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.     

Miniprimer PCR, a New Lens for Viewing the Microbial World

Molecular methods based on the 16S rRNA gene sequence are used widely in microbial ecology to reveal the diversity of microbial populations in environmental samples. Here we show that a new PCR method using an engineered polymerase and 10-nucleotide “miniprimers” expands the scope of detectable sequences beyond those detected by standard methods using longer primers and Taq polymerase. After testing the method in silico to identify divergent ribosomal genes in previously cloned environmental sequences, we applied the method to soil and microbial mat samples, which revealed novel 16S rRNA gene sequences that would not have been detected with standard primers. Deeply divergent sequences were discovered with high frequency and included representatives that define two new division-level taxa, designated CR1 and CR2, suggesting that miniprimer PCR may reveal new dimensions of microbial diversity.     

Assembly-Driven Community Genomics of a Hypersaline Microbial Ecosystem

Microbial populations inhabiting a natural hypersaline lake ecosystem in Lake Tyrrell, Victoria, Australia, have been characterized using deep metagenomic sampling, iterative de novo assembly, and multidimensional phylogenetic binning. Composite genomes representing habitat-specific microbial populations were reconstructed for eleven different archaea and one bacterium, comprising between 0.6 and 14.1% of the planktonic community. Eight of the eleven archaeal genomes were from microbial species without previously cultured representatives. These new genomes provide habitat-specific reference sequences enabling detailed, lineage-specific compartmentalization of predicted functional capabilities and cellular properties associated with both dominant and less abundant community members, including organisms previously known only by their 16S rRNA sequences. Together, these data provide a comprehensive, culture-independent genomic blueprint for ecosystem-wide analysis of protein functions, population structure, and lifestyles of co-existing, co-evolving microbial groups within the same natural habitat. The “assembly-driven” community genomic approach demonstrated in this study advances our ability to push beyond single gene investigations, and promotes genome-scale reconstructions as a tangible goal in the quest to define the metabolic, ecological, and evolutionary dynamics that underpin environmental microbial diversity.     

16S rRNA Gene-Based Oligonucleotide Microarray for Environmental Monitoring of the Betaproteobacterial Order “Rhodocyclales”

For simultaneous identification of members of the betaproteobacterial order “Rhodocyclales” in environmental samples, a 16S rRNA gene-targeted oligonucleotide microarray (RHC-PhyloChip) consisting of 79 probes was developed. Probe design was based on phylogenetic analysis of available 16S rRNA sequences from all cultured and as yet uncultured members of the “Rhodocyclales.” The multiple nested probe set was evaluated for microarray hybridization with 16S rRNA gene PCR amplicons from 29 reference organisms. Subsequently, the RHC-PhyloChip was successfully used for cultivation-independent “Rhodocyclales” diversity analysis in activated sludge from an industrial wastewater treatment plant. The implementation of a newly designed “Rhodocyclales”-selective PCR amplification system prior to microarray hybridization greatly enhanced the sensitivity of the RHC-PhyloChip and thus enabled the detection of “Rhodocyclales” populations with relative abundances of less than 1% of all bacteria (as determined by fluorescence in situ hybridization) in the activated sludge. The presence of as yet uncultured Zoogloea-, Ferribacterium/Dechloromonas-, and Sterolibacterium-related bacteria in the industrial activated sludge, as indicated by the RHC-PhyloChip analysis, was confirmed by retrieval of their 16S rRNA gene sequences and subsequent phylogenetic analysis, demonstrating the suitability of the RHC-PhyloChip as a novel monitoring tool for environmental microbiology.     

Frequency of formation of chimeric molecules as a consequence of PCR coamplification of 16S rRNA genes from mixed bacterial genomes.

PCR is routinely used in amplification and cloning of rRNA genes from environmental DNA samples for studies of microbial community structure and identification of novel organisms. There have been concerns about generation of chimeric sequences as a consequence of PCR coamplification of highly conserved genes, because such sequences may lead to reports of nonexistent organisms. To quantify the frequency of chimeric molecule formation, mixed genomic DNAs from eight actinomycete species whose 16S rRNA sequences had been determined were used for PCR coamplification of 16S rRNA genes. A large number of cloned 16S ribosomal DNAs were examined by sequence analysis, and chimeric molecules were identified by multiple-sequence alignment with reference species. Here, we report that the level of occurrence of chimeric sequences after 30 cycles of PCR amplification was 32%. We also show that PCR-induced chimeras were formed between different rRNA gene copies from the same organism. Because of the wide use of PCR for direct isolation of 16S rRNA sequences from environmental DNA to assess microbial diversity, the extent of chimeric molecule formation deserves serious attention.     

Short-Read Assembly of Full-Length 16S Amplicons Reveals Bacterial Diversity in Subsurface Sediments

In microbial ecology, a fundamental question relates to how community diversity and composition change in response to perturbation. Most studies have had limited ability to deeply sample community structure (e.g. Sanger-sequenced 16S rRNA libraries), or have had limited taxonomic resolution (e.g. studies based on 16S rRNA hypervariable region sequencing). Here, we combine the higher taxonomic resolution of near-full-length 16S rRNA gene amplicons with the economics and sensitivity of short-read sequencing to assay the abundance and identity of organisms that represent as little as 0.01% of sediment bacterial communities. We used a new version of EMIRGE optimized for large data size to reconstruct near-full-length 16S rRNA genes from amplicons sheared and sequenced with Illumina technology. The approach allowed us to differentiate the community composition among samples acquired before perturbation, after acetate amendment shifted the predominant metabolism to iron reduction, and once sulfate reduction began. Results were highly reproducible across technical replicates, and identified specific taxa that responded to the perturbation. All samples contain very high alpha diversity and abundant organisms from phyla without cultivated representatives. Surprisingly, at the time points measured, there was no strong loss of evenness, despite the selective pressure of acetate amendment and change in the terminal electron accepting process. However, community membership was altered significantly. The method allows for sensitive, accurate profiling of the “long tail” of low abundance organisms that exist in many microbial communities, and can resolve population dynamics in response to environmental change.     

Development and Application of a Real-Time PCR Approach for Quantification of Uncultured Bacteria in the Central Baltic Sea

We have developed a highly sensitive approach to assess the abundance of uncultured bacteria in water samples from the central Baltic Sea by using a noncultured member of the “Epsilonproteobacteria” related to Thiomicrospira denitrificans as an example. Environmental seawater samples and samples enriched for the target taxon provided a unique opportunity to test the approach over a broad range of abundances. The approach is based on a combination of taxon- and domain-specific real-time PCR measurements determining the relative T. denitrificans-like 16S rRNA gene and 16S rRNA abundances, as well as the determination of total cell counts and environmental RNA content. It allowed quantification of T. denitrificans-like 16S rRNA molecules or 16S rRNA genes as well as calculation of the number of ribosomes per T. denitrificans-like cell. Every real-time measurement and its specific primer system were calibrated using environmental nucleic acids obtained from the original habitat for external standardization. These standards, as well as the respective samples to be measured, were prepared from the same DNA or RNA extract. Enrichment samples could be analyzed directly, whereas environmental templates had to be preamplified with general bacterial primers before quantification. Preamplification increased the sensitivity of the assay by more than 4 orders of magnitude. Quantification of enrichments with or without a preamplification step yielded comparable results. T. denitrificans-like 16S rRNA molecules ranged from 7.1 × 10³ to 4.4 × 10⁹ copies ml⁻¹ or 0.002 to 49.7% relative abundance. T. denitrificans-like 16S rRNA genes ranged from 9.0 × 10¹ to 2.2 ×10⁶ copies ml⁻¹ or 0.01 to 49.7% relative abundance. Detection limits of this real-time-PCR approach were 20 16S rRNA molecules or 0.2 16S rRNA gene ml⁻¹. The number of ribosomes per T. denitrificans-like cell was estimated to range from 20 to 200 in seawater and reached up to 2,000 in the enrichments. The results indicate that our real-time PCR approach can be used to determine cellular and relative abundances of uncultured marine bacterial taxa and to provide information about their levels of activity in their natural environment.     

Diversity and Distribution in Hypersaline Microbial Mats of Bacteria Related to Chloroflexus spp.

Filamentous bacteria containing bacteriochlorophylls c and a were enriched from hypersaline microbial mats. Based on phylogenetic analyses of 16S rRNA gene sequences, these organisms form a previously undescribed lineage distantly related to Chloroflexus spp. We developed and tested a set of PCR primers for the specific amplification of 16S rRNA genes from filamentous phototrophic bacteria within the kingdom of “green nonsulfur bacteria.” PCR products recovered from microbial mats in a saltern in Guerrero Negro, Mexico, were subjected to cloning or denaturing gradient gel electrophoresis and then sequenced. We found evidence of a high diversity of bacteria related to Chloroflexus which exhibit different distributions along a gradient of salinity from 5.5 to 16%.     

Use of 16S rRNA and rpoB Genes as Molecular Markers for Microbial Ecology Studies

Several characteristics of the 16S rRNA gene, such as its essential function, ubiquity, and evolutionary properties, have allowed it to become the most commonly used molecular marker in microbial ecology. However, one fact that has been overlooked is that multiple copies of this gene are often present in a given bacterium. These intragenomic copies can differ in sequence, leading to identification of multiple ribotypes for a single organism. To evaluate the impact of such intragenomic heterogeneity on the performance of the 16S rRNA gene as a molecular marker, we compared its phylogenetic and evolutionary characteristics to those of the single-copy gene rpoB. Full-length gene sequences and gene fragments commonly used for denaturing gradient gel electrophoresis were compared at various taxonomic levels. Heterogeneity found between intragenomic 16S rRNA gene copies was concentrated in specific regions of rRNA secondary structure. Such “heterogeneity hot spots” occurred within all gene fragments commonly used in molecular microbial ecology. This intragenomic heterogeneity influenced 16S rRNA gene tree topology, phylogenetic resolution, and operational taxonomic unit estimates at the species level or below. rpoB provided comparable phylogenetic resolution to that of the 16S rRNA gene at all taxonomic levels, except between closely related organisms (species and subspecies levels), for which it provided better resolution. This is particularly relevant in the context of a growing number of studies focusing on subspecies diversity, in which single-copy protein-encoding genes such as rpoB could complement the information provided by the 16S rRNA gene.     

Generation of multimillion-sequence 16S rRNA gene libraries from complex microbial communities by assembling paired-end illumina reads

Microbial communities host unparalleled taxonomic diversity. Adequate characterization of environmental and host-associated samples remains a challenge for microbiologists, despite the advent of 16S rRNA gene sequencing. In order to increase the depth of sampling for diverse bacterial communities, we developed a method for sequencing and assembling millions of paired-end reads from the 16S rRNA gene (spanning the V3 region; ～200 nucleotides) by using an Illumina genome analyzer. To confirm reproducibility and to identify a suitable computational pipeline for data analysis, sequence libraries were prepared in duplicate for both a defined mixture of DNAs from known cultured bacterial isolates (>1 million postassembly sequences) and an Arctic tundra soil sample (>6 million postassembly sequences). The Illumina 16S rRNA gene libraries represent a substantial increase in number of sequences over all extant next-generation sequencing approaches (e.g., 454 pyrosequencing), while the assembly of paired-end 125-base reads offers a methodological advantage by incorporating an initial quality control step for each 16S rRNA gene sequence. This method incorporates indexed primers to enable the characterization of multiple microbial communities in a single flow cell lane, may be modified readily to target other variable regions or genes, and demonstrates unprecedented and economical access to DNAs from organisms that exist at low relative abundances.     

沙颍河下游城市黑臭内河沉积物微生物群落季节变化特征

为了研究沙颍河下游城市黑臭内河不同季节沉积物微生物群落特征,对安徽省阜阳市黑臭内河中清河、七渔河表层沉积物进行16S rDNA高通量测序。结果发现:黑臭河流中沉积物的微生物多样性指数均不高,但是表现出一定的变化规律,即春季>冬季≥夏季>秋季;通过冗余性分析发现微生物多样性受季节与沉积物pH影响较显著。分析沉积物门水平上的微生物群落结构发现,季节、温度、TN及SOM对微生物影响较大。变形杆菌、厚壁菌门、绿弯菌门、疣微菌门、拟杆菌门和放线菌门等优势菌门的相对丰度在季节水平上存在差异,春季厚壁菌门、绿弯菌门、放线菌门和酸杆菌门相对丰度较高,其中绿弯菌门和酸杆菌门是已知指示污染的微生物,变形杆菌门相对较少。秋季疣微菌门与拟杆菌门相对丰度显著减小,变形杆菌门相对其他季节显著增加。样品中共发现16个硫酸盐还原菌(SRB)菌属,其中Desulfoprunum是丰度最高的菌属。春季沉积物中SRB的类群最多,相对丰度最大;硫酸盐还原菌群与SO42-、TN、SOM、Cl–等呈显著正相关。上述结果为营养盐控制时机的选择从而有效避免河流中黑臭物质的产生提供了一定参考。     

Threatened Corals Provide Underexplored Microbial Habitats

Contemporary in-depth sequencing of environmental samples has provided novel insights into microbial community structures, revealing that their diversity had been previously underestimated. Communities in marine environments are commonly composed of a few dominant taxa and a high number of taxonomically diverse, low-abundance organisms. However, studying the roles and genomic information of these “rare” organisms remains challenging, because little is known about their ecological niches and the environmental conditions to which they respond. Given the current threat to coral reef ecosystems, we investigated the potential of corals to provide highly specialized habitats for bacterial taxa including those that are rarely detected or absent in surrounding reef waters. The analysis of more than 350,000 small subunit ribosomal RNA (16S rRNA) sequence tags and almost 2,000 nearly full-length 16S rRNA gene sequences revealed that rare seawater biosphere members are highly abundant or even dominant in diverse Caribbean corals. Closely related corals (in the same genus/family) harbored similar bacterial communities. At higher taxonomic levels, however, the similarities of these communities did not correlate with the phylogenetic relationships among corals, opening novel questions about the evolutionary stability of coral-microbial associations. Large proportions of OTUs (28.7–49.1%) were unique to the coral species of origin. Analysis of the most dominant ribotypes suggests that many uncovered bacterial taxa exist in coral habitats and await future exploration. Our results indicate that coral species, and by extension other animal hosts, act as specialized habitats of otherwise rare microbes in marine ecosystems. Here, deep sequencing provided insights into coral microbiota at an unparalleled resolution and revealed that corals harbor many bacterial taxa previously not known. Given that two of the coral species investigated are listed as threatened under the U.S. Endangered Species Act, our results add an important microbial diversity-based perspective to the significance of conserving coral reefs.     

Comparative Analysis of DNA Word Abundances in Four Yeast Genomes Using a Novel Statistical Background Model

Previous studies have shown that the identification and analysis of both abundant and rare k-mers or “DNA words of length k” in genomic sequences using suitable statistical background models can reveal biologically significant sequence elements. Other studies have investigated the uni/multimodal distribution of k-mer abundances or “k-mer spectra” in different DNA sequences. However, the existing background models are affected to varying extents by compositional bias. Moreover, the distribution of k-mer abundances in the context of related genomes has not been studied previously. Here, we present a novel statistical background model for calculating k-mer enrichment in DNA sequences based on the average of the frequencies of the two (k-1) mers for each k-mer. Comparison of our null model with the commonly used ones, including Markov models of different orders and the single mismatch model, shows that our method is more robust to compositional AT-rich bias and detects many additional, repeat-poor over-abundant k-mers that are biologically meaningful. Analysis of overrepresented genomic k-mers (4≤k≤16) from four yeast species using this model showed that the fraction of overrepresented DNA words falls linearly as k increases; however, a significant number of overabundant k-mers exists at higher values of k. Finally, comparative analysis of k-mer abundance scores across four yeast species revealed a mixture of unimodal and multimodal spectra for the various genomic sub-regions analyzed.     

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.     

Digital Genotyping of Macrosatellites and Multicopy Genes Reveals Novel Biological Functions Associated with Copy Number Variation of Large Tandem Repeats

Tandem repeats are common in eukaryotic genomes, but due to difficulties in assaying them remain poorly studied. Here, we demonstrate the utility of Nanostring technology as a targeted approach to perform accurate measurement of tandem repeats even at extremely high copy number, and apply this technology to genotype 165 HapMap samples from three different populations and five species of non-human primates. We observed extreme variability in copy number of tandemly repeated genes, with many loci showing 5–10 fold variation in copy number among humans. Many of these loci show hallmarks of genome assembly errors, and the true copy number of many large tandem repeats is significantly under-represented even in the high quality ‘finished’ human reference assembly. Importantly, we demonstrate that most large tandem repeat variations are not tagged by nearby SNPs, and are therefore essentially invisible to SNP-based GWAS approaches. Using association analysis we identify many cis correlations of large tandem repeat variants with nearby gene expression and DNA methylation levels, indicating that variations of tandem repeat length are associated with functional effects on the local genomic environment. This includes an example where expansion of a macrosatellite repeat is associated with increased DNA methylation and suppression of nearby gene expression, suggesting a mechanism termed “repeat induced gene silencing”, which has previously been observed only in transgenic organisms. We also observed multiple signatures consistent with altered selective pressures at tandemly repeated loci, suggesting important biological functions. Our studies show that tandemly repeated loci represent a highly variable fraction of the genome that have been systematically ignored by most previous studies, copy number variation of which can exert functionally significant effects. We suggest that future studies of tandem repeat loci will lead to many novel insights into their role in modulating both genomic and phenotypic diversity.     

Microbial Diversity in Hummock and Hollow Soils of Three Wetlands on the Qinghai-Tibetan Plateau Revealed by 16S rRNA Pyrosequencing

The wetlands of the Qinghai-Tibetan Plateau are believed to play an important role in global nutrient cycling, but the composition and diversity of microorganisms in this ecosystem are poorly characterized. An understanding of the effects of geography and microtopography on microbial populations will provide clues to the underlying mechanisms that structure microbial communities. In this study, we used pyrosequencing-based analysis of 16S rRNA gene sequences to assess and compare the composition of soil microbial communities present in hummock and hollow soils from three wetlands (Dangxiong, Hongyuan and Maduo) on the Qinghai-Tibetan Plateau, the world’s highest plateau. A total of 36 bacterial phyla were detected. Proteobacteria (34.5% average relative abundance), Actinobacteria (17.3%) and Bacteroidetes (11%) had the highest relative abundances across all sites. Chloroflexi, Acidobacteria, Verrucomicrobia, Firmicutes, and Planctomycetes were also relatively abundant (1–10%). In addition, archaeal sequences belonging to Euryarchaea, Crenarchaea and Thaumarchaea were detected. Alphaproteobacteria sequences, especially of the order Rhodospirillales, were significantly more abundant in Maduo than Hongyuan and Dangxiong wetlands. Compared with Hongyuan soils, Dangxiong and Maduo had significantly higher relative abundances of Gammaproteobacteria sequences (mainly order Xanthomonadales). Hongyuan wetland had a relatively high abundance of methanogens (mainly genera Methanobacterium, Methanosarcina and Methanosaeta) and methanotrophs (mainly Methylocystis) compared with the other two wetlands. Principal coordinate analysis (PCoA) indicated that the microbial community structure differed between locations and microtopographies and canonical correspondence analysis indicated an association between microbial community structure and soil properties or geography. These insights into the microbial community structure and the main controlling factors in wetlands of the Qinghai-Tibetan Plateau provide a valuable background for further studies on biogeochemical processes in this distinct ecosystem.     

Dynamics of Bacterial and Fungal Communities on Decaying Salt Marsh Grass

Both bacteria and fungi play critical roles in decomposition processes in many natural environments, yet only rarely have they been studied as an integrated microbial community. Here we describe the bacterial and fungal assemblages associated with two decomposition stages of Spartina alterniflora detritus in a productive southeastern U.S. salt marsh. 16S rRNA genes and 18S-to-28S internal transcribed spacer (ITS) regions were used to target the bacterial and ascomycete fungal communities, respectively, based on DNA sequence analysis of isolates and environmental clones and by using community fingerprinting based on terminal restriction fragment length polymorphism (T-RFLP) analysis. Seven major bacterial taxa (six affiliated with the α-Proteobacteria and one with the Cytophagales) and four major fungal taxa were identified over five sample dates spanning 13 months. Fungal terminal restriction fragments (T-RFs) were informative at the species level; however, bacterial T-RFs frequently comprised a number of related genera. Amplicon abundances indicated that the salt marsh saprophyte communities have little-to-moderate variability spatially or with decomposition stage, but considerable variability temporally. However, the temporal variability could not be readily explained by either successional shifts or simple relationships with environmental factors. Significant correlations in abundance (both positive and negative) were found among dominant fungal and bacterial taxa that possibly indicate ecological interactions between decomposer organisms. Most associations involved one of four microbial taxa: two groups of bacteria affiliated with the α-Proteobacteria and two ascomycete fungi (Phaeosphaeria spartinicola and environmental isolate “4clt”).