Sequencing and phylogenetic analysis of the coding region of six common rotavirus strains: Evidence for intragenogroup reassortment among co-circulating G1P[8] and G2P[4] strains from the United States

Clinical, epidemiological, and genotypic characteristics of rotavirus infection in hospitalized infants and young children in Yunnan Province

Rotaviruses are the most important pathogenic cause of non-bacterial diarrhea in infants and children. Approximately 60% of hospital admissions for acute diarrhea worldwide are caused by rotavirus infection. Rotavirus infection and hospitalization among children in China are a social burden, resulting in economic loss. The prevalence and geographical distribution of rotavirus genotypes is variable, partially due to population migration. Due to the unique geographical conditions and climate in Yunnan Province, several viruses with new genotypes have emerged, and multiple genotypes have become co-epidemic. In this study, rotavirus infection screening and genetic characterization of epidemic strains were performed in 149,492 infants and children admitted to hospitals in six representative prefectures in Yunnan Province between 2019 and 2021. The prevalence of rotavirus infection was 13.39% and was highest in January and lowest in September. G9P[8] was the main epidemic rotavirus genotype. Other epidemic genotypes included G2P[4], G8P[8], G9P[4], G2P[8], G3P[8], G4P[8], G3P[4], and G4P[6]. Phylogenetic analysis revealed that locally epidemic strains were influenced by importation of strains from neighboring provinces and other Asian countries. These findings provide a scientific basis for rotavirus prevention and control and lay a foundation for preliminary studies to establish a rotavirus surveillance network in Yunnan Province.

Whole-Genome Analyses Identifies Multiple Reassortant Rotavirus Strains in Rwanda Post-Vaccine Introduction

Children in low-and middle-income countries, including Rwanda, experience a greater burden of rotavirus disease relative to developed countries. Evolutionary mechanisms leading to multiple reassortant rotavirus strains have been documented over time which influence the diversity and evolutionary dynamics of novel rotaviruses. Comprehensive rotavirus whole-genome analysis was conducted on 158 rotavirus group A (RVA) samples collected pre- and post-vaccine introduction in children less than five years in Rwanda. Of these RVA positive samples, five strains with the genotype constellations G4P[4]-I1-R2-C2-M2-A2-N2-T1-E1-H2 (n = 1), G9P[4]-I1-R2-C2-M2-A1-N1-T1-E1-H1 (n = 1), G12P[8]-I1-R2-C2-M1-A1-N2-T1-E2-H3 (n = 2) and G12P[8]-I1-R1-C1-M1-A2-N2-T2-E1-H1 (n = 1), with double and triple gene reassortant rotavirus strains were identified. Phylogenetic analysis revealed a close relationship between the Rwandan strains and cognate human RVA strains as well as the RotaTeq^® vaccine strains in the VP1, VP2, NSP2, NSP4 and NSP5 gene segments. Pairwise analyses revealed multiple differences in amino acid residues of the VP7 and VP4 antigenic regions of the RotaTeq^® vaccine strain and representative Rwandan study strains. Although the impact of such amino acid changes on the effectiveness of rotavirus vaccines has not been fully explored, this analysis underlines the potential of rotavirus whole-genome analysis by enhancing knowledge and understanding of intergenogroup reassortant strains circulating in Rwanda post vaccine introduction.

Whole genome and in-silico analyses of G1P[8] rotavirus strains from pre- and post-vaccination periods in Rwanda

Rwanda was the first low-income African country to introduce RotaTeq vaccine into its Expanded Programme on Immunization in May 2012. To gain insights into the overall genetic make-up and evolution of Rwandan G1P[8] strains pre- and post-vaccine introduction, rotavirus positive fecal samples collected between 2011 and 2016 from children under the age of 5 years as part of ongoing surveillance were genotyped with conventional RT-PCR based methods and whole genome sequenced using the Illumina MiSeq platform. From a pool of samples sequenced (n = 158), 36 were identified as G1P[8] strains (10 pre-vaccine and 26 post-vaccine), of which 35 exhibited a typical Wa-like genome constellation. However, one post vaccine strain, RVA/Human-wt/RWA/UFS-NGS:MRC-DPRU442/2012/G1P[8], exhibited a RotaTeq vaccine strain constellation of G1-P[8]-I2-R2-C2-M2-A3-N2-T6-E2-H3, with most of the gene segments having a close relationship with a vaccine derived reassortant strain, previously reported in USA in 2010 and Australia in 2012. The study strains segregated into two lineages, each containing a paraphyletic pre- and post-vaccine introduction sub-lineages. In addition, the study strains demonstrated close relationship amongst each other when compared with globally selected group A rotavirus (RVA) G1P[8] reference strains. For VP7 neutralization epitopes, amino acid substitutions observed at positions T91A/V, S195D and M217T in relation to the RotaTeq vaccine were radical in nature and resulted in a change in polarity from a polar to non-polar molecule, while for the VP4, amino acid differences at position D195G was radical in nature and resulted in a change in polarity from a polar to non-polar molecule. The polarity change at position T91A/V of the neutralizing antigens might play a role in generating vaccine-escape mutants, while substitutions at positions S195D and M217T may be due to natural fluctuation of the RVA. Surveillance of RVA at whole genome level will enhance further assessment of vaccine impact on circulating strains, the frequency of reassortment events under natural conditions and epidemiological fitness generated by such events.

Uncovering the First Atypical DS-1-like G1P[8] Rotavirus Strains That Circulated during Pre-Rotavirus Vaccine Introduction Era in South Africa

Emergence of DS-1-like G1P[8] group A rotavirus (RVA) strains during post-rotavirus vaccination period has recently been reported in several countries. This study demonstrates, for the first time, rare atypical DS-1-like G1P[8] RVA strains that circulated in 2008 during pre-vaccine era in South Africa. Rotavirus positive samples were subjected to whole-genome sequencing. Two G1P[8] strains (RVA/Human-wt/ZAF/UFS-NGS-MRC-DPRU1971/2008/G1P[8] and RVA/Human-wt/ZAF/UFS-NGS-MRC-DPRU1973/2008/G1P[8]) possessed a DS-1-like genome constellation background (I2-R2-C2-M2-A2-N2-T2-E2-H2). The outer VP4 and VP7 capsid genes of the two South African G1P[8] strains had the highest nucleotide (amino acid) nt (aa) identities of 99.6–99.9% (99.1–100%) with the VP4 and the VP7 genes of a locally circulating South African strain, RVA/Human-wt/ZAF/MRC-DPRU1039/2008/G1P[8]. All the internal backbone genes (VP1–VP3, VP6, and NSP1-NSP5) had the highest nt (aa) identities with cognate internal genes of another locally circulating South African strain, RVA/Human-wt/ZAF/MRC-DPRU2344/2008/G2P[6]. The two study strains emerged through reassortment mechanism involving locally circulating South African strains, as they were distinctly unrelated to other reported atypical G1P[8] strains. The identification of these G1P[8] double-gene reassortants during the pre-vaccination period strongly supports natural RVA evolutionary mechanisms of the RVA genome. There is a need to maintain long-term whole-genome surveillance to monitor such atypical strains.

Divergence of group a rotavirus with genetic variations before and after introduction of rotavirus vaccines in northern Taiwan

Despite the development of vaccines in 2006, rotavirus is still a major cause of acute gastroenteritis worldwide. This study was performed to analyze the presence of circulating rotaviruses before and after the introduction of rotavirus vaccines to allow phylogenetic comparisons of vaccine strains in northern Taiwan.Rotavirus genotyping and sequencing of rotavirus VP7 and VP4 PCR products were performed by Reverse Transcriptase Polymerase Chain Reaction and DNA autosequencing. Phylogenies were constructed by the neighbor-joining and maximum-likelihood methods using CLUSTAL W software included in the MEGA software package (version 6.0).Between April 2004 and December 2012, a total of 101 rotavirus specimens from pediatric patients with acute gastroenteritis hospitalized in Chang Gung Children's Hospital were amplified, and their VP4 and VP7 sequences were determined. These 101 specimens consisted of 55 pre-vaccine strains (G1 [13, 23.6%], G2 [12, 21.8%], G3 [16, 29.1%], and G9 [14, 25.5%]) and 46 post-vaccine strains (G1 [25, 54.3%], G2 [12, 26.1%], G3 [5, 10.9%], and G9 [4, 8.7%]). The most common combination of the G and P types was G2P[4], accounting for 36% cases, followed by G9P[8] (25%), G1P[8] (20%), G3P[4] (15%), G3P[8] (10%), G1P[4] (5%), and G2P[8] (5%). Phylogenetic analysis showed that only the G1 and P[8] genotypes clustered in the same lineages with the rotavirus vaccine strains.Based on our results, the inclusion of G9, modified G2 and G3 with target lineages, and the combination G2P[4] and G9P[8] in the rotavirus vaccines in Taiwan is warranted as a vaccination strategy.

Rotavirus Burden, Genetic Diversity and Impact of Vaccine in Children under Five in Tanzania

In Tanzania, rotavirus infections are responsible for 72% of diarrhea deaths in children under five. The Rotarix vaccine was introduced in early 2013 to mitigate rotavirus infections. Understanding the disease burden and virus genotype trends over time is important for assessing the impact of rotavirus vaccine in Tanzania. When assessing the data for this review, we found that deaths of children under five declined after vaccine introduction, from 8171/11,391 (72% of diarrhea deaths) in 2008 to 2552/7087 (36% of diarrhea deaths) in 2013. Prior to vaccination, the prevalence of rotavirus infections in children under five was 18.1–43.4%, 9.8–51%, and 29–41% in Dar es Salaam, Mwanza and Tanga, respectively, and after the introduction of vaccines, these percentages declined to 17.4–23.5%, 16–19%, and 10–29%, respectively. Rotaviruses in Tanzania are highly diverse, and include genotypes of animal origin in children under five. Of the genotypes, 10%, 28%, and 7% of the strains are untypable in Dar es Salaam, Tanga, and Zanzibar, respectively. Mixed rotavirus genotype infection accounts for 31%, 29%, and 12% of genotypes in Mwanza, Tanga and Zanzibar, respectively. The vaccine effectiveness ranges between 53% and 75% in Mwanza, Manyara and Zanzibar. Rotavirus vaccination has successfully reduced the rotavirus burden in Tanzania; however, further studies are needed to better understand the relationship between the wildtype strain and the vaccine strain as well as the zoonotic potential of rotavirus in the post-vaccine era.

Rotavirus strain diversity in Eastern and Southern African countries before and after vaccine introduction

BACKGROUND

The African Rotavirus Surveillance Network has been detecting and documenting rotavirus genotypes in the African sub-continent since 1998 in anticipation of the rollout of rotavirus vaccination in routine Expanded Programme on Immunisation. This paper reports distribution of the rotavirus strains circulating in 15 Eastern and Southern African (ESA) countries from 2010-2015 as part of active World Health Organization (WHO) rotavirus surveillance, and investigates possibility of emergence of non-vaccine or unusual strains in six selected countries post-vaccine introduction.

MATERIAL AND METHODS

Stool samples were collected from children <5 years of age presenting with acute gastroenteritis at sentinel hospitals pre- and post-rotavirus vaccine introduction. Samples were tested for group A rotavirus using an enzyme immunoassay by the national and sentinel laboratories. At the WHO Rotavirus Regional Reference Laboratory in South Africa, molecular characterisation was determined by PAGE (n = 4186), G and P genotyping (n = 6447) and DNA sequencing for both G and P types (n = 400).

RESULTS

The six-year surveillance period demonstrated that 23.8% of the strains were G1P[8], followed by G2P[4] (11.8%), G9P[8] (10.4%), G12P[8] (4.9%), G2P[6] (4.2%) and G3P[6] (3.7%) in 15 ESA countries. There was no difference in circulating strains pre- and post-rotavirus vaccine introduction with yearly fluctuation of strains observed over time. Atypical rotavirus G and P combinations (such as G1P[4], G2P[8], G9P[4] and G12P[4]) that might have arisen through inter-genogroup or inter-genotypes reassortment were detected at low frequency (2%). Close genetic relationship of African strains were reflected on the phylogenetic analysis, strains segregated together to form an African cluster in the same lineages/sub-lineage or monophyletic branch.

CONCLUSION

There has been considerable concern about strain replacement post-vaccine introduction, it was not clear at this early stage whether observed cyclical changes of rotavirus strains were due to vaccine pressure or this was just part of natural annual fluctuations in the six ESA countries, long-term surveillance is required.

Whole genomic analysis of G2P[4] human Rotaviruses in Mymensingh, north-central Bangladesh

Rotavirus A (RVA) is a dominant causative agent of acute gastroenteritis in children worldwide. G2P[4] is one of the most common genotypes among human rotavirus (HRV) strains, and has been persistently prevalent in South Asia including Bangladesh. In the present study, whole genome sequences of a total of 16 G2P[4] HRV strains (8 strains each in 2010 and 2013) detected in Mymensingh, north-central Bangladesh were determined. These strains had typical DS-1-like genotype constellation. Most of gene segments from DS-1 genogroup exhibited high level sequence identities to each other (>98%), while slight diversity was observed for VP1, VP3, and NSP4 genes. By phylogenetic analysis, individual RNA segments were classified into one (V) or two-three lineages (V–VI or V–VII). In terms of lineages (sublineages) of 11 gene segments, the 16 Bangladeshi strains could be further classified into four clades (A-D) containing 8 lineage constellations, revealing the presence of three clades (A-C) with three lineage constellations in 2010, and a single clade (D) with four constellations in 2013. Therefore, co-existence of multiple G2P[4] HRV strains with different lineage constellations, and change in clades for the study period were demonstrated. Although amino acids in the antigenic regions on VP7 and VP4 were mostly identical to those of global G2P[4] strains after 2000, VP4 of clade D RVAs in 2013 had alanine and proline at positions 88 and 114, respectively, which are novel substitutions compared with recent global G2P[4] strains. Replacement of lineage constellations associated with unique amino acid changes in the antigenic region in VP4 suggested continuous genetic evolutionary state for emerging new G2P[4] rotavirus strains in Bangladesh.

Whole genome analyses of G1P[8] rotavirus strains from vaccinated and non-vaccinated South African children presenting with diarrhea

Group A rotaviruses (RVAs) are the leading cause of severe gastroenteritis and eventually death among infants and young children worldwide, and disease prevention and management through vaccination is a public health priority. In August 2009, Rotarix™ was introduced in the South African Expanded Programme on Immunisation. As a result, substantial reductions in RVA disease burden have been reported among children younger than 5 years old. Rotavirus strain surveillance post-vaccination is crucial to, inter alia, monitor and study the evolution of vaccine escape strains. Here, full-genome sequence data for the 11 gene segments from 11 South African G1P[8] rotavirus strains were generated, including 5 strains collected from non-vaccinated children during the 2004-2009 rotavirus seasons and 6 strains collected from vaccinated children during the 2010 rotavirus season. These data were analyzed to gain insights into the overall genetic makeup and evolution of South African G1P[8] rotavirus strains and to compare their genetic backbones with those of common human Wa-like RVAs from other countries, as well as with the Rotarix™ and RotaTeq™ G1P[8] vaccine components. All 11 South African G1P[8] strains revealed a complete Wa-like genotype constellation of G1-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1. On the basis of sequence similarities, the South African G1P[8] strains (with the exception of strain RVA/Human-wt/ZAF/1262/2004/G1P[8]) were closely related to each other (96-100% identity in all gene segments). Comparison to the Rotarix™ and RotaTeq™ G1P[8] vaccine components revealed a moderate nucleotide identity of 89-96% and 93-95%, respectively. The results indicated that none of the gene segments of these 11 South African G1P[8] strains were vaccine-derived. This study illustrates that large-scale next generation sequencing will provide crucial information on the influence of the vaccination program on evolution of rotavirus strains. This is the first report to describe full genomic analyses of G1P[8] RVA strains collected from both non-vaccinated and vaccinated children in South Africa.

Detection of novel rotavirus strain by vaccine postlicensure surveillance

Surveillance for rotavirus-associated diarrhea after implementation of rotavirus vaccination can assess vaccine effectiveness and identify disease-associated genotypes. During active vaccine postlicensure surveillance in the United States, we found a novel rotavirus genotype, G14P[24], in a stool sample from a child who had diarrhea. Unusual rotavirus strains may become more prevalent after vaccine implementation.

Molecular epidemiology of contemporary G2P[4] human rotaviruses cocirculating in a single U.S. community: footprints of a globally transitioning genotype

Group A rotaviruses (RVs) remain a leading cause of childhood gastroenteritis worldwide. Although the G/P types of locally circulating RVs can vary from year to year and differ depending upon geographical location, those with G1P[8], G2P[4], G3P[8], G4P[8], G9P[8], and G12P[8] specificities typically dominate. Little is known about the evolution and diversity of G2P[4] RVs and the possible role that widespread vaccine use has had on their increased frequency of detection. To address these issues, we analyzed the 12 G2P[4] RV isolates associated with a rise in RV gastroenteritis cases at Vanderbilt University Medical Center (VUMC) during the 2010-2011 winter season. Full-genome sequencing revealed that the isolates had genotype 2 constellations typical of DS-1-like viruses (G2P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2). Phylogenetic analyses showed that the genome segments of the isolates were comprised of two or three different subgenotype alleles; this enabled recognition of three distinct clades of G2P[4] viruses that caused disease at VUMC in the 2010-2011 season. Although the three clades cocirculated in the same community, there was no evidence of interclade reassortment. Bayesian analysis of 328 VP7 genes of G2 viruses isolated in the last 39 years indicate that existing G2 VP7 gene lineages continue to evolve and that novel lineages, as represented by the VUMC isolates, are constantly being formed. Moreover, G2 lineages are characteristically shaped by lineage turnover events that introduce new globally dominant strains every 7 years, on average. The ongoing evolution of G2 VP7 lineages may give rise to antigenic changes that undermine vaccine effectiveness in the long term.

IMPORTANCE

Little is known about the diversity of cocirculating G2 rotaviruses and how their evolution may undermine the effectiveness of rotavirus vaccines. To expand our understanding of the potential genetic range exhibited by rotaviruses circulating in postvaccine communities, we analyzed part of a collection of rotaviruses recovered from pediatric patients in the United States from 2010 to 2011. Examining the genetic makeup of these viruses revealed they represented three segregated groups that did not exchange genetic material. The distinction between these three groups may be explained by three separate introductions. By comparing a specific gene, namely, VP7, of the recent rotavirus isolates to those from a collection recovered from U.S. children between 1974 and 1991 and other globally circulating rotaviruses, we were able to reconstruct the timing of events that shaped their ancestry. This analysis indicates that G2 rotaviruses are continuously evolving, accumulating changes in their genetic material as they infect new patients.

Evolution of DS-1-like human G2P[4] rotaviruses assessed by complete genome analyses

Group A rotaviruses (RVAs) are a leading cause of viral gastroenteritis in children, with G2P[4] RVA being one of the most common human strains worldwide. The complete genome sequences of nine G2P[4] RVA strains, selected from a 26-year archival collection (1985-2011) established in Palermo, Italy, were determined. A strain associated with a peak of G2P[4] RVA activity in 1996 resembled a reassortant strain identified in Kenya in 1982 and differed completely in genomic make up from more recent strains that circulated during 2004-2011. Conversely, the 2004-2011 G2P[4] RVAs were genetically more similar to contemporary RVA strains circulating globally. Recent G2P[4] strains possessed either single or multiple genome segments (VP1, VP3 and/or NSP4) likely derived from ruminant viruses through intra-genotype reassortment. Amino acid substitutions were selected and maintained over time in the VP7 and VP8* antigenic proteins, allowing the circulation of two contemporary G2P[4] variants to be distinguished. Altogether, these findings suggest that major changes in the genomic composition of recent G2P[4] RVAs occurred in the early 2000s, leading to the appearance of a novel variant of the DS-1-like genotype constellation. Whether the modifications observed in the neutralizing antigens and in the genome composition of modern G2P[4] RVAs may affect the long-term effectiveness of the vaccination programmes remains to be explored.

Comparative evaluation of commercially available manual and automated nucleic acid extraction methods for rotavirus RNA detection in stools

Rotaviruses are a major cause of viral gastroenteritis in children. For accurate and sensitive detection of rotavirus RNA from stool samples by reverse transcription-polymerase chain reaction (RT-PCR), the extraction process must be robust. However, some extraction methods may not remove the strong RT-PCR inhibitors known to be present in stool samples. The objective of this study was to evaluate and compare the performance of six extraction methods used commonly for extraction of rotavirus RNA from stool, which have never been formally evaluated: the MagNA Pure Compact, KingFisher Flex and NucliSENS easyMAG instruments, the NucliSENS miniMAG semi-automated system, and two manual purification kits, the QIAamp Viral RNA kit and a modified RNaid kit. Using each method, total nucleic acid or RNA was extracted from eight rotavirus-positive stool samples with enzyme immunoassay optical density (EIA OD) values ranging from 0.176 to 3.098. Extracts prepared using the MagNA Pure Compact instrument yielded the most consistent results by qRT-PCR and conventional RT-PCR. When extracts prepared from a dilution series were extracted by the 6 methods and tested, rotavirus RNA was detected in all samples by qRT-PCR but by conventional RT-PCR testing, only the MagNA Pure Compact and KingFisher Flex extracts were positive in all cases. RT-PCR inhibitors were detected in extracts produced with the QIAamp Viral RNA Mini kit. The findings of this study should prove useful for selection of extraction methods to be incorporated into future rotavirus detection and genotyping protocols.

Sensitive and specific quantitative detection of rotavirus A by one-step real-time reverse transcription-PCR assay without antecedent double-stranded-RNA denaturation

A real-time quantitative reverse transcription-PCR (qRT-PCR) assay using the recombinant thermostable Thermus thermophilus (rTth) enzyme was developed to detect and quantify rotavirus A (RVA). By using rTth polymerase, significant improvement was achieved over the existing real-time RT-PCR assays, which require denaturation of the RVA double-stranded RNA (dsRNA) prior to assay setup. Using a dsRNA transcript for segment 7, which encodes the assay target NSP3 gene, the limit of detection for the improved assay was calculated to be approximately 1 genome copy per reaction. The NSP3 qRT-PCR assay was validated using a panel of 1,906 stool samples, 23 reference RVA strains, and 14 nontarget enteric virus samples. The assay detected a diverse number of RVA genotypes and did not detect other enteric viruses, demonstrating analytical sensitivity and specificity for RVA in testing stool samples. A XenoRNA internal process control was introduced and detected in a multiplexed qRT-PCR format. Because it does not require an antecedent dsRNA denaturation step, this assay reduces the possibility of sample cross-contamination and requires less hands-on time than other published qRT-PCR protocols for RVA detection.

Human G9P[8] rotavirus strains circulating in Cameroon, 1999-2000: Genetic relationships with other G9 strains and detection of a new G9 subtype

Group A rotaviruses (RV-A) are the leading cause of viral gastroenteritis in children worldwide and genotype G9P[8] is one of the five most common genotypes detected in humans. In order to gain insight into the degree of genetic variability of G9P[8] strains circulating in Cameroon, stool samples were collected during the 1999-2000 rotavirus season in two different geographic regions in Cameroon (Southwest and Western Regions). By RT-PCR, 15 G9P[8] strains (15/89=16.8%) were identified whose genomic configurations was subsequently determined by complete or partial gene sequencing. In general, all Cameroonian G9 strains clustered into current globally-spread sublineages of the VP7 gene and displayed 86.6-100% nucleotide identity amongst themselves and 81.2-99.5% nucleotide identity with global G9 strains. The full genome classification of all Cameroonian strains was G9-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1 but phylogenetic analysis of each gene revealed that the strains were spread across 4 or more distinct lineages. An unusual strain, RVA/Human-wt/CMR/6788/1999/G9P[8], which shared the genomic constellation of other Cameroonian G9P[8] strains, contained a novel G9 subtype which diverged significantly (18.8% nucleotide and 19% amino acid distance) from previously described G9 strains. Nucleotide and amino acid alignments revealed that the 3' end of this gene is highly divergent from other G9 VP7 genes suggesting that it arose through extensive accumulation of point mutations. The results of this study demonstrate that diverse G9 strains circulated in Cameroon during 1999-2000.

G8 rotaviruses with conserved genotype constellations detected in Malawi over 10 years (1997-2007) display frequent gene reassortment among strains co-circulating in humans

Rotavirus A, the most common cause of severe diarrhoea in children worldwide, occurs in five major VP7 (G) and VP4 (P) genotype combinations, comprising G1P[8], G2P[4], G3P[8], G4P[8] and G9P[8]. However, G8, a common bovine rotavirus genotype, has been reported frequently among children in African countries. Surveillance of rotavirus gastroenteritis conducted in a sentinel hospital in Blantyre, Malawi between 1997 and 2007 provided a rare opportunity to examine the whole genotype constellation of G8 strains and their evolution over time. A sample of 27 (9.0 %) of 299 G8 strains was selected to represent each surveillance year and a range of P genotypes, which shifted in predominance from P[6] to P[4] and P[8] during the study period. Following cell culture adaptation, whole genome sequencing demonstrated that the genetic background of 26 strains possessed the DS-1 genotype constellation. A single G8P[6] strain was a reassortant in which both NSP2 and NSP5 genes from strains with the Wa genotype constellation had been inserted into a strain with the DS-1 genotype background. Phylogenetic analysis suggested frequent reassortment among co-circulating strains with the DS-1 genotype constellation. Little evidence was identified to suggest the introduction of contemporary bovine rotavirus genes into any of the 27 G8 strains examined. In conclusion, Malawian G8 strains are closely related to other human strains with the DS-1 genotype constellation. They have evolved over the last decade through genetic reassortment with other human rotaviruses, changing their VP4 genotypes while maintaining a conserved genotype constellation for the remaining structural and non-structural proteins.

Whole genomic analysis of human G1P[8] rotavirus strains from different age groups in China

G1P[8] rotaviruses are an important cause of diarrhea in humans in China. To date, there are no reports on the whole genomic analysis of the Chinese G1P[8] rotaviruses. To determine the origin and overall genetic makeup of the recent Chinese G1P[8] strains, the whole genomes of three strains, RVA/Human-wt/CHN/E1911/2009/G1P[8], RVA/Human-tc/CHN/R588/2005/G1P[8] and RVA/Human-tc/CHN/Y128/2004/G1P[8], detected in an infant, a child and an adult, respectively, were analyzed. Strains E1911, R588 and Y128 exhibited a typical Wa-like genotype constellation. Except for the NSP3 gene of E1911, the whole genomes of strains E1911, R588 and Y128 were found to be more closely related to those of the recent Wa-like common human strains from different countries than those of the prototype G1P[8] strain, or other old strains. On the other hand, the NSP3 gene of E1911 was genetically distinct from those of Y128, R588, or other Wa-like common human strains, and appeared to share a common origin with those of the porcine-like human G9 strains, providing evidence for intergenotype reassortment events. Comparisons of the amino acid residues defining the VP7 and VP4 antigenic domains revealed several mismatches between these Chinese G1P[8] strains and the G1 and P[8] strains contained in the currently licensed rotavirus vaccines Rotarix(TM )and RotaTeq(TM).

Genotype constellation and evolution of group A rotaviruses infecting humans

Numerous rotavirus group A (RVA) strains with distinct G-genotype and P-genotype combinations have been described infecting humans worldwide. However, the increasing amount of complete RVA genome data which have become available, suggest that only RVA strains with 2 discrete genotype constellations have been successful in sustaining infection of humans worldwide over longer periods of time. Those genotype constellations have been designated I1-R1-C1-M1-A1-N1-T1-E1-H1 and I2-R2-C2-M2-A2-N2-T2-E2-H2 and are also known as Wa-like and DS-1-like, respectively. RVAs of other genotype constellations which were able to spread to a limited extent in the human population are AU-1-related RVA strains (I3-R3-C3-M3-A3/A12-N3-T3-E3-H3/H6) in combination with G3P[9] or G12P[9], and neonatal G10P[11] RVA strains in India (bovine×human Wa-like reassortants). On the basis of the analysis of complete genomes, it is suggested that the overall genetic diversity of epidemiologically widespread human RVA strains is more limited than generally assumed. This conclusion has consequences for how we look at host range restriction and the criteria according to which the effectiveness of RVA universal mass vaccination programs is assessed.

Whole genome analyses of African G2, G8, G9, and G12 rotavirus strains using sequence-independent amplification and 454® pyrosequencing

High mortality rates caused by rotaviruses are associated with several strains such as G2, G8, G9, and G12 rotaviruses. Rotaviruses with G9 and G12 genotypes emerged worldwide in the past two decades. G2 and G8 rotaviruses are however also characterized frequently across Africa. To understand the genetic constellation of African G2, G8, G9, and G12 rotavirus strains and their possible origin, sequence-independent cDNA synthesis, amplification, and 454(®) pyrosequencing of the whole genomes of five human African rotavirus strains were performed. RotaC and phylogenetic analysis were used to assign and confirm the genotypes of the strains. Strains RVA/Human-wt/MWI/1473/2001/G8P[4], RVA/Human-wt/ZAF/3203WC/2009/G2P[4], RVA/Human-wt/ZAF/3133WC/2009/G12P[4], RVA/Human-wt/ZAF/3176WC/2009/G12P[6], and RVA/Human-wt/ZAF/GR10924/1999/G9P[6] were assigned G8-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2, G2-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2, G12-P[4]-I1-R1-C1-M1-A1-N1-T1-E1-H1, G12-P[6]-I1-R1-C1-M1-A1-N1-T1-E1-H1, and G9-P[6]-I2-R2-C2-M2-A2-N2-T2-E2-H2 genotypes, respectively. The detection of both Wa- and DS-1-like genotypes in strain RVA/Human-wt/ZAF/3133WC/2009/G12P[4] and Wa-like, DS-1-like and P[6] genotypes in strain RVA/Human-wt/ZAF/GR10924/1999/G9P[6] implies that these two strains were generated through intergenogroup genome reassortment. The close similarity of the genome segments of strain RVA/Human-wt/MWI/1473/2001/G8P[4] to artiodactyl-like, human-bovine reassortant strains and human rotavirus strains suggests that it originated from or shares a common origin with bovine strains. It is therefore possible that this strain might have emerged through interspecies genome reassortment between human and artiodactyl rotaviruses. This study illustrates the swift characterization of all the 11 rotavirus genome segments by using a single set of universal primers for cDNA synthesis followed by 454(®) pyrosequencing and RotaC analysis.

First reports of human rotavirus G8P[4] gastroenteritis in the United States

In 2009, three children were hospitalized in Rochester, NY, with sequence-confirmed G8P[4] rotavirus gastroenteritis-the first U.S. detection of this uncommon strain more typically found in Africa. Continued monitoring of G8P[4] and other rotavirus genotypes not represented in current vaccines is essential to assess whether vaccination will result in an increase in prevalence of these strains.

Whole genome analysis of multiple rotavirus strains from a single stool specimen using sequence-independent amplification and 454® pyrosequencing reveals evidence of intergenotype genome segment recombination

Infection of a single host cell with two or more different rotavirus strains creates conditions favourable for evolutionary mechanisms like reassortment and recombination that can generate novel strains. Despite numerous reports describing mixed rotavirus infections, whole genome characterisation of rotavirus strains in a mixed infection case has not been reported. Double-stranded RNA, exhibiting a long electropherotype pattern only, was extracted from a single human stool specimen (RVA/Human-wt/ZAF/2371WC/2008/G9P[8]). Both short and long electropherotype profiles were however detected in the sequence-independent amplified cDNA derived from the dsRNA, suggesting infection with more than one rotavirus strain. 454® pyrosequencing of the amplified cDNA revealed co-infection of at least four strains. Both genotype 1 (Wa-like) and genotype 2 (DS-1-like) were assigned to the consensus sequences obtained from the nine genome segments encoding NSP1-NSP5, VP1-VP3 and VP6. Genotypes assigned to the genome segments encoding VP4 were P[4] (DS-1-like), P[6] (ST3-like) and P[8] (Wa-like) genotypes. Since four distinct genotypes [G2 (DS-1-like), G8, G9 (Wa-like) and G12] were assigned to the four consensus nucleotide sequences obtained for genome segment 9 (VP7), it was concluded that at least four distinct rotaviruses were present in the stool. Intergenotype genome recombination events were observed in genome segments encoding NSP2, NSP4 and VP6. The close similarities of some of the genome segments encoding NSP2, VP6 and VP7 to artiodactyl rotaviruses suggest that some of the infecting strains shared common ancestry with animal strains, or that interspecies transmission occurred previously. The sequence-independent genome amplification technology coupled with 454® pyrosequencing used in this study enabled the characterisation of the whole genomes of multiple rotavirus strains in a single stool specimen that was previously assigned single genotypes, i.e. G9P[8], by sequence-dependent RT-PCR.

Whole-genomic analysis of rotavirus strains: current status and future prospects

Studies on genetic diversity of rotaviruses have been primarily based on the genes encoding the antigenically significant VP7 and VP4 proteins. Since the rotavirus genome has 11 segments of RNA that are vulnerable to reassortment events, analyses of the VP7 and VP4 genes may not be sufficient to obtain conclusive data on the overall genetic diversity, or true origin of strains. In the last few years following the advent of the whole-genome-based genotype classification system, the whole genomes of at least 167 human group A rotavirus strains have been analyzed, providing a plethora of new and important information on the complex origin of strains, inter- and intra-genogroup reassortment events, animal–human reassortment events, zoonosis, and genetic linkages involving different group A rotavirus gene segments. In addition, the whole genomes of a limited number of human group B, C and novel group rotavirus strains have been analyzed. This article briefly reviews the available data on whole-genomic analysis of human rotavirus strains. The significance and future prospects of whole-genome-based studies are also discussed.

Whole-genome analysis reveals the complex evolutionary dynamics of Kenyan G2P[4] human rotavirus strains

Although G2P[4] rotaviruses are common causes of acute childhood diarrhoea in Africa, to date there are no reports on whole genomic analysis of African G2P[4] strains. In this study, the nearly complete genome sequences of two Kenyan G2P[4] strains, AK26 and D205, detected in 1982 and 1989, respectively, were analysed. Strain D205 exhibited a DS-1-like genotype constellation, whilst strain AK26 appeared to be an intergenogroup reassortant with a Wa-like NSP2 genotype on the DS-1-like genotype constellation. The VP2-4, VP6-7, NSP1, NSP3 and NSP5 genes of strain AK26 and the VP2, VP4, VP7 and NSP1–5 genes of strain D205 were closely related to those of the prototype or other human G2P[4] strains. In contrast, their remaining genes were distantly related, and, except for NSP2 of AK26, appeared to originate from or share a common origin with rotavirus genes of artiodactyl (ruminant and camelid) origin. These observations highlight the complex evolutionary dynamics of African G2P[4] rotaviruses.

Uniformity of rotavirus strain nomenclature proposed by the Rotavirus Classification Working Group (RCWG)

In April 2008, a nucleotide-sequence-based, complete genome classification system was developed for group A rotaviruses (RVs). This system assigns a specific genotype to each of the 11 genome segments of a particular RV strain according to established nucleotide percent cutoff values. Using this approach, the genome of individual RV strains are given the complete descriptor of Gx-P[x]-Ix-Rx-Cx-Mx-Ax-Nx-Tx-Ex-Hx. The Rotavirus Classification Working Group (RCWG) was formed by scientists in the field to maintain, evaluate and develop the RV genotype classification system, in particular to aid in the designation of new genotypes. Since its conception, the group has ratified 51 new genotypes: as of April 2011, new genotypes for VP7 (G20-G27), VP4 (P[28]-P[35]), VP6 (I12-I16), VP1 (R5-R9), VP2 (C6-C9), VP3 (M7-M8), NSP1 (A15-A16), NSP2 (N6-N9), NSP3 (T8-T12), NSP4 (E12-E14) and NSP5/6 (H7-H11) have been defined for RV strains recovered from humans, cows, pigs, horses, mice, South American camelids (guanaco), chickens, turkeys, pheasants, bats and a sugar glider. With increasing numbers of complete RV genome sequences becoming available, a standardized RV strain nomenclature system is needed, and the RCWG proposes that individual RV strains are named as follows: RV group/species of origin/country of identification/common name/year of identification/G- and P-type. In collaboration with the National Center for Biotechnology Information (NCBI), the RCWG is also working on developing a RV-specific resource for the deposition of nucleotide sequences. This resource will provide useful information regarding RV strains, including, but not limited to, the individual gene genotypes and epidemiological and clinical information. Together, the proposed nomenclature system and the NCBI RV resource will offer highly useful tools for investigators to search for, retrieve, and analyze the ever-growing volume of RV genomic data.