整合单细胞与批量转录组分析揭示UTF1相关精子发生功能障碍并鉴定ACTRT2与AKAP4作为非梗阻性无精子症（NOA）生物标志物

《Frontiers in Genetics》：Integrated single-cell and bulk transcriptomic analyses reveal UTF1-associated spermatogenic dysfunction and identify ACTRT2 and AKAP4 as biomarkers in non-obstructive azoospermia

引言：非梗阻性无精子症（non-obstructive azoospermia, NOA）是一种由精子发生失败（spermatogenic failure）导致的严重男性不育症，UTF1对维持未分化精原细胞（undifferentiated spermatogonia）和调控生殖细胞发育至关重要。本研究旨在探讨UTF1相关精子发生功能障碍在NOA中的作用，并鉴定潜在的诊断生物标志物（diagnostic biomarkers）。

方法：研究人员将NOA转录组数据集与来自Utf1条件敲除和正常小鼠睾丸的单细胞RNA测序（scRNA-seq）数据整合。应用生物信息学分析和机器学习筛选生物标志物，同时开展功能富集、免疫浸润、细胞-细胞通讯（cell-cell communication）和伪时间轨迹分析（pseudotime trajectory analysis），并在Utf1缺陷睾丸中进行组织学验证（histological validation）。

结果：精原细胞（spermatogonia）、支持细胞（Sertoli cells）和精子细胞（spermatids）被鉴定为关键细胞群体。ACTRT2和AKAP4被确认为一致下调的NOA生物标志物，并被整合入一个校准良好的列线图（nomogram）。这两个基因与细胞骨架重塑（cytoskeletal remodeling）和减数分裂调控（meiotic regulation）相关。NOA还以免疫浸润（immune infiltration）和细胞-细胞通讯改变为特征，而Utf1缺陷损害睾丸发育和精子发生。

讨论：总之，UTF1缺陷与NOA中的精子发生功能障碍和微环境重塑（microenvironmental remodeling）密切相关，而ACTRT2和AKAP4可能作为该过程的潜在诊断生物标志物。这些发现为NOA的分子基础提供了新见解，并可能支持未来的疾病分层和靶向干预策略。

研究背景与问题：男性不育是重要公共卫生问题，约7%男性受不育影响，无精子症是最严重临床表型之一；其中非梗阻性无精子症（non-obstructive azoospermia, NOA）占无精子症可达70%，多数病因不明。睾丸精子提取联合卵胞浆内单精子注射（TESE-ICSI）可使部分患者获得生育，但对睾丸内完全缺乏成熟精子者仍无有效治疗。精子发生依赖精原细胞、支持细胞（Sertoli cells）等支持细胞、免疫微环境和细胞间通讯，NOA涉及生殖细胞分化缺陷、支持细胞功能障碍、免疫微环境改变和通讯失衡，且细胞异质性高。UTF1是未分化精原细胞/精原干细胞关键标志和调控因子，NOA中其表达失调并伴空间定位改变；但Utf1全球敲除胚胎致死，缺乏组织特异性条件敲除模型，其调控精子发生及NOA发病机制尚待阐明。为此，研究人员构建睾丸特异性Utf1条件敲除小鼠，结合小鼠单细胞转录组与人NOA批量转录组，筛选并验证关键分子生物标志物，以阐明UTF1依赖的精子发生损伤机制，为NOA精准诊断和靶向干预提供理论框架。

研究概述与意义：研究人员生成Utf1^cKO小鼠，开展从体内表型到单细胞转录组分析，并整合人NOA数据筛选验证ACTRT2和AKAP4。结论提示UTF1缺陷与NOA精子发生功能障碍及微环境重塑密切相关，ACTRT2和AKAP4可能作为潜在诊断生物标志物。该研究为NOA分子基础、疾病分层和靶向干预提供新线索，但临床适用性仍需更大、特征明确队列验证。

方法概述（不超过250字）：训练队列为GSE145467（GPL4133；10例NOA/10例对照），验证队列为GSE45887（GPL6244；16例NOA/4例对照）；小鼠单细胞数据含3只对照和3只Utf1^cKO；Utf1^fl/+小鼠购自GemPharmatech（江苏，中国），与Vasa-Cre工具鼠杂交。关键技术包括Seurat单细胞质控与注释、高维加权基因共表达网络分析（hdWGCNA）、limma差异分析、基因本体/京都基因与基因组百科全书（GO/KEGG）富集与蛋白-蛋白相互作用（PPI）/分子复合物检测（MCODE）、支持向量机-递归特征消除（SVM-RFE）/随机森林（RF）/Boruta机器学习、列线图、基因集富集分析（GSEA）与单样本GSEA（ssGSEA）免疫浸润、miRNA-mRNA-转录因子网络及药物预测、CellChat通讯、Monocle伪时间，以及苏木精-伊红（HE）、免疫荧光（IF）、免疫组织化学（IHC）、TUNEL、计算机辅助精子分析（CASA）和酶联免疫吸附试验（ELISA）验证。

研究结果：
3.1 基于单细胞测序的细胞亚群表征与hdWGCNA：对49,901个细胞和26,608个基因质控后，UMAP获得10个转录簇并注释为8类主要细胞及“Others”；选取精原细胞、支持细胞和精子细胞行hdWGCNA，软阈值7（R²=0.92），获3个共表达模块、300个关键模块相关基因（KMRGs）及240个人直系同源的精子发生相关基因。
3.2 候选基因筛选及潜在生物学机制：GSE145467中获4,688个差异表达基因（DEGs；1,550上调、3,138下调），与240个基因交集得73个候选基因；GO/KEGG显示其富集于精子细胞分化、精子活力、鞭毛组装、精子DNA浓缩、顶体反应、谷胱甘肽代谢、活性氧反应、糖酵解和脂肪酸降解等，PPI网络显示候选基因间功能互联。
3.3 生物标志物鉴定及列线图性能评价：MCODE最高分子复合物子网含13个候选基因；SVM-RFE保留ACTRT2和AKAP4，与RF、Boruta结果交集亦为二者；两队列中二者在NOA均显著下调。训练队列AUC为Boruta=1.000、RF=1.000、SVM-RFE=0.920，验证队列AUC分别为0.750、0.750、0.125，验证受限；UTF1在GSE145467中显著上调，在GSE45887中无显著差异。列线图总分149对应94.6%的NOA预测概率，HL检验P=0.911提示校准良好。
3.4 ACTRT2和AKAP4的染色体与亚细胞定位及通路富集：ACTRT2位于1号染色体，AKAP4位于X染色体；ACTRT2主要定位于细胞骨架、胞质等，AKAP4主要定位于核，并分布于胞质、骨架、线粒体等；GSEA显示ACTRT2涉及52条通路，主要为微管动力学和细胞骨架重塑，AKAP4涉及59条通路，主要为鞭毛组装和减数分裂纺锤体调控，二者共享51条通路，包括微管成核、解聚和动粒-微管连接。
3.5 免疫细胞浸润及其与生物标志物的相关性：28种免疫亚群中23种存在显著差异，22种在NOA升高，包括巨噬细胞、调节性T细胞（Tregs）、髓源性抑制细胞（MDSCs）和肥大细胞，记忆B细胞为唯一显著降低群体；Tregs与γδT细胞、CD56^dim自然杀伤（NK）细胞和MDSCs呈强正相关；ACTRT2和AKAP4与多数T细胞、NK细胞、MDSCs、巨噬细胞、肥大细胞和单核细胞显著负相关，与记忆B细胞显著正相关。
3.6 潜在调控机制与药物预测：构建miRNA-mRNA-转录因子（TF）网络，4个miRNA预测靶向ACTRT2、7个靶向AKAP4，hsa-miR-27a-3p为共享调控因子；24个TF预测调控ACTRT2、55个调控AKAP4，10个共享TF包括AR、FOXA2、SPI1和JUND；药物预测显示3种化合物靶向ACTRT2（bisphenol A、3,3′,4,4′,5-pentachlorobiphenyl、pentadecafluorooctanoic acid），2种靶向AKAP4（clidinium bromide、pivampicillin）。
3.7 关键细胞鉴定与细胞间通讯网络：ACTRT2和AKAP4在精子细胞高表达，ACTRT2在该细胞群Utf1^cKO与对照间差异显著，故精子细胞为关键细胞；Utf1^cKO组细胞间通讯数量和强度增加，涉及支持细胞、生精细胞和免疫细胞。对照中精子细胞主要经App-Cd74与巨噬细胞、经Cadm1-Cadm1与支持细胞和精原细胞作用；Utf1^cKO中新增Dll3-Notch1/3、Grn-Sort1和Cadm1-Cadm1等与周管细胞、内皮细胞和上皮细胞的相互作用。
3.8 精子细胞的功能通路与分化轨迹：精子细胞高激活代谢和离子通道相关通路，包括肠道脂质吸收、组胺受体信号、SLC22转运体有机阴离子转运和钾转运通道，而RNA加工、免疫功能和氧化应激解除相关通路受抑；精子细胞分为3个亚群，伪时间呈单一连续分化轨迹，Utf1^cKO组早期细胞较少、晚期细胞较多；ACTRT2在早期升高、中期达峰后下降，AKAP4在早中期缓慢升高至峰值后持续下降。
3.9 Utf1缺陷导致睾丸发育受损及精子数量和活力降低：IHC和Western blot显示Utf1^cKO小鼠8 dpp和35 dpp睾丸UTF1蛋白显著降低；35 dpp体重和睾丸/体重比降低，精子计数和活力显著下降，血清睾酮显著降低。
3.10 Utf1缺陷导致雄性小鼠不育：60 dpp交配实验中，Utf1^cKO雄鼠与野生型雌鼠的阴道栓数和总产仔数减少，平均窝仔数显著低于野生型配对及Utf1^flox/+Vasa^Cre+雄鼠配对。
3.11 Utf1缺陷导致第一次精子发生波功能障碍及生精小管各阶段生殖细胞丢失：HE显示10 dpp精原细胞减少，20 dpp和35 dpp精母细胞和圆精子细胞不同程度耗竭、细胞分层紊乱，管腔无成熟精子；IF显示野生型8 dpp UTF1主要位于精原细胞核，35 dpp主要位于圆精子细胞核。
3.12 Utf1缺陷降低35 dpp Utf1^cKO雄性小鼠睾丸细胞增殖能力：IHC显示PCNA和MKI67在Utf1^cKO睾丸中显著降低，提示生殖细胞增殖能力下降。
3.13 TUNEL染色显示Utf1^cKO睾丸生殖细胞凋亡增加：35 dpp Utf1^cKO生精小管TUNEL阳性细胞明显增多，阳性信号强度显著升高，提示凋亡增强参与生殖细胞进行性丢失。

讨论总结与研究结论翻译：讨论指出，NOA并非单一生殖细胞缺陷，而是精原细胞维持、支持细胞功能、局部免疫状态和细胞通讯共同失衡的异质性过程。UTF1可能维持未分化精原细胞与支持细胞间功能稳定；Utf1缺陷后早期分化受损、晚期细胞积聚且通讯增强，提示分化启动受损而下游程序异常推进并停滞。ACTRT2与精原细胞铁稳态和线粒体完整性相关，AKAP4为精子鞭毛纤维鞘关键结构蛋白；二者共享微管成核、微管动力学、动粒-微管连接、细胞骨架重塑和鞭毛组装等通路，反映从早期精原细胞失调到晚期精子结构失败的连续病理轴。Utf1^cKO小鼠睾酮降低，提示生精和内分泌功能均受影响；免疫浸润异常及Dll3-Notch、Grn-Sort1、App-Cd74等通讯增强可能代表损伤微环境中的代偿或病理重塑。局限包括样本量较小、跨物种整合偏倚、缺乏因果验证、未对UTF1本身进行PCR/qRT-PCR/Western blot验证、无救援实验、列线图外部验证有限及临床参数相关性不显著。研究结论翻译：为阐明UTF1在睾丸发育和精子发生中的作用，本研究建立了生殖细胞特异性Utf1条件敲除（Utf1^cKO）小鼠模型。结果表明，Utf1缺陷与精子发生功能障碍和微环境重塑密切相关。整合分析进一步鉴定ACTRT2和AKAP4为该过程的潜在诊断生物标志物。总体而言，这些发现揭示了UTF1在维持正常精子发生中的关键作用，并为NOA的分子基础提供了新见解，支持其在疾病分层和靶向干预中的潜在价值。然而，由于这些发现主要来自实验模型和整合分析，其临床适用性仍需在更大、特征明确的患者队列中进一步验证。