Laura Sala-Comorera1
Tristan M. Nolan1
Liam J. Reynolds1
Anjan Venkatesh1
Lily Cheung1
Niamh A. Martin1Jayne H. Stephens1Aurora Gitto1Gregory M. P. O’Hare2
John J. O’Sullivan3
Wim G. Meijer1*- 1UCD School of Biomolecular and Biomedical Science, UCD Earth Institute, UCD Conway Institute, University College Dublin, Dublin, Ireland
- 2UCD School of Computer Science, UCD Earth Institute, University College Dublin, Dublin, Ireland
- 3UCD School of Civil Engineering, UCD Dooge Centre for Water Resources Research, UCD Earth Institute, University College Dublin, Dublin, Ireland
Fecal pollution of surface water may introduce bacteria and bacteriophages harboring antibiotic resistance genes (ARGs) into the aquatic environment. Watercourses discharging into the marine environment, especially close to designated bathing waters, may expose recreational users to fecal pollution and therefore may increase the likelihood that they will be exposed to ARGs. This study compares the bacterial and bacteriophage ARG profiles of two rivers (River Tolka and Liffey) and two small urban streams (Elm Park and Trimleston Streams) that discharge close to two marine bathing waters in Dublin Bay. Despite the potential differences in pollution pressures experienced by these waterways, microbial source tracking analysis showed that the main source of pollution in both rivers and streams in the urban environment is human contamination. All ARGs included in this study, blaTEM, blaSHV, qnrS, and sul1, were present in all four waterways in both the bacterial and bacteriophage fractions, displaying a similar ARG profile. We show that nearshore marine bathing waters are strongly influenced by urban rivers and streams discharging into these, since they shared a similar ARG profile. In comparison to rivers and streams, the levels of bacterial ARGs were significantly reduced in the marine environment. In contrast, the bacteriophage ARG levels in freshwater and the marine were not significantly different. Nearshore marine bathing waters could therefore be a potential reservoir of bacteriophages carrying ARGs. In addition to being considered potential additional fecal indicators organism, bacteriophages may also be viewed as indicators of the spread of antimicrobial resistance.
Introduction
Microbial antibiotic resistance is a severe threat to public health, resulting in failure to treat a range of infections, in extended hospital treatment and in increased healthcare costs (Stewardson et al., 2016; Cassini et al., 2019; Roope et al., 2019; Jit et al., 2020). Furthermore, it is estimated that multidrug resistant pathogens will lead to 10 million deaths by 2050 (OECD, 2018; WHO, 2019). Initially, the problem was mainly approached from a clinical perspective. However, more recently, the importance of the environment as a major contributor to the spread of antimicrobial resistance in the human and animal population has become increasingly clear (Allcock et al., 2017; Leonard et al., 2018; Hernando-Amado et al., 2020). Successful management strategies and policies to combat the increase in antimicrobial resistance will therefore have to be based on a One Health approach, which recognizes the connectivity between animal, human, and environmental health (Hernando-Amado et al., 2019; Van Bruggen et al., 2019).
Antibiotic resistance in bacteria may arise from chromosomal mutations, however, for most types of antimicrobial resistance, the acquisition of antibiotic resistance genes (ARGs) mediated by horizontal gene transfer mechanisms is a more common mechanism (Andersson and Hughes, 2010; Huddleston, 2014; Baquero et al., 2019). Transduction by bacteriophages is an important mechanism in spreading ARGs within a microbial population (Muniesa et al., 2013b; Balcazar, 2014; Brown-Jaque et al., 2015; Balcázar, 2018; Maganha De et al., 2021). The importance of bacteriophages in spreading ARGs and virulence genes, was recently underscored by the discovery of lateral transduction, which is the basis of genomic hypermobility (Chen et al., 2018; Chiang et al., 2019).
Environments characterized by high levels of fecal matter such as sewage, animal slurry, sludge, and effluent of wastewater treatment plants are hotspots of antibiotic resistant bacteria and phages harboring ARGs, forming an ideal environment for horizontal gene transfer to occur (Calero-Cáceres et al., 2014; Quirós et al., 2014; Ross and Topp, 2015; Calero-Cáceres and Muniesa, 2016; Guo et al., 2017; Yang et al., 2021; Zieliński et al., 2021). It has been shown that transfection and transduction of ARGs from environmental phages conferred resistance to the recipient bacteria (Battaglioli et al., 2011; Colomer-Lluch et al., 2011; Gunathilaka et al., 2017; Wang et al., 2018b; Yang et al., 2021).
The effluent of wastewater treatment plants is an important route by which ARGs may enter the aquatic environment (Rizzo et al., 2013; Marti et al., 2014b; Rodriguez-Mozaz et al., 2015; Lekunberri et al., 2017; Zhou et al., 2020; Nguyen et al., 2021). Sewerage misconnections and leaking septic tanks may be a source of untreated sewage entering rivers and other waterbodies (Kay et al., 2008; Hinojosa et al., 2020; Reynolds et al., 2021). Furthermore, agricultural land run-off may also negatively impact water quality and introduce ARGs into waterbodies (Unc and Goss, 2004; Ballesté et al., 2020). Watercourses discharging into the marine environment, especially close to designated bathing waters, may expose the users to fecal pollution and therefore may increase the likelihood that they will be exposed to ARGs (Molina et al., 2014; Leonard et al., 2018; Ahmed et al., 2020; Reynolds et al., 2020; Sala-Comorera et al., 2021b).
This study focuses on the bacterial and bacteriophage ARG profiles of two rivers entering a large urban environment and two small, completely urban, streams that discharge close to two marine bathing waters. We show that the ARG profiles of these streams and rivers in an urban environment are highly similar and have a strong impact on nearshore marine bathing waters. In contrast to bacterial ARGs, ARGs associated with bacteriophages appear to persist in the nearshore marine environment.
Materials and Methods
Site Location
Dublin, the capital of Ireland, is a coastal city on the Irish Sea with 560,000 inhabitants. Near 1,905,000 people live within the greater Dublin area and around sixty watercourses, ranging from major rivers to small streams, which discharge into Dublin Bay. Dublin Bay, which is a UNESCO biosphere, has three designated bathing waters: Dollymount, Sandymount, and Merrion Strands.
In this study, two rivers (River Liffey and River Tolka), two streams (Elm Park Stream and Trimleston Stream), and two bathing waters (Merrion Strand and Sandymount Strand) were selected (Figure 1 and Supplementary Table 1). The River Liffey rises in a pristine area near Kippure in the Wicklow Mountains and flows through agricultural land to reaching Dublin city. The River Liffey has a 125 km course and it is the main river flowing into Dublin Bay (Sweeney, 2017). The River Tolka is 33 km long, and it is the second largest river by flow in Dublin Bay. The river rises near Culmullin Cross Road and flows through agricultural and industrial land into the north of Dublin city (Sweeney, 2017). The Elm Park Stream and the Trimleston Stream catchment is completely urban areas with a population of 40,000 people. The length of the streams are 3.8 and 1.7 km, respectively, and the depth is less than 10 cm. Both streams discharge close to a designated bathing area in Dublin Bay.
Figure 1. Location of the sampling sites selected in this study. The sampling stations are indicated by red squares (rivers and streams) and stars (bathing waters). The inset box shows the location of Dublin Bay in Ireland indicated by a yellow square.
Grab samples were collected at the tidal limit at a depth of 10–20 cm, every 5 weeks over 15 months (from September 2018 to November 2019). Bathing water samples were taken during high tide at Merrion and Sandymount Strands. A total of 85 samples collected and stored at 4°C before being processed within 6 h.
Enumeration of Fecal Indicator Organisms
The levels of Escherichia coli and intestinal enterococci were determined by membrane filtration. Water samples were filtered through 0.45 μm pore size nitrocellulose membranes (Thermo Scientific) and placed on Tryptone Bile X-Glucuronide agar (Sigma-Aldrich) at 37°C for 4 h, followed by an incubation at 44°C for 18 h to enumerate E. coli (ISO, 2001). Intestinal enterococci were enumerated by incubating the membrane on Slanetz and Bartley agar (Oxoid) at 37°C for 48 h. After incubation, membranes were transferred into Bile Aesculin agar at 44°C for 2 h to confirm positive intestinal enterococci colonies (ISO, 2000).
Electron Microscopy
River samples (100 ml) were concentrated by ultrafiltration using 100 kDa Amicon Ultra-15 Centrifugal Filter units (Millipore) and 5 μl of the concentrated samples were pipetted onto a 200-mesh copper grid coated with formvar. Samples were negatively stained with 5 μl of 2% uranyl acetate stain and incubated for 2 min. The grids were imaged using a transmission electron microscope Tecnai G2 (FEI Tecnai) operating at 120 kV. Untreated sewage (50 ml) collected in the influent of a local wastewater treatment plant was used as a positive control.
DNA Extraction From the Bacterial Fraction
DNA of the water samples was extracted after concentrating by filtering 100 ml through 0.22 μm mixed cellulose ester membrane filters. The filters were then transferred in 500 μl of GITC buffer [5 M guanidine thiocyanate, 100 mM EDTA (pH 8), and 0.5% sarkosyl] and stored at −20°C. DNA was extracted using the DNeasy Blood and Tissue kit (Qiagen) with some modifications as reported previously (Gourmelon et al., 2007). The DNA was eluted in a final volume of 70 μl.
DNA Extraction From the Bacteriophage Fraction
DNA of bacteriophages present in water was extracted using the protocol described by Colomer-Lluch et al. (2011). Briefly, 100 ml of water was passed through 0.22 μm low protein binding polyethersulfone filters. The filtrate was concentrated 200-fold using 100 kDa Amicon Ultra-15 Centrifugal Filter units (Millipore), and subjected to chloroform extraction in a 1:1 (v/v) ratio, followed by a DNAse (100 U/ml) treatment at 37°C for 1 h and 10 min at 80°C. At this point, an aliquot of 10 μl was collected as a control to ensure the complete removal of free DNA. Phage particles were subjected to proteinase K digestion (0.2 mg/ml) for 1 h at 56°C followed by phenol/chloroform extraction and ethanol precipitation. The resulting DNA was dissolved into 20 μl water.
Quantification of Gene Target
Primers to amplify the ARGs blaTEM, blaSHV, qnrS, and sul1 and the human (HF183) and ruminant (CF128) microbial source tracking markers, as well as the cycling conditions are described in Supplementary Table 2. The MST markers and ARGs were quantified as previously described (Ballesté et al., 2020; Reynolds et al., 2020). Standard curves were generated using linearized cloned standards between 100 and 106 gene copies to quantify target gene levels in each sample (Ballesté et al., 2020; Reynolds et al., 2020). The limit of detection of each assay was determined as the lowest concentration of DNA target detected in 95% or more of replicates, whereas the limit of quantification was determined as the lowest concentration of DNA quantified within 0.5 SD of the log10 concentration (Supplementary Table 2; Blanchard et al., 2012; AFNOR, 2015). All samples and negative controls were analyzed in duplicate in each 96-well plate. The absence of non-packaged DNA in the bacteriophage extraction protocol was verified with 16S rRNA gene amplification by PCR (AllTaq Master Mix Kit, Qiagen) and ARGs genes by qPCR. Only negative samples were used for the subsequent analysis. MST markers and ARGs concentrations were expressed as gene copies per 100 ml (GC/100 ml). The amplification efficiency of each reaction was calculated using the E = 10(1/slope) − 1 equation (Rutledge, 2003).
Data Analysis
The non-parametric Mann–Whitney paired-test and Kruskal–Wallis test with Dunn’s post hoc analysis was used to assess significant differences between microbial source tracking markers and ARGs between waterbodies. The values of the qPCR targets were log10 transformed and Spearman correlation was used to identify relationships between variables. A significance cut-off of p ≤ 0.05 was used for all analyses. Statistical analysis was carried out using GraphPad Prism 9.1.0. software (GraphPad Software).
Results
Fecal Contamination of Rivers, Urban Streams, and Bathing Waters
All water samples (n = 85) were positive for E. coli and intestinal enterococci (Figure 2). The River Tolka and the River Liffey, which flow through agricultural areas before reaching the city, had similar median levels of fecal indicator organisms as the small urban Elm Park and Trimleston streams. The fecal indicators in the latter varied more than three orders of magnitude. In contrast, the fecal indicator levels were less variable for the River Liffey during this 15-month period, which had the lowest median concentration for the fecal indicators. In general, all rivers and streams received fecal contamination along their course, independently of the catchment area.
Figure 2. Levels of the fecal indicators Escherichia coli, intestinal enterococci (I.E.). The dots show the concentration of each sample and the median are indicated as a horizontal line.
The median fecal indicator levels in Merrion and Sandymount Strands were similar with individual values varying by up to four orders of magnitude. Some of the bathing water samples analyzed exceeded the 90th percentile value for sufficient water quality parameters for coastal waters (≥500 CFU/100 ml for E. coli and ≥185 CFU/100 ml for intestinal enterococci) according to the European bathing water quality Directive (EU, 2006).
Fecal Contamination of Rivers and Streams Is Predominantly Human in Nature
The Bacteroidales human (HF183) and ruminant (CF128) markers were deployed to determine the biological origins of pollution (Figure 3). The River Liffey catchment is larger than that of the River Tolka and has more agricultural land use. Therefore, as expected, more River Liffey (80%) than River Tolka (66%) samples were positive for the ruminant marker. The levels of the ruminant marker reached values of 5.1 × 105 and 3.4 × 105 GC/100 ml in River Liffey and River Tolka, respectively. In contrast the marker was only sporadically detected in the Elm Park stream (40%) and in only one sample in the Trimleston stream. When the CF128 marker tested positive in Elm Park stream samples, concentrations up to 2.8 × 105 GC/100 ml were reported. The presence of this marker in an urban stream is explained by the run-off from a small pasture grazed by a few heads of cattle on a private property near the Elm Park stream. Similar to the streams, the ruminant CF128 marker was present in 47% of the samples in Merrion Strand and in only one sample in Sandymount Strand.
Figure 3. Boxplot representation of the concentrations and percentage of positive samples for the human (HF183) and ruminant (CF128) marker in water samples from rivers (Tolka and Liffey), streams (Elm Park and Trimleston), and bathing waters (Sandymount and Merrion Strands). In the boxplots the lower hinge represents 25% quantile, upper hinge 75% quantile, and center line the median. The whiskers show the maximum and the lowest value. The percentage of samples above the quantification limit is indicated, only values above the quantification limit are plotted.
The human marker was detected in all river and stream samples and in nearly 90% of the bathing water samples. HF183 levels in the streams ranged from 2.3 × 103 to 1.7 × 106 GC/100 ml and from 7.6 × 102 to 4.7 × 104 GC/100 ml in river samples. The median values for rivers and streams differed by less than one order of magnitude. Variations in HF183 levels were observed in both beaches ranging from below the limit of quantification to 6.1 × 104 GC/100 ml for Merrion Strand and 2.2 × 104 GC/100 ml for Sandymount Strand. The median concentrations of the marker in marine samples were lower than those of rivers and streams.
All the rivers and streams were continuously impacted by human pollution. In contrast, the ruminant marker was less prevalent in rivers, streams, and bathing waters. Thus, anthropogenic activities are therefore likely to be the primary driver of fecal pollution at the sampling stations in the different watercourses.
Rivers and Urban Streams Have a Similar ARGs Profile
Bacteriophages may be an important reservoir of ARGs yet are often overlooked. We therefore wanted to analyze water samples for ARGs in both bacteria and bacteriophages. To validate our bacteriophage enrichment procedure, the concentrate was examined using an electron microscope. Bacteriophages with morphologies corresponding to families Myoviridae, Siphoviridae, and Podoviridae (Demuth et al., 1993; Colomer-Lluch et al., 2011) were present (Figure 4).
Figure 4. Electron micrographs showing bacteriophages isolated in concentrated water samples from Elm Park Stream (1, 4, 5, 7, 11), River Tolka (2, 8, 10), Trimleston Stream (6, 9), and raw sewage (3, 12). Images 1–3 show Myoviridae phages, images 4–6 Siphoviridae phages and images 7–9 Podoviridae phages. Scale bars represent 100 nm, whereas the scale bar in micrograph-11 represents 200 nm.
The presence of antibiotic resistance genes in the bacterial and bacteriophage fractions was assessed by selecting four ARGs, blaTEM, blaSHV, qnrS, and sul1. Bacterial ARGs were found in 73–100% of the river samples and the percentage of positive samples was even higher in the stream samples (87–100%). However, only 20–93% of the bacteriophage fractions from the rivers and 42–100% from the streams contained ARGs (Figure 5). In all cases the level of ARGs in the bacteriophage fraction was lower than in the bacterial fraction (Figure 6). The bacteriophage DNA preparations did not contain 16S rDNA.
Figure 5. Frequency of detection (%) for the antibiotic resistance genes in freshwater (River Tolka and Liffey, Elm Park, and Trimleston Streams) and bathing waters (Merrion and Sandymount Strands).
Figure 6. Boxplot representation of the ARGs concentrations in the bacterial (B) and bacteriophage (P) fractions in water samples from (A) rivers (Tolka and Liffey), (B) streams (Elm Park and Trimleston), and (C) bathing waters (Merrion and Sandymount Strands). In the boxplots the lower hinge represents 25% quantile, upper hinge 75% quantile, and center line the median. The whiskers show the maximum and the lowest value. Only values above the quantification limit are plotted.
The blaTEM and sul1 genes the most prevalent in the bacterial fraction, while blaTEM was most abundant in the bacteriophage fraction of the river and stream samples. The blaTEM median levels ranged from 1.1 × 104 to 4.0 × 104 GC/100 ml in the bacterial fraction and from 1.7 × 103 to 1.4 × 104 GC/100 ml in bacteriophage fraction. The sul1 gene was the most abundant in the bacterial fraction of the River Tolka, River Liffey, and the Trimleston stream with levels ranging from 3.2 × 104 to 7.6 × 104 GC/100 ml. Interestingly, the levels of sul1 in the bacteriophage DNA were 10- to 100-fold lower.
The median concentrations for the bacterial blaSHV and qnrS genes were in the same range (3 log10 GC/100 ml) but were one order of magnitude lower than blaTEM and sul1. For the bacteriophage fraction, the median levels for blaSHV and qnrS genes observed in rivers were in the same range (2 log10 GC/100 ml) as those found in the urban streams. However, the qnrS gene was more commonly found in the urban streams (47–64%) in comparison to the largest rivers (20–40%). Interestingly, an opposite pattern was obtained for blaSHV. Overall, the profiles of the bacterial ARGs in the urban streams and the largest rivers were similar, since no statistical differences for the majority of ARGs were reported between rivers and streams (Kruskal–Wallis, p > 0.05, Supplementary Table 3).
Spearman correlation analysis was performed to establish the relationship between the levels of ARGs and the human fecal marker. Since there was no significant difference between the ARG profiles of rivers and streams, all samples were treated as a single dataset. The blaTEM, blaSHV, and qnrS genes in the bacterial DNA fraction correlated moderately with the human fecal marker (Spearman correlation, ρ = 0.292–0.335, p < 0.05, Supplementary Table 4). However, there was no significant correlation between the human marker and ARGs in the bacteriophage fraction.
ARGs in Bathing Waters Have an Urban Profile
The rivers and streams in this study are mostly impacted by human fecal contamination and have the same ARG profile, despite their substantial difference in size. We hypothesized that the discharge of these streams and rivers into Dublin Bay would affect the nearshore marine environment.
All ARGs present in the streams and rivers were also found in the two bathing waters. The detection frequency profile of the ARGs in the rivers, streams, and strands was very similar (Figure 5). The blaTEM and sul1 genes were the most frequently detected genes, followed by blaSHV and qnrS. The levels of individual ARGs in both the bacterial and bacteriophage fractions in Merrion and Sandymount Strand did not differ significantly (Mann–Whitney, p > 0.05, Supplementary Table 5), which is not surprising, considering that these bathing waters are adjacent to each other. As was observed for the rivers and streams, in the nearshore marine environment blaTEM and sul1 were the most abundant ARGs in the bacterial fraction, and blaTEM was the most abundant in the bacteriophage fraction (Figure 6). The median level of blaTEM ranged from 2.6 × 103 to 6.3 × 103 GC/100 ml and from 3.7 × 102 to 4.8 × 103 GC/100 ml for sul1. As was the case for the rivers and streams, the levels of blaSHV and qnrS were around 0.5 log10 lower than the two most abundant genes (Figure 6).
In comparison to the freshwater samples, the median concentration for bacterial ARGs levels (blaTEM, blaSHV and sul1) decreased by significantly two and ninefold in marine bathing waters (Mann–Whitney, p = 0.0001–0.041, Supplementary Table 5 and Supplementary Figure 1), The qnrS levels did not change significantly (Mann–Whitney, p = 0.466). In contrast, there was no significant reduction in ARG levels for the phage fraction (Mann–Whitney, p = 0.093–0.653, Supplementary Table 5 and Supplementary Figure 1).
Discussion
The rivers Liffey and Tolka rise outside Dublin and flow through an agricultural area before entering the city and discharging into Dublin Bay. In contrast, the two streams included in this study are completely urban and are therefore experiencing different pollution pressures than the much larger rivers. Despite these potential differences in pollution pressures, microbial source tracking analysis showed that the main source of pollution in both rivers and streams in the urban environment is human contamination, which presumably enters the rivers and streams through for example combined sewer overflows and sewerage misconnections.
The four clinically relevant ARGs included in this study blaTEM, blaSHV, qnrS, and sul1, were present in all four waterways and occurred in both the bacterial as well as in the bacteriophage fraction. These ARGs selected confer resistance to the most common antibiotics prescribed in Ireland (HSE, 2019). Interestingly, the ARG profiles of the four waterways were very similar, both in terms of the levels of ARGs and the frequency of detection. This is consistent with the urban environment as the main fecal impactor on these waterways. In line with this, a moderate level of correlation between the bacterial blaTEM, blaSHV, and qnrS and the HF183 marker was observed, but not for sul1. In contrast, there was no correlation between the levels of the four ARGs in bacteriophage fraction and HF183. Partial correlation between ARGs and indicators of fecal contamination has been reported previously (Calero-Cáceres et al., 2017).
In the bacterial fraction, the sul1 and blaTEM genes were the most prevalent followed by blaSHV and qnrS. The blaTEM gene was the most prominent ARG in the bacteriophage fraction, followed by sul1, blaSHV, and qnrS. Resistance to β-lactam antimicrobial agents by blaTEM, blaSHV are widely distributed in aquatic ecosystems (Colomer-Lluch et al., 2011; Anand et al., 2016; Calero-Cáceres et al., 2017; Zhang et al., 2019), which might indicate that they are particularly amenable to propagation through transduction (Wang et al., 2018a). Sulfonamides are one of the oldest antimicrobial synthetic or semi-synthetic drug classes that have been also used for the treatment of animals (Dasenaki and Thomaidis, 2017) and is an authorized antibiotic for use in aquaculture in Europe (Santos and Ramos, 2018). The sul1 gene was the second most abundant and variable in concentration in the bacteriophage fraction. This high degree of variability of sul1 in rivers has also been observed in Mediterranean human-impacted rivers (Calero-Cáceres et al., 2017; Lekunberri et al., 2017) and rivers in China (Yang et al., 2018). Sulfonamide resistance is associated with mobile genetic elements, like class 1 integrons (int1) (Jiang et al., 2019), which may explain the presence in the bacteriophage fraction. Resistance to fluoroquinolones have been associated with clinical Enterobacteriaceae isolates as well as in waterborne bacteria, but with lower prevalence in rivers (Poirel et al., 2012; Marti and Balcázar, 2013; Lekunberri et al., 2017).
Aquatic environments are ideally suited for the dispersal of ARGs (Marti et al., 2014a). Rivers and streams may therefore have a lasting effect on the presence of ARGs in the marine environment into which they discharge. In addition, the discharge of wastewater treatment plants will also add to the presence of ARG in the marine environment. The ARG profiles of two bathing waters in Dublin Bay were not significantly different from those of the rivers and streams, although the concentrations and frequency of detection were lower. The blaTEM gene is also the most abundant in the marine environment, which is consistent with the few marine studies in the Mediterranean Sea and the Indian Ocean, where blaTEM gene was the most prevalent and abundant gene in the bacteriophage fraction (Calero-Cáceres and Balcázar, 2019; Blanco-Picazo et al., 2020).
Interestingly, there was a significant two to ninefold reduction in median values of the bacterial blaTEM, blaSHV, and sul1 genes, whereas the levels of qnrS did not change significantly. Freshwater or intestinal bacteria carrying ARGs are likely to die off rapidly in the marine environment, which would account for the decrease of three of the four ARGs. Resistance to fluoroquinolone in marine environments is related to the intrinsic resistance of marine Vibrionaceae and Shewanellaceae family species which possess chromosome-encoded Qnr-like proteins (Poirel et al., 2005; Cattoir et al., 2007). The decrease of bacterial ARGs in the marine environment was not observed for bacteriophage ARGs. In general, decay rates of bacteriophages in seawater are much lower than those of fecally associated bacteria and therefore persist for prolonged periods of time (Mocé-Llivina et al., 2005; Calero-Cáceres and Muniesa, 2016; Wu et al., 2020; Sala-Comorera et al., 2021a).
Once ARGs have been introduced to the phageome, they may be reintroduced into the bacterial metagenome through horizontal gene transfer mechanisms (Muniesa et al., 2013a). Prolonged persistence of bacteriophages carrying ARGs, may result in transfer of these ARGs to marine microbiota, and eventually find their way into human consumers of seafood. In addition, and perhaps more importantly, ARG carrying bacteriophages may accumulate in filter feeders (e.g., oysters) and enter the food chain in this manner. Furthermore, recreational activities taking place in or on the water may expose people to ARG carrying bacteriophages, which once ingested may transfer these ARGs to the microbiota in the intestinal tract. The EU Bathing Water Directive classifies bathing water solely on the presence of E. coli and intestinal enterococci, with a view to prevent gastrointestinal and respiratory disease (EU, 2006). Although the use of bacteriophages as additional or alternative indicators for fecal contamination has and is being discussed, the data presented here make a case for the inclusion of bacteriophages as indicators for the potential spread of antimicrobial resistance during recreational use of bathing waters.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
LS-C: conceptualization, methodology, validation, formal analysis, data curation, writing—original draft, writing—review and editing, and visualization. TN and LR: investigation and validation. AV, LC, NM, JS, and AG: investigation. GO’H and JO’S: writing—review and editing. WM: conceptualization, supervision, writing—review and editing, supervision, project administration, and funding acquisition. All authors contributed to the article and approved the submitted version.
Funding
This research (Acclimatize) was part funded by the European Regional Development Fund through the Ireland Wales Cooperation program.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.718234/full#supplementary-material
References
AFNOR. (2015). Méthodes D’analyse En Santé Animale — Pcr (Réaction De Polymérisation En Chaîne) — Partie 2?: exigences Et Recommandations Pour Le Développement Et La Validation De La Pcr En Santé Animale. France: AFNOR.
Ahmed, W., Payyappat, S., Cassidy, M., Harrison, N., and Besley, C. (2020). Sewage-Associated Marker Genes Illustrate The Impact Of Wet Weather Overflows And Dry Weather Leakage In Urban Estuarine Waters Of Sydney. Australia. Sci. Total Environ. 705:135390. doi: 10.1016/J.Scitotenv.2019.135390
Allcock, S., Young, E. H., Holmes, M., Gurdasani, D., Dougan, G., Sandhu, M. S., et al. (2017). Antimicrobial Resistance In Human Populations: challenges And Opportunities. Glob. Heal. Epidemiol. Genom. 2, 1–7. doi: 10.1017/Gheg.2017.4
Anand, T., Bera, B. C., Vaid, R. K., Barua, S., Riyesh, T., Virmani, N., et al. (2016). Abundance Of Antibiotic Resistance Genes In Environmental Bacteriophages. J. Gen. Virol. 97, 3458–3466. doi: 10.1099/Jgv.0.000639
Andersson, D. I., and Hughes, D. (2010). Antibiotic Resistance And Its Cost: is It Possible To Reverse Resistance? Nat. Rev. Microbiol. 8, 260–271. doi: 10.1038/Nrmicro2319
Balcazar, J. L. (2014). Bacteriophages As Vehicles For Antibiotic Resistance Genes In The Environment. PLoS Pathog. 10:E1004219. doi: 10.1371/Journal.Ppat.1004219
Balcázar, J. L. (2018). How Do Bacteriophages Promote Antibiotic Resistance In The Environment? Clin. Microbiol. Infect. 24, 447–449. doi: 10.1016/J.Cmi.2017.10.010
Ballesté, E., Demeter, K., Masterson, B., Timoneda, N., Sala-Comorera, L., and Meijer, W. G. (2020). Implementation And Integration Of Microbial Source Tracking In A River Watershed Monitoring Plan. Sci. Total Environ. 736:139573. doi: 10.1016/J.Scitotenv.2020.139573
Baquero, F., Coque, T. M., Martínez, J. L., Aracil-Gisbert, S., and Lanza, V. F. (2019). Gene Transmission In The One Health Microbiosphere And The Channels Of Antimicrobial Resistance. Front. Microbiol. 10:2892. doi: 10.3389/Fmicb.2019.02892
Battaglioli, E. J., Baisa, G. A., Weeks, A. E., Schroll, R. A., Hryckowian, A. J., and Welch, R. A. (2011). Isolation Of Generalized Transducing Bacteriophages For Uropathogenic Strains Of Escherichia Coli. Appl. Environ. Microbiol. 77, 6630–6635. doi: 10.1128/Aem.05307-11
Blanchard, P., Regnault, J., Schurr, F., Dubois, E., and Ribière, M. (2012). Intra-Laboratory Validation Of Chronic Bee Paralysis Virus Quantitation Using An Accredited Standardised Real-Time Quantitative Rt-Pcr Method. J. Virol. Methods 180, 26–31. doi: 10.1016/J.Jviromet.2011.12.005
Blanco-Picazo, P., Roscales, G., Toribio-Avedillo, D., Gómez-Gómez, C., Avila, C., Ballesté, E., et al. (2020). Antibiotic Resistance Genes In Phage Particles From Antarctic and Mediterranean Seawater Ecosystems. Microorganisms 8:1293. doi: 10.3390/Microorganisms8091293
Brown-Jaque, M., Calero-Cáceres, W., and Muniesa, M. (2015). Transfer Of Antibiotic-Resistance Genes Via Phage-Related Mobile Elements. Plasmid 79, 1–7. doi: 10.1016/J.Plasmid.2015.01.001
Calero-Cáceres, W., and Balcázar, J. L. (2019). Antibiotic Resistance Genes In Bacteriophages From Diverse Marine Habitats. Sci. Total Environ. 654, 452–455. doi: 10.1016/J.Scitotenv.2018.11.166
Calero-Cáceres, W., Melgarejo, A., Colomer-Lluch, M., Stoll, C., Lucena, F., Jofre, J., et al. (2014). Sludge As A Potential Important Source Of Antibiotic Resistance Genes In Both The Bacterial and Bacteriophage Fractions. Environ. Sci. Technol. 48, 7602–7611. doi: 10.1021/Es501851s
Calero-Cáceres, W., Méndez, J., Martín-Díaz, J., and Muniesa, M. (2017). The Occurrence Of Antibiotic Resistance Genes In A Mediterranean River and Their Persistence In The Riverbed Sediment. Environ. Pollut. 223, 384–394. doi: 10.1016/J.Envpol.2017.01.035
Calero-Cáceres, W., and Muniesa, M. (2016). Persistence Of Naturally Occurring Antibiotic Resistance Genes In The Bacteria and Bacteriophage Fractions Of Wastewater. Water Res. 95, 11–18. doi: 10.1016/J.Watres.2016.03.006
Cassini, A., Högberg, L. D., Plachouras, D., Quattrocchi, A., Hoxha, A., Simonsen, G. S., et al. (2019). Attributable Deaths and Disability-Adjusted Life-Years Caused By Infections With Antibiotic-Resistant Bacteria In The Eu and The European Economic Area In 2015: a Population-Level Modelling Analysis. Lancet Infect. Dis. 19, 56–66. doi: 10.1016/S1473-3099(18)30605-4
Cattoir, V., Poirel, L., Mazel, D., Soussy, C. J., and Nordmann, P. (2007). Vibrio Splendidus As The Source Of Plasmid-Mediated Qnrs-Like Quinolone Resistance Determinants [1]. Antimicrob. Agents Chemother. 51, 2650–2651. doi: 10.1128/Aac.00070-07
Chen, J., Quiles-Puchalt, N., Chiang, Y. N., Bacigalupe, R., Fillol-Salom, A., Chee, M. S. J., et al. (2018). Genome Hypermobility By Lateral Transduction. Science 362, 207–212. doi: 10.1126/Science.Aat5867
Chiang, Y. N., Penadés, J. R., and Chen, J. (2019). Genetic Transduction By Phages and Chromosomal Islands: the New and Noncanonical. PLoS Pathog. 15:E1007878. doi: 10.1371/Journal.Ppat.1007878
Colomer-Lluch, M., Jofre, J., and Muniesa, M. (2011). Antibiotic Resistance Genes In The Bacteriophage Dna Fraction Of Environmental Samples. PLoS One 6:E17549. doi: 10.1371/Journal.Pone.0017549
Dasenaki, M. E., and Thomaidis, N. S. (2017). “Meat Safety: ii Residues and Contaminants” In Lawrie’s Meat Science: eighth Edition. ed. F. Toldrá. (Netherlands: Elsevier). 553–583. doi: 10.1016/B978-0-08-100694-8.00018-2
Demuth, J., Neve, H., and Witzel, K. P. (1993). Direct Electron Microscopy Study On The Morphological Diversity Of Bacteriophage Populations In Lake Plusssee. Appl. Environ. Microbiol. 59, 3378–3384. doi: 10.1128/Aem.59.10.3378-3384.1993
EU. (2006). Directive 2006/7/Ec Of The European Parliament and Of The Council Of 15 February 2006 Concerning The Management Of Bathing Water Quality And Repealing Directive 76/160/Eec. EU: EUR-Lex.
Gourmelon, M., Caprais, M. P., Ségura, R., Le Mennec, C., Lozach, S., Piriou, J. Y., et al. (2007). Evaluation Of Two Library-Independent Microbial Source Tracking Methods To Identify Sources Of Fecal Contamination In French Estuaries. Appl. Environ. Microbiol. 73, 4857–4866. doi: 10.1128/Aem.03003-06
Gunathilaka, G. U., Tahlan, V., Mafiz, A. I., Polur, M., and Zhang, Y. (2017). Phages In Urban Wastewater Have The Potential To Disseminate Antibiotic Resistance. Int. J. Antimicrob. Agents 50, 678–683. doi: 10.1016/J.Ijantimicag.2017.08.013
Guo, J., Li, J., Chen, H., Bond, P. L., and Yuan, Z. (2017). Metagenomic Analysis Reveals Wastewater Treatment Plants As Hotspots Of Antibiotic Resistance Genes and Mobile Genetic Elements. Water Res. 123, 468–478. doi: 10.1016/J.Watres.2017.07.002
Hernando-Amado, S., Coque, T. M., Baquero, F., and Martínez, J. L. (2019). Defining And Combating Antibiotic Resistance From One Health And Global Health Perspectives. Nat. Microbiol. 4, 1432–1442. doi: 10.1038/S41564-019-0503-9
Hernando-Amado, S., Coque, T. M., Baquero, F., and Martínez, J. L. (2020). Antibiotic Resistance: moving From Individual Health Norms To Social Norms In One Health And Global Health. Front. Microbiol. 11:1914. doi: 10.3389/Fmicb.2020.01914
Hinojosa, J., Green, J., Estrada, F., Herrera, J., Mata, T., Phan, D., et al. (2020). Determining The Primary Sources Of Fecal Pollution Using Microbial Source Tracking Assays Combined With Land-Use Information In The Edwards Aquifer. Water Res. 184:116211. doi: 10.1016/J.Watres.2020.116211
Huddleston, J. R. (2014). Horizontal Gene Transfer In The Human Gastrointestinal Tract: potential Spread Of Antibiotic Resistance Genes. Infect. Drug Resist. 7, 167–176. doi: 10.2147/Idr.S48820
ISO. (2000). Iso 7899-2:2000 - Water Quality – Detection And Enumeration Of Intestinal Enterococci – Part 2: membrane Filtration Method. Switzerland: ISO.
ISO. (2001). Iso 16649-1:2001 - Microbiology Of Food And Animal Feeding Stuffs — Horizontal Method For The Enumeration Of Beta-Glucuronidase-Positive Escherichia Coli — Part 1: colony-Count Technique At 44 Degrees C Using Membranes And 5-Bromo-4-Chloro-3-Indolyl beta-D-glucuronide. Switzerland: ISO.
Jiang, H., Cheng, H., Liang, Y., Yu, S., Yu, T., Fang, J., et al. (2019). Diverse Mobile Genetic Elements And Conjugal Transferability Of Sulfonamide Resistance Genes (Sul1, Sul2, And Sul3) In Escherichia Coli Isolates From Penaeus Vannamei And Pork From Large Markets In Zhejiang. China. Front. Microbiol. 10:1787. doi: 10.3389/Fmicb.2019.01787
Jit, M., Ng, D. H. L., Luangasanatip, N., Sandmann, F., Atkins, K. E., Robotham, J. V., et al. (2020). Quantifying The Economic Cost Of Antibiotic Resistance And The Impact Of Related Interventions: rapid Methodological Review, Conceptual Framework And Recommendations For Future Studies. BMC Med. 18:38. doi: 10.1186/S12916-020-1507-2
Kay, D., Crowther, J., Stapleton, C. M., Wyer, M. D., Fewtrell, L., Edwards, A., et al. (2008). Faecal Indicator Organism Concentrations In Sewage And Treated Effluents. Water Res. 42, 442–454. doi: 10.1016/J.Watres.2007.07.036
Lekunberri, I., Villagrasa, M., Balcázar, J. L., and Borrego, C. M. (2017). Contribution Of Bacteriophage And Plasmid Dna To The Mobilization Of Antibiotic Resistance Genes In A River Receiving Treated Wastewater Discharges. Sci. Total Environ. 60, 206–209. doi: 10.1016/J.Scitotenv.2017.05.174
Leonard, A. F. C., Zhang, L., Balfour, A. J., Garside, R., Hawkey, P. M., Murray, A. K., et al. (2018). Exposure To And Colonisation By Antibiotic-Resistant E. Coli In Uk Coastal Water Users: environmental Surveillance, Exposure Assessment, And Epidemiological Study (Beach Bum Survey). Environ. Int. 114, 326–333. doi: 10.1016/J.Envint.2017.11.003
Maganha De, A. K. A. C., Borrego, C. M., and Balcázar, J. L. (2021). Antimicrobial Resistance And Bacteriophages: an Overlooked Intersection In Water Disinfection. Trends Microbiol. 29, 517–527. doi: 10.1016/J.Tim.2020.12.011
Marti, E., and Balcázar, J. L. (2013). Comprehensive Analytical Chemistry. Netherlands: elsevier B.V. 671–684. doi: 10.1016/B978-0-444-62657-8.00019-7
Marti, E., Variatza, E., and Balcazar, J. L. (2014a). The Role Of Aquatic Ecosystems As Reservoirs Of Antibiotic Resistance. Trends Microbiol. 22, 36–41. doi: 10.1016/J.Tim.2013.11.001
Marti, E., Variatza, E., and Balcázar, J. L. (2014b). Bacteriophages As A Reservoir Of Extended-Spectrum B -Lactamase And Fluoroquinolone Resistance Genes In The Environment. Clin. Microbiol. Infect. 20, O456–O459. doi: 10.1111/1469-0691.12446
Mocé-Llivina, L., Lucena, F., and Jofre, J. (2005). Enteroviruses And Bacteriophages In Bathing Waters. Appl. Environ. Microbiol. 71, 6838–6844. doi: 10.1128/Aem.71.11.6838-6844.2005
Molina, M., Hunter, S., Cyterski, M., Peed, L. A., Kelty, C. A., Sivaganesan, M., et al. (2014). Factors Affecting The Presence Of Human-Associated And Fecal Indicator Real-Time Quantitative Pcr Genetic Markers In Urban-Impacted Recreational Beaches. Water Res. 64, 196–208. doi: 10.1016/J.Watres.2014.06.036
Muniesa, M., Colomer-Lluch, M., and Jofre, J. (2013a). Could Bacteriophages Transfer Antibiotic Resistance Genes From Environmental Bacteria To Human-Body Associated Bacterial Populations? Mob. Genet. Elements 3:E25847. doi: 10.4161/Mge.25847
Muniesa, M., Colomer-Lluch, M., and Jofre, J. (2013b). Potential Impact Of Environmental Bacteriophages In Spreading Antibiotic Resistance Genes. Future Microbiol. 8, 739–751. doi: 10.2217/Fmb.13.32
Nguyen, A. Q., Vu, H. P., Nguyen, L. N., Wang, Q., Djordjevic, S. P., Donner, E., et al. (2021). Monitoring Antibiotic Resistance Genes In Wastewater Treatment: current Strategies And Future Challenges. Sci. Total Environ. 783:146964. doi: 10.1016/J.Scitotenv.2021.146964
Poirel, L., Cattoir, V., and Nordmann, P. (2012). Plasmid-Mediated Quinolone Resistance; Interactions Between Human, Animal, and Environmental Ecologies. Front. Microbiol. 3:24. doi: 10.3389/Fmicb.2012.00024
Poirel, L., Rodriguez-Martinez, J. M., Mammeri, H., Liard, A., and Nordmann, P. (2005). Origin Of Plasmid-Mediated Quinolone Resistance Determinant Qnra. Antimicrob. Agents Chemother. 49, 3523–3525. doi: 10.1128/Aac.49.8.3523-3525.2005
Quirós, P., Colomer-Lluch, M., Martínez-Castillo, A., Miró, E., Argente, M., Jofre, J., et al. (2014). Antibiotic Resistance Genes In The Bacteriophage Dna Fraction Of Human Fecal Samples. Antimicrob. Agents Chemother. 58, 606–609. doi: 10.1128/Aac.01684-13
Reynolds, L. J., Martin, N. A., Sala-Comorera, L., Callanan, K., Doyle, P., O’leary, C., et al. (2021). Identifying Sources Of Faecal Contamination In A Small Urban Stream Catchment: a Multiparametric Approach. Front. Microbiol. 12:1580. doi: 10.3389/Fmicb.2021.661954
Reynolds, L. J., Sala-Comorera, L., Martin, N. A., Nolan, T. M., Stephens, J. H., Gitto, A., et al. (2020). Correlation Between Antimicrobial Resistance And Faecal Contamination In Small Urban Streams And Bathing Waters. Sci. Total Environ. 739:140242. doi: 10.1016/J.Scitotenv.2020.140242
Rizzo, L., Manaia, C., Merlin, C., Schwartz, T., Dagot, C., Ploy, M. C., et al. (2013). Urban Wastewater Treatment Plants As Hotspots For Antibiotic Resistant Bacteria And Genes Spread Into The Environment: a Review. Sci. Total Environ. 447, 345–360. doi: 10.1016/J.Scitotenv.2013.01.032
Rodriguez-Mozaz, S., Chamorro, S., Marti, E., Huerta, B., Gros, M., Sànchez-Melsió, A., et al. (2015). Occurrence Of Antibiotics And Antibiotic Resistance Genes In Hospital And Urban Wastewaters And Their Impact On The Receiving River. Water Res. 69, 234–242. doi: 10.1016/J.Watres.2014.11.021
Roope, L. S. J., Smith, R. D., Pouwels, K. B., Buchanan, J., Abel, L., Eibich, P., et al. (2019). The Challenge Of Antimicrobial Resistance: what Economics Can Contribute. Science 364:eaau4679. doi: 10.1126/Science.Aau4679
Ross, J., and Topp, E. (2015). Abundance Of Antibiotic Resistance Genes In Bacteriophage Following Soil Fertilization With Dairy Manure Or Municipal Biosolids, And Evidence For Potential Transduction. Appl. Environ. Microbiol. 81, 7905–7913. doi: 10.1128/Aem.02363-15
Rutledge, R. G. (2003). Mathematics Of Quantitative Kinetic Pcr And The Application Of Standard Curves. Nucleic Acids Res. 31:e93. doi: 10.1093/Nar/Gng093
Sala-Comorera, L., Reynolds, L. J., Martin, N. A., O’sullivan, J. J., Meijer, W. G., and Fletcher, N. F. (2021a). Decay Of Infectious Sars-Cov-2 And Surrogates In Aquatic Environments. Water Res. 201:117090. doi: 10.1016/J.Watres.2021.117090
Sala-Comorera, L., Reynolds, L. J., Martin, N. A., Pascual-Benito, M., Stephens, J. H., Nolan, T. M., et al. (2021b). CrAssphage As Human Molecular Marker To Evaluate Temporal And Spatial Variability In Faecal Contamination Of Urban Marine Bathing Waters. Sci. Total Environ. 789:147828. doi: 10.1016/J.Scitotenv.2021.147828
Santos, L., and Ramos, F. (2018). Antimicrobial Resistance In Aquaculture: current Knowledge and Alternatives To Tackle The Problem. Int. J. Antimicrob. Agents 52, 135–143. doi: 10.1016/J.Ijantimicag.2018.03.010
Stewardson, A. J., Allignol, A., Beyersmann, J., Graves, N., Schumacher, M., Meyer, R., et al. (2016). The Health And Economic Burden Of Bloodstream Infections Caused By Antimicrobial-Susceptible And Non-Susceptible Enterobacteriaceae And Staphylococcus Aureus In European Hospitals, 2010 And 2011: a Multicentre Retrospective Cohort Study. Eurosurveillance 21:30319. doi: 10.2807/1560-7917.Es.2016.21.33.30319
Unc, A., and Goss, M. J. (2004). Transport Of Bacteria From Manure And Protection Of Water Resources. Appl. Soil Ecol. 25, 1–18. doi: 10.1016/J.Apsoil.2003.08.007
Van Bruggen, A. H. C., Goss, E. M., Havelaar, A., Van Diepeningen, A. D., Finckh, M. R., and Morris, J. G. (2019). One Health - Cycling Of Diverse Microbial Communities As A Connecting Force For Soil, Plant, Animal, Human And Ecosystem Health. Sci. Total Environ. 664, 927–937. doi: 10.1016/J.Scitotenv.2019.02.091
Wang, M., Liu, P., Zhou, Q., Tao, W., Sun, Y., and Zeng, Z. (2018a). Estimating The Contribution Of Bacteriophage To The Dissemination Of Antibiotic Resistance Genes In Pig Feces. Environ. Pollut. 238, 291–298. doi: 10.1016/J.Envpol.2018.03.024
Wang, M., Xiong, W., Liu, P., Xie, X., Zeng, J., Sun, Y., et al. (2018b). Metagenomic Insights Into The Contribution Of Phages To Antibiotic Resistance In Water Samples Related To Swine Feedlot Wastewater Treatment. Front. Microbiol. 9:2474. doi: 10.3389/Fmicb.2018.02474
WHO (2019). No Time to Wait: securing The Future From Drug-Resistant Infections. Switzerland: world Health Organization.
Wu, Z., Greaves, J., Arp, L., Stone, D., and Bibby, K. (2020). Comparative Fate Of Crassphage With Culturable And Molecular Fecal Pollution Indicators During Activated Sludge Wastewater Treatment. Environ. Int. 136:105452. doi: 10.1016/J.Envint.2019.105452
Yang, Y., Shi, W., Lu, S. Y., Liu, J., Liang, H., Yang, Y., et al. (2018). Prevalence Of Antibiotic Resistance Genes In Bacteriophage Dna Fraction From Funan River Water In Sichuan. China. Sci. Total Environ. 626, 835–841. doi: 10.1016/J.Scitotenv.2018.01.148
Yang, Y., Xing, S., Chen, Y., Wu, R., Wu, Y., Wang, Y., et al. (2021). Profiles Of Bacteria/Phage-Comediated Args In Pig Farm Wastewater Treatment Plants In China: association With Mobile Genetic Elements, Bacterial Communities And Environmental Factors. J. Hazard. Mater. 404:124149. doi: 10.1016/J.Jhazmat.2020.124149
Zhang, A., Call, D. R., Besser, T. E., Liu, J., Jones, L., Wang, H., et al. (2019). B -Lactam Resistance Genes In Bacteriophage And Bacterial Dna From Wastewater, River Water, And Irrigation Water In Washington State. Water Res. 161, 335–340. doi: 10.1016/J.Watres.2019.06.026
Zhou, P., Yang, X., Lou, Wang, X. G., Hu, B., Zhang, L., et al. (2020). A Pneumonia Outbreak Associated With A New Coronavirus Of Probable Bat Origin. Nature 579, 270–273. doi: 10.1038/S41586-020-2012-7
Keywords: antibiotic resistance genes, microbial source tracking, bacteriophages, fecal pollution, rivers, urban streams, bathing waters
Citation: Sala-Comorera L, Nolan TM, Reynolds LJ, Venkatesh A, Cheung L, Martin NA, Stephens JH, Gitto A, O’Hare GMP, O’Sullivan JJ and Meijer WG (2021) Bacterial and Bacteriophage Antibiotic Resistance in Marine Bathing Waters in Relation to Rivers and Urban Streams. Front. Microbiol. 12:718234. doi: 10.3389/fmicb.2021.718234
Received: 31 May 2021; Accepted: 28 June 2021;
Published: 26 July 2021.
Edited by:
Lorena Rodriguez-Rubio, University of Barcelona, SpainReviewed by:
William Calero-Cáceres, Technical University of Ambato, EcuadorMaite Muniesa, University of Barcelona, Spain
Copyright © 2021 Sala-Comorera, Nolan, Reynolds, Venkatesh, Cheung, Martin, Stephens, Gitto, O’Hare, O’Sullivan and Meijer. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wim G. Meijer, d2ltLm1laWplckB1Y2QuaWU=