样本收集与数据生成
获取研究所需的基因组和甲基化数据，以分析混合群体的遗传和表观遗传变异。
收集256个安博塞利狒狒个体的全基因组重测序数据和RRBS甲基化数据，并额外使用178个个体计算等位基因频率。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
混合影响野生灵长类杂交区DNA甲基化的遗传结构
Admixture influences the genetic architecture of DNA methylation in a wild primate hybrid zone
实验完整度
高该研究包含多个独立证据层级，包括自然杂交群体的基因组和表观基因组分析、种间差异比较、meQTL定位、整合mSTARR-seq功能实验，且包含功能验证和机制验证。
主要模型
肯尼亚安博塞利生态系统野生狒狒杂交群体（黄狒狒×狒狒） 未混合的黄狒狒（来自米库米国家公园） 未混合的狒狒（来自马赛马拉国家保护区及华盛顿国家灵长类研究中心） K562细胞系（用于mSTARR-seq报告分析）
重点核对
样本量：256个安博塞利狒狒个体（具有重测序和甲基化数据），另外有178个个体用于等位基因频率分析 DNA甲基化数据：RRBS，分析了636,955个常染色体CpG位点 局部祖先状态：基于全基因组重测序和超过1000万个祖先信息标记，准确率96% meQTL定位：使用IMAGE方法，分析172,044个SNP-CpG对，FDR<10% mSTARR-seq实验：使用K562细胞，6个重复，检测甲基化依赖的增强子活性
摘要
背景：杂交区在进化生物学中扮演核心角色，因为它们是研究性状和物种分化的自然实验室。基因调控的变化对此过程有重要贡献。然而，混合对杂交群体中基因调控变异的影响程度仍知之甚少。在此，我们结合来自295只杂交狒狒（一个单一、深入研究自然群体的成员）的全基因组重测序和DNA甲基化数据，研究混合如何影响这种重要表观遗传标记的遗传结构。结果：我们发现局部遗传祖先常能预测DNA甲基化水平，并重演亲本物种间的差异。通过进行甲基化数量性状位点定位，我们表明这些差异主要源于等位基因频率的进化差异。因此，杂交群体中的混合增加了DNA甲基化的变异，包括引入原本不变的、影响DNA甲基化的遗传变异。最后，我们整合大规模并行报告分析数据，表明混合来源的DNA甲基化变异改变了增强子活性和基因表达。结论：这些结果共同表明，混合如何在自然杂交区显著改变基因调控性状的遗传结构。它们还提示DNA甲基化的遗传结构在近缘灵长类中是保守的，表明遗传对基因调控的影响可能在数百万年的时间尺度上保持稳定。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 混合如何影响自然杂交群体中DNA甲基化的遗传结构？
- 核心机制
- DNA甲基化差异主要由cis作用meQTL的等位基因频率差异引起，而非效应大小或因果位点的差异；混合通过引入原本不变的调节变异增加DNA甲基化的遗传方差。
- 主要证据
- 局部祖先与100,482个CpG位点的甲基化水平显著相关，且与种间差异一致；meQTL效应大小和等位基因频率差异预测了局部祖先效应（Spearman's rho=0.762, p<10^-300）；混合使meQTL的MAF平均增加2.78%，并增加了受影响的CpG位点的甲基化水平标准差（平均增加19.7%）；mSTARR-seq实验表明少数位点具有甲基化依赖的增强子活性。
- 研究意义
- 表明混合能在自然杂交区显著改变基因调控性状的遗传结构，为理解基因调控进化和杂交的进化后果提供新视角，并支持跨近缘物种遗传效应大小可能部分通用。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
局部祖先与DNA甲基化的关联分析
检验局部遗传祖先是否与DNA甲基化水平相关，并区分局部和整体祖先效应。
使用MACAU模型，对每个CpG位点建模甲基化读段数，以局部和整体祖先为预测变量，并控制年龄、性别、等级、早期逆境、亚硫酸氢盐转化率、测序批次和亲缘关系等协变量。
种间差异分析
验证安博塞利混杂群体中祖先相关的甲基化差异是否重演未混合黄狒狒和狒狒之间的差异。
重新分析来自黄狒狒（n=6）和狒狒（n=9）的RRBS数据，检测差异甲基化CpG位点，并与安博塞利祖先效应比较。
meQTL定位
鉴定影响DNA甲基化水平的遗传变异，并测试等位基因频率差异是否解释祖先效应。
使用IMAGE方法对218,047个SNP-CpG对进行meQTL定位，结合基因型和等位基因特异性甲基化数据。
功能后果分析
探究祖先相关的CpG位点是否具有功能性影响，特别是通过改变增强子活性影响基因表达。
分析祖先相关CpG位点在基因组功能元件中的分布；测试其与祖先相关表达基因的共定位；整合mSTARR-seq数据检测甲基化依赖的增强子活性。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
DNA甲基化谱分析
甲基化数量性状位点(mQTL)定位
差异甲基化分析
基因表达关联分析
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| DNA甲基化谱分析 | 简化代表性亚硫酸氢盐测序 | -- | -- |
| 全基因组重测序 | 全基因组重测序 | -- | -- |
| meQTL定位 | IMAGE | -- | -- |
| mSTARR-seq实验 | Covaris S220聚焦超声仪 | Covaris | -- |
| Msp1限制性内切酶 | -- | -- | |
| pmSTARRseq1载体 | -- | -- | |
| M.Sss1甲基转移酶 | -- | -- | |
| K562细胞系 | -- | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 样本收集 | 安博塞利狒狒个体数量：256个具有重测序和甲基化数据，另有178个用于等位基因频率 阅读提示：Methods 研究群体和样本 |
| DNA甲基化测序 | RRBS数据的处理和分析参数，包括mapping到Panubis1.0参考基因组 阅读提示：Methods 祖先对DNA甲基化的影响 |
| 祖先效应分析 | 局部和整体祖先计算，最小等位基因频率差异20%，准确率96% 阅读提示：Methods 祖先对DNA甲基化的影响 |
| meQTL定位 | 基因型调用标准，至少50%个体有基因型，MAF≥0.05，使用IMAGE工具 阅读提示：Methods 遗传对DNA甲基化的影响 |
| mSTARR-seq实验 | 处理：M.Sss1甲基化与缓冲液对照，转染K562细胞，6个重复，48小时后收获 阅读提示：Methods mSTARR-seq实验和数据 |