The coordination of major events in C4 photosynthesis evolution in the genus Flaveria.

C4 photosynthesis is a remarkable complex trait, elucidations of the evolutionary trajectory of C4 photosynthesis from its ancestral C3 pathway can help us better understand the generic principles of the evolution of complex traits and guide the engineering of C3 crops for higher yields. Here, we used the genus Flaveria that contains C3, C3-C4, C4-like and C4 species as a system to study the evolution of C4 photosynthesis. We first mapped transcript abundance, protein sequence and morphological features onto the phylogenetic tree of the genus Flaveria, and calculated the evolutionary correlation of different features; we then predicted the relative changes of ancestral nodes of those features to illustrate the major events during the evolution of C4 photosynthesis. We found that gene expression and protein sequence showed consistent modification patterns in the phylogenetic tree. High correlation coefficients ranging from 0.46 to 0.9 among gene expression, protein sequence and morphology were observed. The greatest modification of those different features consistently occurred at the transition between C3-C4 species and C4-like species. Our results show highly coordinated changes in gene expression, protein sequence and morphological features, which support evolutionary major events during the evolution of C4 metabolism.

Dynamic changes of genome sizes and gradual gain of cell-specific distribution of C4 enzymes during C4 evolution in genus Flaveria.

C4 plants are believed to have evolved from C3 plants through various C3 -C4 intermediate stages in which a photorespiration-dependent CO2 concentration system known as C2 photosynthesis operates. Genes involved in the C4 cycle were thought to be recruited from orthologs present in C3 species and developed cell-specific expression during C4 evolution. To understand the process of establishing C4 photosynthesis, we performed whole-genome sequencing and investigated expression and mesophyll- or bundle-sheath-cell-specific localization of phosphoenolpyruvate carboxylase (PEPC), NADP-malic enzyme (NADP-ME), pyruvate, orthophosphate dikinase (PPDK) in C3 , C3 -C4 intermediate, C4 -like, and C4 Flaveria species. While genome sizes vary greatly, the number of predicted protein-coding genes was similar among C3 , C3 -C4 intermediate, C4 -like, and C4 Flaveria species. Cell-specific localization of the PEPC, NADP-ME, and PPDK transcripts was insignificant or weak in C3 -C4 intermediate species, whereas these transcripts were expressed cell-type specific in C4 -like species. These results showed that elevation of gene expression and cell-specific control of pre-existing C4 cycle genes in C3 species was involved in C4 evolution. Gene expression was gradually enhanced during C4 evolution, whereas cell-specific control was gained independently of quantitative transcriptional activation during evolution from C3 -C4 intermediate to C4 photosynthesis in genus Flaveria.

Metabolic Labeling of RNAs Uncovers Hidden Features and Dynamics of the Arabidopsis Transcriptome.

Transcriptome analysis by RNA sequencing (RNA-seq) has become an indispensable research tool in modern plant biology. Virtually all RNA-seq studies provide a snapshot of the steady state transcriptome, which contains valuable information about RNA populations at a given time but lacks information about the dynamics of RNA synthesis and degradation. Only a few specialized sequencing techniques, such as global run-on sequencing, have been used to provide information about RNA synthesis rates in plants. Here, we demonstrate that RNA labeling with the modified, nontoxic uridine analog 5-ethynyl uridine (5-EU) in Arabidopsis (Arabidopsis thaliana) seedlings provides insight into plant transcriptome dynamics. Pulse labeling with 5-EU revealed nascent and unstable RNAs, RNA processing intermediates generated by splicing, and chloroplast RNAs. Pulse-chase experiments with 5-EU allowed us to determine RNA stabilities without the need for chemical transcription inhibitors such as actinomycin and cordycepin. Inhibitor-free, genome-wide analysis of polyadenylated RNA stability via 5-EU pulse-chase experiments revealed RNAs with shorter half-lives than those reported after chemical inhibition of transcription. In summary, our results indicate that the Arabidopsis nascent transcriptome contains unstable RNAs and RNA processing intermediates and suggest that polyadenylated RNAs have low stability in plants. Our technique lays the foundation for easy, affordable, nascent transcriptome analysis and inhibitor-free analysis of RNA stability in plants.

The C4 Ppc promoters of many C4 grass species share a common regulatory mechanism for gene expression in the mesophyll cell.

C4 photosynthetic plants have evolved from C3 ancestors and are characterized by differential expression of several hundred genes. Strict compartmentalization of key C4 enzymes either to mesophyll (M) or bundle sheath cells is considered a crucial step towards the evolution of C4 photosynthesis. In this study, we demonstrate that the 5'-flanking sequences of the C4 type phosphoenolpyruvate carboxylase (Ppc) gene from three C4 grass species could drive M-cell-specific expression of a reporter gene in rice. In addition to that, we identified about 450 bp (upstream of their transcription start site) of the analyzed C4 Ppc promoters contain all the essential regulatory elements for driving M-cell-specific expression in rice leaves. Importantly, four motifs of conserved nucleotide sequences (CNSs) were also determined, which are essential for the activity of the promoter. A putative interaction between the CNSs and an unknown upstream element(s) is required for driving M-cell-specific expression. This work identifies the evolutionary conservation of C4 Ppc regulatory mechanisms of multiple closely related C4 grass species.

Reporter-based forward genetic screen to identify bundle sheath anatomy mutants in A. thaliana.

The evolution of C4 photosynthesis proceeded stepwise with each small step increasing the fitness of the plant. An important pre-condition for the introduction of a functional C4 cycle is the photosynthetic activation of the C3 bundle sheath by increasing its volume and organelle number. Therefore, to engineer C4 photosynthesis into existing C3 crops, information about genes that control the bundle sheath cell size and organelle content is needed. However, very little information is known about the genes that could be manipulated to create a more C4 -like bundle sheath. To this end, an ethylmethanesulfonate (EMS)-based forward genetic screen was established in the Brassicaceae C3  species Arabidopsis thaliana. To ensure a high-throughput primary screen, the bundle sheath cells of A. thaliana were labeled using a luciferase (LUC68) or by a chloroplast-targeted green fluorescent protein (sGFP) reporter using a bundle sheath specific promoter. The signal strengths of the reporter genes were used as a proxy to search for mutants with altered bundle sheath anatomy. Here, we show that our genetic screen predominantly identified mutants that were primarily affected in the architecture of the vascular bundle, and led to an increase in bundle sheath volume. By using a mapping-by-sequencing approach the genomic segments that contained mutated candidate genes were identified.

Efficient 2-phosphoglycolate degradation is required to maintain carbon assimilation and allocation in the C4 plant Flaveria bidentis.

Photorespiration is indispensable for oxygenic photosynthesis since it detoxifies and recycles 2-phosphoglycolate (2PG), which is the primary oxygenation product of Rubisco. However, C4 plant species typically display very low rates of photorespiration due to their efficient biochemical carbon-concentrating mechanism. Thus, the broader relevance of photorespiration in these organisms remains unclear. In this study, we assessed the importance of a functional photorespiratory pathway in the C4 plant Flaveria bidentis using knockdown of the first enzymatic step, namely 2PG phosphatase (PGLP). The isolated RNAi lines showed strongly reduced amounts of PGLP protein, but distinct signs of the photorespiratory phenotype only emerged below 5% residual PGLP protein. Lines with this characteristic were stunted in growth, had strongly increased 2PG content, exhibited accelerated leaf senescence, and accumulated high amounts of branched-chain and aromatic amino acids, which are both characteristics of incipient carbon starvation. Oxygen-dependent gas-exchange measurements consistently suggested the cumulative impairment of ribulose-1,5-bisphosphate regeneration with increased photorespiratory pressure. Our results indicate that photorespiration is essential for maintaining high rates of C4 photosynthesis by preventing the 2PG-mediated inhibition of carbon utilization efficiency. However, considerably higher 2PG accumulation can be tolerated compared to equivalent lines of C3 plants due to the differential distribution of specific enzymatic steps between the mesophyll and bundle sheath cells.

Expression of SULTR2;2, encoding a low-affinity sulphur transporter, in the Arabidopsis bundle sheath and vein cells is mediated by a positive regulator.

The bundle sheath provides a conduit linking veins and mesophyll cells. In the C3 plant Arabidopsis thaliana, it also plays important roles in oxidative stress and sulphur metabolism. However, the mechanisms responsible for the patterns of gene expression that underpin these metabolic specializations are poorly understood. Here, we used the Arabidopsis SULTR2;2 gene as a model to better understand mechanisms that restrict expression to the bundle sheath. Deletion analysis indicated that the SULTR2;2 promoter contains a short region necessary for expression in the bundle sheath and veins. This sequence acts as a positive regulator and is tolerant to multiple consecutive deletions indicating considerable redundancy in the cis-elements involved. It is highly conserved in SULTR2;2 genes of the Brassicaceae and is functional in the distantly related C4 species Flaveria bidentis that belongs to the Asteraceae. We conclude that expression of SULTR2;2 in the bundle sheath and veins is underpinned by a highly redundant sequence that likely represents an ancient and conserved mechanism found in families as diverse as the Asteraceae and Brassicaceae.

De novo Transcriptome Assembly and Comparison of C3, C3-C4, and C4 Species of Tribe Salsoleae (Chenopodiaceae).

C4 photosynthesis is a carbon-concentrating mechanism that evolved independently more than 60 times in a wide range of angiosperm lineages. Among other alterations, the evolution of C4 from ancestral C3 photosynthesis requires changes in the expression of a vast number of genes. Differential gene expression analyses between closely related C3 and C4 species have significantly increased our understanding of C4 functioning and evolution. In Chenopodiaceae, a family that is rich in C4 origins and photosynthetic types, the anatomy, physiology and phylogeny of C4, C2, and C3 species of Salsoleae has been studied in great detail, which facilitated the choice of six samples of five representative species with different photosynthetic types for transcriptome comparisons. mRNA from assimilating organs of each species was sequenced in triplicates, and sequence reads were de novo assembled. These novel genetic resources were then analyzed to provide a better understanding of differential gene expression between C3, C2 and C4 species. All three analyzed C4 species belong to the NADP-ME type as most genes encoding core enzymes of this C4 cycle are highly expressed. The abundance of photorespiratory transcripts is decreased compared to the C3 and C2 species. Like in other C4 lineages of Caryophyllales, our results suggest that PEPC1 is the C4-specific isoform in Salsoleae. Two recently identified transporters from the PHT4 protein family may not only be related to the C4 syndrome, but also active in C2 photosynthesis in Salsoleae. In the two populations of the C2 species S. divaricata transcript abundance of several C4 genes are slightly increased, however, a C4 cycle is not detectable in the carbon isotope values. Most of the core enzymes of photorespiration are highly increased in the C2 species compared to both C3 and C4 species, confirming a successful establishment of the C2 photosynthetic pathway. Furthermore, a function of PEP-CK in C2 photosynthesis appears likely, since PEP-CK gene expression is not only increased in S. divaricata but also in C2 species of other groups.

Photosynthesis in C3-C4 intermediate Moricandia species.

Evolution of C4 photosynthesis is not distributed evenly in the plant kingdom. Particularly interesting is the situation in the Brassicaceae, because the family contains no C4 species, but several C3-C4 intermediates, mainly in the genus Moricandia Investigation of leaf anatomy, gas exchange parameters, the metabolome, and the transcriptome of two C3-C4 intermediate Moricandia species, M. arvensis and M. suffruticosa, and their close C3 relative M. moricandioides enabled us to unravel the specific C3-C4 characteristics in these Moricandia lines. Reduced CO2 compensation points in these lines were accompanied by anatomical adjustments, such as centripetal concentration of organelles in the bundle sheath, and metabolic adjustments, such as the balancing of C and N metabolism between mesophyll and bundle sheath cells by multiple pathways. Evolution from C3 to C3-C4 intermediacy was probably facilitated first by loss of one copy of the glycine decarboxylase P-protein, followed by dominant activity of a bundle sheath-specific element in its promoter. In contrast to recent models, installation of the C3-C4 pathway was not accompanied by enhanced activity of the C4 cycle. Our results indicate that metabolic limitations connected to N metabolism or anatomical limitations connected to vein density could have constrained evolution of C4 in Moricandia.

A MEM1-like motif directs mesophyll cell-specific expression of the gene encoding the C4 carbonic anhydrase in Flaveria.

The first two reactions of C4 photosynthesis are catalysed by carbonic anhydrase (CA) and phosphoenolpyruvate carboxylase (PEPC) in the leaf mesophyll (M) cell cytosol. Translatome experiments using a tagged ribosomal protein expressed under the control of M and bundle-sheath (BS) cell-specific promoters showed transcripts encoding CA3 from the C4 species Flaveria bidentis were highly enriched in polysomes from M cells relative to those of the BS. Localisation experiments employing a CA3-green fluorescent protein fusion protein showed F. bidentis CA3 is a cytosolic enzyme. A motif showing high sequence homology to that of the Flaveria M expression module 1 (MEM1) element was identified approximately 2 kb upstream of the F. bidentis and F. trinervia ca3 translation start sites. MEM1 is located in the promoter of C4 Flaveria ppcA genes, which encode the C4-associated PEPC, and is necessary for M-specific expression. No MEM1-like sequence was found in the 4 kb upstream of the C3 species F. pringlei ca3 translation start site. Promoter-reporter fusion experiments demonstrated the region containing the ca3 MEM1-like element also directs M-specific expression. These results support the idea that a common regulatory switch drives the expression of the C4 Flaveria ca3 and ppcA1 genes specifically in M cells.

C3 cotyledons are followed by C4 leaves: intra-individual transcriptome analysis of Salsola soda (Chenopodiaceae).

Some species of Salsoleae (Chenopodiaceae) convert from C3 photosynthesis during the seedling stage to the C4 pathway in adult leaves. This unique developmental transition of photosynthetic pathways offers the exceptional opportunity to follow the development of the derived C4 syndrome from the C3 condition within individual plants, avoiding phylogenetic noise. Here we investigate Salsola soda, a little-studied species from tribe Salsoleae, using an ontogenetic approach. Anatomical sections, carbon isotope (δ13C) values, transcriptome analysis by means of mRNA sequencing, and protein levels of the key C4 enzyme phosphoenolpyruvate carboxylase (PEPC) were examined from seed to adult plant stages. Despite a previous report, our results based on δ13C values, anatomy and transcriptomics clearly indicate a C3 phase during the cotyledon stage. During this stage, the entire transcriptional repertoire of the C4 NADP-malic enzyme type is detected at low levels compared to a significant increase in true leaves. In contrast, abundance of transcripts encoding most of the major photorespiratory enzymes is not significantly decreased in leaves compared to cotyledons. PEPC polypeptide was detected only in leaves, correlating with increased PEPC transcript abundance from the cotyledon to leaf stage.

Glycine decarboxylase in C3, C4 and C3-C4 intermediate species.

The glycine decarboxylase complex (GDC) plays a central role in photorespiration. GDC is localized in the mitochondria and together with serine hydroxymethyltransferase it converts two molecules of glycine to one molecule of serine, CO2 and NH3. Overexpression of GDC subunits in the C3 species Arabidopsis thaliana can increase the metabolic flux through the photorespiratory pathway leading to enhanced photosynthetic efficiency and consequently to an enhanced biomass production of the transgenic plants. Changing the spatial expression patterns of GDC subunits was an important step during the evolution of C3-C4 intermediate and likely also C4 plants. Restriction of the GDC activity to the bundle sheath cells led to the establishment of a photorespiratory CO2 pump.

Most photorespiratory genes are preferentially expressed in the bundle sheath cells of the C4 grass Sorghum bicolor.

One of the hallmarks of C4 plants is the division of labor between two different photosynthetic cell types, the mesophyll and the bundle sheath cells. C4 plants are of polyphyletic origin and, during the evolution of C4 photosynthesis, the expression of thousands of genes was altered and many genes acquired a cell type-specific or preferential expression pattern. Several lines of evidence, including computational modeling and physiological and phylogenetic analyses, indicate that alterations in the expression of a key photorespiration-related gene, encoding the glycine decarboxylase P subunit, was an early and important step during C4 evolution. Restricting the expression of this gene to the bundle sheath led to the establishment of a photorespiratory CO2 pump. We were interested in whether the expression of genes related to photorespiration remains bundle sheath specific in a fully optimized C4 species. Therefore we analyzed the expression of photorespiratory and C4 cycle genes using RNA in situ hybridization and transcriptome analysis of isolated mesophyll and bundle sheath cells in the C4 grass Sorghum bicolor It turns out that the C4 metabolism of Sorghum is based solely on the NADP-dependent malic enzyme pathway. The majority of photorespiratory gene expression, with some important exceptions, is restricted to the bundle sheath.

Photorespiration connects C3 and C4 photosynthesis.

C4 plants evolved independently more than 60 times from C3 ancestors. C4 photosynthesis is a complex trait and its evolution from the ancestral C3 photosynthetic pathway involved the modification of the leaf anatomy and the leaf physiology accompanied by changes in the expression of thousands of genes. Under high temperature, high light, and the current CO2 concentration in the atmosphere, the C4 pathway is more efficient than C3 photosynthesis because it increases the CO2 concentration around the major CO2 fixating enzyme Rubisco. The oxygenase reaction and, accordingly, photorespiration are largely suppressed. In the present review we describe a scenario for C4 evolution that not only includes the avoidance of photorespiration as the major driving force for C4 evolution but also highlights the relevance of changes in the expression of photorespiratory genes in inducing and establishing important phases on the path from C3 to C4.

RNA-Seq based phylogeny recapitulates previous phylogeny of the genus Flaveria (Asteraceae) with some modifications.

BACKGROUND: The genus Flaveria has been extensively used as a model to study the evolution of C4 photosynthesis as it contains C3 and C4 species as well as a number of species that exhibit intermediate types of photosynthesis. The current phylogenetic tree of the genus Flaveria contains 21 of the 23 known Flaveria species and has been previously constructed using a combination of morphological data and three non-coding DNA sequences (nuclear encoded ETS, ITS and chloroplast encoded trnL-F). RESULTS: Here we developed a new strategy to update the phylogenetic tree of 16 Flaveria species based on RNA-Seq data. The updated phylogeny is largely congruent with the previously published tree but with some modifications. We propose that the data collection method provided in this study can be used as a generic method for phylogenetic tree reconstruction if the target species has no genomic information. We also showed that a "F. pringlei" genotype recently used in a number of labs may be a hybrid between F. pringlei (C3) and F. angustifolia (C3-C4). CONCLUSIONS: We propose that the new strategy of obtaining phylogenetic sequences outlined in this study can be used to construct robust trees in a larger number of taxa. The updated Flaveria phylogenetic tree also supports a hypothesis of stepwise and parallel evolution of C4 photosynthesis in the Flavaria clade.

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

C4 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 C2 photosynthesis) is a prerequisite for the evolution of C4. However, a mechanistic model explaining the tight connection between the evolution of C4 and C2 photosynthesis is currently lacking. Here we address this question through comparative transcriptomic and biochemical analyses of closely related C3, C3-C4, and C4 species, combined with Flux Balance Analysis constrained through a mechanistic model of carbon fixation. We show that C2 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 C4 cycle. Our findings thus show how C2 photosynthesis represents a pre-adaptation for the C4 system, where the evolution of the C2 system establishes important C4 components as a side effect.

Evolution of the Phosphoenolpyruvate Carboxylase Protein Kinase Family in C3 and C4 Flaveria spp.

The key enzyme for C4 photosynthesis, Phosphoenolpyruvate Carboxylase (PEPC), evolved from nonphotosynthetic PEPC found in C3 ancestors. In all plants, PEPC is phosphorylated by Phosphoenolpyruvate Carboxylase Protein Kinase (PPCK). However, differences in the phosphorylation pattern exist among plants with these photosynthetic types, and it is still not clear if they are due to interspecies differences or depend on photosynthetic type. The genus Flaveria contains closely related C3, C3-C4 intermediate, and C4 species, which are evolutionarily young and thus well suited for comparative analysis. To characterize the evolutionary differences in PPCK between plants with C3 and C4 photosynthesis, transcriptome libraries from nine Flaveria spp. were used, and a two-member PPCK family (PPCKA and PPCKB) was identified. Sequence analysis identified a number of C3- and C4-specific residues with various occurrences in the intermediates. Quantitative analysis of transcriptome data revealed that PPCKA and PPCKB exhibit inverse diel expression patterns and that C3 and C4 Flaveria spp. differ in the expression levels of these genes. PPCKA has maximal expression levels during the day, whereas PPCKB has maximal expression during the night. Phosphorylation patterns of PEPC varied among C3 and C4 Flaveria spp. too, with PEPC from the C4 species being predominantly phosphorylated throughout the day, while in the C3 species the phosphorylation level was maintained during the entire 24 h. Since C4 Flaveria spp. evolved from C3 ancestors, this work links the evolutionary changes in sequence, PPCK expression, and phosphorylation pattern to an evolutionary phase shift of kinase activity from a C3 to a C4 mode.

Evolution of C4 photosynthesis in the genus flaveria: establishment of a photorespiratory CO2 pump.

C4 photosynthesis is nature's most efficient answer to the dual activity of ribulose-1,5-bisphosphate carboxylase/oxygenase and the resulting loss of CO(2) by photorespiration. Gly decarboxylase (GDC) is the key component of photorespiratory CO(2) release in plants and is active in all photosynthetic tissues of C(3) plants, but only in the bundle sheath cells of C(4) plants. The restriction of GDC to the bundle sheath is assumed to be an essential and early step in the evolution of C(4) photosynthesis, leading to a photorespiratory CO(2) concentrating mechanism. In this study, we analyzed how the P-protein of GDC (GLDP) became restricted to the bundle sheath during the transition from C(3) to C(4) photosynthesis in the genus Flaveria. We found that C(3) Flaveria species already contain a bundle sheath-expressed GLDP gene in addition to a ubiquitously expressed second gene, which became a pseudogene in C(4) Flaveria species. Analyses of C(3)-C(4) intermediate Flaveria species revealed that the photorespiratory CO(2) pump was not established in one single step, but gradually. The knowledge gained by this study sheds light on the early steps in C(4) evolution.

Predicting C4 photosynthesis evolution: modular, individually adaptive steps on a Mount Fuji fitness landscape.

An ultimate goal of evolutionary biology is the prediction and experimental verification of adaptive trajectories on macroevolutionary timescales. This aim has rarely been achieved for complex biological systems, as models usually lack clear correlates of organismal fitness. Here, we simulate the fitness landscape connecting two carbon fixation systems: C3 photosynthesis, used by most plant species, and the C4 system, which is more efficient at ambient CO2 levels and elevated temperatures and which repeatedly evolved from C3. Despite extensive sign epistasis, C4 photosynthesis is evolutionarily accessible through individually adaptive steps from any intermediate state. Simulations show that biochemical subtraits evolve in modules; the order and constitution of modules confirm and extend previous hypotheses based on species comparisons. Plant-species-designated C3-C4 intermediates lie on predicted evolutionary trajectories, indicating that they indeed represent transitory states. Contrary to expectations, we find no slowdown of adaptation and no diminishing fitness gains along evolutionary trajectories.