A Kinesin-14 Motor Activates Neocentromeres to Promote Meiotic Drive in Maize.

Maize abnormal chromosome 10 (Ab10) encodes a classic example of true meiotic drive that converts heterochromatic regions called knobs into motile neocentromeres that are preferentially transmitted to egg cells. Here, we identify a cluster of eight genes on Ab10, called the Kinesin driver (Kindr) complex, that are required for both neocentromere motility and preferential transmission. Two meiotic drive mutants that lack neocentromere activity proved to be kindr epimutants with increased DNA methylation across the entire gene cluster. RNAi of Kindr induced a third epimutant and corresponding loss of meiotic drive. Kinesin gliding assays and immunolocalization revealed that KINDR is a functional minus-end-directed kinesin that localizes specifically to knobs containing 180 bp repeats. Sequence comparisons suggest that Kindr diverged from a Kinesin-14A ancestor ∼12 mya and has driven the accumulation of > 500 Mb of knob repeats and affected the segregation of thousands of genes linked to knobs on all 10 chromosomes.

Fertilizer management for global ammonia emission reduction.

Crop production is a large source of atmospheric ammonia (NH3), which poses risks to air quality, human health and ecosystems1-5. However, estimating global NH3 emissions from croplands is subject to uncertainties because of data limitations, thereby limiting the accurate identification of mitigation options and efficacy4,5. Here we develop a machine learning model for generating crop-specific and spatially explicit NH3 emission factors globally (5-arcmin resolution) based on a compiled dataset of field observations. We show that global NH3 emissions from rice, wheat and maize fields in 2018 were 4.3 ± 1.0 Tg N yr-1, lower than previous estimates that did not fully consider fertilizer management practices6-9. Furthermore, spatially optimizing fertilizer management, as guided by the machine learning model, has the potential to reduce the NH3 emissions by about 38% (1.6 ± 0.4 Tg N yr-1) without altering total fertilizer nitrogen inputs. Specifically, we estimate potential NH3 emissions reductions of 47% (44-56%) for rice, 27% (24-28%) for maize and 26% (20-28%) for wheat cultivation, respectively. Under future climate change scenarios, we estimate that NH3 emissions could increase by 4.0 ± 2.7% under SSP1-2.6 and 5.5 ± 5.7% under SSP5-8.5 by 2030-2060. However, targeted fertilizer management has the potential to mitigate these increases.

Distinct kinesin motors drive two types of maize neocentromeres.

A maize chromosome variant called abnormal chromosome 10 (Ab10) converts knobs on chromosome arms into neocentromeres, causing their preferential segregation to egg cells in a process known as meiotic drive. We previously demonstrated that the gene Kinesin driver (Kindr) on Ab10 encodes a kinesin-14 required to mobilize neocentromeres made up of the major tandem repeat knob180. Here we describe a second kinesin-14 gene, TR-1 kinesin (Trkin), that is required to mobilize neocentromeres made up of the minor tandem repeat TR-1. Trkin lies in a 4-Mb region of Ab10 that is not syntenic with any other region of the maize genome and shows extraordinary sequence divergence from Kindr and other kinesins in plants. Despite its unusual structure, Trkin encodes a functional minus end-directed kinesin that specifically colocalizes with TR-1 in meiosis, forming long drawn out neocentromeres. TRKIN contains a nuclear localization signal and localizes to knobs earlier in prophase than KINDR. The fact that TR-1 repeats often co-occur with knob180 repeats suggests that the current role of the TRKIN/TR-1 system is to facilitate the meiotic drive of the KINDR/knob180 system.

Molecular mechanisms underlying gene regulatory variation of maize metabolic traits.

Variation in gene expression levels is pervasive among individuals and races or varieties, and has substantial agronomic consequences, for example, by contributing to hybrid vigor. Gene expression level variation results from mutations in regulatory sequences (cis) and/or transcription factor (TF) activity (trans), but the mechanisms underlying cis- and/or trans-regulatory variation of complex phenotypes remain largely unknown. Here, we investigated gene expression variation mechanisms underlying the differential accumulation of the insecticidal compounds maysin and chlorogenic acid in silks of widely used maize (Zea mays) inbreds, B73 and A632. By combining transcriptomics and cistromics, we identified 1,338 silk direct targets of the maize R2R3-MYB TF Pericarp color1 (P1), consistent with it being a regulator of maysin and chlorogenic acid biosynthesis. Among these P1 targets, 464 showed allele-specific expression (ASE) between B73 and A632 silks. Allelic DNA-affinity purification sequencing identified 34 examples in which P1 allelic specific binding (ASB) correlated with cis-expression variation. From previous yeast one-hybrid studies, we identified 9 TFs potentially implicated in the control of P1 targets, with ASB to 83 out of 464 ASE genes (cis) and differential expression of 4 out of 9 TFs between B73 and A632 silks (trans). These results provide a molecular framework for understanding universal mechanisms underlying natural variation of gene expression levels, and how the regulation of metabolic diversity is established.

The heat shock factor 20-HSF4-cellulose synthase A2 module regulates heat stress tolerance in maize.

Temperature shapes the geographical distribution and behavior of plants. Understanding the regulatory mechanisms underlying the plant heat stress response is important for developing climate-resilient crops, including maize (Zea mays). To identify transcription factors (TFs) that may contribute to the maize heat stress response, we generated a dataset of short- and long-term transcriptome changes following a heat treatment time course in the inbred line B73. Co-expression network analysis highlighted several TFs, including the class B2a heat shock factor (HSF) ZmHSF20. Zmhsf20 mutant seedlings exhibited enhanced tolerance to heat stress. Furthermore, DNA affinity purification sequencing and Cleavage Under Targets and Tagmentation assays demonstrated that ZmHSF20 binds to the promoters of Cellulose synthase A2 (ZmCesA2) and three class A Hsf genes, including ZmHsf4, repressing their transcription. We showed that ZmCesA2 and ZmHSF4 promote the heat stress response, with ZmHSF4 directly activating ZmCesA2 transcription. In agreement with the transcriptome analysis, ZmHSF20 inhibited cellulose accumulation and repressed the expression of cell wall-related genes. Importantly, the Zmhsf20 Zmhsf4 double mutant exhibited decreased thermotolerance, placing ZmHsf4 downstream of ZmHsf20. We proposed an expanded model of the heat stress response in maize, whereby ZmHSF20 lowers seedling heat tolerance by repressing ZmHsf4 and ZmCesA2, thus balancing seedling growth and defense.

NIN-LIKE PROTEIN3.2 inhibits repressor Aux/IAA14 expression and enhances root biomass in maize seedlings under low nitrogen.

Plants generally enhance their root growth in the form of greater biomass and/or root length to boost nutrient uptake in response to short-term low nitrogen (LN). However, the underlying mechanisms of short-term LN-mediated root growth remain largely elusive. Our genome-wide association study, haplotype analysis, and phenotyping of transgenic plants showed that the crucial nitrate signaling component NIN-LIKE PROTEIN3.2 (ZmNLP3.2), a positive regulator of root biomass, is associated with natural variations in root biomass of maize (Zea mays L.) seedlings under LN. The monocot-specific gene AUXIN/INDOLE-3-ACETIC ACID14 (ZmAux/IAA14) exhibited opposite expression patterns to ZmNLP3.2 in ZmNLP3.2 knockout and overexpression lines, suggesting that ZmNLP3.2 hampers ZmAux/IAA14 transcription. Importantly, ZmAux/IAA14 knockout seedlings showed a greater root dry weight (RDW), whereas ZmAux/IAA14 overexpression reduced RDW under LN compared with wild-type plants, indicating that ZmAux/IAA14 negatively regulates the RDW of LN-grown seedlings. Moreover, in vitro and vivo assays indicated that AUXIN RESPONSE FACTOR19 (ZmARF19) binds to and transcriptionally activates ZmAux/IAA14, which was weakened by the ZmNLP3.2-ZmARF19 interaction. The zmnlp3.2 ZmAux/IAA14-OE seedlings exhibited further reduced RDW compared with ZmAux/IAA14 overexpression lines when subjected to LN treatment, corroborating the ZmNLP3.2-ZmAux/IAA14 interaction. Thus, our study reveals a ZmNLP3.2-ZmARF19-ZmAux/IAA14 module regulating root biomass in response to nitrogen limitation in maize.

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.

The RALF signaling pathway regulates cell wall integrity during pollen tube growth in maize.

Autocrine signaling pathways regulated by RAPID ALKALINIZATION FACTORs (RALFs) control cell wall integrity during pollen tube germination and growth in Arabidopsis (Arabidopsis thaliana). To investigate the role of pollen-specific RALFs in another plant species, we combined gene expression data with phylogenetic and biochemical studies to identify candidate orthologs in maize (Zea mays). We show that Clade IB ZmRALF2/3 mutations, but not Clade III ZmRALF1/5 mutations, cause cell wall instability in the sub-apical region of the growing pollen tube. ZmRALF2/3 are mainly located in the cell wall and are partially able to complement the pollen germination defect of their Arabidopsis orthologs AtRALF4/19. Mutations in ZmRALF2/3 compromise pectin distribution patterns leading to altered cell wall organization and thickness culminating in pollen tube burst. Clade IB, but not Clade III ZmRALFs, strongly interact as ligands with the pollen-specific Catharanthus roseus RLK1-like (CrRLK1L) receptor kinases Z. mays FERONIA-like (ZmFERL) 4/7/9, LORELEI-like glycosylphosphatidylinositol-anchor (LLG) proteins Z. mays LLG 1 and 2 (ZmLLG1/2), and Z. mays pollen extension-like (PEX) cell wall proteins ZmPEX2/4. Notably, ZmFERL4 outcompetes ZmLLG2 and ZmPEX2 outcompetes ZmFERL4 for ZmRALF2 binding. Based on these data, we suggest that Clade IB RALFs act in a dual role as cell wall components and extracellular sensors to regulate cell wall integrity and thickness during pollen tube growth in maize and probably other plants.

Dual binding of chromomethylase domains to H3K9me2-containing nucleosomes directs DNA methylation in plants.

DNA methylation and histone modification exert epigenetic control over gene expression. CHG methylation by CHROMOMETHYLASE3 (CMT3) depends on histone H3K9 dimethylation (H3K9me2), but the mechanism underlying this relationship is poorly understood. Here, we report multiple lines of evidence that CMT3 interacts with H3K9me2-containing nucleosomes. CMT3 genome locations nearly perfectly correlated with H3K9me2, and CMT3 stably associated with H3K9me2-containing nucleosomes. Crystal structures of maize CMT3 homolog ZMET2, in complex with H3K9me2 peptides, showed that ZMET2 binds H3K9me2 via both bromo adjacent homology (BAH) and chromo domains. The structures reveal an aromatic cage within both BAH and chromo domains as interaction interfaces that capture H3K9me2. Mutations that abolish either interaction disrupt CMT3 binding to nucleosomes and show a complete loss of CMT3 activity in vivo. Our study establishes dual recognition of H3K9me2 marks by BAH and chromo domains and reveals a distinct mechanism of interplay between DNA methylation and histone modification.

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.

Premeiotic 24-nt phasiRNAs are present in the Zea genus and unique in biogenesis mechanism and molecular function.

Reproductive phasiRNAs (phased, small interfering RNAs) are broadly present in angiosperms and play crucial roles in sustaining male fertility. While the premeiotic 21-nt (nucleotides) phasiRNAs and meiotic 24-nt phasiRNA pathways have been extensively studied in maize (Zea mays) and rice (Oryza sativa), a third putative category of reproductive phasiRNAs-named premeiotic 24-nt phasiRNAs-have recently been reported in barley (Hordeum vulgare) and wheat (Triticum aestivum). To determine whether premeiotic 24-nt phasiRNAs are also present in maize and related species and begin to characterize their biogenesis and function, we performed a comparative transcriptome and degradome analysis of premeiotic and meiotic anthers from five maize inbred lines and three teosinte species/subspecies. Our data indicate that a substantial subset of the 24-nt phasiRNA loci in maize and teosinte are already highly expressed at the premeiotic phase. The premeiotic 24-nt phasiRNAs are similar to meiotic 24-nt phasiRNAs in genomic origin and dependence on DCL5 (Dicer-like 5) for biogenesis, however, premeiotic 24-nt phasiRNAs are unique in that they are likely i) not triggered by microRNAs, ii) not loaded by AGO18 proteins, and iii) not capable of mediating PHAS precursor cleavage. In addition, we also observed a group of premeiotic 24-nt phasiRNAs in rice using previously published data. Together, our results indicate that the premeiotic 24-nt phasiRNAs constitute a unique class of reproductive phasiRNAs and are present more broadly in the grass family (Poaceae) than previously known.

Root-exuded specialized metabolites reduce arsenic toxicity in maize.

By releasing specialized metabolites, plants modify their environment. Whether and how specialized metabolites protect plants against toxic levels of trace elements is not well understood. We evaluated whether benzoxazinoids, which are released into the soil by major cereals, can confer protection against arsenic toxicity. Benzoxazinoid-producing maize plants performed better in arsenic-contaminated soils than benzoxazinoid-deficient mutants in the greenhouse and the field. Adding benzoxazinoids to the soil restored the protective effect, and the effect persisted to the next crop generation via positive plant-soil feedback. Arsenate levels in the soil and total arsenic levels in the roots were lower in the presence of benzoxazinoids. Thus, the protective effect of benzoxazinoids is likely soil-mediated and includes changes in soil arsenic speciation and root accumulation. We conclude that exuded specialized metabolites can enhance protection against toxic trace elements via soil-mediated processes and may thereby stabilize crop productivity in polluted agroecosystems.

Investigating the cis-regulatory basis of C3 and C4 photosynthesis in grasses at single-cell resolution.

While considerable knowledge exists about the enzymes pivotal for C4 photosynthesis, much less is known about the cis-regulation important for specifying their expression in distinct cell types. Here, we use single-cell-indexed ATAC-seq to identify cell-type-specific accessible chromatin regions (ACRs) associated with C4 enzymes for five different grass species. This study spans four C4 species, covering three distinct photosynthetic subtypes: Zea mays and Sorghum bicolor (NADP-dependent malic enzyme), Panicum miliaceum (NAD-dependent malic enzyme), Urochloa fusca (phosphoenolpyruvate carboxykinase), along with the C3 outgroup Oryza sativa. We studied the cis-regulatory landscape of enzymes essential across all C4 species and those unique to C4 subtypes, measuring cell-type-specific biases for C4 enzymes using chromatin accessibility data. Integrating these data with phylogenetics revealed diverse co-option of gene family members between species, showcasing the various paths of C4 evolution. Besides promoter proximal ACRs, we found that, on average, C4 genes have two to three distal cell-type-specific ACRs, highlighting the complexity and divergent nature of C4 evolution. Examining the evolutionary history of these cell-type-specific ACRs revealed a spectrum of conserved and novel ACRs, even among closely related species, indicating ongoing evolution of cis-regulation at these C4 loci. This study illuminates the dynamic and complex nature of cis-regulatory elements evolution in C4 photosynthesis, particularly highlighting the intricate cis-regulatory evolution of key loci. Our findings offer a valuable resource for future investigations, potentially aiding in the optimization of C3 crop performance under changing climatic conditions.

Mesophyll airspace unsaturation drives C4 plant success under vapor pressure deficit stress.

A fundamental assumption in plant science posits that leaf air spaces remain vapor saturated, leading to the predominant view that stomata alone control leaf water loss. This concept has been pivotal in photosynthesis and water-use efficiency research. However, recent evidence has refuted this longstanding assumption by providing evidence of unsaturation in the leaf air space of C3 plants under relatively mild vapor pressure deficit (VPD) stress. This phenomenon represents a nonstomatal mechanism restricting water loss from the mesophyll. The potential ubiquity and physiological implications of this phenomenon, its driving mechanisms in different plant species and habitats, and its interaction with other ecological adaptations have. In this context, C4 plants spark particular interest for their importance as crops, bundle sheath cells' unique anatomical characteristics and specialized functions, and notably higher water-use efficiency relative to C3 plants. Here, we confirm reduced relative humidities in the substomatal cavity of the C4 plants maize, sorghum, and proso millet down to 80% under mild VPD stress. We demonstrate the critical role of nonstomatal control in these plants, indicating that the role of the CO2 concentration mechanism in CO2 management at a high VPD may have been overestimated. Our findings offer a mechanistic reconciliation between discrepancies in CO2 and VPD responses reported in C4 species. They also reveal that nonstomatal control is integral to maintaining an advantageous microclimate of relatively higher CO2 concentrations in the mesophyll air space of C4 plants for carbon fixation, proving vital when these plants face VPD stress.

Membralin is required for maize development and defines a branch of the endoplasmic reticulum-associated degradation pathway in plants.

Endoplasmic reticulum (ER)-associated degradation (ERAD) plays key roles in controlling protein levels and quality in eukaryotes. The Ring Finger Protein 185 (RNF185)/membralin ubiquitin ligase complex was recently identified as a branch in mammals and is essential for neuronal function, but its function in plant development is unknown. Here, we report the map-based cloning and characterization of Narrow Leaf and Dwarfism 1 (NLD1), which encodes the ER membrane-localized protein membralin and specifically interacts with maize homologs of RNF185 and related components. The nld1 mutant shows defective leaf and root development due to reduced cell number. The defects of nld1 were largely restored by expressing membralin genes from Arabidopsis thaliana and mice, highlighting the conserved roles of membralin proteins in animals and plants. The excessive accumulation of ß-hydroxy ß-methylglutaryl-CoA reductase in nld1 indicates that the enzyme is a membralin-mediated ERAD target. The activation of bZIP60 mRNA splicing-related unfolded protein response signaling and marker gene expression in nld1, as well as DNA fragment and cell viability assays, indicate that membralin deficiency induces ER stress and cell death in maize, thereby affecting organogenesis. Our findings uncover the conserved, indispensable role of the membralin-mediated branch of the ERAD pathway in plants. In addition, ZmNLD1 contributes to plant architecture in a dose-dependent manner, which can serve as a potential target for genetic engineering to shape ideal plant architecture, thereby enhancing high-density maize yields.

De novo TANGLED1 recruitment from the phragmoplast to aberrant cell plate fusion sites in maize.

Division plane positioning is crucial for proper growth and development in many organisms. In plants, the division plane is established before mitosis, by accumulation of a cytoskeletal structure called the preprophase band (PPB). The PPB is thought to be essential for recruitment of division site-localized proteins, which remain at the division site after the PPB disassembles. Here, we show that the division site-localized protein TANGLED1 (TAN1) is recruited independently of the PPB to the cell cortex by the plant cytokinetic machinery, the phragmoplast, from experiments using both the PPB-defective mutant discordia1 (dcd1) and chemical treatments that disrupt the phragmoplast in maize. TAN1 recruitment to de novo sites on the cortex is partially dependent on intact actin filaments and the myosin XI motor protein OPAQUE1 (O1). These data imply a yet unknown role for TAN1 and possibly other division site-localized proteins during the last stages of cell division when the phragmoplast touches the cell cortex to complete cytokinesis.

UB2/UB3/TSH4-anchored transcriptional networks regulate early maize inflorescence development in response to simulated shade.

Increasing planting density has been adopted as an effective means to increase maize (Zea mays) yield. Competition for light from neighbors can trigger plant shade avoidance syndrome, which includes accelerated flowering. However, the regulatory networks of maize inflorescence development in response to high-density planting remain poorly understood. In this study, we showed that shade-mimicking treatments cause precocious development of the tassels and ears. Comparative transcriptome profiling analyses revealed the enrichment of phytohormone-related genes and transcriptional regulators among the genes co-regulated by developmental progression and simulated shade. Network analysis showed that three homologous Squamosa promoter binding protein (SBP)-like (SPL) transcription factors, Unbranched2 (UB2), Unbranched3 (UB3), and Tasselsheath4 (TSH4), individually exhibited connectivity to over 2,400 genes across the V3-to-V9 stages of tassel development. In addition, we showed that the ub2 ub3 double mutant and tsh4 single mutant were almost insensitive to simulated shade treatments. Moreover, we demonstrated that UB2/UB3/TSH4 could directly regulate the expression of Barren inflorescence2 (BIF2) and Zea mays teosinte branched1/cycloidea/proliferating cell factor30 (ZmTCP30). Furthermore, we functionally verified a role of ZmTCP30 in regulating tassel branching and ear development. Our results reveal a UB2/UB3/TSH4-anchored transcriptional regulatory network of maize inflorescence development and provide valuable targets for breeding shade-tolerant maize cultivars.

Maize PPR-E proteins mediate RNA C-to-U editing in mitochondria by recruiting the trans deaminase PCW1.

RNA C-to-U editing in organelles is essential for plant growth and development; however, the underlying mechanism is not fully understood. Here, we report that pentatricopeptide repeat (PPR)-E subclass proteins carry out RNA C-to-U editing by recruiting the trans deaminase PPR motifs, coiled-coil, and DYW domain-containing protein 1 (PCW1) in maize (Zea mays) mitochondria. Loss-of-function of bZIP and coiled-coil domain-containing PPR 1 (bCCP1) or PCW1 arrests seed development in maize. bCCP1 encodes a bZIP and coiled-coil domain-containing PPR protein, and PCW1 encodes an atypical PPR-DYW protein. bCCP1 is required for editing at 66 sites in mitochondria and PCW1 is required for editing at 102 sites, including the 66 sites that require bCCP1. The PCW1-mediated editing sites are exclusively associated with PPR-E proteins. bCCP1 interacts with PCW1 and the PPR-E protein Empty pericarp7 (EMP7). Two multiple organellar RNA editing factor (MORF) proteins, ZmMORF1 and ZmMORF8, interact with PCW1, EMP7, and bCCP1. ZmMORF8 enhanced the EMP7-PCW1 interaction in a yeast three-hybrid assay. C-to-U editing at the ccmFN-1553 site in maize required EMP7, bCCP1, and PCW1. These results suggest that PPR-E proteins function in RNA editing by recruiting the trans deaminase PCW1 and bCCP1, and MORF1/8 assist this recruitment through protein-protein interactions.

The transcription factors ZmNAC128 and ZmNAC130 coordinate with Opaque2 to promote endosperm filling in maize.

Endosperm filling in maize (Zea mays), which involves nutrient uptake and biosynthesis of storage reserves, largely determines grain yield and quality. However, much remains unclear about the synchronization of these processes. Here, we comprehensively investigated the functions of duplicate NAM, ATAF1/2, and CUC2 (NAC)-type transcription factors, namely, ZmNAC128 and ZmNAC130, in endosperm filling. The gene-edited double mutant zmnac128 zmnac130 exhibits a poorly filled kernel phenotype such that the kernels have an inner cavity. RNA sequencing and protein abundance analysis revealed that the expression of many genes involved in the biosynthesis of zein and starch is reduced in the filling endosperm of zmnac128 zmnac130. Further, DNA affinity purification and sequencing combined with chromatin-immunoprecipitation quantitative PCR and promoter transactivation assays demonstrated that ZmNAC128 and ZmNAC130 are direct regulators of 3 (16-, 27-, and 50-kD) γ-zein genes and 6 important starch metabolism genes (Brittle2 [Bt2], pullulanase-type starch debranching enzyme [Zpu1], granule-bound starch synthase 1 [GBSS1], starch synthase 1 [SS1], starch synthase IIa [SSIIa], and sucrose synthase 1 [Sus1]). ZmNAC128 and ZmNAC130 recognize an additional cis-element in the Opaque2 (O2) promoter to regulate its expression. The triple mutant zmnac128 zmnac130 o2 exhibits extremely poor endosperm filling, which results in more than 70% of kernel weight loss. ZmNAC128 and ZmNAC130 regulate the expression of the transporter genes sugars that will eventually be exported transporter 4c (ZmSWEET4c), sucrose and glucose carrier 1 (ZmSUGCAR1), and yellow stripe-like2 (ZmYSL2) and in turn facilitate nutrient uptake, while O2 plays a supporting role. In conclusion, ZmNAC128 and ZmNAC130 cooperate with O2 to facilitate endosperm filling, which involves nutrient uptake in the basal endosperm transfer layer (BETL) and the synthesis of zeins and starch in the starchy endosperm (SE).

NAKED ENDOSPERM1, NAKED ENDOSPERM2, and OPAQUE2 interact to regulate gene networks in maize endosperm development.

NAKED ENDOSPERM1 (NKD1), NKD2, and OPAQUE2 (O2) are transcription factors important for cell patterning and nutrient storage in maize (Zea mays) endosperm. To study the complex regulatory interrelationships among these 3 factors in coregulating gene networks, we developed a set of nkd1, nkd2, and o2 homozygous lines, including all combinations of mutant and wild-type genes. Among the 8 genotypes tested, we observed diverse phenotypes and gene interactions affecting cell patterning, starch content, and storage proteins. From ∼8 to ∼16 d after pollination, maize endosperm undergoes a transition from cellular development to nutrient accumulation for grain filling. Gene network analysis showed that NKD1, NKD2, and O2 dynamically regulate a hierarchical gene network during this period, directing cellular development early and then transitioning to constrain cellular development while promoting the biosynthesis and storage of starch, proteins, and lipids. Genetic interactions regulating this network are also dynamic. The assay for transposase-accessible chromatin using sequencing (ATAC-seq) showed that O2 influences the global regulatory landscape, decreasing NKD1 and NKD2 target site accessibility, while NKD1 and NKD2 increase O2 target site accessibility. In summary, interactions of NKD1, NKD2, and O2 dynamically affect the hierarchical gene network and regulatory landscape during the transition from cellular development to grain filling in maize endosperm.