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

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.

A Machine Learning Framework Identifies Plastid-Encoded Proteins Harboring C3 and C4 Distinguishing Sequence Information

C4 photosynthesis is known to have at least 61 independent origins across plant lineages making it one of the most notable examples of convergent evolution. Of the >60 independent origins, a predicted 22-24 origins, encompassing greater than 50% of all known C4 species, exist within the Panicoideae, Arundinoideae, Chloridoideae, Micrairoideae, Aristidoideae, and Danthonioideae (PACMAD) clade of the Poaceae family. This clade is therefore primed with species ideal for the study of genomic changes associated with the acquisition of the C4 photosynthetic trait. In this study, we take advantage of the growing availability of sequenced plastid genomes and employ a machine learning (ML) approach to screen for plastid genes harboring C3 and C4 distinguishing information in PACMAD species. We demonstrate that certain plastid-encoded protein sequences possess distinguishing and informative sequence information that allows them to train accurate ML C3/C4 classification models. Our RbcL-trained model, for example, informs a C3/C4 classifier with greater than 99% accuracy. Accurate prediction of photosynthetic type from individual sequences suggests biologically relevant, and potentially differing roles of these sequence products in C3 versus C4 metabolism. With this ML framework, we have identified several key sequences and sites that are most predictive of C3/C4 status, including RbcL, subunits of the NAD(P)H dehydrogenase complex, and specific residues within, further highlighting their potential significance in the evolution and/or maintenance of C4 photosynthetic machinery. This general approach can be applied to uncover intricate associations between other similar genotype-phenotype relationships.

Role of C4 photosynthetic enzyme isoforms in C3 plants and their potential applications in improving agronomic traits in crops

As compared to C3, C4 plants have higher photosynthetic rates and better tolerance to high temperature and drought. These traits are highly beneficial in the current scenario of global warming. Interestingly, all the genes of the C4 photosynthetic pathway are present in C3 plants, although they are involved in diverse non-photosynthetic functions. Non-photosynthetic isoforms of carbonic anhydrase (CA), phosphoenolpyruvate carboxylase (PEPC), malate dehydrogenase (MDH), the decarboxylating enzymes NAD/NADP-malic enzyme (NAD/NADP-ME), and phosphoenolpyruvate carboxykinase (PEPCK), and finally pyruvate orthophosphate dikinase (PPDK) catalyze reactions that are essential for major plant metabolism pathways, such as the tricarboxylic acid (TCA) cycle, maintenance of cellular pH, uptake of nutrients and their assimilation. Consistent with this view differential expression pattern of these non-photosynthetic C3 isoforms has been observed in different tissues across the plant developmental stages, such as germination, grain filling, and leaf senescence. Also abundance of these C3 isoforms is increased considerably in response to environmental fluctuations particularly during abiotic stress. Here we review the vital roles played by C3 isoforms of C4 enzymes and the probable mechanisms by which they help plants in acclimation to adverse growth conditions. Further, their potential applications to increase the agronomic trait value of C3 crops is discussed.

Differential timing of gene expression and recruitment in independent origins of CAM in the Agavoideae (Asparagaceae)

Crassulacean acid metabolism (CAM) photosynthesis has evolved repeatedly across the plant tree of life, however our understanding of the genetic convergence across independent origins remains hampered by the lack of comparative studies. Here, we explore gene expression profiles in eight species from the Agavoideae (Asparagaceae) encompassing three independent origins of CAM. Using comparative physiology and transcriptomics, we examined the variable modes of CAM in this subfamily and the changes in gene expression across time of day and between well watered and drought-stressed treatments. We further assessed gene expression and the molecular evolution of genes encoding phosphoenolpyruvate carboxylase (PPC), an enzyme required for primary carbon fixation in CAM. Most time-of-day expression profiles are largely conserved across all eight species and suggest that large perturbations to the central clock are not required for CAM evolution. By contrast, transcriptional response to drought is highly lineage specific. Yucca and Beschorneria have CAM-like expression of PPC2, a copy of PPC that has never been shown to be recruited for CAM in angiosperms. Together the physiological and transcriptomic comparison of closely related C₃ and CAM species reveals similar gene expression profiles, with the notable exception of differential recruitment of carboxylase enzymes for CAM function.

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.

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.

Molecular Cloning of Novel-Type Phosphoenolpyruvate Carboxylase Isoforms in Pitaya (Hylocereus undatus)

Phosphoenolpyruvate carboxylase (PEPC) is an important enzyme involved in the initial CO₂ fixation of crassulacean acid metabolism (CAM) photosynthesis. To understand the cultivation characteristics of a CAM plant pitaya, it is necessary to clarify the characteristics of PEPC in this species. Here, we cloned three PEPC cDNAs in pitaya, HuPPC1, HuPPC2, and HuPPC3, which encode 942, 934, and 966 amino acid residues, respectively. Phylogenetic analysis indicated that these PEPC belonged to plant-type PEPC (PTPC), although HuPPC1 and HuPPC2 have no Ser-phosphorylation motif in N-terminal region, which is a crucial regulation site in PTPC and contributes to CAM periodicity. HuPPC1 and HuPPC2 phylogenetically unique to the Cactaceae family, whereas HuPPC3 was included in a CAM clade. Two isoforms were partially purified at the protein level and were assigned as HuPPC2 and HuPPC3 using MASCOT analysis. The most distinct difference in enzymatic properties between the two was sensitivity to malate and aspartate, both of which are allosteric inhibitors of PEPC. With 2 mM malate, HuPPC3 was inhibited to 10% of the initial activity, whereas HuPPC2 activity was maintained at 70%. Aspartate inhibited HuPPC3 activity by approximately 50% at 5 mM; however, such inhibition was not observed for HuPPC2 at 10 mM. These results suggest that HuPPC3 corresponds to a general CAM-related PEPC, whereas HuPPC1 and HuPPC2 are related to carbon and/or nitrogen metabolism, with a characteristic regulation mechanism similar to those of Cactaceae plants.

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.

Identification and evolution of C4 photosynthetic pathway genes in plants

BACKGROUND

NADP-malic enzyme (NAPD-ME), and pyruvate orthophosphate dikinase (PPDK) are important enzymes that participate in C₄ photosynthesis. However, the evolutionary history and forces driving evolution of these genes in C₄ plants are not completely understood.

RESULTS

We identified 162 NADP-ME and 35 PPDK genes in 25 species and constructed respective phylogenetic trees. We classified NADP-ME genes into four branches, A1, A2, B1 and B2, whereas PPDK was classified into two branches in which monocots were in branch I and dicots were in branch II. Analyses of selective pressure on the NAPD-ME and PPDK gene families identified four positively selected sites, including 94H and 196H in the a5 branch of NADP-ME, and 95A and 559E in the e branch of PPDK at posterior probability thresholds of 95%. The positively selected sites were located in the helix and sheet regions. Quantitative RT-PCR (qRT-PCR) analyses revealed that expression levels of 6 NADP-ME and 2 PPDK genes from foxtail millet were up-regulated after exposure to light.

CONCLUSION

This study revealed that positively selected sites of NADP-ME and PPDK evolution in C₄ plants. It provides information on the classification and positive selection of plant NADP-ME and PPDK genes, and the results should be useful in further research on the evolutionary history of C₄ plants.

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.

Molecular evolution of key metabolic genes during transitions to C4 and CAM photosynthesis

PREMISE OF THE STUDY

Next-generation sequencing facilitates rapid production of well-sampled phylogenies built from very large genetic data sets, which can then be subsequently exploited to examine the molecular evolution of the genes themselves. We present an evolutionary analysis of 83 gene families (19 containing carbon-concentrating mechanism (CCM) genes, 64 containing non-CCM genes) in the portullugo clade (Caryophyllales), a diverse lineage of mostly arid-adapted plants that contains multiple evolutionary origins of all known photosynthesis types in land plants (C₃ , C₄ , CAM, C₄ -CAM, and various intermediates).

METHODS

We inferred a phylogeny of 197 individuals from 167 taxa using coalescent-based approaches and individual gene family trees using maximum likelihood. Positive selection analyses were conducted on individual gene family trees with a mixed effects model of evolution (MEME). We devised new indices to compare levels of convergence and prevalence of particular residues between CCM and non-CCM genes and between species with different photosynthetic pathways.

KEY RESULTS

Contrary to expectations, there were no significant differences in the levels of positive selection detected in CCM versus non-CCM genes. However, we documented a significantly higher level of convergent amino acid substitutions in CCM genes, especially in C₄ taxa.

CONCLUSIONS

Our analyses reveal a new suite of amino acid residues putatively important for C₄ and CAM function. We discuss both the advantages and challenges of using targeted enrichment sequence data for exploratory studies of molecular evolution.

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.

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

The two carboxylation reactions performed by phosphoenolpyruvate carboxylase (PEPC) and ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) are vital in the fixation of inorganic carbon for C4 plants. The abundance of PEPC is substantially elevated in C4 leaves, while the location of Rubisco is restricted to one of two chloroplast types. These differences compared with C3 leaves have been shown to result in convergent enzyme optimization in some C4 species. Investigation into the kinetic properties of PEPC and Rubisco from Kranz C4, single cell C4, and C3 species in Chenopodiaceae s. s. subfamily Suaedoideae showed that these major carboxylases in C4 Suaedoideae species lack the same mutations found in other C4 systems which have been examined; but still have similar convergent kinetic properties. Positive selection analysis on the N-terminus of PEPC identified residues 364 and 368 to be under positive selection with a posterior probability >0.99 using Bayes empirical Bayes. Compared with previous analyses on other C4 species, PEPC from C4 Suaedoideae species have different convergent amino acids that result in a higher K m for PEP and malate tolerance compared with C3 species. Kinetic analysis of Rubisco showed that C4 species have a higher catalytic efficiency of Rubisco (k catc in mol CO2 mol(-1) Rubisco active sites s(-1)), despite lacking convergent substitutions in the rbcL gene. The importance of kinetic changes to the two-carboxylation reactions in C4 leaves related to amino acid selection is discussed.

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.