A cis-regulatory atlas in maize at single-cell resolution.

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

A Nucleosome Bridging Mechanism for Activation of a Maintenance DNA Methyltransferase.

DNA methylation and H3K9me are hallmarks of heterochromatin in plants and mammals, and are successfully maintained across generations. The biochemical and structural basis for this maintenance is poorly understood. The maintenance DNA methyltransferase from Zea mays, ZMET2, recognizes dimethylation of H3K9 via a chromodomain (CD) and a bromo adjacent homology (BAH) domain, which flank the catalytic domain. Here, we show that dinucleosomes are the preferred ZMET2 substrate, with DNA methylation preferentially targeted to linker DNA. Electron microscopy shows one ZMET2 molecule bridging two nucleosomes within a dinucleosome. We find that the CD stabilizes binding, whereas the BAH domain enables allosteric activation by the H3K9me mark. ZMET2 further couples recognition of H3K9me to an increase in the specificity for hemimethylated versus unmethylated DNA. We propose a model in which synergistic coupling between recognition of nucleosome spacing, H3K9 methylation, and DNA modification allows ZMET2 to maintain DNA methylation in heterochromatin with high fidelity.

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.

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.

The maize preligule band is subdivided into distinct domains with contrasting cellular properties prior to ligule outgrowth.

The maize ligule is an epidermis-derived structure that arises from the preligule band (PLB) at a boundary between the blade and sheath. A hinge-like auricle also develops immediately distal to the ligule and contributes to blade angle. Here, we characterize the stages of PLB and early ligule development in terms of topography, cell area, division orientation, cell wall rigidity and auxin response dynamics. Differential thickening of epidermal cells and localized periclinal divisions contributed to the formation of a ridge within the PLB, which ultimately produces the ligule fringe. Patterns in cell wall rigidity were consistent with the subdivision of the PLB into two regions along a distinct line positioned at the nascent ridge. The proximal region produces the ligule, while the distal region contributes to one epidermal face of the auricles. Although the auxin transporter PIN1 accumulated in the PLB, observed differential auxin transcriptional response did not underlie the partitioning of the PLB. Our data demonstrate that two zones with contrasting cellular properties, the preligule and preauricle, are specified within the ligular region before ligule outgrowth.

Transcription factor bHLH121 regulates root cortical aerenchyma formation in maize.

Root anatomical phenotypes present a promising yet underexploited avenue to deliver major improvements in yield and climate resilience of crops by improving water and nutrient uptake. For instance, the formation of root cortical aerenchyma (RCA) significantly increases soil exploration and resource capture by reducing the metabolic costs of root tissue. A key bottleneck in studying such phenotypes has been the lack of robust high-throughput anatomical phenotyping platforms. We exploited a phenotyping approach based on laser ablation tomography, termed Anatomics, to quantify variation in RCA formation of 436 diverse maize lines in the field. Results revealed a significant and heritable variation for RCA formation. Genome-wide association studies identified a single-nucleotide polymorphism mapping to a root cortex-expressed gene-encoding transcription factor bHLH121. Functional studies identified that the bHLH121 Mu transposon mutant line and CRISPR/Cas9 loss-of-function mutant line showed reduced RCA formation, whereas an overexpression line exhibited significantly greater RCA formation when compared to the wild-type line. Characterization of these lines under suboptimal water and nitrogen availability in multiple soil environments revealed that bHLH121 is required for RCA formation developmentally as well as under studied abiotic stress. Overall functional validation of the bHLH121 gene's importance in RCA formation provides a functional marker to select varieties with improved soil exploration and thus yield under suboptimal conditions.

Cross-species predictive modeling reveals conserved drought responses between maize and sorghum.

Drought tolerance is a highly complex trait controlled by numerous interconnected pathways with substantial variation within and across plant species. This complexity makes it difficult to distill individual genetic loci underlying tolerance, and to identify core or conserved drought-responsive pathways. Here, we collected drought physiology and gene expression datasets across diverse genotypes of the C4 cereals sorghum and maize and searched for signatures defining water-deficit responses. Differential gene expression identified few overlapping drought-associated genes across sorghum genotypes, but using a predictive modeling approach, we found a shared core drought response across development, genotype, and stress severity. Our model had similar robustness when applied to datasets in maize, reflecting a conserved drought response between sorghum and maize. The top predictors are enriched in functions associated with various abiotic stress-responsive pathways as well as core cellular functions. These conserved drought response genes were less likely to contain deleterious mutations than other gene sets, suggesting that core drought-responsive genes are under evolutionary and functional constraints. Our findings support a broad evolutionary conservation of drought responses in C4 grasses regardless of innate stress tolerance, which could have important implications for developing climate resilient cereals.

Bacterial tolerance to host-exuded specialized metabolites structures the maize root microbiome.

Plants exude specialized metabolites from their roots, and these compounds are known to structure the root microbiome. However, the underlying mechanisms are poorly understood. We established a representative collection of maize root bacteria and tested their tolerance against benzoxazinoids (BXs), the dominant specialized and bioactive metabolites in the root exudates of maize plants. In vitro experiments revealed that BXs inhibited bacterial growth in a strain- and compound-dependent manner. Tolerance against these selective antimicrobial compounds depended on bacterial cell wall structure. Further, we found that native root bacteria isolated from maize tolerated the BXs better compared to nonhost Arabidopsis bacteria. This finding suggests the adaptation of the root bacteria to the specialized metabolites of their host plant. Bacterial tolerance to 6-methoxy-benzoxazolin-2-one (MBOA), the most abundant and selective antimicrobial metabolite in the maize rhizosphere, correlated significantly with the abundance of these bacteria on BX-exuding maize roots. Thus, strain-dependent tolerance to BXs largely explained the abundance pattern of bacteria on maize roots. Abundant bacteria generally tolerated MBOA, while low abundant root microbiome members were sensitive to this compound. Our findings reveal that tolerance to plant specialized metabolites is an important competence determinant for root colonization. We propose that bacterial tolerance to root-derived antimicrobial compounds is an underlying mechanism determining the structure of host-specific microbial communities.

ZmTDM1 encodes a tetratricopeptide repeat domain protein and is required for meiotic exit in maize.

Elaborate cell-cycle control must be adopted to ensure the continuity of the meiotic second division and termination after that. Despite its importance, however, the genetic controls underlying the meiotic cell cycle have not been reported in maize. Here, we characterized a meiotic cell-cycle controller ZmTDM1, which is a homolog of Arabidopsis TDM1 and encodes a canonical tetratricopeptide repeat domain protein in maize. The Zmtdm1 homozygous plants exhibited complete male sterility and severe female abortion. In Zmtdm1 mutants, cell-cycle progression was almost identical to that of wild type from leptotene to anaphase II. However, chromosomes in the tetrad failed meiotic termination at the end of the second division and underwent additional divisions in succession without DNA replication, reducing the ploidy to less than haploid in the product. In addition, two ZmTDM1-like homologs (ZmTDML1 and ZmTDML2) were not functional in meiotic cell-cycle control. Moreover, ZmTDM1 interacted with RING-type E3 ubiquitin ligase, revealing that it acts as a subunit of the APC/C E3 ubiquitin ligase complex. Overall, our results identified a regulator of meiotic cell cycle in maize and demonstrated that ZmTDM1 is essential for meiotic exit after meiosis II.

GOLDEN2-like1 is sufficient but not necessary for chloroplast biogenesis in mesophyll cells of C4 grasses.

Chloroplasts are the site of photosynthesis. In land plants, chloroplast biogenesis is regulated by a family of transcription factors named GOLDEN2-like (GLK). In C4 grasses, it has been hypothesized that genome duplication events led to the sub-functionalization of GLK paralogs (GLK1 and GLK2) to control chloroplast biogenesis in two distinct cell types: mesophyll and bundle sheath cells. Although previous characterization of golden2 (g2) mutants in maize has demonstrated a role for GLK2 paralogs in regulating chloroplast biogenesis in bundle sheath cells, the function of GLK1 has remained elusive. Here we show that, contrary to expectations, GLK1 is not required for chloroplast biogenesis in mesophyll cells of maize. Comparisons between maize and Setaria viridis, which represent two independent C4 origins within the Poales, further show that the role of GLK paralogs in controlling chloroplast biogenesis in mesophyll and bundle sheath cells differs between species. Despite these differences, complementation analysis revealed that GLK1 and GLK2 genes from maize are both sufficient to restore functional chloroplast development in mesophyll and bundle sheath cells of S. viridis mutants. Collectively our results suggest an evolutionary trajectory in C4 grasses whereby both orthologs retained the ability to induce chloroplast biogenesis but GLK2 adopted a more prominent developmental role, particularly in relation to chloroplast activation in bundle sheath cells.

Widespread changes to the translational landscape in a maize microRNA biogenesis mutant.

MicroRNAs are short, non-coding RNAs that repress gene expression in both plants and animals and have diverse functions related to growth, development, and stress responses. The ribonuclease, DICER-LIKE1 (DCL1) is required for two steps in plant miRNA biogenesis: cleavage of the primary miRNAs (pri-miRNAs) to release a hairpin structure, called the precursor miRNA (pre-miRNA) and cleavage of the pre-miRNA to generate the miRNA/miRNA* duplex. The mature miRNA guides the RNA-induced silencing complex to target RNAs with complementary sequences, resulting in translational repression and/or RNA cleavage of target mRNAs. However, the relative contribution of translational repression versus mRNA degradation by miRNAs remains unknown at the genome-level in crops, especially in maize. The maize fuzzy tassel (fzt) mutant contains a hypomorphic mutation in DCL1 resulting in broad developmental defects. While most miRNAs are reduced in fzt, the levels of miRNA-targeted mRNAs are not dramatically increased, suggesting that translational regulation by miRNAs may be common. To gain insight into the repression mechanism of plant miRNAs, we combined ribosome profiling and RNA-sequencing to globally survey miRNA activities in maize. Our data indicate that translational repression contributes significantly to regulation of most miRNA targets and that approximately one-third of miRNA targets are regulated primarily at the translational level. Surprisingly, ribosomes appear altered in fzt mutants suggesting that DCL1 may also have a role in ribosome biogenesis. Thus, DICER-LIKE1 shapes the translational landscape in plants through both miRNA-dependent and miRNA-independent mechanisms.

Molecular mechanisms underlying low temperature inhibition of grain filling in maize (Zea mays L.): coordination of growth and cold responses.

Low temperature (LT) greatly restricts grain filling in maize (Zea mays L.), but the relevant molecular mechanisms are not fully understood. To better understand the effect of LT on grain development, 17 hybrids were subjected to LT stress in field trials over 3 years, and two hybrids of them with contrasting LT responses were exposed to 30/20°C and 20/10°C for 7 days during grain filling in a greenhouse. At LT, thousand-kernel weight declined, especially in LT-sensitive hybrid FM985, while grain-filling rate was on average about 48% higher in LT-tolerant hybrid DK159 than FM985. LT reduced starch synthesis in kernel mainly by suppression of transcript levels and enzyme activities for sucrose synthase and hexokinase. Brassinolide (BR) was abundant in DK159 kernel, and genes involved in BR and cytokinin signals were inducible by stress. LT downregulated the genes in light-harvesting complex and photosystem I/II subunits, accompanied by reduced photosynthetic rate and Fv/Fm in ear leaf. The LT-tolerant hybrid could maintain a high soluble sugar content and fast interconversion between sucrose and hexose in the stem internode and cob, improving assimilate allocation to kernel at LT stress and paving the way for simultaneous growth and LT stress responses.

Population-level gene expression can repeatedly link genes to functions in maize.

Transcriptome-wide association studies (TWAS) can provide single gene resolution for candidate genes in plants, complementing genome-wide association studies (GWAS) but efforts in plants have been met with, at best, mixed success. We generated expression data from 693 maize genotypes, measured in a common field experiment, sampled over a 2-h period to minimize diurnal and environmental effects, using full-length RNA-seq to maximize the accurate estimation of transcript abundance. TWAS could identify roughly 10 times as many genes likely to play a role in flowering time regulation as GWAS conducted data from the same experiment. TWAS using mature leaf tissue identified known true-positive flowering time genes known to act in the shoot apical meristem, and trait data from a new environment enabled the identification of additional flowering time genes without the need for new expression data. eQTL analysis of TWAS-tagged genes identified at least one additional known maize flowering time gene through trans-eQTL interactions. Collectively these results suggest the gene expression resource described here can link genes to functions across different plant phenotypes expressed in a range of tissues and scored in different experiments.

Identification of anther thermotolerance genes by the integration of linkage and association analysis in maize.

Maize is one of the world's most important staple crops, yet its production is increasingly threatened by the rising frequency of high-temperature stress (HTS). To investigate the genetic basis of anther thermotolerance under field conditions, we performed linkage and association analysis to identify HTS response quantitative trait loci (QTL) using three recombinant inbred line (RIL) populations and an association panel containing 375 diverse maize inbred lines. These analyses resulted in the identification of 16 co-located large QTL intervals. Among the 37 candidate genes identified in these QTL intervals, five have rice or Arabidopsis homologs known to influence pollen and filament development. Notably, one of the candidate genes, ZmDUP707, has been subject to selection pressure during breeding. Its expression is suppressed by HTS, leading to pollen abortion and barren seeds. We also identified several additional candidate genes potentially underly QTL previously reported by other researchers. Taken together, our results provide a pool of valuable candidate genes that could be employed by future breeding programs aiming at enhancing maize HTS tolerance.

Toward unveiling transcriptome dynamics and regulatory modules at the maternal/filial interface of developing maize kernel.

The basal region of maize (Zea mays) kernels, which includes the pedicel, placenta-chalazal, and basal endosperm transfer layers, serves as the maternal/filial interface for nutrient transfer from the mother plant to the developing seed. However, transcriptome dynamics of this maternal/filial interface remain largely unexplored. To address this gap, we conducted high-temporal-resolution RNA sequencing of the basal and upper kernel regions between 4 and 32 days after pollination and deeply analyzed transcriptome dynamics of the maternal/filial interface. Utilizing 790 specifically and highly expressed genes in the basal region, we performed the gene ontology (GO) term and weighted gene co-expression network analyses. In the early-stage basal region, we identified five MADS-box transcription factors (TFs) as hubs. Their homologs have been demonstrated as pivotal regulators at the maternal/filial interface of rice or Arabidopsis, suggesting their potential roles in maize kernel development. In the filling-stage basal region, numerous GO terms associated with transcriptional regulation and transporters are significantly enriched. Furthermore, we investigated the molecular function of three hub TFs. Through genome-wide DNA affinity purification sequencing combined with promoter transactivation assays, we suggested that these three TFs act as regulators of 10 basal-specific transporter genes involved in the transfer of sugars, amino acids, and ions. This study provides insights into transcriptomic dynamic and regulatory modules of the maternal/filial interface. In the future, genetic investigation of these hub regulators must advance our understanding of maternal/filial interface development and function.

Improving resilience to high temperature in drought: water replenishment enhances sucrose and amino acid metabolisms in maize grain.

Heat stress poses a significant threat to maize, especially when combined with drought. Recent research highlights the potential of water replenishment to ameliorate grain weight loss. However, the mitigating mechanisms of heat in drought stress, especially during the crucial early grain-filling stage, remain poorly understood. We investigated the mechanism for mitigating heat in drought stress by water replenishment from the 12th to the 32nd days after silking in a controlled greenhouse experiment (Exp. I) and field trial (Exp. II). A significant reduction in grain weight was observed in heat stress compared to normal conditions. When water replenishment was applied to increase soil water content (SWC) under heat stress, the grain yield exhibited a notable increase ranging from 28.4 to 76.9%. XY335 variety was used for transcriptome sequencing to analyze starch biosynthesis and amino acid metabolisms in Exp. I. With water replenishment, the transcripts of genes responsible for trehalose 6-phosphate phosphates (TPP), alpha-trehalase (TRE), ADP-glcpyrophosphorylase, and starch synthase activity were stimulated. Additionally, the expression of genes encoding TPP and TRE contributed to an enhanced conversion of trehalose to glucose. This led to the conversion of sucrose from glucose-1-phosphate to ADP-glucose and ADP-glucose to amylopectin, ultimately increasing starch production by 45.1%. Water replenishment to boost SWC during heat stress also elevated the levels of essential amino acids in maize, including arginine, serine, tyrosine, leucine, glutamic acid, and methionine, providing valuable support to maize plants in adversity. Field trials further validated the positive impact of water replenishment on SWC, resulting in a notable increase in grain yield ranging from 7.1 to 9.2%. This study highlights the vital importance of adapting to abiotic stress and underscores the necessity of developing strategies to counteract its adverse effects on crop yield.

The additive function of YIGE2 and YIGE1 in regulating maize ear length.

Ear length (EL) is a key trait that greatly contributes to yield in maize. Although dozens of EL quantitative trait loci have been mapped, very few causal genes have been cloned, and the molecular mechanisms remain largely unknown. Our previous study showed that YIGE1 is involved in sugar and auxin pathways to regulate ear inflorescence meristem (IM) development and thus affects EL in maize. Here, we reveal that YIGE2, the paralog of YIGE1, regulates maize ear development and EL through auxin pathway. Knockout of YIGE2 causes a significant decrease of auxin level, IM length, floret number, EL, and grain yield. yige1 yige2 double mutants had even shorter IM and ears implying that these two genes redundantly regulate IM development and EL. The genes controlling auxin levels are differential expressed in yige1 yige2 double mutants, leading to lower auxin level. These results elucidated the critical role of YIGE2 and the redundancy between YIGE2 and YIGE1 in maize ear development, providing a new genetic resource for maize yield improvement.

Dynamic transcriptome landscape of maize pericarp development.

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