1. Introduction
In recent years, emerging antimicrobial resistance (AMR) has been recognized as a significant public health concern that transcends international boundaries [
1]. In 2019, it was estimated there were 4.95 million deaths globally associated with bacterial AMR, including 1.27 million deaths directly attributable to bacterial AMR, with the highest mortality currently seen in South Asian and sub-Saharan African countries [
2]. There is also considerable morbidity and costs associated with AMR [
3,
4,
5,
6]. If AMR is not adequately tackled, this could reduce gross domestic product per country by up to 3.8% [
5,
6]. As a result, AMR is increasingly seen as the next pandemic unless urgent measures are introduced to reverse the rising rate [
7]. The overuse and misuse of antibiotics either prophylactically or therapeutically in either humans or animal husbandry generates selection pressure, increasing the development of AMR in pathogens and other diverse commensal microbial populations [
8,
9,
10,
11]. Whilst AMR is a universal phenomenon, the burden among low- and middle-income countries (LMICs) is appreciably higher due to economic, political, and environmental factors, including poor governance and infrastructures, as well as a limited number of national initiatives [
12,
13,
14,
15]. This is now changing with increasing recognition of the clinical and economic consequences of AMR. Global initiatives include the Global Action Plan (GAP) by the World Health Organization (WHO) to reduce AMR [
1], as well as initiatives from the World Bank and OECD [
6,
16]. The GAP has been translated into the National Action Plans in Bangladesh, with countries at different stages of their implementation due to resource and other issues [
17,
18,
19,
20].
Other important global initiatives include dividing antibiotics into three different categories based on their resistance potential, which includes the ‘Access’, ‘Watch’, and ‘Reserve’ categories [
21,
22]. ‘Access’ antibiotics should typically be prescribed to treat commonly encountered infections with lower resistance rates and include penicillins. Those in the ‘Watch’ group, which includes azithromycin, quinolones, and a number of cephalosporins, should ideally only be prescribed in critical conditions as they have a greater chance of resistance development. Antibiotics in the ‘Reserve’ category, which include fifth-generation cephalosporins, some carbapenems, and linezolid, should only prescribed in multidrug resistance cases, with the aim of curbing rising AMR rates [
21,
23,
24,
25]. Reducing the extent of ‘Watch’ antibiotics being prescribed or dispensed can help to appreciably reduce the extent of multidrug-resistant organisms [
26]. More recently, the WHO has launched the AWaRe book, giving prescribing guidance for approximately 35 clinical infections, including the choice of drug [
24,
27], dose, and duration for both children and adults, to reduce the extent of inappropriate antibiotic use. The target is that at least 60% of antibiotic use in given settings should be ‘Access’ antibiotics [
23,
24].
Colistin is also classified as a ‘Reserve’ antibiotic following the global increase in the prevalence of carbapenem-resistant
Enterobacteriaceae, [
22,
28]. However, there are now concerns with its over use and resultant resistance development [
29,
30]. This includes in animals where there has been appreciable use of colistin as a growth promoter in recent years [
30], with studies now showing that 96% or more of total worldwide use of colistin is still in poultry and pig farming [
31,
32]. Concerns with the overuse of colistin including both animals and humans, and resultant resistance development via zoonotic gene transfers, coupled with its importance in treating resistant gram-negative infections to reduce morbidity and mortality, has resulted in the World Health Organization and others classifying colistin as an antibiotic of very high importance for use in humans with its use reserved [
21,
33,
34,
35,
36,
37,
38,
39]. Alongside this, many countries now ban the use of colistin as a growth promoter in animal feeds and prophylactically to prevent bacterial infections [
40,
41,
42]. Such measures have shown to be effective in reducing resistant strains [
41,
42,
43]. This is important in Bangladesh, given extensive colistin-resistant
Escherichia coli in broiler meat and chicken feces [
44,
45,
46], exacerbating resistance among patients to colistin in Bangladesh [
47,
48,
49,
50]. Over-the-counter dispensing of antibiotics is also common in Bangladesh and a concern, especially when this involves ‘Watch’ and ‘Reserve’ antibiotics [
51,
52,
53,
54]
Overall, activities to enhance the appropriate use of colistin in both animals and humans are essential as colistin still remains the antibiotic of choice for multiple drug-resistant gram-negative bacterial infections (MDR-GNB). This includes carbapenem-resistant
Acinetobacter baumannii (CRAB) as well as other pathogens resistant to the new antimicrobial agents [
55,
56,
57]. However, the use of colistin as a last resort antibiotic is greatly threatened by the rise of plasmid-borne mobile colistin resistance gene [
58,
59,
60], spreading rapidly via horizontal gene transfer [
61]. Resistance to colistin is generated by the chromosomally mediated modification of lipopolysaccharide (LPS) [
62]. Acquisition of colistin resistance by a novel plasmid-mediated gene,
mcr-1, was first described in
Enterobacteriales from both farm-animal products and humans [
63]. Earlier studies have shown the genotypic linkage of the mobile colistin resistance genes,
mcr-1, to the phenotypic colistin resistance [
64,
65]. Since then, the variants of
mcr-carrying multiple species of
Enterobacteriales have been detected in many countries from environments, animals, and humans [
66,
67,
68,
69]. Subsequently, more variants of transferable colistin resistance
mcr gene (
mcr-1 to
mcr-9) have been described in
Enterobacteriaceae [
70,
71].
In general, a substantial part of
Enterobacteriales is normal flora; however, a few of those microbial pathogens can cause systemic bacteremia [
72], community-acquired infections [
73,
74], and healthcare-associated infections (HAI) [
75]. A recent outbreak with colistin-resistant pathogens in China ended a very high case-fatality rate in humans [
76]. Consequently, the identification of the root cause, transmission, and trajectories of colistin-resistant infection is an increasing priority globally. We are aware that
mcr-gene variants can be detected in the environment, animals, human fecal samples, and food products [
77]. However, only a limited number of studies have also showed the dissemination of plasmids carrying these variants in infants with acute diarrhea [
78]. This is important in Bangladesh, with diarrhea being a major cause of childhood mortality in the country, combined with the increasing prevalence of resistant genes in children exacerbated by the overuse of antibiotics [
79,
80,
81].
Consequently, this study was designed to investigate different variants of the mcr gene and their association with colistin resistance among diarrheal pathogens in infants, children and adults in Bangladesh and the demographic factors associated with the identified mcr variants. The findings can be used to suggest future policies and initiatives where there are concerns.
3. Discussion
We believe this is one of the first studies in Bangladesh to investigate colistin-resistant genes in human diarrheal pathogens among infants and children. Our data clearly showed the association of phenotypic colistin resistance and the
mcr genes (
mcr-1 to
mcr-5), similar to other published studies on colistin resistant bacteria isolated from diarrhea [
78,
84,
85]. This builds on earlier studies in Bangladesh including among children and adults [
68,
86].
Isolates harboring
mcr genes have been detected with high MIC values, showcasing disparities between the agar dilution test and the disk diffusion method [
87]. This variation can be attributed to the slow diffusion of colistin disks on agar medium for the complex and large molecular structure of colistin sulfate. In parallel, some
mcr-negative isolates exhibited resistance to colistin. Several potential reasons may account for this phenomenon. Firstly, mutations in the
mgrB gene [
88], responsible for binding polymyxin antibiotics in the gram-negative cell wall, are particularly prevalent in
Klebsiella pneumoniae [
89,
90]. Secondly, the absence of testing for other variants of
mcr genes including
mcr-6, mcr-7, mcr-8, mcr-9, and
mcr-10, which could also attribute phenotypic colistin resistance [
91,
92]. Thirdly, some resistant bacteria may develop capsules, a polysaccharide coating on the outer surface of the cell wall [
93,
94,
95]. Fourthly, overexpression of efflux pump systems could contribute to resistance development [
96,
97]. Fifthly, modulation in the bacterial cell surface, including alterations in the structure of lipopolysaccharides (LPS) of the cell membrane, might affect the binding of polymyxin antibiotics [
98]. Additional investigations will be necessary to uncover potential molecular explanations for the observed differences between phenotypic and genotypic colistin resistance.
In this study, male children had a higher incidence rate of diarrhea than female children, similar to previous studies [
99,
100]. However, the reason for this difference is unclear. We also found that 57.14% of
mcr-positive isolates were resistant to amoxicillin/clavulanic acid, and 57 to 71% of
mcr-carrying isolates resistant to all generations of cephalosporins, although higher generation were more susceptible. This is a concern with clinicians now prescribing more carbapenems due to the decreased potency of cephalosporins, with meropenem showing more susceptibility among the 17 different antibiotics in the eight groups studied. However, recent studies showed carbapenem-resistant
Enterobacteriaceae are now a global threat [
101,
102,
103]. The MDR status of all
mcr-positive isolates is positive, which is a threat to public health in Bangladesh and beyond, given the ensuing rise in untreatable infectious diseases [
104,
105,
106].
This study underscores the clinical significance of establishing comprehensive surveillance systems for priority antibiotics such as colistin. In addition, we urge the Government of Bangladesh to ban the use of colistin as a growth promoter in animal feeds and prophylactically to prevent bacterial infections, similar to other countries [
40,
41,
42]. This is because such measures have been shown to be effective in reducing resistant strains [
41,
42,
43]. Alongside, this instigates antimicrobial stewardship programs (ASPs) in ambulatory care similar to other important antibiotics in Bangladesh where there are concerns about resistance development [
107,
108,
109,
110]. This includes ASPs among community pharmacists and drug sellers building on ASPs in other sectors in Bangladesh [
52,
111].
We are aware that there are several limitations with this study. Firstly, it was challenging to accurately estimate the real-world scenario since adults experiencing diarrheal problems could readily seek treatment at a healthcare center, whereas children relied on their parents’ decisions and assistance to access medical care. Secondly, the small sample size posed significant obstacles to conducting fully powered statistical analyses. However, our research spanned fifteen distinct districts across Bangladesh, and the outcomes are anticipated to be applicable if similar studies are conducted in other districts in the country. Thirdly, while this study examined mcr gene variants up to mcr-5, newer variants such as mcr-6, mcr-7, mcr-8, and mcr-9 were not investigated. Having said this, efforts were made to maintain internal validity by conducting independent trials when necessary. Despite these limitations, we believe our findings are robust, providing guidance to all key stakeholders in Bangladesh to enhance future sensitivity to colistin.
4. Materials and Methods
4.1. Study Design and Sampling
A prospective cross-sectional study was conducted between January 2020 and December 2020 among diarrheic children visiting the outpatient Department of Uttara Adhunik Medical College Hospital, Dhaka, Bangladesh. A total of 179 children and adults having acute diarrhea participated in this study prior to treatment with any prescribed antibiotics.
Acute diarrhea was defined as three or more liquid, loose, mucus, or bloody stools within 24 h, lasting no longer than 14 days. Fever was defined as a temperature of ≥37.5 °C. Demographic data was taken from each child, and informed consent was obtained from the parents or guardians before sample collection. All relevant demographic, clinical, and laboratory data were recorded and transferred to the questionnaire prepared for this study.
The ages of the children and patients were categorized into five groups, namely <1 year, 1-5 years, 6-10 years, 11-15, and >15 years, based on previous studies [
112]. The income of parents was classified into rich, middle class, or poor, and as we are aware, this can make a difference [
113].
A sterilized cotton swab was dipped in the mucus, purulent or bloody part of the freshly passed stool sample, placed immediately in CaryBlair Medium (Oxoid, Hampshire, UK), and transported to the laboratory for further analysis within six hours of collection.
4.2. Isolation and Identification of Bacteria
Collected samples were pre-enriched in buffered peptone water (Oxoid®, Hampshire, UK) at a dilution ratio of 1:10 and were incubated overnight at 37 ◦C. A loopful of each culture was streaked on MacConkey agar (Liofilchem Inc, Italy) and cysteine-, lactose-, and electrolyte-deficient (CLED) agar (Liofilchem Inc, Italy) and subsequently incubated at 37 ◦C for 24 h in aerobic condition simultaneously. MacConkey agar supports gram-negative diarrheal pathogens (
Supplementary Figure S1A), while CLED agar aids in the growth of gram-negative bacteria and gram-positive cocci if present in diarrheal samples. Colony counts of 10
3 or 10
5 CFU/mL were considered for a cut-off value for a probable diarrheal sample [
114].
Gram’s staining and biochemical tests were initially performed to identify growth-positive bacteria. A rapid biochemical-test kit API 20E (BioMe´rieux, Durham, NC), consisting of carbohydrate batteries and enzymatic substrates in a set of chromogenic panels, was used to verify the isolated identity (
Supplementary Figure S1B) [
115]. A part of the bacterial identity was confirmed by the polymerase chain reaction (PCR) amplification and sequencing of the 16S rDNA gene [
116]. In total, 228 different isolates were generated from all the diarrheal samples. Three bacteria (one Proteus and two staphylococci) were excluded from the study for the next level analyses since colistin resistance is a natural phenomenon of the excluded isolates. The remaining 225 isolates were subjected to assessment of colistin susceptibility and
mcr-1 to
mcr-5 carriage. The isolates were preserved in 30% glycerol in Trypticase Soy Broth (TSB) at – 80°C until further use.
4.3. Phenotypic Colistin Susceptibility Testing
The phenotypic antibiogram profiles of diarrheal isolates against colistin were determined primarily using the Kirby–Bauer disk diffusion method according to the European Committee on Antimicrobial Susceptibility Testing (EUCAST) and Clinical & Laboratory Standards Institute (CLSI) guidelines (
Supplementary Figure S1C) [
117,
118]. A 3-hour bacterial suspension in Mueller–Hinton broth was prepared to a concentration of McFarland 0.5 equivalent and then streaked on Mueller–Hinton agar (MHA, Oxoid, Basingstoke, UK) plates using a cotton swab to ensure consistent growth. The susceptibility of the isolate to the following discs of antibiotics (Oxoid, Basingstoke, UK) was evaluated: colistin (25 µg), Amoxycillin + Clavulanic acid (30 µg), Cefuroxime Sodium (30 µg), Cefixime (30 µg), Cefepime (30 µg), Imipenem (10 µg), Meropenem (10 µg), Nalidixic acid (10 µg), Ciprofloxacin (5 µg), Lomefloxacin (10 µg), levofloxacin (5 µg), Gentamicin (30 µg), Amikacin (30 µg), and Netilmicin (30 µg), Tobramycin (10 µg), Nitrofurantoin(300 µg), Trimethoprim-sulfamethoxazole (25 µg); by placing them on the bacterial lawns and incubating at 37°C overnight. A clear zone was developed around the disc for sensitive bacteria, and zone diameter was measured and evaluated to categorize Bacteria as susceptible (S), intermediate (I), and resistant (R) from the CLSI guideline charts for the appropriate antibiotics tested [
118].
The isolates’ phenotypic antimicrobial susceptibilities were further tested by the agar dilution method [
119]. The lowest concentration of colistin adequate to inhibit the visible growth of bacterial test isolates was determined by minimal inhibitory concentration (MIC) measurement by agar dilution method [
119]. For agar dilution MIC, different concentrations of colistin-sulfate powder (Santa Cruz Biotechnology Inc, TX) from .50 µg/mL to 256.0 µg/mL in a two-fold dilution order were incorporated into MHA medium accordingly. One pure culture colony was inoculated into Mueller-Hinton broth to prepare each test inoculum and incubated for three hours at 37°C that develops a density of inoculum equivalent to 10
4 colony-forming units (CFU) per spot to drop on the MHA. The inoculum density was periodically compared to a 0.5 McFarland standard, equivalent to approximately 10
8 CFU/mL. The plates were incubated at 37°C in the air for 18-20 hours. Agar dilution MICs were performed in duplicates. The experiments were repeated when some single colonies or a thin haze growth was observed within the inoculated spot.
The epsilometer test (E-test) was performed parallel partly using a commercial strip containing a predefined gradient of colistin concentrations (Liofilchem Inc, Italy) to validate colistin MIC determination by the agar dilution method [
87,
120]. Concordant results were found in independent MIC assessment and E-test (
Supplementary Figure S1D,E).
E. coli ATCC 25922 strain was used as the quality control strain for disc diffusion and MIC testing. Besides, a control plate without colistin-sulfate was examined for the growth of both test and control isolates. Following EUCAST and CLSI guidelines, isolates were considered susceptible (S) when the MIC values exhibited ≤2 μg/mL and resistant (R) when MICs appeared >2 μg/mL [
117]. Multidrug-resistant (MDR) isolates were described as those isolates that were found to be resistant to at least three different classes of antibiotics [
121].
4.4. Detection of the Colistin Resistance mcr Genes
All 225 isolates were subjected to a single plex polymerase chain reaction (PCR) to detect the
mcr-1 gene, yielding a 309 bp DNA band, using primers described elsewhere [
63], and confirmed by sequencing. Amplicons were visualized under UV light after 1.2% agarose gel electrophoresis. The other four primer pairs to detect
mcr-2, mcr-3, mcr-4, and
mcr-5 gene amplicons were obtained from a recently published original study [
122]. Multiplex polymerase chain reaction (PCR) was conducted to detect the
mcr-1 to
mcr-5 genes in the isolates. In brief, the modified protocol was as follows: prepared bacterial DNA (2.0 μL) was added to a 2× PCR premixture (15 μL, GeneON, Germany), and five pmol of each primer (1 μL), and deionized water was added to obtain a final volume of 30 μL. Reactions went through an initial denaturation at 94 °C for 15 min followed by 25 cycles of amplification (Applied Biosystems 2720 Thermal Cycler, Singapore), consisting of denaturation for 30 s at 94 °C, annealing for 90 s at 55 °C, and extension for 1 min at 72 °C, and a final 10 min elongation at 72 °C. Expected amplicons for
mcr-1 (309 bp),
mcr-2 (715 bp),
mcr-3 (929 bp),
mcr-4 (1116 bp), and
mcr-5 (1644 bp) underwent electrophoresis through 1.2% agarose gel followed by staining with ethidium bromide and were visualized under UV light (
Supplementary Figure S1E). Lastly, the obtained results were validated by separate single plex PCR analyses of the
mcr genes.
4.5. Statistical Analysis
Using IBM SPSS statistics data editor (version 21) and GraphPad prism software (version 9.5), verified data were entered and then examined. The bivariate analysis did not include missing data. The mcr gene variations carried by diarrheal pathogens and their phenotypic traits were described using both descriptive and inferential statistical methods. Any associations between categorical data were examined using Pearson’s chi-square test, with the appropriate use of Yate’s continuity correction. Fisher’s Exact test results of a 2 x 2 contingency table were presented in place of the chi-square test results if the predicted frequency of the test cannot be assumed. Two-tailed p-values were computed, with a significance level of 0.05.
4.6. Ethics Statements
This study was authorized [No. UAMC/ERC/Recommend-62/2018, dated 09.07.2018] by the Ethics Review Committee of Uttara Adhunik Medical College. All research protocols complied with the Declaration of Helsinki regarding the use of human beings in research.
Each adult study participant provided written informed consent before the collection of their urine samples. For patients under the age of 18, parents or legal guardians were separately asked for written informed permission. Patient’s identities were anonymized.
Author Contributions
S.S and R.M.N.: Sample collection, Methodology, Investigation, Formal analysis, and Manuscript drafting; M.B.H: Data acquisition, Investigation, Data validation, Visualization; U.L.U and M.A.A.: Methodology, Data curation, Writing-Reviewing, Validation and Editing; A.S.M.M and M.R.K.K: Clinical and demographic data acquisition, Resources, Conceptualization, Visualization and Validation; SN: Project administration, Resources, Methodology, Validation, and Visualization; B.G.: Conceptualization, Methodology, Re-writing and Editing, Visualization; SI: Conceptualization, Supervision, Resources, Data curation and analysis, Writing-Reviewing, and Editing, Study coordination. All authors have read and agreed to the published version of the manuscript.