摘要：
【目的】通过分析1998—2021年间我国人感染H9N2亚型禽流感病例的发病时间、所在省份、年龄和性别等信息,明确H9N2亚型禽流感病毒的流行病学特征;通过分析人源H9N2亚型禽流感病毒的基因特征,阐明人源H9N2亚型禽流感病毒的遗传演化规律;为H9N2亚型禽流感病毒跨种间传播的预警和防控提供数据支撑。【方法】基于流感基因数据库、病例报道和文献资料,获得1998—2021年我国人感染H9N2亚型禽流感病毒的病例信息和毒株序列数据。从时间、空间、性别和年龄的分布对感染病例进行分析,明确人源H9N2亚型禽流感病毒感染的流行病学特征。通过DNASTAR中的MegAlign软件对人源H9N2病毒的各基因片段的核苷酸序列进行同源性分析,利用MEGA7.0软件构建系统进化树和分析病毒蛋白关键位点,揭示遗传演化趋势和病毒蛋白关键氨基酸位点的变异情况。通过GISAID网站下载2019—2021年间我国H9N2亚型禽流感病毒的核苷酸序列,利用mafft比对后在MEGA7.0中查看人源与禽源H9N2病毒关键氨基酸位点的突变差异,揭示当前人源和禽源H9N2病毒可能引起的风险。【结果】1998—2021年我国人感染H9N2亚型禽流感病毒病例共71例,从空间分布分析,病例分布于16个省市,其中91.55%的病例来自于南方12个省市;从时间分布分析,2013年以后,我国报道的感染病例呈增长趋势,2013—2021年累计感染病例数占总病例数的61.97%;从性别和年龄分布分析,男、女性别比为1:1.68,感染病例主要见于幼儿和少儿,占总病例的74.14%。对人源H9N2病毒进行基因组比对分析,发现这些病毒均属于欧亚分支,但是这些病毒各基因片段的核苷酸序列同源性差异较大,HA、NA、PB2、PB1、PA、NP、M和NS的同源性分别为75.3%—100%、80.1%—100%、78.7%—100%、82.5%—100%、72.6%—100%、74.1%—100%、65.5%—100%、82.0%—100%;22株具有完整基因片段的病毒分为8个基因型,2003、2008和2013年的基因型与1999年的基因型有明显差异。1998—2021年共有42株人源H9N2病毒株上传HA序列,其中有38株病毒的HA蛋白发生Q226L的突变;共有30株人源H9N2病毒株上传PB2序列,其中9株病毒的PB2蛋白发生E627V突变,1株病毒的PB2蛋白发生E627K突变;1株病毒的PB2蛋白的701位点发生D701N突变,共有31株人源H9N2病毒株上传NS与M序列,NS1蛋白的42位点均为S,M1蛋白的30和215位点的氨基酸分别为D和A。2019-2021年人源H9N2病毒的HA蛋白183与190位点、NS1蛋白42位点均发生突变,人源与禽源H9N2病毒的PB2蛋白701位点均未发生突变。【结论】自2013年以来,我国人感染H9N2亚型禽流感病例数量呈增长趋势,且具有显著的地域、年龄和性别差异。1998年至今,人源H9N2病毒的基因同源性差异较大,不同分支间病毒基因重排频繁,形成了复杂的基因型,提示H9N2亚型禽流感病毒在不断地进化。人源H9N2病毒的关键氨基酸位点出现突变,且在2019—2021年人源比禽源H9N2病毒的关键位点突变率高,提示H9N2亚型禽流感病毒的跨种感染人的潜力逐渐增强。该结果丰富了对人源H9N2病毒认知,为H9N2亚型禽流感病毒防控提供参考。
