生育力的表观遗传学见解:奶牛繁殖中免疫细胞甲基化的参与

Epigenetic insights into fertility: involvement of immune cell methylation in dairy cows reproduction†

作者信息Lotfi Bouzeraa, Helene Martin, Pascal Dufour, Jessica C S Marques, Ronaldo Cerri, Marc-Andre Sirard
PMID39908274
发布时间2026-08-18
DOI10.1093/biolre/ioaf020

实验完整度

基于12头奶牛血液样本的全基因组甲基化测序和生物信息学分析,未涉及功能验证或机制实验。

主要模型

荷斯坦奶牛血液样本(可育组与生育力低下组,每组6头)

重点核对

模型:荷斯坦奶牛血液样本,分组为可育(n=6)和生育力低下(n=6) 样本来源:同一农场、相似年龄(3-8岁)的奶牛血液 分组标准:根据一年后繁殖相关淘汰记录分为生育力低下组,无生育问题的仍在群奶牛为可育组 检测方法:全基因组酶促甲基化测序(EM-seq),Illumina NovaSeq平台 数据分析:差异甲基化分析(q<0.01,|甲基化差异|>25%)

摘要

不孕和产后生殖疾病是养牛业面临的重大挑战,母体免疫系统识别和耐受胚胎的能力对于成功妊娠至关重要。造血细胞中的DNA甲基化可能影响产后生育问题的易感性,因此识别表观遗传变化对于可持续动物生产至关重要。本研究旨在表征免疫细胞中与生育力相关的甲基化组,可能有助于早期检测生育力低下。利用全表观基因组测序和酶促甲基化测序,我们分析了疾病发生前来自十二头荷斯坦奶牛的血液DNA甲基化模式。我们的研究结果显示,可育和生育力低下奶牛之间有216,990个差异甲基化胞嘧啶(DMC)。值得注意的是,三个基因——干扰素tau-3(IFNT3)、KIAA0825和RAS相关蛋白2A——在可育和生育力低下奶牛之间的差异甲基化表现出高度显著性。IFNT3对早期胚胎发育至关重要,在生育力低下奶牛的转录起始位点(TSS) shores区域有7个DMC。此外,KLRA1基因(Ly49)在分析的所有五个基因组区域(TSS shores、外显子、内含子、下游和远端基因间)均含有DMC。其广泛的差异甲基化突出了其对生育力的潜在影响。关键白细胞介素基因,包括IL6、IL15、IL22和IL36G,也显示出多个DMC,进一步证实了免疫系统在牛生育力中的作用。这些发现说明了免疫细胞表观遗传学对牛产后生育力的潜在控制。此外,这项研究表明,仅用220个生物标志物就可能估计发生生育力低下的风险,为加强动物健康管理和改善生育治疗铺平了道路。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
免疫细胞的DNA甲基化模式是否与奶牛产后生育力相关,并可作为生育力低下风险的预测生物标志物?
核心机制
生育力低下奶牛中IFNT3、RAP2A等基因的TSS shores区域高甲基化可能抑制其表达,而KIAA0825等基因低甲基化可能促进表达,从而影响妊娠建立和维持。
主要证据
基于12头奶牛(6头可育,6头生育力低下)的全基因组甲基化测序,鉴定出216,990个DMC,其中IFNT3、KIAA0825和RAP2A等基因在TSS shores区域具有显著差异甲基化。
研究意义
本研究首次探索成年奶牛血细胞DNA甲基化对生育问题的易感性,并提出了220个可预测生育力低下风险的甲基化生物标志物,有望用于早期识别高风险奶牛并改善繁殖管理。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

样本选择和分组

根据奶牛繁殖结果(淘汰原因)将血液样本分为可育和生育力低下两组。

从同一农场采集12头奶牛血液,一年后根据繁殖相关淘汰记录选择6头为生育力低下组,其余6头为可育组。

2

全基因组甲基化测序和数据处理

获取奶牛血液样本的全基因组DNA甲基化图谱。

使用酶促甲基化测序(EM-seq)构建文库并进行全基因组测序,经生物信息学流程进行甲基化调用。

3

差异甲基化分析

比较可育与生育力低下奶牛间的差异甲基化位点。

使用MethylKit进行单CpG和区域差异化分析,筛选显著DMCs并进行基因组区域注释。

4

功能富集和基因网络分析

解析差异甲基化基因的生物学功能和潜在相互作用网络。

使用DAVID进行功能富集分析,并利用BioGRID数据库构建基因相互作用网络。

5

预测标志物筛选

识别能区分可育与生育力低下奶牛的最小甲基化标志物集。

基于高度显著DMCs(-log10(q-value)>10),使用sPLS-DA筛选出220个关键标志物。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
Illumina NovaSeq 测序平台Illumina--
酶促甲基化测序----
Génome Québec测序平台Génome Québec--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
样本选择
分组定义:生育力低下组为采样后一年内因繁殖相关问题被淘汰的奶牛,可育组为仍在群且无生育问题的奶牛;年龄范围3-8岁;样本量n=6每组。
阅读提示:Materials and methods: Animal selection and sample collection
测序与文库制备
使用1 ml血液;EM-seq酶处理(TET2、T4-BGT、APOBEC3A);内部对照λDNA和pUC19;测序平台Illumina NovaSeq;目标读段数约450M 150bp双端读段。
阅读提示:Materials and methods: Enzymatic methyl sequencing
甲基化调用
使用fastp进行过滤、修剪和去重;使用Bismark比对到Bos taurus ARS-UCD1.2.104参考基因组;DRAGEN Methylation Pipeline运行。
阅读提示:Materials and methods: Methyl calling
差异甲基化分析
CpG覆盖度:至少6个reads,且每组至少5个样本;区域分析窗口500bp,重叠250bp,至少5个CpG;差异统计:logistic回归,SLIM方法调整q值;协变量:年龄、父母差异、λ转换率;筛选阈值:q<0.01,甲基化差异>25%。
阅读提示:Materials and methods: Differential methylation analysis
标志物筛选
sPLS-DA分析:验证参数Mfold,folds=5,nrepeat=100,dist='max.dist',measure='BER';基于DMCs -log10(q-value)>10筛选。
阅读提示:Materials and methods: Statistical analyses and visualization methods; Results: Top predictive markers associated with fertility