Serial analysis of rRNA genes and the unexpected dominance of rare members of microbial communities

The role of Nucleic Acid Mimics (NAMs) on FISH-based techniques and applications for microbial detection

Fluorescent in situ hybridization (FISH) is a powerful tool that for more than 30 years has allowed to detect and quantify microorganisms as well as to study their spatial distribution in three-dimensional structured environments such as biofilms. Throughout these years, FISH has been improved in order to face some of its earlier limitations and to adapt to new research objectives. One of these improvements is related to the emergence of Nucleic Acid Mimics (NAMs), which are now employed as alternatives to the DNA and RNA probes that have been classically used in FISH. NAMs such as peptide and locked nucleic acids (PNA and LNA) have provided enhanced sensitivity and specificity to the FISH technique, as well as higher flexibility in terms of applications. In this review, we aim to cover the state-of-the-art of the different NAMs and explore their possible applications in FISH, providing a general overview of the technique advancement in the last decades.

Assessing Biotic and Abiotic Interactions of Microorganisms in Amazonia through Co-Occurrence Networks and DNA Metabarcoding

Species may co-occur due to responses to similar environmental conditions, biological associations, or simply because of coincident geographical distributions. Disentangling patterns of co-occurrence and potential biotic and abiotic interactions is crucial to understand ecosystem function. Here, we used DNA metabarcoding data from litter and mineral soils collected from a longitudinal transect in Amazonia to explore patterns of co-occurrence. We compared data from different Amazonian habitat types, each with a characteristic biota and environmental conditions. These included non-flooded rainforests (terra-firme), forests seasonally flooded by fertile white waters (várzeas) or by unfertile black waters (igapós), and open areas associated with white sand soil (campinas). We ran co-occurrence network analyses based on null models and Spearman correlation for all samples and for each habitat separately. We found that one third of all operational taxonomic units (OTUs) were bacteria and two thirds were eukaryotes. The resulting networks were nevertheless mostly composed of bacteria, with fewer fungi, protists, and metazoans. Considering the functional traits of the OTUs, there is a combination of metabolism modes including respiration and fermentation for bacteria, and a high frequency of saprotrophic fungi (those that feed on dead organic matter), indicating a high turnover of organic material. The organic carbon and base saturation indices were important in the co-occurrences in Amazonian networks, whereas several other soil properties were important for the co-exclusion. Different habitats had similar network properties with some variation in terms of modularity, probably associated with flooding pulse. We show that Amazonian microorganism communities form highly interconnected co-occurrence and co-exclusion networks, which highlights the importance of complex biotic and abiotic interactions in explaining the outstanding biodiversity of the region.

A Tripartite Microbial-Environment Network Indicates How Crucial Microbes Influence the Microbial Community Ecology

Current technologies could identify the abundance and functions of specific microbes, and evaluate their individual effects on microbial ecology. However, these microbes interact with each other, as well as environmental factors, in the form of complex network. Determination of their combined ecological influences remains a challenge. In this study, we developed a tripartite microbial-environment network (TMEN) analysis method that integrates microbial abundance, metabolic function, and environmental data as a tripartite network to investigate the combined ecological effects of microbes. Applying TMEN to analyzing the microbial-environment community structure in the sediments of Hangzhou Bay, one of the most seriously polluted coastal areas in China, we found that microbes were well-organized into 4 bacterial communities and 9 archaeal communities. The total organic carbon, sulfate, chemical oxygen demand, salinity, and nitrogen-related indexes were detected as crucial environmental factors in the microbial-environmental network. With close interactions with these environmental factors, Nitrospirales and Methanimicrococcu were identified as hub microbes with connection advantage. Our TMEN method could close the gap between lack of efficient statistical and computational approaches and the booming of large-scale microbial genomic and environmental data. Based on TMEN, we discovered a potential microbial ecological mechanism that crucial species with significant influence on the microbial community ecology would possess one or two of the community advantages for enhancing their ecological status and essentiality, including abundance advantage and connection advantage.

First-tier detection of intragenomic 16S rRNA gene variation in culturable endophytic bacteria from cacao seeds

Background

Intragenomic variability in 16S rDNA is a limiting factor for taxonomic and diversity characterization of Bacteria, and studies on its occurrence in natural/environmental populations are scarce. In this work, direct DNA amplicon sequencing coupled with frequent-cutter restriction analysis allowed detection of intragenomic 16S rDNA variation in culturable endophytic bacteria from cacao seeds in a fast and attractive manner.

Methods

Total genomic DNA from 65 bacterial strains was extracted and the 16S rDNA hyper variable V5-V9 regions were amplified for enzyme digestion and direct Sanger-type sequencing. The resulting electropherograms were visually inspected and compared to the corresponding AluI-restriction profiles, as well as to complete genome sequences in databases. Restriction analysis were employed to substitute the need of amplicon cloning and re-sequencing. A specifically improved polyacrylamide-gradient electrophoresis allowed to resolve 5-bp differences in restriction fragment sizes. Chi-square analysis on 2 × 2 contingency table tested for the independence between the 'number of AluI bands' and 'type of eletropherogram'.

Results

Two types of electropherograms were obtained: unique template, with single peaks per base (clean chromatograms), and heterogeneous template, with various levels of multiple peaks per base (mixed chromatograms). Statistics revealed significant interaction between number of restriction fragments and type of electropherogram for the same amplicons: clean or mixed ones associated to ≤5 or ≥6 bands, respectively. The mixed-template pattern combined with the AluI-restriction profiles indicated a high proportion of 49% of the culturable endophytes from a tropical environment showing evidence of intragenomic 16S rDNA heterogeneity.

Conclusion

The approach presented here was useful for a rapid, first-tier detection of intragenomic variation in culturable isolates, which can be applied in studies of other natural populations; a preliminary view of intragenomic heterogeneity levels can complement culture-dependent and -independent methods. Consequences of these findings in taxonomic and diversity studies in complex bacterial communities are discussed.

Effect of increased load of high-strength food wastewater in thermophilic and mesophilic anaerobic co-digestion of waste activated sludge on bacterial community structure

In recent years, anaerobic co-digestion (AcoD) has been widely used to improve reactor performance, especially methane production. In this study, we applied two different operating temperatures (thermophilic and mesophilic) and gradually increased the load of food wastewater (FWW) to investigate the bacterial communities during the AcoD of waste activated sludge (WAS) and FWW. As the load of FWW was increased, methane production rate (MPR; L CH4/L d) and methane content (%) in both Thermophilic AcoD (TAcoD) and Mesophilic AcoD (MAcoD) increased significantly; the highest MPR and methane content in TAcoD (1.423 L CH4/L d and 68.24%) and MAcoD (1.233 L CH4/L d and 65.21%) were observed when the FWW mixing ratio was 75%. However, MPR and methane yield in both reactors decreased markedly and methane production in TAcoD ceased completely when only FWW was fed into the reactor, resulting from acidification of the reactor caused by accumulation of organic acids. Pyrosequencing analysis revealed a decrease in bacterial diversity in TAcoD and a markedly different composition of bacterial communities between TAcoD and MAcoD with an increase in FWW load. For example, Bacterial members belonging to two genera Petrotoga (assigned to phylum Thermotogae) and Petrimonas (assigned to phylum Bacteroidetes) became dominant in TAcoD and MAcoD with an increase in FWW load, respectively. In addition, quantitative real-time PCR (qPCR) results showed higher bacterial and archaeal populations (expressed as 16S rRNA gene concentration) in TAcoD than MAcoD with an increase in FWW load and showed maximum population when the FWW mixing ratio was 75% in both reactors. Collectively, this study demonstrated the dynamics of key bacterial communities in TAcoD and MAcoD, which were highly affected by the load of FWW.

Contrasting soil bacterial community structure between the phyla Acidobacteria and Proteobacteria in tropical Southeast Asian and temperate Japanese forests

Soil bacterial community structures of six dominant phyla (Acidobacteria, Proteobacteria, Verrucomicrobia, Planctomycetes, Bacteroidetes and Actinobacteria) and unclassified bacteria detected in tropical Sarawakian and temperate Japanese forests were compared based on 16S rRNA gene sequence variation. The class composition in each phylum was similar among the studied forests; however, significant heterogeneities of class frequencies were detected. Acidobacteria and Proteobacteria were the most dominant phyla in all six forests, but differed in the level of bacterial species diversity, pattern of species occurrence and association pattern of species composition with physicochemical properties in soil. Species diversity among Acidobacteria was approximately half that among Proteobacteria, based on the number of clusters and the Chao1 index, even though a similar number of sequence reads were obtained for these two phyla. In contrast, species diversity within Planctomycetes and Bacteroidetes was nearly as high as within Acidobacteria, despite many fewer sequence reads. The density of species (the number of sequence reads per cluster) correlated negatively with species diversity, and species density within Acidobacteria was approximately twice that within Proteobacteria. Although the percentage of forest-specific species was high for all bacterial groups, sampling site-specific species varied among bacterial groups, indicating limited inter-forest migration and differential movement of bacteria in forest soil. For five of the seven bacterial groups, including Acidobacteria, soil pH appeared to strongly influence species composition, but this association was not observed for Proteobacterial species. Topology of UPGMA trees and pattern of NMDS plots among the forests differed among the bacterial groups, suggesting that each bacterial group has adapted and evolved independently in each forest.

Acetic acid effects on methanogens in the second stage of a two-stage anaerobic system

This study reports on biomass tolerance towards high concentrations of acetic acid (HAc) within the system. Biomass from the second stage of a two-stage anaerobic sludge digestion system was used for this study. Microbial community analysis by 454 pyrosequencing highlighted hydrogenotrophic Methanomicrobiales was the predominant archaeal population in the second stage (>99% of the total archaeal community). Second stage biomass degraded HAc up to 4200 mg HAc L(-1) without observable lag phase. However, at HAc-shock loading of 7400 mg HAc L(-1), it showed a one day lag phase associated with decreased biomass activity. After stepwise HAc-acclimation over 27 d, the biomass degraded HAc of up to 8200 mg HAc L(-1) without observable lag phase. The dominance of Methanomicrobiales had remained unchanged in proportion - while the total archaeal population increased during acclimation. This study showed stepwise acclimation could be an approach to accommodate HAc accumulation and hence higher concentrations resulting from an enhanced first stage.

Dynamics of propionic acid degradation in a two-phase anaerobic system

This paper reports on propionic acid (HPr) degradation in a laboratory scale two-phase anaerobic system, where HPr was accumulated in the acidogenic reactor and degraded in the methanogenic reactor. Batch tests using biomass from the two-phase anaerobic system showed HPr degradation was rarely detectable in the acidogenic reactor when HPr concentration ranged from 639 to 4531mgHPrL(-1) and at pH 4.50 to 6.50. Biomass from the methanogenic reactor could, however, successfully degrade HPr at its initial concentration of up to 4585mgHPrL(-1) at pH 6.40-7.30. ATP results showed that differences in the degradation ability of HPr by the acidogenic and methanogenic biomass may be related with their respective different biomass activities. Results from pyrosequencing showed that the predominant propionic acid oxidizing bacteria (POB) in the methanogenic reactor were Smithella (2.68%) and Syntrophobacter (0.35%); while poor degradation of HPr in the acidogenic reactor may be associated with the low abundance of POB (0.02% Desulfacinum and 0.08% Desulfobulbus). This might have been induced by the long-term unfavorable environment for POB growth in the acidogenic reactor.

Efficient and Accurate OTU Clustering with GPU-Based Sequence Alignment and Dynamic Dendrogram Cutting

De novo clustering is a popular technique to perform taxonomic profiling of a microbial community by grouping 16S rRNA amplicon reads into operational taxonomic units (OTUs). In this work, we introduce a new dendrogram-based OTU clustering pipeline called CRiSPy. The key idea used in CRiSPy to improve clustering accuracy is the application of an anomaly detection technique to obtain a dynamic distance cutoff instead of using the de facto value of 97 percent sequence similarity as in most existing OTU clustering pipelines. This technique works by detecting an abrupt change in the merging heights of a dendrogram. To produce the output dendrograms, CRiSPy employs the OTU hierarchical clustering approach that is computed on a genetic distance matrix derived from an all-against-all read comparison by pairwise sequence alignment. However, most existing dendrogram-based tools have difficulty processing datasets larger than 10,000 unique reads due to high computational complexity. We address this difficulty by developing two efficient algorithms for CRiSPy: a compute-efficient GPU-accelerated parallel algorithm for pairwise distance matrix computation and a memory-efficient hierarchical clustering algorithm. Our experiments on various datasets with distinct attributes show that CRiSPy is able to produce more accurate OTU groupings than most OTU clustering applications.

Response of the rare biosphere to environmental stressors in a highly diverse ecosystem (Zodletone spring, OK, USA)

Within highly diverse ecosystems, the majority of bacterial taxa are present in low abundance as members of the rare biosphere. The rationale for the occurrence and maintenance of the rare biosphere, and the putative ecological role(s) and dynamics of its members within a specific ecosystem is currently debated. We hypothesized that in highly diverse ecosystems, a fraction of the rare biosphere acts as a backup system that readily responds to environmental disturbances. We tested this hypothesis by subjecting sediments from Zodletone spring, a sulfide- and sulfur-rich spring in Southwestern OK, to incremental levels of salinity (1, 2, 3, 4, and 10% NaCl), or temperature (28°, 30°, 32°, and 70 °C), and traced the trajectories of rare members of the community in response to these manipulations using 16S rRNA gene analysis. Our results indicate that multiple rare bacterial taxa are promoted from rare to abundant members of the community following such manipulations and that, in general, the magnitude of such recruitment is directly proportional to the severity of the applied manipulation. Rare members that are phylogenetically distinct from abundant taxa in the original sample (unique rare biosphere) played a more important role in the microbial community response to environmental disturbances, compared to rare members that are phylogenetically similar to abundant taxa in the original sample (non-unique rare biosphere). The results emphasize the dynamic nature of the rare biosphere, and highlight its complexity and non-monolithic nature.

Comprehensive microbial analysis of combined mesophilic anaerobic-thermophilic aerobic process treating high-strength food wastewater

A combined mesophilic anaerobic-thermophilic aerobic process was used to treat high-strength food wastewater in this study. During the experimental period, most of solid residue from the mesophilic anaerobic reactor (R1) was separated by centrifugation and introduced into the thermophilic aerobic reactor (R2) for further digestion. Then, thermophilic aerobically-digested sludge was reintroduced into R1 to enhance reactor performance. The combined process was operated with two different Runs: Run I with hydraulic retention time (HRT) = 40 d (corresponding OLR = 3.5 kg COD/m(3) d) and Run II with HRT = 20 d (corresponding OLR = 7 kg COD/m(3)). For a comparison, a single-stage mesophilic anaerobic reactor (R3) was operated concurrently with same OLRs and HRTs as the combined process. During the overall digestion, all reactors showed high stability without pH control. The combined process demonstrated significantly higher organic matter removal efficiencies (over 90%) of TS, VS and COD and methane production than did R3. Quantitative real-time PCR (qPCR) results indicated that higher populations of both bacteria and archaea were maintained in R1 than in R3. Pyrosequencing analysis revealed relatively high abundance of phylum Actinobacteria in both R1 and R2, and a predominance of phyla Synergistetes and Firmicutes in R3 during Run II. Furthermore, R1 and R2 shared genera (Prevotella, Aminobacterium, Geobacillus and Unclassified Actinobacteria), which suggests synergy between mesophilic anaerobic digestion and thermophilic aerobic digestion. For archaea, in R1 methanogenic archaea shifted from genus Methanosaeta to Methanosarcina, whereas genera Methanosaeta, Methanobacterium and Methanoculleus were predominant in R3. The results demonstrated dynamics of key microbial populations that were highly consistent with an enhanced reactor performance of the combined process.

Ecology and exploration of the rare biosphere

The profound influence of microorganisms on human life and global biogeochemical cycles underlines the value of studying the biogeography of microorganisms, exploring microbial genomes and expanding our understanding of most microbial species on Earth: that is, those present at low relative abundance. The detection and subsequent analysis of low-abundance microbial populations—the 'rare biosphere'—have demonstrated the persistence, population dynamics, dispersion and predation of these microbial species. We discuss the ecology of rare microbial populations, and highlight molecular and computational methods for targeting taxonomic 'blind spots' within the rare biosphere of complex microbial communities.

Detection of bacterial DNA in lymph nodes of Crohn's disease patients using high throughput sequencing

OBJECTIVE

Our aim was to determine whether or not specific microorganisms were transported selectively to lymph nodes in Crohn's disease (CD) by comparing node and mucosal microbial communities in patients and controls. We also sought evidence of dysbiosis and bacterial translocation.

DESIGN

Lymph nodes, and involved and uninvolved mucosal samples were obtained from resections of 58 patients (29 CD, eight 'other inflammatory bowel disease' (IBD) and 21 non-IBD). Universal primers targeting V1-V3 regions of bacterial 16S rRNA genes were used to amplify bacterial DNA and amplicons sequenced using high throughput sequencing. 20 patients (eight CD (28%), two other IBD (25%) and 10 non-IBD (48%)) had PCR positive nodes.

RESULTS

All samples from an individual were similar: there was no evidence of selective concentration of any microorganism in nodes. No specific microorganism was present in the nodes of all CD samples. Escherichia/Shigella were common in all patient groups but patients with ileal CD had a greater proportion of Escherichia coli reads in their nodes than other CD patients (p=0.0475). Campylobacter, Helicobacter and Yersinia were uncommon; Mycobacterium and Listeria were not detected. Dysbiosis was present in all groups but shifts were specific and no common pattern emerged.

CONCLUSIONS

It is unlikely that a single bacterium perpetuates inflammation in late stage CD; dysbiosis was common and we found no evidence of increased bacterial translocation. We believe that future studies should focus on early disease and viable bacteria in nodes, aphthous ulcers and granulomas, as they may be more relevant in the initiation of inflammation in CD.

Bacterial and methanogenic archaeal communities during the single-stage anaerobic digestion of high-strength food wastewater

Single-stage anaerobic digestion (AD) was operated to treat high-strength food wastewater (FWW) derived from food waste recycling facilities at two different organic loading rates (OLRs) of 3.5 (Phase I) and 7 (Phase II) kgCOD/m(3)d. Changes in composition of microbial communities were investigated using quantitative real-time PCR (qPCR) and barcoded-pyrosequencing. At the high FWW loading rate, AD showed efficient performance (i.e., organic matter removal and methane production). Bacterial communities were represented by the phyla Bacteroidetes, Firmicutes, Synergistetes and Actinobacteria. During the entire digestion process, the relative abundance phylum Chloroflexi decreased significantly. The qPCR analysis demonstrated that the methanogenic communities shifted from aceticlastic (Methanosarcinales) to hydrogenotrophic methanogens (Methanobacteriales and Methanomicrobiales) with high increase in the proportion of syntrophic bacterial communities. Canonical correspondence analysis revealed a strong relationship between reactor performance and microbial community shifts.

Airborne bacterial populations above desert soils of the McMurdo Dry Valleys, Antarctica

Bacteria are assumed to disperse widely via aerosolized transport due to their small size and resilience. The question of microbial endemicity in isolated populations is directly related to the level of airborne exogenous inputs, yet this has proven hard to identify. The ice-free terrestrial ecosystem of Antarctica, a geographically and climatically isolated continent, was used to interrogate microbial bio-aerosols in relation to the surrounding ecology and climate. High-throughput sequencing of bacterial ribosomal RNA (rRNA) genes was combined with analyses of climate patterns during an austral summer. In general terms, the aerosols were dominated by Firmicutes, whereas surrounding soils supported Actinobacteria-dominated communities. The most abundant taxa were also common to aerosols from other continents, suggesting that a distinct bio-aerosol community is widely dispersed. No evidence for significant marine input to bioaerosols was found at this maritime valley site, instead local influence was largely from nearby volcanic sources. Back trajectory analysis revealed transport of incoming regional air masses across the Antarctic Plateau, and this is envisaged as a strong selective force. It is postulated that local soil microbial dispersal occurs largely via stochastic mobilization of mineral soil particulates.

Proceedings of the 2013 A.S.P.E.N. Research workshop: the interface between nutrition and the gut microbiome: implications and applications for human health [corrected]

The human and earth microbiomes are among the most important biological agents in understanding and preventing disease. Technology is advancing at a fast pace and allowing for high-resolution analysis of the composition and function of our microbial partners across regions, space, and time. Bioinformaticists and biostatisticians are developing ever more elegant displays to understand the generated megadatasets. A virtual cyberinfrastructure of search engines to cross-reference the rapidly developing data is emerging in line with technologic advances. Nutrition science will reap the benefits of this new field, and its role in preserving the earth and the humans who inhabit it will become evidently clear. In this report we highlight some of the topics of an A.S.P.E.N.-sponsored symposium held during Clinical Nutrition Week in 2013 that address the importance of the human microbiome to human health and disease.

Diffuse flow environments within basalt- and sediment-based hydrothermal vent ecosystems harbor specialized microbial communities

Hydrothermal vents differ both in surface input and subsurface geochemistry. The effects of these differences on their microbial communities are not clear. Here, we investigated both alpha and beta diversity of diffuse flow-associated microbial communities emanating from vents at a basalt-based hydrothermal system along the East Pacific Rise (EPR) and a sediment-based hydrothermal system, Guaymas Basin. Both Bacteria and Archaea were targeted using high throughput 16S rRNA gene pyrosequencing analyses. A unique aspect of this study was the use of a universal set of 16S rRNA gene primers to characterize total and diffuse flow-specific microbial communities from varied deep-sea hydrothermal environments. Both surrounding seawater and diffuse flow water samples contained large numbers of Marine Group I (MGI) Thaumarchaea and Gammaproteobacteria taxa previously observed in deep-sea systems. However, these taxa were geographically distinct and segregated according to type of spreading center. Diffuse flow microbial community profiles were highly differentiated. In particular, EPR dominant diffuse flow taxa were most closely associated with chemolithoautotrophs, and off axis water was dominated by heterotrophic-related taxa, whereas the opposite was true for Guaymas Basin. The diversity and richness of diffuse flow-specific microbial communities were strongly correlated to the relative abundance of Epsilonproteobacteria, proximity to macrofauna, and hydrothermal system type. Archaeal diversity was higher than or equivalent to bacterial diversity in about one third of the samples. Most diffuse flow-specific communities were dominated by OTUs associated with Epsilonproteobacteria, but many of the Guaymas Basin diffuse flow samples were dominated by either OTUs within the Planctomycetes or hyperthermophilic Archaea. This study emphasizes the unique microbial communities associated with geochemically and geographically distinct hydrothermal diffuse flow environments.

Impact of colonoscopy bowel preparation on intestinal microbiota

The gut microbiota is important in maintaining human health, but numerous factors have the potential to alter its composition. Our aim was to examine the impact of a standard bowel preparation on the intestinal microbiota using two different techniques. Fifteen subjects undergoing colonoscopy consumed a bowel preparation comprised of 10 mg bisacodyl and 2 L polyethylene glycol. The microbiota of stool samples, collected one month before, one week before (pre-colonoscopy), and one week, one month, and three to six months after colonoscopy (post-colonoscopy) was evaluated. Two samples were taken three to six months apart from five healthy subjects who did not undergo colonoscopy. Universal primers targeting the V2-V3 region of the 16S rRNA gene were used to PCR amplify all samples for denaturing gradient gel electrophoresis (PCR-DGGE). Pre- and post-colonoscopy samples were compared using Dice's similarity coefficients. Three samples from ten subjects who underwent colonoscopy, and both samples from the five subjects who didn't, were used for high-throughput sequencing of the V1-V3 region of the 16S rRNA gene. Samples were curated and analysed in Mothur. Results of the DGGE analyses show that the fecal microbiota of a small number of subjects had short-term changes. High-throughput sequencing results indicated that the variation between the samples of subjects who underwent colonoscopy was no greater than the variation observed between samples from subjects who did not. We conclude that bowel preparation does not have a lasting effect on the composition of the intestinal microbiota for the majority of subjects.

Effects of OTU clustering and PCR artifacts on microbial diversity estimates

Next-generation sequencing has increased the coverage of microbial diversity surveys by orders of magnitude, but differentiating artifacts from rare environmental sequences remains a challenge. Clustering 16S rRNA sequences into operational taxonomic units (OTUs) organizes sequence data into groups of 97 % identity, helping to reduce data volumes and avoid analyzing sequencing artifacts by grouping them with real sequences. Here, we analyze sequence abundance distributions across environmental samples and show that 16S rRNA sequences of >99 % identity can represent functionally distinct microorganisms, rendering OTU clustering problematic when the goal is an accurate analysis of organism distribution. Strict postsequencing quality control (QC) filters eliminated the most prevalent artifacts without clustering. Further experiments proved that DNA polymerase errors in polymerase chain reaction (PCR) generate a significant number of substitution errors, most of which pass QC filters. Based on our findings, we recommend minimizing the number of PCR cycles in DNA library preparation and applying strict postsequencing QC filters to reduce the most prevalent artifacts while maintaining a high level of accuracy in diversity estimates. We further recommend correlating rare and abundant sequences across environmental samples, rather than clustering into OTUs, to identify remaining sequence artifacts without losing the resolution afforded by high-throughput sequencing.

Groundtruthing next-gen sequencing for microbial ecology-biases and errors in community structure estimates from PCR amplicon pyrosequencing

Analysis of microbial communities by high-throughput pyrosequencing of SSU rRNA gene PCR amplicons has transformed microbial ecology research and led to the observation that many communities contain a diverse assortment of rare taxa-a phenomenon termed the Rare Biosphere. Multiple studies have investigated the effect of pyrosequencing read quality on operational taxonomic unit (OTU) richness for contrived communities, yet there is limited information on the fidelity of community structure estimates obtained through this approach. Given that PCR biases are widely recognized, and further unknown biases may arise from the sequencing process itself, a priori assumptions about the neutrality of the data generation process are at best unvalidated. Furthermore, post-sequencing quality control algorithms have not been explicitly evaluated for the accuracy of recovered representative sequences and its impact on downstream analyses, reducing useful discussion on pyrosequencing reads to their diversity and abundances. Here we report on community structures and sequences recovered for in vitro-simulated communities consisting of twenty 16S rRNA gene clones tiered at known proportions. PCR amplicon libraries of the V3-V4 and V6 hypervariable regions from the in vitro-simulated communities were sequenced using the Roche 454 GS FLX Titanium platform. Commonly used quality control protocols resulted in the formation of OTUs with >1% abundance composed entirely of erroneous sequences, while over-aggressive clustering approaches obfuscated real, expected OTUs. The pyrosequencing process itself did not appear to impose significant biases on overall community structure estimates, although the detection limit for rare taxa may be affected by PCR amplicon size and quality control approach employed. Meanwhile, PCR biases associated with the initial amplicon generation may impose greater distortions in the observed community structure.

Microbialite genetic diversity and composition relate to environmental variables

Microbialites have played an important role in the early history of life on Earth. Their fossilized forms represent the oldest evidence of life on our planet dating back to 3500 Ma. Extant microbialites have been suggested to be highly productive and diverse communities with an evident role in the cycling of major elements, and in contributing to carbonate precipitation. Although their ecological and evolutionary importance has been recognized, the study of their genetic diversity is yet scanty. The main goal of this study was to analyse microbial genetic diversity of microbialites living in different types of environments throughout Mexico, including desert ponds, coastal lagoons and a crater-lake. We followed a pyrosequencing approach of hypervariable regions of the 16S rRNA gene. Results showed that microbialite communities were very diverse (H' = 6-7) and showed geographic variation in composition, as well as an environmental effect related to pH and conductivity, which together explained 33% of the genetic variation. All microbialites had similar proportions of major bacterial and archaeal phyla.

Modeling microbial community structure and functional diversity across time and space

Microbial communities exhibit exquisitely complex structure. Many aspects of this complexity, from the number of species to the total number of interactions, are currently very difficult to examine directly. However, extraordinary efforts are being made to make these systems accessible to scientific investigation. While recent advances in high-throughput sequencing technologies have improved accessibility to the taxonomic and functional diversity of complex communities, monitoring the dynamics of these systems over time and space - using appropriate experimental design - is still expensive. Fortunately, modeling can be used as a lens to focus low-resolution observations of community dynamics to enable mathematical abstractions of functional and taxonomic dynamics across space and time. Here, we review the approaches for modeling bacterial diversity at both the very large and the very small scales at which microbial systems interact with their environments. We show that modeling can help to connect biogeochemical processes to specific microbial metabolic pathways.

Quantitative fluorescence in situ hybridization analysis of microbial consortia from a biogenic gas field in Alaska's Cook Inlet basin

Filter-collected production water samples from a methane-rich gas field in the Cook Inlet basin of Alaska were investigated using whole-cell rRNA-targeted fluorescence in situ hybridization (FISH) and 16S rRNA tag pyrosequencing. Both techniques were consistent in determining the microbial community composition, including the archaeal or bacterial dominance of samples. The archaeal community is dominated by the obligate methylotrophic methanogen genus Methanolobus as well as the nutritional generalist methanogen genus Methanosarcina, which is capable of utilizing acetate, CO(2), and methyl-bearing compounds. The most-abundant bacterial groups are Firmicutes, notably of the Acetobacterium genus, and Cytophaga-Flexibacter-Bacteroides species (CFBs) affiliated with the order Bacteroidales. We observed spatial variation among samples in both the percentage of members of Archaea compared to that of members of Bacteria and the dominant members of the bacterial community, differences which could not be explained with the available geochemical data. Based upon the microbial community composition and the isotopic signature of methane associated with the Cook Inlet basin site, we propose a simplified reaction network beginning with the breakdown of coal macromolecules, followed by fermentation and methylotrophic and acetoclastic methane production.

CLASI-FISH: principles of combinatorial labeling and spectral imaging

Although the number of phylotypes present in a microbial community may number in the hundreds or more, until recently, fluorescence in situ hybridization has been used to label, at most, only a handful of different phylotypes in a single sample. We recently developed a technique, CLASI-FISH for combinatorial labeling and spectral imaging - fluorescence in situ hybridization, to greatly expand the number of distinguishable taxa in a single FISH experiment. The CLASI technique involves labeling microbes of interest with combinations of probes coupled with spectral imaging to allow the use of fluorophores with highly overlapping excitation and emission spectra. Here, we present the basic principles and theory of CLASI-FISH along with some guidelines for performing CLASI-FISH experiments. We further include a protocol for creating fluorescence spectral reference standards, a vital component of successful CLASI-FISH.

Modeling microbial communities: current, developing, and future technologies for predicting microbial community interaction

Never has there been a greater opportunity for investigating microbial communities. Not only are the profound effects of microbial ecology on every aspect of Earth's geochemical cycles beginning to be understood, but also the analytical and computational tools for investigating microbial Earth are undergoing a rapid revolution. This environmental microbial interactome, the system of interactions between the microbiome and the environment, has shaped the planet's past and will undoubtedly continue to do so in the future. We review recent approaches for modeling microbial community structures and the interactions of microbial populations with their environments. Different modeling approaches consider the environmental microbial interactome from different aspects, and each provides insights to different facets of microbial ecology. We discuss the challenges and opportunities for the future of microbial modeling and describe recent advances in microbial community modeling that are extending current descriptive technologies into a predictive science.

Rapid microbial response to the presence of an ancient relic in the Antarctic Dry Valleys

The extreme cold and aridity of the Antarctic McMurdo Dry Valleys have led to the longstanding belief that metabolic rates of soil microbiota are negligible, and that ecosystem changes take place over millennia. Here we report the first direct experimental evidence that soil microbial communities undergo rapid and lasting changes in response to contemporary environmental conditions. Mummified seals, curious natural artifacts found scattered throughout Dry Valleys, alter their underlying soil environment by stabilizing temperatures, elevating relative humidity and reducing ultraviolet exposure. In a unique, multi-year mummified seal transplantation experiment, we found that endemic Dry Valley microbial communities responded to these changes within 3 years, resulting in a sevenfold increase in CO(2) flux and a significant reduction in biodiversity. These findings challenge prevailing ideas about Antarctic Dry Valley ecosystems and indicate that current and future environmental conditions may strongly influence the ecology of the dominant biota in the Dry Valleys.

Comparative analysis of Korean human gut microbiota by barcoded pyrosequencing

Human gut microbiota plays important roles in harvesting energy from the diet, stimulating the proliferation of the intestinal epithelium, developing the immune system, and regulating fat storage in the host. Characterization of gut microbiota, however, has been limited to western people and is not sufficiently extensive to fully describe microbial communities. In this study, we investigated the overall composition of the gut microbiota and its host specificity and temporal stability in 20 Koreans using 454-pyrosequencing with barcoded primers targeting the V1 to V3 region of the bacterial 16S rRNA gene. A total of 303,402 high quality reads covered each sample and 8,427 reads were analyzed on average. The results were compared with those of individuals from the USA, China and Japan. In general, microbial communities were dominated by five previously identified phyla: Actinobacteria, Firmicutes, Bacteroidetes, Fusobacteria, and Proteobacteria. UPGMA cluster analysis showed that the species composition of gut microbiota was host-specific and stable over the duration of the test period, but the relative abundance of each member fluctuated. 43 core Korean gut microbiota were identified by comparison of sequences from each individual, of which 15 species level phylotypes were related to previously-reported butyrate-producing bacteria. UniFrac analysis revealed that human gut microbiota differed between countries: Korea, USA, Japan and China, but tended to vary less between individual Koreans, suggesting that gut microbial composition is related to internal and external characteristics of each country member such as host genetics and diet styles.

Global patterns in the biogeography of bacterial taxa

Bacteria control major nutrient cycles and directly influence plant, animal and human health. However, we know relatively little about the forces shaping their large-scale ecological ranges. Here, we reveal patterns in the distribution of individual bacterial taxa at multiple levels of phylogenetic resolution within and between Earth's major habitat types. Our analyses suggest that while macro-scale habitats structure bacterial distribution to some degree, abundant bacteria (i.e. detectable using 16S rRNA gene sequencing methods) are confined to single assemblages. Additionally, we show that the most cosmopolitan taxa are also the most abundant in individual assemblages. These results add to the growing body of data that support that the diversity of the overall bacterial metagenome is tremendous. The mechanisms governing microbial distribution remain poorly understood, but our analyses provide a framework with which to test the importance of macro-ecological environmental gradients, relative abundance, neutral processes and the ecological strategies of individual taxa in structuring microbial communities.

Ancient origins determine global biogeography of hot and cold desert cyanobacteria

Factors governing large-scale spatio-temporal distribution of microorganisms remain unresolved, yet are pivotal to understanding ecosystem value and function. Molecular genetic analyses have focused on the influence of niche and neutral processes in determining spatial patterns without considering the temporal scale. Here, we use temporal phylogenetic analysis calibrated using microfossil data for a globally sampled desert cyanobacterium, Chroococcidiopsis, to investigate spatio-temporal patterns in microbial biogeography and evolution. Multilocus phylogenetic associations were dependent on contemporary climate with no evidence for distance-related patterns. Massively parallel pyrosequencing of environmental samples confirmed that Chroococcidiopsis variants were specific to either hot or cold deserts. Temporally scaled phylogenetic analyses showed no evidence of recent inter-regional gene flow, indicating populations have not shared common ancestry since before the formation of modern continents. These results indicate that global distribution of desert cyanobacteria has not resulted from widespread contemporary dispersal but is an ancient evolutionary legacy. This highlights the importance of considering temporal scales in microbial biogeography.

Bacterial community changes during bioremediation of aliphatic hydrocarbon-contaminated soil

The microbial community response during the oxygen biostimulation process of aged oil-polluted soils is poorly documented and there is no reference for the long-term monitoring of the unsaturated zone. To assess the potential effect of air supply on hydrocarbon fate and microbial community structure, two treatments (0 and 0.056 mol h⁻¹ molar flow rate of oxygen) were performed in fixed bed reactors containing oil-polluted soil. Microbial activity was monitored continuously over 2 years throughout the oxygen biostimulation process. Microbial community structure before and after treatment for 12 and 24 months was determined using a dual rRNA/rRNA gene approach, allowing us to characterize bacteria that were presumably metabolically active and therefore responsible for the functionality of the community in this polluted soil. Clone library analysis revealed that the microbial community contained many rare phylotypes. These have never been observed in other studied ecosystems. The bacterial community shifted from Gammaproteobacteria to Actinobacteria during the treatment. Without aeration, the samples were dominated by a phylotype linked to the Streptomyces. Members belonging to eight dominant phylotypes were well adapted to the aeration process. Aeration stimulated an Actinobacteria phylotype that might be involved in restoring the ecosystem studied. Phylogenetic analyses suggested that this phylotype is a novel, deep-branching member of the Actinobacteria related to the well-studied genus Acidimicrobium.

Exploring microbial diversity and taxonomy using SSU rRNA hypervariable tag sequencing

Massively parallel pyrosequencing of hypervariable regions from small subunit ribosomal RNA (SSU rRNA) genes can sample a microbial community two or three orders of magnitude more deeply per dollar and per hour than capillary sequencing of full-length SSU rRNA. As with full-length rRNA surveys, each sequence read is a tag surrogate for a single microbe. However, rather than assigning taxonomy by creating gene trees de novo that include all experimental sequences and certain reference taxa, we compare the hypervariable region tags to an extensive database of rRNA sequences and assign taxonomy based on the best match in a Global Alignment for Sequence Taxonomy (GAST) process. The resulting taxonomic census provides information on both composition and diversity of the microbial community. To determine the effectiveness of using only hypervariable region tags for assessing microbial community membership, we compared the taxonomy assigned to the V3 and V6 hypervariable regions with the taxonomy assigned to full-length SSU rRNA sequences isolated from both the human gut and a deep-sea hydrothermal vent. The hypervariable region tags and full-length rRNA sequences provided equivalent taxonomy and measures of relative abundance of microbial communities, even for tags up to 15% divergent from their nearest reference match. The greater sampling depth per dollar afforded by massively parallel pyrosequencing reveals many more members of the “rare biosphere” than does capillary sequencing of the full-length gene. In addition, tag sequencing eliminates cloning bias and the sequences are short enough to be completely sequenced in a single read, maximizing the number of organisms sampled in a run while minimizing chimera formation. This technique allows the cost-effective exploration of changes in microbial community structure, including the rare biosphere, over space and time and can be applied immediately to initiatives, such as the Human Microbiome Project.

Microbes play a critical role in both human and environmental health. The more we explore microbial populations, the more complexity and diversity we find. Phylogenetic trees based on 16S ribosomal RNA genes have been used with great success to identify microbial taxonomy from DNA alone. New DNA sequencing technologies, such as massively parallel pyrosequencing, can provide orders of magnitude more DNA sequences than ever before, however, the sequences are much shorter, so new methods are necessary to identify the microbes from short DNA tags. We demonstrate the effectiveness of identifying microbial taxa by comparing short tags from 16S hypervariable regions against a large database of known 16S genes. Using this technique, hypervariable region tags provide equivalent taxonomy and relative abundances of microbial communities as full-length rRNA sequences. The greater sampling depth afforded by tag pyrosequencing uncovers not only the dominant microbial species, but many more members of the “rare biosphere” than does capillary sequencing of the full-length gene. Tag pyrosequencing greatly enhances projects exploring composition, diversity, and distribution of microbial populations, such as the Human Microbiome Initiative. A companion paper in PLoS Biology (see Dethlefsen et al., doi:10.1371/journal.pbio.0060280) successfully uses this technique to characterize the effects of antibiotics on the human gut microbiota.

Novelty and uniqueness patterns of rare members of the soil biosphere

Soil bacterial communities typically exhibit a distribution pattern in which most bacterial species are present in low abundance. Due to the relatively small size of most culture-independent sequencing surveys, a detailed phylogenetic analysis of rare members of the community is lacking. To gain access to the rarely sampled soil biosphere, we analyzed a data set of 13,001 near-full-length 16S rRNA gene clones derived from an undisturbed tall grass prairie soil in central Oklahoma. Rare members of the soil bacterial community (empirically defined at two different abundance cutoffs) represented 18.1 to 37.1% of the total number of clones in the data set and were, on average, less similar to their closest relatives in public databases when compared to more abundant members of the community. Detailed phylogenetic analyses indicated that members of the soil rare biosphere either belonged to novel bacterial lineages (members of five novel bacterial phyla identified in the data set, as well as members of multiple novel lineages within previously described phyla or candidate phyla), to lineages that are prevalent in other environments but rarely encountered in soil, or were close relatives to more abundant taxa in the data set. While a fraction of the rare community was closely related to more abundant taxonomic groups in the data set, a significant portion of the rare biosphere represented evolutionarily distinct lineages at various taxonomic cutoffs. We reason that these novelty and uniqueness patterns provide clues regarding the origins and potential ecological roles of members of the soil's rare biosphere.

Scratching the surface of the rare biosphere with ribosomal sequence tag primers

Increasingly large datasets of 16S rRNA gene sequences reveal new information about the extent of microbial diversity and the surprising extent of the rare biosphere. Currently, many of the largest datasets are represented by short and variable ribosomal sequence tags (RSTs) that are limited in their ability to accurately assign sequences to broad-scale phylogenetic trees. In this study, we selected 30 rare RSTs from existing sequence datasets and designed primers to amplify c. 1400 bases of the 16S rRNA gene to determine whether these sequences were represented by existing databases or if they might reveal new lineages within the Bacteria. Approximately one-third of the RST primers successfully amplified longer portions of these low-abundance 16S rRNA genes in a specific manner. Subsequent phylogenetic analysis demonstrated that most of these sequences were (1) distantly related to existing cultivated microorganisms and (2) closely related to uncultivated clone sequences that were recently deposited in GenBank. The presence of so many recently collected 16S rRNA gene reference sequences in existing databases suggests that progress is being made quickly towards a microbial census, one which has begun scratching the surface of the 'rare biosphere'.

Quantitative analysis of three hydrogenotrophic microbial groups, methanogenic archaea, sulfate-reducing bacteria, and acetogenic bacteria, within plaque biofilms associated with human periodontal disease

Human subgingival plaque biofilms are highly complex microbial ecosystems that may depend on H(2)-metabolizing processes. Here we investigated the ubiquity and proportions of methanogenic archaea, sulfate reducers, and acetogens in plaque samples from 102 periodontitis patients. In contrast to the case for 65 healthy control subjects, hydrogenotrophic groups were almost consistently detected in periodontal pockets, with the proportions of methanogens and sulfate reducers being significantly elevated in severe cases. In addition, antagonistic interactions among the three microbial groups indicated that they may function as alternative syntrophic partners of secondary fermenting periodontal pathogens.

Quantitative analysis of three hydrogenotrophic microbial groups, methanogenic archaea, sulfate-reducing bacteria, and acetogenic bacteria, within plaque biofilms associated with human periodontal disease

Human subgingival plaque biofilms are highly complex microbial ecosystems that may depend on H(2)-metabolizing processes. Here we investigated the ubiquity and proportions of methanogenic archaea, sulfate reducers, and acetogens in plaque samples from 102 periodontitis patients. In contrast to the case for 65 healthy control subjects, hydrogenotrophic groups were almost consistently detected in periodontal pockets, with the proportions of methanogens and sulfate reducers being significantly elevated in severe cases. In addition, antagonistic interactions among the three microbial groups indicated that they may function as alternative syntrophic partners of secondary fermenting periodontal pathogens.