Downloaded fr by kenia.barra on 01 March 2 Occurrence of pathogenic microorganisms in small drinking-water systems in Costa Rica Kenia Barrantes a,*, Luz Chacón a, Eric Morales a, Luis Rivera-Montero a, Macario Pino b, Alejandra Gamboa Jiménezc, Diana Campos Morac, Pablo Salas Jiménez c, Basilio Silvad and Luis G. Romero-Esquivel b a Infection and Nutrition Section, Health Research Institute, University of Costa Rica, P.O. Box 11501-2060, San José, Costa Rica b Environmental Protection Research Center (CIPA), School of Chemistry, Instituto Tecnológico de Costa Rica (ITCR), P.O. Box 159-7050, Cartago, Costa Rica c Environmental Analysis Laboratory, School of Environmental Sciences, National University, P.O. Box 86-3000, Heredia, Costa Rica d Basic Sciences, National Technological University, P.O. Box 1902–4050, Alajuela, Costa Rica *Corresponding author. E-mail: kenia.barrantes@ucr.ac.cr KB, 0000-0002-2673-9220; LC, 0000-0003-2506-0619; EM, 0000-0001-8090-4947; LR-M, 0000-0001-8739-7466; MP, 0000-0001-8446-4723; PSJ, 0000-0002-7759-6741; LGR-E, 0000-0002-7191-3195 © 2022 The Authors Journal of Water and Health Vol 20 No 2, 344 doi: 10.2166/wh.2022.230 ABSTRACT This study describes the quality of drinking water sampled over 2 years (2018 and 2019) from 20 ASADAS (Spanish acronym for Administrative Associations for Water and Sewer Systems) in Costa Rica. The analysis included Rotavirus (RV), somatic coliphages, fecal coliforms, and Escherichia coli. The ASADAS were categorized into three regions as temperate rainy (region 1), tropical rainy (region 2), and tropical rainy and dry (region 3) according to biogeographic classification. The concentrations of fecal coliforms and E. coli were higher in samples from surface water sources from the ASADAS in region 3 compared to regions 1 and 2. RV-positive samples (24/296) were detected in drink- ing-water samples from regions 2 and 3 during dry and transition seasons, with higher concentrations more frequently in the dry season. In addition, somatic coliphages were detected in samples from the three regions, with higher concentrations in region 2. Furthermore, a stat- istically significant relationship was found between somatic coliphages and diarrheal cases, classified as outbreaks or alerts in the region. Thus, the results confirmed that somatic coliphages are a good indicator of the presence of diarrhea cases in a specific region. Key words: Costa Rica, Escherichia coli, fecal coliforms, Rotavirus, somatic coliphages HIGHLIGHTS • Rotaviruses (RVs) are frequent waterborne pathogens. • Fecal coliforms, Escherichia coli, somatic coliphages, and RV were analyzed from drinking-water samples in Costa Rica. • A significant relationship was found between somatic coliphages and diarrheal cases. • Somatic coliphages are considered a good indicator of diarrhea-causing pathogens in a specific region. INTRODUCTION Drinking water is necessary for life and an essential human right. But, unfortunately, the World Health Organization (WHO) estimates that 700 million people worldwide are not receiving their drinking water from improved water sources, with sub- sequent devastating effects on their quality of life (UN-Water 2014; WHO 2017). A variety of pathogen microorganisms are transmitted by contaminated water. Bacteria, protozoa, helminths, and viruses are good examples of pathogens as etiological agents of waterborne diseases. The most significant microbial risks are water contamination with feces from humans or animals (including birds). Thus, feces can be a source of pathogenic microorgan- isms (WHO 2017). Fecal coliforms (FCs) and Escherichia coli are common bacterial indicators of fecal water contamination. They have tra- ditionally been used to monitor drinking-water quality, and they remain important parameters in monitoring undertaken as part of verification or surveillance. However, water intended for human consumption should contain no fecal indicator organisms (WHO 2017). This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/). om http://iwaponline.com/jwh/article-pdf/20/2/344/1013002/jwh0200344.pdf?guestAccessKey=96e7e231-769a-4283-96cd-579a91d25f30 ntes@ucr.ac.cr 022 https://orcid.org/0000-0002-2673-9220 https://orcid.org/0000-0003-2506-0619 https://orcid.org/0000-0001-8090-4947 https://orcid.org/0000-0001-8739-7466 https://orcid.org/0000-0001-8446-4723 https://orcid.org/0000-0002-7759-6741 https://orcid.org/0000-0002-7191-3195 mailto:kenia.barrantes@ucr.ac.cr http://orcid.org/ http://orcid.org/0000-0002-2673-9220 http://orcid.org/0000-0003-2506-0619 http://orcid.org/0000-0001-8090-4947 http://orcid.org/0000-0001-8739-7466 http://orcid.org/0000-0001-8446-4723 http://orcid.org/0000-0002-7759-6741 http://orcid.org/0000-0002-7191-3195 http://creativecommons.org/licenses/by/4.0/ http://creativecommons.org/licenses/by/4.0/ https://crossmark.crossref.org/dialog/?doi=10.2166/wh.2022.230&domain=pdf&date_stamp=2022-02-03 Journal of Water and Health Vol 20 No 2, 345 Downloaded from http by kenia.barrantes@u on 01 March 2022 Viruses are the smallest pathogens and may be less sensitive to water disinfection than bacteria. They usually persist for long periods in water, and infective doses are typically low. Rotaviruses (RVs), enteroviruses, and noroviruses have been identified as frequent waterborne pathogens among viruses (WHO 2017). They contaminate water through human activities such as leaking sewage and septic systems, agricultural and urban runoff, estuarine and marine waters, sewage outfall, and vessel wastewater discharge (Fong & Lipp 2005; WHO 2017). RVs are the most important cause of gastrointestinal infection in children, and this infection can have severe consequences, including hospitalization and death, particularly in low-income regions. Occasional outbreaks of waterborne disease have been recorded (WHO 2017). Before introducing RV vaccines in 2006, approximately 15,000 deaths, 75,000 hospitalizations, and two million medical visits were associated with RV infection in Latin America and the Caribbean annually (Ureña-Castro et al. 2019). The effec- tiveness of RV vaccines against hospitalizations and severe diarrhea in Latin American children has been approximately 70%, with notable reductions in mortality rates in Brazil, Mexico, and Panama. In Costa Rica, since 2019, two live oral RV vaccines have been included in the National Immunization Program. One of these vaccines is derived from an attenuated human strain of RV, and the other combines five bovine-human reassortant strains (Ulloa-Gutierrez & Avila-Aguero 2014; Ureña-Castro et al. 2019). On the other hand, coliphages are a type of bacteriophage, viruses that use only bacteria as hosts for replication. Coliphages use E. coli and closely related species as hosts and can be released into humans and other warm-blooded animals (WHO 2017). Moreover, they share many properties with human viruses, such as composition, morphology, structure, and replica- tion mode (Jofre et al. 2016). As a result, coliphages are applicable models of enteric viruses in water environments and the sensitivity to treatment and disinfection processes. Therefore, the WHO recommends coliphages for inclusion in verification and surveillance monitoring where source waters are affected by human fecal waste (WHO 2017). Drinking-water coverage in Costa Rica Costa Rica has been considered one of the countries in the Central American region with the best access to improved water sources, reducing the inequity of this resource among rural and urban populations (MINAET 2013). According to the National Institute of Statistics and Census (INEC), drinking-water coverage increased in almost 10 years, from 91.1% in 2010 to 94.1% in 2019 (INEC 2021). Regarding these data, coverage in the urban area increased from 96.7% in 2010 to 97.8% in 2019; however, for rural areas, it increased from 76.1 to 84% in the same period (INEC 2021). Since 2000, a legal framework for the co-management of rural aqueducts was created at the national level to address the pattern of inequity among rural and urban populations and achieve complete nationwide coverage (Dobbin 2013). This national policy created a new type of water management organization with the legal status of community water organizations and increased governmental oversight, called Administrative Associations for Water and Sewer Systems (Asociaciones Administradoras de Sistemas de Acueductos y Alcantarillados Sanitarios or ASADAS) (Dobbin 2013). More than 1,400 ASADAS in the country supply 25.5% of the Costa Rican population. The quality of the water and the supported inputs that an ASADA can obtain, such as technology, human, and financial resources, directly affect the drinking-water quality in a specific population (Soto-Córdoba et al. 2016). The degree of pollution of the source used for water supply is also an important parameter that could affect the drinking- water system’s quality. Surface waters in Costa Rica suffer potential contamination, mostly from untreated sewage, runoffs, and agricultural chemicals. This situation limits the value of surface sources for water supply (Bower 2014). Morales et al. (2019) described that the contamination of catchment areas used as resources in one drinking-water treatment plant in Costa Rica influences the water quality in the final distribution network, even if the water was treated. For that reason, inten- sive and sectorized source monitoring is recommended to locate specific contamination points (Morales et al. 2019). Another critical fact to consider is that access to drinking water does not necessarily imply that the resource is of potable quality. Potable means that it complies with the physical–chemical, biological, and microbiological characteristics according to sanitary regulations, and thus guarantees that its consumption does not imply health risks (WHO 2017). However, accord- ing to the National Water Laboratory of the Costa Rican Institute of Aqueducts and Sewerage (Instituto Costarricense de Acueductos y Alcantarillados, ICAA), in 2020, only 89% of the population supplied received potable drinking-water quality (LNA-AyA 2021), which means that 11% of Costa Rica’s population consume water of non-potable quality, thereby increasing the risk of developing waterborne diseases. Even though Costa Rica has good water sources and generally good water quality for consumption, there are many diarrhea cases annually (Chacón Jiménez et al. 2015). The National Water Laboratory of the ICAA analyzed 115 diarrheal outbreaks ://iwaponline.com/jwh/article-pdf/20/2/344/1013002/jwh0200344.pdf?guestAccessKey=96e7e231-769a-4283-96cd-579a91d25f30 cr.ac.cr Journal of Water and Health Vol 20 No 2, 346 Downloaded fr by kenia.barra on 01 March 2 related to waterborne pathogens from 1999 to 2005 (Valiente & Mora 2007). The study showed that in 70% of the outbreaks, bacterial pathogens or FCs were detected (Valiente & Mora 2007). Thirty per cent of outbreak cases’ water samples were tested negative for bacteriological indicators of fecal pollution, such as FCs (Valiente & Mora 2007). In another study, microbiological indicators of water contamination were analyzed in two rural community aqueducts in the Central Valley of Costa Rica for 1 year (Chacón Jiménez et al. 2015). The authors reported an association between viral indicators in thewater, such as somatic coliphages, and the increase in diarrhea cases in these communities. In those aqueducts, the water was treated using unsuitable technologies such as chlorination to eliminate viruses and parasites. The bacteriological indicators (FC and E. coli) were also analyzed, and they were absent in the samples (Chacón Jiménez et al. 2015). Although several studies described the detection of viral pathogens such as RV or Norovirus in surface and treated water (De Deus et al. 2019; Prez et al. 2020), information regarding the frequency of viral pathogens in water is scarce in the Latin American region. The present study describes the quality of drinking water sampled over 2 years (2018 and 2019) from 20 ASADAS in Costa Rica. The analysis included RVs, somatic coliphages, FCs, and E. coli. In addition, the ASADAS were categorized into three regions according to climate and humidity parameters. As a result, there were observed differences regarding microbiological parameters in the ASADAS from region 3 compared to regions 1 and 2. Also, endemic channels were elaborated to compare and discuss the diarrhea cases reported in those locations. METHODS Drinking-water sampling A total of 296 water samples (treated and non-treated) were collected from 20 ASADAS, located in Cartago, Puntarenas, and Alajuela provinces in Costa Rica (Figure 1). Figure 1 | Location of the 20 ASADAS in Costa Rica according to each climatic region. The scale denotes the distance in km. om http://iwaponline.com/jwh/article-pdf/20/2/344/1013002/jwh0200344.pdf?guestAccessKey=96e7e231-769a-4283-96cd-579a91d25f30 ntes@ucr.ac.cr 022 Journal of Water and Health Vol 20 No 2, 347 Downloaded from http by kenia.barrantes@u on 01 March 2022 Figure 1 shows the location of the three regions on the Costa Rica map. Region 1 was categorized as temperate rainy; region 2 was classified as tropical rainy, and finally, region 3 was tropical rainy and dry. Sample collection was performed during the transition and rainy seasons in 2018 (35 samples in April and May, 65 in Octo- ber and November, respectively) and dry, transition, and rainy seasons in 2019 (99 samples in March and April, 40 in May and June, and 57 in August and September, respectively). Each ASADA was visited three times, one per season. During the sampling period five water samples were collected per visit, one from the raw water source and four from the distribution network. Table 1 indicates the types of water source and treatment used for each ASADA. Climatic conditions The 20 ASADAS were categorized according to the Holdridge life zone system (Holdridge 1971) into three regions as tem- perate rainy (region 1), tropical rainy (region 2), and tropical rainy and dry (region 3) according to biogeographic classification. The climatic regions were obtained from the Atlas of Costa Rica (Ortiz-Malavassi 2014), which considers the accumulated annual precipitation range, mean temperature, and vegetation type. Thus, region 1 included nine ASADAS; region 2 included four ASADAS; and region 3, seven ASADAS. FC and E. coli analysis Microorganisms were enumerated using the most probable number technique (9221E; American Public Health Association 2017). Two hundred milliliters of water were collected using a sterile recipient. The sample was inoculated within less than 30 h of the collection into lauryl tryptose broth (Oxoid) and incubated at 35.0 °C for 48 h. All tubes testing positive after the incubation period were inoculated into EC-MUG broth (Oxoid). After a 24-h incubation period at 44.5 °C, tubes with a posi- tive reaction were checked and reported. E. coli ATCC 25922 was used as a positive control and Pseudomonas aeruginosa ATCC 14502 as a negative control. Somatic coliphages Somatic coliphages were enumerated using the single layer agar method 9924B and E Parts (American Public Health Association 2017) with modifications (Solano Barquero et al. 2012). Briefly, 250 mL of water was collected using a sterile recipient. The sample was first filtered using an 80-μm glass fiber filter (Sartorius Stedim Biotech, Goettingen, Germany) pretreated with beef extract pH 7.2 (Oxoid) to prevent losses associated with viruses sticking to the filter. Then, a second sample filtration using 0.2 μm cellulose acetate filters (Sartorius Stedim Table 1 | Types of water source, treatment, and disinfection systems used by each ASADA Raw water source type Water treatment system Water disinfection system Number of ASADAS Surface Conventional DWTPa Chlorine 1 Slow filtration Chlorine 1 Settler Chlorine 2 Desander Electrolysis 1 Multi-stage filtrationb Chlorine 1 Nonec None 1 Surface/groundwater Conventional DWTPa Chlorine 1 Coarse layered filtration Chlorine 2 Settler Chlorine 1 Desander Chlorine 1 Nonec Chlorine 2 Spring Nonec Chlorine 5 Nonec Electrolysis 1 Total of ASADAS 20 aDWTP represents the Conventional Drinking Water Treatment Plant, which uses a water treatment system, including coagulation, flocculation, sedimentation, and filtration. bMulti-stage filtration uses a water treatment system, including two pre-treatment units, such as gravel filters, with a final unit corresponding to slow sand filters. cNone means no treatment is used to treat water before disinfection. ://iwaponline.com/jwh/article-pdf/20/2/344/1013002/jwh0200344.pdf?guestAccessKey=96e7e231-769a-4283-96cd-579a91d25f30 cr.ac.cr Journal of Water and Health Vol 20 No 2, 348 Downloaded fr by kenia.barra on 01 March 2 Biotech, Goettingen, Germany) was pretreated with beef extract pH 7.2. Then it was mixed with 10 mL of a fresh culture of E. coli ATCC 13706 (O.D.¼0.300 nm), CaCl2, to a final concentration of 0.2 M with 2% trypticase soy agar (TSA). The mix was spread into Petri dishes and incubated at 35 °C overnight. After the incubation period, samples were analyzed to count for plaque-forming units (PFU). Rotavirus Five liters of water were collected using a sterile recipient. Then, according to reference, samples were processed and ana- lyzed to obtain viral quantification (number of copies/100 ml) (Chacón et al. 2020). Spatial description and data analysis The locations of the ASADAS were obtained in the field using a GPS Garmin, model Etrex touch 25. The GPX data were transformed into a vector format. The vectorial database was classified and handled using the SQL language. The climatic regions were obtained from the Atlas of Costa Rica as vector format (Ortiz-Malavassi 2014). The spatial description was performed using ArcGIS 10.7.1. (ESRI) with four data frames (one for each of the three regions and the fourth for the general location). Descriptive statistics, boxplot graphics, odds ratio, and logistic regression were performed using the R software from CRAN. In addition, censored data analysis and χ2 were performed to determine differences between regions and microbiological par- ameters. The NADA package of the R language for non-parametric censored data was used (Helsel 2012). Biogeographic regions were established as the analysis group, delimiting regions 1, 2, and 3. The selection of method for estimating descriptive statistics was determined according to the characteristics of the dataset (sample size and proportion of censored data) (Helsel 2012). For FC and E. coli, the Kaplan–Meier and Turnbull method was used. For somatic coliphages and RVs, the maximum likelihood estimation method was used. Region 2 was excluded from the analysis for somatic coli- phages since more than 80% of samples showed censored data. Differences between groups were obtained using the non- parametric Kendal–Tau correlation method for each microbiological variable. According to references, the registered cases of acute diarrheal disease for regions 1, 2, and 3 were analyzed to elaborate endemic channels (Bortman 1999; Chacón et al. 2021). The geometric mean of historical trends and their confidence inter- vals were calculated using MS Excel (Office 365). The endemic ranges were obtained using the time-series records for at least 5 years before the sampling for this study (in this case, the diarrheal registered cases per week from 2012 to 2019). Mean values and quartiles were used for a defined time unit (epidemiological week). Endemic channels were represented as a graphic elaborated with the R language (version 4.1.1) and ggplot2 package (ver- sion 3.3.1). The registered cases of diarrheal diseases were indicated by an epidemiological week of the study year as a bar, and they were compared with historical data on this disease for the last quinquennium. The data on diarrheal cases was obtained from the Ministry of Health Area located in each region. In addition, inhabitants’ populations per county during the period under study (2012–2019) were obtained from the official demographic estimates published by the National Institute of Statistics and Census of Costa Rica (INEC; https://www.inec.cr/anuario-estadistico). A relationship between somatic coliphages and diarrheal cases was established using multiple logistic regression. This test analyzed the relationship between the presence and absence of each microorganism and the variables such as climatic regions, seasonality, categorized rainfall, and the prevalence of the diarrheal disease. The glm (generalized linear models) function from the R language was used to perform the logistic regression analysis. The exponential coefficient of the logistic regression was calculated to obtain the odds ratios. A significance p-value of ,0.05 was determined for all analyses. RESULTS Regions and microbiological analysis Regarding the microbial indicators and quality of the water source, as shown in Figure 2, the concentrations of FC and E. coli were higher in sources from the ASADAS from region 3 compared to regions 1 and 2 (p,0.05). Also, according to Figure 3, they were higher in samples from the distribution networks from region 2. Furthermore, the higher concentrations from water sources in region 3 were detected during all the sampling periods, in 2018 (transition or rainy seasons) and 2019 (dry, transition, and rainy seasons). In the case of samples from the distribution network, the higher concentrations of FC and E.coli were detectable from the ASADAS in region 2. om http://iwaponline.com/jwh/article-pdf/20/2/344/1013002/jwh0200344.pdf?guestAccessKey=96e7e231-769a-4283-96cd-579a91d25f30 ntes@ucr.ac.cr 022 https://www.inec.cr/anuario-estadistico https://www.inec.cr/anuario-estadistico Figure 2 | Concentration (mean value) of FCs or FC (MPN/100 ml), E. coli (MPN/100 ml), coliphages (PFU/100 ml), and RV (number of viral copies/100 ml) from surface water sources from the ASADAS from regions 1, 2, and 3. Journal of Water and Health Vol 20 No 2, 349 Downloaded from http by kenia.barrantes@u on 01 March 2022 The presence of somatic coliphages was positive in 18% of samples (53/296; 14 samples from sources and 39 from distri- bution networks) in the three regions, in a range from 1 to 1,464 PFU/100 ml. Of coliphages positive samples, 55% belong to region 1 (7/53 from water sources and 22/53 from the distribution network), 26% belong to region 2 (2/53 from water source and 12/53 from distribution networks), and 19% to region 3 (5/53 from water source and 5/53 from distribution networks). Interestingly, RV was detectable from water sources and distribution networks, as observed in Figures 2 and 3. Its detection was possible in 9% of samples (28/296), with higher concentrations during the dry season (more than 1,000 copies/100 ml) compared to transition (100–1,000 copies/100 ml), mainly in regions 2 and 3. However, it was detected in just one source water sample from region 1. Endemic channels of diarrhea cases Four zones were established for the graphical analysis of the endemic ranges: (1) Success zone: diarrhea cases below the 25th percentile of the weekly prevalence (the prevalence of diarrheal cases registered every week); (2) Safety zone: diarrhea cases between the 25th and 50th percentiles of the weekly prevalence; (3) Alert zone: diarrhea cases between the 50th and 75th percentiles of the weekly prevalence; and (4) Outbreak zone: area above the 75th percentile of the weekly prevalence of diar- rheal disease (Bortman 1999). Figure 4 shows the endemic channels of diarrhea diseases from regions 1, 2, and 3 during 2018 and 2019. According to Figure 4(a), in region 1, most of the weeks during 2018 were considered ‘success’ weeks. Fifty-one weeks in 2018 were suc- cessful or safe. One week was considered an ‘alarm’ week, which means the diarrhea cases were outside the expected normal seasonal range. ://iwaponline.com/jwh/article-pdf/20/2/344/1013002/jwh0200344.pdf?guestAccessKey=96e7e231-769a-4283-96cd-579a91d25f30 cr.ac.cr Figure 3 | Concentration (mean value) of FCs or FC (MPN/100 ml), E. coli (MPN/100 ml), coliphages (PFU/100 ml), and RV (number of viral copies/100 ml) from water distribution networks from the ASADAS from regions 1, 2, and 3. Journal of Water and Health Vol 20 No 2, 350 Downloaded fr by kenia.barra on 01 March 2 Nevertheless, during 2019 (Figure 4(b)), there were more weeks considered ‘alarm’ (16 weeks). In addition, 11 weeks of that year were reported as outbreaks of diarrhea. Thus, the cases of diarrhea registered during 2018 and 2019 were relatively high at the time of sampling and coincided with the RV detection in samples. Regarding the diarrheal cases in region 2 (Figure 4(c) and 4(d)), it was observed that in 2018, only 2 weeks were considered alarm, no outbreak weeks were reported that year. In 2019, 10 weeks were alarm, and again, no outbreaks were reported. For region 3 (Figure 4(e) and 4(f)), in 2018, 31 weeks of the year (close to 8 months) were considered as ‘alarm’ or ‘out- break’. In 2019, 14 weeks were also considered in these categories. The observation is congruent with the results obtained from the analysis of microbiological indicators in samples from this region. A significant relationship between somatic coliphages and the diarrheal cases reported in each region was observed (Table 2). The detection of somatic coliphages in water samples was related to 2.40 (1.21–4.80, p,0.05) times to outbreaks and alert weeks of diarrhea cases in a particular region. DISCUSSION Regions and microbiological analysis In general terms, the values of bacterial fecal indicators, such as FC and E. coli, sampled from the water source and distri- bution networks from the three regions were higher according to local and international recommendations for drinking water (WHO 2017). Interestingly, 19 out of the 20 ASADAS in treated water used chlorination (17/20 ASADAS) or electro- lysis (2/20 ASADAS). Just one ASADA did not treat the water at all. Therefore, the disinfection process must secure a good quality for water used in human activities (WHO 2017). om http://iwaponline.com/jwh/article-pdf/20/2/344/1013002/jwh0200344.pdf?guestAccessKey=96e7e231-769a-4283-96cd-579a91d25f30 ntes@ucr.ac.cr 022 Table 2 | Odds ratio values of the epidemiological classification (success, safety, alert, or outbreak) of diarrhea cases and microbiological parameters Coefficient β Odds ratio 2.50% 97.50% p-value Intercept �0.78 0.46 0.31 0.67 ,0.001 FCs �1.08 0.34 0.05 1.67 0.223 E. coli 0.80 2.22 0.43 16.72 0.371 Somatic coliphages 0.88 2.40 1.21 4.80 0.012 RV �1.80 0.17 0.01 0.88 0.088 Bold values denotes a p-value less than 5% (p-value,0.05). Figure 4 | Endemic channels of diarrhea cases for 2018 (a, c and e) and 2019 (b, d and f), from regions 1, 2, and 3. The number of acute diarrheal cases by epidemiological week is represented with black bars, and sampling weeks are marked with a star symbol. The x-axis denotes the epidemiological weeks, and the y-axis denotes the number of acute diarrheal cases reported. Journal of Water and Health Vol 20 No 2, 351 Downloaded from http by kenia.barrantes@u on 01 March 2022 In this study, RVs and coliphages were detectable from treated water (samples from the distribution networks). Virus inac- tivation in water differs from bacteria. Therefore, if free chlorine is used as a disinfectant for water, a controlled process with stable pH, temperature, and organic matter, among others, is required. Also, it is necessary to include an adequate pretreat- ment process before chlorination, such as conventional filtration, direct filtration, or sand slow filtration (Shin & Sobsey 2008; Chacón Jiménez et al. 2015). According to these results, the presence of RV in water is a significant risk of waterborne disease. Conforming to the WHO, this pathogen and E. coli are responsible for most moderate-to-severe diarrhea cases in low-income countries. In addition, diarrheal disease is the second leading cause of death among children under 5 years old (Anderson & Weber 2004; WHO 2019). Linhares et al. (2012) reported a frequency of 53% of RV-positive stool samples in costarrican children with gastroen- teritis under 3 years of age. Ureña-Castro et al. (2019) reported that RV cases predominated in the country in the first 5 months of the year; meanwhile, the increase of RV hospitalizations coincided with dry and cold weather conditions with a 2-month lag (Ureña-Castro et al. 2019). ://iwaponline.com/jwh/article-pdf/20/2/344/1013002/jwh0200344.pdf?guestAccessKey=96e7e231-769a-4283-96cd-579a91d25f30 cr.ac.cr Journal of Water and Health Vol 20 No 2, 352 Downloaded fr by kenia.barra on 01 March 2 Most clinicians who care for adults do not recognize the infection caused by typical pediatric pathogens such as RV (Anderson & Weber 2004). This pathogen is considered a directly transmitted disease due to its high infectivity. However, environmental pathways have usually been ignored (Anderson & Weber 2004; CDC 2011). Kraay et al. (2018) demonstrated that water could disseminate RV in communities and amplify transmission (Kraay et al. 2018). In the model described, RV persists in water sources. Its surface water concentrations and infection incidence vary with temperature, being a critical find- ing regarding the transmission of this disease (Kraay et al. 2018). Interestingly, RV concentration was increased in regions 2 and 3 during the dry season. Ureña-Castro et al. (2019) already reported that RV cases in Costa Rica predominated in the first 5 months of the year, which coincide with the winter months in template regions from the northern hemisphere. Climate variations affect the watersheds and the quality and access of the water. Changes in the availability of water through precipitation, for instance, have a significant impact on the development of communities. Water fluctuations that cause scarcity can lead to the use of inappropriate water sources. Therefore, it will increase the risk of contamination for pathogenic microorganisms causing diarrhea and other waterborne diseases (Duran-Encalada et al. 2017; Morales et al. 2019). RV infections increase during the winter months in template regions, decreasing significantly in summer (Anderson & Weber 2004). However, other authors have indicated that the winter seasonality of RV infections is too simple a generaliz- ation. For example, Levy et al. (2008) observed a correlation between months, temperature, rainfall, relative humidity, and the incidence of RV infections in tropical countries, concluding that RV responds to climate changes, with the highest number of cases found at the colder and drier times of the year. Another critical fact to consider is that universal RV vaccination for infants started in Costa Rica in 2019 and long before in the private medical practice (Ulloa-Gutierrez & Avila-Aguero 2014; Ureña-Castro et al. 2019). Furthermore, several studies have indicated the shedding of pentavalent RV and RV3-BB (human neonatal) vaccine in fecal samples from vaccinated infants (Yen et al. 2011; Cowley et al. 2017). Also, Moresco et al. (2016) reported the thermal and length of stability of one RV vaccine (RotaTeq) in the aquatic environment (Moresco et al. 2016). Therefore, the detection of RV from water in this study cannot be attributable to a wild type or a vaccination RV strain. On the other hand, somatic coliphages comprise groups of bacteriophages that have been suggested as valuable indicators of fecal pollution in water and the behavior and survival of human viruses (Jofre 2009). These viruses replicate in E. coli after infecting it through the cell wall. They may also infect bacterial species closely related to E. coli, such as Klebsiella spp. and Shigella sp. (Jofre et al. 2016). Several studies in Costa Rica have visualized the utility of this indicator related to pathogenic virus circulation in water and a higher frequency of diarrheal cases in a particular community (Chacón Jiménez et al. 2015; Chacón et al. 2020). Endemic channels of diarrhea cases The endemic channel is a graphic representation of the number of diarrhea cases (or any other infectious disease) registered in an area during established periods, in this case, epidemiological weeks in 2018 and 2019. This graphic permits visualization if diarrhea cases in a particular zone increase and start to behave like an outbreak (outbreak is defined as excess cases of one disease in a community, beyond response capabilities). The endemic channel is elaborated using rates of the registered cases of diarrhea for at least the last 5 years. It is used as a criterion to define an infectious disease outbreak from the epidemio- logical point of view (Bortman 1999). This study observed that for regions 1 and 2, no outbreak weeks were observed during 2018 and 2019; meanwhile, for region 3, 12 outbreaks weeks were detected in 2018 and 1 week in 2019. Interestingly, a significant relationship was found between microbiological parameters and climatic regions, including seasonality, categorized rainfall, and diarrheal dis- ease prevalence. No significant differences were observed between regions 1 and 2 for FC and E. coli (p.0.05). However, region 3 showed differences for all fecal indicators (p,0.05). For somatic coliphages, differences were observed between regions 2 and 3 (p,0.05). For RV, no differences were found in any case. Considering the statistical analysis of microbiolo- gical parameters and endemic channels of diarrhea cases, region 3 differed compared to other regions. Somatic coliphages are more related to diarrheal cases than the other microbiological indicators Interestingly, somatic coliphages were more related to the diarrheal cases reported in each region (odds ratio¼2.4; p,0.05). Meanwhile, no significant results were obtained for bacterial fecal parameters (such as FC and E. coli) (Table 1). In addition, Chacón et al. (2020) described a statistically significant relationship between somatic coliphage threshold and the detection of om http://iwaponline.com/jwh/article-pdf/20/2/344/1013002/jwh0200344.pdf?guestAccessKey=96e7e231-769a-4283-96cd-579a91d25f30 ntes@ucr.ac.cr 022 Journal of Water and Health Vol 20 No 2, 353 Downloaded from http by kenia.barrantes@u on 01 March 2022 human enteric viruses such as Norovirus, RV, Enterovirus, and Hepatitis A from wastewater (Chacón et al. 2020). Also, the same authors have reported the presence of somatic coliphages in treated drinking water and their relationship with diarrheal cases in one community from Costa Rica in the absence of bacteriological indicators in the same samples (Chacón Jiménez et al. 2015). Because coliphages are viruses that closely mimic human viral pathogens in terms of environmental dissemination and sur- vival, they are being investigated as a potentially more reliable indicator of fecal pollution and pathogens’ presence in the water. Therefore, the results obtained in the present study confirm that somatic coliphage detection in drinking water is a good indicator of the circulation of diarrhea-causing pathogens and the need to increase the screening of the microbiological quality of the water with other types of indicators of contamination that could correlate better with diarrhea cases. CONCLUSIONS Contaminated water is a source of microbiological pathogens, including bacteria, viruses, protozoa, and helminths, causing waterborne diseases. For example, in Costa Rica, diarrhea outbreaks related to drinking water are reported despite good drinking-water coverage indicators. This study reported RV and microbiological indicators of contamination in drinking water from 20 ASADAS, categorized into three regions (temperate rainy or region 1, tropical rainy or region 2, and tropical rainy and dry or region 3). The micro- biological indicators such as FCs, E. coli, and somatic coliphages were higher in the ASADAS from region 3 than regions 1 and 2. In a few samples in regions 2 and 3, RV was detected with higher concentrations and more frequently in the dry season. No significant results were observed for bacterial indicators and diarrheal cases. However, somatic coliphages, detected in samples from the three regions, showed a statistically significant relationship with diarrheal cases. Furthermore, when coli- phages were detected, it was 2.40 more likely that diarrheal cases were classified as outbreaks or alerts in the region. Thus, this study suggests that somatic coliphages could be a good indicator of diarrhea-causing pathogens in a particular region. ACKNOWLEDGEMENTS The authors thank the Ministry of Health: Dr Tatiana Picado, Dr Laura Sánchez, Dr Sheyris Torres, Dr Glorianela Sancho, Dr Milagro Umaña, Dr María Fuentes, Dr Oscar Bermúdez, Dr Melissa Medina, Dr Viviana García, Dr Asdrúbal Díaz, Dr José Alberto Morales, Dr Danilo Granados, Dr Marvin Quesada, and Dr Alberto Rodríguez for their invaluable help provid- ing us with information regarding diarrheal cases. FUNDING This study was funded by Consejo Nacional de Rectores (CONARE), University of Costa Rica, National University, National Technological University, and Tecnológico de Costa Rica. DATA AVAILABILITY STATEMENT Data cannot be made publicly available; readers should contact the corresponding author for details. REFERENCES Anderson, E. J. & Weber, S. G. 2004 Rotavirus infection in adults. The Lancet 4, 91–98. Available from: https://www.ncbi.nlm.nih.gov/pmc/ articles/PMC7135863/pdf/main.pdf. APHA/AWWA/WEF 2017 Standard Methods for the Examination of Water and Wastewater, 23rd edn. American Public Health Association/ American Water Works Association/Water Environment Federation, Washington, DC, USA. Bortman, M. 1999 Elaboración de corredores o canales endémicos mediante planillas de cálculo. (Elaboration of endemic corridors or channels using spreadsheets). Revista Panamericana de Salud Pública 5 (1), 1–8. https://doi.org/10.1590/s1020-49891999000100001. Bower, K. M. 2014 Water supply and sanitation of Costa Rica. Environmental Earth Sciences 71 (1), 107–123. https://doi.org/10.1007/ s12665-013-2416-x. Centers for Disease Control and Prevention (CDC) 2011 From the Field: Outbreaks of Rotavirus Gastroenteritis Among Elderly Adults in Two Retirement Communities – Illinois, 2011. Available from: https://www.cdc.gov/mmwr/preview/mmwrhtml/mm6042a4.htm. Chacón, L., Barrantes, K., Santamaría-Ulloa, C., Solano, M., Reyes, L., Taylor, L., Valiente, C., Symonds, E. & Achí, R. 2020 A somatic coliphage threshold approach to improve the management of activated sludge wastewater treatment plant effluents in resource-limited regions. Applied and Environmental Microbiology 86 (17), 1–16. https://doi.org/10.1128/AEM.00616-20. ://iwaponline.com/jwh/article-pdf/20/2/344/1013002/jwh0200344.pdf?guestAccessKey=96e7e231-769a-4283-96cd-579a91d25f30 cr.ac.cr http://dx.doi.org/10.1016/S1473-3099(04)00928-4 https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7135863/pdf/main.pdf https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7135863/pdf/main.pdf https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7135863/pdf/main.pdf http://dx.doi.org/10.1590/S1020-49891999000100001 http://dx.doi.org/10.1590/S1020-49891999000100001 http://dx.doi.org/10.1007/s12665-013-2416-x https://www.cdc.gov/mmwr/preview/mmwrhtml/mm6042a4.htm https://www.cdc.gov/mmwr/preview/mmwrhtml/mm6042a4.htm http://dx.doi.org/10.1128/AEM.00616-20 http://dx.doi.org/10.1128/AEM.00616-20 http://dx.doi.org/10.1128/AEM.00616-20 Journal of Water and Health Vol 20 No 2, 354 Downloaded fr by kenia.barra on 01 March 2 Chacón, L., Morales, E., Valiente, C., Reyes, L. & Barrantes, K. 2021 Wastewater-based epidemiology of enteric viruses and surveillance of acute gastrointestinal illness outbreaks in a resource-limited region. American Journal of Tropical Medicine and Hygiene 105 (4), 1004–1012. https://doi.org/10.4269/ajtmh.21-0050. Chacón Jiménez, L., Solano Barquero, M., Barrantes Jiménez, K., Valiente Álvarez, C., Mora Alavarado, D., Reyes Lizano, L. & Achí Araya, R. 2015 Relación entre la presencia de colifagos en agua para consumo humano, las lluvias y las diarreas agudas en Costa Rica. (Relation between coliphages presence in drinking water, rain and acute diarrhea in Costa Rica). Revista Costarricense de Salud Pública 24 (2), 160–167. Cowley, D., Boniface, K., Bogdanovic-Sakran, N., Kirkwood, C. D. & Bines, J. E. 2017 Rotavirus shedding following administration of RV3- BB human neonatal rotavirus vaccine. Human Vaccines and Immunotherapeutics 13 (8), 1908–1915. https://doi.org/10.1080/ 21645515.2017.1323591. De Deus, D. R., Teixeira, D. M., Dos Santos Alves, J. C., Smith, V. C., Da Silva Bandeira, R., Siqueira, J. A. M., De Sá Morais, L. L. C., Resque, H. R. & & Gabbay, Y. B. 2019 Occurrence of norovirus genogroups I and II in recreational water from four beaches in Belém city, Brazilian Amazon region. Journal of Water and Health 17 (3), 442–454. https://doi.org/10.2166/wh.2019.304. Dobbin, K. B. 2013 Co-Management and the Fight for Rural Water Justice: Learning from Costa Rican ASADAS. Pitzer Senior Theses, Paper 33, 109. Available from: http://scholarship.claremont.edu/pitzer_theses/33. Duran-Encalada, J. A., Paucar-Caceres, A., Bandala, E. R. & Wright, G. H. 2017 The impact of global climate change on water quantity and quality: a system dynamics approach to the US–Mexican transborder region. European Journal of Operational Research 256, 567–581. https://doi.org/10.1016/j.ejor.2016.06.016. Fong, T. T. & Lipp, E. K. 2005 Enteric viruses of humans and animals in aquatic environments: health risks, detection, and potential water quality assessment tools. Microbiology and Molecular Biology Reviews 69 (2), 357–371. https://doi.org/10.1128/MMBR.69.2.357–371. 2005. Helsel, D. R. 2012 Statistics for Censored Environmental Data Using Minitab® and R. John Wiley & Sons, Inc., Hoboken, NJ, USA. https:// doi.org/10.1002/9781118162729. Holdridge, L. R. 1971 Life Zone Ecology. Tropical Science Center, 206. Available from: http://reddcr.go.cr/sites/default/files/centro-de- documentacion/holdridge_1966__life_zone_ecology.pdf%0Ahttp://www.cabdirect.org/abstracts/19670604180.html; jsessionid=E82760BE633600B0FC2F82BA1E4D71E2. Instituto Nacional de Estadística y Censos 2021 Objetivos del Desarrollo Sostenible (Sustainable Development Goals). Available from www.inec.cr website: www.inec.cr/objetivos-de-desarrollo-sostenible (accessed 19 April 2021). Jofre, J. 2009 Is the replication of somatic coliphages in water environments significant? Journal of Applied Microbiology 106 (4), 1059–1069. https://doi.org/10.1111/j.1365-2672.2008.03957.x. Jofre, J., Lucena, F., Blanch, A. R. & Muniesa, M. 2016 Coliphages as model organisms in the characterization and management of water resources. Water (Switzerland) 8 (5), 1–21. https://doi.org/10.3390/w8050199. Kraay, A. N. M., Brouwer, A. F., Lin, N., Collender, P. A., Remais, J. V. & Eisenberg, J. N. S. 2018 Modeling environmentally mediated rotavirus transmission: the role of temperature and hydrologic factors. Proceedings of the National Academy of Sciences 115 (12), E2782–E2790. https://doi.org/10.1073/pnas.1719579115. Laboratorio Nacional de Aguas (LNA-AyA) 2021 Cobertura y calidad del agua para uso y consumo humano suministrada por Municipalidades y la ESPH en el 2020 (Coverage and Quality of Water for Human use and Consumption Supplied for Municipalities and ESPH in 2020). Available from: https://www.aya.go.cr/Noticias/Documents/Informe%20cobertura%20agua%20potable%20y% 20saneamiento%202020%20-%20Laboratorio%20Nacional%20de%20Aguas.pdf. Levy, K., Hubbard, A. E. & Eisenberg, J. N. S. 2008 Seasonality of rotavirus disease in the tropics: a systematic review and meta-analysis. International Journal of Epidemiology 38 (6), 1487–1496. https://doi.org/10.1093/ije/dyn260. Linhares, A. C., Macias-Parra, M., Sáez-Llorens, X., Vergara, R., Jimenez, E., Raúl Velázquez, F., Cervantes, Y., Abate, H. J., Rivera, L., Ruttimann, R., Rivera-Medina, D. M., Salinas, B., Ortega-Barria, E., Rubio, P. & Breuer, T. B. 2012 Rotavirus gastroenteritis in Latin America: a hospital-based study in children under 3 years of age. Trials in Vaccinology 1, 36–41. https://doi.org/10.1016/j.trivac. 2012.07.002. MINAET 2013 ‘Agenda del Agua’ Costa Rica 2013–2030 (Water Agenda, Costa Rica 2013–2030) (Issue 2008). Ministerio de Ambiente, Energía y Telecomunicaciones. Available from: http://www.inec.go.cr/A/MT/MedioAmbiente/Publicación/C1/2.1. 4AgendadelAgua2013(DirecióndeaguasMINAE)/GestiónIntegradadelRecursoHídrico.pdf. Morales, E., Solano, M., Morales, R., Reyes, L., Barrantes, K., Achí, R. & Chacón, L. 2019 Evaluation of climatic seasonal influence on the quality of water for human consumption in a San José (Costa Rica) supply system, 2017–2018. Revista Costarricense de Salud Pública 28 (1), 48–58. Available from: https://www.scielo.sa.cr/pdf/rcsp/v28n1/1409-1429-rcsp-28-01-48.pdf. Moresco, V., Damazo, N. A. & Barardi, C. R. M. 2016 Rotavirus vaccine stability in the aquatic environment. Journal of Applied Microbiology 120 (2), 321–328. https://doi.org/10.1111/jam.13021. Ortiz-Malavassi, E. 2014 Atlas digital de Costa Rica. (Digital Atlas of Costa Rica). Escuela de Ingeniería Forestal. Instituto Tecnológico de Costa Rica, Costa Rica. Available from: https://repositoriotec.tec.ac.cr/handle/2238/6749. Prez, V. E., Poma, H. R., Giordano, G. G., Victoria, M., Nates, S. V., Rajal, V. B. & Barril, P. A. 2020 Rotavirus contamination of surface waters from the northwest of Argentina. Journal of Water and Health 18 (3), 409–415. https://doi.org/10.2166/wh.2020.005. om http://iwaponline.com/jwh/article-pdf/20/2/344/1013002/jwh0200344.pdf?guestAccessKey=96e7e231-769a-4283-96cd-579a91d25f30 ntes@ucr.ac.cr 022 http://dx.doi.org/10.4269/ajtmh.21-0050 http://dx.doi.org/10.4269/ajtmh.21-0050 http://dx.doi.org/10.1080/21645515.2017.1323591 http://dx.doi.org/10.1080/21645515.2017.1323591 http://dx.doi.org/10.2166/wh.2019.304 http://dx.doi.org/10.2166/wh.2019.304 http://scholarship.claremont.edu/pitzer_theses/33 http://scholarship.claremont.edu/pitzer_theses/33 http://dx.doi.org/10.1016/j.ejor.2016.06.016 http://dx.doi.org/10.1016/j.ejor.2016.06.016 http://dx.doi.org/10.1128/MMBR.69.2.357-371.2005 http://dx.doi.org/10.1128/MMBR.69.2.357-371.2005 http://reddcr.go.cr/sites/default/files/centro-de-documentacion/holdridge_1966__life_zone_ecology.pdf%0Ahttp://www.cabdirect.org/abstracts/19670604180.html;jsessionid=E82760BE633600B0FC2F82BA1E4D71E2 http://reddcr.go.cr/sites/default/files/centro-de-documentacion/holdridge_1966__life_zone_ecology.pdf%0Ahttp://www.cabdirect.org/abstracts/19670604180.html;jsessionid=E82760BE633600B0FC2F82BA1E4D71E2 http://reddcr.go.cr/sites/default/files/centro-de-documentacion/holdridge_1966__life_zone_ecology.pdf%0Ahttp://www.cabdirect.org/abstracts/19670604180.html;jsessionid=E82760BE633600B0FC2F82BA1E4D71E2 http://reddcr.go.cr/sites/default/files/centro-de-documentacion/holdridge_1966__life_zone_ecology.pdf%0Ahttp://www.cabdirect.org/abstracts/19670604180.html;jsessionid=E82760BE633600B0FC2F82BA1E4D71E2 http://reddcr.go.cr/sites/default/files/centro-de-documentacion/holdridge_1966__life_zone_ecology.pdf%0Ahttp://www.cabdirect.org/abstracts/19670604180.html;jsessionid=E82760BE633600B0FC2F82BA1E4D71E2 http://www.inec.cr/objetivos-de-desarrollo-sostenible http://dx.doi.org/10.1111/j.1365-2672.2008.03957.x http://dx.doi.org/10.3390/w8050199 http://dx.doi.org/10.3390/w8050199 http://dx.doi.org/10.1073/pnas.1719579115 http://dx.doi.org/10.1073/pnas.1719579115 https://www.aya.go.cr/Noticias/Documents/Informe%20cobertura%20agua%20potable%20y%20saneamiento%202020%20-%20Laboratorio%20Nacional%20de%20Aguas.pdf https://www.aya.go.cr/Noticias/Documents/Informe%20cobertura%20agua%20potable%20y%20saneamiento%202020%20-%20Laboratorio%20Nacional%20de%20Aguas.pdf https://www.aya.go.cr/Noticias/Documents/Informe%20cobertura%20agua%20potable%20y%20saneamiento%202020%20-%20Laboratorio%20Nacional%20de%20Aguas.pdf http://dx.doi.org/10.1093/ije/dyn260 http://dx.doi.org/10.1016/j.trivac.2012.07.002 http://dx.doi.org/10.1016/j.trivac.2012.07.002 http://www.inec.go.cr/A/MT/MedioAmbiente/Publicación/C1/2.1.4AgendadelAgua2013(DirecióndeaguasMINAE)/GestiónIntegradadelRecursoHídrico.pdf http://www.inec.go.cr/A/MT/MedioAmbiente/Publicación/C1/2.1.4AgendadelAgua2013(DirecióndeaguasMINAE)/GestiónIntegradadelRecursoHídrico.pdf http://www.inec.go.cr/A/MT/MedioAmbiente/Publicación/C1/2.1.4AgendadelAgua2013(DirecióndeaguasMINAE)/GestiónIntegradadelRecursoHídrico.pdf https://www.scielo.sa.cr/pdf/rcsp/v28n1/1409-1429-rcsp-28-01-48.pdf https://www.scielo.sa.cr/pdf/rcsp/v28n1/1409-1429-rcsp-28-01-48.pdf http://dx.doi.org/10.1111/jam.13021 https://repositoriotec.tec.ac.cr/handle/2238/6749 https://repositoriotec.tec.ac.cr/handle/2238/6749 http://dx.doi.org/10.2166/wh.2020.005 http://dx.doi.org/10.2166/wh.2020.005 Journal of Water and Health Vol 20 No 2, 355 Downloaded from http by kenia.barrantes@u on 01 March 2022 Shin, G. A. & Sobsey, M. D. 2008 Inactivation of norovirus by chlorine disinfection of water. Water Research 42 (17), 4562–4568. https://doi. org/10.1016/j.watres.2008.08.001. Solano Barquero, M., Chacón Jiménez, L. M., Barrantes Jiménez, K. & Achí Araya, R. 2012 Implementación de dos métodos de recuento en placa para la detección de colifagos somáticos, aportes a las metodologías estándar (Implementation of two plate count methods for detection of somatic coliphages and contributions to the standard methodologies). Revista Peruana de Biologia 19 (3), 335–340. https:// doi.org/10.15381/rpb.v19i3.1050. Soto-Córdoba, S. M., Gaviria-Montoya, L. & Pino-Gómez, M. 2016 Situación de la gestión del agua potable en las zonas rurales de la provincia de Cartago, Costa Rica (Situation of drinking water management in the rural zones from the province of Cartago, Costa Rica). Revista Tecnología En Marcha 29 (8), 67. https://doi.org/10.18845/tm.v29i8.2986. Ulloa-Gutierrez, R. & Avila-Aguero, M. L. 2014 Rotavirus vaccination in Central American children. Expert Review of Vaccines 13 (6), 687–690. https://doi.org/10.1586/14760584.2014.905747. UN-Water 2014 Global Analysis and Assessment of Sanitation and Drinking-Water (GLAAS) 2014. pp. 1–12. Available from: https://www. unwater.org/publication_categories/glaas/ (accessed 17 February 2021). Ureña-Castro, K., Ávila, S., Gutierrez, M., Naumova, E. N., Ulloa-Gutierrez, R. & Mora-Guevara, A. 2019 Seasonality of rotavirus hospitalizations at Costa Rica’s National Children’s Hospital in 2010–2015. International Journal of Environmental Research and Public Health 16 (13), 1–13. https://doi.org/10.3390/ijerph16132321. Valiente, C. &Mora, D. 2007 Estudio bacteriológico del agua asociado a brotes de diarrea en Costa Rica, 1999–2005 (Bacteriological study of waterborne diarrheal outbreaks in Costa Rica, 1999–2005). AIDIS 1 (3), 8. https://doi.org/10.22201/iingen.0718378xe.2007.1.3.14433. World Health Organization 2017 Guidelines for Drinking-Water Quality: Fourth Edition Incorporating the First Addendum (Fourth). Licence: CC BY-NC-SA 3.0 IGO, Geneva. Available from: https://apps.who.int/iris/handle/10665/254637 (accessed 20 April 2021). World Health Organization 2019 Water, Sanitation, Hygiene and Health. A Primer for Health Professionals. Available from: https://apps. who.int/iris/bitstream/handle/10665/330100/WHO-CED-PHE-WSH-19.149-eng.pdf?ua=1 (accessed 20 April 2021). Yen, C., Jakob, K., Esona, M. D., Peckham, X., Rausch, J., Hull, J. J., Whittier, S., Gentsch, J. R. & LaRussa, P. 2011 Detection of fecal shedding of rotavirus vaccine in infants following their first dose of pentavalent rotavirus vaccine. Vaccine 29 (24), 4151–4155. https:// doi.org/10.1016/j.vaccine.2011.03.074. First received 9 September 2021; accepted in revised form 21 January 2022. Available online 3 February 2022 ://iwaponline.com/jwh/article-pdf/20/2/344/1013002/jwh0200344.pdf?guestAccessKey=96e7e231-769a-4283-96cd-579a91d25f30 cr.ac.cr http://dx.doi.org/10.1016/j.watres.2008.08.001 http://dx.doi.org/10.15381/rpb.v19i3.1050 http://dx.doi.org/10.15381/rpb.v19i3.1050 http://dx.doi.org/10.15381/rpb.v19i3.1050 http://dx.doi.org/10.18845/tm.v29i8.2986 http://dx.doi.org/10.18845/tm.v29i8.2986 http://dx.doi.org/10.1586/14760584.2014.905747 https://www.unwater.org/publication_categories/glaas/ https://www.unwater.org/publication_categories/glaas/ https://www.unwater.org/publication_categories/glaas/ http://dx.doi.org/10.3390/ijerph16132321 http://dx.doi.org/10.3390/ijerph16132321 http://dx.doi.org/10.22201/iingen.0718378xe.2007.1.3.14433 http://dx.doi.org/10.22201/iingen.0718378xe.2007.1.3.14433 https://apps.who.int/iris/handle/10665/254637 https://apps.who.int/iris/handle/10665/254637 https://apps.who.int/iris/bitstream/handle/10665/330100/WHO-CED-PHE-WSH-19.149-eng.pdf?ua=1 https://apps.who.int/iris/bitstream/handle/10665/330100/WHO-CED-PHE-WSH-19.149-eng.pdf?ua=1 https://apps.who.int/iris/bitstream/handle/10665/330100/WHO-CED-PHE-WSH-19.149-eng.pdf?ua=1 http://dx.doi.org/10.1016/j.vaccine.2011.03.074 http://dx.doi.org/10.1016/j.vaccine.2011.03.074 Occurrence of pathogenic microorganisms in small drinking-water systems in Costa Rica INTRODUCTION Drinking-water coverage in Costa Rica METHODS Drinking-water sampling Climatic conditions FC and E. coli analysis Somatic coliphages Rotavirus Spatial description and data analysis RESULTS Regions and microbiological analysis Endemic channels of diarrhea cases DISCUSSION Regions and microbiological analysis Endemic channels of diarrhea cases Somatic coliphages are more related to diarrheal cases than the other microbiological indicators CONCLUSIONS ACKNOWLEDGEMENTS FUNDING DATA AVAILABILITY STATEMENT REFERENCES