- Research note
- Open Access
Diarrheagenic Escherichia coli in Costa Rican children: a 9-year retrospective study
© The Author(s) 2019
- Received: 28 February 2019
- Accepted: 11 May 2019
- Published: 28 May 2019
This study aimed to estimate diarrheagenic Escherichia coli (DEC) prevalence among pediatric patients with diarrhea at the Costa Rican National Children’s Hospital-Social Security Service (Hospital Nacional de Niños-Caja Costarricense del Seguro Social; HNN-CCSS). DEC variations with respect to rainfall, presence of coinfections, and DEC antimicrobial resistance were also investigated.
A retrospective observational study from January 2008 to December 2016 was conducted. A total of 12 247 gastroenteritis records were analyzed. Annual DEC prevalence ranged from 2.7% (2008) to 9.0% (2013). The most prevalent pathotypes were enteroaggregative E. coli (EAEC) [n = 189 (31%)], enteropathogenic E. coli (EPEC) [n = 145 (24%)] and enteroinvasive E. coli (EIEC) [n = 91 (15%)]. A reduction in the probability of EAEC gastroenteritis was detected as rainfall rose above 200 mm/mo. [(Generalized Additive Model (GAM), p = 0.04)]. Coinfections were observed mainly between EPEC and Campylobacter spp. (10%). Antimicrobial resistance occurred in 0.6%, 29%, and 42% of DEC for ciprofloxacin, trimethoprim/sulfamethoxazole, and ampicillin, respectively.
- Diarrheagenic Escherichia coli
Annually, diarrheal diseases account for more than two million deaths worldwide. Of these, more than half a million correspond to children under 5 years of age . Diarrhea is defined as an increase in water content, volume, or frequency of stools, or as the passage of three or more loose or liquid stools per day [1, 2]. The broad spectrum of etiological agents, clinical syndromes and the multiple routes of pathogen transmission, highlight the need for accurate diagnostic strategies, clinical management and public health control programs . However, this approach is often challenging in developing countries.
Since 2005, the Costa Rican National Children´s Hospital –Social Security Service (Hospital Nacional de Niños-Caja Costarricense del Seguro Social; HNN-CCSS) hosts the only Costa Rican state-run laboratory that performs molecular testing for diarrheagenic Escherichia coli (DEC), identifying six well-known pathotypes: Enteropathogenic E. coli (EPEC), enterotoxigenic E. coli (ETEC), shiga-toxin producing E. coli (STEC), enteroinvasive E. coli (EIEC), enteroaggregative E. coli (EAEC), and diffusely adherent E. coli (DAEC). Moreover, infections provoked by mixed pathotypes of E. coli (MPEC) can be identified when genes belonging to different pathotypes are concurrently found in a given strain. In Costa Rica, a previous study showed a predominant circulation of EPEC and EIEC between 2005 and 2007 . The possible interaction between DEC and rainfall remained unexplored in Costa Rica. The aim of this study was to evaluate DEC prevalence and its association with rainfall, co-infections and resistance to antimicrobials in pediatric patients who were treated at the HNN-CCSS over a 9-year period (2008–2016).
Data included in the study
All stool records suggestive of gastroenteritis (loose or liquid samples, with or without leukocytes) between January 2008 and December 2016 were extracted from HNN-CCSS’s laboratory information systems: Labcore® (Rochem Biocare, Bogotá, Colombia) and Copernico® (Biomérieux, Marcy L’étoile, France) and collated with Microsoft Excel 2016. Diarrheic samples from outpatients under 13 years of age living in the central region of the country were included. Only the first sample per diarrheal episode was considered for analysis. Rainfall records for the study period were obtained from manual and automated weather stations located in the central region of Costa Rica and were kindly provided by the Costa Rican National Meteorological Institute.
All stool samples were analyzed at the time of sampling by the laboratory staff of the hospital. The standard laboratory procedure remained, in essence, unaltered throughout the study period. Predominant growth of yellow lactose-fermenting or bluish non-fermenting colonies on Tergitol-7 agar (Thermo Scientific™ Oxoid™) were investigated as a potential case of DEC gastroenteritis. Suspicious colonies were analyzed through an endpoint multiplex-PCR, as previously described . The assay targets genes eaeA, bfpA (EPEC); stx1, stx2 (STEC); ial (EIEC); lt, st (ETEC); aggR, astA (EAEC); and daaE (DAEC). Biochemical identification of isolates and antimicrobial susceptibility tests were performed using Vitek 2® (BioMérieux) and Clinical Laboratory Standards Institute (CLSI) breakpoints.
Descriptive statistics were used to estimate the frequency of pathogens, antimicrobial susceptibility and coinfections for the study period. Differences in the proportion of pathogens between years were ascertained through adjusted loglinear models . Generalized additive models (GAM) with Poisson distribution were fitted to evaluate the relationship between rainfall and the number of cases per pathotype. The impact of rainfall, DEC pathotypes, and their interaction with diarrheal cases were evaluated using a GAM model with a binomial distribution. The models were adjusted if they showed overdispersion and validated through diagnostic graphs. The level of significance was set at 0.05 (p < 0.05). The statistical analyses were carried out with R version 3.4.1 . The GAM models were adjusted with the MGCV  and MASS  packages, and the smoothed functions that were significant with the VISREG package  and GGPLOT2  were plotted.
A total of 46 906 fecal samples from 2008 to 2016 were studied. Of these, 26% (n = 12 247) matched the definition criteria for diarrhea [2008 (n = 1068); 2009, (n = 1054); 2010, (n = 1301); 2011, (n = 1046); 2012, (n = 1319); 2013, (n = 1249); 2014, (n = 1341); 2015, (n = 1964); 2016, (n = 1905)] and were selected for further analysis.
DEC prevalence throughout 2008–2016
DEC cases and prevalence were distributed as follows: [2008 (n = 29; 2.7%); 2009 (n = 31; 2.9%); 2010 (n = 64; 4.9%); 2011 (n = 56; 5.4%); 2012 (n = 98; 7.4%); 2013 (n = 112; 9.0%); 2014 (n = 65; 4.8%); 2015 (n = 86; 4.4%); 2016 (n = 68; 3.6%)]. Significant differences in prevalence were observed throughout the study period (Log-linear regression (LR), χ2 = 84.2; p < 0.0001). Notably, 2013 showed a nearly 3-fold increase in DEC prevalence (9.0%) with respect to 2008, whereupon prevalence decreased to 3.6% (n = 68) by 2016.
Distribution of DEC pathotypes
DEC pathotypes and rainfall records
Coinfections were found across most pathotypes, except for MPEC, and the global prevalence was estimated at 11% (n = 64). Frequency of co-infections ranged between pathotypes [EPEC 18% (n = 26), ETEC 13% (n = 9), EIEC 13% (n = 12), DAEC 11% (n = 6), EAEC 5% (n = 10), STEC 2% (n = 1)]. EPEC-Campylobacter coinfection was the most prevalent (10%; n = 14).
Antimicrobial susceptibility to ampicillin (AMP), trimethoprim/sulfamethoxazole (TMP/STX) and ciprofloxacin (CIP) differed among DEC pathotypes. Approximately 25-30% of EPEC, ETEC, STEC and EIEC isolates showed resistance to AMP, whereas this proportion increased to 50% for MPEC, 54% for EAEC, and 77% for DAEC. TMP/STX resistance was found in 16% of EPEC and ETEC, 23.3% of STEC and EIEC, 33.9% of EAEC, 40% of MPEC, and 60% of DAEC. A generalized susceptibility to CIP was found across all DEC pathotypes. Over a span of 9 years, only one STEC isolate showed resistance to all three antimicrobials.
To estimate the burden of DEC gastroenteritis, we analyzed more than 46 906 laboratory records, of which only 26% (n = 12 247) fulfilled the inclusion criteria. DEC was identified in 609 (5.0%) of these samples. A similar study found DEC in 7.44% of cases in children under 5 years of age , whilst a small study found it in 6/155 (3.9%) of single infections using a molecular gastrointestinal panel . We previously demonstrated the concurrent circulation of all DEC pathotypes and a predominance of EPEC and EIEC in a study from 2005 to 2007 . The changing epidemiology of this pathogen has been described in different studies and geographical regions [12–18]. Importantly, this work provides evidence showing this phenomenon over an extended period of evaluation. This study unveiled a gradual shift in the contribution of specific pathotypes to DEC gastroenteritis, wherein the share of EAEC cases progressively increased as those attributable to EPEC declined. EAEC has been recently acknowledged as the leading DEC pathotype [13, 15, 19–23].
Regarding STEC, its prevalence remained low over the span of this study. Nonetheless, its association with large outbreaks, hemolytic uremic syndrome (HUS) and the high mortality rates emphasize the relevance of this pathotype for public health [24, 25]. To the best of our knowledge, asymptomatic STEC carriage in children  has not been studied thus far in Costa Rica, thereby hampering the formulation of policies to prevent its dissemination.
DEC and rainfall
The ubiquity of E. coli is well-documented and makes this bacterium a classic indicator of fecal contamination in water sources . DEC presence in water used for human activities emphasizes the need to track not only the frequency of diarrheal disease but also the spatio-temporal interaction between the environment and this pathogen [28–30]. By virtue of the high variability of rainfall in tropical countries, we hereby described an innovative approach, analyzing DEC gastroenteritis in terms of monthly rainfall instead of seasonality. Using this approach, we found a marked drop in the EAEC infection probability and rainfall above 200 mm. We hypothesize that sudden rain showers and incipient rains heralding the rainy season might trigger EAEC diarrheal cases. In contrast, all other DEC pathotypes did not seem to interact with rainfall. Additional environmental and human factors might better predict the DEC infection probability by pathotypes other than EAEC.
We found 11% of coinfections between DEC and other enteropathogens. Changes in gut microbiota induced by DEC infection and the potential interaction with other pathogens have been reviewed elsewhere . Campylobacter spp. was the most prevalent pathogen (4%) co-infecting pediatric patients in this study, which were predominantly associated with EPEC (10% of the cases) and ETEC (6% of the cases). The coinfection between EPEC and Campylobacter was also found at similar rates in a multicenter study in Europe, using a multiplex PCR-based approach .
We found different antimicrobial susceptibility patterns among DEC pathotypes. Notably, 54% of EAEC isolates were mono-resistant to AMP, 34% to TMP/STX, and 33% showed resistance to both. Resistance to multiple antimicrobials is a major concern in EAEC , considering its high prevalence in children under 5 years of age and the possibility of asymptomatic carriage [33, 34]. Unlike the previous study  by Pérez et al. we did not find marked susceptibility to AMP in EPEC. Rather, a 25.6% resistance rate was observed, which was similar to ETEC (24%) and suggested increased rates of resistance in these two pathotypes. We hypothesize that this shift may be explained by a higher circulation, in recent years, of DEC strains carrying mobile genetic elements, which confer resistance to AMP. The rapid changes in local antibiotic susceptibility justify the active surveillance of DEC for epidemiological and public health purposes.
For more than a decade, DEC has been investigated in stool specimens via molecular methods at the HNN-CCSS. Through this work, we are now aware of the local distribution, prevalence, antimicrobial resistance, and the interactions between each DEC pathotype and rainfall for most of this period.
Since March 2017, a syndromic multiplex PCR is part of the routine tests performed directly on diarrheic stool specimens for the simultaneous detection of multiple pathogens. This novel molecular approach will improve diagnostic efficiency through reduced turnaround times and an extended coverage of agents targeted in a single assay, thereby enabling a timely management and the instauration of appropriate antimicrobial therapy when warranted. The overall impact of this strategy will require an exhaustive review, for which this study establishes a reference against which future research can draw upon to determine the influence of technical innovation on the epidemiological profile of infectious gastroenteritis at the HNN-CCSS.
It is possible that a minor fraction of data may be missing due to the retrospective nature of this study. Additionally, this is a hospital and laboratory-based study only. Consequently, it does not include all children diagnosed with DEC gastroenteritis in the central region of Costa Rica. It is likely that the actual DEC pediatric disease was greater than depicted in this study.
Special thanks to Drs. Marco Herrera, Álvaro Vargas, Carlos Santamaría and Mario Vargas for the support during the data collection and elaboration of this manuscript.
The authors declare that they have not received any funding for the elaboration of the present research.
CPC participated in the conception, elaboration, execution and analysis of data. Also was a major contributor in the writing of the manuscript. KLS analyzed, depurated databases, and was a major contributor in the writing of the manuscript. Both authors read and approved the final manuscript.
Ethics approval and consent to participate
Ethical approval was obtained from the HNN-CCSS Scientific Ethics Committee of the (Ref. CLOBI-HNN-015-2015).
Consent for publication
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- WHO. Diarrhoeal disease. WHO. 2017. http://www.who.int/news-room/fact-sheets/detail/diarrhoeal-disease. Accessed 19 Aug 2018.
- Guerrant RL, Van Gilder T, Steiner TS, Thielman NM, Slutsker L, Tauxe RV, et al. Practice guidelines for the management of infectious diarrhea. Clin Infect Dis. 2001;32(3):331–51.View ArticleGoogle Scholar
- Pérez C, Gómez-Duarte OG, Arias ML. Diarrheagenic Escherichia coli in Children from Costa Rica. Am J Trop Med Hyg. 2010;83(2):292–7.View ArticleGoogle Scholar
- Bilder CR, Loughin TM. Analyzing a binary response, part 2: regression models. In: Analysis of categorical data with R. Boca Raton, FL: Chapman and Hall/CRC; 2014. p. 121–9.Google Scholar
- R Development Core Team RFFSC. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. 2008.Google Scholar
- Wood SN. GAMs in Practice: mgcv. In: Generalized additive models: an introduction with R. 2nd ed. New York: Chapman and Hall/CRC; 2017. p. 80.Google Scholar
- Venables WN, Ripley BD. Modern applied statistics with S. 4th ed. New York: Springer; 2002.View ArticleGoogle Scholar
- Breheny P, Burchett W. Visualizing regression models using visreg. R J. 2017;9(2):56–71.View ArticleGoogle Scholar
- Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer; 2009.View ArticleGoogle Scholar
- Gómez-Duarte OG, Romero-Herazo YC, Paez-Canro CZ, Eslava-Schmalbach JH, Arzuza O. Enterotoxigenic Escherichia coli associated with childhood diarrhoea in Colombia, South America. J Infect Dev Ctries. 2013;7:372–81.View ArticleGoogle Scholar
- López-Medina E, Parra B, Dávalos DM, López P, Villamarín E, Pelaez M. Acute gastroenteritis in a pediatric population from Cali, Colombia in the post rotavirus vaccine era. Int J Infect Dis. 2018;73:52–9.View ArticleGoogle Scholar
- Torres AG. Escherichia coli diseases in Latin America-a “One Health” multidisciplinary approach. Pathog Dis. 2017;75(2):ftx012.View ArticleGoogle Scholar
- Imdad A, Foster MA, Iqbal J, Fonnesbeck C, Payne DC, Zhang C, et al. Diarrheagenic Escherichia coli and Acute Gastroenteritis in Children in Davidson County, Tennessee, United States: A Case-control Study. Pediatr Infect Dis J. 2018;37(6). https://journals.lww.com/pidj/Fulltext/2018/06000/Diarrheagenic_Escherichia_coli_and_Acute.10.aspx.View ArticleGoogle Scholar
- Weiler N, Orrego M, Alvarez M, Huber C. Detección molecular de Escherichia coli diarreogénica en pacientes pediátricos con síndrome diarreico agudo en Paraguay. Mem Inst Investig Cienc Salud. 2017;15(1):16–21.View ArticleGoogle Scholar
- Canizalez-Roman A, Flores-Villaseñor HM, Gonzalez-Nuñez E, Velazquez-Roman J, Vidal JE, Muro-Amador S, et al. Surveillance of diarrheagenic Escherichia coli strains isolated from diarrhea cases from children, adults and elderly at Northwest of Mexico. Front Microbiol. 2016;7:1924.View ArticleGoogle Scholar
- Rivas M, Padola NL, Lucchesi PM, Masana M. Diarrheagenic Escherichia coli in Argentina. In: Torres AG, editor. Pathogenic Escherichia coli in Latin America. UAE: Bentham Books; 2010. p. 142–61.Google Scholar
- Vilchez S, Reyes D, Paniagua M, Bucardo F, Möllby R, Weintraub A. Prevalence of diarrhoeagenic Escherichia coli in children from León, Nicaragua. J Med Microbiol. 2009;58(5):630–7. https://doi.org/10.1099/jmm.0.007369-0.View ArticlePubMedGoogle Scholar
- Spano LC, da Cunha KF, Monfardini MV, de Cássia Bergamaschi Fonseca R, Scaletsky ICA. High prevalence of diarrheagenic Escherichia coli carrying toxin-encoding genes isolated from children and adults in southeastern Brazil. BMC Infect Dis. 2017;17(1):773. https://doi.org/10.1186/s12879-017-2872-0.View ArticlePubMedPubMed CentralGoogle Scholar
- Lozer DM, Souza TB, Monfardini MV, Vicentini F, Kitagawa SS, Scaletsky ICA, et al. Genotypic and phenotypic analysis of diarrheagenic Escherichia coli strains isolated from Brazilian children living in low socioeconomic level communities. BMC Infect Dis. 2013;13(1):418. https://doi.org/10.1186/1471-2334-13-418.View ArticlePubMedPubMed CentralGoogle Scholar
- Dutta S, Guin S, Ghosh S, Pazhani GP, Rajendran K, Bhattacharya MK, et al. Trends in the prevalence of diarrheagenic Escherichia coli among hospitalized diarrheal patients in Kolkata, India. PLoS One. 2013;8(2):1–5. https://doi.org/10.1371/journal.pone.0056068.View ArticleGoogle Scholar
- Ali MMM, Mohamed ZK, Klena JD, Ahmed SF, Moussa TAA, Ghenghesh KS. Molecular characterization of diarrheagenic Escherichia coli from Libya. Am J Trop Med Hyg. 2012;86:866–71.View ArticleGoogle Scholar
- Garcia PG, Silva VL, Diniz CG. Occurrence and antimicrobial drug susceptibility patterns of commensal and diarrheagenic Escherichia coli in fecal microbiota from children with and without acute diarrhea. J Microbiol. 2011;49(1):46–52. https://doi.org/10.1007/s12275-011-0172-8.View ArticlePubMedGoogle Scholar
- Somda NS, Bonkoungou OJI, Zongo C, Kpoda DS, Tapsoba F, Bassolé IHN, et al. Prevalence of Escherichia coli virulence genes in patients with diarrhoea in Ouagadougou, Burkina Faso. Afr J Clin Exp Microbiol. 2017;18(4):179–85.View ArticleGoogle Scholar
- Paton JC, Paton AW. Pathogenesis and Diagnosis of Shiga Toxin-Producing Escherichia coli Infections. Clin Microbiol Rev. 1998;11:450–79.View ArticleGoogle Scholar
- Nataro JP, Kaper JB. Diarrheagenic Escherichia coli. Clin Microbiol Rev. 1998;11(1):142–201. http://cmr.asm.org/content/11/1/142.abstract.View ArticleGoogle Scholar
- Harries M, Dreesman J, Rettenbacher-Riefler S, Mertens E. Faecal carriage of extended-spectrum β-lactamase-producing Enterobacteriaceae and Shiga toxin-producing Escherichia coli in asymptomatic nursery children in Lower Saxony (Germany), 2014. Epidemiol Infect. 2016;144(16):3540–8.View ArticleGoogle Scholar
- Edberg SC, Rice EW, Karlin RJ, Allen MJ. Escherichia coli: the best biological drinking water indicator for public health protection. J Appl Microbiol. 2000;88:106S–16S.View ArticleGoogle Scholar
- Ramírez Castillo FY, Avelar González FJ, Garneau P, Márquez Díaz F, Guerrero Barrera AL, Harel J. Presence of multi-drug resistant pathogenic Escherichia coli in the San Pedro River located in the State of Aguascalientes, Mexico. Front Microbiol. 2013;4:147.View ArticleGoogle Scholar
- Estrada-Garcia T, Lopez-Saucedo C, Thompson-Bonilla R, Abonce M, Lopez-Hernandez D, Santos JI, et al. Association of Diarrheagenic Escherichia coli Pathotypes with Infection and Diarrhea among Mexican Children and Association of Atypical Enteropathogenic E. coli with Acute Diarrhea. J Clin Microbiol. 2009;47(1):93–8. https://jcm.asm.org/content/47/1/93.View ArticleGoogle Scholar
- Paredes-Paredes M, Okhuysen PC, Flores J, Mohamed JA, Padda RS, Gonzalez-Estrada A, et al. Seasonality of diarrheagenic Escherichia coli pathotypes in the US students acquiring diarrhea in Mexico. J Travel Med. 2011;18(2):121–5. https://doi.org/10.1111/j.1708-8305.2010.00488.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Gallardo P, Izquierdo M, Vidal RM, Chamorro-Veloso N, Rosselló-Móra R, O’Ryan M, et al. Distinctive gut microbiota is associated with diarrheagenic Escherichia coli infections in chilean children. Front Cell Infect Microbiol. 2017;7:424.View ArticleGoogle Scholar
- Spina A, Kerr KG, Cormican M, Barbut F, Eigentler A, Zerva L, et al. Spectrum of enteropathogens detected by the FilmArray GI panel in a multicentre study of community-acquired gastroenteritis. Clin Microbiol Infect. 2015;21(8):719–28.View ArticleGoogle Scholar
- Amaya E, Reyes D, Vilchez S, Paniagua M, Möllby R, Nord CE, et al. Antibiotic resistance patterns of intestinal Escherichia coli isolates from Nicaraguan children. J Med Microbiol. 2011;60(2):216–22. http://jmm.microbiologyresearch.org/content/journal/jmm/10.1099/jmm.0.020842-0.View ArticleGoogle Scholar
- Hebbelstrup Jensen B, Olsen KEP, Struve C, Krogfelt KA, Petersen AM. Epidemiology and clinical manifestations of enteroaggregative Escherichia coli. Clin Microbiol Rev. 2014;27:614–30.View ArticleGoogle Scholar