Identification of OP354-like human rotavirus strains with subtype Pb in Ghanaian children with diarrhoea
Virology Journal volume 13, Article number: 69 (2016)
Rotaviruses with the P genotype have been associated with majority of infections. Recent improvements in molecular diagnostics have delineated the P genotype into Pa and Pb subtypes. Pa is the previously known P genotype which is common whilst Pb 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.
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.
Of the 80 RVA-positive samples, 57 were successfully sequenced and characterized. Forty-eight of these were identified as P strains of which 5 were characterized as the rare Pb subtype. Phylogenetic analysis of the VP8* fragment of the VP4 genes of these Pb strains revealed a close relationship with prototype OP354-like Pb strain and Pb strains of Russian and South African Pb origin.
The study highlights the importance of regularly updating the primers employed for molecular typing of rotaviruses.
Group A rotaviruses are the most important aetiological agent of severe dehydrating diarrhoea in children less than 5 years of age accounting for an estimated 453,000 deaths annually . Rotaviruses with the P genotype have been associated with majority of these infections, accounting for approximately 74 % of global prevalence . In order to control severe disease caused by rotavirus infection, vaccination is considered an essential strategy. The two currently available rotavirus vaccines (Rotarix™ and RotaTeq™) both include the rotavirus P genotype, which has recently been classified into two subtypes- Pa and Pb (OP354-like) [3, 4]. Whilst Pa which corresponds to the previously known P genotype is common, Pb is rarely detected, genetically distinct from Pa with limited reports globally [4–6]. The prototype OP354 strain bearing Pb specificity was detected in the Malawian strain . A recently inferred evolutionary history of OP354-like Pb rotavirus strains globally showed that these strains emerged relatively recently, spreading from South and East Asia to Europe, Sub-Saharan Africa and North America . In previous WHO sponsored rotavirus surveillance studies in selected African countries,  the genotypes of 80 RVA-positive samples with first round RT-PCR amplification products could not be determined with the recommended primers presently available within the WHO Regional Reference Laboratory in Ghana. This study sought to determine genotypes of previously non-typeable rotavirus VP4 genes using sequence-dependent cDNA amplification methods.
RVA-positive samples (n = 80) whose VP4 genes could not be characterised using the pool of genotype specific primers recommended for use in the WHO Regional Reference Laboratory in Ghana were included in the study. Fifty-seven of the RVA-positive samples were successfully sequenced and characterised using the web-based genotyping tool RotaC version 2.2 . Forty-eight (84.2 %) of these were characterised as P strains of which 43 (89.6 %) were identified as Pa and 5 (10.4 %) as the rare OP354-like Pb subtype. The 5 Pb subtypes detected in the study were found in combination with G9 genotype. Sequence mismatches were detected in the primer binding region of the target gene of the 43 Ghanaian Pa rotavirus strains when aligned with the available WHO Pa primers (1 T-1, 1 T-1Wa and 1 T1-VN) which might have resulted in genotyping failure (data not shown). The Ghanaian Pb rotavirus strains (GHA-949/DC/2010, GHA-716/PML/2010, GHA-094/M/2010, GHA-176/M/2010 and GHA-192/M/2010) shared more than 99 % nucleotide sequence identity amongst themselves and with the Russian Pb strain (Rus/Nov06-1486). However, they shared a slightly lower nucleotide identity (96.3–98.9 %) with other published African Pb strains (Table 1). The Ghanaian strains formed a monophyletic cluster within lineage I of the Pb subtype, clustering along with the prototype strain OP354 and other published strains from Africa (Ethiopia, Malawi, Togo South Africa), Asia (India, Bangladesh, Pakistan) and Eastern-Europe (Russia) (Fig. 1).
Rotavirus continues to be a major cause of diarrhoeal disease among children less than 5 years old. Among the five major rotavirus strains (G1P, G2P, G3P, G4P and G9P) causing severe disease in children, G1P has been found to be the predominant strain [2, 10]. Rotaviruses with the P VP4 genotype are a major cause of acute infantile diarrhoea. Recent improvement in molecular characterisation has delineated the P genotype into Pa and Pb subtypes. The Pa subtype is the previously known common Pgenotype, whilst the Pb (OP354-like P) subtype is rare and has been identified in only a few cases . In this study, we report for the first time rotaviruses with the Pb subtype as a cause of diarrhoeal disease in hospitalized Ghanaian children. The prototype Pb rotavirus (strain OP354) was reported in Malawi, Blantyre during a two year study of viral gastroenteritis in children (1997–1999)  and subsequently detected mainly in Asian countries [11–13]. The Ghanaian OP354-like P strains were phylogenetically closely related to strains mostly from Africa and South Asia suggesting common origin. This finding supports the recently inferred evolutionary history of OP354-like Pb rotavirus strains globally . Compared with other VP4 P lineages, strains bearing OP354-like P genes emerged relatively recently, spreading from South and East Asia to Europe, Sub-Saharan Africa and North America in less than 20 years. Considering that all the Ghanaian OP354-like P strains were isolated in the same rotavirus season, it is unclear whether they are continually circulating or their detection resulted from a one-time seeding event. Continuous surveillance of rotavirus strains using updated PCR primer sets capable of detecting OP354-like P strains should clarify this, and in the process, determine whether they will be an important genotype in the human rotavirus population. Though, Nguyen et al.,  speculated that Pb rotaviruses cause more severe disease than Pa subtypes, this study was unable to validate this position or otherwise. This was due to the disproportionately small number of Pb strains as well as insufficient clinical records. In the present study, 10.4 % (5/48) of previously non-typeable P rotaviruses were successfully sequenced and characterized as the rare OP354-like Pb rotavirus strain. Both of the currently available rotavirus vaccines [Rotarix™ and RotaTeq™] include Pa VP4 gene but not Pb [3, 14]. It has been reported that there is a relatively large genetic distance between OP354-like P strains and the Pa strain contained in both vaccines which is translated into multiple differences in antigenic epitopes . Therefore, it remains unclear whether these vaccines will provide efficient protection against Pb rotavirus strains [4, 15]. Full genome characterization of the Ghanaian Pb subtype is presently being conducted to elucidate their origin and evolutionary dynamics.
We have been able to characterize the VP4 genes of previously non-typeable Ghanaian rotavirus genotypes and established their phylogeny. The rare OP354-like Pb subtype of rotavirus strain has been detected and reported for the first time in Ghana. Expansion of Pb subtype data would lead to a better understanding of transmission dynamics of this rare subtype in Ghana in relation to Pa. The study highlights the importance of regularly updating PCR primers employed for molecular typing of rotaviruses. It also indicates the necessity for continuous surveillance especially in the post-vaccine era.
The study was reviewed and approved by the Institutional Review Board, Noguchi Memorial Institute for Medical Research, Legon, Accra, Ghana.
Viral RNA was extracted from 10 % stool suspensions by the phenol-chloroform method  and purified with the RNaid® Kit (Bio 101, Carlsbad, USA). A one step RT-PCR targeting the VP4 genes of RVA was carried out using the PCR primer pairs, Con2/Con3 and VP4-F/VP4-R to generate 876 bp and 663 bp VP8*gene fragments respectively. The success of gene amplification was verified by electrophoresis on 2 % agarose gel stained with ethidium bromide. The VP4-F/VP4-R RT-PCR products (663 bp) were purified with the QIAquick PCR purification kit (Qiagen/Westburg) following the manufacturer’s instructions. Sequencing of the purified products was performed using the di-deoxynucleotide chain termination method with the ABI PRISM™ BigDye Terminator Cycle Sequencing Reaction Kit v.3.1 (Applied Biosystems, Foster City, CA, USA). After post purification, sequences were read in an automated sequencer (ABI PRISM™ 3130).
Consensus sequences determined from forward and reverse sequences were used to query the nucleotide database in GenBank using the Basic Local Alignment Search Tool [BLAST] (http://blast.ncbi.nlm.nih.gov/Blast.cgi). Genotypes were confirmed using the automated genotyping tool, RotaC v2.0 . Sequences were edited, trimmed to the same length to maximise homology and homology table was constructed using Bioedit v. 7. 0. 5 (Table 1). Nucleotide sequences of confirmed rotavirus genotypes were aligned with cognate gene sequences available in GenBank using the ClustalW algorithm , and phylogenetic analysis was performed using the Neighbour-Joining method in MEGA v6.06 software [18, 19]. Phylogenetic tree was statistically supported by bootstrapping with 1000 replicates, and phylogenetic distances calculated using the Kimura-2 parameter model. The Ghanaian Pb and reference Pb strains used in the construction of the phylogenetic tree are listed in Table 2.
Nucleotide sequence accession numbers
The nucleotide sequences of the Pb rotaviruses reported in the study have been deposited in the GenBank database under accession numbers KM379145-KM379149 (GHA-949/DC/2010; GHA-716/PML/2010; GHA-094/M/2010; GHA-176/M/2010 and GHA-192/M/2010 respectively).
Tate JE, Burton AH, Boschi-Pinto C, Steele AD, Duque J, Parashar UD. The 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(2):136–41.
Estes M, Kapikian A. Rotaviruses. In: Fields BN, Knipe DM, Howley PM, Griffin DE, Lamb RA, Martin MA, Roizman B, Straus ES, editors. Fields virology, vol. 2. 5th ed. Philadelphia: Lippincott Williams & Wilkins; 2007.
Greenberg HB, Estes MK. Rotaviruses: from pathogenesis to vaccination. Gastroenterology. 2009;136:1939–51.
Nagashima S, Kobayashi N, Paul SK, Alam MM, Chawla-Sakar M, Krishnan T. Characterization of full-length VP4 genes of OP354-like P human rotavirus strains detected in Bangladesh representing a novel P subtype. Arch Virol. 2009;54:1223–31.
Aung MS, Aung TS, Oo KY, Ne W, Gosh Syamamoto D, Kobayashi N. Phylogenetic analysis of human rotavirus in Myanmar: detection of Indian-Bangladeshi G1/G2lineages, Chinese G3 lineages and OP354-like P lineage (Pb subtype). South-east Asian J Trop Med Public Health. 2010;41:1393–404.
Cunliffe NA, Gondwe JS, Graham SM, Thindwa BD, Dove W, Broadhead RL, Molyneux ME, Hart CA. Rotavirus strain diversity in Blantyre, Malawi from 1997 to 1999. J Clin Microbiol. 2001;39:836–43.
Zeller M, Heylen E, Damanka S, Pictsch C, Donato C, et al. Emerging OP354-like P rotaviruses have rapidly dispersed from Asia to other continents. Mol Biol Evol. 2015;32:2060–71.
Enweronu-Laryea CC, Sagoe KW, Damanka S, Lartey B, Armah EA. Rotavirus genotypes associated with childhood severe acute diarrhoea in Southern Ghana: a cross-sectional study. Virol J. 2013;10:287.
Maes P, Matthijnssens J, Rahman M, Van Ranst M. RotaC: a web-based tool for the complete genome classification of group A rotaviruses. BMC Microbiol. 2009;9:238.
World Health Organisation. Global rotavirus information and surveillance bulletin. Reporting period: January through December 2010, vol. 4, World Health Organization; 2011.
Nguyen TA, Hoang LP, Pham LD, Hoang KT, Okitsu S, Mizuguchi M, Ushijima H. Use of sequence analysis of the VP4 gene to classify recent Vietnamese rotavirus isolates. Clin Microbiol Infect. 2008;14:235–41.
Samajdar S, Ghosh S, Dutta D, Chawla-Sarkar M, Kobayashi N, Naik TN. Human group A rotavirus P Hun9-like and rare OP354-like strains are circulating among diarrhoeic chidren in Eastern India. Arch Virol. 2008;153:1933–6.
Tamura T, Nishikawa M, Anh DD, Suzuki H. Molecular Epidemiological study of rotavirus and norovirus infection among children with acute gastroenteritis in Nha Trang, Vietnam, December 2005-June 2006. Jpn J Infect Dis. 2010;63:405–11.
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.
Nagashima S, Kobayashi N, Paul SK, Ghosh S, Chawla-Sarkar M, Hossain MA, Krishnan T. Identification of Pb subtype in OP354-like human rotavirus strains by a modified RT-PCR method. Jpn J Infect Dis. 2010;63:208–11.
Steele AD, Alexander JJ. Molecular epidemiology of rotaviruses in black infants in South Africa. J Clin Microbiol. 1987;23:992–4.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.
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.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.
The authors acknowledge the staff of the Department of Electron Microscopy and Histopathology, Noguchi Memorial Institute for Medical Research, University of Ghana for their technical assistance.
“The authors declare that there is no conflict of interest regarding the publication of this paper”.
SD and GEA designed the study. SD, FED, BL and CA participated in the genetic analysis of the data and SD coordinated with GEA, FED and BL in drafting the manuscript. Authors of the manuscript read and approved of the final copy.
SD (PhD) is Chief Research Assistant, FED (PhD) is a Research Fellow, CA (PhD candidate) is Principal Research Assistant and BL (PhD candidate) is Principal Research Assistant in the Department of Electron Microscopy and Histopathology (Regional Rotavirus Reference Laboratory). TA (PhD) is the Head of Department of Microbiology, University of Ghana Medical School and KWS (PhD) is the Head of Medical Virology Laboratory, University of Ghana Medical School. KN (PhD) is a lecturer in the Department of Epidemiology and Disease Control, School of Public Health. CEL (MRCPCH) and RO (FRCPCH) are paediatricians in the Department of Child Health, Korle Bu Teaching Hospital. MO (PhD) is the Head of Department of Electron Microscopy and Histopathology and GEA (PhD) is the Head of Regional Rotavirus Reference Laboratory.
About this article
Cite this article
Damanka, S., Dennis, F.E., Agbemabiese, C. et al. Identification of OP354-like human rotavirus strains with subtype Pb in Ghanaian children with diarrhoea. Virol J 13, 69 (2016). https://doi.org/10.1186/s12985-016-0523-5