High-resolution microbial community reconstruction by integrating short reads from multiple 16S rRNA regions

The emergence of massively parallel sequencing technology has revolutionized microbial profiling, allowing the unprecedented comparison of microbial diversity across time and space in a wide range of host-associated and environmental ecosystems. Although the high-throughput nature of such methods enables the detection of low-frequency bacteria, these advances come at the cost of sequencing read length, limiting the phylogenetic resolution possible by current methods. Here, we present a generic approach for integrating short reads from large genomic regions, thus enabling phylogenetic resolution far exceeding current methods. The approach is based on a mapping to a statistical model that is later solved as a constrained optimization problem. We demonstrate the utility of this method by analyzing human saliva and Drosophila samples, using Illumina single-end sequencing of a 750 bp amplicon of the 16S rRNA gene. Phylogenetic resolution is significantly extended while reducing the number of falsely detected bacteria, as compared with standard single-region Roche 454 Pyrosequencing. Our approach can be seamlessly applied to simultaneous sequencing of multiple genes providing a higher resolution view of the composition and activity of complex microbial communities.     

Phylogenetic clustering of soil microbial communities by 16S rRNA but not 16S rRNA genes

We evaluated phylogenetic clustering of bacterial and archaeal communities from redox-dynamic subtropical forest soils that were defined by 16S rRNA and rRNA gene sequences. We observed significant clustering for the RNA-based communities but not the DNA-based communities, as well as increasing clustering over time of the highly active taxa detected by only rRNA.     

Flow cytometry-assisted cloning of specific sequence motifs from complex 16S rRNA gene libraries

A flow cytometry method was developed for rapid screening and recovery of cloned DNA containing common sequence motifs. This approach, termed fluorescence-activated cell sorting-assisted cloning, was used to recover sequences affiliated with a unique lineage within the Bacteroidetes not abundant in a clone library of environmental 16S rRNA genes.     

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, CX(D), and heterologous coverage curves, CXY(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.     

Explorative multivariate analyses of 16S rRNA gene data from microbial communities in modified-atmosphere-packed salmon and coalfish

Modified-atmosphere packaging (MAP) of foods in combination with low-temperature storage extends product shelf life by limiting microbial growth. We investigated the microbial biodiversity of MAP salmon and coalfish by using an explorative approach and analyzing both the total amounts of bacteria and the microbial group composition (both aerobic and anaerobic bacteria). Real-time PCR analyses revealed a surprisingly large difference in the microbial loads for the different fish samples. The microbial composition was determined by examining partial 16S rRNA gene sequences from 180 bacterial isolates, as well as by performing terminal restriction fragment length polymorphism analysis and cloning 92 sequences from PCR products of DNA directly retrieved from the fish matrix. Twenty different bacterial groups were identified. Partial least-squares (PLS) regression was used to relate the major groups of bacteria identified to the fish matrix and storage time. A strong association of coalfish with Photobacterium phosphoreum was observed. Brochothrix spp. and Carnobacterium spp., on the other hand, were associated with salmon. These bacteria dominated the fish matrixes after a storage period. Twelve Carnobacterium isolates were identified as either Carnobacterium piscicola (five isolates) or Carnobacterium divergens (seven isolates), while the eight Brochothrix isolates were identified as Brochothrix thermosphacta by full-length 16S rRNA gene sequencing. Principal-component analyses and PLS analysis of the growth characteristics (with 49 different substrates) showed that C. piscicola had distinct substrate requirements, while the requirements of B. thermosphacta and C. piscicola were quite divergent. In conclusion, our explorative multivariate approach gave a picture of the total microbial biodiversity in MAP fish that was more comprehensive than the picture that could be obtained previously. Such information is crucial in controlled food production when, for example, the hazard analysis of critical control points principle is used.     

Archaeal communities in mangrove soil characterized by 16S rRNA gene clones

An archaeal 16S rRNA gene library was constructed from mangrove soil. Phylogenetic analysis revealed archaea in mangrove soil including the Crenarchaeota (80.4%) and Euryarchaeota (19.6%) phyla. The archaeal community in mangrove soil appears to be a mixture of organisms found in a variety of environments with the majority being of marine origin.     

Alignment-independent comparisons of human gastrointestinal tract microbial communities in a multidimensional 16S rRNA gene evolutionary space

We present a novel approach for comparing 16S rRNA gene clone libraries that is independent of both DNA sequence alignment and definition of bacterial phylogroups. These steps are the major bottlenecks in current microbial comparative analyses. We used direct comparisons of taxon density distributions in an absolute evolutionary coordinate space. The coordinate space was generated by using alignment-independent bilinear multivariate modeling. Statistical analyses for clone library comparisons were based on multivariate analysis of variance, partial least-squares regression, and permutations. Clone libraries from both adult and infant gastrointestinal tract microbial communities were used as biological models. We reanalyzed a library consisting of 11,831 clones covering complete colons from three healthy adults in addition to a smaller 390-clone library from infant feces. We show that it is possible to extract detailed information about microbial community structures using our alignment-independent method. Our density distribution analysis is also very efficient with respect to computer operation time, meeting the future requirements of large-scale screenings to understand the diversity and dynamics of microbial communities.     

RNA microarray for estimating relative abundance of 16S rRNA in microbial communities

We developed an RNA microarray protocol in which total RNA from a microbial community was attached to a slide glass, and rRNA was detected by fluorescently labeled oligonucleotide probes. The RNA microarray requires only 4 h for hybridization and enables double staining and estimating relative abundance of rRNA.     

Composition of Acridid gut bacterial communities as revealed by 16S rRNA gene analysis

The gut bacterial community from four species of feral locusts and grasshoppers was determined by denaturing gradient gel electrophoresis (DGGE) analysis of bacterial 16S rRNA gene fragments. The study revealed an effect of phase polymorphism on gut bacterial diversity in brown locusts from South Africa. A single bacterial phylotype, consistent with Citrobacter sp. dominated the gut microbiota of two sympatric populations of Moroccan and Italian locusts in Spain. There was evidence for Wollbachia sp. in the meadow grasshopper caught locally in the UK. Sequence analysis of DGGE products did not reveal evidence for unculturable bacteria and homologies suggested that bacterial species were principally Gammaproteobacteria from the family Enterobacteriaceae similar to those recorded previously in laboratory reared locusts.     

Effect of phenylurea herbicides on soil microbial communities estimated by analysis of 16S rRNA gene fingerprints and community-level physiological profiles

The effect of three phenyl urea herbicides (diuron, linuron, and chlorotoluron) on soil microbial communities was studied by using soil samples with a 10-year history of treatment. Denaturing gradient gel electrophoresis (DGGE) was used for the analysis of 16S rRNA genes (16S rDNA). The degree of similarity between the 16S rDNA profiles of the communities was quantified by numerically analysing the DGGE band patterns. Similarity dendrograms showed that the microbial community structures of the herbicide-treated and nontreated soils were significantly different. Moreover, the bacterial diversity seemed to decrease in soils treated with urea herbicides, and sequence determination of several DGGE fragments showed that the most affected species in the soils treated with diuron and linuron belonged to an uncultivated bacterial group. As well as the 16S rDNA fingerprints, the substrate utilization patterns of the microbial communities were compared. Principal-component analysis performed on BIOLOG data showed that the functional abilities of the soil microbial communities were altered by the application of the herbicides. In addition, enrichment cultures of the different soils in medium with the urea herbicides as the sole carbon and nitrogen source showed that there was no difference between treated and nontreated soil in the rate of transformation of diuron and chlorotoluron but that there was a strong difference in the case of linuron. In the enrichment cultures with linuron-treated soil, linuron disappeared completely after 1 week whereas no significant transformation was observed in cultures inoculated with nontreated soil even after 4 weeks. In conclusion, this study showed that both the structure and metabolic potential of soil microbial communities were clearly affected by a long-term application of urea herbicides.     

16S rRNA gene profiling of rhizospheric microbial community of Eichhornia crassipes

Molecular Biology Reports - The rhizosphere of a plant is an important interface for the plant-microbe interaction that plays a significant role in the uptake and removal of heavy metal from...     

基于16S rRNA基因测序分析微生物群落多样性

微生物群落多样性的研究对于挖掘微生物资源,探索微生物群落功能,阐明微生物群落与生境间的关系具有重要意义。随着宏基因组概念的提出以及测序技术的快速发展,16S rRNA基因测序在微生物群落多样性的研究中已被广泛应用。文中系统地介绍了16S rRNA基因测序分析流程中的四个重要环节,包括测序平台与扩增区的选择、测序数据预处理以及多样性分析方法,就其面临的问题与挑战进行了探讨并对未来的研究方向进行了展望,以期为微生物群落多样性相关研究提供参考。     

A meta-analysis of 16S rRNA gene clone libraries from the polymicrobial black band disease of corals

Black band disease (BBD) is a polymicrobial disease affecting corals on reefs worldwide. Since 2002, researchers have constructed clone libraries from the BBD consortium using 16S rRNA gene primers targeting a variety of phyla. In the present study, a meta-analysis was conducted of published libraries from 2002 through the present that contain bacterial sequence data associated with individual clones and BBD samples. The libraries analyzed were from 87 BBD samples, collected from 16 species of scleractinian corals in 10 different geographic locations, and included 327 unique operational taxonomic units (OTUs). One OTU (cyanobacterial) was present in 62 (71%) samples, and three (one Cytophaga-Flavobacter-Bacteriodetes and two alphaproteobacterial) were present in 11 (13%) of the samples. The remaining 323 unique OTUs were present in <10% of the samples. The Alphaproteobacteria was the most diversely represented group. Analysis of clone libraries using nonmetric multidimensional scaling indicated strong regional specificity of BBD microbial populations, but limited host coral specificity. The results of this analysis support the hypotheses that: (1) a specific cyanobacterium may be the primary pathogen of BBD; (2) additional functional groups, required for BBD pathobiology, are represented by variable opportunistic species; and (3) opportunistic BBD microorganisms are primarily derived from the environment.     

Bacterial,archaeal and eukaryotic diversity in Arctic sediment as revealed by 16S rRNA and 18S rRNA gene clone libraries analysis

We studied the microbial diversity in the sediment from the Kongsfjorden, Svalbard, Arctic, in the summer of 2005 based on the analysis of 16S rRNA and 18S rRNA gene clone libraries. The sequences of the cloned 16S rRNA and 18S rRNA gene inserts were used to determine the species identity or closest relatives by comparison with sequences of known species. Compared to the other samples acquired in Arctic and Antarctic, which are different from that of ours, the microbial diversity in our sediment is much higher. The bacterial sequences were grouped into 11 major lineages of the domain Bacteria: Proteobacteria (include α-, β-, γ-, δ-, and ε-Proteobacteria); Bacteroidetes; Fusobacteria; Firmicutes; Chloroflexi; Chlamydiae; Acidobacteria; Actinobacteria; Planctomycetes; Verrucomicrobiae and Lentisphaerae. Crenarchaeota were dominant in the archaeal clones containing inserts. In addition, six groups from eukaryotes including Cercozoa, Fungi, Telonema, Stramenopiles, Alveolata, and Metazoa were identified. Remarkably, the novel group Lentisphaerae was reported in Arctic sediment at the first time. Our study suggested that Arctic sediment as a unique habitat may contain substantial microbial diversity and novel species will be discovered.     

Comparative analysis of eukaryotic marine microbial assemblages from 18S rRNA gene and gene transcript clone libraries by using different methods of extraction

Eukaryotic marine microbes play pivotal roles in biogeochemical nutrient cycling and ecosystem function, but studies that focus on the protistan biogeography and genetic diversity lag-behind studies of other microbes. 18S rRNA PCR amplification and clone library sequencing are commonly used to assess diversity that is culture independent. However, molecular methods are not without potential biases and artifacts. In this study, we compare the community composition of clone libraries generated from the same water sample collected at the San Pedro Ocean Time Series (SPOTs) station in the northwest Pacific Ocean. Community composition was assessed using different cell lysis methods (chemical and mechanical) and the extraction of different nucleic acids (DNA and RNA reverse transcribed to cDNA) to build Sanger ABI clone libraries. We describe specific biases for ecologically important phylogenetic groups resulting from differences in nucleic acid extraction methods that will inform future designs of eukaryotic diversity studies, regardless of the target sequencing platform planned.