Why are plastid genomes retained in non-photosynthetic organisms?

Organellar phylogenomics at the epidendroid orchid base, with a focus on the mycoheterotrophic Wullschlaegelia

BACKGROUND AND AIMS

Heterotrophic plants have long been a challenge for systematists, exemplified by the base of the orchid subfamily Epidendroideae, which contains numerous mycoheterotrophic species.

METHODS

Here we address the utility of organellar genomes in resolving relationships at the epidendroid base, specifically employing models of heterotachy, or lineage-specific rate variation over time. We further conduct comparative analyses of plastid genome evolution in heterotrophs and structural variation in matK.

KEY RESULTS

We present the first complete plastid genomes (plastomes) of Wullschlaegelia, the sole genus of the tribe Wullschlaegelieae, revealing a highly reduced genome of 37 kb, which retains a fraction of the genes present in related autotrophs. Plastid phylogenomic analyses recovered a strongly supported clade composed exclusively of mycoheterotrophic species with long branches. We further analysed mitochondrial gene sets, which recovered similar relationships to those in other studies using nuclear data, but the placement of Wullschlaegelia remains uncertain. We conducted comparative plastome analyses among Wullschlaegelia and other heterotrophic orchids, revealing a suite of correlated substitutional and structural changes relative to autotrophic species. Lastly, we investigated evolutionary and structural variation in matK, which is retained in Wullschlaegelia and a few other 'late stage' heterotrophs and found evidence for structural conservation despite rapid substitution rates in both Wullschlaegelia and the leafless Gastrodia.

CONCLUSIONS

Our analyses reveal the limits of what the plastid genome can tell us on orchid relationships in this part of the tree, even when applying parameter-rich heterotachy models. Our study underscores the need for increased taxon sampling across all three genomes at the epidendroid base, and illustrates the need for further research on addressing heterotachy in phylogenomic analyses.

Photosynthetic demands on translational machinery drive retention of redundant tRNA metabolism in plant organelles

Eukaryotic nuclear genomes often encode distinct sets of translation machinery for function in the cytosol vs. organelles (mitochondria and plastids). This raises questions about why multiple translation systems are maintained even though they are capable of comparable functions and whether they evolve differently depending on the compartment where they operate. These questions are particularly interesting in plants because translation machinery, including aminoacyl-transfer RNA (tRNA) synthetases (aaRS), is often dual-targeted to the plastids and mitochondria. These organelles have different functions, with much higher rates of translation in plastids to supply the abundant, rapid-turnover proteins required for photosynthesis. Previous studies have indicated that plant organellar aaRS evolve more slowly compared to mitochondrial aaRS in eukaryotes that lack plastids. Thus, we investigated the evolution of nuclear-encoded organellar and cytosolic aaRS and tRNA maturation enzymes across a broad sampling of angiosperms, including nonphotosynthetic (heterotrophic) plant species with reduced plastid gene expression, to test the hypothesis that translational demands associated with photosynthesis constrain the evolution of enzymes involved in organellar tRNA metabolism. Remarkably, heterotrophic plants exhibited wholesale loss of many organelle-targeted aaRS and other enzymes, even though translation still occurs in their mitochondria and plastids. These losses were often accompanied by apparent retargeting of cytosolic enzymes and tRNAs to the organelles, sometimes preserving aaRS-tRNA charging relationships but other times creating surprising mismatches between cytosolic aaRS and mitochondrial tRNA substrates. Our findings indicate that the presence of a photosynthetic plastid drives the retention of specialized systems for organellar tRNA metabolism.

Evolution of the Genetic Code in the Ascoideales (CUG-Ser2) Yeast Clade: The Ancestral tRNA-Leu(CAG) Gene Is Retained in Most Saccharomycopsis Species but Is Nonessential and Not Used for Translation

In the yeast genera Saccharomycopsis and Ascoidea, which comprise the taxonomic order Ascoideales, nuclear genes use a nonstandard genetic code in which CUG codons are translated as serine instead of leucine, due to a tRNA-Ser with the unusual anticodon CAG. However, some species in this clade also retain an ancestral tRNA-Leu gene with the same anticodon. One of these species, Ascoidea asiatica, has been shown to have a stochastic proteome in which proteins contain ∼50% Ser and 50% Leu at CUG codon sites, whereas previously examined Saccharomycopsis species translate CUG only as Ser. Here, we investigated the presence, conservation, and possible functionality of the tRNA-Leu(CAG) gene in the genus Saccharomycopsis. We sequenced the genomes of 23 strains that, together with previously available data, include almost every known species of this genus. We found that most Saccharomycopsis species have genes for both tRNA-Leu(CAG) and tRNA-Ser(CAG). However, tRNA-Leu(CAG) has been lost in Saccharomycopsis synnaedendra and Saccharomycopsis microspora, and its predicted cloverleaf structure is aberrant in all the other Saccharomycopsis species. We deleted the tRNA-Leu(CAG) gene of Saccharomycopsis capsularis and found that it is not essential. Proteomic analyses in vegetative and sporulating cultures of S. capsularis and Saccharomycopsis fermentans showed only translation of CUG as Ser. Despite its unusual structure, the tRNA-Leu(CAG) gene shows evidence of sequence conservation among Saccharomycopsis species, particularly in its acceptor stem and leucine identity elements, which suggests that it may have been retained in order to carry out an unknown nontranslational function.

Dinoflagellate chloroplasts as a model for extreme genome reduction and fragmentation in organelles - The COCOA principle for gene retention

The genomes of peridinin-containing dinoflagellate chloroplasts have a very unusual organisation. These genomes are highly fragmented and greatly reduced, with most of the usual complement of chloroplast genes relocated to the nucleus. Dinoflagellate chloroplasts highlight evolutionary changes that are found to varying extents in a number of other organelle genomes. These include the chloroplast genome of the green alga Boodlea and other Cladophorales, and the mitochondrial genomes of blood-sucking and chewing lice, the parasitic plant Rhopalocnemis phalloides, the red alga Rhodosorus marinus and other members of the Stylonematophyceae, diplonemid flagellates, and some Cnidaria. Consideration of the coding content of the remnant chloroplast genomes indicates that organelles may preferentially retain genes for proteins important in initiating assembly of complexes, and the same is largely true for mitochondria. We propose a new principle, of CO-location for COntrol of Assembly (COCOA), indicating the importance of retaining these genes in the organelle. This adds to, but does not invalidate, the existing hypotheses of the multisubunit completion principle, CO-location for Redox Regulation (CORR) and Control by Epistasy of Synthesis (CES).

Ecological Predictors of Organelle Genome Evolution: Phylogenetic Correlations with Taxonomically Broad, Sparse, Unsystematized Data

Comparative analysis of variables across phylogenetically linked observations can reveal mechanisms and insights in evolutionary biology. As the taxonomic breadth of the sample of interest increases, challenges of data sparsity, poor phylogenetic resolution, and complicated evolutionary dynamics emerge. Here, we investigate a cross-eukaryotic question where all these problems exist: which organismal ecology features are correlated with gene retention in mitochondrial and chloroplast DNA (organelle DNA or oDNA). Through a wide palette of synthetic control studies, we first characterize the specificity and sensitivity of a collection of parametric and non-parametric phylogenetic comparative approaches to identify relationships in the face of such sparse and awkward datasets. This analysis is not directly focused on oDNA, and so provides generalizable insights into comparative approaches with challenging data. We then combine and curate ecological data coupled to oDNA genome information across eukaryotes, including a new semi-automated approach for gathering data on organismal traits from less systematized open-access resources including encyclopedia articles on species and taxa. The curation process also involved resolving several issues with existing datasets, including enforcing the clade-specificity of several ecological features and fixing incorrect annotations. Combining this unique dataset with our benchmarked comparative approaches, we confirm support for several known links between organismal ecology and organelle gene retention, identify several previously unidentified relationships constituting possible ecological contributors to oDNA genome evolution, and provide support for a recently hypothesized link between environmental demand and oDNA retention. We, with caution, discuss the implications of these findings for organelle evolution and of this pipeline for broad comparative analyses in other fields.

Connecting Species-Specific Extents of Genome Reduction in Mitochondria and Plastids

Mitochondria and plastids have both dramatically reduced their genomes since the endosymbiotic events that created them. The similarities and differences in the evolution of the two organelle genome types have been the target of discussion and investigation for decades. Ongoing work has suggested that similar mechanisms may modulate the reductive evolution of the two organelles in a given species, but quantitative data and statistical analyses exploring this picture remain limited outside of some specific cases like parasitism. Here, we use cross-eukaryote organelle genome data to explore evidence for coevolution of mitochondrial and plastid genome reduction. Controlling for differences between clades and pseudoreplication due to relatedness, we find that extents of mtDNA and ptDNA gene retention are related to each other across taxa, in a generally positive correlation that appears to differ quantitatively across eukaryotes, for example, between algal and nonalgal species. We find limited evidence for coevolution of specific mtDNA and ptDNA gene pairs, suggesting that the similarities between the two organelle types may be due mainly to independent responses to consistent evolutionary drivers.

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: The plasticity of plastids

Many chloroplast-bearing plants and algae lost their photosynthetic activity during evolution but retained their chloroplasts for other functions. A group of dinoflagellate algae apparently lost one half of their photosynthetic machinery but retained the other, providing a novel mechanism for light perception.

Extreme plastomes in holoparasitic Balanophoraceae are not the norm

BACKGROUND

Balanophoraceae plastomes are known for their highly condensed and re-arranged nature alongside the most extreme nucleotide compositional bias known to date, culminating in two independent reconfigurations of their genetic code. Currently, a large portion of the Balanophoraceae diversity remains unexplored, hindering, among others, evolutionary pattern recognition. Here, we explored newly sequenced plastomes of Sarcophyte sanguinea and Thonningia sanguinea. The reconstructed plastomes were analyzed using various methods of comparative genomics based on a representative taxon sampling.

RESULTS

Sarcophyte, recovered sister to the other sampled Balanophoraceae s. str., has plastomes up to 50% larger than those currently published. Its gene set contains five genes lost in any other species, including matK. Five cis-spliced introns are maintained. In contrast, the Thonningia plastome is similarly reduced to published Balanophoraceae and retains only a single cis-spliced intron. Its protein-coding genes show a more biased codon usage compared to Sarcophyte, with an accumulation of in-frame TAG stop codons. Structural plastome comparison revealed multiple, previously unknown, structural rearrangements within Balanophoraceae.

CONCLUSIONS

For the "minimal plastomes" of Thonningia, we propose a genetic code change identical to sister genus Balanophora. Sarcophyte however differs drastically from our current understanding on Balanophoraceae plastomes. With a less-extreme nucleotide composition, there is no evidence for an altered genetic code. Using comparative genomics, we identified a hotspot for plastome reconfiguration in Balanophoraceae. Based on previously published and newly identified structural reconfigurations, we propose an updated model of evolutionary plastome trajectories for Balanophoraceae, illustrating a much greater plastome diversity than previously known.

Plastid phylogenomics and molecular evolution of Thismiaceae (Dioscoreales)

PREMISE

Species in Thismiaceae can no longer photosynthesize and instead obtain carbon from soil fungi. Here we infer Thismiaceae phylogeny using plastid genome data and characterize the molecular evolution of this genome.

METHODS

We assembled five Thismiaceae plastid genomes from genome skimming data, adding to previously published data for phylogenomic inference. We investigated plastid-genome structural changes, considering locally colinear blocks (LCBs). We also characterized possible shifts in selection pressure in retained genes by considering changes in the ratio of nonsynonymous to synonymous changes (ω).

RESULTS

Thismiaceae experienced two major pulses of gene loss around the early diversification of the family, with subsequent scattered gene losses across descendent lineages. In addition to massive size reduction, Thismiaceae plastid genomes experienced occasional inversions, and there were likely two independent losses of the plastid inverted repeat (IR) region. Retained plastid genes remain under generally strong purifying selection (ω << 1), with significant and sporadic weakening or strengthening in several instances. The bifunctional trnE-UUC gene of Thismia huangii may retain a secondary role in heme biosynthesis, despite a probable loss of functionality in protein translation. Several cis-spliced group IIA introns have been retained, despite the loss of the plastid intron maturase, matK.

CONCLUSIONS

We infer that most gene losses in Thismiaceae occurred early and rapidly, following the initial loss of photosynthesis in its stem lineage. As a species-rich, fully mycoheterotrophic lineage, Thismiaceae provide a model system for uncovering the unique and divergent ways in which plastid genomes evolve in heterotrophic plants.

Reconstruction of Plastid Proteomes of Apicomplexans and Close Relatives Reveals the Major Evolutionary Outcomes of Cryptic Plastids

Apicomplexans and related lineages comprise many obligate symbionts of animals; some of which cause notorious diseases such as malaria. They evolved from photosynthetic ancestors and transitioned into a symbiotic lifestyle several times, giving rise to species with diverse non-photosynthetic plastids. Here, we sought to reconstruct the evolution of the cryptic plastids in the apicomplexans, chrompodellids, and squirmids (ACS clade) by generating five new single-cell transcriptomes from understudied gregarine lineages, constructing a robust phylogenomic tree incorporating all ACS clade sequencing datasets available, and using these to examine in detail, the evolutionary distribution of all 162 proteins recently shown to be in the apicoplast by spatial proteomics in Toxoplasma. This expanded homology-based reconstruction of plastid proteins found in the ACS clade confirms earlier work showing convergence in the overall metabolic pathways retained once photosynthesis is lost, but also reveals differences in the degrees of plastid reduction in specific lineages. We show that the loss of the plastid genome is common and unexpectedly find many lineage- and species-specific plastid proteins, suggesting the presence of evolutionary innovations and neofunctionalizations that may confer new functional and metabolic capabilities that are yet to be discovered in these enigmatic organelles.

Phylogenomics and plastome evolution of a Brazilian mycoheterotrophic orchid, Pogoniopsis schenckii

PREMISE

Pogoniopsis likely represents an independent photosynthesis loss in orchids. We use phylogenomic data to better identify the phylogenetic placement of this fully mycoheterotrophic taxon, and investigate its molecular evolution.

METHODS

We performed likelihood analysis of plastid and mitochondrial phylogenomic data to localize the position of Pogoniopsis schenckii in orchid phylogeny, and investigated the evolution of its plastid genome.

RESULTS

All analyses place Pogoniopsis in subfamily Epidendroideae, with strongest support from mitochondrial data, which also place it near tribe Sobralieae with moderately strong support. Extreme rate elevation in Pogoniopsis plastid genes broadly depresses branch support; in contrast, mitochondrial genes are only mildly rate elevated and display very modest and localized reductions in bootstrap support. Despite considerable genome reduction, including loss of photosynthesis genes and multiple translation apparatus genes, gene order in Pogoniopsis plastomes is identical to related autotrophs, apart from moderately shifted inverted repeat (IR) boundaries. All cis-spliced introns have been lost in retained genes. Two plastid genes (accD, rpl2) show significant strengthening of purifying selection. A retained plastid tRNA gene (trnE-UUC) of Pogoniopsis lacks an anticodon; we predict that it no longer functions in translation but retains a secondary role in heme biosynthesis.

CONCLUSIONS

Slowly evolving mitochondrial genes clarify the placement of Pogoniopsis in orchid phylogeny, a strong contrast with analysis of rate-elevated plastome data. We documented the effects of the novel loss of photosynthesis: for example, despite massive gene loss, its plastome is fully colinear with other orchids, and it displays only moderate shifts in selective pressure in retained genes.

Endosymbiotic selective pressure at the origin of eukaryotic cell biology

The dichotomy that separates prokaryotic from eukaryotic cells runs deep. The transition from pro- to eukaryote evolution is poorly understood due to a lack of reliable intermediate forms and definitions regarding the nature of the first host that could no longer be considered a prokaryote, the first eukaryotic common ancestor, FECA. The last eukaryotic common ancestor, LECA, was a complex cell that united all traits characterising eukaryotic biology including a mitochondrion. The role of the endosymbiotic organelle in this radical transition towards complex life forms is, however, sometimes questioned. In particular the discovery of the asgard archaea has stimulated discussions regarding the pre-endosymbiotic complexity of FECA. Here we review differences and similarities among models that view eukaryotic traits as isolated coincidental events in asgard archaeal evolution or, on the contrary, as a result of and in response to endosymbiosis. Inspecting eukaryotic traits from the perspective of the endosymbiont uncovers that eukaryotic cell biology can be explained as having evolved as a solution to housing a semi-autonomous organelle and why the addition of another endosymbiont, the plastid, added no extra compartments. Mitochondria provided the selective pressures for the origin (and continued maintenance) of eukaryotic cell complexity. Moreover, they also provided the energetic benefit throughout eukaryogenesis for evolving thousands of gene families unique to eukaryotes. Hence, a synthesis of the current data lets us conclude that traits such as the Golgi apparatus, the nucleus, autophagosomes, and meiosis and sex evolved as a response to the selective pressures an endosymbiont imposes.

Danxiaorchismangdangshanensis (Orchidaceae, Epidendroideae), a new species from central Fujian Province based on morphological and genomic data

Danxiaorchismangdangshanensis, a new mycoheterotrophic species from Fujian Province, China, is described and illustrated. The new species is morphologically similar to D.singchiana, but its callus of labellum is a less distinctive Y-shape with three auricles on the apex, four pollinia that are narrowly elliptic in shape and equal in size, and it lacks fine roots. The plastome of D.mangdangshanensis is highly degraded. Phylogenetic analyses distinguished D.mangdangshanensis from its congeners, D.singchiana and D.yangii, with strong support based on nrITS + matK and plastomes, respectively.

The convoluted history of haem biosynthesis

The capacity of haem to transfer electrons, bind diatomic gases, and catalyse various biochemical reactions makes it one of the essential biomolecules on Earth and one that was likely used by the earliest forms of cellular life. Since the description of haem biosynthesis, our understanding of this multi-step pathway has been almost exclusively derived from a handful of model organisms from narrow taxonomic contexts. Recent advances in genome sequencing and functional studies of diverse and previously neglected groups have led to discoveries of alternative routes of haem biosynthesis that deviate from the 'classical' pathway. In this review, we take an evolutionarily broad approach to illuminate the remarkable diversity and adaptability of haem synthesis, from prokaryotes to eukaryotes, showing the range of strategies that organisms employ to obtain and utilise haem. In particular, the complex evolutionary histories of eukaryotes that involve multiple endosymbioses and horizontal gene transfers are reflected in the mosaic origin of numerous metabolic pathways with haem biosynthesis being a striking case. We show how different evolutionary trajectories and distinct life strategies resulted in pronounced tensions and differences in the spatial organisation of the haem biosynthesis pathway, in some cases leading to a complete loss of a haem-synthesis capacity and, rarely, even loss of a requirement for haem altogether.

De Novo Transcriptome Meta-Assembly of the Mixotrophic Freshwater Microalga Euglena gracilis

Euglena gracilis is a well-known photosynthetic microeukaryote considered as the product of a secondary endosymbiosis between a green alga and a phagotrophic unicellular belonging to the same eukaryotic phylum as the parasitic trypanosomatids. As its nuclear genome has proven difficult to sequence, reliable transcriptomes are important for functional studies. In this work, we assembled a new consensus transcriptome by combining sequencing reads from five independent studies. Based on a detailed comparison with two previously released transcriptomes, our consensus transcriptome appears to be the most complete so far. Remapping the reads on it allowed us to compare the expression of the transcripts across multiple culture conditions at once and to infer a functionally annotated network of co-expressed genes. Although the emergence of meaningful gene clusters indicates that some biological signal lies in gene expression levels, our analyses confirm that gene regulation in euglenozoans is not primarily controlled at the transcriptional level. Regarding the origin of E. gracilis, we observe a heavily mixed gene ancestry, as previously reported, and rule out sequence contamination as a possible explanation for these observations. Instead, they indicate that this complex alga has evolved through a convoluted process involving much more than two partners.

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.

Comparative Analysis of Plastid Genomes in the Non-photosynthetic Genus Thismia Reveals Ongoing Gene Set Reduction

Heterotrophic plants provide intriguing examples of reductive evolution. This is especially evident in the reduction of their plastid genomes, which can potentially proceed toward complete genome loss. Several milestones at the beginning of this path of degradation have been described; however, little is known about the latest stages of plastome reduction. Here we analyze a diversity of plastid genomes in a set of closely related non-photosynthetic plants. We demonstrate how a gradual loss of genes shapes the miniaturized plastomes of these plants. The subject of our study, the genus Thismia, represents the mycoheterotrophic monocot family Thismiaceae, a group that may have experienced a very ancient (60-80 mya) transition to heterotrophy. In all 18 species examined, the plastome is reduced to 14-18 kb and is highly AT-biased. The most complete observed gene set includes accD, seven ribosomal protein genes, three rRNA, and two tRNA genes. Different clades of Thismia have undergone further gene loss (complete absence or pseudogenization) compared to this set: in particular, we report two independent losses of rps2 and rps18.

Correlation Analysis Reveals an Important Role of GC Content in Accumulation of Deletion Mutations in the Coding Region of Angiosperm Plastomes

Variation in GC content is assumed to correlate with various processes, including mutation biases, recombination, and environmental parameters. To date, most genomic studies exploring the evolution of GC content have focused on nuclear genomes, but relatively few have concentrated on organelle genomes. We explored the mechanisms maintaining the GC content in angiosperm plastomes, with a particular focus on the hypothesis of phylogenetic dependence and the correlation with deletion mutations. We measured three genetic traits, namely, GC content, A/T tracts, and G/C tracts, in the coding region of plastid genomes for 1382 angiosperm species representing 350 families and 64 orders, and tested the phylogenetic signal. Then, we performed correlation analyses and revealed the variation in evolutionary rate of selected traits using RRphylo. The plastid GC content in the coding region varied from 28.10% to 43.20% across angiosperms, with a few non-photosynthetic species showing highly reduced values, highlighting the significance of functional constraints. We found strong phylogenetic signal in A/T tracts, but weak ones in GC content and G/C tracts, indicating adaptive potential. GC content was positively and negatively correlated with G/C and A/T tracts, respectively, suggesting a trade-off between these two deletion events. GC content evolved at various rates across the phylogeny, with significant increases in monocots and Lamiids, and a decrease in Fabids, implying the effects of some other factors. We hypothesize that variation in plastid GC content might be a mixed strategy of species to optimize fitness in fluctuating climates, partly through influencing the trade-off between AT → GC and GC → AT mutations.

Plastome of the mycoheterotrophic eudicot Exacum paucisquama (Gentianaceae) exhibits extensive gene loss and a highly expanded inverted repeat region

Mycoheterotrophic plants are highly specialized species able to acquire organic carbon from symbiotic fungi, with relaxed dependence on photosynthesis for carbon fixation. The relaxation of the functional constraint of photosynthesis and thereby the relaxed selective pressure on functional photosynthetic genes usually lead to substantial gene loss and a highly degraded plastid genome in heterotrophs. In this study, we sequenced and analyzed the plastome of the eudicot Exacum paucisquama, providing the first plastid genome of a mycoheterotroph in the family Gentianaceae to date. The E. paucisquama plastome was 44,028 bp in length, which is much smaller than the plastomes of autotrophic eudicots. Although the E. paucisquama plastome had a quadripartite structure, a distinct boundary shift was observed in comparison with the plastomes of other eudicots. We detected extensive gene loss and only 21 putative functional genes (15 protein-coding genes, four rRNA genes and two tRNA genes). Our results provide valuable information for comparative evolutionary analyses of plastomes of heterotrophic species belonging to different phylogenetic groups.

Severe Plastid Genome Size Reduction in a Mycoheterotrophic Orchid, Danxiaorchis singchiana, Reveals Heavy Gene Loss and Gene Relocations

Danxiaorchis singchiana (Orchidaceae) is a leafless mycoheterotrophic orchid in the subfamily Epidendroideae. We sequenced the complete plastome of D. singchiana. The plastome has a reduced size of 87,931 bp, which includes a pair of inverted repeat (IR) regions of 13,762 bp each that are separated by a large single copy (LSC) region of 42,575 bp and a small single copy (SSC) region of 17,831 bp. When compared to its sister taxa, Cremastra appendiculata and Corallorhiza striata var. involuta, D. singchiana showed an inverted gene block in the LSC and SSC regions. A total of 61 genes were predicted, including 21 tRNA, 4 rRNA, and 36 protein-coding genes. While most of the housekeeping genes were still intact and seem to be protein-coding, only four photosynthesis-related genes appeared presumably intact. The majority of the presumably intact protein-coding genes seem to have undergone purifying selection (dN/dS < 1), and only the psaC gene was positively selected (dN/dS > 1) when compared to that in Cr. appendiculata. Phylogenetic analysis of 26 complete plastome sequences from 24 species of the tribe Epidendreae had revealed that D. singchiana diverged after Cr. appendiculata and is sister to the genus Corallorhiza with strong bootstrap support (100%).

Comprehensive genomic analyses with 115 plastomes from algae to seed plants: structure, gene contents, GC contents, and introns

BACKGROUND

Chloroplasts are a common character in plants. The chloroplasts in each plant lineage have shaped their own genomes, plastomes, by structural changes and transferring many genes to nuclear genomes during plant evolution. Some plastid genes have introns that are mostly group II introns.

OBJECTIVE

This study aimed to get genomic and evolutionary insights on the plastomes from green algae to flowering plants.

METHODS

Plastomes of 115 species from green algae, bryophytes, pteridophytes (spore bearing vascular plants), gymnosperms, and angiosperms were mined from NCBI organelle genome database. Plastome structure, gene contents and GC contents were analyzed by the in-house developed Phyton code. Intronic features including presence/absence, length, intron phases were analyzed by manually in the annotated information in NCBI.

RESULTS

The canonical quadripartite structures were retained in most plastomes except of a few plastomes that had lost an invert repeat (IR). Expansion or reduction or deletion of IRs resulted in the length variation of the plastomes. The number of protein coding genes ranged from 40 to 92 with an average 79.43 ± 5.84 per plastome and gene losses were apparent in specific lineages. The number of trn genes ranged from 13 to 33 with an average 21.19 ± 2.42 per plastome. Ribosomal RNA genes, rrn, were located in the IRs so that they were present in a duplicate except of the species that had lost one of the IR. GC contents were variable from 24.9 to 51.0% with an average 38.21 ± 3.27%, indicating bias to high AT contents. Plastid introns were present in 18 protein coding genes, six trn genes, and one rrn gene. Intron losses occurred among the orthologous genes in different plant lineages. The plastid introns were long compared with the nuclear introns, which might be related with the spliceosome nuclear introns and self-splicing group II plastid introns. The trnK-UUU intron contained the maturase encoding matK gene except in the chlorophyte algae and monilophyte ferns in which the trnK-UUU was lost, but matK retained. There were many annotation artefacts in the intron positions in the NCBI database. In the analysis of intron phases, phase 0 introns were more frequent than those of phase 2 and 3 introns. Phase polymorphism was observed in the introns of clpP which was derived from nucleotide insertion. Plastid trn introns were long compared to the archaeal or eukaryotic nuclear tRNA introns. Of the six plastid trn introns, one was at the D loop and other five were at the anticodon loop. The insertion sites were conserved among the trn genes in archaea, eukaryotic nuclear and plastid tRNA genes.

CONCLUSIONS

Current study refurbrished the previous findings of structural variations, gene contents, and GC contents of the chloroplast genomes from green algae to flowering plants. The study also included some noble findings and discussions on the plastome introns including their length variations and phase variation. We also presented and corrected some false annotations on the introns in protein coding and tRNA genes in the genome database, which might be confirmed by the chloroplast transcriptome analysis in the future.

Dinoflagellates with relic endosymbiont nuclei as models for elucidating organellogenesis

Nucleomorphs are relic endosymbiont nuclei so far found only in two algal groups, cryptophytes and chlorarachniophytes, which have been studied to model the evolutionary process of integrating an endosymbiont alga into a host-governed plastid (organellogenesis). However, past studies suggest that DNA transfer from the endosymbiont to host nuclei had already ceased in both cryptophytes and chlorarachniophytes, implying that the organellogenesis at the genetic level has been completed in the two systems. Moreover, we have yet to pinpoint the closest free-living relative of the endosymbiotic alga engulfed by the ancestral chlorarachniophyte or cryptophyte, making it difficult to infer how organellogenesis altered the endosymbiont genome. To counter the above issues, we need novel nucleomorph-bearing algae, in which endosymbiont-to-host DNA transfer is on-going and for which endosymbiont/plastid origins can be inferred at a fine taxonomic scale. Here, we report two previously undescribed dinoflagellates, strains MGD and TGD, with green algal endosymbionts enclosing plastids as well as relic nuclei (nucleomorphs). We provide evidence for the presence of DNA in the two nucleomorphs and the transfer of endosymbiont genes to the host (dinoflagellate) genomes. Furthermore, DNA transfer between the host and endosymbiont nuclei was found to be in progress in both the MGD and TGD systems. Phylogenetic analyses successfully resolved the origins of the endosymbionts at the genus level. With the combined evidence, we conclude that the host-endosymbiont integration in MGD/TGD is less advanced than that in cryptophytes/chrorarachniophytes, and propose the two dinoflagellates as models for elucidating organellogenesis.

Highly Reduced Plastid Genomes of the Non-photosynthetic Dictyochophyceans Pteridomonas spp. (Ochrophyta, SAR) Are Retained for tRNA-Glu-Based Organellar Heme Biosynthesis

Organisms that have lost their photosynthetic capabilities are present in a variety of eukaryotic lineages, such as plants and disparate algal groups. Most of such non-photosynthetic eukaryotes still carry plastids, as these organelles retain essential biological functions. Most non-photosynthetic plastids possess genomes with varied protein-coding contents. Such remnant plastids are known to be present in the non-photosynthetic, bacteriovorous alga Pteridomonas danica (Dictyochophyceae, Ochrophyta), which, regardless of its obligatory heterotrophic lifestyle, has been reported to retain the typically plastid-encoded gene for ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) large subunit (rbcL). The presence of rbcL without photosynthetic activity suggests that investigating the function of plastids in Pteridomonas spp. would likely bring unique insights into understanding the reductive evolution of plastids, their genomes, and plastid functions retained after the loss of photosynthesis. In this study, we demonstrate that two newly established strains of the non-photosynthetic genus Pteridomonas possess highly reduced plastid genomes lacking rbcL gene, in contrast to the previous report. Interestingly, we discovered that all plastid-encoded proteins in Pteridomonas spp. are involved only in housekeeping processes (e.g., transcription, translation and protein degradation), indicating that all metabolite synthesis pathways in their plastids are supported fully by nuclear genome-encoded proteins. Moreover, through an in-depth survey of the available transcriptomic data of another strain of the genus, we detected no candidate sequences for nuclear-encoded, plastid-directed Fe-S cluster assembly pathway proteins, suggesting complete loss of this pathway in the organelle, despite its widespread conservation in non-photosynthetic plastids. Instead, the transcriptome contains plastid-targeted components of heme biosynthesis, glycolysis, and pentose phosphate pathways. The retention of the plastid genomes in Pteridomonas spp. is not explained by the Suf-mediated constraint against loss of plastid genomes, previously proposed for Alveolates, as they lack Suf genes. Bearing all these findings in mind, we propose the hypothesis that plastid DNA is retained in Pteridomonas spp. for the purpose of providing glutamyl-tRNA, encoded by trnE gene, as a substrate for the heme biosynthesis pathway.

Plastome of mycoheterotrophic Burmannia itoana Mak. (Burmanniaceae) exhibits extensive degradation and distinct rearrangements

Plastomes of heterotrophs went through varying degrees of degradation along with the transition from autotrophic to heterotrophic lifestyle. Here, we identified the plastome of mycoheterotrophic species Burmannia itoana and compared it with those of its reported relatives including three autotrophs and one heterotroph (Thismia tentaculata) in Dioscoreales. B. itoana yields a rampantly degraded plastome reduced in size and gene numbers at the advanced stages of degradation. Its length is 44,463 bp with a quadripartite structure. B. itoana plastome contains 33 tentatively functional genes and six tentative pseudogenes, including several unusually retained genes. These unusual retention suggest that the inverted repeats (IRs) regions and possibility of being compensated may prolong retention of genes in plastome at the advanced stage of degradation. Otherwise, six rearrangements including four inversions (Inv1/Inv2/Inv3/Inv4) and two translocations (Trans1/Trans2) were detected in B. itoana plastome vs. its autotrophic relative B. disticha. We speculate that Inv1 may be mediated by recombination of distinct tRNA genes, while Inv2 is likely consequence of extreme gene losses due to the shift to heterotrophic lifestyle. The other four rearrangements involved in IRs and small single copy region may attribute to multiple waves of IRs and overlapping inversions. Our study fills the gap of knowledge about plastomes of heterotroph in Burmannia and provides a new evidence for the convergent degradation patterns of plastomes en route to heterotrophic lifestyle.

Complete plastome sequencing of Allium paradoxum reveals unusual rearrangements and the loss of the ndh genes as compared to Allium ursinum and other onions

In this work the complete chloroplast DNAs of Allium paradoxum and Allium ursinum, two edible species of Allium subg. Amerallium (the first lineage), were sequenced, assembled, annotated, and compared with complete Allium plastomes of the second and third evolutionary lines from GenBank database. The A. ursinum plastome contains 90 predicted genes (81 unique) including 5 pseudogenes, while A. paradoxum has 88 predicted genes (79 unique) including 19 pseudogenes. The comparative analysis has revealed that the A. paradoxum plastome differs markedly from those of other species. Due to many deletions, the A. paradoxum plastome is the shortest of known for Allium species, being only 145,819 bp long. The most prominent distinctions are (1) a 4825 bp long local inversion that spans from the ndhE to the rpl32 gene in the small single copy region and (2) pseudogenization, or the loss of all NADH-genes. In contrast, the plastome of A. ursinum - a species from the first evolutionary line (as well as A. paradoxum) - resembles the Allium species of the second and third evolutionary lines, showing no large rearrangements or discrepancies in gene content. It is unclear yet whether only A. paradoxum was affected by some evolutionary events or its close relatives from both sect. Briseis and other sections of Amerallium were altered as well. We speculate the sunlight-intolerant, shade-loving nature of A. paradoxum and the impairment of the ndh genes in its plastome could be interrelated phenomena.

Rhopalocnemis phalloides has one of the most reduced and mutated plastid genomes known

Although most plant species are photosynthetic, several hundred species have lost the ability to photosynthesize and instead obtain nutrients via various types of heterotrophic feeding. Their plastid genomes markedly differ from the plastid genomes of photosynthetic plants. In this work, we describe the sequenced plastid genome of the heterotrophic plant Rhopalocnemis phalloides, which belongs to the family Balanophoraceae and feeds by parasitizing other plants. The genome is highly reduced (18,622 base pairs vs. approximately 150 kbp in autotrophic plants) and possesses an extraordinarily high AT content, 86.8%, which is inferior only to AT contents of plastid genomes of Balanophora, a genus from the same family. The gene content of this genome is quite typical of heterotrophic plants, with all of the genes related to photosynthesis having been lost. The remaining genes are notably distorted by a high mutation rate and the aforementioned AT content. The high AT content has led to sequence convergence between some of the remaining genes and their homologs from AT-rich plastid genomes of protists. Overall, the plastid genome of R. phalloides is one of the most unusual plastid genomes known.

Unprecedented Parallel Photosynthetic Losses in a Heterotrophic Orchid Genus

Heterotrophic plants are evolutionary experiments in genomic, morphological, and physiological change. Yet, genomic sampling gaps exist among independently derived heterotrophic lineages, leaving unanswered questions about the process of genome modification. Here, we have sequenced complete plastid genomes for all species of the leafless orchid genus Hexalectris, including multiple individuals for most, and leafy relatives Basiphyllaea and Bletia. Our objectives are to determine the number of independent losses of photosynthesis and to test hypotheses on the process of genome degradation as a result of relaxed selection. We demonstrate four to five independent losses of photosynthesis in Hexalectris based on degradation of the photosynthetic apparatus, with all but two species displaying evidence of losses, and variation in gene loss extending below the species level. Degradation in the atp complex is advanced in Hexalectris warnockii, whereas only minimal degradation (i.e., physical loss) has occurred among some "housekeeping" genes. We find genomic rearrangements, shifts in Inverted Repeat boundaries including complete loss in one accession of H. arizonica, and correlations among substitutional and genomic attributes. Our unprecedented finding of multiple, independent transitions to a fully mycoheterotrophic lifestyle in a single genus reveals that the number of such transitions among land plants is likely underestimated. This study underscores the importance of dense taxon sampling, which is highly informative for advancing models of genome evolution in heterotrophs. Mycoheterotrophs such as Hexalectris provide forward-genetic opportunities to study the consequences of radical genome evolution beyond what is possible with mutational studies in model organisms alone.

There Is Treasure Everywhere: Reductive Plastid Evolution in Apicomplexa in Light of Their Close Relatives

The phylum Apicomplexa (Alveolates) comprises a group of host-associated protists, predominately intracellular parasites, including devastating parasites like Plasmodium falciparum, the causative agent of malaria. One of the more fascinating characteristics of Apicomplexa is their highly reduced (and occasionally lost) remnant plastid, termed the apicoplast. Four core metabolic pathways are retained in the apicoplast: heme synthesis, iron–sulfur cluster synthesis, isoprenoid synthesis, and fatty acid synthesis. It has been suggested that one or more of these pathways are essential for plastid and plastid genome retention. The past decade has witnessed the discovery of several apicomplexan relatives, and next-generation sequencing efforts are revealing that they retain variable plastid metabolic capacities. These data are providing clues about the core genes and pathways of reduced plastids, while at the same time further confounding our view on the evolutionary history of the apicoplast. Here, we examine the evolutionary history of the apicoplast, explore plastid metabolism in Apicomplexa and their close relatives, and propose that the differences among reduced plastids result from a game of endosymbiotic roulette. Continued exploration of the Apicomplexa and their relatives is sure to provide new insights into the evolution of the apicoplast and apicomplexans as a whole.

Small Genomes and Big Data: Adaptation of Plastid Genomics to the High-Throughput Era

Plastid genome sequences are becoming more readily available with the increase in high-throughput sequencing, and whole-organelle genetic data is available for algae and plants from across the diversity of photosynthetic eukaryotes. This has provided incredible opportunities for studying species which may not be amenable to in vivo study or genetic manipulation or may not yet have been cultured. Research into plastid genomes has pushed the limits of what can be deduced from genomic information, and in particular genomic information obtained from public databases. In this Review, we discuss how research into plastid genomes has benefitted enormously from the explosion of publicly available genome sequence. We describe two case studies in how using publicly available gene data has supported previously held hypotheses about plastid traits from lineage-restricted experiments across algal and plant diversity. We propose how this approach could be used across disciplines for inferring functional and biological characteristics from genomic approaches, including integration of new computational and bioinformatic approaches such as machine learning. We argue that the techniques developed to gain the maximum possible insight from plastid genomes can be applied across the eukaryotic tree of life.

Discovering the in vitro potent inhibitors against Babesia and Theileria parasites by repurposing the Malaria Box: A review

There is an innovative approach to discovering and developing novel potent and safe anti-Babesia and anti-Theileria agents for the control of animal piroplasmosis. Large-scale screening of 400 compounds from a Malaria Box (a treasure trove of 400 diverse compounds with antimalarial activity has been established by Medicines for Malaria Venture) against the in vitro growth of bovine Babesia and equine Babesia and Theileria parasites was performed, and the data were published in a brief with complete dataset from 236 screens of the Malaria Box compounds. Therefore, in this review, we explored and discussed in detail the in vitro inhibitory effects of 400 antimalarial compounds (200 drug-like and 200 probe-like) from the Malaria Box against Babesia (B.) bovis, B. bigemina, B. caballi, and Theileria (T.) equi. Seventeen hits were the most interesting with regard to bovine Babesia parasites, with mean selectivity indices (SIs) greater than 300 and half maximal inhibitory concentration (IC₅₀s) ranging from 50 to 410 nM. The most interesting compounds with regard to equine Babesia and Theileria parasites were MMV020490 and MMV020275, with mean SIs > 258.68 and >251.55, respectively, and IC₅₀s ranging from 76 to 480 nM. Ten novel anti-B. bovis, anti-B. bigemina, anti-T. equi, and anti-B. caballi hits, MMV666093, MMV396794, MMV006706, MMV665941, MMV085203, MMV396693, MMV006787, MMV073843, MMV007092, and MMV665875, with nanomole levels of IC₅₀ were identified. The most interesting hits were MMV396693, MMV073843, MMV666093, and MMV665875, with mean SIs greater than 307.8 and IC₅₀s ranging from 43 to 630 nM for both bovine Babesia and equine Babesia and Theileria parasites. Screening the Malaria Box against the in vitro growth of Babesia and Theileria parasites helped with the discovery of new drugs than those traditionally used, diminazene aceturate and imidocarb dipropionate, and indicated the potential of the Malaria Box in finding new, potent antibabesial drugs.

Prediction and large-scale analysis of primary operons in plastids reveals unique genetic features in the evolution of chloroplasts

While bacterial operons have been thoroughly studied, few analyses of chloroplast operons exist, limiting the ability to study fundamental elements of these structures and utilize them for synthetic biology. Here, we describe the creation of a plastome-specific operon database (link provided below) achieved by combining experimental tools and predictive modeling. Using a Reverse-Transcription-PCR based method and published data, we determined the transcription-state of 213 gene pairs from four plastomes of evolutionary distinct organisms. By analyzing sequence-based features computed for our dataset, we were able to highlight fundamental characteristics differentiating between operon pairs and non-operon pairs. These include an interesting tendency toward maintaining similar messenger RNA-folding profiles in operon gene pairs, a feature that failed to yield any informative separation in cyanobacteria, suggesting that it catches unique traits of operon gene expression, which have evolved post-endosymbiosis. Subsequently, we used this feature set to train a random-forest classifier for operon prediction. As our results demonstrate the ability of our predictor to obtain accurate (84%) and robust predictions on unlabeled datasets, we proceeded to building operon maps for 2018 sequenced plastids. Our database may now present new opportunities for promoting metabolic engineering and synthetic biology in chloroplasts.

Novel genetic code and record-setting AT-richness in the highly reduced plastid genome of the holoparasitic plant Balanophora

Plastid genomes (plastomes) vary enormously in size and gene content among the many lineages of nonphotosynthetic plants, but key lineages remain unexplored. We therefore investigated plastome sequence and expression in the holoparasitic and morphologically bizarre Balanophoraceae. The two Balanophora plastomes examined are remarkable, exhibiting features rarely if ever seen before in plastomes or in any other genomes. At 15.5 kb in size and with only 19 genes, they are among the most reduced plastomes known. They have no tRNA genes for protein synthesis, a trait found in only three other plastid lineages, and thus Balanophora plastids must import all tRNAs needed for translation. Balanophora plastomes are exceptionally compact, with numerous overlapping genes, highly reduced spacers, loss of all cis-spliced introns, and shrunken protein genes. With A+T contents of 87.8% and 88.4%, the Balanophora genomes are the most AT-rich genomes known save for a single mitochondrial genome that is merely bloated with AT-rich spacer DNA. Most plastid protein genes in Balanophora consist of ≥90% AT, with several between 95% and 98% AT, resulting in the most biased codon usage in any genome described to date. A potential consequence of its radical compositional evolution is the novel genetic code used by Balanophora plastids, in which TAG has been reassigned from stop to tryptophan. Despite its many exceptional properties, the Balanophora plastome must be functional because all examined genes are transcribed, its only intron is correctly trans-spliced, and its protein genes, although highly divergent, are evolving under various degrees of selective constraint.

Comparative analysis of nucleus-encoded plastid-targeting proteins in Rafflesia cantleyi against photosynthetic and non-photosynthetic representatives reveals orthologous systems with potentially divergent functions

Parasitic plants are known to discard photosynthesis thus leading to the deletion or loss of the plastid genes. Despite plastid genome reduction in non-photosynthetic plants, some nucleus-encoded proteins are transported back to the plastid to carry out specific functions. In this work, we study such proteins in Rafflesia cantleyi, a member of the holoparasitic genus well-known for producing the largest single flower in the world. Our analyses of three transcriptome datasets, two holoparasites (R. cantleyi and Phelipanche aegyptiaca) and one photosynthetic plant (Arabidopsis thaliana), suggest that holoparasites, such as R. cantleyi, retain some common plastid associated processes such as biosynthesis of amino acids and lipids, but are missing photosynthesis components that can be extensions of these pathways. The reconstruction of two selected biosynthetic pathways involving plastids correlates the trend of plastid retention to pathway complexity - transcriptome evidence for R. cantleyi suggests alternate mechanisms in regulating the plastidial heme and terpenoid backbone biosynthesis pathways. The evolution to holoparasitism from autotrophy trends towards devolving the plastid genes to the nuclear genome despite the functional sites remaining in the plastid, or maintaining non-photosynthetic processes in the plastid, before the eventual loss of the plastid and any site dependent functions.

Genome Reports: Contracted Genes and Dwarfed Plastome in Mycoheterotrophic Sciaphila thaidanica (Triuridaceae, Pandanales)

With a reduced need for photosynthesis, the plastome of parasitic and mycoheterotrophic plants degrades. In the tiny, fully mycoheterotrophic plant Sciaphila thaidanica, we find one of the smallest plastomes yet encountered. Its size is just 12,780 bp and it contains only 20 potentially functional housekeeping genes. Thus S. thaidanica fits the proposed model of gene loss in achlorophyllous plants. The most astonishing feature of the plastome is its extremely compact nature, with more than half of the genes having overlapping reading frames. Additionally, intergenic sequences have been reduced to a bare minimum, and the retained genes have been reduced in length both compared with the orthologous genes in another mycoheterotrophic species of Sciaphila and in the autotrophic relative Carludovica.

Cytonuclear integration and co-evolution

The partitioning of genetic material between the nucleus and cytoplasmic (mitochondrial and plastid) genomes within eukaryotic cells necessitates coordinated integration between these genomic compartments, with important evolutionary and biomedical implications. Classic questions persist about the pervasive reduction of cytoplasmic genomes via a combination of gene loss, transfer and functional replacement - and yet why they are almost always retained in some minimal form. One striking consequence of cytonuclear integration is the existence of 'chimeric' enzyme complexes composed of subunits encoded in two different genomes. Advances in structural biology and comparative genomics are yielding important insights into the evolution of such complexes, including correlated sequence changes and recruitment of novel subunits. Thus, chimeric cytonuclear complexes provide a powerful window into the mechanisms of molecular co-evolution.

The monoplastidic bottleneck in algae and plant evolution

Plastids in plants and algae evolved from the endosymbiotic integration of a cyanobacterium by a heterotrophic eukaryote. New plastids can only emerge through fission; thus, the synchronization of bacterial division with the cell cycle of the eukaryotic host was vital to the origin of phototrophic eukaryotes. Most of the sampled algae house a single plastid per cell and basal-branching relatives of polyplastidic lineages are all monoplastidic, as are some non-vascular plants during certain stages of their life cycle. In this Review, we discuss recent advances in our understanding of the molecular components necessary for plastid division, including those of the peptidoglycan wall (of which remnants were recently identified in moss), in a wide range of phototrophic eukaryotes. Our comparison of the phenotype of 131 species harbouring plastids of either primary or secondary origin uncovers that one prerequisite for an algae or plant to house multiple plastids per nucleus appears to be the loss of the bacterial genes minD and minE from the plastid genome. The presence of a single plastid whose division is coupled to host cytokinesis was a prerequisite of plastid emergence. An escape from such a monoplastidic bottleneck succeeded rarely and appears to be coupled to the evolution of additional layers of control over plastid division and a complex morphology. The existence of a quality control checkpoint of plastid transmission remains to be demonstrated and is tied to understanding the monoplastidic bottleneck.

On the brink: the highly reduced plastomes of nonphotosynthetic Ericaceae

Ericaceae (the heather family) is a large and diverse group of plants that forms elaborate symbiotic relationships with mycorrhizal fungi, and includes several nonphotosynthetic lineages. Using an extensive sample of fully mycoheterotrophic (MH) species, we explored inter- and intraspecific variation as well as selective constraints acting on the plastomes of these unusual plants. The plastomes of seven MH genera were analysed in a phylogenetic context with two geographically disparate individuals sequenced for Allotropa, Monotropa, and Pityopus. The plastomes of nonphotosynthetic Ericaceae are highly reduced in size (c. 33-41 kbp) and content, having lost all photosynthesis-related genes, and are reduced to encoding housekeeping genes as well as a protease subunit (clpP)-like and acetyl-CoA carboxylase subunit D (accD)-like open reading frames. Despite an increase in the rate of their nucleotide substitutions, the remaining protein-coding genes are typically under purifying selection in full MHs. We also identified ribosomal proteins under relaxed or neutral selection. These plastomes also exhibit striking structural rearrangements. Intraspecific variation within MH Ericaceae ranges from a few differences (Allotropa) to extensive population divergences (Monotropa, Hypopitys), which indicates that cryptic speciation may be occurring in several lineages. The pattern of gene loss within fully MH Ericaceae plastomes suggests an advanced state of degradation.

Structure of Plastid Genomes of Photosynthetic Eukaryotes

This review presents current views on the plastid genomes of higher plants and summarizes data on the size, structural organization, gene content, and other features of plastid DNAs. Special emphasis is placed on the properties of organization of land plant plastid genomes (nucleoids) that distinguish them from bacterial genomes. The prospects of genetic engineering of chloroplast genomes are discussed.

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.

Plastomes on the edge: the evolutionary breakdown of mycoheterotroph plastid genomes

Contents 48 I. 48 II. 50 III. 53 54 References 54 SUMMARY: We examine recent evidence for ratchet-like genome degradation in mycoheterotrophs, plants that obtain nutrition from fungi. Initial loss of the NADH dehydrogenase-like (NDH) complex may often set off an irreversible evolutionary cascade of photosynthetic gene losses. Genes for plastid-encoded subunits of RNA polymerase and photosynthetic enzymes with secondary functions (Rubisco and ATP synthase) can persist initially, with nonsynchronous and quite broad windows in the relative timing of their loss. Delayed losses of five core nonbioenergetic genes (especially trnE and accD, which respectively code for glutamyl tRNA and a subunit of acetyl-CoA carboxylase) probably explain long-term persistence of heterotrophic plastomes. The observed range of changes of mycoheterotroph plastomes is similar to that of holoparasites, although greater diversity of both probably remains to be discovered. These patterns of gene loss/retention can inform research programs on plastome function.

Degeneration of photosynthetic capacity in mixotrophic plants, Chimaphila japonica and Pyrola decorata (Ericaceae)

The evolution of photosynthesis is an important feature of mixotrophic plants. Previous inferences proposed that mixotrophic taxa tend to retain most genes relating to photosynthetic functions but vary in plastid gene content. However, no sequence data are available to test this hypothesis in Ericaceae. To investigate changes in plastid genomes that may result from a transition from autotrophy to mixotrophy, the plastomes of two mixotrophic plants, Pyrola decorata and Chimaphila japonica, were sequenced at Illumina's Genome Analyzer and compared to the published plastome of the autotrophic plant Rhododendron simsii, which also belongs to Ericaceae. The greatest discrepancy between mixotrophic and autotrophic plants was that ndh genes for both P. decorata and C. japonica plastomes have nearly all become pseudogenes. P. decorata and C. japonica also retained all genes directly involved in photosynthesis under strong selection. The calculated rate of nonsynonymous nucleotide substitutions and synonymous substitutions of protein-coding genes (dN/dS) showed that substitution rates in shade plants were apparently higher than those in sunlight plants. The two mixotrophic plastomes were generally very similar to that of non-parasitic plants, although ndh genes were largely pseudogenized. Photosynthesis genes under strong selection were retained in the two mixotrophs, however, with greatly increased substitution rates. Further research is needed to gain a clearer understanding of the evolution of autotrophy and mixotrophy in Ericaceae.

Effect of light and prey availability on gene expression of the mixotrophic chrysophyte, Ochromonas sp

Background

Ochromonas is a genus of mixotrophic chrysophytes that is found ubiquitously in many aquatic environments. Species in this genus can be important consumers of bacteria but vary in their ability to perform photosynthesis. We studied the effect of light and bacteria on growth and gene expression of a predominantly phagotrophic Ochromonas species. Axenic cultures of Ochromonas sp. were fed with heat-killed bacteria (HKB) and grown in constant light or darkness. RNA was extracted from cultures in the light or in the dark with HKB present (Light + HKB; Dark + HKB), and in the light after HKB were depleted (Light + depleted HKB).

Results

There were no significant differences in the growth or bacterial ingestion rates between algae grown in light or dark conditions. The availability of light led to a differential expression of only 8% of genes in the transcriptome. A number of genes associated with photosynthesis, phagotrophy, and tetrapyrrole synthesis was upregulated in the Light + HKB treatment compared to Dark + HKB. Conversely, the comparison between the Light + HKB and Light + depleted HKB treatments revealed that the presence of HKB led to differential expression of 59% of genes, including the majority of genes involved in major carbon and nitrogen metabolic pathways. Genes coding for unidirectional enzymes for the utilization of glucose were upregulated in the presence of HKB, implying increased glycolytic activities during phagotrophy. Algae without HKB upregulated their expression of genes coding for ammonium transporters, implying uptake of inorganic nitrogen from the culture medium when prey were unavailable.

Conclusions

Transcriptomic results agreed with previous observations that light had minimal effect on the population growth of Ochromonas sp. However, light led to the upregulation of a number of phototrophy- and phagotrophy-related genes, while the availability of bacterial prey led to prominent changes in major carbon and nitrogen metabolic pathways. Our study demonstrated the potential of transcriptomic approaches to improve our understanding of the trophic physiologies of complex mixotrophs, and revealed responses in Ochromonas sp. not apparent from traditional culture studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3549-1) contains supplementary material, which is available to authorized users.

Evolution of the Tetrapyrrole Biosynthetic Pathway in Secondary Algae: Conservation, Redundancy and Replacement

Tetrapyrroles such as chlorophyll and heme are indispensable for life because they are involved in energy fixation and consumption, i.e. photosynthesis and oxidative phosphorylation. In eukaryotes, the tetrapyrrole biosynthetic pathway is shaped by past endosymbioses. We investigated the origins and predicted locations of the enzymes of the heme pathway in the chlorarachniophyte Bigelowiella natans, the cryptophyte Guillardia theta, the “green” dinoflagellate Lepidodinium chlorophorum, and three dinoflagellates with diatom endosymbionts (“dinotoms”): Durinskia baltica, Glenodinium foliaceum and Kryptoperidinium foliaceum. Bigelowiella natans appears to contain two separate heme pathways analogous to those found in Euglena gracilis; one is predicted to be mitochondrial-cytosolic, while the second is predicted to be plastid-located. In the remaining algae, only plastid-type tetrapyrrole synthesis is present, with a single remnant of the mitochondrial-cytosolic pathway, a ferrochelatase of G. theta putatively located in the mitochondrion. The green dinoflagellate contains a single pathway composed of mostly rhodophyte-origin enzymes, and the dinotoms hold two heme pathways of apparently plastidal origin. We suggest that heme pathway enzymes in B. natans and L. chlorophorum share a predominantly rhodophytic origin. This implies the ancient presence of a rhodophyte-derived plastid in the chlorarachniophyte alga, analogous to the green dinoflagellate, or an exceptionally massive horizontal gene transfer.

Understanding the evolution of holoparasitic plants: the complete plastid genome of the holoparasite Cytinus hypocistis (Cytinaceae)

Background and Aims Plant plastid genomes are highly conserved in size, gene content and structure; however, parasitic plants are a noticeable exception to this evolutionary stability. Although the evolution of parasites could help to better understand plastome evolution in general, complete plastomes of parasites have been sequenced only for some lineages so far. Here we contribute to filling this gap by providing and analysing the complete plastome sequence of Cytinus hypocistis, the first parasite sequenced for Malvales and a species suspected to have an extremely small genome. Methods We sequenced and assembled de novo the plastid genome of Cytinus hypocistis using a shotgun approach on genomic DNA. Phylogenomic analyses based on coding regions were performed on Malvidae. For each coding region present in Cytinus, we tested for relaxation or intensification of selective pressures in the Cytinus lineage compared with autotrophic Malvales. Key Results Cytinus hypocistis has an extremely divergent genome that is among the smallest sequenced to date (19·4 kb), with only 23 genes and no inverted repeat regions. Phylogenomic analysis confirmed the position of Cytinus within Malvales. All coding regions of Cytinus plastome presented very high substitution rates compared with non-parasitic Malvales. Conclusions Some regions were inferred to be under relaxed negative selection in Cytinus, suggesting that further plastome reduction is occurring due to relaxed purifying selection associated with the loss of photosynthetic activity. On the other hand, increased selection intensity and strong positive selection were detected for rpl22 in the Cytinus lineage, which might indicate an evolutionary role in the host-parasite arms race, a point that needs further research.

Parallel evolution of highly conserved plastid genome architecture in red seaweeds and seed plants

Background

The red algae (Rhodophyta) diverged from the green algae and plants (Viridiplantae) over one billion years ago within the kingdom Archaeplastida. These photosynthetic lineages provide an ideal model to study plastid genome reduction in deep time. To this end, we assembled a large dataset of the plastid genomes that were available, including 48 from the red algae (17 complete and three partial genomes produced for this analysis) to elucidate the evolutionary history of these organelles.

Results

We found extreme conservation of plastid genome architecture in the major lineages of the multicellular Florideophyceae red algae. Only three minor structural types were detected in this group, which are explained by recombination events of the duplicated rDNA operons. A similar high level of structural conservation (although with different gene content) was found in seed plants. Three major plastid genome architectures were identified in representatives of 46 orders of angiosperms and three orders of gymnosperms.

Conclusions

Our results provide a comprehensive account of plastid gene loss and rearrangement events involving genome architecture within Archaeplastida and lead to one over-arching conclusion: from an ancestral pool of highly rearranged plastid genomes in red and green algae, the aquatic (Florideophyceae) and terrestrial (seed plants) multicellular lineages display high conservation in plastid genome architecture. This phenomenon correlates with, and could be explained by, the independent and widely divergent (separated by >400 million years) origins of complex sexual cycles and reproductive structures that led to the rapid diversification of these lineages.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0299-5) contains supplementary material, which is available to authorized users.

Lineage-Specific Reductions of Plastid Genomes in an Orchid Tribe with Partially and Fully Mycoheterotrophic Species

The plastid genome (plastome) of heterotrophic plants like mycoheterotrophs and parasites shows massive gene losses in consequence to the relaxation of functional constraints on photosynthesis. To understand the patterns of this convergent plastome reduction syndrome in heterotrophic plants, we studied 12 closely related orchids of three different lifeforms from the tribe Neottieae (Orchidaceae). We employ a comparative genomics approach to examine structural and selectional changes in plastomes within Neottieae. Both leafy and leafless heterotrophic species have functionally reduced plastid genome. Our analyses show that genes for the NAD(P)H dehydrogenase complex, the photosystems, and the RNA polymerase have been lost functionally multiple times independently. The physical reduction proceeds in a highly lineage-specific manner, accompanied by structural reconfigurations such as inversions or modifications of the large inverted repeats. Despite significant but minor selectional changes, all retained genes continue to evolve under purifying selection. All leafless Neottia species, including both visibly green and nongreen members, are fully mycoheterotrophic, likely evolved from leafy and partially mycoheterotrophic species. The plastomes of Neottieae span many stages of plastome degradation, including the longest plastome of a mycoheterotroph, providing invaluable insights into the mechanisms of plastome evolution along the transition from autotrophy to full mycoheterotrophy.

Comparative analysis of plastid genomes of non-photosynthetic Ericaceae and their photosynthetic relatives

Although plastid genomes of flowering plants are typically highly conserved regarding their size, gene content and order, there are some exceptions. Ericaceae, a large and diverse family of flowering plants, warrants special attention within the context of plastid genome evolution because it includes both non-photosynthetic and photosynthetic species with rearranged plastomes and putative losses of "essential" genes. We characterized plastid genomes of three species of Ericaceae, non-photosynthetic Monotropa uniflora and Hypopitys monotropa and photosynthetic Pyrola rotundifolia, using high-throughput sequencing. As expected for non-photosynthetic plants, M. uniflora and H. monotropa have small plastid genomes (46 kb and 35 kb, respectively) lacking genes related to photosynthesis, whereas P. rotundifolia has a larger genome (169 kb) with a gene set similar to other photosynthetic plants. The examined genomes contain an unusually high number of repeats and translocations. Comparative analysis of the expanded set of Ericaceae plastomes suggests that the genes clpP and accD that are present in the plastid genomes of almost all plants have not been lost in this family (as was previously thought) but rather persist in these genomes in unusual forms. Also we found a new gene in P. rotundifolia that emerged as a result of duplication of rps4 gene.