Use of sequence analysis of the VP4 gene to classify recent Vietnamese rotavirus isolates

Antigenic epitope analysis of Pakistani G3 and G9 rotavirus strains compared to vaccine strains revealed multiple amino acid differences

Rotaviruses belong to genotype VP4-P[8] are a significant cause of severe loose diarrhea in infants and young children. In the present study, we characterised the complete genome of three of the Pakistani P[8]b RVA strains by Illumina HiSeq sequencing technology to determine the complete genotype constellation providing insight into the evolutionary dynamics of their genes using maximum likelihood analysis. The maximum genomic sequences of our study strains were similar to more recent human Wa-Like G1P[8]a, G3P[8]a, G4P[6], G4P[8], G9P[4], G9P[8]a, G11P[25],G12P[8]a and G12P[6] strains circulating around the world. Therefore, strains PAK274, PAK439 and PAK624 carry natively distinctive VP4 gene with universally common human Wa-Like genetic backbone. Comparing our study P[8]b strains with vaccines strains RotarixTM and RotaTeqTM, multiple amino acid differences were examined between vaccine virus antigenic epitopes and Pakistani isolates. Over time, these differences may result in the selection for strains that will escape the vaccine-induced RVA-neutralizing-antibody effect.

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.

High genetic diversity of species A rotaviruses detected in swine farms in Chile

Rotavirus A is one of the main causative agents of diarrhoea in lactating and weaned pigs worldwide. Its impact in the swine industry is well documented. However, in Chile, the current epidemiological status of rotavirus on porcine farms is unknown. This study evaluated the current epidemiologic status of rotavirus A infection in Chile using on-farm detection techniques, electrophoretic confirmation, genotyping and phylogenetic clustering by analysis of partial sequences of VP4 and VP7 genes. Rotavirus A was detected in four out of five farms with an overall prevalence of 17.7 % in diarrhoeic pigs. The average age of diarrhoea onset was at 32±6.2 days, corresponding to weaning pigs, and rotavirus was not detected in lactating piglets. Molecular characterization indicated that genotypes G5, G3, P[7] and P[13] are currently the most widely represented on these pigs farms. The phylogenetic analysis showed that farms shared similar G types (VP7), which might denote a common origin. Meanwhile, [P] types (VP4) showed considerable genetic diversity, and this might represent a high rate of reassortment of this genetic segment in rotavirus circulating in the researched area. These findings demonstrate the importance of considering both the geographical and production factors to accurately determine rotavirus prevalence status at the national level, and have relevant implications in determining effective strategies for rotavirus infection control on porcine farms.

Analysis of complete genome sequences of G9P[19] rotavirus strains from human and piglet with diarrhea provides evidence for whole-genome interspecies transmission of nonreassorted porcine rotavirus

Whole genomes of G9P[19] human (RVA/Human-wt/THA/CMH-S070-13/2013/G9P[19]) and porcine (RVA/Pig-wt/THA/CMP-015-12/2012/G9P[19]) rotaviruses concurrently detected in the same geographical area in northern Thailand were sequenced and analyzed for their genetic relationships using bioinformatic tools. The complete genome sequence of human rotavirus RVA/Human-wt/THA/CMH-S070-13/2013/G9P[19] was most closely related to those of porcine rotavirus RVA/Pig-wt/THA/CMP-015-12/2012/G9P[19] and to those of porcine-like human and porcine rotaviruses reference strains than to those of human rotavirus reference strains. The genotype constellation of G9P[19] detected in human and piglet were identical and displayed as the G9-P[19]-I5-R1-C1-M1-A8-N1-T1-E1-H1 genotypes with the nucleotide sequence identities of VP7, VP4, VP6, VP1, VP2, VP3, NSP1, NSP2, NSP3, NSP4, and NSP5 at 99.0%, 99.5%, 93.2%, 97.7%, 97.7%, 85.6%, 89.5%, 93.2%, 92.9%, 94.0%, and 98.1%, respectively. The findings indicate that human rotavirus strain RVA/Human-wt/THA/CMH-S070-13/2013/G9P[19] containing the genome segments of porcine genetic backbone is most likely a human rotavirus of porcine origin. Our data provide an evidence of interspecies transmission and whole-genome transmission of nonreassorted G9P[19] porcine RVA to human occurring in nature in northern Thailand.

Identification of OP354-like human rotavirus strains with subtype P[8]b in Ghanaian children with diarrhoea

BACKGROUND

Rotaviruses with the P[8] genotype have been associated with majority of infections. Recent improvements in molecular diagnostics have delineated the P[8] genotype into P[8]a and P[8]b subtypes. P[8]a is the previously known P[8] genotype which is common whilst P[8]b subtype also known as OP354-like strain is genetically distinct, rarely detected and reported from a few countries. In a previous study, the P-types could not be determined for 80 RVA-positive samples by conventional RT-PCR genotyping methods with the recommended pool of P-genotype specific primers used in the WHO Regional Rotavirus Reference Laboratory in Ghana. The present study employed sequence-dependent cDNA amplification method to genotype previously non-typeable P-types.

METHODS

Viral RNAs were extracted and rotavirus VP4 genes amplified by one step RT-PCR using gene specific primers. PCR amplicons were purified, sequenced and sequences aligned with cognate gene sequences available in GenBank using the ClustalW algorithm. Phylogenetic analysis was performed using the Neighbour-Joining method in MEGA v6.06 software. Phylogenetic tree was statistically supported by bootstrapping with 1000 replicates, and distances calculated using the Kimura-2 parameter model.

RESULTS

Of the 80 RVA-positive samples, 57 were successfully sequenced and characterized. Forty-eight of these were identified as P[8] strains of which 5 were characterized as the rare P[8]b subtype. Phylogenetic analysis of the VP8* fragment of the VP4 genes of these P[8]b strains revealed a close relationship with prototype OP354-like P[8]b strain and P[8]b strains of Russian and South African P[8]b origin.

CONCLUSION

The study highlights the importance of regularly updating the primers employed for molecular typing of rotaviruses.

Molecular characterisation of wild-type G1P[8] and G3P[8] rotaviruses isolated in Vietnam 2008 during a vaccine trial

Rotavirus vaccines work better in developed countries than in developing countries, leading to the question of whether the circulating strains are different in these two settings. In 2008, a clinical trial of the pentavalent rotavirus vaccine was performed in Nha Trang, Vietnam, in which the efficacy was reported to be 64 %. Although samples were collected independently from the clinical trial, we examined faecal specimens from children hospitalised for rotavirus diarrhoea and found that G3P[8] and G1P[8] were co-dominant at the time of the clinical trial. The aim of this study was to explore whether they were divergent from the strains circulating in the developed countries where the vaccine efficacy is high. Two G3P[8] and two G1P[8] strains that were regarded as representatives based on their electropherotypes were selected for full-genome sequencing. The genotype constellation was G1/G3-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1. All but the VP4 genes, one of which belonged to the emerging P[8]b genotype (OP354-like VP4), clustered into one or more lineages/alleles with the strains circulating in developed countries, with ≥97.5 % nucleotide sequence identity. Additionally, 10 G1 and 12 G3 VP7 sequences as well as 31 VP4 sequences were determined. No amino acid differences were observed between the Vietnamese strains and strains in the developed countries that were likely to have affected the neutralisation specificity of their VP7 and VP4. In conclusion, apart from prevalent P[8]b VP4, virtually no differences were observed between the predominant strains circulating in Vietnam at the time of the clinical trial and the strains in the developed countries; hence, the lower vaccine efficacy was more likely to be due to factors other than strain divergence.

Emerging OP354-Like P[8] Rotaviruses Have Rapidly Dispersed from Asia to Other Continents

The majority of human group A rotaviruses possess the P[8] VP4 genotype. Recently, a genetically distinct subtype of the P[8] genotype, also known as OP354-like P[8] or lineage P[8]-4, emerged in several countries. However, it is unclear for how long the OP354-like P[8] gene has been circulating in humans and how it has spread. In a global collaborative effort 98 (near-)complete OP354-like P[8] VP4 sequences were obtained and used for phylogeographic analysis to determine the viral migration patterns. During the sampling period, 1988-2012, we found that South and East Asia acted as a source from which strains with the OP354-like P[8] gene were seeded to Africa, Europe, and North America. The time to the most recent common ancestor (TMRCA) of all OP354-like P[8] genes was estimated at 1987. However, most OP354-like P[8] strains were found in three main clusters with TMRCAs estimated between 1996 and 2001. The VP7 gene segment of OP354-like P[8] strains showed evidence of frequent reassortment, even in localized epidemics, suggesting that OP354-like P[8] genes behave in a similar manner on the evolutionary level as other P[8] subtypes. The results of this study suggest that OP354-like P[8] strains have been able to disperse globally in a relatively short time period. This, in combination with a relatively large genetic distance to other P[8] subtypes, might result in a lower vaccine effectiveness, underscoring the need for a continued surveillance of OP354-like P[8] strains, especially in countries where rotavirus vaccination programs are in place.

Dominance of emerging G9 and G12 genotypes and polymorphism of VP7 and VP4 of rotaviruses from Bhutanese children with severe diarrhea prior to the introduction of vaccine

A prospective study was performed to determine the molecular characteristics of rotaviruses circulating among children aged <5 years in Bhutan. Stool samples were collected from February 2010 through January 2011 from children who attended two tertiary care hospitals in the capital Thimphu and the eastern regional headquarters, Mongar. The samples positive for rotavirus was mainly comprised genotype G1, followed by G12 and G9. The VP7 and VP4 genes of all genotypes clustered mainly with those of neighboring countries, thereby indicating that they shared common ancestral strains. The VP7 gene of Bhutanese G1 strains belonged to lineage 1c, which differed from the lineages of vaccine strains. Mutations were also identified in the VP7 gene of G1 strains, which may be responsible for neutralization escape strains. Furthermore, we found that lineage 4 of P[8] genotype differed antigenically from the vaccine strains, and mutations were identified in Bhutanese strains of lineage 3. The distribution of rotavirus genotypes varies among years, therefore further research is required to determine the distribution of rotavirus strain genotypes in Bhutan.

Novel porcine-like human G26P[19] rotavirus identified in hospitalized paediatric diarrhoea patients in Ho Chi Minh City, Vietnam

During a hospital-based diarrhoeal disease study conducted in Ho Chi Minh City, Vietnam from 2009 to 2010, we identified four symptomatic children infected with G26P[19] rotavirus (RV) – an atypical variant that has not previously been reported in human gastroenteritis. To determine the genetic structure and investigate the origin of this G26P[19] strain, the whole genome of a representative example was characterized, revealing a novel genome constellation: G26–P[19]–I5–R1–C1–M1–A8–N1–T1–E1–H1. The genome segments were most closely related to porcine (VP7, VP4, VP6 and NSP1) and Wa-like porcine RVs (VP1–3 and NSP2–5). We proposed that this G26P[19] strain was the product of zoonotic transmission coupled with one or more reassortment events occurring in human and/or animal reservoirs. The identification of such strains has potential implications for vaccine efficacy in south-east Asia, and outlines the utility of whole-genome sequencing for studying RV diversity and zoonotic potential during disease surveillance.

Genetic analyses reveal differences in the VP7 and VP4 antigenic epitopes between human rotaviruses circulating in Belgium and rotaviruses in Rotarix and RotaTeq

Two live-attenuated rotavirus group A (RVA) vaccines, Rotarix (G1P[8]) and RotaTeq (G1-G4, P[8]), have been successfully introduced in many countries worldwide, including Belgium. The parental RVA strains used to generate the vaccines were isolated more than 20 years ago in France (G4 parental strain in RotaTeq) and the United States (all other parental strains). At present, little is known about the relationship between currently circulating human RVAs and the vaccine strains. In this study, we determined sequences for the VP7 and VP4 outer capsid proteins of representative G1P[8], G2P[4], G3P[8], G4P[8], G9P[8], and G12P[8] RVAs circulating in Belgium during 2007 to 2009. The analyses showed that multiple amino acid differences existed between the VP7 and VP4 antigenic epitopes of the vaccine viruses and the Belgian isolates, regardless of their G and P genotypes. However, the highest variability was observed among the circulating G1P[8] RVA strains and the G1 and P[8] components of both RVA vaccines. In particular, RVA strains of the P[8] lineage 4 (OP354-like) showed a significant number of amino acid differences with the P[8] VP4 of both vaccines. In addition, the circulating Belgian G3 RVA strains were found to possibly possess an extra N-linked glycosylation site compared to the G3 RVA vaccine strain of RotaTeq. These results indicate that the antigenic epitopes of RVA strains contained in the vaccines differ substantially from those of the currently circulating RVA strains in Belgium. Over time, these differences might result in selection for strains that escape the RVA neutralizing-antibody pressure induced by vaccines.

The emergence of rotavirus G12 and the prevalence of enteric viruses in hospitalized pediatric diarrheal patients in southern Vietnam

Diarrhea is a major cause of childhood morbidity and mortality in developing countries, and the majority of infections are of viral etiology. We aimed to compare the etiological prevalence of the major enteric viruses in an urban and a rural setting in southern Vietnam. We simultaneously screened fecal specimens from 362 children in Ho Chi Minh City and Dong Thap province that were hospitalized with acute diarrhea over a 1-month-long period for four viral gastrointestinal pathogens. Rotavirus was the most common pathogen identified, but there was a differential prevalence of rotavirus and norovirus between the urban and rural locations. Furthermore, rotavirus genotyping and phylogenetic analysis again differentiated the genotypes by the sampling location. Our data show a disproportional distribution of enteric viral pathogens in urban and rural locations, and we provide evidence of continual importation of new rotavirus strains into southern Vietnam and report the emergence of rotavirus genotype G12.

Prevalence of human rotavirus genotypes in Wuhan, China, during 2008-2011: changing trend of predominant genotypes and emergence of strains with the P[8]b subtype of the VP4 gene

Hospital-based surveillance of rotavirus genotypes was conducted in Wuhan, China, between March 2008 and May 2011. The detection rates of group A rotavirus were 24.6% (458/1859) and 12.1% (96/795) in children and adults, respectively, with diarrhea. Among the 554 positive specimens, the most frequent genotype was G3P[8] (57.9%), followed by G1P[8] (29.4%). Compared with previous studies in Wuhan (2000-2008), the relative frequency of G3P[8] has been decreasing year by year, while the predominant genotype G3 shifted to G1 in 2011. In the present study, a rare P[8]b subtype of the VP4 gene (OP354-like P[8]) was identified in nine strains. Full-length sequences of VP7, VP4, VP6 and NSP4 genes of two G9P[8]b strains (RVA/Human-wt/CHN/E1545/2009/G9P[8]b and RVA/Human-wt/CHN/Z1108/2008/G9P[8]b) were determined for phylogenetic analysis. The four genes of these strains were closely related to one another, and the G9-VP7 genes of these strains belonged to lineage III, which contains globally spreading G9 rotaviruses. The full-length sequence of VP4 gene segments of the P[8]b strains in Wuhan clustered with those of P[8]b strains in Vietnam, Russia and Belgium, while they were distinct from those of the OP354 strain from Malawi and Bangladeshi strains. The VP6 and NSP4 genes of two P[8]b strains belonged to the I1 and E1 genotype, respectively, and clustered with those of strains belonging to Wa-like human rotaviruses from various Asian countries. These findings indicate the changing epidemiologic trend of rotavirus genotypes in Wuhan, i.e., the shift of the predominant type from G3 to G1 and the emergence of P[8]b strains genetically related to those distributed in other Asian countries.

Diverse origin of P[19] rotaviruses in children with acute diarrhea in Taiwan: Detection of novel lineages of the G3, G5, and G9 VP7 genes

We previously reported the detection of genotype P[19] rotavirus strains from children hospitalized with acute dehydrating diarrhea during a 5-year surveillance period in Taiwan. The characterization of five P[19] strains (0.4% of all typed), including three G3P[19], a novel G5P[19], and a unique G9P[19] genotype is described in this study. Phylogenetic analysis of the VP4, VP7, VP6, and NSP4 genes was performed, which demonstrated novel lineages for respective genotypes of the VP4 and the VP7 genes. The sequence similarities of the P[19] VP4 gene among Taiwanese human strains was higher (nt, 91.5-96.2%; aa, 93.7-97.6%) than to other P[19] strains (nt, 83.5-86.6%; aa, 89.4-94.1%) from different regions of the world. The VP7 gene of the three G3P[19] Taiwanese strains shared up to 93.4% nt and 97.5% aa identity to each other but had lower similarity to reference strain sequences available in GenBank (nt, <90.1%; aa, <95.6%). Similarly, the VP7 gene of the novel G5P[19] strain was only moderately related to the VP7 gene of reference G5 strains (nt, 82.2-87.3%; aa, 87.0-93.1%), while the VP7 gene of the single G9P[19] strain was genetically distinct from other known human and animal G9 rotavirus strains (nt, ≤ 92.0%; aa, ≤ 95.7%). Together, these findings suggest that the Taiwanese P[19] strains originated by independent interspecies transmission events. Synchronized surveillance of human and animal rotaviruses in Taiwan should identify possible hosts of these uncommon human rotavirus strains.

Molecular characterization of group A rotaviruses detected in children with gastroenteritis in Ireland in 2006–2009

Community and hospital-acquired cases of human rotavirus are responsible for millions of gastroenteritis cases in children worldwide, chiefly in developing countries, and vaccines are now available. During surveillance activity for human rotavirus infections in Ireland, between 2006 and 2009, a total of 420 rotavirus strains were collected and analysed. Upon either PCR genotyping and sequence analysis, a variety of VP7 (G1–G4 and G9) and VP4 (P[4], P[6], P[8] and P[9]) genotypes were detected. Strains G1P[8] were found to be predominant throughout the period 2006–2008, with slight fluctuations seen in the very limited samples available in 2008–2009. Upon either PCR genotyping and sequence analysis of selected strains, the G1, G3 and G9 viruses were found to contain E1 (Wa-like) NSP4 and I1 VP6 genotypes, while the analysed G2 strains possessed E2 NSP4 and I2 VP6 genotypes, a genetic make-up which is highly conserved in the major human rotavirus genogroups Wa- and Kun-like, respectively. Upon sequence analysis of the most common VP4 genotype, P[8], at least two distinct lineages were identified, both unrelated to P[8] Irish rotaviruses circulating in previous years, and more closely related to recent European humans rotaviruses. Moreover, sequence analysis of the VP7 of G1 rotaviruses revealed the onset of a G1 variant, previously unseen in the Irish population.

Detection and characterisation of group A rotavirus in asymptomatic piglets in southern Ireland

Porcine group A rotaviruses (GARV) are causative agents of enteritis in piglets and are a large reservoir of genetic material for the diversification of human GARVs. Accumulation of information on the genetic heterogeneity of porcine viruses is pivotal for readily characterising unusual human strains. Screening of 292 fecal samples, collected from 4-5- to 8-9-week-old asymptomatic pigs from four herds in Ireland between 2005 and 2007 resulted in 19 (6.5%) samples testing positive by reverse-transcription PCR (RT-PCR) for GARV. The strains were molecularly characterized to collate data on the VP7 and partial VP4 outer capsid genes. By sequence analysis of the VP7 gene, the Irish strains were identified as G2, G4, G5, G9 and G11 viruses. The G11 strains were closely related to other human and porcine G11 strains, while the G2 strains resembled porcine G2 viruses detected recently in Europe and southern Asia. The G4 strains were distantly related to other G4 human and animal strains, constituting a separate G4 VP7 lineage. Analysis of the G5 strains revealed that they were similar to a selection of G5 human and porcine strains, while the G9 strains resembled other porcine G9 viruses. By sequence analysis of the VP8* fragment of the VP4, the Irish viruses were characterised as P[6], P[7], P[13], P[13]/[22], P[26] and P[32].

Characterization of full-length VP4 genes of OP354-like P[8] human rotavirus strains detected in Bangladesh representing a novel P[8] subtype

The G1 and G9 rotavirus strains MMC71 and MMC38 (subgroup II, NSP4 genogroup B), respectively, isolated from children in Bangladesh, were analyzed genetically. Full-length VP4 genes of these strains had 98.9% identity to each other and showed 83.9-89.4% identity to those of the P[4] and P[8] rotaviruses. Phylogenetic analysis of VP4 nucleotide sequences revealed that strains MMC38 and MMC71 were located in a lineage of P[8] strains. However, the cluster was highly divergent from the previously established P[8] strains. The VP8* portions of strains MMC38 and MMC71 showed more than 93.9% nucleotide sequence identity to OP354-like P[8] strains, and these strains were clustered into the same lineage. These findings indicate that the VP4 of these strains should be classified into a subtype of the P[8] genotype (P[8]b) that is distinct from that of common P[8] rotaviruses (P[8]a).

Human group A rotavirus P[8] Hun9-like and rare OP354-like strains are circulating among diarrhoeic children in Eastern India

During 2004-2006, group A rotavirus P[8] strains were the major VP4 genotype (43.2%, n = 317) among diarrhoeic children in Eastern India. Phylogenetic analysis of VP8* amino acid sequences of 16 of these strains with other P[8] strains revealed four distinct lineages. P[8] strains from Eastern India clustered within rare OP354-like and Hun9-like lineages, pointing towards co-prevalence of divergent P[8] strains. Although it is unclear whether the observed genetic diversity might affect to some extent the efficacy of vaccines, the present study emphasized further efforts to address the much lacking information on diversity of P[8] strains across the Indian subcontinent.