Quantitative Comparisons of 16S rRNA Gene Sequence Libraries from Environmental Samples

To determine the significance of differences between clonal libraries of environmental rRNA gene sequences, differences between homologous coverage curves, C_X(D), and heterologous coverage curves, C_XY(D), were calculated by a Cramér-von Mises-type statistic and compared by a Monte Carlo test procedure. This method successfully distinguished rRNA gene sequence libraries from soil and bioreactors and correctly failed to find differences between libraries of the same composition.     

Effects of Transgenic Hybrid Aspen Overexpressing Polyphenol Oxidase on Rhizosphere Diversity

This study assessed the potential effects of transgenic aspen overexpressing a polyphenol oxidase gene on diversity in rhizosphere communities. Cultivation-independent methods were used to better delineate bacterial and fungal populations associated with transgenic and nontransgenic trees. Gene libraries for the bacterial component of the rhizosphere were established using 16S rRNA and chaperonin-60 (CPN-60) gene sequences, while the fungal community was characterized using 18S rRNA gene sequences. The 16S rRNA gene libraries were dominated by alphaproteobacterial sequences, while the CPN-60 gene libraries were dominated by members of the Bacteroidetes/Chlorobi group. In both the CPN-60 and 16S rRNA libraries, there were differences in only minor components of the bacterial community between transgenic and unmodified trees, and no significant differences in species diversity were observed. Compared to the bacterial gene libraries, greater coverage of the underlying population was achieved with the fungal 18S rRNA libraries. Members of the Zygomycota, Chytridiomycota, Ascomycota, and Basidiomycota were recovered from both libraries. The dominant groups of fungi associated with each tree type were very similar, although there were some qualitative differences in the recovery of less-abundant fungi, likely as a result of the underlying heterogeneity of the fungal population. The methods employed revealed only minor differences between the bacterial and fungal communities associated with transgenic and unmodified trees.     

Effect of Condensed Tannins on Bovine Rumen Protist Diversity Based on 18S rRNA Gene Sequences

Molecular diversity of protists from bovine rumen fluid incubated with condensed tannins of Leucaena leucocephala hybrid‐Rendang at 20 mg/500 mg dry matter (treatment) or without condensed tannins (control) was investigated using 18S rRNA gene library. Clones from the control library were distributed within nine genera, but clones from the condensed tannin treatment clone library were related to only six genera. Diversity estimators such as abundance‐based coverage estimation and Chao1 showed significant differences between the two libraries, although no differences were found based on Shannon–Weaver index and Libshuff.     

Low sequencing efforts bias analyses of shared taxa in microbial communities

The potential for comparing microbial community population structures has been greatly enhanced by developments in next generation sequencing methods that can generate hundreds of thousands to millions of reads in a single run. Conversely, many microbial community comparisons have been published with no more than 1,000 sequences per sample. These studies have presented data on levels of shared operational taxonomic units (OTUs) between communities. Due to lack of coverage, that approach might compromise the conclusions about microbial diversity and the degree of difference between environments. In this study, we present data from recent studies that highlight this problem. Also, we analyzed datasets of 16 rRNA sequences with small and high sequence coverage from different environments to demonstrate that the level of sequencing effort used for analyzing microbial communities biases the results. We randomly sampled pyrosequencing-generated 16S rRNA gene libraries with increasing sequence effort. Sequences were used to calculate Good's coverage, the percentage of shared OTUs, and phylogenetic distance measures. Our data showed that simple counts of presence/absence of taxonomic unities do not reflect the real similarity in membership and structure of the bacterial communities and that community comparisons based on phylogenetic tests provide a way to test statistically significant differences between two or more environments without need an exhaustive sampling effort.     

Integration of microbial ecology and statistics: a test to compare gene libraries

Libraries of 16S rRNA genes provide insight into the membership of microbial communities. Statistical methods help to determine whether differences in library composition are artifacts of sampling or are due to underlying differences in the communities from which they are derived. To contribute to a growing statistical framework for comparing 16S rRNA libraries, we present a computer program, integral -LIBSHUFF, which calculates the integral form of the Cramér-von Mises statistic. This implementation builds upon the LIBSHUFF program, which uses an approximation of the statistic and makes a number of modifications that improve precision and accuracy. Once integral -LIBSHUFF calculates the P values, when pairwise comparisons are tested at the 0.05 level, the probability of falsely identifying a significant P value is 0.098 for a study with two libraries, 0.265 for three libraries, and 0.460 for four libraries. The potential negative effects of making the multiple pairwise comparisons necessitate correcting for the increased likelihood that differences between treatments are due to chance and do not reflect biological differences. Using integral -LIBSHUFF, we found that previously published 16S rRNA gene libraries constructed from Scottish and Wisconsin soils contained different bacterial lineages. We also analyzed the published libraries constructed for the zebrafish gut microflora and found statistically significant changes in the community during development of the host. These analyses illustrate the power of integral -LIBSHUFF to detect differences between communities, providing the basis for ecological inference about the association of soil productivity or host gene expression and microbial community composition.     

Phylogenetic analysis of bacteria associated with Laminaria saccharina

Bacterial communities associated with the brown alga Laminaria saccharina from the Baltic Sea and from the North Sea were investigated using denaturing gradient gel electrophoresis and 16S rRNA gene clone libraries. The rhizoid, cauloid, meristem and phyloid revealed different 16S rRNA gene denaturing gradient gel electrophoresis banding patterns indicating a specific association of bacterial communities with different parts of the alga. Associations with cauloid and meristem were more specific, while less specific associations were obtained from the old phyloid. In addition, seasonal and geographical differences in the associated communities were observed. Results from 16S rRNA gene libraries supported these findings. Bacterial phylotypes associated with the alga were affiliated with the Alphaproteobacteria (nine phylotypes), Gammaproteobacteria (nine phylotypes) and the Bacteroidetes group (four phylotypes). A number of bacteria associated with other algae and other marine macroorganisms were among the closest relatives of phylotypes associated with L. saccharina.     

Integration of Microbial Ecology and Statistics: a Test To Compare Gene Libraries

Libraries of 16S rRNA genes provide insight into the membership of microbial communities. Statistical methods help to determine whether differences in library composition are artifacts of sampling or are due to underlying differences in the communities from which they are derived. To contribute to a growing statistical framework for comparing 16S rRNA libraries, we present a computer program, ∫-LIBSHUFF, which calculates the integral form of the Cramér-von Mises statistic. This implementation builds upon the LIBSHUFF program, which uses an approximation of the statistic and makes a number of modifications that improve precision and accuracy. Once ∫-LIBSHUFF calculates the P values, when pairwise comparisons are tested at the 0.05 level, the probability of falsely identifying a significant P value is 0.098 for a study with two libraries, 0.265 for three libraries, and 0.460 for four libraries. The potential negative effects of making the multiple pairwise comparisons necessitate correcting for the increased likelihood that differences between treatments are due to chance and do not reflect biological differences. Using ∫-LIBSHUFF, we found that previously published 16S rRNA gene libraries constructed from Scottish and Wisconsin soils contained different bacterial lineages. We also analyzed the published libraries constructed for the zebrafish gut microflora and found statistically significant changes in the community during development of the host. These analyses illustrate the power of ∫-LIBSHUFF to detect differences between communities, providing the basis for ecological inference about the association of soil productivity or host gene expression and microbial community composition.     

Investigation of the methanogen population structure and activity in a brackish lake sediment

The methanogen community in sediment from the edge of a small brackish lake connected to the Beaulieu Estuary (Hampshire, UK) was investigated by analysis of 16S rRNA gene diversity using new methanogen-specific primers plus Archaea-specific primers. 16S rRNA gene primers previously used for polymerase chain reaction (PCR) detection of methanogenic Archaea from a variety of environments were evaluated by in silico testing. The primers displayed variable coverage of the four main orders of methanogens, highlighting the importance of this type of primer evaluation. Three PCR primer sets were designed using novel reverse primers to facilitate specific amplification of the orders Methanomicrobiales/Methanosarcinales, Methanobacteriales and Methanococcales. Diversity of the methanogen functional gene, methyl coenzyme M reductase (mcrA), was also studied. All gene libraries constructed from this sediment indicated that Methanomicrobiales and Methanosarcinales were the only methanogens detected. There was good agreement between the relative sequence abundances in the methanogen-specific 16S rRNA gene library and terminal restriction fragment length polymorphism (T-RFLP) profiling, suggesting that the population was dominated by putative H2 CO2 utilizing Methanomicrobiales, although acetate-utilizing methanogens were also present. The methanogen population analyses were in agreement with methanogenic activity measurements, which indicated that bicarbonate methanogenesis was higher than acetate methanogenesis at all depths measured and overall there was a significant difference (P = 0.001) between the rates of the two pathways. This study demonstrates the utility of new 16S rRNA gene PCR primers targeting specific methanogenic orders, and the combined results suggest that the CO2 reduction pathway dominates methanogenesis in the brackish sediment investigated.     

Geographical distribution and diversity of bacteria associated with natural populations of Drosophila melanogaster

Drosophila melanogaster is one of the most widely used model systems in biology. However, little is known about its associated bacterial community. As a first step towards understanding these communities, we compared bacterial 16S rRNA gene sequence libraries recovered from 11 natural populations of adult D. melanogaster. Bacteria from these sequence libraries were grouped into 74 distinct taxa, spanning the phyla Proteobacteria, Bacteroidetes, and Firmicutes, which were unevenly spread across host populations. Summed across populations, the distribution of abundance of genera was closely fit by a power law. We observed differences among host population locations both in bacterial community richness and in composition. Despite this significant spatial variation, no relationship was observed between species richness and a variety of abiotic factors, such as temperature and latitude. Overall, bacterial communities associated with adult D. melanogaster hosts are diverse and differ across host populations.     

Molecular characterization of microbial communities in deep coal seam groundwater of northern Japan

We investigated microbial methanogenesis and community structure based on 16S rRNA gene sequences from a coal seam aquifer located 843–907 m below ground level in northern Japan; additionally, we studied the δ¹³C and δ²H (δD) of coal‐bed gases and other physicochemical parameters. Although isotopic analysis suggested a thermocatalytic origin for the gases, the microbial activity and community structure strongly implied the existence of methanogenic microbial communities in situ. Methane was generated in the enrichment cultures of the hydrogenotrophic and methylotrophic microorganisms obtained from coal seam groundwater. Methanogen clones dominated the archaeal 16S rRNA gene libraries and were mostly related to the hydrogenotrophic genus Methanoculleus and the methylotrophic genus Methanolobus. Bacterial 16S rRNA gene libraries were dominated by the clones related to the genera Acetobacterium and Syntrophus which have a symbiotic association with methanogens. LIBSHUFF analysis revealed that N₂ gas injected into the coal seam (for enhanced methane production) does not affect the coverage of archaeal and bacterial populations. However, amova analysis does provide evidence for a change in the genetic diversity of archaeal populations that are dominated by methanogens. Therefore, N₂ injection into the coal seam might affect the cycling of matter by methanogens in situ.     

Bacterial biodiversity from Roopkund Glacier, Himalayan mountain ranges, India

The bacterial diversity of two soil samples collected from the periphery of the Roopkund glacial lake and one soil sample from the surface of the Roopkund Glacier in the Himalayan ranges was determined by constructing three 16S rRNA gene clone libraries. The three clone libraries yielded a total of 798 clones belonging to 25 classes. Actinobacteria was the most predominant class (>10% of the clones) in the three libraries. In the library from the glacial soil, class Betaproteobacteria (24.2%) was the most predominant. The rarefaction analysis indicated coverage of 43.4 and 41.2% in the samples collected from the periphery of the lake thus indicating a limited bacterial diversity covered; at the same time, the coverage of 98.4% in the glacier sample indicated most of the diversity was covered. Further, the bacterial diversity in the Roopkund glacier soil was low, but was comparable with the bacterial diversity of a few other glaciers. The results of principal component analysis based on the 16S rRNA gene clone library data, percentages of OTUs and biogeochemical data revealed that the lake soil samples were different from the glacier soil sample and the biogeochemical properties affected the diversity of microbial communities in the soil samples.     

Microbial sequences retrieved from environmental samples from seasonal Arctic snow and meltwater from Svalbard, Norway

16S rRNA gene (rrs) clone libraries were constructed from two snow samples (May 11, 2007 and June 7, 2007) and two meltwater samples collected during the spring of 2007 in Svalbard, Norway (79°N). The libraries covered 19 different microbial classes, including Betaproteobacteria (21.3%), Sphingobacteria (16.4%), Flavobacteria (9.0%), Acidobacteria (7.7%) and Alphaproteobacteria (6.5%). Significant differences were detected between the two sets of sample libraries. First, the meltwater libraries had the highest community richness (Chao1: 103.2 and 152.2) and Shannon biodiversity indices (between 3.38 and 3.59), when compared with the snow libraries (Chao1: 14.8 and 59.7; Shannon index: 1.93 and 3.01). Second, ∫-LIBSHUFF analyses determined that the bacterial communities in the snow libraries were significantly different from those of the meltwater libraries. Despite these differences, our data also support the theory that a common core group of microbial populations exist within a variety of cryohabitats.     

Molecular profiling of rhizosphere microbial communities associated with healthy and diseased black spruce (Picea mariana) seedlings grown in a nursery

Bacterial and fungal populations associated with the rhizosphere of healthy black spruce (Picea mariana) seedlings and seedlings with symptoms of root rot were characterized by cloned rRNA gene sequence analysis. Triplicate bacterial and fungal rRNA gene libraries were constructed, and 600 clones were analyzed by amplified ribosomal DNA restriction analysis and grouped into operational taxonomical units (OTUs). A total of 84 different bacterial and 31 different fungal OTUs were obtained and sequenced. Phylogenetic analyses indicated that the different OTUs belonged to a wide range of bacterial and fungal taxa. For both groups, pairwise comparisons revealed that there was greater similarity between replicate libraries from each treatment than between libraries from different treatments. Significant differences between pooled triplicate samples from libraries of genes from healthy seedlings and pooled triplicate samples from libraries of genes from diseased seedlings were also obtained for both bacteria and fungi, clearly indicating that the rhizosphere-associated bacterial and fungal communities of healthy and diseased P. mariana seedlings were different. The communities associated with healthy and diseased seedlings also showed distinct ecological parameters as indicated by the calculated diversity, dominance, and evenness indices. Among the main differences observed at the community level, there was a higher proportion of Acidobacteria, Gammaproteobacteria, and Homobasidiomycetes clones associated with healthy seedlings, while the diseased-seedling rhizosphere harbored a higher proportion of Actinobacteria, Sordariomycetes, and environmental clones. The methodological approach described in this study appears promising for targeting potential rhizosphere-competent biological control agents against root rot diseases occurring in conifer nurseries.     

Molecular Profiling of Rhizosphere Microbial Communities Associated with Healthy and Diseased Black Spruce (Picea mariana) Seedlings Grown in a Nursery

Bacterial and fungal populations associated with the rhizosphere of healthy black spruce (Picea mariana) seedlings and seedlings with symptoms of root rot were characterized by cloned rRNA gene sequence analysis. Triplicate bacterial and fungal rRNA gene libraries were constructed, and 600 clones were analyzed by amplified ribosomal DNA restriction analysis and grouped into operational taxonomical units (OTUs). A total of 84 different bacterial and 31 different fungal OTUs were obtained and sequenced. Phylogenetic analyses indicated that the different OTUs belonged to a wide range of bacterial and fungal taxa. For both groups, pairwise comparisons revealed that there was greater similarity between replicate libraries from each treatment than between libraries from different treatments. Significant differences between pooled triplicate samples from libraries of genes from healthy seedlings and pooled triplicate samples from libraries of genes from diseased seedlings were also obtained for both bacteria and fungi, clearly indicating that the rhizosphere-associated bacterial and fungal communities of healthy and diseased P. mariana seedlings were different. The communities associated with healthy and diseased seedlings also showed distinct ecological parameters as indicated by the calculated diversity, dominance, and evenness indices. Among the main differences observed at the community level, there was a higher proportion of Acidobacteria, Gammaproteobacteria, and Homobasidiomycetes clones associated with healthy seedlings, while the diseased-seedling rhizosphere harbored a higher proportion of Actinobacteria, Sordariomycetes, and environmental clones. The methodological approach described in this study appears promising for targeting potential rhizosphere-competent biological control agents against root rot diseases occurring in conifer nurseries.     

Seasonal diversity of planktonic protists in Southwestern Alberta rivers over a 1-year period as revealed by terminal restriction fragment length polymorphism and 18S rRNA gene library analyses

The temporal dynamics of planktonic protists in river water have received limited attention despite their ecological significance and recent studies linking phagotrophic protists to the persistence of human-pathogenic bacteria. Using molecular-based techniques targeting the 18S rRNA gene, we studied the seasonal diversity of planktonic protists in Southwestern Alberta rivers (Oldman River Basin) over a 1-year period. Nonmetric multidimensional scaling analysis of terminal restriction fragment length polymorphism (T-RFLP) data revealed distinct shifts in protistan community profiles that corresponded to season rather than geographical location. Community structures were examined by using clone library analysis; HaeIII restriction profiles of 18S rRNA gene amplicons were used to remove prevalent solanaceous plant clones prior to sequencing. Sanger sequencing of the V1-to-V3 region of the 18S rRNA gene libraries from spring, summer, fall, and winter supported the T-RFLP results and showed marked seasonal differences in the protistan community structure. The spring library was dominated by Chloroplastidae (29.8%), Centrohelida (28.1%), and Alveolata (25.5%), while the summer and fall libraries contained primarily fungal clones (83.0% and 88.0%, respectively). Alveolata (35.6%), Euglenozoa (24.4%), Chloroplastida (15.6%), and Fungi (15.6%) dominated the winter library. These data demonstrate that planktonic protists, including protozoa, are abundant in river water in Southwestern Alberta and that conspicuous seasonal shifts occur in the community structure.     

Diversity of bacterial community in freshwater of Woopo wetland

Diversity of bacterial community in water layer of Woopo wetland was investigated. Cultivable bacterial strains were isolated by the standard dilution plating technique and culture-independent 16S rRNA gene clones were obtained directly from DNA extracts of a water sample. Amplified rDNA restriction analysis (ARDRA) was applied onto both of the isolates and 16S rRNA gene clones. Rarefaction curves, coverage rate and diversity indices of ARDRA patterns were calculated. Representative isolates and clones of all the single isolate/clone phylotype were partially sequenced and analyzed phylogenetically. Sixty-four and 125 phylotypes were obtained from 203 bacterial isolates and 235 culture-independent 16S rRNA gene clones, respectively. Bacterial isolates were composed of 4 phyla, of which Firmicutes (49.8%) and Actinobacteria (32.0%) were predominant. Isolates were affiliated with 58 species. Culture-independent 16S rRNA gene clones were composed of 8 phyla, of which Proteobacteria (62.2%), Actinobacteria (15.5%), and Bacteroidetes (13.7%) were predominant. Diversity of 16S rRNA gene clones originated from cultivation-independent DNA extracts was higher than that of isolated bacteria.     

Spatial distribution of sponge-associated bacteria in the Mediterranean sponge Tethya aurantium

The local distribution of the bacterial community associated with the marine sponge Tethya aurantium Pallas 1766 was studied. Distinct bacterial communities were found to inhabit the endosome and cortex. Clear differences in the associated bacterial populations were demonstrated by denaturing gradient gel electrophoresis (DGGE) and analysis of 16S rRNA gene clone libraries. Specifically associated phylotypes were identified for both regions: a new phylotype of Flexibacteria was recovered only from the sponge cortex, while Synechococcus species were present mainly in the sponge endosome. Light conduction via radiate spicule bundles conceivably facilitates the unusual association of Cyanobacteria with the sponge endosome. Furthermore, a new monophyletic cluster of sponge-derived 16S rRNA gene sequences related to the Betaproteobacteria was identified using analysis of 16S rRNA gene clone libraries. Members of this cluster were specifically associated with both cortex and endosome of T. aurantium.     

Picoeukaryotic Plankton Diversity at the Helgoland Time Series Site as Assessed by Three Molecular Methods

We analyzed picoeukaryote assemblages in the German Bight at the Helgoland time series site by sequencing cloned eukaryotic 18S rRNA genes in six genetic libraries plus one library from the Orkney Islands from a cruise of opportunity. The libraries were constructed from environmental samples collected at different periods of the year. The same samples were also analyzed using a fingerprinting technique, single-strand conformational polymorphism (SSCP), and DNA microarrays with class-level oligonucleotide probes. One hundred unique clones were analyzed from each library, thus insuring over 85% coverage of the library. The V4 region of the 18S rRNA gene was sequenced from each of these clones, thus providing the most discrimination among the clones. The nonphotosynthetic picoeukaryotic component dominated over the photosynthetic one and was represented by the ciliates at 45% and group II alveolates at 42%. Prasinophytes dominated the photosynthetic group at 40%, but other picoplankton groups, such as bolidomonads and chrysophytes, were also present. Totally novel groups were found in the cryptomonads and in the dinoflagellates. A new algal group sister to the cryptophyte nuclear gene and the glaucocystophytes was also found. These three groups have been found in other picoeukaryotic planktonic clone libraries. SSCP analyses at closer time intervals suggest that clone libraries should be made at weekly intervals if succession in the picoeukaryotic plankton community is to be monitored accurately. A comparison of annual samples suggests thatthere appears to be an annual cycle with regard to species composition. Microarray analysis supported the clone library data and offered a faster means of community analysis, which can be performed with similar accuracy and with higher throughput for a more in-depth analysis.     

Development of a real-time PCR method for quantification of the three genera Dehalobacter, Dehalococcoides, and Desulfitobacterium in microbial communities

We developed standard curves based on plasmids containing a 16S rRNA gene of a member of one of the three genera Dehalobacter, Desulfitobacterium, and Dehalococcoides. A large difference in amplification efficiency between the standard curves was observed ranging from 1.5 to 2.0. The total eubacterial 16S rRNA gene copy number determined in a sample DNA by using eubacterial primers and the three standard curves led to differences in the estimated copy numbers of a factor up to 73. However, the amplification efficiencies for one specific standard curve were the same independent of the PCR primer pair used. This allowed the determination of the abundance of a population expressed as fractional number, hence, the percentage of genus-specific copy numbers within the total eubacterial 16S rRNA gene copy numbers. Determination of the fractional numbers in DNA mixtures of known composition showed the accuracy of this approach. The average difference in threshold value between two 10-fold dilutions of DNA of pure cultures, mixtures thereof and of environmental samples was -3.45+/-0.34, corresponding to an average almost optimal amplification efficiency of 1.95. This indicated that the low amplification efficiency of certain standard curves seemed to be mainly a problem of the plasmid DNA used and not of the 16S rRNA gene of the target genera.     

Molecular-phylogenetic characterization of microbial communities imbalances in the small intestine of dogs with inflammatory bowel disease

An association between luminal commensal bacteria and inflammatory bowel disease (IBD) has been suggested in humans, but studies investigating the intestinal microbial communities of dogs with IBD have not been published. The aim of this study was to characterize differences of the small intestinal microbial communities between dogs with IBD and healthy control dogs. Duodenal brush cytology samples were endoscopically collected from 10 dogs with IBD and nine healthy control dogs. DNA was extracted and 16S rRNA gene was amplified using universal bacterial primers. Constructed 16S rRNA gene clone libraries were compared between groups. From a total of 1240 selected clones, 156 unique 16S rRNA gene sequences were identified, belonging to six phyla: Firmicutes (53.4%), Proteobacteria (28.4%), Bacteroidetes (7.0%), Spirochaetes (5.2%), Fusobacteria (3.4%), Actinobacteria (1.1%), and Incertae sedis (1.5%). Species richness was significantly lower in the IBD group (P=0.038). Principal component analysis indicated that the small intestinal microbial communities of IBD and control dogs are composed of distinct microbial communities. The most profound difference involved enrichment of the IBD dogs with members of the Enterobacteriaceae family. However, differences involving members of other families, such as Clostridiaceae, Bacteroidetes and Spirochaetes, were also identified. In conclusion, canine IBD is associated with altered duodenal microbial communities compared with healthy controls.