Molecular evolution of rDNA in early diverging Metazoa: First comparative analysis and phylogenetic application of complete SSU rRNA secondary structures in Porifera

Background  

The cytoplasmic ribosomal small subunit (SSU, 18S) ribosomal RNA (rRNA) is the most frequently-used gene for molecular phylogenetic studies. However, information regarding its secondary structure is neglected in most phylogenetic analyses. Incorporation of this information is essential in order to apply specific rRNA evolutionary models to overcome the problem of co-evolution of paired sites, which violates the basic assumption of the independent evolution of sites made by most phylogenetic methods. Information about secondary structure also supports the process of aligning rRNA sequences across taxa. Both aspects have been shown to increase the accuracy of phylogenetic reconstructions within various taxa.     

Phylogeny of Molting Protostomes (Ecdysozoa) as Inferred from 18S and 28S rRNA Gene Sequences

Phylogenetic relationships within the group of molting protostomes were reconstructed by comparing the sets of 18S and 28S rRNA gene sequences considered either separately or in combination. The reliability of reconstructions was estimated from the bootstrap indices for major phylogenetic tree nodes and from the degree of congruence of phylogenetic trees obtained by different methods. By either criterion, the phylogenetic trees reconstructed on the basis of both 18 and 28S rRNA gene sequences were better than those based on the 18S or 28S sequences alone. The results of reconstruction are consistent with the phylogenetic hypothesis classifying protostomes into two major clades: molting Ecdysozoa (Priapulida + Kinorhyncha, Nematoda + Nematomorpha, Onychophora + Tardigrada, Myriapoda + Chelicerata, and Crustacea + Hexapoda) and nonmolting Lophotrochozoa (Plathelminthes, Nemertini, Annelida, Mollusca, Echiura, and Sipuncula). Nematomorphs (Nematomorpha) do not belong to the clade Cephalorhyncha (Priapulida + Kinorhyncha). It is concluded that combined data on the 18S and 28S rRNA gene sequences provide a more reliable basis for phylogenetic inferences.__________Translated from Molekulyarnaya Biologiya, Vol. 39, No. 4, 2005, pp. 590–601.Original Russian Text Copyright © 2005 by Petrov, Vladychenskaya.     

Trypanosomes are monophyletic: evidence from genes for glyceraldehyde phosphate dehydrogenase and small subunit ribosomal RNA

The genomes of Trypanosoma brucei, Trypanosoma cruzi and Leishmania major have been sequenced, but the phylogenetic relationships of these three protozoa remain uncertain. We have constructed trypanosomatid phylogenies based on genes for glycosomal glyceraldehyde phosphate dehydrogenase (gGAPDH) and small subunit ribosomal RNA (SSU rRNA). Trees based on gGAPDH nucleotide and amino acid sequences (51 taxa) robustly support monophyly of genus Trypanosoma, which is revealed to be a relatively late-evolving lineage of the family Trypanosomatidae. Other trypanosomatids, including genus Leishmania, branch paraphyletically at the base of the trypanosome clade. On the other hand, analysis of the SSU rRNA gene data produced equivocal results, as trees either robustly support or reject monophyly depending on the range of taxa included in the alignment. We conclude that the SSU rRNA gene is not a reliable marker for inferring deep level trypanosome phylogeny. The gGAPDH results support the hypothesis that trypanosomes evolved from an ancestral insect parasite, which adapted to a vertebrate/insect transmission cycle. This implies that the switch from terrestrial insect to aquatic leech vectors for fish and some amphibian trypanosomes was secondary. We conclude that the three sequenced pathogens, T. brucei, T. cruzi and L. major, are only distantly related and have distinct evolutionary histories.     

The evolutionary position of Lestoniidae revealed by molecular autapomorphies in the secondary structure of rRNA besides phylogenetic reconstruction (Insecta: Hemiptera: Heteroptera)

Pentatomoidea (stink bugs and their relatives) is the third largest superfamily in Heteroptera, or the true bugs. The phylogenetic relationships among the families within Pentatomoidea remain controversial. The family Lestoniidae is morphologically highly specialized, currently including only two species endemic to Australia. Previous researchers have suggested a close relationship of Lestoniidae to either Plataspidae or Acanthosomatidae, based on morphological characters. In this study, phylogenetic tree reconstruction revealed that Lestoniidae and Acanthosomatidae form a monophyletic clade. In addition, in comparisons of the secondary structures of 18S and 28S rRNAs representing 15 families of Pentatomoidea, four length‐variable regions in 18S and 28S rRNAs that can serve as autapomorphies for the clade Lestoniidae + Acanthosomatidae were recognized. Among them, E in 18S rRNA and D3‐1 and D5‐1 in 28S rRNA are unique in length in Lestoniidae and Acanthosomatidae. Based on the new molecular evidence and morphological evidence published by previous authors, Lestoniidae is suggested to be a highly specialized group derived from a common ancestor with Acanthosomatidae.     

Molecular phylogeny of Collembola inferred from ribosomal RNA genes

The classification of taxa within Collembola (Springtails, Hexapoda) has been controversial. In this study, we combined complete 18S rRNA gene with partial 28S rRNA gene (D7-D10) sequences to investigate the phylogeny of Collembola. About 2500 aligned sites of thirty species representing 29 genera from14 families of Collembola were analyzed, including one species of Neelipleona from which no sequence has been reported previously.The phylogenetic trees were obtained by different methods (maximum parsimony, maximum likelihood, and Bayesian analysis). Our results supported the monophyly of two of the four taxonomic groups of Collembola summarized by Deharveng [Deharveng, L., 2004. Recent advances in Collembola systematics. Pedobiologia 48, 415–433.], namely of Poduromorpha and of Symphypleona. Within Poduromorpha, Neanuridae was monophyletic with high support, but Hypogastruridae was not. Entomobryomorpha was paraphyletic, as the Tomoceroidea (Tomoceridae and Oncopoduridae) was found to be apart from the other entomobryomorphs. In the latter Isotomoidea and Entomobryoidea joined into a group with moderate support. Within Symphypleona, the phylogenetic relationship [(Sminthuridae + Bourletiellidae) + Sminthurididae] was consistent with traditional morphological studies. Neelipleona grouped with Symphypleona in all trees, with moderate support in the ML and Bayesian analyses.     

Vicariance of Pyrocoelia fireflies (Coleoptera: Lampyridae) in the Ryukyu islands,Japan

We performed molecular phylogenetic analyses based on the mitochondrial COI gene (687 bp) and the nuclear 28S rRNA gene (715 bp) and reconstructed phylogenetic trees of the Pyrocoelia fireflies in the Ryukyu Islands and eastern Asia. Age calibration was done using a robust geological constraint: the Okinawa trough and associated straits began to rift at 1.55 Ma, isolating the Ryukyu Islands from the Chinese continent, Japanese islands, Taiwan island and some of the islands from each other. We suggest that the physical isolation of these islands began to generate the allopatric speciation within these islands, so the timing of this isolation was assigned to an appropriate node. The topology is completely concordant among phylogenetic trees reconstructed using MEGA (maximum‐likelihood), raxmlGUI (maximum‐likelihood) and BEAST (Bayesian inference; including combined analysis of COI and 28S rRNA genes). Two lineages are recognized, related to their emergence time; spring to summer, and autumn. In each lineage, vicariance is inferred to have begun at 1.55 Ma from our phylogenetic and geological analyses. In lineage 1, P. oshimana (Amami), P. matsumurai (Okinawa), P. discicollis (W. Japan), P. fumosa (E. Japan) and P. abdominalis (Yaeyama) were differentiated. In lineage 2, P. rufa (Tsushima and Korea), P. miyako (Miyako‐jima), P. atripennis (Ishigaki‐jima) and P. praetexta (Taiwan and HongKong) were differentiated. Pyrocoelia analis (Taiwan and China) emerges throughout the year except for winter, and constitutes another lineage. We suggest that Pyrocoelia fireflies differentiated at 2 Ma to generate these three lineages. The base substitution rate for the COI gene is estimated as 4.48% Myr^–1 and that for the 28S rRNA gene is 0.394%  Myr^–1, and these rates were used in a combined BEAty analysis in BEAST.     

Morphology and Phylogenetic Analyses of Three Novel Naegleria Isolated from Freshwaters on Jeju Island,Korea, During the Winter Period

The genus Naegleria is one of the best known heterolobosean groups, and is the causative agent of primary amoebic meningoencephalitis. This group is rarely studied in temperate regions during winter. Here, three novel Naegleria were isolated from freshwaters on Jeju Island, Korea, during winter. Two isolates were amoeboflagellates, and one of the three amoebae did not undergo enflagellation. All amoebae had eruptive pseudopodia, and the layer of refractile granules around a large nucleus. They formed a cyst with ~2 pores in the cyst stage. The amoeboflagellate form had two flagella and no division in the flagellate stage, and no cytostome. These features are very similar to typical Naegleria. Furthermore, our isolates were able to grow at > 30 °C, suggesting that they had different thermophilicity from Naegleria in polar regions. All amoebae were largely encysted at 5 or 10 °C, indicating that they were likely encysted during winter. Based on the 18S rRNA gene and the ITS1‐5.8S rRNA gene‐ITS2 sequences, the phylogenetic analyses consistently revealed that the isolates are members of the Naegleria group. However, the isolates differ from other species in both phylogenetic trees. Thus, Naegleria in cold habitats appeared to have a high degree of novelty, but their thermophilicity may be dependent on locality.     

Molecular phylogenetics of Trypanosomatidae: contrasting results from 18S rRNA and protein phylogenies

Phylogenetic analyses of the family Trypanosomatidae have been conducted using both 18S rRNA gene sequences and a variety of protein sequences. Using a variety of phylogenetic methods, 18S rRNA phylogenies indicate that the genus Trypanosoma is not monophyletic. Rather, they suggest that the American and African trypanosomes constitute distinct clades. By contrast, phylogenetic analyses of available sequences in 42 protein families gene generally supported monophyly of the genus Trypanosoma. One possible explanation for these conflicting results is poor taxon sampling in the case of protein coding genes, most of which have been sequenced for only a few species of Trypanosomatidae.     

RNA SEQUENCE DATA IN PHYLOGENETIC RECONSTRUCTION: TESTING THE LIMITS OF ITS RESOLUTION

Abstract— 18S ribosomal RNA sequences from 11 echinoderms are analysed using parsimony to investigate phylogenetic relationships. Their estimated divergence limes range from less than 20 Ma to more than 550 Ma before present. Phylogenies based on 18S rRNA sequence data are compared with well-established morphological phylogenies to discover at what evolutionary distance the two approaches start to produce incongruent results. Three regions of the 18S rRNA molecule are analysed separately and together, and paired and unpaired sites are also treated separately and combined.
Results show that a parsimony analysis of sequence data produces reliable results only when taxa have diverged more recently than about 100 Ma. At greater evolutionary distances (up to 250 Ma), paired nucleotides produce more reliable results than unpaired, while paired and unpaired data combined produce intermediate results. All trees within about 1% of the most parsimonious solution ought to be accepted. Transversions give results almost as reliable as paired regions though there were relatively few informative sites. The relationships of echinoderm classes, which diverged 450–550 Ma ago, are unresolved by 18S rRNA data.     

Phylogeny of protostome moulting animals (Ecdysozoa) inferred from 18 and 28S rRNA gene sequences

Reliability of reconstruction of phylogenetic relationships within a group of protostome moulting animals was evaluated by means of comparison of 18 and 28S rRNA gene sequences sets both taken separately and combined. Reliability of reconstructions was evaluated by values of the bootstrap support of major phylogenetic tree nodes and by degree of congruence of phylogenetic trees inferred by various methods. By both criteria, phylogenetic trees reconstructed from the combined 18 and 28S rRNA gene sequences were better than those inferred from 18 and 28S sequences taken separately. Results obtained are consistent with phylogenetic hypothesis separating protostome animals into two major clades, moulting Ecdysozoa (Priapulida + Kinorhyncha, Nematoda + Nematomorpha, Onychophora + Tardigrada, Myriapoda + Chelicerata, Crustacea + Hexapoda) and unmoulting Lophotrochozoa (Plathelminthes, Nemertini, Annelida, Mollusca, Echiura, Sipuncula). Clade Cephalorhyncha does not include nematomorphs (Nematomorpha). Conclusion was taken that it is necessary to use combined 18 and 28S data in phylogenetic studies.     

The complete mitochondrial genome of Japanese sparrowhawk (Accipiter gularis) and the phylogenetic relationships among some predatory birds

Phylogenetic relationships among raptors, especially various groups are rather complex and controversial. We determined the complete mtDNA of Japanese sparrowhawk, and estimated phylogenetic trees based on the complete mtDNA alignment of it and 36 other raptor species, to clarify raptor phylogenetics. Phylogenetic trees were also estimated using a multiple sequence alignment of 12S rRNA and 16S rRNA from 81 typical species in GenBank, to further clarify the phylogenetic relationships of several groups among the raptors. The new mtDNA is a circular molecule, 17 917 bp in length, containing the 37 typical genes, with a pseudo-control region. ATG is generally the start codon, TAA is the most frequent stop codon. All tRNAs can be folded into canonical cloverleaf secondary structures except for tRNA^{Ser (AGY)} and tRNA^{Leu (CUN)}, which are missing the “DHU” arm. Phylogenetic relationships demonstrate that raptors can be divided into four branches: Accipitriformes, Falconiformes, Strigiformes and Caprimulgiformes in this study. We suggest that Accipitriformes should to be an independent order, Accipitriformes. The results also indicate that Accipitriformes contains three clades: Accipitridae, Pandionidae and Sagittariidae. Strigiformes includes species from Tytonidae and Strigidae. Caprimulgiformes contains Aegothelidae and Caprimulgidae.     

Introducing W.A.T.E.R.S.: a Workflow for the Alignment,Taxonomy, and Ecology of Ribosomal Sequences

Background  

For more than two decades microbiologists have used a highly conserved microbial gene as a phylogenetic marker for bacteria and archaea. The small-subunit ribosomal RNA gene, also known as 16 S rRNA, is encoded by ribosomal DNA, 16 S rDNA, and has provided a powerful comparative tool to microbial ecologists. Over time, the microbial ecology field has matured from small-scale studies in a select number of environments to massive collections of sequence data that are paired with dozens of corresponding collection variables. As the complexity of data and tool sets have grown, the need for flexible automation and maintenance of the core processes of 16 S rDNA sequence analysis has increased correspondingly.     

Morphological and molecular characterization of the name-bearing type species Rimaleptus binucleatus (Kahl, 1931), with a phylogenetic re-analysis of dileptid evolutionary history (Ciliophora: Litostomatea: Rhynchostomatia)

We collected Rimaleptus binucleatus from soil in the surroundings of the town of Ulsan, South Korea. Its morphology and 18S rRNA gene were studied using standard methods. This possibly widely distributed species is characterized by: (i) a size of about 170–400 × 20–65 μm; (ii) a narrowly to cylindrically dileptid body with proboscis occupying about 30–65% of body length; (iii) two dorsal contractile vacuoles; (iv) two size groups of rod-shaped extrusomes; and (v) about 18–29 ciliary rows, 4–6 of them anteriorly differentiated into a staggered dorsal brush. Phylogenetic analyses of five new rhynchostomatian 18S rRNA gene sequences supported monophylies of the orders Tracheliida and Dileptida, but revealed that the genera Rimaleptus and Pseudomonilicaryon are polyphyletic. Monophyly of genera with two macronuclear nodules was consistently rejected, but monophylies of dileptids with many scattered macronuclear nodules and of dileptids with moniliform macronucleus and multi-rowed dorsal brush could not be excluded by statistical topology tests. Nevertheless, phylogenetic network analyses indicated considerable conflict in the phylogenetic signal provided by the 18S rRNA gene to resolve unambiguously relationships among dileptid genera.     

Application of <Emphasis Type="Italic">gyrB</Emphasis> and <Emphasis Type="Italic">parE</Emphasis> sequence similarity analyses for differentiation of species within the genus <Emphasis Type="Italic">Geobacillus</Emphasis>

The primary structures of the genes encoding the β-subunits of a type II topoisomerase (gyrase, gyrB) and a type IV topoisomerase (parE) were determined for 15 strains of thermophilic bacteria of the genus Geobacillus. The obtained sequences were used for analysis of the phylogenetic similarity between members of this genus. Comparison of the phylogenetic trees of geobacilli constructed on the basis of the 16S rRNA, gyrB, and parE gene sequences demonstrated that the level of genetic distance between the sequences of the genes encoding the β-subunits of type II topoisomerases significantly exceeded the values obtained by comparative analysis of the 16S rRNA gene sequences of Geobacillus strains. It was shown that, unlike the 16S rRNA gene analysis, comparative analysis of the gyrB and parE gene sequences provided a more precise determination of the phylogenetic position of bacteria at the species level. The data obtained suggest the possibility of using the genes encoding the β-subunits of type II topoisomerases as phylogenetic markers for determination of the species structure of geobacilli.     

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.     

Morphological and molecular description of Blastocrithidia cyrtomeni sp. nov. (Kinetoplastea: Trypanosomatidae) associated with Cyrtomenus bergi Froeschner (Hemiptera: Cydnidae) from Colombia

A new trypanosomatid species, Blastocrithidia cyrtomeni, is herein described using morphological and molecular data. It was found parasitising the alimentary tract of the insect host Cyrtomenus bergi, a polyphagous pest. The morphology of B. cyrtomeni was investigated using light and transmission microscopy and molecular phylogeny was inferred from the sequences of spliced leader RNA (SL rRNA) - 5S rRNA gene repeats and the 18S small subunit (SSU) rRNA gene. Epimastigotes of variable size with straphanger cysts adhering to the middle of the flagellum were observed in the intestinal tract, hemolymph and Malpighian tubules. Kinetoplasts were always observed anterior to the nucleus. The ultrastructure of longitudinal sections of epimastigotes showed the flagellum arising laterally from a relatively shallow flagellar pocket near the kinetoplast. SL RNA and 5S rRNA gene repeats were positive in all cases, producing a 0.8-kb band. The amplicons were 797-803 bp long with > 98.5% identity, indicating that they originated from the same organism. According to the sequence analysis of the SL-5S rRNA gene repeats and the 18S SSU rRNA gene, B. cyrtomeni is different from all other known species or isolates of Trypanosomatidae. Both analyses indicate that among known species, it is most closely related to Blastocrithidia triatomae.     

Transcription and processing map of the 4.5S-5S rRNA intergenic region (ITS3) from rapeseed (Brassica napus) chloroplasts

Brassica napus   can be applied to other known plant sequences. The conservation of those sequences suggests a functional role for them possibly by providing recognition structures for endogenous RNases involved in the 4.5S rRNA-5S rRNA maturing process. Most of the processing signals detected here are located at single-stranded regions of a proposed secondary structure. Received: 1 June 1999 / Revision received: 14 December 1999 / Accepted: 18 December 1999     

Identification,genetic variation,and structural analysis of 18S rRNA of Theileria orientalis and Theileria velifera-like isolates from Myanmar

Ribosomal RNA genes have been widely used for the identification and phylogenetic analysis of various organisms, including parasitic protozoa. Here, we report nine near full-length Theileria orientalis 18S rRNA gene sequences from cattle from different areas of Myanmar. Phylogenetic analysis of the 18S rRNA genes revealed a considerably close genetic relationship among T. orientalis isolates from Australia, China, Japan, Korea, Myanmar, and Pakistan. We also obtained four Theileria velifera-like (Theileria cf. velifera) 18S rRNA gene sequences from two cattle and two water buffaloes from the northernmost area of Myanmar. The phylogenetic analysis of T. cf. velifera isolates from Myanmar along with T. velifera and T. cf. velifera isolates from African countries suggested an evolutionary lineage of greater complexity in T. velifera-related parasites. DNA alignment analysis indicated the presence of 51 and 55 nucleotide variation positions within the 18S rRNA genes from 15 T. orientalis and 11 T. velifera-related isolates, respectively. Alignment entropy analysis of the 18S rRNA sequences indicated that both T. orientalis and T. velifera-related isolates had three hyper variable regions, corresponding to V2, V4, and V7 regions in eukaryotes. The degree of variation was prominent in the V2 in T. orientalis and V4 in T. velifera-related isolates. The secondary structure analysis of the 18S rRNA predicted using minimum free energy algorism revealed that the structure of V4 region differed most significantly between T. orientalis and T. velifera. These results provide novel insights into common structures, variations and functions of small subunit rRNA in Theileria species.     

Evaluation of sequence alignments and oligonucleotide probes with respect to three-dimensional structure of ribosomal RNA using ARB software package

Background  

Availability of high-resolution RNA crystal structures for the 30S and 50S ribosomal subunits and the subsequent validation of comparative secondary structure models have prompted the biologists to use three-dimensional structure of ribosomal RNA (rRNA) for evaluating sequence alignments of rRNA genes. Furthermore, the secondary and tertiary structural features of rRNA are highly useful and successfully employed in designing rRNA targeted oligonucleotide probes intended for in situ hybridization experiments. RNA3D, a program to combine sequence alignment information with three-dimensional structure of rRNA was developed. Integration into ARB software package, which is used extensively by the scientific community for phylogenetic analysis and molecular probe designing, has substantially extended the functionality of ARB software suite with 3D environment.