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

Human Rotaviruses of Multiple Genotypes Acquire Conserved VP4 Mutations during Serial Passage

Human rotaviruses exhibit limited tropism and replicate poorly in most cell lines. Attachment protein VP4 is a key rotavirus tropism determinant. Previous studies in which human rotaviruses were adapted to cultured cells identified mutations in VP4. However, most such studies were conducted using only a single human rotavirus genotype. In the current study, we serially passaged 50 human rotavirus clinical specimens representing five of the genotypes most frequently associated with severe human disease, each in triplicate, three to five times in primary monkey kidney cells then ten times in the MA104 monkey kidney cell line. From 13 of the 50 specimens, we obtained 25 rotavirus antigen-positive lineages representing all five genotypes, which tended to replicate more efficiently in MA104 cells at late versus early passage. We used Illumina next-generation sequencing and analysis to identify variants that arose during passage. In VP4, variants encoded 28 mutations that were conserved for all P[8] rotaviruses and 12 mutations that were conserved for all five genotypes. These findings suggest there may be a conserved mechanism of human rotavirus adaptation to MA104 cells. In the future, such a conserved adaptation mechanism could be exploited to study human rotavirus biology or efficiently manufacture vaccines.

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.

Genomic epidemiology of the rotavirus G2P[4] strains in coastal Kenya pre- and post-rotavirus vaccine introduction, 2012-8

The introduction of rotavirus vaccines into the national immunization programme in many countries has led to a decline in childhood diarrhoea disease burden. Coincidentally, the incidence of some rotavirus group A (RVA) genotypes has increased, which may result from non-vaccine-type replacement. Here, we investigate the evolutionary genomics of rotavirus G2P[4] which has shown an increase in countries that introduced the monovalent Rotarix® vaccine. We examined sixty-three RVA G2P[4] strains sampled from children (aged below 13 years) admitted to Kilifi County Hospital, coastal Kenya, pre- (2012 to June 2014) and post-(July 2014 to 2018) rotavirus vaccine introduction. All the sixty-three genome sequences showed a typical DS-1-like genome constellation (G2-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2). Pre-vaccine G2 sequences predominantly classified as sub-lineage IVa-3 and co-circulated with low numbers of sub-lineage IVa-1 strains, whereas post-vaccine G2 sequences mainly classified into sub-lineage IVa-3. In addition, in the pre-vaccine period, P[4] sub-lineage IVa strains co-circulated with low numbers of P[4] lineage II strains, but P[4] sub-lineage IVa strains predominated in the post-vaccine period. On the global phylogeny, the Kenyan pre- and post-vaccine G2P[4] strains clustered separately, suggesting that different virus populations circulated in the two periods. However, the strains from both periods exhibited conserved amino acid changes in the known antigenic epitopes, suggesting that replacement of the predominant G2P[4] cluster was unlikely a result of immune escape. Our findings demonstrate that the pre- and post-vaccine G2P[4] strains circulating in Kilifi, coastal Kenya, differed genetically but likely were antigenically similar. This information informs the discussion on the consequences of rotavirus vaccination on rotavirus diversity.

Genomic Analysis of G2P[4] Group A Rotaviruses in Zambia Reveals Positive Selection in Amino Acid Site 7 of Viral Protein 3

The G2P[4] genotype is among the rotavirus strains that circulate commonly in humans. Several countries have reported its immediate upsurge after the introduction of rotavirus vaccination, raising concern about sub-optimal vaccine effectiveness against this genotype in the long term. This study aimed to gain insight into the evolution of post-vaccine Zambian G2P[4] group A rotavirus (RVA) strains and their overall genetic make-up by analysis of sequence alignments at the amino acid (AA) level. Twenty-nine Zambian G2P[4] rotavirus strains were subjected to whole-genome sequencing using the Illumina MiSeq® platform. All the strains exhibited the typical DS-1-like genotype constellation, and the nucleotide sequences of the 11 genome segments showed high nucleotide similarities (>97%). Phylogenetic analyses together with representative global G2P[4] RVA showed that Zambian strains clustered into human lineages IV (for VP2, VP4, VP7, NSP1, and NSP5), V (for VP1, VP3, VP6, NSP2, and NSP3), and XXIII (for NSP4). The AA differences between the lineages where the study strains clustered and lineages of global reference strains were identified and analyzed. Selection pressure analysis revealed that AA site seven in the Viral Protein 3 (VP3) genome segment was under positive selection. This site occurs in the region of intrinsic disorder in the VP3 protein, and Zambian G2P[4] strains could potentially be utilizing this intrinsically disordered region to survive immune pressure. The Zambian G2P[4] strains from 2012 to 2016 comprised the G2P[4] strains that have been circulating globally since the early 2000s, highlighting the epidemiological fitness of these contemporary G2P[4] strains. Continuous whole-genome surveillance of G2P[4] strains remains imperative to understand their evolution during the post-vaccination period.

M. Eljaro Z, Abu-Awwad F, Bozya Y, Hjaija D, Bhat N, Leader T, Ramlawi A, Marzouqa H. Introducing ROTAVAC® to the occupied Palestinian Territories: Impact on diarrhea incidence, rotavirus prevalence and genotype composition

BACKGROUND

Rotavirus infection remains an important cause of morbidity and mortality in children. The introduction of vaccination programs in more than 100 countries has contributed to a decrease in hospitalizations and mortality. This study investigates the epidemiological impact of the rotavirus vaccine ROTAVAC® in the Palestinian Territories, the first country to switch from ROTARIX® to this new vaccine.

METHODS

Clinical surveillance data was collected fromchildren younger than 5attendingoutpatient clinics throughout Gaza withdiarrhea between 2015 and 2020. The incidence of all-cause diarrhea was assessed using an interrupted time-series approach. Rotavirus prevalence was determined at the Caritas Baby Hospital in the West Bank usingELISA on stool specimen of children younger than 5with diarrhea. Genotyping was performed on 325 randomly selected rotavirus-positive samples from January 2015 through December 2020 using multiplex PCR analysis.

RESULTS

Average monthly diarrhea casesdropped by 16.7% annually fromintroduction of rotavirus vaccination in May 2016 to the beginning of the SARS-CoV-2 epidemic in March 2020 for a total of 53%. Case count declines were maintained afterthe switchto ROTAVAC® in October 2018. Rotavirus positivity in stool samples declined by 67.1% over the same period without change followingthe switch to ROTAVAC®. The distribution of predominant genotypes in rotavirus-positive stool samples changed from a pre-vaccination G1P [8] to G9P[8] and G12P[8] during the ROTARIX® period and G2P[4] after the introduction of ROTAVAC®.

CONCLUSION

ROTAVAC® has shown epidemiological impact on par with ROTARIX® after its introduction to the national immunization schedule in the Palestinian Territories. A molecular genotype shift from a pre-vaccination predominance of G1P[8] to a current predominance of G2P[4] requires more long-term surveillance.

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.

Rotavirus Strain Trends in United States, 2009-2016: Results from the National Rotavirus Strain Surveillance System (NRSSS)

Before the introduction of vaccines, group A rotaviruses (RVA) were the leading cause of acute gastroenteritis in children worldwide. The National Rotavirus Strain Surveillance System (NRSSS) was established in 1996 by the Centers for Disease Control and Prevention (CDC) to perform passive RVA surveillance in the USA. We report the distribution of RVA genotypes collected through NRSSS during the 2009-2016 RVA seasons and retrospectively examine the genotypes detected through the NRSSS since 1996. During the 2009-2016 RVA seasons, 2134 RVA-positive fecal specimens were sent to the CDC for analysis of the VP7 and VP4 genes by RT-PCR genotyping assays and sequencing. During 2009-2011, RVA genotype G3P[8] dominated, while G12P[8] was the dominant genotype during 2012-2016. Vaccine strains were detected in 1.7% of specimens and uncommon/unusual strains, including equine-like G3P[8] strains, were found in 1.9%. Phylogenetic analyses showed limited VP7 and VP4 sequence variation within the common genotypes with 1-3 alleles/lineages identified per genotype. A review of 20 years of NRSSS surveillance showed two changes in genotype dominance, from G1P[8] to G3P[8] and then G3P[8] to G12P[8]. A better understanding of the long-term effects of vaccine use on epidemiological and evolutionary dynamics of circulating RVA strains requires continued surveillance.

Evolutionary changes between pre- and post-vaccine South African group A G2P[4] rotavirus strains, 2003-2017

The transient upsurge of G2P[4] group A rotavirus (RVA) after Rotarix vaccine introduction in several countries has been a matter of concern. To gain insight into the diversity and evolution of G2P[4] strains in South Africa pre- and post-RVA vaccination introduction, whole-genome sequencing was performed for RVA positive faecal specimens collected between 2003 and 2017 and samples previously sequenced were obtained from GenBank (n=103; 56 pre- and 47 post-vaccine). Pre-vaccine G2 sequences predominantly clustered within sub-lineage IVa-1. In contrast, post-vaccine G2 sequences clustered mainly within sub-lineage IVa-3, whereby a radical amino acid (AA) substitution, S15F, was observed between the two sub-lineages. Pre-vaccine P[4] sequences predominantly segregated within sub-lineage IVa while post-vaccine sequences clustered mostly within sub-lineage IVb, with a radical AA substitution R162G. Both S15F and R162G occurred outside recognised antigenic sites. The AA residue at position 15 is found within the signal sequence domain of Viral Protein 7 (VP7) involved in translocation of VP7 into endoplasmic reticulum during infection process. The 162 AA residue lies within the hemagglutination domain of Viral Protein 4 (VP4) engaged in interaction with sialic acid-containing structure during attachment to the target cell. Free energy change analysis on VP7 indicated accumulation of stable point mutations in both antigenic and non-antigenic regions. The segregation of South African G2P[4] strains into pre- and post-vaccination sub-lineages is likely due to erstwhile hypothesized stepwise lineage/sub-lineage evolution of G2P[4] strains rather than RVA vaccine introduction. Our findings reinforce the need for continuous whole-genome RVA surveillance and investigation of contribution of AA substitutions in understanding the dynamic G2P[4] epidemiology.

[Molecular monitoring of the rotavirus (Reoviridae: Sedoreovirinae: Rotavirus: Rotavirus A) strains circulating in Nizhny Novgorod (2012-2020): detection of the strains with the new genetic features]

INTRODUCTION

The pentavalent rotavirus vaccine has been registered in Russia, however, the vaccination coverage remains low, and an annual increase in the incidence of rotavirus infection is unavoidable. In this regard, molecular monitoring of rotaviruses in order to search for new variants possessing epidemic potential is an urgent task.

MATERIAL AND METHODS

PCR genotyping and VP4 and VP7 genes sequencing were used to characterize rotaviruses circulating in Nizhny Novgorod in 2012-2020. The phylogenetic analysis of the strains was carried out using the BEAST software package.

RESULTS

The spectrum included 17 genotypes with predominance of G9P[8] (37,4%). Detected in this study genotypes G1P[4], G1P[9], G2P[8], G4P[4], G4P[6], G8P[8], and G9P[4] were not previously identified in Nizhny Novgorod. The circulation of DS-1-like strains possessing genotypes G1P[8], G3P[8], G8P[8], or G9P[8] and a short RNA pattern had been shown. Rotaviruses of the common genotypes were genetically heterogeneous and belonged to different phylogenetic lineages and/or sublineages (P[4]-IV-a; P[4]-IV-b; P[8]-3.1; P[8]-3.3; P[8]-3.4 and P[8]-3.6; G1-I; G1-II; G2-IVa-1; G2-IVa-3; G3-1; G3-3; G4-I-c; G9-III; G9-VI).

DISCUSSION

These results extend the available data on the genotypic structure of rotavirus populations in Russia and show the genetic diversity of viral strains. G3P[8] DS-1-like viruses were representatives of the G3-1 lineage, new for the territory of Russia, and had the largest number of amino acid substitutions in the VP7 antigenic epitopes.

CONCLUSION

The emergence and spread of strains with new genetic features may allow rotavirus to overcome the immunological pressure formed by natural and vaccine-induced immunity, and maintain or increase the incidence of rotavirus infection.

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.

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.

Increasing detection of rotavirus G2P[4] strains in Nizhny Novgorod, Russia, between 2016 and 2019

Rotavirus infection is one of the leading causes of acute gastroenteritis in children in their first years of life. We studied the genotypic diversity of rotavirus A (RVA) strains in Nizhny Novgorod, Russia, during the period 2016-19. In total, 4714 samples of faeces from children admitted to the Nizhny Novgorod Hospital for Infectious Diseases with acute gastroenteritis were examined. The share of rotavirus-positive samples was 31.5% in 2016-17. It decreased to 21.6% in 2018-19. In Nizhny Novgorod, all six global types of RVA were detected (G1P[8], G2P[4], G3P[8], G4P[8], G9P[8] and G12P[8]), as well as sporadic samples with genotypes G9P[4], G3P[9], G9P[9], G8P[8], G2P[8], G4P[4], G3P[9]. The fraction of strains with genotype G2P[4] gradually increased from 5.9% in 2016-17 to 39.1% in 2018-19. Simultaneously, the proportion of G9P[8] strains decreased from 63.2% to 27.7% in the same period. Phylogenetic analysis showed that rotaviruses with the G2P[4] genotype carried ubiquitous alleles of the VP7 and VP4 genes during the period of their prevalence: G2-IVa-1 and G2-IVa-3; P[4]-IVa and P[4]-IVb. As rotavirus vaccination is not widely used in the region because it is not included in the national vaccination calendar in Russia so far, the increase in the number of G2P[4] RVA is likely due to natural strain fluctuations.

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.

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.

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.

Evolution of P[8], P[4], and P[6] VP8* genes of human rotaviruses globally reported during 1974 and 2017: possible implications for rotavirus vaccines in development

Non-replicating parenteral rotavirus (RV) vaccine candidates are in development in an attempt to overcome the lower efficacy and effectiveness of oral RV vaccines in low-income countries. One of the leading candidates is a truncated recombinant VP8* protein, expressed in Escherichia coli from original sequences of the prototype RV genotypes P[8], P[4], or P[6] isolated before 1983. Since VP8* is highly variable, it was considered useful to examine the evolutionary changes of RV strains reported worldwide over time in relation to the three P2-VP8 vaccine strains. Here, we retrieved from the GenBank 6,366 RV VP8* gene sequences of P[8], P[4], or P[6] strains isolated between 1974 and 2017, in 77 countries, and compared them with those of the three P2-VP8 vaccine strains: Wa (USA, 1974, G1P[8]), DS-1 (USA, 1976, G2P[4]), and 1076 (Sweden, 1983, G2P[6]). Phylogenetic analysis showed that 94.9% (4,328/4,560), 99.8% (1,141/1,143), and 100% (663/663) of the P[8], P[4], and P[6] strains, respectively, reported globally between 1974 and 2018 belong to non-vaccine lineages. These P[8], P[4], and P[6] RV strains have a mean of 9%, 5%, and 6% amino acid difference from the corresponding vaccine strains. Additionally, in the USA, the mean percentage difference between all the P[8] RV strains and the original Wa strain increased over time: 4% (during 1974-1980), 5% (1988-1991), and 9% (2005-2013). Our analysis substantiated high evolutionary changes in VP8* of the P[8], P[4], and P[6] major RV strains and their increasing variations from the candidate subunit vaccine strains over time. These findings may have implications for the development of new RV vaccines.

Whole genome analyses of DS-1-like Rotavirus A strains detected in children with acute diarrhoea in southern Mozambique suggest several reassortment events

We report the first whole genome constellations of Mozambican rotavirus A strains detected between 2012 and 2013 in the Mavalane General Hospital in Maputo city and Manhiça District Hospital in the Manhiça district. Consensus sequences for ten DS-1-like strains (G2P[4] and G8P[4]) were identified with an Illumina Miseq platform using cDNA prepared from dsRNA extracted from stool samples, without genome amplification or prior adaptation to cell culture. Comparison of previously reported genotyping results and the consensus sequences described in this study, indicated that the genotype primers specific for G12 and P[4] might require revision. Phylogenetic analyses indicated diversity among the G2P[4] Mozambican strains and suggested reassortment between G2P[4] and G8P[4] Mozambican strains, as well as the intragenogroup reassortment of all the genome segments encoding VP1, 2, 3 and 6 for strain RVA/Human-wt/MOZ/0045/2012G8P[4]. These results highlight the necessity to determine whole genome constellations to confirm surveillance data in Africa and to monitor the growing diversity in DS-1-like strains.

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

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

GENETIC RELATEDNESS OF EPIZOOTIC HEMORRHAGIC DISEASE VIRUS SEROTYPE 2 FROM 2012 OUTBREAK IN THE USA

During summer and early fall of 2012, the US experienced the largest outbreak of hemorrhagic disease (HD) on record; deer (both Odocoileus virginianus and Odocoileus hemionus) in 35 states were affected, including many northern states where HD typically does not occur. Epizootic hemorrhagic disease virus (EHDV) was the predominant virus isolated, with serotype 2 (EHDV-2) representing 66% (135/205) of all isolated viruses. Viruses within the EHDV serogroup are genetically similar, but we hypothesized that subtle genetic distinctions between viruses would exist across the geographic range of the outbreak if multiple EHDV-2 strains were responsible. We examined viral relatedness and molecular epidemiology of the outbreak by sequencing the mammalian binding protein (VP2) gene and the insect vector binding protein (VP7) gene of 34 EHDV-2 isolates from 2012 across 21 states. Nucleotide sequences of VP2 had 99.0% pairwise identity; VP7 nucleotide sequences had 99.1% pairwise identity. Very few changes were observed in either protein at the amino acid level. Despite the high genetic similarity between isolates, subtle nucleotide differences existed. Both VP2 and VP7 gene sequences separated into two distinct clades based on patterns of single-nucleotide polymorphisms after phylogenetic analysis. The clades were divided geographically into eastern and western clades, although those divisions were not identical between VP2 and VP7. There was also an association between percent sequence identity and geographic distance between isolates. We concluded that multiple EHDV-2 strains contributed to this outbreak.

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

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

Emergence of Human G2P[4] Rotaviruses in the Post-vaccination Era in South Korea: Footprints of Multiple Interspecies Re-assortment Events

After the introduction of two global rotavirus vaccines, RotaTeq in 2007 and Rotarix in 2008 in South Korea, G1[P8] rotavirus was the major rotavirus genotype in the country until 2012. However, in this study, an emergence of G2P[4] as the dominant genotype during the 2013 to 2015 season has been reported. Genetic analysis revealed that these viruses had typical DS-1-like genotype constellation and showed evidence of re-assortment in one or more genome segments, including the incorporation of NSP4 genes from strains B-47/2008 from a cow and R4/Haryana/2007 from a buffalo in India, and the VP1 and VP3 genes from strain GO34/1999 from a goat in Bangladesh. Compared to the G2 RotaTeq vaccine strain, 17–24 amino acid changes, specifically A87T, D96N, S213D, and S242N substitutions in G2 epitopes, were observed. These results suggest that multiple interspecies re-assortment events might have contributed to the emergence of G2P[4] rotaviruses in the post-vaccination era in South Korea.

Genetic analysis of rotavirus G2P[4] strains in Pune, Western India: circulation of a novel reassortant bearing E6 NSP4 genotype

In India, G2P[4] strains are known to be the second most predominant group A rotaviruses causing acute gastroenteritis among children. This study was performed to determine the diversity within VP7(G), VP4(P), VP6(I) and NSP4(E) genes of 16 G2P[4] rotavirus strains detected in children hospitalized for acute gastroenteritis in Pune, Western India during 2009-2013. Fourteen strains showed G2-P[4]-I2-E2 and two strains showed G2-P[4]-I2-E6 genotype constellation. Phylogenetic analysis showed their clustering into G2-IV-a3, P[4]-5bi/ii, I2-3ii and E2-4i/ii or E6 genotypes/lineages. These data reveal inter- and/or intra-genotypic variations in a genogroup-2 constellation of G2P[4] rotavirus strains circulating in Pune, Western India, providing evidence of a novel G2P[4] reassortant bearing a rare NSP4 genotype, E6 during 2009-2013.

Genetic determinants restricting the reassortment of heterologous NSP2 genes into the simian rotavirus SA11 genome

Rotaviruses (RVs) can evolve through the process of reassortment, whereby the 11 double-stranded RNA genome segments are exchanged among strains during co-infection. However, reassortment is limited in cases where the genes or encoded proteins of co-infecting strains are functionally incompatible. In this study, we employed a helper virus-based reverse genetics system to identify NSP2 gene regions that correlate with restricted reassortment into simian RV strain SA11. We show that SA11 reassortants with NSP2 genes from human RV strains Wa or DS-1 were efficiently rescued and exhibit no detectable replication defects. However, we could not rescue an SA11 reassortant with a human RV strain AU-1 NSP2 gene, which differs from that of SA11 by 186 nucleotides (36 amino acids). To map restriction determinants, we engineered viruses to contain chimeric NSP2 genes in which specific regions of AU-1 sequence were substituted with SA11 sequence. We show that a region spanning AU-1 NSP2 gene nucleotides 784–820 is critical for the observed restriction; yet additional determinants reside in other gene regions. In silico and in vitro analyses were used to predict how the 784–820 region may impact NSP2 gene/protein function, thereby informing an understanding of the reassortment restriction mechanism.

Emergence of human G2P[4] rotaviruses containing animal derived gene segments in the post-vaccine era

The introduction of Rotarix into the Belgian immunization program in June 2006 coincided with an increase of the relative prevalence of G2P[4] strains. However, the genetic composition of these persistent G2P[4] strains has not been investigated. Therefore, we have investigated the NSP4 gene of 89 Belgian G2P[4] strains detected between 1999 and 2013, covering both pre- and post-vaccination periods. The NSP4 genes were divided over seven separate clusters of which six were more closely related to animal than to human strains. The NSP4 genes that clustered more closely to animal DS-1-like strains were isolated after 2004–2005 and were found throughout multiple seasons. Complete genome sequencing of 28 strains identified several other gene segments that clustered more closely to animal than to human DS-1-like strains. These findings suggest that frequent interspecies reassortments may have played a role in the spread of G2P[4] rotaviruses in the post-vaccination period in Belgium.

Rotavirus vaccine and health-care utilization for rotavirus gastroenteritis in Tsu City, Japan

Background

Rotavirus vaccines were introduced in Japan in November 2011. We evaluated the subsequent reduction of the health-care burden of rotavirus gastroenteritis.

Methods

We conducted active surveillance for rotavirus gastroenteritis among children under 5 years old before and after the vaccine introduction. We surveyed hospitalization rates for rotavirus gastroenteritis in children in Tsu City, Mie Prefecture, Japan, from 2007 to 2015 and surveyed the number of outpatient visits at a Tsu City clinic from 2010 to 2015. Stool samples were obtained for rotavirus testing and genotype investigation. We assessed rotavirus vaccine coverage for infants living in Tsu City.

Results

In the pre-vaccine years (2007–2011), hospitalization rates for rotavirus gastroenteritis in children under 5 years old were 5.5, 4.3, 3.1 and 3.9 cases per 1000 person-years, respectively. In the post-vaccine years (2011–2015), the rates were 3.0, 3.5, 0.8 and 0.6 cases per 1000 person-years, respectively. The hospitalization rate decreased significantly in the 2013–2014 and 2014–2015 seasons compared to the average of the seasons before vaccine introduction (P < 0.0001). In one pre-vaccine year (2010–2011), the number of outpatient visits due to the rotavirus infection was 66. In the post-vaccine years (2011–2015), the numbers for each season was 23, 23, 7 and 5, respectively. The most dominant rotavirus genotype shifted from G3P[8] to G1P[8] and to G2P[4]. The coverage of one dose of rotavirus vaccine in Tsu City was 56.5% in 2014.

Conclusion

After the vaccine introduction, the hospitalization rates and outpatient visits for rotavirus gastroenteritis greatly decreased.

Molecular characterization of emerging G9P[4] rotavirus strains possessing a rare E6 NSP4 or T1 NSP3 genotype on a genogroup-2 backbone using a refined classification framework

Rotavirus infections associated with unusual strains are an emerging concern in rotavirus vaccination programmes. Recently, an increase in circulation of unusual G9P[4] strains was reported from different regions of India, placing this genotype in third position, after G1P[8] and G2P[4], of the most common rotavirus strains. The aim of the present study was to analyse the complete genomic constellation of three G9P[4] strains (RV09, RV10 and RV11), determine their genetic relatedness to other genogroup-2 strains and understand the evolution of a rare E6 and other NSP4 genotypes. All strains revealed the presence of a genogroup-2 backbone, with RV09 constituting the NSP3 T1 genotype and RV10 and RV11 bearing the NSP4 E6 genotype. A refined criterion adopted to classify the nine internal gene segments of G2P[4] and non-G2P[4] strains with the genogroup-2 backbone into lineages and sub-lineages indicated divergence of >8 % (except NSP1: >5.5 %) for lineages and >3 % for sub-lineages. The VP1 and/or VP3 genes of study strains showed close relationships with animal-like human rotaviruses. The estimated evolutionary rate for the NSP4 E6 genotype was marginally higher (3.78×10-3 substitutions per site per year) than that of genotypes E1 (2.6×10-3 substitutions per site per year) and E2 (3.06×10-3 substitutions per site per year), suggesting a step towards adaptation of E6 on a genogroup-2 backbone. The time and origin of the most recent common ancestor of E6 genotype were estimated to be 1981 and South Asia, respectively. Full-genome and evolutionary analyses performed in this study for G9P[4] strains will help better understand the extent of gene reassortment and origin in unusual rotavirus strains that may remain viable and cause infections in humans.

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

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

Genomic constellation and evolution of Ghanaian G2P[4] rotavirus strains from a global perspective

Understanding of the genetic diversity and evolution of Rotavirus A (RVA) strains, a common cause of severe diarrhoea in children, needs to be based on the analysis at the whole genome level in the vaccine era. This study sequenced the whole genomes of six representative G2P[4] strains detected in Ghana from 2008 to 2013, and analysed them phylogenetically with a global collection of G2P[4] strains and African non-G2P[4] DS-1-like strains. The genotype constellation of the study strains was G2-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Strains from the same season were highly identical across the whole genome while strains from different seasons were more divergent from each other. The VP7, VP4, VP2, NSP1, and NSP5 genes belonged to lineage IVa; the VP6, VP1, NSP2, and NSP3 genes belonged to lineage V, and all these genes evolved in the same fashion as the global strains. In the NSP4 gene, lineages V (2008) and X (2009) were replaced by VI (2012/2013) whereas in the VP3 gene, lineage V (2008/2009) was replaced by VII (2012/2013) and these replacements coincided with the vaccine introduction period (2012). The evolutionary rate of the NSP4 gene was 1.2×10^-3 substitutions/site/year and was rather comparable to that of the remaining 10 genes. The multiple NSP4 lineages were explained by intra-genotype reassortment with co-circulating African human DS-1-like strains bearing G2[6], G3P[6], G6[6] and G8. There was no explicit evidence of the contribution of animal RVA strains to the genome of the Ghanaian G2P[4] strains. In summary, this study revealed the dynamic evolution of the G2P[4] strains through intra-genotype reassortment events leading to African specific lineages such IX and X in the NSP4 gene. So far, there was no evidence of a recent direct involvement of animal RVA genes in the genome diversity of African G2P[4] strains.

Reassortment in segmented RNA viruses: mechanisms and outcomes

Segmented RNA viruses are widespread in nature and include important human, animal and plant pathogens, such as influenza viruses and rotaviruses. Although the origin of RNA virus genome segmentation remains elusive, a major consequence of this genome structure is the capacity for reassortment to occur during co-infection, whereby segments are exchanged among different viral strains. Therefore, reassortment can create viral progeny that contain genes that are derived from more than one parent, potentially conferring important fitness advantages or disadvantages to the progeny virus. However, for segmented RNA viruses that package their multiple genome segments into a single virion particle, reassortment also requires genetic compatibility between parental strains, which occurs in the form of conserved packaging signals, and the maintenance of RNA and protein interactions. In this Review, we discuss recent studies that examined the mechanisms and outcomes of reassortment for three well-studied viral families - Cystoviridae, Orthomyxoviridae and Reoviridae - and discuss how these findings provide new perspectives on the replication and evolution of segmented RNA viruses.

Whole genome analysis of Vietnamese G2P[4] rotavirus strains possessing the NSP2 gene sharing an ancestral sequence with Chinese sheep and goat rotavirus strains

Because imminent introduction into Vietnam of a vaccine against Rotavirus A is anticipated, baseline information on the whole genome of representative strains is needed to understand changes in circulating strains that may occur after vaccine introduction. In this study, the whole genomes of two G2P[4] strains detected in Nha Trang, Vietnam in 2008 were sequenced, this being the last period during which virtually no rotavirus vaccine was used in this country. The two strains were found to be >99.9% identical in sequence and had a typical DS-1 like G2-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2 genotype constellation. Analysis of the Vietnamese strains with >184 G2P[4] strains retrieved from GenBank/EMBL/DDBJ DNA databases placed the Vietnamese strains in one of the lineages commonly found among contemporary strains, with the exception of the NSP2 and NSP4 genes. The NSP2 genes were found to belong to a previously undescribed lineage that diverged from Chinese sheep and goat rotavirus strains, including a Chinese rotavirus vaccine strain LLR with 95% nucleotide identity; the time of their most recent common ancestor was 1975. The NSP4 genes were found to belong, together with Thai and USA strains, to an emergent lineage (VIII), adding further diversity to ever diversifying NSP4 lineages. Thus, there is a need to enhance surveillance of locally-circulating strains from both children and animals at the whole genome level to address the effect of rotavirus vaccines on changing strain distribution.

Diversity of group A rotavirus genes detected in the Triângulo Mineiro region, Minas Gerais, Brazil

Group A rotaviruses are the main causative agent of infantile gastroenteritis. The segmented nature of the viral genome allows reassortment of genome segments, which can generate genetic variants. In this study, we characterized the diversity of the VP7, VP4 (VP8*), VP6, NSP4, and NSP5 genes of the rotaviruses that circulated from 2005 to 2011 in the Triângulo Mineiro (TM) region of Brazil. Samples with genotypes G2 (sublineages IVa-1 and IVa-3), G1 (sublineage I-A), G9 (lineage III), G12 (lineages II and III), G8 (lineage II), G3 (lineage III), P[4] (sublineages IVa and IVb), P[8] (sublineages P[8]-3.6, P[8]-3.3, and P[8]-3.1), I2 (lineage VII), E2 (lineages VI, XII, and X), and H2 (lineage III) were identified. The associations found in the samples were G1, G9, or G12 with P[8]-I1-E1-H1; G2 or G8 with P[4]-I2-E2-H2; G12 with I3-E3-H6; and G3 with P[4]-I2-E3-H3 (previously unreported combination). Reassortment events in G2P[4] strains and an apparent pattern of temporal segregation within the lineages were observed. Five TM samples contained genes that exhibited high nucleotide and amino acid identities with strains of animal origin. The present study includes a period of pre- and post-introduction of rotavirus vaccination in all Brazilian territories, thereby serving as a basis for monitoring changes in the genetic constitution of rotaviruses. The results also contribute to the understanding of the diversity and evolution of rotaviruses in a global context.

G2 rotavirus within an emergent VP7 evolutionary lineage circulating in children with acute diarrhea in Guangxi Province of China, 2014

Routine surveillance revealed that the prevalence of P[4] rotaviruses circulating in children with acute diarrhea in Guangxi Province, China, increased in 2014. However, VP7 genotyping for these P[4] rotaviruses was unsuccessful. Exhaustive database searching and sequence analysis indicated that the G genotype of these P[4] rotaviruses was G2, and the VP7 genes clustered with recently emerging G2 strains in several countries within an emergent evolutionary lineage that was distinct from the previously designated lineages I-IV as well as lineage V including porcine rotaviruses. Further studies are essential to monitor the potential global spread of this emerging G2 rotavirus.

Phylogenetic investigation of enteric bovine coronavirus in Ireland reveals partitioning between European and global strains

Background

Bovine coronavirus is a primary cause of neonatal calf diarrhea worldwide, and is also associated with acute diarrhea in adult cattle during the winter season. There are no reports on molecular characterization of bovine coronavirus in Ireland, and little data exists apart from serological studies.

Findings

In this study, 11 neonatal (mean age 9 days) calf BCoV strains from the south of Ireland were collected over a one year period and characterized using molecular methods. The spike gene which encodes a protein involved in viral entry, infectivity and immune response shows the most variability amongst the isolates and was subsequently selected for in depth analysis. Phylogenetic analysis of the spike gene revealed that the Irish strains clustered with novel BCoV strains from Europe in a unique clade, possibly indicating lineage partitioning. Direct analysis of alignments identified amino acid changes in the spike protein unique to the Irish clade.

Conclusion

Thus, monitoring of bovine coronavirus in Ireland is important as the current isolates in circulation in the south of Ireland may be diverging from the available vaccine strain, which may have implications regarding future BCoV vaccine efficacy.

Genome-Wide Evolutionary Analyses of G1P[8] Strains Isolated Before and After Rotavirus Vaccine Introduction

Rotaviruses are the most important etiological agent of acute gastroenteritis in young children worldwide. Among the first countries to introduce rotavirus vaccines into their national immunization programs were Belgium (November 2006) and Australia (July 2007). Surveillance programs in Belgium (since 1999) and Australia (since 1989) offer the opportunity to perform a detailed comparison of rotavirus strains circulating pre- and postvaccine introduction. G1P[8] rotaviruses are the most prominent genotype in humans, and a total of 157 G1P[8] rotaviruses isolated between 1999 and 2011 were selected from Belgium and Australia and their complete genomes were sequenced. Phylogenetic analysis showed evidence of frequent reassortment among Belgian and Australian G1P[8] rotaviruses. Although many different phylogenetic subclusters were present before and after vaccine introduction, some unique clusters were only identified after vaccine introduction, which could be due to natural fluctuation or the first signs of vaccine-driven evolution. The times to the most recent common ancestors for the Belgian and Australian G1P[8] rotaviruses ranged from 1846 to 1955 depending on the gene segment, with VP7 and NSP4 resulting in the most recent estimates. We found no evidence that rotavirus population size was affected after vaccine introduction and only six amino acid sites in VP2, VP3, VP7, and NSP1 were identified to be under positive selective pressure. Continued surveillance of G1P[8] strains is needed to determine long-term effects of vaccine introductions, particularly now rotavirus vaccines are implemented in the national immunization programs of an increasing number of countries worldwide.

Rotavirus genotypes and the indigenous children of Brazilian midwest in the vaccine era, 2008-2012: Footprints of animal genome

World group A rotavirus (RVA) surveillance data provides useful estimates of the disease burden, however, indigenous population might require special consideration. The aim of this study was to describe the results of G- and P-types from Brazilian native children ≤ 3 years. Furthermore, selected strains have been analyzed for the VP7, VP6, VP4, and NSP4 encoding genes in order to gain insight into genetic variability of Brazilian strains. A total of 149 samples, collected during 2008-2012, were tested for RVA using ELISA and PAGE, following by RT-PCR and sequencing. RVA infection was detected in 8.7% of samples (13/149). Genotype G2P[4] was detected in 2008 and 2010, G8P[6] in 2009, and G3P[8] in 2011. The phylogenetic analysis of the VP7 and VP4 genes grouped the Brazilian G2P[4] and G3P[8] strains within the lineages currently circulating in humans worldwide. However, the phylogenetic analysis of the VP6 and NSP4 from the Brazilian G2P[4] strains, and the VP7 and NSP4 from the Brazilian G3P[8] strains suggest a distant common ancestor with different animal strains (bovine, caprine, and porcine). The epidemiological and genetic information obtained in the present study is expected to provide an updated understanding of RVA genotypes circulating in the native infant population, and to formulate policies for the use of RVA vaccines in indigenous Brazilian people. Moreover, these results highlight the great diversity of human RVA strains circulating in Brazil, and an in-depth surveillance of human and animal RVA will lead to a better understanding of the complex dynamics of RVA evolution.

NCBI viral genomes resource

Recent technological innovations have ignited an explosion in virus genome sequencing that promises to fundamentally alter our understanding of viral biology and profoundly impact public health policy. Yet, any potential benefits from the billowing cloud of next generation sequence data hinge upon well implemented reference resources that facilitate the identification of sequences, aid in the assembly of sequence reads and provide reference annotation sources. The NCBI Viral Genomes Resource is a reference resource designed to bring order to this sequence shockwave and improve usability of viral sequence data. The resource can be accessed at http://www.ncbi.nlm.nih.gov/genome/viruses/ and catalogs all publicly available virus genome sequences and curates reference genome sequences. As the number of genome sequences has grown, so too have the difficulties in annotating and maintaining reference sequences. The rapid expansion of the viral sequence universe has forced a recalibration of the data model to better provide extant sequence representation and enhanced reference sequence products to serve the needs of the various viral communities. This, in turn, has placed increased emphasis on leveraging the knowledge of individual scientific communities to identify important viral sequences and develop well annotated reference virus genome sets.

Analysis of rotavirus genotypes in Korea during 2013: an increase in the G2P[4] genotype after the introduction of rotavirus vaccines

BACKGROUND

Group A rotavirus is the leading cause of acute gastroenteritis in children worldwide. We investigated G and P genotypes of group A rotavirus strains isolated from patients during 2013 and investigated which genotypes were identified from vaccinated patients.

METHODS

From January to December 2013, 2235 fecal specimens were tested for rotavirus antigen, of which 374 specimens (16.7%) showed positive results. Strains from 288 rotavirus-positive specimens were genotyped using PCR and sequencing, and individual patients' corresponding vaccine histories were investigated through the Korean Center for Disease Control website.

RESULTS

G2 (22.6%) and P[4] (24.0%) were the most frequently identified G and P genotypes, respectively; accordingly, G2P[4] (19.8%) was the most prevalent G/P genotype observed in this period. G4P[6] (10.1%) was the second most prevalent G/P genotype and was mostly detected in neonates. Other genotypes, G1P[8], G9P[8], G1P[6], and G3P[6], were also detected. Of 288 rotavirus-positive specimens, 48 specimens were obtained from previously vaccinated patients. G2P[4] was also the genotype most frequently isolated from vaccinated patients. VP7 epitope analysis of G1P[8] and G2P[4] strains showed at least one amino acid differences in comparison with Rotarix and RotaTeq vaccine strains. The genotypic distribution of rotavirus strains in Korea has been shown temporal and geographical differences.

CONCLUSION

This study showed that G2P[4] was the genotype most frequently isolated from both vaccinated and unvaccinated patients in Korea during 2013. However, it is unclear whether the change of predominant genotype is due to the effect of vaccination or due to natural variation.

Review of global rotavirus strain prevalence data from six years post vaccine licensure surveillance: is there evidence of strain selection from vaccine pressure?

Comprehensive reviews of pre licensure rotavirus strain prevalence data indicated the global importance of six rotavirus genotypes, G1P[8], G2P[4], G3P[8], G4P[8], G9P[8] and G12P[8]. Since 2006, two vaccines, the monovalent Rotarix (RV1) and the pentavalent RotaTeq (RV5) have been available in over 100 countries worldwide. Of these, 60 countries have already introduced either RV1 or RV5 in their national immunization programs. Post licensure vaccine effectiveness is closely monitored worldwide. This review aimed at describing the global changes in rotavirus strain prevalence over time. The genotype distribution of the nearly 47,000 strains that were characterized during 2007-2012 showed similar picture to that seen in the preceding period. An intriguing finding was the transient predominance of heterotypic strains, mainly in countries using RV1. Unusual and novel antigen combinations continue to emerge, including some causing local outbreaks, even in vaccinated populations. In addition, vaccine strains have been found in both vaccinated infants and their contacts and there is evidence for genetic interaction between vaccine and wild-type strains. In conclusion, the post-vaccine introduction strain prevalence data do not show any consistent pattern indicative of selection pressure resulting from vaccine use, although the increased detection rate of heterotypic G2P[4] strains in some countries following RV1 vaccination is unusual and this issue requires further monitoring.

Characterization of G2P[4] rotavirus strains causing outbreaks of gastroenteritis in the Northern Territory, Australia, in 1999, 2004 and 2009

Outbreaks of rotavirus diarrhea cause a large disease burden in the Alice Springs region of the Northern Territory, Australia. The introduction of the rotavirus vaccine Rotarix® has been associated with an increase in detection of G2P[4] strains in many countries. However, G2P[4] emergence has also been observed in vaccine-naive countries, suggesting a general global increase in the circulation of G2P[4] strains. A G2P[4] rotavirus outbreak occurred in 2009, 28 months after the introduction of the Rotarix® vaccine and 43 children were hospitalized. Pre-vaccine introduction, G2P[4] strains were observed associated with large outbreaks in 1999 and 2004. To determine the genetic relationship between these strains whole genome sequence analysis was conducted on representative strains from each of the G2P[4] outbreaks, in 1999, 2004 and 2009. Phylogenetic analysis revealed the majority of genes from 2009 outbreak strain clustered with contemporary global strains, while the VP7 gene clustered with contemporary and older strains and was antigenically distinct to the majority of contemporary global G2P[4] strains; suggesting the strain was an intragenogroup reassortant. The 1999 and 2009 strains appear to share similar evolutionary origins, and both had a high degree of genetic identity to previously identified Australian and global strains. Conversely, the 2004 outbreak strain was more divergent in comparison to Australian and global strains. The 1999 and 2004 outbreaks likely occurred due to the accumulation of immunologically naïve children in the population following low levels of G2P[4] rotavirus disease in the community in the years prior to each outbreak. The 2009 outbreak was associated with moderate vaccine coverage in the population and vaccine efficacy against the strain was low. The circulation of this unusual strain in the population combined with low vaccine coverage and diminished vaccine efficacy likely contributed to the outbreak occurring in this population.

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

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

Longitudinal analysis of VP7 gene of group A human rotavirus G2P[4] strains circulating in the pre-vaccine era in Sapporo, Japan from 1991 to 2011

Sequence analysis of the VP7 gene in 23 group A human rotavirus G2P[4] strains obtained during 1991-2011, that is, the pre-vaccine era, in Sapporo, Japan showed considerable genetic diversity, mainly in variable regions. Recent G2P[4] epidemic strains were located in sublineage IVa with a distinctive substitution of D96N. This study provides background data on the genetic variability of G2P[4] rotavirus-VP7 gene prior to the widespread use of rotavirus vaccines in Japan.

Absence of genetic differences among G10P[11] rotaviruses associated with asymptomatic and symptomatic neonatal infections in Vellore, India

Rotaviruses (RVs) are leading causes of severe diarrhea and vomiting in infants and young children. RVs with G10P[11] genotype specificity have been associated with symptomatic and asymptomatic neonatal infections in Vellore, India. To identify possible viral genetic determinants responsible for differences in symptomology, the genome sequences of G10P[11] RVs in stool samples of 19 neonates with symptomatic infections and 20 neonates with asymptomatic infections were determined by Sanger and next-generation sequencing. The data showed that all 39 viruses had identical genotype constellations (G10-P[11]-I2-R2-C2-M2-A1-N1-T1-E2-H3), the same as those of the previously characterized symptomatic N155 Vellore isolate. The data also showed that the RNA and deduced protein sequences of all the Vellore G10P[11] viruses were nearly identical; no nucleotide or amino acid differences were found that correlated with symptomatic versus asymptomatic infection. Next-generation sequencing data revealed that some stool samples, both from neonates with symptomatic infections and from neonates with asymptomatic infections, also contained one or more positive-strand RNA viruses (Aichi virus, astrovirus, or salivirus/klassevirus) suspected of being potential causes of pediatric gastroenteritis. However, none of the positive-strand RNA viruses could be causally associated with the development of symptoms. These results indicate that the diversity of clinical symptoms in Vellore neonates does not result from genetic differences among G10P[11] RVs; instead, other undefined factors appear to influence whether neonates develop gastrointestinal disease symptoms.

IMPORTANCE

Rotavirus (RV) strains have been identified that preferentially replicate in neonates, in some cases, without causing gastrointestinal disease. Surveillance studies have established that G10P[11] RVs are a major cause of neonatal infection in Vellore, India, with half of infected neonates exhibiting symptoms. We used Sanger and next-generation sequencing technologies to contrast G10P[11] RVs recovered from symptomatic and asymptomatic neonates. Remarkably, the data showed that the RNA genomes of the viruses were virtually indistinguishable and lacked any differences that could explain the diversity of clinical outcomes among infected Vellore neonates. The sequencing results also indicated that some symptomatic and some asymptomatic Vellore neonates were infected with other enteric viruses (Aichi virus, astrovirus, salvirus/klassevirus); however, none could be correlated with the presence of symptoms in neonates. Together, our findings suggest that other poorly defined factors, not connected to the genetic makeup of the Vellore G10P[11] viruses, influence whether neonates develop gastrointestinal disease symptoms.