Eukaryotic richness in the abyss: insights from pyrotag sequencing

Background

The deep sea floor is considered one of the most diverse ecosystems on Earth. Recent environmental DNA surveys based on clone libraries of rRNA genes confirm this observation and reveal a high diversity of eukaryotes present in deep-sea sediment samples. However, environmental clone-library surveys yield only a modest number of sequences with which to evaluate the diversity of abyssal eukaryotes.

Methodology/Principal Findings

Here, we examined the richness of eukaryotic DNA in deep Arctic and Southern Ocean samples using massively parallel sequencing of the 18S ribosomal RNA (rRNA) V9 hypervariable region. In very small volumes of sediments, ranging from 0.35 to 0.7 g, we recovered up to 7,499 unique sequences per sample. By clustering sequences having up to 3 differences, we observed from 942 to 1756 Operational Taxonomic Units (OTUs) per sample. Taxonomic analyses of these OTUs showed that DNA of all major groups of eukaryotes is represented at the deep-sea floor. The dinoflagellates, cercozoans, ciliates, and euglenozoans predominate, contributing to 17%, 16%, 10%, and 8% of all assigned OTUs, respectively. Interestingly, many sequences represent photosynthetic taxa or are similar to those reported from the environmental surveys of surface waters. Moreover, each sample contained from 31 to 71 different metazoan OTUs despite the small sample volume collected. This indicates that a significant faction of the eukaryotic DNA sequences likely do not belong to living organisms, but represent either free, extracellular DNA or remains and resting stages of planktonic species.

Conclusions/Significance

In view of our study, the deep-sea floor appears as a global DNA repository, which preserves genetic information about organisms living in the sediment, as well as in the water column above it. This information can be used for future monitoring of past and present environmental changes.     

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.     

Comparative Analysis of Pyrosequencing and a Phylogenetic Microarray for Exploring Microbial Community Structures in the Human Distal Intestine

Background

Variations in the composition of the human intestinal microbiota are linked to diverse health conditions. High-throughput molecular technologies have recently elucidated microbial community structure at much higher resolution than was previously possible. Here we compare two such methods, pyrosequencing and a phylogenetic array, and evaluate classifications based on two variable 16S rRNA gene regions.

Methods and Findings

Over 1.75 million amplicon sequences were generated from the V4 and V6 regions of 16S rRNA genes in bacterial DNA extracted from four fecal samples of elderly individuals. The phylotype richness, for individual samples, was 1,400–1,800 for V4 reads and 12,500 for V6 reads, and 5,200 unique phylotypes when combining V4 reads from all samples. The RDP-classifier was more efficient for the V4 than for the far less conserved and shorter V6 region, but differences in community structure also affected efficiency. Even when analyzing only 20% of the reads, the majority of the microbial diversity was captured in two samples tested. DNA from the four samples was hybridized against the Human Intestinal Tract (HIT) Chip, a phylogenetic microarray for community profiling. Comparison of clustering of genus counts from pyrosequencing and HITChip data revealed highly similar profiles. Furthermore, correlations of sequence abundance and hybridization signal intensities were very high for lower-order ranks, but lower at family-level, which was probably due to ambiguous taxonomic groupings.

Conclusions

The RDP-classifier consistently assigned most V4 sequences from human intestinal samples down to genus-level with good accuracy and speed. This is the deepest sequencing of single gastrointestinal samples reported to date, but microbial richness levels have still not leveled out. A majority of these diversities can also be captured with five times lower sampling-depth. HITChip hybridizations and resulting community profiles correlate well with pyrosequencing-based compositions, especially for lower-order ranks, indicating high robustness of both approaches. However, incompatible grouping schemes make exact comparison difficult.     

Effect of Cry1Ab protein on rhizobacterial communities of Bt-maize over a four-year cultivation period

Background

Bt-maize is a transgenic variety of maize expressing the Cry toxin from Bacillus turingiensis. The potential accumulation of the relative effect of the transgenic modification and the cry toxin on the rhizobacterial communities of Bt-maize has been monitored over a period of four years.

Methodology/Principal Findings

The accumulative effects of the cultivation of this transgenic plant have been monitored by means of high throughput DNA pyrosequencing of the bacterial DNA coding for the 16S rRNA hypervariable V6 region from rhizobacterial communities. The obtained sequences were subjected to taxonomic, phylogenetic and taxonomic-independent diversity studies. The results obtained were consistent, indicating that variations detected in the rhizobacterial community structure were possibly due to climatic factors rather than to the presence of the Bt-gene. No variations were observed in the diversity estimates between non-Bt and Bt-maize.

Conclusions/Significance

The cultivation of Bt-maize during the four-year period did not change the maize rhizobacterial communities when compared to those of the non-Bt maize. This is the first study to be conducted with Bt-maize during such a long cultivation period and the first evaluation of rhizobacterial communities to be performed in this transgenic plant using Next Generation Sequencing.     

Development of an Analysis Pipeline Characterizing Multiple Hypervariable Regions of 16S rRNA Using Mock Samples

Objectives

There is much speculation on which hypervariable region provides the highest bacterial specificity in 16S rRNA sequencing. The optimum solution to prevent bias and to obtain a comprehensive view of complex bacterial communities would be to sequence the entire 16S rRNA gene; however, this is not possible with second generation standard library design and short-read next-generation sequencing technology.

Methods

This paper examines a new process using seven hypervariable or V regions of the 16S rRNA (six amplicons: V2, V3, V4, V6-7, V8, and V9) processed simultaneously on the Ion Torrent Personal Genome Machine (Life Technologies, Grand Island, NY). Four mock samples were amplified using the 16S Ion Metagenomics Kit^™ (Life Technologies) and their sequencing data is subjected to a novel analytical pipeline.

Results

Results are presented at family and genus level. The Kullback-Leibler divergence (D_KL), a measure of the departure of the computed from the nominal bacterial distribution in the mock samples, was used to infer which region performed best at the family and genus levels. Three different hypervariable regions, V2, V4, and V6-7, produced the lowest divergence compared to the known mock sample. The V9 region gave the highest (worst) average D_KL while the V4 gave the lowest (best) average D_KL. In addition to having a high D_KL, the V9 region in both the forward and reverse directions performed the worst finding only 17% and 53% of the known family level and 12% and 47% of the genus level bacteria, while results from the forward and reverse V4 region identified all 17 family level bacteria.

Conclusions

The results of our analysis have shown that our sequencing methods using 6 hypervariable regions of the 16S rRNA and subsequent analysis is valid. This method also allowed for the assessment of how well each of the variable regions might perform simultaneously. Our findings will provide the basis for future work intended to assess microbial abundance at different time points throughout a clinical protocol.     

Comparison of Species Richness Estimates Obtained Using Nearly Complete Fragments and Simulated Pyrosequencing-Generated Fragments in 16S rRNA Gene-Based Environmental Surveys

Pyrosequencing-based 16S rRNA gene surveys are increasingly utilized to study highly diverse bacterial communities, with special emphasis on utilizing the large number of sequences obtained (tens to hundreds of thousands) for species richness estimation. However, it is not yet clear how the number of operational taxonomic units (OTUs) and, hence, species richness estimates determined using shorter fragments at different taxonomic cutoffs correlates with the number of OTUs assigned using longer, nearly complete 16S rRNA gene fragments. We constructed a 16S rRNA clone library from an undisturbed tallgrass prairie soil (1,132 clones) and used it to compare species richness estimates obtained using eight pyrosequencing candidate fragments (99 to 361 bp in length) and the nearly full-length fragment. Fragments encompassing the V1 and V2 (V1+V2) region and the V6 region (generated using primer pairs 8F-338R and 967F-1046R) overestimated species richness; fragments encompassing the V3, V7, and V7+V8 hypervariable regions (generated using primer pairs 338F-530R, 1046F-1220R, and 1046F-1392R) underestimated species richness; and fragments encompassing the V4, V5+V6, and V6+V7 regions (generated using primer pairs 530F-805R, 805F-1046R, and 967F-1220R) provided estimates comparable to those obtained with the nearly full-length fragment. These patterns were observed regardless of the alignment method utilized or the parameter used to gauge comparative levels of species richness (number of OTUs observed, slope of scatter plots of pairwise distance values for short and nearly complete fragments, and nonparametric and parametric species richness estimates). Similar results were obtained when analyzing three other datasets derived from soil, adult Zebrafish gut, and basaltic formations in the East Pacific Rise. Regression analysis indicated that these observed discrepancies in species richness estimates within various regions could readily be explained by the proportions of hypervariable, variable, and conserved base pairs within an examined fragment.Culture-independent 16S rRNA gene surveys are now routinely utilized to examine the microbial diversity in various environmental habitats. However, in surveys of highly diverse ecosystems, the size of clone libraries typically constructed (100 to 500 clones) allows for the identification only of members of the community that are present in high abundance (2, 13, 14, 17, 24, 51). In addition to the failure to detect the rare members of the ecosystem, these relatively small datasets provide inaccurate estimates when used for computing species richness within an ecosystem. Regardless of the approach utilized to estimate species richness, the estimates obtained are highly dependent on sample size, and smaller datasets typically result in the underestimation of species richness (14, 44, 47, 55).The use of a pyrosequencing-based approach (40) in 16S gene-based diversity surveys promises to overcome both of the above-mentioned problems associated with inadequate sampling. The large number of 16S rRNA gene sequences produced (hundreds of thousands) allows access to rare members of the community (25; J. M. Tiedje, presented at the 108th General Meeting of the American Society for Microbiology, Boston, MA, 2008), as well as a relatively more accurate estimation of species richness. However, with the introduction of this new technology, it is necessary to correlate the results obtained from newer pyrosequencing-based surveys to the extensive collection of longer, capillary sequence-generated 16S rRNA gene sequences that has been deposited in public databases during the last 2 decades. Several recent studies have examined the utility of pyrosequencing fragments in providing an accurate survey of overall community structure (36) and investigated the ability of various fragments spanning the 16S rRNA gene to accurately predict the phylogenetic affiliation of pyrosequencing-generated fragments at various taxonomic cutoffs (35, 54). As such, these admirable efforts gave useful insights into the advantages and limitations of the pyrosequencing approach in 16S-based community surveys, pinpointed specific regions that provide better phylogenetic resolution than other pyrosequencing-generated regions, and provided a quantitative assessment of binning accuracy at various empirical cutoffs.However, while issues regarding correlating phylogenies of shorter and longer fragments are actively being addressed, efforts to calibrate species richness data obtained from various pyrosequencing fragments at various taxonomic cutoffs to estimates obtained using longer 16S rRNA gene fragments are still lacking. It is unclear how pairwise distances and, hence, operational taxonomic unit (OTU) assignments and species richness estimates computed using various shorter fragments spanning various regions of the 16S rRNA gene will correlate to pairwise distances computed using the nearly complete 16S rRNA gene. Elucidating such differences between shorter and nearly complete fragments, as well as between shorter fragments representing different regions in the 16S rRNA gene, is absolutely necessary for accurate meta-analysis of species richness in previously published and future datasets constructed using various sequencing approaches.Here, we constructed, sequenced, and analyzed a 16S rRNA library of 1,132 clones generated from an undisturbed tallgrass prairie soil in central Oklahoma and compared the numbers of OTUs and species richness values obtained using the full-length data sets (with and without the application of the Lane mask filter that excludes hypervariable regions from the phylogenetic analysis) (32) and fragments simulating pyrosequencing output generated by clipping where known conserved bacterial primers are encountered in the 16S rRNA gene. The lengths of the chosen simulated-pyrosequencing fragments represent amplicons that have been generated using the original GS20 pyrosequencing platform (≈100 bp) (25, 44, 48), similar to those currently being generated using the GS FLX pyrosequencing platform (≈250 bp) (1, 20, 35) or amplicons produced using the anticipated increase in the new GS XLR pyrosequencing platform (>250 bp). We show that the choice of the pyrosequenced fragment could indeed impact the number of OTUs calculated at different taxonomic cutoffs, with some fragments underestimating and others overestimating such parameters compared to the results with longer, nearly complete 16S rRNA gene fragments. We also show that even more marked differences could be encountered when comparing two pyrosequencing fragments within the same molecule. Further, we established a regression analysis that explains the nature of the observed discrepancies using the proportions of the hypervariable, variable, and conserved bases within fragments.     

Potential accumulative effect of the herbicide glyphosate on glyphosate-tolerant maize rhizobacterial communities over a three-year cultivation period

Background

Glyphosate is a herbicide that is liable to be used in the extensive cultivation of glyphosate-tolerant cultivars. The potential accumulation of the relative effect of glyphosate on the rhizobacterial communities of glyphosate-tolerant maize has been monitored over a period of three years.

Methodology/Principal Findings

The composition of rhizobacterial communities is known to vary with soil texture, hence, the analyses have been performed in two agricultural fields with a different soil texture. The accumulative effects of glyphosate have been monitored by means of high throughput DNA pyrosequencing of the bacterial DNA coding for the 16S rRNA hypervariable V6 region from rhizobacterial communities. The relative composition of the rhizobacterial communities does vary in each field over the three-year period. The overall distribution of the bacterial phyla seems to change from one year to the next similarly in the untreated and glyphosate-treated soils in both fields. The two methods used to estimate bacterial diversity offered consistent results and are equally suitable for diversity assessment.

Conclusions/Significance

The glyphosate treatment during the three-year period of seasonal cultivation in two different fields did not seem to significantly change the maize rhizobacterial communities when compared to those of the untreated soil. This may be particularly relevant with respect to a potential authorisation to cultivate glyphosate-tolerant maize in the European Union.     

Oral microbiome profiles: 16S rRNA pyrosequencing and microarray assay comparison

Objectives

The human oral microbiome is potentially related to diverse health conditions and high-throughput technology provides the possibility of surveying microbial community structure at high resolution. We compared two oral microbiome survey methods: broad-based microbiome identification by 16S rRNA gene sequencing and targeted characterization of microbes by custom DNA microarray.

Methods

Oral wash samples were collected from 20 individuals at Memorial Sloan-Kettering Cancer Center. 16S rRNA gene survey was performed by 454 pyrosequencing of the V3–V5 region (450 bp). Targeted identification by DNA microarray was carried out with the Human Oral Microbe Identification Microarray (HOMIM). Correlations and relative abundance were compared at phylum and genus level, between 16S rRNA sequence read ratio and HOMIM hybridization intensity.

Results

The major phyla, Firmicutes, Proteobacteria, Bacteroidetes, Actinobacteria, and Fusobacteria were identified with high correlation by the two methods (r = 0.70∼0.86). 16S rRNA gene pyrosequencing identified 77 genera and HOMIM identified 49, with 37 genera detected by both methods; more than 98% of classified bacteria were assigned in these 37 genera. Concordance by the two assays (presence/absence) and correlations were high for common genera (Streptococcus, Veillonella, Leptotrichia, Prevotella, and Haemophilus; Correlation = 0.70–0.84).

Conclusion

Microbiome community profiles assessed by 16S rRNA pyrosequencing and HOMIM were highly correlated at the phylum level and, when comparing the more commonly detected taxa, also at the genus level. Both methods are currently suitable for high-throughput epidemiologic investigations relating identified and more common oral microbial taxa to disease risk; yet, pyrosequencing may provide a broader spectrum of taxa identification, a distinct sequence-read record, and greater detection sensitivity.     

Community Analysis of Chronic Wound Bacteria Using 16S rRNA Gene-Based Pyrosequencing: Impact of Diabetes and Antibiotics on Chronic Wound Microbiota

Background

Bacterial colonization is hypothesized to play a pathogenic role in the non-healing state of chronic wounds. We characterized wound bacteria from a cohort of chronic wound patients using a 16S rRNA gene-based pyrosequencing approach and assessed the impact of diabetes and antibiotics on chronic wound microbiota.

Methodology/Principal Findings

We prospectively enrolled 24 patients at a referral wound center in Baltimore, MD; sampled patients'' wounds by curette; cultured samples under aerobic and anaerobic conditions; and pyrosequenced the 16S rRNA V3 hypervariable region. The 16S rRNA gene-based analyses revealed an average of 10 different bacterial families in wounds—approximately 4 times more than estimated by culture-based analyses. Fastidious anaerobic bacteria belonging to the Clostridiales family XI were among the most prevalent bacteria identified exclusively by 16S rRNA gene-based analyses. Community-scale analyses showed that wound microbiota from antibiotic treated patients were significantly different from untreated patients (p = 0.007) and were characterized by increased Pseudomonadaceae abundance. These analyses also revealed that antibiotic use was associated with decreased Streptococcaceae among diabetics and that Streptococcaceae was more abundant among diabetics as compared to non-diabetics.

Conclusions/Significance

The 16S rRNA gene-based analyses revealed complex bacterial communities including anaerobic bacteria that may play causative roles in the non-healing state of some chronic wounds. Our data suggest that antimicrobial therapy alters community structure—reducing some bacteria while selecting for others.     

Assessment of Species Diversity and Distribution of an Ancient Diatom Lineage Using a DNA Metabarcoding Approach

Background

Continuous efforts to estimate actual diversity and to trace the species distribution and ranges in the natural environments have gone in equal pace with advancements of the technologies in the study of microbial species diversity from microscopic observations to DNA-based barcoding. DNA metabarcoding based on Next Generation Sequencing (NGS) constitutes the latest advancement in these efforts. Here we use NGS data from different sites to investigate the geographic range of six species of the diatom family Leptocylindraceae and to identify possible new taxa within the family.

Methodology/Principal Findings

We analysed the V4 and V9 regions of the nuclear-encoded SSU rDNA gene region in the NGS database of the European ERA-Biodiversa project BioMarKs, collected in plankton and sediments at six coastal sites in European coastal waters, as well as environmental sequences from the NCBI database. All species known in the family Leptocylindraceae were detected in both datasets, but the much larger Illumina V9 dataset showed a higher species coverage at the various sites than the 454 V4 dataset. Sequences identical or similar to the references of Leptocylindrus aporus, L. convexus, L. danicus/hargravesii and Tenuicylindrus belgicus were found in the Mediterranean Sea, North Atlantic Ocean and Black Sea as well as at locations outside Europe. Instead, sequences identical or close to that of L. minimus were found in the North Atlantic Ocean and the Black Sea but not in the Mediterranean Sea, while sequences belonging to a yet undescribed taxon were encountered only in Oslo Fjord and Baffin Bay.

Conclusions/Significance

Identification of Leptocylindraceae species in NGS datasets has expanded our knowledge of the species biogeographic distribution and of the overall diversity of this diatom family. Individual species appear to be widespread, but not all of them are found everywhere. Despite the sequencing depth allowed by NGS and the wide geographic area covered by this study, the diversity of this ancient diatom family appears to be low, at least at the level of the marker used in this study.     

The Skin Microbiome in Healthy and Allergic Dogs

Background

Changes in the microbial populations on the skin of animals have traditionally been evaluated using conventional microbiology techniques. The sequencing of bacterial 16S rRNA genes has revealed that the human skin is inhabited by a highly diverse and variable microbiome that had previously not been demonstrated by culture-based methods. The goals of this study were to describe the microbiome inhabiting different areas of the canine skin, and to compare the skin microbiome of healthy and allergic dogs.

Methodology/Principal Findings

DNA extracted from superficial skin swabs from healthy (n = 12) and allergic dogs (n = 6) from different regions of haired skin and mucosal surfaces were used for 454-pyrosequencing of the 16S rRNA gene. Principal coordinates analysis revealed clustering for the different skin sites across all dogs, with some mucosal sites and the perianal regions clustering separately from the haired skin sites. The rarefaction analysis revealed high individual variability between samples collected from healthy dogs and between the different skin sites. Higher species richness and microbial diversity were observed in the samples from haired skin when compared to mucosal surfaces or mucocutaneous junctions. In all examined regions, the most abundant phylum and family identified in the different regions of skin and mucosal surfaces were Proteobacteria and Oxalobacteriaceae. The skin of allergic dogs had lower species richness when compared to the healthy dogs. The allergic dogs had lower proportions of the Betaproteobacteria Ralstonia spp. when compared to the healthy dogs.

Conclusions/Significance

The study demonstrates that the skin of dogs is inhabited by much more rich and diverse microbial communities than previously thought using culture-based methods. Our sequence data reveal high individual variability between samples collected from different patients. Differences in species richness was also seen between healthy and allergic dogs, with allergic dogs having lower species richness when compared to healthy dogs.     

Diversity of 23S rRNA Genes within Individual Prokaryotic Genomes

Background

The concept of ribosomal constraints on rRNA genes is deduced primarily based on the comparison of consensus rRNA sequences between closely related species, but recent advances in whole-genome sequencing allow evaluation of this concept within organisms with multiple rRNA operons.

Methodology/Principal Findings

Using the 23S rRNA gene as an example, we analyzed the diversity among individual rRNA genes within a genome. Of 184 prokaryotic species containing multiple 23S rRNA genes, diversity was observed in 113 (61.4%) genomes (mean 0.40%, range 0.01%–4.04%). Significant (1.17%–4.04%) intragenomic variation was found in 8 species. In 5 of the 8 species, the diversity in the primary structure had only minimal effect on the secondary structure (stem versus loop transition). In the remaining 3 species, the diversity significantly altered local secondary structure, but the alteration appears minimized through complex rearrangement. Intervening sequences (IVS), ranging between 9 and 1471 nt in size, were found in 7 species. IVS in Deinococcus radiodurans and Nostoc sp. encode transposases. T. tengcongensis was the only species in which intragenomic diversity >3% was observed among 4 paralogous 23S rRNA genes.

Conclusions/Significance

These findings indicate tight ribosomal constraints on individual 23S rRNA genes within a genome. Although classification using primary 23S rRNA sequences could be erroneous, significant diversity among paralogous 23S rRNA genes was observed only once in the 184 species analyzed, indicating little overall impact on the mainstream of 23S rRNA gene-based prokaryotic taxonomy.     

Composition of the Vaginal Microbiota in Women of Reproductive Age – Sensitive and Specific Molecular Diagnosis of Bacterial Vaginosis Is Possible?

Background and Objective

Bacterial vaginosis (BV) is the most common vaginal disorder, characterized by depletion of the normal lactobacillus-dominant microbiota and overgrowth of commensal anaerobic bacteria. This study aimed to investigate the composition of the vaginal microbiota in women of reproductive age (healthy women and women with BV), with the view of developing molecular criteria for BV diagnosis.

Materials and Methods

Vaginal samples from 163 women (79 control, 73 BV and 11 intermediate (Lactobacillary grade II flora) cases) were analyzed using 454 pyrosequencing of the hypervariable regions V3–V4 of the 16S rRNA gene and 16 quantitative bacterial species/genus-specific real-time PCR assays. Sensitivities and specificities of potential BV markers were computed using the Amsel criteria as reference standard for BV. The use of quantitative thresholds for prediction of BV, determined for both relative abundance measured with 454 pyrosequencing and bacterial load measured with qPCR, was evaluated.

Results

Relative to the healthy women, the BV patients had in their vaginal microbiota significantly higher prevalence, loads and relative abundances of the majority of BV associated bacteria. However, only Gardnerella vaginalis, Atopobium vaginae, Eggerthella, Prevotella, BVAB2 and Megasphaera type 1 detected at or above optimal thresholds were highly predictable for BV, with the best diagnostic accuracy shown for A. vaginae. The depletion of Lactobacillus species combined with the presence of either G. vaginalis or A. vaginae at diagnostic levels was a highly accurate BV predictor.

Conclusions

Quantitative determination of the presence of G. vaginalis, A. vaginae, Eggerthella, Prevotella, BVAB2 and Megasphaera type 1 as well as the depletion of Lactobacillus was highly accurate for BV diagnosis. Measurements of abundance of normal and BV microbiota relative to total bacteria in vaginal fluid may provide more accurate BV diagnosis, and be used for test-of-cure, rather than qualitative detection or absolute counts of BV related microorganisms.     

Molecular detection of eukaryotes in a single human stool sample from Senegal

Background

Microbial eukaryotes represent an important component of the human gut microbiome, with different beneficial or harmful roles; some species are commensal or mutualistic, whereas others are opportunistic or parasitic. The diversity of eukaryotes inhabiting humans remains relatively unexplored because of either the low abundance of these organisms in human gut or because they have received limited attention from a whole-community perspective.

Methodology/Principal Finding

In this study, a single fecal sample from a healthy African male was studied using both culture-dependent methods and extended molecular methods targeting the 18S rRNA and ITS sequences. Our results revealed that very few fungi, including Candida spp., Galactomyces spp., and Trichosporon asahii, could be isolated using culture-based methods. In contrast, a relatively a high number of eukaryotic species could be identified in this fecal sample when culture-independent methods based on various primer sets were used. A total of 27 species from one sample were found among the 977 analyzed clones. The clone libraries were dominated by fungi (716 clones/977, 73.3%), corresponding to 16 different species. In addition, 187 sequences out of 977 (19.2%) corresponded to 9 different species of plants; 59 sequences (6%) belonged to other micro-eukaryotes in the gut, including Entamoeba hartmanni and Blastocystis sp; and only 15 clones/977 (1.5%) were related to human 18S rRNA sequences.

Conclusion

Our results revealed a complex eukaryotic community in the volunteer’s gut, with fungi being the most abundant species in the stool sample. Larger investigations are needed to assess the generality of these results and to understand their roles in human health and disease.     

Landscape Patterns in Rainforest Phylogenetic Signal: Isolated Islands of Refugia or Structured Continental Distributions?

Objectives

Identify patterns of change in species distributions, diversity, concentrations of evolutionary history, and assembly of Australian rainforests.

Methods

We used the distribution records of all known rainforest woody species in Australia across their full continental extent. These were analysed using measures of species richness, phylogenetic diversity (PD), phylogenetic endemism (PE) and phylogenetic structure (net relatedness index; NRI). Phylogenetic structure was assessed using both continental and regional species pools. To test the influence of growth-form, freestanding and climbing plants were analysed independently, and in combination.

Results

Species richness decreased along two generally orthogonal continental axes, corresponding with wet to seasonally dry and tropical to temperate habitats. The PE analyses identified four main areas of substantially restricted phylogenetic diversity, including parts of Cape York, Wet Tropics, Border Ranges, and Tasmania. The continental pool NRI results showed evenness (species less related than expected by chance) in groups of grid cells in coastally aligned areas of species rich tropical and sub-tropical rainforest, and in low diversity moist forest areas in the south-east of the Great Dividing Range and in Tasmania. Monsoon and drier vine forests, and moist forests inland from upland refugia showed phylogenetic clustering, reflecting lower diversity and more relatedness. Signals for evenness in Tasmania and clustering in northern monsoon forests weakened in analyses using regional species pools. For climbing plants, values for NRI by grid cell showed strong spatial structuring, with high diversity and PE concentrated in moist tropical and subtropical regions.

Conclusions/Significance

Concentrations of rainforest evolutionary history (phylo-diversity) were patchily distributed within a continuum of species distributions. Contrasting with previous concepts of rainforest community distribution, our findings of continuous distributions and continental connectivity have significant implications for interpreting rainforest evolutionary history and current day ecological processes, and for managing rainforest diversity in changing circumstances.     

Analysis of 16S rRNA amplicon sequencing options on the Roche/454 next-generation titanium sequencing platform

Background

16S rRNA gene pyrosequencing approach has revolutionized studies in microbial ecology. While primer selection and short read length can affect the resulting microbial community profile, little is known about the influence of pyrosequencing methods on the sequencing throughput and the outcome of microbial community analyses. The aim of this study is to compare differences in output, ease, and cost among three different amplicon pyrosequencing methods for the Roche/454 Titanium platform

Methodology/Principal Findings

The following three pyrosequencing methods for 16S rRNA genes were selected in this study: Method-1 (standard method) is the recommended method for bi-directional sequencing using the LIB-A kit; Method-2 is a new option designed in this study for unidirectional sequencing with the LIB-A kit; and Method-3 uses the LIB-L kit for unidirectional sequencing. In our comparison among these three methods using 10 different environmental samples, Method-2 and Method-3 produced 1.5–1.6 times more useable reads than the standard method (Method-1), after quality-based trimming, and did not compromise the outcome of microbial community analyses. Specifically, Method-3 is the most cost-effective unidirectional amplicon sequencing method as it provided the most reads and required the least effort in consumables management.

Conclusions

Our findings clearly demonstrated that alternative pyrosequencing methods for 16S rRNA genes could drastically affect sequencing output (e.g. number of reads before and after trimming) but have little effect on the outcomes of microbial community analysis. This finding is important for both researchers and sequencing facilities utilizing 16S rRNA gene pyrosequencing for microbial ecological studies.     

Fine-Scale Bacterial Beta Diversity within a Complex Ecosystem (Zodletone Spring,OK, USA): The Role of the Rare Biosphere

Background

The adaptation of pyrosequencing technologies for use in culture-independent diversity surveys allowed for deeper sampling of ecosystems of interest. One extremely well suited area of interest for pyrosequencing-based diversity surveys that has received surprisingly little attention so far, is examining fine scale (e.g. micrometer to millimeter) beta diversity in complex microbial ecosystems.

Methodology/Principal Findings

We examined the patterns of fine scale Beta diversity in four adjacent sediment samples (1mm apart) from the source of an anaerobic sulfide and sulfur rich spring (Zodletone spring) in southwestern Oklahoma, USA. Using pyrosequencing, a total of 292,130 16S rRNA gene sequences were obtained. The beta diversity patterns within the four datasets were examined using various qualitative and quantitative similarity indices. Low levels of Beta diversity (high similarity indices) were observed between the four samples at the phylum-level. However, at a putative species (OTU_0.03) level, higher levels of beta diversity (lower similarity indices) were observed. Further examination of beta diversity patterns within dominant and rare members of the community indicated that at the putative species level, beta diversity is much higher within rare members of the community. Finally, sub-classification of rare members of Zodletone spring community based on patterns of novelty and uniqueness, and further examination of fine scale beta diversity of each of these subgroups indicated that members of the community that are unique, but non novel showed the highest beta diversity within these subgroups of the rare biosphere.

Conclusions/Significance

The results demonstrate the occurrence of high inter-sample diversity within seemingly identical samples from a complex habitat. We reason that such unexpected diversity should be taken into consideration when exploring gamma diversity of various ecosystems, as well as planning for sequencing-intensive metagenomic surveys of highly complex ecosystems.     

Prospecting Environmental Mycobacteria: Combined Molecular Approaches Reveal Unprecedented Diversity

Background

Environmental mycobacteria (EM) include species commonly found in various terrestrial and aquatic environments, encompassing animal and human pathogens in addition to saprophytes. Approximately 150 EM species can be separated into fast and slow growers based on sequence and copy number differences of their 16S rRNA genes. Cultivation methods are not appropriate for diversity studies; few studies have investigated EM diversity in soil despite their importance as potential reservoirs of pathogens and their hypothesized role in masking or blocking M. bovis BCG vaccine.

Methods

We report here the development, optimization and validation of molecular assays targeting the 16S rRNA gene to assess diversity and prevalence of fast and slow growing EM in representative soils from semi tropical and temperate areas. New primer sets were designed also to target uniquely slow growing mycobacteria and used with PCR-DGGE, tag-encoded Titanium amplicon pyrosequencing and quantitative PCR.

Results

PCR-DGGE and pyrosequencing provided a consensus of EM diversity; for example, a high abundance of pyrosequencing reads and DGGE bands corresponded to M. moriokaense, M. colombiense and M. riyadhense. As expected pyrosequencing provided more comprehensive information; additional prevalent species included M. chlorophenolicum, M. neglectum, M. gordonae, M. aemonae. Prevalence of the total Mycobacterium genus in the soil samples ranged from 2.3×10⁷ to 2.7×10⁸ gene targets g⁻¹; slow growers prevalence from 2.9×10⁵ to 1.2×10⁷ cells g⁻¹.

Conclusions

This combined molecular approach enabled an unprecedented qualitative and quantitative assessment of EM across soil samples. Good concordance was found between methods and the bioinformatics analysis was validated by random resampling. Sequences from most pathogenic groups associated with slow growth were identified in extenso in all soils tested with a specific assay, allowing to unmask them from the Mycobacterium whole genus, in which, as minority members, they would have remained undetected.     

Increased Throughput by Parallelization of Library Preparation for Massive Sequencing

Background

Massively parallel sequencing systems continue to improve on data output, while leaving labor-intensive library preparations a potential bottleneck. Efforts are currently under way to relieve the crucial and time-consuming work to prepare DNA for high-throughput sequencing.

Methodology/Principal Findings

In this study, we demonstrate an automated parallel library preparation protocol using generic carboxylic acid-coated superparamagnetic beads and polyethylene glycol precipitation as a reproducible and flexible method for DNA fragment length separation. With this approach the library preparation for DNA sequencing can easily be adjusted to a desired fragment length. The automated protocol, here demonstrated using the GS FLX Titanium instrument, was compared to the standard manual library preparation, showing higher yield, throughput and great reproducibility. In addition, 12 libraries were prepared and uniquely tagged in parallel, and the distribution of sequence reads between these indexed samples could be improved using quantitative PCR-assisted pooling.

Conclusions/Significance

We present a novel automated procedure that makes it possible to prepare 36 indexed libraries per person and day, which can be increased to up to 96 libraries processed simultaneously. The yield, speed and robust performance of the protocol constitute a substantial improvement to present manual methods, without the need of extensive equipment investments. The described procedure enables a considerable efficiency increase for small to midsize sequencing centers.