Aetiological, seasonal and antibiotic susceptibility patterns of diarrhoeal diseases in Bhutan (2016–2022): a retrospective study of surveillance data ======================================================================================================================================================== * Kinley Gyem * Sonam Pelden * Dorji Tshering * Kinley Penjor * Rinzin Wangchuk * Sangay Dorji * Jigme Tenzin * Birdi Lal Phuyel ## Abstract **Objectives** This study aimed to identify the aetiological spectrum, seasonal distribution and antimicrobial resistance patterns of diarrhoeal diseases in Bhutan. **Study design and setting** The study used a cross-sectional, retrospective analysis of secondary data gathered through a passive, hospital-based sentinel surveillance for diarrhoeal disease across 12 hospitals, representing Bhutan’s demographically diverse regions. **Participants** A total of 3429 participants’ data of all age groups who presented with diarrhoea at sentinel hospitals between 1 January 1 2016 and 31 December 2022 were analysed. **Results** Diarrhoeagenic *Escherichia coli* (DEC), *Shigella*, *Salmonella* and *Aeromonas* spp. were predominant bacterial pathogens, while *Rotavirus, Astrovirus* and *Norovirus* were the leading viral pathogens. Coinfections were observed in 195 cases. Children under nine were significantly affected than the other age groups. Seasonal trends revealed that bacterial pathogen incidence peaked during the summer/monsoon season, viral pathogens were more common in winter and spring, and parasites persisted year-round. Among the antibiotics tested, gentamicin, chloramphenicol, ceftriaxone and tetracycline exhibited high efficacy, with susceptibility rates of 93.4%, 87.2%, 81.5% and 69.5%, respectively. Conversely, high resistance rates were observed for amoxicillin (80.3%), ampicillin (77.4%) and nalidixic acid (69.5%). Multidrug resistance was prevalent, with β-lactamase production contributing to resistance rates of 80.7% to penicillin and 65.4% to fluoroquinolones groups. Cephalosporin resistance was also notable, with rates of 34.4% for cephalexin, 40.0% for cefazolin and 16.9% for ceftriaxone. **Conclusions** DEC and *Rotavirus* were identified as the leading causes of diarrhoea, with significant resistance patterns observed in common bacterial isolates. These findings underscore the need for DEC screening in paediatric cases and emphasise the need for sustained antimicrobial resistance surveillance. * EPIDEMIOLOGY * Gastrointestinal infections * Neglected Diseases * Observational Study * PUBLIC HEALTH ### STRENGTHS AND LIMITATIONS OF THIS STUDY * Leverages passive sentinel surveillance to assess the aetiology, seasonality and antimicrobial susceptibility of diarrhoeal diseases. * Provides a hospital-based, real-world perspective on diarrhoeal diseases. * Includes antimicrobial susceptibility, helping to understand resistance trends in enteric pathogens. * Provides potential under-representation of the community cases due to reliance on hospital-based surveillance. ## Background Diarrhoeal disease is a major global health challenge, exerting a substantial burden on public health systems and disproportionately affecting vulnerable populations, particularly children. Despite advances in sanitation and increased public health awareness and prevention efforts, diarrhoeal diseases continue to be a leading cause of morbidity and mortality, predominantly stemming from contaminated food and water sources worldwide. Diarrhoeal disease accounts for nearly 1.7 billion cases annually, resulting in approximately 800 000 fatalities, with a disproportionate impact on children under five in developing regions such as South Asia and Sub-Saharan Africa. In these areas, inadequate sanitation, poor hygiene practices and limited access to clean water contribute to the persistence of these diseases. Viral and bacterial enteropathogens remain the primary causes of acute gastroenteritis, affecting populations in both developed and developing countries.1 2 While diarrhoeal disease manifests across all age groups, the aetiology and clinical trajectory vary depending on age and specific causative agents. Most studies focus predominantly on children under five, resulting in limited data on the aetiology and susceptibility trends among adults. Although diarrhoea accounts for only 2% of deaths in adults, it can still contribute to the spread of enteric infections to other vulnerable populations.3 Furthermore, research on the aetiology of diarrhoeal disease in the Bhutanese population is limited. A deeper understanding of pathogen aetiology, seasonal patterns and antimicrobial susceptibility is essential for the development of targeted preventive measures and effective antimicrobial therapies. In Bhutan, diarrhoea is the leading cause of morbidity and the second most frequently reported disease after respiratory infections, presenting an ongoing public health challenge. Bhutan, a low-middle-income country in South Asia situated in the Eastern Himalayas between China and India, has made considerable advancements in healthcare and sanitation. Nevertheless, in 2022, there were 28 179 reported cases of diarrhoea, including 23 272 cases of acute watery diarrhoea and 2112 cases of acute bloody diarrhoea. Children under 9 years were particularly affected with 4232 cases and 2 mortalities reported that year.4 5 Although diarrhoea-associated mortality remains relatively low, disease morbidity has not shown a significant decline. Rapid detection of pathogenic organisms is crucial for timely patient management and identifying potential outbreak sources. Additionally, the increased resistance to commonly used antimicrobials highlights the need for a comprehensive understanding of local factors contributing to these resistance patterns. Since 2016, the Enteric, Zoonotic and Vector-borne Disease Laboratory within the Royal Center for Disease Control (RCDC) has implemented a hospital-based sentinel surveillance system to monitor diarrhoeal diseases. This system, relying on data from sentinel hospitals, focuses on patients presenting with diarrhoea. Over 6 years, surveillance efforts have aimed to unravel the prevalence and seasonal distribution of major agents responsible for diarrhoeal disease. This study determined the aetiology, seasonal patterns and antimicrobial resistance of pathogens in individuals who sought medical attention and were admitted to these sentinel sites from 2016 to 2022. ## Methods ### Bhutan Diarrhoeal Surveillance and Information System (BDSIS) The National Diarrhoeal Disease Surveillance programme, based at the RCDC, operates as a passive hospital-based sentinel surveillance system, which facilitates the rapid identification of pathogens causing diarrhoea. This programme ensures timely patient management and aids in identifying potential outbreak sources. It provides crucial data for evidence-based decision-making and effective disease control strategies, ultimately improving health outcomes. This study uses previously collected stool specimens and associated data from 12 strategically selected sentinel hospitals: one national referral hospital, two regional referral hospitals and nine district hospitals, as illustrated in figure 1. These hospitals served as critical nodes for a comprehensive data collection network, focusing on individuals with diarrhoeal diseases. The RCDC epidemiology unit, EZVDL, and sentinel hospitals collaborate to maintain and share surveillance data through the in-house-developed BDSIS real-time database system.4 ![Figure 1](http://bmjopen.bmj.com/https://bmjopen.bmj.com/content/bmjopen/15/1/e086332/F1.medium.gif) [Figure 1](http://bmjopen.bmj.com/content/15/1/e086332/F1) Figure 1 Diarrhoea surveillance sentinel site. Participants presenting with diarrhoea, defined as three or more loose stools in 24 hours, were documented using a standardised Case Investigation Form (CIF). Cases linked epidemiologically to outbreaks were excluded from this analysis.6 ### Patient and public involvement Patients and the public were not directly involved in the research process in a formalised or participatory manner. The study relies on secondary data from a hospital-based sentinel diarrhoeal surveillance programme, where the data were collected as part of routine healthcare or surveillance efforts. ### Sample preparation and laboratory methods Fresh stool samples were used to detect the causative agents. Patients who had taken antibiotics prior to sample collection were excluded. Samples from the district hospital without microbiology capacity were delivered to the microbiology referral laboratory in Carry-Blair transport media and cryovials maintaining the cold chain following triple packaging and International Air Transport Association) regulations.7 All specimens were tested for a comprehensive panel of enteric pathogens, including eight bacterial species: *Shigella* spp., *Salmonella* spp., *Campylobacter* spp., *Yersinia enterocolitica*, *Aeromonas* spp., *Enterococcus* spp., *Plesiomonas* spp. and Diarrhoeagenic *Escherichia coli* (DEC), which includes *Enteropathogenic E. coli* (EPEC), *Enterotoxigenic E. coli* (ETEC), *Enterohemorrhagic E. coli* (EHEC), *Enteroaggregative E. coli* (EAEC) and *Enteroinvasive E. coli* (EIEC). Viral pathogens screened included *Norovirus, Rotavirus, Astrovirus* and *Adenovirus,* and two parasitic pathogens, *Giardia lamblia* and *Cryptosporidium*, were also examined. Standard microbiological procedures and techniques were employed throughout the testing. ELISA was used for viral detection, and multiplex PCR assays were employed to pathotype DEC, allowing for the simultaneous detection of multiple *E. coli* strains. Serum agglutination tests were used for subtyping *Shigella* and *Salmonella* spp. On-site microscopic examinations of stool samples were performed to identify intestinal ova and parasites, ensuring a thorough assessment of potential pathogens. Antibiotic susceptibility testing was conducted using the Kirby-Bauer disc diffusion method, adhering to the guidelines established by the Clinical and Laboratory Standards Institute (CLSI).8 Bacterial isolates were tested against a comprehensive panel of antibiotics, including amoxicillin (AMX 30 µg), ampicillin (AMP 10 µg), cephazolin (CZO 30 µg), cephalexin (LEX 30 µg), chloramphenicol (CHL 30 µg), ciprofloxacin (CIP 5 µg), ceftriaxone (CRO 30 µg), nalidixic acid (NAL 30 µg), trimethoprim-sulfamethoxazole (SXT 25 µg), gentamicin (GEN 10 µg) and tetracycline (TCY 30 µg). Susceptibility was interpreted based on the CLSI-defined breakpoints, categorising results as susceptible, intermediate or resistant. Multidrug resistance (MDR) was identified as resistance to at least one antimicrobial agent in three or more different antimicrobial classes, namely β-lactam or β-lactamase inhibitors (penicillin and cephalosporins groups), providing a detailed understanding of the resistance patterns observed in the bacterial isolates.9–11 Coinfection was defined as the presence of two or more pathogens in a single sample, indicating simultaneous infections by multiple pathogens. This approach ensures a comprehensive assessment of both susceptibility and the presence of various pathogens, facilitating more accurate diagnosis and treatment planning. ### Data management and statistical analysis Demographic and laboratory data were systematically extracted from the BDSIS which facilitated the structured and consistent collection of clinical and laboratory information for each patient. The integration with BDSIS enabled seamless data management from extraction to analysis, supporting a comprehensive retrospective review of the data. Incomplete or missing data were flagged and excluded from specific analyses where appropriate.4 Statistical analysis was conducted using the STATA statistical package (V.13.1, StataCorp LP). Descriptive statistics, including frequencies and percentages, were used to summarise demographic characteristics (age, sex), pathogen prevalence, seasonal variations, coinfections and antimicrobial resistance patterns. Stratified analysis was performed to assess pathogen incidence and resistance patterns by age group and season. Two-tailed p values <0.05 were considered statistically significant. Seasons were defined based on their climatic characteristics, dividing the calendar year into four distinct categories: winter (December to February), characterised by colder temperatures; spring (March to May), marked by blossoming flora; summer (June to September), associated with warmer temperatures and monsoon rains; and autumn (October and November), characterised by milder temperatures and changing foliage.12 The seasonality of diarrhoeal pathogens was assessed by calculating the monthly incidence rates and identifying peaks for bacterial and viral infections relative to the seasonal changes. The term ‘suspicious food’ refers to self-reported dietary items by individuals that healthcare professionals deemed potentially linked to the onset of diarrhoea. ### Quality control procedures To ensure the reliability and accuracy of the data and results throughout the study, the stool samples were transported to the referral laboratory under cold chain conditions, following strict packaging protocols to maintain the cold chain and minimise the risk of sample degradation. A standardised CIF was used and the data were collected by the trained clinicians and laboratory officials across the 12 sentinel hospitals involved in the study, ensuring consistent documentation of participant data. Laboratory testing protocols followed established microbiological standards, where all participating laboratories were trained to reduce variability. For antibiotic susceptibility testing, the Kirby-Bauer disc diffusion method was employed, using *E. coli* (ATCC 25922) as a quality control reference to ensure accuracy. A comprehensive range of diagnostic methods was implemented, including ELISA for viral detection, multiplex PCR assays for pathotyping DEC and serum agglutination tests for subtyping *Shigella* and *Salmonella*, which enhanced both sensitivity and specificity in pathogen detection. Demographic and laboratory data were systematically extracted from the BDSIS, facilitating structured data collection. Any incomplete or missing data were flagged and excluded from analyses as appropriate. Data analysis was performed using the STATA statistical package, applying descriptive statistics to assess variations in pathogen incidence and resistance patterns by age group and season. ## Results ### Sociodemographic characteristics A total of 3429 subjects were enrolled since the initiation of the surveillance in 2016. Among these, 1810 (53.0%) were males and 1619 (47.0%) females. Participants’ age spectrum spanned from 1 month to 91 years, with a mean age of 20.4 years. A significant proportion, 1765 (51.5%), were under 9 years old, primarily categorised as dependents or students. Most cases, 2888 (84.2%), received treatment in the outpatient department, while 541 (16.0%) required admission. The sample characteristics predominantly constituted loose stool 2472 (72.0%), while 918 (27.0%) presented with watery stool and 39 (1.0%) had bloody stools. Mucus was positive in 595 (23.0%) samples, and blood was detected in 140 (6.0%) cases, as detailed in table 1. Among all cases, 4.0% (143/3429) of diarrhoea incidents were linked with the consumption of suspected junk food, and 1.5% (53/3429) had a history of travel within or outside the country prior to the onset of illness. View this table: [Table 1](http://bmjopen.bmj.com/content/15/1/e086332/T1) Table 1 Sociodemographic and clinical characteristics of the sample (n=3429) ### Aetiology of enteric bacteria and seasonality Enteric pathogens, including *E. coli,* were identified in 3099 samples, representing a positive rate of 90.4% from a total of 3429 stool samples analysed. Bacterial agents accounted for 550 (16.0%), viral for 490 (14.3%), parasitic for 174 (5.1%) and 195 (5.7%) coinfections. Diarrhoeal disease was predominantly caused by enteric bacteria such as DEC, *Shigella, Salmonella* and *Aeromonas* spp., with positive rates of 8.4%, 4.3%, 3.0% and 2.0%, respectively. *Rotavirus* emerged as the predominant enteric virus pathogen, with a rate of 9.5% (239/2513). *Norovirus* and *Astrovirus* followed as the second aetiological agents, accounting for 5.2% (94/1804) and 4.9% (112/2292) of samples, respectively. Among parasites, *G. lamblia* and *Cryptosporidium* exhibited positive rates of 4.0% (86/2402) and 2.2% (53/2402), as illustrated in figure 2. The DEC type identified included 167 EAEC, 39 ETEC, 24 EPEC, 16 EIEC and 14 EHEC. ![Figure 2](http://bmjopen.bmj.com/https://bmjopen.bmj.com/content/bmjopen/15/1/e086332/F2.medium.gif) [Figure 2](http://bmjopen.bmj.com/content/15/1/e086332/F2) Figure 2 Aetiology of enteric pathogens (n=1214). *E. histolytica*, *Entamoeba histolytica*; EAEC, *Enteroaggregative Escherichia coli*; EHEC, *Enterohemorrhagic Escherichia coli*; EIEC, *Enteroinvasive Escherichia coli*; EPEC, *Enteropathogenic Escherichia coli*; ETEC, *Enterotoxigenic Escherichia coli*; LF *E. coli*, *Lactose fermenter Escherichia coli*; NLF *E. coli*, *Non-lactose fermenter Escherichia coli*; Spp, species; STEC, *Shiga toxin-producing Escherichia coli*. The most prevalent coinfections were virus*-*virus (*rotavirus–astrovirus*) and virus*-*bacteria (*rotavirus*–DEC), each occurring in 65 cases. Additionally, 13 cases of bacteria*-*parasite and 10 cases of virus*-*parasite and parasite*-*parasite coinfections were observed, as detailed in online supplemental table 1. Notably, a significant number of enteric pathogens (n=737) were detected in individuals under the age of 9 years (n=1765), with the agent distribution by age group illustrated in figure 3. ### Supplementary data [[bmjopen-2024-086332supp003.pdf]](pending:yes) ![Figure 3](http://bmjopen.bmj.com/https://bmjopen.bmj.com/content/bmjopen/15/1/e086332/F3.medium.gif) [Figure 3](http://bmjopen.bmj.com/content/15/1/e086332/F3) Figure 3 Agent distribution by age group. The number of participants in each age group is as follows: 0–9 years (n=344), 10–19 years (n=16), 20–29 years (n=44), 30–39 years (n=24), 40–49 years (n=10), 50–59 years (n=25) and 60+ years (n=23). DEC, diarrhoeagenic *Escherichia coli*; Spp, species. The examination of enteric pathogens across various months revealed a consistent presence of enteric bacteria causing diarrhoea throughout the year, with the highest prevalence of bacterial agents recorded during the spring and monsoon seasons. DEC, *Shigella* spp. and *Salmonella* spp. were more common in spring, summer and early fall compared with late fall and winter. *Aeromonas* spp. were notably more prevalent from late spring to the end of summer, indicating a seasonal trend. Enteric viruses were predominant from winter through early spring (January to March), while parasites were detected consistently across all seasons, with peaks during the warm and humid summer months (June and July) and the colder months (January). *Rotavirus* and *Astrovirus* had a longer peak period from January to July, whereas *Norovirus* exhibited a less distinct seasonality, clustering around winter and monsoon seasons. Bacterial pathogens were more frequent during the summer, accounting for 232/550 (42.2%) of all cases, while viral pathogens were more common in winter and spring, accounting for 281/490 (57.3%) of all cases. Although enteric pathogens were detected year-round, cases were significantly increased during the monsoon seasons, as illustrated in figure 4. ![Figure 4](http://bmjopen.bmj.com/https://bmjopen.bmj.com/content/bmjopen/15/1/e086332/F4.medium.gif) [Figure 4](http://bmjopen.bmj.com/content/15/1/e086332/F4) Figure 4 Seasonal prevalence of enteric pathogens. Crypto, *Cryptosporidium*; DEC, diarrhoeagenic *Escherichia coli;* Spp, species. ### Antibiotic susceptibility patterns The susceptibility of isolated enteric bacterial pathogens was assessed against the 11 most commonly used antibiotics. GEN, CHL, CRO and TCY emerged as the most effective antibiotics, displaying overall pathogen susceptibility rates of 93.4%, 87.2%, 81.5% and 69.5%, respectively. Conversely, high resistance rates were observed for AMX, AMP and NAL, with prevalence rates of 80.3%, 77.4% and 69.5%, as illustrated in figure 5. ![Figure 5](http://bmjopen.bmj.com/https://bmjopen.bmj.com/content/bmjopen/15/1/e086332/F5.medium.gif) [Figure 5](http://bmjopen.bmj.com/content/15/1/e086332/F5) Figure 5 Overall antibiogram of the bacterial agents (n=556). AMX, amoxicillin; AMP, ampicillin; CZO, cephazolin; CRO, ceftriaxone; LEX, cephalexin; CHL, chloramphenicol; CIP, ciprofloxacin; GEN, gentamicin; NAL, nalidixic acid; TCY, tetracycline; SXT, trimethoprim-sulfamethoxazole. *Shigella* spp. exhibited notable susceptibility to various antibiotics, with high susceptibility observed for GEN (99.0%), CRO (84.3%), CHL (78.3%), CZO (77.6%) and LEX (56.7%) with similar patterns in *Salmonella*, *Aeromonas* and DEC. However, CIP demonstrated more than 50.0% intermediate results in DEC and 26.3% in other bacterial pathogens, indicating an emerging resistance to this antibiotic. β-Lactamase production was observed in several isolates, predominantly within DEC, significantly contributing to their MDR profiles and indicating the presence of extended-spectrum β-lactamases (ESBLs). These isolates exhibited resistance to multiple classes of antibiotics, with 80.7% resistant to penicillin and 65.4% to fluoroquinolones. First-generation and third-generation cephalosporin resistance was also noted, with 34.4% LEX, 40.0% CZO and 16.9% CRO. Additionally, some isolates demonstrated resistance to non-β-lactam antibiotics, with 40.0% resistance to SXT and 33.9% to TCY. Detailed antimicrobial susceptibility patterns for the different bacterial species are presented in online supplemental figure 1A and 1B. ### Supplementary data [[bmjopen-2024-086332supp001.pdf]](pending:yes) ### Supplementary data [[bmjopen-2024-086332supp002.pdf]](pending:yes) ## Discussion This surveillance study offers valuable insights into the epidemiology of diarrhoea-causing pathogens, identifying the most prevalent pathogens across all age groups, observing that children under nine are disproportionately affected, with males showing a slightly higher prevalence (57.0% compared with 43.0%). Statistical analysis revealed a strong correlation (p<0.001) between pathogen prevalence, gender and age group, underscoring a particularly high risk among young children. This finding highlights the increased risk of diarrhoeal diseases among vulnerable groups, particularly males, and underscores the need for targeted public health interventions.13–15 The study also highlights distinct seasonal patterns in diarrhoeal disease, offering valuable data that could enhance the ability to predict and manage outbreaks. Additionally, the study provides crucial data on the antimicrobial susceptibility patterns of these pathogens, highlighting emerging resistance profiles that were not the focus of earlier studies, which concentrated on case-control data or hospitalised patients.16–19 The findings have significant implications for public health, particularly in the context of managing diarrhoeal diseases and combating antimicrobial resistance in low-income countries like Bhutan. In such settings, diarrhoea remains a widespread issue caused by a growing range of enteric pathogens. Accurately determining the aetiology is essential for effective treatment and improving patient outcomes.1 *Rotavirus* and DEC collectively represent the primary causes of moderate-to-severe diarrhoea in low-income countries.20 In our study, EAEC, *Shigella* spp., *Salmonella* spp. and *Rotavirus* were identified as the most prevalent pathogens associated with diarrhoeal disease, highlighting their significant role in the local burden of gastrointestinal infections. These pathogens exhibit similar and variable prevalence rates across different countries, as reported in the literature.1 21–23 Notably, over 16.0% (550/3429) of stool samples in our study contained established diarrhoeagenic pathogens, with DEC showing a high prevalence, especially in cases of childhood diarrhoea, a significant health issue in developing nations. Our findings align with those from comparable studies.24–26 Besides its prevalence in low-income settings, DEC is also emerging as antimicrobial-resistant enteropathogens in developed nations.2 Previous studies have highlighted the emergence of DEC accounting for up to 25.0% of all diarrhoeal diseases in developing countries.21 24 While our findings are consistent with this trend, contrasting data from a study in China identified *Norovirus* as the most common pathogen,3 and research in Bangladesh highlighted *Campylobacter, Shigella* and *Vibrio cholera* as major bacterial pathogens linked to diarrhoea.24 25 This variability emphasises the importance of continuous surveillance to track emerging pathogens and monitor seasonal trends specific to each region. The high prevalence of DEC and *Rotavirus* in Bhutan highlights the need for targeted public health measures, such as routine screening, especially in vulnerable paediatric populations. Identifying these pathogens enables more effective disease management strategies, including vaccination, improved sanitation and timely clinical interventions. *Rotavirus* is associated with substantial hospitalisations and paediatric mortality as established in previous studies.1 27 It emerged as the second most prevalent diarrhoeal pathogen in our study, which is consistent with prior research.20 Despite being widely recognised as the primary cause of severe diarrhoea in children,28 29 its burden and impact on adults have often been underestimated.30 However, our findings challenge this perspective as 6.2% (71/1151) of the adults in our surveillance cohort tested positive for *Rotavirus*. This rate of infection in adults suggests that they may play a potential role in the transmission dynamics and propagation of *Rotavirus* outbreak, especially in households or areas with close contact between age groups, and requires further exploration as reported.31 Our findings revealed numerous coinfections in 16.0% (195/1124) of the positive cases, with virus+bacteria and virus+virus combinations being the most common. Similar results have been observed in other studies that coinfections with enteric pathogens can intensify diarrhoeal symptoms, increasing the frequency of episodes and raising the risk of hospitalisation, and in severe cases, fatalities.23 32 These findings underscore the need for vigilant monitoring and management of coinfections, as they pose a higher risk for severe disease outcomes. The study also reveals distinct seasonal patterns varied among detected pathogens, with viral pathogens most prevalent during colder months, particularly peaking in January through March, while bacterial pathogens were more frequent in warmer months (spring and summer), and parasitic pathogens observed throughout the year with peaks in June, July and January. These patterns are generally consistent with previously documented seasonality,1 27 though some studies report high *Rotavirus* incidence in dry and hot seasons.1 33 Environmental factors such as temperature, precipitation and humidity likely contributed to these patterns by impacting exposure, host immunity and pathogenicity as reported elsewhere.34 Public health strategies should therefore incorporate environmental and climate-related data to anticipate seasonal outbreaks and implement timely interventions. By leveraging this information, health systems can enhance outbreak preparedness and response efforts, ultimately reducing the burden of diarrhoeal diseases in the population. The antimicrobial susceptibility analysis observed in our study reveals significant resistance among various enteric pathogens. *Shigella* spp. demonstrated high susceptibility to GEN, CRO, CHL and ceftazidime consistent with findings from previous studies.35 36 Conversely, *Salmonella* spp. exhibited high resistance to AMX, AMP and NAL, differing from earlier findings.22 35 37 DEC and *Aeromonas* species showed resistance to AMP, LEX and CZO, which aligns with another research.38 However, CRO, CHL and GEN were the most effective antibiotics in our findings, as similarly reported elsewhere.22 Fluoroquinolone resistance was noted in 50.0% of bacterial pathogens, demonstrating emerging resistance to this critical class of antibiotics. Overall, our report underscores the complex and diverse antimicrobial resistance patterns among enteric pathogens, with 80.7% of bacterial isolates exhibiting ESBL-producing MDR. This trend mirrors a global public health challenge, complicating treatment options, prolonging illness durations, escalating healthcare costs and increasing mortality rates.39 40 The study identified a low incidence of *Campylobacter* within the cohort, which may be attributable to methodological limitations as its detection often relies on culture-based methods, which are known for their limited sensitivity as previously documented by sources.41 42 Inconsistencies in diagnostic practices across the 12 participating hospitals and potential detection inaccuracies likely influenced the detection rates.40 41 Implementing active surveillance and adopting more advanced, sensitive diagnostic methods are essential to gain a more accurate understanding of *Campylobacter* prevalence. Our study has limitations. Supply disruption during the COVID-19 pandemic impacted ELISA testing for enteric viruses, potentially underestimating their prevalence and inconsistencies in the denominator for analysis. Additionally, due to laboratory constraints, *Rotavirus* genotyping was not performed and reliance on hospital-based surveillance may have under-represented cases that did not seek medical attention. ## Conclusions Our study provides a valuable understanding of diarrhoeal diseases in Bhutan, emphasising the significant roles of DEC and *Rotavirus* as the primary contributors to diarrhoea. Children under 9 years of age are particularly vulnerable to bacterial, viral and parasitic infections, a susceptibility shared by the elderly. Social determinants of health, including poverty, limited access to clean water and living conditions, have a notable impact on infection rates, especially in children under nine in low-resource settings. Seasonal patterns were observed with viral pathogens peaking from winter to summer, bacterial pathogens more prevalent in spring and summer, and parasites persisting year-round. Research on coinfections and secondary infections remains crucial to understanding these complex dynamics. The high prevalence of zoonotic pathogens such as *Salmonella* and *Aeromonas* in diarrhoeal cases suggests that contaminated food and water sources, as well as animal-human transmission, play a role in spreading these infections. Strengthening food safety measures, improving water sanitation and monitoring environmental reservoirs of enteric pathogens are essential to reduce the incidence of diarrhoea. Additionally, the significant presence of β-lactamase-producing MDR pathogens and emerging resistance to fluoroquinolones underline the need for coordinated interventions. Enhanced surveillance of diarrhoeal diseases and antimicrobial resistance is crucial for monitoring resistant enteric bacteria, vital for informing treatment guidelines. Adopting a One Health approach with integrated surveillance of ESBL-producing bacteria across human, animal and environmental health sectors could be instrumental in tackling this growing public health threat. Our findings emphasise the importance of including DEC in routine screening protocols for childhood diarrhoea, given its rising trends in MDR as it evolves from a gut flora into a significant health threat. Continued monitoring of *Rotavirus* serotypes is necessary to track genotypic trends critical for effective disease management. Moreover, the integration of advanced molecular tools, such as next-generation sequencing, is vital for detecting novel antimicrobial resistance mechanisms and monitoring the emergence of MDR pathogens. ## Data availability statement Data are available upon reasonable request. Data may be obtained from a third party and are not publicly available. All data relevant to the study are included in the article or uploaded as supplementary information. Not applicable. ## Ethics statements ### Patient consent for publication Not applicable. ### Ethics approval The study protocol was approved by the Research Ethics Board of Health, Ministry of Health, Thimphu, Bhutan (REBH/Approval/2023/029). ## Acknowledgments The authors thank the sentinel sites for their unwavering support of the diarrhoeal surveillance activities. ## Footnotes * Contributors KG led the conception, design, data acquisition, analysis and interpretation, drafted the manuscript, revised it critically for intellectual content and has approved the final version for publication. SP contributed to the study’s conception and design and reviewed the final manuscript draft. DT, KP, SD, RW, JT and BLP supported data management and critically reviewed the manuscript. All authors contributed to the data acquisition and have read and approved the final version. All authors agree that KG will be the corresponding author for this work. She has the authority to make necessary revisions per the journal’s requirements and will handle all correspondence. As the guarantor for this manuscript, she assumes full responsibility for the data’s integrity and the analysis’s accuracy. * Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not-for-profit sectors. * Map disclaimer The inclusion of any map (including the depiction of any boundaries therein), or of any geographic or locational reference, does not imply the expression of any opinion whatsoever on the part of *BMJ* concerning the legal status of any country, territory, jurisdiction or area or of its authorities. Any such expression remains solely that of the relevant source and is not endorsed by *BMJ*. Maps are provided without any warranty of any kind, either express or implied. * Competing interests None declared. * Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research. * Provenance and peer review Not commissioned; externally peer reviewed. * Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise. [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/) This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/). ## References 1. Bicer S , Col D , Erdag GC , et al . A retrospective analysis of acute gastroenteritis agents in children admitted to a university hospital pediatric emergency unit. Jundishapur J Microbiol 2014;7:e9148. [doi:10.5812/jjm.9148](http://dx.doi.org/10.5812/jjm.9148) 2. Afum T , Asandem DA , Asare P , et al . Diarrhea-Causing Bacteria and Their Antibiotic Resistance Patterns Among Diarrhea Patients From Ghana. Front Microbiol 2022;13:894319. [doi:10.3389/fmicb.2022.894319](http://dx.doi.org/10.3389/fmicb.2022.894319) 3. Gong X-H , Wu H-Y , Li J , et al . Epidemiology, aetiology and seasonality of infectious diarrhoea in adult outpatients through active surveillance in Shanghai, China, 2012–2016: a cross-sectional study. BMJ Open 2018;8:e019699. [doi:10.1136/bmjopen-2017-019699](http://dx.doi.org/10.1136/bmjopen-2017-019699) 4. RCDC surveillance system. Available: [http://192.168.100.25/COMMONNEW/dashboard/dashboard.php](http://192.168.100.25/COMMONNEW/dashboard/dashboard.php) [Accessed 08 Mar 2023]. 5. Ministry of health. 2022. Available: [www.health.gov.bt](http://www.health.gov.bt) 6. Diarrhoeal disease. Available: [https://www.who.int/news-room/fact-sheets/detail/diarrhoeal-disease](https://www.who.int/news-room/fact-sheets/detail/diarrhoeal-disease) [Accessed 14 Mar 2024]. 7. IATA. Dangerous goods regulations (DGR). Available: [https://www.iata.org/en/publications/dgr/](https://www.iata.org/en/publications/dgr/) [Accessed 08 Mar 2023]. 8. Methods for dilution antimicrobial susceptibility tests for bacteria that grow aerobically; Approved standard-ninth edition. n.d. Available: [www.clsi.org](http://www.clsi.org) 9. Huang I-F , Lee W-Y , Wang J-L , et al . Fecal carriage of multidrug-resistant Escherichia coli by community children in southern Taiwan. BMC Gastroenterol 2018;18:86. [doi:10.1186/s12876-018-0807-x](http://dx.doi.org/10.1186/s12876-018-0807-x) 10. Tanwar J , Das S , Fatima Z , et al . Multidrug Resistance: An Emerging Crisis. Hindawi Limited, 2014. 11. Magiorakos A-P , Srinivasan A , Carey RB , et al . Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin Microbiol Infect 2012;18:268–81. [doi:10.1111/j.1469-0691.2011.03570.x](http://dx.doi.org/10.1111/j.1469-0691.2011.03570.x) [CrossRef](http://bmjopen.bmj.com/lookup/external-ref?access_num=10.1111/j.1469-0691.2011.03570.x&link_type=DOI) [PubMed](http://bmjopen.bmj.com/lookup/external-ref?access_num=21793988&link_type=MED&atom=%2Fbmjopen%2F15%2F1%2Fe086332.atom) 12. National center for hydrology and meteorology. Available: [https://www.nchm.gov.bt/home/pageMenu/781](https://www.nchm.gov.bt/home/pageMenu/781) [Accessed 16 Jan 2024]. 13. Kombat MY , Kushitor SB , Sutherland EK , et al . Prevalence and predictors of diarrhea among children under five in Ghana. BMC Public Health 2024;24:154. [doi:10.1186/s12889-023-17575-7](http://dx.doi.org/10.1186/s12889-023-17575-7) 14. Sevilimedu V , Pressley KD , Snook KR , et al . Gender-based differences in water, sanitation and hygiene-related diarrheal disease and helminthic infections: a systematic review and meta-analysis. Trans R Soc Trop Med Hyg 2016;110:637–48. [doi:10.1093/trstmh/trw080](http://dx.doi.org/10.1093/trstmh/trw080) [PubMed](http://bmjopen.bmj.com/lookup/external-ref?access_num=28115686&link_type=MED&atom=%2Fbmjopen%2F15%2F1%2Fe086332.atom) 15. Azanaw J , Malede A , Yalew HF , et al . Determinants of diarrhoeal diseases among under-five children in Africa (2013-2023): a comprehensive systematic review highlighting geographic variances, socioeconomic influences, and environmental factors. BMC Public Health 2024;24:2399. [doi:10.1186/s12889-024-19962-0](http://dx.doi.org/10.1186/s12889-024-19962-0) 16. Levine MM , Nasrin D , Acácio S , et al . Diarrhoeal disease and subsequent risk of death in infants and children residing in low-income and middle-income countries: analysis of the GEMS case-control study and 12-month GEMS-1A follow-on study. Lancet Glob Health 2020;8:e204–14. [doi:10.1016/S2214-109X(19)30541-8](http://dx.doi.org/10.1016/S2214-109X(19)30541-8) [CrossRef](http://bmjopen.bmj.com/lookup/external-ref?access_num=10.1016/S2214-109X(19)30541-8&link_type=DOI) 17. Polage CR , Solnick JV , Cohen SH . Nosocomial diarrhea: evaluation and treatment of causes other than Clostridium difficile. Clin Infect Dis 2012;55:982–9. [doi:10.1093/cid/cis551](http://dx.doi.org/10.1093/cid/cis551) [CrossRef](http://bmjopen.bmj.com/lookup/external-ref?access_num=10.1093/cid/cis551&link_type=DOI) [PubMed](http://bmjopen.bmj.com/lookup/external-ref?access_num=22700831&link_type=MED&atom=%2Fbmjopen%2F15%2F1%2Fe086332.atom) 18. Farfán-García AE , Zhang C , Imdad A , et al . Case-Control Pilot Study on Acute Diarrheal Disease in a Geographically Defined Pediatric Population in a Middle Income Country. Int J Pediatr 2017;2017:6357597. [doi:10.1155/2017/6357597](http://dx.doi.org/10.1155/2017/6357597) 19. Jabak SJ , Kawam L , El Mokahal A , et al . Management of acute diarrhea in the emergency department of a tertiary care university medical center. J Int Med Res 2022;50:3000605221115385. [doi:10.1177/03000605221115385](http://dx.doi.org/10.1177/03000605221115385) [PubMed](http://bmjopen.bmj.com/lookup/external-ref?access_num=35942579&link_type=MED&atom=%2Fbmjopen%2F15%2F1%2Fe086332.atom) 20. Wangchuk S , Dorji T , Tshering K , et al . A Prospective Hospital-based Surveillance to Estimate Rotavirus Disease Burden in Bhutanese Children under 5 Years of Age. Trop Med Health 2015;43:63–8. [doi:10.2149/tmh.2014-22](http://dx.doi.org/10.2149/tmh.2014-22) [PubMed](http://bmjopen.bmj.com/lookup/external-ref?access_num=25859154&link_type=MED&atom=%2Fbmjopen%2F15%2F1%2Fe086332.atom) 21. Vu Nguyen T , Le Van P , Le Huy C , et al . Etiology and epidemiology of diarrhea in children in Hanoi, Vietnam. Int J Infect Dis 2006;10:298–308. [doi:10.1016/j.ijid.2005.05.009](http://dx.doi.org/10.1016/j.ijid.2005.05.009) [CrossRef](http://bmjopen.bmj.com/lookup/external-ref?access_num=10.1016/j.ijid.2005.05.009&link_type=DOI) [PubMed](http://bmjopen.bmj.com/lookup/external-ref?access_num=16458564&link_type=MED&atom=%2Fbmjopen%2F15%2F1%2Fe086332.atom) [Web of Science](http://bmjopen.bmj.com/lookup/external-ref?access_num=000239341400010&link_type=ISI) 22. Getie M , Abebe W , Tessema B . Prevalence of enteric bacteria and their antimicrobial susceptibility patterns among food handlers in Gondar town, Northwest Ethiopia. Antimicrob Resist Infect Control 2019;8:111. [doi:10.1186/s13756-019-0566-7](http://dx.doi.org/10.1186/s13756-019-0566-7) 23. Potgieter N , Heine L , Ngandu JPK , et al . High Burden of Co-Infection with Multiple Enteric Pathogens in Children Suffering with Diarrhoea from Rural and Peri-Urban Communities in South Africa. Pathogens 2023;12:315. [doi:10.3390/pathogens12020315](http://dx.doi.org/10.3390/pathogens12020315) 24. Shrivastava AK , Kumar S , Mohakud NK , et al . Multiple etiologies of infectious diarrhea and concurrent infections in a pediatric outpatient-based screening study in Odisha, India. Gut Pathog 2017;9:16. [doi:10.1186/s13099-017-0166-0](http://dx.doi.org/10.1186/s13099-017-0166-0) 25. Albert MJ , Faruque ASG , Faruque SM , et al . Case-control study of enteropathogens associated with childhood diarrhea in Dhaka, Bangladesh. J Clin Microbiol 1999;37:3458–64. [doi:10.1128/JCM.37.11.3458-3464.1999](http://dx.doi.org/10.1128/JCM.37.11.3458-3464.1999) [Abstract/FREE Full Text](http://bmjopen.bmj.com/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjEwOiIzNy8xMS8zNDU4IjtzOjQ6ImF0b20iO3M6MjY6Ii9ibWpvcGVuLzE1LzEvZTA4NjMzMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 26. Vilchez S , Reyes D , Paniagua M , et al . Prevalence of diarrhoeagenic Escherichia coli in children from León, Nicaragua. J Med Microbiol 2009;58:630–7. [doi:10.1099/jmm.0.007369-0](http://dx.doi.org/10.1099/jmm.0.007369-0) [CrossRef](http://bmjopen.bmj.com/lookup/external-ref?access_num=10.1099/jmm.0.007369-0&link_type=DOI) [PubMed](http://bmjopen.bmj.com/lookup/external-ref?access_num=19369525&link_type=MED&atom=%2Fbmjopen%2F15%2F1%2Fe086332.atom) 27. Isenbarger DW , Hien BT , Ha HT , et al . Prospective study of the incidence of diarrhoea and prevalence of bacterial pathogens in a cohort of Vietnamese children along the Red River. Epidemiol Infect 2001;127:229–36. [doi:10.1017/S0950268801005933](http://dx.doi.org/10.1017/S0950268801005933) [PubMed](http://bmjopen.bmj.com/lookup/external-ref?access_num=11693500&link_type=MED&atom=%2Fbmjopen%2F15%2F1%2Fe086332.atom) 28. Crawford SE , Ramani S , Tate JE , et al . Rotavirus infection. Nat Rev Dis Primers 2017;3:17083. [doi:10.1038/nrdp.2017.83](http://dx.doi.org/10.1038/nrdp.2017.83) [PubMed](http://bmjopen.bmj.com/lookup/external-ref?access_num=29119972&link_type=MED&atom=%2Fbmjopen%2F15%2F1%2Fe086332.atom) 29. Omatola CA , Olaniran AO . Rotaviruses: From Pathogenesis to Disease Control-A Critical Review. Viruses 2022;14:875. [doi:10.3390/v14050875](http://dx.doi.org/10.3390/v14050875) 30. Dumic I , Nordin T , Jecmenica M , et al . Gastrointestinal Tract Disorders in Older Age. Can J Gastroenterol Hepatol 2019;2019:6757524. [doi:10.1155/2019/6757524](http://dx.doi.org/10.1155/2019/6757524) 31. For personal use. Only reproduce with permission from the Lancet. n.d. Available: [http://infection.thelancet.com](http://infection.thelancet.com) 32. Zhang S-X , Zhou Y-M , Xu W , et al . Impact of co-infections with enteric pathogens on children suffering from acute diarrhea in southwest China. Infect Dis Poverty 2016;5:64. [doi:10.1186/s40249-016-0157-2](http://dx.doi.org/10.1186/s40249-016-0157-2) 33. Asadi L , Pourlak T , Ahmadi B , et al . Etiological Agents of Pediatric Diarrhea in Ardebil, Northwestern Iran. Arch Pediatr Infect Dis 2018;In Press:1–5. [doi:10.5812/pedinfect.11771](http://dx.doi.org/10.5812/pedinfect.11771) 34. Chao DL , Roose A , Roh M , et al . The seasonality of diarrheal pathogens: A retrospective study of seven sites over three years. PLoS Negl Trop Dis 2019;13:e0007211. [doi:10.1371/journal.pntd.0007211](http://dx.doi.org/10.1371/journal.pntd.0007211) 35. Mekonnen GK , Mengistie B , Sahilu G , et al . Etiologies of diarrhea and drug susceptibility patterns of bacterial isolates among under-five year children in refugee camps in Gambella Region, Ethiopia: a case control study. BMC Infect Dis 2019;19:1008. [doi:10.1186/s12879-019-4599-6](http://dx.doi.org/10.1186/s12879-019-4599-6) 36. Shrestha SK , Shrestha J , Mason CJ , et al . Etiology of Acute Diarrheal Disease and Antimicrobial Susceptibility Pattern in Children Younger Than 5 Years Old in Nepal. Am J Trop Med Hyg 2023;108:174–80. [doi:10.4269/ajtmh.21-1219](http://dx.doi.org/10.4269/ajtmh.21-1219) [PubMed](http://bmjopen.bmj.com/lookup/external-ref?access_num=36509064&link_type=MED&atom=%2Fbmjopen%2F15%2F1%2Fe086332.atom) 37. Mama M , Alemu G . Prevalence, antimicrobial susceptibility patterns and associated risk factors of Shigella and Salmonella among food handlers in Arba Minch University, South Ethiopia. BMC Infect Dis 2016;16:686. [doi:10.1186/s12879-016-2035-8](http://dx.doi.org/10.1186/s12879-016-2035-8) 38. San Joaquin VH , Scribner RK , Pickett DA , et al . Antimicrobial susceptibility of Aeromonas species isolated from patients with diarrhea. Antimicrob Agents Chemother 1986;30:794–5. [doi:10.1128/AAC.30.5.794](http://dx.doi.org/10.1128/AAC.30.5.794) [Abstract/FREE Full Text](http://bmjopen.bmj.com/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjg6IjMwLzUvNzk0IjtzOjQ6ImF0b20iO3M6MjY6Ii9ibWpvcGVuLzE1LzEvZTA4NjMzMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 39. Husna A , Rahman MM , Badruzzaman ATM , et al . Extended-Spectrum β-Lactamases (ESBL): Challenges and Opportunities. Biomedicines 2023;11:2937. [doi:10.3390/biomedicines11112937](http://dx.doi.org/10.3390/biomedicines11112937) 40. Doi Y , Iovleva A , Bonomo RA . The ecology of extended-spectrum β-lactamases (ESBLs) in the developed world. J Travel Med 2017;24:S44–51. [doi:10.1093/jtm/taw102](http://dx.doi.org/10.1093/jtm/taw102) [CrossRef](http://bmjopen.bmj.com/lookup/external-ref?access_num=10.1093/jtm/taw102&link_type=DOI) [PubMed](http://bmjopen.bmj.com/lookup/external-ref?access_num=28521000&link_type=MED&atom=%2Fbmjopen%2F15%2F1%2Fe086332.atom) 41. Buss JE , Cresse M , Doyle S , et al . Campylobacter culture fails to correctly detect Campylobacter in 30% of positive patient stool specimens compared to non-cultural methods. Eur J Clin Microbiol Infect Dis 2019;38:1087–93. [doi:10.1007/s10096-019-03499-x](http://dx.doi.org/10.1007/s10096-019-03499-x) [PubMed](http://bmjopen.bmj.com/lookup/external-ref?access_num=30783889&link_type=MED&atom=%2Fbmjopen%2F15%2F1%2Fe086332.atom) 42. Platts-Mills JA , Liu J , Gratz J , et al . Detection of Campylobacter in stool and determination of significance by culture, enzyme immunoassay, and PCR in developing countries. J Clin Microbiol 2014;52:1074–80. [doi:10.1128/JCM.02935-13](http://dx.doi.org/10.1128/JCM.02935-13) [Abstract/FREE Full Text](http://bmjopen.bmj.com/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjk6IjUyLzQvMTA3NCI7czo0OiJhdG9tIjtzOjI2OiIvYm1qb3Blbi8xNS8xL2UwODYzMzIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9)