Monophyly of Primary Photosynthetic Eukaryotes: Green Plants, Red Algae, and Glaucophytes

Phylogenomic analysis identifies red algal genes of endosymbiotic origin in the chromalveolates

Endosymbiosis has spread photosynthesis to many branches of the eukaryotic tree; however, the history of photosynthetic organelle (plastid) gain and loss remains controversial. Fortuitously, endosymbiosis may leave a genomic footprint through the transfer of endosymbiont genes to the "host" nucleus (endosymbiotic gene transfer, EGT). EGT can be detected through comparison of host genomes to uncover the history of past plastid acquisitions. Here we focus on a lineage of chlorophyll c-containing algae and protists ("chromalveolates") that are postulated to share a common red algal secondary endosymbiont. This plastid is originally of cyanobacterial origin through primary endosymbiosis and is closely related among the Plantae (i.e., red, green, and glaucophyte algae). To test these ideas, an automated phylogenomics pipeline was used with a novel unigene data set of 5,081 expressed sequence tags (ESTs) from the haptophyte alga Emiliania huxleyi and genome or EST data from other chromalveolates, red algae, plants, animals, fungi, and bacteria. We focused on nuclear-encoded proteins that are targeted to the plastid to express their function because this group of genes is expected to have phylogenies that are relatively easy to interpret. A total of 708 genes were identified in E. huxleyi that had a significant Blast hit to at least one other taxon in our data set. Forty-six of the alignments that were derived from the 708 genes contained at least one other chromalveolate (i.e., besides E. huxleyi), red and/or green algae (or land plants), and one or more cyanobacteria, whereas 15 alignments contained E. huxleyi, one or more other chromalveolates, and only cyanobacteria. Detailed phylogenetic analyses of these data sets turned up 19 cases of EGT that did not contain significant paralogy and had strong bootstrap support at the internal nodes, allowing us to confidently identify the source of the plastid-targeted gene in E. huxleyi. A total of 17 genes originated from the red algal lineage, whereas 2 genes were of green algal origin. Our data demonstrate the existence of multiple red algal genes that are shared among different chromalveolates, suggesting that at least a subset of this group may share a common origin.

The largest subunit of RNA polymerase II from the Glaucocystophyta: functional constraint and short-branch exclusion in deep eukaryotic phylogeny

BACKGROUND

Evolutionary analyses of the largest subunit of RNA polymerase II (RPB1) have yielded important and at times provocative results. One particularly troublesome outcome is the consistent inference of independent origins of red algae and green plants, at odds with the more widely accepted view of a monophyletic Plantae comprising all eukaryotes with primary plastids. If the hypothesis of a broader kingdom Plantae is correct, then RPB1 trees likely reflect a persistent phylogenetic artifact. To gain a better understanding of RNAP II evolution, and the presumed artifact relating to green plants and red algae, we isolated and analyzed RPB1 from representatives of Glaucocystophyta, the third eukaryotic group with primary plastids.

RESULTS

Phylogenetic analyses incorporating glaucocystophytes do not recover a monophyletic Plantae; rather they result in additional conflicts with the most widely held views on eukaryotic relationships. In particular, glaucocystophytes are recovered as sister to several amoebozoans with strong support. A detailed investigation shows that this clade can be explained by what we call "short-branch exclusion," a phylogenetic artifact integrally associated with "long-branch attraction." Other systematic discrepancies observed in RPB1 trees can be explained as phylogenetic artifacts; however, these apparent artifacts also appear in regions of the tree that support widely held views of eukaryotic evolution. In fact, most of the RPB1 tree is consistent with artifacts of rate variation among sequences and co-variation due to functional constraints related to C-terminal domain based RNAP II transcription.

CONCLUSION

Our results reveal how subtle and easily overlooked biases can dominate the overall results of molecular phylogenetic analyses of ancient eukaryotic relationships. Sources of potential phylogenetic artifact should be investigated routinely, not just when obvious "long-branch attraction" is encountered.

The peculiar distribution of class I and class II aldolases in diatoms and in red algae

Diatom plastids probably evolved by secondary endocytobiosis from a red alga that was up by a eukaryotic host cell. Apparently, this process increased the complexity of the intracellular distribution of metabolic enzymes. We identified genes encoding fructose-bisphosphate aldolases (FBA) in two centric (Odontella sinensis, Thalassiosira pseudonana) and one pennate (Phaeodactylum tricornutum) diatoms and found that four different aldolases are present in both groups: two plastid targeted class II enzymes (FBAC1 and FBAC2), one cytosolic class II (FBA3) and one cytosolic class I (FBA4) enzyme. The pennate Phaeodactylum possesses an additional plastidic class I enzyme (FBAC5). We verified the classification of the different aldolases in the diatoms by enzymatic characterization of isolated plastids and whole cell extracts. Interestingly, our results imply that in plastids of centric and pennate diatoms mainly either class I or class II aldolases are active. We also identified genes for both class I and class II aldolases in red algal EST databases, thus presenting a fascinating example of the reutilization and recompartmentalization of different aldolase isoenzymes during secondary endocytobiosis but as well demonstrating the limited use of metabolic enzymes as markers for the interpretation of phylogenetic histories in algae.

Homologous protein import machineries in chloroplasts and cyanelles

The cyanelles of the glaucocystophyte alga Cyanophora paradoxa resemble endosymbiotic cyanobacteria, especially in the presence of a peptidoglycan wall between the inner and outer envelope membranes. However, it is now clear that cyanelles are in fact primitive plastids. Phylogenetic analyses of plastid, nuclear and mitochondrial genes support a single primary endosymbiotic event. In this scenario, cyanelles and all other plastid types are derived from an ancestral photosynthetic organelle combining the high gene content of rhodoplasts and the peptidoglycan wall of cyanelles. This means that the import apparatuses of all primary plastids, i.e. those from glaucocystophytes, red algae, green algae and higher plants, should be homologous. If this is the case, then transit sequences should be similar and heterologous import experiments feasible. Thus far, heterologous in vitro import has been shown in one direction only: precursors from C. paradoxa were imported into isolated pea or spinach chloroplasts. Cyanelle transit sequences differ from chloroplast stroma targeting peptides in containing in their N-terminal domain an invariant phenylalanine residue which is shown here to be crucial for import. In addition, we now demonstrate that heterologous precursors are readily imported into isolated cyanelles, provided that the essential phenylalanine residue is engineered into the N-terminal part of chloroplast transit peptides. The cyanelle and likely also the rhodoplast import apparatus can be envisaged as prototypes with a single receptor/channel showing this requirement for N-terminal phenylalanine. In chloroplasts, multiple receptors with overlapping and less stringent specificities have evolved, explaining the efficient heterologous import of native precursors from C. paradoxa.

Dense taxonomic EST sampling and its applications for molecular systematics of the Coleoptera (beetles)

Expressed sequence tag (EST) sequences can provide a wealth of data for phylogenetic and genomic studies, but the utility of these resources is restricted by poor taxonomic sampling. Here, we use small EST libraries (<1,000 clones) to generate phylogenetic markers across a broad sample of insects, focusing on the species-rich Coleoptera (beetles). We sequenced over 23,000 ESTs from 34 taxa, which produced 8,728 unique sequences after clustering nonredundant sequences. Between taxa, the sequences could be grouped into 731 gene clusters, with the largest corresponding to mitochondrial DNA transcripts and gene families chymotrypsin, actin, troponin, and tubulin. While levels of paralogy were high in most gene clusters, several midsized clusters including many ribosomal protein (RP) genes appeared to be free of expressed paralogs. To evaluate the utility of EST data for molecular systematics, we curated available transcripts for 66 RP genes from representatives of the major groups of Coleoptera. Using supertree and supermatrix approaches for phylogenetic analysis, the results were consistent with the emerging phylogenetic conclusions about basal relationships in Coleoptera. Numerous small EST libraries from a taxonomically densely sampled lineage can provide a core set of genes that together act as a scaffold in phylogenetic reconstruction, comparative genomics, and studies of gene evolution.

A plastid in the making: evidence for a second primary endosymbiosis

One of the major steps in the evolution of life was the origin of photosynthesis in nucleated cells underpinning the evolution of plants. It is well accepted that this evolutionary process was initiated when a photosynthetic bacterium (a cyanobacterium) was taken up by a colorless host cell, probably more than a billion years ago, and transformed into a photosynthetic organelle (a plastid) during a process known as primary endosymbiosis. Here, we use sequence comparisons and phylogenetic analyses of the prokaryotic rDNA operon to show that the thecate, filose amoeba Paulinella chromatophora Lauterborn obtained its photosynthetic organelles by a similar but more recent process, which involved a different cyanobacterium, indicating that the evolution of photosynthetic organelles from cyanobacteria was not a unique event, as is commonly believed, but may be an ongoing process.

The tree of eukaryotes

Recent advances in resolving the tree of eukaryotes are converging on a model composed of a few large hypothetical 'supergroups', each comprising a diversity of primarily microbial eukaryotes (protists, or protozoa and algae). The process of resolving the tree involves the synthesis of many kinds of data, including single-gene trees, multigene analyses, and other kinds of molecular and structural characters. Here, we review the recent progress in assembling the tree of eukaryotes, describing the major evidence for each supergroup, and where gaps in our knowledge remain. We also consider other factors emerging from phylogenetic analyses and comparative genomics, in particular lateral gene transfer, and whether such factors confound our understanding of the eukaryotic tree.

Heterotachy and long-branch attraction in phylogenetics

BACKGROUND

Probabilistic methods have progressively supplanted the Maximum Parsimony (MP) method for inferring phylogenetic trees. One of the major reasons for this shift was that MP is much more sensitive to the Long Branch Attraction (LBA) artefact than is Maximum Likelihood (ML). However, recent work by Kolaczkowski and Thornton suggested, on the basis of simulations, that MP is less sensitive than ML to tree reconstruction artefacts generated by heterotachy, a phenomenon that corresponds to shifts in site-specific evolutionary rates over time. These results led these authors to recommend that the results of ML and MP analyses should be both reported and interpreted with the same caution. This specific conclusion revived the debate on the choice of the most accurate phylogenetic method for analysing real data in which various types of heterogeneities occur. However, variation of evolutionary rates across species was not explicitly incorporated in the original study of Kolaczkowski and Thornton, and in most of the subsequent heterotachous simulations published to date, where all terminal branch lengths were kept equal, an assumption that is biologically unrealistic.

RESULTS

In this report, we performed more realistic simulations to evaluate the relative performance of MP and ML methods when two kinds of heterogeneities are considered: (i) within-site rate variation (heterotachy), and (ii) rate variation across lineages. Using a similar protocol as Kolaczkowski and Thornton to generate heterotachous datasets, we found that heterotachy, which constitutes a serious violation of existing models, decreases the accuracy of ML whatever the level of rate variation across lineages. In contrast, the accuracy of MP can either increase or decrease when the level of heterotachy increases, depending on the relative branch lengths. This result demonstrates that MP is not insensitive to heterotachy, contrary to the report of Kolaczkowski and Thornton. Finally, in the case of LBA (i.e. when two non-sister lineages evolved faster than the others), ML outperforms MP over a wide range of conditions, except for unrealistic levels of heterotachy.

CONCLUSION

For realistic combinations of both heterotachy and variation of evolutionary rates across lineages, ML is always more accurate than MP. Therefore, ML should be preferred over MP for analysing real data, all the more so since parametric methods also allow one to handle other types of biological heterogeneities much better, such as among sites rate variation. The confounding effects of heterotachy on tree reconstruction methods do exist, but can be eschewed by the development of mixture models in a probabilistic framework, as proposed by Kolaczkowski and Thornton themselves.

An Empirical Assessment of Long-Branch Attraction Artefacts in Deep Eukaryotic Phylogenomics

In the context of exponential growing molecular databases, it becomes increasingly easy to assemble large multigene data sets for phylogenomic studies. The expected increase of resolution due to the reduction of the sampling (stochastic) error is becoming a reality. However, the impact of systematic biases will also become more apparent or even dominant. We have chosen to study the case of the long-branch attraction artefact (LBA) using real instead of simulated sequences. Two fast-evolving eukaryotic lineages, whose evolutionary positions are well established, microsporidia and the nucleomorph of cryptophytes, were chosen as model species. A large data set was assembled (44 species, 133 genes, and 24,294 amino acid positions) and the resulting rooted eukaryotic phylogeny (using a distant archaeal outgroup) is positively misled by an LBA artefact despite the use of a maximum likelihood-based tree reconstruction method with a complex model of sequence evolution. When the fastest evolving proteins from the fast lineages are progressively removed (up to 90%), the bootstrap support for the apparently artefactual basal placement decreases to virtually 0%, and conversely only the expected placement, among all the possible locations of the fast-evolving species, receives increasing support that eventually converges to 100%. The percentage of removal of the fastest evolving proteins constitutes a reliable estimate of the sensitivity of phylogenetic inference to LBA. This protocol confirms that both a rich species sampling (especially the presence of a species that is closely related to the fast-evolving lineage) and a probabilistic method with a complex model are important to overcome the LBA artefact. Finally, we observed that phylogenetic inference methods perform strikingly better with simulated as opposed to real data, and suggest that testing the reliability of phylogenetic inference methods with simulated data leads to overconfidence in their performance. Although phylogenomic studies can be affected by systematic biases, the possibility of discarding a large amount of data containing most of the nonphylogenetic signal allows recovering a phylogeny that is less affected by systematic biases, while maintaining a high statistical support.