Differentiating between point and non-point source nutrient loadings of wastewater in an agriculturally impacted area using a hybrid statistical model.

Nutrient contamination from point and non-point sources can lead to harmful consequences, such as algal blooms. Point and non-point nutrient loading estimation is determined using modeling approaches and often require an abundance of variables and observations for calibration. Small rural streams that lack water use designations often lack available data to utilize current modeling strategies. This study proposes the use of a 3-phase hybrid stepwise statistical modeling approach using generalized linear mixed models (GLMM) and a reference stream. Two streams in Central Texas were sampled for 13 months between February 2020 and February 2021, one being impacted by a wastewater treatment plant (WWTP). Dissolved phosphorus (PO4-P), ammonia (NH3-N), nitrite/nitrate (NO2 + NO3-N), total nitrogen (TN), and total phosphorus (TP) were sampled in both streams for each month. Non-point sources of contamination, such as land use/land cover and geomorphology composition, were quantified for both sub-basin drainage areas. Phase I models predicted nutrient concentrations in the reference stream using non-point source variables along with discharge and temporal variables. Best fit models were carried forward to phase II and leveraged a point-source variable, which is a naïve estimate of effluent nutrient concentration in the absence of assimilation. Phase II model coefficients highlight the significance of point-source contamination in predicting nutrient concentration, but overall lacked the ability to make future predictions under new hydrologic regimes from WWTP intensification. Phase III models included deterministically calculating an uptake variable using the relationship between discharge and wetted widths, predicting background non-point concentrations by leveraging phase I models, and calculating future nutrient loadings from WWTP intensification. This approach predicted significant increases in nutrient concentrations under planned WWTP intensification scenarios and decreased uptake efficiencies under the new hydrologic regimes.

Diazotrophy modulates cyanobacteria stoichiometry through functional traits that determine bloom magnitude and toxin production.

Harmful cyanobacterial blooms are an increasing threat to water quality. The interactions between two eco-physiological functional traits of cyanobacteria, diazotrophy (nitrogen (N)-fixation) and N-rich cyanotoxin synthesis, have never been examined in a stoichiometric explicit manner. We explored how a gradient of resource N:phosphorus (P) affects the biomass, N, P stoichiometry, light-harvesting pigments, and cylindrospermopsin production in a N-fixing cyanobacterium, Aphanizomenon. Low N:P Aphanizomenon cultures produced the same biomass as populations grown in high N:P cultures. The biomass accumulation determined by carbon, indicated low N:P Aphanizomenon cultures did not have a N-fixation growth tradeoff, in contrast to some other diazotrophs that maintain stoichiometric N homeostasis at the expense of growth. However, N-fixing Aphanizomenon populations produced less particulate cylindrospermopsin and had undetectable dissolved cylindrospermopsin compared to non-N-fixing populations. The pattern of low to high cyanotoxin cell quotas across an N:P gradient in the diazotrophic cylindrospermopsin producer is similar to the cyanotoxin cell quota response in non-diazotrophic cyanobacteria. We suggest that diazotrophic cyanobacteria may be characterized into two broad functional groups, the N-storage-strategists and the growth-strategists, which use N-fixation differently and may determine patterns of bloom magnitude and toxin production in nature.

Comparative Analysis of RNA-Extraction Approaches and Associated Influences on RT-qPCR of the SARS-CoV-2 RNA in a University Residence Hall and Quarantine Location.

Wastewater-based epidemiology (WBE) provides an early warning and trend analysis approach for determining the presence of COVID-19 in a community and complements clinical testing in assessing the population level, even as viral loads fluctuate. Here, we evaluate combinations of two wastewater concentration methods (i.e., ultrafiltration and composite supernatant-solid), four pre-RNA extraction modifications, and three nucleic acid extraction kits using two different wastewater sampling locations. These consisted of a quarantine facility containing clinically confirmed COVID-19-positive inhabitants and a university residence hall. Of the combinations examined, composite supernatant-solid with pre-RNA extraction consisting of water concentration and RNA/DNA shield performed the best in terms of speed and sensitivity. Further, of the three nucleic acid extraction kits examined, the most variability was associated with the Qiagen kit. Focusing on the quarantine facility, viral concentrations measured in wastewater were generally significantly related to positive clinical cases, with the relationship dependent on method, modification, kit, target, and normalization, although results were variable-dependent on individual time points (Kendall's Tau-b (τ) = 0.17 to 0.6) or cumulatively (Kendall's Tau-b (τ) = -0.048 to 1). These observations can support laboratories establishing protocols to perform wastewater surveillance and monitoring efforts for COVID-19.

Compensatory dynamics of lotic algae break down nonlinearly with increasing nutrient enrichment.

One important mechanism governing the temporal maintenance of biodiversity is asynchrony in co-occurring competitors due to fluctuating environments (i.e., compensatory dynamics). Temporal niche partitioning has evolved in response to predictable oscillations in environmental conditions so that species may offset competition, but we do not yet have a clear understanding of how novel anthropogenic stressors alter seasonal patterns of succession. Many primary producers are nutrient limited, and enrichment may decrease the importance of environmental fluctuations that govern which species are effective competitors under naturally low nutrient regimes. Consequently, elevated nutrient concentrations may synchronize species responses to seasonality. By studying benthic algal assemblages over 2 years from 35 streams that spanned a wide gradient of nutrient enrichment, we found that compensatory dynamics characterizing seasonal succession under natural nutrient regimes broke down at relatively low levels of total phosphorus (P) enrichment (~ 25 µg/L). With increasing P more species were able to coexist at any given time, and seasonal variation in assemblage composition was characterized by synchronous swings in species biovolumes. We also observed much higher instability in assemblage biovolumes with declines in compensatory dynamics, which indicates that anthropogenic alteration of nutrient regimes can affect community stability by changing the dominant mode of seasonal succession. Our findings indicate that compensatory fluctuations of stream algae are driven by seasonality and provide insight about how nutrient enrichment alters evolved drivers of species coexistence.

A comparison of four commercially available RNA extraction kits for wastewater surveillance of SARS-CoV-2 in a college population.

Localized wastewater surveillance has allowed for public health officials to gain a broader understanding of SARS-CoV-2 viral prevalence in the community allowing public health officials time to prepare for impending outbreaks. Given variable levels of virus in the population through public health interventions, proper concentration and extraction of viral RNA is a key step in ensuring accurate detections. With many commercial RNA extraction kits and methodologies available, the performance of 4 different kits were evaluated for SARS-CoV-2 RNA detection in wastewater, specifically focusing on their applicability to lower population densities such as those at university campus dorms. Raw wastewater samples were collected at 4 sites on a college campus over a 24 hour period as a composite sample. Included in these sites was an isolation site that housed students that tested positive for Covid-19 via nasopharyngeal swabs. These samples were analyzed using the following kits: Qiagen All Prep PowerViral DNA/RNA kit, New England BioLabs Monarch RNA MiniPrep Kit, and Zymo Quick RNA-Viral Kit, and the Zymo Quick-RNA Fecal/Soil Microbe MicroPrep Kit. All four sites were processed according to the manufacturer's guidelines. Extractions were then quantified with RT-qPCR one-step reactions using an N2 primer and a linearized plasmid standard. While the Zymo Quick-RNA Fecal/Soil Microbe MicroPrep Kit (also known as the Zymo Environ Water RNA Kit) only recovered approximately 73% (±38%) SARS-CoV-2 RNA compared to the Zymo Quick-RNA Viral kit, it was the most time efficient kit to yield comparable results. This extraction kit had a cumulative processing time of approximately 5 h compared, while the other three kits had processing times between approximately 9 and 9.5 h. Based on the current research, the most effective kits for smaller population densities are pellet based and include a homogenization, inhibitor removal, and RNA preservation step.

Copper and Gold Nanoparticles Increase Nutrient Excretion Rates of Primary Consumers.

Freshwater ecosystems are exposed to engineered nanoparticles through municipal and industrial wastewater-effluent discharges and agricultural nonpoint source runoff. Because previous work has shown that engineered nanoparticles from these sources can accumulate in freshwater algal assemblages, we hypothesized that nanoparticles may affect the biology of primary consumers by altering the processing of two critical nutrients associated with growth and survivorship, nitrogen and phosphorus. We tested this hypothesis by measuring the excretion rates of nitrogen and phosphorus of Physella acuta, a ubiquitous pulmonate snail that grazes heavily on periphyton, exposed to either copper or gold engineered nanoparticles for 6 months in an outdoor wetland mesocosm experiment. Chronic nanoparticle exposure doubled nutrient excretion when compared to the control. Gold nanoparticles increased nitrogen and phosphorus excretion rates more than copper nanoparticles, but overall, both nanoparticles led to higher consumer excretion, despite contrasting particle stability and physiochemical properties. Snails in mesocosms enriched with nitrogen and phosphorus had overall higher excretion rates than ones in ambient (no nutrients added) mesocosms. Stimulation patterns were different between nitrogen and phosphorus excretion, which could have implications for the resulting nutrient ratio in the water column. These results suggest that low concentrations of engineered nanoparticles could alter the metabolism of consumers and increase consumer-mediated nutrient recycling rates, potentially intensifying eutrophication in aquatic systems, for example, the increased persistence of algal blooms as observed in our mesocosm experiment.

A Metagenome-Based Investigation of Gene Relationships for Non-Substrate-Associated Microbial Phosphorus Cycling in the Water Column of Streams and Rivers.

Phosphorus (P) is a nutrient of primary importance in all living systems, and it is especially important in streams and rivers which are sensitive to anthropogenic P inputs and eutrophication. Microbes are accepted as the primary mineralizers and solubilizers of P improving bioavailability for organisms at all trophic levels. Here, we use a genomics approach with metagenome sequencing of 24 temperate streams and rivers representing a total P (TP) gradient to identify relationships between functional genes, functional gene groupings, P, and organisms within the P biogeochemical cycle. Combining information from network analyses, functional groupings, and system P levels, we have constructed a System Relational Overview of Gene Groupings (SROGG) which is a cohesive system level representation of P cycle gene and nutrient relationships. Using SROGG analysis in concert with other statistical approaches, we found that the compositional makeup of P cycle genes is strongly correlated to environmental P whereas absolute abundance of P genes shows no significant correlation to environmental P. We also found orthophosphate (PO43-) to be the dominant factor correlating with system P cycle gene composition with little evidence of a strong organic phosphorous correlation present even in more oligotrophic streams.

Polychlorinated biphenyl (PCB) contamination in Galveston Bay, Texas: Comparing concentrations and profiles in sediments, passive samplers, and fish.

The industrialized portion of the Houston Ship Channel (HSC) is heavily contaminated with anthropogenic contaminants, most prominent of which are the polychlorinated biphenyls (PCBs). This contamination has driven adaptive evolution in a keystone species for Galveston Bay, the Gulf killifish (Fundulus grandis). We investigated the geographical extent of PCB impacts by sampling 12 sites, ranging from the heavily industrialized upper portion of the HSC to Galveston Island. At each site, PCB concentrations and profiles were determined in three environmental compartments: sediment, water (polyethylene passive samplers), and fish tissue (resident Gulf killifish). We observed a steep gradient of PCB contamination, ranging from 4.00 to 100,000 ng/g organic carbon in sediment, 290-110,000 ng/g lipid in fish, and 4.5-2300 ng/g polyethylene in passive samplers. The PCB congener profiles in Gulf killifish at the most heavily contaminated sites were shifted toward the higher chlorinated PCBs and were highly similar to the sediment contamination profiles. In addition, while magnitude of total PCB concentrations in sediment and total fish contamination levels were highly correlated between sites, the relative PCB congener profiles in fish and passive samplers were more alike. This strong correlation, along with a lack of dependency of biota-sediment accumulation factors with total contamination rates, confirm the likely non-migratory nature of Gulf killifish and suggest their contamination levels are a good site-specific indicator of contamination in the Galveston Bay area. The spatial gradient of PCB contamination in Galveston Bay was evident in all three matrices studied and was observed effectively using Gulf killifish contamination as an environmentally relevant bioindicator of localized contamination in this environment.

Microbial Community Structure and Function Decoupling Across a Phosphorus Gradient in Streams.

Phosphorus (P) is a key biological element with important and unique biogeochemical cycling in natural ecosystems. Anthropogenic phosphorus inputs have been shown to greatly affect natural ecosystems, and this has been shown to be especially true of freshwater systems. While the importance of microbial communities in the P cycle is widely accepted, the role, composition, and relationship to P of these communities in freshwater systems still hold many secrets. Here, we investigated combined bacterial and archaeal communities utilizing 16S ribosomal RNA (rRNA) gene sequencing and computationally predicted functional metagenomes (PFMs) in 25 streams representing a strong P gradient. We discovered that 16S rRNA community structure and PFMs demonstrate a degree of decoupling between structure and function in the system. While we found that total phosphorus (TP) was correlated to the structure and functional capability of bacterial and archaeal communities in the system, turbidity had a stronger, but largely independent, correlation. At TP levels of approximately 55 µg/L, we see sharp differences in the abundance of numerous ecologically important taxa related to vegetation, agriculture, sediment, and other ecosystem inhabitants.

Freshwater eutrophication drives sharp reductions in temporal beta diversity.

Eutrophication has become one of the most widespread anthropogenic forces impacting freshwater biological diversity. One potentially important mechanism driving biodiversity changes in response to eutrophication is the alteration of seasonal patterns of succession, particularly among species with short, synchronous, life cycles. We tested the hypothesis that eutrophication reduces seasonally driven variation in species assemblages by focusing on an understudied aspect of biodiversity: temporal beta diversity (ßt ). We estimated the effect of eutrophication on ßt by sampling benthic macroinvertebrate assemblages bimonthly for two years across 35 streams spanning a steep gradient of total phosphorus (P) and benthic algal biomass (as chlorophyll a [chl a]). Two widely used metrics of ß diversity both declined sharply in response to increasing P and chl a, regardless of covariates. The most parsimonious explanatory model for ßt included an interaction between P and macroinvertebrate biomass, which revealed that ßt was lower when macroinvertebrate biomass was relatively high. Macroinvertebrate biomass explained a greater amount of deviance in ßt at lower to moderate concentrations of P, providing additional explanatory power where P concentration alone was unable to fully explain declines in ßt . Chl a explained similar amounts of deviance in ßt in comparison to the best P model, but only when temperature variability, which was positively related to ßt , also was included in the model. Declines in ßt suggest that nutrient enrichment decreases the competitive advantage that specialists gain by occupying particular temporal niches, which leads to assemblages dominated by generalists that exhibit little seasonal turnover. The collapse of seasonal variation in assemblage composition we observed in our study suggests that treating dynamic communities as static assemblages is a simplification that may fail to detect the full impact of anthropogenic stressors. Our results show that eutrophication leads to more temporally homogenous communities and therefore degrades a fundamental facet of biodiversity.

Effects of pulsed atrazine exposures on autotrophic community structure, biomass, and production in field-based stream mesocosms.

The authors performed a multiple-pulsed atrazine experiment to measure responses of autotrophic endpoints in outdoor stream mesocosms. The experiment was designed to synthetically simulate worst-case atrazine chemographs from streams in agricultural catchments to achieve 60-d mean concentrations of 0 µg/L (control), 10 µg/L, 20 µg/L, and 30 µg/L. The authors dosed triplicate streams with pulses of 0 µg/L, 50 µg/L, 100 µg/L, and 150 µg/L atrazine for 4 d, followed by 7 d without dosing. This 11-d cycle occurred 3 times, followed by a recovery (untreated) period from day 34 to day 60. Mean ± standard error 60-d atrazine concentrations were 0.07 ± 0.03 µg/L, 10.7 ± 0.05 µg/L, 20.9 ± 0.24 µg/L, and 31.0 ± 0.17 µg/L for the control, 10-µg/L, 20-µg/L, and 30-µg/L treatments, respectively. Multivariate analyses revealed that periphyton and phytoplankton community structure did not differ among treatments on any day of the experiment, including during the atrazine pulses. Control periphyton biomass in riffles was higher immediately following the peak of the first atrazine pulse and remained slightly higher than some of the atrazine treatments on most days through the peak of the last pulse. However, periphyton biomass was not different among treatments at the end of the present study. Phytoplankton biomass was not affected by atrazine. Metaphyton biomass in pools was higher in the controls near the midpoint of the present study and remained higher on most days for the remainder of the study. Ceratophyllum demersum, a submersed macrophyte, biomass was higher in controls than in 20-µg/L and 30-µg/L treatments before pulse 3 but was not different subsequent to pulse 3 through the end of the present study. Maximum daily dissolved oxygen (DO, percentage of saturation) declined during each pulse in approximate proportion to magnitude of dose but rapidly converged among treatments after the third pulse. However, DO increased in controls relative to all atrazine treatments during the last 17 d of the experiment, likely a result of metaphyton cover in the pools. Finally, atrazine significantly limited uptake of PO4(3-) and uptake and/or denitrification of NO3(-) but only during pulses; percentage of dose removed from the water column was >85% for P and >95% for N after pulse 3 through the end of the present study. Collectively, only DO and metaphyton biomass differed at the end of the present study and only slightly. Some other endpoints were affected but only during pulses, if at all. The high levels of primary production and accumulation of algal biomass in all streams suggest that effects of pulses of atrazine at the concentrations used in the present study appear transient and likely do not represent ecologically significant adverse outcomes to periphyton, phytoplankton, and aquatic macrophytes, particularly in agricultural streams subjected to high nutrient loads.

Influence of drought and total phosphorus on diel pH in wadeable streams: implications for ecological risk assessment of ionizable contaminants.

Climatological influences on site-specific ecohydrology are particularly germane in semiarid regions where instream flows are strongly influenced by effluent discharges. Because many traditional and emerging aquatic contaminants, such as pharmaceuticals, are ionizable, we examined diel surface water pH patterns (i.e., change in pH over a 24-h period) at 23 wadeable streams in central Texas, USA, representing a gradient of nutrient enrichment during consecutive summers of 2006 and 2007. The years of our study were characterized by decidedly different instream flows, which likely affected production:respiration dynamics and led to distinctions in diel pH patterns between 2006 and 2007. Site-specific ambient water quality criteria for NH(3) and the aquatic toxicity of the model weak base pharmaceutical sertraline were predicted using continuous water quality monitoring data from the sites. Drought conditions of 2006 significantly increased (p<0.05) diel pH changes compared to high instream flows of 2007,and the magnitude of diel pH variability was most pronounced at nutrient-enriched sites in 2006. Differences in diel pH change patterns between 2006 and 2007 affected predictions of the environmental fate and effects for model weak base pharmaceuticals and NH(3). Overall, site-specific diel pH was more variable at some sites than the difference in mean surface water pH between the 2 summers. Diel pH variability affected regulatory criteria, because 20% of the study sites in 2006 experienced greater than 5-fold differences in National Ambient Water Quality Criteria for NH(3) over 24-h periods. Our study emphasizes the potential uncertainty that diel pH variability may introduce in site-specific assessments and provides recommendations for environmental assessment of ionizable contaminants.

Exploring Lemna gibba thresholds to nutrient and chemical stressors: differential effects of triclosan on internal stoichiometry and nitrate uptake across a nitrogen:phosphorus gradient.

Nutrient enrichment often co-occurs with chemical stressors in aquatic ecosystems, but the impacts of these multiple stressors across nutrient gradients is poorly understood and not typically addressed in ecotoxicity studies of lower trophic level models. Moreover, laboratory assays performed to determine threshold responses of aquatic macrophytes to contaminants typically use growth and morphometric endpoints to establish threshold effects and seldom report other important functional responses of lower trophic levels. Using the aquatic macrophyte Lemna gibba, we examined influences of varying nitrogen (N) and phosphorus (P) levels in combination with triclosan, a widely used antimicrobial agent in consumer care products, on internal carbon (C):N:P and NO(3) (-) uptake kinetics. Triclosan modulated L. gibba tissue N and P content, and these stoichiometric responses for P-limited plants to triclosan exposure were more sensitive than growth endpoints employed in standardized phytotoxicity assays. Nitrate uptake capacities were also differentially inhibited by triclosan exposure according to external nutrient levels. Uptake rates for plants cultured and exposed under saturating N-levels were inhibited by more than threefold compared with N-limited plants. The results suggest that stoichiometric and nutrient uptake responses to chemical stressors provide useful information regarding adverse ecological thresholds not defined in standardized phytotoxicity assays with aquatic macrophytes. Our findings further indicate that site-specific impacts of chemicals associated with the wide ambient ranges of N and P typical of surface waters may be anticipated in lower trophic levels. Future studies should examine adverse effects of other stressors to these ecologically relevant endpoints, which may be useful in environmental assessment and management.

Influence of nitrogen and phosphorus concentrations and ratios on Lemna gibba growth responses to triclosan in laboratory and stream mesocosm experiments.

The effects of co-occurring nutrient and contaminant stressors are very likely to interact in aquatic systems, particularly at the level of primary producers. Site-specific nitrogen (N) and phosphorus (P) concentrations are often much lower and differ in relative availability than those used in nutrient-saturated laboratory assays for aquatic plants, which can introduce uncertainty in prospective ecological hazard and risk assessments. Because triclosan, an antimicrobial agent included in personal care products, potentially presents high relative risk among antimicrobial agents to aquatic plants and algae, we performed laboratory experiments with the model aquatic macrophyte Lemna gibba across a gradient of environmentally relevant N:P levels with and without triclosan co-exposure. Frond numbers (7 d) were significantly higher in N:P treatments of 16 and 23 but were lower in N:P of 937 and 2,500 treatments relative to standardized control media (N:P=3). When triclosan co-exposure occurred at high nutrient concentrations, frond number median effective concentration values at N:P 0.75, 3, and 16 were more than twofold lower than triclosan median effective concentration values in low nutrient media N:P ratios. However, a triclosan median effective concentration for frond number was twofold lower at N:P of 2,500 than at other N:P ratios in low concentration media. Influences of P enrichment on triclosan toxicity to L. gibba were further explored during a 14-d outdoor experimental stream mesocosm study. Effects of 2.6 and 20.8 microg L(-1) triclosan on L. gibba growth rates were more pronounced with increasing P treatment levels, which was generally consistent with our laboratory observations. Findings from these laboratory and field studies indicate that site-specific nutrient concentrations and ratios should be considered during assessments of primary producer responses to chemical stressors.