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.

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 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.

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.

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.

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.

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.

Polycomb-dependent differential chromatin compartmentalization determines gene coregulation in Arabidopsis

In animals, distant H3K27me3-marked Polycomb targets can establish physical interactions forming repressive chromatin hubs. In plants, growing evidence suggests that H3K27me3 acts directly or indirectly to regulate chromatin interactions, although how this histone modification modulates 3D chromatin architecture remains elusive. To decipher the impact of the dynamic deposition of H3K27me3 on the Arabidopsis thaliana nuclear interactome, we combined genetics, transcriptomics, and several 3D epigenomic approaches. By analyzing mutants defective for histone H3K27 methylation or demethylation, we uncovered the crucial role of this chromatin mark in short- and previously unnoticed long-range chromatin loop formation. We found that a reduction in H3K27me3 levels led to a decrease in the interactions within Polycomb-associated repressive domains. Regions with lower H3K27me3 levels in the H3K27 methyltransferase clf mutant established new interactions with regions marked with H3K9ac, a histone modification associated with active transcription, indicating that a reduction in H3K27me3 levels induces a global reconfiguration of chromatin architecture. Altogether, our results reveal that the 3D genome organization is tightly linked to reversible histone modifications that govern chromatin interactions. Consequently, nuclear organization dynamics shapes the transcriptional reprogramming during plant development and places H3K27me3 as a key feature in the coregulation of distant genes.

The lncRNA APOLO interacts with the transcription factor WRKY42 to trigger root hair cell expansion in response to cold

Plant long noncoding RNAs (lncRNAs) have emerged as important regulators of chromatin dynamics, impacting on transcriptional programs leading to different developmental outputs. The lncRNA AUXIN-REGULATED PROMOTER LOOP (APOLO) directly recognizes multiple independent loci across the Arabidopsis genome and modulates their three-dimensional chromatin conformation, leading to transcriptional shifts. Here, we show that APOLO recognizes the locus encoding the root hair (RH) master regulator ROOT HAIR DEFECTIVE 6 (RHD6) and controls RHD6 transcriptional activity, leading to cold-enhanced RH elongation through the consequent activation of the transcription factor gene RHD6-like RSL4. Furthermore, we demonstrate that APOLO interacts with the transcription factor WRKY42 and modulates its binding to the RHD6 promoter. WRKY42 is required for the activation of RHD6 by low temperatures and WRKY42 deregulation impairs cold-induced RH expansion. Collectively, our results indicate that a novel ribonucleoprotein complex with APOLO and WRKY42 forms a regulatory hub to activate RHD6 by shaping its epigenetic environment and integrate signals governing RH growth and development.

Overlapping roles of spliceosomal components SF3B1 and PHF5A in rice splicing regulation

The SF3B complex, a multiprotein component of the U2 snRNP of the spliceosome, plays a crucial role in recognizing branch point sequence and facilitates spliceosome assembly and activation. Several chemicals that bind SF3B1 and PHF5A subunits of the SF3B complex inhibit splicing. We recently generated a splicing inhibitor-resistant SF3B1 mutant named SF3B1GEX1ARESISTANT 4 (SGR4) using CRISPR-mediated directed evolution, whereas splicing inhibitor-resistant mutant of PHF5A (Overexpression-PHF5A GEX1A Resistance, OGR) was generated by expressing an engineered version PHF5A-Y36C. Global analysis of splicing in wild type and these two mutants revealed the role of SF3B1 and PHF5A in splicing regulation. This analysis uncovered a set of genes whose intron retention is regulated by both proteins. Further analysis of these retained introns revealed that they are shorter, have a higher GC content, and contain shorter and weaker polypyrimidine tracts. Furthermore, splicing inhibition increased seedlings sensitivity to salt stress, consistent with emerging roles of splicing regulation in stress responses. In summary, we uncovered the functions of two members of the plant branch point recognition complex. The novel strategies described here should be broadly applicable in elucidating functions of splicing regulators, especially in studying the functions of redundant paralogs in plants.

Butt et al. used CRISPR-mediated directed evolution to generate rice mutants for the spliceosome components SF3B1 and PHF5A. They demonstrate that these mutants have different levels of sensitivity to salt treatments and suggest that the strategies they employed can be used in the future to study functions of redundant paralogs in plants.

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.

Plant Long Noncoding RNAs: New Players in the Field of Post-Transcriptional Regulations

The first reference to the "C-value paradox" reported an apparent imbalance between organismal genome size and morphological complexity. Since then, next-generation sequencing has revolutionized genomic research and revealed that eukaryotic transcriptomes contain a large fraction of non-protein-coding components. Eukaryotic genomes are pervasively transcribed and noncoding regions give rise to a plethora of noncoding RNAs with undeniable biological functions. Among them, long noncoding RNAs (lncRNAs) seem to represent a new layer of gene expression regulation, participating in a wide range of molecular mechanisms at the transcriptional and post-transcriptional levels. In addition to their role in epigenetic regulation, plant lncRNAs have been associated with the degradation of complementary RNAs, the regulation of alternative splicing, protein sub-cellular localization, the promotion of translation and protein post-translational modifications. In this review, we report and integrate numerous and complex mechanisms through which long noncoding transcripts regulate post-transcriptional gene expression in plants.

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.

Long noncoding RNAs shape transcription in plants

The advent of novel high-throughput sequencing techniques has revealed that eukaryotic genomes are massively transcribed although only a small fraction of RNAs exhibits protein-coding capacity. In the last years, long noncoding RNAs (lncRNAs) have emerged as regulators of eukaryotic gene expression in a wide range of molecular mechanisms. Plant lncRNAs can be transcribed by alternative RNA polymerases, acting directly as long transcripts or can be processed into active small RNAs. Several lncRNAs have been recently shown to interact with chromatin, DNA or nuclear proteins to condition the epigenetic environment of target genes or modulate the activity of transcriptional complexes. In this review, we will summarize the recent discoveries about the actions of plant lncRNAs in the regulation of gene expression at the transcriptional level.

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.

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.

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.

White lupin is an annual crop cultivated for protein rich seeds and can produce cluster roots for efficient phosphate acquisition. Here, the authors generate high quality genome assemblies of a cultivated accession, a landrace, and a wild relative and provides insight into soil exploration and seed quality.

Alternative Splicing in the Regulation of Plant-Microbe Interactions

As sessile organisms, plants are continuously exposed to a wide range of biotic interactions. While some biotic interactions are beneficial or even essential for the plant (e.g. rhizobia and mycorrhiza), others such as pathogens are detrimental and require fast adaptation. Plants partially achieve this growth and developmental plasticity by modulating the repertoire of genes they express. In the past few years, high-throughput transcriptome sequencing have revealed that, in addition to transcriptional control of gene expression, post-transcriptional processes, notably alternative splicing (AS), emerged as a key mechanism for gene regulation during plant adaptation to the environment. AS not only can increase proteome diversity by generating multiple transcripts from a single gene but also can reduce gene expression by yielding isoforms degraded by mechanisms such as nonsense-mediated mRNA decay. In this review, we will summarize recent discoveries detailing the contribution of AS to the regulation of plant-microbe interactions, with an emphasis on the modulation of immunity receptor function and other components of the signaling pathways that deal with pathogen responses. We will also discuss emerging evidences that AS could contribute to dynamic reprogramming of the plant transcriptome during beneficial interactions, such as the legume-symbiotic interaction.

Assessing the Response of Small RNA Populations to Allopolyploidy Using Resynthesized Brassica napus Allotetraploids

Allopolyploidy, combining interspecific hybridization with whole genome duplication, has had significant impact on plant evolution. Its evolutionary success is related to the rapid and profound genome reorganizations that allow neoallopolyploids to form and adapt. Nevertheless, how neoallopolyploid genomes adapt to regulate their expression remains poorly understood. The hypothesis of a major role for small noncoding RNAs (sRNAs) in mediating the transcriptional response of neoallopolyploid genomes has progressively emerged. Generally, 21-nt sRNAs mediate posttranscriptional gene silencing by mRNA cleavage, whereas 24-nt sRNAs repress transcription (transcriptional gene silencing) through epigenetic modifications. Here, we characterize the global response of sRNAs to allopolyploidy in Brassica, using three independently resynthesized Brassica napus allotetraploids originating from crosses between diploid Brassica oleracea and Brassica rapa accessions, surveyed at two different generations in comparison with their diploid progenitors. Our results suggest an immediate but transient response of specific sRNA populations to allopolyploidy. These sRNA populations mainly target noncoding components of the genome but also target the transcriptional regulation of genes involved in response to stresses and in metabolism; this suggests a broad role in adapting to allopolyploidy. We finally identify the early accumulation of both 21- and 24-nt sRNAs involved in regulating the same targets, supporting a posttranscriptional gene silencing to transcriptional gene silencing shift at the first stages of the neoallopolyploid formation. We propose that reorganization of sRNA production is an early response to allopolyploidy in order to control the transcriptional reactivation of various noncoding elements and stress-related genes, thus ensuring genome stability during the first steps of neoallopolyploid formation.

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 studies¹, 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.

Thermopriming triggers splicing memory in Arabidopsis

Abiotic and biotic stresses limit crop productivity. Exposure to a non-lethal stress, referred to as priming, can allow plants to survive subsequent and otherwise lethal conditions; the priming effect persists even after a prolonged stress-free period. However, the molecular mechanisms underlying priming are not fully understood. Here, we investigated the molecular basis of heat-shock memory and the role of priming in Arabidopsis thaliana. Comprehensive analysis of transcriptome-wide changes in gene expression and alternative splicing in primed and non-primed plants revealed that alternative splicing functions as a novel component of heat-shock memory. We show that priming of plants with a non-lethal heat stress results in de-repression of splicing after a second exposure to heat stress. By contrast, non-primed plants showed significant repression of splicing. These observations link 'splicing memory' to the ability of plants to survive subsequent and otherwise lethal heat stress. This newly discovered priming-induced splicing memory may represent a general feature of heat-stress responses in plants and other organisms as many of the key components are conserved among eukaryotes. Furthermore, this finding could facilitate the development of novel approaches to improve plant survival under extreme heat stress.

Cancer Risk Related to Uncommon Migration within Italy

Aims and background

Given the industrialization in Italy over this past century much migration has occurred within the country especially from southern to northern regions. Following repeated drainings of the pre-existing marsh area (Pontina plain) during the 1930s the Latina province received an unusual north-south immigration from the regions of Veneto, Friuli and Emilia Romagna. This consisted principally of manual workers, farmers and their families. Four new towns developed after a few years (Littoria, later renamed Latina, Sabaudia, Pontinia and Aprilia), and the whole province quickly reached a population of 60,000. The availability of a population-based Cancer Registry in the Latina province allowed us to assess the cancer risk in this migrant population.

Methods

Standardized Incidence Ratios (SIRs) according to cancer site and sex were computed for residents over the age of 55 years, born in northern Italy. Population data, by sex, age and region of birth were based on the 1981 census. The age-sex-site specific incidence rates for the 1983-1987 period for the entire population of the Latina province over 55 years of age were used as standard.

Results

A significant excess of cancer risk for subjects of both sexes born in northern Italy was found. In addition, a statistically significant higher risk was observed for the cancers of the lung, skin (non-melanomas) and prostate in males, and of the mouth-pharynx, lung and skin (non-melanomas) in females.

Conclusions

SIRs for all sites confirm the findings from other studies on migrants in Italy and strongly support the hypothesis that the place of birth has an important influence on the frequency of cancer. Some possible etiological factors are suggested for cancer sites with higher frequencies in northern-born subjects.

Splicing regulation by long noncoding RNAs

Massive high-throughput sequencing techniques allowed the identification of thousands of noncoding RNAs (ncRNAs) and a plethora of different mRNA processing events occurring in higher organisms. Long ncRNAs can act directly as long transcripts or can be processed into active small si/miRNAs. They can modulate mRNA cleavage, translational repression or the epigenetic landscape of their target genes. Recently, certain long ncRNAs have been shown to play a crucial role in the regulation of alternative splicing in response to several stimuli or during disease. In this review, we focus on recent discoveries linking gene regulation by alternative splicing and its modulation by long and small ncRNAs.

Root Development in Medicago truncatula: Lessons from Genetics to Functional Genomics

This decade introduced "omics" approaches, such as genomics, transcriptomics, proteomics, and metabolomics in association with reverse and forward genetic approaches, developed earlier, to try to identify molecular pathways involved in the development or in the response to environmental conditions as well as in animals and plants. This review summarizes studies that utilized "omics" strategies to unravel the root development in the model legume Medicago truncatula and how external factors such as soil mineral status or the presence of bacteria and fungi affect root system architecture in this species. We also compare these "omics" data to the knowledges concerning the Arabidopsis thaliana root development, nowadays considered as the model of allorhiz root systems. However, unlike legumes, this species is unable to interact with soil nitrogen-fixing rhizobia and arbuscular-mycorrhizal (AM) fungi to develop novel root-derived symbiotic structures. Differences in root organization, development, and regulatory pathways between these two model species have been highlighted.

Stable Inactivation of MicroRNAs in Medicago truncatula Roots

MicroRNAs are key regulators in the development processes or stress responses in plants. In the last decade, several conserved or non-conserved microRNAs have been identified in Medicago truncatula. Different strategies leading to the inactivation of microRNAs in plants have been described. Here, we propose a protocol for an effective inactivation of microRNAs using a STTM strategy in M. truncatula transgenic roots.

One Gene, Many Proteins: Mapping Cell-Specific Alternative Splicing in Plants

Pre-mRNA alternative splicing (AS) generates protein variants from a single gene that can create novel regulatory opportunities. In this issue of Developmental Cell, Li et al. (2016) present a high-resolution expression map of AS events in Arabidopsis root tissues, giving insight into cell-type- and stage-specific AS mechanisms in plants.

The Arabidopsis SWI/SNF protein BAF60 mediates seedling growth control by modulating DNA accessibility

BACKGROUND

Plant adaptive responses to changing environments involve complex molecular interplays between intrinsic and external signals. Whilst much is known on the signaling components mediating diurnal, light, and temperature controls on plant development, their influence on chromatin-based transcriptional controls remains poorly explored.

RESULTS

In this study we show that a SWI/SNF chromatin remodeler subunit, BAF60, represses seedling growth by modulating DNA accessibility of hypocotyl cell size regulatory genes. BAF60 binds nucleosome-free regions of multiple G box-containing genes, opposing in cis the promoting effect of the photomorphogenic and thermomorphogenic regulator Phytochrome Interacting Factor 4 (PIF4) on hypocotyl elongation. Furthermore, BAF60 expression level is regulated in response to light and daily rhythms.

CONCLUSIONS

These results unveil a short path between a chromatin remodeler and a signaling component to fine-tune plant morphogenesis in response to environmental conditions.

TMV induces RNA decay pathways to modulate gene silencing and disease symptoms

RNA decay pathways comprise a combination of RNA degradation mechanisms that are implicated in gene expression, development and defense responses in eukaryotes. These mechanisms are known as the RNA Quality Control or RQC pathways. In plants, another important RNA degradation mechanism is the post-transcriptional gene silencing (PTGS) mediated by small RNAs (siRNAs). Notably, the RQC pathway antagonizes PTGS by preventing the entry of dysfunctional mRNAs into the silencing pathway to avoid global degradation of mRNA by siRNAs. Viral transcripts must evade RNA degrading mechanisms, thus viruses encode PTGS suppressor proteins to counteract viral RNA silencing. Here, we demonstrate that tobacco plants infected with TMV and transgenic lines expressing TMV MP and CP (coat protein) proteins (which are not linked to the suppression of silencing) display increased transcriptional levels of RNA decay genes. These plants also showed accumulation of cytoplasmic RNA granules with altered structure, increased rates of RNA decay for transgenes and defective transgene PTGS amplification. Furthermore, knockdown of RRP41 or RRP43 RNA exosome components led to lower levels of TMV accumulation with milder symptoms after infection, several developmental defects and miRNA deregulation. Thus, we propose that TMV proteins induce RNA decay pathways (in particular exosome components) to impair antiviral PTGS and this defensive mechanism would constitute an additional counter-defense strategy that lead to disease symptoms.

Small RNA profiles in soybean primary root tips under water deficit

BACKGROUND

Soybean (Glycine max) production is significantly hampered by frequent droughts in many regions of the world including the United States. Identifying microRNA (miRNA)-controlled posttranscriptional gene regulation under drought will enhance our understanding of molecular basis of drought tolerance in this important cash crop. Indeed, miRNA profiles in soybean exposed to drought were studied but not from the primary root tips, which is not only a main zone of water uptake but also critical for water stress sensing and signaling.

METHODS

Here we report miRNA profiles specifically from well-watered and water-stressed primary root tips (0 to 8 mm from the root apex) of soybean. Small RNA sequencing confirmed the expression of vastly diverse miRNA (303 individual miRNAs) population, and, importantly several conserved miRNAs were abundantly expressed in primary root tips.

RESULTS

Notably, 12 highly conserved miRNA families were differentially regulated in response to water-deficit; six were upregulated while six others were downregulated at least by one fold (log2) change. Differentially regulated soybean miRNAs are targeting genes include auxin response factors, Cu/Zn Superoxide dismutases, laccases and plantacyanin and several others.

CONCLUSIONS

These results highlighted the importance of miRNAs in primary root tips both under control and water-deficit conditions; under control conditions, miRNAs could be important for cell division, cell elongation and maintenance of the root apical meristem activity including quiescent centre whereas under water stress differentially regulated miRNAs could decrease auxin signaling and oxidative stress as well as other metabolic processes that save energy and water.

Detection of generic differential RNA processing events from RNA-seq data

RNA-seq data analysis has revealed abundant alternative splicing in eukaryotic mRNAs. However, splicing is only one of many processing events that transcripts may undergo during their lifetime. We present here RNAprof (RNA profile analysis), a program for the detection of differential processing events from the comparison of RNA-seq experiments. RNAprof implements a specific gene-level normalization procedure and compares RNA-seq coverage profiles at nucleotide resolution to detect regions of significant coverage differences, independently of splice sites or other gene features. We used RNAprof to analyze the effect of alternative-splicing regulators NSRa and NSRb on the Arabidopsis thaliana transcriptome. A number of intron retention events and alternative transcript structures were specifically detected by RNAprof and confirmed by qRT-PCR. Further tests using a public Mus musculus RNA-seq dataset and comparisons with other RNA isoform predictors showed that RNAprof uniquely identified sets of highly significant processing events as well as other relevant library-specific differences in RNA-seq profiles. This highlights an important layer of variation that remains undetected by current protocols for RNA-seq analysis.

LHP1 Regulates H3K27me3 Spreading and Shapes the Three-Dimensional Conformation of the Arabidopsis Genome

Precise expression patterns of genes in time and space are essential for proper development of multicellular organisms. Dynamic chromatin conformation and spatial organization of the genome constitute a major step in this regulation to modulate developmental outputs. Polycomb repressive complexes (PRCs) mediate stable or flexible gene repression in response to internal and environmental cues. In Arabidopsis thaliana, LHP1 co-localizes with H3K27me3 epigenetic marks throughout the genome and interacts with PRC1 and PRC2 members as well as with a long noncoding RNA. Here, we show that LHP1 is responsible for the spreading of H3K27me3 towards the 3' end of the gene body. We also identified a subset of LHP1-activated genes and demonstrated that LHP1 shapes local chromatin topology in order to control transcriptional co-regulation. Our work reveals a general role of LHP1 from local to higher conformation levels of chromatin configuration to determine its accessibility to define gene expression patterns.

Arabidopsis CLAVATA1 and CLAVATA2 receptors contribute to Ralstonia solanacearum pathogenicity through a miR169-dependent pathway

Bacterial wilt caused by Ralstonia solanacearum is one of the most destructive bacterial plant diseases. Although many molecular determinants involved in R. solanacearum adaptation to hosts and pathogenesis have been described, host components required for disease establishment remain poorly characterized. Phenotypical analysis of Arabidopsis mutants for leucine-rich repeat (LRR)-receptor-like proteins revealed that mutations in the CLAVATA1 (CLV1) and CLAVATA2 (CLV2) genes confer enhanced disease resistance to bacterial wilt. We further investigated the underlying mechanisms using genetic, transcriptomic and molecular approaches. The enhanced resistance of both clv1 and clv2 mutants to the bacteria did not require the well characterized CLV signalling modules involved in shoot meristem homeostasis, and was conditioned by neither salicylic acid nor ethylene defence-related hormones. Gene expression microarray analysis performed on clv1 and clv2 revealed deregulation of genes encoding nuclear transcription factor Y subunit alpha (NF-YA) transcription factors whose post-transcriptional regulation is known to involve microRNAs from the miR169 family. Both clv mutants showed a defect in miR169 accumulation. Conversely, overexpression of miR169 abrogated the resistance phenotype of clv mutants. We propose that CLV1 and CLV2, two receptors involved in CLV3 perception during plant development, contribute to bacterial wilt through a signalling pathway involving the miR169/NF-YA module.

EXPLOITING SENTINEL-1 AMPLITUDE DATA FOR GLACIER SURFACE VELOCITY FIELD MEASUREMENTS: FEASIBILITY DEMONSTRATION ON BALTORO GLACIER

The leading idea of this work is to continuously retrieve glaciers surface velocity through SAR imagery, in particular using the amplitude data from the new ESA satellite sensor Sentinel-1 imagery. These imagery key aspects are the free access policy, the very short revisit time (down to 6 days with the launch of the Sentinel-1B satellite) and the high amplitude resolution (up to 5 m). In order to verify the reliability of the proposed approach, a first experiment has been performed using Sentinel-1 imagery acquired over the Karakoram mountain range (North Pakistan) and Baltoro and other three glaciers have been investigated. During this study, a stack of 11 images acquired in the period from October 2014 to September 2015 has been used in order to investigate the potentialities of the Sentinel-1 SAR sensor to retrieve the glacier surface velocity every month. The aim of this test was to measure the glacier surface velocity between each subsequent pair, in order to produce a time series of the surface velocity fields along the investigated period. The necessary coregistration procedure between the images has been performed and subsequently the glaciers areas have been sampled using a regular grid with a 250 × 250 meters posting. Finally the surface velocity field has been estimated, for each image pair, using a template matching procedure, and an outlier filtering procedure based on the signal to noise ratio values has been applied, in order to exclude from the analysis unreliable points. The achieved velocity values range from 10 to 25 meters/month and they are coherent to those obtained in previous studies carried out on the same glaciers and the results highlight that it is possible to have a continuous update of the glacier surface velocity field through free Sentinel-1 imagery, that could be very useful to investigate the seasonal effects on the glaciers fluid-dynamics.