banner

Blog

Jan 19, 2024

Mobile resistome of microbial communities and antimicrobial residues from drinking water supply systems in Rio de Janeiro, Brazil

Scientific Reports volume 12, Article number: 19050 (2022) Cite this article

1319 Accesses

1 Altmetric

Metrics details

Antibiotic resistance genes (ARGs) are widespread in the environment due to the overuse of antibiotics and other pollutants, posing a threat to human and animal health. In this study, we evaluated antimicrobial residues, bacterial diversity and ARGs in two important watersheds, Guandu and São João, that supply drinking water to Rio de Janeiro city, Brazil. In addition, tap water samples were collected from three different cities in Rio de Janeiro State, including the metropolitan area of Rio de Janeiro city. Clarithromycin, sulfamethoxazole and azithromycin were found in untreated water and drinking water in all samples. A greater abundance of Proteobacteria was observed in Guandu and São João watersheds, with most of the sequences belonging to the Gammaproteobacteria class. A plasmidome-focused metagenomics approach revealed 4881 (Guandu), 3705 (São João) and 3385 (drinking water) ARGs mainly associated with efflux systems. The genes encoding metallo-β-lactamase enzymes (blaAIM, blaGIM, blaIMP, and blaVIM) were detected in the two watersheds and in drinking water samples. Moreover, we demonstrated the presence of the colistin resistance genes mcr-3 and mcr-4 (both watersheds) and mcr-9 (drinking water and Guandu) for the first time in Brazil. Our data emphasize the importance of introducing measures to reduce the disposal of antibiotics and other pollutants capable of promoting the occurrence and spread of the microbial resistome on aquatic environments and predicting possible negative impacts on human health.

The environmental impacts that most affect the quality of aquatic ecosystems and, consequently, public health are strongly associated with inadequately treated or untreated wastewater1,2. Water pollution can occur due to a lack of sanitation and/or waste discharge without treatment by point or diffuse sources3,4.

Several substances have been considered emerging contaminants, including new pesticides, antimicrobials, personal care products, some by-products from water disinfection processes, sweeteners such as sucralose, nanomaterials, and some microorganisms5,6. Recent estimates indicate that antimicrobials are the main classes of drugs capable of causing some of the greatest environmental impacts7.

Antimicrobials have been widely used in human and veterinary medicine. However, about 70 to 80% of ingested doses are excreted unchanged and discharged to water bodies, mainly through wastewater generated from hospitals and pharmaceutical industries8. These drugs are only partially removed by wastewater treatment and, depending on the compound, they can still be found at levels ranging between 10 and 1000 ng L−1 in effluents9,10.

The antimicrobials carried by effluent disposal in the environment, even at low levels, are a key signal that promotes gene dissemination and consequently increased resistance11. Many antimicrobials are naturally biodegradable compounds, but synthetic drugs such as quinolones are more resistant to biodegradation in the environment. This leads to prolonged effects on bacterial communities and to a substantial impact on increased resistance. Even when antimicrobial contamination is eliminated, the resistance determinants can be maintained and disseminated within and between microbial populations12,13.

Additionally, the disposal of antimicrobial residues in aquatic environments can not only cause impacts on the biodiversity and function of ecosystems but may also select for antibiotic resistant bacteria (ARB) and stimulate the dissemination of antimicrobial resistance genes (ARG)14. Mobile genetic elements (MGE), including phages, plasmids, and transposons, among others, mediate this spread15. Plasmids, in particular, are rapidly disseminated in the environment and play a major role in microbial evolution and adaptation as vehicles of gene transfer16.

The plasmidome is defined as the total plasmid populations within a given community17. Plasmidome analysis provides information on the composition and structure of the mobile resistome. Therefore, it is considered a promising approach that provides information about the types of plasmids present in the studied microbial community, and the MGE contained in these plasmids18.

The role of non-clinical environments in increasing the spread of ARBs has not been fully elucidated. In general, ARGs are not easily removed from polluted areas, even when the selective pressure exerted by pollutants is gone. This may also explain why ARGs are often found in antimicrobial-free environments19,20. Resistance to antimicrobials, initially confined to hospitals, has also been observed in the natural environment, which raises great concern regarding the impacts on human health. ARBs and ARGs can be dispersed into raw sources of drinking water, mainly through discharge of human and animal waste, wastewater treatment plant, hospital sewage, and agricultural practices such as manure application21,22.

The present study evaluated the presence, distribution, and abundance of antimicrobial residues and ARGs in two important watersheds for drinking water supply to South Central regions of Rio de Janeiro state, Brazil, including the metropolitan region of Rio de Janeiro city, Brazil, using high-performance liquid chromatography coupled to tandem mass spectrometry and a culture-independent approach, respectively. Our study may provide relevant information about the structure, complexity, and content of the plasmidome of these waters, which can pose serious threats to human health.

The samples were grouped as São João watershed, Guandu watershed, and Drinking water. Clarithromycin was the most frequently detected antimicrobial, being observed in 80% (8/10) of the Guandu samples, 40% (4/10) of the São João samples and 36% (4/11) of drinking water samples. In the first collection at São João River mouth, cefoperazone concentrations > 500 ng L−1 were found. The sulfamethoxazole concentration, belonging to the sulfonamide class, ranged from 47.4 to 340.5 ng L−1 in the second collection in the Macacos and Queimados rivers, respectively. The Unamar drinking water sample presented levels of 12.5 ng L−1 of this antimicrobial. Azithromycin was found in the drinking water samples from Itaguaí, in Macacos river and the São João river mouth at a concentration below 10 ng L−1, and 49.9 ng L−1 in the second collection in the Queimados river. Troleandomycin and roxithromycin were detected only in the Leblon drinking water sample at concentrations < 10 ng L−1 (Fig. 1).

Detection frequency and concentration levels of antimicrobials in grouped samples of São João watershed, Guandu watershed, and Drinking water.

Nine bacterial phyla were observed, with Proteobacteria predominating in 93.5% (29/31) of the samples, followed by Actinobacteria and Bacteroidetes. Within the phylum Proteobacteria, the predominant class was Gammaproteobacteria (36 to 46%), followed by Alphaproteobacteria (11–13%). In the Bacteroidetes phylum, the Bacteroidia class was the most abundant (12–17%). In the Actinobacteria phylum, the class Actinobacteria was the most abundant (12–14%) (Fig. 2). There were no significant differences (p > 0.05) in the alpha diversity of microbial communities associated with the seasonality of sample collection. Therefore, the samples were grouped as São João watershed, Guandu watershed, and Drinking water (Fig. 3a).

Relative abundance of bacterial composition at the class level in grouped samples of São João watershed, Guandu watershed, and Drinking water.

Diversity indices (Averages). (A) Comparison of the alpha diversity of the Operational taxonomic unit (OTUs) between the collection of untreated and treated water samples. Shannon coverage analysis; (B) Comparison of alpha diversity of bacterial OTUs among the analyzed water samples. Chao1 coverage analysis; (C) Principal coordinate analysis (PCoA) between bacterial communities present in the Guandu watershed (green), São João watershed (blue), and drinking water (red).

The effect of sample types on alpha bacterial diversity was assessed based on the richness of OTUs (absolute number of taxa), diversity, and uniformity. The bacterial communities in the drinking water samples presented lower diversity, compared to untreated water samples from the watersheds (p < 0.02 with Shannon Index) (Fig. 3A and B). There were no statistically significant differences between these samples so they were pooled. The water samples also showed significant variations (p < 0.001) concerning the beta diversity of bacterial communities. Meanwhile, the samples did not show significant variations (p < 0.103) by the Jaccard method (Fig. 3C). We believe the wide variations observed between the samples reflect the environmental conditions that are very dynamic in these aquatic habitats, such as rainfall, and sewage discharge, for example. As our main goal was to determine antimicrobial concentrations and the presence of ARGs in the plasmidome in these samples, we believe there is no need to collect more samples. More detailed Alpha diversity measures are presented in supplementary material Table S1, revealing that bacterial community richness (Chao1), diversity (Shannon and Simpson), and evenness (Shannon even) varied widely among the samples.

A total of 3,490,453 paired-end reads were generated for São João watershed, 2,719,506 for Guandu and 3,302,359 drinking water samples. After quality control and assembly, 6197 contigs from the São João watershed (mean sequence length 2.043 bp), 4866 from Guandu watershed (mean sequence length 4776 bp), and 5185 contigs from drinking water (mean sequence length 3255 bp) were analyzed. Eighteen distinct subsystems, containing genes attributed to resistance and adaptation to antimicrobials, metals, and other environmental pollutants, were distributed among all samples, according to analyzes at the MG-RAST database. While the Guandu watershed showed the greatest diversity of subsystems (n = 16), São João exhibited the least diversity (n = 4) while samples of drinking water revealed 12 subsystems.

Fifty-seven percent of São João watershed sequences and 4% of Guandu sequences were attributed to Cadmium_resistance; this system was not found in drinking water samples. The Cobalt-zinc-cadmium_resistance subsystem, including zinc, cobalt, and cadmium efflux systems encoded mainly by the cusA gene and the czc operon, was found in 14% of the São João watershed, 16% of Guandu watershed, and 13% of drinking water sequences. The presence of genetic determinants of the Multidrug_Resistance_Efflux_Pumps subsystem was also revealed in the three sample groups (14% São João watershed, 17% Guandu watershed, 10% drinking water) mainly composed of members of the Multidrug And Toxic compound Extrusion (MATE) family in the São João and Guandu watersheds samples (> 90%). In the drinking water samples, the MATE family (37.5%), the superfamily pumps of efflux resistance-nodulation-cell division (RND) encoded by the cmeA gene (12.5%), and the macA/macB macrolide efflux system (25%) were found (Fig. 4a).

(A) Relative abundance of sequences related to antimicrobial resistance subsystems and toxic compounds in grouped samples of São João watershed, Guandu watershed, and Drinking water; (B) Relative abundance of ARGs concerning antimicrobial classes, in the three groups of samples.

Through the CARD database searches, 4881 antimicrobial resistance genes were annotated from the Guandu watershed, 3705 from the São João watershed, and 3,385 from the drinking water samples. This analysis revealed the prevalence of three types of resistance mechanisms: antibiotic efflux, antibiotic target alteration/protection, and antibiotic inactivation, evenly distributed among the samples (Fig. 4b).

In general, ARGs found in the samples encode many types of proteins and enzymes capable of conferring resistance to antimicrobials. The genes capable of conferring resistance to 4 or more classes of antimicrobials were grouped as Multidrug, while the less abundant antimicrobial classes (< 1%) were grouped as Others. Resistance genes to macrolides were prevalent in all three sampling sites (15.2% São João, 14.9% Guandu, and 15.3% drinking water), followed by ARGs against glycopeptides, tetracyclines, fluoroquinolone, β-lactams, and others. The most abundant genes in both watersheds and drinking water samples were macB and tetA(58), both associated with efflux systems. It is worth noting the presence of the genes blaAIM (Guandu and São João watersheds), blaGIM (drinking water, Guandu and São João watersheds), blaIMP (Guandu watershed), and blaVIM (drinking water, Guandu, and São João watersheds) encoding metallo-β-lactamase enzymes. Even more surprising, the presence of the mcr-3 (Drinking water, Guandu and São João watersheds), mcr-4 (Guandu watershed), and mcr-9 (Drinking water and Guandu watershed) genes that can provide resistance to colistin and are currently considered a serious public health issue, was also observed.

The consumption of contaminated drinking water is a major pathway for environmental ARB to enter into the human gut22. The São João and Guandu watersheds are frequently impacted by elevated levels of chemical and fecal pollution, due to a lack of sewage treatment in the surrounding cities. Therefore, this study aimed to determine the presence of antimicrobial residues in water resources and to evaluate their possible impacts on the microbial resistome. Antimicrobial residues in the aquatic environment may originate from hospital, domestic, rural effluents (aquaculture, livestock)23,24, and the pharmaceutical industry25.

In this study, substances from the β-lactam, macrolide, and sulfonamide antimicrobial classes were detected in samples of untreated waters, being the Guandu watershed the one with the highest rates of antimicrobials. Clarithromycin, azithromycin, and sulfamethoxazole have also been found in drinking water samples. It is worthy to noting that clarithromycin residues were observed in all three environments, as well as a high proportion of resistance genes to this antibiotic, suggesting an impact of the drug on the spread of resistance genes26. In addition, troleandomycin and roxithromycin were detected only in drinking water samples. Macrolides, such as clarithromycin and azithromycin, are widely administered in human and animal medicine and can be carried to water resources through agricultural soil and the application of sewage sludges or fertilizers27.

Many antimicrobials have already been detected in drinking water in several developed countries at levels generally < 100 ng L−128. High levels of carbamazepine, clofibric acid, and sulfamethoxazole have also been observed in drinking water in countries in Europe and North America29. However, our data showed that although cephalexin was found in 60% (6/10) of the samples from Guandu and 20% (2/10) from São João in levels ranging between < 10 ng L−1 and > 500 ng L−1 it was not detected in drinking waters. Some antimicrobials can be eliminated via abiotic or biotic degradation, but their continued introduction can make them pseudo persistent in aquatic environments29. The presence of antimicrobials in drinking water is due to their incomplete removal during conventional treatment steps in WWTPs. In addition, antimicrobial residues can accelerate the emergence and evolution of ARB and ARGs in the environment30.

Water also represents an important way of spreading bacteria between different aquatic environments, including freshwater habitats that harbor the richest bacterial diversity31. It is well known that heterotrophic prokaryotes play important roles in the structure and dynamics of trophic networks and the remineralization of organic matter32. We find that most of the environments analyzed (93.5%; 29/31) are dominated by the phyla Proteobacteria, Actinobacteria, and Bacteroidetes. Proteobacteria have great metabolic diversity, a fact that allows their dissemination in the most varied environments33.

A relevant finding of our study was the high abundance of Actinobacteria in the Juturnaíba Dam, which is one of the groups best known for containing organisms that produce antimicrobials and carriers of MDR profiles, and one of the most prevalent sources of ARGs34. Actinobacteria have many acetyltransferases and phosphotransferases, which represent the greatest resistance mechanism to aminoglycosides34.

Physical, chemical, and biological pollution can influence the microbial composition, with potential effects on water quality and safety. Although maintaining a safe and reliable supply of drinking water is of critical importance, few potentially pathogenic microorganisms are recognized and even less are regulated35. The homogeneity between the bacterial communities showed in our study can be explained by the urbanization process which leads to the discharge of high loads of untreated waste and sewage containing feces and xenobiotic compounds into bodies of water36,37.Environments that are not impacted have a higher prevalence of phyla Acidobacteria and Verrucomicrobia36. Also, the treatment of drinking water aims to reduce the microbial load, which explains why alpha and beta diversities of the microbial communities of untreated watershed waters were higher compared to the drinking water communities38,39.

Our plasmidome analyses revealed a high abundance of the czc (cobalt-zinc-cadmium) efflux system in all environments analyzed. This system is strongly associated with locations impacted by oil, sludge, metals, and other urban waste40. There is great concern about the relationship between this system and the increase in antimicrobial resistance. P. aeruginosa strains isolated from urinary catheters, susceptible to carbapenems and carrying the czc operon, have demonstrated resistance to imipenem, when exposed to zinc. Also, analysis of cross-resistance mechanisms revealed co-regulation of czcR overexpression and a decrease in oprD expression, that encodes a channel associated with resistance to carbapenems, especially imipenem41,42.

It is already known that ARBs can survive the selective pressures that occur during the water treatment process43. Meanwhile, the removal of ARGs varies depending on the water treatment scheme. Chlorine disinfection can eliminate many ARBs but does not destroy ARGs resulting in their discharge into aquatic environments44,45.

Several ARGs detected in our study have been previously reported in environmental waters, sediments and soils, drinking water, and WWTP46,47. The blaNDM and blaCTX-M-type genes often found on mobile genetic elements such as plasmids, have already been reported in drinking water globally22,48. We also show the presence of several genes encoding MBLs (metallo-β-lactamases), carbapenemases, and the mcr (mobile colistin resistance) genes, which confers resistance to colistin, associated with plasmids.

It is noteworthy that, as far as we know, this is the first study reporting the presence of the mcr-3 and mcr-9 genes in drinking water samples. Some mcr-like genes have already been described in water systems, such as mcr-1 which, although it was first described in Enterobacteriaceae isolated from animals, food, and humans in China49, has already been revealed in China's water systems46. So far, the presence of the mcr-9 gene has not been described in Brazil, either in bacterial isolates or through metagenomic studies. Although colistin residues were not revealed in our study, the occurrence of mcr-like genes in the Guandu and São João watersheds and drinking water samples could be related to the low levels of colistin and/or other drugs important for spreading of mcr resistance genes in those environments. Stanton et al.50, demonstrated evidence that antibiotics in low concentrations (below minimum inhibitory concentrations) promote the emergence and persistence of antibiotic resistance in natural environments. In fact, colistin has been heavily added to animal feed, as a growth promoter in cattle, pigs, and poultry, in Brazil51. The same was observed in some Asian countries, including China, India, Japan, and Vietnam where colistin is widely used to improve weight gain in animals52. In Europe, it is mainly used to treat infections caused by Enterobacteriaceae in pigs, chickens, cows, sheep, and goats53.

The β-lactamase enzyme production is the most common mechanism of bacterial resistance to β-lactam antimicrobials53. In this study, the genes encoding carbapenemases were found only in one of the evaluated watersheds (Guandu). However, Gram-negative bacteria carrying genes for resistance to carbapenems have already been isolated from samples of rivers, wastewater, and drinking water, emphasizing their high potential for dissemination in the environment54,55.

In addition to other carbapenem resistance genes, the blaGIM and blaVIM genes were found in the drinking water samples in the present study. In general, strains carrying MBL genes are MDR which are a serious therapeutic problem in clinical isolates. The description of genes encoding MBLs associated with MGEs has considerably increased attention to these enzymes, including them among the main threats to human health for the twenty-first century56,57.

Widespread antimicrobial resistance represents a serious threat to human health since it is associated with the loss of the therapeutic potential of antibiotics and the consequent morbidity and mortality58. Currently, studies suggest that chemical compounds that are not antimicrobial can also select and stimulate antimicrobial resistance, such as heavy metals59, disinfectants60, disinfection by-products61, and nanomaterials62.

Treatment of drinking water in treatment plants is not aimed at removing disinfectants, and quite often, a residual amount is maintained in the supply system to prolong water quality during distribution63. However, the consequences and selective pressures of such residuals are not usually considered regarding the presence of ARGs, ARBs, and MGEs. Besides that, these plants do not efficiently remove antimicrobials and metals64,65. Hence, this selection pressure caused by disinfection, antimicrobials and metallic agents could continue throughout distribution, and bacteria carrying resistance determinants, or capable of acquiring them, may persist in the drinking water11,66.

Current regulations do not establish the monitoring and control of ARBs, ARGs, and MGEs in drinking water and wastewater. Our data emphasize the importance of introducing measures to reduce the disposal of antibiotics and other pollutants capable of promoting the enrichment and maintenance of the microbial resistome. Besides that, our data indicate that mitigation strategies should be put in place to reduce the risk of AMR and to prolong the efficacy of the currently available antimicrobial agents for use in animals and humans.

Five collection sites were selected in Guandu Watershed (22°50′22.11″ S and 43°36′36.70″ O) (Queimados, Guandu, Piraí and Macacos rivers and Guandu Dam), and five sites in São João Watershed (22°37′36.60″ S and 42°17′54.36″ O) (Capivari, Bacaxá and São João rivers, São João river mouth and Juturnaíba Dam). Samples (5 L from each site) were collected in sterile bottles, six months apart (January and June/2015). In addition, in January/2015, 11 household tap drinking water samples (5 L from each point) were collected in different neighborhoods (one sample per neighborhood) from Rio de Janeiro city (Centro, Copacabana, Ilha do Governador, Jacarepaguá, Jardim Botânico, Leblon, Realengo, Santa Teresa, Vista Alegre) and from the small towns of Unamar and Itaguaí (one sample per city) in Rio de Janeiro State. All samples were collected in three replicates and refrigerated until processed in the laboratory within 24 h. All experiments were performed using kits and internal controls. Metadata for all samples is included in supplementary Table 2.

All drinking water is provided by Guandu Water Treatment Plant (GWTP). The GWTP is in the Guinness Book as the world's largest drinking water treatment plant in continuous production with a flow rate of about 45,000 L per second67. Upon reaching the GWTP, a chemical coagulant is added to the water, followed by a polyelectrolyte. With the coagulant adequately dispersed, the water passes through hydraulic flocculators, whose controlled agitation promotes the collision of the particles and consequently the agglutination, forming the flocs. The water then enters the sedimentation tanks (decanters), where the speed is reduced, and the flakes already formed and with greater weight sink to the bottom. The clarified water is collected through channels on the blade's surface and distributed to the filtration system. The filters are composed of layers of sand with a granulometry capable of retaining the finest particles that are still present in the clarified water. After being filtered, the water flows to the contact reservoirs, where disinfection occurs with the addition of chlorine. After being disinfected, the water is fed through underground channels to the high-pressure lifts. In these channels, the pH correction occurs with the addition of quicklime. Fluoride is also applied to treated water as an auxiliary agent in the fight against dental caries68.

Amoxicillin tryhidrate (AMOX), ampicillin (AMPI), cefaclor (CFCL), cefadroxil (CFDX), cefalexin hydrate (CFLX), cefazolin (CFZL), clarithromycin (CLA), ciprofloxacin hydrochloride (CPF), norfloxacin (NOR), tetracycline hydrochloride (TC) and sulfamethoxazole (SMZ) were chemical reference substances from the Brazilian Pharmacopeial Convention (Santa Maria, RS, Brazil). Azithromycin dehydrate (AZI), roxithromycin (ROX), spiramycin (SPI), oleandomycin (OLE), tilmicosin (TILM), and cefquinome sulphate salt (CFQN) were obtained from Dr. Ehrenstorfer (Augsburg, Germany). Oxytetracycline (OTC), doxycycline hyclate (DC), hydrochloride salts of chlortetracycline (CTC) and demeclocycline (DMC), dapsone (DAP), sulfacetamide (SCT), sulfadimethoxine (SDM), sulfamerazine (SFM), sulfamethazine (SMT), sulphaquinoxaline (SQN), sulfathiazole (STZ), tylosin tartarate (TYL), troleandomycin (TRO), erythromycin (ERY), cephapirin sodium salt (CPPN), ceftiofur (CFTF), cefoperazone (CFPZ), benzylpenicillin sodium salt (PENG), oxacillin sodium salt hydrate (OXA), moxifloxacin (MXF) and ofloxacin (OFX) were supplied from US Pharmacopeial Convention (Rockville, MD, USA). Phenoxymethylpenicillin potassium salt (PENV), cloxacillin sodium salt hydrate (CLOX), dicloxacillin sodium salt hydrate (DCLOX) and nafcillin sodium salt (NAFC) were supplied from WHO Collaborating Centre for Chemical Reference Substances (Stockholm, Sweden). Methacycline (MTC), 4-epioxytetracycline (4-EOTC), 4-epitetracycline (4-ETC) and 4-epichlortetracycline hydrochloride (4-ECTC) were acquired from Acros (Pittsburgh, PA, USA). Ampicillin-d5 (AMPID5) was purchased from Purity Grade Standards (San Francisco, CA, USA). Desacetylcephapirin (DESAC) was supplied from Bristol-Myers Squibb (New York, USA).

Methanol (MeOH) and acetonitrile (ACN) HPLC grade, hydrochloric acid (HCl) and formic acid (FOA) analytical grade were purchased from Merck (Darmstadt, Germany). Sodium hydroxide (NaOH), acetone (ACE) and ascorbic acid (ASA) were purchased from Merck (Darmstadt, Germany). Ethylenediaminetetracetic acid disodium dihydrate (EDTA) was acquired from Calbiochem (Gibbstown, NJ, USA). Ultrapure water was obtained from a Milli-Q purification system (Millipore, Bedford, MA, USA).

Solid-phase extraction (SPE) was performed with 60 mg Oasis® HLB cartridges from Waters Corp. (Milford, MA, USA). Membrane filters of polyvinylidene fluoride (PVDF) with pore size 0.22 µm were purchased from Millipore (Billerica, MA, USA).

The stock standard solutions were prepared to obtain a concentration of approximately 1000 μg mL−1. Stock solutions of β-lactams (BL) were prepared in water while those of fluoroquinolones (FQ) in a 0.03 mol L−1 NaOH. Finally, macrolides (MC), sulfonamides (SF) and tetracyclines (TC) solutions were prepared in MeOH. The amount weighed for each standard was calculated considering purity, water content and free acid/basic corrections. The solutions were transferred to microtubes and stored in a freezer at − 70 °C or below. DMC and AMPID5 were used as internal standards.

Intermediate and working standard solutions were freshly prepared at several concentrations by appropriate dilution of stock standard solutions.

The extraction methodology for antimicrobial residues was based on the standard method from the United States Environmental Protection Agency (US EPA)—Method 169469, with modifications described by Monteiro et al.70.

Samples were previously filtered through filter paper and 0.22 µm PVDF membrane filter. A 50 mL aliquot of each sample was spiked with 100 ng L−1 of the internal standards, acidified to pH 2.5 with HCl, and 2 mL of 25 mg L−1 EDTA stock solution was added. For drinking water samples, 2 mL of 625 mg L−1 ASA was added to reduce any residual chlorine. This solution was applied to an Oasis® HLB cartridge that had been previously conditioned in sequence with 3 mL of MeOH, 3 mL of ultrapure water and 3 mL of ultrapure water acidified to pH 2.5 with HCl. After being washed twice with 2 mL of water, SPE cartridges were vacuum-dried (− 35 kPa) for 2 min. Antimicrobials were eluted with three portions of 2 mL MeOH and one portion of 2 mL ACE, using gravity flow only. 4 mL aliquots of the eluate were transferred to two centrifuge tubes and evaporated to dryness with N2 in a temperature up to 47 °C. The residues were reconstituted with 1 mL of 0.1% FOA:MeOH (80:20, v/v) for TC and SF analysis and 1 mL of MeOH:H2O (65:35, v/v) for BL, MC and FQ analysis, vortexed for 30 s and filtered through a 0.22 µm polyvinylidene fluoride (PVDF) syringe filter into amber auto-sampler vials.

The chromatographic analysis was performed on a Shimadzu Prominence HPLC (Kyoto, Japan) equipped with a quaternary pump (LC-20AD), a membrane degasser (DGU-20A5), an auto-sampler (SIL-20AC), a column oven (CTO-20AC) and a system controller (CBM-20A) interfaced to a triple quadrupole mass spectrometer (API5000, Applied Biosystems/MDS Sciex, Foster City, CA, USA) with the TurboIonSpray® source. Analyst® V1.4.2 LC/MS control software was used. The analytical column was a Pursuit™ C18 RS (100 mm × 2 mm i.d., 3 µm particle size, 200 Å), with a respective guard column (Varian, Lake Forest, CA, USA). Mobile phases A, B and C were prepared using water, ACN and MeOH, respectively, all with 0.1% FOA. A gradient elution program for TC and SF method was used with a flow rate of 0.15 mL min−1 at 25 °C and for BL, MC and FQ another gradient elution was used with a flow rate of 0.30 mL min−1 at 35 °C. The injection volume was 25 µL for both methods. The auto-sampler was set at 4 °C. Positive electrospray ionization technique (ESI +) in Multiple Reaction Monitoring (MRM) acquisition mode was used to monitor two ions for each substance.

A six-point calibration set was freshly prepared by spiking varying levels of working standard solutions in ultrapure water. The analytical curves for all analytes in the concentration range from 25 to 1000 ng L−1 were constructed in order to quantify the analytes in samples.

The chromatographic peaks were integrated with the IntelliQuan algorithm of the Analyst® software. A signal-to-noise ratio of the peaks equal or greater than 3:1 for at least 2 transitions was required for detection. Relative retention times and relative abundances between quantification and confirmation MRM transitions in both matrix-fortified standards and samples were used as confirmation criteria according to recommended tolerances in 2002/657/EC Commission Decision, that was in place when this work was carried out71. Samples were considered contaminated when analytes were detected according to identification criteria by liquid chromatography tandem mass spectrometry and the concentration values exceeded the limits of detection (LOD).

Water samples were filtered through 0.22 µm cellulose acetate membranes (Millipore, USA) under aseptic conditions. All experiments were performed using kits and internal controls. DNA was extracted from the filters using the PowerWater DNA Isolation Kit (Qiagen Science, USA). For the preparation of the amplicon library, the DNA was quantified using a Qubit 2.0 Fluorometer (ThermoFisher Scientific, USA), and samples were diluted to achieve the concentration of 5 ng μL−1 per sample. The V4 hypervariable region of the 16S rRNA gene was amplified by PCR using the primers 16Sf (5′- GTGCCAGCMGCCGCGGTAA-3′) and 16Sr (5′- GGACTACHVGGGTWTCTAAT-3′) with the appropriate barcodes and adapters72. PCR products were purified using the ChargeSwitch™ PCR Clean-Up Kit (ThermoFisher Scientific, USA). Each individual sample library was diluted to 4 nM and then pooled and paired-end sequenced on a MiSeq system (Illumina Inc. USA), using the 500 cycles MiSeq Reagent v2 Kit.

Quality analysis of raw reads was carried out with FastQC software73 and the filtering of sequences with an average quality equal to or greater than 20 was performed by the PRINSEQ program74. Data analysis was performed using QIIME (Quantitative Insights Into Microbial Ecology) 1.9.175. The data were compared with the SILVA Ribosomal RNA database (non-redundant) 132 release72 with a maximum e-value of 1e-5, and a minimum identity of 99%, which generated a table with taxonomic groups. Statistical analyzes such as Alpha and Beta diversity were calculated using the MicrobiomeAnalyst web platform (https://www.microbiomeanalyst.ca/)76,77. The diversity of the bacterial communities was assessed using the Chao1 index and Simpson's index calculated for OTUs with the evolutionary distance of 0.01 (or 99% 16S rRNA gene sequence similarity). Principal coordinate analysis (PCoA) between bacterial communities present in the Guandu, São João watersheds and in drinking water was built using the Jaccard method with PERMANOVA and using the bacterial OTUs.

Plasmid DNA (pDNA) was extracted from the filters by alkaline lysis using the Plasmid Mini Kit (Qiagen Science, USA) according to the manufacturer's protocol. pDNA was precipitated with isopropanol and washed with 70% ethanol. To remove possible traces of genomic DNA the precipitate was treated with ATP-dependent Plasmid Safe DNase (Epicentre, USA) according to manufacturer's instructions46. The pDNA was quantified using a Qubit 2.0 Fluorometer (ThermoFisher ScientificTM, USA) according to the manufacturer's manual.

A pDNA sequencing library was prepared using the Nextera XT DNA Library Prep Kit (Illumina Inc. USA) following the manufacturer's recommendations. Paired-end sequencing was performed with the 600 cycles Miseq Reagent Kit v.3 on the MiSeq platform (Illumina Inc. USA). Sequence quality checks were performed with the FastQC software73 and sequence filtering with an average quality of 20 or higher was performed by PRINSEQ78. Contigs were assembled with SPAdes version 3.1379,80 using the metaSPAdes pipeline.

The sequences were analyzed by the MG-RAST (Meta Genome—Rapid Annotation using Subsystem Technology) platform81, where the annotation can be viewed in several different categories, including subsystems. A subsystem can be understood as a set of functional roles that implement a certain biological or structural process82. The subsystems are classified into hierarchical levels, so that level 1 includes general catabolic and anabolic functions (for example, DNA metabolism), and levels 2 and 3 contain more specific pathways, such as antimicrobial resistance and other compounds83.

In addition, the sequences were compared against the Comprehensive Antibiotic Resistance Database (CARD) database84,85 with DIAMOND86. Only alignments with an e-value < 1e5, coverage > 60% and amino acid identity > 30% were considered87.

This article does not contain any studies with human or animals performed by any of the authors.

All data employed in this report are available in GenBank under BioProject PRJNA812588 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA812588). In addition, metagenome sequence data are available on MG-RAST under accession numbers mgm4919709.3, mgm4919786.3, mgm4919818.3.

Ilyas, M. et al. Environmental and health impacts of industrial wastewater effluents in Pakistan: A review. Rev. Environ. Health 34, 171–186 (2019).

Article MathSciNet PubMed CAS Google Scholar

Roy, M. A. et al. A Metagenomic approach to evaluating surface water quality in Haiti. Int. J. Environ. Res. Public Health 15, 2211 (2018).

Article PubMed Central CAS Google Scholar

Liu, M. & Lu, J. Support vector machine-an alternative to artificial neuron network for water quality forecasting in an agricultural nonpoint source polluted river?. Environ. Sci. Pollut. Res. 21, 11036–11053 (2014).

Article CAS Google Scholar

Yao, H., Qian, X., Yin, H., Gao, H. & Wang, Y. Regional risk assessment for point source pollution based on a water quality model of the Taipu River, China. Risk Anal. 35, 265–277 (2015).

Article PubMed Google Scholar

Dey, S., Bano, F. & Malik, A. Pharmaceuticals and personal care product (PPCP) contamination—a global discharge inventory. in Pharmaceuticals and Personal Care Products: Waste Management and Treatment Technology 1–26 (Elsevier, 2019). doi:https://doi.org/10.1016/B978-0-12-816189-0.00001-9.

Richardson, S. D. & Ternes, T. A. Water analysis: Emerging contaminants and current issues. Anal. Chem. 86, 2813–2848 (2014).

Article PubMed CAS Google Scholar

Felis, E. et al. Antimicrobial pharmaceuticals in the aquatic environment—occurrence and environmental implications. Eur. J. Pharmacol. 866, 172813 (2020).

Article PubMed CAS Google Scholar

Nassiri Koopaei, N. & Abdollahi, M. Health risks associated with the pharmaceuticals in wastewater. Daru 25, 9 (2017).

Article PubMed PubMed Central Google Scholar

Gillings, M. R. et al. Using the class 1 integron-integrase gene as a proxy for anthropogenic pollution. ISME J. 9, 1269–1279 (2015).

Article PubMed CAS Google Scholar

Baquero, F., Martínez, J.-L. & Cantón, R. Antibiotics and antibiotic resistance in water environments. Curr. Opin. Biotechnol. 19, 260–265 (2008).

Article PubMed CAS Google Scholar

Martínez, J. L. Antibiotics and antibiotic resistance genes in natural environments. Science 321, 365–367 (2008).

Article ADS PubMed Google Scholar

Grenni, P., Ancona, V. & Barra Caracciolo, A. Ecological effects of antibiotics on natural ecosystems: A review. Microchem. J. 136, 25–39 (2018).

Article CAS Google Scholar

Kraemer, S. A., Ramachandran, A. & Perron, G. G. Antibiotic pollution in the environment: From microbial ecology to public policy. Microorganisms 7, 180 (2019).

Article PubMed Central CAS Google Scholar

Subirats, J. et al. Emerging contaminants and nutrients synergistically affect the spread of class 1 integron-integrase (intI1) and sul1 genes within stable streambed bacterial communities. Water Res. 138, 77–85 (2018).

Article PubMed CAS Google Scholar

Partridge, S. R., Kwong, S. M., Firth, N. & Jensen, S. O. Mobile genetic elements associated with antimicrobial resistance. Clin. Microbiol. Rev. https://doi.org/10.1128/CMR.00088-17 (2018).

Article PubMed PubMed Central Google Scholar

Wein, T., Hülter, N. F., Mizrahi, I. & Dagan, T. Emergence of plasmid stability under non-selective conditions maintains antibiotic resistance. Nat. Commun. 10, 2595 (2019).

Article ADS PubMed PubMed Central Google Scholar

Walker, A. Welcome to the plasmidome. Nat. Rev. Microbiol. 10, 379 (2012).

Article PubMed CAS Google Scholar

Norman, A. et al. An improved method for including upper size range plasmids in metamobilomes. PLoS ONE 9, e104405 (2014).

Article ADS PubMed PubMed Central Google Scholar

Aminov, R. I. & Mackie, R. I. Evolution and ecology of antibiotic resistance genes. FEMS Microbiol. Lett. 271, 147–161 (2007).

Article PubMed CAS Google Scholar

Rahman, M. H., Nonaka, L., Tago, R. & Suzuki, S. Occurrence of two genotypes of tetracycline (TC) resistance gene tet(M) in the TC-resistant bacteria in marine sediments of Japan. Environ. Sci. Technol. 42, 5055–5061 (2008).

Article ADS PubMed CAS Google Scholar

Unc, A. & Goss, M. J. Transport of bacteria from manure and protection of water resources. Appl. Soil Ecol. 25, 1–18 (2004).

Article Google Scholar

Alawi, M., Torrijos, T. V. & Walsh, F. Plasmid-mediated antimicrobial resistance in drinking water. Environ. Adv. https://doi.org/10.1016/j.envadv.2022.100191 (2022).

Article Google Scholar

Kümmerer, K. Antibiotics in the aquatic environment–a review–part I. Chemosphere 75, 417–434 (2009).

Article ADS PubMed Google Scholar

Kümmerer, K. Antibiotics in the aquatic environment—A review—Part II. Chemosphere https://doi.org/10.1016/j.chemosphere.2008.12.006 (2009).

Article PubMed Google Scholar

Cardoso, O., Porcher, J.-M. & Sanchez, W. Factory-discharged pharmaceuticals could be a relevant source of aquatic environment contamination: Review of evidence and need for knowledge. Chemosphere 115, 20–30 (2014).

Article ADS PubMed CAS Google Scholar

He, Y. et al. Antibiotic resistance genes from livestock waste: occurrence, dissemination, and treatment. npj Clean Water 3, 4 (2020).

Article Google Scholar

Topp, E., Renaud, J., Sumarah, M. & Sabourin, L. Reduced persistence of the macrolide antibiotics erythromycin, clarithromycin and azithromycin in agricultural soil following several years of exposure in the field. Sci. Total Environ. https://doi.org/10.1016/j.scitotenv.2016.03.210 (2016).

Article PubMed Google Scholar

Segura, P. A., François, M., Gagnon, C. & Sauvé, S. Review of the occurrence of anti-infectives in contaminated wastewaters and natural and drinking waters. Environ. Health Perspect. 117, 675–684 (2009).

Article PubMed PubMed Central CAS Google Scholar

Benotti, M. J. et al. Pharmaceuticals and endocrine disrupting compounds in U.S. drinking water. Environ. Sci. Technol. 43, 597–603 (2009).

Article ADS PubMed CAS Google Scholar

Ben, Y. et al. Human health risk assessment of antibiotic resistance associated with antibiotic residues in the environment: A review. Environ. Res. 169, 483–493 (2018).

Article PubMed Google Scholar

Vaz-Moreira, I., Nunes, O. C. & Manaia, C. M. Bacterial diversity and antibiotic resistance in water habitats: Searching the links with the human microbiome. FEMS Microbiol. Rev. 38, 761–778 (2014).

Article PubMed CAS Google Scholar

Feng, B.-W. et al. Bacterial diversity of water and sediment in the Changjiang estuary and coastal area of the East China Sea. FEMS Microbiol. Ecol. 70, 80–92 (2009).

Article PubMed Google Scholar

Qin, Y. et al. Bacterial abundance and diversity in pond water supplied with different feeds. Sci. Rep. https://doi.org/10.1038/srep35232 (2016).

Article PubMed PubMed Central Google Scholar

Forsberg, K. J. et al. Bacterial phylogeny structures soil resistomes across habitats. Nature https://doi.org/10.1038/nature13377 (2014).

Article PubMed PubMed Central Google Scholar

Stanish, L. F. et al. Factors influencing bacterial diversity and community composition in municipal drinking waters in the Ohio River Basin, USA. PLoS ONE 11, e0157966 (2016).

Article PubMed PubMed Central Google Scholar

Medeiros, J. D. et al. Comparative metagenome of a stream impacted by the urbanization phenomenon. Braz. J. Microbiol. https://doi.org/10.1016/j.bjm.2016.06.011 (2016).

Article PubMed PubMed Central Google Scholar

Paul, Michael J. & Meyer, Judy L. Streams in the urban landscape. In Urban Ecology (eds Marzluff, John M. et al.) (Springer, 2008).

Google Scholar

Nogueira, I. S., Nabout, J. C., Oliveira, J. E. & Silva, K. D. Diversidade (alfa, beta e gama) da comunidade fitoplanctônica de quatro lagos artificiais urbanos do município de Goiânia, GO. Hoehnea 35, 219–233 (2008).

Article Google Scholar

Tuomisto, H. A consistent terminology for quantifying species diversity? Yes, it does exist. Oecologia 164, 853–860 (2010).

Article ADS PubMed Google Scholar

Cabral, L. et al. Anthropogenic impact on mangrove sediments triggers differential responses in the heavy metals and antibiotic resistomes of microbial communities. Environ. Pollut. https://doi.org/10.1016/j.envpol.2016.05.078 (2016).

Article PubMed Google Scholar

Conejo, M. C., García, I., Martínez-Martínez, L., Picabea, L. & Pascual, A. Zinc eluted from siliconized latex urinary catheters decreases OprD expression, causing carbapenem resistance in Pseudomonas aeruginosa. Antimicrob. Agents Chemother. 47, 2313–2315 (2003).

Article PubMed PubMed Central CAS Google Scholar

Perron, K. et al. CzcR-CzcS, a two-component system involved in heavy metal and carbapenem resistance in Pseudomonas aeruginosa. J. Biol. Chem. 279, 8761–8768 (2004).

Article PubMed CAS Google Scholar

Guo, X., Li, J., Yang, F., Yang, J. & Yin, D. Prevalence of sulfonamide and tetracycline resistance genes in drinking water treatment plants in the Yangtze River Delta, China. Sci. Total Environ. https://doi.org/10.1016/j.scitotenv.2014.06.035 (2014).

Article PubMed Google Scholar

Furukawa, T., Jikumaru, A., Ueno, T. & Sei, K. Inactivation effect of antibiotic-resistant gene using chlorine disinfection. Water (Switzerland) https://doi.org/10.3390/w9070547 (2017).

Article Google Scholar

Khan, H., Miao, X., Liu, M., Ahmad, S. & Bai, X. Behavior of last resort antibiotic resistance genes (mcr-1 and blaNDM-1) in a drinking water supply system and their possible acquisition by the mouse gut flora. Environ. Pollut. https://doi.org/10.1016/j.envpol.2019.113818 (2020).

Article PubMed PubMed Central Google Scholar

Li, A.-D., Li, L.-G. & Zhang, T. Exploring antibiotic resistance genes and metal resistance genes in plasmid metagenomes from wastewater treatment plants. Front. Microbiol. 6, 1025 (2015).

Article PubMed PubMed Central Google Scholar

Nesme, J. et al. Large-scale metagenomic-based study of antibiotic resistance in the environment. Curr. Biol. https://doi.org/10.1016/j.cub.2014.03.036 (2014).

Article PubMed Google Scholar

Taggar, G., Attiq Rehman, M., Boerlin, P. & Diarra, M. Molecular epidemiology of carbapenemases in enterobacteriales from humans, animals, food and the environment. Antibiotics 9, 693 (2020).

Article PubMed Central CAS Google Scholar

Liu, Y.-Y. et al. Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: A microbiological and molecular biological study. Lancet. Infect. Dis. 16, 161–168 (2016).

Article PubMed Google Scholar

Fernandes, M. R. et al. Silent dissemination of colistin-resistant Escherichia coli in South America could contribute to the global spread of the mcr-1 gene. Eurosurveillance https://doi.org/10.2807/1560-7917.ES.2016.21.17.30214 (2016).

Article PubMed Google Scholar

Kempf, I., Jouy, E. & Chauvin, C. Colistin use and colistin resistance in bacteria from animals. Int. J. Antimicrob. Agents https://doi.org/10.1016/j.ijantimicag.2016.09.016 (2016).

Article PubMed Google Scholar

Catry, B. et al. Use of colistin-containing products within the European Union and European Economic Area (EU/EEA): Development of resistance in animals and possible impact on human and animal health. Int. J. Antimicrob. Agents 46, 297–306 (2015).

Article PubMed CAS Google Scholar

Rossolini, G. M. et al. Metallo-β-lactamase producers in environmental microbiota: New molecular class B enzyme in Janthinobacterium lividum. Antimicrob. Agents Chemother. https://doi.org/10.1128/AAC.45.3.837-844.2001 (2001).

Article PubMed PubMed Central Google Scholar

Marathe, N. P. et al. Untreated urban waste contaminates Indian river sediments with resistance genes to last resort antibiotics. Water Res. https://doi.org/10.1016/j.watres.2017.07.060 (2017).

Article PubMed Google Scholar

Piedra-Carrasco, N. et al. Carbapenemase-producing enterobacteriaceae recovered from a Spanish river ecosystem. PLoS ONE https://doi.org/10.1371/journal.pone.0175246 (2017).

Article PubMed PubMed Central Google Scholar

Islam, M. A. et al. Environmental spread of New Delhi metallo-β- lactamase-1-producing multidrug-resistant bacteria in Dhaka, Bangladesh. Appl. Environ. Microbiol. https://doi.org/10.1128/AEM.00793-17 (2014).

Article PubMed Google Scholar

Mombini, S., Rezatofighi, S. E., Kiyani, L. & Motamedi, H. Diversity and metallo-β-lactamase-producing genes in Pseudomonas aeruginosa strains isolated from filters of household water treatment systems. J. Environ. Manage. https://doi.org/10.1016/j.jenvman.2018.10.068 (2019).

Article PubMed Google Scholar

Ashbolt, N. J. et al. Human health risk assessment (HHRA) for environmental development and transfer of antibiotic resistance. Environ. Health Perspect. https://doi.org/10.1289/ehp.1206316 (2013).

Article PubMed PubMed Central Google Scholar

Seiler, C. & Berendonk, T. U. Heavy metal driven co-selection of antibiotic resistance in soil and water bodies impacted by agriculture and aquaculture. Front. Microbiol. 3, 399 (2012).

Article PubMed PubMed Central Google Scholar

Zhang, Y., Gu, A. Z., He, M., Li, D. & Chen, J. Subinhibitory concentrations of disinfectants promote the horizontal transfer of multidrug resistance genes within and across genera. Environ. Sci. Technol. 51, 570–580 (2017).

Article ADS PubMed CAS Google Scholar

Li, D., Zeng, S., He, M. & Gu, A. Z. Water disinfection byproducts induce antibiotic resistance-role of environmental pollutants in resistance phenomena. Environ. Sci. Technol. 50, 3193–3201 (2016).

Article ADS PubMed PubMed Central CAS Google Scholar

Ding, C. et al. Enhanced uptake of antibiotic resistance genes in the presence of nanoalumina. Nanotoxicology 10, 1051–1060 (2016).

Article PubMed CAS Google Scholar

Fish, K. E., Reeves-McLaren, N., Husband, S. & Boxall, J. Author Correction: Uncharted waters: The unintended impacts of residual chlorine on water quality and biofilms. NPJ Biofilms Microbiomes 8, 55 (2022).

Article PubMed PubMed Central Google Scholar

Bergeron, S., Boopathy, R., Nathaniel, R., Corbin, A. & LaFleur, G. Presence of antibiotic resistant bacteria and antibiotic resistance genes in raw source water and treated drinking water. Int. Biodeterior. Biodegrad. 102, 370–374 (2015).

Article CAS Google Scholar

Xi, C. et al. Prevalence of antibiotic resistance in drinking water treatment and distribution systems. Appl. Environ. Microbiol. 75, 5714–5718 (2009).

Article ADS PubMed PubMed Central CAS Google Scholar

Colomer-Lluch, M., Jofre, J. & Muniesa, M. Antibiotic resistance genes in the bacteriophage DNA fraction of environmental samples. PLoS ONE 6, e17549 (2011).

Article ADS PubMed PubMed Central CAS Google Scholar

Britto, A. L., Maiello, A. & Quintslr, S. Water supply system in the Rio de Janeiro Metropolitan Region: Open issues, contradictions, and challenges for water access in an emerging megacity. J. Hydrol. 573, 1007–1020 (2019).

Article Google Scholar

Bacha, L. et al. Risk of Collapse in water quality in the Guandu River (Rio de Janeiro, Brazil). Microb. Ecol. 84, 314–324 (2022).

Article PubMed CAS Google Scholar

U. S. Environmental Protection Agency. Method 1694 : Pharmaceuticals and Personal Care Products in Water, Soil, Sediment, and Biosolids by HPLC/MS/MS. EPA Method (2007) doi:https://doi.org/10.1002/etc.3451.

Monteiro, M. et al. Development and validation of liquid chromatography-tandem mass spectrometry methods for determination of beta-lactams, macrolides, fluoroquinolones, sulfonamides and tetracyclines in surface and drinking water from Rio de Janeiro, Brazil. J. Braz. Chem. Soc. https://doi.org/10.21577/0103-5053.20170203 (2017).

Article Google Scholar

Communities European. Commission Decision of 12 August 2002 implementing Council Directive 96/23/EC concerning the performance of analytical methods and the interpretation of results. 8–36 (2002).

Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K. & Schloss, P. D. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl. Environ. Microbiol. 79, 5112–5120 (2013).

Article ADS PubMed PubMed Central CAS Google Scholar

Andrews, S. FastQC: A quality control tool for high throughput sequence data. Babraham Bioinformatics http://www.bioinformatics.babraham.ac.uk/projects/ (2010) doi:citeulike-article-id:11583827.

Schmieder, R. & Edwards, R. Insights into antibiotic resistance through metagenomic approaches. Future Microbiol. 7, 73–89 (2012).

Article PubMed CAS Google Scholar

Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336 (2010).

Article PubMed PubMed Central CAS Google Scholar

Dhariwal, A. et al. MicrobiomeAnalyst: A web-based tool for comprehensive statistical, visual and meta-analysis of microbiome data. Nucleic Acids Res. 45, W180–W188 (2017).

Article PubMed PubMed Central CAS Google Scholar

Chong, J., Liu, P., Zhou, G. & Xia, J. Using MicrobiomeAnalyst for comprehensive statistical, functional, and meta-analysis of microbiome data. Nat. Protoc. https://doi.org/10.1038/s41596-019-0264-1 (2020).

Article PubMed Google Scholar

Schmieder, R. & Edwards, R. Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863–864 (2011).

Article PubMed PubMed Central CAS Google Scholar

Bankevich, A. et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477 (2012).

Article MathSciNet PubMed PubMed Central CAS Google Scholar

Nurk, S., Meleshko, D., Korobeynikov, A. & Pevzner, P. A. metaSPAdes: A new versatile metagenomic assembler. Genome Res. 27, 824–834 (2017).

Article PubMed PubMed Central CAS Google Scholar

Keegan, K. P., Glass, E. M. & Meyer, F. MG-RAST, a metagenomics service for analysis of microbial community structure and function. Methods Mol. Biol. 1399, 207–233 (2016).

Article PubMed CAS Google Scholar

Overbeek, R. et al. The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Res. https://doi.org/10.1093/nar/gki866 (2005).

Article PubMed PubMed Central Google Scholar

Dinsdale, E. A. et al. Functional metagenomic profiling of nine biomes. Nature https://doi.org/10.1038/nature06810 (2008).

Article PubMed Google Scholar

Jia, B. et al. CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 45, D566–D573 (2017).

Article PubMed CAS Google Scholar

Alcock, B. P. et al. CARD 2020: Antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. https://doi.org/10.1093/nar/gkz935 (2020).

Article PubMed Google Scholar

Buchfink, B., Reuter, K. & Drost, H.-G. Sensitive protein alignments at tree-of-life scale using DIAMOND. Nat. Methods 18, 366–368 (2021).

Article PubMed PubMed Central CAS Google Scholar

Dias, M. F. et al. Exploring the resistome, virulome and microbiome of drinking water in environmental and clinical settings. Water Res. 174, 115630 (2020).

Article PubMed CAS Google Scholar

Download references

This work was funded by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Fundação Oswaldo Cruz, Conselho Nacional de Desenvolvimento Científico e Tecnológico.

Instituto Nacional de Controle de Qualidade Em Saúde INCQS/FIOCRUZ, Oswaldo Cruz Foundation, Rio de Janeiro, RJ, 4365, Brazil

Kayo Bianco, Beatriz Oliveira de Farias, Andressa Silva Gonçalves-Brito, Ana Paula Alves do Nascimento, Mariana Magaldi, Kaylanne Montenegro, Claudia Flores, Samara Oliveira, Mychelle Alves Monteiro, Bernardete Ferraz Spisso, Mararlene Ulberg Pereira, Rosana Gomes Ferreira & Maysa Mandetta Clementino

Rio de Janeiro State University, Rio de Janeiro, Brazil

Rodolpho Mattos Albano & Alexander Machado Cardoso

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

K.B. collected the samples, performed experimental works and data analysis, wrote the manuscript; B.O.F. performed experimental works and data analysis; A.S.G.B. performed experimental works and data analysis; A.P.A.N. contributed to the data analysis, and the writing of the paper (review and editing); M.M. performed experimental works and data analysis; K.M. contributed to the data analysis and the writing of the paper; C.F. performed experimental works and data analysis; S.O. performed experimental works and data analysis; M.A.M. contributed to the experimental work; B.F.S. data analysis, and the writing of the paper (review and editing); R.M.A., B.F.S., M.U.P., R.G.F. contributed to the experimental work, data analysis, and the writing of the paper (review and editing); A.M.C. contributed to the data analysis and the writing of the paper (review and editing); M.M.C. was responsible for funding acquisition, project administration, supervision, and contributed to the writing of the paper (review and editing). All authors read and approved the manuscript.

Correspondence to Kayo Bianco.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Bianco, K., de Farias, B.O., Gonçalves-Brito, A.S. et al. Mobile resistome of microbial communities and antimicrobial residues from drinking water supply systems in Rio de Janeiro, Brazil. Sci Rep 12, 19050 (2022). https://doi.org/10.1038/s41598-022-21040-7

Download citation

Received: 29 March 2022

Accepted: 22 September 2022

Published: 09 November 2022

DOI: https://doi.org/10.1038/s41598-022-21040-7

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

SHARE