Intra-Genomic Variation in the Ribosomal Repeats of Nematodes

Ribosomal loci represent a major tool for investigating environmental diversity and community structure via high-throughput marker gene studies of eukaryotes (e.g. 18S rRNA). Since the estimation of species’ abundance is a major goal of environmental studies (by counting numbers of sequences), understanding the patterns of rRNA copy number across species will be critical for informing such high-throughput approaches. Such knowledge is critical, given that ribosomal RNA genes exist within multi-copy repeated arrays in a genome. Here we measured the repeat copy number for six nematode species by mapping the sequences from whole genome shotgun libraries against reference sequences for their rRNA repeat. This revealed a 6-fold variation in repeat copy number amongst taxa investigated, with levels of intragenomic variation ranging from 56 to 323 copies of the rRNA array. By applying the same approach to four C. elegans mutation accumulation lines propagated by repeated bottlenecking for an average of ~400 generations, we find on average a 2-fold increase in repeat copy number (rate of increase in rRNA estimated at 0.0285-0.3414 copies per generation), suggesting that rRNA repeat copy number is subject to selection. Within each Caenorhabditis species, the majority of intragenomic variation found across the rRNA repeat was observed within gene regions (18S, 28S, 5.8S), suggesting that such intragenomic variation is not a product of selection for rRNA coding function. We find that the dramatic variation in repeat copy number among these six nematode genomes would limit the use of rRNA in estimates of organismal abundance. In addition, the unique pattern of variation within a single genome was uncorrelated with patterns of divergence between species, reflecting a strong signature of natural selection for rRNA function. A better understanding of the factors that control or affect copy number in these arrays, as well as their rates and patterns of evolution, will be critical for informing estimates of global biodiversity.     

RiboFR-Seq: a novel approach to linking 16S rRNA amplicon profiles to metagenomes

16S rRNA amplicon analysis and shotgun metagenome sequencing are two main culture-independent strategies to explore the genetic landscape of various microbial communities. Recently, numerous studies have employed these two approaches together, but downstream data analyses were performed separately, which always generated incongruent or conflict signals on both taxonomic and functional classifications. Here we propose a novel approach, RiboFR-Seq (Ribosomal RNA gene flanking region sequencing), for capturing both ribosomal RNA variable regions and their flanking protein-coding genes simultaneously. Through extensive testing on clonal bacterial strain, salivary microbiome and bacterial epibionts of marine kelp, we demonstrated that RiboFR-Seq could detect the vast majority of bacteria not only in well-studied microbiomes but also in novel communities with limited reference genomes. Combined with classical amplicon sequencing and shotgun metagenome sequencing, RiboFR-Seq can link the annotations of 16S rRNA and metagenomic contigs to make a consensus classification. By recognizing almost all 16S rRNA copies, the RiboFR-seq approach can effectively reduce the taxonomic abundance bias resulted from 16S rRNA copy number variation. We believe that RiboFR-Seq, which provides an integrated view of 16S rRNA profiles and metagenomes, will help us better understand diverse microbial communities.     

Phylogenetic Profiling and Diversity of Bacterial Communities in the Death Valley,an Extreme Habitat in the Atacama Desert

The Atacama Desert, one of the driest deserts in the world, represents a unique extreme environmental ecosystem to explore the bacterial diversity as it is considered to be at the dry limit for life. A 16S rRNA gene (spanning the hyper variable V3 region) library was constructed from an alkaline sample of unvegetated soil at the hyperarid margin in the Atacama Desert. A total of 244 clone sequences were used for MOTHUR analysis, which revealed 20 unique phylotypes or operational taxonomic units (OTUs). V3 region amplicons of the 16S rRNA were suitable for distinguishing the bacterial community to the genus and specie level. We found that all OTUs were affiliated with taxa representative of the Firmicutes phylum. The extremely high abundance of Firmicutes indicated that most bacteria in the soil were spore-forming survivors. In this study we detected a narrower diversity as compared to other ecological studies performed in other areas of the Atacama Desert. The reported genera were Oceanobacillus (representing the 69.5 % of the clones sequenced), Bacillus, Thalassobacillus and Virgibacillus. The present work shows physical and chemical parameters have a prominent impact on the microbial community structure. It constitutes an example of the communities adapted to live in extreme conditions caused by dryness and metal concentrations .

Electronic supplementary material

The online version of this article (doi:10.1007/s12088-015-0539-3) contains supplementary material, which is available to authorized users.     

Discordant phylogenies within the rrn loci of Rhizobia

It is evident from complete genome sequencing results that lateral gene transfer and recombination are essential components in the evolutionary process of bacterial genomes. Since this has important implications for bacterial systematics, the primary objective of this study was to compare estimated evolutionary relationships among a representative set of alpha-Proteobacteria by sequencing analysis of three loci within their rrn operons. Tree topologies generated with 16S rRNA gene sequences were significantly different from corresponding trees assembled with 23S rRNA gene and internally transcribed space region sequences. Besides the incongruence in tree topologies, evidence that distinct segments along the 16S rRNA gene sequences of bacteria currently classified within the genera Bradyrhizobium, Mesorhizobium and Sinorhizobium have a reticulate evolutionary history was also obtained. Our data have important implications for bacterial taxonomy, because currently most taxonomic decisions are based on comparative 16S rRNA gene sequence analysis. Since phylogenetic placement based on 16S rRNA gene sequence divergence perhaps is questionable, we suggest that the proposals of bacterial nomenclature or changes in their taxonomy that have been made may not necessarily be warranted. Accordingly, a more conservative approach should be taken in the future, in which taxonomic decisions are based on the analysis of a wider variety of loci and comparative analytical methods are used to estimate phylogenetic relationships among the genomes under consideration.     

Divergence and redundancy of 16S rRNA sequences in genomes with multiple rrn operons

The level of sequence heterogeneity among rrn operons within genomes determines the accuracy of diversity estimation by 16S rRNA-based methods. Furthermore, the occurrence of widespread horizontal gene transfer (HGT) between distantly related rrn operons casts doubt on reconstructions of phylogenetic relationships. For this study, patterns of distribution of rrn copy numbers, interoperonic divergence, and redundancy of 16S rRNA sequences were evaluated. Bacterial genomes display up to 15 operons and operon numbers up to 7 are commonly found, but ~40% of the organisms analyzed have either one or two operons. Among the Archaea, a single operon appears to dominate and the highest number of operons is five. About 40% of sequences among 380 operons in 76 bacterial genomes with multiple operons were identical to at least one other 16S rRNA sequence in the same genome, and in 38% of the genomes all 16S rRNAs were invariant. For Archaea, the number of identical operons was only 25%, but only five genomes with 21 operons are currently available. These considerations suggest an upper bound of roughly threefold overestimation of bacterial diversity resulting from cloning and sequencing of 16S rRNA genes from the environment; however, the inclusion of genomes with a single rrn operon may lower this correction factor to ~2.5. Divergence among operons appears to be small overall for both Bacteria and Archaea, with the vast majority of 16S rRNA sequences showing <1% nucleotide differences. Only five genomes with operons with a higher level of nucleotide divergence were detected, and Thermoanaerobacter tengcongensis exhibited the highest level of divergence (11.6%) noted to date. Overall, four of the five extreme cases of operon differences occurred among thermophilic bacteria, suggesting a much higher incidence of HGT in these bacteria than in other groups.     

Bacterial Diversity in Linglong Gold Mine,China

Bacteria have been actively regulating cycles of various elements in the environment. To explore the potential bacterial role in gold biogeochemical cycling, this study analyzed the bacterial diversity of mine rock (MR) and surface soil (SS) samples from Linglong gold mine using 16S rRNA gene clone library analysis and cultivation method. From MR, 24 operational taxonomic units (OTUs) were identified from MR, covering 3 phyla and 18 genera. Meanwhile, 24 OTUs were identified from SS, including 4 phyla and 18 genera. Compared with 16S rRNA gene clone library analysis, 28 aerobic and 34 anaerobic isolates were obtained, whereas 26 aerobic and 71 anaerobic strains were isolated from SS. The cultivable bacteria were affiliated with Firmicutes, Proteobacteria and Actinobacteria phyla, and dominated by Firmicutes. These results underscore the high level of bacterial diversity in the gold mine. Our study provides information on the microbial diversity in Linglong gold mine and sheds light on the existence and potential function of bacteria in the gold biogeochemical cycling.     

Bacterial diversity,community structure and potential growth rates along an estuarine salinity gradient

Very little is known about growth rates of individual bacterial taxa and how they respond to environmental flux. Here, we characterized bacterial community diversity, structure and the relative abundance of 16S rRNA and 16S rRNA genes (rDNA) using pyrosequencing along the salinity gradient in the Delaware Bay. Indices of diversity, evenness, structure and growth rates of the surface bacterial community significantly varied along the transect, reflecting active mixing between the freshwater and marine ends of the estuary. There was no positive correlation between relative abundances of 16S rRNA and rDNA for the entire bacterial community, suggesting that abundance of bacteria does not necessarily reflect potential growth rate or activity. However, for almost half of the individual taxa, 16S rRNA positively correlated with rDNA, suggesting that activity did follow abundance in these cases. The positive relationship between 16S rRNA and rDNA was less in the whole water community than for free-living taxa, indicating that the two communities differed in activity. The 16S rRNA:rDNA ratios of some typically marine taxa reflected differences in light, nutrient concentrations and other environmental factors along the estuarine gradient. The ratios of individual freshwater taxa declined as salinity increased, whereas the 16S rRNA:rDNA ratios of only some typical marine bacteria increased as salinity increased. These data suggest that physical and other bottom-up factors differentially affect growth rates, but not necessarily abundance of individual taxa in this highly variable environment.     

Distribution of High Bacterial Taxa Across the Chronosequence of Two Alpine Glacier Forelands

Little is known about the changes in abundance of microbial taxa in relation to the chronosequence of receding glaciers. This study investigated how the abundances of ten bacterial phyla or classes varied along successional gradients in two glaciers, Ödenwinkelkees and Rotmoosferner, in the central Alps. Quantitative PCR was used to estimate the abundance of the different bacterial taxa in extended glacier chronosequences, including 10- to 160-year-old successional stages, the surface of the glacier, and a fully established soil. Actinobacteria (15–30%) was the dominant group within the chronosequences. Several taxa showed significant differences in the number of taxa-specific 16S rRNA gene copies per nanogram of DNA and/or in the ratio of taxa-specific to the total bacterial 16S rRNA gene copies (i.e., the relative abundance of the different taxa within the bacterial community) between the established soils or the glacier surface and the 10- to 160-year-old successional stages. A significantly higher proportion of Βetaproteobacteria (20%) was observed on the surface of both glaciers. However, no differences were observed between the 10- to 160-year-old successional stages in the number of taxa-specific 16S rRNA gene copies per nanogram of DNA or in the ratio of taxa-specific to the total bacterial 16S rRNA gene copies for the different taxa. Nevertheless, when the relative abundance data from all the studied taxa were combined and analyzed altogether, most of the sites could be distinguished from one other. This indicates that the overall composition of the bacterial community was more affected than the abundance of the targeted taxa by changes in environmental conditions along the chronosequences.     

Incorporating 16S Gene Copy Number Information Improves Estimates of Microbial Diversity and Abundance

The abundance of different SSU rRNA (“16S”) gene sequences in environmental samples is widely used in studies of microbial ecology as a measure of microbial community structure and diversity. However, the genomic copy number of the 16S gene varies greatly – from one in many species to up to 15 in some bacteria and to hundreds in some microbial eukaryotes. As a result of this variation the relative abundance of 16S genes in environmental samples can be attributed both to variation in the relative abundance of different organisms, and to variation in genomic 16S copy number among those organisms. Despite this fact, many studies assume that the abundance of 16S gene sequences is a surrogate measure of the relative abundance of the organisms containing those sequences. Here we present a method that uses data on sequences and genomic copy number of 16S genes along with phylogenetic placement and ancestral state estimation to estimate organismal abundances from environmental DNA sequence data. We use theory and simulations to demonstrate that 16S genomic copy number can be accurately estimated from the short reads typically obtained from high-throughput environmental sequencing of the 16S gene, and that organismal abundances in microbial communities are more strongly correlated with estimated abundances obtained from our method than with gene abundances. We re-analyze several published empirical data sets and demonstrate that the use of gene abundance versus estimated organismal abundance can lead to different inferences about community diversity and structure and the identity of the dominant taxa in microbial communities. Our approach will allow microbial ecologists to make more accurate inferences about microbial diversity and abundance based on 16S sequence data.     

Bacterial diversity in biofilms formed on condenser tube surfaces in a nuclear power plant

To elucidate the bacterial diversity in biofilms formed on a condenser tube from a nuclear power plant, 16S rRNA gene sequences were examined using a PCR-cloning-sequencing approach. Twelve operational taxonomic units were retrieved in the clone library, and the estimated species richness was low (13.2). Most of the clones (94.7%) were affiliated with α-Proteobacteria; Planctomycetes and γ-Proteobacteria were much rarer. Interestingly, except for one clone belonging to Pseudoalteromonas, most of the sequences displayed sequence similarities <97% of those of the closest type strains. Based on 16S rRNA phylogenetic analysis, most bacteria were assigned to novel taxa above the species level. The low species richness and unusual bacterial composition may be attributable to selective pressure from chlorine in the cooling water. To prevent or control bacterial biofilms in cooling circuits, additional studies of the physiology and ecology of these species will be essential.     

Bacterial diversity in rock varnish of extreme arid region of Turpan Basin

The bacterial community composition and diversity in rock varnish of Turpan Basin were investigated by restriction fragment length polymorphism (RFLP) and clone library of the 16S rRNA gene. 114 positive clones were screened, which could be grouped into 28 phylotypes and then further divided into 23 different operational taxonomic units (OTUs). These were affiliated into 5 phyla (Acidobacteria, Proteobacteria, Chloroflexi, Firmicutes and Cyanobacteria). Clones from actinobacteria were the dominant, accounting for 67.5% of total clones in the library, followed by Proteobacteria (15.8%), Chloroflexi (13.2%), Firmicutes (2.6%) and Cyanobacteria (0.9%). Rubrobacter (accounts for 35%) in the phylum Actinobacteria was the dominant genus and contained many species which might be resistant to gamma radiation. A 70% of the library clone sequences showed less 97% similarity to 16S rRNA gene sequences of standard strains obtained by pure culture. Shannon–Wiener index value of this study is 2.52 and is lower than deep-sea sediments, soils, lakes and other environments. Results of this study showed that bacterial diversity in rock varnishes of Turpan Basin was low, but maybe exist a large number of new unknown taxons, especially species that could well adapted to drought and resist radiation.     

Bacterial Microbiota Profiling in Gastritis without Helicobacter pylori Infection or Non-Steroidal Anti-Inflammatory Drug Use

Recent 16S ribosomal RNA gene (rRNA) molecular profiling of the stomach mucosa revealed a surprising complexity of microbiota. Helicobacter pylori infection and non-steroidal anti-inflammatory drug (NSAID) use are two main contributors to gastritis and peptic ulcer. However, little is known about the association between other members of the stomach microbiota and gastric diseases. In this study, cloning and sequencing of the 16S rRNA was used to profile the stomach microbiota from normal and gastritis patients. One hundred and thirty three phylotypes from eight bacterial phyla were identified. The stomach microbiota was found to be closely adhered to the mucosa. Eleven Streptococcus phylotypes were successfully cultivated from the biopsies. One to two genera represented a majority of clones within any of the identified phyla. We further developed two real-time quantitative PCR assays to quantify the relative abundance of the Firmicutes phylum and the Streptococcus genus. Significantly higher abundance of the Firmicutes phylum and the Streptococcus genus within the Firmicutes phylum was observed in patients with antral gastritis, compared with normal controls. This study suggests that the genus taxon level can largely represent much higher taxa such as the phylum. The clinical relevance and the mechanism underlying the altered microbiota composition in gastritis require further functional studies.     

Analysis of 23S rRNA genes in metagenomes - a case study from the Global Ocean Sampling Expedition

As an evolutionary marker, 23S ribosomal RNA (rRNA) offers more diagnostic sequence stretches and greater sequence variation than 16S rRNA. However, 23S rRNA is still not as widely used. Based on 80 metagenome samples from the Global Ocean Sampling (GOS) Expedition, the usefulness and taxonomic resolution of 23S rRNA were compared to those of 16S rRNA. Since 23S rRNA is approximately twice as large as 16S rRNA, twice as many 23S rRNA gene fragments were retrieved from the GOS reads than 16S rRNA gene fragments, with 23S rRNA gene fragments being generally about 100 bp longer. Datasets for 16S and 23S rRNA sequences revealed similar relative abundances for major marine bacterial and archaeal taxa. However, 16S rRNA sequences had a better taxonomic resolution due to their significantly larger reference database.Reevaluation of the specificity of previously published PCR amplification primers and group specific fluorescence in situ hybridization probes on this metagenomic set of non-amplified 23S rRNA sequences revealed that out of 16 primers investigated, only two had more than 90% target group coverage. Evaluations of two probes, BET42a and GAM42a, were in accordance with previous evaluations, with a discrepancy in the target group coverage of the GAM42a probe when evaluated against the GOS metagenomic dataset.     

Changes in the microbial community structure of bacteria,archaea and fungi in response to elevated CO2 and warming in an Australian native grassland soil

The microbial community structure of bacteria, archaea and fungi is described in an Australian native grassland soil after more than 5 years exposure to different atmospheric CO₂ concentrations ([CO₂]) (ambient, + 550 ppm) and temperatures (ambient, + 2°C) under different plant functional types (C ₃ and C ₄ grasses) and at two soil depths (0–5 cm and 5–10 cm). Archaeal community diversity was influenced by elevated [CO₂], while under warming archaeal 16S rRNA gene copy numbers increased for C ₄ plant Themeda triandra and decreased for the C ₃ plant community (P < 0.05). Fungal community diversity resulted in three groups based upon elevated [CO₂], elevated [CO₂] plus warming and ambient [CO₂]. Overall bacterial community diversity was influenced primarily by depth. Specific bacterial taxa changed in richness and relative abundance in response to climate change factors when assessed by a high‐resolution 16S rRNA microarray (PhyloChip). Operational taxonomic unit signal intensities increased under elevated [CO₂] for both Firmicutes and Bacteroidetes, and increased under warming for Actinobacteria and Alphaproteobacteria. For the interaction of elevated [CO₂] and warming there were 103 significant operational taxonomic units (P < 0.01) representing 15 phyla and 30 classes. The majority of these operational taxonomic units increased in abundance for elevated [CO₂] plus warming plots, while abundance declined in warmed or elevated [CO₂] plots. Bacterial abundance (16S rRNA gene copy number) was significantly different for the interaction of elevated [CO₂] and depth (P < 0.05) with decreased abundance under elevated [CO₂] at 5–10 cm, and for Firmicutes under elevated [CO₂] (P < 0.05). Bacteria, archaea and fungi in soil responded differently to elevated [CO₂], warming and their interaction. Taxa identified as significantly climate‐responsive could show differing trends in the direction of response (‘+’ or ‘?’) under elevated CO₂ or warming, which could then not be used to predict their interactive effects supporting the need to investigate interactive effects for climate change. The approach of focusing on specific taxonomic groups provides greater potential for understanding complex microbial community changes in ecosystems under climate change.     

Bacterial diversity in biofilms formed on condenser tube surfaces in a nuclear power plant

To elucidate the bacterial diversity in biofilms formed on a condenser tube from a nuclear power plant, 16S rRNA gene sequences were examined using a PCR-cloning-sequencing approach. Twelve operational taxonomic units were retrieved in the clone library, and the estimated species richness was low (13.2). Most of the clones (94.7%) were affiliated with α-Proteobacteria; Planctomycetes and γ-Proteobacteria were much rarer. Interestingly, except for one clone belonging to Pseudoalteromonas, most of the sequences displayed sequence similarities <97% of those of the closest type strains. Based on 16S rRNA phylogenetic analysis, most bacteria were assigned to novel taxa above the species level. The low species richness and unusual bacterial composition may be attributable to selective pressure from chlorine in the cooling water. To prevent or control bacterial biofilms in cooling circuits, additional studies of the physiology and ecology of these species will be essential.     

Archaeal abundance across a pH gradient in an arable soil and its relationship to bacterial and fungal growth rates

Soil pH is one of the most influential factors for the composition of bacterial and fungal communities, but the influence of soil pH on the distribution and composition of soil archaeal communities has yet to be systematically addressed. The primary aim of this study was to determine how total archaeal abundance (quantitative PCR [qPCR]-based estimates of 16S rRNA gene copy numbers) is related to soil pH across a pH gradient (pH 4.0 to 8.3). Secondarily, we wanted to assess how archaeal abundance related to bacterial and fungal growth rates across the same pH gradient. We identified two distinct and opposite effects of pH on the archaeal abundance. In the lowest pH range (pH 4.0 to 4.7), the abundance of archaea did not seem to correspond to pH. Above this pH range, there was a sharp, almost 4-fold decrease in archaeal abundance, reaching a minimum at pH 5.1 to 5.2. The low abundance of archaeal 16S rRNA gene copy numbers at this pH range then sharply increased almost 150-fold with pH, resulting in an increase in the ratio between archaeal and bacterial copy numbers from a minimum of 0.002 to more than 0.07 at pH 8. The nonuniform archaeal response to pH could reflect variation in the archaeal community composition along the gradient, with some archaea adapted to acidic conditions and others to neutral to slightly alkaline conditions. This suggestion is reinforced by observations of contrasting outcomes of the (competitive) interactions between archaea, bacteria, and fungi toward the lower and higher ends of the examined pH gradient.     

Effect of the Earthworms Lumbricus terrestris and Aporrectodea caliginosa on Bacterial Diversity in Soil

Earthworms ingest large amounts of soil and have the potential to radically alter the biomass, activity, and structure of the soil microbial community. In this study, the diversity of eight bacterial groups from fresh soil, gut, and casts of the earthworms Lumbricus terrestris and Aporrectodea caliginosa were studied by single-strand conformation polymorphism (SSCP) analysis using both newly designed 16S rRNA gene-specific primer sets targeting Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria, Deltaproteobacteria, Bacteroidetes, Verrucomicrobia, Planctomycetes, and Firmicutes and a conventional universal primer set for SSCP, with RNA and DNA as templates. In parallel, the study of the relative abundance of these taxonomic groups in the same samples was performed using fluorescence in situ hybridization. Bacteroidetes, Alphaproteobacteria, and Betaproteobacteria were predominant in communities from the soil and worm cast samples. Representatives of classes Flavobacteria and Sphingobacteria (Bacteroidetes) and Pseudomonas spp. (low-abundant Gammaproteobacteria) were detected in soil and worm cast samples with conventional and taxon-targeting SSCP and through the sequence analysis of 16S rRNA clone libraries. Physiologically active unclassified Sphingomonadaceae (Alphaproteobacteria) and Alcaligenes spp. (Betaproteobacteria) also maintained their diversities during transit through the earthworm intestine and were found on taxon-targeting SSCP profiles from the soil and worm cast samples. In conclusion, our results suggest that some specific bacterial taxonomic groups maintain their diversity and even increase their relative numbers during transit through the gastrointestinal tract of earthworms.     

Shotgun metagenomics and microscopy indicate diverse cyanophytes,other bacteria,and microeukaryotes in the epimicrobiota of a northern Chilean wetland Nostoc (Cyanobacteria)

Prokaryotic Nostoc, one of the world's most conspicuous and widespread algal genera (similar to eukaryotic algae, plants, and animals) is known to support a microbiome that influences host ecological roles. Past taxonomic characterizations of surface microbiota (epimicrobiota) of free‐living Nostoc sampled from freshwater systems employed 16S rRNA genes, typically amplicons. We compared taxa identified from 16S, 18S, 23S, and 28S rRNA gene sequences filtered from shotgun metagenomic sequence and used microscopy to illuminate epimicrobiota diversity for Nostoc sampled from a wetland in the northern Chilean Altiplano. Phylogenetic analysis and rRNA gene sequence abundance estimates indicated that the host was related to Nostoc punctiforme PCC 73102. Epimicrobiota were inferred to include 18 epicyanobacterial genera or uncultured taxa, six epieukaryotic algal genera, and 66 anoxygenic bacterial genera, all having average genomic coverage ≥90X. The epicyanobacteria Geitlerinemia, Oscillatoria, Phormidium, and an uncultured taxon were detected only by 16S rRNA gene; Gloeobacter and Pseudanabaena were detected using 16S and 23S; and Phormididesmis, Neosynechococcus, Symphothece, Aphanizomenon, Nodularia, Spirulina, Nodosilinea, Synechococcus, Cyanobium, and Anabaena (the latter corroborated by microscopy), plus two uncultured cyanobacterial taxa (JSC12, O77) were detected only by 23S rRNA gene sequences. Three chlamydomonad and two heterotrophic stramenopiles genera were inferred from 18S; the streptophyte green alga Chaetosphaeridium globosum was detected by microscopy and 28S rRNA genes, but not 18S rRNA genes. Overall, >60% of epimicrobial taxa were detected by markers other than 16S rRNA genes. Some algal taxa observed microscopically were not detected from sequence data. Results indicate that multiple taxonomic markers derived from metagenomic sequence data and microscopy increase epimicrobiota detection.     

Taxonomic Classification of Bacterial 16S rRNA Genes Using Short Sequencing Reads: Evaluation of Effective Study Designs

Massively parallel high throughput sequencing technologies allow us to interrogate the microbial composition of biological samples at unprecedented resolution. The typical approach is to perform high-throughout sequencing of 16S rRNA genes, which are then taxonomically classified based on similarity to known sequences in existing databases. Current technologies cause a predicament though, because although they enable deep coverage of samples, they are limited in the length of sequence they can produce. As a result, high-throughout studies of microbial communities often do not sequence the entire 16S rRNA gene. The challenge is to obtain reliable representation of bacterial communities through taxonomic classification of short 16S rRNA gene sequences. In this study we explored properties of different study designs and developed specific recommendations for effective use of short-read sequencing technologies for the purpose of interrogating bacterial communities, with a focus on classification using naïve Bayesian classifiers. To assess precision and coverage of each design, we used a collection of ∼8,500 manually curated 16S rRNA gene sequences from cultured bacteria and a set of over one million bacterial 16S rRNA gene sequences retrieved from environmental samples, respectively. We also tested different configurations of taxonomic classification approaches using short read sequencing data, and provide recommendations for optimal choice of the relevant parameters. We conclude that with a judicious selection of the sequenced region and the corresponding choice of a suitable training set for taxonomic classification, it is possible to explore bacterial communities at great depth using current technologies, with only a minimal loss of taxonomic resolution.     

Multiple copies of 16S rRNA gene affect the restriction patterns and DGGE profile revealed by analysis of genome database

The use of 16S rRNA gene has been a “golden” method to determine the diversity of microbial communities in environmental samples, phylogenetic relationships of prokaryotes and taxonomic position of newly isolated organisms. However due to the presence of multiple heterogeneous 16S rRNA gene copies in many strains, the interpretation of microbial ecology via 16S rRNA sequences is complicated. Purpose of present paper is to demonstrate the extent to which the multiple heterogeneous 16S rRNA gene copies affect RFLP patterns and DGGE profiles by using the genome database. In present genome database, there are 782 bacterial strains in total whose genomes have been completely sequenced and annotated. Among the total strains, 639 strains (82%) possess multiple 16S rRNA gene copies, 415 strains (53%) whose multiple copies are heterogeneous in sequences as revealed by alignment, 236 strains (30%) whose multiple copies show different restrict patterns by CSP6I+HinfI, MspI+RsaI or HhaI as analyzed in silico. Polymorphisms of the multiple copies in certain strains were further characterized by G+C% and phy-logentic distances based on the sequences of V3 region, which are linked to DGGE patters. Polymorphisms of a few strains were shown as examples. Using artificial communities, it is demonstrated that the presence of multiple heterogeneous 16S rRNA gene copies potentially leads to over-estimation of the diversity of a community. It is suggested that care must be taken when interpreting 16S rRNA-based RFLP and DGGE data and profiling an environmental community.