VITCOMIC2: visualization tool for the phylogenetic composition of microbial communities based on 16S rRNA gene amplicons and metagenomic shotgun sequencing

Background

The 16S rRNA gene-based amplicon sequencing analysis is widely used to determine the taxonomic composition of microbial communities. Once the taxonomic composition of each community is obtained, evolutionary relationships among taxa are inferred by a phylogenetic tree. Thus, the combined representation of taxonomic composition and phylogenetic relationships among taxa is a powerful method for understanding microbial community structure; however, applying phylogenetic tree-based representation with information on the abundance of thousands or more taxa in each community is a difficult task. For this purpose, we previously developed the tool VITCOMIC (VIsualization tool for Taxonomic COmpositions of MIcrobial Community), which is based on the genome-sequenced microbes’ phylogenetic information. Here, we introduce VITCOMIC2, which incorporates substantive improvements over VITCOMIC that were necessary to address several issues associated with 16S rRNA gene-based analysis of microbial communities.

Results

We developed VITCOMIC2 to provide (i) sequence identity searches against broad reference taxa including uncultured taxa; (ii) normalization of 16S rRNA gene copy number differences among taxa; (iii) rapid sequence identity searches by applying the graphics processing unit-based sequence identity search tool CLAST; (iv) accurate taxonomic composition inference and nearly full-length 16S rRNA gene sequence reconstructions for metagenomic shotgun sequencing; and (v) an interactive user interface for simultaneous representation of the taxonomic composition of microbial communities and phylogenetic relationships among taxa. We validated the accuracy of processes (ii) and (iv) by using metagenomic shotgun sequencing data from a mock microbial community.

Conclusions

The improvements incorporated into VITCOMIC2 enable users to acquire an intuitive understanding of microbial community composition based on the 16S rRNA gene sequence data obtained from both metagenomic shotgun and amplicon sequencing.

     

Incorporating 16S Gene Copy Number Information Improves Estimates of Microbial Diversity and Abundance

The abundance of different SSU rRNA (“16S”) gene sequences in environmental samples is widely used in studies of microbial ecology as a measure of microbial community structure and diversity. However, the genomic copy number of the 16S gene varies greatly – from one in many species to up to 15 in some bacteria and to hundreds in some microbial eukaryotes. As a result of this variation the relative abundance of 16S genes in environmental samples can be attributed both to variation in the relative abundance of different organisms, and to variation in genomic 16S copy number among those organisms. Despite this fact, many studies assume that the abundance of 16S gene sequences is a surrogate measure of the relative abundance of the organisms containing those sequences. Here we present a method that uses data on sequences and genomic copy number of 16S genes along with phylogenetic placement and ancestral state estimation to estimate organismal abundances from environmental DNA sequence data. We use theory and simulations to demonstrate that 16S genomic copy number can be accurately estimated from the short reads typically obtained from high-throughput environmental sequencing of the 16S gene, and that organismal abundances in microbial communities are more strongly correlated with estimated abundances obtained from our method than with gene abundances. We re-analyze several published empirical data sets and demonstrate that the use of gene abundance versus estimated organismal abundance can lead to different inferences about community diversity and structure and the identity of the dominant taxa in microbial communities. Our approach will allow microbial ecologists to make more accurate inferences about microbial diversity and abundance based on 16S sequence data.     

Taxonomic composition of Lake Baikal bacterioneuston communities

The taxonomic composition of microbial communities of Lake Baikal surface microlayer was studied by pyrosequencing of the 16S rDNA amplicons. Statistically reliable differences were found between bacterioneuston of the shallow and deep-water stations. The shallow station community was characterized by higher diversity than the deep-water one. While bacterioneuston communities were shown to be less diverse than the water column communities, their diversity was comparable to that of other biofilm associations. Microbial communities of Lake Baikal surface microlayer were shown to be similar to those of the water column in the composition of predominant phyla, while differing considerably at the genus level. Bacterioneuston of Lake Baikal was comparable to microbial communities of the surface microlayer of other freshwater basins, although it was characterized by high abundance of the Alphaproteobacteria and Verrucomicrobia. High abundance of photoheterotrophs compared to the water column communities of other freshwater basins was another distinctive feature of Lake Baikal bacterioneuston. Our results showed the Lake Baikal surface microlayer to be a specific microbial community with low species diversity and relatively high abundance of photoheterotrophic microorganisms.     

Pyrosequencing of the Chaperonin-60 Universal Target as a Tool for Determining Microbial Community Composition

We compared dideoxy sequencing of cloned chaperonin-60 universal target (cpn60 UT) amplicons to pyrosequencing of amplicons derived from vaginal microbial communities. In samples pooled from a number of individuals, the pyrosequencing method produced a data set that included virtually all of the sequences that were found within the clone library and revealed an additional level of taxonomic richness. However, the relative abundances of the sequences were different in the two datasets. These observations were expanded and confirmed by the analysis of paired clone library and pyrosequencing datasets from vaginal swabs taken from four individuals. Both for individuals with a normal vaginal microbiota and for those with bacterial vaginosis, the pyrosequencing method revealed a large number of low-abundance taxa that were missed by the clone library approach. In addition, we showed that the pyrosequencing method generates a reproducible profile of microbial community structure in replicate amplifications from the same community. We also compared the taxonomic composition of a vaginal microbial community determined by pyrosequencing of 16S rRNA amplicons to that obtained using cpn60 universal primers. We found that the profiles generated by the two molecular targets were highly similar, with slight differences in the proportional representation of the taxa detected. However, the number of operational taxonomic units was significantly higher in the cpn60 data set, suggesting that the protein-encoding gene provides improved species resolution over the 16S rRNA target. These observations demonstrate that pyrosequencing of cpn60 UT amplicons provides a robust, reliable method for deep sequencing of microbial communities.Scientific interest in human microbial communities is growing, and basic concepts about the “human microbiome” are evolving rapidly (3, 34). Molecular phylogenetic analysis of 16S rRNA-encoding DNA sequences has revealed a vast diversity of uncultured microbial symbionts that influence animal physiology in ways only beginning to be understood. In particular, microbial species inhabiting the human vagina are thought to play an important role in host health (10). A shift in the composition of the vaginal microbiota from “normal” (Lactobacillus dominated) to a state defined as bacterial vaginosis (BV; increased abundance of gram-negative organisms) is associated with a range of negative outcomes, including pelvic inflammatory disease, preterm births, and the acquisition of sexually transmitted diseases (21, 22, 37). This observation has led to an increased interest in determining the composition of the vaginal microbiota by culture-independent methods (8, 11, 17, 25, 30, 35, 36). However, established cloning and sequencing techniques remain time- and labor-intensive, severely limiting the reach of phylogenetic or functional surveys of microbial communities across body sites, individuals, geographic areas, and scales of time.The advent of next-generation ultra-high-throughput sequencing technologies, in particular, the GS FLX (454 Life Sciences, Branford, CT), has removed an important quantitative barrier in molecular analysis by increasing the number of reads from a gene or genome by orders of magnitude in a single run (20). Unfortunately, the short average length of pyrosequencing reads (∼200 bp compared to ∼700 bp using dideoxy sequencing) presents a new set of problems. The results of recent application of this technology to analysis of 16S rRNA gene sequences from microbes in vaginal samples have demonstrated that short reads are more likely to generate matches to multiple sequences in the rRNA sequence database and that taxonomic and phylogenetic resolution was limited due to strong similarities between 16S rRNA sequences from closely related species (32).An alternative molecular target for microbial identification and phylogenetic analysis is cpn60, a gene that encodes the 60-kDa chaperonin or heat shock protein (HSP60/GroEL) (13). The cpn60 gene is universal in eubacteria and eukaryotes and an extensive, curated reference database is available (13) (http://cpndb.cbr.nrc.ca). The cpn60 universal target (UT) offers key advantages, including short target length (549 to 567 bp), sufficient resolving power to distinguish closely related species and subspecies, and a relatively uniform distribution of variability across the entire length of the target (9, 12). The use of the cpn60 UT has been well established for phylogenetic analysis of complex samples (4, 14) and has recently been applied to vaginal microbial communities (11). In the present study, we examined the feasibility of pyrosequencing for determining the composition of the vaginal microbiota using the cpn60 UT. We compared the microbial community structure generated by pyrosequencing of cpn60 amplicons using the GS FLX with dideoxy sequencing based on clone libraries generated from the same samples. In addition, we evaluated the microbial community profiles generated by pyrosequencing of cpn60 UT amplicons and 16S rRNA amplicons from the same vaginal samples.     

Short pyrosequencing reads suffice for accurate microbial community analysis

Pyrosequencing technology allows us to characterize microbial communities using 16S ribosomal RNA (rRNA) sequences orders of magnitude faster and more cheaply than has previously been possible. However, results from different studies using pyrosequencing and traditional sequencing are often difficult to compare, because amplicons covering different regions of the rRNA might yield different conclusions. We used sequences from over 200 globally dispersed environments to test whether studies that used similar primers clustered together mistakenly, without regard to environment. We then tested whether primer choice affects sequence-based community analyses using UniFrac, our recently-developed method for comparing microbial communities. We performed three tests of primer effects. We tested whether different simulated amplicons generated the same UniFrac clustering results as near-full-length sequences for three recent large-scale studies of microbial communities in the mouse and human gut, and the Guerrero Negro microbial mat. We then repeated this analysis for short sequences (100-, 150-, 200- and 250-base reads) resembling those produced by pyrosequencing. The results show that sequencing effort is best focused on gathering more short sequences rather than fewer longer ones, provided that the primers are chosen wisely, and that community comparison methods such as UniFrac are surprisingly robust to variation in the region sequenced.     

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.     

Design and Experimental Application of a Novel Non-Degenerate Universal Primer Set that Amplifies Prokaryotic 16S rRNA Genes with a Low Possibility to Amplify Eukaryotic rRNA Genes

The deep sequencing of 16S rRNA genes amplified by universal primers has revolutionized our understanding of microbial communities by allowing the characterization of the diversity of the uncultured majority. However, some universal primers also amplify eukaryotic rRNA genes, leading to a decrease in the efficiency of sequencing of prokaryotic 16S rRNA genes with possible mischaracterization of the diversity in the microbial community. In this study, we compared 16S rRNA gene sequences from genome-sequenced strains and identified candidates for non-degenerate universal primers that could be used for the amplification of prokaryotic 16S rRNA genes. The 50 identified candidates were investigated to calculate their coverage for prokaryotic and eukaryotic rRNA genes, including those from uncultured taxa and eukaryotic organelles, and a novel universal primer set, 342F-806R, covering many prokaryotic, but not eukaryotic, rRNA genes was identified. This primer set was validated by the amplification of 16S rRNA genes from a soil metagenomic sample and subsequent pyrosequencing using the Roche 454 platform. The same sample was also used for pyrosequencing of the amplicons by employing a commonly used primer set, 338F-533R, and for shotgun metagenomic sequencing using the Illumina platform. Our comparison of the taxonomic compositions inferred by the three sequencing experiments indicated that the non-degenerate 342F-806R primer set can characterize the taxonomic composition of the microbial community without substantial bias, and is highly expected to be applicable to the analysis of a wide variety of microbial communities.     

Impact of biochar application to soil on the root-associated bacterial community structure of fully developed greenhouse pepper plants

Adding biochar to soil has environmental and agricultural potential due to its long-term carbon sequestration capacity and its ability to improve crop productivity. Recent studies have demonstrated that soil-applied biochar promotes the systemic resistance of plants to several prominent foliar pathogens. One potential mechanism for this phenomenon is root-associated microbial elicitors whose presence is somehow augmented in the biochar-amended soils. The objective of this study was to assess the effect of biochar amendment on the root-associated bacterial community composition of mature sweet pepper (Capsicum annuum L.) plants. Molecular fingerprinting (denaturing gradient gel electrophoresis and terminal restriction fragment length polymorphism) of 16S rRNA gene fragments showed a clear differentiation between the root-associated bacterial community structures of biochar-amended and control plants. The pyrosequencing of 16S rRNA amplicons from the rhizoplane of both treatments generated a total of 20,142 sequences, 92 to 95% of which were affiliated with the Proteobacteria, Bacteroidetes, Actinobacteria, and Firmicutes phyla. The relative abundance of members of the Bacteroidetes phylum increased from 12 to 30% as a result of biochar amendment, while that of the Proteobacteria decreased from 71 to 47%. The Bacteroidetes-affiliated Flavobacterium was the strongest biochar-induced genus. The relative abundance of this group increased from 4.2% of total root-associated operational taxonomic units (OTUs) in control samples to 19.6% in biochar-amended samples. Additional biochar-induced genera included chitin and cellulose degraders (Chitinophaga and Cellvibrio, respectively) and aromatic compound degraders (Hydrogenophaga and Dechloromonas). We hypothesize that these biochar-augmented genera may be at least partially responsible for the beneficial effect of biochar amendment on plant growth and viability.     

Impact of Substratum Surface on Microbial Community Structure and Treatment Performance in Biological Aerated Filters

The impact of substratum surface property change on biofilm community structure was investigated using laboratory biological aerated filter (BAF) reactors and molecular microbial community analysis. Two substratum surfaces that differed in surface properties were created via surface coating and used to develop biofilms in test (modified surface) and control (original surface) BAF reactors. Microbial community analysis by 16S rRNA gene-based PCR-denaturing gradient gel electrophoresis (DGGE) showed that the surface property change consistently resulted in distinct profiles of microbial populations during replicate reactor start-ups. Pyrosequencing of the bar-coded 16S rRNA gene amplicons surveyed more than 90% of the microbial diversity in the microbial communities and identified 72 unique bacterial species within 19 bacterial orders. Among the 19 orders of bacteria detected, Burkholderiales and Rhodocyclales of the Betaproteobacteria class were numerically dominant and accounted for 90.5 to 97.4% of the sequence reads, and their relative abundances in the test and control BAF reactors were different in consistent patterns during the two reactor start-ups. Three of the five dominant bacterial species also showed consistent relative abundance changes between the test and control BAF reactors. The different biofilm microbial communities led to different treatment efficiencies, with consistently higher total organic carbon (TOC) removal in the test reactor than in the control reactor. Further understanding of how surface properties affect biofilm microbial communities and functional performance would enable the rational design of new generations of substrata for the improvement of biofilm-based biological treatment processes.     

Rare but active taxa contribute to community dynamics of benthic biofilms in glacier‐fed streams

Glaciers harbour diverse microorganisms, which upon ice melt can be released downstream. In glacier‐fed streams microorganisms can attach to stones or sediments to form benthic biofilms. We used 454‐pyrosequencing to explore the bulk (16S rDNA) and putatively active (16S rRNA) microbial communities of stone and sediment biofilms across 26 glacier‐fed streams. We found differences in community composition between bulk and active communities among streams and a stronger congruence between biofilm types. Relative abundances of rRNA and rDNA were positively correlated across different taxa and taxonomic levels, but at lower taxonomic levels, the higher abundance in either the active or the bulk communities became more apparent. Here, environmental variables played a minor role in structuring active communities. However, we found a large number of rare taxa with higher relative abundances in rRNA compared with rDNA. This suggests that rare taxa contribute disproportionately to microbial community dynamics in glacier‐fed streams. Our findings propose that high community turnover, where taxa repeatedly enter and leave the ‘seed bank’, contributes to the maintenance of microbial biodiversity in harsh ecosystems with continuous environmental perturbations, such as glacier‐fed streams.     

Multiple copies of 16s rRNA gene affect the restriction patterns and DGGE profile as revealed by analysis of genome database

The use of 16S rRNA gene has been a "golden" method to determine the diversity of microbial communities in environmental samples, phylogenetic relationships of prokaryotes and taxonomic position of newly isolated organisms. However due to the presence of multiple heterogeneous 16S rRNA gene copies in many strains, the interpretation of microbial ecology via 16S rRNA sequences is complicated. Purpose of present paper is to demonstrate the extent to which the multiple heterogeneous 16S rRNA gene copies affect RFLP patterns and DGG E profiles by using the genome database. In present genome database, there are 782 bacterial strains in total whose genomes have been completely sequenced and annotated. Among the total strains, 639 strains (82%) possess multiple 16S rRNA gene copies, 415 strains (53%) whose multiple copies are heterogeneous in sequences as revealed by alignment, 236 strains (30%) whose multiple copies show different restrict patterns by CSP61 + Hinfl, MspI + Rsal or HhaI as analyzed in silico. Polymorphisms of the multiple copies in certain strains were further characterized by G + C% and phylogentic distances based on the sequences of V3 region, which are linked to DGGE patters. Polymorphisms of a few strains were shown as examples. Using artificial communities, it is demonstrated that the presence of multiple heterogeneous 16S rRNA gene copies potentially leads to over-estimation of the diversity of a community. It is suggested that care must be taken when interpreting 16S rRNA-based RFLP and DGGE data and profiling an environmental community.     

Phylogenetic Profiling and Diversity of Bacterial Communities in the Death Valley,an Extreme Habitat in the Atacama Desert

The Atacama Desert, one of the driest deserts in the world, represents a unique extreme environmental ecosystem to explore the bacterial diversity as it is considered to be at the dry limit for life. A 16S rRNA gene (spanning the hyper variable V3 region) library was constructed from an alkaline sample of unvegetated soil at the hyperarid margin in the Atacama Desert. A total of 244 clone sequences were used for MOTHUR analysis, which revealed 20 unique phylotypes or operational taxonomic units (OTUs). V3 region amplicons of the 16S rRNA were suitable for distinguishing the bacterial community to the genus and specie level. We found that all OTUs were affiliated with taxa representative of the Firmicutes phylum. The extremely high abundance of Firmicutes indicated that most bacteria in the soil were spore-forming survivors. In this study we detected a narrower diversity as compared to other ecological studies performed in other areas of the Atacama Desert. The reported genera were Oceanobacillus (representing the 69.5 % of the clones sequenced), Bacillus, Thalassobacillus and Virgibacillus. The present work shows physical and chemical parameters have a prominent impact on the microbial community structure. It constitutes an example of the communities adapted to live in extreme conditions caused by dryness and metal concentrations .

Electronic supplementary material

The online version of this article (doi:10.1007/s12088-015-0539-3) contains supplementary material, which is available to authorized users.     

Microbial rRNA gene expression and co‐occurrence profiles associate with biokinetics and elemental composition in full‐scale anaerobic digesters

This study examined whether the abundance and expression of microbial 16S rRNA genes were associated with elemental concentrations and substrate conversion biokinetics in 20 full‐scale anaerobic digesters, including seven municipal sewage sludge (SS) digesters and 13 industrial codigesters. SS digester contents had higher methane production rates from acetate, propionate and phenyl acetate compared to industrial codigesters. SS digesters and industrial codigesters were distinctly clustered based on their elemental concentrations, with higher concentrations of NH₃‐N, Cl, K and Na observed in codigesters. Amplicon sequencing of 16S rRNA genes and reverse‐transcribed 16S rRNA revealed divergent grouping of microbial communities between mesophilic SS digesters, mesophilic codigesters and thermophilic digesters. Higher intradigester distances between Archaea 16S rRNA and rRNA gene profiles were observed in mesophilic codigesters, which also had the lowest acetate utilization biokinetics. Constrained ordination showed that microbial rRNA and rRNA gene profiles were significantly associated with maximum methane production rates from acetate, propionate, oleate and phenyl acetate, as well as concentrations of NH₃‐N, Fe, S, Mo and Ni. A co‐occurrence network of rRNA gene expression confirmed the three main clusters of anaerobic digester communities based on active populations. Syntrophic and methanogenic taxa were highly represented within the subnetworks, indicating that obligate energy‐sharing partnerships play critical roles in stabilizing the digester microbiome. Overall, these results provide new evidence showing that different feed substrates associate with different micronutrient compositions in anaerobic digesters, which in turn may influence microbial abundance, activity and function.     

Microbial Communities Can Be Described by Metabolic Structure: A General Framework and Application to a Seasonally Variable,Depth-Stratified Microbial Community from the Coastal West Antarctic Peninsula

Taxonomic marker gene studies, such as the 16S rRNA gene, have been used to successfully explore microbial diversity in a variety of marine, terrestrial, and host environments. For some of these environments long term sampling programs are beginning to build a historical record of microbial community structure. Although these 16S rRNA gene datasets do not intrinsically provide information on microbial metabolism or ecosystem function, this information can be developed by identifying metabolisms associated with related, phenotyped strains. Here we introduce the concept of metabolic inference; the systematic prediction of metabolism from phylogeny, and describe a complete pipeline for predicting the metabolic pathways likely to be found in a collection of 16S rRNA gene phylotypes. This framework includes a mechanism for assigning confidence to each metabolic inference that is based on a novel method for evaluating genomic plasticity. We applied this framework to 16S rRNA gene libraries from the West Antarctic Peninsula marine environment, including surface and deep summer samples and surface winter samples. Using statistical methods commonly applied to community ecology data we found that metabolic structure differed between summer surface and winter and deep samples, comparable to an analysis of community structure by 16S rRNA gene phylotypes. While taxonomic variance between samples was primarily driven by low abundance taxa, metabolic variance was attributable to both high and low abundance pathways. This suggests that clades with a high degree of functional redundancy can occupy distinct adjacent niches. Overall our findings demonstrate that inferred metabolism can be used in place of taxonomy to describe the structure of microbial communities. Coupling metabolic inference with targeted metagenomics and an improved collection of completed genomes could be a powerful way to analyze microbial communities in a high-throughput manner that provides direct access to metabolic and ecosystem function.     

Different response of bacteria,archaea and fungi to process parameters in nine full-scale anaerobic digesters

Biogas production is a biotechnological process realized by complex bacterial, archaeal and likely fungal communities. Their composition was assessed in nine full-scale biogas plants with distinctly differing feedstock input and process parameters. This study investigated the actually active microbial community members by using a comprehensive sequencing approach based on ribosomal 16S and 28S rRNA fragments. The prevailing taxonomical units of each respective community were subsequently linked to process parameters. Ribosomal rRNA of bacteria, archaea and fungi, respectively, showed different compositions with respect to process parameters and supplied feedstocks: (i) bacterial communities were affected by the key factors temperature and ammonium concentration; (ii) composition of archaea was mainly related to process temperature; and (iii) relative abundance of fungi was linked to feedstocks supplied to the digesters. Anaerobic digesters with a high methane yield showed remarkably similar bacterial communities regarding identified taxonomic families. Although archaeal communities differed strongly on genus level from each other, the respective digesters still showed high methane yields. Functional redundancy of the archaeal communities may explain this effect. 28S rRNA sequences of fungi in all nine full-scale anaerobic digesters were primarily classified as facultative anaerobic Ascomycota and Basidiomycota. Since the presence of ribosomal 28S rRNA indicates that fungi may be active in the biogas digesters, further research should be carried out to examine to which extent they are important players in anaerobic digestion processes.     

Environmental factors influencing the distribution of rRNA from Verrucomicrobia in soil

The Verrucomicrobia constitute a newly discovered division of the Bacteria identified as a numerically abundant component of soil microbial communities in numerous sites around the world. The relative abundance of rRNA from Verrucomicrobia was investigated in the soil to examine the influence of specific environmental factors on the distribution of Verrucomicrobia and to better understand the distribution of this group in terrestrial ecosystems. The abundance of the verrucomicrobial rRNA was determined by using a novel oligonucleotide probe that is specific for verrucomicrobial 16S rRNA. The abundance of verrucomicrobial 16S rRNA in soil microbial communities was determined in relation to plant community composition and soil management history over a period of 2 years. Additional samples were analyzed to determine if verrucomicrobial rRNA relative abundance changes in relation to either soil depth or soil moisture content. The Verrucomicrobia composed 1.9+/-0.2% of the microbial community rRNA present in the 85 soil samples examined. The distribution of verrucomicrobial rRNA in the soil reveals that Verrucomicrobia are significantly affected by environmental characteristics that change in relation to time, soil history, and soil depth, and reveals that a statistically significant amount of the variation in verrucomicrobial rRNA abundance can be explained by changes in soil moisture content.     

Diversity of Kenyan soda lake alkaliphiles assessed by molecular methods

DNA was extracted from water and sediment samples taken from soda lakes of the Kenyan-Tanzanian Rift Valley. DNA was also extracted from microbial enrichment cultures of sediment samples. 16S rRNA genes were amplified by the polymerase chain reaction and microbial diversity was studied using denaturing gradient gel electrophoresis (DGGE) of 16S rDNA amplicons. Cloning and sequencing of single DGGE bands showed that they usually contained mixed amplicons. Several of the amplicon sequences had high identities, up to 99%, with 16S rRNA genes of organisms previously isolated from soda lakes, while others were only distantly related, with identities as low as 82%. Phylogenetic analysis of the sequenced amplicons indicated that sequences were related to the haloarchaeal, Bacillus/Clostridium, Rhodobacterium/Thioalcalovibrio/ Methylobacter, and Cytophaga/Flavobacterium/Bacteroides (CFB) groups and the enterobacteria/Aeromonas/Vibrio part of the 3 subdivision of the Proteobacteria.Communicated by K. Horikoshi     

MetaGeniE: Characterizing Human Clinical Samples Using Deep Metagenomic Sequencing

With the decreasing cost of next-generation sequencing, deep sequencing of clinical samples provides unique opportunities to understand host-associated microbial communities. Among the primary challenges of clinical metagenomic sequencing is the rapid filtering of human reads to survey for pathogens with high specificity and sensitivity. Metagenomes are inherently variable due to different microbes in the samples and their relative abundance, the size and architecture of genomes, and factors such as target DNA amounts in tissue samples (i.e. human DNA versus pathogen DNA concentration). This variation in metagenomes typically manifests in sequencing datasets as low pathogen abundance, a high number of host reads, and the presence of close relatives and complex microbial communities. In addition to these challenges posed by the composition of metagenomes, high numbers of reads generated from high-throughput deep sequencing pose immense computational challenges. Accurate identification of pathogens is confounded by individual reads mapping to multiple different reference genomes due to gene similarity in different taxa present in the community or close relatives in the reference database. Available global and local sequence aligners also vary in sensitivity, specificity, and speed of detection. The efficiency of detection of pathogens in clinical samples is largely dependent on the desired taxonomic resolution of the organisms. We have developed an efficient strategy that identifies “all against all” relationships between sequencing reads and reference genomes. Our approach allows for scaling to large reference databases and then genome reconstruction by aggregating global and local alignments, thus allowing genetic characterization of pathogens at higher taxonomic resolution. These results were consistent with strain level SNP genotyping and bacterial identification from laboratory culture.     

Nutrient and acetate amendment leads to acetoclastic methane production and microbial community change in a non‐producing Australian coal well

Coal mining is responsible for 11% of total anthropogenic methane emission thereby contributing considerably to climate change. Attempts to harvest coalbed methane for energy production are challenged by relatively low methane concentrations. In this study, we investigated whether nutrient and acetate amendment of a non‐producing sub‐bituminous coal well could transform the system to a methane source. We tracked cell counts, methane production, acetate concentration and geochemical parameters for 25 months in one amended and one unamended coal well in Australia. Additionally, the microbial community was analysed with 16S rRNA gene amplicon sequencing at 17 and 25 months after amendment and complemented by metagenome sequencing at 25 months. We found that cell numbers increased rapidly from 3.0 × 10⁴ cells ml^?1 to 9.9 × 10⁷ in the first 7 months after amendment. However, acetate depletion with concomitant methane production started only after 12–19 months. The microbial community was dominated by complex organic compound degraders (Anaerolineaceae, Rhodocyclaceae and Geobacter spp.), acetoclastic methanogens (Methanothrix spp.) and fungi (Agaricomycetes). Even though the microbial community had the functional potential to convert coal to methane, we observed no indication that coal was actually converted within the time frame of the study. Our results suggest that even though nutrient and acetate amendment stimulated relevant microbial species, it is not a sustainable way to transform non‐producing coal wells into bioenergy factories.     

Microbial community composition in sediments resists perturbation by nutrient enrichment

Functional redundancy in bacterial communities is expected to allow microbial assemblages to survive perturbation by allowing continuity in function despite compositional changes in communities. Recent evidence suggests, however, that microbial communities change both composition and function as a result of disturbance. We present evidence for a third response: resistance. We examined microbial community response to perturbation caused by nutrient enrichment in salt marsh sediments using deep pyrosequencing of 16S rRNA and functional gene microarrays targeting the nirS gene. Composition of the microbial community, as demonstrated by both genes, was unaffected by significant variations in external nutrient supply in our sampling locations, despite demonstrable and diverse nutrient-induced changes in many aspects of marsh ecology. The lack of response to external forcing demonstrates a remarkable uncoupling between microbial composition and ecosystem-level biogeochemical processes and suggests that sediment microbial communities are able to resist some forms of perturbation.