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

Размещено на http://www.allbest.ru/

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

Ипатов Юрий Аркадьевич

Новиков Петр Сергеевич

Шургин Алексей Иванович

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

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

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

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

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

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

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

Для биологического анализа необходимо обработать изображение (рис. 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а.

Библиография

электрофорез полимеразный цепной реакция

1. SaikiR.K., GelfandD.H., StoffelS., HiguchiR. // Science. 1988. V. 239. N 2. Р.487.

2. П.С. Новиков, О.В. Шейкина, Т.Н. Милютина Изменчивость плюсовых деревьев сосны обыкновенной на архиве клонов по ISSR-маркерам // Вестник МарГТУ. Радиотехнические и инфокоммуникационные системы. - Йошкар-Ола: МарГТУ- № 3, 2011. - С. 82-87.

3. Самарский, А.А. Численные методы / А.А. Самарский, А.В. Гулин - М.: Наука, 1989. -432 с.

4. Тихонов, В.И. Оптимальный прием сигналов/ В.И. Тихонов - М.: Радио и связь, 1983. - 320 с.

5. Прэтт У. Цифровая обработка изображений: в 2 кн. М.: Мир, 1982.- 790 с.

6. Ипатов Ю.А., Кревецкий А.В., Шмакин В.О. Проектирование распределенной наземной системы мониторинга за лесными пожарами // NB: Кибернетика и программирование. - 2013. - 2. - C. 20-28. URL: http://www.e-notabene.ru/kp/article_8309.html.

7. Ю.А. Ипатов, А.В. Кревецкий Обнаружение текстурных переходов на изображениях препаратов клеточной структуры древесины // Программные системы и вычислительные методы. - 2013. - 1. - C. 106-115. DOI: 10.7256/2305-6061.2013.01.9.

8. А.Г. Коробейников, С.С. Кувшинов, С.Ю. Блинов, А.В. Лейман Анализ принципов создания и работы стеганографических алгоритмов // Программные системы и вычислительные методы. - 2012. - 1. - C. 28-36.

Размещено на Allbest.ru