key: cord-011105-or9azf1g authors: Huang, Zheng; Yao, Dong; Xiao, Shan; Yang, Dong; Ou, Xinhua title: Full-genome sequences of GII.13[P21] recombinant norovirus strains from an outbreak in Changsha, China date: 2020-04-30 journal: Arch Virol DOI: 10.1007/s00705-020-04643-1 sha: doc_id: 11105 cord_uid: or9azf1g On 31 March 2019, 68 school students suffered from vomiting, diarrhea, and abdominal pain after participating in a group activity at a commercial park. In this outbreak, multiple norovirus genotypes were observed, including GII.2[P16], GII.17[P17], and GII.13[P21]. Further, we determined the full-genome sequences of two strains of GII.13[P21] recombinant noroviruses, which were 7434 nt long. Phylogenetic analysis based on open reading frames (ORFs) 1 and 2 revealed that these recombinants were related to stains of different genotypes from different countries. The full genome nucleotide sequences of the two isolates were 97.0% and 98.0% identical to those of strains from London and Thailand, respectively. Simplot analysis revealed the presence of a break point at nt 5059 in the ORF1 region. The histo-blood group antigen binding sites were conserved in both recombinant viruses. Our findings not only provide valuable genetic information about a recombinant norovirus but also contribute to our general understanding of the evolution, genetic diversity, and distribution of noroviruses. Noroviruses are a major causative agent of acute gastroenteritis, with high prevalence across the globe [1] . Norovirus outbreaks in public spaces, such as kindergartens and primary or secondary schools, are generally associated with low hygiene levels and contaminated food or water [2] [3] [4] . Noroviruses belong to the family Caliciviridae and have three open reading frames (ORFs) in their genome. ORF1 encodes a large nonstructural protein (Ntpase, p22, VPg, 3CLpro and RdRp), ORF2 encodes the major structural capsid protein (VP1), and ORF3 encodes a minor structural protein (VP2) [5] . Based on the amino acid sequence of the capsid protein, noroviruses have been classified into 10 genogroups (GI~GX) and approximately 53 genotypes so far [6, 7] . Of the genogroups, GI, GII and GIV can be found in human infections, with genogroup GII being the most common [8] . The GII. 4 [P4] genotype has been reported in many countries and is the prevalent strain in most outbreaks and human infections [9] [10] [11] . Previous reports have shown that GII. 17 [P17] became the main epidemic genotype of norovirus in China, in 2014 [11] [12] [13] [14] . GII. 21 [P21] and GII. 13 [P13] have also been found in the coastal environment and in seafood in Korea. These two genotypes share a close phylogenetic relationship with GII.17 [P17] . However, these genotypes are found infrequently in human infections, and this might be associated with a unique histo-blood group antigen (HBGA) binding site involved in host susceptibility and the presence of decoy glycan receptors in the human gastrointestinal tract that prevents binding of the virus [15, 16] . Due to the high frequency of genetic recombination at the ORF1 and ORF2 junction, a number of recombinant strains have emerged under natural selection [1, 17, 18] [20] [21] [22] [23] . Although those genotypes have not caused severe outbreaks, studying these strains has helped us to obtain more information on the genetic diversity and gene constellation of noroviruses. Here, we report two GII. 13 [P21] recombinant strains from an outbreak in the city of Changsha. In this outbreak, there were also two cases of coinfection with GI and GII. To better understand the genetic background and molecular characteristics of these viruses, we carried out a comprehensive analysis of the complete genome sequences of these noroviruses obtained by next-generation sequencing (NGS). The stool sample in this study was provided by the YueLu District Center for Disease Prevention and Control and tested by the Changsha Center for Disease Prevention and Control (CSCDC). The epidemiological and clinical information on laboratory-confirmed human cases of norovirus infection were collected by the CSCDC staff and medical doctors at the hospital. The Institutional Review Board reviewed and approved the use of those samples for this research (no. CSCDC-2019-008). Six stool samples were obtained from patients who suffered from fever, vomiting and diarrhea on 1 April 2019. The samples were stored at −70°C until RNA extraction. Viral RNA was extracted using a QIAamp RNA Mini Kit (QIAGEN), and detected using a norovirus real-time PCR kit (Jiangsu Bioperfectus Technologies, China) according to the manufacturer's instructions. Then, the region corresponding to the ORF1 and ORF2 junction was amplified from norovirus-positive samples by reverse transcription polymerase chain reaction (RT-PCR) using norovirus type I and II ORF(1+2) junction region gene amplification kits (BioGerm, China). ORF1 and 2 junction sequences were obtained from the company (BioGerm, China). The genotypes of the norovirus sequences were determined using the Norovirus Genotyping Tool website (http://www.rivm.nl/ mpf/norov irus/typin gtool ) [7] . The full genome was amplified from of norovirus-positive samples by RT-PCR using a SuperScript TM III One- Step RT-PCR System with a Platinum TM Taq High Fidelity DNA Polymerase Kit (Thermo Fisher, USA), performed in a GeneAmp PCR System 9700 (Thermo Fisher, USA) [24, 25] . PCR products were purified using AMPure XP beads (Beckman, USA) according to the manufacturer's instructions and eluted using 45 μl of nucleic-acid-free water. The purified nucleic acid sequences were quantitated by Qubit 2.0 using a dsDNA HS (High Sensitivity) Assay Kit. A Nextera XT DNA Library Prep Kit (Illumina, USA) was used to construct a DNA library with 1 ng of input DNA. Then, the samples were sequenced using a Miseq v2 Reagent Kit (Illumina, USA) on a Miseq platform (Illumina, USA). The sequence data were analyzed using Fastqc, Cutadapt and Virus Identification Pipeline (VIP) software. Sequences were assembled using SPAdes-3.13.0 software. Sequence alignments were performed with Clustal W using Molecular Evolutionary Genetic Analysis software version 6 (MEGA 6). A phylogenetic tree based on fullgenome sequences was constructed by the maximumlikelihood method in MEGA 6 (1000 bootstrap replicates). Phylogenetic trees based on RdRp and VP1 sequences were constructed by the neighbor-joining method with the Kimura two-parameter model in MEGA 6 (1000 bootstrap replicates). The other complete and partial genome sequences were downloaded from NCBI. To identify break points in the genomes of the recombinant strains, their sequences were analyzed using SimPlot 3.5.1 software. The analysis was conducted using 1500-bp sequences from the ORF1 and ORF2 regions with a window size of 200 nt and a step size of 20 nt [24] . Amino acid sequences and HBGA binding sites were analyzed using Biological Sequence Alignment Editor (BioEdit 7.0.5). The outbreak occurred in a senior high school after 786 students participated in a group activity at a commercial park. On 29 March 2019, the students had lunch and dinner at the commercial park, and the first case of illness with nausea and vomiting was reported on 30 March. In this outbreak, 68 cases were reported, 31 males and 37 females, including one teacher. The infection cases were distributed in 13 classes. The symptoms in students who suffered from gastroenteritis mainly included dizziness (48.53%), nausea (75.00%), vomiting (83.82%), diarrhea (57.35%) and abdominal pain (47.06%). In all cases, the symptoms were resolved within 48-72 h, with no severe outcome. To identify the pathogens causing acute gastroenteritis in this outbreak, six stool specimens from patients, 10 anal swab samples from the park staff, and food and drinking water samples were collected for norovirus testing. The results showed that the six stool specimens were GII positive, two of which exhibited a mixed infection with GI and GII. However, other samples from the commercial park were negative. The ORF1 and ORF2 junction was amplified from six samples by RT-PCR and sequenced, and the following genotypes were found: GII. 13 To investigate the genetic relationship between the recombinants and other sequences available in the GenBank database, full-genome sequences and ORF1-2 fragments of strains were analyzed by constructing a phylogenetic tree (Fig. 1A) . A phylogenetic tree based on a fragment of the RdRp gene revealed that the recombinant strains belonged to the same branch as viruses from the United Kingdom (MH218651) and Bhutan (MH702263) and shared 98.4%-98.5% sequence identity with MH702263 ( Fig. 2A) . The ORF2 fragment of these strains belonged to the GII.13 genotype branch and shared 98.2% sequence identity with MG892908 (Fig. 2B) . The complete genome sequences of the GII. 13 [15] . This site consists of eight residues located at the top of each P domain. The eight residues are as follows: N297 and W298 of the B loop; S357, T359, and S360 of the N loop; and N395, N397, and T398 of the T loop. None of these residues were mutated in the GII. 13 [P21] recombinant strains in this study. However, other substitutions that might affect the structure of the binding site were also observed -N294S, N309S, and V394Q -as well as a V541A substitution in the C-terminal amino acid sequence. The HBGA binding sites were conserved in the GII. 13 [P21] recombinant strains. T135A and K139R substitutions were found in the RdRp segment of these viruses. GI and GII noroviruses are the most important viral cause of non-bacterial gastroenteritis in humans globally. A large proportion of the population is infected with noroviruses via contact with contaminated food, water, and environment every year [17, 18, 26] . Various [15] . Hence, GII.21[P21] and GII. 13 [P13] genotypes are still evolving and spreading at a lower rate. Although these genotypes might infect humans only sporadically, it is still useful to study the genetic diversity of these norovirus genotypes. Due to the lack of in vitro cell culture systems and in vivo animal models for human noroviruses, there is still a lack of information on its pathogenesis and epidemiology. Whole-genome sequencing of recombinant strains is important for vaccine development strategies and studies of viral evolution. Hence, we determined the genome sequences of GII.13[P21] recombinant strains isolated in the city of Changsha, China. From the epidemiological information, although only visitors of the commercial park suffered from gastroenteritis, all food and water samples from the park were free of contamination, indicating that the source of the contamination might have been direct contact with infected persons or contaminated environments. Based on the genetic and sequence data, we conclude that the outbreak involved multiple genotypes, including GII.4[P4], GII.2[P16], and GII.17 [P17] , which have been reported frequently in previous outbreaks [10] . These results reveal that the prevalent genotypes continue to infect humans. Phylogenetic analysis of the ORF1 and ORF2 segments showed that the recombinant strains belong to different branches. These results suggest that sequences in these recombinant strains might have originated in other countries neighboring China or that the strains might have evolved from local strains that have not been detected in the environment before. The high degree of sequence similarity between the two GII. 13 [13] . The break point identified in this study was located in ORF1 segment, at the same position as reported previously in other recombinant strains [23] . Based on these findings, it is suggested that recombination between these genotypes might occur more easily than with other genotypes. Although the genotypes GII. 21 [P21] and GII. 13 [P13] represent a new evolutionary lineage of norovirus selected by HBGAs, the binding sites of the recombinant strains are still conserved [16] . The role of substitutions in the B, N and T loops of the P domain, including N294S, N309S and V394Q, which are associated with adaptation of the HBGA binding site, needs to be clarified in future research. TheT135A and K139R substitutions in the ORF1 sequence and the V541A substitution in ORF2 sequences may have no influence on the structure of the binding pocket. Although the sporadic norovirus genotypes do not cause serious epidemic worldwide, unlike GII.2[P16] and GII.4[P4], due to these special binding sites, human infections are associated with severe vomiting and diarrhea. Hence, it is important for us to enrich the database of norovirus genome sequences. Although only a few samples were collected in this outbreak, we obtained the full genome of GII.13[P21] recombinant strains that have rarely been reported in China. Due to the high rate of genetic exchange in the norovirus genome, the virus can escape immune monitoring in individuals and acquire new host specificity [31] . Our research is important for understanding the diversity and wide distribution of noroviruses. Evolution of norovirus Genotypic and epidemiologic trends of norovirus outbreaks in the United States Environmental Surveillance for Noroviruses in Selected South African Wastewaters Outbreaks of norovirus and acute gastroenteritis associated with British Columbia Oysters Structure(s), function(s), and inhibition of the RNA-dependent RNA polymerase of noroviruses Advances in laboratory methods for detection and typing of norovirus Updated classification of norovirus genogroups and genotypes Proposal for a unified norovirus nomenclature and genotyping Mechanisms of GII.4 norovirus evolution Emergence and predominance of norovirus GII Emergence of a new GII.17 norovirus variant in patients with acute gastroenteritis in Gastroenteritis Outbreaks Caused by Norovirus GII An outbreak of gastroenteritis associated with GII.17 Norovirus-contaminated secondary water supply system in Wuhan Noroviruses Recognize Glycans with a Terminal beta-Galactose via an Unconventional Glycan Binding Site A Unique Human Norovirus Lineage with a Distinct HBGA Binding Interface Human norovirus transmission and evolution in a changing world Epidemiology of Norovirus outbreaks reported to the public health emergency event surveillance system, China A norovirus intervariant GII.4 recombinant in Victoria Emergence of norovirus GII.P16-GII.2 strains in patients with acute gastroenteritis in Huzhou Emergence of Norovirus GII.17 variants among children with acute gastroenteritis in South Korea Analysis of uncommon norovirus recombinants from Manaus, Amazon region, Brazil: GII.P22/GII.5, GII.P7/GII.6 and GII Novel recombinant GII.P16_GII.13 and GII.P16_GII.3 norovirus strains in Italy Full-genomic analysis of a human norovirus recombinant GII.12/13 novel strain isolated from South Korea Full-genome sequence analysis of an uncommon norovirus genotype, GII.21, from South Korea Detection of norovirus and rotavirus present in suspended and dissolved forms in drinking water sources Emergence of GII.4 Sydney norovirus in South Korea during the winter of 2012-2013 Distribution of norovirus and sapovirus genotypes with emergence of NoV GII.P16/GII.2 recombinant strains in Chiang Mai, Thailand GII.21 in children with Diarrhea Occurrence of novel GII.17 and GII.21 norovirus variants in the coastal environment of South Korea in 2015 Norovirus-host interaction: implications for disease control and prevention Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations Acknowledgements We acknowledge all of the clinicians who collected data and samples.Funding This study was supported by the Hunan Provincial Health Commission Foundation (No. B2017220).Norovirus strains from an outbreak in Changsha, China The authors declare that they have no conflict of interest.