Citation: Yuanyun Ao, Xiaolu Xie, Xiaogen Dong, Miao Jin, Zhaojun Duan. Genetic Analysis of an Emerging GII.P2–GII.2 Norovirus Associated with a 2016 Outbreak of Acute Gastroenteritis in China[J]. VIROLOGICA SINICA, 2019, 34 (1): 111-114 https://doi.org/10.1007/s12250-019-00084-6
Received: 28 September, 2018; Accepted: 19 December 2018; Published: 6 February 2019
Copyright: © Wuhan Institute of Virology, CAS 2019
Data Availability: All relevant data are within the paper and its Supporting Information files.
Corresponding author: Zhaojun Duan, E-mail: firstname.lastname@example.org.
Yuanyun Ao and Xiaolu Xie have contributed equally to this work.
Electronic supplementary material The online version of this article (https://doi.org/10.1007/s12250-019-00084-6) contains supplementary material, which is available to authorized users
Noroviruses are positive-sense, single-stranded RNA viruses belonging to Caliciviridae and account for more than 50% of all acute gastroenteritis (AGE) outbreaks worldwide and cause an estimated 200, 000 deaths per year among children < 5 years of age, primarily in developing countries (Hall et al.2012; Glass et al.2009).The norovirus genome contains three open reading frames (ORFs). ORF1 encodes nonstructural proteins, including an RNAdependent RNA polymerase (RdRp), while ORF2 and ORF3 encode the major (VP1) and minor (VP2) structural proteins, respectively.Based on the sequences of their VP1 genes, noroviruses are classified into at least seven genogroups (GI-GVII) and more than 30 genotypes (Vinjeé 2015).GI and GII viruses are responsible for most human infections.
Over the past 2 decades, GII.4 viruses have been the most common genotype to cause norovirus outbreaks worldwide, and new GII.4 variants have emerged approximately every 2-3 years through accumulations of mutations and intra-genotype recombination (Lindesmith et al. 2012;Eden et al.2013).During the 2014-2015 winter season, a new GII.17 variant emerged and became the predominant virus in China and parts of Asia, through changes in blockade antibody epitopes and the possible acquisition of a mutated polymerase (Chan et al.2015). However, during 2016-2017, a GII.P16-GII.2 strain reemerged as the major cause of norovirus outbreaks in China and other countries (Ao et al.2017; Bidalot et al.2017; Thongprachum et al.2017; Niendorf et al.2017; Cannon et al.2017), possibly driven by mutations in its nonstructural proteins, particularly the polymerases (Tohma et al. 2017;Ao et al.2018).This necessitates extensive epidemiological studies on GII.2 noroviruses causing both outbreaks and sporadic cases.
In this study, a total of 190 fecal samples associated with 37 norovirus outbreaks from December 2016 to July 2017 in the Fengtai District, Beijing City, were analyzed, using the process developed by the US CDC (Cannon et al. 2017).Thirty-six of 37 outbreaks were caused by the reemerging GII.P16-GII.2 norovirus, while three samples from one outbreak contained the GⅡ.P2-GII.2 norovirus. The analysis of the partially genotyped VP1 sequences (275 nt) of three GII.P2-GII.2 strains (BJFTJYX, BJFTZTX, and BJFTNJY) demonstrated 100% nt identity to the sequence of the reemerging GII.P16-GII.2 norovirus predominant in 2016-2017(Ao et al.2018), but the sequences were distinct from previous GII.P2-GII.2 strains, suggesting the presence of a new GII.P2-GII.2 strain.Therefore, we investigated the evolutionary and genetic characteristics of this strain.The capsid-encoding genes from three GII.P2-GII.2 strains were amplified from three samples using nested RT-PCR, yielding an approximately 2.5-kb amplicon, as described previously (Ao et al. 2018).The complete genome of one representative GII.P2- GII.2 strain BJFTJYX was amplified and sequenced.All sequences obtained in this study have been deposited in GenBank (Accession Numbers MH158635 and MH671553-MH671554).
The time-scale evolutionary trees of 148 GII.2 complete VP1 sequences and 79 partial GII.P2 RdRp sequences were estimated using the strict clock model, GTR+G/HKY+G substitution model, and Bayesian skyline coalescent model implemented in BEAST software (version 1.8.2).Markov chain Monte Carlo (MCMC) sample chains were run for 4×108 and 2.5×108 steps for the VP1 and RdRp genes, respectively, sampling both every 10, 000 generations.The VP1 gene phylogenetic tree could be grouped into three major clusters, according to their polymerase genotypes and sequence divergences, one each for the reemerging GII.P16-GII.2(2016-2017), GII.P16-GII.2(2008-2015), and GII.P2-GII.2(2004-2010) strains.BJFTJYX, BJFTZTX, and BJFTNJY formed peripheral branches in the cluster for the re-emerging GII.P16-GII.2 (2016-2017), which were closely related to the lineage of the 2011-2012 GII.P16-GII.2 strains, including the HS255/USA/2011(KJ407074.2) and Miyagi1/Japan/2012 (LC145787.1) strains.However, they were separated from the previous cluster of the GII.P2-GII.2 (2004-2010) strains (Fig. 1A).The divergence time analysis showed that the most recent common ancestor (MRCA) of the emerging GII.P2-GII.2 strain appeared in 2011-2012, whereas that of the GII.P16-GII.2 strains reemerging in 2016-2017 appeared in 2013(Fig. 1A), indicating that the appearance of the emerging GII.P2- GII.2 strain was earlier than that of the reemerging GII.P16-GII.2 viruses predominant in 2016-2017.These data suggest that this virus did not originate from the previous GII.P2-GII.2 strains but probably resulted from a novel recombination event.The RdRp gene phylogenetic tree could be divided into four lineages based on their evolutionary distances and detected time (Fig. 1B). Lineage 1 can be further split into two sub-lineages.The emerging GII.P2-GII.2 strain BJFTJYX belonged to lineage 1 and sub-lineage 2, which also included three GII.P2-GII.2 strains identified from July to October 2016 in Russia (Fig. 1B).Despite the absence of complete VP1 nt sequences, partial VP1 sequences (600 nucleotides) of the GII.P2-GII.2 strains from Russia in the sub-lineage 2 of lineage 1 were also found to be most related to those of BJFTJYX, BJFTZTX, and BJFTNJY, also belonging to the emerging GII.P2-GII.2 norovirus.This suggests that the emerging GII.P2-GII.2 strain in China may have been introduced from the strains in Russia and that the newly emerged viruses had spread across the countries.The strains in the sub-lineage 2 of lineage 1 were revealed to be relatively closely related to the 2009-2011 GII.P2- GII.2 strains, including the Rotterdam-5-709/NL/2011 (KX446510.1) and Rotterdam-5-0/NL/2009 (KX446499.1) strains.The MRCA of the emerging GII.P2-GII.2 strain appeared in 2012, which is consistent with the results of the analysis of the GII.2 VP1 genes (Fig. 1B).Hence, we conclude that the GII.P2-GII.2 strain emerging in 2016 is a novel recombinant virus that was generated from the 2011-2012 GII.P16-GII.2 and 2011 GII.P2-GII.2 strains in approximately 2011-2012.
The complete genome sequences of the reemerging GII.P16-GII.2 strain BJSMQ (KY421122.1) and GII.P2- GII.2 strain Hokkaido-16(LC209465.1) were used as reference sequences to further predict the recombination nucleotide breakpoint of the emerging GII.P2-GII.2 strain BJFTJYX by the Simplot software v.3.5.1.The SimPlot analysis showed that BJFTJYX shared the highest nucleotide sequence identity with the ORF1 of the GII.P2- GII.2 strain HenrytonSP17, and with the ORF2 and ORF3 of the reemerging GII.P16-GII.2 strain BJSMQ (Supplementary Figure S1).The recombination breakpoint was predicted to be at nucleotide 5177, near the ORF1/ORF2 overlap region (Supplementary Figure S1).
Pairwise comparison of all GII.2 VP1 sequences over 4 decades showed that the GII.2 prototype Snow Mountain Virus (SMV)(AY134748.1) showed an aa divergence of ＜3%, although the root-to-tip divergence plot of the VP1 nt sequences showed a strong clock-like evolution with a coefficient of determination (R2) value of 0.96.The analysis of VP1 sequences revealed that the emerging GII.P2- GII.2 strains showed maximum nt identities of approximately 99% to the 2016-2017 reemerging GII.P16-GII.2 strain CQ25(KY421149.1), while they showed the highest aa identity to the GII.P2-GII.2 strain OC08079 (BAL60765.1), with a divergence of only 1 aa.The alignment of GII.2 VP1 sequences also showed that although the residues adjacent to the HBGA-binding site I appeared to change over time compared to SMV, the major residues of three sites constituting the HBGA-binding interface remained conserved in the emerging GII.P2-GII.2 strain.These data indicate that GII.2 noroviruses have remained very stable over 4 decades, supporting the model presented by Parra et al.(2017).RdRp sequence analysis revealed that the RdRp proteins from these GII.P2 strains showed an aa divergence of less than 1%, although the root-to-tip divergence plot of the GII.P2 partial RdRp nt sequences revealed a moderate clock-like evolution over time (R2=0.60).The RdRp gene of the emerging GII.P2- GII.2 strain BJFTJYX shared 93.5%-95% nt and 99.0%- 99.6% aa identities with those of the previous GII.P2-GII.2 strains.The sequence analyses indicate that the RdRp proteins from GII.P2 strains are highly stable at the aa level.
Our analyses suggest that the GII.P2-GII.2 strain that emerged in 2016 harbors a capsid sequence identical to that in the reemerging GII.P16-GII.2 strain.However, the reemerging GII.P16-GII.2 strain, but not the emerging GII.P2-GII.2 strain, became predominant in 2016-2017 in Fengtai, implying that factors other than capsid proteins contribute to the reemerging GII.P16-GII.2 virus epidemic (Tohma et al.2017; Ao et al.2018).Epidemiological data also found that the GII.P2-GII.2 strains were identified during 1989-2010 and were gradually replaced by the GII.P16-GII.2 strains, which started to circulate in 2008 (Mizukoshi et al.2017).Thus, a possible explanation for the predominance of the reemerging GII.P16-GII.2 viruses is that their GII.P16 polymerase could have a bigger impact on fitness than the GII.P2 polymerase.Enhanced global surveillance of GII.2 genotype norovirus epidemiology is necessary to further understand its evolution and significance in public health.
This work was supported by the Special National Project on Research and Development of Key Biosafety Technologies (2016YFC120 1900) and the National Natural Science Foundation of China (31500139).
Animal and Human Rights Statement
This study was approved by the Institutional Review Board of the China CDC for human subject protection.
- . Ao Y, Wang J, Ling H, He Y, Dong X, Wang X, Peng J, Zhang H, Jin M, Duan Z. 2017. Norovirus GII. P16/GII.2-associated gastroenteritis, China, 23: 1172-1175.
- . Ao Y, Cong X, Jin M, Sun X, Wei X, Wang J, Zhang Q, Song J, Yu J, Cui J, Qi J, Tan M, Duan Z. 2018. Genetic analysis of reemerging GII. P16-GII.2 noroviruses in 2016-2017 in China. J Infect Dis, 218: 133-143.
- . Bidalot M, Thery L, Kaplon J, De Rougemont A, Ambert-Balay K. 2016. Emergence of new recombinant noroviruses GII.2, France, winter 2016 to 2017. Euro Surveill, 22.pii: 30508.
- . Cannon JL, Barclay L, Collins NR, Wikswo ME, Castro CJ, Magana LC, Gregoricus N, Marine RL, Chhabra P, Vinje J. 2017. Genetic and epidemiologic trends of norovirus outbreaks in the United States from 2013 to 2016 demonstrated emergence of novel GII. 4 recombinant viruses. J Clin Microbiol, 55: 2208-2221.
- . Chan MC, Lee N, Hung TN, Kwok K, Cheung K, Tin EK, Lai RW, Nelson EA, Leung TF, Chan PK. 2015. Rapid emergence and predominance of a broadly recognizing and fast-evolving norovirus GII. 17 variant in late 2014. Nat Commun, 6: 10061.
- . Eden JS, Tanaka MM, Boni MF, Rawlinson WD, White PA. 2013. Recombination within the pandemic norovirus GII. 4 lineage. J Virol, 87: 6270-6282.
- . Glass RI, Parashar UD, Estes MK. 2009. Norovirus gastroenteritis. N Engl J Med, 361: 1776-1785.
- . Hall AJ, Eisenbart VG, Etingue AL, Gould LH, Lopman BA, Parashar UD. 2012. Epidemiology of foodborne norovirus outbreaks, United States, 2001-2008. Emerg Infect Dis, 18: 1566-1573.
- . Lindesmith LC, Beltramello M, Donaldson EF, Corti D, Swanstrom J, Debbink K, Lanzavecchia A, Baric RS. 2012. Immunogenetic mechanisms driving norovirus GII.4 antigenic variation. PLoS Pathog, 8: e1002705.
- . Mizukoshi F, Nagasawa K, Doan YH, Haga K, Yoshizumi S, Ueki Y, Shinohara M, Ishikawa M, Sakon N, Shigemoto N, Okamoto-Nakagawa R, Ochi A, Murakami K, Ryo A, Suzuki Y, Katayama K, Kimura H. 2017. Molecular evolution of the RNA-dependent RNA polymerase and capsid genes of human norovirus genotype GII. 2 in Japan during 2004-2015. Front Microbiol, 8: 705.
- . Niendorf S, Jacobsen S, Faber M, Eis-Hubinger AM, Hofmann J, Zimmermann O, Hohne M, Bock CT. 2017. Steep rise in norovirus cases and emergence of a new recombinant strain GII.2, Germany, winter 2016. Euro Surveill, 22.pii: 30447.
- . Parra GI, Squires RB, Karangwa CK, Johnson JA, Lepore CJ, Sosnovtsev SV, Green KY. 2017. Static and evolving norovirus genotypes: implications for epidemiology and immunity. PLoS Pathog, 13: e1006136.
- . Thongprachum A, Okitsu S, Khamrin P, Maneekarn N, Hayakawa S, Ushijima H. 2017. Emergence of norovirus GII. 2 and its novel recombination during the gastroenteritis outbreak in Japanese children in mid-2016. Infect Genet Evol, 51: 86-88.
- . Tohma K, Lepore CJ, Ford-Siltz LA, Parra GI. 2017. Phylogenetic analyses suggest that factors other than the capsid protein play a role in the epidemic potential of GII. 2 norovirus. mSphere, 2.pii: e00187-17.
- . Vinjé J. 2015. Advances in laboratory methods for detection and typing of norovirus. J Clin Microbio, 53: 373-381.