Characterizing regulatory and functional differentiation between maize mesophyll and bundle sheath cells by transcriptomic analysis

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.

Fluctuating light induces a significant photoinhibition of photosystem I in maize

In nature, light intensity usually fluctuates and a sudden shade-sun transition can induce photodamage to photosystem I (PSI) in many angiosperms. Photosynthetic regulation in fluctuating light (FL) has been studied extensively in C3 plants; however, little is known about how C4 plants cope FL to prevent PSI photoinhibition. We here compared photosynthetic responses to FL between maize (Zea mays, C4) and tomato (Solanum lycopersicum, C3) grown under full sunlight. Maize leaves had significantly higher cyclic electron flow (CEF) activity and lower photorespiration activity than tomato. Upon a sudden shade-sun transition, maize showed a significant stronger transient PSI over-reduction than tomato, resulting in a significant greater PSI photoinhibition in maize after FL treatment. During the first seconds upon shade-sun transition, CEF was stimulated in maize at a much higher extent than tomato, favoring the rapid formation of trans-thylakoid proton gradient (ΔpH), which was helped by a transient down-regulation of chloroplast ATP synthase activity. Therefore, modulation of ΔpH by regulation of CEF and chloroplast ATP synthase adjusted PSI redox state at donor side, which partially compensated for the deficiency of photorespiration. We propose that C4 plants use different photosynthetic strategies for coping with FL as compared with C3 plants.

Maize GOLDEN2-LIKE proteins enhance drought tolerance in rice by promoting stomatal closure

Drought has become one of the most severe abiotic stresses experienced in agricultural production across the world. Plants respond to water deficit via stomatal movements in the leaves, which are mainly regulated by abscisic acid (ABA). A previous study from our lab showed that constitutive expression of maize (Zea mays L.) GOLDEN2-LIKE (GLK) transcription factors in rice (Oryza sativa L.) can improve stomatal conductance and plant photosynthetic capacity under field conditions. In the present study, we uncovered a function of ZmGLK regulation of stomatal movement in rice during drought stress. We found that elevated drought tolerance in rice plants overexpressing ZmGLK1 or GOLDEN2 (ZmG2) was conferred by rapid ABA-mediated stomatal closure. Comparative analysis of RNA-sequencing (RNA-seq) data from the rice leaves and DNA affinity purification sequencing (DAP-seq) results obtained in vitro revealed that ZmGLKs played roles in regulating ABA-related and stress-responsive pathways. Four upregulated genes closely functioning in abiotic stress tolerance with strong binding peaks in the DAP-seq data were identified as putative target genes of ZmGLK1 and ZmG2 in rice. These results demonstrated that maize GLKs play an important role in regulating stomatal movements to coordinate photosynthesis and stress tolerance. This trait is a valuable target for breeding drought-tolerant crop plants without compromising photosynthetic capacity.

Regulatory NADH dehydrogenase-like complex optimizes C4 photosynthetic carbon flow and cellular redox in maize

C4 plants typically operate a CO2 concentration mechanism from mesophyll (M) cells into bundle sheath (BS) cells. NADH dehydrogenase-like (NDH) complex is enriched in the BS cells of many NADP-malic enzyme (ME) type C4 plants and is more abundant in C4 than in C3 plants, but to what extent it is involved in the CO2 concentration mechanism remains to be experimentally investigated. We created maize and rice mutants deficient in NDH function and then used a combination of transcriptomic, proteomic, and metabolomic approaches for comparative analysis. Considerable decreases in growth, photosynthetic activities, and levels of key photosynthetic proteins were observed in maize but not rice mutants. However, transcript abundance for many cyclic electron transport (CET) and Calvin-Benson cycle components, as well as BS-specific C4 enzymes, was increased in maize mutants. Metabolite analysis of the maize ndh mutants revealed an increased NADPH : NADP ratio, as well as malate, ribulose 1,5-bisphosphate (RuBP), fructose 1,6-bisphosphate (FBP), and photorespiration intermediates. We suggest that by optimizing NADPH and malate levels and adjusting NADP-ME activity, NDH functions to balance metabolic and redox states in the BS cells of maize (in addition to ATP supply), coordinating photosynthetic transcript abundance and protein content, thus directly regulating the carbon flow in the two-celled C4 system of maize.

GOLDEN2-like1 is sufficient but not necessary for chloroplast biogenesis in mesophyll cells of C4 grasses

Chloroplasts are the site of photosynthesis. In land plants, chloroplast biogenesis is regulated by a family of transcription factors named GOLDEN2-like (GLK). In C4 grasses, it has been hypothesized that genome duplication events led to the sub-functionalization of GLK paralogs (GLK1 and GLK2) to control chloroplast biogenesis in two distinct cell types: mesophyll and bundle sheath cells. Although previous characterization of golden2 (g2) mutants in maize has demonstrated a role for GLK2 paralogs in regulating chloroplast biogenesis in bundle sheath cells, the function of GLK1 has remained elusive. Here we show that, contrary to expectations, GLK1 is not required for chloroplast biogenesis in mesophyll cells of maize. Comparisons between maize and Setaria viridis, which represent two independent C4 origins within the Poales, further show that the role of GLK paralogs in controlling chloroplast biogenesis in mesophyll and bundle sheath cells differs between species. Despite these differences, complementation analysis revealed that GLK1 and GLK2 genes from maize are both sufficient to restore functional chloroplast development in mesophyll and bundle sheath cells of S. viridis mutants. Collectively our results suggest an evolutionary trajectory in C4 grasses whereby both orthologs retained the ability to induce chloroplast biogenesis but GLK2 adopted a more prominent developmental role, particularly in relation to chloroplast activation in bundle sheath cells.

C4 leaf development and evolution

C4 photosynthesis is more efficient than C3 photosynthesis for two reasons. First, C4 plants have evolved efficient C4 enzymes to suppress wasteful photorespiration and enhance CO2 fixation. Second, C4 leaves have Kranz anatomy in which the veins are surrounded by one layer of bundle sheath (BS) cells and one layer of mesophyll (M) cells. The BS and M cells are functionally well differentiated and also well coordinated for rapid assimilation of atmospheric CO2 and transport of photo-assimilates between the two types of cells. Recent comparative transcriptomics of developing M and BS cells in young maize embryonic leaves revealed not only potential regulators of BS and M cell differentiation but also rapid early BS cell differentiation whereas slower, more prolonged M cell differentiation, contrary to the traditional view of a far simpler process of M cell development. Moreover, new upstream regulators of Kranz anatomy development have been identified and a number of gene co-expression modules for early vascular development have been inferred. Also, a candidate gene regulatory network associated with Kranz anatomy and vascular development has been constructed. Additionally, how whole genome duplication (WGD) may facilitate C4 evolution has been studied and the reasons for why the same WGD event led to successful C4 evolution in Gynandropsis gynandra but not in the sister species Tarenaya hassleriana have been proposed. Finally, new future research directions are suggested.

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.

Photosynthesis and food security: the evolving story of C4 rice

Traditional "Green Revolution" cereal breeding strategies to improve yield are now reaching a plateau in our principal global food crop rice. Photosynthesis has now become a major target of international consortia to increase yield potential. Synthetic biology is being used across multiple large projects to improve photosynthetic efficiency. This review follows the genesis and progress of one of the first of these consortia projects, now in its 13th year; the Bill and Melinda Gates funded C4 Rice Project. This project seeks to install the biochemical and anatomical attributes necessary to support C4 photosynthesis in the C3 crop rice. Here we address the advances made thus far in installing the biochemical pathway and some of the key targets yet to be reached.

Identification of Two Flip-Over Genes in Grass Family as Potential Signature of C4 Photosynthesis Evolution

In flowering plants, C4 photosynthesis is superior to C3 type in carbon fixation efficiency and adaptation to extreme environmental conditions, but the mechanisms behind the assembly of C4 machinery remain elusive. This study attempts to dissect the evolutionary divergence from C3 to C4 photosynthesis in five photosynthetic model plants from the grass family, using a combined comparative transcriptomics and deep learning technology. By examining and comparing gene expression levels in bundle sheath and mesophyll cells of five model plants, we identified 16 differentially expressed signature genes showing cell-specific expression patterns in C3 and C4 plants. Among them, two showed distinctively opposite cell-specific expression patterns in C3 vs. C4 plants (named as FOGs). The in silico physicochemical analysis of the two FOGs illustrated that C3 homologous proteins of LHCA6 had low and stable pI values of ~6, while the pI values of LHCA6 homologs increased drastically in C4 plants Setaria viridis (7), Zea mays (8), and Sorghum bicolor (over 9), suggesting this protein may have different functions in C3 and C4 plants. Interestingly, based on pairwise protein sequence/structure similarities between each homologous FOG protein, one FOG PGRL1A showed local inconsistency between sequence similarity and structure similarity. To find more examples of the evolutionary characteristics of FOG proteins, we investigated the protein sequence/structure similarities of other FOGs (transcription factors) and found that FOG proteins have diversified incompatibility between sequence and structure similarities during grass family evolution. This raised an interesting question as to whether the sequence similarity is related to structure similarity during C4 photosynthesis evolution.

Mutations in NAKED-ENDOSPERM IDD genes reveal functional interactions with SCARECROW during leaf patterning in C4 grasses

Leaves comprise a number of different cell-types that are patterned in the context of either the epidermal or inner cell layers. In grass leaves, two distinct anatomies develop in the inner leaf tissues depending on whether the leaf carries out C3 or C4 photosynthesis. In both cases a series of parallel veins develops that extends from the leaf base to the tip but in ancestral C3 species veins are separated by a greater number of intervening mesophyll cells than in derived C4 species. We have previously demonstrated that the GRAS transcription factor SCARECROW (SCR) regulates the number of photosynthetic mesophyll cells that form between veins in the leaves of the C4 species maize, whereas it regulates the formation of stomata in the epidermal leaf layer in the C3 species rice. Here we show that SCR is required for inner leaf patterning in the C4 species Setaria viridis but in this species the presumed ancestral stomatal patterning role is also retained. Through a comparative mutant analysis between maize, setaria and rice we further demonstrate that loss of NAKED-ENDOSPERM (NKD) INDETERMINATE DOMAIN (IDD) protein function exacerbates loss of function scr phenotypes in the inner leaf tissues of maize and setaria but not rice. Specifically, in both setaria and maize, scr;nkd mutants exhibit an increased proportion of fused veins with no intervening mesophyll cells. Thus, combined action of SCR and NKD may control how many mesophyll cells are specified between veins in the leaves of C4 but not C3 grasses. Together our results provide insight into the evolution of cell patterning in grass leaves and demonstrate a novel patterning role for IDD genes in C4 leaves.

Simultaneous detection of miRNA and mRNA at the single-cell level in plant tissues

Detecting the simultaneous presence of a microRNA (miRNA) and a mRNA in a specific tissue can provide support for the prediction that the miRNA regulates the mRNA. Although two such methods have been developed for mammalian tissues, they have a low signal-noise ratio and/or poor resolution at the single-cell level. To overcome these drawbacks, we develop a method that uses sequence-specific miRNA-locked nucleic acid (LNA) and mRNA-LNA probes. Moreover, it augments the detection signal by rolling circle amplification, achieving a high signal-noise ratio at the single-cell level. Dot signals are counted for determining the expression levels of mRNA and miRNA molecules in specific cells. We show a high sequence specificity of our miRNA-LNA probe, revealing that it can discriminate single-base mismatches. Numerical quantification by our method is tested in transgenic rice lines with different gene expression levels. We conduct several applications. First, the spatial expression profiling of osa-miR156 and OsSPL12 in rice leaves reveals their specific expression in mesophyll cells. Second, studying rice and its mutant lines with our method reveals opposite expression patterns of miRNA and its target mRNA in tissues. Third, the dynamic expression profiles of ZmGRF8 and zma-miR396 during maize leaf development provide evidence that zma-miR396 regulates the preferential spatial expression of ZmGRF8 in bundle sheath cells. Finally, our method can be scaled up to simultaneously detect multiple miRNAs and mRNAs in a tissue. Thus, it is a sensitive and versatile technique for studying miRNA regulation of plant tissue development.

Gene family evolution and natural selection signatures in Datura spp. (Solanaceae)

Elucidating the diversification process of congeneric species makes it necessary to identify the factors promoting species variation and diversification. Comparative gene family analysis allows us to elucidate the evolutionary history of species by identifying common genetic/genomic mechanisms underlying species responses to biotic and abiotic environments at the genomic level. In this study, we analyzed the high-quality transcriptomes of four Datura species, D. inoxia, D. pruinosa, D. stramonium, and D. wrightii. We performed a thorough comparative gene family analysis to infer the role of selection in molecular variation, changes in protein physicochemical properties, and gain/loss of genes during their diversification processes. The results revealed common and species-specific signals of positive selection, physicochemical divergence and/or expansion of metabolic genes (e.g., transferases and oxidoreductases) associated with terpene and tropane metabolism and some resistance genes (R genes). The gene family analysis presented here is a valuable tool for understanding the genome evolution of economically and ecologically significant taxa such as the Solanaceae family.

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.

Network and epigenetic characterization of subsets of genes specifically expressed in maize bundle sheath cells

Bundle sheath (BS) cells exhibit dramatically structural differences and functional variations at physiological, biochemical and epigenetic levels as compared to mesophyll (M) cells in maize. The regulatory mechanisms controlling functional divergences between M and BS have been extensively investigated. However, BS cell-related regulatory networks are still not completely characterized. To address this, we herein conducted WGCNA-related co-expression assays using bulk M and BS cell RNA-seq data sets and identified a module containing 384 genes highly expressed in BS cells (including 20 hub TFs) instead of M cells. According to the hub TF centered regulatory network, we found that Dof22 and Dof30 might act as key regulators in the regulation of expression of BS-specific genes, and several MYB TFs exhibited a high collaboration with Dof TFs. By comparing the enrichment levels of histone modifications, we found that genes in the aforementioned module were more enriched with histone acetylation as compared to other BS-enriched DEGs with similar expression levels. Moreover, we found that a subset of genes functioning in photosynthesis, protein auto processing and enzymatic activities were significantly enriched with broad H3K4me3. Thus, our study provides evidence showing that regulatory network and histone modifications may play vital roles in the regulation of a subset of genes with important functions in BS cells.

Gene co-expression reveals the modularity and integration of C4 and CAM in Portulaca

C4 photosynthesis and Crassulacean acid metabolism (CAM) have been considered as largely independent adaptations despite sharing key biochemical modules. Portulaca is a geographically widespread clade of over 100 annual and perennial angiosperm species that primarily use C4 but facultatively exhibit CAM when drought stressed, a photosynthetic system known as C4 + CAM. It has been hypothesized that C4 + CAM is rare because of pleiotropic constraints, but these have not been deeply explored. We generated a chromosome-level genome assembly of Portulaca amilis and sampled mRNA from P. amilis and Portulaca oleracea during CAM induction. Gene co-expression network analyses identified C4 and CAM gene modules shared and unique to both Portulaca species. A conserved CAM module linked phosphoenolpyruvate carboxylase to starch turnover during the day-night transition and was enriched in circadian clock regulatory motifs in the P. amilis genome. Preservation of this co-expression module regardless of water status suggests that Portulaca constitutively operate a weak CAM cycle that is transcriptionally and posttranscriptionally upregulated during drought. C4 and CAM mostly used mutually exclusive genes for primary carbon fixation, and it is likely that nocturnal CAM malate stores are shuttled into diurnal C4 decarboxylation pathways, but we found evidence that metabolite cycling may occur at low levels. C4 likely evolved in Portulaca through co-option of redundant genes and integration of the diurnal portion of CAM. Thus, the ancestral CAM system did not strongly constrain C4 evolution because photosynthetic gene networks are not co-regulated for both daytime and nighttime functions.

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.

Maize Golden2-like transcription factors boost rice chloroplast development, photosynthesis, and grain yield

Chloroplasts are the sites for photosynthesis, and two Golden2-like factors act as transcriptional activators of chloroplast development in rice (Oryza sativa L.) and maize (Zea mays L.). Rice OsGLK1 and OsGLK2 are orthologous to maize ZmGLK1 (ZmG1) and ZmGLK2 (ZmG2), respectively. However, while rice OsGLK1 and OsGLK2 act redundantly to regulate chloroplast development in mesophyll cells, maize ZmG1 and ZmG2 are functionally specialized and expressed in different cell-specific manners. To boost rice chloroplast development and photosynthesis, we generated transgenic rice plants overexpressing ZmG1 and ZmG2, individually or simultaneously, with constitutive promoters (pZmUbi::ZmG1 and p35S::ZmG2) or maize promoters (pZmG1::ZmG1, pZmG2::ZmG2, and pZmG1::ZmG1/pZmG2::ZmG2). Both ZmG1 and ZmG2 genes were highly expressed in transgenic rice leaves. Moreover, ZmG1 and ZmG2 showed coordinated expression in pZmG1::ZmG1/pZmG2::ZmG2 plants. All Golden2-like (GLK) transgenic plants had higher chlorophyll and protein contents, Rubisco activities and photosynthetic rates per unit leaf area in flag leaves. However, the highest grain yields occurred when maize promoters were used; pZmG1::ZmG1, pZmG2::ZmG2, and pZmG1::ZmG1/pZmG2::ZmG2 transgenic plants showed increases in grain yield by 51%, 47%, and 70%, respectively. In contrast, the pZmUbi::ZmG1 plant produced smaller seeds without yield increases. Transcriptome analysis indicated that maize GLKs act as master regulators promoting the expression of both photosynthesis-related and stress-responsive regulatory genes in both rice shoot and root. Thus, by promoting these important functions under the control of their own promoters, maize GLK1 and GLK2 genes together dramatically improved rice photosynthetic performance and productivity. A similar approach can potentially improve the productivity of many other crops.

The C4 cycle and beyond: diverse metabolic adaptations accompany dual-cell photosynthetic functions in Setaria

C4 photosynthesis is typically characterized by the spatial compartmentalization of the photosynthetic reactions into mesophyll (M) and bundle sheath (BS) cells. Initial carbon fixation within M cells gives rise to C4 acids, which are transported to the BS cells. There, C4 acids are decarboxylated so that the resulting CO2 is incorporated into the Calvin cycle. This work is focused on the study of Setaria viridis, a C4 model plant, closely related to several major feed and bioenergy grasses. First, we performed the heterologous expression and biochemical characterization of Setaria isoforms for chloroplastic NADP-malic enzyme (NADP-ME) and mitochondrial NAD-malic enzyme (NAD-ME). The kinetic parameters obtained agree with a major role for NADP-ME in the decarboxylation of the C4 acid malate in the chloroplasts of BS cells. In addition, mitochondria-located NAD-ME showed regulatory properties that could be important in the context of the operation of the C4 carbon shuttle. Secondly, we compared the proteomes of M and BS compartments and found 825 differentially accumulated proteins that could support different metabolic scenarios. Most interestingly, we found evidence of metabolic strategies to insulate the C4 core avoiding the leakage of intermediates by either up-regulation or down-regulation of chloroplastic, mitochondrial, and peroxisomal proteins. Overall, the results presented in this work provide novel data concerning the complexity of C4 metabolism, uncovering future lines of research that will undoubtedly contribute to the expansion of knowledge on this topic.

Ground tissue circuitry regulates organ complexity in maize and Setaria

[Figure: see text].

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.

Identification of the Golden-2-like transcription factors gene family in Gossypium hirsutum

Background

Golden2-Like (GLK) transcription factors are a type of transcriptional regulator in plants. They play a pivotal role in the plant physiological activity process and abiotic stress response.

Methods

In this study, the potential function of GLK family genes in Gossypium hirsutum was studied based on genomic identification, phylogenetic analysis, chromosome mapping and cis-regulatory elements prediction. Gene expression of nine key genes were analyzed by qRT-PCR experiments.

Results

Herein, we identified a total of 146 GhGLK genes in Gossypium hirsutum, which were unevenly distributed on each of the chromosomes. There were significant differences in the number and location of genes between the At sub-genome and the Dt sub-genome. According to the phylogenetic analysis, they were divided into ten subgroups, each of which had very similar number and structure of exons and introns. Some cis-regulatory elements were identified through promoter analysis, including five types of elements related to abiotic stress response, five types of elements related to phytohormone and five types of elements involved in growth and development. Based on public transcriptome data analysis, we identified nine key GhGLKs involved in salt, cold, and drought stress. The qRT-PCR results showed that these genes had different expression patterns under these stress conditions, suggesting that GhGLK genes played an important role in abiotic stress response. This study laid a theoretical foundation for the screening and functional verification of genes related to stress resistance of GLK gene family in cotton.

Genomic insights into the fast growth of paulownias and the formation of Paulownia witches' broom

Paulownias are among the fastest growing trees in the world, but they often suffer tremendous loss of wood production due to infection by Paulownia witches' broom (PaWB) phytoplasmas. In this study, we have sequenced and assembled a high-quality nuclear genome of Paulownia fortunei, a commonly cultivated paulownia species. The assembled genome of P. fortunei is 511.6 Mb in size, with 93.2% of its sequences anchored to 20 pseudo-chromosomes, and it contains 31 985 protein-coding genes. Phylogenomic analyses show that the family Paulowniaceae is sister to a clade composed of Phrymaceae and Orobanchaceae. Higher photosynthetic efficiency is achieved by integrating C3 photosynthesis and the crassulacean acid metabolism pathway, which may contribute to the extremely fast growth habit of paulownia trees. Comparative transcriptome analyses reveal modules related to cambial growth and development, photosynthesis, and defense responses. Additional genome sequencing of PaWB phytoplasma, combined with functional analyses, indicates that the effector PaWB-SAP54 interacts directly with Paulownia PfSPLa, which in turn causes the degradation of PfSPLa by the ubiquitin-mediated pathway and leads to the formation of witches' broom. Taken together, these results provide significant insights into the biology of paulownias and the regulatory mechanism for the formation of PaWB.

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.

Ribosome profiling elucidates differential gene expression in bundle sheath and mesophyll cells in maize

The efficiencies offered by C4 photosynthesis have motivated efforts to understand its biochemical, genetic, and developmental basis. Reactions underlying C4 traits in most C4 plants are partitioned between two cell types, bundle sheath (BS), and mesophyll (M) cells. RNA-seq has been used to catalog differential gene expression in BS and M cells in maize (Zea mays) and several other C4 species. However, the contribution of translational control to maintaining the distinct proteomes of BS and M cells has not been addressed. In this study, we used ribosome profiling and RNA-seq to describe translatomes, translational efficiencies, and microRNA abundance in BS- and M-enriched fractions of maize seedling leaves. A conservative interpretation of our data revealed 182 genes exhibiting cell type-dependent differences in translational efficiency, 31 of which encode proteins with core roles in C4 photosynthesis. Our results suggest that non-AUG start codons are used preferentially in upstream open reading frames of BS cells, revealed mRNA sequence motifs that correlate with cell type-dependent translation, and identified potential translational regulators that are differentially expressed. In addition, our data expand the set of genes known to be differentially expressed in BS and M cells, including genes encoding transcription factors and microRNAs. These data add to the resources for understanding the evolutionary and developmental basis of C4 photosynthesis and for its engineering into C3 crops.

Unoccupied aerial systems discovered overlooked loci capturing the variation of entire growing period in maize

Traditional phenotyping methods, coupled with genetic mapping in segregating populations, have identified loci governing complex traits in many crops. Unoccupied aerial systems (UAS)-based phenotyping has helped to reveal a more novel and dynamic relationship between time-specific associated loci with complex traits previously unable to be evaluated. Over 1,500 maize (Zea mays L.) hybrid row plots containing 280 different replicated maize hybrids from the Genomes to Fields (G2F) project were evaluated agronomically and using UAS in 2017. Weekly UAS flights captured variation in plant heights during the growing season under three different management conditions each year: optimal planting with irrigation (G2FI), optimal dryland planting without irrigation (G2FD), and a stressed late planting (G2LA). Plant height of different flights were ranked based on importance for yield using a random forest (RF) algorithm. Plant heights captured by early flights in G2FI trials had higher importance (based on Gini scores) for predicting maize grain yield (GY) but also higher accuracies in genomic predictions which fluctuated for G2FD (-0.06∼0.73), G2FI (0.33∼0.76), and G2LA (0.26∼0.78) trials. A genome-wide association analysis discovered 52 significant single nucleotide polymorphisms (SNPs), seven were found consistently in more than one flights or trial; 45 were flight or trial specific. Total cumulative marker effects for each chromosome's contributions to plant height also changed depending on flight. Using UAS phenotyping, this study showed that many candidate genes putatively play a role in the regulation of plant architecture even in relatively early stages of maize growth and development.

The bundle sheath of rice is conditioned to play an active role in water transport as well as sulfur assimilation and jasmonic acid synthesis

Leaves comprise multiple cell types but our knowledge of the patterns of gene expression that underpin their functional specialization is fragmentary. Our understanding and ability to undertake the rational redesign of these cells is therefore limited. We aimed to identify genes associated with the incompletely understood bundle sheath of C₃ plants, which represents a key target associated with engineering traits such as C₄ photosynthesis into Oryza sativa (rice). To better understand the veins, bundle sheath and mesophyll cells of rice, we used laser capture microdissection followed by deep sequencing. Gene expression of the mesophyll is conditioned to allow coenzyme metabolism and redox homeostasis, as well as photosynthesis. In contrast, the bundle sheath is specialized in water transport, sulphur assimilation and jasmonic acid biosynthesis. Despite the small chloroplast compartment of bundle sheath cells, substantial photosynthesis gene expression was detected. These patterns of gene expression were not associated with the presence or absence of specific transcription factors in each cell type, but were instead associated with gradients in expression across the leaf. Comparative analysis with C₃ Arabidopsis identified a small gene set preferentially expressed in the bundle sheath cells of both species. This gene set included genes encoding transcription factors from 14 orthogroups and proteins allowing water transport, sulphate assimilation and jasmonic acid synthesis. The most parsimonious explanation for our findings is that bundle sheath cells from the last common ancestor of rice and Arabidopsis were specialized in this manner, and as the species diverged these patterns of gene expression have been maintained.

Comparison of Transcriptional Response of C3 and C4 Plants to Drought Stress Using Meta-Analysis and Systems Biology Approach

Drought stress affects a range of plant processes. It is still not well-known how C3 and C4 plants respond to drought. Here, we used a combination of meta-analysis and network analysis to compare the transcriptional responses of Oryza sativa (rice), a C3 plant, and Zea mays (maize), a C4 plant, to drought stress. The findings showed that drought stress changes the expression of genes and affects different mechanisms in the C3 and C4 plants. We identified several genes that were differentially expressed genes (DEGs) under stress conditions in both species, most of which are associated with photosynthesis, molecule metabolic process, and response to stress. Additionally, we observed that many DEGs physically located within the quantitative trait locus regions are associated with C isotope signature (d13C), photosynthetic gas exchange, and root characteristics traits. Through the gene co-expression and differential co-expression network methods, we identified sets of genes with similar and different behaviors among C3 and C4 plants during drought stress. This result indicates that mitogen-activated protein kinases (MAPK) signaling pathway plays an important part in the differences between the C3 and C4 species. The present study provides a better understanding of the mechanisms underlying the response of C3 and C4 plants to drought stress, which may useful for engineering drought tolerance in plants.

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

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

A cis-regulatory atlas in maize at single-cell resolution

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

A cis-regulatory atlas in maize at single-cell resolution

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

Evidence for phloem loading via the abaxial bundle sheath cells in maize leaves

Leaves are asymmetric, with different functions for adaxial and abaxial tissue. The bundle sheath (BS) of C3 barley (Hordeum vulgare) is dorsoventrally differentiated into three types of cells: adaxial structural, lateral S-type, and abaxial L-type BS cells. Based on plasmodesmatal connections between S-type cells and mestome sheath (parenchymatous cell layer below bundle sheath), S-type cells likely transfer assimilates toward the phloem. Here, we used single-cell RNA sequencing to investigate BS differentiation in C4 maize (Zea mays L.) plants. Abaxial BS (abBS) cells of rank-2 intermediate veins specifically expressed three SWEET sucrose uniporters (SWEET13a, b, and c) and UmamiT amino acid efflux transporters. SWEET13a, b, c mRNAs were also detected in the phloem parenchyma (PP). We show that maize has acquired a mechanism for phloem loading in which abBS cells provide the main route for apoplasmic sucrose transfer toward the phloem. This putative route predominates in veins responsible for phloem loading (rank-2 intermediate), whereas rank-1 intermediate and major veins export sucrose from the PP adjacent to the sieve element companion cell complex, as in Arabidopsis thaliana. We surmise that abBS identity is subject to dorsoventral patterning and has components of PP identity. These observations provide insights into the unique transport-specific properties of abBS cells and support a modification to the canonical phloem loading pathway in maize.

Single-cell RNA sequencing of developing maize ears facilitates functional analysis and trait candidate gene discovery

Crop productivity depends on activity of meristems that produce optimized plant architectures, including that of the maize ear. A comprehensive understanding of development requires insight into the full diversity of cell types and developmental domains and the gene networks required to specify them. Until now, these were identified primarily by morphology and insights from classical genetics, which are limited by genetic redundancy and pleiotropy. Here, we investigated the transcriptional profiles of 12,525 single cells from developing maize ears. The resulting developmental atlas provides a single-cell RNA sequencing (scRNA-seq) map of an inflorescence. We validated our results by mRNA in situ hybridization and by fluorescence-activated cell sorting (FACS) RNA-seq, and we show how these data may facilitate genetic studies by predicting genetic redundancy, integrating transcriptional networks, and identifying candidate genes associated with crop yield traits.

Cereal asparagine synthetase genes

Asparagine synthetase catalyses the transfer of an amino group from glutamine to aspartate to form glutamate and asparagine. The accumulation of free (nonprotein) asparagine in crops has implications for food safety because free asparagine is the precursor for acrylamide, a carcinogenic contaminant that forms during high-temperature cooking and processing. Here we review publicly available genome data for asparagine synthetase genes from species of the Pooideae subfamily, including bread wheat and related wheat species (Triticum and Aegilops spp.), barley (Hordeum vulgare) and rye (Secale cereale) of the Triticeae tribe. Also from the Pooideae subfamily: brachypodium (Brachypodium dIstachyon) of the Brachypodiae tribe. More diverse species are also included, comprising sorghum (Sorghum bicolor) and maize (Zea mays) of the Panicoideae subfamily and rice (Oryza sativa) of the Ehrhartoideae subfamily. The asparagine synthetase gene families of the Triticeae species each comprise five genes per genome, with the genes assigned to four groups: 1, 2, 3 (subdivided into 3.1 and 3.2) and 4. Each species has a single gene per genome in each group, except that some bread wheat varieties (genomes AABBDD) and emmer wheat (Triticum dicoccoides; genomes AABB) lack a group 2 gene in the B genome. This raises questions about the ancestry of cultivated pasta wheat and the B genome donor of bread wheat, suggesting that the hybridisation event that gave rise to hexaploid bread wheat occurred more than once. In phylogenetic analyses, genes from the other species cluster with the Triticeae genes, but brachypodium, sorghum and maize lack a group 2 gene, while rice has only two genes, one group 3 and one group 4. This means that TaASN2, the most highly expressed asparagine synthetase gene in wheat grain, has no equivalent in maize, rice, sorghum or brachypodium. An evolutionary pathway is proposed in which a series of gene duplications gave rise to the five genes found in modern Triticeae species.

Transcriptome analysis of a novel maize bsd C4 mutant using RNA-seq

The C4 plants like maize are at an advantage because they exhibit higher carbon conversion efficiency than C3 during photosynthesis. Using the evaluation of photosynthetic phenotypes and subcellular structure, and high-quality transcriptome analysis for four types of leaves from a bundle sheath defective maize mutant (bsd) and 501 wild line, the key target genes, important transcription factors, and specific pathways were obtained, which may regulate the C-concentrating mechanisms and antioxidant protection of the photosynthetic system, gibberellin signaling, ribosome editing, glycolysis, and chlorophyll biosynthesis. Based on these target genes, a novel network with photosynthetic transformation efficiency with oxidative decarboxylation and ribosome regulation was filtered innovatively by Cytoscape, which adds to our understanding for high-efficiency C-fixation and its genetic improvement in C3 and C4 plants.

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

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

Regulation and Evolution of C4 Photosynthesis

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

Subdivision of Light Signaling Networks Contributes to Partitioning of C4 Photosynthesis

Plants coordinate the expression of photosynthesis-related genes in response to growth and environmental changes. In species that conduct two-cell C4 photosynthesis, expression of photosynthesis genes is partitioned such that leaf mesophyll and bundle sheath cells accumulate different components of the photosynthetic pathway. The identities of the regulatory networks that facilitate this partitioning are unknown. Here, we show that differences in light perception between mesophyll and bundle sheath cells facilitate differential regulation and accumulation of photosynthesis gene transcripts in the C4 crop maize (Zea mays). Key components of the photosynthesis gene regulatory network differentially accumulated between mesophyll and bundle sheath cells, indicative of differential network activity across cell types. We further show that blue (but not red) light is necessary and sufficient to activate photosystem II assembly in mesophyll cells in etiolated maize. Finally, we demonstrate that 61% of all light-induced mesophyll and bundle sheath genes were induced only by blue light or only by red light, but not both. These findings provide evidence that subdivision of light signaling networks is a component of cellular partitioning of C4 photosynthesis in maize.

Insights Into the Regulation of the Expression Pattern of Calvin-Benson-Bassham Cycle Enzymes in C3 and C4 Grasses

C₄ photosynthesis is characterized by the compartmentalization of the processes of atmospheric uptake of CO₂ and its conversion into carbohydrate between mesophyll and bundle-sheath cells. As a result, most of the enzymes participating in the Calvin-Benson-Bassham (CBB) cycle, including RubisCO, are highly expressed in bundle-sheath cells. There is evidence that changes in the regulatory sequences of RubisCO contribute to its bundle-sheath-specific expression, however, little is known about how the spatial-expression pattern of other CBB cycle enzymes is regulated. In this study, we use a computational approach to scan for transcription factor binding sites in the regulatory regions of the genes encoding CBB cycle enzymes, SBPase, FBPase, PRK, and GAPDH-B, of C₃ and C₄ grasses. We identified potential cis-regulatory elements present in each of the genes studied here, regardless of the photosynthetic path used by the plant. The trans-acting factors that bind these elements have been validated in A. thaliana and might regulate the expression of the genes encoding CBB cycle enzymes. In addition, we also found C₄-specific transcription factor binding sites in the genes encoding CBB cycle enzymes that could potentially contribute to the pathway-specific regulation of gene expression. These results provide a foundation for the functional analysis of the differences in regulation of genes encoding CBB cycle enzymes between C₃ and C₄ grasses.

Expression profiles of cell-wall related genes vary broadly between two common maize inbreds during stem development

Background

The cellular machinery for cell wall synthesis and metabolism is encoded by members of large multi-gene families. Maize is both a genetic model for grass species and a potential source of lignocellulosic biomass from crop residues. Genetic improvement of maize for its utility as a bioenergy feedstock depends on identification of the specific gene family members expressed during secondary wall development in stems.

Results

High-throughput sequencing of transcripts expressed in developing rind tissues of stem internodes provided a comprehensive inventory of cell wall-related genes in maize (Zea mays, cultivar B73). Of 1239 of these genes, 854 were expressed among the internodes at ≥95 reads per 20 M, and 693 of them at ≥500 reads per 20 M. Grasses have cell wall compositions distinct from non-commelinid species; only one-quarter of maize cell wall-related genes expressed in stems were putatively orthologous with those of the eudicot Arabidopsis. Using a slope-metric algorithm, five distinct patterns for sub-sets of co-expressed genes were defined across a time course of stem development. For the subset of genes associated with secondary wall formation, fifteen sequence motifs were found in promoter regions. The same members of gene families were often expressed in two maize inbreds, B73 and Mo17, but levels of gene expression between them varied, with 30% of all genes exhibiting at least a 5-fold difference at any stage. Although presence-absence and copy-number variation might account for much of these differences, fold-changes of expression of a CADa and a FLA11 gene were attributed to polymorphisms in promoter response elements.

Conclusions

Large genetic variation in maize as a species precludes the extrapolation of cell wall-related gene expression networks even from one common inbred line to another. Elucidation of genotype-specific expression patterns and their regulatory controls will be needed for association panels of inbreds and landraces to fully exploit genetic variation in maize and other bioenergy grass species.

Genome-Wide Transcription Factor Binding in Leaves from C3 and C4 Grasses

The majority of plants use C3 photosynthesis, but over 60 independent lineages of angiosperms have evolved the C4 pathway. In most C4 species, photosynthesis gene expression is compartmented between mesophyll and bundle-sheath cells. We performed DNaseI sequencing to identify genome-wide profiles of transcription factor binding in leaves of the C4 grasses Zea mays, Sorghum bicolor, and Setaria italica as well as C3 Brachypodium distachyon In C4 species, while bundle-sheath strands and whole leaves shared similarity in the broad regions of DNA accessible to transcription factors, the short sequences bound varied. Transcription factor binding was prevalent in gene bodies as well as promoters, and many of these sites could represent duons that influence gene regulation in addition to amino acid sequence. Although globally there was little correlation between any individual DNaseI footprint and cell-specific gene expression, within individual species transcription factor binding to the same motifs in multiple genes provided evidence for shared mechanisms governing C4 photosynthesis gene expression. Furthermore, interspecific comparisons identified a small number of highly conserved transcription factor binding sites associated with leaves from species that diverged around 60 million years ago. These data therefore provide insight into the architecture associated with C4 photosynthesis gene expression in particular and characteristics of transcription factor binding in cereal crops in general.

Redundant SCARECROW genes pattern distinct cell layers in roots and leaves of maize

The highly efficient C₄ photosynthetic pathway is facilitated by ‘Kranz’ leaf anatomy. In Kranz leaves, closely spaced veins are encircled by concentric layers of photosynthetic bundle sheath (inner) and mesophyll (outer) cells. Here, we demonstrate that, in the C₄ monocot maize, Kranz patterning is regulated by redundant function of SCARECROW 1 (ZmSCR1) and a previously uncharacterized homeologue: ZmSCR1h. ZmSCR1 and ZmSCR1h transcripts accumulate in ground meristem cells of developing leaf primordia and in Zmscr1;Zmscr1h mutant leaves, most veins are separated by one rather than two mesophyll cells; many veins have sclerenchyma above and/or below instead of mesophyll cells; and supernumerary bundle sheath cells develop. The mutant defects are unified by compromised mesophyll cell development. In addition to Kranz defects, Zmscr1;Zmscr1h mutants fail to form an organized endodermal layer in the root. Collectively, these data indicate that ZmSCR1 and ZmSCR1h redundantly regulate cell-type patterning in both the leaves and roots of maize. Leaf and root pathways are distinguished, however, by the cell layer in which they operate – mesophyll at a two-cell distance from leaf veins versus endodermis immediately adjacent to root vasculature.

Summary: Two duplicated maize SCARECROW genes control the development of the endodermis in roots and the mesophyll in leaves.

Decarboxylation mechanisms of C4 photosynthesis in Saccharum spp.: increased PEPCK activity under water-limiting conditions

BACKGROUND

C4 plants have been classified into three subtypes based on the enzymes used to decarboxylate C4 acids in the bundle sheath cells (NADP-ME, NAD-ME and PEPCK pathways). Evidences indicate that, depending on environmental factors, C4 plants may exhibit a certain degree of flexibility in the use of the decarboxylation mechanisms. In this context, the objective was to extend the knowledge on the degree of flexibility between the pathways of decarboxylation in sugarcane, a NADP-ME species, at different levels of water deficit.

RESULTS

An experiment was carried out with two cultivars - RB92579 (tolerant to water deficit) and SP80-3280 (susceptible to water deficit) subjected to moderate level (- 1.5 to - 1.8 MPa), severe level (below - 2.0 MPa) and recovery (48 h after rehydration) and changes in the activities of the enzymes involved in the three C4 mechanisms and in gene expression were investigated. Our results showed that sugarcane uses the PEPCK pathway as a decarboxylation mechanism in addition to the NADP-ME, which was more evident under water deficit conditions for both cultivars.

CONCLUSIONS

The results obtained here, show that sugarcane increases the use of the PEPCK pathway as a decarboxylation mechanism, in addition to the NADP-ME pathway, under conditions of water deficit, particularly in the tolerant cultivar.

Comparative transcriptomics method to infer gene coexpression networks and its applications to maize and rice leaf transcriptomes

Time-series transcriptomes of a biological process obtained under different conditions are useful for identifying the regulators of the process and their regulatory networks. However, such data are 3D (gene expression, time, and condition), and there is currently no method that can deal with their full complexity. Here, we developed a method that avoids time-point alignment and normalization between conditions. We applied it to analyze time-series transcriptomes of developing maize leaves under light-dark cycles and under total darkness and obtained eight time-ordered gene coexpression networks (TO-GCNs), which can be used to predict upstream regulators of any genes in the GCNs. One of the eight TO-GCNs is light-independent and likely includes all genes involved in the development of Kranz anatomy, which is a structure crucial for the high efficiency of photosynthesis in C₄ plants. Using this TO-GCN, we predicted and experimentally validated a regulatory cascade upstream of SHORTROOT1, a key Kranz anatomy regulator. Moreover, we applied the method to compare transcriptomes from maize and rice leaf segments and identified regulators of maize C₄ enzyme genes and RUBISCO SMALL SUBUNIT2 Our study provides not only a powerful method but also novel insights into the regulatory networks underlying Kranz anatomy development and C₄ photosynthesis.

The role of alanine and aspartate aminotransferases in C4 photosynthesis

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

Evolution of the Cell Wall Gene Families of Grasses

Grasses and related commelinid monocot species synthesize cell walls distinct in composition from other angiosperm species. With few exceptions, the genomes of all angiosperms contain the genes that encode the enzymes for synthesis of all cell-wall polysaccharide, phenylpropanoid, and protein constituents known in vascular plants. RNA-seq analysis of transcripts expressed during development of the upper and lower internodes of maize (Zea mays) stem captured the expression of cell-wall-related genes associated with primary or secondary wall formation. High levels of transcript abundances were not confined to genes associated with the distinct walls of grasses but also of those associated with xyloglucan and pectin synthesis. Combined with proteomics data to confirm that expressed genes are translated, we propose that the distinctive cell-wall composition of grasses results from sorting downstream from their sites of synthesis in the Golgi apparatus and hydrolysis of the uncharacteristic polysaccharides and not from differential expression of synthases of grass-specific polysaccharides.

Understanding the Genetic Basis of C4 Kranz Anatomy with a View to Engineering C3 Crops

One of the most remarkable examples of convergent evolution is the transition from C₃ to C₄ photosynthesis, an event that occurred on over 60 independent occasions. The evolution of C₄ is particularly noteworthy because of the complexity of the developmental and metabolic changes that took place. In most cases, compartmentalized metabolic reactions were facilitated by the development of a distinct leaf anatomy known as Kranz. C₄ Kranz anatomy differs from ancestral C₃ anatomy with respect to vein spacing patterns across the leaf, cell-type specification around veins, and cell-specific organelle function. Here we review our current understanding of how Kranz anatomy evolved and how it develops, with a focus on studies that are dissecting the underlying genetic mechanisms. This research field has gained prominence in recent years because understanding the genetic regulation of Kranz may enable the C₃-to-C₄ transition to be engineered, an endeavor that would significantly enhance crop productivity.

Strategies and tools to improve crop productivity by targeting photosynthesis

Crop productivity needs to substantially increase to meet global food and feed demand for a rapidly growing world population. Agricultural technology developers are pursuing a variety of approaches based on both traditional technologies such as genetic improvement, pest control and mechanization as well as new technologies such as genomics, gene manipulation and environmental modelling to develop crops that are capable of meeting growing demand. Photosynthesis is a key biochemical process that, many suggest, is not yet optimized for industrial agriculture or the modern global environment. We are interested in identifying control points in maize photoassimilation that are amenable to gene manipulation to improve overall productivity. Our approach encompasses: developing and using novel gene discovery techniques, translating our discoveries into traits and evaluating each trait in a stepwise manner that reflects a modern production environment. Our aim is to provide step change advancement in overall crop productivity and deliver this new technology into the hands of growers.This article is part of the themed issue 'Enhancing photosynthesis in crop plants: targets for improvement'.