研究队列建立与样本采集
建立儿童起病多发性硬化患者与健康对照的外周血DNA样本队列，用于全基因组甲基化分析。
多中心回顾性研究纳入175名意大利受试者（122名PedMS患者，53名健康对照），从外周血提取gDNA。PedMS患者入组标准为起病年龄<18岁、符合Krupp标准诊断、疾病持续时间<4年。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
儿童多发性硬化的表观遗传学洞察：DNA甲基化突显EB病毒感染与宿主免疫网络的参与
Epigenetic Insights in Pediatric Multiple Sclerosis: DNA Methylation Highlights an Involvement of Epstein-Barr Virus Infection and Host Immune Networks
实验完整度
高全血DNA甲基化芯片分析175例样本，含敏感性亚组和内部验证，并行网络与通路富集验证。
主要模型
122例儿童起病多发性硬化患者外周血DNA 53例儿童健康对照外周血DNA 未治疗PedMS亚组（n=21） LETA及其他DMT治疗亚组（n=54）
重点核对
DMRs检出标准：最少7个连续DMPs，DMR间最小间隔（lambda值）1,000 bp，FDR调整p值0.05，mean Δβ阈值0.02 内部验证采用1/3–2/3分区，发现样本n=59，验证样本n=116，DMRcate中单个CpG位点检测阈值放宽至FDR<0.1 敏感性分析1：仅未治疗PedMS患者（n=21）与HC比较 敏感性分析2：未治疗+LETA+其他DMT治疗患者（n=54）与HC比较 线性模型协变量：免疫细胞估计比例矩阵的前3个主成分（捕获86%细胞方差）和性别 网络分析保留交互证据评分>0.95的链接，使用1步邻域法，子网络165个节点和680条边
摘要
背景与目标：儿童起病多发性硬化（PedMS）因环境暴露与疾病起病之间间隔较短，为研究MS发病机制提供了独特机会。本研究旨在识别在疾病早期介导环境触发因素与基因组因素相互作用的差异甲基化区域（DMRs）。方法：在这项多中心研究中，我们使用Infinium Methylation EPIC Array v2分析了175名意大利受试者（122名PedMS患者，53名健康对照）的外周血DNA。所有患者均在起病4年内入组。生物信息学流程包括批次效应校正、归一化以及针对性别和推断细胞类型比例校正的差异甲基化分析。使用DMRcate检测DMRs。通过NDEx工具将DMR重叠基因投射到STRING相互作用组进行网络分析，随后使用gProfiler进行通路富集分析。为确保稳健性，我们对未治疗和低效治疗（LETA）亚组进行了内部1/3–2/3验证和敏感性分析。结果：我们的分析揭示了55个DMRs，其中28个影响启动子区域。主要信号涉及CHI3L2、LPCAT1、CD19、PARP1、CPT1B-CHKB和IGF2R等基因，包括若干MS相关位点。映射到CHI3L2/DENND2D和LPCAT1的信号高度稳健，在所有验证和敏感性亚组中持续存在。网络分析识别出一个由165个节点组成的子网络，其中Epstein-Barr病毒（EBV）感染成为排名最高的富集通路（p = 3.23 *10−40）。高中心性调控枢纽包括PARP1和IGF2R，后者与已知印记控制区域重叠。讨论：我们的发现识别出早期PedMS中一种独特的表观遗传特征，提示DNA甲基化变化并非随机，而是汇聚于高中心性调控枢纽。与EBV相关通路和印记位点的稳定关联强化了这一假说：表观遗传重塑是MS起病中早期环境暴露与遗传易感性之间的关键界面。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 识别儿童起病多发性硬化早期阶段中介环境触发因素与基因组因素相互作用的差异甲基化区域（DMRs）。
- 核心机制
- DNA甲基化变化并非随机，而是汇聚于高中心性调控枢纽（如PARP1和IGF2R），表观遗传重塑作为早期环境暴露（如EBV感染）与遗传易感性之间的关键界面参与PedMS起病。
- 主要证据
- 175名意大利受试者（122 PedMS，53 HC）外周血DNA的Infinium Methylation EPIC Array v2分析识别55个DMRs；CHI3L2/DENND2D和LPCAT1信号在所有验证和敏感性亚组中持续存在；网络分析识别EBV感染为排名最高富集通路（p = 3.23 *10−40），PARP1和IGF2R为高中心性枢纽。
- 研究意义
- 研究提示区域DNA甲基化分析能够在疾病最早阶段揭示可重复的表观遗传模式，并强调未来需整合靶向测序及暴露数据以建立表观遗传特征与环境暴露之间的功能联系。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
全基因组DNA甲基化芯片检测
获取全基因组DNA甲基化谱，用于识别差异甲基化探针和区域。
亚硫酸氢盐转化后的DNA杂交到Infinium Methylation EPIC Array v2.0，覆盖超过930,000个甲基化位点（参考基因组hg38），使用Illumina NextSeq 550测序仪测量信号强度。
数据预处理与细胞类型估计
校正技术偏差和细胞异质性，为差异甲基化分析提供校正后的甲基化值。
使用ChAMP v2.21.1进行预处理，β-mixture quantile normalization校正II型探针偏差，ComBat方法校正载玻片水平批次效应。使用DNAmAge calculator估计免疫细胞比例（CD4+ T细胞、CD8+ T细胞、B细胞、NK细胞、单核细胞、粒细胞），并生物信息学推断性别。
差异甲基化区域识别
识别PedMS患者与健康对照之间协调且统计学显著的DNA甲基化差异区域。
在甲基化M值矩阵上拟合多元线性回归模型，使用RnBeads v2.21.2中的limma方法。使用DMRcate方法（针对EPICv2更新）检测DMRs，重新映射4,130个交叉杂交探针，过滤重复探针。
网络分析与通路富集
探索DMR相关基因的相互作用网络及富集的生物学通路，识别关键调控枢纽。
使用NDEx工具将DMR重叠基因投射到STRING v12.0参考相互作用组，保留交互证据评分>0.95的链接（来自蛋白-蛋白相互作用、共表达和数据库三个证据域），过滤后相互作用组包含10,430个节点和58,559条边。使用1步邻域法探索网络，使用CentiScaPe估计节点中心性，使用gProfiler进行通路富集分析。
功能基因组注释
将顶级DMRs整合到功能基因组数据中，评估其调控潜力。
将顶级DMRs与ENCODE项目（Phase 4）功能基因组数据整合，包括Broad Institute生成的H3K27ac组蛋白乙酰化轨迹（使用GM12878淋巴母细胞样细胞系作为B淋巴细胞免疫分析的金标准参考），并使用SCREEN portal访问候选顺式调控元件多细胞注册库（cCREs，v3）。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| DNA提取与定量 | NanoDrop 2000分光光度计 | Thermo Scientific | -- |
| Quant-iT PicoGreen双链DNA检测试剂盒 | ThermoFisher Scientific | -- | |
| Qubit双链DNA高灵敏度检测试剂盒 | Life Technologies | -- | |
| Qubit双链DNA宽范围检测试剂盒 | Life Technologies | -- | |
| 亚硫酸氢盐转化 | EZ DNA甲基化试剂盒 | -- | -- |
| 全基因组DNA甲基化分析 | Infinium甲基化EPIC芯片v2.0 | Illumina | -- |
| Illumina NextSeq 550测序仪 | Illumina | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 研究队列与样本 | PedMS患者122例、健康对照53例；PedMS入组标准为起病年龄<18岁、符合Krupp标准、疾病持续时间<4年；对照组排除系统性自身免疫病。 阅读提示：请回看Methods的Study Population小节，核对入组/排除标准及样本量。 |
| DNA甲基化芯片检测 | Infinium Methylation EPIC Array v2.0；500 ng全血gDNA经EZ DNA Methylation Kit转化；探针过滤标准（检测p>0.01、bead count<3、非CpG、SNP重叠）。 阅读提示：请回看Methods的DNA Extraction and Bisulphite Treatment of DNA和Genome-Wide DNA Methylation Analysis小节。 |
| 差异甲基化分析 | 线性模型协变量为免疫细胞比例前3个主成分和性别；DMRcate参数：最少7个连续DMPs、lambda=1,000 bp、FDR<0.05、mean Δβ=0.02。 阅读提示：请回看Methods的Differential Methylation Analysis and Functional Annotation小节。 |
| 内部验证与敏感性分析 | 内部验证1/3–2/3分区，发现n=59、验证n=116，DMRcate阈值FDR<0.1；敏感性分析1（未治疗n=21）和敏感性分析2（未治疗+LETA+其他DMT n=54）。 阅读提示：请回看Results中关于敏感性分析和内部验证的描述，以及Methods中DMR检测参数。 |
| 网络分析与通路富集 | STRING v12.0、交互证据评分>0.95、1步邻域法、CentiScaPe、gProfiler；子网络165个节点和680条边。 阅读提示：请回看Methods的Network Analysis小节。 |