《Middle East Fertility Society Journal》：Elucidating the role of transcription factors in molecular pathways underlying infertility in endometriosis: a bioinformatics approach
编辑推荐：
为阐明子宫内膜异位症导致不孕的分子机制，本研究通过整合分析多个人群的转录组数据，鉴定出20个保守的核心基因（Hub genes），并揭示其中12个为转录因子。研究发现这些基因在转录调控、发育生物学通路（如WNT信号）及上皮屏障功能中发挥关键作用，为探索其作为潜在生物标志物或治疗靶点提供了新方向。
子宫内膜异位症，这个困扰着全球约10%育龄期女性的慢性炎症性疾病，不仅是长期剧烈盆腔痛的元凶，更是导致女性不孕的重要原因。虽然经血逆流理论等经典假说被广泛接受，但其导致不孕的具体分子机制仍如雾里看花。目前，子宫内膜异位症的诊断依赖于有创的腹腔镜检查，确诊平均被延迟4到11年之久，迫切需要非侵入性的生物标志物。现有的标志物如CA-125，其特异性与敏感性尚不足以独当一面。因此，寻找稳定、可靠、跨越人群差异的分子靶点，成为解开子宫内膜异位症导致不孕谜团的关键突破口，也是开发新疗法的希望所在。一篇发表在《Middle East Fertility Society Journal》上的研究，正是通过生物信息学方法，试图从纷繁复杂的基因数据中，揪出那批“核心”分子。
研究人员为攻克这一难题，巧妙地整合了来自不同人群的四个独立基因芯片数据集（GSE7305、GSE120103、GSE23339、GSE51981）。他们利用Limma包进行差异表达基因（DEGs）分析，随后在不同数据集间取交集，以找到不受人群异质性影响的、最为保守的差异基因。接着，利用STRING数据库构建蛋白-蛋白相互作用（PPI）网络，并在Cytoscape软件中运用cytoHubba插件的多种算法（MCC、Degree、EPC、DMNC）筛选出最关键的“核心基因”。最后，通过功能富集分析（GO和通路分析）以及GeneMANIA网络分析，深入探索这些核心基因的生物学功能和相互作用。
差异基因鉴定与交叉分析
研究首先分析了四个数据集，分别识别出不同数量的差异表达基因。通过绘制韦恩图并取交集，最终在所有四个数据集中鉴定出203个稳定的重叠差异表达基因，包括7个共同上调基因和196个共同下调基因。这一分析过程可通过下图直观展示：
核心基因识别与转录因子占比
在203个重叠差异表达基因的基础上，研究构建了包含142个节点和307条边的蛋白相互作用网络。通过综合四种算法评分，最终锁定了一个包含20个基因的“共识核心基因”集，它们是：ESR1、EPCAM、MUC1、MMP9、CLDN7、MET、RAB25、PRSS8、PGR、DKK1、FOXA2、MSX1、HOXB4、HOXA9、HOXA10、PAX8、GATA2、SOX17、RUNX1和 MAL2。对这些基因的分析揭示了一个关键发现：在这20个核心基因中，有多达12个（占60%）是转录因子。下图展示了这20个核心基因的PPI网络，节点大小和颜色反映了其连接度：
功能富集分析揭示关键通路与生物学过程
研究对20个核心基因进行了功能富集分析。京都基因与基因组百科全书（KEGG）通路分析显示，这些基因显著富集在“癌症通路”、“癌症中的蛋白聚糖”和“紧密连接”等通路。而Reactome通路分析则凸显了它们在“核受体信号传导”、“RUNX1调控的转录”、“WNT信号”、“原肠胚形成”和“雌激素依赖的基因表达”等发育生物学通路中的重要作用。基因本体（Gene Ontology, GO）分析进一步细化了其功能。在生物学过程层面，这些基因主要与“组织发育”、“模式规范”、“胚胎形态发生”以及“RNA聚合酶II介导的DNA模板转录正调控”相关。在分子功能层面，它们主要作为“转录调控区核酸结合蛋白”和“DNA结合转录激活因子”。在细胞组分层面，它们富集于“转录调节复合体”、“蛋白质-DNA复合体”、“染色质”以及“紧密连接”等结构。下图汇总了通路和GO富集分析的主要发现：
核心基因的功能相互作用网络
利用GeneMANIA工具构建的相互作用网络显示，这些核心基因之间存在高度的共表达和物理相互作用，特别是涉及FOXA2、SOX17、PAX8、HOXA9、HOXA10、HOXB4、GATA2、RUNX1和DKK1等转录因子。这个网络在“模式规范”、“区域化”、“胚胎形态发生”和“胚胎器官发育”等发育过程中显著富集。下图直观展示了这个相互作用网络：
结论与讨论
本研究的整合性分析揭示，尽管研究人群存在异质性，但存在一套在子宫内膜异位症中持续发生改变的核心基因集合。其中转录因子占比高达60%，这强烈暗示转录失调是子宫内膜异位症发病的核心机制之一。这些转录因子，如FOXA2、SOX17、PAX8以及HOX家族成员（HOXA10、HOXA9、HOXB4），是调控子宫腺体发育、子宫内膜容受性和上皮稳态的关键“主控开关”。功能富集分析表明，这些核心基因协同作用于转录调控、组织发育和上皮屏障功能。例如，GO分析中“紧密连接”的富集提示，子宫内膜上皮屏障完整性的破坏可能促进了异位病灶的定植和侵袭。同时，通路分析中WNT信号、核受体（如ESR1和PGR）信号等发育相关通路的富集，将子宫内膜异位症与异常的发育编程和激素反应失衡联系起来。例如，已知HOXA10的表达异常会影响子宫内膜容受性，而孕激素抵抗可导致WNT/β-连环蛋白（WNT/β-catenin）信号通路异常激活，这都与本研究发现的信号通路相互印证。总而言之，该研究通过跨越人群的数据整合，筛选出一组在子宫内膜异位症中稳定失调的核心基因，其中转录因子扮演了枢纽角色。这些基因通过扰乱关键的转录网络、发育通路和上皮屏障功能，共同促成了子宫内膜异位症的发生发展及其相关的不孕结局。这些发现不仅深化了对疾病分子机制的理解，更重要的是，这20个核心基因，特别是那12个转录因子，为后续开发诊断生物标志物和新型靶向疗法（例如针对特定的转录因子或WNT信号通路）提供了极具潜力的候选靶点，为饱受子宫内膜异位症不孕困扰的女性带来了新的希望。
