关于人类睾丸生殖寿命的单细胞转录组图谱研究的学术报告
本文报告了一项由Lina Cui, Xichen Nie, Yixuan Guo, Pengcheng Ren, Yifei Guo, Xiaoyan Wang, Ran Li, James M. Hotaling, Bradley R. Cairns & Jingtao Guo等研究者共同完成的原创性研究工作。该研究团队主要来自中国科学院动物研究所（中国北京）、北京干细胞与再生医学研究所以及美国犹他大学医学院（包括霍华德·休斯医学研究所、亨茨曼癌症研究所、泌尿外科等）。这项研究成果发表于Nature Aging期刊，于2025年4月正式出版（卷5，期号，页码658-674）。
该研究领域属于生殖生物学、衰老生物学与单细胞组学的交叉领域。尽管男性生殖衰老与生育力下降、后代健康风险增加相关，但其背后的细胞类型特异性分子机制尚不完全清楚。既往研究多集中在特定年龄段（如青年与老年）的比较，或分析精子基因组、蛋白质组等，缺乏在整个男性生殖寿命（成年后至老年）跨度内，对睾丸内所有细胞类型进行系统、高分辨率解析的研究。同时，肥胖（以身体质量指数BMI衡量）作为重要的生理影响因素，其与衰老如何协同损害男性生育力，也亟待阐明。
因此，本研究旨在构建一个覆盖男性整个生殖寿命的、全面的人类睾丸单细胞转录组图谱，以系统性回答以下几个核心科学问题：1）睾丸衰老对不同细胞类型的影响是同步的还是异步的？2）哪些细胞类型对衰老最敏感，其分子特征是什么？3）衰老对男性生育力的影响是累积性的还是即时性的？4）衰老与BMI如何相互作用，共同影响睾丸功能？最终目标是为理解睾丸衰老的分子和细胞机制奠定基础，并寻找潜在诊断标志物和治疗靶点。
本研究流程复杂且系统，主要包括样本收集与处理、单细胞测序与数据生成、生物信息学与机器学习分析、以及实验验证等多个环节。
第一，样本收集与单细胞图谱构建。 研究团队从35名年龄介于21至69岁的男性遗体器官捐献者（通过DonorConnect机构获取）处获取了睾丸样本。所有捐献者均有后代或证实具备正常精子发生能力，以确保研究聚焦于生理性衰老而非疾病状态。样本经快速处理、冷冻保存或固定后，用于后续分析。研究采用10x Genomics Chromium平台对睾丸组织进行单细胞RNA测序（scRNA-seq）。经过严格的质量控制（过滤低质量细胞和基因），最终获得了来自35个个体的共计214,369个高质量细胞的转录组数据。通过无监督聚类和UMAP降维分析，成功鉴定了14个主要的睾丸细胞亚群。这些亚群包括生殖细胞系的各个阶段：未分化精原细胞、分化中精原细胞、早期初级精母细胞、晚期初级精母细胞、圆形精子细胞和伸长精子细胞；以及多种体细胞：支持细胞、睾丸管周细胞、间质细胞、平滑肌细胞、内皮细胞、巨噬细胞、淋巴细胞和B细胞。这构建了一个覆盖整个生殖寿命的健康人类睾丸单细胞转录组图谱。
第二，机器学习构建睾丸衰老时钟与细胞特异性分析。 为了量化不同细胞类型在预测生物学年龄中的重要性，研究团队采用了监督式机器学习方法。他们将供者按年龄分为20多岁至60多岁五个组，并基于每个细胞类型的基因表达谱，测试了包括极端梯度提升分类器（XGBoost）在内的八种模型。结果显示，XGBoost模型在所有细胞类型中均表现出最佳性能。分析发现，体细胞（如内皮细胞、管周细胞、巨噬细胞、间质细胞）在年龄预测中的准确性（>0.8）显著高于大部分生殖细胞，提示体细胞对衰老的反应更强烈，可能是指示睾丸衰老状态的关键。通过分析模型的特征重要性，研究者在不同体细胞类型中均发现了免疫反应相关基因的突出贡献，表明免疫激活是睾丸衰老时钟的核心特征，并可能作为预测年龄相关生育力下降的标志物。
第三，体细胞衰老反应的双波次模式解析。 通过对不同年龄段体细胞进行配对差异基因表达分析，研究者发现了衰老相关的两波显著变化。第一波变化出现在30多岁的供者中，主要体现在睾丸管周细胞。与20多岁的供者相比，30多岁供者的TPCs表现出大量差异表达基因，功能富集分析显示，与细胞外基质合成、胶原纤维组织和氧化应激反应相关的通路显著下调，而与核糖体生物发生、细胞质翻译相关的通路上调。这一发现通过PAS染色（显示基底膜）和胶原蛋白I免疫荧光染色得到实验验证：在30多岁时，睾丸生精小管的基底膜厚度显著增加，胶原蛋白I的沉积也发生明显变化。这表明TPCs在30多岁时已进入一种衰老的“启动”状态，其特征是ECM重塑和增厚。第二波变化则集中在50多岁的供者中，涉及多种体细胞类型的广泛功能改变。支持细胞的比例下降，细胞间通讯增强，其表面与营养感知相关的受体（如胰岛素受体INSR、胰岛素样生长因子1受体IGF1R）表达上调。间质细胞的差异基因则富集在类固醇生物合成与代谢、激素反应等通路，提示其睾酮合成功能发生改变。此外，巨噬细胞的免疫反应通路和平滑肌细胞的肌肉收缩通路也发生显著变化。这些结果表明，在50岁左右，多种睾丸体细胞的功能状态发生协同转变，共同构成了衰老微环境。
第四，生殖细胞衰老的累积性效应分析。 对生殖细胞的分析揭示了不同的衰老模式。研究发现，精原干细胞的数量在50多岁时显著减少，但其转录组在整个生命周期中相对稳定，表现出对衰老微环境的渐进式反应。然而，单倍体生殖细胞（圆形和伸长精子细胞）的转录组在40-50岁时发生了显著变化。差异基因表达分析显示，与精子鞭毛运动、精子细胞分化、DNA损伤修复等相关的基因表达发生动态改变。此外，研究还利用单细胞转座子元件分析（scTE）发现，某些转座子家族（如LINE1， Alu）的表达在老年供者的单倍体生殖细胞中增加，这可能影响基因组稳定性。通过SCENIC基因调控网络分析，研究者进一步识别出在老年精子细胞中活跃的关键转录因子调控子，这些调控子参与调控纤毛组装、自噬、精子运动等功能通路。这些结果共同表明，单倍体生殖细胞的功能变化具有累积性，在中年后期（40-50岁）变得尤为明显，可能直接影响精子质量和功能。
第五，体重指数与衰老的协同效应研究。 为了探究环境生理因素BMI的影响，研究者整合了供者的BMI信息进行分析。结果发现，高BMI（特别是≥30）与圆形和伸长精子细胞的频率降低相关，提示其导致亚生育状态。更重要的是，年龄与BMI存在显著的交互作用：在40岁以后，高BMI对精子发生能力的负面影响更为突出。通过构建预测模型，研究确定年龄超过45岁或BMI超过30（尤其是BMI>40的III级肥胖）是男性亚生育的重要风险阈值。对体细胞的差异基因分析进一步显示，高BMI组（BMI≥30）在30多岁和50多岁时表现出的基因表达变化模式，与低BMI组在衰老过程中发生的变化相似，表明高BMI可能“加速”了睾丸的衰老进程。
上述流程环环相扣，逻辑紧密。首先，构建高质量的单细胞图谱是基础，为后续所有分析提供了数据资源。其次，机器学习衰老时钟的分析结果（体细胞对年龄预测更准确）直接引导研究者将重点放在体细胞的衰老反应上，并发现了TPCs在30多岁的早期变化。这一发现通过差异表达分析和组织学染色（PAS和胶原蛋白I）得到双重验证，构成了“第一波衰老”的有力证据。随后，对50多岁体细胞的深入分析，揭示了支持细胞、间质细胞等功能的协同改变，构成了“第二波衰老”。与此同时，对生殖细胞的分析则从另一个维度补充了全貌：精原干细胞数量的减少和单倍体生殖细胞转录组的累积性变化，解释了衰老过程中精子产量和质量可能下降的细胞基础。最后，整合BMI数据的分析，将生理性衰老研究与重要的环境风险因素联系起来，阐明了肥胖如何与年龄协同加剧男性生育力下降，提升了研究的临床与现实意义。所有结果相互支撑，共同绘制了一幅睾丸衰老的动态、多细胞、多因素的分子图谱。
本研究的核心结论是：人类睾丸衰老是一个异步、多阶段的过程，体细胞和生殖细胞以不同的模式和节奏参与其中。具体而言，TPCs是衰老的“先锋”响应细胞，在30多岁时即发生ECM相关的显著重塑；而在50多岁时，多种体细胞（如支持细胞、间质细胞）发生功能转变，与单倍体生殖细胞的累积性转录组变化同步，共同标志着睾丸功能的显著转折点。此外，高BMI会加剧年龄相关的睾丸功能衰退，尤其是在40岁以后。
本研究的科学价值极高：1）资源价值：创建了首个跨越整个男性生殖寿命的、大规模的健康人类睾丸单细胞转录组图谱，为生殖和衰老领域提供了宝贵的公共数据资源。2）机制洞察：首次系统揭示了睾丸衰老的“双波次”细胞动力学模型，明确了关键细胞类型和分子通路（如ECM、免疫、类固醇代谢）在特定时间点的变化，深化了对睾丸衰老机制的认知。3）方法学创新：成功将机器学习模型与单细胞组学深度结合，用于识别对衰老最敏感的细胞类型和生物标志物，为衰老研究提供了新范式。4）临床与公共卫生意义：明确了年龄（>45岁）和肥胖（BMI>30，尤其是>40）作为男性亚生育的关键风险阈值，为男性生殖健康管理、生育咨询以及开发延缓睾丸衰老或改善生育力的干预策略（如靶向ECM、调控局部免疫）提供了潜在靶点和理论依据。