H. pylori外膜蛋白抗体反应通过MHC/HLA-DMA甲基化影响胃癌风险的多层面孟德尔随机化证据整合研究

《Genes》：H. pylori Outer Membrane Protein Antibody Response, MHC/HLA-DMA Methylation, and Gastric Cancer Risk: A Multi-Layer Evidence Integration Study Using Mendelian Randomization

　　

背景：幽门螺杆菌（Helicobacter pylori，H. pylori）感染是胃癌的明确危险因素；然而，关于抗原特异性抗体性状、宿主甲基化与胃癌风险之间关联的证据仍然有限。方法：采用正向双样本孟德尔随机化（Mendelian Randomization，MR）设计，研究人员整合了甲基化数量性状位点（methylation Quantitative Trait Loci，mQTLs）、外部队列、癌症基因组图谱-胃腺癌（The Cancer Genome Atlas–Stomach Adenocarcinoma，TCGA-STAD）数据及功能注释，通过三阶段框架（抗体性状筛选、CpG位点优先排序和外部验证）评估了从宿主抗体反应经DNA甲基化到胃癌的两步通路。结果：在六种H. pylori抗原中，针对外膜蛋白（Outer Membrane Protein，OMP）的抗体反应与胃癌风险的关联最强（OR = 1.186，95% CI：1.079–1.304，原始p = 4.08 × 10^?4）。CpG水平优先排序鉴定出cg24290947（注释至HLA-DMA）为与胃癌易感性相关的甲基化位点，MR分析呈逆关联（beta = ?0.2347，95% CI：?0.3980至?0.0713，原始p = 0.00487）。外部探索显示，cg24290947甲基化在H. pylori阳性胃黏膜中升高；然而，在12对配对队列中其直接肿瘤-正常组织差异并不显著，尽管观察到区域性甲基化扰动。在三个队列中，HLA-DMA和MHC-II表达在肿瘤组织中均一致高于正常组织，而TCGA-STAD分析不支持cg24290947甲基化与HLA-DMA/MHC-II表达之间存在稳定的线性关系；六个免疫组成校正模型中三个显示名义显著性，但存在模型依赖性。这些信号富集于抗原加工与呈递通路。结论：cg24290947/HLA-DMA代表一个有待进一步研究的、产生假说的区域性甲基化候选位点，而非一个已确认的线性转录调控因子。MR证据基于欧洲血统汇总数据；其普适性和任何临床预测价值有待独立评估。

**论文解读**

**一、研究背景与问题**

胃癌是全球发病率和死亡率最高的恶性肿瘤之一。根据GLOBOCAN 2022数据，全球每年新发胃癌超96.8万例，死亡约66万例。该病起病隐匿，早期常无明显症状，多数患者确诊时已属晚期，五年生存率显著降低，而传统血清标志物如CEA、CA19-9和CA72-4在早期阶段敏感性和特异性有限，构成早期检测的重大技术瓶颈。

H. pylori是胃癌最重要的可控危险因素，被国际癌症研究机构列为I类致癌物。H. pylori抗原包括细胞毒素相关抗原A（Cytotoxin-associated antigen A，CagA）、空泡细胞毒素（Vacuolating cytotoxin，VacA）、伴侣蛋白GroEL（Chaperonin GroEL）和OMP等，其致癌潜能存在差异且具有人群特异性效应。研究表明，在东亚人群中OMP对非贲门胃癌风险的影响最强（相对风险增加2.66倍），超过传统毒力因子CagA（2.20倍）；而在欧洲人群中仅CagA和GroEL具有显著关联，OMP无显著效应，凸显了宿主-病原体相互作用的复杂性。

异常DNA甲基化，特别是启动子CpG岛（CpG island）高甲基化，在H. pylori诱导的胃黏膜细胞恶性转化中发挥关键作用。H. pylori感染是诱导胃黏膜异常DNA甲基化最重要的环境因素之一，且H. pylori诱导的异常DNA甲基化可在胃黏膜中形成"表观遗传场缺陷"（epigenetic field defect），某些位点的甲基化水平在H. pylori根除后仍持续升高。

尽管观察性研究已揭示H. pylori感染、DNA甲基化和胃癌之间的密切关联，但现有证据存在两个主要局限：其一，多数研究将H. pylori感染作为同质性暴露处理，未能区分特定毒力抗原对宿主表观遗传修饰的差异化效应；其二，从抗原特异性免疫反应经DNA甲基化至胃癌的完整中介链尚未被系统建立。为此，本研究采用正向双样本MR方法，以宿主对H. pylori抗原的抗体反应相关遗传变异作为暴露工具变量，以mQTLs作为中介工具变量，系统评估抗体性状→DNA甲基化→胃癌的两步通路。

**二、主要技术方法（约250字）**

研究人员构建了四阶段分析框架：第一阶段对六种H. pylori抗体性状与胃癌风险进行MR分析，确定OMP为主导抗体性状。第二阶段以OMP为锚点，整合GoDMC mQTL数据（基于32,851名参与者外周血DNA甲基化数据，覆盖420,509个甲基化位点）进行CpG发现，经探察性层级筛选保留132个候选CpG位点。第三阶段对候选CpG进行CpG-胃癌MR分析以评估其遗传关联。第四阶段利用外部基因表达综合数据库（Gene Expression Omnibus，GEO）甲基化数据集（GSE99553、GSE164988、GSE30601）和表达数据集（GSE29272含134对肿瘤-正常样本、GSE54129、GSE66229）及TCGA-STAD 354例配对甲基化-表达样本，结合定向功能注释进行多层面证据整合。MR分析以逆方差加权法（Inverse Variance Weighted，IVW）为主估计量，辅以MR-Egger、加权中位数法等敏感性分析，并以Cochran's Q检验评估异质性，MR-Egger截距检验评估方向性水平多效性。

**三、主要研究结果**

**3.1 OMP抗体被优先确定为主导H. pylori抗体性状**

在六种H. pylori抗体性状中，OMP与胃癌风险的关联最为稳健（IVW OR = 1.186，95% CI：1.079–1.304，原始p = 4.08 × 10^?4），经BH-FDR（Benjamini-Hochberg假发现率）校正后仍显著（校正p = 0.0020）。其余抗体性状VacA、GroEL、过氧化氢酶（Catalase）和CagA均未观察到显著关联（原始p > 0.05）；UREA因协调后无可用的单核苷酸多态性（Single Nucleotide Polymorphism，SNP）而未能获得有效IVW估计值。敏感性分析未检测到显著异质性（Cochran's Q p = 0.2159）或方向性水平多效性（MR-Egger截距p = 0.2247；MR-PRESSO全局检验p = 0.208）。

**3.2 OMP相关CpG发现鉴定出MHC/HLA区域候选位点**

在全局CpG水平共检出652个OMP相关CpG，其中519个经多重校正后达到FDR < 0.05。经探察性层级标准筛选后保留132个候选位点，进入两个非互斥评估分支：严格的HLA/MHC来源注释CpG分支（3个）和CpG-胃癌MR评估分支（最终整合为7个CpG展示），两组仅在cg24290947/HLA-DMA处重叠。

**3.3 CpG-胃癌MR优先确定cg24290947/HLA-DMA**

在7个候选CpG中，cg24290947/HLA-DMA显示最强的遗传关联（beta = ?0.2347，95% CI：?0.3980至?0.0713，原始p = 0.0049），敏感性分析未发现显著异质性（Cochran's Q p = 0.5054）或方向性水平多效性（MR-Egger截距p = 0.7320）。其余6个CpG均未达到FDR < 0.05的显著性水平。

**3.4 外部甲基化数据集提供H. pylori暴露、直接位点和区域性甲基化证据**

GSE99553显示H. pylori阳性胃黏膜中cg24290947甲基化升高（delta beta = 0.0623，p = 0.0052）。GSE164988的12对配对样本中该位点直接肿瘤-正常差异不显著（delta beta = 0.0121，p = 0.6772）。GSE30601的平台未覆盖cg24290947，但在HLA-DMA/HLA-DMB/BRD2区域，6个可检测CpG中有3个达到FDR < 0.05。

**3.5 外部表达数据集支持胃癌中MHC-II抗原呈递背景**

三个外部GEO表达队列中，HLA-DMA表达在肿瘤中均高于正常组织（GSE29272：delta = 0.225，FDR = 0.0280；GSE54129：delta = 0.443，FDR = 0.0015；GSE66229：delta = 0.052，FDR = 0.0453）。MHC-II平均z评分标志也呈现一致的肿瘤高于正常趋势（GSE29272：delta = 0.404；GSE54129：delta = 1.047；GSE66229：delta = 0.153），均达到FDR < 0.05。

**3.6 TCGA-STAD联合分析不支持线性甲基化-表达偶联**

在354例TCGA-STAD原发肿瘤样本中，未校正分析显示cg24290947甲基化与HLA-DMA表达无显著相关（Spearman rho = ?0.017，p = 0.7485），与MHC-II标志亦无显著相关（rho = 0.043，p = 0.4146）。在六个组织组成校正模型中，ImmuneScore校正后甲基化项与HLA-DMA表达呈名义逆关联（beta = ?1.5743，p = 0.0015），但该结果未在不同组成协变量间一致复制：StromalScore校正后不显著（p = 0.156），TumorPurity校正后的HLA-DMA关联在排除高影响样本后不再达到名义显著性（p = 0.0882）。

**3.7 整合证据分类支持cg24290947/HLA-DMA为区域性甲基化锚点**

七种OMP相关候选CpG（cg24290947、cg17239008、cg04317588、cg13047157、cg21042276、cg09474017、cg08491668）及其注释基因（HLA-DMA、HLA-DQB1、TNXB、BRD2）在五层证据体系（"++"、"+"、"±"、"–"、"NA"）中接受评估，cg24290947/HLA-DMA在多个层级获得一致支持，但其TCGA-STAD甲基化-表达偶联证据仅为模型依赖的临界支持。

**3.8 定向功能注释支持抗原呈递免疫背景**

以MHC/HLA区域候选基因和MHC-II相关基因为输入的GO/KEGG/Reactome富集分析显示，该基因集显著富集于抗原加工与呈递相关通路：KEGG抗原加工与呈递（FDR = 7.05 × 10^?19）、肠道IgA产生免疫网络（FDR = 8.13 × 10^?16）、Reactome MHC II类抗原呈递（FDR = 1.56 × 10^?12）以及GO外源性肽抗原经MHC II类加工呈递（FDR = 7.03 × 10^?14）。

**四、讨论总结**

研究通过MR分析首次确定OMP为与胃癌风险关联最强的H. pylori抗体靶标（OR = 1.186），并通过多阶段分析鉴定出注释至HLA-DMA的cg24290947为潜在胃癌风险相关甲基化位点。外部甲基化探索显示该位点在H. pylori阳性胃黏膜中甲基化升高，且HLA-DMA/MHC邻近区域在胃癌中存在甲基化扰动。三个独立胃癌队列中HLA-DMA和MHC-II标志均呈肿瘤高表达方向。然而，TCGA-STAD未校正和组成校正分析均不支持cg24290947甲基化与HLA-DMA表达之间存在简单线性调控关系，且OMP抗体性状→胃癌（正向）、OMP抗体性状→cg24290947甲基化（正向）和cg24290947→胃癌（逆关联）的方向模式不支持传统正向中介通路，这可能反映该区域连锁不平衡（Linkage Disequilibrium，LD）对OMP遗传工具变量的影响。因此，cg24290947更适合被解读为区域性甲基化信号而非直接中介因子。HLA-DMA编码HLA-DM复合物的α亚基，其核心功能是作为"肽加载催化剂"移除MHC II类分子上的II类相关恒定链肽（Class II-associated Invariant Chain Peptide，CLIP），促进高亲和力抗原肽的高效加载。区域性甲基化扰动和MHC II类通路基因在肿瘤组织中的整体表达升高提示，cg24290947甲基化可能反映该位点响应H. pylori相关免疫压力的更广泛表观基因组重塑。研究存在若干局限：MR和mQTL分析依赖公开汇总数据，mQTL工具变量来自外周血而非胃黏膜，研究结果主要基于欧洲人群，外部GEO数据集来自不同平台且无法合并，HLA/MHC区域未进行正式共定位分析等。

**五、研究结论**

本研究在cg24290947/HLA-DMA处鉴定出一个OMP抗体相关甲基化候选位点，并发现其与胃癌风险呈逆遗传关联。外部组织证据支持区域性背景，但未建立与HLA-DMA的因果调控关系。遗传学结果基于欧洲血统汇总数据，仍受限于工具变量有效性、组织特异性、多重性和区域LD等局限。该候选位点有待独立的生物学和跨血统研究验证；其临床预测价值尚未得到证实。