基于遗传锚点的哮喘气道上皮DNA甲基化调控基因表达的因果效应研究

《SCIENCE ADVANCES》：Estimating in silico causal effects of DNA methylation on gene expression through genetic anchors in airway epithelium in asthma

本研究针对DNA甲基化（DNAm）与基因表达关联（eQTM）中因果方向不明的问题，利用SNP作为遗传工具变量，在哮喘患儿鼻上皮中开展因果中介分析。结果显示73%的eQTM对中DNAm是基因表达的“因”，且显著富集于免疫通路，提示DNAm是哮喘防治的潜在靶点。

　　

背景：表观遗传学中的“鸡生蛋”难题

在哮喘等复杂疾病的研究中，科学家们经常观察到DNA甲基化（DNA methylation, DNAm）与基因表达之间存在强烈的统计关联（即eQTM, expression quantitative trait methylation）。然而，一个根本性的问题长期悬而未决：是DNA甲基化的变化导致了基因表达的改变，还是基因表达反过来影响了甲基化水平？ 厘清这一因果关系，对于将表观遗传标记从“生物标志物”升级为“可干预靶点”至关重要。

传统的孟德尔随机化（MR）等方法在大样本中效果显著，但在中小规模的人群队列中往往因统计效能不足而“失灵”。此外，既往研究多局限于顺式（cis）调控或血液样本，对于哮喘发病关键场所——气道上皮——中的跨染色体（trans）调控网络，其因果结构仍是一片迷雾。

技术路线概览

本研究利用EVA-PR（Epigenetic Variation and Childhood Asthma in Puerto Ricans）队列的455例（含219例哮喘患儿）鼻上皮样本，构建了SNP-DNAm-基因表达的三元因果链。研究采用了JS-AB（Joint Significance - Adaptive Bootstrap）因果中介分析模型，在控制假阳性的同时大幅降低了对样本量的要求。通过整合cis/trans-eQTM、meQTL（甲基化数量性状位点）和eQTL（表达数量性状位点）分析，系统绘制了气道上皮中表观遗传调控的因果图谱。

研究结果解析

Far more SNPs influence DNAm than gene expression（遗传变异更倾向于调控甲基化）

在对约900万对eQTM关联的筛查中，研究人员发现了一个关键现象：SNP对DNA甲基化的影响（meQTL）远多于其对基因表达的直接作用（eQTL）。在约338万对与SNP相关的eQTM中，61%仅表现为meQTL，而仅12%同时具有meQTL和eQTL信号。这表明，在气道上皮中，遗传变异更倾向于通过“先改变甲基化，再间接影响表达”的级联路径发挥作用，而非直接调控转录。

Genetic coregulations of DNAm and gene expression are uncommon（遗传共定位并不常见）

进一步分析发现，meQTL和eQTL信号共享同一因果遗传变异的比例很低（仅约14%）。这意味着DNA甲基化和基因表达在很大程度上受不同的遗传机制调控，而非由同一个SNP“双管齐下”。这一发现强调了在分析中区分这两种遗传效应的必要性。

DNAm regulates gene expression for most eQTM pairs（甲基化是基因表达的“因”）

这是本研究最核心的结论。在通过严格检验的34,891对显著eQTM中，73%的案例支持DNA甲基化是基因表达的“因”（SNP → DNAm → Expression），而反向路径（基因表达影响甲基化）仅占28%。例如，位于CST1（Cystatin SN）基因附近的SNP rs134559，其A等位基因通过降低cg13586696位点的甲基化水平，进而导致CST1表达上调，最终与哮喘表型关联。这一具体案例验证了“遗传→表观遗传→转录”的完整因果链。

DNAm remains predominantly causal after cell-type adjustment（细胞组分调整后的稳健性）

为了排除鼻腔样本中细胞类型比例混杂的影响，研究团队引入了去卷积算法估计上皮细胞和免疫细胞比例并进行校正。结果显示，即便在调整细胞组分后，DNAm作为因果主导的比例依然高达81%。这证明了“甲基化驱动表达”这一现象是生物学上的真实效应，而非细胞异质性造成的统计假象。

Plausible causal meQTL trios significantly affect asthma-related diseases（因果链富集于免疫通路与哮喘）

功能富集分析显示，由DNA甲基化介导的因果基因（meQTL trios）显著富集于抗原处理与呈递、自身免疫疾病等免疫相关通路，且与哮喘GWAS（全基因组关联分析）信号高度重合。相比之下，直接调控表达的eQTL基因则缺乏这种特异性。这表明，表观遗传介导的调控网络在哮喘的免疫失调中扮演了枢纽角色。

结论与展望

这项发表于《Science Advances》的研究，通过创新的统计遗传学方法，有力证明了在哮喘患儿的气道上皮中，DNA甲基化主要是基因表达的“上游调控者”而非“下游结果”。这一结论具有重要的转化医学价值：

1. 1.
   靶点验证：它确认了DNA甲基化（如cg13586696）作为潜在治疗靶点的合理性，为开发针对特定甲基化位点的表观遗传编辑疗法提供了理论依据。
2. 2.
   机制洞察：揭示了哮喘遗传风险变异（如17q21区域的SNP）很可能主要是通过“挟持”甲基化机制来发挥作用的，这为理解“遗传-环境”互作提供了新视角。
3. 3.
   方法学突破：JS-AB模型的应用为小样本队列中进行高维因果推断树立了范例，尤其适用于珍贵的临床组织样本研究。

未来，基于这些因果CpG位点开发的血浆或鼻腔拭子DNA甲基化标志物，有望成为预测哮喘风险或评估治疗响应的新型“表观遗传生物标志物”。