Fast statistical tests for detecting heterotachy in protein evolution

The w statistic introduced by Lockhart et al. (1998. A covariotide model explains apparent phylogenetic structure of oxygenic photosynthetic lineages. Mol Biol Evol. 15:1183-1188) is a simple and easily calculated statistic intended to detect heterotachy by comparing amino acid substitution patterns between two monophyletic groups of protein sequences. It is defined as the difference between the fraction of varied sites in both groups and the fraction of varied sites in each group. The w test has been used to distinguish a covarion process from equal rates and rates variation across sites processes. Using simulation we show that the w test is effective for small data sets and for data sets that have low substitution rates in the groups but can have difficulties when these conditions are not met. Using site entropy as a measure of variability of a sequence site, we modify the w statistic to a w' statistic by assigning as varied in one group those sites that are actually varied in both groups but have a large entropy difference. We show that the w' test has more power to detect two kinds of heterotachy processes (covarion and bivariate rate shifts) in large and variable data. We also show that a test of Pearson's correlation of the site entropies between two monophyletic groups can be used to detect heterotachy and has more power than the w' test. Furthermore, we demonstrate that there are settings where the correlation test as well as w and w' tests do not detect heterotachy signals in data simulated under a branch length mixture model. In such cases, it is sometimes possible to detect heterotachy through subselection of appropriate taxa. Finally, we discuss the abilities of the three statistical tests to detect a fourth mode of heterotachy: lineage-specific changes in proportion of variable sites.     

Evolutionary Relationships Among the Eukaryotic Crown Taxa Taking into Account Site-to-Site Rate Variation in 18S rRNA

In this study we constructed a bootstrapped distance tree of 500 small subunit ribosomal RNA sequences from organisms belonging to the so-called crown of eukaryote evolution. Taking into account the substitution rate of the individual nucleotides of the rRNA sequence alignment, our results suggest that (1) animals, true fungi, and choanoflagellates share a common origin: The branch joining these taxa is highly supported by bootstrap analysis (bootstrap support [BS] > 90%), (2) stramenopiles and alveolates are sister groups (BS = 75%), (3) within the alveolates, dinoflagellates and apicomplexans share a common ancestor BS > 95%), while in turn they both share a common origin with the ciliates (BS > 80%), and (4) within the stramenopiles, heterokont algae, hyphochytriomycetes, and oomycetes form a monophyletic grouping well supported by bootstrap analysis (BS > 85%), preceded by the well-supported successive divergence of labyrinthulomycetes and bicosoecids. On the other hand, many evolutionary relationships between crown taxa are still obscure on the basis of 18S rRNA. The branching order between the animal-fungal-choanoflagellates clade and the chlorobionts, the alveolates and stramenopiles, red algae, and several smaller groups of organisms remains largely unresolved. When among-site rate variation is not considered, the inferred tree topologies are inferior to those where the substitution rate spectrum for the 18S rRNA is taken into account. This is primarily indicated by the erroneous branching of fast-evolving sequences. Moreover, when different substitution rates among sites are not considered, the animals no longer appear as a monophyletic grouping in most distance trees. Received: 11 June 1997 / Accepted: 21 July 1997     

Heterotachy, an important process of protein evolution.

Because of functional constraints, substitution rates vary among the positions of a protein but are usually assumed to be constant at a given site during evolution. The distribution of the rates across the sequence positions generally fits a Gamma distribution. Models of sequence evolution were accordingly designed and led to improved phylogenetic reconstruction. However, it has been convincingly demonstrated that the evolutionary rate of a given position is not always constant throughout time. We called such within-site rate variations heterotachy (for "different speed" in Greek). Yet, heterotachy was found among homologous sequences of distantly related organisms, often with different functions. In such cases, the functional constraints are likely different, which would explain the different distribution of variable sites. To evaluate the importance of heterotachy, we focused on amino acid sequences of mitochondrial cytochrome b, for which the function is likely the same in all vertebrates. Using 2,038 sequences, we demonstrate that 95% of the variable positions are heterotachous, i.e., underwent dramatic variations of substitution rate among vertebrate lineages. Heterotachy even occurs at small evolutionary scale, and in these cases it is very unlikely to be related to functional changes. Since a large number of sequences are required to efficiently detect heterotachy, the extent of this phenomenon could not be estimated for all proteins yet. It could be as large as for cytochrome b, since this protein is not a peculiar case. The observations made here open several new avenues of research, such as the understanding of the evolution of functional constraints or the improvement of phylogenetic reconstruction methods.     

Artifactual phylogenies caused by correlated distribution of substitution rates among sites and lineages: the good, the bad, and the ugly

Despite the advances in understanding molecular evolution, current phylogenetic methods barely take account of a fraction of the complexity of evolution. We are chiefly constrained by our incomplete knowledge of molecular evolutionary processes and the limits of computational power. These limitations lead to the establishment of either biologically simplistic models that rarely account for a fraction of the complexity involved or overfitting models that add little resolution to the problem. Such oversimplified models may lead us to assign high confidence to an incorrect tree (inconsistency). Rate-across-site (RAS) models are commonly used evolutionary models in phylogenetic studies. These account for heterogeneity in the evolutionary rates among sites but do not account for changing within-site rates across lineages (heterotachy). If heterotachy is common, using RAS models may lead to systematic errors in tree inference. In this work we show possible misleading effects in tree inference when the assumption of constant within-site rates across lineages is violated using maximum likelihood. Using a simulation study, we explore the ways in which gamma stationary models can lead to wrong topology or to deceptive bootstrap support values when the within-site rates change across lineages. More precisely, we show that different degrees of heterotachy mislead phylogenetic inference when the model assumed is stationary. Finally, we propose a geometry-based approach to visualize and to test for the possible existence of bias due to heterotachy.     

Combined mitochondrial and nuclear sequences support the monophyly of forcipulatacean sea stars

Previous molecular phylogenetic analyses of forcipulatacean sea stars (Echinodermata: Asteroidea) have reconstructed a non-monophyletic order Forcipulatida, provided that two or more forcipulate families are included. This result could mean that one or more assumptions of the reconstruction method was violated, or else the traditional classification could be erroneous. The present molecular phylogenetic analysis included 12 non-forcipulatacean and 39 forcipulatacean sea stars, with multiple representatives of all but one of the forcipulate families and/or subfamilies. Bayesian analysis of approximately 4.2kb of sequence data representing seven partitions (nuclear 18S rRNA and 28S rRNA, mitochondrial 12S rRNA, 16S rRNA, 5 tRNAs and cytochrome oxidase I with first and second codon positions analyzed separately from third codon positions) recovered a consensus tree with three well-supported clades (78%-100% bootstrap support) that corresponded at least approximately to traditional taxonomic ranks: the superorder Forcipulatacea (Forcipulatida + Brisingida) + Pteraster, the Brisingida/Brisingidae and Asteriidae + Rathbunaster + Pycnopodia. When a molecular clock was enforced, the partitioned Bayesian analysis recovered the traditional Forcipulatacea. Five of six genera represented by two or more species were monophyletic with 100% bootstrap support. Most of the traditional subfamilial and familial groupings within the Forcipulatida were either unresolved or non-monophyletic. The separate partitions differed considerably in estimates of model parameters, mainly between nuclear sequences (with high GC content, low rates of sequence substitution and high transition/transversion rate ratios) and mitochondrial sequences.     

Ribosomal DNA and resolution of branching order among the ascomycota: how many nucleotides are enough?

Molecular phylogenies for the fungi in the Ascomycota rely heavily on 18S rRNA gene sequences but this gene alone does not answer all questions about relationships. Particularly problematical are the relationships among the first ascomycetes to diverge, the Archiascomycetes, and the branching order among the basal filamentous ascomycetes, the Euascomycetes. Would more data resolve branching order? We used the jackknife and bootstrapping resampling approach that constitutes the "pattern of resolved nodes" method to address the relationship between number of variable sites in a DNA sequence alignment and support for taxonomic clusters. We graphed the effect of increasing sizes of subsamples of the 18S rRNA gene sequences on bootstrap support for nodes in the Ascomycota tree. Nodes responded differently to increasing data. Some nodes, those uniting the filamentous ascomycetes for example, would still have been well supported with only two thirds of the 18S rRNA gene. Other nodes, like the one uniting the Archiascomycetes as a monophyletic group, would require about double the number of variable sites available in the 18S gene for 95% neighbor-joining bootstrap support. Of the several groups emerging at the base of the filamentous ascomycetes, the Pezizales receive the most support as the first to diverge. Our analysis suggests that we would also need almost three times as much sequence data as that provided by the 18S gene to confirm the basal position for the Pezizales and more than seven times as much data to resolve the next group to diverge. If more data from other genes show the same pattern, the lack of resolution for the filamentous ascomycetes may indicate rapid radiation within this clade.     

Exploring among-site rate variation models in a maximum likelihood framework using empirical data: effects of model assumptions on estimates of topology,branch lengths,and bootstrap support

We have investigated the effects of different among-site rate variation models on the estimation of substitution model parameters, branch lengths, topology, and bootstrap proportions under minimum evolution (ME) and maximum likelihood (ML). Specifically, we examined equal rates, invariable sites, gamma-distributed rates, and site-specific rates (SSR) models, using mitochondrial DNA sequence data from three protein-coding genes and one tRNA gene from species of the New Zealand cicada genus Maoricicada. Estimates of topology were relatively insensitive to the substitution model used; however, estimates of bootstrap support, branch lengths, and R-matrices (underlying relative substitution rate matrix) were strongly influenced by the assumptions of the substitution model. We identified one situation where ME and ML tree building became inaccurate when implemented with an inappropriate among-site rate variation model. Despite the fact the SSR models often have a better fit to the data than do invariable sites and gamma rates models, SSR models have some serious weaknesses. First, SSR rate parameters are not comparable across data sets, unlike the proportion of invariable sites or the alpha shape parameter of the gamma distribution. Second, the extreme among-site rate variation within codon positions is problematic for SSR models, which explicitly assume rate homogeneity within each rate class. Third, the SSR models appear to give severe underestimates of R-matrices and branch lengths relative to invariable sites and gamma rates models in this example. We recommend performing phylogenetic analyses under a range of substitution models to test the effects of model assumptions not only on estimates of topology but also on estimates of branch length and nodal support.     

Heterotachy and functional shift in protein evolution

Study of structure/function relationships constitutes an important field of research, especially for modification of protein function and drug design. However, the fact that rational design (i.e. the modification of amino acid sequences by means of directed mutagenesis, based on knowledge of the three-dimensional structure) appears to be much less efficient than irrational design (i.e. random mutagenesis followed by in vitro selection) clearly indicates that we understand little about the relationships between primary sequence, three-dimensional structure and function. The use of evolutionary approaches and concepts will bring insights to this difficult question. The increasing availability of multigene family sequences that has resulted from genome projects has inspired the creation of novel in silico evolutionary methods to predict details of protein function in duplicated (paralogous) proteins. The underlying principle of all such approaches is to compare the evolutionary properties of homologous sequence positions in paralogs. It has been proposed that the positions that show switches in substitution rate over time--i.e., 'heterotachous sites'--are good indicators of functional divergence. However, it appears that heterotachy is a much more general process, since most variable sites of homologous proteins with no evidence of functional shift are heterotachous. Similarly, it appears that switches in substitution rate are as frequent when paralogous sequences are compared as when orthologous sequences are compared. Heterotachy, instead of being indicative of functional shift, may more generally reflect a less specific process related to the many intra- and inter-molecular interactions compatible with a range of more or less equally viable protein conformations. These interactions will lead to different constraints on the nature of the primary sequences, consistently with theories suggesting the non-independence of substitutions in proteins. However, a specific type of amino acid variation might constitute a good indicator of functional divergence: substitutions occurring at positions that are generally slowly evolving. Such substitutions at constrained sites are indeed much more frequent soon after gene duplication. The identification and analysis of these sites by complementing structural information with evolutionary data may represent a promising direction to future studies dealing with the functional characterization of an ever increasing number of multi-gene families identified by complete genome analysis.     

Noncoding sequences from the slowly evolving chloroplast inverted repeat in addition to rbcL data do not support gnetalean affinities of angiosperms

We developed PCR primers against highly conserved regions of the rRNA operon located within the inverted repeat of the chloroplast genome and used these to amplify the region spanning from the 3' terminus of the 23S rRNA gene to the 5' terminus of the 5S rRNA gene. The sequence of this roughly 500-bp region, which includes the 4.5S rRNA gene and two chloroplast intergenic transcribed spacer regions (cpITS2 and cpITS3), was determined from 20 angiosperms, 7 gymnosperms, and 16 ferns (21,700 bp). Sequences for the large subunit of ribulose bisphosphate carboxylase/oxygenase (rbcL) from the same or confamilial genera were analyzed in both separate and combined data sets. Due to the low substitution rate in the inverted repeat region, noncoding sequences in the cpITS region are not saturated with substitutions, in contrast to synonymous sites in rbcL, which are shown to evolve roughly six times faster than noncoding cpITS sequences. Several length polymorphisms with very clear phylogenetic distributions were detected in the data set. Results of phylogenetic analyses provide very strong bootstrap support for monophyly of both spermatophytes and angiosperms. No support for a sister group relationship between Gnetales and angiosperms in either cpITS or rbcL data was found. Rather, weak bootstrap support for monophyly of gymnosperms studied and for a basal position for the aquatic angiosperm Nymphaea among angiosperms studied was observed. Noncoding sequences from the inverted repeat region of chloroplast DNA appear suitable for study of land plant evolution.      

Functional divergence prediction from evolutionary analysis: a case study of vertebrate hemoglobin

It is a central assumption of evolution that gene duplications provide the genetic raw material from which to create proteins with new functions. The increasing availability in multigene family sequences that has resulted from genome projects has inspired the creation of novel in silico approaches to predict details of protein function. The underlying principle of all such approaches is to compare the evolutionary properties of homologous sequence positions in paralogous proteins. It has been proposed that the positions that show switches in substitution rate over time-i.e., "heterotachous sites," are good indicators of functional divergence. Here, we analyzed the alpha and beta paralogous subunits of hemoglobin in search for such signatures. We found as many heterotachous sites in comparisons between groups of paralogous subunits (alpha/beta) as between orthologous ones (alpha/alpha, beta/beta). Thus, the importance of substitution rate shifts as predictors of specialization between protein subfamilies might be reconsidered. Instead, such shifts may reflect a more general process of protein evolution, consistent with the fact that they can be compatible with function conservation. As an alternative, we focused on those residues showing highly constrained states in two sequence groups, but different in each group, and we named them CBD (for "constant but different"). As opposed to heterotachous positions, CBD sites were markedly overrepresented in paralogous (alpha/beta) comparisons, as opposed to orthologous ones (alpha/alpha, beta/beta), identifying them as likely signatures of functional specialization between the two subunits. When superimposed onto the three-dimensional structure of hemoglobin, CBD positions consistently appeared to cluster preferentially on inter-subunit surfaces, two contact areas crucial to function in vertebrate tetrameric hemoglobin. The identification and analysis of CBD sites by complementing structural information with evolutionary data may represent a promising direction for future studies dealing with the functional characterization of a growing number of multigene families identified by complete genome analyses.     

Investigation of molluscan phylogeny using large-subunit and small-subunit nuclear rRNA sequences

The Mollusca represent one of the most morphologically diverse animal phyla, prompting a variety of hypotheses on relationships between the major lineages within the phylum based upon morphological, developmental, and paleontological data. Analyses of small-ribosomal RNA (SSU rRNA) gene sequence have provided limited resolution of higher-level relationships within the Mollusca. Recent analyses suggest large-subunit (LSU) rRNA gene sequences are useful in resolving deep-level metazoan relationships, particularly when combined with SSU sequence. To this end, LSU (approximately 3.5 kb in length) and SSU (approximately 2 kb) sequences were collected for 33 taxa representing the major lineages within the Mollusca to improve resolution of intraphyletic relationships. Although the LSU and combined LSU+SSU datasets appear to hold potential for resolving branching order within the recognized molluscan classes, low bootstrap support was found for relationships between the major lineages within the Mollusca. LSU+SSU sequences also showed significant levels of rate heterogeneity between molluscan lineages. The Polyplacophora, Gastropoda, and Cephalopoda were each recovered as monophyletic clades with the LSU+SSU dataset. While the Bivalvia were not recovered as monophyletic clade in analyses of the SSU, LSU, or LSU+SSU, the Shimodaira-Hasegawa test showed that likelihood scores for these results did not differ significantly from topologies where the Bivalvia were monophyletic. Analyses of LSU sequences strongly contradict the widely accepted Diasoma hypotheses that bivalves and scaphopods are closely related to one another. The data are consistent with recent morphological and SSU analyses suggesting scaphopods are more closely related to gastropods and cephalopods than to bivalves. The dataset also presents the first published DNA sequences from a neomeniomorph aplacophoran, a group considered critical to our understanding of the origin and early radiation of the Mollusca.     

Seed plant phylogeny: gnetophytes are derived conifers and a sister group to Pinaceae

The phylogenetic position of gnetophytes has long been controversial. We sequenced parts of the genes coding for the largest subunit of nuclear RNA polymerase I, II, and III and combined these sequences with those of four chloroplast genes, two mitochondrial genes, and 18S rRNA genes to address this issue. Both maximum likelihood and maximum parsimony analyses of the sites not affected by high substitution levels strongly support a phylogeny where gymnosperms and angiosperms are monophyletic, where cycads are at the base of gymnosperm tree and are followed by ginkgos, and where gnetophytes are grouped within conifers as the sister group of pines. The evolution of several morphological and molecular characters of gnetophytes and conifers will therefore need to be reinterpreted.     

A phylogeny of anisopterous dragonflies (Insecta,Odonata) using mtRNA genes and mixed nucleotide/doublet models

The application of mixed nucleotide/doublet substitution models has recently received attention in RNA‐based phylogenetics. Within a Bayesian approach, it was shown that mixed models outperformed analyses relying on simple nucleotide models. We analysed an mt RNA data set of dragonflies representing all major lineages of Anisoptera plus outgroups, using a mixed model in a Bayesian and parsimony (MP) approach. We used a published mt 16S rRNA secondary consensus structure model and inferred consensus models for the mt 12S rRNA and tRNA valine. Secondary structure information was used to set data partitions for paired and unpaired sites on which doublet or nucleotide models were applied, respectively. Several different doublet models are currently available of which we chose the most appropriate one by a Bayes factor test. The MP reconstructions relied on recoded data for paired sites in order to account for character covariance and an application of the ratchet strategy to find most parsimonious trees. Bayesian and parsimony reconstructions are partly differently resolved, indicating sensitivity of the reconstructions to model specification. Our analyses depict a tree in which the damselfly family Lestidae is sister group to a monophyletic clade Epiophlebia + Anisoptera, contradicting recent morphological and molecular work. In Bayesian analyses, we found a deep split between Libelluloidea and a clade ‘Aeshnoidea’ within Anisoptera largely congruent with Tillyard’s early ideas of anisopteran evolution, which had been based on evidently plesiomorphic character states. However, parsimony analysis did not support a clade ‘Aeshnoidea’, but instead, placed Gomphidae as sister taxon to Libelluloidea. Monophyly of Libelluloidea is only modestly supported, and many inter‐family relationships within Libelluloidea do not receive substantial support in Bayesian and parsimony analyses. We checked whether high Bayesian node support was inflated owing to either: (i) wrong secondary consensus structures; (ii) under‐sampling of the MCMC process, thereby missing other local maxima; or (iii) unrealistic prior assumptions on topologies or branch lengths. We found that different consensus structure models exert strong influence on the reconstruction, which demonstrates the importance of taxon‐specific realistic secondary structure models in RNA phylogenetics.     

Relative rates of evolution in the coding and control regions of African mtDNAs

Reduced median networks of African haplogroup L mitochondrial DNA (mtDNA) sequences were analyzed to determine the pattern of substitutions in both the noncoding control and coding regions. In particular, we attempted to determine the causes of the previously reported (Howell et al. 2004) violation of the molecular clock during the evolution of these sequences. In the coding region, there was a significantly higher rate of substitution at synonymous sites than at nonsynonymous sites as well as in the tRNA and rRNA genes. This is further evidence for the operation of purifying selection during human mtDNA evolution. For most sites in the control region, the relative rate of substitution was similar to the rate of neutral evolution (assumed to be most closely approximated by the substitution rate at 4-fold degenerate sites). However, there are a number of mutational hot spots in the control region, approximately 3% of the total sites, that have a rate of substitution greater than the neutral rate, at some sites by more than an order of magnitude. It is possible either that these sites are evolving under conditions of positive selection or that the substitution rate at some sites in the control region is strongly dependent upon sequence context. Finally, we obtained preliminary evidence for "nonideal" evolution in the control region, including haplogroup-specific substitution patterns and a decoupling between relative rates of substitution in the control and coding regions.     

The effects of nucleotide substitution model assumptions on estimates of nonparametric bootstrap support

The use of parameter-rich substitution models in molecular phylogenetics has been criticized on the basis that these models can cause a reduction both in accuracy and in the ability to discriminate among competing topologies. We have explored the relationship between nucleotide substitution model complexity and nonparametric bootstrap support under maximum likelihood (ML) for six data sets for which the true relationships are known with a high degree of certainty. We also performed equally weighted maximum parsimony analyses in order to assess the effects of ignoring branch length information during tree selection. We observed that maximum parsimony gave the lowest mean estimate of bootstrap support for the correct set of nodes relative to the ML models for every data set except one. For several data sets, we established that the exact distribution used to model among-site rate variation was critical for a successful phylogenetic analysis. Site-specific rate models were shown to perform very poorly relative to gamma and invariable sites models for several of the data sets most likely because of the gross underestimation of branch lengths. The invariable sites model also performed poorly for several data sets where this model had a poor fit to the data, suggesting that addition of the gamma distribution can be critical. Estimates of bootstrap support for the correct nodes often increased under gamma and invariable sites models relative to equal rates models. Our observations are contrary to the prediction that such models cause reduced confidence in phylogenetic hypotheses. Our results raise several issues regarding the process of model selection, and we briefly discuss model selection uncertainty and the role of sensitivity analyses in molecular phylogenetics.     

Deuterostome phylogeny and the sister group of the chordates: evidence from molecules and morphology

Complete coding regions of the 18S rRNA gene of an enteropneust hemichordate and an echinoid and ophiuroid echinoderm were obtained and aligned with 18S rRNA gene sequences of all major chordate clades and four outgroups. Gene sequences were analyzed to test morphological character phylogenies and to assess the strength of the signal. Maximum- parsimony analysis of the sequences fails to support a monophyletic Chordata; the urochordates form the sister taxon to the hemichordates, and together this clade plus the echinoderms forms the sister taxon to the cephalochordates plus craniates. Decay, bootstrap, and tree-length distribution analyses suggest that the signal for inference of dueterostome phylogeny is weak in this molecule. Parsimony analysis of morphological plus molecular characters supports both monophyly of echinoderms plus enteropneust hemichordates and a sister group relationship of this clade to chordates. Evolutionary parsimony does not support chordate monophyly. Neighbor-joining, Fitch-Margoliash, and maximum-likelihood analyses support a chordate lineage that is the sister group to an echinoderm-plus-hemichordate lineage. The results illustrate both the limitations of the 18S rRNA molecule alone for high- level phylogeny inference and the importance of considering both molecular and morphological data in phylogeny reconstruction.      

Evolutionary models accounting for layers of selection in protein-coding genes and their impact on the inference of positive selection

The selective forces acting on a protein-coding gene are commonly inferred using evolutionary codon models by contrasting the rate of nonsynonymous substitutions to the rate of synonymous substitutions. These models usually assume that the synonymous substitution rate, Ks, is homogenous across all sites, which is justified if synonymous sites are free from selection. However, a growing body of evidence indicates that the DNA and RNA levels of protein-coding genes are subject to varying degrees of selective constraints due to various biological functions encoded at these levels. In this paper, we develop evolutionary models that account for these layers of selection by allowing for both among-site variability of substitution rates at the DNA/RNA level (which leads to Ks variability among protein-coding sites) and among-site variability of substitution rates at the protein level (Ka variability). These models are constructed so that positive selection is either allowed or not. This enables statistical testing of positive selection when variability at the DNA/RNA substitution rate is accounted for. Using this methodology, we show that variability of the baseline DNA/RNA substitution rate is a widespread phenomenon in coding sequence data of mammalian genomes, most likely reflecting varying degrees of selection at the DNA and RNA levels. Additionally, we use simulations to examine the impact that accounting for the variability of the baseline DNA/RNA substitution rate has on the inference of positive selection. Our results show that ignoring this variability results in a high rate of erroneous positive-selection inference. Our newly developed model, which accounts for this variability, does not suffer from this problem and hence provides a likelihood framework for the inference of positive selection on a background of variability in the baseline DNA/RNA substitution rate.     

Evaluation of Morphological Variation in the Lichen Diploschistes ocellatus (Ascomycota, Ostropales): Evidence from Nuclear rDNA ITS Sequence Data

Abstract: A morphologically distinct population (var. almeriensis ) of the lichen species Diploschistes ocellatus from southern Spain is examined using morphological and molecular methods. While the morphological examination revealed some unique features of var. almeriensis, the ITS sequence data did not support a distinction from D. ocellatus s. str. The Kishino-Hasegawa test did not rule out trees distinguishing between the two morphotypes, although spectral analytical data support an inclusion of var. almeriensis in D. ocellatus. The influence of potentially erronously aligned ambiguous sites is evaluated by analysis of different alignments. However, all alignments agree in placing all specimens of D. ocellatus together and, consequently, in the consensus tree D. ocellatus s. lat. forms a monophyletic clade with high bootstrap support.     

Phylogenetic relationships of the Raphidophyceae and Xanthophyceae as inferred from nucleotide sequences of the 18S ribosomal RNA gene

Some earlier studies suggested an evolutionary relationship between the Raphidophyceae (chloromonads) and Xanthophyceae (yellow-green algae), whereas other studies suggested relationships with different algal classes or the öomycete fungi. To evaluate the relationships, we determined the complete nucleotide sequences of the 18S ribosomal RNA gene from the raphidophytes Vacuolaria virescens, Chattonella subsalsa, and Heterosigma carterae, and the xanthophytes Vaucheria bursata, Botrydium stoloniferum, Botrydiopsis intercedens, and Xanthonema debile. The results showed that the Xanthophyceae were most closely related to the Phaeophyceae. A cladistic analysis of combined data sets (nucleotide sequences, ultrastructure, and pigments) suggested the Raphidophyceae are the sister taxon to the Phaeophyceae-Xanthophyceae clade, but the bootstrap value was low (40%). The raphidophyte genera were united with high (100%) bootstrap values, supporting a hypothesis based upon ultrastructural features that marine and freshwater raphidophytes form a monophyletic group. We examined the relationship between Vaucheria, a siphoneous xanthophyte alga, and the öomycetes, and we confirmed that Vaucheria is a member of the class Xanthophyceae. Partial nucleotide sequences of the 18S rRNA gene from eight xanthophytes (including Bumillariopsis filiformis, Heterococcus caespitiosus, and Mischococcus sphaerocephalus) produce a phylogeny that is not congruent with the current morphology-based classification scheme.