Improved serial analysis of V1 ribosomal sequence tags (SARST-V1) provides a rapid, comprehensive, sequence-based characterization of bacterial diversity and community composition

Serial analysis of ribosomal sequence tags (SARST) is a recently developed technology that can generate large 16S rRNA gene (rrs) sequence data sets from microbiomes, but there are numerous enzymatic and purification steps required to construct the ribosomal sequence tag (RST) clone libraries. We report here an improved SARST method, which still targets the V1 hypervariable region of rrs genes, but reduces the number of enzymes, oligonucleotides, reagents, and technical steps needed to produce the RST clone libraries. The new method, hereafter referred to as SARST-V1, was used to examine the eubacterial diversity present in community DNA recovered from the microbiome resident in the ovine rumen. The 190 sequenced clones contained 1055 RSTs and no less than 236 unique phylotypes (based on > or = 95% sequence identity) that were assigned to eight different eubacterial phyla. Rarefaction and monomolecular curve analyses predicted that the complete RST clone library contains 99% of the 353 unique phylotypes predicted to exist in this microbiome. When compared with ribosomal intergenic spacer analysis (RISA) of the same community DNA sample, as well as a compilation of nine previously published conventional rrs clone libraries prepared from the same type of samples, the RST clone library provided a more comprehensive characterization of the eubacterial diversity present in rumen microbiomes. As such, SARST-V1 should be a useful tool applicable to comprehensive examination of diversity and composition in microbiomes and offers an affordable, sequence-based method for diversity analysis.     

Unexpectedly High Bacterial Diversity in Arctic Tundra Relative to Boreal Forest Soils, Revealed by Serial Analysis of Ribosomal Sequence Tags

Arctic tundra and boreal forest soils have globally relevant functions that affect atmospheric chemistry and climate, yet the bacterial composition and diversity of these soils have received little study. Serial analysis of ribosomal sequence tags (SARST) and denaturing gradient gel electrophoresis (DGGE) were used to compare composite soil samples taken from boreal and arctic biomes. This study comprises an extensive comparison of geographically distant soil bacterial communities, involving the analysis of 12,850 ribosomal sequence tags from six composite soil samples. Bacterial diversity estimates were greater for undisturbed arctic tundra soil samples than for boreal forest soil samples, with the highest diversity associated with a sample from an extreme northern location (82^oN). The lowest diversity estimate was obtained from an arctic soil sample that was disturbed by compaction and sampled from a greater depth. Since samples from the two biomes did not form distinct clusters on the basis of SARST data and DGGE fingerprints, factors other than latitude likely influenced the phylogenetic compositions of these communities. The high number of ribosomal sequences analyzed enabled the identification of possible cosmopolitan and endemic bacterial distributions in particular soils.     

Unexpectedly high bacterial diversity in arctic tundra relative to boreal forest soils, revealed by serial analysis of ribosomal sequence tags

Arctic tundra and boreal forest soils have globally relevant functions that affect atmospheric chemistry and climate, yet the bacterial composition and diversity of these soils have received little study. Serial analysis of ribosomal sequence tags (SARST) and denaturing gradient gel electrophoresis (DGGE) were used to compare composite soil samples taken from boreal and arctic biomes. This study comprises an extensive comparison of geographically distant soil bacterial communities, involving the analysis of 12,850 ribosomal sequence tags from six composite soil samples. Bacterial diversity estimates were greater for undisturbed arctic tundra soil samples than for boreal forest soil samples, with the highest diversity associated with a sample from an extreme northern location (82(o)N). The lowest diversity estimate was obtained from an arctic soil sample that was disturbed by compaction and sampled from a greater depth. Since samples from the two biomes did not form distinct clusters on the basis of SARST data and DGGE fingerprints, factors other than latitude likely influenced the phylogenetic compositions of these communities. The high number of ribosomal sequences analyzed enabled the identification of possible cosmopolitan and endemic bacterial distributions in particular soils.     

Composition of microbial communities in hexachlorocyclohexane (HCH) contaminated soils from Spain revealed with a habitat-specific microarray

Microarray technology was used to characterize and compare hexachlorocyclohexane (HCH) contaminated soils from Spain. A library of 2,290 hypervariable 16S rRNA gene sequences was prepared with serial analysis of ribosomal sequence tags (SARST) from a composite of contaminated and uncontaminated soils. By designing hybridization probes specific to the 100 most abundant ribosomal sequence tags (RSTs) in the composite library, the RST array was designed to be habitat-specific and predicted to monitor the most abundant polymerase chain reaction (PCR)-amplified phylotypes in the individual samples. The sensitivity and specificity of the RST array was tested with a series of pure culture-specific probes and hybridized with labelled soil PCR products to generate hybridization patterns for each soil. Sequencing of prominent bands in denaturing gradient gel electrophoresis (DGGE) fingerprints derived from these soils provided a means by which we successfully confirmed the habitat-specific array design and validated the bulk of the probe signals. Non-metric multidimensional scaling revealed correlations between probe signals and soil physicochemical parameters. Among the strongest correlations to total HCH contamination were probe signals corresponding to unknown Gamma Proteobacteria, potential pollutant-degrading phylotypes, and several organisms with acid-tolerant phenotypes. The strongest correlations to alpha-HCH were probe signals corresponding to the genus Sphingomonas, which contains known HCH degraders. This suggests that the population detected was enriched in situ by HCH contamination and may play a role in HCH degradation. Other environmental parameters were also likely instrumental in shaping community composition in these soils. The results highlight the power of habitat-specific microarrays for comparing complex microbial communities.     

Serial analysis of V6 ribosomal sequence tags (SARST-V6): a method for efficient, high-throughput analysis of microbial community composition

Serial analysis of ribosomal sequence tags (SARST) is a novel technique for characterizing microbial community composition. The SARST method captures sequence information from concatemers of short 16S rDNA polymerase chain reaction (PCR) amplicons from complex populations of DNA. Here, we describe a similar method, serial analysis of V6 ribosomal sequence tags (SARST-V6), which targets the V6 hypervariable region of bacterial 16S rRNA genes. The SARST-V6 technique exploits internal primer sequences to generate compatible restriction digest overhangs, thereby improving upon the efficiency of SARST. Serial analysis of V6 ribosomal sequence tags of bacterial community composition in hydrothermal marine sediments from Guaymas Basin resembled results of cloning and sequencing of single, full-length PCR products from ribosomal RNA genes of the same microbial community. Both methods identified the same major bacterial groups, but only SARST-V6 recovered thermodesulfobacteria and gamma-proteobacteria sequences, while only full-length PCR product cloning recovered candidate division OP11 se-quences. There were differences in the relative frequencies of some phylotypes. The disparities reflect differences in the amplicon pool obtained during initial amplification that may result from different primer affinities or DNA degradation. These results demonstrate the utility of SARST-V6 in collecting taxonomically informative data for high-throughput analysis of microbial communities.     

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'.     

Phylogenetic analysis of ribosomal RNA operons from uncultivated coastal marine bacterioplankton

Analyses of small subunit ribosomal RNA genes (SSU rDNAs) have significantly influenced our understanding of the composition of aquatic microbial assemblages. Unfortunately, SSU rDNA sequences often do not have sufficient resolving power to differentiate closely related species. To address this general problem for uncultivated bacterioplankton taxa, we analysed and compared sequences of polymerase chain reaction (PCR)-generated and bacterial artificial chromosome (BAC)-derived clones that contained most of the SSU rDNAs, the internal transcribed spacer (ITS) and the large subunit ribosomal RNA gene (LSU rDNA). The phylogenetic representation in the rRNA operon PCR library was similar to that reported previously in coastal bacterioplankton SSU rDNA libraries. We observed good concordance between the phylogenetic relationships among coastal bacterioplankton inferred from SSU or LSU rDNA sequences. ITS sequences confirmed the close intragroup relationships among members of the SAR11, SAR116 and SAR86 clades that were predicted by SSU and LSU rDNA sequence analyses. We also found strong support for homologous recombination between the ITS regions of operons from the SAR11 clade.     

Genetic Diversity of Microbial Eukaryotes in Anoxic Sediment of the Saline Meromictic Lake Namako-ike (Japan): On the Detection of Anaerobic or Anoxic-tolerant Lineages of Eukaryotes

Available sequence data on eukaryotic small-subunit ribosomal DNA (SSU rDNA) directly retrieved from various environments have increased recently, and the diversity of microbial eukaryotes (protists) has been shown to be much greater than previously expected. However, the molecular information accumulated to date does still not thoroughly reveal ecological distribution patterns of microbial eukaryotes. In the ongoing challenge to detect anaerobic or anoxic-tolerant lineages of eukaryotes, we directly extracted DNA from the anoxic sediment of a saline meromictic lake, constructed genetic libraries of PCR-amplified SSU rDNA, and performed phylogenetic analyses with the cloned SSU rDNA sequences. Although a few sequences could not be confidently assigned to any major eukaryotic groups in the analyses and are debatable regarding their taxonomic positions, most sequences obtained have affiliations with known major lineages of eukaryotes (Cercozoa, Alveolata, Stramenopiles, and Opisthokonta). Among these sequences, some branched with lineages predominantly composed of uncultured environmental clones retrieved from other anoxic environments, while others were closely related to those of eukaryotic parasites (e.g. Phytomyxea of Cercozoa, Gregarinea of Alveolata, and Ichthyosporea of Opisthokonta).     

Microbial composition of near-boiling silica-depositing thermal springs throughout Yellowstone National Park

The extent of hyperthermophilic microbial diversity associated with siliceous sinter (geyserite) was characterized in seven near-boiling silica-depositing springs throughout Yellowstone National Park using environmental PCR amplification of small-subunit rRNA genes (SSU rDNA), large-subunit rDNA, and the internal transcribed spacer (ITS). We found that Thermocrinis ruber, a member of the order Aquificales, is ubiquitous, an indication that primary production in these springs is driven by hydrogen oxidation. Several other lineages with no known close relatives were identified that branch among the hyperthermophilic bacteria. Although they all branch deep in the bacterial tree, the precise phylogenetic placement of many of these lineages is unresolved at this time. While some springs contained a fair amount of phylogenetic diversity, others did not. Within the same spring, communities in the subaqueous environment were not appreciably different than those in the splash zone at the edge of the pool, although a greater number of phylotypes was found along the pool's edge. Also, microbial community composition appeared to have little correlation with the type of sinter morphology. The number of cell morphotypes identified by fluorescence in situ hybridization and scanning electron microscopy was greater than the number of phylotypes in SSU clone libraries. Despite little variation in Thermocrinis ruber SSU sequences, abundant variation was found in the hypervariable ITS region. The distribution of ITS sequence types appeared to be correlated with distinct morphotypes of Thermocrinis ruber in different pools. Therefore, species- or subspecies-level divergences are present but not detectable in highly conserved SSU sequences.     

Diversity of picoplanktonic prasinophytes assessed by direct nuclear SSU rDNA sequencing of environmental samples and novel isolates retrieved from oceanic and coastal marine ecosystems

Picoplanktonic prasinophytes are well represented in culture collections and marine samples. In order to better characterize this ecologically important group, we compared the phylogenetic diversity of picoplanktonic prasinophyte strains available at the Roscoff Culture Collection (RCC) and that of nuclear SSU rDNA sequences from environmental clone libraries obtained from oceanic and coastal ecosystems. Among the 570 strains avalaible, 91 belonged to prasinophytes, 65 were partially sequenced, and we obtained the entire SSU rDNA sequence for a selection of 14 strains. Within the 18 available environmental clone libraries, the prasinophytes accounted for 12% of the total number of clones retrieved (142 partial sequences in total), and we selected 9 clones to obtain entire SSU rDNA sequence. Using this approach, we obtained a subsequent genetic database that revealed the presence of seven independent lineages among prasinophytes, including a novel clade (clade VII). This new clade groups the genus Picocystis, two unidentified coccoid strains, and 4 environmental sequences. For each of these seven lineages, at least one representative is available in culture. The three picoplanktonic genera Ostreococcus, Micromonas, and Bathycoccus (order Mamiellales), were the best represented prasinophytes both in cultures and genetic libraries. SSU rDNA phylogenetic analyses suggest that the genus Bathycoccus forms a very homogeneous group. In contrast, the genera Micromonas and Ostreococcus turned out to be quite complex, consisting of three and four independent lineages, respectively. This report of the overall diversity of picoeukaryotic prasinophytes reveals a group of ecologically important and diverse marine microorganims that are well represented by isolated cultures.     

Comparative microbial diversity in the gastrointestinal tracts of food animal species

Molecular tools based on small subunit (SSU) rDNA gene sequencesoffer a powerful and rapid tool for the analysis of complexmicrobial communities found in the gastrointestinal tracts (GIT)of food animal species. Extensive comparative sequence analysisof SSU rRNA molecules representing a wide diversity of organismsshows that different regions of the molecule vary in sequenceconservation. Oligonucleotides complementing regions of universallyconversed SSU rRNA sequences are used as universal probes, whilethose complementing more variable regions of sequence are usefulas selective probes targeting species, genus, or phylogeneticgroups. Different approaches derive different information andthis is highly dependent on the type of target nucleic acidemployed and the conceptual and technical basis used for nucleicacid probe design. Generally these approaches can be dividedinto DNA-based methods employing empirically characterized probesand rRNA-based methods based on comparative sequence analysisfor design and interpretation of "rational" probes. Polymerasechain reaction (PCR) based techniques can also be applied tothe analysis of microbial communities in the GIT. Direct cloningof SSU rDNA genes amplified from these complex communities canbe used to determine the extent of diversity in these GIT communities.Denaturing gradient gel electrophoresis (DGGE) is another powerfultool for profiling microbial diversity of microbial communitiesin GI tracts. Sequence analysis of the excised DGGE ampliconscan then be used to presumptively identify predominant bacterialspecies. Examples of how these molecular approaches are beingused to study the microbial diversity of communities from steersfed different diets, swine fed probiotics, and Atlantic salmonfed aquaculture diets are presented.     

Microbial Composition of Near-Boiling Silica-Depositing Thermal Springs throughout Yellowstone National Park

The extent of hyperthermophilic microbial diversity associated with siliceous sinter (geyserite) was characterized in seven near-boiling silica-depositing springs throughout Yellowstone National Park using environmental PCR amplification of small-subunit rRNA genes (SSU rDNA), large-subunit rDNA, and the internal transcribed spacer (ITS). We found that Thermocrinis ruber, a member of the order Aquificales, is ubiquitous, an indication that primary production in these springs is driven by hydrogen oxidation. Several other lineages with no known close relatives were identified that branch among the hyperthermophilic bacteria. Although they all branch deep in the bacterial tree, the precise phylogenetic placement of many of these lineages is unresolved at this time. While some springs contained a fair amount of phylogenetic diversity, others did not. Within the same spring, communities in the subaqueous environment were not appreciably different than those in the splash zone at the edge of the pool, although a greater number of phylotypes was found along the pool's edge. Also, microbial community composition appeared to have little correlation with the type of sinter morphology. The number of cell morphotypes identified by fluorescence in situ hybridization and scanning electron microscopy was greater than the number of phylotypes in SSU clone libraries. Despite little variation in Thermocrinis ruber SSU sequences, abundant variation was found in the hypervariable ITS region. The distribution of ITS sequence types appeared to be correlated with distinct morphotypes of Thermocrinis ruber in different pools. Therefore, species- or subspecies-level divergences are present but not detectable in highly conserved SSU sequences.     

Molecular analysis of bacterial community structure and diversity in unimproved and improved upland grass pastures

Bacterial community structure and diversity in rhizospheres in two types of grassland, distinguished by both plant species and fertilization regimen, were assessed by performing a 16S ribosomal DNA (rDNA) sequence analysis of DNAs extracted from triplicate soil plots. PCR products were cloned, and 45 to 48 clones from each of the six libraries were partially sequenced. Phylogenetic analysis of the resultant 275 clone sequences indicated that there was considerable variation in abundance in replicate unfertilized, unimproved soil samples and fertilized, improved soil samples but that there were no significant differences in the abundance of any phylogenetic group. Several clone sequences were identical in the 16S rDNA region analyzed, and the clones comprised eight pairs of duplicate clones and two sets of triplicate clones. Many clones were found to be most closely related to environmental clones obtained in other studies, although three clones were found to be identical to culturable species in databases. The clones were clustered into operational taxonomic units at a level of sequence similarity of >97% in order to quantify diversity. In all, 34 clusters containing two or more sequences were identified, and the largest group contained nine clones. A number of diversity, dominance, and evenness indices were calculated, and they all indicated that diversity was high, reflecting the low coverage of rDNA libraries achieved. Differences in diversity between sample types were not observed. Collector's curves, however, indicated that there were differences in the underlying community structures; in particular, there was reduced diversity of organisms of the alpha subdivision of the class Proteobacteria (alpha-proteobacteria) in improved soils.     

Effect of PCR amplicon size on assessments of clone library microbial diversity and community structure

PCR-based surveys of microbial communities commonly use regions of the small-subunit ribosomal RNA (SSU rRNA) gene to determine taxonomic membership and estimate total diversity. Here we show that the length of the target amplicon has a significant effect on assessments of microbial richness and community membership. Using operational taxonomic unit (OTU)- and taxonomy-based tools, we compared the V6 hypervariable region of the bacterial SSU rRNA gene of three amplicon libraries of c. 100, 400 and 1000 base pairs (bp) from each of two hydrothermal vent fluid samples. We found that the smallest amplicon libraries contained more unique sequences, higher diversity estimates and a different community structure than the other two libraries from each sample. We hypothesize that a combination of polymerase dissociation, cloning bias and mispriming due to secondary structure accounts for the differences. While this relationship is not linear, it is clear that the smallest amplicon libraries contained more different types of sequences, and accordingly, more diverse members of the community. Because divergent and lower abundant taxa can be more readily detected with smaller amplicons, they may provide better assessments of total community diversity and taxonomic membership than longer amplicons in molecular studies of microbial communities.     

Soil bacterial community shift correlated with change from forest to pasture vegetation in a tropical soil.

The change in vegetative cover of a Hawaiian soil from forest to pasture led to significant changes in the composition of the soil bacterial community. DNAs were extracted from both soil habitats and compared for the abundance of guanine-plus-cytosine (G+C) content, by analysis of abundance of phylotypes of small-subunit ribosomal DNA (SSU rDNA) amplified from fractions with 63 and 35% G+C contents, and by phylogenetic analysis of the dominant rDNA clones in the 63% G+C content fraction. All three methods showed differences between the forest and pasture habitats, providing evidence that vegetation had a strong influence on microbial community composition at three levels of taxon resolution. The forest soil DNA had a peak in G+C content of 61%, while the DNA of the pasture soil had a peak in G+C content of 67%. None of the dominant phylotypes found in the forest soil were detected in the pasture soil. For the 63% G+C fraction SSU rDNA sequence analysis of the three most dominant members revealed that their phyla changed from Fibrobacter and Syntrophomonas assemblages in the forest soil to Burkholderia and Rhizobium-Agrobacterium assemblages in the pasture soil.     

An advanced molecular strategy to identify bacterial communities on art objects

The application of culture-independent techniques based on molecular biological methods, especially on the PCR amplification of 16S rRNA genes, attempts to overcome some shortcomings of conventional cultivation methods and reveals far more complex bacterial communities on art objects than can be shown by cultivation methods. One of the major challenges of investigating microbial growth on art objects by molecular means is the extraction of DNA, due to small sample amounts and PCR inhibitors. In the present study, we introduce a DNA extraction protocol, which allowed the extraction of PCR-amplifiable DNA from samples derived from lime wall paintings and loamy soil underground. The DNA extracts were used to amplify 16S ribosomal fragments, which were subsequently analyzed by denaturing gradient gel electrophoresis (DGGE). In parallel with the DGGE analysis, clone libraries containing PCR fragments of the ribosomal gene were constructed and clones were screened by DGGE. Clone libraries allow the inclusion of the entire 16S rDNA sequence in the phylogenetic analyses of microorganisms, providing a more reliable phylogenetic identification of microorganisms than is obtained from sequence analyses of excised and directly sequenced DGGE bands.     

Illumina-based analysis of microbial community diversity

Microbes commonly exist in milieus of varying complexity and diversity. Although cultivation-based techniques have been unable to accurately capture the true diversity within microbial communities, these deficiencies have been overcome by applying molecular approaches that target the universally conserved 16S ribosomal RNA gene. The recent application of 454 pyrosequencing to simultaneously sequence thousands of 16S rDNA sequences (pyrotags) has revolutionized the characterization of complex microbial communities. To date, studies based on 454 pyrotags have dominated the field, but sequencing platforms that generate many more sequence reads at much lower costs have been developed. Here, we use the Illumina sequencing platform to design a strategy for 16S amplicon analysis (iTags), and assess its generality, practicality and potential complications. We fabricated and sequenced paired-end libraries of amplified hyper-variable 16S rDNA fragments from sets of samples that varied in their contents, ranging from a single bacterium to highly complex communities. We adopted an approach that allowed us to evaluate several potential sources of errors, including sequencing artifacts, amplification biases, non-corresponding paired-end reads and mistakes in taxonomic classification. By considering each source of error, we delineate ways to make biologically relevant and robust conclusions from the millions of sequencing reads that can be readily generated by this technology.