Evolutionary dynamics of human rotaviruses: balancing reassortment with preferred genome constellations

Isolation and identification of BRV G6P[1] strain in Heilongjiang province, Northeast China

Bovine rotavirus (BRV) is the main cause of acute gastroenteritis in calves, resulting in significant economic losses to the cattle industry worldwide. Additionally, BRV has multiple genotypes, which could enable cross-species transmission, thereby posing a significant risk to public health. However, there is a problem of multiple genotypes coexisting in BRV, and the cross-protection effect between different genotypes of rotavirus strains is not effective enough. Therefore, mastering clinical epidemic genotypes and using epidemic genotype strains for vaccine preparation is an effective means of preventing and controlling BRV. In this study, BRV strain DQ2020 in MA104 cells was identified by transmission electron microscopy (TEM), reverse transcription polymerase chain reaction (RT-PCR), and colloidal gold immunochromatographic test strips. The whole genome of BRV strain DQ2020 was sequenced and pathogenicity in suckling mice was assessed. The results showed that after 10 passages in MA104 cells, BRV strain DQ2020 induced cytopathic effects. Wheel-shaped virus particles (diameter, ~80 nm) were observed by TEM. A target band of 382 bp was detected by RT-PCR, a positive band was detected with the colloidal gold immunochromatographic test strips, and significant green fluorescence was observed by indirect immunofluorescence (IFA). The highest median tissue culture infectious dose of strain DQ2020 after 9 passages in MA104 cells was 10-4.81 viral particles/0.1 mL. Based on phylogenetic analysis of 11 gene fragments, the genotype of BRV strain DQ2020 was G6-P[1]-I2-R2-C2-M2-A11-N2-T6-E2-H3, confirming transmission of the G6-P[1] genotype in Chinese cattle herds. Further analysis showed that the isolated strain was a reassortant of bovine (VP7, VP6, NSP3, and NSP5), human (VP4, VP1, VP2, VP3, NSP2, and NSP4), and ovine (NSP1) rotaviruses. BRV strain DQ2020 caused damage to the intestinal villi of suckling mice and diarrhea, confirming pathogenicity. In summary, this study identified a reassortant strain of bovine, human, and ovine rotavirus that is pathogenic to lactating mice, and conducted whole genome sequence analysis, providing valuable insights for the genetic evolution of the virus and the development of vaccines.

Genetic and antigenic characterization of two diarrhoeicdominant rotavirus A genotypes G3P[12] and G14P[12] circulating in the global equine population

Equine rotavirus species A (ERVA) G3P[12] and G14P[12] are two dominant genotypes that cause foal diarrhoea with a significant economic impact on the global equine industry. ERVA can also serve as a source of novel (equine-like) rotavirus species A (RVA) reassortants with zoonotic potential as those identified previously in 2013-2019 when equine G3-like RVA was responsible for worldwide outbreaks of severe gastroenteritis and hospitalizations in children. One hurdle to ERVA research is that the standard cell culture system optimized for human rotavirus replication is not efficient for isolating ERVA. Here, using an engineered cell line defective in antiviral innate immunity, we showed that both equine G3P[12] and G14P[12] strains can be rapidly isolated from diarrhoeic foals. The genome sequence analysis revealed that both G3P[12] and G14P[12] strains share the identical genotypic constellation except for VP7 and VP6 segments in which G3P[12] possessed VP7 of genotype G3 and VP6 of genotype I6 and G14P[12] had the combination of VP7 of genotype G14 and VP6 of genotype I2. Further characterization demonstrated that two ERVA genotypes have a limited cross-neutralization. The lack of an in vitro broad cross-protection between both genotypes supported the increased recent diarrhoea outbreaks due to equine G14P[12] in foals born to dams immunized with the inactivated monovalent equine G3P[12] vaccine. Finally, using the structural modelling approach, we provided the genetic basis of the antigenic divergence between ERVA G3P[12] and G14P[12] strains. The results of this study will provide a framework for further investigation of infection biology, pathogenesis and cross-protection of equine rotaviruses.

Emergence of a novel pathogenic porcine G1P[7] rotavirus in China

Among group A rotaviruses (RVAs), the G1 genotype is the main genotype causing diarrhea in children, but it has rarely been reported in pigs. During our epidemiological investigation, we detected G1P[7] rotavirus infection in piglets across several provinces in China and then isolated a porcine G1P[7] rotavirus strain (CN1P7). Sequencing revealed that the virus constellation was G1-P[7]-I5-R1-C1-M1-A8-N1-T1-E1-H1. Phylogenetic analyses revealed that CN1P7 most likely emerged due to genetic reassortment among porcine, human, giant panda and dog rotavirus strains. In vivo experiments were conducted on two-day-old piglets, which revealed that the CN1P7 strain was pathogenic to piglets. The virus was shed through the digestive tract and respiratory tract. In addition to the intestine, the CN1P7 strain displayed extraintestinal tropisms in piglets. Histopathological analysis revealed that the lung and small intestine were the targets of CN1P7. This study is the first to explore the molecular and pathogenic characterization of a pig-origin G1P[7] rotavirus.

Whole-genome analysis of human group A rotaviruses in 1980s Japan and evolutionary assessment of global Wa-like strains across half a century

Historically, the Wa-like strains of human group A rotavirus (RVA) have been major causes of gastroenteritis. However, since the 2010s, the circulation of non-Wa-like strains has been increasingly reported, indicating a shift in the molecular epidemiology of RVA. Although understanding RVA evolution requires the analysis of both current and historical strains, comprehensive pre-1980's sequencing data are scarce globally. We determined the whole-genome sequences of representative strains from six RVA gastroenteritis outbreaks observed at an infant home in Sapporo, Japan, between 1981 and 1989. These outbreaks were mainly caused by G1 or G3 Wa-like strains, resembling strains from the United States in the 1970s-1980s and from Malawi in the 1990s. Phylogenetic analysis of these infant home strains, together with Wa-like strains collected worldwide from the 1970s to 2020, revealed a notable trend: pre-2010 strains diverged into multiple lineages in many genomic segments, whereas post-2010 strains tended to converge into a single lineage. However, Bayesian skyline plot indicated near-constant effective population sizes from the 1970s to 2020, and selection pressure analysis identified positive selection only at amino acid 75 of NSP2. These results suggest that evidence supporting the influence of rotavirus vaccines, introduced globally since 2006, on Wa-like RVA molecular evolution is lacking at present, and phylogenetic analysis may simply reflect natural fluctuations in RVA molecular evolution. Evaluating the long-term impact of RV vaccines on the molecular evolution of RVA requires sustained surveillance.

The Number and Pattern of Viral Genomic Reassortments are not Necessarily Identifiable from Segment Trees

Reassortment is an evolutionary process common in viruses with segmented genomes. These viruses can swap whole genomic segments during cellular co-infection, giving rise to novel progeny formed from the mixture of parental segments. Since large-scale genome rearrangements have the potential to generate new phenotypes, reassortment is important to both evolutionary biology and public health research. However, statistical inference of the pattern of reassortment events from phylogenetic data is exceptionally difficult, potentially involving inference of general graphs in which individual segment trees are embedded. In this paper, we argue that, in general, the number and pattern of reassortment events are not identifiable from segment trees alone, even with theoretically ideal data. We call this fact the fundamental problem of reassortment, which we illustrate using the concept of the "first-infection tree," a potentially counterfactual genealogy that would have been observed in the segment trees had no reassortment occurred. Further, we illustrate four additional problems that can arise logically in the inference of reassortment events and show, using simulated data, that these problems are not rare and can potentially distort our observation of reassortment even in small data sets. Finally, we discuss how existing methods can be augmented or adapted to account for not only the fundamental problem of reassortment, but also the four additional situations that can complicate the inference of reassortment.

VP4 Mutation Boosts Replication of Recombinant Human/Simian Rotavirus in Cell Culture

Rotavirus A (RVA) is the leading cause of diarrhea requiring hospitalization in children and causes over 100,000 annual deaths in Sub-Saharan Africa. In order to generate next-generation vaccines against African RVA genotypes, a reverse genetics system based on a simian rotavirus strain was utilized here to exchange the antigenic capsid proteins VP4, VP7 and VP6 with those of African human rotavirus field strains. One VP4/VP7/VP6 (genotypes G9-P[6]-I2) triple-reassortant was successfully rescued, but it replicated poorly in the first cell culture passages. However, the viral titer was enhanced upon further passaging. Whole genome sequencing of the passaged virus revealed a single point mutation (A797G), resulting in an amino acid exchange (E263G) in VP4. After introducing this mutation into the VP4-encoding plasmid, a VP4 mono-reassortant as well as the VP4/VP7/VP6 triple-reassortant replicated to high titers already in the first cell culture passage. However, the introduction of the same mutation into the VP4 of other human RVA strains did not improve the rescue of those reassortants, indicating strain specificity. The results show that specific point mutations in VP4 can substantially improve the rescue and replication of recombinant RVA reassortants in cell culture, which may be useful for the development of novel vaccine strains.

Phylodynamic characteristics of reassortant DS-1-like G3P[8]-strains of rotavirus type A isolated in Nizhny Novgorod (Russia)

Since 2013, there has been an increase in reports of the spread of a double intergroup reassortant strain of rotavirus type A (RVA) with the genotype G3P[8] and other genes belonging to the second genogroup I2-R2-C2-M2-A2-N2-T2-E2-H2. In our study, we provide a molecular genetic characterization of rotaviruses with genotype G3P[8]-I2 isolated in Nizhny Novgorod. In our study, we used RT-PCR, Sanger sequencing, RNA-PAGE methods. Phylogenetic and phylodynamic analysis were performed using the Bayesian approach. According to our study, there was a significant increase in the proportion of G3P[8] from 15% during the period of 2020-2021 to 53% during the period of 2021-2022 in Nizhny Novgorod, Russia. Phylogenetic analysis based on the VP4 gene revealed that DS-1-like RVAs isolated in Nizhny Novgorod belong to different clusters of the P[8]-3.1 lineage, with a level of variation ranging from 1.1% to 1.3%. Based on the VP6 gene, the equine-like RVAs identified by us carry genetic variants belonging to three distinct clusters of the lineage I2-V, with a variation level ranging from 2.0% to 4.5%. These data indicate the genotypic diversity of circulating DS-1-like G3 RVAs. Phylogenetic analysis of the VP7 gene allowed us to assign the isolates identified in our study to the G3-1 lineage. We estimated that the circulation of the most recent common ancestor of the spreading strains dates back to 2002. Additionally, we determined the typical level of mutations in the VP7 gene, which amounted to 2.14*10-3 substitutions/per site/per year.

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.

Modeling the evolution of segment trees reveals deficiencies in current inferential methods for genomic reassortment

Reassortment is an evolutionary process common in viruses with segmented genomes. These viruses can swap whole genomic segments during cellular co-infection, giving rise to new viral variants. Large-scale genome rearrangements, such as reassortment, have the potential to quickly generate new phenotypes, making the understanding of viral reassortment important to both evolutionary biology and public health research. In this paper, we argue that reassortment cannot be reliably inferred from incongruities between segment phylogenies using the established remove-and-rejoin or coalescent approaches. We instead show that reassortment must be considered in the context of a broader population process that includes the dynamics of the infected hosts. Using illustrative examples and simulation we identify four types of evolutionary events that are difficult or impossible to reconstruct with incongruence-based methods. Further, we show that these specific situations are very common and will likely occur even in small samples. Finally, we argue that existing methods can be augmented or modified to account for all the problematic situations that we identify in this paper. Robust assessment of the role of reassortment in viral evolution is difficult, and we hope to provide conceptual clarity on some important methodological issues that can arise in the development of the next generation of tools for studying reassortment.

Surveillance of Human Rotaviruses in Wuhan, China (2019-2022): Whole-Genome Analysis of Emerging DS-1-like G8P[8] Rotavirus

Group A rotaviruses (RVAs) are major etiologic agents of gastroenteritis in infants and young children worldwide. To study the prevalence and genetic characteristics of RVAs, a hospital-based surveillance study was conducted in Wuhan, China from June 2019 through May 2022. The detection rates of RVAs were 19.40% (142/732) and 3.51% (8/228) in children and adults, respectively. G9P[8] was the predominant genotype, followed by G8P[8] and G3P[8]. G8P[8] emerged and was dominant in the 2021-2022 epidemic season. The genome constellation of six G8P[8] strains was assigned to G8-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Phylogenetic analysis revealed that the VP7, VP4, VP2, VP3, NSP1, NSP2, NSP3, and NSP5 genes of these G8P[8] strains clustered closely with those of the G8P[8] strains in Asia and were distant from those of the P[8] and G2P[4] strains simultaneously detected in Wuhan. In contrast, the VP1, VP6, and NSP4 genes were closely related to the typical G2P[4] rotavirus, including those of G2P[4] strains simultaneously detected in Wuhan. The detection rate of RVAs decreased in the COVID-19 pandemic era. It was deduced that the G8P[8] rotaviruses that emerged in China may be reassortants, carrying the VP6, VP1, and NSP4 genes derived from the G2P[4] rotavirus in the backbone of the neighboring DS-1-like G8P[8] strains represented by CAU17L-103.

The Evolution of Post-Vaccine G8P[4] Group a Rotavirus Strains in Rwanda; Notable Variance at the Neutralization Epitope Sites

Africa has a high level of genetic diversity of rotavirus strains, which is suggested to be a possible reason contributing to the suboptimal effectiveness of rotavirus vaccines in this region. One strain that contributes to this rotavirus diversity in Africa is the G8P[4]. This study aimed to elucidate the entire genome and evolution of Rwandan G8P[4] strains. Illumina sequencing was performed for twenty-one Rwandan G8P[4] rotavirus strains. Twenty of the Rwandan G8P[4] strains had a pure DS-1-like genotype constellation, and one strain had a reassortant genotype constellation. Notable radical amino acid differences were observed at the neutralization sites when compared with cognate regions in vaccine strains potentially playing a role in neutralization escape. Phylogenetic analysis revealed that the closest relationship was with East African human group A rotavirus (RVA) strains for five of the genome segments. Two genome sequences of the NSP4 genome segment were closely related to bovine members of the DS-1-like family. Fourteen VP1 and eleven VP3 sequences had the closest relationships with the RotaTeq™ vaccine WC3 bovine genes. These findings suggest that the evolution of VP1 and VP3 might have resulted from reassortment events with RotaTeq™ vaccine WC3 bovine genes. The close phylogenetic relationship with East African G8P[4] strains from Kenya and Uganda suggests co-circulation in these countries. These findings highlight the need for continued whole-genomic surveillance to elucidate the evolution of G8P[4] strains, especially after the introduction of rotavirus vaccination.

Strain-Specific Interactions between the Viral Capsid Proteins VP4, VP7 and VP6 Influence Rescue of Rotavirus Reassortants by Reverse Genetics

Rotavirus A (RVA) genome segments can reassort upon co-infection of target cells with two different RVA strains. However, not all reassortants are viable, which limits the ability to generate customized viruses for basic and applied research. To gain insight into the factors that restrict reassortment, we utilized reverse genetics and tested the generation of simian RVA strain SA11 reassortants carrying the human RVA strain Wa capsid proteins VP4, VP7, and VP6 in all possible combinations. VP7-Wa, VP6-Wa, and VP7/VP6-Wa reassortants were effectively rescued, but the VP4-Wa, VP4/VP7-Wa, and VP4/VP6-Wa reassortants were not viable, suggesting a limiting effect of VP4-Wa. However, a VP4/VP7/VP6-Wa triple-reassortant was successfully generated, indicating that the presence of homologous VP7 and VP6 enabled the incorporation of VP4-Wa into the SA11 backbone. The replication kinetics of the triple-reassortant and its parent strain Wa were comparable, while the replication of all other rescued reassortants was similar to SA11. Analysis of the predicted structural protein interfaces identified amino acid residues, which might influence protein interactions. Restoring the natural VP4/VP7/VP6 interactions may therefore improve the rescue of RVA reassortants by reverse genetics, which could be useful for the development of next generation RVA vaccines.

Mammalian orthoreoviruses exhibit rare genotype variability in genome constellations

Mammalian orthoreoviruses (reoviruses) are currently classified based on properties of the attachment protein, σ1. Four reovirus serotypes have been identified, three of which are represented by well-studied prototype human reovirus strains. Reoviruses contain ten segments of double-stranded RNA that encode 12 proteins and can reassort during coinfection. To understand the breadth of reovirus genetic diversity and its potential influence on reassortment, the sequence of the entire genome should be considered. While much is known about the prototype strains, a thorough analysis of all ten reovirus genome segment sequences has not previously been conducted. We analyzed phylogenetic relationships and nucleotide sequence conservation for each of the ten segments of more than 60 complete or nearly complete reovirus genome sequences, including those of the prototype strains. Using these relationships, we defined genotypes for each segment, with minimum nucleotide identities of 77-88% for most genotypes that contain several representative sequences. We applied segment genotypes to determine reovirus genome constellations, and we propose implementation of an updated reovirus genome classification system that incorporates genotype information for each segment. For most sequenced reoviruses, segments other than S1, which encodes σ1, cluster into a small number of genotypes and a limited array of genome constellations that do not differ greatly over time or based on animal host. However, a small number of reoviruses, including prototype strain Jones, have constellations in which segment genotypes differ from those of most other sequenced reoviruses. For these reoviruses, there is little evidence of reassortment with the major genotype. Future basic research studies that focus on the most genetically divergent reoviruses may provide new insights into reovirus biology. Analysis of available partial sequences and additional complete reovirus genome sequencing may also reveal reassortment biases, host preferences, or infection outcomes that are based on reovirus genotype.

Complete genome analyses of G12P[8] rotavirus strains from hospitalized children in Surabaya, Indonesia, 2017-2018

Rotavirus A (RVA) is a major viral cause of acute gastroenteritis (AGE) worldwide. G12 RVA strains have emerged globally since 2007. There has been no report of the whole genome sequences of G12 RVAs in Indonesia. We performed the complete genome analysis by the next-generation sequencing of five G12 strains from hospitalized children with AGE in Surabaya from 2017 to 2018. All five G12 strains were Wa-like strains (G12-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1) and were clustered into lineage-III of VP7 gene phylogenetic tree. STM430 sample was observed as a mixed-infection between G12 and G1 strains: G12/G1-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1. A phylogenetic tree analysis revealed that all five Indonesian G12 strains (SOEP379, STM371, STM413, STM430, and STM433) were genetically close to each other in all 11 genome segments with 98.0%-100% nucleotide identities, except VP3 and NSP4 of STM430, suggesting that these strains have originated from a similar ancestral G12 RVA. The VP3 and NSP4 genome segments of STM430-G12P[8] were separated phylogenetically from those of the other four G12 strains, probably due to intra-genotype reassortment between the G12 and G1 Wa-like strains. The change from G12P[6] lineage-II in 2007 to G12P[8] lineage-III 2017-2018 suggests the evolution and diversity of G12 RVAs in Indonesia over the past approximately 10 years.

Full genotype constellations analysis of unusual DS-1-like G12P[6] and G6P[8] rotavirus strains detected in Brazil, 2019

Rotavirus A (RVA) is a major cause of acute gastroenteritis (AGE) in children worldwide. We report unusual RVA G12P[6] and G6P[8] strains isolated from fecal samples from Brazilian children hospitalized for AGE. The characterized RVA have genome segments backbone: G12-P[6]/ G6-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2 of DS-1-like genogroup. Our study describes the first identification of G6P[8], a DS-1-like genogroup strain. Nucleotide analysis of VP7 and VP4 genes revealed that all G12 Brazilian strains clustered into the sub-lineages IIIB, mostly associated with P[6] lineage I. Additionally, our G6 lineage I strains were closely related to German G6 genotypes, bound with P[8] lineage III, differing from both vaccine strains. The comparative sequence analysis of our strains with vaccine strains revealed amino acid substitutions located in immunodominant regions of VP7 and VP4 proteins. Continuous monitoring of RVA genotypes is essential to evaluate the impact of vaccination on the dynamic nature of RVA evolution.

Genome Analysis of the G6P6 Genotype of Porcine Group C Rotavirus in China

Swine enteric disease is the predominant cause of morbidity and mortality, and viral species involved in swine enteric disease include rotaviruses and coronaviruses, among others. Awareness of the circulating porcine rotavirus group C (PoRVC) in pig herds is critical to evaluate the potential impact of infection. At present, due to the lack of disease awareness and molecular diagnostic means, the research on RVC infection in China is not well-studied. In this study, diarrhea samples collected from pig farms were detected positive for RVC by PCR, and the full-length RVC was not previously reported for Chinese pig farms. This rotavirus strain was designated as RVC/Pig/CHN/JS02/2018/G6P6. A natural recombination event was observed with breakpoints at nucleotides (nt) 2509 to 2748 of the VP2 gene. Phylogenetic analysis based on nsp1 revealed that a new branch A10 formed. Collectively, our data suggest a potentially novel gene recombination event of RVC in the VP2 gene. These findings provide a new insight into the evolution of the rotavirus.

Genomic and evolutionary characteristics of G9P[8], the dominant group a rotavirus in China (2016–2018)

G9P[8] became the predominant rotavirus A (RVA) genotype in China in 2012. To evaluate its genetic composition at the whole-genome level, 115 G9P[8] RVA strains isolated from children under 5 years old were sequenced and characterized. All 13 strains in 2016 and 2017 and an additional 54 strains in 2018 were genotyped as G9-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1. The other 48 strains in 2018 were all genotyped as G9-P[8]-I1-R1-C1-M1-A1-N1-T1-E2-H1, with the NSP4 gene characterized as a DS-1-like genotype. The time of the most recent common ancestor (tMRCA) and evolution rates of the VP7, VP4, and NSP4 (E1 and E2) genes of these strains were estimated by Bayesian evolutionary dynamics analysis. We estimated the evolution rates (nt substitutions per site per year) as 1.38 × 10^–3 [the 95% highest posterior density (HPD) was 1.09–1.72 × 10^–3] for VP7, 0.87 × 10^–3 (95% HPD: 0.75–1.00 × 10^–3) for VP4, 0.56 × 10^–3 (95% HPD: 0.41–0.73 × 10^–3) for NSP4-E1, and 1.35 × 10^–3 (95% HPD: 0.92–1.86 × 10^–3) for NSP4-E2. The tMRCA was estimated to be 1935.4 (95% HPD: 1892.4–1961.3) for VP7, 1894.3 (95% HPD: 1850.5–1937.8) for VP4, 1929.4 (95% HPD: 1892.4–1961.3) for NSP4-E1, and 1969.2 (95% HPD: 1942.2–1985.3) for NSP4-E2. The baseline genetic information in this study is expected to improve our understanding of the genomic and evolutionary characteristics of the rotavirus genome. Furthermore, it will provide a basis for the development of next-generation rotavirus vaccines for humans.

Rotavirus vaccine-derived cases in Belgium: Evidence for reversion of attenuating mutations and alternative causes of gastroenteritis

Since the introduction of live-attenuated rotavirus vaccines in Belgium in 2006, surveillance has routinely detected rotavirus vaccine-derived strains. However, their genomic landscape and potential role in gastroenteritis have not been thoroughly investigated. We compared VP7 and VP4 nucleotide sequences obtained from rotavirus surveillance with the Rotarix vaccine sequence. As a result, we identified 80 vaccine-derived strains in 5125 rotavirus-positive infants with gastroenteritis from 2007 to 2018. Using both viral metagenomics and reverse transcription qPCR, we evaluated the vaccine strains and screened for co-infecting enteropathogens. Among the 45 patients with known vaccination status, 39 were vaccinated and 87% received the vaccine less than a month before the gastroenteritis episode. Reconstruction of 30 near complete vaccine-derived genomes revealed 0-11 mutations per genome, with 88% of them being non-synonymous. This, in combination with several shared amino acid changes among strains, pointed at selection of minor variant(s) present in the vaccine. We also found that some of these substitutions were true revertants (e.g., F167L on VP4, and I45T on NSP4). Finally, co-infections with known (e.g., Clostridioides difficile and norovirus) and divergent or emerging (e.g., human parechovirus A1, salivirus A2) pathogens were detected, and we estimated that 35% of the infants likely had gastroenteritis due to a 'non-rotavirus' cause. Conversely, we could not rule out the vaccine-derived gastroenteritis in over half of the cases. Continued studies inspecting reversion to pathogenicity should monitor the long-time safety of live-attenuated rotavirus vaccines. All in all, the complementary approach with NGS and qPCR provided a better understanding of rotavirus vaccine strain evolution in the Belgian population and epidemiology of co-infecting enteropathogens in suspected rotavirus vaccine-derived gastroenteritis cases.

Occurrence of Rotavirus A Genotypes and Other Enteric Pathogens in Diarrheic Suckling Piglets from Spanish Swine Farms

Simple Summary

Neonatal diarrhea is a major cause of economic losses in the swine industry worldwide and has significant impact in Spain, which is one of the biggest pork producers globally. Multiple infectious agents can contribute to this condition, with some viruses such as species A rotavirus (RVA) playing a major role. Studies on their occurrence and genetic diversity are essential for development of RVA vaccines. In this study, fecal samples from diarrheic suckling piglets originating from farms distributed throughout Spain were analyzed for RVA and four other common enteric pathogens using molecular methods. The individual prevalence was 89.4%, 64.4%, 44.9%, 33.7% and 4.4% for Clostridiumperfringens, Clostridioides (formerly Clostridium) difficile, species A rotavirus, species C rotavirus and porcine epidemic diarrhea virus, respectively. Most specimens (96.9%) were positive for at least one of the target pathogens and concurrent infections were common. The molecular characterization of RVA positive specimens of specific genes used for genotyping revealed the extensive genetic diversity of RVA strains circulating in swine herds in Spain. Comparison with genotypes contained in the commercial vaccine available in Spain showed differences in the identity of the predominant RVA genotypes from diarrheic piglets in the sampled pig farms. These findings contribute to the surveillance of RVA strains circulating in swine herds in Spain and may help optimize target vaccine design.

Abstract

Species A rotavirus (RVA) is a major viral pathogen causing diarrhea in suckling piglets. Studies on its genetic heterogeneity have implications for vaccine efficacy in the field. In this study, fecal samples (n = 866) from diarrheic piglets younger than 28 days were analyzed over a two-year period (2018–2019). Samples were submitted from 426 farms located in 36 provinces throughout Spain and were tested using real-time PCR (qPCR) and reverse transcription real-time PCR (RT-qPCR) for five enteric pathogens. The individual prevalence was 89.4%, 64.4%, 44.9%, 33.7% and 4.4% for Clostridiumperfringens, Clostridioides (formerly Clostridium) difficile, species A rotavirus, species C rotavirus and porcine epidemic diarrhea virus, respectively. Most specimens (96.9%) were positive for at least one of the target pathogens, and more than 80% of samples harbored mixed infections. Nucleotide sequencing of 70 specimens positive for RVA revealed the presence of the VP7 genotypes G4, G9, G3, G5, G11 and the VP4 genotypes P7, P23, P6 and P13, with the combinations G4P7 and G9P23 being the most prevalent, and especially in the areas with the highest pig population. The study shows the extensive genetic diversity of RVA strains as well as discrepancies with the genotypes contained in the vaccine available in Spain, and multiple amino acid differences in antigenic epitopes of different G- and P- genotypes with the vaccine strains. Further investigations are needed to determine the efficacy of the vaccine to confer clinical protection against heterologous strains.

High Prevalence of Coinfecting Enteropathogens in Suspected Rotavirus Vaccine Breakthrough Cases

Despite the global use of rotavirus vaccines, vaccine breakthrough cases remain a pediatric health problem. In this study, we investigated suspected rotavirus vaccine breakthrough cases using next-generation sequencing (NGS)-based viral metagenomics (n = 102) and a panel of semiquantitative reverse transcription-PCR (RT-qPCR) (n = 92) targeting known enteric pathogens. Overall, we identified coinfections in 80% of the cases. Enteropathogens such as adenovirus (32%), enterovirus (15%), diarrheagenic Escherichia coli (1 to 14%), astrovirus (10%), Blastocystis spp. (10%), parechovirus (9%), norovirus (9%), Clostridioides (formerly Clostridium) difficile (9%), Dientamoeba fragilis (9%), sapovirus (8%), Campylobacter jejuni (4%), and Giardia lamblia (4%) were detected. Except for a few reassortant rotavirus strains, unusual genotypes or genotype combinations were not present. However, in addition to well-known enteric viruses, divergent variants of enteroviruses and nonclassic astroviruses were identified using NGS. We estimated that in 31.5% of the patients, rotavirus was likely not the cause of gastroenteritis, and in 14.1% of the patients, it contributed together with another pathogen(s) to disease. The remaining 54.4% of the patients likely had a true vaccine breakthrough infection. The high prevalence of alternative enteropathogens in the suspected rotavirus vaccine breakthrough cases suggests that gastroenteritis is often the result of a coinfection and that rotavirus vaccine effectiveness might be underestimated in clinical and epidemiological studies.

Development of a Real-Time Reverse Transcription-PCR Assay To Detect and Quantify Group A Rotavirus Equine-Like G3 Strains

Since 2013, group A rotavirus strains characterized as novel DS-1-like intergenogroup reassortant "equine-like G3" strains have emerged and spread across 5 continents among human populations in at least 14 countries. Here, we report a novel one-step TaqMan quantitative real-time reverse transcription-PCR assay developed to genotype and quantify the viral load for samples containing rotavirus equine-like G3 strains. Using a universal G forward primer and a newly designed reverse primer and TaqMan probe, we developed and validated an assay with a linear dynamic range of 227 to 2.3 × 10⁹ copies per reaction and a limit of detection of 227 copies. The percent positive agreement, percent negative agreement, and precision of our assay were 100.00%, 99.63%, and 100.00%, respectively. This assay can simultaneously detect and quantify the viral load for samples containing DS-1-like intergenogroup reassortant equine-like G3 strains with high sensitivity and specificity, faster turnaround time, and decreased cost. It will be valuable for high-throughput screening of stool samples collected to monitor equine-like G3 strain prevalence and circulation among human populations throughout the world.

FAST Proteins: Development and Use of Reverse Genetics Systems for Reoviridae Viruses

Reverse genetics systems for viruses, the technology used to generate gene-engineered recombinant viruses from artificial genes, enable the study of the roles of the individual nucleotides and amino acids of viral genes and proteins in infectivity, replication, and pathogenicity. The successful development of a reverse genetics system for poliovirus in 1981 accelerated the establishment of protocols for other RNA viruses important for human health. Despite multiple efforts, rotavirus (RV), which causes severe gastroenteritis in infants, was refractory to reverse genetics analysis, and the first complete reverse genetics system for RV was established in 2017. This novel technique involves use of the fusogenic protein FAST (fusion-associated small transmembrane) derived from the bat-borne Nelson Bay orthoreovirus, which induces massive syncytium formation. Co-transfection of a FAST-expressing plasmid with complementary DNAs encoding RV genes enables rescue of recombinant RV. This review focuses on methodological insights into the reverse genetics system for RV and discusses applications and potential improvements to this system.

Whole Genome Analysis of Human Rotaviruses Reveals Single Gene Reassortant Rotavirus Strains in Zambia

Rotarix® vaccine was implemented nationwide in Zambia in 2013. In this study, four unusual strains collected in the post-vaccine period were subjected to whole genome sequencing and analysis. The four strains possessed atypical genotype constellations, with at least one reassortant genome segment within the constellation. One of the strains (UFS-NGS-MRC-DPRU4749) was genetically and phylogenetically distinct in the VP4 and VP1 gene segments. Pairwise analyses demonstrated several amino acid disparities in the VP4 antigenic sites of this strain compared to that of Rotarix®. Although the impact of these amino acid disparities remains to be determined, this study adds to our understanding of the whole genomes of reassortant strains circulating in Zambia following Rotarix® vaccine introduction.

Human Rotavirus Reverse Genetics Systems to Study Viral Replication and Pathogenesis

Human rotaviruses (HuRVAs) are highly important causes of acute gastroenteritis in infants and young children worldwide. A lack of reliable and reproducible reverse genetics systems for HuRVAs has limited a proper understanding of HuRVA biology and also the rational design of live-attenuated vaccines. Since the development of the first reverse genetics system for RVAs (partially plasmid-based reverse genetics system) in 2006, there have been many efforts with the goal of generating infectious recombinant HuRVAs entirely from cloned cDNAs. However, the establishment of a HuRVA reverse genetics system was very challenging until 2019. This review article provides an overview of the historical background of the recent development of long-awaited HuRVA reverse genetics systems, beginning with the generation of recombinant human-simian reassortant RVAs with the aid of a helper virus in 2006 and the generation of recombinant animal (simian) RVAs in a helper virus-free manner in 2017, and culminating in the generation of recombinant HuRVAs entirely from plasmid cDNAs in 2019. Notably, the original HuRVA reverse genetics system has already been optimized to increase the efficiency of virus generation. Although the application of HuRVA reverse genetics systems has only just been initiated, these technologies will help to answer HuRVA research questions regarding viral replication and pathogenicity that could not be addressed before, and to develop next-generation vaccines and intestine-specific rotaviral vectors.

Detection and molecular characteristics of bovine rotavirus A in dairy calves in China

Background

Bovine group A rotavirus (BoRVA) is a major cause of severe gastroenteritis in newborn dairy calves. Only one study has investigated the G and P genotypes among dairy calves in a few regions of China, which were G6 and P[5]. Therefore, data on the prevalence and molecular characteristics of BoRVA in dairy calves in China remains limited.

Objectives

The purpose of this study was to investigate the prevalence and molecular characteristics of BoRVA in dairy calves in China.

Methods

269 dairy calves diarrheic samples from 23 farms in six provinces in China were collected to detect BoRVA using reverse transcription polymerase chain reaction.

Results

71% of samples were determined to be BoRVA-positive. Two G genotypes (G6, G10) and two P genotypes (P[1], P[5]) were identified, and G6P[1] BoRVA was the predominant strain. Moreover, the VP7 and VP4 gene sequences of these dairy calf BoRVA strains revealed abundant genetic diversity. Interestingly, eight out of 17 complete G6 VP7 sequences were clustered into G6 lineage VI and analysis showed the strains were closely related to Chinese yak BoRVA strains.

Conclusions

The results of this study show that BoRVA circulates widely among dairy calves in China, and the dominant genotype in circulation is G6P[1], first report on molecular characteristics of complete P[5] VP4 genes in chinese dairy calves. These results will help us to further understand the prevalence and genetic evolution of BoRVA among dairy calves in China and, thus, prevent the disease more effectively.

Rotavirus A Genome Segments Show Distinct Segregation and Codon Usage Patterns

Reassortment of the Rotavirus A (RVA) 11-segment dsRNA genome may generate new genome constellations that allow RVA to expand its host range or evade immune responses. Reassortment may also produce phylogenetic incongruities and weakly linked evolutionary histories across the 11 segments, obscuring reassortment-specific epistasis and changes in substitution rates. To determine the co-segregation patterns of RVA segments, we generated time-scaled phylogenetic trees for each of the 11 segments of 789 complete RVA genomes isolated from mammalian hosts and compared the segments’ geodesic distances. We found that segments 4 (VP4) and 9 (VP7) occupied significantly different tree spaces from each other and from the rest of the genome. By contrast, segments 10 and 11 (NSP4 and NSP5/6) occupied nearly indistinguishable tree spaces, suggesting strong co-segregation. Host-species barriers appeared to vary by segment, with segment 9 (VP7) presenting the weakest association with host species. Bayesian Skyride plots were generated for each segment to compare relative genetic diversity among segments over time. All segments showed a dramatic decrease in diversity around 2007 coinciding with the introduction of RVA vaccines. To assess selection pressures, codon adaptation indices and relative codon deoptimization indices were calculated with respect to different host genomes. Codon usage varied by segment with segment 11 (NSP5) exhibiting significantly higher adaptation to host genomes. Furthermore, RVA codon usage patterns appeared optimized for expression in humans and birds relative to the other hosts examined, suggesting that translational efficiency is not a barrier in RVA zoonosis.

Genome-wide analyses of human noroviruses provide insights on evolutionary dynamics and evidence of coexisting viral populations evolving under recombination constraints

Norovirus is a major cause of acute gastroenteritis worldwide. Over 30 different genotypes, mostly from genogroup I (GI) and II (GII), have been shown to infect humans. Despite three decades of genome sequencing, our understanding of the role of genomic diversification across continents and time is incomplete. To close the spatiotemporal gap of genomic information of human noroviruses, we conducted a large-scale genome-wide analyses that included the nearly full-length sequencing of 281 archival viruses circulating since the 1970s in over 10 countries from four continents, with a major emphasis on norovirus genotypes that are currently underrepresented in public genome databases. We provided new genome information for 24 distinct genotypes, including the oldest genome information from 12 norovirus genotypes. Analyses of this new genomic information, together with those publicly available, showed that (i) noroviruses evolve at similar rates across genomic regions and genotypes; (ii) emerging viruses evolved from transiently-circulating intermediate viruses; (iii) diversifying selection on the VP1 protein was recorded in genotypes with multiple variants; (iv) non-structural proteins showed a similar branching on their phylogenetic trees; and (v) contrary to the current understanding, there are restrictions on the ability to recombine different genomic regions, which results in co-circulating populations of viruses evolving independently in human communities. This study provides a comprehensive genetic analysis of diverse norovirus genotypes and the role of non-structural proteins on viral diversification, shedding new light on the mechanisms of norovirus evolution and transmission.

Norovirus is a highly diverse enteric pathogen. The large genomic database accumulated in the last three decades advanced our understanding of norovirus diversity; however, this information is limited by geographical bias, sporadic times of collection, and missing or incomplete genome sequences. In this multinational collaborative study, we mined archival samples collected since the 1970s and sequenced nearly full-length new genomes from 281 historical noroviruses, including the first full-length genomic sequences for three genotypes. Using this novel dataset, we found evidence for restrictions in the recombination of genetically disparate viruses and that diversifying selection results in new variants with different epidemiological profiles. These new insights on the diversification of noroviruses could provide baseline information for the study of future epidemics and ultimately the prevention of norovirus infections.

Wa-1 Equine-Like G3P[8] Rotavirus from a Child with Diarrhea in Colombia

Rotavirus A (RVA) has been considered the main cause of diarrheal disease in children under five years in emergency services in both developed and developing countries. RVA belongs to the Reoviridae family, which comprises 11 segments of double-stranded RNA (dsRNA) as a genomic constellation that encodes for six structural and five to six nonstructural proteins. RVA has been classified in a binary system with Gx[Px] based on the spike protein (VP4) and the major outer capsid glycoprotein (VP7), respectively. The emerging equine-like G3P[8] DS-1-like strains reported worldwide in humans have arisen an important concern. Here, we carry out the complete genome characterization of a previously reported G3P[8] strain in order to recognize the genetic diversity of RVA circulating among infants in Colombia. A near-full genome phylogenetic analysis was done, confirming the presence of the novel equine-like G3P[8] with a Wa-like backbone for the first time in Colombia. This study demonstrated the importance of surveillance of emerging viruses in the Colombian population; furthermore, additional studies must focus on the understanding of the spread and transmission dynamic of this important RVA strain in different areas of the country.

Emergence in 2017-2019 of novel reassortant equine-like G3 rotavirus strains in Palermo, Sicily

Rotavirus A (RVA) is a major etiologic agent of gastroenteritis in children worldwide. Hospital-based surveillance of viral gastroenteritis in paediatric population in Palermo (Italy) from 2017 onwards revealed a sharp increase in G3P[8] RVAs, accounting for 71% of all the RVAs detected in 2019. This pattern had not been observed before in Italy, with G3 RVA usually being detected at rates lower than 3%. In order to investigate this unique epidemiological pattern, the genetic diversity of G3 RVAs identified during a 16-year long surveillance (2004-2019) was explored by systematic sequencing of the VP7 and VP4 genes and by whole genome sequencing of selected G3 strains, representative of the various RVA seasons. Sequence and phylogenetic analyses of the VP7 and VP4 genes revealed the emergence, in 2017 of reassortant equine-like G3P[8], which gradually replaced former G3P[8] strains. The G3P[8] circulating before 2017 showed a Wa-like constellation of genome segments while the G3P[8] that emerged in 2017 had a DS-1-like backbone. On direct inspection of the VP7 and VP4 antigenic epitopes, the equine-like G3P[8] strains possessed several amino acid variations in neutralizing regions compared with vaccine strains. The equine-like G3P[8] RVAs are a further example of the zoonotic impact of animal viruses on human health.

Rotavirus in Calves and Its Zoonotic Importance

Rotavirus is a major pathogen responsible for diarrheal disease in calves, resulting in loss of productivity and economy of farmers. However, various facets of diarrheal disease caused by rotavirus in calves in the world are inadequately understood, considering that diarrheal disease caused by rotavirus is a vital health problem in calves that interrupts production benefits with reduced weight gain and increased mortality, and its potential for zoonotic spread. The pathological changes made by rotavirus are almost exclusively limited to the small intestine that leads to diarrhea. It is environmentally distributed worldwide and was extensively studied. Reassortment is one of the important mechanisms for generating genetic diversity of rotaviruses and eventually for viral evolution. So, the primary strategy is to reduce the burden of rotavirus infections by practicing early colostrum's feeding in newborn calves, using vaccine, and improving livestock management. Rotaviruses have a wide host range, infecting many animal species as well as humans. As it was found that certain animal rotavirus strains had antigenic similarities to some human strains, this may be an indication for an animal to play a role as a source of rotavirus infection in humans. Groups A to C have been shown to infect both humans and animals. The most commonly detected strains in both human and animals are G2, G3, G4, and G9, P [6]. Therefore, this review was made to get overview epidemiology status and zoonotic importance of bovine rotavirus.

Rescue of Infectious Rotavirus Reassortants by a Reverse Genetics System Is Restricted by the Receptor-Binding Region of VP4

The rotavirus species A (RVA) capsid contains the spike protein VP4, which interacts with VP6 and VP7 and is involved in cellular receptor binding. The capsid encloses the genome consisting of eleven dsRNA segments. Reassortment events can result in novel strains with changed properties. Using a plasmid-based reverse genetics system based on simian RVA strain SA11, we previously showed that the rescue of viable reassortants containing a heterologous VP4-encoding genome segment was strain-dependent. In order to unravel the reasons for the reassortment restrictions, we designed here a series of plasmids encoding chimeric VP4s. Exchange of the VP4 domains interacting with VP6 and VP7 was not sufficient for rescue of viable viruses. In contrast, the exchange of fragments encoding the receptor-binding region of VP4 resulted in virus rescue. All parent strains and the rescued reassortants replicated efficiently in MA-104 cells used for virus propagation. In contrast, replication in BSR T7/5 cells used for plasmid transfection was only efficient for the SA11 strain, whereas the rescued reassortants replicated slowly, and the parent strains failing to produce reassortants did not replicate. While future research in this area is necessary, replication in BSR T7/5 cells may be one factor that affects the rescue of RVAs.

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

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

Comparative genomic analysis of genogroup 1 and genogroup 2 rotaviruses circulating in seven US cities, 2014-2016

For over a decade, the New Vaccine Surveillance Network (NVSN) has conducted active rotavirus (RVA) strain surveillance in the USA. The evolution of RVA in the post-vaccine introduction era and the possible effects of vaccine pressure on contemporary circulating strains in the USA are still under investigation. Here, we report the whole-gene characterization (eleven ORFs) for 157 RVA strains collected at seven NVSN sites during the 2014 through 2016 seasons. The sequenced strains included 52 G1P[8], 47 G12P[8], 18 G9P[8], 24 G2P[4], 5 G3P[6], as well as 7 vaccine strains, a single mixed strain (G9G12P[8]), and 3 less common strains. The majority of the single and mixed strains possessed a Wa-like backbone with consensus genotype constellation of G1/G3/G9/G12-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1, while the G2P[4], G3P[6], and G2P[8] strains displayed a DS-1-like genetic backbone with consensus constellation of G2/G3-P[4]/P[6]/P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Two intergenogroup reassortant G1P[8] strains were detected that appear to be progenies of reassortment events between Wa-like G1P[8] and DS-1-like G2P[4] strains. Two Rotarix^® vaccine (RV1) and two RV5 derived (vd) reassortant strains were detected. Phylogenetic and similarity matrices analysis revealed 2-11 sub-genotypic allelic clusters among the genes of Wa- and DS-1-like strains. Most study strains clustered into previously defined alleles. Amino acid (AA) substitutions occurring in the neutralization epitopes of the VP7 and VP4 proteins characterized in this study were mostly neutral in nature, suggesting that these RVA proteins were possibly under strong negative or purifying selection in order to maintain competent and actual functionality, but fourteen radical (AA changes that occur between groups) AA substitutions were noted that may allow RVA strains to gain a selective advantage through immune escape. The tracking of RVA strains at the sub-genotypic allele constellation level will enhance our understanding of RVA evolution under vaccine pressure, help identify possible mechanisms of immune escape, and provide valuable information for formulation of future RVA vaccines.

Reverse Genetics Approach for Developing Rotavirus Vaccine Candidates Carrying VP4 and VP7 Genes Cloned from Clinical Isolates of Human Rotavirus

Species A rotaviruses (RVs) are a leading cause of severe acute gastroenteritis in infants and children younger than 5 years. Currently available RV vaccines were adapted from wild-type RV strains by serial passage of cultured cells or by reassortment between human and animal RV strains. These traditional methods require large-scale screening and genotyping to obtain vaccine candidates. Reverse genetics is a tractable, rapid, and reproducible approach to generating recombinant RV vaccine candidates carrying any VP4 and VP7 genes that provide selected antigenicity. Here, we developed a vaccine platform by generating recombinant RVs carrying VP4 (P[4] and P[8]), VP7 (G1, G2, G3, G8, and G9), and/or VP6 genes cloned from human RV clinical samples using the simian RV SA11 strain (G3P[2]) as a backbone. Neutralization assays using monoclonal antibodies and murine antisera revealed that recombinant VP4 and VP7 monoreassortant viruses exhibited altered antigenicity. However, replication of VP4 monoreassortant viruses was severely impaired. Generation of recombinant RVs harboring a chimeric VP4 protein for SA11 and human RV gene components revealed that the VP8* fragment was responsible for efficient infectivity of recombinant RVs. Although this system must be improved because the yield of vaccine viruses directly affects vaccine manufacturing costs, reverse genetics requires less time than traditional methods and enables rapid production of safe and effective vaccine candidates.IMPORTANCE Although vaccines have reduced global RV-associated hospitalization and mortality over the past decade, the multisegmented genome of RVs allows reassortment of VP4 and VP7 genes from different RV species and strains. The evolutionary dynamics of novel RV genotypes and their constellations have led to great genomic and antigenic diversity. The reverse genetics system is a powerful tool for manipulating RV genes, thereby controlling viral antigenicity, growth capacity, and pathogenicity. Here, we generated recombinant simian RVs (strain SA11) carrying heterologous VP4 and VP7 genes cloned from clinical isolates and showed that VP4- or VP7-substituted chimeric viruses can be used for antigenic characterization of RV outer capsid proteins and as improved seed viruses for vaccine production.

Phylogeographic diversity and hybrid zone of Hantaan orthohantavirus collected in Gangwon Province, Republic of Korea

Background

Hantaan orthohantavirus (Hantaan virus, HTNV), harbored by Apodemus agrarius (the striped field mouse), causes hemorrhagic fever with renal syndrome (HFRS) in humans. Viral genome-based surveillance at new expansion sites to identify HFRS risks plays a critical role in tracking the infection source of orthohantavirus outbreak. In the Republic of Korea (ROK), most studies demonstrated the serological prevalence and genetic diversity of orthohantaviruses collected from HFRS patients or rodents in Gyeonggi Province. Gangwon Province is a HFRS-endemic area with a high incidence of patients and prevalence of infected rodents, ROK. However, the continued epidemiology and surveillance of orthohantavirus remain to be investigated.

Methodology/Principal findings

Whole-genome sequencing of HTNV was accomplished in small mammals collected in Gangwon Province during 2015–2018 by multiplex PCR-based next-generation sequencing. To elucidate the geographic distribution and molecular diversity of viruses, we conducted phylogenetic analyses of HTNV tripartite genomes. We inferred the hybrid zone using cline analysis to estimate the geographic contact between two different HTNV lineages in the ROK. The graph incompatibility based reassortment finder performed reassortment analysis. A total of 12 HTNV genome sequences were completely obtained from A. agrarius newly collected in Gangwon Province. The phylogenetic and cline analyses demonstrated the genetic diversity and hybrid zone of HTNV in the ROK. Genetic exchange analysis suggested the possibility of reassortments in Cheorwon-gun, a highly HFRS-endemic area.

Conclusions/Significance

The prevalence and distribution of HTNV in HFRS-endemic areas of Gangwon Province enhanced the phylogeographic map for orthohantavirus outbreak monitoring in ROK. This study revealed the hybrid zone reflecting the genetic diversity and evolutionary dynamics of HTNV circulating in Gangwon Province. The results arise awareness of rodent-borne orthohantavirus diseases for physicians in the endemic area of ROK.

The genetic and molecular epidemiological studies on small mammals derived from hemorrhagic fever with renal syndrome (HFRS)-endemic areas have consistently conducted for the public health surveillance and mitigation of orthohantavirus outbreak in the Republic of Korea (ROK). Implementing viral genome-based surveillance at new expansion sites that may identify HFRS risks is critical for tracking the location of orthohantavirus infections and diagnosing HFRS. In the present study, whole-genome sequences of Hantaan virus (HTNV) from small mammals in Gangwon Province were recovered using multiplex PCR-based next-generation sequencing during 2015–2018. In HFRS-endemic regions including Cheorwon-gun, Chuncheon-si, and Hwacheon-gun, additional HTNV genome sequences contribute to establish a high-resolution phylogeographic map for tracking emerging orthohantavirus-induced diseases. The cline analysis revealed a remarkable hybrid zone by showing spatial contact regions of HTNV at two sites (Cheorwon-gun and Hwacheon-gun) and the spatial separation and sequence divergence across genome segments of HTNV in Gangwon Province. These results demonstrate the genetic diversity and hybrid zone of HTNV circulating in Gangwon Province. These findings increase an awareness raising about HFRS in the endemic area of ROK.

Genotype constellations of the rotavirus A strains circulating in Nizhny Novgorod, Russia, 2017-2018

Currently, the full-genome-based classification is widely used to investigate rotavirus A (RVA) strains found in different countries around the world. However, the information on the full genotypes of rotaviruses circulating in Russia is limited. Using partial sequencing, this study determined the full genotype constellations of 15 RVA strains in total commonly detected in Nizhny Novgorod (European part of Russia) in 2017-2018, three from each of the following genotypes G1P[8], G4P[8], and G9P[8] and six from G2P[4]. There were two intergenogroup mono-reassortants possessing an identical genotype constellation of G4-P[8]-I1-R1-C1-M1-A1-N1-T1-E2-H1 with the DS-1-like NSP4 gene of probably local origin. A variety of subgenotype lineages and their combinations of Wa-like rotaviruses and genetic heterogeneity among G9P[8] and G1P[8] strains were shown on the basis of phylogenetic analysis of each gene. Moreover, two distinct co-circulating variants that differed in all 11 genome segments were found among DS-1-like rotaviruses.

Rapid generation of rotavirus single-gene reassortants by means of eleven plasmid-only based reverse genetics

Reassortment is an important mechanism in the evolution of group A rotaviruses (RVAs), yielding viruses with novel genetic and phenotypic traits. The classical methods for generating RVA reassortants with the desired genetic combinations are laborious and time-consuming because of the screening and selection processes required to isolate a desired reassortant. Taking advantage of a recently developed RVA reverse genetics system based on just 11 cloned cDNAs encoding the RVA genome (11 plasmid-only system), we prepared a panel of simian SA11-L2 virus-based single-gene reassortants, each containing 1 segment derived from human KU virus of the G1P[8] genotype. It was shown that there was no gene-specific restriction of the reassortment potential. In addition to these 11 single-gene reassortants, a triple-gene reassortant with KU-derived core-encoding VP1–3 gene segments with the SA11-L2 genetic background, which make up a virion composed of the KU-based core, and SA11-L2-based intermediate and outer layers, could also be prepared with the 11 plasmid-only system. Finally, for possible clinical application of this system, we generated a series of VP7 reassortants representing all the major human RVA G genotypes (G1–4, G9 and G12) efficiently. The preparation of each of these single-gene reassortants was achieved within just 2 weeks. Our results demonstrate that the 11 plasmid-only system allows the rapid and reliable generation of RVA single-gene reassortants, which will be useful for basic research and clinical applications.

Leveraging a Validated in silico Approach to Elucidate Genotype-Specific VP7 Epitopes and Antigenic Relationships of Porcine Rotavirus A

Rotavirus A (RVA) remains one of the most widespread causes of diarrheal disease and mortality in piglets despite decades of research and efforts to boost lactogenic immunity for passive protection. Genetic changes at B cell epitopes (BCEs) may be driving failure of lactogenic immunity, which relies on production of IgA antibodies to passively neutralize RVA within the piglet gut, yet little research has mapped epitopes to swine-specific strains of RVA. Here we describe a bioinformatic approach to predict BCEs on the VP7 outer capsid protein using sequence data alone. We first validated the approach using a previously published dataset of VP7-specific cross-neutralization titers, and found that amino acid changes at predicted BCEs on the VP7 protein allowed for accurate recapitulation of antigenic relationships among the strains. Applying the approach to a dataset of swine RVA sequences identified 9 of the 11 known BCEs previously mapped to swine strains, indicating that epitope prediction can identify sites that are known to drive neutralization escape in vitro. Additional genotype-specific BCEs were also predicted that may be the cause of antigenic differences among strains of RVA on farms and should be targeted for further confirmatory work. The results of this work lay the groundwork for high throughput, immunologically-relevant analysis of swine RVA sequence data, and provide potential sites that can be targeted with vaccines to reduce piglet mortality and support farm health.

Genetic Characterizations and Molecular Evolution of VP7 Gene in Human Group A Rotavirus G1

Rotavirus group A (RVA) G1 is one leading genotype circulating in humans worldwide, and related molecular information from a global perspective is still limited. Here, we present a comprehensive description of the genetic characterizations and molecular evolution of the RVA G1 VP7 gene. Our results show that RVA G1 can be divided into two lineages and multiple sub-lineages with a relatively high genetic diversity. Vaccine strains are phylogenetic, closer to lineage I. The evolutionary rate of the RVA G1 VP7 gene is 8.869 × 10^-4 substitutions/site/year, and its most recent common ancestor was in 1933. The RVA G1 VP7 gene shows a linear evolution at the nucleotide level and a linear accumulation of difference at the amino acid level. Sub-lineage replacement of G1 VP7 gene is also observed and the effective population size of the G1 VP7 gene has had great change in the past decades and has remained stable in recent years. Altogether, the RVA G1 VP7 gene constantly evolves and there is no clear evidence that the evolution of the RVA G1 VP7 gene was influenced by vaccines. Continuous surveillance is still indispensable to evaluate the molecular epidemiology of RVA, especially in the post-vaccination era.

First isolation and whole-genome characterization of a G9P[14] rotavirus strain from a diarrheic child in Egypt

An unusual group A rotavirus (RVA) strain (RVA/Human-tc/EGY/AS997/2012/G9[14]) was isolated for the first time in a faecal sample from a 6-month-old child who was hospitalized for treatment of acute gastroenteritis in Egypt in 2012. Whole-genome analysis showed that the strain AS997 had a unique genotype constellation: G9-P[14]-I2-R2-C2-M2-A11-N2-T1-E2-H1. Phylogenetic analysis indicated that the strain AS997 had the consensus P[14] genotype constellation with the G9, T1 and H1 reassortment. This suggests either a mixed gene configuration originated from a human Wa-like strain and a P[14]-containing animal virus, or that this P[14] could have been acquired via reassortment of human strains only. The study shows the possible roles of interspecies transmission and multiple reassortment events leading to the generation of novel rotavirus genotypes and underlines the importance of whole-genome characterization of rotavirus strains in surveillance studies.

Comparative Sequence Analysis of Historic and Current Porcine Rotavirus C Strains and Their Pathogenesis in 3-Day-Old and 3-Week-Old Piglets

The increased prevalence of porcine group C rotavirus (PRVC) in suckling piglets and the emergence of new genetically distinct PRVC strains are concerning due to the associated significant economic losses they cause to the swine industry. We sequenced and analyzed two new PRVC strains, RV0104 (G3), and RV0143 (G6) and compared their pathogenesis with that of the historic strain Cowden (G1) in gnotobiotic (Gn) pigs. Near complete genome sequence analysis confirmed that these two strains were distinct from one another and the Cowden strain. VP1, VP2, VP6, NSP1-NSP3, and NSP5 genes were more similar between Cowden and RV0143, whereas VP3, VP7, and NSP4 shared higher nucleotide identity between Cowden and RV0104. Three-day-old and 3-week-old Gn piglets were inoculated with 10⁵ FFU/piglet of Cowden, RV0104 or RV0143, or mock. All 3-day-old piglets developed severe diarrhea, anorexia, and lethargy, with mean PRVC fecal shedding titers peaking and numerically higher in RV0104 and RV0143 piglets on post infection day (PID) 2. Histopathological examination of the small intestine revealed that the 3-day-old Cowden and RV0104 inoculated piglets were mildly affected, while significant destruction of small intestinal villi was observed in the RV0143 inoculated piglets. Consistent with the highest degree of pathological changes in the small intestines, the RV0143 inoculated piglets had numerically higher levels of serum IL-17 and IFN-α cytokines and numerically lower PRVC IgA geometric mean antibody titers. Milder pathological changes and overall higher titers of PRVC IgA antibodies were observed in 3-week-old vs. 3-day-old piglets. Additionally, diarrhea was only observed in RV0104 and RV0143 (but not Cowden) inoculated 3-week-old piglets, while levels of serum IL-10 and PRVC IgA antibodies were higher in Cowden inoculated pigs, consistent with the lack of diarrhea. Thus, we confirmed that these current, genetically heterogeneous PRVC strains possess distinct pathobiological characteristics that may contribute to the increased prevalence of PRVC diarrhea in neonatal suckling piglets.

Full genome characterization of novel DS-1-like G9P[8] rotavirus strains that have emerged in Thailand

The emergence and rapid spread of unusual DS-1-like intergenogroup reassortant rotaviruses having G1/3/8 genotypes have been recently reported from major parts of the world (Africa, Asia, Australia, Europe, and the Americas). During rotavirus surveillance in Thailand, three novel intergenogroup reassortant strains possessing the G9P[8] genotype (DBM2017-016, DBM2017-203, and DBM2018-291) were identified in three stool specimens from diarrheic children. In the present study, we determined and analyzed the full genomes of these three strains. On full-genomic analysis, all three strains were found to share a unique genotype constellation comprising both genogroup 1 and 2 genes: G9-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Phylogenetic analysis demonstrated that each of the 11 genes of the three strains was closely related to that of emerging DS-1-like intergenogroup reassortant, human, and/or locally circulating human strains. Thus, the three strains were suggested to be multiple reassortants that had acquired the G9-VP7 genes from co-circulating Wa-like G9P[8] rotaviruses in the genetic background of DS-1-like intergenogroup reassortant (likely equine-like G3P[8]) strains. To our knowledge, this is the first description of emerging DS-1-like intergenogroup reassortant strains having the G9P[8] genotype. Our observations will add to the growing insights into the dynamic evolution of emerging DS-1-like intergenogroup reassortant rotaviruses through reassortment.

Reverse Genetics Reveals a Role of Rotavirus VP3 Phosphodiesterase Activity in Inhibiting RNase L Signaling and Contributing to Intestinal Viral Replication In Vivo

Our understanding of how rotavirus (RV) subverts host innate immune signaling has greatly increased over the past decade. However, the relative contribution of each virus-encoded innate immune antagonist has not been fully studied in the context of RV infection in vivo Here, we present both in vitro and in vivo evidence that the host interferon (IFN)-inducible 2'-5'-oligoadenylate synthetase (OAS) and RNase L pathway effectively suppresses the replication of heterologous RV strains. VP3 from homologous RVs relies on its 2'-5'-phosphodiesterase (PDE) domain to counteract RNase L-mediated antiviral signaling. Using an RV reverse-genetics system, we show that compared to the parental strain, VP3 PDE mutant RVs replicated at low levels in the small intestine and were shed less in the feces of wild-type mice, and such defects were rescued in Rnasel-/- suckling mice. Collectively, these findings highlight an important role of VP3 in promoting viral replication and pathogenesis in vivo in addition to its well-characterized function as the viral RNA-capping enzyme.IMPORTANCE Rotaviruses are significant human pathogens that result in diarrhea, dehydration, and deaths in many children around the world. Rotavirus vaccines have suboptimal efficacy in low- to middle-income countries, where the burden of the diseases is the most severe. With the ultimate goal of improving current vaccines, we aim to better understand how rotavirus interacts with the host innate immune system in the small intestine. Here, we demonstrate that interferon-activated RNase L signaling blocks rotavirus replication in a strain-specific manner. In addition, virus-encoded VP3 antagonizes RNase L activity both in vitro and in vivo These studies highlight an ever-evolving arms race between antiviral factors and viral pathogens and provide a new means of targeted attenuation for next-generation rotavirus vaccine design.

Genetic Diversity of Human Rotavirus A Among Hospitalized Children Under-5 Years in Lebanon

Human rotavirus remains a major cause of gastroenteritis worldwide despite the availability of effective vaccines. In this study, we investigated the genetic diversity of rotaviruses circulating in Lebanon. We genetically characterized the VP4 and VP7 genes encoding the outer capsid proteins of 132 rotavirus-associated gastroenteritis specimens, previously identified in hospitalized children (<5 years) from 2011 to 2013 in Lebanon. These included 43 vaccine-breakthrough specimens and the remainder were from non-vaccinated subjects. Phylogenetic analysis of VP4 and VP7 genes revealed distinct clustering compared to the vaccine strains, and several substitutions were identified in the antigenic epitopes of Lebanese specimens. No unique changes were identified in the breakthrough specimens compared to non-breakthroughs that could explain the occurrence of infection in vaccinated children. Further, we report the emergence of a rare P[8] OP354-like strain with a G9 VP7 in Lebanon, possessing high genetic variability in their VP4 compared to vaccine strains. Therefore, human rotavirus strains circulating in Lebanon and globally have accumulated numerous substitutions in their antigenic sites compared to those currently used in the licensed vaccines. The successful spread and continued genetic drift of these strains over time might undermine the effectiveness of the vaccines. The effect of such changes in the antigenic sites on vaccine efficacy remains to be assessed.

Genetic diversity of group A rotaviruses in Moscow in 2018-2019

Here, we present the results of a study in which 639 samples obtained between October 2018 and April 2019 from patients with symptoms of acute gastroenteritis were tested for the presence of a rotavirus infection. The antigen of group A rotavirus was detected in 160 samples (25% of those tested). To study the genetic diversity of group A rotavirus, RNA was isolated from the samples, and polymerase chain reaction combined with reverse transcription (RT-PCR) with primers specific for the VP4, VP6, and VP7 genes of group A rotaviruses was performed. At least one fragment of the group A rotavirus genome was found in 101 samples (15.8%). These fragments were sequenced, and their G and P genotypes-as well as their combinations-were determined. The predominant G genotypes were G9 (35.8% of all genotyped samples) and G4 (28.4%), but the rare G12 genotype was also found (3.0%). The dominant P genotype was P[8]. The spectrum of certain G/P combinations of genotypes included seven variants. The most common variants were G9P[8] (37.2%) and G4P[8] (30.2%).

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.

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.