Study on hereditary susceptibility genetic markers to anti-tuberculosis drug induced liver injury in Chinese population
Chenxi Zhou, Mo Li, Cong Huai, Lin He, Shengying Qin,Bio-X Institute, Shanghai Jiaotong University, Shanghai 200030, China通讯作者: 秦胜营,研究员,研究方向:药物基因组学与个体化医学。E-mail:chinsir@sjtu.edu.cn
编委: 张天宇
收稿日期:2019-11-24修回日期:2019-12-26网络出版日期:2020-04-20
基金资助: |
Editorial board:
Received:2019-11-24Revised:2019-12-26Online:2020-04-20
Fund supported: |
作者简介 About authors
周晨希,硕士,专业方向:生物学。E-mail:
摘要
为系统性研究中国人群中抗结核药物引发肝损伤的易感基因标记,本研究以41例抗结核药物引发肝损伤的病人和39例健康对照为研究对象,采用Haloplex捕获测序的方法对其基因组中药物代谢、转运和免疫相关通路的109个基因进行靶向测序。用Plink软件对DNA突变位点与肝损伤的发生进行关联分析,以千人基因组计划东亚人群作为对照组,对显著性位点进行验证,并用SIFT和Polyphen2软件对预测显著关联的位点进行功能预测。结果发现UGT1A4 rs2011404 (χ 2 = 4.6809, P =0.0305)是抗结核药物引发的肝损伤的易感基因标记,且rs2011404突变可能引起UGT1A4蛋白的功能障碍。本研究为临床上对抗结核药物的合理用药提供了有益的参考。
关键词:
Abstract
To systematically study the susceptible genetic markers for liver injury induced by anti-tuberculosis drugs in the Chinese population, 109 genes related to drug metabolism, transport and immunity were captured by Haloplex capture technique from DNA samples of 41 patients with liver injury induced by anti-tuberculosis drugs and 39 healthy controls, and sequenced completely. Association study was conducted using Plink software. To verify the significant candidate SNPs, the χ 2 study was expanded to the control group from the 1000-person Genome Project of the East Asian population. SIFT and Polyphen2 software were used to predict the functional significance of the associated SNPs. Our results identified the UGT1A4 rs2011404 (χ 2 = 4.6809, P = 0.0305) as a susceptible genetic marker for liver injury induced by anti-tuberculosis drugs, and rs2011404 mutation might contribute to UGT1A4 protein dysfunction. This study has provided a potentially useful reference for establishing the precision medicine in rational uses of anti-tuberculosis drugs in the clinic.
Keywords:
PDF (473KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
周晨希, 李沫, 怀聪, 贺林, 秦胜营. 中国人群中抗结核药物引发肝损伤的易感基因标记研究. 遗传[J], 2020, 42(4): 374-379 doi:10.16288/j.yczz.19-271
Chenxi Zhou.
21世纪以来,我国的结核病发病率已降低为本世纪初的一半,但每年新增结核病患者仍有100万之多[1,2]。目前抗结核治疗的标准用药为异烟肼、利福平、吡嗪酰胺联合乙胺丁醇或链霉素等[3]。其中,异烟肼、利福平和吡嗪酰胺均具有潜在的肝毒性[3],其引发的肝脏不良反应发生率约为2.55%~11.9%[4,5],严重时可能会引发肝衰竭导致死亡。抗结核药物所致肝损伤还是导致用药失败、继而诱发结核病耐药的重要原因之一。而抗结核药物引发的肝损伤(anti-tuberculosis drug-induced liver injury, ATLI)具有个体差异明显,难以预测的特点,给病人带来极大的健康危害和经济损失,增加了医疗负担。
现有研究发现遗传因素是抗结核药物引发肝损伤个体间差异的关键因素之一[6]。结核药物导致肝损伤的机制可能有两种:一是药物的异常代谢导致毒性产物的积累,引起胞内的氧化应激反应和线粒体通透性的改变,继而导致肝细胞凋亡或坏死[7];二是药物引起肝脏的免疫调控和炎症反应[8]。因此,本研究系统性地选择了与药物代谢、转运、炎症和免疫相关的109个基因,在抗结核病四联用药(利福平、异烟肼、乙胺丁醇和吡嗪酰胺)引发肝损伤的病人与健康人中进行了关联分析,旨在找出抗结核药物导致肝损伤的易感基因位点,作为潜在的临床预测分子标记。
1 对象与方法
1.1 实验对象
病例组收集自上海肺科医院服用抗结核病四联用药(利福平、异烟肼、乙胺丁醇和吡嗪酰胺)后出现药物性肝损伤(drug-induced liver injury, DILI)的结核病患者。纳入标准:(1)有完整的用药信息和服药前后肝功能信息;(2)按照国际严重药物不良反应联盟(International Serious Adverse Events Consortium, iSAEC)标准评判肝损伤:①丙氨酸氨基转移酶或谷草转氨酶活性高于5倍正常值上限;或②碱性磷酸酶活性高于2倍正常值上限(当没有骨源性原因引起碱性磷酸酶值上升时);或③丙氨酸氨基转移酶活性高于3倍正常值上限并且总胆红素超过2倍正常值上限。(3)因果关系评分量表(roussel uclaf causality assessment method, RUCAM)得分≥7 (极可能或确定为药物引起的肝损伤)[9]。排除标准:同时服用中草药或中成药的患者。病例组最终纳入41例患者,男25例,女16例,年龄19~70岁,平均年龄44.88 (±17.63)岁。对照组为来自上海新华医院的健康体检者39例。其中,男23例,女15例,年龄22~77岁,平均年龄35.66 (±12.83)岁。
本研究经过上海交通大学Bio-X研究院伦理委员会批准(批准号: M2011003),获得所有参与者的同意并签署知情同意书。
1.2 基因组DNA抽提
采用AxyPrep Blood Genomic DNA Miniprep Kit (Axygen,美国)提取受试者基因组DNA。1.3 Haloplex捕获测序
针对109个药物相关代谢基因(表1)的外显子、3ʹUTR和5ʹUTR区域,使用Agilent Technologies的在线设计平台SureDesign (www.agilent.com/genomics/ suredesign) 设计Haloplex捕获探针。使用Haloplex Target Enrichment System Kits (Agilent Technologies, 美国) 酶切基因组DNA,与Haloplex探针杂交。经过扩增,纯化后的杂交片段由2100 Bioanalyzer (Agilent,美国)和Qubit dsDNA HS Assay Kit/Qubit 2.0定量分析仪(Life Technology,美国)测定片段大小及浓度。长度在175~625 bp之间的杂交片段为合格目的片段。使用Illumina MiSeq测序平台对合格片段进行测序。使用Agilent SureCall软件对测序结果进行分析,得到突变SNP位点的信息。Table 1
表1
表1抗结核药物导致肝损伤的109个候选基因
Table 1
基因类型 | 基因家族 | 目标捕获基因 |
---|---|---|
一相药物 代谢基因 | CYP家族基因 | CYP1A1、CYP1A2、CYP1B1、CYP2A13、CYP2A6、CYP2A7、CYP2B6、CYP2C18、CYP2C19 CYP2C8、CYP2C9、CYP2D6、CYP2E1、CYP2F1、CYP2J2、CYP2R1、CYP2S1、CYP2U1、 CYP2W1、CYP3A4、CYP3A43、CYP3A5、CYP3A、CYP4V2、CYP4X1和CYP4Z1 |
超氧化物歧化酶 | SOD1、SOD2和SOD3 | |
含黄素单加氧化酶 | FMO1、FMO2、FMO3、FMO4和FMO5 | |
其他 | TPMT | |
二相药物 代谢基因 | UGT家族 | UGT1A1、UGT1A3、UGT1A4、UGT1A6、UGT1A9、UGT2B15和UGT2B7 |
GST家族 | GSTM1和GSTT1 | |
NAT家族 | NAT1和NAT2 | |
三相药物 转运蛋白 | ATP结合盒式 转运蛋白 | ABCB1、ABCB11、ABCB4、ABCC1、ABCC10、ABCC1、ABCC12、ABCC2、ABCC3、ABCC4、ABCC5和ABCC6 |
其他转运蛋白 | SLCO1B1、POU2F1和POU5F1 | |
调控受体 | 调控受体 | AHR、ARNT、ESR1、ESR2、FOXA1、FOXA2、FOXA3、HNF1A、HNF1B、HNF4A、HNF4G、NR1H2、NR1H3、NR1H4、NR1I、NR1I3、PPARA、PPARD、PPARG、RARA、RARB、RARG和VDR |
固有免疫 | 白细胞介素 | IL10、IL12A、IL12B、IL13、IL18、IL1A、IL1B、IL2、IL4、IL5、IL6、IL7和IL9 |
其他细胞激素 | CCL2、IFNG和TNF | |
适应性免疫 | 人类白细胞抗原 | HLA-A、HLA-B、HLA-C、HLA-DPA1、HLA-DPB1、HLA-DQA1、HLA-DQB1、HLA-DRA和HLA-DRB1 |
新窗口打开|下载CSV
1.4 统计学分析
首先对SNP位点进行质量控制,删除缺失率> 10%,哈迪-温伯格遗传平衡检验(Hardy-Weinberg equilibrium, HWE) P<10-4的SNP位点。针对每个SNP位点病例组和对照组的等位基因频数进行卡方检验,P<0.05为有统计学意义。在验证研究中,将病例与对照组中存在的显著差异位点与其在千人基因组计划东亚人群中的等位基因频数进行卡方检验。统计学分析由PLINK (version1.07,1.5 突变位点功能预测
使用SNP-nexus在线软件(2 结果与分析
2.1 捕获测序结果
在80例样本中,使用Haloplex测序试剂盒和Miseq平台对109个药物代谢相关基因(表1)进行捕获测序。突变的平均测序深度大于100×,检测到突变位点共5130个。质量控制后,位于57个基因上的228个SNP位点最终纳入关联分析。2.2 与抗结核药物导致肝损伤相关的风险SNP/SNV位点
在41例患者和39例健康人样本中对SNP/SNV和抗结核药物引发的肝损伤进行关联分析,结果见表2。其中3个位点显示出显著性差异,分别为UGT1A4 rs2011404 (OR = 0.1939, 95%CI = 0.0740- 0.5083, P = 0.0004), CYP2D6 rs16947 (OR = 6.5620, 95%CI = 1.4140-30.4500, P = 0.0069)和CYP2S1 rs338599 (OR = 0.3944, 95%CI = 0.1646-0.9446, P = 0.0330)。为提高统计准确性,将对照组扩大为千人基因组计划中的东亚人群(n = 504)对显著性位点进行了验证。结果发现,UGT1A4 rs2011404 (χ2= 4.6809, P = 0.0305)依然具有显著性差异,与抗结核药物导致肝损伤显著相关(表3)。Table 2
表2
表2与抗结核药物导致的肝损伤显著相关的SNP/SNV
Table 2
基因 | 编号 | 频数 REF/ALT | 比值比(OR) | 95%置信区间(95%CI) | P值 | |
---|---|---|---|---|---|---|
结核病人 | 健康人 | |||||
UGT1A4 | rs2011404 | 6/74 | 23/55 | 0.1939 | 0.0740~0.5083 | 0.0004 |
CYP2D6 | rs16947 | 64/12 | 70/2 | 6.5620 | 1.4140~30.4500 | 0.0069 |
CYP2S1 | rs338599 | 71/9 | 56/18 | 0.3944 | 0.1646~0.9446 | 0.0330 |
新窗口打开|下载CSV
Table 3
表3
表3验证与抗结核药物导致肝损伤显著相关的SNP/SNV
Table 3
基因 | 编号 | 等位基因 REF/ALT | 等位基因频率 | 卡方(χ2) | P值 | |
---|---|---|---|---|---|---|
结核病人 | 东亚健康人群数据库 | |||||
UGT1A4 | rs2011404 | T/G | 0.0750/0.9250 | 0.0258/0.9742 | 4.6809 | 0.0305 |
CYP2D6 | rs16947 | G/A | 0.8421/0.1579 | 0.8601/0.1399 | 0.0698 | 0.7917 |
CYP2S1 | rs338599 | C/G | 0.8875/0.1125 | 0.7897/0.2103 | 3.7978 | 0.0513 |
新窗口打开|下载CSV
2.3 显著性SNP/SNV突变位点蛋白功能影响预测
用SIFT和PolyPhen2对rs2011404突变的导致的蛋白功能变化进行预测,结果显示rs2011404位置上的T突变为G,可能导致157位的半胱氨酸突变为色氨酸,很有可能为有害的突变(表4)。Table 4
表4
表4显著性SNP/SNV突变位点蛋白功能影响预测
Table 4
基因 | 编号 | 位置 | 等位基因 | 氨基酸改变 | SIFT预测 | PolyPhen2预测 |
---|---|---|---|---|---|---|
UGT1A4 | rs2011404 | 编码区 | T>G | C157W或C157C | deleterious | Probably damaging |
CYP2D6 | rs16947 | 编码区 | G>A | R296C或R274C或R245C | tolerated | Benign |
CYP2S1 | rs338599 | 编码区 | C>G | P74P | tolerated | - |
新窗口打开|下载CSV
用SWISS-MODEL和PyMol对UGT1A4蛋白结构进行同源模拟(图1A),蛋白结构显示rs2011404 (C157W)位于α螺旋上。野生型UGT1A4 157位的半胱氨酸可能会形成二硫键(图1B),而其突变为色氨酸(图1C)后,可能会影响二硫键的形成,进而影响α螺旋的稳定性。
图1
新窗口打开|下载原图ZIP|生成PPT图1突变蛋白结构分析
A:野生型UGT1A4蛋白结构模拟(粉色部分为α螺旋,青色部分为β折叠,蓝色部位为157位半胱氨酸);B:野生型UGT1A4的157位半胱氨酸局部放大示意图:C:rs2011404突变型UGT1A4的157位色氨酸局部放大示意图。
Fig. 1Structure analysis of mutant protein
3 讨论
药物性肝损伤是在药物使用过程中,因药物或其代谢产物产生超敏反应导致的患者肝细胞毒性损伤及肝功能异常,为最常见和严重的药物不良反应(adverse drug reaction, ADR)之一。临床上,药物性肝损伤导致了3%~5%的黄疸症,也是造成急性肝功能衰竭的最主要原因,重者可致死亡[10]。现已知超过1100种药物具有潜在的肝毒性,包括中草药、抗结核药、抗感染药、解热镇痛药,抗肿瘤药等[11]。中国人群中,抗结核药物以异烟肼,利福平和吡嗪酰胺为主引发的肝损伤占所有药物性肝损伤的21.56%,是仅次于中草药的第二大造成药物性肝损伤的原因[12]。非遗传因素(如性别、年龄、饮酒、肝病史、伴随感染和营养状况等)是抗结核药物引发肝损伤的危险因素,但其在临床上对肝损伤发生的预测作用有限[6]。国内外研究发现一相药物代谢酶CYP2E1[13]、二相药物代谢酶NAT2[14]、GSTM1[15]、三相药物转运体ABCB1[16]、SLCO1B1[17]、免疫调节基因HLA- DQA1和HLA-DQB1[18]等与抗结核药物引发的肝损伤显著相关。这些研究验证了遗传因素对抗结核药物引发肝损伤的贡献。现有的中国人群中的研究多针对少数候选基因,可能会遗漏潜在的SNP/SNV位点,缺少全面的研究。本研究利用Haloplex捕获测序的方法,捕获与药物代谢、转运、炎症和免疫相关的一、二和三相代谢酶,调控受体,以及固有免疫和适应性免疫家族的共109个基因上的突变位点,旨在系统性地寻找抗结核药物引发肝损伤的潜在遗传标记。
本研究首次发现rs2011404 (471T>G)与抗结核药物引发的肝损伤显著相关。当对照组人群扩大为千人基因组计划中的东亚人群时,此SNP位点依然呈阳性。rs2011404位于UGT1A4基因区段上,此位点上T突变为C或G。其中,T>C不会造成氨基酸的改变,而T>G造成157位的半胱氨酸(C)突变为色氨酸(W)。后续的SIFT和PolyPhen2预测指出,C157W突变极有可能为有害的突变,影响UGT1A4代谢酶的正常功能。UGT1A4 (UDP glucuronosyltransferase family 1 member A4)属于UGT (UDP glucuronosyltransferase)家族,编码葡萄糖醛酸转移酶,负责催化一系列外源或内源化合物氨基的葡萄糖醛酸化(N-glucuronidation),属于II相药物代谢 酶[19]。Chang等[20]在中国台湾人群中发现UGT1A1* 27,UGT1A1*28会增加抗结核药物引发的肝损伤的风险,而尚未有研究报道UGT1A4与抗结核药物引发的肝损伤的相关性。Gufford 等[21]和Lee 等[22]证明利福平会诱导UGT1A4的表达,继而影响药物的代谢。与之相反,Cao等[23]发现利福平会抑制UGT1A4的活性。目前尚未发现异烟肼,吡嗪酰胺和乙胺丁醇对UGT酶有抑制作用。本研究发现UGT1A4与抗结核药物引发的肝损伤显著相关,而对于UGT1A4在抗结核药物代谢中的作用机制仍需进一步的实验验证。
CYP2D6 rs16947和CYP2S1 rs338599在80例患者的关联研究中呈阳性,而在扩大对照组为千人基因组计划东亚人群后未达到统计学意义。这一结果可能是由于本次研究的样本收集困难,样本量较少,而捕获测序后纳入的位点较多导致的,这可能是本研究的一个局限。
综上所述,本研究通过对109个药物代谢、转运、炎症和免疫相关的基因突变位点的系统性研究,发现UGT1A4 rs2011404作为抗结核药物引发的肝损伤的潜在基因标记,为临床上对抗结核药物的合理用药提供了有益的参考。
参考文献 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 2]
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 2]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]