- Research note
- Open Access
Molecular characteristics of rotavirus genotypes circulating in the south of Benin, 2016–2018
BMC Research Notes volume 13, Article number: 485 (2020)
Rotavirus remains the main causative agent of gastroenteritis in young children in countries that have not yet introduced the vaccine. In Benin, rotavirus vaccine was introduced late December 2019 into the EPI. This study aims to provide pre-vaccination era rotavirus genotyping data in Benin. These data can supplement data from the surveillance system of Ministry of Health of Benin which is supported by the World Health Organization (WHO).
Of the 420 diarrheal stool samples, actively collected in southern Benin from July 2016 through November 2018 from children under 5 years old and suffering from gastroenteritis, 167 (39.8%) samples were rotavirus EIA positive. 186 (44.3%) samples contained amplifiable rotavirus RNA detected by qRT-PCR method and were genotyped using one-step RT-PCR multiplex genotyping method. G1P represents the predominant genotype (32%) followed by the G2P (26%), G3P (16%), G12P (13%) and mixed G and P types (1%). Four samples (2%) could not be assigned both G and P type specificity.
Diarrhea ranks as the fifth leading cause of mortality among children under 5 years old . Infectious diarrhea can be caused by microbial agents of bacterial, parasitic, viral or mycotic origin . Viruses, especially rotavirus, are predominant as causative agents with more than 80% of cases [3, 4]. Globally, deaths associated with rotavirus in 2013 were estimated at 122,000–215,000 in children under 5 years of age and the largest number of rotavirus deaths occurred in sub-Saharan Africa .
Rotaviruses belong to the Reoviridae family, and possess a triple-layered icosahedral capsid enclosing a genome of 11 segments of double-stranded RNA encoding six structural and five or six non-structural proteins [6, 7]. The outer capsid proteins VP7 and VP4 define the G and P genotypes, respectively [8, 9] and rotaviruses are currently classified into 36 G and 51 P genotypes [RCWG, https://rega.kuleuven.be/cev/viralmetagenomics/virus-classification]. Although a large number of G/P genotype combinations have been reported [10, 11], the rotavirus genotypes G1P, G2P, G3P, G4P and G9P together comprise up to three quarters of human rotavirus infections worldwide [12, 13].
Rotarix and RotaTeq, have been approved and licensed in many countries. Recently, two vaccines, ROTAVAC® and ROTASIL® have been prequalified by World Health Organization (WHO) . By the end of 2017, 32 (68%) of 47 countries in the WHO's African Region had introduced rotavirus vaccine into their national immunization programs (NIP) .
In Benin, rotavirus vaccine was introduced into the NIP late December 2019. It is therefore necessary to continue rotavirus strain surveillance in order to monitor rotavirus genotype post-vaccination era. This study aims to provide data for the pre-vaccine introduction period on rotavirus genotypes circulating in Benin. These data can supplement data from the surveillance system of the Ministry of Health of Benin which is supported by the WHO.
Study design, specimen handling and transport
This study is a descriptive investigation of the molecular epidemiology of rotavirus from July 2016 through November 2018 from children under 5 years with diarrhea. Diarrhea stool samples were collected from the Mênontin, Anastasie, and St. Vincent Health Centres in Cotonou and from the Central Polyclinic of Abomey-Calavi (Additional file 1: Figure S1a). Children with bloody diarrhea and those older than 5 years were excluded from this study.
All stool samples were stored at 4 °C in the laboratory or − 20 °C to facilitate further testing. All stool samples were shipped on dry ice to the CDC, USA for genotyping analysis.
Stool processing and nucleic acid extraction
A 10% stool suspension was prepared for each sample using phosphate-buffered saline and RNA was extracted from the suspension using either the MagMAX™-96 Viral RNA Isolation Kit (ThermoFisher Scientific, Vilnius, Lithuania) on the KingFisher™ Flex Purification System (ThermoFisher Scientific, Vantaa, Finland) or the MagNA Pure Compact RNA extraction kit on the MagNA Pure Compact instrument (Roche Applied Science, Indianapolis, IN, USA) following the manufacturer’s instructions. Prior to each of the above extraction procedures, 2 µL of 109 unit/µL of MS2 bacteriophage RNA (ZeptoMetrix, Buffalo, NY, USA), were spiked into a 48 µL or 98 µL volume of 10% stool suspension to serve as an internal process control.
VP7 and VP4 genotyping
Genotyping was performed using reverse transcription-polymerase chain reaction (RT-PCR) to determine the G and P-genotypes. VP7 and VP4 genotyping RT-PCR was performed using a conventional multiplexed one-step amplification process with slight modifications . In brief, the genotype G2 (G2-R4), G4 (G4-R2) and G9 (G9-R2) specific primers were replaced with an updated versions G2 (G2-R1: TAT GTA GTC CAT YGT ATT AGT), G4 (G4-R1: GAG CAT TCG MTA ATA MTG ATA ATA C), and G9 (G9-R3: CAG AGT ATY YTT CCA TTC HGT ATC TCC) primers. The VP7 and VP4 conventional multiplexed one-step RT-PCR genotyping product was electrophoresed on 3% agarose gels containing GelRed (Biotium, Heyward, CA, USA) for 2–3 h at 100 V and products were detected under UV transillumination or were analyzed on the LabChip®GX instrument (Caliper, Life Sciences, MA, USA) using the HT DNA 1 K or 5 K reagent kit (Dual protocol DNA Analysis and Quantitation) with the HT Extended Range LabChip (Caliper Life Sciences, MA, USA) as described previously .
Demographic characteristics and genotyping results
Distribution of the study population by age and sex
In this study, children aged from 7 to 12 months were the most represented with a percentage of 41.9% followed by those under 6 months of age (Additional file 1: Table S1). The age average was 11.3 months and the sex ratio M/F was 1.12.
All the samples were tested by EIA and the EIA positivity rate (prevalence) was 39.8% (167/420).
Distribution of genotypes in Benin
Out of the 420 samples, 186 (44.3%) samples contained rotavirus RNA detected by RT-PCR methods. Four samples (2%) could not be assigned both G and P type specificity (non-typeable samples). G1P represents the predominant genotype (32%) circulating in southern Benin, followed by the G2P (26%), G3P (16%) and G12P (13%). Other uncommon genotypes in Benin such as G9P (3%), G1P (3%), G3P (2%), G2P (1%), G2P (1%), G12P (1%), G9P (1%) and mixed genotypes (1%) were detected, albeit at low rate (Fig. 1a).
Distribution of genotypes by gender and age in Benin: There was greatest variability of genotypes in children aged from 0 to 18 months. This variability was less after 18 months. Uncommon genotypes were most observed in children under 24 months of age (Fig. 1b). However, there were no significant differences in genotype distribution when genders were compared (p > 0.05).
Distribution of genotypes in Cotonou
From the 330 samples collected, 147 (44.6%) contained detectable rotavirus RNA and 144 were genotyped for G and P. The common genotype combinations were G1P (30%) followed by G2P (23%), G3P (16%), and G12P (15%) (Fig. 1c). Three (2%) stool samples were non-typeable for G and P and only one (1%) sample showed mixed genotype. Less common genotypes such as G1P (4%), G9P (3%), G3P (2%), G2P (1%) G2P (1%), G12P (1%) and G9P (1%) were detected.
There is great variability of genotypes in children aged from 0 to 18 months. This variability is less after 18 months (Additional file 1: Figure S1b).
Distribution of genotypes in Abomey-Calavi
In Abomey-Calavi, 39 (43.3%) stool samples were positive by RT-PCR. G1P and G2P genotype combinations were the most common (38% and 36% respectively), followed by G3P (15%), G12P (5%). No mixed genotypes were found and one sample was untypeable (Additional file 1: Figure S1c).
The greatest variability of genotypes was more observed in children less than 18 months in Abomey-Calavi. (Additional file 1: Figure S1d).
In the two localities that hosted the study, G1P was the prevalent genotype, accounting for 38% in Abomey-Calavi and 30.0% in Cotonou or 32.0% across the two sites. According to a systematic review  describing the epidemiological situation in Africa, G1P similarly was the predominant genotype with an overall prevalence of 22.64%. Taking into account the five regions of Africa (West, East, Central, North and South), the rate obtained varied from one region to another with a higher prevalence in North Africa (37.10%), while it was 14.35% in the West Africa .
At country level in West Africa, before vaccine introduction, Ghana reported G1P (20%) as the predominant genotype , but others genotypes such as G6P was found to be medically important in Burkina Faso , G12P in Côte d’Ivoire  and G4P and G12P in Nigeria [21, 22]. During the post vaccine period, Ghana reported G1P as the fourth common genotype detected with a rate of 8.0% . Thus, there was a decrease in the prevalence rate and a replacement of strains after the introduction of the vaccine. Nearly similar observations are made in Zambia where G1P initially reported as the predominant genotype (49%) in 2008 during the pre-vaccine period remained predominant the first year after vaccine introduction at 25% but changed in subsequent years . On the other hand, according to the study of Seheri et al. , in Eastern and Southern African countries before and after vaccine introduction, data from six countries showed no difference in strains circulation during the pre- and post-vaccine introduction eras. So, it is not conclusive that after vaccine introduction there is a switch of strains due to vaccine pressure or this is just an annual fluctuation of strains.
In systematic review by Ouermi et al.  showed that G1P was the most predominant genotype, followed by G2P, G9P and G2P. Similarly, in our study, G2P proved to be the second-most prevalent genotype circulating in the south of the country where the study was conducted. However, in Niger and the Democratic Republic of Congo, G2P was the predominant genotype before vaccine introduction [25, 26]. However, G2P was found circulating in the Central Africa Republic at a rate of 13%, in Cameroon at 5.9%, Burkina Faso at < 1% and Gabon at < 2% [27,28,29,30].
G12P is the second emerging genotype, detected for the first time in Philippines in 1987, in children under 2 years of age . Mainly limited to Asian countries, G12P has been reported worldwide over the last 20 years .
Commonly detected at low rates in Africa, over the past decade, G12P has been isolated as a predominant genotype in some countries such as Côte d'Ivoire and Nigeria in the pre-vaccination period [20, 22], and Ghana in the post-vaccination period . Also, according to some studies, G12P may come from a reassortment between human and pig rotavirus strains . The presence of this genotype in our study without any pre-existing vaccination context, strongly suggests that G12P appearance in Benin, is due to natural rotavirus genotype occurrence, not vaccine pressure.
On the other hand, G9P recognized as being the fifth-most prevailing genotype [10, 12] was found in our study at a rate of 3%. This genotype seems to be more present in East and North Africa [33, 34]. Even at this relatively low level, it deserves to be monitored in Benin since it was one of the emerging genotypes of the last 20 years .
During the study 2% (4/186) of the samples were untypeable for G or P, which is less than the rate of untypeable samples previously reported from Sub-Saharan countries (8.6–14.6%) . This low level of untypeable G/P genotypes could be explained by the technique used which incorporates more recently designed primers.
Unlike reported by the African Rotavirus Surveillance Network where mixed genotypes accounted for 12–14% , in our study, a single sample (1%) showed a mixed genotype. Indeed, close association of humans with domesticated animals in most countries of the developing region lead to gene reassortment events within commonly circulating animal rotavirus strains and thus give rise to a large genomic diversity and frequent occurrence of mixed infections . However, studies from some Africans settings have shown no mixed genotypes  or a very low rate (1%) of mixed infections  similar to our findings. Also, as suggested by comments from the works of Boula and colleagues, the incorporation and use of new and updated one step multiplex genotyping assays in this study could strongly explain the low rate of mixed genotypes.
Atypical rotavirus G and P combinations were detected at low frequency (2%) similar to the rate described previously by Seheri et al. . No correlation between age group and genotype was observed in this study, although it has been noted a higher variability of genotypes in children less than 18 months. This observation could probably due to the fact that more samples came from the youngest age groups.
It is important that further studies be conducted to determine the true intrinsic determinants of occurrence of rotavirus gastroenteritis in children. Studies on the genetic differences in histo-blood group antigens in our populations (secretor, non-secretor, Lewis antigens) could help to learn more about rotavirus epidemiology.
The limitations of this study are due to the fact that samples came from the southern part of Benin where the distribution of genotypes could differ from that of the northern region due to the difference in seasons and the living behavior of the populations between the two regions.
Availability of data and materials
The raw genotyping data including those in the supplementary table and figures are kept by the Centers for Disease Control and Prevention (CDC) and can be requested from Dr. Mathew Esona using email address: firstname.lastname@example.org.
Enzyme linked immunosorbent assay
National Immunization Programs
Quantitative reverse transcription-polymerase chain reaction
Reverse transcription-polymerase chain reaction
World Health Organization
GDD Collaborators. Estimates of the global, regional, and national morbidity, mortality, and aetiologies of diarrhoea in 195 countries: a systematic analysis for the Global Burden of Disease Study 2016. Lancet Infect Dis. 2018;11:1211–28.
Liu L, Johnson HL, Cousens S, et al. Global, regional, and national causes of child mortality: an updated systematic analysis for 2010 with time trends since 2000. Lancet. 2012;379:2151–61.
Tate JE, Burton AH, Boschi-Pinto C, et al. 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.
WHO. Rotavirus vaccines. https://www.who.int/wer. 2013;5:49–64.
Tate JE, Burton AH, Boschi-Pinto C, Parashar UD, World Health Organization-Coordinated Global Rotavirus Surveillance N. Global, regional, and national estimates of rotavirus mortality in children < 5 years of age, 2000–2013. Clin Infect Dis. 2016;62(Suppl 2):S96–105.
Desselberger U. Rotaviruses. Virus Res. 2014;190:75–96.
Mihalov-Kovacs E, Gellert A, Marton S, et al. Candidate new rotavirus species in sheltered dogs, Hungary. Emerg Infect Dis. 2015;21:660–3.
Gentsch JR, Glass RI, Woods P, et al. Identification of group A rotavirus gene 4 types by polymerase chain reaction. J Clin Microbiol. 1992;30:1365–73.
Gouvea V, Glass RI, Woods P, et al. Polymerase chain reaction amplification and typing of rotavirus nucleic acid from stool specimens. J Clin Microbiol. 1990;28:276–82.
Iturriza-Gomara M, Dallman T, Banyai K, et al. Rotavirus genotypes co-circulating in Europe between 2006 and 2009 as determined by EuroRotaNet, a pan-European collaborative strain surveillance network. Epidemiol Infect. 2011;139:895–909.
Matthijnssens J, Van Ranst M. Genotype constellation and evolution of group A rotaviruses infecting humans. Curr Opin Virol. 2012;2:426–33.
Banyai K, Laszlo B, Duque J, et al. 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;30(Suppl 1):A122-130.
Esona MD, Gautam R. Rotavirus. Clin Lab Med. 2015;35:363–91.
Kirkwood CD, Steele AD. Rotavirus vaccines in China: improvement still required. JAMA Netw Open. 2018;1:e181579.
Mwenda JM, Parashar UD, Cohen AL, Tate JE. Impact of rotavirus vaccines in Sub-Saharan African countries. Vaccine. 2018;36:7119–23.
Esona MD, Gautam R, Tam KI, Williams A, Mijatovic-Rustempasic S, Bowen MD. Multiplexed one-step RT-PCR VP7 and VP4 genotyping assays for rotaviruses using updated primers. J Virol Methods. 2015;223:96–104.
Ouermi D, Soubeiga D, Nadembega WMC, et al. Molecular epidemiology of rotavirus in children under five in Africa (2006–2016): a systematic review. Pak J Biol Sci. 2017;20:59–69.
Lartey BL, Damanka S, Dennis FE, et al. Rotavirus strain distribution in Ghana pre- and post-rotavirus vaccine introduction. Vaccine. 2018;36:7238–42.
Nordgren J, Bonkoungou IJ, Nitiema LW, et al. Rotavirus in diarrheal children in rural Burkina Faso: high prevalence of genotype G6P. Infect Genet Evol. 2012;12:1892–8.
Boni-Cisse C, Meite S, Mlan AB, et al. Genotypic characterization of rotavirus in children under 5 years circulating in Cote D’Ivoire from 2010 to 2013. Virol J. 2018;15:78.
Alkali BR, Daneji AI, Magaji AA, Bilbis LS, Bande F. Molecular characterization of human rotavirus from children with diarrhoeal disease in Sokoto State, Nigeria. Mol Biol Int. 2016;2016:1876065.
Japhet MO, Famurewa O, Iturriza-Gomara M, et al. Group A rotaviruses circulating prior to a national immunization programme in Nigeria: cinical manifestations, high G12P frequency, intra-genotypic divergence of VP4 and VP7. J Med Virol. 2018;90:239–49.
Simwaka JC, Mpabalwani EM, Seheri M, et al. Diversity of rotavirus strains circulating in children under five years of age who presented with acute gastroenteritis before and after rotavirus vaccine introduction, University Teaching Hospital, Lusaka, Zambia, 2008–2015. Vaccine. 2018;36:7243–7.
Seheri LM, Magagula NB, Peenze I, et al. Rotavirus strain diversity in Eastern and Southern African countries before and after vaccine introduction. Vaccine. 2018;36:7222–30.
Page AL, Jusot V, Mamaty AA, et al. Rotavirus surveillance in urban and rural areas of Niger, April 2010–March 2012. Emerg Infect Dis. 2014;20:573–80.
Pukuta ES, Esona MD, Nkongolo A, et al. Molecular surveillance of rotavirus infection in the Democratic Republic of the Congo August 2009 to June 2012. Pediatr Infect Dis J. 2014;33:355–9.
Banga-Mingo V, Waku-Kouomou D, Gody JC, et al. Molecular surveillance of rotavirus infection in Bangui, Central African Republic, October 2011–September 2013. Infect Genet Evol. 2014;28:476–9.
Bonkoungou IJO, Ouedraogo N, Tamini L, et al. Rotavirus and norovirus in children with severe diarrhea in Burkina Faso before rotavirus vaccine introduction. J Med Virol. 2018;90:1453–60.
Lekana-Douki SE, Kombila-Koumavor C, Nkoghe D, Drosten C, Drexler JF, Leroy EM. Molecular epidemiology of enteric viruses and genotyping of rotavirus A, adenovirus and astrovirus among children under 5 years old in Gabon. Int J Infect Dis. 2015;34:90–5.
Ndze VN, Papp H, Achidi EA, et al. One year survey of human rotavirus strains suggests the emergence of genotype G12 in Cameroon. J Med Virol. 2013;85:1485–90.
Taniguchi K, Urasawa T, Kobayashi N, Gorziglia M, Urasawa S. Nucleotide sequence of VP4 and VP7 genes of human rotaviruses with subgroup I specificity and long RNA pattern: implication for new G serotype specificity. J Virol. 1990;64:5640–4.
Rahman M, Matthijnssens J, Yang X, et al. Evolutionary history and global spread of the emerging g12 human rotaviruses. J Virol. 2007;81:2382–90.
Abugalia M, Cuevas L, Kirby A, et al. Clinical features and molecular epidemiology of rotavirus and norovirus infections in Libyan children. J Med Virol. 2011;83:1849–56.
Kiulia NM, Peenze I, Dewar J, et al. Molecular characterisation of the rotavirus strains prevalent in Maua, Meru North, Kenya. East Afr Med J. 2006;83:360–5.
Shoja Z, Tagliamonte M, Jalilvand S, et al. Molecular characterization analysis of the outer protein layer (VP7) from human rotavirus A genotype G1 isolate identified in Iran: implications for vaccine development. New Microbiol. 2012;35:415–27.
Mwenda JM, Tate JE, Parashar UD, et al. African rotavirus surveillance network: a brief overview. Pediatr Infect Dis J. 2014;33(Suppl 1):S6-8.
Moure UAE, Banga-Mingo V, Gody JC, Mwenda JM, Fandema J, Waku-Kouomou D, et al. Emergence of G12 and G9 rotavirus genotypes in the Central African Republic, January 2014 to February 2016. BMC Res Notes. 2018;11(1):5.
We thank the staff of CDC, Atlanta for the technical support, all the children and their parents who participated in this study. We express our gratitude to Dr. Leonard Peruski, CDC/Lead, International Laboratory Operations, for his inputs in the writing of this paper.
The authors have no funding to disclose. Funding for genotyping of specimens was provided by the CDC.
Ethics approval and consent to participate
This study was authorized by the ethics committee of the Ministry of Health of Benin. We obtained a valid administrative authorization no. 188/16/MS/DDSATL-LIT/ZS-AS/SAAE/DGAP/SGA dated 15th June 2016, before collecting samples from hospitals in Cotonou (Menontin/Anastasie/Ste Vincent de Paul) and Abomey-Calavi (Central Clinic of Abomey-Calavi) in South Benin (West Africa). Stool samples were taken after written informed consent from parents or guardians of children.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Distribution of the study population by age and gender. Figure S1a. Sample collection sites. (1) Abomey Calavi in the Atlantic Region (a = Central Clinic of Abomey Calavi) and (2) Cotonou in the Central Region (b = Anastasis Hospital, c = Polyclinic St Vincent de Paul and d = Mênontin Hospital). Figure S1b. Distribution of genotypes by age in Cotonou. Figure S1c. Distribution of genotypes from 90 stool samples in Abomey-Calavi. Figure S1d. Distribution of genotypes according to age in Abomey-Calavi.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Agbla, J.M., Esona, M.D., Agbankpe, A.J. et al. Molecular characteristics of rotavirus genotypes circulating in the south of Benin, 2016–2018. BMC Res Notes 13, 485 (2020). https://doi.org/10.1186/s13104-020-05332-7