Rotavirus genotyping is performed by using reverse transcription PCR with type-specific-primers. Because the high rotavirus mutation rate generates an extensive genomic variation, different G-type-specific primer sets are applied in different geographical locations. In Bangladesh, a significant proportion (36.9%) of the rotavirus strains isolated in 2002 could not be G-typed using the routinely used primer set. To investigate the reason why the strains were untypeable, nucleotide sequencing of the VP7 genes was performed.
Results
Four nucleotide substitutions at the G1 primer-binding site of the VP7 gene of Bangladeshi G1 rotaviruses rendered a major proportion of circulating strains untypeable using the routine primer set. Using an alternative primer set, we could identify G1 rotaviruses as the most prevalent genotype (44.8%), followed by G9 (21.7%), G2 (15.0%) and G4 (13.8%).
Conclusion
Because of the natural variation in the rotaviral gene sequences, close monitoring of rotavirus genotyping methods is important.
Background
Rotaviruses remain the most common cause of acute gastroenteritis worldwide and cause an estimated 600,000 deaths in children less than 5 years of age [20]. The high disease burden motivated major efforts to develop a suitable rotavirus vaccine. However, the vaccine efficacy is being challenged by the extensive strain diversity of the rotaviruses [3, 7–9, 13, 14].
Rotaviruses belong to the Reoviridae, and their genome consists of 11 segments of double stranded RNA. The gene segment coding for the VP7 glycoprotein is the basis for genotyping group A rotaviruses into at least 15 G-genotypes. Among them, G1, G2, G3, G4 and G9 are the most common G-types in humans [5, 15, 16, 19, 21, 23]. The importance of type-specific immunological protection against rotavirus disease is still under discussion [13].
G-genotyping is performed using type-specific-primer-based RT-PCR. Two common primer sets introduced by Gouvea et al. [6] and Das et al. [2] are currently used in rotavirus G-typing surveillance programs [22]. A failure to genotype or mistyping has already been reported from different parts of the world. These reports showed that nucleotide sequence differences between the target region of the respective genes and the primer sequences used for typing led to the genotyping failure [1, 10, 11, 17].
The Dhaka hospital of ICDDR,B, situated in the central Bangladesh, and the Matlab hospital, located 45 km south east of Dhaka respectively treat about 100,000 and 15,000 diarrhoeal patients each year. A hospital surveillance system has been established in these hospitals by ICDDR,B to collect information on clinical, epidemiological and demographic characteristics of the patients attending the hospital since 1978. In Bangladesh, rotavirus strains have previously been typed using a variety of techniques. Serotyping was introduced by Ward et al. [28] with specimens collected during 1985–1986 in Dhaka using neutralization with hyperimmune antisera against prototype rotavirus strains G1, G2, G3 and G4. They concluded that epitopic variations between rotavirus strains influenced the sensitivity of serotyping. Fun et al. [4] detected the major rotavirus types (G1 to G4) by RNA hybridization with serotype-specific synthetic oligonucleotide probes, but this method could not type 33.3% of the Bangladeshi rotaviruses. Likewise, RT-PCR depending on type-specific oligonucleotide primers failed to type a significant portion of rotaviruses. Rotavirus surveillance studies in Bangladesh between 1987 and 1997 reported that 1,095 (43.7%) samples out of 2,515 were G-untypeable [25–27].
In this study, we characterized rotavirus positive stool specimens collected in the Dhaka and Matlab hospitals during 2002 by using RT-PCR based on the primer set described by Das et al. [2]. We found that a major proportion of the specimens were untypeable. Nucleotide sequences of VP7 genes were performed to investigate the reason why they were untypeable with the routine primer set. The untypeable specimens were further characterized by using a different primer set described by Gouvea et al. [6].
Results and Discussion
Detection of rotavirus strains
In 2002, a total of 3,803 patients with history of watery diarrhea were included in the hospital surveillance system. In Dhaka and Matlab, 535 (27.2%) and 358 (19.4%) specimens were positive for group A rotavirus antigens by enzyme immunoassay.
G typing using the Das primer set
Rotavirus G-typing was carried out for all rotavirus positive specimens from Matlab and for every fourth of the rotavirus-positive specimens from Dhaka. Some samples were excluded from this study due to unavailability of sufficient amount of stool specimens for testing. G-typing was performed on 433 rotavirus ELISA-positive stool samples by RT-PCR using the primer set described by Das et al. [4], which was routinely used in our laboratory. The most prevalent G types were G9 (20.5%); G2 (14.6%), and G4 (13.8%). G1 comprised only 11.6% of the isolates and 36.9% of the rotavirus-positive samples were untypeable.
VP7 gene sequence analysis
We amplified the VP7 genes of five randomly selected untypeable strains (Dhaka162-02, Dhaka18-02, Dhaka164-02 Dhaka165-02 and Matlab26-02) using the VP7 consensus primers Beg9-End9 as described by Gouva et al. [6] and sequenced their complete open reading frame [GenBank:AY631050, GenBank:AY631054]. They were typed as G1 rotaviruses by using BLAST homology searches (99–100% nucleotide and amino acid identities with the Indian G1 rotavirus strain, ISO-4). To compare them with the typeable G1 sequences, the VP7 genes of two typeable G1 strains, Dhaka8-02 [GenBank:AY631049] and Matlab159-02 [GenBank:AY631055] were sequenced. We found that the nucleotide sequences of the typeable and untypeable G1 strains were 100% identical at the G1 primer-binding sites. We aligned the target G1 VP7 sequence with the Das G1 primer sequence (reverse primer, 9T1-1; 5'-TCTTGTCAAAGCAAATAATG-3'; nt 176–195, prototype strain Wa [GenBank:M21843]) to determine if there was any mismatch between them. Four mismatches were found in the Das G1 primer, 9T1-1, at the 5' end (Fig. 2). Due to these mismatches, the Das G1 primer failed to detect most (75%) of the G1 strains. Since, the primer set had perfect matches at the 3' end, it could detect 25% of the G1 rotaviruses. When we compared the target sequence with the Gouvea G1 primer sequence (forward primer, aBT1; 5'-CAAGTACTCAAATCAATGATGG-3'; nt 314–335, prototype strain Wa), we found only one mismatch (Fig. 2). Therefore, the Gouvea G1 primer was found to be more suitable for typing our G1 strains.
Figure 2
Nucleotide mismatches in the primers. The target sequence is the VP7 gene of G1 Bangladeshi strains. The G1 rotavirus VP7 gene specific primers were described by Das et al. [2] and Gouvea et al. [6]. Mismatches are in red.
Distribution of G types using Gouvea primer set
The untypeable specimens were typed using the primer set described by Gouvea et al. [6]. After typing with the Gouvea primer set, the distribution of rotavirus G-types changed dramatically (Fig. 1). Type G1 now accounted for 44.8% of the isolates and became the most prevalent genotype, and the number of untypeable strains was reduced from 36.9 to 2.1 %. The other common G types were G9 (21.7%), G2 (15.0%), and G4 (13.8%). The previous studies in Bangladesh reported that G4 strains were the most prevalent strains during 1992–1997 and a significant number of rotavirus strains were untypeable using the Das primer set [26]. It is likely that the Das primer set could not detect most of the G1 rotaviruses in the previous years and that a majority of the untypeable rotaviruses were G1 strains.
Figure 1
Effect of untypeable strains on G typing of Bangladeshi rotavirus strains isolated in 2002 (n = 433).
Conclusion
Because of the natural variation in the rotaviral gene sequences, close monitoring of rotavirus genotyping methods is important. The findings described in this paper will be important for genotyping strategies in the rotavirus surveillance studies.
Materials and methods
Sample collection
Stool specimens were collected from patients who presented with diarrhea to the Dhaka and Matlab hospitals of ICDDR,B in 2002. In the Dhaka hospital, stool specimens are routinely collected from every 50th patient and in Matlab hospital, every patient with diarrhea submits a stool specimen for testing.
Rotavirus antigen detection
Rotavirus antigens (group A-specific VP6 proteins) were detected in the stool specimen using a solid phase sandwich type enzyme immunoassay modelled after Dakopatts commercial kit incorporating rabbit hyperimmune antisera produced at ICDDR,B and an anti-human rotavirus-horseradish peroxidase conjugate (Dakopatts, Copenhagen, Denmark) using the same criteria for determination of positivity as those used by the Dakopatts kit [26].
RNA extraction
The QIAamp Viral RNA mini kit (Qiagen/Westburg, Leusden, The Netherlands) was used according to the manufacturer's instructions for the extraction of rotavirus RNA from the stool samples.
RT-PCR
A reverse transcriptase-polymerase chain reaction (RT-PCR) was carried out using the Qiagen OneStep RT-PCR Kit (Qiagen/Westburg) as previously described by Das et al. [2] and Gouvea et al. [6] for rotavirus G-types (G1, G2, G3, G4 and G9) using type-specific oligonucleotide primers. The reaction was carried out with an initial reverse transcription step at 45°C for 30 min, followed by 35 cycles of amplification (30 sec at 94°C, 30 sec at 50°C, 1 min at 72°C), and a final extension of 7 min at 72°C in a thermal cycler (Eppendorf, Hamburg, AG). PCR products were run on a 2% agarose gel, and stained with ethidium bromide. Specific segment sizes for the different G types were visualized under UV-light.
Nucleotide sequencing
The amplified PCR products were purified with the QIA quick PCR purification kit (Qiagen/Westburg), and sequenced in both directions using the dideoxy-nucleotide chain termination method with the ABI PRISM® BigDye Terminator Cycle Sequencing Reaction kit (Perkin-Elmer Applied Biosystems, Foster City, California) on an automated sequencer (ABI PRISM™ 3100). The Beg9 and End9 RT-PCR primers were used as sequencing primers.
Sequence analysis
The chromatogram sequencing files were inspected using Chromas 2.2 (Technelysium, Queensland, Australia), and consensus sequences were prepared using SeqMan II (DNASTAR, Madison, WI). Multiple sequence alignments were performed using CLUSTALX 1.81 [24]. Sequences were manually edited in the GeneDoc version 2.6.002 alignment editor [18].
Sequence submission
The nucleotide sequence data were deposited in GenBank using the National Center for Biotechnology Information (NCBI, Bethesda, MD) Sequin 5.15 submission tool http://www3.ncbi.nlm.nih.gov/ under accession numbers AY631049-AY631055.
Declarations
Acknowledgements
This study was funded by the Program for Appropriate Technology in Health (PATH), grant number GAT 770-790-01451-SPS. ICDDR,B acknowledges with gratitude the commitment of PATH to the Centre's research efforts.
Authors’ Affiliations
(1)
ICDDR,B: Centre for Health and Population Research
(2)
Laboratory of Clinical and Epidemiological Virology, Rega Institute for Medical Research, University of Leuven
References
Adah MI, Rohwedder A, Olaleyle OD, Werchau H: Nigerian rotavirus serotype G8 could not be typed by PCR due to nucleotide mutation at the 3' end of the primer binding site.Arch Virol 1997, 142: 1881-1887. 10.1007/s007050050206View ArticlePubMedGoogle Scholar
Das BK, Gentsch JR, Cicirello HG, Woods PA, Gupta A, Ramachandran M, Kumar R, Bhan MK, Glass RI: Characterization of rotavirus strains from newborns in New Delhi, India.J Clin Microbiol 1994, 32: 1820-1822.PubMed CentralPubMedGoogle Scholar
Estes MK: Rotaviruses and their replication. In Fields virology. Edited by: Howley PM. Philadelphia, Lippincott Williams & Wilkins; 2001:1747-1786.Google Scholar
Fun BN, Unicomb L, Rahim Z, Banu NN, Podder G, Clemens J, Van Loon FP, Rao MR, Malek A, Tzipori S: Rotavirus-associated diarrhea in rural Bangladesh: two-year study of incidence and serotype distribution.J Clin Microbiol 1991, 29: 1359-1363.PubMed CentralPubMedGoogle Scholar
Gentsch JR, Woods PA, Ramachandran M, Das BK, Leite JP, Alfieri A, Kumar R, Bhan MK, Glass RI: Review of G and P typing results from a global collection of rotavirus strains: implications for vaccine development.J Infect Dis 1996, 174: S30-S36.View ArticlePubMedGoogle Scholar
Gouvea V, Glass RI, Woods P, Taniguchi K, Clark HF, Forrester B: Polymerase chain reaction amplification and typing of rotavirus nucleic acid from stool specimens.J Clin Microbiol 1990, 28: 276-282.PubMed CentralPubMedGoogle Scholar
Green KY, Sears JF, Taniguchi K, Midthun K, Hoshino Y, Gorziglia M, Nishikawa K, Urasawa S, Kapikian AZ, Chanock RM: Prediction of human rotavirus serotype by nucleotide sequence analysis of the VP7 protein gene.J Virol 1988, 62: 1819-1823.PubMed CentralPubMedGoogle Scholar
Green KY, Midthun K, Gorziglia M, Hoshino Y, Kapikian AZ, Chanock RM, Flores J: Comparison of the amino acid sequences of the major neutralization protein of four human rotavirus serotypes.Virology 1987, 161: 153-159. 10.1016/0042-6822(87)90181-4View ArticlePubMedGoogle Scholar
Hoshino Y, Nishikawa K, Benfield DA, Gorziglia M: Mapping of antigenic sites involved in serotype-cross-reactive neutralization on group A rotavirus outercapsid glycoprotein VP7.Virology 1994, 199: 233-237. 10.1006/viro.1994.1117View ArticlePubMedGoogle Scholar
Iturriza-Gomara M, Kang G, Mammen A, Jana AK, Abraham M, Desselberger U, Brown D, Gray J: Characterization of G10P[11] rotaviruses causing acute gastroenteritis in neonates and infants in Vellore, India.J Clin Microbiol 2004, 42: 2541-2547. 10.1128/JCM.42.6.2541-2547.2004View ArticlePubMedGoogle Scholar
Iturriza-Gomara M, Green J, Brown DW, Desselberger U, Gray JJ: Diversity within the VP4 gene of rotavirus P[8] strains: implications for reverse transcription-PCR genotyping.J Clin Microbiol 2000, 38: 898-901.PubMed CentralPubMedGoogle Scholar
Jin Q, Ward RL, Knowlton DR, Gabbay YB, Linhares AC, Rappaport R, Woods PA, Glass RI, Gentsch JR: Divergence of VP7 genes of G1 rotaviruses isolated from infants vaccinated with reassortant rhesus rotaviruses.Arch Virol 1996, 141: 2057-2076.View ArticlePubMedGoogle Scholar
Kapikian AZ, Hoshino Y, Chanock RM: Rotaviruses. In Fields virology. Edited by: Howley PM. Philadelphia, Lippincott Williams & Wilkins; 2001:1787-1833.Google Scholar
Linhares AC, Lanata CF, Hausdorff WP, Gabbay WP, Black RE: Reappraisal of the Peruvian and Brazilian lower titer tetravalent rhesus-human reassortant rotavirus vaccine efficacy trials: analysis by severity of diarrhea.Pediatr Infect Dis J 1999, 18: 1001-1016. 10.1097/00006454-199911000-00013View ArticlePubMedGoogle Scholar
Liprandi F, Gerder M, Bastidas Z, Lopez JA, Pujol FH, Ludert JH, Joelsson DB, Ciarlet MA: Novel Type of VP4 Carried by a Porcine Rotavirus Strain.Virology 2003, 315: 373-380. 10.1016/S0042-6822(03)00534-8View ArticlePubMedGoogle Scholar
Martella V, Ciarlet M, Camarda A, Pratelli A, Tempesta M, Greco G, Cavalli A, Elia G, Decaro N, Terio V, Bozzo G, Camero M, Buonavoglia C: Molecular characterization of the VP4, VP6, VP7, and NSP4 genes of lapine rotaviruses identified in Italy: emergence of a novel VP4 genotype.Virology 2003, 314: 358-370. 10.1016/S0042-6822(03)00418-5View ArticlePubMedGoogle Scholar
Maunula L, von Bonsdorff C-H: Short sequences define genetic lineages: phylogenetic analysis of group A rotaviruses based on partial sequences of genome segments 4 and 9.J Gen Virol 1998, 79: 321-332.View ArticlePubMedGoogle Scholar
Nicholas KB, Nicholas HB, Deerfield DW: GeneDoc: analysis and visualization of genetic variation.Embnet News 1997, 4: 14.Google Scholar
Okada J, Urasawa T, Kobayashi N, Taniguchi K, Hasegawa A, Mise K, Urasawa S: New P serotype of group A human rotavirus closely related to that of a porcine rotavirus.J Med Virol 2000, 60: 63-69. 10.1002/(SICI)1096-9071(200001)60:1<63::AID-JMV11>3.0.CO;2-4View ArticlePubMedGoogle Scholar
Parashar UD, Hummelman EG, Bresee JS, Miller MA, Glass RI: Global illness and deaths caused by rotavirus disease in children.Emerg Infect Dis 2003, 9: 565-572.PubMed CentralView ArticlePubMedGoogle Scholar
Rao CD, Gowda K, Reddy BS: Sequence analysis of VP4 and VP7 genes of nontypeable strains identifies a new pair of outer capsid proteins representing novel P and G genotypes in bovine rotaviruses.Virology 2000, 276: 104-113. 10.1006/viro.2000.0472View ArticlePubMedGoogle Scholar
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. 10.1002/rmv.448View ArticlePubMedGoogle Scholar
Sereno MM, Gorziglia MI: The outer capsid protein VP4 of murine rotavirus strain Eb represents a tentative new P type.Virology 1994, 199: 500-504. 10.1006/viro.1994.1153View ArticlePubMedGoogle Scholar
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools.Nucleic Acids Res 1997, 25: 4876-4882. 10.1093/nar/25.24.4876PubMed CentralView ArticlePubMedGoogle Scholar
Unicomb LE, Bingnan F, Rahim Z, Banu NN, Gomes JG, Podder G, Munshi MH, Tzipori SR: A one-year survey of rotavirus strains from three locations in Bangladesh.Arch Virol 1993, 132: 201-208. 10.1007/BF01309854View ArticlePubMedGoogle Scholar
Unicomb LE, Podder G, Gentsch JR, Woods PA, Hasan KZ, Faruque AS, Albert MJ, Glass RI: Evidence of high-frequency genomic reassortment of group A rotavirus strains in Bangladesh: emergence of type G9 in 1995.J Clin Microbiol 1999, 37: 1885-1891.PubMed CentralPubMedGoogle Scholar
Unicomb LE, Kilgore PE, Faruque ASG, Hamadani JD, Fuchs GJ, Albert MJ, Glass RI: Anticipating rotavirus vaccines: hospital-based surveillance for rotavirus diarrhea and estimates of disease burden in Bangladesh.Pediatr Infect Dis J 1997, 16: 947-951. 10.1097/00006454-199710000-00008View ArticlePubMedGoogle Scholar
Ward RL, McNeal MM, Clemens JD, Sack DA, Rao M, Huda N, Green KY, Kapikian AZ, Coulson BS, Bishop RF: Reactivities of serotyping monoclonal antibodies with culture-adapted human rotaviruses.J Clin Microbiol 1991, 29: 449-456.PubMed CentralPubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
