Invited Review Beyond parasitic convergence: unravelling the evolution of the organellar genomes in holoparasites

BACKGROUND

The molecular evolution of organellar genomes in angiosperms has been studied extensively, with some lineages, such as parasitic ones, displaying unique characteristics. Parasitism has emerged 12 times independently in angiosperm evolution. Holoparasitism is the most severe form of parasitism, and is found in ~10 % of parasitic angiosperms. Although a few holoparasitic species have been examined at the molecular level, most reports involve plastomes instead of mitogenomes. Parasitic plants establish vascular connections with their hosts through haustoria to obtain water and nutrients, which facilitates the exchange of genetic information, making them more susceptible to horizontal gene transfer (HGT). HGT is more prevalent in the mitochondria than in the chloroplast or nuclear compartments.

SCOPE

This review summarizes current knowledge on the plastid and mitochondrial genomes of holoparasitic angiosperms, compares the genomic features across the different lineages, and discusses their convergent evolutionary trajectories and distinctive features. We focused on Balanophoraceae (Santalales), which exhibits extraordinary traits in both their organelles.

CONCLUSIONS

Apart from morphological similarities, plastid genomes of holoparasitic plants also display other convergent features, such as rampant gene loss, biased nucleotide composition and accelerated evolutionary rates. In addition, the plastomes of Balanophoraceae have extremely low GC and gene content, and two unexpected changes in the genetic code. Limited data on the mitochondrial genomes of holoparasitic plants preclude thorough comparisons. Nonetheless, no obvious genomic features distinguish them from the mitochondria of free-living angiosperms, except for a higher incidence of HGT. HGT appears to be predominant in holoparasitic angiosperms with a long-lasting endophytic stage. Among the Balanophoraceae, mitochondrial genomes exhibit disparate evolutionary paths with notable levels of heteroplasmy in Rhopalocnemis and unprecedented levels of HGT in Lophophytum. Despite their differences, these Balanophoraceae share a multichromosomal mitogenome, a feature also found in a few free-living angiosperms.

Evolution of cox2 introns in angiosperm mitochondria and efficient splicing of an elongated cox2i691 intron

In angiosperms, the mitochondrial cox2 gene harbors up to two introns, commonly referred to as cox2i373 and cox2i691. We studied the cox2 from 222 fully-sequenced mitogenomes from 30 angiosperm orders and analyzed the evolution of their introns. Unlike cox2i373, cox2i691 shows a distribution among plants that is shaped by frequent intron loss events driven by localized retroprocessing. In addition, cox2i691 exhibits sporadic elongations, frequently in domain IV of introns. Such elongations are poorly related to repeat content and two of them showed the presence of LINE transposons, suggesting that increasing intron size is very likely due to nuclear intracelular DNA transfer followed by incorporation into the mitochondrial DNA. Surprisingly, we found that cox2i691 is erroneously annotated as absent in 30 mitogenomes deposited in public databases. Although each of the cox2 introns is ∼1.5 kb in length, a cox2i691 of 4.2 kb has been reported in Acacia ligulata (Fabaceae). It is still unclear whether its unusual length is due to a trans-splicing arrangement or the loss of functionality of the interrupted cox2. Through analyzing short-read RNA sequencing of Acacia with a multi-step computational strategy, we found that the Acacia cox2 is functional and its long intron is spliced in cis in a very efficient manner despite its length.

Break-induced replication is the primary recombination pathway in plant somatic hybrid mitochondria: a model for mt-HGT

Somatic hybrids between distant species offer a remarkable model to study genomic recombination events after mitochondrial fusion. Recently, our lab described highly chimeric mitogenomes in two somatic hybrids between the Solanaceae Nicotiana tabacum and Hyoscyamus niger resulting from interparental homologous recombination. To better examine the recombination map in somatic hybrid mitochondria, we developed a more sensitive bioinformatic strategy to detect recombination activity based on high-throughput sequencing without assembling the hybrid mitogenome. We generated a new intergeneric somatic hybrid between N. tabacum and Physochlaina orientalis and re-analyzed the somatic hybrids previously generated in our lab. We inferred 213 homologous recombination events across repeats of 2.1 kb on average. Most of them (~80%) were asymmetrical, consistent with the break-induced replication (BIR) pathway. Only rare (2.74%) non-homologous events were detected. Interestingly, independent events frequently occurred in the same regions within and across somatic hybrids, suggesting the existence of recombination hotspots in plant mitogenomes. BIR is the main pathway of interparental recombination in somatic hybrid mitochondria. Findings of this study are relevant to mitogenome editing assays and to mechanistic aspects of DNA integration following mitochondrial DNA horizontal transfer events.

Cytonuclear coevolution in a holoparasitic plant with highly disparate organellar genomes

Contrasting substitution rates in the organellar genomes of Lophophytum agree with the DNA repair, replication, and recombination gene content. Plastid and nuclear genes whose products form multisubunit complexes co-evolve. The organellar genomes of the holoparasitic plant Lophophytum (Balanophoraceae) show disparate evolution. In the plastid, the genome has been severely reduced and presents a > 85% AT content, while in the mitochondria most protein-coding genes have been replaced by homologs acquired by horizontal gene transfer (HGT) from their hosts (Fabaceae). Both genomes carry genes whose products form multisubunit complexes with those of nuclear genes, creating a possible hotspot of cytonuclear coevolution. In this study, we assessed the evolutionary rates of plastid, mitochondrial and nuclear genes, and their impact on cytonuclear evolution of genes involved in multisubunit complexes related to lipid biosynthesis and proteolysis in the plastid and those in charge of the oxidative phosphorylation in the mitochondria. Genes from the plastid and the mitochondria (both native and foreign) of Lophophytum showed extremely high and ordinary substitution rates, respectively. These results agree with the biased loss of plastid-targeted proteins involved in angiosperm organellar repair, replication, and recombination machinery. Consistent with the high rate of evolution of plastid genes, nuclear-encoded subunits of plastid complexes showed disproportionate increases in non-synonymous substitution rates, while those of the mitochondrial complexes did not show different rates than the control (i.e. non-organellar nuclear genes). Moreover, the increases in the nuclear-encoded subunits of plastid complexes were positively correlated with the level of physical interaction they possess with the plastid-encoded ones. Overall, these results suggest that a structurally-mediated compensatory factor may be driving plastid-nuclear coevolution in Lophophytum, and that mito-nuclear coevolution was not altered by HGT.

Native and foreign mitochondrial and nuclear encoded proteins conform the OXPHOS complexes of a holoparasitic plant

The intimate contact between the holoparasitic plant Lophophytum mirabile (Balanophoraceae) and its host plant (Fabaceae) facilitates the exchange of genetic information, increasing the frequency of horizontal gene transfer (HGT). Lophophytum stands out because it acquired a large number of mitochondrial genes (greater than 20) from its legume host that replaced the majority of the native homologs. These foreign genes code for proteins that form multisubunit enzyme complexes, such as those in the oxidative phosphorylation system (OXPHOS) and cytochrome c maturation (ccm) system, together with dozens of nuclear-encoded subunits. However, the existence and the origin of the nuclear subunits that form the major part of the OXPHOS and ccm system in Lophophytum remain unknown. It was proposed that nuclear-encoding genes whose products interact with foreign mitochondrial proteins are also foreign, minimizing the incompatibilities that could arise in the assembly and functioning of these multiprotein complexes. We identified a nearly complete set of OXPHOS and ccm system subunits evolving under selective constraints in the transcriptome of Lophophytum, indicating that OXPHOS is functional and resembles that of free-living angiosperms. Maximum Likelihood phylogenetic analyses revealed a single case of HGT in the nuclear genes, which results in mosaic OXPHOS and ccm system in Lophophytum. These observations raise new questions about the evolution and physiology of this parasitic plant. A putative case of cooperation between two foreign (one mitochondrial and one nuclear) genes is presented.

Deepred-Mt: Deep representation learning for predicting C-to-U RNA editing in plant mitochondria

In land plant mitochondria, C-to-U RNA editing converts cytidines into uridines at highly specific RNA positions called editing sites. This editing step is essential for the correct functioning of mitochondrial proteins. When using sequence homology information, edited positions can be computationally predicted with high precision. However, predictions based on the sequence contexts of such edited positions often result in lower precision, which is limiting further advances on novel genetic engineering techniques for RNA regulation. Here, a deep convolutional neural network called Deepred-Mt is proposed. It predicts C-to-U editing events based on the 40 nucleotides flanking a given cytidine. Unlike existing methods, Deepred-Mt was optimized by using editing extent information, novel strategies of data augmentation, and a large-scale training dataset, constructed with deep RNA sequencing data of 21 plant mitochondrial genomes. In comparison to predictive methods based on sequence homology, Deepred-Mt attains significantly better predictive performance, in terms of average precision as well as F1 score. In addition, our approach is able to recognize well-known sequence motifs linked to RNA editing, and shows that the local RNA structure surrounding editing sites may be a relevant factor regulating their editing. These results demonstrate that Deepred-Mt is an effective tool for predicting C-to-U RNA editing in plant mitochondria. Source code, datasets, and detailed use cases are freely available at https://github.com/aedera/deepredmt.

Plastomes in the holoparasitic family Balanophoraceae: Extremely high AT content, severe gene content reduction, and two independent genetic code changes

The transition to a heterotrophic lifestyle in angiosperms is characterized by convergent evolutionary changes. Plastid genome remodeling includes dramatic functional and physical reductions with the highest degrees observed in fully heterotrophic plants. Genes related to photosynthesis are generally absent or pseudogenized, while a few genes related to other metabolic processes that take place within the plastid are almost invariably maintained. The family Balanophoraceae consists of root holoparasites that present reduced plastid genomes with an extraordinarily elevated AT content and the single genetic code change ever documented in land plant plastomes (the stop codon TAG now codes for tryptophan). Here, we studied the plastomes of Lophophytum leandri and Ombrophytum subterraneum (Balanophoraceae) that showed the remarkable absence of the gene trnE, a highly biased nucleotide composition, and an independent genetic code change (the standard stop codon TGA codes for tryptophan). This is the second genetic code change identified in land plant plastomes. Analysis of the transcriptome of Lophophytum indicated that the entire C5 pathway typical of plants is conserved despite the lack of trnE in its plastome. A hypothetical model of plastome evolution in the Balanophoraceae is presented.

Horizontal gene transfers dominate the functional mitochondrial gene space of a holoparasitic plant

Although horizontal gene transfer (HGT) is common in angiosperm mitochondrial DNAs (mtDNAs), few cases of functional foreign genes have been identified. The one outstanding candidate for large-scale functional HGT is the holoparasite Lophophytum mirabile, whose mtDNA has lost most native genes but contains intact foreign homologs acquired from legume host plants. To investigate the extent to which this situation results from functional replacement of native by foreign genes, functional mitochondrial gene transfer to the nucleus, and/or loss of mitochondrial biochemical function in the context of extreme parasitism, we examined the Lophophytum mitochondrial and nuclear transcriptomes by deep paired-end RNA sequencing. Most foreign mitochondrial genes in Lophophytum are highly transcribed, accurately spliced, and efficiently RNA edited. By contrast, we found no evidence for functional gene transfer to the nucleus or loss of mitochondrial functions in Lophophytum. Many functional replacements occurred via the physical replacement of native genes by foreign genes. Some of these events probably occurred as the final act of HGT itself. Lophophytum mtDNA has experienced an unprecedented level of functional replacement of native genes by foreign copies. This raises important questions concerning population-genetic and molecular regimes that underlie such a high level of foreign gene takeover.

Computational Detection of Plant RNA Editing Events

Computers are able to systematically exploit RNA-seq data allowing us to efficiently detect RNA editing sites in a genome-wide scale. This chapter introduces a very flexible computational framework for detecting RNA editing sites in plant organelles. This framework comprises three major steps: RNA-seq data processing, RNA read alignment, and RNA editing site detection. Each step is discussed in sufficient detail to be implemented by the reader. As a study case, the framework will be used with publicly available sequencing data to detect C-to-U RNA editing sites in the coding sequences of the mitochondrial genome of Nicotiana tabacum.

Multichromosomal structure and foreign tracts in the Ombrophytum subterraneum (Balanophoraceae) mitochondrial genome

Horizontal gene transfer (HGT) is frequent in parasitic plant mitochondria as a result of vascular connections established in host-parasite relationships. Recent studies of the holoparasitic plant Lophophytum mirabile (Balanophoraceae) revealed the unprecedented acquisition of a large amount of mitochondrial sequences from its legume host. We focused on a close relative, the generalist holoparasite Ombrophytum subterraneum, to examine the incidence of HGT events in the mitochondrial genome (mtDNA). The mtDNA of O. subterraneum assembles into 54 circular chromosomes, only 34 of which contain the 51 full-length coding regions. Numerous foreign tracts (totaling almost 100 kb, ~ 14% of the mtDNA), including 12 intact genes, were acquired by HGT from the Asteraceae hosts. Nine chromosomes concentrate most of those regions and eight are almost entirely foreign. Native homologs of each foreign gene coexist in the mtDNA and are potentially functional. A large proportion of shorter regions were related to the Fabaceae (a total of ~ 110 kb, 15.4%), some of which were shared with L. mirabile. We also found evidence of foreign sequences donated by angiosperm lineages not reported as hosts (Apocynaceae, Euphorbiaceae, Lamiaceae, and Malvales). We propose an evolutionary hypothesis that involves ancient transfers from legume hosts in the common ancestor of Ombrophytum and Lophophytum followed by more recent transfer events in L. mirabile. Besides, the O. subterraneum mtDNA was also subjected to additional HGT events from diverse angiosperm lineages, including large and recent transfers from the Asteraceae, and also from Lamiaceae.

Elucidating genomic patterns and recombination events in plant cybrid mitochondria

KEY MESSAGE

Cybrid plant mitochondria undergo homologous recombination, mainly BIR, keep a single allele for each gene, and maintain exclusive sequences of each parent and a single copy of the homologous regions. The maintenance of a dynamic equilibrium between the mitochondrial and nuclear genomes requires continuous communication and a high level of compatibility between them, so that alterations in one genetic compartment need adjustments in the other. The co-evolution of nuclear and mitochondrial genomes has been poorly studied, even though the consequences and effects of this interaction are highly relevant for human health, as well as for crop improvement programs and for genetic engineering. The mitochondria of plants represent an excellent system to understand the mechanisms of genomic rearrangements, chimeric gene formation, incompatibility between nucleus and cytoplasm, and horizontal gene transfer. We carried out detailed analyses of the mtDNA of a repeated cybrid between the solanaceae Nicotiana tabacum and Hyoscyamus niger. The mtDNA of the cybrid was intermediate between the size of the parental mtDNAs and the sum of them. Noticeably, most of the homologous sequences inherited from both parents were lost. In contrast, the majority of the sequences exclusive of a single parent were maintained. The mitochondrial gene content included a majority of N. tabacum derived genes, but also chimeric, two-parent derived, and H. niger-derived genes in a tobacco nuclear background. Any of these alterations in the gene content could be the cause of CMS in the cybrid. The parental mtDNAs interacted through 28 homologous recombination events and a single case of illegitimate recombination. Three main homologous recombination mechanisms were recognized in the cybrid mitochondria. Break induced replication (BIR) pathway was the most frequent. We propose that BIR could be one of the mechanisms responsible for the loss of the majority of the repeated regions derived from H. niger.

Towards a comprehensive picture of C-to-U RNA editing sites in angiosperm mitochondria

Our understanding of the dynamic and evolution of RNA editing in angiosperms is in part limited by the few editing sites identified to date. This study identified 10,217 editing sites from 17 diverse angiosperms. Our analyses confirmed the universality of certain features of RNA editing, and offer new evidence behind the loss of editing sites in angiosperms. RNA editing is a post-transcriptional process that substitutes cytidines (C) for uridines (U) in organellar transcripts of angiosperms. These substitutions mostly take place in mitochondrial messenger RNAs at specific positions called editing sites. By means of publicly available RNA-seq data, this study identified 10,217 editing sites in mitochondrial protein-coding genes of 17 diverse angiosperms. Even though other types of mismatches were also identified, we did not find evidence of non-canonical editing processes. The results showed an uneven distribution of editing sites among species, genes, and codon positions. The analyses revealed that editing sites were conserved across angiosperms but there were some species-specific sites. Non-synonymous editing sites were particularly highly conserved (~ 80%) across the plant species and were efficiently edited (80% editing extent). In contrast, editing sites at third codon positions were poorly conserved (~ 30%) and only partially edited (~ 40% editing extent). We found that the loss of editing sites along angiosperm evolution is mainly occurring by replacing editing sites with thymidines, instead of a degradation of the editing recognition motif around editing sites. Consecutive and highly conserved editing sites had been replaced by thymidines as result of retroprocessing, by which edited transcripts are reverse transcribed to cDNA and then integrated into the genome by homologous recombination. This phenomenon was more pronounced in eudicots, and in the gene cox1. These results suggest that retroprocessing is a widespread driving force underlying the loss of editing sites in angiosperm mitochondria.

Metazoan Ribosome Inactivating Protein encoding genes acquired by Horizontal Gene Transfer

Ribosome inactivating proteins (RIPs) are RNA N-glycosidases that depurinate a specific adenine residue in the conserved sarcin/ricin loop of 28S rRNA. These enzymes are widely distributed among plants and their presence has also been confirmed in several bacterial species. Recently, we reported for the first time in silico evidence of RIP encoding genes in metazoans, in two closely related species of insects: Aedes aegypti and Culex quinquefasciatus. Here, we have experimentally confirmed the presence of these genes in mosquitoes and attempted to unveil their evolutionary history. A detailed study was conducted, including evaluation of taxonomic distribution, phylogenetic inferences and microsynteny analyses, indicating that mosquito RIP genes derived from a single Horizontal Gene Transfer (HGT) event, probably from a cyanobacterial donor species. Moreover, evolutionary analyses show that, after the HGT event, these genes evolved under purifying selection, strongly suggesting they play functional roles in these organisms.

Unparalleled replacement of native mitochondrial genes by foreign homologs in a holoparasitic plant

Horizontal gene transfer (HGT) among flowering plant mitochondria occurs frequently and, in most cases, leads to nonfunctional transgenes in the recipient genome. Parasitic plants are particularly prone to this phenomenon, but their mitochondrial genomes (mtDNA) have been largely unexplored. We undertook a large-scale mitochondrial genomic study of the holoparasitic plant Lophophytum mirabile (Balanophoraceae). Comprehensive phylogenetic analyses were performed to address the frequency, origin, and impact of HGT. The sequencing of the complete mtDNA of L. mirabile revealed the unprecedented acquisition of host-derived mitochondrial genes, representing 80% of the protein-coding gene content. All but two of these foreign genes replaced the native homologs and are probably functional in energy metabolism. The genome consists of 54 circular-mapping chromosomes, 25 of which carry no intact genes. The likely functional replacement of up to 26 genes in L. mirabile represents a stunning example of the potential effect of rampant HGT on plant mitochondria. The use of host-derived genes may have a positive effect on the host-parasite relationship, but could also be the result of nonadaptive forces.

Dissimilar evolutionary histories of two resistance gene families in the genus Solanum

Genomic analyses have shown that most genes in eukaryotic lineages belong to families. Gene families vary in terms of number of members, nucleotide similarity, gene integrity, expression, and function. Often, the members of gene families are arranged in clusters, which contribute to maintaining similarity among gene copies and also to generate duplicates through replication errors. Gene families offer us an opportunity to examine the forces involved in the evolution of the genomes and to study recombination events and genomic rearrangements. In this work, we focused on the evolution of two plant resistance gene families, Sw5 and Mi-1, and analyzed the completely sequenced nuclear genomes of potato and tomato. We first noticed that the potato genome carries larger resistance gene families than tomato, but all gene copies are pseudogenes. Second, phylogenetic analyses indicated that Sw5 and Mi-1 gene families had dissimilar evolutionary histories. In contrast to Sw5, Mi-1 homologues suffered repeated gene conversion events among the gene copies, particularly in the tomato genome.

Homologous recombination and retention of a single form of most genes shape the highly chimeric mitochondrial genome of a cybrid plant

The structure and evolution of angiosperm mitochondrial genomes are driven by extremely high rates of recombination and rearrangement. An excellent experimental system for studying these events is offered by cybrid plants, in which parental mitochondria usually fuse and their genomes recombine. Little is known about the extent, nature and consequences of mitochondrial recombination in these plants. We conducted the first study in which the organellar genomes of a cybrid - between Nicotiana tabacum and Hyoscyamus niger - were sequenced and compared to those of its parents. This cybrid mitochondrial genome is highly recombinant, reflecting at least 30 crossovers and five gene conversions between its parental genomes. It is also surprisingly large (41% and 64% larger than the parental genomes), yet contains single alleles for 90% of mitochondrial genes. Recombination produced a remarkably chimeric cybrid mitochondrial genome and occurred entirely via homologous mechanisms involving the double-strand break repair and/or break-induced replication pathways. Retention of a single form of most genes could be advantageous to minimize intracellular incompatibilities and/or reflect neutral forces that preferentially eliminate duplicated regions. We discuss the relevance of these findings to the surprisingly frequent occurrence of horizontal gene - and genome - transfer in angiosperm mitochondrial DNAs.

The chloroplast genome of Hyoscyamus niger and a phylogenetic study of the tribe Hyoscyameae (Solanaceae)

The tribe Hyoscyameae (Solanaceae) is restricted to Eurasia and includes the genera Archihyoscyamus, Anisodus, Atropa, Atropanthe, Hyoscyamus, Physochlaina, Przewalskia and Scopolia. Even though the monophyly of Hyoscyameae is strongly supported, the relationships of the taxa within the tribe remain unclear. Chloroplast markers have been widely used to elucidate plant relationships at low taxonomic levels. Identification of variable chloroplast intergenic regions has been developed based on comparative genomics of chloroplast genomes, but these regions have a narrow phylogenetic utility. In this study, we present the chloroplast genome sequence of Hyoscyamus niger and make comparisons to other solanaceous plastid genomes in terms of gene order, gene and intron content, editing sites, origins of replication, repeats, and hypothetical open reading frames. We developed and sequenced three variable plastid markers from eight species to elucidate relationships within the tribe Hyoscyameae. The presence of a horizontally transferred intron in the mitochondrial cox1 gene of some species of the tribe is considered here a likely synapomorphy uniting five genera of the Hyoscyameae. Alternatively, the cox1 intron could be a homoplasious character acquired twice within the tribe. A homoplasious inversion in the intergenic plastid spacer trnC-psbM was recognized as a source of bias and removed from the data set used in the phylogenetic analyses. Almost 12 kb of plastid sequence data were not sufficient to completely resolve relationships among genera of Hyoscyameae but some clades were identified. Two alternative hypotheses of the evolution of the genera within the tribe are proposed.

Horizontal transfer of entire genomes via mitochondrial fusion in the angiosperm Amborella

We report the complete mitochondrial genome sequence of the flowering plant Amborella trichopoda. This enormous, 3.9-megabase genome contains six genome equivalents of foreign mitochondrial DNA, acquired from green algae, mosses, and other angiosperms. Many of these horizontal transfers were large, including acquisition of entire mitochondrial genomes from three green algae and one moss. We propose a fusion-compatibility model to explain these findings, with Amborella capturing whole mitochondria from diverse eukaryotes, followed by mitochondrial fusion (limited mechanistically to green plant mitochondria) and then genome recombination. Amborella's epiphyte load, propensity to produce suckers from wounds, and low rate of mitochondrial DNA loss probably all contribute to the high level of foreign DNA in its mitochondrial genome.

A HYPOTHESIS FOR PLASTID EVOLUTION IN CHROMALVEOLATES(1)

Four eukaryotic lineages, namely, haptophytes, alveolates, cryptophytes, and heterokonts, contain in most cases photosynthetic and nonphotosynthetic members-the photosynthetic ones with secondary plastids with chl c as the main photosynthetic pigment. These four photosynthetic lineages were grouped together on the basis of their pigmentation and called chromalveolates, which is usually understood to imply loss of plastids in the nonphotosynthetic members. Despite the ecological and economic importance of this group of organisms, the phylogenetic relationships among these algae are only partially understood, and the so-called chromalveolate hypothesis is very controversial. This review evaluates the evidence for and against this grouping and summarizes the present understanding of chromalveolate evolution. We also describe a testable hypothesis that is intended to accommodate current knowledge based on plastid and nuclear genomic data, discuss the implications of this model, and comment on areas that require further examination.

Frequent, phylogenetically local horizontal transfer of the cox1 group I Intron in flowering plant mitochondria

Horizontal gene transfer is surprisingly common among plant mitochondrial genomes. The first well-established case involves a homing group I intron in the mitochondrial cox1 gene shown to have been frequently acquired via horizontal transfer in angiosperms. Here, we report extensive additional sampling of angiosperms, including 85 newly sequenced introns from 30 families. Analysis of all available data leads us to conclude that, among the 640 angiosperms (from 212 families) whose cox1 intron status has been characterized thus far, the intron has been acquired via roughly 70 separate horizontal transfer events. We propose that the intron was originally seeded into angiosperms by a single transfer from fungi, with all subsequent inferred transfers occurring from one angiosperm to another. The pattern of angiosperm-to-angiosperm transfer is biased toward exchanges between plants belonging to the same family. Illegitimate pollination is proposed as one potential factor responsible for this pattern, given that aberrant, cross-species pollination is more likely between close relatives. Other potential factors include shared vectoring agents or common geographic locations. We report the first apparent cases of loss of the cox1 intron; losses are accompanied by retention of the exonic coconversion tract, which is located immediately downstream of the intron and which is a product of the intron's self-insertion mechanism. We discuss the many reasons why the cox1 intron is so frequently and detectably transferred, and rarely lost, and conclude that it should be regarded as the "canary in the coal mine" with respect to horizontal transfer in angiosperm mitochondria.

Sorting wheat from chaff in multi-gene analyses of chlorophyll c-containing plastids

Photosynthetic eukaryotes contain primary, secondary or tertiary plastids, depending on the source of the organelle (a cyanobacterium or a photosynthetic eukaryote). Plastid phylogeny is relatively well investigated, but molecular phylogenies have conflicted as a function of gene choice, taxon-representations, and analytical method. To better understand the influences of these variables, we performed analyses of a multi-gene data set based on 62 plastid-associated genes of 15 taxa representing the major plastid lineages. In an attempt to distinguish phylogenetic signal from non-phylogenetic patterns, we analyzed the data using a wide range of phylogenetic methods and examined the effect of covarion evolution and compositional bias. The data suggest that the chlorophyll c-containing plastids are monophyletic and acquired their plastids from the red algae after the emergence of the Cyanidiales. The relationships among chl c-containing plastids are particularly hard to resolve. This is the largest data set used for this purpose; the analyses show that cryptophyte plastids are sister to other chl c-containing plastids, and haptophyte and peridinin-containing dinoflagellate plastids are closely related.

Plastid genes in a non-photosynthetic dinoflagellate

Dinoflagellates are a diverse group of protists, comprising photosynthetic and heterotrophic free-living species, as well as parasitic ones. About half of them are photosynthetic with peridinin-containing plastids being the most common. It is uncertain whether non-photosynthetic dinoflagellates are primitively so, or have lost photosynthesis. Studies of heterotrophic species from this lineage may increase our understanding of plastid evolution. We analyzed an EST project of the early-diverging heterotrophic dinoflagellate Crypthecodinium cohnii looking for evidence of past endosymbiosis. A large number of putative genes of cyanobacterial or algal origin were identified using BLAST, and later screened by metabolic function. Phylogenetic analyses suggest that several proteins could have been acquired from a photosynthetic endosymbiont, arguing for an earlier plastid acquisition in dinoflagellates. In addition, intact N-terminal plastid-targeting peptides were detected, indicating that C. cohnii may contain a reduced plastid and that some of these proteins are imported into this organelle. A number of metabolic pathways, such as heme and isoprenoid biosynthesis, seem to take place in the plastid. Overall, these data indicate that C. cohnii is derived from a photosynthetic ancestor and provide a model for loss of photosynthesis in dinoflagellates and their relatives. This represents the first extensive genomic analysis of a heterotrophic dinoflagellate.

Rate variation as a function of gene origin in plastid-derived genes of peridinin-containing dinoflagellates

Peridinin-pigmented dinoflagellates contain secondary plastids that seem to have undergone more nearly complete plastid genome reduction than other eukaryotes. Many typically plastid-encoded genes appear to have been transferred to the nucleus, with a few remaining genes found on minicircles. To understand better the evolution of the dinoflagellate plastid, four categories of plastid-associated genes in dinoflagellates were defined based on their history of transfer and evaluated for rate of sequence evolution, including minicircle genes (presumably plastid-encoded), genes probably transferred from the plastid to the nucleus (plastid-transferred), and genes that were likely acquired directly from the nucleus of the previous plastid host (nuclear-transferred). The fourth category, lateral-transferred genes, are plastid-associated genes that do not appear to have a cyanobacterial origin. The evolutionary rates of these gene categories were compared using relative rate tests and likelihood ratio tests. For comparison with other secondary plastid-containing organisms, rates were calculated for the homologous sequences from the haptophyte Emiliania huxleyi. The evolutionary rate of minicircle and plastid-transferred genes in the dinoflagellate was strikingly higher than that of nuclear-transferred and lateral-transferred genes and, also, substantially higher than that of all plastid-associated genes in the haptophyte. Plastid-transferred genes in the dinoflagellate had an accelerated rate of evolution that was variable but, in most cases, not as extreme as the minicircle genes. Furthermore, the nuclear-transferred and lateral-transferred genes showed rates of evolution that are similar to those of other taxa. Thus, nucleus-to-nucleus transferred genes have a more typical rate of sequence evolution, while those whose history was wholly or partially within the dinoflagellate plastid genome have a markedly accelerated rate of evolution.