心房颤动患者口服盐酸贝普地尔校正QT间期敏感性的遗传学研究

《Journal of the American Heart Association》：Genetic Investigation of Corrected QT Interval Sensitivity to Oral Bepridil Hydrochloride Hydrate in Patients With Atrial Fibrillation

　　

盐酸贝普地尔（Bepridil Hydrochloride Hydrate）是一种用于治疗阵发性心房颤动（Paroxysmal Atrial Fibrillation, PAF）的抗心律失常药物，但其应用受到校正QT间期（QTc）延长及心室性心律失常风险的限制。本研究旨在识别与导管消融术后贝普地尔相关QTc延长相关的遗传变异。研究人员纳入了523例2013年11月至2023年3月期间于广岛大学接受贝普地尔治疗的病例，其中445例纳入全基因组关联研究（Genome‐Wide Association Study, GWAS），78例作为独立复制队列。QTc于治疗启动后1个月测定，QTc延长定义为男性≥450 ms、女性≥460 ms。采用加性（Additive）、显性（Dominant）及隐性（Recessive）遗传模型进行逻辑回归分析。GWAS中显示提示性关联（P<1.00×10^?6）的变异在复制队列中评估，随后进行随机效应荟萃（Meta）分析，全基因组显著性定义为P<5.00×10^?8。质量控制后共443例纳入GWAS分析，其中122例（28%）出现QTc延长。隐性模型下，FSHR（Follicle‐Stimulating Hormone Receptor，卵泡刺激素受体）基因内含子变异rs12622919显示提示性关联。在合并发现队列与复制队列的随机效应荟萃分析中，rs12622919达到全基因组显著性。探索性性别分层分析提示女性中存在潜在关联（P=0.026），但未观察到统计学显著的基因型‐性别交互作用。纳入rs12622919可改善QTc延长的判别能力（曲线下面积[Area Under The Curve, AUC] 0.77 [95% CI, 0.70–0.84]；Welch's t检验，P<0.001），且经DeLong检验改善仍具显著性（AUC 0.78 [95% CI, 0.74–0.84]；P=0.0078）。研究结论认为，FSHR变异rs12622919与贝普地尔相关QTc延长相关，可能作为个体风险评估的遗传生物标志物，但由于发现基于荟萃分析，尚需独立验证。

## 研究背景与目的

获得性长QT综合征以QT间期延长和潜在致死性室性心律失常为特征，可由药物（如抗心律失常药、抗生素）、电解质紊乱及心动过缓等外部因素触发。先天性长QT综合征由基因突变导致，钾通道功能丧失或钠/钙通道功能增益突变均可致心脏复极延迟。近年研究表明，遗传易感性不仅参与先天性长QT综合征，亦涉及获得性长QT综合征，患者可携带致病基因变异，凸显了遗传背景在获得性疾病中的重要性。多中心研究揭示，药物诱发尖端扭转型室性心动过速患者可携带主要候选基因（KCNQ1、KCNH2、SCN5A、KCNE1及KCNE2）变异，国际指南已建议对40岁以下药物诱发尖端扭转型室性心动过速患者进行基因检测。此外，常见遗传变异对QTc延长的影响亦被广泛研究，提示药物诱导QTc延长的遗传背景具有异质性。

盐酸贝普地尔是一种多离子通道阻滞剂，用于阵发性心房颤动的预防与治疗。尽管临床疗效确切，其已知风险包括QT间期延长及尖端扭转型室性心动过速等致死性室性心律失常。鉴于个体间QTc延长程度差异显著，识别易感患者至关重要。本研究旨在探究接受导管消融术后口服贝普地尔治疗的阵发性心房颤动患者中，贝普地尔诱导QTc延长的遗传背景，重点关注常见遗传变异即单核苷酸多态性（Single Nucleotide Polymorphism, SNP），通过全基因组关联研究识别影响QTc延长易感性的遗传因素。

传统GWAS多采用加性模型，但隐性模型亦具价值。例如，2型糖尿病研究中已发现若干仅隐性模型可检测的变异；药物代谢领域，CYP2C19多态性为典型范例，双等位基因功能缺失变异携带者酶活性显著降低、血栓事件风险增加。因此，本研究同时采用加性、显性及隐性模型进行全面分析，以期为个体化风险预测及更安全的抗心律失常策略提供依据。

## 主要技术方法

研究队列来源于广岛大学心血管内科2013年11月至2023年3月期间523例首次接受射频导管消融（Radiofrequency Catheter Ablation, RFCA）的阵发性心房颤动患者，术后次日启动口服盐酸贝普地尔（100 mg/d）治疗。基线QTc于RFCA入院时心电图测定，治疗后QTc于用药1个月门诊随访时测定，采用Bazett公式校正，由至少两名心律失常专科医师独立手动测量V5导联，差异<5%方采纳。排除标准涵盖年龄/QTc/体重指数（Body Mass Index, BMI）极端值（均值±4标准差）、非窦性心律、先天性长QT综合征、心源性猝死家族史及其他高危表型。

GWAS采用亚洲筛查芯片基因分型，经Minimac4软件及5.6K自建参考面板进行基因型推断，Info评分≥0.3的推断变异纳入分析。样本质控标准包括：样本检出率≥99%、基因型缺失率<1%、性别一致、排除近亲关系（PLINK软件计算PI_HAT>0.25）及基于1000基因组计划第三期数据主成分分析的东亚人群聚类异常值。遗传标记质控包括：基因分型效率≥99%、最小等位基因频率>1%及对照组Hardy-Weinberg平衡P值>1.00×10^?6。合格常染色体变异采用PLINK软件进行逻辑回归分析，调整年龄、性别、BMI及前3个主成分，分别采用加性、显性及隐性模型。结果以分位数-分位数图及曼哈顿图可视化。

复制队列78例于2022至2023年招募，候选变异采用多重聚合酶链反应-侵入试验（Invader assay）及QuantStudio 7 Flex实时聚合酶链反应系统分型。发现队列与复制队列采用METASOFT软件进行Han和Eskin改良随机效应荟萃分析，显著性阈值P<5.00×10^?8。

表达数量性状位点（Expression Quantitative Trait Loci, eQTL）分析利用基因型-组织表达（Genotype‐Tissue Expression, GTEx）数据库第8版，检测心房附属器、左心室、肝脏及肾脏组织。功能编码变异关联分析基于国家老年医学与老年学中心生物库4344例全基因组测序数据，采用ANNOVAR软件注释±1 Mb范围内变异，PLINK评估连锁不平衡（连锁不平衡系数r²）。预测模型构建采用75%训练集与25%测试集随机分割，逻辑回归模型比较含（年龄、性别、BMI、基线QTc）与不含rs12622919的模型性能，1000次Bootstrap重采样评估稳健性，AUC比较采用DeLong检验。

## 研究结果

**GWAS参与者人口学数据**：445例患者中，122例发生QTc延长，323例未发生。QTc延长组年龄较高、BMI较低，基线及治疗后QTc均较长。血液学指标组间无显著差异。Ⅰa类抗心律失常药物使用在非QTc延长组更常见，其他药物使用无显著差异。

**GWAS结果**：443例通过严格质控，3个遗传模型量化-量化图显示无基因组膨胀（基因组膨胀因子1.00–1.01）。3个变异达提示性显著：加性模型rs2185830、加性及显性模型rs957185、隐性模型rs12622919，比值比（Odds Ratio, OR）2.32–7.92，均与风险增加相关。

**复制研究与荟萃分析**：3个提示性变异于独立复制队列分型。rs12622919显示名义显著性（P=0.024），但未达Bonferroni校正阈值。发现队列与复制队列荟萃分析中，rs12622919达全基因组显著性（OR 3.88 [95% CI, 2.44–6.11]；P=8.52×10^?9）。

**rs12622919的eQTL效应**：GTEx数据库分析显示，rs12622919风险等位基因与左心室PPP1R21（Protein Phosphatase 1 Regulatory Subunit 21）表达降低显著相关（P<0.05），肝脏及肾脏无显著关联。

**rs12622919与功能变异的关联**：±1 Mb范围内141个功能变异中，无编码变异与rs12622919存在强连锁不平衡（r²>0.8）。

**rs12622919与临床因素的关联**：线性回归显示rs12622919与QTc延长显著相关（β=0.25 [95% CI, 0.072–0.43]；P<0.05）。性别分层分析仅女性显著（β=0.39 [95% CI, 0.046–0.74]；P<0.05），但GWAS中无显著基因型-性别交互作用（隐性模型交互项P=0.77），该发现属探索性。

**rs12622919预测QTc延长的潜力**：含rs12622919的模型AUC达0.77（95% CI, 0.70–0.84），显著高于不含该变异模型（0.73 [95% CI, 0.65–0.81]；Welch's t检验P<0.001）。DeLong检验验证差异显著（AUC 0.78 [95% CI, 0.74–0.84] vs. 0.75 [95% CI, 0.70–0.79]；P=0.0078）。性别分层模型中，男性纳入rs12622919后AUC为0.76（95% CI, 0.66–0.84），女性为0.80（95% CI, 0.66–0.93）。

## 讨论与结论

本研究采用加性、显性及隐性模型进行GWAS，通过隐性模型发现FSHR基因内含子变异rs12622919与贝普地尔治疗期间QTc延长相关。性别分层线性回归及预测模型提示该变异与女性口服贝普地尔后QTc延长显著相关，具潜在生物标志物价值，但男性未观察到显著关联。

FSHR基因编码卵泡刺激素受体，该激素由垂体前叶促性腺细胞分泌，在哺乳动物生殖中起核心作用。女性中FSH促进卵泡发育，男性中调节支持细胞功能。绝经后循环FSH水平显著升高，FSHR表达于卵巢外组织（骨、脂肪组织、血管平滑肌），高FSH水平与绝经后女性心血管风险增加相关。近期研究还显示FSH促进绝经女性心房纤维化，故FSHR遗传变异可能通过调控激素信号间接影响心血管生理。尖端扭转型室性心动过速患者存在FSH降低、雌二醇升高的激素特征，绝经后女性雌二醇相对优势可能通过快速延迟整流钾电流（I_Kr）抑制促进QTc延长。复极储备由缓慢延迟整流钾电流及其他K⁺、Ca²⁺电流维持，雌二醇不仅抑制I_Kr，还增加缓慢延迟整流钾电流（I_Ks）并降低L型钙电流，但此调节不足以完全抵消I<sub>Kr抑制致QT延长。孕酮通过增强IKs及抑制L型钙电流对QT延长具保护作用，绝经后孕酮消失而雌二醇经脂肪组织外腺腺持续低水平产生，此激素环境改变可能影响复极储备。研究表明IKr持续抑制导致缓慢延迟整流钾电流代偿性增加，由miR-133a及miR-133b介导；miR-133b参与FSH依赖性雌二醇产生，提示FSH-FSHR轴与miRNA介导的激素调控可能影响心肌离子通道表达，为FSHR内含子变异与女性QT延长关联提供支持。

虽然rs12622919功能影响尚不明确，可能通过调控FSH信号参与QTc延长。基于GTEx数据库，该变异与左心室PPP1R21表达降低显著相关。PPP1R21具RNA结合活性，定位于早期内涵体，双等位基因功能缺失变异导致神经发育综合征伴内涵体功能障碍。其在发育中表达于小鼠皮层，参与内涵体分选成熟及膜蛋白（包括离子通道）循环降解。PPP1R21功能受损可能破坏离子通道循环，导致动作电位时程改变及QTc延长，尽管直接证据尚缺，但本研究提示影响PPP1R21表达的遗传变异可能影响药物诱导QTc延长易感性。

尽管rs12622919功能待阐明，线性回归及AUC性能评估显示其作为预测生物标志物的潜力。本研究为日本人群贝普地尔相关QTc延长调查提供有用框架，rs12622919或为女性患者QTc延长的预测标记。

研究局限性包括：队列规模有限影响统计效能，发现队列仅443例、复制队列78例，低频率变异效应估计不稳定；未对3种遗传模型进行Bonferroni校正，多重检验可能增加假阳性，但rs12622919具稳健全基因组显著性；QTc延长定义阈值的选择可能影响结果；基线QTc测量时胺碘酮等药物停药洗脱期可能不足，RFCA本身可能影响治疗后QTc；预测模型候选SNP自完整数据集识别，内部训练-测试分割无法完全避免过拟合；研究仅限日本人群，推广性不确定；GTEx数据库以欧洲裔为主、亚洲代表不足，需谨慎解读。

发表在《Journal of the American Heart Association》的此项研究通过GWAS发现rs12622919为贝普地尔治疗期间QTc延长的相关遗传变异。探索性性别分层分析提示女性潜在关联，但无显著基因型-性别交互作用，该发现应谨慎解读。这些结果表明遗传分析可能有助于阵发性心房颤动患者接受贝普地尔治疗时QTc延长的风险分层。鉴于主要统计证据基于荟萃分析，需在更大队列中进行独立验证以确认该关联的稳健性及临床实用性。</sub>