- 1Department of Biology and Ecology Center, Utah State University, Logan, UT, United States
- 2Cary Institute of Ecosystem Studies, Millbrook, NY, United States
- 3Department of Plant and Wildlife Sciences, Brigham Young University, Provo, UT, United States
- 4Department of Biology, Southern Utah University, Cedar City, UT, United States
Bacteria in stream biofilms contribute to stream biogeochemical processes and are potentially sensitive to the substantial levels of pollution entering urban streams. To examine the effects of contaminants on stream biofilm bacteria in situ, we exposed growing biofilms to experimental additions of nutrients [nitrogen (N), phosphorus (P), and iron (Fe)], pharmaceuticals (caffeine and diphenhydramine), nutrients plus pharmaceuticals, or no contaminants using contaminant exposure substrates (CES) in three catchments in northern Utah. We performed our study at montane and urban sites to examine the influence of existing pollution on biofilm response. We identified bacterial core communities (core) for each contaminant treatment at each land-use type (e.g., nutrient addition montane bacterial core, nutrient addition urban bacterial core, pharmaceutical addition montane bacterial core) by selecting all taxa found in at least 75% of the samples belonging to each specific grouping. Montane and urban land-use distinguished bacterial cores, while nutrients and pharmaceuticals had subtle, but nonetheless distinct effects. Nutrients enhanced the dominance of already abundant copiotrophs [i.e., Pseudomonadaceae (Gammaproteobacteria) and Comamonadaceae (Betaproteobacteria)] within bacterial cores at montane and urban sites. In contrast, pharmaceuticals fostered species-rich bacterial cores containing unique contaminant-degrading taxa within Pseudomonadaceae and Anaerolineaceae (Chloroflexi). Surprisingly, even at urban sites containing ambient pharmaceutical pollution, pharmaceutical additions increased bacterial core richness, specifically within DR-16 (Betaproteobacteria), WCHB1-32 (Bacteroidetes), and Leptotrichiaceae (Fusobacteria). Nutrients exerted greater selective force than pharmaceuticals in nutrient plus pharmaceutical addition treatments, creating bacterial cores more closely resembling those under nutrient rather than pharmaceutical addition, and promoting unique Oscillatoriales (Cyanobacteria) taxa in urban streams. Our results show that additions of N, P, and Fe intensified the dominance of already abundant copiotrophs, while additions of caffeine and diphenhydramine enabled unique taxa associated with contaminant degradation to participate in bacterial cores. Further, biofilm bacteria at urban sites remained sensitive to pharmaceuticals commonly present in waters, suggesting a dynamic interplay among pharmaceutical pollution, bacterial diversity, and contaminant degradation.
Introduction
Biofilms house the primary form of bacterial life in freshwater streams (Costerton et al., 1978; Geesey et al., 1978; Lock et al., 1984; Battin et al., 2016). Biofilms form from single or co-aggregated microbial cells attaching to stream benthic surfaces, forming chain-like microcolonies that coalesce into larger and more complex communities of bacteria, archaea, algae, fungi, protozoa, and viruses embedded in a polysaccharide extracellular matrix (Battin et al., 2007, 2016; Besemer et al., 2012; Besemer, 2015). Upstream biological propagules of dispersing bacteria enter the water column from soils and tributaries within the surrounding catchment and contribute to burgeoning biofilms. Typically, stream biofilms are dominated by taxa from three main phyla: Proteobacteria, Bacteroidetes, and Cyanobacteria (Besemer, 2015) and natural physiochemical water conditions such as pH, water temperature, and dissolved organic carbon concentrations influence biofilm bacterial composition (Besemer, 2015). For example, stream bacterial communities were strongly correlated with stream water pH across the 3000-ha Hubbard Brook watershed in New Hampshire, United States, where an increase in stream water pH was associated with an increase in the relative abundance of Proteobacteria and decrease in Acidobacteria (Fierer et al., 2007b). In addition to natural river conditions, contaminants associated with human activities may dramatically alter biofilm bacterial biomass and composition (Tank and Dodds, 2003; Rosi-Marshall et al., 2013; Wang et al., 2018).
Nutrients entering streams and rivers through wastewater, fertilizers, and atmospheric deposition (Carey et al., 2013) are major freshwater contaminants with potential to alter bacterial communities. Human activities have dramatically enhanced global nitrogen (N) and phosphorus (P) cycles (Falkowski et al., 2000), with global riverine exports of dissolved inorganic N and P estimated to have increased by about 30% between 1970 and 2000 (Seitzinger et al., 2010). Nutrient pollution is recognized as a leading contributor to water quality impairment (Withers et al., 2014) and can have important consequences for stream bacteria. The growth of aquatic heterotrophs and oxygenic phototrophs is often limited by the availability of N and P, and to a lesser extent iron (Fe) required for photosynthesis (Van Horn et al., 2011; Larson et al., 2018). Further, nutrients may alter bacterial composition as taxa differ in their ability to thrive under various nutrient conditions (Fierer et al., 2007a). Among Proteobacteria, one of the major constituents of stream biofilms, Betaproteobacteria and Gammaproteobacteria contain many opportunistic or copiotroph species that thrive in nutrient-rich conditions (Shi et al., 1999; Fierer et al., 2007a), while Alphaproteobacteria may thrive in oligotrophic waters (Pinhassi and Berman, 2003). In heterotrophic biofilms grown in stream mesocosms, a 10-fold increase in nutrient concentrations led to a 20-fold rise in biomass and enhanced the abundance of copiotrophic heterotrophs such as Betaproteobacteria and Gammaproteobacteria taxa (Van Horn et al., 2011). Urban nutrient run-off in the cities surrounding the Jialing River, China, has dramatically altered biofilms with elevated levels of P, nitrate, and Fe explaining 63.0% of the variation within bacterial communities (Wang et al., 2018).
Like nutrients, pharmaceuticals are ubiquitous in human impacted surface waters (Kolpin et al., 2002; Pal et al., 2010) and can alter biofilm bacterial biomass and composition (Lawrence et al., 2005; Rosi-Marshall et al., 2013; Rosi et al., 2018). Pharmaceutical pollution is also a global issue: pharmaceuticals or their transformation products are at least detectable in many aquatic ecosystems of seventy-one countries on all continents (aus der Beek et al., 2016). Although typically detected at trace concentrations in freshwaters (ng L–1 to μg L–1; Kolpin et al., 2002; Pal et al., 2010), many pharmaceuticals are designed to be biologically active at even low doses and consequently, are expected to affect non-target organisms, such as bacteria. Caffeine, a stimulant found in a variety of beverages (e.g., coffee, tea, soft drinks, and energy drinks) and medications and diphenhydramine, an antihistamine, are both particularly prevalent pharmaceutical pollutants in freshwaters (Kolpin et al., 2002) with both altering biofilm bacteria. Exposure to high experimental caffeine and diphenhydramine additions decreased respiration of heterotrophic biofilms by 53 and 63%, respectively (Rosi-Marshall et al., 2013). The same study also found that diphenhydramine reduced the relative abundance of Flavobacterium, a common constituent of stream biofilms (Besemer et al., 2012), but increased the abundance of durable Pseudomonas species (Gammaproteobacteria), a genus known for its ability to degrade high molecular weight organic compounds (Dastidar et al., 1976; Palleroni, 2010; Rosi-Marshall et al., 2013), suggesting that pharmaceuticals may have negative effects on sensitive taxa but positive effects on taxa capable of degrading recalcitrant compounds. Indeed, caffeine, which may serve as the sole carbon (C) and N source for certain bacteria, increased bacterial biomass and biofilm thickness by over 30% in biofilms from the Saskatchewan River, Canada (Lawrence et al., 2005). Stream bacterial communities may even respond to pharmaceutical pollution by demonstrating enhanced contaminant-degradation potential; the abundance of xenobiotic metabolic genes within bacterial communities increased with the highly-urbanized Jialing River, China (Wang et al., 2018). Given the ubiquitous nature and potential influence of pharmaceuticals on stream biofilms, research to understand the effects of pharmaceuticals such as caffeine and diphenhydramine can enhance our understanding of anthropogenic effects on stream ecosystems.
Contaminant exposure substrates (CES; Rosi-Marshall et al., 2013; Costello et al., 2016) can be used to examine the effects of nutrients, pharmaceuticals, and other contaminants on stream biofilm bacteria in situ. CES employ a similar design as nutrient diffusing substrates (NDS) which are widely used in stream ecology to examine nutrient limitation (Tank et al., 2007; Capps et al., 2011). These devices consist of a small container filled with unamended (control) or contaminant-amended agar (contaminant treatment) capped with a porous surface such as a fritted glass disk or cellulose sponge (Tank et al., 2007; Capps et al., 2011). CES are incubated in the stream for 2–3 weeks and biofilms which colonize the porous surfaces are exposed to contaminants diffusing out of the agar. Contaminant effects on biofilm biomass, activity, or community composition can be assessed by comparing biofilms on control and contaminant treatment CES. CES provide ecologically-relevant results as studies are conducted in situ. Further, the method is cost-effective and easily replicated, enabling the inclusion of multiple contaminant treatments and study sites. Importantly, contaminant diffusion rates typically decline throughout the in-stream incubation and may be influenced multiple factors such as temperature, stream flow, and water column contaminant concentrations (Shaw et al., 2015). In addition, the composition of biofilms grown on CES may differ from natural biofilms, as CES biofilms develop over only 2–3 weeks on porous surfaces which are structurally homogenous in comparison to natural stream surfaces (Hoellein et al., 2010). Thus, CES may be less appropriate for determining specific ecotoxicology thresholds but can play an important role in detecting contaminant effects on stream biofilm growth, activity, and composition.
To more fully understand the impact of human activities on stream biofilm bacterial community composition, we used CES to assess in situ responses of stream biofilm bacteria to nutrient (N, P, Fe), pharmaceutical (caffeine, diphenhydramine) and nutrient plus pharmaceutical additions in montane and urban streams in three catchments in northern Utah, United States. We included N, P, and Fe as these nutrients may limit biofilm growth and activity at their natural levels (Van Horn et al., 2011; Larson et al., 2018) and human activities have dramatically increased N and P levels in streams worldwide (Seitzinger et al., 2010). We included caffeine and diphenhydramine in our pharmaceutical addition treatment as both are widespread contaminants in human-impacted streams (Kolpin et al., 2002) and have previously been shown to influence biofilm growth, function, and composition (Lawrence et al., 2005; Rosi-Marshall et al., 2013). We incorporated a nutrient plus pharmaceutical addition treatment to examine the combined effects of the two contaminant classes which may enter streams through similar sources (i.e., wastewater effluent, livestock) and thus co-occur (Rosi-Marshall and Royer, 2012). We identified and compared bacterial core communities (core) for each contaminant treatment by land-use combination, where bacterial core taxa were defined as all taxa that occurred in at least 75% of the samples in a specific grouping. We measured ambient nutrient and pharmaceutical concentrations to examine the potential for existing pollution to influence bacterial responses to experimental nutrient and pharmaceutical additions. We hypothesized that nutrients would select for bacterial cores dominated by copiotrophic taxa, while pharmaceuticals would select for cores composed of contaminant-degrading taxa. Further, we hypothesized that existing nutrient and pharmaceutical pollution in urban streams would decrease the sensitivity of biofilm bacteria to experimental additions of these pollutants.
Materials and Methods
Study Areas
We conducted nutrient and pharmaceutical addition bioassays at sites located along a gradient of montane to urban land-use in three catchments, Logan River, Red Butte Creek, and Middle Provo River, located along the Wasatch Range in northern Utah, United States. These catchments contain different mixtures of upland forest/sage-steppe, urban, sub-urban, and agricultural land-use. The Logan River originates in the Bear River Range, passes through a series of three hydroelectric dams located near the base of the mountains, and then flows through a mixture of urban and agricultural land-uses. Red Butte Creek flows from the foothills of the Wasatch Range where it passes through a water supply reservoir and then makes a rapid transition into the built-out, urban environment of the University of Utah campus and Salt Lake City. Located between the Jordanelle and Deer Creek Reservoirs, the Middle Provo River flows through the Heber Valley, an area undergoing a rapid change from rural to urban land-use. Our research was conducted during the summer of 2015.
Stream Biofilm Contaminant Exposure Substrates
We examined the effects of nutrients and pharmaceuticals on stream biofilm bacterial community composition using CES (photograph provided in Supplementary Figure 1; Rosi-Marshall et al., 2013; Costello et al., 2016). To construct CES, we filled 30-mL plastic cups (Polycon, Madan Plastics) with agar amended with one of four contaminant treatments: no additions (control), 0.5M N + 0.5M P + 0.5 mM Fe [nutrient (Nu) addition; 7000 mg N L–1 + 15500 mg P L–1 + 28 mg Fe L–1], 2.5 mM caffeine + 2.5 mM diphenhydramine [pharmaceutical (Ph) addition; 490 mg caffeine L–1 + 640 mg diphenhydramine L–1], and 0.5M N + 0.5M phosphorus + 0.5 mM Fe + 2.5 mM caffeine + 2.5 mM diphenhydramine [combined nutrient and pharmaceutical (NuPh) addition]. Ethylenediaminetetraacetic acid was added in all treatments with Fe. We chose these contaminant treatment concentrations to match N and P concentrations typically used in NDS studies (Tank et al., 2007; Capps et al., 2011). Similarly, previous CES studies have employed caffeine and diphenhydramine concentrations between 2.5 mM (Shaw et al., 2015) and 13–15 mM (Rosi-Marshall et al., 2013). We capped the agar with a porous glass disk and incubated three replicates of each contaminant treatment in situ at montane and urban sites for 18–26 days in the Logan River, Red Butte Creek, and Middle Provo River catchments (4 contaminant treatments × 2 land-use types × 3 catchments × 3 replicates = 72). During the in-stream incubation, biofilms which colonized the porous disks capping the CES were exposed to contaminants diffusing out of the agar, which likely declined over time (Costello et al., 2016). In a laboratory study of NDS containing 0.5M N or 0.5M P, initial N and P release rates were near 1 mg h–1 and declined log-linearly to 0.1 mg h–1 after 21 days (Rugenski et al., 2008). Diphenhydramine and caffeine release rates likely exceeded environmental concentrations by approximately four orders of magnitude. A laboratory experiment of CES containing pharmaceuticals at 2.5 mM concentrations measured average release rates of 49 ng diphenhydramine min–1 cm–2 and 41 ng caffeine min–1 cm–2 over 21 days (Shaw et al., 2015). In comparison, a typical wastewater-impacted stream with diphenhydramine concentrations of 100 ng L–1 would expose biofilms to 0.002 ng diphenhydramine min–1 cm–2 (Costello et al., 2016).
At the end of the in-stream incubation, we collected the biofilm-colonized disks in sterile plastic bags (Whirl-Pak®) and transported them on ice to the laboratory where they were stored at 4°C. Due to the loss of several CES that occurred during the in-stream incubations, only two replicates of the control and nutrient plus pharmaceutical addition treatments were recovered at the urban site in the Logan River catchment and all bacterial community inferences were based on 70 samples.
Effects on Bacterial Community Composition
We performed a 16S rRNA gene metabarcoding analysis to examine bacterial community composition in montane and urban stream biofilms exposed to contaminant treatments. We extracted genomic DNA from biofilms using the DNeasy Powersoil Kit (Qiagen, Inc., Germantown, MD, United States). We then targeted the V4 hypervariable region of the 16S rRNA gene using primers 16Sf (5′-GTGCCAGCMGCCGCGGTAA-3′) and 16Sr (5′-GGACTACHVGGGTWTCTAAT-3′). The primers contained a series of repeating 8-bp barcodes, a forward or reverse Illumina primer, linker region, and primer pad to facilitate a dual-indexed Illumina sequencing approach (Caporaso et al., 2012; Kozich et al., 2013; Aanderud et al., 2018). The InvitrogenTM AccuPrimeTM Pfx SuperMix was used for the generation of 16S rRNA gene amplicons. The thermocycler settings were: an initial denaturation step of 94°C for 3 min, followed by 35 cycles of denaturation at 94°C for 45 s, annealing at 55°C for 60 s, and elongation at 72°C for 90 s. A final elongation step was set for 72°C for 10 min and samples were then held at 4°C. Samples were then normalized using the SequalPrepTM Normalization Plate (96) Kit (Invitrogen, Carlsbad, CA, United States) and quantified on a Qubit® 2.0 Fluorometer (Thermo Fisher Scientific, Waltham, MA, United States). The samples were sequenced on an Illumina HiSeq 2500 platform (2 × 250; Illumina Biotechnology, San Diego, CA, United States) at the Brigham Young University DNA Sequencing Center1. Primer dimers were removed via Pippin Prep (Sage Science, Beverly, MA, United States). Illumina sequence reads were analyzed in mothur (v. 1.35.1; Schloss et al., 2009) following the MiSeq SOP protocol (Kozich et al., 2013). We removed chimeras with UCHIME (Edgar et al., 2011) and eliminated chloroplast, mitochondrial, archaeal, and eukaryotic gene sequences based on reference sequences from the Ribosomal Database Project (Cole et al., 2009). We then aligned sequences against the SILVA 128 database (silva.nr_v128.align; Pruesse et al., 2007) with the SEED aligner and created operational taxonomic units (OTUs) based on uncorrected pairwise distances using a minimum coverage of 99% and minimum pairwise sequence similarity of 97%. We used the ‘phyloseq’ package (McMurdie and Holmes, 2013) in R (R Core Team, 2017) to combine taxonomy files and OTU tables generated in mothur. Samples contained a total of 1,856,706 sequences and 16,438 unique OTUs with a mean sequencing coverage of 98 ± 0.1%. We rarefied data to the smallest sample size (27,704 sequences, 18,984 unique OTUs after rarefication) and then assigned OTUs from the total bacterial community into two abundance categories: abundant (≥0.1% relative abundance) and rare (<0.1% relative abundance; Pedrós-Alió, 2012; Aanderud et al., 2015). The resulting rarefied dataset was used to perform all of the following analyses.
Nutrient and Pharmaceutical Effects on Biofilm Bacterial Core Communities in Montane and Urban Streams
We identified bacterial core communities (cores) for each contaminant treatment by land-use combination (i.e., control montane, control urban, nutrient montane, etc.) by selecting all OTUs that occurred in at least 75% of the samples in a specific grouping using the ‘microbiome’ package (Lahti et al., 2017) in R. We defined bacterial cores using a 75% persistence cutoff to identify bacterial taxa that consistently responded to contaminant treatments in montane and urban biofilms and thus characterize common bacterial responses to contaminant additions at each type of land-use. Our 75% persistence cutoff also ensured sufficient number of taxa above the cutoff threshold in all bacterial cores (>200 OTUs). We visualized the composition of bacterial cores with Krona plots (Ondov et al., 2011) generated in R. Krona plots displayed the percentages of bacterial core taxa, calculated as a taxon’s total sequence count divided by the total sequence count of all taxa across all samples contributing to a bacterial core. We employed the same method to calculate and display the contribution of abundant compared with rare OTUs.
We examined the effects of land-use, nutrient addition, and pharmaceutical addition on the composition of bacterial cores. We visualized community differences among cores by performing a principal coordinate analysis using a Bray–Curtis abundance-based distance matrix from the relative recovery of OTUs with the vegdist function in the ‘vegan’ package (Oksanen et al., 2014). We then performed a PERMANOVA analysis (Anderson, 2001) using the adonis function to test for the effects of land-use, nutrient addition, and pharmaceutical addition with catchment nested within land-use, nutrient addition, and pharmaceutical addition.
Next, we focused on the effects of nutrients and pharmaceuticals. We calculated bacterial core diversity as Shannon’s diversity index, and richness as the number of OTUs in each bacterial core. We performed separate two-way ANOVAs for bacterial cores at montane and urban sites, where nutrient addition and pharmaceutical addition were included as fixed effects with catchment nested within nutrient addition and pharmaceutical addition. Treatment differences were compared using Tukey’s post hoc tests. Last, to visualize the distinctiveness of bacterial cores exposed to contaminant treatments at each type of land-use, we generated area-proportional Euler diagrams based on counts of the number of shared and unique taxa in bacterial cores of each contaminant treatment at each type of land-use using the ‘eulerr’ package (Larsson, 2019). We visualized the composition of unique taxa in bacterial cores with Krona plots created using the methods described previously.
Effects on Biofilm Biomass
We characterized the effects of nutrients and pharmaceuticals on stream biofilm biomass to assess contaminant effects on biofilm growth and facilitate comparisons to existing CES studies where biomass is frequently measured. We measured ash-free dry mass (AFDM) and chlorophyll a concentrations of biofilms grown on CES to estimate total organic biomass and photoautotrophic biomass, respectively. One-half of each fritted glass disk was used for chlorophyll determination and the other for AFDM. The disk area was estimated with the Image J software (Schneider et al., 2012). Chlorophyll was extracted in hot ethanol, refrigerated overnight, and measured on a Turner Designs TD-700 Fluorometer with pheophytin correction after acidification (Hauer and Lamberti, 2011). AFDM samples were dried for 2 days at 60°C, combusted at 450°C for 2 h, and the mass lost during combustion used to calculate AFDM (APHA, 1998). We examined the effects of nutrient addition and pharmaceutical addition on chlorophyll and AFDM in montane and urban streams using two-way ANOVAs as described above for core comparisons.
Ambient Nutrient and Pharmaceutical Concentrations
We measured ambient nutrient concentrations to examine the potential for existing nutrient pollution to influence biofilm bacteria. We collected grab samples for total nitrogen (TN) and total phosphorus (TP) analyses, samples of stream water filtered through pre-combusted Whatman GF/F filters for nitrite + nitrate (hereafter ), ammonium (), and soluble reactive phosphorus (SRP) analyses, and 1.5-mL samples of stream water filtered through a 0.2-μm filter for dissolved total Fe and dissolved ferrous Fe analyses. We immediately added 40-μl of Ferrozine to the dissolved ferrous Fe samples. All samples were collected in duplicate. TN was quantified using a potassium persulfate digestion (Nydahl, 1978) followed by a cadmium reduction for (APHA, 1998, EPA method 353.2). TP was quantified using a potassium persulfate digestion followed by an ascorbic acid molybdenum reaction for SRP (Murphy and Riley, 1962, EPA method 365.1). and SRP were measured using the methods described above without digestion. was quantified with an automated alkaline phenolhypochlorite reaction followed by spectrophotometric analysis (EPA method 350.1, Solorzano, 1969; APHA, 1998). Dissolved total Fe and dissolved ferrous Fe were measured colorimetrically at a wavelength of 562 nm (Stookey, 1970). The method detection limits were 8.0 μg L–1 for TN, 12.0 μg L–1 for TP, 6.0 μg L–1 for , 4.0 μg L–1 for , 1.0 μg L–1 for SRP, and 17.1 μg L–1 for dissolved total Fe and dissolved ferrous Fe. We performed one-way ANOVAs with land-use as a fixed effect and catchment nested within land-use to test for differences in nutrients concentrations between montane and urban streams.
To examine the potential for pharmaceutical pollution to influence biofilm bacteria, we measured ambient concentrations of pharmaceuticals by deploying Passive Organic Chemical Integrative Samplers (POCIS) in streams at the urban site in each catchment (Alvarez et al., 2004). Pharmaceutical POCIS containing Waters Oasis HLBTM polymer sorbent were obtained from Environmental Sampling Technologies (EST, Inc., St. Joseph, MO, United States) and deployed at each location in a stainless-steel strainer (Zicome Collapsible Vegetable Steamer 7.5 × 7.3 × 2.3 inches) except at one location where a protective covering was fashioned out of plastic landscape cloth. POCIS are typically deployed in the field for an average of 30 days, although deployments have ranged from weeks to months (Alvarez, 2010). Accordingly, we deployed our POCIS for 20–26 days during the summer and fall of 2015. The concentrations of 19 pharmaceuticals were measured with high performance liquid chromatography combined with tandem mass spectrometry (Method 1694; EPA, 2007) by the University of Nebraska Water Sciences Laboratory.
Results
Bacterial Core Communities
Overall, biofilm bacterial cores were composed of mostly a few abundant taxa (Figures 1, 2). Bacterial cores at montane and urban sites contained a mean (± SEM) of 385 ± 37 and 576 ± 73 of OTUs, which represented only 2 and 3% of OTUs in the total bacterial community, respectively (Supplementary Table 1). Abundant taxa comprised 34 ± 3.4 and 24 ± 2.4% of bacterial core OTUs at montane and urban sites, respectively (Supplementary Table 1). In comparison, abundant taxa comprised only 0.81% of OTUs in the total bacterial community. The most abundant bacterial core phyla included Proteobacteria (mean relative abundance ± SEM: montane = 62 ± 3.8%; urban = 60 ± 2.9%), Bacteroidetes (montane = 27 ± 3.2%; urban = 16 ± 1.5%), Verrucomicrobia (montane = 4.4 ± 0.34%; urban = 7.2 ± 1.10%), and Cyanobacteria (montane = 3.9 ± 0.75%; urban = 6.5 ± 0.54%; Figures 1, 2). Control, nutrient, pharmaceutical, and nutrient plus pharmaceutical bacterial cores shared 35 and 32% of their OTUs at montane and urban sites, respectively (Figure 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/526545/fmicb-11-526545-HTML/image_m/fmicb-11-526545-g001.jpg)
Figure 1. Krona plots of biofilm bacterial cores in control, nutrient (Nu), pharmaceutical (Ph), and nutrient plus pharmaceutical (NuPh) treatments in montane streams. Bacterial cores defined as OTUs present in more than 75% of samples of a given contaminant treatment in montane streams across the three catchments. Inner rings, in order from smallest to largest, show percentages of individual taxa at phylum, class, order, and family levels and colors indicate different phyla. The outer ring shows the percentage of individual taxa at the OTU level and colors indicate abundant (≥0.1% relative abundance of total bacterial community) and rare OTUs (<0.1% relative abundance of total bacterial community). To simplify Krona plots, we only included taxa belonging to eight most abundant phyla across all bacterial cores, which comprised at least 99% of the abundance of bacterial cores at montane sites. Bacterial cores based on 16S rRNA gene community libraries (97% similarity cut-off). Top families are labeled with the following abbreviations: Comam. = Comamonadaceae, Flavo. = Flavobacteriaceae, Oscilla. = Oscillatoriales, and Pseudo. = Pseudomonadaceae.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/526545/fmicb-11-526545-HTML/image_m/fmicb-11-526545-g002.jpg)
Figure 2. Krona plots of biofilm bacterial cores on control, nutrient (Nu), pharmaceutical (Ph), and nutrient plus pharmaceutical (NuPh) treatments at urban sites. Krona plots are described in Figure 1. The eight phyla comprised at least 99% of the abundance of bacterial cores at urban sites. Top families are labeled with the following abbreviations: Comam. = Comamonadaceae, Flavo. = Flavobacteriaceae, Oscilla. = Oscillatoriales, and Pseudo. = Pseudomonadaceae.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/526545/fmicb-11-526545-HTML/image_m/fmicb-11-526545-g003.jpg)
Figure 3. Venn diagrams of shared and unique bacterial core OTUs in control, nutrient (Nu), pharmaceutical (Ph), and nutrient plus pharmaceutical (NuPh) treatments in montane (left diagram) and urban streams (right diagram). Values are the number of OTUs in each section and section areas are proportional to the number of OTUs. Bacterial cores in montane streams contained a total of 611 OTUs and bacterial cores in urban streams contained a total of 959 OTUs.
Land-Use and Contaminant Effects on Bacterial Core Structure
Bacterial cores were primarily structured by urban and montane land-use and secondarily by nutrient and pharmaceutical additions. In the PCoA ordination, the two different land-uses clearly separated bacterial cores along the primary axis (Figure 4). In general, bacterial cores at montane sites in the three catchments were compositionally distinct, while bacterial cores at urban sites were more similar to each other. PERMANOVA analysis confirmed the separation in ordination space, with land-use explaining 27% of the variation (land-use: F1,66 = 27, p < 0.001, R2 = 0.27). In comparison to land-use, nutrients and pharmaceutical additions had less of an impact but still helped separate bacterial cores along axis 2 and contributed to the formation of distinct bacterial cores based on PERMANOVA results (Nu: F1,66 = 5.9, p < 0.001, R2 = 0.06, Ph: F1,66 = 2.4, p = 0.02, R2 = 0.02).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/526545/fmicb-11-526545-HTML/image_m/fmicb-11-526545-g004.jpg)
Figure 4. Principle coordinate analysis (PCoA) generated using Bray–Curtis distances on a sample × OTU matrix of 16S rRNA gene community libraries. Bacterial cores in control, nutrient addition, pharmaceutical addition, and nutrient plus pharmaceutical addition treatments are represented as different shapes with montane and urban streams in the Logan, Red Butte (RB), and Middle Provo (Provo) catchments represented as different colors.
Nutrient Effects on Bacterial Cores
Nutrients enhanced the dominance of the Pseudomonadaceae and Comamonadaceae in both nutrient and nutrient plus pharmaceutical bacterial cores (Figures 1, 2). The abundance of the Pseudomonadaceae, a family containing seven core taxa, was at least 5- and 3-fold higher in nutrient and nutrient plus pharmaceutical cores, respectively, than in control cores (Figures 1, 2). Comamonadaceae, which contained the highest abundance of OTUs in control bacterial cores and contained 44 core taxa, was 1.7- and 1.5-more abundant in nutrient and nutrient plus pharmaceutical cores as compared to controls at urban sites, respectively (Figure 1). Nutrients depressed bacterial core richness by at least 19 and 26% at montane and urban sites, respectively (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/526545/fmicb-11-526545-HTML/image_m/fmicb-11-526545-g005.jpg)
Figure 5. Richness of bacterial cores in control, nutrient addition, pharmaceutical addition, and nutrient plus pharmaceutical addition treatments in montane (left panel) and urban streams (right panel). Values are means (n = 8–9) shown with standard error of the mean. In two-way ANOVAs, nutrients and pharmaceuticals had interactive effects on bacterial core richness, with nutrients decreasing and pharmaceuticals increasing richness in montane (Nu addition × Ph addition, F1,24 = 11.7, p = 0.002; Nu addition, F1,24 = 677, p < 0.001; Ph addition, F1,24 = 70.6, p < 0.001) and urban streams (Nu addition × Ph addition, F1,22 = 188, p < 0.001; Nu addition, F1,22 = 846, p < 0.001; Ph addition, F1,22 = 64.7, p < 0.001). All treatments were significantly different within each type of land-use.
Pharmaceutical Effects on Bacterial Cores
Pharmaceutical effects generally contrasted with those of nutrients, as pharmaceuticals increased bacterial core richness (Figure 5) but had few effects on abundant families (Figures 1, 2). Pharmaceuticals had the most dramatic effect on richness in urban streams, where pharmaceutical exposure increased bacterial core richness 38% above controls (Figure 5B). In montane streams, pharmaceuticals only increased bacterial core richness by 6% above controls (Figure 5A). Pharmaceuticals had a minimal effect on abundant families with the exception of Pseudomonadaceae, which was 4.8- and 1.5-times more abundant in pharmaceutical as compared to control bacterial cores at urban and montane sites, respectively (Figures 1, 2).
Pharmaceutical bacterial cores were distinct, containing many unique taxa associated with contaminant-degradation. Based on the Venn diagrams of the OTUs found in contaminant treatment bacterial cores at urban sites, 21% of OTUs were unique to the pharmaceutical bacterial cores, while only 7.8, 3.7, and 1.5% were unique to the control, nutrient, and nutrient plus pharmaceutical bacterial cores, respectively (Figure 3B). The pharmaceutical bacterial core was less unique at montane sites, where 14, 11, 2.5, and 4.0% of OTUs were unique to pharmaceutical, control, nutrient, and nutrient plus pharmaceutical bacterial cores, respectively (Figure 3A). Top families of unique pharmaceutical bacterial core taxa at urban sites (Figure 6) included DR-16 (Betaproteobacteria; 25% of unique pharmaceutical bacterial core OTUs), WCHB1-32 (Bacteroidetes; 6.8%), and Leptotrichiaceae (Fusobacteria; 6.8%). At montane sites, top families of unique pharmaceutical bacterial core taxa included unclassified Bacteroidetes (10.0%), Oscillatoriales (9.8%), and Sphingobacteriales env.ops 17 (Bacteroidetes; 6.9%).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/526545/fmicb-11-526545-HTML/image_m/fmicb-11-526545-g006.jpg)
Figure 6. Krona plots of unique taxa in bacterial cores in control, nutrient (Nu), pharmaceutical (Ph), and nutrient plus pharmaceutical (NuPh) treatments in urban streams. Rings, in order from smallest to largest, show percentages of individual taxa at phylum, class, order, family, and OTU levels. Colors indicate different phyla. Note that the phyla and color key are different from those in Figures 1, 2. To simplify Krona plots, we only included taxa belonging to eight most abundant phyla of unique taxa, which comprised at least 96% of the abundance of unique taxa. Top families of unique taxa in each bacterial core are labeled using the following abbreviations: Aero. = Aeromonadaceae, Chitin. = Chitinophagaceae, Cyto. = Cytophagaceae, DR-16 = DR-16 family, env.OPS 17 = Sphingobacteriales env.OPS 17, Lepto. = Leptotrichiaceae, Morax. = Moraxellaceae, OM190 = OM190 family, Oscilla. = Oscillatoriales, Plancto. = Planctomycetaceae, WCHB1-32 = WCHB1-32 family, Xantho. = Xanthomonadaceae.
Nutrient Plus Pharmaceutical Effects on Bacterial Cores
Although the composition of nutrient plus pharmaceutical bacterial cores closely resembled that of nutrient bacterial cores, nutrient plus pharmaceutical bacterial cores also contained taxa not found in other contaminant treatment bacterial cores and the unique taxa in the nutrient plus pharmaceutical core were distinct from unique taxa in the control, nutrient, and pharmaceutical cores at urban sites (Figure 6). Oscillatoriales comprised 80 and 11% of the OTUs unique to nutrient plus pharmaceutical bacterial cores at urban and montane sites, respectively (Figure 6). Other major families of unique nutrient plus pharmaceutical bacterial core taxa at urban sites included Sphingobacteriales env.OPS 17 (Bacteroidetes; 3.7% of unique nutrient plus pharmaceutical bacterial core OTUs), Cytophagaceae (Bacteroidetes; 2.5%), and Nitrosomonadaceae (Betaproteobacteria; 2.5%), while major families at montane sites included Verrucomicrobiaceae (Verrucomicrobia; 19%), Rhodobacteraceae (Alphaproteobacteria; 17%), and unclassified Burkholderiales (Betaproteobacteria; 6.6%).
Ambient Nutrients and Pharmaceuticals
Urban biofilms were exposed to distinct signatures of nutrient and pharmaceutical pollution associated with the different types of land-use in the catchment. Although variable, ambient nutrient concentrations were generally higher in urban compared to montane streams (Table 1): dissolved and total N and P and total dissolved Fe concentrations were significantly greater in urban as compared to montane streams (all models one-way ANOVAs with land-use as the predictor variable: p < 0.02). The increase in ambient nutrient concentrations was likely due to human land-use, as differences in the form and ratio of N and P across the three catchments reflected differences in the intensity of urbanization in the three watersheds. The urban stream in the Logan catchment, which is impacted by a mixture of agricultural and urban land-use, had greater concentrations and DIN:SRP ratios than urban streams in the Red Butte and Middle Provo catchments, potentially reflecting chemical fertilizer inputs associated with agriculture (Kaushal et al., 2011). In contrast, P concentrations increased between montane and urban streams in Red Butte Creek and Middle Provo, potentially reflecting the influence of urban land-use, which generally increases P inputs due to sewage and septic inputs (Paul and Meyer, 2001; Walsh et al., 2005; Kaushal et al., 2011). In addition, the Middle Provo urban stream is also the only stream in our study located downstream of a wastewater treatment plant and the increase in P concentrations between montane and urban streams in this catchment may reflect the influence of treated wastewater inputs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/526545/fmicb-11-526545-HTML/image_m/fmicb-11-526545-t001.jpg)
Table 1. Ambient nutrient concentrations in montane and urban streams in the Logan River, Red Butte Creek, and Middle Provo River catchments during 2015 CES bioassays.
Similarly, urban biofilms were also exposed to different forms of pharmaceutical pollution associated with human land-use in the three catchments (Table 2). Caffeine, one of the contaminants in our pharmaceutical treatment, and 1,7-dimethylzanthine, a metabolite of caffeine, were present in urban streams in all three catchments. Diphenhydramine, the other contaminant in our pharmaceutical treatment, was only detected in the Logan River and Red Butte Creek catchments. The different types of human land-use created distinct chemical signatures of human- and livestock-associated pharmaceuticals. Sulfa antibiotics (sulfamethazine and sulfamethoxazole), pharmaceuticals primarily associated with livestock (Mathew et al., 2007), were most frequently detected in the Logan and Middle Provo catchments, which are both impacted by agricultural land-use. Pharmaceuticals associated with human consumption, including 1,7-dimethylxanthine (a metabolite of caffeine), acetaminophen, caffeine, carbamazepine, and cotinine (a metabolite of nicotine) were detected in all three catchments, reflecting anthropogenic impacts found in both rural and urban environments.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/526545/fmicb-11-526545-HTML/image_m/fmicb-11-526545-t002.jpg)
Table 2. Concentrations of pharmaceuticals that accumulated in Polar Organic Chemical Integrated Samplers (POCIS) in urban streams in the Logan River, Red Butte Creek, Middle Provo River catchments during summer (June–July) and fall [September (Sept.)–October (Oct.)].
Contaminant Effects on Biofilm Biomass
Urban biofilm biomass was nutrient-limited. Nutrient additions increased chlorophyll a concentrations by approximately 250% (Figure 7A; Nu: F1,47 = 36, p < 0.001), but increased AFDM by only 25% (Figure 7B; Nu: F1,49 = 5.2, p = 0.03). Pharmaceuticals in urban streams and nutrients and pharmaceuticals in montane streams did not influence biofilm biomass.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/526545/fmicb-11-526545-HTML/image_m/fmicb-11-526545-g007.jpg)
Figure 7. Chlorophyll a concentrations (A) and ash-free dry mass (B) of biofilms grown on control, nutrient addition, pharmaceutical addition, and nutrient plus pharmaceutical addition treatments in urban streams. Values are means ± SEM (n = 3). Treatments with different lowercase letters were significantly different from each other (all p ≤ 0.03).
Discussion
Nutrient and pharmaceutical additions created distinct bacterial cores in stream biofilms at montane and urban sites across the three catchments. Nutrients enhanced the dominance of Pseudomonadaceae and Comamonadaceae in bacterial cores and depressed bacterial core richness. Pseudomonadaceae and Comamonadaceae both contain many opportunistic copiotrophs often dominating nutrient-rich environments (Pinhassi and Berman, 2003; Yu et al., 2017; Jani et al., 2018; Roberto et al., 2018) and are members of Gammaproteobacteria and Betaproteobacteria, respectively, both copiotrophic Proteobacterial classes (Shi et al., 1999; Fierer et al., 2007a). Additions of N, P, and Fe, the nutrients included in our treatment, likely led to the dominance of these copiotrophs, which were able to outcompete other taxa for these nutrients, ultimately reducing richness (Mittelbach et al., 2001; Besemer, 2015). Nutrient pollution has previously demonstrated a negative effect on bacterial richness within a nutrient-rich urban stream, where nitrate concentrations were four-times greater than pristine streams and biofilm bacterial richness and evenness estimates were visibly lower and under the dominance of several species (Lear et al., 2009). While the results of our CES experiment suggest nutrient enrichment may depress bacterial core richness, the composition of biofilms grown on CES likely differed from natural biofilms at our study sites. A previous CES study which employed cellulose sponges a biofilm colonization substrates found that bacterial and fungal diversity were on the low end of value reported for stream leaf litter and sediments, potentially due to the structural homogeneity of the CES substrate, short incubation time, and/or microbial identification methodology (Hoellein et al., 2010). Research examining the effects of nutrient enrichment on the bacterial core of natural biofilms is an important avenue of future research. Diminished bacterial richness associated with nutrient additions may have important implications for stream ecosystems, as the loss of bacterial taxa may also reduce genetic diversity and, consequently, biofilm resilience to disturbance (Sogin et al., 2006; Pedrós-Alió, 2012).
The effects of pharmaceuticals largely contrasted those of nutrients, as pharmaceuticals increased bacterial core richness and fostered unique taxa associated with contaminant tolerance and/or degradation. The pharmaceuticals used in our study, caffeine and diphenhydramine, may have served as a C and/or energy source for taxa capable of degrading one or both compounds. Indeed, pharmaceuticals increased the relative abundance of Pseudomonadaceae, a family which contains many taxa known for their ability to degrade toxic organic molecules (Palleroni, 2010) and includes the majority of bacteria capable of degrading caffeine (Summers et al., 2015). Experimental diphenhydramine additions have previously been shown to increase the abundance Pseudomonas sp. in stream biofilms (Rosi-Marshall et al., 2013). Further, a small percentage (0.9%) of unique pharmaceutical bacterial core taxa at montane sites belonged to the Anaerolineaceae (Chloroflexi), a family associated with diphenhydramine degradation in anaerobic digester sludge (Wolfson et al., 2018). Many taxa that were unique to pharmaceutical bacterial cores are associated with the breakdown of complex C substrates or are commonly found in highly contaminated environments. DR-16, which contained the highest abundance of unique pharmaceutical bacterial core taxa at urban sites, resides in a quinoline degrading microbial consortium and may help degrade aromatic compounds (Wu et al., 2015). Other highly abundant families of unique pharmaceutical bacterial core taxa at urban sites included WCHB1-32, a family found in a hydrocarbon- and chlorinated-solvent-contaminated aquifer undergoing bioremediation (Dojka et al., 1998) and Leptotrichiaceae, a family containing pathogenic bacteria (Lory, 2014) associated with sewage infrastructure (Shanks et al., 2013; Gonzalez-Martinez et al., 2018) that often contains pharmaceuticals. However, pharmaceutical release rates from CES likely exposed CES biofilms to pharmaceutical concentrations that exceeded environmental conditions by approximately four orders of magnitude (Shaw et al., 2015; Costello et al., 2016), likely resulting in more pronounced pharmaceutical effects. Future CES studies that includes environmentally-relevant pharmaceutical concentrations and dose-response curves (Costello et al., 2016) could provide insight into the ecological effects of pharmaceutical pollution and how these effects change with concentration.
Our results highlight the potential interplay between pharmaceutical pollution, diversity, and pharmaceutical degradation. The increase in bacterial core richness associated with pharmaceutical addition is notable, as contaminant degradation rates may increase with microbial community diversity (Dell’Anno et al., 2012; Hernandez-Raquet et al., 2013; Delgado-Baquerizo et al., 2016). Bacterial biodiversity may provide a rich set of taxa that can act as a consortium supporting specialized metabolic pathways required to degrade complex organic pollutants (Fuentes et al., 2014; Jousset et al., 2017). However, we are only able to hypothesize that bacteria degraded caffeine and diphenhydramine. The pharmaceutical-induced compositional shifts that we found may also result from toxicity and/or tolerance to these contaminants. Stable isotope probing, which tracks the incorporation of an isotopically-labeled substrate of interest into specific taxa within a microbial community (Neufeld et al., 2007), may identify potential pharmaceutical-degrading taxa. Future studies that examine the potential role of biofilm bacterial taxa in pollutant degradation may reveal the importance of this group in stream ecosystem function and resilience to anthropogenic pollution (Jousset et al., 2017).
Nutrients had much stronger effects than pharmaceuticals on bacterial cores when the two contaminant classes were added in combination. The composition of the top bacterial taxa, richness, and diversity in the nutrient plus pharmaceutical bacterial cores largely resembled those in the nutrient bacterial cores. Unexpectedly, nutrient plus pharmaceutical bacterial cores contained taxa not found in nutrient or pharmaceutical bacterial cores. The vast majority of unique nutrient plus pharmaceutical bacterial core taxa in urban streams belonged to the order Oscillatoriales (Cyanobacteria), suggesting the taxa may possess superior competitive abilities under combined nutrient and pharmaceutical enriched conditions. Indeed, Oscillatoriales outcompeted other bacteria for energy and nutrients under N-rich conditions in salt marsh communities (Kearns et al., 2016) and increased in response to pharmaceutical pollution in a wastewater-impacted stream (Corcoll et al., 2014). Further research to identify bacterial taxa that thrive in multi-pollutant conditions may help identify biofilm responses to the variety of contaminants associated with urbanization (Walsh et al., 2005).
Last, we were most surprised to find that in contrast to our hypothesis, prior exposure to nutrient and pharmaceutical pollution in urban streams did not dampen the effects of our contaminant treatments on biofilm bacteria. Instead, pharmaceutical additions elevated bacterial core richness by 37% at urban sites and only 6% at montane sites. Existing nutrient and pharmaceutical pollution may have acted as an environmental filter selecting for taxa ready to capitalize on caffeine and diphenhydramine. Similarly, experimental additions of the antibiotic ciprofloxacin had an immense effect on community composition of urban biofilms in streams along a suburban to urban gradient in Maryland, United States (Rosi et al., 2018), with ciprofloxacin increasing the abundance of Pedobacter spp., a genus known for its ability to degrade unusual compounds (Kämpfer, 2010), and an unclassified genus from the family Bradyrhizobiaceae, which included several strains that are able to metabolize halogenated aromatic compounds (Kamal and Wyndham, 1990; Van der Woude et al., 1994; Egland et al., 2001). Thus, urbanization may create a predictable occurrence of pharmaceuticals entering streams that allows contaminant-tolerant bacteria to colonize biofilms.
Conclusion
Both nutrients and pharmaceuticals structured biofilm bacterial cores at montane and urban sites, albeit in distinctly different ways. While N, P, and Fe additions enhanced the dominance of already abundant, copiotrophic heterotrophs, such as Pseudomonadaceae and Comamonadaceae, caffeine and diphenhydramine additions supported rich bacterial cores containing unique taxa associated with contaminant-degradation. In addition, nutrient plus pharmaceutical bacterial cores contained unique taxa, many of which belonged to Oscillatoriales, suggesting that these taxa possess superior competitive abilities to scavenge nutrients and tolerate and/or capitalize on multiple contaminants at a single time. Surprisingly, exposure to ambient nutrient and pharmaceutical pollution failed to reduce the effects of these contaminants on bacterial cores at urban sites, highlighting the potential for pollution to act as an environmental filter selecting for taxa that may readily respond to urban water conditions.
Data Availability Statement
The datasets generated for this study can be found in the datasets generated and analyzed for this study are publicly available online at the Hydroshare database: http://www.hydroshare.org/resource/1f9c2b8d4fbd4c849e8e492cb0fd9198.
Author Contributions
EO, MB, ER, DL, and ZA contributed to the study design. EO and DL carried out the experiment and lab work. EO and TS carried out the analyses. EO, MB, and ZA wrote the manuscript. All authors substantially contributed to commenting and revising it.
Funding
This project was supported by the Utah State University Presidential Doctoral Research Fellowship Program, the NSF-funded iUTAH project (NSF EPSCoR award IIA 1208732), and the Ecology Center at Utah State University.
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.
Acknowledgments
We thank members of the Utah State University Aquatic Biogeochemistry and Brigham Young University Microbial and Ecosystem Ecology Laboratories.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.526545/full#supplementary-material
Footnotes
References
Aanderud, Z., Jones, S., Fierer, N., and Lennon, J. T. (2015). Resuscitation of the rare biosphere contributes to pulses of ecosystem activity. Front. Microbiol. 6:24. doi: 10.3389/fmicb.2015.00024
Aanderud, Z. T., Saurey, S., Ball, B. A., Wall, D. H., Barrett, J. E., Muscarella, M. E., et al. (2018). Stoichiometric shifts in soil C:N:P promote bacterial taxa dominance, maintain biodiversity, and deconstruct community assemblages. Front. Microbiol. 9:1401. doi: 10.3389/fmicb.2018.01401
Alvarez, D. A. (2010). Guidelines for the Use of the Semipermeable Membrane Device (SPMD) and the Polar Organic Chemical Integrative Sampler (POCIS) in Environmental Monitoring Studies, Vol. 1. Reston, VA: US Geological Survey, 1–28. doi: 10.3133/tm1D4
Alvarez, D. A., Petty, J. D., Huckins, J. N., Jones-Lepp, T. L., Getting, D. T., Goddard, J. P., et al. (2004). Development of a passive, in situ, integrative sampler for hydrophilic organic contaminants in aquatic environments. Environ. Toxicol. Chem. 23, 1640–1648. doi: 10.1897/03-603
Anderson, M. J. (2001). A new method for non-parametric multivariate analysis of variance. Austral Ecol. 26, 32–46. doi: 10.1111/j.1442-9993.2001.01070.pp.x
APHA (1998). Standard Methods for the Examination of Water and Wastewater, 20th Edn. Washington, DC: American Public Health Association.
aus der Beek, T., Weber, F.-A., Bergmann, A., Hickmann, S., Ebert, I., Hein, A., et al. (2016). Pharmaceuticals in the environment—global occurrences and perspectives. Environ. Toxicol. Chem. 35, 823–835. doi: 10.1002/etc.3339
Battin, T. J., Besemer, K., Bengtsson, M. M., Romani, A. M., and Packmann, A. I. (2016). The ecology and biogeochemistry of stream biofilms. Nat. Rev. Microbiol. 14, 251–263. doi: 10.1038/nrmicro.2016.15
Battin, T. J., Sloan, W. T., Kjelleberg, S., Daims, H., Head, I. M., Curtis, T. P., et al. (2007). Microbial landscapes: new paths to biofilm research. Nat. Rev. Microbiol. 5, 76–81. doi: 10.1038/nrmicro1556
Besemer, K. (2015). Biodiversity, community structure and function of biofilms in stream ecosystems. Res. Microbiol. 166, 774–781. doi: 10.1016/j.resmic.2015.05.006
Besemer, K., Peter, H., Logue, J. B., Langenheder, S., Lindström, E. S., Tranvik, L. J., et al. (2012). Unraveling assembly of stream biofilm communities. ISME J. 6, 1459–1468. doi: 10.1038/ismej.2011.205
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Huntley, J., Fierer, N., et al. (2012). Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 6, 1621–1624. doi: 10.1038/ismej.2012.8
Capps, K. A., Booth, M. T., Collins, S. M., Davison, M. A., Moslemi, J. M., El-Sabaawi, R. W., et al. (2011). Nutrient diffusing substrata: a field comparison of commonly used methods to assess nutrient limitation. J. North Am. Benthol. Soc. 30, 522–532. doi: 10.1899/10-146.1
Carey, R. O., Hochmuth, G. J., Martinez, C. J., Boyer, T. H., Dukes, M. D., Toor, G. S., et al. (2013). Evaluating nutrient impacts in urban watersheds: challenges and research opportunities. Environ. Pollut. 173, 138–149. doi: 10.1016/j.envpol.2012.10.004
Cole, J. R., Wang, Q., Cardenas, E., Fish, J., Chai, B., Farris, R. J., et al. (2009). The ribosomal database project: improved alignments and new tools for rrna analysis. Nucleic Acids Res. 37, D141–D145. doi: 10.1093/nar/gkn879
Corcoll, N., Acuña, V., Barceló, D., Casellas, M., Guasch, H., Huerta, B., et al. (2014). Pollution-induced community tolerance to non-steroidal anti-inflammatory drugs (NSAIDs) in fluvial biofilm communities affected by WWTP effluents. Chemosphere 112, 185–193. doi: 10.1016/j.chemosphere.2014.03.128
Costello, D. M., Rosi-Marshall, E. J., Shaw, L. E., Grace, M. R., and Kelly, J. J. (2016). A novel method to assess effects of chemical stressors on natural biofilm structure and function. Freshw. Biol. 61, 2129–2140. doi: 10.1111/fwb.12641
Costerton, J. W., Geesey, G. G., and Cheng, K.-J. (1978). How bacteria stick. Sci. Am. 238, 86–95. doi: 10.1038/scientificamerican0178-86
Dastidar, S. G., Saha, P. K., Sanyamat, B., and Chakrabarty, A. N. (1976). Antibacterial activity of ambodryl and benadryl. J. Appl. Bacteriol. 41, 209–214. doi: 10.1111/j.1365-2672.1976.tb00621.x
Delgado-Baquerizo, M., Giaramida, L., Reich, P. B., Khachane, A. N., Hamonts, K., Edwards, C., et al. (2016). Lack of functional redundancy in the relationship between microbial diversity and ecosystem functioning. J. Ecol. 104, 936–946. doi: 10.1111/1365-2745.12585
Dell’Anno, A., Beolchini, F., Rocchetti, L., Luna, G. M., and Danovaro, R. (2012). High bacterial biodiversity increases degradation performance of hydrocarbons during bioremediation of contaminated harbor marine sediments. Environ. Pollut. 167, 85–92. doi: 10.1016/j.envpol.2012.03.043
Dojka, M. A., Hugenholtz, P., Haack, S. K., and Pace, N. R. (1998). Microbial diversity in a hydrocarbon-and chlorinated-solvent-contaminated aquifer undergoing intrinsic bioremediation. Appl. Environ. Microbiol. 64, 3869–3877. doi: 10.1128/AEM.64.10.3869-3877.1998
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200.
Egland, P. G., Gibson, J., and Harwood, C. S. (2001). Reductive, coenzyme A-mediated pathway for 3-Chlorobenzoate degradation in the phototrophic Bacterium Rhodopseudomonas palustris. Appl. Environ. Microbiol. 67, 1396–1399. doi: 10.1128/AEM.67.3.1396-1399.2001
EPA (2007). Method 1694: Pharmaceuticals and Personal Care Products in Water, Soil, Sediments, and Biosolids by HPLC/MS/MS. EPA-821-R-08-002. Washington, DC: United States Environmental Protection Agency.
Falkowski, P., Scholes, R. J., Boyle, E. E. A., Canadell, J., Canfield, D., Elser, J., et al. (2000). The global carbon cycle: a test of our knowledge of earth as a system. Science 290, 291–296. doi: 10.1126/science.290.5490.291
Fierer, N., Bradford, M. A., and Jackson, R. B. (2007a). Toward an ecological classification of soil bacteria. Ecology 88, 1354–1364. doi: 10.1890/05-1839
Fierer, N., Morse, J. L., Berthrong, S. T., Bernhardt, E. S., and Jackson, R. B. (2007b). Environmental controls on the landscape-scale biogeography of stream bacterial communities. Ecology 88, 2162–2173. doi: 10.1890/06-1746.1
Fuentes, S., Méndez, V., Aguila, P., and Seeger, M. (2014). Bioremediation of petroleum hydrocarbons: catabolic genes, microbial communities, and applications. Appl. Microbiol. Biotechnol. 98, 4781–4794. doi: 10.1007/s00253-014-5684-9
Geesey, G. G., Mutch, R., Costerton, J. T., and Green, R. B. (1978). Sessile bacteria: an important component of the microbial population in small mountain streams. Limnol. Ocean. 23, 1214–1223. doi: 10.4319/lo.1978.23.6.1214
Gonzalez-Martinez, A., Sihvonen, M., Munoz-Palazon, B., Rodriguez-Sanchez, A., Mikola, A., and Vahala, R. (2018). Microbial ecology of full-scale wastewater treatment systems in the polar arctic circle: Archaea, Bacteria and Fungi. Sci. Rep. 8:2208. doi: 10.1038/s41598-018-20633-5
Hernandez-Raquet, G., Durand, E., Braun, F., Cravo-Laureau, C., and Godon, J.-J. (2013). Impact of microbial diversity depletion on xenobiotic degradation by sewage-activated sludge. Environ. Microbiol. Rep. 5, 588–594. doi: 10.1111/1758-2229.12053
Hoellein, T. J., Tank, J. L., Kelly, J. J., and Rosi-Marshall, E. J. (2010). Seasonal variation in nutrient limitation of microbial biofilms colonizing organic and inorganic substrata in streams. Hydrobiologia 649, 331–345. doi: 10.1007/s10750-010-0276-x
Jani, K., Ghattargi, V., Pawar, S., Inamdar, M., Shouche, Y., and Sharma, A. (2018). Anthropogenic activities induce depletion in microbial communities at urban sites of the river Ganges. Curr. Microbiol. 75, 79–83. doi: 10.1007/s00284-017-1352-5
Jousset, A., Bienhold, C., Chatzinotas, A., Gallien, L., Gobet, A., Kurm, V., et al. (2017). Where less may be more: how the rare biosphere pulls ecosystems strings. ISME J. 11, 853–862. doi: 10.1038/ismej.2016.174
Kamal, V. S., and Wyndham, R. C. (1990). Anaerobic phototrophic metabolism of 3-chlorobenzoate by Rhodopseudomonas palustris WS17. Appl. Environ. Microbiol. 56, 3871–3873. doi: 10.1128/AEM.56.12.3871-3873.1990
Kaushal, S. S., Groffman, P. M., Band, L. E., Elliott, E. M., Shields, C. A., and Kendall, C. (2011). Tracking nonpoint source nitrogen pollution in human-impacted watersheds. Environ. Sci. Technol. 45, 8225–8232. doi: 10.1021/es200779e
Kearns, P. J., Angell, J. H., Howard, E. M., Deegan, L. A., Stanley, R. H., and Bowen, J. L. (2016). Nutrient enrichment induces dormancy and decreases diversity of active bacteria in salt marsh sediments. Nat. Commun. 7:12881. doi: 10.1038/ncomms12881
Kolpin, D. W., Furlong, E. T., Meyer, M. T., Thurman, E. M., Zaugg, S. D., Barber, L. B., et al. (2002). Pharmaceuticals, hormones, and other organic wastewater contaminants in U.S. streams, 1999-2000:? a national reconnaissance. Environ. Sci. Technol. 36, 1202–1211. doi: 10.1021/es011055j
Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K., and Schloss, P. D. (2013). 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. doi: 10.1128/AEM.01043-13
Lahti, L., Shetty, S., and Blake, T. (2017). Tools for Microbiome Analysis in R. Microbiome Package Version. Available online at: http://microbiome.github.io/microbiome (accessed December 10, 2018).
Larson, C. A., Mirza, B., Rodrigues, J. L. M., and Passy, S. I. (2018). Iron limitation effects on nitrogen-fixing organisms with possible implications for cyanobacterial blooms. FEMS Microbiol. Ecol. 94:fiy046. doi: 10.1093/femsec/fiy046
Larsson, J. (2019). eulerr: Area-Proportional Euler and Venn Diagrams with Ellipses. R package version 5.1.0.
Lawrence, J. R., Swerhone, G. D., Wassenaar, L. I., and Neu, T. R. (2005). Effects of selected pharmaceuticals on riverine biofilm communities. Can. J. Microbiol. 51, 655–669. doi: 10.1139/w05-047
Lear, G., Boothroyd, I. K. G., Turner, S. J., Roberts, K., and Lewis, G. D. (2009). A comparison of bacteria and benthic invertebrates as indicators of ecological health in streams. Freshw. Biol. 54, 1532–1543. doi: 10.1111/j.1365-2427.2009.02190.x
Lock, M. A., Wallace, R. R., Costerton, J. W., Ventullo, R. M., and Charlton, S. E. (1984). River epilithon: toward a structural-functional model. Oikos 42, 10–22. doi: 10.2307/3544604
Lory, S. (2014). “The family Leptotrichiaceae,” in The Prokaryotes: Firmicutes and Tenericutes, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin: Springer), 213–214. doi: 10.1007/978-3-642-30120-9_357
Mathew, A. G., Cissell, R., and Liamthong, S. (2007). Antibiotic resistance in bacteria associated with food animals: a United States perspective of livestock production. Foodborne Pathog. Dis. 4, 115–133. doi: 10.1089/fpd.2006.0066
McMurdie, P. J., and Holmes, S. (2013). phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217
Mittelbach, G. G., Steiner, C. F., Scheiner, S. M., Gross, K. L., Reynolds, H. L., Waide, R. B., et al. (2001). What is the observed relationship between species richness and productivity? Ecology 82, 2381–2396. doi: 10.1890/0012-9658(2001)082[2381:WITORB]2.0.CO;2
Murphy, J., and Riley, J. P. (1962). A modified single solution method for the determination of phosphate in natural waters. Anal. Chim. Acta 27, 31–36. doi: 10.1016/S0003-2670(00)88444-5
Neufeld, J. D., Wagner, M., and Murrell, J. C. (2007). Who eats what, where and when? Isotope-labelling experiments are coming of age. ISME J. 1, 103–110. doi: 10.1038/ismej.2007.30
Nydahl, F. (1978). On the peroxodisulphate oxidation of total nitrogen in waters to nitrate. Water Res. 12, 1123–1130. doi: 10.1016/0043-1354(78)90060-X
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., O’Hara, R. B., Simpson, G. L., et al. (2014). vegan: Community Ecology Package. 2014. R package version 2.0–0.
Ondov, B. D., Bergman, N. H., and Phillippy, A. M. (2011). Interactive metagenomic visualization in a Web browser. BMC Bioinformatics 12:385. doi: 10.1186/1471-2105-12-385
Pal, A., Gin, K. Y.-H., Lin, A. Y.-C., and Reinhard, M. (2010). Impacts of emerging organic contaminants on freshwater resources: review of recent occurrences, sources, fate and effects. Sci. Total Environ. 408, 6062–6069. doi: 10.1016/j.scitotenv.2010.09.026
Palleroni, N. J. (2010). The Pseudomonas story. Environ. Microbiol. 12, 1377–1383. doi: 10.1111/j.1462-2920.2009.02041.x
Paul, M. J., and Meyer, J. L. (2001). Streams in the urban landscape. Annu. Rev. Ecol. Syst. 32, 333–365. doi: 10.1146/annurev.ecolsys.32.081501.114040
Pedrós-Alió, C. (2012). The rare bacterial biosphere. Annu. Rev. Mar. Sci. 4, 449–466. doi: 10.1146/annurev-marine-120710-100948
Pinhassi, J., and Berman, T. (2003). Differential growth response of colony-forming α-and γ-proteobacteria in dilution culture and nutrient addition experiments from Lake Kinneret (Israel), the Eastern Mediterranean Sea, and the Gulf of Eilat. Appl. Environ. Microbiol. 69, 199–211. doi: 10.1128/AEM.69.1.199-211.2003
Pruesse, E., Quast, C., Knittel, K., Fuchs, B. M., Ludwig, W., Peplies, J., et al. (2007). SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 35, 7188–7196. doi: 10.1093/nar/gkm864
R Core Team (2017). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Roberto, A. A., Van Gray, J. B., and Leff, L. G. (2018). Sediment bacteria in an urban stream: spatiotemporal patterns in community composition. Water Res. 134, 353–369. doi: 10.1016/j.watres.2018.01.045
Rosi, E. J., Bechtold, H. A., Snow, D., Rojas, M., Reisinger, A. J., and Kelly, J. J. (2018). Urban stream microbial communities show resistance to pharmaceutical exposure. Ecosphere 9:e02041. doi: 10.1002/ecs2.2041
Rosi-Marshall, E. J., Kincaid, D. W., Bechtold, H. A., Royer, T. V., Rojas, M., and Kelly, J. J. (2013). Pharmaceuticals suppress algal growth and microbial respiration and alter bacterial communities in stream biofilms. Ecol. Appl. 23, 583–593. doi: 10.1890/12-0491.1
Rosi-Marshall, E. J., and Royer, T. V. (2012). Pharmaceutical compounds and ecosystem function: an emerging research challenge for aquatic ecologists. Ecosystems 15, 867–880. doi: 10.1007/s10021-012-9553-z
Rugenski, A. T., Marcarelli, A. M., and Bechtold, H. A. (2008). Effects of temperature and concentration on nutrient release rates from nutrient diffusing substrates. J. North Am. Benthol. Soc. 27, 52–57. doi: 10.1899/07-046R1.1
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Schneider, C. A., Rasband, W. S., and Eliceiri, K. W. (2012). NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675. doi: 10.1038/nmeth.2089
Seitzinger, S. P., Mayorga, E., Bouwman, A. F., Kroeze, C., Beusen, A. H. W., Billen, G., et al. (2010). Global river nutrient export: a scenario analysis of past and future trends. Glob. Biogeochem. Cycles 24:GB0A08. doi: 10.1029/2009GB003587
Shanks, O. C., Newton, R. J., Kelty, C. A., Huse, S. M., Sogin, M. L., and McLellan, S. L. (2013). Comparison of the microbial community structures of untreated wastewaters from different geographic locales. Appl. Environ. Microbiol. 79, 2906–2913. doi: 10.1128/AEM.03448-12
Shaw, L., Phung, C., and Grace, M. (2015). Pharmaceuticals and personal care products alter growth and function in lentic biofilms. Environ. Chem. 12, 301–306. doi: 10.1071/EN14141
Shi, Y., Zwolinski, M. D., Schreiber, M. E., Bahr, J. M., Sewell, G. W., and Hickey, W. J. (1999). Molecular analysis of microbial community structures in pristine and contaminated aquifers: field and laboratory microcosm experiments. Appl. Environ. Microbiol. 65, 2143–2150. doi: 10.1128/AEM.65.5.2143-2150.1999
Sogin, M. L., Morrison, H. G., Huber, J. A., Welch, D. M., Huse, S. M., and Neal, P. R., et al. (2006). Microbial diversity in the deep sea and the underexplored “rare biosphere.” Proc. Natl. Acad. Sci. U.S.A. 103, 12115–12120. doi: 10.1073/pnas.0605127103
Solorzano, L. (1969). Determination of ammonia in natural waters by the phenolhypochlorite method. Limnol. Ocean. 14, 799–801. doi: 10.4319/lo.1969.14.5.0799
Stookey, L. L. (1970). Ferrozine—a new spectrophotometric reagent for iron. Anal. Chem. 42, 779–781. doi: 10.1021/ac60289a016
Summers, R. M., Mohanty, S. K., Gopishetty, S., and Subramanian, M. (2015). Genetic characterization of caffeine degradation by bacteria and its potential applications. Microb. Biotechnol. 8, 369–378. doi: 10.1111/1751-7915.12262
Tank, J. L., Bernot, M. J., and Rosi-Marshall, E. J. (2007). “Nitrogen limitation and uptake,” in Methods in Stream Ecology, 2nd Edn, eds F. R. Hauer and G. A. Lamberti (San Diego, CA: Academic Press), 213–238. doi: 10.1016/B978-012332908-0.50012-7
Tank, J. L., and Dodds, W. K. (2003). Nutrient limitation of epilithic and epixylic biofilms in ten North American streams. Freshw. Biol. 48, 1031–1049. doi: 10.1046/j.1365-2427.2003.01067.x
Van der Woude, B. J., De Boer, M., Van der Put, N. M. J., Van der Geld, F. M., Prins, R. A., and Gottschal, J. C. (1994). Anaerobic degradation of halogenated benzoic acids by photoheterotrophic bacteria. FEMS Microbiol. Lett. 119, 199–207. doi: 10.1111/j.1574-6968.1994.tb06889.x
Van Horn, D. J., Sinsabaugh, R. L., Takacs-Vesbach, C. D., Mitchell, K. R., and Dahm, C. N. (2011). Response of heterotrophic stream biofilm communities to a gradient of resources. Aquat. Microb. Ecol. 64, 149–161. doi: 10.3354/ame01515
Walsh, C. J., Roy, A. H., Feminella, J. W., Cottingham, P. D., Groffman, P. M., and Morgan, R. P. (2005). The urban stream syndrome: current knowledge and the search for a cure. J. North Am. Benthol. Soc. 24, 706–723. doi: 10.1899/04-028.1
Wang, L., Zhang, J., Li, H., Yang, H., Peng, C., Peng, Z., et al. (2018). Shift in the microbial community composition of surface water and sediment along an urban river. Sci. Total Environ. 627, 600–612. doi: 10.1016/j.scitotenv.2018.01.203
Withers, P. J. A., Neal, C., Jarvie, H. P., and Doody, D. G. (2014). Agriculture and eutrophication: where do we go from here? Sustainability 6, 5853–5875. doi: 10.3390/su6095853
Wolfson, S. J., Porter, A. W., Villani, T. S., Simon, J. E., and Young, L. Y. (2018). The antihistamine diphenhydramine is demethylated by anaerobic wastewater microorganisms. Chemosphere 202, 460–466. doi: 10.1016/j.chemosphere.2018.03.093
Wu, G., Zhang, M., Wei, C., Wang, Y., Yao, X., Zhao, L., et al. (2015). Draft genome sequence of Thauera sp. DTG from a denitrifying quinoline degrading microbial consortium. Appl. Environ. Biotechnol. 1, 38–43. doi: 10.18063/AEB.2016.01.012
Keywords: water quality, stream biofilms, nutrients, pharmaceuticals, bacteria, contaminant diffusing substrate
Citation: Ogata EM, Baker MA, Rosi EJ, Smart TB, Long D Jr and Aanderud ZT (2020) Nutrients and Pharmaceuticals Structure Bacterial Core Communities in Urban and Montane Stream Biofilms. Front. Microbiol. 11:526545. doi: 10.3389/fmicb.2020.526545
Received: 14 January 2020; Accepted: 22 September 2020;
Published: 15 October 2020.
Edited by:
Marcelino T. Suzuki, Sorbonne Université, FranceReviewed by:
Natàlia Corcoll, University of Gothenburg, SwedenSarah J. Fischer, University of Colorado Boulder, United States
Copyright © 2020 Ogata, Baker, Rosi, Smart, Long and Aanderud. 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: Elizabeth M. Ogata, bethogata@gmail.com