重新审视DNA甲基化微阵列上的遗传伪影揭示新的生物学意义

Revisiting genetic artifacts on DNA methylation microarrays exposes novel biological implications

作者信息Benjamin Planterose Jiménez, Manfred Kayser, Athina Vidaki
PMID34548083
发布时间2021-09
DOI10.1186/s13059-021-02484-y

实验完整度

高

研究涉及多个独立证据层级,包括基于MZ双胞胎的遗传控制数据集、多种细胞组织样本、WGBS验证、多种工具验证及机制模型,且包含功能验证和机制验证。

主要模型

E-risk双胞胎队列(450K,全血) MZ-Adipose(MuTHER,皮下脂肪组织) Twins WGBS(全血和脂肪组织)

重点核对

遗传控制:使用MZ双胞胎配对一致性进行验证 SNP MAF估计:与1000 Genomes EUR数据比对 共甲基化窗口:±1kb WGBS验证:MZ/DZ双胞胎配对一致性及组织差异

摘要

背景:Illumina DNA甲基化微阵列实现了表观基因组水平的分析,广泛应用于发现健康和疾病中新的DNA甲基化变异。然而,微阵列的探针设计无法充分考虑人类巨大的遗传多样性,导致遗传伪影的产生。在DNA甲基化遗传力和甲基化数量性状位点的研究中,区分真实的遗传影响与伪影尤为重要。尽管这很重要,但由于对这些伪影的机制理解有限,当前应对遗传伪影的策略滞后。结果:为解决这一问题,我们开发并评估了UMtools,一个R包,包含基于荧光强度信号量化并鉴定遗传伪影的新方法。利用我们的方法,我们在一个遗传控制的单卵双胞胎数据集上建模并验证了已知的SNP/indel,并从DNA甲基化数据中估计了次要等位基因频率,经验性检测了dbSNP中未包含的变异。此外,我们识别了遗传伪影彼此之间或与印记、X染色体失活或组织特异性调控相互作用的实例。最后,我们提出了一种基于共甲基化的新策略,可以区分遗传伪影和真实的基因组影响。结论:我们提供了一个图谱,以导航DNA甲基化微阵列上遇到的遗传伪影的巨大多样性。总体而言,我们的研究为DNA甲基化微阵列中表观遗传变异的遗传组分研究范式转变奠定了基础。补充信息:在线版本包含补充材料,可访问10.1186/s13059-021-02484-y。关键词:DNA甲基化微阵列,遗传伪影,单卵双胞胎,meQTL

实验结论

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

研究问题
DNA甲基化微阵列中遗传伪影如何运作,以及如何将其与真实的遗传影响区分开来。
核心机制
遗传变异(如SNP/indel)位于探针的CpG/SBE位点或结合位点,可导致探针失败或伪装成甲基化等位基因,从而产生伪影。此外,共甲基化模式可用于区分meQTL与遗传伪影。
主要证据
基于450K芯片的E-risk双胞胎队列、多种组织样本及WGBS数据,通过UMtools工具量化荧光信号,识别伪影并与共甲基化关联。
研究意义
挑战当前基于探针排除列表的做法,提出数据驱动的策略,为DNA甲基化遗传组分研究提供新方法。

研究路径

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

1

开发并基准化UMtools工具

建立分析DNA甲基化微阵列原始荧光信号的工具,用于检测遗传伪影。

开发了U/M plots、bGMM、CVlogT、BC(CVlogT)、corMZ(CVlogT)和K-caller等工具,并在E-risk双胞胎队列上评估其性能。

2

注释450K探针相关遗传变异

根据dbSNP151,系统注释450K和EPIC平台上与探针相关的SNP和indel。

将遗传变异分为CpG位点、SBE位点(I型探针)及其他结合位点,并分类CpG/SBE-SNP为16种类型。

3

评估CpG/SBE-SNPs的伪影表现

系统评估CpG/SBE位点SNP如何导致遗传伪影。

使用UMtools的工具分析CpG/SBE-SNP类别,发现大部分符合预期,但某些类别(如type I (+) SNP1: C↔T)出现意外探针失败,并提出了包含内部CpG的模型解释。

4

评估探针结合位点SNP的影响

评估远离CpG/SBE位点的SNP如何影响探针性能。

分析探针结合位点SNP与3'端距离对探针失败的影响,发现SNP距离3'端15bp以上无影响。

5

检测indels和未注释变异

评估indels和未注册变异对荧光信号的影响。

发现indel可导致探针失败或伪装等位基因,并检测到未注释的大indel影响六个探针。

6

分析高阶遗传变异及与生物变异交互

识别复杂遗传伪影,如三等位SNP、双SNP、与印记、X失活、组织特异性等交互。

发现三等位SNP、CpG位点双SNP形成四簇、遗传伪影与印记/X失活/组织特异性交互等案例。

7

评估探针排除列表的局限性

评估现有探针排除列表的假阳性和假阴性。

分析van Dongen研究的排除列表,发现其高度富集遗传伪影探针,但仍有少量伪影泄漏,且排除列表包含大量健康探针。

8

验证NINJ2 meQTL

通过共甲基化区分meQTL和遗传伪影,并验证NINJ2基因内meQTL。

使用共甲基化分析、MAF估计、SNP imputation、WGBS验证等方法,确定rs34038797为可能的顺式因果变异,并构建机制模型。

研究方法

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

产品清单

实验环节名称品牌货号
preprocessCore----
minfi----
illuminaio----
EMCluster----
dbscan----
GEOquery----
Sushi----
ggplot2----
plotly----
bedtools----
Impute2----

关键环节

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

环节核对要点
数据来源
所使用的队列及组织类型
阅读提示:Methods部分:Datasets
遗传控制
MZ双胞胎配对一致性
阅读提示:Methods部分:UMtools benchmarking
共甲基化分析
共甲基化窗口大小
阅读提示:Methods部分:Verification of the NINJ2 meQTL
统计分析
MAF阈值
阅读提示:Methods部分:Statistical analysis of genetic artifacts