Intragenic recombination influences rotavirus diversity and evolution

Isolation, identification and genetic variation analysis of avian orthoreovirus in commercial broilers in China from 2016 to 2021

In the last decade, the emergence of variant strains of avian orthoreovirus (ARV) has caused an enormous economic impact on the poultry industry across China and other countries. This study aimed to evaluate the molecular evolution of the ARV lineages detected in Chinese commercial broiler farms. Firstly, ARV isolation and identification of commercial broiler arthritis cases from different provinces in China from 2016 to 2021 were conducted. A total of 51 pure ARV isolates were obtained. Sequencing results showed that there were five genotypes of the strains isolated in this study, of which genotype 1 ARV predominated, accounting for 56.9% (29/51). The whole gene sequences of 19 ARV representative isolates were successfully obtained. The genetic evolution analysis of 10 genome segments of 19 ARV isolates showed that the σC-encoding gene had evolved into six different lineages, while the other genome segments only differentiated into two to four different lineages. The results of recombination analysis showed that recombination events were present in the L3, M1 and S1 genome segments. Analysis of the variation of the key factor σC protein showed that the nucleotide and amino acid homologies of the σC were low among the different genotypes. Three-dimensional structural visualization analysis showed that all the structural changes of σC protein were concentrated in the spherical domain at the C-terminal, which is associated with host receptor binding.

Epidemiological, molecular, and evolutionary characteristics of G1P[8] rotavirus in China on the eve of RotaTeq application

Introduction

This study, conducted in China prior to RotaTeq's launch, examined the epidemiological, molecular, and evolutionary features of the G1P[8] genotype RVA in children admitted with diarrhea, to aid in evaluating its efficacy and impact on G1P[8] RVA in China.

Methods

Data from the Chinese viral diarrhea surveillance network were collected from January 2016 to December 2018. RVA strains identified as the G1P[8] genotype were subjected to whole-genome sequencing. Neutralizing epitope, amino acid selection pressure, and evolution dynamics analyses on VP7 and VP4 were performed using BioEdit v.7.0.9.0 and PyMOL v.2.5.2, four algorithms (MEME, SLAC, FEL, and FUBAR) in the Datamonkey online software, and the MCMC model in BEAST v. 1.10.4, respectively. Phylogenetic and identity features of 11 genes were assessed by DNAStar and MEGA v.7.

Results

Results showed that the detection rate of G1P[8] in China from 2016 to 2018 was generally low with significant seasonality. The whole genome of G1P[8] of four 2016 childhood diarrhea specimens was successfully sequenced. Phylogenetic and neutralizing epitope analysis showed that Rotavin-M1 might have better protection on G1P[8] prevalent in China than Rotarix and RotaTeq. Two conserved N-glycosylation sites on VP7 of Chinese G1P[8] might affect the protective effect of the vaccine. Evolution rate and selection pressure analysis identified the possibility of rapidly evolving and adapting to the new environment introduced by vaccines of G1P[8], whereas positive selection specific to VP4 indicated the potential tendency to select for dominant traits. Identity and phylogeny analysis showed that Chinese G1P[8] from before 2018 was generally stable with possible genetic recombination among local strains.

Discussion

These findings not only are of great significance for predicting the prevalence of G1P [8] in China, but also provide data reference for evaluating rotavirus vaccine efficacy.

Full genome sequence analysis of the predominant and uncommon G9P[4] rotavirus strains circulating in Tehran, Iran, 2021-2022: Evidence for inter and intra-genotype recombination

Group A rotaviruses (RVAs) are a major cause of acute gastroenteritis in children under 5 years of age worldwide. Herein, the genetic sequences of 11 RNA segments from three uncommon G9P[4] RVA strains found in the stool samples of children under 5 years of age in Iran were analyzed using next-generation sequencing (NGS) technology. The genomic constellations of these three uncommon G9P[4] strains indicated the presence of the double and quadruple reassortants of two G9P[4] strains, containing the VP7/NSP2 and VP7/VP2/NSP2/NSP4 genes on a DS-1-like genetic background, respectively. The genome of one strain indicated a Wa-like genetic backbone in a single-reassortant with the VP4 of the DS1-like human strains. With the exception of VP1, VP2, VP7, NSP2, NSP3, and NSP4 genes, which clustered with RVA of human origins belonging to cognate gene sequences of genogroup 1/2 genotypes/lineages, the remaining five genes (VP8/VP4, VP3, VP6, NSP1, NSP5) displayed direct evidence of recombination. It is presumed that the presence of uncommon G9P[4] strains in Iran is not linked to vaccination pressure, but rather to the high prevalence of RVA co-infection or the direct import of these uncommon RVA reassortants strains from other countries (especially those that have implemented RV vaccination).

Isolation, Genomic Characterization and Evolution of Six Porcine Rotavirus A Strains in a Pig Farming Group

Porcine rotavirus (PoRV) is a significant enteric pathogen causing gastroenteritis in piglets, which causes huge economic loss to the Chinese pig industry. In this study, six porcine rotavirus A strains were isolated from three adjacent sow farms belonging to the same company within one year, which suffered severe diarrhea outbreaks. AHBZ2303 (G11P[7]) and AHBZ2305 (G9P[23]), AHBZ2304 (G9P[23]) and AHBZ2312 (G4P[6]), AHBZ2310 (G9P[23]) and AHBZ2402 (G5P[23]) were isolated from Farm A, B and C, respectively. All six isolates were related to human rotavirus through complete genome analysis, suggesting the potential cross-species infection between humans and pigs. Evolutionary analysis revealed that AHBZ2303 and AHBZ2304 likely emerged simultaneously in Farm A and B, and then AHBZ2304 was introduced to Farm A and C, leading to the emergence of AHBZ2305 and AHBZ2310. Recombination and large variation were identified for AHBZ2312 and AHBZ2402. These findings provided insights into the transmission and evolution of PoRV among farms and underscored the need for enhanced monitoring to mitigate the risk of outbreaks from novel variants.

Complete genome constellation of a dominant Bovine rotavirus genotype circulating in Bangladesh reveals NSP4 intragenic recombination with human strains

Rotavirus A is a leading cause of non-bacterial gastroenteritis in humans and domesticated animals. Despite the vast diversity of bovine Rotavirus A strains documented in South Asian countries, there are very few whole genomes available for phylogenetic study. A cross-sectional study identified a high prevalence of the G6P[11] genotype of bovine Rotavirus A circulating in the commercial cattle population in Bangladesh. Next-generation sequencing and downstream phylogenetic analysis unveiled all 11 complete gene segments of this strain (BD_ROTA_CVASU), classifying it under the genomic constellation G6P[11]-I2-R2-C2-M2-A13-N2-T6-E2-H3, which belongs to a classical DS-1-like genomic backbone. We found strong evidence of intragenic recombination between human and bovine strains in the Non-structural protein 4 (NSP4) gene, which encodes a multifunctional enterotoxin. Our analyses highlight frequent zoonotic transmissions of rotaviruses in diverse human-animal interfaces, which might have contributed to the evolution and pathogenesis of this dominant genotype circulating in the commercial cattle population in Bangladesh.

Evolution of DS-1-like G8P[8] rotavirus A strains from Vietnamese children with acute gastroenteritis (2014-21): Adaptation and loss of animal rotavirus-derived genes during human-to-human spread

Animal rotaviruses A (RVAs) are considered the source of emerging, novel RVA strains that have the potential to cause global spread in humans. A case in point was the emergence of G8 bovine RVA consisting of the P[8] VP4 gene and the DS-1-like backbone genes that appeared to have jumped into humans recently. However, it was not well documented what evolutionary changes occurred on the animal RVA-derived genes during circulation in humans. Rotavirus surveillance in Vietnam found that DS-1-like G8P[8] strains emerged in 2014, circulated in two prevalent waves, and disappeared in 2021. This surveillance provided us with a unique opportunity to investigate the whole process of evolutionary changes, which occurred in an animal RVA that had jumped the host species barrier. Of the 843 G8P[8] samples collected from children with acute diarrhoea in Vietnam between 2014 and 2021, fifty-eight strains were selected based on their distinctive electropherotypes of the genomic RNA identified using polyacrylamide gel electrophoresis. Whole-genome sequence analysis of those fifty-eight strains showed that the strains dominant during the first wave of prevalence (2014-17) carried animal RVA-derived VP1, NSP2, and NSP4 genes. However, the strains from the second wave of prevalence (2018-21) lost these genes, which were replaced with cognate human RVA-derived genes, thus creating strain with G8P[8] on a fully DS-1-like human RVA gene backbone. The G8 VP7 and P[8] VP4 genes underwent some point mutations but the phylogenetic lineages to which they belonged remained unchanged. We, therefore, propose a hypothesis regarding the tendency for the animal RVA-derived genes to be expelled from the backbone genes of the progeny strains after crossing the host species barrier. This study underlines the importance of long-term surveillance of circulating wild-type strains in order to better understand the adaptation process and the fate of newly emerging, animal-derived RVA among the human population. Further studies are warranted to disclose the molecular mechanisms by which spillover animal RVAs become readily transmissible among humans, and the roles played by the expulsion of animal-derived genes and herd immunity formed in the local population.

Genetic characterization and evidence for multiple reassortments of rotavirus A G3P[3] in dogs and cats in Thailand

Rotavirus A (RVA) causes gastroenteritis in humans and animals. The zoonotic potential of RVA has been reported and raises major concerns, especially in animal-human interface settings. The study aimed to characterize and investigate the genetic diversity among RVAs in dogs and cats in Thailand. We collected 572 rectal swab samples from dogs and cats in Bangkok animal hospitals from January 2020 to June 2021. The one-step RT-PCR assay detected RVAs in 1.92% (11/572) of the samples, with 2.75% (8/290) in dogs and 1.06% (3/282) in cats. Two canine RVA and one feline RVA were subjected to whole genome sequencing. Our results showed that all three viruses were identified as RVA genotype G3P[3]. The genetic constellation of RVAs is unique for different species. For canine RVAs is G3-P [3]-I3-R3-C3-M3-A9-N2-T3-E3-H6, while Feline RVA is G3-P [3]-I8-R3-C3-M3-A9-N3-T3-E3-H6. Notably, both canine and feline RVAs contained the AU-1 genetic constellation with multiple reassortments. The results of phylogenetic, genetic, and bootscan analyses showed that canine RVAs may have reassorted from dog, human, and cat RVAs. While feline RVA was closely related to RVAs in humans, bats, and simians. This study provided genetic characteristics and diversity of RVAs in dogs and cats and suggested possible multiple reassortments, suggesting the zoonotic potential of the viruses. Thus, public health awareness should be raised regarding the zoonotic potential of RVAs in dogs and cats. Further studies on RVAs on a larger scale in dogs and cats in Thailand are needed.

Cryo-EM structure of rotavirus B NSP2 reveals its unique tertiary architecture

Rotavirus (RV) NSP2 is a multifunctional RNA chaperone that exhibits numerous activities that are essential for replication and viral genome packaging. We performed an in silico analysis that highlighted a distant relationship of NSP2 from rotavirus B (RVB) to proteins from other human RVs. We solved a cryo-electron microscopy structure of RVB NSP2 that shows structural differences with corresponding proteins from other human RVs. Based on the structure, we identified amino acid residues that are involved in RNA interactions. Anisotropy titration experiments showed that these residues are important for nucleic acid binding. We also identified structural motifs that are conserved in all RV species. Collectively, our data complete the structural characterization of rotaviral NSP2 protein and demonstrate its structural diversity among RV species.IMPORTANCERotavirus B (RVB), also known as adult diarrhea rotavirus, has caused epidemics of severe diarrhea in China, India, and Bangladesh. Thousands of people are infected in a single RVB epidemic. However, information on this group of rotaviruses remains limited. As NSP2 is an essential protein in the viral life cycle, including its role in the formation of replication factories, it may be a target for future antiviral strategy against viruses with similar mechanisms.

Spatio-temporal spread and evolution of Lassa virus in West Africa

BACKGROUND

Lassa fever is a hemorrhagic disease caused by Lassa virus (LASV), which has been classified by the World Health Organization as one of the top infectious diseases requiring prioritized research. Previous studies have provided insights into the classification and geographic characteristics of LASV lineages. However, the factor of the distribution and evolution characteristics and phylodynamics of the virus was still limited.

METHODS

To enhance comprehensive understanding of LASV, we employed phylogenetic analysis, reassortment and recombination detection, and variation evaluation utilizing publicly available viral genome sequences.

RESULTS

The results showed the estimated the root of time of the most recent common ancestor (TMRCA) for large (L) segment was approximately 634 (95% HPD: [385879]), whereas the TMRCA for small (S) segment was around 1224 (95% HPD: [10301401]). LASV primarily spread from east to west in West Africa through two routes, and in route 2, the virus independently spread to surrounding countries through Liberia, resulting in a wider spread of LASV. From 1969 to 2018, the effective population size experienced two significant increased, indicating the enhanced genetic diversity of LASV. We also found the evolution rate of L segment was faster than S segment, further results showed zinc-binding protein had the fastest evolution rate. Reassortment events were detected in multiple lineages including sub-lineage IIg, while recombination events were observed within lineage V. Significant amino acid changes in the glycoprotein precursor of LASV were identified, demonstrating sequence diversity among lineages in LASV.

CONCLUSION

This study comprehensively elucidated the transmission and evolution of LASV in West Africa, providing detailed insights into reassortment events, recombination events, and amino acid variations.

Phylogenetic and Recombination Analysis of Clinical Vitreous Humor-Derived Adenovirus Isolates Reveals Discordance Between Serotype and Phylogeny

Purpose

To sequence, identify, and perform phylogenetic and recombination analysis on three clinical adenovirus samples taken from the vitreous humor at the Bascom Palmer Eye Institute.

Methods

The PacBio Sequel II was used to sequence the genomes of the three clinical adenovirus isolates. To identify the isolates, a full genome-based multiple sequence alignment (MSA) of 722 mastadenoviruses was generated using multiple alignment using fast Fourier transform (MAFFT). MAFFT was also used to generate genome-based human adenovirus B (HAdV-B) MSAs, as well as HAdV-B fiber, hexon, and penton protein-based MSAs. To examine recombination within HAdV-B, RF-Net 2 and Bootscan software programs were used.

Results

In the course of classifying three new atypical ocular adenovirus samples, taken from the vitreous humor, we found that all three isolates were HAdV-B species. The three Bascom Palmer HAdV-B genomes were then combined with over 300 HAdV-B genome sequences, including nine ocular HAdV-B genome sequences. Attempts to categorize the penton, hexon, and fiber serotypes using phylogeny of the three Bascom Palmer samples were inconclusive due to incongruence between serotype and phylogeny in the dataset. Recombination analysis using a subset of HAdV-B strains to generate a hybridization network detected recombination between nonhuman primate and human-derived strains, recombination between one HAdV-B strain and the HAdV-E outgroup, and limited recombination between the B1 and B2 clades.

Conclusions

The discordance between serotype and phylogeny detected in this study suggests that the current classification system does not accurately describe the natural history and phylogenetic relationships among adenoviruses.

Interspecies recombination in NSP3 gene in the first porcine rotavirus H in Russia identified using nanopore-based metagenomic sequencing

During the last decade, porcine rotavirus H was detected in the USA, Asian regions, South Africa, Brazil, and a couple of European countries. In the presented study, the virus was identified in piglets on a farrow-to-finish farm in Russia during metagenomic surveillance. Currently, it is the first identification of this species in the country. As a diagnostic method, nanopore-based metagenomic sequencing was applied. The obtained nanopore reads allowed for the assembly of 10 genome segments out of 11. The phylogenetic analysis revealed the virus belonged to the porcine cluster and had GX-P3-I3-R3-C3-M8-A7-N1-T5-E3-H3 genome constellation. Moreover, three potential new genotype groups for VP3, NSP1, and NSP3 genes were determined. Additionally, a recombination between RVH and RVC in the NSP3 gene was detected. The study provides significant information about a novel RVH strain.

The impact and complete genome characterisation of viruses involved in outbreaks of gastroenteritis in a farrow-to-finish holding

Viral enteric pathogens continuously burden intensive pig farming, causing gastrointestinal diseases of epidemic and endemic nature. The present study investigated two diarrhoea outbreaks on a large farrow-to-finish holding and subsequent circulation of outbreak-related enteric viruses. These viruses were characterised by whole genome sequencing, and statistical evaluation of the impact on specific production metrics was performed. The results provided evidence that the Porcine epidemic diarrhoea virus-swine enteric coronavirus (PEDV-SeCoV) S gene recombinant strain was responsible for the first outbreak, whilst Rotavirus A (RVA) in a mixed infection with Rotavirus B (RVB) and porcine kobuvirus (PKV) probably caused the second diarrhoea outbreak. Whole genome characterisation revealed a porcine origin of all viruses involved and significant heterogeneity of RVB strain, proposing four novel genotypes and changes in RVB VP1 genotype classification. The statistical evaluation confirmed only a minor disturbance in the number of weaned pigs per sow, with statistical forecasting showing positive trends. A follow-up study corroborated the endemicity of RVA and PKV, in contrast to PEDV-SeCoV. Punctual, comprehensive and timely investigation of diarrhoea outbreaks is a prerequisite for applying adequate pig health and biosecurity management. Calculating such outbreaks' impact on production metrics can potentially shape future decisions on management improvements.

Whole-genome characterization of common rotavirus strains circulating in Vellore, India from 2002 to 2017: emergence of non-classical genomic constellations

Rotaviruses (RVs) are the most common etiological agent of acute gastroenteritis among young children, even after vaccine introduction in low-income countries. A whole-genome classification representing the 11 RV genes, was introduced for surveillance and characterization of RVs. This study characterized the common circulating strains in Vellore, India from 2002 to 2017 to understand rotavirus strain diversity and evolution using Whole genome sequencing (WGS) carried out on Illumina MiSeq. The 89% (92% of Wa-like, 86% of DS-1-like) of strains had classical constellations, while reassortant constellations were seen in 11% (8% of Wa-like, 14% of DS-1-like) of the strains. The rare E6-NSP4 in combination with DS-1 like G1P[8] and the emergence of the OP-354 subtype of P[8] were identified. Phylogenetics of RV strains revealed multiple subtypes circulating in the past 15 years, with strong evidence of animal to human gene transmission among several strains.

Genome analyses of species A rotavirus isolated from various mammalian hosts in Northern Ireland during 2013-2016

Rotavirus group A (RVA) is the most important cause of acute diarrhoea and severe dehydration in young mammals. Infection in livestock is associated with significant mortality and economic losses and, together with wildlife reservoirs, acts as a potential source of zoonotic transmission. Therefore, molecular surveillance of circulating RVA strains in animal species is necessary to assess the risks posed to humans and their livestock. An RVA molecular epidemiological surveillance study on clinically diseased livestock species revealed high prevalence in cattle and pigs (31 per cent and 18 per cent, respectively) with significant phylogenetic diversity including a novel and divergent ovine artiodactyl DS-1-like constellation G10-P[15]-I2-R2-C2-M2-A11-N2-T6-E2-H3. An RVA gene reassortment occurred in an RVA asymptomatic pig and identified as a G5-P[13] strain, and a non-structural protein (NSP)2 gene had intergenomically reassorted with a human RVA strain (reverse zoonosis) and possessed a novel NSP4 enterotoxin E9 which may relate to the asymptomatic RVA infection. Analysis of a novel sheep G10-P[15] strain viral protein 4 gene imparts a putative homologous intergenic and interspecies recombination event, subsequently creating the new P[15] divergent lineage. While surveillance across a wider range of wildlife and exotic species identified generally negative or low prevalence, a novel RVA interspecies transmission in a non-indigenous pudu deer (zoo origin) with the constellation of G6-P[11]12-R2-C2-M2-A3-N2-T6-E2-H3 was detected at a viral load of 11.1 log10 copies/gram. The detection of novel emerging strains, interspecies reassortment, interspecies infection, and recombination of RVA circulating in animal livestock and wildlife reservoirs is of paramount importance to the RVA epidemiology and evolution for the One Health approach and post-human vaccine introduction era where highly virulent animal RVA genotypes have the potential to be zoonotically transmitted.

Evolutionary analysis of all eleven genes of species C rotaviruses circulating in humans and domestic animals

Species C rotaviruses (RVC) are the second most common rotavirus species known to cause gastroenteritis in humans and pigs and with occurrence documented in cattle, dogs, ferrets, and sloth bears. Despite the host-specific nature of RVC genotypes, cross-species transmission, reassortment, and recombination events are also documented. In the present study, we inferred the evolutionary history of globally circulating RVC strains, including time scale stasis, the most probable ancestral country, and the most probable source host using Bayesian methods implemented in BEAST v.1.8.4. The human-derived RVC strains were majorly monophyletic and further grouped into two lineages. The RVC strains derived from pigs were monophyletic for the VP1 and the remaining genes were classified into 2 to 4 groups based on the high posterior support. The root mean age for all the genes indicated the circulation of RVC for over 800 years. Overall, the time to Most Recent Common Ancestor of human RVC strains dated back to the beginning of the 20th century. The VP7 and NSP2 genes had the lowest rates of evolution compared to other genes. The majority of the genes of RVC showed their origin in Japan except for VP7 and VP4 genes in South Korea. The phylogeographic analysis with the country as a trait showed the role of Japan, China, and India in the dispersion of the virus. In the current study, significant transmission links between different hosts were analyzed for the first time using the host as a trait. Significant transmission links between pigs and other animal species as well as humans indicate possible transmission from the pig as a source host and suggest monitoring of proximity with animals.

Interspecies transmission of porcine-originated G4P[6] rotavirus A between pigs and humans: a synchronized spatiotemporal approach

As a leading viral cause of acute gastroenteritis in both humans and pigs, rotavirus A (RVA) poses a potential public health concern. Although zoonotic spillover of porcine RVA strains to humans is sporadic, it has been detected worldwide. The origin of chimeric human-animal strains of RVA is closely linked to the crucial role of mixed genotypes in driving reassortment and homologous recombination, which play a major role in shaping the genetic diversity of RVA. To better understand how genetically intertwined porcine and zoonotic human-derived G4P[6] RVA strains are, the present study employed a spatiotemporal approach to whole-genome characterization of RVA strains collected during three consecutive RVA seasons in Croatia (2018-2021). Notably, sampled children under 2 years of age and weanling piglets with diarrhea were included in the study. In addition to samples tested by real-time RT-PCR, genotyping of VP7 and VP4 gene segments was conducted. The unusual genotype combinations detected in the initial screening, including three human and three porcine G4P[6] strains, were subjected to next-generation sequencing, followed by phylogenetic analysis of all gene segments, and intragenic recombination analysis. Results showed a porcine or porcine-like origin for each of the eleven gene segments in all six RVA strains. The G4P[6] RVA strains detected in children most likely resulted from porcine-to-human interspecies transmission. Furthermore, the genetic diversity of Croatian porcine and porcine-like human G4P[6] strains was propelled by reassortment events between porcine and porcine-like human G4P[6] RVA strains, along with homologous intragenotype and intergenotype recombinations in VP4, NSP1, and NSP3 segments. Described concurrent spatiotemporal approach in investigating autochthonous human and animal RVA strains is essential in drawing relevant conclusions about their phylogeographical relationship. Therefore, continuous surveillance of RVA, following the One Health principles, may provide relevant data for assessing the impact on the protectiveness of currently available vaccines.

Genetic heterogeneity of group A rotaviruses: a review of the evolutionary dynamics and implication on vaccination

INTRODUCTION

Human rotavirus remains a major etiology of acute gastroenteritis among under 5-year children worldwide despite the availability of oral vaccines. The genetic instability of rotavirus and the ability to form different combinations from the different G- and P-types reshapes the antigenic landscape of emerging strains which often display limited or no antigen identities with the vaccine strain. As evidence also suggests, the selection of the antigenically distinct novel or rare strains and their successful spread in the human population has raised concerns regarding undermining the effectiveness of vaccination programs.

AREAS COVERED

We review aspects related to current knowledge about genetic and antigenic heterogeneity of rotavirus, the mechanism of genetic diversity and evolution, and the implication of genetic change on vaccination.

EXPERT OPINION

Genetic changes in the segmented genome of rotavirus can alter the antigenic landscape on the virion capsid and further promote viral fitness in a fully vaccinated population. Against this background, the potential risk of the appearance of new rotavirus strains over the long term would be better predicted by a continued and increased close monitoring of the variants across the globe to identify any change associated with disease dynamics.

Genetic diversity, reassortment, and recombination of mammalian orthoreoviruses from Japanese porcine fecal samples

Mammalian orthoreoviruses (MRVs) are non-enveloped double-stranded RNA viruses with a broad host range. MRVs are prevalent worldwide, and in Japan, they have been isolated from various hosts, including humans, dogs, cats, wild boars, and pigs, and they have also been found in sewage. However, Japanese porcine MRVs have not been genetically characterized. While investigating porcine enteric viruses including MRV, five MRVs were isolated from the feces of Japanese pigs using MA104 cell culture. Genetic analysis of the S1 gene revealed that the Japanese porcine MRV isolates could be classified as MRV-2 and MRV-3. Whole genome analysis showed that Japanese porcine MRVs exhibited genetic diversity, although they shared sequence similarity with porcine MRV sequences in the DDBJ/EMBL/GenBank database. Several potential intragenetic reassortment events were detected among MRV strains from pigs, sewage, and humans in Japan, suggesting zoonotic transmission. Furthermore, homologous recombination events were identified in the M1 and S1 genes of Japanese porcine MRV. These findings imply that different strains of Japanese porcine MRV share a porcine MRV genomic backbone and have evolved through intragenetic reassortment and homologous recombination events.

Complete Genome Sequencing Reveals Unusual Equine Rotavirus A of Bat Origin from India

Rotaviruses are the most common viral agents associated with foal diarrhea. Between 2014 and 2017, the annual prevalence of rotavirus in diarrheic foals ranged between 18 and 28% in Haryana (India). Whole-genome sequencing of two equine rotavirus A (ERVA) isolates (RVA/Horse-wt/IND/ERV4/2017 and RVA/Horse-wt/IND/ERV6/2017) was carried out to determine the genotypic constellations (GCs) of ERVAs. The GCs of both the isolates were G3-P[3]-I8-R3-C3-M3-A9-N3-T3-E3-H6, a unique combination reported for ERVAs so far. Both the isolates carried VP6 of genotype I8, previously unreported from equines. Upon comparison with RVAs of other species, the GC of both isolates was identical to that of a bat rotavirus strain, MSLH14, isolated from China in 2012. The nucleotide sequences of the genes encoding VP3, NSP2, and NSP3 shared >95.81% identity with bat RVA strains isolated from Africa (Gabon). The genes encoding VP1, VP2, VP7, NSP1, and NSP4 shared 94.82% to 97.12% nucleotide identities with the human strains which have zoonotic links to bats (RCH272 and MS2015-1-0001). The VP6 genes of both strains were distinct and had the highest similarity of only 87.08% with that of CMH222, a human strain of bat origin. The phylogenetic analysis and lineage studies revealed that VP7 of both isolates clustered in a new lineage (lineage X) of the G3 genotype with bat, human, and alpaca strains. Similarly, VP4 clustered in a distinct P[3] lineage. These unusual findings highlight the terra incognita of the genomic diversity of equine rotaviruses and support the need for the surveillance of RVAs in animals and humans with a "one health" approach. IMPORTANCE Rotaviruses are globally prevalent diarrheal pathogens in young animals including foals, piglets, calves, goats, sheep, cats, and dogs along with humans. The genome of rotaviruses consists of 11 segments, which enables them to undergo reshuffling by reassortment of segments from multiple species during mixed infections. In this study, the prevalence of equine rotaviruses was 32.11% in organized equine farms of North India. The complete genome analysis of two ERVA isolates revealed an unusual genomic constellation, which was previously reported only in a bat RVA strain. A segment-wise phylogenetic analysis revealed that most segments of both isolates were highly similar either to bat or to bat-like human rotaviruses. The occurrence of unusual bat-like rotaviruses in equines emphasizes the need of extensive surveillance of complete genomes of both animal and human rotaviruses with a "one health" approach.

Genomic diversity and intragenic recombination of species C rotaviruses

Rotavirus C (RVC) is a major cause of diarrhoea in swine, cattle, and humans worldwide. RVC exhibits sequence diversity in all 11 genes, especially in VP4 and VP7, and all segment-based genotyping has been performed similar to rotavirus A. To date, recombination events have been reported in rotavirus A and B. However, there are no reports describing gene recombination of RVC, except for recombination in NSP3 between RVC and rotavirus H. In this study, nine porcine RVC strains identified in Japanese pigs were completely sequenced and analysed together with RVC sequences from the GenBank database. The analyses showed that sequences of the VP4, VP2, and NSP1 of several porcine RVC strains did not branch with any of those of the RVC strains in the GenBank database, suggesting new genotypes. Several homologous recombination events, between or within genotypes, were identified in the VP4, VP7, VP2, NSP1, and NSP3 genes. Of these, nine, one, and one intergenotypic recombination events in the VP4, VP2, and NSP3 genes, respectively, were supported with sufficient statistical values. Although these findings suggest occurrences of the intragenic recombination events in the RVC genome, potential sequence errors and poor sequence assemblies in the databases should be watched with care. The results in this study present data about the important recombination events of the RVCs, which influence evolution of the virus by aiding them to gain genetic diversity and plasticity, although further sequence data will be necessary to obtain more comprehensive understanding of such mechanisms.

Effect of rotavirus genetic diversity on vaccine impact

Group A rotaviruses (RVAs) are the leading cause of gastroenteritis, causing 0.2 million deaths and several million hospitalisations globally each year. Four rotavirus vaccines (Rotarix^TM , RotaTeq^TM , Rotavac^® and ROTASIIL^® ) have been pre-qualified by the World Health Organization (WHO), but the two newly pre-qualified vaccines (Rotavac^® and ROTASIIL^® ) are currently only in use in Palestine and India, respectively. In 2009, WHO strongly proposed that rotavirus vaccines be included in the routine vaccination schedule of all countries around the world. By the end of 2019, a total of 108 countries had administered rotavirus vaccines, and 10 countries have currently been approved by Gavi for the introduction of rotavirus vaccine in the near future. With 39% of global coverage, rotavirus vaccines have had a substantial effect on diarrhoeal morbidity and mortality in different geographical areas, although efficacy appears to be higher in high income settings. Due to the segmented RNA genome, the pattern of RVA genotypes in the human population is evolving through interspecies transmission and/or reassortment events for which the vaccine might be less effective in the future. However, despite the relative increase in some particular genotypes after rotavirus vaccine use, the overall efficacy of rotavirus mass vaccination worldwide has not been affected. Some of the challenges to improve the effect of current rotavirus vaccines can be solved in the future by new rotavirus vaccines and by vaccines currently in progress.

Rotavirus research: 2014-2020

Rotaviruses are major causes of acute gastroenteritis in infants and young children worldwide and also cause disease in the young of many other mammalian and of avian species. During the recent 5-6 years rotavirus research has benefitted in a major way from the establishment of plasmid only-based reverse genetics systems, the creation of human and other mammalian intestinal enteroids, and from the wide application of structural biology (cryo-electron microscopy, cryo-EM tomography) and complementary biophysical approaches. All of these have permitted to gain new insights into structure-function relationships of rotaviruses and their interactions with the host. This review follows different stages of the viral replication cycle and summarizes highlights of structure-function studies of rotavirus-encoded proteins (both structural and non-structural), molecular mechanisms of viral replication including involvement of cellular proteins and lipids, the spectrum of viral genomic and antigenic diversity, progress in understanding of innate and acquired immune responses, and further developments of prevention of rotavirus-associated disease.

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.

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.

Unexpected Role of Rotavirus G3P[8] Infection in Causing Severe Diarrhea in a Major Tertiary Referral Hospital in the Prevaccine Era

Background: Rotavirus A species is associated with severe gastroenteritis in children. Rotavirus G1P[8] was the most prevalent genotype found in Kuwait in a study conducted between 2005 and 2006. The RotaTeq vaccine was included in the Kuwait national immunization program at the end of 2017. Objectives: Since there is no available data on the rotavirus genotypes circulating before the introduction of the vaccine, we conducted a study to investigate the role of rotaviruses in causing severe diarrhea in children hospitalized in a major tertiary referral hospital in Kuwait during the year 2016. Methods: Viral RNA was isolated from the stool samples of 101 children under five years of age, hospitalized for acute gastroenteritis. Rotavirus VP4 and VP7 dsRNA were detected by RT-PCR, and their partial sequences were analyzed by phylogenic analysis. Results: Rotavirus dsRNA was detected in 24.7% of children with median age of 1 year. The genotype G3P[8] accounted for 47% of cases, followed by G1P[8] (26%), G9P[8] (10.5%), G4P[8] (10.5%), and G9P[4] (5%). Only VP7 nucleotide sequences of rotavirus G3 or G4 type clustered in the same lineage as RotaTeq vaccine, while most VP4 nucleotide sequences of rotavirus P[8] type clustered in a different lineage compared to Rotarix and RotaTeq vaccines. Conclusions: Our findings highlight the role of rotavirus G3P[8] in causing severe diarrhea and invites future investigations to know whether the recent introduction of RotaTeq vaccine in Kuwait selects certain genotypes and subgenomic lineages.

Genomic analysis of group A rotavirus G12P[8] including a new Japanese strain revealed evidence for intergenotypic recombination in VP7 and VP4 genes

Group A rotavirus is a leading cause of severe acute gastroenteritis worldwide. In this study, the first complete coding sequences of 11 RNA segments of human group A rotavirus G12P[8] in Japan were determined by an unbiased viral metagenomics. Its genomic constellation (VP7-VP4-VP6-VP1-VP2-VP3-NSP1-NSP2-NSP3-NSP4-NSP5 genes) was identified as G12-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1. When performing the genetic analysis, we discovered an intergenotypic recombination event in the pig group A rotavirus G12P[8] strain BUW-14-A008. The novel recombination was found between two different genotypes G12 and G3 in the VP7 gene, and P[8] and P[13] in the VP4 gene.

Recombination and Positive Selection Differentially Shaped the Diversity of Betacoronavirus Subgenera

The Betacoronavirus genus of mammal-infecting viruses includes three subgenera (Sarbecovirus, Embecovirus, and Merbecovirus), in which most known human coronaviruses, including SARS-CoV-2, cluster. Coronaviruses are prone to host shifts, with recombination and positive selection possibly contributing to their high zoonotic potential. We analyzed the role of these two forces in the evolution of viruses belonging to the Betacoronavirus genus. The results showed that recombination has been pervasive during sarbecovirus evolution, and it is more widespread in this subgenus compared to the other two. In both sarbecoviruses and merbecoviruses, recombination hotspots are clearly observed. Conversely, positive selection was a less prominent force in sarbecoviruses compared to embecoviruses and merbecoviruses and targeted distinct genomic regions in the three subgenera, with S being the major target in sarbecoviruses alone. Overall, the results herein indicate that Betacoronavirus subgenera evolved along different trajectories, which might recapitulate their host preferences or reflect the origins of the presently available coronavirus sequences.

Whole Genome In-Silico Analysis of South African G1P[8] Rotavirus Strains Before and After Vaccine Introduction Over A Period of 14 Years

Rotavirus G1P[8] strains account for more than half of the group A rotavirus (RVA) infections in children under five years of age, globally. A total of 103 stool samples previously characterized as G1P[8] and collected seven years before and seven years after introducing the Rotarix® vaccine in South Africa were processed for whole-genome sequencing. All the strains analyzed had a Wa-like constellation (G1-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1). South African pre- and post-vaccine G1 strains were clustered in G1 lineage-I and II while the majority (84.2%) of the P[8] strains were grouped in P[8] lineage-III. Several amino acid sites across ten gene segments with the exception of VP7 were under positive selective pressure. Except for the N147D substitution in the antigenic site of eight post-vaccine G1 strains when compared to both Rotarix® and pre-vaccine strains, most of the amino acid substitutions in the antigenic regions of post-vaccine G1P[8] strains were already present during the pre-vaccine period. Therefore, Rotarix® did not appear to have an impact on the amino acid differences in the antigenic regions of South African post-vaccine G1P[8] strains. However, continued whole-genome surveillance of RVA strains to decipher genetic changes in the post-vaccine period remains imperative.