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.

Plant Noncoding RNAs: Hidden Players in Development and Stress Responses.

A large and significant portion of eukaryotic transcriptomes consists of noncoding RNAs (ncRNAs) that have minimal or no protein-coding capacity but are functional. Diverse ncRNAs, including both small RNAs and long ncRNAs (lncRNAs), play essential regulatory roles in almost all biological processes by modulating gene expression at the transcriptional and posttranscriptional levels. In this review, we summarize the current knowledge of plant small RNAs and lncRNAs, with a focus on their biogenesis, modes of action, local and systemic movement, and functions at the nexus of plant development and environmental responses. The complex connections among small RNAs, lncRNAs, and small peptides in plants are also discussed, along with the challenges of identifying and investigating new classes of ncRNAs.

Natural temperature fluctuations promote COOLAIR regulation of FLC.

Plants monitor many aspects of their fluctuating environments to help align their development with seasons. Molecular understanding of how noisy temperature cues are registered has emerged from dissection of vernalization in Arabidopsis, which involves a multiphase cold-dependent silencing of the floral repressor locus FLOWERING LOCUS C (FLC). Cold-induced transcriptional silencing precedes a low probability PRC2 epigenetic switching mechanism. The epigenetic switch requires the absence of warm temperatures as well as long-term cold exposure. However, the natural temperature inputs into the earlier transcriptional silencing phase are less well understood. Here, through investigation of Arabidopsis accessions in natural and climatically distinct field sites, we show that the first seasonal frost strongly induces expression of COOLAIR, the antisense transcripts at FLC Chamber experiments delivering a constant mean temperature with different fluctuations showed the freezing induction of COOLAIR correlates with stronger repression of FLC mRNA. Identification of a mutant that ectopically activates COOLAIR revealed how COOLAIR up-regulation can directly reduce FLC expression. Consistent with this, transgenes designed to knockout COOLAIR perturbed the early phase of FLC silencing. However, all transgenes designed to remove COOLAIR resulted in increased production of novel convergent FLC antisense transcripts. Our study reveals how natural temperature fluctuations promote COOLAIR regulation of FLC, with the first autumn frost acting as a key indicator of autumn/winter arrival.

Reciprocal Regulation of the TOR Kinase and ABA Receptor Balances Plant Growth and Stress Response.

As sessile organisms, plants must adapt to variations in the environment. Environmental stress triggers various responses, including growth inhibition, mediated by the plant hormone abscisic acid (ABA). The mechanisms that integrate stress responses with growth are poorly understood. Here, we discovered that the Target of Rapamycin (TOR) kinase phosphorylates PYL ABA receptors at a conserved serine residue to prevent activation of the stress response in unstressed plants. This phosphorylation disrupts PYL association with ABA and with PP2C phosphatase effectors, leading to inactivation of SnRK2 kinases. Under stress, ABA-activated SnRK2s phosphorylate Raptor, a component of the TOR complex, triggering TOR complex dissociation and inhibition. Thus, TOR signaling represses ABA signaling and stress responses in unstressed conditions, whereas ABA signaling represses TOR signaling and growth during times of stress. Plants utilize this conserved phospho-regulatory feedback mechanism to optimize the balance of growth and stress responses.

Flavonol-mediated stabilization of PIN efflux complexes regulates polar auxin transport.

The transport of auxin controls the rate, direction and localization of plant growth and development. The course of auxin transport is defined by the polar subcellular localization of the PIN proteins, a family of auxin efflux transporters. However, little is known about the composition and regulation of the PIN protein complex. Here, using blue-native PAGE and quantitative mass spectrometry, we identify native PIN core transport units as homo- and heteromers assembled from PIN1, PIN2, PIN3, PIN4 and PIN7 subunits only. Furthermore, we show that endogenous flavonols stabilize PIN dimers to regulate auxin efflux in the same way as does the auxin transport inhibitor 1-naphthylphthalamic acid (NPA). This inhibitory mechanism is counteracted both by the natural auxin indole-3-acetic acid and by phosphomimetic amino acids introduced into the PIN1 cytoplasmic domain. Our results lend mechanistic insights into an endogenous control mechanism which regulates PIN function and opens the way for a deeper understanding of the protein environment and regulation of the polar auxin transport complex.

The calcium transporter ANNEXIN1 mediates cold-induced calcium signaling and freezing tolerance in plants.

The transient elevation of cytosolic free calcium concentration ([Ca2+ ]cyt ) induced by cold stress is a well-established phenomenon; however, the underlying mechanism remains elusive. Here, we report that the Ca2+ -permeable transporter ANNEXIN1 (AtANN1) mediates cold-triggered Ca2+ influx and freezing tolerance in Arabidopsis thaliana. The loss of function of AtANN1 substantially impaired freezing tolerance, reducing the cold-induced [Ca2+ ]cyt increase and upregulation of the cold-responsive CBF and COR genes. Further analysis showed that the OST1/SnRK2.6 kinase interacted with and phosphorylated AtANN1, which consequently enhanced its Ca2+ transport activity, thereby potentiating Ca2+ signaling. Consistent with these results and freezing sensitivity of ost1 mutants, the cold-induced [Ca2+ ]cyt elevation in the ost1-3 mutant was reduced. Genetic analysis indicated that AtANN1 acts downstream of OST1 in responses to cold stress. Our data thus uncover a cascade linking OST1-AtANN1 to cold-induced Ca2+ signal generation, which activates the cold response and consequently enhances freezing tolerance in Arabidopsis.

Prediction of conserved and variable heat and cold stress response in maize using cis-regulatory information.

Changes in gene expression are important for responses to abiotic stress. Transcriptome profiling of heat- or cold-stressed maize genotypes identifies many changes in transcript abundance. We used comparisons of expression responses in multiple genotypes to identify alleles with variable responses to heat or cold stress and to distinguish examples of cis- or trans-regulatory variation for stress-responsive expression changes. We used motifs enriched near the transcription start sites (TSSs) for thermal stress-responsive genes to develop predictive models of gene expression responses. Prediction accuracies can be improved by focusing only on motifs within unmethylated regions near the TSS and vary for genes with different dynamic responses to stress. Models trained on expression responses in a single genotype and promoter sequences provided lower performance when applied to other genotypes but this could be improved by using models trained on data from all three genotypes tested. The analysis of genes with cis-regulatory variation provides evidence for structural variants that result in presence/absence of transcription factor binding sites in creating variable responses. This study provides insights into cis-regulatory motifs for heat- and cold-responsive gene expression and defines a framework for developing models to predict expression responses across multiple genotypes.

Misregulation of MYB16 expression causes stomatal cluster formation by disrupting polarity during asymmetric cell divisions.

Stomatal pores and the leaf cuticle regulate evaporation from the plant body and balance the tradeoff between photosynthesis and water loss. MYB16, encoding a transcription factor involved in cutin biosynthesis, is expressed in stomatal lineage ground cells, suggesting a link between cutin biosynthesis and stomatal development. Here, we show that the downregulation of MYB16 in meristemoids is directly mediated by the stomatal master transcription factor SPEECHLESS (SPCH) in Arabidopsis thaliana. The suppression of MYB16 before an asymmetric division is crucial for stomatal patterning, as its overexpression or ectopic expression in meristemoids increased stomatal density and resulted in the formation of stomatal clusters, as well as affecting the outer cell wall structure. Expressing a cutinase gene in plants ectopically expressing MYB16 reduced stomatal clustering, suggesting that cutin affects stomatal signaling or the polarity setup in asymmetrically dividing cells. The clustered stomatal phenotype was rescued by overexpressing EPIDERMAL PATTERNING FACTOR2, suggesting that stomatal signaling was still functional in these plants. Growing seedlings ectopically expressing MYB16 on high-percentage agar plates to modulate tensile strength rescued the polarity and stomatal cluster defects of these seedlings. Therefore, the inhibition of MYB16 expression by SPCH in the early stomatal lineage is required to correctly place the polarity protein needed for stomatal patterning during leaf morphogenesis.

Gene regulatory networks for compatible versus incompatible grafts identify a role for SlWOX4 during junction formation.

Grafting has been adopted for a wide range of crops to enhance productivity and resilience; for example, grafting of Solanaceous crops couples disease-resistant rootstocks with scions that produce high-quality fruit. However, incompatibility severely limits the application of grafting and graft incompatibility remains poorly understood. In grafts, immediate incompatibility results in rapid death, but delayed incompatibility can take months or even years to manifest, creating a significant economic burden for perennial crop production. To gain insight into the genetic mechanisms underlying this phenomenon, we developed a model system using heterografting of tomato (Solanum lycopersicum) and pepper (Capsicum annuum). These grafted plants express signs of anatomical junction failure within the first week of grafting. By generating a detailed timeline for junction formation, we were able to pinpoint the cellular basis for this delayed incompatibility. Furthermore, we inferred gene regulatory networks for compatible self-grafts and incompatible heterografts based on these key anatomical events, which predict core regulators for grafting. Finally, we examined the role of vascular development in graft formation and uncovered SlWOX4 as a potential regulator of graft compatibility. Following this predicted regulator up with functional analysis, we show that Slwox4 homografts fail to form xylem bridges across the junction, demonstrating that indeed, SlWOX4 is essential for vascular reconnection during grafting, and may function as an early indicator of graft failure.

Metabolomic selection for enhanced fruit flavor.

Although they are staple foods in cuisines globally, many commercial fruit varieties have become progressively less flavorful over time. Due to the cost and difficulty associated with flavor phenotyping, breeding programs have long been challenged in selecting for this complex trait. To address this issue, we leveraged targeted metabolomics of diverse tomato and blueberry accessions and their corresponding consumer panel ratings to create statistical and machine learning models that can predict sensory perceptions of fruit flavor. Using these models, a breeding program can assess flavor ratings for a large number of genotypes, previously limited by the low throughput of consumer sensory panels. The ability to predict consumer ratings of liking, sweet, sour, umami, and flavor intensity was evaluated by a 10-fold cross-validation, and the accuracies of 18 different models were assessed. The prediction accuracies were high for most attributes and ranged from 0.87 for sourness intensity in blueberry using XGBoost to 0.46 for overall liking in tomato using linear regression. Further, the best-performing models were used to infer the flavor compounds (sugars, acids, and volatiles) that contribute most to each flavor attribute. We found that the variance decomposition of overall liking score estimates that 42% and 56% of the variance was explained by volatile organic compounds in tomato and blueberry, respectively. We expect that these models will enable an earlier incorporation of flavor as breeding targets and encourage selection and release of more flavorful fruit varieties.

HYL1-CLEAVAGE SUBTILASE 1 (HCS1) suppresses miRNA biogenesis in response to light-to-dark transition.

The core plant microprocessor consists of DICER-LIKE 1 (DCL1), SERRATE (SE), and HYPONASTIC LEAVES 1 (HYL1) and plays a pivotal role in microRNA (miRNA) biogenesis. However, the proteolytic regulation of each component remains elusive. Here, we show that HYL1-CLEAVAGE SUBTILASE 1 (HCS1) is a cytoplasmic protease for HYL1-destabilization. HCS1-excessiveness reduces HYL1 that disrupts miRNA biogenesis, while HCS1-deficiency accumulates HYL1. Consistently, we identified the HYL1K154A mutant that is insensitive to the proteolytic activity of HCS1, confirming the importance of HCS1 in HYL1 proteostasis. Moreover, HCS1-activity is regulated by light/dark transition. Under light, cytoplasmic CONSTITUTIVE PHOTOMORPHOGENIC 1 (COP1) E3 ligase suppresses HCS1-activity. COP1 sterically inhibits HCS1 by obstructing HYL1 access into the catalytic sites of HCS1. In contrast, darkness unshackles HCS1-activity for HYL1-destabilization due to nuclear COP1 relocation. Overall, the COP1-HYL1-HCS1 network may integrate two essential cellular pathways: the miRNA-biogenetic pathway and light signaling pathway.

Arabidopsis TBP-ASSOCIATED FACTOR 12 ortholog NOBIRO6 controls root elongation with unfolded protein response cofactor activity.

Plant root growth is indeterminate but continuously responds to environmental changes. We previously reported on the severe root growth defect of a double mutant in bZIP17 and bZIP28 (bz1728) modulating the unfolded protein response (UPR). To elucidate the mechanism by which bz1728 seedlings develop a short root, we obtained a series of bz1728 suppressor mutants, called nobiro, for rescued root growth. We focused here on nobiro6, which is defective in the general transcription factor component TBP-ASSOCIATED FACTOR 12b (TAF12b). The expression of hundreds of genes, including the bZIP60-UPR regulon, was induced in the bz1728 mutant, but these inductions were markedly attenuated in the bz1728nobiro6 mutant. In view of this, we assigned transcriptional cofactor activity via physical interaction with bZIP60 to NOBIRO6/TAF12b. The single nobiro6/taf12b mutant also showed an altered sensitivity to endoplasmic reticulum stress for both UPR and root growth responses, demonstrating that NOBIRO6/TAF12b contributes to environment-responsive root growth control through UPR.

Stable establishment of organ polarity occurs several plastochrons before primordium outgrowth in Arabidopsis.

In many species, leaves are initiated at the flanks of shoot meristems. Subsequent growth usually occurs mainly in the plane of the leaf blade, which leads to the formation of a bifacial leaf with dorsoventral identities. In a classical set of surgical experiments in potato meristems, Sussex provided evidence that dorsoventrality depends on a signal emanating from the meristem center. Although these results could be reproduced in tomato, this concept has been debated. We revisited these experiments in Arabidopsis, in which a range of markers are available to target the precise site of ablation. Using specific markers for organ founder cells and dorsoventral identity, we were unable to perturb the polarity of leaves and sepals long before organ outgrowth. Although results in Solanaceae suggested that dorsoventral patterning was unstable during early development, we found that, in Arabidopsis, the local information contained within and around the primordium is able to withstand major invasive perturbations, long before polarity is fully established.

PIF4 negatively modulates cold tolerance in tomato anthers via temperature-dependent regulation of tapetal cell death.

Extreme temperature conditions seriously impair male reproductive development in plants; however, the molecular mechanisms underlying the response of anthers to extreme temperatures remain poorly described. The transcription factor phytochrome-interacting factor4 (PIF4) acts as a hub that integrates multiple signaling pathways to regulate thermosensory growth and architectural adaptation in plants. Here, we report that SlPIF4 in tomato (Solanum lycopersicum) plays a pivotal role in regulating cold tolerance in anthers. CRISPR (clustered regularly interspaced short palindromic repeats)-associated nuclease Cas9-generated SlPIF4 knockout mutants showed enhanced cold tolerance in pollen due to reduced temperature sensitivity of the tapetum, while overexpressing SlPIF4 conferred pollen abortion by delaying tapetal programmed cell death (PCD). SlPIF4 directly interacts with SlDYT1, a direct upstream regulator of SlTDF1, both of which (SlDYT1 and SlTDF1) play important roles in regulating tapetum development and tapetal PCD. Moderately low temperature (MLT) promotes the transcriptional activation of SlTDF1 by the SlPIF4-SlDYT1 complex, resulting in pollen abortion, while knocking out SlPIF4 blocked the MLT-induced activation of SlTDF1. Furthermore, SlPIF4 directly binds to the canonical E-box sequence in the SlDYT1 promoter. Collectively, these findings suggest that SlPIF4 negatively regulates cold tolerance in anthers by directly interacting with the tapetal regulatory module in a temperature-dependent manner. Our results shed light on the molecular mechanisms underlying the adaptation of anthers to low temperatures.

Arabidopsis NF-YCs play dual roles in repressing brassinosteroid biosynthesis and signaling during light-regulated hypocotyl elongation.

Light functions as the primary environmental stimulus and brassinosteroids (BRs) as important endogenous growth regulators throughout the plant lifecycle. Photomorphogenesis involves a series of vital developmental processes that require the suppression of BR-mediated seedling growth, but the mechanism underlying the light-controlled regulation of the BR pathway remains unclear. Here, we reveal that nuclear factor YC proteins (NF-YCs) function as essential repressors of the BR pathway during light-controlled hypocotyl growth in Arabidopsis thaliana. In the light, NF-YCs inhibit BR biosynthesis by directly targeting the promoter of the BR biosynthesis gene BR6ox2 and repressing its transcription. NF-YCs also interact with BIN2, a critical repressor of BR signaling, and facilitate its stabilization by promoting its Tyr200 autophosphorylation, thus inhibiting the BR signaling pathway. Consistently, loss-of-function mutants of NF-YCs show etiolated growth and constitutive BR responses, even in the light. Our findings uncover a dual role of NF-YCs in repressing BR biosynthesis and signaling, providing mechanistic insights into how light antagonizes the BR pathway to ensure photomorphogenic growth in Arabidopsis.

Nonspecific phospholipase C4 hydrolyzes phosphosphingolipids and sustains plant root growth during phosphate deficiency.

Phosphate is a vital macronutrient for plant growth, and its availability in soil is critical for agricultural sustainability and productivity. A substantial amount of cellular phosphate is used to synthesize phospholipids for cell membranes. Here, we identify a key enzyme, nonspecific phospholipase C4 (NPC4) that is involved in phosphosphingolipid hydrolysis and remodeling in Arabidopsis during phosphate starvation. The level of glycosylinositolphosphorylceramide (GIPC), the most abundant sphingolipid in Arabidopsis thaliana, decreased upon phosphate starvation. NPC4 was highly induced by phosphate deficiency, and NPC4 knockouts in Arabidopsis decreased the loss of GIPC and impeded root growth during phosphate starvation. Enzymatic analysis showed that NPC4 hydrolyzed GIPC and displayed a higher activity toward GIPC as a substrate than toward the common glycerophospholipid phosphatidylcholine. NPC4 was associated with the plasma membrane lipid rafts in which GIPC is highly enriched. These results indicate that NPC4 uses GIPC as a substrate in planta and the NPC4-mediated sphingolipid remodeling plays a positive role in root growth in Arabidopsis response to phosphate deficiency.

The B3 domain-containing transcription factor ZmABI19 coordinates expression of key factors required for maize seed development and grain filling.

Grain filling in maize (Zea mays) is regulated by a group of spatiotemporally synchronized transcription factors (TFs), but the factors that coordinate their expression remain unknown. We used the promoter of the grain filling-specific TF gene Opaque2 (O2) to screen upstream regulatory factors and identified a B3 domain TF, ZmABI19, that directly binds to the O2 promoter for transactivation. zmabi19 mutants displayed developmental defects in the endosperm and embryo, and mature kernels were opaque and reduced in size. The accumulation of zeins, starch and lipids dramatically decreased in zmabi19 mutants. RNA sequencing revealed an alteration of the nutrient reservoir activity and starch and sucrose metabolism in zmabi19 endosperms, and plant phytohormone signal transduction and lipid metabolism in zmabi19 embryos. Chromatin immunoprecipitation followed by sequencing coupled with differential expression analysis identified 106 high-confidence direct ZmABI19 targets. ZmABI19 directly regulates multiple key grain filling TFs including O2, Prolamine-box binding factor 1, ZmbZIP22, NAC130, and Opaque11 in the endosperm and Viviparous1 in the embryo. A number of phytohormone-related genes were also bound and regulated by ZmABI19. Our results demonstrate that ZmABI19 functions as a grain filling initiation regulator. ZmABI19 roles in coupling early endosperm and embryo development are also discussed.

Seed-coat protective neolignans are produced by the dirigent protein AtDP1 and the laccase AtLAC5 in Arabidopsis.

Lignans/neolignans are generally synthesized from coniferyl alcohol (CA) in the cinnamate/monolignol pathway by oxidation to generate the corresponding radicals with subsequent stereoselective dimerization aided by dirigent proteins (DIRs). Genes encoding oxidases and DIRs for neolignan biosynthesis have not been identified previously. In Arabidopsis thaliana, the DIR AtDP1/AtDIR12 plays an essential role in the 8-O-4' coupling in neolignan biosynthesis by unequivocal structural determination of the compound missing in the atdp1 mutant as a sinapoylcholine (SC)-conjugated neolignan, erythro-3-{4-[2-hydroxy-2-(4-hydroxy-3-methoxyphenyl)-1-hydroxymethylethoxy]-3,5-dimethoxyphenyl}acryloylcholine. Phylogenetic analyses showed that AtDP1/AtDIR12 belongs to the DIR-a subfamily composed of DIRs for 8-8' coupling of monolignol radicals. AtDP1/AtDIR12 is specifically expressed in outer integument 1 cells in developing seeds. As a putative oxidase for neolignan biosynthesis, we focused on AtLAC5, a laccase gene coexpressed with AtDP1/AtDIR12. In lac5 mutants, the abundance of feruloylcholine (FC)-conjugated neolignans decreased to a level comparable to those in the atdp1 mutant. In addition, SC/FC-conjugated neolignans were missing in the seeds of mutants defective in SCT/SCPL19, an enzyme that synthesizes SC. These results strongly suggest that AtDP1/AtDIR12 and AtLAC5 are involved in neolignan biosynthesis via SC/FC. A tetrazolium penetration assay showed that seed coat permeability increased in atdp1 mutants, suggesting a protective role of neolignans in A. thaliana seeds.

Ectopic expression of Triticum polonicum VRT-A2 underlies elongated glumes and grains in hexaploid wheat in a dosage-dependent manner.

Flower development is an important determinant of grain yield in crops. In wheat (Triticum spp.), natural variation for the size of spikelet and floral organs is particularly evident in Triticum turgidum ssp. polonicum (also termed Triticum polonicum), a tetraploid subspecies of wheat with long glumes, lemmas, and grains. Using map-based cloning, we identified VEGETATIVE TO REPRODUCTIVE TRANSITION 2 (VRT2), which encodes a MADS-box transcription factor belonging to the SHORT VEGETATIVE PHASE family, as the gene underlying the T. polonicum long-glume (P1) locus. The causal P1 mutation is a sequence rearrangement in intron-1 that results in ectopic expression of the T. polonicum VRT-A2 allele. Based on allelic variation studies, we propose that the intron-1 mutation in VRT-A2 is the unique T. polonicum subspecies-defining polymorphism, which was later introduced into hexaploid wheat via natural hybridizations. Near-isogenic lines differing for the P1 locus revealed a gradient effect of P1 across spikelets and within florets. Transgenic lines of hexaploid wheat carrying the T. polonicum VRT-A2 allele show that expression levels of VRT-A2 are highly correlated with spike, glume, grain, and floral organ length. These results highlight how changes in expression profiles, through variation in cis-regulation, can affect agronomic traits in a dosage-dependent manner in polyploid crops.

Time of the day prioritizes the pool of translating mRNAs in response to heat stress.

The circadian clock helps organisms to anticipate and coordinate gene regulatory responses to changes in environmental stimuli. Under growth limiting temperatures, the time of the day modulates the accumulation of polyadenylated mRNAs. In response to heat stress, plants will conserve energy and selectively translate mRNAs. How the clock and/or the time of the day regulates polyadenylated mRNAs bound by ribosomes in response to heat stress is unknown. In-depth analysis of Arabidopsis thaliana translating mRNAs found that the time of the day gates the response of approximately one-third of the circadian-regulated heat-responsive translatome. Specifically, the time of the day and heat stress interact to prioritize the pool of mRNAs in cue to be translated. For a subset of mRNAs, we observed a stronger gated response during the day, and preferentially before the peak of expression. We propose previously overlooked transcription factors (TFs) as regulatory nodes and show that the clock plays a role in the temperature response for select TFs. When the stress was removed, the redefined priorities for translation recovered within 1 h, though slower recovery was observed for abiotic stress regulators. Through hierarchical network connections between clock genes and prioritized TFs, our work provides a framework to target key nodes underlying heat stress tolerance throughout the day.