Impact of harmful algal bloom severity on bacterial communities in a full-scale biological filtration system for drinking water treatment

The occurrence of harmful algal blooms (HABs) in freshwater environments has been expanded worldwide with growing frequency and severity. HABs can pose a threat to public water supplies, raising concerns about safety of treated water. Many studies have provided valuable information about the impacts of HABs and management strategies on the early-stage treatment processes (e.g., pre-oxidation and coagulation/flocculation) in conventional drinking water treatment plants (DWTPs). However, the potential effect of HAB-impacted water in the granular media filtration has not been well studied. Biologically-active filters (BAFs), which are used in drinking water treatment and rely largely on bacterial community interactions, have not been examined during HABs in full-scale DWTPs. In this study, we assessed the bacterial community structure of BAFs, functional profiles, assembly processes, and bio-interactions in the community during both severe and mild HABs. Our findings indicate that bacterial diversity in BAFs significantly decreases during severe HABs due to the predominance of bloom-associated bacteria (e.g., Spingopyxis, Porphyrobacter, and Sphingomonas). The excitation-emission matrix combined with parallel factor analysis (EEM-PARAFAC) confirmed that filter influent affected by the severe HAB contained a higher portion of protein-like substances than filter influent samples during a mild bloom. In addition, BAF community functions showed increases in metabolisms associated with intracellular algal organic matter (AOM), such as lipids and amino acids, during severe HABs. Further ecological process and network analyses revealed that severe HAB, accompanied by the abundance of bloom-associated taxa and increased nutrient availability, led to not only strong stochastic processes in the assembly process, but also a bacterial community with lower complexity in BAFs. Overall, this study provides deeper insights into BAF bacterial community structure, function, and assembly in response to HABs.

The Lake Erie HABs Grab: A binational collaboration to characterize the western basin cyanobacterial harmful algal blooms at an unprecedented high-resolution spatial scale

Monitoring of cyanobacterial bloom biomass in large lakes at high resolution is made possible by remote sensing. However, monitoring cyanobacterial toxins is only feasible with grab samples, which, with only sporadic sampling, results in uncertainties in the spatial distribution of toxins. To address this issue, we conducted two intensive "HABs Grabs" of microcystin (MC)-producing Microcystis blooms in the western basin of Lake Erie. These were one-day sampling events during August of 2018 and 2019 in which 100 and 172 grab samples were collected, respectively, within a six-hour window covering up to 2,270 km² and analyzed using consistent methods to estimate the total mass of MC. The samples were analyzed for 57 parameters, including toxins, nutrients, chlorophyll, and genomics. There were an estimated 11,513 kg and 30,691 kg of MCs in the western basin during the 2018 and 2019 HABs Grabs, respectively. The bloom boundary poses substantial issues for spatial assessments because MC concentration varied by nearly two orders of magnitude over very short distances. The MC to chlorophyll ratio (MC:chl) varied by a factor up to 5.3 throughout the basin, which creates challenges for using MC:chl to predict MC concentrations. Many of the biomass metrics strongly correlated (r > 0.70) with each other except chlorophyll fluorescence and phycocyanin concentration. While MC and chlorophyll correlated well with total phosphorus and nitrogen concentrations, MC:chl correlated with dissolved inorganic nitrogen. More frequent MC data collection can overcome these issues, and models need to account for the MC:chl spatial heterogeneity when forecasting MCs.

Highlighting the promise of qPCR-based environmental monitoring: response of the ribosomal RNA:DNA ratio of calanoid copepods to toxic cyanobacteria

Calanoid copepods are integral to aquatic food webs and may drive the bioaccumulation of toxins and heavy metals, spread of infectious diseases, and occurrence of toxic cyanobacterial harmful algal blooms (HABs) in freshwater aquatic systems. However, interrelationships between copepod and cyanobacterial population dynamics and ecophysiology remain unclear. Insights into these relationships are important to aquatic resource management, as they may help guide mitigation efforts. We developed a calanoid copepod qPCR assay to investigate how copepod abundance and physiological status relate to the abundance of cyanobacteria and the concentration of total microcystin in a HAB-prone freshwater multi-use eutrophic lake. Through in silico and in vitro validation of primers and analyses of time series, we demonstrate that our assay can be used as a reliable tool for environmental monitoring. Importantly, copepod RNA:DNA ratios on and shortly after the day when microcystin concentration was at its highest within the lake were not significantly lower (or higher) than before or after this period, suggesting that copepods may have been tolerant of microcystin levels observed and capable of perpetuating bloom events by consuming competitors of toxic cyanobacteria.

Ecological insights into assembly processes and network structures of bacterial biofilms in full-scale biologically active carbon filters under ozone implementation

To address the adverse effects of harmful algal blooms, there are increased demands over the implementation of ozone coupled with biologically active carbon (BAC) filters in the drinking water treatment plants. Although the microbial biofilms are vital elements to support the proper performance of BAC filters, except for taxonomic affiliations, little is known about the assembly mechanisms of microbial communities in the full-scale BAC filters. This study aimed to examine how the assembly processes and their associated factors (e.g., influent characteristics, biological interactions) drive the temporal dynamics of bacterial communities in full-scale BAC filters, which underwent ozone implementation (five consecutive seasons from 2017 to 2018). The results revealed that along with the increase of bacterial taxonomic richness and evenness, stochastic processes became more crucial to determine the bacterial community assembly in the summer and autumn after ozone implementation (relative contribution: 61.23% and 83.75%, respectively). Moreover, their corresponding networks possessed simple network structures with lower modularity than other seasons, which implied lesser biological interactions among bacterial populations. The correlation between taxonomic and predicted functional diversities using functional redundancy index indicated that relatively high levels of bacterial functional redundancy (>0.83) were generally present in BAC filters. However, compared to other seasons, significantly higher degrees of functional redundancy existed in the summer and autumn after ozone implementation (0.85 ± 0.01 and 0.86 ± 0.01, respectively). Overall, this work improves our understanding of the microbial ecology of full-scale BAC filters by providing a conceptual framework that characterizes bacterial biofilm assembly processes relevant to performance optimization of full-scale BAC filters.

Impact of algal organic matter on the performance, cyanotoxin removal, and biofilms of biologically-active filtration systems

The occurrence of harmful algal blooms dominated by toxic cyanobacteria has induced continuous loadings of algal organic matter (AOM) and toxins in drinking water treatment plants. However, the impact of AOM on the active biofilms and microbial community structures of biologically-active filtration (BAF), which directly affects the contaminant removal, is not well understood. In this study, we systematically examined the effects of AOM on BAF performance and bacterial biofilm formation over 240 days, tracing the removal of specific AOM components, a cyanotoxin [microcystin-LR (MC-LR)], and microbial community responses. The component analysis (excitation and emission matrix analysis) results for AOM revealed that terrestrial humic-like substances showed the highest removal among all the identified components and were strongly correlated to MC-LR removal. In addition, reduced empty bed contact time and deactivation of biofilms significantly decreased BAF performances for both AOM and MC-LR. The active biofilm, bacterial community structure, and mlrA gene (involved in microcystin degradation) abundance demonstrated that bacterial biofilm composition responded to AOM and MC-LR, in which Rhodocyclaceae, Saprospiraceae, and Comamonadaceae were dominant. In addition, MC-LR biodegradation appeared to be more active at the top than at the bottom layer in BAF. Overall, this study provides deeper insights into the role of biofilms and filter operation on the fate of AOM and MC-LR in BAF.

Assessing the chemical compositions and disinfection byproduct formation of biofilms: Application of fluorescence excitation-emission spectroscopy coupled with parallel factor analysis

There are increased concerns over the contributions of biofilms to disinfection byproduct (DBP) formation in engineered water systems (EWS). However, monitoring the biomolecular characteristics of biofilms to understand their impacts on DBP formation has been a great challenge as it requires complex analytical techniques. This study aimed to examine the applicability of fluorescence excitation-emission matrices (EEMs) coupled with parallel factor analysis (PARAFAC) to assess the chemical compositions and DBP formation of biofilms. Biofilms were collected from reactors grown on R2A media, as well as two drinking water-related organic substrates such as humic substances and algal organic matter. The chemical composition and formation of carbonaceous and nitrogenous DBPs of biofilms were continuously monitored every 21 days for 168 days and correlated with the derived EEM-PARAFAC components. Results indicated that all biofilm samples comprised mostly of protein-like components (∼90%), and to a lesser extent, humic-like components (∼10%). Strong correlations were generally found between tryptophan-like substances and the studied DBP formation (R2min ≥ 0.76, P < 0.05), indicating that they play a major role in producing biofilm-derived DBPs upon chlorination. Moreover, significant discrepancies between the chemical compositions and DBP formation of biofilms and their corresponding feed solutions were observed, likely due to biotransformation and biosorption processes. Overall, this work highlights that EEM-PARAFAC analysis is a promising tool to monitor the biomolecular characteristics of biofilm components and to predict the subsequent DBP formation in optimizing disinfection protocols for EWS.

Composition and Drivers of Gut Microbial Communities in Arctic-Breeding Shorebirds

Gut microbiota can have important effects on host health, but explanatory factors and pathways that determine gut microbial composition can differ among host lineages. In mammals, host phylogeny is one of the main drivers of gut microbiota, a result of vertical transfer of microbiota during birth. In birds, it is less clear what the drivers might be, but both phylogeny and environmental factors may play a role. We investigated host and environmental factors that underlie variation in gut microbiota composition in eight species of migratory shorebirds. We characterized bacterial communities from 375 fecal samples collected from adults of eight shorebird species captured at a network of nine breeding sites in the Arctic and sub-Arctic ecoregions of North America, by sequencing the V4 region of the bacterial 16S ribosomal RNA gene. Firmicutes (55.4%), Proteobacteria (13.8%), Fusobacteria (10.2%), and Bacteroidetes (8.1%) dominated the gut microbiota of adult shorebirds. Breeding location was the main driver of variation in gut microbiota of breeding shorebirds (R² = 11.6%), followed by shorebird host species (R² = 1.8%), and sampling year (R² = 0.9%), but most variation remained unexplained. Site variation resulted from differences in the core bacterial taxa, whereas rare, low-abundance bacteria drove host species variation. Our study is the first to highlight a greater importance of local environment than phylogeny as a driver of gut microbiota composition in wild, migratory birds under natural conditions.

Dynamics of the physiochemical and community structures of biofilms under the influence of algal organic matter and humic substances

Increased loading of algal organic matter (AOM) during harmful algal blooms not only burdens water treatment processes but also challenges safe drinking water delivery. While organic constituents promote biofilm growth in drinking water distribution systems (DWDS), the effects of AOM on biofilm formation in DWDS are not well understood. Herein, three parallel biofilm reactors were used to assess and compare how treated AOM- and humic substance (HS)-impacted bulk water, and R2A medium (a control) affect biofilm development for 168 days. The 16S rRNA gene sequencing analysis revealed that the bacterial communities in biofilms were clustered with the organic matter types in bulk water, where Family Comamonadaceae was the most dominant but showed different temporal dynamics depending on the organic matter characteristics in bulk water. Higher diversity was observed in the biofilms grown in AOM-impacted bulk water (BFAOM) than biofilms grown in HS-impacted (BFHS) and R2A-impacted bulk water (BFR2A) as the biofilms matured. In addition, some taxa (e.g., Rhodobacteraceae and Sphingomonadaceae) were enriched in BFAOM compared to BFHS and BFR2A. The biofilm image analysis results indicated that compared to BFHS, BFAOM and BFR2A had relatively thinner and heterogeneous physical structures with lower amounts of cell biomass, extracellular polymeric substances (EPS), and higher EPS protein/polysaccharide ratios. Overall, this study revealed how AOM- and HS-impacted bulk water shape the physiochemical and community structures of biofilms, which can provide insights into assessing biofilm-associated risks and optimizing disinfection practices for biofilm control in DWDS.

Categorical performance characteristics of method ISO 7899-2 and indicator value of intestinal enterococci for bathing water quality monitoring

Intestinal enterococci indicate the fecal contamination of bathing waters. This study defines the performance characteristics of the reference method ISO 7899-2:2000 with water samples collected from inland and coastal bathing areas in Finland. From a total of 341 bacterial isolates grown on Slanetz and Bartley medium, 63.6% were confirmed as intestinal enterococci on bile aesculin agar. The partial 16S rRNA gene sequences showed that Enterococcus faecium and Enterococcus faecalis clades accounted for 93.1% of the confirmed isolates. The range of the false positive and false negative rate of the ISO 7899-2 was 0.0-18.5% and 5.6-57.1%, respectively, being affected by the presumptive colony count on the membrane. The analysis of multiple sample volumes is proposed to reach 10-100 colonies per membrane when 47 mm diameter membranes are used to prevent overestimation of low counts and underestimation of the high counts.

Microbial degradation of Cold Lake Blend and Western Canadian select dilbits by freshwater enrichments

Treatability experiments were conducted to determine the biodegradation of diluted bitumen (dilbit) at 5 and 25 °C for 72 and 60 days, respectively. Microbial consortia obtained from the Kalamazoo River Enbridge Energy spill site were enriched on dilbit at both 5 (cryo) and 25 (meso) ºC. On every sampling day, triplicates were sacrificed and residual hydrocarbon concentrations (alkanes and polycyclic aromatic hydrocarbons) were determined by GCMS/MS. The composition and relative abundance of different bacterial groups were identified by 16S rRNA gene sequencing analysis. While some physicochemical differences were observed between the two dilbits, their biodegradation profiles were similar. The rates and extent of degradation were greater at 25 °C. Both consortia metabolized 99.9% of alkanes; however, the meso consortium was more effective at removing aromatics than the cryo consortium (97.5 vs 70%). Known hydrocarbon-degrading bacteria were present in both consortia (Pseudomonas, Rhodococcus, Hydrogenophaga, Parvibaculum, Arthrobacter, Acidovorax), although their relative abundances depended on the temperatures at which they were enriched. Regardless of the dilbit type, the microbial community structure significantly changed as a response to the diminishing hydrocarbon load. Our results demonstrate that dilbit can be effectively degraded by autochthonous microbial consortia from sites with recent exposure to dilbit contamination.

Microbial activity influences electrical conductivity of biofilm anode

This study assessed the conductivity of a Geobacter-enriched biofilm anode in a microbial electrochemical cell (MxC) equipped with two gold anodes (25 mM acetate medium), as different proton gradients were built throughout the biofilm. There was no pH gradient across the biofilm anode at 100 mM phosphate buffer (current density 2.38 A/m2) and biofilm conductivity (Kbio) was as high as 0.87 mS/cm. In comparison, an inner biofilm became acidic at 2.5 mM phosphate buffer in which dead cells were accumulated at ∼80 μm of the inner biofilm anode. At this low phosphate buffer, Kbio significantly decreased by 0.27 mS/cm, together with declined current density of 0.64 A/m2. This work demonstrates that biofilm conductivity depends on the composition of live and dead cells in the conductive biofilm anode.

Whole-Community Metagenomics in Two Different Anammox Configurations: Process Performance and Community Structure

Anaerobic ammonia oxidation (anammox) combined with partial nitritation (PN) is an innovative treatment process for energy-efficient nitrogen removal from wastewater. In this study, we used genome-based metagenomics to investigate the overall community structure and anammox species enriched in suspended growth (SGR) and attached growth packed-bed (AGR) anammox reactors after 220 days of operation. Both reactors removed more than 85% of the total inorganic nitrogen. Metagenomic binning and phylogenetic analysis revealed that two anammox population genomes, affiliated with the genus Candidatus Brocadia, were differentially abundant between the SGR and AGR. Both of the genomes shared an average nucleotide identify of 83%, suggesting the presence of two different species enriched in both of the reactors. Metabolic reconstruction of both population genomes revealed key aspects of their metabolism in comparison to known anammox species. The community composition of both the reactors was also investigated to identify the presence of flanking community members. Metagenomics and 16S rRNA gene amplicon sequencing revealed the dominant flanking community members in both reactors were affiliated with the phyla Anaerolinea, Ignavibacteria, and Proteobacteria. Findings from this research adds two new species, Ca. Brocadia sp. 1 and Ca. Brocadia sp. 2, to the genus Ca. Brocadia and sheds light on their metabolism in engineered ecosystems.

Analysis of human mitochondrial DNA sequences from fecally polluted environmental waters as a tool to study population diversity

Mitochondrial signature sequences have frequently been used to study human population diversity around the world. Traditionally, this requires obtaining samples directly from individuals which is cumbersome, time consuming and limited to the number of individuals that participated in these types of surveys. Here, we used environmental DNA extracts to determine the presence and sequence variability of human mitochondrial sequences as a means to study the diversity of populations inhabiting in areas nearby a tropical watershed impacted with human fecal pollution. We used high-throughput sequencing (Illumina) and barcoding to obtain thousands of sequences from the mitochondrial hypervariable region 2 (HVR2) and determined the different haplotypes present in 10 different water samples. Sequence analyses indicated a total of 19 distinct variants with frequency greater than 5%. The HVR2 sequences were associated with haplogroups of West Eurasian (57.6%), Sub-Saharan African (23.9%), and American Indian (11%) ancestry. This was in relative accordance with population census data from the watershed sites. The results from this study demonstrates the potential value of mitochondrial sequence data retrieved from fecally impacted environmental waters to study the population diversity of local municipalities. This environmental DNA approach may also have other public health implications such as tracking background levels of human mitochondrial genes associated with diseases. It may be possible to expand this approach to other animal species inhabiting or using natural water systems.

High Biofilm Conductivity Maintained Despite Anode Potential Changes in a Geobacter-Enriched Biofilm

This study systematically assessed intracellular electron transfer (IET) and extracellular electron transfer (EET) kinetics with respect to anode potential (E_anode ) in a mixed-culture biofilm anode enriched with Geobacter spp. High biofilm conductivity (0.96-1.24 mS cm^-1 ) was maintained during E_anode changes from -0.2 to +0.2 V versus the standard hydrogen electrode (SHE), although the steady-state current density significantly decreased from 2.05 to 0.35 A m^-2 in a microbial electrochemical cell. Substantial increase of the Treponema population was observed in the biofilm anode at E_anode =+0.2 V, which reduced intracellular electron-transfer kinetics associated with the maximum specific substrate-utilization rate by a factor of ten. This result suggests that fast EET kinetics can be maintained under dynamic E_anode conditions in a highly conductive biofilm anode as a result of shift of main EET players in the biofilm anode, although E_anode changes can influence IET kinetics.

The Roles of Biofilm Conductivity and Donor Substrate Kinetics in a Mixed-Culture Biofilm Anode

We experimentally assessed the kinetics and thermodynamics of electron transfer (ET) from the donor substrate (acetate) to the anode for a mixed-culture biofilm anode. We interpreted the results with a modified biofilm-conduction model consisting of three ET steps in series: (1) intracellular ET, (2) non-Ohmic extracellular ET (EET) from an outer membrane protein to an extracellular cofactor (EC), and (3) ET from the EC to the anode by Ohmic-conduction in the biofilm matrix. The steady-state current density was 0.82 ± 0.03 A/m² in a miniature microbial electrochemical cell operated at fixed anode potential of -0.15 V versus the standard hydrogen electrode. Illumina 16S-rDNA and -rRNA sequences showed that the Geobacter genus was less than 30% of the community of the biofilm anode. Biofilm conductivity was high at 2.44 ± 0.42 mS/cm, indicating that the maximum current density could be as high as 270 A/m² if only Ohmic-conduction EET was limiting. Due to the high biofilm conductivity, the maximum energy loss for Ohmic-conduction EET was negligible, 0.085 mV. The energy loss in the second ET step also was small, only 20 mV, and the potential for the EC involved in the second ET was -0.15 V, a value documenting that >99% of the EC was in the oxidized state. Monod kinetics for utilization of acetate were relatively slow, and at least 87% of the energy loss was in the intracellular step. Thus, intracellular ET was the main kinetic and thermodynamic bottleneck to ET from donor substrate to the anode for a highly conductive biofilm.

Ohmic resistance affects microbial community and electrochemical kinetics in a multi-anode microbial electrochemical cell

Multi-anode microbial electrochemical cells (MxCs) are considered as one of the most promising configurations for scale-up of MxCs, but understanding of anode kinetics in multiple anodes is limited in the MxCs. In this study we assessed microbial community and electrochemical kinetic parameters for biofilms on individual anodes in a multi-anode MxC to better comprehend anode fundamentals. Microbial community analysis targeting 16S rRNA Illumina sequencing showed that Geobacter genus was abundant (87%) only on the biofilm anode closest to a reference electrode (low ohmic energy loss) in which current density was the highest among three anodes. In comparison, Geobacter populations were less than 1% for biofilms on other two anodes distant from the reference electrode (high ohmic energy loss), generating small current density. Half-saturation anode potential (E_KA) was the lowest at -0.251 to -0.242 V (vs. standard hydrogen electrode) for the closest biofilm anode to the reference electrode, while E_KA was as high as -0.134 V for the farthest anode. Our study proves that electric potential of individual anodes changed by ohmic energy loss shifts biofilm communities on individual anodes and consequently influences electron transfer kinetics on each anode in the multi-anode MxC.

Use of functional gene expression and respirometry to study wastewater nitrification activity after exposure to low doses of copper

Autotrophic nitrification in biological nitrogen removal systems has been shown to be sensitive to the presence of heavy metals in wastewater treatment plants. Using transcriptase-quantitative polymerase chain reaction (RT-qPCR) data, we examined the effect of copper on the relative expression of functional genes (i.e., amoA, hao, nirK, and norB) involved in redox nitrogen transformation in batch enrichment cultures obtained from a nitrifying bioreactor operated as a continuous reactor (24-h hydraulic retention time). 16S ribosomal RNA (rRNA) gene next-generation sequencing showed that Nitrosomonas-like populations represented 60-70% of the bacterial community, while other nitrifiers represented <5%. We observed a strong correspondence between the relative expression of amoA and hao and ammonia removal in the bioreactor. There were no considerable changes in the transcript levels of amoA, hao, nirK, and norB for nitrifying samples exposed to copper dosages ranging from 0.01 to 10 mg/L for a period of 12 h. Similar results were obtained when ammonia oxidation activity was measured via specific oxygen uptake rate (sOUR). The lack of nitrification inhibition by copper at doses lower than 10 mg/L may be attributed to the role of copper as cofactor for ammonia monooxygenase or to the sub-inhibitory concentrations of copper used in this study. Overall, these results demonstrate the use of molecular methods combined with conventional respirometry assays to better understand the response of wastewater nitrifying systems to the presence of copper.

Effects of Cr(III) and Cr(VI) on nitrification inhibition as determined by SOUR, function-specific gene expression and 16S rRNA sequence analysis of wastewater nitrifying enrichments

The effect of Cr(III) and Cr(VI) on nitrification was examined with samples from nitrifying enrichment cultures using three different approaches: by measuring substrate (ammonia) specific oxygen uptake rates (SOUR), by using RT-qPCR to quantify the transcripts of functional genes involved in nitrification, and by analysis of 16S rRNA sequences to determine changes in structure and activity of the microbial communities. The nitrifying bioreactor was operated as a continuous reactor with a 24 h hydraulic retention time. The samples were exposed in batch vessels to Cr(III) (10-300 mg/L) and Cr(VI) (1-30 mg/L) for a period of 12 h. There was considerable decrease in SOUR with increasing dosages for both Cr(III) and Cr(VI), however Cr(VI) was more inhibitory than Cr(III). Based on the RT-qPCR data, there was reduction in the transcript levels of amoA and hao for increasing Cr(III) dosage, which corresponded well with the ammonia oxidation activity measured via SOUR. For Cr(VI) exposure, there was comparatively little reduction in amoA expression while hao expression decreased for 1-3 mg/L Cr(VI) and increased at 30 mg/L Cr(VI). While Nitrosomonas spp. were the dominant bacteria in the bioreactor, based on 16S rRNA sequencing, there was a considerable reduction in Nitrosomonas activity upon exposure to 300 mg/L Cr(III). In contrast, a relatively small reduction in activity was observed at 30 mg/L Cr(VI) loading. Our data that suggest that both Cr(III) and Cr(VI) were inhibitory to nitrification at concentrations near the high end of industrial effluent concentrations.

Impact of Heavy Metals on Transcriptional and Physiological Activity of Nitrifying Bacteria

Heavy metals can inhibit nitrification, a key process for nitrogen removal in wastewater treatment. The transcriptional responses of amoA, hao, nirK, and norB were measured in conjunction with specific oxygen uptake rate (sOUR) for nitrifying enrichment cultures exposed to different metals (Ni(II), Zn(II), Cd(II), and Pb(II)). There was significant decrease in sOUR with increasing concentrations for Ni(II) (0.03-3 mg/L), Zn(II) (0.1-10 mg/L), and Cd(II) (0.03-1 mg/L) (p < 0.05). However, no considerable changes in sOUR were observed with Pb(II) (1-100 mg/L), except at a dosage of 1000 mg/L causing 84% inhibition. Based on RT-qPCR data, the transcript levels of amoA and hao decreased when exposed to Ni(II) dosages. Slight up-regulation of amoA, hao, and nirK (0.5-1.5-fold) occurred after exposure to 0.3-3 mg/L Zn(II), although their expression decreased for 10 mg/L Zn(II). With the exception of 1000 mg/L Pb(II), stimulation of all genes occurred on Cd(II) and Pb(II) exposure. While overall the results show that RNA-based function-specific assays can be used as potential surrogates for measuring nitrification activity, the degree of inhibition inferred from sOUR and gene transcription is different. We suggest that variations in transcription of functional genes may supplement sOUR based assays as early warning indicators of upsets in nitrification.

Metabolic and genomic analysis elucidates strain-level variation in Microbacterium spp. isolated from chromate contaminated sediment

Hexavalent chromium [Cr(VI)] is a soluble carcinogen that has caused widespread contamination of soil and water in many industrial nations. Bacteria have the potential to aid remediation as certain strains can catalyze the reduction of Cr(VI) to insoluble and less toxic Cr(III). Here, we examine Cr(VI) reducing Microbacterium spp. (Cr-K1W, Cr-K20, Cr-K29, and Cr-K32) isolated from contaminated sediment (Seymore, Indiana) and show varying chromate responses despite the isolates' phylogenetic similarity (i.e., identical 16S rRNA gene sequences). Detailed analysis identified differences based on genomic metabolic potential, growth and general metabolic capabilities, and capacity to resist and reduce Cr(VI). Taken together, the discrepancies between the isolates demonstrate the complexity inter-strain variation can have on microbial physiology and related biogeochemical processes.

New perspectives on microbial community distortion after whole-genome amplification

Whole-genome amplification (WGA) has become an important tool to explore the genomic information of microorganisms in an environmental sample with limited biomass, however potential selective biases during the amplification processes are poorly understood. Here, we describe the effects of WGA on 31 different microbial communities from five biotopes that also included low-biomass samples from drinking water and groundwater. Our findings provide evidence that microbiome segregation by biotope was possible despite WGA treatment. Nevertheless, samples from different biotopes revealed different levels of distortion, with genomic GC content significantly correlated with WGA perturbation. Certain phylogenetic clades revealed a homogenous trend across various sample types, for instance Alpha- and Betaproteobacteria showed a decrease in their abundance after WGA treatment. On the other hand, Enterobacteriaceae, an important biomarker group for fecal contamination in groundwater and drinking water, were strongly affected by WGA treatment without a predictable pattern. These novel results describe the impact of WGA on low-biomass samples and may highlight issues to be aware of when designing future metagenomic studies that necessitate preceding WGA treatment.

Bacterial composition in a metropolitan drinking water distribution system utilizing different source waters

We investigated the bacterial composition of water samples from two service areas within a drinking water distribution system (DWDS), each associated with a different primary source of water (groundwater, GW; surface water, SW) and different treatment process. Community analysis based on 16S rRNA gene clone libraries indicated that Actinobacteria (Mycobacterium spp.) and α-Proteobacteria represented nearly 43 and 38% of the total sequences, respectively. Sequences closely related to Legionella, Pseudomonas, and Vibrio spp. were also identified. In spite of the high number of sequences (71%) shared in both areas, multivariable analysis revealed significant differences between the GW and SW areas. While the dominant phylotypes where not significantly contributing in the ordination of samples, the populations associated with the core of phylotypes (1-10% in each sample) significantly contributed to the differences between both service areas. Diversity indices indicate that the microbial community inhabiting the SW area is more diverse and contains more distantly related species coexisting with local assemblages as compared with the GW area. The bacterial community structure of SW and GW service areas were dissimilar, suggesting that their respective source water and/or water quality parameters shaped by the treatment processes may contribute to the differences in community structure observed.

Microbial community response to chlorine conversion in a chloraminated drinking water distribution system

Temporary conversion to chlorine (i.e., "chlorine burn") is a common approach to controlling nitrification in chloraminated drinking water distribution systems, yet its effectiveness and mode(s) of action are not fully understood. This study characterized occurrence of nitrifying populations before, during and after a chlorine burn at 46 sites in a chloraminated distribution system with varying pipe materials and levels of observed nitrification. Quantitative polymerase chain reaction analysis of gene markers present in nitrifying populations indicated higher frequency of detection of ammonia oxidizing bacteria (AOB) (72% of samples) relative to ammonia oxidizing archaea (AOA) (28% of samples). Nitrospira nitrite oxidizing bacteria (NOB) were detected at 45% of samples, while presence of Nitrobacter NOB could not be confirmed at any of the samples. During the chlorine burn, the numbers of AOA, AOB, and Nitrospira greatly reduced (i.e., 0.8-2.4 log). However, rapid and continued regrowth of AOB and Nitrospira were observed along with nitrite production in the bulk water within four months after the chlorine burn, and nitrification outbreaks appeared to worsen 6-12 months later, even after adopting a twice annual burn program. Although high throughput sequencing of 16S rRNA genes revealed a distinct community shift and higher diversity index during the chlorine burn, it steadily returned towards a condition more similar to pre-burn than burn stage. Significant factors associated with nitrifier and microbial community composition included water age and sampling location type, but not pipe material. Overall, these results indicate that there is limited long-term effect of chlorine burns on nitrifying populations and the broader microbial community.

Genome sequencing reveals the environmental origin of enterococci and potential biomarkers for water quality monitoring

Enterococci are common members of the gut microbiome and their ease of culturing has facilitated worldwide use as indicators of fecal pollution of waters. However, enterococci were recently shown to persist in environmental habitats, often in the absence of fecal input, potentially confounding water quality assays. Toward resolving this issue and providing a more complete picture of natural enterococci diversity, 11 isolates of Enterococcus faecalis recovered from freshwater watersheds (environmental) were sequenced and compared to 59 available enteric genomes. Phenotypically and phylogenetically the environmental E. faecalis were indistinguishable from their enteric counterparts. However, distinct environmental- and enteric-associated gene signatures, encoding mostly accessory nutrient utilization pathways, were detected among the variable genes. Specifically, a nickel uptake operon was over-represented in environmental genomes, while genes for utilization of sugars thought to be abundant in the gut such as xylose were over-represented in enteric genomes. The distribution and phylogeny of these identified signatures suggest that ancestors of E. faecalis resided in extra-enteric habitats, challenging the prevailing commensal view of enterococci ecology. Thus, habitat-associated gene content changes faster than core genome phylogeny and may include biomarkers for reliably detecting fecal contaminants for improved microbial water quality monitoring.

Syntrophic interactions between H2-scavenging and anode-respiring bacteria can improve current density in microbial electrochemical cells

High current density of 10.0-14.6A/m(2) and COD removal up to 96% were obtained in a microbial electrochemical cell (MEC) fed with digestate at hydraulic retention time (HRT) of 4d and 8d. Volatile fatty acids became undetectable in MEC effluent (HRT 8d), except for trivial acetate (4.16±1.86mgCOD/L). Accumulated methane only accounted for 3.42% of ΔCOD. Pyrosequencing analyses showed abundant fermenters (Kosmotoga spp.) and homoacetogens (Treponema spp.) in anolytes. In anode biofilm, propionate fermenters (Kosmotoga, and Syntrophobacter spp.), homoacetogens (Treponema spp.), and anode-respiring bacteria (ARB) (Geobacter spp. and Dysgonomonas spp.) were dominant. These results imply that syntrophic interactions among fermenters, homoacetogens and ARB would allow MECs to maintain high current density and coulombic efficiency.

Multi-laboratory evaluations of the performance of Catellicoccus marimammalium PCR assays developed to target gull fecal sources

Here we report results from a multi-laboratory (n = 11) evaluation of four different PCR methods targeting the 16S rRNA gene of Catellicoccus marimammalium originally developed to detect gull fecal contamination in coastal environments. The methods included a conventional end-point PCR method, a SYBR(®) Green qPCR method, and two TaqMan(®) qPCR methods. Different techniques for data normalization and analysis were tested. Data analysis methods had a pronounced impact on assay sensitivity and specificity calculations. Across-laboratory standardization of metrics including the lower limit of quantification (LLOQ), target detected but not quantifiable (DNQ), and target not detected (ND) significantly improved results compared to results submitted by individual laboratories prior to definition standardization. The unit of measure used for data normalization also had a pronounced effect on measured assay performance. Data normalization to DNA mass improved quantitative method performance as compared to enterococcus normalization. The MST methods tested here were originally designed for gulls but were found in this study to also detect feces from other birds, particularly feces composited from pigeons. Sequencing efforts showed that some pigeon feces from California contained sequences similar to C. marimammalium found in gull feces. These data suggest that the prevalence, geographic scope, and ecology of C. marimammalium in host birds other than gulls require further investigation. This study represents an important first step in the multi-laboratory assessment of these methods and highlights the need to broaden and standardize additional evaluations, including environmentally relevant target concentrations in ambient waters from diverse geographic regions.

Correlative assessment of fecal indicators using human mitochondrial DNA as a direct marker

Identifying the source of surface water fecal contamination is paramount to mitigating pollution and risk to human health. Fecal bacteria such as E. coli have been staple indicator organisms for over a century, however there remains uncertainty with E. coli-based metrics since these bacteria are abundant in the environment. The relationships between the presence of direct indicator of human waste (human mitochondrial DNA), human-specific Bacteroidales, and E. coli were studied for water samples taken from an urban creek system (Duck Creek Watershed, Cincinnati, OH) impacted by combined sewer overflows. Logistic regression analysis shows that human-specific Bacteroidales correlates much more closely to human mitochondrial DNA (R = 0.62) relative to E. coli (R = 0.33). We also examine the speciation of Bacteroidales within the Duck Creek Watershed using next-generation sequencing technology (Ion Torrent) and show the most numerous populations to be associated with sewage. Here we demonstrate that human-specific Bacteroidales closely follow the dynamics of human mitochondrial DNA concentration changes, indicating that these obligate anaerobes are more accurate than E. coli for fecal source tracking, lending further support to risk overestimation using coliforms, especially fecal coliforms and E. coli.

Establishment and early succession of bacterial communities in monochloramine-treated drinking water biofilms

Monochloramine is an increasingly used drinking water disinfectant and has been shown to increase nitrifying bacteria and mycobacteria in drinking waters. The potential successions and development of these bacteria were examined by 16S rRNA gene clone libraries generated from various biofilms within a water distribution system simulator. Biofilms were obtained from in-line and off-line devices using borosilicate glass beads, along with polycarbonate coupons from annular reactors incubated for up to 8 months in monochloramine-treated drinking water. No significant difference in community structures was observed between biofilm devices and coupon material; however, all biofilm communities that developed on different devices underwent similar successions over time. Early stages of biofilm formation were dominated by Serratia (29%), Cloacibacterium (23%), Diaphorobacter (16%), and Pseudomonas (7%), while Mycobacterium-like phylotypes were the most predominant populations (> 27%) in subsequent months. The development of members of the nontuberculous mycobacteria (NTM) after 3 months may impact individuals with predisposing conditions, while nitrifiers (related to Nitrospira moscoviensis and Nitrosospira multiformis) could impact water quality. Overall, 90% of the diversity in all the clone library samples was associated with the phyla Proteobacteria, Actinobacteria, and Bacteroidetes. These results provide an ecological insight into biofilm bacterial successions in monochloramine-treated drinking water.

Functional microbial diversity explains groundwater chemistry in a pristine aquifer

BACKGROUND

The diverse microbial populations that inhabit pristine aquifers are known to catalyze critical in situ biogeochemical reactions, yet little is known about how the structure and diversity of this subsurface community correlates with and impacts upon groundwater chemistry. Herein we examine 8,786 bacterial and 8,166 archaeal 16S rRNA gene sequences from an array of monitoring wells in the Mahomet aquifer of east-central Illinois. Using multivariate statistical analyses we provide a comparative analysis of the relationship between groundwater chemistry and the microbial communities attached to aquifer sediment along with those suspended in groundwater.

RESULTS

Statistical analyses of 16S rRNA gene sequences showed a clear distinction between attached and suspended communities; with iron-reducing bacteria far more abundant in attached samples than suspended, while archaeal clones related to groups associated with anaerobic methane oxidation and deep subsurface gold mines (ANME-2D and SAGMEG-1, respectively) distinguished the suspended community from the attached. Within the attached bacterial community, cloned sequences most closely related to the sulfate-reducing Desulfobacter and Desulfobulbus genera represented 20% of the bacterial community in wells where the concentration of sulfate in groundwater was high (> 0.2 mM), compared to only 3% in wells with less sulfate. Sequences related to the genus Geobacter, a genus containing ferric-iron reducers, were of nearly equal abundance (15%) to the sulfate reducers under high sulfate conditions, however their relative abundance increased to 34% when sulfate concentrations were < 0.03 mM. Also, in areas where sulfate concentrations were <0.03 mM, archaeal 16S rRNA gene sequences similar to those found in methanogens such as Methanosarcina and Methanosaeta comprised 73-80% of the community, and dissolved CH4 ranged between 220 and 1240 μM in these groundwaters. In contrast, methanogens (and their product, CH4) were nearly absent in samples collected from groundwater samples with > 0.2 mM sulfate. In the suspended fraction of wells where the concentration of sulfate was between 0.03 and 0.2 mM, the archaeal community was dominated by sequences most closely related to the ANME-2D, a group of archaea known for anaerobically oxidizing methane. Based on available energy (∆GA) estimations, results varied little for both sulfate reduction and methanogenesis throughout all wells studied, but could favor anaerobic oxidation of methane (AOM) in wells containing minimal sulfate and dihydrogen, suggesting AOM coupled with H2-oxidizing organisms such as sulfate or iron reducers could be an important pathway occurring in the Mahomet aquifer.

CONCLUSIONS

Overall, the results show several distinct factors control the composition of microbial communities in the Mahomet aquifer. Bacteria that respire insoluble substrates such as iron oxides, i.e. Geobacter, comprise a greater abundance of the attached community than the suspended regardless of groundwater chemistry. Differences in community structure driven by the concentration of sulfate point to a clear link between the availability of substrate and the abundance of certain functional groups, particularly iron reducers, sulfate reducers, methanogens, and methanotrophs. Integrating both geochemical and microbiological observations suggest that the relationships between these functional groups could be driven in part by mutualism, especially between ferric-iron and sulfate reducers.

Determining hot spots of fecal contamination in a tropical watershed by combining land-use information and meteorological data with source-specific assays

The objective of this study was to combine knowledge of environmental, topographical, meteorological, and anthropologic factors in the Río Grande de Arecibo (RGA) watershed in Puerto Rico with information provided by microbial source tracking (MST) to map hot spots (i.e., likely sources) of fecal contamination. Water samples were tested for the presence of human and bovine fecal contamination in addition to fecal indicator bacteria and correlated against several land uses and the density of septic tanks, sewers, and latrines. Specifically, human sources were positively correlated with developed (r = 0.68), barren land uses (r = 0.84), density of septic tanks (r = 0.78), slope (r = 0.63), and the proximity to wastewater treatment plants (WWTPs) (r = 0.82). Agricultural land, the number of upstream National Pollution Discharge Elimination System (NPDES) facilities, and density of latrines were positively associated with the bovine marker (r = 0.71; r = 0.74; and r = 0.68, respectively). Using this information, we provided a hot spot map, which shows areas that should be closely monitored for fecal contamination in the RGA watershed. The results indicated that additional bovine assays are needed in tropical regions. We concluded that meteorological, topographical, anthropogenic, and land cover data are needed to evaluate and verify the performance of MST assays and, therefore, to identify important sources of fecal contamination in environmental waters.

Detection of fecal bacteria and source tracking identifiers in environmental waters using rRNA-based RT-qPCR and rDNA-based qPCR assays

In this study, we evaluated the use of RT-qPCR assays targeting rRNA gene sequences for the detection of fecal bacteria in water samples. We challenged the RT-qPCR assays against RNA extracted from sewage effluent (n = 14), surface water (n = 30), and treated source water (n = 15) samples. Additionally, we applied the same assays using DNA as the qPCR template. The targeted fecal bacteria were present in most of the samples tested, although in several cases, the detection frequency increased when RNA was used as the template. For example, the majority of samples that tested positive for E. coli and Campylobacter spp. in surface waters, and for human-specific Bacteroidales, E. coli, and Enterococcus spp. in treated source waters were only detected when rRNA was used as the original template. The difference in detection frequency using rRNA or rDNA (rRNA gene) was sample- and assay-dependent, suggesting that the abundance of active and nonactive populations differed between samples. Statistical analyses for each population exhibiting multiple quantifiable results showed that the rRNA copy numbers were significantly higher than the rDNA counterparts (p < 0.05). Moreover, the detection frequency of rRNA-based assays were in better agreement with the culture-based results of E. coli, intestinal enterococci, and thermotolerant Campylobacter spp. in surface waters than that of rDNA-based assays, suggesting that rRNA signals were associated to active bacterial populations. Our data show that using rRNA-based approaches significantly increases detection sensitivity for common fecal bacteria in environmental waters. These findings have important implications for microbial water quality monitoring and public health risk assessments.

Metagenome analyses of corroded concrete wastewater pipe biofilms reveal a complex microbial system

Background

Concrete corrosion of wastewater collection systems is a significant cause of deterioration and premature collapse. Failure to adequately address the deteriorating infrastructure networks threatens our environment, public health, and safety. Analysis of whole-metagenome pyrosequencing data and 16S rRNA gene clone libraries was used to determine microbial composition and functional genes associated with biomass harvested from crown (top) and invert (bottom) sections of a corroded wastewater pipe.

Results

Taxonomic and functional analysis demonstrated that approximately 90% of the total diversity was associated with the phyla Actinobacteria, Bacteroidetes, Firmicutes and Proteobacteria. The top (TP) and bottom pipe (BP) communities were different in composition, with some of the differences attributed to the abundance of sulfide-oxidizing and sulfate-reducing bacteria. Additionally, human fecal bacteria were more abundant in the BP communities. Among the functional categories, proteins involved in sulfur and nitrogen metabolism showed the most significant differences between biofilms. There was also an enrichment of genes associated with heavy metal resistance, virulence (protein secretion systems) and stress response in the TP biofilm, while a higher number of genes related to motility and chemotaxis were identified in the BP biofilm. Both biofilms contain a high number of genes associated with resistance to antibiotics and toxic compounds subsystems.

Conclusions

The function potential of wastewater biofilms was highly diverse with level of COG diversity similar to that described for soil. On the basis of the metagenomic data, some factors that may contribute to niche differentiation were pH, aerobic conditions and availability of substrate, such as nitrogen and sulfur. The results from this study will help us better understand the genetic network and functional capability of microbial members of wastewater concrete biofilms.

Molecular survey of concrete sewer biofilm microbial communities

The microbial composition of concrete biofilms within wastewater collection systems was studied using molecular assays. SSU rDNA clone libraries were generated from 16 concrete surfaces of manholes, a combined sewer overflow, and sections of a corroded sewer pipe. Of the 2457 sequences analyzed, α-, β-, γ-, and δ-Proteobacteria represented 15%, 22%, 11%, and 4% of the clones, respectively. β-Proteobacteria (47%) sequences were more abundant in the pipe crown than any of the other concrete surfaces. While 178 to 493 Operational Taxonomic Units (OTUs) were associated with the different concrete samples, only four sequences were shared among the different clone libraries. Bacteria implicated in concrete corrosion were found in the clone libraries while archaea, fungi, and several bacterial groups were also detected using group-specific assays. The results showed that concrete sewer biofilms are more diverse than previously reported. A more comprehensive molecular database will be needed to better study the dynamics of concrete biofilms.

Comparative fecal metagenomics unveils unique functional capacity of the swine gut

Background

Uncovering the taxonomic composition and functional capacity within the swine gut microbial consortia is of great importance to animal physiology and health as well as to food and water safety due to the presence of human pathogens in pig feces. Nonetheless, limited information on the functional diversity of the swine gut microbiome is available.

Results

Analysis of 637, 722 pyrosequencing reads (130 megabases) generated from Yorkshire pig fecal DNA extracts was performed to help better understand the microbial diversity and largely unknown functional capacity of the swine gut microbiome. Swine fecal metagenomic sequences were annotated using both MG-RAST and JGI IMG/M-ER pipelines. Taxonomic analysis of metagenomic reads indicated that swine fecal microbiomes were dominated by Firmicutes and Bacteroidetes phyla. At a finer phylogenetic resolution, Prevotella spp. dominated the swine fecal metagenome, while some genes associated with Treponema and Anareovibrio species were found to be exclusively within the pig fecal metagenomic sequences analyzed. Functional analysis revealed that carbohydrate metabolism was the most abundant SEED subsystem, representing 13% of the swine metagenome. Genes associated with stress, virulence, cell wall and cell capsule were also abundant. Virulence factors associated with antibiotic resistance genes with highest sequence homology to genes in Bacteroidetes, Clostridia, and Methanosarcina were numerous within the gene families unique to the swine fecal metagenomes. Other abundant proteins unique to the distal swine gut shared high sequence homology to putative carbohydrate membrane transporters.

Conclusions

The results from this metagenomic survey demonstrated the presence of genes associated with resistance to antibiotics and carbohydrate metabolism suggesting that the swine gut microbiome may be shaped by husbandry practices.

Microbial community structure during nitrate and perchlorate reduction in ion-exchange brine using the hydrogen-based membrane biofilm reactor (MBfR)

Detoxification of perchlorate by microbial communities under denitrifying conditions has been recently reported, although the identity of the mixed populations involved in perchlorate reduction is not well understood. In order to address this, the bacterial diversity of membrane biofilm reactors (MBfR) set up under autotrophic denitrifying and perchlorate-reducing conditions were examined by analyses of the 16S rRNA gene sequences of clone libraries. Inocula from diverse locations were tested for their ability to reduce nitrate and perchlorate in synthetic ion exchange spent brine (45g/l NaCl) using H(2)-based MBfRs. Phylogenetic analysis of 16S rRNA gene sequences showed that proteobacterial species dominated the biofilm communities, particularly nitrate-reducing gamma-proteobacteria. Even though the inocula to the MBfRs came from different sources, clones closely related to Marinobacter hydrocarbonoclasticus represented 53% of all clones in the MBfR biofilms. The clone libraries contained no known perchlorate-reducing bacteria, which suggest that denitrifiers carried out perchlorate reduction, probably by secondary-utilization.

Identification of bacterial populations in drinking water using 16S rRNA-based sequence analyses

Intracellular RNA is rapidly degraded in stressed cells and is more unstable outside of the cell than DNA. As a result, RNA-based methods have been suggested to study the active microbial fraction in environmental matrices. The aim of this study was to identify bacterial populations in drinking water by analyzing 16S rRNA-based clone libraries. Hollow-fiber ultrafiltration was used to concentrate bacterial communities from 40l of tap water collected at 12 different times during three different summer months from a single point-of-use. Total RNA was extracted from the microbial concentrates and used to develop 16S rRNA-based clone libraries. Phylogenetic analyses of 1231 partial 16S rRNA gene sequences showed that difficult-to-classify bacterial sequences were the most predominant clones, representing 57.6% of the sequences analyzed. Within these unclassified clades, most sequences were closely related to sequences retrieved from previous DNA- and RNA-based drinking water studies. Other bacterial groups represented in this study included Proteobacteria, cyanobacteria, Actinobacteria, Bacteroidetes, and Planctomycetes. Overall, the results suggest that these bacterial groups are amongst potentially active bacteria in drinking water. Diversity analyses of clones generated show that while overall diversity is similar amongst the different months, membership changes with respect to time. The results from this study further improve our understanding of the molecular diversity and bacterial population dynamics of drinking water microbial communities. Moreover, these results provide the sequence foundation for the development of molecular assays that target active drinking water bacteria.

Quo vadis source tracking? Towards a strategic framework for environmental monitoring of fecal pollution

Advances in microbial source tracking (MST) have largely been driven by the need to comply with water quality standards based on traditional indicator bacteria. Recently, a number of culture-independent, and library-independent methods based on polymerase chain reaction (PCR) have been gaining popularity among source trackers. However, only a limited number of these methods have been successfully used in field applications, primarily due to the fact that many of them are still being developed. In this critical outlook, we examine different viewpoints associated with the practical use of MST to identify critical research gaps, propose a priority-based timeline to address them, and outline emerging technologies that will likely impact the future of source tracking. We propose that it is necessary to consider each of these aspects in order to advance towards a unifying framework in source identification, so that fecal pollution monitoring can be reliably used for comprehensive environmental microbial monitoring, to develop risk assessment models, and to implement and validate adequate management practices.

Identifying fecal sources in a selected catchment reach using multiple source-tracking tools

Given known limitations of current microbial source-tracking (MST) tools, emphasis on small, simple study areas may enhance interpretations of fecal contamination sources in streams. In this study, three MST tools-Escherichia coli repetitive element polymerase chain reaction (rep-PCR), coliphage typing, and Bacteroidales 16S rDNA host-associated markers-were evaluated in a selected reach of Plum Creek in south-central Nebraska. Water-quality samples were collected from six sites. One reach was selected for MST evaluation based on observed patterns of E. coli contamination. Despite high E. coli concentrations, coliphages were detected only once among water samples, precluding their use as a MST tool in this setting. Rep-PCR classification of E. coli isolates from both water and sediment samples supported the hypothesis that cattle and wildlife were dominant sources of fecal contamination, with minor contributions by horses and humans. Conversely, neither ruminant nor human sources were detected by Bacteroidales markers in most water samples. In bed sediment, ruminant- and human-associated Bacteroidales markers were detected throughout the interval from 0 to 0.3 m, with detections independent of E. coli concentrations in the sediment. Although results by E. coli-based and Bacteroidales-based MST methods led to similar interpretations, detection of Bacteroidales markers in sediment more commonly than in water indicates that different tools to track fecal contamination (in this case, tools based on Bacteroidales DNA and E. coli isolates) may have varying relevance to the more specific goal of tracking the sources of E. coli in watersheds. This is the first report of simultaneous, toolbox approach application of a library-based and marker-based MST analyses to flowing surface water.

Identification of bacterial DNA markers for the detection of human fecal pollution in water

We used genome fragment enrichment and bioinformatics to identify several microbial DNA sequences with high potential for use as markers in PCR assays for detection of human fecal contamination in water. Following competitive solution-phase hybridization of total DNA from human and pig fecal samples, 351 plasmid clones were sequenced and were determined to define 289 different genomic DNA regions. These putative human-specific fecal bacterial DNA sequences were then analyzed by dot blot hybridization, which confirmed that 98% were present in the source human fecal microbial community and absent from the original pig fecal DNA extract. Comparative sequence analyses of these sequences suggested that a large number (43.5%) were predicted to encode bacterial secreted or surface-associated proteins. Deoxyoligonucleotide primers capable of annealing to a subset of 26 of the candidate sequences predicted to encode factors involved in interactions with host cells were then used in the PCR and did not amplify markers in DNA from any additional pig fecal specimens. These 26 PCR assays exhibited a range of specificity in tests with 11 other animal sources, with more than half amplifying markers only in specimens from dogs or cats. Four assays were more specific, detecting markers only in specimens from humans, including those from 18 different human populations examined. We then demonstrated the potential utility of these assays by using them to detect human fecal contamination in several impacted watersheds.

Assessment of fecal pollution sources in a small northern-plains watershed using PCR and phylogenetic analyses of Bacteroidetes 16S rRNA gene

We evaluated the efficacy, sensitivity, host-specificity, and spatial/temporal dynamics of human- and ruminant-specific 16S rRNA gene Bacteroidetes markers used to assess the sources of fecal pollution in a fecally impacted watershed. Phylogenetic analyses of 1271 fecal and environmental 16S rRNA gene clones were also performed to study the diversity of Bacteroidetes in this watershed. The host-specific assays indicated that ruminant feces were present in 28-54% of the water samples and in all sampling seasons, with increasing frequency in downstream sites. The human-targeted assays indicated that only 3-5% of the water samples were positive for human fecal signals, although a higher percentage of human-associated signals (19-24%) were detected in sediment samples. Phylogenetic analysis indicated that 57% of all water clones clustered with yet-to-be-cultured Bacteroidetes species associated with sequences obtained from ruminant feces, further supporting the prevalence of ruminant contamination in this watershed. However, since several clusters contained sequences from multiple sources, future studies need to consider the potential cosmopolitan nature of these bacterial populations when assessing fecal pollution sources using Bacteroidetes markers. Moreover, additional data is needed in order to understand the distribution of Bacteroidetes host-specific markers and their relationship to water quality regulatory standards.

Competitive metagenomic DNA hybridization identifies host-specific microbial genetic markers in cow fecal samples

Several PCR methods have recently been developed to identify fecal contamination in surface waters. In all cases, researchers have relied on one gene or one microorganism for selection of host-specific markers. Here we describe the application of a genome fragment enrichment (GFE) method to identify host-specific genetic markers from fecal microbial community DNA. As a proof of concept, bovine fecal DNA was challenged against a porcine fecal DNA background to select for bovine-specific DNA sequences. Bioinformatic analyses of 380 bovine enriched metagenomic sequences indicated a preponderance of Bacteroidales-like regions predicted to encode membrane-associated and secreted proteins. Oligonucleotide primers capable of annealing to select Bacteroidales-like bovine GFE sequences exhibited extremely high specificity (>99%) in PCR assays with total fecal DNAs from 279 different animal sources. These primers also demonstrated a broad distribution of corresponding genetic markers (81% positive) among 148 different bovine sources. These data demonstrate that direct metagenomic DNA analysis by the competitive solution hybridization approach described is an efficient method for identifying potentially useful fecal genetic markers and for characterizing differences between environmental microbial communities.

Use of competitive DNA hybridization to identify differences in the genomes of bacteria

Although recent technological advances in DNA sequencing and computational biology now allow scientists to compare entire microbial genomes, comparisons of closely related bacterial species and individual isolates by whole-genome sequencing approaches remains prohibitively expensive for most laboratories. Here we report the development and testing of a biochemical approach for targeted sequencing of only those chromosomal regions that differ between two DNA preparations. The method, designated GFE (genome fragment enrichment) uses competitive solution hybridization and positive selection to obtain genomic DNA fragments that are present in one pool of fragments but not another. Repeated comparisons of the genomes of Enterococcus faecalis and E. faecium led to the identification of 225 putative genome-specific DNA fragments. Species and strain variations within these fragments were confirmed by both experimental and bioinformatic analyses. The E. faecalis genome-specific sequences identified included both a preponderance of those predicted to encode surface-exposed proteins, as well as several previously described unique marker regions embedded within highly conserved rrn operons. The GFE strategy we describe efficiently identified genomic differences between two enterococcal genomes, and will be widely applicable for studying genetic variation among closely related bacterial species.

Development and applications of microbial ecogenomic indicators for monitoring water quality: report of a workshop assessing the state of the science, research needs and future directions

This article brings forth recommendations from a workshop sponsored by the U.S. Environmental Protection Agency's Science to Achieve Results (STAR) and Environmental Monitoring and Assessment (EMAP) Programs and by the Council of State Governments, held during May 2002 in Kansas City, Kansas. The workshop assembled microbial ecologists and environmental scientists to determine what research and science is needed to bring existing molecular biological approaches and newer technologies arising from microbial genomic research into environmental monitoring and water quality assessments. Development of genomics and proteomics technologies for environmental science is a very new area having potential to improve environmental water quality assessments. The workshop participants noted that microbial ecologists are already using molecular biological methods well suited for monitoring and water quality assessments and anticipate that genomics-enabled technologies could be made available for monitoring within a decade. Recommendations arising from the workshop include needs for (i) identification of informative microbial gene sequences, (ii) improved understandings of linkages between indicator taxa, gene expression and environmental condition, (iii) technological advancements towards field application, and (iv) development of the appropriate databases.

Identification of active bacterial communities in a model drinking water biofilm system using 16S rRNA-based clone libraries

Recent phylogenetic studies have used DNA as the target molecule for the development of environmental 16S rRNA gene clone libraries. As DNA may persist in the environment, DNA-based libraries cannot be used to identify metabolically active bacteria in water systems. In this study, an annular reactor was used to generate model drinking water biofilms grown on polycarbonate slides. High-quality RNA was extracted from 2-month-old biofilms and used to generate 16S rRNA-based clones. Sequencing analyses of 16S rRNA-based clones suggested that the active bacterial fraction consisted of a few dominant bacterial groups related to Nevskia ramosa and to uncultured bacteria. Several of these bacterial groups were closely related to clones characterized in a DNA-based clone library also generated in this study. Altogether, these results suggest that some of the predominant drinking water bacteria identified using DNA-based techniques are indeed active.

Host distributions of uncultivated fecal Bacteroidales bacteria reveal genetic markers for fecal source identification

The purpose of this study was to examine host distribution patterns among fecal bacteria in the order Bacteroidales, with the goal of using endemic sequences as markers for fecal source identification in aquatic environments. We analyzed Bacteroidales 16S rRNA gene sequences from the feces of eight hosts: human, bovine, pig, horse, dog, cat, gull, and elk. Recovered sequences did not match database sequences, indicating high levels of uncultivated diversity. The analysis revealed both endemic and cosmopolitan distributions among the eight hosts. Ruminant, pig, and horse sequences tended to form host- or host group-specific clusters in a phylogenetic tree, while human, dog, cat, and gull sequences clustered together almost exclusively. Many of the human, dog, cat, and gull sequences fell within a large branch containing cultivated species from the genus Bacteroides. Most of the cultivated Bacteroides species had very close matches with multiple hosts and thus may not be useful targets for fecal source identification. A large branch containing cultivated members of the genus Prevotella included cloned sequences that were not closely related to cultivated Prevotella species. Most ruminant sequences formed clusters separate from the branches containing Bacteroides and Prevotella species. Host-specific sequences were identified for pigs and horses and were used to design PCR primers to identify pig and horse sources of fecal pollution in water. The primers successfully amplified fecal DNAs from their target hosts and did not amplify fecal DNAs from other species. Fecal bacteria endemic to the host species may result from evolution in different types of digestive systems.