An atlas of the tomato epigenome reveals that KRYPTONITE shapes TAD-like boundaries through the control of H3K9ac distribution.

In recent years, the exploration of genome three-dimensional (3D) conformation has yielded profound insights into the regulation of gene expression and cellular functions in both animals and plants. While animals exhibit a characteristic genome topology defined by topologically associating domains (TADs), plants display similar features with a more diverse conformation across species. Employing advanced high-throughput sequencing and microscopy techniques, we investigated the landscape of 26 histone modifications and RNA polymerase II distribution in tomato (Solanum lycopersicum). Our study unveiled a rich and nuanced epigenetic landscape, shedding light on distinct chromatin states associated with heterochromatin formation and gene silencing. Moreover, we elucidated the intricate interplay between these chromatin states and the overall topology of the genome. Employing a genetic approach, we delved into the role of the histone modification H3K9ac in genome topology. Notably, our investigation revealed that the ectopic deposition of this chromatin mark triggered a reorganization of the 3D chromatin structure, defining different TAD-like borders. Our work emphasizes the critical role of H3K9ac in shaping the topology of the tomato genome, providing valuable insights into the epigenetic landscape of this agriculturally significant crop species.

Epigenetic control during root development and symbiosis.

The roots of plants play multiples functions that are essential for growth and development, including anchoring to the soil and water and nutrient acquisition. These underground organs exhibit the plasticity to modify their root system architecture in response to environmental cues allowing adaptation to change in water and nutrient availability. In addition, roots enter in mutualistic interactions with soil microorganisms, e.g. the root nodule symbiosis established between a limited group of plants and nitrogen fixing soil bacteria and the arbuscular mycorrhiza symbiosis involving most land plants and fungi of the Glomeromycetes phylum. In the past 20 years, genetic approaches allowed the identification and functional characterization of genes required for the specific programs of root development, root nodule and arbuscular mycorrhiza symbioses. These genetic studies provided evidence that the program of root nodule symbiosis recruited components of the arbuscular mycorrhiza symbiosis and the root developmental programs. The execution of these programs is strongly influenced by epigenetic changes -DNA methylation and histone post-translational modifications- that alter chromatin conformation modifying the expression of key genes. In this review, we summarize recent advances that highlighted how DNA methylation and histone post-translational modifications, as well as chromatin remodeling factors and long non-coding RNAs, shape the root system architecture and allow the successful establishment of both root nodule and arbuscular mycorrhiza symbioses. We anticipate that the analysis of dynamic epigenetic changes and chromatin 3D structure in specific single-cells or tissue types of root organs will illuminate our understanding of how root developmental and symbiotic programs are orchestrated, opening exciting questions and new perspectives to modulate agronomical and ecological traits linked to nutrient acquisition.

A lateral organ boundaries domain transcription factor acts downstream of the auxin response factor 2 to control nodulation and root architecture in Medicago truncatula.

Legume plants develop two types of root postembryonic organs, lateral roots and symbiotic nodules, using shared regulatory components. The module composed by the microRNA390, the Trans-Acting SIRNA3 (TAS3) RNA and the Auxin Response Factors (ARF)2, ARF3, and ARF4 (miR390/TAS3/ARFs) mediates the control of both lateral roots and symbiotic nodules in legumes. Here, a transcriptomic approach identified a member of the Lateral Organ Boundaries Domain (LBD) family of transcription factors in Medicago truncatula, designated MtLBD17/29a, which is regulated by the miR390/TAS3/ARFs module. ChIP-PCR experiments evidenced that MtARF2 binds to an Auxin Response Element present in the MtLBD17/29a promoter. MtLBD17/29a is expressed in root meristems, lateral root primordia, and noninfected cells of symbiotic nodules. Knockdown of MtLBD17/29a reduced the length of primary and lateral roots and enhanced lateral root formation, whereas overexpression of MtLBD17/29a produced the opposite phenotype. Interestingly, both knockdown and overexpression of MtLBD17/29a reduced nodule number and infection events and impaired the induction of the symbiotic genes Nodulation Signaling Pathway (NSP) 1 and 2. Our results demonstrate that MtLBD17/29a is regulated by the miR390/TAS3/ARFs module and a direct target of MtARF2, revealing a new lateral root regulatory hub recruited by legumes to act in the root nodule symbiotic program.

Long noncoding RNA-mediated epigenetic regulation of auxin-related genes controls shade avoidance syndrome in Arabidopsis.

The long noncoding RNA (lncRNA) AUXIN-REGULATED PROMOTER LOOP (APOLO) recognizes a subset of target loci across the Arabidopsis thaliana genome by forming RNA-DNA hybrids (R-loops) and modulating local three-dimensional chromatin conformation. Here, we show that APOLO regulates shade avoidance syndrome by dynamically modulating expression of key factors. In response to far-red (FR) light, expression of APOLO anti-correlates with that of its target BRANCHED1 (BRC1), a master regulator of shoot branching in Arabidopsis thaliana. APOLO deregulation results in BRC1 transcriptional repression and an increase in the number of branches. Accumulation of APOLO transcription fine-tunes the formation of a repressive chromatin loop encompassing the BRC1 promoter, which normally occurs only in leaves and in a late response to far-red light treatment in axillary buds. In addition, our data reveal that APOLO participates in leaf hyponasty, in agreement with its previously reported role in the control of auxin homeostasis through direct modulation of auxin synthesis gene YUCCA2, and auxin efflux genes PID and WAG2. We show that direct application of APOLO RNA to leaves results in a rapid increase in auxin signaling that is associated with changes in the plant response to far-red light. Collectively, our data support the view that lncRNAs coordinate shade avoidance syndrome in A. thaliana, and reveal their potential as exogenous bioactive molecules. Deploying exogenous RNAs that modulate plant-environment interactions may therefore become a new tool for sustainable agriculture.

Dynamics of ribosome composition and ribosomal protein phosphorylation in immune signaling in Arabidopsis thaliana.

In plants, the detection of microbe-associated molecular patterns (MAMPs) induces primary innate immunity by the activation of mitogen-activated protein kinases (MAPKs). We show here that the MAMP-activated MAPK MPK6 not only modulates defense through transcriptional regulation but also via the ribosomal protein translation machinery. To understand the effects of MPK6 on ribosomes and their constituent ribosomal proteins (RPs), polysomes, monosomes and the phosphorylation status of the RPs, MAMP-treated WT and mpk6 mutant plants were analysed. MAMP-activation induced rapid changes in RP composition of monosomes, polysomes and in the 60S ribosomal subunit in an MPK6-specific manner. Phosphoproteome analysis showed that MAMP-activation of MPK6 regulates the phosphorylation status of the P-stalk ribosomal proteins by phosphorylation of RPP0 and the concomitant dephosphorylation of RPP1 and RPP2. These events coincide with a significant decrease in the abundance of ribosome-bound RPP0s, RPP1s and RPP3s in polysomes. The P-stalk is essential in regulating protein translation by recruiting elongation factors. Accordingly, we found that RPP0C mutant plants are compromised in basal resistance to Pseudomonas syringae infection. These data suggest that MAMP-induced defense also involves MPK6-induced regulation of P-stalk proteins, highlighting a new role of ribosomal regulation in plant innate immunity.

Pathogen-induced m6A dynamics affect plant immunity.

Posttranscriptional regulation of mRNA mediated by methylation at the N6 position of adenine (N6-methyladenosine [m6A]) has profound effects on transcriptome regulation in plants. Focused studies across eukaryotes offer glimpses into the processes governed by m6A throughout developmental and disease states. However, we lack an understanding of the dynamics and the regulatory potential of m6A during biotic stress in plants. Here, we provide a comprehensive look into the effects of m6A on both the short-term and long-term responses to pathogen signaling in Arabidopsis (Arabidopsis thaliana). We demonstrate that m6A-deficient plants are more resistant to bacterial and fungal pathogen infections and have altered immune responses. Furthermore, m6A deposition is specifically coordinated on transcripts involved in defense and immunity prior to and proceeding the pathogen signal flagellin. Consequently, the dynamic modulation of m6A on specific stress-responsive transcripts is correlated with changes in abundance and cleavage of these transcripts. Overall, we show that the m6A methylome is regulated prior to and during simulated and active pathogen stress and functions in the coordination and balancing of normal growth and pathogen responses.

The long intergenic noncoding RNA ARES modulates root architecture in Arabidopsis.

Long noncoding RNAs (lncRNAs) have emerged as important regulators of gene expression in plants. They have been linked to a wide range of molecular mechanisms, including epigenetics, miRNA activity, RNA processing and translation, and protein localization or stability. In Arabidopsis, characterized lncRNAs have been implicated in several physiological contexts, including plant development and the response to the environment. Here we searched for lncRNA loci located nearby key genes involved in root development and identified the lncRNA ARES (AUXIN REGULATOR ELEMENT DOWNSTREAM SOLITARYROOT) downstream of the lateral root master gene IAA14/SOLITARYROOT (SLR). Although ARES and IAA14 are co-regulated during development, the knockdown and knockout of ARES did not affect IAA14 expression. However, in response to exogenous auxin, ARES knockdown impairs the induction of its other neighboring gene encoding the transcription factor NF-YB3. Furthermore, knockdown/out of ARES results in a root developmental phenotype in control conditions. Accordingly, a transcriptomic analysis revealed that a subset of ARF7-dependent genes is deregulated. Altogether, our results hint at the lncRNA ARES as a novel regulator of the auxin response governing lateral root development, likely by modulating gene expression in trans.

Transcription factors KANADI 1, MYB DOMAIN PROTEIN 44, and PHYTOCHROME INTERACTING FACTOR 4 regulate long intergenic noncoding RNAs expressed in Arabidopsis roots.

Thousands of long intergenic noncoding RNAs (lincRNAs) have been identified in plant genomes. While some lincRNAs have been characterized as important regulators in different biological processes, little is known about the transcriptional regulation for most plant lincRNAs. Through the integration of 8 annotation resources, we defined 6,599 high-confidence lincRNA loci in Arabidopsis (Arabidopsis thaliana). For lincRNAs belonging to different evolutionary age categories, we identified major differences in sequence and chromatin features, as well as in the level of conservation and purifying selection acting during evolution. Spatiotemporal gene expression profiles combined with transcription factor (TF) chromatin immunoprecipitation (ChIP) data were used to construct a TF-lincRNA regulatory network containing 2,659 lincRNAs and 15,686 interactions. We found that properties characterizing lincRNA expression, conservation, and regulation differ between plants and animals. Experimental validation confirmed the role of 3 TFs, KANADI 1, MYB DOMAIN PROTEIN 44, and PHYTOCHROME INTERACTING FACTOR 4, as key regulators controlling root-specific lincRNA expression, demonstrating the predictive power of our network. Furthermore, we identified 58 lincRNAs, regulated by these TFs, showing strong root cell type-specific expression or chromatin accessibility, which are linked with genome-wide association studies genetic associations related to root system development and growth. The multilevel genome-wide characterization covering chromatin state information, promoter conservation, and chromatin immunoprecipitation-based TF binding, for all detectable lincRNAs across 769 expression samples, permits rapidly defining the biological context and relevance of Arabidopsis lincRNAs through regulatory networks.

Synergistic action of the Arabidopsis spliceosome components PRP39a and SmD1b in promoting posttranscriptional transgene silencing.

Besides regulating splicing, the conserved spliceosome component SmD1 (Small nuclear ribonucleoprotein D1)b promotes posttranscriptional silencing of sense transgenes (S-PTGS [post-transcriptional genesilencing]). Here, we show that the conserved spliceosome component PRP39 (Pre-mRNA-processing factor 39)a also plays a role in S-PTGS in Arabidopsis thaliana. However, PRP39a and SmD1b actions appear distinct in both splicing and S-PTGS. Indeed, RNAseq-based analysis of expression level and alternative splicing in prp39a and smd1b mutants identified different sets of deregulated transcripts and noncoding RNAs. Moreover, double mutant analyses involving prp39a or smd1b and RNA quality control (RQC) mutants revealed distinct genetic interactions for SmD1b and PRP39a with nuclear RQC machineries, suggesting nonredundant roles in the RQC/PTGS interplay. Supporting this hypothesis, a prp39a smd1b double mutant exhibited enhanced suppression of S-PTGS compared to the single mutants. Because the prp39a and smd1b mutants (i) showed no major changes in the expression of PTGS or RQC components or in small RNA production and (ii) do not alter PTGS triggered by inverted-repeat transgenes directly producing dsRNA (IR-PTGS), PRP39a, and SmD1b appear to synergistically promote a step specific to S-PTGS. We propose that, independently from their specific roles in splicing, PRP39a and SmD1b limit 3'-to-5' and/or 5'-to-3' degradation of transgene-derived aberrant RNAs in the nucleus, thus favoring the export of aberrant RNAs to the cytoplasm where their conversion into double-stranded RNA (dsRNA) initiates S-PTGS.

Long non-coding RNAs reveal new regulatory mechanisms controlling gene expression.

A plethora of non-coding RNAs have been found in eukaryotes, notably with the advent of modern sequencing technologies to analyze the transcriptome. Apart from the well-known housekeeping RNA genes (such as the ribosomal RNA or the transfer RNA), many thousands of transcripts detected are not evidently linked to a protein-coding gene. These, so called non-coding RNAs, may code for crucial regulators of gene expression, the small si/miRNAs, for small peptides (translated under specific conditions) or may act as long RNA molecules (antisense, intronic or intergenic long non-coding RNAs or lncRNAs). The lncRNAs interact with members of multiple machineries involved in gene regulation. In this review, we discussed about how plant lncRNAs permitted to discover new regulatory mechanisms acting in epigenetic control, chromatin 3D structure and alternative splicing. These novel regulations diversified the expression patterns and protein variants of target protein-coding genes and are an important element of the response of plants to environmental stresses and their adaptation to changing conditions.

Une pléthore d'ARN non codants a été découverte chez les eucaryotes, notamment avec l'avènement des technologies modernes de séquençage pour analyser le transcriptome. En dehors des gènes ARN domestiques bien connus (tels que l'ARN ribosomal ou l'ARN de transfert), plusieurs milliers de transcrits détectés ne sont pas manifestement liés à un gène codant pour une protéine. Ces ARN, appelés ARN non codants, peuvent coder pour des régulateurs cruciaux de l'expression des gènes, les petits si/miARN, pour de petits peptides (traduits dans des conditions spécifiques) ou peuvent agir comme de longues molécules d'ARN (ARN non codants antisens, introniques ou intergéniques ou lncRNA). Les lncRNAs interagissent avec les membres des multiples machineries impliquées dans la régulation des gènes. Dans cette revue, nous avons discuté de la façon dont les lncRNAs des plantes ont permis de découvrir de nouveaux mécanismes de régulation agissant dans le contrôle épigénétique, la structure 3D de la chromatine et l'épissage alternatif. Ces nouvelles régulations ont diversifié les profils d'expression des gènes codants pour les protéines cibles et constituent un élément important de la réponse des plantes aux stress environnementaux et de leur adaptation aux conditions changeantes.

HSFA1a modulates plant heat stress responses and alters the 3D chromatin organization of enhancer-promoter interactions.

The complex and dynamic three-dimensional organization of chromatin within the nucleus makes understanding the control of gene expression challenging, but also opens up possible ways to epigenetically modulate gene expression. Because plants are sessile, they evolved sophisticated ways to rapidly modulate gene expression in response to environmental stress, that are thought to be coordinated by changes in chromatin conformation to mediate specific cellular and physiological responses. However, to what extent and how stress induces dynamic changes in chromatin reorganization remains poorly understood. Here, we comprehensively investigated genome-wide chromatin changes associated with transcriptional reprogramming response to heat stress in tomato. Our data show that heat stress induces rapid changes in chromatin architecture, leading to the transient formation of promoter-enhancer contacts, likely driving the expression of heat-stress responsive genes. Furthermore, we demonstrate that chromatin spatial reorganization requires HSFA1a, a transcription factor (TF) essential for heat stress tolerance in tomato. In light of our findings, we propose that TFs play a key role in controlling dynamic transcriptional responses through 3D reconfiguration of promoter-enhancer contacts.

Beyond transcription: compelling open questions in plant RNA biology.

The study of RNAs has become one of the most influential research fields in contemporary biology and biomedicine. In the last few years, new sequencing technologies have produced an explosion of new and exciting discoveries in the field but have also given rise to many open questions. Defining these questions, together with old, long-standing gaps in our knowledge, is the spirit of this article. The breadth of topics within RNA biology research is vast, and every aspect of the biology of these molecules contains countless exciting open questions. Here, we asked 12 groups to discuss their most compelling question among some plant RNA biology topics. The following vignettes cover RNA alternative splicing; RNA dynamics; RNA translation; RNA structures; R-loops; epitranscriptomics; long non-coding RNAs; small RNA production and their functions in crops; small RNAs during gametogenesis and in cross-kingdom RNA interference; and RNA-directed DNA methylation. In each section, we will present the current state-of-the-art in plant RNA biology research before asking the questions that will surely motivate future discoveries in the field. We hope this article will spark a debate about the future perspective on RNA biology and provoke novel reflections in the reader.

A conserved HSF:miR169:NF-YA loop involved in tomato and Arabidopsis heat stress tolerance.

Heat stress transcription factors (HSFs) and microRNAs (miRNAs) regulate different stress and developmental networks in plants. Regulatory feedback mechanisms are at the basis of these networks. Here, we report that plants improve their heat stress tolerance through HSF-mediated transcriptional regulation of MIR169 and post-transcriptional regulation of Nuclear Factor-YA (NF-YA) transcription factors. We show that HSFs recognize tomato (Solanum lycopersicum) and Arabidopsis MIR169 promoters using yeast one-hybrid/chromatin immunoprecipitation-quantitative PCR. Silencing tomato HSFs using virus-induced gene silencing (VIGS) reduced Sly-MIR169 levels and enhanced Sly-NF-YA9/A10 target expression. Further, Sly-NF-YA9/A10 VIGS knockdown tomato plants and Arabidopsis plants overexpressing At-MIR169d or At-nf-ya2 mutants showed a link with increased heat tolerance. In contrast, Arabidopsis plants overexpressing At-NF-YA2 and those expressing a non-cleavable At-NF-YA2 form (miR169d-resistant At-NF-YA2) as well as plants in which At-miR169d regulation is inhibited (miR169d mimic plants) were more sensitive to heat stress, highlighting NF-YA as a negative regulator of heat tolerance. Furthermore, post-transcriptional cleavage of NF-YA by elevated miR169 levels resulted in alleviation of the repression of the heat stress effector HSFA7 in tomato and Arabidopsis, revealing a retroactive control of HSFs by the miR169:NF-YA node. Hence, a regulatory feedback loop involving HSFs, miR169s and NF-YAs plays a critical role in the regulation of the heat stress response in tomato and Arabidopsis plants.

The Arabidopsis APOLO and human UPAT sequence-unrelated long noncoding RNAs can modulate DNA and histone methylation machineries in plants.

BACKGROUND: RNA-DNA hybrid (R-loop)-associated long noncoding RNAs (lncRNAs), including the Arabidopsis lncRNA AUXIN-REGULATED PROMOTER LOOP (APOLO), are emerging as important regulators of three-dimensional chromatin conformation and gene transcriptional activity. RESULTS: Here, we show that in addition to the PRC1-component LIKE HETEROCHROMATIN PROTEIN 1 (LHP1), APOLO interacts with the methylcytosine-binding protein VARIANT IN METHYLATION 1 (VIM1), a conserved homolog of the mammalian DNA methylation regulator UBIQUITIN-LIKE CONTAINING PHD AND RING FINGER DOMAINS 1 (UHRF1). The APOLO-VIM1-LHP1 complex directly regulates the transcription of the auxin biosynthesis gene YUCCA2 by dynamically determining DNA methylation and H3K27me3 deposition over its promoter during the plant thermomorphogenic response. Strikingly, we demonstrate that the lncRNA UHRF1 Protein Associated Transcript (UPAT), a direct interactor of UHRF1 in humans, can be recognized by VIM1 and LHP1 in plant cells, despite the lack of sequence homology between UPAT and APOLO. In addition, we show that increased levels of APOLO or UPAT hamper VIM1 and LHP1 binding to YUCCA2 promoter and globally alter the Arabidopsis transcriptome in a similar manner. CONCLUSIONS: Collectively, our results uncover a new mechanism in which a plant lncRNA coordinates Polycomb action and DNA methylation through the interaction with VIM1, and indicates that evolutionary unrelated lncRNAs with potentially conserved structures may exert similar functions by interacting with homolog partners.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis.

BACKGROUND: Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis. RESULTS: We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage. CONCLUSIONS: AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

The Rice Serine/Arginine Splicing Factor RS33 Regulates Pre-mRNA Splicing during Abiotic Stress Responses.

Abiotic stresses profoundly affect plant growth and development and limit crop productivity. Pre-mRNA splicing is a major form of gene regulation that helps plants cope with various stresses. Serine/arginine (SR)-rich splicing factors play a key role in pre-mRNA splicing to regulate different biological processes under stress conditions. Alternative splicing (AS) of SR transcripts and other transcripts of stress-responsive genes generates multiple splice isoforms that contribute to protein diversity, modulate gene expression, and affect plant stress tolerance. Here, we investigated the function of the plant-specific SR protein RS33 in regulating pre-mRNA splicing and abiotic stress responses in rice. The loss-of-function mutant rs33 showed increased sensitivity to salt and low-temperature stresses. Genome-wide analyses of gene expression and splicing in wild-type and rs33 seedlings subjected to these stresses identified multiple splice isoforms of stress-responsive genes whose AS are regulated by RS33. The number of RS33-regulated genes was much higher under low-temperature stress than under salt stress. Our results suggest that the plant-specific splicing factor RS33 plays a crucial role during plant responses to abiotic stresses.

Silencing the conserved small nuclear ribonucleoprotein SmD1 target gene alters susceptibility to root-knot nematodes in plants.

Root-knot nematodes (RKNs) are among the most damaging pests of agricultural crops. Meloidogyne is an extremely polyphagous genus of nematodes that can infect thousands of plant species. A few genes for resistance (R-genes) to RKN suitable for use in crop breeding have been identified, but virulent strains and species of RKN have emerged that render these R-genes ineffective. Secretion of RKN effectors targeting plant functions mediates the reprogramming of root cells into specialized feeding cells, the giant cells, essential for RKN development and reproduction. Conserved targets among plant species define the more relevant strategies for controlling nematode infection. The EFFECTOR18 (EFF18) protein from M. incognita interacts with the spliceosomal small nuclear ribonucleoprotein D1 (SmD1) in Arabidopsis (Arabidopsis thaliana), disrupting its function in alternative splicing regulation and modulating the giant cell transcriptome. We show here that EFF18 is a conserved RKN-specific effector that targets this conserved spliceosomal SmD1 protein in Solanaceae. This interaction modulates alternative splicing events produced by tomato (Solanum lycopersicum) in response to M. incognita infection. The alteration of SmD1 expression by virus-induced gene silencing in Solanaceae affects giant cell formation and nematode development. Thus, our work defines a promising conserved SmD1 target gene to develop broad resistance for the control of Meloidogyne spp. in plants.

Non-B DNA structures emerging from plant genomes.

Noncanonical (non-B) DNA structures, including G-quadruplexes, R-loops, and various DNA alternate conformations, affect genome organization, stability, and activity. Very few non-B structures have been mapped across plant genomes. Recently, Ma et al. identified intercalated-motifs (cytosine-rich, four-stranded quadruplex DNA) in rice, expanding how non-B DNA may contribute to plant gene regulation.

The lncRNA MARS modulates the epigenetic reprogramming of the marneral cluster in response to ABA.

Clustered organization of biosynthetic non-homologous genes is emerging as a characteristic feature of plant genomes. The co-regulation of clustered genes seems to largely depend on epigenetic reprogramming and three-dimensional chromatin conformation. In this study, we identified the long non-coding RNA (lncRNA) MARneral Silencing (MARS), localized inside the Arabidopsis marneral cluster, which controls the local epigenetic activation of its surrounding region in response to abscisic acid (ABA). MARS modulates the POLYCOMB REPRESSIVE COMPLEX 1 (PRC1) component LIKE HETEROCHROMATIN PROTEIN 1 (LHP1) binding throughout the cluster in a dose-dependent manner, determining H3K27me3 deposition and chromatin condensation. In response to ABA, MARS decoys LHP1 away from the cluster and promotes the formation of a chromatin loop bringing together the MARNERAL SYNTHASE 1 (MRN1) locus and a distal ABA-responsive enhancer. The enrichment of co-regulated lncRNAs in clustered metabolic genes in Arabidopsis suggests that the acquisition of novel non-coding transcriptional units may constitute an additional regulatory layer driving the evolution of biosynthetic pathways.