The genome of Spironucleus salmonicida highlights a fish pathogen adapted to fluctuating environments.

Spironucleus salmonicida causes systemic infections in salmonid fish. It belongs to the group diplomonads, binucleated heterotrophic flagellates adapted to micro-aerobic environments. Recently we identified energy-producing hydrogenosomes in S. salmonicida. Here we present a genome analysis of the fish parasite with a focus on the comparison to the more studied diplomonad Giardia intestinalis. We annotated 8067 protein coding genes in the ∼12.9 Mbp S. salmonicida genome. Unlike G. intestinalis, promoter-like motifs were found upstream of genes which are correlated with gene expression, suggesting a more elaborate transcriptional regulation. S. salmonicida can utilise more carbohydrates as energy sources, has an extended amino acid and sulfur metabolism, and more enzymes involved in scavenging of reactive oxygen species compared to G. intestinalis. Both genomes have large families of cysteine-rich membrane proteins. A cluster analysis indicated large divergence of these families in the two diplomonads. Nevertheless, one of S. salmonicida cysteine-rich proteins was localised to the plasma membrane similar to G. intestinalis variant-surface proteins. We identified S. salmonicida homologs to cyst wall proteins and showed that one of these is functional when expressed in Giardia. This suggests that the fish parasite is transmitted as a cyst between hosts. The extended metabolic repertoire and more extensive gene regulation compared to G. intestinalis suggest that the fish parasite is more adapted to cope with environmental fluctuations. Our genome analyses indicate that S. salmonicida is a well-adapted pathogen that can colonize different sites in the host.

On the reversibility of parasitism: adaptation to a free-living lifestyle via gene acquisitions in the diplomonad Trepomonas sp. PC1.

BACKGROUND: It is generally thought that the evolutionary transition to parasitism is irreversible because it is associated with the loss of functions needed for a free-living lifestyle. Nevertheless, free-living taxa are sometimes nested within parasite clades in phylogenetic trees, which could indicate that they are secondarily free-living. Herein, we test this hypothesis by studying the genomic basis for evolutionary transitions between lifestyles in diplomonads, a group of anaerobic eukaryotes. Most described diplomonads are intestinal parasites or commensals of various animals, but there are also free-living diplomonads found in oxygen-poor environments such as marine and freshwater sediments. All these nest well within groups of parasitic diplomonads in phylogenetic trees, suggesting that they could be secondarily free-living. RESULTS: We present a transcriptome study of Trepomonas sp. PC1, a diplomonad isolated from marine sediment. Analysis of the metabolic genes revealed a number of proteins involved in degradation of the bacterial membrane and cell wall, as well as an extended set of enzymes involved in carbohydrate degradation and nucleotide metabolism. Phylogenetic analyses showed that most of the differences in metabolic capacity between free-living Trepomonas and the parasitic diplomonads are due to recent acquisitions of bacterial genes via gene transfer. Interestingly, one of the acquired genes encodes a ribonucleotide reductase, which frees Trepomonas from the need to scavenge deoxyribonucleosides. The transcriptome included a gene encoding squalene-tetrahymanol cyclase. This enzyme synthesizes the sterol substitute tetrahymanol in the absence of oxygen, potentially allowing Trepomonas to thrive under anaerobic conditions as a free-living bacterivore, without depending on sterols from other eukaryotes. CONCLUSIONS: Our findings are consistent with the phylogenetic evidence that the last common ancestor of diplomonads was dependent on a host and that Trepomonas has adapted secondarily to a free-living lifestyle. We believe that similar studies of other groups where free-living taxa are nested within parasites could reveal more examples of secondarily free-living eukaryotes.

Gene transfer and diversification of microbial eukaryotes.

The importance of lateral gene transfer in genome evolution of microbial eukaryotes is slowly being appreciated. Acquisitions of genes have led to metabolic adaptation in diverse eukaryotic lineages. In most cases the metabolic genes have originated from prokaryotes, often followed by sequential transfers between eukaryotes. However, the knowledge of gene transfer in eukaryotes is still mainly based on anecdotal evidence. Some of the observed patterns may be biases in experimental approaches and sequence databases rather than evolutionary trends. Rigorous systematic studies of gene acquisitions that allow for the possibility of exchanges of all categories of genes from all sources are needed to get a more objective view of gene transfer in eukaryote evolution. It may be that the role of gene transfer in the diversification process of microbial eukaryotes currently is underestimated.

Genome-wide analyses of recombination suggest that Giardia intestinalis assemblages represent different species.

Giardia intestinalis is a major cause of waterborne enteric disease in humans. The species is divided into eight assemblages suggested to represent separate Giardia species based on host specificities and the genetic divergence of marker genes. We have investigated whether genome-wide recombination occurs between assemblages using the three available G. intestinalis genomes. First, the relative nonsynonymous substitution rates of the homologs were compared for 4,009 positional homologs. The vast majority of these comparisons indicate genetic isolation without interassemblage recombinations. Only a region of 6 kbp suggests genetic exchange between assemblages A and E, followed by gene conversion events. Second, recombination-detecting software fails to identify within-gene recombination between the different assemblages for most of the homologs. Our results indicate very low frequency of recombination between the syntenic core genes, suggesting that G. intestinalis assemblages are genetically isolated lineages and thus should be viewed as separated Giardia species.

CD4⁺ T-cell immunity in the female genital tract is critically dependent on local mucosal immunization.

Immunizations via the i.n. and intravaginal (ivag) routes effectively generate strong genital tract antibody-mediated immunity. To what extent the same is true for T-cell responses is incompletely known. Therefore, we set out to investigate optimal conditions for stimulation of genital tract CD4(+) T-cell responses, using adoptive transfer of mouse DO11.10 TCR transgenic T cells specific for OVA and OVA conjugated to cholera toxin (CT) as an immunogen. We observed that progesterone was required for a T-cell response following ivag immunization, whereas estradiol prevented a response. Although i.n. immunization stimulated OVA-specific CD4(+) T-cell responses in the draining LNs, it was substantially less effective compared to ivag. More importantly, an ivag booster immunization was absolutely required to attract T cells to the genital tract mucosa itself. While clinical use of CT is precluded because of its toxicity, we developed a combined adjuvant vector based on a non-toxic derivative of CT and immune-stimulating complexes. The CTA1-DD/immune-stimulating complexes (ISCOMs) adjuvant together with major outer membrane protein was effective at stimulating genital tract CD4(+) T-cell immunity and protection against a live chlamydial infection, which holds promise for the development of mucosal vaccines against sexually transmitted infections.

A chromosome-scale reference genome for Spironucleus salmonicida.

Spironucleus salmonicida is a diplomonad causing systemic infection in salmon. The first S. salmonicida genome assembly was published 2014 and has been a valuable reference genome in protist research. However, the genome assembly is fragmented without assignment of the sequences to chromosomes. In our previous Giardia genome study, we have shown how a fragmented genome assembly can be improved with long-read sequencing technology complemented with optical maps. Combining Pacbio long-read sequencing technology and optical maps, we are presenting here this new S. salmonicida genome assembly in nine near-complete chromosomes with only three internal gaps at long repeats. This new genome assembly is not only more complete sequence-wise but also more complete at annotation level, providing more details into gene families, gene organizations and chromosomal structure. This near-complete reference genome will aid comparative genomics at chromosomal level, and serve as a valuable resource for the diplomonad community and protist research.

Draft genome sequencing of giardia intestinalis assemblage B isolate GS: is human giardiasis caused by two different species?

Giardia intestinalis is a major cause of diarrheal disease worldwide and two major Giardia genotypes, assemblages A and B, infect humans. The genome of assemblage A parasite WB was recently sequenced, and the structurally compact 11.7 Mbp genome contains simplified basic cellular machineries and metabolism. We here performed 454 sequencing to 16x coverage of the assemblage B isolate GS, the only Giardia isolate successfully used to experimentally infect animals and humans. The two genomes show 77% nucleotide and 78% amino-acid identity in protein coding regions. Comparative analysis identified 28 unique GS and 3 unique WB protein coding genes, and the variable surface protein (VSP) repertoires of the two isolates are completely different. The promoters of several enzymes involved in the synthesis of the cyst-wall lack binding sites for encystation-specific transcription factors in GS. Several synteny-breaks were detected and verified. The tetraploid GS genome shows higher levels of overall allelic sequence polymorphism (0.5 versus <0.01% in WB). The genomic differences between WB and GS may explain some of the observed biological and clinical differences between the two isolates, and it suggests that assemblage A and B Giardia can be two different species.

Genome analysis and comparative genomics of a Giardia intestinalis assemblage E isolate.

BACKGROUND: Giardia intestinalis is a protozoan parasite that causes diarrhea in a wide range of mammalian species. To further understand the genetic diversity between the Giardia intestinalis species, we have performed genome sequencing and analysis of a wild-type Giardia intestinalis sample from the assemblage E group, isolated from a pig. RESULTS: We identified 5012 protein coding genes, the majority of which are conserved compared to the previously sequenced genomes of the WB and GS strains in terms of microsynteny and sequence identity. Despite this, there is an unexpectedly large number of chromosomal rearrangements and several smaller structural changes that are present in all chromosomes. Novel members of the VSP, NEK Kinase and HCMP gene families were identified, which may reveal possible mechanisms for host specificity and new avenues for antigenic variation. We used comparative genomics of the three diverse Giardia intestinalis isolates P15, GS and WB to define a core proteome for this species complex and to identify lineage-specific genes. Extensive analyses of polymorphisms in the core proteome of Giardia revealed differential rates of divergence among cellular processes. CONCLUSIONS: Our results indicate that despite a well conserved core of genes there is significant genome variation between Giardia isolates, both in terms of gene content, gene polymorphisms, structural chromosomal variations and surface molecule repertoires. This study improves the annotation of the Giardia genomes and enables the identification of functionally important variation.

Large genomic differences between the morphologically indistinguishable diplomonads Spironucleus barkhanus and Spironucleus salmonicida.

BACKGROUND: Microbial eukaryotes show large variations in genome structure and content between lineages, indicating extensive flexibility over evolutionary timescales. Here we address the tempo and mode of such changes within diplomonads, flagellated protists with two nuclei found in oxygen-poor environments. Approximately 5,000 expressed sequence tag (EST) sequences were generated from the fish commensal Spironucleus barkhanus and compared to sequences from the morphologically indistinguishable fish parasite Spironucleus salmonicida, and other diplomonads. The ESTs were complemented with sequence variation studies in selected genes and genome size determinations. RESULTS: Many genes detected in S. barkhanus and S. salmonicida are absent in the human parasite Giardia intestinalis, the most intensively studied diplomonad. For example, these fish diplomonads show an extended metabolic repertoire and are able to incorporate selenocysteine into proteins. The codon usage is altered in S. barkhanus compared to S. salmonicida. Sequence variations were found between individual S. barkhanus ESTs for many, but not all, protein coding genes. Conversely, no allelic variation was found in a previous genome survey of S. salmonicida. This difference was confirmed by sequencing of genomic DNA. Up to five alleles were identified for the cloned S. barkhanus genes, and at least nineteen highly expressed S. barkhanus genes are represented by more than four alleles in the EST dataset. This could be explained by the presence of a non-clonal S. barkhanus population in the culture, by a ploidy above four, or by duplications of parts of the genome. Indeed, genome size estimations using flow cytometry indicated similar haploid genome sizes in S. salmonicida and G. intestinalis (approximately 12 Mb), whereas the S. barkhanus genome is larger (approximately 18 Mb). CONCLUSIONS: This study indicates extensive divergent genome evolution within diplomonads. Genomic traits such as codon usage, frequency of allelic sequence variation, and genome size have changed considerably between S. barkhanus and S. salmonicida. These observations suggest that large genomic differences may accumulate in morphologically indistinguishable eukaryotic microbes.

Metabolic Reconstruction Elucidates the Lifestyle of the Last Diplomonadida Common Ancestor.

The identification of ancestral traits is essential to understanding the evolution of any group. In the case of parasitic groups, this helps us understand the adaptation to this lifestyle and a particular host. Most diplomonads are parasites, but there are free-living members of the group nested among the host-associated diplomonads. Furthermore, most of the close relatives within Fornicata are free-living organisms. This leaves the lifestyle of the ancestor unclear. Here, we present metabolic maps of four different diplomonad species. We identified 853 metabolic reactions and 147 pathways present in at least one of the analyzed diplomonads. Our study suggests that diplomonads represent a metabolically diverse group in which differences correlate with different environments (e.g., the detoxification of arsenic). Using a parsimonious analysis, we also provide a description of the putative metabolism of the last Diplomonadida common ancestor. Our results show that the acquisition and loss of reactions have shaped metabolism since this common ancestor. There is a net loss of reaction in all branches leading to parasitic diplomonads, suggesting an ongoing reduction in the metabolic capacity. Important traits present in host-associated diplomonads (e.g., virulence factors and the synthesis of UDP-N-acetyl-d-galactosamine) are shared with free-living relatives. The last Diplomonadida common ancestor most likely already had acquired important enzymes for the salvage of nucleotides and had a reduced capacity to synthesize nucleotides, lipids, and amino acids de novo, suggesting that it was an obligate host-associated organism.IMPORTANCE Diplomonads are a group of microbial eukaryotes found in oxygen-poor environments. There are both parasitic (e.g., Giardia intestinalis) and free-living (e.g., Trepomonas) members in the group. Diplomonads are well known for their anaerobic metabolism, which has been studied for many years. Here, we reconstructed whole metabolic networks of four extant diplomonad species as well as their ancestors, using a bioinformatics approach. We show that the metabolism within the group is under constant change throughout evolutionary time, in response to the environments that the different lineages explore. Both gene losses and gains are responsible for the adaptation processes. Interestingly, it appears that the last Diplomonadida common ancestor had a metabolism that is more similar to extant parasitic than free-living diplomonads. This suggests that the host-associated lifestyle of parasitic diplomonads, such as the human parasite G. intestinalis, is an old evolutionary adaptation.

The compact genome of Giardia muris reveals important steps in the evolution of intestinal protozoan parasites.

Diplomonad parasites of the genus Giardia have adapted to colonizing different hosts, most notably the intestinal tract of mammals. The human-pathogenic Giardia species, Giardia intestinalis, has been extensively studied at the genome and gene expression level, but no such information is available for other Giardia species. Comparative data would be particularly valuable for Giardia muris, which colonizes mice and is commonly used as a prototypic in vivo model for investigating host responses to intestinal parasitic infection. Here we report the draft-genome of G. muris. We discovered a highly streamlined genome, amongst the most densely encoded ever described for a nuclear eukaryotic genome. G. muris and G. intestinalis share many known or predicted virulence factors, including cysteine proteases and a large repertoire of cysteine-rich surface proteins involved in antigenic variation. Different to G. intestinalis, G. muris maintains tandem arrays of pseudogenized surface antigens at the telomeres, whereas intact surface antigens are present centrally in the chromosomes. The two classes of surface antigens engage in genetic exchange. Reconstruction of metabolic pathways from the G. muris genome suggest significant metabolic differences to G. intestinalis. Additionally, G. muris encodes proteins that might be used to modulate the prokaryotic microbiota. The responsible genes have been introduced in the Giardia genus via lateral gene transfer from prokaryotic sources. Our findings point to important evolutionary steps in the Giardia genus as it adapted to different hosts and it provides a powerful foundation for mechanistic exploration of host-pathogen interaction in the G. muris-mouse pathosystem.

Convergent evolution: gene sharing by eukaryotic plant pathogens.

Oomycetes and filamentous parasitic fungi are plant pathogens that have undergone convergent evolution. A recent study has shown that these microbial eukaryotes have exchanged metabolic genes, which might explain some of their phenotypic similarities.

Horizontal gene transfer between microbial eukaryotes.

Comparative genomics have identified two loosely defined classes of genes: widely distributed core genes that encode proteins for central functions in the cell and accessory genes that are patchily distributed across lineages and encode taxa-specific functions. Studies of microbial eukaryotes show that both categories undergo horizontal gene transfer (HGT) from prokaryotes, but also between eukaryotic organisms. Intra-domain gene transfers of most core genes seem to be relatively infrequent and therefore comparatively easy to detect using phylogenetic methods. In contrast, phylogenies of accessory genes often have complex topologies with little or no resemblance of organismal relationships typically with eukaryotes and prokaryotes intermingled, making detailed evolutionary histories difficult to interpret. Nevertheless, this suggests significant rates of gene transfer between and among the three domains of life for many of these genes, affecting a considerably diversity of eukaryotic microbes, although the current depth of taxonomic sampling usually is insufficient to pin down individual transfer events. The occurrence of intra-domain transfer among microbial eukaryotes has important implications for studies of organismal phylogeny as well as eukaryote genome evolution in general.

Lateral Acquisitions Repeatedly Remodel the Oxygen Detoxification Pathway in Diplomonads and Relatives.

Oxygen and reactive oxygen species (ROS) are important stress factors for cells because they can oxidize many large molecules. Fornicata, a group of flagellated protists that includes diplomonads, have anaerobic metabolism but are still able to tolerate fluctuating levels of oxygen. We identified 25 protein families putatively involved in detoxification of oxygen and ROS in this group using a bioinformatics approach and propose how these interact in an oxygen detoxification pathway. These protein families were divided into a central oxygen detoxification pathway and accessory pathways for the synthesis of nonprotein thiols. We then used a phylogenetic approach to investigate the evolutionary origin of the components of this putative pathway in Diplomonadida and other Fornicata species. Our analyses suggested that the diplomonad ancestor was adapted to low-oxygen levels, was able to reduce O2 to H2O in a manner similar to extant diplomonads, and was able to synthesize glutathione and l-cysteine. Several genes involved in the pathway have complex evolutionary histories and have apparently been repeatedly acquired through lateral gene transfer and subsequently lost. At least seven genes were acquired independently in different Fornicata lineages, leading to evolutionary convergences. It is likely that acquiring these oxygen detoxification proteins helped anaerobic organisms (like the parasitic Giardia intestinalis) adapt to low-oxygen environments (such as the digestive tract of aerobic hosts).

Evolution of the cutinase gene family: evidence for lateral gene transfer of a candidate Phytophthora virulence factor.

Lateral gene transfer (LGT) can facilitate the acquisition of new functions in recipient lineages, which may enable them to colonize new environments. Several recent publications have shown that gene transfer between prokaryotes and eukaryotes occurs with appreciable frequency. Here we present a study of interdomain gene transfer of cutinases -- well documented virulence factors in fungi -- between eukaryotic plant pathogens Phytophthora species and prokaryotic bacterial lineages. Two putative cutinase genes were cloned from Phytophthora brassicae and Northern blotting experiments showed that these genes are expressed early during the infection of the host Arabidopsis thaliana and induced during cyst germination of the pathogen. Analysis of the gene organisation of this gene family in Phytophthora ramorum and P. sojae showed three and ten copies in tight succession within a region of 5 and 25 kb, respectively, probably indicating a recent expansion in Phytophthora lineages by gene duplications. Bioinformatic analyses identified orthologues only in three genera of Actinobacteria, and in two distantly related eukaryotic groups: oomycetes and fungi. Together with phylogenetic analyses this limited distribution of the gene in the tree of life strongly support a scenario where cutinase genes originated after the origin of land plants in a microbial lineage living in proximity of plants and subsequently were transferred between distantly related plant-degrading microbes. More precisely, a cutinase gene was likely acquired by an ancestor of P. brassicae, P. sojae, P. infestans and P. ramorum, possibly from an actinobacterial source, suggesting that gene transfer might be an important mechanism in the evolution of their virulence. These findings could indeed provide an interesting model system to study acquisition of virulence factors in these important plant pathogens.

Synchronisation of Giardia lamblia: identification of cell cycle stage-specific genes and a differentiation restriction point.

The intestinal parasite Giardia lamblia undergoes cell differentiations that entail entry into and departure from the replicative cell cycle. The pathophysiology of giardiasis depends directly upon the ability of the trophozoite form to replicate in the host upper small intestine. Thus, cell proliferation is tightly linked to disease. However, studies of cell cycle regulation in Giardia have been hampered by the inability to synchronise cultures. Here we report that Giardia isolates of the major human genotypes A and B can be synchronised using aphidicolin, a mycotoxin that reversibly inhibits replicative DNA polymerases in eukaryotic cells. Aphidicolin arrests Giardia trophozoites in the early DNA synthesis (S) phase of the cell cycle. We identified a set of cell cycle orthologues in the Giardia genome using bioinformatic analyses and showed that synchronised parasites express these genes in a cell cycle stage-specific manner. The synchronisation method also showed that during encystation, exit from the ordinary cell cycle occurs preferentially in G(2) and defines a restriction point for differentiation. Synchronisation opens up possibilities for further molecular and cell biological studies of chromosome replication, mitosis and segregation of the complex cytoskeleton in Giardia.

Dominance of Giardia assemblage B in León, Nicaragua.

Giardiasis is a major problem in León, Nicaragua, yet despite this no data are available regarding the prevalence of different Giardia genotypes in this area. To address this question, a molecular analysis of Giardia isolates from humans and dogs living in the same area in León, Nicaragua was performed. Giardia isolates from 119 Nicaraguan patients and 8 dogs were successfully genotyped using single and/or nested beta-giardin PCR with subsequent restriction length fragment polymorphism (RFLP) analysis. The analyses of human samples yielded 94 (79%) assemblage B isolates and 25 (21%) assemblage A isolates. Only the non-human-associated assemblages C and D were found in the dog samples. Sixteen isolates with assemblage A pattern, 26 isolates with assemblage B pattern and all dog isolates were further characterized by sequencing the nested beta-giardin PCR product and by molecular analyses of the glutamate dehydrogenase (gdh) gene. Within the study area the assemblage A isolates were highly genetically homogenous, showing only sub-genotypes A2 (n=3) or A3 (n=13) at the beta-giardin locus and AII only at the gdh locus while assemblage B showed a high genetic polymorphism at both loci. Seven different sub-genotypes were identified within 13 of the sequenced assemblage B beta-giardin isolates. The remaining 13 sequenced assemblage B-isolates appeared to contain several different variants of the beta-giardin gene since the chromatograms displayed one to seven double peaks. The gdh sequences showed an even higher polymorphism since only 2 of 26 assemblage B isolates were without double peaks. Two mixed infections between assemblage A and B were found when the gdh gene was analyzed. Polymorphisms were also observed in the dog-associated assemblages C and D, but to a lesser extent than in assemblage B.

A genomic survey of the fish parasite Spironucleus salmonicida indicates genomic plasticity among diplomonads and significant lateral gene transfer in eukaryote genome evolution.

BACKGROUND: Comparative genomic studies of the mitochondrion-lacking protist group Diplomonadida (diplomonads) has been lacking, although Giardia lamblia has been intensively studied. We have performed a sequence survey project resulting in 2341 expressed sequence tags (EST) corresponding to 853 unique clones, 5275 genome survey sequences (GSS), and eleven finished contigs from the diplomonad fish parasite Spironucleus salmonicida (previously described as S. barkhanus). RESULTS: The analyses revealed a compact genome with few, if any, introns and very short 3' untranslated regions. Strikingly different patterns of codon usage were observed in genes corresponding to frequently sampled ESTs versus genes poorly sampled, indicating that translational selection is influencing the codon usage of highly expressed genes. Rigorous phylogenomic analyses identified 84 genes--mostly encoding metabolic proteins--that have been acquired by diplomonads or their relatively close ancestors via lateral gene transfer (LGT). Although most acquisitions were from prokaryotes, more than a dozen represent likely transfers of genes between eukaryotic lineages. Many genes that provide novel insights into the genetic basis of the biology and pathogenicity of this parasitic protist were identified including 149 that putatively encode variant-surface cysteine-rich proteins which are candidate virulence factors. A number of genomic properties that distinguish S. salmonicida from its human parasitic relative G. lamblia were identified such as nineteen putative lineage-specific gene acquisitions, distinct mutational biases and codon usage and distinct polyadenylation signals. CONCLUSION: Our results highlight the power of comparative genomic studies to yield insights into the biology of parasitic protists and the evolution of their genomes, and suggest that genetic exchange between distantly-related protist lineages may be occurring at an appreciable rate in eukaryote genome evolution.

A cyanobacterial gene in nonphotosynthetic protists--an early chloroplast acquisition in eukaryotes?

Since the incorporation of mitochondria and chloroplasts (plastids) into the eukaryotic cell by endosymbiosis, genes have been transferred from the organellar genomes to the nucleus of the host, via an ongoing process known as endosymbiotic gene transfer. Accordingly, in photosynthetic eukaryotes, nuclear genes with cyanobacterial affinity are believed to have originated from endosymbiotic gene transfer from chloroplasts. Analysis of the Arabidopsis thaliana genome has shown that a significant fraction (2%-9%) of the nuclear genes have such an endosymbiotic origin. Recently, it was argued that 6-phosphogluconate dehydrogenase (gnd)-the second enzyme in the oxidative pentose phosphate pathway-was one such example. Here we show that gnd genes with cyanobacterial affinity also are present in several nonphotosynthetic protistan lineages, such as Heterolobosea, Apicomplexa, and parasitic Heterokonta. Current data cannot definitively resolve whether these groups acquired the gnd gene by primary and/or secondary endosymbiosis or via an independent lateral gene transfer event. Nevertheless, our data suggest that chloroplasts were introduced into eukaryotes much earlier than previously thought and that several major groups of heterotrophic eukaryotes have secondarily lost photosynthetic plastids.