Confident phylogenetic identification of uncultured prokaryotes through long read amplicon sequencing of the 16S-ITS-23S rRNA operon

Amplicon sequencing of the 16S rRNA gene is the predominant method to quantify microbial compositions and to discover novel lineages. However, traditional short amplicons often do not contain enough information to confidently resolve their phylogeny. Here we present a cost-effective protocol that amplifies a large part of the rRNA operon and sequences the amplicons with PacBio technology. We tested our method on a mock community and developed a read-curation pipeline that reduces the overall read error rate to 0.18%. Applying our method on four environmental samples, we captured near full-length rRNA operon amplicons from a large diversity of prokaryotes. The method operated at moderately high-throughput (22286–37,850 raw ccs reads) and generated a large amount of putative novel archaeal 23S rRNA gene sequences compared to the archaeal SILVA database. These long amplicons allowed for higher resolution during taxonomic classification by means of long (∼1000 bp) 16S rRNA gene fragments and for substantially more confident phylogenies by means of combined near full-length 16S and 23S rRNA gene sequences, compared to shorter traditional amplicons (250 bp of the 16S rRNA gene). We recommend our method to those who wish to cost-effectively and confidently estimate the phylogenetic diversity of prokaryotes in environmental samples at high throughput.     

Identification of microorganisms by PCR amplification and sequencing of a universal amplified ribosomal region present in both prokaryotes and eukaryotes

The small ribosomal subunit contains 16S rRNA in prokaryotes and 18S rRNA in eukaryotes. Even though it has been known that some small ribosomal sequences are conserved in 16S rRNA and 18S rRNA molecules, they have been used separately for taxonomic and phylogenetic studies. Here, we report the existence of two highly conserved ribosomal sequences in all organisms that allow the amplification of a zone containing approximately 495 bp in prokaryotes and 508 bp in eukaryotes which we have named the "Universal Amplified Ribosomal Region" (UARR). Amplification and sequencing of this zone is possible using the same two universal primers (U1F and U1R) designed on the basis of two highly conserved ribosomal sequences. The UARR encompasses the V6, V7 and V8 domains from SSU rRNA in both prokaryotes and eukaryotes. The internal sequence of this zone in prokaryotes and eukaryotes is variable and the differences become less marked on descent from phyla to species. Nevertheless, UARR sequence allows species from the same genus to be differentiated. Thus, by UARR sequence analysis the construction of universal phylogenetic trees is possible, including species of prokaryotic and eukaryotic microorganisms together. Single isolates of prokaryotic and eukaryotic microorganisms from different sources can be identified by amplification and sequencing of UARR using the same pair of primers and the same PCR and sequencing conditions.     

Copy number of ribosomal operons in prokaryotes and its effect on phylogenic analyses

Different aspects of the presence of multiple copies of ribosomal operons in prokaryotic genomes are reviewed. Structure of prokaryotic ribosomal operons is briefly described. The available data are summarized regarding the copy number of ribosomal genes in various prokaryotic genomes, the degree of polymorphism of their individual copies, physiological and evolutionary aspects of the presence of the multiple copies of ribosomal genes. The review also considers the influence of the presence of multiple copies of ribosomal genes on the results of identification of prokaryotic isolates and of the studies of prokaryotic diversity in environmental samples based on phylogenetic analysis of 16S rRNA gene sequences.     

Copy Number of Ribosomal Operons in Prokaryotes and Its Effect on Phylogenetic Analyses

Different aspects of the presence of multiple copies of ribosomal operons in prokaryotic genomes are reviewed. The structure of prokaryotic ribosomal operons is briefly described. The available data are summarized regarding the copy number of ribosomal genes in various prokaryotic genomes, the degree of polymorphism of their individual copies, and physiological and evolutional aspects of the presence of the multiple copies of ribosomal genes. The review also considers the influence of the presence of multiple copies of ribosomal genes on the results of identification of prokaryotic isolates and of the studies of prokaryotic diversity in environmental samples based on phylogenetic analysis of 16S rRNA gene sequences.     

Nucleotide sequence of the single ribosomal RNA operon of pea chloroplast DNA

The nucleotide sequence of an 8 kbp region of pea ( Pisum sativum L.) chloroplast DNA containing the rRNA operon and putative promoter sites has been determined and compared to the corresponding sequences from maize, tobacco and the liverwort Marchantia polymorpha . The chloroplast DNA species of all vascular plants investigated, with the exception of a few legumes including pea, and of Marchantia contain an inverted repeat with an rRNA operon. The pea rRNA operon is the first sequenced rRNA operon from a plant with only one copy of the rRNA genes per molecule of chloroplast DNA. The organization of the operon is the same as for maize, tobacco and Marchantia . i.e. tRNA-Val gene/16S rRNA gene/spacer with intron-containing genes for tRNA-Ile and tRNA-Ala/23S rRNA gene/4.5S rRNA gene/5S rRNA gene. Current evidence suggests that the tRNA-Val gene may not be contranscribed with the other genes. For pea 16S, 23S, 4.5S and 5S rRNA have 1488, 2813, 105 and 121 nucleotides, respectively. The homologies of the entire operon (the tRNA-Val gene - 5S rRNA region) to those from tobacco, maize and Marchantia are 88, 82 and 79%, respectively. The corresponding homologies for tobacco/maize, tobacco/ Marchantia and maize/ Marchantia have similar values. The 16S and 23S rRNA genes from pea are more than 90% homologous to those from the 3 other species. We conclude that the fact that pea only has one set of rRNA genes per molecule of chloroplast DNA is apparently not correlated with any significant difference between the pea operon and the rRNA operons from tobacco, maize and Marchantia .     

RibAlign: a software tool and database for eubacterial phylogeny based on concatenated ribosomal protein subunits

Background  

Until today, analysis of 16S ribosomal RNA (rRNA) sequences has been the de-facto gold standard for the assessment of phylogenetic relationships among prokaryotes. However, the branching order of the individual phlya is not well-resolved in 16S rRNA-based trees. In search of an improvement, new phylogenetic methods have been developed alongside with the growing availability of complete genome sequences. Unfortunately, only a few genes in prokaryotic genomes qualify as universal phylogenetic markers and almost all of them have a lower information content than the 16S rRNA gene. Therefore, emphasis has been placed on methods that are based on multiple genes or even entire genomes. The concatenation of ribosomal protein sequences is one method which has been ascribed an improved resolution. Since there is neither a comprehensive database for ribosomal protein sequences nor a tool that assists in sequence retrieval and generation of respective input files for phylogenetic reconstruction programs, RibAlign has been developed to fill this gap.     

原核生物全基因组中16S rRNA基因的识别

【目的】识别原核生物全基因组中的16S rRNA基因。【方法】本文依据基因序列的GC碱基含量、碱基3-周期性和马尔可夫链3个方面的特性,构建了识别原核生物全基因组中16S rRNA基因的三层过滤模型。【结果】经检验,模型的特异性、敏感性和马修斯相关系数分别为99.58%、91.60%和91.49%。【结论】结果表明,本文所提出的方法可以高效、准确地识别出16S rRNA基因。     

Gene copy number variation and its significance in cyanobacterial phylogeny

ABSTRACT: BACKGROUND: In eukaryotes, variation in gene copy numbers is often associated with deleterious effects, but may also have positive effects. For prokaryotes, studies on gene copy number variation are rare. Previous studies have suggested that high numbers of rRNA gene copies can be advantageous in environments with changing resource availability, but further association of gene copies and phenotypic traits are not documented. We used one of the morphologically most diverse prokaryotic phyla to test whether numbers of gene copies are associated with levels of cell differentiation. RESULTS: We implemented a search algorithm that identified 44 genes with highly conserved copies across 22 fully sequenced cyanobacterial taxa. For two very basal cyanobacterial species, Gloeobacter violaceus and a thermophilic Synechococcus species, distinct phylogenetic positions previously found were supported by identical protein coding gene copy numbers. Furthermore, we found that increased ribosomal gene copy numbers showed a strong correlation to cyanobacteria capable of terminal cell differentiation. Additionally, we detected extremely low variation of 16S rRNA sequence copies within the cyanobacteria. We compared our results for 16S rRNA to three other eubacterial phyla (Chroroflexi, Spirochaetes and Bacteroidetes). Based on Bayesian phylogenetic inference and the comparisons of genetic istances, we could confirm that cyanobacterial 16S rRNA paralogs and orthologs show significantly stronger conservation than found in other eubacterial phyla. Conclusions: A higher number of ribosomal operons could potentially provide an advantage to terminally differentiated cyanobacteria. Furthermore, we suggest that 16S rRNA gene copies in cyanobacteria are homogenized by both concerted evolution and purifying selection. In addition, the small ribosomal subunit in cyanobacteria appears to evolve at extraordinary slow evolutionary rates, an observation that has been made previously for morphological characteristics of cyanobacteria.     

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.     

Quantification of microbial communities in near-surface and deeply buried marine sediments on the Peru continental margin using real-time PCR

Deeply buried marine sediments harbour a large fraction of all prokaryotes on Earth but it is still unknown which phylogenetic and physiological microbial groups dominate the deep biosphere. In this study real-time PCR allowed a comparative quantitative microbial community analysis in near-surface and deeply buried marine sediments from the Peru continental margin. The 16S rRNA gene copy numbers of prokaryotes and Bacteria were almost identical with a maximum of 10(8)-10(10) copies cm(-3) in the near-surface sediments. Archaea exhibited one to three orders of magnitude lower 16S rRNA gene copy numbers. The 18S rRNA gene of Eukarya was always at least three orders of magnitude less abundant than the 16S rRNA gene of prokaryotes. The 16S rRNA gene of the Fe(III)- and Mn(IV)-reducing bacterial family Geobacteraceae and the dissimilatory (bi)sulfite reductase gene (dsrA) of sulfate-reducing prokaryotes were abundant with 10(6)-10(8) copies cm(-3) in near-surface sediments but showed lower numbers and an irregular distribution in the deep sediments. The copy numbers of all genes decreased with sediment depth exponentially. The depth gradients were steeper for the gene copy numbers than for numbers of total prokaryotes (acridine orange direct counts), which reflects the ongoing degradation of the high-molecular-weight DNA with sediment age and depth. The occurrence of eukaryotic DNA also suggests DNA preservation in the deeply buried sediments.     

Discordant 16S and 23S rRNA gene phylogenies for the genus Helicobacter: implications for phylogenetic inference and systematics

Analysis of 16S rRNA gene sequences has become the primary method for determining prokaryotic phylogeny. Phylogeny is currently the basis for prokaryotic systematics. Therefore, the validity of 16S rRNA gene-based phylogenetic analyses is of fundamental importance for prokaryotic systematics. Discrepancies between 16S rRNA gene analyses and DNA-DNA hybridization and phenotypic analyses have been noted in the genus Helicobacter. To clarify these discrepancies, we sequenced the 23S rRNA genes for 55 helicobacter strains representing 41 taxa (>2,700 bases per sequence). Phylogenetic-tree construction using neighbor-joining, parsimony, and maximum likelihood methods for 23S rRNA gene sequence data yielded stable trees which were consistent with other phenotypic and genotypic methods. The 16S rRNA gene sequence-derived trees were discordant with the 23S rRNA gene trees and other data. Discrepant 16S rRNA gene sequence data for the helicobacters are consistent with the horizontal transfer of 16S rRNA gene fragments and the creation of mosaic molecules with loss of phylogenetic information. These results suggest that taxonomic decisions must be supported by other phylogenetically informative macromolecules, such as the 23S rRNA gene, when 16S rRNA gene-derived phylogeny is discordant with other credible phenotypic and genotypic methods. This study found Wolinella succinogenes to branch with the unsheathed-flagellum cluster of helicobacters by 23S rRNA gene analyses and whole-genome comparisons. This study also found intervening sequences (IVSs) in the 23S rRNA genes of strains of 12 Helicobacter species. IVSs were found in helices 10, 25, and 45, as well as between helices 31' and 27'. Simultaneous insertion of IVSs at three sites was found in H. mesocricetorum.     

Restriction enzyme analysis of Bacillus subtilis ribosomal ribonucleic acid genes.

The organization of the ribosomal ribonucleic acid (rRNA) genes (rDNA) of Bacillus subtilis was examined by cleaving the genome with several restriction endonucleases. The rDNA sequences were assayed by hybridization with purified radioactive rRNA's. Our interpretation of the resulting electrophoretic patterns is strengthened by an analysis of a fragment of B. subtilis rDNA cloned in Escherichia coli. The results indicated that there are eight rRNA operons in B. subtilis. Each operon contains one copy of the sequences coding for 16S, 23S, and 5S rRNA. The sequences coding for 5S rRNA were shown to be more closely linked to the 23S rRNA genes than to the 16S rRNA genes.     

rRNA operon copy number reflects ecological strategies of bacteria

Although natural selection appears to favor the elimination of gene redundancy in prokaryotes, multiple copies of each rRNA-encoding gene are common on bacterial chromosomes. Despite this conspicuous deviation from single-copy genes, no phenotype has been consistently associated with rRNA gene copy number. We found that the number of rRNA genes correlates with the rate at which phylogenetically diverse bacteria respond to resource availability. Soil bacteria that formed colonies rapidly upon exposure to a nutritionally complex medium contained an average of 5.5 copies of the small subunit rRNA gene, whereas bacteria that responded slowly contained an average of 1.4 copies. In soil microcosms pulsed with the herbicide 2,4-dichlorophenoxyacetic acid (2,4-D), indigenous populations of 2,4-D-degrading bacteria with multiple rRNA genes ( = 5.4) became dominant, whereas populations with fewer rRNA genes ( = 2.7) were favored in unamended controls. These findings demonstrate phenotypic effects associated with rRNA gene copy number that are indicative of ecological strategies influencing the structure of natural microbial communities.     

检测肠道细菌Feacalibacterium prausnitzii的三对PCR引物的特异性比较

【目的】比较3对基于16S rRNA基因、用于检测人肠道中重要细菌Feacalibacterium prausnitzii的引物(FPR-1/FPR-2、FPR-2F/Fprau645R和Fprau223F/Fprau420R)的特异性。【方法】用Clustal X比对每个引物与F.prausnitzii和其他细菌的16S rRNA基因的序列。在Ribosomal Database Project(RDP)数据库中使用Probe Match工具比较每个引物匹配的Faecalibacterium spp.序列数目。利用本课题组建立的中国人粪便菌群的16S rRNA基因全长文库的7255个克隆序列,用Simulated PCR(SPCR)预测每对引物检测到的F.prausnitzii和其他细菌的克隆数;用3对引物分别对代表克隆进行PCR扩增。用3对引物分别对14个健康人的粪便样品进行实时定量PCR。【结果】引物Fprau645R的3端最后一个碱基与非F.prausnitzii序列的错配度高于其它引物,它在RDP中匹配的Faecalibacterium spp.序列数占其匹配的细菌序列数的百分比(97.6%)显著高于其他引物。SPCR预测,3对引物检测到的F.prausnitzii克隆数均为1171左右;在检测到的非Faecalibacterium spp.克隆中,FPR-2F/Fprau645R主要是Subdoligranulum spp.,而FPR-1/FPR-2和Fprau223F/Fprau420R还有Oscillibacter spp.、Ruminococcus spp.和unclassified Ruminococcaceae等。真实PCR与SPCR的结果吻合。实时定量PCR中,FPR-1/FPR-2和Fprau223F/Fprau420R检测到的细菌数量高于FPR-2F/Fprau645R。【结论】3对引物能检测到F.prausnitzii和Subdoligranulum spp.,FPR-2F/Fprau645R的特异性优于FPR-1/FPR-2和Fprau223F/Fprau420R。     

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.     

A robust approach to estimate relative phytoplankton cell abundances from metagenomes

Phytoplankton account for >45% of global primary production, and have an enormous impact on aquatic food webs and on the entire Earth System. Their members are found among prokaryotes (cyanobacteria) and multiple eukaryotic lineages containing chloroplasts. Genetic surveys of phytoplankton communities generally consist of PCR amplification of bacterial (16S), nuclear (18S) and/or chloroplastic (16S) rRNA marker genes from DNA extracted from environmental samples. However, our appreciation of phytoplankton abundance or biomass is limited by PCR-amplification biases, rRNA gene copy number variations across taxa, and the fact that rRNA genes do not provide insights into metabolic traits such as photosynthesis. Here, we targeted the photosynthetic gene psbO from metagenomes to circumvent these limitations: the method is PCR-free, and the gene is universally and exclusively present in photosynthetic prokaryotes and eukaryotes, mainly in one copy per genome. We applied and validated this new strategy with the size-fractionated marine samples collected by Tara Oceans, and showed improved correlations with flow cytometry and microscopy than when based on rRNA genes. Furthermore, we revealed unexpected features of the ecology of these ecosystems, such as the high abundance of picocyanobacterial aggregates and symbionts in the ocean, and the decrease in relative abundance of phototrophs towards the larger size classes of marine dinoflagellates. To facilitate the incorporation of psbO in molecular-based surveys, we compiled a curated database of >18,000 unique sequences. Overall, psbO appears to be a promising new gene marker for molecular-based evaluations of entire phytoplankton communities.