Sequence analysis of the whole genomes of five African human G9 rotavirus strains

Genomic Analysis of Rwandan G9P[8] Rotavirus Strains Pre- and Post-RotaTeq® Vaccine Reveals Significant Distinct Sub-Clustering in a Post-Vaccination Cohort

Although the introduction of rotavirus vaccines has substantially contributed to the reduction in rotavirus morbidity and mortality, concerns persist about the re-emergence of variant strains that might alter vaccine effectiveness in the long term. The G9 strains re-emerged in Africa during the mid-1990s and have more recently become predominant in some countries, such as Ghana and Zambia. In Rwanda, during the 2011 to 2015 routine surveillance period, G9P[8] persisted during both the pre- and post-vaccine periods. The pre-vaccination cohort was based on the surveillance period of 2011 to 2012, and the post-vaccination cohort was based on the period of 2013 to 2015, excluding 2014. The RotaTeq® vaccine that was first introduced in Rwanda in 2012 is genotypically heterologous to Viral Protein 7 (VP7) G9. This study elucidated the whole genome of Rwandan G9P[8] rotavirus strains pre- and post-RotaTeq® vaccine introduction. Fecal samples from Rwandan children under the age of five years (pre-vaccine n = 23; post-vaccine n = 7), conventionally genotyped and identified as G9P[8], were included. Whole-genome sequencing was then performed using the Illumina® MiSeq platform. Phylogenetic analysis and pair-wise sequence analysis were performed using MEGA6 software. Distinct clustering of three post-vaccination study strains was observed in all 11 gene segments, compared to the other Rwandan G9P[8] study strains. Specific amino acid differences were identified across the gene segments of these three 2015 post-vaccine strains. Important amino acid differences were identified at position N242S in the VP7 genome segment of the three post-vaccine G9 strains compared to the other G9 strains. This substitution occurs at a neutralization epitope site and may slightly affect protein interaction at that position. These findings indicate that the Rwandan G9P[8] strains revealed a distinct sub-clustering pattern among post-vaccination study strains circulating in Rwanda, with changes at neutralization epitopes, which may play a role in neutralization escape from vaccine candidates. This emphasizes the need for continuous whole-genome surveillance to better understand the evolution and epidemiology of the G9P[8] strains post-vaccination.

Whole Genome Analysis of African G12P[6] and G12P[8] Rotaviruses Provides Evidence of Porcine-Human Reassortment at NSP2, NSP3, and NSP4

Group A rotaviruses (RVA) represent the most common cause of pediatric gastroenteritis in children <5 years, worldwide. There has been an increase in global detection and reported cases of acute gastroenteritis caused by RVA genotype G12 strains, particularly in Africa. This study sought to characterize the genomic relationship between African G12 strains and determine the possible origin of these strains. Whole genome sequencing of 34 RVA G12P[6] and G12P[8] strains detected from the continent including southern (South Africa, Zambia, Zimbabwe), eastern (Ethiopia, Uganda), central (Cameroon), and western (Togo) African regions, were sequenced using the Ion Torrent PGM method. The majority of the strains possessed a Wa-like backbone with consensus genotype constellation of G12-P[6]/P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1, while a single strain from Ethiopia displayed a DS-1-like genetic constellation of G12-P[6]-I2-R2-C2-M2-A2-N2-T2-E2-H2. In addition, three Ethiopian and one South African strains exhibited a genotype 2 reassortment of the NSP3 gene, with genetic constellation of G12-P[8]-I1-R1-C1-M1-A1-N1-T2-E1-H1. Overall, 10 gene segments (VP1–VP4, VP6, and NSP1–NSP5) of African G12 strains were determined to be genetically related to cognate gene sequences from globally circulating human Wa-like G12, G9, and G1 strains with nucleotide (amino acid) identities in the range of 94.1–99.9% (96.5–100%), 88.5–98.5% (93–99.1%), and 89.8–99.0% (88.7–100%), respectively. Phylogenetic analysis showed that the Ethiopian G12P[6] possessing a DS-1-like backbone consistently clustered with G2P[4] strains from Senegal and G3P[6] from Ethiopia with the VP1, VP2, VP6, and NSP1–NSP4 genes. Notably, the NSP2, NSP3, and NSP4 of most of the study strains exhibited the closest relationship with porcine strains suggesting the occurrence of reassortment between human and porcine strains. Our results add to the understanding of potential roles that interspecies transmission play in generating human rotavirus diversity through reassortment events and provide insights into the evolutionary dynamics of G12 strains spreading across selected sub-Saharan Africa regions.

Report of the 1st African Enteric Viruses Genome Initiative (AEVGI) Data and Bioinformatics Workshop on whole-genome analysis of some African rotavirus strains held in Bloemfontein, South Africa

The University of the Free State - Next Generation Sequencing (NGS) Unit, Bloemfontein, South Africa, hosted a data and bioinformatics workshop from 19 to 22 June 2018. The workshop was coordinated by the African Enteric Viruses Genome Initiative (AEVGI) with support from the Bill & Melinda Gates Foundation. The event introduced technologies in NGS and data analysis with focus on the rotavirus (RV) genome. The workshop fostered interactions and networking between professionals, scientific experts, technicians and students. The courses provided an overview of RV diarrhoea and its burden in Africa, while highlighting the key resources and methodologies in NGS and advanced bioinformatics in deciphering vaccine impact. It was concluded that, despite the reported significant decline in RV associated-diarrhoea mortality and morbidity in Africa due to RV vaccine impact, the need for continuous surveillance and genomic characterization to better understand the ever-changing dynamics of RV strains is imperative.

Rotavirus G9P[4], G9P[6] and G1P[6] strains isolated from children with acute gastroenteritis in Pune, western India, 2013-2015: evidence for recombination in genes encoding VP3, VP4 and NSP1

Species A rotaviruses (RVAs) are genetically diverse pathogens. These are the most evolutionarily adaptable organisms, with a multitude of mechanisms for evolutionary change. To date, full-genome classification has been proved to be an excellent tool for studying the evolution of unusual rotavirus strains. As limited data are available from Pune (Maharashtra), western India, the current study was undertaken with the aim of understanding the genetic diversity in three (G1P[6], G9P[4] and G9P[4]) unusual RVA strains circulating in Pune, India during 2013-2015. Full-genome analysis of these strains classified them as G1-P[6]-I1-R1-C1-M1-A1-N1-T1-E1-H1, G9-P[4]-I2-R2-C2-[M1-M2_R]-[A1-A2_R]-N2-T2-E6-H2 and G9-[P4-P6_R]-I1-R1-C1-M1-A1-N1-T1-E1-H1. Sequencing and phylogenetic analysis of the structural and non-structural genes of these unusual RVA strains showed nucleotide/amino acid identities of 82.3-98.5 %/77.3-99.8 % and 86.6-97.6 %/89.6-97.8 % between the strains of the study. Evidence of recombination events was found within the genes encoding VP3, VP4 and NSP1, which showed a combination of genetic information for genogroup 1 [M1/P[6]/A1] and genogroup 2 [M2/P[4]/A2] strains. This study will facilitate future investigations into the molecular pathogenesis of such RVAs as the exchange of whole or partial genetic material between rotaviruses through recombination contributes directly to their diversification, adaptation and evolution.

Multiple Introductions and Antigenic Mismatch with Vaccines May Contribute to Increased Predominance of G12P[8] Rotaviruses in the United States

Rotavirus is the leading global cause of diarrheal mortality for unvaccinated children under 5 years of age. The outer capsid of rotavirus virions consists of VP7 and VP4 proteins, which determine viral G and P types, respectively, and are primary targets of neutralizing antibodies. Successful vaccination depends upon generating broadly protective immune responses following exposure to rotaviruses presenting a limited number of G- and P-type antigens. Vaccine introduction resulted in decreased rotavirus disease burden but also coincided with the emergence of uncommon G and P genotypes, including G12. To gain insight into the recent predominance of G12P[8] rotaviruses in the United States, we evaluated 142 complete rotavirus genome sequences and metadata from 151 clinical specimens collected in Nashville, TN, from 2011 to 2013 through the New Vaccine Surveillance Network. Circulating G12P[8] strains were found to share many segments with other locally circulating strains but to have distinct constellations. Phylogenetic analyses of G12 sequences and their geographic sources provided evidence for multiple separate introductions of G12 segments into Nashville, TN. Antigenic epitopes of VP7 proteins of G12P[8] strains circulating in Nashville, TN, differ markedly from those of vaccine strains. Fully vaccinated children were found to be infected with G12P[8] strains more frequently than with other rotavirus genotypes. Multiple introductions and significant antigenic mismatch may in part explain the recent predominance of G12P[8] strains in the United States and emphasize the need for continued monitoring of rotavirus vaccine efficacy against emerging rotavirus genotypes.IMPORTANCE Rotavirus is an important cause of childhood diarrheal disease worldwide. Two immunodominant proteins of rotavirus, VP7 and VP4, determine G and P genotypes, respectively. Recently, G12P[8] rotaviruses have become increasingly predominant. By analyzing rotavirus genome sequences from stool specimens obtained in Nashville, TN, from 2011 to 2013 and globally circulating rotaviruses, we found evidence of multiple introductions of G12 genes into the area. Based on sequence polymorphisms, VP7 proteins of these viruses are predicted to present themselves to the immune system very differently than those of vaccine strains. Many of the sick children with G12P[8] rotavirus in their diarrheal stools also were fully vaccinated. Our findings emphasize the need for continued monitoring of circulating rotaviruses and the effectiveness of the vaccines against strains with emerging G and P genotypes.

Phylogenetic analysis of human G9P[8] rotavirus strains circulating in Jiangsu, China between 2010 and 2016

Rotavirus A (RVA) is the leading cause of acute viral gastroenteritis in children under 5 years of age worldwide. G9P[8] is a common RVA genotype that has been persistently prevalent in Jiangsu, China. To determine the genetic diversity of G9P[8] RVAs, 7 representative G9P[8] strains collected from Suzhou Children's Hospital between 2010 and 2016 (named JS2010-JS2016) were analyzed through whole-genome sequencing. All evaluated strains showed the Wa-like constellation G9-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1. Furthermore, phylogenetic analysis revealed that the VP7 genes of all strains clustered into lineage G9-III and G9-VI. With the exception of strain JS2012 (P[8]-4), the VP4 sequences of all strains belonged to the P[8]-3 lineage. Sequencing further revealed that amino acid substitutions were present in the antigenic regions of the VP7 and VP4 genes of all strains. Moreover, there were multiple substitutions in antigenic sites I and II of the nonstructural protein 4 (NSP4) genes, whereas the other NSP genes were relatively conserved. In conclusion, our phylogenetic analysis of these 7 G9P[8] strains suggests that RVA varied across regions and time. Therefore, our findings suggest that continued surveillance is necessary to explore the molecular evolutionary characteristics of RVA for better prevention and treatment of acute viral gastroenteritis.

Whole-genome sequencing and analyses identify high genetic heterogeneity, diversity and endemicity of rotavirus genotype P[6] strains circulating in Africa

Rotavirus A (RVA) exhibits a wide genotype diversity globally. Little is known about the genetic composition of genotype P[6] from Africa. This study investigated possible evolutionary mechanisms leading to genetic diversity of genotype P[6] VP4 sequences. Phylogenetic analyses on 167 P[6] VP4 full-length sequences were conducted, which included six porcine-origin sequences. Of the 167 sequences, 57 were newly acquired through whole genome sequencing as part of this study. The other 110 sequences were all publicly-available global P[6] VP4 full-length sequences downloaded from GenBank. The strength of association between the phenotypic features and the phylogeny was also determined. A number of reassortment and mixed infections of RVA genotype P[6] strains were observed in this study. Phylogenetic analyses demostrated the extensive genetic diversity that exists among human P[6] strains, porcine-like strains, their concomitant clades/subclades and estimated that P[6] VP4 gene has a higher substitution rate with the mean of 1.05E-3 substitutions/site/year. Further, the phylogenetic analyses indicated that genotype P[6] strains were endemic in Africa, characterised by an extensive genetic diversity and long-time local evolution of the viruses. This was also supported by phylogeographic clustering and G-genotype clustering of the P[6] strains when Bayesian Tip-association Significance testing (BaTS) was applied, clearly supporting that the viruses evolved locally in Africa instead of spatial mixing among different regions. Overall, the results demonstrated that multiple mechanisms such as reassortment events, various mutations and possibly interspecies transmission account for the enormous diversity of genotype P[6] strains in Africa. These findings highlight the need for continued global surveillance of rotavirus diversity.

Completely genomic and evolutionary characteristics of human-dominant G9P[8] group A rotavirus strains in Yunnan, China

G9P[8] rotavirus A (RVA) has been identified as the predominant genotype circulating in Yunnan, China. To elucidate the molecular characteristics of its genetic composition at the whole-genome level, the genomes of 12 strains isolated from paediatric patients with diarrhoea were fully sequenced and characterized. Eleven of the 12 strains were genotyped as G9-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1, which is closely related to the Wa-like genotype 1 constellation. In contrast, one strain was genotyped as G9-P[8]-I1-R1-C1-M1-A1-N2-T1-E1-H1, with the NSP2 gene characterized as a DS-1 like genotype. Bayesian phylogenetic analysis indicated that G9 strains had emerged in 1932 with an estimated average evolutionary rate of 1.63×10-3 substitutions/site/year. Considering the high prevalence and fast evolutionary rate of G9P[8] rotaviruses, our results suggest that G9P[8] RVA should be strictly monitored in China.

Distinguishing the genotype 1 genes and proteins of human Wa-like rotaviruses vs. porcine rotaviruses

Group A rotaviruses (RVAs) are 11-segmented, double-stranded RNA viruses and important causes of gastroenteritis in the young of many animal species. Previous studies have suggested that human Wa-like RVAs share a close evolutionary relationship with porcine RVAs. Specifically, the VP1-VP3 and NSP2-5/6 genes of these viruses are usually classified as genotype 1 with >81% nucleotide sequence identity. Yet, it remains unknown whether the genotype 1 genes and proteins of human Wa-like strains are distinguishable from those of porcine strains. To investigate this, we performed comprehensive bioinformatic analyses using all known genotype 1 gene sequences. The RVAs analyzed represent wildtype strains isolated from humans or pigs at various geographical locations during the years of 2004-2013, including 11 newly-sequenced porcine RVAs from Brazil. We also analyzed archival strains that were isolated during the years of 1977-1992 as well as atypical strains involved in inter-species transmission between humans and pigs. We found that, in general, the genotype 1 genes of typical modern human Wa-like RVAs clustered together in phylogenetic trees and were separate from those of typical modern porcine RVAs. The only exception was for the NSP5/6 gene, which showed no host-specific phylogenetic clustering. Using amino acid sequence alignments, we identified 34 positions that differentiated the VP1-VP3, NSP2, and NSP3 genotype 1 proteins of typical modern human Wa-like RVAs versus typical modern porcine RVAs and documented how these positions vary in the archival/unusual isolates. No host-specific amino acid positions were identified for NSP4, NSP5, or NSP6. Altogether, the results of this study support the notion that human Wa-like RVAs and porcine RVAs are evolutionarily related, but indicate that some of their genotype 1 genes and proteins have diverged over time possibly as a reflection of sequestered replication and protein co-adaptation in their respective hosts.

Whole genome detection of rotavirus mixed infections in human, porcine and bovine samples co-infected with various rotavirus strains collected from sub-Saharan Africa

Group A rotaviruses (RVA) are among the main global causes of severe diarrhea in children under the age of 5years. Strain diversity, mixed infections and untypeable RVA strains are frequently reported in Africa. We analysed rotavirus-positive human stool samples (n=13) obtained from hospitalised children under the age of 5years who presented with acute gastroenteritis at sentinel hospital sites in six African countries, as well as bovine and porcine stool samples (n=1 each), to gain insights into rotavirus diversity and evolution. Polyacrylamide gel electrophoresis (PAGE) analysis and genotyping with G-(VP7) and P-specific (VP4) typing primers suggested that 13 of the 15 samples contained more than 11 segments and/or mixed G/P genotypes. Full-length amplicons for each segment were generated using RVA-specific primers and sequenced using the Ion Torrent and/or Illumina MiSeq next-generation sequencing platforms. Sequencing detected at least one segment in each sample for which duplicate sequences, often having distinct genotypes, existed. This supported and extended the PAGE and RT-PCR genotyping findings that suggested these samples were collected from individuals that had mixed rotavirus infections. The study reports the first porcine (MRC-DPRU1567) and bovine (MRC-DPRU3010) mixed infections. We also report a unique genome segment 9 (VP7), whose G9 genotype belongs to lineage VI and clusters with porcine reference strains. Previously, African G9 strains have all been in lineage III. Furthermore, additional RVA segments isolated from humans have a clear evolutionary relationship with porcine, bovine and ovine rotavirus sequences, indicating relatively recent interspecies transmission and reassortment. Thus, multiple RVA strains from sub-Saharan Africa are infecting mammalian hosts with unpredictable variations in their gene segment combinations. Whole-genome sequence analyses of mixed RVA strains underscore the considerable diversity of rotavirus sequences and genome segment combinations that result from a complex evolutionary history involving multiple host species.

Whole genomic analysis of human G12P[6] and G12P[8] rotavirus strains that have emerged in Kenya: identification of porcine-like NSP4 genes

G12 rotaviruses are globally emerging rotavirus strains causing severe childhood diarrhea. However, the whole genomes of only a few G12 strains have been fully sequenced and analyzed, of which only one G12P[4] and one G12P[6] are from Africa. In this study, we sequenced and characterized the complete genomes of three G12 strains (RVA/Human-tc/KEN/KDH633/2010/G12P[6], RVA/Human-tc/KEN/KDH651/2010/G12P[8], and RVA/Human-tc/KEN/KDH684/2010/G12P[6]) identified in three stool specimens from children with acute diarrhea in Kenya, Africa. On whole genomic analysis, all three Kenyan G12 strains were found to have a Wa-like genetic backbone: G12-P[6]-I1-R1-C1-M1-A1-N1-T1-E1-H1 (strains KDH633 and KDH684) and G12-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1 (strain KDH651). Phylogenetic analysis showed that most genes of the three strains examined in this study were genetically related to globally circulating human G1, G9, and G12 strains. Of note is that the NSP4 genes of strains KDH633 and KDH684 appeared to be of porcine origin, suggesting the occurrence of reassortment between human and porcine strains. Furthermore, strains KDH633 and KDH684 were very closely related to each other in all the 11 gene segments, indicating derivation of the two strains from a common origin. On the other hand, strain KDH651 consistently formed distinct clusters of 10 of the 11 gene segments (VP1-2, VP4, VP6-7, and NSP1-5), indicating a distinct origin of strain KDH651 from that of strains KDH633 and KDH684. To our knowledge, this is the first report on whole genome-based characterization of G12 strains that have emerged in Kenya. Our observations will provide important insights into the evolutionary dynamics of emerging G12 rotaviruses in Africa.

Genetic diversity of porcine group A rotavirus strains in the UK

•

This is the first study of rotavirus genotypes circulating in UK pigs.

•

Rotavirus transmission between pigs and humans is not thought to be common in the UK.

•

Human rotavirus genotype P[8] found in a UK pig.

•

The uncommon rotavirus genotype P[32] is widespread in UK pig herds.

Rotavirus is endemic in pig farms where it causes a loss in production. This study is the first to characterise porcine rotavirus circulating in UK pigs. Samples from diarrheic pigs with rotavirus enteritis obtained between 2010 and 2012 were genotyped in order to determine the diversity of group A rotavirus (GARV) in UK pigs. A wide range of rotavirus genotypes were identified in UK pigs: six G types (VP7); G2, G3, G4, G5, G9 and G11 and six P types (VP4); P[6], P[7], P[8], P[13], P[23], and P[32]. With the exception of a single P[8] isolate, there was less than 95% nucleotide identity between sequences from this study and any available rotavirus sequences.

The G9 and P[6] genotypes are capable of infecting both humans and pigs, but showed no species cross-over within the UK as they were shown to be genetically distinct, which suggested zoonotic transmission is rare within the UK. We identified the P[8] genotype in one isolate, this genotype is almost exclusively found in humans. The P[8] was linked to a human Irish rotavirus isolate in the same year. The discovery of human genotype P[8] rotavirus in a UK pig confirms this common human genotype can infect pigs and also highlights the necessity of surveillance of porcine rotavirus genotypes to safeguard human as well as porcine health.

Rotavirus surveillance in Kisangani, the Democratic Republic of the Congo, reveals a high number of unusual genotypes and gene segments of animal origin in non-vaccinated symptomatic children

Group A rotavirus (RVA) infections form a major public health problem, especially in low-income countries like the Democratic Republic of the Congo (COD). However, limited data on RVA diversity is available from sub-Saharan Africa in general and the COD in particular. Therefore, the first aim of this study was to determine the genetic diversity of 99 RVAs detected during 2007–2010 in Kisangani, COD. The predominant G-type was G1 (39%) and the most predominant P-type was P[6] (53%). A total of eight different G/P-combinations were found: G1P[8] (28%), G8P[6] (26%), G2P[4] (14%), G12P[6] (13%), G1P[6] (11%), G9P[8] (4%), G4P[6] (2%) and G8P[4] (1%). The second aim of this study was to gain insight into the diversity of P[6] RVA strains in the COD. Therefore, we selected five P[6] RVA strains in combination with the G1, G4, G8 (2x) or G12 genotype for complete genome analysis. Complete genome analysis showed that the genetic background of the G1P[6] and G12P[6] strains was entirely composed of genotype 1 (Wa-like), while the segments of the two G8P[6] strains were identified as genotype 2 (DS-1-like). Interestingly, all four strains possessed a NSP4 gene of animal origin. The analyzed G4P[6] RVA strain was found to possess the unusual G4-P[6]-I1-R1-C1-M1-A1-N1-T7-E1-H1 constellation. Although the majority of its genes (if not all), were presumably of porcine origin, this strain was able to cause gastro-enteritis in humans. The high prevalence of unusual RVA strains in the COD highlights the need for continued surveillance of RVA diversity in the COD. These results also underline the importance of complete genetic characterization of RVA strains and indicate that reassortments and interspecies transmission among human and animal RVAs strains occur regularly. Based on these data, RVA vaccines will be challenged with a wide variety of different RVA strain types in the COD.

Whole-genome analyses of DS-1-like human G2P[4] and G8P[4] rotavirus strains from Eastern, Western and Southern Africa

Group A rotaviruses (RVAs) with distinct G and P genotype combinations have been reported globally. We report the genome composition and possible origin of seven G8P[4] and five G2P[4] human RVA strains based on the genetic evolution of all 11 genome segments at the nucleotide level. Twelve RVA ELISA positive stool samples collected in the representative countries of Eastern, Southern and West Africa during the 2007-2012 surveillance seasons were subjected to sequencing using the Ion Torrent PGM and Illumina MiSeq platforms. A reference-based assembly was performed using CLC Bio's clc_ref_assemble_long program, and full-genome consensus sequences were obtained. With the exception of the neutralising antigen, VP7, all study strains exhibited the DS-1-like genome constellation (P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2) and clustered phylogenetically with reference strains having a DS-1-like genetic backbone. Comparison of the nucleotide and amino acid sequences with selected global cognate genome segments revealed nucleotide and amino acid sequence identities of 81.7-100 % and 90.6-100 %, respectively, with NSP4 gene segment showing the most diversity among the strains. Bayesian analyses of all gene sequences to estimate the time of divergence of the lineage indicated that divergence times ranged from 16 to 44 years, except for the NSP4 gene where the lineage seemed to arise in the more distant past at an estimated 203 years ago. However, the long-term effects of changes found within the NSP4 genome segment should be further explored, and thus we recommend continued whole-genome analyses from larger sample sets to determine the evolutionary mechanisms of the DS-1-like strains collected in Africa.

Analysis of human rotaviruses from a single location over an 18-year time span suggests that protein coadaption influences gene constellations

Rotaviruses (RVs) are 11-segmented, double-stranded RNA viruses that cause severe gastroenteritis in children. In addition to an error-prone genome replication mechanism, RVs can increase their genetic diversity by reassorting genes during host coinfection. Such exchanges allow RVs to acquire advantageous genes and adapt in the face of selective pressures. However, reassortment may also impose fitness costs if it unlinks genes/proteins that have accumulated compensatory, coadaptive mutations and that operate best when kept together. To better understand human RV evolutionary dynamics, we analyzed the genome sequences of 135 strains (genotype G1/G3/G4-P[8]-I1-C1-R1-A1-N1-T1-E1-H1) that were collected at a single location in Washington, DC, during the years 1974 to 1991. Intragenotypic phylogenetic trees were constructed for each viral gene using the nucleotide sequences, thereby defining novel allele level gene constellations (GCs) and illuminating putative reassortment events. The results showed that RVs with distinct GCs cocirculated during the vast majority of the collection years and that some of these GCs persisted in the community unchanged by reassortment. To investigate the influence of protein coadaptation on GC maintenance, we performed a mutual information-based analysis of the concatenated amino acid sequences and identified an extensive covariance network. Unexpectedly, amino acid covariation was highest between VP4 and VP2, which are structural components of the RV virion that are not thought to directly interact. These results suggest that GCs may be influenced by the selective constraints placed on functionally coadapted, albeit noninteracting, viral proteins. This work raises important questions about mutation-reassortment interplay and its impact on human RV evolution.

IMPORTANCE

Rotaviruses are devastating human pathogens that cause severe diarrhea and kill >450,000 children each year. The virus can evolve by accumulating mutations and by acquiring new genes from other strains via a process called reassortment. However, little is known about the relationship between mutation accumulation and gene reassortment for rotaviruses and how it impacts viral evolution. In this study, we analyzed the genome sequences of human strains found in clinical fecal specimens that were collected at a single hospital over an 18-year time span. We found that many rotaviruses did not reassort their genes but instead maintained them as specific sets (i.e., constellations). By analyzing the encoded proteins, we discovered concurrent amino acid changes among them, which suggests that they are functionally coadapted to operate best when kept together. This study increases our understanding of how rotaviruses evolve over time in the human population.

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.

Genetic diversity of rotavirus genome segment 6 (encoding VP6) in Pretoria, South Africa

BACKGROUND

Rotavirus viral protein 6 (VP6), encoded by genome segment (GS) 6, is the primary target for rotavirus diagnosis by serological and some molecular techniques. Selected full length nucleotide sequences of GS 6 of rotavirus strains from South Africa were sequenced and analysed to determine genetic diversity and variations within the circulating rotaviruses.

FINDINGS

The VP6 amplicons were sequenced using the Sanger ABI 3130xl. Phylogenetic and pairwise analysis revealed that the VP6 genes of the study strains belonged to two different VP6 [I] genotypes. Five sequences were assigned genotype I1 and seven as genotype I2. Comparison of the group specific antigenic regions of the South African strains to the reference strains, shows that the South African VP6 sequences belonging to the VP6 genotype I2 were highly conserved, with only two amino acids changes at positions 239 (T›N) and 261(I›V). On the other hand, South African VP6 sequences belonging to I1 genotypes revealed several amino acid variations mostly within the antigenic region III.

CONCLUSIONS

Rotavirus strains with I1 and I2 genotype are predominantly circulating within the South African communities of which the later seems to be more conserved within the antigenic regions. The observed genetic variations observed within GS 6 of rotaviruses analysed in the current study are unlikely to impact negatively on the performance of the current VP6-based detection methods. Nevertheless, investigators should continually consider this diversity and adapt the primer design for the detection and characterization of the VP6 gene accordingly.