Статья: Влияние аберрантной экспрессии микроРНК на эффективность лучевой терапии опухолей прямой кишки

Министерство здравоохранения Российской Федерации

ФГБУ «Национальный медицинский исследовательский центр онкологии»

Влияние аберрантной экспрессии микроРНК на эффективность лучевой терапии опухолей прямой кишки

Кутилин Д.С., Гусарева М.А.,

Кошелева Н.Г., Габричидзе П.Н.,

Донцов В.А., Легостаев В.М.,

Шляхова О.В., Лиман Н.А.,

Солнцева А.А., Васильева Е.О.

Аннотация

Лучевая терапия - один из главных методов лечения опухолей прямой кишки. В современной радиотерапевтической практике отсутствие положительной реакции у больных на предоперационное облучение достаточно распространено, что обусловлено формированием радиорезистентности опухолевых клеток, ассоциированной с аномалиями на генетическом и эпигенетическом уровне. Так, ранее нами была выявлена связь транскрипционной активности генетических локусов BRCA2, H2AX, RBBP8, CASP9 и BCL2 с эффективностью лучевой терапии у больных раком прямой кишки. Поэтому целью исследования стала идентификация с помощью биоинформационного анализа микроРНК, таргетирующих гены BRCA2, H2AX, RBBP8, CASP9 и BCL2, и изучение влияния изменения экспрессии этих микроРНК на эффективность лучевой терапии опухолей прямой кишки. Биоинформационный поиск микроРНК, таргетирующих гены BRCA2, H2AX, RBBP8, CASP9 и BCL2, осуществляли с использованием модифицированного алгоритма TarPmiR и базы данных mirDB. Относительную экспрессию микроРНК оценивали методом RT-qPCR. Проведенный биоинформационный анализ позволил установить, что гены BRCA2, H2AX, RBBP8, CASP9 и BCL2 таргетируются 1338 микроРНК, из которых валидированы в базе данных miRDB и обладают минимальной свободной энергией взаимодействия с мРНК 26 микроРНК (9 для BCL2, 2 для BRCA2, 12 для CASP9, 2 для H2AX и 1 для RBBP8). Профилирование экспрессии 26 микроРНК в нормальных и опухолевых тканях у больных раком прямой кишки с полным регрессом опухоли (группа 1, n=17) и с незначительным регрессом опухоли или отсутствием динамики (группа 2, n=15) позволило выявить 13 дифференциально экспрессирующихся микроРНК (miR-1249-5p, miR-6874-5p, miR-4728-5p, miR-6808-5p, miR-3202, miR- 5195-3p, miR-6820-3p, miR-557, miR-6757-3p, miR-1273h-5p,miR-6737-5p, miR-130b и miR-6819-5p), регулирующих системы апоптоза и репарации ДНК, и ассоциированных с эффективность/неэффективностью лучевой терапии опухолей прямой кишки.

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

Annotation

Influence of aberrant microRNA expression on radiotherapy efficiency for rectal tumors

Kutilin D.S., Gusareva M.A., Kosheleva N.G., Gabrichidze P.N., Dontsov V.A., Legostaev V.M., Shlyakhova O.V., Liman N.A., Solntseva A.A., Vasileva E.O., National medical research oncology center, Russian Federation

Radiation therapy is one of the main methods in the treatment of rectal cancer, however, in oncological practice, there are often cases of non-response patients to preoperative radiation therapy, which is associated with the radioresistance of tumor cells, which depends on their genetic and epigenetic characteristics. Previously, it was shown that the expression of BRCA2, H2AX, CASP9, RBBP8 and BCL2 genes affects the effectiveness of radiation therapy for rectal tumors. Therefore, the aim of the study was to identify, using bioinformatic analysis, microRNAs targeting genes BRCA2, H2AX, CASP9, RBBP8 and BCL2, and to study the effect of changes in the expression of these microRNAs on the effectiveness of rectal tumors radiation therapy. Bioinformatic search for microRNAs targeting genes BRCA2, H2AX, CASP9, RBBP8 and BCL2 was carried out using a modified TarPmiR algorithm and the mirDB database. The relative expression of microRNA was assessed by the RT- qPCR method. The bioinformatic analysis performed made it possible to establish that genes BRCA2, H2AX, CASP9, RBBP8 and BCL2 targeted by 1338 microRNAs, 26 microRNAs of which are validated in the miRDB database and have the minimum free energy of interaction with mRNA (9 for BCL2, 2 for BRCA2, 12 for CASP9 , 2 for H2AX and 1 for RBBP8). Expression profiling of 26 miRNAs in normal and tumor tissues in patients with rectal cancer with complete tumor regression (group 1, n = 17) and with insignificant tumor regression or no dynamics (group 2, n = 15) made it possible to identify 13 differentially expressed miRNAs (miR-1249-5p, miR- 6874-5p, miR-4728-5p, miR-6808-5p, miR-3202, miR-5195-3p, miR-6820-3p, miR-557, miR-6757-3p, miR-1273h- 5p,miR-6737-5p, miR-130b и miR-6819-5p), providing regulation of the DNA repair system and apoptosis, and associated with the efficacy / ineffectiveness of rectal tumors radiation therapy.

Keywords: rectal cancer, microRNA, radiation therapy, signaling pathways, apoptosis, DNA repair, radioresistance.

Постановка проблемы

Среди онкологических заболеваний по числу летальных исходов 4-е место занимают злокачественные новообразования кишечника. В данной группе наибольший процент составляют опухоли прямой кишки [1; 2]. В настоящее время для лечения рака прямой кишки применяют комбинированные протоколы, отличающиеся от протоколов для других отделов кишечника, так как сочетают предоперационную лучевую терапию (в стандартном варианте - разовая очаговая доза (РОД) 2,4 Гр до суммарной очаговой дозы (СОД) 54 Гр) и последующее хирургическое вмешательство [3; 4]. Анализ результатов отечественной и мировой радиотерапевтической практики показывает наличие большого числа (до 30-40%) случаев с радиорезистентностью опухолей, обусловленной молекулярными (генетическими/эпигенетическими) аномалиями в опухолевых клетках, в частности аберрантной экспрессией определенных групп генов. Ранее нами было показано, что транскрипционная активность генетических локусов, регулирующих системы репарации (RBBP8, H2AX, BRCA2), пролиферации и апоптоза (CASP9, BCL2), связана с эффективностью проводимой лучевой терапии у больных раком прямой кишки [1]. Как известно, важнейшими регуляторами транскрипционной активности генов являются микроРНК, представляющие собой некодирующие одноцепочечные РНК, содержащие приблизительно два десятка нуклеотидов, осуществляющие регуляцию трансляции и деградации матричной РНК посредством связывания с комплементарными сайтами в нетранслируемых участках молекул последних, служащих их мишенями [5].

Поэтому целью исследования стало выявление с помощью биоинформационного анализа микроРНК, таргетирующих гены RBBP8, H2AX, BRCA2, CASP9 и BCL2, и изучение влияния изменения экспрессии этих микроРНК на эффективность лучевой терапии опухолей прямой кишки.

Материалы и методы исследования

Биоинформационный поиск микроРНК, таргетирующих гены RBBP8, H2AX, BRCA2, CASP9 и BCL2, осуществляли с использованием модифицированного алгоритма TarPmiR и базы данных mirDB. TarPmiR для прогнозирования сайта связывания мРНК и микроРНК использует алгоритм машинного обучения random forest («случайный лес»), сочетающий в себе метод случайных подпространств и метод бэггинга Бреймана. Результатом работы алгоритма является значение вероятности того, что обнаруженный сайт-мишень кандидат и есть истинный сайт-мишень [6; 7].

В исследовании участвовали 32 пациента (возраст от 45 до 65 лет), госпитализированных в ФГБУ «НМИЦ онкологии» с диагнозом аденокарцинома прямой кишки (G1-2). До лучевой терапии (ЛТ) от этих пациентов при проведении видеоколоноскопии (ВКС) получали препараты биопсии - фрагменты немалигнизированных (нормальных) и опухолевых тканей прямой кишки. Фрагменты ткани, извлеченные во время биопсии, мгновенно замораживали при минус 195 °C и в последующем использовали для выделения суммарной РНК. ЛТ проводили по стандартной схеме - РОД 2,4 Гр, СОД 54,0 Гр - на линейном ускорителе частиц Novalis TX.

Препараты суммарной РНК выделяли методом фенол-хлороформной экстракции в модификации [7]. Для выявления микроРНК и малой РНК U6 препараты суммарной РНК подвергали реакции обратной транскрипции (ОТ), которая проводилась одновременно с полиаденилированием РНК, с использованием специфичных RT-праймеров. Олигонуклеотидные праймеры были разработаны нами с использованием алгоритма предложенного Balcells I. [8]. Реакцию ОТ проводили в одном повторе, отдельно для каждой микроРНК, для этого использовали лиофилизированную смесь «Мастер-микс ОТ» («Вектор- Бест», Россия). ОТ проводили при 16 °С (15 минут), далее при 42 °С (15 минут), обратную транскриптазу инактивировали инкубацией при 95 °С в течение 2 минут [7].

Относительную экспрессию микроРНК оценивали методом RT-qPCR. Постановку RT-qPCR каждого образца проводили в 3 повторах. Для подбора референсных локусов оценивали стабильность экспрессии микроРНК по алгоритму geNorm [9]. Первоначальный список предполагаемых референсных локусов для микро-РНК включал: miR-191; miR-23a и U6. С помощью geNorm для нормализации данных по экспрессии микроРНК был выбран U6. Нормализацию проводили по референсному локусу и уровню экспрессии соответствующих микроРНК в образцах нормальной ткани. Относительную экспрессию микроРНК (REmir) определяли по формуле REmir = E'AACt, где Е - вычисленная эффективность реакции амплификации [5].

Статистический анализ данных проводили с использованием языка программирования Python (библиотека SciPy). Различия между группами оценивались с использованием U-критерия Манна-Уитни, для корректировки множественного сравнения была применена поправка Бонферрони [7]. Также проводили анализ избыточной представленности микроРНК в сигнальных путях (ORA, Over-Representation Analysis). Этот метод определяет долю участия дифференциально экспрессирующихся микроРНК в определенных сигнальных путях, его результатом является получение списка наиболее значимых сигнальных путей. Статистическая значимость в ORA рассчитывалась с применением точного критерия Фишера [10].

радиорезистентность опухолевый кишка генетический лучевой

Результаты исследования и их обсуждение

Биоинформационный поиск микроРНК, таргетирующих гены RBBP8, H2AX, BRCA2, CASP9 и BCL2, с использованием модифицированного алгоритма TarPmiR, выявил 1338 микроРНК (рис. 1А), из которых валидированы в базе данных miRDB 86 микроРНК (рис. 1Б), в том числе 26 с минимальным значением термодинамического потенциала (свободной энергии) взаимодействия микроРНК и матричной РНК (табл. 1): 9 микроРНК, взаимодействующих с BCL2, 2 микроРНК, взаимодействующие с BRCA2, 12 микроРНК, взаимодействующих с CASP9, 2 микроРНК, взаимодействующие с H2AX, и 1 микроРНК, взаимодействующая с RBBP8.

Таблица 1

МикроРНК, таргетирующие гены RBBP8, H2AX, BRCA2, CASP9 и BCL2, выявленные с использованием алгоритма TarPmiR*

* - представлены только микроРНК, валидированные в базе данных miRDB и с минимальной свободной энергией взаимодействия микроРНК и матричной РНК.

Рис. 1. А - все микроРНК, таргетирующие гены RBBP8, H2AX, BRCA2, CASP9 и BCL2, выявленные с использованием алгоритма TarPmiR. Б - микроРНК, выявленные с использованием алгоритма TarPmiR и валидированные в базе данных miRDB

Для 26 микроРНК, выявленных с использованием биоинформационных подходов, проведено профилирование экспрессии в нормальных и опухолевых тканях у больных раком прямой кишки.

Анализ результатов ЛТ опухолей прямой кишки у 32 больных позволил разделить их на 2 группы. У 17 больных после ЛТ был зафиксирован полный регресс опухоли (группа 1), у 8 больных был зафиксирован незначительный регресс опухоли, и у 7 больных отсутствовала динамика изменения размера опухоли (группа 2). В этих двух группах пациентов наблюдалась дифференциальная экспрессия ряда микроРНК (как относительно нормальной ткани, так и относительно каждой из групп). Так, у больных с незначительным регрессом опухоли или отсутствием динамики после лучевой терапии (n = 15) выявлено статистически значимое (p<0,05) снижение экспрессии miR-1249 в 2,5 раза относительно экспрессии в нормальной ткани и в 2,3 раза относительно экспрессии в опухолевой ткани больных с полным регрессом опухоли, снижение экспрессии miR-6820 в 1,7 раза относительно экспрессии в нормальной ткани ив 1,6 раза относительно экспрессии в опухолевой ткани у больных с полным регрессом опухоли, снижение miR-4717 в 2,5 раза относительно нормальной ткани, снижение hsa-miR-3943 в 3,3 раза по сравнению с уровнем экспрессии в нормальной ткани, снижение miR-557 в 2,5 раза относительно нормальной ткани и в 2,3 раза относительно пациентов с полным регрессом опухоли, при этом наблюдается увеличение экспрессии hsa-miR-4690-5p в 2,1 раза относительно нормальной ткани (рис. 2).

Также у больных с полным регрессом опухоли была снижена экспрессия miR-4717 в 1,8 раза (р<0,05) относительно уровня экспрессии в нормальной ткани прямой кишки. Данные микроРНК (miR-557, miR-1249, miR-6820, miR-4717, miR-3943), за исключением miR-4717, в опухолевой ткани у больных с незначительным регрессом опухоли или отсутствием динамики обладают преимущественно сниженной экспрессией как относительно нормальной ткани, так и относительно опухолевой ткани пациентов с полным регрессом опухоли после лучевой терапии. Эти микроРНК таргетируют антиапоптозный ген BCL2, и их сниженная активность может способствовать сохранению на повышенном уровне транскрипционной активности этого гена-мишени. BCL2 кодирует белок, участвующий в подавлении апоптоза двумя путями: во-первых, этот белок изменяет проницаемость мембран митохондрий и предотвращает выход из них цитохрома C, тем самым останавливая активацию каспазного каскада; во-вторых, он инактивирует белок CED4 (фактор активации апоптоза) [11]. Соответствующий эффект увеличения экспрессии гена BCL2 у пациентов с незначительным регрессом опухоли или отсутствием динамики после лучевой терапии мы и обнаружили в предыдущем исследовании [1].