Phylogenetic specificity and reproducibility and new method for analysis of terminal restriction fragment profiles of 16S rRNA genes from bacterial communities

Terminal restriction fragment (TRF) analysis of 16S rRNA genes is an increasingly popular method for rapid comparison of microbial communities, but analysis of the data is still in a developmental stage. We assessed the phylogenetic resolution and reproducibility of TRF profiles in order to evaluate the limitations of the method, and we developed an essential analysis technique to improve the interpretation of TRF data. The theoretical phylogenetic resolution of TRF profiles was determined based on the specificity of TRFs predicted from 3,908 16S rRNA gene sequences. With sequences from the Proteobacteria or gram-positive division, as much as 73% of the TRFs were phylogenetically specific (representing strains from at most two genera). However, the fraction decreased when sequences from the two divisions were combined. The data show that phylogenetic inference will be most effective if TRF profiles represent only a single bacterial division or smaller group. The analytical precision of the TRF method was assessed by comparing nine replicate profiles of a single soil DNA sample. Despite meticulous care in producing the replicates, numerous small, irreproducible peaks were observed. As many as 85% of the 169 distinct TRFs found among the profiles were irreproducible (i.e., not present in all nine replicates). Substantial variation also occurred in the height of synonymous peaks. To make comparisons of microbial communities more reliable, we developed an analytical procedure that reduces variation and extracts a reproducible subset of data from replicate TRF profiles. The procedure can also be used with other DNA fingerprinting techniques for microbial communities or microbial genomes.     

T-Align, a web-based tool for comparison of multiple terminal restriction fragment length polymorphism profiles

Terminal restriction fragment length polymorphism (tRFLP) is a potentially high-throughput method for the analysis of complex microbial communities. Comparison of multiple tRFLP profiles to identify shared and unique components of microbial communities however, is done manually, which is both time consuming and error prone. This paper describes a freely accessible web-based program, T-Align (http://inismor.ucd.ie/~talign/), which addresses this problem. Initially replicate profiles are compared and used to generate a single consensus profile containing only terminal restriction fragments that occur in all replicate profiles. Subsequently consensus profiles representing different communities are compared to produce a list showing whether a terminal restriction fragment (TRF) is present in a particular sample and its relative fluorescence intensity. The use of T-Align thus allows rapid comparison of numerous tRFLP profiles. T-Align is demonstrated by alignment of tRFLP profiles generated from bacterioplankton communities collected from the Irish and Celtic Seas in November 2000. Ubiquitous TRFs and site-specific TRFs were identified using T-Align.     

Reproducibility and quantitation of amplicon sequencing-based detection

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

Interpreting Ecological Diversity Indices Applied to Terminal Restriction Fragment Length Polymorphism Data: Insights from Simulated Microbial Communities

Ecological diversity indices are frequently applied to molecular profiling methods, such as terminal restriction fragment length polymorphism (T-RFLP), in order to compare diversity among microbial communities. We performed simulations to determine whether diversity indices calculated from T-RFLP profiles could reflect the true diversity of the underlying communities despite potential analytical artifacts. These include multiple taxa generating the same terminal restriction fragment (TRF) and rare TRFs being excluded by a relative abundance (fluorescence) threshold. True community diversity was simulated using the lognormal species abundance distribution. Simulated T-RFLP profiles were generated by assigning each species a TRF size based on an empirical or modeled TRF size distribution. With a typical threshold (1%), the only consistently useful relationship was between Smith and Wilson evenness applied to T-RFLP data (TRF-E_var) and true Shannon diversity (H′), with correlations between 0.71 and 0.81. TRF-H′ and true H′ were well correlated in the simulations using the lowest number of species, but this correlation declined substantially in simulations using greater numbers of species, to the point where TRF-H′ cannot be considered a useful statistic. The relationships between TRF diversity indices and true indices were sensitive to the relative abundance threshold, with greatly improved correlations observed using a 0.1% threshold, which was investigated for comparative purposes but is not possible to consistently achieve with current technology. In general, the use of diversity indices on T-RFLP data provides inaccurate estimates of true diversity in microbial communities (with the possible exception of TRF-E_var). We suggest that, where significant differences in T-RFLP diversity indices were found in previous work, these should be reinterpreted as a reflection of differences in community composition rather than a true difference in community diversity.     

Phylogenetic diversity, composition and distribution of bacterioplankton community in the Dongjiang River, China

Bacterioplankton community compositions in the Dongjiang River were characterized using denaturing gradient gel electrophoresis (DGGE) and 16S rRNA gene clone library construction. Water samples in nine different sites were taken along the mainstem and three tributaries. In total, 24 bands from DGGE gels and 406 clones from the libraries were selected and sequenced, subsequently analyzed for the bacterial diversity and composition of those microbial communities. Bacterial 16S rRNA gene sequences from freshwater bacteria exhibited board phylogenetic diversity, including sequences representing the Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria, Actinobacteria, Bacteriodetes, Verrucomicrobia, and candidate division TM7. Members of Betaproteobacteria group were the most dominant in all sampling sites, followed by Gammaproteobacteria, Alphaproteobacteria, and Actinobacteria. DGGE profiles and the ∫-LIBSHUFF analysis revealed similar patterns of bacterial diversity among most sampling sites, while spatial distribution variances existed in all sites along the river basin. Statistical analysis showed that bacterial species distribution strongly correlated with environmental variables, such as nitrate and ammonia, suggesting that nitrogen nutrients may shape the microbial community structure and composition in the Dongjiang River. This study had important implications for the comparison with other rivers elsewhere and contributed to the growing data set on the factors that structure bacterial communities in freshwater ecosystems.     

Assessment of microbial diversity in four southwestern United States soils by 16S rRNA gene terminal restriction fragment analysis

The ability of terminal restriction fragment (T-RFLP or TRF) profiles of 16S rRNA genes to provide useful information about the relative diversity of complex microbial communities was investigated by comparison with other methods. Four soil communities representing two pinyon rhizosphere and two between-tree (interspace) soil environments were compared by analysis of 16S rRNA gene clone libraries and culture collections (Dunbar et al., Appl. Environ. Microbiol. 65:1662-1669, 1998) and by analysis of 16S rDNA TRF profiles of community DNA. The TRF method was able to differentiate the four communities in a manner consistent with previous comparisons of the communities by analysis of 16S rDNA clone libraries. TRF profiles were not useful for calculating and comparing traditional community richness or evenness values among the four soil environments. Statistics calculated from RsaI, HhaI, HaeIII, and MspI profiles of each community were inconsistent, and the combined data were not significantly different between samples. The detection sensitivity of the method was tested. In standard PCRs, a seeded population comprising 0.1 to 1% of the total community could be detected. The combined results demonstrate that TRF analysis is an excellent method for rapidly comparing the relationships between bacterial communities in environmental samples. However, for highly complex communities, the method appears unable to provide classical measures of relative community diversity.     

Sample pooling obscures diversity patterns in intertidal ciliate community composition and structure

The aim of this study was to assess the effect of sample pooling on the portrayal of ciliate community structure and composition in intertidal sediment samples. Molecular ciliate community profiles were obtained from nine biological replicates distributed in three discrete sampling plots and from samples that were pooled prior to RNA extraction using terminal restriction fragment polymorphism (T-RFLP) analyses of SSU rRNA. Comparing the individual replicates of one sampling plot with each other, we found a differential variability among the individual biological replicates. T-RFLP profiles of pooled samples displayed a significantly different community composition compared with the cumulative individual biological replicate samples. We conclude that sample pooling obscures diversity patterns in ciliate and possibly also other microbial eukaryote studies. However, differences between pooled samples and replicates were less pronounced when community structure was analyzed. We found that the most abundant T-RFLP peaks were generally shared between biological replicates and pooled samples. Assuming that the most abundant taxa in an ecosystem under study are also the ones driving ecosystem processes, sample pooling may still be effective for the analyses of ecological key players.     

Interpreting ecological diversity indices applied to terminal restriction fragment length polymorphism data: insights from simulated microbial communities

Ecological diversity indices are frequently applied to molecular profiling methods, such as terminal restriction fragment length polymorphism (T-RFLP), in order to compare diversity among microbial communities. We performed simulations to determine whether diversity indices calculated from T-RFLP profiles could reflect the true diversity of the underlying communities despite potential analytical artifacts. These include multiple taxa generating the same terminal restriction fragment (TRF) and rare TRFs being excluded by a relative abundance (fluorescence) threshold. True community diversity was simulated using the lognormal species abundance distribution. Simulated T-RFLP profiles were generated by assigning each species a TRF size based on an empirical or modeled TRF size distribution. With a typical threshold (1%), the only consistently useful relationship was between Smith and Wilson evenness applied to T-RFLP data (TRF-E(var)) and true Shannon diversity (H'), with correlations between 0.71 and 0.81. TRF-H' and true H' were well correlated in the simulations using the lowest number of species, but this correlation declined substantially in simulations using greater numbers of species, to the point where TRF-H' cannot be considered a useful statistic. The relationships between TRF diversity indices and true indices were sensitive to the relative abundance threshold, with greatly improved correlations observed using a 0.1% threshold, which was investigated for comparative purposes but is not possible to consistently achieve with current technology. In general, the use of diversity indices on T-RFLP data provides inaccurate estimates of true diversity in microbial communities (with the possible exception of TRF-E(var)). We suggest that, where significant differences in T-RFLP diversity indices were found in previous work, these should be reinterpreted as a reflection of differences in community composition rather than a true difference in community diversity.     

Web-based phylogenetic assignment tool for analysis of terminal restriction fragment length polymorphism profiles of microbial communities

Culture-independent DNA fingerprints are commonly used to assess the diversity of a microbial community. However, relating species composition to community profiles produced by community fingerprint methods is not straightforward. Terminal restriction fragment length polymorphism (T-RFLP) is a community fingerprint method in which phylogenetic assignments may be inferred from the terminal restriction fragment (T-RF) sizes through the use of web-based resources that predict T-RF sizes for known bacteria. The process quickly becomes computationally intensive due to the need to analyze profiles produced by multiple restriction digests and the complexity of profiles generated by natural microbial communities. A web-based tool is described here that rapidly generates phylogenetic assignments from submitted community T-RFLP profiles based on a database of fragments produced by known 16S rRNA gene sequences. Users have the option of submitting a customized database generated from unpublished sequences or from a gene other than the 16S rRNA gene. This phylogenetic assignment tool allows users to employ T-RFLP to simultaneously analyze microbial community diversity and species composition. An analysis of the variability of bacterial species composition throughout the water column in a humic lake was carried out to demonstrate the functionality of the phylogenetic assignment tool. This method was validated by comparing the results generated by this program with results from a 16S rRNA gene clone library.     

Web-Based Phylogenetic Assignment Tool for Analysis of Terminal Restriction Fragment Length Polymorphism Profiles of Microbial Communities

Culture-independent DNA fingerprints are commonly used to assess the diversity of a microbial community. However, relating species composition to community profiles produced by community fingerprint methods is not straightforward. Terminal restriction fragment length polymorphism (T-RFLP) is a community fingerprint method in which phylogenetic assignments may be inferred from the terminal restriction fragment (T-RF) sizes through the use of web-based resources that predict T-RF sizes for known bacteria. The process quickly becomes computationally intensive due to the need to analyze profiles produced by multiple restriction digests and the complexity of profiles generated by natural microbial communities. A web-based tool is described here that rapidly generates phylogenetic assignments from submitted community T-RFLP profiles based on a database of fragments produced by known 16S rRNA gene sequences. Users have the option of submitting a customized database generated from unpublished sequences or from a gene other than the 16S rRNA gene. This phylogenetic assignment tool allows users to employ T-RFLP to simultaneously analyze microbial community diversity and species composition. An analysis of the variability of bacterial species composition throughout the water column in a humic lake was carried out to demonstrate the functionality of the phylogenetic assignment tool. This method was validated by comparing the results generated by this program with results from a 16S rRNA gene clone library.     

Short-Read Assembly of Full-Length 16S Amplicons Reveals Bacterial Diversity in Subsurface Sediments

In microbial ecology, a fundamental question relates to how community diversity and composition change in response to perturbation. Most studies have had limited ability to deeply sample community structure (e.g. Sanger-sequenced 16S rRNA libraries), or have had limited taxonomic resolution (e.g. studies based on 16S rRNA hypervariable region sequencing). Here, we combine the higher taxonomic resolution of near-full-length 16S rRNA gene amplicons with the economics and sensitivity of short-read sequencing to assay the abundance and identity of organisms that represent as little as 0.01% of sediment bacterial communities. We used a new version of EMIRGE optimized for large data size to reconstruct near-full-length 16S rRNA genes from amplicons sheared and sequenced with Illumina technology. The approach allowed us to differentiate the community composition among samples acquired before perturbation, after acetate amendment shifted the predominant metabolism to iron reduction, and once sulfate reduction began. Results were highly reproducible across technical replicates, and identified specific taxa that responded to the perturbation. All samples contain very high alpha diversity and abundant organisms from phyla without cultivated representatives. Surprisingly, at the time points measured, there was no strong loss of evenness, despite the selective pressure of acetate amendment and change in the terminal electron accepting process. However, community membership was altered significantly. The method allows for sensitive, accurate profiling of the “long tail” of low abundance organisms that exist in many microbial communities, and can resolve population dynamics in response to environmental change.     

Anaerobic microbial communities in Lake Pavin, a unique meromictic lake in France

The Bacteria and Archaea from the meromictic Lake Pavin were analyzed in samples collected along a vertical profile in the anoxic monimolimnion and were compared to those in samples from the oxic mixolimnion. Nine targeted 16S rRNA oligonucleotide probes were used to assess the distribution of Bacteria and Archaea and to investigate the in situ occurrence of sulfate-reducing bacteria and methane-producing Archaea involved in the terminal steps of the anaerobic degradation of organic material. The diversity of the complex microbial communities was assessed from the 16S rRNA polymorphisms present in terminal restriction fragment (TRF) depth patterns. The densities of the microbial community increased in the anoxic layer, and Archaea detected with probe ARCH915 represented the largest microbial group in the water column, with a mean Archaea/Eubacteria ratio of 1.5. Terminal restriction fragment length polymorphism (TRFLP) analysis revealed an elevated archaeal and bacterial phylotype richness in anoxic bottom-water samples. The structure of the Archaea community remained rather homogeneous, while TRFLP patterns for the eubacterial community revealed a heterogeneous distribution of eubacterial TRFs.     

Effect of hydrogen on soil bacterial community structure in two soils as determined by terminal restriction fragment length polymorphism

The metabolism of hydrogen evolved from HUP^? legume nodules can alter bacterial community structures in the rhizosphere. Our earlier experiments demonstrated increased hydrogen uptake and appearance of white spots within bacterial colonies in H₂-treated soil. We were also able to isolate hydrogen-oxidizing bacteria from soil samples exposed to hydrogen, but not from samples exposed to air. To further understand the effect of hydrogen metabolism on soil microbial communities, in this study 16S rRNA terminal restriction fragment (TRF) profiles of different soil samples exposed to hydrogen gas under laboratory, greenhouse, and field conditions were analyzed. Relationships between soil bacterial community structures from hydrogen-treated soil samples and controls, illustrated by UPGMA (unpaired group mathematical averages) dendrograms, indicated a significant contribution of hydrogen metabolism to the variation in bacterial community. The intensity variation of TRF peaks includes both hydrogen-utilizing bacteria, whose growth were stimulated by hydrogen exposure, and other bacterial species whose growth was inhibited. Comparison of TRF profiles between laboratory and greenhouse samples showed that T-RFLP is a useful technique in the detection of root-related effects on soil bacterial community structure.     

A novel method for determining microflora composition using dynamic phylogenetic analysis of 16S ribosomal RNA deep sequencing data

Deep sequencing of the 16S rRNA gene provides a comprehensive view of bacterial communities in a particular environment and has expanded our ability to study the impact of the microflora on human health and disease. Current analysis methods rely on comparisons of the sequences generated with an expanding but limited set of annotated 16S rRNA sequences or phylogenic clustering of sequences based on arbitrary similarity cutoffs. We describe a novel approach to characterize bacterial composition using deep sequencing of 16S rRNA gene. Our method defines operational taxonomic units based on phylogenetic tree reconstruction and dynamic clustering of sequences using solely sequencing data. These OTUs can be used to identify differences in bacteria abundance between environments. This approach can perform better than previous phylogenetic methods and will significantly improve our understanding of the microfloral role on human diseases by providing a comprehensive analysis of the microbial composition from various bacterial communities.     

Biogeography of Actinomycete Communities and Type II Polyketide Synthase Genes in Soils Collected in New Jersey and Central Asia

Soil microbial communities are believed to be comprised of thousands of different bacterial species. One prevailing idea is that “everything is everywhere, and the environment selects,” implying that all types of bacteria are present in all environments where their growth requirements are met. We tested this hypothesis using actinomycete communities and type II polyketide synthase (PKS) genes found in soils collected from New Jersey and Uzbekistan (n = 91). Terminal restriction fragment length polymorphism analysis using actinomycete 16S rRNA and type II PKS genes was employed to determine community profiles. The terminal fragment frequencies in soil samples had a lognormal distribution, indicating that the majority of actinomycete phylotypes and PKS pathways are present infrequently in the environment. Less than 1% of peaks were detected in more than 50% of samples, and as many as 18% of the fragments were unique and detected in only one sample. Actinomycete 16S rRNA fingerprints clustered by country of origin, indicating that unique populations are present in North America and Central Asia. Sequence analysis of type II PKS gene fragments cloned from Uzbek soil revealed 35 novel sequence clades whose levels of identity to genes in the GenBank database ranged from 68 to 92%. The data indicate that actinomycetes are patchily distributed but that distinct populations are present in North American and Central Asia. These results have implications for microbial bioprospecting and indicate that the cosmopolitan actinomycete species and PKS pathways may account for only a small proportion of the total diversity in soil.     

Identification and phylogeny of Arabian snakes: Comparison of venom chromatographic profiles versus 16S rRNA gene sequences

Identification of snake species is important for various reasons including the emergency treatment of snake bite victims. We present a simple method for identification of six snake species using the gel filtration chromatographic profiles of their venoms. The venoms of Echis coloratus, Echis pyramidum, Cerastes gasperettii, Bitis arietans, Naja arabica, and Walterinnesia aegyptia were milked, lyophilized, diluted and centrifuged to separate the mucus from the venom. The clear supernatants were filtered and chromatographed on fast protein liquid chromatography (FPLC). We obtained the 16S rRNA gene sequences of the above species and performed phylogenetic analysis using the neighbor-joining method. The chromatograms of venoms from different snake species showed peculiar patterns based on the number and location of peaks. The dendrograms generated from similarity matrix based on the presence/absence of particular chromatographic peaks clearly differentiated Elapids from Viperids. Molecular cladistics using 16S rRNA gene sequences resulted in jumping clades while separating the members of these two families. These findings suggest that chromatographic profiles of snake venoms may provide a simple and reproducible chemical fingerprinting method for quick identification of snake species. However, the validation of this methodology requires further studies on large number of specimens from within and across species.     

The phylogenetic diversity of metagenomes

Phylogenetic diversity--patterns of phylogenetic relatedness among organisms in ecological communities--provides important insights into the mechanisms underlying community assembly. Studies that measure phylogenetic diversity in microbial communities have primarily been limited to a single marker gene approach, using the small subunit of the rRNA gene (SSU-rRNA) to quantify phylogenetic relationships among microbial taxa. In this study, we present an approach for inferring phylogenetic relationships among microorganisms based on the random metagenomic sequencing of DNA fragments. To overcome challenges caused by the fragmentary nature of metagenomic data, we leveraged fully sequenced bacterial genomes as a scaffold to enable inference of phylogenetic relationships among metagenomic sequences from multiple phylogenetic marker gene families. The resulting metagenomic phylogeny can be used to quantify the phylogenetic diversity of microbial communities based on metagenomic data sets. We applied this method to understand patterns of microbial phylogenetic diversity and community assembly along an oceanic depth gradient, and compared our findings to previous studies of this gradient using SSU-rRNA gene and metagenomic analyses. Bacterial phylogenetic diversity was highest at intermediate depths beneath the ocean surface, whereas taxonomic diversity (diversity measured by binning sequences into taxonomically similar groups) showed no relationship with depth. Phylogenetic diversity estimates based on the SSU-rRNA gene and the multi-gene metagenomic phylogeny were broadly concordant, suggesting that our approach will be applicable to other metagenomic data sets for which corresponding SSU-rRNA gene sequences are unavailable. Our approach opens up the possibility of using metagenomic data to study microbial diversity in a phylogenetic context.     

Chitinase genes revealed and compared in bacterial isolates, DNA extracts and a metagenomic library from a phytopathogen-suppressive soil

Soil that is suppressive to disease caused by fungal pathogens is an interesting source to target for novel chitinases that might be contributing towards disease suppression. In this study, we screened for chitinase genes, in a phytopathogen-suppressive soil in three ways: (1) from a metagenomic library constructed from microbial cells extracted from soil, (2) from directly extracted DNA and (3) from bacterial isolates with antifungal and chitinase activities. Terminal restriction fragment length polymorphism (T-RFLP) of chitinase genes revealed differences in amplified chitinase genes from the metagenomic library and the directly extracted DNA, but approximately 40% of the identified chitinase terminal restriction fragments (TRFs) were found in both sources. All of the chitinase TRFs from the isolates were matched to TRFs in the directly extracted DNA and the metagenomic library. The most abundant chitinase TRF in the soil DNA and the metagenomic library corresponded to the TRF¹⁰³ of the isolate Streptomyces mutomycini and/or Streptomyces clavifer . There were good matches between T-RFLP profiles of chitinase gene fragments obtained from different sources of DNA. However, there were also differences in both the chitinase and the 16S rRNA gene T-RFLP patterns depending on the source of DNA, emphasizing the lack of complete coverage of the gene diversity by any of the approaches used.     

Metagenomics: Read Length Matters

Obtaining an unbiased view of the phylogenetic composition and functional diversity within a microbial community is one central objective of metagenomic analysis. New technologies, such as 454 pyrosequencing, have dramatically reduced sequencing costs, to a level where metagenomic analysis may become a viable alternative to more-focused assessments of the phylogenetic (e.g., 16S rRNA genes) and functional diversity of microbial communities. To determine whether the short (~100 to 200 bp) sequence reads obtained from pyrosequencing are appropriate for the phylogenetic and functional characterization of microbial communities, the results of BLAST and COG analyses were compared for long (~750 bp) and randomly derived short reads from each of two microbial and one virioplankton metagenome libraries. Overall, BLASTX searches against the GenBank nr database found far fewer homologs within the short-sequence libraries. This was especially pronounced for a Chesapeake Bay virioplankton metagenome library. Increasing the short-read sampling depth or the length of derived short reads (up to 400 bp) did not completely resolve the discrepancy in BLASTX homolog detection. Only in cases where the long-read sequence had a close homolog (low BLAST E-score) did the derived short-read sequence also find a significant homolog. Thus, more-distant homologs of microbial and viral genes are not detected by short-read sequences. Among COG hits, derived short reads sampled at a depth of two short reads per long read missed up to 72% of the COG hits found using long reads. Noting the current limitation in computational approaches for the analysis of short sequences, the use of short-read-length libraries does not appear to be an appropriate tool for the metagenomic characterization of microbial communities.