《Journal of Animal Science and Biotechnology》：Integrated methylome–transcriptome profiling reveals epigenetic regulation of immune activation pathways and CSN3-associated lactation repression in bovine subclinical mastitis
编辑推荐：
背景：亚临床乳腺炎(SCM)是奶牛生产的主要制约因素，由复杂的宿主-病原体相互作用驱动。尽管SCM相关的转录反应已被广泛研究，但稳定调控这些程序的表观遗传机制仍缺乏深入表征，特别是在杂交牛群体中。本研究旨在通过整合Vrindavani(Bos taurus ×
背景：亚临床乳腺炎(SCM)是奶牛生产的主要制约因素，由复杂的宿主-病原体相互作用驱动。尽管SCM相关的转录反应已被广泛研究，但稳定调控这些程序的表观遗传机制仍缺乏深入表征，特别是在杂交牛群体中。本研究旨在通过整合Vrindavani(Bos taurus × Bos indicus)牛乳体细胞的全基因组甲基化和转录组数据，表征基于DNA甲基化的调控网络。方法：对SCM感染和健康对照奶牛的乳体细胞进行了全基因组甲基化(n=6)和相应的转录组谱分析(n=6)。结果：差异甲基化分析(q-value < 0.05)鉴定出62,940个差异甲基化胞嘧啶(DMCs)、7,706个差异甲基化区域(DMRs)和6,203个差异甲基化基因(DMGs)，SCM中主要表现为低甲基化偏向。利用严格的甲基化(≥10%)和表达变化(|log? fold change| ≥ 1; P of GMM < 0.001)阈值进行的整合分析鉴定出1,407个差异甲基化且表达的基因(DMEGs)。功能富集分析揭示了47条KEGG通路和30个基因本体(GO)生物过程术语(FDR < 0.05)，主要与免疫信号传导和炎症反应相关。相比之下，一部分DMEGs显示出与泌乳和代谢相关基因的甲基化相关抑制。通过qPCR实验验证了选定的基因，包括炎性介质S100A8的上调和关键乳蛋白基因κ-酪蛋白(CSN3)的下调。结论：这些发现提供了乳体细胞中SCM的DNA甲基化和转录组景观的综合视图，并证明表观遗传重塑与免疫通路的协调激活以及泌乳相关功能的抑制相关。结果为理解亚临床乳腺炎的分子基础做出了贡献，并可能为未来的生物标志物开发和改善奶牛疾病抗性的表观遗传策略提供信息。
论文解读：牛亚临床乳腺炎表观遗传调控机制的整合多组学研究
研究背景
亚临床乳腺炎(SCM)是全球乳业面临的一项重大挑战，因其缺乏可见临床症状，常导致慢性感染，严重影响牛奶质量和产量。尽管现有的基因组关联研究(GWAS)和转录组学(RNA-seq)已识别出参与免疫反应的关键基因，但这些技术仅能提供基因表达动态的瞬时快照，往往无法捕捉支撑SCM慢性炎症状态的稳定调控机制。表观遗传学，特别是DNA甲基化，作为一种稳定且动态的“细胞记忆”形式，能够在不改变DNA序列的情况下响应病原体刺激调节基因表达。然而，不同乳腺炎致病病原体诱导的DNA甲基化特征及其在杂交牛群体中的调控网络仍未得到充分探索。为此，研究人员针对印度重要的杂交牛品种Vrindavani(Bos taurus × Bos indicus)开展了全基因组甲基化与转录组的整合分析，旨在阐明SCM发病过程中的表观遗传重编程机制。该研究成果发表于《Journal of Animal Science and Biotechnology》。
关键技术方法
研究人员选取了Vrindavani杂交奶牛群，通过加州乳腺炎测试(CMT)和直接显微镜体细胞计数(DMSCC)筛选并建立了SCM组(SCC ≥ 200,000 cells/mL)和健康对照组(HC, SCC ≤ 150,000 cells/mL)队列，收集其乳体细胞。采用全基因组甲基化测序(WGMS)和RNA测序(RNA-seq)技术分别获取DNA甲基化和转录组图谱。利用MethGET程序和Gaussian Mixture Model(GMM)算法对甲基化与表达数据进行整合分析，识别差异甲基化且表达的基因(DMEGs)。随后利用clusterProfiler等工具进行功能富集分析，并通过实时荧光定量PCR(qPCR)对关键候选基因进行实验验证。
研究结果
1. 亚临床乳腺炎的全基因组DNA低甲基化特征
研究人员通过对乳体细胞的全基因组甲基化测序分析发现，SCM样本表现出显著的全局DNA低甲基化现象。与健康对照组相比，SCM组中平均甲基化水平在启动子、基因体及CpG岛(CGI)结构中均显著下降。差异甲基化分析鉴定出62,940个差异甲基化胞嘧啶(DMCs)，其中93.8%为低甲基化位点，仅有6.2%为高甲基化。在区域层面，7,706个差异甲基化区域(DMRs)中98.8%表现为低甲基化。染色体分布分析显示，这种低甲基化偏向在所有染色体上普遍存在，且热点区域主要集中在第18、19和11号染色体上。
2. 差异甲基化胞嘧啶(DMCs)和区域(DMRs)的鉴定
研究人员进一步分析了DMCs和DMRs的基因组分布特征。结果显示，DMCs主要富集在CpG岛的岸边(CGI shores)和架区(CGI shelves)，这表明炎症应答可能涉及对这些关键调控区域的靶向去甲基化。功能注释显示，低甲基化基因显著富集于免疫信号传导以及神经发生等与发育重塑相关的生物学过程。相比之下，高甲基化基因集未检测到显著的功能富集通路，表明SCM的表观遗传特征以激活型调控为主。
3. 差异甲基化且表达基因(DMEGs)的识别
为了厘清DNA甲基化与转录调控之间的相互作用，研究人员应用高斯混合模型(GMM)对多组学数据进行了整合。研究共鉴定出1,407个独特的DMEGs。相关性分析显示，尽管相关系数绝对值较小，但启动子甲基化百分比与基因表达值之间存在显著的负相关关系。这表明经典的表观遗传调控范式——启动子低甲基化伴随转录激活——在SCM中依然成立。热图分析清晰地展示了健康组与SCM组在甲基化与表达模式上的分离。
4. 整合表观基因组和转录组分析揭示协调的免疫反应
对上调的DMEGs进行功能富集分析，研究人员发现47条KEGG通路和30个GO生物过程术语显著富集。这些通路主要涉及利什曼病、类风湿性关节炎、哮喘等感染性疾病以及自身免疫和炎症反应。基因本体分析强调了免疫反应、中性粒细胞趋化作用以及对脂多糖的细胞反应是核心过程。此外，MHC II类蛋白复合物、细胞外空间和分泌蛋白等细胞组分显著富集，提示了抗原呈递和细胞外通讯在宿主防御中的重要性。相反，下调的DMEGs功能富集较弱，主要涉及G蛋白偶联受体(GPCR)信号通路、嗅觉受体活性以及激素反应，特别是与乳蛋白基因和激素响应过程相关，这反映了泌乳相关通路的抑制。
5. qPCR分析验证启动子甲基化介导的基因调控
为了验证生物信息学分析的结果，研究人员选取了代表性基因进行qPCR验证。结果显示，RNA-Seq中观察到的表达变化趋势在所有八个选定基因中均得到了证实。具体而言，促炎基因S100A8、TLR4和CASP1表现出启动子低甲基化与转录上调的关联；而关键乳蛋白基因CSN3以及PRKAR2B、OXCT1等代谢相关基因则表现出启动子高甲基化与转录下调的强烈一致性。其中，CSN3的qPCR验证显示极显著差异(P = 0.011)，证实了表观遗传沉默是乳腺炎期间乳蛋白合成受抑的保守机制。
讨论与结论
本研究表明，SCM的病理特征不仅在于持续的免疫激活，还在于整个牛基因组的深刻表观遗传重塑。研究发现超过90%的差异甲基化胞嘧啶表现为启动子低甲基化，这与之前在牛乳腺炎及其他慢性炎症疾病中的报道一致，支持了低甲基化是乳腺炎特征性表观遗传标志的观点。这种靶向性的低甲基化导致了免疫相关基因(如S100A8、CCL3、CASP1)的协调激活，同时伴随着关键泌乳基因(如CSN3)的启动子高甲基化及其表达抑制。这种表观遗传层面的资源重新分配，使得乳腺上皮细胞和免疫细胞在应对感染时优先保障防御机制，而牺牲了产奶性能。值得注意的是，本研究鉴定出的部分DMEGs与外部数据集存在重叠，证明了这些甲基化-表达耦合可能是SCM发病机制中的保守分子特征。综上所述，该研究将SCM视为一种乳腺细胞特性的表观遗传重编程，确立了DNA甲基化作为连接免疫信号通路与代谢及泌乳程序之间的重要调控开关，为理解亚临床乳腺炎的分子基础及开发基于表观遗传的疾病监测工具提供了理论基础。
