Introduction 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 [ 1 ]. 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 [ 2 ]. 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 [ 3 ]. 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). Main text Methods 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. Laboratory workup 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 [ 3 ]. 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. Statistical analyses 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 [ 4 ]. 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 &lt; 0.05). The statistical analyses were carried out with R version 3.4.1 [ 5 ]. The GAM models were adjusted with the MGCV [ 6 ] and MASS [ 7 ] packages, and the smoothed functions that were significant with the VISREG package [ 8 ] and GGPLOT2 [ 9 ] were plotted. Results Database records 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 &lt; 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 The proportion of DEC pathotypes differed significantly over the years ( LR, χ 2 = 195.6; p &lt; 0.0001 ). EAEC was the most frequent pathotype identified [31% (n = 189)], followed by EPEC [24% (n = 145)], EIEC [15% (n = 91)], ETEC [11% (n = 70)], DAEC [9% (n = 54)], STEC [8% (n = 46)], and MPEC [2% (n = 14)] (see Fig. 1 ). Fig. 1 Distribution of DEC pathotypes from 2008 to 2016 DEC pathotypes and rainfall records Rainfall records depicted two well-defined seasons (dry and rainy), as well as noticeable variations in monthly rainfall averages between the years (Fig. 2 a). The interaction between monthly rainfall and DEC cases was significant (GAM; χ 2 = 33.22; p &lt; 0.05). This interaction had a sigmoidal pattern, as shown in Fig. 2 b. Fig. 2 a Radial plot showing monthly average rainfall (mm) throughout the study period (2008–2016). b Relationship between cases of DEC diarrheal disease and monthly rainfall (2008–2016). The line represents the predicted smoothed function for DEC and the grey-shaded area corresponds to 95% point-wise confidence bands We also estimated the infection probability for each pathotype. Among the pathotypes investigated, only the EAEC infection probability and monthly rainfall was significant (Fig. 3 ). Fig. 3 DEC infection probability changes and monthly rainfall (mm). The lines represent the predicted smoothed functions of DEC pathotypes and the grey-shaded area the 95% point-wise confidence bands Coinfections 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 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. Discussion 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 [ 10 ], whilst a small study found it in 6/155 (3.9%) of single infections using a molecular gastrointestinal panel [ 11 ]. We previously demonstrated the concurrent circulation of all DEC pathotypes and a predominance of EPEC and EIEC in a study from 2005 to 2007 [ 3 ]. 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 [ 26 ] 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 [ 27 ]. 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. Coinfections 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 [ 31 ]. 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 [ 32 ]. Antimicrobial susceptibility 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 [ 15 ], considering its high prevalence in children under 5 years of age and the possibility of asymptomatic carriage [ 33 , 34 ]. Unlike the previous study [ 3 ] 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. Conclusions 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. Limitations 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.
Original work titled: 'Diarrheagenic Escherichia coli in Costa Rican children: a 9-year retrospective study' (BMC Research Notes/BioMed Central) by the following authors: Cristian Pérez-Corrales; Kevin Leandro-Sandí DOI: 10.1186/s13104-019-4313-1