通过EWMR分析识别因果CpG位点
评估DNA甲基化对IPF风险的因果效应。
使用gIVW、gEgger和Wald ratio方法，基于mQTL和IPF GWAS数据，对420,509个CpG位点进行EWMR分析，并进行敏感性分析。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
表观基因组-wide孟德尔随机化结合多组学验证鉴定了特发性肺纤维化的表观遗传驱动因素。
Epigenome-wide Mendelian randomization with multi-omics validation identifies epigenetic drivers of idiopathic pulmonary fibrosis.
实验完整度
高包含EWMR、共定位、DMA、eQTM、转录组分析、TF预测、细胞模型和RNA-seq等多层次验证，明确了功能与机制。
主要模型
A549肺泡上皮细胞系（TGFβ1诱导模型） IPF患者肺组织样本 IPF患者血液样本（mQTL数据） 对照肺组织样本
重点核对
EWMR分析中，CpG位点与IPF的因果效应（Bonferroni校正p<0.05/134,316） 共定位分析中，PP.H4>0.80的CpG位点 DMA中，肺组织DMPs（FDR校正p<0.05） eQTM分析中，CpG位点与基因表达的Pearson相关性（p<0.05） TGFβ1诱导的A549细胞模型中，MAN2A2和ZNF384的表达变化（n=3，p<0.05）
摘要
特发性肺纤维化（IPF）是一种病因不明且缺乏有效治疗的复杂疾病。虽然DNA甲基化已被认为与IPF发病机制有关，但表观遗传因素对IPF的组织特异性因果效应仍未确定。在这里，我们利用基于血液的420,509个CpG位点的甲基化数量性状位点和IPF的全基因组关联研究进行表观基因组-wide孟德尔随机化，以阐明CpG位点对IPF的因果效应。经Bonferroni校正后，共有452个CpG位点显示出对IPF风险的推定因果效应。其中，13个CpG位点与IPF相关的遗传因素显示出强烈的共定位证据。具体而言，MAN2A2和TRIM27基因内CpG位点的DNA甲基化在IPF肺和对照组之间显示出显著差异，并与这些基因在肺组织中的mRNA表达改变相关。根据转录因子数据库，MAN2A2中的CpG位点是ZNF384的结合位点。TGFβ1诱导的肺泡上皮细胞中的RNA测序证实，与对照组相比，MAN2A2和ZNF384的表达显著降低。总的来说，我们的研究提示MAN2A2内的DNA甲基化与IPF风险之间存在推定因果关联，其中肺特异性DNA甲基化可能扰乱ZNF384与MAN2A2之间的相互作用，揭示了这些基因在IPF发病机制中的新作用。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- DNA甲基化是否对特发性肺纤维化（IPF）风险具有因果效应，并识别潜在的因果CpG位点及其组织特异性功能。
- 核心机制
- MAN2A2基因内的DNA甲基化（cg01482407）可能通过破坏转录因子ZNF384与其启动子区域的结合，导致MAN2A2表达下调，从而促进IPF发病。
- 主要证据
- EWMR和共定位分析识别出13个与IPF相关的CpG位点；肺组织DMA和eQTM分析证实MAN2A2和TRIM27的甲基化与表达改变；TGFβ1诱导的A549细胞模型RNA-seq验证了MAN2A2和ZNF384表达下降。
- 研究意义
- 揭示了MAN2A2和ZNF384在IPF发病中的新作用，为IPF的预防和治疗提供了潜在靶点。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
共定位分析验证共享遗传信号
评估因果CpG位点与IPF风险位点是否共享遗传变异。
使用coloc.abf()函数，基于全部cis-mQTLs（无LD clumping）和IPF GWAS数据进行共定位分析。
肺组织差异甲基化分析（DMA）
验证因果CpG位点在IPF肺组织中的甲基化差异。
使用GEO数据集GSE175458的肺组织DNA甲基化数据（345例IPF，202例对照），通过ChAMP包进行DMA分析。
eQTM分析关联甲基化与基因表达
评估因果CpG位点的甲基化水平与靶基因表达的相关性。
使用GSE175458的甲基化数据和GSE175457的转录组数据，对DMA验证的CpG位点进行eQTM分析。
转录组分析验证基因表达差异
独立验证候选基因在IPF肺组织中的差异表达及其与肺功能的相关性。
使用GEO数据集GSE47460，比较IPF患者和对照肺组织中候选基因的表达，并分析其与肺功能（FVC, DLCO）的相关性。
转录因子结合预测和实验证据
探索因果CpG位点调控基因表达的潜在转录调控机制。
使用AnimalTFDB4.0预测CpG位点附近的TF结合位点，通过hTFtarget数据库查询ChIP-seq实验证据。
TGFβ1诱导的AEC细胞模型验证
在体外模型中验证MAN2A2和ZNF384的表达变化及其潜在功能。
使用TGFβ1处理A549细胞，通过RNA-seq评估MAN2A2、ZNF384和TRIM27的表达变化。
独立GWAS数据中的复制分析
验证发现结果在不同IPF病例定义中的稳健性。
使用独立临床诊断的IPF GWAS数据（Allen et al.），对发现的765个CpG位点进行复制EWMR分析。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
验证样本或模型特征
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 细胞培养 | DMEM/F-12培养基 | Gibco | C11330500BT |
| 胎牛血清 | Gibco | 1027016 | |
| 青霉素/链霉素 | Gibco | 15140-122 | |
| 药物处理 | 转化生长因子β1 | Sino Biological | 10804-HNAC |
| RNA测序 | NovaSeq6000测序平台 | Illumina | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| EWMR分析 | cis-mQTL的筛选条件（p<1e-8，r2<0.3，F>10）；IPF GWAS数据来源和样本量；Bonferroni校正的阈值。 阅读提示：Methods：Data sources，Epigenome-wide Mendelian randomization (EWMR) |
| 共定位分析 | 用于共定位的SNP列表（所有cis-mQTLs，无LD clumping）；先验概率设定；PP.H4的阈值。 阅读提示：Methods：Colocalization analysis |
| DMA | 使用的GEO数据集和样本量；DMA方法（ChAMP包）；DMP的FDR阈值。 阅读提示：Methods：Differential methylation analysis (DMA) |
| eQTM分析 | 配对的甲基化和基因表达数据集；用于分析的CpG位点筛选标准；Pearson相关的显著性阈值。 阅读提示：Methods：Expression quantitative trait methylation (eQTM) analysis |
| 细胞实验 | 细胞系（A549）；TGFβ1处理浓度（10 ng/mL）和时间（48小时）；样品重复次数（3次）。 阅读提示：Methods：Cell culture and treatment |