ISSN 0300-9092 (Print)
ISSN 2412-5679 (Online)

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

Тимофеева А.В., Федоров И.С., Гохберг Я.А., Калинина Е.А.

ФГБУ «Национальный медицинский исследовательский центр акушерства, гинекологии и перинатологии имени академика В.И. Кулакова» Минздрава России, Москва, Россия
Современные методы диагностики рецептивности эндометрия, введенные в клиническую практику (ERA, Igenomix; ER-Map, iGLS; ERPeak, CooperGenomics), основаны на транскриптомном анализе ткани эндометрия при инвазивном способе его получения, что отменяет перенос эмбриона в том же цикле, что и выполнение биопсии эндометрия. Цель: Оценка рецептивности эндометрия путем количественного определения малых некодирующих (мнк)РНК в маточной жидкости (МЖ) в день переноса криоконсервированного эмбриона (КЭ) на фоне циклической гормональной терапии (ЦГТ) и построение модели логистической регрессии расчета вероятности готовности эндометрия к имплантации эмбриона путем сопоставления транскриптома МЖ от пациенток с положительным и отрицательным исходом программы ВРТ. Материалы и методы: В исследование были включены 54 женщины, которым проводили аспирацию МЖ в объеме 5–50 мкл в зависимости от уровня ее секреции непосредственно перед переносом КЭ с использованием катетера (COOK, Австралия). мнкРНК, выделенные из МЖ набором miRNeasy Serum/Plasma Kit (Qiagen), были проанализированы методом глубокого секвенирования на платформе NextSeq 500/550 (Illumina, USA) c последующей валидацией полученных данных методом количественной ПЦР в реальном времени наборами miScript II RT Kit и miScript SYBR Green PCR Kit (Qiagen, Hilden, Germany). Результаты: Образцы МЖ были классифицированы на 2 группы – «рецептивный эндометрий» и «нерецептивный эндометрий» – в зависимости от результатов проведенной программы ВРТ (наличие и отсутствие имплантации соответственно). По профилю мнкРНК в образцах МЖ и толщине эндометрия на момент переноса эмбриона в полость матки были построены семь моделей логистической регрессии, среди которых наиболее точной является комбинация содержания miR-1180-3p в МЖ и толщины эндометрия (71% чувствительность, 88% специфичность) ввиду отсутствия зависимости включенных в нее переменных по данным корреляционного анализа по Спирмену (r=0,02; p=0,9) и статистической значимости всех входящих в модель коэффициентов (р<0,05). Заключение: Индивидуализация подхода в определении «окна имплантации» при проведении ЭКО на фоне ЦГТ в соответствии с уровнем мнкРНК в МЖ может играть важную роль в повышении эффективности программ ВРТ. Ограничением использования в клинической практике модели логистической регрессии, построенной в настоящем исследовании на обучающей выборке пациенток, является отсутствие информации об имплантационном потенциале переносимой в полость матки бластоцисты. Отрицательный исход программы ВРТ может быть обусловлен не отсутствием рецептивного эндометрия, а сниженным качеством самого эмбриона. Одновременное определение рецептивности эндометрия и имплантационного потенциала эмбриона по профилю мнкРНК приведет к снижению процента ложноотрицательных результатов и улучшению качества модели. Необходимо расширить обучающую выборку и проверить точность построенных моделей на независимой тестовой выборке. Вклад авторов: Тимофеева А.В. – получение экспериментальных данных, написание статьи; Федоров И.С. – получение экспериментальных данных, статистическая обработка данных; Гохберг Я.А. – сбор образцов маточной жидкости, создание клинической базы пациенток; Калинина Е.А. – редактирование статьи. Конфликт интересов: Авторы заявляют об отсутствии конфликта интересов. Финансирование: Работа выполнена при финансовой поддержке Российского научного фонда в рамках гранта № 22-15-00363 «Эпигенетические и биохимические аспекты патологии беременности при нарушениях инвазивных свойств трофобласта: от ранней диагностики к профилактике материнской и перинатальной заболеваемости» в соответствии с соглашением № 22-15-00363 между Российским научным фондом, руководителем проекта Тимофеевой А.В. и НМИЦ АГП им. В. И. Кулакова о предоставлении гранта на проведение фундаментальных научных исследований и поисковых научных исследований от 13.05.2022 г. Одобрение Этического комитета: Исследование было одобрено локальным Этическим комитетом ФГБУ «НМИЦ АГП им. В.И. Кулакова» Минздрава России. Согласие пациентов на публикацию: Пациенты подписали информированное согласие на публикацию своих данных. Обмен исследовательскими данными: Данные, подтверждающие выводы этого исследования, доступны по запросу у автора, ответственного за переписку, после одобрения ведущим исследователем. Для цитирования: Тимофеева А.В., Федоров И.С., Гохберг Я.А., Калинина Е.А. Оценка рецептивности эндометрия по уровню малых некодирующих РНК в маточном аспирате у женщин на фоне циклической гормональной терапии. Акушерство и гинекология. 2023; 7: 90-102 https://dx.doi.org/10.18565/aig.2023.110

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

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

Бесплодие представляет собой глобальную проблему, затрагивающую 10–15% супружеских пар репродуктивного возраста в большинстве стран мира [1, 2]. За последние десятилетия в рамках проведения программ экстракорпорального оплодотворения (ЭКО) в лечении бесплодия достигнут огромный прогресс по оптимизации условий культивирования эмбрионов [3, 4], оценке морфологических и молекулярно-биологических критериев их качества [5–8], выбору стадии развития эмбриона, подходящей для переноса в полость матки [9, 10]. Но не только жизнеспособность и качество эмбриона играют первостепенную роль в успешной имплантации. Рецептивность эндометрия является решающим фактором в обеспечении имплантации эмбриона в течение определенного периода менструального цикла, называемого «окном имплантации» [11, 12]. Эта средняя стадия секреторной фазы наступает через 6–10 дней после пика лютеинизирующего гормона и регулируется цитокинами, факторами роста и другими молекулами, ключевыми для развития эндометрия и эмбриона [13, 14]. Выявлено, что около 60% повторных неудач вспомогательных репродуктивных технологий (ВРТ) обусловлено смещением «окна имплантации» [12, 15, 16]; поэтому определение готовности эндометрия к имплантации эмбриона является важной задачей в области репродукции.

Для определения рецептивности эндометрия были разработаны различные методы, включая инструментальные и инвазивные методы исследования [17]. Согласно данным метаанализа, специфичность и/или чувствительность инструментальных методов исследования по анализу толщины, объема, трехслойности эндометрия, кровотока в эндометрии, гистероскопии не позволяют достоверно оценить рецептивность эндометрия с целью прогнозирования наступления беременности в программах ВРТ [12]. Ультраструктурными маркерами рецептивности эндометрия являются пиноподии, считающиеся предпочтительными местами взаимодействия эмбриона и эндометрия [18, 19]. Кроме того, появление пиноподий совпадает с увеличением экспрессии других биомаркеров «окна имплантации», таких как муцин 1, ингибирующий лейкемию фактор (LIF) и его рецептор (интегрин αVβ3) [20], гомеодомен-содержащий транскрипционный фактор 10 (HOX10) [21]. Оценка числа и формы пиноподий при гистологическом анализе биоптата эндометрия была использована для разработки метода персонализированного переноса эмбрионов [22]. Но во избежание ошибок в подсчете числа пиноподий в образце люминальной поверхности эндометрия необходим микроскопический просмотр не менее 60 областей образца ввиду высокой вариабельности числа пиноподий (от 5 до 20%) [23], что является трудоемким процессом и исключает его использование в рутинной клинической практике.

Использование метода транскриптомного анализа биоптата эндометрия (мРНК и микроРНК) позволило выявлять пациентов со смещением окна имплантации и повысить частоту наступления беременности, особенно в группе с повторными неудачами имплантации в программе ЭКО [24–26]. Несмотря на большую информативность и точность метода транскриптомного анализа в сравнении с гистологическим исследованием биоптата эндометрия, его дороговизна и инвазивность являются препятствиями для рутинного применения в клинической практике.

Альтернативным методом оценки состояния эндометрия является исследование маточной жидкости (МЖ), аспирация которой является минимально инвазивной и безопасной и может быть выполнена непосредственно перед переносом эмбрионов без отрицательного влияния на последующую имплантацию, в отличие от общепринятой биопсии эндометрия [27, 28]. Было выявлено изменение паттерна экспрессии малых некодирующих (мнк)РНК экзосом МЖ в рецептивную стадию в сравнении с прорецептивной стадией секреторной фазы естественного менструального цикла и при контролируемой овариальной стимуляции [29], но связи с наступлением беременности выявить не удалось.

Для минимизации влияния на эндометрий контролируемой овариальной стимуляции [30] нами была выбрана модель переноса эмбриона в криоцикле.

Цель исследования: оценка рецептивности эндометрия путем количественного определения мнкРНК в МЖ в день переноса криоконсервированного эмбриона (КЭ) на фоне циклической гормональной терапии (ЦГТ) и построение модели логистической регрессии расчета вероятности готовности эндометрия к имплантации эмбриона путем сопоставления транскриптома МЖ от пациенток с положительным и отрицательным исходом программы ВРТ.

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

Пациенты

В отделение вспомогательных репродуктивных технологий в лечении бесплодия им. проф. Б.В. Леонова ФГБУ «НМИЦ АГП им. В.И. Кулакова» обратились 54 супружеские пары для проведения программы переноса КЭ на фоне ЦГТ. Критериями включения в исследование были возраст пациентов от 29 до 39 лет, нормальный овариальный резерв у женщин с трубно-перитонеальным фактором бесплодия, нормозооспермия или мужской фактор бесплодия без выраженной патозо­оспермии. Критериями невключения в исследование были: 1) наличие верифицированного наружного генитального эндометриоза III–IV стадии; 2) онкологические заболевания; 3) синдром поликистозных яичников; 4) патология эндометрия; 5) интерстициальная и/или субсерозная миома матки более 4 см; 6) субмукозная миома, деформирующая полость матки; 7) генетические аномалии; 8) пороки развития половых органов; 9) перенесенные оперативные вмешательства на яичниках; 10) тяжелые формы мужского бесплодия.

Всем пациентам было проведено полное клинико-лабораторное обследование согласно Клиническим рекомендациям «Женское бесплодие», а также в соответствии с Приказом Минздрава России от 31.07.2020 N 803н «О порядке использования вспомогательных репродуктивных технологий, противопоказаниях и ограничениях к их применению». У всех пациентов, включенных в исследование, было получено информированное добровольное согласие на осуществление аспирации МЖ в день переноса эмбриона с помощью одноразовых гибких катетеров.

Протокол подготовки эндометрия к переносу размороженного эмбриона на фоне ЦГТ

На 2–4-й день менструального цикла пациентам был назначен препарат эстрадиола валерат в дозировке 2 мг/сут, и после достижения нормальной толщины эндометрия (>7 мм) на 14–15-й день цикла вводили гестагены (200 мг/сут микронизированного прогестерона или 10 мг/сут дидрогестерона). День переноса КЭ в полость матки рассчитывали путем прибавления числа дней культивирования эмбриона с момента оплодотворения яйцеклетки в предыдущем цикле стимуляции суперовуляции к дате начала введения гестагенов в текущем криоцикле, как показано на рисунке 1.

93-1.jpg (180 KB)

Наступление беременности диагностировали при получении положительного результата β-хорионического гонадотропина человека через 14 дней после переноса эмбриона и при визуализации плодного яйца в полости матки на 21-й день.

Аспирация МЖ

Аспирация МЖ проводилась непосредственно перед проведением переноса КЭ с использованием гибкого катетера (COOK, Австралия), который соединялся со стерильным шприцом (1 мл). Полученный материал МЖ объемом 5–50 мкл помещали в стерильный эппендорф и доводили объем физиологическим раствором 0,9% NaCl до 200 мкл. Далее материал транспортировали в лабораторию прикладной транскриптомики отдела системной биологии в репродукции ФГБУ «НМИЦ АГП им. В.И. Кулакова» для хранения при -80°С и дальнейшего молекулярно-биологического анализа МЖ.

Выделение мнкРНК из образцов МЖ

Из собранных образцов МЖ были выделены РНК колоночным способом с использованием набора miRNeasy Serum/Plasma Kit (Qiagen).

Глубокое секвенирование мнкРНК

Из 14 мкл элюата колонки miRNeasy Serum/Plasma Kit (Qiagen), содержащего мнкРНК из МЖ, 7 мкл были использованы для синтеза кДНК-библиотек набором NEBNext Multiplex Small RNA Library Prep Set for Illumina (Set11 and Set2, New England Biolab®), амплифицированных в течение 20 циклов в ходе полимеразной цепной реакции (ПЦР), очищенных в 6% полиакриламидном геле и секвенированных на платформе NextSeq 500/550 (Illumina). Секвенированные последовательности в диапазоне от 16 до 50 п.н. были картированы на базы данных человека GRCh38.p15, miRBase v21 и piRNABase с использованием алгоритма Bowtie [31]. Анализ дифференциальной экспрессии мнкРНК проводили с помощью программного пакета DESeq2 [32].

Обратная транскрипция и количественная ПЦР в реальном времени

Из 14 мкл элюата колонки miRNeasy Serum/Plasma Kit (Qiagen), содержащего мнкРНК из МЖ, 5 мкл были использованы для синтеза кДНК набором miScript II RT Kit (Qiagen) по протоколу фирмы-производителя. Количественная ПЦР в реальном времени была проведена с использованием набора miScript SYBR Green PCR Kit (Qiagen, Hilden, Germany) и смысловых праймеров, специфичных для определенных мкРНК и пивиРНК (табл. 1). Программа ПЦР: (1) 15 мин при 95°С и (2) 40 циклов: 94°С в течение 15 с, оптимизированная температура отжига (45–61,6°С) в течение 30 с и 70°С в течение 30 с; (3) нагревание реакционной смеси от 65 до 95°C с шагом 0,1°C для построения кривой плавления продукта ПЦР в термоциклере StepOnePlus (Applied Biosystems). Относительную экспрессию мнкРНК в МЖ определяли методом ∆Ct, используя hsa_piR_017716 в качестве референсной РНК.

94-1.jpg (162 KB)

Статистический анализ

Для статистической обработки результатов использовали скрипты, написанные на языке R [33], и программу RStudio [34]. Соответствие анализируемых параметров закону нормального распределения оценивали по значениям теста Шапиро–Уилка. Статистический анализ проводили с помощью теста Манна–Уитни при парном сравнении в случае, когда распределение не соответствовало закону нормального распределения. При распределении признаков с распределением, отличающимся от нормального, их описывали в виде медианы (Me) и квартилей Q1 и Q3 в формате Me (Q1; Q3). Величину порогового уровня значимости p принимали равной 0,05. Если значение p было меньше 0,001, то p указывали в формате p<0,001.

Модели логистической регрессии разрабатывали с использованием программы RStudio путем поэтапного включения и исключения мнкРНК-предикторов рецептивного эндометрия в соответствии с их вкладом в модель. Прогностическую способность модели оценивали методом ROC-анализа (Receiver operating characteristic) по величине AUC (Area Under Curve), статистической значимости, уровню специфичности и чувствительности.

Результаты

На первом этапе исследования был проведен ретроспективный анализ анамнестических и лабораторных данных супружеских пар до вступления в программу ЭКО, эмбрионального этапа в ходе стимуляции суперовуляции и последующего криопротокола на фоне ЦГТ в двух сформированных группах в зависимости от наличия или отсутствия имплантации эмбриона (табл. 2).

96-1.jpg (229 KB)

Из таблицы 2 следует, что в группе с отрицательным результатом программы ВРТ при переносе КЭ наблюдали статистически значимое повышение уровня тестостерона в 1,5 раза (р=0,02) и понижение уровня прогестерона в 2,1 раза (р<0,001) относительно группы с наступившей беременностью до вступления в программу ЭКО. При этом у женщин с положительным результатом программы ВРТ выявлены более высокие значения толщины эндометрия в день переноса эмбриона, по сравнению с группой женщин, у которых беременность не наступила (р<0,001). По анамнестическим данным и характеристикам эмбрионального этапа в ходе программы стимуляции суперовуляции сопоставляемые группы не различались.

На втором этапе исследования методом глубокого секвенирования была проведена количественная оценка мнкРНК (мкРНК и пивиРНК) в 31 из 54 образцов МЖ, собранных непосредственно перед проведением переноса КЭ. В 15 образцах МЖ женщин с отсутствием имплантации эмбриона и в 16 образцах МЖ женщин с наступившей беременностью было идентифицировано 289 мкРНК и 488 пивиРНК. Из данного списка мнкРНК в программе RStudio путем поэтапного включения и исключения каждой молекулы были найдены оптимальные комбинации мнкРНК-маркеров рецептивного эндометрия в соответствии с их вкладом в построение моделей логистической регрессии (рис. 2), где в качестве зависимой переменной (переменной отклика) выступала степень готовности эндометрия к имплантации (0 – имплантация, 1 – отсутствие имплантации). Все модели, построенные по содержанию мкРНК (рис. 2А) и пивиРНК (рис. 2Б) в МЖ, были статистически значимы и обладали высокой специфичностью (84–100%), а значит, высокой диагностической ценностью наличия рецептивного эндометрия в день переноса КЭ. Более низкие значения чувствительности построенных моделей (65–87%) могут быть обусловлены отсутствием информации об имплантационном потенциале эмбриона, переносимого в полость матки, которую можно получить до этапа криоконсервации на стадии морулы или бластоцисты [7], так как причиной отсутствия имплантации может быть не только нарушение формирования рецептивного эндометрия, но и качество самого эмбриона.

95-1.jpg (57 KB)

Данные секвенирования были валидированы методом количественной ПЦР в реальном времени на всей выборке пациентов, указанных в таблице 2. По значениям «-ΔСt» для 18 молекул мнкРНК, перечисленных в таблице 1, и величине толщины эндометрия были построены модели логистической регрессии в программе RStudio путем нахождения оптимальной комбинации объясняющих переменных при поэтапном их включении и исключении в соответствии с вкладом в построение модели, где в качестве зависимой переменной (переменной отклика) выступала готовность эндометрия к имплантации эмбриона. Была найдена комбинация из объясняющих переменных (модель 1 рис. 3), имеющая наибольшую площадь под кривой (AUC=0,89) и состоящая из семи мнкРНК и толщины эндометрия. Методом обратного исключения из модели 1 рисунка 3 пошагово удаляли по одной самой незначимой объясняющей переменной с получением моделей 2, 3, 4, 5, 6, 7, где каждая последующая модель является производной от предыдущей (табл. 3). В результирующей модели 7 рисунка 3 все объясняющие переменные оказались статистически значимыми (табл. 3) и независимыми по данным корреляционного анализа по Спирмену (r=0,02; p=0,9). Cпецифичность и чувствительность модели 7 рисунка 3 составила 88 и 71% соответственно, что говорит о высокой диагностической значимости сочетанного измерения уровня miR-1180-3p в МЖ и толщины эндометрия в определении готовности эндометрия к имплантации КЭ в день его переноса. Формула расчета вероятности рецептивного эндометрия для модели 7 рисунка 3 представлена ниже:

96-2.jpg (4 KB)

где: х1 – ΔCt(miR-1180-3p), х2­ – толщина эндо­метрия (мм).

97-1.jpg (229 KB)

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

98-1.jpg (156 KB)

Важно отметить, что однофакторная модель, включающая лишь толщину эндометрия, обладает меньшей специфичностью для оценки его репептивности (модель 8 рис. 3, Sp=76%), как указано в таблице 3.

При анализе функциональной значимости мнкРНК, участвующих в моделях рисунка 3, были идентифицированы потенциальные гены-мишени пивиРНК с использованием алгоритма miRanda [8], а гены-мишени мкРНК были найдены в программе Funrich (http://www.funrich.org/). При анализе обогащения сигнальных путей, в которых задействованы выявленные гены-мишени мкРНК и пивиРНК, обнаружено их участие в процессах клеточной адгезии и проведения сигнала внутрь клеток через рецепторы ростовых факторов и интерлейкинов, необходимых для имплантации эмбриона и его последующего роста и развития (рис. 4).

99-1.jpg (130 KB)

Обсуждение

Существенным недостатком современных методов диагностики рецептивности эндометрия, введенных в клиническую практику (ERA, Igenomix; ER-Map, iGLS; ERPeak, CooperGenomics), является инвазивный способ получения ткани эндометрия, что отменяет перенос эмбриона в том же цикле, что и выполнение биопсии ткани. Известно, что «окно имплантации» открывается в момент резкого изменения транскрипционного профиля в безреснитчатых, но не в реснитчатых эпителиоцитах, и в стромальных фибробластах, трансформирующихся в децидуальные клетки [35]. Поэтому высока вероятность различий соотношения реснитчатых, безреснитчатых эпителиоцитов и стромальных клеток в разных образцах биоптата эндометрия, взятых у одной и той же женщины. Более того, при различных патологических состояниях эндометрия, в том числе хроническом эндометрите, показатель отношения плотности реснитчатых и безреснитчатых эпителиоцитов значимо ниже контрольных значений [36]. Вышеуказанные обстоятельства могут привести к неправильной интерпретации получаемых данных транскриптомного анализа биоптата эндометрия. Кроме того, «окно имплантации» в менструальном цикле во время проведения инвазивного теста может быть сдвинуто относительно «окна имплантации» в последующих циклах переноса эмбриона из-за физиологической межцикловой вариабельности. Поэтому требуется использовать неинвазивный способ анализа рецептивности эндометрия, позволяющий предоставить комплексную оценку трансформации его клеток для обеспечения имплантации эмбриона.

В настоящем исследовании при сочетании методов глубокого секвенирования и количественной ПЦР в реальном времени были идентифицированы внеклеточные мнкРНК в МЖ в день переноса эмбриона, являющиеся маркерами рецептивного эндометрия. Построены модели логистической регрессии в виде комбинации количественных характеристик мнкРНК и толщины эндометрия. Наиболее точной для диагностики рецептивного эндометрия у женщин в день переноса КЭ на фоне ЦГТ оказалась модель сочетанного измерения уровня miR-1180-3p в МЖ и толщины эндометрия ввиду отсутствия зависимости включенных в нее переменных по данным корреляционного анализа по Спирмену (r=0,02; p=0,9) и статистической значимости всех входящих в модель коэффициентов (р<0,05). Ограничением настоящего исследования является малый размер выборки, исходя из формулы N=50+8m, где N – требуемый объем выборки, m – количество независимых переменных [37]. При использовании двух объясняющих переменных (уровень miR-1180-3p в МЖ и толщина эндометрия в день переноса эмбриона, модель 7 рисунка 3) объем выборки приближен к допустимому и отклоняется от референсного показателя меньше чем на 20%. Но для того, чтобы проверить диагностическую значимость исходной модели 1 рисунка 3 с использованием восьми объясняющих переменных, необходимо расширить обучающую выборку до 114 образцов, что мы и планируем сделать в будущем.

При анализе генов-мишеней мнкРНК-маркеров рецептивного эндометрия выявлено, что в основном они участвуют в процессах: 1) адгезии, опосредованные синдеканом-1, нектином, глипиканом, семейством интегринов; 2) протеолиза внеклеточного матрикса посредством активатора плазминогена урокиназного типа; 3) дифференцировки гранулоцитов и моноцитов-макрофагов под действием интерлейкинов-3 и -5, гранулоцитарно-макрофагального колониестимулирующего фактора; 4) иммунных реакциях при участии интерферона гамма; 5) клеточного роста и пролиферации, опосредованных действием тромбоцитарного фактора роста, эпидермального фактора роста, сосудистого ростового фактора, инсулиноподобного фактора роста 1, эстрогена, сфингозин-1-фосфат сигнальной системы; 6) везикулярного транспорта, контролируемого фактором АДФ-рибозилирования. Полученные данные согласуются с результатами секвенирования мРНК эпителиальных и стромальных клеток эндометрия человека в рецептивную стадию относительно прорецептивной стадии секреторной фазы менструального цикла женщин [38]. Авторы работы подчеркивают ключевую роль белков клеточной адгезии (в том числе, галектинов, интегрина бета-1, базигина и остеопонтина), регуляторов секреторных процессов и организации внеклеточного матрикса на начальных этапах имплантации эмбриона, таких как аппозиция, адгезия и инвазия. Построены сети белок-белковых взаимодействий бластоцисты и эндометрия [38–40], где одним из основных белков является интегрин бета-1. В составе мембраны клеток эндометрия интегрин бета-1 взаимодействует с гликопротеином BSG на поверхности трофэктодермы, в которой, в свою очередь, уровень интегрина бета-1 регулируется галектином-3, секретирующимся стромальными и эпителиальными клетками эндометрия, что в конечном итоге и обеспечивает миграцию и инвазию клеток трофобласта [41]. Ввиду участия целого каскада взаимодействующих друг с другом белков в процессе имплантации эмбриона, необходимо до его переноса в полость матки не только оценивать рецептивность эндометрия, но и заранее знать имплантационный потенциал самого эмбриона, ориентируясь не только на его морфологические параметры, но и на молекулярно-биологические характеристики. Важно отметить, что модель логистической регрессии для определения готовности эндометрия к имплантации эмбриона по уровню мнкРНК в день переноса КЭ на фоне ЦГТ была разработана в настоящем исследовании без учета имплантационного потенциала эмбриона, что сказалось на числе ложноотрицательных результатов (29%). В предыдущих исследованиях нами были разработаны модели логистической регрессии для определения имплантационного потенциала эмбриона по уровню мнкРНК в среде его культивирования на стадии морулы и бластоцисты в цикле контролируемой овариальной стимуляции [7] с целью дальнейшего переноса в полость матки или криоконсервации качественного эмбриона. Поэтому в дальнейшем для улучшения диагностической способности построенных моделей планируется одновременная оценка рецептивности эндометрия и имплантационного потенциала эмбриона по профилю мнкРНК.

Значимость разработанных нами моделей подтверждается литературными данными о потенциальной роли входящих в них мкРНК (miR-30d-5p, miR-425-5p, miR-34c-5p, miR-1180-3p) в процессах, ассоциированных с имплантацией эмбриона. Роль семейства miR-30 в формировании рецептивного эндометрия выявлена в различных исследованиях [42–44]; в частности, наблюдали повышение экспрессии miR-30d-5p в период окна имплантации [43]. Напротив, при повторных неудачах имплантации обнаружено резкое снижение уровня экспрессии miR-30d-5p и повышение уровня экспрессии гена-мишени супрессора передачи сигнала цитокинов 1 (SOCS1) в эндометрии на 7-е сутки после пика лютеинизирующего гормона [45]. На мышиных моделях, нокаутных по miR-30d-5p, продемонстрировано влияние дефицита miR-30d-5p в эндометрии и у эмбриона на имплантацию и последующее фетальное развитие [46]. В других исследованиях доказана роль miR-425-5p в регуляции клеточной пролиферации, миграции и инвазии и выявлен повышенный уровень ее экспрессии в тканях различных опухолей, в том числе карциноме почек [47], гепатоклеточной карциноме [48], раке шейки матки [49], раке желудка [50]. miR-34c-5p, являясь ингибитором эпителиально-мезенхимального перехода, клеточной инвазии и миграции путем подавления уровня экспрессии гена-мишени Notch1, вовлечена в патогенез эндометриоза при снижении уровня ее экспрессии в эктопическом эндометрии [51]. Значительное подавление экспрессии miR-34c-5p в эутопическом эндометрии отмечают в период «окна имплантации», в течение которого резко повышается уровень белка гена-мишени GAS, необходимого для имплантации эмбриона; в то время как дефицит экспрессии GAS1 обнаружен в эндометрии женщин с неоднократными неудачами имплантации в сравнении с рецептивным эндометрием [52]. miR-1180-3p рассматривают в качестве супрессора многих эпителиальных опухолей, в том числе рака мочевого пузыря [53], рака поджелудочной железы [54], колоректального рака [55], и повышенный уровень экспрессии miR-1180-3p обусловливает снижение пролиферации и подвижности раковых клеток путем воздействия на ген-мишень COL12A1 [55].

Заключение

Индивидуализация подхода в определении «окна имплантации» при проведении ЭКО на фоне ЦГТ в соответствии с уровнем мнкРНК в МЖ может играть важную роль в повышении эффективности программ ВРТ. Ограничением использования в клинической практике модели логистической регрессии, построенной в настоящем исследовании на обучающей выборке пациенток, является отсутствие информации об имплантационном потенциале переносимой в полость матки бластоцисты. Отрицательный исход программы ВРТ может быть обусловлен не отсутствием рецептивного эндометрия, а сниженным качеством самого эмбриона. Одновременное определение рецептивности эндометрия и имплантационного потенциала эмбриона по профилю мнкРНК приведет к снижению процента ложноотрицательных результатов и улучшению качества модели. Точность построенных моделей необходимо проверять на независимой тестовой выборке.

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

  1. Gerrits T., Van Rooij F., Esho T., Ndegwa W., Goossens J., Bilajbegovic A. et al. Infertility in the Global South: Raising awareness and generating insights for policy and practice. Facts Views Vis Obgyn. 2017; 9(1): 39-44.
  2. Sun H., Gong T.T., Jiang Y.T., Zhang S., Zhao Y. H.,Wu Q.J. Global, regional, and national prevalence and disability-adjusted life-years for infertility in 195 countries and territories, 1990-2017: results from a global burden of disease study. Aging. 2017; 11(23): 10952-91. https://dx.doi.org/10.18632/aging.102497.
  3. Mani S., Mainigi M. Embryo culture conditions and the epigenome. Semin. Reprod. Med. 2018; 36(3-04): 211-20. https://dx.doi.org/10.1055/s-0038-1675777.
  4. Simopoulou M., Sfakianoudis K., Rapani A., Giannelou P., Anifandis G., Bolaris S. et al. Considerations regarding embryo culture conditions: from media to epigenetics. In Vivo. 2018; 32(3): 451-60. https://dx.doi.org/10.21873/invivo.11261.
  5. Kirkegaard K., Agerholm I.E., Ingerslev H.J. Time-lapse monitoring as a tool for clinical embryo assessment. Hum. Reprod. 2012; 27(5): 1277-85.https://dx.doi.org/10.1093/humrep/des079.
  6. Gardner D.K., Balaban B. Assessment of human embryo development using morphological criteria in an era of time-lapse, algorithms and 'OMICS': is looking good still important? Mol. Hum. Reprod. 2016; 22(10): 704-18.https://dx.doi.org/10.1093/molehr/gaw057.
  7. Timofeeva A.V., Fedorov I.S., Shamina M. A., Chagovets V.V., Makarova N.P., Kalinina E.A. et al. Clinical relevance of secreted small noncoding RNAs in an embryo implantation potential prediction at morula and blastocyst development stages. life (Basel). 2021; 11(12): 1328. https://dx.doi.org/10.3390/life11121328.
  8. Timofeeva A.V., Drapkina Y.S., Fedorov I.S., Chagovets V.V., Makarova N.P., Shamina M.A. et al. Small noncoding RNA signatures for determining the developmental potential of an embryo at the morula ssage. Int. J. Mol. Sci. 2020; 21(24): 9399. https://dx.doi.org/10.3390/ijms21249399.
  9. Simopoulou M., Sfakianoudis K., Tsioulou P., Rapani A., Maziotis E., Giannelou P. et al. Should the flexibility enabled by performing a day-4 embryo transfer remain as a valid option in the IVF laboratory? A systematic review and network meta-analysis. J. Assist. Reprod. Genet. 2019; 36(6): 1049-61.https://dx.doi.org/10.1007/s10815-019-01475-0.
  10. Li Y.X., Wang J., Sun T.Z., Lv M.Q., Ge P., Li H.N. et al. Pregnancy outcomes after day 5 versus day 6 blastocyst-stage embryo transfer: A systematic review and meta-analysis. J. Obstet. Gynaecol. Res. 2020; 46(4): 595-605.https://dx.doi.org/10.1111/jog.14188.
  11. Luddi A., Pavone V., Semplici B., Governini L., Criscuoli M., Paccagnini E. et al. Organoids of human endometrium: a powerful in vitro model for the endometrium-embryo cross-talk at the implantation site. Cells. 2020 ;9(5): 1121. https://dx.doi.org/10.3390/cells9051121.
  12. Craciunas L., Gallos I., Chu J., Bourne T., Quenby S., Brosens J.J. et al. Conventional and modern markers of endometrial receptivity: a systematic review and meta-analysis. Hum. Reprod. Update. 2019; 25(2): 202-23.https://dx.doi.org/10.1093/humupd/dmy044.
  13. Massimiani M., Lacconi V., La Civita F., Ticconi C., Rago R. et al. Molecular signaling regulating endometrium-blastocyst crosstalk. Int. J. Mo.l Sci. 2019; 21(1): 23. https:/dx./doi.org/10.3390/ijms21010023.
  14. Kieu V., Lantsberg D., Mizrachi Y., Stern C., Polyakov A., Teh W.T. A survey study of endometrial receptivity tests and immunological treatments in in vitro fertilisation (IVF). Aust. N. Z. J. Obstet. Gynaecol. 2022; 62(2): 306-11.https://dx.doi.org/10.1111/ajo.13466.
  15. Sebastian-Leon P., Garrido N., Remohí J., Pellicer A., Diaz-Gimeno P. Asynchronous and pathological windows of implantation: two causes of recurrent implantation failure. Hum. Reprod. 2018; 33(4): 626-35.https://dx.doi.org/10.1093/humrep/dey023.
  16. Valdes C.T., Schutt A., Simon C. Implantation failure of endometrial origin: it is not pathology, but our failure to synchronize the developing embryo with a receptive endometrium. Fertil. Steril. 2017; 108(1): 15-8.https://dx.doi.org/10.1016/j.fertnstert.2017.05.033.
  17. Aghajanova L., Hamilton A.E., Giudice L.C. Uterine receptivity to human embryonic implantation: histology, biomarkers, and transcriptomics. Semin. Cell Dev. Biol. 2008; 19(2): 204-11. https://dx.doi.org/10.1016/j.semcdb.2007.10.008.
  18. Nejatbakhsh R., Kabir-Salmani M., Dimitriadis E., Hosseini A., Taheripanah R., Sadeghi Y. et al. Subcellular localization of L-selectin ligand in the endometrium implies a novel function for pinopodes in endometrial receptivity. Reprod. Biol. Endocrinol. 2012; 10: 46. Erratum in: Reprod. Biol. Endocrinol. 2021;19(1): 62. https://dx.doi.org/10.1186/1477-7827-10-46.
  19. Quinn K.E., Matson B.C., Wetendorf M., Caron K.M. Pinopodes: recent advancements, current perspectives, and future directions. Mol. Cell. Endocrinol. 2020; 501: 110644. https://dx.doi.org/10.1016/j.mce.2019.110644.
  20. Xu B., Sun X., Li L., Wu L., Zhang A., Feng Y. Pinopodes, leukemia inhibitory factor, integrin-β3, and mucin-1 expression in the peri-implantation endometrium of women with unexplained recurrent pregnancy loss. Fertil. Steril. 2012; 98(2): 389-95. https://dx.doi.org/10.1016/j.fertnstert.2012.04.032.
  21. Li F., Zhang M., Zhang Y., Liu T., Qu X. GnRH analogues may increase endometrial Hoxa10 promoter methylation and affect endometrial receptivity. Mol. Med. Rep. 2015; 11(1): 509-14. https://dx.doi.org/10.3892/mmr.2014.2680.
  22. Qiong Z., Jie H., Yonggang W., Bin X., Jing Z., Yanping L. Clinical validation of pinopode as a marker of endometrial receptivity: a randomized controlled trial. Fertil. Steril. 2017; 108(3): 513-7.e2. https://dx.doi.org/10.1016/j.fertnstert.2017.07.006.
  23. Jin X.Y., Zhao L.J., Luo D.H., Liu L., Dai Y.D., Hu X.X. et. al. Pinopode score around the time of implantation is predictive of successful implantation following frozen embryo transfer in hormone replacement cycles. Hum. Reprod. 2017; 32(12): 2394-403. https://dx.doi.org/10.1093/humrep/dex312.
  24. Hashimoto T., Koizumi M., Doshida M., Toya M., Sagara E., Oka N. et. al. Efficacy of the endometrial receptivity array for repeated implantation failure in Japan: a retrospective, two-centers study. Reprod. Med. Biol. 2017; 16(3):290-6. https://dx.doi.org/10.1002/rmb2.12041.
  25. Tan J., Kan A., Hitkari J., Taylor B., Tallon N., Warraich G. et. al. The role of the endometrial receptivity array (ERA) in patients who have failed euploid embryo transfers. J. Assist. Reprod. Genet. 2018; 35(4): 683-92.https://dx.doi.org/10.1007/s10815-017-1112-2.
  26. Altmäe S., Koel M., Võsa U., Adler P., Suhorutšenko M., Laisk-Podar T. et. al. Meta-signature of human endometrial receptivity: a meta-analysis and validation study of transcriptomic biomarkers. Sci. Rep. 2017; 7(1): 10077. https://dx.doi.org/10.1038/s41598-017-10098-3.
  27. van der Gaast M.H., Beier-Hellwig K., Fauser B.C., Beier H.M., Macklon N.S. Endometrial secretion aspiration prior to embryo transfer does not reduce implantation rates. Reprod. Biomed. Online. 2003; 7(1): 105-9.https://dx.doi.org/10.1016/s1472-6483(10)61737-3.
  28. Matorras R., Quevedo S., Corral B., Prieto B., Exposito A., Mendoza R. et. al. Proteomic pattern of implantative human endometrial fluid in in vitro fertilization cycles. Arch. Gynecol. Obstet. 2018; 297(6): 1577-86.https://dx.doi.org/10.1007/s00404-018-4753-1.
  29. Li T., Greenblatt E.M., Shin M.E., Brown T.J., Chan C. Cargo small non-coding RNAs of extracellular vesicles isolated from uterine fluid associate with endometrial receptivity and implantation success. Fertil. Steril. 2021; 115(5): 1327-36. https://dx.doi.org/10.1016/j.fertnstert.2020.10.046.
  30. Bergenheim S.J., Saupstad M., Pistoljevic N., Andersen A.N., Forman J.L., Løssl K. et al. Immediate versus postponed frozen embryo transfer after IVF/ICSI: a systematic review and meta-analysis. Hum. Reprod. Update. 2021; 27(4): 623-42. https://dx.doi.org/10.1093/humupd/dmab002.
  31. Langmead B., Trapnell C., Pop M., Salzberg S.L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009; 10: R25. https://dx.doi.org/10.1186/gb-2009-10-3-r25.
  32. 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.
  33. Team R.C. A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria. Available at:https://www.R-project.org Accessed 10.03. 2021.
  34. Team Rs. RStudio: Integrated Development for R. RStudio. Available at: http://www.rstudio.com/ Accessed 23.03.2021.
  35. Wang W., Vilella F., Alama P., Moreno I., Mignardi M., Isakova A. et. al. Single-cell transcriptomic atlas of the human endometrium during the menstrual cycle. Nat. Med. 2020; 26(10): 1644-53. https://dx.doi.org/10.1038/s41591-020-1040-z.
  36. Казачков Е.Л., Воропаева Е.Е., Казачкова Э.А., Затворницкая А.В., Дуб А.А., Мирошниченко Л.Е. Морфологическая характеристика эндометрия у пациенток с миомой матки и хроническим эндометритом при бесплодии. Архив патологии. 2019; 81(6): 41-8.
  37. Green S.B. How many subjects does it take to do a regression analysis. Multivariate Behav. Res. 1991; 26(3): 499-510. https://dx.doi.org/10.1207/s15327906mbr2603_7.
  38. Koel M., Krjutškov K., Saare M., Samuel K., Lubenets D., Katayama S. et. al. Human endometrial cell-type-specific RNA sequencing provides new insights into the embryo-endometrium interplay. Hum. Reprod. Open. 2022; 2022(4): hoac043. https://dx.doi.org/10.1093/hropen/hoac043.
  39. Evans J., Hutchison J., Salamonsen L.A., Greening D.W. Proteomic insights into endometrial receptivity and embryo-endometrial epithelium interaction for implantation reveal critical determinants of fertility. Proteomics. 2020; 20(1): e1900250. https://dx.doi.org/10.1002/pmic.201900250.
  40. Ruane P.T., Garner T., Parsons L., Babbington P.A., Wangsaputra I., Kimber S.J. et. al. Trophectoderm differentiation to invasive syncytiotrophoblast is promoted by endometrial epithelial cells during human embryo implantation. Hum. Reprod. 2022; 37(4): 777-92. https://dx.doi.org/10.1093/humrep/deac008.
  41. Bojić-Trbojević Ž., Jovanović Krivokuća M., Vilotić A., Kolundžić N., Stefanoska I., Zetterberg F. et. al. Human trophoblast requires galectin-3 for cell migration and invasion. Sci. Rep. 2019; 9(1): 2136. https://dx.doi.org/10.1038/s41598-018-38374-w.
  42. Vilella F., Moreno-Moya J.M., Balaguer N., Grasso A., Herrero M., Martínez S. et. al. Hsa-miR-30d, secreted by the human endometrium, is taken up by the pre-implantation embryo and might modify its transcriptome. Development. 2015; 142(18): 3210-21. https://dx.doi.org/10.1242/dev.124289.
  43. Sha A.G., Liu J.L., Jiang X.M., Ren J.Z., Ma C.H., Lei W. et. al. Genome-wide identification of micro-ribonucleic acids associated with human endometrial receptivity in natural and stimulated cycles by deep sequencing. Fertil. Steril. 2011; 96(1): 150-5.e5. https://dx.doi.org/10.1016/j.fertnstert.2011.04.072.
  44. Kuokkanen S., Chen B., Ojalvo L., Benard L., Santoro N., Pollard J.W. Genomic profiling of microRNAs and messenger RNAs reveals hormonal regulation in microRNA expression in human endometrium. Biol. Reprod. 2010; 82(4): 791-801. https://dx.doi.org/10.1095/biolreprod.109.081059.
  45. Zhao Y., He D., Zeng H., Luo J., Yang S., Chen J. et. al. Expression and significance of miR-30d-5p and SOCS1 in patients with recurrent implantation failure during implantation window. Reprod. Biol. Endocrinol. 2021; 19(1): 138. https://dx.doi.org/10.1186/s12958-021-00820-2.
  46. Balaguer N., Moreno I., Herrero M., Gonzáléz-Monfort M., Vilella F., Simón C. MicroRNA-30d deficiency during preconception affects endometrial receptivity by decreasing implantation rates and impairing fetal growth. Am. J. Obstet. Gynecol. 2019; 221(1): 46.e1-46.e16. https://dx.doi.org/10.1016/j.ajog.2019.02.047.
  47. Quan J., Li Y., Pan X., Lai Y., He T., Lin C. et. al. Oncogenic miR-425-5p is associated with cellular migration, proliferation and apoptosis in renal cell carcinoma. Oncol. Lett. 2018; 16(2): 2175-84. https://dx.doi.org/10.3892/ol.2018.8948.
  48. Fang F., Song T., Zhang T., Cui Y., Zhang G., Xiong Q. MiR-425-5p promotes invasion and metastasis of hepatocellular carcinoma cells through SCAI-mediated dysregulation of multiple signaling pathways. Oncotarget. 2017; 8(19): 31745-57. https://dx.doi.org/10.18632/oncotarget.15958.
  49. Sun L., Jiang R., Li J., Wang B., Ma C., Lv Y. et. al. MicoRNA-425-5p is a potential prognostic biomarker for cervical cancer. Ann. Clin. Biochem. 2017; 54(1): 127-33. https://dx.doi.org/10.1177/0004563216649377.
  50. Zhang Z., Li Y., Fan L., Zhao Q., Tan B., Li Z. et. al. microRNA-425-5pis upregulated in human gastric cancer and contributes to invasion and metastasis in vitro and in vivo. Exp. Ther. Med. 2015; 9(5): 1617-22.https://dx.doi.org/10.3892/etm.2015.2318.
  51. Luo Y., Wang D., Chen S., Yang, Q. The role of miR-34c-5p/Notch in epithelial-mesenchymal transition (EMT) in endometriosis. Cell. Signal. 2020; 72: 109666. https://dx.doi.org/10.1016/j.cellsig.2020.109666.
  52. Tan Q., Shi S., Liang J., Zhang X., Cao D., Wang Z. MicroRNAs in small extracellular vesicles indicate successful embryo implantation during early pregnancy. Cells. 2020; 9(3): 645. https://dx.doi.org/10.3390/cells9030645.
  53. Ge Q., Wang C., Chen Z., Li F., Hu J., Ye Z. The suppressive effects of miR-1180-5p on the proliferation and tumorigenicity of bladder cancer cells. Histol. Histopathol. 2017; 32(1): 77-86. https://dx.doi.org/10.14670/HH-11-772.
  54. Gu L., Zhang J., Shi M., Peng C. The effects of miRNA-1180 on suppression of pancreatic cancer. Am. J. Transl. Res. 2017; 9(6): 2798-806.
  55. Li C., Jin W., Zhang D., Tian S. Clinical significance of microRNA-1180-3p for colorectal cancer and effect of its alteration on cell function. Bioengineered. 2021; 12(2): 10491-500. https://dx.doi.org/10.1080/21655979.2021.1997694.

Поступила 03.05.2023

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

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

Тимофеева Анжелика Владимировна, к.б.н., заведующая лабораторией прикладной транскриптомики отдела системной биологии в репродукции, Национальный медицинский исследовательский центр акушерства, гинекологии и перинатологии имени академика В.И. Кулакова Минздрава России, v_timofeeva@oparina4.ru,
https://orcid.org/0000-0003-2324-9653, 117997, Россия, Москва, ул. Академика Опарина, д. 4.
Федоров Иван Сергеевич, м.н.с. лаборатории прикладной транскриптомики отдела системной биологии в репродукции, Национальный медицинский исследовательский центр акушерства, гинекологии и перинатологии имени академика В.И. Кулакова Минздрава России, is_fedorov@oparina4.ru, https://orcid.org/0000-0002-2104-5887, 117997, Россия, Москва, ул. Академика Опарина, д. 4.
Гохберг Яэль Александровна, аспирант отделения вспомогательных технологий в лечении бесплодия им. профессора Б.В. Леонова, Национальный медицинский исследовательский центр акушерства, гинекологии и перинатологии имени академика В.И. Кулакова Минздрава России, dr.yaelgokhberg@gmail.com,
https://orcid.org/0000-0003-3637-6096, 117997, Россия, Москва, ул. Академика Опарина, д. 4.
Калинина Елена Анатольевна, д.м.н, профессор, заведующая отделением вспомогательных технологий в лечении бесплодия им. профессора Б.В. Леонова, Национальный медицинский исследовательский центр акушерства, гинекологии и перинатологии имени академика В.И. Кулакова Минздрава России,
e_kalinina@oparina4.ru, https://orcid.org/0000-0002-8922-2878, 117997, Россия, Москва, ул. Академика Опарина, 4.

Также по теме

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