Skip to main content

Differences in epidemiology of enteropathogens in children pre- and post-rotavirus vaccine introduction in Kilifi, coastal Kenya



Kenya introduced Rotarix® (GlaxoSmithKline Biologicals, Rixensart, Belgium) vaccination into its national immunization programme beginning July 2014. The impact of this vaccination program on the local epidemiology of various known enteropathogens is not fully understood.


We used a custom TaqMan Array Card (TAC) to screen for 28 different enteropathogens in 718 stools from children aged less than 13 years admitted to Kilifi County Hospital, coastal Kenya, following presentation with diarrhea in 2013 (before vaccine introduction) and in 2016–2018 (after vaccine introduction). Pathogen positivity rate differences between pre- and post-Rotarix® vaccination introduction were examined using both univariate and multivariable logistic regression models.


In 665 specimens (92.6%), one or more enteropathogen was detected, while in 323 specimens (48.6%) three or more enteropathogens were detected. The top six detected enteropathogens were: enteroaggregative Escherichia coli (EAggEC; 42.1%), enteropathogenic Escherichia coli (EPEC; 30.2%), enterovirus (26.9%), rotavirus group A (RVA; 24.8%), parechovirus (16.6%) and norovirus GI/GII (14.4%). Post-rotavirus vaccine introduction, there was a significant increase in the proportion of samples testing positive for EAggEC (35.7% vs. 45.3%, p = 0.014), cytomegalovirus (4.2% vs. 9.9%, p = 0.008), Vibrio cholerae (0.0% vs. 2.3%, p = 0.019), Strongyloides species (0.8% vs. 3.6%, p = 0.048) and Dientamoeba fragilis (2.1% vs. 7.8%, p = 0.004). Although not reaching statistical significance, the positivity rate of adenovirus 40/41 (5.8% vs. 7.3%, p = 0.444), norovirus GI/GII (11.2% vs. 15.9%, p = 0.089), Shigella species (8.7% vs. 13.0%, p = 0.092) and Cryptosporidium spp. (11.6% vs. 14.7%, p = 0.261) appeared to increase post-vaccine introduction. Conversely, the positivity rate of sapovirus decreased significantly post-vaccine introduction (7.8% vs. 4.0%, p = 0.030) while that of RVA appeared not to change (27.4% vs. 23.5%, p = 0.253). More enteropathogen coinfections were detected per child post-vaccine introduction compared to before (mean: 2.7 vs. 2.3; p = 0.0025).


In this rural Coastal Kenya setting, childhood enteropathogen infection burden was high both pre- and post-rotavirus vaccination introduction. Children who had diarrheal admissions post-vaccination showed an increase in coinfections and changes in specific enteropathogen positivity rates. This study highlights the utility of multipathogen detection platforms such as TAC in understanding etiology of childhood acute gastroenteritis in resource-limited regions.


In 2016, there were approximately 446,000 deaths worldwide caused by diarrheal illnesses among children aged < 5 years, with the majority occurring in low-income countries [1, 2]. Rotavirus group A (RVA) was estimated to be responsible for approximately 128,000 of these deaths.

The World Health Organization (WHO) recommended inclusion of RVA vaccines into national immunization programs (NIPs) of all countries in 2009 [3]. As of November 2021, 114 countries had included rotavirus vaccination into their NIPs [4]. Kenya included Rotarix® (GlaxoSmithKline Biologicals, Rixensart, Belgium), one of the WHO pre-qualified rotavirus vaccines, in its NIP in July 2014 [3], with doses given at 6 and 10 weeks of age. Based on antigen testing (enzyme immunoassay; EIA), Kilifi County Hospital (KCH), Coastal Kenya, reported a 57% reduction in rotavirus hospitalizations in the first year after vaccine introduction and an 80% reduction in the second year among children < 5-year-olds [4]. Further, KCH found a 64% Rotarix® vaccine effectiveness in children aged < 5 years [5], which is similar to other low-income settings in sub-Saharan Africa [5, 6], but lower than high-income settings [7, 8].

Over 30 potential enteropathogens can cause diarrhea. In addition to a reduction in RVA-induced diarrhea, a referral peri-urban hospital in Central Kenya observed a 31% decrease in all-cause diarrheal hospital admissions in the first year and a 58% decrease in the second year [9]. While RVA has clearly been responsible for a large burden of severe diarrheal cases among children < 5 years-old, there are many other enteropathogens that contribute to severe diarrheal illness in Kenya [10,11,12,13] but their epidemiology, both pre- and post-rotavirus vaccine introduction, is under-studied. Previously, we examined change in prevalence of 5 enteric viruses, including RVA, pre- and post-rotavirus vaccination demonstrating an increase in norovirus GII and a decrease in RVA, however this study did not examine many other potential infectious causes of diarrheal illness [14].

Using a custom TaqMan Array Card (TAC), this study aimed to compare the epidemiological patterns of 28 different enteropathogens in children < 13 years-old admitted to KCH with diarrhea, pre- and post-rotavirus vaccine introduction.


Study design and population

KCH is a referral facility that primarily serves residents of Kilifi County located on the Kenyan Coast [15]. From September 2009, a surveillance study of rotavirus was established at KCH pediatric ward targeting children < 13 years old who presented with diarrhea (defined as passing three or more loose stools in a 24-h period with or without visible blood) [16]. At KCH, pediatric admission is determined by the clinician’s assessment of the severity of the illness and children must stay overnight once admitted. A single stool specimen was collected within 48 h of admission and immediately put into − 4 °C before transferring to the adjoining KEMRI-Wellcome Trust Research Programme (KWTRP) laboratories for long-term storage at −  80 °C. The current retrospective analysis targeted samples collected in the year 2013 (pre-vaccine period) and 2016–2018 (post-vaccine period). During these periods, sample collection was interrupted 6 times by healthcare worker strikes; 3 times each pre- and post-vaccine introduction [17].

Sample processing

Total nucleic acid (TNA) was extracted from stool using two different approaches. First, the cador Pathogen 96 QIAcube HT Kit (QIAGEN, Hilden, Germany) with 200 µL liquid or 200 mg solid stool was used on samples from 2013. Second, the QIAamp Fast DNA Stool Mini Kit (QIAGEN, Hilden, Germany) with 200 µL or 200 mg stool after bead-beating and centrifugation to remove debris was used on samples from 2016 to 2018. An extraction blank was included in every processing batch for quality control. TNA extracts were stored at − 80 °C for up to 130 weeks and defrosted once, only for the enteropathogen diagnostics.

TaqMan array card (TAC) analysis

A custom TAC spotted with lyophilized target-specific primers and probes was used to detect enteropathogens from the TNA extracts. The GASTRO v4.0K TAC was designed by MC at Public Health England (PHE), Cambridge, England [18] and manufactured by ThermoFisher, USA. GASTRO v4.0K screened for 28 enteropathogens: 12 viruses, 10 bacteria, 5 protozoa and 1 helminth (Fig. 1 and Additional file 1: Table S1). All TAC targets were validated previously by PHE and include 16S bacterial and 18S RNA as internal controls, and bacteriophage MS2 as an external control. A Rotarix® vaccine-specific assay targeting the NSP2 gene was included to identify RVA vaccine shedding due to recent vaccination. The molecular probes were labeled at the 5’ end with 6-carboxyfluorescein (FAM) reporter dye and NFQ-MGB quencher dye at the 3’ end [18].

Fig. 1
figure 1

The configuration of the GASTRO v4.0 K TaqMan Array Card used for screening 28 enteropathogens. Some of the enteropathogens had multiple spots (targets) on the card

For each TAC, 8 samples were tested (48-wells per sample; Fig. 1). A mix of 60 µL of nuclease-free water, 25 µL of TaqMan Fast Virus 1-Step MasterMix (ThermoFisher, USA) containing ROX reference dye and 15 µL of TNA was prepared and loaded onto the TAC sample portal. TACs were spun twice at 1200 rpm for two minutes and sealed. TACs were processed using the Quantistudio-7-Flex Real-time PCR System (Thermofisher, USA). The thermocycling conditions were 50 °C for 5 min followed by 45 cycles of 95 °C for 1 s and 60 °C for 20 s. Extraction blanks and no template controls were included once in every 10 runs for quality control.

Amplification data were analyzed by the cycle threshold (Ct) method (QuantiStudio™Real-Time PCR software v.1.1). A threshold value of 0.2 florescence units was applied to all our analyses and baseline range set to automatic. The amplification curves for presumptive positive samples were all visually inspected for confirmation. Pathogen prevalence was evaluated at Ct cut-off levels; <40.0, < 35.0 and < 30.0.

Clinical and vaccination data

Clinical data were collected using a standard clinical research form (CRF) entered into an electronic database [19]. The CRF included admission date, discharge date, presenting signs and symptoms, co-morbidities, and discharge diagnosis and outcome. Diarrhea was categorized as acute (< 7 days), prolonged (7–13 days), persistent (14–28 days), and chronic (> 28 days). Diarrheal severity was assessed by Vesikari Clinical Severity Score (VCSS) [20]. All samples were tested for RVA using a commercial EIA antigen detection kit (ProsPeCT, Oxoid, Thermo Scientific). All pediatric admissions were offered HIV testing using two rapid antibody tests according to national guidelines [21].

Rotavirus vaccination status for children was obtained from the Kilifi electronic vaccine registry [22]. The registry captures the child’s vaccination data real-time in clinics within the Kilifi health and demographic surveillance system (KHDSS) operated by KWTRP. This is then linked with KCH surveillance data using unique KHDSS identification numbers.


We compared clinical and demographic characteristics of (i) included and excluded children and (ii) those recruited pre- and post-rotavirus vaccination using chi-squared analysis and Fisher’s exact test for categorical variables and Wilcoxon-rank sum and t-tests for continuous variables. As this was the first use of this TAC in Kenya, we descriptively explored different lower limits of detection (Ct ≤ 30, 35, 40) for a positive result. Using Ct ≤ 35 as the lower limit of detection, pathogen prevalence was calculated as the proportion positive over the total tested with 95% confidence intervals (CI) calculated using the exact method. We assessed differences in pathogen presence in pre-post-rotavirus vaccination periods by univariate logistic regression. Any enteropathogen that was marginally different (p < 0.1) in univariate analyses, as well as rotavirus (due to epidemiological interest), were included in multivariable logistic regression analyses adjusted for age, sex, Vesikari score, HIV-status, hospital length of stay (LoS), and fatality. Separate multivariable logistic regression analyses per pathogen were conducted due to the high prevalence of co-infections. The analyses were performed in Stata version 16.1 and R version 3.6.3.

Ethical considerations

Before samplecollection, informed written consent was gained from each child’s parent/guardian. The Scientific and Ethics Review Unit (SERU) at Kenya Medical Research Institute, Nairobi, approved the study protocol (Protocol #3624).


Participant characteristics

A total of 1317 KCH pediatric admissions with diarrhea were eligible for inclusion during study periods (2013, 2016–2018). Of these, 1250 (95%) were unique patients; 718 (57.4%) provided an adequate stool sample for TAC processing. Children who did and did not have a stool sample were similar except those with a stool sample: (i) were more likely to have reported vomiting at admission (p = 0.001); (ii) were hospitalized for longer (mean 5.8 days versus 4.3 days, p < 0.001); (iii) were more likely to have been discharged alive (p < 0.001); (iv) were more likely to have or have been exposed to HIV (p = 0.049); and (v) were less likely to be included in 2017 and 2018 than in 2013 and 2016 (p < 0.001) (Table 1).

Table 1 Demographic and clinical characteristics of the eligible, included and excluded participants and, among participants, pre-vaccine versus the post-vaccine introduction period comparisons

The children who had provided sufficient stool had a median age of 13.3 months (IQR: 8.0–23.0). Clinical and demographic characteristics of children from pre- and post-vaccine introduction periods (n = 718) were similar except that those from the pre-vaccine period were (i) significantly younger (11.5 vs. 14.1 months, p = 0.004); (ii) less likely to have known HIV exposure, but more likely to have unknown HIV status (p < 0.001); and (iii) had a lower median hospital stay (4 vs. 5 days, p = 0.005) than those from the post-vaccine period (Table 1).

Pathogen positivity rate across different ct cut-offs

Of the 718 specimens analyzed by TAC, using a cycle threshold cut-off of ≤ 40.0 Ct (which was considered an appropriate cut-off based on clinical samples from the UK [18, 23]) one or more enteropathogen was detected in 94.0% of the children (n = 675). Overall, the top five enteropathogens detected at the ≤ 40.0 Ct cut off value were: enteroaggregative Escherichia coli (EAggEC) 44.6%; enterovirus 33.0%; enteropathogenic Escherichia coli (EPEC) 32.3%; RVA 25.1%; and parechovirus 20.8% (Fig. 2A). When the limit of detection was lowered to ≤ 30 Ct (high pathogen burden), 12/26 of detected pathogens decreased in prevalence by a factor of 2 or more (Fig. 2B). The largest declines were observed with verotoxogenic E. coli (VTEC) (11 fold), Dientamoeba fragilis (8.3 fold), cytomegalovirus (7.4 fold), Strongyloides spp. (5.0 fold), Salmonella spp. (4.0 fold), Clostridium perfringens/difficile (3.8 fold) and hepatitis E virus (3.0 fold). Only three targets (D. fragilis, hepatitis E virus, and VTEC) showed a reduction by factor 2 or more when comparing a Ct value of ≤ 40 versus ≤ 35 as the limit of detection. On average, the prevalence of the detected pathogens decreased by 1.3 fold when the Ct cut-off was lowered from 40.0 to 35.0 (Fig. 2C). Notably, even with the most stringent Ct cut-off (≤ 30), the prevalence of RVA was unaffected: 25.1% (21.9–28.4) at Ct ≤ 40 versus 22.3% (19.3–25.5) at Ct ≤ 30. No samples tested positive for Vibrio parahaemolyticus, Entamoeba histolytica or the Rotarix® vaccine strain at any Ct detection limit. All subsequent analyses have used the Ct cut-off of ≤ 35.0 for all targets. Detection above this threshold has previously been considered clinically insignificant or irreproducible in studies from low-income settings [24, 25].

Fig. 2
figure 2

Prevalence of enteropathogens in stool samples from 718 children hospitalised at KCH at different PCR cycle threshold (Ct) cut-offs. The bars represent organism type and include 95% confidence interval errors bars for the proportions. Panel (a) prevalence of the detected enteropathogens when applying a Ct cut-off of ≤ 40 to define positives. Panel (b) prevalence of the detected enteropathogens when applying a Ct cut-off ≤ 35 to define positives. Panel (c) prevalence of detected the enteropathogens when applying a Ct cut-off ≤ 30 to define positives

Pathogen positivity rate pre-post-rotavirus vaccine introduction

Among the 718 samples, 241 (33.6%) were from the pre-vaccine period (2013) and 477 (66.4%) were from the post-vaccine introduction period (2016–2018; Table 2). In univariate logistic regression analyses a significant increase in positivity rate was observed for cytomegalovirus (4.2% vs. 9.9%, OR = 2.53, p = 0.010), EAggEC (35.7% vs. 45.3%, OR = 1.49, p = 0.014), D. fragilis (2.0% vs. 7.8%, OR = 3.97, p = 0.004) and Strongyloides spp. (0.8% vs. 3.6%, OR = 4,42, p = 0.048). Also, an increase in the positivity rate was observed for adenovirus 40/41 (5.8% vs. 7.3%, OR = 1.28, p = 0.444), norovirus GI/GII (11.2% vs. 15.9%, OR = 1.5, p = 0.089), Shigella spp. (8.7% vs. 13.0%, OR = 1.57, p = 0.092) and Cryptosporidium spp. (11.6% vs. 14.7%, OR = 1.31, p = 0.261) post-vaccine introduction but this did not reach statistical significance (Table 2). A significant decrease in test positivity rate was observed for sapovirus (7.8% vs. 4.0%, OR = 0.49, p = 0.030).

Table 2 Univariate and multivariate logistic regression comparing enteropathogen infection in the pre- vs post-rotavirus vaccine introduction periods (n = 718)

In multivariable models adjusting for age, Vesikari score, fatality, gender, LOS, and HIV-status, cytomegalovirus (OR = 2.81, 95% CI 1.36, 5.81), EAggEC (OR = 1.58, 95% CI 1.13, 2.20), D. fragilis (OR = 3.90, 95% CI 1.49, 10.21), and Strongyloides spp. (OR = 4.60, 95% CI 1.03, 20.55) were all detected among significantly more patients in the post-vaccine than in the pre-vaccine period. Sapovirus no longer demonstrated a significant difference between the two time periods (Table 2).

No differences were observed for RVA pre-post-RVA vaccination (27.4% vs. 23.5%, p = 0.253). Under the multivariable models the risk of RVA positive status did not change after adjusting for age and receipt of recent RVA vaccination (data not shown). The absence of reduction in RVA prevalence and absence of association with vaccine status remained even when we restricted our analyses to children < 5 years. No significant changes were observed for the other remaining pathogens.

Enteropathogen co-infections

Only 53 children (7.4%) had not even a single enteropathogen detected. Among the rest (n = 665), 148 (22.3%) were mono-infected, while 323 (48.6%) had ≥ 3 pathogens detected. Including all children, 72% (n = 517) had ≥ 2 enteropathogens detected. Five or more enteropathogens were detected in 84 children (11.7%). A greater number of different enteropathogens were detected per child in the post-vaccine compared to pre-vaccine period (mean: 2.7 vs. 2.3 respectively, p = 0.0025).

Characteristics of rotavirus positive samples

The distribution of RVA TAC assay Ct values for positive samples stratified by various potential influencing factors is shown in Fig. 3. The median Ct value was statistically similar between (a) RVA EIA test positive and negative samples (p = 0.214; Fig. 3A), (b) TAC positive samples pre- and post-vaccine introduction (Fig. 3B; p = 0.098), (c) across individuals who had received various numbers of rotavirus vaccine doses including the unvaccinated group (Fig. 3C; p = 0.197). However we noted that the median Ct values among older children were significantly lower than those among younger children, regardless of the sampling year (Fig. 3D; p < 0.001).

Fig. 3
figure 3

Boxplots examining the relationship between the observed RVA detection cycle threshold (Ct) from TAC assay and (a) previous enzyme-linked immunosorbent assay (ELISA) result on the same samples, (b) period of sample collection, (c) rotavirus vaccination status, and (d) age group in stool samples from 718 children at KCH. P-values are derived from an independent t-test, one-way ANOVA, and where applicable Kruskal Wallis was used


Our study shows very high positivity rate of enteropathogens among hospitalized children with diarrhea, extensive co-infections, and differences in pathogen positivity before and after introduction of rotavirus vaccination in this rural coastal setting of Kenya. Further, despite introduction of childhood rotavirus vaccination, RVA was the fourth most detected and we did not observe a significant difference in RVA positivity rate between pre- and post-vaccination introduction periods following TAC analysis. There have been few studies focused on illuminating on changes in multiple enteropathogen presentation post-rotavirus vaccination [26,27,28,29,30,31]. Our findings here align with previous research in low-middle income countries (LMIC) settings, such as Global Enteric Multicenter Study (GEMS) [25] and Malnutrition and Enteric Disease Study (MAL-ED) study [32], which also found high rates of enteropathogen co-infections. The changes in pathogen positivity rate pre-post rotavirus vaccination and increase in co-infection rate suggests a changing landscape of diarrhea etiology and epidemiology in Kenya. We hypothesize that the demographics of hospitalized children with severe diarrhea after rotavirus vaccination is changing to somewhat older children and those with additional risk factors such as malnutrition and HIV exposure.

The top nine enteropathogens detected in this study (in descending order) were EAggEC, EPEC, enterovirus, RVA, parechovirus, norovirus GI/GII, Cryptosporidium spp., Shigella spp. and Giardia lambia. A study of bacterial diarrheal infections in < 5 year-olds in Western Kenya revealed that E. coli and Shigella species were the most common bacterial infections [33]. Our previous study examining 5 viral enteropathogens among children < 5 years at KCH found RVA, and adenovirus as the commonest enteric viruses [14]. In the GEMS study which included children < 5 years-old and community controls, the top nine enteropathogens associated with diarrheal illness were RVA, enterotoxigenic E. coli, Cryptosporidium spp., H. pylori, norovirus GII, Aeromonas spp., adenovirus 40/41, EPEC, Shigella spp., C. difficile, and sapovirus [34]. A post-vaccine introduction study of hospitalized children aged < 5 years in Tanzania found RVA, ETEC, Shigella spp., Cryptosporidium spp. and astrovirus as the leading causes of diarrhea [26]. Although all the above studies took place in different geographic regions and/or used different methods of pathogen detection from the current study, there is still considerable consistency in the major enteropathogens observed.

While the enteropathogens detected in our study could cause diarrhea, several of them have been associated with asymptomatic shedding or non-diarrheic syndromes [35]. EAggEC was the most common enteropathogen detected here but several previous studies e.g., one in Malawi found EAggEC frequently in both cases (51.8%) and controls (47.8%) [28]. Generally, the role of EAggEC in acute diarrhea in sub-Saharan Africa settings remains uncertain but prolonged carriage has been associated with child growth faltering [36]. EPEC, the second most prevalent pathogen we detected here has been linked with diarrhea, for example in a Malawi study, where it was found in 18.0% of cases versus 8.3% of controls [28]. Enterovirus, the third top enteropathogen we detected has been previously associated with a wide spectrum of clinically distinct syndromes including respiratory and neurological illnesses [37]. Thus, some of the detections here could be long-term carriage unassociated with the current diarrhea episode. To better understand changes in enteropathogen epidemiology pre-post-rotavirus vaccine introduction, repeated measures studies on the same individuals including community controls are necessary [25, 32]. Such studies will also be important in understanding the longer-term impact of vaccination, carriage, and infection on children in LMIC settings.

RVA positivity rate in this study post vaccine introduction was 23.5%. RVA has been reported in several studies in sub-Sarahan African settings be the leading cause severe acute gastroenteritis (AGE) even post-vaccine introduction [26, 28]. In Tanzania, RVA was found to be a leading diarrhea hospitalization (attributable fraction (AF), 25.8%) [26] and in Malawi (AF, 34.5%) [28]. Nonetheless, the Global Rotavirus Surveillance Network showed that RVA prevalence decreased in < 5-year-olds between 2008 and 2016 following national introduction of the RVA vaccine across 82 countries, although to a lesser extent among low-income countries [24]. In Kenya, using controlled interrupted time series analysis, a > 50% decline in RVA cases were observed in two sites following Rotarix® introduction [38]. One of these sites was KCH. Further in a previous analysis, using RT-PCR, we also observed a reduction in RVA among children < 5-year-olds (23.3% vs. 13.8%) [14]. The current study found no reduction in RVA prevalence among KCH children; however, there were some key differences from the previous studies. First, to detect RVA infection, we used TAC, which is considerably more sensitive than the conventional RT-PCR [39] or EIA used in other studies [14, 40]. Second, we included children up to 13 years of age, whereas the other studies included only < 5-year-olds. Third, previous Kenyan studies, samples were collected between 2014 and 2017 [38] and 2003, 2013, 2016, and 2019 [14], which does not overlap entirely with the current study. The inclusion of older children and the use of TAC provide evidence of the overall prevalence of RVA in children, as opposed to only the vaccine eligible children [41].

Apart from the methodological differences with previous studies, the absence of a significant decline of RVA positivity post-vaccination observed here could be further due to local circulation of distinct RVA genotypes for which Rotarix ® has limited efficacy. We have previously reported a temporal increase in prevalence of Rotarix ® vaccine heterotypic genotypes in Kenya post-vaccination (e.g. G2P[4] and G3P[8]) [44]. It has also been reported elsewhere that RVA prevalence increased with age following introduction of vaccines [45], which might have occurred in this study that included children up to 13-years of age. Indeed overall, children from the post-rotavirus vaccine introduction period were older compared to pre-vaccine.

Studies have shown that RVA load in stool samples during an episode, as determined by RT-PCR Ct values, is significantly correlated with disease severity [46, 47]. It was found that children with more severe diarrhea excreted more virus, as evidenced by lower RT-PCR Ct values. We found that older children had significantly lower median Ct values compared with younger children. There may be a link between RVA being more severe in older children and either the non-vaccination status or the waning immune response after vaccination [48, 49]. Similar Ct values were seen in samples that were TAC positive but EIA negative, as well as in samples that were pre-vaccine versus post-vaccine. It is necessary to investigate why some TAC positives were missed by EIA to determine if there is a sensitivity issue with available EIA kits post-vaccine introduction.

The reanalysis of GEMS case-control study using TAC found Shigella spp./EIEC to be the topmost attributable pathogen for moderate-to-severe diarrhea [25]. Here, we observed an increase in prevalence of Shigella spp./EIEC pre-post rotavirus vaccine introduction although this was not statistically significant; 8.7% (95% CI 5.5–13.0%) versus 13.0 (95% CI 10.2–16.4%). The prevalence observed in this study was close to that observed in a study in neighboring Tanzania post-vaccine introduction study that found a Shigella spp./EIEC prevalence of 14.5% (95% CI 10.2–22.8%) [26].

Norovirus has been increasingly reported as a most important cause of sporadic and epidemic AGE post-rotavirus vaccine introduction [30, 50]. Although norovirus GI/GII was sixth as the most frequently detected enteropathogen in our study, its contribution to diarrhea etiology is likely to be significant given the suspected carriage of some of the top detected enteropathogens. Congruent with our previous findings [14], we observed its trend towards increased positivity rate pre-post-rotavirus vaccination introduction (11.2% vs. 15.9%) although this did not reach statistical significance.

We detected 11 V. cholerae positive cases, all occurring post-rotavirus vaccine introduction. This finding is consistent with several small outbreaks that were reported in Kenya during the period including in refugee camps [51]. The sources of this ongoing V. cholerae transmission in Kenya e.g. whether is a result of sustained cryptic transmission or repeated introductions into the country is unclear. In contrast to other studies in other low-resource settings, no samples in our study tested positive for Vibrio parahaemolyticus or Entamoeba histolytica [26, 32]. This was possibly due to different epidemiology for these pathogens in our setting or failure of the molecular assay we used due to primer/probe sequence mismatch [52].

This study had some limitations. Samples were collected more than two years prior to analysis, which could compromise their integrity [53]. The study setting might have experienced a further change in enteropathogen burden since 2018 in part due to the coronavirus disease 2019 (COVID-19) pandemic as has been observed elsewhere [54, 55]. However, this does not invalidate the hypotheses generated by the findings regarding enteropathogen dynamics, but rather suggests that continuous surveillance is critical for adapting and maximizing diarrheal illness prevention. We did not analyze healthy controls to adjust for the background prevalence of the detected enteropathogens in this population. Of eligible children, 56.9% provided samples, which limited our ability to assess the entire population of hospitalized children with diarrhea. As a result of health worker industrial action in Kenya, enrollment was interrupted five times [17]. In the pre-vaccine and post-vaccine introduction years, the total nucleic acid extraction procedure varied, but we did not expect this to have significantly impacted enteropathogen prevalence. Despite previous studies finding an association between pathogen quantities and diarrhea attribution, we did not attempt to model these associations in this study [25]. Given the unbalanced time series nature of this study it is difficult to draw conclusions about enteropathogen dynamics following rotavirus vaccine introduction. We reported pathogen positivity rate for study eligible children who provided a sample, but more robust analysis would require a population-based denominator to calculate population-wide incidence changes.


We provide a snapshot of enteropathogens infection burden among children with diarrheal illness in a low-income setting of coastal Kenya pre-post rotavirus vaccine introduction. Our findings highlight the utility of multiplex pathogen detection arrays in enteropathogen infection surveillance to improve understanding on the causes of differential global rotavirus vaccine efficacy. Our study revealed that most children presenting with diarrheal illness in this setting were coinfected with multiple enteropathogens emphasizing the need for local improved available enteropathogen vaccine coverage and emphasis of non-pharmaceutical interventions such as WASH (water, sanitation, and hygiene) to reduce carriage and transmission of these enteropathogens [56]. Overall, the TAC platform provides a fast, broad, and accurate diagnostic tool for broad routine surveillance of enteropathogens and quick screening of outbreak samples to identify potential etiological pathogens [57]. This platform has potential to offer insight into the patient infection profiles in a clinical setting to support high-quality interventions and individual patient care. The clinical implications of the frequent coinfections and if this should influence patient care during and post admission requires follow-up investigation.

Data availability

The data described in this manuscript can be accessed by via the link Further information can be accessed submitting a request form to our Data Governance Committee (


  1. Troeger C, Blacker BF, Khalil IA, et al. 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. The Lancet Infectious Diseases 2018; 18:1211–28.

    Article  Google Scholar 

  2. Rogawski ET, Liu J, Platts-Mills JA, et al. Use of quantitative molecular diagnostic methods to investigate the effect of enteropathogen infections on linear growth in children in low-resource settings: longitudinal analysis of results from the MAL-ED cohort study. Lancet Glob Health 2018; 6:e1319-e28.

    Article  Google Scholar 

  3. WHO. Rotavirus vaccines: WHO position paper- July 2021. Weekly Epidemiological Record. 2021;96:301–219.

    Google Scholar 

  4. Varghese T, Kang G, Steele AD. Understanding rotavirus vaccine efficacy and effectiveness in countries with high child mortality. Vaccines (Basel). 2022;10:346.

    Article  CAS  Google Scholar 

  5. Steele AD, Victor JC, Carey ME, et al. Experiences with rotavirus vaccines: can we improve rotavirus vaccine impact in developing countries? Human Vaccines & Immunotherapeutics 2019; 15:1215–27.

    Article  CAS  Google Scholar 

  6. Murunga N, G PO, Maia M, C NA. Effectiveness of Rotarix (®) vaccine in Africa in the first decade of progressive introduction, 2009–2019: systematic review and meta-analysis. Wellcome Open Res 2020; 5:187.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Willame C, Vonk Noordegraaf-Schouten M, Gvozdenovic E, et al. Effectiveness of the oral human attenuated rotavirus vaccine: a systematic review and meta-analysis—2006–2016. Open Forum Infect Dis 2018; 5:ofy292.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Walker JL, Andrews NJ, Atchison CJ, et al. Effectiveness of oral rotavirus vaccination in England against rotavirus-confirmed and all-cause acute gastroenteritis. Vaccine X 2019; 1:100005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Wandera EA, Mohammad S, Bundi M, et al. Impact of rotavirus vaccination on rotavirus and all-cause gastroenteritis in peri-urban Kenyan children. Vaccine 2017; 35:5217–23.

    Article  CAS  PubMed  Google Scholar 

  10. Omore R, Khagayi S, Ogwel B, et al. Rates of hospitalization and death for all-cause and rotavirus acute gastroenteritis before rotavirus vaccine introduction in Kenya, 2010–2013. BMC Infectious Diseases 2019; 19:47.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Swierczewski BE, Odundo EA, Koech MC, et al. Surveillance for enteric pathogens in a case–control study of acute diarrhea in Western Kenya. Trans R Soc Trop Med Hyg 2013; 107:83–90.

    Article  CAS  PubMed  Google Scholar 

  12. Njuguna C, Njeru I, Mgamb E, et al. Enteric pathogens and factors associated with acute bloody diarrhoea, Kenya. BMC Infect Dis 2016; 16:477.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Florez ID, Niño-Serna LF, Beltrán-Arroyave CP. Acute infectious diarrhea and gastroenteritis in children. Current Infectious Disease Rep. 2020;22:4.

    Article  PubMed  Google Scholar 

  14. Lambisia AW, Onchaga S, Murunga N, Lewa CS, Nyanjom SG, Agoti CN. Trends of five common diarrhea-associated enteric viruses pre- and post-rotavirus vaccine introduction in Coastal Kenya. Pathogens. 2020;9:660.

    Article  PubMed Central  Google Scholar 

  15. Scott JA, Bauni E, Moisi JC, et al. Profile: The Kilifi Health and Demographic Surveillance System (KHDSS). Int J Epidemiol 2012; 41:650–7.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Khagayi S, Omore R, Otieno GP, et al. Effectiveness of monovalent rotavirus vaccine against hospitalization with acute rotavirus gastroenteritis in Kenyan Children. Clin Infect Dis 2020; 70:2298–305.

    Article  CAS  PubMed  Google Scholar 

  17. Ong’ayo G, Ooko M, Wang’ondu R, et al. Effect of strikes by health workers on mortality between 2010 and 2016 in Kilifi, Kenya: a population-based cohort analysis. Lancet Glob Health 2019; 7:e961-e7.

    Google Scholar 

  18. Marongiu L, Shain E, Drumright L, Lillestøl R, Somasunderam D, Curran MD. Analysis of TaqMan Array Cards Data by an Assumption-Free Improvement of the maxRatio Algorithm Is More Accurate than the Cycle-Threshold Method. PLOS ONE 2016; 11:e0165282.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Talbert A, Ngari M, Bauni E, et al. Mortality after inpatient treatment for diarrhea in children: a cohort study. BMC Medicine. 2019;17:20.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Ruuska T, Vesikari T. Rotavirus disease in Finnish children: use of numerical scores for clinical severity of diarrhoeal episodes. Scand J Infect Dis 1990; 22:259–67.

    Article  CAS  PubMed  Google Scholar 

  21. (Kenya) NAaSCPN. The Kenya HIV testing services guidelines 3rd edition. ed: Ministry of Health 2015, 2015.

  22. Adetifa IMO, Bwanaali T, Wafula J, et al. Cohort Profile: The Kilifi Vaccine Monitoring Study. Int J Epidemiol 2017; 46:792-h.

    Google Scholar 

  23. Heaney J, Rolfe K, Gleadall NS, Greatorex JS, Curran MD. Low-Density TaqMan® Array Cards for the Detection of Pathogens. Methods in Microbiology 2015; 42:199–218.

    Article  PubMed Central  Google Scholar 

  24. Iturriza-Gómara M, Jere KC, Hungerford D, et al. Etiology of diarrhea among hospitalized children in Blantyre, Malawi, following rotavirus vaccine introduction: a case-control study. J Infectious Dis 2019; 220:213–8.

    Article  Google Scholar 

  25. Liu J, Platts-Mills JA, Juma J, et al. Use of quantitative molecular diagnostic methods to identify causes of diarrhoea in children: a reanalysis of the GEMS case-control study. Lancet 2016; 388:1291–301.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Platts-Mills JA, Amour C, Gratz J, et al. Impact of Rotavirus vaccine introduction and postintroduction etiology of diarrhea requiring hospital admission in in Haydom, Tanzania, a Rural African Setting. Clin Infect Dis 2017; 65:1144–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Operario DJ, Platts-Mills JA, Nadan S, et al. Etiology of severe acute watery diarrhea in children in the global rotavirus surveillance network using quantitative polymerase chain reaction. J Infect Dis 2017; 216:220–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Iturriza-Gómara M, Jere KC, Hungerford D, et al. Etiology of diarrhea among hospitalized children in Blantyre, Malawi, following rotavirus vaccine introduction: a case-control study. J Infect Dis 2019; 220:213–8.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Alsuwaidi AR, Al Dhaheri K, Al Hamad S, et al. Etiology of diarrhea by multiplex polymerase chain reaction among young children in the United Arab Emirates: a case-control study. BMC Infect Dis 2021; 21:7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Ballard SB, Requena D, Mayta H, et al. Enteropathogen changes after rotavirus vaccine scale-up. Pediatrics. 2022;149.

    Article  PubMed  Google Scholar 

  31. Harrison CJ, Hassan F, Lee B, et al. Multiplex PCR pathogen detection in acute gastroenteritis among hospitalized us children compared with healthy controls duringin the post-rotavirus vaccine era. Open Forum Infectious Diseases. 2021;8:ofab592.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Platts-Mills JA, Liu J, Rogawski ET, et al. Use of quantitative molecular diagnostic methods to assess the aetiology, burden, and clinical characteristics of diarrhoea in children in low-resource settings: a reanalysis of the MAL-ED cohort study. Lancet Glob Health 2018; 6:e1309-e18.

    Google Scholar 

  33. RonoSalinah J, Rose K, Fabian E. Prevalence and clinico-demographic characteristics associated with bacterial diarrhea among HIV positive and negative children aged below five years at Moi Teaching and Referral Hospital, Kenya. American Journal of Life Sciences 2014; 2:1–8.

    Google Scholar 

  34. Kotloff KL, Nasrin D, Blackwelder WC, et al. The incidence, aetiology, and adverse clinical consequences of less severe diarrhoeal episodes among infants and children residing in low-income and middle-income countries: a 12-month case-control study as a follow-on to the Global Enteric Multicenter Study (GEMS). Lancet Glob Health 2019; 7:e568-e84.

    Article  Google Scholar 

  35. Horwood PF, Soli KW, Maure T, et al. A high burden of asymptomatic gastrointestinal infections in traditional communities in Papua New Guinea. Am J Trop Med Hyg 2017; 97:1872–5.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Rogawski ET, Guerrant RL, Havt A, et al. Epidemiology of enteroaggregative Escherichia coli infections and associated outcomes in the MAL-ED birth cohort. PLoS Negl Trop Dis 2017; 11:e0005798.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Lugo D, Krogstad P. Enteroviruses in the early 21st century: new manifestations and challenges. Curr Opin Pediatr 2016; 28:107–13.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Otieno GP, Bottomley C, Khagayi S, et al. Impact of the introduction of rotavirus vaccine on hospital admissions for diarrhea among children in Kenya: a controlled interrupted time-series analysis. Clin Infect Dis 2020; 70:2306–13.

    Article  CAS  PubMed  Google Scholar 

  39. Zhou HL, Bessey T, Wang SM, et al. Burden and etiology of moderate and severe diarrhea in children less than 5 years of age living in north and south of China: prospective, population-based surveillance. Gut Pathog 2021; 13:33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Moutelíková R, Dvořáková Heroldová M, Holá V, Sauer P, Prodělalová J. Human rotavirus A detection: comparison of enzymatic immunoassay and rapid chromatographic test with two quantitative RT-PCR assays. Epidemiol Mikrobiol Imunol 2018; 67:110–3.

    PubMed  Google Scholar 

  41. Soares-Weiser K, Maclehose H, Bergman H, et al. Vaccines for preventing rotavirus diarrhoea: vaccines in use. Cochrane Database Syst Rev 2012; 11:Cd008521.

    PubMed  Google Scholar 

  42. Gurgel RQ, Correia JB, Cuevas LE. Effect of rotavirus vaccination on circulating virus strains. Lancet 2008; 371:301–2.

    Article  PubMed  Google Scholar 

  43. Esona MD, Gautam R, Katz E, et al. Comparative genomic analysis of genogroup 1 and genogroup 2 rotaviruses circulating in seven US cities, 2014–2016. Virus Evol 2021; 7:veab023.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Mwanga MJ, Owor BE, Ochieng JB, et al. Rotavirus group A genotype circulation patterns across Kenya before and after nationwide vaccine introduction, 2010–2018. BMC Infect Dis 2020; 20:504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Bar-Zeev N, Jere KC, Bennett A, et al. Population impact and effectiveness of monovalent rotavirus vaccination in urban malawian children 3 years after vaccine introduction: ecological and case-control analyses. Clin Infect Dis 2016; 62 Suppl 2:S213-9.

    PubMed  Google Scholar 

  46. Kang G, Iturriza-Gomara M, Wheeler JG, et al. Quantitation of group A rotavirus by real-time reverse-transcription-polymerase chain reaction: correlation with clinical severity in children in South India. J Med Virol 2004; 73:118–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Bennett A, Pollock L, Jere KC, et al. Duration and density of fecal rotavirus shedding in vaccinated Malawian children with rotavirus gastroenteritis. J Infect Dis 2020; 222:2035–40.

    Article  CAS  PubMed  Google Scholar 

  48. Kyo K, Takano C, Kasuga Y, et al. Severe rotavirus gastroenteritis in children older than 5 years after vaccine introduction. J Infect Chemother 2021; 27:598–603.

    Article  CAS  PubMed  Google Scholar 

  49. Shim JO, Chang JY, Shin S, Moon JS, Ko JS. Changing distribution of age, clinical severity, and genotypes of rotavirus gastroenteritis in hospitalized children after the introduction of vaccination: a single center study in Seoul between 2011 and 2014. BMC Infect Dis 2016; 16:287.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. Lucero Y, Matson DO, Ashkenazi S, George S, O’Ryan M. Norovirus: facts and reflections from past, present, and future. Viruses. 2021;13:2399.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Golicha Q, Shetty S, Nasiblov O, et al. Cholera Outbreak in Dadaab Refugee Camp, Kenya—November 2015–June 2016. MMWR Morb Mortal Wkly Rep 2018; 67:958–61.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Liu J, Gratz J, Amour C, et al. A laboratory-developed TaqMan Array Card for simultaneous detection of 19 enteropathogens. J Clin Microbiol 2013; 51:472–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Schrader C, Schielke A, Ellerbroek L, Johne R. PCR inhibitors—occurrence, properties and removal. J Appl Microbiol 2012; 113:1014–26.

    Article  CAS  PubMed  Google Scholar 

  54. Maldonado-Barrueco A, García-Rodríguez J, Yániz-Ramirez J, et al. Impact of the SARS-CoV-2 pandemic on the prevalence and incidence of gastrointestinal viruses in children up to five years old: a retrospective cohort study. Microbiol Spectr. 2022;10:e0266921.

    Article  PubMed  CAS  Google Scholar 

  55. Ahn SY, Park JY, Lim IS, et al. Changes in the occurrence of gastrointestinal infections after COVID-19 in Korea. J Korean Med Sci 2021; 36:e180.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Apondi Wandera E, Muriithi B, Kathiiko C, et al. Impact of integrated water, sanitation, hygiene, health and nutritional interventions on diarrhoea disease epidemiology and microbial quality of water in a resource-constrained setting in Kenya: a controlled intervention study. Trop Med Int Health. 2022.

    Article  Google Scholar 

  57. Lappan R, Jirapanjawat T, Williamson DA, Lange S, Chown SL, Greening C. Simultaneous detection of multiple pathogens with the TaqMan Array Card. MethodsX 2022; 9:101707.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We thank the study participants who provided the samples and members of the Virus Epidemiology and Control (VEC) group at KEMRI-Wellcome Trust Programme who collected the samples and performed the laboratory processing. We thank Prof. Philip Bejon for stimulating discussions about the work and comments provided on the initial draft manuscript. This study is published with permission of the from Director of KEMRI.


The views expressed in this article are solely those of the authors and not of the funders or affiliated institutions.


This study was funded by the Cambridge-Africa ALBORADA Research Fund. Dr Charles Agoti was supported by the Initiative to Develop African Research Leaders (IDeAL) through the DELTAS Africa Initiative [DEL-15-003]. The DELTAS Africa Initiative is an independent funding scheme of the African Academy of Sciences (AAS)’s Alliance for Accelerating Excellence in Science in Africa (AESA) and supported by the New Partnership for Africa’s Development Planning and Coordinating Agency (NEPAD Agency). Dr Lydia Drumright was supported by the following funding sources during this study: Medical Research Council (award # MR/S013164/1), the Alan Turing Institute through a Research Fellowship, the National Institutes for Health Research (award # CDF-2011-04-017) and the Cambridge Biomedical Research Center, the Economic and Social Science Research Council (award # AH/R001952/1), and an Isaac Newton/Wellcome Trust ISSF award to understand gastrointestinal infection epidemiology through use of this TAC. The views expressed in this report are those of the authors and not necessarily those of AAS, NEPAD Agency, The Wellcome, the UK government.

Author information

Authors and Affiliations



CNA, BB, SDWF, DJN, and LND sourced the study funding. CNA, EM, DJN MDC designed the study laboratory assay (the GASTRO IVK Taqman Array Card). CAN. BB, EM, BB and LND wrote the study research protocol. AWL and CNA did the laboratory experiments. NM and MN managed the study data and did the data analysis. CNA, NM and LND wrote first draft of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Charles N. Agoti.

Ethics declarations

Competing interests

All the authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Details of the targets on the Gastro v4.0K

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Agoti, C.N., Curran, M.D., Murunga, N. et al. Differences in epidemiology of enteropathogens in children pre- and post-rotavirus vaccine introduction in Kilifi, coastal Kenya. Gut Pathog 14, 32 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: