《Hypertension in Pregnancy》：A preeclampsia prediction model incorporating methylenetetrahydrofolate reductase, vitamin D receptor, and uterine artery pulsatility index
编辑推荐：
目的：本研究评估了亚甲基四氢叶酸还原酶（MTHFR）rs1801133和维生素D受体（VDR）FokI rs2228570多态性联合子宫动脉搏动指数（PI）对子痫前期（PE）的预测价值，并探讨了遗传风险背后的中介通路。
方法：开展了一项两阶段前瞻性研究。第一阶
目的：本研究评估了亚甲基四氢叶酸还原酶（MTHFR）rs1801133和维生素D受体（VDR）FokI rs2228570多态性联合子宫动脉搏动指数（PI）对子痫前期（PE）的预测价值，并探讨了遗传风险背后的中介通路。
方法：开展了一项两阶段前瞻性研究。第一阶段（n=73）通过靶向测序筛选PE相关单核苷酸多态性（SNPs）；第二阶段（n=195）验证候选位点并收集临床和超声指标。通过多因素logistic回归确定独立预测因子以构建列线图，并采用受试者工作特征（ROC）曲线、校准曲线和决策曲线分析评估列线图。通过中介分析量化血清生物标志物的贡献。
结果：MTHFR rs1801133的T等位基因在显性模型下与PE风险呈边缘性关联（OR=1.82，95% CI：0.96–3.44，P=0.067），而VDR rs2228570无显著关联（OR=1.74，95% CI：0.89–3.39，P=0.104）。孕前体质指数（BMI）≥25 kg/m2（OR=2.75，95% CI：1.45–5.21，P=0.002）和子宫动脉PI升高（OR=4.75，95% CI：1.38–16.41，P=0.014）是独立危险因素。列线图的曲线下面积（AUC）为0.711（95% CI：0.640–0.783），灵敏度为75.53%，特异度为59.41%。不对称二甲基精氨酸（ADMA）和8-异前列腺素F2α分别介导了MTHFR对PE风险效应的77.5%和26.0%。
结论：孕前BMI和子宫动脉PI是PE的独立预测因子。整合模型具有中等预测价值。ADMA介导的内皮功能障碍可能参与MTHFR相关的PE风险。
论文解读
一、研究背景与意义
子痫前期（preeclampsia，PE）是妊娠期最严重的高血压疾病之一，以妊娠20周后新发高血压伴多器官受累为特征。系统综述和荟萃分析报告中国妊娠期高血压疾病总体患病率为7.30%，其中PE占4.50%。PE可导致胎盘早剥、胎儿生长受限和早产等不良结局，同时增加母体远期心血管风险，是全球孕产妇和围产期死亡的重要贡献因素。因此，准确早期识别高危个体并及时干预，对降低母婴发病率至关重要。
当前PE的临床预测主要依赖常规因素如孕前体质指数（body mass index，BMI）、尿蛋白和肌酐，以及子宫动脉多普勒参数等影像学指标。子宫动脉搏动指数（pulsatility index，PI）反映中孕期子宫胎盘灌注状态，当胎盘血管重塑受损时升高，是PE的重要预测指标。然而，现有模型主要依赖常规临床变量，在很大程度上忽视了遗传易感性，导致预测性能欠佳，在精准医学中的应用价值有限。
PE的发病机制复杂，涉及遗传因素与环境因素的相互作用。亚甲基四氢叶酸还原酶（methylenetetrahydrofolate reductase，MTHFR）基因是叶酸代谢的关键调控因子，其多态性可通过影响酶活性改变同型半胱氨酸（homocysteine，Hcy）的代谢水平，从而诱导血管内皮功能损伤，已被证明与PE的遗传易感性密切相关。维生素D受体（vitamin D receptor，VDR）基因的FokI多态性通过调控维生素D信号通路影响血管平滑肌功能平衡和炎症反应，也可能参与PE的发病。尽管已有研究探讨了单一遗传位点或超声指标与PE的关联，但将遗传多态性与子宫动脉PI结合构建多因素整合PE风险预测模型的研究仍较匮乏，且缺乏探索性分析，难以明确基于多因素logistic回归的联合预测可行性及核心变量筛选方向。
为此，研究人员设计了一项两阶段前瞻性研究，旨在：（1）探讨MTHFR rs1801133和VDR FokI rs2228570与PE的关联；（2）检验其与子宫动脉PI的联合预测价值并探索构建列线图的可行性；（3）通过中介分析探索潜在的中介通路，为未来研究提供初步证据。该论文发表在《Hypertension in Pregnancy》期刊。
二、主要关键技术方法
研究人员采用两阶段前瞻性研究设计。第一阶段（探索队列，2022年10月至2023年4月）纳入自贡市妇幼保健院42例PE患者和31例对照，通过Illumina NovaSeq平台高通量靶向测序对8个候选SNP位点进行初步基因分型和关联分析。第二阶段（验证与建模队列，2023年5月至2025年7月）纳入101例PE患者和94例对照，采用竞争性等位基因特异性PCR（KASP）技术对候选位点进行验证性基因分型，并收集临床、实验室和超声指标。子宫动脉PI于孕20–24周由两名统一培训且不知晓分组的超声医师使用Samsung WS80A彩色多普勒超声诊断仪测量，取双侧PI平均值，通过组内相关系数（ICC）评估测量一致性。采用多因素logistic回归识别独立危险因素并构建列线图，通过Bootstrap重采样（1000次）进行内部验证，使用ROC曲线、校准曲线和决策曲线评估模型性能。中介分析采用Bootstrap抽样（5000次）量化血清生物标志物在遗传风险与疾病关联中的中介作用。
三、研究结果
基线特征：两阶段队列中PE组孕前BMI均显著高于同期健康对照组（第一阶段p=0.012；第二阶段p<0.001）。第二阶段队列中，PE组PI和纤维蛋白原（FIB）水平显著高于对照组（PI：p=0.044；FIB：p=0.037），血清25-羟基维生素D（25(OH)D）水平显著低于对照组（p<0.001），血清尿酸（UA）水平略高（p=0.011）。两组在年龄、初产妇比例、甘油三酯（TG）、低密度脂蛋白胆固醇（LDL-C）和Hcy水平方面无显著差异。
第一阶段遗传关联分析：在73例样本的探索队列中，MTHFR基因rs1801133（C677T）与PE风险显著相关，携带至少一个T变异等位基因（C/T或T/T）的孕妇与野生型C/C纯合子相比，患病风险显著增加（显性模型：OR=2.52，95% CI：1.00–6.40，p=0.050）。VDR基因rs2228570（FokI）显示潜在关联趋势（等位基因模型：OR=1.89，p=0.060）。其余6个SNP在各遗传模型中均无显著关联。所有SNP在对照组中的基因型分布均符合哈迪-温伯格平衡（Hardy-Weinberg equilibrium，HWE）。
KASP与高通量测序一致性验证：以高通量测序结果为“金标准”，对第一阶段73例样本的8个SNP位点进行平行检测和一致性分析，所有位点一致性百分比均高于97%，Cohen's Kappa值均大于0.95（范围：0.958–1.000），表明KASP分型技术具有极高的重复性和准确性。
第二阶段独立队列遗传关联分析：MTHFR rs1801133的突变T等位基因在病例组频率（46.0%）显著高于对照组（34.0%，p=0.016），显性模型差异显著（p=0.030），加性模型支持该关联（p=0.016），隐性模型未达显著水平（p=0.087）。VDR rs2228570的野生型C等位基因在病例组频率（62.4%）显著高于对照组（52.1%，p=0.041），隐性模型（C/C vs. C/T+T/T）差异显著（p=0.048），加性模型支持该关联（p=0.041），显性模型无显著关联（p=0.167）。
独立性与联合风险分析：单因素分析显示孕前BMI升高、子宫动脉PI升高、MTHFR rs1801133 T/T纯合突变基因型和VDR rs2228570 C/C野生型纯合基因型均与PE风险显著相关。所有候选预测因子的方差膨胀因子均低于2，适合进行多因素logistic建模。多因素logistic回归分析显示，孕前BMI≥25 kg/m2（OR=2.75，p=0.002）和子宫动脉PI升高（OR=4.75，p=0.014）是PE的独立危险因素。两个遗传位点虽未达统计学显著，但均显示风险增加趋势（rs1801133：OR=1.82，p=0.067；rs2228570：OR=1.74，p=0.104）。
预测模型性能验证与可视化：基于最终模型构建了整合孕前BMI、子宫动脉PI和两个遗传风险位点的列线图。模型的ROC曲线下面积为0.711（95% CI：0.640–0.783），灵敏度为75.53%，特异度为59.41%。经Bootstrap自抽样（1000次）内部验证后，校正后的一致性指数（Harrell's C-index）为0.659。校准曲线显示偏差校正后的预测曲线与理想对角线拟合良好，表明模型预测概率与实际观察概率基本一致。临床决策曲线分析显示，在阈值概率0.1至0.5范围内，使用该模型指导干预决策可获得高于“全部干预”或“全部不干预”策略的临床净收益。
MTHFR rs1801133基因型影响PE风险的机制初探：PE患者血清中反映内皮功能障碍的ADMA、反映氧化应激的8-异前列腺素F2α（8-iso-PGF2α）和炎症标志物超敏C反应蛋白（hs-CRP）水平均显著高于健康孕妇（均p<0.05）。在健康对照组中，携带MTHFR rs1801133 T等位基因的个体基线血清ADMA水平显著高于C/C基因型者（0.66 vs. 0.57 μmol/L，p=0.025）。中介分析显示，血清ADMA介导了遗传风险总效应的77.5%（间接效应=0.086，p=0.020），8-iso-PGF2α也发挥显著的部分中介作用（中介比例为26.0%，p=0.020），而hs-CRP和传统公认的Hcy均未显示显著中介效应。
四、讨论与结论
讨论总结：本研究通过两阶段设计，初步探索了整合遗传易感性（MTHFR rs1801133、VDR FokI rs2228570）和胎盘血流动力学指标（子宫动脉PI）对PE的预测价值，并构建了相应的临床列线图模型。更重要的是，通过中介分析提出了MTHFR基因风险效应中独立于Hcy的新机制通路。MTHFR rs1801133 T等位基因在两个独立队列中均显示与PE风险增加相关，符合显性遗传模式；VDR rs2228570 C等位基因（尤其是C/C野生型纯合基因型）与PE风险相关，倾向于隐性遗传模式。这是首次在同一PE人群中同时验证分别参与代谢（叶酸/Hcy）和免疫/发育（维生素D信号）通路的两个基因位点的风险效应，提示PE的遗传背景具有高度异质性和通路交叉性。中介分析未证实血清Hcy是该遗传变异与PE发病之间的关键中介因子，而ADMA的高中介比例提示MTHFR rs1801133 T等位基因可能通过影响体内ADMA水平间接干扰一氧化氮（NO）依赖性血管稳态，从而与PE风险增加相关。8-iso-PGF2α的中介效应提示氧化应激增强也参与其中，氧化应激与内皮功能障碍在PE发病中常相互强化。子宫动脉PI升高作为PE的独立危险因素，与现有病理生理学理解和大量临床证据高度一致，其中孕期（20–24周）PI评估窗口恰是具备临床干预潜力的“警示期”。研究局限性包括：探索阶段纳入的8个候选SNP中6个未显示显著关联，可能确实与PE无关或样本量不足未能检测弱效应；作为单中心研究可能存在选择偏倚和结果外推受限；模型验证仅采用内部Bootstrap重采样，缺乏外部独立队列验证；机制探索仅检测有限标志物，未全面覆盖MTHFR基因相关代谢通路；未考虑基因-基因和基因-环境交互作用。
结论翻译：本研究初步证实了MTHFR rs1801133和VDR FokI rs2228570基因多态性与PE风险的关联，明确了孕前BMI和子宫动脉PI的独立预测价值，所构建的多因素logistic联合预测模型为PE的早期风险筛查提供了新的探索方向。同时，本研究初步揭示了ADMA和氧化应激在MTHFR基因多态性所致PE风险中的中介作用，为理解PE的发病机制提供了新线索。
