中国农业科学 ›› 2022, Vol. 55 ›› Issue (20): 4075-4090.doi: 10.3864/j.issn.0578-1752.2022.20.017

• 畜牧·兽医 • 上一篇    

1998-2021年我国人感染H9N2亚型禽流感病毒的遗传演化规律

王艳文(),王梦静,张虹,高鑫鑫,郭晶(),李旭勇()   

  1. 聊城大学农学院,山东聊城 252000
  • 收稿日期:2021-08-06 接受日期:2022-04-07 出版日期:2022-10-16 发布日期:2022-10-24
  • 通讯作者: 郭晶,李旭勇
  • 作者简介:王艳文,E-mail: 15290068938@163.com
  • 基金资助:
    山东省自然科学基金(ZR2021MC087);兽医生物技术国家重点实验室开放基金(SKLVBF201906);兽医生物技术国家重点实验室开放基金(SKLVBF202005)

Evolution of Human H9N2 Avian Influenza Virus in China from 1998 to 2021

WANG YanWen(),WANG MengJing,ZHANG Hong,GAO XinXin,GUO Jing(),LI XuYong()   

  1. College of Agronomy, Liaocheng University, Liaocheng 252000, Shandong
  • Received:2021-08-06 Accepted:2022-04-07 Online:2022-10-16 Published:2022-10-24
  • Contact: Jing GUO,XuYong LI

摘要:

【目的】 通过分析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亚型禽流感病毒防控提供参考。

关键词: 禽流感病毒, H9N2亚型, 流行病学, 遗传进化, 跨种间传播

Abstract:

【Objective】 The epidemiological characteristics of H9N2 avian influenza virus were clarified by analyzing the onset time, province, age, gender and other information of human cases infected with H9N2 avian influenza from 1998 to 2021 in China. By analyzing the genetic characteristics of human H9N2 subtype avian influenza virus, the genetic evolution law of human H9N2 subtype avian influenza virus was elucidated. This study could provide data support for the early warning and control of H9N2 subtype avian influenza virus cross-species transmission. 【Method】 According to the case report, literature data and gene bank, the case information and strain sequence data of human infection with H9N2 subtype avian influenza virus in China from 1998 to 2021 were obtained. The infection cases were analyzed from the time distribution, space distribution, gender and age distribution, and the epidemiological characteristics of human-derived H9N2 subtype avian influenza virus infection were found. The MegAlign software in DNASTAR were used to analyze Nucleotide sequence of the homology of various gene fragments of human H9N2 isolates, a phylogenetic tree was constructed by MEGA7.0 software, and the key sites of viral proteins were analyzed. The genetic evolution trend of viral proteins and the variation of key amino acid sites were clarified. The 2019-2021 H9N2 subtype avian influenza virus nucleotide sequence in China was downloaded from the GISAID website. After mafft alignment, the key amino acid site mutation differences between human and avian H9N2 viruses were examined in MEGA7.0, the potential risks brought by current human and avian H9N2 viruses were revealed. 【Result】 From1998 to 2021, a total of 71 cases of human infection with the H9N2 subtype avian influenza virus occurred in China. From the spatial distribution analysis, the cases were distributed in 16 provinces and cities, of which 91.55% of the cases were from 12 provinces and cities in the south. In terms of time, the number of reported infections has been on the rise after 2013, and the cumulative number of infections from 2013 to 2021 accounted for 61.97% of the total number of infections. From the analysis of gender and age distribution, the ratio of male to female was 1:1.68. Human infection with H9N2 virus was mainly seen in infants and young children, accounting for 74.14% of the total number of cases. Genome comparison analysis of human H9N2 viruses found that these viruses belonged to the Eurasian branch, but the nucleotide sequence of homology of these viruses gene fragments was quite different. The nucleotide homology of HA, NA, PB2, PB1, PA, NP, M, and NS were 75.3%-100%, 80.1%-100%, 78.7%-100%, 82.5%-100%, 72.6%-100%, 74.1 -100%, 65.5%-100%, and 82.0%-100%, respectively. According to the genealogy of human isolates, 22 human isolates with complete gene fragments could be divided into 8 genotypes, genotypes in 2003, 2008, and 2013 were significantly different from genotypes in 1999. A total of 42 human H9N2 virus strains uploaded HA sequences from 1998 to 2021, and 38 of them had the mutation of Q226L in the HA protein; a total of 30 human virus strains uploaded PB2 sequence from 1998 to 2021, and 9 strains of which had the mutation of E627V in the PB2 protein, 1 of which had the mutation of E627K in the PB2 protein; the 701 site of the PB2 protein of one virus had the mutation of D701N, a total of 31 virus strains uploaded NS and M sequence from 1998 to 2021, the 42nd positions of the NS1 protein were all S, and the amino acids at positions 30th and 215th of the M1 protein were D and A, respectively. 【Conclusion】 Since 2013, the number of human infections with H9N2 subtype avian influenza reported in China has been on the rise, with significant differences in geographic, age and gender distribution. Since 1998, the genetic homology between human H9N2 isolates was quite different, and the viral gene rearrangement between different branches was frequent, forming a complex genotype. It was suggested that the H9N2 subtype avian influenza virus was constantly evolving. Mutations in key amino acid sites of human H9N2 virus and mutation rate of human viruses was higher than avian H9N2 viruses from 2019 to 2021, suggesting that the potential of H9N2 subtype avian influenza virus to infect humans was gradually increasing. This result enriched the knowledge of human H9N2 virus and provided an important reference for the prevention and control of H9N2 subtype avian influenza virus.

Key words: avian influenza virus, H9N2 subtype, epidemiology, genetic evolution: cross-species transmission