An Expanded Maize Gene Expression Atlas based on RNA Sequencing and its Use to Explore Root Development

The genetic architecture of cell type-specific cis regulation in maize

Gene expression and complex phenotypes are determined by the activity of cis-regulatory elements. However, an understanding of how extant genetic variants affect cis regulation remains limited. Here, we investigated the consequences of cis-regulatory diversity using single-cell genomics of more than 0.7 million nuclei across 172 Zea mays (maize) inbreds. Our analyses pinpointed cis-regulatory elements distinct to domesticated maize and revealed how historical transposon activity has shaped the cis-regulatory landscape. Leveraging population genetics principles, we fine-mapped about 22,000 chromatin accessibility-associated genetic variants with widespread cell type-specific effects. Variants in TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR-binding sites were the most prevalent determinants of chromatin accessibility. Finally, integrating chromatin accessibility-associated variants, organismal trait variation, and population differentiation revealed how local adaptation has rewired regulatory networks in unique cellular contexts to alter maize flowering.

Functional characterization of a maize UDP-glucosyltransferase involved in the biosynthesis of flavonoid 7-O-glucosides and di-O-glucosides

Flavonoids are glycosylated in the final steps of their synthesis by UDP-dependent glycosyltransferase enzymes. We present the functional characterization of the first maize flavonoid O-glucosyltransferase enzyme from family 85, which exhibits properties not previously described. ZmUGT85W2 catalyzes the O-glucosylation of flavonols, flavones and flavanones, to form primarily 7-O-glucosides, but also flavonol O-glucoside positional isomers, flavones and flavonol di-O-glucosides. ZmUGT85W2 exhibited a differential kinetic behavior depending on the flavonoid acceptor, showing hyperbolic dependence for flavonols and sigmoidal response for flavanones and flavones. Structural and molecular docking analyses predicted conserved residues interacting with the sugar donor, with close contact with the 7-hydroxyl of the flavonoid acceptors, consistent with enzymatic activity results. In addition, ZmUGT85W2 is induced by UV-B radiation, and its expression is controlled by the B and PL1 transcription factors. Consistently, higher levels of flavone and flavonol O-glycosides are accumulated in leaves of plants exposed to solar UV-B compared to control plants, suggesting that ZmUGT85W2 is involved in the biosynthesis of these metabolites in maize leaves, contributing to UV-B tolerance. The activity of ZmUGT85W2, along with its elevated expression in silks and pericarps expressing the R2R3-MYB transcription factor P1, highlights its critical role in the accumulation of flavonoid O-glucosides in these tissues. Together, our findings reveal a key step in maize flavonoid O-glycosides biosynthesis, with the observed positive cooperative behaviors suggesting that ZmUGT85W2 plays a crucial role in finely regulating metabolic flux towards these compounds in planta.

Quantifying Nitrogen Uptake Rates of Maize Roots Using Stable Isotopes

Nitrogen is an essential element for plant growth and development; however, application of nitrogen (N)-based fertilizers comes with a high environmental cost. This includes the energy required for production, volatilization from fields, and runoff or leaching to waterways triggering algal blooms. As such, a key goal in plant breeding programs is to develop varieties that maintain yield while requiring less fertilization. Central to this goal is understanding how roots take up nitrogen and finding traits that represent improvements in the net uptake. Maize, one of the most widely produced crops in the world, has seminal, crown, and brace root types, each under independent developmental control. Recent evidence suggests that these independent developmental patterns may result in different nutrient uptake characteristics. As such, understanding the uptake dynamics of each root type under different environmental conditions is an essential aspect for the selection of new maize varieties. A key method for tracking nitrogen uptake is the use of the 15N stable isotope, which is naturally less abundant than the main 14N isotope. This method involves replacing the 14N in nutrient solutions with 15N, exogenously providing it to the plant tissues (roots in this case), and then measuring the 15N content of the tissues after a fixed amount of time. Here, we provide a brief overview of nitrogen uptake and remobilization in maize, and discuss current techniques for measuring nutrient uptake, with a focus on methods using stable isotopes of nitrogen.

Single-cell RNA sequencing reveals a key regulator ZmEREB14 affecting shoot apex development and yield formation in maize

Shoot apical meristem (SAM) is the origin of aerial structure formation in the plant life cycle. However, the mechanisms underlying the maize SAM development are still obscure. Here, approximately 12 700 cells were captured from the 5-day-old shoot apex of maize using a high-throughput single-cell transcriptome sequencing. According to the gene expression patterns, we partitioned the cells into 8 cell types with 13 transcriptionally distinct cell clusters and traced the developmental trajectory of shoot apex. Regulatory network analysis of transcription factors (TFs) showed that three core TFs, AP2-EREBP-transcription factor 14 (ZmEREB14, Zm00001d052087), MYB histone 4 (ZmMYB4) and HSF-transcription factor 8 (ZmHSF8) potentially regulated the SAM development. Functional validation revealed that ZmEREB14 affected the SAM development and thereby regulated the maize yield formation. Our results characterised the inherent heterogeneity of SAM at single-cell resolution and provided new insights into the mechanisms of SAM development.

ZmEREB25 transcription factor mediates transactivation of core starch synthetic genes in maize endosperm via interaction with ZmARF27

Starch, as the primary storage material in maize endosperm, is essential in determining yield and quality. Although the starch biosynthetic pathway in maize has been well-documented, the transcriptional network underlying endosperm starch synthesis remains elusive. Through a comprehensive co-expression analysis, we screened an endosperm-preferential AP2/ERF transcription factor ZmEREB25, which exhibited a strong correlation with the expression pattern of starch biosynthetic genes in maize endosperm. ZmEREB25 enhanced the promoter activities of the core starch biosynthetic genes, namely Sh2, SSIIIa and SSI, through specific binding to the GCCGAC-containing elements present in their promoters. Given that ZmEREB25 lacked the transactivation capacity, we further identified an ARF transcription factor, ZmARF27, that interacted with ZmEREB25 to coordinately transactivate the promoters of Sh2, SSIIIa and SSI genes via direct binding to these promoters. Our present study demonstrated that the ZmEREB25-ZmARF27 complex is crucial for transactivating core starch synthetic genes in maize endosperm and uncovered a novel regulatory pathway for starch synthesis in maize endosperm.

The regulatory potential of transposable elements in maize

The genomes of flowering plants consist largely of transposable elements (TEs), some of which modulate gene regulation and function. However, the repetitive nature of TEs and difficulty of mapping individual TEs by short-read-sequencing have hindered our understanding of their regulatory potential. We demonstrate that long-read chromatin fiber sequencing (Fiber-seq) comprehensively identifies accessible chromatin regions (ACRs) and CpG methylation across the maize genome. We uncover stereotypical ACR patterns at young TEs that degenerate with evolutionary age, resulting in TE-enhancers preferentially marked by a novel plant-specific epigenetic feature: simultaneous hyper-CpG methylation and chromatin accessibility. We show that TE ACRs are co-opted as gene promoters and that ACR-containing TEs can facilitate gene amplification. Lastly, we uncover a pervasive epigenetic signature - hypo-5mCpG methylation and diffuse chromatin accessibility - directing TEs to specific loci, including the loci that sparked McClintock's discovery of TEs.

Characterization and functional analysis of conserved non-coding sequences among poaceae: insights into gene regulation and phenotypic variation in maize

BACKGROUND

Conserved non-coding sequences (CNS) are islands of non-coding sequences conserved across species and play an important role in regulating the spatiotemporal expression of genes. Identification of CNS provides valuable information about potentially functional genomic elements, regulatory regions, and helps to gain insights into the genetic basis of crop agronomic traits.

RESULTS

Here, we comprehensively analyze CNS in maize, by comparing the genomes of maize inbred line B73 (Zea mays ssp. mays), its close wild relative Zea mays spp. mexicana, and other grasses in Poaceae, including sorghum (Sorghum bicolor), foxtail millet (Setaria italica) and two adlay (Coix lacryma) cultivars. There were 289,931 CNS found in two syntenic gene pairs, while 51,701 CNS were conserved within at least three species. To explore the regulatory characteristics of the CNS identified, the flanking regions of CNS were compared with the peaks called using both transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and chromatin immunoprecipitation with high-throughput sequencing (ChIP-Seq) data of histone modifications. It was found that CNS in maize were enriched in open chromatin regions compared with randomly selected non-coding regions of similar length. A significant enrichment of transcription factor binding sites was found within CNS sequences, including different transcription factors involved in abiotic stress response, such as OBP (OBF-BINDING PROTEIN) family and Adof1 (Encodes dof zinc finger protein). To investigate the epigenetic modification patterns in CNS, ChIP-Seq data for histone modifications H3K9ac, H3K4me3, H3K36me3, H3K9me3, and H3K27ac were further analyzed to depict the changes along CNS. Our findings revealed significantly elevated levels of transcription-promoting histone modifications in the CNS regions compared to randomly selected non-coding sequences with an equal number and similar length. Notably, CNS were also identified on both Vgt1 (Vegetative to generative transition 1) and ZmCCT10. In addition, CNS with potential functions were identified based on SNPs within CNS significantly associated with various agronomic traits in maize, which holds potential utility in molecular breeding for maize.

CONCLUSIONS

In summary, we identified and characterized CNS in maize through genomic comparative analysis, which provides valuable insights into their potential regulatory effects on gene expression and phenotypic variation.

Global Transcriptomic Analysis of Inbred Lines Reveal Candidate Genes for Response to Maize Lethal Necrosis

Maize lethal necrosis (MLN) is a significant threat to food security in Sub-Saharan Africa (SSA), with limited commercial inbred lines displaying tolerance. This study analyzed the transcriptomes of four commercially used maize inbred lines and a non-adapted inbred line, all with varying response levels to MLN. RNA-Seq revealed differentially expressed genes in response to infection by maize chlorotic mottle virus (MCMV) and sugarcane mosaic virus (SCMV), the causative agents of MLN. Key findings included the identification of components of the plant innate immune system, such as differentially regulated R genes (mainly LRRs), and activation/deactivation of virus resistance pathways, including RNA interference (RNAi) via Argonaute (AGO), Dicer-like proteins, and the ubiquitin-proteasome system (UPS) via RING/U-box and ubiquitin ligases. Genes associated with redox signaling, WRKY transcription factors, and cell modification were also differentially expressed. Additionally, the expression of translation initiation and elongation factors, eIF4E and eIF4G, correlated with the presence of MLN viruses. These findings provide valuable insights into the molecular mechanisms of MLN resistance and highlight potential gene candidates for engineering or selecting MLN-resistant maize germplasm for SSA.

A Zea genus-specific micropeptide controls kernel dehydration in maize

Kernel dehydration rate (KDR) is a crucial production trait that affects mechanized harvesting and kernel quality in maize; however, the underlying mechanisms remain unclear. Here, we identified a quantitative trait locus (QTL), qKDR1, as a non-coding sequence that regulates the expression of qKDR1 REGULATED PEPTIDE GENE (RPG). RPG encodes a 31 amino acid micropeptide, microRPG1, which controls KDR by precisely modulating the expression of two genes, ZmETHYLENE-INSENSITIVE3-like 1 and 3, in the ethylene signaling pathway in the kernels after filling. microRPG1 is a Zea genus-specific micropeptide and originated de novo from a non-coding sequence. Knockouts of microRPG1 result in faster KDR in maize. By contrast, overexpression or exogenous application of the micropeptide shows the opposite effect both in maize and Arabidopsis. Our findings reveal the molecular mechanism of microRPG1 in kernel dehydration and provide an important tool for future crop breeding.

Expression divergence of BAG gene family in maize under heat stress

Heat stress poses a significant challenge for maize production, especially during the spring when high temperatures disrupt cellular processes, impeding plant growth and development. The B-cell lymphoma-2 (Bcl-2) associated athanogene (BAG) gene family is known to be relatively conserved across various species. It plays a crucial role as molecular chaperone cofactors that are responsible for programmed cell death and tumorigenesis. Once the plant is under heat stress, the BAG genes act as co-chaperones and modulate the molecular functions of HSP70/HSC70 saving the plant from the damage of high temperature stress. The study was planned to identify and characterize the BAG genes for heat stress responsiveness in maize. Twenty-one (21) BAG genes were identified in the latest maize genome. The evolutionary relationship of Zea mays BAGs (ZmBAGs) with Arabidopsis thaliana, Solanum lycopersicum, Theobroma cacao, Sorghum bicolor, Ananas comosus, Physcomitrium patens, Oryza sativa and Populus trichocarpa were represented by the phylogenetic analysis. Differential expressions of BAG gene family in leaf, endosperm, anther, silk, seed and developing embryo depict their contribution to the growth and development. The in-silico gene expression analysis indicated ZmBAG-8 (Zm00001eb170080), and ZmBAG-11 (Zm00001eb237960) showed higher expression under abiotic stresses (cold, heat and salinity). The RT-qPCR further confirmed the expression of ZmBAG-8 and ZmBAG-11 in plant leaf tissue across the contrasting inbred lines and their F1 hybrid (DR-139, UML-1 and DR-139 × UML-1) when exposed to heat stress. Furthermore, the protein-protein interaction networks of ZmBAG-8 and ZmBAG-11 further elucidated their role in stress tolerance related pathways. This research offers a roadmap to plan functional research and utilize ZmBAG genes to enhance heat tolerance in grasses.

Cold mediates maize root hair developmental plasticity via epidermis-specific transcriptomic responses

Cold stress during early development limits maize (Zea mays L.) production in temperate zones. Low temperatures restrict root growth and reprogram gene expression. Here, we provide a systematic transcriptomic landscape of maize primary roots, their tissues, and cell types in response to cold stress. The epidermis exhibited a unique transcriptomic cold response, and genes involved in root hair formation were dynamically regulated in this cell type by cold. Consequently, activation of genes involved in root hair tip growth contributed to root hair recovery under moderate cold conditions. The maize root hair defective mutants roothair defective 5 (rth5) and roothair defective 6 (rth6) displayed enhanced cold tolerance with respect to primary root elongation. Furthermore, DEHYDRATION RESPONSE ELEMENT-BINDING PROTEIN 2.1 (DREB2.1) was the only member of the dreb subfamily of AP2/EREB transcription factor genes upregulated in primary root tissues and cell types but exclusively downregulated in root hairs upon cold stress. Plants overexpressing dreb2.1 significantly suppressed root hair elongation after moderate cold stress. Finally, the expression of rth3 was regulated by dreb2.1 under cold conditions, while rth6 transcription was regulated by DREB2.1 irrespective of the temperature regime. We demonstrated that dreb2.1 negatively regulates root hair plasticity at low temperatures by coordinating the expression of root hair defective genes in maize.

ZmPHR1 contributes to drought resistance by modulating phosphate homeostasis in maize

As an essential macronutrient, phosphorus (P) is often a limiting nutrient because of its low availability and mobility in soils. Drought is a major environmental stress that reduces crop yield. How plants balance and combine P-starvation responses (PSRs) and drought resistance is unclear. In this study, we identified the transcription factor ZmPHR1 as a major regulator of PSRs that modulates phosphate (Pi) signaling and homeostasis. We found that maize zmphr1 mutants had reduced P concentration and were sensitive to Pi starvation, whereas ZmPHR1-OE lines displayed elevated Pi concentration and yields. In addition, 57% of PSR genes and nearly 70% of ZmPHR1-regulated PSR genes in leaves were transcriptionally responsive to drought. Under moderate and early drought conditions, the Pi concentration of maize decreased, and PSR genes were up-regulated before drought-responsive genes. The ZmPHR1-OE lines exhibited drought-resistant phenotypes and reduced stomatal apertures, whereas the opposite was true of the zmphr1 mutants. ZmPT7-OE lines and zmspx3 mutants, which had elevated Pi concentration, also exhibited drought resistance, but zmpt7 mutants were sensitive to drought. Our results suggest that ZmPHR1 plays a central role in integrating Pi and drought signals and that Pi homeostasis improves the ability of maize to combat drought.

Genome-wide association study reveals genetic basis and candidate genes for chlorophyll content of leaves in maize (Zea mays L.)

The chlorophyll content (CC) directly affects photosynthesis, growth, and yield. However, the genetic basis of CC is still unclear in maize (Zea mays L.). Here, we conducted a genome-wide association study using mixed linear model for CC of the fifth leaves at seedling stage (CCFSS) and the ear leaves at filling stage (CCEFS) for 334 maize inbred lines. The heritability estimates for CCFSS and CCEFS, obtained via variance components analysis using the lme4 package in R, were 70.84% and 78.99%, respectively, indicating that the CC of leaves is primarily controlled by genetic factors. A total of 15 CC-related SNPs and 177 candidate genes were identified with a p-value < 4.49 × 10-5, which explained 4.98-7.59% of the phenotypic variation. Lines with more favorable gene variants showed higher CC. Meanwhile, Gene Ontology (GO) analysis implied that these candidate genes were probably related to chlorophyll biosynthesis. In addition, gene-based association analyses revealed that six variants in GRMZM2G037152, GRMZM5G816561, GRMZM2G324462, and GRMZM2G064657 genes were significantly (p-value < 0.01) correlated with CC, of which GRMZM2G064657 (encodes a phosphate transporter protein) and GRMZM5G816561 (encodes a cytochrome P450 protein) were specifically highly expressed in leaves tissues. Interestingly, these candidate genes were previously reported to involve in the regulation of the contents of chlorophyll in plants or Chlamydomonas. These results may contribute to the understanding of genetic basis and molecular mechanisms of maize CC and the selection of maize varieties with improved CC.

Functional Optimization in Distinct Tissues and Conditions Constrains the Rate of Protein Evolution

Understanding the main determinants of protein evolution is a fundamental challenge in biology. Despite many decades of active research, the molecular and cellular mechanisms underlying the substantial variability of evolutionary rates across cellular proteins are not currently well understood. It also remains unclear how protein molecular function is optimized in the context of multicellular species and why many proteins, such as enzymes, are only moderately efficient on average. Our analysis of genomics and functional datasets reveals in multiple organisms a strong inverse relationship between the optimality of protein molecular function and the rate of protein evolution. Furthermore, we find that highly expressed proteins tend to be substantially more functionally optimized. These results suggest that cellular expression costs lead to more pronounced functional optimization of abundant proteins and that the purifying selection to maintain high levels of functional optimality significantly slows protein evolution. We observe that in multicellular species both the rate of protein evolution and the degree of protein functional efficiency are primarily affected by expression in several distinct cell types and tissues, specifically, in developed neurons with upregulated synaptic processes in animals and in young and fast-growing tissues in plants. Overall, our analysis reveals how various constraints from the molecular, cellular, and species' levels of biological organization jointly affect the rate of protein evolution and the level of protein functional adaptation.

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.

GWAS and Meta-QTL Analysis of Kernel Quality-Related Traits in Maize

The quality of corn kernels is crucial for their nutritional value, making the enhancement of kernel quality a primary objective of contemporary corn breeding efforts. This study utilized 260 corn inbred lines as research materials and assessed three traits associated with grain quality. A genome-wide association study (GWAS) was conducted using the best linear unbiased estimator (BLUE) for quality traits, resulting in the identification of 23 significant single nucleotide polymorphisms (SNPs). Additionally, nine genes associated with grain quality traits were identified through gene function annotation and prediction. Furthermore, a total of 697 quantitative trait loci (QTL) related to quality traits were compiled from 27 documents, followed by a meta-QTL analysis that revealed 40 meta-QTL associated with these traits. Among these, 19 functional genes and reported candidate genes related to quality traits were detected. Three significant SNPs identified by GWAS were located within the intervals of these QTL, while the remaining eight significant SNPs were situated within 2 Mb of the QTL. In summary, the findings of this study provide a theoretical framework for analyzing the genetic basis of corn grain quality-related traits and for enhancing corn quality.

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.

An atlas of Brachypodium distachyon lateral root development

The root system of plants is a vital part for successful development and adaptation to different soil types and environments. A major determinant of the shape of a plant root system is the formation of lateral roots, allowing for expansion of the root system. Arabidopsis thaliana, with its simple root anatomy, has been extensively studied to reveal the genetic program underlying root branching. However, to get a more general understanding of lateral root development, comparative studies in species with a more complex root anatomy are required. Here, by combining optimized clearing methods and histology, we describe an atlas of lateral root development in Brachypodium distachyon, a wild, temperate grass species. We show that lateral roots initiate from enlarged phloem pole pericycle cells and that the overlying endodermis reactivates its cell cycle and eventually forms the root cap. In addition, auxin signaling reported by the DR5 reporter was not detected in the phloem pole pericycle cells or young primordia. In contrast, auxin signaling was activated in the overlying cortical cell layers, including the exodermis. Thus, Brachypodium is a valuable model to investigate how signaling pathways and cellular responses have been repurposed to facilitate lateral root organogenesis.

Brassinosteroid biosynthesis and signaling: Conserved and diversified functions of core genes across multiple plant species

Brassinosteroids (BRs) are important regulators that control myriad aspects of plant growth and development, including biotic and abiotic stress responses, such that modulating BR homeostasis and signaling presents abundant opportunities for plant breeding and crop improvement. Enzymes and other proteins involved in the biosynthesis and signaling of BRs are well understood from molecular genetics and phenotypic analysis in Arabidopsis thaliana; however, knowledge of the molecular functions of these genes in other plant species, especially cereal crop plants, is minimal. In this manuscript, we comprehensively review functional studies of BR genes in Arabidopsis, maize, rice, Setaria, Brachypodium, and soybean to identify conserved and diversified functions across plant species and to highlight cases for which additional research is in order. We performed phylogenetic analysis of gene families involved in the biosynthesis and signaling of BRs and re-analyzed publicly available transcriptomic data. Gene trees coupled with expression data provide a valuable guide to supplement future research on BRs in these important crop species, enabling researchers to identify gene-editing targets for BR-related functional studies.

Exploring the Roles of the Swi2/Snf2 Gene Family in Maize Abiotic Stress Responses

The maize Snf2 gene family plays a crucial role in chromatin remodeling and response to environmental stresses. In this study, we identified and analyzed 35 members of the maize Snf2 gene family (ZmCHR1 to ZmCHR35) using the Ensembl Plants database. Each protein contained conserved SNF2-N and Helicase-C domains. Phylogenetic analysis revealed six groups among the Snf2 proteins, with an uneven distribution across subfamilies. Physicochemical analysis indicated that the Snf2 proteins are hydrophilic, with varied amino acid lengths, isoelectric points, and molecular weights, and are predominantly localized in the nucleus. Chromosomal mapping showed that these genes are distributed across all ten maize chromosomes. Gene structure analysis revealed diverse exon-intron arrangements, while motif analysis identified 20 conserved motifs. Collinearity analysis highlighted gene duplication events, suggesting purifying selection. Cis-regulatory element analysis suggested involvement in abiotic and biotic stress responses. Expression analysis indicated tissue-specific expression patterns and differential expression under various stress conditions. Specifically, qRT-PCR validation under drought stress showed that certain Snf2 genes were upregulated at 12 h and downregulated at 24 h, revealing potential roles in drought tolerance. These findings provide a foundation for further exploration of the functional roles of the maize Snf2 gene family in development and stress responses.

Nonphotochemical quenching kinetics GWAS in sorghum identifies genes that may play conserved roles in maize and Arabidopsis thaliana photoprotection

Photosynthetic organisms must cope with rapid fluctuations in light intensity. Nonphotochemical quenching (NPQ) enables the dissipation of excess light energy as heat under high light conditions, whereas its relaxation under low light maximizes photosynthetic productivity. We quantified variation in NPQ kinetics across a large sorghum (Sorghum bicolor) association panel in four environments, uncovering significant genetic control for NPQ. A genome-wide association study (GWAS) confidently identified three unique regions in the sorghum genome associated with NPQ and suggestive associations in an additional 61 regions. We detected strong signals from the sorghum ortholog of Arabidopsis thaliana Suppressor Of Variegation 3 (SVR3) involved in plastid-nucleus signaling. By integrating GWAS results for NPQ across maize (Zea mays) and sorghum-association panels, we identified a second gene, Non-yellowing 1 (NYE1), originally studied by Gregor Mendel in pea (Pisum sativum) and involved in the degradation of photosynthetic pigments in light-harvesting complexes. Analysis of nye1 insertion alleles in A. thaliana confirmed the effect of this gene on NPQ kinetics in eudicots. We extended our comparative genomics GWAS framework across the entire maize and sorghum genomes, identifying four additional loci involved in NPQ kinetics. These results provide a baseline for increasing the accuracy and speed of candidate gene identification for GWAS in species with high linkage disequilibrium.

Integrating GWAS with a gene co-expression network better prioritizes candidate genes associated with root metaxylem phenes in maize

Root metaxylems are phenotypically diverse structures whose function is particularly important under drought stress. Significant research has dissected the genetic machinery underlying metaxylem phenotypes in dicots, but that of monocots are relatively underexplored. In maize (Zea mays), a robust pipeline integrated a genome-wide association study (GWAS) of root metaxylem phenes under well-watered and water-stress conditions with a gene co-expression network to prioritize the strongest gene candidates. We identified 244 candidate genes by GWAS, of which 103 reside in gene co-expression modules most relevant to xylem development. Several candidate genes may be involved in biosynthetic processes related to the cell wall, hormone signaling, oxidative stress responses, and drought responses. Of those, six gene candidates were detected in multiple root metaxylem phenes in both well-watered and water-stress conditions. We posit that candidate genes that are more essential to network function based on gene co-expression (i.e., hubs or bottlenecks) should be prioritized and classify 33 essential genes for further investigation. Our study demonstrates a new strategy for identifying promising gene candidates and presents several gene candidates that may enhance our understanding of vascular development and responses to drought in cereals.

Effect of linolenic acid on triterpenoids production by the liquid fermentation of Antrodia cinnamomea

Liquid state fermentation is now a commonly used route to obtain triterpenoids from Antrodia cinnamomea, and linolenic acid can significantly promote triterpenoids synthesis, whereas its action mechanism has not been studied. Here, we comprehensively performed an investigation on the mechanism of linolenic acid to promote triterpenoids production in liquid-state fermentation of A. cinnamomea. Results showed that the addition of linolenic acid increased the unsaturated fatty acid index, fluidity, and permeability in the cell membrane of A. cinnamomea mycelia, favored the absorption of nutrients in the medium by the mycelium, enhanced the material exchange inside and outside, and thus promoted mycelial growth and triterpenoids synthesis. Moreover, 767 significantly differentially expressed genes were detected by adding linolenic acid, including 212 upregulated genes and 555 downregulated genes. The upregulated genes were mainly enriched in metabolism, glycolytic pathway, TCA cycle, and pyruvate metabolism. It was seen that the addition of linolenic acid improved the cell metabolic activity and promoted the synthesis of secondary metabolites, proving that the addition of linolenic acid improved the metabolic viability of cells and promoted secondary metabolite synthesis.

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

Solving the Jigsaw puzzle of phytosterol diversity by a novel sterol methyltransferase from Zea mays

Phytosterols are vital structural and regulatory components in plants. Zea mays produces a series of phytosterols that are specific to corn. However, the underline biosynthetic mechanism remains elusive. In this study, we identified a novel sterol methyltransferase from Z. mays (ZmSMT1-2) which showed a unique feature compared with documented plant SMTs. ZmSMT1-2 showed a substrate preference for cycloartenol. Using S-adenosyl-L-methionine (AdoMet) as a donor, ZmSMT1-2 converted cycloartenol into alkylated sterols with unique side-chain architectures, including Δ25(27) (i.e., cyclolaudenol and cycloneolitsol) and Δ24(25) (i.e., cyclobranol) sterols. Cycloneolitsol is identified as a product of SMTs for the first time. Our discovery provides a previously untapped mechanism for phytosterol biosynthesis and adds another layer of diversity of sterol biosynthesis.

Genome-Wide Association Study Identifies Quantitative Trait Loci and Candidate Genes Involved in Deep-Sowing Tolerance in Maize (Zea mays L.)

Deep sowing is an efficient strategy for maize to ensure the seedling emergence rate under adverse conditions such as drought or low temperatures. However, the genetic basis of deep-sowing tolerance-related traits in maize remains largely unknown. In this study, we performed a genome-wide association study on traits related to deep-sowing tolerance, including mesocotyl length (ML), coleoptile length (CL), plumule length (PL), shoot length (SL), and primary root length (PRL), using 255 maize inbred lines grown in three different environments. We identified 23, 6, 4, and 4 quantitative trait loci (QTLs) associated with ML, CL, PL, and SL, respectively. By analyzing candidate genes within these QTLs, we found a γ-tubulin-containing complex protein, ZmGCP2, which was significantly associated with ML, PL, and SL. Loss of function of ZmGCP2 resulted in decreased PL, possibly by affecting the cell elongation, thus affecting SL. Additionally, we identified superior haplotypes and allelic variations of ZmGCP2 with a longer PL and SL, which may be useful for breeding varieties with deep-sowing tolerance to improve maize cultivation.

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

Meta-Quantitative Trait Loci Analysis and Candidate Gene Mining for Drought Tolerance-Associated Traits in Maize (Zea mays L.)

Drought is one of the major abiotic stresses with a severe negative impact on maize production globally. Understanding the genetic architecture of drought tolerance in maize is a crucial step towards the breeding of drought-tolerant varieties and a targeted exploitation of genetic resources. In this study, 511 quantitative trait loci (QTL) related to grain yield components, flowering time, and plant morphology under drought conditions, as well as drought tolerance index were collected from 27 published studies and then projected on the IBM2 2008 Neighbors reference map for meta-analysis. In total, 83 meta-QTL (MQTL) associated with drought tolerance in maize were identified, of which 20 were determined as core MQTL. The average confidence interval of MQTL was strongly reduced compared to that of the previously published QTL. Nearly half of the MQTL were confirmed by co-localized marker-trait associations from genome-wide association studies. Based on the alignment of rice proteins related to drought tolerance, 63 orthologous genes were identified near the maize MQTL. Furthermore, 583 candidate genes were identified within the 20 core MQTL regions and maize-rice homologous genes. Based on KEGG analysis of candidate genes, plant hormone signaling pathways were found to be significantly enriched. The signaling pathways can have direct or indirect effects on drought tolerance and also interact with other pathways. In conclusion, this study provides novel insights into the genetic and molecular mechanisms of drought tolerance in maize towards a more targeted improvement of this important trait in breeding.

Comparison of the transcriptome and metabolome of wheat (Triticum aestivum L.) proteins content during grain formation provides insight

Introduction

Wheat is a food crop with a large global cultivation area, and the content and quality of wheat glutenin accumulation are important indicators of the quality of wheat flour.

Methods

To elucidate the gene expression regulation and metabolic characteristics related to the gluten content during wheat grain formation, transcriptomic and metabolomic analyses were performed for the high gluten content of the Xinchun 26 cultivar and the low proteins content of the Xinchun 34 cultivar at three periods (7 d, 14 d and 21 d) after flowering.

Results

Transcriptomic analysis revealed that 5573 unique differentially expressed genes (DEGs) were divided into two categories according to their expression patterns during the three periods. The metabolites detected were mainly divided into 12 classes. Lipid and lipid-like molecule levels and phenylpropanoid and polyketide levels were the highest, and the difference analysis revealed a total of 10 differentially regulated metabolites (DRMs) over the three periods. Joint analysis revealed that the DEGs and DRMs were significantly enriched in starch and sucrose metabolism; the citrate cycle; carbon fixation in photosynthetic organisms; and alanine, aspartate and glutamate metabolism pathways. The genes and contents of the sucrose and gluten synthesis pathways were analysed, and the correlation between gluten content and its related genes was calculated. Based on weighted correlation network analysis (WGCNA), by constructing a coexpression network, a total of 5 specific modules and 8 candidate genes that were strongly correlated with the three developmental stages of wheat grain were identified.

Discussion

This study provides new insights into the role of glutenin content in wheat grain formation and reveals potential regulatory pathways and candidate genes involved in this developmental process.

Genetic analysis of pericarp pigmentation variation in Corn Belt dent maize

The US standard for maize commercially grown for grain specifies that yellow corn can contain at maximum 5% corn of other colors. Inbred parents of commercial hybrids typically have clear pericarp, but transgressive segregants in breeding populations can display variation in pericarp pigmentation. We identified 10 doubled haploid biparental populations segregating for pigmented pericarp and evaluated qualitative genetic models using chi-square tests of observed and expected frequencies. Pigmentation ranged from light to dark brown color, and pigmentation intensity was quantitatively measured across 1,327 inbred lines using hue calculated from RGB pixel values. Genetic mapping was used to identify loci associated with pigmentation intensity. For 9 populations, pigmentation inheritance best fit a hypothesis of a 2- or 3-gene epistatic model. Significant differences in pigment intensity were observed across populations. W606S-derived inbred lines with the darkest pericarp often had clear glumes, suggesting the presence of a novel P1-rw allele, a hypothesis supported by a significant quantitative trait locus peak at P1. A separate quantitative trait locus region on chromosome 2 between 221.64 and 226.66 Mbp was identified in LH82-derived populations, and the peak near p1 was absent. A genome-wide association study using 416 inbred lines from the Wisconsin Diversity panel with full genome resequencing revealed 4 significant associations including the region near P1. This study supports that pericarp pigmentation among dent maize inbreds can arise by transgressive segregation when pigmentation in the parental generation is absent and is partially explained by functional allelic variation at the P1 locus.

Maize Dek407 Encodes the Nitrate Transporter 1.5 and Is Required for Kernel Development

The kernel serves as the storage organ and harvestable component of maize, and it plays a crucial role in determining crop yield and quality. Understanding the molecular and genetic mechanisms of kernel development is of considerable importance for maize production. In this study, we obtained a mutant, which we designated defective kernel 407 (dek407), through ethyl methanesulfonate mutagenesis. The dek407 mutant exhibited reduced kernel size and kernel weight, as well as delayed grain filling compared with those of the wild type. Positional cloning and an allelism test revealed that Dek407 encodes a nitrate transporter 1/peptide transporter family (NPF) protein and is the allele of miniature 2 (mn2) that was responsible for a poorly filled defective kernel phenotype. A transcriptome analysis of the developing kernels showed that the mutation of Dek407 altered the expression of phytohormone-related genes, especially those genes associated with indole-3-acetic acid synthesis and signaling. Phytohormone measurements and analysis indicated that the endogenous indole-3-acetic acid content was significantly reduced by 66% in the dek407 kernels, which may be the primary cause of the defective phenotype. We further demonstrated that natural variation in Dek407 is associated with kernel weight and kernel size. Therefore, Dek407 is a potential target gene for improvement of maize yield.

A role for fermentation in aerobic conditions as revealed by computational analysis of maize root metabolism during growth by cell elongation

The root is a well-studied example of cell specialisation, yet little is known about the metabolism that supports the transport functions and growth of different root cell types. To address this, we used computational modelling to study metabolism in the elongation zone of a maize lateral root. A functional-structural model captured the cell-anatomical features of the root and modelled how they changed as the root elongated. From these data, we derived constraints for a flux balance analysis model that predicted metabolic fluxes of the 11 concentric rings of cells in the root. We discovered a distinct metabolic flux pattern in the cortical cell rings, endodermis and pericycle (but absent in the epidermis) that involved a high rate of glycolysis and production of the fermentation end-products lactate and ethanol. This aerobic fermentation was confirmed experimentally by metabolite analysis. The use of fermentation in the model was not obligatory but was the most efficient way to meet the specific demands for energy, reducing power and carbon skeletons of expanding cells. Cytosolic acidification was avoided in the fermentative mode due to the substantial consumption of protons by lipid synthesis. These results expand our understanding of fermentative metabolism beyond that of hypoxic niches and suggest that fermentation could play an important role in the metabolism of aerobic tissues.

GWAS and Meta-QTL Analysis of Yield-Related Ear Traits in Maize

Maize ear traits are an important component of yield, and the genetic basis of ear traits facilitates further yield improvement. In this study, a panel of 580 maize inbred lines were used as the study material, eight ear-related traits were measured through three years of planting, and whole genome sequencing was performed using the maize 40 K breeding chip based on genotyping by targeted sequencing (GBTS) technology. Five models were used to conduct a genome-wide association study (GWAS) on best linear unbiased estimate (BLUE) of ear traits to find the best model. The FarmCPU (Fixed and random model Circulating Probability Unification) model was the best model for this study; a total of 104 significant single nucleotide polymorphisms (SNPs) were detected, and 10 co-location SNPs were detected simultaneously in more than two environments. Through gene function annotation and prediction, a total of nine genes were identified as potentially associated with ear traits. Moreover, a total of 760 quantitative trait loci (QTL) associated with yield-related traits reported in 37 different articles were collected. Using the collected 760 QTL for meta-QTL analysis, a total of 41 MQTL (meta-QTL) associated with yield-related traits were identified, and 19 MQTL detected yield-related ear trait functional genes and candidate genes that have been reported in maize. Five significant SNPs detected by GWAS were located within these MQTL intervals, and another three significant SNPs were close to MQTL (less than 1 Mb). The results provide a theoretical reference for the analysis of the genetic basis of ear-related traits and the improvement of maize yield.

The DEAD-box RNA helicase ZmRH48 is required for the splicing of multiple mitochondrial introns, mitochondrial complex biosynthesis, and seed development in maize

RNA helicases participate in nearly all aspects of RNA metabolism by rearranging RNAs or RNA-protein complexes in an adenosine triphosphate-dependent manner. Due to the large RNA helicase families in plants, the precise roles of many RNA helicases in plant physiology and development remain to be clarified. Here, we show that mutations in maize (Zea mays) DEAD-box RNA helicase 48 (ZmRH48) impair the splicing of mitochondrial introns, mitochondrial complex biosynthesis, and seed development. Loss of ZmRH48 function severely arrested embryogenesis and endosperm development, leading to defective kernel formation. ZmRH48 is targeted to mitochondria, where its deficiency dramatically reduced the splicing efficiency of five cis-introns (nad5 intron 1; nad7 introns 1, 2, and 3; and ccmFc intron 1) and one trans-intron (nad2 intron 2), leading to lower levels of mitochondrial complexes I and III. ZmRH48 interacts with two unique pentatricopeptide repeat (PPR) proteins, PPR-SMR1 and SPR2, which are required for the splicing of over half of all mitochondrial introns. PPR-SMR1 interacts with SPR2, and both proteins interact with P-type PPR proteins and Zm-mCSF1 to facilitate intron splicing. These results suggest that ZmRH48 is likely a component of a splicing complex and is critical for mitochondrial complex biosynthesis and seed development.

Genetic basis of maize maternal haploid induction beyond MATRILINEAL and ZmDMP

In maize, doubled haploid (DH) lines are created in vivo through crosses with maternal haploid inducers. Their induction ability, usually expressed as haploid induction rate (HIR), is known to be under polygenic control. Although two major genes (MTL and ZmDMP) affecting this trait were recently described, many others remain unknown. To identify them, we designed and performed a SNP based (~9007) genome-wide association study using a large and diverse panel of 159 maternal haploid inducers. Our analyses identified a major gene near MTL, which is present in all inducers and necessary to disrupt haploid induction. We also found a significant quantitative trait loci (QTL) on chromosome 10 using a case-control mapping approach, in which 793 noninducers were used as controls. This QTL harbors a kokopelli ortholog, whose role in maternal haploid induction was recently described in Arabidopsis. QTL with smaller effects were identified on six of the ten maize chromosomes, confirming the polygenic nature of this trait. These QTL could be incorporated into inducer breeding programs through marker-assisted selection approaches. Further improving HIR is important to reduce the cost of DH line production.

Investigating the genetic basis of maize ear characteristics: a comprehensive genome-wide study utilizing high-throughput phenotypic measurement method and system

The morphology of maize ears plays a critical role in the breeding of new varieties and increasing yield. However, the study of traditional ear-related traits alone can no longer meet the requirements of breeding. In this study, 20 ear-related traits, including size, shape, number, and color, were obtained in 407 maize inbred lines at two sites using a high-throughput phenotypic measurement method and system. Significant correlations were found among these traits, particularly the novel trait ear shape (ES), which was correlated with traditional traits: kernel number per row and kernel number per ear. Pairwise comparison tests revealed that the inbred lines of tropical-subtropical were significantly different from other subpopulations in row numbers per ear, kernel numbers per ear, and ear color. A genome-wide association study identified 275, 434, and 362 Single nucleotide polymorphisms (SNPs) for Beijing, Sanya, and best linear unbiased prediction scenarios, respectively, explaining 3.78% to 24.17% of the phenotypic variance. Furthermore, 58 candidate genes with detailed functional descriptions common to more than two scenarios were discovered, with 40 genes being associated with color traits on chromosome 1. After analysis of haplotypes, gene expression, and annotated information, several candidate genes with high reliability were identified, including Zm00001d051328 for ear perimeter and width, zma-MIR159f for ear shape, Zm00001d053080 for kernel width and row number per ear, and Zm00001d048373 for the blue color channel of maize kernels in the red-green-blue color model. This study emphasizes the importance of researching novel phenotypic traits in maize by utilizing high-throughput phenotypic measurements. The identified genetic loci enrich the existing genetic studies related to maize ears.

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.

Combining multiplex gene editing and doubled haploid technology in maize

A major advantage of using CRISPR/Cas9 for gene editing is multiplexing, that is, the simultaneous targeting of many genes. However, primary transformants typically contain hetero-allelic mutations or are genetic mosaic, while genetically stable lines that are homozygous are desired for functional analysis. Currently, a dedicated and labor-intensive effort is required to obtain such higher-order mutants through several generations of genetic crosses and genotyping. We describe the design and validation of a rapid and efficient strategy to produce lines of genetically identical plants carrying various combinations of homozygous edits, suitable for replicated analysis of phenotypical differences. This approach was achieved by combining highly multiplex gene editing in Zea mays (maize) with in vivo haploid induction and efficient in vitro generation of doubled haploid plants using embryo rescue doubling. By combining three CRISPR/Cas9 constructs that target in total 36 genes potentially involved in leaf growth, we generated an array of homozygous lines with various combinations of edits within three generations. Several genotypes show a reproducible 10% increase in leaf size, including a septuple mutant combination. We anticipate that our strategy will facilitate the study of gene families via multiplex CRISPR mutagenesis and the identification of allele combinations to improve quantitative crop traits.

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.

A non-JA producing oxophytodienoate reductase functions in salicylic acid-mediated antagonism with jasmonic acid during pathogen attack

Peroxisome-localized oxo-phytodienoic acid (OPDA) reductases (OPR) are enzymes converting 12-OPDA into jasmonic acid (JA). However, the biochemical and physiological functions of the cytoplasmic non-JA producing OPRs remain largely unknown. Here, we generated Mutator-insertional mutants of the maize OPR2 gene and tested its role in resistance to pathogens with distinct lifestyles. Functional analyses showed that the opr2 mutants were more susceptible to the (hemi)biotrophic pathogens Colletotrichum graminicola and Ustilago maydis, but were more resistant to the necrotrophic fungus Cochliobolus heterostrophus. Hormone profiling revealed that increased susceptibility to C. graminicola was associated with decreased salicylic acid (SA) but increased JA levels. Mutation of the JA-producing lipoxygenase 10 (LOX10) reversed this phenotype in the opr2 mutant background, corroborating the notion that JA promotes susceptibility to this pathogen. Exogenous SA did not rescue normal resistance levels in opr2 mutants, suggesting that this SA-inducible gene is the key downstream component of the SA-mediated defences against C. graminicola. Disease assays of the single and double opr2 and lox10 mutants and the JA-deficient opr7opr8 mutants showed that OPR2 negatively regulates JA biosynthesis, and that JA is required for resistance against C. heterostrophus. Overall, this study uncovers a novel function of a non-JA producing OPR as a major negative regulator of JA biosynthesis during pathogen infection, a function that leads to its contrasting contribution to either resistance or susceptibility depending on pathogen lifestyle.

A dolabralexin-deficient mutant provides insight into specialized diterpenoid metabolism in maize

Two major groups of specialized metabolites in maize (Zea mays), termed kauralexins and dolabralexins, serve as known or predicted diterpenoid defenses against pathogens, herbivores, and other environmental stressors. To consider the physiological roles of the recently discovered dolabralexin pathway, we examined dolabralexin structural diversity, tissue-specificity, and stress-elicited production in a defined biosynthetic pathway mutant. Metabolomics analyses support a larger number of dolabralexin pathway products than previously known. We identified dolabradienol as a previously undetected pathway metabolite and characterized its enzymatic production. Transcript and metabolite profiling showed that dolabralexin biosynthesis and accumulation predominantly occur in primary roots and show quantitative variation across genetically diverse inbred lines. Generation and analysis of CRISPR-Cas9-derived loss-of-function Kaurene Synthase-Like 4 (Zmksl4) mutants demonstrated dolabralexin production deficiency, thus supporting ZmKSL4 as the diterpene synthase responsible for the conversion of geranylgeranyl pyrophosphate precursors into dolabradiene and downstream pathway products. Zmksl4 mutants further display altered root-to-shoot ratios and root architecture in response to water deficit. Collectively, these results demonstrate dolabralexin biosynthesis via ZmKSL4 as a committed pathway node biochemically separating kauralexin and dolabralexin metabolism, and suggest an interactive role of maize dolabralexins in plant vigor during abiotic stress.

Causes and consequences of endogenous hypoxia on growth and metabolism of developing maize kernels

Maize (Zea mays) kernels are the largest cereal grains, and their endosperm is severely oxygen deficient during grain fill. The causes, dynamics, and mechanisms of acclimation to hypoxia are minimally understood. Here, we demonstrate that hypoxia develops in the small, growing endosperm, but not the nucellus, and becomes the standard state, regardless of diverse structural and genetic perturbations in modern maize (B73, popcorn, sweet corn), mutants (sweet4c, glossy6, waxy), and non-domesticated wild relatives (teosintes and Tripsacum species). We also uncovered an interconnected void space at the chalazal pericarp, providing superior oxygen supply to the placental tissues and basal endosperm transfer layer. Modeling indicated a very high diffusion resistance inside the endosperm, which, together with internal oxygen consumption, could generate steep oxygen gradients at the endosperm surface. Manipulation of oxygen supply induced reciprocal shifts in gene expression implicated in controlling mitochondrial functions (23.6 kDa Heat-Shock Protein, Voltage-Dependent Anion Channel 2) and multiple signaling pathways (core hypoxia genes, cyclic nucleotide metabolism, ethylene synthesis). Metabolite profiling revealed oxygen-dependent shifts in mitochondrial pathways, ascorbate metabolism, starch synthesis, and auxin degradation. Long-term elevated oxygen supply enhanced the rate of kernel development. Altogether, evidence here supports a mechanistic framework for the establishment of and acclimation to hypoxia in the maize endosperm.

Mining genes regulating root system architecture in maize based on data integration analysis

KEY MESSAGE

A new strategy that integrated multiple public data resources was established to construct root gene co-expression network and mine genes regulating root system architecture in maize. A root gene co-expression network, containing 13,874 genes, was constructed. A total of 53 root hub genes and 16 priority root candidate genes were identified. One priority root candidate was further functionally verified using overexpression transgenic maize lines. Root system architecture (RSA) is crucial for crops productivity and stress tolerance. In maize, few RSA genes are functionally cloned, and effective discovery of RSA genes remains a great of challenge. In this work, we established a strategy to mine maize RSA genes by integrating functionally characterized root genes, root transcriptome, weighted gene co-expression network analysis (WGCNA) and genome-wide association analysis (GWAS) of RSA traits based on public data resources. A total of 589 maize root genes were collected by searching well-characterized root genes in maize or homologous genes of other species. We performed WGCNA to construct a maize root gene co-expression network containing 13874 genes based on public available root transcriptome data, and further discovered the 53 hub genes related to root traits. In addition, by the prediction function of obtained root gene co-expression network, a total of 1082 new root candidate genes were explored. By further overlapping the obtained new root candidate gene with the root-related GWAS of RSA candidate genes, 16 priority root candidate genes were identified. Finally, a priority root candidate gene, Zm00001d023379 (encodes pyruvate kinase 2), was validated to modulate root open angle and shoot-borne roots number using its overexpression transgenic lines. Our results develop an integration analysis method for effectively exploring regulatory genes of RSA in maize and open a new avenue to mine the candidate genes underlying complex traits.

Integrative transcriptome and metabolome analysis reveals the mechanism of exogenous melatonin alleviating drought stress in maize roots

Melatonin (MT) is essential for plant development and drought adaptation. However, the molecular and metabolic mechanisms underlying MT-induced drought tolerance in maize roots remain largely unclear. Herein, we investigated the effects of MT on drought tolerance in maize roots using integrated transcriptomic and metabolomic analyses, and identified MT-induced genes and metabolites associated with drought resistance. Compared with the untreated control plants, MT application alleviated the deleterious effects of drought on roots, by decreasing the malondialdehyde level and increasing the solute potential, eventually promoting root growth. Transcriptome and metabolome analysis demonstrated that MT significantly upregulates the expression of genes related to flavonoid biosynthesis (PAL, C4H, 4CL, HCT, CHS, CHI, F3'5'H, and DFR), activates drought-responsive transcription factors (ERFs, NACs, MYBs, and bHLHs), and regulates hormone signaling-related genes, especially ethylene response factors (ERF4, ERF81, and ERF110). Moreover, MT increased the accumulation of flavonoid metabolites, particularly apigenin, luteolin, and quercetin, under drought-stress conditions. These findings were further supported by quantitative real-time polymerase chain reaction analysis and total flavonoid measurements. Altogether, our findings suggest that MT promotes maize root growth during drought by regulating flavonoid synthesis pathways, transcription factors, and plant hormone signals. This study provides new insights into the complex mechanisms by which MT enhances crop resistance to drought damage.

Inference of differentially expressed genes using generalized linear mixed models in a pairwise fashion

Background

Technological advances involving RNA-Seq and Bioinformatics allow quantifying the transcriptional levels of genes in cells, tissues, and cell lines, permitting the identification of Differentially Expressed Genes (DEGs). DESeq2 and edgeR are well-established computational tools used for this purpose and they are based upon generalized linear models (GLMs) that consider only fixed effects in modeling. However, the inclusion of random effects reduces the risk of missing potential DEGs that may be essential in the context of the biological phenomenon under investigation. The generalized linear mixed models (GLMM) can be used to include both effects.

Methods

We present DEGRE (Differentially Expressed Genes with Random Effects), a user-friendly tool capable of inferring DEGs where fixed and random effects on individuals are considered in the experimental design of RNA-Seq research. DEGRE preprocesses the raw matrices before fitting GLMMs on the genes and the derived regression coefficients are analyzed using the Wald statistical test. DEGRE offers the Benjamini-Hochberg or Bonferroni techniques for P-value adjustment.

Results

The datasets used for DEGRE assessment were simulated with known identification of DEGs. These have fixed effects, and the random effects were estimated and inserted to measure the impact of experimental designs with high biological variability. For DEGs’ inference, preprocessing effectively prepares the data and retains overdispersed genes. The biological coefficient of variation is inferred from the counting matrices to assess variability before and after the preprocessing. The DEGRE is computationally validated through its performance by the simulation of counting matrices, which have biological variability related to fixed and random effects. DEGRE also provides improved assessment measures for detecting DEGs in cases with higher biological variability. We show that the preprocessing established here effectively removes technical variation from those matrices. This tool also detects new potential candidate DEGs in the transcriptome data of patients with bipolar disorder, presenting a promising tool to detect more relevant genes.

Conclusions

DEGRE provides data preprocessing and applies GLMMs for DEGs’ inference. The preprocessing allows efficient remotion of genes that could impact the inference. Also, the computational and biological validation of DEGRE has shown to be promising in identifying possible DEGs in experiments derived from complex experimental designs. This tool may help handle random effects on individuals in the inference of DEGs and presents a potential for discovering new interesting DEGs for further biological investigation.

Local auxin biosynthesis regulates brace root angle and lodging resistance in maize

Root lodging poses a major threat to maize production, resulting in reduced grain yield and quality, and increased harvest costs. Here, we combined expressional, genetic, and cytological studies to demonstrate a role of ZmYUC2 and ZmYUC4 in regulating gravitropic response of the brace root and lodging resistance in maize. We show that both ZmYUC2 and ZmYUC4 are preferentially expressed in root tips with partially overlapping expression patterns, and the protein products of ZmYUC2 and ZmYUC4 are localized in the cytoplasm and endoplasmic reticulum, respectively. The Zmyuc4 single mutant and Zmyuc2/4 double mutant exhibit enlarged brace root angle compared with the wild-type plants, with larger brace root angle being observed in the Zmyuc2/4 double mutant. Consistently, the brace root tips of the Zmyuc4 single mutant and Zmyuc2/4 double mutant accumulate less auxin and are defective in proper reallocation of auxin in response to gravi-stimuli. Furthermore, we show that the Zmyuc4 single mutant and the Zmyuc2/4 double mutant display obviously enhanced root lodging resistance. Our combined results demonstrate that ZmYUC2- and ZmYUC4-mediated local auxin biosynthesis is required for normal gravity response of the brace roots and provide effective targets for breeding root lodging resistant maize cultivars.

Comprehensive Analysis of the INDETERMINATE DOMAIN (IDD) Gene Family and Their Response to Abiotic Stress in Zea mays

Transcription factors (TFs) are important regulators of numerous gene expressions due to their ability to recognize and combine cis-elements in the promoters of target genes. The INDETERMINATE DOMAIN (IDD) gene family belongs to a subfamily of C2H2 zinc finger proteins and has been identified only in terrestrial plants. Nevertheless, little study has been reported concerning the genome-wide analysis of the IDD gene family in maize. In total, 22 ZmIDD genes were identified, which can be distributed on 8 chromosomes in maize. On the basis of evolutionary relationships and conserved motif analysis, ZmIDDs were categorized into three clades (1, 2, and 3), each owning 4, 6, and 12 genes, respectively. We analyzed the characteristics of gene structure and found that 3 of the 22 ZmIDD genes do not contain an intron. Cis-element analysis of the ZmIDD promoter showed that most ZmIDD genes possessed at least one ABRE or MBS cis-element, and some ZmIDD genes owned the AuxRR-core, TCA-element, TC-rich repeats, and LTR cis-element. The Ka:Ks ratio of eight segmentally duplicated gene pairs demonstrated that the ZmIDD gene families had undergone a purifying selection. Then, the transcription levels of ZmIDDs were analyzed, and they showed great differences in diverse tissues as well as abiotic stresses. Furthermore, regulatory networks were constructed through the prediction of ZmIDD-targeted genes and miRNAs, which can inhibit the transcription of ZmIDDs. In total, 6 ZmIDDs and 22 miRNAs were discovered, which can target 180 genes and depress the expression of 9 ZmIDDs, respectively. Taken together, the results give us valuable information for studying the function of ZmIDDs involved in plant development and climate resilience in maize.

Genome-Wide Meta-Analysis of QTLs Associated with Root Traits and Implications for Maize Breeding

Root system architecture (RSA), also known as root morphology, is critical in plant acquisition of soil resources, plant growth, and yield formation. Many QTLs associated with RSA or root traits in maize have been identified using several bi-parental populations, particularly in response to various environmental factors. In the present study, a meta-analysis of QTLs associated with root traits was performed in maize using 917 QTLs retrieved from 43 mapping studies published from 1998 to 2020. A total of 631 QTLs were projected onto a consensus map involving 19,714 markers, which led to the prediction of 68 meta-QTLs (MQTLs). Among these 68 MQTLs, 36 MQTLs were validated with the marker-trait associations available from previous genome-wide association studies for root traits. The use of comparative genomics approaches revealed several gene models conserved among the maize, sorghum, and rice genomes. Among the conserved genomic regions, the ortho-MQTL analysis uncovered 20 maize MQTLs syntenic to 27 rice MQTLs for root traits. Functional analysis of some high-confidence MQTL regions revealed 442 gene models, which were then subjected to in silico expression analysis, yielding 235 gene models with significant expression in various tissues. Furthermore, 16 known genes viz., DXS2, PHT, RTP1, TUA4, YUC3, YUC6, RTCS1, NSA1, EIN2, NHX1, CPPS4, BIGE1, RCP1, SKUS13, YUC5, and AW330564 associated with various root traits were present within or near the MQTL regions. These results could aid in QTL cloning and pyramiding in developing new maize varieties with specific root architecture for proper plant growth and development under optimum and abiotic stress conditions.

An ARF24-ZmArf2 module influences kernel size in different maize haplotypes

Members of the ADP-ribosylation factor family, which are GTP-binding proteins, are involved in metabolite transport, cell division, and expansion. Although there has been a significant amount of research on small GTP-binding proteins, their roles and functions in regulating maize kernel size remain elusive. Here, we identified ZmArf2 as a maize ADP-ribosylation factor-like family member that is highly conserved during evolution. Maize zmarf2 mutants showed a characteristic smaller kernel size. Conversely, ZmArf2 overexpression increased maize kernel size. Furthermore, heterologous expression of ZmArf2 dramatically elevated Arabidopsis and yeast growth by promoting cell division. Using expression quantitative trait loci (eQTL) analysis, we determined that ZmArf2 expression levels in various lines were mainly associated with variation at the gene locus. The promoters of ZmArf2 genes could be divided into two types, pS and pL, that were significantly associated with both ZmArf2 expression levels and kernel size. In yeast-one-hybrid screening, maize Auxin Response Factor 24 (ARF24) is directly bound to the ZmArf2 promoter region and negatively regulated ZmArf2 expression. Notably, the pS and pL promoter types each contained an ARF24 binding element: an auxin response element (AuxRE) in pS and an auxin response region (AuxRR) in pL, respectively. ARF24 binding affinity to AuxRR was much higher compared with AuxRE. Overall, our results establish that the small G-protein ZmArf2 positively regulates maize kernel size and reveals the mechanism of its expression regulation.