- Open Access
Human rotavirus strains circulating in Venezuela after vaccine introduction: predominance of G2P and reemergence of G1P
Virology Journalvolume 14, Article number: 58 (2017)
Rotavirus (RV) is the most common cause of severe childhood diarrhea worldwide. Despite Venezuela was among the first developing countries to introduce RV vaccines into their national immunization schedules, RV is still contributing to the burden of diarrhea. Concerns exist about the selective pressure that RV vaccines could exert on the predominant types and/or emergence of new strains.
To assess the impact of RV vaccines on the genotype distribution 1 year after the vaccination was implemented, a total of 912 fecal specimens, collected from children with acute gastroenteritis in Caracas from February 2007 to April 2008, were screened, of which 169 (18.5%) were confirmed to be RV positive by PAGE. Rotavirus-associated diarrhea occurred all year-round, although prevailed during the coolest and driest months among unvaccinated children under 24 months old. Of 165 RV strains genotyped for G (VP7) and P (VP4) by seminested multiplex RT-PCR, 77 (46.7%) were G2P and 63 (38.2%) G1P. G9P, G3P and G2P were found in a lower proportion (7.3%). Remarkable was also the detection of <5% of uncommon combinations (G8P, G8P, G1P and G4P) and 3.6% of mixed infections. A changing pattern of G/P-type distribution was observed during the season studied, with complete predominance of G2P from February to June 2007 followed by its gradual decline and the reemergence of G1P, predominant since January 2008. Phylogenetic analysis of VP7 and VP4 genes revealed a high similarity among G2P and global strains belonging to G2-II and P-V lineages. The amino acid substitution 96D → N, related with reemergence of the G2 genotype elsewhere, was observed. The G1P strains from Caracas were grouped into the lineages G1-I and P-III, along with geographically remote G1P rotaviruses, but they were rather distant from Rotarix® vaccine and pre-vaccine strains. Unique amino acid substitutions observed on neutralization domains of the VP7 sequence from Venezuelan post-vaccine G1P could have conditioned their re-emergence and a more efficient dissemination into susceptible population.
The results suggest that natural fluctuations of genotypes in combination with forces driving the genetic evolution could determine the spread of novel strains, whose long-term effect on the efficacy of available vaccines should be determined.
Human rotavirus (RV) is the most important cause of severe, acute gastroenteritis in infants and young children in both developed and developing countries. Until 2008, rotavirus infections were estimated to cause approximately 453,000 deaths per year worldwide in children aged <5 years, most of them living in low-income countries . Two live RV vaccines (oral and attenuated), the monovalent human RV vaccine Rotarix® (GlaxoSmithKline Biologicals, Rixensart, Belgium), containing a single human G1P RV strain, and the pentavalent bovine–human, reassortant vaccine RotaTeq® (Merck & Co., Inc., Whitestation, NJ, United States), were licensed in 2006. Both have demonstrated very good safety and efficacy profiles in large clinical trials mainly in the United States, Europe and Latin America [2–4]. The inclusion of rotavirus vaccine in the immunization programs was recommended by World Health Organization since 2007 in regions where a significant public health impact was anticipated , and later in 2009, the introduction in developing countries was proposed worldwide. The implementation in many countries of the world has demonstrated to reduce the disease burden of rotavirus-specific diarrhea and death related in several regions [6–8].
Venezuela was among the first developing countries to introduce RV vaccines into the national immunization schedules in 2006. Both vaccines have been available: the Rotarix® is the RV vaccine of choice in the public health care sector, and therefore the most used, while the RotaTeq® is obtainable only in private facilities. Although some data have shown decline in hospitalizations and deaths related to severe diarrhea for any causes in children under 5-years of age after RV vaccine introduction , a report from the Center for Disease Control and Prevention (CDC) indicated no substantial changes in the percentage of RV diarrhea cases in 2010 (31%) compared with 2006 (32%) in Venezuela, where coverage has fluctuated around 49% in the last decade . On the other hand, vaccine efficacy has been shown to be remarkably lower in developing countries . Strain surveillance is needed to assess the impact of immunization on the RV diarrhea incidence and the variability of the circulating viruses.
In the last two decades, RV molecular genotyping has provided valuable information about the diversity of RV outer capside (VP7 or “G”, and VP4 or “P”) and some nonstructural proteins from strains circulating throughout the world. Currently, 27 G genotypes and 37 P genotypes have been described [12, 13]. Extensive molecular epidemiological studies globally have indicated that only a small number of genotypes have prevailed: G1P, G2P, G3P and G4P have been the most important in humans worldwide. Since 1995, G9P increased dramatically, being now considered the fifth globally important RV genotype [14–16], and G12P is currently also being increasingly detected around the world [17, 18]. Other G types of animal origin such as G5, G6 and G8 have acquired epidemiological relevance in some geographical areas of Africa, Asia and South America [19–21].
Temporal and geographical changes in genotype prevalence patterns have been seen with periodic emergence of novel strains, particularly in developing countries [14, 22] where the evolutionary dynamics of RV are complicated by a greater diversity. Surveillance of RV infection conducted in the past two decades in Venezuela has showed also a broad diversity and temporal variations for the G and P types circulating, with alternating predominance of G1P, G3P or G4P [23–26] and moderate rates of circulation of G2P strains, in addition to the emergence of G9P and appearance of G8P [25, 26]. The effects of such variability on the efficacy of the vaccines have to be elucidated, and more whether vaccination may lead to the replacement of vaccine-type strains. Reports from Brazil, Australia, Portugal and USA have described changes in RV genotype prevalence following vaccine introduction [27–31], but whether the pattern and distribution of the most prevalent RV strains could be modified by the immunologic pressure exerted by the vaccines remains uncertain. The emergence of RV genotypes that are not efficiently controlled by the immune response induced by the vaccine is a possibility , but additional studies are needed.
The main aim of this study was to estimate the prevalence of RV gastroenteritis and distribution of circulating G (VP7) and P (VP4) genotypes from clinical isolates causing symptomatic infections in children with diarrhea, living in the metropolitan area of Caracas, one year after the RV vaccination started in Venezuela. In addition, phylogenetic analysis based on VP7 and VP4 gene of some isolates collected during 2007-2008 was performed in comparison with pre-vaccination RV strains collected in 2003 and global reference strains. In order to examine the potential modifications under vaccine selective pressure, VP7 and VP4 amino acid sequences were also deduced and analyzed to verify the occurrence of eventual substitutions.
Stool collection, RV testing and controls
Between February 2007 and April 2008, fecal samples were collected from 912 children younger than 10 years of age, who were attended for acute diarrhea, defined as three or more liquid stools over a 24-h period, in a private clinical setting of a medical center in Caracas, Venezuela. Out of 912, 69% proceeded from children over 16 months of age at the starting time of the study, therefore not eligible for rotavirus vaccination.
Stool specimens were immediately screened for the presence of rotavirus by use of an immunochromatographic rapid test, Rota-Strip (Coris BioConcept, Gembloux, Belgium), according to the manufacturer’s instructions. The samples were stored at -20°C until successively tested.
Cell culture-adapted RV strains grown in MA-104 cells were used as control strains for G- and P-genotyping assays. G/P type designations are as recommended by the Rotavirus Classification Working Group  and Reoviridae Study Group of the International Committee on Taxonomy of Viruses. The study was cleared by the IVIC ethic committee. A written informed consent was obtained from the parents or legal guardians at enrollment.
RNA extraction, gel electrophoresis and silver staining
Rotavirus positive samples detected by immunochromatographic test were further analyzed by polyacrylamide gel electrophoresis (PAGE) in a 7% gel following viral RNA genome extraction by phenol and chloroform-isoamyl alcohol treatment, and ethanol precipitation. Double strand RNA (dsRNA) segments of RV were separated by PAGE at room temperature for 3 h at 90 V and stained with silver nitrate as described previously . The electrophoretic migration patterns (electropherotypes) of the RNA segments were also analyzed for further genetic characterization of the strains.
G/P genotyping of rotavirus
Rotavirus G (VP7) and P (VP4) genotypes were determined as described previously . For this purpose, nucleic acids were extracted from 5% fecal suspensions clarified supernatant using the QIAamp® Viral RNA Mini kit (QIAGEN®, Hilden, Germany) according to the manufacturer’s instructions. The extracted RNA was reverse transcribed and G and P genotyping was performed by semi-nested multiplex polymerase chain reaction method (RT-PCR) in two rounds, using type-specific primers for VP7 gene (G1-4, G8-10 types) and for VP4 gene (P, P, P, P and P types) separately, as previously described . PCR products were analysed by agarose gel electrophoresis and ethidium bromide staining.
The partial first-round PCR-derived product of the VP4 and VP7 genes from selected four RV strains circulating in Venezuela during the years 2007/2008 and strains obtained during the year 2003 before vaccine introduction in Valencia (Carabobo state, Venezuela) , representative of the main RV genotypes found, were purified using a commercial column (QIAquick PCR purification kit, QIAGEN®, Hilden, Germany) and automatically sequenced in both directions using BigDye Terminator cycle chemistry and a 3130XL DNA analyzer (Applied Biosystems, Foster City, CA, USA). A commercially obtained lyophilized dose (Lot n. A41CA419A) of Rotarix® vaccine (GlaxoSmithKline Biologicals, Rixensart, Belgium), used in Venezuela during the 2007-2008 season, was reconstituted according to the manufacturer’s instructions, the viral genome was extracted and amplified as described above, and purified VP4/VP7 first-round PCR amplicons also sequenced.
Phylogenetic analysis of VP7 and VP4
Alignment and comparison of nucleotide and deduced amino acid sequences of VP7 and VP4 genes from the selected Venezuelan RV isolates were performed by using BLAST® (Basic Local Alignment Search Tool) web service of U.S. National Library of Medicine (blast.ncbi.nlm.nih.gov), and ClustalW method. Phylogenetic and molecular evolutionary analyses were conducted using MEGA version 6.06 . The phylogenetic trees were generated by neighbor-joining method and with the Kimura’s two-parameter as substitution model, from nucleotide sequences of the RV strains in this study together with global reference strains available in the NCBI (National Center for Biotechnology Information) GenBank database (https://www.ncbi.nlm.nih.gov/nucleotide/) and Rotarix® vaccine strain. The statistical significance of the branches was assessed by bootstrap resampling analysis (1000 replicates).
Protein structure analysis of VP7
A high-quality homology-based model from VP7 deduced amino acid sequences from G1P strains from this study was generated automatically from the ModBase database by the ModWeb server version r181 from the University of California (San Francisco, USA) . The models were built by using the crystal structure (Protein Data Bank identifier [PDB ID], 3FMG) of RV outer capsid protein VP7 trimer in complex with a neutralizing Fab as template . Protein structural analysis and modeling were performed using the UCSF Chimera-Molecular Modeling System v1.11 (University of California, San Francisco) , and the Adobe Illustrator CS6 software was utilized as a design tool of the image.
Nucleotide sequence accession numbers
The VP7 and VP4 gene sequences of the Venezuelan RV isolates reported in this work were deposited in the GenBank nucleotide sequence database under the accession numbers KY039372-KY039373 and KY053848-KY053851 for VP7 gene, and KY056539-KY056544 for VP4 gene.
Data were analysed for the comparisons of prevalence rates using 2x2 tables with χ2 test, or Fisher’s exact test (two-tailed, 95% confidence intervals) when the size sample was less than 5 [Epi Info™ 22.214.171.124, CDC, Atlanta, GA, USA). Student’s test was applied for comparisons of variable values. Tests were considered significant when p <0.05.
RV detection, PAGE and G/P genotyping
Of a total of 912 stool specimens, 206 (22.6%) rotavirus-positive were detected by immunochromatographic test, of which 169 (82%) were confirmed by PAGE, showing a prevalence rate of RV infection of 18.5%. Based on the electrophoretic migration pattern by PAGE, all the isolates exhibited a profile of RV group A (4-2-3-2) (Fig. 1). Eighty-two (49%) of them showed a long (faster-moving gene segment 11) and 87 (51%) a short (slower-moving gene segment 11) pattern.
A variable detection rate of RV diarrhea was observed during the entire period studied, ranging from 2.1 to 50.1%, with the highest frequency observed during the coolest and driest months, from February to May 2007 and from February to March 2008, and a gradual decline of RV detection from March until August 2007 (Fig. 2, left axis). The comparison between the detection average rate (24%) of 2007 RV peak season (February, March and April) and that of 2008 (36%) revealed a significant increase (p = 0.008) of RV diarrhea cases between the two period.
Rotavirus caused diarrhea in children from 7 to 84 months of age, but the most of cases occurred during the first 2 years of life (median age 22 months). The male to female ratio was approximately 1:1. The majority of RV strains (84%) were detected from children who had not received any RV vaccine.
A total of 165 (97.6%) strains were characterized for the G (VP7) and P (VP4) genotype; the remaining four (2.4%) stool specimens were not genotyped because a PCR product could not be generated. The RT-PCR assays revealed that 77 (46.7%) of the 165 strains, were G2P, 63 (38.2%) G1P, four (2.4%) G8P, two (1.2%) G9P and two (1.2%) G1P throughout the study. Additionally, it was found one (0.6%) strain for each one of the following combinations: G3P, G2P, G4P and G8P. Six (3.6%) samples showed a mixed infection (G1-G2/P-P) and seven (4.2%) were G or P untypeable (NT).
The analysis of genome by PAGE revealed that most (96%) of the G2P strains had a short RNA pattern, among which could be differentiated nine electropherotypes (data not shown) overall very similar, differing only in the migration of dsRNA segments 7 to 9 upon co-electrophoresis. The G1P strains showed predominantly a long RNA pattern, and at least five slightly different electropherotypes were recognized among them (data not shown), evident for the electrophoretic mobility of the dsRNA segments 3, 7 to 9. All the G8, G9, G4 and G3 strains studied showed a long pattern of RNA migration, while the strain G2P had a short profile.
Remarkably, a changing pattern of the G/P type distribution was demonstrated throughout the RV season 2007-2008. Although G2P represented nearly half of the RV strains during the whole study, this genotype was detected almost during the entire period studied and showed a complete predominance over the others during the 2007 RV peak season. The proportion of G2P RV infected children fluctuated around an average of 80% during the first half-year studied, before falling below 10% in the remaining period (Fig. 2, right axis). This significant reduction of their prevalence coincided with an increase in the number of children infected with G1P RVs during the following months. Even though G1P was scarcely represented for most of the year 2007, it emerged to become the dominant genotype since January until March 2008 (Fig. 2, right axis), when prevailed over the others by approximately 67%.
The median age of the children infected with RV G1P and G2P was identical (21 months, p = 0.7). The unusual G8P strains detected between April and May 2007, and G8P in October, were collected mainly from children younger than 24 months of age. The other genotypes appeared sporadically (Fig. 2, right axis).
Nucleotide and deduced amino acid sequence analysis
VP7 sequence analysis
Phylogenetic analysis were based on partial sequences of VP7 gene encoding region derived from a PCR amplicon of 884-bp, comprising the hypervariable regions of VP7 protein and surface-exposed amino acids that show intergenotypic variability among prevalent human G and P genotypes (36, 41, 42). Approximately 820 nucleotides of VP7 gene from three randomly selected G1P RV isolates (strains CCS-96/2007, CCS-174/2008 and CCS-193/2008) and one G2P (strain CCS-42/2007) circulating in Caracas during the years 2007-2008 could be analyzed. Additionally two Venezuelan strains, VCE19063/2003 and VCE15377/2003, both detected during 2003 before vaccine introduction, and defined as G1P and G2P respectively in previous genotyping studies , were also included for the comparison, as shown in Fig. 3a. Genotype-specific lineages were assigned as previously described [38, 39].
In the phylogenetic tree of VP7 nucleotide sequences showed in Fig. 3a, the G1P RV strains from Caracas (2007-2008) were grouped in a unique branch, revealing more than 99.8% identities to each other at nucleotide level, and 100% at amino acid level (data not shown). These strains were clustered into the lineage I together with G1 reference strains from India reported in the years 2002, 2007 and 2008 (respectively ISO-4, NIV-0716581 and NIV-088325), for which showed the highest (≥99.5%) nucleotide identities.
Amino acid differences on the VP7 and VP4 proteins were investigated. A comparison of residues that constituted the epitopes defining the neutralization domains on VP7 and VP4 of the Venezuelan strains G1P and G2P, Rotarix® vaccine and other reference strains belonging to the same genotypes, was shown in Figs. 4 and 5. The G1P RV strains from Caracas shared the same amino acid substitutions on the VP7 with Indian strains and other global G1 RVs belonging to the same lineage I and described in geographically distant countries. On the other hand, they showed a lower identity (93%) at both nucleotide and amino acid level to the Rotarix® vaccine strain, which instead was more similar (98%) to the Venezuelan pre-vaccine strain VCE19063/2003 clustered into the lineage II (Fig. 3a). Several amino acid substitutions were observed along the deduced amino acid VP7 sequence of the three 2007/2008 G1P RVs from Caracas in comparison to Rotarix vaccine strain (data not shown), three of them (94N → S, 123S → N and 217M → T) located into the 7-1a and 7-2 neutralization domains within the hypervariable regions of VP7 protein, which were absent in most of the global reference strains belonging to G1 lineages not-I (Fig. 4a). A potential N-linked glycosylation site created by the substitution 123S → N was only observed in the strains of G1-lineage I, including the RV from Caracas, and IX (Fig. 4a). In addition, two substitutions at positions 68A → S and 72Q → R, adjacent to the conserved glycosylation motif Asn-X-Thr (residues 69–71), and others like the mutation 41Y → F, were found on the VP7 of the local strains from Caracas, but not on the Rotarix® and VCE19063/2003 strain (data not shown). No changes were visualized in Ca++ binding sites described by Aoki et al.  (data not shown).
On the other hand, the strain CCS-42/2007 of G2P genotype detected in Caracas grouped into the lineage G2-II of VP7 gene (Fig. 3a), which showed the highest nucleotide identity to the strains CMH028/2007 (99.2%) reported from Thailand, Dhaka4/2003 (98.6%) from Bangladesh, PSAL3273-F/2008 (98.1%) from Brazil and IND/KOL-32-08 from India (97.9%). A lower identity at nucleotide (92.6%) and at amino acid (95%) level was shown to the Venezuelan pre-vaccine strain VCE15377/2003 (data not shown) that clustered into the lineage III. The deduced amino acid VP7 sequence of the strain CCS-42/2007 from Caracas was very similar to most of the reference strains analyzed belonging to the lineage II, reported before and after vaccine licensure (Fig. 4b). Amino acid differences located into the 7-1a and 7-1b neutralization domains were found at residues 96, 213 and 242 in comparison with the VP7 sequence of the VCE15377/2003 (Fig. 4b).
VP4 sequence analysis
Based on the analysis of partial sequences of VP4 gene (more or less 890 nucleotides) encoding for the VP8* and part of the VP5* subunits of VP4 protein, the 2007/2008 G1P RV strains from Caracas exhibited an extremely high genetic similarity among themselves ranging from 99 to 99.8%. They were grouped within the P-III lineage together with the strain VCE19063, detected in Valencia in 2003 before vaccine introduction in Venezuela (Fig. 3b), against which demonstrated a slightly lower nucleotide and amino acid identity, 97 and 98% respectively. Furthermore, VP4 gene of the Venezuelan RV strains showed high percentages of identity (>95.3%) to reference P strains of the same lineage, isolated in different continents during the last two decades (data not shown). A lower identity (around 90%) to the Rotarix® vaccine strain, belonging to the lineage P-I, was found. In fact, the deduced amino acid VP4 sequence analysis revealed several amino acid substitutions along the surface-exposed antigenic epitopes of the VP8* portion defining neutralization domains , in comparison to the VP4 of Rotarix® vaccine strain: the substitutions 150E → D and 195N → G within the antigenic region 8-1, and 113N → D, 125S → N, 131S → R and 135N → D into the antigenic region 8-3 (Fig. 5a). Two of the three Venezuelan 2007/2008 G1P strains showed an amino acid change (I → T) at position 196, not related with escape neutralization sites (Fig. 5a) and a unique residue motif K-I-L-V at position 346–349 (data not shown).
The Venezuelan G2P strain CCS-42/2007 isolated from Caracas in 2007 segregated into the P-V lineage (Fig. 3b), sharing a nucleotide identity >98.4% with P global strains of the same period, such as the strains MMC6/2005 and Dhaka4/2003 from Bangladesh, and Py1138ASR07/2007 from Paraguay. Nucleotide and amino acid identity values, respectively of 97.4 and 99%, were found when it was compared to the VCE15377/2003 from Valencia belonging to the same lineage P-V. The Fig. 5b shows the amino acid differences located into the neutralization domains of the VP8* portion of the analyzed strains. The potential trypsin cleavage sites at arginine 240 and 246 were both conserved in all the VP4 sequences studied from Venezuelan RV G1P and G2P strains, as well as the proline at positions 68, 71, 224 and 225, and the cysteine at position 215 (data not shown), residues that have been described as highly conserved in the VP8* gene portion of human RV strains .
VP7 protein structural analysis and modeling
In order to investigate the effect of the mentioned mutations on the conformation of VP7 protein, an structural analysis comparing the VP7 model from 2007/2008 G1P post-vaccine strains from Caracas and G1P Rotarix® vaccine strain was performed. The analysis showed that both VP7 sequences were perfectly superimposable, where the amino acid substitutions were conservative and did not affect the molecular conformation (Fig. 6). Particularly, no structural modification appeared to have been generated by the three substitutions (94N → S, 123S → N and 217M → T) located into the neutralization domains along the surface-exposed antigenic epitopes on the VP7 of the Venezuelan 2007/2008 G1P strains (Fig. 6).
RV is associated with significant morbidity and mortality rates, particularly in developing countries. Studies performed in Brazil, Mexico, El Salvador and Panama have reported reduction in gastroenteritis and RV disease burden following the introduction of RV vaccines [28, 42, 43], illustrating the health benefits of these strategies. A substantial drop in deaths and hospital admissions after vaccine introduction has been also revealed by some previous works in Bolivia, Honduras and Venezuela . The RV-associated diarrhea rate showed in the present study is remarkably close (>18%) to that reported in Valencia city, not far away from Caracas, some years ago before the introduction of RV vaccines into the national immunization program [45, 46]. This result could be explained by the high percentage of unvaccinated children included in the present study. However, even though the vaccines have been applied in Venezuela since 2006, some authors reported that RV continued causing a considerable number of cases of diarrhea during 2010 . The data showed here might not reflect national and regional trends in rotavirus detection, but the results raise concerns about the factors that could compromise the success of the immunization programs in use in this context.
The efficacy of vaccines in a region can be largely influenced by multiple factors ranging from concurrent enteric infections, malnutrition, immune status, health care access, vaccine coverage rates of the population but also by the RV epidemiology . Although the present study included a relatively short period of time to capture the effects of the seasonality or vaccine introduction, it revealed the occurrence of RV peaks in cooler and drier months. This result contrasts with a previous study performed in the same geographical setting years ago where minimal or absent seasonality was found, but it agrees with the pattern of RV infection observed in Valencia during the years 1998–2002 . The median age of RV infected children included in the present study was higher than that of the children involved in the mentioned study , which could explain this difference. It agrees with the association described by other authors between the seasonality differences of RV infection in different regions of a same country and the age of children , and with a RV epidemiology systematic review/metaanalysis that reveals a trend of infection to occur in the cool, dry seasons in tropical countries, similar to that observed in temperate zones .
On the other hand, concerns exist also about the effects of the vaccination on the RV diversity, considering the selective pressure that RV vaccines could exert on viral populations, leading to changes in the distribution of the genotypes that would compromise the success of the immunization. The assessment of changing pattern of RV infections in children with diarrhea in Caracas during the RV season 2007–2008 after vaccine introduction revealed a distribution of G/P genotypes prevalence similar to that found elsewhere: G1P and G2P, two of the five most common RV genotypes found globally during the past three decades [14, 18], were the strains predominantly detected in Caracas. Furthermore, the low proportion of G9P, G3P, G1P detected ratifies the diversity of the RV genotypes circulating in this geographic region as previously described [23–26]. The occurrence of a few unusual genotype combinations as G8P, G1P, G4P and G8P was also demonstrated. G8P strains has been found sporadically in humans. In Latin America this genotype has been previously described among RV of animal origin  and in humans, specifically in Venezuela from an infant with diarrhea in 2003 , and Guatemala from an adult in 2009 . The G8P genotype is believed to have originated from animal RVs that were introduced into human populations through interspecies transmission and/or gene reassortment . Although the G8P strains isolated in Caracas were not able to spread efficiently in the population, as only few strains were detected, further studies are required to understand the transmission potentiality and the origin of these genotype combination that could take part of the complex dynamics of RV evolution.
Extensive global epidemiological surveys have demonstrated that the incidence of rotavirus individual genotypes can show a yearly fluctuation, and multiple rotavirus types can cocirculate within the same region [14, 22, 52]. A temporal variation of the predominant RV genotype occurred in Caracas, where G2P was the RV genotype leading during the 2007 following vaccine introduction in Venezuela, and G1P reemerged during the subsequent months until prevailing in 2008. Predominance of G2P field strains has been repeatedly described almost worldwide concurrently with the introduction of a universal mass vaccination program with the monovalent G1P Rotarix® vaccine [28, 52, 53], but this occurrence has been also seen in the past in countries when had not yet introduced the RV vaccination into the public sector, such as Honduras, Paraguay, Guatemala, Argentina and Bangladesh [14, 52, 54]. Studies performed in Valencia (Venezuela) during the 2003 showed a pattern of fluctuation of RV genotypes similar to that described in the present study, where the higher prevalence of G2P strains was unrelated to RV vaccination . Some authors have found findings suggesting that the relative high frequency of G2P reported in several Latin American countries can reflect a regional phenomenon or a cyclic pattern of RV strains instead of the selective pressure created by the implementation of RV monovalent vaccine [54, 55].
Comparison of the VP7 gene sequence of the Venezuelan G2P strains revealed that RV from Caracas was more similar to global strains of the lineage G2-II, described in old and recent times, than to the pre-vaccine strain from Valencia. The data do not allow determining the precise evolutionary relationship among RV isolated before and after vaccine. Nevertheless, some variations observed after analysis of the deduced amino acid sequence of VP7 protein, and to a lesser extent of VP4, of RV G2P detected after vaccine introduction in Venezuela, arouse interest, particularly two substitutions located into the 7-1a and 7-1b neutralization domains: the 96D → N substitution, which implied a radical change from a negatively charged (aspartic acid) to an uncharged (asparagine) amino acid, and the 213N → D, capable to confer a negatively charged residue (aspartic acid). Moreover, other two amino acid changes observed on the VP4 of this isolate, the 150N → D and the 89N → D, substituted both a polar/uncharged with a polar/negatively charged amino acidic residue. It remains uncertain if such electrical changes on the outer capsid viral protein could have conditioned somehow the spread of the G2P strains into the population in Caracas. The substitution 96D → N on the VP7 has been strongly related in the last decade by other authors to an abrupt increase or reemergence of G2 strains in different European, African and Asian countries, as United Kingdom, Nepal, South Africa, Taiwan, Thailand, Bangladesh and Japan [39, 56–60].
The predominance of G2P genotype in Caracas during the 2007 and its consecutive decline (from 30.1% to <10%) to be overcome by G1P during 2008, is a trend that was also described in other countries as Nicaragua some years ago before vaccine introduction . This event could have been result of differential virus fitness among susceptible and immunological protected hosts. Interestingly, no significant difference in the median age of the G1P and G2P infected children was observed in Caracas, suggesting that other factors beyond the age would be conditioning the host susceptibility to the infection. Following the introduction of the vaccine in a region, different selective pressures can be exerted on the viral populations by the homotypic and heterotypic immunity vaccine-induced. It has been extensively shown that the monovalent G1P Rotarix® vaccine used worldwide confers homotypic and heterotypic protection, providing good and sustained protection against all encountered non-G1 strains (G3, G4, G9) sharing the same P type (P) [2, 4, 62, 63]. Although it may be less efficacious against the G2 strains, an immune response can be achieved through cross-reactive epitopes on VP7 and VP4 proteins, together with VP6 or non-structural proteins epitopes . Thus, the results suggest that the combination of the weaker natural or vaccine-induced immunity against G2P and the natural fluctuations of this genotype could have favored the subsequent reemergence of G1P strains observed in Caracas at the end of 2007, as proposed by other authors .
The complete predominance of G1P strains here described over the other types during the following months (early 2008) was unexpected. Since a large proportion of children included in the present study was unvaccinated, it could be assumed they were highly susceptible to be infected by G1P strains circulating at that time. On the other hand, these strains could have acquired some selective advantage favoring their transmission. In fact, the phylogenetic analysis of the VP7 revealed that the G1P strains circulating in Caracas following RV vaccine introduction were grouped into a single genetically homogeneous clade of lineage G1-I, very close to contemporary strains described in geographically remote countries, but they were rather distantly related to Rotarix® vaccine and the Venezuelan pre-vaccine 2003 G1P strain, belonging both to the lineage G1-II. Thus, the post-vaccine G1P strains could have accumulated mutations and evolved over time into a variant that might escape from vaccine induced antibodies. It is difficult to establish when this variant arose, considering that data on the RV strains circulating in nearest pre-vaccination period are not available, but it cannot be excluded that it might have been introduced in Venezuela from other countries or were vaccine-derived.
Studies with animal and human RVs have demonstrated that neutralizing antibodies against VP7 protein play a critical role in vaccine-mediated immunological protection against RV disease, as well as the VP8* subunit of VP4 protein participates in viral infectivity and neutralization . Amino acid substitutions in outer capsid proteins could thus result in escape of viruses from neutralizing antibodies, affect viral fitness, and/or change receptor preference. The Venezuelan G1P post-vaccine strains showed amino acid differences on the 7-1a and 7-2 antigenic epitopes of the VP7 protein that have been also described by other authors on current global G1 strains . Here, the modifications observed were: i) the mutation 217M → T into the domain 7-2 that changed a methionine (non-polar/hydrophobic) with a tyrosine (polar/uncharged), probably defining the lineage G1-I; ii) the 68A → S, which determined a change from an alanine (non-polar, small) to a serine (polar/neutral), adjacent to a conserved glycosylation motif located on the residues 69 to 71; and, iii) a potential N-linked glycosylation site created by the mutation 123S → N, which was absent on the VP7 of Rotarix® vaccine strain and other G1 strains not belonging to the lineage I or IX, whose far ranging effect on the antigenicity of this epitope is unknown. Previous analysis has revealed that the amino acid residues located in the 7-1a and 7-2 antigenic epitopes are distributed fairly heterogeneously across the face of the VP7 molecule, defining the neutralization domains . Modeling structural analysis of the VP7 from Venezuelan G1P strains did not showed differences in the molecular conformation of these epitopes when compared with the VP7 from the Rotarix® vaccine strain, suggesting that the mentioned amino acid mutations did not alter the conformation of the neutralization domains. Nevertheless, although they were structurally conservative mutations, other types of studies are needed to understand whether the changes in electric charges described into the immunodominant regions could have affected the protein-antibody binding and led to loss of vaccine-induced protection.
VP4 analysis revealed that the RVs G1P analyzed, circulating in Caracas after vaccine implementation, were very similar to the G1P pre-vaccine strain VCE19063/2003 from Valencia, which was grouped in the same genetic lineage P-III, but differed from Rotarix® vaccine strain in some amino acid substitutions located into neutralization domains 8-1 and 8-3. Among them, the 196I → T was exclusively present on the VP4 of the G1P strains isolated from Caracas in 2008 (Fig. 5a) and replaced a non-polar/hydrophobic with a polar/uncharged amino acid, determining a change of electrical charge.
Antigenic variants of viruses are expected to emerge as the population immunity evolves. The emergence and later predominance of G1P strains observed in Caracas during the 2007-2008 RV season and described in this study suggests that these RV could have escaped from the immune response evocated against strains previously circulating in the region. The observed amino acid changes on the outer capsid proteins VP7 and/or VP4, which determined alterations in the electrical charges of residues located onto the neutralization domains, could have subtly affected the binding of neutralizing antibodies and conferred a selective pressure influencing the viral fitness, and favoring the transmission of the viruses. Thus, the result of forces and balances that drive rotavirus natural evolution would determine the spread of novel strains. The impact that the widespread use of RV vaccines will have on the diversity and evolution of human RVs is hardly predictable. The vaccination continues being the best control strategy against the RV diarrhea, thus it is important to achieve higher vaccine coverage rates in countries as Venezuela. Although official data are not available, it is presumed that the current coverage for oral RV vaccines in this country would surpass 75%. Therefore it is expected that vaccine impact may intensify in the next years.
This study highlights the importance of monitoring the transitions in the prevalence of genotypes and understanding of their effect on the efficacy of currently available vaccines. The genotype prevalence changes described after vaccine introduction in Caracas raise concerns, but such variations should be interpreted cautiously in the global context over time of a determined geographic area. Additionally, a continued surveillance of the genetic characteristics of RVs circulating will continue to be needed to obtain a better view of the long-term effects of vaccine introductions, to assess intra-genotype evolution that may lead to selection for strains that could escape homotypic immunity from the vaccine or that are antigenically different from those included in the vaccine formulations, and to determine the potential of their global spread. Since the immunogenicity and efficacy of the RV vaccines may be challenged by evolution of the viral genome of RV circulating, it is also important to accomplish the full genome analysis of strains collected in different time or geographic regions of a same country and know their evolutionary profile during post-licensure surveillance.
Basic Local Alignment Search Tool
Double strand RNA
Molecular Evolutionary Genetics Analysis
National Center for Biotechnology Information
Polyacrylamide gel electrophoresis
Protein Data Bank
Reverse transcriptase polymerase chain reaction
Tate JE, Burton AH, Boschi-Pinto C, Steele AD, Duque J, Parashar UD, WHO-coordinated Global Rotavirus Surveillance Network. 2008 estimate of worldwide rotavirus-associated mortality in children younger than 5 years before the introduction of universal rotavirus vaccination programmes: a systematic review and meta-analysis. Lancet Infect Dis. 2012;12:136–41.
Ruiz-Palacios GM, Pérez-Schael I, Velázquez FR, Abate H, Breuer T, Clemens SC, Cheuvart B, Espinoza F, Gillard P, Innis BL, Cervantes Y, Linhares AC, López P, Macías-Parra M, Ortega-Barría E, Richardson V, Rivera-Medina DM, Rivera L, Salinas B, Pavía-Ruz N, Salmerón J, Rüttimann R, Tinoco JC, Rubio P, Nuñez E, Guerrero ML, Yarzábal JP, Damaso S, Tornieporth N, Sáez-Llorens X, Vergara RF, Vesikari T, Bouckenooghe A, Clemens R, De Vos B, O’Ryan M, Human Rotavirus Vaccine Study Group. Safety and efficacy of an attenuated vaccine against severe rotavirus gastroenteritis. N Engl J Med. 2006;354:11–22.
Vesikari T, Matson DO, Dennehy P, Van Damme P, Santosham M, Rodriguez Z, Dallas MJ, Heyse JF, Goveia MG, Black SB, Shinefield HR, Christie CD, Ylitalo S, Itzler RF, Coia ML, Onorato MT, Adeyi BA, Marshall GS, Gothefors L, Campens D, Karvonen A, Watt JP, O’Brien KL, DiNubile MJ, Clark HF, Boslego JW, Offit PA, Heaton PM, Rotavirus Efficacy and Safety Trial (REST) Study Team. Safety and efficacy of a pentavalent human-bovine (WC3) reassortant rotavirus vaccine. N Engl J Med. 2006;354:23–33.
Salinas B, Pérez Schael I, Linhares AC, Ruiz Palacios GM, Guerrero ML, Yarzábal JP, Cervantes Y, Costa Clemens S, Damaso S, Hardt K, De Vos B. Evaluation of safety, immunogenicity and efficacy of an attenuated rotavirus vaccine, RIX4414: A randomized, placebo-controlled trial in Latin American infants. Pediatr Infect Dis J. 2005;24:807–16.
World Health Organization. Rotavirus vaccines. Wkly Epidemiol Rec. 2007;82(32):285–96.
Richardson V, Hernandez-Pichardo J, Quintanar-Solares M, Esparza-Aguilar M, Johnson B, Gomez-Altamirano CM, Parashar U, Patel R. Effect of rotavirus vaccination on death from childhood diarrhea in Mexico. N Engl J Med. 2010;362:300–5.
Paternina-Caicedo A, Parashar UD, Alvis-Guzmán N, De Oliveira LH, Castaño-Zuluaga A, Cotes-Cantillo K, Gamboa-Garay O, Coronell-Rodríguez W, De la Hoz-Restrepo F. Effect of rotavirus vaccine on childhood diarrhea mortality in five Latin American countries. Vaccine. 2015;doi: 10.1016/j.vaccine.2015.06.058.
Munos MK, Walker CL, Black RE. The effect of rotavirus vaccine on diarrhoea mortality. Int J Epidemiol. 2010;39 Suppl 1:i56–62. doi:10.1093/ije/dyq022.
Franco-Soto V, Rísquez-Parra A, Colmenares R, Casanova M, Celis D, Calderon N, Márquez L. Impact of the vaccination against rotavirus on morbidity and mortality in diarrhea. Arch Venez Pueric Pediatr. 2013;76:93–8.
Centers for Disease Control and Prevention (CDC). Progress in the Introduction of Rotavirus Vaccine - Latin America and the Caribbean, 2006–2010. Morb Mortal Wkly Rep (MMWR). 2011;60:1611–4. https://www.cdc.gov/mmwr/preview/mmwrhtml/mm6047a2.htm.
Glass RI, Parashar U, Patel M, Gentsch J, Jiang B. 2014. Rotavirus vaccines: successes and challenges. J Infect. 2014; doi:10.1016/j.jinf.2013.09.010.
Matthijnssens J, Ciarlet M, Rahman M, Attoui H, Bányai K, Estes MK, Gentsch JR, Iturriza-Gómara M, Kirkwood CD, Martella V, Mertens P, Nakagomi O, Patton JT, Ruggeri F, Saif LF, Santos N, Steyer A, Taniguchi K, Desselberger U, Van Ranst M. Recommendations for the classification of group A rotaviruses using all 11 genomic RNA segments. Arch Virol. 2008;153:1621–9.
Trojnar E, Sachsenröder J, Twardziok S, Reetz J, Otto PH, Johne R. Identification of an avian group A rotavirus containing a novel VP4 gene with a close relationship to those of mammalian rotaviruses. J Gen Virol. 2013; doi:10.1099/vir.0.047381-0.
Santos N, Hoshino Y. Global distribution of rotavirus serotypes/genotypes and its implication for the development and implementation of an effective rotavirus vaccine. Rev Med Virol. 2005;15:29–56.
World Health Organization. Global Rotavirus Information and Surveillance Bulletin. Reporting period: January to December 2010. WHO Rotavirus Surveillance Websites.2011;4. http://www.who.int/immunization/sage/3_Final_RV_bulletin_Jan_Dec_2010_Data_nov11.pdf?ua=1.
Iturriza-Gómara M, Dallman T, Bányai K, Böttiger B, Buesa J, Diedrich S, Fiore L, Johansen K, Koopmans M, Korsun N, Koukou D, Kroneman A, László B, Lappalainen M, Maunula L, Marques AM, Matthijnssens J, Midgley S, Mladenova Z, Nawaz S, Poljsak-Prijatelj M, Pothier P, Ruggeri FM, Sanchez-Fauquier A, Steyer A, Sidaraviciute-Ivaskeviciene I, Syriopoulou V, Tran AN, Usonis V, VAN Ranst M, DE Rougemont A, Gray J. Rotavirus genotypes co-circulating in Europe between 2006 and 2009 as determined by EuroRotaNet, a pan-European collaborative strain surveillance network. Epidemiol Infect. 2011; doi: 10.1017/S0950268810001810.
Rahman M, Matthijnssens J, Yang X, Delbeke T, Arijs I, Taniguchi K, et al. Evolutionary history and global spread of the emerging G12 human rotaviruses. J Virol. 2007;81:2382–90.
Bányai K, László B, Duque J, Steele AD, Nelson EA, Gentsch JR, Parashar UD. Systematic review of regional and temporal trends in global rotavirus strain diversity in the pre rotavirus vaccine era: insights for understanding the impact of rotavirus vaccination programs. Vaccine. 2012; doi:10.1016/j.vaccine.2011.09.111.
Luchs A, Timenetsky Mdo C. G8P rotaviruses isolated from Amerindian children in Mato Grosso do Sul, Brazil, during 2009: close relationship of the G and P genes with those of bovine and bat strains. J Gen Virol. 2014;95:627–41.
Da Silva MF, Tort LF, Goméz MM, Assis RM, Volotão Ede M, de Mendonc¸ a, MC, Bello G, Leite JP. VP7 gene of human rotavirus A genotype G5: phylo-genetic analysis reveals the existence of three different lineages worldwide. JMed Virol. 2011;83:357–66.
Kang G, Desai R, Arora R, Chitamabar S, Naik TN, Krishnan T, Deshpande J, Gupte MD, Venkatasubramaniam S, Gentsch JR, Parashar UD, Indian Rotavirus Strain Surveillance Network, Mathew A, Anita Sr, Ramani S, Sowmynarayanan TV, Moses PD, Agarwal I, Simon A, Bose A, Arora R, Chhabra P, Fadnis P, Bhatt J, Shetty SJ, Saxena VK, Mathur M, Jadhav A, Roy S, Mukherjee A, Singh NB. Diversity of circulating rotavirus strains in children hospitalized with diarrhea in India, 2005–2009. Vaccine. 2013;31:2879–83.
Linhares AC, Stupka JA, Ciapponi A, Bardach AE, Glujovsky D, Aruj PK, Mazzoni A, Rodriguez JA, Rearte A, Lanzieri TM, Ortega-Barria E, Colindres R. Burden and typing of rotavirus group A in Latin America and the Caribbean: systematic review and meta-analysis. Rev Med Virol. 2011; doi:10.1002/rmv.682.
Pérez-Schael I, González R, Fernández R, Alfonso E, Inaty D, Boher Y, Sarmiento L. Epidemiological features of rotavirus infection in Caracas, Venezuela: Implications for rotavirus immunization programs. J Med Virol. 1999;59:520–6.
White L, García D, Boher Y, Blanco M, Pérez M, Romer H, Flores J, Pérez-Schael I. Temporal distribution of human rotavirus serotypes 1, 2, 3 and 4 in Venezuelan children with gastroenteritis during 1979–1989. J Med Virol. 1991;34:79–84.
González R, Rivero L. Genetic diversity of rotavirus group a: correlation between G3 type and severity of the infection. Valencia, Venezuela. Invest Clin. 2013;54:34–46.
Vizzi E, Piñeros O, González GG, Zambrano JL, Ludert JE, Liprandi F. Genotyping of human rotaviruses circulating among children with diarrhea in Valencia, Venezuela. J Med Virol. 2011; doi:10.1002/jmv.22211.
Gurgel RQ, Cuevas LE, Vieira SC, Barros VC, Fontes PB, Salustino EF, Nakagomi O, Nakagomi T, Dove W, Cunliffe N, Hart CA. Predominance of rotavirus PG2 in a vaccinated population, Brazil. Emerg Infect Dis. 2007;13:1571–3.
Nakagomi T, Cuevas LE, Gurgel RG, Elrokhsi SH, Belkhir YA, Abugalia M, Dove W, Montenegro FM, Correia JB, Nakagomi O, Cunliffe NA, Hart CA. Apparent extinction of non-G2 rotavirus strains from circulation in Recife, Brazil, after the introduction of rotavirus vaccine. 2008. Arch Virol. 2008; doi:10.1007/s00705-007-0028-z.
Kirkwood CD, Boniface K, Barnes GL, Bishop RF. Distribution of rotavirus genotypes after introduction of rotavirus vaccines, Rotarix® and RotaTeq®, into the National Immunization Program of Australia. Pediatr Infect Dis J. 2011; doi:10.1097/INF.0b013e3181fefd90.
Rodrigues F, Iturriza-Gómara M, Marlow R, Gray J, Nawaz S, Januário L, Finn A. The evolving epidemiology of rotavirus gastroenteritis in central Portugal with modest vaccine coverage. J Clin Virol. 2013; doi:10.1016/j.jcv.2012.10.016.
Hull JJ, Teel EN, Kerin TK, Freeman MM, Esona MD, Gentsch JR, Cortese MM, Parashar UD, Glass RI, Bowen MD. United States rotavirus strain surveillance from 2005 to 2008: genotype prevalence before and after vaccine introduction. 2011. National Rotavirus Strain Surveillance System. Pediatr Infect Dis J. 2011; doi:10.1097/INF.0b013e3181fefd.
Pitzer VE, Patel MM, Lopman BA, Viboud C, Parashar UD, Grenfell BT. Modeling rotavirus strain dynamics in developed countries to understand the potential impact of vaccination on genotype distributions. Proc Natl Acad Sci USA. 2011; doi:10.1073/pnas.1110507108.
Herring AJ, Inglis NF, Ojeh CK, Snodgrass DR, Menzies JD. Rapid diagnosis of rotavirus infection by direct detection of viral nucleic acid in silver-stained polyacrylamide gels. J Clin Microbiol. 1982;16:473–7.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Pieper U, Webb BM, Dong GQ, Schneidman-Duhovny D, Fan H, Kim SJ, Khuri N, Spill YG, Weinkam P, Hammel M, Tainer JA, Nilges M, Sali A. MODBASE, a database of annotated comparative protein structure models and associated resources. Nucleic Acids Res. 2014;42:D336–346.
Aoki ST, Settembre EC, Trask SD, Greenberg HB, Harrison SC, Dormitzer PR. Structure of rotavirus outer-layer protein VP7 bound with a neutralizing Fab. Science. 2009;324:1444–7.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE. UCSF Chimera--a visualization system for exploratory research and analysis. J Comput Chem. 2004;25:1605–12.
Parra GI, Bok K, Martinez V, Russomando G, Gomez J. Molecular characterization and genetic variation of the VP7 gene of human rotaviruses isolated in Paraguay. J Med Virol. 2005;77:579–86.
Page NA, Steele AD. Antigenic and genetic characterization of serotype G2 human rotavirus strains from South Africa from 1984 to 1998. J Med Virol. 2004;72:320–7.
Dormitzer PR, Sun ZY, Wagner G, Harrison SC. The rhesus rotavirus VP4 sialic acid binding domain has a galectin fold with a novel carbohydrate binding site. EMBO J. 2002;21:885–97.
Silva MF, Rose TL, Gómez MM, Carvalho-Costa FA, Fialho AM, Assis RM, Andrade Jda S, Volotão Ede M, Leite JP. G1P species A rotavirus over 27 years--pre- and post-vaccination eras--in Brazil: full genomic constellation analysis and no evidence for selection pressure by Rotarix® vaccine. Infect Genet Evol. 2015; doi:10.1016/j.meegid.2014.12.030.
Gurgel RG1, Bohland AK, Vieira SC, Oliveira DM, Fontes PB, Barros VF, Ramos MF, Dove W, Nakagomi T, Nakagomi O, Correia JB, Cunliffe N, Cuevas LE. Incidence of rotavirus and all-cause diarrhea in northeast Brazil following the introduction of a national vaccination program. Gastroenterology. 2009; doi:10.1053/j.gastro.2009.07.046.
Desai R, Oliveira LH, Parashar UD, Lopman B, Tate JE, Patel MM. Reduction in morbidity and mortality from childhood diarrhoeal disease after species A rotavirus vaccine introduction in Latin America - a review. Mem Inst Oswaldo Cruz. 2011;106:907–11.
De Oliveira LH, Giglio N, Ciapponi A, García Martí S, Kuperman M, Sanwogou NJ, Ruiz-Matus C, Marinho de Sousa MF. Temporal trends in diarrhea-related hospitalizations and deaths in children under age 5 before and after the introduction of the rotavirus vaccine in four Latin American countries. Vaccine. 2013; doi:10.1016/j.vaccine.2013.05.065.
Salinas B, González G, González R, Escalona M, Materán M, Schael IP. Epidemiologic and clinical characteristics of rotavirus disease during five years of surveillance in Venezuela. Pediatr Infect Dis J. 2004;23:S161–7.
González GG, Liprandi F, Ludert JE. Molecular epidemiology of enteric viruses in children with sporadic gastroenteritis in Valencia, Venezuela. J Med Virol. 2011; doi:10.1002/jmv.22185.
Patel M, Shane AL, Parashar UD, Jiang B, Gentsch JR, Glass RI. Oral rotavirus vaccines: how well will they work where they are needed most? J Infect Dis. 2009; doi:10.1086/605035.
Pérez-Schael I., González R, Salinas B. Severity and age of rotavirus diarrhea, but not socioeconomic conditions, are associated with rotavirus seasonality in Venezuela. J Med Virol. 2009; doi:10.1002/jmv.21420.
Levy K, Hubbard AE, Eisenberg JN. 2009. Seasonality of rotavirus disease in the tropics: a systematic review and meta-analysis. Int J Epidemiol. 2009; doi:10.1093/ije/dyn260.
Matthijnssens J, Potgieter CA, Ciarlet M, Parreño V, Martella V, Bányai K, Garaicoechea L, Palombo EA, Novo L, Zeller M, Arista S, Gerna G, Rahman M, Van Ranst M. Are human P rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla? J Virol. 2009; doi:10.1128/JVI.02246-08.
Gautam R, Mijatovic-Rustempasic S, Roy S, Esona MD, Lopez B, Mencos Y, Rey-Benito G, Bowen MD. Full genomic characterization and phylogenetic analysis of a zoonotic human G8P rotavirus strain detected in a sample from Guatemala. Infect Genet Evol. 2015; doi:10.1016/j.meegid.2015.05.004.
Leite JP, Carvalho-Costa FA, Linhares AC. Group A rotavirus genotypes and the ongoing Brazilian experience - A Review. Mem Inst Oswaldo Cruz. 2008;103:745–53.
Matthijnssens J, Zeller M, Heylen E, De Coster S, Vercauteren J, Braeckman T, Van Herck K, Meyer N, Pirçon JY, Soriano-Gabarro M, Azou M, Capiau H, De Koster J, Maernoudt AS, Raes M, Verdonck L, Verghote M, Vergison A, Van Damme P, Van Ranst M; RotaBel study group. Higher proportion of G2P rotaviruses in vaccinated hospitalized cases compared with unvaccinated hospitalized cases, despite high vaccine effectiveness against heterotypic G2P rotaviruses. Clin. Microbiol. Infect. 2014; doi:10.1111/1469-0691.12612.
Esteban LE, Rota RP, Gentsch JR, Jiang B, Esona M, Glass RI, Glikmann G, Castello AA. 2010. Molecular epidemiology of group A rotavirus in Buenos Aires, Argentina 2004-2007: reemergence of G2P and emergence of G9P strains. J Med Virol. 2010; doi:10.1002/jmv.21745.
Luchs A, Cilli A, Morillo SG, de Cassia Compagnoli Carmona R, do Carmo SampaioTavares Timenetsky M. Rotavirus in adults Brazil, 2004–2011: G2P domi-nance and potential impact on vaccination. Braz J Infect Dis. 2014;18:59–69.
Iturriza-Gomara M, Cubitt D, Desselberger U, Gray J. Amino acid substitution within the VP7 protein of G2 rotavirus strains associated with failure to serotype. J Clin Microbiol. 2001;39:3796–8.
Ahmed K, Ahmed S, Mitui MT, Rahman A, Kabir L, Hannan A, Nishizono A, Nakagomi O. Molecular characterization of VP7 gene of human rotaviruses from Bangladesh. Virus Genes. 2010;40:347–56.
Trinh QD, Pham NT, Nguyen TA, Phan TG, Yan H, le Hoang P, Khamrin P, Maneekarn N, Li Y, Okitsu S, Mizuguchi M, Ushijima H. Sequence analysis of the VP7 gene of human rotaviruses G2 and G4 isolated in Japan, China, Thailand, and Vietnam during 2001-2003. J Med Virol. 2010;82:878–85.
Doan YH, Nakagomi T, Cunliffe NA, Pandey BD, Sherchand JB, Nakagomi O. The occurrence of amino acid substitutions D96N and S242N in VP7 of emergent G2P rotaviruses in Nepal in 2004–2005: a global and evolutionary perspective. Arch Virol. 2011;156:1969–78.
Do LP, Nakagomi T, Doan YH, Kitahori Y, Nakagomi O. 2014. Molecular evolution of the VP7 gene of Japanese G2 rotaviruses before vaccine introduction. Arch Virol. 2014; doi:10.1007/s00705-013-1804-6.
Espinoza F, Bucardo F, Paniagua M, Svensson L, Hallander HO, Bondeson K. Shifts of rotavirus G and P types in Nicaragua–2001-2003. Pediatr Infect Dis J. 2006;25:1078–80.
Chandran A, Fitzwater S, Zhen A, Santosham M. Prevention of rotavirus gastroenteritis in infants and children: rotavirus vaccine safety, efficacy, and potential impact of vaccines. Biologics. 2010;9:213–29.
Vesikari T, Karvonen A, Prymula R, Schuster V, Tejedor JC, Cohen R, Meurice F, Han HH, Damaso S, Bouckenooghe A. Efficacy of human rotavirus vaccine against rotavirus gastroenteritis during the first 2 years of life in European infants: randomised, double-blind controlled study. Lancet. 2007;370:1757–63.
Grimwood K, Kirkwood CD. Human rotavirus vaccines: too early for the strain to tell. Lancet. 2008;371:1144–5.
Pitzer VE, Bilcke J, Heylen E, Crawford FW, Callens M, De Smet F, Van Ranst M, Zeller M, Matthijnssens J. Did Large-Scale Vaccination Drive Changes in the Circulating Rotavirus Population in Belgium? Sci Rep. 2015; doi: 10.1038/srep18585.
Mackow ER, Shaw RD, Matsui SM, Vo PT, Dang MN, Greenberg HB. The rhesus rotavirus gene encoding protein VP3: location of amino acids involved in homologous and heterologous rotavirus neutralization and identification of a putative fusion region. Proc Natl Acad Sci U S A. 1988;85:645–9.
Zeller M, Patton JT, Heylen E, De Coster S, Ciarlet M, Van Ranst M, Matthijnssens J. Genetic analyses reveal differences in the VP7 and VP4 antigenic epitopes between human rotaviruses circulating in Belgium and rotaviruses in Rotarix and RotaTeq. J Clin Microbiol. 2012; doi:10.1128/JCM.05590-11.
Matthijnssens J, Ciarlet M, McDonald S, Attoui H, Bányai K, Brister J, Buesa J, Esona M, Estes M, Gentsch J, Iturriza-Gómara M, Johne R, Kirkwood K, Martella V, Mertens P, Nakagomi O, Parreño V, Rahman M, Ruggeri F, Saif L, Santos N, Steyer A, Taniguchi K, Patton J, Desselberger U, Van Rans M. Uniformity of rotavirus strain nomenclature proposed by The Rotavirus Classification Working Group (RCWG). Arch Virol. 2011;156:1397–413.
McDonald SM, Matthijnssens J, McAllen JK, Hine E, Overton L, Wang S, Lemey P, Zeller M, Van Ranst M, Spiro DJ, Patton JT. Evolutionary dynamics of human rotaviruses: balancing reassortment with preferred genome constellations. PLoS Pathog. 2009;5:e1000634.
The authors are in debt with the infants and their families who participated in this study, and they thank all professional staff from the Policlínica Metropolitana’s Laboratory from Caracas, particularly to Mrs. Larissa Deninzon for her assistance during sample and data collection.
This study was partially supported by GRANT LOCTI (Centro Médico Loira, C.A.) and by regular funds n. 250 from the Instituto Venezolano de Investigaciones Científicas (IVIC), Venezuela.
Availability of data and materials
The authors have confirmed that the nucleotide sequence data for the RV isolates reported in this study have been deposited in the GenBank nucleotide sequence database under the accession numbers KY039372-KY039373 and KY053848-KY053851 for VP7 gene, and KY056539-KY056544 for VP4 gene.
EV, LN and JAS conceived the study. LN, EV and MDO were engaged in sample and data collection. EV participated in the design of the study and experimental protocols, adapted protocols for the molecular assays, supervised the laboratory procedures and PCR quality control, realized the data collection from DNA sequencing assay, the analysis results, and wrote the initial draft of the manuscript. MDO and OAP performed the molecular analysis. EV, RF, AC and FL contributed to data interpretation. JLZ, AC, RF, LN, JAS and FL helped check the draft manuscript. All authors read and approved the final manuscript.
LN is an employee of GlaxoSmithKline since 2010. The others authors have no conflict of interest to declare.
Consent for publication
All the participants are consent for publication.
Ethics approval and consent to participate
All the assays performed in this research complied with the current laws of Venezuela and were approved by the IVIC Ethics Committee.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.