- Research article
- Open Access
- Open Peer Review
Changes in norovirus genotype diversity in gastroenteritis outbreaks in Alberta, Canada: 2012–2018
BMC Infectious Diseasesvolume 19, Article number: 177 (2019)
The emergence of norovirus genotype GII.4 variants has been associated with gastroenteritis pandemics worldwide, prompting molecular surveillance for early detection of novel strains. In this study, we aimed to analyze the outbreak activity of norovirus and characterize the norovirus strains circulating in Alberta between July 2012 and February 2018.
Stool samples from gastroenteritis outbreaks in Alberta were tested for norovirus at the Provincial Laboratory for Public Health using a multiplex real time-RT PCR assay. The ORF1 and ORF2-genotypes of norovirus positive samples were assigned based on phylogenetic analyses of partial polymerase and capsid sequences, respectively.
A total of 530 norovirus outbreaks were identified. During July 2012 and June 2017 there was a gradual decrease in the annual number of GII.4 outbreaks, however, outbreak numbers increased from June 2017–February 2018. Four novel strains emerged: GII.17 Kawasaki in July 2014–June 2015, GII.P16/GII.4 Sydney in July 2015–June 2016, GII.P16/GII.2 and GII.P4 New Orleans/GII.4 Sydney in July 2016–June 2017. GII.Pe/GII.4 Sydney was the single predominant strain responsible for the majority (over 50%) of all norovirus outbreaks up to June 2015. Between June 2017 and February 2018, GII.P16/GII.4 Sydney was the leading strain causing 63% of all norovirus outbreaks.
GII.4 stands as the predominant capsid genotype causing a large majority of the norovirus outbreaks in early 2018. An increase in genotype diversity was observed in the last years, characterized by a high circulation of non-GII.4 strains and GII.4 recombinants.
With advancement in molecular testing, Norovirus (NoV) is recognized as the leading cause of morbidity and mortality from diarrhoeal disease across all ages . Approximately 18% of all acute gastroenteritis around the world is caused by NoV with an estimated economic burden of 60 billion US dollars per year [2, 3]. Norovirus has a single-stranded positive-sense RNA genome of ~ 7.5 Kb encompassing three open reading frames (ORFs). The three ORFs encode a non-structural polyprotein, the major capsid gene (VP1) and the minor capsid gene (VP2).
The NoV genome has high genetic variability at the P2 subdomain of VP1. P2 forms the outermost part of the virion and contains binding sites for histo-blood group antigens (HBGAs), attachment factors for human NoV strains. Norovirus can evade host immune responses by antigenic drift, through a process similar to influenza’s epochal evolution . Recombination is another mechanism of evolution for NoV that often occurs at the ORF1/ORF2 junction, further increasing the genetic diversity of the virus . Seven NoV genogroups (GI to GVII) have been described to date of which GI, GII and GIV can cause gastroenteritis in humans, with GIV being less common . NoV genogroups are sub classified into genotypes; at least 14 ORF1-based and 9 ORF2-based genotypes have been described for GI and 27 ORF1-based and 22 ORF2-based genotypes for GII .
NoV GII.4 has been the most common genotype in circulation worldwide since the mid-1990s. New genetic clusters or variants of NoV GII.4 have emerged periodically and caused pandemics including: US95/96 in 1996, Farmington Hills in 2002, Hunter in 2004, Den Haag in 2006, New Orleans in 2010 and Sydney in 2012 [7,8,9]. No global pandemic GII.4 strain has yet emerged after GII.4 Sydney. A novel cluster of GII.17 named Kawasaki became predominant in Asia during the 2014–2015 winter season  raising concern over a possible new global pandemic NoV genotype. Since its emergence in Asia, GII.17 Kawasaki has also been reported in other continents although with lower prevalence .
It is yet unknown why GII.4 has remained as the predominant genotype for over two decades but several features of this genotype support its enhanced circulation, including higher rates of evolution and a progressive accumulation of mutations that help evade host immune responses [11, 12]. Moreover, pandemic GII.4 variants can broadly bind to a wide set of HBGAs, a feature that favours virus transmission by providing NoV a large pool of individuals genetically susceptible to infection .
The aim of the present study is to describe the outbreak activity of NoV from July 2012 to February 2018 in Alberta, Canada. We provide information on the norovirus genotypes in circulation and their relevance in outbreak settings. This study provides important data for vaccine development and enhances our understanding of norovirus disease burden.
Gastroenteritis outbreak investigations in Alberta were managed by public health officials in collaboration with the Provincial Laboratory for Public Health (ProvLab) . Stool samples collected between July 2012 and February 2018 during outbreak investigations were tested at ProvLab for NoV genogroup I and genogroup II using a real time RT-PCR assay . A NoV- confirmed outbreak was defined as ≥2 epidemiologically linked cases with gastroenteritis and at least one sample tested positive for NoV. Data was analyzed from July to June of the following year considering the winter seasonality of NoV.
Outbreak settings were classified into 6 different groups: 1) community long-term care, hospital long-term care, supportive living, and group homes facilities; 2) hospital acute care; 3) food establishments, catering events, food shops, community functions, conferences and hotels; 4) day care centers; 5) other types of group residences, e.g. camp, dormitory; and 6) other settings including community shelters, community services, household, schools and cruise ships. Groups 1 and 2 were classified as healthcare-related settings, whereas groups 3 to 6 were classified as non-healthcare-related.
For each NoV-positive outbreak, one NoV positive sample was selected for genotyping. Briefly, the nucleic acid extract from each stool sample was subjected to reverse transcription (RT) with random primers and the resulting cDNA was PCR-amplified. Samples collected up to February 2017, were amplified in region C (ORF2) using primer pair G2SKF/G2SKR for genogroup II strains or primer pair G1SKF/G1SKR for genogroup I strains . Samples with emergent or unassigned genotypes based on region C sequence analysis were further genotyped targeting the 3’end of the polymerase gene using primers LV4282-99F  and COG2R . Samples collected from March 2017 onwards, were genotyped using a dual polymerase-capsid genotyping protocol based on a single PCR amplicon obtained with primer pair MON432/G1SKR for genogroup I strains and primer pair MON431/G2SKR for genogroup II strains . All PCR products were subjected to Sanger sequencing and genotypes were assigned using the Norovirus Genotyping tool . A large majority of strains from outbreaks occurring between July 2012 and June 2015 were left uncharacterized at ORF1; retrospective characterization was not attempted based on observations that 69% of norovirus GII outbreaks in Alberta were caused by a single ORF2 genotype, GII.4 Sydney, and reports from North America  and diverse countries from different continents  suggest that GII.Pe/GII.4 Sydney was the major strain in circulation worldwide during that time frame.
Phylogenetic analyses of GII.17, GII.P16/GII.4 Sydney, GII.P4 New Orleans/GII.4 Sydney and GII.P16/GII.2 sequences obtained in our study were performed with MEGA 6.06 . Maximum likelihood trees were constructed using the substitution model producing the lowest Bayesian Information Criterion scores, as calculated by the software. For all trees, branch significance was estimated based on 1000 bootstrap replicates.
Cloning, expression and purification of recombinant P domain proteins
The capsid P domain of two GII.P16/GII.4 Sydney outbreak strains (AB-2016-26 and AB-2016-190) were amplified by RT-PCR using forward primer ACGCGGATCCTCAAGAACTAAACCATTCTCTGTCC and reverse primer ATAAGAATGCGGCCGCTTAGCAAAAGCAATCGCCACGGCAATCGCATACTGCACGTCTACGCCCCGTTCC and cloned into pGEX-4 T-1 vector (GST Gene fusion System, GE Healthcare Life Sciences) between the Bam HI and Not I sites. An RGD4C tag (CDCRGDCFC) was linked to the C terminus of the P domain for P particle formation . The recombinant P domain protein was expressed in E. coli (BL21, DE3) with induction by 0.25 mM isopropyl-β-D-thiogalactopyranoside (IPTG) at room temperature (~ 21 °C) overnight as described elsewhere . Purification of the glutathione S-transferase (GST)-P domain fusion protein was performed using resin of Glutathione Sepharose 4 Fast Flow (GE Healthcare Life Sciences) according to the manufacturer’s instruction. GST was removed from the target proteins by thrombin (GE Healthcare Life Sciences) cleavage either on beads or in solution (phosphate buffer saline, PBS, pH 7.4).
Saliva binding assay of P-domain proteins
The saliva-based binding assays were performed as previously described . Briefly, boiled human saliva with known HBGA phenotypes collected from Cincinnati, OH, United States, were diluted 1000-fold and used to coat 96-well microtiter plates (Dynex Immulon; Dynatech, Franklin, MA). After blocking with 5% non-fat milk in PBS, different concentrations of P-domain protein (15, 7.5, 3.75 ng/μl) were added to the wells. The bound P proteins were detected using a guinea pig anti-NoV antiserum (1:3000), followed by the addition of HRP-conjugated goat anti-guinea pig IgG. The HRP activity was then measured with TMB kit (Kierkegaard & Perry Laboratories, Gaithersburg, MD) and the OD450 values were read with an ELISA spectrum reader (Tecan, Durham, NC).
The proportion of NoV GI and GII outbreaks by settings were compared using the Chi-square exact test. The annual numbers of NoV positive outbreaks occurring between July 2012 and June 2017 were compared to those occurring in the previous 5 years, from July 2007 to June 2012, using a one tailed t-test, and a significance of p < 0.05.
Norovirus outbreaks: Annual activity and circulating genotypes
A total of 1572 gastroenteritis outbreak investigations were performed in Alberta between July 1st 2012 and February 30th 2018, of which 859 (54.6%) had specimens submitted to the ProvLab for laboratory testing. Norovirus was identified in 530 (61.7%) of all tested outbreaks. The monthly distribution of NoV-positive outbreaks peaked in the winter months (Fig. 1), however, spring peaks with higher activity than that of winter months occurred in March 2014 and May 2016. Compared to historical data, the annual numbers of NoV outbreaks between July 2012 to June 2017 were lower than those observed in the previous 5 years (July 2007 to June 2012 vs. July 2012 to June 2017, p = 0.0489, one-tailed t-test).
Genogroup II strains were responsible for 440 out of 530 (83.0%) of laboratory confirmed NoV outbreaks (Table 1), while genogroup I and mixed genogroup I and genogroup II strains were responsible for 83 (15.7%) and 7 (1.3%) NoV outbreaks, respectively. Twenty ORF2-based genotypes were identified throughout the study period: GI.1, GI.2, GI.3, GI.4, GI.5, GI.6, GI.7, GI.9, GII.1, GII.2, GII.3, GII.4, GII.5, GII.6, GII.7, GII.8, GII.13, GII.14, GII.16 and GII.17 (Table 2). Overall, GII.4 was the most common genotype and was responsible for at least 319 (60.2%) out of 530 NoV-positive outbreaks. Strains carrying an ORF2 of variant Sydney were accountable for the majority of GII.4 outbreaks (297/319, 93.1%). GII.4 variants Den Haag and New Orleans, which caused pandemics in 2006 and 2009, respectively, disappeared after June 2014. Genotypes GI.6 and GI.7 had an increased presence during July 2012 to June 2013 whereas GI.3 was the predominant GI strain from July 2015 to June 2016. Four novel NoV strains emerged during the last three annual periods of study (Fig. 2): 1) GII.17; 2) GII.P16/GII.4 Sydney, a recombinant strain carrying a polymerase gene of GII.16 and a GII.4 Sydney capsid gene; 3) GII.P16/GII.2, a recombinant strain carrying a GII.16 polymerase gene and a GII.2 capsid gene and 4) GII.P4 New Orleans/GII.4 Sydney, a recombinant strain with a GII.4 New Orleans polymerase gene and a GII.4 Sydney capsid gene.
NoV GII.17 outbreaks in Alberta
The first GII.17 outbreak in Alberta occurred in September 2014. Phylogenetic analysis of complete capsid sequences identified all GII.17 strains as part of the Kawasaki cluster (Additional file 1: Figure S1). Representatives of two sub-clusters of GII.17 Kawasaki previously defined by Chan et al.  were identified: strains from 2014 grouped into the Kawasaki323-like sub-cluster (0.7–1% nucleotide and 0.7–0.9% amino acid p-distances with AB983218), and strains from 2015 and 2016 grouped into the Kawasaki308-like sub-cluster (0.4–0.9% nucleotide and 0.2–0.6% amino acid p-distances with LC037415). GII.17 was responsible for 4.3, 11.5 and 25.6 (%) of all NoV outbreaks during the periods of July 2014–June 2015, July 2015–June 2016 and July 2016–June 2017, respectively. No GII.17 outbreaks were observed after April 2017.
Novel GII.P16/GII.4 Sydney recombinant
Several GII.4 strains causing outbreaks in early 2016 could not be classified at the variant level by the genotyping tool based solely on region C sequences. BLAST analysis with additional sequence data from the 3’end of ORF1 (~ 740 nt) revealed a high nucleotide identity (99%) with Kawasaki194-like sequences (GenBank accession number LC175468), previously characterized as GII.P16/GII.4 Sydney recombinant strains . Phylogenetic trees of partial ORF1 and partial ORF2 sequences of this recombinant are shown on Fig. 3. The complete capsid sequences of four GII.P16/GII.4 Sydney recombinants collected between 2016 and 2018, AlbertaEI84/2016, AlbertaEI277/2016, AlbertaEI487/2017 and AlbertaEI231/2018, had 3.3–3.9% nucleotide and 1.5–2.2% amino acid p-distances with the GII.Pe/GII.4 Sydney 2012 reference strain NSW0514/AU/2012 (accession number JX459908) (Additional file 2: Figure S2). Compared with GII.Pe/GII.4 Sydney, the GII.P16/GII.4 Sydney strains did not present distinctive mutations at the P2 domain, however, all GII.P16/GII.4 Sydney sequences in our data set shared two isoleucines at residues 119 and 145, both located at the S domain (Additional file 3: Figure S3). The first GII.P16/GII.4 Sydney outbreak in Alberta occurred in February 2016. Overall, the GII.P16/GII.4 Sydney recombinant was responsible for 21 (27.0%) of 78 NoV outbreaks from July 2015–June 2016, 3 (7.7%) out of 39 NoV outbreaks from July 2016–June 2017 and 44 (62.9%) out of 70 NoV outbreaks from July 2017–February 2018.
Saliva binding profile of the novel GII.P16/GII.4 Sydney recombinant strain
In order to investigate the HBGA binding properties of the novel GII.P16/GII.4 Sydney recombinant, we performed binding assays using p-particles of two representative strains, AB-2016-26 and AB-2016-190 and saliva of individuals with known HBGA phenotypes. Both strains demonstrated ability to bind saliva of A, B and O secretors and did not bind saliva of non-secretors (Additional file 4: Figure S4). AB-2016-190 showed stronger binding to secretors than AB-2016-26, but both strains showed no binding to non-secretors, similar to the binding pattern of Syd9-2B, a GII.Pe/GII.4 Sydney strain.
Novel GII.P16/GII.2 recombinant
Based on retrospective ORF1-genotyping analysis of GII.2 outbreaks back to November 2011, we identified that the first GII.P16/GII.2 outbreak in Alberta occurred in August 2016. In BLAST analysis, the GII.P16/GII.2 outbreak strain shared highest identity (99%) with Volvic-E15317, a strain collected in France in late 2016 (Fig. 4). GII.P16/GII.2 was the predominant NoV strain between July 2016–June 2017 (13 of 39 NoV outbreaks, 33.3%).
GII.P4 New Orleans/GII.4 Sydney recombinant
The recombinant strain GII.P4 New Orleans/GII.4 Sydney was identified during the period July 2016–June 2017 (Fig. 3). The strain was responsible for four out of five GII.4 outbreaks from that period and was not observed after June 2017.
NoV outbreak settings
Between July 2012 and February 2018, the majority of outbreaks occurred in community long-term care, hospital long-term care, supportive living, and group homes facilities (63/83, 75.9% of all NoV GI and 346/440, 78.6% of all NoV GII outbreaks, respectively) and hospital acute care (5/83, 6% of GI and 53/440,12% of GII outbreaks) (Additional file 5: Figure S5). A larger proportion of GI outbreaks occurred in non-health care settings compared to GII (15/83, 18.1% vs. 41/440, 9.3%, p = 0.0180, Chi-square test).
During July 2012 to February 2018 the NoV outbreak activity in Alberta presented important changes in terms of overall disease burden and NoV genotype distribution as compared to historical data. Between July 2012 and June 2017, the annual numbers of NoV outbreaks had a gradual decline. Genotype GII.4 outbreaks, which have previously represented over 50% of all NoV outbreaks, decreased to 23% in July 2016–June 2017. The steep decline in NoV outbreaks was followed by rebound in July 2017–February 2018, when GII.4 outbreaks were accountable for 76% of all NoV outbreaks.
GII.4 Sydney remained as the single predominant strain responsible for a large majority (> 50%) of all NoV outbreaks from July 2012 to June 2015. In contrast, the following 2 years, July 2015 to June 2017, had an atypical high proportion of non-GII.4 outbreaks. Moreover, four novel NoV strains emerged in the province: GII.17 Kawasaki, GII.P16/GII.4 Sydney, GII.P16/GII.2, and GII.P4 New Orleans/GII.4 Sydney.
The GII.17 Kawasaki strain was first reported in Asia in September 2014 where it quickly became the predominant genotype during the 2014–2015-winter season raising concern of a global pandemic . Since its emergence, GII.17 Kawasaki has undergone diversification into at least three sub-clusters . In Alberta, we identified two sub-clusters: the Kawasaki 308-like sub-cluster and the Kawasaki 323-like sub-cluster. Although both clusters started to circulate in the province at similar timeframes as reported in Asia , the GII.17 Kawasaki had limited prevalence in Alberta, as also observed in other regions outside Asia, including Europe and the United States . The factors limiting the transmission of GII.17 Kawasaki outside Asia are largely unknown. The GII.17 Kawasaki strains display a broad saliva HBGA-binding profile comparable to pandemic GII.4 strains , which emphasizes the potential of GII.17 as a pandemic strain.
The novel GII.P16/GII.4 Sydney recombinant was the most common strain among NoV outbreaks from July 2015–June 2016. Although the prevalence of this recombinant decreased temporarily in July 2016–June 2017, its circulation increased sharply in the more recent period, July 2017 and February 2018, when GII.P16/GII.4 Sydney was again responsible for a large majority of outbreaks. The GII.P16/GII.4 Sydney recombinant also circulated in the United States and Germany during the 2015–2016 winter season [18, 26]. Our analysis of capsid sequences from this recombinant collected at three different time points (2016, 2017 and 2018) demonstrated no unique amino acid substitutions compared to GII.Pe/GII.4 Sydney, its ORF2 parental strain. Furthermore, we observed that GII.P16/GII.4 Sydney can bind to saliva of secretor HBGA phenotypes but not to that of non-secretors, similarly as GII.Pe/GII.4 Sydney strains. Our findings agree with a recent phylogenetic study  and suggest that amino acid substitutions outside ORF2, rather than antigenic change at the capsid gene, conferred an evolutionary advantage to GII.P16/GII.4 Sydney over its parental strain GII.Pe/GII.4 Sydney.
The GII.P4 New Orleans/GII.4 Sydney recombinant had low circulation between July 2016 and June 2017 and was no longer observed in the subsequent period. We hypothesize that the transmission of this recombinant in Alberta might have been restricted by its antigenic resemblance to two highly circulated strains, GII.4 New Orleans and GII.4 Sydney. GII.4 has been the predominant genotype worldwide among humans since the 1990s but had very limited circulation in the 1970s and 1980s [28, 29]. It is possible that the decline in GII.4 NoV outbreaks observed between July 2015 and June 2017 is a consequence of herd immunity accumulated in the population after two decades of high prevalence of various GII.4 variants and exhaustion of mutational sites in GII.4 Sydney. The GII.P16/GII.4 Sydney recombinant might have partially circumvented such factors by acquiring the ORF1 genes from the less circulating, and presumably more fit, GII.P16. Further surveillance data is still required to understand the turnover of NoV genotypes in humans.
The last recombinant that emerged in Alberta in the time period of the present study is GII.P16/GII.2. Although GII.P16/GII.2 outbreaks have been reported previously in Asia in 2009, phylogenetic time-scale analyses performed by others suggest that the polymerase of the recent GII.P16/GII.2 strain is rather closely related to that of GII.P16/GII.4 Sydney and carries amino acid substitutions that could have conferred the novel GII.P16/GII.2 an evolutionary advantage . GII.P16/GII.2 strains have also been reported in Germany, France, Japan, China and the United States since mid-2016 [18, 30]. Notably, GII.P16/GII.2 was the predominant strain in the province briefly in July 2016–June 2017, but was replaced thereafter by GII.P16/GII.4 Sydney.
A limitation of our study is the limited collection of ORF1 data for the period between July 2012 and June 2015, which could have lead us to miss novel or unusual recombinants. However, we believe the likelihood of missing important strains is low since 69% of the norovirus GII outbreaks in Alberta within those years were caused by a single ORF2 genotype, GII.4 Sydney. Reports from North America  and diverse countries around different continents  suggest that GII.Pe/GII.4 Sydney was the major strain in circulation worldwide during that time frame and thus, we believe that most, if not all of GII.4 Sydney from that period, carried a GII.Pe polymerase.
In summary, we provided 6 years of systematic molecular surveillance data of norovirus outbreaks in Alberta and identified a GII.4 recombinant, GII.P16/GII.4 Sydney, as the most prevalent strain causing outbreaks in the province in early 2018. An important shift and increase in genotype diversity was noticed in recent years, which should be considered for vaccine development. Ongoing surveillance of the molecular epidemiology of NoV using a dual ORF1/ORF2 genotyping scheme is indispensable to better identify the disease burden of new emerging strains and understand their evolutionary pathways.
Histo-blood group antigens
Open reading frame
Provincial Laboratory for Public Health
Pires SM, Fischer-Walker CL, Lanata CF, Devleesschauwer B, Hall AJ, Kirk MD, Duarte AS, Black RE, Angulo FJ. Aetiology-specific estimates of the global and regional incidence and mortality of Diarrhoeal diseases commonly transmitted through food. PLoS One. 2015;12:e0142927.
Bartsch SM, Lopman BA, Ozawa S, Hall AJ, Lee BY. Global economic burden of norovirus gastroenteritis. PLoS One. 2016;4:e0151219.
Ahmed SM, Hall AJ, Robinson AE, Verhoef L, Premkumar P, Parashar UD, Koopmans M, Lopman BA. Global prevalence of norovirus in cases of gastroenteritis: a systematic review and meta-analysis. Lancet Infect Dis. 2014;8:725–30.
Siebenga JJ, Vennema H, Renckens B, de Bruin E, van der Veer B, Siezen RJ, Koopmans M. Epochal evolution of GGII.4 norovirus capsid proteins from 1995 to 2006. J Virol. 2007;18:9932–41.
Eden JS, Tanaka MM, Boni MF, Rawlinson WD, White PA. Recombination within the pandemic norovirus GII.4 lineage. J Virol. 2013;87:6270–82.
Vinje J. Advances in laboratory methods for detection and typing of norovirus. J Clin Microbiol. 2015;2:373–81.
van Beek J, Ambert-Balay K, Botteldoorn N, Eden JS, Fonager J, Hewitt J, Iritani N, Kroneman A, Vennema H, Vinje J, et al. Indications for worldwide increased norovirus activity associated with emergence of a new variant of genotype II.4, late 2012. Euro Surveill. 2013;18:8–9.
Hasing ME, Lee BE, Preiksaitis JK, Tellier R, Honish L, Senthilselvan A, Pang XL. Emergence of a new norovirus GII.4 variant and changes in the historical biennial pattern of norovirus outbreak activity in Alberta, Canada, from 2008 to 2013. J Clin Microbiol. 2013;7:2204–11.
Siebenga JJ, Vennema H, Zheng DP, Vinje J, Lee BE, Pang XL, Ho EC, Lim W, Choudekar A, Broor S, et al. Norovirus illness is a global problem: emergence and spread of norovirus GII.4 variants, 2001-2007. J Infect Dis. 2009;5:802–12.
Chan MCW, Hu Y, Chen H, Podkolzin AT, Zaytseva EV, Komano J, Sakon N, Poovorawan Y, Vongpunsawad S, Thanusuwannasak T, et al. Global spread of norovirus GII.17 Kawasaki 308, 2014-2016. Emerg Infect Dis. 2017;8:1359–4.
Debbink K, Lindesmith LC, Donaldson EF, Costantini V, Beltramello M, Corti D, Swanstrom J, Lanzavecchia A, Vinje J, Baric RS. Emergence of new pandemic GII.4 Sydney norovirus strain correlates with escape from herd immunity. J Infect Dis. 2013;11:1877–87.
Parra GI, Squires RB, Karangwa CK, Johnson JA, Lepore CJ, Sosnovtsev SV, Green KY. Static and evolving norovirus genotypes: implications for epidemiology and immunity. PLoS Pathog. 2017;1:e1006136.
Yang Y, Xia M, Tan M, Huang P, Zhong W, Pang XL, Lee BE, Meller J, Wang T, Jiang X. Genetic and phenotypic characterization of GII-4 noroviruses that circulated during 1987 to 2008. J Virol. 2010;18:9595–607.
Pang XL, Preiksaitis JK, Lee B. Multiplex real time RT-PCR for the detection and quantitation of norovirus genogroups I and II in patients with acute gastroenteritis. J Clin Virol. 2005;2:168–71.
Kojima S, Kageyama T, Fukushi S, Hoshino FB, Shinohara M, Uchida K, Natori K, Takeda N, Katayama K. Genogroup-specific PCR primers for detection of Norwalk-like viruses. J Virol Methods. 2002;1–2:107–14.
Iritani N, Kaida A, Kubo H, Abe N, Murakami T, Vennema H, Koopmans M, Takeda N, Ogura H, Seto Y. Epidemic of genotype GII.2 noroviruses during spring 2004 in Osaka City, Japan. J Clin Microbiol. 2008;7:2406–9.
Kageyama T, Kojima S, Shinohara M, Uchida K, Fukushi S, Hoshino FB, Takeda N, Katayama K. Broadly reactive and highly sensitive assay for Norwalk-like viruses based on real-time quantitative reverse transcription-PCR. J Clin Microbiol. 2003;4:1548–57.
Cannon JL, Barclay L, Collins NR, Wikswo ME, Castro CJ, Magana LC, Gregoricus N, Marine RL, Chhabra P, Vinje J. 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. 2017;7:2208–21.
Kroneman A, Vennema H, Deforche K, v d Avoort H, Penaranda S, Oberste MS, Vinje J, Koopmans M An automated genotyping tool for enteroviruses and noroviruses. J Clin Virol 2011;2:121–125.
van Beek J, de Graaf M, Al-Hello H, Allen DJ, Ambert-Balay K, Botteldoorn N, Brytting M, Buesa J, Cabrerizo M, Chan M, et al. Molecular surveillance of norovirus, 2005-16: an epidemiological analysis of data collected from the NoroNet network. Lancet Infect Dis. 2018;18:545–53.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;12:2725–9.
Tan M, Jiang X. The p domain of norovirus capsid protein forms a subviral particle that binds to histo-blood group antigen receptors. J Virol. 2005;22:14017–30.
Chan MC, Lee N, Hung TN, Kwok K, Cheung K, Tin EK, Lai RW, Nelson EA, Leung TF, Chan PK. Rapid emergence and predominance of a broadly recognizing and fast-evolving norovirus GII.17 variant in late 2014. Nat Commun. 2015;6:10061.
Matsushima Y, Shimizu T, Ishikawa M, Komane A, Okabe N, Ryo A, Kimura H, Katayama K, Shimizu H. Complete Genome Sequence of a Recombinant GII.P16-GII.4 Norovirus Detected in Kawasaki City, Japan, in 2016. Genome Announc. 2016;4:5.
Zhang XF, Huang Q, Long Y, Jiang X, Zhang T, Tan M, Zhang QL, Huang ZY, Li YH, Ding YQ et al. An outbreak caused by GII.17 norovirus with a wide spectrum of HBGA-associated susceptibility. Sci Rep. 2015;17687.
Niendorf S, Jacobsen S, Faber M, Eis-Hubinger AM, Hofmann J, Zimmermann O, Hohne M, Bock CT. Steep rise in norovirus cases and emergence of a new recombinant strain GII.P16-GII.2, Germany, winter 2016. Euro Surveill. 2017;22:4.
Ruis C, Roy S, Brown JR, Allen DJ, Goldstein RA, Breuer J. The emerging GII.P16-GII.4 Sydney 2012 norovirus lineage is circulating worldwide, arose by late-2014 and contains polymerase changes that may increase virus transmission. PLoS One. 2017;6:e0179572.
Bok K, Abente EJ, Realpe-Quintero M, Mitra T, Sosnovtsev SV, Kapikian AZ, Green KY. Evolutionary dynamics of GII.4 noroviruses over a 34-year period. J Virol. 2009;22:11890–901.
van Beek J, de Graaf M, Xia M, Jiang X, Vinje J, Beersma M, de Bruin E, van de Vijver D, Holwerda M, van Houten M, et al. J Gen Virol. 2016.
Ao Y, Cong X, Jin M, Sun X, Wei X, Wang J, Zhang Q, Song J, Yu J, Cui J, et al. Genetic analysis of reemerging GII.P16-GII.2 noroviruses in 2016-2017 in China. J Infect Dis. 2018;1:133–43.
Pang XL, Preiksaitis JK, Wong S, Li V, Lee BE. Influence of novel norovirus GII.4 variants on gastroenteritis outbreak dynamics in Alberta and the northern territories, Canada between 2000 and 2008. PLoS One. 2010;7:e11599.
Dang Thanh H, Than VT, Nguyen TH, Lim I, Kim W. Emergence of Norovirus GII.17 Variants among Children with Acute Gastroenteritis in South Korea. PLoS One. 2016 May 5;11(5):e0154284.
Lindesmith LC, Brewer-Jensen PD, Mallory ML, Debbink K, Swann EW, Vinje J, Baric RS. Antigenic characterization of a novel recombinant GII.P16-GII.4 Sydney norovirus strain with minor sequence variation leading to antibody escape. J Infect Dis. 2018;7:1145–52.
We thank Andy Wang, Jay Gamma, Joanne Cho and Oliver Jin for their assistance with norovirus sequencing. We also thank the Provincial Laboratory for Public Health staff, especially the virology and molecular departments performing testing for norovirus and the collaborative teams participating in outbreak investigations in Alberta, Northwest Territories, Yukon, and Nunavut, and the members of the Public Health (ProvLab) Outbreak Investigation Committee (OINC).
This work was supported by the Provincial Laboratory for Public Health, Alberta, Canada, and Alberta Health, Canada. The Provincial Laboratory for Public Health (ProvLab) runs the norovirus outbreak molecular surveillance program in Alberta. ProvLab scientists participated in the study design, collection, analysis and interpretation of data as well as writing the manuscript. Alberta Health had no role in the study design, data collection and analysis or manuscript preparation.
Availability of data and materials
The norovirus complete capsid sequences obtained during the current study are available in GenBank under accession numbers MH922873 through MH922876.
Ethics approval and consent to participate
This study was approved by the University of Alberta Biomedical Health Research Ethics Board. The board waived the need of participant consent because samples were anonymized (study ID Pro00037093).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Maximum likelihood phylogeny of GII.17 strains based on complete ORF2 sequences. GII.17 strains from outbreaks in Alberta are shown with triangles. A) Tree based on nucleotide sequences constructed using the Tamura-Nei substitution model assuming gamma-distributed rates of evolution among sites B) Tree based on amino acid sequences constructed using the Jones-Taylor-Thornton (JTT) substitution model assuming gamma-distributed rates of evolution among sites. Both trees were rooted using a GII.13 outgroup reference sequence. Branch significance was estimated based on 1000 bootstrap replicates. Clusters and subclusters are shown as defined in previous studies [31, 32] (PDF 92 kb)
Figure S2. Maximum likelihood phylogeny of complete GII.P16/GII.4 Sydney capsid sequences. The maximum-likelihood tree was constructed using the Tamura-Nei substitution model, assuming gamma-distributed rates of evolution among sites. The analysis included strains representative of different countries and different years. (PDF 14 kb)
Figure S3. Description of data: Analysis of variable sites in capsid sequences of GII.P16/GII.4 Sydney recombinants compared to GII.Pe/GII.4 Sydney strains. Tables depict all variable amino acid residues identified in the alignment used to construct the tree shown in Additional file 2. Evolving amino acid positions previously reported by Lindesmith et al.  to play a role in the evasion of antibody immune responses are shown in red. (PDF 80 kb)
Figure S4. Saliva-binding of NoV P-domain proteins. The P-domain proteins of (A) AB-2016-26 (GII.P16/GII.4 Sydney), (B) AB-2016-190 (GII.P16/GII.4 Sydney) and (C) Syd9-2B (GII.Pe/GII.4 Sydney) were tested in their ability to bind saliva from individuals with different HBGAs profiles. (PDF 255 kb)
Figure S5. Norovirus outbreak settings in Alberta by genogroup. Outbreaks from mixed GI and GII strains (n = 7) were excluded from the analysis. (PDF 23 kb)