Evolution of c4 phosphoenolpyruvate carboxylase. Genes and proteins: a case study with the genus Flaveria

Testing the kinetic tradeoff between bicarbonate versus phosphoenolpyruvate affinity and glucose-6 phosphate response of phosphoenolpyruvate carboxylase from two C4 grasses

Phosphoenolpyruvate (PEP) carboxylase (PEPC) has an anaplerotic role in central plant metabolism but also initiates the carbon concentrating mechanism during C4 photosynthesis. The C4 PEPC has different binding affinities (Km) for PEP (K0.5PEP) and HCO3- (K0.5HCO3), and allosteric regulation by glucose-6-phosphate (G6-P) compared to non-photosynthetic isoforms. These differences are linked to specific changes in amino acids within PEPC. For example, region II (residues 302-433, Zea mays numbering) has been identified as important for G6-P regulation and within this region residue 353 may be conserved in C4 PEPC enzymes. Additionally, residue 780 influences the C4 PEPC kinetic properties and may interact with region II as well as residue 353 to influence G6-P regulation. We test the hypothesis that variation within region II, including residue 353, and their interactions with residue 780 influence the kinetic and allosteric regulation by G6-P of two C4 PEPC isozymes from two C4 grasses. The data does not support a kinetic tradeoff between K0.5HCO3 and K0.5PEP in these PEPC isozymes. Additionally, these enzymes had different response to G6-P that was only partially attributed to region II, residue 353 and residue 780. This data provides new insights into factors influencing the kinetic variation of C4 PEPC isozymes.

Agrobacterium-mediated transient transformation of Flaveria bidentis leaves: a novel method to examine the evolution of C4 photosynthesis

The genus Flaveria has been studied extensively as a model for the evolution of C4 photosynthesis. Thus far, molecular analyses in this genus have been limited due to a dearth of genomic information and the lack of a rapid and efficient transformation protocol. Since their development, Agrobacterium-mediated transient transformation protocols have been instrumental in understanding many biological processes in a range of plant species. However, this technique has not been applied to the genus Flaveria. Here, an efficient protocol for the Agrobacterium-mediated transient transformation of the leaves of the C4 species Flaveria bidentis is presented. This technique has the distinct advantages of rapid turnaround, the ability to co-transform with multiple constructs, and the capacity to assay coding and non-coding regions of Flaveria genomes in a homologous context. To illustrate the utility of this protocol, the quantitative transcriptional regulation of phosphoenolpyruvate carboxylase, the primary carboxylase of C4 plants, was investigated. A 24 bp region in the ppcA1 proximal promoter was found to elicit high levels of reporter gene expression. The Agrobacterium-mediated transient transformation of F. bidentis leaves will accelerate the understanding of the biology and evolution of C4 photosynthesis in the genus Flaveria as well as in other C4 lineages.

C4 Phosphoenolpyruvate Carboxylase: Evolution and transcriptional regulation

Photosynthetic phosphoenolpyruvate carboxylase (PEPC) catalyses the irreversible carboxylation of phosphoenolpyruvate (PEP), producing oxaloacetate (OAA). This enzyme catalyses the first step of carbon fixation in C4 photosynthesis, contributing to the high photosynthetic efficiency of C4 plants. PEPC is also involved in replenishing tricarboxylic acid cycle intermediates, such as OAA, being involved in the C/N balance. In plants, PEPCs are classified in two types: bacterial type (BTPC) and plant-type (PTPC), which includes photosynthetic and non-photosynthetic PEPCs. During C4 evolution, photosynthetic PEPCs evolved independently. C4 PEPCs evolved to be highly expressed and active in a spatial-specific manner. Their gene expression pattern is also regulated by developmental cues, light, circadian clock as well as adverse environmental conditions. However, the gene regulatory networks controlling C4 PEPC gene expression, namely its cell-specificity, are largely unknown. Therefore, after an introduction to the evolution of PEPCs, this review aims to discuss the current knowledge regarding the transcriptional regulation of C4 PEPCs, focusing on cell-specific and developmental expression dynamics, light and circadian regulation, as well as response to abiotic stress. In conclusion, this review aims to highlight the evolution, transcriptional regulation by different signals and importance of PEPC in C4 photosynthesis and its potential as tool for crop improvement.

Investigating the Unique Ability of Trichodesmium To Fix Carbon and Nitrogen Simultaneously Using MiMoSA

The open ocean is an extremely competitive environment, partially due to the dearth of nutrients. Trichodesmium erythraeum, a marine diazotrophic cyanobacterium, is a keystone species in the ocean due to its ability to fix nitrogen and leak 30 to 50% into the surrounding environment, providing a valuable source of a necessary macronutrient to other species. While there are other diazotrophic cyanobacteria that play an important role in the marine nitrogen cycle, Trichodesmium is unique in its ability to fix both carbon and nitrogen simultaneously during the day without the use of specialized cells called heterocysts to protect nitrogenase from oxygen. Here, we use the advanced modeling framework called multiscale multiobjective systems analysis (MiMoSA) to investigate how Trichodesmium erythraeum can reduce dimolecular nitrogen to ammonium in the presence of oxygen. Our simulations indicate that nitrogenase inhibition is best modeled as Michealis-Menten competitive inhibition and that cells along the filament maintain microaerobia using high flux through Mehler reactions in order to protect nitrogenase from oxygen. We also examined the effect of location on metabolic flux and found that cells at the end of filaments operate in distinctly different metabolic modes than internal cells despite both operating in a photoautotrophic mode. These results give us important insight into how this species is able to operate photosynthesis and nitrogen fixation simultaneously, giving it a distinct advantage over other diazotrophic cyanobacteria because they can harvest light directly to fuel the energy demand of nitrogen fixation. IMPORTANCE Trichodesmium erythraeum is a marine cyanobacterium responsible for approximately half of all biologically fixed nitrogen, making it an integral part of the global nitrogen cycle. Interestingly, unlike other nitrogen-fixing cyanobacteria, Trichodesmium does not use temporal or spatial separation to protect nitrogenase from oxygen poisoning; instead, it operates photosynthesis and nitrogen fixation reactions simultaneously during the day. Unfortunately, the exact mechanism the cells utilize to operate carbon and nitrogen fixation simultaneously is unknown. Here, we use an advanced metabolic modeling framework to investigate and identify the most likely mechanisms Trichodesmium uses to protect nitrogenase from oxygen. The model predicts that cells operate in a microaerobic mode, using both respiratory and Mehler reactions to dramatically reduce intracellular oxygen concentrations.

High-quality ice plant reference genome analysis provides insights into genome evolution and allows exploration of genes involved in the transition from C3 to CAM pathways

Ice plant (Mesembryanthemum crystallinum), a member of the Aizoaceae family, is a typical halophyte crop and a model plant for studying the mechanism of transition from C3 photosynthesis to crassulacean acid metabolism (CAM). Here, we report a high-quality chromosome-level ice plant genome sequence. This 98.05% genome sequence is anchored to nine chromosomes, with a total length of 377.97 Mb and an N50 scaffold of 40.45 Mb. Almost half of the genome (48.04%) is composed of repetitive sequences, and 24 234 genes have been annotated. Subsequent to the ancient whole-genome triplication (WGT) that occurred in eudicots, there has been no recent whole-genome duplication (WGD) or WGT in ice plants. However, we detected a novel WGT event that occurred in the same order in Simmondsia chinensis, which was previously overlooked. Our findings revealed that ice plants have undergone chromosome rearrangements and gene removal during evolution. Combined with transcriptome and comparative genomic data and expression verification, we identified several key genes involved in the CAM pathway and constructed a comprehensive network. As the first genome of the Aizoaceae family to be released, this report will provide a rich data resource for comparative and functional genomic studies of Aizoaceae, especially for studies on salt tolerance and C3-to-CAM transitions to improve crop yield and resistance.

Conservation and Divergence of Phosphoenolpyruvate Carboxylase Gene Family in Cotton

Phosphoenolpyruvate carboxylase (PEPC) is an important enzyme in plants, which regulates carbon flow through the TCA cycle and controls protein and oil biosynthesis. Although it is important, there is little research on PEPC in cotton, the most important fiber crop in the world. In this study, a total of 125 PEPCs were identified in 15 Gossypium genomes. All PEPC genes in cotton are divided into six groups and each group generally contains one PEPC member in each diploid cotton and two in each tetraploid cotton. This suggests that PEPC genes already existed in cotton before their divergence. There are additional PEPC sub-groups in other plant species, suggesting the different evolution and natural selection during different plant evolution. PEPC genes were independently evolved in each cotton sub-genome. During cotton domestication and evolution, certain PEPC genes were lost and new ones were born to face the new environmental changes and human being needs. The comprehensive analysis of collinearity events and selection pressure shows that genome-wide duplication and fragment duplication are the main methods for the expansion of the PEPC family, and they continue to undergo purification selection during the evolutionary process. PEPC genes were widely expressed with temporal and spatial patterns. The expression patterns of PEPC genes were similar in G. hirsutum and G. barbadense with a slight difference. PEPC2A and 2D were highly expressed in cotton reproductive tissues, including ovule and fiber at all tested developmental stages in both cultivated cottons. However, PEPC1A and 1D were dominantly expressed in vegetative tissues. Abiotic stress also induced the aberrant expression of PEPC genes, in which PEPC1 was induced by both chilling and salinity stresses while PEPC5 was induced by chilling and drought stresses. Each pair (A and D) of PEPC genes showed the similar response to cotton development and different abiotic stress, suggesting the similar function of these PEPCs no matter their origination from A or D sub-genome. However, some divergence was also observed among their origination, such as PEPC5D was induced but PEPC5A was inhibited in G. barbadense during drought treatment, suggesting that a different organized PEPC gene may evolve different functions during cotton evolution. During cotton polyploidization, the homologues genes may refunction and play different roles in different situations.

Mining for allelic gold: finding genetic variation in photosynthetic traits in crops and wild relatives

Improvement of photosynthetic traits in crops to increase yield potential and crop resilience has recently become a major breeding target. Synthetic biology and genetic technologies offer unparalleled opportunities to create new genetics for photosynthetic traits driven by existing fundamental knowledge. However, large 'gene bank' collections of germplasm comprising historical collections of crop species and their relatives offer a wealth of opportunities to find novel allelic variation in the key steps of photosynthesis, to identify new mechanisms and to accelerate genetic progress in crop breeding programmes. Here we explore the available genetic resources in food and fibre crops, strategies to selectively target allelic variation in genes underpinning key photosynthetic processes, and deployment of this variation via gene editing in modern elite material.

Global analysis of switchgrass (Panicum virgatum L.) transcriptomes in response to interactive effects of drought and heat stresses

BACKGROUND

Sustainable production of high-quality feedstock has been of great interest in bioenergy research. Despite the economic importance, high temperatures and water deficit are limiting factors for the successful cultivation of switchgrass in semi-arid areas. There are limited reports on the molecular basis of combined abiotic stress tolerance in switchgrass, particularly the combination of drought and heat stress. We used transcriptomic approaches to elucidate the changes in the response of switchgrass to drought and high temperature simultaneously.

RESULTS

We conducted solely drought treatment in switchgrass plant Alamo AP13 by withholding water after 45 days of growing. For the combination of drought and heat effect, heat treatment (35 °C/25 °C day/night) was imposed after 72 h of the initiation of drought. Samples were collected at 0 h, 72 h, 96 h, 120 h, 144 h, and 168 h after treatment imposition, total RNA was extracted, and RNA-Seq conducted. Out of a total of 32,190 genes, we identified 3912, as drought (DT) responsive genes, 2339 and 4635 as, heat (HT) and drought and heat (DTHT) responsive genes, respectively. There were 209, 106, and 220 transcription factors (TFs) differentially expressed under DT, HT and DTHT respectively. Gene ontology annotation identified the metabolic process as the significant term enriched in DTHT genes. Other biological processes identified in DTHT responsive genes included: response to water, photosynthesis, oxidation-reduction processes, and response to stress. KEGG pathway enrichment analysis on DT and DTHT responsive genes revealed that TFs and genes controlling phenylpropanoid pathways were important for individual as well as combined stress response. For example, hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase (HCT) from the phenylpropanoid pathway was induced by single DT and combinations of DTHT stress.

CONCLUSION

Through RNA-Seq analysis, we have identified unique and overlapping genes in response to DT and combined DTHT stress in switchgrass. The combination of DT and HT stress may affect the photosynthetic machinery and phenylpropanoid pathway of switchgrass which negatively impacts lignin synthesis and biomass production of switchgrass. The biological function of genes identified particularly in response to DTHT stress could further be confirmed by techniques such as single point mutation or RNAi.

Metabolic profiles in C3, C3-C4 intermediate, C4-like, and C4 species in the genus Flaveria

C4 photosynthesis concentrates CO2 around Rubisco in the bundle sheath, favouring carboxylation over oxygenation and decreasing photorespiration. This complex trait evolved independently in >60 angiosperm lineages. Its evolution can be investigated in genera such as Flaveria (Asteraceae) that contain species representing intermediate stages between C3 and C4 photosynthesis. Previous studies have indicated that the first major change in metabolism probably involved relocation of glycine decarboxylase and photorespiratory CO2 release to the bundle sheath and establishment of intercellular shuttles to maintain nitrogen stoichiometry. This was followed by selection for a CO2-concentrating cycle between phosphoenolpyruvate carboxylase in the mesophyll and decarboxylases in the bundle sheath, and relocation of Rubisco to the latter. We have profiled 52 metabolites in nine Flaveria species and analysed 13CO2 labelling patterns for four species. Our results point to operation of multiple shuttles, including movement of aspartate in C3-C4 intermediates and a switch towards a malate/pyruvate shuttle in C4-like species. The malate/pyruvate shuttle increases from C4-like to complete C4 species, accompanied by a rise in ancillary organic acid pools. Our findings support current models and uncover further modifications of metabolism along the evolutionary path to C4 photosynthesis in the genus Flaveria.

The limiting factors and regulatory processes that control the environmental responses of C3, C3-C4 intermediate, and C4 photosynthesis

Here, we describe a model of C3, C3-C4 intermediate, and C4 photosynthesis that is designed to facilitate quantitative analysis of physiological measurements. The model relates the factors limiting electron transport and carbon metabolism, the regulatory processes that coordinate these metabolic domains, and the responses to light, carbon dioxide, and temperature. It has three unique features. First, mechanistic expressions describe how the cytochrome b6f complex controls electron transport in mesophyll and bundle sheath chloroplasts. Second, the coupling between the mesophyll and bundle sheath expressions represents how feedback regulation of Cyt b6f coordinates electron transport and carbon metabolism. Third, the temperature sensitivity of Cyt b6f is differentiated from that of the coupling between NADPH, Fd, and ATP production. Using this model, we present simulations demonstrating that the light dependence of the carbon dioxide compensation point in C3-C4 leaves can be explained by co-occurrence of light saturation in the mesophyll and light limitation in the bundle sheath. We also present inversions demonstrating that population-level variation in the carbon dioxide compensation point in a Type I C3-C4 plant, Flaveria chloraefolia, can be explained by variable allocation of photosynthetic capacity to the bundle sheath. These results suggest that Type I C3-C4 intermediate plants adjust pigment and protein distributions to optimize the glycine shuttle under different light and temperature regimes, and that the malate and aspartate shuttles may have originally functioned to smooth out the energy supply and demand associated with the glycine shuttle. This model has a wide range of potential applications to physiological, ecological, and evolutionary questions.

Differential Expression, Tissue-Specific Distribution, and Posttranslational Controls of Phosphoenolpyruvate Carboxylase

Phosphoenolpyruvate carboxylase (PEPC) is a ubiquitous cytosolic enzyme, which is crucial for plant carbon metabolism. PEPC participates in photosynthesis by catalyzing the initial fixation of atmospheric CO2 and is abundant in both C4 and crassulacean acid metabolism leaves. PEPC is differentially expressed at different stages of plant development, mostly in leaves, but also in developing seeds. PEPC is known to show tissue-specific distribution in leaves and in other plant organs, such as roots, stems, and flowers. Plant PEPC undergoes reversible phosphorylation and monoubiquitination, which are posttranslational modifications playing important roles in regulatory processes and in protein localization. Phosphorylation activates the PEPC enzyme, making it more sensitive to glucose-6-phosphate and less sensitive to malate or aspartate. PEPC phosphorylation is known to be diurnally regulated and delicately changed in response to various environmental stimuli, in addition to light. PEPCs belong to a small gene family encoding several plant-type and distantly related bacterial-type PEPCs. This paper provides a minireview of the general information on PEPCs in both C4 and C3 plants.

Genome-Wide Identification and Analysis of the Phosphoenolpyruvate Carboxylase Gene Family in Suaeda aralocaspica, an Annual Halophyte With Single-Cellular C4 Anatomy

Phosphoenolpyruvate carboxylase (PEPC) plays pivotal roles in the carbon fixation of photosynthesis and a variety of metabolic and stress pathways. Suaeda aralocaspica belongs to a single-cellular C₄ species and carries out a photosynthetic pathway in an unusually elongated chlorenchyma cell, which is expected to have PEPCs with different characteristics. To identify the different isoforms of PEPC genes in S. aralocaspica and comparatively analyze their expression and regulation patterns as well as the biochemical and enzymatic properties in this study, we characterized a bacterial-type PEPC (BTPC; SaPEPC-4) in addition to the two plant-type PEPCs (PTPCs; SaPEPC-1 and SaPEPC-2) using a genome-wide identification. SaPEPC-4 presented a lower expression level in all test combinations with an unknown function; two SaPTPCs showed distinct subcellular localizations and different spatiotemporal expression patterns but positively responded to abiotic stresses. Compared to SaPEPC-2, the expression of SaPEPC-1 specifically in chlorenchyma cell tissues was much more active with the progression of development and under various stresses, particularly sensitive to light, implying the involvement of SaPEPC-1 in a C₄ photosynthetic pathway. In contrast, SaPEPC-2 was more like a non-photosynthetic PEPC. The expression trends of two SaPTPCs in response to light, development, and abiotic stresses were also matched with the changes in PEPC activity in vivo (native) or in vitro (recombinant), and the biochemical properties of the two recombinant SaPTPCs were similar in response to various effectors while the catalytic efficiency, substrate affinity, and enzyme activity of SaPEPC-2 were higher than that of SaPEPC-1 in vitro. All the different properties between these two SaPTPCs might be involved in transcriptional (e.g., specific cis-elements), posttranscriptional [e.g., 5'-untranslated region (5'-UTR) secondary structure], or translational (e.g., PEPC phosphorylation/dephosphorylation) regulatory events. The comparative studies on the different isoforms of the PEPC gene family in S. aralocaspica may help to decipher their exact role in C₄ photosynthesis, plant growth/development, and stress resistance.

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

C₄ plants are believed to have evolved from C₃ plants through various C₃ -C₄ intermediate stages in which a photorespiration-dependent CO₂ concentration system known as C₂ photosynthesis operates. Genes involved in the C₄ cycle were thought to be recruited from orthologs present in C₃ species and developed cell-specific expression during C₄ evolution. To understand the process of establishing C₄ 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 C₃ , C₃ -C₄ intermediate, C₄ -like, and C₄ Flaveria species. While genome sizes vary greatly, the number of predicted protein-coding genes was similar among C₃ , C₃ -C₄ intermediate, C₄ -like, and C₄ Flaveria species. Cell-specific localization of the PEPC, NADP-ME, and PPDK transcripts was insignificant or weak in C₃ -C₄ intermediate species, whereas these transcripts were expressed cell-type specific in C₄ -like species. These results showed that elevation of gene expression and cell-specific control of pre-existing C₄ cycle genes in C₃ species was involved in C₄ evolution. Gene expression was gradually enhanced during C₄ evolution, whereas cell-specific control was gained independently of quantitative transcriptional activation during evolution from C₃ -C₄ intermediate to C₄ photosynthesis in genus Flaveria.

Developing Portulaca oleracea as a model system for functional genomics analysis of C4/CAM photosynthesis

Previously regarded as an intriguing photosynthetic curiosity, the occurrence of C4 and Crassulacean acid metabolism (CAM) photosynthesis within a single organism has recently emerged as a source of information for future biotechnological use. Among C4/CAM facultative species, Portulaca oleracea L. has been used as a model for biochemical and gene expression analysis of C4/CAM under field and laboratory conditions. In the present work, we focussed on developing molecular tools to facilitate functional genomics studies in this species, from the optimisation of RNA isolation protocols to a method for stable genetic transformation. Eleven variations of RNA extraction procedures were tested and compared for RNA quantity and quality. Also, 7 sample sets comprising total RNA from hormonal and abiotic stress treatments, distinct plant organs, leaf developmental stages, and subspecies were used to select, among 12 reference genes, the most stable reference genes for RT-qPCR analysis of each experimental condition. Furthermore, different explant sources, Agrobacterium tumefaciens strains, and regeneration and antibiotic selection media were tested in various combinations to optimise a protocol for stable genetic transformation of P. oleracea. Altogether, we provide essential tools for functional gene analysis in the context of C4/CAM photosynthesis, including an efficient RNA isolation method, preferred reference genes for RT-qPCR normalisation for a range of experimental conditions, and a protocol to produce P. oleracea stable transformants using A. tumefaciens.

Kinetic variation in grass phosphoenolpyruvate carboxylases provides opportunity to enhance C4 photosynthetic efficiency

The high rates of photosynthesis and the carbon-concentrating mechanism (CCM) in C₄ plants are initiated by the enzyme phosphoenolpyruvate (PEP) carboxylase (PEPC). The flow of inorganic carbon into the CCM of C₄ plants is driven by PEPC's affinity for bicarbonate (K_HCO3 ), which can be rate limiting when atmospheric CO₂ availability is restricted due to low stomatal conductance. We hypothesize that natural variation in K_HCO3 across C₄ plants is driven by specific amino acid substitutions to impact rates of C₄ photosynthesis under environments such as drought that restrict stomatal conductance. To test this hypothesis, we measured K_HCO3 from 20 C₄ grasses to compare kinetic properties with specific amino acid substitutions. There was nearly a twofold range in K_HCO3 across these C₄ grasses (24.3 ± 1.5 to 46.3 ± 2.4 μm), which significantly impacts modeled rates of C₄ photosynthesis. Additionally, molecular engineering of a low-HCO₃^- affinity PEPC identified key domains that confer variation in K_HCO3 . This study advances our understanding of PEPC kinetics and builds the foundation for engineering increased-HCO₃^- affinity and C₄ photosynthetic efficiency in important C₄ crops.

Physiological and Proteomic Signatures Reveal Mechanisms of Superior Drought Resilience in Pearl Millet Compared to Wheat

Presently, pearl millet and wheat are belonging to highly important cereal crops. Pearl millet, however, is an under-utilized crop, despite its superior resilience to drought and heat stress in contrast to wheat. To investigate this in more detail, we performed comparative physiological screening and large scale proteomics of drought stress responses in drought-tolerant and susceptible genotypes of pearl millet and wheat. These chosen genotypes are widely used in breeding and farming practices. The physiological responses demonstrated large differences in the regulation of root morphology and photosynthetic machinery, revealing a stay-green phenotype in pearl millet. Subsequent tissue-specific proteome analysis of leaves, roots and seeds led to the identification of 12,558 proteins in pearl millet and wheat under well-watered and stress conditions. To allow for this comparative proteome analysis and to provide a platform for future functional proteomics studies we performed a systematic phylogenetic analysis of all orthologues in pearl millet, wheat, foxtail millet, sorghum, barley, brachypodium, rice, maize, Arabidopsis, and soybean. In summary, we define (i) a stay-green proteome signature in the drought-tolerant pearl millet phenotype and (ii) differential senescence proteome signatures in contrasting wheat phenotypes not capable of coping with similar drought stress. These different responses have a significant effect on yield and grain filling processes reflected by the harvest index. Proteome signatures related to root morphology and seed yield demonstrated the unexpected intra- and interspecies-specific biochemical plasticity for stress adaptation for both pearl millet and wheat genotypes. These quantitative reference data provide tissue- and phenotype-specific marker proteins of stress defense mechanisms which are not predictable from the genome sequence itself and have potential value for marker-assisted breeding beyond genome assisted breeding.

Recurrent sequence evolution after independent gene duplication

Background

Convergent and parallel evolution provide unique insights into the mechanisms of natural selection. Some of the most striking convergent and parallel (collectively recurrent) amino acid substitutions in proteins are adaptive, but there are also many that are selectively neutral. Accordingly, genome-wide assessment has shown that recurrent sequence evolution in orthologs is chiefly explained by nearly neutral evolution. For paralogs, more frequent functional change is expected because additional copies are generally not retained if they do not acquire their own niche. Yet, it is unknown to what extent recurrent sequence differentiation is discernible after independent gene duplications in different eukaryotic taxa.

Results

We develop a framework that detects patterns of recurrent sequence evolution in duplicated genes. This is used to analyze the genomes of 90 diverse eukaryotes. We find a remarkable number of families with a potentially predictable functional differentiation following gene duplication. In some protein families, more than ten independent duplications show a similar sequence-level differentiation between paralogs. Based on further analysis, the sequence divergence is found to be generally asymmetric. Moreover, about 6% of the recurrent sequence evolution between paralog pairs can be attributed to recurrent differentiation of subcellular localization. Finally, we reveal the specific recurrent patterns for the gene families Hint1/Hint2, Sco1/Sco2 and vma11/vma3.

Conclusions

The presented methodology provides a means to study the biochemical underpinning of functional differentiation between paralogs. For instance, two abundantly repeated substitutions are identified between independently derived Sco1 and Sco2 paralogs. Such identified substitutions allow direct experimental testing of the biological role of these residues for the repeated functional differentiation. We also uncover a diverse set of families with recurrent sequence evolution and reveal trends in the functional and evolutionary trajectories of this hitherto understudied phenomenon.

What Matters for C4 Transporters: Evolutionary Changes of Phosphoenolpyruvate Transporter for C4 Photosynthesis

C4 photosynthesis is a complex trait that evolved from its ancestral C3 photosynthesis by recruiting pre-existing genes. These co-opted genes were changed in many aspects compared to their counterparts in C3 species. Most of the evolutionary changes of the C4 shuttle enzymes are well characterized, however, evolutionary changes for the recruited metabolite transporters are less studied. Here we analyzed the evolutionary changes of the shuttle enzyme phosphoenolpyruvate (PEP) transporter (PPT) during its recruitment from C3 to C4 photosynthesis. Our analysis showed that among the two PPT paralogs PPT1 and PPT2, PPT1 was the copy recruited for C4 photosynthesis in multiple C4 lineages. During C4 evolution, PPT1 gained increased transcript abundance, shifted its expression from predominantly in root to in leaf and from bundle sheath cell to mesophyll cell, and gained more rapid and long-lasting responsiveness to light. Modifications occurred in both regulatory and coding regions in C4 PPT1 as compared to C3 PPT1, however, the PEP transporting function of PPT1 remained. We found that PPT1 of a Flaveria C4 species recruited a MEM1 B submodule in the promoter region, which might be related to the increased transcript abundance of PPT1 in C4 mesophyll cells. The case study of PPT further suggested that high transcript abundance in a proper location is of high priority for PPT to support C4 function.

Regulation and Evolution of C4 Photosynthesis

C₄ photosynthesis evolved multiple times independently from ancestral C₃ photosynthesis in a broad range of flowering land plant families and in both monocots and dicots. The evolution of C₄ photosynthesis entails the recruitment of enzyme activities that are not involved in photosynthetic carbon fixation in C₃ plants to photosynthesis. This requires a different regulation of gene expression as well as a different regulation of enzyme activities in comparison to the C₃ context. Further, C₄ photosynthesis relies on a distinct leaf anatomy that differs from that of C₃, requiring a differential regulation of leaf development in C₄. We summarize recent progress in the understanding of C₄-specific features in evolution and metabolic regulation in the context of C₄ photosynthesis.

Physiological and Proteomic Signatures Reveal Mechanisms of Superior Drought Resilience in Pearl Millet Compared to Wheat

Presently, pearl millet and wheat are belonging to highly important cereal crops. Pearl millet, however, is an under-utilized crop, despite its superior resilience to drought and heat stress in contrast to wheat. To investigate this in more detail, we performed comparative physiological screening and large scale proteomics of drought stress responses in drought-tolerant and susceptible genotypes of pearl millet and wheat. These chosen genotypes are widely used in breeding and farming practices. The physiological responses demonstrated large differences in the regulation of root morphology and photosynthetic machinery, revealing a stay-green phenotype in pearl millet. Subsequent tissue-specific proteome analysis of leaves, roots and seeds led to the identification of 12,558 proteins in pearl millet and wheat under well-watered and stress conditions. To allow for this comparative proteome analysis and to provide a platform for future functional proteomics studies we performed a systematic phylogenetic analysis of all orthologues in pearl millet, wheat, foxtail millet, sorghum, barley, brachypodium, rice, maize, Arabidopsis, and soybean. In summary, we define (i) a stay-green proteome signature in the drought-tolerant pearl millet phenotype and (ii) differential senescence proteome signatures in contrasting wheat phenotypes not capable of coping with similar drought stress. These different responses have a significant effect on yield and grain filling processes reflected by the harvest index. Proteome signatures related to root morphology and seed yield demonstrated the unexpected intra- and interspecies-specific biochemical plasticity for stress adaptation for both pearl millet and wheat genotypes. These quantitative reference data provide tissue- and phenotype-specific marker proteins of stress defense mechanisms which are not predictable from the genome sequence itself and have potential value for marker-assisted breeding beyond genome assisted breeding.

A single serine to alanine substitution decreases bicarbonate affinity of phosphoenolpyruvate carboxylase in C4Flaveria trinervia

Phosphoenolpyruvate (PEP) carboxylase (PEPc) catalyzes the first committed step of C4 photosynthesis generating oxaloacetate from bicarbonate (HCO3-) and PEP. It is hypothesized that PEPc affinity for HCO3- has undergone selective pressure for a lower KHCO3 (Km for HCO3-) to increase the carbon flux entering the C4 cycle, particularly during conditions that limit CO2 availability. However, the decrease in KHCO3 has been hypothesized to cause an unavoidable increase in KPEP (Km for PEP). Therefore, the amino acid residue S774 in the C4 enzyme, which has been shown to increase KPEP, should lead to a decrease in KHCO3. Several studies reported the effect S774 has on KPEP; however, the influence of this amino acid substitution on KHCO3 has not been tested. To test these hypotheses, membrane-inlet mass spectrometry (MIMS) was used to measure the KHCO3 of the photosynthetic PEPc from the C4Flaveria trinervia and the non-photosynthetic PEPc from the C3F. pringlei. The cDNAs for these enzymes were overexpressed and purified from the PEPc-less PCR1 Escherichia coli strain. Our work in comparison with previous reports suggests that KHCO3 and KPEP are linked by specific amino acids, such as S774; however, these kinetic parameters respond differently to the tested allosteric regulators, malate and glucose-6-phosphate.

The role of alanine and aspartate aminotransferases in C4 photosynthesis

Alanine and aspartate are essential transfer metabolites for C₄ species of the NAD-malic enzyme and phosphoenolpyruvate carboxykinase subtype. To some degree both amino acids are also part of the metabolite shuttle in NADP-malic enzyme plants. In comparison with C₃ species, the majority of C₄ species are therefore characterised by enhanced expression and activity of alanine and aspartate aminotransferases (AT) in the photosynthetically active tissue. Both enzymes exist in multiple copies and have been found in different subcellular compartments. We tested whether different C₄ species show preferential recruitment of enzymes from specific lineages and subcellular compartments. Phylogenetic analysis of alanine and aspartate ATs from a variety of monocot and eudicot C₄ species and their C₃ relatives was combined with subcellular prediction tools and analysis of the subsequent transcript amounts in mature leaves. Recruitment of aspartate AT from a specific subcellular compartment was strongly connected to the biochemical subtype. Deviation from the main model was however observed in Gynandropsis gynandra. The configuration of alanine AT generally differed in monocot and eudicot species. C₄ monocots recruited an alanine AT from a specific cytosolic branch, but eudicots use alanine AT copies from a mitochondrial branch. Generally, plants display high plasticity in the setup of the C₄ pathway. Beside the common models for the different C₄ subtypes, individual solutions were found for plant groups or lineages.

Rapid and pervasive development- and tissue-specific homeolog expression partitioning in newly formed inter-subspecific rice segmental allotetraploids

Background

In diverse plant taxa, whole-genome duplication (WGD) events are major sources of phenotypic novelty. Studies of gene expression in synthetic polyploids have shown immediate expression and functional partitioning of duplicated genes among different tissues. Many studies of the tissue-specific homeolog expression partitioning have focused on allopolyploids that have very different parental genomes, while few studies have focused on autopolyploids or allopolyploids that have similar parental genomes.

Results

In this study, we used a set of reciprocal F1 hybrids and synthetic tetraploids constructed from subspecies (japonica and indica) of Asian rice (Oryza sativa L.) as a model to gain insights into the expression partitioning of homeologs among tissues in a developmental context. We assayed the tissue-specific silencing (TSS) of the parental homeologs of 30 key genes in the hybrids and tetraploids relative to the in vitro “hybrids” (parental mixes) using Sequenom MassARRAY. We found that the parental mix and synthetic tetraploids had higher frequencies of homeolog TSS than the F1, revealing an instantaneous role of WGD on homeolog expression partitioning.

Conclusions

Our observations contradicted those of previous studies in which newly formed allopolyploids had a low TSS frequency, similar to that of F1 hybrids, suggesting that the impact of WGD on homeolog expression requires a longer time to manifest. In addition, we found that the TSS frequency in the tetraploids varied at different growth stages and that roots had a much higher frequency of TSS than leaves, which indicated that developmental and metabolic traits may influence the expression states of duplicated genes in newly formed plant polyploids.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5150-7) contains supplementary material, which is available to authorized users.

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.

Loss of the Chloroplast Transit Peptide from an Ancestral C3 Carbonic Anhydrase Is Associated with C4 Evolution in the Grass Genus Neurachne

Neurachne is the only known grass lineage containing closely related C3, C3-C4 intermediate, and C4 species, making it an ideal taxon with which to study the evolution of C4 photosynthesis in the grasses. To begin dissecting the molecular changes that led to the evolution of C4 photosynthesis in this group, the complementary DNAs encoding four distinct β-carbonic anhydrase (CA) isoforms were characterized from leaf tissue of Neurachne munroi (C4), Neurachne minor (C3-C4), and Neurachne alopecuroidea (C3). Two genes (CA1 and CA2) each encode two different isoforms: CA1a/CA1b and CA2a/CA2b. Transcript analyses found that CA1 messenger RNAs were significantly more abundant than transcripts from the CA2 gene in the leaves of each species examined, constituting ∼99% of all β-CA transcripts measured. Localization experiments using green fluorescent protein fusion constructs showed that, while CA1b is a cytosolic CA in all three species, the CA1a proteins are differentially localized. The N. alopecuroidea and N. minor CA1a isoforms were imported into chloroplasts of Nicotiana benthamiana leaf cells, whereas N. munroi CA1a localized to the cytosol. Sequence analysis indicated an 11-amino acid deletion in the amino terminus of N. munroi CA1a relative to the C3 and C3-C4 proteins, suggesting that chloroplast targeting of CA1a is the ancestral state and that loss of a functional chloroplast transit peptide in N. munroi CA1a is associated with the evolution of C4 photosynthesis in Neurachne spp. Remarkably, this mechanism is homoplastic with the evolution of the C4-associated CA in the dicotyledonous genus Flaveria, although the actual mutations in the two lineages differ.

Regulatory gateways for cell-specific gene expression in C4 leaves with Kranz anatomy

C₄ photosynthesis is a carbon-concentrating mechanism that increases delivery of carbon dioxide to RuBisCO and as a consequence reduces photorespiration. The C₄ pathway is therefore beneficial in environments that promote high photorespiration. This pathway has evolved many times, and involves restricting gene expression to either mesophyll or bundle sheath cells. Here we review the regulatory mechanisms that control cell-preferential expression of genes in the C₄ cycle. From this analysis, it is clear that the C₄ pathway has a complex regulatory framework, with control operating at epigenetic, transcriptional, post-transcriptional, translational, and post-translational levels. Some genes of the C₄ pathway are regulated at multiple levels, and we propose that this ensures robust expression in each cell type. Accumulating evidence suggests that multiple genes of the C₄ pathway may share the same regulatory mechanism. The control systems for C₄ photosynthesis gene expression appear to operate in C₃ plants, and so it appears that pre-existing mechanisms form the basis of C₄ photosynthesis gene expression.

Histone Acetylation Modifications Affect Tissue-Dependent Expression of Poplar Homologs of C4 Photosynthetic Enzyme Genes

Histone modifications play important roles in regulating the expression of C₄ photosynthetic genes. Given that all enzymes required for the C₄ photosynthesis pathway are present in C₃ plants, it has been hypothesized that this expression regulatory mechanism has been conserved. However, the relationship between histone modification and the expression of homologs of C₄ photosynthetic enzyme genes has not been well determined in C₃ plants. In the present study, we cloned nine hybrid poplar (Populus simonii × Populus nigra) homologs of maize (Zea mays) C₄ photosynthetic enzyme genes, carbonic anhydrase (CA), pyruvate orthophosphate dikinase (PPDK), phosphoenolpyruvate carboxykinase (PCK), and phosphoenolpyruvate carboxylase (PEPC), and investigated the correlation between the expression levels of these genes and the levels of promoter histone acetylation modifications in four vegetative tissues. We found that poplar homologs of C₄ homologous genes had tissue-dependent expression patterns that were mostly well-correlated with the level of histone acetylation modification (H3K9ac and H4K5ac) determined by chromatin immunoprecipitation assays. Treatment with the histone deacetylase inhibitor trichostatin A further confirmed the role of histone acetylation in the regulation of the nine target genes. Collectively, these results suggest that both H3K9ac and H4K5ac positively regulate the tissue-dependent expression pattern of the PsnCAs, PsnPPDKs, PsnPCKs, and PsnPEPCs genes and that this regulatory mechanism seems to be conserved among the C₃ and C₄ species. Our findings provide new insight that will aid efforts to modify the expression pattern of these homologs of C₄ genes to engineer C₄ plants from C₃ plants.

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

The first two reactions of C₄ 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 C₄ 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 C₄ Flaveria ppcA genes, which encode the C₄-associated PEPC, and is necessary for M-specific expression. No MEM1-like sequence was found in the 4 kb upstream of the C₃ 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 C₄ Flaveria ca3 and ppcA1 genes specifically in M cells.

Cloning of PEPC-1 from a C4 halophyte Suaeda aralocaspica without Kranz anatomy and its recombinant enzymatic activity in responses to abiotic stresses

Phosphoenolpyruvate carboxylase (PEPC) is a key enzyme of C4 photosynthetic pathway and plays an important biochemical role in higher plants and micro organisms. To gain understanding of the role of PEPC in stress adaptation in plant, we cloned PEPC gene from Suaeda aralocaspica, a C4 species without Kranz anatomy, and performed a series of experiments with PEPC gene expressed in Escherichia coli under various abiotic stresses. Results showed that, based on the homology cloning and 5'-RACE technique, the full-length cDNA sequence of PEPC (2901 bp) from S. aralocaspica was obtained, which shares the typical conserved domains to documented PEPCs and was identified as PEPC-1 in accord to the reported partial sequence (ppc-1) in S. aralocaspica. qRT-PCR analysis revealed the expression patterns of PEPC-1 and PEPC-2 (known as ppc-2, another plant type of PEPC) in S. aralocaspica, suggesting that PEPC-1 was up-regulated during seed germination and under NaCl stress, and presented higher level in chlorenchyma than other tissues, which were significantly different with PEPC-2. Afterwards, PEPC-1 was recombinant in E. coli (pET-28a-PEPC) and expressed as an approximate 110 kDa protein. Under various abiotic stresses, the recombinant E. coli strain harboring with PEPC-1 showed significant advantage in growth at 400-800 mmol L(-1) NaCl, 10-20% PEG6000, 25 and 30 °C lower temperature, 50-200 μmol L(-1) methyl viologen, and pH 5.0 and 9.0 condition, compared to control. Further analysis of the enzymatic characteristics of the recombinant PEPC-1 suggests that it was the higher enzyme activity of PEPC-1 which might confer the stress tolerance to E. coli. We speculate that over expression of PEPC-1 is probably related to regulation of oxaloacetate (OAA) in tricarboxylic acid (TCA) cycle in E. coli, which may contribute to further understanding of the physiological function of PEPC in S. aralocaspica.

Posttranslational Protein Modifications in Plant Metabolism

Posttranslational modifications (PTMs) of proteins greatly expand proteome diversity, increase functionality, and allow for rapid responses, all at relatively low costs for the cell. PTMs play key roles in plants through their impact on signaling, gene expression, protein stability and interactions, and enzyme kinetics. Following a brief discussion of the experimental and bioinformatics challenges of PTM identification, localization, and quantification (occupancy), a concise overview is provided of the major PTMs and their (potential) functional consequences in plants, with emphasis on plant metabolism. Classic examples that illustrate the regulation of plant metabolic enzymes and pathways by PTMs and their cross talk are summarized. Recent large-scale proteomics studies mapped many PTMs to a wide range of metabolic functions. Unraveling of the PTM code, i.e. a predictive understanding of the (combinatorial) consequences of PTMs, is needed to convert this growing wealth of data into an understanding of plant metabolic regulation.

Comparison of plant-type phosphoenolpyruvate carboxylases from rice: identification of two plant-specific regulatory regions of the allosteric enzyme

Phosphoenolpyruvate carboxylase (PEPC) is a key enzyme of primary metabolism in bacteria, algae and vascular plants, and it undergoes allosteric regulation by various metabolic effectors. Rice (Oryza sativa) has five plant-type PEPCs, four cytosolic and one chloroplastic. We investigated their kinetic properties using recombinant proteins and found that, like most plant-type PEPCs, rice cytosolic isozymes were activated by glucose 6-phosphate and by alkaline pH. In contrast, no such activation was observed for the chloroplastic isozyme, Osppc4. In addition, Osppc4 showed low affinity for the substrate phosphoenolpyruvate (PEP) and very low sensitivities to allosteric inhibitors aspartate and glutamate. By comparing the isozyme amino acid sequences and three-dimensional structures simulated on the basis of the reported crystal structures, we identified two regions where Osppc4 has unique features that can be expected to affect its kinetic properties. One is the N-terminal extension; replacement of the extension of Osppc2a (cytosolic) with that from Osppc4 reduced the aspartate and glutamate sensitivities to about one-tenth of the wild-type values but left the PEP affinity unaffected. The other is the N-terminal loop, in which a conserved lysine at the N-terminal end is replaced with a glutamate-alanine pair in Osppc4. Replacement of the lysine of Osppc2a with glutamate-alanine lowered the PEP affinity to a quarter of the wild-type level (down to the Osppc4 level), without affecting inhibitor sensitivity. Both the N-terminal extension and the N-terminal loop are specific to plant-type PEPCs, suggesting that plant-type isozymes acquired these regions so that their activity could be regulated properly at the sites where they function.

Phylogeny of C4-photosynthesis enzymes based on algal transcriptomic and genomic data supports an archaeal/proteobacterial origin and multiple duplication for most C4-related genes

Both Calvin-Benson-Bassham (C₃) and Hatch-Slack (C₄) cycles are most important autotrophic CO₂ fixation pathways on today’s Earth. C₃ cycle is believed to be originated from cyanobacterial endosymbiosis. However, studies on evolution of different biochemical variants of C₄ photosynthesis are limited to tracheophytes and origins of C₄-cycle genes are not clear till now. Our comprehensive analyses on bioinformatics and phylogenetics of novel transcriptomic sequencing data of 21 rhodophytes and 19 Phaeophyceae marine species and public genomic data of more algae, tracheophytes, cyanobacteria, proteobacteria and archaea revealed the origin and evolution of C₄ cycle-related genes. Almost all of C₄-related genes were annotated in extensive algal lineages with proteobacterial or archaeal origins, except for phosphoenolpyruvate carboxykinase (PCK) and aspartate aminotransferase (AST) with both cyanobacterial and archaeal/proteobacterial origin. Notably, cyanobacteria may not possess complete C₄ pathway because of the flawed annotation of pyruvate orthophosphate dikinase (PPDK) genes in public data. Most C₄ cycle-related genes endured duplication and gave rise to functional differentiation and adaptation in different algal lineages. C₄-related genes of NAD-ME (NAD-malic enzyme) and PCK subtypes exist in most algae and may be primitive ones, while NADP-ME (NADP-malic enzyme) subtype genes might evolve from NAD-ME subtype by gene duplication in chlorophytes and tracheophytes.

The draft genome and transcriptome of Amaranthus hypochondriacus: a C4 dicot producing high-lysine edible pseudo-cereal

Grain amaranths, edible C4 dicots, produce pseudo-cereals high in lysine. Lysine being one of the most limiting essential amino acids in cereals and C4 photosynthesis being one of the most sought-after phenotypes in protein-rich legume crops, the genome of one of the grain amaranths is likely to play a critical role in crop research. We have sequenced the genome and transcriptome of Amaranthus hypochondriacus, a diploid (2n = 32) belonging to the order Caryophyllales with an estimated genome size of 466 Mb. Of the 411 linkage single-nucleotide polymorphisms (SNPs) reported for grain amaranths, 355 SNPs (86%) are represented in the scaffolds and 74% of the 8.6 billion bases of the sequenced transcriptome map to the genomic scaffolds. The genome of A. hypochondriacus, codes for at least 24,829 proteins, shares the paleohexaploidy event with species under the superorders Rosids and Asterids, harbours 1 SNP in 1,000 bases, and contains 13.76% of repeat elements. Annotation of all the genes in the lysine biosynthetic pathway using comparative genomics and expression analysis offers insights into the high-lysine phenotype. As the first grain species under Caryophyllales and the first C4 dicot genome reported, the work presented here will be beneficial in improving crops and in expanding our understanding of angiosperm evolution.

Positive selection of Kranz and non-Kranz C4 phosphoenolpyruvate carboxylase amino acids in Suaedoideae (Chenopodiaceae)

In subfamily Suaedoideae, four independent gains of C4 photosynthesis are proposed, which includes two parallel origins of Kranz anatomy (sections Salsina and Schoberia) and two independent origins of single-cell C4 anatomy (Bienertia and Suaeda aralocaspica). Additional phylogenetic support for this hypothesis was generated from sequence data of the C-terminal portion of the phosphoenolpyruvate carboxylase (PEPC) gene used in C4 photosynthesis (ppc-1) in combination with previous sequence data. ppc-1 sequence was generated for 20 species in Suaedoideae and two outgroup Salsola species that included all types of C4 anatomies as well as two types of C3 anatomies. A branch-site test for positively selected codons was performed using the software package PAML. From labelling of the four branches where C4 is hypothesized to have developed (foreground branches), residue 733 (maize numbering) was identified to be under positive selection with a posterior probability >0.99 and residue 868 at the >0.95 interval using Bayes empirical Bayes (BEB). When labelling all the branches within C4 clades, the branch-site test identified 13 codons to be under selection with a posterior probability >0.95 by BEB; this is discussed considering current information on functional residues. The signature C4 substitution of an alanine for a serine at position 780 in the C-terminal end (which is considered a major determinant of affinity for PEP) was only found in four of the C4 species sampled, while eight of the C4 species and all the C3 species have an alanine residue; indicating that this substitution is not a requirement for C4 function.

Multiple isoforms of phosphoenolpyruvate carboxylase in the Orchidaceae (subtribe Oncidiinae): implications for the evolution of crassulacean acid metabolism

Phosphoenolpyruvate carboxylase (PEPC) catalyses the initial fixation of atmospheric CO2 into oxaloacetate and subsequently malate. Nocturnal accumulation of malic acid within the vacuole of photosynthetic cells is a typical feature of plants that perform crassulacean acid metabolism (CAM). PEPC is a ubiquitous plant enzyme encoded by a small gene family, and each member encodes an isoform with specialized function. CAM-specific PEPC isoforms probably evolved from ancestral non-photosynthetic isoforms by gene duplication events and subsequent acquisition of transcriptional control elements that mediate increased leaf-specific or photosynthetic-tissue-specific mRNA expression. To understand the patterns of functional diversification related to the expression of CAM, ppc gene families and photosynthetic patterns were characterized in 11 closely related orchid species from the subtribe Oncidiinae with a range of photosynthetic pathways from C3 photosynthesis (Oncidium cheirophorum, Oncidium maduroi, Rossioglossum krameri, and Oncidium sotoanum) to weak CAM (Oncidium panamense, Oncidium sphacelatum, Gomesa flexuosa and Rossioglossum insleayi) and strong CAM (Rossioglossum ampliatum, Trichocentrum nanum, and Trichocentrum carthagenense). Phylogenetic analysis revealed the existence of two main ppc lineages in flowering plants, two main ppc lineages within the eudicots, and three ppc lineages within the Orchidaceae. Our results indicate that ppc gene family expansion within the Orchidaceae is likely to be the result of gene duplication events followed by adaptive sequence divergence. CAM-associated PEPC isoforms in the Orchidaceae probably evolved from several independent origins.

Towards an integrative model of C4 photosynthetic subtypes: insights from comparative transcriptome analysis of NAD-ME, NADP-ME, and PEP-CK C4 species

C4 photosynthesis affords higher photosynthetic carbon conversion efficiency than C3 photosynthesis and it therefore represents an attractive target for engineering efforts aiming to improve crop productivity. To this end, blueprints are required that reflect C4 metabolism as closely as possible. Such blueprints have been derived from comparative transcriptome analyses of C3 species with related C4 species belonging to the NAD-malic enzyme (NAD-ME) and NADP-ME subgroups of C4 photosynthesis. However, a comparison between C3 and the phosphoenolpyruvate carboxykinase (PEP-CK) subtype of C4 photosynthesis is still missing. An integrative analysis of all three C4 subtypes has also not been possible to date, since no comparison has been available for closely related C3 and PEP-CK C4 species. To generate the data, the guinea grass Megathyrsus maximus, which represents a PEP-CK species, was analysed in comparison with a closely related C3 sister species, Dichanthelium clandestinum, and with publicly available sets of RNA-Seq data from C4 species belonging to the NAD-ME and NADP-ME subgroups. The data indicate that the core C4 cycle of the PEP-CK grass M. maximus is quite similar to that of NAD-ME species with only a few exceptions, such as the subcellular location of transfer acid production and the degree and pattern of up-regulation of genes encoding C4 enzymes. One additional mitochondrial transporter protein was associated with the core cycle. The broad comparison identified sucrose and starch synthesis, as well as the prevention of leakage of C4 cycle intermediates to other metabolic pathways, as critical components of C4 metabolism. Estimation of intercellular transport fluxes indicated that flux between cells is increased by at least two orders of magnitude in C4 species compared with C3 species. In contrast to NAD-ME and NADP-ME species, the transcription of photosynthetic electron transfer proteins was unchanged in PEP-CK. In summary, the PEP-CK blueprint of M. maximus appears to be simpler than those of NAD-ME and NADP-ME plants.

Getting the most out of natural variation in C4 photosynthesis

C4 photosynthesis is a complex trait that has a high degree of natural variation, involving anatomical and biochemical changes relative to the ancestral C3 state. It has evolved at least 66 times across a variety of lineages and the evolutionary route from C3 to C4 is likely conserved but not necessarily genetically identical. As such, a variety of C4 species are needed to identify what is fundamental to the C4 evolutionary process in a global context. In order to identify the genetic components of C4 form and function, a number of species are used as genetic models. These include Zea mays (maize), Sorghum bicolor (sorghum), Setaria viridis (Setaria), Flaveria bidentis, and Cleome gynandra. Each of these species has different benefits and challenges associated with its use as a model organism. Here, we propose that RNA profiling of a large sampling of C4, C3-C4, and C3 species, from as many lineages as possible, will allow identification of candidate genes necessary and sufficient to confer C4 anatomy and/or biochemistry. Furthermore, C4 model species will play a critical role in the functional characterization of these candidate genes and identification of their regulatory elements, by providing a platform for transformation and through the use of gene expression profiles in mesophyll and bundle sheath cells and along the leaf developmental gradient. Efforts should be made to sequence the genomes of F. bidentis and C. gynandra and to develop congeneric C3 species as genetic models for comparative studies. In combination, such resources would facilitate discovery of common and unique C4 regulatory mechanisms across genera.

Phosphoenolpyruvate carboxylase identified as a key enzyme in erythrocytic Plasmodium falciparum carbon metabolism

Phospoenolpyruvate carboxylase (PEPC) is absent from humans but encoded in the Plasmodium falciparum genome, suggesting that PEPC has a parasite-specific function. To investigate its importance in P. falciparum, we generated a pepc null mutant (D10^Δpepc), which was only achievable when malate, a reduction product of oxaloacetate, was added to the growth medium. D10^Δpepc had a severe growth defect in vitro, which was partially reversed by addition of malate or fumarate, suggesting that pepc may be essential in vivo. Targeted metabolomics using ¹³C-U-D-glucose and ¹³C-bicarbonate showed that the conversion of glycolytically-derived PEP into malate, fumarate, aspartate and citrate was abolished in D10^Δpepc and that pentose phosphate pathway metabolites and glycerol 3-phosphate were present at increased levels. In contrast, metabolism of the carbon skeleton of ¹³C,¹⁵N-U-glutamine was similar in both parasite lines, although the flux was lower in D10^Δpepc; it also confirmed the operation of a complete forward TCA cycle in the wild type parasite. Overall, these data confirm the CO₂ fixing activity of PEPC and suggest that it provides metabolites essential for TCA cycle anaplerosis and the maintenance of cytosolic and mitochondrial redox balance. Moreover, these findings imply that PEPC may be an exploitable target for future drug discovery.

The genome of the human malaria parasite Plasmodium falciparum encodes a protein called phosphoenolpyruvate carboxylase (PEPC) absent from the human host. PEPC is known to fix CO₂ to generate metabolites used for energy metabolism in plants and bacteria, but its function in malaria parasites remained an enigma. Our study aimed to elucidate the role and importance of PEPC in P. falciparum in its host red blood cell by generating a gene deletion mutant in P. falciparum. This was only achievable in the presence of high concentrations of malate were added to the culture medium. The mutant generated (D10^Δpepc) had a severe growth defect, which was rescued partially by malate or fumarate (but not any other downstream metabolites), suggesting that they feed into the same metabolic pathway. Using heavy isotope labelled ¹³C-U-D-glucose and ¹³C-bicarbonate we showed that PECP has an important role in intermediary carbon metabolism and is vital for the maintenance of cytosolic and mitochondrial redox balance. Together these findings imply that PEPC may be an exploitable target for future drug discovery.

Photosynthetic gene expression in higher plants

Within the chloroplasts of higher plants and algae, photosynthesis converts light into biological energy, fueling the assimilation of atmospheric carbon dioxide into biologically useful molecules. Two major steps, photosynthetic electron transport and the Calvin-Benson cycle, require many gene products encoded from chloroplast as well as nuclear genomes. The expression of genes in both cellular compartments is highly dynamic and influenced by a diverse range of factors. Light is the primary environmental determinant of photosynthetic gene expression. Working through photoreceptors such as phytochrome, light regulates photosynthetic genes at transcriptional and posttranscriptional levels. Other processes that affect photosynthetic gene expression include photosynthetic activity, development, and biotic and abiotic stress. Anterograde (from nucleus to chloroplast) and retrograde (from chloroplast to nucleus) signaling insures the highly coordinated expression of the many photosynthetic genes between these different compartments. Anterograde signaling incorporates nuclear-encoded transcriptional and posttranscriptional regulators, such as sigma factors and RNA-binding proteins, respectively. Retrograde signaling utilizes photosynthetic processes such as photosynthetic electron transport and redox signaling to influence the expression of photosynthetic genes in the nucleus. The basic C3 photosynthetic pathway serves as the default form used by most of the plant species on earth. High temperature and water stress associated with arid environments have led to the development of specialized C4 and CAM photosynthesis, which evolved as modifications of the basic default expression program. The goal of this article is to explain and summarize the many gene expression and regulatory processes that work together to support photosynthetic function in plants.

Greater efficiency of photosynthetic carbon fixation due to single amino-acid substitution

The C4-photosynthetic carbon cycle is an elaborated addition to the classical C3-photosynthetic pathway, which improves solar conversion efficiency. The key enzyme in this pathway, phosphoenolpyruvate carboxylase, has evolved from an ancestral non-photosynthetic C3 phosphoenolpyruvate carboxylase. During evolution, C4 phosphoenolpyruvate carboxylase has increased its kinetic efficiency and reduced its sensitivity towards the feedback inhibitors malate and aspartate. An open question is the molecular basis of the shift in inhibitor tolerance. Here we show that a single-point mutation is sufficient to account for the drastic differences between the inhibitor tolerances of C3 and C4 phosphoenolpyruvate carboxylases. We solved high-resolution X-ray crystal structures of a C3 phosphoenolpyruvate carboxylase and a closely related C4 phosphoenolpyruvate carboxylase. The comparison of both structures revealed that Arg884 supports tight inhibitor binding in the C3-type enzyme. In the C4 phosphoenolpyruvate carboxylase isoform, this arginine is replaced by glycine. The substitution reduces inhibitor affinity and enables the enzyme to participate in the C4 photosynthesis pathway.

A possible role for C4 photosynthetic enzymes in tolerance of Zea mays to NaCl

Treatment of 14-day-old maize cultivars (Hybrid351 and Giza2) with 250 mM NaCl significantly reduced shoot fresh and dry weights and protein content during the subsequent 12 days. The magnitude of reduction was more pronounced in Giza than Hybrid. Both cultivars contained converging levels of protein for the enzymes phosphoenolpyruvate carboxylase (PEPC), malate dehydrogenase (MDH), pyruvate phosphate dikinase (PPDK) and ribulose 1,5-bisphosphate carboxylase/oxygenase (Rubisco) under normal conditions; however, NaCl led to increase these levels in Hybrid and decrease them in Giza. Moreover, NaCl significantly inhibited the activities of PEPC, MDH and PPDK in both cultivars during the first 2 days, thereafter the inhibition nullified only in Hybrid; nonetheless, Rubisco was the least affected enzyme in both cultivars. In addition, NaCl slightly increased V (max) of PEPC, MDH and PPDK in Hybrid with no change in K (m); nevertheless V (max) dropped in Giza with an increase in K (m) of only PEPC and MDH. Also K (cat), K (cat)/K (m) and V (max)/K (m) of all enzymes were lower in treated Giza than in treated Hybrid. The increased V (max) of all enzymes in only Hybrid by NaCl confirms that they were synthesised more in Hybrid than in Giza. However, the decreased V (max) in Giza concomitant with the increased K (m) points to an interference of salinity with synthesis of enzymes and their structural integrity. This would lead to a noncompetitive inhibition for the enzymes. These findings declare that maize tolerance to NaCl was larger in Hybrid compared to Giza due to a role for C4 enzymes.

Molecular evolution of genes recruited into C₄ photosynthesis

The C₄ pathway is found in 62 lineages of land plants. We assess evidence for parallel versus convergent evolution of C₄ photosynthesis from three approaches: (i) studies of specific genes and cis-elements controlling their expression; (ii) phylogenetic analyses of mRNAs and inferred amino acid sequences; and (iii) analysis of C₃ and C₄ genomes and transcriptomes. Evidence suggests that although convergent evolution is common, parallel evolution can underlie both changes to gene expression and amino acid sequence. cis-elements that direct cell specificity in C₄ leaves are present in C₃ orthologues of genes recruited into C₄, probably facilitating this parallel evolution. From this, and genomic data, we propose that gene duplication followed by neofunctionalisation is not necessarily important in the evolution of C₄ biochemistry.

Multiple Arabidopsis genes primed for recruitment into C₄ photosynthesis

C(4) photosynthesis occurs in the most productive crops and vegetation on the planet, and has become widespread because it allows increased rates of photosynthesis compared with the ancestral C(3) pathway. Leaves of C(4) plants typically possess complicated alterations to photosynthesis, such that its reactions are compartmented between mesophyll and bundle sheath cells. Despite its complexity, the C(4) pathway has arisen independently in 62 separate lineages of land plants, and so represents one of the most striking examples of convergent evolution known. We demonstrate that elements in untranslated regions (UTRs) of multiple genes important for C(4) photosynthesis contribute to the metabolic compartmentalization characteristic of a C(4) leaf. Either the 5' or the 3' UTR is sufficient for cell specificity, indicating that functional redundancy underlies this key aspect of C(4) gene expression. Furthermore, we show that orthologous PPDK and CA genes from the C(3) plant Arabidopsis thaliana are primed for recruitment into the C(4) pathway. Elements sufficient for M-cell specificity in C(4) leaves are also present in both the 5' and 3' UTRs of these C(3) A. thaliana genes. These data indicate functional latency within the UTRs of genes from C(3) species that have been recruited into the C(4) pathway. The repeated recruitment of pre-existing cis-elements in C(3) genes may have facilitated the evolution of C(4) photosynthesis. These data also highlight the importance of alterations in trans in producing a functional C(4) leaf, and so provide insight into both the evolution and molecular basis of this important type of photosynthesis.

Regulation of the photorespiratory GLDPA gene in C(4) flaveria: an intricate interplay of transcriptional and posttranscriptional processes

The mitochondrial Gly decarboxylase complex (GDC) is a key component of the photorespiratory pathway that occurs in all photosynthetically active tissues of C(3) plants but is restricted to bundle sheath cells in C(4) species. GDC is also required for general cellular C(1) metabolism. In the Asteracean C(4) species Flaveria trinervia, a single functional GLDP gene, GLDPA, encodes the P-subunit of GDC, a decarboxylating Gly dehydrogenase. GLDPA promoter reporter gene fusion studies revealed that this promoter is active in bundle sheath cells and the vasculature of transgenic Flaveria bidentis (C(4)) and the Brassicacean C(3) species Arabidopsis thaliana, suggesting the existence of an evolutionarily conserved gene regulatory system in the bundle sheath. Here, we demonstrate that GLDPA gene regulation is achieved by an intricate interplay of transcriptional and posttranscriptional mechanisms. The GLDPA promoter is composed of two tandem promoters, P(R2) and P(R7), that together ensure a strong bundle sheath expression. While the proximal promoter (P(R7)) is active in the bundle sheath and vasculature, the distal promoter (P(R2)) drives uniform expression in all leaf chlorenchyma cells and the vasculature. An intron in the 5' untranslated leader of P(R2)-derived transcripts is inefficiently spliced and apparently suppresses the output of P(R2) by eliciting RNA decay.

The remarkable diversity of plant PEPC (phosphoenolpyruvate carboxylase): recent insights into the physiological functions and post-translational controls of non-photosynthetic PEPCs

PEPC [PEP (phosphoenolpyruvate) carboxylase] is a tightly controlled enzyme located at the core of plant C-metabolism that catalyses the irreversible β-carboxylation of PEP to form oxaloacetate and Pi. The critical role of PEPC in assimilating atmospheric CO2 during C4 and Crassulacean acid metabolism photosynthesis has been studied extensively. PEPC also fulfils a broad spectrum of non-photosynthetic functions, particularly the anaplerotic replenishment of tricarboxylic acid cycle intermediates consumed during biosynthesis and nitrogen assimilation. An impressive array of strategies has evolved to co-ordinate in vivo PEPC activity with cellular demands for C4–C6 carboxylic acids. To achieve its diverse roles and complex regulation, PEPC belongs to a small multigene family encoding several closely related PTPCs (plant-type PEPCs), along with a distantly related BTPC (bacterial-type PEPC). PTPC genes encode ~110-kDa polypeptides containing conserved serine-phosphorylation and lysine-mono-ubiquitination sites, and typically exist as homotetrameric Class-1 PEPCs. In contrast, BTPC genes encode larger ~117-kDa polypeptides owing to a unique intrinsically disordered domain that mediates BTPC's tight interaction with co-expressed PTPC subunits. This association results in the formation of unusual ~900-kDa Class-2 PEPC hetero-octameric complexes that are desensitized to allosteric effectors. BTPC is a catalytic and regulatory subunit of Class-2 PEPC that is subject to multi-site regulatory phosphorylation in vivo. The interaction between divergent PEPC polypeptides within Class-2 PEPCs adds another layer of complexity to the evolution, physiological functions and metabolic control of this essential CO2-fixing plant enzyme. The present review summarizes exciting developments concerning the functions, post-translational controls and subcellular location of plant PTPC and BTPC isoenzymes.

Structural and metabolic transitions of C4 leaf development and differentiation defined by microscopy and quantitative proteomics in maize

C(4) grasses, such as maize (Zea mays), have high photosynthetic efficiency through combined biochemical and structural adaptations. C(4) photosynthesis is established along the developmental axis of the leaf blade, leading from an undifferentiated leaf base just above the ligule into highly specialized mesophyll cells (MCs) and bundle sheath cells (BSCs) at the tip. To resolve the kinetics of maize leaf development and C(4) differentiation and to obtain a systems-level understanding of maize leaf formation, the accumulation profiles of proteomes of the leaf and the isolated BSCs with their vascular bundle along the developmental gradient were determined using large-scale mass spectrometry. This was complemented by extensive qualitative and quantitative microscopy analysis of structural features (e.g., Kranz anatomy, plasmodesmata, cell wall, and organelles). More than 4300 proteins were identified and functionally annotated. Developmental protein accumulation profiles and hierarchical cluster analysis then determined the kinetics of organelle biogenesis, formation of cellular structures, metabolism, and coexpression patterns. Two main expression clusters were observed, each divided in subclusters, suggesting that a limited number of developmental regulatory networks organize concerted protein accumulation along the leaf gradient. The coexpression with BSC and MC markers provided strong candidates for further analysis of C(4) specialization, in particular transporters and biogenesis factors. Based on the integrated information, we describe five developmental transitions that provide a conceptual and practical template for further analysis. An online protein expression viewer is provided through the Plant Proteome Database.