Экспрессия экзосомальных микроРНК при раке яичников высокой степени злокачественности и эндометриоидных кистах яичников

Юрова М.В., Эльдаров Ч.М., Бобров М.Ю., Хабас Г.Н., Павлович С.В.

1) ФГБУ «Национальный медицинский исследовательский центр акушерства, гинекологии и перинатологии имени академика В.И. Кулакова» Минздрава России, Москва, Россия; 2) ГАОУ ВО «Первый Московский государственный медицинский университет им. И.М. Сеченова» Минздрава России (Сеченовский Университет), Москва, Россия
Экзосомальные микроРНК рассматриваются в качестве перспективных неинвазивных маркеров различных пролиферативных процессов, в частности, затрагивающих органы женской репродуктивной системы.
Цель: Оценить экспрессию экзосомальных микроРНК в плазме крови пациентов с серозным раком яичников (РЯ) и пациентов с наружным генитальным эндометриозом и эндометриоидными кистами яичников (ЭКЯ) для выявления потенциальных маркеров данных заболеваний.
Материалы и методы: В исследование были включены 7 пациенток с гистологически верифицированным серозным РЯ высокой степени злокачественности и 6 женщин контрольной группы без данной патологии («Контроль-РЯ»), а также 6 пациентов с ЭКЯ и 5 женщинрепродуктивного возраста без данной патологии («Контроль-ЭКЯ»). Экзосомы выделяли из плазмы крови, состав экзосомальных микроРНК определяли методом высокопроизводительного секвенирования нового поколения.
Результаты: При попарном сравнении представленности микроРНК в группах «РЯ» — «Контроль- РЯ» и «ЭКЯ» — «Контроль-ЭКЯ» было выявлено соответственно 22 и 13 дифференциально экспрессирующихся микроРНК (ДЭМ), статистически значимо отличающихся в 2 раза и более. По данным биоинформатического анализа, были предложены потенциальные маркеры исследуемых заболеваний: для РЯ - miR-141-3p, -199a-5p, -200b-3p, -203a-3p, -224-5p, -4488; для ЭКЯ - miR-92b-5p, -486-5p, -3184-3p, -4732-5p, -423-5p. Для ДЭМ в каждой паре сравнения были определены потенциальные гены- мишени, на основании списков которых был проведен поиск путей внутриклеточной сигнализации, наиболее вероятно задействованных в патогенезе РЯ и ЭКЯ. Были идентифицированы следующие потенциальные гены-мишени нескольких сигнальных путей, регулируемых одной или несколькими ДЭМ при РЯ и ЭКЯ: AKT1, ATM, BARD1, BAX, BCL2, BRCA1, CASP3, CDK4, CHEK1, CHEK2, JAK1, MDM2, PLK1, PTEN, RB1, SMAD2, SMAD3, TP53. Функциональная кластеризация генов-мишеней показала их участие в процессах регуляции сигнальных путей факторов роста, клеточного цикла, апоптоза и репарации ДНК.
Заключение: Выявление уникальных ДЭМ в экзосомах плазмы крови свидетельствует о наличии специфичных изменений профиля микроРНК, характерного для РЯ и ЭКЯ. Перечисленные экзосомальные микроРНК могут рассматриваться в качестве кандидатов в маркеры изученных пролиферативных процессов. Многочисленные взаимосвязи между выявленными генами-мишенями микроРНК указывают на их существенный вклад и совместную вовлеченность в процессы пролиферативной активности при РЯ и ЭКЯ.

Ключевые слова

рак яичников
экзосомы
экзосомальные микроРНК
эндометриоз

Диагностика патологических пролиферативных процессов органов женской репродуктивной системы является актуальной проблемой современной гинекологии. В структуре всех онкогинекологических причин смертности у женщин злокачественные новообразования яичников занимают пятую позицию в структуре летальности, обусловленной опухолями всех локализаций у женщин [1]. Эндометриоз характеризуется формированием эндометриоидной ткани вне полости матки и поражает примерно 10-15% женщин репродуктивного возраста [2, 3], вызывая бесплодие и существенно ухудшая качество жизни пациентов [4]. Эндометриоз считается доброкачественным заболеванием; тем не менее, клетки эндометриоидных гетеротопий проявляют черты, сходные с онкотрансформированными клетками: нарушение регуляции пролиферации, апоптоза, миграции, способность стимулировать ангиогенез, резистентность к воздействию имму-нокомпетентных клеток приводят к росту и распро-странению эндометриоидных очагов [5].

Циркулирующие микроРНК представляют собой подкласс малых некодирующих РНК, которые были обнаружены практически во всех жидкостях организма, где, как предполагается, участвуют в межклеточных коммуникациях [6]. Биогенез микроРНК в клетке происходит сначала в ядре, а затем в цитоплазме, где они обеспечивают пост-транскрипционную регуляцию экспрессии целевых генов [7]. Изменение экспрессии микроРНК в опу-холевых тканях яичника и при эндометриозе яичника было продемонстрировано в ряде работ, что свидетельствует об их существенной роли в патогенезе данных заболеваний [8, 9]. МикроРНК могут секретироваться в комплексе со специфическими белками, либо в составе экзосом в межклеточное пространство или циркуляторное русло [10]. Экзосомы являются разновидностью внеклеточных везикул, содержащих специфический молекулярный репертуар, благодаря которому они способны влиять на активность клеток-мишеней [10, 11]. Изменение состава циркулирующих экзосом, в частности содержания микроРНК, может отражать активность патологического процесса, и в этой связи экзосомальные микроРНК рассматриваются в качестве перспективных неинвазивных маркеров [12]. В ряде исследований был изучен диагностический потенциал циркулирующих микроРНК при различных типах рака яичника (РЯ) [13]. Результаты оказались противоречивыми, что может быть обусловлено несколькими причинами: различиями в методах отбора и группировки пациентов, использовании сыворотки или плазмы крови, а также в методах детекции и идентификации микроРНК. Кроме этого, в подавляющем большинстве исследований анализировался общий состав микроРНК плазмы или сыворотки крови без индивидуального анализа экзосомальных фракций. Необходимо отметить, что состав циркулирующих экзосомальных микроРНК при эндометриозе остается практически не изученным. Еще одним ограничением исследований в данной области является применение в подавляющем большинстве работ метода полимеразной цепной реакции, который, обладая высокой чувствительностью в определении конкретных микроРНК, не позволяет оценить весь спектр микроРНК в исследуемых образцах. В этой связи, в данной работе был использован метод секвенирования нового поколения, позволяющий проводить наиболее полное профилирование экзо- сомальных микроРНК.

Целью проведенного исследования была оценка дифференциальной экспрессии экзосомальных микроРНК в плазме крови пациентов с распро-страненными стадиями серозного РЯ и наружного генитального эндометриоза, а также добровольцев соответствующих групп контроля для выявления потенциальных маркеров данных заболеваний.

Материалы и методы

В одномоментное (поперечное) исследование, выполненное в НМИЦ АГП им. В.И. Кулакова (Москва, Россия), были включены пациенты с серозным РЯ высокой степени злокачественности III—IV стадий по классификации FIGO (the International Federation of Gynecology and Obstetrics) (n=7) и с эндометриоидной кистой яичника (ЭКЯ, n=6) при распространенном наружном генитальном эндометриозе III-IV стадий по классификации rASRM (the revised American Society for Reproductive Medicine). Диагноз был верифицирован гистологически. Для образцов крови пациентов исследуемых групп были подобраны образцы крови женщин сопоставимого постменопаузального (n=6, 65 (10) лет в группе контроля и 63 (5) года в группе пациентов с РЯ, р=0,59) и репродуктивного (n=5, 34 (2) года в группе контроля и 31 (1) год у пациентов с ЭКЯ, р=0,06) возраста соответственно.

Критерии включения: информированное добровольное согласие на участие в исследовании; данные о сопутствующей экстрагенитальной патологии (в том числе острых инфекционно-воспалительных и аутоиммунных процессах), доброкачественных и злокачественных заболеваниях любой локализации на момент включения в исследование и перенесенных ранее; фолликулярная (пролиферативная) фаза менструального цикла для женщин репродуктивного возраста; постменопауза для пациентов с РЯ и женщин, составивших их группу контроля; отсутствие приема гормональных препаратов. Отсутствие заболеваний репродуктивных органов и экстрагенитальной патологии у наблюдаемых групп контроля исключено на основании данных анамнеза, ультразвукового исследования органов малого таза, клинического анализа крови и онкомаркеров СА125 и НЕ4).

Образцы крови пациентов в объеме 10 мл были получены до оперативного вмешательства. Транспортировка крови в лабораторию осуществ-лялась в пробирках с этилендиаминтетрауксусной кислотой в течение 30 минут от момента получения материала. Кровь центрифугировалась в течение 20 минут при 300 g, при 4°C. После центрифугирования надосадочная жидкость отбиралась и повторно центрифугировалась 15 минут при 3000 g. Полученную плазму вносили в криопробирки порциями по 1,5—2 мл. Промаркированные образцы замораживались и хранились при температуре -80°С. Пробоподготовка плазмы и хранение образцов осуществлялись в условиях Лаборатории по сбору и хранению биоматериалов НМИЦ АГП им. В.И. Кулакова.

Для поиска потенциальных маркеров РЯ и ЭКЯ было произведено исследование представленности микроРНК в экзосомах различных групп пациентов методом высокопроизводительного секвенирования. Перед выделением РНК, плазму крови размораживали и центрифугировали 16000 g при 4°С. Экзосомальные РНК выделяли из 3 мл полученной после центрифугирования плазмы крови при помощи набора реагентов ExoRNeasy Maxi kit (Qiagen), и далее использовали для создания библиотек ком-плементарных ДНК малых РНК при помощи набора NEBNext Multiplex Small RNA Library Prep Set for Illumina, строго в соответствии с методикой про-изводителя. Анализ микроРНК проводили методом высокопроизводительного секвенирования на платформе Illumina NextSeq, с использованием набора реагентов Illumina NextSeq 500/500 High Output Kit ver. 2,5 (75 cycles).

В результате секвенирования были получены прочтения нуклеотидных последовательностей для образцов, которые были картированы на базу данных miRbase v21 (miRbase.org) для идентификации зрелых микроРНК. Для нормализации данных и анализа дифференциальной экспрессии микроРНК в группах использовался алгоритм DESeq2 [14]. Для дальнейшего исследования были отобраны микроРНК с количеством прочтений >10, учитывалась также кратность изменений >1,5. При отборе микроРНК статистически значимыми считались результаты при рadj<0,05, где рadj — p-значение с учетом поправки средней доли ложных отклонений гипотез при проведении множественных сравнений (false discovery rate, FDR). Для отобранных микроРНК был выполнен поиск возможных мишеней в базе данных MicroRNA-Target Interactions (miRTarBase [15]), содержащей информацию более, чем о пятидесяти тысячах взаимодействий микроРНК с мишенями, данные о которых регулярно обновляются согласно систематическому анализу релевантных исследований. Были отобраны только те взаимодействия «микроРНК - ген», которые были выявлены одним или несколькими релевантными («сильными», по классификации авторов базы данных) методами: полимеразная цепная реакция в реальном времени (qPCR), Вестерн-блот (WesternBlot), методы высокопроизводительного репортерного анализа (reporter assay kit). Полученные списки потенциальных мишеней использовались для поиска сигнальных путей, регулируемых этими мишенями и, как следствие, изменяющих свою активность. Использовалась модель известных и прогнозируемых белок-белковых взаимодействий (база https://string-db.org/) и проводился анализ представленности генов- мишеней в различных сигнальных путях из базы данных WikiPatways (https://www.wikipathways.org). В результате «обогащения» сигнальных путей для каждой пары сравниваемых групп были отобраны те из них, которые продемонстрировали наибольшую статистическую значимость (XD-score, q-value), определяемую с помощью точного теста Фишера перекрывающихся с заданным списком мишеней. При выборе сигнальных путей использовалась поправка на тканеспецифичность - были отобраны пути с наиболее высокими показателями представленности, свидетельствующими об активности при патологических процессах в яичниках.

Результаты

По данным секвенирования в экзосомах было выявлено 480 микроРНК с различной представлен-ностью в группах сравнения. Списки дифференциально экспрессирующихся микроРНК (ДЭМ), а также кратность и направление изменений представлены в таблице 1 и на рисунке 1.

71-1.jpg (68 KB)

Среди выявленных ДЭМ, дифференцирующих пациентов от добровольцев репродуктивного возраста группы контроля (рис. 1А), 22 демонстрировали кратность изменений в 2 раза и более при ЭКЯ: 14 ДЭМ - с повышенной и 8 ДЭМ - с пониженной экспрессией. При РЯ были получены статистически значимые различия экспрессии 13 микроРНК (повышение - 12 ДЭМ, снижение - 1 ДЭМ), которые дифференцировали пациентов от добровольцев постменопаузального возраста группы контроля (рис. 1Б); наибольшая кратность изменений наблюдалась для miR-200b-3p, miR-4488, miR-141-3p (увеличение в 21,5, 16,7 и 11,7 раз, соответственно). Также была проведена оценка различий представленности микроРНК в экзосомах крови больных РЯ и ЭКЯ. Было выявлено 27 ДЭМ, из которых 18 - с повышенной и 9 - с пониженной экспрессией в 2 и более раз, соответственно (табл. 1, рис. 1В).

72-1.jpg (430 KB)

В связи с тем, что исследуемые группы представлены пациентами разных возрастов, для оценки возможного влияния возрастного фактора также было проведено сравнение между контрольными группами. Были установлены достоверные различия представленности 17 экзосомальных микроРНК у женщин групп контроля репродуктивного и пост-менопаузального возрастов: в постменопаузе выявлено повышение 5 ДЭМ и снижение 12 ДЭМ. Однако пересечение списков ДЭМ групп сравнения «РЯ - ЭКЯ» и «Контроль-РЯ» - «Контроль-ЭКЯ» показало присутствие всего двух общих ДЭМ (miR- 484 и miR-4508); кроме того, их экспрессия при РЯ и ЭКЯ была снижена по сравнению с женщинами групп контроля. Это свидетельствует об отсутствии влияния возрастного фактора на представленность остальных микроРНК в группе сравнения «РЯ - ЭКЯ» (рис. 2).

73-1.jpg (162 KB)

Пересечение списков ДЭМ, полученных при сравнениях «ЭКЯ» - «Контроль-ЭКЯ», «РЯ» - «Контроль-РЯ», не выявило общих микроРНК (рис. 2), что свидетельствует об уникальности набора микроРНК для каждой патологии. Списки микроРНК по каждому пересечению представлены в таблице 1. Также в результате данного пересечения для каждой группы были установлены уникальные ДЭМ с кратностью изменений в 2 раза и более: 18 микроРНК - при ЭКЯ, 2 микроРНК - при РЯ и 26 микроРНК - при «РЯ - ЭКЯ» (выделены курсивом в табл. 1). Наличие уникальных ДЭМ в экзосомах плазмы крови может свидетельствовать о специфическом изменении их состава в зависимости от исследуемой патологии. Отобранные в качестве потенциальных маркеров ДЭМ представлены в таблице 2.

Кроме уникальных ДЭМ, в представленные списки были добавлены ДЭМ, которые встречались в парах сравнения «РЯ - ЭКЯ» и «Контроль-РЯ» - «Контроль-ЭКЯ». Так, для РЯ маркерные ДЭМ miR-199a-5p и miR-224-5p выявлялись также при сравнении контрольных групп разного возраста. С учетом распространенности РЯ преимущественно в постменопаузальном возрасте, использование данных ДЭМ в качестве маркеров не будет приводить к ложноположительным результатам. Тем не менее, возрастной фактор необходимо учитывать при их потенциальном использовании.

Полученные в ходе исследования списки ДЭМ были использованы для поиска потенциальных генов-мишеней в базе данных MirTarBase v. 7.0. Поскольку микроРНК теоретически могут реализовывать взаимодействия с большим количеством матричных РНК (мРНК), по результатам поиска отбирались только наиболее вероятные пары «микроРНК - мРНК», взаимодействие которых было установлено с высокой степенью валидации (подтвержденные одновременно методами qPCR, WesternBlot, reporter assay kit и др. по данным MirTarBase). Полученные таким образом набо-ры генов-мишеней были использованы для анализа представленности функциональных групп генов (gene set enrichment analysis) по базе данных WikiPathways. На основании статистических параметров значимости и достоверности выбора наборов генов-мишеней для того или иного процесса (XD, q-value [16]) и их представленности для каждой группы сравнения были отобраны пути, из которых были выделены списки наиболее представленных в двух группах сравнения. Параметр XD указывает, насколько значимо изучаемый список генов интегрируется в тот или иной процесс, реализуемый определенной сетью молекулярных взаимодействий. Параметр q представляет собой скорректированное значение достоверности (p) для параметра XD. Отобранные пути для групп сравнения представлены в таблице 3.

74-1.jpg (289 KB)

Среди клеточных процессов и сигнальных каскадов особенного внимания заслуживают пути, акти-вируемые эстрогенами и андрогенами, факторами роста и морфогенными регуляторами. Также была выявлена группа путей, контролируемых интерлей-кинами (IL-2, -3, -4, -5, -6, -7, -9), которые могут участвовать в регуляции воспалительных процессов при формировании очагов патологического процесса при раке и эндометриозе. Также были установлены внутриклеточные процессы, участвующие в регуляции репарации ДНК при ее повреждении. Обращает на себя внимание обогащение с высоким значением параметра XD по регуляторным путям, ассоциированным с повреждением ДНК, реализующимся при участии микроРНК (miRNA involved in DDR). Активация указанных путей, а также сигнализация, опосредованная интерлейкинами, TP53 и факторами роста может свидетельствовать о развитии воспалительного процесса и окислительного стресса в местах формирования исследованных экзосом. Повреждение ДНК, в свою очередь, может быть фактором, индуцирующим процессы апоптоза, аутофагии и клеточного старения, которые также были выявлены в результате обогащения списков генов-мишеней ДЭМ.

Также по генам-мишеням микроРНК в обеих парах сравнения достоверно выявлялся интегральный опухолевый путь (Intedrated Cancer Pathway). В обеих группах были выявлены следующие гены- мишени, участвующие в данном пути: AKT1, ATM, BARD1, BAX, BCL2, BRCA1, CASP3, CDK4, CHEK1, CHEK2, JAK1, MDM2, PLK1, PTEN, RB1, SMAD2, SMAD3, TP53. Анализ функциональных взаимодействий показал многочисленные взаимосвязи между выявленными генами-мишенями микроРНК, что указывает на их существенный вклад и совместную вовлеченность в процессы опухолевого роста. Функциональная кластеризация генов-мишеней показала их участие в процессах регуляции клеточного цикла, сигнальных путей репарации ДНК, путях, активируемых эстрогенами и андрогенами, факторами роста и интерлейкинами. На рисунке 3 показана кластеризация и взаимосвязь генов- мишеней, а также их участие в клеточных процессах и путях внутриклеточной сигнализации.

75-1.jpg (325 KB)

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

Обсуждение

Изучение экзосом является перспективным направлением молекулярной медицины в связи с входящими в их состав нуклеиновыми кислотами, белками, гликоконъюгатами, липидами и даже органеллами — данные компоненты динамично отображают патофизиологические условия [17], участвуют в межклеточных коммуникациях [10, 11] и при этом защищены от внешнего ферментативного воздействия липидным бислоем везикулярной мембраны [18]. В представленном исследовании был проведен транскриптомный анализ экспрессии циркулирующих экзосомальных микроРНК плазмы крови у больных с серозным РЯ III—IV стадий и пациентов с распространенным наружным генитальным эндометриозом III—IV стадий и ЭКЯ. Было показано, что перечень микроРНК изменен в крови пациентов с РЯ и ЭКЯ, по сравнению с образцами крови наблюдаемых групп контроля.

Количество прочтений при секвенировании является параметром относительной представленности микроРНК в исследуемом образце. Слишком низкие значения количества прочтений существенно снижают вероятность детекции и использования в сравнительном анализе кандидатных микроРНК. Поэтому в качестве потенциальных маркеров необ-ходимо отбирать ДЭМ с максимально возможными количеством прочтений и кратностью изменений. На основании данных критериев был проведен отбор 11 микроРНК (табл. 2).

Изменение представленности микроРНК может быть обусловлено не только конкретным заболеванием, но и гормональными изменениями, характерными для фаз менструального цикла [19, 20] и периодов активности репродуктивной системы [21] у пациентов разных возрастов. Однако данные остаются противоречивыми. В исследовании S. Moustafa et al. определение перечня микроРНК (miR-125b-5p, miR-150-5p, miR-342-3p, miR-451a, miR-3613-5p, let-7b) позволило высокоточно дифференцировать образцы пациентов с легкими и тяжелыми формами эндометриоза (AUC=0,94) вне зависимости от фазы менструального цикла и приема гормональных препаратов (комбинированные оральные контрацептивы, прогестерон, и др.) [22]. N. Zafari et al. также были получены данные об отсутствии влияния фазы менструального цикла на экспрессию miR-199b- 3p и miR-224-5p, тогда как изменения let-7d-3p у пациентов с эндометриозом были взаимосвязаны с фазой цикла, а у наблюдаемых группы контроля — нет [20]. Для проверки гипотезы о конфаундерном влиянии возраста нами было проведено сравнение экспрессии микроРНК у здоровых добровольцев репродуктивного и постменопаузального возрастов. Несмотря на то, что были получены статистически значимые различия экспрессии 26 экзосомальных предположительно «возраст-ассоциированных», микроРНК, сравнение их экспрессии у здоровых добровольцев и пациентов с пролиферативными процессами не обнаружило совпадений, что свидетельствует об отсутствии влияния возрастного фактора на представленность данных микроРНК при РЯ и ЭКЯ. Это означает, что использование данных ДЭМ в качестве маркеров не будет приводить к ложноположительным результатам. С одной стороны, это позволяет считать применение предложенных микроРНК (табл. 2) релевантным для выявления РЯ и ЭКЯ в различных возрастных группах и, в том числе, для использования в качестве маркеров при сложностях дифференциальной диагностики данных пролиферативных процессов на распространенных стадиях. Однако, с другой стороны, учитывая ограничения проведенного пилотного исследования, мы считаем, что при потенциальном использовании микроРНК возрастной фактор по- прежнему необходимо учитывать.

Полученные нами данные по некоторым маркерам сопоставлены с результатами опубликованных исследований и приведены в таблице 4.

76-1.jpg (260 KB)

Согласно данным литературы, к наиболее изученным при заболеваниях яичников относятся микроРНК семейств miR-200, let-7, miR-199 [19, 20, 37]. В работе Q. Cao et al. [38] была показана прогностическая значимость определения микроРНК в тканях серозного РЯ высокой степени злокачественности: повышение экспрессии miR-200a, miR-200b и miR-200c ассоциировано с сокращением общей выживаемости. Нами было обнаружено повышение экспрессии представителей данного семейства (miR-203a-3p, miR-200b-5p и miR-141-3p), что можно интерпретировать в качестве отображения процессов, происходящих в опухолевой ткани. Преимуществом является обнаружение изменений при помощи неинвазивного метода. В дальнейшем может быть также изучена роль данного маркера в прогнозировании выживаемости пациентов с распространенными стадиями РЯ. Что касается обнаруженного повышения miR-141-3p при РЯ, S. Masoumi-Dehghi et al. было показано, что miR-141-3p индуцирует активность рецептора 2 типа эндотелиального фактора роста (VEGFR-2), что интенсифицирует опухолевый ангиогенез [23]. Также при исследовании активности сигнального пути JAK-STAT3 в культурах клеток РЯ авторами было показано, что секретируемые опухолью экзосомальные микроРНК miR-141-3p, выступая в качестве важных посредников межклеточной коммуникации, подавляют активность цитокин-индуцируемых супрессоров (SOCS5), что приводит к стимуляции активности JAK-STAT3 (The Janus kinase/signal transducer and activator of transcription) в эндотелиальных клетках. Нами были получены данные о повышении экспрессии miR-203a-3p (padj=0,011) при распространенных стадиях серозного РЯ; аналогичные данные были опубликованы D. Taylor et al. в 2008 г. [30]. Авторы также показали эквивалентные изменения в экзо- сомах плазмы крови и непосредственно в клетках опухоли. Полученные изменения экспрессии данной микроРНК могут быть объяснены механизмами эпигенетической регуляции, в том числе гипометилированием.

miR-486-5p отобран нами, как маркер распростра-ненного эндометриоза с ЭКЯ (padj=0,025), что подтвер-ждает воспроизводимость исследования микроРНК. N. Nakamura et al. было показано повышение экспрессии данной микроРНК в сыворотке крови и перитонеальной жидкости при РЯ, по сравнению с изменениями при эндометриозе. В связи с зависимостью экспрессии miR-486-5p от степени распространения заболевания предложено использовать определение данной микроРНК в качестве маркера тяжести течения эндометриоза. Также ими была показана роль данной микроРНК в индукции пролиферации при ЭКЯ и РЯ, ассоциированном с ЭКЯ [32].

Наряду с сопоставимыми данными нами были получены результаты, не опубликованные ранее при изученных нозологиях (табл. 4) или отличающиеся от опубликованных ранее данных. Ранее J. Liu et al. было показано, что при дифференцировке образцов крови пациентов с РЯ от образцов пациентов группы контроля определение miR- 4732-5p в плазме характеризуется 85,7% чувствительностью и 82,4% специфичностью (AUC=0,889, р<0,0001) [36]. Также авторами было показано, что данная микроРНК является информативным маркером для мониторинга пациентов с ранними стадиями для выявления прогрессирования заболевания (р=0,018). В отличие от упомянутого исследования, нами не было обнаружено статистически значимого изменения miR-4732-5p у пациентов с серозным РЯ, по сравнению с пациентами с ЭКЯ и наблюдаемыми женщинами группы контроля. Следует отметить, что изменения данной микроРНК нами были отмечены при ЭКЯ — у пациентов наблюдалось повышение ее экспрессии, по сравнению с наблюдаемыми женщинами группы контроля (padj=0,015).

Для более углубленного изучения биологической роли выявленных ДЭМ нами был проведен анализ регуляторной функции потенциальных 1288 мишеней в нескольких сигнальных путях, наиболее вероятно задействованных в патогенезе и дифференцировке изучаемых пролиферативных процессов. На основании выявленных наиболее статистически значимых различий были идентифицированы гены- мишени нескольких сигнальных путей: AKT1, ATM, BARD1, BAX, BCL2, BRCA1, CASP3, CDK4, CHEK1, CHEK2, JAK1, MDM2, PLK1, PTEN, RB1, SMAD2, SMAD3, TP53 (табл. 2). Функциональная кластеризация обозначенных генов-мишеней показала их участие в процессах регуляции клеточного цикла, сигнальных путей факторов роста, апоптоза и репарации ДНК, что может служить подтверждением их специфичности для патогенетических процессов, ассоциированных с исследуемыми пролиферативными патологиями яичников. Было установлено сходство анализируемых групп по активации нескольких сигнальных путей — многочисленные взаимосвязи между выявленными генами-мишенями микроРНК указывают на их существенный вклад и совместную вовлеченность в процессы пролиферативной активности при РЯ и ЭКЯ. Подробно механизмы регуляции изложены в ряде работ, в частности в обзоре 2021 г. A. Gajek et al. [39], а также в оригинальных исследованиях [40, 41].

Заключение

Благодаря стабильности экзосом плазмы крови изучение микроРНК, входящих в их состав, является перспективным направлением поиска маркеров пролиферативных процессов. Наличие ДЭМ в экзо- сомах плазмы крови свидетельствует о специфических изменениях при РЯ и ЭКЯ. В качестве маркерных для РЯ и ЭКЯ были выбраны следующие ДЭМ: miR-141-3p, -199a-5p, -200b-3p, -203a-3p, -224-5p, -4488 (повышены при РЯ), miR-486-5p, -3184-3p, -4732-5p, -423-5p, -92b-5p (повышены при ЭКЯ).

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

Список литературы

1. Siegel R.L., Miller K.D., Jemal A. Cancer statistics, 2017. CA Cancer J. Clin. 2017; 67(1):7-30. https://dx.doi.org/10.3322/caac.21387.2017; 67(1): 7-30.

2. Serdar E. Bulun M.D. Endometriosis. N. Engl. J. Med. 2009; 360: 268-79. https://dx.doi.org/10.1056/NEJMra0804690.

3. Павлович С.В., Юрова М.В., Мелкумян А.Г., Франкевич В.Е., Хабас Г.Н., Чаговец В.В. Биомаркеры при новообразованиях яичников: возможности, ограничения и перспективы применения у женщин репродуктивного возраста. Акушерство и гинекология. 2019; 11: 65-73.

4. Fassbender A., Burney R.O., O D.F., D'Hooghe T., Giudice L. Update on biomarkers for the detection of endometriosis. Biomed. Res. Int. 2015; 2015: 130854. https://dx.doi.org/10.1155/2015/130854.

5. Aznaurova Y.B., Zhumataev M.B., Roberts T.K., Aliper A.M., Zhavoronkov A.A. Molecular aspects of development and regulation of endometriosis. Reprod. Biol. Endocrinol. 2014 Jun 13; 12:50. https://dx.doi.org/10.1186/ 1477-7827-12-50.

6. Cortez M.A., Bueso-Ramos C., Ferdin J., Lopez-Berestein G., Sood A.K., Calin G.A. MicroRNAs in body fluids-the mix of hormones and biomarkers. Nat. Rev. Clin. Oncol. 2011; 8(8): 467-77. https://dx.doi.org/10.1038/nrclinonc.2011.76.

7. Mirzaei H., Gholamin S., Shahidsales S., Sahebkar A., Jaafari M.R., Mirzaei H.R. et al. MicroRNAs as potential diagnostic and prognostic biomarkers in melanoma. Eur. J. Cancer. 2016; 53: 25-32. https://dx.doi.org/10.1016/ j.ejca.2015.10.009.

8. Iorio M.V., Visone R., Di Leva G., Donati V., Petrocca F., Casalini P. et al. MicroRNA signatures in human ovarian cancer. Cancer Res. 2007; 67(18): 8699-707. https://dx.doi.org/10.1158/0008-5472.CAN-07-1936.

9. Filigheddu N., Gregnanin I., Porporato P.E., Surico D., Perego B., Galli L. et al. Differential expression of micrornas between eutopic and ectopic endometrium in ovarian endometriosis. J. Biomed. Biotechnol. 2010; 2010: 369549. https://dx.doi.org/10.1155/2010/369549.

10. Anfossi S., Babayan A., Pantel K., Calin G.A. Clinical utility of circulating non-coding RNAs — an update. Nat. Rev. Clin. Oncol. 2018; 15(9): 541-63. https://dx.doi.org/10.1038/s41571-018-0035-x.

11. Farooqi A.A., Desai N.N., Qureshi M.Z., Librelotto D.R.N., Gasparri M.L., Bishayee A. et al. Exosome biogenesis, bioactivities and functions as new delivery systems of natural compounds. Biotechnol. Adv. 2018; 36(1): 328-34. https://dx.doi.org/10.1016/j.biotechadv.2017.12.010.

12. Giannopoulou L., Zavridou M., Kasimir-Bauer S., Lianidou E.S. Liquid biopsy in ovarian cancer: the potential of circulating miRNAs and exosomes. Transl. Res. 2019; 205: 77-91. https://dx.doi.org/10.1016/j.trsl.2018.10.003.

13. Aboutalebi H., Bahrami A., Soleimani A., Saeedi N., Rahmani F., Khazaei M. et al. The diagnostic, prognostic and therapeutic potential of circulating microRNAs in ovarian cancer. Int. J. Biochem. Cell Biol. 2020; 124: 105765. https://dx.doi.org/10.1016/j.biocel.2020.105765.

14. Love M.I., Huber W., Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15(12): 550. https://dx.doi.org/10.1186/s13059-014-0550-8.

15. Licursi V., Conte F., Fiscon G., Paci P. MIENTURNET: An interactive web tool for microRNA-target enrichment and network-based analysis. BMC Bioinformatics. 2019; 20(1):545. https://dx.doi.org/10.1186/ s12859-019-3105-x.

16. Glaab E., Baudot A., Krasnogor N., Schneider R., Valencia A. EnrichNet: network-based gene set enrichment analysis. Bioinformatics. 2012; 28(18): i451- 7. https://dx.doi.org/10.1093/bioinformatics/bts389.

17. Pegtel D.M., Gould S.J. Exosomes. Annu. Rev. Biochem. 2019; 88: 487-514. https://dx.doi.org/10.1146/annurev-biochem-013118-111902.

18. Kim K.M., Abdelmohsen K., Mustapic M., Kapogiannis D., Gorospe M. RNA in extracellular vesicles. Wiley Interdiscip. Rev. RNA. 2018; 8(4): 10.1002/ wrna.1413. https://dx.doi.org/10.1002/wrna.1413.

19. Cho S., Mutlu L., Grechukhina O., Taylor H.C. Circulating microRNAs as potential biomarkers for endometriosis. Fertil. Steril. 2015; 103(5): 1252-60. e1. 103(5):1252- 103(5):1252-60. e1. https://dx.doi.org/10.1016/ j.fertnstert.2015.02.013.

20. Zafari N., Tarafdari A.M., Izadi P. A panel of plasma miRNAs 199b-3p, 224-5p and Let-7d-3p as non-invasive diagnostic biomarkers for endometriosis. Reprod. Sci. 2021; 28(4): 991-9. https://dx.doi.org/10.1007/ s43032-020-00415-z.

21. Vanhie A., O D., Peterse D., Beckers A., Cuellar A., Fassbender A. et al. Plasma miRNAs as biomarkers for endometriosis. Hum. Reprod. 2019; 34(9): 1650-60. https://dx.doi.org/10.1093/humrep/dez116.

22. Moustafa S., Burn M., Mamillapalli R., Nematian S., Flores V., Taylor H.S. Accurate diagnosis of endometriosis using serum microRNAs. Am. J. Obstet. Gynecol. 2020; 223(4): 557.e1-557.e11. https://dx.doi.org/10.1016/ j. ajog. 2020.02.050.

23. Masoumi-Dehghi S., Babashah S. microRNA-141-3p-containing small extracellular vesicles derived from epithelial ovarian cancer cells promote endothelial cell angiogenesis through activating the JAK / STAT3 and NF-к B signaling pathways. J. Cell Commun. Signal. 2020; 14(2): 233-44. https://dx.doi.org/10.1007/s12079-020-00548-5.

24. Leonova A., Turpin V.E., Agarwal S.K., Leonardi M., Foster W.G. A critical appraisal of the circulating levels of differentially expressed microRNA in endometriosis. Biol. Reprod. 2021; 105(5): 1075-85. https://dx.doi.org/10.1093/ biolre/ioab134.

25. Truong G., Guanzon D., Kinhal V., Elfeky O., Lai A., Longo S. et al. Oxygen tension regulates the miRNA profile and bioactivity of exosomes released from extravillous trophoblast cells — Liquid biopsies for monitoring complications of pregnancy. PLoS One. 2017; 12(3): e0174514. https://dx.doi.org/10.1371/ journal.pone.0174514.

26. Hu J., Tang T., Zeng Z., Wu J., Tan X., Yan J. The expression of small RNAs in exosomes of follicular fluid altered in human polycystic ovarian syndrome. PeerJ. 2020; 8: e8640. https://dx.doi.org/10.7717/peerj.8640.

27. Papari E., Noruzinia M., Kashani L., Foster W.G. Identification of candidate microRNA markers of endometriosis with the use of nextgeneration sequencing and quantitative real-time polymerase chain reaction. Fertil. Steril. 2020; 113(6): 1232-41. https://dx.doi.org/10.1016/ j.fertnstert.2020.01.026.

28. Bjorkman S., Taylor H.S. MicroRNAs in endometriosis: biological function and emerging biomarker candidates. Biol. Reprod. 2019; 100(5): 1135-46. https://dx.doi.org/10.1093/biolre/ioz014.

29. Roman-Canal B., Moiola C.P, Gatius S., Bonnin S., Ruiz-Miro M., Gonzalez E. et al. EV-associated miRNAs from peritoneal lavage are a source of biomarkers in endometrial cancer. Cancers (Basel). 2019; 11(6): 839. https://dx.doi.org/10.3390/cancers11060839.

30. Taylor D.D., Gercel-Taylor C. MicroRNA signatures of tumor-derived exosomes as diagnostic biomarkers of ovarian cancer. Gynecol. Oncol. 2008; 110(1): 13-21. https://dx.doi.org/10.1016/j.ygyno.2008.04.033.

31. Chen Y., Wang K., Xu Y., Guo P., Hong B., Cao Y. et al. Alteration of myeloid-derived suppressor cells, chronic inflammatory cytokines, and exosomal miRNA contribute to the peritoneal immune disorder of patients with endometriosis. Reprod. Sci. 2018; 26(8): 1130-8. https://dx.doi.org/10.1177/1933719118808923.

32. Nakamura N., Terai Y., Nunode M., Kokunai K., Konishi H., Taga S. et al. The differential expression of miRNAs between ovarian endometrioma and endometriosis-associated ovarian cancer. J. Ovarian Res. 2020; 13(1): 51. https://dx.doi.org/10.1186/s13048-020-00652-5.

33. Bj0rnetr0 T., Redalen K.R., Meltzer S., Thusy anthan N.S., Samiappan R., Jegerschold C. et al. An experimental strategy unveiling exosomal microRNAs 486-5p, 181a-5p and 30d-5p from hypoxic tumour cells as circulating indicators of high-risk rectal cancer. J. Extracell. Vesicles. 2019; 8(1): 1567219. https://dx.doi.org/10.1080/20013078.2019.1567219.

34. Liu C., Li M., Hu Y., Shi N., Yu H., Liu H., Lian H. miR-486-5p attenuates tumor growth and lymphangiogenesis by targeting neuropilin-2 in colorectal carcinoma. Onco Targets Ther. 2016; 9: 2865-71. https://dx.doi.org/10.2147/ OTT.S103460.

35. Narasimhan A., Ghosh S., Stretch C., Greiner R., Bathe O.F., Baracos V., Damaraju S. Small RNAome pro fi ling from human skeletal muscle: novel miRNAs and their targets associated with cancer cachexia. J. Cachexia Sarcopenia Muscle. 2017; 8(3): 405-16. https://dx.doi.org/10.1002/ jcsm.12168.

36. Liu J., Yoo J., Ho J.Y., Jung Y., Lee S., Hur S.Y., Choi Y.J. Plasma-derived exosomal miR-4732-5p is a promising noninvasive diagnostic biomarker for epithelial ovarian cancer. J. Ovarian Res. 2021; 14(1): 59. https://dx.doi.org/10.1186/s13048-021-00814-z.

37. Cochrane D.R., Howe E.N., Spoelstra N.S., Richer J.K. Loss of miR-200c: a marker of aggressiveness and chemoresistance in female reproductive cancers. J. Oncol. 2010; 2010: 821717. https://dx.doi.org/10.1155/2010/821717.

38. Cao Q., Lu K., Dai S., Hu Y., Fan W. Clinicopathological and prognostic implications of the miR-200 family in patients with epithelial ovarian cancer. Int. J. Clin. Exp. Pathol. 2014; 7(5): 2392-401.

39. Gajek A., Gralewska P., Marczak A., Rogalska A. Current implications of microRNAs in genome stability and stress responses of ovarian cancer. Cancers (Basel). 2021; 13(11): 2690. https://dx.doi.org/10.3390/cancers13112690.

40. Бобров М.Ю., Балашов И.С., Филиппова Е.С., Альмова И.К., Тимофеева А.В., Гусар В.А., Боровиков П.И., Хилькевич Е.Г., Чупрынин В.Д., Павлович С.В. Оценка экспрессии микроРНК в очагах ретроцервикаль- ного эндометриоза. Акушерство и гинекология. 2018; 6: 55-61.

41. Бобров М.Ю., Балашов И.С., Филиппова Е.С., Альмова И.К., Хилькевич Е.Г., Павлович С.В., Наумов В.А., Боровиков П.И., Сухих Г.Т. Использование транскриптомных баз данных для анализа патогенетических факторов эндометриоза. Акушерство и гинекология. 2017; 4: 34-44.

Поступила 11.01.2022

Принята в печать 25.02.2022

Об авторах / Для корреспонденции

Юрова Мария Владимировна, аспирант кафедры акушерства, гинекологии, перинатологии и репродуктивного здоровья ИПО, Первый МГМУ им. И.М. Сеченова Минздрава России (Сеченовский Университет), 119991, Россия, Москва, ул. Трубецкая, д. 8, стр. 2; специалист, Национальный медицинский исследовательский центр акушерства, гинекологии и перинатологии имени академика В.И. Кулакова Минздрава России, 117997, Россия, Москва, ул. Академика Опарина, д. 4, +7(495)609-14-00, m_yurova@oparina4.ru, https://orcid.org/0000-0002-0179-7635
Эльдаров Чупалав Максудович, к.б.н., с.н.с. лаборатории молекулярной патофизиологии, Национальный медицинский исследовательский центр акушерства, гинекологии и перинатологии им. академика В.И. Кулакова Минздрава России, 117997, Россия, Москва, ул. Академика Опарина, д. 4, +7(495)438-77-00, ch_eldarov@oparina4.ru
Бобров Михаил Юрьевич, к.х.н., руководитель лаборатории молекулярной патофизиологии, Национальный медицинский исследовательский центр акушерства, гинекологии и перинатологии им. академика В.И. Кулакова Минздрава России, 117997, Россия, Москва, ул. Академика Опарина, д. 4, mbobr@mail.ru
Хабас Григорий Николаевич, к.м.н., руководитель отделения инновационной онкологии и гинекологии, Национальный медицинский исследовательский центр акушерства, гинекологии и перинатологии им. академика В.И. Кулакова» Минздрава России, 117997, Россия, Москва, ул. Академика Опарина, д. 4., +7(495)438-78-33, g_khabas@oparina4.ru, https://orcid.org/0000-0002-5011-9152
Павлович Станислав Владиславович, к.м.н., Ученый секретарь, Национальный медицинский исследовательский центр акушерства, гинекологии и перинатологии им. академика В.И. Кулакова Минздрава России, 117997, Россия, Москва, ул. Академика Опарина, д. 4, +7(495)438-20-88, s_pavlovich@oparina4.ru, https://orcid.org/0000-0002-1313-7079
Автор, ответственный за переписку: Мария Владимировна Юрова, +7(926)385-45-48, hi5melisa@gmail.com, m_yurova@oparina4.com

Вклад авторов: Юрова М.В. — разработка дизайна представленного исследования, анализ данных доступной литературы, создание и ведение базы данных, статистическая обработка и интерпретация клинических данных, написание текста, редактирование, подготовка документации; Бобров М.Ю. — разработка дизайна представленного исследования, проведение всех этапов практической части исследования, редактирование, биологическая интерпретация полученных данных, окончательное утверждение публикуемой версии рукописи, подготовка документации; Эльдаров Ч.М. — разработка дизайна представленного исследования, анализ данных доступной литературы, работа с базами данных, статистическая обработка и интерпретация данных, написание текста, редактирование, подготовка документации; Хабас Г.Н. — оперативное лечение пациентов, клиническая интерпретация полученных данных, критический пересмотр написанного текста, внесение правок и дополнений, окончательное утверждение публикуемой версии рукописи; Павлович С.В. — разработка дизайна представленного исследования, редактирование, клиническая интерпретация полученных данных, окончательное утверждение публикуемой версии рукописи, подготовка документации.
Конфликт интересов: Авторы заявляют об отсутствии конфликта интересов.
Финансирование: Исследование выполнено в рамках государственного задания Минздрава РФ «Разработка подходов для неинвазивной диагностики эндометриоза на основе омиксных технологий»; № госрегистрации: АААА-А19-119021490132-9 на период 01.01.19—31.12.21 гг.
Одобрение Этического комитета: Данное исследование одобрено комиссиями по этике биомедицинских исследований — локальным Этическим комитетом ФГАОУ ВО «Первый МГМУ им. И.М. Сеченова» Минздрава России (Сеченовский Университет) (протокол № 16-19 от 05.12.2019) и комиссией по вопросам медицинской этики ФГБУ «НМИЦ АГП им. В.И. Кулакова» Минздрава России (протокол № 10 от 05 декабря 2019 г.).
Согласие пациентов на публикацию: Пациенты подписали информированное согласие на публикацию своих данных.
Обмен исследовательскими данными: Данные, подтверждающие выводы этого исследования, доступны по запросу у автора, ответственного за переписку, после одобрения ведущим исследователем.
Для цитирования: Юрова М.В., Эльдаров Ч.М., Бобров М.Ю., Хабас Г.Н., Павлович С.В. Экспрессия экзосомальных микроРНК при раке яичников высокой степени злокачественности и эндометриоидных кистах яичников.
Акушерство и гинекология. 2022; 3: 68-79
https://dx.doi.org/10.18565/aig.2022.3.68-79

Также по теме

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