Photosynthesis moved into eukaryotic organisms through the endosymbiotic uptake of cyanobacteria that, over time, were integrated as organelles[@b1]. Primary plastids, direct descendants of this process, are found in three evolutionary lineages: red algae, green plants and glaucophytes[@b2]. Photosynthesis later spread horizontally through 'secondary' endosymbioses, resulting in most other eukaryotic algae[@b3][@b4]. These organisms are genetically complex and have been likened to 'Matryoshka' (Russian nesting dolls)[@b5], with a cyanobacterium nested within a red alga that is further nested within a secondary host cell. As most genes that are retained from endosymbionts end up in the secondary host's nuclear genome[@b6], they present difficulties in molecular phylogenetic investigations. Although the secondary movement of green plastids appears well established[@b7], the evolution of red-lineage plastids remains among the most controversial problems in broad-scale eukaryotic phylogenomics[@b8][@b9][@b10].

Four major algal taxa contain plastids descended secondarily from red algal primary plastids, the Ochrophyta, Haptophyta, Cryptophyta and Dinophyta (dinoflagellates)[@b4][@b9]. Possible evolutionary or systematic connections among these groups have been suggested for over 60 years[@b11]. In 1981, Cavalier-Smith[@b12] proposed the former three as members of a distinct kingdom, the Chromista, connecting their plastids through a single secondary endosymbiotic uptake of a red alga. This interpretation later was expanded as the chromalveolate hypothesis to include the Dinophyta, moving the proposed single endosymbiosis back to a common ancestor of chromists and the Alveolata, a group of protists that includes dinoflagellates, ciliates and apicomplexan parasites[@b13].

Both alveolates and chromists include numerous heterotrophic forms that lack plastids altogether[@b14]. Ochrophytes, in particular, have long been recognized as related most closely to various heterotrophic organisms (for example, oomycetes) based on shared cytological features[@b15][@b16]; collectively they are grouped into the Heterokonta[@b12] (or Stramenopiles[@b16]). A monophyletic Heterokonta also has been supported consistently by molecular phylogenetic analyses[@b17][@b18][@b19], and cryptophytes and haptophytes fall outside this group based on all available evidence (both cytological and molecular). Therefore, if chromalveolate plastids are related by direct descent from a single endosymbiosis, they must have been lost from aplastidial heterokont and alveolate groups. More recently, molecular studies also have identified heterotrophic relatives of both cryptophytes and haptophytes[@b20], although these relationships are more tentative because they are not, as yet, supported by the kinds of shared-derived cytological features that unite the Heterokonta.

Combined with evidence that chromist and dinoflagellate plastids all can be traced to a single red algal secondary endosymbiosis[@b5][@b8][@b21], the increasing number of plastid losses implied by the chromalveolate model have prompted various proposals of multiple or serial endosymbioses[@b8][@b9][@b22][@b23][@b24]. To date, compelling evidence in support of a single, unambiguous evolutionary pathway has not emerged[@b8], but it appears likely that the history of these plastids has included some level of reticulate rather than vertical inheritance. Additional endosymbioses would mean additional genetic layers in the genomes of these 'Matryoshka algae', further complicating their analyses through typical phylogenomic methods[@b25]. The problem is exacerbated by the generally weak tree-building signal in genomes for resolving deeper nodes of the eukaryotic tree of life[@b17], as well as biases in sequence data that result in widespread phylogenetic artifacts[@b26][@b27][@b28][@b29]. Consequently, although phylogenomic investigations easily recover most classically defined taxa, and often cluster them into larger assemblies, they have difficulty producing consistent, strongly supported relationships among major eukaryotic lineages.

The relatively weak phylogenetic associations across major eukaryotic groups can allow other statistical patterns to emerge, which we have explored to assess impacts from past endosymbioses on chromist genomes. In an earlier study of heterokonts, we uncovered a strong relationship between the overall similarity of a given query genome to other eukaryotic groups, as measured by top matches in BLAST searches, and the number of sequences in the databases from each target group[@b30]. This relationship likely reflects the combined effects of expanded paralogous gene families and a reduced chance that ancestral genes and gene families have been lost from larger and/or multiple genomes (see ref. [@b30] for expanded discussion). Basically, the more potential matches present in a database, the greater the probability that a given query sequence will produce a match simply by chance. Later, we expanded those analyses to show that the relationship holds consistently across a diversity of eukaryotes for which complete genome data were available ([Supplementary Fig. 1](#S1){ref-type="supplementary-material"}).

The stochastic properties of this relationship allow the use of regression analyses to generate null models for the expected similarity of any individual genome to those of other eukaryotic groups. Outliers, that is, taxa significantly more similar to a query than expected, can be identified from studentized residuals[@b31]. Such strong deviations rise above other sources of variation, for example, phylogenetic noise or individual horizontal gene transfers[@b32], which are captured in the typical residual variance around the regression line. Therefore, they likely reflect either a close phylogenetic relationship or the impact of endosymbiotic gene transfer (EGT)[@b6][@b33] between two more divergent taxa. If the deviation reflects a relatively recent divergence between nuclear genomes, then it is representative of consistent phylogenetic signal across the genome and the relationship is expected to be recovered using typical phylogenomic analyses[@b30] (and see [Supplementary Fig. 1](#S1){ref-type="supplementary-material"}). In the case of the three chromist algal groups it is more likely to reflect EGT, given that they have not been linked closely through phylogenetic analyses of nuclear genomes, but they clearly are connected by the common presence of red secondary plastids. Here, we present results that reject direct monophyletic inheritance of chromist plastids with plastid losses accounting for non-photosynthetic chromists and, further, provide strong evidence for a specific pathway of serial endosymbioses that accounts for the disjunct phylogenetic distribution of chromist plastids.

Results
=======

Linear regressions and detection of outliers
--------------------------------------------

We performed linear regressions and calculated studentized residuals for the number of most similar BLAST hits to all taxa using complete genomes from three of the four major chromalveolate algal lineages. Because no complete, well-annotated genome is available from dinoflagellates, we could not include them in our analyses. All regressions showed strong and significant relationships and revealed an interesting pattern with respect to outlier taxa. Although there appears to be no strong evidence to suggest a close relationship between host cells of haptophytes and heterokonts[@b19], studentized residuals show that *Emiliania huxleyi* (Haptophyta) is significantly more similar to heterokonts (\>3 s.d.) than predicted by the regression model ([Fig. 1](#f1){ref-type="fig"}). An analysis using *Guillardia theta* (Cryptophyta) as the query genome produced similar results; in this case, the association with heterokonts was not significant, although at a level (\>2 s.d.) that generally is considered to be potentially significant in analyses of studentized residuals, and worthy of further scrutiny (see ref. [@b31] for full discussion). In contrast, there was no indication that *Emiliania* and *Guillardia* genomes are more similar to each other than predicted by the regression model, regardless of which genome was used as the query ([Fig. 1](#f1){ref-type="fig"}). This latter result is interesting because it is at odds with the widely discussed Hacrobia hypothesis, which argues that cryptophytes and haptophytes are the most closely related of the three chromist algal groups[@b34]. When a randomly chosen ochrophyte genome (the diatom *Phaeodactylum tricornutum*) was used as a reciprocal query, cryptophytes and haptophytes were the two largest outliers, although neither deviated significantly from their predicted values ([Supplementary Fig. 2](#S1){ref-type="supplementary-material"} and [Supplementary Table 1](#S1){ref-type="supplementary-material"}). Thus, regressions provide strong support for an endosymbiotic association between heterokonts and haptophytes, suggest one between heterokonts and cryptophytes and give no indication that cryptophytes or haptophytes are related either phylogenetically or through a direct endosymbiosis.

The overall results of regression analyses suggest a pattern of serial endosymbioses, including the tertiary adoption of a cryptophyte by the ancestor of ochrophytes, and then a quaternary uptake of an ochrophyte into the ancestor of haptophytes. [Figure 2a](#f2){ref-type="fig"} depicts these serial transfers mapped onto a tree of eukaryotic relationships adapted from the recent phylogenetic investigation of cryptophyte and haptophyte evolution by Burki *et al*.[@b19] Because the regressions alone indicate only greater similarities than expected between query genomes to target taxa, they do not polarize plastid transfers in this direction; however, cryptophytes retain a residual red algal nucleus from the original secondary endosymbiosis[@b35], which is absent from either of the other two groups. Thus, a series of endosymbioses beginning with haptophytes is unreasonable because it requires a highly unlikely scenario of retention and loss of the nucleomorph, or the even less likely assumption that the nucleomorph was regenerated in cryptophytes. A more recent quaternary endosymbiosis between haptophytes and ochrophytes also is consistent with strong similarities between the two in plastid pigmentation and ultrastructure, which often led to the inclusion of haptophytes within ochrophyte classes in historical classification systems (see ref. [@b36] for review).

Contingency tests of model predictions
--------------------------------------

The model of plastid evolution that we propose makes explicit predictions as to where specific footprints of endosymbiosis (that is, significant evidence of EGT) should be found. First, if cryptophytes donated their plastid to establish ochrophytes, the transfer presumably occurred after the ochrophyte ancestor diverged from other heterokonts that remained heterotrophic. Thus, there should be evidence of EGT from cryptophytes into ochrophytes, but not into the heterotrophic relatives of ochrophytes. To test this prediction, we carried out Fisher's exact tests on the numbers of most similar BLAST matches to four heterokonts, two photosynthetic and two heterotrophic, to determine whether evidence of EGT from cryptophytes was significantly greater than in control genomes with no putative genealogic or endosymbiotic relationship to either group. To avoid complications from shared photosynthetic pathways, which clearly make algal genomes more similar to each other under any circumstances, only genes presumably unrelated to plastid functions were included in these analyses. Both photosynthetic heterokonts tested, the diatom *Phaeodactylum* and pelagophyte *Aureococcus*, contained highly significant evidence of EGT compared with control genomes ([Table 1](#t1){ref-type="table"}). In contrast, no significant evidence was found for EGT from cryptophytes in the heterotrophic oomycete *Phytophthora* or protist *Blastocystis*. In fact, in three of four tests the non-photosynthetic heterokont showed a weaker putative impact from cryptophyte genes than did the paired control genome ([Table 1](#t1){ref-type="table"}). Thus, Fisher's exact tests offer consistent and highly significant support for the hypothesized tertiary endosymbiosis between cryptophytes and ochrophytes.

We next tested predictions from the second hypothesized endosymbiosis, the quaternary adoption of an ochrophyte by the ancestor of haptophytes. In this case, because heterotrophic and photosynthetic heterokonts are closely related, most genes unrelated to plastid function should be relatively similar between the two. Consequently, evidence of EGT into haptophytes should be detectable from all heterokonts regardless of metabolic lifestyle. This prediction was strongly supported in all Fisher's exact tests ([Table 1](#t1){ref-type="table"}). We note that if a plastid transfer occurred in the opposite direction, from a haptophyte to the ancestor of ochrophytes, then no significant association would be expected between non-photosynthetic heterokonts and haptophytes (see cryptophyte example above). Therefore, not only do results from contingency tests support specific predictions of the model inferred objectively from regression analyses, they also provide significant evidence against alternative hypotheses with mutually exclusive predictions (see further discussion below).

Phylogenetic analyses of plastid genes
--------------------------------------

We performed one additional test of the proposed model of serial endosymbioses. Because genes retained within plastids have not moved between genomes, and were not included in our regressions and contingency tests, they represent independent data that should track the pathway of plastids through the respective chromalveolate lineages. Phylogenetic analyses using 5,818 inferred amino-acid positions from plastid genes were fully consistent with our proposed model ([Fig. 2b](#f2){ref-type="fig"}). Cryptophyte plastids emerged from within red algae and as the sister group to an ochrophyte--haptophyte clade. Although this tree is not compelling evidence by itself, particularly given the small number of chromist plastid genomes included, it takes on greater significance as an *a posteriori* test of an explicitly defined phylogenetic hypothesis.

To avoid introducing additional variables, plastid genes were sampled only from those photosynthetic organisms included in our primary investigation of nuclear genomes. Nevertheless, our results agree with inferences from broader phylogenetic analyses carried out by Baurain *et al*.[@b18] in their previous investigation of the chromalveolate model. Although neither study showed haptophyte plastids originating from within ochrophytes, as might be expected from the endosymbiotic transfer proposed, this could require much broader sampling of the global diversity of ochrophyte plastid genomes. Ochrophytes and haptophytes are believed to have emerged at around the same time during the late Proterozoic[@b37], meaning the proposed plastid transfer must have originated from an early branch of ochrophytes, certainly before the broad radiation of most extant forms. Presumably, a substantial fraction of early diverging ochrophytes have been lost to the multiple mass extinctions that have punctuated more than 500 million years of intervening evolution[@b38]. Thus, it would be surprising if initial sampling had recovered haptophyte plastids nested within ochrophytes, rather than as their sister group. These factors also could explain why the horizontally transferred bacterial *rpl36* gene, which has been cited as a diagnostic feature linking cryptophyte and haptophyte plastids[@b39], has not been found in ochrophyte plastid genomes to date.

Discussion
==========

For most of the last 15 years, the chromalveolate hypothesis has been the dominant model of red-lineage plastid evolution; it assumes a single secondary origin of plastids in the common ancestor of chromist algae, dinoflagellates and related heterotrophic taxa. Consequently, plastids must have been lost from the ancestors of chromalveolate groups where they are not currently present, including oomycetes and other aplastidial heterokonts. If the hypothesis of a continuous vertical descent of chromist plastids is correct, then what we have identified as a footprint of cryptophyte-EGT in ochrophyte genomes would reflect a phylogenetic rather than an endosymbiotic relationship. If so, that phylogenetic signal (as measured in genes unrelated to plastid function) also should be present in *Phytophthora* and *Blastocystis*, which share a more recent common ancestor with ochrophytes than do other chromist algae. Instead, our results ([Table 1](#t1){ref-type="table"}) provide significant evidence that the greater than expected similarity between cryptophyte and ochrophyte genomes does not predate the divergence between ochrophytes and aplastidial heterokonts. This is most reasonably explained by EGT from a tertiary cryptophyte endosymbiont, and is strong evidence against the vertical descent of chromist plastids from a single ancestral secondary endosymbiosis. As such, it is an explicit falsification of the core assumption of the chromalveolate hypothesis, and is in line with several other recent investigations that found evidence inconsistent with various aspects of that model[@b8][@b18][@b19][@b30]. It also emphasizes the importance of understanding how EGT has shaped photosynthetic eukaryotes, particularly 'Matryoshka algae' with complex, chimeric genomes.

Most previous inferences of EGT as evidence of a past endosymbiosis have been based on some indeterminate number of discrepant gene phylogenies suggesting the presence of foreign 'algal' genes in a given genome. Such trees can be interpreted as endosymbiotic footprints, but also could be examples of more general horizontal gene transfer or phylogenetic/sampling artifacts[@b32][@b40]. Our use of objective regression models and explicit negative controls in contingency tests shows that the footprints of EGT proposed in this study rise significantly above what can be expected from more random sources of variation within eukaryotic genomes. Nevertheless, even such highly significant support for greater than expected similarities among chromist algal genomes does not provide rigorous evidence for EGT on a gene-by-gene basis. This will require further investigations of our proposed model of serial endosymbioses, including careful and rigorous analyses of individual genes and gene sets that underlie the genome-wide signals that we uncovered.

Beyond supporting a specific pathway of endosymbiosis through chromist algae, our results also provide strong evidence against alternative models of higher-order red-lineage plastid evolution. For example, if all three, or even any two chromist algal lineages had been founded by independent secondary endosymbioses involving a red alga, we should not find significant endosymbiotic connections both between cryptophytes and ochrophytes, and between ochrophytes and haptophytes. Independent secondary endosymbioses predict that at least one of the chromist algae would show a significant endosymbiotic association with red algae only, and not with either of the other two chromist groups. We also note that, unlike independent secondary events, connecting chromist plastids through serial endosymbioses is consistent with evidence that all chromalveolate plastids trace back to a single secondary red algal endosymbiont[@b8].

Likewise, an independent tertiary endosymbiosis between cryptophytes and haptophytes is inconsistent with the absence of an association between their genomes in regression analyses ([Fig. 1](#f1){ref-type="fig"}). More importantly, it is contradicted by the highly significant support that we find for EGT from ochrophytes into haptophytes ([Fig. 1](#f1){ref-type="fig"}, [Table 1](#t1){ref-type="table"}); this would not be present had haptophytes adopted a cryptophyte endosymbiont directly. The highly significant footprint of genes from non-photosynthetic heterokonts in the haptophyte genome also rules out the endosymbiotic movement of a haptophyte into ochrophytes, as noted above.

The specific pattern of serial endosymbioses that we propose also explains conflicting phylogenetic signals uncovered previously across the complete *Emiliania huxleyi* genome[@b41]. Although the dominant tree-building signal present recovered haptophytes as sister to the large SAR[@b42] complex (stramenopiles (haptophytes), alveolates, rhizarians), \~30% of individual nuclear genes showed a specific phylogenetic affinity with heterokonts[@b41]. In contrast, no strong, concerted affinities were found for either cryptophyte or red algal genes. This is fully consistent with the quaternary endosymbiosis that we propose between ochrophytes and haptophytes, but is not easily reconciled with most alternative models of chromist plastid evolution. In conclusion, the results of our investigation represent a qualitative advance in understanding how the major chromist lineages acquired photosynthesis. In addition, they highlight the utility of new and alternative approaches for untangling the complex and reticulate nature of plastid evolution.

Methods
=======

Database creation and sequence similarity searches
--------------------------------------------------

The general approach to database construction, BLAST searches and parsing of BLAST results followed those described previously[@b30]. Briefly, to investigate signal from EGT in genomes of chromists and related taxa, a sequence database was constructed including all inferred protein-encoding genes from the complete genomes of 50 organisms distributed across 14 well-established eukaryotic lineages ([Supplementary Table 2](#S1){ref-type="supplementary-material"}). BLASTP searches were performed against this database using all inferred protein sequences as indicated from each specific query genome described in our results. Each sequence from a given query genome was used in only one search against the complete database. To identify homologues from the databases most reliably, a stringent *E*-value cut-off of 1e−20 was used in all searches ([Supplementary Data 1](#S1){ref-type="supplementary-material"}). All sequences returned using a given query genome were retained in a data set in ranked order for use in computational pipelines described below. Only complete, well-annotated genomes were used for two primary reasons. First, our statistical treatments comparing potential signals from EGT across genomes are based on the total number of genes present, which cannot be determined accurately from transcriptome data alone. In addition, because investigations of EGT seek to measure the contribution of foreign genes present in the organisms assessed, it is essential that such genes can be placed in a genomic context, and are not simply present as mRNAs that could represent external contamination.

Sorting pipelines
-----------------

We wrote a custom JAVA program ([Supplementary Data 2](#S1){ref-type="supplementary-material"}) to parse output files from BLAST searches. Relevant data, including all sequences producing significant alignments (≤e−20), their source genomes, group designation and *E*-values, were extracted and formatted along with the query sequences. We then created a database using MySQL to store these entries to facilitate multiple further computations, including various complicated SQL queries generated on the basis of the research questions arising from regression analyses. Information from BLAST output files for all sequences returned from each given query genome was stored in a specific, individual table.

The JAVA program that we developed used two inputs. The first was a category file, indicating to which of the 14 eukaryotic groups each of the 50 genomes belonged. The 14 groups were labelled with group ID numbers. The second input was the BLAST output file described above. The program generated output files in the following format:

Input sequence/similar sequence found/hit\#/similar sequence/group ID/group name/genome/score/e-value, where 'similar sequence' indicates whether any sequence above our cut-off value was returned for the query sequence. 'Hit\#' (hit-number) stands for the sequential rank number of a given sequence returned in a BLAST search. For the results reported here, we focused largely on best matches, after first screening out matches to members of the target group to which the query genome belonged, including all self-hits to the query genome itself. Output tables constructed in this format could be loaded directly to our database and sorted in whatever manner required for a given statistical test. We then constructed various SQL queries to conduct counting and parsing functions against the data file. As a screen for genes unrelated to photosynthesis, we included a command that the gene must be present in at least two genomes from the broadly sampled heterotrophic taxa metazoa and/or fungi.

Statistical analyses
--------------------

Linear regressions were carried out in SPSS (version 20) without a constant, given that a data set with zero genes must yield zero BLAST matches. For each of the individual query genomes described in our results, regressions were run on the number of BLASTP hits returned from each major eukaryotic group ([Supplementary Table 2](#S1){ref-type="supplementary-material"}) versus the number of genes in each group's database. Studentized residuals were calculated in each case to determine whether significant outliers were present.

Fisher's exact tests were carried out using online software as described previously[@b15] to test specific predictions of the model of serial plastid evolution that emerged from regression analyses. Two control genomes were chosen for each test from among those contained within our large database using the following criteria; (1) they were from organisms with no inferred history of photosynthesis, (2) they had no inferred phylogenetic or endosymbiotic associations with any of the algal lineages tested and (3) they were the closest in size (smaller and larger, respectively), among the genomes in our database, to the algal genome tested. Significance of one-tailed *P* values for multiple tests were determined using a sequential Bonferroni test[@b43] with an *a priori ∝*=0.001.

Phylogenetic analyses
---------------------

A tree of plastid relationships was generated using inferred amino-acid sequences from the 11 largest genes retained in all plastid genomes from the 16 photosynthetic species included in our analyses of EGT in nuclear genomes. To avoid complications from additional experimental variables, we only analysed plastid genes from the photosynthetic organisms already present in our large database, so that results were directly comparable to nuclear gene analyses. Sequences were aligned using MUSCLE in MEGA 5.2.2 (ref. [@b44]). To minimize potential phylogenetic artifacts, we took a conservative approach, removing all regions with inferred gaps even if they appeared to align well otherwise. Tests for the most appropriate model of sequence evolution were performed for all individual alignments in MEGA 5.2.2, with a cpREV+G+I model inferred in all cases. A single concatenated alignment of 5,818 positions was created and used for tree reconstruction. Trees were inferred using maximum-likelihood in PhyML[@b45] with 1,000 bootstrap repetitions, and by Bayesian inference using MrBayes[@b46] run for 1 million generations with trees sampled every 100 generations. On the basis of standard deviations among split frequencies, the final 5,000 trees sampled were used to construct a Bayesian consensus tree and infer posterior probabilities. To determine effects (potential long-branch artifacts) from rooting the tree with the green plastid clade, we repeated the phylogenetic analyses (as described above) with a sub-alignment containing sequences only from red algal plastids and their chromist descendants.

Author contributions
====================

J.W.S. was involved in the overall study design and data interpretation, and was the primary author of the paper; J.S. performed statistical tests and phylogenetic analyses; J.Y. assembled the data sets and carried out BLAST analyses; H.G. wrote and implemented computational programs for parsing BLAST results; Q.D. was involved in design of the overall study and computational programs; J.H. participated in the overall study design and bioinformatic analyses. All the authors contributed editorially to the paper.

Additional information
======================

**How to cite this article:** Stiller, J. W. *et al*. The evolution of photosynthesis in chromist algae through serial endosymbioses. *Nat. Commun.* 5:5764 doi: 10.1038/ncomms6764 (2014).

Supplementary Material {#S1}
======================

###### Supplementary Information

Supplementary Figures 1-2, Supplementary Tables 1-2, and Supplementary References

###### Supplementary Data 1

Tables of top Blast matches among chromist groups used in the regressions depicted in Figure 1 and Supplementary Figure 2, and highlighted in Supplementary Table 1. These lists contain sequences responsible for the greater than expected similarities among chromist nuclear genomes; however, individual sequences should not be considered specific evidence of endosymbiotic gene transfer without explicit gene-by-gene analyses.

###### Supplementary Data 2

Java program files, along with explanatory readme file, contain scripts used to extract Blast matches used in statistical analyses. The Java program parses Blast output files, extracts the data needed and stores them in a database, so that various queries can be performed to assemble data for statistical analyses.

This research effort received support from NSF Research Collaboration Network grant 0741907, NSF Assembling the Tree of Life grant (DEB 0830024), NSFC Oversea, Hong Kong, Macao collaborative grant (31328003) and the CAS/SAFEA International Partnership Program for Creative Research Teams.

![Evidence of EGT from regression analyses with significant outliers.\
Linear regressions on relationships between the number of most similar sequences (measured as top BLAST matches) from 13 other eukaryotic lineages to all inferred protein-encoding genes in haptophyte (solid line, black triangles) and cryptophyte (dashed line, grey circles) genomes. The largest outliers are highlighted with studentized residuals (residual divided by its standard deviation; 's.d.' on figure) and images of the taxa that share more genes than expected based on the model, heterokonts in both cases. For comparison, the data points that indicate measures of similarity between cryptophytes and haptophytes genome also are highlighted. No studentized residual is provided for the data point indicating the number of top BLAST hits to haptophytes using the cryptophyte genome as query, because it falls at the predicted value from the regression line. Arrows indicate the proposed direction of EGT based on our overall model, except for the double-pointed arrow (the X indicating no evidence of EGT in either direction), which shows the reciprocal results between haptophytes and cryptophytes. A query's matches to its own group (for example, haptophyte to haptophyte) are not counted, meaning these reciprocal data points are both unpaired on the figure.](ncomms6764-f1){#f1}

![Model of serial plastid endosymbioses and a test using genes from plastid genomes.\
(**a**) Model of serial plastid endosymbioses suggested by regression analyses. The relationships depicted agree with general inferences from eukaryotic phylogenomics that heterokonts are not closely related to cryptophytes or haptophytes, and that each of the three groups emerge from mutually exclusive clades containing heterotrophic relatives. With respect to our proposed model of serial endosymbioses, the specific topology of the tree is not important, only that the three chromist algal groups do not form an exclusively monophyletic grouping that excludes aplastidial heterokonts. (**b**) Tree of plastid relationships based on an alignment of 5,818 amino-acid positions from genes inherited directly through the plastid genome, an independent data set for testing the model of plastid transfer inferred from EGT. The tree shown was recovered using both Bayesian and maximum-likelihood (ML) approaches. Bayesian posterior probabilities and ML bootstrap support values are provided for each node, with a star indicating 1.0 and 100% support, respectively. A second bootstrap value on nodes in the red plastid clade is from ML analyses performed in the absence of the green plastid lineage as outgroup. Bayesian probabilities are based on 5,000 sampled trees and bootstrap support values are from 1,000 replicates in each case.](ncomms6764-f2){#f2}

###### Fisher's exact tests for significant impacts of EGT from non-photosynthetic genes in heterokont and haptophyte genomes.

  **Query**                              **Top BLAST search matches to:**    **Ratio (adjusted to genome size) target:control**   ***P***   
  ----------------- -------------------- ----------------------------------- ---------------------------------------------------- --------- ----------
  *Phaeodactylum*   Photosynthetic       276                                 48[\*](#t1-fn1){ref-type="fn"}                       3.55       1.7e−20
                                                                             185[†](#t1-fn2){ref-type="fn"}                       2.13       3.4e−16
  *Aureococcus*     Photosynthetic       298                                 59[\*](#t1-fn1){ref-type="fn"}                       3.13       2.1e−19
                                                                             225[†](#t1-fn2){ref-type="fn"}                       1.9        1.7e−13
  *Phytophthora*    Non-photosynthetic   248                                 103[\*](#t1-fn1){ref-type="fn"}                      1.49        0.0003
                                                                             554[†](#t1-fn2){ref-type="fn"}                       0.62        0.999
  *Blastocystis*    Non-photosynthetic   75                                  50[\*](#t1-fn1){ref-type="fn"}                       0.93        0.695
                                                                             255[†](#t1-fn2){ref-type="fn"}                       0.44        0.999
  **Haptophyte**                         **Heterokonts**                                                                                         
  *Emiliania*       Photosynthetic       1,088[‡](#t1-fn3){ref-type="fn"}    436[§](#t1-fn4){ref-type="fn"}                       5.0        2.7e−204
                                                                             402[†](#t1-fn2){ref-type="fn"}                       2.0        8.3e−36
  *Emiliania*       Photosynthetic       658[\|\|](#t1-fn5){ref-type="fn"}   456[§](#t1-fn4){ref-type="fn"}                       3.85       1.8e−110
                                                                             396[†](#t1-fn2){ref-type="fn"}                       1.62       1.9e−15

^\*^Heterolobosea.

^†^Amoebozoa.

^‡^Matches to photosynthetic heterokonts.

^§^Fungi.

^\|\|^Matches to non-photosynthetic heterokonts.
