Создание автоматизированной системы анализа изображений полимеразной цепной реакции

Ипатов Юрий Аркадьевич кандидат технических наук доцент, кафедра информатики, Поволжский государственный технологический университет 424000, Россия, Республика Марий Эл, г. Йошкар-Ола, пл. Ленина, 3 Ipatov Yurii Arkad'evich PhD in Technical Science Associate Professor of the Department of Informatics at the Volga State University of Technology 424000, Russia, Mari El, Yoshkar-Ola, pl. Lenina, d. 3 philsilver@mail.ru Другие публикации этого автора     Новиков Петр Сергеевич ассистент, кафедра лесной селекции, недревесных ресурсов и биотехнологии, Поволжский государственный технологический университет 424000, Республика Марий Эл, г. Йошкар-Ола, пл. Ленина, дом 3. Novikov Petr Sergeevich Assistant of the Department of Forest Breeding, Biotechnology and Non-timber Resources at the Volga State University of Technology 424000, Russia, Mari El, Yoshkar-Ola, pl. Lenina, d. 3 philsilver@mail.ru

Шургин Алексей Иванович доцент, кафедра лесной селекции, недревесных ресурсов и биотехнологии, Поволжский государственный технологический университет 424000, Республика Марий Эл, г. Йошкар-Ола, пл. Ленина, дом 3. Shurgin Aleksei Ivanovich Associate Professor of the Department of Forest Breeding, Biotechnology and Non-timber Resources at the Volga State University of Technology 424000, Russia, Mari El, Yoshkar-Ola, pl. Lenina, d. 3 philsilver@mail.ru

Введение

Современные подходы и методы ДНК-анализа позволяют получать точные и достоверные результаты генетической изменчивости и разнообразия видов. Использование ДНК-маркеров дает возможность провести такие исследования и сделать вывод о генетическом разнообразии. Однако некоторые методы является дорогими и не эффективными, в то время как метод ДНК-маркеров показывает хорошие результаты, и доказывают свою эффективность для решения такого рода задач.

Использование полимеразной цепной реакции (ПЦР) позволяет увеличивать малые концентрации определённых фрагментов нуклеиновой кислоты (ДНК) в биологическом материале. При этом, один из ресурсоемких и трудоемких этапов формирования первичных данных, является анализ цифровых изображений результатов химической реакции^[1].

Разработке необходимых теоретических методов и практических реализаций автоматизированного анализа данных будет посвящена настоящая работа.

Анализ изображения гелей, получаемых в результате электрофореза образцов содержащих продукты ПЦР с использованием праймеров, представляет собой трудно формализуемую задачу. На рис.1. представлено изображение продуктов ЦПР. Это связанно не только с большим числом полос выявляемых при электрофорезе, но и с рядом проблем возникающих в результате не возможности стандартизации всех условий проведения электрофореза. Сочетанное использование таких коммерческих программ как RAD-BioMed, FastPCR, Sigma позволяет решить ряд задач, возникающие при анализе изображений гелей, но не все ^[2]. При этом, этап формирования первичных данных, является одним из ресурсоемких и трудоемких, а возможность её автоматизации, на сегодняшний не решается в полной мере.

Рис. 1. Изображение электрофореграммы ПЦР-продуктов
ДНК сосны с праймером (CA)₆AGG.

Стандартные реализации разработчиков программного обеспечения, как отечественных, так и зарубежных фирм, не всегда учитывают специфику решаемых задач. Предлагаемая разработка автоматизированной системы на базе методов распознавания и анализа цифровых изображений гелей, позволит минимизировать трудоемкие операции выполняемых человеком и повысит качество проводимых работ.

Синтез алгоритма

Для биологического анализа необходимо обработать изображение (рис.1.) и определить положения объектов, являющихся продуктами ПЦР. Данный тип электрофореграммы был получен с помощью оборудования и программного комплекса RAD-BioMed. Анализ изображений показывает, что праймеры и продукты ПЦР, расположены в вертикальных пространственных областях, которые не пересекаются. Такое расположение дает возможность, с помощью стандартного программного обеспечения автоматически распознать и назначить уникальные линии, определяющие эти области.

Алгоритм будут состоять из следующих шагов:

1. Выполняем операцию накопления яркостных отсчетов вдоль вертикали кадра, а затем усредняем полученный сигнал. Получившееся изображение сигнала представлено на рис.1.

Рис. 2. Изображение сигнала накопленной и усредненной яркости изображения вдоль вертикали кадра (рис.1.), а также функция интерполяции

Форма полученного сигнала свидетельствует о том, что исходное изображение содержит мультипликативную помеху. Для её устранения необходимо, получить оценку данного типа помехи. Экспериментальный анализ показывает, что для этих задач можно использовать полином второй степени ^[3], чтобы получить интерполирующий сигнал.

2. Получаем разностный сигнал (рис.3) для которого определяем локальные максимумы по определенному пороговому значению.

Рис. 3. Разностный сигнал для функций по рис.2.

Результирующий массив данных, характеризует пространственные координаты вертикальных линий на изображении вдоль горизонтали кадра. Недостающие вертикальные координаты для линий будет определяться двумя точками, задаваемыми вручную. Таким образом, удается определить пространственные координаты начала и конца всех линий на изображении.

3. Далее используя принцип согласованной фильтрации ^[4] и имея априорно известную импульсную характеристику двухмерного фильтра ^[5] находим пространственные координаты для продуктов ПЦР.

Рис.4. Участок изображения (рис.1) с найденными продуктами ПЦР

Результат обнаружения можно посмотреть на рис.4.

Выводы

Предлагаемая разработка автоматизированной системы на базе методов распознавания и анализа цифровых изображений гелей ПЦР, позволит минимизировать трудоемкие операции выполняемых человеком и повысит качество проводимых работ. Сравнительный анализ ручных методов и созданного подхода выигрывает по времени на полтора порядка. При этом учитывая показатели массовости и точности, повышают экономический эффект от использования данной разработки. Синтезированный алгоритм имеют свидетельство об официальной регистрации программы для ЭВМ Федеральная службы по интеллектуальной собственности, патентам и товарным знакам РФ №2012619172 «BioImageGeles PCR Analysis» от 11.10.2012 года.

Работа выполнена при финансовой поддержке РФФИ, проект 13-01-00427а.