Sequence analysis of pol gene of HIV-1 strains in Kunming, Yunnan Province
-
摘要:
目的 了解云南省昆明市人类免疫缺陷病毒1型(human immunodeficiency virus type 1,HIV-1)的亚型分布特征及其pol结构基因的变异特征。 方法 随机抽取2015—2019年昆明市抗病毒治疗半年以上的样本390例,收集其流行病学信息,血浆提纯HIV-1 RNA,采用巢式聚合酶链式反应扩增HIV-1 pol基因,应用Genotyping Tool、BLAST在线分型工具及MEGA 6.06软件进行基因分型,运用MEGA Distance程序计算基因距离,Entropy程序分析氨基酸多态性,进而对HIV-1流行株的pol基因序列特征进行分析。 结果 成功扩增获得pol基因序列318例,云南省昆明市的HIV-1亚型主要为B亚型占比3.1%,C亚型占比4.4%,CRF01_AE占比26.1%,CRF07_BC占比20.4%,CRF08_BC占比42.8%,其他亚型占比3.1%。各个亚型在性别、年龄及感染途径的构成上差异均有统计学意义(均有P < 0.001)。CRF01_AE、CRF07_BC、CRF08_BC的基因距离分别为(0.018±0.008)、(0.015±0.006)和(0.006±0.004)bp,差异均有统计学意义(均有P < 0.01)。CRF01_AE亚型中静脉药瘾感染者与性接触感染者中共有45个位点的氨基酸差异存在统计学意义(均有P < 0.05),异性性传播与同性性传播中共156个位点有差异。 结论 昆明市HIV-1亚型分布复杂多样,以CRF08_BC为主,传播途径以异性性传播为主,CRF01_AE亚型为昆明地区变异程度较高的流行株,不同感染途径毒株的氨基酸突变位点存在差异,应密切监视其变化。 -
关键词:
- 人类免疫缺陷病毒1型 /
- 亚型分布 /
- 基因距离 /
- 氨基酸变异
Abstract:Objective To investigate the distribution of human immunodeficiency virus type 1(HIV-1) subtypes and the variation of pol gene in Kunming, Yunnan Province. Methods Randomly selected 390 samples with HAART for more than half a year in Kunming from 2015 to 2019, and collected the epidemiological information. Purified the HIV-1 RNA from plasma and amplified the HIV-1 pol gene by nested polymerase chain reaction. Online analysis tools like Genotyping, BLAST and MEGA 6.06 software were used to determine strain subtypes. Distance program of MEGA was used to calculate the gene distance, and Entropy program was used to analyze the amino acid polymorphism, so as to analyze the pol gene sequence characteristics of HIV-1 strains. Results A total of 318 pol gene sequences were successfully amplified. The major subtypes of HIV-1 strains in Kunming were B subtype(3.1%), C subtype(4.4%), CRF01_AE(26.1%), CRF07_BC(20.4%), CRF08_ BC(42.8%) and other subtypes (3.1%). There were statistically significant differences in the distribution of HIV-1 genotypes in sex, age and infection routes (all P < 0.001).The genetic distances of pol gene for CRF01_AE, CRF07_BC and CRF08_BC were (0.018±0.008), (0.015±0.006) and (0.006±0.004) bp (all P < 0.01). In CRF01_ AE, there were 45 different sites in amino acid sequences between intravenous drug addicts and sexual contacts and 156 different sites between hetrosexual contacts and homosexual contacts, and the differences were sitatistically significant (P < 0.05). Conclusion There are complex and diverse HIV-1 subtypes in Kunming, and CRF08_BC is the predominant epidemic subtype. The main route of infection is heterosexual transmission. CRF01_ AE subtype is an epidemic strain with high degree of variation in Kunming. There are differences in amino acid mutation sites of strains from different infection routes, which should be closely monitored. -
表 1 2015—2019年昆明市318例HIV/AIDS患者人口学特征及基因亚型分布情况[n(%)]
Table 1. Demographic characteristics and gene subtype distribution of 318 HIV/AIDS patients in Kunming City from 2015 to 2019 [n(%)]
变量 例数(n=318) B C CRF01_AE CRF07_BC CRF08_BC 其他 χ2值 P值 性别 29.381 < 0.01 男 224(70.4) 9(4.0) 8(3.6) 75(33.5) 42(18.8) 82(36.6) 8(3.6) 女 94(29.6) 1(1.0) 6(6.4) 8(8.5) 23(24.4) 54(57.4) 2(2.1) 年龄(岁) 42.344 < 0.01 < 20 20(6.3) - 2(10.0) 1(5.0) 5(25.0) 12(60.0) - 20~ < 40 153(48.1) 6(3.9) 8(5.2) 57(37.3) 32(20.9) 42(27.5) 8(5.2) 40~ < 60 135(42.5) 4(3.0) 3(2.2) 25(18.5) 26(19.3) 75(55.6) 2(1.5) ≥60 10(3.1) - 1(10.0) - 2(20.0) 7(70.0) - 感染途径 113.846 < 0.01 静脉吸毒 86(27.0) - 3(3.5) 7(8.1) 19(22.1) 57(66.3) - 异性性传播 120(37.3) 4(3.3) 7(5.8) 25(20.8) 25(20.8) 56(46.7) 3(2.5) 同性性传播 62(19.5) 4(6.5) - 38(61.3) 12(19.4) 2(3.2) 6(9.7) 母婴传播 19(6.0) - 2(10.5) 1(5.3) 5(26.3) 11(57.9) - 其他 31(9.7) 2(6.5) 2(6.5) 12(38.7) 4(12.9) 10(32.3) 1(3.2) 合计 318(100.0) 10(3.1) 14(4.4) 83(26.1) 65(20.4) 136(42.8) 10(3.1) 注:其他亚型(例数)依次为A(2)、CRF55_01B(6)、CRF78_cpx(1)以及CRF86_BC(1)。 表 2 2015—2019年昆明市318例HIV/AIDS患者HIV-1 pol基因的三种主要亚型基因距离
Table 2. Gene distances of three main subtypes of HIV-1 pol gene in 318 HIV/AIDS patients in Kunming from 2015 to 2019
HIV-1亚型 例数(n=276) 基因距离(x±s, bp) F值 P值 CRF01_AE 75 0.018±0.008 10 739.856 < 0.001 CRF07_BC 58 0.015±0.006 CRF08_BC 143 0.006±0.004 -
[1] Preston BD, Dougherty JP. Mechanisms of retroviral mutation[J]. Trends Microbiol, 1996, 4(1): 16-21. DOI: 10.1016/0966-842x(96)81500-9. [2] Ho DD, Neumann AU, Perelson AS, et al. Rapid turnover of plasma virions and CD4 lymphocytes in HIV-1 infection[J]. Nature, 1995, 373(6510): 123-126. DOI: 10.1038/373123a0. [3] Wei X, Ghosh SK, Taylor ME, et al. Viral dynamics in human immunodeficiency virus type 1 infection[J]. Nature, 1995, 373(6510): 117-122. DOI: 10.1038/373117a0. [4] Robertson DL, Hahn BH, Sharp PM. Recombination in AIDS viruses[J]. J Mol Evol, 1995, 40(3): 249-259. DOI: 10.1007/BF00163230. [5] Hué S, Clewley JP, Cane PA, et al. HIV-1 pol gene variation is sufficient for reconstruction of transmissions in the era of antiretroviral therapy[J]. AIDS, 2004, 18(5): 719-728. DOI: 10.1097/00002030-200403260-00002. [6] 伦玉华, 罗皓, 方晓君, 等. 东莞市HIV-1流行株pol基因序列特征和亚型分析[J]. 中华疾病控制杂志, 2018, 22(12): 1297-1299. DOI: 10.16462/j.cnki.zhjbkz.2018.12.022.Lun YH, Luo H, Fang XJ, et al. Subtype and sequence analysis of pol gene among HIV-1 strains in Dongguan[J]. Chin J Dis Control Prev, 2018, 22(12): 1297-1299. DOI: 10.16462/j.cnki.zhjbkz.2018.12.022. [7] Santos AF, Soares MA. HIV Genetic Diversity and Drug Resistance[J]. Viruses, 2010, 2(2): 503-531. DOI: 10.3390/v2020503. [8] Njouom R, Pasquier C, Sandres-Sauné K, et al. Assessment of HIV-1 subtyping for Cameroon strains using phylogenetic analysis of pol gene sequences[J]. J Virol Methods, 2003, 110(1): 1-8. DOI: 10.1016/s0166-0934(03)00080-6. [9] Kessler HH, Deuretzbacher D, Stelzl E, et al. Determination of human immunodeficiency virus type1 subtypes by a rapid method useful for the routine diagnostic laboratory[J]. Clin Diagn Lab Immunol, 2001, 8(5): 1018-1020. DOI: 10.1128/CDLI.8.5.1018-1020.2001. [10] Yahi N, Fantini J, Tourres C, et al. Use of drug resistance sequence data for the systematic detection of non-B human immunodeficiency virus type 1 (HIV-1) subtypes: how to create a sentinel site for monitoring the geneticdiversity of HIV-1 at a country scale[J]. J Infect Dis, 2001, 183(9): 1311-1317. DOI: 10.1086/319859. [11] Xu HT, Colby-Germinario SP, Asahchop EL, et al. Effect of mutations at position E138 in HIV-1 reverse transcriptase and their interactions with the M184I mutation on defining patterns of resistance to nonnucleoside reverse transcriptase inhibitors rilpivirine and etravirine[J]. Antimicrob Agents Chemother, 2013, 57(7): 3100-3109. DOI: 10.1128/AAC.00348-13. [12] 霍青青, 李庆海, 刘思宇, 等. 中国部分地区HIV-1 CRF01_AE流行情况[J]. 中国病毒病杂志, 2019, 9(3): 236-241. DOI: 10.16505/j.2095-0136.2019.0017.Huo QQ, Li QH, Liu SY, et al. Prevalence of HIV-1 CRF01_AE strain in partial area of China[J]. Chin J Viral Dis, 2019, 9(3): 236-241. DOI: 10.16505/j.2095-0136.2019.0017. [13] 钟敏, 杨壁珲, 杨绍敏. 云南地区123例艾滋病抗病毒治疗失败病人的HIV-1基因耐药突变位点分析[J]. 中华检验医学杂志, 2013, 36(1): 53-58. DOI: 10.3760/cma.j.issn.1009-9158.2013.01.014.Zhong M, Yang BH, Yang SM. Analysis of genotype resistance mutations sites among 123 antiretroviral-therapy failure in HIV/AIDS patients in Yunnan Province[J]. Chin J Lab Med, 2013, 36(1): 53-58. DOI: 10.3760/cma.j.issn.1009-9158.2013.01.014. [14] Li X, Li W, Zhong P, et al. Nationwide trends in molecular epidemiology of HIV-1 in China[J]. AIDS Res Hum Retroviruses, 2016, 32(9): 851-859. DOI: 10.1089/AID.2016.0029. [15] 杨绍敏, 樊移山, 李惠琴, 等. 云南省德宏和昆明地区高效抗逆转录病毒治疗后HIV-1流行毒株pol区基因变异分析[J]. 中华检验医学杂志, 2011, 34(4): 315-320. DOI: 10.3760/cma.j.issn.1009-9158.2011.04.006.Yang SM, Fan YS, Li HQ, et al. Investigation of pol gene variation of HIV-1 epidemic strains after treatment with HARRT at Dehong prefecture and Kunming in Yunnan Province[J]. Chin J Lab Med, 2011, 34(4): 315-320. DOI: 10.3760/cma.j.issn.1009-9158.2011.04.006. [16] Chen X, Zhou YH, Ye M, et al. Burmese injecting drug users in Yunnan play a pivotal role in the cross-border transmission of HIV-1 in the China-Myanmar border region[J]. Virulence. 2018, 9(1): 1195-1204. DOI: 10.1080/21505594.2018.1496777. [17] Feng Y, Takebe Y, Wei H, et al. Geographic origin and evolutionary history of China's two predominant HIV-1 circulating recombinant forms, CRF07_BC and CRF08_BC[J]. Sci Rep, 2016, 6(1): 19279. DOI: 10.1038/srep19279. [18] Zai JJ, Liu HZ, Lu ZZ, et al. Tracing the transmissiondynamics of HIV-1 CRF55_01B[J]. Sci Rep. 2020, 10(1): 5098. DOI: 10.1038/s41598-020-61870-x. [19] 柳兴凤, 李志坚, 方星, 等. 2013—2017年贵州省艾滋病患者HIV亚型分布[J]. 中华疾病控制杂志, 2019, 23(12): 1523-1526. DOI: 10.16462/j.cnki.zhjbkz.2019.12.017.Liu XF, Li ZJ, Fang X, et al. Distribution of HIV subtypes in HIV/AIDS patients in Guizhou from 2013 to 2017[J]. Chin J Dis Control Prev, 2019, 23(12): 1523-1526. DOI: 10.16462/j.cnki.zhjbkz.2019.12.017. [20] Zheng X, Tian C, Choi KH, et al. Injecting drug use and HIV infection in southwest China[J]. AIDS, 1994, 8(8): 1141-1147. DOI: 10.1097/00002030-199408000-00017. [21] Zhang L, Wang YJ, Wang BX, et al. Prevalence of HIV-1 subtypes among men who have sex with men in China: a systematic review[J]. Int J STD AIDS, 2015, 26(5): 291-305. DOI: 10.1177/0956462414543841. [22] 李敬云. 艾滋病病毒1型CRF01_AE毒株的过去、现在和将来[J]. 中华流行病学杂志, 2016, 37(4): 443-449. DOI: 10.3760/cma.j.issn.0254-6450.2016.04.001.Li JY. Vision on the past, present and future of HIV-1 CRF01_AE[J]. Chinese Journal of Epidemiology, 2016, 37(4): 443-449. DOI: 10.3760/cma.j.issn.0254-6450.2016.04.001. [23] Kandathil AJ, Ramalingam S, Kannangai R, et al. Molecular epidemiology of HIV[J]. Indian J Med Res, 2005, 121(4): 333-344. [24] 邢辉, 梁浩, 万卓越, 等. 中国CRF01-AE亚型人类免疫缺陷病毒毒株的分子流行病学研究[J]. 中华预防医学杂志, 2004, 38(5): 12-16. DOI: 10.3760/j:issn:0253-9624.2004.05.004.Xing H, Liang H, Wan ZY, et al. Distribution of recombinant human immunodeficiency virus type-1 CRF01_AE strains in China and its sequence variations in the env V3-C3 region[J]. Chin J Prev Med, 2004, 38(5): 12-16. DOI: 10.3760/j:issn:0253-9624.2004.05.004. [25] 杨洪, 叶黎, 苏玲, 等. 2011—2015年四川省男男性行为人群HIV-1新发感染率及流行趋势分析[J]. 中华预防医学杂志, 2019, 53(3): 327-329. DOI: 10.3760/cma.j.issn.0253-9624.2019.03.018.Yang H, Ye L, Su L, et al. Ananalysis on incidence of HIV-1epidemics among men who have sex with men in Sichuan Province during 2011-2015 [J]. Chin J Prev Med, 2019, 53(3): 327-329. DOI: 10.3760/cma.j.issn.0253-9624.2019.03.018.