国内少数民族VKORC1 基因多态性研究的荟萃分析

普健，薛媚，代松源，陈浙涵，钱孟佼

1.首都儿科研究所附属儿童医院普通（新生儿）外科/北京协和医学院研究生院，北京 100020；2.红河哈尼族彝族自治州第三人民医院甲状腺外科，云南个旧 661199；3.昆明医科大学第一附属医院心内科，云南昆明 650032；4.浙江中医药大学研究生院第二临床医学院，浙江杭州 310053；5.云南省滇南中心医院（红河哈尼族彝族自治州第一人民医院）心胸外科，云南个旧 661199

华法林（Warfarin）作为临床常用抗凝药物，广泛运用于血栓性疾病的治疗和预防[1]。因为其狭窄的治疗窗及个体差异，用药者会有出血或血栓形成的风险。大量的因素能影响华法林维持剂量[2]，其中华法林相关基因：维生素K 环氧化物还原酶复合体亚单位1（VKORC1）的影响较为突出。VKORC1基因可解释近46.9% 的中国人群个体差异[3]。Wang L 等[4]研究显示VKORC1 基因可解释40%左右华法林个体剂量差异。

国外研究发现不同人种VKORC1 基因突变频率的差异能影响华法林使用剂量[5]，国内少数民族群体VKORC1 基因的研究开展较少。Perera MA等[6]发现：为了减少混合基因带来的误差，很多研究把少数民族基因排除在外，但少数民族基因的数据不可或缺。所以中国少数民族是否具有不同于汉族人群的基因多态性？本文就此问题对近年来国内少数民族及汉族华法林VKORC1 基因相关研究进行荟萃分析。

本研究收集自2010 年1 月—2021 年12 月发表的中国少数民族VKORC1 基因研究的文章，仅检索到哈萨克族，维吾尔族，蒙古族，白族和藏族，其他少数民族的研究暂未发现有相关研究报道。本文严格按照纳入、排除标准筛选,在排除重复及不合格研究后,最终纳入9 篇文献。

1 材料与方法

1.1 文献检索

检索PubMed、Cochrane、中国知网（CNKI）、万方和维普数据库,检索年限为2010 年1 月—2021 年12月。中文检索词：VKORC1 基因、华法林、民族、中国；英文检索词:warfarin, vkrc1, gene, minority,China。

1.2 纳入与排除标准

纳入标准：VKORC1 基因研究；中国少数民族人群；中国汉族人群；华法林相关研究。排除标准：缺乏VKORC1 基因数据的文献；非中国地区研究的文献；文献时间较早。

1.3 统计方法

应用SPSS 25.0 统计学软件进行分析，计数资料以[n(%)]表示，采用χ2检验或Fisher 确切概率法；符合正态分布的计量资料以（±s）表示，采用t检验；差异性分析：采用卡方检验时，若P＜0.0125，提示两组基因分布及华法林剂量存在统计学差异。采用多重比较t检验时，P＜0.0167 表示有统计学意义。

2 结果

2.1 VKORC1 多态性的分布数据

近年来国内各少数民族VKORC1 1639G＞A 基因的研究数据汇总，见表1。Zeng WT 等[7]对于云南省大理白族和香格里拉藏族的VKORC1-3673G＞A基因研究的数据，见表2。

表 1 各少数民族及汉族的VKORC1-1639G＞A 基因型分布数据

2.2 各少数民族与汉族VKORC1 基因型比较

VKORC1-1639G＞A 基因型中，汉族与维吾尔族、蒙古族、哈萨克族VKORC1 基因型分布比较，差异有统计学意义（P＜0.05）。而VKORC1-3673G＞A基因型中，汉族与白族基因型分布比较，差异无统计学意义（P＞0.05）。见表3、表4。

表 3 各少数民族与汉族的VKORC1-1639G＞A 基因分布差异性比较

表 4 各少数民族与汉族的VKORC1-3673G＞A 基因分布差异性比较

2.3 华法林剂量比较

对少数民族不用基因型间华法林维持剂量进行统计分析后可得：对于维吾尔族、蒙古族，基因型AA、AG、GG 的华法林日均剂量比较，差异有统计学意义（P＜0.05），按剂量要求GG 型＞AG 型＞AA 型。见表5、表6。

3 讨论

3.1 解释结果

华法林通过抑制VKORC1 引发其抗凝作用，阻碍维生素K 由环氧化物转化成氢醌，从而阻断凝血因子Ⅱ、Ⅶ、Ⅸ、Ⅹ及抗凝蛋白C 和S 的活化达到抗凝的目的。

本研究选取的文献中包含了VKORC1-1639G＞A 和VKORC1-3673G＞A 两个VKORC1 基因的SNP位点。VKORC1-1639G＞A 是VKORC1 基因较常见的位点，有研究表明VKORC1 基因1639G＞A 多态性位点与华法林抗凝治疗关系密切，是种族与民族之间华法林维持剂量出现差异的决定因素。VKORC1-1639GG 被称作非变异或野生基因型，VKORC1-1639GA 是杂合突变基因型，与野生基因型相比，携带杂合突变基因型的患者需要华法林剂量相对较低的；VKORC1-1639AA 是纯合子突变基因型，华法林剂量更低[16-17]。有研究发现携带有VKORC1-1639A 基因的个体对华法林极为敏感[18-20]。鲁晓春等[21]的研究显示中国人群杂合AG型华法林维持剂量较AA 增加40.8%，纯合GG 型较AA 型增加84.2%，AG、GG 型较AA 型华法林剂量要求明显大，按剂量要求GG 型＞AG 型＞AA 型。所以拥有更高VKORC1 基因GG 型频率的维吾尔族、蒙古族、哈萨克族会比汉族人群需要更高的华法林治疗剂量。VKORC1-3673G＞A 是另外一处基因位点，其可以解释20%～40%的个体间华法林剂量差异。

表 2 各少数民族及汉族的VKORC1-3673G＞A 基因型分布数据

表 5 各少数民族VKORC1-1639G＞A 基因华法林日均剂量数据

通过统计学分析可知：汉族VKORC1-1639G＞A基因的基因型分布与维吾尔族、蒙古族、哈萨克族比较，差异有统计学意义（P＜0.001、P=0.007、P＜0.001），可见我国少数民族确实具有不同于汉族人群的VKORC1 基因型及基因频率。对于华法林日均剂量：维吾尔族在VKORC1-1639G＞A 基因AA 型和AG 型，AA 型和GG 型，AG 型和GG 型比较，差异有统计学意义（P=0.001、P＜0.001、P=0.003），蒙古族在VKORC1-1639G＞A 基因AA 型和AG 型，AA 型和GG 型的剂量上比较，差异有统计学意义（P=0.001、P＜0.001），AG 型和GG 型比较，差异无统计学意义（P=0.160）。

3.2 非预期结果

对于VKORC1-3673G＞A 基因中藏族、白族人群与汉族人群不能认为存在统计学差异的结果，分析可能是因为：本研究选取的云南省少数民族与汉族人群分布呈大杂居趋势，使得少数民族人群基因受到到汉族人群基因影响。蒙古族在VKORC1-1639G＞A 基因AG 型和GG 型比较，差异无统计学意义（P＞0.05），分析可能是因为：原文献所取样本量太少，存在选择偏倚。

3.3 研究不足

所收集的研究中采集的样本数量不够大，所研究的少数民族种类也不够丰富。在笔者所构建的统计模型中，若增加其他少数民族的研究样本量可以更进一步深化探讨VKORC1 基因型分布规律。

3.4 结论

综上可见，国内少数民族确实存在其独特的VKORC1 基因多态性,而这些基因多态性又会影响患者华法林的日均剂量，进而对抗凝治疗产生影响。

3.5 研究意义

目前暂未发现有对国内少数民族华法林VKORC1 基因及华法林维持剂量研究数据的综述文献报道，本文有助于进一步推进完善华法林个体化用药。目前已有通过检测基因多态性及临床数据（年龄、体质量等）来构建华法林预测剂量公式来指导华法林的临床应用。当前国际华法林遗传药理协会（International Warfarin Pharmacogenetics Consortium, IWPC）制定剂量计算公式的数据主要来源于欧洲人群，缺乏亚洲人群的基因数据[22]。近年来针对亚洲人群的华法林预测剂量公式的研究也开展了很多，但Asiimwe IG 等[23]的研究观点认为：当前研究集中于主流人种，对于一些少数民族的研究数据缺乏，以此构建的华法林预测公式会存在有偏倚的风险。所以对于少数民族华法林相关基因多态性及华法林剂量的研究有助于完善华法林预测剂量公式，后续研究可以进一步完善各个少数民族与汉族人群的VKORC1 基因多态性数据，构建少数民族相关华法林预测剂量公式，借助华法林预测剂量公式提高患者抗凝治疗的有效性和安全性。

3.6 展望

中国是一个少数民族大国，民族众多，民族地区多而分散，他们拥有较高频率的突变基因型。面对这种情况，完善各地区，特别是少数民族聚集地区的华法林相关基因的检测并建立相适应的剂量预测公式是有必要的，以此能有效降低患者出血及血栓形成的风险。而国内各民族所受环境，饮食，海拔，特殊风俗等因素对个体华法林代谢及华法林相关基因是否会产生影响，这也有待各位学者完善研究。