Molecular characteristics of rotavirus genotypes circulating in the south of Benin, 2016-2018

doi:10.21203/rs.3.rs-51311/v4

PDF

Research note

Molecular characteristics of rotavirus genotypes circulating in the south of Benin, 2016-2018

https://doi.org/10.21203/rs.3.rs-51311/v4

This work is licensed under a CC BY 4.0 License

Journal Publication

published 19 Oct, 2020

Read the published version in BMC Research Notes →

You are reading this latest preprint version

Objective

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).

Results

Of the 420 diarrheal stool samples, actively collected in southern Benin from July 2016 through November 2018 from children under five 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[8] represents the predominant genotype (32%) followed by the G2P[4] (26%), G3P[6] (16%), G12P[8] (13%) and mixed G and P types (1%). Four samples (2%) could not be assigned both G and P type speciﬁcity.

Health Economics & Outcomes Research

Health Policy

pediatric

rotavirus

surveillance

Genotypes

Benin

Diarrhea ranks as the fifth leading cause of mortality among children under five years old (1). Infectious diarrhea can be caused by microbial agents of bacterial, parasitic, viral or mycotic origin (2). 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 five years of age and the largest number of rotavirus deaths occurred in sub-Saharan Africa (5).

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[8], G2P[4], G3P[8], G4P[8] and G9P[8] 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) (14). 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) (15).

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 five 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 (Figure S1a; Additional file). Children with bloody diarrhea and those older than five 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 10⁹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 (16). 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-3hrs 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 1K or 5K reagent kit (Dual protocol DNA Analysis and Quantitation) with the HT Extended Range LabChip (Caliper Life Sciences, MA, USA) as described previously (16).

Demographic characteristics and genotyping results.

Distribution of the study population by age and sex.

In this study, children aged from 7-12 months were the most represented with a percentage of 41.9% followed by those under six months of age (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).

Genotyping Results.

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 speciﬁcity (non-typeable samples). G1P[8] represents the predominant genotype (32%) circulating in southern Benin, followed by the G2P[4] (26%), G3P[6] (16%) and G12P[8] (13%). Other uncommon genotypes in Benin such as G9P[8] (3%), G1P[6] (3%), G3P[8] (2%), G2P[8] (1%), G2P[6] (1%), G12P[6] (1%), G9P[4] (1%) and mixed genotypes (1%) were detected, albeit at low rate (Figure 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 (Figure 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[8] (30%) followed by G2P[4] (23%), G3P[6] (16%), and G12P[8] (15%) (Figure 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[6] (4%), G9P[8] (3%), G3P[8] (2%), G2P[8] (1%) G2P[6] (1%), G12P[6] (1%) and G9P[4] (1%) were detected.

There is great variability of genotypes in children aged from 0-18 months. This variability is less after 18 months (Figure S1b; Additional file).

Distribution of genotypes in Abomey-Calavi.

In Abomey-Calavi, 39 (43.3%) stool samples were positive by RT-PCR. G1P[8] and G2P[4] genotype combinations were the most common (38% and 36% respectively), followed by G3P[6] (15%), G12P[8] (5%). No mixed genotypes were found and one sample was untypeable (Figure S1c; Additional file).

The greatest variability of genotypes was more observed in children less than 18 months in Abomey-Calavi. (Figure S1d; Additional file).

In the two localities that hosted the study, G1P[8] 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 (17) describing the epidemiological situation in Africa, G1P[8] 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 (17).

At country level in West Africa, before vaccine introduction, Ghana reported G1P[8] (20%) as the predominant genotype (18), but others genotypes such as G6P[6] was found to be medically important in Burkina Faso (19), G12P[8] in Côte d’Ivoire (20) and G4P[8] and G12P[8] in Nigeria (21, 22). During the post vaccine period, Ghana reported G1P[8] as the fourth common genotype detected with a rate of 8.0% (18). 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[8] 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 (23). On the other hand, according to the study of Seheri and colleagues, 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 (24). 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 and colleagues showed that G1P[8] was the most predominant genotype, followed by G2P[4], G9P[8] and G2P[6] (17). Similarly, in our study, G2P[4] 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[4] was the predominant genotype before vaccine introduction (25, 26). However, G2P[4] 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-30).

G12P[8] is the second emerging genotype, detected for the first time in Philippines in 1987, in children under 2 years of age (31). Mainly limited to Asian countries, G12P[8] has been reported worldwide over the last 20 years (32).

Commonly detected at low rates in Africa, over the past decade, G12P[8] 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 (18). Also, according to some studies, G12P[8] may come from a reassortment between human and pig rotavirus strains (32). The presence of this genotype in our study without any pre-existing vaccination context, strongly suggests that G12P[8] appearance in Benin, is due to natural rotavirus genotype occurrence, not vaccine pressure.

On the other hand, G9P[8] 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 (35).

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%) (36). 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% (36), 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 (20). However, studies from some Africans settings have shown no mixed genotypes (37) or a very low rate (1%) of mixed infections (20) 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 and colleagues (24). 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.

EIA: Enzyme linked Immunosorbent Assay
NIP: National Immunization Programs
qRT-PCR: quantitative Reverse Transcription-Polymerase Chain Reaction
RT-PCR: Reverse Transcription-Polymerase Chain Reaction
WHO: World Health Organization.

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 15^th 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

Not applicable

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: [email protected] and materials deposition are not applicable in this study

COMPETING INTERESTS

The authors declare that they have no competing interests

FUNDING

The authors have no funding to disclose. Funding for genotyping of specimens was provided by the CDC.

AUTHORS' CONTRIBUTIONS

All authors read and approved the manuscript before submission. JMA, ACC, AJA, TVD, OR were involved in samples collection process, JMA, AJA, TVD, OR performed the EIA and MDE, RG, MDB performed the genotyping. JMA, AJA, MDE, RG, MDB contributed in data analysis and writing of the first draft of the article; HSB, MDE, ACC, TVD, AJA, OR review the article and JMA, MDE, AJA, MDB, HSB were major contributors in writing the manuscript. MDE, RG, MBD performed the RT-PCR quality control and HSB, MDB, MDE supervised the work. All authors read and approved the final manuscript.

ACKNOWLEDGEMENTS:

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.

CONTRIBUTOR’ INFORMATIONS

Jijoho Michel Agbla, [email protected]
Mathew D. Esona, [email protected]
Alidéhou Jerrold Agbankpé, [email protected]
Annick Capo-Chichi, [email protected]
Rashi Gautam, [email protected]
Tamègnon Victorien Dougnon, [email protected]
Osseni Razack, [email protected]
Michael D. Bowen, [email protected]
Honoré Sourou Bankolé, [email protected]

Collaborators. GDD. 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-1228.
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-2161.
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-141.
(WHO) WHO. Rotavirus vaccines. http://wwwwhoint/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-S105.
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-663.
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-1373.
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-282.
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-433.
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-391.
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-7123.
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-7242.
Nordgren J, Bonkoungou IJ, Nitiema LW, et al. Rotavirus in diarrheal children in rural Burkina Faso: high prevalence of genotype G6P[6]. Infect Genet Evol. 2012;12:1892-1898.
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: Clinical manifestations, high G12P[8] frequency, intra-genotypic divergence of VP4 and VP7. J Med Virol. 2018;90:239-249.
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-7247.
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-7230.
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-580.
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-359.
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-479.
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-1460.
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-95.
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-1490.
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-5644.
Rahman M, Matthijnssens J, Yang X, et al. Evolutionary history and global spread of the emerging g12 human rotaviruses. J Virol. 2007;81:2382-2390.
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-1856.
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-365.
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-427.
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 Jan 5;11(1):5.
Boula A, Waku-Kouomou D, Njiki Kinkela M, Esona MD, Kemajou G, Mekontso D, et al. Molecular surveillance of rotavirus strains circulating in Yaoundé, Cameroon, September 2007-December 2012. Infect Genet Evol. 2014 Dec;28:470–5.

BMCSupplementarymaterialRevised.pptx

PDF

Journal Publication

published 19 Oct, 2020

Read the published version in BMC Research Notes →

Submission checks completed at journal
07 Oct, 2020
Editorial decision: Accept
06 Oct, 2020

You are reading this latest preprint version

Molecular characteristics of rotavirus genotypes circulating in the south of Benin, 2016-2018

Status:

Journal Publication

Version 4

Abstract

Figures

Introduction

Methods

Results

Discussion

Limitations

List Of Abbreviations

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 4