Developmental dynamics of Kranz cell transcriptional specificity in maize leaf reveals early onset of C4-related processes

The WRKY transcription factor family and senescence in switchgrass

BACKGROUND

Early aerial senescence in switchgrass (Panicum virgatum) can significantly limit biomass yields. WRKY transcription factors that can regulate senescence could be used to reprogram senescence and enhance biomass yields.

METHODS

All potential WRKY genes present in the version 1.0 of the switchgrass genome were identified and curated using manual and bioinformatic methods. Expression profiles of WRKY genes in switchgrass flag leaf RNA-Seq datasets were analyzed using clustering and network analyses tools to identify both WRKY and WRKY-associated gene co-expression networks during leaf development and senescence onset.

RESULTS

We identified 240 switchgrass WRKY genes including members of the RW5 and RW6 families of resistance proteins. Weighted gene co-expression network analysis of the flag leaf transcriptomes across development readily separated clusters of co-expressed genes into thirteen modules. A visualization highlighted separation of modules associated with the early and senescence-onset phases of flag leaf growth. The senescence-associated module contained 3000 genes including 23 WRKYs. Putative promoter regions of senescence-associated WRKY genes contained several cis-element-like sequences suggestive of responsiveness to both senescence and stress signaling pathways. A phylogenetic comparison of senescence-associated WRKY genes from switchgrass flag leaf with senescence-associated WRKY genes from other plants revealed notable hotspots in Group I, IIb, and IIe of the phylogenetic tree.

CONCLUSIONS

We have identified and named 240 WRKY genes in the switchgrass genome. Twenty three of these genes show elevated mRNA levels during the onset of flag leaf senescence. Eleven of the WRKY genes were found in hotspots of related senescence-associated genes from multiple species and thus represent promising targets for future switchgrass genetic improvement. Overall, individual WRKY gene expression profiles could be readily linked to developmental stages of flag leaves.

Identification of Photosynthesis-Associated C4 Candidate Genes through Comparative Leaf Gradient Transcriptome in Multiple Lineages of C3 and C4 Species

Leaves of C4 crops usually have higher radiation, water and nitrogen use efficiencies compared to the C3 species. Engineering C4 traits into C3 crops has been proposed as one of the most promising ways to repeal the biomass yield ceiling. To better understand the function of C4 photosynthesis, and to identify candidate genes that are associated with the C4 pathways, a comparative transcription network analysis was conducted on leaf developmental gradients of three C4 species including maize, green foxtail and sorghum and one C3 species, rice. By combining the methods of gene co-expression and differentially co-expression networks, we identified a total of 128 C4 specific genes. Besides the classic C4 shuttle genes, a new set of genes associated with light reaction, starch and sucrose metabolism, metabolites transportation, as well as transcription regulation, were identified as involved in C4 photosynthesis. These findings will provide important insights into the differential gene regulation between C3 and C4 species, and a good genetic resource for establishing C4 pathways in C3 crops.

A Semi-parametric Bayesian Approach for Differential Expression Analysis of RNA-seq Data

RNA-sequencing (RNA-seq) technologies have revolutionized the way agricultural biologists study gene expression as well as generated a tremendous amount of data waiting for analysis. Detecting differentially expressed genes is one of the fundamental steps in RNA-seq data analysis. In this paper, we model the count data from RNA-seq experiments with a Poisson-Gamma hierarchical model, or equivalently, a negative binomial (NB) model. We derive a semi-parametric Bayesian approach with a Dirichlet process as the prior model for the distribution of fold changes between the two treatment means. An inference strategy using Gibbs algorithm is developed for differential expression analysis. The results of several simulation studies show that our proposed method outperforms other methods including the popularly applied edgeR and DESeq methods. We also discuss an application of our method to a dataset that compares gene expression between bundle sheath and mesophyll cells in maize leaves.

Photosynthetic Genes and Genes Associated with the C4 Trait in Maize Are Characterized by a Unique Class of Highly Regulated Histone Acetylation Peaks on Upstream Promoters

Histone modifications contribute to gene regulation in eukaryotes. We analyzed genome-wide histone H3 Lysine (Lys) 4 trimethylation and histone H3 Lys 9 acetylation (two modifications typically associated with active genes) in meristematic cells at the base and expanded cells in the blade of the maize (Zea mays) leaf. These data were compared with transcript levels of associated genes. For individual genes, regulations (fold changes) of histone modifications and transcript levels were much better correlated than absolute intensities. When focusing on regulated histone modification sites, we identified highly regulated secondary H3 Lys 9 acetylation peaks on upstream promoters (regulated secondary upstream peaks [R-SUPs]) on 10% of all genes. R-SUPs were more often found on genes that were up-regulated toward the blade than on down-regulated genes and specifically, photosynthetic genes. Among those genes, we identified six genes encoding enzymes of the C4 cycle and a significant enrichment of genes associated with the C4 trait derived from transcriptomic studies. On the DNA level, R-SUPs are frequently associated with ethylene-responsive elements. Based on these data, we suggest coevolution of epigenetic promoter elements during the establishment of C4 photosynthesis.

Analysis of transcriptional regulation and tissue-specific expression of Avicennia marina Plasma Membrane Protein 3 suggests it contributes to Na(+) transport and homoeostasis in A. marina

Plasma membrane proteins (PMP3) play a role in cation homoeostasis. The 5' flanking sequence of stress inducible, Avicennia marina PMP3 (AmPMP3prom) was transcriptionally fused to (a) GUS or (b) GFP-AmPMP3 and analyzed in transgenic tobacco. Tissue-histochemical GUS and GFP:AmPMP3 localization are co-incident under basal and stress conditions. AmPMP3prom directed GUS activity is highest in roots. Basal transcription is conferred by a 388bp segment upstream of the translation start site. A 463bp distal enhancer in the AmPMP3prom confers enhanced expression under salinity in all tissues and also responds to increases in salinity. The effect of a central, stem-specific negative regulatory region is suppressed by the distal enhancer. The A. marina rhizosphere encounters dynamic changes in salinity at the inter-tidal interface. The complex, tissue-specific transcriptional responsiveness of AmPMP3 to salinity appears to have evolved in response to these changes. Under salinity, guard cell and phloem-specific expression of GFP:AmPMP3 is highly enhanced. Mesophyll, trichomes, bundle sheath, parenchymatous cortex and xylem parenchyma also show GFP:AmPMP3 expression. Cis-elements conferring stress, root and vascular-specific expression are enriched in the AmPMP3 promoter. Pronounced vascular-specific AmPMP3 expression suggests a role in salinity induced Na(+) transport, storage, and secretion in A. marina.

Insights into C4 metabolism from comparative deep sequencing

C4 photosynthesis suppresses the oxygenation activity of Ribulose Bisphosphate Carboxylase Oxygenase and so limits photorespiration. Although highly complex, it is estimated to have evolved in 66 plant lineages, with the vast majority lacking sequenced genomes. Transcriptomics has recently initiated assessments of the degree to which transcript abundance differs between C3 and C4 leaves, identified novel components of C4 metabolism, and also led to mathematical models explaining the repeated evolution of this complex phenotype. Evidence is accumulating that this complex and convergent phenotype is partly underpinned by parallel evolution of structural genes, but also regulatory elements in both cis and trans. Furthermore, it appears that initial events associated with acquisition of C4 traits likely represent evolutionary exaptations related to non-photosynthetic processes.

MET1 is a thylakoid-associated TPR protein involved in photosystem II supercomplex formation and repair in Arabidopsis

Photosystem II (PSII) requires constant disassembly and reassembly to accommodate replacement of the D1 protein. Here, we characterize Arabidopsis thaliana MET1, a PSII assembly factor with PDZ and TPR domains. The maize (Zea mays) MET1 homolog is enriched in mesophyll chloroplasts compared with bundle sheath chloroplasts, and MET1 mRNA and protein levels increase during leaf development concomitant with the thylakoid machinery. MET1 is conserved in C3 and C4 plants and green algae but is not found in prokaryotes. Arabidopsis MET1 is a peripheral thylakoid protein enriched in stroma lamellae and is also present in grana. Split-ubiquitin assays and coimmunoprecipitations showed interaction of MET1 with stromal loops of PSII core components CP43 and CP47. From native gels, we inferred that MET1 associates with PSII subcomplexes formed during the PSII repair cycle. When grown under fluctuating light intensities, the Arabidopsis MET1 null mutant (met1) showed conditional reduced growth, near complete blockage in PSII supercomplex formation, and concomitant increase of unassembled CP43. Growth of met1 in high light resulted in loss of PSII supercomplexes and accelerated D1 degradation. We propose that MET1 functions as a CP43/CP47 chaperone on the stromal side of the membrane during PSII assembly and repair. This function is consistent with the observed differential MET1 accumulation across dimorphic maize chloroplasts.

The evolutionary ecology of C4 plants

C4 photosynthesis is a physiological syndrome resulting from multiple anatomical and biochemical components, which function together to increase the CO2 concentration around Rubisco and reduce photorespiration. It evolved independently multiple times and C4 plants now dominate many biomes, especially in the tropics and subtropics. The C4 syndrome comes in many flavours, with numerous phenotypic realizations of C4 physiology and diverse ecological strategies. In this work, we analyse the events that happened in a C3 context and enabled C4 physiology in the descendants, those that generated the C4 physiology, and those that happened in a C4 background and opened novel ecological niches. Throughout the manuscript, we evaluate the biochemical and physiological evidence in a phylogenetic context, which demonstrates the importance of contingency in evolutionary trajectories and shows how these constrained the realized phenotype. We then discuss the physiological innovations that allowed C4 plants to escape these constraints for two important dimensions of the ecological niche--growth rates and distribution along climatic gradients. This review shows that a comprehensive understanding of C4 plant ecology can be achieved by accounting for evolutionary processes spread over millions of years, including the ancestral condition, functional convergence via independent evolutionary trajectories, and physiological diversification.

Maize and millet transcription factors annotated using comparative genomic and transcriptomic data

Background

Transcription factors (TFs) contain DNA-binding domains (DBDs) and regulate gene expression by binding to specific DNA sequences. In addition, there are proteins, called transcription coregulators (TCs), which lack DBDs but can alter gene expression through interaction with TFs or RNA Polymerase II. Therefore, it is interesting to identify and classify the TFs and TCs in a genome. In this study, maize (Zea mays) and foxtail millet (Setaria italica), two important species for the study of C4 photosynthesis and kranz anatomy, were selected.

Result

We conducted a comprehensive genome-wide annotation of TFs and TCs in maize B73 and in two strains of foxtail millet, Zhang gu and Yugu1, and classified them into families. To gain additional support for our predictions, we searched for their homologous genes in Arabidopsis or rice and studied their gene expression level using RNA-seq and microarray data. We identified many new TF and TC families in these two species, and described some evolutionary and functional aspects of the 9 new maize TF families. Moreover, we detected many pseudogenes and transposable elements in current databases. In addition, we examined tissue expression preferences of TF and TC families and identified tissue/condition-specific TFs and TCs in maize and millet. Finally, we identified potential C4-related TF and TC genes in maize and millet.

Conclusions

Our results significantly expand current TF and TC annotations in maize and millet. We provided supporting evidence for our annotation from genomic and gene expression data and identified TF and TC genes with tissue preference in expression. Our study may facilitate the study of regulation of gene expression, tissue morphogenesis, and C4 photosynthesis in maize and millet. The data we generated in this study are available at http://sites.google.com/site/jjlmmtf.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-818) contains supplementary material, which is available to authorized users.