Characterization of Coastal Urban Watershed Bacterial Communities Leads to Alternative Community-Based Indicators

doi:10.1371/journal.pone.0011285

PLoS One. 2010; 5(6): e11285.

Published online 2010 June 23. doi: 10.1371/journal.pone.0011285

PMCID: PMC2890573

Copyright Wu et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Characterization of Coastal Urban Watershed Bacterial Communities Leads to Alternative Community-Based Indicators

Cindy H. Wu,¹ Bram Sercu,² Laurie C. Van De Werfhorst,² Jakk Wong,¹ Todd Z. DeSantis,¹ Eoin L. Brodie,¹ Terry C. Hazen,¹ Patricia A. Holden,² and Gary L. Andersen¹^*

¹Ecology Department, Earth Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California, United States of America

²Donald Bren of School of Environmental Science and Management, University of California Santa Barbara, Santa Barbara, California, United States of America

Ching-Hong Yang, Editor

University of Wisconsin-Milwaukee, United States of America

* E-mail: GLAndersen/at/lbl.gov

Conceived and designed the experiments: CHW PAH GLA. Performed the experiments: CHW BS LCVDW. Analyzed the data: CHW BS JW TZD ELB TCH GLA. Contributed reagents/materials/analysis tools: PAH GLA. Wrote the paper: CHW BS ELB TCH PAH GLA.

Received February 10, 2010; Accepted May 1, 2010.

This article has been cited by other articles in PMC.

Abstract

Background

Microbial communities in aquatic environments are spatially and temporally dynamic due to environmental fluctuations and varied external input sources. A large percentage of the urban watersheds in the United States are affected by fecal pollution, including human pathogens, thus warranting comprehensive monitoring.

Methodology/Principal Findings

Using a high-density microarray (PhyloChip), we examined water column bacterial community DNA extracted from two connecting urban watersheds, elucidating variable and stable bacterial subpopulations over a 3-day period and community composition profiles that were distinct to fecal and non-fecal sources. Two approaches were used for indication of fecal influence. The first approach utilized similarity of 503 operational taxonomic units (OTUs) common to all fecal samples analyzed in this study with the watershed samples as an index of fecal pollution. A majority of the 503 OTUs were found in the phyla Firmicutes, Proteobacteria, Bacteroidetes, and Actinobacteria. The second approach incorporated relative richness of 4 bacterial classes (Bacilli, Bacteroidetes, Clostridia and α-proteobacteria) found to have the highest variance in fecal and non-fecal samples. The ratio of these 4 classes (BBC [ratio]

A) from the watershed samples demonstrated a trend where bacterial communities from gut and sewage sources had higher ratios than from sources not impacted by fecal material. This trend was also observed in the 124 bacterial communities from previously published and unpublished sequencing or PhyloChip- analyzed studies.

Conclusions/Significance

This study provided a detailed characterization of bacterial community variability during dry weather across a 3-day period in two urban watersheds. The comparative analysis of watershed community composition resulted in alternative community-based indicators that could be useful for assessing ecosystem health.

Introduction

Given that water sustains life, it is not surprising that a large percentage of the world's population lives near coastal regions [1], [2]. Coastal urban watersheds in the United States offer aesthetics and recreational value, serve as catchments for storm runoff, establish biological corridors for movements of wildlife, and provide buffers between developed areas and downstream waterways. As human populations increase, so does urbanization and lasting anthropogenic affects on creeks and coastal ecosystems [3]. According to a USEPA report (2007), 45% of streams and rivers, and 32% of bays and estuaries are impaired in the United States. Sources of impairment include pathogens and sewage discharges [4]. The presence of bacterial pollutants warrants comprehensive bacteriological characterization of these water bodies in order for us to understand their fate and transport in the environment.

Since pathogens often come from fecal sources, regulatory agencies require monitoring fecal indicator bacteria (FIB) for water quality assessments. Culture-dependent assays such as total coliform, fecal coliform and enterococci, and culture-independent assays such as quantitative PCR (qPCR) for Bacteroides and Bifidobacterium spp. [5] have been used as proxies for fecal pollution. However, enumeration of these indicator organisms often does not accurately represent the health of the ecosystem or associated risk [6] as these indicators are ubiquitous, persistent, regenerative [7], [8] and have low correlations with pathogen survival [9], [10] in the environment. Reliance upon single, even source-specific, markers of fecal pollution can be ineffective if they are labile or persistent relative to pathogens. The use of multiple indicators for tracking fecal contamination could circumvent the problem of single marker absence or presence and strengthen overall diagnoses of microbiological water quality [6], [7], [8], [9], [11].

With the advent of high throughput culture-independent characterization of microbial communities, such as microarray and sequencing approaches [12], [13], [14], [15], [16], detailed studies of bacterial community fluctuations due to physical, chemical and biological influences are now feasible. One such phylogenetic microarray, the PhyloChip, targets much of the known diversity within Bacteria and Archaea, and has been employed in a number of complex environments and conditions [17], [18], [19], [20], [21], [22], [23], [24], [25]. The current version (G2) of the PhyloChip provides the capability of identifying up to 8,741 Bacterial and Archaeal OTUs simultaneously [17], and allows for relative quantification of individual OTUs over a wide dynamic range [18], [26]. The highly parallel and reproducible nature of this array allows tracking community dynamics over time and treatment.

Bacterial communities in urban watersheds are sensitive to environmental perturbations and could provide information on impacts of fecal influence and overall ecosystem health. It is important to monitor the conditions of these watersheds because they are intricately tied in with downstream waterways, which could have public health risk and economic implications. Previous studies monitoring FIB most probable numbers (MPN) in urban creeks have found high temporal variability even during dry weather [27], [28], [29]. In Santa Barbara, California, exfiltration from sewer lines into the storm drain systems has been suspected to cause the observed high densities of FIB and human-specific Bacteroides markers (HBM) in urban watersheds that discharge into a recreational beach [29]. Here we analyze whole bacterial communities from the same Sercu et al. [29] samples in order to gain insights regarding the temporal and spatial dynamics of urban watershed bacterial community composition relevant to fecal pollution. Amplified 16S rRNA gene sequences from creek (including storm drains), lagoon and ocean sites in the Lower Mission Creek and Laguna watersheds in Santa Barbara, CA, along with 3 samples of fecal origin, were hybridized onto the PhyloChip for a complete microbial community analysis. Characterization of the whole bacterial community is crucial for understanding fluctuations of various bacterial groups, and could lead to more robust health risk indication by integrating data from multiple bacteria taxa. This work represents the first application of a comprehensive phylogenetic array for the purpose of characterizing urban watershed bacterial communities. Findings from this work suggest that such an approach could be useful for complementing multiple individual tests that are now typically employed to diagnose microbiological water quality related to public health.

Results

Resolving community differences by habitats

Samples were categorized into 4 habitat types: fecal, ocean, lagoon, and creek (Figure 1). Comparisons of Bray-Curtis distances of the communities, using Multi-Response Permutation Procedure (MRPP) [30], indicated significant differences between the samples from the different habitat types. Non-metric multidimensional scaling (NMDS) ordination illustrated that the bacterial communities were separated by habitat types for most of the samples, except for M2a and M2b (Figure 2). Salinity measurements at one of the lagoon sites (M2) were low, at ~1 ppt, on days 1 (M2a) and 2 (M2b) (Table S1). On day 3 (M2c), the salinity increased to 5.3 ppt, and a corresponding community composition shift was observed (Figure 2). The bacterial communities of M2a and M2b were more similar to creek samples with low salinity and M2c was more similar to the M4a and M4b lagoon samples, which had higher salinity measurements of 7.3–9.5 ppt. Lagoon sample M4c had lower salinity measurements and the community was more similar to creek samples than to M4a and M4b.

Figure 1

Sampling sites along Mission (M4–M9) and Laguna Channel (M2 and M3) watersheds.

Figure 2

NMDS plot of PhyloChip community distances.

Distributions of detected operational taxonomic units (OTUs) at the class level were compared among all habitat types, shown as relative richness (Figure 3A). The relative richness was normalized to the total number of OTUs detected in all of the samples from the same habitat type. We focused on classes that exhibited high variability of relative richness across the 4 habitats. The top 10 classes with the highest standard deviations were (in descending order): Clostridia, α-proteobacteria, Bacilli, γ-proteobacteria, β-proteobacteria, Actinobacteria, Flavobacteria, Bacteroidetes, Cyanobacteria and ε-proteobacteria. Of those classes, only Clostridia, Bacilli, and Bacteroidetes had higher relative richness in fecal samples than in creek, lagoon and ocean samples (Figure 3B). Only α-proteobacteria had lower richness in fecal samples than in creek, lagoon, and ocean samples (Figure 3C). The characteristics and potential of these 4 classes as indicators of fecal influence will be discussed further.

Figure 3

Bacterial community composition comparison across fecal, lagoon, creek and ocean samples.

Fecal sample-associated OTUs

In order to define bacteria that were common to all 3 fecal samples used in this study, a set of 503 OTUs, found in all fecal samples but not ubiquitous in the 27 watershed samples, were characterized and defined as fecal sample-associated OTUs (FSAO). The FSAO subpopulation consisted of 43% Firmicutes (out of the 503 OTUs), 28% Proteobacteria, 9% Bacteroidetes and 5% Actinobacteria (Figure S1). Of the Firmicutes (218 OTUs), 56% were from the order Clostridiales including the families Lachnospiraceae, Peptostreptococcaceae, Peptococcaceae, Acidaminococcaceae and Clostridiaceae; 17% were from the order Bacillales including Bacillaceae, Halobacillaceae, and Staphylococcaceae; and 17% were from Lactobacillales which included the families of Lactobacillaceae, Enterococcaceae and Streptococcaceae. In the Proteobacteria phylum (141 OTUs), 30% were from Enterobacteriales including Enterobacteriaceae; 7% were from Alteromonadales including Alteromonadaceae, and Shewanellaceae; 8% of the OTUs were from the order Burkholderiales including Burkholderiaceae, Comamonadaceae, Alcaligenaceae, Oxalobacteraceae, and Ralstoniaceae. The counts of FSAO for each of the three days are shown in Figure S2. The FSAO counts were highest at M9, M8, M6, M3 and M2 and lowest at M4 and M1. The 3-day average FSAO counts for sites M9, M6, M3, and M2 were significantly different (t-test, p-value<0.0001) from counts of M4, and M1.

Variable and stable subpopulations

PhyloChip analysis of subpopulations from each site for which the fluorescence intensities fluctuated the most (variable) and the least (stable) were examined over the course of the three-day sampling period. These variable and stable subpopulations consisted of OTUs from the top and bottom deciles after sorting based on variance of fluorescence intensity over the 3 days. A similarity metric, from the UniFrac [31] distance measure, was illustrated with boxplots for comparison of the median, upper and lower quartiles. Variable subpopulations of M6 were the most similar to the FSAO composition in comparison to the other sites (Figure 4A). Sites M9 and M3 were the second and third most similar to the FSAO. However, the similarity to FSAO for site M9 was not significantly different from that of M6 or M3. A pattern of decreasing similarity from M9, M6 and M3 to immediate downstream sites was illustrated. The majority of FSAO detected in the variable subpopulations was in the orders of Enterobacteriales (39 out of 58 FSAO detected in the variable subpopulation) for M6, Campylobacterales (6 out of 44) for M9, and Flavobacteriales (4 out of 31) for M3. The M9 stable subpopulation was the most similar to the FSAO, and was significantly different from the similarity to FSAO of all other sites (Figure 4B). Many of the FSAO in the M9 stable subpopulation were in the order of Bacillales (17 out of 47).

Figure 4

Boxplots of UniFrac similarity metrics between water and fecal-sample-associated OTUs (FSAO).

Ratio of Bacilli, Bacteroidetes and Clostridia to α-proteobacteria

Four bacterial classes, which exhibited highly fluctuating relative richness across the habitat types, were further explored as representatives of the fecal bacterial community (Figure 3A). The combined percentage of Bacilli, Bacteroidetes and Clostridia relative richness was 28.5% of total detected in the fecal samples, whereas in creek, lagoon and ocean they were less than 13.5% (Figure 3B). Almost 15% of the relative richness in creek water, lagoon and ocean samples were α-proteobacteria, while the percentage of α-proteobacteria found in fecal samples was 7% (Figure 3C). The relative richness ratio of Bacilli, Bacteroidetes and Clostridia to α-proteobacteria (BBC [ratio]

A) for fecal samples was more than 4-fold higher than the ratios of the other habitat types (Figure 3D). The BBC [ratio]

A ratio was calculated for each of the samples from the different sites (Figure 5). Site M6 exhibited the highest BBC [ratio]

A, and sites M1 and M4 had low BBC [ratio]

A ratios compared to the rest of the sites.

Figure 5

Bacillus, Bacteroidetes, Clostridia to α-proteobacteria ratios (BBC [ratio]

A) from each site.

Retrospective comparison of BBC [ratio]

A ratios from 16S rRNA gene clone library sequencing- and PhyloChip-analyzed samples

The BBC

A ratios of 124 communities characterized by clone-library sequencing and PhyloChip were compared (Figure 6). Detailed descriptions of the communities are included in Table S2. From published sequencing studies, we calculated the BBC [ratio]

A ratios of bacterial communities from 54 mammalian intestines [32], 5 sewage-associated samples [33], [34], [35], [36], [37], and 19 non-fecal samples [23], [38], [39], [40], [41], [42], [43], [44], [45], [46], [47], [48], [49], [50], [51], [52], [53]. Likewise, from PhyloChip-analyzed samples, we determined the BBC [ratio]

A ratios from communities of 11 gut [Brodie et al., unpublished; Marchesi et al., unpublished; Nguyen et al., unpublished; This study] [54], 17 sewage-associated [Conrad et al., unpublished; Sercu et al. unpublished; Wu et al., unpublished], and 18 non-fecal samples [Sercu et al., unpublished; This study] [55]. Anoxic non-fecal samples were included in this comparison as well. For both PhyloChip- and library sequencing-analyzed bacterial communities, gut and sewage-associated samples generally had higher BBC [ratio]

A ratios than non-fecal samples, except for anoxic non-fecal samples, which had an overlapping range with sewage-associated samples. There were also a few communities that did not follow the general BBC [ratio]

A ratio trend. The community of a nitrifying-denitrifying activated sludge [35] had much lower BBC [ratio]

A ratio than the rest of the sequenced sewage-associated communities. Also, beetle posterior hindgut and midgut communities had lower BBC [ratio]

A ratios than beetle anterior hindgut communities and the other PhyloChip-analyzed gut samples [Nguyen et al., unpublished].

Figure 6

Bacillus, Bacteroidetes, Clostridia to α-proteobacteria ratios (BBC [ratio]

A) of communities analyzed by sequencing or PhyloChip.

Discussion

Microbial communities in surface waters are highly responsive to perturbation, shifting with tidal cycles [56], salinity gradients [57], [58], dissolved organic matter concentration [59], and chemical stress [60], [61], [62]. The detection of short-term fluctuations in community composition suggests changes in environmental conditions, nutrients or bacterial sources. An effect of increased salinity due to tidal influence on bacterial composition was observed in this study where the coastal lagoon communities were more similar to creek communities with comparable salinity measurements (Figure 2). Salinity was more strongly correlated to community composition than the other environmental variables measured based on canonical correspondence analysis (data not shown). This result corroborated observations by others [63], [64], [65]. In addition to being highly sensitive to environmental fluctuations, the response time of community composition shift was within a 24-hour period.

The detection of this rapid community response could be useful for indication of external bacterial inputs, such as from fecal sources. FSAO, derived from the human fecal and untreated sewage samples, were used to represent fecal communities. One caveat is that the OTUs in the FSAO list are specific to the 3 fecal samples used in this study, and do not represent all fecal communities in all environments. However, the prevalent bacterial phyla found in the FSAO are the same as those observed in published studies of human gastrointestinal tract samples [66], [67], [68], [69] and turkey cecal samples [70]. Therefore, community similarity to FSAO could potentially indicate the presence of fecal bacteria. This hypothesis was tested by comparing the community distances between FSAO and variable/stable subpopulations at each of the site (Figure 4A and 4B).

Examination of the variable and stable subpopulations brings to light the bacterial temporal fluctuations across the 3 days. The variable subpopulation represents OTUs with highly fluctuating relative abundances, perhaps due to rapid growth, decay or large sporadic influx of bacterial sources. The stable subpopulation represents OTUs with constant relative abundances. These stable subpopulation OTUs are likely associated with endemic bacteria that are able to grow and persist under the in situ environmental conditions or are from consistent external sources.

UniFrac analysis showed that the variable subpopulation of M6 was the most similar to the FSAO (Figure 4A). This suggested intermittent exposure to fecal sources at this site, which was supported by elevated but numerically variable HBM densities and FIB MPN (Figure S3). The prevalence of Enterobacteriales in the variable subpopulation falls in line with the high FIB MPN observed at site M6, and further supports the use of similarity of the variable subpopulation with FSAO for demonstrating fecal pollution. Similarity of M9 variable subpopulation to FSAO was not significantly different from that of the M6 (Figure 4A). This indicated that there were OTUs in the M9 variable subpopulation that were also found in the FSAO, but they were mostly from the order of Campylobacterales, and not represented by FIB or HBM detection. The similarity to FSAO decreased gradually from drains to downstream sites (i.e. M9 to M7 and M6 to M4), illustrating possible fecal community presence at the drains and die-off or dilution effects as the communities flow downstream.

Interestingly, the stable subpopulation at M9 was most similar to FSAO out of all the sites, even though the FIB densities met the California water quality standards on 2 out of the 3 days and no HBM was detected (Figure 4B and Figure S3). The non-detection of HBM at M9 could be due to Bacteroides DNA concentration being below the quantitative PCR detection limit of 0.5×10³–10⁴ targets L⁻¹ [29] or that the fecal source was non-human. The top three families present in the M9 stable subpopulation were Bacillaceae, Staphylococcaceae and Lachnospiraceae. While Bacillaceae and Staphylococcaceae have been observed in non-aquatic environments [22], [26], Lachnospiraceae are primarily associated with cow rumen [71], human bowel [67] and anaerobic digesters [72]. Therefore, the data suggested that some of the OTUs detected at M9 could have a fecal, but non-human, origin. However, further confirmatory work is needed to distinguish between a consistent fecal source or bacterial re-growth as the cause for the similarity between M9 stable subpopulation and FSAO.

The FSAO includes OTUs that contain fecal coliforms, which have been demonstrated to re-grow and persist in the environment leading to false-positive water quality diagnoses [6], [8], [73]. This study further explores the potential of using alternative organisms that are independent of coliforms as fecal indicators by introducing the BBC [ratio]

A ratio. The ratio excludes coliform bacteria, thus, potentially avoids false-positive results associated with coliforms, and integrates counts for organisms widespread in non-fecal “pristine” environments to assess ecosystem health.

Bacteroidetes and Clostridia are enriched within the gut microbiota of many mammals [32], [66], [67], [68], [69], [70], and specific species within these 2 classes have been proposed as fecal indicators [5], [10], [74]. However, they are also found in anoxic saline aquatic environments [40], [45], [49], estuaries [38], the deep ocean [41], and high elevation lakes [59]. The class of Bacilli, which includes the indicator species Enterococcus, is commonly found in fecal samples such as the human gastrointestinal tract [69], turkey intestines [69], [70] and aerobic thermophilic swine wastewater bioreactors [75]. All 3 classes are dominant groups found in a chicken fecal metagenomic study [76]. α-proteobacteria, have been found as primary surface colonizers in coastal marine waters [77] and have the ability to thrive under low-nutrient conditions [56]. The BBC [ratio]

A ratio incorporates the relative richness of OTUs prevalent in these 4 bacterial classes associated with fecal and non-fecal samples to reflect possible fecal inputs, rather than the use of single organism presence or absence. Previous studies have suggested the use of ratios for indicating human or non-human fecal pollution [78], determining fecal age and enteric viral content [79], [80], representing the nutrient status of soil ecosystems [81], [82], identifying land use in wetland soils [83], and eutrophy in aquatic systems [84].

In order to assess the applicability of the observations from our watersheds to other samples, we calculated the BBC [ratio]

A ratio from previously published and unpublished studies (Table S2). BBC [ratio]

A ratios of gut samples analyzed by DNA sequencing or PhyloChip are not completely comparable, mainly due to differences in sample processing including primers used, PCR conditions and coverage differences across phylogenetic groups on the PhyloChip. However, within communities analyzed by sequencing from different research groups employing varying protocols, the gut, sewage-associated and non-fecal samples exhibited the same BBC [ratio]

A ratio trend as those communities analyzed by PhyloChip processed with a consistent standardized protocol. The distribution of BBC [ratio]

A ratios from these studies illustrates that gut and sewage-associated samples have higher BBC [ratio]

A ratio than non-fecal samples regardless of analysis methods (Figure 6). Anoxic non-fecal polluted environments also have similar ratios of BBC [ratio]

A as sewage-associated samples (Figure 6). This is most likely an attribute of similar growth conditions favoring both anaerobic and fecal bacteria. The indication of anoxic non-fecal environments is often times pertinent for determining public health risks. Anoxic conditions could lead to eutrophication in both fresh and salt water environments, which changes nutrient cycling, water quality and biodiversity [84]. Eutrophication has led to toxic algal blooms that adversely affect human and wildlife health [85], [86].

Kendall rank correlation of FIB, HBM, FSAO and BBC [ratio]

A ratios from all sites indicated significant positive correlations of BBC [ratio]

A ratios with HBM, total coliform, enterococcus and FSAO counts, but not with E. coli (Table S3). However, many of the samples had reached the total coliform measurement maximum detection limit of 24,196 MPN, therefore, the correlation of total coliform with BBC [ratio]

A ratio might be misleading. The result also illustrated that even though the BBC [ratio]

A ratio did not contain fecal coliforms, the fecal pollution pattern was similar to that indicated by the FSAO where coliforms were included. The drain site M6 was the only site where all lines of evidence, i.e. similarity of variable subpopulations to FSAO, FIB, HBM, and BBC [ratio]

A ratios, pointed to the presence of fecal contamination. At site M1 (ocean), all data indicated a community with the least fecal influence. The data for the rest of the sites (M2, M3, M4, M5, M7, M8 and M9) indicated varying degrees of influence by fecal sources. Also, communities from drains (M6 and M9) were the most similar to organisms found in the fecal samples, although different fecal organisms were detected in the two drains.

Knowledge of who is there and how they change over time and location is the hallmark of an ecosystems approach to studying urban watersheds. We used this concept to track the microbial community dynamics over a three-day period at a location with a history of frequent fecal contamination. In spite of the confounding effect of the movement of water through this watershed, several patterns that correlated with the presence of human fecal contamination were observed. By using the PhyloChip we are able to identify a significantly greater number of bacterial OTUs than is typically examined in coastal watersheds. Comparison of the microbial inventory of the watershed samples with local sewage samples and a human fecal sample led to the identification of specific organisms that were associated with either potential human fecal sources or with the watershed. From this information we observed 503 OTUs that were common to the three fecal samples (FSAO) and the ratios of observed classes of organisms that demonstrated the largest differences between human fecal sources and the receiving waters (BBC [ratio]

A ratio). Whereas most research for measuring fecal influences on coastal watersheds uses a bottom-up approach to hypothesize that a specific organism is representative of the source, we employed a top-down approach that looked at a large number of potential bacterial contaminants from a majority of the known bacterial diversity to identify a diverse collection of organisms associated with fecal pollution. The advantage of this approach is that we can use the findings of the BBC [ratio]

A ratio and the FSAO as the basis for additional bottom-up, controlled experiments to examine their applicability at other locations and with other human fecal sources. Using this more detailed microbial community characterization, it may be possible to move away from generic, single indicators to a community-indicator approach for assessing fecal contamination or environments conducive to pathogen growth.

Materials and Methods

Ethics statement

The Human Subjects Committee of University of California, Santa Barbara was informed of the anonymous human sample used in this study, and declared that the sample did not meet the definition of a human subject sample, therefore, no approval was necessary for it's use.

Sample description, collection and extraction

Mission Creek and Laguna Channel flow through an urbanized area of downtown Santa Barbara and discharge at a popular bathing beach. As described previously [29], water column samples from 3 consecutive days (a

day 1, b

day 2, c

day 3), during the dry season (June 2005), were collected from 9 locations (M1–M9) within the Mission Creek and Laguna watersheds in Santa Barbara, California (Figure 1). Samples were delineated into different habitat types: creek (M3, M5–M9, where M6 and M9 were from drains), lagoon (M2 and M4), and ocean (M1). One sample per day was collected at approximately the same time on each of the 3 days. No rain occurred at least 48 hours prior to or during the sampling. The creek flow rate, taken at M5, was 0.016 m³s⁻¹. Both watersheds discharged into the same lagoon at M2 and M4. Surface water flowed from the lagoon into the ocean (M1) at the time of sampling. Three fecal samples, 1 human feces (H), from Santa Barbara, and 2 raw sewage, from the influent at El Estero Wastewater Treatment plant (Santa Barbara, CA) (S1, S2), were also collected. Dissolved oxygen (DO), pH, temperature and salinity were measured along with each sampling [29]. Water samples were filtered in the lab onto 0.22 µm filters on the day of the sampling and stored at −20°C until nucleic acid extractions. DNA was extracted using the UltraClean Water DNA kit (MoBio Laboratories, Inc. Carlsbad, CA, USA), and archived at −20°C. Concentrations of fecal indicator bacteria (FIB) which includes total coliforms, E. coli, and Enterococcus spp., and quantitative PCR (qPCR) measurements of Human-specific Bacteroides Marker (HBM) were reported previously [29].

16S rRNA gene amplification for microarray analysis

Genes encoding 16S rRNA were amplified from the gDNA using non-degenerate Bacterial primers 27F and 1492R [87]. Polymerase chain reaction (PCR) was carried out using the TaKaRa Ex Taq system (Takara Bio Inc, Otsu, Japan). The amplification protocol was previously described [17].

Microarray processing, and image data analysis

Microarray analysis was performed using the PhyloChip, an Affymetrix-platform microarray. The protocols were previously reported [17]. Briefly, amplicons were concentrated to a volume less than 40 µl by isopropanol precipitation. The DNA amplicons were then fragmented with DNAse (Invitrogen, Carlsbad, CA, USA), biotin labeled, denatured, and hybridized to the DNA microarray at 48°C overnight (>16 hr). The arrays were subsequently washed and stained. Reagents, conditions, and equipments involved are detailed elsewhere [88]. Arrays were scanned using a GeneArray Scanner (Affymetrix, Santa Clara, CA, USA).

The CEL files obtained from the Affymetrix software that produced information about the fluorescence intensity of each probe were analyzed. The detailed criteria for scoring the probe fluorescence intensities were described elsewhere [17], [18], [89]. Briefly, a probe set consisted of 11 or more specific 25-mers (probes) that were prevalent in members of a given OTU but were dissimilar from sequences outside the given OTU. Probes with sequences complementing all 25 base pairs of the target sequences were termed perfect match (PM) probes. Each PM probe was matched with a control 25-mer, identical in all positions except the 13^th base, termed mismatch (MM) probe. The PM and MM constituted a probe pair that were analyzed together. The probe pairs were scored as positive if the following two criteria were met: 1) the intensity of fluorescence from the PM probe was greater than 1.3 times the intensity from the MM probe, and 2) the difference in intensity (PM minus MM), was at least 500 times greater than the squared noise value. The CEL files from this study are available upon request.

The taxonomic position of each OTU as well as the accompanying NCBI accession numbers of the sequences composing each OTU can be viewed in outline format at: http://greengenes.lbl.gov/Download/Taxonomic_Outlines/G2_chip_SeqDescByOTU_tax_outline.txt.

PhyloChip data normalization

PhyloChip data normalization was performed using R [90]. To correct for variation associated with quantification of amplicon target (quantification variation), and downstream variation associated with target fragmentation, labeling, hybridization, washing, staining and scanning (microarray technical variation) a two-step normalization procedure was developed. First, for each PhyloChip experiment, a scaling factor best explaining the intensities of the spiked control probes under a multiplicative error model was estimated using a maximum-likelihood procedure [54]. The intensities in each experiment were multiplied with its corresponding optimal scaling factor. Second, the intensities for each experiment were corrected for the variation in total array intensity by dividing the intensities with its corresponding total array intensity separately for Bacteria and Archea. The normalized data is available in Table S4.

Statistical Analysis

All statistical analyses were carried out in R [90], except for the canonical correspondence analysis (CCA). Bray-Curtis distances were calculated using normalized fluorescence intensity with the ecodist package [91]. Non-metric multidimensional scaling (NMDS) and multi response permutation procedure (MRPP) was performed using the vegan package. Student t-test and Kendall rank correlation from the stats package were used to compare samples. A relaxed neighbor-joining tree was generated using Clearcut [92] and used for UniFrac analysis [31]. Unweighted UniFrac distances, converted to similarity metrics, were calculated for FSAO, variable and stable subpopulations. CCA was carried out using PCOrd [93]. There were no DO, pH and salinity data for sampling days 1 and 2 for site 6, and all 3 days of sampling for site 8. No environmental variables were measured for fecal sample data. Therefore, best-estimate values were inserted based on values measured from the nearest sites on the same day for the CCA. Fecal sample environmental variables were estimated based on reported values in literature.

PhyloChip derived parameters

Unless otherwise stated, an OTU was considered present when at least 90% of its assigned probe pairs for its corresponding probe set were positive (positive fraction ≥0.9). For example, if 10 out of 11 probe pairs are positive, the positive fraction is 0.909 and the OTU is considered present.

Fecal-sample associated OTUs (FSAO) - OTUs that were present in all 3 fecal samples, and in all 27 water samples were tabulated separately. The list of 503 FSAO was derived by removing those OTUs found in all 27 water samples from the OTUs that were present in the fecal samples. The OTUs in each sample which were also found on the list of 503 FSAO were tallied and presented as the FSAO count.

Variable and stable subpopulations - OTUs that were present in at least one of the 3 samples from each site were tabulated and variances of the fluorescence intensities across the 3 days for those OTUs were generated. The OTUs were sorted by variance in descending order. The OTUs in the top deciles (90^th percentile) were defined as the variable subpopulation, and OTUs in the bottom deciles (10^th percentile) were defined as the stable subpopulation.

The BBC

A ratio of phyloChip samples - The number of OTUs in the classes of Bacilli (Bac), Bacteroidetes (Bct), Clostridia (Cls), and α-proteobacteria (A) where the positive fraction equal to 1 were tallied. The ratio was calculated using the following formula:

A mathematical equation, expression, or formula.
Object name is pone.0011285.e001.jpg

where

A mathematical equation, expression, or formula.
Object name is pone.0011285.e002.jpg

A mathematical equation, expression, or formula.
Object name is pone.0011285.e003.jpg

A mathematical equation, expression, or formula.
Object name is pone.0011285.e004.jpg

A mathematical equation, expression, or formula.
Object name is pone.0011285.e005.jpg

The count for unique OTUs in each of the class was normalized by dividing by the total number of OTUs in each class detectable by the G2 PhyloChip. The denominators were predetermined based on the number of OTUs assigned for each bacterial class on the G2 PhyloChip design.

The BBC

A ratio of published 16S rRNA gene clone library sequencing samples

Aligned sequences in the Greengenes [94] database were downloaded and re-classified using the PhyloChip (G2) taxonomy on the Greengenes website (http://greengene.lbl.gov). Aligned DNA sequences of various environmental communities were also obtained from [63]. The counts of unique OTUs were tallied for each bacterial class. The BBC [ratio]

A ratios were calculated using the formulas mentioned above. If no OTU was detected for that class, the count was set to 0.5.

Supporting Information

Figure S1

Phylum level profile of 503 fecal sample-associated OTUs (FSAO), and order level profiles of Firmicutes and Proteobacteria. Pie chart illustrates that the FSAO consist of 43% Firmicutes and 28% Proteobacteria. Most of the Firmicutes OTUs are in the order of Clostridiales, and most of the Proteobacteria OTUs are in Enterobacteriales.

(0.31 MB TIF)

Click here for additional data file.^{(302K, tif)}

Figure S2

Counts of fecal-sample-associated OTUs (FSAO) at each site. Each bar represents one sample from each day. OTUs in each sample which were also found on the list of the 503 FSAO were tallied and presented as the FSAO count.

(0.15 MB TIF)

Click here for additional data file.^{(144K, tif)}

Figure S3

Measurements of Human-specific Bacteroides Marker (HBM), Total Coliform (TC), E. coli (EC), and Enterococcus (ENT) counts. Bars represent HBM values. Lines represent TC, EC and ENT most probable number (MPN).

(0.49 MB TIF)

Click here for additional data file.^{(475K, tif)}

Table S1

Environmental variables measured concurrently with the bacterial community samples. Dissolved oxygen, temperature, salinity and pH were measured at the time of sampling and reported here.

(0.42 MB TIF)

Click here for additional data file.^{(410K, tif)}

Table S2

Description of bacterial communities analyzed by sequencing and PhyloChip used in Figure 6. Gut, sewage-associated and non-fecal samples analyzed by clone-library sequencing and PhyloChip used for the Bacilli, Bacteroidetes, Clostridia to α-proteobacteria ratio (BBC [ratio]

A ratio) are described. All DNA sequences from sequencing samples had a minimum length of 1250 base pairs, except for those with the (*) symbol where the minimum sequence length was 200 base pairs.

(0.90 MB TIF)

Click here for additional data file.^{(874K, tif)}

Table S3

Kendall rank correlation tau coefficient and p-values (in parenthesis). Measurements from all 27 water samples were used. The (*) symbol denotes statistical significance (p-value<0.05) differences. Abbreviations: Human Bacteroides Marker (HBM); total coliform (TC); E. coli (EC); enterococcus (ENT); fecal-sample associated OTUs (FSAO); Bacilli, Bacteroidetes, and Clostridia to α-proteobacteria ratio (BBC [ratio]

A).

(0.09 MB TIF)

Click here for additional data file.^{(83K, tif)}

Table S4

Total OTUs detected by PhyloChip for all 30 samples. Positive fraction and normalized fluorescence intensity values are reported.

(3.87 MB XLS)

Click here for additional data file.^{(3.6M, xls)}

Acknowledgments

We are grateful to Yvette M. Piceno (LBNL) for PhyloChip processing help, Ulas Karaoz (LBNL) for advice on R scripts, and Sherry Seybold (LBNL) for making Figure 1. We would also like to thank Catherine Lozupone for providing aligned DNA sequence files, as well as Julian Marchesi, Nhu Nguyen and Mark Conrad for providing PhyloChip data files. Additional assistance from the City of Santa Barbara and from UCSB researchers was as per Sercu et al., 2009.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: Part of this work was performed at Lawrence Berkeley National Laboratory under Department of Energy contract number DE-AC02-05CH11231 and funded by the Rathmann Family Foundation and the California State Water Resources Control Board Prop. 50 Clean Beaches Initiative grant with additional assistance by the Department of Public Health, Marin County California (CHW, JW, TZD, ELB, TCH, GLA). Additional funding was provided by the City of Santa Barbara through Measure B funding, the California State Water Resources Control Board Prop. 50 Clean Beaches Initiative grant and by the Switzer Foundation through a Leadership Grant. Flow data were provided through the NSF-funded Santa Barbara Long Term Ecological Research project, NSF OCE 9982105 and OCE 0620276 (BS, LCVDW, PAH). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Stewart JR, Gast RJ, Fujioka RS, Solo-Gabriele HM, Meschke JS, et al. The coastal environment and human health: microbial indicators, pathogens, sentinels and reservoirs. Environ Health. 2008;7:1476–1069X. [PMC free article] [PubMed]

2. UNEP. Marine and coastal ecosystems and human well-being: A synthesis report based on the findings of the Millennium Ecosystem Assessment. Nairobi, Kenya: United Nations Environmental Program; 2006. pp. 1–76.

3. Platt RH. Urban watershed management: Sustainability, one stream at a time. Environment. 2006;48:38–38.

4. USEPA. National Water Quality Inventory: Report to Congress, 2002 Reporting Cycle. 2007. US Environmental Protection Agency.

5. Savichtcheva O, Okabe S. Alternative indicators of fecal pollution: Relations with pathogens and conventional indicators, current methodologies for direct pathogen monitoring and future application perspectives. Water Res. 2006;40:2463–2476. [PubMed]

6. Hazen TC, Toranzos GA. In: Tropical source water. McFeters GA, editor. New York: Springer-Verlag; 1990.

7. Field KG, Bernhard AE, Brodeur TJ. Molecular approaches to microbiological monitoring: Fecal source detection. Environl Monit Assess. 2003;81:313–326. [PubMed]

8. Harwood VJ, Levine AD, Scott TM, Chivukula V, Lukasik J, et al. Validity of the indicator organism paradigm for pathogen reduction in reclaimed water and public health protection. Appl Environ Microbiol. 2005;71:3163–3170. [PMC free article] [PubMed]

9. Savichtcheva O, Okayama N, Okabe S. Relationships between Bacteroides 16S rRNA genetic markers and presence of bacterial enteric pathogens and conventional fecal indicators. Water Res. 2007;41:3615–3628. [PubMed]

10. Wery N, Lhoutellier C, Ducray F, Delgenes JP, Godon JJ. Behaviour of pathogenic and indicator bacteria during urban wastewater treatment and sludge composting, as revealed by quantitative PCR. Water Res. 2008;42:53–62. [PubMed]

11. Shibata T, Solo-Gabriele HM, Fleming LE, Elmir S. Monitoring marine recreational water quality using multiple microbial indicators in an urban tropical environment. Water Res. 2004;38:3119–3131. [PMC free article] [PubMed]

12. Andersen GL, He Z, DeSantis TZ, Brodie EL, Zhou J. The Use of Microarrays in Microbial Ecology. Environ Mol Microbiol. 2010:87–109.

13. Cardenas E, Tiedje JM. New tools for discovering and characterizing microbial diversity. Curr Opin Biotechnol. 2008;19:544–549. [PubMed]

14. Hamady M, Knight R. Microbial community profiling for human microbiome projects: Tools, techniques, and challenges. Genome Res. 2009;19:1141–1152. [PubMed]

15. Nelson KE, Bryan PA, White BA. Genomics and Metagenomics: History and Progress. Environ Mol Microbiol. 2010:21–36.

16. Raes J, Bork P. Systems microbiology - Timeline - Molecular eco-systems biology: towards an understanding of community function. Nat Rev Microbiol. 2008;6:693–699. [PubMed]

17. Brodie EL, DeSantis TZ, Joyner DC, Baek SM, Larsen JT, et al. Application of a high-density oligonucleotide microarray approach to study bacterial population dynamics during uranium reduction and reoxidation. Appl Environ Microbiol. 2006;72:6288–6298. [PMC free article] [PubMed]

18. Brodie EL, DeSantis TZ, Parker JPM, Zubietta IX, Piceno YM, et al. Urban aerosols harbor diverse and dynamic bacterial populations. Proc Nat Acad Sci. 2007;104:299–304. [PubMed]

19. Cruz-Martinez K, Suttle KB, Brodie EL, Power ME, Andersen GL, et al. Despite strong seasonal responses, soil microbial consortia are more resilient to long-term changes in rainfall than overlying grassland. ISME J. 2009;3:738–744. [PubMed]

20. DeAngelis KM, Brodie EL, DeSantis TZ, Andersen GL, Lindow SE, et al. Selective progressive response of soil microbial community to wild oat roots. ISME J. 2009;3:168–178. [PubMed]

21. Flanagan JL, Brodie EL, Weng L, Lynch SV, Garcia O, et al. Loss of bacterial diversity during antibiotic treatment of intubated patients colonized with Pseudomonas aeruginosa. J Clin Microbiol. 2007;45:1954–1962. [PMC free article] [PubMed]

22. Sagaram US, DeAngelis KM, Trivedi P, Andersen GL, Lu SE, et al. Bacterial Diversity Analysis of Huanglongbing Pathogen-Infected Citrus, Using PhyloChip Arrays and 16S rRNA Gene Clone Library Sequencing. Appl Environ Microbiol. 2009;75:1566–1574. [PMC free article] [PubMed]

23. West NJ, Obernosterer I, Zemb O, Lebaron P. Major differences of bacterial diversity and activity inside and outside of a natural iron-fertilized phytoplankton bloom in the Southern Ocean. Environ Microbiol. 2008;10:738–756. [PubMed]

24. Wrighton KC, Agbo P, Warnecke F, Weber KA, Brodie EL, et al. A novel ecological role of the Firmicutes identified in thermophilic microbial fuel cells. ISME J. 2008;2:1146–1156. [PubMed]

25. Cox MJ, Huang YJ, Fujimura KE, Liu JT, McKean M, et al. Lactobacillus casei Abundance Is Associated with Profound Shifts in the Infant Gut Microbiome. PLoS One. 2010;5:e8745. [PMC free article] [PubMed]

26. La Duc MT, Dekas A, Osman S, Moissl C, Newcombe D, et al. Isolation and characterization of bacteria capable of tolerating the extreme conditions of clean room environments. Appl Environ Microbiol. 2007;73:2600–2611. [PMC free article] [PubMed]

27. Isovitsch Parks SL, VanBriesen JM. Evaluating Temporal Variability in Bacterial Indicator Samples for an Urban Watershed. Journal of Environmental Engineering. 2009;32:1294–1303.

28. Petersen TM, Rifai HS, Suarez MP, Stein AR. Bacteria loads from point and nonpoint sources in an urban watershed. J Environ Eng-ASCE. 2005;131:1414–1425.

29. Sercu B, Van De Werfhorst LC, Murray J, Holden PA. Storm Drains are Sources of Human Fecal Pollution during Dry Weather in Three Urban Southern California Watersheds. Environ Sci Technol. 2009;43:293–298. [PubMed]

30. Mielke PW, Berry KJ. Permutation tests for common locations among samples with unequal variances. J Educ Behav Stat. 1994;19:217–236.

31. Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71:8228–8235. [PMC free article] [PubMed]

32. Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, et al. Evolution of mammals and their gut microbes. Science. 2008;320:1647–1651. [PMC free article] [PubMed]

33. Chouari R, Le Paslier D, Daegelen P, Ginestet P, Weissenbach J, et al. Novel predominant archaeal and bacterial groups revealed by molecular analysis of an anaerobic sludge digester. Environ Microbiol. 2005;7:1104–1115. [PubMed]

34. Chouari R, Le Paslier D, Dauga C, Daegelen P, Weissenbach J, et al. Novel major bacterial candidate division within a municipal anaerobic sludge digester. Appl Environ Microbiol. 2005;71:2145–2153. [PMC free article] [PubMed]

35. Juretschko S, Loy A, Lehner A, Wagner M. The microbial community composition of a nitrifying-denitrifying activated sludge from an industrial sewage treatment plant analyzed by the full-cycle rRNA approach. Systematic and Applied Microbiology. 2002;25:84–99. [PubMed]

36. LaPara TM, Nakatsu CH, Pantea L, Alleman JE. Phylogenetic analysis of bacterial communities in mesophilic and thermophilic bioreactors treating pharmaceutical wastewater. Appl Environ Microbiol. 2000;66:3951–3959. [PMC free article] [PubMed]

37. Roest K, Heilig H, Smidt H, de Vos WM, Stams AJM, et al. Community analysis of a full-scale anaerobic bioreactor treating paper mill wastewater. Systematic and Applied Microbiology. 2005;28:175–185. [PubMed]

38. Acinas SG, Klepac-Ceraj V, Hunt DE, Pharino C, Ceraj I, et al. Fine-scale phylogenetic architecture of a complex bacterial community. Nature. 2004;430:551–554. [PubMed]

39. Benlloch S, RodriguezValera F, MartinezMurcia AJ. Bacterial diversity in two coastal lagoons deduced from 16S rDNA PCR amplification and partial sequencing. FEMS Microbiol Ecol. 1995;18:267–279.

40. Daffonchio D, Borin S, Brusa T, Brusetti L, van der Wielen P, et al. Stratified prokaryote network in the oxic-anoxic transition of a deep-sea halocline. Nature. 2006;440:203–207. [PubMed]

41. DeLong EF, Preston CM, Mincer T, Rich V, Hallam SJ, et al. Community genomics among stratified microbial assemblages in the ocean's interior. Science. 2006;311:496–503. [PubMed]

42. Franklin MP, McDonald IR, Bourne DG, Owens NJP, Upstill-Goddard RC, et al. Bacterial diversity in the bacterioneuston (sea surface microlayer): the bacterioneuston through the looking glass. Environ Microbiology. 2005;7:723–736. [PubMed]

43. Fuhrman JA, McCallum K, Davis AA. Phylogenetic diversity of subsurface marine microbial communities from the Atlantic and Pacific oceans. Appl Environ Microbiol. 1993;59:1294–1302. [PMC free article] [PubMed]

44. Glockner FO, Zaichikov E, Belkova N, Denissova L, Pernthaler J, et al. Comparative 16S rRNA analysis of lake bacterioplankton reveals globally distributed phylogenetic clusters including an abundant group of actinobacteria. Appl Environ Microbiol. 2000;66:5053–+. [PMC free article] [PubMed]

45. Humayoun SB, Bano N, Hollibaugh JT. Depth distribution of microbial diversity in Mono Lake, a meromictic soda lake in California. Appl Environ Microbiol. 2003;69:1030–1042. [PMC free article] [PubMed]

46. Madrid VM, Taylor GT, Scranton MI, Chistoserdov AY. Phylogenetic diversity of bacterial and archaeal communities in the anoxic zone of the Cariaco Basin. Appl Environ Microbiol. 2001;67:1663–1674. [PMC free article] [PubMed]

47. Mosier AC, Murray AE, Fritsen CH. Microbiota within the perennial ice cover of Lake Vida, Antarctica. FEMS Microb Ecol. 2007;59:274–288. [PubMed]

48. Urbach E, Vergin KL, Young L, Morse A, Larson GL, et al. Unusual bacterioplankton community structure in ultra-oligotrophic Crater Lake. Limnology and Oceanography. 2001;46:557–572.

49. van der Wielen P, Bolhuis H, Borin S, Daffonchio D, Corselli C, et al. The enigma of prokaryotic life in deep hypersaline anoxic basins. Science. 2005;307:121–123. [PubMed]

50. Williams MM, Domingo JWS, Meckes MC, Kelty CA, Rochon HS. Phylogenetic diversity of drinking water bacteria in a distribution system simulator. J Appl Microbiol. 2004;96:954–964. [PubMed]

51. Woebken D, Fuchs BA, Kuypers MAA, Amann R. Potential interactions of particle-associated anammox bacteria with bacterial and archaeal partners in the Namibian upwelling system. Appl Environ Microbiol. 2007;73:4648–4657. [PMC free article] [PubMed]

52. Zaballos M, Lopez-Lopez A, Ovreas L, Bartual SG, D'Auria G, et al. Comparison of prokaryotic diversity at offshore oceanic locations. FEMS Microbiol Ecol. 2006;56:389–405. [PubMed]

53. Zrafi-Nouira I, Guermazi S, Chouari R, Safi NMD, Pelletier E, et al. Molecular diversity analysis and bacterial population dynamics of an adapted seawater microbiota during the degradation of Tunisian zarzatine oil. Biodegradation. 2009;20:467–486. [PubMed]

54. Ivanov II, Atarashi K, Manel N, Brodie EL, Shima T, et al. Induction of Intestinal Th17 Cells by Segmented Filamentous Bacteria. Cell. 2009;139:485–498. [PMC free article] [PubMed]

55. DeSantis TZ, Brodie EL, Moberg JP, Zubieta IX, Piceno YM, et al. High-density universal 16S rRNA microarray analysis reveals broader diversity than typical clone library when sampling the environment. Microb Ecol. 2007;53:371–383. [PubMed]

56. Chauhan A, Cherrier J, Williams HN. Impact of sideways and bottom-up control factors on bacterial community succession over a tidal cycle. Proc Nat Acad Sci. 2009;106:4301–4306. [PubMed]

57. Crump BC, Hopkinson CS, Sogin ML, Hobbie JE. Microbial biogeography along an estuarine salinity gradient: Combined influences of bacterial growth and residence time. App Environ Microbiol. 2004;70:1494–1505. [PMC free article] [PubMed]

58. Kaartokallio H, Laamanen M, Sivonen K. Responses of Baltic Sea ice and open-water natural bacterial communities to salinity change. Appl Environ Microbiol. 2005;71:4364–4371. [PMC free article] [PubMed]

59. Nelson CE. Phenology of high-elevation pelagic bacteria: the roles of meteorologic variability, catchment inputs and thermal stratification in structuring communities. ISME J. 2009;3:13–30. [PubMed]

60. Bodtker G, Thorstenson T, Lillebo BLP, Thorbjornsen BE, Ulvoen RH, et al. The effect of long-term nitrate treatment on SRB activity, corrosion rate and bacterial community composition in offshore water injection systems. J Ind Microbiol Biotechnol. 2008;35:1625–1636. [PubMed]

61. Hirayama H, Takai K, Inagaki F, Yamato Y, Suzuki M, et al. Bacterial community shift along a subsurface geothermal water stream in a Japanese gold mine. Extremophiles. 2005;9:169–184. [PubMed]

62. Wassel RA, Mills AL. Changes in water and sediment bacterial community structure in a lake receiving acid-mine drainage. Microbial Ecology. 1983;9:155–169.

63. Lozupone CA, Knight R. Global patterns in bacterial diversity. Proc Nat Acad Sci. 2007;104:11436–11440. [PubMed]

64. Hawkins RJ, Purdy KJ. 2007. pp. 187–194. Genotypic distribution of an indigenous model microorganism along an estuarine gradient. [PubMed]

65. Troussellier M, Schafer H, Batailler N, Bernard L, Courties C, et al. Bacterial activity and genetic richness along an estuarine gradient (Rhone River plume, France). Aquatic Microbial Ecology. 2002;28:13–24.

66. Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L, et al. Diversity of the human intestinal microbial flora. Science. 2005;308:1635–1638. [PMC free article] [PubMed]

67. Frank DN, Amand ALS, Feldman RA, Boedeker EC, Harpaz N, et al. Molecular-phylogenetic characterization of microbial community imbalances in human inflammatory bowel diseases. Proc Nat Acad Sci. 2007;104:13780–13785. [PubMed]

68. Mahowald MA, Rey FE, Seedorf H, Turnbaugh PJ, Fulton RS, et al. Characterizing a model human gut microbiota composed of members of its two dominant bacterial phyla. Proc Nat Acad Sci. 2009;106:5859–5864. [PubMed]

69. Rajilic-Stojanovic M, Heilig H, Molenaar D, Kajander K, Surakka A, et al. Development and application of the human intestinal tract chip, a phylogenetic microarray: analysis of universally conserved phylotypes in the abundant microbiota of young and elderly adults. Environ Microbiol. 2009;11:1736–1751. [PMC free article] [PubMed]

70. Scupham AJ, Patton TG, Bent E, Bayles DO. Comparison of the cecal microbiota of domestic and wild turkeys. Microbial Ecology. 2008;56:322–331. [PubMed]

71. Dowd SE, Callaway TR, Wolcott RD, Sun Y, McKeehan T, et al. Evaluation of the bacterial diversity in the feces of cattle using 16S rDNA bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP). BMC Microbiol. 2008;8 [PMC free article] [PubMed]

72. Cotta M, Forster R. The family Lachnospiraceae, including the genera Butyrivibrio, Lachnospira and Roseburia. In: Dworkin M, Schleifer K-H, Rosenberg E, Falkow S, editors. The Prokaryotes - A handbook on the biology of Bacteria: Bacteria: Firmicutes, Cyanobacteria. 3 ed. New York, NY: Springer-Verlag New York, LLC; 2006. pp. 1002–1021.

73. Ishii S, Ksoll WB, Hicks RE, Sadowsky MJ. Presence and growth of naturalized Escherichia coli in temperate soils from lake superior watersheds. Appl Environ Microbiol. 2006;72:612–621. [PMC free article] [PubMed]

74. Fogarty LR, Voytek MA. Comparison of Bacteroides-Prevotella 16S rRNA genetic markers for fecal samples from different animal species. Appl Environl Microbiol. 2005;71:5999–6007. [PMC free article] [PubMed]

75. Juteau P, Tremblay D, Villemur R, Bisaillon JG, Beaudet R. Analysis of the bacterial community inhabiting an aerobic thermophilic sequencing batch reactor (AT-SBR) treating swine waste. Appl Microbiol Biotechnol. 2005;66:115–122. [PubMed]

76. Lu JR, Domingo JS, Shanks OC. Identification of chicken-specific fecal microbial sequences using a metagenomic approach. Water Res. 2007;41:3561–3574. [PubMed]

77. Dang HY, Li TG, Chen MN, Huang GQ. Cross-Ocean distribution of Rhodobacterales bacteria as primary surface colonizers in temperate coastal marine waters. Appl Environ Microbiol. 2008;74:52–60. [PMC free article] [PubMed]

78. Geldreich EE. Fecal coliform and fecal streptococcus density relationships in waste discharges and receiving waters. CRC Crit Rev Environ Contr. 1976;6:349–369.

79. Black LE, Brion GM, Freitas SJ. Multivariate logistic regression for predicting total culturable virus presence at the intake of a potable-water treatment plant: Novel application of the atypical coliform/total coliform ratio. Appl Environ Microbiol. 2007;73:3965–3974. [PMC free article] [PubMed]

80. Nieman J, Brion GM. 2003. pp. 45–49. Novel bacterial ratio for predicting faecal age. [PubMed]

81. Smit E, Leeflang P, Gommans S, van den Broek J, van Mil S, et al. Diversity and seasonal fluctuations of the dominant members of the bacterial soil community in a wheat field as determined by cultivation and molecular methods. Appl Environ Microbiol. 2001;67:2284–2291. [PMC free article] [PubMed]

82. Sun HY, Deng SP, Raun WR. Bacterial community structure and diversity in a century-old manure-treated agroecosystem. Appl Environ Microbiol. 2004;70:5868–5874. [PMC free article] [PubMed]

83. Hartman WH, Richardson CJ, Vilgalys R, Bruland GL. Environmental and anthropogenic controls over bacterial communities in wetland soils. Proc Nat Acad Sci. 2008;105:17842–17847. [PubMed]

84. Paerl HW, Dyble J, Moisander PH, Noble RT, Piehler MF, et al. Microbial indicators of aquatic ecosystem change: current applications to eutrophication studies. FEMS Microbiol Ecol. 2003;46:233–246. [PubMed]

85. Graneli E, Weberg M, Salomon PS. Harmful algal blooms of allelopathic microalgal species: The role of eutrophication. Harmful Algae. 2008;8:94–102.

86. Moustafa A, Loram JE, Hackett JD, Anderson DM, Plumley FG, et al. Origin of Saxitoxin Biosynthetic Genes in Cyanobacteria. Plos One. 2009;4:1932–6203. [PMC free article] [PubMed]

87. Wilson KH, Blitchington RB, Greene RC. Amplification of Bacterial 16S ribosomal DNA with polymerase chain reaction. App Environ Microbiol. 1990;28:1942–1946. [PMC free article] [PubMed]

88. Masuda N, Church GM. Escherichia coli gene expression responsive to levels of the response regulator EvgA. J Bacteriol. 2002;184:6225–6234. [PMC free article] [PubMed]

89. Mei R, Hubbell E, Bekiranov S, Mittmann M, Christians FC, et al. Probe selection for high-density oligonucleotide arrays. Proc Nat Acad Sci. 2003;100:11237–11242. [PubMed]

90. Team RCD. 2008. R: A language and environment for statistical computing.

91. Goslee SC, Urban DL. The ecodist package for dissimilarity-based analysis of ecological data. J Stat Softw. 2007;22:1–19.

92. Evans J, Sheneman L, Foster JA. Relaxed neighbor-joining: a fast distance-based phylogenetic tree. Construction method. J Mol Evol. 2006;62:785–792. [PubMed]

93. McCune B, Mefford MJ. PC-ORD. Multivariate Analysis of Ecological Data. Version 5. Gleneden Beach, Oregon, USA: MjM Software; 2006.

94. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–5072. [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of
Public Library of Science