Cistrome and Epicistrome Features Shape the Regulatory DNA Landscape

STOP1 regulates CCX1-mediated Ca2+ homeostasis for plant adaptation to Ca2+ deprivation

Calcium (Ca) is essential for plant growth and stress adaptation, yet its availability is often limited in acidic soils, posing a major threat to crop production. Understanding the intricate mechanisms orchestrating plant adaptation to Ca deficiency remains elusive. Here, we show that the Ca deficiency-enhanced nuclear accumulation of the transcription factor SENSITIVE TO PROTON RHIZOTOXICITY 1 (STOP1) in Arabidopsis thaliana confers tolerance to Ca deprivation, with the global transcriptional responses triggered by Ca deprivation largely impaired in the stop1 mutant. Notably, STOP1 activates the Ca deprivation-induced expression of CATION/Ca2+ EXCHANGER 1 (CCX1) by directly binding to its promoter region, which facilitates Ca2+ efflux from endoplasmic reticulum to cytosol to maintain Ca homeostasis. Consequently, the constitutive expression of CCX1 in the stop1 mutant partially rescues the Ca deficiency phenotype by increasing Ca content in the shoots. These findings uncover the pivotal role of the STOP1-CCX1 axis in plant adaptation to low Ca, offering alternative manipulating strategies to improve plant Ca nutrition in acidic soils and extending our understanding of the multifaceted role of STOP1.

VIP1 and its close homologs confer mechanical stress tolerance in Arabidopsis leaves

VIP1, an Arabidopsis thaliana basic leucine zipper transcription factor, and its close homologs are imported from the cytoplasm to the nucleus when cells are exposed to mechanical stress. They bind to AGCTG (G/T) and regulate mechanical stress responses in roots. However, their role in leaves is unclear. To clarify this, mutant lines (QM1 and QM2) that lack the functions of VIP1 and its close homologs (bZIP29, bZIP30 and PosF21) were generated. Brushing more severely damaged QM1 and QM2 leaves than wild-type leaves. Genes regulating stress responses and cell wall properties were downregulated in brushed QM2 leaves and upregulated in brushed VIP1-GFP-overexpressing (VIP1-GFPox) leaves compared to wild-type leaves in a transcriptome analysis. The VIP1-binding sequence AGCTG (G/T) was enriched in the promoters of genes downregulated in brushed QM2 leaves compared to wild-type leaves and in those upregulated in brushed VIP1-GFPox leaves. Calmodulin-binding transcription activators (CAMTAs) are known regulators of mechanical stress responses, and the CAMTA-binding sequence CGCGT was enriched in the promoters of genes upregulated in the brushed QM2 leaves and in those downregulated in the brushed VIP1-GFPox leaves. These findings suggest that VIP1 and its homologs upregulate genes via AGCTG (G/T) and influence CAMTA-dependent gene expression to enhance mechanical stress tolerance in leaves.

Overexpression of a pearl millet WRKY transcription factor gene, PgWRKY74, in Arabidopsis retards shoot growth under dehydration and salinity-stressed conditions

Pearl millet (Cenchrus americanus) is a cereal crop that can tolerate high temperatures, drought, and low-fertility conditions where other crops lose productivity. However, genes regulating this ability are largely unknown. Transcription factors (TFs) regulate transcription of their target genes, regulate downstream biological processes, and thus are candidates for regulators of such tolerance of pearl millet. PgWRKY74 encodes a group IIc WRKY TF in pearl millet and is downregulated by drought. PgWRKY74 may have a role in drought tolerance. The objective of this study was to gain insights into the physiological and biochemical functions of PgWRKY74. Yeast one-hybrid and gel shift assays were performed to examine transcriptional activation potential and deoxyribonucleic acid (DNA)-binding ability, respectively. Transgenic Arabidopsis thaliana plants overexpressing PgWRKY74-green fluorescent protein (GFP) fusion gene were generated and tested for growth and stress-responsive gene expression under mannitol and NaCl-stressed conditions. A construct with PgWRKY74 enabled yeast reporter cells to survive on test media in the yeast one-hybrid assays. The electrophoretic mobility of DNA with putative WRKY TF-binding motifs was lower in the presence of a recombinant PgWRKY74 protein than its absence. The PgWRKY74-GFP-overexpressing Arabidopsis plants exhibited smaller rosette areas than did wild-type plants under mannitol-stressed and NaCl-stressed conditions, and exhibited weaker expression of RD29B, which is induced by the stress-related phytohormone abscisic acid (ABA), under the mannitol-stressed condition. PgWRKY74 have transcriptional activation potential and DNA-binding ability, and can negatively regulate plant responses to mannitol and NaCl stresses, possibly by decreasing ABA levels or ABA sensitivity.

Dynamics of DNA methylation and its impact on plant embryogenesis

Flowering plants exhibit unique DNA methylation dynamics during development. Particular attention can be focused on seed development and the embryo, which represents the starting point of the sporophytic life cycle. A build-up of CHH methylation is now recognized as highly characteristic of embryo development. This process is thought to occur in order to silence potentially harmful transposable element expression, though roles in promoting seed dormancy and dessication tolerance have also been revealed. Recent studies show that increased CHH methylation in embryos inhabits both novel loci, unmethylated elsewhere in the plant, as well as shared loci, exhibiting more dense methylation. The role of DNA methylation in cis-regulatory gene regulation in plants is less well established compared to mammals, and here we discuss both transposable element regulation and the potential role of DNA methylation in dynamic gene expression.

Widespread position-dependent transcriptional regulatory sequences in plants

Much of what we know about eukaryotic transcription stems from animals and yeast; however, plants evolved separately for over a billion years, leaving ample time for divergence in transcriptional regulation. Here we set out to elucidate fundamental properties of cis-regulatory sequences in plants. Using massively parallel reporter assays across four plant species, we demonstrate the central role of sequences downstream of the transcription start site (TSS) in transcriptional regulation. Unlike animal enhancers that are position independent, plant regulatory elements depend on their position, as altering their location relative to the TSS significantly affects transcription. We highlight the importance of the region downstream of the TSS in regulating transcription by identifying a DNA motif that is conserved across vascular plants and is sufficient to enhance gene expression in a dose-dependent manner. The identification of a large number of position-dependent enhancers points to fundamental differences in gene regulation between plants and animals.

Genomes of diverse Actinidia species provide insights into cis-regulatory motifs and genes associated with critical traits

BACKGROUND

Kiwifruit, belonging to the genus Actinidia, represents a unique fruit crop characterized by its modern cultivars being genetically diverse and exhibiting remarkable variations in morphological traits and adaptability to harsh environments. However, the genetic mechanisms underlying such morphological diversity remain largely elusive.

RESULTS

We report the high-quality genomes of five Actinidia species, including Actinidia longicarpa, A. macrosperma, A. polygama, A. reticulata, and A. rufa. Through comparative genomics analyses, we identified three whole genome duplication events shared by the Actinidia genus and uncovered rapidly evolving gene families implicated in the development of characteristic kiwifruit traits, including vitamin C (VC) content and fruit hairiness. A range of structural variations were identified, potentially contributing to the phenotypic diversity in kiwifruit. Notably, phylogenomic analyses revealed 76 cis-regulatory elements within the Actinidia genus, predominantly associated with stress responses, metabolic processes, and development. Among these, five motifs did not exhibit similarity to known plant motifs, suggesting the presence of possible novel cis-regulatory elements in kiwifruit. Construction of a pan-genome encompassing the nine Actinidia species facilitated the identification of gene DTZ79_23g14810 specific to species exhibiting extraordinarily high VC content. Expression of DTZ79_23g14810 is significantly correlated with the dynamics of VC concentration, and its overexpression in the transgenic roots of kiwifruit plants resulted in increased VC content.

CONCLUSIONS

Collectively, the genomes and pan-genome of diverse Actinidia species not only enhance our understanding of fruit development but also provide a valuable genomic resource for facilitating the genome-based breeding of kiwifruit.

A Bamboo HD-Zip Transcription Factor PeHDZ72 Conferred Drought Tolerance by Promoting Sugar and Water Transport

Drought drastically affects plant growth, development and productivity. Plants respond to drought stress by enhancing sugar accumulation and water transport. Homeodomain-leucine zipper (HD-Zip) transcription factors (TFs) participate in various aspects of plant growth and stress response. However, the internal regulatory mechanism of HD-Zips in moso bamboo (Phyllostachys edulis) remains largely unknown. In this study, we identified an HD-Zip member, PeHDZ72, which was highly expressed in bamboo shoots and roots and was induced by drought. Furthermore, PeSTP_46019, PeSWEET_23178 and PeTIP4-3 were identified as downstream genes of PeHDZ72 in moso bamboo by DAP-seq. The expressions of these three genes were all induced by drought stress. Y1H, DLR and GUS activity assays demonstrated that PeHDZ72 could bind to three types of HD-motifs in the promoters of these three genes. Overexpression of PeHDZ72 led to a remarkable enhancement in drought tolerance in transgenic rice, with significantly improved soluble sugar and sucrose contents. Meanwhile, the expressions of OsSTPs, OsSWEETs and OsTIP were all upregulated in transgenic rice under drought stress. Overall, our results indicate that drought stress might induce the expression of PeHDZ72, which in turn activated downstream genes PeSTP_46019, PeSWEET_23178 and PeTIP4-3, contributing to the improvement of cellular osmotic potential in moso bamboo in response to drought stress.

OsMBF1a Facilitates Seed Germination by Regulating Biosynthesis of Gibberellic Acid and Abscisic Acid in Rice

Seed germination is a pivotal stage in the plant life cycle, orchestrated by a myriad of internal and external factors, notably plant hormones. The underlying molecular mechanisms governing rice seed germination remain largely unelucidated. Herein, we uncover OsMBF1a as a crucial regulatory factor that employs a dual strategy to promote seed germination: positively activating genes involved in gibberellin (GA) biosynthesis pathways, while negatively regulating key genes responsible for abscisic acid (ABA) synthesis. Furthermore, OsMBF1a modulates the endogenous levels of ABA and GA in rice seeds, reinforcing its central role in the germination process. The expression of ZmMBF1a and ZmMBF1b, the homologous genes in maize, in rice seeds similarly affects germination, indicating the conserved functionality of MBF1 family genes in regulating seed germination. This study provides novel insights into the molecular mechanisms underlying rice seed germination and underscores the significance of MBF1 family genes in plant growth and development.

MYB-related transcription factors control chloroplast biogenesis

Chloroplast biogenesis is dependent on master regulators from the GOLDEN2-LIKE (GLK) family of transcription factors. However, glk mutants contain residual chlorophyll, indicating that other proteins must be involved. Here, we identify MYB-related transcription factors as regulators of chloroplast biogenesis in the liverwort Marchantia polymorpha and angiosperm Arabidopsis thaliana. In both species, double-mutant alleles in MYB-related genes show very limited chloroplast development, and photosynthesis gene expression is perturbed to a greater extent than in GLK mutants. Genes encoding enzymes of chlorophyll biosynthesis are controlled by MYB-related and GLK proteins, whereas those allowing CO2 fixation, photorespiration, and photosystem assembly and repair require MYB-related proteins. Regulation between the MYB-related and GLK transcription factors appears more extensive in A. thaliana than in M. polymorpha. Thus, MYB-related and GLK genes have overlapping as well as distinct targets. We conclude that MYB-related and GLK transcription factors orchestrate chloroplast development in land plants.

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.

Architecture and dynamics of the abscisic acid gene regulatory network

The plant hormone abscisic acid (ABA) regulates essential processes in plant development and responsiveness to abiotic and biotic stresses. ABA perception triggers a post-translational signaling cascade that elicits the ABA gene regulatory network (GRN), encompassing hundreds of transcription factors (TFs) and thousands of transcribed genes. To further our knowledge of this GRN, we performed an RNA-seq time series experiment consisting of 14 time points in the 16 h following a one-time ABA treatment of 5-week-old Arabidopsis rosettes. During this time course, ABA rapidly changed transcription levels of 7151 genes, which were partitioned into 44 coexpressed modules that carry out diverse biological functions. We integrated our time-series data with publicly available TF-binding site data, motif data, and RNA-seq data of plants inhibited in translation, and predicted (i) which TFs regulate the different coexpression clusters, (ii) which TFs contribute the most to target gene amplitude, (iii) timing of engagement of different TFs in the ABA GRN, and (iv) hierarchical position of TFs and their targets in the multi-tiered ABA GRN. The ABA GRN was found to be highly interconnected and regulated at different amplitudes and timing by a wide variety of TFs, of which the bZIP family was most prominent, and upregulation of genes encompassed more TFs than downregulation. We validated our network models in silico with additional public TF-binding site data and transcription data of selected TF mutants. Finally, using a drought assay we found that the Trihelix TF GT3a is likely an ABA-induced positive regulator of drought tolerance.

Establishment of single-cell transcriptional states during seed germination

Germination involves highly dynamic transcriptional programs as the cells of seeds reactivate and express the functions necessary for establishment in the environment. Individual cell types have distinct roles within the embryo, so must therefore have cell type-specific gene expression and gene regulatory networks. We can better understand how the functions of different cell types are established and contribute to the embryo by determining how cell type-specific transcription begins and changes through germination. Here we describe a temporal analysis of the germinating Arabidopsis thaliana embryo at single-cell resolution. We define the highly dynamic cell type-specific patterns of gene expression and how these relate to changing cellular function as germination progresses. Underlying these are unique gene regulatory networks and transcription factor activity. We unexpectedly discover that most embryo cells transition through the same initial transcriptional state early in germination, even though cell identity has already been established during embryogenesis. Cells later transition to cell type-specific gene expression patterns. Furthermore, our analyses support previous findings that the earliest events leading to the induction of seed germination take place in the vasculature. Overall, our study constitutes a general framework with which to characterize Arabidopsis cell transcriptional states through seed germination, allowing investigation of different genotypes and other plant species whose seed strategies may differ.

The methyl jasmonate-responsive transcription factor SmERF106 promotes tanshinone accumulation in Salvia miltiorrhiza

Understanding the regulatory biosynthesis mechanisms of active compounds in herbs is vital for the preservation and sustainable use of natural medicine resources. Diterpenoids, which play a key role in plant growth and resistance, also serve as practical products for humans. Tanshinone, a class of abietane-type diterpenes unique to the Salvia genus, such as Salvia miltiorrhiza, is an excellent model for studying diterpenoids. In this study, we discovered that a transcription factor, SmERF106, responds to MeJA induction and is located in the nucleus. It exhibits a positive correlation with the expression of SmKSL1 and SmIDI1, which are associated with tanshinone biosynthesis. We performed DNA affinity purification sequencing (DAP-seq) to predict genes that may be transcriptionally regulated by SmERF106. Our cis-elements analysis suggested that SmERF106 might bind to GCC-boxes in the promoters of SmKSL1 and SmIDI1. This indicates that SmKSL1 and SmIDI1 could be potential target genes regulated by SmERF106 in the tanshinone biosynthesis pathway. Their interaction was then demonstrated through a series of in vitro and in vivo binding experiments, including Y1H, EMSA, and Dual-LUC. Overexpression of SmERF106 in the hairy root of S. miltiorrhiza led to a significant increase in tanshinone content and the transcriptional levels of SmKSL1 and SmIDI1. In summary, we found that SmERF106 can activate the transcription of SmKSL1 and SmIDI1 in response to MeJA induction, thereby promoting tanshinone biosynthesis. This discovery provides new insights into the regulatory mechanisms of tanshinones in response to JA and offers a potential gene tool for tanshinone metabolic engineering strategy.

GenomicLinks: deep learning predictions of 3D chromatin interactions in the maize genome

Gene regulation in eukaryotes is partly shaped by the 3D organization of chromatin within the cell nucleus. Distal interactions between cis-regulatory elements and their target genes are widespread, and many causal loci underlying heritable agricultural traits have been mapped to distal non-coding elements. The biology underlying chromatin loop formation in plants is poorly understood. Dissecting the sequence features that mediate distal interactions is an important step toward identifying putative molecular mechanisms. Here, we trained GenomicLinks, a deep learning model, to identify DNA sequence features predictive of 3D chromatin interactions in maize. We found that the presence of binding motifs of specific transcription factor classes, especially bHLH, is predictive of chromatin interaction specificities. Using an in silico mutagenesis approach we show the removal of these motifs from loop anchors leads to reduced interaction probabilities. We were able to validate these predictions with single-cell co-accessibility data from different maize genotypes that harbor natural substitutions in these TF binding motifs. GenomicLinks is currently implemented as an open-source web tool, which should facilitate its wider use in the plant research community.

A NAC transcription factor MNAC3-centered regulatory network negatively modulates rice immunity against blast disease

NAC transcription factors (TFs) are pivotal in plant immunity against diverse pathogens. Here, we report the functional and regulatory network of MNAC3, a novel NAC TF, in rice immunity. MNAC3, a transcriptional activator, negatively modulates rice immunity against blast and bacterial leaf blight diseases and pathogen-associated molecular pattern (PAMP)-triggered immune responses. MNAC3 binds to a CACG cis-element and activates the transcription of immune-negative target genes OsINO80, OsJAZ10, and OsJAZ11. The negative function of MNAC3 in rice immunity depends on its transcription of downstream genes such as OsINO80 and OsJAZ10. MNAC3 interacts with immunity-related OsPP2C41 (a protein phosphatase), ONAC066 (a NAC TF), and OsDjA6 (a DnaJ chaperone). ONAC066 and OsPP2C41 attenuate MNAC3 transcriptional activity, while OsDjA6 promotes it. Phosphorylation of MNAC3 at S163 is critical for its negative functions in rice immunity. OsPP2C41, which plays positive roles in rice blast resistance and chitin-triggered immune responses, dephosphorylates MNAC3, suppressing its transcriptional activity on the target genes OsINO80, OsJAZ10, and OsJAZ11 and promoting the translocation of MNAC3 from nucleus to cytoplasm. These results establish a MNAC3-centered regulatory network in which OsPP2C41 dephosphorylates MNAC3, attenuating its transcriptional activity on downstream immune-negative target genes in rice. Together, these findings deepen our understanding of molecular mechanisms in rice immunity and offer a novel strategy for genetic improvement of rice disease resistance.

Functional Characterization of Accessible Chromatin in Common Wheat

Eukaryotic gene transcription is fine-tuned by precise spatiotemporal interactions between cis-regulatory elements (CREs) and trans-acting factors. However, how CREs individually or coordinated with epigenetic marks function in regulating homoeolog bias expression is still largely unknown in wheat. In this study, through comprehensively characterizing open chromatin coupled with DNA methylation in the seedling and spikelet of common wheat, we observed that differential chromatin openness occurred between the seedling and spikelet, which plays important roles in tissue development through regulating the expression of related genes or through the transcription factor (TF)-centered regulatory network. Moreover, we found that CHH methylation may act as a key determinant affecting the differential binding of TFs, thereby resulting in differential expression of target genes. In addition, we found that sequence variations in MNase hypersensitive sites (MHSs) result in the differential expression of key genes responsible for important agronomic traits. Thus, our study provides new insights into the roles of CREs in regulating tissue or homoeolog bias expression, and controlling important agronomic traits in common wheat. It also provides potential CREs for genetic and epigenetic manipulation toward improving desirable traits for wheat molecule breeding.

Voice from both sides: a molecular dialogue between transcriptional activators and repressors in seed-to-seedling transition and crop adaptation

High-quality seeds provide valuable nutrients to human society and ensure successful seedling establishment. During maturation, seeds accumulate storage compounds that are required to sustain seedling growth during germination. This review focuses on the epigenetic repression of the embryonic and seed maturation programs in seedlings. We begin with an extensive overview of mutants affecting these processes, illustrating the roles of core proteins and accessory components in the epigenetic machinery by comparing mutants at both phenotypic and molecular levels. We highlight how omics assays help uncover target-specific functional specialization and coordination among various epigenetic mechanisms. Furthermore, we provide an in-depth discussion on the Seed dormancy 4 (Sdr4) transcriptional corepressor family, comparing and contrasting their regulation of seed germination in the dicotyledonous species Arabidopsis and two monocotyledonous crops, rice and wheat. Finally, we compare the similarities in the activation and repression of the embryonic and seed maturation programs through a shared set of cis-regulatory elements and discuss the challenges in applying knowledge largely gained in model species to crops.

Syntrichia ruralis: emerging model moss genome reveals a conserved and previously unknown regulator of desiccation in flowering plants

Water scarcity, resulting from climate change, poses a significant threat to ecosystems. Syntrichia ruralis, a dryland desiccation-tolerant moss, provides valuable insights into survival of water-limited conditions. We sequenced the genome of S. ruralis, conducted transcriptomic analyses, and performed comparative genomic and transcriptomic analyses with existing genomes and transcriptomes, including with the close relative S. caninervis. We took a genetic approach to characterize the role of an S. ruralis transcription factor, identified in transcriptomic analyses, in Arabidopsis thaliana. The genome was assembled into 12 chromosomes encompassing 21 169 protein-coding genes. Comparative analysis revealed copy number and transcript abundance differences in known desiccation-associated gene families, and highlighted genome-level variation among species that may reflect adaptation to different habitats. A significant number of abscisic acid (ABA)-responsive genes were found to be negatively regulated by a MYB transcription factor (MYB55) that was upstream of the S. ruralis ortholog of ABA-insensitive 3 (ABI3). We determined that this conserved MYB transcription factor, uncharacterized in Arabidopsis, acts as a negative regulator of an ABA-dependent stress response in Arabidopsis. The new genomic resources from this emerging model moss offer novel insights into how plants regulate their responses to water deprivation.

Identification of plant transcriptional activation domains

Gene expression in Arabidopsis is regulated by more than 1,900 transcription factors (TFs), which have been identified genome-wide by the presence of well-conserved DNA-binding domains. Activator TFs contain activation domains (ADs) that recruit coactivator complexes; however, for nearly all Arabidopsis TFs, we lack knowledge about the presence, location and transcriptional strength of their ADs1. To address this gap, here we use a yeast library approach to experimentally identify Arabidopsis ADs on a proteome-wide scale, and find that more than half of the Arabidopsis TFs contain an AD. We annotate 1,553 ADs, the vast majority of which are, to our knowledge, previously unknown. Using the dataset generated, we develop a neural network to accurately predict ADs and to identify sequence features that are necessary to recruit coactivator complexes. We uncover six distinct combinations of sequence features that result in activation activity, providing a framework to interrogate the subfunctionalization of ADs. Furthermore, we identify ADs in the ancient AUXIN RESPONSE FACTOR family of TFs, revealing that AD positioning is conserved in distinct clades. Our findings provide a deep resource for understanding transcriptional activation, a framework for examining function in intrinsically disordered regions and a predictive model of ADs.

Plant pangenomes for crop improvement, biodiversity and evolution

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

Rewiring gene circuitry for plant improvement

Aspirations for high crop growth and yield, nutritional quality and bioproduction of materials are challenged by climate change and limited adoption of new technologies. Here, we review recent advances in approaches to profile and model gene regulatory activity over developmental and response time in specific cells, which have revealed the basis of variation in plant phenotypes: both redeployment of key regulators to new contexts and their repurposing to control different slates of genes. New synthetic biology tools allow tunable, spatiotemporal regulation of transgenes, while recent gene-editing technologies enable manipulation of the regulation of native genes. Ultimately, understanding how gene circuitry is wired to control form and function across varied plant species, combined with advanced technology to rewire that circuitry, will unlock solutions to our greatest challenges in agriculture, energy and the environment.

Expression of Iron Metabolism Genes Is Potentially Regulated by DOF Transcription Factors in Dendrocalamus latiflorus Leaves

Transcription factors (TFs) are crucial pre-transcriptional regulatory mechanisms that can modulate the expression of downstream genes by binding to their promoter regions. DOF (DNA binding with One Finger) proteins are a unique class of TFs with extensive roles in plant growth and development. Our previous research indicated that iron content varies among bamboo leaves of different colors. However, to our knowledge, genes related to iron metabolism pathways in bamboo species have not yet been studied. Therefore, in the current study, we identified iron metabolism related (IMR) genes in bamboo and determined the TFs that significantly influence them. Among these, DOFs were found to have widespread effects and potentially significant impacts on their expression. We identified specific DOF members in Dendrocalamus latiflorus with binding abilities through homology with Arabidopsis DOF proteins, and established connections between some of these members and IMR genes using RNA-seq data. Additionally, molecular docking confirmed the binding interactions between these DlDOFs and the DOF binding sites in the promoter regions of IMR genes. The co-expression relationship between the two gene sets was further validated using q-PCR experiments. This study paves the way for research into iron metabolism pathways in bamboo and lays the foundation for understanding the role of DOF TFs in D. latiflorus.

Prioritizing Maize Metabolic Gene Regulators through Multi-Omic Network Integration

Elucidating gene regulatory networks is a major area of study within plant systems biology. Phenotypic traits are intricately linked to specific gene expression profiles. These expression patterns arise primarily from regulatory connections between sets of transcription factors (TFs) and their target genes. Here, we integrated 46 co-expression networks, 283 protein-DNA interaction (PDI) assays, and 16 million SNPs used to identify expression quantitative trait loci (eQTL) to construct TF-target networks. In total, we analyzed ∼4.6M interactions to generate four distinct types of TF-target networks: co-expression, PDI, trans -eQTL, and cis -eQTL combined with PDIs. To functionally annotate TFs based on their target genes, we implemented three different network integration strategies. We evaluated the effectiveness of each strategy through TF loss-of function mutant inspection and random network analyses. The multi-network integration allowed us to identify transcriptional regulators of several biological processes. Using the topological properties of the fully integrated network, we identified potential functionally redundant TF paralogs. Our findings retrieved functions previously documented for numerous TFs and revealed novel functions that are crucial for informing the design of future experiments. The approach here-described lays the foundation for the integration of multi-omic datasets in maize and other plant systems.

GRAPHICAL ABSTRACT

Transcription factor binding specificities of the oomycete Phytophthora infestans reflect conserved and divergent evolutionary patterns and predict function

BACKGROUND

Identifying the DNA-binding specificities of transcription factors (TF) is central to understanding gene networks that regulate growth and development. Such knowledge is lacking in oomycetes, a microbial eukaryotic lineage within the stramenopile group. Oomycetes include many important plant and animal pathogens such as the potato and tomato blight agent Phytophthora infestans, which is a tractable model for studying life-stage differentiation within the group.

RESULTS

Mining of the P. infestans genome identified 197 genes encoding proteins belonging to 22 TF families. Their chromosomal distribution was consistent with family expansions through unequal crossing-over, which were likely ancient since each family had similar sizes in most oomycetes. Most TFs exhibited dynamic changes in RNA levels through the P. infestans life cycle. The DNA-binding preferences of 123 proteins were assayed using protein-binding oligonucleotide microarrays, which succeeded with 73 proteins from 14 families. Binding sites predicted for representatives of the families were validated by electrophoretic mobility shift or chromatin immunoprecipitation assays. Consistent with the substantial evolutionary distance of oomycetes from traditional model organisms, only a subset of the DNA-binding preferences resembled those of human or plant orthologs. Phylogenetic analyses of the TF families within P. infestans often discriminated clades with canonical and novel DNA targets. Paralogs with similar binding preferences frequently had distinct patterns of expression suggestive of functional divergence. TFs were predicted to either drive life stage-specific expression or serve as general activators based on the representation of their binding sites within total or developmentally-regulated promoters. This projection was confirmed for one TF using synthetic and mutated promoters fused to reporter genes in vivo.

CONCLUSIONS

We established a large dataset of binding specificities for P. infestans TFs, representing the first in the stramenopile group. This resource provides a basis for understanding transcriptional regulation by linking TFs with their targets, which should help delineate the molecular components of processes such as sporulation and host infection. Our work also yielded insight into TF evolution during the eukaryotic radiation, revealing both functional conservation as well as diversification across kingdoms.

Hyperspectral imaging of liverwort Marchantia polymorpha identifies MpWRKY10 as a key regulator defining Foliar pigmentation patterns

Foliar pigmentation patterns vary among plant species and growth conditions. In this study, we utilize hyperspectral imaging to assess foliar pigmentation in the bryophyte Marchantia polymorpha under nutrient stress and identify associated genetic factors. Using singular value decomposition (SVD) for feature selection, we quantitate color variations induced by deficiencies in phosphate, nitrate, magnesium, calcium, and iron. Pseudo-colored thallus images show that disrupting MpWRKY10 causes irregular pigmentation with auronidin accumulation. Transcriptomic profiling shows that MpWRKY10 regulates phenylpropanoid pathway enzymes and R2R3-MYB transcription factors during phosphate deficiency, with MpMYB14 upregulation preceding pigment accumulation. MpWRKY10 is downregulated in older, pigmented thalli under phosphate deficiency but maintained in young thalli, where it suppresses pigmentation genes. This downregulation is absent in pigmented thalli due to aging. Comparative transcriptome analysis suggests similar WRKY and MYB roles in nutrient response and pigmentation in red-leaf lettuce, alluding to conserved genetic factors controlling foliar pigmentation patterns under nutrient deficiency.

MpANT regulates meristem development in Marchantia polymorpha

Meristems are crucial for organ formation, but our knowledge of their molecular evolution is limited. Here, we show that AINTEGUMENTA (MpANT) in the euANT branch of the APETALA2-like transcription factor family is essential for meristem development in the nonvascular plant Marchantia polymorpha. MpANT is expressed in the thallus meristem. Mpant mutants show defects to maintain meristem identity and undergo meristem duplication, while MpANT overexpressers show ectopic thallus growth. MpANT directly upregulates MpGRAS9 in the SHORT-ROOT (SHR) branch of the GRAS family. In the vascular plant Arabidopsis thaliana, the euANT-branch genes PLETHORAs (AtPLTs) and AtANT are involved in the formation and maintenance of root/shoot apical meristems and lateral organ primordia, and AtPLTs directly target SHR-branch genes. In addition, euANTs bind through a similar DNA-binding motif to many conserved homologous genes in M. polymorpha and A. thaliana. Overall, the euANT pathway has an evolutionarily conserved role in meristem development.

Systematic identification of transcriptional activation domains from non-transcription factor proteins in plants and yeast

Transcription factors can promote gene expression through activation domains. Whole-genome screens have systematically mapped activation domains in transcription factors but not in non-transcription factor proteins (e.g., chromatin regulators and coactivators). To fill this knowledge gap, we employed the activation domain predictor PADDLE to analyze the proteomes of Arabidopsis thaliana and Saccharomyces cerevisiae. We screened 18,000 predicted activation domains from >800 non-transcription factor genes in both species, confirming that 89% of candidate proteins contain active fragments. Our work enables the annotation of hundreds of nuclear proteins as putative coactivators, many of which have never been ascribed any function in plants. Analysis of peptide sequence compositions reveals how the distribution of key amino acids dictates activity. Finally, we validated short, "universal" activation domains with comparable performance to state-of-the-art activation domains used for genome engineering. Our approach enables the genome-wide discovery and annotation of activation domains that can function across diverse eukaryotes.

Computational Reconstruction of the Transcription Factor Regulatory Network Induced by Auxin in Arabidopsis thaliana L

In plant hormone signaling, transcription factor regulatory networks (TFRNs), which link the master transcription factors to the biological processes under their control, remain insufficiently characterized despite their crucial function. Here, we identify a TFRN involved in the response to the key plant hormone auxin and define its impact on auxin-driven biological processes. To reconstruct the TFRN, we developed a three-step procedure, which is based on the integrated analysis of differentially expressed gene lists and a representative collection of transcription factor binding profiles. Its implementation is available as a part of the CisCross web server. With the new method, we distinguished two transcription factor subnetworks. The first operates before auxin treatment and is switched off upon hormone application, the second is switched on by the hormone. Moreover, we characterized the functioning of the auxin-regulated TFRN in control of chlorophyll and lignin biosynthesis, abscisic acid signaling, and ribosome biogenesis.

Beneficial Effects of Phosphite in Arabidopsis thaliana Mediated by Activation of ABA, SA, and JA Biosynthesis and Signaling Pathways

Phosphite (Phi) has gained attention in agriculture due to its biostimulant effect on crops. This molecule has been found to benefit plant performance by providing protection against pathogens, improving yield and fruit quality as well as nutrient and water use efficiency. It is still unclear how Phi enhances plant growth and protects against multiple stresses. It has been hypothesized that Phi acts by directly affecting the pathogens and interacting with the plant cellular components and molecular machinery to elicit defense responses. This study elucidates the mechanisms underlying Phi's beneficial effects on plants, revealing their complex interplay with fundamental signaling pathways. An RNA-seq study of Arabidopsis seedlings under optimal and limiting phosphate conditions helped us unveil Phi's role in promoting plant growth by activating the expression of the genes involved in the biosynthesis and signaling pathways associated with abscisic acid (ABA), salicylic acid (SA), and jasmonic acid (JA). The expression of ABA-related genes, known for their involvement in stress response and development regulation, is triggered by Phi treatment, contributing to enhanced resilience and growth. Simultaneously, the activation of the SA pathway, associated with defense responses, suggests Phi's potential in bolstering plant immunity. Moreover, Phi influences JA biosynthesis and signaling, which are crucial for defense against herbivores and pathogens, thereby strengthening plants' defenses. Our findings reveal a multifaceted mechanism through which Phi benefits Arabidopsis development. Understanding its intricate interplay with key signaling pathways opens avenues for leveraging Phi as a strategic tool to enhance plant resilience, immunity, and growth in agricultural and ecological contexts.

A rapid aureochrome opto-switch enables diatom acclimation to dynamic light

Diatoms often outnumber other eukaryotic algae in the oceans, especially in coastal environments characterized by frequent fluctuations in light intensity. The identities and operational mechanisms of regulatory factors governing diatom acclimation to high light stress remain largely elusive. Here, we identified the AUREO1c protein from the coastal diatom Phaeodactylum tricornutum as a crucial regulator of non-photochemical quenching (NPQ), a photoprotective mechanism that dissipates excess energy as heat. AUREO1c detects light stress using a light-oxygen-voltage (LOV) domain and directly activates the expression of target genes, including LI818 genes that encode NPQ effector proteins, via its bZIP DNA-binding domain. In comparison to a kinase-mediated pathway reported in the freshwater green alga Chlamydomonas reinhardtii, the AUREO1c pathway exhibits a faster response and enables accumulation of LI818 transcript and protein levels to comparable degrees between continuous high-light and fluctuating-light treatments. We propose that the AUREO1c-LI818 pathway contributes to the resilience of diatoms under dynamic light conditions.

Mapping of flumioxazin tolerance in a snap bean diversity panel leads to the discovery of a master genomic region controlling multiple stress resistance genes

Introduction

Effective weed management tools are crucial for maintaining the profitable production of snap bean (Phaseolus vulgaris L.). Preemergence herbicides help the crop to gain a size advantage over the weeds, but the few preemergence herbicides registered in snap bean have poor waterhemp (Amaranthus tuberculatus) control, a major pest in snap bean production. Waterhemp and other difficult-to-control weeds can be managed by flumioxazin, an herbicide that inhibits protoporphyrinogen oxidase (PPO). However, there is limited knowledge about crop tolerance to this herbicide. We aimed to quantify the degree of snap bean tolerance to flumioxazin and explore the underlying mechanisms.

Methods

We investigated the genetic basis of herbicide tolerance using genome-wide association mapping approach utilizing field-collected data from a snap bean diversity panel, combined with gene expression data of cultivars with contrasting response. The response to a preemergence application of flumioxazin was measured by assessing plant population density and shoot biomass variables.

Results

Snap bean tolerance to flumioxazin is associated with a single genomic location in chromosome 02. Tolerance is influenced by several factors, including those that are indirectly affected by seed size/weight and those that directly impact the herbicide's metabolism and protect the cell from reactive oxygen species-induced damage. Transcriptional profiling and co-expression network analysis identified biological pathways likely involved in flumioxazin tolerance, including oxidoreductase processes and programmed cell death. Transcriptional regulation of genes involved in those processes is possibly orchestrated by a transcription factor located in the region identified in the GWAS analysis. Several entries belonging to the Romano class, including Bush Romano 350, Roma II, and Romano Purpiat presented high levels of tolerance in this study. The alleles identified in the diversity panel that condition snap bean tolerance to flumioxazin shed light on a novel mechanism of herbicide tolerance and can be used in crop improvement.

Plant enhancers exhibit both cooperative and additive interactions among their functional elements

Enhancers are cis-regulatory elements that shape gene expression in response to numerous developmental and environmental cues. In animals, several models have been proposed to explain how enhancers integrate the activity of multiple transcription factors. However, it remains largely unclear how plant enhancers integrate transcription factor activity. Here, we use Plant STARR-seq to characterize 3 light-responsive plant enhancers-AB80, Cab-1, and rbcS-E9-derived from genes associated with photosynthesis. Saturation mutagenesis revealed mutations, many of which clustered in short regions, that strongly reduced enhancer activity in the light, in the dark, or in both conditions. When tested in the light, these mutation-sensitive regions did not function on their own; rather, cooperative interactions with other such regions were required for full activity. Epistatic interactions occurred between mutations in adjacent mutation-sensitive regions, and the spacing and order of mutation-sensitive regions in synthetic enhancers affected enhancer activity. In contrast, when tested in the dark, mutation-sensitive regions acted independently and additively in conferring enhancer activity. Taken together, this work demonstrates that plant enhancers show evidence for both cooperative and additive interactions among their functional elements. This knowledge can be harnessed to design strong, condition-specific synthetic enhancers.

The dimorphic diaspore model Aethionema arabicum (Brassicaceae): Distinct molecular and morphological control of responses to parental and germination temperatures

Plants in habitats with unpredictable conditions often have diversified bet-hedging strategies that ensure fitness over a wider range of variable environmental factors. A striking example is the diaspore (seed and fruit) heteromorphism that evolved to maximize species survival in Aethionema arabicum (Brassicaceae) in which external and endogenous triggers allow the production of two distinct diaspores on the same plant. Using this dimorphic diaspore model, we identified contrasting molecular, biophysical, and ecophysiological mechanisms in the germination responses to different temperatures of the mucilaginous seeds (M+ seed morphs), the dispersed indehiscent fruits (IND fruit morphs), and the bare non-mucilaginous M- seeds obtained by pericarp (fruit coat) removal from IND fruits. Large-scale comparative transcriptome and hormone analyses of M+ seeds, IND fruits, and M- seeds provided comprehensive datasets for their distinct thermal responses. Morph-specific differences in co-expressed gene modules in seeds, as well as in seed and pericarp hormone contents, identified a role of the IND pericarp in imposing coat dormancy by generating hypoxia affecting abscisic acid (ABA) sensitivity. This involved expression of morph-specific transcription factors, hypoxia response, and cell wall remodeling genes, as well as altered ABA metabolism, transport, and signaling. Parental temperature affected ABA contents and ABA-related gene expression and altered IND pericarp biomechanical properties. Elucidating the molecular framework underlying the diaspore heteromorphism can provide insight into developmental responses to globally changing temperatures.

LAFL Factors in Seed Development and Phase Transitions

Development is a chain reaction in which one event leads to another until the completion of a life cycle. Phase transitions are milestone events in the cycle of life. LEAFY COTYLEDON1 (LEC1), ABA INSENSITIVE3 (ABI3), FUSCA3 (FUS3), and LEC2 proteins, collectively known as LAFL, are master transcription factors (TFs) regulating seed and other developmental processes. Since the initial characterization of the LAFL genes, more than three decades of active research has generated tremendous amounts of knowledge about these TFs, whose roles in seed development and germination have been comprehensively reviewed. Recent advances in cell biology with genetic and genomic tools have allowed the characterization of the LAFL regulatory networks in previously challenging tissues at a higher throughput and resolution in reference species and crops. In this review, we provide a holistic perspective by integrating advances at the epigenetic, transcriptional, posttranscriptional, and protein levels to exemplify the spatiotemporal regulation of the LAFL networks in Arabidopsis seed development and phase transitions, and we briefly discuss the evolution of these TF networks.

The miR159a-DUO1 module regulates pollen development by modulating auxin biosynthesis and starch metabolism in citrus

Achieving seedlessness in citrus varieties is one of the important objectives of citrus breeding. Male sterility associated with abnormal pollen development is an important factor in seedlessness. However, our understanding of the regulatory mechanism underlying the seedlessness phenotype in citrus is still limited. Here, we determined that the miR159a-DUO1 module played an important role in regulating pollen development in citrus, which further indirectly modulated seed development and fruit size. Both the overexpression of csi-miR159a and the knocking out of DUO1 in Hong Kong kumquat (Fortunella hindsii) resulted in small and seedless fruit phenotypes. Moreover, pollen was severely aborted in both transgenic lines, with arrested pollen mitotic I and abnormal pollen starch metabolism. Through additional cross-pollination experiments, DUO1 was proven to be the key target gene for miR159a to regulate male sterility in citrus. Based on DNA affinity purification sequencing (DAP-seq), RNA-seq, and verified interaction assays, YUC2/YUC6, SS4 and STP8 were identified as downstream target genes of DUO1, those were all positively regulated by DUO1. In transgenic F. hindsii lines, the miR159a-DUO1 module down-regulated the expression of YUC2/YUC6, which decreased indoleacetic acid (IAA) levels and modulated auxin signaling to repress pollen mitotic I. The miR159a-DUO1 module reduced the expression of the starch synthesis gene SS4 and sugar transport gene STP8 to disrupt starch metabolism in pollen. Overall, this work reveals a new mechanism by which the miR159a-DUO1 module regulates pollen development and elucidates the molecular regulatory network underlying male sterility in citrus.

Wheat NAC transcription factor NAC5-1 is a positive regulator of senescence

Wheat (Triticum aestivum L.) is an important source of both calories and protein in global diets, but there is a trade-off between grain yield and protein content. The timing of leaf senescence could mediate this trade-off as it is associated with both declines in photosynthesis and nitrogen remobilization from leaves to grain. NAC transcription factors play key roles in regulating senescence timing. In rice, OsNAC5 expression is correlated with increased protein content and upregulated in senescing leaves, but the role of the wheat ortholog in senescence had not been characterized. We verified that NAC5-1 is the ortholog of OsNAC5 and that it is expressed in senescing flag leaves in wheat. To characterize NAC5-1, we combined missense mutations in NAC5-A1 and NAC5-B1 from a TILLING mutant population and overexpressed NAC5-A1 in wheat. Mutation in NAC5-1 was associated with delayed onset of flag leaf senescence, while overexpression of NAC5-A1 was associated with slightly earlier onset of leaf senescence. DAP-seq was performed to locate transcription factor binding sites of NAC5-1. Analysis of DAP-seq and comparison with other studies identified putative downstream target genes of NAC5-1 which could be associated with senescence. This work showed that NAC5-1 is a positive transcriptional regulator of leaf senescence in wheat. Further research is needed to test the effect of NAC5-1 on yield and protein content in field trials, to assess the potential to exploit this senescence regulator to develop high-yielding wheat while maintaining grain protein content.

Optimizing data integration improves gene regulatory network inference in Arabidopsis thaliana

MOTIVATIONS

Gene regulatory networks (GRNs) are traditionally inferred from gene expression profiles monitoring a specific condition or treatment. In the last decade, integrative strategies have successfully emerged to guide GRN inference from gene expression with complementary prior data. However, datasets used as prior information and validation gold standards are often related and limited to a subset of genes. This lack of complete and independent evaluation calls for new criteria to robustly estimate the optimal intensity of prior data integration in the inference process.

RESULTS

We address this issue for two regression-based GRN inference models, a weighted random forest (weigthedRF) and a generalized linear model estimated under a weighted LASSO penalty with stability selection (weightedLASSO). These approaches are applied to data from the root response to nitrate induction in Arabidopsis thaliana. For each gene, we measure how the integration of transcription factor binding motifs influences model prediction. We propose a new approach, DIOgene, that uses model prediction error and a simulated null hypothesis in order to optimize data integration strength in a hypothesis-driven, gene-specific manner. This integration scheme reveals a strong diversity of optimal integration intensities between genes, and offers good performance in minimizing prediction error as well as retrieving experimental interactions. Experimental results show that DIOgene compares favorably against state-of-the-art approaches and allows to recover master regulators of nitrate induction.

AVAILABILITY AND IMPLEMENTATION

The R code and notebooks demonstrating the use of the proposed approaches are available in the repository https://github.com/OceaneCsn/integrative_GRN_N_induction.

AraLeTA: An Arabidopsis leaf expression atlas across diurnal and developmental scales

Mature plant leaves are a composite of distinct cell types, including epidermal, mesophyll, and vascular cells. Notably, the proportion of these cells and the relative transcript concentrations within different cell types may change over time. While gene expression data at a single-cell level can provide cell-type-specific expression values, it is often too expensive to obtain these data for high-resolution time series. Although bulk RNA-seq can be performed in a high-resolution time series, RNA-seq using whole leaves measures average gene expression values across all cell types in each sample. In this study, we combined single-cell RNA-seq data with time-series data from whole leaves to assemble an atlas of cell-type-specific changes in gene expression over time for Arabidopsis (Arabidopsis thaliana). We inferred how the relative transcript concentrations of different cell types vary across diurnal and developmental timescales. Importantly, this analysis revealed 3 subgroups of mesophyll cells with distinct temporal profiles of expression. Finally, we developed tissue-specific gene networks that form a community resource: an Arabidopsis Leaf Time-dependent Atlas (AraLeTa). This allows users to extract gene networks that are confirmed by transcription factor-binding data and specific to certain cell types at certain times of day and at certain developmental stages. AraLeTa is available at https://regulatorynet.shinyapps.io/araleta/.

A CUC1/auxin genetic module links cell polarity to patterned tissue growth and leaf shape diversity in crucifer plants

How tissue-level information encoded by fields of regulatory gene activity is translated into the patterns of cell polarity and growth that generate the diverse shapes of different species remains poorly understood. Here, we investigate this problem in the case of leaf shape differences between Arabidopsis thaliana, which has simple leaves, and its relative Cardamine hirsuta that has complex leaves divided into leaflets. We show that patterned expression of the transcription factor CUP-SHAPED COTYLEDON1 in C. hirsuta (ChCUC1) is a key determinant of leaf shape differences between the two species. Through inducible genetic perturbations, time-lapse imaging of growth, and computational modeling, we find that ChCUC1 provides instructive input into auxin-based leaf margin patterning. This input arises via transcriptional regulation of multiple auxin homeostasis components, including direct activation of WAG kinases that are known to regulate the polarity of PIN-FORMED auxin transporters. Thus, we have uncovered a mechanism that bridges biological scales by linking spatially distributed and species-specific transcription factor expression to cell-level polarity and growth, to shape diverse leaf forms.

Genome-wide transcript expression analysis reveals major chickpea and lentil genes associated with plant branching

The search for elite cultivars with better architecture has been a demand by farmers of the chickpea and lentil crops, which aims to systematize their mechanized planting and harvesting on a large scale. Therefore, the identification of genes associated with the regulation of the branching and architecture of these plants has currently gained great importance. Herein, this work aimed to gain insight into transcriptomic changes of two contrasting chickpea and lentil cultivars in terms of branching pattern (little versus highly branched cultivars). In addition, we aimed to identify candidate genes involved in the regulation of shoot branching that could be used as future targets for molecular breeding. The axillary and apical buds of chickpea cultivars Blanco lechoso and FLIP07-318C, and lentil cultivars Castellana and Campisi, considered as little and highly branched, respectively, were harvested. A total of 1,624 and 2,512 transcripts were identified as differentially expressed among different tissues and contrasting cultivars of chickpea and lentil, respectively. Several gene categories were significantly modulated such as cell cycle, DNA transcription, energy metabolism, hormonal biosynthesis and signaling, proteolysis, and vegetative development between apical and axillary tissues and contrasting cultivars of chickpea and lentil. Based on differential expression and branching-associated biological function, ten chickpea genes and seven lentil genes were considered the main players involved in differentially regulating the plant branching between contrasting cultivars. These collective data putatively revealed the general mechanism and high-effect genes associated with the regulation of branching in chickpea and lentil, which are potential targets for manipulation through genome editing and transgenesis aiming to improve plant architecture.

PLETHORA transcription factors promote early embryo development through induction of meristematic potential

Plants are dependent on divisions of stem cells to establish cell lineages required for growth. During embryogenesis, early division products are considered to be stem cells, whereas during post-embryonic development, stem cells are present in meristems at the root and shoot apex. PLETHORA/AINTEGUMENTA-LIKE (PLT/AIL) transcription factors are regulators of post-embryonic meristem function and are required to maintain stem cell pools. Despite the parallels between embryonic and post-embryonic stem cells, the role of PLTs during early embryogenesis has not been thoroughly investigated. Here, we demonstrate that the PLT regulome in the zygote, and apical and basal cells is in strong congruence with that of post-embryonic meristematic cells. We reveal that out of all six PLTs, only PLT2 and PLT4/BABY BOOM (BBM) are expressed in the zygote, and that these two factors are essential for progression of embryogenesis beyond the zygote stage and first divisions. Finally, we show that other PLTs can rescue plt2 bbm defects when expressed from the PLT2 and BBM promoters, establishing upstream regulation as a key factor in early embryogenesis. Our data indicate that generic PLT factors facilitate early embryo development in Arabidopsis by induction of meristematic potential.

Strigolactone induces D14-dependent large-scale changes in gene expression requiring SWI/SNF chromatin remodellers

Strigolactones (SL) function as plant hormones in control of multiple aspects of plant development, mostly via the regulation of gene expression. Immediate early-gene regulation by SL remains unexplored due to difficulty in dissecting early from late gene expression responses to SL. We used synthetic SL, rac-GR24 treatment of protoplasts and RNA-seq to explore early SL-induced changes in gene expression over time (5-180 minutes) and discovered rapid, dynamic and SL receptor D14-dependent regulation of gene expression in response to rac-GR24. Importantly, we discovered a significant dependence of SL signalling on chromatin remodelling processes, as the induction of a key SL-induced transcription factor BRANCHED1 requires the SWI/SNF chromatin remodelling ATPase SPLAYED (SYD) and leads to upregulation of a homologue SWI/SNF ATPase BRAHMA. ATAC-seq profiling of genome-wide changes in chromatin accessibility in response to rac-GR24 identified large-scale changes, with over 1400 differentially accessible regions. These changes in chromatin accessibility often precede transcriptional changes and are likely to harbour SL cis-regulatory elements. Importantly, we discovered that this early and extensive modification of the chromatin landscape also requires SYD. This study, therefore, provides evidence that SL signalling requires regulation of chromatin accessibility, and it identifies genomic locations harbouring likely SL cis-regulatory sequences.

Recent advances in exploring transcriptional regulatory landscape of crops

Crop breeding entails developing and selecting plant varieties with improved agronomic traits. Modern molecular techniques, such as genome editing, enable more efficient manipulation of plant phenotype by altering the expression of particular regulatory or functional genes. Hence, it is essential to thoroughly comprehend the transcriptional regulatory mechanisms that underpin these traits. In the multi-omics era, a large amount of omics data has been generated for diverse crop species, including genomics, epigenomics, transcriptomics, proteomics, and single-cell omics. The abundant data resources and the emergence of advanced computational tools offer unprecedented opportunities for obtaining a holistic view and profound understanding of the regulatory processes linked to desirable traits. This review focuses on integrated network approaches that utilize multi-omics data to investigate gene expression regulation. Various types of regulatory networks and their inference methods are discussed, focusing on recent advancements in crop plants. The integration of multi-omics data has been proven to be crucial for the construction of high-confidence regulatory networks. With the refinement of these methodologies, they will significantly enhance crop breeding efforts and contribute to global food security.

The transcription factor bZIP44 cooperates with MYB10 and MYB72 to regulate the response of Arabidopsis thaliana to iron deficiency stress

Nicotianamine (NA) plays a crucial role in transporting metal ions, including iron (Fe), in plants; therefore, NICOTIANAMINE SYNTHASE (NAS) genes, which control NA synthesis, are tightly regulated at the transcriptional level. However, the transcriptional regulatory mechanisms of NAS genes require further investigations. In this study, we determined the role of bZIP44 in mediating plant response to Fe deficiency stress by conducting transformation experiments and assays. bZIP44 positively regulated the response of Arabidopsis to Fe deficiency stress by interacting with MYB10 and MYB72 to enhance their abilities to bind at NAS2 and NAS4 promoters, thereby increasing NAS2 and NAS4 transcriptional levels and promote NA synthesis. In summary, the transcription activities of bZIP44, MYB10, and MYB72 were induced in response to Fe deficiency stress, which enhanced the interaction between bZIP44 and MYB10 or MYB72 proteins, synergistically activated the transcriptional activity of NAS2 and NAS4, promoted NA synthesis, and improved Fe transport, thereby enhancing plant tolerance to Fe deficiency stress.

Regulatory mechanisms used by ZmMYB39 to enhance drought tolerance in maize (Zea mays) seedlings

Drought is a significant abiotic stressor that limits maize (Zea mays L.) growth and development. Thus, enhancing drought tolerance is critical for promoting maize production. Our findings demonstrated that ZmMYB39 is an MYB transcription factor with transcriptional activation activity. Drought stress experiments involving ZmMYB39 overexpression and knockout lines indicated that ZmMYB39 positively regulated drought stress tolerance in maize. DAP-Seq, EMSA, dual-LUC, and RT-qPCR provided initial insights into the molecular regulatory mechanisms by which ZmMYB39 enhances drought tolerance in maize. ZmMYB39 directly promoted the expression of ZmP5CS1, ZmPOX1, ZmSOD2, ZmRD22, ZmNAC49, and ZmDREB2A, which are involved in stress resistance. ZmMYB39 enhanced drought tolerance by interacting with and promoting the expression of ZmFNR1, ZmHSP20, and ZmDOF6. Our study offers a theoretical basis for understanding the molecular regulatory networks involved in maize drought stress response. Furthermore, ZmMYB39 serves as a valuable genetic resource for breeding drought-resistant maize.

Genomic analyses provide insights into sex differentiation of tetraploid strawberry (Fragaria moupinensis)

The strawberry genus, Fragaria, exhibits a wide range of sexual systems and natural ploidy variation. Nearly, all polyploid strawberry species exhibit separate sexes (dioecy). Research has identified the sex-determining sequences as roughly conserved but with repeatedly changed genomic locations across octoploid strawberries. However, it remains unclear whether tetraploid wild strawberries evolved dioecy independently or shared a common origin with octoploid strawberries. In this study, we investigated the sex determinants of F. moupinensis, a dioecious plant with heterogametic females (ZW). Utilizing a combination of haplotype-resolved genome sequencing of the female F. moupinensis, k-mer-based and coverage-based genome-wide association studies (GWAS), and transcriptomic analysis, we discovered a non-recombining, approximately 33.6 kb W-specific region on chromosome 2a. Within this region, only one candidate sex-determining gene (FmoAFT) was identified. Furthermore, an extensive resequencing of the entire Fragaria genus indicated that the W-specific region displays conservative female specificity across all tetraploid species. This observation suggests that dioecy evolved independently in tetraploid and octoploid strawberries. Moreover, employing virus-induced gene silencing (VIGS), we knocked down the expression of the FmoAFT homologue transcript in cultivated strawberries, revealing its potential role in promoting female functions during early carpel development. We also applied DNA affinity purification sequencing (DAP-seq) and yeast one-hybrid assays to identify potential direct targets of FmoAFT. These insights shed new light on the genetic basis and evolutionary history of sex determination in strawberries, thereby facilitating the formulation of strategies to manipulate sex determination in breeding programs.

Insights into multilevel spatial regulation within the root stem cell niche

All differentiated root cells derive from stem cells spatially organized within the stem cell niche (SCN), a microenvironment located within the root tip. Here, we compiled recent advances in the understanding of how the SCN drives the establishment and maintenance of cell types. The quiescent center (QC) is widely recognized as the primary driver of cell fate determination, but it is recently considered a convergence center of multiple signals. Cell identity of the cortex endodermis initials is mainly driven by the regulatory feedback loops between transcription factors (TFs), acting as mobile signals between neighboring cells, including the QC. As exemplified in the vascular initials, the precise spatial expression of these regulatory TFs is connected with a dynamic hormonal interplay. Thus, stem cell maintenance and cell differentiation are regulated by a plethora of signals forming a complex, multilevel regulatory network. Integrating the transcriptional and post-translational regulations, protein-protein interactions, and mobile signals into models will be fundamental for the comprehensive understanding of SCN maintenance and differentiation.

Transcription factor OsWRKY72 controls rice leaf angle by regulating LAZY1-mediated shoot gravitropism

Leaf angle is a major trait of ideal architecture, which is considered to influence rice (Oryza sativa) cultivation and grain yield. Although a few mutants with altered rice leaf inclination angles have been reported, the underlying molecular mechanism remains unclear. In this study, we showed that a WRKY transcription factor gene, OsWRKY72, was highly expressed in the leaf sheath and lamina joint. Phenotypic analyses showed that oswrky72 mutants had smaller leaf angles than the wild type, while OsWRKY72 overexpression lines exhibited an increased leaf angle. This observation suggests that OsWRKY72 functions as a positive regulator, promoting the enlargement of the leaf angle. Our bioinformatics analysis identified LAZY1 as the downstream gene of OsWRKY72. Electrophoretic mobility shift assays and dual-luciferase analysis revealed that OsWRKY72 directly inhibited LAZY1 by binding to its promoter. Moreover, knocking out OsWRKY72 enhanced shoot gravitropism, which contrasted with the phenotype of lazy1 plants. These results imply that OsWRKY72 regulates the leaf angle through gravitropism by reducing the expression of LAZY1. In addition, OsWRKY72 could directly regulate the expression of other leaf angle-related genes such as FLOWERING LOCUS T-LIKE 12 (OsFTL12) and WALL-ASSOCIATED KINASE 11 (OsWAK11). Our study indicates that OsWRKY72 contributes positively to the expansion of the leaf angle by interfering with shoot gravitropism in rice.

ETHYLENE RESPONSE FACTOR 070 inhibits flowering in Pak-choi by indirectly impairing BcLEAFY expression

APETALA2/ethylene responsive factors respond to ethylene and participate in many biological and physiological processes, such as plant morphogenesis, stress resistance, and hormone signal transduction. Ethylene responsive factor 070 (BcERF070) is important in flowering. However, the underlying molecular mechanisms of BcERF070 in floral transition in response to ethylene signaling have not been fully characterized. Herein, we explored the function of BcERF070 in Pak-choi [Brassica campestris (syn. Brassica rapa) ssp. chinensis]. Ethylene treatment induced BcERF070 expression and delayed flowering in Pak-choi. Silencing of BcERF070 induced flowering in Pak-choi. BcERF070 interacted with major latex protein-like 328 (BcMLP328), which forms a complex with helix-loop-helix protein 30 (BcbHLH30) to enhance the transcriptional activity of BcbHLH30 on LEAFY (BcLFY), ultimately promoting flowering. However, BcERF070 impaired the BcMLP328-BcbHLH30 complex activation of LEAFY (BcLFY), ultimately inhibiting flowering in Pak-choi. BcERF070 directly promoted the expression of the flowering inhibitor gene B-box 29 (BcBBX29) and delayed flowering by reducing FLOWERING LOCUS T (BcFT) expression. These results suggest that BcERF070 mediates ethylene-reduced flowering by impairing the BcMLP328-BcbHLH30 complex activation of BcLFY and by directly promoting the gene expression of the flowering inhibition factor BcBBX29 to repress BcFT expression. The findings contribute to understanding the molecular mechanisms underlying floral transition in response to ethylene in plants.

Identification of Eleutherococcus senticosus NAC transcription factors and their mechanisms in mediating DNA methylation of EsFPS, EsSS, and EsSE promoters to regulate saponin synthesis

BACKGROUND

The formation of pharmacologically active components in medicinal plants is significantly impacted by DNA methylation. However, the exact mechanisms through which DNA methylation regulates secondary metabolism remain incompletely understood. Research in model species has demonstrated that DNA methylation at the transcription factor binding site within functional gene promoters can impact the binding of transcription factors to target DNA, subsequently influencing gene expression. These findings suggest that the interaction between transcription factors and target DNA could be a significant mechanism through which DNA methylation regulates secondary metabolism in medicinal plants.

RESULTS

This research conducted a comprehensive analysis of the NAC family in E. senticosus, encompassing genome-wide characterization and functional analysis. A total of 117 EsNAC genes were identified and phylogenetically divided into 15 subfamilies. Tandem duplications and chromosome segment duplications were found to be the primary replication modes of these genes. Motif 2 was identified as the core conserved motif of the genes, and the cis-acting elements, gene structures, and expression patterns of each EsNAC gene were different. EsJUB1, EsNAC047, EsNAC098, and EsNAC005 were significantly associated with the DNA methylation ratio in E. senticosus. These four genes were located in the nucleus or cytoplasm and exhibited transcriptional self-activation activity. DNA methylation in EsFPS, EsSS, and EsSE promoters significantly reduced their activity. The methyl groups added to cytosine directly hindered the binding of the promoters to EsJUB1, EsNAC047, EsNAC098, and EsNAC005 and altered the expression of EsFPS, EsSS, and EsSE genes, eventually leading to changes in saponin synthesis in E. senticosus.

CONCLUSIONS

NAC transcription factors that are hindered from binding by methylated DNA are found in E. senticosus. The incapacity of these NACs to bind to the promoter of the methylated saponin synthase gene leads to subsequent alterations in gene expression and saponin synthesis. This research is the initial evidence showcasing the involvement of EsNAC in governing the impact of DNA methylation on saponin production in E. senticosus.