Полнотранскриптомный анализ событий альтернативного сплайсинга РНК при выраженной и умеренной формах задержки роста плода
Обоснование. Задержка роста плода (ЗРП) остается одной из ведущих причин перинатальной заболеваемости, а также серьезным фактором риска неблагоприятных исходов для здоровья ребенка в долгосрочной перспективе, включая повышенную вероятность неврологических, метаболических и сердечно-сосудистых заболеваний. Несмотря на высокий интерес к этой проблеме, молекулярные механизмы, лежащие в основе ЗРП, изучены недостаточно. Особенно мало данных о роли посттранскрипционной регуляции, в частности альтернативного сплайсинга (АС), в развитии этого заболевания, хотя именно этот процесс детерминирует разнообразие изоформ РНК и широту диапазона функциональных свойств клеток, определяя способности адаптироваться к патологическим воздействиям и степень подверженности к заболеваниям, включая акушерские патологии. Цель исследования: характеристика профилей АС децидуальных клеток (ДК) плаценты, определяющих тяжесть течения ЗРП. Материал и методы. Исследование выполнено на образцах плацентарной ткани пациенток с умеренной и выраженной формами ЗРП. Осуществлен полнотранскриптомный анализ ДК плаценты, которые были получены с использованием технологии лазерной микродиссекции препаратов тонких окрашенных срезов. Полнотранскриптомное секвенирование рибонуклеиновой кислоты (РНК) выполнено с использованием SMARTer Stranded Total RNA-Seq kit v2 («Takara BIO»). Анализ событий АС проведен с помощью пакета «MAJIQ» в операционной системе Linux. Для оценки вероятности дифференциальных событий применялись стандартные параметры MAJIQ, включая порог значимости p < 0,05 и минимальное изменение индекса сплайсинга |ΔPSI| > 0,20 между группами, что соответствует рекомендациям разработчиков алгоритма для выявления биологически значимых событий. Результаты. В изученных выборках было обнаружено 13 688 событий АС в 4 002 генах, экспрессирующихся в ДК. Более 52% событий являются идентичными для обеих групп. Как среди аннотированных, так и de novo событий при выраженной форме ЗРП выявлено статистически значимое снижение частоты альтернативного первого экзона (χ2 = 8,48; p = 0,004; χ2 = 6,15; p = 0,014 соответственно). Альтернативно сплайсированные гены, специфичные для выраженной ЗРП, вовлечены в следующие биологические процессы: каталитическая активность, действующая на нуклеиновые кислоты (pFDR = 0,020); регуляция активности ГТФаз (pFDR = 0,021); регуляция активности нуклеозидтрифосфатазы (pFDR = 0,021) и активность пептидной N-ацетилтрансферазы (pFDR = 0,028). При сравнении групп с умеренной и выраженной ЗРП идентифицированы 84 дифференциально сплайсированных гена (0,200 < deltaPSI < 0,648; p < 0,05), статистически значимо ассоциированных с такими биологическими и сигнальными путями, как различные виды репарации ДНК, сигнальный путь лиганд-управляемых ионных каналов, везикулярный транспорт к плазматической мембране, регуляция метаболизма мРНК, организация пероксисом, организация лизосом, морфогенез, сигнальный путь SMAD, метаболизм серы и АТФ-зависимое ремоделирование хроматина. Заключение. Полученные данные указывают на существование определенного набора молекулярных изменений на уровне АС, характерных для ЗРП независимо от степени ее тяжести. Паттерны АС, специфичные для выраженной ЗРП, ассоциированы с нарушением базовых регуляторных систем клетки. Результаты функциональной аннотации дифференциально сплайсированных генов свидетельствуют о влиянии АС на посттранскрипционный контроль, клеточную архитектуру и межклеточные сигнальные взаимодействия при выраженной ЗРП.