Анализ дифференциального метилирования геномов методами непредвзятого скрининга
Скачать 396.16 Kb.
|
Экстракция и анализ геномной ДНК. Экстракцию ДНК клеточных линий проводили стандартным методом (Sambrook J. et al., 1989). Рестрикцию ДНК, лигирование с адаптерами, полимеразную цепную реакцию и электрофорез проводили стандартными методами (Sambrook J. et al., 1989). Использовали рестриктазы SmaI, XmaI (Сибэнзим, Россия), ДНК лигазу T4 (Fermentas, Литва), ДНК-полимеразу Taq (Сибэнзим, Россия). Валидацию статуса метилирования изучаемых локусов проводили метилспецифическим секвенированием (Hajkova P., 2002) и метилспецифическим анализом конформации однонитевых фрагментов (Bianco T., 1999). Реакцию автоматического прямого секвенирования проводили на приборе ABI Prism 3100, с использованием BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, США) по протоколам производителя. Разработка алгоритмов и программного обеспечения непредвзятого скрининга дифференциального метилирования геномов. При разработке программ использованы методы модульного, объектно-ориентированного и динамического программирования, цифровой обработки сигналов, численной оптимизации, разработки пользовательских интерфейсов, а также реляционная методология организации хранилища данных. Использован портал организации совместной работы http://www.assembla.com, в состав которого интегрированы репозиторий с системой управления версиями, система управления задачами, пополняемая участниками проекта база знаний. Целевой платформой разработки выбрана программная платформа .NET Framework для Microsoft Windows, позволяющая вести разработку с использованием различных специализированных языков программирования. Программы написаны на языках программирования C# и С++. Для хранения данных продуктов обработки исходной последовательности генома эндонуклеазой рестрикции используется реляционная система управления базами данных Microsoft SQL Server. Возможность задания индексов по предварительно вычисленным фланкирующим последовательностям фрагментов позволяет значительно ускорить выборку данных. Реляционная методология организации хранилища данных характеризуется простотой структуры данных, удобным табличным представлением и возможностью использования формального аппарата алгебры отношений и реляционного исчисления для обработки данных. Программное обеспечение планирования и анализа результатов экспериментов. В качестве инструмента для дизайна эксперимента по амплификации интерметилированных сайтов (АИМС) путем предварительного компьютерного моделирования ожидаемых результатов использована собственная программа компьютерной симуляции AIMS in silico. В  модифицированном методе АИМС (схема оригинального метода представлена на Error: Reference source not found) разделение флуоресцентно меченых продуктов проводили капиллярным электрофорезом на приборе ABI Prism 3100 (Applied Biosystems, США). Результаты фрагментного анализа ДНК обрабатывали с использованием программы GeneMapper 2.0 (Applied Biosystems, США) и собственного программного обеспечения PeakPick. Геномную идентификацию дифференциально метилированных участков ДНК, выявляемых методом АИМС, проводили методом рестриктазного картирования на основе моделей, сформированных программой AIMS in silico. В качестве материала для компьютерного моделирования эпигеномных экспериментов использовали последовательность генома человека hg19, GRCh37 Genome Reference Consortium http://www.ncbi.nlm.nih.gov/projects/genome/assembly/ grc/ в файлах формата FASTA (http://www.ncbi.nlm.nih.gov/BLAST/fasta.shtml). Данные электрофореграмм получены из выходных файлов автоматического генетического анализатора ABI Prism 3100 (Applied Biosystems, США) формата ABIF (http://www.appliedbiosystems.com). III. РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ. 3.1. Современный алгоритм непредвзятого скрининга дифференциального метилированния геномов Идеологическую и методологическую основу проведенного исследования составил специально разработанный алгоритм скрининга дифференциального метилирования геномов. Алгоритм состоит из трех основных блоков: 1) генноинженерный метод формирования репрезентаций эпигеномов; 2) метод высокоразрешающего разделения отдельных элементов репрезентации (исследуемых фрагментов ДНК); 3) методы компьютерного моделирования результатов экспериментов, планирования in silico, анализа результатов исследований in vitro. В качестве основного подхода к формированию репрезентаций эпигеномов был принят метод, в наибольшей степени удовлетворяющий сформулированным требованиям – непредвзятому характеру скрининга, возможности стандартизации процедуры эксперимента, простоте и экономичности протоколов, возможности быстрой и технически несложной идентификации геномной принадлежности дифференциально метилированных участков ДНК, и возможности тестирования большого количества (десятков) образцов в одном эксперименте. Из опубликованных к настоящему времени методов таким требованиям в наибольшей степени удовлетворяет АИМС - амплификация интерметилированных сайтов (Frigola J. et al., 2002). Серьезным недостатком метода АИМС является техническая сложность геномной идентификации идентификации продуктов реакции. Это относится как к методу визуализации (радиоавтография низкопроцентного денатурирующего полиакриламидного геля большого размера с последующим совмещением изображения с реальным гелем), так и к методу определения геномной принадлежности интересующих фрагментов (клонирование амплификатов, содержащихся в элюатах интересующих областей, с последующим секвенированием множества клонов). В идеале метод, генерирующий такое значительное количество целевых фрагментов ДНК, как АИМС (в чем его несомненный плюс) должен быть реализован на платформе, обеспечивающей разрешение фрагментов ДНК по длине с точностью до одного нуклеотида. Такая платформа представлена на сегодняшний день капиллярным электрофорезом в формате фрагментного анализа. Примененный способ дискриминации продуктов АИМС, синтезированных с флуоресцентно меченых праймеров, с помощью многоканального капиллярного электрофореза значительно повышает разрешающую способность метода и избавляет от необходимости использования радиоактивно меченых материалов. Высокая разрешающая способность капиллярного электрофореза в сочетании с компьютерным обеспечением анализа результатов позволяют с высокой точностью определять длины выявляемых дифференциально метилированных фрагментов генома. Такая информация может использоваться для идентификации геномной локализации интересующих фрагментов in silico, что исключает этапы физической изоляции их из геля, клонирования и секвенирования, снижая, тем самым, временные затраты, трудоёмкость и ресурсоёмкость исследования. В  целом, сочетание преимуществ АИМС, капиллярного электрофореза и компьютерного анализа позволяют разработать оригинальный современный алгоритм непредвзятого скрининга дифференциального метилирования геномов (Error: Reference source not found). 3.2. Алгоритм и компьютерная программа моделирования результатов экспериментов АИМС с использованием баз данных нуклеотидных последовательностей В классическом варианте реализация АИМС приводит к формированию сложных фингерпринтов, количественный и качественный состав которых практически невозможно предугадать. С точки зрения идеологии непредвзятости скрининга метилирования такая ситуация оптимальна. В то же время получение четких информативных наборов продуктов АИМС возможно далеко не при любых условиях экспериментов. Технология АИМС исключительно консервативна. На этапе лигирования фланкирующие участки всех продуктов рестрикции образцов геномной ДНК унифицируются с точки зрения последующей ПЦР, проводимой с универсальных праймеров. Нуклеотидный состав праймеров на 90-100% определяется составом адаптеров. Жесткие условия ПЦР, дающие высокую специфичность реакции и отсутствие неспецифичных продуктов, обеспечивают высокий уровень стандартизации метода и повторяемости результатов. Во всем протоколе АИМС существуют только два относительно гибких сегмента. Это выбор изошизомеров рестриктаз, используемых при подготовке репрезентаций геномов, и вариации 3’-концов праймеров, которые могут различаться по длине на 1-4 нуклеотидов. Эти вариабельные участки праймеров («удлинители») могут различаться не только по длине, но и по нуклеотидному составу. Таким образом, дизайн эксперимента включает в себя лишь выбор рестриктаз, протяженности и нуклеотидного состава удлинителей. От этих факторов зависит количество возможных результирующих фрагментов АИМС и их геномная принадлежность. Тем не менее, разнообразие картин АИМС, возникающее при разных сочетаниях указанных переменных параметров, исключительно велико. Адекватный дизайн экспериментов АИМС можно обеспечить путем предварительного математического моделирования ожидаемых результатов. Такая возможность обеспечивается разработанной программой компьютерной симуляции AIMS in silico. Алгоритм работы программы включает моделирование полного набора фрагментов ДНК, получаемых при обработке исследуемого генома рестриктазой, выбранной для формирования геномной библиотеки (традиционно для АИМС это рестриктаза XmaI). Полученный набор фрагментов представляет собой полную репрезентацию АИМС, служащую основой для моделирования возможных исходов АИМС in vitro при различных условиях. Полученная полная репрезентация виртуально подвергается лигированию с предложенным пользователем адаптером и ПЦР с универсальными праймерами, содержащими заданные удлинители на 3'-конце. Получаемые полные наборы возможных продуктов АИМС для заданных условий можно характеризовать с количественной и качественной точек зрения, причем в последнем случае реальное клонирование и секвенирование отдельных фрагментов ДНК заменяется сравнением с уже существующими геномными базами данных. Таким образом, алгоритм моделирования результатов экспериментов воспроизводит лабораторный протокол АИМС. Интерактивный интерфейс разработанной на основе этого алгоритма программы AIMS in silico позволяет осуществлять ввод файлов, содержащих нуклеотидные последовательности, предназначенные для анализа, и выбор пользователем исходных условий эксперимента. К вводимым условиям эксперимента относятся: сайт узнавания рестриктазы, используемой для скрининга дифференциального метилирования, нуклеотидный состав короткой и длинной частей адаптера, нуклеотидный состав специфического удлинителя универсального праймера. Предусмотрена возможность ограничения длин результирующих фрагментов АИМС, в зависимости от способа их физического разделения. Выходные данные: общее количество продуктов ПЦР, которые могут быть получены при проведении АИМС с заданными условиями, их распределение по длинам, нуклеотидные последовательности каждого из ПЦР-продуктов. В графическом виде информация представляется как виртуальная хроматограмма капиллярного электрофореза, построенная в виде суперпозиции пиков Гаусса, соответствующих индивидуальным продуктам АИМС. Выбор отдельного пика при помощи курсора вызывает информацию о нем в специальном окне (Рис. 1). При составлении программы использованы преимущества технологии баз данных. Собственно анализ целевых последовательностей начинается с создания пользователем базы данных, содержащей информацию о каждом участке изучаемого генома, возникающем после обработки геномной ДНК эндонуклеазой рестрикции: длину продукта рестрикции и его нуклеотидный состав. После того, как такая база данных создана и сохранена, возможность быстрого обращения к ее содержимому обеспечивает высокую скорость симуляции результатов АИМС для различных удлинителей универсального праймера.  Рис. 1. Пример окна программы AIMS in silico с выходными данными. 3.3. количественная и качественная характеристика продуктов АИМС, получаемых при различных экпериментальных условиях Характеристика распределения сайтов узнавания эндонуклеазы рестрикции XmaI в геноме человека. В классическом варианте АИМС используется эндонуклеаза рестрикции XmaI с сайтом узнавания C^CCGGG. Программа AIMS in silico позволила оценить количественный и качественный состав продуктов гидролиза геномной ДНК человека этой рестриктазой. Общее число продуктов гидролиза составило 374146. Распределение последовательности CCCGGG по геному оказалось в целом случайным. Характер распределения сайтов узнавания XmaI одинаков на всех хромосомах. Митохондриальный геном последовательностей CCCGGG не содержит. Случайное пространственное распределение сайтов CCCGGG отличает рестриктазу XmaI от часто используемой в эпигенетических исследованиях HpaII, сайты узнавания которой преимущественно концентрируются в областях CpG-островков. Случайное распределение сайтов CCCGGG дало основание предположить, что использование XmaI в экспериментах АИМС позволит проводить дифференциальный анализ метилирования самых разнообразных, с точки зрения структуры и функции, геномных элементов. Это предположение полностью подтверждается более детальным анализом продуктов рестрикции XmaI. Характер распределения XmaI-фрагментов, ограниченных порогом длин эффективного разделения капиллярным электрофорезом, показан на Рис. 2, А. Общее число продуктов АИМС длиной до 600 нуклеотидов составило 73102.  Рис. 2. Виртуальные электорофореграммы АИМС, полученные с помощью программы AIMS in silico: А - характер распределения продуктов АИМС без удлинителей, отражающий распределение в геноме XmaI-фрагментов длиной до 600 н.п. Б, В – виртуальные электорофореграммы АИМС с удлинителями CG и CCG, соответственно. Качественный анализ полученной репрезентации генома показал, что она обогащена повторяющимися последовательностями с длинами коровых единиц от 100 до 1000 нуклеотидов, принадлежащими самым разнообразным классам рассеянных и тандемных повторов: Alu, LTR, прителомерным повторам, неклассифицированным повторам низкой сложности и др. Некоторые близкородственные, или, по крайней мере, в значительной степени гомологичные высококопийные повторяющиеся последовательности образовали выраженные кластеры продуктов АИМС (несколько высоких, до 1500 копий, пиков на хроматограмме, в области длин до 80 нуклеотидов, Рис. 2, А). Присутствуют в разных соотношениях и другие классы геномных последовательностей (уникальные межгенные и интронные участки, CpG-островки и т.д.). Количественная иерархия продуктов АИМС генома человека Для ограничения сложности картин АИМС с целью облегчения анализа электрофореграмм при проведении ПЦР используются праймеры, удлиненные на несколько нуклеотидов. Обычно выбор удлинителей праймеров проводится случайным образом и утверждается либо отвергается по результатам экспериментов. Авторы метода АИМС изначально постулировали, что при прочих равных условиях добавление каждого дополнительного нуклеотида к универсальным праймерам предполагает снижение сложности результирующей картины (т.е. снижение количества продуктов АИМС) примерно в 5 раз для C,G и в 3 раза для A,T. Такое предположение сделано, исходя из представленности нуклеотидов в геноме человека. Случайный выбор удлинителей с использованием такой приблизительной эмпирической формулы для оценки сложности репрезентаций геномов снижает ценность АИМС как универсального метода скрининга, необходимой характеристикой которого должен стать высокий уровень стандартизации. Кроме того, использование различных наборов нуклеотидов в качестве удлинителей должно сдвигать спектральный состав результирующих репрезентаций либо в А,Т-богатые области, либо в области CpG-островков, и это влияние должно, несомненно, учитываться не только при анализе результатов, но и на самом первом этапе исследования – в процессе определения дизайна эксперимента. Для предоставления исследователям возможности обоснованного выбора удлинителей с помощью программы AIMS in silico сформировано графическое представление количественной иерархии продуктов АИМС генома человека (Рис. 3). Приведенная гистограмма ясно показывает, что для подавляющего числа семейств удлинителей достаточно трех нуклеотидов, чтобы снизить сложность картины АИМС таким образом, чтобы продукты ПЦР можно было четко дифференцировать после разделения в полиакриламидном геле. Речь идет о количестве продуктов, не превышающем второго порядка в пределах длин до одной тысячи нуклеотидов. Предоставляемая программой симуляции АИМС in silico возможность построения иерархического дерева конечных продуктов в зависимости от длин и нуклеотидного состава удлинителей не только позволяет уточнить функцию снижения сложности для каждого конкретного случая, но и демонстрирует удивительные результаты. Так, характер первого же нуклеотида в удлинителе практически не влияет на сложность картины: присоединение нуклеотида А, Т, С или G приводит к ограничению количества продуктов соответственно до 8.5, 5.5, 7 и 7 тысяч. При этом разница результатов присоединения А и Т значительно более выражена, чем разница между А/Т и C/G. Чрезвычайно выраженные различия наблюдаются на уровне удлинителей из двух нуклеотидов. Так, в семействе удлинителей с первым нуклеотидом Т значения результирующих количеств продуктов АИМС для ТА, ТС, ТG и ТТ составляют 20, 152, 191 и 2547 соответственно, т.е. разница в степени сложности репрезентаций составляет целых два порядка. Снижение сложности до состояния практически полной неинформативности для некоторых удлинителей наступает уже на уровне трех нуклеотидов. В частности, АИМС in silico с одним из удлинителей ААТ, ТАТ, ТСА и GТА приводит к возникновению лишь одного продукта реакции. В то же время, использование некоторых других тринуклеотидных удлинителей приводит к едва заметному снижению сложности исходной картины: например, для АGG количество фрагментов составляет 2829, для ТТС – 2161. Экспериментальное определение этих параметров потребовало бы невероятных затрат труда, времени и материалов, а неопределение наверняка привело бы к ложной интерпретации опытных данных. В  Рис. 3. Иерархическая диаграмма продуктов АИМС до 1 т.п.н., предсказанных для классических условий эксперимента. Синим цветом обозначено общее количество первичных фрагментов (при отсутствии удлинителей) – около 100 тыс. Красным цветом указаны количества продуктов, получаемых при использовании однонуклеотидных удлинителей, желтым – двухнуклеотидных, зеленым – трехнуклеотидных. Диапазон предсказанных количеств фрагментов составляет от 1 до 100000, поэтому ось абсцисс представлена логарифмической шкалой. ажным качественным параметром, определяющим дизайн АИМС, является характер распределения продуктов ПЦР по длинам в анализируемой области. Распределение должно быть в целом равномерным, длины отдельных фрагментов в идеале должны различаться не менее чем на два нуклеотида. В этом случае будет обеспечиваться четкая идентификация пиков на хроматограмме капиллярного электрофореза и однозначность результатов анализа. Не менее важным фактором, влияющим на выбор удлинителей для проведения АИМС, является их полная гомология с теми или иными структурно-функциональными элементами генома, определяющая конечный качественный состав анализируемых фрагментов ДНК. Качественный состав продуктов АИМС генома человека и дизайн экспериментов АИМС при помощи программы AIMS in silico Нуклеотидные последовательности продуктов АИМС, предоставляемые программой AIMS in silico, позволяют однозначно локализовать анализируемые фрагменты и, пользуясь доступными базами данных, определить принадлежность каждого из них к тому или иному классу геномных последовательностей. Используя эту возможность, можно целиком осуществить обоснованный дизайн эксперимента АИМС в зависимости от целей исследования и доступной приборной базы, ни разу не прибегая к планировочным экспериментам. Ниже приводится конкретный пример, демонстрирующий возможности дизайна экспериментов АИМС in silico. Поставим задачу определения условий АИМС для изучения характера метилирования не менее 100 CpG-островков в образцах условно нормальной и опухолевой тканей молочной железы. Планирование такого эксперимента будет включать следующие этапы. Прежде всего, предложим рабочую гипотезу, ограничивавшую качественный состав потенциальных удлинителей. Гипотеза заключается в простом и очевидном предположении, что для обогащения фракции продуктов АИМС фрагментами классических CpG-островков необходимо использовать удлинители, содержащие по крайней мере один динуклеотид CG. Среди двухнуклеотидных такой удлинитель, очевидно, лишь один – CG, и результаты компьютерной симуляции АИМС (Рис. 2, Б) показывают, что его использование в эксперименте нецелесообразно: в пределах длин до 600 нуклеотидов он создает 308 фрагментов, причем их распределение не соответствует сформулированным выше требованиям. Среди тринуклеотидных удлинителей последовательностей, содержащих CG, восемь: ACG, CCG, CGA, CGC, CGT, CGG, GCG, и TCG. Количество продуктов АИМС длиной до 600 нуклеотидов для каждого из них - соответственно 10, 71, 3, 73, 3, 43, 65, и 21. Ни один из этих удлинителей не дает возможности анализа 100 фрагментов АИМС, поэтому выберем из них три, обещающих достаточно высокую информативность – CCG, CGC и GCG. Рассмотрение результатов виртуального электрофореза для удлинителя CGC (Error: Reference source not found, Б) показывает, что его высокая информативность, предполагаемая исходя из количества продуктов АИМС, обманчива. Из 73 фрагментов 60 группируются в области до 250 нуклеотидов; формируются 8 кластеров фрагментов одной длины, относящихся к 2-10 геномным локусам. Кластер фрагментов с длинами 222 н.п. включает в себя 5 продуктов АИМС, представляющих участки сателлитных повторов и не имеющих отношения к CpG-островкам. Одиннадцать фрагментов с общей длиной 44 н.п. и три фрагмента с общей длиной 46 н.п. представлены высокогомологичными участками непромоторных CpG-островков. Что же касается предсказанных продуктов АИМС, получаемых при использовании удлинителей CCG и GCG, то их распределение более равномерно, расстояние между большинством пиков превышает два нуклеотида, обе репрезентации содержат лишь по одному кластеру из трех продуктов АИМС (Error: Reference source not found, А и В). По этим признакам удлинители CCG и GCG в большей степени соответствуют оптимальным параметрам репрезентации, чем CGC. На следующем этапе необходимо определить качественный состав продуктов АИМС, получение которых возможно с выбранными удлинителями CCG и GCG. Используя нуклеотидные последовательности продуктов АИМС, с помощью программы BLAT определяем их геномную локализацию и принадлежность к классу CpG-островков. Р  ассмотрим результаты этого этапа более подробно на примере удлинителя CCG. Пример соответствующей виртуальной хроматограммы представлен на Error: Reference source not found, А. Состав предсказанной репрезентации генома отражен на Error: Reference source not found. Анализ показывает, что абсолютное большинство результирующих продуктов АИМС при удлинителе CCG принадлежит классическим протяженным CpG-островкам с высоким содержанием CpG-динуклеотидов: от 25 до 400, в среднем 125. Аналогичные результаты дает количественный и качественный анализ продуктов АИМС при удлинителе GCG. Таким образом, совместный анализ образца биологического материала с использованием удлинителей CCG и GCG позволяет определять состояние метилирования 136 геномных участков, 122 из которых представлены CpG-островками. Дизайн эксперимента АИМС с заранее заданными параметрами завершен.  3.4. Алгоритм и компьютерная программа геномной идентификации выявляемых дифференциально метилированных участков ДНК Особенности первичной нуклеотидной последовательности продуктов АИМС не позволяют проводить их точную идентификацию по молекулярной массе в соответствии с распределением, предсказанным in silico. Исключительно обогащенные G,C-нуклеотидами фрагменты ДНК мигрируют в среде разделения со скоростями, отличными от таковых, характерных для фрагментов с усредненным нуклеотидным составом. Кроме того, использование денатурирующей среды для проведения капиллярного электрофореза зачастую вызывает различия электрофоретической подвижности G,C-обогащенных комплементарных фрагментов ДНК, приводя к двоению соответствующих сигналов на электрофореграмме. Традиционный дизайн АИМС предполагает клонирование и секвенирование интересующих фрагментов ДНК для определения их геномной принадлежности. Пользуясь тем, что определение нуклеотидной последовательности генома человека уже практически полностью завершено и, соответственно, последовательности возможных продуктов АИМС также известны, этапы клонирования и секвенирования заменены последовательными этапами рестриктазного картирования in silico и in vitro. Для осуществления рестриктазного картирования in silico был введен дополнительный модуль в программу AIMS in silico. Моделирование рестриктазного картирования проводится на основе полной репрезентации АИМС, предсказанной для конкретных условий эксперимента, с использованием рестриктаз, выбираемых из заранее сформированного списка. Сигналы виртуальной электрофореграммы, соответствующие «гидролизованным» продуктам АИМС при выборе конкретной рестриктазы, определяются по изменению цвета. Виртуальная электрофореграмма результатов рестриктазного картирования продуктов АИМС интерактивна; предусмотрена возможность копирования полной нуклеотидной последовательности выбранного продукта АИМС при помощи курсора (Рис. 4).   Разработан алгоритм интеграции рестриктазного картирования in silico и in vitro. Картирование реально получаемых продуктов АИМС проводится, как и в случае компьютерного моделирования, на основе полной репрезентации АИМС, получаемой для конкретных условий эксперимента. Для создания полной репрезентации in vitro из лабораторного протокола АИМС исключается этап гидролиза ДНК метилчувствительной рестриктазой SmaI. Таким образом моделируется состояние полного метилирования изучаемого генома. 3.5. Анализ дифференциального метилирования ДНК клеточных линий рака молочной железы При сравнении электрофореграмм продуктов АИМС геномов клеточных линий с удлинителем CCG выявлены дифференциально метилированные локусы в диапазонах длин 150-167, 167-176 и 176-190 п.н. (Рис. 5).  Рис. 5. Сравнение электрофореграмм продуктов АИМС с удлинителем CCG клеточных линий РМЖ. Моделирование с помощью программы AIMS in silico предсказывает получение в эксперименте АИМС с удлинителем CCG в указанном диапазоне трех групп продуктов: 1)ccg164, ccg165, ccg168; 2) ccg176, ccg177.1, ccg177.2, и 3) ccg189. Краткая характеристика геномной принадлежности локусов продуктов приведена в таблице 1. Виртуальная электрофореграмма продуктов представлена на Рис. 6. С помощью программы AIMS in silico получен план рестриктазного картирования (табл. 2) Таблица 1. Характеристика геномной принадлежности локусов.
* расстояние до TSS (transcription start site – сайт инициации транскрипции)  Рис. 6. Результаты моделирования продуктов АИМС с удлинителем CCG в диапазоне длин 160-190 п.н. Таблица 2. План рестриктазного картирования
Гидролиз продуктов АИМС эндонуклеазами рестрикции и сравнение электрофореграмм позволили однозначно определить положение пиков продуктов на электрофореграммах (Рис. 7).  Рис. 7. Пример позиционирования пиков продуктов на электрофореграмме продуктов АИМС клеточной линии MCF7 до (черным) и после гидролиза BssHII (красным). Зеленой стрелкой обозначен пик продукта, соответствующий ccg177.2, синей стрелкой – продукт гидролиза. С использованием метода метилспецифического секвенирования (Рис. 8) построены карты метилирования (Рис. 9) выявленных локусов дифференциального метилирования.  Рис. 8. Пример электрофореграммы метилспецифического секвенирования локуса ccg177.2 в клеточной линии MCF7. Стрелками указаны сайты узнавания рестриктазы SmaI в полностью метилированном состоянии (после бисульфитной обработки – TTCGGG).  Рис. 9. Карты метилирования локуса генома человека, содержащего фрагмент АИМС ccg177.2. Белыми кружками показаны неметилированные CpG-динуклеотиды, черными – метилированные, серыми – частично метилированные, желтыми – CpG-динуклеотиды, состояние метилирования которых определить не удалось. Сравнение результатов, полученных методами АИМС и метилчувствительного секвенирования, показало совпадение оценок статуса метилирования сайтов SmaI по всем исследуемым локусам во всех образцах, кроме локуса ccg164 (PURB) в линии HBL и локусов ccg168 и ccg177.1 в линии T 47D, для которых метилспецифическое секвенирование не подтвердило наличия метилирования. Возможно, наблюдаемые расхождения объясняются более высокой чувствительностью метода АИМС по сравнению с секвенированием по Сэнгеру. Сравнение электрофореграмм полной репрезентации и продуктов АИМС клеточной линии HBL 100 (Рис. 10) показало, что уровень сигнала продукта АИМС, соответствующего локусу ccg168, примерно вдвое ниже, чем на электрофореграмме полной репрезентации, что согласуется с данными карты метилирования по локусу ccg168 для клеточной линии HBL 100 – один сайт SmaI полностью метилирован, второй в полуметилированном состоянии (Рис. 11). Уровень сигнала продукта АИМС, соответствующего локусу ccg164, значительно ниже, чем на электрофореграмме полной репрезентации, что позволяет предположить наличие в образце минорной субпопуляции геномов (мозаичность), в которых соответствующие CpG-динуклеотиды находятся в метилированном состоянии.  Рис. 10. Сравнение электрофореграмм полной репрезентации (черным) и продуктов АИМС (красным) клеточной линии HBL-100.  Рис. 11. Карты метилирования локуса генома человека, содержащего фрагмент АИМС ccg168 (ANK3). Здесь и далее белыми кружками показаны неметилированные CpG-динуклеотиды, черными – метилированные, серыми – в полуметилированном состоянии, желтым – CpG-динуклеотиды, состояние метилирования которых определить не удалось. Сравнение электрофореграмм полной репрезентации и продуктов АИМС клеточной линии T 47D (Рис. 12) показывает, что уровень сигналов продуктов АИМС, соответствующих локусам ccg168 и ccg177.1 значительно ниже, чем на электрофореграмме полной репрезентации, что может быть объяснено тем же явлением.  Рис. 12. Сравнение электрофореграмм полной репрезентации (черным) и продуктов АИМС (красным) клеточной линии T 47D. Метилспецифический анализ конформации однонитевых фрагментов АИМС ccg168 показал в образце клеточной линии HBL 100 наличие метилированных молекул, обладающих измененной электрофоретической подвижностью относительно неметилированной ДНК (). Для образца клеточной линии T 47D показано наличие только фракций молекул с измененной электрофоретической подвижностью, в то время как результаты метилспецифического секвенирования (Рис. 11) не показывают наличия метилирования ни по одному CpG-динуклеотиду, что может быть объяснено высокой гетерогенностью клеточной линии по статусу метилирования в связи с повреждением механизма поддержания метилирования в этой клеточной линии (Ting A.H. et al., 2006).  Рис. 13. Метилспецифический анализ конформации однонитевых фрагментов по локусу АИМС ccg168. 1, 3, 4, 5, 7 – образцы неметилированной ДНК. 2 – образец ДНК клеточной линии HBL 100. 6 – образец ДНК клеточной линии T 47D. Стрелками указаны сигналы, соответствующие неметилированным молекулам. Проведенное исследование говорит о высокой чувствительности разработанного метода анализа дифференциального метилирования геномов. Анализ геномов клеточных линий рака молочной железы ZR 75 1, HBL 100, HS 578 T, BT 474, MCF7 и T 47D выявил 6 дифференциально метилированных геномных локусов, которые ранее не изучались с точки зрения дифференциального метилирования: участки CpG-островков, расположенных в первых интронах генов ANK3 и MAFK, в промоторной области гена C2CD2, в межгенных областях на хромосомах 12q13.13 и 13q32.1, а также участок первого экзона гена PURB, не являющийся CpG-островком. Выявление дифференциального метилирования неканонической локализации подтверждает непредвзятость разработанного алгоритма скрининга. |
Похожие:
 | Об утверждении порядков проведения пренатальной (дородовой) диагностики нарушений развития ребенка, универсального аудиологического... |  | Анализ правовой базы для уточнения понятия социального предприятия и связанных с ним терминов, а также инструментария государственного... |
| Как правило, задачи эти со многими неизвестными, со сложным и вариативным составом исходных данных и возможных решений. Чтобы уверенно... | | Клиентом признается нп огр и всеми её членами добровольной инициативой и стимулируется внутри нп огр управленческими методами. Предпосылкой... |
| Изучение стереотипов женского поведения методами психосемантики на материале России и сша1 | | Целью практических занятий является ознакомление студентов с теоретическими положениями дисциплины, с признаками и методами кодирования... |
| Учебно-методическое пособие предназначено для студентов, изучающих курсы «Прикладная социология», «Статистический анализ социальной... | | Широкополосный анализатор поляризации лазерного излучения на основе нанографитовой плёнки |
| Способы диагностики перикардитов; Принципы про ведения дифференциального диагноза перикардита с другими заболева- ниями сердца; Вопросы... | | Анализ потенциальныхисточников информации для мониторинга и оценки эффективности налоговых льгот (в том числе анализ судебной практики;... |