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.

Cistrome and Epicistrome Features Shape the Regulatory DNA Landscape.

The cistrome is the complete set of transcription factor (TF) binding sites (cis-elements) in an organism, while an epicistrome incorporates tissue-specific DNA chemical modifications and TF-specific chemical sensitivities into these binding profiles. Robust methods to construct comprehensive cistrome and epicistrome maps are critical for elucidating complex transcriptional networks that underlie growth, behavior, and disease. Here, we describe DNA affinity purification sequencing (DAP-seq), a high-throughput TF binding site discovery method that interrogates genomic DNA with in-vitro-expressed TFs. Using DAP-seq, we defined the Arabidopsis cistrome by resolving motifs and peaks for 529 TFs. Because genomic DNA used in DAP-seq retains 5-methylcytosines, we determined that >75% (248/327) of Arabidopsis TFs surveyed were methylation sensitive, a property that strongly impacts the epicistrome landscape. DAP-seq datasets also yielded insight into the biology and binding site architecture of numerous TFs, demonstrating the value of DAP-seq for cost-effective cistromic and epicistromic annotation in any organism.

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.

Genetic dissection of cis-regulatory control of ZmWUSCHEL1 expression by type B RESPONSE REGULATORS.

Mutations in cis-regulatory regions play an important role in the domestication and improvement of crops by altering gene expression. However, assessing the in vivo impact of cis-regulatory elements (CREs) on transcriptional regulation and phenotypic outcomes remains challenging. Previously, we showed that the dominant Barren inflorescence3 (Bif3) mutant of maize (Zea mays) contains a duplicated copy of the homeobox transcription factor gene ZmWUSCHEL1 (ZmWUS1), named ZmWUS1-B. ZmWUS1-B is controlled by a spontaneously generated novel promoter region that dramatically increases its expression and alters patterning and development of young ears. Overexpression of ZmWUS1-B is caused by a unique enhancer region containing multimerized binding sites for type B RESPONSE REGULATORs (RRs), key transcription factors in cytokinin signaling. To better understand how the enhancer increases the expression of ZmWUS1 in vivo, we specifically targeted the ZmWUS1-B enhancer region by CRISPR-Cas9-mediated editing. A series of deletion events with different numbers of type B RR DNA binding motifs (AGATAT) enabled us to determine how the number of AGATAT motifs impacts in vivo expression of ZmWUS1-B and consequently ear development. In combination with dual-luciferase assays in maize protoplasts, our analysis reveals that AGATAT motifs have an additive effect on ZmWUS1-B expression, while the distance separating AGATAT motifs does not appear to have a meaningful impact, indicating that the enhancer activity derives from the sum of individual CREs. These results also suggest that in maize inflorescence development, there is a threshold of buffering capacity for ZmWUS1 overexpression.

Necrotic upper tips1 mimics heat and drought stress and encodes a protoxylem-specific transcription factor in maize.

Maintaining sufficient water transport during flowering is essential for proper organ growth, fertilization, and yield. Water deficits that coincide with flowering result in leaf wilting, necrosis, tassel browning, and sterility, a stress condition known as "tassel blasting." We identified a mutant, necrotic upper tips1 (nut1), that mimics tassel blasting and drought stress and reveals the genetic mechanisms underlying these processes. The nut1 phenotype is evident only after the floral transition, and the mutants have difficulty moving water as shown by dye uptake and movement assays. These defects are correlated with reduced protoxylem vessel thickness that indirectly affects metaxylem cell wall integrity and function in the mutant. nut1 is caused by an Ac transposon insertion into the coding region of a unique NAC transcription factor within the VND clade of Arabidopsis NUT1 localizes to the developing protoxylem of root, stem, and leaf sheath, but not metaxylem, and its expression is induced by flowering. NUT1 downstream target genes function in cell wall biosynthesis, apoptosis, and maintenance of xylem cell wall thickness and strength. These results show that maintaining protoxylem vessel integrity during periods of high water movement requires the expression of specialized, dynamically regulated transcription factors within the vasculature.

NEEDLE1 encodes a mitochondria localized ATP-dependent metalloprotease required for thermotolerant maize growth.

Meristems are highly regulated structures ultimately responsible for the formation of branches, lateral organs, and stems, and thus directly affect plant architecture and crop yield. In meristems, genetic networks, hormones, and signaling molecules are tightly integrated to establish robust systems that can adapt growth to continuous inputs from the environment. Here we characterized needle1 (ndl1), a temperature-sensitive maize mutant that displays severe reproductive defects and strong genetic interactions with known mutants affected in the regulation of the plant hormone auxin. NDL1 encodes a mitochondria-localized ATP-dependent metalloprotease belonging to the FILAMENTATION TEMPERATURE-SENSITIVE H (FTSH) family. Together with the hyperaccumulation of reactive oxygen species (ROS), ndl1 inflorescences show up-regulation of a plethora of stress-response genes. We provide evidence that these conditions alter endogenous auxin levels and disrupt primordia initiation in meristems. These findings connect meristem redox status and auxin in the control of maize growth.

OsFD4 promotes the rice floral transition via florigen activation complex formation in the shoot apical meristem.

In rice, the florigens Heading Date 3a (Hd3a) and Rice Flowering Locus T 1 (RFT1), OsFD-like basic leucine zipper (bZIP) transcription factors, and Gf14 proteins assemble into florigen activation/repressor complexes (FACs/FRCs), which regulate transition to flowering in leaves and apical meristem. Only OsFD1 has been described as part of complexes promoting flowering at the meristem, and little is known about the role of other bZIP transcription factors, the combinatorial complexity of FAC formation, and their DNA-binding properties. Here, we used mutant analysis, protein-protein interaction assays and DNA affinity purification (DAP) sequencing coupled to in silico prediction of binding syntaxes to study several bZIP proteins that assemble into FACs or FRCs. We identified OsFD4 as a component of a FAC promoting flowering at the shoot apical meristem, downstream of OsFD1. The osfd4 mutants are late flowering and delay expression of genes promoting inflorescence development. Protein-protein interactions indicate an extensive network of contacts between several bZIPs and Gf14 proteins. Finally, we identified genomic regions bound by bZIPs with promotive and repressive effects on flowering. We conclude that distinct bZIPs orchestrate floral induction at the meristem and that FAC formation is largely combinatorial. While binding to the same consensus motif, their DNA-binding syntax is different, suggesting discriminatory functions.

The FUSED LEAVES1-ADHERENT1 regulatory module is required for maize cuticle development and organ separation.

All aerial epidermal cells in land plants are covered by the cuticle, an extracellular hydrophobic layer that provides protection against abiotic and biotic stresses and prevents organ fusion during development. Genetic and morphological analysis of the classic maize adherent1 (ad1) mutant was combined with genome-wide binding analysis of the maize MYB transcription factor FUSED LEAVES1 (FDL1), coupled with transcriptional profiling of fdl1 mutants. We show that AD1 encodes an epidermally-expressed 3-KETOACYL-CoA SYNTHASE (KCS) belonging to a functionally uncharacterized clade of KCS enzymes involved in cuticular wax biosynthesis. Wax analysis in ad1 mutants indicates that AD1 functions in the formation of very-long-chain wax components. We demonstrate that FDL1 directly binds to CCAACC core motifs present in AD1 regulatory regions to activate its expression. Over 2000 additional target genes of FDL1, including many involved in cuticle formation, drought response and cell wall organization, were also identified. Our results identify a regulatory module of cuticle biosynthesis in maize that is conserved across monocots and eudicots, and highlight previously undescribed factors in lipid metabolism, transport and signaling that coordinate organ development and cuticle formation.

VviNAC33 promotes organ de-greening and represses vegetative growth during the vegetative-to-mature phase transition in grapevine.

Plants undergo several developmental transitions during their life cycle. In grapevine, a perennial woody fruit crop, the transition from vegetative/green-to-mature/woody growth involves transcriptomic reprogramming orchestrated by a small group of genes encoding regulators, but the underlying molecular mechanisms are not fully understood. We investigated the function of the transcriptional regulator VviNAC33 by generating and characterizing transgenic overexpressing grapevine lines and a chimeric repressor, and by exploring its putative targets through a DNA affinity purification sequencing (DAP-seq) approach combined with transcriptomic data. We demonstrated that VviNAC33 induces leaf de-greening, inhibits organ growth and directly activates the expression of STAY-GREEN PROTEIN 1 (SGR1), which is involved in Chl and photosystem degradation, and AUTOPHAGY 8f (ATG8f), which is involved in the maturation of autophagosomes. Furthermore, we show that VviNAC33 directly inhibits AUXIN EFFLUX FACILITATOR PIN1, RopGEF1 and ATP SYNTHASE GAMMA CHAIN 1T (ATPC1), which are involved in photosystem II integrity and activity. Our results show that VviNAC33 plays a major role in terminating photosynthetic activity and organ growth as part of a regulatory network governing the vegetative-to-mature phase transition.

A Synthetic Approach Allows Rapid Characterization of the Maize Nuclear Auxin Response Circuit.

Auxin plays a key role across all land plants in growth and developmental processes. Although auxin signaling function has diverged and expanded, differences in the molecular functions of signaling components have largely been characterized in Arabidopsis (Arabidopsis thaliana). Here, we used the nuclear Auxin Response Circuit recapitulated in yeast (Saccharomyces cerevisiae) system to functionally annotate maize (Zea mays) auxin signaling components, focusing on genes expressed during the development of ear and tassel inflorescences. All 16 maize auxin/indole-3-acetic acid repressor proteins were degraded in response to auxin with rates that depended on both receptor and repressor identities. When fused to the maize TOPLESS homolog RAMOSA1 ENHANCER LOCUS2, maize auxin/indole-3-acetic acids were able to repress AUXIN RESPONSE FACTOR transcriptional activity. A complete auxin response circuit comprising all maize components, including the ZmAFB2/3 b1 maize AUXIN SIGNALING F-BOX (AFB) receptor, was fully functional. The ZmAFB2/3 b1 auxin receptor was more sensitive to hormone than AtAFB2 and allowed for rapid circuit activation upon auxin addition. These results validate the conserved role of predicted auxin response genes in maize as well as provide evidence that a synthetic approach can facilitate broader comparative studies across the wide range of species with sequenced genomes.

Improving architectural traits of maize inflorescences.

The domestication and improvement of maize resulted in radical changes in shoot architecture relative to its wild progenitor teosinte. In particular, critical modifications involved a reduction of branching and an increase in inflorescence size to meet the needs for human consumption and modern agricultural practices. Maize is a major contributor to global agricultural production by providing large and inexpensive quantities of food, animal feed, and ethanol. Maize is also a classic system for studying the genetic regulation of inflorescence formation and its enlarged female inflorescences directly influence seed production and yield. Studies on the molecular and genetic networks regulating meristem proliferation and maintenance, including receptor-ligand interactions, transcription factor regulation, and hormonal control, provide important insights into maize inflorescence development and reveal potential avenues for the targeted modification of specific architectural traits. In this review, we summarize recent findings on the molecular mechanisms controlling inflorescence formation and discuss how this knowledge can be applied to improve maize productivity in the face of present and future environmental challenges.

RAMOSA1 ENHANCER LOCUS2-Mediated Transcriptional Repression Regulates Vegetative and Reproductive Architecture.

Transcriptional repression in multicellular organisms orchestrates dynamic and precise gene expression changes that enable complex developmental patterns. Here, we present phenotypic and molecular characterization of the maize (Zea mays) transcriptional corepressor RAMOSA1 ENHANCER LOCUS2 (REL2), a unique member of the highly conserved TOPLESS (TPL) family. Analysis of single recessive mutations in rel2 revealed an array of vegetative and reproductive phenotypes, many related to defects in meristem initiation and maintenance. To better understand how REL2-mediated transcriptional complexes relate to rel2 phenotypes, we performed protein interaction assays and transcriptional profiling of mutant inflorescences, leading to the identification of different maize transcription factors and regulatory pathways that employ REL2 repression to control traits directly impacting maize yields. In addition, we used our REL2 interaction data to catalog conserved repression motifs present on REL2 interactors and showed that two of these, RLFGV- and DLN-type motifs, interact with the C-terminal WD40 domain of REL2 rather than the N terminus, which is known to bind LxLxL EAR motifs. These findings establish that the WD40 domain of TPL family proteins is an independent protein interaction surface that may work together with the N-terminal domain to allow the formation of large macromolecular complexes of functionally related transcription factors.

Expanding the Regulatory Network for Meristem Size in Plants.

The remarkable plasticity of post-embryonic plant development is due to groups of stem-cell-containing structures called meristems. In the shoot, meristems continuously produce organs such as leaves, flowers, and stems. Nearly two decades ago the WUSCHEL/CLAVATA (WUS/CLV) negative feedback loop was established as being essential for regulating the size of shoot meristems by maintaining a delicate balance between stem cell proliferation and cell recruitment for the differentiation of lateral primordia. Recent research in various model species (Arabidopsis, tomato, maize, and rice) has led to discoveries of additional components that further refine and improve the current model of meristem regulation, adding new complexity to a vital network for plant growth and productivity.

Auxin signaling modules regulate maize inflorescence architecture.

In plants, small groups of pluripotent stem cells called axillary meristems are required for the formation of the branches and flowers that eventually establish shoot architecture and drive reproductive success. To ensure the proper formation of new axillary meristems, the specification of boundary regions is required for coordinating their development. We have identified two maize genes, BARREN INFLORESCENCE1 and BARREN INFLORESCENCE4 (BIF1 and BIF4), that regulate the early steps required for inflorescence formation. BIF1 and BIF4 encode AUXIN/INDOLE-3-ACETIC ACID (Aux/IAA) proteins, which are key components of the auxin hormone signaling pathway that is essential for organogenesis. Here we show that BIF1 and BIF4 are integral to auxin signaling modules that dynamically regulate the expression of BARREN STALK1 (BA1), a basic helix-loop-helix (bHLH) transcriptional regulator necessary for axillary meristem formation that shows a striking boundary expression pattern. These findings suggest that auxin signaling directly controls boundary domains during axillary meristem formation and define a fundamental mechanism that regulates inflorescence architecture in one of the most widely grown crop species.

The boron efflux transporter ROTTEN EAR is required for maize inflorescence development and fertility.

Although boron has a relatively low natural abundance, it is an essential plant micronutrient. Boron deficiencies cause major crop losses in several areas of the world, affecting reproduction and yield in diverse plant species. Despite the importance of boron in crop productivity, surprisingly little is known about its effects on developing reproductive organs. We isolated a maize (Zea mays) mutant, called rotten ear (rte), that shows distinct defects in vegetative and reproductive development, eventually causing widespread sterility in its inflorescences, the tassel and the ear. Positional cloning revealed that rte encodes a membrane-localized boron efflux transporter, co-orthologous to the Arabidopsis thaliana BOR1 protein. Depending on the availability of boron in the soil, rte plants show a wide range of phenotypic defects that can be fully rescued by supplementing the soil with exogenous boric acid, indicating that rte is crucial for boron transport into aerial tissues. rte is expressed in cells surrounding the xylem in both vegetative and reproductive tissues and is required for meristem activity and organ development. We show that low boron supply to the inflorescences results in widespread defects in cell and cell wall integrity, highlighting the structural importance of boron in the formation of fully fertile reproductive organs.

Transport of boron by the tassel-less1 aquaporin is critical for vegetative and reproductive development in maize.

The element boron (B) is an essential plant micronutrient, and B deficiency results in significant crop losses worldwide. The maize (Zea mays) tassel-less1 (tls1) mutant has defects in vegetative and inflorescence development, comparable to the effects of B deficiency. Positional cloning revealed that tls1 encodes a protein in the aquaporin family co-orthologous to known B channel proteins in other species. Transport assays show that the TLS1 protein facilitates the movement of B and water into Xenopus laevis oocytes. B content is reduced in tls1 mutants, and application of B rescues the mutant phenotype, indicating that the TLS1 protein facilitates the movement of B in planta. B is required to cross-link the pectic polysaccharide rhamnogalacturonan II (RG-II) in the cell wall, and the percentage of RG-II dimers is reduced in tls1 inflorescences, indicating that the defects may result from altered cell wall properties. Plants heterozygous for both tls1 and rotten ear (rte), the proposed B efflux transporter, exhibit a dosage-dependent defect in inflorescence development under B-limited conditions, indicating that both TLS1 and RTE function in the same biological processes. Together, our data provide evidence that TLS1 is a B transport facilitator in maize, highlighting the importance of B homeostasis in meristem function.

BARREN STALK FASTIGIATE1 is an AT-hook protein required for the formation of maize ears.

Ears are the seed-bearing inflorescences of maize (Zea mays) plants and represent a crucial component of maize yield. The first step in the formation of ears is the initiation of axillary meristems in the axils of developing leaves. In the classic maize mutant barren stalk fastigiate1 (baf1), first discovered in the 1950s, ears either do not form or, if they do, are partially fused to the main stalk. We positionally cloned Baf1 and found that it encodes a transcriptional regulator containing an AT-hook DNA binding motif. Single coorthologs of Baf1 are found in syntenic regions of brachypodium (Brachypodium distachyon), rice (Oryza sativa), and sorghum (Sorghum bicolor), suggesting that the gene is likely present in all cereal species. Protein-protein interaction assays suggest that BAF1 is capable of forming homodimers and heterodimers with other members of the AT-hook family. Another transcriptional regulator required for ear initiation is the basic helix-loop-helix protein BARREN STALK1 (BA1). Genetic and expression analyses suggest that Baf1 is required to reach a threshold level of Ba1 expression for the initiation of maize ears. We propose that Baf1 functions in the demarcation of a boundary region essential for the specification of a stem cell niche.

Large-scale single-cell profiling of stem cells uncovers redundant regulators of shoot development and yield trait variation.

Stem cells in plant shoots are a rare population of cells that produce leaves, fruits and seeds, vital sources for food and bioethanol. Uncovering regulators expressed in these stem cells will inform crop engineering to boost productivity. Single-cell analysis is a powerful tool for identifying regulators expressed in specific groups of cells. However, accessing plant shoot stem cells is challenging. Recent single-cell analyses of plant shoots have not captured these cells, and failed to detect stem cell regulators like CLAVATA3 and WUSCHEL . In this study, we finely dissected stem cell-enriched shoot tissues from both maize and arabidopsis for single-cell RNA-seq profiling. We optimized protocols to efficiently recover thousands of CLAVATA3 and WUSCHEL expressed cells. A cross-species comparison identified conserved stem cell regulators between maize and arabidopsis. We also performed single-cell RNA-seq on maize stem cell overproliferation mutants to find additional candidate regulators. Expression of candidate stem cell genes was validated using spatial transcriptomics, and we functionally confirmed roles in shoot development. These candidates include a family of ribosome-associated RNA-binding proteins, and two families of sugar kinase genes related to hypoxia signaling and cytokinin hormone homeostasis. These large-scale single-cell profiling of stem cells provide a resource for mining stem cell regulators, which show significant association with yield traits. Overall, our discoveries advance the understanding of shoot development and open avenues for manipulating diverse crops to enhance food and energy security.

Transcription factor binding site divergence across maize inbred lines drives transcriptional and phenotypic variation.

Regulatory elements are important constituents of plant genomes that have shaped ancient and modern crops. Their identification, function, and diversity in crop genomes however are poorly characterized, thus limiting our ability to harness their power for further agricultural advances using induced or natural variation. Here, we use DNA affinity purification-sequencing (DAP-seq) to map transcription factor (TF) binding events for 200 maize TFs belonging to 30 distinct families and heterodimer pairs in two distinct inbred lines historically used for maize hybrid plant production, providing empirical binding site annotation for 5.3% of the maize genome. TF binding site comparison in B73 and Mo17 inbreds reveals widespread differences, driven largely by structural variation, that correlate with gene expression changes. TF binding site presence-absence variation helps clarify complex QTL such as vgt1, an important determinant of maize flowering time, and DICE, a distal enhancer involved in herbivore resistance. Modification of TF binding regions via CRISPR-Cas9 mediated editing alters target gene expression and phenotype. Our functional catalog of maize TF binding events enables collective and comparative TF binding analysis, and highlights its value for agricultural improvement.