《Neurobiology of Pain》：Next-generation DNA methylation sequencing: loci and regions differentially methylated in adolescents with chronic postsurgical pain
编辑推荐：
DNA甲基化（DNAm）可能将环境暴露与疼痛易感性联系起来。研究人员（使用新一代酶促甲基化测序（NG EM-seq））检测了与慢性术后疼痛（CPSP）相关的术前DNAm模式，CPSP是接受肌肉骨骼手术的青少年中一个主要的生物心理社会问题，招募自六个研究中心（n
DNA甲基化（DNAm）可能将环境暴露与疼痛易感性联系起来。研究人员（使用新一代酶促甲基化测序（NG EM-seq））检测了与慢性术后疼痛（CPSP）相关的术前DNAm模式，CPSP是接受肌肉骨骼手术的青少年中一个主要的生物心理社会问题，招募自六个研究中心（n=166；年龄15.5岁（IQR 14.3–16.8）；51%男性）。疼痛和功能测量被纵向追踪。CPSP定义为术后两个月后疼痛>3/10。研究人员识别出289个差异甲基化位点（DML；HB p<0.1），这些位点在自助法后满足稳健性标准，以及20个与CPSP相关的区域（DMR）（FDR p<0.05）。使用CPSP特异性背景，DML注释到57个基因，富集于生物过程（细胞对cAMP的反应（p=0.001）、神经元动作电位（p=0.01）、钠离子跨膜转运（p=0.03）、钾离子跨膜转运（p=0.035）和突触后膜电位调节（p=0.041））、分子功能（细胞内cAMP激活的阳离子通道活性（p=0.001）、cAMP结合（p=0.003）、电压门控钠通道活性（p=0.01）和电压门控钾通道活性（p=0.02））以及细胞组分（HCN通道复合体（p=0.002）、突触前活性区膜（p=0.006）和突触后膜（p=0.049））。DMR未富集基因通路，但显示了转录因子基序富集于应激反应（TGF-β/SMAD）、神经元分化/可塑性（NR4A/TEAD）、炎症/免疫（STATs）和昼夜节律信号（RORA/RORB/NPAS2）。研究结果表明与CPSP相关的潜在DNAm特征，并对其在CPSP发展中的可能作用提出了假设。
**论文解读：新一代DNA甲基化测序揭示青少年慢性术后疼痛的表观遗传学特征**
**研究背景与问题**
慢性术后疼痛（CPSP）是儿童和青少年术后常见且严重的并发症，影响11%至60%的儿科手术患者，显著影响其情绪、身体和社会功能，并带来巨大经济负担。尽管CPSP具有遗传调控基础，但遗传关联研究缺乏可重复性，凸显了对易感性理解的不足。环境因素（如父母照顾、社会经济地位和心理社会因素）可通过表观遗传机制（如DNA甲基化（DNAm））修饰基因表达，从而影响疼痛反应和CPSP的发展。由于全基因组DNAm覆盖的局限性，既往研究可能遗漏了与CPSP相关的关键位点，特别是基因间区。因此，有必要采用更全面的测序方法，在更大的多中心青少年队列中，系统探索术前DNAm模式与CPSP的关联，并评估其调控作用。
**研究目的**
本研究旨在通过新一代酶促甲基化测序（NG EM-seq）技术，在多中心、接受不同肌肉骨骼手术的青少年队列中，寻找与CPSP相关的术前全基因组DNAm差异位点和区域，并揭示其潜在的生物学功能。
**主要技术方法**
1. **样本队列与临床评估**：在美国六个儿科中心招募343名接受脊柱融合术或Nuss手术的青少年（8岁以上），排除自身免疫病、恶性肿瘤等。术前采集全血样本，纵向追踪术后2-6个月和10-12个月的疼痛、功能和生活质量，CPSP定义为术后2个月后疼痛评分>3/10。最终纳入166名有完整甲基化数据和随访数据的受试者（中位年龄15.5岁，51%男性，64%为CPSP）。
2. **DNA甲基化测序与数据处理**：使用NEB Next酶促甲基化测序试剂盒和Illumina NovaSeq 6000平台进行全基因组测序（150 bp双端 reads）。通过Bismark比对到GRCh38/hg38基因组并调用甲基化。质量控制包括：排除覆盖度<7×、位于性染色体或与已知SNP重叠的CpG位点，保留在≥90%样本中存在的位点。
3. **差异甲基化分析**：采用两步法：首先用LIMMA线性模型筛选与CPSP相关的CpG位点（p<0.05），然后对候选位点进行逻辑回归，校正临床协变量（批次、术前疼痛等），以FDR<0.1为显著阈值。通过1000次非参数自助法验证稳健性（方向一致性>99%，p值稳定性>70%，置信区间排除零值）。使用DMRcate识别差异甲基化区域（DMR）（FDR<0.05）。
4. **功能富集分析**：使用topGO进行基因本体（GO）富集分析，使用monaLisa进行转录因子（TF）基序富集分析，背景基因集为CPSP特异性基因集（1305个候选基因和31个训练基因）。
**研究结果**
### 2.1. 受试者特征
共166名受试者纳入分析。CPSP组（n=108）与non-CPSP组（n=58）在年龄、性别、手术类型上可比，但CPSP组术前疼痛更强、功能残疾更严重（p<0.1）。术后2-6个月和10-12个月，CPSP组功能残疾指数（FDI）更高，生活质量（PedsQL）更低（p<0.001）。
### 2.2. 差异甲基化分析
#### 2.2.1. DNA甲基化预处理
成功提取所有术前血样DNA，测序质量高（Phred评分中位数36）。经过质量控制，保留44,255个CpG位点用于下游分析。通过EpiDISH估算并校正个体细胞类型比例。
#### 2.2.3. CpG甲基化谱与质量控制
技术重复显示高重现性（7×覆盖度下平均R2=0.9155）。排除与SNP重叠、性染色体及低覆盖度位点后，最终分析位点数为44,255。
### 2.3. 与CPSP相关的差异甲基化位点（DMLs）
#### 2.3.1. CpG位点主要位于基因间/内含子区域
第一阶段筛选出2210个CpG位点（p<0.05），其中37%位于基因间区，30%位于内含子区，约69%为低甲基化。
#### 2.3.2. 校正协变量后与CPSP相关的差异甲基化CpG位点主要为低甲基化
第二阶段逻辑回归并校正批次和术前疼痛后，646个CpG位点显著（FDR<0.1）。自助法验证后，289个CpG位点满足稳健性标准，其中70.2%为低甲基化（平均β差异=-0.04），主要位于基因间/内含子区。
### 2.4. 使用CPSP特异性背景集的基因富集分析
289个CpG位点注释到57个独特基因，包括32%非编码RNA、33%蛋白编码基因和32%假基因。GO富集分析显示显著富集于：
- **生物过程（BP）**：细胞对cAMP的反应（p=0.001）、神经元动作电位（p=0.01）、钠离子跨膜转运（p=0.03）、钾离子跨膜转运（p=0.035）、突触后膜电位调节（p=0.041）。
- **分子功能（MF）**：细胞内cAMP激活的阳离子通道活性（p=0.001）、cAMP结合（p=0.003）、电压门控钠通道活性（p=0.01）、电压门控钾通道活性（p=0.02）。
- **细胞组分（CC）**：HCN通道复合体（p=0.002）、突触前活性区膜（p=0.006）、突触后膜（p=0.049）。
KEGG和Reactome通路未发现显著富集。
### 2.5. 差异甲基化区域（DMR）分析
#### 2.5.1. 显著DMR的基因和基序富集分析
DMRcate识别出20个显著DMR（HMFDR<0.05），注释到7个独特基因（如ZNF407、CWH43、EMBP1等）。基因富集分析无显著结果，但转录因子（TF）基序富集分析显示显著过表达，涉及：
- **应激反应**：SMAD3、SMAD4、SMAD5（TGF-β/SMAD通路）。
- **神经元分化/可塑性**：NR4A2::RXRA、TEAD2、TEAD4。
- **炎症/免疫**：STAT1、STAT4、STAT5A、STAT6。
- **昼夜节律信号**：RORA、RORB、NPAS2。
**讨论与结论**
**讨论总结**：本研究首次在多中心青少年CPSP队列中使用全基因组酶促甲基化测序，识别出289个主要位于调控区域的差异甲基化位点（DML）和20个差异甲基化区域（DMR），揭示CPSP与广泛的低甲基化模式相关。DML富集于神经元和离子通道相关功能，提示表观遗传调控可能通过影响神经元兴奋性、突触可塑性和中枢敏化参与急性疼痛向慢性的转变。DMR虽未显示基因通路富集，但转录因子基序富集于应激、炎症、昼夜节律和神经元可塑性通路，为CPSP的机制研究提供了新假设。研究局限性包括：样本量有限、测序深度较低、丢失随访的偏倚、以及DML注释主要位于非编码区域，限制了直接基因水平解读。未来需在独立队列中验证，并进行多组学整合和功能实验。
**结论翻译**：这项DNA甲基化测序研究在多中心儿科队列中首次提供了全基因组DNA甲基化差异与CPSP相关的证据。这些发现支持DNAm作为CPSP生物学贡献因素，并为探索表观遗传调控如何促进儿童疼痛持续性的机制研究打开了大门，尽管在建立任何临床或风险分层作用之前，需要重复和正式评估预测性能。
