- Research article
- Open Access
A non-photosynthetic green alga illuminates the reductive evolution of plastid electron transport systems
BMC Biology volume 18, Article number: 126 (2020)
Plastid electron transport systems are essential not only for photosynthesis but also for dissipating excess reducing power and sinking excess electrons generated by various redox reactions. Although numerous organisms with plastids have lost their photoautotrophic lifestyles, there is a spectrum of known functions of remnant plastids in non-photosynthetic algal/plant lineages; some of non-photosynthetic plastids still retain diverse metabolic pathways involving redox reactions while others, such as apicoplasts of apicomplexan parasites, possess highly reduced sets of functions. However, little is known about underlying mechanisms for redox homeostasis in functionally versatile non-photosynthetic plastids and thus about the reductive evolution of plastid electron transport systems.
Here we demonstrated that the central component for plastid electron transport systems, plastoquinone/plastoquinol pool, is still retained in a novel strain of an obligate heterotrophic green alga lacking the photosynthesis-related thylakoid membrane complexes. Microscopic and genome analyses revealed that the Volvocales green alga, chlamydomonad sp. strain NrCl902, has non-photosynthetic plastids and a plastid DNA that carries no genes for the photosynthetic electron transport system. Transcriptome-based in silico prediction of the metabolic map followed by liquid chromatography analyses demonstrated carotenoid and plastoquinol synthesis, but no trace of chlorophyll pigments in the non-photosynthetic green alga. Transient RNA interference knockdown leads to suppression of plastoquinone/plastoquinol synthesis. The alga appears to possess genes for an electron sink system mediated by plastid terminal oxidase, plastoquinone/plastoquinol, and type II NADH dehydrogenase. Other non-photosynthetic algae/land plants also possess key genes for this system, suggesting a broad distribution of an electron sink system in non-photosynthetic plastids.
The plastoquinone/plastoquinol pool and thus the involved electron transport systems reported herein might be retained for redox homeostasis and might represent an intermediate step towards a more reduced set of the electron transport system in many non-photosynthetic plastids. Our findings illuminate a broadly distributed but previously hidden step of reductive evolution of plastid electron transport systems after the loss of photosynthesis.
Multiple lineages possess plastids acquired through independent endosymbioses in eukaryotes [1, 2]. An engulfed cyanobacterium by a heterotrophic eukaryote evolved to form the first plastid, an evolutionary event that occurred in the last common ancestor of Archaeplastida, which comprises land plants, green algae, red algae, and glaucophytes. Other algal lineages possess red alga-derived or green alga-derived complex plastids acquired by secondary or tertiary endosymbiotic events. The fact that plastids have been laterally transferred multiple times through eukaryote-eukaryote endosymbioses has spawned much interest in untangling the evolutionary events around the diversity of extant plastids.
Photosynthesis provides algae and plants with biochemical energy in forms of ATP and NADPH by conversion of solar energy through the plastid ATP synthase complex following photosynthetic electron transport accomplished by the cytochrome b6/f complex, plastoquinone/plastoquinol (PQ/PQH2), plastocyanin, ferredoxin (Fd), Fd:NADP+ reductase (FNR), and photosynthetic pigment systems, i.e., photosystems (PS) I and II [3,4,5,6]. The generated NADPH and ATP are utilized for carbon fixation through the Calvin Benson cycle as well as for various other metabolic pathways that are essential in the photosynthetic plastid. Regardless of the beneficial aspects of photosynthesis that allows autotrophic lifestyles, some species of algal and plant lineages have lost their photoautotrophic lifestyles secondarily. Almost all photosynthetic lineages appear to include such secondary heterotrophs, indicating multiple losses of photoautotrophic lifestyles independently in eukaryotic evolution .
The most well-studied species among such secondary heterotrophs is the malaria parasite Plasmodium falciparum (Apicomplexa), which bears a non-photosynthetic plastid called an apicoplast, with only a few metabolic functions such as biosynthesis of heme, Fe-S cluster, fatty acids, and isopentenyl pyrophosphates . It also retains the most reduced system of photosynthetic electron transport comprising only Fd and FNR [9, 10]. However, recent investigations have greatly expanded the spectrum of known functions of remnant plastids in non-photosynthetic algal/plant lineages [11,12,13,14,15,16,17,18,19] and have been reviewed in recent papers [20, 21]. For example, the non-photosynthetic diatom plastid in Nitzschia sp. NIES-3581 still retains multiple redox reactions for glycolysis, pentose phosphate pathway, and biosynthesis of a variety of amino acids, in addition to the functions found in the apicoplast, except for isopentenyl pyrophosphate synthesis . The volvocales green algae Polytomella spp. and the trebouxiophycean green alga Helicosporidium sp. also possess similarly complex metabolisms in non-photosynthetic plastids [11, 12]. In contrast, biosynthetic pathways for certain amino acids and fatty acids are reported to be lost in non-photosynthetic plastids of the chrysophycean “Spumella” sp. NIES-1846 bearing only glycolysis and biosynthesis of heme and Fe-S cluster . The newly discovered sister lineage of red algae, Rhodelphidia, possesses a non-photosynthetic plastid metabolically functioning only for the synthesis of heme and Fe-S clusters . Nuclear-encoded proteins imported across plastid membranes after translation in the cytosol are responsible for all the functions introduced above [8,9,10,11,12,13,14,15,16,17,18,19].
In contrast to the metabolic functions, few studies have focused on the evolution of the plastid electron transport system in non-photosynthetic plastids. In addition to photosynthetic linear electron transport, multiple branched pathways of electron transport mediated by PQ/PQH2 and plastid terminal oxidase (PTOX) are equipped in photosynthetic plastids as particular components to avoid the fatal photodamage caused by inefficient linear electron transport due to an imbalanced ratio of reducing power and ATP [22,23,24,25]. In particular, chlororespiration is a system in which type II NADH dehydrogenase (NDH2) re-oxidizes the reducing power NADPH from which electrons are discarded through PQ/PQH2 and PTOX in photosynthetic green algae [22,23,24,25]. Chlororespiration is also likely to contribute to redox homeostasis as a safety valve when excess NADPH is generated by plastid biochemical reactions . The system mediated by PQ/PQH2 and PTOX also contributes as a sink for electrons generated through the carotenoid biosynthetic pathway [22,23,24,25].
It is noteworthy that some non-photosynthetic plastids still retain numerous biochemical functions involved in redox reactions using NADP+/NADPH and involved in the carotenoid biosynthetic pathway (e.g., [13, 16]). In turn, this suggests that the electron transport systems for redox homeostasis that dissipate excess reducing power and sink excess electrons might also be dedicated to certain functions of non-photosynthetic plastids. However, underlying mechanisms by which metabolically versatile non-photosynthetic plastids still conduct numerous redox reactions remain to be elucidated.
In this study, we demonstrated that a novel strain of a heterotrophic green alga, chlamydomonad sp. NrCl902, retains biosynthesis of PQ/PQH2. A transient RNA interference (RNAi) knockdown experiment strongly suggests that the product of homogentisate transferase gene detected in this study is involved in the synthesis of the PQ/PQH2 pool of this alga. In addition to sequences for Fd-FNR, the transcriptome analysis detected sequences for the plastidal, PQ/PQH2-mediated electron sink systems such as PTOX, which are likely to be required for carotenoid biosynthesis and redox homeostasis in the non-photosynthetic plastid. Our in-depth survey further indicates a broad distribution of the gene set for plastid electron transport systems that are simpler than that of photosynthetic plastids but more complex than that of apicoplasts. Our study thus unveiled a previously hidden step for reductive evolution of photosynthetic electron transport systems, along with the evolution of photoautotrophic algae to heterotrophic protists.
Results and discussion
A novel lineage of plastid-bearing, non-photosynthetic Volvocales species
We established and maintained the axenic, clonal culture of chlamydomonad sp. strain NrCl902, a colorless Volvocales green alga (Fig. 1a), in a medium containing sodium acetate as the sole carbon/energy source (see “Methods” section for details). The cell possesses an orange eyespot (Fig. 1b), located in plastids of photosynthetic Volvocales green algae . Indeed, a membrane-bound structure containing starch granules but lacking accumulation of thylakoid membranes was observed via transmission electron microscopy (Fig. 1c). In light of these findings, we conclude that the chlamydomonad sp. possesses a non-photosynthetic plastid that is capable of synthesizing starch granules.
After SYBR Green staining, we could observe organellar DNAs in addition to the nuclear DNA (Fig. 1d, e). Most of organellar DNA signals were derived from mitochondrial DNAs as they colocalized with the Mitotracker fluorescence (Fig. 1d–g). However, some organellar DNA signals were outside the Mitotracker fluorescence (Fig. 1g), strongly suggesting some of them are derived from a plastid DNA in the new non-photosynthetic green alga.
In 18S rRNA gene phylogeny, chlamydomonad sp. NrCl902 was branched with Chlamydomonas pseudoplanoconvexa (Genbank accession number: AB602849), a photosynthetic species, with the highest bootstrap support and the highest PhyloBayes posterior probability (Fig. 2a). Other non-photosynthetic Volvocales green algae, Polytoma spp. and Polytomella spp., were distantly related to chlamydomonad sp. NrCl902. We removed Polytoma oviformis U22936  from our analysis as the nucleotide sequence had been reported to be a chimera . Thus, the chlamydomonad sp. NrCl902 is the third independent lineage of non-photosynthetic Volvocales species.
For the first time, we succeeded in assembling the complete plastid genome in non-photosynthetic Volvocales green algae, given a fragmented assembly of the previously published plastid genome of the non-photosynthetic Volvocales green alga Polytomella uvella . The 176-kb-long, circularly mapping plastid genome of chlamydomonad sp. NrCl902 does not carry the genes or pseudogenes for photosynthesis-related thylakoid membrane complexes, the carbon fixation pathway, and chlorophyll biosynthesis (Additional file 1: Fig. S1). The proteins encoded in the chlamydomonad sp. NrCl902 plastid are related to proteolysis, transcription, translation, protein transport, and plastid division. Notably, the protein-coding gene repertoire of NrCl902 is completely identical with that of P. uvella (highlighted by gray back in Fig. 2b; Additional file 1: Fig. S1). As P. uvella and our strain NrCl902 have lost photosynthesis independently (Fig. 2a), the shared gene repertoire represents the convergent reductive plastid genome evolution after loss of photosynthesis.
Plastoquinone/plastoquinol pool in the functionally versatile non-photosynthetic plastid
Our prediction of potential plastid metabolic functions in the non-photosynthetic green alga reveals a previously undiscovered function of non-photosynthetic plastids in this organism. We reconstructed the metabolic pathways of plastids in chlamydomonad sp. by transcriptome analyses with C. reinhardtii chloroplast functions as references; these included nitrite assimilation, sulfur assimilation, biosynthesis of isoprenoids, starches, fatty acids, glycerolipids, heme, chlorophylls, carotenoids, quinones, Fe-S clusters, and various amino acids, as well as carbon fixation through the Calvin Benson cycle . By reciprocal blast analyses, we detected contigs encoding plastid-targeted proteins that are involved in most of the above functions, except for photosynthesis-related functions such as photosynthetic thylakoid membrane complexes, chlorophyll biosynthesis, and carbon fixation (Additional file 2: Fig. S2; Additional file 3: Table S1). It is possible that the plastid pentose phosphate pathway found in chlamydomonad sp. NrCl902 works for supply of ATP, NADPH, and erythrose 4-phosphate, the latter of which is then utilized in the shikimate pathway for aromatic amino acid biosynthesis (Additional file 2: Fig. S2), as proposed in the non-photosynthetic diatom plastid . These complex metabolic pathways in the non-photosynthetic plastid are most likely fueled by cytosolic/mitochondrial metabolisms through plastid transporters such as triose phosphate transporters and glucose phosphate transporters (Additional file 4: Fig. S3; Additional file 5: Fig. S4; Additional file 6: Table S2). In addition, the plastid triose phosphate transporter would play roles for not only import of sugar phosphates but also export of glycerate 3-phosphate (Additional file 4: Fig. S3). Interestingly, sequences for the PQ/PQH2 biosynthetic pathway that follows isoprenoid synthesis were present in the transcriptome data (Fig. 3a in detail). To support this finding, we conducted liquid chromatography (LC)-tandem mass spectrometry (MS/MS) analyses for the detection of PQ/PQH2. In addition to ubiquinone and ubiquinol, which are the electron acceptor and donor, respectively, in mitochondrial oxidative phosphorylation (Additional file 4: Fig. S3; Additional file 6: Table S2), LC analysis detected an explicit, candidate peak of PQH2 and a faint, candidate peak of PQ (Fig. 3b). The MS spectrum confirmed that the explicit peak was PQH2-9 (Fig. 3c). After treatment of cell extracts with FeCl3 for oxidization, the PQH2-9 peak disappeared and the candidate peak of PQ became explicit, indicating that the faint peak detected by the initial analysis was the oxidized form of PQH2-9, i.e., PQ-9 (Fig. 3b). Indeed, the MS/MS spectrum indicated that the peak was most likely PQ-9 (Fig. 3d). In line with the above analyses, the non-photosynthetic plastid of chlamydomonad sp. NrCl902 was demonstrated to possess PQ-9 as the electron carrier. This is the first report of PQ/PQH2 pool in an organism with non-photosynthetic plastids that lack photosynthesis-related thylakoid membrane complexes.
To further confirm involvement of the detected sequences, such as homogentisate solanesyltransferase (HST), in biosynthesis of the PQ/PQH2 pool of this alga, we performed a transient RNAi knockdown experiment with electroporation of the double-strand RNAs (HST1; Fig. 4a) for the gene. After the electroporation, HST transcripts were undetectable by the reverse transcriptase PCR (RT-PCR) assay while transcripts of actin, a housekeeping gene for Actin cytoskeleton, were detected as a positive control (Fig. 4b). This is in stark contrast to the control experiment for which the entire experimental procedure was followed but except for electroporation without double-strand RNAs; transcripts of both genes were detected by the RT-PCR assays (Fig. 4b). The plastid PQ/PQH2 pool size relative to the UQ/UQH2 pool size of presumably mitochondria was significantly decreased after the double-strand HST1 RNA electroporation to ca. 60% of that in the control experiments (Fig. 4c), which was confirmed by the LC-MS/MS analyses. Given these, the transient knockdown of the transcripts for HST likely caused suppression of PQH2 synthesis. Interestingly, cell growth has been suppressed in the knockdown samples during 2 days after the electroporation (Fig. 4d), suggesting PQ/PQH2 deficiency might affect cell viability. To investigate a possibility of off-target effects in the cell growth, two other distinct double-strand RNA molecules (HST2 and HST3; Fig. 4a) within the same target region of HST1 were designed for independent RNAi experiments, and suppression of both HST mRNA expression and cell growth was observed as well (Fig. 4e, f). It is noteworthy that the number of cells in the knockdown samples became comparable with that of the control cells by 3 days after the electroporation (Fig. 4f) and that the RT-PCR assays for cells 3 days after the electroporation detected both HST and actin transcripts (Fig. 4g). These observations indicate that the RNAi effect in the current condition has lasted during a few days after the electroporation.
Possible roles of PQ/PQH2 pool in the non-photosynthetic plastid
PQ is known as the electron acceptor from the carotenoid biosynthetic pathway in which phytoene desaturase (PDS or CrtP) and ζ-carotene desaturase (ZDS or CrtQ) pull out electrons from carotenoid precursors . Transcriptome analyses of chlamydomonad sp. NrCl902 revealed the presence of contigs encoding carotenoid biosynthesis proteins that could convert an isoprenoid to zeaxanthin (Fig. 2a; Additional file 7: Fig. S5). Indeed, high-performance liquid chromatography (HPLC) analyses demonstrated that the chlamydomonad sp. possesses several different carotenoids including β-carotene, γ-carotene, 3-hydroxy-γ-carotene, and zeaxanthin (Additional file 7: Fig. S5). Those carotenoids might be either components of the eyespot (Fig. 1b ) or scavengers of reactive oxygen species . The non-photosynthetic green alga Polytomella magna (Fig. 2a) is known to have an eyespot and genes for the carotenoid biosynthesis pathway, although what carotenoids are synthesized remains unknown . In addition, certain carotenoids were also detected and were thought to be antioxidants in the malaria P. falciparum, although the exact localization of the carotenoid biosynthesis remained unclear in the parasite .
If PQ is utilized as an electron acceptor in the biosynthesis of carotenoids in NrCl902, a certain electron sink system should be functional in the plastid to re-oxidize PQH2 to PQ. Plastid terminal oxidase (PTOX) functions in the oxidation of PQH2 as a key component of the plastid sink of excess electrons in photosynthetic plastids [22, 24, 25] (Fig. 5a), and the transcriptome data also contain a contig encoding PTOX with a plastid transit peptide (Additional file 3: Table S1). We thus inferred that PTOX in chlamydomonad sp. NrCl902 is involved in the oxidation of PQH2 generated in association with carotenoid biosynthesis (Fig. 5b).
The PQ pool is also reduced by cyclic electron transport and by chlororespiration mediated by type II NADH dehydrogenase (NDH2) in the photosynthetic model alga C. reinhardtii [23,24,25] (Fig. 5a), in addition to carotenoid biosynthesis. We further investigated the candidate roles of PQ/PQH2 in the non-photosynthetic green alga chlamydomonad sp. NrCl902. We detected a sequence for NDH2 with a plastid transit peptide (Additional file 3: Table S1). Because the presumably plastidal NDH2 likely catalyzes the redox exchange between NADPH and PQ , PQ/PQH2 pool may also be involved in a part of the electron transport system to the terminal electron sink PTOX from NADPH. If correct, NDH2 together with PQ/PQH2 may be a component of an electron transport system for chlororespiration in the non-photosynthetic plastid of our strain NrCl902 (Fig. 5b), as chlororespiration has been proposed to be involved in NADP+/NADPH homeostasis previously .
In addition to the possible roles of the PQ/PQH2 pool as the electron sink systems for redox homeostasis in the non-photosynthetic green algal plastid, it cannot be completely ruled out the PQ/PQH2 pool has other roles. For example, it is also possible that nuclear gene expression is regulated by monitoring the redox state of the PQ/PQH2 pool as proposed in photosynthetic organisms .
Broad distribution of electron transport systems in eukaryotes with non-photosynthetic plastids
Intriguingly, the chlamydomonad sp. NrCl902 transcriptome data revealed two additional sequences related to the plastid electron transport system, plastid-targeted Fd and FNR. Although their exact roles in the non-photosynthetic green alga NrCl902 are unclear, the FNR of the non-photosynthetic plastid may reversely function to reduce Fd by NADPH-derived electrons as demonstrated in the cyanobacterial heterocyst where linear electron transport is absent  and in non-photosynthetic plastids of land plants with root-type FNR . The reverse electron transport from NADPH to Fd could explain how the non-photosynthetic plastid forms the reduced Fd required for reduced Fd-dependent nitrite and sulfate assimilation, reduction of oxidized thioredoxin by ferredoxin:thioredoxin reductase, and for Fe-S cluster assembly (Additional file 2: Fig. S2; Additional file 3: Table S1) [8, 38,39,40,41]. Otherwise, the plastid Fd-FNR system might also contribute to NADP+/NADPH homeostasis by re-oxidizing excess NADPHs  in the non-photosynthetic green algal plastid.
It was unclear how broadly genes for the plastid redox homeostasis distribute in non-photosynthetic algae/land plants, while only the plastid Fd-FNR system in apicoplasts and Rhodelphidia [9, 10, 18] and PTOX in a non-photosynthetic chrysophyte  have been reported. Our additional in-depth survey suggests that the NADPH-Fd electron transport system and the PQ/PQH2-mediated electron transport system are still retained in various non-photosynthetic algae/plants. Plastid Fd and FNR are conserved in all the examined algae/land plants with non-photosynthetic plastids, though we detected plastid flavodoxin (Fld), a functional analog of ferredoxin, in a non-photosynthetic diatom  (Fig. 5c). Therefore, the electron transport between Fd/Fld and NADP(H) mediated by FNR is highly likely one of the indispensable functions for maintaining non-photosynthetic plastids. In contrast, the sequences for NDH2, PTOX, and PQH2 synthesis show a rather punctate distribution (Fig. 5c; Additional file 8: Table S3; Additional file 9: Table S4); NDH2 sequences are detected only from green algae, sequences for PTOX and PQH2 synthesis from a diatom, a chrysophyte, a land plant and green algae. The organisms used in this comparative analysis have lost photosynthesis independently [8,9,10,11,12,13,14,15,16,17,18,19], indicating convergent evolution towards retention of genes for NDH2, PTOX, and PQH2 synthesis. Note that some of the sequences for PQ/PQH2 biosynthesis obtained from the transcriptome data lack N-terminal plastid-targeting sequences (Fig. 5c; Additional file 8: Table S3). It remains unclear whether the lack of plastid-targeting sequence in some sequences is caused by their incomplete, 5′ truncate sequences in the transcriptome data or reflects their localization outside plastids. Phylogenetic analyses of the key proteins, FNR, PTOX, and methyl-6-solanyl-1,4-benzoquinol methyltransferase (MMT) for the last step of PQH2 synthesis, show no clear evidence of lateral transfer, but rather many of the sequences are suggested to be of vertical inheritance (Additional file 10: Fig S6, Additional file 11: Fig S7, Additional file 12: Fig. S8). In addition to chlamydomonad sp. NrCl902, through vertical inheritance, multiple distinct lineages with non-photosynthetic plastids might retain or might have retained the functional, abovementioned PTOX- and PQ/PQH2-mediated electron transport system until recently (Fig. 5c). Whether these species bearing the genes are indeed capable of synthesizing PQ/PQH2 should be investigated biochemically.
Conclusion: evolutionary principle of the electron transport system
Plastids are widely known as light-dependent cellular powerhouse organelles supplying sugar through photosynthesis to the algal and plant cells [1,2,3,4,5,6]. However, research over the past decade has shown that plastids exist on a functional and evolutionary continuum that includes a variety of non-photosynthetically functioning organelles [7, 8, 11,12,13,14,15,16,17,18,19,20,21]. Nevertheless, the evolutionary transitions of the electron transport system and the branched pathways for the electron sink between photosynthetic and non-photosynthetic plastids have remained unclear. Our findings clearly demonstrate that the entire photosynthetic electron transport system (Fig. 5a) is not always straightforwardly evolved to the simplest form comprising only Fd and FNR, soon after loss of photosynthesis (Fig. 5d). Instead, some non-photosynthetic plastids have the PQ/PQH2-mediated electron transport system, and they might reflect an intermediate step for reductive evolution of the plastid electron transport system (Fig. 5a, b, and d). First, associated with loss of photosynthesis, components for the photosynthetic electron transport system are lost, i.e., PSI, PSII, cytochrome b6/f, and plastocyanin, while the PQ/PQH2-mediated electron transport system and the Fd-FNR system are retained (Fig. 5b). Subsequently, only Fd-FNR system is left after loss of PQ/PQH2-mediated electron transport system as seen in a non-photosynthetic plastid of the trebouxiophycean green alga Helicosporidium sp. (Fig. 5d). Given multiple non-photosynthetic plastid-bearing species with genes to synthesize PQH2 (Fig. 5c), plastids in the green alga NrCl902 and some other non-photosynthetic algae/plant might belong to the evolutionary stages shown in Fig. 5b. Metabolic pathways that require redox reactions such as glycolysis and/or the pentose phosphate pathway might be a key constraint against the evolutionary loss of PQ/PQH2 pool. In this point of view, the Helicosporidium plastid might have an alternative way for redox homeostasis as its genome lacks PTOX and PQ/PQH2 biosynthesis (Fig. 5c) regardless of plastid-targeted proteins for a variety of biochemical reactions . Otherwise, the complex metabolic pathways retained in the Helicosporidium plastid might be no longer functional efficiently and thus might require no efficient system for redox homeostasis. In contrast, “Spumella” sp. likely lacks various pathways but possesses genes for PTOX and some homologs for PQ/PQH2 synthesis (Fig. 5c ), suggesting multiple evolutionary routes towards loss of PTOX and PQ/PQH2 biosynthesis. Additional systematic analyses exploring the links between reductive evolution of the electron transport system and metabolisms with redox reactions in the non-photosynthetic plastids of diverse eukaryotes may provide further insights into the principle that would govern the functional reduction of plastids.
Cultivation of strain
Chlamydomonad sp. NrCl902 was isolated by a single-cell isolation technique with a glass pipet, from a freshwater sediment sample collected from a paddy field in Niinomi, Sammu, Chiba, Japan (35° 37′ 01″ N, 140° 24′ 59″ E). The cell was grown and maintained in the AFAC medium  at 25 °C under the dark condition. The culture was deposited to Microbial Culture Collection of National Institute for Environmental Studies (NIES ), Japan, as NIES-4405.
Cell shape was observed by the CCD Camera DP74-CU (Olympus) equipped with a light microscope (Zeiss). To stain DNA, SYBR Green I nucleic acid stain (Thermo Fisher Scientific) was added to the cell suspension to give a final dilution of 1:1000. To stain mitochondrial membranes, MitoTracker CMTM Ros (Thermo Fisher Scientific) were added to the cell culture to give a final concentration of 200 nmol L-1. After staining, cells were observed using BX51 fluorescence/differential interference microscope (Olympus) connected to DP72 charge-coupled device camera. Ultrastructure was observed by transmission electron microscopy as follows. Cells were frozen in liquid propane at − 175 °C and substituted with 2% glutaraldehyde, 0.5% tannic acid in acetone, and 2% distilled water at − 80 °C for 2 days. The sample was incubated at 4 °C for 2 h, following to incubation at − 20 °C for 2 h. The sample was rinsed with acetone 4 times for 15 min each and fixed with 2% osmium tetroxide (OsO4) in acetone at 4 °C for 60 min. Dehydration was performed with ethanol 3 times for 30 min each, followed by additional dehydration with ethanol at room temperature overnight. The sample was infiltrated with propylene oxide 2 times for 30 min each and put into a 70:30 mixture of propylene oxide and resin (Quetol-651; Nisshin EM Co.) for 1 h. After volatilization of propylene oxide, the sample was transferred into a fresh 100% resin, followed by polymerization at 60 °C for 48 h. The resin block was ultrathin-sectioned at 70 nm with a diamond knife using a ultramicrotome (ULTRACUT, Leica), and then the sections were placed on the copper grids. The sections were stained with 2% uranyl acetate at room temperature for 15 min, rinsed with distilled water, and secondarily stained with lead stain solution (Sigma-Aldrich Co.) at room temperature for 3 min. The sections were observed by a transmission electron microscope (JEN-1400Plus; JEOL Ltd.) at an acceleration voltage of 80 kV. Digital images (3296 × 2472 pixels) were taken with a CCD camera (EM-14830RUBY2; JEOL Ltd.).
Phylogenetic position of chlamydomonad sp. NrCl902
Nuclear small subunit rRNA gene sequence of NrCl902 was retrieved from the nuclear DNA sequencing data obtained with Illumina (see below). The rRNA sequence of NrCl902 was aligned with those of Volvocales species by MAFFT  and ambiguously aligned sites were removed with Bioedit . The dataset comprised of 54 taxa and 1726 sites was subjected to phylogenetic analyses with PhyloBayes 4.1  and IQtree 1.6.5  for Bayesian and maximum likelihood frameworks, respectively. The best-fitting available model based on the Bayesian Information Criterion was the TIM3e + I + Γ model, which was used for estimation of the maximum likelihood tree and for a bootstrap analysis with the 100 pseudoreplicates. PhyloBayes analyses were performed under the CAT-GTR + Γ model with two independent Markov chain Monte Carlo chains (MCMC) were run for 50,000 trees, sampling every 100 trees, with burnin of 12,500 trees. Two chains converged with maxdiff = 0.09. Subsequently, the consensus tree with branch lengths and Bayesian posterior probabilities (BPPs) were calculated from the rest of the sampled trees. The 18S rRNA gene sequence of NrCl902, and the dataset used for the analysis are attached in Additional file 13: Datasets.
Organellar DNA sequencing
Total DNA was extracted with Plant DNA extraction kit (Jena Biosciences) according to the manufacturer’s instruction. Total DNA was sent to Hokkaido System Science Co. for subjecting to library construction with v2 / TruSeq DNA PCR-free Sample Prep Kit (Illumina) and HiSeq2500 sequencing, resulting in 45.9 million paired-end reads. Adapter trimming and quality filtering were performed with fastX toolkit . In quality filtering, reads with quality scores > 20 for at least 75% of their length were retained, resulting in 36.9 million paired-end reads. The filtered short reads were subjected to KmerGenie 1.7044  to predict an assembled genome size which is the required parameter for the following HGAP-based assembling. DNA was also sequenced by PacBio RSII, with SMRT cell 8Pac V3 and DNA Polymerase Binding Kit P6 v2, in Macrogen, and the resultant 1.6 Gb subreads were subjected to assembling by HGAP v3  through DDBJ Pipeline . By the tblastN search  with plastid-encoded protein sequences of Chlamydomonas reinhardtii (GenBank no. BK000554), a single contig derived from a plastid genome was detected. For error correction, we mapped the filtered Illumina short reads onto the contig with Bowtie2  with default settings, and detected errors on the PacBio contig were manually corrected. By using PCR assay followed by the Sanger sequencing to fulfill the gap between the termini of the contig, we obtained the complete sequence of plastid genome with 176,432 bp in length. Protein-coding genes were identified with Mfannot  and blastX-based homology search to the nr database of GenBank . Transfer RNA genes were identified with Mfannot  and tRNAscan . Similarly, the mitochondrial contig was detected and annotated by comparison with the mitochondrial DNA sequence of C. reinhardtii (GenBank no. EU306622 ) and its close relative C. leiostraca (GenBank no. KP696389 ).
Nuclear-encoded protein sequences
Total RNA was extracted with Trizol (Thermo Fisher Scientific) according to the manufacturer’s instruction. RNA was subjected to library construction with TruSeq RNA Sample Prep Kit v2 (Illumina) and HiSeq2500 sequencing, resulting in 43.7 million paired-end reads. Adapter trimming and quality filtering were performed with fastX toolkit . In quality filtering, reads with quality scores > 20 for at least 75% of their length ≥ 50 bp long were retained, resulting in 21.8 million paired-end reads. Assembling was performed by trinity 2.4.0 with default settings. The C. reinhardtii chloroplast functions for photosynthesis; nitrite assimilation; sulfur assimilation; biosynthesis of isoprenoids, starches, fatty acids, glycerolipids, heme, chlorophylls, carotenoids, Fe-S cluster, and various amino acids; and carbon fixation through the Calvin Benson cycle were partly comprised of glycolysis and the pentose phosphate pathway . By using tblastN  with the C. reinhardtii sequences involved in the above plastid functions as queries with a cut off criterion of e− 10, we searched for homologous sequences in the transcriptome data of chlamydomonad sp. NrCl902. We then confirmed that the detected sequences are not distant paralogues by using blastP with deduced amino acid sequences as queries against the non-redundant protein database . To check whether encoded proteins in the detected sequences possess plastid-targeting transit peptides at the N-termini, ChloroP 1.1  was subjected with a cutoff value of 0.5 (Additional file 3: Table S1). We also surveyed sequences for the cytosolic glycolysis/gluconeogenesis, the acetate metabolisms, and the mitochondrial tricarboxylic cycle as described above. Mitochondrial targeting sequences were surveyed by Mitofates .
Survey of plastid-targeted proteins for electron transport systems in other non-photosynthetic algae/plant lineages
Plastid-targeted proteins for electron transport systems (Fig. 5c) were surveyed as described above for transcriptome data of Polytomella spp. [12, 13], Monotropa hypopitys , Nitzschia sp. , and “Spumella” sp. , and genome data of Helicosporidium sp. , Rhodelphis limneticus , and Plasmodium falciparum . Before the analyses, we checked the quality of the transcriptome data with gVolante . The BUSCO v2 “complete + partial” scores with 303 eukaryote conserved proteins by gVolante were 94.06% in NrCl902, 93.73% and 92.74% in Polytomella parva and Polytomella magna, respectively, and 93.4% in Monotropa. The scores for the transcriptome data of Nitzschia and Spumella used in this study were already reported to be higher than 80% [16, 17]. Detected sequences are shown in Additional file 8: Table S3 and 9: Tables S4, respectively. Plastid-targeting sequences of detected homologs were investigated as performed in previous studies [11,12,13,14,15,16,17,18]. Phylogenetic analyses of FNR, PTOX, and MMT for the last step of PQH2 synthesis were performed with IQtree 1.6.5  with 100 bootstrap analyses. Details of the datasets and the used models are described in legends for Figs. S6-S8 (Additional file 10: Fig. S6, Additional file11: Fig. S7, Additional file 12: Fig. S8). Datasets used for the analyses are attached in Additional file 13: Datasets.
Detection of carotenoids and quinones/quinols
The pigments including carotenoids were extracted with acetone/methanol (7:2, v/v) using an ultrasonicator. Each pigment was separated using a C8-HPLC column  equipped in the Separations Module Waters 2695 (Waters) and then detected by the Photodiode Array Detector Waters 2996 (Waters) and the Multi λ Fluorescence Detector Waters 2475 (Waters). Carotenoids were identified on the basis of their retention time and characteristic absorption spectra. We also confirmed the identification of each carotene with a method using a C18-HPLC column .
For quinone/quinol extraction, an aliquot of 2-propanol (LC-MS grade, Kanto Chemical Co., Inc., Tokyo, Japan) were added to a microtube containing the pelleted fresh cells of chlamydomonad sp. NrC1902 and then placed in an ice-cooled ultrasonication bath for extraction for 1 min. The supernatant was immediately separated from suspends by centrifugation. A half of the supernatant was directly injected into the HPLC apparatus for analysis. The other half was treated with ferric chloride (final concentration, 1.2 mM) before the analysis for oxidation of total quinones and quinols. The LC-MS/MS instrument was composed of a Shimadzu Nexera X2 liquid chromatography system, comprising a CBM-20A communication bus module, two DGU-20A3R/5R HPLC degassing units, three LC-30AD solvent delivery units constituting a ternary pumping system, an SIL-30 AC autosampler, a CTO-20 AC column oven, and a LC-MS-8030 triple quadrupole mass spectrometer connected through an atmospheric pressure chemical ionization (APCI) interface (Shimadzu, Kyoto, Japan). The system was coupled to a personal computer configured to run the Shimadzu LabSolution software. Reverse-phase HPLC was performed on a Zorbax Eclipse Plus C18 column (Rapid Resolution HT, 3.0 × 100 mm, 1.8 μm silica particle size; Agilent Technologies, Santa Clara, USA) with the binary gradient (flow rate of 0.5 mL min-1) of ethyl acetate (LC-MS grade, Honeywell, Seelze, Germany) with 0.1% [v/v] formic acid (LC-MS grade, Wako Pure Chemical Industries, Ltd., Osaka, Japan) in methanol (LC-MS grade, Kanto Chemical Co., Inc., Tokyo, Japan) with 0.1% [v/v] formic acid as follows (all v/v); 30% for 1.0 min, 30–80% in 4.0 min, 80% for 3 min, 80–30% in 0.1 min, and 30% for 3.9 min. All the mobile phases were degassed in vacuo with ultrasonication. The mobile-phase reservoir bottles were designed to prevent any contact between the mobile phases and air during analysis. The APCI was set as following conditions: nebulizer gas flow, 3.0 L min− 1; interface temperature, 350 °C; desolvation line temperature, 200 °C; heat block temperature, 200 °C; drying gas flow, 5 L min− 1. Parameters of the multiple reaction monitoring (MRM) in the positive ion mode of the mass spectrometer were summarized in Additional file 14 (Table S5). Parameters of the Q3 product ion scan in the positive ion mode of the mass spectrometer were summarized in Additional file 15 (Table S6).
RNAi knockdown of homogentisate solanesyltransferase
In silico prediction of suitable regions in the homogentisate solanesyltransferase (HST) transcript for RNAi was performed with siDirect 2.0 , unveiling accumulation of the suitable regions at ca. 600 bp of the 5′ terminal region (data not shown). Therefore, we decided to prepare three double-strand RNA molecules for the 5′ terminal region (Fig. 4a) as follows. One microgram of total RNA extracted from cells of chlamydomonad sp. NrCl902 with Trizol (Thermo Fisher Scientific) was subjected to cDNA synthesis with 3′ RACE System for Rapid Amplification of cDNA Ends (Thermo Fisher Scientific) according to the manufacturer’s instructions. PCR assay with the primer set, 5′-AGCCTGAATAATGGCGCAAG-3′ and 5′-TGACGAAGGCGGTGATGAAG-3′, was performed in order to obtain DNA fragments of homogentisate solanesyltransferase (HST). Using the HST DNA fragment as the template, the nested PCR was performed with primers (HST1F: 5′-CTAATACGACTCACTATAGGGAGAATAATGGCGCAAGATCAGCTTC-3′ and HST1R: 5′-CTAATACGACTCACTATAGGGAGACTTGTTCACCACGTCAATGTCC-3′) in order to prepare the DNA fragment attached with the T7 promoter at the 5′ end, which is subsequently used as the template for double-strand RNA synthesis with MEGA script RNAi Kit (Thermo Fisher Scientific) equipped with T7 RNA polymerase. The double-strand RNA molecule prepared as describe above is called HST1. To reduce the possibility of off-target effect, we additionally prepared two distinct double-strand RNA molecules, called HST2 and HST3, targeting the HST transcript (Fig. 4a). RT-PCR assays were performed using the following primer sets: HST2F: 5′-CTAATACGACTCACTATAGGGAGAAGCCTGAATAATGGCGCAAGAT-3′ and HST2R: 5′-CTAATACGACTCACTATAGGGAGACACCTGGCTGTCATTTGTGGA-3′, and HST3F: 5′-CTAATACGACTCACTATAGGGAGAAAATTCAGCCATGCGTTTTGG-3′ and HST3R: 5′-CTAATACGACTCACTATAGGGAGACCACACCGGTTGACATCTCG-3′. The synthesized PCR products were used as the templates for double-strand RNA synthesis as described above.
Chlamydomonad sp. NrCl902 were preliminarily cultivated as described above for 4 days. The cultivated cells were further incubated in a fresh AFAC medium for 24 h. Cells were collected by centrifugation at 2500 rpm for 5 min and washed with the TAP medium (Thermo Fisher Scientific) containing 40 mM sucrose. The 8.0 × 105 cells were suspended in 40 μL of the sucrose-containing TAP medium. The cell suspension was placed into an electroporation cuvette with a 2-mm gap (NEPAGENE). By using NEPA21 Super Electroporator (NEPAGENE), electroporation was performed with no double-strand RNAs or 4 μg of HST1 double-strand RNAs for nine cell suspension samples, each containing 8.0 × 105 cells in the sucrose-containing TAP medium under the following electroporation conditions reported in Yamano et al.  with some modifications. For details, parameters of poring pulse were one polarity-exchange pulse of 300 V with 8 ms pulse length, 50 ms pulse interval, and a 40% decay rate, while those of transfer pulse were a ten polarity-exchanged pulse of 20 V with 50 ms pulse length, 50 ms pulse interval, and a 40% decay rate. In this case, the measured value of electrical impedance was within 440–500 Ω in the cell conditions described above. Of nine cell suspension samples, three were subjected to cell counting, three others to the RT-PCR assays, and the others to the PQ/PQH2 detection/quantification. Electroporated cells in the TAP medium were transferred to the AFAC medium and cultivated as described above, under the dark condition. Cells were counted every 12 h under the light microscopy until 48 h after electroporation. Total RNA was extracted from cells 1.5 days after the electroporation as described above. RT-PCR assays with the total RNA as the template were conducted with the primer set for HST (5′-AGCCTGAATAATGGCGCAAG-3′ and 5′-TGACGAAGGCGGTGATGAAG-3′) or with that for actin (5′-ACTCATACGTCGGTGATGAG-3′ and 5′-GCTCCATCAAGATCTTCATC-3′). Cells 2 days after electroporation were subjected to quinone/quinol extraction as described above. The quinols in the total extract was oxidized into quinones with ferric chloride before analyses described above.
To reduce the possibility of off-target effect, six cell suspensions were subjected to the electroporation with no double-strand RNAs or each of HST1, HST2, and HST3 double-strand RNAs as described above: of six cell suspension samples, three were for the RT-PCR assays at 1.5 days after the electroporation and the others for cell counting for 3 days after the electroporation. Note that after the cell counting, the cells were used for the RT-PCR assays to check the RNAi effect after 3 days from the electroporation.
Availability of data and materials
DNA and RNA sequencing data of chlamydomonad sp. NrCl902 were deposited to DDBJ (BioProject Accessions, PRJDB10134  and PRJDB9052 , respectively). DDBJ accession numbers of the plastid and mitochondrial genomes are LC516060  and LC516061 , respectively. The culture strain NrCl902 used in this study is available in National Institute for Environmental Studies (NIES ), Japan, as NIES-4405.
Keeling PJ. The number, speed, and impact of plastid endosymbioses in eukaryotic evolution. Annu Rev Plant Biol. 2013;64:583–607.
Archibald JM. Genomic perspectives on the birth and spread of plastids. Proc Natl Acad Sci U S A. 2015;112:10147–53.
Niyogi KK. Safety valves for photosynthesis. Curr Opin Plant Biol. 2000;3:455–60.
Rochaix JD. Regulation of photosynthetic electron transport. Biochim Biophys Acta. 1807;2011:375–83.
Hohmann-Marriott MF, Blankenship RE. The photosynthetic world. In: Eaton-Rye JJ, Tripathy BC, Sharkey TD, Editors. Photosynthesis. Berlin: Springer; 2012. p. 3–32.
Cornic G, Baker NR. Electron transport in leaves: a physiological perspective. In: Eaton-Rye JJ, Tripathy BC, Sharkey TD, Editors. Photosynthesis. Berlin: Springer; 2012. p. 591-605.
Hadariová L, Vesteg M, Hampl V, Krajčovič J. Reductive evolution of chloroplasts in non-photosynthetic plants, algae and protists. Curr Genet. 2018;64:365–87.
Ralph SA, van Dooren GG, Waller RF, Crawford MJ, Fraunholz MJ, Foth BJ, et al. Tropical infectious diseases: metabolic maps and functions of the Plasmodium falciparum apicoplast. Nature Rev Microbiol. 2004;2:203–16.
Röhrich RC, Englert N, Troschke K, Reichenberg A, Hintz M, Seeber F, et al. Reconstitution of an apicoplast-localised electron transfer pathway involved in the isoprenoid biosynthesis of Plasmodium falciparum. FEBS Lett. 2005;579:6433–8.
Balconi E, Pennati A, Crobu D, Pandini V, Cerutti R, Zanetti G, et al. The ferredoxin-NADP+ reductase/ferredoxin electron transfer system of Plasmodium falciparum. FEBS J. 2009;276:3825–36.
Pombert JF, Blouin NA, Lane C, Boucias D, Keeling PJ. A lack of parasitic reduction in the obligate parasitic green alga Helicosporidium. PLoS Genet. 2014;10:e1004355.
Smith DR, Lee RW. A plastid without a genome: evidence from the nonphotosynthetic green algal genus Polytomella. Plant Physiol. 2014;164:1812–9.
Asmail SR, Smith DR. Retention, erosion, and loss of the carotenoid biosynthetic pathway in the nonphotosynthetic green algal genus Polytomella. New Phytol. 2016;209:899–903.
Ravin NV, Gruzdev EV, Beletsky AV, Mazur AM, Prokhortchouk EB, Filyushin MA, et al. The loss of photosynthetic pathways in the plastid and nuclear genomes of the non-photosynthetic mycoheterotrophic eudicot Monotropa hypopitys. BMC Plant Biol. 2016;16:238.
Janouškovec J, Gavelis GS, Burki F, Dinh D, Bachvaroff TR, Gornik SG, et al. Major transitions in dinoflagellate evolution unveiled by phylotranscriptomics. Proc Natl Acad Sci U S A. 2017;114:E171–80.
Kamikawa R, Moog D, Zauner S, Tanifuji G, Ishida KI, Miyashita H, et al. A non-photosynthetic diatom reveals early steps of reductive evolution in plastids. Mol Biol Evol. 2017;34:2355–66.
Dorrell RG, Azuma T, Nomura M, Audren de Kerdrel G, Paoli L, Yang S, et al. Principles of plastid reductive evolution illuminated by nonphotosynthetic chrysophytes. Proc Natl Acad Sci U S A. 2019;116:6914–23.
Gawryluk RMR, Tikhonenkov DV, Hehenberger E, Husnik F, Mylnikov AP, Keeling PJ. Non-photosynthetic predators are sister to red algae. Nature. 2019;572:240–3.
Janouškovec J, Paskerova GG, Miroliubova TS, Mikhailov KV, Birley T, Aleoshin VV, et al. Apicomplexan-like parasites are polyphyletic and widely but selectively dependent on cryptic plastid organelles. eLife. 2019;8:e49662.
Maciszewski K, Karnkowska A. Should I stay or should I go? Retention and loss of components in vestigial endosymbiotic organelles. Curr Opin Genet Dev. 2019;58–59:33–9.
Sibbald SJ, Archibald JM. Genomic insights into plastid evolution. Genome Biol Evol. 2020;12:evaa096.
Carol P, Stevenson D, Bisanz C, Breitenbach J, Sandmann G, Mache R, et al. Mutations in the Arabidopsis gene IMMUTANS cause a variegated phenotype by inactivating a chloroplast terminal oxidase associated with phytoene desaturation. Plant Cell. 1999;11:57–68.
Shikanai T. Cyclic electron transport around photosystem I: genetic approaches. Annu Rev Plant Biol. 2007;58:199–217.
Nawrocki WJ, Tourasse NJ, Taly A, Rappaport F, Wollman FA. The plastid terminal oxidase: its elusive function points to multiple contributions to plastid physiology. Annu Rev Plant Biol. 2015;66:49–74.
Peltier G, Aro EM, Shikanai T. NDH-1 and NDH-2 plastoquinone reductases in oxygenic photosynthesis. Annu Rev Plant Biol. 2015;67:55–80.
Peltier G, Cournac L. Chlororespiration. Annu Rev Plant Biol. 2002;53:523–50.
Kreimer G. The green algal eyespot apparatus: a primordial visual system and more? Curr Genet. 2009;55:19–43.
Figueroa-Martinez F, Nedelcu AM, Smith DR, Adrian RP. When the lights go out: the evolutionary fate of free-living colorless green algae. New Phytol. 2015;206:972–82.
Nakada T, Misawa K, Nozaki H. Molecular systematics of Volvocales (Chlorophyceae, Chlorophyta) based on exhaustive 18S rRNA phylogenetic analyses. Mol Phylogenet Evol. 2008;48:281–91.
Figueroa-Martinez F, Nedelcu AM, Smith DR, Reyes-Prieto A. The plastid genome of Polytoma uvella is the largest known among colorless algae and plants and reflects contrasting evolutionary paths to nonphotosynthetic lifestyles. Plant Physiol. 2017;173:932–43.
Stern D. The Chlamydomonas sourcebook: organellar and metabolic processes. Cambridge, Massachusetts: Academic Press; 2009.
Alboresi A, Dall’Osto L, Aprile A, Carillo P, Roncaglia E, Cattivelli L, Bassi R. Reactive oxygen species and transcript analysis upon excess light treatment in wild-type Arabidopsis thaliana vs a photosensitive mutant lacking zeaxanthin and lutein. BMC Plant Biol. 2011;11:62.
Tonhosolo R, D’Alexandri FL, de Rosso VV, Gazarini ML, Matsumura MY, Peres VJ, et al. Carotenoid biosynthesis in intraerythrocytic stages of Plasmodium falciparum. J Biol Chem. 2008;284:9974–85.
Cheng Z, Sattler S, Maeda H, Sakuragi Y, Bryant DA, DellaPenna D. Highly divergent methyltransferases catalyze a conserved reaction in tocopherol and plastoquinone synthesis in cyanobacteria and photosynthetic eukaryotes. Plant Cell. 2003;15:2343–56.
Puthiyaveetil S, Ibrahim IM, Allen JF. Oxidation-reduction signalling components in regulatory pathways of state transitions and photosystem stoichiometry adjustment in chloroplasts. Plant Cell Environ. 2012;35:347–59.
Razquin P, Fillat MF, Schmitz S, Stricker O, Böhme H, Gómez-Moreno C, et al. Expression of ferredoxin–NADP+ reductase in heterocysts from Anabaena sp. Biochem J. 1996;316:157–60.
Lange H, Kaut A, Kispal G, Lill R. A mitochondrial ferredoxin is essential for biogenesis of cellular iron-sulfur proteins. Proc Natl Acad Sci U S A. 2000;97:1050–5.
Mulo P, Medina M. Interaction and electron transfer between ferredoxin–NADP+ oxidoreductase and its partners: structural, functional, and physiological implications. Photosynth Res. 2017;134:265–80.
Karlusich JJP, Carrillo N. Evolution of the acceptor side of photosystem I: ferredoxin, flavodoxin, and ferredoxin-NADP+ oxidoreductase. Photosynth Res. 2017;134:235–50.
Seeber F, Soldati-Favre D. Metabolic pathways in the apicoplast of apicomplexa. Int Rev Cell Mol Biol. 2010;281:161–228.
Yoshida K, Hisabori T. Distinct electron transfer from ferredoxin–thioredoxin reductase to multiple thioredoxin isoforms in chloroplasts. Biochem J. 2017;474:1347–60.
La Roche J, Boyd PW, McKay RML, Geider RJ. Flavodoxin as an in situ marker for iron stress in phytoplankton. Nature. 1996;382:802–5.
NIES Microbial Culture Collection Media List. https://mcc.nies.go.jp/02medium-e.html Accessed 10 Jul 2020.
NIES Microbial Culture Collection. https://mcc.nies.go.jp/ Accessed 8 Aug 2020.
Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2003;30:3059–66.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.
Lartillot N, Philippe H. A Bayesian mixture model for across-site heterogeneities in the amino-acid replacement process. Mol Biol Evol. 2004;21:1095–109.
Nguyen LT, Schmidt HA, von Haeseler A, Minh BQ. IQ-tree: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32:268–74.
fastX toolkit. http://hannonlab.cshl.edu/fastx_toolkit/ Accessed 10 Jul 2020.
Chikhi R, Medvedev P. Informed and automated k-mer size selection for genome assembly. Bioinformatics. 2014;30:31–7.
Chin CS, Alexander DH, Marks P, Klammer AA, Drake J, Heiner C, et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods. 2013;10:563–9.
Nagasaki H, Mochizuki T, Kodama Y, Saruhashi S, Morizaki S, Sugawara H, et al. DDBJ read annotation pipeline: a cloud computing-based pipeline for high-throughput analysis of next-generation sequencing data. DNA Res. 2013;20:383–90.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
Mfannot. http://megasun.bch.umontreal.ca/cgi-bin/mfannot/mfannotInterface.pl Accessed 10 Jul 2020.
Lowe TM, Chan PP. tRNAscan-SE on-line: integrating search and context for analysis of transfer RNA genes. Nucl Acids Res. 2016;44:W54–7.
Chlamydomonas reinhardtii strain CC-277 cw15 mt+ mitochondrion, complete genome. https://www.ncbi.nlm.nih.gov/nuccore/EU306622 Accessed 13 Aug 2020.
Chlamydomonas leiostraca strain SAG 11–49 mitochondrion, complete genome. https://www.ncbi.nlm.nih.gov/nuccore/KP696389 Accessed 13 Aug 2020.
Emanuelsson O, Nielsen H, von Heijne G. ChloroP, a neural network-based method for predicting chloroplast transit peptides and their cleavage sites. Protein Sci. 1999;8:978–84.
Fukasawa Y, Tsuji J, Fu SC, Tomii K, Horton P, Imai K. MitoFates: improved prediction of mitochondrial targeting sequences and their cleavage sites. Mol Cell Proteomics. 2015;14:1113–26.
Nishimura O, Hara Y, Kuraku S. gVolante for standardizing completeness assessment of genome and transcriptome assemblies. Bioinformatics. 2017;33:3635–7.
Zapata M, Rodríguez F, Garrido JL. Separation of chlorophylls and carotenoids from marine phytoplankton: a new HPLC method using a reversed phase C8 column and pyridine-containing mobile phases. Mar Ecol Prog Ser. 2000;195:29–45.
Takaichi S. Characterization of carotenes in a combination of a C18 HPLC column with isocratic elution and absorption spectra with a photodiode-array detector. Photosyn Res. 2000;65:93–9.
Naito Y, Yoshimura J, Morishita S, Ui-Tei K. siDirect 2.0: updated software for designing functional siRNA with reduced seed-dependent off-target effect. BMC Bioinformatics. 2009;10:392.
Yamano T, Iguchi H, Fukuzawa H. Rapid transformation of Chlamydomonas reinhardtii without cell-wall removal. J Biosci Bioeng. 2013;115:691–4.
Genome sequencing of a non-photosynthetic green alga. https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJDB10134 Accessed 8 Aug 2020.
Transcriptome of a non-photosynthetic green alga. https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJDB9052 Accessed 8 Aug 2020.
Volvocales sp. NrCl902 leucoplast DNA, complete genome. https://www.ncbi.nlm.nih.gov/nuccore/LC516060 Accessed 13 Aug 2020.
Volvocales sp. NrCl902 mitochondrial DNA, complete genome. https://www.ncbi.nlm.nih.gov/nuccore/LC516061 Accessed 13 Aug 2020.
Computations were partially performed on the NIG supercomputer at ROIS National Institute of Genetics.
This work was in part supported by JSPS grants awarded to RK (no. 19H03274) and to YK (no. 18H03743). RK was supported by The Yanmar Environmental Sustainability Support Association (KI019045).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Plastid genome of chlamydomonad sp. NrCl902. Dark gray boxes show canonical plastid genes. Closed boxes show genes with no homolog in other organisms and intronic ORFs. Light gray boxes between dark gray ones show introns. Thin bars show tRNA genes.
Predicted plastid metabolic map of chlamydomonad sp. NrCl902. 1. Heme synthesis, 2. branched chain amino acid synthesis, 3. sulfate assimilation, 4. Fe-S cluster synthesis, 5. aspartate-to-lysine conversion, 6. starch metabolism, 7. pentose phosphate pathway, 8. nitrite assimilation, 9. aromatic amino acid biosynthesis, 10. fatty acid synthesis, and 11. glycerolipid synthesis. Sequences in which plastid-targeting transit peptides were detected are shown by green circles with solid lines, while those lacking explicit plastid-targeting transit peptides are shown by green circle with no line. Details for isoprenoid (IPP), carotenoid, and plastoquinone syntheses are depicted in Fig. 3a. Substrate possibly imported from the cytosol and/or mitochondria are highlighted in blue, while those possibly exported to contribute to the cytosolic and mitochondrial functions are in red. Abbreviations of proteins are explained in Additional file 3 (Table S1).
Plastid-targeted protein sequences detected in the transcriptome data of chlamydomonad sp. NrCl902.
Interaction between the non-photosynthetic plastid and other compartments. Protein sequences with clear plastid-targeting sequences are shown by light green circles enclosed by solid lines. Protein sequences with clear mitochondrial targeting sequences are shown by orange circles enclosed by red lines. Protein sequences with neither targeting sequence are shown by light blue circles. “?”: Aconitase that catalyzes the conversion from citrate to isocitrate outside mitochondria was not detected in the transcriptome data of chlamydomonad sp. NrCl902. Translation initiation from the 2nd methionine of the mitochondrial aconitase gene might express Aconitase functioning outside mitochondria. Abbreviations are explained in Tables S1 and S2.
Mitochondrial genome of chlamydomonad sp. NrCl902. Conserved mitochondrial genes are shown by closed boxes, while intronic open reading frames are shown in gray. Intron regions are shown as open boxes. L1-L6 and S1-S3 show large subunit and small subunit rRNA gene fragments. Transfer RNA genes are shown by their amino acids and anticodons in parentheses. Given a variety of structures of mitochondrial genomes in Volvocales, i.e., circular genomes and tandem repeats of linear genome, it remains unclear whether the mitochondrial genome is a circularly mapping molecule or a linear, tandemly repeated molecule.
Mitochondrial and cytosolic protein sequences for carbon and energy metabolisms detected in the transcriptome data of chlamydomonad sp. NrCl902.
Carotenoid biosynthesis in chlamydomonad sp. NrCl902. A. The detailed pathway for carotenoid biosynthesis and structures of carotenoids predicted to be synthesized in this pathway. B. HPLC profile for carotenoid detection by absorbance at 450 nm in chlamydomonad sp. NrCl902.
Detected sequences for electron transport systems in non-photosynthetic species of algal and plant lineages.
Detected sequences for representative plastid metabolic functions in the non-photosynthetic eudicot plant Monotropa hypopytis.
Maximum likelihood tree of plastid terminal oxidase in eukaryotes. Non-photosynthetic algae are highlighted in black. Numbers on branches are bootstrap values equal to or higher than 80%. The dataset comprised of 67 taxa and 236 sites was analyzed with IQtree under the LG + I + Γ model selected with Bayesian Information Criterion.
Maximum likelihood tree of Ferredoxin:NADP+ oxidoreductase in eukaryotes. Non-photosynthetic algae are highlighted in black. Numbers on branches are bootstrap values equal to or higher than 80%. The dataset comprised of 96 taxa and 303 sites was analyzed with IQtree under the LG + I + Γ model selected with Bayesian Information Criterion.
Maximum likelihood tree of MPBQ/MSBQ methyltransferase in eukaryotes. A. Plant-type MMT. The dataset comprised of 17 taxa and 250 sites was analyzed with IQtree under the LG + Γ model selected with Bayesian Information Criterion. B. Divergent type MMT. The dataset comprised of 35 taxa and 292 sites was analyzed with IQtree under the LG + I + Γ model selected with Bayesian Information Criterion. Non-photosynthetic algae are highlighted in black. Numbers on branches are bootstrap values equal to or higher than 80%.
Datasets for phylogenetic analyses.
Parameters of the multiple reaction monitoring (MRM) in the positive ion mode of the mass spectrometer used in the present study.
Parameters of the Q3 product ion scan in the positive ion mode of the mass spectrometer used in the present study.
About this article
Cite this article
Kayama, M., Chen, JF., Nakada, T. et al. A non-photosynthetic green alga illuminates the reductive evolution of plastid electron transport systems. BMC Biol 18, 126 (2020). https://doi.org/10.1186/s12915-020-00853-w
- Non-photosynthetic plastid
- Plastid terminal oxidase
- Transient RNA interference knockdown