The role of photorespiration during the evolution of C4 photosynthesis in the genus Flaveria

Abstract
eLife digest
Introduction
Results
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

C₄ photosynthesis represents a most remarkable case of convergent evolution of a complex trait, which includes the reprogramming of the expression patterns of thousands of genes. Anatomical, physiological, and phylogenetic and analyses as well as computational modeling indicate that the establishment of a photorespiratory carbon pump (termed C₂ photosynthesis) is a prerequisite for the evolution of C₄. However, a mechanistic model explaining the tight connection between the evolution of C₄ and C₂ photosynthesis is currently lacking. Here we address this question through comparative transcriptomic and biochemical analyses of closely related C₃, C₃–C₄, and C₄ species, combined with Flux Balance Analysis constrained through a mechanistic model of carbon fixation. We show that C₂ photosynthesis creates a misbalance in nitrogen metabolism between bundle sheath and mesophyll cells. Rebalancing nitrogen metabolism requires anaplerotic reactions that resemble at least parts of a basic C₄ cycle. Our findings thus show how C₂ photosynthesis represents a pre-adaptation for the C₄ system, where the evolution of the C₂ system establishes important C₄ components as a side effect.

https://doi.org/10.7554/eLife.02478.001

eLife digest

Environmental pressures sometimes cause different organisms to independently evolve the same traits. A dramatic example of this phenomenon, which is called convergent evolution, can be seen in the modes used by plants to convert carbon dioxide from the air into starch during photosynthesis.

Early plants existed in an environment with high levels of carbon dioxide in the air. Over time, carbon dioxide levels decreased, so plants evolved more efficient types of photosynthesis to cope. A very efficient type of photosynthesis, called C₄ photosynthesis essentially represents a carbon dioxide concentration mechanism. It has evolved at least 62 times independently in 19 different families of flowering plants.

Scientists have shown that a less advanced, low-efficiency version of photosynthetic carbon dioxide concentration, called C₂ photosynthesis, is a stepping-stone to C₄ photosynthesis. It is also known that the evolution of C₄ photosynthesis required changes to the expression patterns of thousands of genes, but the exact mechanism that leads from C₂ photosynthesis to C₄ photosynthesis is not clear.

To explore this in greater detail, Mallmann, Heckmann et al. studied plants from the genus Flaveria, which belongs to the same family as sunflowers and asters. Under identical greenhouse conditions, plants that use three different photosynthetic pathways—C₃ photosynthesis, C₄ photosynthesis, or an intermediate between the two—were grown and their gene expression patterns were compared. Computer simulations were used to model the metabolism of plants that relied on C₂ photosynthesis.

Based on the modeling, it appears that C₂ photosynthesis shifts the balance of nitrogen metabolism between two types of cell that are critical to photosynthesis. To rebalance the nitrogen, several genes are expressed to trigger an ammonia recycling mechanism. The same genes are turned on during C₄ photosynthesis, and this recycling mechanism include parts of the C₄ process.

The findings of Mallmann, Heckmann et al. suggest that the initial steps in C₄ photosynthesis evolved to prevent nitrogen imbalance. Over time, this mechanism was co-opted to become part of a more efficient form of photosynthesis, which may explain why so many different plants evolved from C₂ to C₄ photosynthesis.

https://doi.org/10.7554/eLife.02478.002

Introduction

The dual-specific enzyme ribulose 1,5-bisphosphate carboxylase/oxygenase (Rubisco) catalyzes two opposing reactions—the carboxylation and the oxygenation of ribulose 1,5-bisphosphate. The former reaction yields 3-phosphoglycerate (3-PGA), whereas the latter produces 2-phosphoglycolate (2-PG). 3-PGA is reduced to carbohydrates in the Calvin–Benson cycle and incorporated into biomass. However, 2-PG is toxic, which requires its removal by a metabolic repair pathway called photorespiration (Anderson, 1971; Bowes et al., 1971; Ogren, 1984; Leegood et al., 1995). In the photorespiratory cycle, 2-PG is regenerated to 3-PGA, but it involves the release of formerly assimilated CO₂ and NH₃, entails energy costs for the plants and reduces the efficiency of photosynthesis by up to 30% (Ehleringer et al., 1991; Bauwe et al., 2010; Raines, 2011; Fernie et al., 2013). Eight core enzymes are required for photorespiration, which in higher plants are located in the chloroplast, the peroxisome, and the mitochondrion (Bauwe et al., 2010; Figure 1A). The pathway rescues ¾ of the carbon, which would otherwise be lost through the oxygenase activity of Rubisco (Peterhansel et al., 2010; Fernie et al., 2013). Ammonia refixation in the chloroplast by the combined activities of glutamine synthase (GS) and glutamine oxoglutarate aminotransferase (GOGAT) is an integral part of photorespiration.

Figure 1

Download asset Open asset

The genus *Flaveria* as a model organism to study C₄ evolution.

Schematic view of the photorespiratory pathway (A), the NADP-ME type C₄ pathway as it can be found in C₄ Flaveria species (B) and the C₂ photosynthesis pathway (C). (D) Phylogeny and physiological properties of selected Flaveria species. The phylogeny was redrawn according to McKown et al. (2005), CO₂ compensation points are taken from Ku et al. (1991), incorporation of ¹⁴CO₂ is from Moore et al. (1987) and the ratios of GLDP B (expressed in all chlorenchyma cells) and GLDP A (expressed in bundle sheath cells only) are from Schulze et al. (2013). (Abbreviations: AGT: serine glyoxylate aminotransferase; AlaAT: alanine aminotransferase; AspAT: aspartate aminotransferase; GDC: glycine decarboxylase complex; GGT: glutamate, glyoxylate-aminotransferase; GLYK: D-glycerate 3-kinase; GOX: glycolate oxidase; HPR: hydroxypyruvate reductase; MDH: malate dehydrogenase; NADP-ME: NADP dependent malic enzyme; PEPC: phosphoenolpyruvate carboxylase; PGLP: 2-phosphoglycerate phosphatase; PPDK pyruvate, phosphate-dikinase; RUBISCO: Ribulose-1,5-bisphosphat-carboxylase/-oxygenase; SHM: serine hydroxymethyltransferase; 2-OG: oxoglutarate; 2-PG 2-phosphoglycolate; 3-PGA: 3-phosphoglycerate; Gln: glutamine; Glu: glutamate; OAA: oxaloacetate; PEP: phosphoenolpyruvate; TP: triosephosphate).

https://doi.org/10.7554/eLife.02478.003

In hot and dry environments and under low atmospheric CO₂ conditions, when the oxygenation activity of Rubisco is increased, the high rate of photorespiration becomes unfavorable for the plants (Sage, 2001, 2013). C₄ plants possess a mechanism that minimizes the oxygenase function of Rubisco and thereby reduces photorespiration and decreases the loss of carbon. C₄ photosynthesis is based on a division of labor between two different cell types, mesophyll and bundle sheath cells, which are organized in a wreath-like structure called ‘Kranz Anatomy’ (Haberlandt, 1904; Dengler and Nelson, 1999). Atmospheric CO₂ is initially fixed in the mesophyll by phosphoenolpyruvate carboxylase (PEPC), and the resulting four-carbon compound is transported to the bundle sheath cells and decarboxylated by NADP/NAD malic enzyme or phosphoenolpyruvate carboxykinase (Hatch et al., 1975). Thereby CO₂ is concentrated at the site of the Rubisco in the bundle sheath cells (Hatch, 1987), outcompeting the molecular oxygen. As a consequence, photorespiration is drastically reduced as compared to C₃ plants, and C₄ plants are characterized by a high photosynthetic efficiency (Figure 1B).

C₄ plants have evolved multiple times independently from C₃ ancestors. The evolution of C₄ photosynthesis occurred at least 62 times in 19 different families of the angiosperms (Sage et al., 2011), implying a low evolutionary barrier towards expression of this trait. The analysis of recent intermediate species (Bauwe and Kolukisaoglu, 2003; Sage, 2004; Bauwe, 2011; Sage et al., 2012, 2013; Schulze et al., 2013) indicates that establishing a photorespiratory CO₂ pump was an early and important step in the evolution towards C₄ photosynthesis (Figure 1C). Since the two-carbon compound glycine serves as a transport metabolite, this photorespiratory CO₂ concentrating mechanism is also termed C₂ photosynthesis. Computational modeling of the evolutionary trajectory from C₃ to C₄ photosynthesis indicated C₂ photosynthesis represented an evolutionary intermediate state (Heckmann et al., 2013; Williams et al., 2013) as well suggesting that C₂ photosynthesis is a prerequisite for the evolution of C_4. However, it remained unclear if the evolution of C₂ photosynthesis fosters the evolution of C₄ photosynthesis beyond providing a selection pressure to reallocate Rubisco to the bundle sheath.

In the present study, we have used the genus Flaveria as a model system for investigating the transition from C₂ to C₄ photosynthesis. To this end, we study a phylogenetic framework consisting of C₃, C₃–C₄ intermediate, and C₄ species (Powell, 1978; Edwards and Ku, 1987; Ku et al., 1991) of this genus which rather recently evolved C₄ (Christin et al., 2011), focusing on genes encoding photorespiratory enzymes and other components of C₂ photosynthesis. The genus Flaveria contains three main phylogenetic groups, of which the first diverging group includes all C₃ Flaveria. Clade B contains seven C₃–C₄ intermediate species and the C₄-like species F. brownii. All C₄ Flaveria species belong to clade A, which also contains several C₄-like species and the C₃–C₄ intermediate F. ramosissima (McKown et al., 2005; Figure 1D). We hypothesized that the analysis of species in the genus Flaveria combined with in silico modeling elucidates the evolutionary changes accompanying and following the establishment of the C₂ pathway. To this end we simulated the metabolism of C₂ plants by coupling a mechanistic model of C₃–C₄ intermediate photosynthesis (von Caemmerer, 2000; Heckmann et al., 2013) with a detailed modified stoichiometric model of C₄ photosynthesis (Dal'Molin et al., 2010), and investigated the evolution of C₄ photosynthesis and photorespiration by following the changes in mRNA and protein abundance along the evolutionary path.

RNA and protein amounts of the majority of the photorespiratory enzymes were reduced in C₄ as compared to C₃ species. In contrast, photorespiratory mRNA and protein amounts did not decrease in the C₃–C₄ intermediate species but were mostly equal or even higher than in the C₃ species, demonstrating that the establishment of the photorespiratory CO₂ pump in the genus Flaveria relies on coordinated changes in the expression of all core photorespiratory enzymes. Metabolic modeling in combination with comparisons of transcript abundances in the different Flaveria species strongly indicates that introduction of C₂ photosynthesis has a direct impact on the nitrogen metabolism of the leaf. Its implementation necessitates the parallel establishment of components of the C₄ cycle to cope with these changes in refixation of photorespiratory nitrogen. Based on these results, we predict a mechanistic interaction between C₄ and C₂ photosynthesis.

Results

Selection of Flaveria species, cultivation of plant material and experimental design

To study the evolution of the expression of photorespiratory and C₄ cycle genes during the transition from C₃ to C₄ photosynthesis in the genus Flaveria, nine species reflecting the evolutionary trajectory taken were selected, including two C₃ (F. robusta and F. pringlei), two C₄ (F. bidentis and F. trinervia), and five C₃–C₄ intermediate species (Figure 1D). According to their CO₂ compensation points and the percentage of carbon initially fixed into malate and aspartate, F. chloraefolia and F. pubescens were earlier classified as type I C₃–C₄ intermediates. F. anomala and F. ramosissima belong to the type II C₃–C₄ intermediates and F. brownii is classified as a C₄-like species (Edwards and Ku, 1987; Moore et al., 1987; Cheng et al., 1988; Ku et al., 1991). Type I C₃–C₄ intermediates are defined as solely relying on the photorespiratory CO₂ concentration cycle whereas a basal C₄ cycle activity is present in type II C₃–C₄ intermediates species. C₄-like species exhibit much higher C₄ cycle activities but lack complete bundle sheath compartmentation of Rubisco activity (Edwards and Ku, 1987).

Four independent experiments with plants grown during different seasons were performed to identify differences between the species that are dependent on their different modes of photosynthesis and independent of environmental influences. For each experiment the plants were seeded concurrently and grown side-by-side under greenhouse conditions. The second and fourth visible leaves from the top of all nine species were harvested at noon on the same day for transcript and protein analysis. Plants for experiment one were harvested in September 2009, for experiment two in June 2010, for experiment three in October 2010 and for experiment four in April 2011. The amounts of the core photorespiratory and C₄ enzymes were assessed by immunoblotting using specific antibodies raised against synthetic peptides or recombinant proteins. The abundances of the corresponding RNAs as well of C₄ cycle associated transcripts were quantified by total transcriptome sequencing.

The transcript profiles of the individual Flaveria species were comparable throughout all four experiments

The transcriptomes of the different Flaveria species were sequenced by Illumina technology following standard procedures. In total, close to 200 Gb of raw sequence data were produced. After filtering of low quality reads 30 to 58 million reads per species and experiment were quantified (Figure 2—source data 1). In a cross species approach, we mapped the sequences onto the minimal set of Arabidopsis thaliana coding sequences using the BLAST-like alignment tool BLAT (Kent, 2002) as described previously (Gowik et al., 2011) (Figure 2—source data 2, data available from the Dryad Digital Repository: http://dx.doi.org/10.5061/dryad.q827h). We were able to align approx. 50% of our reads to the Arabidopsis transcripts. This is lower as compared to a similar approach using 454 sequencing (Gowik et al., 2011) and likely due to the shorter read length of the Illumina compared to the 454 reads. To overcome the low mapping efficiency, the leaf transcriptomes of Flaveria species were assembled de novo based on 454 (Gowik et al., 2011) and Illumina reads (this study). Among the contigs from F. robusta, we identified full-length transcripts for all photorespiratory and C₄ genes in the focus of the present study and used these for further read mapping and detailed analysis.

To evaluate the variation between the four independent experiments, we performed hierarchical sample clustering and a principal component analysis of the transcript profiles derived from read mapping on the minimal set of Arabidopsis coding sequences. Hierarchical sample clustering using Pearson correlation and average linkage clustering shows that the transcript profiles of all Flaveria species were quite similar in all four experiments since the samples cluster strictly species-wise (Figure 2A). The transcriptome patterns are influenced by the photosynthesis type and the phylogenetic relationships of the different species. The two C₄ species, both belonging to clade A, cluster together as do the two C₃ species that belong to the basal Flaveria species. Within the C₃–C₄ intermediates the two more advanced intermediates F. ramosissima and F. anomala cluster together, the only pattern which contradicts phylogenetic proximity since F. ramosissima belongs to clade A and F. anomala belongs to clade B. The last cluster consists out of the C₃–C₄ intermediates F. chloraefolia and F. pubescens, and the C₄-like species F. brownii.

Figure 2

Download asset Open asset

Variation of transcript profiles of the individual *Flaveria* species between the four experiments.

(A) Hierarchical sample clustering of all expressed transcripts. The tree was calculated with the MEV program using the HCL module with Pearson correlation and the average linkage method. (B) Principal component analysis of transcript levels. The first three components explain 27% of the total variance.

https://doi.org/10.7554/eLife.02478.004

Figure 2—source data 1 Results of the Illumina sequencing and cross species read mapping.: https://doi.org/10.7554/eLife.02478.005
Download elife-02478-fig2-data1-v2.xlsx
Figure 2—source data 2 Quantitative information for all reads mapped in a cross species approach onto the reference transcriptome from Arabidopsis thaliana.: https://doi.org/10.7554/eLife.02478.006
Download elife-02478-fig2-data2-v2.xlsx

Principle component analysis supports the results of the hierarchical clustering. The samples are mainly separated by photosynthesis type and phylogenetic relationships with the two intermediate species from different phylogenetic trajectories again forming a tight cluster (Figure 2B). The first three components, shown in Figure 2B, explain only 27% of the total variance. This is in good accordance with earlier results where it was shown that about 16% of all analyzed genes showed photosynthesis type related expression changes when the transcriptomes of the C₄ species F. trinervia and F. bidentis, the C₃ species F. robusta and F. pringlei and the C₃–C₄ intermediate species F. ramosissima were compared (Gowik et al., 2011).

Amounts of photorespiratory transcripts and proteins indicate that the C₂ pathway was established early during C₄ evolution in Flaveria and is present also in the C₄-like species F. brownii

Photorespiratory genes are expressed in all species and photorespiratory proteins are detected in all species. To visualize the differences in transcript and protein abundance heat maps were plotted (Figure 3). The transcription of all photorespiratory genes except the transport proteins DIT1 and DIT2 and one isoform of GLDH was downregulated in the C₄ species F. bidentis and F. trinervia compared to the C₃ species F. pringlei und F. robusta (Figure 3A, Figure 3—source data 1). Both dicarboxylate transporters play an important role in generating the transfer acids in the C₄ pathway of NADP-ME plants such as F. trinervia and F. bidentis (Renne et al., 2003; Gowik et al., 2011; Kinoshita et al., 2011). This may explain why their expression pattern is more similar to the C₄ genes than to the other photorespiratory genes.

Figure 3 with 2 supplements see all

Download asset Open asset

Abundance of photorespiratory transcripts and proteins in leaves of individual *Flaveria* species.

Normalized transcript (A) and protein (B) amounts are plotted as heat maps. Transcript amounts were determined by Illumina sequencing of the leaf transcriptomes and read mapping on selected *F. robusta* full length transcript sequences. Protein amounts were determined by protein gel blots. See Figure 3—source data 1 for absolute transcript levels, Figure 3—source data 2 for protein quantification and Figure 3—figure supplements 1 and 2 for immunoblots. Fp: *F. pringlei* (C₃); Fro: *F. robusta* (C₃); Fc: *F. chloraefolia* (C₃–C₄); Fpu: *F. pubescens* (C₃–C₄); Fa: *F. anomala* (C₃–C₄); Fra: *F. ramosissima* (C₃–C₄); Fbr: *F. brownii* (C₄-like); Fb: *F. bidentis* (C₄); Ft: *F. trinervia* (C₄).

https://doi.org/10.7554/eLife.02478.007

Figure 3—source data 1 Transcript abundance of photorespiratory genes determined by read mapping on F. robusta full length transcript sequences.: https://doi.org/10.7554/eLife.02478.008
Download elife-02478-fig3-data1-v2.xlsx
Figure 3—source data 2 Quantification of photorespiratory proteins by protein gel blot.: https://doi.org/10.7554/eLife.02478.009
Download elife-02478-fig3-data2-v2.xlsx

The amounts of photorespiratory transcripts did not decrease gradually from C₃ to C₄ but the expression levels in the C₃–C₄ intermediate species F. chloraefolia, F. pubescens, F. anomala and F. ramosissima were mostly equal or higher than in the C₃ species. An exception are the transcripts of one GLDP, one GLDH and one SHM isoform which are drastically down-regulated also in the C₃–C₄ intermediate species. It was shown earlier that the down-regulation of this GLDP isoform is tightly associated with the establishment of the C₂ pathway in Flaveria (Schulze et al., 2013). The down-regulation of the GLDH and SHM isogenes might have similar reasons since both enzymes are also involved in glycine decarboxylation. Only the C₄-like species F. brownii is intermediate with respect to photorespiratory transcripts. 19 of 27 transcripts are reduced compared to the C₃–C₄ intermediate and C₃ species but have higher levels than the true C₄ species F. bidentis and F. trinervia. Exceptions are the components of the glycine decarboxylase complex as the respective transcripts levels are equal to these in the C₃ and C₃–C₄ intermediate species (Figure 3A).

The expression patterns described above were not only found for the genes encoding the core enzymes of photorespiration but also for the genes responsible for recycling of ammonia set free during photorespiration, GS/GOGAT. Also the genes of recently discovered transporters associated with photorespiration, PLGG1 and BOU (Eisenhut et al., 2013; Pick et al., 2013), behave accordingly.

To test whether transcript abundance reflects protein abundance, amounts of core photorespiratory proteins in the leaves of all nine species were quantified by protein gel blots. To this end we generated antibodies against conserved peptides from Flaveria GLDP, GLDT, GLDL, SHM, HPR, PGLP and GLYK proteins. Total leaf proteins were extracted from plant material harvested together with the material used for RNA isolation and equal amounts of protein were separated via SDS gel-electrophoresis prior to blotting (Figure 3—figure supplement 1). The changes of protein amounts essentially reflected the changes of the amounts of the corresponding transcripts (Figure 3B, Figure 3—figure supplement 2, Figure 3—source data 2). The amounts of core photorespiratory proteins in the C₃–C₄ intermediates were equal to the amounts in the C₃ species. A clear reduction of these proteins can be observed only for the true C₄ species and the C₄–like species. F. brownii exhibits intermediate amounts of most photorespiratory proteins. This indicates that the regulation of photorespiratory genes mainly occurs on the transcriptional level and that our approach to analyze the photorespiratory activity by comparative transcriptomics is reasonable.

While the overall patterns remain similar between all independent experiments, individual proteins and transcripts vary between the four experiments. This likely reflects adjustments of photorespiratory gene expression to the different light and temperature conditions in our green house in the different seasons of the year.

We conclude that the four experiments support the establishment of a photorespiratory C₂ cycle early during C₄ evolution in Flaveria and that this C₂ cycle was maintained until Rubisco activity was constricted to the bundle sheath cells in the true C₄ Flaveria species.

An integrated model of C₂ photosynthesis

While the principal physiological differences between C₃ and C₄ leaves are widely understood, knowledge about the metabolic reconfiguration required to implement a functional C₂ pathway into a C₃ leaf is incomplete. In particular, moving glycine from mesophyll to bundle sheath cells (Hylton et al., 1988; Morgan et al., 1993) does not only translocate carbon, it also transports one nitrogen atom per two carbon atoms. Evidently, implementing the C₂ carbon pump requires balancing of metabolic routes to maintain homeostasis of both carbon and nitrogen metabolism (Monson and Rawsthorne, 2000). How this can be achieved is non-intuitive and it thus requires a systematic analysis by metabolic modeling. To this end, we simulated the leaf metabolism of a C₂ plant using an integrated model. We coupled a mechanistic model of C₃–C₄ intermediate photosynthesis (von Caemmerer, 2000; Heckmann et al., 2013) with a modified genome-scale stoichiometric model of C₄ photosynthesis that was designed to describe the entire metabolic interactions of mesophyll and bundle sheath cells in C₄ leaves (Dal'Molin et al., 2010).

We used the mechanistic model to predict constraints for the stoichiometric model. It provided values for net CO₂ uptake, Rubisco carboxylation as well as oxygenation in mesophyll and bundle sheath, CO₂ leakage from the bundle sheath, PEPC activity in the mesophyll, activity of NADP-ME in the bundle sheath, plasmodesmatal flux of glycine and serine, and decarboxylation by the GDC. Given specific activities of the C₂ and C₄ cycles in the mechanistic model, we used flux balance analysis (FBA) to predict detailed flux distributions that follow biologically realistic optimality criteria (Varma and Palsson, 1994). We employed a maximization of leaf biomass production, followed by a minimization of the sum of absolute fluxes including transport processes. In the minimization of total flux, we allocated higher weights to plasmodesmatal fluxes in order to account for the trade-off between CO₂ leakage and diffusion of metabolites between the cells. This framework allows us to investigate the most parsimonious implementation of C₂ and C₄ cycles, given a hypothesis about which metabolites are suitable for plasmodesmatal transport.

The first outcome of simulating the photorespiratory CO₂ pump was that the establishment of the C₂ pathway has indeed a direct impact on the nitrogen metabolism of the leaf. It transports two molecules of glycine from the mesophyll to the bundle sheath, where one molecule each of serine, CO₂, and ammonium are produced. CO₂ is fixed by bundle sheath Rubisco and serine is transferred back to the mesophyll, where it is used for the regeneration of phosphoglycerate and photorespiratory glycine. This results in a net transport of CO₂ but also ammonia from the mesophyll to the bundle sheath. To create a noticeable CO₂ enrichment in the bundle sheath, the C₂ cycle must run with an appreciable capacity; indeed, the mechanistic model of C₃–C₄ intermediate photosynthesis predicted an oxygenation rate of Rubisco of about one third of its carboxylation rate. Running at such rates, the C₂ cycle will create a massive nitrogen imbalance between mesophyll and bundle sheath cells, as was also predicted earlier by Monson and Rawsthorne (2000). Within the stoichiometric model, the free diffusion of ammonia between the two cell types was not allowed, since ammonia is toxic and known to effectively uncouple electrochemical gradients (Krogmann et al., 1959). Thus, ammonia must be refixed in the bundle sheath cells and transferred back to the mesophyll in the form of amino acids. According to the intergrated model, ammonia is fixed by glutamine synthetase and glutamine oxoglutarate aminotransferase (GS/GOGAT) in the bundle sheath cells (Figure 4). Consistent with this prediction, we found that GS/GOGAT transcripts were upregulated in the C₃–C₄ intermediate species (Figure 3).

Figure 4

Download asset Open asset

Flux Balance Analysis of the C₂ photosynthetic pathway.

Predicted fluxes if (A) major amino acids and the corresponding oxoacids and dicarbonic acids are allowed to freely diffuse between cells, (B) the α-ketoglutarate and glutamate transfer between mesophyll and bundle sheath was constrained (C) additionally the transfer of alanine and pyruvate between mesophyll and bundle sheath was constrained (D) transfer of all nitrogen containing compounds except for glycine and serine, which are used by the C₂ cycle were constrained. Fluxes are given in µmol s⁻¹ m⁻². Values in brackets show minimum and maximum of flux resulting from flux variability analysis. Flux of dissolved gasses, sucrose, inorganic compounds and processes that carry flux below 1 µmol s⁻¹ m⁻² are not shown. The sums of absolute fluxes over the plasmodesmata for the different variants were (A): 17.8 µmol s⁻¹ m⁻²; (B): 18.4 µmol s⁻¹ m⁻²; (C): 19.0 µmol s⁻¹ m⁻²; (D): 22.1 µmol s⁻¹ m⁻². See Figure 4—source data 1 for plasmodesmatal fluxes.

https://doi.org/10.7554/eLife.02478.012

Figure 4—source data 1 Fluxes over plasmodesmata depending on the weight on plasmodesmatal fluxes including flux variability analysis.: https://doi.org/10.7554/eLife.02478.013
Download elife-02478-fig4-data1-v2.xlsx

Estimating whether a certain metabolite is suitable for maintaining a diffusional gradient between mesophyll and bundle sheath is an unsolved problem. The impact on regulatory mechanisms and homeostasis of the C₃ leaf may render some metabolites unsuitable to serve as transport metabolites. We address this problem by modeling multiple scenarios that assume different transport metabolites.

If major amino acids and the corresponding oxoacids and dicarbonic acids are allowed to freely diffuse between cells in an integrated model representing a C₂ cycle, glutamate is predicted to be transferred to the mesophyll, where it is deaminated by GGT, regenerating the photorespiratory glycine. The resulting 2-oxoglutarate is transferred back to the bundle sheath cells (Figure 4A). The model preference for glutamate/2-oxoglutarate reflects the minimization of total flux in the FBA model, as this effectively minimizes the number of active enzymatic reactions and holds the plasmodesmatal flux for ammonia balance at one acceptor and one transport metabolite.

To elucidate if alternative solutions exist that contain more steps but retain the same biomass output, the 2-oxoglutarate transfer between mesophyll and bundle sheath was constrained to prevent the glutamate/2-oxoglutarate exchange. The integrated model then predicts an alanine/pyruvate shuttle (Figure 4B). The glutamate produced by GS/GOGAT activity in the bundle sheath cells is used by alanine aminotransferase (Ala-AT) to aminate pyruvate. The resulting alanine is transferred to the mesophyll and trans-aminated by Ala-AT resulting in pyruvate and glutamate. The glutamate is used to regenerate photorespiratory glycine and pyruvate is transferred back to the bundle sheath.

If alanine and pyruvate transfer are also constrained, the model predicts an aspartate/malate shuttle (Figure 4C). This includes the oxidation of malate in the bundle sheath. The resulting oxaloacetate (OAA) is aminated by aspartate aminotransferase (Asp-AT) and aspartate moves to the mesophyll. Here aspartate is trans-aminated by Asp-AT and malate is regenerated by reduction of the resulting OAA and transferred to the bundle sheath.

In all these scenarios, further increasing the weights on plasmodesmatal flux leads to transporter metabolites with increased N carrying capacity such as asparagine (Figure 4—source data 1).

In a restrictive scenario, all nitrogen containing compounds were excluded from plasmodesmatal transport, except for glycine and serine, which are used by the C₂ cycle itself. In this case, the model predicts that bundle sheath derived ammonia is transferred from glutamate to phosphohydroxy-pyruvate by phosphoserine aminotransferase to yield phosphoserine; phosphoserine is then converted to serine by phosphoserine phosphatase. Finally, the serine moves to the mesophyll. This variant includes the transfer of 3-phosphoglycerate from the mesophyll to the bundle sheath, where it is converted to phosphohydroxy pyruvate by 3-phosphoglycerate dehydrogenase (Figure 4D).

The model predicts a mechanistic interaction between C₂ and C₄ cycle

In C₃ plants, basal activities of the typical C₄ cycle enzymes are present (Aubry et al., 2011). When our integrated model is parameterized to include an active C₄ cycle, it predicts that a contingent of the bundle sheath ammonia will be transferred to the mesophyll cells by the C₄ cycle as a biomass neutral alternative to the 2-OG/Glu shuttle or as the unique solution when additional weight on plasmodesmatal fluxes is applied (Figure 5—source data 1). In this solution malate is decarboxylated in the bundle sheath cells. CO₂ is refixed by Rubisco, and the resulting pyruvate is aminated by Ala-AT. Alanine moves to the mesophyll cells, where ammonia is fed into the photorespiratory cycle by Ala-AT and GGT. The resulting pyruvate is converted back to malate by PPDK, PEPC, and NADPH-dependent MDH (Figure 5A). Flux variability analysis shows that only marginal variability in the fluxes of the shuttle is possible (Figure 5—source data 1). According to our model predictions, the cycle is active even at low PEPC activities, such as those measured in C₃ Flaveria species (Gowik et al., 2011; Heckmann et al., 2013). When the C₄ cycle runs with low capacity, according the model, the surplus of bundle sheath ammonia is transferred back to the mesophyll by the glutamate/2-oxoglutarate shuttle. Once the capacity of the C₄ cycle gradually increases, the recirculation of nitrogen is shifted from the glutamate/2-oxoglutarate shuttle towards the C₄ cycle (Figure 5B). The predicted biomass production increases linearly with C₄ cycle activity (Figure 5C). Thus, our model predicts a strong interaction between C₂ and C₄ photosynthesis.

Figure 5

Download asset Open asset

Mechanistic interaction between C₂ and C₄ cycle.

(A) Predicted fluxes when the model is parameterized to include activity of the C₄ cycle enzymes. Fluxes are given in µmol s⁻¹ m⁻². Values in brackets show minimum and maximum of flux resulting from flux variability analysis. The sum of absolute flux over plasmodesmata was 21.9 µmol s⁻¹ m⁻². Flux of dissolved gasses, sucrose, inorganic compounds and processes that carry flux below 1 µmol s⁻¹ m⁻² are not shown. See Figure 5—source data 1 for plasmodesmatal fluxes. (B) Predicted activities of Ala-AT in mesophyll (black line) and bundle sheath (gray line) cells and predicted transfer of α-ketoglutarate from mesophyll to bundle sheath cells (black dashed line) and glutamate from bundle sheath to mesophyll cells (gray dashed line) at low C₄ cycle activities. (C) Changes in biomass production with varying (low) activity of the C₄ cycle in a C₂ plant.

https://doi.org/10.7554/eLife.02478.014

Figure 5—source data 1 Fluxes over plasmodesmata depending on the weight on plasmodesmatal fluxes including flux variability analysis.: https://doi.org/10.7554/eLife.02478.015
Download elife-02478-fig5-data1-v2.xlsx

Analysis of C₄ cycle gene expression in C₃–C₄ intermediate Flaveria species

When the C₂ cycle is running with high capacity, our integrated modeling approach predicts the necessity of auxiliary metabolite fluxes between mesophyll and bundle sheath cells to prevent a massive nitrogen imbalance. Among those auxiliary fluxes were the pyruvate/alanine and the malate/aspartate exchanges. The metabolites used in these shuttles also serve as transport metabolites in C₄ photosynthesis. Furthermore, the model highlights the possibility that a low capacity C₄ cycle balances part of the C₂ cycle ammonia production. Therefore we analyzed in detail the expression of C₄ cycle related genes in our dataset. True C₄ Flaverias, such as F. bidentis or F. trinervia, are believed to use a NADP-ME type C₄ cycle (Moore et al., 1984; Ku et al., 1991; Meister et al., 1996; Gowik et al., 2011). All genes associated with this type of C₄ photosynthesis are gradually upregulated in the analyzed C₃–C₄ intermediate species in line with their degree of ‘C₄-ness’. This is true for the typical C₄ enzymes like PEPC, PPDK, MDH, NADP-ME, Ala-AT and a plastidic aspartate aminotransferase (Asp-AT), as well as for several C₄ associated transporters, such as the pyruvate transporter BASS2, the H⁺/Na⁺ exchanger NHD, the PEP translocator CUE1 and the putative malate and aspartate transporters DIT1 and DIT2 (Weber and von Caemmerer, 2010; Brautigam et al., 2011; Furumoto et al., 2011; Gowik et al., 2011). The regulators of the C₄ enzymes (like PEPC kinase or the PPDK regulatory protein) and enzymes with auxiliary functions of C₄ enzymes (like pyrophosphatases or adenosinmonophosphatases) show a similar pattern (Figure 6, Figure 6—source data 1). To corroborate the results of the transcript abundance measurements, selected C₄ cycle enzymes (PEPC, PPDK, and NADP-ME) were measured by immunoblotting. The protein abundance correlates well with the transcript abundance (Figure 6, Figure 6—source data 2).

Figure 6 with 1 supplement see all

Download asset Open asset

Abundance of C₄ related transcripts and proteins in leaves of individual Flaveria species.

Normalized transcript (A) and protein (B) levels are plotted as heat maps. Transcript amounts were determined by Illumina sequencing of the leaf transcriptomes and read mapping on selected *F. robusta* full length transcript sequences. Protein amounts were determined by protein gel blots. See Figure 6—source data 2 for absolute transcript level, Figure 6—source data 2 for protein quantification and Figure 3—figure supplement 1 for immunoblots. (C) Mean values of transcript levels from all four experiments were clustered by hierarchical using the HCL module of MEV program with Pearson correlation and the average linkage method. The relative transcript abundance for PEPC, PPDK, NADP-ME and Ala-AT (mean values from all four experiments) are plotted for all nine species. Fp: *F. pringlei* (C₃); Fro: *F. robusta* (C₃); Fc: *F. chloraefolia* (C₃–C₄); Fpu: *F. pubescens* (C₃–C₄); Fa: *F. anomala* (C₃–C₄); Fra: *F. ramosissima* (C₃–C₄); Fbr: *F. brownii* (C₄–like); Fb: *F. bidentis* (C₄); Ft: *F. trinervia* (C₄).

https://doi.org/10.7554/eLife.02478.016

Figure 6—source data 1 Transcript abundance of C₄ cycle genes determined by read mapping on F. robusta full length transcript sequences.: https://doi.org/10.7554/eLife.02478.017
Download elife-02478-fig6-data1-v2.xlsx
Figure 6—source data 2 Quantification of C₄ proteins by protein gel blots.: https://doi.org/10.7554/eLife.02478.018
Download elife-02478-fig6-data2-v2.xlsx

The expression changes of C₄ cycle genes do not all follow the same quantitative pattern (Figure 6C). Although all of these genes gradually increase in expression when plants gain C₄ properties, as judged, for example, by the percentage of ¹⁴CO₂ directly fixed into C₄ acids (Vogan and Sage, 2011), the quantitative changes in gene expression are quite different. PEPC and PPDK transcript amounts increase slowly in the C₃–C₄ intermediates F. chloraefolia, F. pubescens, and F. anomala, more steeply in the advanced C₃–C₄ intermediate F. ramosissima and the C₄-like species F. brownii before reaching the highest transcript abundances in the true C₄ species (Figure 6C). In contrast, NADP-ME and Ala-AT gene expression already increase in expression in the more C₃-like intermediate species. Their expression rises more linearly in the further advanced intermediates and plateaus in the C₄-like and C₄ species. If one uses the different Flaveria species as evolutionary proxies as suggested by the results of Heckmann et al. (2013), these results suggest that NADP-ME and Ala-AT are strongly upregulated earlier in evolution than other C₄ core enzymes like PEPC or PPDK.

F. chloraefolia is classified as a type I C₃–C₄ intermediate species, and no enhanced C₄ cycle activity should be present in this species based on the classification. We detected upregulation of all NADP-ME type associated C₄ genes, with some of the genes showing comparable small increases in expression (Figure 6). This is in line with the results of ¹⁴CO₂ uptake studies that indicate about 14% of CO₂ is directly incorporated into C₄ acids in F. chloraefolia, whereas only 6% goes into C₄ acids directly in the C₃ species F. pringlei (Moore et al., 1987). We think therefore that a basal C₄ cycle activity is present in F. chloraefolia and its classification as type I C₃–C₄ intermediate is questionable.

A gene encoding a mitochondrial NAD dependent malate dehydrogenase as well as several cytosolic and especially one mitochondrial Asp-AT were upregulated exclusively in the C₃–C₄ intermediate species and the C₄-like F. brownii (Figure 6). Often, high activities of these genes are associated with the NAD-ME or PEP-CK type of C₄ photosynthesis. NAD dependent malic enzyme and PEP carboxykinase genes were only very lowly expressed in all analyzed Flaverias, and no obvious differences between the C₃–C₄ intermediates and the other species could be found (Figure 2—source data 2, data available from the Dryad Digital Repository: http://dx.doi.org/10.5061/dryad.q827h).

We found no transcriptomic evidence that ammonia is recirculated by the phosphoserine pathway predicted by the model that restricts the free diffusion of all amino acids except serine and glycine. The amounts of transcripts for all three enzymes of this pathway, i.e., phosphoserine aminotransferase, phosphoserine phosphatase, and 3-phosphoglycerate dehydrogenase, were found to be very low in all analyzed Flaveria species (Figure 2—source data 2, data available from the Dryad Digital Repository: http://dx.doi.org/10.5061/dryad.q827h) (Mallmann et al., 2014).

Taken together, these data imply that the anaplerotic ammonia shuttle, required to maintain the nitrogen homeostasis in mesophyll and bundles sheath cells of plants performing C₂ photosynthesis, is active in all analyzed C₃–C₄ Flaveria species, as predicted by the computer simulations. Furthermore, it appears that even the most C₃-like C₃–C₄ intermediate species analyzed within the present study, F. chloraefolia, exhibits low level C₄ cycle activity. This activity is again in accordance with the in silico model, which predicts the C₄ cycle to be a highly efficient ammonia recirculation pathway.

Discussion

Photorespiration is mainly seen as a wasteful process, which arises from a malfunction of Rubisco and reduces photosynthetic efficiency (Ogren, 1984). In a high CO₂ atmosphere, Rubisco can operate efficiently. But the current atmospheric CO₂ concentration, combined with heat and drought, leads to an enhanced oxygenase activity and thereby the photosynthetic efficiency decreases (Raines, 2011). Up to 30% of the initially fixed CO₂ may be lost by photorespiration (Bauwe et al., 2010). C₄ plants avoid this problem by enriching CO₂ at the site of Rubisco. CO₂ is prefixed in the mesophyll and released in the bundle sheath cells, where Rubisco is operating (Hatch, 1987). The establishment of the photorespiratory CO₂ pump, which relocates the release of photorespiratory CO₂ to the bundle sheath cells, appears to be an important intermediate step towards the C₄ cycle and our detailed study of Flaveria intermediate species suggests that genes associated with C₄ photosynthesis also played a role in the C₂ cycle.

Implementation of the C₂ pathway leads to high expression of photorespiratory genes in C₃-C₄ intermediate Flaveria species

The expression of photorespiratory genes, including all genes encoding the core enzymes of the pathway, most of the transporters, and the enzymes involved in ammonia refixation, is not downregulated in the analyzed intermediate species; the transcript and protein amounts remain constant or in some cases are even higher compared to C₃ species. A significant drop in photorespiratory gene expression is only observed in the C₄-like species F. brownii and is decreased further in the C₄ species. Together with earlier results (Schulze et al., 2013), this indicates that indeed a C₂ photosynthetic cycle is active in all these C₃-C₄ intermediate Flaveria species and that a reduction in photorespiratory transcripts and proteins only occurs once the amounts of Rubisco have been reduced in the mesophyll as was described for the C₄-like species F. brownii (Bauwe, 1984; Holaday et al., 1988). Rubisco reduction in the mesophyll is thus a late step of C₄ evolution, which in the Flaveria series appears to not occur gradually but rather abruptly towards the end of the evolutionary trajectory. It is followed by a strong increase of C₄ cycle activity, as can be deduced from the upregulation of PEPC and PPDK genes in the real C₄ species (Figure 6C), when the primary CO₂ fixation is completely taken over by PEPC. In the intermediate species C₂ and C₄ cycles operate in parallel leading to similar or higher photorespiratory gene expression compared with the C₃ species.

Analysis of C₄ cycle gene expression supports the predictions of the C₂ model for C₃–C₄ intermediate Flaveria species and implies the early establishment of a complete C₄ pathway

The model of the C₂ cycle and the underlying metabolism proposes GS/GOGAT, Ala-AT, and Asp-AT to be involved in balancing the amino groups during C₂ cycle operation (Figure 4). The transcriptome data from the C₃-C₄ intermediate Flaveria species largely support the results of our integrated model for the C₂ pathway (Figure 6). In these species we found an upregulation of genes involved in the three most likely mechanisms for the recovery of ammonia predicted by the model. GS/GOGAT, which catalyzes the primary refixation of ammonia in the bundle sheath cells, is important for all three versions of ammonia shuttles (Figure 4) and is upregulated in the intermediate species. Transcripts for the glutamate/2-oxoglutarate shuttle, the alanine/pyruvate shuttle, and the aspartate/malate shuttle are enriched in all C₃–C₄ intermediates compared to the C₃ and C₄ Flaverias. For the alanine/pyruvate shuttle, Ala-AT is needed in the bundle sheath and the mesophyll cells. Ala-AT is upregulated already in the least advanced C₃–C₄ intermediates F. chloraefolia and F. pubescens, but also in all the other C₃–C₄ intermediates. Ala-AT transcripts are also highly abundant in the true C₄ Flaverias since Ala-AT is directly involved in the C₄ cycle when alanine is used as transport metabolite.

We found several Asp-AT and two MDH genes upregulated in the C₃–C₄ intermediate species (Figure 6). The chloroplast-located MDH and Asp-AT genes are involved in the C₄ cycle of C₄ Flaverias, in which malate and aspartate are used concurrently as C₄ transport metabolites (Meister et al., 1996). Two further Asp-AT genes and another MDH gene were found to be upregulated exclusively in the C₃–C₄ intermediates, including the C₄-like species F. brownii. The most likely reason for upregulation of these genes is their involvement in the recirculation of photorespiratory ammonia by a malate/aspartate shuttle.

The pathways of ammonia recirculation between mesophyll and bundle sheath foreshadow the establishment of a true C₄ cycle (Figure 4). All variants described above need the establishment of inter- and intra-cellular transport capacities for amino acids and small organic acids, which are also needed for a functional C₄ cycle (Weber and von Caemmerer, 2010). The existence of an aspartate/malate and an alanine/pyruvate shuttle anticipates important components of a functional C₄ pathway. Our transcript data imply that both of these shuttles are active in C₃–C₄ intermediate Flaverias. Only a few additions would be required to convert these pathways of ammonia recirculation into a C₄-like CO₂ concentration mechanism, that is, malate would have to be decarboxylated in the bundle sheath cells and pyruvate would have to be converted to malate in the mesophyll. Our transcript data implies that this conversion of the photorespiratory ammonia recirculation pathway into a C₄-like CO₂ concentrating pump must have been an early event in C₄ evolution of Flaveria since already in the least advanced intermediates such as F. chloraefolia and F. pubescens, NADP-ME transcripts are elevated and their amounts increase in parallel with Ala-AT and Asp-AT transcript levels.

To extend the pathways of ammonia recirculation into a rudimentary C₄ cycle, a capacity to regenerate malate from pyruvate in the mesophyll is required. As deduced from the transcriptome data, the enzymatic functions required are also already enhanced in the least advanced C₃–C₄ intermediates, since we observe a low but consistent upregulation of PEPC and PPDK genes in these species compared to the C₃ Flaverias. Measurements of radiolabeled CO₂ incorporation support the view that a rudimentary C₄ cycle is already operating in intermediate Flaveria species (Rumpho et al., 1984; Monson et al., 1986; Moore et al., 1987; Chastain and Chollet, 1989). F. chloraefolia as well as F. pubescens incorporate a higher percentage of ¹⁴CO₂ into the C₄ compounds malate and aspartate (11.3% and 24.9%) than the C₃ species F. pringlei and F. cronquistii (4.1% and 7.7%) (Vogan and Sage, 2011). Thus even the least advanced intermediates analyzed in this study run already a low-level C₄ cycle, which assists in recycling the ammonia liberated by GDC in the bundle sheath cells.

The question arises whether amino group transfer initially exclusively happened via amino acid/oxoacid pairs or whether the enzymatic content of C₃ plants immediately supported a shuttle that also involved decarboxylation and carboxylation reactions. C₃ plants have considerable capacity for the decarboxylation of four-carbon organic acids in their bundle sheath cell (Hibberd and Quick, 2002; Brown et al., 2010) and measurements of total leaf NAD-ME and NADP-ME activity in C₃ plants repeatedly demonstrated basal activities for various C₃ species (Wheeler et al., 2005; Aubry et al., 2011; Maier et al., 2011). C₃ plants also accumulate high amounts of organic C₄ acids like malate or fumarate during the day (Zell et al., 2010), which are produced by PEPC, the only enzyme capable of producing C₄ acids de novo. It is tempting to hypothesize that plants use a malate/alanine shuttle to recycle parts of the ammonia liberated by glycine decarboxylation from the very beginning of the C₂ cycle.

Elevating the C₄ cycle activity in a C₂ plant enhances the CO₂ fixation capacity

If the C₄ cycle is superimposed onto a C₂ cycle operating in a C₃–C₄ intermediate plant, the C₂ photosynthesis model predicts a mechanistic interaction between the C₂ and C₄ cycles (Figure 5). When the C₄ cycle is running, the photorespiratory ammonia is recirculated from the bundle sheath to the mesophyll cells by moving malate from the mesophyll to the bundle sheath and transferring alanine back to the mesophyll. This malate/alanine cycling leads to a net transport of ammonia from the bundle sheath into the mesophyll cells. In contrast to the other mechanisms of ammonia recirculation described above, the C₄ cycle does not only lead to a net transport of ammonia from the bundle sheath to the mesophyll but additionally also to a net transport of CO₂ in the opposite direction. Thus CO₂ is transferred from the mesophyll to the bundle sheath without increasing the number of transport processes between the cells. By elevating the CO₂ concentration in the bundle sheath cells the C₄ cycle acts cooperatively with the C₂ cycle. The bundle sheath Rubisco would work under a more elevated CO₂ concentration and thus operate more effectively compared to a pure C₂ plant, leading to an increased biomass production. The C₄ cycle thus has a dual beneficial effect: an efficient nitrogen shuttle is combined with a CO₂ concentrating pump.

To investigate the possible interaction with regard to biomass, a C₄ cycle at the enzyme capacities of C₃ plants was allowed and tested for biomass changes (Figure 5 C). When the C₄ cycle is running with PEPC activities comparable to those found in C₃ Flaveria species, the model already predicts a gain in biomass production compared to the C₂ cycle on its own. Under these conditions, the bulk of photorespiratory ammonia is recycled through a rudimentary C₄ cycle limited by the C₄ cycle flux capacity. The model predicts that biomass production will be further enhanced with higher activity of the C₄ cycle. Consequently, there is permanent positive selection on enhancing the activity of the currently rate limiting enzyme once a C₄ cycle is running.

The evolutionary scenario described above is in good agreement with the Flaveria transcriptome data. We observe gradual increases in the amounts of C₄ transcript with increasing ‘C₄-ness’ of the C₃–C₄ intermediates until the most advanced species F. brownii. The abundance of NADP-ME and Ala-AT transcripts increases faster than the transcript abundance of the other core C₄ genes like PEPC, PPDK, MDH or Asp-AT. This implies that these evolutionary changes were driven by selection on high bundle sheath decarboxylation capacity, consistent with the idea that the C₄ cycle began as an auxiliary pathway to the C₂ cycle to recirculate photorespiratory ammonia. Hence, in this early phase, the main purpose of the C₄ cycle was to provide the ammonia acceptor pyruvate. The C₂ model and its evolutionary implications are consistent with the properties of the C₃–C₄ intermediate Flaveria species including F. brownii, which possess mesophyll Rubisco activity and consequently the C₂ photosynthetic pathway. The next iteration during C₄ evolution in Flaveria must have been the restriction of Rubisco activity to the bundle sheath, making the C₂ cycle obsolete, as observed for the true C₄ Flaveria species.

A scenario for C₄ evolution in the genus Flaveria—a general blueprint for the evolution of C₄ photosynthesis?

The establishment of a photorespiratory CO₂ pump, termed C₂ photosynthesis, is thought to be an important step in C₄ evolution. Recent work has shown how C₃ Flaverias were preconditioned for the evolution of the C₂ pathway and how the C₂ cycle was implemented on the molecular level (Sage et al., 2013; Schulze et al., 2013). Together with the present work, this gives us a detailed picture of what happened in the early and intermediate stages during C₄ evolution in Flaveria.

We have argued that the establishment of the C₂ cycle requires the implementation of at least components of the C₄ pathway, if not the whole pathway. This fact might be a partial explanation for the polyphyletic evolution of C₄ photosynthesis. Only the C₂ cycle has to evolve to set a system on a slippery slope towards C₄ photosynthesis. Nature seems to confirm this idea. So far, 66 independent origins of C₄ photosynthesis could be identified. In contrast, there are only seven known groups with independent origins of C₂ plants and no direct ancestry to C₄ species (Sage et al., 2012). If one assumes that all recent C₄ lineages evolved via C₂ intermediates, which appears likely (Sage et al., 2012; Heckmann et al., 2013; Williams et al., 2013), this would mean that the C₂ pathway evolved 73 times independently and that over 90% of these C₂ plant containing lineages proceeded to C₄ photosynthesis. This indicates that the C₂ photosynthetic pathway must indeed be a strong enabler of C₄ photosynthesis. It will be highly enlightening to analyze these C₂ groups without ancestry to C₄ species, like Moricandia, Steinchisma or Mollugo, to find out in how far they differ from groups that evolved the C₄ pathway and why C₄ evolution may have been hampered in these groups.

The close evolutionary interconnection of the C₂ and the C₄ pathway could be seen as an example of metabolic exaptation (Barve and Wagner, 2013). Exaptation or pre-adaptation was defined as an adaptation involving the co-option of traits that originally evolved for a different purpose (Gould and Vrba, 1982). While both C₂ and C₄ act as carbon shuttles to the bundle sheath cells, the two systems achieve this goal through different biochemical processes. In particular, the amino acid shuttle in the C₂ system evolved to transport nitrogen, and its later use in C₄ photosynthesis to shuttle carbon thus represents a molecular exaptation. Our findings therefore corroborate the general idea that the evolution of complex traits may be accelerated through exaptations (Darwin, 1872; Gould and Vrba, 1982; Barve and Wagner, 2013).

We do not know if the scenario on the early and intermediate stages of evolution described above is limited to the genus Flaveria or if it is valid for C₄ evolution in general. Our prediction of the C₂ pathway being a strong facilitator of C₄ evolution should apply to all C₄ origins, as the integrated model is not specific to Flaveria.

Materials and methods

Plant material

Request a detailed protocol

F. pringlei, F. robusta, F. chloraefolia, F. pubescens, F. anomala, F. ramosissima, F. brownii, F. bidentis and F. trinervia plants were grown in the green house at University of Duesseldorf side-by-side and harvested at four different points of time over the year. The plants were grown in 17-cm pots on soil (C-400 with Cocopor [Stender Erden, Schermbeck Germany] fertilized with 3 g/l Osmocote exact standard 3 to 4 M [Marysville, USA]) with additional light for 16 hr per day until 50 to 60 cm height and before the onset of flowering.

Plants for experiment one were harvested in September, for experiment two in June, for experiment three in October and for experiment four in April. The plant material was immediately frozen in liquid nitrogen, stored at −80°C and used for the following analyses.

RNA isolation, transcriptome sequencing and analysis

Request a detailed protocol

Total RNA was isolated from the second and fourth leaves according to (Westhoff et al., 1991) followed by a DNAse treatment. After phenol/chloroform extraction and precipitation with NaAc and isopropyl alcohol the RNA was dissolved in H₂O. The RNA quality was tested with the Agilent 2100 bioanalyzer. 1 µg of total RNA was used for cDNA library generation, which was accomplished with the TruSeq RNA Sample Preparation Kit (Illumina Inc., San Diego, USA) via the Low-Throughput Protocol (TruSeq RNA Sample Preparation Guide, Illumina Proprietary Catalog # RS-930-2001, Part # 15008136 Rev. A, November 2010). Clusters were generated with the TruSeq SR Cluster Kit v2 according to the Reagent Preparation Guide with the Illumina cBot device. The single read sequencing was performed with the Illumina HiSeq2000.

Sequences of transcripts from genes involved in photorespiraton, C₄ photosynthesis and refixation and recirculation of photorespiratory ammonia were identified among de novo assembled transcripts of F. robusta. De novo assembly was performed with either CLC Genomics Workbench (CLC-Bio, Aaarhus, Denmark) or the Velvet/Oases software package (Schulz et al., 2012) using F. robusta 454 (Gowik et al., 2011) and Illumina reads (this study).

After quality control and processing, Illumina reads were aligned to the F. robusta transcript sequences with the CLC Genomics Workbench using standard parameters. Read mapping against a minimal set of coding sequences (Brautigam et al., 2011) of the TAIR 9 release of the Arabidopsis thaliana genome (http://www.Arabidopsis.org/) was performed using BLAT (Kent, 2002) as described in (Gowik et al., 2011).

The MEV software package (http://www.tm4.org/mev.html) was used for plotting heat maps, hierarchical clustering and principal component analysis.

Protein isolation and quantification

Request a detailed protocol

Total proteins were isolated from plant material harvested together with the material for RNA isolation according to Shen et al. (2007) and quantified using the RC-DC protocol (Bio-Rad Laboratories, Hercules, USA). 30 µg of total protein was electrophoresed on polyacrylamide-SDS gels (Schägger and von Jagow, 1987) and electrophoretically transferred to nitrocellulose membranes (Protran BA85, 0.45 μm; Schleicher & Schuell, Dassel, Germany) for 1 hr with 0.8 mA per cm². Specific primary antibodies were raised against conserved Flaveria peptides (Agrisera Vännäs, Sweden). For the detection of specific proteins the nitrocellulose membranes were incubated with the primary antibodies and a Horseradish peroxidase-conjugated secondary antibody (Sigma-Aldrich, St. Louis, USA). An enhanced chemiluminescent Horseradish peroxide substrate was added and signals were recorded using a Fuji LAS-4000 mini CCD camera system. The signals were quantified with the Multi Gage analysis software (Fujifilm, Tokyo, Japan). As loading control a gel was stained for 45 min with 0.25% Coomassie blue, 50% methanol, 7% acetic acid, and destained in 50% methanol, 7% acetic acid.

Coupling a mechanistic model with a genome-scale metabolic reconstruction

Request a detailed protocol

In order to model the metabolic integration of C₂ and C₄ cycle in the context of leaf metabolism, we conducted Flux Balance Analysis (FBA) based on a genome-scale metabolic reconstruction of C₄ metabolism, C4GEM (Dal'Molin et al., 2010). This reconstruction contains a complex biomass reaction including carbohydrates, cell wall components, amino acids and nucleotides (Dal'Molin et al., 2010).

FBA is a powerful tool to understand the adaptation of metabolism on a genomic scale. Since metabolite concentrations are not modeled explicitly, fluxes related to carbon concentration mechanisms (CCMs) cannot be captured by this constraint-based approach alone. To account for this issue, we coupled the FBA model with a mechanistic model of C₃–C₄ photosynthesis (von Caemmerer, 2000; Heckmann et al., 2013).

C4GEM representing NADP-ME types was provided by the authors and FBA was conducted using this model:

Maximize c^{T} v

subject to $S v = 0$ .

v_{\min, i} \leq v_{i} \leq v_{\max, i}

where c is the vector of coefficients in the objective function, here the leaf biomass production. v is the vector of fluxes through the network reactions, S is the stoichiometric matrix of the metabolic network, and v_min and v_max represent constraints on the respective fluxes.

In order to test hypotheses concerning nitrogen metabolism in C₃–C₄ intermediate plants, S had to be modified. The plasmosdesmatal transport reactions in the original C4GEM model include malate, pyruvate, 3-phosphoglycerate, trioses, phosphates, sucrose, aspartate, alanine, phosphoenolpyruvate, CO₂, and O₂. Reactions were added to S in order to include transport of serine, glycine, glutamate, glutamine, asparagine, threonine, 2-oxoglutarate and water over the mesophyll/bundle sheath interface. Furthermore, the lack of photosystem II in the bundle sheath of certain C₄ plants does not hold in our scenario (Nakamura et al., 2013) and we added a reaction for linear electron transport to the bundle sheath. C4GEM does not contain a reaction for a plastidal NADP-dependent malate dehydrogenase in the bundle sheath; we added this reaction to S.

In addition to the stoichiometric matrix S, the constraints used in C4GEM were modified:

The original constraint on leaf sucrose production was changed to result in an output ratio of sucrose to amino acids of about 5 (Riens et al., 1991). Fixed constraints on production of starch and fatty acids are not appropriate in the coupled framework. Since we are not aware of data that explains how these fluxes scale with net CO₂ assimilation rate, the constraints were removed from the model. Reactions belonging to the GS/GOGAT system were assumed to be irreversible. Nitrogen is available in the form of nitrate as opposed to NH₃ in the original model. Since there is no evidence suggesting mesophyll specificity of PEPC in intermediate Flaveria species, we unconstrained PEPC flux in the bundle sheath.

To couple the genome-scale FBA model with the mechanistic model of carbon fixation, the following reactions were constrained using the values predicted by the mechanistic model: net CO₂ uptake, Rubisco carboxylation and oxygenation in mesophyll and bundle sheath, CO₂ leakage from the bundle sheath, PEPC activity in the mesophyll, activity of NADP-ME in the bundle sheath, plasmodesmatal flux of glycine and serine and decarboxylation by the GDC complex. The lower bound on glycine diffusion (V_min,Gly), serine diffusion (V_min,Ser), and GDC reaction (V_min,GDC) can be obtained from the rate of Rubisco oxygenation in the mesophyll (V_om) and the fraction of photorespiratory CO₂ in the bundle sheath derived from mesophyll oxygenations (ξ):

v_{\min, G l y} = ξ V_{o m}, v_{\min, S e r} = 0.5 ξ V_{o m}, v_{\min, G D C} = 0.5 ξ V_{o m}

The mechanistic model was parameterized to the C₃ state as given in Heckmann et al. (2013), with the exception of the parameter ξ, which was set to a value of 0.98 (i.e., the majority of GDC activity was restricted to the bundle sheath. Derivation from transcriptome data is given in Heckmann et al. (2013)). These constraints on the reactions of the photorespiratory pump are necessary to adequately predict C₂ photosynthesis because of the inability of FBA alone to model CCMs (see discussion above).

In the FBA part of the model, a minimization of total flux (MTF) analysis was conducted in order to narrow down the space of optimal solutions:

Minimize \sum_{i = 1}^{n} w_{i} | v_{i} |

subject to: $S v = 0$ .

v_{\min, i} \leq v_{i} \leq v_{\max, i}

c^{T} v = c^{T} v_{F B A}

where v_FBA is the flux distribution of the FBA optimization described above. w denotes a vector of weights, where plasmodesmatal flux received a higher weighting factor (1.1 for plasmodesmatal exchange, 1 for the remaining reactions). This method implements a simple minimization of protein costs for a given optimal biomass production. The higher weights on plasmodesmatal fluxes account for the trade-off between CO₂ containment in the bundle sheath and metabolite diffusion between the cells. Since this trade-off is difficult to quantify, we conducted a sensitivity analysis by varying the weight on plasmodesmatal transport reactions.

In order to investigate the possible range that fluxes can take while yielding an optimal solution, flux variability analysis was conducted:

For each v_i:

Maximize or Minimize $v_{i}$ .

Subject to: $S v = 0$ .

v_{\min, i} \leq v_{i} \leq v_{\max, i}

c^{T} v = c^{T} v_{F B A}

\sum_{i = 1}^{n} w_{i} | v_{i} | = s_{o p t}

Where s_opt is the minimum for the weighted sum of absolute flux found in the MTF optimization.

All simulations were conducted in the R environment for statistical computing (R Core Team, 2013) using the sybil library (Gelius-Dietrich et al., 2013).

Accession numbers

Request a detailed protocol

The read data have been submitted to the National Center for Biotechnology Information Short Read Archive under accession numbers SRP036880 (F. bidentis), SRP036881 (F. anomala), SRP036883 (F. brownii), SRP036884 (F. chloraefolia), SRP036885 (F. pringlei), SRP037526 (F. pubescens), SRP037527 (F. ramosissima), SRP037528 (F. robusta) and SRP037529 (F. trinervia).

Data availability

The following data sets were generated

1. Mallmann J
2. Heckmann D
3. Bräutigam A
4. Lercher MJ
5. Weber APM
6. Westhoff P
7. Gowik U
(2014) Flaveria bidentis leaf transcriptomes
Publicly available at NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/sra).

http://www.ncbi.nlm.nih.gov/sra/?term=SRP036880
1. Mallmann J
2. Heckmann D
3. Bräutigam A
4. Lercher MJ
5. Weber APM
6. Westhoff P
7. Gowik U
(2014) Flaveria anomala leaf transcriptomes
Publicly available at NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/sra).

http://www.ncbi.nlm.nih.gov/sra/?term=SRP036881
1. Mallmann J
2. Heckmann D
3. Bräutigam A
4. Lercher MJ
5. Weber APM
6. Westhoff P
7. Gowik U
(2014) Flaveria brownii leaf transcriptomes
Publicly available at NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/sra).

http://www.ncbi.nlm.nih.gov/sra/?term=SRP036883
1. Mallmann J
2. Heckmann D
3. Bräutigam A
4. Lercher MJ
5. Weber APM
6. Westhoff P
7. Gowik U
(2014) Flaveria chloreafolia leaf transcriptomes
Publicly available at NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/sra).

http://www.ncbi.nlm.nih.gov/sra/?term=SRP036884
1. Mallmann J
2. Heckmann D
3. Bräutigam A
4. Lercher MJ
5. Weber APM
6. Westhoff P
7. Gowik U
(2014) Flaveria pringlei leaf transcriptomes
Publicly available at NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/sra).

http://www.ncbi.nlm.nih.gov/sra/?term=SRP036885
1. Mallmann J
2. Heckmann D
3. Bräutigam A
4. Lercher MJ
5. Weber APM
6. Westhoff P
7. Gowik U
(2014) Flaveria pubescens leaf transcriptomes
Publicly available at NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/sra).

http://www.ncbi.nlm.nih.gov/sra/?term=SRP037526
1. Mallmann J
2. Heckmann D
3. Bräutigam A
4. Lercher MJ
5. Weber APM
6. Westhoff P
7. Gowik U
(2014) Flaveria ramosissima leaf transcriptomes
Publicly available at NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/sra).

http://www.ncbi.nlm.nih.gov/sra/?term=SRP037527
1. Mallmann J
2. Heckmann D
3. Bräutigam A
4. Lercher MJ
5. Weber APM
6. Westhoff P
7. Gowik U
(2014) Flaveria robusta leaf transcriptomes
Publicly available at NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/sra).

http://www.ncbi.nlm.nih.gov/sra/?term=SRP037528
1. Mallmann J
2. Heckmann D
3. Bräutigam A
4. Lercher MJ
5. Weber APM
6. Westhoff P
7. Gowik U
(2014) Flaveria trinervia leaf transcriptomes
Publicly available at NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/sra).

http://www.ncbi.nlm.nih.gov/sra/?term=SRP037529
1. Mallmann J
2. Heckmann D
3. Bräutigam A
4. Lercher MJ
5. Weber APM
6. Westhoff P
7. Gowik U
(2014) Data from: the role of photorespiration during the evolution of C4 photosynthesis in the genus Flaveria
Available at Dryad Digital Repository under a CC0 Public Domain Dedication.

http://dx.doi.org/10.5061/dryad.q827h

References

1. Anderson LE
(1971) Chloroplast and cytoplasmic enzymes. II. Pea leaf triose phosphate isomerases
Biochimica et Biophysica Acta 235:237–244.

https://doi.org/10.1016/0005-2744(71)90051-9
- Google Scholar
(2011) The role of proteins in C(3) plants prior to their recruitment into the C(4) pathway
Journal of Experimental Botany 62:3049–3059.

https://doi.org/10.1093/jxb/err012
- Google Scholar
1. Barve A
2. Wagner A
(2013) A latent capacity for evolutionary innovation through exaptation in metabolic systems
Nature 500:203–206.

https://doi.org/10.1038/nature12301
- Google Scholar
1. Bauwe H
(1984) Photosynthetic enzyme activities and immunofluorescence studies on the localization of ribulose-1, 5-bisphosphate carboxylase/oxygenase in leaves of C3, C4, and C3−C4 intermediate species of Flaveria (Asteraceae)
Biochemie und Physiologie der Pflanzen 179:253–268.

https://doi.org/10.1016/S0015-3796(84)80041-4
- Google Scholar
Book
1. Bauwe H
(2011)
Photorespiration: the bridge to C4 photosynthesis

In: Raghavendra AS, Sage RF, editors. C4 Photosynthesis and related CO2 concentrating mechanisms. Dordrecht, Netherlands: Springer. pp. 81–108.
- Google Scholar
(2010) Photorespiration: players, partners and origin
Trends in Plant Science 15:330–336.

https://doi.org/10.1016/j.tplants.2010.03.006
- Google Scholar
1. Bauwe H
2. Kolukisaoglu Ü
(2003) Genetic manipulation of glycine decarboxylation
Journal of Experimental Botany 54:1523–1535.

https://doi.org/10.1093/jxb/erg171
- Google Scholar
(1971) Phosphoglycolate production catalyzed by ribulose diphosphate carboxylase
Biochemical and Biophysical Research Communications 45:716–722.

https://doi.org/10.1016/0006-291X(71)90475-X
- Google Scholar
1. Brautigam A
2. Kajala K
3. Wullenweber J
4. Sommer M
5. Gagneul D
6. Weber KL
7. Carr KM
8. Gowik U
9. Mass J
10. Lercher MJ
11. Westhoff P
12. Hibberd JM
13. Weber AP
(2011) An mRNA blueprint for C4 photosynthesis derived from comparative transcriptomics of closely related C3 and C4 species
Plant Physiology 155:142–156.

https://doi.org/10.1104/pp.110.159442
- Google Scholar
1. Brown NJ
2. Palmer BG
3. Stanley S
4. Hajaji H
5. Janacek SH
6. Astley HM
7. Parsley K
8. Kajala K
9. Quick WP
10. Trenkamp S
11. Fernie AR
12. Maurino VG
13. Hibberd JM
(2010) C4 acid decarboxylases required for C4 photosynthesis are active in the mid-vein of the C3 species Arabidopsis thaliana, and are important in sugar and amino acid metabolism
The Plant Journal 61:122–133.

https://doi.org/10.1111/j.1365-313X.2009.04040.x
- Google Scholar
1. Chastain CJ
2. Chollet R
(1989) Interspecific variation in assimilation of (14)CO2 into C4 acids by leaves of C3, C4 and C3-C4 intermediate Flaveria species near the CO2 compensation concentration
Planta 179:81–88.

https://doi.org/10.1007/BF00395774
- Google Scholar
1. Cheng SH
2. Moore BD
3. Edwards GE
4. Ku MSB
(1988) Photosynthesis in Flaveria brownii, a C4-like species. Leaf anatomy, characteristics of CO2 exchange, compartmentation of photosynthetic enzymes, and metabolism of 14CO2
Plant Physiology 87:867–873.

https://doi.org/10.1104/pp.87.4.867
- Google Scholar
(2011) C4 eudicots are not younger than C4 monocots
Journal of Experimental Botany 62:3171–3181.

https://doi.org/10.1093/jxb/err041
- Google Scholar
(2010) C4GEM, a genome-scale metabolic model to study C4 plant metabolism
Plant Physiology 154:1871–1885.

https://doi.org/10.1104/pp.110.166488
- Google Scholar
Book
1. Darwin C
(1872)
On the origin of species by means of natural selection, or the preservation of favoured races in the struggle for life (6th Edn)

London: John Murray.
- Google Scholar
Book
1. Dengler NG
2. Nelson T
(1999)
Leaf structure and development in C4 plants

In: Rowan FS, Russell KM, editors. C4 plant biology. San Diego: Academic Press. pp. 133–172.
- Google Scholar
Book
1. Edwards GE
2. Ku MS
(1987)
Biochemistry of C3-C4 intermediates

In: Hatch MD, Boardman NK, editors. The biochemistry of plants. New York: Academic Press. pp. 275–325.
- Google Scholar
(1991) Climate change and the evolution of C4 photosynthesis
Trends in Ecology & Evolution 6:95–99.

https://doi.org/10.1016/0169-5347(91)90183-X
- Google Scholar
1. Eisenhut M
2. Planchais S
3. Cabassa C
4. Guivarc'h A
5. Justin AM
6. Taconnat L
7. Renou JP
8. Linka M
9. Gagneul D
10. Timm S
11. Bauwe H
12. Carol P
13. Weber AP
(2013) Arabidopsis A BOUT DE SOUFFLE is a putative mitochondrial transporter involved in photorespiratory metabolism and is required for meristem growth at ambient CO(2) levels
The Plant Journal 73:836–849.

https://doi.org/10.1111/tpj.12082
- Google Scholar
1. Fernie AR
2. Bauwe H
3. Eisenhut M
4. Florian A
5. Hanson DT
6. Hagemann M
7. Keech O
8. Mielewczik M
9. Nikoloski Z
10. Peterhänsel C
11. Roje S
12. Sage R
13. Timm S
14. von Cammerer S
15. Weber APM
16. Westhoff P
(2013) Perspectives on plant photorespiratory metabolism
Plant Biology 15:748–753.

https://doi.org/10.1111/j.1438-8677.2012.00693.x
- Google Scholar
1. Furumoto T
2. Yamaguchi T
3. Ohshima-Ichie Y
4. Nakamura M
5. Tsuchida-Iwata Y
6. Shimamura M
7. Ohnishi J
8. Hata S
9. Gowik U
10. Westhoff P
11. Brautigam A
12. Weber AP
13. Izui K
(2011) A plastidial sodium-dependent pyruvate transporter
Nature 476:472–475.

https://doi.org/10.1038/nature10250
- Google Scholar
(2013) sybil - efficient constraint-based modelling in R
BMC Systems Biology 7:125.

https://doi.org/10.1186/1752-0509-7-125
- Google Scholar
1. Gould SJ
2. Vrba ES
(1982)
Exaptation-a missing term in the science of form

Paleobiology 8:4–15.
- Google Scholar
1. Gowik U
2. Brautigam A
3. Weber KL
4. Weber APM
5. Westhoff P
(2011) Evolution of C4 photosynthesis in the genus Flaveria: how many and which genes does it take to make C4?
The Plant Cell 23:2087–2105.

https://doi.org/10.1105/tpc.111.086264
- Google Scholar
Book
1. Haberlandt G
(1904)
Physiologische Pflanzenanatomie (3rd edn)

Leipzig, Germany: Verlag von Wilhelm Engelmann.
- Google Scholar
1. Hatch MD
(1987) C4 photosynthesis: a unique blend of modified biochemistry, anatomy and ultrastructure
Biochimica et Biophysica Acta 895:81–106.

https://doi.org/10.1016/S0304-4173(87)80009-5
- Google Scholar
1. Hatch M
2. Kagawa T
3. Craig S
(1975) Subdivision of C4-pathway species based on differing C4 acid decarboxylating systems and ultrastructural features
Functional Plant Biology 2:111–128.

https://doi.org/10.1071/PP9750111
- Google Scholar
1. Heckmann D
2. Schulze S
3. Denton A
4. Gowik U
5. Westhoff P
6. Weber APM
7. Lercher MJ
(2013) Predicting C4 photosynthesis evolution: modular, individually adaptive steps on a Mount fuji fitness landscape
Cell 153:1579–1588.

https://doi.org/10.1016/j.cell.2013.04.058
- Google Scholar
1. Hibberd JM
2. Quick WP
(2002) Characteristics of C4 photosynthesis in stems and petioles of C3 flowering plants
Nature 415:451–454.

https://doi.org/10.1038/415451a
- Google Scholar
(1988) Enzymic and photosynthetic characteristics of reciprocal F(1) hybrids of Flaveria pringlei (C(3)) and Flaveria brownii (C(4)-Like species)
Plant Physiology 87:484–490.

https://doi.org/10.1104/pp.87.2.484
- Google Scholar
(1988) Glycine decarboxylase is confined to the bundle-sheath cells of leaves of C3‚àíC4 intermediate species
Planta 175:452–459.

https://doi.org/10.1007/bf00393064
- Google Scholar
1. Kent WJ
(2002) BLAT–the BLAST-like alignment tool
Genome Research 12:656–664.

https://doi.org/10.1101/gr.229202
- Google Scholar
1. Kinoshita H
2. Nagasaki J
3. Yoshikawa N
4. Yamamoto A
5. Takito S
6. Kawasaki M
7. Sugiyama T
8. Miyake H
9. Weber AP
10. Taniguchi M
(2011) The chloroplastic 2-oxoglutarate/malate transporter has dual function as the malate valve and in carbon/nitrogen metabolism
The Plant Journal 65:15–26.

https://doi.org/10.1111/j.1365-313X.2010.04397.x
- Google Scholar
(1959) Uncouplers of spinach chloroplast photosynthetic Phosphorylation
Plant Physiology 34:272–277.

https://doi.org/10.1104/pp.34.3.272
- Google Scholar
1. Ku MSB
2. Wu J
3. Dai Z
4. Scott RA
5. Chu C
6. Edwards GE
(1991) Photosynthetic and photorespiratory characteristics of Flaveria species
Plant Physiology 96:518–528.

https://doi.org/10.1104/pp.96.2.518
- Google Scholar
(1995) The regulation and control of photorespiration
Journal of Experimental Botany 46:1397–1414.

https://doi.org/10.1093/jxb/46.special_issue.1397
- Google Scholar
(2011) Malate decarboxylases: evolution and roles of NAD (P)-ME isoforms in species performing C4 and C3 photosynthesis
Journal of Experimental Botany 62:3061–3069.

https://doi.org/10.1093/jxb/err024
- Google Scholar
1. Mallmann J
2. Heckmann D
3. Bräutigam A
4. Lercher MJ
5. Weber APM
6. Westhoff P
7. Gowik U
(2014) Data from: the role of photorespiration during the evolution of C4 photosynthesis in the genus Flaveria
Dryad Digital Repository.

https://doi.org/10.5061/dryad.q827h
- Google Scholar
(2005) Phylogeny of Flaveria (Asteraceae) and inference of C4 photosynthesis evolution
American Journal of Botany 92:1911–1928.

https://doi.org/10.3732/ajb.92.11.1911
- Google Scholar
(1996) The roles of malate and aspartate in C4 photosynthetic metabolism of Flaveria bidentis (L.)
Planta 199:262–269.

https://doi.org/10.1007/BF00196567
- Google Scholar
Book
1. Monson R
2. Rawsthorne S
(2000)
CO2 assimilation in C3-C4 intermediate plants

In: Leegood R, Sharkey T, Caemmerer S, editors. Photosynthesis, Vol 9. Netherlands: Springer. pp. 533–550.
- Google Scholar
1. Monson RK
2. Moore BD
3. Ku MS
4. Edwards GE
(1986) Co-function of C3-and C4-photosynthetic pathways in C3, C 4 and C3-C4 intermediate Flaveria species
Planta 168:493–502.

https://doi.org/10.1007/BF00392268
- Google Scholar
1. Moore B
2. Ku M
3. Edwards G
(1987)
C4 photosynthesis and light-dependent accumulation of inorganic carbon in leaves of C3−C4 and C4 Flaveria species

Australian Journal of Plant Physiology 14:657–668.
- Google Scholar
(1984) Isolation of leaf bundle sheath protoplasts from C4 dicot species and intracellular localization of selected enzymes
Plant Science Letters 35:127–138.

https://doi.org/10.1016/0304-4211(84)90186-X
- Google Scholar
(1993) Coordination of the cell-specific distribution of the four subunits of glycine decarboxylase and of serine hydroxymethyltransferase in leaves of C3-C4 intermediate species from different genera
Planta 190:468–473.

https://doi.org/10.1007/bf00224785
- Google Scholar
1. Nakamura N
2. Iwano M
3. Havaux M
4. Yokota A
5. Munekage YN
(2013) Promotion of cyclic electron transport around photosystem I during the evolution of NADP–malic enzyme-type C4 photosynthesis in the genus Flaveria
New Phytologist 199:832–842.

https://doi.org/10.1111/nph.12296
- Google Scholar
1. Ogren WL
(1984) Photorespiration: pathways, regulation, and Modification
Annual Review of Plant Physiology 35:415–442.

https://doi.org/10.1146/annurev.pp.35.060184.002215
- Google Scholar
(2010) Photorespiration
The Arabidopsis Book 8:e0130.

https://doi.org/10.1199/tab.0130
- Google Scholar
1. Pick TR
2. Brautigam A
3. Schulz MA
4. Obata T
5. Fernie AR
6. Weber AP
(2013) PLGG1, a plastidic glycolate glycerate transporter, is required for photorespiration and defines a unique class of metabolite transporters
Proceedings of the National Academy of Sciences of the United States of America 110:3185–3190.

https://doi.org/10.1073/pnas.1215142110
- Google Scholar
1. Powell M
(1978) Systematics of Flaveria (Flaveriinae-Asteraceae)
Annals of the Missouri Botanical Garden 65:590–636.

https://doi.org/10.2307/2398862
- Google Scholar
1. R Core Team
(2013) R foundation for statistical computing
R foundation for statistical computing, Vienna, Austria, URL http://www.R-project.org/.

http://www.R-project.org/
- Google Scholar
1. Raines CA
(2011) Increasing photosynthetic carbon assimilation in C3 plants to improve crop yield: current and future strategies
Plant Physiology 155:36–42.

https://doi.org/10.1104/pp.110.168559
- Google Scholar
1. Renne P
2. Dressen U
3. Hebbeker U
4. Hille D
5. Flugge UI
6. Westhoff P
7. Weber AP
(2003) The Arabidopsis mutant dct is deficient in the plastidic glutamate/malate translocator DiT2
The Plant Journal 35:316–331.

https://doi.org/10.1046/j.1365-313X.2003.01806.x
- Google Scholar
1. Riens B
2. Lohaus G
3. Heineke D
4. Heldt HW
(1991) Amino acid and sucrose content determined in the cytosolic, chloroplastic, and vacuolar compartments and in the Phloem sap of spinach leaves
Plant Physiology 97:227–233.

https://doi.org/10.1104/pp.97.1.227
- Google Scholar
1. Rumpho ME
2. Ku MS
3. Cheng SH
4. Edwards GE
(1984) Photosynthetic characteristics of C(3)-C(4) intermediate Flaveria species : III. Reduction of photorespiration by a limited C(4) pathway of photosynthesis in Flaveria ramosissima
Plant Physiology 75:993–996.

https://doi.org/10.1104/pp.75.4.993
- Google Scholar
1. Sage RF
(2001) Environmental and evolutionary preconditionsfor the origin and diversification of the C4 photosyntheticsyndrome
Plant Biology 3:202–213.

https://doi.org/10.1055/s-2001-15206
- Google Scholar
1. Sage RF
(2004) The evolution of C4 photosynthesis
New Phytologist 161:341–370.

https://doi.org/10.1111/j.1469-8137.2004.00974.x
- Google Scholar
1. Sage RF
(2013) Photorespiratory compensation: a driver for biological diversity
Plant Biology 15:624–638.

https://doi.org/10.1111/plb.12024
- Google Scholar
(2011) The C4 plant lineages of planet Earth
Journal of Experimental Botany 62:3155–3169.

https://doi.org/10.1093/jxb/err048
- Google Scholar
(2012) Photorespiration and the evolution of C4 Photosynthesis
Annual Review of Plant Biology 63:19–47.

https://doi.org/10.1146/annurev-arplant-042811-105511
- Google Scholar
1. Sage TL
2. Busch FA
3. Johnson DC
4. Friesen PC
5. Stinson CR
6. Stata M
7. Sultmanis S
8. Rahman BA
9. Rawsthorne S
10. Sage RF
(2013) Initial events during the evolution of C4 photosynthesis in C3 species of Flaveria
Plant Physiology 163:1266–1276.

https://doi.org/10.1104/pp.113.221119
- Google Scholar
1. Schägger H
2. von Jagow G
(1987) Tricine-sodium dodecyl sulfate-polyacrylamide gel electrophoresis for the separation of proteins in the range from 1 to 100 kDa
Analytical Biochemistry 166:368–379.

https://doi.org/10.1016/0003-2697(87)90587-2
- Google Scholar
(2012) Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels
Bioinformatics 28:1086–1092.

https://doi.org/10.1093/bioinformatics/bts094
- Google Scholar
1. Schulze S
2. Mallmann J
3. Burscheidt J
4. Koczor M
5. Streubel M
6. Bauwe H
7. Gowik U
8. Westhoff P
(2013) Evolution of C4 photosynthesis in the genus flaveria: establishment of a photorespiratory CO2 pump
The Plant Cell 25:2522–2535.

https://doi.org/10.1105/tpc.113.114520
- Google Scholar
1. Shen Y
2. Khanna R
3. Carle CM
4. Quail PH
(2007) Phytochrome induces rapid PIF5 phosphorylation and degradation in response to red-light activation
Plant Physiology 145:1043–1051.

https://doi.org/10.1104/pp.107.105601
- Google Scholar
1. Varma A
2. Palsson BO
(1994)
Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110

Applied and Environmental Microbiology 60:3724–3731.
- Google Scholar
1. Vogan PJ
2. Sage RF
(2011) Water-use efficiency and nitrogen-use efficiency of C3-C4 intermediate species of Flaveria Juss. (Asteraceae)
Plant, Cell & Environment 34:1415–1430.

https://doi.org/10.1111/j.1365-3040.2011.02340.x
- Google Scholar
Book
1. von Caemmerer S
(2000)
Biochemical models of leaf photosynthesis

Collingwood, Australia: Csiro Publishing.
- Google Scholar
1. Weber AP
2. von Caemmerer S
(2010) Plastid transport and metabolism of C3 and C4 plants–comparative analysis and possible biotechnological exploitation
Current Opinion in Plant Biology 13:257–265.

https://doi.org/10.1016/j.pbi.2010.01.007
- Google Scholar
(1991) Differential accumulation of plastid transcripts encoding photosystem II components in the mesophyll and bundle-sheath cells of monocotyledonous NADP-malic enzyme-type C4 plants
Planta 184:377–388.

https://doi.org/10.1007/bf00195340
- Google Scholar
(2005) A comprehensive analysis of the NADP-malic enzyme gene family of Arabidopsis
Plant Physiology 139:39–51.

https://doi.org/10.1104/pp.105.065953
- Google Scholar
(2013) Phenotypic landscape inference reveals multiple evolutionary paths to C4 photosynthesis
eLife 2:e00961.

https://doi.org/10.7554/eLife.00961
- Google Scholar
1. Zell MB
2. Fahnenstich H
3. Maier A
4. Saigo M
5. Voznesenskaya EV
6. Edwards GE
7. Andreo C
8. Schleifenbaum F
9. Zell C
10. Drincovich MF
11. Maurino VG
(2010) Analysis of Arabidopsis with highly reduced levels of malate and fumarate sheds light on the role of these organic acids as storage carbon molecules
Plant Physiology 152:1251–1262.

https://doi.org/10.1104/pp.109.151795
- Google Scholar

Article and author information

Author details

Julia Mallmann

Institute for Plant Molecular and Developmental Biology, Heinrich-Heine-Universität, Düsseldorf, Germany

Contribution
JM, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article

Contributed equally with
David Heckmann

Competing interests
The authors declare that no competing interests exist.
David Heckmann

Institute for Computer Science, Heinrich-Heine-Universität, Düsseldorf, Germany

Contribution
DH, In silico modeling, Conception and design, Analysis and interpretation of data, Drafting or revising the article

Contributed equally with
Julia Mallmann

Competing interests
The authors declare that no competing interests exist.
Andrea Bräutigam

Institute of Plant Biochemistry, Heinrich-Heine-Universität, Düsseldorf, Germany

Contribution
AB, Conception and design, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Martin J Lercher
1. Institute for Computer Science, Heinrich-Heine-Universität, Düsseldorf, Germany
2. Cluster of Excellence on Plant Sciences (CEPLAS), Heinrich-Heine-Universität, Düsseldorf, Germany
Contribution
MJL, Conception and design, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Andreas PM Weber
1. Institute of Plant Biochemistry, Heinrich-Heine-Universität, Düsseldorf, Germany
2. Cluster of Excellence on Plant Sciences (CEPLAS), Heinrich-Heine-Universität, Düsseldorf, Germany
Contribution
APMW, Conception and design, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Peter Westhoff
1. Institute for Plant Molecular and Developmental Biology, Heinrich-Heine-Universität, Düsseldorf, Germany
2. Cluster of Excellence on Plant Sciences (CEPLAS), Heinrich-Heine-Universität, Düsseldorf, Germany
Contribution
PW, Conception and design, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Udo Gowik

Institute for Plant Molecular and Developmental Biology, Heinrich-Heine-Universität, Düsseldorf, Germany

Contribution
UG, Conception and design, Analysis and interpretation of data, Drafting or revising the article

For correspondence
gowik@uni-duesseldorf.de

Competing interests
The authors declare that no competing interests exist.

Funding

Deutsche Forschungsgemeinschaft (FOR 1186)

Andreas PM Weber
Peter Westhoff

Deutsche Forschungsgemeinschaft (IRTG 1525)

David Heckmann

Deutsche Forschungsgemeinschaft (CRC 680)

Martin J Lercher

Deutsche Forschungsgemeinschaft (EXC 1028)

Martin J Lercher
Andreas PM Weber
Peter Westhoff

European Union (EU Project 3to4)

Andreas PM Weber
Peter Westhoff

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank the ‘Genomics and Transcriptomics laboratory’ of the ‘Biologisch-Medizinischen Forschungszentrum’ (BMFZ) at the Heinrich-Heine-University of Duesseldorf (Germany) for technical support and conducting the Illumina sequencing.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.