Microdiversity characterizes prevalent phylogenetic clades in the glacier-fed stream microbiome

Glacier-fed streams (GFSs) are extreme and rapidly vanishing ecosystems, and yet they harbor diverse microbial communities. Although our understanding of the GFS microbiome has recently increased, we do not know which microbial clades are ecologically successful in these ecosystems, nor do we understand potentially underlying mechanisms. Ecologically successful clades should be more prevalent across GFSs compared to other clades, which should be reflected as clade-wise distinctly low phylogenetic turnover. However, methods to assess such patterns are currently missing. Here we developed and applied a novel analytical framework, “phyloscore analysis”, to identify clades with lower spatial phylogenetic turnover than other clades in the sediment microbiome across twenty GFSs in New Zealand. These clades constituted up to 44% and 64% of community α-diversity and abundance, respectively. Furthermore, both their α-diversity and abundance increased as sediment chlorophyll a decreased, corroborating their ecological success in GFS habitats largely devoid of primary production. These clades also contained elevated levels of putative microdiversity than others, which could potentially explain their high prevalence in GFSs. This hitherto unknown microdiversity may be threatened as glaciers shrink, urging towards further genomic and functional exploration of the GFS microbiome.Subject terms: Climate-change ecology, Microbial ecology, Community ecology, Theoretical ecology     

Macroscopic Streamer Growths in Acidic, Metal-Rich Mine Waters in North Wales Consist of Novel and Remarkably Simple Bacterial Communities

The microbial composition of acid streamers (macroscopic biofilms) in acidic, metal-rich waters in two locations (an abandoned copper mine and a chalybeate spa) in north Wales was studied using cultivation-based and biomolecular techniques. Known chemolithotrophic and heterotrophic acidophiles were readily isolated from disrupted streamers, but they accounted for only <1 to 7% of the total microorganisms present. Fluorescent in situ hybridization (FISH) revealed that 80 to 90% of the microbes in both types of streamers were β-Proteobacteria. Terminal restriction fragment length polymorphism analysis of the streamers suggested that a single bacterial species was dominant in the copper mine streamers, while two distinct bacteria (one of which was identical to the bacterium found in the copper mine streamers) accounted for about 90% of the streamers in the spa water. 16S rRNA gene clone libraries showed that the β-proteobacterium found in both locations was closely related to a clone detected previously in acid mine drainage in California and that its closest characterized relatives were neutrophilic ammonium oxidizers. Using a modified isolation technique, this bacterium was isolated from the copper mine streamers and shown to be a novel acidophilic autotrophic iron oxidizer. The β-proteobacterium found only in the spa streamers was closely related to the neutrophilic iron oxidizer Gallionella ferruginea. FISH analysis using oligonucleotide probes that targeted the two β-proteobacteria confirmed that the biodiversity of the streamers in both locations was very limited. The microbial compositions of the acid streamers found at the two north Wales sites are very different from the microbial compositions of the previously described acid streamers found at Iron Mountain, California, and the Rio Tinto, Spain.     

An all-taxon microbial inventory of the Moorea coral reef ecosystem

The Moorea Coral Reef Long Term Ecological Research (LTER) Site (17.50°S, 149.83°W) comprises the fringe of coral reefs and lagoons surrounding the volcanic island of Moorea in the Society Islands of French Polynesia. As part of our Microbial Inventory Research Across Diverse Aquatic LTERS biodiversity inventory project, we characterized microbial community composition across all three domains of life using amplicon pyrosequencing of the V6 (bacterial and archaeal) and V9 (eukaryotic) hypervariable regions of small-subunit ribosomal RNA genes. Our survey spanned eight locations along a 130-km transect from the reef lagoon to the open ocean to examine changes in communities along inshore to offshore gradients. Our results illustrate consistent community differentiation between inshore and offshore ecosystems across all three domains, with greater richness in all domains in the reef-associated habitats. Bacterial communities were more homogenous among open ocean sites spanning >100 km than among inshore sites separated by <1 km, whereas eukaryotic communities varied more offshore than inshore, and archaea showed more equal levels of dissimilarity among subhabitats. We identified signature communities representative of specific geographic and geochemical milieu, and characterized co-occurrence patterns of specific microbial taxa within the inshore ecosystem including several bacterial groups that persist in geographical niches across time. Bacterial and archaeal communities were dominated by few abundant taxa but spatial patterning was consistent through time and space in both rare and abundant communities. This is the first in-depth inventory analysis of biogeographic variation of all three microbial domains within a coral reef ecosystem.     

Microbial biodiversity in glacier-fed streams

While glaciers become increasingly recognised as a habitat for diverse and active microbial communities, effects of their climate change-induced retreat on the microbial ecology of glacier-fed streams remain elusive. Understanding the effect of climate change on microorganisms in these ecosystems is crucial given that microbial biofilms control numerous stream ecosystem processes with potential implications for downstream biodiversity and biogeochemistry. Here, using a space-for-time substitution approach across 26 Alpine glaciers, we show how microbial community composition and diversity, based on 454-pyrosequencing of the 16S rRNA gene, in biofilms of glacier-fed streams may change as glaciers recede. Variations in streamwater geochemistry correlated with biofilm community composition, even at the phylum level. The most dominant phyla detected in glacial habitats were Proteobacteria, Bacteroidetes, Actinobacteria and Cyanobacteria/chloroplasts. Microorganisms from ice had the lowest α diversity and contributed marginally to biofilm and streamwater community composition. Rather, streamwater apparently collected microorganisms from various glacial and non-glacial sources forming the upstream metacommunity, thereby achieving the highest α diversity. Biofilms in the glacier-fed streams had intermediate α diversity and species sorting by local environmental conditions likely shaped their community composition. α diversity of streamwater and biofilm communities decreased with elevation, possibly reflecting less diverse sources of microorganisms upstream in the catchment. In contrast, β diversity of biofilms decreased with increasing streamwater temperature, suggesting that glacier retreat may contribute to the homogenisation of microbial communities among glacier-fed streams.     

Functional genomic analysis of an uncultured δ-proteobacterium in the sponge Cymbastela concentrica

Marine sponges are ancient, sessile, filter-feeding metazoans, which represent a significant component of the benthic communities throughout the world. Sponges harbor a remarkable diversity of bacteria, however, little is known about the functional properties of such bacterial symbionts. In this study, we present the genomic and functional characterization of an uncultured δ-proteobacterium associated with the sponge Cymbastela concentrica. We show that this organism represents a novel phylogenetic clade and propose that it lives in association with a cyanobacterium. We also provide an overview of the predicted functional and ecological properties of this δ-proteobacterium, and discuss its complex interactions with surrounding cells and milieu, including traits of cell attachment, nutrient transport and protein–protein interactions.     

Evidence of Physiological Decoupling from Grassland Ecosystem Drivers by an Encroaching Woody Shrub

Shrub encroachment of grasslands is a transformative ecological process by which native woody species increase in cover and frequency and replace the herbaceous community. Mechanisms of encroachment are typically assessed using temporal data or experimental manipulations, with few large spatial assessments of shrub physiology. In a mesic grassland in North America, we measured inter- and intra-annual variability in leaf δ¹³C in Cornus drummondii across a grassland landscape with varying fire frequency, presence of large grazers and topographic variability. This assessment of changes in individual shrub physiology is the largest spatial and temporal assessment recorded to date. Despite a doubling of annual rainfall (in 2008 versus 2011), leaf δ¹³C was statistically similar among and within years from 2008-11 (range of −28 to −27‰). A topography*grazing interaction was present, with higher leaf δ¹³C in locations that typically have more bare soil and higher sensible heat in the growing season (upland topographic positions and grazed grasslands). Leaf δ¹³C from slopes varied among grazing contrasts, with upland and slope leaf δ¹³C more similar in ungrazed locations, while slopes and lowlands were more similar in grazed locations. In 2011, canopy greenness (normalized difference vegetation index – NDVI) was assessed at the centroid of individual shrubs using high-resolution hyperspectral imagery. Canopy greenness was highest mid-summer, likely reflecting temporal periods when C assimilation rates were highest. Similar to patterns seen in leaf δ¹³C, NDVI was highest in locations that typically experience lowest sensible heat (lowlands and ungrazed). The ability of Cornus drummondii to decouple leaf physiological responses from climate variability and fire frequency is a likely contributor to the increase in cover and frequency of this shrub species in mesic grassland and may be generalizable to other grasslands undergoing woody encroachment.     

Partitioning of beta‐diversity reveals distinct assembly mechanisms of plant and soil microbial communities in response to nitrogen enrichment

Nitrogen (N) deposition poses a serious threat to terrestrial biodiversity and alters plant and soil microbial community composition. Species turnover and nestedness reflect the underlying mechanisms of variations in community composition. However, it remains unclear how species turnover and nestedness contribute to different responses of taxonomic groups (plants and soil microbes) to N enrichment. Here, based on a 13‐year consecutive multi‐level N addition experiment in a semiarid steppe, we partitioned community β‐diversity into species turnover and nestedness components and explored how and why plant and microbial communities reorganize via these two processes following N enrichment. We found that plant, soil bacterial, and fungal β‐diversity increased, but their two components showed different patterns with increasing N input. Plant β‐diversity was mainly driven by species turnover under lower N input but by nestedness under higher N input, which may be due to a reduction in forb species, with low tolerance to soil Mn²⁺, with increasing N input. However, turnover was the main contributor to differences in soil bacterial and fungal communities with increasing N input, indicating the phenomenon of microbial taxa replacement. The turnover of bacteria increased greatly whereas that of fungi remained within a narrow range with increasing N input. We further found that the increased soil Mn²⁺ concentration was the best predictor for increasing nestedness of plant communities under higher N input, whereas increasing N availability and acidification together contributed to the turnover of bacterial communities. However, environmental factors could explain neither fungal turnover nor nestedness. Our findings reflect two different pathways of community changes in plants, soil bacteria, and fungi, as well as their distinct community assembly in response to N enrichment. Disentangling the turnover and nestedness of plant and microbial β‐diversity would have important implications for understanding plant–soil microbe interactions and seeking conservation strategies for maintaining regional diversity.     

Anaerobic oxidation of methane at different temperature regimes in Guaymas Basin hydrothermal sediments

Anaerobic oxidation of methane (AOM) was investigated in hydrothermal sediments of Guaymas Basin based on δ¹³C signatures of CH₄, dissolved inorganic carbon and porewater concentration profiles of CH₄ and sulfate. Cool, warm and hot in-situ temperature regimes (15–20 °C, 30–35 °C and 70–95 °C) were selected from hydrothermal locations in Guaymas Basin to compare AOM geochemistry and 16S ribosomal RNA (rRNA), mcrA and dsrAB genes of the microbial communities. 16S rRNA gene clone libraries from the cool and hot AOM cores yielded similar archaeal types such as Miscellaneous Crenarchaeotal Group, Thermoproteales and anaerobic methane-oxidizing archaea (ANME)-1; some of the ANME-1 archaea formed a separate 16S rRNA lineage that at present seems to be limited to Guaymas Basin. Congruent results were obtained by mcrA gene analysis. The warm AOM core, chemically distinct by lower porewater sulfide concentrations, hosted a different archaeal community dominated by the two deep subsurface archaeal lineages Marine Benthic Group D and Marine Benthic Group B, and by members of the Methanosarcinales including ANME-2 archaea. This distinct composition of the methane-cycling archaeal community in the warm AOM core was confirmed by mcrA gene analysis. Functional genes of sulfate-reducing bacteria and archaea, dsrAB, showed more overlap between all cores, regardless of the core temperature. 16S rRNA gene clone libraries with Euryarchaeota-specific primers detected members of the Archaeoglobus clade in the cool and hot cores. A V6-tag high-throughput sequencing survey generally supported the clone library results while providing high-resolution detail on archaeal and bacterial community structure. These results indicate that AOM and the responsible archaeal communities persist over a wide temperature range.     

Linking microbial community structure to β-glucosidic function in soil aggregates

To link microbial community 16S structure to a measured function in a natural soil, we have scaled both DNA and β-glucosidase assays down to a volume of soil that may approach a unique microbial community. β-Glucosidase activity was assayed in 450 individual aggregates, which were then sorted into classes of high or low activities, from which groups of 10 or 11 aggregates were identified and grouped for DNA extraction and pyrosequencing. Tandem assays of ATP were conducted for each aggregate in order to normalize these small groups of aggregates for biomass size. In spite of there being no significant differences in the richness or diversity of the microbial communities associated with high β-glucosidase activities compared with the communities associated with low β-glucosidase communities, several analyses of variance clearly show that the communities of these two groups differ. The separation of these groups is partially driven by the differential abundances of members of the Chitinophagaceae family. It may be observed that functional differences in otherwise similar soil aggregates can be largely attributed to differences in resource availability, rather than to the presence or absence of particular taxonomic groups.     

Effect of Primers Hybridizing to Different Evolutionarily Conserved Regions of the Small-Subunit rRNA Gene in PCR-Based Microbial Community Analyses and Genetic Profiling

Genetic profiling techniques of microbial communities based on PCR-amplified signature genes, such as denaturing gradient gel electrophoresis or single-strand-conformation polymorphism (SSCP) analysis, are normally done with PCR products of less than 500-bp. The most common target for diversity analysis, the small-subunit rRNA genes, however, are larger, and thus, only partial sequences can be analyzed. Here, we compared the results obtained by PCR targeting different variable (V) regions (V2 and V3, V4 and V5, and V6 to V8) of the bacterial 16S rRNA gene with primers hybridizing to evolutionarily conserved flanking regions. SSCP analysis of single-stranded PCR products generated from 13 different bacterial species showed fewer bands with products containing V4-V5 (average, 1.7 bands per organism) than with V2-V3 (2.2 bands) and V6-V8 (2.3 bands). We found that the additional bands (>1 per organism) were caused by intraspecies operon heterogeneities or by more than one conformation of the same sequence. Community profiles, generated by PCR-SSCP from bacterial-cell consortia extracted from rhizospheres of field-grown maize (Zea mays), were analyzed by cloning and sequencing of the dominant bands. A total of 48 sequences could be attributed to 34 different strains from 10 taxonomical groups. Independent of the primer pairs, we found proteobacteria (α, β, and γ subgroups) and members of the genus Paenibacillus (low G+C gram-positive) to be the dominant organisms. Other groups, however, were only detected with single primer pairs. This study gives an example of how much the selection of different variable regions combined with different specificities of the flanking “universal” primers can affect a PCR-based microbial community analysis.     

Seasonal Changes in an Alpine Soil Bacterial Community in the Colorado Rocky Mountains

The period when the snowpack melts in late spring is a dynamic time for alpine ecosystems. The large winter microbial community begins to turn over rapidly, releasing nutrients to plants. Past studies have shown that the soil microbial community in alpine dry meadows of the Colorado Rocky Mountains changes in biomass, function, broad-level structure, and fungal diversity between winter and early summer. However, little specific information exists on the diversity of the alpine bacterial community or how it changes during this ecologically important period. We constructed clone libraries of 16S ribosomal DNA from alpine soil collected in winter, spring, and summer. We also cultivated bacteria from the alpine soil and measured the seasonal abundance of selected cultured isolates in hybridization experiments. The uncultured bacterial communities changed between seasons in diversity and abundance within taxa. The Acidobacterium division was most abundant in the spring. The winter community had the highest proportion of Actinobacteria and members of the Cytophaga/Flexibacter/Bacteroides (CFB) division. The summer community had the highest proportion of the Verrucomicrobium division and of β-Proteobacteria. As a whole, α-Proteobacteria were equally abundant in all seasons, although seasonal changes may have occurred within this group. A number of sequences from currently uncultivated divisions were found, including two novel candidate divisions. The cultured isolates belonged to the α-, β-, and γ-Proteobacteria, the Actinobacteria, and the CFB groups. The only uncultured sequences that were closely related to the isolates were from winter and spring libraries. Hybridization experiments showed that actinobacterial and β-proteobacterial isolates were most abundant during winter, while the α- and γ-proteobacterial isolates tested did not vary significantly. While the cultures and clone libraries produced generally distinct groups of organisms, the two approaches gave consistent accounts of seasonal changes in microbial diversity.     

Metagenomic Profiling of Microbial Composition and Antibiotic Resistance Determinants in Puget Sound

Human-health relevant impacts on marine ecosystems are increasing on both spatial and temporal scales. Traditional indicators for environmental health monitoring and microbial risk assessment have relied primarily on single species analyses and have provided only limited spatial and temporal information. More high-throughput, broad-scale approaches to evaluate these impacts are therefore needed to provide a platform for informing public health. This study uses shotgun metagenomics to survey the taxonomic composition and antibiotic resistance determinant content of surface water bacterial communities in the Puget Sound estuary. Metagenomic DNA was collected at six sites in Puget Sound in addition to one wastewater treatment plant (WWTP) that discharges into the Sound and pyrosequenced. A total of ∼550 Mbp (1.4 million reads) were obtained, 22 Mbp of which could be assembled into contigs. While the taxonomic and resistance determinant profiles across the open Sound samples were similar, unique signatures were identified when comparing these profiles across the open Sound, a nearshore marina and WWTP effluent. The open Sound was dominated by α-Proteobacteria (in particular Rhodobacterales sp.), γ-Proteobacteria and Bacteroidetes while the marina and effluent had increased abundances of Actinobacteria, β-Proteobacteria and Firmicutes. There was a significant increase in the antibiotic resistance gene signal from the open Sound to marina to WWTP effluent, suggestive of a potential link to human impacts. Mobile genetic elements associated with environmental and pathogenic bacteria were also differentially abundant across the samples. This study is the first comparative metagenomic survey of Puget Sound and provides baseline data for further assessments of community composition and antibiotic resistance determinants in the environment using next generation sequencing technologies. In addition, these genomic signals of potential human impact can be used to guide initial public health monitoring as well as more targeted and functionally-based investigations.     

Influence of an Oyster Reef on Development of the Microbial Heterotrophic Community of an Estuarine Biofilm

We characterized microbial biofilm communities developed over two very closely located but distinct benthic habitats in the Pensacola Bay estuary using two complementary cultivation-independent molecular techniques. Biofilms were grown for 7 days on glass slides held in racks 10 to 15 cm over an oyster reef and an adjacent muddy sand bottom. Total biomass and optical densities of dried biofilms showed dramatic differences for oyster reef versus non-oyster reef biofilms. This study assessed whether the observed spatial variation was reflected in the heterotrophic prokaryotic species composition. Genomic biofilm DNA from both locations was isolated and served as a template to amplify 16S rRNA genes with universal eubacterial primers. Fluorescently labeled PCR products were analyzed by terminal restriction fragment length polymorphism, creating a genetic fingerprint of the composition of the microbial communities. Unlabeled PCR products were cloned in order to construct a clone library of 16S rRNA genes. Amplified ribosomal DNA restriction analysis was used to screen and define ribotypes. Partial sequences from unique ribotypes were compared with existing database entries to identify species and to construct phylogenetic trees representative of community structures. A pronounced difference in species richness and evenness was observed at the two sites. The biofilm community structure from the oyster reef setting had greater evenness and species richness than the one from the muddy sand bottom. The vast majority of the bacteria in the oyster reef biofilm were related to members of the γ- and δ-subdivisions of Proteobacteria, the Cytophaga-Flavobacterium -Bacteroides cluster, and the phyla Planctomyces and Holophaga-Acidobacterium. The same groups were also present in the biofilm harvested at the muddy sand bottom, with the difference that nearly half of the community consisted of representatives of the Planctomyces phylum. Total species richness was estimated to be 417 for the oyster reef and 60 for the muddy sand bottom, with 10.5% of the total unique species identified being shared between habitats. The results suggest dramatic differences in habitat-specific microbial diversity that have implications for overall microbial diversity within estuaries.     

Spatial and Temporal Variation of Archaeal,Bacterial and Fungal Communities in Agricultural Soils

Background

Soil microbial communities are in constant change at many different temporal and spatial scales. However, the importance of these changes to the turnover of the soil microbial communities has been rarely studied simultaneously in space and time.

Methodology/Principal Findings

In this study, we explored the temporal and spatial responses of soil bacterial, archaeal and fungal β-diversities to abiotic parameters. Taking into account data from a 3-year sampling period, we analyzed the abundances and community structures of Archaea, Bacteria and Fungi along with key soil chemical parameters. We questioned how these abiotic variables influence the turnover of bacterial, archaeal and fungal communities and how they impact the long-term patterns of changes of the aforementioned soil communities. Interestingly, we found that the bacterial and fungal β-diversities are quite stable over time, whereas archaeal diversity showed significantly higher fluctuations. These fluctuations were reflected in temporal turnover caused by soil management through addition of N-fertilizers.

Conclusions

Our study showed that management practices applied to agricultural soils might not significantly affect the bacterial and fungal communities, but cause slow and long-term changes in the abundance and structure of the archaeal community. Moreover, the results suggest that, to different extents, abiotic and biotic factors determine the community assembly of archaeal, bacterial and fungal communities.     

Supragingival Microbial Profiles of Permanent and Deciduous Teeth in Children with Mixed Dentition

Objectives

The present study was designed to investigate the microbial profiles of teeth in different locations in mixed-dentition-stage children, and to compare the microbiomes of permanent and deciduous teeth in the same healthy oral cavity.

Methods

Supragingival plaque samples of teeth in various locations—the first permanent molars, deciduous molars, deciduous canines and incisors and permanent incisors—were collected from 20 healthy mixed-dentition-stage children with 10–12 permanent teeth erupted. Plaque DNA was extracted, and the V3–V4 hypervariable region of the bacterial 16S rRNA gene was amplified and subjected to sequencing.

Results

On average, 18,051 high-quality sequences per sample were generated. Permanent tooth sites tended to host more diverse bacterial communities than those of deciduous tooth sites. A total of 12 phyla, 21 classes, 38 orders, 66 families, 74 genera were detected ultimately. Five predominant phyla (Proteobacteria, Firmicutes, Bacteroidetes, Fusobacteria and Actinobacteria) were highly variable among sites. Of 26 genera with a mean relative abundance of >0.1%, 16 showed significant differences in relative abundance among the groups. More than 20% of the total operational taxonomical units were detected only in permanent or deciduous teeth. The variation in the microbial community composition was due mainly to permanent teeth being enriched in Actinomyces and deciduous teeth in Treponema. The core microbiome of supragingival plaque in mixed dentition comprised 19 genera with complex correlationships.

Conclusion

Our results suggest differences in microbial diversity and composition between permanent and deciduous teeth sites in mixed dentition. Moreover, the core microbiome of these sites was determined. These findings enhance our understanding of the development of the native oral microbiota with age.     

Stochastic Assembly of Bacteria in Microwell Arrays Reveals the Importance of Confinement in Community Development

The structure and function of microbial communities is deeply influenced by the physical and chemical architecture of the local microenvironment and the abundance of its community members. The complexity of this natural parameter space has made characterization of the key drivers of community development difficult. In order to facilitate these characterizations, we have developed a microwell platform designed to screen microbial growth and interactions across a wide variety of physical and initial conditions. Assembly of microbial communities into microwells was achieved using a novel biofabrication method that exploits well feature sizes for control of innoculum levels. Wells with incrementally smaller size features created populations with increasingly larger variations in inoculum levels. This allowed for reproducible growth measurement in large (20 μm diameter) wells, and screening for favorable growth conditions in small (5, 10 μm diameter) wells. We demonstrate the utility of this approach for screening and discovery using 5 μm wells to assemble P. aeruginosa colonies across a broad distribution of innoculum levels, and identify those conditions that promote the highest probability of survivial and growth under spatial confinement. Multi-member community assembly was also characterized to demonstrate the broad potential of this platform for studying the role of member abundance on microbial competition, mutualism and community succession.     

Reproducibility and quantitation of amplicon sequencing-based detection

To determine the reproducibility and quantitation of the amplicon sequencing-based detection approach for analyzing microbial community structure, a total of 24 microbial communities from a long-term global change experimental site were examined. Genomic DNA obtained from each community was used to amplify 16S rRNA genes with two or three barcode tags as technical replicates in the presence of a small quantity (0.1% wt/wt) of genomic DNA from Shewanella oneidensis MR-1 as the control. The technical reproducibility of the amplicon sequencing-based detection approach is quite low, with an average operational taxonomic unit (OTU) overlap of 17.2%±2.3% between two technical replicates, and 8.2%±2.3% among three technical replicates, which is most likely due to problems associated with random sampling processes. Such variations in technical replicates could have substantial effects on estimating β-diversity but less on α-diversity. A high variation was also observed in the control across different samples (for example, 66.7-fold for the forward primer), suggesting that the amplicon sequencing-based detection approach could not be quantitative. In addition, various strategies were examined to improve the comparability of amplicon sequencing data, such as increasing biological replicates, and removing singleton sequences and less-representative OTUs across biological replicates. Finally, as expected, various statistical analyses with preprocessed experimental data revealed clear differences in the composition and structure of microbial communities between warming and non-warming, or between clipping and non-clipping. Taken together, these results suggest that amplicon sequencing-based detection is useful in analyzing microbial community structure even though it is not reproducible and quantitative. However, great caution should be taken in experimental design and data interpretation when the amplicon sequencing-based detection approach is used for quantitative analysis of the β-diversity of microbial communities.     

Denitrifier Community Composition along a Nitrate and Salinity Gradient in a Coastal Aquifer

Nitrogen flux into the coastal environment via submarine groundwater discharge may be modulated by microbial processes such as denitrification, but the spatial scales at which microbial communities act and vary are not well understood. In this study, we examined the denitrifying community within the beach aquifer at Huntington Beach, California, where high-nitrate groundwater is a persistent feature. Nitrite reductase-encoding gene fragments (nirK and nirS), responsible for the key step in the denitrification pathway, were PCR amplified, cloned, and sequenced from DNAs extracted from aquifer sediments collected along a cross-shore transect, where groundwater ranged in salinity from 8 to 34 practical salinity units and in nitrate concentration from 0.5 to 330 μM. We found taxonomically rich and novel communities, with all nirK clones exhibiting <85% identity and nirS clones exhibiting <92% identity at the amino acid level to those of cultivated denitrifiers and other environmental clones in the database. Unique communities were found at each site, despite being located within 40 m of each other, suggesting that the spatial scale at which denitrifier diversity and community composition vary is small. Statistical analyses of nir sequences using the Monte Carlo-based program ∫-Libshuff confirmed that some populations were indeed distinct, although further sequencing would be required to fully characterize the highly diverse denitrifying communities at this site.     

Seasonal Dynamics of Shallow-Hyporheic-Zone Microbial Community Structure along a Heavy-Metal Contamination Gradient

Heavy metals contaminate numerous freshwater streams and rivers worldwide. Previous work by this group demonstrated a relationship between the structure of hyporheic microbial communities and the fluvial deposition of heavy metals along a contamination gradient during the fall season. Seasonal variation has been documented in microbial communities in numerous terrestrial and aquatic environments, including the hyporheic zone. The current study was designed to assess whether relationships between hyporheic microbial community structure and heavy-metal contamination vary seasonally by monitoring community structure along a heavy-metal contamination gradient for more than a year. No relationship between total bacterial abundance and heavy metals was observed (R² = 0.02, P = 0.83). However, denaturing gradient gel electrophoresis pattern analysis indicated a strong and consistent linear relationship between the difference in microbial community composition (populations present) and the difference in the heavy metal content of hyporheic sediments throughout the year (R² = 0.58, P < 0.001). Correlations between heavy-metal contamination and the abundance of four specific phylogenetic groups (most closely related to the α, β, and γ-proteobacteria and cyanobacteria) were apparent only during the fall and early winter, when the majority of organic matter is deposited into regional streams. These seasonal data suggest that the abundance of susceptible populations responds to heavy metals primarily during seasons when the potential for growth is highest.