R-Loop Mediated trans Action of the APOLO Long Noncoding RNA.

In eukaryotes, three-dimensional genome organization is critical for transcriptional regulation of gene expression. Long noncoding RNAs (lncRNAs) can modulate chromatin conformation of spatially related genomic locations within the nucleus. Here, we show that the lncRNA APOLO (AUXIN-REGULATED PROMOTER LOOP) recognizes multiple distant independent loci in the Arabidopsis thaliana genome. We found that APOLO targets are not spatially associated in the nucleus and that APOLO recognizes its targets by short sequence complementarity and the formation of DNA-RNA duplexes (R-loops). The invasion of APOLO to the target DNA decoys the plant Polycomb Repressive Complex 1 component LHP1, modulating local chromatin 3D conformation. APOLO lncRNA coordinates the expression of distal unrelated auxin-responsive genes during lateral root development in Arabidopsis. Hence, R-loop formation and chromatin protein decoy mediate trans action of lncRNAs on distant loci. VIDEO ABSTRACT.

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.

From environmental responses to adaptation: the roles of plant lncRNAs.

As sessile organisms, plants are continuously exposed to heterogeneous and changing environments and constantly need to adapt their growth strategies. They have evolved complex mechanisms to recognize various stress factors, activate appropriate signaling pathways, and respond accordingly by reprogramming the expression of multiple genes at the transcriptional, post-transcriptional, and even epigenome levels to tolerate stressful conditions such as drought, high temperature, nutrient deficiency, and pathogenic interactions. Apart from protein-coding genes, long non-coding RNAs (lncRNAs) have emerged as key players in plant adaptation to environmental stresses. They are transcripts larger than 200 nucleotides without protein-coding potential. Still, they appear to regulate a wide range of processes, including epigenetic modifications and chromatin reorganization, as well as transcriptional and post-transcriptional modulation of gene expression, allowing plant adaptation to various environmental stresses. LncRNAs can positively or negatively modulate stress responses, affecting processes such as hormone signaling, temperature tolerance, and nutrient deficiency adaptation. Moreover, they also seem to play a role in stress memory, wherein prior exposure to mild stress enhances plant ability to adapt to subsequent stressful conditions. In this review, we summarize the contribution of lncRNAs in plant adaptation to biotic and abiotic stresses, as well as stress memory. The complex evolutionary conservation of lncRNAs is also discussed and provides insights into future research directions in this field.

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.

The Compact Root Architecture 2 systemic pathway is required for the induction of cytokinins and of the miR399 in Medicago truncatula N-satisfied plants.

Plants use a combination of sophisticated local and systemic pathways to optimize growth depending on heterogeneous nutrient availability in the soil. Legume plants can acquire mineral nitrogen (N) either through their roots or via a symbiotic interaction with N-fixing rhizobia bacteria housed in so-called root nodules. To identify shoot-to-root systemic signals acting in Medicago truncatula plants at N-deficit or N-satiety, plants were grown in a split-root experimental design, in which either high or low N was provided to a half of the root system, allowing the analysis of systemic pathways independently of any local N response. Among the plant hormone families analyzed, the cytokinin trans-Zeatin accumulated in plants at N-satiety. Cytokinin application by petiole feeding led to an inhibition of both root growth and nodulation. In addition, an exhaustive analysis of miRNAs revealed that miR2111 accumulates systemically under N-deficit in both shoots and non-treated distant roots, whereas a miRNA related to inorganic Phosphate (Pi)-acquisition, the miR399, does so in plants grown at N-satiety. These two accumulation patterns are dependent on CRA2 (Compact Root Architecture 2), a receptor required for CEP (C-terminally Encoded Peptide) signaling. Constitutive ectopic expression of the miR399 reduced nodule numbers and root biomass depending on Pi availability, suggesting that the miR399-dependent Pi-acquisition regulatory module controlled by N-availability affects the development of the whole legume plant root system.

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.

Role of MPK4 in pathogen-associated molecular pattern-triggered alternative splicing in Arabidopsis.

Alternative splicing (AS) of pre-mRNAs in plants is an important mechanism of gene regulation in environmental stress tolerance but plant signals involved are essentially unknown. Pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) is mediated by mitogen-activated protein kinases and the majority of PTI defense genes are regulated by MPK3, MPK4 and MPK6. These responses have been mainly analyzed at the transcriptional level, however many splicing factors are direct targets of MAPKs. Here, we studied alternative splicing induced by the PAMP flagellin in Arabidopsis. We identified 506 PAMP-induced differentially alternatively spliced (DAS) genes. Importantly, of the 506 PAMP-induced DAS genes, only 89 overlap with the set of 1950 PAMP-induced differentially expressed genes (DEG), indicating that transcriptome analysis does not identify most DAS events. Global DAS analysis of mpk3, mpk4, and mpk6 mutants in the absence of PAMP treatment showed no major splicing changes. However, in contrast to MPK3 and MPK6, MPK4 was found to be a key regulator of PAMP-induced DAS events as the AS of a number of splicing factors and immunity-related protein kinases is affected, such as the calcium-dependent protein kinase CPK28, the cysteine-rich receptor like kinases CRK13 and CRK29 or the FLS2 co-receptor SERK4/BKK1. Although MPK4 is guarded by SUMM2 and consequently, the mpk4 dwarf and DEG phenotypes are suppressed in mpk4 summ2 mutants, MPK4-dependent DAS is not suppressed by SUMM2, supporting the notion that PAMP-triggered MPK4 activation mediates regulation of alternative splicing.

Regulatory long non-coding RNAs in root growth and development.

As sessile organisms, plants have evolved sophisticated mechanisms of gene regulation to cope with changing environments. Among them, long non-coding RNAs (lncRNAs) are a class of RNAs regulating gene expression at both transcriptional and post-transcriptional levels. They are highly responsive to environmental cues or developmental processes and are generally involved in fine-tuning plant responses to these signals. Roots, in addition to anchoring the plant to the soil, allow it to absorb the major part of its mineral nutrients and water. Furthermore, roots directly sense environmental constraints such as mineral nutrient availability and abiotic or biotic stresses and dynamically adapt their growth and architecture. Here, we review the role of lncRNAs in the control of root growth and development. In particular, we highlight their action in fine-tuning primary root growth and the development of root lateral organs, such as lateral roots and symbiotic nodules. Lastly, we report their involvement in plant response to stresses and the regulation of nutrient assimilation and homeostasis, two processes leading to the modification of root architecture. LncRNAs could become interesting targets in plant breeding programs to subtly acclimate crops to coming environmental changes.

The Arabidopsis lncRNA ASCO modulates the transcriptome through interaction with splicing factors.

Alternative splicing (AS) is a major source of transcriptome diversity. Long noncoding RNAs (lncRNAs) have emerged as regulators of AS through different molecular mechanisms. In Arabidopsis thaliana, the AS regulators NSRs interact with the ALTERNATIVE SPLICING COMPETITOR (ASCO) lncRNA. Here, we analyze the effect of the knock-down and overexpression of ASCO at the genome-wide level and find a large number of deregulated and differentially spliced genes related to flagellin responses and biotic stress. In agreement, ASCO-silenced plants are more sensitive to flagellin. However, only a minor subset of deregulated genes overlaps with the AS defects of the nsra/b double mutant, suggesting an alternative way of action for ASCO. Using biotin-labeled oligonucleotides for RNA-mediated ribonucleoprotein purification, we show that ASCO binds to the highly conserved spliceosome component PRP8a. ASCO overaccumulation impairs the recognition of specific flagellin-related transcripts by PRP8a. We further show that ASCO also binds to another spliceosome component, SmD1b, indicating that it interacts with multiple splicing factors. Hence, lncRNAs may integrate a dynamic network including spliceosome core proteins, to modulate transcriptome reprogramming in eukaryotes.

GCN5 modulates salicylic acid homeostasis by regulating H3K14ac levels at the 5' and 3' ends of its target genes.

The modification of histones by acetyl groups has a key role in the regulation of chromatin structure and transcription. The Arabidopsis thaliana histone acetyltransferase GCN5 regulates histone modifications as part of the Spt-Ada-Gcn5 Acetyltransferase (SAGA) transcriptional coactivator complex. GCN5 was previously shown to acetylate lysine 14 of histone 3 (H3K14ac) in the promoter regions of its target genes even though GCN5 binding did not systematically correlate with gene activation. Here, we explored the mechanism through which GCN5 controls transcription. First, we fine-mapped its GCN5 binding sites genome-wide and then used several global methodologies (ATAC-seq, ChIP-seq and RNA-seq) to assess the effect of GCN5 loss-of-function on the expression and epigenetic regulation of its target genes. These analyses provided evidence that GCN5 has a dual role in the regulation of H3K14ac levels in their 5' and 3' ends of its target genes. While the gcn5 mutation led to a genome-wide decrease of H3K14ac in the 5' end of the GCN5 down-regulated targets, it also led to an increase of H3K14ac in the 3' ends of GCN5 up-regulated targets. Furthermore, genome-wide changes in H3K14ac levels in the gcn5 mutant correlated with changes in H3K9ac at both 5' and 3' ends, providing evidence for a molecular link between the depositions of these two histone modifications. To understand the biological relevance of these regulations, we showed that GCN5 participates in the responses to biotic stress by repressing salicylic acid (SA) accumulation and SA-mediated immunity, highlighting the role of this protein in the regulation of the crosstalk between diverse developmental and stress-responsive physiological programs. Hence, our results demonstrate that GCN5, through the modulation of H3K14ac levels on its targets, controls the balance between biotic and abiotic stress responses and is a master regulator of plant-environmental interactions.

Landscape of the Noncoding Transcriptome Response of Two Arabidopsis Ecotypes to Phosphate Starvation.

Root architecture varies widely between species; it even varies between ecotypes of the same species, despite strong conservation of the coding portion of their genomes. By contrast, noncoding RNAs evolve rapidly between ecotypes and may control their differential responses to the environment, since several long noncoding RNAs (lncRNAs) are known to quantitatively regulate gene expression. Roots from ecotypes Columbia and Landsberg erecta of Arabidopsis (Arabidopsis thaliana) respond differently to phosphate starvation. Here, we compared transcriptomes (mRNAs, lncRNAs, and small RNAs) of root tips from these two ecotypes during early phosphate starvation. We identified thousands of lncRNAs that were largely conserved at the DNA level in these ecotypes. In contrast to coding genes, many lncRNAs were specifically transcribed in one ecotype and/or differentially expressed between ecotypes independent of phosphate availability. We further characterized these ecotype-related lncRNAs and studied their link with small interfering RNAs. Our analysis identified 675 lncRNAs differentially expressed between the two ecotypes, including antisense RNAs targeting key regulators of root-growth responses. Misregulation of several lincRNAs showed that at least two ecotype-related lncRNAs regulate primary root growth in ecotype Columbia. RNA-sequencing analysis following deregulation of lncRNA NPC48 revealed a potential link with root growth and transport functions. This exploration of the noncoding transcriptome identified ecotype-specific lncRNA-mediated regulation in root apexes. The noncoding genome may harbor further mechanisms involved in ecotype adaptation of roots to different soil environments.

Full Length Transcriptome Highlights the Coordination of Plastid Transcript Processing.

Plastid gene expression involves many post-transcriptional maturation steps resulting in a complex transcriptome composed of multiple isoforms. Although short-read RNA-Seq has considerably improved our understanding of the molecular mechanisms controlling these processes, it is unable to sequence full-length transcripts. This information is crucial, however, when it comes to understanding the interplay between the various steps of plastid gene expression. Here, we describe a protocol to study the plastid transcriptome using nanopore sequencing. In the leaf of Arabidopsis thaliana, with about 1.5 million strand-specific reads mapped to the chloroplast genome, we could recapitulate most of the complexity of the plastid transcriptome (polygenic transcripts, multiple isoforms associated with post-transcriptional processing) using virtual Northern blots. Even if the transcripts longer than about 2500 nucleotides were missing, the study of the co-occurrence of editing and splicing events identified 42 pairs of events that were not occurring independently. This study also highlighted a preferential chronology of maturation events with splicing happening after most sites were edited.

Arabidopsis HEAT SHOCK TRANSCRIPTION FACTORA1b regulates multiple developmental genes under benign and stress conditions.

In Arabidopsis thaliana, HEAT SHOCK TRANSCRIPTION FACTORA1b (HSFA1b) controls resistance to environmental stress and is a determinant of reproductive fitness by influencing seed yield. To understand how HSFA1b achieves this, we surveyed its genome-wide targets (ChIP-seq) and its impact on the transcriptome (RNA-seq) under non-stress (NS), heat stress (HS) in the wild type, and in HSFA1b-overexpressing plants under NS. A total of 952 differentially expressed HSFA1b-targeted genes were identified, of which at least 85 are development associated and were bound predominantly under NS. A further 1780 genes were differentially expressed but not bound by HSFA1b, of which 281 were classified as having development-associated functions. These genes are indirectly regulated through a hierarchical network of 27 transcription factors (TFs). Furthermore, we identified 480 natural antisense non-coding RNA (cisNAT) genes bound by HSFA1b, defining a further mode of indirect regulation. Finally, HSFA1b-targeted genomic features not only harboured heat shock elements, but also MADS box, LEAFY, and G-Box promoter motifs. This revealed that HSFA1b is one of eight TFs that target a common group of stress defence and developmental genes. We propose that HSFA1b transduces environmental cues to many stress tolerance and developmental genes to allow plants to adjust their growth and development continually in a varying environment.

GDP-L-fucose is required for boundary definition in plants.

The CUP-SHAPED COTYLEDON (CUC) transcription factors control plant boundary formation, thus allowing the emergence of novel growth axes. While the developmental roles of the CUC genes in different organs and across species are well characterized, upstream and downstream events that contribute to their function are still poorly understood. To identify new players in this network, we performed a suppressor screen of CUC2g-m4, a line overexpressing CUC2 that has highly serrated leaves. We identified a mutation that simplifies leaf shape and affects MURUS1 (MUR1), which is responsible for GDP-L-fucose production. Using detailed morphometric analysis, we show that GDP-L-fucose has an essential role in leaf shape acquisition by sustaining differential growth at the leaf margins. Accordingly, reduced CUC2 expression levels are observed in mur1 leaves. Furthermore, genetic analyses reveal a conserved role for GDP-L-fucose in different developmental contexts where it contributes to organ separation in the same pathway as CUC2. Taken together, our results reveal that GDP-L-fucose is necessary for proper establishment of boundary domains in various developmental contexts.

The iRoCS Toolbox--3D analysis of the plant root apical meristem at cellular resolution.

To achieve a detailed understanding of processes in biological systems, cellular features must be quantified in the three-dimensional (3D) context of cells and organs. We described use of the intrinsic root coordinate system (iRoCS) as a reference model for the root apical meristem of plants. iRoCS enables direct and quantitative comparison between the root tips of plant populations at single-cell resolution. The iRoCS Toolbox automatically fits standardized coordinates to raw 3D image data. It detects nuclei or segments cells, automatically fits the coordinate system, and groups the nuclei/cells into the root's tissue layers. The division status of each nucleus may also be determined. The only manual step required is to mark the quiescent centre. All intermediate outputs may be refined if necessary. The ability to learn the visual appearance of nuclei by example allows the iRoCS Toolbox to be easily adapted to various phenotypes. The iRoCS Toolbox is provided as an open-source software package, licensed under the GNU General Public License, to make it accessible to a broad community. To demonstrate the power of the technique, we measured subtle changes in cell division patterns caused by modified auxin flux within the Arabidopsis thaliana root apical meristem.

Modification of plant Rac/Rop GTPase signalling using bacterial toxin transgenes.

Bacterial protein toxins which modify Rho GTPase are useful for the analysis of Rho signalling in animal cells, but these toxins cannot be taken up by plant cells. We demonstrate in vitro deamidation of Arabidopsis Rop4 by Escherichia coli Cytotoxic Necrotizing Factor 1 (CNF1) and glucosylation by Clostridium difficile toxin B. Expression of the catalytic domain of CNF1 caused modification and activation of co-expressed Arabidopsis Rop4 GTPase in tobacco leaves, resulting in hypersensitive-like cell death. By contrast, the catalytic domain of toxin B modified and inactivated co-expressed constitutively active Rop4, blocking the hypersensitive response caused by over-expression of active Rops. In transgenic Arabidopsis, both CNF1 and toxin B inhibited Rop-dependent polar morphogenesis of leaf epidermal cells. Toxin B expression also inhibited Rop-dependent morphogenesis of root hairs and trichome branching, and resulted in root meristem enlargement and dwarf growth. Our results show that CNF1 and toxin B transgenes are effective tools in Rop GTPase signalling studies.

A miR169 isoform regulates specific NF-YA targets and root architecture in Arabidopsis.

In plants, roots are essential for water and nutrient acquisition. MicroRNAs (miRNAs) regulate their target mRNAs by transcript cleavage and/or inhibition of protein translation and are known as major post-transcriptional regulators of various developmental pathways and stress responses. In Arabidopsis thaliana, four isoforms of miR169 are encoded by 14 different genes and target diverse mRNAs, encoding subunits A of the NF-Y transcription factor complex. These miRNA isoforms and their targets have previously been linked to nutrient signalling in plants. By using mimicry constructs against different isoforms of miR169 and miR-resistant versions of NF-YA genes we analysed the role of specific miR169 isoforms in root growth and branching. We identified a regulatory node involving the particular miR169defg isoform and NF-YA2 and NF-YA10 genes that acts in the control of primary root growth. The specific expression of MIM169defg constructs altered specific cell type numbers and dimensions in the root meristem. Preventing miR169defg-regulation of NF-YA2 indirectly affected laterial root initiation. We also showed that the miR169defg isoform affects NF-YA2 transcripts both at mRNA stability and translation levels. We propose that a specific miR169 isoform and the NF-YA2 target control root architecture in Arabidopsis.

Evolution and diverse roles of the CUP-SHAPED COTYLEDON genes in Arabidopsis leaf development.

CUP-SHAPED COTYLEDON2 (CUC2) and the interacting microRNA miR164 regulate leaf margin dissection. Here, we further investigate the evolution and the specific roles of the CUC1 to CUC3 genes during Arabidopsis thaliana leaf serration. We show that CUC2 is essential for dissecting the leaves of a wide range of lobed/serrated Arabidopsis lines. Inactivation of CUC3 leads to a partial suppression of the serrations, indicating a role for this gene in leaf shaping. Morphometric analysis of leaf development and genetic analysis provide evidence for different temporal contributions of CUC2 and CUC3. Chimeric constructs mixing CUC regulatory sequences with different coding sequences reveal both redundant and specific roles for the three CUC genes that could be traced back to changes in their expression pattern or protein activity. In particular, we show that CUC1 triggers the formation of leaflets when ectopically expressed instead of CUC2 in the developing leaves. These divergent fates of the CUC1 and CUC2 genes after their formation by the duplication of a common ancestor is consistent with the signature of positive selection detected on the ancestral branch to CUC1. Combining experimental observations with the retraced origin of the CUC genes in the Brassicales, we propose an evolutionary scenario for the CUC genes.

Variational attenuation correction in two-view confocal microscopy.

BACKGROUND: Absorption and refraction induced signal attenuation can seriously hinder the extraction of quantitative information from confocal microscopic data. This signal attenuation can be estimated and corrected by algorithms that use physical image formation models. Especially in thick heterogeneous samples, current single view based models are unable to solve the underdetermined problem of estimating the attenuation-free intensities. RESULTS: We present a variational approach to estimate both, the real intensities and the spatially variant attenuation from two views of the same sample from opposite sides. Assuming noise-free measurements throughout the whole volume and pure absorption, this would in theory allow a perfect reconstruction without further assumptions. To cope with real world data, our approach respects photon noise, estimates apparent bleaching between the two recordings, and constrains the attenuation field to be smooth and sparse to avoid spurious attenuation estimates in regions lacking valid measurements. CONCLUSIONS: We quantify the reconstruction quality on simulated data and compare it to the state-of-the art two-view approach and commonly used one-factor-per-slice approaches like the exponential decay model. Additionally we show its real-world applicability on model organisms from zoology (zebrafish) and botany (Arabidopsis). The results from these experiments show that the proposed approach improves the quantification of confocal microscopic data of thick specimen.

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.