Evidence for Transitional Stages in the Evolution of Euglenid Group II Introns and Twintrons in the Monomorphina aenigmatica Plastid Genome

Background

Photosynthetic euglenids acquired their plastid by secondary endosymbiosis of a prasinophyte-like green alga. But unlike its prasinophyte counterparts, the plastid genome of the euglenid Euglena gracilis is riddled with introns that interrupt almost every protein-encoding gene. The atypical group II introns and twintrons (introns-within-introns) found in the E. gracilis plastid have been hypothesized to have been acquired late in the evolution of euglenids, implying that massive numbers of introns may be lacking in other taxa. This late emergence was recently corroborated by the plastid genome sequences of the two basal euglenids, Eutreptiella gymnastica and Eutreptia viridis, which were found to contain fewer introns.

Methodology/Principal Findings

To gain further insights into the proliferation of introns in euglenid plastids, we have characterized the complete plastid genome sequence of Monomorphina aenigmatica, a freshwater species occupying an intermediate phylogenetic position between early and late branching euglenids. The M. aenigmatica UTEX 1284 plastid genome (74,746 bp, 70.6% A+T, 87 genes) contains 53 intron insertion sites, of which 41 were found to be shared with other euglenids including 12 of the 15 twintron insertion sites reported in E. gracilis.

Conclusions

The pattern of insertion sites suggests an ongoing but uneven process of intron gain in the lineage, with perhaps a minimum of two bursts of rapid intron proliferation. We also identified several sites that represent intermediates in the process of twintron evolution, where the external intron is in place, but not the internal one, offering a glimpse into how these convoluted molecular contraptions originate.     

The complete chloroplast genome of the chlorarachniophyte Bigelowiella natans: evidence for independent origins of chlorarachniophyte and euglenid secondary endosymbionts

Chlorarachniophytes are amoeboflagellate cercozoans that acquired a plastid by secondary endosymbiosis. Chlorarachniophytes are the last major group of algae for which there is no completely sequenced plastid genome. Here we describe the 69.2-kbp chloroplast genome of the model chlorarachniophyte Bigelowiella natans. The genome is highly reduced in size compared with plastids of other photosynthetic algae and is closer in size to genomes of several nonphotosynthetic plastids. Unlike nonphotosynthetic plastids, however, the B. natans chloroplast genome has not sustained a massive loss of genes, and it retains nearly all of the functional photosynthesis-related genes represented in the genomes of other green algae. Instead, the genome is highly compacted and gene dense. The genes are organized with a strong strand bias, and several unusual rearrangements and inversions also characterize the genome; notably, an inversion in the small-subunit rRNA gene, a translocation of 3 genes in the major ribosomal protein operon, and the fragmentation of the cluster encoding the large photosystem proteins PsaA and PsaB. The chloroplast endosymbiont is known to be a green alga, but its evolutionary origin and relationship to other primary and secondary green plastids has been much debated. A recent hypothesis proposes that the endosymbionts of chlorarachniophytes and euglenids share a common origin (the Cabozoa hypothesis). We inferred phylogenies using individual and concatenated gene sequences for all genes in the genome. Concatenated gene phylogenies show a relationship between the B. natans plastid and the ulvophyte-trebouxiophyte-chlorophyte clade of green algae to the exclusion of Euglena. The B. natans plastid is thus not closely related to that of Euglena, which suggests that plastids originated independently in these 2 groups and the Cabozoa hypothesis is false.     

Phylogeny of phagotrophic euglenids (Euglenozoa) as inferred from hsp90 gene sequences

Molecular phylogenies of euglenids are usually based on ribosomal RNA genes that do not resolve the branching order among the deeper lineages. We addressed deep euglenid phylogeny using the cytosolic form of the heat-shock protein 90 gene (hsp90), which has already been employed with some success in other groups of euglenozoans and eukaryotes in general. Hsp90 sequences were generated from three taxa of euglenids representing different degrees of ultrastructural complexity, namely Petalomonas cantuscygni and wild isolates of Entosiphon sulcatum, and Peranema trichophorum. The hsp90 gene sequence of P. trichophorum contained three short introns (ranging from 27 to 31 bp), two of which had non-canonical borders GG-GG and GG-TG and two 10-bp inverted repeats, suggesting a structure similar to that of the non-canonical introns described in Euglena gracilis. Phylogenetic analyses confirmed a closer relationship between kinetoplastids and diplonemids than to euglenids, and supported previous views regarding the branching order among primarily bacteriovorous, primarily eukaryovorous, and photosynthetic euglenids. The position of P. cantuscygni within Euglenozoa, as well as the relative support for the nodes including it were strongly dependent on outgroup selection. The results were most consistent when the jakobid Reclinomonas americana was used as the outgroup. The most robust phylogenies place P. cantuscygni as the most basal branch within the euglenid clade. However, the presence of a kinetoplast-like mitochondrial inclusion in P. cantuscygni deviates from the currently accepted apomorphy-based definition of the kinetoplastid clade and highlights the necessity of detailed studies addressing the molecular nature of the euglenid and diplonemid mitochondrial genome.     

Evolutionary History of the Enzymes Involved in the Calvin–Benson Cycle in Euglenids

Euglenids are an ancient lineage that may have existed as early as 2 billion years ago. A mere 65 years ago, Melvin Calvin and Andrew A. Benson performed experiments on Euglena gracilis and elucidated the series of reactions by which carbon was fixed and reduced during photosynthesis. However, the evolutionary history of this pathway (Calvin–Benson cycle) in euglenids was more complex than Calvin and Benson could have imagined. The chloroplast present today in euglenophytes arose from a secondary endosymbiosis between a phagotrophic euglenid and a prasinophyte green alga. A long period of evolutionary time existed before this secondary endosymbiotic event took place, which allowed for other endosymbiotic events or gene transfers to occur prior to the establishment of the green chloroplast. This research revealed the evolutionary history of the major enzymes of the Calvin–Benson cycle throughout the euglenid lineage and showed that the majority of genes for Calvin–Benson cycle enzymes shared an ancestry with red algae and/or chromophytes suggesting they may have been transferred to the nucleus prior to the acquisition of the green chloroplast.     

MOLECULAR EVOLUTION OF EUGLENOZOAN PARAXONEMAL ROD GENES par1 AND par2 COINCIDES WITH PHYLOGENETIC RECONSTRUCTION BASED ON SMALL SUBUNIT rDNA DATA1

Emergent flagella of Euglenozoa consist of two prominent structural elements: the axoneme built by microtubules with motor proteins to enable the movement of the flagellum and a highly organized protein structure of unknown function, called the paraxonemal rod (PAR), which consists of two major proteins paralleling the axoneme of euglenid and kinetoplastid emergent flagella. These flagellar structures are considered apomorphic characters of Euglenozoa. We examined the evolution of the genes par1 and par2 encoding the two major proteins, where we could show that these proteins are encoded by two very similar genes found in kinetoplastids and euglenids. The branching pattern indicated a gene duplication before the diversification into euglenids and kinetoplastids. In the clades of the genes, subtrees of euglenid and kinetoplastid monophyla arose. Both genes showed strong genetic diversity with biased GC content at taxon rather than at gene level. We also examined phylogenies inferred from PAR genes that are well in agreement with established small subunit rDNA analyses. Both showed further separation of the euglenid subtree into primary osmotrophs and a phototrophic clade, including secondarily derived osmotrophs.     

Functional gene losses occur with minimal size reduction in the plastid genome of the parasitic liverwort Aneura mirabilis

Aneura mirabilis is a parasitic liverwort that exploits an existing mycorrhizal association between a basidiomycete and a host tree. This unusual liverwort is the only known parasitic seedless land plant with a completely nonphotosynthetic life history. The complete plastid genome of A. mirabilis was sequenced to examine the effect of its nonphotosynthetic life history on plastid genome content. Using a partial genomic fosmid library approach, the genome was sequenced and shown to be 108,007 bp with a structure typical of green plant plastids. Comparisons were made with the plastid genome of Marchantia polymorpha, the only other liverwort plastid sequence available. All ndh genes are either absent or pseudogenes. Five of 15 psb genes are pseudogenes, as are 2 of 6 psa genes and 2 of 6 pet genes. Pseudogenes of cysA, cysT, ccsA, and ycf3 were also detected. The remaining complement of genes present in M. polymorpha is present in the plastid of A. mirabilis with intact open reading frames. All pseudogenes and gene losses co-occur with losses detected in the plastid of the parasitic angiosperm Epifagus virginiana, though the latter has functional gene losses not found in A. mirabilis. The plastid genome sequence of A. mirabilis represents only the second liverwort, and first mycoheterotroph, to have its plastid genome sequenced. We observed a pattern of genome evolution congruent with functional gene losses in parasitic angiosperms but suggest that its plastid genome represents a genome in the early stages of decay following the relaxation of selection pressures.     

Ins and outs of plastid genome evolution

Recent findings have established cracks in the straight-laced image of the plastid genome as a molecule whose sole function is photosynthesis and whose gene content is highly conserved. Genes for numerous non-photosynthetic functions have been identified. Algal plastid genomes contain many genes with no homologs in angiosperms, and the recent transfer of genes from the plastid to the nuclear genome has been described. Wholesale abandonment of genes encoding photosynthetic and gene-expression functions has occurred in the plastid genomes of a non-green plant and alga. The origins of plastid DNA, its use in phylogenetic studies, and the origins of plastid introns are also reviewed.     

Complete plastid genome sequence of the chickpea (Cicer arietinum) and the phylogenetic distribution of rps12 and clpP intron losses among legumes (Leguminosae)

Chickpea (Cicerarietinum, Leguminosae), an important grain legume, is widely used for food and fodder throughout the world. We sequenced the complete plastid genome of chickpea, which is 125,319bp in size, and contains only one copy of the inverted repeat (IR). The genome encodes 108 genes, including 4 rRNAs, 29 tRNAs, and 75 proteins. The genes rps16, infA, and ycf4 are absent in the chickpea plastid genome, and ndhB has an internal stop codon in the 5'exon, similar to other legumes. Two genes have lost their introns, one in the 3'exon of the transpliced gene rps12, and the one between exons 1 and 2 of clpP; this represents the first documented case of the loss of introns from both of these genes in the same plastid genome. An extensive phylogenetic survey of these intron losses was performed on 302 taxa across legumes and the related family Polygalaceae. The clpP intron has been lost exclusively in taxa from the temperate "IR-lacking clade" (IRLC), whereas the rps12 intron has been lost in most members of the IRLC (with the exception of Wisteria, Callerya, Afgekia, and certain species of Millettia, which represent the earliest diverging lineages of this clade), and in the tribe Desmodieae, which is closely related to the tribes Phaseoleae and Psoraleeae. Data provided here suggest that the loss of the rps12 intron occurred after the loss of the IR. The two new genomic changes identified in the present study provide additional support of the monophyly of the IR-loss clade, and resolution of the pattern of the earliest-branching lineages in this clade. The availability of the complete chickpea plastid genome sequence also provides valuable information on intergenic spacer regions among legumes and endogenous regulatory sequences for plastid genetic engineering.     

Complete sequence and analysis of the plastid genome of the unicellular red alga Cyanidioschyzon merolae.

The complete nucleotide sequence of the plastid genome of the unicellular primitive red alga Cyanidioschyzon merolae 10D (Cyanidiophyceae) was determined. The genome is a circular DNA composed of 149,987 bp with no inverted repeats. The G + C content of this plastid genome is 37.6%. The C. merolae plastid genome contains 243 genes, which are distributed on both strands and consist of 36 RNA genes (3 rRNAs, 31 tRNAs, tmRNA, and a ribonuclease P RNA component) and 207 protein genes, including unidentified open reading frames. The striking feature of this genome is the high degree of gene compaction; it has very short intergenic distances (approximately 40% of the protein genes were overlapped) and no genes have introns. This genome encodes several genes that are rarely found in other plastid genomes. A gene encoding a subunit of sulfate transporter (cysW) is the first to be identified in a plastid genome. The cysT and cysW genes are located in the C. merolae plastid genome in series, and they probably function together with other nuclear-encoded components of the sulfate transport system. Our phylogenetic results suggest that the Cyanidiophyceae, including C. merolae, are a basal clade within the red lineage plastids.     

Evolution of the Chloroplast Genome in Photosynthetic Euglenoids: A Comparison of Eutreptia viridis and Euglena gracilis (Euglenophyta)

The chloroplast genome of Eutreptia viridis Perty, a basal taxon in the photosynthetic euglenoid lineage, was sequenced and compared with that of Euglena gracilis Ehrenberg, a crown species. Several common gene clusters were identified and gene order, conservation, and sequence similarity was assessed through comparisons with Euglena gracilis. Significant gene rearrangements were present between Eutreptia viridis and Euglena gracilis chloroplast genomes. In addition, major expansion has occurred in the Euglena gracilis chloroplast accounting for its larger size. However, the key chloroplast genes are present and differ only in the absence of psaM and roaA in Eutreptia viridis, and psaI in Euglena gracilis, suggesting a high level of gene conservation within the euglenoid lineage. Further comparisons with the plastid genomes of closely related green algal taxa have provided additional support for the hypothesis that a Pyramimonas-like alga was the euglenoid chloroplast donor via secondary endosymbiosis.     

The evolution of the plastid chromosome in land plants: gene content, gene order, gene function

This review bridges functional and evolutionary aspects of plastid chromosome architecture in land plants and their putative ancestors. We provide an overview on the structure and composition of the plastid genome of land plants as well as the functions of its genes in an explicit phylogenetic and evolutionary context. We will discuss the architecture of land plant plastid chromosomes, including gene content and synteny across land plants. Moreover, we will explore the functions and roles of plastid encoded genes in metabolism and their evolutionary importance regarding gene retention and conservation. We suggest that the slow mode at which the plastome typically evolves is likely to be influenced by a combination of different molecular mechanisms. These include the organization of plastid genes in operons, the usually uniparental mode of plastid inheritance, the activity of highly effective repair mechanisms as well as the rarity of plastid fusion. Nevertheless, structurally rearranged plastomes can be found in several unrelated lineages (e.g. ferns, Pinaceae, multiple angiosperm families). Rearrangements and gene losses seem to correlate with an unusual mode of plastid transmission, abundance of repeats, or a heterotrophic lifestyle (parasites or myco-heterotrophs). While only a few functional gene gains and more frequent gene losses have been inferred for land plants, the plastid Ndh complex is one example of multiple independent gene losses and will be discussed in detail. Patterns of ndh-gene loss and functional analyses indicate that these losses are usually found in plant groups with a certain degree of heterotrophy, might rendering plastid encoded Ndh1 subunits dispensable.     

Migration of the plastid genome to the nucleus in a peridinin dinoflagellate

Dinoflagellate algae are important primary producers and of significant ecological and economic impact because of their ability to form "red tides". They are also models for evolutionary research because of an unparalleled ability to capture photosynthetic organelles (plastids) through endosymbiosis. The nature and extent of the plastid genome in the dominant perdinin-containing dinoflagellates remain, however, two of the most intriguing issues in plastid evolution. The plastid genome in these taxa is reduced to single-gene minicircles encoding an incomplete (until now 15) set of plastid proteins. The location of the remaining photosynthetic genes is unknown. We generated a data set of 6,480 unique expressed sequence tags (ESTs) from the toxic dinoflagellate Alexandrium tamarense (for details, see the Experimental Procedures in the Supplemental Data) to find the missing plastid genes and to understand the impact of endosymbiosis on genome evolution. Here we identify 48 of the non-minicircle-encoded photosynthetic genes in the nuclear genome of A. tamarense, accounting for the majority of the photosystem. Fifteen genes that are always found on the plastid genome of other algae and plants have been transferred to the nucleus in A. tamarense. The plastid-targeted genes have red and green algal origins. These results highlight the unique position of dinoflagellates as the champions of plastid gene transfer to the nucleus among photosynthetic eukaryotes.     

Genome evolution of a tertiary dinoflagellate plastid

The dinoflagellates have repeatedly replaced their ancestral peridinin-plastid by plastids derived from a variety of algal lineages ranging from green algae to diatoms. Here, we have characterized the genome of a dinoflagellate plastid of tertiary origin in order to understand the evolutionary processes that have shaped the organelle since it was acquired as a symbiont cell. To address this, the genome of the haptophyte-derived plastid in Karlodinium veneficum was analyzed by Sanger sequencing of library clones and 454 pyrosequencing of plastid enriched DNA fractions. The sequences were assembled into a single contig of 143 kb, encoding 70 proteins, 3 rRNAs and a nearly full set of tRNAs. Comparative genomics revealed massive rearrangements and gene losses compared to the haptophyte plastid; only a small fraction of the gene clusters usually found in haptophytes as well as other types of plastids are present in K. veneficum. Despite the reduced number of genes, the K. veneficum plastid genome has retained a large size due to expanded intergenic regions. Some of the plastid genes are highly diverged and may be pseudogenes or subject to RNA editing. Gene losses and rearrangements are also features of the genomes of the peridinin-containing plastids, apicomplexa and Chromera, suggesting that the evolutionary processes that once shaped these plastids have occurred at multiple independent occasions over the history of the Alveolata.     

Tobacco plastid ribosomal protein S18 is essential for cell survival

Plastid genomes contain a conserved set of genes most of which are involved in either photosynthesis or gene expression. Among the ribosomal protein genes present in higher plant plastid genomes, rps18 is special in that it is absent from the plastid genomes of several non-green unicellular organisms, including Euglena longa and Toxoplasma gondii. Here we have tested whether the ribosomal protein S18 is required for translation by deleting the rps18 gene from the tobacco plastid genome. We report that, while deletion of the rps18 gene was readily obtained, no homoplasmic Δrps18 plants or leaf sectors could be isolated. Instead, segregation into homoplasmy led to severe defects in leaf development suggesting that the knockout of rps18 is lethal and the S18 protein is required for cell survival. Our data demonstrate that S18 is indispensable for plastid ribosome function in tobacco and support an essential role for plastid translation in plant development. Moreover, we demonstrate the occurrence of flip-flop recombination on short inverted repeat sequences which generates different isoforms of the transformed plastid genome that differ in the orientation a 70 kb segment in the large single-copy region. However, infrequent occurrence of flip-flop recombination and random segregation of plastid genomes result in the predominant presence of only one of the isoforms in many tissue samples. Implications for the interpretation of chloroplast transformation experiments and vector design are discussed.     

PHYLOGENY OF PHAGOTROPHIC EUGLENIDS (EUGLENOZOA): A MOLECULAR APPROACH BASED ON CULTURE MATERIAL AND ENVIRONMENTAL SAMPLES1

Molecular studies based on small subunit (SSU) rDNA sequences addressing euglenid phylogeny hitherto suffered from the lack of available data about phagotrophic species. To extend the taxon sampling, SSU rRNA genes from species of seven genera of phagotrophic euglenids were investigated. Sequence analyses revealed an increasing genetic diversity among euglenid SSU rDNA sequences compared with other well‐known eukaryotic groups, reflecting an equally broad diversity of morphological characters among euglenid phagotrophs. Phylogenetic inference using standard parsimony and likelihood approaches as well as Bayesian inference and spectral analyses revealed no clear support for euglenid monophyly. Among phagotrophs, monophyly of Petalomonas cantuscygni and Notosolenus ostium, both comprising simple ingestion apparatuses, is strongly supported. A moderately supported clade comprises phototrophic euglenids and primary osmotrophic euglenids together with phagotrophs, exhibiting a primarily flexible pellicle composed of numerous helically arranged strips and a complex ingestion apparatus with two supporting rods and four curved vanes. Comparison of molecular and morphological data is used to demonstrate the difficulties to formulate a hypothesis about how the ingestion apparatus evolved in this group.     

In situ plastid and mitochondrial DNA determination: Implication of the observed MINIMAL PLASTID GENOME NUMBER

The total loss of plastid DNA has never been reported for any alga or plant cell line, with the sole exception of the protozoan Euglena, yet plastid distribution at mitosis is apparently stochastric (Birky and Skavaril, Journal of Theoretical Biology, vol. 106, pp. 441–447, 1984) and accidental loss might be expected. It is not obvious how stem cells of photosynthetic eukaryotes avoid this problem. The chrysophyte alga Ochromonas danica, described as having but one or two plastids, can proliferate indefinitely without the benefit of photosynthesis. Under such conditions its plastid genome copy number per cell might drop to the absolute minimum compatible with maintaining its inheritance. In situ quantitation of Ochromonas plastid DNA in both photosynthetic and enriched mixotrophic growth, and in heterotrophic growth in prolonged darkness, suggests that plastids are capable of very wide variation (7 to >;200 genomes/plastid) in their DNA content, and likewise, cells can vary from one to >;8 plastids per cell, with total genomes numbers from 7 to >;1,000 per cell. Among many growth conditions tested, the smallest plastids were found in rapidly dividing cells grown in the dark, many of which contained but one plastid. The inability to find plastids with fewer than seven plastid genome equivalents of DNA, even in these rapidly multiplying cells grown in total darkness for months, suggests that multiple copies of the plastid genome may be very carefully maintained, even in the prolonged absence of photosynthesis. This implies that multiple copies are important for reasons other than photosynthetic capability; two possibilities are the biosynthetic steps necessary for eukaryote cell survival known to occur solely within a plastid, and/or the potential that multiple plastid genome copies provide to escape the effects of Muller's ratchet.     

Hypothesis: Gene‐rich plastid genomes in red algae may be an outcome of nuclear genome reduction

Red algae (Rhodophyta) putatively diverged from the eukaryote tree of life >1.2 billion years ago and are the source of plastids in the ecologically important diatoms, haptophytes, and dinoflagellates. In general, red algae contain the largest plastid gene inventory among all such organelles derived from primary, secondary, or additional rounds of endosymbiosis. In contrast, their nuclear gene inventory is reduced when compared to their putative sister lineage, the Viridiplantae, and other photosynthetic lineages. The latter is thought to have resulted from a phase of genome reduction that occurred in the stem lineage of Rhodophyta. A recent comparative analysis of a taxonomically broad collection of red algal and Viridiplantae plastid genomes demonstrates that the red algal ancestor encoded ~1.5× more plastid genes than Viridiplantae. This difference is primarily explained by more extensive endosymbiotic gene transfer (EGT) in the stem lineage of Viridiplantae, when compared to red algae. We postulate that limited EGT in Rhodophytes resulted from the countervailing force of ancient, and likely recurrent, nuclear genome reduction. In other words, the propensity for nuclear gene loss led to the retention of red algal plastid genes that would otherwise have undergone intracellular gene transfer to the nucleus. This hypothesis recognizes the primacy of nuclear genome evolution over that of plastids, which have no inherent control of their gene inventory and can change dramatically (e.g., secondarily non‐photosynthetic eukaryotes, dinoflagellates) in response to selection acting on the host lineage.