Способ определения метилирования сайтов pucgpy регуляторных областей генов-онкомаркеров колоректального рака методом glad-пцр-анализа и набор олигонуклеотидных праймеров и флуоресцентно-меченых зондов для осуществления указанного способа

Способ определения метилирования сайтов pucgpy регуляторных областей генов-онкомаркеров колоректального рака методом glad-пцр-анализа и набор олигонуклеотидных праймеров и флуоресцентно-меченых зондов для осуществления указанного способа

Иллюстрации

Показать все

Реферат

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

Перспективным инструментом молекулярно-генетической диагностики и мониторинга онкологических заболеваний является эпигенетическая диагностика, то есть идентификация болезни путем выявления эпигенетических маркеров. Эпигенетические маркеры представляют собой химические группы, которые присоединяются к геномной ДНК человека (метилирование цитозиновых оснований с образованием 5-метилцитозина, 5mC) либо к ассоциированным с ней гистоновым комплексам и оказывают влияние на функциональную активность генов. Наиболее распространенным эпигенетическим нарушением при канцерогенезе является локальное гиперметилирование отдельных участков ДНК, так называемых CpG-островков, в регуляторных областях генов-онкосупрессоров, которые в норме лишены метальных групп. В настоящее время показано, что подобное аберрантное метилирование характерно для начальных стадий большинства ненаследственных (спорадических) форм рака, которые в среднем составляют более 90% всех случаев злокачественных новообразований.

Для осуществления эпигенетической диагностики злокачественных новообразований необходимо решить два ключевых вопроса:

а) какой метод детекции следует использовать для определения статуса метилирования генов;

б) метилирование регуляторных участков каких генов ассоциировано с конкретным онкологическим заболеванием?

Для решения первого вопроса, а именно для анализа статуса метилирования генов, применяется ряд методов, которые можно разделить на четыре группы.

Первая группа методов основана на бисульфитной конверсии ДНК, влекущей изменения ее нуклеотидного состава. В дальнейшем нужный участок конвертированного генома может быть амплифицирован в ПЦР и секвенирован для определения всех метилированных цитозиновых остатков. Если исследователя интересует статус метилирования конкретного CG-динуклеотида, то используются более простые способы анализа. Так, возможно проведение метилчувствительной ПЦР (другой употребляемый термин - «метилспецифическая ПЦР»), в которой проводятся две реакции с использованием праймеров, подобранных для случаев прохождения или непрохождения конверсии в исследуемом динуклеотиде. Другим широко используемым способом определения метилирования является обработка конвертированной ДНК эндонуклеазами рестрикции (метод COBRA, основанный на появлении или исчезновении сайтов узнавания ферментов в измененной ДНК). К недостаткам описанных методов относятся существенная деградация и значительные потери нуклеиновых кислот при проведении бисульфитной конверсии, что делает проблематичным анализ малого количества материала, трудоемкость, дороговизна исследования и значительные затраты времени [Umer M, Herceg Z. Deciphering the epigenetic code: an overview of DNA methylation analysis methods // Antioxid Redox Signal. - 2013. - V. 18. - P. 1972-1986].

Вторая группа методов включает в себя процедуры гидролиза ДНК метилчувствительными эндонуклеазами. Анализ статуса метилирования генов основан на том, что данные рестриктазы не способны распознавать и расщеплять последовательности, содержащие метилированные CG-сайты (например, фермент HpaII, узнающий и гидролизующий лишь неметилированный тетрануклеотид CCGG). Результаты реакции гидролиза могут быть проанализированы с помощью ПЦР с фланкирующими праймерами. Существенным ограничением для широкого применения этих методов является относительно небольшое число известных эндонуклеаз рестрикции, распознающих последовательности, содержащие CG-динуклеотиды, и при этом чувствительных к метилированию цитозиновых остатков в них [Shen L, Waterland RA. Methods of DNA methylation analysis // Curr Opin Clin Nutr Metab Care. - 2007. - V. 10. - P. 576-581]. Данный подход также не позволяет выявлять частично метилированную ДНК, тогда как именно такой паттерн метилирования CpG-островков характерен для самых ранних стадий канцерогенеза. Кроме того, использование данных методов предъявляет более высокие требования как к чистоте анализируемой ДНК, так и к качеству проведения исследования на всех стадиях, поскольку вывод о наличии метилирования в этих методах делается при отсутствии гидролиза ДНК метилчувствительными эндонуклеазами рестрикции. Однако недогидролиз ДНК может происходить в случаях наличия загрязнений или нарушений условий анализа. Это приводит к ложным выводам о наличии метилирования и как следствие к ложно-положительному диагнозу.

Третья группа методов базируется на использовании метил-связывающих белков (methyl binding domain containing proteins, MBD-белки) для изоляции метилированной ДНК и ее последующей амплификации. Известна, например, технология MethylMeter, разработанная специалистами американской компании RiboMed (патент Канады №2724343, опубл. 15.05.2008 г.), позволяющая достигать высоких показателей чувствительности и специфичности. Данная технология используется производителями некоторых тест-систем, например, G-CIMP DecisionDX американской компании Castle Biosciences. Для защиты данной тест-системы подана заявка США №2014272986, опубликованная 18.09.2014 г. Тем не менее, подобные методы остаются весьма трудоемкими, а для четкого количественного анализа результатов требуется сложное дорогостоящее оборудование.

Четвертая группа методов базируется на использовании сайт-специфических метилзависимых ДНК эндонуклеаз. После их открытия и начала коммерческого производства их препаратов появилась возможность использования этого нового класса ферментов в эпигенетических исследованиях. В отличие от обычных эндонуклеаз рестрикции, метилзависимые ферменты узнают только метилированные последовательности ДНК, что делает их еще одним удобным инструментом изучения статуса метилирования генов. В число ферментов данного класса входят производимые в России компанией ООО «СибЭнзайм» эндонуклеазы GlaI, BisI, BlsI, MteI и ряд других [http://russia.sibenzyme.com/products/m2_type].

Благодаря совпадению последовательности ДНК, являющейся продуктом метилирования de novo, и субстратной специфичности фермента GlaI исследователями ООО «СибЭнзайм» был разработан метод, названный GLAD-ПЦР (GlaI hydrolysis and Ligation Adapter Dependent PCR), который позволяет определять наличие сайтов R(5mC)GY в исследуемом участке ДНК [Кузнецов В. В. с соавт. Способ определения нуклеотидной последовательности Pu(5mC)GPy в заданном положении протяженной ДНК // Патент РФ №2525710 С1, опубл. 20.08.2014, Бюл. № 23]. Метод включает гидролиз ДНК ферментом GlaI, присоединение специфического адаптера к концам гидролизованной ДНК и последующую реакцию ПЦР в режиме реального времени с праймеров, окружающих целевой участок ДНК. Метод позволяет определить наличие одного или нескольких сайтов Pu(5mC)GPy внутри изучаемого фрагмента ДНК. Таким образом, метод GLAD-ПЦР может быть использован для определения статуса метилирования регуляторных участков генов.

Данный метод обладает рядом преимуществ, заключающихся в простоте выполнения, способности выявлять лишь метилированные сайты в сложных смесях ДНК и возможности количественно оценивать их число. Следует также отметить, что метод обладает значительной универсальностью и может быть применен для анализа любых сложных смесей ДНК, например, ДНК из опухолей, содержащей нуклеиновые кислоты не только из малигнантных, но и из нормальных клеток (клетки сосудов, соединительной ткани и т.д.). Апробация GLAD-ПЦР показала, что метод обладает высокой чувствительностью, позволяя детектировать метилированные сайты в сильноразбавленных растворах.

Указанный метод является универсальным методом детекции метилирования сайтов R(5mC)GY, позволяющим сформировать и описать конкретную панель генов, специфичную, например, для колоректального рака. Но в самом описании метода специфичные для различных видов онкологии панели генов отсутствуют.

Рак толстого кишечника и прямой кишки (колоректальный рак) - один из наиболее часто встречающихся видов злокачественных новообразований. В развитых странах мира ежегодно диагностируется более миллиона случаев этого заболевания, а летальность от рака этой локализации составляет более 33% даже в странах с высоким уровнем онкологической помощи [Walther A. et al. Genetic prognostic and predictive markers in colorectal cancer // Nat. Rev. Cancer. - 2009. - V. 9. - P. 489-499]. Как и в случае большинства других онкологических заболеваний, выявление данного вида рака на ранних стадиях во многом определяет успешность лечения, ввиду чего развитие чувствительных молекулярно-генетических, в том числе эпигенетических, методов диагностики колоректального рака несомненно является актуальной задачей. Однако начальные стадии патологии внутренних органов зачастую не сопровождаются определяемыми при наружном осмотре аномалиями или ощущаемыми пациентами нарушениями функций организма, что делает раннюю диагностику в этих случаях крайне затруднительной. Таким образом, использование биохимических и генетических анализов для диагностики колоректального рака является предпочтительным способом выявления заболевания, тем более, что такая диагностика возможна еще на бессимптомной стадии [Gyparaki MT, Basdra EK, Papavassiliou AG. DNA methylation biomarkers as diagnostic and prognostic tools in colorectal cancer // J Mol Med (Berl). - 2013. - V. 91. - P. 1249-1256].

Для разработки соответствующей диагностической системы необходимо найти ответ на второй вопрос: метилирование регуляторных участков каких именно генов ассоциировано с развитием колоректального рака?

Одна из первых масштабных работ по поиску маркерных генов для эпигенетической диагностики колоректального рака была проведена в 2008 году сотрудниками компании Epigenomics AG (Германия) [Lofton-Day C et al. DNA methylation biomarkers for blood-based colorectal cancer screening // Clin Chem. - 2008. - V. 54. - P. 414-423]. На первом этапе работы были использованы метилчувствительные ферменты рестрикции для анализа более 600 кандидатных фрагментов геномной ДНК, выделенной из опухолей и нормальной ткани. Были выбраны 56 участков ДНК, которые достоверно различались по степени метилирования в здоровых и малигнантных клетках. Далее эти участки дополнительно проверялись на микрочипах и в ПЦР в реальном времени. Это позволило сократить число генов, наиболее пригодных для использования в качестве эпигенетических маркеров колоректального рака, до семи (SEPT9, NGFR, TMEFF2, ALX4, EYA4, FOXL2, SIX6).

В дальнейшем большое число авторов опубликовало результаты своих исследований по поиску генов-онкомаркеров для эпигенетической диагностики колоректального рака, проведенных различными методами, включая такие современные подходы, как применение ДНК-микрочипов (например, серии Infinium HumanMethylation, Illumina, США), ограниченное (RRBS, Reduced Representation Bisulfite Sequencing) и полногеномное бисульфитное секвенирование (WGBS, Whole Genome Bisulphite Sequencing) секвенирующими системами нового поколения (Next Generation Sequencing, NGS) [Reinders J. et al. // Epigenomics. - 2010. - V. 2. - P. 209-220; Kaneda A., Yagi K. // Cancer Sci. - 2011. - V. 102. - P. 18-24; Kibriya M.G. et al. // BMC Med Genomics. - 2011. - 4:50; Kim Y.H. et al. // Ann Surg Oncol. - 2011. - V. 18. - P. 2338-2347; Gu H. et al. // Nat Protoc. - 2011. - V. 6. - P. 468-481; Yi J.M. et al. // Tumour Biol. - 2012. - V. 33. - P. 363-372; Hinoue T. et al. // Genome Res. - 2012. - V. 22. - P. 271-282; Lange C.P. et al. // PLoS One. - 2012. - 7:e50266; Li H. et al. // Oncol Rep. - 2012. - V. 28. - P. 99-104; Roperch J.P. et al. // BMC Cancer. - 2013. - P. 13-566; Naumov V.A. // Epigenetics. - 2013. - V. 8. - P. 921-934; Chen Y.A. et al. // Epigenetics. - 2013. - V. 8. - P. 203-209; Harper K.N. et al. // Cancer Epidemiol Biomarkers Prev. - 2013. - V. 22. - P. 1052-1060; Mitchell S.M. et al. // BMC Cancer. - 2014. - P. 14-54; Melson J. et al. // Int J Cancer. - 2014. - V. 134. - P. 2656-2662; Papadia C. et al. // Inflamm Bowel Dis. - 2014. - V. 20. - P. 271-277].

Некоторые из обнаруженных генов-онкомаркеров уже применяются в коммерческих тест-системах для диагностики колоректального рака. В 2008-2012 годах в компании Epigenomics AG (Германия) была разработана и запатентована тест-система «Epi proColon kit 2.0» для диагностики колоректального рака (международная заявка №WO2005001142, опубл. 06.01.2005 г.; заявка США №2014179563, опубл. 26.06.2014 г.) на основе выявления специфического эпигенетического маркера SEPT9 в периферической крови [Heichman KA. Blood-based testing for colorectal cancer screening // Mol Diagn Ther. 2014. - V. 18. - P. 127-135]. Метилирование другого гена (VIM, виментин) в малигнантных клетках послужило основой для создания тест-системы «ColoGuard assay» от компании LabCorp (США) [Lao VV, Grady WM. Epigenetics and colorectal cancer // Nat Rev Gastroenterol Hepatol. - 2011. - V. 8. - P. 686-700]. Однако подобные моно-тест-системы не предназначены для постановки окончательного диагноза, и их показатели могут служить лишь дополнительными индикаторами развития или отсутствия заболевания.

В совокупности полученные исследователями данные свидетельствуют о ненадежности диагностики, построенной на анализе одного гена и, соответственно, об актуальности разработки и внедрения тест-систем, основанных на определении статуса метилирования ДНК регуляторных областей ряда генов-онкомаркеров, позволяющих проводить диагностику более надежно и точно. Известна работа Lind с соавт., создавших диагностический набор для анализа метилирования шести генов: CNRIP1, FBN1, INA, MAL, SNCA и SPG20 [Lind G.E. et al. Identification of an epigenetic biomarker panel with high sensitivity and specificity for colorectal cancer and adenomas // Mol Cancer 2011. - V 10. - P. 85].

Наиболее близким аналогом (прототипом) является способ определения диагноза рака у индивидуума с подозрением на колоректальный рак (Международная заявка №WO2014062218, МПК C12Q1/68, опубл. 24.04.2014 г.), включающий: получение образца от индивидуума с подозрением на рак; гидролиз высокоочищенной геномной ДНК из исследуемого биоматериала с использованием бисульфитной конверсии ДНК; определения наличия или отсутствия высокого уровня экспрессии генов (определение уровня метилирования регуляторных областей генов-онкомаркеров с использованием цифровой ПЦР) по отношению к нормальному базовому стандарту для одной диагностической панели, включающей следующие биомаркеры: TFPI2, THBD, C9ORF50, ADHFEl, FGF12, PTPRT, ZNF568, KIAA1026, SFMBT2, и / или AUTS2; и диагностики рака при наличии высокого уровня экспрессии по отношению к нормальному базовому стандарту по меньшей мере одного биомаркера.

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

Техническим результатом предлагаемого изобретения является упрощение и повышение точности способа определения метилирования сайтов RCGY регуляторных областей генов-онкомаркеров колоректального рака и расширение его функциональных возможностей.

Указанный технический результат достигается разработанным набором олигонуклеотидных праймеров и флуоресцентно-меченых зондов для идентификации злокачественного поражения толстого кишечника и прямой кишки посредством анализа в образцах ДНК метилирования сайтов RCGY в заданных положениях регуляторных областей генов-онкомаркеров методом GLAD-ПЦР, имеющих структуру, приведенную в приложении (список последовательностей) и таблице 1:

Таблица 1

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера ADHFE1 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера ALX4 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера CNRIP1 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера EID3 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера ELMO1 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера ESR1 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера FBN1(1) (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера FBN1(2) (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера FBN1(3.1) (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера FBN1(3.3) (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера HLTF (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера LAMA1 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера NEUROG1 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера NGFR (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера RARB (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера RXRG (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера RYR2 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера SDC2 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера SEPT9 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера SFRP2 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера SOCS3 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера SOX17 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера THBD (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера TMEFF2 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера UCHL1 (5'→3'):

- для специфических праймеров и зонда на участки CpG-островков регуляторных областей гена-онкомаркера VIM (5'→3'):

где g - геномный праймер (genome); р - флуоресцентно-меченый зонд (zond); h - гибридный праймер (hybrid); FAM - флуоресцеин; BHQ1 - гаситель флуоресценции (Black Hole Quencher 1).

Учет результатов анализа проводится методом гибридизационно-флуоресцентной детекции продуктов ПЦР-амплификации в режиме реального времени.

Указанный технический результат достигается также созданием способа идентификации злокачественного поражения толстого кишечника и прямой кишки посредством анализа в образцах ДНК метилирования сайтов RCGY в заданных положениях регуляторных областей генов-онкомаркеров методом GLAD-ПЦР, включающем формирование панели генов-онкомаркеров колоректального рака: ADHFE1, ALX4, CNRIP1, EID3, ELMO1, ESR1, FBN1, HLTF, LAMA1, NEUROG1, NGFR, RARB, RXRG, RYR2, SDC2, SEPT9, SFRP2, SOCS3, SOX17, THBD, TMEFF2, UCHL1 и VIM путем отбора генов-онкомаркеров колоректального рака при биоинформационном анализе и рассмотрения нуклеотидных последовательностей предположительных регуляторных участков этих генов, выявление располагающихся в них CpG-островков, анализ первичной структуры CpG-островков и их соответствие определенным критериям выбора. В регуляторной области (промоторе или первом экзоне) этих генов обычно присутствуют несколько сайтов RCGY, метилирование которых может вызывать подавление экспрессии гена при колоректальном раке. Поэтому далее на образцах ДНК из опухолевого материала или клеточных линий осуществляют отбор наиболее метилированных сайтов с помощью GLAD-ПЦР-анализа, которые в могут использоваться для диагностики и/или прогнозирования течения колоректального рака.

Далее осуществляют гидролиз высокоочищенной геномной ДНК из исследуемого биоматериала метилзависимой сайт-специфической ДНК-эндонуклеазой GlaI, лигирование универсального олигонуклеотидного адаптера , где р - фосфат с последующей амплификацией в реальном времени с использованием геномных праймеров и зондов, комплементарных последовательностям регуляторных областей указанных генов в исследуемой ДНК, и гибридных праймеров, 3'-концы которых комплементарны 5'-концам ДНК в заданных местах гидролиза GlaI, а оставшаяся часть комплементарнаадаптерной последовательности, и делают заключение о наличии последовательности R(5mC)GY, где 5mC - 5-метилцитозин, R-А или G, Y-Т или С в регуляторных областях указанных генов по появлению флуоресцентного сигнала для определения статуса метилирования сайтов в клинических образцах методом GLAD-ПЦР-анализа и на основе информации о статусе метилирования одного или нескольких сайтов делают вывод о наличии или отсутствия колоректального рака.

Структура адаптера подбирается таким образом, чтобы его нуклеотидная последовательность была уникальна для генома человека, оптимальная по GC составу и защищена от образования самокомплементарных структур.

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

Изобретение иллюстрируется следующим графическим материалом, представленным на фигурах 1-3. На фиг. 1 представлен пример поиска наиболее часто метилированного сайта RCGY на выбранном участке регуляторной области гена FBN1(3.3). Как видно из фиг. 1, таким оказался третий сайт в ампликоне, которому соответствует гибридный праймер .

На панели А представлен фрагмент нуклеотидной последовательности; показаны участки связывания геномного праймера, зонда и геном-специфичных 3'-тетрануклеотидов гибридных праймеров. На панели Б приведены кривые накопления флуоресценции в ходе ПЦР с использованием различных гибридных праймеров. Заглавными буквами «GCGG» обозначен выбранный гибридный праймер, соответствующий точке гидролиза ДНК по наиболее часто метилированному третьему сайту RCGY в ампликоне.

На фиг. 2: представлен результат GLAD-ПЦР-анализа отобранного сайта RCGY в регуляторной области гена SEPT9 в клинических образцах. Погоризонтальной оси отображены номера клинических образцов, где T – соответствует опухолевому образцу, а N – здоровому. По вертикальной оси приведено усредненное значение Cq для каждого образца с указанием доверительного интервала.

На фиг. 3 приведена ROC (Receiver Operating Characteristic) кривая, полученная в результате ROC анализа результатов GLAD-ПЦР-анализа статуса метилирования RCGY сайта-маркера колоректального рака гена SEPT9. Параметры ROC-анализа приведены в таблице 4. ROC-анализ является стандартным методом оценки диагностической ценности тестов.

Сущность изобретения поясняется следующими примерами 1-3, реализующими способ определения метилирования сайтов RCGY регуляторных областей генов-онкомаркеров колоректального рака методом GLAD-ПЦР-анализа, а также состав олигонуклеотидных праймеров и флуоресцентно-меченых зондов для осуществления указанного способа.

Пример 1. Разработка специфических праймеров и флуоресцентно-меченых зондов на участки CpG-островков регуляторных областей генов-онкомаркеров колоректального рака.

Изучение опубликованных в литературе данных по метилированию генов при онкопатологиях нижних отделов кишечника, данных RRBS-секвенирования (Reduced Representation Bisulfite Sequencing, бисульфитное секвенирование ограниченных выборок локусов) в проекте ENCODE/HudsonAlpha и последующий биоинформационный анализ нуклеотидных последовательностей регуляторных участков описанных ранее генов позволил выбрать кандидатную комбинацию участков CpG-островков регуляторных областей генов-онкомаркеров колоректального рака, удовлетворяющих критериям анализа методом GLAD-ПЦР (табл. 2).

На первом этапе работы был проведен анализ литературных данных последних лет по метилированию генов при онкопатологиях нижних отделов кишечника, данных RRBS-секвенирования (Reduced Representation Bisulfite Sequencing, бисульфитное секвенирование ограниченных выборок локусов) в проекте ENCODE/HudsonAlpha и сформирован список из 57 генов-онкомаркеров, метилирование которых происходит с высокой частотой в ДНК клеток колоректального рака. Этот список включает в себя следующие гены: SEPT9, SYNE1, FGF12, GAD2, HAND1, CACNA1G, CNRIP1, TLX2, LOX, LAMA1, GAS7, SLC6A15, BOLL, TFPI2, SOCS3, TRH, MDFI, C1orf70, RUNX3, TMEFF2, SFRP2, WNT5A, IGFBP7, TFPI2, FOXE1, NPY, WIF1, CIDEB, RIC3, KCNQ5, COL4A1, SPG20, SND1, PENK, CDO1, THBD, EDIL3, DUSP26, UCHL1, STOX2, ELMO1, SLC30A10, IGF2, ESR1, EYA4, NALCN, FLI1, RXRg, RYR2, PRKAR1B, BEND5, ITGA4, ADHFE1, NGFR, GATA4, NEYROG1, FBN1. При выборе генов для данного списка учитывалось также количество проанализированных проб в описанных авторами статей экспериментах, так как достоверные данные могут быть получены лишь при большом числе образцов ДНК в выборке. Большое значение также придавалось такому описанному авторами работ свойству, как возможность выявления метилирования гена на ранних стадиях заболевания, хотя, к сожалению, число публикаций, в которых рассматривалась зависимость метилирования гена от стадии онкопатологии, очень незначительно.

На втором этапе работы были рассмотрены нуклеотидные последовательности предположительных регуляторных участков этих 57 генов и выявлены CpG-островки, располагающиеся в них. Первичные структуры CpG-островков были проанализированы на их соответствие критериям выбора:

1) Наличие данных из опубликованных работ, показывающих высокую частоту метилирования регуляторного участка гена в исследованных выборках образцов ДНК от больных колоректальным раком при низкой частоте метилирования данного участка у здоровых доноров или в прилегающих нормальных тканях толстой и прямой кишки онкопациентов.

2) Присутствие СpG-островка в области начала гена (или в районе точки старта кодирования альтернативной изоформы белка - продукта гена), что в большинстве случаев служит свидетельством регуляции активности данного гена путем метилирования этой области.

3) Присутствие потенциального сайта узнавания эндонуклеазы GlaI (RCGY) в нуклеотидной последовательности на выбранном участке.

4) Достаточное расстояние (50 п.н. и более) между двумя потенциальными сайтами узнавания эндонуклеазы GlaI для одновременного отжига на этом участке флуоресцентно-меченого зонда и геномного праймера и возможность подбора для этого участка праймера и зонда, не подверженных теоретической кросс-гибридизации между собой и с последовательностью адаптера.

5) Значение G+C состава на кандидатном участке, не превышающее 80%, что вызвано ограничениями системы амплификации ДНК in vitro, в которой используется HotStart Taq ДНК-полимераза.

6) Уникальность нуклеотидной последовательности кандидатного участка анализа (невхождение его в число повторяющихся элементов генома) во избежание неспецифического отжига подобранного геномного праймера в других областях генома, что может привести к значительному снижению эффективности амплификации и затруднению детекции результатов при проведении ПЦР в реальном времени.

В результате проведенного биоинформационного анализа регуляторных участков каждого из перечисленных генов для проведения экспериментального исследования были отобраны 26 областей из данных участков (табл.2), которые удовлетворяют всем требуемым критериям проведения анализа с помощью метода GLAD-ПЦР, и подобраны специфические праймеры и зонды на выбранные участки регуляторных областей генов-онкомаркеров (табл.3).

Таблица 2.

Участки регуляторных областей генов-онкомаркеров перспективные для GLAD-ПЦР-анализа.

Ген (регион)	Название гена	Расположение в хромосоме	Координаты участка в хромосоме (версия генома человека GRCh38.p7)
ADHFE1	alcohol dehydrogenase, iron containing 1	8q12.3	chr8:66432540-66432621
ALX4	ALX homeobox 4	11p11.2	chr11:44304627-44304798
CNRIP1	cannabinoid receptor interacting protein 1	2p13	chr2:68319340-68319419
EID3	EP300 interacting inhibitor of differentiation 3	12q23.3	chr12:104303544-104303697
ELMO1	engulfment and cell motility 1	7p14.1	chr7:37448622-37448735
ESR1	estrogen receptor 1	6q24-q27	chr6:151807784-151807876
FBN1(1)	fibrillin 1	15q21.1	chr15:48644977-48645057
FBN1(2)			chr15:48645128-48645191
FBN1(3.1)			chr15:48645489-48645638
FBN1(3.3)			chr15:48645378-48645537
HLTF	helicase-like transcription factor	3q25.1-q26.1	chr3:149086693-149086780
LAMA1	laminin subunit alpha 1	18p11.3	chr18:7117909-7118044
NEUROG1	neurogenin 1	5q23-q31	chr5:135536157-135536283
NGFR	nerve growth factor receptor	17q21-q22	chr17:49495262-49495389
RARB	retinoic acid receptor beta	3p24	chr3:25428290-25428491
RXRG	retinoid X receptor gamma	1q22-q23	chr1:165445168-165445307
RYR2	ryanodine receptor 2	1q43	chr1:237043053-237043209
SDC2	syndecan 2	8q22-q23	chr8:96494024-96494159
SEPT9	septin 9	17q25.3	chr17:77372872-77372956
SFRP2	secreted frizzled-related protein 2	4q31.3	chr4:153789082-153789314