隐睾睾丸m6A表观转录组重编程与天华肉羊精子发生障碍的关联解析

《Animals》：Neonatal Colostrum Intake in Dairy Heifers Can Affect Lactation Performance

　　

隐睾（cryptorchidism）是哺乳动物常见先天性生殖障碍，表现为单侧或双侧睾丸未正常降入阴囊而滞留于腹腔或腹股沟管。研究人员旨在探究隐睾如何改变RNA上最丰富的内部修饰N6-甲基腺苷（m6A）的调控格局。研究以9月龄天华肉羊为对象，采集正常与隐睾睾丸组织，组织学观察显示隐睾组支持细胞（Sertoli cell）与间质细胞（Leydig cell）数量减少、生精小管直径变窄。转录组测序（RNA-seq）表明隐睾睾丸中m6A调控因子mRNA表达失调；甲基化RNA免疫沉淀测序（MeRIP-seq）共鉴定2080个差异m6A峰，其中1449个高甲基化、631个低甲基化，主要富集于编码区（CDS）、终止密码子及3′非翻译区（3′ UTR）。整合分析发现差异修饰与差异表达基因富集于凋亡（apoptosis）、AMP活化蛋白激酶（AMPK）、叉头框O（FoxO）、肿瘤坏死因子（TNF）及丝裂原活化蛋白激酶（MAPK）通路。研究人员筛选出RRM2、GTSE1、ATM、GORAB、INSR、HSPA1L、MAPK9等7个候选基因，其功能涉及凋亡、细胞周期、DNA修复、应激反应与精子发生障碍。该工作深化了对雄性不育生物学基础的认识，并为家畜生殖健康研究提供表观转录组学依据。

研究背景与问题提出

隐睾是雄性哺乳动物常见的先天性生殖发育异常，其核心病理在于睾丸滞留于体内高温环境，造成生精上皮损伤与精子发生失败。在天华肉羊等家畜中，隐睾虽发病率不高，但可显著削弱种公羊繁殖性能。既往研究多集中于解剖与内分泌层面，对隐睾所致精子发生障碍的表观转录组机制认识有限。N⁶-甲基腺苷（m⁶A，真核mRNA最常见内部修饰）由甲基转移酶（如METTL3-METTL14复合体）、去甲基化酶（如FTO、ALKBH5）与阅读蛋白（如YTHDF家族、IGF2BP家族）共同调控，参与RNA剪接、稳定性与翻译。已有证据提示m^{66A全景图谱与转录组协同变化尚未被刻画。因此，研究人员开展本研究，以明确隐睾环境下m6A区域富集与基因表达的改变及其潜在生物学意义。}

研究人员开展的工作与结论

研究人员利用正常与先天性双侧隐睾9月龄天华肉羊睾丸组织，联合组织形态学、液相色谱-串联质谱（LC-MS/MS）、MeRIP-seq与RNA-seq，构建睾丸m⁶A甲基化图谱，并整合多维数据筛选候选基因。结论显示：隐睾睾丸呈现生精细胞耗竭与细胞组成改变；全局m⁶A含量未变，但区域修饰峰显著重排；m⁶A调控因子表达失调；差异修饰基因富集于凋亡与应激相关通路，7个候选基因可作为后续机制研究切入点。该研究为隐睾相关雄性不育提供了表观转录组学观测基础。

主要关键技术方法

研究样本来自甘肃民勤德福农业有限公司饲喂的同龄同管理9月龄天华肉羊，包含3例正常与3例先天性双侧隐睾睾丸组织。核心方法包括：苏木精-伊红（HE）组织形态计量；LC-MS/MS定量总RNA中m⁶A核苷；MeRIP-seq检测m⁶A峰分布与差异甲基化；RNA-seq定量转录本表达；exomePeak2进行峰识别与差异分析；DESeq2分析差异表达；ClusterProfiler完成GO与KEGG功能富集；qRT-PCR验证m⁶A调控因子表达趋势。

研究结果

3.1 组织学观察

HE染色与计量分析表明，隐睾组睾丸重量降低，生精小管管腔精子缺失，支持细胞与间质细胞计数及小管直径均显著小于正常组，提示生精停滞与细胞微环境退化。

3.2 全局m⁶A含量

LC-MS/MS显示两组总RNA中m⁶A摩尔比无显著差异，说明隐睾并未改变整体修饰总量，而可能影响局部修饰布局。

3.3 测序质控与比对

MeRIP-seq与RNA-seq数据经过滤后比对至绵羊参考基因组（Ovis_aries Ensembl v109），Q20、Q30与GC分布良好，满足后续峰值与表达分析要求。

3.4 甲基化峰检测与注释

共检出隐睾组16308个、正常组18826个m⁶A峰，主要位于CDS、终止密码子邻近区与3′ UTR；经典RRACH基序显著富集，证实图谱可靠。隐睾组峰强度分布与正常组总体相近但区域偏好改变。

3.5 差异甲基化基因与功能富集

差异峰分析获2080个峰（1449高甲基化、631低甲基化）。差异甲基化基因（DMGs）GO富集于转录调控、磷酸化等过程；KEGG涉及胰岛素抵抗、AMPK、Hippo与FoxO通路，提示代谢与应激相关转录后调控紊乱。

3.6 转录组概览

RNA-seq鉴定出5201个上调与5471个下调基因。差异表达基因（DEGs）富集于p53信号通路，与高温诱导生精细胞凋亡的已知机制相呼应。

3.7 MeRIP与RNA整合分析

1297个基因同时具备差异m⁶A峰与差异表达，分入高甲基上调、高甲基下调、低甲基上调、低甲基下调四类。通路富集于凋亡、AMPK、FoxO、TNF与MAPK，指向DNA损伤、代谢应激与炎症响应网络。

3.8 m⁶A调控因子表达

转录组与qRT-PCR一致显示，隐睾组METTL3与YTHDF3上调，ALKBH5、IGF2BP1、IGF2BP3下调，反映writer-reader-eraser系统失衡可能参与生精损伤。

讨论部分总结

讨论指出，隐睾睾丸的高温、缺氧与内分泌扰动环境引发生精细胞耗竭，而m⁶A区域重排并非全局总量变化，更可能通过阅读蛋白特异性影响靶RNA稳定性与翻译。p53通路相关基因（RRM2、GTSE1、ATM、GORAB）的差异甲基化提示m⁶A可能参与DNA损伤应答与凋亡阈值设定；INSR、HSPA1L、MAPK9则关联代谢应激与氧化损伤。由于采用批量组织测序，细胞组成改变对信号有贡献，因果机制尚需单细胞或条件敲除验证。

结论翻译

本研究比较了正常与隐睾绵羊睾丸的区域性m⁶A富集与转录丰度，分析揭示了二者协同差异图谱并鉴定出候选基因。由于测量基于批量组织，细胞组成变化可能贡献于观测模式。这些发现为未来研究提供了观察性资源，但m⁶A调控在生精损伤中的因果作用仍有待确立。