Integration of omic networks in a developmental atlas of maize

Population genomics uncovers loci for trait improvement in the indigenous African cereal tef (Eragrostis tef)

Tef (Eragrostis tef) is an indigenous African cereal that is gaining global attention as a gluten-free "superfood" with high protein, mineral, and fibre contents. However, tef yields are limited by lodging and by losses during harvest owing to its small grain size (150× lighter than wheat). Breeders must also consider a strong cultural preference for white-grained over brown-grained varieties. Tef is relatively understudied with limited "omics" resources. Here, we resequence 220 tef accessions from an Ethiopian diversity collection and also perform multi-locational phenotyping for 25 agronomic and grain traits. Grain metabolome profiling reveals differential accumulation of fatty acids and flavonoids between white and brown grains. k-mer and SNP-based genome-wide association uncover important marker-trait associations, including a significant 70 kb peak for panicle morphology containing the tef orthologue of rice qSH1-a transcription factor regulating inflorescence morphology in cereals. We also observe a previously unknown relationship between grain size, colour, and fatty acids. These traits are highly associated with retrotransposon insertions in homoeologues of TRANSPARENT TESTA 2, a known regulator of grain colour. Our study provides valuable resources for tef research and breeding, facilitating the development of improved cultivars with desirable agronomic and nutritional properties.

A candidate association study of transcription factors in maize revealed the ZmPLATZ15-ZmEREB200 module as a key regulator of waterlogging tolerance at the seedling stage

Soil waterlogging is a major abiotic stress that severely impairs maize growth and development by inducing hypoxic conditions that disrupt essential physiological processes. Transcription factors (TFs) play crucial roles in modulating plant responses to waterlogging stress by regulating the expression of stress-related genes that enhance or diminish stress tolerance. In this study, we conducted an association analysis to identify 11 TFs closely associated with waterlogging stress in maize. Notably, the PLATZ family emerged as a novel and significant contributor to waterlogging stress. Overexpression of ZmPLATZ15 resulted in increased sensitivity to waterlogging at maize seedlings. Conversely, ZmEREB200, a member of the maize Group VII ERF (ZmERFVII) family, was significantly downregulated in the ZmPLATZ15 overexpression lines under waterlogging stress. Promoter analysis revealed that ZmPLATZ15 regulates ZmEREB200 by binding to the A/T-rich motifs in the ZmEREB200 promoter. Interestingly, overexpression of ZmEREB200 was found to enhance waterlogging tolerance at maize seedlings. To further elucidate their roles, we analyzed the transcriptomic profiles of ZmPLATZ15 and ZmEREB200 overexpression lines under waterlogging stress. The overlapping differentially expressed genes in both ZmPLATZ15 and ZmEREB200 overexpression lines were significantly enriched in pathways associated with redox balance and salicylic acid metabolism, both of which are crucial for modulating waterlogging tolerance at maize seedlings. Metabolomic analysis revealed that antioxidant enzyme activity, salicylic acid, and glutathione levels were decreased in OE-ZmPLATZ15, while these metabolites were significantly increased in OE-ZmEREB200. These contrasting metabolic responses in overexpression lines may underlie their different tolerances to waterlogging stress. Our findings provide valuable insights into the regulatory mechanisms underlying maize's response to waterlogging stress and highlight the potential of TFs as tools for developing maize varieties with enhanced waterlogging tolerance.

Mitochondrial Energy Homeostasis and Membrane Interaction Regulate the Rapid Growth of Moso Bamboo

The rapid growth of moso bamboo is primarily attributed to the swift elongation of its internodes. While mitochondria are known to provide energy for various cellular processes, the specific mechanisms by which they facilitate rapid growth in bamboo remain elusive. In this study, we optimised the procedures for mitochondria isolation and performed a comprehensive analysis of mitochondrial dynamics and proteomics from internodes at various growth stages, including the initial growth (IG) stage, the starting of cell division (SD), and the rapid elongation (RE). Confocal observation demonstrated that cells in the RE stage have a higher mitochondrial density and increased mitochondrial motility compared to other stages. Proteomic analysis of isolated mitochondria revealed an upregulation of the tricarboxylic acid cycle, along with a synchronous increase in both mitochondrial- and nuclear-encoded components of oxidative phosphorylation in RE cells. Moreover, the upregulation of various mitochondrial membrane transporters in RE cells suggests an enhanced exchange of metabolic intermediates and inorganic ions with the cytosol. Intriguingly, ultrastructural analysis and pharmacological treatments revealed membrane interactions between the endoplasmic reticulum (ER) and mitochondria in RE cells. In conclusion, our study provides novel insights into mitochondrial function and the intracellular dynamics that regulate the rapid growth of moso bamboo.

The combined effect of decreased stomatal density and aperture increases water use efficiency in maize

Stomata play a crucial role in balancing carbon dioxide uptake and water vapor loss, thereby regulating plant water use efficiency (WUE). Enhancing WUE is important for sustainable agriculture and food security, particularly for crops such as maize (Zea mays L.), as climate change and growing global food demand exacerbate limitations on water availability. Genetic factors controlling stomatal density and levels of the plant hormone abscisic acid (ABA) in leaves, which affect stomatal aperture, are key determinants of stomatal conductance (gs) and intrinsic WUE (iWUE). In this study, we demonstrate that stomatal density and stomatal aperture have a combined effect on gs and iWUE in maize. Using near-isogenic lines (NILs) and CRISPR/Cas9 mutants, we show that combining reduced stomatal density and reduced stomatal aperture can improve iWUE without compromising photosynthesis. This effect is pronounced at both, optimal and high temperatures. These findings highlight the potential of targeting multiple stomatal traits through genetic stacking to enhance WUE, offering a promising strategy for crop adaptation to water-limited environments.

Fertilization-induced synergid cell death by RALF12-triggered ROS production and ethylene signaling

Fertilization-dependent elimination of the persistent synergid cell is essential to block supernumerary pollen tubes and thus to avoid polyspermy in flowering plants. Little is known about the molecular mechanisms ensuring timely induction and execution of synergid cell death. We analyzed manually isolated maize synergid cells along their degeneration and show that they are gland cells expressing batteries of genes encoding small secreted proteins under control of the MYB98 transcription factor. This network is down-regulated after fertilization, while genes involved in reactive oxygen species (ROS) production, ethylene biosynthesis and response, senescence, and oxidative stress regulation are induced before synergid elimination and its ultimate fusion with the endosperm. We further show that the fertilization-induced RALF12 peptide specifically triggers mitochondrial ROS and apoptosis, while ethylene promotes synergid degeneration. In conclusion, this study sheds light on developmental programmed cell death (dPCD) in plants and provides a unique resource to discover unknown PCD regulators.

Endothelial transcriptomic, epigenomic and proteomic data challenge the proposed role for TSAd in vascular permeability

The vascular endothelial growth factor VEGF drives excessive vascular permeability to cause tissue-damaging oedema in neovascular and inflammatory diseases across multiple organs. Several molecular pathways have been implicated in VEGF-induced hyperpermeability, including binding of the VEGF-activated tyrosine kinase receptor VEGFR2 by the T-cell specific adaptor (TSAd) to recruit a SRC family kinase to induce junction opening between vascular endothelial cells (ECs). Inconsistent with a universal role for TSAd in permeability signalling, immunostaining approaches previously reported TSAd only in dermal and kidney vasculature. To address this discrepancy, we have mined publicly available omics data for expression of TSAd and other permeability-relevant signal transducers in multiple organs affected by VEGF-induced vascular permeability. Unexpectedly, TSAd transcripts were largely absent from EC single cell RNAseq data, whereas transcripts for other permeability-relevant signal transducers were detected readily. TSAd transcripts were also lacking from half of the EC bulk RNAseq datasets examined, and in the remaining datasets appeared at low levels concordant with models of leaky transcription. Epigenomic EC data located the TSAd promoter to closed chromatin in ECs, and mass spectrometry-derived EC proteomes typically lacked TSAd. By suggesting that TSAd is not actively expressed in ECs, our findings imply that TSAd is likely not critical for linking VEGFR2 to downstream signal transducers for EC junction opening.

Integration of multi-omics data accelerates molecular analysis of common wheat traits

Integration of multi-omics data can provide information on biomolecules from different layers to illustrate the complex biology systematically. Here, we build a multi-omics atlas containing 132,570 transcripts, 44,473 proteins, 19,970 phosphoproteins, and 12,427 acetylproteins across wheat vegetative and reproductive phases. Using this atlas, we elucidate transcriptional regulation network, contributions of post-translational modification (PTM) and transcript level to protein abundance, and biased homoeolog expression and PTM in wheat. The genes/proteins related to wheat development and disease resistance are systematically analyzed, thus identifying phosphorylation and/or acetylation modifications for the seed proteins controlling wheat grain quality and the disease resistance-related genes. Lastly, a unique protein module TaHDA9-TaP5CS1, specifying de-acetylation of TaP5CS1 by TaHDA9, is discovered, which regulates wheat resistance to Fusarium crown rot via increasing proline content. Our atlas holds great promise for fast-tracking molecular biology and breeding studies in wheat and related crops.

Ferredoxin-mediated mechanism for efficient nitrogen utilization in maize

Nitrogen (N) is an essential macronutrient for plant development and, ultimately, yield. Identifying the genetic components and mechanisms underlying N use efficiency in maize (Zea mays L.) is thus of great importance. Nitrate (NO3-) is the preferred inorganic N source in maize. Here we performed a genome-wide association study of shoot NO3- accumulation in maize seedlings grown under low-NO3- conditions, identifying the ferredoxin family gene ZmFd4 as a major contributor to this trait. ZmFd4 interacts and co-localizes with nitrite reductases (ZmNiRs) in chloroplasts to promote their enzymatic activity. Furthermore, ZmFd4 forms a high-affinity heterodimer with its closest paralogue, ZmFd9, in a NO3--sensitive manner. Although ZmFd4 exerts similar biochemical functions as ZmFd9, ZmFd4 and ZmFd9 interaction limits their ability to associate with ZmNiRs and stimulate their activity. Knockout lines for ZmFd4 with decreased NO3- contents exhibit more efficient NO3- assimilation, and field experiments show consistently improved N utilization and grain yield under N-deficient conditions. Our work thus provides molecular and mechanistic insights into the natural variation in N utilization, instrumental for genetic improvement of yield in maize and, potentially, in other crops.

Maize mutants in miR394-regulated genes show improved drought tolerance

Water limitation represents one of the major threats to agricultural production, which often leads to drought stress and results in compromised growth, development and yield of crop species. Drought tolerance has been intensively studied in search of potential targets for molecular approaches to crop improvement. However, drought adaptive traits are complex, and our understanding of the physiological and genetic basis of drought tolerance is still incomplete. The miR394-LCR pathway is a conserved regulatory module shown to participate in several aspects of plant growth and development, including stress response. Here, we characterized the miR394 pathway in maize, which harbours two genetic loci producing an evolutionarily conserved mature zma-miR394 targeting two transcripts coding for F-Box proteins, named hereby ZmLCR1 and ZmLCR2. Arabidopsis plants overexpressing the zma-MIR394B gene showed high tolerance to drought conditions compared to control plants. Moreover, analysis of the growth and development of single and double maize mutant plants in ZmLCR genes indicate that these mutations do not affect plant fitness when they grow in normal watering conditions, but mutants showed better survival than wild-type plants under water deprivation conditions. This increased drought tolerance is based on more efficient intrinsic water use, changes in root architecture and increased epicuticular wax content under water-limiting conditions. Our results indicate that the miR394-regulated ZmLCR genes are involved in drought stress tolerance and are remarkable candidates for maize crop improvement.

Structural Variants Contribute to Phenotypic Variation in Maize

Comprehensively identifying the loci shaping trait variation has been challenging, in part because standard approaches often miss many types of genetic variants. Structural variants (SVs), especially transposable elements (TEs), are likely to affect phenotypic variation but we lack methods that can detect polymorphic SVs and TEs using short-read sequencing data. Here, we used a whole genome alignment between two maize genotypes to identify polymorphic SVs and then genotyped a large maize diversity panel for these variants using short-read sequencing data. After characterising SV variation in the panel, we identified SV polymorphisms that are associated with life history traits and genotype-by-environment (GxE) interactions. While most of the SVs associated with traits contained TEs, only two of the SVs had boundaries that clearly matched TE breakpoints indicative of a TE insertion, while the other polymorphisms were likely caused by deletions. One of the SVs that appeared to be caused by a TE insertion had the most associations with gene expression compared to other trait-associated SVs. All of the SVs associated with traits were in linkage disequilibrium with nearby single nucleotide polymorphisms (SNPs), suggesting that the approach used here did not identify unique associations that would have been missed in a SNP association study. Overall, we have (1) created a technique to genotype SV polymorphisms across a large diversity panel using support from genomic short-read sequencing alignments and (2) connected this presence/absence SV variation to diverse traits and GxE interactions.

Pentyl leaf volatiles promote insect and pathogen resistance via enhancing ketol-mediated defense responses

Plants emit an array of volatile organic compounds in response to stresses. Six-carbon green leaf volatiles (GLVs) and five-carbon pentyl leaf volatiles (PLVs) are fatty acid-derived compounds involved in intra- and inter-species communications. Unlike extensively studied GLVs, the biological activities of PLVs remain understudied. Maize (Zea mays L.) contains a unique monocot-specific lipoxygenase, ZmLOX6, that is unable to oxidize fatty acids and instead possesses a hydroperoxide lyase-like activity to specifically produce PLVs. Here, we show that disruption of ZmLOX6 reduced resistance to fall armyworm (FAW; Spodoptera frugiperda) and fungal pathogens Colletotrichum graminicola and Cochliobolus heterostrophus. Metabolite profiling revealed that reduced resistance to insects and pathogens was associated with decreased production of PLVs and ketols, including the better studied α-ketol, 9,10-KODA (9-hydroxy-10-oxo-12(Z)-octadecadienoic acid). Exogenous PLV and 9,10-KODA treatments rescued the resistance of lox6 mutants to FAW and the pathogens. Surprisingly, the susceptible-to-herbivory lox6 mutants produced greater levels of wound-induced jasmonates, suggesting potential substrate competition between JA and PLV pathway branches and highlighting a strong role of PLVs in defense against insects. Similarly, likely due to substrate competition between GLV and PLV synthesis pathways, in response to C. graminicola infection, lox6 mutants accumulated elevated levels of GLVs, which promote susceptibility to this pathogen. Mutation of the GLV-producing ZmLOX10 in the lox6 mutant background reversed the susceptibility to C. graminicola, unveiling the contrasting roles of PLVs and GLVs in resistance to this pathogen. Overall, this study uncovered a potent signaling role of PLVs in defense against insect herbivory and fungal pathogens with distinct lifestyles.

Maize unstable factor for orange1 encodes a nuclear protein that affects redox accumulation during kernel development

The basal endosperm transfer layer (BETL) of the maize (Zea mays L.) kernel is composed of transfer cells for nutrient transport to nourish the developing kernel. To understand the spatiotemporal processes required for BETL development, we characterized 2 unstable factor for orange1 (Zmufo1) mutant alleles. The BETL defects in these mutants were associated with high levels of reactive oxygen species, oxidative DNA damage, and cell death. Interestingly, antioxidant supplementation in in vitro cultured kernels alleviated the cellular defects in mutants. Transcriptome analysis of the loss-of-function Zmufo1 allele showed differential expression of tricarboxylic acid cycle, redox homeostasis, and BETL-related genes. The basal endosperms of the mutant alleles had high levels of acetyl-CoA and elevated histone acetyltransferase activity. The BETL cell nuclei showed reduced electron-dense regions, indicating sparse heterochromatin distribution in the mutants compared with wild-type. Zmufo1 overexpression further reduced histone methylation marks in the enhancer and gene body regions of the pericarp color1 (Zmp1) reporter gene. Zmufo1 encodes an intrinsically disordered nuclear protein with very low sequence similarity to known proteins. Yeast two-hybrid and luciferase complementation assays established that ZmUFO1 interacts with proteins that play a role in chromatin remodeling, nuclear transport, and transcriptional regulation. This study establishes the critical function of Zmufo1 during basal endosperm development in maize kernels.

Structural variants contribute to phenotypic variation in maize

Comprehensively identifying the loci shaping trait variation has been challenging, in part because standard approaches often miss many types of genetic variants. Structural variants (SVs), especially transposable elements (TEs), are likely to affect phenotypic variation but we lack methods that can detect polymorphic structural variants and TEs using short-read sequencing data. Here, we used a whole genome alignment between two maize genotypes to identify polymorphic structural variants and then genotyped a large maize diversity panel for these variants using short-read sequencing data. After characterizing SV variation in the panel, we identified SV polymorphisms that are associated with life history traits and genotype-by-environment (GxE) interactions. While most of the SVs associated with traits contained TEs, only two of the SVs had boundaries that clearly matched TE breakpoints indicative of a TE insertion, while the other polymorphisms were likely caused by deletions. One of the SVs that appeared to be caused by a TE insertion had the most associations with gene expression compared to other trait-associated SVs. All of the SVs associated with traits were in linkage disequilibrium with nearby single nucleotide polymorphisms (SNPs), suggesting that the approach used here did not identify unique associations that would have been missed in a SNP association study. Overall, we have created a technique to genotype SV polymorphisms across a large diversity panel using support from genomic short-read sequencing alignments and connecting this presence/absence SV variation to diverse traits and GxE interactions.

The C2H2-type zinc finger transcription factor ZmDi19-7 regulates plant height and organ size by promoting cell size in maize

The drought-induced protein 19 (Di19) gene family encodes a Cys2/His2 zinc-finger protein implicated in responses to diverse plant stressors. To date, potential roles of these proteins as transcription factors remain largely elusive in maize. Here, we show that ZmDi19-7 gene exerts pivotal functions in regulation of plant height and organ growth by modulating the cell size in maize. ZmDi19-7 physically interacts with ubiquitin receptor protein ZmDAR1b, which is indispensable in ubiquitination of ZmDi19-7 and affects its protein stability. Further genetic analysis demonstrated that ZmDAR1b act in a common pathway with ZmDi19-7 to regulate cell size in maize. ZmDi19-7, severing as a transcriptional factor, is significantly enriched in conserved DiBS element in the promoter region of ZmHSP22, ZmHSP18c, ZmSAUR25, ZmSAUR55, ZmSAUR7 and ZmXTH23 and orchestrates the expression of these genes involving in auxin-mediated cell expansion and protein processing in the endoplasmic reticulum. Thus, our findings demonstrate that ZmDi19-7 is an important newfound component of the ubiquitin-proteasome pathway in regulation of plant height and organ size in maize. These discoveries highlight potential targets for the genetic improvement of maize in the future.

A comprehensive atlas of endogenous peptides in maize

In this study, we present a comprehensive peptidomic atlas of 13 maize tissues, covering both vegetative and reproductive phases. Using a three-frame translation of canonical coding sequences, we identified 6100 nonredundant endogenous peptides, significantly expanding the known plant peptide repertoire. By integrating peptidomic coexpression profiles with previously reported proteomic profiles, we found that the peptide abundance did not consistently correlate with the abundance of their source proteins, suggesting the presence of complex regulatory mechanisms. This integrated peptidomic and proteomic map can serve as a valuable resource for exploring the functional roles of endogenous peptides in maize development and facilitates the investigation of the functional relationship among genes, peptides, and proteins across various biological contexts.

Single-cell proteomics differentiates Arabidopsis root cell types

Single-cell proteomics (SCP) is an emerging approach to resolve cellular heterogeneity within complex tissues of multi-cellular organisms. Here, we demonstrate the feasibility of SCP on plant samples using the model plant Arabidopsis thaliana. Specifically, we focused on examining isolated single cells from the cortex and endodermis, which are two adjacent root cell types derived from a common stem cell lineage. From 756 root cells, we identified 3763 proteins and 1118 proteins/cell. Ultimately, we focus on 3217 proteins quantified following stringent filtering. Of these, we identified 596 proteins whose expression is enriched in either the cortex or endodermis and are able to differentiate these closely related plant cell types. Collectivity, this study demonstrates that SCP can resolve neighboring cell types with distinct functions, thereby facilitating the identification of biomarkers and candidate proteins to enable functional genomics.

Transcriptomic Characterization of Genes Harboring Markers Linked to Maize Yield

BACKGROUND

It is currently believed that breeding priorities, including maize breeding, should focus on introducing varieties with greater utility value, specifically higher yields, into production. Global modern maize breeding relies on various molecular genetics techniques. Using the above mentioned technologies, we can identify regions of the genome that are associated with various phenotypic traits, including yield, which is of fundamental importance for understanding and manipulating these regions.

OBJECTIVES

The aim of the study was to analyze the expression of candidate genes associated with maize yield. To better understand the function of the analyzed genes in increasing maize yield, their expression in different organs and tissues was also assessed using publicly available transcriptome data.

METHODS

RT-qPCR analyses were performed using iTaq Universal SYBR Green Supermix (Bio-Rad, Hercules, CA, USA) and CFX96 Touch Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA). Each of the performed RT-qPCR experiments consisted of three biological replicates and three technical replicates, the results of which were averaged.

RESULTS

The research results allowed us to select three out of six candidate genes (cinnamoyl-CoA reductase 1-CCR1, aspartate aminotransferase-AAT and sucrose transporter 1-SUT1), which can significantly affect grain yield in maize. Not only our studies but also literature reports clearly indicate the participation of CCR1, AAT and SUT1 in the formation of yield. Identified molecular markers located within these genes can be used in breeding programs to select high yielding maize genotypes.

Data-driven identification of environmental variables influencing phenotypic plasticity to facilitate breeding for future climates

Phenotypic plasticity describes a genotype's ability to produce different phenotypes in response to different environments. Breeding crops that exhibit appropriate levels of plasticity for future climates will be crucial to meeting global demand, but knowledge of the critical environmental factors is limited to a handful of well-studied major crops. Using 727 maize (Zea mays L.) hybrids phenotyped for grain yield in 45 environments, we investigated the ability of a genetic algorithm and two other methods to identify environmental determinants of grain yield from a large set of candidate environmental variables constructed using minimal assumptions. The genetic algorithm identified pre- and postanthesis maximum temperature, mid-season solar radiation, and whole season net evapotranspiration as the four most important variables from a candidate set of 9150. Importantly, these four variables are supported by previous literature. After calculating reaction norms for each environmental variable, candidate genes were identified and gene annotations investigated to demonstrate how this method can generate insights into phenotypic plasticity. The genetic algorithm successfully identified known environmental determinants of hybrid maize grain yield. This demonstrates that the methodology could be applied to other less well-studied phenotypes and crops to improve understanding of phenotypic plasticity and facilitate breeding crops for future climates.

Genetic architecture of kernel-related traits in sweet and waxy maize revealed by genome-wide association analysis

Introduction

Maize (Zea mays L.) is one of the most important crops worldwide, the kernel size-related traits are the major components of maize grain yield.

Methods

To dissect the genetic architecture of four kernel-related traits of 100-kernel weight, kernel length, kernel width, and kernel diameter, a genome-wide association study (GWAS) was conducted in the waxy and sweet maize panel comprising of 447 maize inbred lines re-sequenced at the 5× coverage depth. GWAS analysis was carried out with the mixed linear model using 1,684,029 high-quality SNP markers.

Results

In total, 49 SNPs significantly associated with the four kernel-related traits were identified, including 46 SNPs on chromosome 3, two SNPs on chromosome 4, and one SNP on chromosome 7. Haplotype regression analysis identified 338 haplotypes that significantly affected these four kernel-related traits. Genomic selection (GS) results revealed that a set of 10,000 SNPs and a training population size of 30% are sufficient for the application of GS in waxy and sweet maize breeding for kernel weight and kernel size. Forty candidate genes associated with the four kernel-related traits were identified, including both Zm00001d000707 and Zm00001d044139 expressed in the kernel development tissues and stages with unknown functions.

Discussion

These significant SNPs and important haplotypes provide valuable information for developing functional markers for the implementation of marker-assisted selection in breeding. The molecular mechanism of Zm00001d000707 and Zm00001d044139 regulating these kernel-related traits needs to be investigated further.

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.

Genome-Wide Identification and Analysis of Maize DnaJ Family Genes in Response to Salt, Heat, and Cold at the Seedling Stage

DnaJ proteins, also known as HSP40s, play a key role in plant growth and development, and response to environmental stress. However, little comprehensive research has been conducted on the DnaJ gene family in maize. Here, we identify 91 ZmDnaJ genes from maize, which are likely distributed in the chloroplast, nucleus, and cytoplasm. Our analysis revealed that ZmDnaJs were classified into three types, with conserved protein motifs and gene structures within the same type, particularly among members of the same subfamily. Gene duplication events have likely contributed to the expansion of the ZmDnaJ family in maize. Analysis of cis-regulatory elements in ZmDnaJ promoters suggested involvement in stress responses, growth and development, and phytohormone sensitivity in maize. Specifically, four cis-acting regulatory elements associated with stress responses and phytohormone regulation indicated a role in adaptation. RNA-seq analysis showed constitutive expression of most ZmDnaJ genes, some specifically in pollen and endosperm. More importantly, certain genes also responded to salt, heat, and cold stresses, indicating potential interaction between stress regulatory networks. Furthermore, early responses to heat stress varied among five inbred lines, with upregulation of almost tested ZmDnaJ genes in B73 and B104 after 6 h, and fewer genes upregulated in QB1314, MD108, and Zheng58. After 72 h, most ZmDnaJ genes in the heat-sensitive inbred lines (B73 and B104) returned to normal levels, while many genes, including ZmDnaJ55, 79, 88, 90, and 91, remained upregulated in the heat-tolerant inbred lines (QB1314, MD108, and Zheng58) suggesting a synergistic function for prolonged protection against heat stress. In conclusion, our study provides a comprehensive analysis of the ZmDnaJ family in maize and demonstrates a correlation between heat stress tolerance and the regulation of gene expression within this family. These offer a theoretical basis for future functional validation of these genes.

COI1 F-box proteins regulate DELLA protein levels, growth, and photosynthetic efficiency in maize

The F-box protein Coronatine Insensitive (COI) is a receptor for the jasmonic acid signaling pathway in plants. To investigate the functions of the 6 maize (Zea mays) COI proteins (COI1a, COI1b, COI1c, COI1d, COI2a, and COI2b), we generated single, double, and quadruple loss-of-function mutants. The pollen of the coi2a coi2b double mutant was inviable. The coi1 quadruple mutant (coi1-4x) exhibited shorter internodes, decreased photosynthesis, leaf discoloration, microelement deficiencies, and accumulation of DWARF8 and/or DWARF9, 2 DELLA family proteins that repress the gibberellic acid (GA) signaling pathway. Coexpression of COI and DELLA in Nicotiana benthamiana showed that the COI proteins trigger proteasome-dependent DELLA degradation. Many genes that are downregulated in the coi1-4x mutant are GA-inducible. In addition, most of the proteins encoded by the downregulated genes are predicted to be bundle sheath- or mesophyll-enriched, including those encoding C4-specific photosynthetic enzymes. Heterologous expression of maize Coi genes in N. benthamiana showed that COI2a is nucleus-localized and interacts with maize jasmonate zinc-finger inflorescence meristem domain (JAZ) proteins, the canonical COI repressor partners. However, maize COI1a and COI1c showed only partial nuclear localization and reduced binding efficiency to the tested JAZ proteins. Together, these results show the divergent functions of the 6 COI proteins in regulating maize growth and defense pathways.

Genome-wide characterization and expression profiling of FARL (FHY3/FAR1) family genes in Zea mays

A significant role of the plant is played by the transcription factor FARL, which is light signal transduction as well as plant growth and development. Despite being transposases, FARL has developed a variety of dominant biological actions in evolution and speciation. On the other hand, little is known about the Zea mays FARL protein family. This study identifies and characterizes fifteen ZmFARL genes genome-wide, and RNA sequencing data was used to profile their expression. 105 FARL proteins from five plant species were classified into five groups based on sequence alignment and phylogeny. The ZmFARL genes' exon-intron and motif distribution were conserved based on their evolutionary group. The fifteen ZmFARL genes were distributed over seven of the ten Z. mays chromosomes, although no duplication was discovered. Cis-element analysis reveals that ZmFARL genes play a variety of activities, including tissue-specific, stress- and hormone-responsive expressions. Furthermore, the results of the RNA sequencing used to profile expression showed that the genes ZmFARL2 and ZmFARL5 were much more expressed than other genes in various tissues, particularly in leaf characteristics. The identification of likely genes involved in cellular activity in Z. mays and related species will be aided by the characterization of the FARL genes.

Teosinte Pollen Drive guides maize diversification and domestication by RNAi

Selfish genetic elements contribute to hybrid incompatibility and bias or 'drive' their own transmission1,2. Chromosomal drive typically functions in asymmetric female meiosis, whereas gene drive is normally post-meiotic and typically found in males. Here, using single-molecule and single-pollen genome sequencing, we describe Teosinte Pollen Drive, an instance of gene drive in hybrids between maize (Zea mays ssp. mays) and teosinte mexicana (Z. mays ssp. mexicana) that depends on RNA interference (RNAi). 22-nucleotide small RNAs from a non-coding RNA hairpin in mexicana depend on Dicer-like 2 (Dcl2) and target Teosinte Drive Responder 1 (Tdr1), which encodes a lipase required for pollen viability. Dcl2, Tdr1 and the hairpin are in tight pseudolinkage on chromosome 5, but only when transmitted through the male. Introgression of mexicana into early cultivated maize is thought to have been critical to its geographical dispersal throughout the Americas3, and a tightly linked inversion in mexicana spans a major domestication sweep in modern maize4. A survey of maize traditional varieties and sympatric populations of teosinte mexicana reveals correlated patterns of admixture among unlinked genes required for RNAi on at least four chromosomes that are also subject to gene drive in pollen from synthetic hybrids. Teosinte Pollen Drive probably had a major role in maize domestication and diversification, and offers an explanation for the widespread abundance of 'self' small RNAs in the germ lines of plants and animals.

Insights into maydis leaf blight resistance in maize: a comprehensive genome-wide association study in sub-tropics of India

BACKGROUND

In the face of contemporary climatic vulnerabilities and escalating global temperatures, the prevalence of maydis leaf blight (MLB) poses a potential threat to maize production. This study endeavours to discern marker-trait associations and elucidate the candidate genes that underlie resistance to MLB in maize by employing a diverse panel comprising 336 lines. The panel was screening for MLB across four environments, employing standard artificial inoculation techniques. Genome-wide association studies (GWAS) and haplotype analysis were conducted utilizing a total of 128,490 SNPs obtained from genotyping-by-sequencing (GBS).

RESULTS

GWAS identified 26 highly significant SNPs associated with MLB resistance, among the markers examined. Seven of these SNPs, reported in novel chromosomal bins (9.06, 5.01, 9.01, 7.04, 4.06, 1.04, and 6.05) were associated with genes: bzip23, NAGS1, CDPK7, aspartic proteinase NEP-2, VQ4, and Wun1, which were characterized for their roles in diminishing fungal activity, fortifying defence mechanisms against necrotrophic pathogens, modulating phyto-hormone signalling, and orchestrating oxidative burst responses. Gene mining approach identified 22 potential candidate genes associated with SNPs due to their functional relevance to resistance against necrotrophic pathogens. Notably, bin 8.06, which hosts five SNPs, showed a connection to defense-regulating genes against MLB, indicating the potential formation of a functional gene cluster that triggers a cascade of reactions against MLB. In silico studies revealed gene expression levels exceeding ten fragments per kilobase million (FPKM) for most genes and demonstrated coexpression among all candidate genes in the coexpression network. Haplotype regression analysis revealed the association of 13 common significant haplotypes at Bonferroni ≤ 0.05. The phenotypic variance explained by these significant haplotypes ranged from low to moderate, suggesting a breeding strategy that combines multiple resistance alleles to enhance resistance to MLB. Additionally, one particular haplotype block (Hap_8.3) was found to consist of two SNPs (S8_152715134, S8_152460815) identified in GWAS with 9.45% variation explained (PVE).

CONCLUSION

The identified SNPs/ haplotypes associated with the trait of interest contribute to the enrichment of allelic diversity and hold direct applicability in Genomics Assisted Breeding for enhancing MLB resistance in maize.

Engineering the future cereal crops with big biological data: toward intelligence-driven breeding by design

How to feed 10 billion human populations is one of the challenges that need to be addressed in the following decades, especially under an unpredicted climate change. Crop breeding, initiating from the phenotype-based selection by local farmers and developing into current biotechnology-based breeding, has played a critical role in securing the global food supply. However, regarding the changing environment and ever-increasing human population, can we breed outstanding crop varieties fast enough to achieve high productivity, good quality, and widespread adaptability? This review outlines the recent achievements in understanding cereal crop breeding, including the current knowledge about crop agronomic traits, newly developed techniques, crop big biological data research, and the possibility of integrating them for intelligence-driven breeding by design, which ushers in a new era of crop breeding practice and shapes the novel architecture of future crops. This review focuses on the major cereal crops, including rice, maize, and wheat, to explain how intelligence-driven breeding by design is becoming a reality.

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

A foundational large language model for edible plant genomes

Significant progress has been made in the field of plant genomics, as demonstrated by the increased use of high-throughput methodologies that enable the characterization of multiple genome-wide molecular phenotypes. These findings have provided valuable insights into plant traits and their underlying genetic mechanisms, particularly in model plant species. Nonetheless, effectively leveraging them to make accurate predictions represents a critical step in crop genomic improvement. We present AgroNT, a foundational large language model trained on genomes from 48 plant species with a predominant focus on crop species. We show that AgroNT can obtain state-of-the-art predictions for regulatory annotations, promoter/terminator strength, tissue-specific gene expression, and prioritize functional variants. We conduct a large-scale in silico saturation mutagenesis analysis on cassava to evaluate the regulatory impact of over 10 million mutations and provide their predicted effects as a resource for variant characterization. Finally, we propose the use of the diverse datasets compiled here as the Plants Genomic Benchmark (PGB), providing a comprehensive benchmark for deep learning-based methods in plant genomic research. The pre-trained AgroNT model is publicly available on HuggingFace at https://huggingface.co/InstaDeepAI/agro-nucleotide-transformer-1b  for future research purposes.

MTP family analysis and association study reveal the role of ZmMTP11 in lead (Pb) accumulation

The metal tolerance protein (MTP) gene family plays an essential role in the transport of heavy metals, however the function of the MTP family in transporting lead (Pb) was still unclear in plants. In this study, we identified and characterized 12 ZmMTPs in the whole genome of maize. These ZmMTP genes were divided into three subfamilies in evolution, namely Zn-CDF, Zn/Fe-CDF, Mn-CDF subfamilies, which showed diverse expression patterns in different tissues of maize. Using gene-based association analyses, we identified a Pb accumulation-related MTP member in maize, ZmMTP11, which was located in plasma membrane and had the potential of transporting Pb ion. Under the Pb treatment, ZmMTP11 showed a generally decreased expression relative to the normal conditions. Heterologous expressions of ZmMTP11 in yeast, Arabidopsis, and rice demonstrated that ZmMTP11 enhanced Pb accumulation in the cells without affecting yeast and plant growth under Pb stress. Remarkably, the increased Pb concentration in the plant roots did not cause changes in Pb content in the shoots. Our study provides new insights into the genetic improvement of heavy metal tolerance in plants and contributes to bioremediation of Pb-contaminant soils.

Comprehensive identification and expression analysis of FAR1/FHY3 genes under drought stress in maize (Zea mays L.)

Background

FAR1/FHY3 transcription factors are derived from transposase, which play important roles in light signal transduction, growth and development, and response to stress by regulating downstream gene expression. Although many FAR1/FHY3 members have been identified in various species, the FAR1/FHY3 genes in maize are not well characterized and their function in drought are unknown.

Method

The FAR1/FHY3 family in the maize genome was identified using PlantTFDB, Pfam, Smart, and NCBI-CDD websites. In order to investigate the evolution and functions of FAR1 genes in maize, the information of protein sequences, chromosome localization, subcellular localization, conserved motifs, evolutionary relationships and tissue expression patterns were analyzed by bioinformatics, and the expression patterns under drought stress were detected by quantitative real-time polymerase chain reaction (qRT-PCR).

Results

A total of 24 ZmFAR members in maize genome, which can be divided into five subfamilies, with large differences in protein and gene structures among subfamilies. The promoter regions of ZmFARs contain abundant abiotic stress-responsive and hormone-respovensive cis-elements. Among them, drought-responsive cis-elements are quite abundant. ZmFARs were expressed in all tissues detected, but the expression level varies widely. The expression of ZmFARs were mostly down-regulated in primary roots, seminal roots, lateral roots, and mesocotyls under water deficit. Most ZmFARs were down-regulated in root after PEG-simulated drought stress.

Conclusions

We performed a genome-wide and systematic identification of FAR1/FHY3 genes in maize. And most ZmFARs were down-regulated in root after drought stress. These results indicate that FAR1/FHY3 transcription factors have important roles in drought stress response, which can lay a foundation for further analysis of the functions of ZmFARs in response to drought stress.

ZmNF-YA1 Contributes to Maize Thermotolerance by Regulating Heat Shock Response

Zea mays (maize) is a staple food, feed, and industrial crop. Heat stress is one of the major stresses affecting maize production and is usually accompanied by other stresses, such as drought. Our previous study identified a heterotrimer complex, ZmNF-YA1-YB16-YC17, in maize. ZmNF-YA1 and ZmNF-YB16 were positive regulators of the drought stress response and were involved in maize root development. In this study, we investigated whether ZmNF-YA1 confers heat stress tolerance in maize. The nf-ya1 mutant and overexpression lines were used to test the role of ZmNF-YA1 in maize thermotolerance. The nf-ya1 mutant was more temperature-sensitive than the wild-type (WT), while the ZmNF-YA1 overexpression lines showed a thermotolerant phenotype. Higher malondialdehyde (MDA) content and reactive oxygen species (ROS) accumulation were observed in the mutant, followed by WT and overexpression lines after heat stress treatment, while an opposite trend was observed for chlorophyll content. RNA-seq was used to analyze transcriptome changes in nf-ya1 and its wild-type control W22 in response to heat stress. Based on their expression profiles, the heat stress response-related differentially expressed genes (DEGs) in nf-ya1 compared to WT were grouped into seven clusters via k-means clustering. Gene Ontology (GO) enrichment analysis of the DEGs in different clades was performed to elucidate the roles of ZmNF-YA1-mediated transcriptional regulation and their contribution to maize thermotolerance. The loss function of ZmNF-YA1 led to the failure induction of DEGs in GO terms of protein refolding, protein stabilization, and GO terms for various stress responses. Thus, the contribution of ZmNF-YA1 to protein stabilization, refolding, and regulation of abscisic acid (ABA), ROS, and heat/temperature signaling may be the major reason why ZmNF-YA1 overexpression enhanced heat tolerance, and the mutant showed a heat-sensitive phenotype.

Recent advances in exploring transcriptional regulatory landscape of crops

Crop breeding entails developing and selecting plant varieties with improved agronomic traits. Modern molecular techniques, such as genome editing, enable more efficient manipulation of plant phenotype by altering the expression of particular regulatory or functional genes. Hence, it is essential to thoroughly comprehend the transcriptional regulatory mechanisms that underpin these traits. In the multi-omics era, a large amount of omics data has been generated for diverse crop species, including genomics, epigenomics, transcriptomics, proteomics, and single-cell omics. The abundant data resources and the emergence of advanced computational tools offer unprecedented opportunities for obtaining a holistic view and profound understanding of the regulatory processes linked to desirable traits. This review focuses on integrated network approaches that utilize multi-omics data to investigate gene expression regulation. Various types of regulatory networks and their inference methods are discussed, focusing on recent advancements in crop plants. The integration of multi-omics data has been proven to be crucial for the construction of high-confidence regulatory networks. With the refinement of these methodologies, they will significantly enhance crop breeding efforts and contribute to global food security.

Dual-Model GWAS Analysis and Genomic Selection of Maize Flowering Time-Related Traits

An appropriate flowering period is an important selection criterion in maize breeding. It plays a crucial role in the ecological adaptability of maize varieties. To explore the genetic basis of flowering time, GWAS and GS analyses were conducted using an associating panel consisting of 379 multi-parent DH lines. The DH population was phenotyped for days to tasseling (DTT), days to pollen-shedding (DTP), and days to silking (DTS) in different environments. The heritability was 82.75%, 86.09%, and 85.26% for DTT, DTP, and DTS, respectively. The GWAS analysis with the FarmCPU model identified 10 single-nucleotide polymorphisms (SNPs) distributed on chromosomes 3, 8, 9, and 10 that were significantly associated with flowering time-related traits. The GWAS analysis with the BLINK model identified seven SNPs distributed on chromosomes 1, 3, 8, 9, and 10 that were significantly associated with flowering time-related traits. Three SNPs 3_198946071, 9_146646966, and 9_152140631 showed a pleiotropic effect, indicating a significant genetic correlation between DTT, DTP, and DTS. A total of 24 candidate genes were detected. A relatively high prediction accuracy was achieved with 100 significantly associated SNPs detected from GWAS, and the optimal training population size was 70%. This study provides a better understanding of the genetic architecture of flowering time-related traits and provides an optimal strategy for GS.

Factors specifying sex determination in maize

Plant architecture is an important feature for agronomic performance in crops. In maize, which is a monoecious plant, separation of floral organs to produce specific gametes has been studied from different perspectives including genetic, biochemical and physiological. Maize mutants affected in floral organ development have been key to identifying genes, hormones and other factors like miRNAs important for sex determination. In this review, we describe floral organ formation in maize, representative mutants and genes identified with a function in establishing sexual identity either classified as feminizing or masculinizing, and its relationship with hormones associated with sexual organ identity as jasmonic acid, brassinosteroid and gibberellin. Finally, we discuss the challenges and scopes of future research in maize sex determination.

Perfluorobutanoic acid triggers metabolic and transcriptional reprogramming in wheat seedlings

The environmental risks of fluorinated alternatives are of great concern with the phasing out of perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate. Here, multi-omics (i.e., metabolomics and transcriptomics) coupled with physiological and biochemical analyses were employed to investigate the stress responses of wheat seedings (Triticum aestivum L.) to perfluorobutanoic acid (PFBA), one of the short-chain per- and polyfluoroalkyl substances (PFAS) and PFOA alternatives, at environmentally relevant concentrations (0.1-100 ng/g). After 28 days of soil exposure, PFBA boosted the generation of OH and O2- in wheat seedlings, resulting in lipid peroxidation, protein perturbation and impaired photosynthesis. Non-enzymatic antioxidant defense systems (e.g., glutathione, phenolics, and vitamin C) and enzymatic antioxidant copper/zinc superoxide dismutase were strikingly activated (p < 0.05). PFBA-triggered oxidative stress induced metabolic and transcriptional reprogramming, including carbon and nitrogen metabolisms, lipid metabolisms, immune responses, signal transduction processes, and antioxidant defense-related pathways. Down-regulation of genes related to plant-pathogen interaction suggested suppression of the immune-response, offering a novel understanding on the production of reactive oxygen species in plants under the exposure to PFAS. The identified MAPK signaling pathway illuminated a novel signal transduction mechanism in plant cells in response to PFAS. These findings provide comprehensive understandings on the phytotoxicity of PFBA to wheat seedlings and new insights into the impacts of PFAS on plants.

Establishing an optimized ATAC-seq protocol for the maize

The advent of next-generation sequencing in crop improvement offers unprecedented insights into the chromatin landscape closely linked to gene activity governing key traits in plant development and adaptation. Particularly in maize, its dynamic chromatin structure is found to collaborate with massive transcriptional variations across tissues and developmental stages, implying intricate regulatory mechanisms, which highlights the importance of integrating chromatin information into breeding strategies for precise gene controls. The depiction of maize chromatin architecture using Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) provides great opportunities to investigate cis-regulatory elements, which is crucial for crop improvement. In this context, we developed an easy-to-implement ATAC-seq protocol for maize with fewer nuclei and simple equipment. We demonstrate a streamlined ATAC-seq protocol with four key steps for maize in which nuclei purification can be achieved without cell sorting and using only a standard bench-top centrifuge. Our protocol, coupled with the bioinformatic analysis, including validation by read length periodicity, key metrics, and correlation with transcript abundance, provides a precise and efficient assessment of the maize chromatin landscape. Beyond its application to maize, our testing design holds the potential to be applied to other crops or other tissues, especially for those with limited size and amount, establishing a robust foundation for chromatin structure studies in diverse crop species.

ZmPILS6 is an auxin efflux carrier required for maize root morphogenesis

Plant root systems play a pivotal role in plant physiology and exhibit diverse phenotypic traits. Understanding the genetic mechanisms governing root growth and development in model plants like maize is crucial for enhancing crop resilience to drought and nutrient limitations. This study focused on identifying and characterizing ZmPILS6, an annotated auxin efflux carrier, as a key regulator of various crown root traits in maize. ZmPILS6-modified roots displayed reduced network area and suppressed lateral root formation, which are desirable traits for the "steep, cheap, and deep" ideotype. The research revealed that ZmPILS6 localizes to the endoplasmic reticulum and plays a vital role in controlling the spatial distribution of indole-3-acetic acid (IAA or "auxin") in primary roots. The study also demonstrated that ZmPILS6 can actively efflux IAA when expressed in yeast. Furthermore, the loss of ZmPILS6 resulted in significant proteome remodeling in maize roots, particularly affecting hormone signaling pathways. To identify potential interacting partners of ZmPILS6, a weighted gene coexpression analysis was performed. Altogether, this research contributes to the growing knowledge of essential genetic determinants governing maize root morphogenesis, which is crucial for guiding agricultural improvement strategies.

Dissection of Common Rust Resistance in Tropical Maize Multiparent Population through GWAS and Linkage Studies

Common rust (CR), caused by Puccina sorghi, is a major foliar disease in maize that leads to quality deterioration and yield losses. To dissect the genetic architecture of CR resistance in maize, this study utilized the susceptible temperate inbred line Ye107 as the male parent crossed with three resistant tropical maize inbred lines (CML312, D39, and Y32) to generate 627 F7 recombinant inbred lines (RILs), with the aim of identifying maize disease-resistant loci and candidate genes for common rust. Phenotypic data showed good segregation between resistance and susceptibility, with varying degrees of resistance observed across different subpopulations. Significant genotype effects and genotype × environment interactions were observed, with heritability ranging from 85.7% to 92.2%. Linkage and genome-wide association analyses across the three environments identified 20 QTLs and 62 significant SNPs. Among these, seven major QTLs explained 66% of the phenotypic variance. Comparison with six SNPs repeatedly identified across different environments revealed overlap between qRUST3-3 and Snp-203,116,453, and Snp-204,202,469. Haplotype analysis indicated two different haplotypes for CR resistance for both the SNPs. Based on LD decay plots, three co-located candidate genes, Zm00001d043536, Zm00001d043566, and Zm00001d043569, were identified within 20 kb upstream and downstream of these two SNPs. Zm00001d043536 regulates hormone regulation, Zm00001d043566 controls stomatal opening and closure, related to trichome, and Zm00001d043569 is associated with plant disease immune responses. Additionally, we performed candidate gene screening for five additional SNPs that were repeatedly detected across different environments, resulting in the identification of five candidate genes. These findings contribute to the development of genetic resources for common rust resistance in maize breeding programs.

Multiomics unravels potential molecular switches in the C3 to CAM transition of Mesembryanthemum crystallinum

Mesembryanthemum crystallinum (common ice plant), a facultative CAM plant, shifts from C3 to CAM photosynthesis under salt stress, enhancing water use efficiency. Here we used transcriptomics, proteomics, and targeted metabolomics to profile molecular changes during the diel cycle of C3 to CAM transition. The results confirmed expected changes associated with CAM photosynthesis, starch biosynthesis and degradation, and glycolysis/gluconeogenesis. Importantly, they yielded new discoveries: 1) Transcripts displayed greater circadian regulation than proteins. 2) Oxidative phosphorylation and inositol methylation may play important roles in initiating the transition. 3) V-type H+-ATPases showed consistent transcriptional regulation, aiding in vacuolar malate uptake. 4) A protein phosphatase 2C, a major component in the ABA signaling pathway, may trigger the C3 to CAM transition. Our work highlights the potential molecular switches in the C3 to CAM transition, including the potential role of ABA signaling. SIGNIFICANCE: The common ice plant is a model facultative CAM plant, and under stress conditions it can shift from C3 to CAM photosynthesis within a three-day period. However, knowledge about the molecular changes during the transition and the molecular switches enabling the transition is lacking. Multi-omic analyses not only revealed the molecular changes during the transition, but also highlighted the importance of ABA signaling, inositol methylation, V-type H+-ATPase in initiating the shift. The findings may explain physiological changes and nocturnal stomatal opening, and inform future synthetic biology effort in improving crop water use efficiency and stress resilience.

Genome-wide expression quantitative trait locus analysis reveals silk-preferential gene regulatory network in maize

Silk of maize (Zea mays L.) contains diverse metabolites with complicated structures and functions, making it a great challenge to explore the mechanisms of metabolic regulation. Genome-wide identification of silk-preferential genes and investigation of their expression regulation provide an opportunity to reveal the regulatory networks of metabolism. Here, we applied the expression quantitative trait locus (eQTL) mapping on a maize natural population to explore the regulation of gene expression in unpollinated silk of maize. We obtained 3,985 silk-preferential genes that were specifically or preferentially expressed in silk using our population. Silk-preferential genes showed more obvious expression variations compared with broadly expressed genes that were ubiquitously expressed in most tissues. We found that trans-eQTL regulation played a more important role for silk-preferential genes compared to the broadly expressed genes. The relationship between 38 transcription factors and 85 target genes, including silk-preferential genes, were detected. Finally, we constructed a transcriptional regulatory network around the silk-preferential gene Bx10, which was proposed to be associated with response to abiotic stress and biotic stress. Taken together, this study deepened our understanding of transcriptome variation in maize silk and the expression regulation of silk-preferential genes, enhancing the investigation of regulatory networks on metabolic pathways.

Should I stay or should I go? Trafficking of plant extra-nuclear transcription factors

At the heart of all biological processes lies the control of nuclear gene expression, which is primarily achieved through the action of transcription factors (TFs) that generally contain a nuclear localization signal (NLS) to facilitate their transport into the nucleus. However, some TFs reside in the cytoplasm in a transcriptionally inactive state and only enter the nucleus in response to specific signals, which in plants include biotic or abiotic stresses. These extra-nuclear TFs can be found in the cytosol or associated with various membrane systems, including the endoplasmic reticulum and plasma membrane. They may be integral proteins with transmembrane domains or associate peripherally with the lipid bilayer via acylation or membrane-binding domains. Although over 30 plant TFs, most of them involved in stress responses, have been experimentally shown to reside outside the nucleus, computational predictions suggest that this number is much larger. Understanding how extra-nuclear TFs are trafficked into the nucleus is essential for reconstructing transcriptional regulatory networks that govern major cellular pathways in response to biotic and abiotic signals. Here, we provide a perspective on what is known on plant extranuclear-nuclear TF retention, nuclear trafficking, and the post-translational modifications that ultimately enable them to regulate gene expression upon entering the nucleus.

Genome-wide investigation of the LARP gene family: focus on functional identification and transcriptome profiling of ZmLARP6c1 in maize pollen

BACKGROUND

The La-related proteins (LARPs) are a superfamily of RNA-binding proteins associated with regulation of gene expression. Evidence points to an important role for post-transcriptional control of gene expression in germinating pollen tubes, which could be aided by RNA-binding proteins.

RESULTS

In this study, a genome-wide investigation of the LARP proteins in eight plant species was performed. The LARP proteins were classified into three families based on a phylogenetic analysis. The gene structure, conserved motifs, cis-acting elements in the promoter, and gene expression profiles were investigated to provide a comprehensive overview of the evolutionary history and potential functions of ZmLARP genes in maize. Moreover, ZmLARP6c1 was specifically expressed in pollen and ZmLARP6c1 was localized to the nucleus and cytoplasm in maize protoplasts. Overexpression of ZmLARP6c1 enhanced the percentage pollen germination compared with that of wild-type pollen. In addition, transcriptome profiling analysis revealed that differentially expressed genes included PABP homologous genes and genes involved in jasmonic acid and abscisic acid biosynthesis, metabolism, signaling pathways and response in a Zmlarp6c1::Ds mutant and ZmLARP6c1-overexpression line compared with the corresponding wild type.

CONCLUSIONS

The findings provide a basis for further evolutionary and functional analyses, and provide insight into the critical regulatory function of ZmLARP6c1 in maize pollen germination.

Sugar transporters spatially organize microbiota colonization along the longitudinal root axis of Arabidopsis

Plant roots are functionally heterogeneous in cellular architecture, transcriptome profile, metabolic state, and microbial immunity. We hypothesized that axial differentiation may also impact spatial colonization by root microbiota along the root axis. We developed two growth systems, ArtSoil and CD-Rhizotron, to grow and then dissect Arabidopsis thaliana roots into three segments. We demonstrate that distinct endospheric and rhizosphere bacterial communities colonize the segments, supporting the hypothesis of microbiota differentiation along the axis. Root metabolite profiling of each segment reveals differential metabolite enrichment and specificity. Bioinformatic analyses and GUS histochemistry indicate microbe-induced accumulation of SWEET2, 4, and 12 sugar uniporters. Profiling of root segments from sweet mutants shows altered spatial metabolic profiles and reorganization of endospheric root microbiota. This work reveals the interdependency between root metabolites and microbial colonization and the contribution of SWEETs to spatial diversity and stability of microbial ecosystem.

Maize actin depolymerizing factor 1 (ZmADF1) negatively regulates pollen development

The normal development of pollen grains and the completion of double fertilization in embryos are crucial for both the sexual reproduction of angiosperms and grain production. Actin depolymerizing factor (ADF) regulates growth, development, and responses to biotic and abiotic stress by binding to actin in plants. In this study, the function of the ZmADF1 gene was validated through bioinformatic analysis, subcellular localization, overexpression in maize and Arabidopsis, and knockout via CRISPR/Cas9. The amino acid sequence of ZmADF1 exhibited high conservation and a similar tertiary structure to that of ADF homologs. Subcellular localization analysis revealed that ZmADF1 is localized mainly to the nucleus and cytoplasm. The ZmADF1 gene was specifically expressed in maize pollen, and overexpression of the ZmADF1 gene decreased the number of pollen grains in the anthers of transgenic Arabidopsis plants. The germination rate of pollen and the empty seed shell rate in the fruit pods of the overexpressing plants were significantly greater than those in the wild-type (WT) plants. In maize, the pollen viability of the knockout lines was significantly greater than that of both the WT and the overexpressing lines. Our results confirmed that the ZmADF1 gene was specifically expressed in pollen and negatively regulated pollen quantity, vigor, germination rate, and seed setting rate. This study provides insights into ADF gene function and possible pathways for improving high-yield maize breeding.

The maize PLASTID TERMINAL OXIDASE (PTOX) locus controls the carotenoid content of kernels

Carotenoids perform a broad range of important functions in humans; therefore, carotenoid biofortification of maize (Zea mays L.), one of the most highly produced cereal crops worldwide, would have a global impact on human health. PLASTID TERMINAL OXIDASE (PTOX) genes play an important role in carotenoid metabolism; however, the possible function of PTOX in carotenoid biosynthesis in maize has not yet been explored. In this study, we characterized the maize PTOX locus by forward- and reverse-genetic analyses. While most higher plant species possess a single copy of the PTOX gene, maize carries two tandemly duplicated copies. Characterization of mutants revealed that disruption of either copy resulted in a carotenoid-deficient phenotype. We identified mutations in the PTOX genes as being causal of the classic maize mutant, albescent1. Remarkably, overexpression of ZmPTOX1 significantly improved the content of carotenoids, especially β-carotene (provitamin A), which was increased by ~threefold, in maize kernels. Overall, our study shows that maize PTOX locus plays an important role in carotenoid biosynthesis in maize kernels and suggests that fine-tuning the expression of this gene could improve the nutritional value of cereal grains.

Maize multi-omics reveal leaf water status controlling of differential transcriptomes, proteomes and hormones as mechanisms of age-dependent osmotic stress response in leaves

Drought-induced osmotic stress severely affects the growth and yield of maize. However, the mechanisms underlying the different responses of young and old maize leaves to osmotic stress remain unclear. To gain a systematic understanding of age-related stress responses, we compared osmotic-stress-induced changes in maize leaves of different ages using multi-omics approaches. After short-term osmotic stress, old leaves suffered more severe water deficits than young leaves. The adjustments of transcriptomes, proteomes, and hormones in response to osmotic stress were more dynamic in old leaves. Metabolic activities, stress signaling pathways, and hormones (especially abscisic acid) responded to osmotic stress in an age-dependent manner. We identified multiple functional clusters of genes and proteins with potential roles in stress adaptation. Old leaves significantly accumulated stress proteins such as dehydrin, aquaporin, and chaperones to cope with osmotic stress, accompanied by senescence-like cellular events, whereas young leaves exhibited an effective water conservation strategy mainly by hydrolyzing transitory starch and increasing proline production. The stress responses of individual leaves are primarily determined by their intracellular water status, resulting in differential transcriptomes, proteomes, and hormones. This study extends our understanding of the mechanisms underlying plant responses to osmotic stress.

The BulkECexplorer compiles endothelial bulk transcriptomes to predict functional versus leaky transcription

Transcriptomic data can be mined to understand the molecular activity of cell types. Yet, functional genes may remain undetected in RNA sequencing (RNA-seq) experiments for technical reasons, such as insufficient read depth or gene dropout. Conversely, RNA-seq experiments may detect lowly expressed mRNAs thought to be biologically irrelevant products of leaky transcription. To represent a cell type's functional transcriptome more accurately, we propose compiling many bulk RNA-seq datasets into a compendium and applying established classification models to predict whether detected transcripts are likely products of active or leaky transcription. Here, we present the BulkECexplorer (bulk RNA-seq endothelial cell explorer) compendium of 240 bulk RNA-seq datasets from five vascular endothelial cell subtypes. This resource reports transcript counts for genes of interest and predicts whether detected transcripts are likely the products of active or leaky gene expression. Beyond its usefulness for vascular biology research, this resource provides a blueprint for developing analogous tools for other cell types.

A review of artificial intelligence-assisted omics techniques in plant defense: current trends and future directions

Plants intricately deploy defense systems to counter diverse biotic and abiotic stresses. Omics technologies, spanning genomics, transcriptomics, proteomics, and metabolomics, have revolutionized the exploration of plant defense mechanisms, unraveling molecular intricacies in response to various stressors. However, the complexity and scale of omics data necessitate sophisticated analytical tools for meaningful insights. This review delves into the application of artificial intelligence algorithms, particularly machine learning and deep learning, as promising approaches for deciphering complex omics data in plant defense research. The overview encompasses key omics techniques and addresses the challenges and limitations inherent in current AI-assisted omics approaches. Moreover, it contemplates potential future directions in this dynamic field. In summary, AI-assisted omics techniques present a robust toolkit, enabling a profound understanding of the molecular foundations of plant defense and paving the way for more effective crop protection strategies amidst climate change and emerging diseases.

Evidence that variation in root anatomy contributes to local adaptation in Mexican native maize

Mexican native maize (Zea mays ssp. mays) is adapted to a wide range of climatic and edaphic conditions. Here, we focus specifically on the potential role of root anatomical variation in this adaptation. Given the investment required to characterize root anatomy, we present a machine-learning approach using environmental descriptors to project trait variation from a relatively small training panel onto a larger panel of genotyped and georeferenced Mexican maize accessions. The resulting models defined potential biologically relevant clines across a complex environment that we used subsequently for genotype-environment association. We found evidence of systematic variation in maize root anatomy across Mexico, notably a prevalence of trait combinations favoring a reduction in axial hydraulic conductance in varieties sourced from cooler, drier highland areas. We discuss our results in the context of previously described water-banking strategies and present candidate genes that are associated with both root anatomical and environmental variation. Our strategy is a refinement of standard environmental genome-wide association analysis that is applicable whenever a training set of georeferenced phenotypic data is available.

Genetic basis of maize stalk strength decoded via linkage and association mapping

Stalk lodging is a severe problem that limits maize production worldwide, although little attention has been given to its genetic basis. Here we measured rind penetrometer resistance (RPR), an effective index for stalk lodging, in a multi-parent population of 1948 recombinant inbred lines (RILs) and an association population of 508 inbred lines (AMP508). Linkage and association mapping identified 53 and 29 single quantitative trait loci (QTLs) and 50 and 19 pairs of epistatic interactions for RPR in the multi-parent population and AMP508 population, respectively. Phenotypic variation explained by all identified epistatic QTLs (up to ~5%) was much less than that explained by all single additive QTLs (up to ~33% in the multi-parent population and ~ 60% in the AMP508 population). Among all detected QTLs, only eight single QTLs explained >10% of phenotypic variation in single RIL populations. Alleles that increased RPR were enriched in tropical/subtropical (TST) groups from the AMP508 population. Based on genome-wide association studies in both populations, we identified 137 candidate genes affecting RPR, which were assigned to multiple biological processes, such as the biosynthesis of cell wall components. Sixty-six candidate genes were cross-validated by multiple methods or populations. Most importantly, 23 candidate genes were upregulated or downregulated in high-RPR lines relative to low-RPR lines, supporting the associations between candidate genes and RPR. These findings reveal the complex nature of the genetic basis underlying RPR and provide loci or candidate genes for developing elite varieties that are resistant to stalk lodging via molecular breeding.