Shared and genetically distinct Zea mays transcriptome responses to ongoing and past low temperature exposure

Cellular dynamics in the maize leaf growth zone during recovery from chilling depends on the leaf developmental stage

KEY MESSAGE

A novel non-steady-state kinematic analysis shows differences in cell division and expansion determining a better recovery from a 3-day cold spell in emerged compared to non-emerged maize leaves. Zea mays is highly sensitive to chilling which frequently occurs during its seedling stage. Although the direct effect of chilling is well studied, the mechanisms determining the subsequent recovery are still unknown. Our goal is to determine the cellular basis of the leaf growth response to chilling and during recovery of leaves exposed before or after their emergence. We first studied the effect of a 3-day cold spell on leaf growth at the plant level. Then, we performed a kinematic analysis to analyse the dynamics of cell division and elongation during recovery of the 4th leaf after exposure to cold before or after emergence. Our results demonstrated cold more strongly reduced the final length of non-emerged than emerged leaves (- 13 vs. - 18%). This was not related to growth differences during cold, but a faster and more complete recovery of the growth of emerged leaves. This difference was due to a higher cell division rate on the 1st and a higher cell elongation rate on the 2nd day of recovery, respectively. The dynamics of cell division and expansion during recovery determines developmental stage-specific differences in cold tolerance of maize leaves.

Identifying yield-related genes in maize based on ear trait plasticity

BACKGROUND

Phenotypic plasticity is defined as the phenotypic variation of a trait when an organism is exposed to different environments, and it is closely related to genotype. Exploring the genetic basis behind the phenotypic plasticity of ear traits in maize is critical to achieve climate-stable yields, particularly given the unpredictable effects of climate change. Performing genetic field studies in maize requires development of a fast, reliable, and automated system for phenotyping large numbers of samples.

RESULTS

Here, we develop MAIZTRO as an automated maize ear phenotyping platform for high-throughput measurements in the field. Using this platform, we analyze 15 common ear phenotypes and their phenotypic plasticity variation in 3819 transgenic maize inbred lines targeting 717 genes, along with the wild type lines of the same genetic background, in multiple field environments in two consecutive years. Kernel number is chosen as the primary target phenotype because it is a key trait for improving the grain yield and ensuring yield stability. We analyze the phenotypic plasticity of the transgenic lines in different environments and identify 34 candidate genes that may regulate the phenotypic plasticity of kernel number.

CONCLUSIONS

Our results suggest that as an integrated and efficient phenotyping platform for measuring maize ear traits, MAIZTRO can help to explore new traits that are important for improving and stabilizing the yield. This study indicates that genes and alleles related with ear trait plasticity can be identified using transgenic maize inbred populations.

Transcriptomic Analysis Reveals the Correlation between End-of-Day Far Red Light and Chilling Stress in Setaria viridis

Low temperature and end-of-day far-red (EOD-FR) light signaling are two key factors limiting plant production and geographical location worldwide. However, the transcriptional dynamics of EOD-FR light conditions during chilling stress remain poorly understood. Here, we performed a comparative RNA-Seq-based approach to identify differentially expressed genes (DEGs) related to EOD-FR and chilling stress in Setaria viridis. A total of 7911, 324, and 13431 DEGs that responded to low temperature, EOD-FR and these two stresses were detected, respectively. Further DEGs analysis revealed that EOD-FR may enhance cold tolerance in plants by regulating the expression of genes related to cold tolerance. The result of weighted gene coexpression network analysis (WGCNA) using 13431 nonredundant DEGs exhibited 15 different gene network modules. Interestingly, a CO-like transcription factor named BBX2 was highly expressed under EOD-FR or chilling conditions. Furthermore, we could detect more expression levels when EOD-FR and chilling stress co-existed. Our dataset provides a valuable resource for the regulatory network involved in EOD-FR signaling and chilling tolerance in C₄ plants.

Understanding and Comprehensive Evaluation of Cold Resistance in the Seedlings of Multiple Maize Genotypes

Maize is a cold-sensitive crop, and it exhibits severe retardation of growth and development when exposed to cold snaps during and right after seedling emergence. Although different agronomic, physiological, and molecular approaches have been tried to overcome the problems related to cold stress in recent years, the mechanisms causing cold resistance in maize are still unclear. Screening and breeding of varieties for cold resistance may be a sustainable option to boost maize production under low-temperature environments. Herein, seedlings of 39 different maize genotypes were treated under both 10 °C low temperature and 22 °C normal temperature conditions for 7 days, to assess the changes in seven growth parameters, two membrane characteristics, two reactive oxygen species (ROS) levels, and four antioxidant enzymes activities. The changes in ten photosynthetic performances, one osmotic substance accumulation, and three polyamines (PAs) metabolisms were also measured. Results indicated that significant differences among genotypes, temperature treatments, and their interactions were found in 29 studied traits, and cold–stressed seedlings were capable to enhance their cold resistance by maintaining high levels of membrane stability index (66.07%); antioxidant enzymes activities including the activity of superoxide dismutase (2.44 Unit g⁻¹ protein), peroxidase (1.65 Unit g⁻¹ protein), catalase (0.65 μM min⁻¹ g⁻¹ protein), and ascorbate peroxidase (5.45 μM min⁻¹ g⁻¹ protein); chlorophyll (Chl) content, i.e., Chl a (0.36 mg g⁻¹ FW) and Chl b (0.40 mg g⁻¹ FW); photosynthetic capacity such as net photosynthetic rate (5.52 μM m⁻² s⁻¹) and ribulose 1,5–biphosphate carboxylase activity (6.57 M m⁻² s⁻¹); PAs concentration, mainly putrescine (274.89 nM g⁻¹ FW), spermidine (52.69 nM g⁻¹ FW), and spermine (45.81 nM g⁻¹ FW), particularly under extended cold stress. Importantly, 16 traits can be good indicators for screening of cold–resistant genotypes of maize. Gene expression analysis showed that GRMZM2G059991, GRMZM2G089982, GRMZM2G088212, GRMZM2G396553, GRMZM2G120578, and GRMZM2G396856 involved in antioxidant enzymes activity and PAs metabolism, and these genes may be used for genetic modification to improve maize cold resistance. Moreover, seven strong cold–resistant genotypes were identified, and they can be used as parents in maize breeding programs to develop new varieties.

Can we improve the chilling tolerance of maize photosynthesis through breeding?

Chilling tolerance is necessary for crops to thrive in temperate regions where cold snaps and lower baseline temperatures place limits on life processes; this is particularly true for crops of tropical origin such as maize. Photosynthesis is often adversely affected by chilling stress, yet the maintenance of photosynthesis is essential for healthy growth and development, and most crucially for yield. In this review, we describe the physiological basis for enhancing chilling tolerance of photosynthesis in maize by examining nine key responses to chilling stress. We synthesize current knowledge of genetic variation for photosynthetic chilling tolerance in maize with respect to each of these traits and summarize the extent to which genetic mapping and candidate genes have been used to understand the genomic regions underpinning chilling tolerance. Finally, we provide perspectives on the future of breeding for photosynthetic chilling tolerance in maize. We advocate for holistic and high-throughput approaches to screen for chilling tolerance of photosynthesis in research and breeding programmes in order to develop resilient crops for the future.

Prediction of conserved and variable heat and cold stress response in maize using cis-regulatory information

Changes in gene expression are important for responses to abiotic stress. Transcriptome profiling of heat- or cold-stressed maize genotypes identifies many changes in transcript abundance. We used comparisons of expression responses in multiple genotypes to identify alleles with variable responses to heat or cold stress and to distinguish examples of cis- or trans-regulatory variation for stress-responsive expression changes. We used motifs enriched near the transcription start sites (TSSs) for thermal stress-responsive genes to develop predictive models of gene expression responses. Prediction accuracies can be improved by focusing only on motifs within unmethylated regions near the TSS and vary for genes with different dynamic responses to stress. Models trained on expression responses in a single genotype and promoter sequences provided lower performance when applied to other genotypes but this could be improved by using models trained on data from all three genotypes tested. The analysis of genes with cis-regulatory variation provides evidence for structural variants that result in presence/absence of transcription factor binding sites in creating variable responses. This study provides insights into cis-regulatory motifs for heat- and cold-responsive gene expression and defines a framework for developing models to predict expression responses across multiple genotypes.

Local adaptation contributes to gene expression divergence in maize

Gene expression links genotypes to phenotypes, so identifying genes whose expression is shaped by selection will be important for understanding the traits and processes underlying local adaptation. However, detecting local adaptation for gene expression will require distinguishing between divergence due to selection and divergence due to genetic drift. Here, we adapt a QST−FST framework to detect local adaptation for transcriptome-wide gene expression levels in a population of diverse maize genotypes. We compare the number and types of selected genes across a wide range of maize populations and tissues, as well as selection on cold-response genes, drought-response genes, and coexpression clusters. We identify a number of genes whose expression levels are consistent with local adaptation and show that genes involved in stress response show enrichment for selection. Due to its history of intense selective breeding and domestication, maize evolution has long been of interest to researchers, and our study provides insight into the genes and processes important for in local adaptation of maize.

Zea mays RNA-seq estimated transcript abundances are strongly affected by read mapping bias

Background

Genetic variation for gene expression is a source of phenotypic variation for natural and agricultural species. The common approach to map and to quantify gene expression from genetically distinct individuals is to assign their RNA-seq reads to a single reference genome. However, RNA-seq reads from alleles dissimilar to this reference genome may fail to map correctly, causing transcript levels to be underestimated. Presently, the extent of this mapping problem is not clear, particularly in highly diverse species. We investigated if mapping bias occurred and if chromosomal features associated with mapping bias. Zea mays presents a model species to assess these questions, given it has genotypically distinct and well-studied genetic lines.

Results

In Zea mays, the inbred B73 genome is the standard reference genome and template for RNA-seq read assignments. In the absence of mapping bias, B73 and a second inbred line, Mo17, would each have an approximately equal number of regulatory alleles that increase gene expression. Remarkably, Mo17 had 2–4 times fewer such positively acting alleles than did B73 when RNA-seq reads were aligned to the B73 reference genome. Reciprocally, over one-half of the B73 alleles that increased gene expression were not detected when reads were aligned to the Mo17 genome template. Genes at dissimilar chromosomal ends were strongly affected by mapping bias, and genes at more similar pericentromeric regions were less affected. Biased transcript estimates were higher in untranslated regions and lower in splice junctions. Bias occurred across software and alignment parameters.

Conclusions

Mapping bias very strongly affects gene transcript abundance estimates in maize, and bias varies across chromosomal features. Individual genome or transcriptome templates are likely necessary for accurate transcript estimation across genetically variable individuals in maize and other species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07577-3.

A temporal hierarchy underpins the transcription factor-DNA interactome of the maize UPR

Adverse environmental conditions reduce crop productivity and often increase the load of unfolded or misfolded proteins in the endoplasmic reticulum (ER). This potentially lethal condition, known as ER stress, is buffered by the unfolded protein response (UPR), a set of signaling pathways designed to either recover ER functionality or ignite programmed cell death. Despite the biological significance of the UPR to the life of the organism, the regulatory transcriptional landscape underpinning ER stress management is largely unmapped, especially in crops. To fill this significant knowledge gap, we performed a large-scale systems-level analysis of the protein-DNA interaction (PDI) network in maize (Zea mays). Using 23 promoter fragments of six UPR marker genes in a high-throughput enhanced yeast one-hybrid assay, we identified a highly interconnected network of 262 transcription factors (TFs) associated with significant biological traits and 831 PDIs underlying the UPR. We established a temporal hierarchy of TF binding to gene promoters within the same family as well as across different families of TFs. Cistrome analysis revealed the dynamic activities of a variety of cis-regulatory elements (CREs) in ER stress-responsive gene promoters. By integrating the cistrome results into a TF network analysis, we mapped a subnetwork of TFs associated with a CRE that may contribute to UPR management. Finally, we validated the role of a predicted network hub gene using the Arabidopsis system. The PDIs, TF networks, and CREs identified in our work are foundational resources for understanding transcription-regulatory mechanisms in the stress responses and crop improvement.

Changes in photosynthesis and respiratory metabolism of maize seedlings growing under low temperature stress may be regulated by arbuscular mycorrhizal fungi

Arbuscular mycorrhizal fungi as an important soil microbe have been demonstrated to mitigate the harmful effects of stress on plants. However, little is known about the molecular mechanisms underlying the AM symbiotic response to low temperature. Here, differentially expressed genes (DEGs) in the maize seedlings were identified after inoculating AMF under low temperature conditions. A total of 10,400 DEGs were obtained among four treatments, including non-inoculated AMF under ambient temperature (NMA), inoculated AMF under ambient temperature (MA), non-inoculated with low temperature stress (NML), and inoculated with low temperature stress (ML). The relative expression of 858 genes increased and that of 497 genes decreased in AM plants under low temperature stress. 24 DEGs were identified related to photosynthesis and respiratory metabolism. Among these DEGs, 10 genes were upregulated, and 14 genes were downregulated. The results show that inoculating AMF might decrease the production and transmission of electrons under low temperature, and the cyclic electron flow process in chloroplasts was stimulated to protect plants against low temperature. The fungi also influenced transmission of electrons and production of phosphoric acid in mitochondria in response to low temperature. CO₂ assimilation capacity was affected and the tricarboxylic acid cycle was promoted by the adjustments in the glycolysis, pentose phosphate pathway, gamma-aminobutyric acid shunt pathway, and glyoxylic acid cycle to produce more adenosine triphosphate and raw materials for other metabolic pathways under low temperature. These findings provide new insight into low temperature tolerance induced by AMF, and help identify genes for further investigation and functional analyses.

Elucidating the regulatory roles of microRNAs in maize (Zea mays L.) leaf growth response to chilling stress

MAIN CONCLUSION: miRNAs control leaf size of maize crop during chilling stress tolerance by regulating developmentally important transcriptional factors and sustaining redox homeostasis of cells. Chilling temperature (0-15 °C) is a major constraint for the cultivation of maize (Zea mays) which inhibits the early growth of maize leading to reduction in leaf size. Growth and development take place in meristem, elongation, and mature zones that are linearly located along the leaf base to tip. To prevent shortening of leaf caused by chilling, this study aims to elucidate the regulatory roles of microRNA (miRNA) genes in the controlling process switching between growth and developmental stages. In this respect, hybrid maize ADA313 seedlings were treated to the chilling temperature which caused 26% and 29% reduction in the final leaf length and a decline in cell production of the fourth leaf. The flow cytometry data integrated with the expression analysis of cell cycle genes indicated that the reason for the decline was a failure proceeding from G2/M rather than G1/S. Through an miRNome analysis of 321 known maize miRNAs, 24, 6, and 20 miRNAs were assigned to putative meristem, elongation, and mature zones, respectively according to their chilling response. To gain deeper insight into decreased cell production, in silico, target prediction analysis was performed for meristem specific miRNAs. Among the miRNAs, miR160, miR319, miR395, miR396, miR408, miR528, and miR1432 were selected for confirming the potential of negative regulation with their predicted targets by qRT-PCR. These findings indicated evidence for improvement of growth and yield under chilling stress of the maize.

Maize Response to Low Temperatures at the Gene Expression Level: A Critical Survey of Transcriptomic Studies

Maize is a cold-sensitive plant whose physiological reactions to sub-optimal temperatures are well understood, but their molecular foundations are only beginning to be deciphered. In an attempt to identify key genes involved in these reactions, we surveyed several independent transcriptomic studies addressing the response of juvenile maize to moderate or severe cold. Among the tens of thousands of genes found to change expression upon cold treatment less than 500 were reported in more than one study, indicating an astonishing variability of the expression changes, likely depending on the experimental design and plant material used. Nearly all these "common" genes were specific to either moderate or to severe cold and formed distinct interaction networks, indicating fundamentally different responses. Moreover, down-regulation of gene expression dominated strongly in moderate cold and up-regulation prevailed in severe cold. Very few of these genes have ever been mentioned in the literature as cold-stress-related, indicating that most response pathways remain poorly known at the molecular level. We posit that the genes identified by the present analysis are attractive candidates for further functional studies and their arrangement in complex interaction networks indicates that a re-interpretation of the present state of knowledge on the maize cold-response is justified.