Synergistic Binding of bHLH Transcription Factors to the Promoter of the Maize NADP-ME Gene Used in C4 Photosynthesis Is Based on an Ancient Code Found in the Ancestral C3 State

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.

Is winter coming? Impact of the changing climate on plant responses to cold temperature

Climate change is causing alterations in annual temperature regimes worldwide. Important aspects of this include the reduction of winter chilling temperatures as well as the occurrence of unpredicted frosts, both significantly affecting plant growth and yields. Recent studies advanced the knowledge of the mechanisms underlying cold responses and tolerance in the model plant Arabidopsis thaliana. However, how these cold-responsive pathways will readjust to ongoing seasonal temperature variation caused by global warming remains an open question. In this review, we highlight the plant developmental programmes that depend on cold temperature. We focus on the molecular mechanisms that plants have evolved to adjust their development and stress responses upon exposure to cold. Covering both genetic and epigenetic aspects, we present the latest insights into how alternative splicing, noncoding RNAs and the formation of biomolecular condensates play key roles in the regulation of cold responses. We conclude by commenting on attractive targets to accelerate the breeding of increased cold tolerance, bringing up biotechnological tools that might assist in overcoming current limitations. Our aim is to guide the reflection on the current agricultural challenges imposed by a changing climate and to provide useful information for improving plant resilience to unpredictable cold regimes.

Differential chromatin binding preference is the result of the neo-functionalization of the TB1 clade of TCP transcription factors in grasses

The understanding of neo-functionalization of plant transcription factors (TFs) after gene duplication has been extensively focused on changes in protein-protein interactions, the expression pattern of TFs, or the variation of cis-elements bound by TFs. Yet, the main molecular role of a TF, that is, its specific chromatin binding for the direct regulation of target gene expression, continues to be mostly overlooked. Here, we studied the TB1 clade of the TEOSINTE BRANCHED 1, CYCLOIDEA, PROLIFERATING CELL FACTORS (TCP) TF family within the grasses (Poaceae). We identified an Asp/Gly amino acid replacement within the TCP domain, originated within a paralog TIG1 clade exclusive for grasses. The heterologous expression of Zea mays TB1 and its two paralogs BAD1 and TIG1 in Arabidopsis mutant plants lacking the TB1 ortholog BRC1 revealed distinct functions in plant development. Notably, the Gly acquired in the TIG1 clade does not impair TF homodimerization and heterodimerization, while it modulates chromatin binding preferences. We found that in vivo TF recognition of target promoters depends on this Asp/Gly mutation and directly impacts downstream gene expression and subsequent plant development. These results provided new insights into how natural selection fine-tunes gene expression regulation after duplication of TFs to define plant architecture.

Epi-Brassinolide Regulates ZmC4 NADP-ME Expression through the Transcription Factors ZmbHLH157 and ZmNF-YC2

Maize is a main food and feed crop with great production potential and high economic benefits. Improving its photosynthesis efficiency is crucial for increasing yield. Maize photosynthesis occurs mainly through the C4 pathway, and NADP-ME (NADP-malic enzyme) is a key enzyme in the photosynthetic carbon assimilation pathway of C4 plants. ZmC4-NADP-ME catalyzes the release of CO₂ from oxaloacetate into the Calvin cycle in the maize bundle sheath. Brassinosteroid (BL) can improve photosynthesis; however, its molecular mechanism of action remains unclear. In this study, transcriptome sequencing of maize seedlings treated with epi-brassinolide (EBL) showed that differentially expressed genes (DEGs) were significantly enriched in photosynthetic antenna proteins, porphyrin and chlorophyll metabolism, and photosynthesis pathways. The DEGs of C4-NADP-ME and pyruvate phosphate dikinase in the C4 pathway were significantly enriched in EBL treatment. Co-expression analysis showed that the transcription level of ZmNF-YC2 and ZmbHLH157 transcription factors was increased under EBL treatment and moderately positively correlated with ZmC4-NADP-ME. Transient overexpression of protoplasts revealed that ZmNF-YC2 and ZmbHLH157 activate C4-NADP-ME promoters. Further experiments showed ZmNF-YC2 and ZmbHLH157 transcription factor binding sites on the -1616 bp and -1118 bp ZmC4 NADP-ME promoter. ZmNF-YC2 and ZmbHLH157 were screened as candidate transcription factors mediating brassinosteroid hormone regulation of the ZmC4 NADP-ME gene. The results provide a theoretical basis for improving maize yield using BR hormones.

Overexpression of cytoplasmic C4 Flaveria bidentis carbonic anhydrase in C3 Arabidopsis thaliana increases amino acids, photosynthetic potential, and biomass

An important method to improve photosynthesis in C₃ crops, such as rice and wheat, is to transfer efficient C₄ characters to them. Here, cytosolic carbonic anhydrase (CA: βCA3) of the C₄ Flaveria bidentis (Fb) was overexpressed under the control of ³⁵ S promoter in Arabidopsis thaliana, a C₃ plant, to enhance its photosynthetic efficiency. Overexpression of CA resulted in a better supply of the substrate HCO3- for the endogenous phosphoenolpyruvate carboxylase in the cytosol of the overexpressers, and increased its activity for generating malate that feeds into the tricarboxylic acid cycle. This provided additional carbon skeleton for increased synthesis of amino acids aspartate, asparagine, glutamate, and glutamine. Increased amino acids contributed to higher protein content in the transgenics. Furthermore, expression of FbβCA3 in Arabidopsis led to a better growth due to expression of several genes leading to higher chlorophyll content, electron transport, and photosynthetic carbon assimilation in the transformants. Enhanced CO₂ assimilation resulted in increased sugar and starch content, and plant dry weight. In addition, transgenic plants had lower stomatal conductance, reduced transpiration rate, and higher water-use efficiency. These results, taken together, show that expression of C₄ CA in the cytosol of a C₃ plant can indeed improve its photosynthetic capacity with enhanced water-use efficiency.

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.

Wheat-Thinopyrum Substitution Lines Imprint Compensation Both From Recipients and Donors

Even frequently used in wheat breeding, we still have an insufficient understanding of the biology of the products via distant hybridization. In this study, a transcriptomic analysis was performed for six Triticum aestivum-Thinopyrum elongatum substitution lines in comparison with the host plants. All the six disomic substitution lines showed much stronger "transcriptomic-shock" occurred on alien genomes with 57.43-69.22% genes changed expression level but less on the recipient genome (2.19-8.97%). Genome-wide suppression of alien genes along chromosomes was observed with a high proportion of downregulated genes (39.69-48.21%). Oppositely, the wheat recipient showed genome-wide compensation with more upregulated genes, occurring on all chromosomes but not limited to the homeologous groups. Moreover, strong co-upregulation of the orthologs between wheat and Thinopyrum sub-genomes was enriched in photosynthesis with predicted chloroplastic localization, which indicates that the compensation happened not only on wheat host genomes but also on alien genomes.

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.

Insights into Regulation of C2 and C4 Photosynthesis in Amaranthaceae/Chenopodiaceae Using RNA-Seq

Amaranthaceae (incl. Chenopodiaceae) shows an immense diversity of C4 syndromes. More than 15 independent origins of C4 photosynthesis, and the largest number of C4 species in eudicots signify the importance of this angiosperm lineage in C4 evolution. Here, we conduct RNA-Seq followed by comparative transcriptome analysis of three species from Camphorosmeae representing related clades with different photosynthetic types: Threlkeldia diffusa (C3), Sedobassia sedoides (C2), and Bassia prostrata (C4). Results show that B. prostrata belongs to the NADP-ME type and core genes encoding for C4 cycle are significantly upregulated when compared with Sed. sedoides and T. diffusa. Sedobassia sedoides and B. prostrata share a number of upregulated C4-related genes; however, two C4 transporters (DIT and TPT) are found significantly upregulated only in Sed. sedoides. Combined analysis of transcription factors (TFs) of the closely related lineages (Camphorosmeae and Salsoleae) revealed that no C3-specific TFs are higher in C2 species compared with C4 species; instead, the C2 species show their own set of upregulated TFs. Taken together, our study indicates that the hypothesis of the C2 photosynthesis as a proxy towards C4 photosynthesis is questionable in Sed. sedoides and more in favour of an independent evolutionary stable state.

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

Spatial separation of the photosynthetic reactions is a key feature of C₄ metabolism. In most C₄ 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.

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.]

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.

Regulation and Evolution of C4 Photosynthesis

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

Physiological and Molecular Mechanism Involved in Cold Stress Tolerance in Plants

Previous studies have reported that low temperature (LT) constrains plant growth and restricts productivity in temperate regions. However, the underlying mechanisms are complex and not well understood. Over the past ten years, research on the process of adaptation and tolerance of plants during cold stress has been carried out. In molecular terms, researchers prioritize research into the field of the ICE-CBF-COR signaling pathway which is believed to be the important key to the cold acclimation process. Inducer of CBF Expression (ICE) is a pioneer of cold acclimation and plays a central role in C-repeat binding (CBF) cold induction. CBFs activate the expression of COR genes via binding to cis-elements in the promoter of COR genes. An ICE-CBF-COR signaling pathway activates the appropriate expression of downstream genes, which encodes osmoregulation substances. In this review, we summarize the recent progress of cold stress tolerance in plants from molecular and physiological perspectives and other factors, such as hormones, light, and circadian clock. Understanding the process of cold stress tolerance and the genes involved in the signaling network for cold stress is essential for improving plants, especially crops.

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.

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.

Key changes in gene expression identified for different stages of C4 evolution in Alloteropsis semialata

C4 photosynthesis is a complex trait that boosts productivity in tropical conditions. Compared with C3 species, the C4 state seems to require numerous novelties, but species comparisons can be confounded by long divergence times. Here, we exploit the photosynthetic diversity that exists within a single species, the grass Alloteropsis semialata, to detect changes in gene expression associated with different photosynthetic phenotypes. Phylogenetically informed comparative transcriptomics show that intermediates with a weak C4 cycle are separated from the C3 phenotype by increases in the expression of 58 genes (0.22% of genes expressed in the leaves), including those encoding just three core C4 enzymes: aspartate aminotransferase, phosphoenolpyruvate carboxykinase, and phosphoenolpyruvate carboxylase. The subsequent transition to full C4 physiology was accompanied by increases in another 15 genes (0.06%), including only the core C4 enzyme pyruvate orthophosphate dikinase. These changes probably created a rudimentary C4 physiology, and isolated populations subsequently improved this emerging C4 physiology, resulting in a patchwork of expression for some C4 accessory genes. Our work shows how C4 assembly in A. semialata happened in incremental steps, each requiring few alterations over the previous step. These create short bridges across adaptive landscapes that probably facilitated the recurrent origins of C4 photosynthesis through a gradual process of evolution.

The molecular evolution of C4 photosynthesis: opportunities for understanding and improving the world's most productive plants

C4 photosynthesis is a convergent evolutionary trait that enhances photosynthetic efficiency in a variety of environmental conditions. It has evolved repeatedly following a fall in atmospheric CO2 concentration such that there is up to a 30 million year difference in the amount of time that natural selection has had to improve C4 function between the oldest and youngest C4 lineages. This large difference in time, coupled with the phylogenetic distance between lineages, has resulted in a large disparity in anatomy, physiology, and biochemistry between extant C4 species. This review summarizes the myriad of molecular sequence changes that have been linked to the evolution of C4 photosynthesis. These range from single nucleotide changes to duplication of entire genes, and provide a roadmap for how natural selection has adapted enzymes and pathways for enhanced C4 function. Finally, this review discusses how this molecular diversity can provide opportunities for understanding and improving photosynthesis for multiple important C4 food, feed, and bioenergy crops.

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.

Genome-wide analysis of the basic Helix-Loop-Helix (bHLH) transcription factor family in maize

BACKGROUND

In plants, the basic helix-loop-helix (bHLH) transcription factors play key roles in diverse biological processes. Genome-wide comprehensive and systematic analyses of bHLH proteins have been well conducted in Arabidopsis, rice, tomato and other plant species. However, only few of bHLH family genes have been functional characterized in maize.

RESULTS

In this study, our genome-wide analysis identified 208 putative bHLH family proteins (ZmbHLH proteins) in maize (Zea mays). We classified these proteins into 18 subfamilies by comparing the ZmbHLHs with Arabidopsis thaliana bHLH proteins. Phylogenetic analysis, conserved protein motifs, and exon-intron patterns further supported the evolutionary relationships among these bHLH proteins. Genome distribution analysis found that the 208 ZmbHLH loci were located non-randomly on the ten maize chromosomes. Further, analysis of conserved cis-elements in the promoter regions, protein interaction networks, and expression patterns in roots, leaves, and seeds across developmental stages, suggested that bHLH family proteins in maize are probably involved in multiple physiological processes in plant growth and development.

CONCLUSION

We performed a genome-wide, systematic analysis of bHLH proteins in maize. This comprehensive analysis provides a useful resource that enables further investigation of the physiological roles and molecular functions of the ZmbHLH transcription factors.