《Journal of Genetic Engineering and Biotechnology》：Integrative dual-track transcriptomics reveals stage-specific coordination, regulatory divergence, and HSP90AA1-associated remodeling in human folliculogenesis
编辑推荐：
人卵泡发生依赖于卵母细胞与其周围颗粒细胞之间协调而非完全相同的发育重塑。然而这两大区室何时保持同步、何时分化为谱系特异性的调控状态仍不明确。本研究对人RNA-seq数据集GSE107746进行了整合双轨重新分析，将卵母细胞和颗粒细胞分别建模为卵泡发育进程中不同
人卵泡发生依赖于卵母细胞与其周围颗粒细胞之间协调而非完全相同的发育重塑。然而这两大区室何时保持同步、何时分化为谱系特异性的调控状态仍不明确。本研究对人RNA-seq数据集GSE107746进行了整合双轨重新分析，将卵母细胞和颗粒细胞分别建模为卵泡发育进程中不同但发育关联的区室。对148个测序文库的分析显示，区室身份是转录组变异的主导来源，支持进行区室感知的下游解读。在此框架下，卵母细胞遵循相对连续的发展轨迹，相邻阶段间已可见显著的转录重塑；而颗粒细胞早期阶段对比较弱，但在窦前卵泡及排卵前过渡期出现明显更强的晚期重组，尤其围绕窦状卵泡和预排卵转变。功能富集表明卵母细胞成熟与RNA加工及更广泛的基因组调控重塑相关，而颗粒细胞成熟则以渐进的线粒体和生物能量激活为主。共表达分析显示两区室均含有强晚期程序及反向的早期模块，表明成熟具有共享的系统层面架构，尽管这些程序的hub基因组成和生物学内容大体为区室特异性。机器学习验证强化了这种不对称性：卵母细胞阶段分类最好由紧凑的特征基因(eigengene)表征恢复，而颗粒细胞阶段区分则由更广的差异表达衍生特征集更好解析。在基因水平，HSP90AA1作为一个阶段关联标志物表现出区室特异性行为，在卵母细胞发育中呈进行性衰减、归属于选定的卵母细胞blue模块，并在颗粒细胞中显示出更急剧的过渡动态。综上，这些发现支持如下模型：人卵泡发生通过协调但非等同的转录组重塑推进，系统层面具有共享的发育逻辑，但生殖细胞与体细胞区室中具有不同的分子执行方式。
论文解读：整合双轨转录组学揭示人卵泡发生过程中阶段特异性协调、调控分化及HSP90AA1相关的重塑
该研究由Jing Wang、Zhimin Xu与Jun Ning完成，拟发表于《Journal of Genetic Engineering and Biotechnology》。人卵泡发生（human folliculogenesis）是卵母细胞与周围颗粒细胞（granulosa cells, GCs）通过双向通讯、代谢支持和信号交换共同成熟的长期发育过程，而非孤立的细胞自主程序。现有研究虽提供了阶段关联基因目录及单细胞图谱，但未阐明卵母细胞—颗粒细胞程序在哪些发育阶段保持协调、在哪些节点开始反映生殖系与体细胞区室间的调控递接（regulatory handoff）即出现调控分化（regulatory divergence）。公共数据集GSE107746提供了同一卵泡中卵母细胞与对应颗粒细胞跨五个卵泡阶段的配对RNA-seq数据，是探究此问题的理想资源。研究人员假设人卵泡发生并非卵母细胞与颗粒细胞的均匀平行成熟，而是存在阶段特异性区室间协调偏移，关键发育转折伴有调控分化及区室主导的分子程序。为此研究人员对GSE107746进行了整合双轨（dual-track）转录组再分析，旨在识别阶段特异性转录重编程、共享与区室特异性生物学过程及与发育协调或分化相关的候选调控模块。
主要关键技术方法：
研究人员采用GEO数据库人卵泡发生配对RNA-seq数据集GSE107746（148个双端测序文库，含77个卵母细胞及71个颗粒细胞文库，覆盖原始/初级/次级/窦/排卵前阶段），经fastp进行接头修剪与质控，Salmon v1.10.2基于GRCh38脱靶序列（decoy-aware）索引进行转录本定量，tximport转为基因水平计数矩阵后经edgeR过滤低表达基因并做TMM（trimmed mean of M values）标准化得24,684个基因的分析矩阵。分别按卵母细胞与颗粒细胞区室用limma-voom进行阶段感知差异表达分析（differential expression analysis, DE）；取显著且|log2FC|≥1基因用g:Profiler2做GO BP/KEGG/Reactome功能富集分析（functional enrichment analysis）；分区室构建加权基因共表达网络（weighted gene co-expression network analysis, WGCNA），软阈值选取满足无标度拓扑标准，动态剪切识别模块并计算模块特征基因（module eigengene, ME）与有序阶段进展的Pearson相关性以筛选阶段关联模块及hub基因；最后分区室用Elastic Net多项式逻辑回归结合重复分层5折交叉验证，比较特征基因、差异表达基因集、hub基因集等不同特征空间对发育阶段的多分类预测性能以验证区室特异性发育程序的可区分性。
研究结果
3.1. Final cohort structure and stage distribution across compartments
最终分析集含148个文库，卵母细胞区室77例（原始17、初级25、次级12、窦23），颗粒细胞区室71例（原始8、初级15、次级6、窦24、排卵前18）。原始/初级/次级/窦为两区室共有阶段可供直接跨区室比较，排卵前仅存于颗粒细胞区室。确立了双轨分析中共享阶段比对与颗粒细胞独有终末延伸并行的分析框架。
3.2. Raw-read quality control and trimming outcomes / 3.3. Preprocessing improved technical read quality…
原始FASTQ普遍存在3′接头污染及末端质量衰减，fastp修剪后Per Base Sequence Quality通过数由116升至244，Adapter Content通过由1升至119，Q30碱基占比明显上升。残留GC含量及序列组成警告属转录组文库生物学特征，修剪后数据被锁定为Salmon定量的输入。
3.4. Salmon quantification completed successfully across all libraries
所有148个文库均成功完成GRCh38脱靶感知Salmon转录本丰度估计，产生完整quant.sf文件供下游基因水平汇总。
3.5. Expression filtering reduced the imported gene set to an analysis-ready matrix
tximport初始81,412基因→去零计数63,607→edgeR表达感知过滤后保留24,684基因，TMM标准化生成log-CPM矩阵用于后续全部分析，样本标识全程锁定可追溯。
3.6. Oocyte-only PCA revealed progressive stage-related restructuring…
全队列PCA显示PC1(28.25%)按区室分离；卵母细胞单独PCA(PC1=18.01%, PC2=9.60%)示原始与初级重叠偏负PC1侧，次级向正PC1偏移，窦卵泡最正且沿PC2分散，支持卵母细胞呈渐进而非离散的阶段连续体。
3.7. Granulosa-only PCA showed stronger stage resolution…
颗粒细胞单独PCA(PC1=15.18%, PC2=10.11%)示早期间有重叠但窦及排卵前形成紧密独立簇，尤其排卵前最鲜明，说明去除区室效应后颗粒细胞晚期阶段转录分辨率更高。探索性筛查标记5个潜在离群颗粒细胞样本（多为原始阶段）但因符合早阶段异质性生物学预期予以保留。
3.8. Stage-wise differential expression was extensive in both oocytes and granulosa cells
卵母细胞各相邻阶段对比均大量显著（Primary vs Primordial 10,551; Secondary vs Primary 10,553; Antral vs Secondary 9,498基因，FDR<0.05），Antral vs Primordial最强(15,589)，方向不对称——Primary vs Primordial以下调为主，Secondary vs Primary以上调为主，最大|log2FC|=12.90，体现渐进转录重编程。颗粒细胞早期对比弱（Secondary vs Primordial仅578基因显著），但窦及排卵前相关对比极强烈（Preovulatory vs Primordial 17,488基因显著），以晚期大幅下调基因为主，最大|log2FC|=11.90，与PCA晚期分辨率增强一致。Omnibus阶段效应检验列出生阶段强关联基因（卵母细胞含UCHL1、ZFAND2A、BOD1、NLRP9等；颗粒细胞含ACTB、RPLP0、EEF1A1、TPT1、RBP1等）。
3.9. HSP90AA1 showed distinct stage-dependent expression dynamics in oocytes and granulosa cells
HSP90AA1在卵母细胞中随发育逐步衰减（logCPM：原始11.53→初级11.50→次级11.25→窦10.59，Antral vs Primordial log2FC=-0.94, adj.P=4.39×10-3）；在颗粒细胞中呈二次卵泡阶段峰值（原始11.46、初级11.37、次级14.18、窦11.23、排卵前9.34），Secondary vs Primary log2FC=+2.81(adj.P=1.17×10-4)及排卵前显著下调，说明其虽为两区室共有阶段关联基因，但颗粒细胞中发育重塑幅度与转折远更显著。
3.10. Functional enrichment resolves distinct developmental programs in oocytes and granulosa cells
卵母细胞Antral vs Primordial上调基因富集DNA重组(recombination)、染色体分离(segregation)、DNA复制与双链断裂修复等基因组调控程序；Primary vs Secondary中初级卵母细胞高表达基因富集RNA剪接(splicing)、mRNA加工(processing)、核糖核蛋白复合体生物合成及翻译起始，次级高表达侧无凝聚性富集信号——提示卵母细胞成熟含晚期基因组维护/染色质调控激活及中间阶段RNA加工架构重组。颗粒细胞Antral vs Secondary上调基因富集线粒体翻译、氧化磷酸化(OXPHOS)、电子传递偶联ATP合成，该代谢程序在Preovulatory vs Antral进一步强化，而Antral相对较高基因富集DNA复制/染色体分离——说明进入排卵前状态时颗粒细胞从较增殖/染色体动力学程序转向代谢特化的晚熟程序，整体颗粒细胞成熟以渐进线粒体与生物能量重编程为核心。
3.11. Co-expression programs underlying follicular progression in the oocyte and granulosa compartments
卵母细胞WGCNA鉴定到最大turquoise模块(r=+0.909, FDR=5.57×10-28)为正关联晚期成熟程序（hub含GEMIN5、BOD1、BTG4等，契合基因组调控重塑），royalblue(r=-0.755)代表早期翻译/生物合成程序（hub含COX6A1、RPS14、RPL19等），blue模块(r=-0.745)代表早期调控/结构程序（hub含GPM6B、RELN、SATB1、FGF13等，HSP90AA1归属此早期衰减blue模块），两反向模块随成熟进展递减。颗粒细胞WGCNA turquoise模块(r=+0.918, FDR=7.54×10-27)为正关联晚熟生物能量程序（hub含MRO、GRAMD1B、CNPY2等，与线粒体及OXPHOS富集吻合），magenta(r=-0.718)为广域早颗粒信号/应激/胞外基质样程序（hub含HSPB8、COL14A1、IER3等），orangered4(r=-0.578)为早颗粒调控/黏附相关程序（hub含DSP、CDH3、JAG1等）。两区室共有"一个主导晚期模块+一个或多个反向早期模块"的系统架构，但hub基因几无重叠，证实平行发育逻辑下分子驱动具区室特异性。
3.12. Machine-learning validation of compartment-specific developmental programs
卵母细胞阶段分类以选定3个WGCNA模块特征基因（selected eigengenes）为输入时平衡准确率最高（0.914±0.005），优于差异表达基因集(0.775)，说明卵母细胞阶段信息最有效地被紧凑模块级摘要捕获；颗粒细胞最佳为特征基因集衍生的差异表达面板（DE panel, 平衡准确率0.839±0.036），紧凑特征基因表现较差(0.695)，说明颗粒细胞阶段判别信息分布于更广转录本水平特征空间、难压缩至少数特征基因。混淆矩阵显示卵母细胞残差误分集中于相邻早期间期，颗粒细胞晚期阶段（窦、排卵前）分辨强而早/中期重叠大，与PCA及DE模式呼应。
3.13. HSP90AA1 links an early oocyte WGCNA program to broader integrative evidence
HSP90AA1归入卵母细胞blue模块（早期状态反向模块），模块成员度(module membership)及早衰阶段基因显著性（gene significance for stage progression为负）与DE中渐进衰减一致，虽非top hub但通过模块结构间接联系ML优选特征基因空间，体现基因水平行为嵌于 broader co-expression program。
讨论与结论总结
讨论部分指出：本研究支持人卵泡发生呈层级结构——区室身份为转录组首要变异轴，卵母细胞呈相对连续渐进重塑（相邻阶段即有大幅DE，PCA重叠多），颗粒细胞早期变化温和而窦及排卵前出现强烈晚期重组（侧重代谢/线粒体激活及脱离增殖染色质程序）。HSP90AA1具区室特异动态，佐证协同但非等同的发育逻辑。WGCNA揭示两区室共有"晚期模块兴起+早期模块衰减"架构，但hub基因基本无交叠，说明系统层面平行而在分子执行上层层有别。机器学习进一步量化此不对称性。研究局限含单队列推断、未证因果及部分模块可能含复合生物学，但提供了区室感知的卵泡成熟转录程序视图。
结论（翻译）：本对人卵泡发生的整合双轨再分析表明，卵母细胞与颗粒细胞在整体协调框架下通过明显不同的转录组及调控轨迹推进发育。区室身份是数据集变异的主导轴，证实卵泡成熟不可解读为单一均一过程。在此结构中，卵母细胞表现更广泛且连续的阶段关联转录重塑，而颗粒细胞早期阶段变化较小随后出现更强晚期重组。功能富集、共表达架构及预测建模均指向同一核心模式：两区室虽经历有序发育进程，却启动不同生物学程序且阶段信息在其转录组中编码强度不同。卵母细胞成熟更多关联调控与发育重塑，颗粒细胞成熟更紧密联结代谢与线粒体激活。HSP90AA1的行为进一步阐释此种区室特异性分化，支持其作为本研究所识别广泛发育框架中有意义的标志物。综上，这些发现通过显示发育协调与调控分化共存于卵泡各阶段，精炼了人卵泡发生的转录组图谱。此区室感知视角为未来卵巢生物学机制研究及女性生殖功能相关生物标志物发现和转化研究提供更佳基础。
