论文部分内容阅读
目的:分析乙型肝炎病毒(HBV)CD重组体的遗传变异情况;探寻藏族人群HBV感染者血清中表面抗原(HBsAg)和抗表面抗原抗体(HBsAb)双阳性发生率及原因;分析HBV/CD重组体的时空动态传播过程和进化遗传学特征(包括基因序列起源、核苷酸替换率和居群动力学等)。方法:利用多阶段随机抽样的方法,在西藏自治区和青海省海南州进行社区藏族人群血清样本采集,西藏七个地市中共收集HBsAg阳性血清样本852例,青海省收集HBsAg 阳性血清样本411例。按照采样地区藏族人群人口分布情况,选取具有代表性的血清样本进行DNA提取,分别对HBV全长和BCP区采用巢式PCR方法进行扩增、测序和拼接,以获得HBV全序列。对获得的全序列进行重组分析,对不同重组型的基线情况、遗传变异、血清型进行比较;对重组体内不同基因型片段进行系统发育分析。以青藏高原藏族人群HBsAg和HBsAb双阳性血清样本做为观察组,并以1:4配比选择年龄、性别、病毒DNA水平和HBeAg状态同质的HBsAg 阳性/HBsAb阴性血清样本作为对照组,对两组间基因序列S区、PreS区氨基酸突变率、缺失和全基因组核苷酸突变位点分布情况进行比较。使用自NCBI中检索、收集的HBV C2亚基因型、D基因型和CD重组体全长序列,合并本研究中获得的185株全长序列,构建数据库。使用BEAST软件包对CD重组体中的D基因型片段和D1-5亚基因型;CD重组体中的C基因型片段和C2亚基因型分别进行系统发生生物地理学分析。同时,对CD1和CD2全长序列单独进行溯祖分析,并使用贝叶斯skyline重建方法推断种群数量的历史变化。结果:共获得185条HBV全长序列,其基因分型主要为CD重组基因型(181/185,97.84%),包括:132株CD1型,其nt10-799为D基因型,其余部分为C基因型;49株CD2型,其nt10-1499为D基因型,其余部分为C基因型;除CD重组型,还获得了 4株C2亚基因型序列。185株测序样本中,HBeAg阳性组的病毒DNA水平高于HBeAg阴性组(P=0.001);HBeAg阴性组年龄高于HBeAg阳性组(P=0.006)。CD1和CD2的C基因型区域均与C2亚基因型最为接近;CD重组体的D基因型段(nt10-799)与D4亚基因型最为接近,而CD2型nt800-1499序列在系统发育树中与D1-D3更为接近。CD2型重组体主要分布于与印度接壤的山南地区和日喀则市(P=0)。CD重组体在S区氨基酸aa76、aa120和aa129与D基因型存在差异。HBV/CD重组体的主要血清学型为ayw2型(175/181,96.7%)。与CD1型相比,CD2型存在较多的S207N突变。HBV/CD重组体的T1762/A1764双突变发生率低于C基因型,与D基因型接近。HBV/CD病毒株的 G1896A(11.1%vs 25.6%)和A2189C(11.1%vs 50.0%)突变频率低于 C2基因型,G1613A(8.9%vs 10.3%)突变水平与C2基因型相近。青藏高原1263份HBsAg阳性样本的双阳性率为2.14%(27/1263)。在S蛋白全长、N端,MHR(aa100-169),“a”决定簇(aa124-147),茎环结构第一环(aa124-137),第二环(aa139-147)和C端的氨基酸突变率,双阳性组均显著高于对照组;双阳性组中PreS区缺失显著高于对照组(以上结果中,均为P<0.001)。PreSI和PreS2氨基酸突变率在两组间无显著性差异(P=0.19;P=0.89)。PreC/C区C2002T、A2159G、A2189C、C2198A在两组间的分布均存在显著差异。双阳性组和对照组的核苷酸点突变比较,PreS/S区C3189A和T825C、P区C930A在两组间分布存在显著性差异。其中,PreS区C129T(S155L)仅存在于对照组;PreS区C3189A(D103E)和S区T825C(V224A)突变仅存在于双阳性组。HBVD基因型和CD重组体中的D型区域(nt10-799)可溯祖至271年前(95%HPD 138-477)的亚洲南部,其平均进化速率约为1.26×104(95%HPD 5.92×10-5-1.93×10-4)/位点/年。D4亚基因型是CD1和CD2中最为接近的分支,于1850年左右(95%HPD 1743-1934)进入青藏高原;1911年左右(95%HPD 1845-1971)分化为CD1和CD2内相应片段。CD1型较CD2型出现更早,表明先形成D基因型区段为nt10-799的CD1,后续形成D基因型区段为10-1499的CD2。HBV/CD1和HBV/CD2重组体分别形成三个独立的进化群。HBV/CD1重组体自1980年代开始至2000年,其种群大小保持稳定;HBV/CD2重组体自1997年首次报道以来,其种群大小经过短暂上升之后,随即出现了明显的下降趋势。全球C2亚型片段的传播最先起源于亚洲东部,首先传播到中国西北和亚洲北部,随后东亚分离株又进一步播散到美洲。结论:HBV/CD重组体是高原地区流行的主要基因型。HBeAg阴性组患者体内重组体DNA水平更高,病毒复制速度较快,病情处于进展期。年龄增长是导致HBeAg在感染HBV/CD的藏族人群血清中发生转阴的主要原因。HBV/CD1重组来源为D4和C2亚基因型;HBV/CD2可能经过重组来源为D1-3、D4和C2亚基因型的多次重组。基于各基因型的地理分布,形成CD重组体的主要重组来源可能为中国(C2亚基因型)和印度(D4亚基因型)。与D基因型相比,HBV/CD重组体的主要血清型无明显变化。HBV/CD重组体感染者均在高海拔地区生活,造成这种现象的可能原因是高原居民具有相近的遗传背景特征或长期的高海拔生活导致其具有相同或相似的生活习惯。青藏高原HBV表面抗原阳性人群中HBsAg和HBsAb双阳性的发生率接近大多数不同基因型的既往研究结果。与对照组相比,双阳性组的氨基酸突变率在S区全长及其内的各个区段均显著增加。S区高突并不是血清双阳性病的充分必要条件,并不能完全解释HBsAg/HBsAb共存现象,PreS区的改变可能是导致HBsAg/HBsAb双阳性的因素之一。S区T825C(V224A)、PreS区C3189A(D103E)仅存在于双阳性组中,这些氨基酸变异导致蛋白结构发生变化,可能是双阳性发生的重要条件。HBVD基因型最可能的起源地是南亚地区。D基因型原型株在18世纪中叶逐步从印度传播到中亚进而传播到欧洲、地中海地区和非洲。在此过程中,19世纪中叶先后分化出D1、D3和D2亚基因型。同时19世纪中叶自印度向东北方向传播,逐步形成D4亚基因型并在1910年左右进入青藏高原及周边地区与C2型发生重组,这与英军两次入侵西藏时间相吻合。欧洲主要HBV D亚型的传播可追溯到二十世纪早期,包括第一次世界大战特别是第二次世界大战在内的一段时期。此过程可能与战争中产生的创伤、输血和疫苗注射,以及有组织的大规模人口迁移导致HBV的快速传播和演化有关。C2基因型向美洲等地的播散则与亚洲在二十世纪八九十年代的欧美移民潮有关,随着大量东亚地区亚裔人口向欧美等发达国家移民导致C2亚基因型的世界范围传播。