Kranz and single-cell forms of C4 plants in the subfamily Suaedoideae show kinetic C4 convergence for PEPC and Rubisco with divergent amino acid substitutions

Multiple highly expressed phosphoenolpyruvate carboxylase genes have divergent enzyme kinetic properties in two C4 grasses

BACKGROUND AND AIMS

Phosphoenolpyruvate (PEP) carboxylase (PEPC) catalyses the irreversible carboxylation of PEP with bicarbonate to produce oxaloacetate. This reaction powers the carbon-concentrating mechanism (CCM) in plants that perform C4 photosynthesis. This CCM is generally driven by a single PEPC gene product that is highly expressed in the cytosol of mesophyll cells. We found two C4 grasses, Panicum miliaceum and Echinochloa colona, that each have two highly expressed PEPC genes. We characterized the kinetic properties of the two most abundant PEPCs in E. colona and P. miliaceum to better understand how the enzyme's amino acid structure influences its function.

METHODS

Coding sequences of the two most abundant PEPC proteins in E. colona and P. miliaceum were synthesized by GenScript and were inserted into bacteria expression plasmids. Point mutations resulting in substitutions at conserved amino acid residues (e.g. N-terminal serine and residue 890) were created via site-directed PCR mutagenesis. The kinetic properties of semi-purified plant PEPCs from Escherichia coli were analysed using membrane-inlet mass spectrometry and a spectrophotometric enzyme-coupled reaction.

KEY RESULTS

The two most abundant P. miliaceum PEPCs (PmPPC1 and PmPPC2) have similar sequence identities (>95 %), and as a result had similar kinetic properties. The two most abundant E. colona PEPCs (EcPPC1 and EcPPC2) had identities of ~78 % and had significantly different kinetic properties. The PmPPCs and EcPPCs had different responses to allosteric inhibitors and activators, and substitutions at the conserved N-terminal serine and residue 890 resulted in significantly altered responses to allosteric regulators.

CONCLUSIONS

The two, significantly expressed C4Ppc genes in P. miliaceum were probably the result of genomes combining from two closely related C4Panicum species. We found natural variation in PEPC's sensitivity to allosteric inhibition that seems to bypass the conserved 890 residue, suggesting alternative evolutionary pathways for increased malate tolerance and other kinetic properties.

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.

Lateral Gene Transfer Acts As an Evolutionary Shortcut to Efficient C4 Biochemistry

The adaptation of proteins for novel functions often requires changes in their kinetics via amino acid replacement. This process can require multiple mutations, and therefore extended periods of selection. The transfer of genes among distinct species might speed up the process, by providing proteins already adapted for the novel function. However, this hypothesis remains untested in multicellular eukaryotes. The grass Alloteropsis is an ideal system to test this hypothesis due to its diversity of genes encoding phosphoenolpyruvate carboxylase, an enzyme that catalyzes one of the key reactions in the C₄ pathway. Different accessions of Alloteropsis either use native isoforms relatively recently co-opted from other functions or isoforms that were laterally acquired from distantly related species that evolved the C₄ trait much earlier. By comparing the enzyme kinetics, we show that native isoforms with few amino acid replacements have substrate K_M values similar to the non-C₄ ancestral form, but exhibit marked increases in catalytic efficiency. The co-option of native isoforms was therefore followed by rapid catalytic improvements, which appear to rely on standing genetic variation observed within one species. Native C₄ isoforms with more amino acid replacements exhibit additional changes in affinities, suggesting that the initial catalytic improvements are followed by gradual modifications. Finally, laterally acquired genes show both strong increases in catalytic efficiency and important changes in substrate handling. We conclude that the transfer of genes among distant species sharing the same physiological novelty creates an evolutionary shortcut toward more efficient enzymes, effectively accelerating evolution.

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.

Evolutionary Convergence of C4 Photosynthesis: A Case Study in the Nyctaginaceae

C4 photosynthesis evolved over 65 times, with around 24 origins in the eudicot order Caryophyllales. In the Caryophyllales family Nyctaginaceae, the C4 pathway is known in three genera of the tribe Nyctagineae: Allionia, Okenia and Boerhavia. Phylogenetically, Allionia and Boerhavia/Okenia are separated by three genera whose photosynthetic pathway is uncertain. To clarify the distribution of photosynthetic pathways in the Nyctaginaceae, we surveyed carbon isotope ratios of 159 species of the Nyctaginaceae, along with bundle sheath (BS) cell ultrastructure, leaf gas exchange, and C4 pathway biochemistry in five species from the two C4 clades and closely related C3 genera. All species in Allionia, Okenia and Boerhavia are C4, while no C4 species occur in any other genera of the family, including three that branch between Allionia and Boerhavia. This demonstrates that C4 photosynthesis evolved twice in Nyctaginaceae. Boerhavia species use the NADP-malic enzyme (NADP-ME) subtype of C4 photosynthesis, while Allionia species use the NAD-malic enzyme (NAD-ME) subtype. The BS cells of Allionia have many more mitochondria than the BS of Boerhavia. Bundle sheath mitochondria are closely associated with chloroplasts in Allionia which facilitates CO2 refixation following decarboxylation by mitochondrial NAD-ME. The close relationship between Allionia and Boerhavia could provide insights into why NADP-ME versus NAD-ME subtypes evolve, particularly when coupled to analysis of their respective genomes. As such, the group is an excellent system to dissect the organizational hierarchy of convergent versus divergent traits produced by C4 evolution, enabling us to understand when convergence is favored versus when divergent modifications can result in a common phenotype.

Kinetic Modifications of C4 PEPC Are Qualitatively Convergent, but Larger in Panicum Than in Flaveria

C4 photosynthesis results from a set of anatomical features and biochemical components that act together to concentrate CO2 within the leaf and boost productivity. This complex trait evolved independently many times, resulting in various realizations of the phenotype, but in all C4 plants the primary fixation of atmospheric carbon is catalyzed by phosphoenolpyruvate carboxylase. Comparisons of C4 and non-C4 PEPC from a few closely related species suggested that the enzyme was modified to meet the demands of the C4 cycle. However, very few C4 groups have been investigated, hampering general conclusions. To test the hypothesis that distant C4 lineages underwent convergent biochemical changes, we compare the kinetic variation between C4 and non-C4 PEPC from a previously assessed young lineage (Flaveria, Asteraceae) with those from an older lineage found within the distantly related grass family (Panicum). Despite the evolutionary distance, the kinetic changes between the non-C4 and C4 PEPC are qualitatively similar, with a decrease in sensitivity for inhibitors, an increased specificity (k cat/K m) for bicarbonate, and a decreased specificity (k cat/K m) for PEP. The differences are more pronounced in the older lineage Panicum, which might indicate that optimization of PEPC for the C4 context increases with evolutionary time.

Exploring molecular evolution of Rubisco in C3 and CAM Orchidaceae and Bromeliaceae

BACKGROUND

The CO2-concentrating mechanism associated to Crassulacean acid metabolism (CAM) alters the catalytic context for Rubisco by increasing CO2 availability and provides an advantage in particular ecological conditions. We hypothesized about the existence of molecular changes linked to these particular adaptations in CAM Rubisco. We investigated molecular evolution of the Rubisco large (L-) subunit in 78 orchids and 144 bromeliads with C3 and CAM photosynthetic pathways. The sequence analyses were complemented with measurements of Rubisco kinetics in some species with contrasting photosynthetic mechanism and differing in the L-subunit sequence.

RESULTS

We identified potential positively selected sites and residues with signatures of co-adaptation. The implementation of a decision tree model related Rubisco specific variable sites to the leaf carbon isotopic composition of the species. Differences in the Rubisco catalytic traits found among C3 orchids and between strong CAM and C3 bromeliads suggested Rubisco had evolved in response to differing CO2 concentration.

CONCLUSIONS

The results revealed that the variability in the Rubisco L-subunit sequence in orchids and bromeliads is composed of coevolving sites under potential positive adaptive signal. The sequence variability was related to δ13C in orchids and bromeliads, however it could not be linked to the variability found in the kinetic properties of the studied species.

Cloning, Expression and Characterization of Two Beta Carbonic Anhydrases from a Newly Isolated CO2 Fixer, Serratia marcescens Wy064

Bacterial strains from karst landform soil were enriched via chemostat culture in the presence of sodium bicarbonate. Two chemolithotrophic strains were isolated and identified as Serratia marcescens Wy064 and Bacillus sp. Wy065. Both strains could grow using sodium bicarbonate as the sole carbon source. Furthermore, the supplement of the medium with three electron donors (Na₂S, NaNO₂, and Na₂S₂O₃) improved the growth of both strains. The activities of carbonic anhydrase (CA) and ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) could be detected in the crude enzyme of strain Wy064, implying that the strain Wy064 might employ Calvin cycle to fix CO₂. S. marcescens genome mining revealed four potential CA genes designated CA1-CA4. The proteins encoded by genes CA1-3 were cloned and expressed in Escherichia coli. The purified recombinant enzymes of CA1 and CA3 exhibited CO₂ hydration activities, whereas enzyme CA2 was expressed in inclusion bodies. A CO₂ hydration assay demonstrated that the specific activity of CA3 was significantly higher than that of CA1. The maximum CO₂ hydration activities for CA1 and CA3 were observed at pH 7.5 and 40 °C. The activities of CA1 and CA3 were significantly enhanced by several metal ions, especially Zn²⁺, which resulted in 21.1-fold and 26.1-fold increases of CO₂ hydration activities, respectively. The apparent K _m and V _max for CO₂ as substrate were 27 mM and 179 WAU/mg for CA1, and 14 mM and 247 WAU/mg for CA3, respectively. Structure modeling combined with sequence analysis indicated that CA1 and CA3 should belong to the Type II β-CA.

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.

Efficient In Vivo Screening Method for the Identification of C4 Photosynthesis Inhibitors Based on Cell Suspensions of the Single-Cell C4 Plant Bienertia sinuspersici

The identification of novel herbicides is of crucial importance to modern agriculture. We developed an efficient in vivo assay based on oxygen evolution measurements using suspensions of chlorenchyma cells isolated from the single-cell C₄ plant Bienertia sinuspersici to identify and characterize inhibitors of C₄ photosynthesis. This novel approach fills the gap between conventional in vitro assays for inhibitors targeting C₄ key enzymes and in vivo experiments on whole plants. The assay addresses inhibition of the target enzymes in a plant context thereby taking care of any reduced target inhibition due to metabolization or inadequate uptake of small molecule inhibitors across plant cell walls and membranes. Known small molecule inhibitors targeting C₄ photosynthesis were used to validate the approach. To this end, we tested pyruvate phosphate dikinase inhibitor bisindolylmaleimide IV and phosphoenolpyruvate carboxylase inhibitor okanin. Both inhibitors show inhibition of plant photosynthesis at half-maximal inhibitory concentrations in the sub-mM range and confirm their potential to act as a new class of C₄ selective inhibitors.

Identification of the allosteric site for neutral amino acids in the maize C4 isozyme of phosphoenolpyruvate carboxylase: The critical role of Ser-100

The isozymes of photosynthetic phosphoenolpyruvate carboxylase from C4 plants (PEPC-C4) play a critical role in their atmospheric CO2 assimilation and productivity. They are allosterically activated by phosphorylated trioses or hexoses, such as d-glucose 6-phosphate, and inhibited by l-malate or l-aspartate. Additionally, PEPC-C4 isozymes from grasses are activated by glycine, serine, or alanine, but the allosteric site for these compounds remains unknown. Here, we report a new crystal structure of the isozyme from Zea mays (ZmPEPC-C4) with glycine bound at the monomer-monomer interfaces of the two dimers of the tetramer, making interactions with residues of both monomers. This binding site is close to, but different from, the one proposed to bind glucose 6-phosphate. Docking experiments indicated that d/l-serine or d/l-alanine could also bind to this site, which does not exist in the PEPC-C4 isozyme from the eudicot plant Flaveria, mainly because of a lysyl residue at the equivalent position of Ser-100 in ZmPEPC-C4 Accordingly, the ZmPEPC-C4 S100K mutant is not activated by glycine, serine, or alanine. Amino acid sequence alignments showed that PEPC-C4 isozymes from the monocot family Poaceae have either serine or glycine at this position, whereas those from Cyperaceae and eudicot families have lysine. The size and charge of the residue equivalent to Ser-100 are not only crucial for the activation of PEPC-C4 isozymes by neutral amino acids but also affect their affinity for the substrate phosphoenolpyruvate and their allosteric regulation by glucose 6-phosphate and malate, accounting for the reported kinetic differences between PEPC-C4 isozymes from monocot and eudicot plants.

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.

Surveying Rubisco Diversity and Temperature Response to Improve Crop Photosynthetic Efficiency

The threat to global food security of stagnating yields and population growth makes increasing crop productivity a critical goal over the coming decades. One key target for improving crop productivity and yields is increasing the efficiency of photosynthesis. Central to photosynthesis is Rubisco, which is a critical but often rate-limiting component. Here, we present full Rubisco catalytic properties measured at three temperatures for 75 plants species representing both crops and undomesticated plants from diverse climates. Some newly characterized Rubiscos were naturally "better" compared to crop enzymes and have the potential to improve crop photosynthetic efficiency. The temperature response of the various catalytic parameters was largely consistent across the diverse range of species, though absolute values showed significant variation in Rubisco catalysis, even between closely related species. An analysis of residue differences among the species characterized identified a number of candidate amino acid substitutions that will aid in advancing engineering of improved Rubisco in crop systems. This study provides new insights on the range of Rubisco catalysis and temperature response present in nature, and provides new information to include in models from leaf to canopy and ecosystem scale.

Rubisco Catalytic Properties and Temperature Response in Crops

Rubisco catalytic traits and their thermal dependence are two major factors limiting the CO2 assimilation potential of plants. In this study, we present the profile of Rubisco kinetics for 20 crop species at three different temperatures. The results largely confirmed the existence of significant variation in the Rubisco kinetics among species. Although some of the species tended to present Rubisco with higher thermal sensitivity (e.g. Oryza sativa) than others (e.g. Lactuca sativa), interspecific differences depended on the kinetic parameter. Comparing the temperature response of the different kinetic parameters, the Rubisco Km for CO2 presented higher energy of activation than the maximum carboxylation rate and the CO2 compensation point in the absence of mitochondrial respiration. The analysis of the Rubisco large subunit sequence revealed the existence of some sites under adaptive evolution in branches with specific kinetic traits. Because Rubisco kinetics and their temperature dependency were species specific, they largely affected the assimilation potential of Rubisco from the different crops, especially under those conditions (i.e. low CO2 availability at the site of carboxylation and high temperature) inducing Rubisco-limited photosynthesis. As an example, at 25°C, Rubisco from Hordeum vulgare and Glycine max presented, respectively, the highest and lowest potential for CO2 assimilation at both high and low chloroplastic CO2 concentrations. In our opinion, this information is relevant to improve photosynthesis models and should be considered in future attempts to design more efficient Rubiscos.