Germline genomic and methylomic dynamics following three generations of early-life metabolic challenges

技术摘要：三代早期生命代谢挑战后的生殖细胞基因组与甲基化组动态研究

问题陈述
环境和饮食因素，特别是早期生命阶段的代谢挑战（如营养过剩），已知会对健康和发育产生多代效应。虽然先前的研究已经证实，代谢紊乱可以通过表观遗传机制（如 miRNA、DNA 甲基化）通过父系生殖细胞进行传递，但关于多代（连续性）暴露对基因组产生的长期影响仍不为人所知。具体而言，目前尚不清楚持续三代的早期生命代谢压力如何影响生殖细胞基因组（SNP 和 CNV），以及这些基因组变化是否由表观基因组改变（DNA 甲基化）驱动或与之相关。大多数现有研究侧重于单次发育暴露的跨代传递（跨代效应），而非连续的多代暴露，而后者更符合环境致肥胖素（obesogens）的真实情况。

方法论
本研究采用 ICR (CD1) 杂合系小鼠模型，通过减少窝仔数来模拟早期生命代谢挑战。

• 实验设计： 追踪了三代（F0–F2）。过营养（ON）组的窝仔数被减少至每窝 4 只，而对照（CT）组为每窝 8 只。该模型已被证实可诱导新生儿肥胖及随后的代谢综合征。每一代的雄性个体均与外部健康的雌性交配，以隔离父系谱系效应。
• 表型分析： 在 4 月龄时评估生理和代谢参数（体重、血糖、胰岛素、肝脏和脂肪组织重量）。
• 测序策略： 提取精子 DNA，并采用并行分析方法：利用基因型测序（GBS）进行基因组变异分析，以及结合甲基化免疫沉淀的 GBS-MeDIP 进行甲基化谱分析。
• 基因组分析：
  • SNPs： 遵循 GATK 最佳实践进行识别。使用主成分分析（PCA）和层次聚类（身份状态，Identity by State）来评估遗传相似性和处理效应。
  • CNVs： 使用 GATK 的 GermlineCNVCaller 进行生殖细胞拷贝数变异（CNV）检测。应用了严格的过滤流程，包括随机抽样对照以消除假阳性。
  • 注释： 使用置换检验确定富集程度，从而对 CNV 进行转座元件（TEs）和重复序列（REs）的注释。
• 表观基因组分析：
  • DMRs： 通过代内（CT vs. ON）和代间（F1 vs. F2）比较，利用 Mann-Whitney 检验识别差异甲基化区域（DMRs）。
  • 通路分析： 对与甲基化窗口相关的基因进行 KEGG 和 PathBIX 通路富集分析。
• 整合： 计算 Spearman 秩相关，以调查父系甲基化水平与后代中新发 SNP 或 CNV 出现之间的关系。

关键结果

1. 基因组影响 (SNPs)： SNP 数据的 PCA 显示，不相关的 ON 家族聚集在一起，并与严格按亲缘关系聚类的 CT 家族区分开来。这表明代谢挑战诱导了一种共同的基因组特征。处理效应解释了约 4.6% 的遗传方差。此外，与 CT 组相比，ON 组在后代中表现出新发 SNP 的减少，表明存在遗传多样性的约束。
2. 基因组不稳定性 (CNVs)： ON 组出现了特定的 CNV 事件（有 7 个推定的事件从 F1 维持至 F2）。这些事件在长散布核元件（LINEs）和长末端重复序列（LTRs）中显著富集。缺失事件在世代间保持稳定，而重复事件则呈现扩张趋势。
3. 表观基因组影响： 甲基化组的全局 PCA 显示，ON 和 CT 组之间没有明显的分离，表明不存在广泛的全局甲基化偏移。然而，通路富集分析显示，与 CT 组相比，ON 组在多糖代谢、脂肪分解、P53、mTOR 和 hedgehog 信号传导相关的通路上表现出较弱或缺失的富集信号。
4. 区域甲基化变化： DMR 分析识别了 F2 代中的特定变化，包括 LINEs 和 SINEs 的低甲基化，以及 Fhod3 和 Etl4 基因的高甲基化。代间比较也强调了 Fhod3 和 $\gamma$-卫星重复序列中的 DMRs。
5. 基因组-表观基因组动态：
   • SNPs： 在 CT 组中，较高的父系甲基化与后代中较高的 SNP 出现率相关（这与甲基化胞嘧啶的可变性已知特性一致）。相反，在 ON 组中观察到了反向相关，这表明在代谢压力下，甲基化的 CpG 可能更不易发生突变，这可能是由于针对氧化损伤的保护机制所致。
   • CNVs： CNV 中 LINEs/LTRs 的富集，结合这些元件的差异甲基化，表明重复序列甲基化的失调可能驱动了 CNV 的形成。

意义与主张
作者声称，本研究提供了多代代谢挑战可以诱导生殖细胞基因组和表观基因组改变的证据。主要主张包括：

• 约束 vs. 不稳定性： 代谢挑战似乎在 SNP 水平上约束了遗传多样性（较少的新发突变），同时在结构水平上诱发了基因组不稳定性（CNVs）。
• 机制： 观察到的基因组不稳定性（CNVs）很可能是由转座元件（特别是 LINEs 和 LTRs）的失调介导的，而非由广泛的全局 DNA 甲基化变化驱动。
• 现实性： 通过追踪连续暴露的三代，本研究提供了比单次暴露跨代模型更真实的视角，展示了环境致肥胖素如何影响进化潜力和基因组架构。
• 新颖性： 这是首个专门探讨在多代暴露情境下，基因组如何受精子 DNA 甲基化变化影响的受控研究，强调了表观遗传扰动可能导致结构性基因组变异的复杂相互作用。

论文得出结论，早期生命阶段的代谢挑战具有长期的基因组架构影响，这种影响可能由转座元件活性介导，并强调了在现实暴露场景下同时研究基因组和表观基因组动态的必要性。