Full genome-based classification of rotaviruses reveals a common origin between human Wa-Like and porcine rotavirus strains and human DS-1-like and bovine rotavirus strains

The performance of biological and tertiary wastewater treatment procedures for rotaviruses A removal

Enteric viruses, generally found in sewage, are recognized as the main cause of waterborne and foodborne public health outbreaks. Among leading enteric viruses, the Rotavirus A (RVA) detection in wastewater appeared to be a novel approach to monitor the emergence of these viruses in some countries where the viral gastroenteritis surveillance is almost absent such as in Tunisia. The RVA detection and quantification in an industrial sewage purification plant of Charguia I (Tunis, Tunisia) were achieved to evaluate the performance of activated sludge procedures coupled to a macrofiltration monolamp ultraviolet irradiation type C (UV-C₂₅₄) disinfection reactor. This UV-C₂₅₄ system was preceded by a fiberglass cartridge filter system with an average porosity of 45 μm to clarify the water and thus increase its UV transmittance. A total of 140 composite sewage samples was collected from this line of treatment and analyzed for RVA detection. The detection and the viral load quantification of RVA were performed using real-time reverse transcription polymerase chain reaction (RT-PCR). The virological results showed in general that RVA were detected at high frequency of 98% (137/140). In fact, the RVA detection rates at the exit of the two studied wastewater treatment were about 100% at the exit of the activated sludge procedure. It means that all wastewater sampled at this last step of treatment was positive for RVA detection. On the other hand, 92.5% of the wastewater samples taken at the exit of the monolamp UV-C₂₅₄ reactor were positive for the RVA. However, the RVA quantification results expressed as viral load showed a significant reduction in the means of RVA viral loads at the exit of the biological activated sludge procedure and the tertiary UV-C₂₅₄ treatment, showing in general an improved treated wastewater virological quality. Therefore, the RVA load removal rates recorded at the two successive stages of treatment, the activated sludge and the UV-C₂₅₄ treatment, were around 85% and 73%, respectively, as compared to the one with 100% registered for the raw wastewater. In addition, good physical-chemical and bacteriological qualities of the treated sewage were found at the exit of the two considered wastewater treatment procedures. The present investigation represents the first Tunisian environmental report showing the good effectiveness and performance of the biological and the tertiary treatments for RVA removal. Therefore, an improved and an optimized tertiary disinfection treatment was needed since it could be a good means for getting better viral water quality and for minimizing the transmission and dissemination of human infectious viral diseases.

Reference Human Rotavirus A Genome Sequence from a Previously Vaccinated Child with Diarrhea in Nigeria

In 2018, a 26-month-old girl, fully vaccinated with Rotarix in 2016, presented with fever, diarrhea, and vomiting. A rapid test showed that her feces contained rotavirus A (RVA). VP7 reverse transcription-PCR (RT-PCR) and Illumina sequencing showed that a G1P[8] strain with a Wa-like genotype constellation was the etiologic agent. This is the first near-complete RVA genome sequence from Nigeria.

In 2018, a 26-month-old girl, fully vaccinated with Rotarix in 2016, presented with fever, diarrhea, and vomiting. A rapid test showed that her feces contained rotavirus A (RVA). VP7 reverse transcription-PCR (RT-PCR) and Illumina sequencing showed that a G1P[8] strain with a Wa-like genotype constellation was the etiologic agent. This is the first near-complete RVA genome sequence from Nigeria.

Vaccine Effectiveness against DS-1-Like Rotavirus Strains in Infants with Acute Gastroenteritis, Malawi, 2013-2015

Atypical DS-1-like G1P[8] rotaviruses emerged in 2013 in Malawi after rotavirus vaccine introduction. Vaccine effectiveness among infants hospitalized with acute DS-1-like G1P[8] rotavirus gastroenteritis was 85.6% (95% CI 34.4%-96.8%). These findings suggest that vaccine provides protection against these strains despite their emergence coinciding with vaccine introduction.

Rotavirus Replication: Gaps of Knowledge on Virus Entry and Morphogenesis

In 1973, rotaviruses A (RVAs) were discovered as major causative agents of acute gastroenteritis in infants and young children worldwide. The infectious RV virion is an icosahedral particle composed of three concentric protein layers surrounding the 11 double-stranded (dsRNA) segments. An in vitro replication system for RVs in permanent cell lines was developed in 1982 and expanded to replication in intestinal organoids in 2015. However, the details of rotavirus (RV) entry into cells and particle maturation mechanisms at the molecular level remain incompletely understood. Slowing down human RVA replication in cell culture on ice allowed morphological visualization of virus particle entry and the assembly of triple-layered particles (virion). Although RVAs are non-enveloped viruses, after virus attachment to the cell membrane, the virus enters the cell by perforating the plasma membrane by a fusion mechanism involving VP5* of the cleaved VP4 protein, as the alternative virus entry route besides the receptor-mediated endocytosis which is generally accepted. After assembling double-layered particles (DLPs) in viroplasm or cytoplasm, they appear to be connected with the endoplasmic reticulum (ER) membrane and become coated with outer capsid proteins (VP4 and VP7) in a coating process. The perforation of the ER membrane is caused by an unknown mechanism following interaction between non-structural protein 4 (NSP4) and the inner capsid protein VP6 of the DLPs. The coating process is closely related to the formation of a hetero-oligomeric complex (NSP4, VP4 and VP7). These lines of evidence suggest the existence of novel mechanisms of RV morphogenesis.

Intragenic recombination influences rotavirus diversity and evolution

Because of their replication mode and segmented dsRNA genome, homologous recombination is assumed to be rare in the rotaviruses. We analyzed 23,627 complete rotavirus genome sequences available in the NCBI Virus Variation database, and found 109 instances of homologous recombination, at least eleven of which prevailed across multiple sequenced isolates. In one case, recombination may have generated a novel rotavirus VP1 lineage. We also found strong evidence for intergenotypic recombination in which more than one sequence strongly supported the same event, particularly between different genotypes of segment 9, which encodes the glycoprotein, VP7. The recombined regions of many putative recombinants showed amino acid substitutions differentiating them from their major and minor parents. This finding suggests that these recombination events were not overly deleterious, since presumably these recombinants proliferated long enough to acquire adaptive mutations in their recombined regions. Protein structural predictions indicated that, despite the sometimes substantial amino acid replacements resulting from recombination, the overall protein structures remained relatively unaffected. Notably, recombination junctions appear to occur nonrandomly with hot spots corresponding to secondary RNA structures, a pattern seen consistently across segments. In total, we found strong evidence for recombination in nine of eleven rotavirus A segments. Only segments 7 (NSP3) and 11 (NSP5) did not show strong evidence of recombination. Collectively, the results of our computational analyses suggest that, contrary to the prevailing sentiment, recombination may be a significant driver of rotavirus evolution and may influence circulating strain diversity.

Full genomic characterization of a porcine rotavirus strain detected in an asymptomatic piglet in Accra, Ghana

Background

The introduction of rotavirus A vaccination across the developing world has not proved to be as efficacious as first hoped. One cause of vaccine failure may be infection by zoonotic rotaviruses that are very variable antigenically from the vaccine strain. However, there is a lack of genomic information about the circulating rotavirus A strains in farm animals in the developing world that may be a source of infection for humans. We therefore screened farms close to Accra, Ghana for animals sub-clinically infected with rotavirus A and then sequenced the virus found in one of these samples.

Results

6.1% of clinically normal cows and pigs tested were found to be Rotavirus A virus antigen positive in the faeces. A subset of these (33.3%) were also positive for virus RNA. The most consistently positive pig sample was taken forward for metagenomic sequencing. This gave full sequence for all open reading frames except segment 5 (NSP1), which is missing a single base at the 5′ end. The virus infecting this pig had genome constellation G5-P[7]-I5-R1-C1-M1-A8-N1-T7-E1-H1, a known porcine genotype constellation.

Conclusions

Farm animals carry rotavirus A infection sub-clinically at low frequency. Although the rotavirus A genotype discovered here has a pig-like genome constellation, a number of the segments most closely resembled those isolated from humans in suspected cases of zoonotic transmission. Therefore, such viruses may be a source of variable gene segments for re-assortment with other viruses to cause vaccine breakdown. It is recommended that further human and pig strains are characterized in West Africa, to better understand this dynamic.

Pre-vaccine rotavirus surveillance in Buenos Aires, Argentina. Characterization of an emergent G1P[8] strain associated to fatal cases in 2014

Group A rotaviruses (RVA) are the most frequent etiological agents causing severe diarrhea in infants and surveillance of genotype, and genetic characteristics of circulating strains are necessary in order to evaluate vaccine programs. The objectives of this work were to describe G and P genotype from 2012 through 2014 in Buenos Aires, Argentina completing an overview of 19 years of genotype surveillance in our region and to characterize an emerging G1P[8] strain associated with severe cases and five fatalities in 2014. We performed genotyping by RT-PCR. The sequencing of several genes, phylogenetic analyses, and comparative epidemiological data were used to know the origin and phylogenetic relationships of the emerging G1P[8] strain. Along with this report, 19 years of continuous RVA genotype surveillance in Argentina in the pre-vaccine era was covered. During the last year of this surveillance, 2014, a significantly increased incidence of RVA associated gastroenteritis was related to the reemergence of G1P[8] strains, being these ones detected in low frequency in the last nine years. Interestingly, the patients affected were significantly older when compared with those from the last six seasons. Additionally, phylogenetic analysis of several genes infer that these G1P[8] strains were closely related to Asian strains circulating during 2012 and 2013. In addition to this, the suggested extra continental origin for the 2014 G1P[8] strains and the very low circulation of G1 type during nine years probably explain the increased incidence and severity in the gastroenteritis cases and the particular epidemiologic characteristics. In conclusion, this work gives us a whole panorama of the pre-vaccine era of the RVA molecular epidemiology in the most populated region of Argentina. In this way, this work inspires us to continue with this type of studies in the post-vaccination era.

Reverse Genetics System for a Human Group A Rotavirus

Group A rotavirus (RV) is a major cause of acute gastroenteritis in infants and young children worldwide. Recently, we established an entirely plasmid-based reverse genetics system for simian RV strain SA11. Although that system was robust enough to generate reassortant RVs, including human RV gene segments, and enabled better understanding of the biological differences between animal and human RV strains, a complete reverse genetics system for human RV strains is desirable. Here, we established a plasmid-based reverse genetics system for G4P[8] human RV strain Odelia. This technology was used to generate a panel of monoreassortant viruses between human and simian RV strains for all of the 11 gene segments demonstrating full compatibility between human and simian RV strains. Furthermore, we generated recombinant viruses lacking the C-terminal region of the viral nonstructural protein NSP1 and used it to define the biological function of the interaction between NSP1 and its target protein β-transducin repeat-containing protein (β-TrCP) during viral replication. While the NSP1 truncation mutant lacking the C-terminal 13 amino acids displayed lower β-TrCP degradation activity, it replicated as efficiently as the wild-type virus. In contrast, the truncation mutant lacking the C-terminal 166 amino acids of NSP1 replicated poorly, suggesting that the C-terminal region of NSP1 plays critical roles in viral replication. The system reported here will allow generation of engineered recombinant virus harboring desired mutations, increase our understanding of the molecular biology of human RV, and facilitate development of novel therapeutics and vaccines.IMPORTANCE Reverse genetics, an approach used to generate viruses from cloned cDNA, has increased our understanding of virus biology. Worldwide research led to the development of an entirely plasmid-based reverse genetics system for the simian RV laboratory strain. Although the technique allows generation of gene-modified recombinant RVs, biological differences between animal and human RVs mean that reverse genetics systems for human RV strains are still needed. Here, we describe a reverse genetics system for the high-yield human RV strain Odelia, which replicates efficiently and is suitable for in vitro molecular studies. Monoreassortant viruses between simian and human RV strains and NSP1 mutant viruses generated by the rescue system enabled study of the biological functions of viral gene segments. This human RV reverse genetics system will facilitate study of RV biology and development of vaccines and vectors.

Coding-Gene Coevolution Analysis of Rotavirus Proteins: A Bioinformatics and Statistical Approach

Rotavirus remains a major cause of diarrhea in infants and young children worldwide. The permanent emergence of new genotypes puts the potential effectiveness of vaccines under serious question. The distribution of unusual genotypes subject to viral fitness is influenced by interactions among viral proteins. The present work aimed at analyzing the genetic constellation and the coevolution of rotavirus coding genes for the available rotavirus genotypes. Seventy-two full genome sequences of different genetic constellations were analyzed using a genetic algorithm. The results revealed an extensive genome-wide covariance network among the 12 viral proteins. Altogether, the emergence of new genotypes represents a challenge to the outcome and success of vaccination and the coevolutionary analysis of rotavirus proteins may boost efforts to better understand the interaction networks of proteins during viral replication/transcription.

Genetic variability of pig and human rotavirus group A isolates from Slovakia

The aim of this work was the genetic typing of RVA isolates originating from pigs and human patients in Slovakia. Seventy-eight rectal swabs from domestic pigs and 30 stool samples from humans were collected. The whole VP7 (G genotypes), VP6 (I genotypes) and partial VP4 (P genotypes) ORFs were amplified by RT-PCR. Genetic variability was higher amongst porcine sequences, where four G genotypes (G3, G4, G5, G11), two P genotypes (P[6], P[13]) and one I5 genotype were detected. Human RVA strains were represented by two G genotypes (G1, G3), two I genotypes (I1, I2), and one P genotype (P[8]). Genetic analysis did not show a relationship between Slovakian porcine and human RVA strains, but phylogenetic grouping of some Slovakian porcine sequences with Hungarian human sequences in both G and P genotypes was observed.

Characterization of rotavirus possessing a DS-1-like VP3 gene from pigs in Brazil: Evidence for zooanthroponotic transmission

Porcine group A rotavirus (RVA) strains SUI15A and SUI24A are suggested to have VP3 genes of human origin possessing DS-1-like backbone. The aim of the present study was to analyse the genome of two strains (SUI15A and SUI24A) and understand the evolution of a rare human-like M2 genotype in pigs. On partial genomic analysis, strains SUI24A (G3-P[13]-I5-R1-C1-M2-A8-N1-T7-E1-H1) and SUI15A (G3-P[x]-Ix-R1-C1-M2-Ax-Nx-T7-E1-H1) were found to have VP3 gene RVA different from those of typical porcine RVA strains described in Brazil and worldwide. This genotypic constellation was a novel constellation that has not been reported previously in both humans and pigs. Furthermore, on phylogenetic analysis, VP3 gene of strains appeared to be of human origin. Therefore, suggested to have evidence for human-to-porcine zooanthroponotic transmission.

Shedding of oral pentavalent bovine-human reassortant rotavirus vaccine indicates high uptake rate of vaccine and prominence of G-type G1

BACKGROUND

Live oral pentavalent bovine-human reassortant rotavirus (RV) vaccine, RotaTeq®, contains bovine rotaviruses reassorted with human G-types G1, G2, G3 and G4, and P-type P[8]. Shedding of RotaTeq® vaccine, as studied by RT-PCR, has been shown to be more common than initially reported, and may include formation of vaccine-derived double-reassortant G1P[8] RVs. We studied the extent and duration of RotaTeq® vaccine virus shedding, genotypes shed, and clinical symptoms associated with shedding.

MATERIAL AND METHODS

We enrolled a total of 301 infants who received RotaTeq® vaccine according to Finnish schedule at 2, 3 and 5 months of age. Stool samples were collected 5-10 days after the first and 0-7 days before the third dose of the vaccine. Additional stool samples 6 and 12 weeks later were collected if the second stool sample was positive. All stools were studied with RT-PCR for RV VP7, VP4 and VP6. Parents filled a symptom diary for a week after each vaccine dose.

RESULTS

We found that 93% of the vaccinees shed vaccine related viral particles in one sample taken 5-10 days after the first dose, indicating that stool shedding is very common and may be regarded as a marker of successful vaccination. Genotype G1 was the predominant genotype in shedding, often in association with P[8], and the only genotype found in long-term shedding. Also G4 was commonly detected whereas other vaccine G-types and bovine-type P[5] were not.

CONCLUSIONS

Shedding of RotaTeq® vaccine-derived viruses is a sign for successful vaccination. Intense shedding of G1 with or without P[8]reflects effective multiplication and may be an important factor in the induction of protective immunity. Shedding of G1 containing vaccine viruses may be prolonged up to 8 months of age. These results suggest that the pentavalent vaccine functions largely like a monovalent G1 vaccine. Eudra-CT: 2014-004252-60.

Molecular characterization of group A rotavirus from rhesus macaques (Macaca mulatta) at human-wildlife interfaces in Bangladesh

Group A rotavirus (RVA) is an important cause of diarrhoea in people, especially children, and animals globally. Due to the segmented nature of the RVA genome, animal RVA strains have the potential to adapt to the human host through reassortment with other co-infecting human viruses. Macaques share food and habitat with people, resulting in close interaction between these two species. This study aimed to detect and characterize RVA in rhesus macaques (Macaca mulatta) in Bangladesh. Faecal samples (N = 454) were collected from apparently healthy rhesus macaques from nine different sites in Bangladesh between February and March 2013. The samples were tested by one-step, real-time, reverse transcriptase-polymerase chain reaction (PCR). Four percent of samples (n = 20; 95% CI 2.7%-6.7%) were positive for RVA. RVA positive samples were further characterized by nucleotide sequence analysis of two structural protein gene fragments, VP4 (P genotype) and VP7 (G genotype). G3, G10, P[3] and P[15] genotypes were identified and were associated as G3P[3], G3P[15] and G10P[15]. The phylogenetic relationship between macaque RVA strains from this study and previously reported human strains indicates possible transmission between humans and macaques in Bangladesh. To our knowledge, this is the first report of detection and characterization of rotaviruses in rhesus macaques in Bangladesh. These data will not only aid in identifying viral sharing between macaques, human and other animals, but will also improve the development of mitigation measures for the prevention of future rotavirus outbreaks.

VP7 and VP4 genotypes of rotaviruses cocirculating in Iran, 2015 to 2017: Comparison with cogent sequences of Rotarix and RotaTeq vaccine strains before their use for universal mass vaccination

The present study was conducted to analyze the genotypic diversity of circulating species A rotavirus (RVA) strains in Iran and also to investigate comparative analysis between the genotypes of VP4 and VP7 of cocirculating RVA and vaccine strains before the vaccine is introduced in the national immunization program. The G3-lineage I was found in this study as the most common G genotype which was followed by G9-lineage III, G1-lineages I, II, G12-lineage III, G2-lineage IV, and G4-lineage I. Also, P[8]-lineages III, IV was found as the predominant P genotype which was followed by P[4]-lineage V, and P[6]-lineage I. Overally, G3P[8] was determined as the most common combination. Moreover, the analysis of the VP7 antigenic epitopes showed that several amino acid differences existed between circulating Iranian and the vaccine strains. The comparison of genotype G1 of Iranian and vaccine strains (RotaTeq and Rotarix), and genotypes G2, G3, and G4 of Iranian and RotaTeq vaccine strains revealed three to five amino acids differences on the VP7 antigenic epitopes. Furthermore, analyzing of the VP8* epitopes of Iranian P[8] strains indicated that they contained up to 11 and 14 amino acid differences with Rotarix and RotaTeq, respectively. Based on different patterns of amino acid substitutions in circulating and vaccine strains, the emergence of antibody escaping mutants and potentially the decrease of immune protection might ensue in vaccinated children. However, considering the broad cross-protective activity of RVA vaccines, their efficacy should be monitored after the introduction in Iran.

Full genome characterization of human G3P[6] and G3P[9] rotavirus strains in Lebanon

Rotaviruses are the most common infectious agents causing severe diarrheal diseases in young children globally. Three rare human rotavirus strains, two G3P[9] and one G3P[6], were detected in stool samples of children under 5 years of age hospitalized for gastroenteritis in Lebanon during the course of a surveillance study. Complete genomes of these strains were sequenced using VirCapSeq-VERT, a capture based high-throughput sequencing method. Genomic sequences were further characterized by using phylogenetic analyses with global RVA G3P[6]/P[9] strains, other vaccine and reference strains. Genetic analysis revealed that the G3P[6] strain emerged as a DS-1/Wa-like mono-reassortant strain with a potential Ethiopian origin. The two G3P[9] strains possessed a mixed DS-1/Wa/AU-1-like origin indicating that these may have evolved via multiple reassortment events involving feline, human and bovine rotaviruses. Furthermore, analysis of these strains revealed high antigenic variability compared to the vaccine strains. Additional studies are essential to fully understand the evolutionary dynamics of G3P[6]/P[9] strains spreading worldwide and their implications on vaccine effectiveness.

Genotype distribution of Group A rotavirus in children before and after massive vaccination in Latin America and the Caribbean: Systematic review

BACKGROUND

During the last decade, most of Latin American and the Caribbean (LAC) countries have implemented oral live rotavirus vaccines in their national vaccination programs with remarkable results. However, it has been suggested that massive vaccination could lead to the replacement of circulating genotypes or the emergence of new variants or neutralizing antibodies escape mutants, which may reduce the effectiveness of the vaccine. The objective was to analyze the genetic diversity of Group A rotavirus before and after the introduction of universal vaccination in LAC.

METHODS

We conducted a systematic review of studies published in PubMed, Scielo and LILACS. There were considered only LAC countries with rotavirus massive vaccination strategy which had described circulating genotypes data in children under 5 years of age, either for surveillance or vaccine effectiveness purposes, from 2001 to 2017. Systematic review stages were carried out following the recommendations of PRISMA.

RESULTS

Of the 18 countries that included any of the two licensed rotavirus vaccines in their national schedules since 2006, only 7 (~39%) presented studies of RVA genetic diversity before and after implementation, and met the inclusion criteria. Four of them (Argentina, Brazil, Colombia and Nicaragua) experienced a rapid switch from Wa-like to DS-1-like strains. Also, G1P[8] association, considered the most predominant worldwide in the pre-vaccination era, decreased significantly and was only frequently detected in Venezuela and Nicaragua. No defined pattern of emergence at high frequencies of unusual associations was observed in the post vaccination period, except for some evidence of G9P[4] in Colombia, G3P[6] and G1P[4] in Nicaragua.

CONCLUSIONS

Even though the evidence shows a DS-1-like change trend, data from studies conducted in Latin America and the Caribbean are diverse and still not sufficient to assess the impact of vaccines on viral ecology or if genetic diversity is influenced by natural mechanisms of fluctuation.

Whole Genome Analysis of Selected Human Group A Rotavirus Strains Revealed Evolution of DS-1-Like Single- and Double-Gene Reassortant Rotavirus Strains in Pakistan During 2015-2016

Acute gastroenteritis due to group A rotaviruses (RVAs) is the leading cause of infant and childhood morbidity and mortality particularly in developing countries including Pakistan. In this study we have characterized the whole genomes of five RVA strains (PAK56, PAK419, PAK585, PAK622, and PAK663) using the Illumina HiSeq platform. The strains PAK56 and PAK622 exhibited a typical Wa-like genotype constellation (G9-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1 and G3-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1, respectively), whereas PAK419, PAK585, and PAK663 exhibited distinct DS-1-like genotype constellations (G3P[4]-I2-R2-C2-M2-A2-N2-T1-E2-H2, G1P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2, and G3P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2, respectively). Despite their DS-1-like genotype constellation, strain PAK585 possessed the typical Wa-like G1P[8] genotypes, whereas both PAK419 and PAK663 possessed the G3 genotype. In addition, PAK419 also possessed the Wa-like NSP3 genotype T1, suggesting that multiple reassortments have occurred. On Phylogenetic analysis, all of the gene segments of the five strains examined in this study were genetically related to globally circulating human G1, G2, G3, G6, G8, G9, and G12 strains. Interestingly, the NSP2 gene of strain PAK419 showed closest relationship with Indian bovine strain (India/HR/B91), suggesting the occurrence of reassortment between human and bovine RVA strains. Furthermore, strains PAK419, PAK585, and PAK663 were closely related to one another in most of their gene segments, indicating that these strains might have been derived from a common ancestor. To our knowledge this is the first whole genome-based molecular characterization of human rotavirus strains in Pakistan. The results of our study will enhance our existing knowledge on the diversity and evolutionary dynamics of novel RVA strains including DS-1-like intergenogroup reassortant strains spreading in Asian countries including Pakistan, in the pre-vaccine era. Therefore, continuous surveillance is recommended to monitor the evolution, spread and genetic stability of novel reassortant rotavirus strains derived from such events.

Emergence of Rare Bovine-Human Reassortant DS-1-Like Rotavirus A Strains with G8P[8] Genotype in Human Patients in the Czech Republic

Group A Rotaviruses (RVA) are the leading cause of acute gastroenteritis in children and a major cause of childhood mortality in low-income countries. RVAs are mostly host-specific, but interspecies transmission and reassortment between human and animal RVAs significantly contribute to their genetic diversity. We investigated the VP7 and VP4 genotypes of RVA isolated from 225 stool specimens collected from Czech patients with gastroenteritis during 2016–2019. The most abundant genotypes were G1P[8] (42.7%), G3P[8] (11.1%), G9P[8] (9.8%), G2P[4] (4.4%), G4P[8] (1.3%), G12P[8] (1.3%), and, surprisingly, G8P[8] (9.3%). Sequence analysis of G8P[8] strains revealed the highest nucleotide similarity of all Czech G8 sequences to the G8P[8] rotavirus strains that were isolated in Vietnam in 2014/2015. The whole-genome backbone of the Czech G8 strains was determined with the use of next-generation sequencing as DS-1-like. Phylogenetic analysis of all segments clustered the Czech isolates with RVA strains that were formerly described in Southeast Asia, which had emerged following genetic reassortment between bovine and human RVAs. This is the first time that bovine–human DS-1-like G8P[8] strains were detected at a high rate in human patients in Central Europe. Whether the emergence of this unusual genotype reflects the establishment of a new RVA strain in the population requires the continuous monitoring of rotavirus epidemiology.

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.

Development of a Colloidal Gold Immunochromatographic Strip Assay for Rapid Detection of Bovine Rotavirus

Bovine rotavirus (BRV) is one of main pathogens responsible for diarrhea, fever, and vomiting. In this study, we developed a colloidal gold immunochromatographic test strip for detecting BRV according to the principle of double-antibody sandwich. The monoclonal antibodies (mAbs) and polyclonal antibodies (pAbs) were prepared and purified. On the strip, the purified mAbs labeled with the colloidal gold were used as the detector, and the goat anti-mouse antibodies and purified pAbs were coated on the nitrocellulose membranes as the control line and the test line, respectively. We optimized different reaction conditions, including the amount of mAbs, the pH of colloidal gold solution, coating solution, blocking solution, sample pad treatment solution, antibody concentration in control line, and antibody concentration in detection line. In specificity assay, the strip had high specificity in detecting BRV. No cross-reaction was observed in detecting other viruses. The detection sensitivity of the strip was found to be 1 × 10³ TCID₅₀/0.1 mL. Two hundred twenty clinical samples were detected with the strip compared to reverse transcription-polymerase chain reaction. No false-negative or false-positive results were found, and the results obtained by the two methods were similar. In conclusion, we developed a novel immunochromatographic strip to rapidly detect BRV. The strip developed exhibited high sensitivity and specificity for BRV detection. It could be a rapid, convenient, and effective method for the rapid diagnosis of BRV infection in the fields.

Molecular Characteristics of Novel Mono-Reassortant G9P[8] Rotavirus A Strains Possessing the NSP4 Gene of the E2 Genotype Detected in Tokyo, Japan

Rotavirus A (RVA) has been detected in patients with gastroenteritis even after vaccine introduction in Japan. To investigate circulating RVA strains, RVA-positive stool specimens obtained in Tokyo in 2017 and 2018 were analyzed using next-generation sequencing. A total of 50 and 21 RVA samples were obtained in 2017 and 2018, respectively. In 2017, G2P[4] (40.0%) was the most prevalent strain, followed by G3P[8] (DS-1-like) (28.0%), G8P[8] (10.0%), G3P[8] (Wa-like) (8.0%), G9P[8]-E1 (8.0%), and mixed infection (6.0%). In 2018, G3P[8] (DS-1-like) (28.6%) and G9P[8]-E2 (28.6%) were the most prevalent strains, followed by G9P[8]-E1 (19.0%), G2P[4] (9.5%), G8P[8] (9.5%), and mixed infection (4.8%). Six G9P[8]-E2 strains detected in 2018 showed an atypical genotype constellation (G9P[8]-I1-R1-C1-M1-A1-N1-T1-E2-H1), which had not been reported previously. Phylogenetic analyses suggested that the RVA virus was generated by inter-genogroup reassortment between commonly circulating G9P[8] and G2P[4] strains in Japan. The G9P[8] strain seemed to be reassorted with only the NSP4 gene of the E2 genotype of the G2P[4] strain. Since this newly-emerged G9P[8]-E2 virus was detected in different locations in Tokyo, the virus appears to have already begun to spread to a wider area.

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.

Dual Recognition of Sialic Acid and αGal Epitopes by the VP8* Domains of the Bovine Rotavirus G6P[5] WC3 and of Its Mono-reassortant G4P[5] RotaTeq Vaccine Strains

Group A rotaviruses initiate infection through the binding of the VP8* domain of the VP4 protein to sialic acids (SAs) or histo-blood group antigens (HBGAs). Although the bovine G6P[5] WC3 strain is an important animal pathogen and is used as the backbone in the bovine-human reassortant RotaTeq vaccine, the receptor(s) for their P[5] VP8* domain has remained elusive. Using a variety of approaches, we demonstrated that the WC3 and bovine-human mono-reassortant G4P[5] vaccine strains recognize both α2,6-linked SA and αGal HBGA as ligands. Neither ligand is expressed on human small intestinal epithelial cells, explaining the absence of natural human infection by P[5]-bearing strains. However, we observed that the P[5]-bearing WC3 and G4P[5] RotaTeq vaccine strains could still infect human intestinal epithelial cells. Thus, the four P[5] RotaTeq vaccine strains potentially binding to additional alternative receptors may be efficient and effective in providing protection against severe rotavirus disease in human.

Group A rotaviruses, an important cause of severe diarrhea in children and young animals, initiate infection via interactions of the VP8* domain of the VP4 spike protein with cell surface sialic acids (SAs) or histo-blood group antigens (HBGAs). Although the bovine G6P[5] WC3 strain is an important animal pathogen and is also used in the bovine-human reassortant RotaTeq vaccine, the receptor(s) for the VP8* domain of WC3 and its reassortant strains have not yet been identified. In the present study, HBGA- and saliva-binding assays showed that both G6P[5] WC3 and mono-reassortant G4P[5] strains recognized the αGal HBGA. The infectivity of both P[5]-bearing strains was significantly reduced in αGal-free MA-104 cells by pretreatment with a broadly specific neuraminidase or by coincubation with the α2,6-linked SA-specific Sambucus nigra lectin, but not by the α2,3-linked specific sialidase or by Maackia amurensis lectin. Free NeuAc and the αGal trisaccharide also prevented the infectivity of both strains. This indicated that both P[5]-bearing strains utilize α2,6-linked SA as a ligand on MA104 cells. However, the two strains replicated in differentiated bovine small intestinal enteroids and in their human counterparts that lack α2,6-linked SA or αGal HBGA, suggesting that additional or alternative receptors such as integrins, hsp70, and tight-junction proteins bound directly to the VP5* domain can be used by the P[5]-bearing strains to initiate the infection of human cells. In addition, these data also suggested that P[5]-bearing strains have potential for cross-species transmission.

IMPORTANCE Group A rotaviruses initiate infection through the binding of the VP8* domain of the VP4 protein to sialic acids (SAs) or histo-blood group antigens (HBGAs). Although the bovine G6P[5] WC3 strain is an important animal pathogen and is used as the backbone in the bovine-human reassortant RotaTeq vaccine, the receptor(s) for their P[5] VP8* domain has remained elusive. Using a variety of approaches, we demonstrated that the WC3 and bovine-human mono-reassortant G4P[5] vaccine strains recognize both α2,6-linked SA and αGal HBGA as ligands. Neither ligand is expressed on human small intestinal epithelial cells, explaining the absence of natural human infection by P[5]-bearing strains. However, we observed that the P[5]-bearing WC3 and G4P[5] RotaTeq vaccine strains could still infect human intestinal epithelial cells. Thus, the four P[5] RotaTeq vaccine strains potentially binding to additional alternative receptors may be efficient and effective in providing protection against severe rotavirus disease in human.

Natural selection contributed to immunological differences between hunter-gatherers and agriculturalists

The shift from a hunter-gatherer to an agricultural mode of subsistence is believed to have been associated with profound changes in the burden and diversity of pathogens across human populations. Yet, the extent to which the advent of agriculture affected the evolution of the human immune system remains unknown. Here we present a comparative study of variation in the transcriptional responses of peripheral blood mononuclear cells to bacterial and viral stimuli between Batwa rainforest hunter-gatherers and Bakiga agriculturalists from Uganda. We observed increased divergence between hunter-gatherers and agriculturalists in the early transcriptional response to viruses compared with that for bacterial stimuli. We demonstrate that a significant fraction of these transcriptional differences are under genetic control and we show that positive natural selection has helped to shape population differences in immune regulation. Across the set of genetic variants underlying inter-population immune-response differences, however, the signatures of positive selection were disproportionately observed in the rainforest hunter-gatherers. This result is counter to expectations on the basis of the popularized notion that shifts in pathogen exposure due to the advent of agriculture imposed radically heightened selective pressures in agriculturalist populations.

Molecular characterization of human group A rotavirus genotypes circulating in Rawalpindi, Islamabad, Pakistan during 2015-2016

Group A rotaviruses (RVA) are one of the major causes of acute gastroenteritis (AGE) in young children worldwide. Owing to lack of proper surveillance programs and health facilities, developing countries of Asia and Africa carry a disproportionately heavy share of the RVA disease burden. The aim of this hospital-based study was to investigate the circulation of RVA genotypes in Rawalpindi and Islamabad, Pakistan in 2015 and 2016, prior to the implementation of RVA vaccine. 639 faecal samples collected from children under 10 years of age hospitalized with AGE were tested for RVA antigen by ELISA. Among 171 ELISA positive samples, 143 were successfully screened for RT-PCR and sequencing. The prevalence of RVA was found to be 26.8% with the highest frequency (34.9%) found among children of age group 6-11 months. The most predominant circulating genotypes were G3P[8] (22.4%) followed by G12P[6] (20.3%), G2P[4] (12.6%), G1P[8] (11.9%), G9P[6] (11.9%), G3P[4] (9.1%), G1P[6] (4.2%), G9P[8] (4.2%), and G3P[6] (0.7%). A single mixed genotype G1G3P[8] was also detected. The findings of this study provide baseline data, that will help to assess if future vaccination campaigns using currently available RVA vaccine will reduce RVA disease burden and instigate evolutionary changes in the overall RVA biology. The high prevalence of RVA infections in Pakistan require to improve and strengthen the surveillance and monitoring system for RVA. This will provide useful information for health authorities in planning public health care strategies to mitigate the disease burden caused by RVA.

Whole-genome analysis of rotavirus G4P[6] strains isolated from Korean neonates: association of Korean neonates and rotavirus P[6] genotypes

Background

Group A rotaviruses are the major causative agents of pediatric gastroenteritis worldwide. Several studies have reported the predominance of G4P[6] rotavirus genotypes in Korean neonates, which is uncommon in other countries. Therefore, the purposes of this study were to determine the genotype constellations of complete genomes of G4P[6] rotavirus strains isolated from Korean neonates using next-generation sequencing, to compare these sequences with other G4P[6] strains in other countries, and to determine the reason for the predominance of G4P[6] genotypes in Korean neonates.

Results

Twenty rotavirus G4P[6] strains, isolated from January 2013 to January 2016, were selected for whole-genome sequencing. Eleven rotavirus genes were amplified using specific primer sets, and sequencing was carried out using an Ion S5 XL next-generation sequencing platform. Genotypes of each gene were determined, and phylogenetic analyses were performed to investigate genetic distances between genes of rotaviruses in this study and those of other rotavirus G4P[6] strains whose whole-genome sequences were previously published. All 20 rotavirus strains in this study had the same genotype: G4-P[6]-I1-R1-C1-M1-A1-N1-T1-E1-H1, representing the Wa-like genotype constellation. BLAST searches of 20 G4P[6] rotavirus strains revealed that all G4 sequences in this study showed the highest nucleotide identity to G4 sequences of G4P[6] rotavirus strains isolated in Korea in 2008 (GenBank accession number: FJ603447). Additionally, P[6] gene sequences in this study showed the highest nucleotide identity to P[6] sequences of G4P[6] strains detected in Korea in 2002 (AY158093). Phylogenetic and nucleotide sequence analyses showed that G4P[6] strains in this study and previously reported G4P[6] strains in Korea were mostly detected in neonates and had similar G4 and P[6] sequences compared with other G4P[6] strains detected in other countries.

Conclusions

This study showed that the whole-genome constellation of rotavirus G4P[6] strains from Korean neonates resembled a Wa-like genotype constellation. Additionally, rotavirus genotypes detected in Korean neonates had unique P[6] sequences, which may be the cause of Korean neonatal rotavirus infection.

Electronic supplementary material

The online version of this article (10.1186/s13099-019-0318-5) contains supplementary material, which is available to authorized users.

Whole genome characterization and evolutionary analysis of OP354-like P[8] Rotavirus A strains isolated from Ghanaian children with diarrhoea

In 2010, the rare OP354-like P[8]b rotavirus subtype was detected in children less than 2 years old in Ghana. In this follow-up study, to provide insight into the evolutionary history of the genome of Ghanaian P[8]b strains RVA/Human-wt/GHA/GHDC949/2010/G9P[8] and RVA/Human-wt/GHA/GHM0094/2010/G9P[8] detected in an infant and a 7-month old child hospitalised for acute gastroenteritis, we sequenced the complete genome using both Sanger sequencing and Illumina MiSeq technology followed by phylogenetic analysis of the near-full length sequences. Both strains possessed the Wa-like/genotype 1 constellation G9P[8]b-I1-R1-C1-M1-A1-N1-T1-E1-H1. Sequence comparison and phylogenetic inference showed that both strains were identical at the lineage level throughout the 11 genome segments. Their VP7 sequences belonged to the major sub-lineage of the G9-lineage III whereas their VP4 sequences belonged to P[8]b cluster I. The VP7 and VP4 genes of the study strains were closely related to a Senegalese G9P[8]b strain detected in 2009. In the remaining nine genome segments, both strains consistently clustered together with Wa-like RVA strains possessing either P[8]a or P[8]b mostly of African RVA origin. The introduction of a P[8]b subtype VP4 gene into the stable Wa-like strain backbone may result in strains that might propagate easily in the human population, with a potential to become an important public health concern, especially because it is not certain if the monovalent rotavirus vaccine (Rotarix) used in Ghana will be efficacious against such strains. Our analysis of the full genomes of GHM0094 and GHDC949 adds to knowledge of the genetic make-up and evolutionary dynamics of P[8]b rotavirus strains.

Full-length genome analysis of the first human G8P[14] rotavirus strain from Morocco suggests evidence of zoonotic transmission

An unusual group A rotavirus (RVA) strain MAR/ma31/2011/G8P[14] was detected for the first time in Morocco in a stool sample from hospitalized child aged 18 months suffering from acute gastroenteritis and fever in 2011. Complete genome sequencing of the ma31 strain was done using the capillary sequencing technology. The analysis revealed the G8-P[14]-I2-R2-C2-M2-A11-N2-T6-E2-H3 constellation and the backbone genes: I2-R2-C2-M2-A11-N2-T6-E2-H3 are commonly found in RVA strains from artiodactyls such as cattle. The constellation was shared with another Italian zoonotic G8P[14] strains (BA01 and BA02), two Hungarian human strains (182-02 and BP1062) and a sheep RVA strain OVR762. Phylogenetic analysis of each genome segment of ma31 revealed a mixed gene configuration originated from animals and human. Comparison of the antigenic regions of VP7 and VP4 amino acid sequences between ma31 strain and selected animal and human strains bearing G8 and or P[14], showed a high level of conservation, while many substitutions was observed in comparison with RotaTeq™ and Rotarix™ vaccine strains. In contrast, alignment analysis of the four antigenic sites of VP6 revealed a high degree of conservation. These findings reveal a typical zoonotic origin of the strain and confirm a high potential for RVA zoonotic transmission between bovine and humans, allowing the generation of novel rotavirus genotypes.

Whole-gene analysis of inter-genogroup reassortant rotaviruses from the Dominican Republic: Emergence of equine-like G3 strains and evidence of their reassortment with locally-circulating strains

Inter-genogroup reassortant group A rotavirus (RVA) strains possessing a G3 VP7 gene of putative equine origin (EQL-G3) have been detected in humans since 2013. Here we report detection of EQL-G3P[8] RVA strains from the Dominican Republic collected in 2014-16. Whole-gene analysis of RVA in stool specimens revealed 16 EQL-G3P[8] strains, 3 of which appear to have acquired an N1 NSP1 gene from locally-circulating G9P[8] strains and a novel G2P[8] reassortant possessing 7 EQL-G3-associated genes and 3 genes from a locally-circulating G2P[4] strain. Phylogenetic/genetic analyses of VP7 gene sequences revealed nine G3 lineages (I-IX) with newly-assigned lineage IX encompassing all reported human EQL-G3 strains along with the ancestral equine strain. VP1 and NSP2 gene phylogenies suggest that EQL-G3P[8] strains were introduced into the Dominican Republic from Thailand. The emergence of EQL-G3P[8] strains in the Dominican Republic and their reassortment with locally-circulating RVA could have implications for current vaccination strategies.

Sub-genotype phylogeny of the non-G, non-P genes of genotype 2 Rotavirus A strains

Recent increase in the detection of unusual G1P[8], G3P[8], G8P[8], and G9P[4] Rotavirus A (RVA) strains bearing the DS-1-like constellation of the non-G, non-P genes (hereafter referred to as the genotype 2 backbone) requires better understanding of their evolutionary relationship. However, within a genotype, there is lack of a consensus lineage designation framework and a set of common sequences that can serve as references. Phylogenetic analyses were carried out on over 8,500 RVA genotype 2 genes systematically retrieved from the rotavirus database within the NCBI Virus Variation Resource. In line with previous designations, using pairwise comparison of cogent nucleotide sequences and stringent bootstrap support, reference lineages were defined. This study proposes a lineage framework and provides a dataset ranging from 34 to 145 sequences for each genotype 2 gene for orderly lineage designation of global genotype 2 genes of RVAs detected in human and animals. The framework identified five to 31 lineages depending on the gene. The least number of lineages (five to seven) were observed in genotypes A2 (NSP1), T2 (NSP3) and H2 (NSP5) which are limited to human RVA whereas the most number of lineages (31) was observed in genotype E2 (NSP4). Sharing of the same lineage constellations of the genotype 2 backbone genes between recently-emerging, unusual G1P[8], G3P[8], G8P[8] and G9P[4] reassortants and many contemporary G2P[4] strains provided strong support to the hypothesis that unusual genotype 2 strains originated primarily from reassortment events in the recent past involving contemporary G2P[4] strains as one parent and ordinary genotype 1 strains or animal RVA strains as the other. The lineage framework with selected reference sequences will help researchers to identify the lineage to which a given genotype 2 strain belongs, and trace the evolutionary history of common and unusual genotype 2 strains in circulation.

Whole-genome characterization of G12 rotavirus strains detected in Mozambique reveals a co-infection with a GXP[14] strain of possible animal origin

A high prevalence of G12 rotavirus strains has previously been reported in southern Mozambique. In this study, the full genomes of five Mozambican G12 strains were determined directly from stool using an Illumina Miseq platform. One sample (0060) contained an intergenogroup co-infection of a G12P[8] Wa-like strain and a GXP[14] DS-1-like strain. The sequences of seven genome segments, detected for the GXP[14] strain, clustered with a diverse group of mostly animal strains, suggesting co-infection with a strain of possible animal origin. The stool samples contained G12P[6] rotavirus strains with Wa-like backbones. Phylogenetic analyses of the VP4 and VP7 encoding segments of the G12P[6] strains suggested that they were reassortants containing backbones that are similar to that of the G12P[8] strain. The study confirms previous observations of interspecies transmission and emphasizes the importance of whole-genome sequencing in order to evaluate rotavirus co-infections and reassortants.

Evidence of zoonotic transmission of VP6 and NSP4 genes into human species A rotaviruses isolated in Pakistan in 2010

Introduction of animal group A rotavirus (RVA) gene segments into the human RVA population is a major factor shaping the genetic landscape of human RVA strains. The VP6 and NSP4 genes of 74 G/P-genotyped RVA isolates collected in Rawalpindi during 2010 were analyzed, revealing the presence of VP6 genotypes I1 (60.8%) and I2 (39.2%) and NSP4 genotypes E1 (60.8%), E2 (28.3%) and E-untypable (10.8%) among the circulating human RVA strains. The typical human RVA combinations I1E1 and I2E2 were found in 59.4% and 24.3% of the cases, respectively, whereas 5.4% of the RVA strains were reassortants, i.e., either I1E2 or I2E1. The phylogeny of the NSP4 gene showed that one G2P[4] and two G1P[6] RVA strains clustered with porcine E1 RVA strains or RVA strains that were considered to be (partially) of porcine origin. In addition, the NSP4 gene segment of the unusual human G6P[1] RVA strains clustered closely with bovine E2 RVA strains, further strengthening the hypothesis of an interspecies transmission event. The study further demonstrates the role of genomic re-assortment and the involvement of interspecies transmission in the evolution of human RVA strains. The VP6 and NSP4 nucleotide sequences analyzed in the study received the GenBank accession numbers KC846908- KC846971 and KC846972-KC847037, respectively.

Molecular Epidemiology and Clinical Features of Rotavirus Infection Among Pediatric Patients in East Java, Indonesia During 2015-2018: Dynamic Changes in Rotavirus Genotypes From Equine-Like G3 to Typical Human G1/G3

Group A rotavirus (RVA) is the most important cause of severe gastroenteritis among children worldwide, and effective RVA vaccines have been introduced in many countries. Here we performed a molecular epidemiological analysis of RVA infection among pediatric patients in East Java, Indonesia, during 2015-2018. A total of 432 stool samples were collected from hospitalized pediatric patients with acute gastroenteritis. None of the patients in this cohort had been immunized with an RVA vaccine. The overall prevalence of RVA infection was 31.7% (137/432), and RVA infection was significantly more prevalent in the 6- to 11-month age group than in the other age groups (P < 0.05). Multiplex reverse transcription-PCR (RT-PCR) revealed that the most common G-P combination was equine-like G3P[8] (70.8%), followed by equine-like G3P[6] (12.4%), human G1P[8] (8.8%), G3P[6] (1.5%), and G1P[6] (0.7%). Interestingly, the equine-like strains were exclusively detected until May 2017, but in July 2017 they were completely replaced by a typical human genotype (G1 and G3), suggesting that the dynamic changes in RVA genotypes from equine-like G3 to typical human G1/G3 in Indonesia can occur even in the country with low RVA vaccine coverage rate. The mechanism of the dynamic changes in RVA genotypes needs to be explored. Infants and children with RVA-associated gastroenteritis presented more frequently with some dehydration, vomiting, and watery diarrhea, indicating a greater severity of RVA infection compared to those with non-RVA gastroenteritis. In conclusion, a dynamic change was found in the RVA genotype from equine-like G3 to a typical human genotype. Since severe cases of RVA infection were prevalent, especially in children aged 6 to 11 months or more generally in those less than 2 years old, RVA vaccination should be included in Indonesia's national immunization program.

Development and evaluation of a one-step multiplex real-time TaqMan® RT-qPCR assay for the detection and genotyping of equine G3 and G14 rotaviruses in fecal samples

BACKGROUND

Equine rotavirus A (ERVA) is the leading cause of diarrhea in neonatal foals and has a negative impact on equine breeding enterprises worldwide. Among ERVA strains infecting foals, the genotypes G3P[12] and G14P[12] are the most prevalent, while infections by strains with other genomic arrangements are infrequent. The identification of circulating strains of ERVA is critical for diagnostic and surveillance purposes, as well as to understand their molecular epidemiology. Current genotyping methods available for ERVA and rotaviruses affecting other animal species rely on Sanger sequencing and are significantly time-consuming, costly and labor intensive. Here, we developed the first one-step multiplex TaqMan® real-time reverse transcription polymerase chain reaction (RT-qPCR) assay targeting the NSP3 and VP7 genes of ERVA G3 and G14 genotypes for the rapid detection and G-typing directly from fecal specimens.

METHODS

A one-step multiplex TaqMan® RT-qPCR assay targeting the NSP3 and VP7 genes of ERVA G3 and G14 genotypes was designed. The analytical sensitivity was assessed using serial dilutions of in vitro transcribed RNA containing the target sequences while the analytical specificity was determined using RNA and DNA derived from a panel of group A rotaviruses along with other equine viruses and bacteria. The clinical performance of this multiplex assay was evaluated using a panel of 177 fecal samples and compared to a VP7-specific standard RT-PCR assay and Sanger sequencing. Limits of detection (LOD), sensitivity, specificity, and agreement were determined.

RESULTS

The multiplex G3 and G14 VP7 assays demonstrated high specificity and efficiency, with perfect linearity. A 100-fold difference in their analytical sensitivity was observed when compared to the singleplex assays; however, this difference did not have an impact on the clinical performance. Clinical performance of the multiplex RT-qPCR assay demonstrated that this assay had a high sensitivity/specificity for every target (100% for NSP3, > 90% for G3 VP7 and > 99% for G14 VP7, respectively) and high overall agreement (> 98%) compared to conventional RT-PCR and sequencing.

CONCLUSIONS

This new multiplex RT-qPCR assay constitutes a useful, very reliable tool that could significantly aid in the rapid detection and G-typing of ERVA strains circulating in the field.

Whole genome characterization of feline-like G3P[8] reassortant rotavirus A strains bearing the DS-1-like backbone genes detected in Vietnam, 2016

While conducting rotavirus gastroenteritis surveillance in Vietnam, two G3P[8] rotavirus A specimens possessing an identical short RNA electropherotype were detected. They were RVA/Human-wt/VNM/0232/2016/G3P[8] and RVA/Human-wt/VNM/0248/2016/G3P[8], and recovered from 9 and 23 months old boys, respectively. The patients developed diarrhoea within one-week interval in March 2016 but in places >100 km apart in northern Vietnam. Whole genome sequencing of the two G3P[8] rotavirus A strains revealed that their genomic RNA sequences were identical across the 11 genome segments, suggesting that they derived from a single clone. The backbone gene constellation was I2-R2-C2-M2-A2-N2-T2-E2-H2. The backbone genes and the VP4 gene had a virtually identical nucleotide sequences with identities ranging from 99.2 to 100% to the corresponding genes of RVA/Human-wt/VNM/1149/2014/G8P[8]; the prototype of recently-emerging bovine-like G8P[8] reassortant strains in Vietnam. On the other hand, the VP7 gene was 98.8% identical with that of RVA/Human-wt/CHN/E2451/2011/G3P[9], and they were clustered together in the lineage represented by RVA/Cat-tc/JPN/FRV-1/1986/G3P[9]. The observations led us to hypothesise that one of the bovine-like G8P[8] strains bearing the DS-1-like backbone genes reassorted with a locally circulating FRV-1-like strain to gain the G3 VP7 gene and to emerge as a thus-far undescribed feline-like G3P[8] reassortant strain. The identification of feline-like G3P[8] strains bearing the DS-1-like backbone genes exemplifies the strength and necessity of the whole genome sequencing approach in monitoring, describing and understanding the evolutionary changes that are occurring in emerging strains and their interactions with co-circulating strains.

Molecular chracterisation of porcine group A rotaviruses: Studies on the age-related occurrence and spatial distribution of circulating virus genotypes in Poland

Rotaviruses of group A (RVAs) commonly occur in farm animals. In pigs, they cause acute gastrointestinal disease which is considered as significant factor of economic losses in pig farming. The aim of the study was an assessment of the prevalence of rotavirus (RV) infections in farmed pigs in Poland, genotype identification of the virus strains in conjunction with their age-related occurrence and regional (province) distribution pattern in pig herds. In total, 920 pig faecal samples were collected from pigs between the ages of one week and two years old from 131 farms. RVAs were detected using ELISA and molecular methods followed by a sequence-based identification of G (VP7) and P (VP4) virus genotypes. RV antigen was found in 377 (41%) of pig faecal samples. The correlation between pig age and frequency of RV infections was shown. In the Polish pig population, 145 RVA strains representing 33 GP genotypes were identified. Subsequent molecular analysis revealed an age-dependent and regional diversity in distribution of genotypes and virus strains. Besides typical pig RVA strains, novel strains such as G5P [34], G9P[34], and human G1P[8] were identified in this animal host. Findings from this study showed a change over time in the genotype occurrence of circulating pig RVAs in Poland. The high genetic variability of RV strains and acquisition of new virus genotypes have led to the emergence of novel, genetically distinct RVAs. The changes in the genotype occurrence of RVA strains in pigs indicate the need for their continuous epidemiological surveillance.

Improvement of Rotavirus Genotyping Method by Using the Semi-Nested Multiplex-PCR With New Primer Set

Rotavirus A (RVA) is a major cause of gastroenteritis in infants and young children. After vaccine introduction, RVA surveillance has become more important for monitoring changes in genotype distribution, and the semi-nested multiplex-PCR is a popular method for RVA genotyping. In particular, the VP7 primer set reported by Gouvea and colleagues in 1990 is still widely used worldwide as the recommended WHO primer set in regional and national reference RVA surveillance laboratories. However, this primer set yielded some mistakes with recent epidemic strains. The newly emerged equine-like G3 strains were mistyped as G1, G8 strains were mistyped as G3, the G9 lineage 3 strains showed very weak band, and the G9 lineage 6 strains showed a G9-specific band and a non-specific band. Gouvea’s standard protocol has become relatively unreliable for identifying genotypes correctly. To overcome this limitation, we redesigned the primer set to include recent epidemic strains. Our new primer set enabled us to correctly identify the VP7 genotypes of representative epidemic strains by agarose gel electrophoresis (G1, G2, human typical G3, equine-like G3, G4, G8, G9, and G12). We believe that the multiplex-PCR method with our new primer set is a useful and valuable tool for surveillance of RVA epidemics.

Impact of rotavirus and hepatitis A virus by worldwide climatic changes during the period between 2000 and 2013

Enteric viruses are present in the environment as a result of the discharge of poorly or untreated wastewater. The spread of enteric viruses in the environment depend to human activities like stools of infected individuals ejected in the external environment can be transmitted by water sources and back to susceptible individuals for other cycles of illness. Among the enteric viruses Rotaviruses (RV) and Hepatitis A viruses (HAV) is the most detected in wastewater causing gastroenteritis and acute hepatitis. Therefore, it is of interest to climate change, mainly temperature and carbon Dioxide (CO2) variations, on Rotavirus and Hepatitis A as a model of enteric viruses present in the aquatic environment using computational modelling tools. The results of genetic ratio showed a negative correlation between the epidemiological data and the mutation rate. However, the correlation was positive between the temperature, CO2 increase, and the rate of mutation. The positive correlation is explained by the adaptation of the viruses to the climatic changes, the RNA polymerase of the RV induces errors to adapt to the environmental conditions. The simultaneous increase in number of infections and temperature in 2010 has been demonstrated in previous studies deducing that viral pathogenicity increase with temperature increase.

Characterization of an Unusual DS-1-Like G8P[8] Rotavirus Strain from Japan in 2017: Evolution of Emerging DS-1-Like G8P[8] Strains through Reassortment

The emergence of unusual DS-1-like intergenogroup reassortant rotaviruses with a bovine-like G8 genotype (DS-1-like G8P[8] strains) has been reported in several Asian countries. During the rotavirus surveillance program in Japan in 2017, a DS-1-like G8P[8] strain (RVA/Human-wt/JPN/SO1162/2017/G8P[8]) was identified in 43 rotavirus-positive stool samples. Strain SO1162 was shown to have a unique genotype constellation, including genes from both genogroup 1 and 2: G8-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Phylogenetic analysis revealed that the VP1 gene of strain SO1162 appeared to have originated from DS-1-like G1P[8] strains from Thailand and Vietnam, while the remaining 10 genes were closely related to those of previously reported DS-1-like G8P[8] strains. Thus, SO1162 was suggested to be a reassortant strain that acquired the VP1 gene from Southeast Asian DS-1-like G1P[8] strains on the genetic background of co-circulating DS-1-like G8P[8] strains. Our findings provide important insights into the evolutionary dynamics of emerging DS-1-like G8P[8] strains.

Genome analysis of a G9P[23] group A rotavirus isolated from a dog with diarrhea in China

Genotype G9 is an emerging genotype among species A rotavirus (RVA) circulating in humans and pigs worldwide. In this study, an RVA strain designated RVA/Dog-tc/CHN/SCCD-A/2017/G9P[23] was isolated in cell culture from a pet dog stool sample with acute diarrhea, and its whole genome was sequenced. The genotype constellation of SCCD-A was G9-P[23]-I5-R1-C1-M1-A8-N1-T1-E1-H1. All genome segments except the VP1 gene were closely related to the genes from porcine RVA strains or porcine-like human RVA strains. On the other hand, the VP1 gene clustered in a distinct lineage only with that of a G5P[6] porcine-like human RVA, preventing the identification of the exact host species origin, but very unlikely to be originated from human RVA. In addition, phylogenetic analysis showed that the G9 VP7 gene of SCCD-A clustered into a novel sublineage within the lineage III of G9. This first isolation of a G9P[23] RVA from a pet dog may justify the exploration of the role dogs play in the interaction of RVA circulating in pigs and humans.

•

First identified G9P[23] group A rotavirus from dog and the genome of RVA/Dog-tc/CHN/SCCD-A/2017/G9P[23]was determined.

•

The strain’s genotype constellation as G9-P[23]-I5-R1-C1-M1-A8-N1-T1-E1-H1.

•

G9P[23] from dog may justify the exploration of the role dogs play in the interaction of RVA circulating in pigs and humans.

Rotavirus Infection Alters Splicing of the Stress-Related Transcription Factor XBP1

XBP1 is a stress-regulated transcription factor also involved in mammalian host defenses and innate immune response. Our investigation of XBP1 RNA splicing during rotavirus infection revealed that an additional XBP1 RNA (XBP1es) that corresponded to exon skipping in the XBP1 pre-RNA is induced depending on the rotavirus strain used. We show that the translation product of XBP1es (XBP1es) has trans-activation properties similar to those of XBP1 on ER stress response element (ERSE) containing promoters. Using monoreassortant between ES⁺ ("skipping") and ES^- ("nonskipping") strains of rotavirus, we show that gene 7 encoding the viral translation enhancer NSP3 is involved in this phenomenon and that exon skipping parallels the nuclear relocalization of cytoplasmic PABP. We further show, using recombinant rotaviruses carrying chimeric gene 7, that the ES⁺ phenotype is linked to the eIF4G-binding domain of NSP3. Because the XBP1 transcription factor is involved in stress and immunological responses, our results suggest an alternative way to activate XBP1 upon viral infection or nuclear localization of PABP.IMPORTANCE Rotavirus is one of the most important pathogens causing severe gastroenteritis in young children worldwide. Here we show that infection with several rotavirus strains induces an alternative splicing of the RNA encoding the stressed-induced transcription factor XBP1. The genetic determinant of XBP1 splicing is the viral RNA translation enhancer NSP3. Since XBP1 is involved in cellular stress and immune responses and since the XBP1 protein made from the alternatively spliced RNA is an active transcription factor, our observations raise the question of whether alternative splicing is a cellular response to rotavirus infection.

Study of Complete Genome Sequences of Rotavirus A Epidemics and Evolution in Japan in 2012-2014

A comprehensive molecular epidemiological study using next-generation sequencing technology was conducted on 333 rotavirus A (RVA)-positive specimens collected from six sentinel hospitals across Japan over three consecutive seasons (2012–2014). The majority of the RVA isolates were grouped into five genotype constellations: Wa-like G1P[8], DS-1-like G1P[8], G2P[4], G3P[8] and G9P[8]. Phylogenetic analysis showed that the distribution of strains varied by geographical locations and epidemic seasons. The VP7 genes of different G types were estimated to evolve at 7.26 × 10^-4–1.04 × 10^-3 nucleotide substitutions per site per year. The Bayesian time-scaled tree of VP7 showed that the time to the most recent common ancestor of epidemic strains within a region was 1–3 years, whereas that of the epidemic strains across the country was 2–6 years. This study provided, for the first time, the timeframe during which an epidemic strain spread locally and within the country and baseline information needed to predict how rapidly RVAs spread.

The evolving epidemiology of rotavirus A infection in Brazil a decade after the introduction of universal vaccination with Rotarix®

Background

Brazil introduced the monovalent rotavirus vaccine (Rotarix®) in 2006. This study aimed to assess the epidemiology and genotype distribution of species-A rotavirus (RVA) in Brazil, comparing the pre- and post-vaccination periods.

Methods

Laboratory-based RVA surveillance included 866 municipalities in 22 Brazilian states, over a 21-year period. A total of 16,185 children with diarrheal diseases (DD) aged up to 12 years between 1996 and 2005 (pre-vaccination period, n = 7030) and from 2006 to 2017 (post-vaccination period, n = 9155) were enrolled. RVA was detected using ELISA immune assay and/or polyacrylamide gel electrophoresis and genotyped using nested PCR and/or nucleotide sequencing. RVA-positivity and genotypes detection rates were compared in distinct periods and age groups and Rotarix vaccination status.

Results

RVA-positivity in pre- and post-vaccination periods was, respectively: 4–11 months bracket, 33.3% (668/2006) and 16.3% (415/2547) (p <  0.001); 12–24 months, 28.2% (607/2154) and 22.2% (680/3068) (p <  0.001); 25–48 months, 17.4% (215/1235) and 29.4% (505/1720) (p <  0.001). Genotypes distribution in the pre- and post-vaccination periods was, respectively: G1P [8]/G1P[Not Typed], 417/855 (48.8%) and 118/1835 (6.4%) (p <  0.001); G2P [4]/G2P[NT], 47/855 (5.5%) and 838/1835 (45.7%) (p <  0.001); G3P [8]/G3P[NT], 55/855 (6.4%) and 253/1835 (13.8%) (p <  0.001); G9P [8]/G9P[NT], 238/855 (27.8%) and 152/1835 (8.3%) (p <  0.001); G12P [8]/G129P[NT], 0/871 (0%) and 249/1835(13.6%) (p <  0.001). Concerning infants aged 4–11 months, RVA frequency in fully vaccinated and non-vaccinated individuals was 11.9% (125/1052) and 24.5% (58/237) (p <  0.001), respectively. In children aged 12–24 months, RVA detection rate was 18.1% (253/1395) and 29.6% (77/260) (p <  0.001), for the vaccinated and non-vaccinated individuals, respectively (p <  0.001).

Conclusions

RVA infection was significantly less frequent in children aged ≤2 years with DD after implementing vaccination, mainly among vaccinated children. It was also observed a decrease of P [8] circulation and emergence of G2P[4] in 2005, and afterwards in the post-vaccine era, with spreading of G12P[8] in 2014–2015 and of G3P[8] in 2017. Continuous RVA surveillance must be carried out in this scenario.

Genome constellations of 24 porcine rotavirus group A strains circulating on commercial Thai swine farms between 2011 and 2016

Rotavirus A (RVA) infection is a major cause of diarrhea-related illness in young children. RVA is also one of the most common enteric viruses detected on pig farms and contributes to substantial morbidity and mortality in piglets. Long-term multi-site surveillance of RVA on Thai swine farms to determine the diversity of RVA strains in circulation is currently lacking. In this study, we characterized the 11 segments of the RVA genome from 24 Thai porcine RVA strains circulating between 2011 and 2016. We identified G9 (15/24) and P[13] (12/24) as the dominant genotypes. The dominant G and P combinations were G9P[13] (n = 6), G9P[23] (n = 6), G3P[13] (n = 5), G9P[19] (n = 3), G4P[6] (n = 2), G4P[19] (n = 1), and G5P[13] (n = 1). Genome constellation of the Thai strains showed the predominance of Wa-like genotype (Gx-P[x]-I1/I5-R1-C1-M1-A8-N1-T1/T7-E1/E9-H1) with evidence of reassortment between the porcine and human RVA strains (e.g., G4-P[6]-I1-R1-C1-M1-A8-N1-T1-E1-H1 and G9-P[19]-I5-R1-C1-M1-A8-N1-T7-E9-H1). To assess the potential effectiveness of rotavirus vaccination, the Thai RVA strains were compared to the RVA strains represented in the swine rotavirus vaccine, which showed residue variations in the antigenic epitope on VP7 and shared amino acid identity below 90% for G4 and G5 strain. Several previous studies suggested these variations might effect on virus neutralization specificity and vaccine efficacy. Our study illustrates the importance of RVA surveillance beyond the G/P genotyping on commercial swine farms, which is crucial for controlling viral transmission.

Development of a live attenuated trivalent porcine rotavirus A vaccine against disease caused by recent strains most prevalent in South Korea

Porcine rotaviruses cause severe economic losses in the Korean swine industry due to G- and P-genotype mismatches between the predominant field and vaccine strains. Here, we developed a live attenuated trivalent porcine group A rotavirus vaccine using 80 cell culture passages of the representative Korean predominant strains G8P[7] 174-1, G9P[23] PRG942, and G5P[7] K71. Vaccination with the trivalent vaccine or its individual components induced no diarrhea during the first 2 weeks post-vaccination, i.e., the vaccines were attenuated. Challenge of trivalent-vaccinated or component-vaccinated piglets with homologous virulent strain(s) did not induce diarrhea for 2 weeks post-challenge. Immunization with the trivalent vaccine or its individual components also alleviated the histopathological lesions in the small intestines caused by challenge with the corresponding original virulent strain(s). Fecal secretory IgAs specific for each of vaccine strains were detected starting at 14 days post-vaccination (dpv), and IgA levels gradually increased up to 28 dpv. Oral immunization with the trivalent vaccine or its individual components induced high levels of serum virus-neutralizing antibody by 7 dpv. No diarrhea was observed in any experimental piglets during five consecutive passages of each vaccine strain. Our data indicated that the live attenuated trivalent vaccine was safe and effective at protecting piglets from diarrhea induced by challenge exposure of homologous virulent strains. This trivalent vaccine will potentially contribute toward controlling porcine rotavirus disease in South Korea and other countries where rotavirus infections with similar G and P genotypes are problematic.

Group A Rotavirus VP1 Polymerase and VP2 Core Shell Proteins: Intergenotypic Sequence Variation and In Vitro Functional Compatibility

Group A rotaviruses (RVAs) are classified according to a nucleotide sequence-based system that assigns a genotype to each of the 11 double-stranded RNA (dsRNA) genome segments. For the segment encoding the VP1 polymerase, 22 genotypes (R1 to R22) are defined with an 83% nucleotide identity cutoff value. For the segment encoding the VP2 core shell protein, which is a functional VP1-binding partner, 20 genotypes (C1 to C20) are defined with an 84% nucleotide identity cutoff value. However, the extent to which the VP1 and VP2 proteins encoded by these genotypes differ in their sequences or interactions has not been described. Here, we sought to (i) delineate the relationships and sites of variation for VP1 and VP2 proteins belonging to the known RVA genotypes and (ii) correlate intergenotypic sequence diversity with functional VP1-VP2 interaction(s) during dsRNA synthesis. Using bioinformatic approaches, we revealed which VP1 and VP2 genotypes encode divergent proteins and identified the positional locations of amino acid changes in the context of known structural domains/subdomains. We then employed an in vitro dsRNA synthesis assay to test whether genotype R1, R2, R4, and R7 VP1 polymerases could be enzymatically activated by genotype C1, C2, C4, C5, and C7 VP2 core shell proteins. Genotype combinations that were incompatible informed the rational design and in vitro testing of chimeric mutant VP1 and VP2 proteins. The results of this study connect VP1 and VP2 nucleotide-level diversity to protein-level diversity for the first time, and they provide new insights into regions/residues critical for VP1-VP2 interaction(s) during viral genome replication.IMPORTANCE Group A rotaviruses (RVAs) are widespread in nature, infecting numerous mammalian and avian hosts and causing severe gastroenteritis in human children. RVAs are classified using a system that assigns a genotype to each viral gene according to its nucleotide sequence. To date, 22 genotypes have been described for the gene encoding the viral polymerase (VP1), and 20 genotypes have been described for the gene encoding the core shell protein (VP2). Here, we analyzed if/how the VP1 and VP2 proteins encoded by the known RVA genotypes differ from each other in their sequences. We also used a biochemical approach to test whether the intergenotypic sequence differences influenced how VP1 and VP2 functionally engage each other to mediate RNA synthesis in a test tube. This work is important because it increases our understanding of RVA protein-level diversity and raises new ideas about the VP1-VP2 binding interface(s) that is important for viral replication.