Heterotachy processes in rhodophyte-derived secondhand plastid genes: Implications for addressing the origin and evolution of dinoflagellate plastids

Serial transfer of plastids from one eukaryotic host to another is the key process involved in evolution of secondhand plastids. Such transfers drastically change the environment of the plastids and hence the selection regimes, presumably leading to changes over time in the characteristics of plastid gene evolution and to misleading phylogenetic inferences. About half of the dinoflagellate protists species are photosynthetic and unique in harboring a diversity of plastids acquired from a wide range of eukaryotic algae. They are therefore ideal for studying evolutionary processes of plastids gained through secondary and tertiary endosymbioses. In the light of these processes, we have evaluated the origin of 2 types of dinoflagellate plastids, containing the peridinin or 19'-hexanoyloxyfucoxanthin (19'-HNOF) pigments, by inferring the phylogeny using "covarion" evolutionary models allowing the pattern of among-site rate variation to change over time. Our investigations of genes from secondary and tertiary plastids derived from the rhodophyte plastid lineage clearly reveal "heterotachy" processes characterized as stationary covarion substitution patterns and changes in proportion of variable sites across sequences. Failure to accommodate covarion-like substitution patterns can have strong effects on the plastid tree topology. Importantly, multigene analyses performed with probabilistic methods using among-site rate and covarion models of evolution conflict with proposed single origin of the peridinin- and 19'-HNOF-containing plastids, suggesting that analysis of secondhand plastids can be hampered by convergence in the evolutionary signature of the plastid DNA sequences. Another type of sequence convergence was detected at protein level involving the psaA gene. Excluding the psaA sequence from a concatenated protein alignment grouped the peridinin plastid with haptophytes, congruent with all DNA trees. Altogether, taking account of complex processes involved in the evolution of dinoflagellate plastid sequences (both at the DNA and amino acid level), we demonstrate the difficulty of excluding independent, tertiary origin for both the peridinin and 19'-HNOF plastids involving engulfment of haptophyte-like algae. In addition, the refined topologies suggest the red algal order, Porphyridales, as the endosymbiont ancestor of the secondary plastids in cryptophytes, haptophytes, and heterokonts.     

Heterotachy and tree building: a case study with plastids and eubacteria

The nature of heterotachy at the center of recent controversy over the relative performance of tree-building methods is different from the form of heterotachy that has been inferred in empirical studies. The latter have suggested that proportions of variable sites (p(var)) vary among orthologues and among paralogues. However, the strength of this inference, describing what may be one of the most important evolutionary properties of sequence data, has remained weak. Consequently, other models of sequence evolution have been proposed to explain some long-branch attraction (LBA) problems that could be attributed to differences in p(var). For an empirical case with plastid and eubacterial RNA polymerase sequences, we confirm using capture-recapture estimates and simulations that p(var) can differ among orthologues in anciently diverged evolutionary lineages. We find that parsimony and a least squares distance method that implements an overly simple model of sequence evolution are susceptible to LBA induced by this form of heterotachy. Although homogeneous maximum likelihood inference was found to be robust to model misspecification in our specific example, we caution against assuming that it will always be so.     

Combined heat shock protein 90 and ribosomal RNA sequence phylogeny supports multiple replacements of dinoflagellate plastids

Dinoflagellates harbour diverse plastids obtained from several algal groups, including haptophytes, diatoms, cryptophytes, and prasinophytes. Their major plastid type with the accessory pigment peridinin is found in the vast majority of photosynthetic species. Some species of dinoflagellates have other aberrantly pigmented plastids. We sequenced the nuclear small subunit (SSU) ribosomal RNA (rRNA) gene of the "green" dinoflagellate Gymnodinium chlorophorum and show that it is sister to Lepidodinium viride, indicating that their common ancestor obtained the prasinophyte (or other green alga) plastid in one event. As the placement of dinoflagellate species that acquired green algal or haptophyte plastids is unclear from small and large subunit (LSU) rRNA trees, we tested the usefulness of the heat shock protein (Hsp) 90 gene for dinoflagellate phylogeny by sequencing it from four species with aberrant plastids (G. chlorophorum, Karlodinium micrum, Karenia brevis, and Karenia mikimotoi) plus Alexandrium tamarense, and constructing phylogenetic trees for Hsp90 and rRNAs, separately and together. Analyses of the Hsp90 and concatenated data suggest an ancestral origin of the peridinin-containing plastid, and two independent replacements of the peridinin plastid soon after the early radiation of the dinoflagellates. Thus, the Hsp90 gene seems to be a promising phylogenetic marker for dinoflagellate phylogeny.     

Modelling heterotachy in phylogenetic inference by reversible-jump Markov chain Monte Carlo

The rate at which a given site in a gene sequence alignment evolves over time may vary. This phenomenon--known as heterotachy--can bias or distort phylogenetic trees inferred from models of sequence evolution that assume rates of evolution are constant. Here, we describe a phylogenetic mixture model designed to accommodate heterotachy. The method sums the likelihood of the data at each site over more than one set of branch lengths on the same tree topology. A branch-length set that is best for one site may differ from the branch-length set that is best for some other site, thereby allowing different sites to have different rates of change throughout the tree. Because rate variation may not be present in all branches, we use a reversible-jump Markov chain Monte Carlo algorithm to identify those branches in which reliable amounts of heterotachy occur. We implement the method in combination with our 'pattern-heterogeneity' mixture model, applying it to simulated data and five published datasets. We find that complex evolutionary signals of heterotachy are routinely present over and above variation in the rate or pattern of evolution across sites, that the reversible-jump method requires far fewer parameters than conventional mixture models to describe it, and serves to identify the regions of the tree in which heterotachy is most pronounced. The reversible-jump procedure also removes the need for a posteriori tests of 'significance' such as the Akaike or Bayesian information criterion tests, or Bayes factors. Heterotachy has important consequences for the correct reconstruction of phylogenies as well as for tests of hypotheses that rely on accurate branch-length information. These include molecular clocks, analyses of tempo and mode of evolution, comparative studies and ancestral state reconstruction. The model is available from the authors' website, and can be used for the analysis of both nucleotide and morphological data.     

Functional Correlation in Amino Acid Residue Mutations of Yeast Iso-2-Cytochrome c that Is Consistent with the Prediction of the Concomitantly Variable Codon Theory in Cytochrome c Evolution

Fitch and Markowitz' theory of concomitantly variable codons (covarions) in evolution predicted the existence of functional correlation in amino acid residue mutations among present-day cytochromes c. Mutational analysis was carried out on yeast iso-2-cytochrome c, where hydrophobic core residues I20, M64, L85, and M98 and surface residue L9 were mutated, in selected combinations, to those found in mammalian and bird cytochromes c. The functionality assay is based upon the ability of yeast cells to grow in YPGE medium. Furthermore, experiments on the single M64L and M98L mutations as well as the double M64L/M98L mutation using NMR showed that the effects of these mutations are to perturb the structural integrity of the protein. We identified functional correlation in two cases of a pair of residue mutations, the I20 V and M98 L pair and the L9 I and L85 I pair. In both cases, only one of the two alternative, putative evolutionary pathways leads to a functional protein and the corresponding pairs of residue mutations are among those found in present-day cytochromes c. Since valine is predicted to be at position 20 in the ancestral form of cytochrome c, the present data provide an explanation for the ancient requirement of leucine rather than methionine in position 98. The present data provide further evidence for the role of those specific atom–atom interactions in directing a pathway in the evolutionary changes of the amino acid sequence that have taken place in cytochrome c, in accordance with Fitch and Markowitz.     

Covarion structure in plastid genome evolution: a new statistical test

Covarion models of molecular evolution allow the rate of evolution of a site to vary through time. There are few simple and effective tests for covarion evolution, and consequently, little is known about the presence of covarion processes in molecular evolution. We describe two new tests for covarion evolution and demonstrate with simulations that they perform well under a wide range of conditions. A survey of covarion evolution in sequenced plastid genomes found evidence of covarion drift in at least 26 out of 57 genes. Covarion evolution is most evident in first and second codon positions of the plastid genes, and there is no evidence of covarion evolution in third codon positions. Therefore, the significant covarion tests are likely due to changes in the selective constraints of amino acids. The frequency of covarion evolution within the plastid genome suggests that covarion processes of evolution were important in generating the observed patterns of sequence variation among plastid genomes.     

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.     

An improved statistical method for detecting heterotachy in nucleotide sequences

The principle of heterotachy states that the substitution rate of sites in a gene can change through time. In this article, we propose a powerful statistical test to detect sites that evolve according to the process of heterotachy. We apply this test to an alignment of 1289 eukaryotic rRNA molecules to 1) determine how widespread the phenomenon of heterotachy is in ribosomal RNA, 2) to test whether these heterotachous sites are nonrandomly distributed, that is, linked to secondary structure features of ribosomal RNA, and 3) to determine the impact of heterotachous sites on the bootstrap support of monophyletic groupings. Our study revealed that with 21 monophyletic taxa, approximately two-thirds of the sites in the considered set of sequences is heterotachous. Although the detected heterotachous sites do not appear bound to specific structural features of the small subunit rRNA, their presence is shown to have a large beneficial influence on the bootstrap support of monophyletic groups. Using extensive testing, we show that this may not be due to heterotachy itself but merely due to the increased substitution rate at the detected heterotachous sites.