Bacterial phylogeny based on 16S and 23S rRNA sequence analysis

Abstract: Molecular phylogeny increasingly supports the understanding of organismal relationships and provides the basis for the classification of microorganisms according to their natural affiliations. Comparative sequence analysis of ribosomal RNAs or the corresponding genes currently is the most widely used approach for the reconstruction of microbial phylogeny. The highly and less conserved primary and higher order structure elements of rRNAs document the history of microbial evolution and are informative for definite phylogenetic levels. An optimal alignment of the primary structures and a careful data selection are prerequisites for reliable phylogenetic conclusions. rRNA based phylogenetic trees can be reconstructed and the significance of their topologies evaluated by applying distance, maximum parsimony and maximum likelihood methods of phylogeny inference in comparison, and by fortuitous or directed resampling of the data set. Phylogenetic trees based on almost equivalent data sets of bacterial 23S and 16S rRNAs are in good agreement and their overall topologies are supported by alternative phylogenetic markers such as elongation factors and ATPase subunits. Besides their phylogenetic information content, the differently conserved primary structure regions of rRNAs provide target sites for specific hybridization probes which have been proven to be powerful tools for the identification of microbes on the basis of their phylogenetic relationships.     

Aligned 18S and insect phylogeny

The nuclear small subunit rRNA (18S) has played a dominant role in the estimation of relationships among insect orders from molecular data. In previous studies, 18S sequences have been aligned by unadjusted automated approaches (computer alignments that are not manually readjusted), most recently with direct optimization (simultaneous alignment and tree building using a program called "POY"). Parsimony has been the principal optimality criterion. Given the problems associated with the alignment of rRNA, and the recent availability of the doublet model for the analysis of covarying sites using Bayesian MCMC analysis, a different approach is called for in the analysis of these data. In this paper, nucleotide sequence data from the 18S small subunit rRNA gene of insects are aligned manually with reference to secondary structure, and analyzed under Bayesian phylogenetic methods with both GTR+I+G and doublet models in MrBayes. A credible phylogeny of Insecta is recovered that is independent of the morphological data and (unlike many other analyses of 18S in insects) not contradictory to traditional ideas of insect ordinal relationships based on morphology. Hexapoda, including Collembola, are monophyletic. Paraneoptera are the sister taxon to a monophyletic Holometabola but weakly supported. Ephemeroptera are supported as the sister taxon of Neoptera, and this result is interpreted with respect to the evolution of direct sperm transfer and the evolution of flight. Many other relationships are well-supported but several taxa remain problematic, e.g., there is virtually no support for relationships among orthopteroid orders. A website is made available that provides aligned 18S data in formats that include structural symbols and Nexus formats.     

16S rRNA及rpoC1基因用于螺旋藻、节旋藻系统发育的比较北大核心CSCD

【目的】利用16S rRNA和rpoC1基因分子标记研究螺旋藻、节旋藻的系统发育关系,并对其区分能力进行比较。【方法】以84株螺旋藻、节旋藻为研究对象,对其进行16S rRNA、rpoC1基因序列的扩增、测序及分析,并对构建的系统发育树进行对比。【结果】rpoC1基因序列保守位点所占比例49.7%、平均G+C百分含量47.7%和序列相似度76%–100%明显低于16S rRNA基因序列的79.4%、55.6%和91%–100%,其变异程度高于16S rRNA基因;基于16S rRNA、rpoC1基因构建的系统发育NJ树拓扑结构基本一致,84株实验藻株分为2个属3个类群,其中仅F-351、F-904-2、F-1070和TJBC14-1藻株为螺旋藻,其余均为节旋藻;虽然2个基因都不能区分形态种和地理种,但rpoC1基因NJ树的置信度(100%)高于16S rRNA基因(99%),属内分群效果也明显优于16S rRNA基因。【结论】支持了螺旋藻、节旋藻为两个不同属的结论,且在属内分类时rpoC1基因比16S rRNA基因具有更高的区分度。     

基于16S rRNA和Cyt b基因序列的鰶亚科系统发育研究

对鰶亚科4属5种鱼类的线粒体基因组16S rRNA和Cyt b基因片段序列进行序列比较和系统发育关系分析。结果显示:5种鰶亚科鱼类的16S rRNA和Cy tb基因片段同源序列长度分别为525 bp和402 bp,序列联合后的序列总长度为927 bp,其中多态位点178个,简约信息位点123个。选取太平洋鲱Clupea pallasii和大西洋鲱C.harengus为外类群,采用邻接法(NJ)、最大简约法(MP)、最大似然法(ML)和贝叶斯法(BI)分别对2个基因片段序列进行了聚类分析,并联合2个基因片段利用邻接法、最大简约法和贝叶斯法进行分析。系统发育分析显示:斑鰶Konosirus punctatus与花鰶Clupanodon thriss亲缘关系最近,分布于美洲大陆的真鰶属Dorosoma鱼类与印度洋、太平洋分布的斑鰶属、花鰶属和海鰶属Nematalosa鱼类亲缘关系较远。     

Contributions to molecular systematics of water scavenger beetles (Hydrophilidae,Coleoptera)

Phylogenetic relationships within Hydrophilidae were examined by analyses of separate and combined nuclear and mitochondrial markers (28S rRNA, 18S rRNA, 16S rRNA, 12S rRNA, COI and COII genes). The preferred (Bayesian) tree topology suggests a sister group relationship between Spercheidae and Hydrophilidae, supporting the ‘hydrophilid lineage’; Epimetopidae are placed on the base of the ‘helophorid branch’, the monophyly of Sphaeridiinae is highly supported, nested deeply within Hydrophilidae closest to Enochrus, making Hydrophilinae and Acidocerini paraphyletic; Hydrobius appears as sister taxon to (Hydrochara + Hydrophilus) without a closer relationship to Acidocerini; the hydrophiloid–histeroid sister group relationship is confirmed. The topology of several taxa remains contradictory, and requires further investigations with a larger taxon sampling and additional molecular markers.     

Phylogenetic position of Salinibacter ruber based on concatenated protein alignments

A total of 22 genes from the genome of Salinibacter ruber strain M31 were selected in order to study the phylogenetic position of this species based on protein alignments. The selection of the genes was based on their essential function for the organism, dispersion within the genome, and sufficient informative length of the final alignment. For each gene, an individual phylogenetic analysis was performed and compared with the resulting tree based on the concatenation of the 22 genes, which rendered a single alignment of 10,757 homologous positions. In addition to the manually chosen genes, an automatically selected data set of 74 orthologous genes was used to reconstruct a tree based on 17,149 homologous positions. Although single genes supported different topologies, the tree topology of both concatenated data sets was shown to be identical to that previously observed based on small subunit (SSU) rRNA gene analysis, in which S. ruber was placed together with Bacteroidetes. In both concatenated data sets the bootstrap was very high, but an analysis with a gradually lower number of genes indicated that the bootstrap was greatly reduced with less than 12 genes. The results indicate that tree reconstructions based on concatenating large numbers of protein coding genes seem to produce tree topologies with similar resolution to that of the single 16S rRNA gene trees. For classification purposes, 16S rRNA gene analysis may remain as the most pragmatic approach to infer genealogic relationships.     

Corresponding 16S rRNA gene segments in Rhizobiaceae and Aeromonas yield discordant phylogenies

Previous evidence has indicated that the 16S rRNA genes in certain species of Aeromonas may have a history of lateral transfer and recombination. A comparative analysis of patterns of 16S nucleotide sequence polymorphism among species of Rhizobium and Agrobacterium was conducted to determine if there is similar evidence for chimeric 16S genes in members of the Rhizobiaceae. Results from phylogenetic analyses and comparison of patterns of nucleotide sequence polymorphism in portions of rhizobial 16S genes revealed the same type of segment-dependent polymorphic site partitioning that was previously reported for Aeromonas. These results support the hypothesis that certain 16S segments in rhizobia may have a history of lateral transfer and recombination.Abbreviations 16S rRNA 16S ribosomal ribonucleic acid - 16S the 16S rRNA gene     

Evaluation of a novel 16S rRNA/tRNA mitochondrial marker for the identification and phylogenetic analysis of shrimp species belonging to the superfamily Penaeoidea

In this study, we infer the phylogenetic relationships within commercial shrimp using sequence data from a novel mitochondrial marker consisting of an approximately 530-bp region of the 16S ribosomal RNA (rRNA)/transfer RNA (tRNA)^Val genes compared with two other mitochondrial genes: 16S rRNA and cytochrome c oxidase I (COI). All three mitochondrial markers were considerably AT rich, exhibiting values up to 78.2% for the species Penaeus monodon in the 16S rRNA/tRNA^Val genes, notably higher than the average among other Malacostracan mitochondrial genomes. Unlike the 16S rRNA and COI genes, the 16S rRNA/tRNA^Val marker evidenced that Parapenaeus is more closely related to Metapenaeus than to Solenocera, a result that seems to be more in agreement with the taxonomic status of these genera. To our knowledge, our study using the 16S rRNA/tRNA^Val gene as a marker for phylogenetic analysis offers the first genetic evidence to confirm that Pleoticus muelleri and Solenocera agassizi constitute a separate group and that they are more related to each other than to genera belonging to the family Penaeidae. The 16S rRNA/tRNA^Val region was also found to contain more variable sites (56%) than the other two regions studied (33.4% for the 16S rRNA region and 42.7% for the COI region). The presence of more variable sites in the 16S rRNA/tRNA^Val marker allowed the interspecific differentiation of all 19 species examined. This is especially useful at the commercial level for the identification of a large number of shrimp species, particularly when the lack of morphological characteristics prevents their differentiation.     

A comparison of primer sets for detecting 16S rRNA and hydrazine oxidoreductase genes of anaerobic ammonium-oxidizing bacteria in marine sediments

Published polymerase chain reaction primer sets for detecting the genes encoding 16S rRNA gene and hydrazine oxidoreductase (hzo) in anammox bacteria were compared by using the same coastal marine sediment samples. While four previously reported primer sets developed to detect the 16S rRNA gene showed varying specificities between 12% and 77%, an optimized primer combination resulted in up to 98% specificity, and the recovered anammox 16S rRNA gene sequences were >95% sequence identical to published sequences from anammox bacteria in the Candidatus “Scalindua” group. Furthermore, four primer sets used in detecting the hzo gene of anammox bacteria were highly specific (up to 92%) and efficient, and the newly designed primer set in this study amplified longer hzo gene segments suitable for phylogenetic analysis. The optimized primer set for the 16S rRNA gene and the newly designed primer set for the hzo gene were successfully applied to identify anammox bacteria from marine sediments of aquaculture zone, coastal wetland, and deep ocean where the three ecosystems form a gradient of anthropogenic impact. Results indicated a broad distribution of anammox bacteria with high niche-specific community structure within each marine ecosystem.     

Unusual intragenomic and interspecific variability of <Emphasis Type="Italic">Geobacillus</Emphasis> 16S-23S rRNA internal transcribed spacers

The aim of this study was to evaluate the inter-and intraspecific as well as intragenomic variability of Geobacillus 16S–23S rRNA internal transcribed spacers without tRNA genes and to compare these sequences with sequences bearing tRNA genes. In this study the structural analysis was performed in a unique way because the length and the sequence of the structural blocks were adjusted to fit the structure of 16S–23S rRNA internal transcribed spacers of five different Geobacillus species. Our study demonstrated the mosaic-like structure of 16S–23S rRNA internal transcribed spacers in Geobacillus. Some characteristics of these spacers of geobacilli were not previously reported for other bacteria: unusually short conserved sequence in the 5′ end region, some identical conserved blocks in both 5′ and 3′ regions of 16S–23S rRNA internal transcribed spacers, the same sequence blocks in both 16S–23S and 23S–5S rRNA intergenic spacers. Our study demonstrated quite uniform arrangement of the sequence blocks in Geobacillus thermodenitrificans. This species diverged early in the phylogenetic tree of the genus Geobacillus. For the phylogenetically recent species Geobacillus kaustophilus and Geobacillus lituanicus the low inter-and intraspecific, but high intragenomic variability, as a consequence of recent phylogenetic events, was established.     

Phylogenetic analysis based on 28S rRNA of Babesia spp. in ruminants in China

Molecular phylogenetic analyses are mainly based on the small ribosomal RNA subunit (18S rRNA), internal transcribed spacer regions, and other molecular markers. We compared the phylogenetic relationships of Babesia spp. using large subunit ribosomal RNA, i.e., 28S rRNA, and the united 28S + 18S rRNA sequence fragments from 11 isolates of Babesia spp. collected in China. Due to sequence length and variability, the 28S rRNA gene contained more information than the 18S rRNA gene and could be used to elucidate the phlyogenetic relationships of B. motasi, B. major, and B. bovis. Thus, 28S rRNA is another candidate marker that can be used for the phylogenetic analysis of Babesia spp. However, the united fragment (28S + 18S) analysis provided better supported phylogenetic relationships than single genes for Babesia spp. in China.     

Among-site rate variation and phylogenetic analysis of 12S rRNA in sigmodontine rodents

We analyze sequences from two mitochondrial genes, cytochrome b (cyt b) and 12S rRNA (12S), for a group of sigmodontine rodents among which phylogenetic relationships are well understood based on concordance of morphological, chromosomal, allozyme, and other DNA data sets. Because these two genes are physically linked on the nonrecombining mitochondrial genome, they necessarily share the same history. Phylogenetic analysis of the cyt b gene recovers the well-corroborated relationships, generally with strong support. None of the methods that we employed, including variously weighted parsimony, neighbor joining on both single-rate and gamma-corrected distances, and maximum likelihood, were able to recover these relationships for the 12S gene. Parsimony analyses of the 12S data resulted in a relatively strongly supported placement of Peromyscus eremicus that conflicts with that suggested by cyt b and all other data. There is extreme among-site rate variation in the 12S sequences and moderate levels in the cyt b sequences. This highly skewed distribution of rates in the 12S gene makes phylogenetic analyses of these sequences particularly susceptible to the misleading effects of nonindependence and other nonrandom noise, suggesting that phylogenetic analyses of data sets that contain a great deal of among-site rate variation be interpreted with caution.      

Molecular identification and phylogenetic relationships of seven Indian Sciaenids (Pisces: Perciformes,Sciaenidae) based on 16S rRNA and cytochrome c oxidase subunit I mitochondrial genes

The partial sequences of 16S rRNA and cytochrome c oxidase subunit I (COI) mitochondrial genes were analyzed for species identification and phylogenetic relationships among the commercially important Indian sciaenids (Otolithes cuvieri, Otolithes ruber, Johnius dussumieri, Johnius elongatus, Johnieops vogleri, Otolithoides biauritus and Protonibea diacanthus). Sequence analysis of both genes revealed that the seven species fell into three distinct groups, which were genetically distant from each other and exhibited identical phylogenetic resolution. Partial sequences of both the genes provided sufficient phylogenetic information to distinguish the seven sciaenids indicating the usefulness of mtDNA-based approach in species identification.     

Multiple Substitutions Affect the Phylogenetic Utility of Cytochrome b and 12S rDNA Data: Examining a Rapid Radiation in Leporid (Lagomorpha) Evolution

Partial sequences of two mitochondrial genes, the 12S ribosomal gene (739 bp) and the cytochrome b gene (672 bp), were analyzed in hopes of reconstructing the evolutionary relationships of 11 leporid species, representative of seven genera. However, partial cytochrome b sequences were of little phylogenetic value in this study. A suite of pairwise comparisons between taxa revealed that at the intergeneric level, the cytochrome b gene is saturated at synonymous coding positions due to multiple substitution events. Furthermore, variation at the nonsynonymous positions is limited, rendering the cytochrome b gene of little phylogenetic value for assessing the relationships between leporid genera. If the cytochrome b data are analyzed without accounting for these two classes of nucleotides (i.e., synonymous and nonsynonymous sites), one may incorrectly conclude that signal exists in the cytochrome b data. The mitochondrial 12S rRNA gene, on the other hand, has not experienced excessive saturation at either stem or loop positions. Phylogenies reconstructed from the 12S rDNA data support hypotheses based on fossil evidence that African rock rabbits (Pronolagus) are outside of the main leporid stock and that leporids experienced a rapid radiation. However, the molecular data suggest that this radiation event occurred in the mid-Miocene several millions of years earlier than the Pleistocene dates suggested by paleontological evidence. Received: 23 April 1998 / Accepted: 14 May 1998     

Borrelia tanukii Sp. Nov. and Borrelia turdae Sp. Nov. Found from Ixodid Ticks in Japan: Rapid Species Identification by 16S rRNA Gene-Targeted PCR Analysis

Based on the results of RFLP-ribotyping, whole DNA/DNA hybridization and phylogenetic analysis of the 16S rRNA gene, we previously defined two genomic groups of spirochetes closely related to Borrelia burgdorferi sensu lato: group Hk501 for strains isolated from Ixodes tanuki ticks and group Ya501 for strains isolated from Ixodes turdus ticks. In this report, we propose that group Hk501 should be classified as Borrelia tanukii sp. nov. and group Ya501 as Borrelia turdae sp. nov. The alignment of previously published Borrelia 16S rRNA gene sequences led us to design species-specific PCR primer sets. The primers allowed the rapid identification of B. tanukii and B. turdae.     

Networks and groups within the genus Neisseria: analysis of argF, recA, rho, and 16S rRNA sequences from human Neisseria species.

To understand the pattern of nucleotide sequence variation among bacteria that frequently exchange chromosomal genes, we analyzed sequences of the recA, argF, and rho genes, as well as part of the small-subunit (16S) rRNA gene, from about 50 isolates of human commensal Neisseria species and the pathogenic N. meningitidis and N. gonorrhoeae. Almost all isolates of these species could be assigned to five phylogenetic groups that are found for all genes examined and generally are supported by high bootstrap values. In contrast, the phylogenetic relationships among groups varied according to the gene analyzed with notable incongruences involving N. cinerea and N. lactamica. Further analysis using split decomposition showed that for each gene, including 16S rRNA, the patterns of sequence divergence within N. meningitidis and closely related species were inconsistent with a bifurcating treelike phylogeny and better represented by an interconnected network. These data indicate that the human commensal Neisseria species can be separated into discrete groups of related species but that the relationships both within and among these groups, including those reconstructed using 16S rRNA, have been distorted by interspecies recombination events.     

Phylogenetic analysis, structural evolution and functional divergence of the 12-oxo-phytodienoate acid reductase gene family in plants

Background

Whenever different data sets arrive at conflicting phylogenetic hypotheses, only testable causal explanations of sources of errors in at least one of the data sets allow us to critically choose among the conflicting hypotheses of relationships. The large (28S) and small (18S) subunit rRNAs are among the most popular markers for studies of deep phylogenies. However, some nodes supported by this data are suspected of being artifacts caused by peculiarities of the evolution of these molecules. Arthropod phylogeny is an especially controversial subject dotted with conflicting hypotheses which are dependent on data set and method of reconstruction. We assume that phylogenetic analyses based on these genes can be improved further i) by enlarging the taxon sample and ii) employing more realistic models of sequence evolution incorporating non-stationary substitution processes and iii) considering covariation and pairing of sites in rRNA-genes.

Results

We analyzed a large set of arthropod sequences, applied new tools for quality control of data prior to tree reconstruction, and increased the biological realism of substitution models. Although the split-decomposition network indicated a high noise content in the data set, our measures were able to both improve the analyses and give causal explanations for some incongruities mentioned from analyses of rRNA sequences. However, misleading effects did not completely disappear.

Conclusion

Analyses of data sets that result in ambiguous phylogenetic hypotheses demand for methods, which do not only filter stochastic noise, but likewise allow to differentiate phylogenetic signal from systematic biases. Such methods can only rely on our findings regarding the evolution of the analyzed data. Analyses on independent data sets then are crucial to test the plausibility of the results. Our approach can easily be extended to genomic data, as well, whereby layers of quality assessment are set up applicable to phylogenetic reconstructions in general.     

16S rRNA phylogenetic analysis and quantification of <Emphasis Type="Italic">Korarchaeota</Emphasis> indigenous to the hot springs of Kamchatka,Russia

The candidate archaeal division Korarchaeota is known primarily from deeply branching sequences of 16S rRNA genes PCR-amplified from hydrothermal springs. Parallels between the phylogeny of these genes and the geographic locations where they were identified suggested that Korarchaeota exhibit a high level of endemism. In this study, the influence of geographic isolation and select environmental factors on the diversification of the Korarchaeota was investigated. Fourteen hot springs from three different regions of Kamchatka, Russia were screened by PCR using Korarchaeota-specific and general Archaea 16S rRNA gene-targeting primers, cloning, and sequencing. Phylogenetic analyses of these sequences with Korarchaeota 16S rRNA sequences previously identified from around the world suggested that all Kamchatka sequences cluster together in a unique clade that subdivides by region within the peninsula. Consistent with endemism, 16S rRNA gene group-specific quantitative PCR of all Kamchatka samples detected only the single clade of Korarchaeota that was found by the non-quantitative PCR screening. In addition, their genes were measured in only low numbers; small Korarchaeota populations would present fewer chances for dispersal to and colonization of other sites. Across the entire division of Korarchaeota, common geographic locations, temperatures, or salinities of identification sites united sequence clusters at different phylogenetic levels, suggesting varied roles of these factors in the diversification of Korarchaeota.     

Taxonomic re-evaluation of Aphanizomenon flos-aquae NH-5 based on morphology and 16S rRNA gene sequences

The taxonomy of Aphanizomenon flos-aquae strain NH-5, a producer of cyanotoxins, was re-evaluated by comparison with six other Aphanizomenon strains using morphological characteristics and 16S rRNA gene sequences. Strain NH-5 was concluded to be improperly identified as Aph. flos-aquae based upon (1) lack of bundle formation in the trichomes, (2) location of akinetes next to heterocytes, (3) lower similarities (less than 97.5%) in the 16S rRNA gene sequences relative to Aph. flos-aquae strains, and (4) comparison within a phylogenetic tree constructed from 16S rRNA gene sequences. The Aphanizomenon strains investigated in this study are classified to four morphological groups as described by the classical taxonomy of Komárek & Kovácik (1989). This classification was supported from the phylogenetic results of 16S rRNA gene sequences. This study also discusses the generic boundaries between Aphanizomenon and Anabaena.     

Molecular evolution of the mitochondrial 12S rRNA in Ungulata (mammalia)

The complete 12S rRNA gene has been sequenced in 4 Ungulata (hoofed eutherians) and 1 marsupial and compared to 38 available mammalian sequences in order to investigate the molecular evolution of the mitochondrial small-subunit ribosomal RNA molecule. Ungulata were represented by one artiodactyl (the collared peccary, Tayassu tajacu, suborder Suiformes), two perissodactyls (the Grevy's zebra, Equus grevyi, suborder Hippomorpha; the white rhinoceros, Ceratotherium simum, suborder Ceratomorpha), and one hyracoid (the tree hyrax, Dendrohyrax dorsalis). The fifth species was a marsupial, the eastern gray kangaroo (Macropus giganteus). Several transition/transversion biases characterized the pattern of changes between mammalian 12S rRNA molecules. A bias toward transitions was found among 12S rRNA sequences of Ungulata, illustrating the general bias exhibited by ribosomal and protein-encoding genes of the mitochondrial genome. The derivation of a mammalian 12S rRNA secondary structure model from the comparison of 43 eutherian and marsupial sequences evidenced a pronounced bias against transversions in stems. Moreover, transversional compensatory changes were rare events within double-stranded regions of the ribosomal RNA. Evolutionary characteristics of the 12S rRNA were compared with those of the nuclear 18S and 28S rRNAs. From a phylogenetic point of view, transitions, transversions and indels in stems as well as transversional and indels events in loops gave congruent results for comparisons within orders. Some compensatory changes in double-stranded regions and some indels in single-stranded regions also constituted diagnostic events. The 12S rRNA molecule confirmed the monophyly of infraorder Pecora and order Cetacea and demonstrated the monophyly of suborder Suiformes. However, the monophyly of the suborder Ruminantia was not supported, and the branching pattern between Cetacea and the artiodactyl suborders Ruminantia and Suiformes was not established. The monophyly of the order Perissodactyla was evidenced, but the relationships between Artiodactyla, Cetacea, and Perissodactyla remained unresolved. Nevertheless, we found no support for a Perissodactyla + Hyracoidea clade, neither with distance approach, nor with parsimony reconstruction. The 12S rRNA was useful to solve intraordinal relationships among Ungulata, but it seemed to harbor too few informative positions to decipher the bushlike radiation of some Ungulata orders, an event which has most probably occurred in a short span of time between 55 and 70 MYA. Correspondence to: E. Douzery