An updated gene atlas for maize reveals organ-specific and stress-induced genes

20 years of the Bio-Analytic Resource for Plant Biology

The Bio-Analytic Resource for Plant Biology ('the BAR', at https://bar.utoronto.ca) is celebrating its 20th year in operation in 2025. The BAR encompasses and provides visualization tools for large 'omics data sets from plants. The BAR covers data from Arabidopsis, tomato, wheat, barley and 29 other plant species (with data for 2 others to be released soon). These data include nucleotide and protein sequence data, gene expression data, protein-protein and protein-DNA interactions, protein structures, subcellular localizations, and polymorphisms. The data are stored in more than 200 relational databases holding 186 GB of data and are presented to the researchers via web apps. These web apps provide data analysis and visualization tools. Some of the most popular tools are eFP ('electronic fluorescent pictograph') Browsers, ePlants and ThaleMine (an Arabidopsis-specific instance of InterMine). The BAR was designated a Global Core Biodata Resource in 2023. Like other GCBRs, the BAR has excellent operational stability, provides access without login requirement, and provides an API for researchers to be able to access BAR data programmatically. We present in this update a new overarching search tool called Gaia that permits easy access to all BAR data, powered by machine learning and artificial intelligence.

The Zea mays PeptideAtlas: A New Maize Community Resource

This study presents the Maize PeptideAtlas resource (www.peptideatlas.org/builds/maize) to help solve questions about the maize proteome. Publicly available raw tandem mass spectrometry (MS/MS) data for maize collected from ProteomeXchange were reanalyzed through a uniform processing and metadata annotation pipeline. These data are from a wide range of genetic backgrounds and many sample types and experimental conditions. The protein search space included different maize genome annotations for the B73 inbred line from MaizeGDB, UniProtKB, NCBI RefSeq, and for the W22 inbred line. 445 million MS/MS spectra were searched, of which 120 million were matched to 0.37 million distinct peptides. Peptides were matched to 66.2% of proteins in the most recent B73 nuclear genome annotation. Furthermore, most conserved plastid- and mitochondrial-encoded proteins (NCBI RefSeq annotations) were identified. Peptides and proteins identified in the other B73 genome annotations will improve maize genome annotation. We also illustrate the high-confidence detection of unique W22 proteins. N-terminal acetylation, phosphorylation, ubiquitination, and three lysine acylations (K-acetyl, K-malonyl, and K-hydroxyisobutyryl) were identified and can be inspected through a PTM viewer in PeptideAtlas. All matched MS/MS-derived peptide data are linked to spectral, technical, and biological metadata. This new PeptideAtlas is integrated in MaizeGDB with a peptide track in JBrowse.

A maize enzyme from the 2-oxoglutarate-dependent oxygenase family with unique kinetic properties, mediates resistance against pathogens and regulates senescence

In plants, salicylic acid (SA) hydroxylation regulates SA homoeostasis, playing an essential role during plant development and response to pathogens. This reaction is catalysed by SA hydroxylase enzymes, which hydroxylate SA producing 2,3-dihydroxybenzoic acid (2,3-DHBA) and/or 2,5-dihydroxybenzoic acid (2,5-DHBA). Several SA hydroxylases have recently been identified and characterised from different plant species, but no such activity has yet been reported in maize. In this work, we describe the identification and characterisation of a new SA hydroxylase in maize plants. This enzyme, with high sequence similarity to previously described SA hydroxylases from Arabidopsis and rice, converts SA into 2,5-DHBA; however, it has different kinetic properties to those of previously characterised enzymes, and it also catalysers the conversion of the flavonoid dihydroquercetin into quercetin in in vitro activity assays, suggesting that the maize enzyme may have different roles in vivo to those previously reported from other species. Despite this, ZmS5H can complement the pathogen resistance and the early senescence phenotypes of Arabidopsis s3h mutant plants. Finally, we characterised a maize mutant in the S5H gene (s5hMu) that has altered growth, senescence and increased resistance against Colletotrichum graminicola infection, showing not only alterations in SA and 2,5-DHBA but also in flavonol levels. Together, the results presented here provide evidence that SA hydroxylases in different plant species have evolved to show differences in catalytic properties that may be important to fine tune SA levels and other phenolic compounds such as flavonols, to regulate different aspects of plant development and pathogen defence.

Prioritizing Maize Metabolic Gene Regulators through Multi-Omic Network Integration

Elucidating gene regulatory networks is a major area of study within plant systems biology. Phenotypic traits are intricately linked to specific gene expression profiles. These expression patterns arise primarily from regulatory connections between sets of transcription factors (TFs) and their target genes. Here, we integrated 46 co-expression networks, 283 protein-DNA interaction (PDI) assays, and 16 million SNPs used to identify expression quantitative trait loci (eQTL) to construct TF-target networks. In total, we analyzed ∼4.6M interactions to generate four distinct types of TF-target networks: co-expression, PDI, trans -eQTL, and cis -eQTL combined with PDIs. To functionally annotate TFs based on their target genes, we implemented three different network integration strategies. We evaluated the effectiveness of each strategy through TF loss-of function mutant inspection and random network analyses. The multi-network integration allowed us to identify transcriptional regulators of several biological processes. Using the topological properties of the fully integrated network, we identified potential functionally redundant TF paralogs. Our findings retrieved functions previously documented for numerous TFs and revealed novel functions that are crucial for informing the design of future experiments. The approach here-described lays the foundation for the integration of multi-omic datasets in maize and other plant systems.

GRAPHICAL ABSTRACT

Sugar sensing in C4 source leaves: a gap that needs to be filled

Plant growth depends on sugar production and export by photosynthesizing source leaves and sugar allocation and import by sink tissues (grains, roots, stems, and young leaves). Photosynthesis and sink demand are tightly coordinated through metabolic (substrate, allosteric) feedback and signalling (sugar, hormones) mechanisms. Sugar signalling integrates sugar production with plant development and environmental cues. In C3 plants (e.g. wheat and rice), it is well documented that sugar accumulation in source leaves, due to source-sink imbalance, negatively feeds back on photosynthesis and plant productivity. However, we have a limited understanding about the molecular mechanisms underlying those feedback regulations, especially in C4 plants (e.g. maize, sorghum, and sugarcane). Recent work with the C4 model plant Setaria viridis suggested that C4 leaves have different sugar sensing thresholds and behaviours relative to C3 counterparts. Addressing this research priority is critical because improving crop yield requires a better understanding of how plants coordinate source activity with sink demand. Here we review the literature, present a model of action for sugar sensing in C4 source leaves, and suggest ways forward.

Identification of candidate gene associated with maize northern leaf blight resistance in a multi-parent population

KEY MESSAGE

QTL mapping combined with genome-wide association studies, revealed a potential candidate gene for  resistance to northern leaf blight in the tropical CATETO-related maize line YML226, providing a basis for marker-assisted selection of maize varieties Northern leaf blight (NLB) is a foliar disease that can cause severe yield losses in maize. Identifying and utilizing NLB-resistant genes is the most effective way to prevent and control this disease. In this study, five important inbred lines of maize were used as parental lines to construct a multi-parent population for the identification of NLB-resistant loci. QTL mapping and GWAS analysis revealed that QTL qtl_YML226_1, which had the largest phenotypic variance explanation (PVE) of 9.28%, and SNP 5-49,193,921 were co-located in the CATETO-related line YML226. This locus was associated with the candidate gene Zm00001d014471, which encodes a pentatricopeptide repeat (PPR) protein. In the coding region of Zm00001d014471, YML226 had more specific SNPs than the other parental lines. qRT-PCR showed that the relative expressions of Zm00001d014471 in inoculated and uninoculated leaves of YML226 were significantly higher, indicating that the expression of the candidate gene was correlated with NLB resistance. The analysis showed that the higher expression level in YML226 might be caused by SNP mutations. This study identified NLB resistance candidate loci and genes in the tropical maize inbred line YML226 derived from the CATETO germplasm, thereby providing a theoretical basis for using modern marker-assisted breeding techniques to select genetic resources resistant to NLB.

Genome-wide identification of ZmMYC2 binding sites and target genes in maize

BACKGROUND

Jasmonate (JA) is the important phytohormone to regulate plant growth and adaption to stress signals. MYC2, an bHLH transcription factor, is the master regulator of JA signaling. Although MYC2 in maize has been identified, its function remains to be clarified.

RESULTS

To understand the function and regulatory mechanism of MYC2 in maize, the joint analysis of DAP-seq and RNA-seq is conducted to identify the binding sites and target genes of ZmMYC2. A total of 3183 genes are detected both in DAP-seq and RNA-seq data, potentially as the directly regulating genes of ZmMYC2. These genes are involved in various biological processes including plant growth and stress response. Besides the classic cis-elements like the G-box and E-box that are bound by MYC2, some new motifs are also revealed to be recognized by ZmMYC2, such as nGCATGCAnn, AAAAAAAA, CACGTGCGTGCG. The binding sites of many ZmMYC2 regulating genes are identified by IGV-sRNA.

CONCLUSIONS

All together, abundant target genes of ZmMYC2 are characterized with their binding sites, providing the basis to construct the regulatory network of ZmMYC2 and better understanding for JA signaling in maize.

Dynamic transcriptome landscape of maize pericarp development

As a maternal tissue, the pericarp supports and protects for other components of seed, such as embryo and endosperm. Despite the importance of maize pericarp in seed, the genome-wide transcriptome pattern throughout maize pericarp development has not been well characterized. Here, we developed RNA-seq transcriptome atlas of B73 maize pericarp development based on 21 samples from 5 days before fertilization (DBP5) to 32 days after fertilization (DAP32). A total of 25 346 genes were detected in programming pericarp development, including 1887 transcription factors (TFs). Together with pericarp morphological changes, the global clustering of gene expression revealed four developmental stages: undeveloped, thickening, expansion and strengthening. Coexpression analysis provided further insights on key regulators in functional transition of four developmental stages. Combined with non-seed, embryo, endosperm, and nucellus transcriptome data, we identified 598 pericarp-specific genes, including 75 TFs, which could elucidate key mechanisms and regulatory networks of pericarp development. Cell wall related genes were identified that reflected their crucial role in the maize pericarp structure building. In addition, key maternal proteases or TFs related with programmed cell death (PCD) were proposed, suggesting PCD in the maize pericarp was mediated by vacuolar processing enzymes (VPE), and jasmonic acid (JA) and ethylene-related pathways. The dynamic transcriptome atlas provides a valuable resource for unraveling the genetic control of maize pericarp development.

Insights into ZmWAKL in maize kernel development: genome-wide investigation and GA-mediated transcription

BACKGROUND

The functional roles of the Wall Associated Kinase (WAK) and Wall Associated Kinase Like (WAKL) families in cellular expansion and developmental processes have been well-established. However, the molecular regulation of these kinases in maize development is limited due to the absence of comprehensive genome-wide studies.

RESULTS

Through an in-depth analysis, we identified 58 maize WAKL genes, and classified them into three distinct phylogenetic clusters. Moreover, structural prediction analysis showed functional conservation among WAKLs across maize. Promoter analysis uncovered the existence of cis-acting elements associated with the transcriptional regulation of ZmWAKL genes by Gibberellic acid (GA). To further elucidate the role of WAKL genes in maize kernels, we focused on three highly expressed genes, viz ZmWAKL38, ZmWAKL42 and ZmWAKL52. Co-expression analyses revealed that their expression patterns exhibited a remarkable correlation with GA-responsive transcription factors (TF) TF5, TF6, and TF8, which displayed preferential expression in kernels. RT-qPCR analysis validated the upregulation of ZmWAKL38, ZmWAKL42, ZmWAKL52, TF5, TF6, and TF8 following GA treatment. Additionally, ZmWAKL52 showed significant increase of transcription in the present of TF8, with ZmWAKL52 localizing in both the plasma membrane and cell wall. TF5 positively regulated ZmWAKL38, while TF6 positively regulated ZmWAKL42.

CONCLUSIONS

Collectively, these findings provide novel insights into the characterization and regulatory mechanisms of specific ZmWAKL genes involved in maize kernel development, offering prospects for their utilization in maize breeding programs.

Biofortification of iron content by regulating a NAC transcription factor in maize

Iron (Fe) deficiency remains widespread among people in developing countries. To help solve this problem, breeders have been attempting to develop maize cultivars with high yields and high Fe concentrations in the kernels. We conducted a genome-wide association study and identified a gene, ZmNAC78 (NAM/ATAF/CUC DOMAIN TRANSCRIPTION FACTOR 78), that regulates Fe concentrations in maize kernels. We cultivated maize varieties with both high yield and high Fe concentrations in their kernels by using a molecular marker developed from a 42-base pair insertion or deletion (indel) in the promoter of ZmNAC78. ZmNAC78 expression is enriched in the basal endosperm transfer layer of kernels, and the ZmNAC78 protein directly regulates messenger RNA abundance of Fe transporters. Our results thus provide an approach to develop maize varieties with Fe-enriched kernels.

Identification of novel loci associated with starch content in maize kernels by a genome-wide association study using an enlarged SNP panel

Starch is a major component of cereals, comprising over 70% of dry weight. It serves as a primary carbon source for humans and animals. In addition, starch is an indispensable industrial raw material. While maize (Zea mays) is a key crop and the primary source of starch, the genetic basis for starch content in maize kernels remains poorly understood. In this study, using an enlarged panel, we conducted a genome-wide association study (GWAS) based on best linear unbiased prediction (BLUP) value for starch content of 261 inbred lines across three environments. Compared with previous study, we identified 14 additional significant quantitative trait loci (QTL), encompassed a total of 42 genes, and indicated that increased marker density contributes to improved statistical power. By integrating gene expression profiling, Gene Ontology (GO) enrichment and haplotype analysis, several potential target genes that may play a role in regulating starch content in maize kernels have been identified. Notably, we found that ZmAPC4, associated with the significant SNP chr4.S_175584318, which encodes a WD40 repeat-like superfamily protein and is highly expressed in maize endosperm, might be a crucial regulator of maize kernel starch synthesis. Out of the 261 inbred lines analyzed, they were categorized into four haplotypes. Remarkably, it was observed that the inbred lines harboring hap4 demonstrated the highest starch content compared to the other haplotypes. Additionally, as a significant achievement, we have developed molecular markers that effectively differentiate maize inbred lines based on their starch content. Overall, our study provides valuable insights into the genetic basis of starch content and the molecular markers can be useful in breeding programs aimed at developing maize varieties with high starch content, thereby improving breeding efficiency.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01437-6.

RNA pseudouridine modification in plants

Pseudouridine is one of the well-known chemical modifications in various RNA species. Current advances to detect pseudouridine show that the pseudouridine landscape is dynamic and affects multiple cellular processes. Although our understanding of this post-transcriptional modification mainly depends on yeast and human models, the recent findings provide strong evidence for the critical role of pseudouridine in plants. Here, we review the current knowledge of pseudouridine in plant RNAs, including its synthesis, degradation, regulatory mechanisms, and functions. Moreover, we propose future areas of research on pseudouridine modification in plants.

VAMP726 from maize and Arabidopsis confers pollen resistance to heat and UV radiation by influencing lignin content of sporopollenin

Sporopollenin in the pollen cell wall protects male gametophytes from stresses. Phenylpropanoid derivatives, including guaiacyl (G) lignin units, are known to be structural components of sporopollenin, but the exact composition of sporopollenin remains to be fully resolved. We analyzed the phenylpropanoid derivatives in sporopollenin from maize and Arabidopsis by thioacidolysis coupled with nuclear magnetic resonance (NMR) and gas chromatography-mass spectrometry (GC-MS). The NMR and GC-MS results confirmed the presence of p-hydroxyphenyl (H), G, and syringyl (S) lignin units in sporopollenin from maize and Arabidopsis. Strikingly, H units account for the majority of lignin monomers in sporopollenin from these species. We next performed a genome-wide association study to explore the genetic basis of maize sporopollenin composition and identified a vesicle-associated membrane protein (ZmVAMP726) that is strongly associated with lignin monomer composition of maize sporopollenin. Genetic manipulation of VAMP726 affected not only lignin monomer composition in sporopollenin but also pollen resistance to heat and UV radiation in maize and Arabidopsis, indicating that VAMP726 is functionally conserved in monocot and dicot plants. Our work provides new insight into the lignin monomers that serve as structural components of sporopollenin and characterizes VAMP726, which affects sporopollenin composition and stress resistance in pollen.

Maize Feature Store: A centralized resource to manage and analyze curated maize multi-omics features for machine learning applications

The big-data analysis of complex data associated with maize genomes accelerates genetic research and improves agronomic traits. As a result, efforts have increased to integrate diverse datasets and extract meaning from these measurements. Machine learning models are a powerful tool for gaining knowledge from large and complex datasets. However, these models must be trained on high-quality features to succeed. Currently, there are no solutions to host maize multi-omics datasets with end-to-end solutions for evaluating and linking features to target gene annotations. Our work presents the Maize Feature Store (MFS), a versatile application that combines features built on complex data to facilitate exploration, modeling and analysis. Feature stores allow researchers to rapidly deploy machine learning applications by managing and providing access to frequently used features. We populated the MFS for the maize reference genome with over 14 000 gene-based features based on published genomic, transcriptomic, epigenomic, variomic and proteomics datasets. Using the MFS, we created an accurate pan-genome classification model with an AUC-ROC score of 0.87. The MFS is publicly available through the maize genetics and genomics database. Database URL  https://mfs.maizegdb.org/.

Characterization of the Isocitrate Dehydrogenase Gene Family and Their Response to Drought Stress in Maize

Isocitrate dehydrogenase (IDH) is a key rate-limiting enzyme in the tricarboxylic acid cycle and acts in glutamine synthesis. IDH also participates in plant growth and development and in response to abiotic stresses. We identified 11 maize IDH genes (ZmIDH) and classified these genes into ZmNAD-IDH and ZmNADP-IDH groups based on their different coenzymes (NAD+ or NADP+). The ZmNAD-IDH group was further divided into two subgroups according to their catalytic and non-catalytic subunits, as in Arabidopsis. The ZmIDHs significantly differed in physicochemical properties, gene structure, conserved motifs, and protein tertiary structure. Promoter prediction analysis revealed that the promoters of these ZmIDHs contain cis-acting elements associated with light response, abscisic acid, phytohormones, and abiotic stresses. ZmIDH is predicted to interact with proteins involved in development and stress resistance. Expression analysis of public data revealed that most ZmIDHs are specifically expressed in anthers. Different types of ZmIDHs responded to abiotic stresses with different expression patterns, but all exhibited responses to abiotic stresses to some extent. In addition, analysis of the public sequence from transcription data in an association panel suggested that natural variation in ZmIDH1.4 will be associated with drought tolerance in maize. These results suggested that ZmIDHs respond differently and/or redundantly to abiotic stresses during plant growth and development, and this analysis provides a foundation to understand how ZmIDHs respond to drought stress in maize.

Identification of genes associated with abiotic stress tolerance in sweetpotato using weighted gene co-expression network analysis

Sweetpotato, Ipomoea batatas (L.), a key food security crop, is negatively impacted by heat, drought, and salinity stress. The orange-fleshed sweetpotato cultivar "Beauregard" was exposed to heat, salt, and drought treatments for 24 and 48 h to identify genes responding to each stress condition in leaves. Analysis revealed both common (35 up regulated, 259 down regulated genes in the three stress conditions) and unique sets of up regulated (1337 genes by drought, 516 genes by heat, and 97 genes by salt stress) and down regulated (2445 genes by drought, 678 genes by heat, and 204 genes by salt stress) differentially expressed genes (DEGs) suggesting common, yet stress-specific transcriptional responses to these three abiotic stressors. Gene Ontology analysis of down regulated DEGs common to both heat and salt stress revealed enrichment of terms associated with "cell population proliferation" suggestive of an impact on the cell cycle by the two stress conditions. To identify shared and unique gene co-expression networks under multiple abiotic stress conditions, weighted gene co-expression network analysis was performed using gene expression profiles from heat, salt, and drought stress treated 'Beauregard' leaves yielding 18 co-expression modules. One module was enriched for "response to water deprivation," "response to abscisic acid," and "nitrate transport" indicating synergetic crosstalk between nitrogen, water, and phytohormones with genes encoding osmotin, cell expansion, and cell wall modification proteins present as key hub genes in this drought-associated module. This research lays the groundwork for exploring to a further degree, mechanisms for abiotic stress tolerance in sweetpotato.

The Genome-Wide Identification, Characterization, and Expression Analysis of the Strictosidine Synthase-like Family in Maize (Zea mays L.)

Maize is often subjected to various environmental stresses. The strictosidine synthase-like (SSL) family is thought to catalyze the key step in the monoterpene alkaloids synthesis pathway in response to environmental stresses. However, the role of ZmSSL genes in maize growth and development and its response to stresses is unknown. Herein, we undertook the systematic identification and analysis of maize SSL genes. Twenty SSL genes were identified in the maize genome. Except for chromosomes 3, 5, 6, and 10, they were unevenly distributed on the remaining 6 chromosomes. A total of 105 SSL genes from maize, sorghum, rice, Aegilops tauschii, and Arabidopsis were divided into five evolutionary groups, and ZmSSL gene structures and conserved protein motifs in the same group were similar. A collinearity analysis showed that tandem duplication plays an important role in the evolution of the SSL family in maize, and ZmSSL genes share more collinear genes in crops (maize, sorghum, rice, and Ae. tauschii) than in Arabidopsis. Cis-element analysis in the ZmSSL gene promoter region revealed that most genes contained many development and stress response elements. We evaluated the expression levels of ZmSSL genes under normal conditions and stress treatments. ZmSSL4-9 were widely expressed in different tissues and were positively or negatively regulated by heat, cold, and infection stress from Colletotrichum graminicola and Cercospora zeina. Moreover, ZmSSL4 and ZmSSL5 were localized in the chloroplast. Taken together, we provide insight into the evolutionary relationships of the ZmSSL genes, which would be useful to further identify the potential functions of ZmSSLs in maize.

Plant stem cells under low oxygen: metabolic rewiring by phytoglobin underlies stem cell functionality

Root growth in maize (Zea mays L.) is regulated by the activity of the quiescent center (QC) stem cells located within the root apical meristem. Here, we show that despite being highly hypoxic under normal oxygen tension, QC stem cells are vulnerable to hypoxic stress, which causes their degradation with subsequent inhibition of root growth. Under low oxygen, QC stem cells became depleted of starch and soluble sugars and exhibited reliance on glycolytic fermentation with the impairment of the TCA cycle through the depressed activity of several enzymes, including pyruvate dehydrogenase (PDH). This finding suggests that carbohydrate delivery from the shoot might be insufficient to meet the metabolic demand of QC stem cells during stress. Some metabolic changes characteristic of the hypoxic response in mature root cells were not observed in the QC. Hypoxia-responsive genes, such as PYRUVATE DECARBOXYLASE (PDC) and ALCOHOL DEHYDROGENASE (ADH), were not activated in response to hypoxia, despite an increase in ADH activity. Increases in phosphoenolpyruvate (PEP) with little change in steady-state levels of succinate were also atypical responses to low-oxygen tensions. Overexpression of PHYTOGLOBIN 1 (ZmPgb1.1) preserved the functionality of the QC stem cells during stress. The QC stem cell preservation was underpinned by extensive metabolic rewiring centered around activation of the TCA cycle and retention of carbohydrate storage products, denoting a more efficient energy production and diminished demand for carbohydrates under conditions where nutrient transport may be limiting. Overall, this study provides an overview of metabolic responses occurring in plant stem cells during oxygen deficiency.

Identification of southern corn rust resistance QTNs in Chinese summer maize germplasm via multi-locus GWAS and post-GWAS analysis

Southern corn rust (SCR) caused by Puccinia polysora Underw is a major disease leading to severe yield losses in China Summer Corn Belt. Using six multi-locus GWAS methods, we identified a set of SCR resistance QTNs from a diversity panel of 140 inbred lines collected from China Summer Corn Belt. Thirteen QTNs on chromosomes 1, 2, 4, 5, 6, and 8 were grouped into three types of allele effects and their associations with SCR phenotypes were verified by post-GWAS case-control sampling, allele/haplotype effect analysis. Relative resistance (RRR) and relative susceptibility (RRs) catering to its inbred carrier were estimated from single QTN and QTN-QTN combos and epistatitic effects were estimated for QTN-QTN combos. By transcriptomic annotation, a set of candidate genes were predicted to be involved in transcriptional regulation (S5_145, Zm00001d01613, transcription factor GTE4), phosphorylation (S8_123, Zm00001d010672, Pgk2- phosphoglycerate kinase 2), and temperature stress response (S6_164a/S6_164b, Zm00001d038806, hsp101, and S5_211, Zm00001d017978, cellulase25). The breeding implications of the above findings were discussed.

Deciphering shared attributes of plant long non-coding RNAs through a comparative computational approach

Over the past decade, long non-coding RNA (lncRNA), which lacks protein-coding potential, has emerged as an essential regulator of the genome. The present study examined 13,599 lncRNAs in Arabidopsis thaliana, 11,565 in Oryza sativa, and 32,397 in Zea mays for their characteristic features and explored the associated genomic and epigenomic features. We found lncRNAs were distributed throughout the chromosomes and the Helitron family of transposable elements (TEs) enriched, while the terminal inverted repeat depleted in lncRNA transcribing regions. Our analyses determined that lncRNA transcribing regions show rare or weak signals for most epigenetic marks except for H3K9me2 and cytosine methylation in all three plant species. LncRNAs showed preferential localization in the nucleus and cytoplasm; however, the distribution ratio in the cytoplasm and nucleus varies among the studied plant species. We identified several conserved endogenous target mimic sites in the lncRNAs among the studied plants. We found 233, 301, and 273 unique miRNAs, potentially targeting the lncRNAs of A. thaliana, O. sativa, and Z. mays, respectively. Our study has revealed that miRNAs, which interact with lncRNAs, target genes that are involved in a diverse array of biological and molecular processes. The miRNA-targeted lncRNAs displayed a strong affinity for several transcription factors, including ERF and BBR-BPC, mutually present in all three plants, advocating their conserved functions. Overall, the present study showed that plant lncRNAs exhibit conserved genomic and epigenomic characteristics and potentially govern the growth and development of plants.

Hybrid performance evaluation and genome-wide association analysis of root system architecture in a maize association population

KEY MESSAGE

The genetic architecture of RSA traits was dissected by GWAS and coexpression networks analysis in a maize association population. Root system architecture (RSA) is a crucial determinant of water and nutrient uptake efficiency in crops. However, the maize genetic architecture of RSA is still poorly understood due to the challenges in quantifying root traits and the lack of dense molecular markers. Here, an association mapping panel including 356 inbred lines were crossed with a common tester, Zheng58, and the test crosses were phenotyped for 12 RSA traits in three locations. We observed a 1.3 ~ sixfold phenotypic variation for measured RSA in the association panel. The association panel consisted of four subpopulations, non-stiff stalk (NSS) lines, stiff stalk (SS), tropical/subtropical (TST), and mixed. Zheng58 × TST has a 2.1% higher crown root number (CRN) and 8.6% less brace root number (BRN) than Zheng58 × NSS and Zheng58 × SS, respectively. Using a genome-wide association study (GWAS) with 1.25 million SNPs and correction for population structure, 191 significant SNPs were identified for root traits. Ninety (47%) of the significant SNPs showed positive allelic effects, and 101 (53%) showed negative effects. Each locus could explain 0.39% to 11.8% of phenotypic variation. By integrating GWAS results and comparing coexpression networks, 26 high-priority candidate genes were identified. Gene GRMZM2G377215, which belongs to the COBRA-like gene family, affected root growth and development. Gene GRMZM2G468657 encodes the aspartic proteinase nepenthesin-1, related to root development and N-deficient response. Collectively, our research provides progress in the genetic dissection of root system architecture. These findings present the further possibility for the genetic improvement of root traits in maize.

The genetic architecture of prolificacy in maize revealed by association mapping and bulk segregant analysis

KEY MESSAGE

Here, we revealed maize prolificacy highly correlated with domestication and identified a causal gene ZmEN1 located in one novel QTL qGEN261 that regulating maize prolificacy by using multiple-mapping methods. The development of maize prolificacy (EN) is crucial for enhancing yield and breeding specialty varieties. To achieve this goal, we employed a genome-wide association study (GWAS) to analyze the genetic architecture of EN in maize. Using 492 inbred lines with a wide range of EN variability, our results demonstrated significant differences in genetic, environmental, and interaction effects. The broad-sense heritability (H2) of EN was 0.60. Through GWAS, we identified 527 significant single nucleotide polymorphisms (SNPs), involved 290 quantitative trait loci (QTL) and 806 genes. Of these SNPs, 18 and 509 were classified as major effect loci and minor loci, respectively. In addition, we performed a bulk segregant analysis (BSA) in an F2 population constructed by a few-ears line Zheng58 and a multi-ears line 647. Our BSA results identified one significant QTL, qBEN1. Importantly, combining the GWAS and BSA, four co-located QTL, involving six genes, were identified. Three of them were expressed in vegetative meristem, shoot tip, internode and tip of ear primordium, with ZmEN1, encodes an unknown auxin-like protein, having the highest expression level in these tissues. It suggested that ZmEN1 plays a crucial role in promoting axillary bud and tillering to encourage the formation of prolificacy. Haplotype analysis of ZmEN1 revealed significant differences between different haplotypes, with inbred lines carrying hap6 having more EN. Overall, this is the first report about using GWAS and BSA to dissect the genetic architecture of EN in maize, which can be valuable for breeding specialty maize varieties and improving maize yield.

Identification and Characterization of VDAC Family in Maize

The voltage-dependent anion channel (VDAC) is the most abundant protein in the outer mitochondrial membrane (OMM) of all eukaryotes, having an important role in the communication between mitochondria and cytosol. The plant VDAC family consists of a wide variety of members that may participate in cell responses to several environmental stresses. However, there is no experimental information about the members comprising the maize VDAC (ZmVDAC) family. In this study, the ZmVDAC family was identified, and described, and its gene transcription profile was explored during the first six days of germination and under different biotic stress stimuli. Nine members were proposed as bona fide VDAC genes with a high potential to code functional VDAC proteins. Each member of the ZmVDAC family was characterized in silico, and nomenclature was proposed according to phylogenetic relationships. Transcript levels in coleoptiles showed a different pattern of expression for each ZmVDAC gene, suggesting specific roles for each one during seedling development. This expression profile changed under Fusarium verticillioides infection and salicylic acid, methyl jasmonate, and gibberellic acid treatments, suggesting no redundancy for the nine ZmVDAC genes and, thus, probably specific and diverse functions according to plant needs and environmental conditions. Nevertheless, ZmVDAC4b was significantly upregulated upon biotic stress signals, suggesting this gene's potential role during the biotic stress response.

Genes and genetic mechanisms contributing to fall armyworm resistance in maize

Maize (Zea mays L.) is a crop of major economic and food security importance globally. The fall armyworm (FAW), Spodoptera frugiperda, can devastate entire maize crops, especially in countries or markets that do not allow the use of transgenic crops. Host-plant insect resistance is an economical and environmentally benign way to control FAW, and this study sought to identify maize lines, genes, and pathways that contribute to resistance to FAW. Of the 289 maize lines phenotyped for FAW damage in artificially infested, replicated field trials over 3 years, 31 were identified with good levels of resistance that could donate FAW resistance into elite but susceptible hybrid parents. The 289 lines were genotyped by sequencing to provide single nucleotide polymorphism (SNP) markers for a genome-wide association study (GWAS), followed by a metabolic pathway analysis using the Pathway Association Study Tool (PAST). GWAS identified 15 SNPs linked to 7 genes, and PAST identified multiple pathways, associated with FAW damage. Top pathways, and thus useful resistance mechanisms for further study, include hormone signaling pathways and the biosynthesis of carotenoids (particularly zeaxanthin), chlorophyll compounds, cuticular wax, known antibiosis agents, and 1,4-dihydroxy-2-naphthoate. Targeted metabolite analysis confirmed that maize genotypes with lower levels of FAW damage tend to have higher levels of chlorophyll a than genotypes with high FAW damage, which tend to have lower levels of pheophytin, lutein, chlorophyll b and β-carotene. The list of resistant genotypes, and the results from the genetic, pathway, and metabolic study, can all contribute to efficient creation of FAW resistant cultivars.

Salt Tolerant Gene 1 contributes to salt tolerance by maintaining photosystem II activity in maize

Salt stress is a major environmental factor limiting crop growth and productivity. Here, we show that Salt-Tolerant Gene 1 (ZmSTG1) contributes to salt tolerance by maintaining photosystem activity in maize. ZmSTG1 encodes an endoplasmic reticulum localized protein and retrotransposon insertion in the promoter region causes differential expression levels in maize inbred lines. Overexpression of ZmSTG1 improved plant growth vigor, and knockout of ZmSTG1 weakened plant growth under normal and salt stress conditions. Transcriptome and metabolome analyses indicated that ZmSTG1 might regulate the expression of lipid trafficking-related genes dependent on the abscisic acid (ABA) signaling pathway, thereby increasing the galactolipids and phospholipid concentrations in the photosynthetic membrane under salt stress. Chlorophyll fluorescence parameters showed that the knockout of ZmSTG1 led to significant impairment of plant photosystem II (PSII) activity under normal and salt stress conditions, whereas overexpression of ZmSTG1 dramatically improved plant PSII activity under salt stress conditions. We also demonstrated that the application of the salt-tolerant locus could enhance salt tolerance in hybrid maize plants. Taken together, we propose that ZmSTG1 may modulate the lipid composition in the photosynthetic membrane by affecting the expression of lipid trafficking-related genes to maintain the photosynthetic activity of plants under salt stress.

MicroRNA408 negatively regulates salt tolerance by affecting secondary cell wall development in maize

Although microRNA408 (miR408) is a highly conserved miRNA, the miR408 response to salt stress differs among plant species. Here, we show that miR408 transcripts are strongly repressed by salt stress and methyl viologen treatment in maize (Zea mays). Application of N, N1-dimethylthiourea partly relieved the NaCl-induced down-regulation of miR408. Transgenic maize overexpressing MIR408b is hypersensitive to salt stress. Overexpression of MIR408b enhanced the rate of net Na+ efflux, caused Na+ to locate in the inter-cellular space, reduced lignin accumulation, and reduced the number of cells in vascular bundles under salt stress. We further demonstrated that miR408 targets ZmLACCASE9 (ZmLAC9). Knockout of MIR408a or MIR408b or overexpression of ZmLAC9 increased the accumulation of lignin, thickened the walls of pavement cells, and improved salt tolerance of maize. Transcriptome profiles of the wild-type and MIR408b-overexpressing transgenic maize with or without salt stress indicated that miR408 negatively regulates the expression of cell wall biogenesis genes under salt conditions. These results indicate that miR408 negatively regulates salt tolerance by regulating secondary cell wall development in maize.

A leucine-rich repeat receptor kinase gene confers quantitative susceptibility to maize southern leaf blight

Southern leaf blight (SLB), caused by the necrotrophic fungal pathogen Cochliobolus heterostrophus (anamorph Bipolaris maydis), is a major foliar disease which causes significant yield losses in maize worldwide. A major quantitative trait locus, qSLB_3.04 , conferring recessive resistance to SLB was previously mapped on maize chromosome 3. Using a combination of map-based cloning, association analysis, ethyl methanesulfonate and transposon mutagenesis, and CRISPR-Cas9 editing, we demonstrate that a leucine-rich repeat receptor-like kinase gene which we have called ChSK1 (Cochliobolus heterostrophus Susceptibility Kinase 1) at qSLB_3.04 causes increased susceptibility to SLB. Genes of this type have generally been associated with the defense response. We present evidence that ChSK1 may be associated with suppression of the basal immune response. These findings contribute to our understanding of plant disease susceptibility genes and the potential to use them for engineering durable disease resistance.

Maize Transcription Factor ZmHsf28 Positively Regulates Plant Drought Tolerance

Identification of central genes governing plant drought tolerance is fundamental to molecular breeding and crop improvement. Here, maize transcription factor ZmHsf28 is identified as a positive regulator of plant drought responses. ZmHsf28 exhibited inducible gene expression in response to drought and other abiotic stresses. Overexpression of ZmHsf28 diminished drought effects in Arabidopsis and maize. Gene silencing of ZmHsf28 via the technology of virus-induced gene silencing (VIGS) impaired maize drought tolerance. Overexpression of ZmHsf28 increased jasmonate (JA) and abscisic acid (ABA) production in transgenic maize and Arabidopsis by more than two times compared to wild-type plants under drought conditions, while it decreased reactive oxygen species (ROS) accumulation and elevated stomatal sensitivity significantly. Transcriptomic analysis revealed extensive gene regulation by ZmHsf28 with upregulation of JA and ABA biosynthesis genes, ROS scavenging genes, and other drought related genes. ABA treatment promoted ZmHsf28 regulation of downstream target genes. Specifically, electrophoretic mobility shift assays (EMSA) and yeast one-hybrid (Y1H) assay indicated that ZmHsf28 directly bound to the target gene promoters to regulate their gene expression. Taken together, our work provided new and solid evidence that ZmHsf28 improves drought tolerance both in the monocot maize and the dicot Arabidopsis through the implication of JA and ABA signaling and other signaling pathways, shedding light on molecular breeding for drought tolerance in maize and other crops.

Local diversity of drought resistance and resilience in Populus tomentosa correlates with the variation of DNA methylation

Little information is known about DNA methylation variation in shaping environment-specific drought resistance and resilience for tree adaptation. In this study, we leveraged RNA sequencing and whole-genome bisulfite sequencing data to dissect the distinction of epigenetic regulation under drought stress and rewater condition of Populus tomentosa accessions from three geographical regions. We demonstrated low resistance and high resilience for accessions from South. Non-CG methylation levels in promoter regions of Southern accessions were lower than accessions from higher latitudes and negatively regulated gene expression. CHH context methylation was more sensitive to drought stress, and the geographical-specific differentially methylated regions were scarcely changed by environmental fluctuation. We identified 60 conserved hub genes within the co-expression networks that correlate with photosynthetic and stomatal morphological traits. Epigenome-wide association studies and genome-wide association studies of these 60 hub genes revealed the interdependency between genetic and epigenetic variation in GATA9 and LECRK-VIII.2, which was associated with stomatal morphology and chlorophyll content. The natural epigenetic variation in GATA9 was also faithfully transmitted to progenies in two family-based F1 populations. This study indicates a functional relationship of DNA methylation diversity with drought resistance and resilience which offers new insights into plants' local adaptation to a stressful environment.

High temporal-resolution transcriptome landscapes of maize embryo sac and ovule during early seed development

Here we provided a high temporal-resolution transcriptome atlas of maize embryo sac and ovule to reveal the gene activity dynamic during early seed development. The early maize (Zea mays) seed development is initiated from double fertilization in the embryo sac and needs to undergo a highly dynamic and complex development process to form the differentiated embryo and endosperm. Despite the importance of maize seed for food, feed, and biofuel, many regulators responsible for controlling its early development are not known yet. Here, we reported a high temporal-resolution transcriptome atlas of embryo sac and ovule based on 44 time point samples collected within the first four days of seed development. A total of 25,187 genes including 1598 transcription factors (TFs) involved in early seed development were detected. Global comparisons of the expressions of these genes revealed five distinct development stages of early seed, which are mainly related to double fertilization, asymmetric cell division of the zygote, as well as coenocyte formation, cellularization and differentiation in endosperm. We identified 3327 seed-specific genes, which more than one thousand seed-specific genes with main expressions during early seed development were newly identified here, including 859 and 186 genes predominantly expressed in the embryo sac and ovule, respectively. Combined with the published transcriptome data of seed, we uncovered the dominant auxin biosynthesis, transport and signaling related genes at different development stages and subregions of seed. These results are helpful for understanding the genetic control of early seed development.

Combined genome-wide association study and gene co-expression network analysis identified ZmAKINβγ1 involved in lead tolerance and accumulation in maize seedlings

Lead (Pb) contamination has become an important abiotic stress that negatively influences crop biomass and yield, threatening human health via food chains. The excavation of causal genes for Pb tolerance in maize will contribute to the breeding of Pb-tolerant maize germplasms. This study aimed to demonstrate the effects of AKINbetagamma-1 protein kinase (ZmAKINβγ1) on maize tolerance to Pb and reveal its molecular mechanisms underlying Pb tolerance. ZmAKINβγ1 was identified using genome-wide association study and weighted gene co-expression network analysis for shoot dry weight (SDW) and root dry weight (RDW) under Pb treatment. The OE and RNAi experiments showed that ZmAKINβγ1 negatively regulated maize tolerance to Pb by reducing SDW and RDW and increasing Pb accumulation in maize. Comparative transcriptome analysis between the OE/RNAi and wild-type lines revealed that ZmAKINβγ1 participated in the pectin metabolism process and nitrogen compound response. Gene-based association analyses revealed that three variants located in ZmAKINβγ1 promoter induced changes in its expression and Pb tolerance among maize lines. The dual-luciferase reporter system verified that the two genotypes (AAT and CGG) of ZmAKINβγ1 promoter had contrasting transcriptional activities. Collectively, ZmAKINβγ1-mediated Pb tolerance provided new insights into the cultivation of Pb-tolerant maize varieties and phytoremediation of Pb-polluted soils.

The In Silico Characterization of Monocotyledonous α-l-Arabinofuranosidases on the Example of Maize

Plant α-l-arabinofuranosidases remove terminal arabinose from arabinose-containing substrates such as plant cell wall polysaccharides, including arabinoxylans, arabinogalactans, and arabinans. In plants, de-arabinosylation of cell wall polysaccharides accompanies different physiological processes such as fruit ripening and elongation growth. In this report, we address the diversity of plant α-l-arabinofuranosidases of the glycoside hydrolase (GH) family 51 through their phylogenetic analysis as well as their structural features. The CBM4-like domain at N-terminus was found to exist only in GH51 family proteins and was detected in almost 90% of plant sequences. This domain is similar to bacterial CBM4, but due to substitutions of key amino acid residues, it does not appear to be able to bind carbohydrates. Despite isoenzymes of GH51 being abundant, in particular in cereals, almost half of the GH51 proteins in Poales have a mutation of the acid/base residue in the catalytic site, making them potentially inactive. Open-source data on the transcription and translation of GH51 isoforms in maize were analyzed to discuss possible functions of individual isoenzymes. The results of homology modeling and molecular docking showed that the substrate binding site can accurately accommodate terminal arabinofuranose and that arabinoxylan is a more favorable ligand for all maize GH51 enzymes than arabinan.

Unleashing the power within short-read RNA-seq for plant research: Beyond differential expression analysis and toward regulomics

RNA-seq has become a state-of-the-art technique for transcriptomic studies. Advances in both RNA-seq techniques and the corresponding analysis tools and pipelines have unprecedently shaped our understanding in almost every aspects of plant sciences. Notably, the integration of huge amount of RNA-seq with other omic data sets in the model plants and major crop species have facilitated plant regulomics, while the RNA-seq analysis has still been primarily used for differential expression analysis in many less-studied plant species. To unleash the analytical power of RNA-seq in plant species, especially less-studied species and biomass crops, we summarize recent achievements of RNA-seq analysis in the major plant species and representative tools in the four types of application: (1) transcriptome assembly, (2) construction of expression atlas, (3) network analysis, and (4) structural alteration. We emphasize the importance of expression atlas, coexpression networks and predictions of gene regulatory relationships in moving plant transcriptomes toward regulomics, an omic view of genome-wide transcription regulation. We highlight what can be achieved in plant research with RNA-seq by introducing a list of representative RNA-seq analysis tools and resources that are developed for certain minor species or suitable for the analysis without species limitation. In summary, we provide an updated digest on RNA-seq tools, resources and the diverse applications for plant research, and our perspective on the power and challenges of short-read RNA-seq analysis from a regulomic point view. A full utilization of these fruitful RNA-seq resources will promote plant omic research to a higher level, especially in those less studied species.

A spatiotemporal transcriptomic network dynamically modulates stalk development in maize

Maize (Zea mays) is an important cereal crop with suitable stalk formation which is beneficial for acquiring an ideal agronomic trait to resist lodging and higher planting density. The elongation pattern of stalks arises from the variable growth of individual internodes driven by cell division and cell expansion comprising the maize stalk. However, the spatiotemporal dynamics and regulatory network of the maize stalk development and differentiation process remain unclear. Here, we report spatiotemporally resolved transcriptomes using all internodes of the whole stalks from developing maize at the elongation and maturation stages. We identified four distinct groups corresponding to four developmental zones and nine specific clusters with diverse spatiotemporal expression patterns among individual internodes of the stalk. Through weighted gene coexpression network analysis, we constructed transcriptional regulatory networks at a fine spatiotemporal resolution and uncovered key modules and candidate genes involved in internode maintenance, elongation, and division that determine stalk length and thickness in maize. Further CRISPR/Cas9-mediated knockout validated the function of a cytochrome P450 gene, ZmD1, in the regulation of stalk length and thickness as predicted by the WGCN. Collectively, these results provide insights into the high genetic complexity of stalk development and the potentially valuable resources with ideal stalk lengths and widths for genetic improvements in maize.

Phenotypic and transcriptomic responses of cultivated sunflower seedlings (Helianthus annuus L.) to four abiotic stresses

Plants encounter and respond to numerous abiotic stresses during their lifetimes. These stresses are often related and could therefore elicit related responses. There are, however, relatively few detailed comparisons between multiple different stresses at the molecular level. Here, we investigated the phenotypic and transcriptomic response of cultivated sunflower (Helianthus annuus L.) seedlings to three water-related stresses (i.e., dry-down, an osmotic challenge, and salt stress), as well as a generalized low-nutrient stress. All four stresses negatively impacted seedling growth, with the nutrient stress having a more divergent response from control as compared to the water-related stresses. Phenotypic responses were consistent with expectations for growth in low-resource environments, including increased (i.e., less negative) carbon fractionation values and leaf C:N ratios, as well as increased belowground biomass allocation. The number of differentially expressed genes (DEGs) under stress was greater in leaf tissue, but roots exhibited a higher proportion of DEGs unique to individual stresses. Overall, the three water-related stresses had a more similar transcriptomic response to each other vs. nutrient stress, though this pattern was more pronounced in root vs. leaf tissue. In contrast to our DEG analyses, co-expression network analysis revealed that there was little indication of a shared response between the four stresses in despite the majority of DEGs being shared between multiple stresses. Importantly, osmotic stress, which is often used to simulate drought stress in experimental settings, had little transcriptomic resemblance to true water limitation (i.e., dry-down) in our study, calling into question its utility as a means for simulating drought.

Smart breeding approaches in post-genomics era for developing climate-resilient food crops

Improving the crop traits is highly required for the development of superior crop varieties to deal with climate change and the associated abiotic and biotic stress challenges. Climate change-driven global warming can trigger higher insect pest pressures and plant diseases thus affecting crop production sternly. The traits controlling genes for stress or disease tolerance are economically imperative in crop plants. In this scenario, the extensive exploration of available wild, resistant or susceptible germplasms and unraveling the genetic diversity remains vital for breeding programs. The dawn of next-generation sequencing technologies and omics approaches has accelerated plant breeding by providing the genome sequences and transcriptomes of several plants. The availability of decoded plant genomes offers an opportunity at a glance to identify candidate genes, quantitative trait loci (QTLs), molecular markers, and genome-wide association studies that can potentially aid in high throughput marker-assisted breeding. In recent years genomics is coupled with marker-assisted breeding to unravel the mechanisms to harness better better crop yield and quality. In this review, we discuss the aspects of marker-assisted breeding and recent perspectives of breeding approaches in the era of genomics, bioinformatics, high-tech phonemics, genome editing, and new plant breeding technologies for crop improvement. In nutshell, the smart breeding toolkit in the post-genomics era can steadily help in developing climate-smart future food crops.

Comparative transcriptome meta-analysis reveals a set of genes involved in the responses to multiple pathogens in maize

Maize production is constantly threatened by the presence of different fungal pathogens worldwide. Genetic resistance is the most favorable approach to reducing yield losses resulted from fungal diseases. The molecular mechanism underlying disease resistance in maize remains largely unknown. The objective of this study was to identify key genes/pathways that are consistently associated with multiple fungal pathogen infections in maize. Here, we conducted a meta-analysis of gene expression profiles from seven publicly available RNA-seq datasets of different fungal pathogen infections in maize. We identified 267 common differentially expressed genes (co-DEGs) in the four maize leaf infection experiments and 115 co-DEGs in all the seven experiments. Functional enrichment analysis showed that the co-DEGs were mainly involved in the biosynthesis of diterpenoid and phenylpropanoid. Further investigation revealed a set of genes associated with terpenoid phytoalexin and lignin biosynthesis, as well as potential pattern recognition receptors and nutrient transporter genes, which were consistently up-regulated after inoculation with different pathogens. In addition, we constructed a weighted gene co-expression network and identified several hub genes encoding transcription factors and protein kinases. Our results provide valuable insights into the pathways and genes influenced by different fungal pathogens, which might facilitate mining multiple disease resistance genes in maize.

Variation in leaf transcriptome responses to elevated ozone corresponds with physiological sensitivity to ozone across maize inbred lines

We examine the impact of sustained elevated ozone concentration on the leaf transcriptome of 5 diverse maize inbred genotypes, which vary in physiological sensitivity to ozone (B73, Mo17, Hp301, C123, and NC338), using long reads to assemble transcripts and short reads to quantify expression of these transcripts. More than 99% of the long reads, 99% of the assembled transcripts, and 97% of the short reads map to both B73 and Mo17 reference genomes. Approximately 95% of the genes with assembled transcripts belong to known B73-Mo17 syntenic loci and 94% of genes with assembled transcripts are present in all temperate lines in the nested association mapping pan-genome. While there is limited evidence for alternative splicing in response to ozone stress, there is a difference in the magnitude of differential expression among the 5 genotypes. The transcriptional response to sustained ozone stress in the ozone resistant B73 genotype (151 genes) was modest, while more than 3,300 genes were significantly differentially expressed in the more sensitive NC338 genotype. There is the potential for tandem duplication in 30% of genes with assembled transcripts, but there is no obvious association between potential tandem duplication and differential expression. Genes with a common response across the 5 genotypes (83 genes) were associated with photosynthesis, in particular photosystem I. The functional annotation of genes not differentially expressed in B73 but responsive in the other 4 genotypes (789) identifies reactive oxygen species. This suggests that B73 has a different response to long-term ozone exposure than the other 4 genotypes. The relative magnitude of the genotypic response to ozone, and the enrichment analyses are consistent regardless of whether aligning short reads to: long read assembled transcripts; the B73 reference; the Mo17 reference. We find that prolonged ozone exposure directly impacts the photosynthetic machinery of the leaf.

COMPILE: a GWAS computational pipeline for gene discovery in complex genomes

BACKGROUND

Genome-Wide Association Studies (GWAS) are used to identify genes and alleles that contribute to quantitative traits in large and genetically diverse populations. However, traits with complex genetic architectures create an enormous computational load for discovery of candidate genes with acceptable statistical certainty. We developed a streamlined computational pipeline for GWAS (COMPILE) to accelerate identification and annotation of candidate maize genes associated with a quantitative trait, and then matches maize genes to their closest rice and Arabidopsis homologs by sequence similarity.

RESULTS

COMPILE executed GWAS using a Mixed Linear Model that incorporated, without compression, recent advancements in population structure control, then linked significant Quantitative Trait Loci (QTL) to candidate genes and RNA regulatory elements contained in any genome. COMPILE was validated using published data to identify QTL associated with the traits of α-tocopherol biosynthesis and flowering time, and identified published candidate genes as well as additional genes and non-coding RNAs. We then applied COMPILE to 274 genotypes of the maize Goodman Association Panel to identify candidate loci contributing to resistance of maize stems to penetration by larvae of the European Corn Borer (Ostrinia nubilalis). Candidate genes included those that encode a gene of unknown function, WRKY and MYB-like transcriptional factors, receptor-kinase signaling, riboflavin synthesis, nucleotide-sugar interconversion, and prolyl hydroxylation. Expression of the gene of unknown function has been associated with pathogen stress in maize and in rice homologs closest in sequence identity.

CONCLUSIONS

The relative speed of data analysis using COMPILE allowed comparison of population size and compression. Limitations in population size and diversity are major constraints for a trait and are not overcome by increasing marker density. COMPILE is customizable and is readily adaptable for application to species with robust genomic and proteome databases.

Dissecting the Genetic Structure of Maize Leaf Sheaths at Seedling Stage by Image-Based High-Throughput Phenotypic Acquisition and Characterization

The rapid development of high-throughput phenotypic detection techniques makes it possible to obtain a large number of crop phenotypic information quickly, efficiently, and accurately. Among them, image-based phenotypic acquisition method has been widely used in crop phenotypic identification and characteristic research due to its characteristics of automation, non-invasive, non-destructive and high throughput. In this study, we proposed a method to define and analyze the traits related to leaf sheaths including morphology-related, color-related and biomass-related traits at V6 stage. Next, we analyzed the phenotypic variation of leaf sheaths of 418 maize inbred lines based on 87 leaf sheath-related phenotypic traits. In order to further analyze the mechanism of leaf sheath phenotype formation, 25 key traits (2 biomass-related, 19 morphology-related and 4 color-related traits) with heritability greater than 0.3 were analyzed by genome-wide association studies (GWAS). And 1816 candidate genes of 17 whole plant leaf sheath traits and 1,297 candidate genes of 8 sixth leaf sheath traits were obtained, respectively. Among them, 46 genes with clear functional descriptions were annotated by single nucleotide polymorphism (SNPs) that both Top1 and multi-method validated. Functional enrichment analysis results showed that candidate genes of leaf sheath traits were enriched into multiple pathways related to cellular component assembly and organization, cell proliferation and epidermal cell differentiation, and response to hunger, nutrition and extracellular stimulation. The results presented here are helpful to further understand phenotypic traits of maize leaf sheath and provide a reference for revealing the genetic mechanism of maize leaf sheath phenotype formation.

Underlying mechanism of accelerated cell death and multiple disease resistance in a maize lethal leaf spot 1 allele

Multiple disease resistance (MDR) in maize has attracted increasing attention. However, the interplay between cell death and metabolite changes and their contributions to MDR remains elusive in maize. In this study, we identified a mutant named as lesion mimic 30 (les30) that showed 'suicidal' lesion formation in the absence of disease and had enhanced resistance to the fungal pathogen Curvularia lunata. Using map-based cloning, we identified the causal gene encoding pheophorbide a oxidase (PAO), which is known to be involved in chlorophyll degradation and MDR, and is encoded by LETHAL LEAF SPOT1 (LLS1). LLS1 was found to be induced by both biotic and abiotic stresses. Transcriptomics analysis showed that genes involved in defense responses and secondary metabolite biosynthesis were mildly activated in leaves of the les30 mutant without lesions, whilst they were strongly activated in leaves with lesions. In addition, in les30 leaves with lesions, there was overaccumulation of defense-associated phytohormones including jasmonic acid and salicylic acid, and of phytoalexins including phenylpropanoids, lignin, and flavonoids, suggesting that their biosynthesis was activated in a lesion-dependent manner. Taken together, our study implies the existence of an interactive amplification loop of interrupted chlorophyll degradation, cell death, expression of defense-related genes, and metabolite changes that results in suicidal lesion formation and MDR, and this has the potential to be exploited by genetic manipulation to improve maize disease resistance.

Meta-QTL analysis and candidate genes identification for various abiotic stresses in maize (Zea mays L.) and their implications in breeding programs

Global climate change leads to the concurrence of a number of abiotic stresses including moisture stress (drought, waterlogging), temperature stress (heat, cold), and salinity stress, which are the major factors affecting maize production. To develop abiotic stress tolerance in maize, many quantitative trait loci (QTL) have been identified, but very few of them have been utilized successfully in breeding programs. In this context, the meta-QTL analysis of the reported QTL will enable the identification of stable/real QTL which will pave a reliable way to introgress these QTL into elite cultivars through marker-assisted selection. In this study, a total of 542 QTL were summarized from 33 published studies for tolerance to different abiotic stresses in maize to conduct meta-QTL analysis using BiomercatorV4.2.3. Among those, only 244 major QTL with more than 10% phenotypic variance were preferably utilised to carry out meta-QTL analysis. In total, 32 meta-QTL possessing 1907 candidate genes were detected for different abiotic stresses over diverse genetic and environmental backgrounds. The MQTL2.1, 5.1, 5.2, 5.6, 7.1, 9.1, and 9.2 control different stress-related traits for combined abiotic stress tolerance. The candidate genes for important transcription factor families such as ERF, MYB, bZIP, bHLH, NAC, LRR, ZF, MAPK, HSP, peroxidase, and WRKY have been detected for different stress tolerances. The identified meta-QTL are valuable for future climate-resilient maize breeding programs and functional validation of candidate genes studies, which will help to deepen our understanding of the complexity of these abiotic stresses.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01294-9.

The Combination of Conventional QTL Analysis, Bulked-Segregant Analysis, and RNA-Sequencing Provide New Genetic Insights into Maize Mesocotyl Elongation under Multiple Deep-Seeding Environments

Mesocotyl length (MES) is an important trait that affects the emergence of maize seedlings after deep-seeding and is closely associated with abiotic stress. The elucidation of constitutive-QTLs (cQTLs) and candidate genes for MES and tightly molecular markers are thus of great importance in marker-assisted selection (MAS) breeding. Therefore, the objective of this study was to perform detailed genetic analysis of maize MES across 346 F_2:3 families, 30/30 extreme bulks of an F₂ population, and two parents by conventional QTL analysis, bulked-segregation analysis (BSA), and RNA-sequencing when maize was sown at the depths of 3, 15, and 20 cm, respectively. QTL analysis identified four major QTLs in Bin 1.09, Bin 3.04, Bin 4.06–4.07, and Bin 6.01 under two or more environments, which explained 2.89–13.97% of the phenotypic variance within a single environment. BSA results revealed the presence of seven significantly linked SNP/InDel regions on chromosomes 1 and 4, and six SNP/InDel regions and the major QTL of qMES4-1 overlapped and formed a cQTL, cQMES4, within the 160.98–176.22 Mb region. In total, 18,001 differentially expressed genes (DEGs) were identified across two parents by RNA-sequencing, and 24 of these genes were conserved core DEGs. Finally, we validated 15 candidate genes in cQMES4 to involve in cell wall structure, lignin biosyntheis, phytohormones (auxin, abscisic acid, brassinosteroid) signal transduction, circadian clock, and plant organ formation and development. Our findings provide a basis for MAS breeding and enhance our understanding of the deep-seeding tolerance of maize.

Impact of arbuscular mycorrhiza on maize P1B-ATPases gene expression and ionome in copper-contaminated soils

Arbuscular mycorrhizal (AM) fungi, symbionts of most land plants, increase plant fitness in metal contaminated soils. To further understand the mechanisms of metal tolerance in the AM symbiosis, the expression patterns of the maize Heavy Metal ATPase (HMA) family members and the ionomes of non-mycorrhizal and mycorrhizal plants grown under different Cu supplies were examined. Expression of ZmHMA5a and ZmHMA5b, whose encoded proteins were predicted to be localized at the plasma membrane, was up-regulated by Cu in non-mycorrhizal roots and to a lower extent in mycorrhizal roots. Gene expression of the tonoplast ZmHMA3a and ZmHMA4 isoforms was up-regulated by Cu-toxicity in shoots and roots of mycorrhizal plants. AM mitigates the changes induced by Cu toxicity on the maize ionome, specially at the highest Cu soil concentration. Altogether these data suggest that in Cu-contaminated soils, AM increases expression of the HMA genes putatively encoding proteins involved in Cu detoxification and balances mineral nutrient uptake improving the nutritional status of the maize plants.

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.

Making a Pangenome Using the Iterative Mapping Approach

Pangenomes have replaced single reference genomes as genetic references, as they contain a better scope of the diversity found in a single species. This protocol outlines the iterative mapping approach in constructing a pangenome, including how to check the raw data, align the data to a reference, how to assemble the data, and how to remove potential contaminants from the final assembly.

Trihelix Transcription Factor ZmThx20 Is Required for Kernel Development in Maize

Maize kernels are the harvested portion of the plant and are related to the yield and quality of maize. The endosperm of maize is a large storage organ that constitutes 80–90% of the dry weight of mature kernels. Maize kernels have long been the study of cereal grain development to increase yield. In this study, a natural mutation that causes abnormal kernel development, and displays a shrunken kernel phenotype, was identified and named “shrunken 2008 (sh2008)”. The starch grains in sh2008 are loose and have a less proteinaceous matrix surrounding them. The total storage protein and the major storage protein zeins are ~70% of that in the wild-type control (WT); in particular, the 19 kDa and 22 kDa α-zeins. Map-based cloning revealed that sh2008 encodes a GT-2 trihelix transcription factor, ZmThx20. Using CRISPR/Cas9, two other alleles with mutated ZmThx20 were found to have the same abnormal kernel. Shrunken kernels can be rescued by overexpressing normal ZmThx20. Comparative transcriptome analysis of the kernels from sh2008 and WT showed that the GO terms of translation, ribosome, and nutrient reservoir activity were enriched in the down-regulated genes (sh2008/WT). In short, these changes can lead to defects in endosperm development and storage reserve filling in seeds.

Transcriptome Reveals Allele Contribution to Heterosis in Maize

Heterosis, which has greatly increased maize yields, is associated with gene expression patterns during key developmental stages that enhance hybrid phenotypes relative to parental phenotypes. Before heterosis can be more effectively used for crop improvement, hybrid maize developmental gene expression patterns must be better understood. Here, six maize hybrids, including the popular hybrid Zhengdan958 (ZC) from China, were studied. Maize hybrids created in-house were generated using an incomplete diallel cross (NCII)-based strategy from four elite inbred parental lines. Differential gene expression (DEG) profiles corresponding to three developmental stages revealed that hybrid partial expression patterns exhibited complementarity of expression of certain parental genes, with parental allelic expression patterns varying both qualitatively and quantitatively in hybrids. Single-parent expression (SPE) and parent-specific expression (PSE) types of qualitative variation were most prevalent, 43.73 and 41.07% of variation, respectively. Meanwhile, negative super-dominance (NSD) and positive super-dominance (PSD) types of quantitative variation were most prevalent, 31.06 and 24.30% of variation, respectively. During the early reproductive growth stage, the gene expression pattern differed markedly from other developmental stage patterns, with allelic expression patterns during seed development skewed toward low-value parental alleles in hybrid seeds exhibiting significant quantitative variation-associated superiority. Comparisons of qualitative gene expression variation rates between ZC and other hybrids revealed proportions of SPE-DEGs (41.36%) in ZC seed DEGs that significantly exceeded the average proportion of SPE-DEGs found in seeds of other hybrids (28.36%). Importantly, quantitative gene expression variation rate comparisons between ZC and hybrids, except for transgressive expression, revealed that the ZC rate exceeded the average rate for other hybrids, highlighting the importance of partial gene expression in heterosis. Moreover, enriched ZC DEGs exhibiting distinct tissue-specific expression patterns belonged to four biological pathways, including photosynthesis, plant hormone signal transduction, biology metabolism and biosynthesis. These results provide valuable technical insights for creating hybrids exhibiting strong heterosis.

Maize unstable factor for orange1 is essential for endosperm development and carbohydrate accumulation

Maize (Zea mays L.) Ufo1-1 is a spontaneous dominant mutation of the unstable factor for orange1 (ufo1). We recently cloned ufo1, which is a Poaceae-specific gene highly expressed during seed development in maize. Here, we have characterized Ufo1-1 and a loss-of-function Ds insertion allele (ufo1-Dsg) to decipher the role of ufo1 in maize. We found that both ufo1 mutant alleles impact sugars and hormones, and have defects in the basal endosperm transfer layer (BETL) and adjacent cell types. The Ufo1-1 BETL had reduced cell elongation and cell wall ingrowth, resulting in cuboidal shaped transfer cells. In contrast, the ufo1-Dsg BETL cells showed a reduced overall size with abnormal wall ingrowth. Expression analysis identified the impact of ufo1 on several genes essential for BETL development. The overexpression of Ufo1-1 in various tissues leads to ectopic phenotypes, including abnormal cell organization and stomata subsidiary cell defects. Interestingly, pericarp and leaf transcriptomes also showed that as compared with wild type, Ufo1-1 had ectopic expression of endosperm development-specific genes. This study shows that Ufo1-1 impacts the expression patterns of a wide range of genes involved in various developmental processes.