The regulation of gene expression required for C4 photosynthesis

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.

Engineering carbon assimilation in plants

Carbon assimilation is a crucial part of the photosynthetic process, wherein inorganic carbon, typically in the form of CO2, is converted into organic compounds by living organisms, including plants, algae, and a subset of bacteria. Although several carbon fixation pathways have been elucidated, the Calvin-Benson-Bassham (CBB) cycle remains fundamental to carbon metabolism, playing a pivotal role in the biosynthesis of starch and sucrose in plants, algae, and cyanobacteria. However, Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO), the key carboxylase enzyme of the CBB cycle, exhibits low kinetic efficiency, low substrate specificity, and high temperature sensitivity, all of which have the potential to limit flux through this pathway. Consequently, RuBisCO needs to be present at very high concentrations, which is one of the factors contributing to its status as the most prevalent protein on Earth. Numerous attempts have been made to optimize the catalytic efficiency of RuBisCO and thereby promote plant growth. Furthermore, the limitations of this process highlight the potential benefits of engineering or discovering more efficient carbon fixation mechanisms, either by improving RuBisCO itself or by introducing alternative pathways. Here, we review advances in artificial carbon assimilation engineering, including the integration of synthetic biology, genetic engineering, metabolic pathway optimization, and artificial intelligence in order to create plants capable of performing more efficient photosynthesis. We additionally provide a perspective of current challenges and potential solutions alongside a personal opinion of the most promising future directions of this emerging field.

Prospects for synthetic biology in 21st Century agriculture

Plant synthetic biology has emerged as a transformative field in agriculture, offering innovative solutions to enhance food security, provide resilience to climate change, and transition to sustainable farming practices. By integrating advanced genetic tools, computational modeling, and systems biology, researchers can precisely modify plant genomes to enhance traits such as yield, stress tolerance, and nutrient use efficiency. The ability to design plants with specific characteristics tailored to diverse environmental conditions and agricultural needs holds great potential to address global food security challenges. Here we highlight recent advancements and applications of plant synthetic biology in agriculture, focusing on key areas such as photosynthetic efficiency, nitrogen fixation, drought tolerance, pathogen resistance, nutrient use efficiency, biofortification, climate resilience, microbiology engineering, synthetic plant genomes, and the integration of artificial intelligence (AI) with synthetic biology. These innovations aim to maximize resource use efficiency, reduce reliance on external inputs, and mitigate environmental impacts associated with conventional agricultural practices. Despite challenges related to regulatory approval and public acceptance, the integration of synthetic biology in agriculture holds immense promise for creating more resilient and sustainable agricultural systems, contributing to global food security and environmental sustainability. Rigorous multi-field testing of these approaches will undoubtedly be required to ensure reproducibility.

Leaf transcriptomes from C3, C3-C4 intermediate, and C4Neurachne species give insights into C4 photosynthesis evolution

The C4 photosynthetic pathway is hypothesized to have evolved from the ancestral C3 pathway through progressive changes in leaf anatomy and biochemistry with extant C3-C4 photosynthetic intermediate species representing phenotypes between species demonstrating full C3 and full C4 states. The Australian endemic genus Neurachne is the only known grass group that contains distinct, closely related species that carry out C3, C3-C4 intermediate, or C4 photosynthesis. To explore and understand the molecular mechanisms underlying C4 photosynthesis evolution in this genus, leaf transcriptomes were generated from two C3, three photosynthetic intermediate (proto-Kranz, C2-like, and C2), and two C4Neurachne species. The data were used to reconstruct phylogenetic relationships in Neurachne, which confirmed two independent C4 origins in the genus. Relative transcript abundances substantiated the photosynthetic phenotypes of individual species and highlighted transcriptional investment differences between species, including between the two C4 species. The data also revealed proteins potentially involved in C4 cycle intermediate transport and identified molecular mechanisms responsible for the evolution of C4-associated proteins in the genus.

Exaptation of ancestral cell-identity networks enables C4 photosynthesis

C4 photosynthesis is used by the most productive plants on the planet, and compared with the ancestral C3 pathway, it confers a 50% increase in efficiency1. In more than 60 C4 lineages, CO2 fixation is compartmentalized between tissues, and bundle-sheath cells become photosynthetically activated2. How the bundle sheath acquires this alternate identity that allows efficient photosynthesis is unclear. Here we show that changes to bundle-sheath gene expression in C4 leaves are associated with the gain of a pre-existing cis-code found in the C3 leaf. From single-nucleus gene-expression and chromatin-accessibility atlases, we uncover DNA binding with one finger (DOF) motifs that define bundle-sheath identity in the major crops C3 rice and C4 sorghum. Photosynthesis genes that are rewired to be strongly expressed in the bundle-sheath cells of C4 sorghum acquire cis-elements that are recognized by DOFs. Our findings are consistent with a simple model in which C4 photosynthesis is based on the recruitment of an ancestral cis-code associated with bundle-sheath identity. Gain of such elements harnessed a stable patterning of transcription factors between cell types that are found in both C3 and C4 leaves to activate photosynthesis in the bundle sheath. Our findings provide molecular insights into the evolution of the complex C4 pathway, and might also guide the rational engineering of C4 photosynthesis in C3 crops to improve crop productivity and resilience3,4.

Genome-wide comparative analysis of photosynthetic enzymatic genes provides novel insights into foxtail millet and other cereals

C4 crops have more efficient photosynthetic pathways that enable their higher photosynthetic capacities as well as nitrogen and water use efficiencies than C3 crops. Previous research has demonstrated that the genomes of C3 species include and express every gene needed for the C4 photosynthesis pathway. However, very little is known about the dynamics and evolutionary history of such genetic evolution in C4 plants. In this study, the genes encoding five key photosynthetic pathway enzymes in the genomes of C3 (rice), C4 (maize, sorghum, and foxtail millet), and CAM (pineapple) crops were identified and compared systematically. The numbers of genes in these photosynthetic enzymes were highest in the C4 crops like sorghum and foxtail millet, while only eight genes were identified in the CAM plant. However, 16 genes were identified in the C3 crop rice. Furthermore, we performed physical, chemical, gene structure and, cis-element analyses to obtain complete insights into these key genes. Tissue-specific expressions showed that most of the photosynthetic genes are expressed in the leaf tissues. Comparisons of the expression characteristics confirmed that the expression patterns of non-photosynthetic gene copies were relatively conserved among the species, while the C4 gene copies in the C4 species acquired new tissue expression patterns during evolution. Additionally, multiple sequence features that could affect C4 gene expressions and subcellular localization were found in the coding and promoter regions. Our research also highlights the variations in how different genes have evolved within the C4 photosynthetic pathway, and we confirmed that specific high expressions in the leaves and right distribution within the cells were crucial for the development of the C4 photosynthetic abilities. The findings of this study are expected to aid in understanding the evolutionary process of the C4 photosynthetic pathway in grasses as well as offer insights for modifying the C4 photosynthetic pathways in wheat, rice, and other significant C3 cereal crops.

Identifying genomic regions associated with C4 photosynthetic activity and leaf anatomy in Alloteropsis semialata

C4 photosynthesis is a complex trait requiring multiple developmental and metabolic alterations. Despite this complexity, it has independently evolved over 60 times. However, our understanding of the transition to C4 is complicated by the fact that variation in photosynthetic type is usually segregated between species that diverged a long time ago. Here, we perform a genome-wide association study (GWAS) using the grass Alloteropsis semialata, the only known species to have C3, intermediate, and C4 accessions that recently diverged. We aimed to identify genomic regions associated with the strength of the C4 cycle (measured using δ13C), and the development of C4 leaf anatomy. Genomic regions correlated with δ13C include regulators of C4 decarboxylation enzymes (RIPK), nonphotochemical quenching (SOQ1), and the development of Kranz anatomy (SCARECROW-LIKE). Regions associated with the development of C4 leaf anatomy in the intermediate individuals contain additional leaf anatomy regulators, including those responsible for vein patterning (GSL8) and meristem determinacy (GIF1). The parallel recruitment of paralogous leaf anatomy regulators between A. semialata and other C4 lineages implies the co-option of these genes is context-dependent, which likely has implications for the engineering of the C4 trait into C3 species.

Single-cell RNA-sequencing provides new insights into the cell-specific expression patterns and transcriptional regulation of photosynthetic genes in bermudagrass leaf blades

As an important warm-season turfgrass species, bermudagrass (Cynodon dactylon L.) flourishes in warm areas around the world due to the existence of the C4 photosynthetic pathway. However, how C4 photosynthesis operates in bermudagrass leaves is still poorly understood. In this study, we performed single-cell RNA-sequencing on 5296 cells from bermudagrass leaf blades. Eight cell clusters corresponding to mesophyll, bundle sheath, epidermis and vascular bundle cells were successfully identified using known cell marker genes. Expression profiling indicated that genes encoding NADP-dependent malic enzymes (NADP-MEs) were highly expressed in bundle sheath cells, whereas NAD-ME genes were weakly expressed in all cell types, suggesting C4 photosynthesis of bermudagrass leaf blades might be NADP-ME type rather than NAD-ME type. The results also indicated that starch synthesis-related genes showed preferential expression in bundle sheath cells, whereas starch degradation-related genes were highly expressed in mesophyll cells, which agrees with the observed accumulation of starch-filled chloroplasts in bundle sheath cells. Gene co-expression analysis further revealed that different families of transcription factors were co-expressed with multiple C4 photosynthesis-related genes, suggesting a complex transcription regulatory network of C4 photosynthesis might exist in bermudagrass leaf blades. These findings collectively provided new insights into the cell-specific expression patterns and transcriptional regulation of photosynthetic genes in bermudagrass.

Transcriptome and small RNA analysis unveils novel insights into the C4 gene regulation in sugarcane

MAIN CONCLUSION

This study reveals miRNA indirect regulation of C4 genes in sugarcane through transcription factors, highlighting potential key regulators like SsHAM3a. C4 photosynthesis is crucial for the high productivity and biomass of sugarcane, however, the miRNA regulation of C4 genes in sugarcane remains elusive. We have identified 384 miRNAs along the leaf gradients, including 293 known miRNAs and 91 novel miRNAs. Among these, 86 unique miRNAs exhibited differential expression patterns, and we identified 3511 potential expressed targets of these differentially expressed miRNAs (DEmiRNAs). Analyses using Pearson correlation coefficient (PCC) and Gene Ontology (GO) enrichment revealed that targets of miRNAs with positive correlations are integral to chlorophyll-related photosynthetic processes. In contrast, negatively correlated pairs are primarily associated with metabolic functions. It is worth noting that no C4 genes were predicted as targets of DEmiRNAs. Our application of weighted gene co-expression network analysis (WGCNA) led to a gene regulatory network (GRN) suggesting miRNAs might indirectly regulate C4 genes via transcription factors (TFs). The GRAS TF SsHAM3a emerged as a potential regulator of C4 genes, targeted by miR171y and miR171am, and exhibiting a negative correlation with miRNA expression along the leaf gradient. This study sheds light on the complex involvement of miRNAs in regulating C4 genes, offering a foundation for future research into enhancing sugarcane's photosynthetic efficiency.

The bHLH transcription factor OsPRI1 activates the Setaria viridis PEPC1 promoter in rice

Photosynthetic efficiency is reduced by the dual role of Rubisco, which acts either as a carboxylase or as an oxygenase, the latter leading to photorespiration. C4 photosynthesis evolved as a carbon-concentrating mechanism to reduce photorespiration. To engineer C4 into a C3 plant, it is essential to understand how C4 genes, such as phosphoenolpyruvate carboxylase (PEPC1), are regulated to be expressed at high levels and in a cell-specific manner. Yeast one-hybrid screening was used to show that OsPRI1, a rice bHLH transcription factor involved in iron homeostasis, binds to the Setaria viridis PEPC1 promoter. This promoter drives mesophyll-specific gene expression in rice. The role of OsPRI1 in planta was characterized using a rice line harbouring SvPEPC1pro ::GUS. We show that OsPRI1 activates the S. viridis PEPC1 promoter by binding to an N-box in the proximal promoter, and that GUS activity is highly reduced in SvPEPC1pro ::GUS lines when OsPRI1 is mutated. Cross-species comparisons showed that the SvPRI1 homolog binds to the SvPEPC1 promoter but the maize ZmPRI1 does not bind to the ZmPEPC1 promoter. Our results suggest that elements of the iron homeostasis pathway were co-opted to regulate PEPC1 gene expression during the evolution of some but not all C4 species.

A Specific Protective Mechanism Against Chloroplast Photo-Reactive Oxygen Species in Phosphate-Starved Rice Plants

Phosphorus (Pi) starvation prevents a good match between light energy absorption and photosynthetic carbon metabolism, generating photo-reactive oxygen species (photo-ROS) in chloroplasts. Plants have evolved to withstand photo-oxidative stress, but the key regulatory mechanism underlying it remains unclear. In rice (Oryza sativa), DEEP GREEN PANICLE1 (DGP1) is robustly up-regulated in response to Pi deficiency. DGP1 decreases the DNA-binding capacities of the transcriptional activators GLK1/2 on the photosynthetic genes involved in chlorophyll biosynthesis, light harvesting, and electron transport. This Pi-starvation-induced mechanism dampens both electron transport rates through photosystem I and II (ETRI and ETRII) and thus mitigates the electron-excessive stress in mesophyll cells. Meanwhile, DGP1 hijacks glycolytic enzymes GAPC1/2/3, redirecting glucose metabolism toward the pentose phosphate pathway with superfluous NADPH production. Phenotypically, light irradiation induces O2 - production in Pi-starved WT leaves but is observably accelerated in dgp1 mutant and impaired in GAPCsRNAi and glk1glk2 lines. Interestingly, overexpressed DGP1 in rice caused hyposensitivity to ROS-inducers (catechin and methyl viologen), but the dgp1 mutant shows a similar inhibitory phenotype with the WT seedlings. Overall, the DGP1 gene serves as a specific antagonizer against photo-ROS in Pi-starved rice plants, which coordinates light-absorbing and anti-oxidative systems by orchestrating transcriptional and metabolic regulations, respectively.

Compartmentation of photosynthesis gene expression in C4 maize depends on time of day

Compared with the ancestral C3 state, C4 photosynthesis occurs at higher rates with improved water and nitrogen use efficiencies. In both C3 and C4 plants, rates of photosynthesis increase with light intensity and are maximal around midday. We determined that in the absence of light or temperature fluctuations, photosynthesis in maize (Zea mays) peaks in the middle of the subjective photoperiod. To investigate the molecular processes associated with these temporal changes, we performed RNA sequencing of maize mesophyll and bundle sheath strands over a 24-h time course. Preferential expression of C4 cycle genes in these cell types was strongest between 6 and 10 h after dawn when rates of photosynthesis were highest. For the bundle sheath, DNA motif enrichment and gene coexpression analyses suggested members of the DNA binding with one finger (DOF) and MADS (MINICHROMOSOME MAINTENANCE FACTOR 1/AGAMOUS/DEFICIENS/Serum Response Factor)-domain transcription factor families mediate diurnal fluctuations in C4 gene expression, while trans-activation assays in planta confirmed their ability to activate promoter fragments from bundle sheath expressed genes. The work thus identifies transcriptional regulators and peaks in cell-specific C4 gene expression coincident with maximum rates of photosynthesis in the maize leaf at midday.

Leaf photosynthetic characteristics of waxy maize in response to different degrees of heat stress during grain filling

BACKGROUND

In the context of climate change, maize is facing unprecedented heat stress (HS) threats during grain filling. Understanding how HS affects yield is the key to reducing the impact of climate change on maize production. Suyunuo5 (SYN5) and Yunuo7 (YN7) were used as materials, and four temperature gradients of 28℃ (day)/20℃ (night; T0, control), 32 °C/24°C (T1, mild HS), 36 °C/28°C (T2, moderate HS), and 40 °C/32°C (T3, severe HS) were set up during grain filling to explore the physiological mechanism of different degrees HS affecting photosynthetic characteristics of leaves in this study.

RESULTS

Results showed that HS accelerated the degradation of chlorophyll, disturbed the metabolism of reactive oxygen species, reduced the activity of antioxidant enzymes, and caused leaf damage. Heat stress induced the down-regulation of photosynthesis-related genes, which results in the decrease of enzymatic activities involved in photosynthesis, thereby inhibiting photosynthesis and reducing yield. Integrated analysis showed that the degree of the negative influence of three HS types during grain filling on leaves and yield was T3 > T2 > T1. The increase in HS disturbed leaf physiological activities and grain filling. Meanwhile, this study observed that the YN7 was more heat tolerance than SYN5 and thus it was recommended to use YN7 in waxy maize planting areas with frequent high temperatures.

CONCLUSIONS

Heat stress during grain filling caused premature senescence of the leaves by inhibiting the ability of leaves to photosynthesize and accelerating the oxidative damage of cells, thereby affecting the waxy maize yield. Our study helped to simulate the productivity of waxy maize under high temperatures and provided assistance for a stable yield of waxy maize under future climate warming.

Impact of pod and seed photosynthesis on seed filling and canopy carbon gain in soybean

There is a limited understanding of the carbon assimilation capacity of nonfoliar green tissues and its impact on yield and seed quality since most photosynthesis research focuses on leaf photosynthesis. In this study, we investigate the photosynthetic efficiency of soybean (Glycine max) pods and seeds in a field setting and evaluate its effect on mature seed weight and composition. We demonstrate that soybean pod and seed photosynthesis contributes 13% to 14% of the mature seed weight. Carbon assimilation by soybean pod and seed photosynthesis can compensate for 81% of carbon loss through the respiration of the same tissues, and our model predicts that soybean pod and seed photosynthesis contributes up to 9% of the total daily carbon gain of the canopy. Chlorophyll fluorescence (CF) shows that the operating efficiency of photosystem II in immature soybean seeds peaks at the 10 to 100 mg seed weight stage, while that of immature pods peaks at the 75 to 100 mg stage. This study provides quantitative information about the efficiency of soybean pod and seed photosynthesis during tissue development and its impact on yield.

Regulation of photosynthesis by mitogen-activated protein kinase in rice: antagonistic adjustment by OsMPK3 and OsMPK6

Photosynthesis is the basis of almost all life on earth and is the main component of crop yield that contributes to the carbohydrate partitioning to the grains. Maintaining the photosynthetic efficiency of plants in challenging environmental conditions by regulating the associated factors is a potential research arena which will help in the improvement of crop yield. Phosphorylation is known to play a pivotal role in the regulation of photosynthesis. Mitogen Activated Protein Kinases (MAPKs) cascade although known to regulate a diverse range of processes does not have any exact reported function in the regulation of photosynthesis. To elucidate the regulatory role of MAPKs in photosynthesis we investigated the changes in net photosynthesis rate and related parameters in DEX inducible over-expressing (OE) lines of two members of MAPK gene family namely, OsMPK3 and OsMPK6 in rice. Interestingly, significant changes were found in net photosynthesis rate and related physiological parameters in OsMPK3 and OsMPK6-OE lines compared to its wild-type relatives. OsMPK3 and OsMPK6 have regulatory effects on nuclear-encoded photosynthetic genes. Untargeted metabolite profiling reveals a higher accumulation of sugars and their derivatives in MPK6 overexpressing plants and a lower accumulation of sugars and organic acids in MPK3 overexpressing plants. The accumulation of amino acids was found in abundance in both MPK3 and MPK6 overexpressing plants. Understanding the effects of MPK3 and MPK6 on the CO2 assimilation of rice plants under normal growth conditions, will help in devising strategies that can be extended for crop improvement.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01383-9.

A rapid method to quantify vein density in C4 plants using starch staining

C4 photosynthesis has evolved multiple times in the angiosperms and typically involves alterations to the biochemistry, cell biology and development of leaves. One common modification found in C4 plants compared with the ancestral C3 state is an increase in vein density such that the leaf contains a larger proportion of bundle sheath cells. Recent findings indicate that there may be significant intraspecific variation in traits such as vein density in C4 plants but to use such natural variation for trait-mapping, rapid phenotyping would be required. Here we report a high-throughput method to quantify vein density that leverages the bundle sheath-specific accumulation of starch found in C4 species. Starch staining allowed high-contrast images to be acquired permitting image analysis with MATLAB- and Python-based programmes. The method works for dicotyledons and monocotolydons. We applied this method to Gynandropsis gynandra where significant variation in vein density was detected between natural accessions, and Zea mays where no variation was apparent in the genotypically diverse lines assessed. We anticipate this approach will be useful to map genes controlling vein density in C4 species demonstrating natural variation for this trait.

Differential subgenome expression underlies biomass accumulation in allotetraploid Pennisetum giganteum

BACKGROUND

Pennisetum giganteum (AABB, 2n = 4x = 28) is a C4 plant in the genus Pennisetum with origin in Africa but currently also grown in Asia and America. It is a crucial forage and potential energy grass with significant advantages in yield, stress resistance, and environmental adaptation. However, the mechanisms underlying these advantageous traits remain largely unexplored. Here, we present a high-quality genome assembly of the allotetraploid P. giganteum aiming at providing insights into biomass accumulation.

RESULTS

Our assembly has a genome size 2.03 Gb and contig N50 of 88.47 Mb that was further divided into A and B subgenomes. Genome evolution analysis revealed the evolutionary relationships across the Panicoideae subfamily lineages and identified numerous genome rearrangements that had occurred in P. giganteum. Comparative genomic analysis showed functional differentiation between the subgenomes. Transcriptome analysis found no subgenome dominance at the overall gene expression level; however, differentially expressed homoeologous genes and homoeolog-specific expressed genes between the two subgenomes were identified, suggesting that complementary effects between the A and B subgenomes contributed to biomass accumulation of P. giganteum. Besides, C4 photosynthesis-related genes were significantly expanded in P. giganteum and their sequences and expression patterns were highly conserved between the two subgenomes, implying that both subgenomes contributed greatly and almost equally to the highly efficient C4 photosynthesis in P. giganteum. We also identified key candidate genes in the C4 photosynthesis pathway that showed sustained high expression across all developmental stages of P. giganteum.

CONCLUSIONS

Our study provides important genomic resources for elucidating the genetic basis of advantageous traits in polyploid species, and facilitates further functional genomics research and genetic improvement of P. giganteum.

Spatial expression patterns of genes encoding sugar sensors in leaves of C4 and C3 grasses

BACKGROUND AND AIMS

The mechanisms of sugar sensing in grasses remain elusive, especially those using C4 photosynthesis even though a large proportion of the world's agricultural crops utilize this pathway. We addressed this gap by comparing the expression of genes encoding components of sugar sensors in C3 and C4 grasses, with a focus on source tissues of C4 grasses. Given C4 plants evolved into a two-cell carbon fixation system, it was hypothesized this may have also changed how sugars were sensed.

METHODS

For six C3 and eight C4 grasses, putative sugar sensor genes were identified for target of rapamycin (TOR), SNF1-related kinase 1 (SnRK1), hexokinase (HXK) and those involved in the metabolism of the sugar sensing metabolite trehalose-6-phosphate (T6P) using publicly available RNA deep sequencing data. For several of these grasses, expression was compared in three ways: source (leaf) versus sink (seed), along the gradient of the leaf, and bundle sheath versus mesophyll cells.

KEY RESULTS

No positive selection of codons associated with the evolution of C4 photosynthesis was identified in sugar sensor proteins here. Expressions of genes encoding sugar sensors were relatively ubiquitous between source and sink tissues as well as along the leaf gradient of both C4 and C3 grasses. Across C4 grasses, SnRK1β1 and TPS1 were preferentially expressed in the mesophyll and bundle sheath cells, respectively. Species-specific differences of gene expression between the two cell types were also apparent.

CONCLUSIONS

This comprehensive transcriptomic study provides an initial foundation for elucidating sugar-sensing genes within major C4 and C3 crops. This study provides some evidence that C4 and C3 grasses do not differ in how sugars are sensed. While sugar sensor gene expression has a degree of stability along the leaf, there are some contrasts between the mesophyll and bundle sheath cells.

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.

Plastids: diving into their diversity, their functions, and their role in plant development

Plastids are a group of essential, heterogenous semi-autonomous organelles characteristic of plants that perform photosynthesis and a diversity of metabolic pathways that impact growth and development. Plastids are remarkably dynamic and can interconvert in response to specific developmental and environmental cues, functioning as a central metabolic hub in plant cells. By far the best studied plastid is the chloroplast, but in recent years the combination of modern techniques and genetic analyses has expanded our current understanding of plastid morphological and functional diversity in both model and non-model plants. These studies have provided evidence of an unexpected diversity of plastid subtypes with specific characteristics. In this review, we describe recent findings that provide insights into the characteristics of these specialized plastids and their functions. We concentrate on the emerging evidence that supports the model that signals derived from particular plastid types play pivotal roles in plant development, environmental, and defense responses. Furthermore, we provide examples of how new technologies are illuminating the functions of these specialized plastids and the overall complexity of their differentiation processes. Finally, we discuss future research directions such as the use of ectopic plastid differentiation as a valuable tool to characterize factors involved in plastid differentiation. Collectively, we highlight important advances in the field that can also impact future agricultural and biotechnological improvement in plants.

C4 gene induction during de-etiolation evolved through changes in cis to allow integration with ancestral C3 gene regulatory networks

C₄ photosynthesis has evolved by repurposing enzymes found in C₃ plants. Compared with the ancestral C₃ state, accumulation of C₄ cycle proteins is enhanced. We used de-etiolation of C₄ Gynandropsis gynandra and C₃ Arabidopsis thaliana to understand this process. C₄ gene expression and chloroplast biogenesis in G. gynandra were tightly coordinated. Although C₃ and C₄ photosynthesis genes showed similar induction patterns, in G. gynandra, C₄ genes were more strongly induced than orthologs from A. thaliana. In vivo binding of TGA and homeodomain as well as light-responsive elements such as G- and I-box motifs were associated with the rapid increase in transcripts of C₄ genes. Deletion analysis confirmed that regions containing G- and I-boxes were necessary for high expression. The data support a model in which accumulation of transcripts derived from C₄ photosynthesis genes in C₄ leaves is enhanced because modifications in cis allowed integration into ancestral transcriptional networks.

Genome-wide identification and comparative analyses of key genes involved in C4 photosynthesis in five main gramineous crops

Compared to C₃ species, C₄ plants showed higher photosynthetic capacity as well as water and nitrogen use efficiency due to the presence of the C₄ photosynthetic pathway. Previous studies have shown that all genes required for the C₄ photosynthetic pathway exist in the genomes of C₃ species and are expressed. In this study, the genes encoding six key C₄ photosynthetic pathway enzymes (β-CA, PEPC, ME, MDH, RbcS, and PPDK) in the genomes of five important gramineous crops (C₄: maize, foxtail millet, and sorghum; C₃: rice and wheat) were systematically identified and compared. Based on sequence characteristics and evolutionary relationships, their C₄ functional gene copies were distinguished from non-photosynthetic functional gene copies. Furthermore, multiple sequence alignment revealed important sites affecting the activities of PEPC and RbcS between the C₃ and C₄ species. Comparisons of expression characteristics confirmed that the expression patterns of non-photosynthetic gene copies were relatively conserved among species, while C₄ gene copies in C₄ species acquired new tissue expression patterns during evolution. Additionally, multiple sequence features that may affect C₄ gene expression and subcellular localization were found in the coding and promoter regions. Our work emphasized the diversity of the evolution of different genes in the C₄ photosynthetic pathway and confirmed that the specific high expression in the leaf and appropriate intracellular distribution were the keys to the evolution of C₄ photosynthesis. The results of this study will help determine the evolutionary mechanism of the C₄ photosynthetic pathway in Gramineae and provide references for the transformation of C₄ photosynthetic pathways in wheat, rice, and other major C₃ cereal crops.

Evolution of gene regulatory network of C4 photosynthesis in the genus Flaveria reveals the evolutionary status of C3-C4 intermediate species

C₄ photosynthesis evolved from ancestral C₃ photosynthesis by recruiting pre-existing genes to fulfill new functions. The enzymes and transporters required for the C₄ metabolic pathway have been intensively studied and well documented; however, the transcription factors (TFs) that regulate these C₄ metabolic genes are not yet well understood. In particular, how the TF regulatory network of C₄ metabolic genes was rewired during the evolutionary process is unclear. Here, we constructed gene regulatory networks (GRNs) for four closely evolutionarily related species from the genus Flaveria, which represent four different evolutionary stages of C₄ photosynthesis: C₃ (F. robusta), type I C₃-C₄ (F. sonorensis), type II C₃-C₄ (F. ramosissima), and C₄ (F. trinervia). Our results show that more than half of the co-regulatory relationships between TFs and core C₄ metabolic genes are species specific. The counterparts of the C₄ genes in C₃ species were already co-regulated with photosynthesis-related genes, whereas the required TFs for C₄ photosynthesis were recruited later. The TFs involved in C₄ photosynthesis were widely recruited in the type I C₃-C₄ species; nevertheless, type II C₃-C₄ species showed a divergent GRN from C₄ species. In line with these findings, a ¹³CO₂ pulse-labeling experiment showed that the CO₂ initially fixed into C₄ acid was not directly released to the Calvin-Benson-Bassham cycle in the type II C₃-C₄ species. Therefore, our study uncovered dynamic changes in C₄ genes and TF co-regulation during the evolutionary process; furthermore, we showed that the metabolic pathway of the type II C₃-C₄ species F. ramosissima represents an alternative evolutionary solution to the ammonia imbalance in C₃-C₄ intermediate species.

Encoded C4 homologue enzymes genes function under abiotic stresses in C3 plant

Plant organisms assimilate CO₂ through the photosynthetic pathway, which facilitates in the synthesis of sugar for plant development. As environmental elements including water level, CO₂ concentration, temperature and soil characteristics change, the plants may recruit series of genes to help adapt the hostile environments and challenges. C4 photosynthesis plants are an excellent example of plant evolutionary adaptation to diverse condition. Compared with C3 photosynthesis plants, C4 photosynthesis plants have altered leaf anatomy and new metabolism for CO₂ capture, with multiple related enzymes such as phosphoenolpyruvate carboxylase (PEPCase), pyruvate orthophosphate dikinase (PPDK), NAD(P)-malic enzyme (NAD(P)-ME), NAD(P) - malate dehydrogenase (NAD(P)-MDH) and carbonic anhydrases (CA), identified to participate in the carbon concentrating mechanism (CCM) pathway. Recently, great achievements about C4 CCM-related genes have been made in the dissection of C3 plant development processes involving various stresses. In this review, we describe the functions of C4 CCM-related homologous genes in carbon and nitrogen metabolism in C3 plants. We further summarize C4 CCM-related homologous genes' functions in response to stresses in C3 plants. The understanding of C4 CCM-related genes' function in response to abiotic stress in plant is important to modify the crop plants for climate diversification.

Development of C4 Biochemistry and Change in Expression of Markers for Photosystems I and II in the Single-Cell C4 Species, Bienertia sinuspersici

Bienertia sinuspersici is one of four identified terrestrial plants that perform C₄ photosynthesis within a single chlorenchyma cell via the compartmentation of organelles and photosynthetic enzymes. The patterns of accumulation of key photosynthetic enzymes and transcripts in developing leaves were examined using immunolocalization and in situ hybridization. The polypeptides of Rubisco large subunit (RbcL) and pyruvate Pi dikinase (PPDK) accumulated equally in all chloroplasts before the formation of two intracellular cytoplasmic compartments: the central (CCC) and peripheral (PCC) cytoplasmic compartments. The differential accumulation of these enzymes was not completed until the leaf had reached maturity, indicating that the transition from C₃ to C₄ photosynthesis occurred during leaf maturation. In mature chlorenchyma cells, RbcL accumulated 20-fold higher in the CCC than in the PCC, while PPDK exhibited a concentration gradient that was the lowest in the chloroplasts in the central region of the CCC and the highest in PCC chloroplasts. The pattern of rbcL transcript accumulation followed that of its polypeptides in developing leaves, suggesting that the expression of this gene was likely controlled by transcriptional and/or post-transcriptional processes. Immunocytochemical results examining the distribution of photosystems I and II in the chloroplasts of chlorenchyma cells from mature leaves showed that PSII is more abundant in chloroplasts of the central compartment, whereas PSI is higher in those of the peripheral compartment. The quantitative real-time PCR results of rbcL, psbA, and psaB transcripts from the isolated chloroplasts of each compartment further supported this observation. Our results suggest that multiple levels of regulation play a role in controlling the differential accumulation of photosynthetic gene expression in the dimorphic chloroplasts of single-cell C₄ species during leaf development.

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.

Progress and prospects of C4 trait engineering in plants

Incorporating C₄ photosynthetic traits into C₃ crops is a rational approach for sustaining future demands for crop productivity. Using classical plant breeding, engineering this complex trait is unlikely to achieve its target. Therefore, it is critical and timely to implement novel biotechnological crop improvement strategies to accomplish this goal. However, a fundamental understanding of C₃ , C₄ , and C₃ -C₄ intermediate metabolism is crucial for the targeted use of biotechnological tools. This review assesses recent progress towards engineering C₄ photosynthetic traits in C₃ crops. We also discuss lessons learned from successes and failures of recent genetic engineering attempts in C₃ crops, highlighting the pros and cons of using rice as a model plant for short-, medium- and long-term goals of genetic engineering. This review provides an integrated approach towards engineering improved photosynthetic efficiency in C₃ crops for sustaining food, fibre and fuel production around the globe.

Photosynthetic mechanism of high yield under an improved wide-narrow row planting pattern in maize

Wide-narrow row maize planting patterns are a popular way to enhance maize yield via improving canopy PAR. To further optimize canopy PAR, we designed an improved wide-narrow row planting pattern (R2) based on the principle of the shortest projection length and the longest illumination of objects on the ground. Compared to the traditional wide-narrow row planting pattern (R1), maize yield increased by about 10% in R2. R2 maize had higher PAR, leaf area index, chlorophyll content, and photosynthetic rates than maize grown in R1. Moreover, compared to maize leaves in R1, the carbon assimilation enzymatic activities were also significantly higher in R2. The higher carbon assimilation enzymatic activity in R2 could account for the increased photosynthetic rate. Thus, the improved wide-narrow row planting pattern could improve photosynthetic performance by enhancing the PAR of the plant canopy, which further promotes the ear number and yield in northeast China.

Evolution of an intermediate C4 photosynthesis in the non-foliar tissues of the Poaceae

Carbon concentrating mechanisms (CCMs) in plants are abaptive features that have evolved to sustain plant growth in unfavorable environments, especially at low atmospheric carbon levels and high temperatures. Uptake of CO₂ and its storage in the aerenchyma tissues of Lycopsids and diurnal acidity fluctuation in aquatic plants during the Palaeozoic era (ca. 300 Ma.) would represent the earliest evolution of a CCM. The CCM parts of the dark reactions of photosynthesis have evolved many times, while the light reactions are conserved across plant lineages. A C₄ type CCM, leaf C₄ photosynthesis is evolved in the PACMAD clade of the Poaceae family. The evolution of C₄ photosynthesis from C₃ photosynthesis was an abaptation. Photosynthesis in reproductive tissues of sorghum and maize (PACMAD clade) has been shown to be of a weaker C₄ type (high CO₂ compensation point, low carbon isotope discrimination, and lack of Rubisco compartmentalization, when compared to the normal C₄ types) than that in the leaves (normal C₄ type). However, this does not fit well with the character polarity concept from an evolutionary perspective. In a recent model proposed for CCM evolution, the development of a rudimentary CCM prior to the evolution of a more efficient CCM (features contrasting to a weaker C₄ type, leading to greater biomass production rate) has been suggested. An intermediate crassulacean acid metabolism (CAM) type of CCM (rudimentary) was reported in the genera, Brassia, Coryanthes, Eriopsis, Peristeria, of the orchids (well-known group of plants that display the CAM pathway). Similarly, we propose here the evolution of a rudimentary CCM (C₄-like type pathway) in the non-foliar tissues of the Poaceae, prior to the evolution of the C₄ pathway as identified in the leaves of the C₄ species of the PACMAD clade.

Plants and water in a changing world: a physiological and ecological perspective

The reduction of greenhouse gases (GHGs) emission by replacing fossil energy stocks with carbon–neutral fuels is a major topic of the political and scientific debate on environmental sustainability. Such shift in energy sources is expected to curtail the accumulation rate of atmospheric CO₂, which is a strong infrared absorber and thus contributes to the global warming effect. Although such change would produce desirable outputs, the consequences of a drastic decrease in atmospheric CO₂ (the substrate of photosynthesis) should be carefully considered in the light of its potential impact on ecosystems stability and agricultural productivity. Indeed, plants regulate CO₂ uptake and water loss through the same anatomical structure: the leaf stomata. A reduced CO₂ availability is thus expected to enhance transpiration rate in plants decreasing their water use efficiency and imposing an increased water demand for both agricultural and wild ecosystems. We suggest that this largely underestimated issue should be duly considered when implementing policies that aim at the mitigation of global environmental changes and, at the same time, promote sustainable agricultural practices, include the preservation of biodiversity. Also, we underlie the important role(s) that modern biotechnology could play to tackle these global challenges by introducing new traits aimed at creating crop varieties with enhanced CO₂ capture and water- and light-use efficiency.

A matter of time: regulatory events behind the synchronization of C4 and crassulacean acid metabolism in Portulaca oleracea

Portulaca species can switch between C4 and crassulacean acid metabolism (CAM) depending on environmental conditions. However, the regulatory mechanisms behind this rare photosynthetic adaptation remain elusive. Using Portulaca oleracea as a model system, here we investigated the involvement of the circadian clock, plant hormones, and transcription factors in coordinating C4 and CAM gene expression. Free-running experiments in constant conditions suggested that C4 and CAM gene expression are intrinsically connected to the circadian clock. Detailed time-course, drought, and rewatering experiments revealed distinct time frames for CAM induction and reversion (days versus hours, respectively), which were accompanied by changes in abscisic acid (ABA) and cytokinin metabolism and signaling. Exogenous ABA and cytokinins were shown to promote and repress CAM expression in P. oleracea, respectively. Moreover, the drought-induced decline in C4 transcript levels was completely recovered upon cytokinin treatment. The ABA-regulated transcription factor genes HB7, NFYA7, NFYC9, TT8, and ARR12 were identified as likely candidate regulators of CAM induction following this approach, whereas NFYC4 and ARR9 were connected to C4 expression patterns. Therefore, we provide insights into the signaling events controlling C4-CAM transitions in response to water availability and over the day/night cycle, highlighting candidate genes for future functional studies in the context of facultative C4-CAM photosynthesis.

Single-Cell Transcriptome and Network Analyses Unveil Key Transcription Factors Regulating Mesophyll Cell Development in Maize

Background: Maize mesophyll (M) cells play important roles in various biological processes such as photosynthesis II and secondary metabolism. Functional differentiation occurs during M-cell development, but the underlying mechanisms for regulating M-cell development are largely unknown. Results: We conducted single-cell RNA sequencing (scRNA-seq) to profile transcripts in maize leaves. We then identified coregulated modules by analyzing the resulting pseudo-time-series data through gene regulatory network analyses. WRKY, ERF, NAC, MYB and Heat stress transcription factor (HSF) families were highly expressed in the early stage, whereas CONSTANS (CO)-like (COL) and ERF families were highly expressed in the late stage of M-cell development. Construction of regulatory networks revealed that these transcript factor (TF) families, especially HSF and COL, were the major players in the early and later stages of M-cell development, respectively. Integration of scRNA expression matrix with TF ChIP-seq and Hi-C further revealed regulatory interactions between these TFs and their targets. HSF1 and COL8 were primarily expressed in the leaf bases and tips, respectively, and their targets were validated with protoplast-based ChIP-qPCR, with the binding sites of HSF1 being experimentally confirmed. Conclusions: Our study provides evidence that several TF families, with the involvement of epigenetic regulation, play vital roles in the regulation of M-cell development in maize.

Chromatin and regulatory differentiation between bundle sheath and mesophyll cells in maize

C₄ plants partition photosynthesis enzymes between the bundle sheath (BS) and the mesophyll (M) cells for the better delivery of CO₂ to RuBisCO and to reduce photorespiration. To better understand how C₄ photosynthesis is regulated at the transcriptional level, we performed RNA-seq, ATAC-seq, ChIP-seq and Bisulfite-seq (BS-seq) on BS and M cells isolated from maize leaves. By integrating differentially expressed genes with chromatin features, we found that chromatin accessibility coordinates with epigenetic features, especially H3K27me3 modification and CHH methylation, to regulate cell type-preferentially enriched gene expression. Not only the chromatin-accessible regions (ACRs) proximal to the genes (pACRs) but also the distal ACRs (dACRs) are determinants of cell type-preferentially enriched expression. We further identified cell type-preferentially enriched motifs, e.g. AAAG for BS cells and TGACC/T for M cells, and determined their corresponding transcription factors: DOFs and WRKYs. The complex interaction between cis and trans factors in the preferential expression of C₄ genes was also observed. Interestingly, cell type-preferentially enriched gene expression can be fine-tuned by the coordination of multiple chromatin features. Such coordination may be critical in ensuring the cell type-specific function of key C₄ genes. Based on the observed cell type-preferentially enriched expression pattern and coordinated chromatin features, we predicted a set of functionally unknown genes, e.g. Zm00001d042050 and Zm00001d040659, to be potential key C₄ genes. Our findings provide deep insight into the architectures associated with C₄ gene expression and could serve as a valuable resource to further identify the regulatory mechanisms present in C₄ species.

Proto Kranz-like leaf traits and cellular ionic regulation are associated with salinity tolerance in a halophytic wild rice

Species of wild rice (Oryza spp.) possess a wide range of stress tolerance traits that can be potentially utilized in breeding climate-resilient cultivated rice cultivars (Oryza sativa) thereby aiding global food security. In this study, we conducted a greenhouse trial to evaluate the salinity tolerance of six wild rice species, one cultivated rice cultivar (IR64) and one landrace (Pokkali) using a range of electrophysiological, imaging, and whole-plant physiological techniques. Three wild species (O. latifolia, O. officinalis and O. coarctata) were found to possess superior salinity stress tolerance. The underlying mechanisms, however, were strikingly different. Na+ accumulation in leaves of O. latifolia, O. officinalis and O. coarctata were significantly higher than the tolerant landrace, Pokkali. Na+ accumulation in mesophyll cells was only observed in O. coarctata, suggesting that O. officinalis and O. latifolia avoid Na+ accumulation in mesophyll by allocating Na+ to other parts of the leaf. The finding also suggests that O. coarctata might be able to employ Na+ as osmolyte without affecting its growth. Further study of Na+ allocation in leaves will be helpful to understand the mechanisms of Na+ accumulation in these species. In addition, O. coarctata showed Proto Kranz-like leaf anatomy (enlarged bundle sheath cells and lower numbers of mesophyll cells), and higher expression of C4-related genes (e.g., NADPME, PPDK) and was a clear outlier with respect to salinity tolerance among the studied wild and cultivated Oryza species. The unique phylogenetic relationship of O. coarctata with C4 grasses suggests the potential of this species for breeding rice with high photosynthetic rate under salinity stress in the future.

Respiratory and C4-photosynthetic NAD-malic enzyme coexist in bundle sheath cell mitochondria and evolved via association of differentially adapted subunits

In plant mitochondria, nicotinamide adenine dinucleotide-malic enzyme (NAD-ME) has a housekeeping function in malate respiration. In different plant lineages, NAD-ME was independently co-opted in C4 photosynthesis. In the C4 Cleome species, Gynandropsis gynandra and Cleome angustifolia, all NAD-ME genes (NAD-MEα, NAD-MEβ1, and NAD-MEβ2) were affected by C4 evolution and are expressed at higher levels than their orthologs in the C3 species Tarenaya hassleriana. In T. hassleriana, the NAD-ME housekeeping function is performed by two heteromers, NAD-MEα/β1 and NAD-MEα/β2, with similar biochemical properties. In both C4 species, this role is restricted to NAD-MEα/β2. In the C4 species, NAD-MEα/β1 is exclusively present in the leaves, where it accounts for most of the enzymatic activity. Gynandropsis gynandra NAD-MEα/β1 (GgNAD-MEα/β1) exhibits high catalytic efficiency and is differentially activated by the C4 intermediate aspartate, confirming its role as the C4-decarboxylase. During C4 evolution, NAD-MEβ1 lost its catalytic activity; its contribution to the enzymatic activity results from a stabilizing effect on the associated α-subunit and the acquisition of regulatory properties. We conclude that in bundle sheath cell mitochondria of C4 species, the functions of NAD-ME as C4 photosynthetic decarboxylase and as a housekeeping enzyme coexist and are performed by isoforms that combine the same α-subunit with differentially adapted β-subunits.

Distinct C4 sub-types and C3 bundle sheath isolation in the Paniceae grasses

In C4 plants, the enzymatic machinery underpinning photosynthesis can vary, with, for example, three distinct C4 acid decarboxylases being used to release CO2 in the vicinity of RuBisCO. For decades, these decarboxylases have been used to classify C4 species into three biochemical sub-types. However, more recently, the notion that C4 species mix and match C4 acid decarboxylases has increased in popularity, and as a consequence, the validity of specific biochemical sub-types has been questioned. Using five species from the grass tribe Paniceae, we show that, although in some species transcripts and enzymes involved in multiple C4 acid decarboxylases accumulate, in others, transcript abundance and enzyme activity is almost entirely from one decarboxylase. In addition, the development of a bundle sheath isolation procedure for a close C3 species in the Paniceae enables the preliminary exploration of C4 sub-type evolution.

Engineering chloroplast development in rice through cell-specific control of endogenous genetic circuits

The engineering of C₄ photosynthetic activity into the C₃ plant rice has the potential to nearly double rice yields. To engineer a two-cell photosynthetic system in rice, the rice bundle sheath (BS) must be rewired to enhance photosynthetic capacity. Here, we show that BS chloroplast biogenesis is enhanced when the transcriptional activator, Oryza sativa Cytokinin GATA transcription factor 1 (OsCGA1), is driven by a vascular specific promoter. Ectopic expression of OsCGA1 resulted in increased BS chloroplast planar area and increased expression of photosynthesis-associated nuclear genes (PhANG), required for the biogenesis of photosynthetically active chloroplasts in BS cells of rice. A further refinement using a DNAse dead Cas9 (dCas9) activation module driven by the same cell-type specific promoter, directed enhanced chloroplast development of the BS cells when gRNA sequences were delivered by the dCas9 module to the promoter of the endogenous OsCGA1 gene. Single gRNA expression was sufficient to mediate the transactivation of both the endogenous gene and a transgenic GUS reporter fused with OsCGA1 promoter. Our results illustrate the potential for tissue-specific dCas9-activation and the co-regulation of genes needed for multistep engineering of C₄ rice.

Change in expression levels of NAD kinase-encoding genes in Flaveria species

Nicotinamide adenine dinucleotides (NAD(H)) and NAD phosphates (NADP(H)) are electron carriers involved in redox reactions and metabolic processes in all organisms. NAD kinase (NADK) is the only enzyme that phosphorylates NAD⁺ into NADP⁺, using ATP as a phosphate donor. In NADP-dependent malic enzyme (NADP-ME)-type C₄ photosynthesis, NADP(H) are required for dehydrogenation by NADP-dependent malate dehydrogenase (NADP-MDH) in mesophyll cells, and decarboxylation by NADP-ME in bundle sheath cells. In this study, we identified five NADK genes (FbNADK1a, 1b, 2a, 2b, and 3) from the C₄ model species Flaveria bidentis. RNA-Seq database analysis revealed higher transcript abundance in one of the chloroplast-type NADK2 genes of C₄F. bidentis (FbNADK2a). Comparative analysis of NADK activity in leaves of C₃, C₃-C₄, and C₄Flaveria showed that C₄Flaveria (F. bidentis and F. trinervia) had higher NADK activity than the other photosynthetic-types of Flaveria. Taken together, our results suggest that chloroplastic NAD kinase appeared to increase in importance as C₃ plants evolved into C₄ plants in the genus Flaveria.

Challenges and Approaches to Crop Improvement Through C3-to-C4 Engineering

With a rapidly growing world population and dwindling natural resources, we are now facing the enormous challenge of increasing crop yields while simultaneously improving the efficiency of resource utilization. Introduction of C4 photosynthesis into C3 crops is widely accepted as a key strategy to meet this challenge because C4 plants are more efficient than C3 plants in photosynthesis and resource usage, particularly in hot climates, where the potential for productivity is high. Lending support to the feasibility of this C3-to-C4 engineering, evidence indicates that C4 photosynthesis has evolved from C3 photosynthesis in multiple lineages. Nevertheless, C3-to-C4 engineering is not an easy task, as several features essential to C4 photosynthesis must be introduced into C3 plants. One such feature is the spatial separation of the two phases of photosynthesis (CO₂ fixation and carbohydrate synthesis) into the mesophyll and bundle sheath cells, respectively. Another feature is the Kranz anatomy, characterized by a close association between the mesophyll and bundle sheath (BS) cells (1:1 ratio). These anatomical features, along with a C4-specific carbon fixation enzyme (PEPC), form a CO₂-concentration mechanism that ensures a high photosynthetic efficiency. Much effort has been taken in the past to introduce the C4 mechanism into C3 plants, but none of these attempts has met with success, which is in my opinion due to a lack of system-level understanding and manipulation of the C3 and C4 pathways. As a prerequisite for the C3-to-C4 engineering, I propose that not only the mechanisms that control the Kranz anatomy and cell-type-specific expression in C3 and C4 plants must be elucidated, but also a good understanding of the gene regulatory network underlying C3 and C4 photosynthesis must be achieved. In this review, I first describe the past and current efforts to increase photosynthetic efficiency in C3 plants and their limitations; I then discuss a systems approach to tackling down this challenge, some practical issues, and recent technical innovations that would help us to solve these problems.

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.

ZmOrphan94 Transcription Factor Downregulates ZmPEPC1 Gene Expression in Maize Bundle Sheath Cells

Spatial separation of the photosynthetic reactions is a key feature of C4 metabolism. In most C4 plants, this separation requires compartmentation of photosynthetic enzymes between mesophyll (M) and bundle sheath (BS) cells. The upstream region of the gene encoding the maize PHOSPHOENOLPYRUVATE CARBOXYLASE 1 (ZmPEPC1) has been shown sufficient to drive M-specific ZmPEPC1 gene expression. Although this region has been well characterized, to date, only few trans-factors involved in the ZmPEPC1 gene regulation were identified. Here, using a yeast one-hybrid approach, we have identified three novel maize transcription factors ZmHB87, ZmCPP8, and ZmOrphan94 as binding to the ZmPEPC1 upstream region. Bimolecular fluorescence complementation assays in maize M protoplasts unveiled that ZmOrphan94 forms homodimers and interacts with ZmCPP8 and with two other ZmPEPC1 regulators previously reported, ZmbHLH80 and ZmbHLH90. Trans-activation assays in maize M protoplasts unveiled that ZmHB87 does not have a clear transcriptional activity, whereas ZmCPP8 and ZmOrphan94 act as activator and repressor, respectively. Moreover, we observed that ZmOrphan94 reduces the trans-activation activity of both activators ZmCPP8 and ZmbHLH90. Using the electromobility shift assay, we showed that ZmOrphan94 binds to several cis-elements present in the ZmPEPC1 upstream region and one of these cis-elements overlaps with the ZmbHLH90 binding site. Gene expression analysis revealed that ZmOrphan94 is preferentially expressed in the BS cells, suggesting that ZmOrphan94 is part of a transcriptional regulatory network downregulating ZmPEPC1 transcript level in the BS cells. Based on both this and our previous work, we propose a model underpinning the importance of a regulatory mechanism within BS cells that contributes to the M-specific ZmPEPC1 gene expression.

Influence of salt stress on C4 photosynthesis in Miscanthus sinensis Anderss

Miscanthus sinensis Anderss. is a good candidate for C₄ bioenergy crop development for marginal lands. As one of the characteristics of marginal lands, salinization is a major limitation to agricultural production. The present work aimed to investigate the possible factors involved in the tolerance of M. sinensis C₄ photosynthesis to salinity stress. Seedlings of two accessions (salt-tolerant 'JM0119' and salt-sensitive 'JM0099') were subjected to 0 mm NaCl (control) or 250 mm NaCl (salt stress treatment) for 2 weeks. The chlorophyll content, parameters of photosynthesis and chlorophyll a fluorescence, activity of C₄ enzymes and expression of C₄ genes were measured. The results showed that photosynthesis rate, transpiration rate, chlorophyll content, PSII operating efficiency, coefficient of photochemical quenching, activity of phosphoenolpyruvate carboxylase (PEPC) and pyruvate, orthophosphate dikinase (PPDK) and gene expression of PEPC and PPDK under salinity were higher after long-term salinity exposure in 'JM0119' than in 'JM0099', while activity of NADP-malate dehydrogenase (NADP-MDH) and NADP-malic enzyme (NADP-ME), together with expression of NADP-MDH and NADP-ME, were much higher in 'JM0099' than in 'JM0119'. In conclusion, the increased photosynthetic capacity under long-term salt stress in the salt-tolerant relative to the salt-sensitive M. sinensis accession was mainly associated with non-stomatal factors, such as reduced chlorophyll loss, higher PSII operating efficiency, enhanced activity of PEPC and PPDK and relatively lower activity of NADP-ME.

Continued Adaptation of C4 Photosynthesis After an Initial Burst of Changes in the Andropogoneae Grasses

C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} photosynthesis is a complex trait that sustains fast growth and high productivity in tropical and subtropical conditions and evolved repeatedly in flowering plants. One of the major C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} lineages is Andropogoneae, a group of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}1200 grass species that includes some of the world’s most important crops and species dominating tropical and some temperate grasslands. Previous efforts to understand C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} evolution in the group have compared a few model C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} plants to distantly related C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} species so that changes directly responsible for the transition to C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} could not be distinguished from those that preceded or followed it. In this study, we analyze the genomes of 66 grass species, capturing the earliest diversification within Andropogoneae as well as their C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} relatives. Phylogenomics combined with molecular dating and analyses of protein evolution show that many changes linked to the evolution of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} photosynthesis in Andropogoneae happened in the Early Miocene, between 21 and 18 Ma, after the split from its C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} sister lineage, and before the diversification of the group. This initial burst of changes was followed by an extended period of modifications to leaf anatomy and biochemistry during the diversification of Andropogoneae, so that a single C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} origin gave birth to a diversity of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} phenotypes during 18 million years of speciation events and migration across geographic and ecological spaces. Our comprehensive approach and broad sampling of the diversity in the group reveals that one key transition can lead to a plethora of phenotypes following sustained adaptation of the ancestral state. [Adaptive evolution; complex traits; herbarium genomics; Jansenelleae; leaf anatomy; Poaceae; phylogenomics.]

A bipartite transcription factor module controlling expression in the bundle sheath of Arabidopsis thaliana

C₄ photosynthesis evolved repeatedly from the ancestral C₃ state, improving photosynthetic efficiency by ~50%. In most C₄ lineages, photosynthesis is compartmented between mesophyll and bundle sheath cells, but how gene expression is restricted to these cell types is poorly understood. Using the C₃ model Arabidopsis thaliana, we identified cis-elements and transcription factors driving expression in bundle sheath strands. Upstream of the bundle sheath preferentially expressed MYB76 gene, we identified a region necessary and sufficient for expression containing two cis-elements associated with the MYC and MYB families of transcription factors. MYB76 expression is reduced in mutant alleles for these transcription factors. Moreover, downregulated genes shared by both mutants are preferentially expressed in the bundle sheath. Our findings are broadly relevant for understanding the spatial patterning of gene expression, provide specific insights into mechanisms associated with the evolution of C₄ photosynthesis and identify a short tuneable sequence for manipulating gene expression in the bundle sheath.

Biosystems Design to Accelerate C3-to-CAM Progression

Global demand for food and bioenergy production has increased rapidly, while the area of arable land has been declining for decades due to damage caused by erosion, pollution, sea level rise, urban development, soil salinization, and water scarcity driven by global climate change. In order to overcome this conflict, there is an urgent need to adapt conventional agriculture to water-limited and hotter conditions with plant crop systems that display higher water-use efficiency (WUE). Crassulacean acid metabolism (CAM) species have substantially higher WUE than species performing C3 or C4 photosynthesis. CAM plants are derived from C3 photosynthesis ancestors. However, it is extremely unlikely that the C3 or C4 crop plants would evolve rapidly into CAM photosynthesis without human intervention. Currently, there is growing interest in improving WUE through transferring CAM into C3 crops. However, engineering a major metabolic plant pathway, like CAM, is challenging and requires a comprehensive deep understanding of the enzymatic reactions and regulatory networks in both C3 and CAM photosynthesis, as well as overcoming physiometabolic limitations such as diurnal stomatal regulation. Recent advances in CAM evolutionary genomics research, genome editing, and synthetic biology have increased the likelihood of successful acceleration of C3-to-CAM progression. Here, we first summarize the systems biology-level understanding of the molecular processes in the CAM pathway. Then, we review the principles of CAM engineering in an evolutionary context. Lastly, we discuss the technical approaches to accelerate the C3-to-CAM transition in plants using synthetic biology toolboxes.

Methods of analysis of chloroplast genomes of C3, Kranz type C4 and Single Cell C4 photosynthetic members of Chenopodiaceae

BACKGROUND

Chloroplast genome information is critical to understanding forms of photosynthesis in the plant kingdom. During the evolutionary process, plants have developed different photosynthetic strategies that are accompanied by complementary biochemical and anatomical features. Members of family Chenopodiaceae have species with C3 photosynthesis, and variations of C4 photosynthesis in which photorespiration is reduced by concentrating CO2 around Rubisco through dual coordinated functioning of dimorphic chloroplasts. Among dicots, the family has the largest number of C4 species, and greatest structural and biochemical diversity in forms of C4 including the canonical dual-cell Kranz anatomy, and the recently identified single cell C4 with the presence of dimorphic chloroplasts separated by a vacuole. This is the first comparative analysis of chloroplast genomes in species representative of photosynthetic types in the family.

RESULTS

Methodology with high throughput sequencing complemented with Sanger sequencing of selected loci provided high quality and complete chloroplast genomes of seven species in the family and one species in the closely related Amaranthaceae family, representing C3, Kranz type C4 and single cell C4 (SSC4) photosynthesis six of the eight chloroplast genomes are new, while two are improved versions of previously published genomes. The depth of coverage obtained using high-throughput sequencing complemented with targeted resequencing of certain loci enabled superior resolution of the border junctions, directionality and repeat region sequences. Comparison of the chloroplast genomes with previously sequenced plastid genomes revealed similar genome organization, gene order and content with a few revisions. High-quality complete chloroplast genome sequences resulted in correcting the orientation the LSC region of the published Bienertia sinuspersici chloroplast genome, identification of stop codons in the rpl23 gene in B. sinuspersici and B. cycloptera, and identifying an instance of IR expansion in the Haloxylon ammodendron inverted repeat sequence. The rare observation of a mitochondria-to-chloroplast inter-organellar gene transfer event was identified in family Chenopodiaceae.

CONCLUSIONS

This study reports complete chloroplast genomes from seven Chenopodiaceae and one Amaranthaceae species. The depth of coverage obtained using high-throughput sequencing complemented with targeted resequencing of certain loci enabled superior resolution of the border junctions, directionality, and repeat region sequences. Therefore, the use of high throughput and Sanger sequencing, in a hybrid method, reaffirms to be rapid, efficient, and reliable for chloroplast genome sequencing.

Overexpression of maize transcription factor mEmBP-1 increases photosynthesis, biomass, and yield in rice

Identifying new options to improve photosynthetic capacity is a major approach to improve crop yield potential. Here we report that overexpression of the gene encoding the transcription factor mEmBP-1 led to simultaneously increased expression of many genes in photosynthesis, including genes encoding Chl a,b-binding proteins (Lhca and Lhcb), PSII (PsbR3 and PsbW) and PSI reaction center subunits (PsaK and PsaN), chloroplast ATP synthase subunit, electron transport reaction components (Fd1 and PC), and also major genes in the Calvin-Benson-Bassham cycle, including those encoding Rubisco, glyceraldehyde phosphate dehydrogenase, fructose bisphosphate aldolase, transketolase, and phosphoribulokinase. These increased expression of photosynthesis genes resulted in increased leaf chlorophyll pigment, photosynthetic rate, biomass growth, and grain yield both in the greenhouse and in the field. Using EMSA experiments, we showed that mEmBP-1a protein can directly bind to the promoter region of photosynthesis genes, suggesting that the direct binding of mEmBP-1a to the G-box domain of photosynthetic genes up-regulates expression of these genes. Altogether, our results show that mEmBP-1a is a major regulator of photosynthesis, which can be used to increase rice photosynthesis and yield in the field.

Genetic Diversity of C4 Photosynthesis Pathway Genes in Sorghum bicolor (L.)

C4 photosynthesis has evolved in over 60 different plant taxa and is an excellent example of convergent evolution. Plants using the C4 photosynthetic pathway have an efficiency advantage, particularly in hot and dry environments. They account for 23% of global primary production and include some of our most productive cereals. While previous genetic studies comparing phylogenetically related C3 and C4 species have elucidated the genetic diversity underpinning the C4 photosynthetic pathway, no previous studies have described the genetic diversity of the genes involved in this pathway within a C4 crop species. Enhanced understanding of the allelic diversity and selection signatures of genes in this pathway may present opportunities to improve photosynthetic efficiency, and ultimately yield, by exploiting natural variation. Here, we present the first genetic diversity survey of 8 known C4 gene families in an important C4 crop, Sorghum bicolor (L.) Moench, using sequence data of 48 genotypes covering wild and domesticated sorghum accessions. Average nucleotide diversity of C4 gene families varied more than 20-fold from the NADP-malate dehydrogenase (MDH) gene family (θπ = 0.2 × 10-3) to the pyruvate orthophosphate dikinase (PPDK) gene family (θπ = 5.21 × 10-3). Genetic diversity of C4 genes was reduced by 22.43% in cultivated sorghum compared to wild and weedy sorghum, indicating that the group of wild and weedy sorghum may constitute an untapped reservoir for alleles related to the C4 photosynthetic pathway. A SNP-level analysis identified purifying selection signals on C4 PPDK and carbonic anhydrase (CA) genes, and balancing selection signals on C4 PPDK-regulatory protein (RP) and phosphoenolpyruvate carboxylase (PEPC) genes. Allelic distribution of these C4 genes was consistent with selection signals detected. A better understanding of the genetic diversity of C4 pathway in sorghum paves the way for mining the natural allelic variation for the improvement of photosynthesis.

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.

Interspecific analysis of diurnal gene regulation in panicoid grasses identifies known and novel regulatory motifs

Background

The circadian clock drives endogenous 24-h rhythms that allow organisms to adapt and prepare for predictable and repeated changes in their environment throughout the day-night (diurnal) cycle. Many components of the circadian clock in Arabidopsis thaliana have been functionally characterized, but comparatively little is known about circadian clocks in grass species including major crops like maize and sorghum.

Results

Comparative research based on protein homology and diurnal gene expression patterns suggests the function of some predicted clock components in grasses is conserved with their Arabidopsis counterparts, while others have diverged in function. Our analysis of diurnal gene expression in three panicoid grasses sorghum, maize, and foxtail millet revealed conserved and divergent evolution of expression for core circadian clock genes and for the overall transcriptome. We find that several classes of core circadian clock genes in these grasses differ in copy number compared to Arabidopsis, but mostly exhibit conservation of both protein sequence and diurnal expression pattern with the notable exception of maize paralogous genes. We predict conserved cis-regulatory motifs shared between maize, sorghum, and foxtail millet through identification of diurnal co-expression clusters for a subset of 27,196 orthologous syntenic genes. In this analysis, a Cochran–Mantel–Haenszel based method to control for background variation identified significant enrichment for both expected and novel 6–8 nucleotide motifs in the promoter regions of genes with shared diurnal regulation predicted to function in common physiological activities.

Conclusions

This study illustrates the divergence and conservation of circadian clocks and diurnal regulatory networks across syntenic orthologous genes in panacoid grass species. Further, conserved local regulatory sequences contribute to the architecture of these diurnal regulatory networks that produce conserved patterns of diurnal gene expression.