Estimating Substitution Rates in Ribosomal RNA Genes

A model is introduced describing nucleotide substitution in ribosomal RNA (rRNA) genes. In this model, substitution in the stem and loop regions of rRNA is modeled with 16- and four-state continuous time Markov chains, respectively. The mean substitution rates at nucleotide sites are assumed to follow gamma distributions that are different for the two types of regions. The simplest formulation of the model allows for explicit expressions for transition probabilities of the Markov processes to be found. These expressions were used to analyze several 16S-like rRNA genes from higher eukaryotes with the maximum likelihood method. Although the observed proportion of invariable sites was only slightly higher in the stem regions, the estimated average substitution rates in the stem regions were almost two times as high as in the loop regions. Therefore, the degree of site heterogeneity of substitution rates in the stem regions seems to be higher than in the loop regions of animal 16S-like rRNAs due to presence of a few rapidly evolving sites. The model appears to be helpful in understanding the regularities of nucleotide substitution in rRNAs and probably minimizing errors in recovering phylogeny for distantly related taxa from these genes.     

Distribution of substitution rates and location of insertion sites in the tertiary structure of ribosomal RNA

The relative substitution rate of each nucleotide site in bacterial small subunit rRNA, large subunit rRNA and 5S rRNA was calculated from sequence alignments for each molecule. Two-dimensional and three-dimensional variability maps of the rRNAs were obtained by plotting the substitution rates on secondary structure models and on the tertiary structure of the rRNAs available from X-ray diffraction results. This showed that the substitution rates are generally low near the centre of the ribosome, where the nucleotides essential for its function are situated, and that they increase towards the surface. An inventory was made of insertions characteristic of the Archaea, Bacteria and Eucarya domains, and for additional insertions present in specific eukaryotic taxa. All these insertions occur at the ribosome surface. The taxon-specific insertions seem to arise randomly in the eukaryotic evolutionary tree, without any phylogenetic relatedness between the taxa possessing them.     

The 23S ribosomal RNA higher-order structure of Pseudomonas cepacia and other prokaryotes

A 23S ribosomal RNA gene of Pseudomonas cepacia has been cloned and sequenced. A general higher-order structure model based on earlier published models has been derived from comparative analysis of 23S-like rRNAs of eubacteria, archaebacteria, organelles and eukaryotes. Differences between the previous models were carefully analyzed and controversial regions evaluated. Moderately large insertions and deletions have been found at new points in the secondary structure. The analysis of 50 published as well as unpublished 23S rRNA sequences provide additional proof for six of the seven previously suggested tertiary interactions within the 23S rRNA. P. cepacia is the first representative of the beta subgroup of the Proteobacteria phylum whose 23S rRNA has been sequenced. A tree reflecting evolutionary relationships of prokaryotes was constructed. The topology of this tree is in good agreement with the 16S rRNA tree.     

Molecular phylogenetic evidence that the phylum Haplosporidia has an alveolate ancestry

The phylogenetic position of the phylum Haplosporidia among other protists was investigated with the complete 16S-like rRNA gene sequences from two species in the phylum: Haplosporidium nelsoni, a parasite of oysters, and Minchinia teredinis, a parasite of shipworms. Because the lack of obvious morphological homologies with other protists hampered decisions regarding taxonomic composition for sequence alignment and phylogenetic analysis, the complete sequences for these two haplosporidians were directed as search queries to the blast/ncbi.nlm.nih.gov electronic mail server. The results of this heuristic similarity search provided a basis for constructing a preliminary higher-taxonomic-level analysis comparing the haplosporidians with species from the slime molds, fungi, algae, amoebae, ciliates, dinoflagellates, and apicomplexans. Maximum parsimony yielded equivocal results, whereas transversionally weighted parsimony suggested an affinity with the alveolates (i.e., the ciliates, dinoflagellates, and apicomplexans). Multiple alignment of the two haplosporidian sequences against 17 taxa in a secondary analysis focusing on the alveolates and subsequent parsimony analysis placed the phylum Haplosporidia as a monophyletic group within the Alveolata and as a taxon of equal rank with the other three alveolate phyla. The precise placement within the Alveolata was sensitive to weighting.      

The unusually long small subunit ribosomal RNA gene found in amitochondriate amoeboflagellate Pelomyxa palustris: its rRNA predicted secondary structure and phylogenetic implication

In order to ascertain a phylogenetic position of the freshwater amitochondriate amoeboflagellate Pelomyxa palustris its small subunit (SSU) rRNA gene was amplified and sequenced. It was shown to be 3502 bp long. The predicted secondary structure of its rRNA includes at least 16 separate expansion zones located in all the variable regions (V1-V9), as well as in some conservative gene regions. Most insertions are represented by sequences of low complexity that have presumably arisen by a slippage mechanism. Relatively conservative, uniformly positioned motifs contained in regions V4 and V7, as well as in some others, made it possible to perform folding. In maximum likelihood, maximum parsimony, and neighbor-joining trees, P. palustris tends to cluster with amitochondriate and secondary lost mitochondria amoebae and amoeboflagellates Entamoeba, Endolimax nana, and Phreatamoeba balamuthi, comprising together with them and aerobic lobose amoebae Vannella, Acanthamoeba, Balamuthia, and Hartmannella a monophyletic cluster. Another pelobiont, Mastigamoeba invertens, does not belong to this cluster. No specific similarity was discovered between the SSU rRNA of P. palustris and amitochondriate taxa of 'Archezoa': Diplomonada, Parabasalia, Microsporidia. Pelomyxa palustris SSU rRNA does not occupy a basal position in the phylogenetic trees and could be ascribed to the so-called eukaryotic 'crown' group if the composition of the latter were not so sensitive to the methods of tree building. Thus, molecular and morphological data suggest that P. palustris represents a secondarily modified eukaryotic lineage.     

The characterization of enzymatically amplified eukaryotic 16S-like rRNA-coding regions

Polymerase chain reaction conditions were established for the in vitro amplification of eukaryotic small subunit ribosomal (16S-like) rRNA genes. Coding regions from algae, fungi, and protozoa were amplified from nanogram quantities of genomic DNA or recombinant plasmids containing rDNA genes. Oligodeoxynucleotides that are complementary to conserved regions at the 5' and 3' termini of eukaryotic 16S-like rRNAs were used to prime DNA synthesis in repetitive cycles of denaturation, reannealing, and DNA synthesis. The fidelity of synthesis for the amplification products was evaluated by comparisons with sequences of previously reported rRNA genes or with primer extension analyses of rRNAs. Fewer than one error per 2000 positions were observed in the amplified rRNA coding region sequences. The primary structure of the 16S-like rRNA from the marine diatom, Skeletonema costatum, was inferred from the sequence of its in vitro amplified coding region.     

Interrelationships among major protistan groups based on a parsimony network of 5S rRNA sequences

To test the validity of the maximum parsimony approach to discern protistan interrelationships, we have derived an optimal network of 16S-like rRNA sequences using our parsimony algorithm and compared it with those reported using the distance matrix method. We have also derived an optimal network topology of 50 5S rRNA sequences through an interactive search using our algorithm. In both these networks, the kinetoplastids and euglenoids form a linkage group with Dictyostelium emerging from its neighbourhood. The cryptophytes, dinoflagellates and chromophytes and green algae emerge as independent lines suggesting that plastids arose more than once during protistan evolution. The large 5S rRNA tree further indicates independent origins of mesozoa and metazoa; kinetoplastids and ciliates; and diphyletic origin of fungi. Comparatively close positions of charales and land plants, chytrids and Zygomycetes, Physarum and amoeba, and red algae and green algae are also seen in this network.     

Members of the Cytophaga–Flavobacterium–Bacteroides phylum as intracellular bacteria of acanthamoebae: proposal of 'Candidatus Amoebophilus asiaticus'

Three Gram-negative, rod-shaped bacteria that were found intracellularly in two environmental and one clinical Acanthamoeba sp. isolates were analysed. Two endocytobionts showing a parasitic behaviour were propagated successfully outside their amoebal host cells and were identified subsequently by comparative 16S rRNA sequence analysis as being most closely affiliated with Flavobacterium succinicans (99% 16S rRNA sequence similarity) or Flavobacterium johnsoniae (98% 16S rRNA sequence similarity). One endocytobiont could neither be cultivated outside its original Acanthamoeba host ( Acanthamoeba sp. TUMSJ-321) nor transferred into other amoebae. Electron microscopy revealed that the amoebal trophozoites and cysts were almost completely filled with cells of this endosymbiont which are surrounded by a host-derived membrane. According to 16S rRNA sequence analysis, this endosymbiont could also be assigned to the Cytophaga – Flavobacterium – Bacteroides (CFB) phylum, but was not closely affiliated to any recognized species within this phylogenetic group (less than 82% 16S rRNA sequence similarity). Identity and intracellular localization of this endosymbiont were confirmed by application of a specific fluorescently labelled 16S rRNA-targeted probe. Based on these findings, we propose classification of this obligate Acanthamoeba endosymbiont as ' Candidatus Amoebophilus asiaticus'. Comparative 18S rRNA sequence analysis of the host of ' Candidatus Amoebophilus asiaticus' revealed its membership with Acanthamoeba 18S rDNA sequence type T4 that comprises the majority of all Acanthamoeba isolates.     

Secondary structure models of the nuclear internal transcribed spacer regions and 5.8S rRNA in Calciodinelloideae (Peridiniaceae) and other dinoflagellates

Secondary structure models of the 5.8S rRNA and both internal transcribed spacers (ITS1 and ITS2) are proposed for Calciodinelloideae (Peridiniaceae) and are also plausible for other dinoflagellates. The secondary structure of the 5.8S rRNA corresponds to previously developed models, with two internal paired regions and at least one 5.8S rRNA–28S rRNA interaction. A general secondary structure model of ITS1 for Calciodinelloideae (and other dinoflagellates), consisting of an open multibranch loop with three major helices, is proposed. The homology of these paired regions with those found in other taxa, published in previous studies (e.g. yeast, green algae and Platyhelmithes) remains to be determined. Finally, a general secondary structure model of ITS2 for Calciodinelloideae (and other dinoflagellates) is reconstructed. Based on the 5.8S rRNA–28S rRNA interaction, it consists of a closed multibranch loop, with four major helices. At least helix III and IV have homology with paired regions found in other eukaryotic taxa (e.g. yeast, green algae and vertebrates). Since the secondary structures of both ITS regions are more conserved than the nucleotide sequences, their analysis helps in understanding molecular evolution and increases the number of structural characters. Thus, the structure models developed in this study may be generally useful for future phylogenetic analyses.     

Molecular phylogenetic profiling of prokaryotic communities in guts of termites with different feeding habits

Termites are an important group of terrestrial insects that harbor an abundant gut microbiota, many of which contribute to digestion, termite nutrition and gas (CH(4), CO(2) and H(2)) emission. With 2200 described species, termites also provide a good model to study relationships between host diet and gut microbial community structure and function. We examined the relationship between diet and gut prokaryotic community profiles in 24 taxonomically and nutritionally diverse species of termites by using nucleic acid probes targeting 16S-like ribosomal RNAs. The relative abundance of domain-specific 16S-like rRNAs recovered from gut extracts varied considerably (ranges: Archaea (0-3%); Bacteria (15-118%)). Although Bacteria were always detectable and the most abundant, differences in domain-level profiles were correlated with termite diet, as evidenced by higher relative abundances of Archaea in guts of soil-feeding termites, compared to those of wood-feeding species in the same family. The oligonucleotide probes also readily distinguished gut communities of wood-feeding taxa in the family Termitidae (higher termites) from those of other wood-feeding termite families (lower termites). The relative abundances of 16S-like archaeal rRNA in guts were positively correlated with rates of methane emission by live termites, and were consistent with previous work linking high relative rates of methanogenesis with the soil (humus)-feeding habit. Probes for methanogenic Archaea detected members of only two families (Methanobacteriaceae and Methanosarcinaceae) in termite guts, and these typically accounted for 60% of the all archaeal probe signal. In four species of termites, Methanosarcinaceae were dominant, a novel observation for animal gut microbial communities, but no clear relationship was apparent between methanogen family profiles and termite diet or taxonomy.     

Higher order structure in ribosomal RNA.

The only reliable general method currently available for determining precise higher order structure in the large ribosomal RNAs is comparative sequence analysis. The method is here applied to reveal 'tertiary' structure in the 16S-like rRNAs, i.e. structure more complex than simple double-helical, secondary structure. From a list of computer-generated potential higher order interactions within 16S rRNA one such interaction considered likely was selected for experimental test. The putative interaction involves a Watson-Crick one to one correspondence between positions 570 and 866 in the molecule (E. coli numbering). Using existing oligonucleotide catalog information several organisms were selected whose 16S rRNA sequences might test the proposed co-variation. In all of the (phylogenetically independent) cases selected, full sequence evidence confirms the predicted one to one (Watson-Crick) correspondence. An interaction between positions 570 and 866 is, therefore, considered proven phylogenetically.     

Novel bacterial endosymbionts of Acanthamoeba spp. related to the Paramecium caudatum symbiont Caedibacter caryophilus

Acanthamoebae are increasingly being recognized as hosts for obligate bacterial endosymbionts, most of which are presently uncharacterized. In this study, the phylogeny of three Gram-negative, rod-shaped endosymbionts and their Acanthamoeba host cells was analysed by the rRNA approach. Comparative analyses of 16S rDNA sequences retrieved from amoebic cell lysates revealed that the endosymbionts of Acanthamoeba polyphaga HN-3, Acanthamoeba sp. UWC9 and Acanthamoeba sp. UWE39 are related to the Paramecium caudatum endosymbionts Caedibacter caryophilus, Holospora elegans a n d Holospora obtusa . With overall 16S rRNA sequence similarities to their closest relative, C. caryophilus , of between 87% and 93%, these endosymbionts represent three distinct new species. In situ hybridization with fluorescently labelled endosymbiont-specific 16S rRNA-targeted probes demonstrated that the retrieved 16S rDNA sequences originated from the endosymbionts and confirmed their intracellular localization. We propose to classify provisionally the endosymbiont of Acanthamoeba polyphaga HN-3 as ' Candidatus Caedibacter acanthamoebae', the endosymbiont of Acanthamoeba sp. strain UWC9 as ' Candidatus Paracaedibacter acanthamoebae' and the endosymbiont of Acanthamoeba sp. strain UWE39 as ' Candidatus Paracaedibacter symbiosus'. The phylogeny of the Acanthamoeba host cells was analysed by comparative sequence analyses of their 18S rRNA. Although Acanthamoeba polyphaga HN-3 clearly groups together with most of the known Acanthamoeba isolates (18S rRNA sequence type 4), Acanthamoeba sp. UWC9 and UWE39 exhibit < 92% 18S rRNA sequence similarity to each other and to other Acanthamoeba isolates. Therefore, we propose two new sequence types (T13 and T14) within the genus Acanthamoeba containing, respectively, Acanthamoeba sp. UWC9 and Acanthamoeba sp. UWE39.     

Algae or Protozoa: Phylogenetic Position of Euglenophytes and Dinoflagellates as Inferred from Mitochondrial Sequences

The chloroplasts of euglenophytes and dinoflagellates have been suggested to be the vestiges of endosymbiotic algae acquired during the process of evolution. However, the evolutionary positions of these organisms are still inconclusive, and they have been tentatively classified as both algae and protozoa. A representative gene of the mitochondrial genome, cytochrome oxidase subunit I (coxI), was chosen and sequenced to clarify the phylogenetic positions of four dinoflagellates, two euglenophytes and one apicomplexan protist. This is the first report of mitochondrial DNA sequences for dinoflagellates and euglenophytes. Our COXI tree shows clearly that dinoflagellates are closely linked to apicomplexan parasites but not with algae. Euglenophytes and algae appear to be only remotely related, with euglenophytes sharing a possible evolutionary link with kinetoplastids. The COXI tree is in general agreement with the tree based on the nuclear encoded small subunit of ribosomal RNA (SSU rRNA) genes, but conflicts with that based on plastid genes. These results support the interpretation that chloroplasts present in euglenophytes and dinoflagellates were captured from algae through endosymbioses, while their mitochondria were inherited from the host cell. We suggest that dinoflagellates and euglenophytes were originally heterotrophic protists and that their chloroplasts are remnants of endosymbiotic algae. Received: 24 March 1997 / Accepted: 21 April 1997     

Phylogenetic position of someChlorella species within the chlorococcales based upon complete small-subunit ribosomal RNA sequences

Summary Complete small-subunit rRNA (16S-like rRNA) coding region sequences were determined for eight species of the Chlorococcales (Chlorophyceae). The genera investigated includePrototheca, Ankistrodesmus, Scenedesmus, and fiveChlorella species. Distance matrix methods were used to infer a phylogenetic tree that describes evolutionary relationships between several plant and green algal groups. The tree exhibits a bifurcation within the Chlorococcales consistent with the division into Oocystaceae and Scenedesmaceae, but three of the fiveChlorella species are more similar to other algae than toChlorella vulgaris. All of the sequences contain primary and secondary structural features that are characteristic of 16S-like rRNAs of chlorophytes and higher plants.Anikstrodesmus stipitatus, however, contains a 394-bp group I intervening sequence in its 16S-like rRNA coding region.     

艾丁嗜盐小盒菌B2菌株16S rRNA核苷酸序列

艾丁嗜盐小盒菌B2菌株(Haloarcula aidinensis, strain B2)16Sr RNA的核苷酸序列已以双脱氧核苷酸链终止法确定。该菌16Sr RNA显示出了典型的古生物类(Archaea)特性。虽然艾丁嗜盐小盒菌B2菌株在序列方面更接近细菌类(Bacteria)的16SrRNA,但它的序列也显示出与真核生物类(Eucarya)的某些特殊的相似性。在序列和结构方面,该菌与细菌类或真核生物类之间的相似程度要高于细菌类与真核生物类之间的相似程度。另外,该菌16SrRNA的序列与其它嗜盐菌序列相比较支持了以前的结论,即艾丁嗜盐小盒菌B2菌株应属于嗜盐小盒菌属(Haloarcula)的一新种。