A phylogenetic analysis of the 18S ribosomal RNA sequence of the coelacanth Latimeria chalumnae

Synopsis Approximately 98% of the sequence of the 18S ribosomal RNA (rRNA) of the coelacanth Latimeria chalumnae was determined by a combination of direct RNA sequencing and sequencing of rRNA genes amplified by the polymerase chain reaction. This sequence was compared with 18S rRNA sequences of similar length from seven other vertebrate species, representing the taxa Petromyzontiformes, Holocephali, Elasmobranchii, Actinopterygii, Dipnoi, Amphibia, and Amniota, in order to determine the most likely sister group of the coelacanth. Maximum parsimony analysis of these sequences resulted in a single most parsimonious tree containing a number of anomalous relationships among these groups. A bootstrap analysis showed that none of the relationships in this tree was significantly supported at the 95% level, however. Addition of data from 15 other vertebrates (providing multiple representatives of most of the higher taxa) resulted in similar ambiguous groupings, as did a number of methods of editing the sites compared (designed to eliminate rapidly evolving positions). These results may be due to a relatively rapid radiation of the major lineages of osteichthyans, the resolution of which will require molecular information from a larger portion of the coelacanth genome.     

18S rRNA alignments derived from different secondary structure models can produce alternative phylogenies

Molecular sequence data are often aligned on the basis of secondary and/or tertiary structure models. However, these models are regularly updated and sometimes differ depending on the way in which they were constructed. We examined whether the choice of a particular 18S rRNA secondary structure model as alignment basis influences phylogeny inference. We therefore compared 18S rRNA phylogenies derived from alignments based on different models. We used: 1. Maximum parsimony; 2. The neighbour-joining method; 3. The maximum-likelihood approach; and 4. Evolutionary parsimony. This demonstrated that the secondary structure model on which an alignment is based may influence: 1. The tree topologies found by these four methods; 2. The numbers of most parsimonious trees found; and 3. The statistical values calculated by the evolutionary parsimony method.     

The phylogenetic status of arthropods, as inferred from 18S rRNA sequences

Partial 18S rRNA sequences of five chelicerate arthropods plus a crustacean, myriapod, insect, chordate, echinoderm, annelid, and platyhelminth were compared. The sequence data were used to infer phylogeny by using a maximum-parsimony method, an evolutionary-distance method, and the evolutionary-parsimony method. The phylogenetic inferences generated by maximum-parsimony and distance methods support both monophyly of the Arthropoda and monophyly of the Chelicerata within the Arthropoda. These results are congruent with phylogenies based on rigorous cladistic analyses of morphological characters. Results support the inclusion of the Arthropoda within a spiralian or protostome coelomate clade that is the sister group of a deuterostome clade, refuting the hypothesis that the arthropods represent the "primitive" sister group of a protostome coelomate clade. Bootstrap analyses and consideration of all trees within 1% of the length of the most parsimonious tree suggest that relationships between the nonchelicerate arthropods and relationships within the chelicerate clade cannot be reliably inferred with the partial 18S rRNA sequence data. With the evolutionary-parsimony method, support for monophyly of the Arthropoda is found in the majority of the combinations analyzed if the coelomates are used as "outgroups." Monophyly of the Chelicerata is supported in most combinations assessed. Our analyses also indicate that the evolutionary-parsimony method, like distance and parsimony, may be biased by taxa with long branches. We suggest that a previous study's inference of the Arthropoda as paraphyletic may be the result of (a) having two few arthropod taxa available for analysis and (b) including long-branched taxa.      

Diversity and evolution of 5S rRNA gene family organization in Pythium

The 5S rRNA gene family organization among 87 species and varieties of Pythium was investigated to assess evolutionary stability of the two patterns detected and to determine which pattern is likely the ancestral state in the genus. Species with filamentous sporangia (Groups A-C according to the ITS phylogenetic tree for Pythium) had 5S genes linked to the rDNA repeat that were predominantly coded for on the DNA strand opposite to the one with the other rRNA genes (‘inverted’ orientation). A small group of species with contiguous sporangia (Group D) is related to Groups A-C but had unlinked 5S genes. The main group of species with spherical zoosporangia (Groups E-J) generally had unlinked 5S genes in tandem arrays. The six species in Group K, although they also have spherical sporangia, had linked genes on the same strand as the other rRNA genes ‘non-inverted’ and most of them had pairs of tandem 5S genes. The evolutionary stability of 5S sequence organization was compared with the stability of morphological characters as interpreted from a phylogeny based on ITS sequence analysis. Features of 5S sequence organization were found to be just as consistent within groups as were the morphological characters. To determine the ancestral type of 5S family organization, a survey of Phytophthora strains was conducted to supply an outgroup reference. The most parsimonious interpretation of the data in this survey yielded the tentative conclusion that the linked condition of the 5S sequences was ancestral.     

The phylogeny of land plants inferred from 18S rDNA sequences: pushing the limits of rDNA signal?

Previous studies of the phylogeny of land plants based on analysis of 18S ribosomal DNA (rDNA) sequences have generally found weak support for the relationships recovered and at least some obviously spurious relationships, resulting in equivocal inferences of land plant phylogeny. We hypothesized that greater sampling of both characters and taxa would improve inferences of land plant phylogeny based on 18S rDNA sequences. We therefore conducted a phylogenetic analysis of complete (or nearly complete) 18S rDNA sequences for 93 species of land plants and 7 green algal relatives. Parsimony analyses with equal weighting of characters and characters state changes and parsimony analyses weighting (1) stem bases half as much as loop bases and (2) transitions half as much as transversions did not produce substantially different topologies. Although the general structure of the shortest trees is consistent with most hypotheses of land plant phylogeny, several relationships, particularly among major groups of land plants, appear spurious. Increased character and taxon sampling did not substantially improve the performance of 18S rDNA in phylogenetic analyses of land plants, nor did analyses designed to accommodate variation in evolutionary rates among sites. The rate and pattern of 18S rDNA evolution across land plants may limit the usefulness of this gene for phylogeny reconstruction at deep levels of plant phylogeny. We conclude that the mosaic structure of 18S rDNA, consisting of highly conserved and highly variable regions, may contain historical signal at two levels. Rapidly evolving regions are informative for relatively recent divergences (e.g., within angiosperms, seed plants, and ferns), but homoplasy at these sites makes it difficult to resolve relationships among these groups. At deeper levels, changes in the highly conserved regions of small-subunit rDNAs provide signal across all of life. Because constraints imposed by the secondary structure of the rRNA may affect the phylogenetic information content of 18S rDNA, we suggest that 18S rDNA sequences be combined with other data and that methods of analysis be employed to accommodate these differences in evolutionary patterns, particularly across deep divergences in the tree of life.     

Convergence in ascospore discharge mechanism among pyrenomycete fungi based on 18S ribosomal RNA gene sequence.

Fungi of the class Pyrenomycetes (Ascomycotina) form a morphological series ranging from those that shoot ascospores (sexual spores) forcibly from the ascus (spore sac) to fungi that ooze ascospores or have no obvious mechanism for ascospore release. Did forcible ascospore discharge evolve within these pyrenomycetes, or has it been lost in the group? We determined the sequences of the 18S ribosomal RNA gene from three fungi and used these, along with six sequences from our previous work and three sequences from GenBank, to infer the phylogeny of 12 ascomycetes with various ascospore discharge mechanisms. The 1720 base pairs of sequence data per fungus yielded 361 variable sites, 198 phylogenetically informative sites, and a single most parsimonious tree requiring 562 nucleotide changes. The tree shows that the capacity to shoot ascospores into the air has been lost or, less probably, gained repeatedly and independently. Species lacking forcible ascospore discharge are intercalated among three lineages of species with forcible discharge. In this tree, seven of the nine internal branches appeared in 95% or more of 500 bootstrap replicates. A tree uniting the fungi with forcible ascospore discharge into a monophyletic group required 45 additional steps and fit significantly less well with the data than the most parsimonious tree, based on a maximum likelihood test. Two of the fungi whose sequence we determined, Pseudallescheria boydii and Sporothrix schenckii, are not closely related to one another, even though both are human pathogens and both are from pyrenomycete lineages lacking forcible ascospore discharge. Using the well-resolved, most parsimonious tree, we inferred base substitution patterns in the 18S rRNA. The transition-to-transversion ratio was 1.9. Of all 12 possible substitutions, 29% were from U to C. At sites corresponding to yeast stem positions, A to G transitions were frequent, perhaps compensating for some of the U to C changes, and maintaining secondary structure base pairing (A to G:U to C = 3:4). In loop or bulge positions without secondary structure base pairing, U to C transitions were still frequent, but A to G transitions were rare (A to G:U to C = 1:5).     

The phylogenetic significance of peptidoglycan types: Molecular analysis of the genera Microbacterium and Aureobacterium based upon sequence comparison of gyrB, rpoB, recA and ppk and 16SrRNA genes

The type strains of 27 species of the genus Microbacterium, family Microbacteriaceae, were analyzed with respect to the phylogeny of the housekeeping genes coding for DNA gyrase subunit B (gyrB), RNA-polymerase subunit B (rpoB), recombinase A (recA) and polyphosphate kinase (ppk). The resulting gene trees were compared to the 16S rRNA gene phylogeny of the same species. The topology of neighbour-joining and maximum parsimony phylogenetic trees based upon nucleic acid sequences and protein sequences of housekeeping genes differed among each other and no gene tree was identical to that of the 16S rRNA gene tree. Only some species showed consistent clustering by all genes analyzed, but the majority of species branched with different neighbours in most gene trees. The failure to phylogenetically cluster type strains into two groups based upon differences in the amino acid composition of peptidoglycan on the basis of 16S rRNA gene sequence similarity, once leading to the union of the genera Microbacterium and Aureobacterium, was also seen in the analysis of recA, rpoB and gyrB gene and protein phylogenies. Analysis of the pkk gene and protein as well as of a concatenate tree, combining sequences of all five genes (total of 3.700 nucleotides), sees members of the former genus Aureobacterium and other type strains with lysine as diagnostic diamino acid to form a coherent cluster that branches within the radiation of Microbacterium species with ornithine in the peptidoglycan.     

Molecular phylogeny of the Gasterosteidae: the importance of using multiple genes

The Gasterosteidae is an important model system in evolutionary biology. Phylogenetic relationships have previously been constructed based upon morphological and behavioral data, but to date no one has investigated those relationships using molecular characters. This paper reports the results of an analysis using sequences from five mitochondrial genes (12S rRNA, 16S rRNA, cytochrome b, ATPase 6, and control region). Phylogenetic analysis of 2879 bp produced a single most parsimonious tree with a consistency index of 72.6%. That tree agrees with the behavior+morphology topology, with one exception: Apeltes is placed as the sister group to Pungitius+Culaea, rather than as the sister-group of (Pungitius+Culaea)+Gasterosteus. This study highlights the importance of using multiple mitochondrial genes in a phylogenetic analysis. Separately, the five genes produced four significantly different topologies, and might have given different versions of gasterosteid relationships had only one or two genes been sequenced. It is thus imperative that comparative biologists choose only trees that contain multiple mitochondrial genes as the basis for studies of evolutionary patterns and processes.     

红豆杉科及相关类群rbcL基因PCR—RFLP分析

运用RFLP方法对红豆杉科及相关类群14种植物叶绿体rbcL基因PCR产物进行限制酶酶切分析,共获29个酶切变异位点。采用PHYLIP软件包对限制位点变异数据进行极大简约法分析得到18个步长为6的最简约树并求得一致树,结果显示：⑴红豆杉科和三尖杉科属单系群;⑵穗花杉属Amentotaxus以置于红豆杉科内为宜,不支持将穗花杉属独立成科的处理方式;⑶白豆杉应为红豆杉科内一个属Pseudotaxus;⑷三尖杉属内篦子三尖杉地位特殊,可设篦子三尖杉组;⑸不赞同将竹柏类从罗汉松属中分离出去成立新科;⑹红豆杉科、三尖杉科和罗汉松科三者间,前两者的关系更为接近。     

红豆杉科及相关类群rbc L基因PCR-RFLP分析

运用RFLP方法对红豆杉科及相关类群14种植物叶绿体rbcL基因PCR产物进行限制酶酶切分析,共获29个酶切变异位点。采用PHYLIP 软件包对限制位点变异数据进行极大简约法分析得到18个步长为6的最简约树并求得一致树,结果显示：（1）红豆杉科和三尖杉科属单系群;(2) 穗花杉属Amentotaxus以置于红豆杉科内为宜,不支持将穗花杉属独立成科的处理方式;(3)白豆杉应为红豆杉科内一个属Pseudotaxus;(4) 三尖杉属内篦子三尖杉地位特殊,可设篦子三尖杉组;(5) 不赞同将竹柏类从罗汉松属中分离出去成立新科;(6) 红豆杉科、三尖杉科和罗汉松科三者间,前两者的关系更为接近。     

Phylogeny and Speciation of Felids

The phylogeny of the Felidae is reconstructed using a total evidence approach combining sequences from 12S rRNA, 16S rRNA, NADH-5, and cytochrome b genes with morphological and karyological characters. The 1504-character data set generated two equally parsimonious trees (CI = 0.413, 1795 steps) of which a strict consensus revealed one polytomy in the placement of the bay cat group. The tree supports several traditional groupings such as the genera Panthera and Lynx and the ocelot group of small South American felids, and it provides additional resolution of relationships within and among the major felid lineages. Combining phylogenetic, distributional, and ecological data indicates that vicariant speciation has played a relatively minor role in the diversification of the felids (approximately 26% of events), while sympatric speciation has been more important than expected on theoretical grounds (approximately 51.8% of events), although postspeciation dispersal may have blurred the boundaries between sympatric, parapatric, and peripheral isolate modes. An examination of ecological changes on the felid tree shows repeated patterns of resource partitioning in time (activity patterns), space (preferred habitat type), and food (as measured by body size) among closely related species. The rapid diversification of the cats thus appears to have been associated more with ecological than with geological separation.     

Investigation of causes of the conflict between taxonomy and molecular phylogeny of trypanosomatids by the example of Leptomonas nabiculae podlipaev, 1987

Results of study of Leptomonas nabiculae using various molecular markers and different material (cultures D2 et Nfm2) contradicted each other and taxonomic position of this species. We investigated morphology of the cells in these cultures as well as in hapantotype of L. nabiculae and those of L. peterhoffi and L. occidentalis that had been described from the same host species. We also reconstructed 18S rRNA gene phylogeny using sequences from both cultures. The D2 culture according to its morphology and phylogenetic position revealed to be a Crithidia that had accompanied L. nabiculae in a mixed infection. We named it Crithidia dedva. The cells in the hapantotypes of the three Leptomonas species and those of the Nfm2 culture represent a single species that is a Herpetomonas (H. nabiculae) judging by morphology and molecular phylogeny. We also showed that the sequence of 18S rRNA gene that had been formerly determined represents a chimaera. This had resulted in the wrong position of this species on the phylogenetic tree that had contradicted results of the analysis of 5s rRNA gene.     

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.     

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 Phylogeny of Fulgoromorpha (Insecta, Hemiptera, Archaeorrhyncha). The Enigmatic Tettigometridae: Evolutionary Affiliations and Historical Biogeography

Previous morphologically based studies by several taxonomists disagree on whether the Tettigometridae represent a sister group to all other fulgoromorphans or whether they are a relatively derived family within the Fulgoromorpha. In this study, several parsimony-based analyses using data sets composed of nucleotide sequences of 18S rDNAs (genes encoding 18S rRNAs) support a monophyletic Fulgoromorpha. All analyses depict Tettigometridae as a relatively derived lineage in a monophyletic relationship with Tropiduchidae. Unscored and unweighted data sets position the tettigometrid + tropiduchid clade as sister to Flatidae. The tettigometrid + tropiduchid clade is supported by four synapomorphic sites, two deletions and two transversions. Constraining tettigometrids to a basal fulgoromorphan lineage significantly reduces parsimony, with several hundreds to thousand of trees being more parsimonious. An analysis employing the Barriel method for scoring potential phylogenetic information in regions containing deletions also supports a non-basal tettigometrid + tropiduchid clade. However, this method results in three equally most parsimonious trees and a sister relationship of the tettigometrid + tropiduchid clade to Flatidae becomes ambiguous.
The molecular-based results showing a derived Tettigometridae agree with previous morphological interpretations of Bourgoin. The current biogeographical distribution of tettigometrids and morphological features supporting the 18S rDNA-based phylogeny are discussed.     

Towards an 18S phylogeny of hexapods: accounting for group-specific character covariance in optimized mixed nucleotide/doublet models

The phylogenetic diversification of Hexapoda is still not fully understood. Morphological and molecular analyses have resulted in partly contradicting hypotheses. In molecular analyses, 18S sequences are the most frequently employed, but it appears that 18S sequences do not contain enough phylogenetic signals to resolve basal relationships of hexapod lineages. Until recently, character interdependence in these data has never been treated seriously, though possibly accounting for the occurrence of biased results. However, software packages are readily available which can incorporate information on character interdependence within a Bayesian approach. Accounting for character covariation derived from a hexapod consensus secondary structure model and applying mixed DNA/RNA substitution models, our Bayesian analysis of 321 hexapod sequences yielded a partly robust tree that depicts many hexapod relationships congruent with morphological considerations. It appears that the application of mixed DNA/RNA models removes many of the anomalies seen in previous studies. We focus on basal hexapod relationships for which unambiguous results are missing. In particular, the strong support for a “Chiastomyaria” clade (Ephemeroptera+Neoptera) obtained in Kjer's [2004. Aligned 18S and insect phylogeny. Syst. Biol. 53, 1–9] study of 18S sequences could not be confirmed by our analysis. The hexapod tree can be rooted with monophyletic Entognatha but not with a clade Ellipura (Collembola+Protura). Compared to previously published contributions, accounting for character interdependence in analyses of rRNA data presents an improvement of phylogenetic resolution. We suggest that an integration of explicit clade-specific rRNA structural refinements is not only possible but an important step in the optimization of substitution models dealing with rRNA data.     

Molecular systematics of Volvocales (Chlorophyceae, Chlorophyta) based on exhaustive 18S rRNA phylogenetic analyses

The taxonomy of Volvocales (Chlorophyceae, Chlorophyta) was traditionally based solely on morphological characteristics. However, because recent molecular phylogeny largely contradicts the traditional subordinal and familial classifications, no classification system has yet been established that describes the subdivision of Volvocales in a manner consistent with the phylogenetic relationships. Towards development of a natural classification system at and above the generic level, identification and sorting of hundreds of sequences based on subjective phylogenetic definitions is a significant step. We constructed an 18S rRNA gene phylogeny based on 449 volvocalean sequences collected using exhaustive BLAST searches of the GenBank database. Many chimeric sequences, which can cause fallacious phylogenetic trees, were detected and excluded during data collection. The results revealed 21 strongly supported primary clades within phylogenetically redefined Volvocales. Phylogenetic classification following PhyloCode was proposed based on the presented 18S rRNA gene phylogeny along with the results of previous combined 18S and 26S rRNA and chloroplast multigene analyses.     

Identification of the most informative regions of the mitochondrial genome for phylogenetic and coalescent analyses

Analysis of complete mitochondrial genome sequences is becoming increasingly common in genetic studies. The availability of full genome datasets enables an analysis of the information content distributed throughout the mitochondrial genome in order to optimize the research design of future evolutionary studies. The goal of our study was to identify informative regions of the human mitochondrial genome using two criteria: (1) accurate reconstruction of a phylogeny and (2) consistent estimates of time to most recent common ancestor (TMRCA). We created two series of datasets by deleting individual genes of varied length and by deleting 10 equal-size fragments throughout the coding region. Phylogenies were statistically compared to the full-coding-region tree, while coalescent methods were used to estimate the TMRCA and associated credible intervals. Individual fragments important for maintaining a phylogeny similar to the full-coding-region tree encompassed bp 577-2122 and 11,399-16,023, including all or part of 12S rRNA, 16S rRNA, ND4, ND5, ND6, and cytb. The control region only tree was the most poorly resolved with the majority of the tree manifest as an unresolved polytomy. Coalescent estimates of TMRCA were less sensitive to removal of any particular fragment(s) than reconstruction of a consistent phylogeny. Overall, we discovered that half the genome, i.e., bp 3669-11,398, could be removed with no significant change in the phylogeny (p(AU)=0.077) while still maintaining overlap of TMRCA 95% credible intervals. Thus, sequencing a contiguous fragment from bp 11,399 through the control region to bp 3668 would create a dataset that optimizes the information necessary for phylogenetic and coalescent analyses and also takes advantage of the wealth of data already available on the control region.     

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.     

Critical analysis of the topology and rooting of the parabasalian 16S rRNA tree

The morphological classification of the protozoan phylum Parabasala is not in absolute agreement with the 16S rRNA phylogeny. However, there are strong indications that tree-construction artifacts play a considerable role in the shaping of the 16S rRNA tree. We have performed rigorous analyses designed to minimize such artifacts using the slow-fast and taxa-exclusion methods. The analyses, which included new sequences from the genera Monocercomonas and Hexamastix, in most respects confirmed the previously suggested tree topology and polyphyly of Hypermastigida and Monocercomonadidae but detected one artificial cluster of long branches (Trichonymphidae, Pseudotrichonymphidae, Hexamastix, and Tricercomitus). They also indicated that the rooting of the phylum on the trichonymphid branch is probably wrong and that reliable rooting on the basis of current data is likely impossible. We discuss the tree topology in the view of anagenesis of cytoskeletal and motility organelles and suggest that a robust taxonomic revision requires extensive analysis of other gene sequences.