Leaf rolling allows quantification of mRNA abundance in mesophyll cells of sorghum

Cell-type-specific transcriptomics uncovers spatial regulatory networks in bioenergy sorghum stems

Bioenergy sorghum is a low-input, drought-resilient, deep-rooting annual crop that has high biomass yield potential enabling the sustainable production of biofuels, biopower, and bioproducts. Bioenergy sorghum's 4-5 m stems account for ~80% of the harvested biomass. Stems accumulate high levels of sucrose that could be used to synthesize bioethanol and useful biopolymers if information about cell-type gene expression and regulation in stems was available to enable engineering. To obtain this information, laser capture microdissection was used to isolate and collect transcriptome profiles from five major cell types that are present in stems of the sweet sorghum Wray. Transcriptome analysis identified genes with cell-type-specific and cell-preferred expression patterns that reflect the distinct metabolic, transport, and regulatory functions of each cell type. Analysis of cell-type-specific gene regulatory networks (GRNs) revealed that unique transcription factor families contribute to distinct regulatory landscapes, where regulation is organized through various modes and identifiable network motifs. Cell-specific transcriptome data was combined with known secondary cell wall (SCW) networks to identify the GRNs that differentially activate SCW formation in vascular sclerenchyma and epidermal cells. The spatial transcriptomic dataset provides a valuable source of information about the function of different sorghum cell types and GRNs that will enable the engineering of bioenergy sorghum stems, and an interactive web application developed during this project will allow easy access and exploration of the data (https://mc-lab.shinyapps.io/lcm-dataset/).

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.

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.

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.

Upregulation of bundle sheath electron transport capacity under limiting light in C4 Setaria viridis

C₄ photosynthesis is a biochemical pathway that operates across mesophyll and bundle sheath (BS) cells to increase CO₂ concentration at the site of CO₂ fixation. C₄ plants benefit from high irradiance but their efficiency decreases under shade, causing a loss of productivity in crop canopies. We investigated shade acclimation responses of Setaria viridis, a model monocot of NADP-dependent malic enzyme subtype, focussing on cell-specific electron transport capacity. Plants grown under low light (LL) maintained CO₂ assimilation rates similar to high light plants but had an increased chlorophyll and light-harvesting-protein content, predominantly in BS cells. Photosystem II (PSII) protein abundance, oxygen-evolving activity and the PSII/PSI ratio were enhanced in LL BS cells, indicating a higher capacity for linear electron flow. Abundances of PSI, ATP synthase, Cytochrome b₆ f and the chloroplast NAD(P)H dehydrogenase complex, which constitute the BS cyclic electron flow machinery, were also increased in LL plants. A decline in PEP carboxylase activity in mesophyll cells and a consequent shortage of reducing power in BS chloroplasts were associated with a more oxidised plastoquinone pool in LL plants and the formation of PSII - light-harvesting complex II supercomplexes with an increased oxygen evolution rate. Our results suggest that the supramolecular composition of PSII in BS cells is adjusted according to the redox state of the plastoquinone pool. This discovery contributes to the understanding of the acclimation of PSII activity in C₄ plants and will support the development of strategies for crop improvement, including the engineering of C₄ photosynthesis into C₃ plants.

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.

An optimized protocol for isolation of high-quality RNA through laser capture microdissection of leaf material

Laser Capture Microdissection is a powerful tool that allows thin slices of specific cell types to be separated from one another. However, the most commonly used protocol, which involves embedding tissue in paraffin wax, results in severely degraded RNA. Yields from low abundance cell types of leaves are particularly compromised. We reasoned that the relatively high temperature used for sample embedding, and aqueous conditions associated with sample preparation prior to microdissection contribute to RNA degradation. Here, we describe an optimized procedure to limit RNA degradation that is based on the use of low-melting-point wax as well as modifications to sample preparation prior to dissection, and isolation of paradermal, rather than transverse sections. Using this approach, high-quality RNA suitable for down-stream applications such as quantitative reverse transcriptase-polymerase chain reactions or RNA-sequencing is recovered from microdissected bundle sheath strands and mesophyll cells of leaf tissue.

ZmbHLH80 and ZmbHLH90 transcription factors act antagonistically and contribute to regulate PEPC1 cell-specific gene expression in maize

Compartmentation of photosynthetic reactions between mesophyll and bundle sheath cells is a key feature of C₄ photosynthesis and depends on the cell-specific accumulation of major C₄ enzymes, such as phosphoenolpyruvate carboxylase 1. The ZmPEPC1 upstream region, which drives light-inducible and mesophyll-specific gene expression in maize, has been shown to keep the same properties when introduced into rice (C₃ plant), indicating that rice has the transcription factors (TFs) needed to confer C₄ -like gene expression. Using a yeast one-hybrid approach, we identified OsbHLH112, a rice basic Helix-Loop-Helix (bHLH) TF that interacts with the maize ZmPEPC1 upstream region. Moreover, we found that maize OsbHLH112 homologues, ZmbHLH80, and ZmbHLH90, also interact with the ZmPEPC1 upstream region, suggesting that these C₄ regulators were co-opted from C₃ plants. A transactivation assay in maize mesophyll protoplasts revealed that ZmbHLH80 represses, whereas ZmbHLH90 activates, ZmPEPC1 expression. In addition, ZmbHLH80 was shown to impair the ZmPEPC1 promoter activation caused by ZmbHLH90. We showed that ZmbHLH80 and ZmbHLH90 bind to the same cis-element within the ZmPEPC1 upstream region either as homodimers or heterodimers. The formation of homo- and heterodimers with higher oligomeric forms promoted by ZmbHLH80 may explain its negative effect on gene transcription. Gene expression analysis revealed that ZmbHLH80 is preferentially expressed in bundle sheath cells, whereas ZmbHLH90 does not show a clear cell-specific expression pattern. Altogether, our results led us to propose a model in which ZmbHLH80 contributes to mesophyll-specific ZmPEPC1 gene expression by impairing ZmbHLH90-mediated ZmPEPC1 activation in the bundle sheath cells.

Maize leaf PPDK regulatory protein isoform-2 is specific to bundle sheath chloroplasts and paradoxically lacks a Pi-dependent PPDK activation activity

In C4 plants, the pyruvate phosphate dikinase regulatory protein (PDRP) regulates the C4 pathway enzyme pyruvate phosphate dikinase (PPDK) in response to changes in incident light intensity. In maize (Zea mays) leaves, two distinct isoforms of PDRP are expressed, ZmPDRP1 and ZmPDRP2. The properties and C4 function of the ZmPDRP1 isoform are well understood. However, the PDRP2 isoform has only recently been identified and its properties and function(s) in maize leaves are unknown. We therefore initiated an investigation into the maize PDRP2 isoform by performing a side by side comparison of its enzyme properties and cell-specific distribution with PDRP1. In terms of enzyme functionality, PDRP2 was found to possess the same protein kinase-specific activity as PDRP1. However, the PDRP2 isoform was found to lack the phosphotransferase activity of the bifunctional PDRP1 isoform except when PDRP2 in the assays is elevated 5- to 10-fold. A primarily immuno-based approach was used to show that PDRP1 is strictly expressed in mesophyll cells and PDRP2 is strictly expressed in bundle sheath strand cells (BSCs). Additionally, using in situ immunolocalization, we establish a regulatory target for PDRP2 by showing a significant presence of C4 PPDK in BSC chloroplasts. However, a metabolic role for PPDK in this compartment is obscure, assuming PPDK accumulating in this compartment would be irreversibly inactivated each dark cycle by a monofunctional PDRP2.

Freeze-quenched maize mesophyll and bundle sheath separation uncovers bias in previous tissue-specific RNA-Seq data

The high efficiency of C4 photosynthesis relies on spatial division of labor, classically with initial carbon fixation in the mesophyll and carbon reduction in the bundle sheath. By employing grinding and serial filtration over liquid nitrogen, we enriched C4 tissues along a developing leaf gradient. This method treats both C4 tissues in an integrity-preserving and consistent manner, while allowing complementary measurements of metabolite abundance and enzyme activity, thus providing a comprehensive data set. Meta-analysis of this and the previous studies highlights the strengths and weaknesses of different C4 tissue separation techniques. While the method reported here achieves the least enrichment, it is the only one that shows neither strong 3' (degradation) bias, nor different severity of 3' bias between samples. The meta-analysis highlighted previously unappreciated observations, such as an accumulation of evidence that aspartate aminotransferase is more mesophyll specific than expected from the current NADP-ME C4 cycle model, and a shift in enrichment of protein synthesis genes from bundle sheath to mesophyll during development. The full comparative dataset is available for download, and a web visualization tool (available at http://www.plant-biochemistry.hhu.de/resources.html) facilitates comparison of the the Z. mays bundle sheath and mesophyll studies, their consistencies and their conflicts.

Independent and Parallel Evolution of New Genes by Gene Duplication in Two Origins of C4 Photosynthesis Provides New Insight into the Mechanism of Phloem Loading in C4 Species

C4 photosynthesis is considered one of the most remarkable examples of evolutionary convergence in eukaryotes. However, it is unknown whether the evolution of C4 photosynthesis required the evolution of new genes. Genome-wide gene-tree species-tree reconciliation of seven monocot species that span two origins of C4 photosynthesis revealed that there was significant parallelism in the duplication and retention of genes coincident with the evolution of C4 photosynthesis in these lineages. Specifically, 21 orthologous genes were duplicated and retained independently in parallel at both C4 origins. Analysis of this gene cohort revealed that the set of parallel duplicated and retained genes is enriched for genes that are preferentially expressed in bundle sheath cells, the cell type in which photosynthesis was activated during C4 evolution. Furthermore, functional analysis of the cohort of parallel duplicated genes identified SWEET-13 as a potential key transporter in the evolution of C4 photosynthesis in grasses, and provides new insight into the mechanism of phloem loading in these C4 species.

KEY WORDS

C4 photosynthesis, gene duplication, gene families, parallel evolution.

From proto-Kranz to C4 Kranz: building the bridge to C4 photosynthesis

In this review, we examine how the specialized "Kranz" anatomy of C4 photosynthesis evolved from C3 ancestors. Kranz anatomy refers to the wreath-like structural traits that compartmentalize the biochemistry of C4 photosynthesis and enables the concentration of CO2 around Rubisco. A simplified version of Kranz anatomy is also present in the species that utilize C2 photosynthesis, where a photorespiratory glycine shuttle concentrates CO2 into an inner bundle-sheath-like compartment surrounding the vascular tissue. C2 Kranz is considered to be an intermediate stage in the evolutionary development of C4 Kranz, based on the intermediate branching position of C2 species in 14 evolutionary lineages of C4 photosynthesis. In the best-supported model of C4 evolution, Kranz anatomy in C2 species evolved from C3 ancestors with enlarged bundle sheath cells and high vein density. Four independent lineages have been identified where C3 sister species of C2 plants exhibit an increase in organelle numbers in the bundle sheath and enlarged bundle sheath cells. Notably, in all of these species, there is a pronounced shift of mitochondria to the inner bundle sheath wall, forming an incipient version of the C2 type of Kranz anatomy. This incipient version of C2 Kranz anatomy is termed proto-Kranz, and is proposed to scavenge photorespiratory CO2. By doing so, it may provide fitness benefits in hot environments, and thus represent a critical first stage of the evolution of both the C2 and C4 forms of Kranz anatomy.

Evolutionary convergence of cell-specific gene expression in independent lineages of C4 grasses

Leaves of almost all C4 lineages separate the reactions of photosynthesis into the mesophyll (M) and bundle sheath (BS). The extent to which messenger RNA profiles of M and BS cells from independent C4 lineages resemble each other is not known. To address this, we conducted deep sequencing of RNA isolated from the M and BS of Setaria viridis and compared these data with publicly available information from maize (Zea mays). This revealed a high correlation (r=0.89) between the relative abundance of transcripts encoding proteins of the core C4 pathway in M and BS cells in these species, indicating significant convergence in transcript accumulation in these evolutionarily independent C4 lineages. We also found that the vast majority of genes encoding proteins of the C4 cycle in S. viridis are syntenic to homologs used by maize. In both lineages, 122 and 212 homologous transcription factors were preferentially expressed in the M and BS, respectively. Sixteen shared regulators of chloroplast biogenesis were identified, 14 of which were syntenic homologs in maize and S. viridis. In sorghum (Sorghum bicolor), a third C4 grass, we found that 82% of these trans-factors were also differentially expressed in either M or BS cells. Taken together, these data provide, to our knowledge, the first quantification of convergence in transcript abundance in the M and BS cells from independent lineages of C4 grasses. Furthermore, the repeated recruitment of syntenic homologs from large gene families strongly implies that parallel evolution of both structural genes and trans-factors underpins the polyphyletic evolution of this highly complex trait in the monocotyledons.

Transcript residency on ribosomes reveals a key role for the Arabidopsis thaliana bundle sheath in sulfur and glucosinolate metabolism

Leaves of angiosperms are made up of multiple distinct cell types. While the function of mesophyll cells, guard cells, phloem companion cells and sieve elements are clearly described, this is not the case for the bundle sheath (BS). To provide insight into the role of the BS in the C3 species Arabidopsis thaliana, we labelled ribosomes in this cell type with a FLAG tag. We then used immunocapture to isolate these ribosomes, followed by sequencing of resident mRNAs. This showed that 5% of genes showed specific splice forms in the BS, and that 15% of genes were preferentially expressed in these cells. The BS translatome strongly implies that the BS plays specific roles in sulfur transport and metabolism, glucosinolate biosynthesis and trehalose metabolism. Much of the C4 cycle is differentially expressed between the C3 BS and the rest of the leaf. Furthermore, the global patterns of transcript residency on BS ribosomes overlap to a greater extent with cells of the root pericycle than any other cell type. This analysis provides the first insight into the molecular function of this cell type in C3 species, and also identifies characteristics of BS cells that are probably ancestral to both C3 and C4 plants.

Getting the most out of natural variation in C4 photosynthesis

C4 photosynthesis is a complex trait that has a high degree of natural variation, involving anatomical and biochemical changes relative to the ancestral C3 state. It has evolved at least 66 times across a variety of lineages and the evolutionary route from C3 to C4 is likely conserved but not necessarily genetically identical. As such, a variety of C4 species are needed to identify what is fundamental to the C4 evolutionary process in a global context. In order to identify the genetic components of C4 form and function, a number of species are used as genetic models. These include Zea mays (maize), Sorghum bicolor (sorghum), Setaria viridis (Setaria), Flaveria bidentis, and Cleome gynandra. Each of these species has different benefits and challenges associated with its use as a model organism. Here, we propose that RNA profiling of a large sampling of C4, C3-C4, and C3 species, from as many lineages as possible, will allow identification of candidate genes necessary and sufficient to confer C4 anatomy and/or biochemistry. Furthermore, C4 model species will play a critical role in the functional characterization of these candidate genes and identification of their regulatory elements, by providing a platform for transformation and through the use of gene expression profiles in mesophyll and bundle sheath cells and along the leaf developmental gradient. Efforts should be made to sequence the genomes of F. bidentis and C. gynandra and to develop congeneric C3 species as genetic models for comparative studies. In combination, such resources would facilitate discovery of common and unique C4 regulatory mechanisms across genera.

A novel RNA binding protein affects rbcL gene expression and is specific to bundle sheath chloroplasts in C4 plants

BACKGROUND

Plants that utilize the highly efficient C4 pathway of photosynthesis typically possess kranz-type leaf anatomy that consists of two morphologically and functionally distinct photosynthetic cell types, the bundle sheath (BS) and mesophyll (M) cells. These two cell types differentially express many genes that are required for C4 capability and function. In mature C4 leaves, the plastidic rbcL gene, encoding the large subunit of the primary CO2 fixation enzyme Rubisco, is expressed specifically within BS cells. Numerous studies have demonstrated that BS-specific rbcL gene expression is regulated predominantly at post-transcriptional levels, through the control of translation and mRNA stability. The identification of regulatory factors associated with C4 patterns of rbcL gene expression has been an elusive goal for many years.

RESULTS

RLSB, encoded by the nuclear RLSB gene, is an S1-domain RNA binding protein purified from C4 chloroplasts based on its specific binding to plastid-encoded rbcL mRNA in vitro. Co-localized with LSU to chloroplasts, RLSB is highly conserved across many plant species. Most significantly, RLSB localizes specifically to leaf bundle sheath (BS) cells in C4 plants. Comparative analysis using maize (C4) and Arabidopsis (C3) reveals its tight association with rbcL gene expression in both plants. Reduced RLSB expression (through insertion mutation or RNA silencing, respectively) led to reductions in rbcL mRNA accumulation and LSU production. Additional developmental effects, such as virescent/yellow leaves, were likely associated with decreased photosynthetic function and disruption of associated signaling networks.

CONCLUSIONS

Reductions in RLSB expression, due to insertion mutation or gene silencing, are strictly correlated with reductions in rbcL gene expression in both maize and Arabidopsis. In both plants, accumulation of rbcL mRNA as well as synthesis of LSU protein were affected. These findings suggest that specific accumulation and binding of the RLSB binding protein to rbcL mRNA within BS chloroplasts may be one determinant leading to the characteristic cell type-specific localization of Rubisco in C4 plants. Evolutionary modification of RLSB expression, from a C3 "default" state to BS cell-specificity, could represent one mechanism by which rbcL expression has become restricted to only one cell type in C4 plants.