High Prevalence of Drug Resistance and Class 1 Integrons in Escherichia coli Isolated From River Yamuna, India: A Serious Public Health Risk

Globally, urban water bodies have emerged as an environmental reservoir of antimicrobial resistance (AMR) genes because resistant bacteria residing here might easily disseminate these traits to other waterborne pathogens. In the present study, we have investigated the AMR phenotypes, prevalent plasmid-mediated AMR genes, and integrons in commensal strains of Escherichia coli, the predominant fecal indicator bacteria isolated from a major urban river of northern India Yamuna. The genetic environment of blaCTX–M–15 was also investigated. Our results indicated that 57.5% of the E. coli strains were resistant to at least two antibiotic classes and 20% strains were multidrug resistant, i.e., resistant to three or more antibiotic classes. The multiple antibiotic resistance index of about one-third of the E. coli strains was quite high (>0.2), reflecting high contamination of river Yamuna with antibiotics. With regard to plasmid-mediated AMR genes, blaTEM–1 was present in 95% of the strains, followed by qnrS1 and armA (17% each), blaCTX–M–15 (15%), strA-strB (12%), and tetA (7%). Contrary to the earlier reports where blaCTX–M–15 was mostly associated with pathogenic phylogroup B2, our study revealed that the CTX-M-15 type extended-spectrum β-lactamases (ESBLs) were present in the commensal phylogroups A and B1, also. The genetic organization of blaCTX–M–15 was similar to that reported for E. coli, isolated from other parts of the world; and ISEcp1 was present upstream of blaCTX–M–15. The integrons of classes 2 and 3 were absent, but class 1 integron gene intI1 was present in 75% of the isolates, denoting its high prevalence in E. coli of river Yamuna. These evidences indicate that due to high prevalence of plasmid-mediated AMR genes and intI1, commensal E. coli can become vehicles for widespread dissemination of AMR in the environment. Thus, regular surveillance and management of urban rivers is necessary to curtail the spread of AMR and associated health risks.


INTRODUCTION
The gastrointestinal tract of humans and animals is regarded as the primary/natural habitat of Escherichia coli. Besides its natural habitat, E. coli is also found in secondary habitats like aquatic and terrestrial reservoirs (Méric et al., 2013). Aquatic environments, especially urban water bodies, harbor a heterogeneous collection of microorganisms originating from fecal, hospital, agricultural, and veterinary sources. Moreover, several studies have suggested that aquatic environments serve as genetic reactors that promote transfer of antimicrobial resistance (AMR) and virulence genes among bacteria (Baquero et al., 2008;Hasegawa et al., 2018).
The phylogrouping methods commonly used for population/clonal studies of E. coli include multilocus sequence typing (MLST) (Aanensen and Spratt, 2005), multilocus enzyme electrophoresis (MLEE), triplex PCR, etc. Triplex PCR is a widely used rapid and simple technique for phylotyping E. coli. In triplex PCR, three genes (chuA, yjaA, and a gene encoding a fragment of a putative lipase esterase) are PCR-amplified (Clermont et al., 2000). Based on the presence or absence of these three genetic elements, a strain can be assigned to belong to any of the four phylogroups-A, B1, B2, and D (Clermont et al., 2000). Several researchers have studied AMR in E. coli isolated from aquatic environments and have reported that these traits are easily transmissible among bacterial species with the help of several mobile genetic elements like integrons, insertion sequences (ISs), plasmids, and transposons (Su et al., 2012;Koczura et al., 2013;Liebana et al., 2013;Pereira et al., 2013;Kaushik et al., 2018). Integrons are regarded as the primary vehicles that disseminate AMR genes among bacterial species (Gillings, 2014) because they can be located on conjugative plasmids, which enhance their horizontal spread.
River Yamuna is a major river of northern India, which is associated with several anthropogenic activities of the population residing in the National Capital Region of India. It gets contaminated with effluents originating from hospital and municipal wastewaters; discharge from livestock, poultry, and agriculture production plants; industries; etc. The high levels of pollutants are expected to provide a positive selection pressure for increasing AMR in bacterial population residing there (Kümmerer, 2009;Tacão et al., 2012). Thus, it is expected to be a crucial reservoir of a diverse E. coli populations and an ideal ecological niche for studying strains with diverse phenotypes, genotypes, and AMR. In an earlier study published from our laboratory, we had described the β-lactam susceptibilities and β-lactamase genes in 61 E. coli strains of all phylogroups (A, B1, B2, and D) isolated from river Yamuna (Bajaj et al., 2015). Though AMR resistance has been investigated for commensal phylogroups of E. coli isolated from veterinary or clinical sources, only a few studies have investigated the AMR phenotypes, genes, and integrons in commensal E. coli isolated from urban rivers. The commensal strains of E. coli in urban rivers can easily disseminate AMR determinants to pathogenic E. coli or other waterborne pathogens via mobile genetic elements. Thus, it is important to study the AMR determinants and integrons in commensal strains residing in water bodies also. Our study is the first report on AMR phenotypes, plasmid-mediated AMR genes, and integrons in the strains of commensal phylogroups (A and B1) of E. coli from river Yamuna, India.

Sample Processing and Isolation of Escherichia coli
Two hundred water samples were collected from different sites along the entire stretch of the river Yamuna, which flows through the National Capital Region of India in sterile screw-capped bottles, transported to the laboratory on ice, and processed within 6 h of the sample collection. A schematic figure showing the details of the sampling sites has been published earlier (Bajaj et al., 2015). Enrichment of the samples for isolation of Escherichia coli was performed using a published method (Ram et al., 2008). Briefly, 100 ml of water sample was filtered through a 0.45 µm membrane filter (Millipore, MA, United States). The membrane filter was cut into four pieces, and each piece was incubated in 50 ml of MacConkey broth at 37 • C, 220 rpm, overnight. The next day, a loopful of the broth culture was streaked on the surface of MacConkey agar plates and incubated at 35 • C for 18-20 h. One hundred sixty-two typical E. coli colonies were selected and maintained as pure cultures on Luria-Bertani (LB) agar slants at 4 • C. Of these, 126 isolates were presumptively identified as E. coli using API 20E strips (bioMérieux, France). API 20E is a standardized kit of biochemical tests used to identify members of the family Enterobacteriaceae and other non-fastidious Gramnegative rods.
Isolation of Genomic DNA, and PCR Amplification of Gene Encoding 16s rRNA and Phylogrouping Based on Triplex PCR DNA was extracted from the E. coli strains using the boiling lysis method (Rodríguez-Baño et al., 2004). The gene encoding 16S rRNA was PCR-amplified using universal eubacterial forward primer 27F (5 AGAGTTTGATCCTGGCTCAG3 ) and reverse primer 1492R (5 ACGGCTACCTTGTTACGACTT3 ). The contents of the PCR mixture were 1 × PCR buffer (1.5 MgCl 2 , 1.5 mM of KCl, 10 mM of Tris-HCl, and 0.1% Triton X-100), 200 µM of the four dNTPs, 1 U of Taq DNA polymerase (New England Biolabs, Ipswich, MA, United States), 10 pmol of forward and reverse primers, and 1 ng of genomic DNA in a final volume of 25 µl. The PCR conditions and methods for purification of PCR amplicons and sequencing have been described earlier (Singhal et al., 2019). Briefly, the PCR amplicons were purified by HiYield TM extraction kit (RBC Bioscience, New Taipei City, Taiwan) and sequenced at a commercial facility using Sanger sequencing (Invitrogen BioServices India Pvt. Ltd., Bangalore, India). The nucleotide sequence homology was analyzed using the nucleotide BLAST (BLASTn) algorithm available at the National Center for Biotechnology Information (NCBI).
The phylogenetic profiles of all the isolates were determined by triplex PCR (Clermont et al., 2000); and 40 strains representing the commensal phylogroups (A and B1) were selected for studying the AMR phenotypes, plasmid-mediated AMR genes, and integrons.

Detection and Analysis of Antimicrobial Resistance Genes
Detection of Genes Encoding β-Lactamases and Genetic Environment of bla CTX−M−15 Genes encoding β-lactamases and ESBLs, viz., bla TEM , bla SHV , and bla CTX−M , were detected by PCR amplification of the AMR genes using group-specific primers that amplified the internal coding regions of the genes (Dhanji et al., 2011;Bajaj et al., 2015). The presence of plasmid-encoded AmpC enzymes of the CMY types was determined by PCR amplification of the AMR gene using published primers (Pérez-Pérez and Hanson, 2002).
The promoter region and genetic environment of bla CTX−M−15 were studied in bla CTX−M−15 -positive E. coli strains by PCR amplification of the corresponding regions using the primers and methods described earlier (Saladin et al., 2002;Dhanji et al., 2011). The primers and the annealing temperatures for amplification of the bla TEM , bla SHV , and bla CTX−M and genetic environment of bla CTX−M are described in Table 1. The contents of the PCR mixture and methods for purification of the PCR amplicons and sequencing were the same as used for 16s rRNA gene sequencing. The nucleotide sequence homology was analyzed using the nucleotide BLAST (BLASTn) available at NCBI.
The presence of aminoglycoside resistance genes was determined by PCR amplification of the genes encoding for linked strA-strB genes and four types of plasmidmediated 16S rRNA methylases-armA, rmtA, rmtB, and rmtC-using published primers and annealing temperatures described in Table 1 (Sunde and Norstrom, 2005;Yamane et al., 2007).
The presence of tetracycline resistance genes was determined by PCR amplification of the tetracycline efflux gene tetA using self-designed primers and genes encoding ribosome protective proteins tetM and tetW using published primers and annealing temperatures (Aminov et al., 2002; Table 1). The contents of the PCR mixture and protocols for purification of the PCR amplicons, sequencing of PMQR, and aminoglycoside-and tetracycline-resistance genes and homology search were the same as used for the 16S rRNA genes.

Detection and Analysis of Integrons and Gene Cassettes
The presence and distribution of integrase genes intI1, intI2, intI3, and integron class 1 gene cassette were determined by PCR amplification using published primers (Kraft et al., 1986;Goldstein et al., 2001;White et al., 2001). The variable regions (VRs) of class 1 integrons, which mainly contain an array of gene cassettes, are flanked at 3 by a conserved segment containing qacE 1 and genes coding for quaternary ammonium and sulfonamide, respectively. VRs were investigated as previously reported in all isolates containing class 1 integrons (Guo et al., 2011). The PCR conditions for amplifying the VRs of integrons were the same as for amplifying the 16S rRNA gene except for the primers and annealing temperatures, which have been summarized in Table 1. The PCR amplicons were purified and sequenced as  qacE1-Fsul1-R 5 -AAGTAATCGCAACATCCG-3 5 -GGGTTTCCGAGAAGGTGATTGC-3 qacE 1, sul1 878 57 Bass et al., 1999;Nandi et al., 2004 described for the 16S rRNA gene, and nucleotide sequence homology was analyzed using the nucleotide BLAST (BLASTn) available at NCBI.

Accession Numbers
Gene sequencing revealed that the bla CTX−M−15 genes of the seven bla CTX−M−15 -positive strains were identical to each other; hence, the partial coding sequence (CDS) of a representative strain (KP20) was submitted to GenBank (NCBI) with the accession number KF040057. The partial CDS of bla TEM−1 was also identical to each other; hence, the sequence of one representative strain (IP1N) was submitted to GenBank (NCBI) with the accession number KF055435. The partial CDS of qnrS was also identical; hence, the CDS of a representative strain (KP20) was submitted to NCBI GenBank under accession number KF055436. The partial CDS of tetA gene of all the three strains (IS47, WB3, and KKC) was submitted under the accession numbers KJ409940-KJ409942.

Molecular Identification and Phylogrouping Based on Triplex PCR
The results of 16S rRNA gene sequencing and homology search using BLAST confirmed that the strains presumptively identified using API 20E strips ( Earlier studies have reported that E. coli strains of all phylogroups were present in river Yamuna (Bajaj et al., 2015;Kaushik et al., 2018). Phylogroups A and B1 of E. coli represent commensal strains, while phylogroups B2 and D represent pathogenic strains (Herzer et al., 1990;Bingen et al., 1998;Lecointre et al., 1998;Picard et al., 1999). Several studies have indicated that the prevalence of virulence genes in commensal strains of E. coli was lesser than in pathogenic strains (Johnson, 1991;Boyd and Hartl, 1998;Lecointre et al., 1998;Picard et al., 1999) but that the commensal strains can easily disseminate AMR determinants to pathogenic E. coli or other waterborne pathogens via mobile genetic elements. Thus, AMR determinants and integrons were investigated in these 40 commensal E. coli strains.

Phenotypic Testing of Antimicrobial Susceptibilities and Extended-Spectrum β-Lactamase Production
Antibiotic susceptibility testing revealed that 95% (n = 38) of the commensal E. coli strains were resistant to ampicillin, while 32% (n = 13) of the strains were resistant to piperacillin. Among the cephalosporins, 42.5% (n = 17) strains were resistant to cefazolin (first-generation cephalosporin), 17.5% (n = 7) to cefuroxime (second-generation cephalosporin), 22.5% (n = 9) to cefotaxime (third-generation cephalosporin), and 15% (n = 6) to cefepime (fourth-generation cephalosporin). The fact that commensal waterborne E. coli were less resistant to newgeneration cephalosporins than ampicillin is normal because ampicillin was a widely prescribed broad-spectrum penicillin, and over time, bacteria might have developed resistance to this antibiotic. Earlier studies have also reported that ampicillin resistance was highly prevalent in commensal strains of E. coli isolated from India and from other parts of the globe like Vietnam, China, Sudan, and Thailand (Dyar et al., 2012;Abdelgader et al., 2018;Lugsomya et al., 2018;Singh A. K. et al., 2018;Purohit et al., 2019). With regard to ESBL production, 17.5% (n = 7) of the commensal E. coli strains tested positive, while none of the strain tested positive for AmpC production ( Table 2). All the ESBL-producing strains were resistant to four or more β-lactam antibiotics, and 57% (n = 4) of the ESBL producers were resistant to ciprofloxacin. Since, penicillins and cephalosporins are the most frequently used antibiotics in India, it is normal that all the seven ESBL-producing E. coli strains were resistant to many antibiotics of these classes. The fact that some ESBL producers were also resistant to ciprofloxacin and streptomycin/kanamycin suggests that besides β-lactam antibiotics, resistance to other antibiotic classes also exhibited coselection.
With regard to quinolone resistance, 35% (n = 14) strains were resistant to older quinolones like nalidixic acid, while 20% (n = 8) and 15% (n = 6) strains were resistant to newer quinolones like ciprofloxacin and ofloxacin, respectively. In this regard, our results are similar to those of other studies that reported lower ciprofloxacin and ofloxacin resistance in waterborne E. coli isolated from other parts of the world (Odonkor and Addo, 2018). Aminoglycoside resistance was observably less prevalent with 12.5% (n = 5) strains resistant to older aminoglycosides like streptomycin and less than 7% strains resistant to new aminoglycosides like kanamycin, tobramycin, netilmicin, and amikacin. An earlier study also reported that a low level of aminoglycoside resistance was present in E. coli strains isolated from aquatic environments of Kuala Lumpur, Malaysia (Hara et al., 2018). With regard to tetracycline, 27.5% (n = 11) of the strains exhibited resistance. The frequent use/misuse of ampicillin, streptomycin, and tetracycline due to frequent prescription, availability, and affordability might be a probable reason for higher bacterial resistance to these antibiotics (Shakya et al., 2013). Additionally, 57.5% (n = 23) E. coli strains were resistant to at least two antibiotic classes, and 20% (n = 8) of the strains were multidrug resistant (MDR), i.e., resistant to three or more antibiotic classes. MDR E. coli were defined as bacteria resistant to antibiotics belonging to three or more antimicrobial classes (Magiorakos et al., 2012). An analysis of MAR index revealed that the MAR indexes of 13 MDR E. coli strains were quite high (>0.2). Tambekar et al. (2006) reported that bacteria isolated from environments where several antibiotics are used usually show MAR index >0.2. The high MAR index of the E. coli strains observed in this study is not surprising because river Yamuna is highly contaminated with effluents originating from hospital and municipal wastewaters; discharge from livestock, poultry, and agriculture production plants; etc.

β-Lactam Resistance and Extended-Spectrum β-Lactamase Encoding Genes
With regard to β-lactam resistance genes, bla TEM −1 was present in 95% of the strains (n = 38) followed by bla CTX−M −15 , which was present in 15% (n = 7) of the E. coli strains. The plasmid-encoded AmpC enzymes (CMY types) were not found in any strain. The presence of bla TEM−1 correlated well with ampicillin resistance in all the strains. Similarly, all the bla CTX−M−15 -positive strains showed phenotypic production of ESBLs. Earlier studies had reported that bla TEM −1 was widely present in E. coli strains isolated from water bodies of India and other countries like Spain, Australia, France, China, and Poland (Lartigue et al., 2002;Tristram and Nichols, 2006;Garcia-Cobos et al., 2008;Ortega et al., 2012;Liu et al., 2014;Ojdana et al., 2014;Bajaj et al., 2015;. CTX-M enzymes belong to the family of ESBLs and are the most widely disseminated ESBLs among Enterobacteriaceae all over the globe . Among these, bla CTX−M−15 is the most widely globally disseminated CTX-M type, which was  first reported from the Indian isolates in 2001 (Karim et al., 2001;Poirel et al., 2002). Later, several studies also reported the prevalence of bla CTX−M−15 in aquatic E. coli isolated from India (Bajaj et al., 2015;Kaushik et al., 2019). Previous studies have associated bla CTX−M−15 in aquatic E. coli with the pathogenic phylogroups B2 (especially those belonging to the genetic lineage ST131) and D (Nicolas-Chanoine et al., 2008;Coque et al., 2008). However, our study revealed that CTX-M-15 type ESBLs were present in the commensal phylogroups A and B1, also.

Aminoglycoside Resistance Genes
The linked strA-strB genes are the most widely prevalent streptomycin resistance genes in E. coli worldwide and encode for phosphotransferases (Poirel et al., 2018). However, strA-strB genes were present in only 12% (n = 5) of the waterborne E. coli. Of these, four strains exhibited phenotypic resistance to streptomycin, while one strain (ISF), despite harboring strA-strB, was phenotypically susceptible for streptomycin. The 16S rRNA methylases methylate certain amino acid residues of the 16S RNA, resulting in resistance to amikacin, tobramycin, gentamicin, and netilmicin (Griffey et al., 1999). Of the four types of plasmid-mediated 16S rRNA methylase investigated, only armA was found to be present in 17.5% (n = 7) of the strains. Earlier studies also reported that a low level of aminoglycoside resistance was present in E. coli strains isolated from aquatic environments in Kuala Lumpur, Malaysia (Hara et al., 2018). No correlation was observed in aminoglycoside resistance and presence of armA, except in the E. coli strain IS45, which also exhibited phenotypic resistance to amikacin, tobramycin, kanamycin, and netilmicin.

Plasmid-Mediated Quinolone Resistance and Tetracycline Resistance Genes
Of the several PMQR genes (qnrA, qnrB, qnrC, qnrD, qnrS, qep, and aac) tested, only qnrS was detected in the E. coli strains isolated from river Yamuna. The qnrS1 was detected as the predominant PMQR gene in about 17% of the aquatic strains (n = 7). Earlier studies have also reported that the qnrS type gene was the most frequently detected PMQR gene in E. coli isolated from environmental E. coli worldwide (Bonemann et al., 2006;Cattoir et al., 2008;Rodriguez-Mozaz et al., 2015;Varela et al., 2016;Hara et al., 2018). The presence of qnrS1 did not correlate with fluoroquinolone resistance, except in E. coli strain NG35. This suggests that the presence of the qnrS gene alone might not be a true indicator of fluoroquinolone resistance and the isolate despite that the presence of qnrS might exhibit phenotypic susceptibility for fluoroquinolones (Mahmud et al., 2020). Of the three tetracycline resistance genes, tetM and tetW were absent and only tetA was present in 7.5% (n = 3) of the strains. Earlier studies have also indicated that tetracycline efflux-related genes like tetA, tetB, and tetC were more prevalent than ribosomal protectionrelated genes (like tetM and tetW) in waterborne E. coli (Zhang et al., 2015;Stange et al., 2016). The presence of tetA correlated well with phenotypic resistance because the three strains that harbored tetA also exhibited tetracycline resistance.

Genetic Environment of bla CTX−M−15
Of the 40 E. coli strains investigated in this study, only seven strains (15%) harbored the bla CTX−M−15 gene. The upstream region of the bla CTX−M−15 was analyzed in these seven strains by PCR amplification and sequencing. Gene sequencing revealed that the ISEcp1 was present in the upstream region of bla CTX−M−15 and orf477 was present in the downstream region of all the seven strains ( Table 2). Several investigators have also reported the presence of orf477 in the downstream region of bla CTX−M−15 (Eckert et al., 2006;Dhanji et al., 2011;Wang et al., 2014;Ben Said et al., 2016). ISEcp1 is the most common and widely reported IS element (Dhanji et al., 2011;Liu et al., 2014;Upadhyay et al., 2015;Ben Said et al., 2016;. In the present study, IS sequence ISEcp1 was found to be present at 48-bp upstream region of bla CTX−M−15 . The 42to 266-bp upstream region has been reported as the preferred insertion site of ISEcp1 for different bla CTX−M genes like CTX-M-1, CTX-M-2, and CTX-M-9. An analysis of the −35 and −10 promoter regions of the bla CTX−M−15 gene of the seven strains revealed that the −35 (TTGAAA) and −10 (TACAAT) regions were present within 3 terminus end of ISEcp1 and 48 bp away from the bla CTX−M−15 start codon. The same organization was previously reported from E. coli strains isolated from different countries of the world (Saladin et al., 2002;Boyd et al., 2004;Canton and Coque, 2006;Lavollay et al., 2006;Dhanji et al., 2011; Frontiers in Microbiology | www.frontiersin.org Liu et al., 2014). The presence of ISEcp1 along with bla CTX−M has been reported from E. coli strains isolated from different parts of the world, indicating that ISEcp1 might be evolutionary associated with bla CTX−M (Karim et al., 2001;Saladin et al., 2002;Dhanji et al., 2011;Liu et al., 2014;Wang et al., 2014). ISEcp1 can mobilize an adjacent gene as a part of transposition units of varying sizes (Zong et al., 2010). It has also been reported that ISEcp1 helps in improving the expression of bla CTX−M in enteric bacteria (Poirel et al., 2003). This is a matter of great concern because the subsequent transfer of ESBL genes from these commensal E. coli to pathogenic E. coli or other bacteria in aquatic water bodies might pose a serious health challenge (Figueira et al., 2011).

Detection and Analysis of Integrons
The class 1 integron gene intI1 was detected in 75% (n = 30) of the isolates ( Table 3). None of the strains harbored class 2 and 3 integrase genes intI2 and intI3. Though some strains of E. coli reportedly harbored class 2 integrons, mostly class 1 integrons have been reported from E. coli isolated from India (Kaushik et al., 2018). Several studies have indicated that integrons of the class 3 were absent in E. coli isolated from water bodies across the globe (Laroche et al., 2009;Su et al., 2012;Pereira et al., 2013). Of the 30 intI1 harboring strains, gene cassette arrays were detected in only 10% (n = 4) of the intI1-positive E. coli strains. Three different types of gene cassette arrays of class 1 were present in the downstream region of intI1 whose size ranged from 1.7 to 2.8 kb ( Table 2). Gene sequencing revealed that five gene cassettes of the dihydrofolate reductase (dfr) resistance gene family (dhfr12 and dfrA1), aminoglycoside (aad) resistance gene family (aadA1, aadA2, and aacA4), and chloramphenicol (CHL) resistance gene family (catB3) were present ( Table 3). The prevalence of class 1 integron intI1 in Indian aquatic isolates was quite high (75%) and alarmingly more than that reported for E. coli isolated from global aquatic environments (Dolejska et al., 2009;Laroche et al., 2009;Pereira et al., 2013;Ghaderpour et al., 2015;Sidhu et al., 2017). E. coli isolates harboring class 1 integrons have been associated with a significantly higher probability for multidrug resistance than those devoid of class 1 integrons (Chen et al., 2011). Moreover, class 1 integron genes intI1 are accompanied by resistance genes for disinfectants and heavy metals (Partridge et al., 2001) and can also easily horizontally transfer between strains originating from different sources (Nagachinta and Chen, 2008;Zhang et al., 2009). Thus, due to high prevalence of class 1 integron gene intI1, the commensal strains of E. coli can become vehicles for widespread dissemination of antibiotic resistance to pathogenic E. coli and other waterborne bacterial pathogens.

CONCLUSION
Our results indicated a high prevalence of drug resistance in Escherichia coli strains of river Yamuna. With regard to plasmid-mediated AMR genes, bla TEM−1 was present in 95% strains followed by qnrS1 and armA (17% each), bla CTX−M−15 (15%), strA-strB (12%), and tetA (7%). Though most of the earlier studies have reported that bla CTX−M−15 in waterborne E. coli was mostly present in pathogenic phylogroup B2, our study revealed that CTX-M-15 type ESBLs were present in the commensal phylogroups A and B1, also. The genetic organization of bla CTX−M−15 was similar to that reported for E. coli globally, and ISEcp1 was present in the upstream region of bla CTX−M−15 . Though integrons of classes 2 and 3 were absent, class 1 integron gene intI was detected in 75% of the isolates, which indicates its high prevalence in the E. coli isolates of the river Yamuna than that reported, globally. The presence of MDR phenotypes, plasmid-mediated AMR genes, and class 1 integron gene intI1 in E. coli is a serious public health risk, because these commensal strains can become potent vehicles for widespread dissemination of AMR determinants to pathogenic E. coli and other waterborne pathogens. Thus, our study suggests an urgent need for regular surveillance and management of natural water bodies to curtail the spread of antibiotic resistance in microorganisms.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.

AUTHOR CONTRIBUTIONS
NSS, NS, MK, and JV analyzed the data. NSS and NS wrote the manuscript. JV conceptualized and supervised. All authors contributed to the article and approved the submitted version.