妊娠期糖尿病应激下胎盘与胎儿基因表达的比较生物信息学分析

《In Silico Pharmacology》：A bioinformatic approach to investigate placental and fetal gene expression during diabetic stress in pregnancy

母体糖尿病会扰乱宫内环境，并对胎儿发育构成显著风险，然而胎盘和胎儿组织对糖尿病应激的分子反应仍未被完全阐明。研究人员采用一种比较生物信息学框架（comparative bioinformatic framework），处理两个公开可用的微阵列数据集（microarray datasets），分别代表胚胎第10.5天（embryonic day 10.5）采集的胎盘/蜕膜（placental/decidua）及胎儿组织。两个数据集均通过RMA归一化（Robust Multi-array Average normalization）、差异基因表达分析（differential gene expression analysis）、主成分分析（principal component analysis, PCA）、功能富集分析（functional enrichment analysis）以及R软件中的Pathview通路可视化（Pathview pathway visualization）进行处理。胎盘/蜕膜组织显示出相对受限的差异表达谱，仅有2个探针达到BH校正后p<0.05，1个探针同时满足BH校正后p<0.05且|log₂FC|>1。研究人员采用基于秩次的基因集富集分析（gene set enrichment analysis, GSEA）覆盖完整基因排序分布，以识别协调的通路水平变化。该分析发现涉及核糖体生物发生（ribosome biogenesis）、RNA加工与剪接（RNA processing and splicing）、细胞周期与DNA复制过程（cell-cycle and DNA-replication processes）以及氧化磷酸化（oxidative phosphorylation）的基因程序富集。相比之下，胎儿组织表现出更为广泛的差异表达，共有2,306个探针同时满足BH校正后p<0.05和|log₂FC|>1。KEGG过度代表性分析（KEGG over-representation analysis）鉴定出62条显著富集通路，包括与轴突导向（axon guidance）、Wnt信号（Wnt signaling）、黏着斑（focal adhesion）、MAPK信号（MAPK signaling）和PI3K-AKT信号（PI3K-AKT signaling）相关的通路。Pathview映射进一步显示，与胎盘/蜕膜组织中相对局部的变化相比，胎儿信号组分的转录水平改变分布更为广泛。这些发现表明，在糖尿病应激下，胎盘/蜕膜与胎儿组织表现出不同的转录反应。胎盘/蜕膜相对受限的差异表达谱与胎儿更广泛的转录改变相一致，支持母胎界面（maternal-fetal interface）可能调节糖尿病应激反应的模型。总的来说，这一比较分析鉴定了与糖尿病妊娠相关的协调分子过程，并为糖尿病妊娠的早期生物标志物开发提供了框架。

论文解读

妊娠期糖尿病是威胁母儿健康的重要妊娠并发症。胎盘作为母体与胎儿之间的关键界面，在营养和氧气交换、内分泌、代谢、血管及免疫调节中发挥核心作用。母体糖尿病会扰乱宫内环境，并可能通过改变胎盘结构和功能、影响胎儿发育编程及基因表达、改变母体代谢和炎症状态等机制对胎儿产生不良影响。然而，胎盘/蜕膜组织与胎儿组织在糖尿病应激下的分子反应差异尚未被充分揭示。本研究基于发育健康与疾病起源假说（Developmental Origins of Health and Disease, DOHaD），旨在利用公共转录组数据，比较糖尿病妊娠中胎盘/蜕膜和胎儿组织的基因表达变化，从而加深对母胎界面响应机制的理解。

研究人员使用两个公开的微阵列数据集：胎盘/蜕膜数据集GSE28277（n=10）和胎儿组织数据集GSE41095（n=6），二者均来自Salbaum等公开的链脲佐菌素（streptozotocin, STZ）诱导糖尿病小鼠模型，采样于胚胎第10.5天。所有数据使用Affymetrix Mouse Genome 430 2.0 Array平台。分析流程包括RMA归一化、limma差异表达分析、主成分分析、火山图、热图、功能富集以及Pathview通路可视化。结果显示，胎盘/蜕膜组织在探针水平仅有1个基因满足差异表达阈值，但基于完整基因排序的GSEA鉴定出RNA剪接、核糖体生物发生、细胞周期和DNA复制等多条通路发生协调变化；胎儿组织则呈现广泛转录改变，2,306个探针满足显著阈值，62条KEGG通路显著富集，其中涉及轴突导向、Wnt信号、黏着斑、MAPK和PI3K-AKT等信号通路。Pathview进一步显示，胎儿组织中MAPK、PI3K-AKT及神经退行相关通路的转录变化分布更为广泛。这些发现提示胎盘/蜕膜和胎儿组织在糖尿病应激下具有不同的转录反应，母胎界面可能对高血糖应激起到一定调节作用，而胎儿组织表现出更强的应激敏感性。该研究为糖尿病妊娠的早期生物标志物发现提供了比较转录组框架。

关键技术方法简要概括如下：研究人员从基因表达综合数据库（Gene Expression Omnibus, GEO）获取GSE28277（胎盘/蜕膜，10个样本）和GSE41095（胎儿组织，6个样本）两个公共芯片数据集，平台均为Affymetrix Mouse Genome 430 2.0 Array。数据采用R语言及Bioconductor包处理，首先进行RMA归一化，随后使用limma进行差异基因表达分析，运用PCA评估样本聚类。针对胎盘/蜕膜数据集差异基因极少的特点，采用基于秩次的基因集富集分析（GSEA）进行通路分析；针对胎儿数据集则使用KEGG过度代表性分析（ORA）鉴定富集通路，并通过Pathview将log₂FC映射至MAPK、PI3K-AKT及神经退行相关通路进行可视化。样本来源于STZ诱导的糖尿病小鼠模型E10.5组织。

研究结果部分包含以下主要内容。

1. 胎盘/蜕膜组织的转录响应：PCA显示对照组与糖尿病组明显分离，糖尿病样本离散度更大。差异表达分析发现45,101个探针中仅2个探针达到BH校正p<0.05，其中Atp1a3满足|log₂FC|>1（log₂FC=1.78，BH p=0.043），Rnf182未满足倍数变化阈值。热图显示多数排名靠前基因在糖尿病样本中表达升高。

2. 胎盘/蜕膜组织的GSEA通路变化：基于完整基因排序的GSEA显示，GO生物过程负富集于核糖体小亚基生物发生（NES=-2.68）、核糖体生物发生（NES=-2.61）、核糖核蛋白复合物生物发生（NES=-2.60）、RNA剪接（NES=-2.47）、rRNA加工（NES=-2.45）和mRNA加工（NES=-2.43）；正富集包括轴突导向、神经元投射导向、细胞外基质组织。KEGG GSEA负富集剪接体（NES=-2.58）和核糖体（NES=-2.54），正富集TNF信号、松弛素信号、胰岛素分泌和Wnt信号。

3. 胎儿组织的差异表达与KEGG富集：胎儿组织有15,844个探针达到BH p<0.05，其中2,306个探针满足双阈值（148个正变化、2,158个负变化），Prl8a2为最大负变化（log₂FC=-6.45）。KEGG过度代表性分析共获得62条显著富集通路，包括轴突导向（GeneRatio=53/921，BH p=4.66×10^-8）、Wnt信号（50/921，p=1.61×10^-7）、黏着斑（52/921，p=2.28×10^-6）以及MAPK信号（60/921，p=4.05×10^-4）、PI3K-AKT信号（64/921，p=2.75×10^-3）等。

4. 跨数据集比较：共享基因热图显示胎盘与胎儿样本整体分离，Lcn2、Kcnq2、Scd1等基因贡献显著。联合PCA进一步表明胎盘与胎儿数据集沿PC1明显分开，提示两个区室具有不同的转录组特征。Pathview映射显示，胎盘组织中PI3K-AKT、MAPK和神经退行通路的转录变化局限于少量节点；胎儿组织则广泛涉及上游受体、信号级联、下游效应器以及线粒体、自噬、凋亡等组分。在神经退行通路（mmu05022）中，胎儿组织多个线粒体组分下调，泛素-蛋白酶体和自噬相关基因上调，CASP3和BAX差异表达。

5. 表达分布差异：胎盘组织log₂FC中位数约0.0015（IQR：-0.077~0.079），范围-1.51~2.39；胎儿组织中位数约0.107（IQR：-0.128~0.284），范围-6.45~2.87。胎儿组织表现出更大的转录变异、更宽的分布范围和更强的负向尾部。

讨论部分指出，胎盘/蜕膜组织在探针水平差异表达受限，但GSEA揭示了RNA加工、翻译和细胞增殖相关通路的协调变化，说明有限差异基因数量并不代表该组织缺乏转录反应。胎儿组织广泛的转录改变及多条信号通路的富集，与糖尿病对增殖、分化、存活和代谢调控的干扰相一致。Pathview显示胎盘与胎儿组织共有线粒体功能障碍和蛋白质加工破坏的信号，但胎儿组织中信号更显著。这些结果与母胎界面可能部分调节母体糖尿病影响的模型一致，但研究尚不能直接证明缓冲效应。局限性包括样本量小、胎盘组织包含蜕膜成分、胎儿样本来自混合胚胎、缺乏独立验证等。研究结论为：这些发现表明，在糖尿病应激下，胎盘/蜕膜与胎儿组织表现出不同的转录反应；胎盘/蜕膜相对受限的差异表达谱与胎儿更广泛的转录改变相一致，支持母胎界面可能调节糖尿病应激反应的模型；总体而言，这一比较分析鉴定了与糖尿病妊娠相关的协调分子过程，并为糖尿病妊娠的早期生物标志物开发提供了框架。