- Research article
- Open Access
Rotavirus group A genotype circulation patterns across Kenya before and after nationwide vaccine introduction, 2010–2018
BMC Infectious Diseases volume 20, Article number: 504 (2020)
Kenya introduced the monovalent G1P  Rotarix® vaccine into the infant immunization schedule in July 2014. We examined trends in rotavirus group A (RVA) genotype distribution pre- (January 2010–June 2014) and post- (July 2014–December 2018) RVA vaccine introduction.
Stool samples were collected from children aged < 13 years from four surveillance sites across Kenya: Kilifi County Hospital, Tabitha Clinic Nairobi, Lwak Mission Hospital, and Siaya County Referral Hospital (children aged < 5 years only). Samples were screened for RVA using enzyme linked immunosorbent assay (ELISA) and VP7 and VP4 genes sequenced to infer genotypes.
We genotyped 614 samples in pre-vaccine and 261 in post-vaccine introduction periods. During the pre-vaccine introduction period, the most frequent RVA genotypes were G1P  (45.8%), G8P  (15.8%), G9P  (13.2%), G2P  (7.0%) and G3P  (3.1%). In the post-vaccine introduction period, the most frequent genotypes were G1P  (52.1%), G2P  (20.7%) and G3P  (16.1%). Predominant genotypes varied by year and site in both pre and post-vaccine periods. Temporal genotype patterns showed an increase in prevalence of vaccine heterotypic genotypes, such as the commonly DS-1-like G2P  (7.0 to 20.7%, P < .001) and G3P  (1.3 to 16.1%, P < .001) genotypes in the post-vaccine introduction period. Additionally, we observed a decline in prevalence of genotypes G8P  (15.8 to 0.4%, P < .001) and G9P  (13.2 to 5.4%, P < .001) in the post-vaccine introduction period. Phylogenetic analysis of genotype G1P , revealed circulation of strains of lineages G1-I, G1-II and P -1, P -III and P -IV. Considerable genetic diversity was observed between the pre and post-vaccine strains, evidenced by distinct clusters.
Genotype prevalence varied from before to after vaccine introduction. Such observations emphasize the need for long-term surveillance to monitor vaccine impact. These changes may represent natural secular variation or possible immuno-epidemiological changes arising from the introduction of the vaccine. Full genome sequencing could provide insights into post-vaccine evolutionary pressures and antigenic diversity.
Childhood diarrhea caused by rotavirus group A (RVA) infection remains a leading cause of morbidity and mortality in young children globally . In 2016, RVA infections were estimated to be responsible for 1,537,000 hospitalization cases, 128,500 deaths globally and over 80% of these deaths occurred in developing countries . In Kenya alone, it is estimated that rotavirus infection accounts for over 3000 deaths annually in children under 5 years of age .
Upon infection by the virus, immune response to RVA by the host is directed to the highly variable VP7 and VP4 genes found on two separate segments of the double-stranded RNA genome . RVA G and P genotypes exist as multiple variants in nature, few of which have been found to infect humans . Up to 36 G and 51 P genotypes have been detected globally in both humans and animals, with multiple G-P combinations . Molecular studies have characterized circulating genotypes worldwide with predominance of genotypes G1P , G2P , G3P , G4P , G9P  and G12P  (in decreasing order of prevalence) [6, 7]. Although the distribution of these genotypes varies from region to region and from one season to another, genotype G1P  has remained the most dominant genotype globally . In Africa, there is a high diversity of genotypes, most commonly G1P , G2P , G9P , G2P , G12P  and G3P , with G1P  and G2P  predominant [7, 8]. Significant RVA infections are also caused by strains of uncommon genotypes including G1P , G2P , G9P , G12P , G8P , G8P  and G12P  [7, 8]. Such uncommon strains also show a wide variation from one region to the other, and can spread globally to become common strains. For instance, genotypes G9P  and G12P  emerged and contributed to a larger proportion of global RVA infections . An understanding of these genetic diversity after vaccine introduction is necessary for design and implementation of effective control programs.
In 2009, the World Health Organization (WHO) recommended the inclusion of RotaTeq® (Merck Vaccines, Whitehouse Station, New Jersey) or Rotarix® (GlaxoSmithKline Biologicals, Rixensart, Belgium) vaccines in the national immunization programs of countries that experience high diarrhea morbidity and mortality burden due to RVA disease . The introduction and increased use of these vaccines have reduced up to 76% of rotavirus hospitalizations in children < 5 years  and averted up to 28,800 deaths globally, including 84% of the deaths in sub-Saharan Africa . Furthermore, these vaccines showed overall good clinical protection against multiple homotypic and heterotypic RVA strains in humans . Kenya incorporated the monovalent G1P  Rotarix® vaccine into the national immunization program in July 2014, administered in two oral doses offered at weeks 6 and 10 of age. The Rotavirus Immunization Programme Evaluation in Kenya (RIPEK) was established as a collaboration among existing rotavirus surveillance platforms across Kenya to monitor the impact of Rotarix vaccine introduction against rotavirus disease and circulating RVA genotypes. Substantial effectiveness of the vaccine in Kenya (vaccine coverage of 72% ) and the entire sub-Saharan region (where disease burden is high) has been recorded, and the decline in incidence of all-cause and rotavirus associated diarrhea admissions has been attributed to vaccine implementation [14,15,16,17,18]. However, there are limited data on RVA diversity in post-vaccine introduction periods in this region. The current report describes the distribution and temporal patterns of RVA genotypes observed before and after Rotarix vaccine introduction in Kenya.
RVA surveillance was carried out in four health facilities in coastal, western and central regions of Kenya. These surveillance sites were: Kilifi County Hospital (KCH) in Kilifi County, Tabitha Clinic (TC), in Kibera, Nairobi County, Saint Elizabeth Lwak Mission Hospital (LMH) and Siaya County Referral Hospital (SCRH) in Siaya County (Supplementary Figure 1). KCH is a government hospital located on the Kenyan coast serving a rural and semi-rural population. Rotavirus surveillance at KCH started in September 2009, and is implemented by KEMRI Wellcome Trust Research Program (KWTRP) a partnership between Kenya Medical Research Institute and Wellcome Trust, UK [19, 20]. SCRH is a government hospital which also serves a rural and semi-rural population in Siaya County, Western Kenya. SCRH has had an active rotavirus surveillance since 2010 and is implemented by KEMRI Center for Global Health Research, KEMRI-CGHR in collaboration with US Centers for Disease Control and Prevention (CDC) . LMH is a private health facility serving a rural population in Asembo, Siaya County where rotavirus surveillance is carried out as part of the Population Based Infectious Disease Surveillance (PBIDS) platform under KEMRI-CGHR and CDC . PBIDS also supports rotavirus surveillance at Tabitha Clinic, a private health facility in Kibera, Nairobi County.
Surveillance was conducted during January 2010 – December 2018 for all the sites except for SCRH where surveillance ended in December 2016. Stool samples were collected from children aged < 13 years (in KCH, TC, LMH) and < 5 years (in SCRH) of age presenting with acute gastroenteritis (AGE). AGE was defined as ≥3 watery stools passed within a 24-h period during the illness for KCH, TC and LMH, while for SCRH, AGE was defined as ≥3 loose stools and/or ≥ 1 episode of unexplained vomiting followed by loose stool within a 24-h period beginning no more than 7 days before the visit to SCRH.
RVA was tested by use of commercially available enzyme immunoassays. The ProSpecT™ Rotavirus Kit (Oxoid, Basingstoke UK) was used to test samples collected from KCH while the Rotaclone® kit (Premier™ Meridian Bioscience, Cincinnati, Ohio, USA) was used to test samples collected from LMH, SCRH, and TC. For samples collected from LMH, KCH and TC, partial fragments of the segments encoding the outer capsid proteins, VP4 (660 bp) and VP7 (881 bp), were amplified in a One-step Reverse Transcriptase Polymerase Chain Reaction (RT-PCR) using previously described primer pairs [23, 24]. Successful amplification was visualized by electrophoresis of the PCR product in a 2% agarose gel. PCR products of confirmed positives were purified using GFX DNA purification kit (GFX-Amersham, Amersham, UK), according to the manufacturer’s instructions. Confirmed positives were then sequenced using Big Dye Terminator 3.1 (Applied Biosystems, Foster City, California, USA) with the same primers as in PCR amplification on an ABI Prism 3130xl Genetic Analyzer (Applied Biosystems, Foster City, California, USA).
Reads from the sequencer were trimmed (removing regions including primer sequences) and assembled into contigs (consensus sequence formed from aligning the forward and reverse reads) using Sequencher version 5.4.6 (Gene Codes Corp Inc., Ann Arbor, MI, USA). The cleaned sequences for the VP7 and VP4 genes have been deposited to GenBank database: https://www.ncbi.nlm.nih.gov/genbank/. Sequence accession numbers for these genes are provided in supplementary Table 2 and supplementary Table 3, respectively. G and P genotypes were determined by submitting cleaned sequences to the online automated Virus Pathogen Resource genotyping tool . Using this method only a single genotype per specimen was identified. However, 8.6% (n = 113) of the samples were typed for only one of the genes due to failure in sequencing and/or contig assembly, supplementary Table 1, and were excluded from the main analysis. For positive samples collected from SCRH, sample processing and genotyping was performed as previously described . Briefly, full length VP7 gene (1062 bp) was amplified using sBeg9 and End9 primers. The resulting cDNA was used for G typing using primers End9, aAT8v, aBT1, aCT2, aDT4, mG3, mG9, mG10 and G12 for G1, G2, G3, G4, G8, G9, G10 and G12 genotypes, respectively [27, 28]. Alternatively, partial-length amplification of the VP4 gene (876 bp) was achieved by use of Con2 and Con3 primers. P typing was performed using a mixture of primers consisting of Con3 and primers 1 T-1, 2 T-1, 3 T-1, 4 T-1, 5 T-1, mP  and P  for genotypes P , P , P , P , P , P  and P , respectively [23, 29]. In this instance, untypable G (8.9%, n = 15/168) and P (22.6%, n = 38/168) genotypes and mixed infections (18.5%, n = 31/168) were observed and subsequently excluded from analysis.
Phylogenetic analysis of the Kenyan G1 and P  vaccine strains
Maximum likelihood (ML) method was used to determine the phylogenetic relationship of the partial sequences of genotypes G1 and P  observed during the pre- and post-vaccine periods in Kenya, and further compared to RVA strains circulating globally. Global contemporaneous sequences were retrieved from GenBank (as of May 2020) and compared to local strains. Sequences were aligned using MAFFT v7.2  and visualized in AliView v1.8 . Best-fit volutionary models were tested and selected in IQTREE v1.6  using the Bayesian Information Criteria . ML trees were inferred using IQTREE with 1000 bootstrap replicates. The resulting trees were visualized and edited in FigTree v1.4.3 (http://tree.bio.ed.ac.uk/software/figtree/). Nucleotide distances matrixes were prepared using the p-distance algorithm inferred in MEGA v10 . The trees were drawn to scale indicating nucleotide substitutions rates per site.
Data collating and analyses were performed in Microsoft Excel and R version 3.5. Genotype data collected from LMH and SCRH were merged to constitute data from Siaya County. The summary of proportions was conducted for each genotype by site and period of detection. Frequency distribution and temporal pattern graphs were generated. Two-sided proportional tests of the most common genotypes during the pre-vaccine (January 2010 thru June 2014) and post-vaccine (July 2014 thru December 2018) eras were performed in R, and P-values < 0.05 were considered statistically significant.
A total of 6562 stool samples were collected from health facilities in LMH, SCRH, TC and KCH. Of the 1312 (20.0%) RVA positive samples, individual G and P genotypes were identified in 928 (70.7%) and 904 (68.9%) samples, respectively, while 875 (66.7%) samples were successfully genotyped for both G and P genes [LMH (n = 38), SCRH (n = 161), TC (n = 84), and KCH (n = 592)]. Of these, 614 (70.2%) were observed in the pre-vaccine (January 2010 – June 2014) and 261 (29.8%) in the post-vaccine (July-2014 – December 2018) introduction periods, respectively.
G and P genotypes circulating in Kenya
Overall, nine different G types (G1, G2, G3,G4, G8, G9, G10, G12, G29) and four P types (P , P , P  and P ) were observed in Kenya. The most common G type was G1 (49.5%) followed by G9 (12.7%), G8 (12.0%), G2 (11.1%) and G3 (10.1%), Table 1. Genotypes G4, G10, G12 and G29 were detected at low frequencies (< 4%). The most common P types were P  (68.4%) and P  (23.6%), while P  (7.7%) and the less common P  (0.3%) were also observed, Table 1.
Rotavirus genotype distribution in pre- and post- vaccination eras
In the entire period, 22 different G/P combinations were identified, Table 2. G1P  (47.7%) was the most common genotype followed by G8P  (11.2%), G2P  (11.1%), G9P  (10.9%), G3P  (5.7%) and G3P  (2.9%). In addition to these common genotypes, multiple other rare genotypes; G1P , G1P , G1P , G2P , G3P , G3P , G4P , G8P , G8P , G9P , G9P , G10P , G12P , G12P , G12P  were observed in low frequency (< 2%), Table 2.
Figure 1 a, b, c and d show the temporal distribution of RVA genotypes in Kilifi, Siaya (i.e, pooled SCRH and LMH data) and Nairobi counties, and Kenya (pooled countrywide), respectively. Genotypes G1P , G2P , G8P , G3P , G3P  and G9P  were observed across all three sites. G1P  was observed in all periodss, in all sites, except in 2011 and 2016 in Nairobi. In contrast, genotypes G2P , G3P , G3P , G8P  and G9P  showed a fluctuating pattern in all three sites. Notably, genotype G2P  was observed in pre- and post-vaccine periods in Kilifi and Nairobi, and only in the post- vaccine period in Siaya. Furthermore, genotype G3P  was observed in Kilifi and Nairobi counties in pre- and post-vaccine periods and only in pre-vaccine period in Siaya, while genotype G3P  was detected in all three sites, although in moderate proportions. Additionally, genotype G9P  was observed in pre- and post- vaccine periods in Kilifi, unlike in Nairobi and Siaya where it was only observed in the pre-vaccine period. Genotype G8P  occurred in pre- and post-vaccine periods in Nairobi county and only in the pre-vaccine period in Kilifi and Siaya counties.
During the period before vaccine introduction in Kenya, G1P  (45.8%) was the predominant genotype observed in this population, followed by G8P  (15.8%), G9P (13.2%), G2P  (7.0%) and G3P  (3.1%). After vaccine introduction, G1P  remained the dominant genotype (52.1%), followed by G2P  (20.7%), G3P  (16.1%), G9P  (5.4%) and G3P  (2.7%). Multiple other genotypes were also observed in pre- and post-vaccine periods, although in low proportions (< 2%), Tables 2, 3.
Overall, at these surveillance sites, the first four RVA periods, saw an alternating pattern of dominance between G8P  and G1P  (Fig. 1d). G8P  was the predominant genotype in 2010 (54/153; 35.3%) and in 2012 (42/150;28.0%), while G1P  predominated in 2011 (117/160;73.1%) and 2013 (60/104;57.7%) through to 2014 (vaccine introduction periods) breaking the cyclic pattern. During this periods, G2P  and G9P  circulated in moderate proportions (1–20%) and (2–14%), respectively (Table 2). After vaccine introduction, G1P  dominated immediately in the first year of vaccine introduction (2015) (91/101, 90.1%). This phenomena was however, short lived and there was a re-emergence of G2P  (35/44; 79.5%) in the second year of vaccine period (2016), which was among the common genotypes in 2017 together with genotype G1P (7/16;43.6%). G3P  (41/66;62.1%) dominated in the fourth year post-vaccine introduction period (2018), whereas G1P  (17/66;25.8%) and G2P  (4/66;6.1%) continuously circulated, although in reduced proportions. Notably, a separate analysis of the post-vaccine G3P  (data not shown), revealed that 4.8% (n = 2/41) were of the equine-like genotype while the rest were wild-type genotype.
In the post-vaccine introduction period, relative to the pre-vaccine period, we observed a significant increase in prevalence of genotypes G2P  (7.0% vs. 20.7%, P < .001) and G3P  (1.3% vs. 16.1%, P < .001) and significant decline in prevalence of genotypes G8P  (15.8% vs. 0.4%, P < .001) and G9P  (13.2% vs. 5.4%, P < .001), Fig. 2. No significant difference in the prevalence of the vaccine genotype G1P  (45.8 vs. 52.1%,P = .35) was measured before and after vaccine introduction. A decline in genotype diversity was observed after vaccine introduction (11 genotypes) as compared to pre-vaccine period (21 genotypes).
Phylogenetic analysis of the Kenya G1 and P  vaccine strains
As shown in Fig. 3a, phylogenetic analysis of Kenyan G1 strains revealed circulation of lineages G1-I and G1-II with representative sequences from all the three sites (Kilifi, Kibera and Lwak). In both pre and post-vaccine periods, detected strains were predominantly found within lineage G1-I with a few within lineage G1-II. Interestingly, within lineage G1-I, pre and post vaccine strains formed separate clusters. However, three post-vaccine strains (2015; n = 2, 2014; n = 1) clustered within pre-vaccine strains, which showed high nucleotide sequence identities of up to 99% to pre-vaccine strains. Additionally, a single pre-vaccine strain, observed in 2014, grouped within the post-vaccine cluster, which showed high nucleotide identities of up to 100% to other post-vaccine strains. The lineage G1-I strains also clustered with other strains isolated from Kenya (2009), Malawi (2012), Togo (2011), Russia (2013), India (2013), Kuwait (2016), USA (2011) and Italy (2015), with high nucleotide identities ranging between 98 and 99%. Lineage II cluster contained strains observed in pre-vaccine period (n = 14) which clustered with other globally isolated strains from Pakistan (2013), India (2012), USA (2013) and Hungary (2012), with nucleotide identities of up to 99%. In contrast, a single post-vaccine strain, isolated in Kibera (KBR337), was observed within lineage G1-II cluster, which showed high nucleotide and amino acid identity of up to 100% to the Rotarix vaccine and other vaccine-derived strains from Japan, Brazil and Korea.
As shown in Fig. 3b, analysis of the Kenyan P  genotype revealed groupings into lineage P -III (predominant lineage) with a few others belonging in lineage P -I and P -IV. Within lineage P -III, distinct clusters of pre and post-vaccine strains were evident with clusters containing genotypes from all the three study sites (Kibera, Lwak and Kilifi). However, three post-vaccine strains (observed in 2015) clustered with pre-vaccine strains, while one pre-vaccine strain (observed in 2014) clustered within post-vaccine genotypes. The strains in lineage III clustered with other strains isolated from Kenya (2009), Togo (2011), Congo (2010), Tanzania (2011), India (2011), USA (2011), Pakistan (2010) and Kuwait (2016) showing nucleotide identities ranging from 97 to 98%. The P -I cluster contained a single post-vaccine strain (KBR337), from Kibera, which showed 100% nucleotide and amino acid identity to the Rotarix vaccine and other vaccine-derived strains observed in Vietnam and Japan. Lineage P -IV cluster contained strains observed in pre (2012 and 2014, from Kibera and Lwak) and post (2018, from Kibera) vaccine periods and clustered with isolates from Russia (2010).
Prior to vaccine introduction, G1P  was the predominant genotype. However, this genotype’s prevalence varied substantially from year to year (dominant in 2011, 2013, 2014 and 2015; less so in 2010 and 2012), emphasizing the limits of short-term surveillance and the potential for misrepresentation of patterns. During 2010 and 2012, a large proportion of children were infected with genotype G8P , which was not detected in any years following vaccination (except during 2017). Genotypes G1P , G2P  and G9P  have also been identified as dominant genotypes in a study involving six countries from Eastern and Southern Africa .
Genotype prevalence varied year to year and from before to after vaccine introduction. Interestingly, while genotype G1P  showed no overall change in prevalence compared to the pre-vaccine period, it predominated only in the first year after vaccine introduction (2015), and declined in prevalence thereafter, in particular to be replaced by G2P  in 2016, and by G3P  in 2018. Genotypes G8P  and G9P , which were the second and third common genotypes in pre-vaccine period, significantly declined in the post-vaccine period. Predominance of G2P  and G3P  has been reported worldwide following universal Rotarix vaccine introduction. For instance, G2P  was observed in Australia  (in states using Rotarix) and Belgium  after vaccine introduction, while data from Ethiopia  and Madagascar  depicted an increase in prevalence of G3P  in post-vaccine introduction period. Studies on the G3P  have revealed the emergence and spread of strains of equine-like G3 genotype co-circulating with wild-type G3 strains [39,40,41]. In here, 4.8% (2/41) were detected to be of equine-like type (data not shown). Further analysis on this genotype will illuminate of the genetic diversity and prevalence of the emerging equine-like in the post-vaccine introduction period. G2P  is one of the genotypes of concern since this fully heterotypic genotype has a different genomic constellation (DS-1-like) to the genotype in the monovalent Rotarix vaccine (Wa-like) . There is no conclusive evidence associating the increase in prevalence of these genotypes to selective vaccine pressure created by implementation of Rotarix monovalent vaccine. Whilst the vaccine offers both homotypic and heterotypic protection [43, 44], the emergence of the fully heterotypic G2P  and partially heterotypic G3P  genotypes, and persistence in circulation of the homotypic G1P  genotype after vaccine introduction emphasizes the need for continued monitoring of impact of vaccine on genotypes.
We observed multiple uncommon G/P combinations including G1P , G1P , G2P ,G4P , G8P ,G8P , G10P , G12P  and G12P  at frequency of < 2%. Genotype P  has been described sporadically in humans and is believed to have originated from animal rotavirus and introduced into humans through interspecies transmission and/or reassortment events . The circulation of these uncommon genotypes demonstrates the high diversity of RVA genotypes in Kenya and concurs with findings from other African countries . However, genotype diversity was seen to decline in post-vaccine period (only 11 GP combinations, compared to 21 in pre-vaccine period), mirroring the experience of other countries, including Brazil and Zimbabwe, which also indicated a decline in genotype variation after vaccine introduction [46, 47]. This is in contrast to other countries such Australia which observed increase in genotypes after vaccine introduction . Although diversity seems to decrease following vaccine introduction, it is unknown whether the observed trends will be sustained in the long-term, especially in African settings where pre-vaccine genotype diversity is high.
Phylogenetic analysis of genotype G1 revealed circulation of strains of lineages G1-I and G1-II while the P  strains formed three distinct clusters revealing circulation of strains of lineages P –I, P -III and P -IV. The Kenyan strains were closely related to other global strains as they segregated together. G1-I and P -III were the predominant lineages for both pre and post-vaccine strains. These lineages were distinct from the Rotarix vaccine lineages showing that the strains were distantly related to the vaccine strain. Similar findings were observed in data collected before and after vaccine introduction in Eastern and Southern Africa countries . Although there were no lineage replacements after vaccine introduction, the genetic diversity between pre and post-vaccine strains, may suggest normal genetic fluctuations or an effect of the vaccine. Interestingly, a single case of vaccine-derived strain was observed in a sample isolated from a child in Kibera, who had received the first dose of the Rotarix vaccine. Phylogenetically, the post-vaccine strain clustered with the Rotarix vaccine, in G1-I and P -I lineages, and showed sequence identities of 100% to the vaccine strain. These clusters also contained previously vaccine-derived strains isolated in Brazil, Korea, and Japan. Acute gastroenteritis, caused by vaccine shedding or horizontal transmissions of the vaccine strain have been reported at a higher rate in other countries [48,49,50]. Further analysis will be necessary to determine whether the observed vaccine strain was as a result of vaccine shedding or horizontal transmission event.
This study provides substantial epidemiological information on changes in distribution and genetic diversity of RVA genotypes in Kenya. Some limitations of this study included fewer surveillance sites and unequal distribution of the number of samples among the sites, potentially underestimating the type of genotypes circulating in the post-vaccine period. However, it is important to consider that there was an overall decline in rotavirus circulation in the post-vaccine era, hence the fewer specimens and less diversity could be a reflection of less rotavirus circulation [14, 15]. Additionally, since this is an ecological study, the changes in distribution and diversity of genotypes in the post-vaccine era cannot directly be attributed to vaccine introduction. Furthermore, due to unsuccessful sequencing and/or contig assembly only two thirds of the total RVA positive samples were fully genotyped. Lastly these findings may not be generalized to the whole country because analysis was based on rotavirus cases observed in health facilities in three counties only.
In conclusion, we highlight the importance of monitoring the transition in the prevalence of genotypes for a better understanding of the performance of the currently available vaccines. The emergence of the fully heterotypic G2P  and partially heterotypic G3P  genotypes after vaccine introduction raises questions about the epidemiological dynamics following vaccine introduction. Previous analysis in Kenya showed that the vaccine had a significant impact on G1P  and non-significant G2P  (although with limited statistical power) , hence, continuous monitoring of the circulating genotypes in the post-vaccine era is necessary. Our findings also highlight existence of considerable variation and genetic diversity within and between Kenyan pre and post-vaccine strains. One strain identified in Kenya, was closely related to the Rotarix® vaccine strain, likely representing shedding or horizontal transmission of the vaccine strain. Additionally, continued surveillance of the genetic characteristics of circulating RVA strains is recommended to obtain a full view of the long-term effects of vaccine introduction. Since immunity to RVA involves immune responses conferred by genes other than the commonly studied P and G genes, vaccine effectiveness might be challenged by changes occurring on non-capsid genes. It is therefore recommended that full genome analysis of genotypes collected in different time or geographic regions be conducted to improve understanding of their evolutionary profile during the post-vaccine introduction period.
Availability of data and materials
Partial sequences for the VP7 and VP4 genes reported in this work have been deposited to GenBank database: https://www.ncbi.nlm.nih.gov/genbank/. Sequence accession numbers for the VP7 and VP4 genes are indicated in supplementary Table 2 and supplementary Table 3, respectively. The datasets used and analyzed during this study are available from the Harvard Dataverse: https://doi.org/10.7910/DVN/0VQ2OK. Public access to the data is restricted. Users who wish to access and use the data should send a request to the KEMRI Wellcome Trust Research Programme data governance committee, which can be contacted by emailing: email@example.com.
Rotavirus group A
World Health Organization
Rotavirus Immunization Programme Evaluation in Kenya
Troeger C, Forouzanfar M, Rao PC, Khalil I, Brown A, Reiner RC, et al. Estimates of global, regional, and national morbidity, mortality, and aetiologies of diarrhoeal diseases: a systematic analysis for the global burden of disease study 2015. Lancet Infect Dis [Internet]. 2017;17(9):909–48 Available from: https://www.thelancet.com/journals/laninf/article/PIIS1473-3099(17)30276-1/fulltext.
Troeger C, Khalil IA, Rao PC, Cao S, Blacker BF, Ahmed T, et al. Rotavirus vaccination and the global burden of rotavirus diarrhea among children younger than 5 years. JAMA Pediatr. 2018;98121:1–8 Available from: http://archpedi.jamanetwork.com/article.aspx?doi=10.1001/jamapediatrics.2018.1960.
Tate JE, Burton AH, Boschi-Pinto C, Parashar UD, Agocs M, Serhan F, et al. 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 Available from: https://academic.oup.com/cid/article/62/suppl_2/S96/2478843.
Desselberger U. Rotaviruses. Virus Res. 2014;190:75–96 Available from: https://www.sciencedirect.com/science/article/pii/S0168170214002640?via%3Dihub.
Metagenomics V. Virus classification [internet]. 2018. Available from: https://rega.kuleuven.be/cev/viralmetagenomics/virus-classification.
Bányai K, László B, Duque J, Steele AD, Nelson EAS, Gentsch JR, 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–30 Available from: https://linkinghub.elsevier.com/retrieve/pii/S0264410X11015532.
Seheri LM, Magagula NB, Peenze I, Rakau K, Ndadza A, Mwenda JM, et al. Rotavirus strain diversity in Eastern and Southern African countries before and after vaccine introduction. Vaccine. 2017; Available from: https://doi.org/10.1016/j.vaccine.2017.11.068.
Mwenda JM, Ntoto KM, Abebe A, Enweronu-Laryea C, Amina I, Mchomvu J, et al. Burden and epidemiology of rotavirus diarrhea in selected African countries: preliminary results from the African Rotavirus Surveillance Network. J Infect Dis. 2010;202(Suppl 1):S5–11 Available from: https://academic.oup.com/jid/article/202/Supplement_1/S5/849136.
Matthijnssens J, Bilcke J, Ciarlet M, Martella V, Bányai K, Rahman M, et al. Rotavirus disease and vaccination: impact on genotype diversity. Future Microbiol. 2009;4(10):1303–16 Available from: https://academic.oup.com/jid/article/218/4/546/5001444.
World Health Organization. Rotavirus vaccines WHO position paper - January 2013. Wkly Epidemiol Rec. 2013;88(5):49–64 Available from: http://www.who.int/wer/2013/wer8805.pdf.
Mwenda JM, Parashar UD, Cohen AL, Tate JE. Impact of rotavirus vaccines in Sub-Saharan African countries. Vaccine. 2018; Available from: https://doi.org/10.1016/j.vaccine.2018.06.026.
De Vos B, Han HH, Bouckenooghe A, Debrus S, Gillard P, Ward R, et al. Live attenuated human rotavirus vaccine, RIX4414, provides clinical protection in infants against rotavirus strains with and without shared G and P genotypes: integrated analysis of randomized controlled trials. Pediatr Infect Dis J. 2009;28(4):261–6 Available from: https://journals.lww.com/pidj/Abstract/2009/04000/Live_Attenuated_Human_Rotavirus_Vaccine,_RIX4414,.1.aspx.
Wandera EA, Mohammad S, Bundi M, Nyangao J, Galata A, Kathiiko C, et al. Impact of rotavirus vaccination on rotavirus hospitalization rates among a resource-limited rural population in Mbita, Western Kenya. Trop Med Int Health. 2018;00(00):1–8 Available from: http://doi.wiley.com/10.1111/tmi.13040%0A. http://www.ncbi.nlm.nih.gov/pubmed/29432666.
Khagayi S, Omore R, Otieno GP, Ogwel B, Ochieng JB, Juma J, et al. Effectiveness of Monovalent Rotavirus Vaccine Against Hospitalization With Acute Rotavirus Gastroenteritis in Kenyan Children. Clin Infect Dis. 2019; Available from: https://academic.oup.com/cid/advance-article/doi/10.1093/cid/ciz664/5536640.
Otieno GP, Bottomley C, Khagayi S, Adetifa I, Ngama M, Omore R, et al. Impact of the introduction of rotavirus vaccine on hospital admissions for diarrhoea among children in Kenya: A controlled interrupted time series analysis. Clin Infect Dis [Internet]. 2019;(Cdc):1–8. Available from: https://academic.oup.com/cid/advance-article/doi/10.1093/cid/ciz912/5572563.
Tate JE, Ngabo F, Donnen P, Gatera M, Uwimana J, Rugambwa C, et al. Effectiveness of Pentavalent rotavirus vaccine under conditions of routine use in Rwanda. Clin Infect Dis [Internet]. 2016;62(suppl 2):S208–12 Available from: https://academic.oup.com/cid/article-lookup/doi/10.1093/cid/civ1016.
Bar-Zeev N, Kapanda L, Tate JE, Jere KC, Iturriza-Gomara M, Nakagomi O, et al. Effectiveness of a monovalent rotavirus vaccine in infants in Malawi after programmatic roll-out: an observational and case-control study. Lancet Infect Dis. 2015;15(4):422–8 Available from: https://www.thelancet.com/journals/laninf/article/PIIS1473-3099(14)71060-6/fulltext.
Armah G, Pringle K, Enweronu-Laryea CC, Ansong D, Mwenda JM, Diamenu SK, et al. Impact and effectiveness of monovalent rotavirus vaccine against severe rotavirus diarrhea in Ghana. Clin Infect Dis. 2016;62(April):S200–7 Available from: https://academic.oup.com/cid/article/62/suppl_2/S200/2478863.
Scott JAG, Bauni E, Moisi JC, Ojal J, Gatakaa H, Nyundo C, et al. Profile: the Kilifi health and demographic surveillance system (KHDSS). Int J Epidemiol. 2012;41(3):650–7 Available from: https://academic.oup.com/ije/article/41/3/650/837001.
Nokes DJ, Abwao J, Pamba A, Peenze I, Dewar J, Maghenda JK, et al. Incidence and clinical characteristics of group a rotavirus infections among children admitted to hospital in Kilifi, Kenya. PLoS Med. 2008;5(7):1154–62 Available from: https://journals.plos.org/plosmedicine/article?id=10.1371/journal.pmed.0050153.
Odhiambo FO, Laserson KF, Sewe M, Hamel MJ, Feikin DR, Adazu K, et al. Profile: the KEMRI/CDC health and demographic surveillance system-Western Kenya. Int J Epidemiol. 2012;41(4):977–87 Available from: https://academic.oup.com/ije/article/41/4/977/689203.
Feikin DR, Olack B, Bigogo GM, Audi A, Cosmas L, Aura B, et al. The burden of common infectious disease syndromes at the clinic and household level from population-based surveillance in rural and urban Kenya. PLoS One. 2011;6(1):1–10 Available from: https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0016085.
Gentsch JR, Glass RI, Woods P, Gouvea V, Gorziglia M, Flores J, et al. Identification of group a rotavirus gene 4 types by polymerase chain reaction. J Clin Microbiol. 1992;30(6):1365–73 Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC265294.
Gómara MI, Cubitt D, Desselberger U, Gray J. Amino acid substitution within the VP7 protein of G2 rotavirus strains associated with failure to serotype. J Clin Microbiol. 2001;39(10):3796–8 Available from: https://jcm.asm.org/content/39/10/3796.
Pickett BE, Greer DS, Zhang Y, Stewart L, Zhou L, Sun G, et al. Virus pathogen database and analysis resource (ViPR): a comprehensive bioinformatics database and analysis resource for the coronavirus research community. Viruses. 2012;4(11):3209–26 Available from: https://www.mdpi.com/1999-4915/4/11/3209.
John B O, Sammy K, Billy O, Reuben O, Evans A, Jane J, et al. Distribution of rotavirus genotypes among children with diarrhea prior to vaccine introduction in Western Kenya. J Infect Dis Epidemiol [Internet]. 2019 Mar 31;5(1). Available from: https://www.clinmedjournals.org/articles/jide/journal-of-infectious-diseases-and-epidemiology-jide-5-070.php?jid=jide.
Gouvea V, Glass RI, Woods P, Taniguchi K, Clark HF, Forrester B, et al. Polymerase chain reaction amplification and typing of rotavirus nucleic acid from stool specimens. J Clin Microbiol. 1990;28(2):276–82 Available from: https://jcm.asm.org/content/28/2/276.
Iturriza-Gómara M, Kang G, Gray J. Rotavirus genotyping: keeping up with an evolving population of human rotaviruses. J Clin Virol. 2004;31(4):259–65 Available from: https://www.sciencedirect.com/science/article/abs/pii/S1386653204001520?via%3Dihub.
Iturriza-gómara M, Green J, Brown DWG, Desselberger U, Gray JJ, Green JON, et al. Diversity within the VP4 gene of rotavirus P  strains : implications for reverse transcription-PCR genotyping. J Clin Microbiol. 2000;38(2):1–5 Available from: https://jcm.asm.org/content/38/2/898.
Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2017;20(4):1–7 Available from: http://academic.oup.com/bib/article/doi/10.1093/bib/bbx108/4106928/MAFFT-online-service-multiple-sequence-alignment.
Larsson A. AliView : a fast and lightweight alignment viewer and editor for large datasets. Oxford Univesity [Internet]. 2014;30(22):3276–8. Available from: https://academic.oup.com/bioinformatics/article/30/22/3276/2391211.
Nguyen L, Schmidt HA, Haeseler A Von, Minh BQ. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Oxford Univesity [Internet]. 2014;32(1):268–74. Available from: https://doi.org/10.1093/molbev/msu300.
Cavanaugh JE. Model selection: Bayesian information criterion. Wiley StatsRef Stat Ref Online. 2016;2(20):1–3 Available from: http://doi.wiley.com/10.1002/9781118445112.stat00247.pub2.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9. Available from. https://doi.org/10.1093/molbev/msy096.
Roczo-Farkas S, Kirkwood CD, Cowley D, Barnes GL, Bishop RF, Bogdanovic-Sakran N, et al. The impact of rotavirus vaccines on genotype diversity: a comprehensive analysis of 2 decades of Australian surveillance data. J Infect Dis. 2018;218(4):546–54 Available from: https://academic.oup.com/jid/article/218/4/546/5001444.
Zeller M, Rahman M, Heylen E, De Coster S, De Vos S, Arijs I, et al. Rotavirus incidence and genotype distribution before and after national rotavirus vaccine introduction in Belgium. Vaccine. 2010;28(47):7507–13. Available from. https://doi.org/10.1016/j.vaccine.2010.09.004.
Gelaw A, Pietsch C, Liebert UG. Molecular epidemiology of rotaviruses in Northwest Ethiopia after national vaccine introduction. Infect Genet Evol. 2018;65(June):300–7. Available from. https://doi.org/10.1016/j.meegid.2018.08.016.
Rahajamanana VL, Raboba JL, Rakotozanany A, Razafindraibe NJ, Andriatahirintsoa EJPR, Razafindrakoto AC, et al. Impact of rotavirus vaccine on all-cause diarrhea and rotavirus hospitalizations in Madagascar. Vaccine. 2018;36(47):7198–204. Available from. https://doi.org/10.1016/j.vaccine.2017.08.091.
Utsumi T, Wahyuni RM, Doan YH, Dinana Z, Soegijanto S, Fujii Y, et al. Equine-like G3 rotavirus strains as predominant strains among children in Indonesia in 2015–2016. Infect Genet Evol. 2018;61(June):224–8. Available from. https://doi.org/10.1016/j.meegid.2018.03.027.
Steele AD, Victor JC, Carey ME, Tate JE, Atherly DE, Pecenka C, et al. Experiences with rotavirus vaccines: can we improve rotavirus vaccine impact in developing countries? Hum Vaccin Immunother. 2019;5515:1–13. Available from. https://doi.org/10.1080/21645515.2018.1553593.
Dóró R, Marton S, Bartókné AH, Lengyel G, Agócs Z, Jakab F, et al. Equine-like G3 rotavirus in Hungary, 2015 -is it a novel intergenogroup reassortant pandemic strain? Acta Microbiol Immunol Hung. 2016;63(2):243–55. Available from. https://doi.org/10.1556/030.63.2016.2.8.
Matthijnssens J, Ciarlet M, Heiman E, Arijs I, Delbeke T, McDonald SM, et al. Full genome-based classification of rotaviruses reveals a common origin between human Wa-like and porcine rotavirus strains and human DS-1-like and bovine rotavirus strains. J Virol. 2008;82(7):3204–19 Available from: http://jvi.asm.org/cgi/doi/10.1128/JVI.02257-07.
Yen C, Figueroa JR, Uribe ES, del Carmen-Hernández L, Tate JE, Parashar UD, et al. Monovalent rotavirus vaccine provides protection against an emerging fully heterotypic G9P rotavirus strain in Mexico. J Infect Dis. 2011;204(5):783–6 Available from: https://academic.oup.com/jid/article-lookup/doi/10.1093/infdis/jir390.
Steele AD, Neuzil KM, Cunliffe NA, Madhi SA, Bos P, Ngwira B, et al. Human rotavirus vaccine Rotarix™ provides protection against diverse circulating rotavirus strains in African infants: a randomized controlled trial. BMC Infect Dis. 2012;12(1):213 Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3462149&tool=pmcentrez&rendertype=abstract.
Matthijnssens J, Potgieter CA, Ciarlet M, Parreno V, Martella V, Banyai K, et al. Are human P  rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla? J Virol. 2009;83(7):2917–29 Available from: http://jvi.asm.org/cgi/doi/10.1128/JVI.02246-08.
Carvalho-Costa FA, Araújo IT, Santos De Assis RM, Fialho AM, De Assis Martins CMM, Bóia MN, et al. Rotavirus genotype distribution after vaccine introduction, Rio de Janeiro, Brazil. Emerg Infect Dis. 2009;15(1):95–7 Available from: https://wwwnc.cdc.gov/eid/article/15/1/07-1136_article.
Mukaratirwa A, Berejena C, Nziramasanga P, Ticklay I, Gonah A, Nathoo K, et al. Distribution of rotavirus genotypes associated with acute diarrhoea in Zimbabwean children less than five years old before and after rotavirus vaccine introduction. Vaccine. 2018; Available from. https://doi.org/10.1016/j.vaccine.2018.03.069.
Kaneko M, Takanashi S, Thongprachum A, Hanaoka N, Fujimoto T, Nagasawa K, et al. Identification of vaccine-derived rotavirus strains in children with acute gastroenteritis in Japan, 2012-2015. PLoS One. 2017;12(9):2012–5. Available from. https://doi.org/10.1371/journal.pone.0184067.
Payne DC, Edwards KM, Bowen MD, Keckley E, Peters J, Esona MD, et al. Sibling transmission of vaccine-derived rotavirus (RotaTeq) associated with rotavirus gastroenteritis. Pediatrics. 2010;125(2) Available from. https://doi.org/10.1542/peds.2009-1901.
Bennett A, Pollock L, Jere KC, Pitzer VE, Lopman B, Parashar U, et al. Infrequent transmission of monovalent human rotavirus vaccine virus to household contacts of vaccinated infants in Malawi. J Infect Dis. 2019;219(11):1730–4. Available from. https://doi.org/10.1093/infdis/jiz002.
The authors are indebted to the participants who provided samples for analysis. We thank Laboratory staff at Virus Epidemiology and Control group, KEMRI Wellcome Trust Research Programme and KEMRI Center for Global Health Research for processing samples and initial analysis.
The findings and conclusions in this report are those of the authors and do not necessarily represent the views of involved affiliated institutions.
This work was funded by Global Alliance for Vaccines (GAVI) to the Emory University Global Health Institute, which supported the Rotavirus Immunization Program and Evaluation in Kenya (RIPEK). The funders had no role in design of the study, collection, analysis and interpretation of the data and in writing of the manuscript.
Ethics approval and consent to participate
The study protocol was submitted for review and approval by institutional ethics review boards of Kenya Medical Research Institute (Scientific Ethics Review Unit protocol no. #3049) and United States Centers for Disease Control and Prevention (CDC protocol no. #6968). CDC’s human subjects research office relied on KEMRI for IRB oversight. Written informed consent for enrolment into the study was obtained from parents or guardians of children admitted with diarrhea symptoms. Permission to access and use the data was provided by the KEMRI Wellcome Trust Research Programme data governance committee.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
; Geographical boundaries and location of the health facilities participating in the rotavirus genotype surveillance program represented by the diamond shapes. Green – Siaya County Referral Hospital, blue – Lwak Mission Hospital, red – Tabitha Clinic and orange – Kilifi County Hospital.
Frequency of partially typed G and P genotypes. Gx and P [x] were unsuccessfully typed due to failure in sequencing and/or contig assembly.
GenBank accession numbers of all VP7 G gene sequences.
GenBank accession numbers of all VP4 P gene sequences.
About this article
Cite this article
Mwanga, M.J., Owor, B.E., Ochieng, J.B. et al. Rotavirus group A genotype circulation patterns across Kenya before and after nationwide vaccine introduction, 2010–2018. BMC Infect Dis 20, 504 (2020). https://doi.org/10.1186/s12879-020-05230-0