Actin Coding Regions: Gene Family Evolution and Use as a Phylogenetic Marker

Selection and paucity of phylogenetic signal challenge the utility of alpha-tubulin in reconstruction of evolutionary history of free-living litostomateans (Protista, Ciliophora)

The class Litostomatea represents a highly diverse but monophyletic group, uniting both free-living and endosymbiotic ciliates. Ribosomal RNA genes and ITS-region sequences helped to recognize and define the main litostomatean lineages, but did not provide enough phylogenetic signal to unambiguously resolve their interrelationships. In this study, we attempted to improve the resolution among main free-living predatory lineages by adding the gene coding for alpha-tubulin. However, our phylogenetic analyses challenged the performance of alpha-tubulin in reconstruction of evolutionary history of free-living litostomateans. We identified several mutually interconnected problems associated with the ciliate alpha-tubulin gene: the paucity of phylogenetic signal, molecular homoplasies and non-neutral evolution. Positive selection may generate molecular homoplasies (parallel evolution), while negative selection may cause a small number of changes and hence little phylogenetic informativness. Both problems were encountered in nucleotide and amino acid alpha-tubulin alignments, indicating an action of various selective pressures. Taking into account the involvement of alpha-tubulin in many essential biological processes, this protein could be so strongly affected by purifying selection that it even might have become an inappropriate molecular marker for reconstruction of phylogenetic relationships. Therefore, a great caution should be paid when tubulin genes are included in phylogenetic and/or phylogenomic analyses.

Actin evolution in ciliates (Protist, Alveolata) is characterized by high diversity and three duplication events

Ciliates possess two distinct nuclear genomes and unique genomic features, including highly fragmented chromosomes and extensive chromosomal rearrangements. Recent transcriptomic surveys have revealed that ciliates have several multi-copy genes providing an ideal template to study gene family evolution. Nonetheless, this process remains little studied in ciliated protozoa and consequently, the evolutionary patterns that govern it are not well understood. In this study, we focused on obtaining fine-scale information relative to ciliate species divergence for the first time. A total of 230 actin gene sequences were derived from this study, among which 217 were from four closely related Pseudokeronopsis species and 13 from other hypotrichous ciliates. Our investigation shows that: (1) At least three duplication events occurred in ciliates: diversification of three actin genes (Actin I, II, III) happened after the divergence of ciliate classes but before that of subclasses. And several recent and genus-specific duplications were followed within Actin I (Sterkiella, Oxytricha, Uroleptus, etc.), Actin II (Sterkiella), respectively. (2) Within the genus Pseudokeronopsis, Actin I gene duplication events happened after P. carnea and P. erythrina diverged. In contrast, in the morphologically similar species P. flava and P. rubra, the duplication event preceded diversification of the two species. The Actin II gene duplication events preceded divergence of the genus Pseudokeronopsis. (3) Phylogenetic analyses revealed that actin is suitable for resolving ciliate classes, but may not be used to infer lower taxon relationships.

Unexpected dynamic gene family evolution in algal actins

Actin is a conserved cytoskeletal protein that is well studied in model organisms although much less is known about actin molecular evolution in taxonomically diverse algae. Here, we analyzed 107 novel partial algal actin sequences and report some unexpected results. First, monophyletic actin gene families in multiple, phylogenetically distantly related algal taxa contain two distinct clades of sequences. One of these clades contains highly conserved sequences, whereas the second has multiple members with a significantly elevated substitution rate. This rate difference is associated with an excess of synonymous substitutions, strongly suggesting that both isoforms are active. These results paint a novel picture of actin gene evolution in algae showing it to be a remarkably dynamic system with duplication, homogenization, and potential functional diversification occurring independently in distantly related lineages.

Genome architecture drives protein evolution in ciliates

Studies of microbial eukaryotes have been pivotal in the discovery of biological phenomena, including RNA editing, self-splicing RNA, and telomere addition. Here we extend this list by demonstrating that genome architecture, namely the extensive processing of somatic (macronuclear) genomes in some ciliate lineages, is associated with elevated rates of protein evolution. Using newly developed likelihood-based procedures for studying molecular evolution, we investigate 6 genes to compare 1) ciliate protein evolution to that of 3 other clades of eukaryotes (plants, animals, and fungi) and 2) protein evolution in ciliates with extensively processed macronuclear genomes to that of other ciliate lineages. In 5 of the 6 genes, ciliates are estimated to have a higher ratio of nonsynonymous/synonymous substitution rates, consistent with an increase in the rate of protein diversification in ciliates relative to other eukaryotes. Even more striking, there is a significant effect of genome architecture within ciliates as the most divergent proteins are consistently found in those lineages with the most highly processed macronuclear genomes. We propose a model whereby genome architecture-specifically chromosomal processing, amitosis within macronuclei, and epigenetics-allows ciliates to explore protein space in a novel manner. Further, we predict that examination of diverse eukaryotes will reveal additional evidence of the impact of genome architecture on molecular evolution.

Highly divergent actins from karyorelictean, heterotrich, and litostome ciliates

We have cloned, sequenced, and characterized cDNA of actins from five ciliate species of three different classes of the phylum Ciliophora: Karyorelictea (Loxodes striatus), Heterotrichea (Blepharisma japonicum, Blepharisma musculus), and Litostomatea (Didinium nasutum, Dileptus margaritifer). Loxodes striatus uses UGA as the stop codon and has numerous in-frame UAA and UAG, which are translated into glutamine. The other four species use UAA as the stop codon and have no in-frame UAG nor UGA. The putative amino acid sequences of the newly determined actin genes were found to be highly divergent as expected from previous findings of other ciliate actins. These sequences were also highly divergent from other ciliate actins, indicating that actin genes are highly diverse even within the phylum Ciliophora. Phylogenetic analysis showed high evolutionary rate of ciliate actins. Our results suggest that the evolutionary rate was accelerated because of the differences in molecular interactions.

Dramatic diversity of ciliate histone H4 genes revealed by comparisons of patterns of substitutions and paralog divergences among eukaryotes

The accumulation of divergent histone H4 amino acid sequences within and between ciliate lineages challenges traditional views of the evolution of this essential eukaryotic protein. We analyzed histone H4 sequences from 13 species of ciliates and compared these data with sequences from well-sampled eukaryotic clades. Ciliate histone H4s differ from one another at as many as 46% of their amino acids, in contrast with the highly conserved character of this protein in most other eukaryotes. Equally striking, we find paralogs of histone H4 within ciliate genomes that differ by up to 25% of their amino acids, whereas paralogs in other eukaryotes share identical or nearly identical amino acid sequences. Moreover, the most divergent H4 proteins within ciliates are found in the lineages with highly processed macronuclear genomes. Our analyses demonstrate that the dual nature of ciliate genomes-the presence of a "germline" micronucleus and a "somatic" macronucleus within each cell-allowed the dramatic variation in ciliate histone genes by altering functional constraints or enabling adaptive evolution of the histone H4 protein, or both.

A kingdom-level phylogeny of eukaryotes based on combined protein data

Current understanding of the higher order systematics of eukaryotes relies largely on analyses of the small ribosomal subunit RNA (SSU rRNA). Independent testing of these results is still limited. We have combined the sequences of four of the most broadly taxonomically sampled proteins available to create a roughly parallel data set to that of SSU rRNA. The resulting phylogenetic tree shows a number of striking differences from SSU rRNA phylogeny, including strong support for most major groups and several major supergroups.

Symbiotic origin of a novel actin gene in the cryptophyte Pyrenomonas helgolandii

Cryptophytes are photosynthetic protists that have acquired their plastids through the secondary symbiotic uptake of a red alga. A remarkable feature of cryptophytes is that they maintain a reduced form of the red algal nucleus, the nucleomorph, between the second and third plastid membranes (periplastidial compartment; PC). The nucleomorph is thought to be a transition state in the evolution of secondary plastids, with this genome ultimately being lost when photosynthesis comes under full control of the "host" nucleus (e.g., as in heterokonts, haptophytes, and euglenophytes). Genes presently found in the nucleomorph seem to be restricted to those involved in its own maintenance and to that of the plastid; other genes were lost as the endosymbiont was progressively reduced to its present state. Surprisingly, we found that the cryptophyte Pyrenomonas helgolandii possesses a novel type of actin gene that originated from the nucleomorph genome of the symbiont. Our results demonstrate for the first time that secondary symbionts can contribute genes to the host lineage which are unrelated to plastid function. These genes are akin to the products of gene duplication or lateral transfer and provide a source of evolutionary novelty that can significantly increase the genetic diversity of the host lineage. We postulate that this may be a common phenomenon in algae containing secondary plastids that has yet to be fully appreciated due to a dearth of evolutionary studies of nuclear genes in these taxa.

Actin gene family evolution and the phylogeny of coleoid cephalopods (Mollusca: Cephalopoda)

Phylogenetic analysis conducted on a 784-bp fragment of 82 actin gene sequences of 44 coleoid cephalopod taxa, along with results obtained from genomic Southern blot analysis, confirmed the presence of at least three distinct actin loci in coleoids. Actin isoforms were characteri zed through phylogenetic analysis of representative cephalopod sequences from each of the three isoforms, along with translated actin cDNA sequences from a diverse array of metazoan taxa downloaded from GenBank. One of the three isoforms found in cephalopods was closely related to actin sequences expressed in the muscular tissues of other molluscs. A second isoform was most similar to cytoplasmic-specific actin amino acid sequences. The muscle type actins of molluscs were found to be distinct from those of arthropods, suggesting at least two independent derivations of muscle actins in the protostome lineage, although statistical support for this conclusion was lacking. Parsimony and maximum-likelihood analyses of two of the isoforms from which >30 orthologous coleoid sequences had been obtained (one of the cytoplasmic actins and the muscle actin) supported the monophyly of several higher-level coleoid taxa. These included the superorders Octopodiformes and Decapodiformes, the order Octopoda, the octopod suborder Incirrata, and the teuthoid suborder Myopsida. The monophyly of several taxonomic groups within the Decapodiformes was not supported, including the orders Teuthoidea and Sepioidea and the teuthoid suborder Oegopsida. Parametric bootstrap analysis conducted on the simulated cytoplasmic actin data set provided statistical support to reject the monophyly of the Sepioidea. Although parametric bootstrap analysis of the muscle actin isoform did not reject sepioid monophyly at the 5% level, the results (rejection at P: = 0.068) were certainly suggestive of sepioid nonmonophyly.

A Search for the Origins of Animals and Fungi: Comparing and Combining Molecular Data

Green plants, animals, and fungi have long held our interest as complex, largely multicellular eukaryotes of indeterminate origin. Considerable progress has now been made toward understanding the evolutionary relationships among these taxa as well as identifying their closest protistan relatives. An exclusive animal-fungal clade (the Opisthokonta) is now widely accepted based on an insertion in the protein synthesis elongation factor 1alpha (EF-1alpha) and molecular phylogenies of ribosomal RNAs and the conservative proteins actin, alpha-tubulin, beta-tubulin, and EF-1alpha. Protein data also suggest that the cellular (dictyostelid) and acellular (myxogastrid) slime molds are a close outgroup to the animal-fungal clade. Subsequent sequencing and phylogenetic analysis of EF-1alpha sequences very strongly support a monophyletic slime mold clade (the Mycetozoa or Eumycetozoa), which also includes the lesser-known protostelid slime molds. Monophyly of the opisthokont and mycetozoan clades, exclusive of green plants, is suggested by individual analyses of EF-1alpha and actin and given strong support by concatenated protein data. Neither the monophyly of the slime molds nor their close relationship to animals and fungi are consistently supported by ribosomal RNA data. Thus, it appears unlikely that any single molecule will accurately reconstruct all higher-order taxonomy.

Actin gene sequence from Euglena gracilis

The full length coding sequence of the Euglena gracilis actin gene was determined by RT-PCR of Euglena gracilis mRNA. Conserved regions in the actin amino acid sequence were used as guides for the synthesis of degenerate primers. Sequence was obtained for 1,238 nucleotides, of which 1,131 were coding for 377 amino acids. Sequence comparisons showed a similarity with other actins of 56% to 80%. Even though most of the actin amino acid sequence was conserved, some regions showed high divergence, i.e. the DNase I-binding loop at the N-terminal region. The construction of a phylogenetic tree based on actin sequences from different organisms placed Euglena gracilis in a cluster with Trypanosoma brucei and Leishmania major.

The tigA gene is a transcriptional fusion of glycolytic genes encoding triose-phosphate isomerase and glyceraldehyde-3-phosphate dehydrogenase in oomycota

Genes encoding triose-phosphate isomerase (TPI) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) are fused and form a single transcriptional unit (tigA) in Phytophthora species, members of the order Pythiales in the phylum Oomycota. This is the first demonstration of glycolytic gene fusion in eukaryotes and the first case of a TPI-GAPDH fusion in any organism. The tigA gene from Phytophthora infestans has a typical Oomycota transcriptional start point consensus sequence and, in common with most Phytophthora genes, has no introns. Furthermore, Southern and PCR analyses suggest that the same organization exists in other closely related genera, such as Pythium, from the same order (Oomycota), as well as more distantly related genera, Saprolegnia and Achlya, in the order Saprolegniales. Evidence is provided that in P. infestans, there is at least one other discrete copy of a GAPDH-encoding gene but not of a TPI-encoding gene. Finally, a phylogenetic analysis of TPI does not place Phytophthora within the assemblage of crown eukaryotes and suggests TPI may not be particularly useful for resolving relationships among major eukaryotic groups.

Origin and evolution of the slime molds (Mycetozoa)

The Mycetozoa include the cellular (dictyostelid), acellular (myxogastrid), and protostelid slime molds. However, available molecular data are in disagreement on both the monophyly and phylogenetic position of the group. Ribosomal RNA trees show the myxogastrid and dictyostelid slime molds as unrelated early branching lineages, but actin and beta-tubulin trees place them together as a single coherent (monophyletic) group, closely related to the animal-fungal clade. We have sequenced the elongation factor-1alpha genes from one member of each division of the Mycetozoa, including Dictyostelium discoideum, for which cDNA sequences were previously available. Phylogenetic analyses of these sequences strongly support a monophyletic Mycetozoa, with the myxogastrid and dictyostelid slime molds most closely related to each other. All phylogenetic methods used also place this coherent Mycetozoan assemblage as emerging among the multicellular eukaryotes, tentatively supported as more closely related to animals + fungi than are green plants. With our data there are now three proteins that consistently support a monophyletic Mycetozoa and at least four that place these taxa within the "crown" of the eukaryote tree. We suggest that ribosomal RNA data should be more closely examined with regard to these questions, and we emphasize the importance of developing multiple sequence data sets.

Evidence for the Heterolobosea from phylogenetic analysis of genes encoding glyceraldehyde-3-phosphate dehydrogenase

The phylogenetic relationships between major slime mould groups and the identification of their unicellular relatives has been a subject of controversy for many years. Traditionally, it has been assumed that two slime mould groups, the acrasids and the dictyostelids were related by virtue of their cellular slime mould habit; a view still endorsed by at least one current classification scheme. However, a decade ago, on the basis of detailed ultrastructural resemblances it was proposed that acrasids of the family Acrasidae were not relatives of other slime moulds but instead related to a group of mostly free-living unicellular amoebae, the Schizopyrenida. The class Heterolobosea was created to contain these organisms and has since figured in many discussions of protist evolution. We sought to test the validity of Heterolobosea by characterizing homologs of the highly conserved glycolytic enzyme glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from an acrasid, Acrasis rosea; a dictyostelid, Dictyostelium discoideum; and the schizopyrenid Naegleria andersoni. Phylogenetic analysis of these and other GAPDH sequences, using maximum parsimony, neighbour-joining distance and maximum likelihood methods strongly supports the Heterolobosea hypothesis and discredits the concept of a cellular slime mould grouping. Moreover, all of our analyses place Dictyostelium discoideum as a relatively recently originating lineage, most closely related to the Metazoa, similar to other recently published phylogenies of protein-coding genes. However, GAPDH phylogenies do not show robust branching orders for most of the relationships between major groups. We propose that several of the incongruencies observed between GAPDH and other molecular phylogenies are artifacts resulting from substitutional saturation of this enzyme.