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.

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.

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.

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.

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.

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.

ChronoRoot: High-throughput phenotyping by deep segmentation networks reveals novel temporal parameters of plant root system architecture.

BACKGROUND: Deep learning methods have outperformed previous techniques in most computer vision tasks, including image-based plant phenotyping. However, massive data collection of root traits and the development of associated artificial intelligence approaches have been hampered by the inaccessibility of the rhizosphere. Here we present ChronoRoot, a system that combines 3D-printed open-hardware with deep segmentation networks for high temporal resolution phenotyping of plant roots in agarized medium. RESULTS: We developed a novel deep learning-based root extraction method that leverages the latest advances in convolutional neural networks for image segmentation and incorporates temporal consistency into the root system architecture reconstruction process. Automatic extraction of phenotypic parameters from sequences of images allowed a comprehensive characterization of the root system growth dynamics. Furthermore, novel time-associated parameters emerged from the analysis of spectral features derived from temporal signals. CONCLUSIONS: Our work shows that the combination of machine intelligence methods and a 3D-printed device expands the possibilities of root high-throughput phenotyping for genetics and natural variation studies, as well as the screening of clock-related mutants, revealing novel root traits.

Insights into long non-coding RNA regulation of anthocyanin carrot root pigmentation.

Carrot (Daucus carota L.) is one of the most cultivated vegetable in the world and of great importance in the human diet. Its storage organs can accumulate large quantities of anthocyanins, metabolites that confer the purple pigmentation to carrot tissues and whose biosynthesis is well characterized. Long non-coding RNAs (lncRNAs) play critical roles in regulating gene expression of various biological processes in plants. In this study, we used a high throughput stranded RNA-seq to identify and analyze the expression profiles of lncRNAs in phloem and xylem root samples using two genotypes with a strong difference in anthocyanin production. We discovered and annotated 8484 new genes, including 2095 new protein-coding and 6373 non-coding transcripts. Moreover, we identified 639 differentially expressed lncRNAs between the phenotypically contrasted genotypes, including certain only detected in a particular tissue. We then established correlations between lncRNAs and anthocyanin biosynthesis genes in order to identify a molecular framework for the differential expression of the pathway between genotypes. A specific natural antisense transcript linked to the DcMYB7 key anthocyanin biosynthetic transcription factor suggested how the regulation of this pathway may have evolved between genotypes.

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.

Wheat chromatin architecture is organized in genome territories and transcription factories.

BACKGROUND: Polyploidy is ubiquitous in eukaryotic plant and fungal lineages, and it leads to the co-existence of several copies of similar or related genomes in one nucleus. In plants, polyploidy is considered a major factor in successful domestication. However, polyploidy challenges chromosome folding architecture in the nucleus to establish functional structures. RESULTS: We examine the hexaploid wheat nuclear architecture by integrating RNA-seq, ChIP-seq, ATAC-seq, Hi-C, and Hi-ChIP data. Our results highlight the presence of three levels of large-scale spatial organization: the arrangement into genome territories, the diametrical separation between facultative and constitutive heterochromatin, and the organization of RNA polymerase II around transcription factories. We demonstrate the micro-compartmentalization of transcriptionally active genes determined by physical interactions between genes with specific euchromatic histone modifications. Both intra- and interchromosomal RNA polymerase-associated contacts involve multiple genes displaying similar expression levels. CONCLUSIONS: Our results provide new insights into the physical chromosome organization of a polyploid genome, as well as on the relationship between epigenetic marks and chromosome conformation to determine a 3D spatial organization of gene expression, a key factor governing gene transcription in polyploids.

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.

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.

High-quality genome sequence of white lupin provides insight into soil exploration and seed quality.

White lupin (Lupinus albus L.) is an annual crop cultivated for its protein-rich seeds. It is adapted to poor soils due to the production of cluster roots, which are made of dozens of determinate lateral roots that drastically improve soil exploration and nutrient acquisition (mostly phosphate). Using long-read sequencing technologies, we provide a high-quality genome sequence of a cultivated accession of white lupin (2n = 50, 451 Mb), as well as de novo assemblies of a landrace and a wild relative. We describe a modern accession displaying increased soil exploration capacity through early establishment of lateral and cluster roots. We also show how seed quality may have been impacted by domestication in term of protein profiles and alkaloid content. The availability of a high-quality genome assembly together with companion genomic and transcriptomic resources will enable the development of modern breeding strategies to increase and stabilize white lupin yield.

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.

Whole-genome landscape of Medicago truncatula symbiotic genes.

Advances in deciphering the functional architecture of eukaryotic genomes have been facilitated by recent breakthroughs in sequencing technologies, enabling a more comprehensive representation of genes and repeat elements in genome sequence assemblies, as well as more sensitive and tissue-specific analyses of gene expression. Here we show that PacBio sequencing has led to a substantially improved genome assembly of Medicago truncatula A17, a legume model species notable for endosymbiosis studies1, and has enabled the identification of genome rearrangements between genotypes at a near-base-pair resolution. Annotation of the new M. truncatula genome sequence has allowed for a thorough analysis of transposable elements and their dynamics, as well as the identification of new players involved in symbiotic nodule development, in particular 1,037 upregulated long non-coding RNAs (lncRNAs). We have also discovered that a substantial proportion (~35% and 38%, respectively) of the genes upregulated in nodules or expressed in the nodule differentiation zone colocalize in genomic clusters (270 and 211, respectively), here termed symbiotic islands. These islands contain numerous expressed lncRNA genes and display differentially both DNA methylation and histone marks. Epigenetic regulations and lncRNAs are therefore attractive candidate elements for the orchestration of symbiotic gene expression in the M. truncatula genome.

Nuclear Speckle RNA Binding Proteins Remodel Alternative Splicing and the Non-coding Arabidopsis Transcriptome to Regulate a Cross-Talk Between Auxin and Immune Responses.

Nuclear speckle RNA binding proteins (NSRs) act as regulators of alternative splicing (AS) and auxin-regulated developmental processes such as lateral root formation in Arabidopsis thaliana. These proteins were shown to interact with specific alternatively spliced mRNA targets and at least with one structured lncRNA, named Alternative Splicing Competitor RNA. Here, we used genome-wide analysis of RNAseq to monitor the NSR global role on multiple tiers of gene expression, including RNA processing and AS. NSRs affect AS of 100s of genes as well as the abundance of lncRNAs particularly in response to auxin. Among them, the FPA floral regulator displayed alternative polyadenylation and differential expression of antisense COOLAIR lncRNAs in nsra/b mutants. This may explains the early flowering phenotype observed in nsra and nsra/b mutants. GO enrichment analysis of affected lines revealed a novel link of NSRs with the immune response pathway. A RIP-seq approach on an NSRa fusion protein in mutant background identified that lncRNAs are privileged direct targets of NSRs in addition to specific AS mRNAs. The interplay of lncRNAs and AS mRNAs in NSR-containing complexes may control the crosstalk between auxin and the immune response pathway.

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.