Транскриптомный анализ здоровых пациентов и пациентов с ПЭ
Для анализа различий в экспрессии генов между контрольной группой (Normal) и пациентами с преэклампсией (PE) использовали пакет R limma для выявления дифференциально экспрессируемых генов в объединенных наборах данных. Всего было выявлено 89 дифференциально экспрессируемых генов (DEG) (p < 0.05 и |log2FC| > 0.5), из которых 69 генов имели повышенную экспрессию, а 20 — пониженную для GSE75010 (Рисунок 1A). Было обнаружено 2 097 дифференциально экспрессируемых генов, из которых 1 121 с повышенной экспрессией и 976 с пониженной экспрессией (Рисунок 1B). Для набора данных GSE60438 было выявлено 7 633 DEG, включая 3 371 ген с повышенной экспрессией и 4 262 гена с пониженной экспрессией (Рисунок 1C). Значительные различия в количестве DEG в трех наборах данных, вероятно, отражают различия в размере выборки, происхождении тканей (ворсины плаценты в сравнении с децидуальной оболочкой базиса) и классификации степени тяжести PE, что дополнительно подчеркивает ценность поиска пересечений между наборами данных для выявления стабильных общих DEG.
Функциональный анализ обогащения дифференциально экспрессируемых генов (ДЕГ)
Для систематической идентификации ключевых дифференциально экспрессируемых генов (DEG), связанных с преэклампсией (PE), был проведен анализ диаграмм Венна для DEG из наборов данных PE GSE75010, GSE10588 и GSE60438, в результате чего было выявлено 56 общих DEG (Рисунок 2A). Для этих 56 общих DEG была построена сеть белок-белковых взаимодействий (PPI) с использованием онлайн-платформы STRING (показатель взаимодействия ≥0.15), после чего 20 основных хаб-генов были определены на основе степени их связности в сети и ранжированы с помощью программного обеспечения Cytoscape (Рисунок 2B). На основе этих генов также был проведен анализ обогащения GO для изучения потенциальных функций в категориях биологических процессов (BP), клеточных компонентов (CC) и молекулярных функций (MF) (Рисунок 2C). DEG, вероятнее всего, были обогащены в процессах секреции гормонов и регуляции сигналинга трансмембранных серин/треониновых киназ-рецепторов. Что касается клеточной локализации, гены в основном были сосредоточены в фокальных адгезиях, контактах клетка-субстрат, а также в ряде генов, связанных с везикулами и лизосомами. Среди молекулярных функций обогащенными терминами оказались пролилгидроксилазная активность, гормональная активность и трансмембранный транспорт липид-связанных молекул. Эти результаты позволяют предположить, что DEG могут играть важную роль в сигналинге клеток плаценты, взаимодействии клеток с матриксом и поддержании функций, связанных с эндокринной системой. В то же время было обнаружено, что анализ обогащения путей KEGG показал наибольшее обогащение DEG во взаимодействиях цитокин-цитокиновый рецептор, взаимодействиях молекул клеточной адгезии, а также в сигнальных путях TGF-β и HIF-1 (Рисунок 2D). Обогащение определенными вирус-ассоциированными терминами (например, вирусы гепатита, вирус Эбола), вероятно, отражает неспецифическое перекрытие с общими иммунорегуляторными генами, а не прямую биологическую связь с PE; поэтому основной интерпретационный акцент сделан на путях, связанных с воспалением, адгезией и гипоксией. Эти пути связаны с воспалительным ответом, адгезивно-опосредованной передачей сигналов и сигналингом, связанным с гипоксией, и могут влиять на прогрессирование PE, изменяя иммуновоспалительный статус и межклеточную коммуникацию в микроокружении плаценты.
Машинное обучение идентифицирует SASH1 как ключевой ген при преэклампсии
Для идентификации ключевых прогностических генов для отбора признаков были применены четыре метода машинного обучения. Регрессия LASSO с перекрестной проверкой показала, что модель достигает минимальной ошибки перекрестной проверки при log(λ.min) = -3.8075. Гены с ненулевыми коэффициентами при данном значении λ были выбраны в качестве ключевых признаков, определенных с помощью LASSO (Рисунок 3A–B). Алгоритм Boruta оценивает значимость признаков путем сравнения с теневыми признаками (shadow features). Гены, значения значимости которых превышали shadowMax, были помечены как Confirmed («Подтвержденные») и рассмотрены как ключевые гены-кандидаты, тогда как гены со значениями значимости ниже shadowMin были помечены как Rejected («Отклоненные») и признаны не вносящими существенного вклада. Гены между этими двумя порогами были помечены как Tentative («Предположительные»), что указывает на потенциальную, но неопределенную значимость (Рисунок 3C). Анализ SVM-RFE продемонстрировал, что при сохранении 40–50 признаков точность перекрестной проверки достигала самого высокого или близкого к нему уровня, что свидетельствует об оптимальной прогностической способности при таком масштабе признаков (Рисунок 3D). Умеренная точность перекрестной проверки (~0.6) указывает на то, что ни один алгоритм не показал оптимальных результатов для данного набора данных; поэтому для повышения надежности отбора признаков была принята стратегия пересечения результатов нескольких алгоритмов. XGBoost количественно определил относительную значимость генов с использованием значений Gain, выявив таким образом признаки с наибольшим вкладом в прогноз модели (Рисунок 3E). Анализ пересечений по диаграмме Венна результатов четырех методов выявил SASH1 как единственный перекрывающийся ген. Несмотря на то, что алгоритмом Boruta ген SASH1 был классифицирован как Tentative, его последовательный выбор тремя другими алгоритмами подтверждает его назначение в качестве консенсусного гена-кандидата (Рисунок 3F).
Одноклеточный анализ
Были исследованы клеточная гетерогенность и молекулярные процессы в тканях плаценты при преэклампсии (PE). Клетки проходили фильтрацию по качеству на основе количества детектированных генов, количества прочтений РНК и процента экспрессии митохондриальных генов. Низкокачественные клетки (менее 200 детектированных генов или более 20% митохондриального генетического контента) были удалены для получения высококачественной одноклеточной выборки для последующего анализа. Метрики контроля качества и идентификация высоковариабельных генов представлены на Supplementary Figure 1A–B. После нормализации и определения высоковариабельных генов был проведен анализ главных компонент (PCA) с использованием 2000 наиболее вариабельных генов; значимые главные компоненты были определены с помощью анализа JackStraw и анализа ElbowPlot (Supplementary Figure 2A–B).
На основе главных компонентов был использован метод однообразного аппроксимированного многообразия и проекции (UMAP) для снижения размерности и проведения неконтролируемой кластеризации, в результате чего все клетки были разделены на 22 кластера, что подчеркивает высокую клеточную гетерогенность плацентарной ткани (Рисунок 4A). Путем изучения экспрессии специфических для кластеров маркерных генов и визуализации дифференциально экспрессируемых генов (DEGs) на тепловой карте (Рисунок 4B), каждый кластер был аннотирован на основе канонического маркера и снова спроецирован на вложение UMAP (Рисунок 4C). Таким образом были идентифицированы основные популяции клеток плаценты, включая B-клетки/плазматические клетки, эндотелиальные клетки, фибробласты/стромальные клетки, макрофаги/дендритные клетки, натуральные киллеры, T-клетки и клетки трофобласта. Сравнение контрольной группы и группы с преэклампсией (PE) показывает, что несколько типов клеток имеют значительное изменение относительного распределения в плацентах при PE, что указывает на ремоделирование микроокружения плаценты (Рисунок 4D).
Происхождение ткани и статус заболевания использовались для отслеживания клеток на одном и том же UMAP-вложении в зависимости от источника ткани (децидуальная оболочка по сравнению с ворсинами) и группы заболевания (контроль по сравнению с PE) (Рисунок 4E). Распределение количества клеток было частично разделено по происхождению ткани, в то время как клетки из разных групп заболеваний были сильно взаимосвязаны. Это позволяет предположить, что PE не вызывает существенной перестройки всей транскриптомной структуры, но может оказывать влияние на отдельные популяции клеток или на определенные молекулярные признаки. Экспрессия SASH1 была нанесена на UMAP для изучения ее распределения среди популяций клеток плаценты (Рисунок 4F–G). SASH1 демонстрирует выраженную клеточно-специфическую экспрессию с заметными различиями между популяциями децидуальных клеток и клеток ворсин. Статистическое сравнение экспрессии SASH1 между группами PE и контроля в различных типах клеток (критерий суммы рангов Вилкоксона) выявило значимые различия в нескольких популяциях клеток плаценты, при этом наиболее выраженное повышение наблюдалось в клетках трофобласта из образцов с PE (Рисунок 4H). Этот результат предоставляет доказательства на уровне отдельных клеток того, что дисрегуляция SASH1 особенно выражена в компартменте трофобласта, что обосновывает использование клеточной линии трофобласта HTR-8/SVneo для последующих функциональных экспериментов.
Нокдаун SASH1 ингибирует индуцированный гипоксией апоптоз трофобласта и воспалительные реакции
Для подтверждения результатов биоинформатического прогнозирования была создана модель клеток трофобласта, связанная с ПЭ, в условиях гипоксии (1% O2). Для оценки экспрессии SASH1 были проведены анализы методом кПЦР и вестерн-блоттинга (Рисунок 5A–B). По сравнению с клетками HTR-8/SVneo в условиях нормоксии, в клетках, подвергнутых воздействию гипоксии, наблюдалось значительное повышение уровней экспрессии мРНК и белка SASH1, что обеспечило предварительное экспериментальное подтверждение биоинформатических данных. Чтобы оценить влияние нокдауна SASH1 на патологические процессы, связанные с ПЭ, в клетках HTR-8/SVneo, подвергнутых гипоксии, была осуществлена сайленсинг-супрессия SASH1. Количественная ПЦР и вестерн-блоттинг (Рисунок 5C–D) показали, что экспрессия SASH1 была значительно снижена во всех группах sh-SASH1 по сравнению с группой sh-NC. Среди них sh-SASH1#3 продемонстрировала наивысшую эффективность нокдауна и, следовательно, была выбрана для последующих экспериментов. Окрашивание TUNEL и ИФА проводились в трех условиях: нормоксический контроль (Normal), гипоксия с отрицательным контролем shRNA (Hypoxia + sh-NC) и гипоксия с нокдауном SASH1 (Hypoxia + sh-SASH1). По сравнению с группой Normal, в группе Hypoxia + sh-NC наблюдалось значительное усиление апоптоза и повышение секреции провоспалительных цитокинов IL-1β, IL-6 и TNF-α, что подтвердило повреждение трофобласта, индуцированное гипоксией. Нокдаун SASH1 заметно ослабил эти эффекты, снизив уровень апоптоза (Рисунок 5E) и уровень цитокинов (Рисунок 5F) до значений, близких к группе Normal, что указывает на то, что нокдаун SASH1 может в значительной степени нивелировать индуцированный гипоксией апоптоз трофобласта и воспалительные реакции. Для дальнейшего установления причинно-следственной роли SASH1 были проведены эксперименты по усилению функции. кПЦР и вестерн-блоттинг подтвердили успешную гиперэкспрессию SASH1 в клетках HTR-8/SVneo (Рисунок 5G–H). В условиях нормоксии гиперэкспрессия SASH1 (OE-SASH1) значительно увеличила апоптоз трофобласта и секрецию провоспалительных цитокинов по сравнению с группой OE-NC, воспроизводя фенотип, индуцированный гипоксией (Рисунок 5I–J). В совокупности данные о потере и усилении функции доказывают, что SASH1 является и необходимым, и достаточным фактором для запуска апоптоза трофобласта и воспалительных реакций, что подтверждает его причинную роль в патологии, связанной с ПЭ.
P-EXOS подавляет индуцированный гипоксией апоптоз трофобласта и воспалительные реакции путем модуляции SASH1
Считается, что P-EXOS обладают потенциальной терапевтической ценностью при ПЭ. В данном исследовании биоинформатический анализ показал, что SASH1 связан с сигнальными путями, относящимися к экзосомам (согласно анализу обогащения GO и KEGG, Рисунок 2C–D), что позволило нам выдвинуть гипотезу о том, что P-EXOS могут оказывать защитный эффект при ПЭ путем регуляции SASH1. Перед функциональными исследованиями была проведена характеристика изолированных P-EXOS. Вестерн-блот-анализ подтвердил наличие маркерных белков экзосом PLAP, CD63 и TSG101 при отсутствии детектируемой экспрессии маркера аппарата Гольджи GM130, что свидетельствует о высокой чистоте изолированных везикул. Просвечивающая электронная микроскопия дополнительно выявила типичные чашеобразные везикулярные структуры (Рисунок 6A–B). Чтобы проверить, могут ли P-EXOS эффективно интернализироваться клетками-мишенями, было проведено иммунофлуоресцентное окрашивание для оценки клеточного захвата P-EXOS. По сравнению с контрольной группой, клетки HTR-8/SVneo, обработанные P-EXOS, демонстрировали выраженные сигналы колокализации на совмещенных изображениях, что указывает на эффективную интернализацию P-EXOS клетками HTR-8/SVneo (Рисунок 6C).
Впоследствии клетки HTR-8/SVneo подвергали совместному культивированию с P-EXOS в условиях гипоксии. Как количественная ПЦР, так и вестерн-блоттинг показали выраженное снижение уровней мРНК и белка SASH1 в группе, обработанной P-EXOS, по сравнению с контролем (Рисунок 6D–E), что стало первичным доказательством регуляторного взаимодействия между P-EXOS и SASH1. В соответствии с этими молекулярными изменениями, окрашивание методом TUNEL выявило существенное ослабление индуцированного гипоксией апоптоза в клетках HTR-8/SVneo после обработки P-EXOS (Рисунок 6F). В то же время анализ методом ELISA продемонстрировал заметное снижение секреции провоспалительных цитокинов, включая IL-6, IL-1β и TNF-α (Рисунок 6G).
Чтобы определить, зависят ли защитные эффекты P-EXOS механистически от супрессии SASH1, в условиях гипоксии был проведен эксперимент по рескью (восстановлению функции). Клетки с гиперэкспрессией SASH1 (OE-SASH1 + P-EXOS) сравнивали с клетками отрицательного контроля (OE-NC + P-EXOS). Гиперэкспрессия SASH1 значительно нивелировала защитные эффекты P-EXOS, о чем свидетельствовали усиление апоптоза (Рисунок 6H>) и повышенная секреция провоспалительных цитокинов (Рисунок 6I>) по сравнению с группой OE-NC + P-EXOS. Эти результаты демонстрируют, что защитный эффект P-EXOS опосредован конкретно через супрессию SASH1, а не через параллельные пути, что подчеркивает терапевтический потенциал P-EXOS при преэклампсии (PE).
ДОСТУПНОСТЬ ДАННЫХ:
Наборы данных, подтверждающие результаты данного исследования, находятся в открытом доступе и были получены из базы данных Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/) под номерами доступа GSE75010, GSE10588, GSE60438 и GSE183338. Аннотирование типов клеток проводилось с использованием базы данных Cell Classification Database (https://ngdc.cncb.ac.cn/celltaxonomy/). Необрезанные изображения вестерн-блоттинга и все остальные необработанные экспериментальные данные, подтверждающие результаты данного исследования, представлены в папке Supplementary Raw Data.

Рисунок 1: Идентификация дифференциально экспрессируемых генов (ДЭГ) в тканях плаценты и децидуальной оболочке при преэклампсии (ПЭ) в трех независимых наборах данных. (A) Диаграмма вулкан ДЭГ между 157 образцами плаценты при ПЭ и 173 образцами плаценты без ПЭ из объединенного набора данных GSE75010. (B) Диаграмма вулкан ДЭГ между образцами плаценты при тяжелой ПЭ (n = 17) и нормальными образцами плаценты (n = 26) из набора данных GSE10588. (C) Диаграмма вулкан ДЭГ между образцами децидуальной оболочки при ПЭ и контрольными образцами децидуальной оболочки при нормальном артериальном давлении из набора данных GSE60438. ДЭГ были идентифицированы при p < 0.05 и |log2FC| > 0.5. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 2: Функциональный обогатительный анализ и анализ PPI общих DEG, связанных с преэклампсией. (A) Диаграмма Венна DEG, выявленных из наборов данных GSE75010, GSE10588 и GSE60438, связанных с PE. (B) Сеть PPI 20 основных хабовых генов, отобранных из 56 общих DEG на основе степени связности и ранжированных с помощью программного обеспечения Cytoscape (показатель взаимодействия STRING ≥0.15). (C) Результаты анализа обогащения GO и KEGG для общих DEG, представленные в виде столбчатых диаграмм. (D) Результаты анализа обогащения GO и KEGG для общих DEG, представленные в виде хордовых диаграмм. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 3: Скрининг потенциальных хаб-генов, связанных с преэклампсией, на основе машинного обучения. (A) Регрессионный анализ LASSO, демонстрирующий профили коэффициентов потенциальных генов. (B) Выбор оптимального параметра регуляризации (λ) при помощи регрессии LASSO с использованием кросс-валидации. (C) Алгоритм Boruta. (D) Анализ SVM-RFE. (E) Алгоритм XGBoost. (F) Диаграмма Венна, показывающая пересечение потенциальных генов, выявленных алгоритмами LASSO, Boruta, SVM-RFE и XGBoost. Сокращения: LASSO = оператор наименьшего абсолютного сжатия и отбора; SVM-RFE = рекурсивное исключение признаков методом опорных векторов. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 4: Транскриптомный анализ отдельных клеток выявляет специфичные для типов клеток паттерны экспрессии SASH1 при преэклампсии. (A) Кластеризация UMAP набора данных отдельных клеток. (B) Анализ аннотации пяти основных маркеров набора данных транскриптома отдельных клеток при ПЭ GSE183338. (C) Аннотация субпопуляций клеток. (D) Столбчатая диаграмма, показывающая распределение пропорций клеток в контрольных образцах и образцах с ПЭ. (E) График UMAP, раскрашенный по происхождению ткани (децидуальная оболочка и ворсины) и статусу заболевания (контроль и ПЭ). (F) График UMAP, показывающий экспрессию SASH1 в популяциях клеток плаценты, сгруппированных по происхождению ткани (децидуальная оболочка и ворсины). (G) График UMAP, показывающий экспрессию SASH1 в популяциях клеток плаценты, сгруппированных по статусу заболевания (контроль и ПЭ). (H) Скрипичный график, показывающий распределение экспрессии SASH1 в различных типах клеток плаценты. ns p > 0.05; *p < 0.05; **p < 0.01; ***p < 0.001. Сокращения: UMAP = равномерное многомерное аппроксимирование и проекция; ns = статистически не значимо. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 5: SASH1 регулирует индуцированный гипоксией апоптоз и воспаление трофобласта. (A) Анализ экспрессии мРНК SASH1 методом кПЦР и денситометрическое определение уровней белка SASH1 в клетках HTR-8/SVneo в нормальных условиях и при гипоксии. (B) Репрезентативные изображения вестерн-блоттинга, соответствующие денситометрическим данным, показанным на (A). (C) Анализ экспрессии мРНК SASH1 методом кПЦР и денситометрическое определение уровней белка SASH1 в клетках HTR-8/SVneo, подвергнутых воздействию гипоксии и трансфецированных конструктами sh-NC или sh-SASH1 (sh-SASH1#1, sh-SASH1#2 и sh-SASH1#3). (D) Репрезентативные изображения вестерн-блоттинга, соответствующие денситометрическим данным, показанным на (C). (E) Окрашивание методом TUNEL клеток HTR-8/SVneo в группах Normal, Hypoxia + sh-NC и Hypoxia + sh-SASH1(#3). Масштабная линейка = 20 µm. (F) Анализ провоспалительных цитокинов (IL-1β, IL-6 и TNF-α) методом ИФА в группах Normal, Hypoxia + sh-NC и Hypoxia + sh-SASH1(#3). (G) Анализ экспрессии мРНК SASH1 методом кПЦР и денситометрическое определение уровней белка SASH1 в клетках HTR-8/SVneo, трансфецированных конструктами OE-NC или OE-SASH1. (H) Репрезентативные изображения вестерн-блоттинга, соответствующие денситометрическим данным, показанным на (G). (I) Окрашивание методом TUNEL клеток HTR-8/SVneo в группах Normal + OE-NC и Normal + OE-SASH1. Масштабная линейка = 20 µm. (J) Анализ провоспалительных цитокинов (IL-1β, IL-6 и TNF-α) методом ИФА в группах Normal + OE-NC и Normal + OE-SASH1. *p < 0.05; **p < 0.01; ***p < 0.001. Данные представлены как среднее значение ± SD по трем независимым биологическим повторностям (n = 3). Нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 6: P-EXOS облегчают вызванное гипоксией повреждение трофобласта путем подавления SASH1. (A) Вестерн-блот анализ маркерных белков P-EXOS PLAP, CD63 и TSG101, где GM130 служит отрицательным контролем. (B) Анализ морфологии P-EXOS с помощью просвечивающей электронной микроскопии (ПЭМ). Масштабный отрезок = 1.0 µm. (C) Изображения, полученные с помощью конфокальной лазерной сканирующей микроскопии, демонстрирующие поглощение меченных PKH67 P-EXOS (зеленый) клетками HTR-8/SVneo (красный). Ядра были контр-окрашены DAPI (синий). Масштабный отрезок = 20 µm. (D) Анализ экспрессии мРНК SASH1 методом кПЦР и денситометрическая количественная оценка уровней белка SASH1 в клетках HTR-8/SVneo в группах Control и P-EXOS. (E) Репрезентативные изображения вестерн-блоттинга, соответствующие денситометрическим данным, показанным на (D). (F) TUNEL-окрашивание клеток HTR-8/SVneo в группах Control и P-EXOS. Масштабный отрезок = 20 µm. (G) ELISA-анализ воспалительных цитокинов (IL-1β, IL-6 и TNF-α) в клетках HTR-8/SVneo, подвергнутых воздействию гипоксии, в группах Control и P-EXOS. (H) TUNEL-окрашивание клеток HTR-8/SVneo в группах OE-NC + P-EXOS и OE-SASH1 + P-EXOS. Масштабный отрезок = 20 µm. (I) ELISA-анализ воспалительных цитокинов (IL-1β, IL-6 и TNF-α) в клетках HTR-8/SVneo, подвергнутых воздействию гипоксии, в группах OE-NC + P-EXOS и OE-SASH1 + P-EXOS. Сокращения: P-EXOS = экзосомы плацентарного происхождения; TEM = просвечивающая электронная микроскопия; PLAP = плацентарная щелочная фосфатаза; TSG101 = ген предрасположенности к опухолям 101. *p < 0.05; **p < 0.01. Данные представлены как среднее значение ±SD по трем независимым биологическим повторностям (n = 3). Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Дополнительный рисунок 1: Контроль качества данных scRNA-seq. (A) Скрипичный график уровней экспрессии генов после контроля качества. (B) График вариабельных признаков, демонстрирующий высоковариабельные гены, выявленные после контроля качества.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительный рисунок 2: Определение статистически значимых главных компонент для кластеризации единичных клеток. (A) Точечный график JackStraw, используемый для идентификации статистически значимых главных компонент. (B) Точечный график «локоть» (Elbow plot), показывающий долю дисперсии, объяснимую главными компонентами.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительная таблица 1: Мишеневые последовательности shRNA, использованные для нокдауна SASH1. Мишеневые последовательности shRNA (5′–3′) для отрицательного контроля (sh-NC) и трех конструктов, нацеленных на SASH1 (sh-SASH1#1, sh-SASH1#2 и sh-SASH1#3), использованных для сайленсинга гена.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительная таблица 2: Последовательности праймеров, использованных для количественной ПЦР. Последовательности прямых и обратных праймеров (5′–3′) для SASH1 и β-actin (внутренний контрольный образец), использованные в анализе qPCR.Пожалуйста, нажмите здесь, чтобы скачать этот файл.