摘要

目的调查云南省4个艾滋病主要流行地区人类免疫缺陷病毒1型(human immunodeficiency virus type 1, HIV-1)毒株gag基因序列的变异特征。方法利用完全随机抽样法抽取2019年1月至2020年12月云南省昆明市、德宏傣族景颇族自治州、红河哈尼族彝族自治州和临沧市抗病毒治疗定点机构进行随访治疗的480例人类免疫缺陷病毒(human immunodeficiency virus, HIV)感染/艾滋病患者, 收集其流行病学信息, 采集血浆标本后送至云南省传染病医院进行分析, 采用巢式聚合酶链反应扩增HIV-1gag基因, 并进行基因分型, 采用MEGA 6.06软件构建系统进化树, 采用VESPA在线分析工具分析特征性氨基酸, 采用Distance程序计算gag, 以及gag蛋白不同区段的基质蛋白p17、衣壳蛋白p24、核衣壳蛋白p7及p6蛋白的基因距离。采用SNAP程序分析同义突变频率与非同义突变频率比值(Ks/Ka值)。多组间统计学比较采用方差分析。结果 404例患者成功获得gag基因序列, 其中以男性居多(250例, 61.9%), 年龄为40~59岁者占59.7%(241/404), 传播途径主要为异性性传播(61.4%, 248/404)。主要流行的HIV-1亚型依次为流行重组型(circulating recombinant form, CRF)08BC 155例(38.4%), CRF01AE 74例(18.3%), 独特重组型(unique rebombinant form, URF)52例(12.9%)、CRF07BC 39例(9.7%), C亚型34例(8.4%), 其他亚型28例(6.9%), B亚型22例(5.4%)。构建系统进化树发现, CRF08BC亚型形成了2个主要传播簇, 2个传播簇间共有8个位点的特征性氨基酸分布存在显著差异, 分别为第79、93、121、122、151、363、395和396位点。CRF01AE、CRF07BC和CRF08BC的gag基因距离分别为0.090±0.004、0.088±0.004和0.078±0.002, 差异有统计学意义(F=118.33, P<0.001)。3种主要亚型中, gag区段的Ks/Ka值分别为4.003±1.309、4.141±0.860和4.514±1.215, 差异有统计学意义(F=1.35, P<0.001);p17区段的Ks/Ka值分别为2.590±0.186、2.831±0.496和2.936±0.475, 差异有统计学意义(F=1.59, P<0.001);p24区段的Ks/Ka值分别为12.579±1.116、10.185±0.494和8.522±0.595, 差异有统计学意义(F=1.61, P<0.001);p7区段的Ks/Ka值分别为10.850±0.711、9.717±0.932和8.522±0.026, 差异有统计学意义(F=4.24, P<0.001);p6区段的Ks/Ka值分别为3.122±0.134、3.040±1.498和4.841±0.353, 差异有统计学意义(F=10.68, P<0.001)。结论云南省4个地区中主要流行的亚型为CRF08BC亚型, CRF08BC的不同传播簇形成了不同的氨基酸组成模式, 不同亚型gag基因不同区段的变异程度不同, 应加强对HIV变异毒株的监测, 控制其流行蔓延。