Multi-transcriptomics identifies targets of the endoribonuclease DNE1 and highlights its coordination with decapping.

Decapping is a crucial step in mRNA degradation in eucaryotes and requires the formation of a holoenzyme complex between the decapping enzyme DECAPPING 2 (DCP2) and the decapping enhancer DCP1. In Arabidopsis (Arabidopsis thaliana), DCP1-ASSOCIATED NYN ENDORIBONUCLEASE 1 (DNE1) is a direct protein partner of DCP1. The function of both DNE1 and decapping are necessary to maintain phyllotaxis, the regularity of organ emergence in the apex. In this study, we combined in vivo mRNA editing, RNA degradome sequencing, transcriptomics and small RNA-omics to identify targets of DNE1 and study how DNE1 and DCP2 cooperate in controlling mRNA fate. Our data reveal that DNE1 mainly contacts and cleaves mRNAs in the coding sequence and has sequence cleavage preferences. DNE1 targets are also degraded through decapping, and both RNA degradation pathways influence the production of mRNA-derived small interfering RNAs. Finally, we detected mRNA features enriched in DNE1 targets including RNA G-quadruplexes and translated upstream open reading frames. Combining these four complementary high-throughput sequencing strategies greatly expands the range of DNE1 targets and allowed us to build a conceptual framework describing the influence of DNE1 and decapping on mRNA fate. These data will be crucial to unveil the specificity of DNE1 action and understand its importance for developmental patterning.

Genome evolution across 1,011 Saccharomyces cerevisiae isolates.

Large-scale population genomic surveys are essential to explore the phenotypic diversity of natural populations. Here we report the whole-genome sequencing and phenotyping of 1,011 Saccharomyces cerevisiae isolates, which together provide an accurate evolutionary picture of the genomic variants that shape the species-wide phenotypic landscape of this yeast. Genomic analyses support a single 'out-of-China' origin for this species, followed by several independent domestication events. Although domesticated isolates exhibit high variation in ploidy, aneuploidy and genome content, genome evolution in wild isolates is mainly driven by the accumulation of single nucleotide polymorphisms. A common feature is the extensive loss of heterozygosity, which represents an essential source of inter-individual variation in this mainly asexual species. Most of the single nucleotide polymorphisms, including experimentally identified functional polymorphisms, are present at very low frequencies. The largest numbers of variants identified by genome-wide association are copy-number changes, which have a greater phenotypic effect than do single nucleotide polymorphisms. This resource will guide future population genomics and genotype-phenotype studies in this classic model system.

PlantRNA 2.0: an updated database dedicated to tRNAs of photosynthetic eukaryotes.

PlantRNA (http://plantrna.ibmp.cnrs.fr/) is a comprehensive database of transfer RNA (tRNA) gene sequences retrieved from fully annotated nuclear, plastidial and mitochondrial genomes of photosynthetic organisms. In the first release (PlantRNA 1.0), tRNA genes from 11 organisms were annotated. In this second version, the annotation was implemented to 51 photosynthetic species covering the whole phylogenetic tree of photosynthetic organisms, from the most basal group of Archeplastida, the glaucophyte Cyanophora paradoxa, to various land plants. tRNA genes from lower photosynthetic organisms such as streptophyte algae or lycophytes as well as extremophile photosynthetic species such as Eutrema parvulum were incorporated in the database. As a whole, about 37 000 tRNA genes were accurately annotated. In the frame of the tRNA genes annotation from the genome of the Rhodophyte Chondrus crispus, non-canonical splicing sites in the D- or T-regions of tRNA molecules were identified and experimentally validated. As for PlantRNA 1.0, comprehensive biological information including 5'- and 3'-flanking sequences, A and B box sequences, region of transcription initiation and poly(T) transcription termination stretches, tRNA intron sequences and tRNA mitochondrial import are included.

Jasmonate signaling controls negative and positive effectors of salt stress tolerance in rice.

Plant responses to salt exposure involve large reconfigurations of hormonal pathways that orchestrate physiological changes towards tolerance. Jasmonate (JA) hormones are essential to withstand biotic and abiotic assaults, but their roles in salt tolerance remain unclear. Here we describe the dynamics of JA metabolism and signaling in root and leaf tissue of rice, a plant species that is highly exposed and sensitive to salt. Roots activate the JA pathway in an early pulse, while the second leaf displays a biphasic JA response with peaks at 1 h and 3 d post-exposure. Based on higher salt tolerance of a rice JA-deficient mutant (aoc), we examined, through kinetic transcriptome and physiological analysis, the salt-triggered processes that are under JA control. Profound genotype-differential features emerged that could underlie the observed phenotypes. Abscisic acid (ABA) content and ABA-dependent water deprivation responses were impaired in aoc shoots. Moreover, aoc accumulated more Na+ in roots, and less in leaves, with reduced ion translocation correlating with root derepression of the HAK4 Na+ transporter gene. Distinct reactive oxygen species scavengers were also stronger in aoc leaves, along with reduced senescence and chlorophyll catabolism markers. Collectively, our results identify contrasted contributions of JA signaling to different sectors of the salt stress response in rice.

Epigenetic silencing of clustered tRNA genes in Arabidopsis.

Beyond their key role in translation, cytosolic transfer RNAs (tRNAs) are involved in a wide range of other biological processes. Nuclear tRNA genes (tDNAs) are transcribed by the RNA polymerase III (RNAP III) and cis-elements, trans-factors as well as genomic features are known to influence their expression. In Arabidopsis, besides a predominant population of dispersed tDNAs spread along the 5 chromosomes, some clustered tDNAs have been identified. Here, we demonstrate that these tDNA clusters are transcriptionally silent and that pathways involved in the maintenance of DNA methylation play a predominant role in their repression. Moreover, we show that clustered tDNAs exhibit repressive chromatin features whilst their dispersed counterparts contain permissive euchromatic marks. This work demonstrates that both genomic and epigenomic contexts are key players in the regulation of tDNAs transcription. The conservation of most of these regulatory processes suggests that this pioneering work in Arabidopsis can provide new insights into the regulation of RNA Pol III transcription in other organisms, including vertebrates.

Pervasive Phenotypic Impact of a Large Nonrecombining Introgressed Region in Yeast.

To explore the origin of the diversity observed in natural populations, many studies have investigated the relationship between genotype and phenotype. In yeast species, especially in Saccharomyces cerevisiae, these studies are mainly conducted using recombinant offspring derived from two genetically diverse isolates, allowing to define the phenotypic effect of genetic variants. However, large genomic variants such as interspecies introgressions are usually overlooked even if they are known to modify the genotype-phenotype relationship. To have a better insight into the overall phenotypic impact of introgressions, we took advantage of the presence of a 1-Mb introgressed region, which lacks recombination and contains the mating-type determinant in the Lachancea kluyveri budding yeast. By performing linkage mapping analyses in this species, we identified a total of 89 loci affecting growth fitness in a large number of conditions and 2,187 loci affecting gene expression mostly grouped into two major hotspots, one being the introgressed region carrying the mating-type locus. Because of the absence of recombination, our results highlight the presence of a sexual dimorphism in a budding yeast for the first time. Overall, by describing the phenotype-genotype relationship in the Lachancea kluyveri species, we expanded our knowledge on how genetic characteristics of large introgression events can affect the phenotypic landscape.

The NRPD1 N-terminus contains a Pol IV-specific motif that is critical for genome surveillance in Arabidopsis.

RNA-guided surveillance systems constrain the activity of transposable elements (TEs) in host genomes. In plants, RNA polymerase IV (Pol IV) transcribes TEs into primary transcripts from which RDR2 synthesizes double-stranded RNA precursors for small interfering RNAs (siRNAs) that guide TE methylation and silencing. How the core subunits of Pol IV, homologs of RNA polymerase II subunits, diverged to support siRNA biogenesis in a TE-rich, repressive chromatin context is not well understood. Here we studied the N-terminus of Pol IV's largest subunit, NRPD1. Arabidopsis lines harboring missense mutations in this N-terminus produce wild-type (WT) levels of NRPD1, which co-purifies with other Pol IV subunits and RDR2. Our in vitro transcription and genomic analyses reveal that the NRPD1 N-terminus is critical for robust Pol IV-dependent transcription, siRNA production and DNA methylation. However, residual RNA-directed DNA methylation observed in one mutant genotype indicates that Pol IV can operate uncoupled from the high siRNA levels typically observed in WT plants. This mutation disrupts a motif uniquely conserved in Pol IV, crippling the enzyme's ability to inhibit retrotransposon mobilization. We propose that the NRPD1 N-terminus motif evolved to regulate Pol IV function in genome surveillance.

Hybrid incompatibility caused by an epiallele.

Hybrid incompatibility resulting from deleterious gene combinations is thought to be an important step toward reproductive isolation and speciation. Here, we demonstrate involvement of a silent epiallele in hybrid incompatibility. In Arabidopsis thaliana accession Cvi-0, one of the two copies of a duplicated histidine biosynthesis gene, HISN6A, is mutated, making HISN6B essential. In contrast, in accession Col-0, HISN6A is essential because HISN6B is not expressed. Owing to these differences, Cvi-0 × Col-0 hybrid progeny that are homozygous for both Cvi-0 HISN6A and Col-0 HISN6B do not survive. We show that HISN6B of Col-0 is not a defective pseudogene, but a stably silenced epiallele. Mutating HISTONE DEACETYLASE 6 (HDA6), or the cytosine methyltransferase genes MET1 or CMT3, erases HISN6B's silent locus identity, reanimating the gene to circumvent hisn6a lethality and hybrid incompatibility. These results show that HISN6-dependent hybrid lethality is a revertible epigenetic phenomenon and provide additional evidence that epigenetic variation has the potential to limit gene flow between diverging populations of a species.

Variation of the meiotic recombination landscape and properties over a broad evolutionary distance in yeasts.

Meiotic recombination is a major factor of genome evolution, deeply characterized in only a few model species, notably the yeast Saccharomyces cerevisiae. Consequently, little is known about variations of its properties across species. In this respect, we explored the recombination landscape of Lachancea kluyveri, a protoploid yeast species that diverged from the Saccharomyces genus more than 100 million years ago and we found striking differences with S. cerevisiae. These variations include a lower recombination rate, a higher frequency of chromosomes segregating without any crossover and the absence of recombination on the chromosome arm containing the sex locus. In addition, although well conserved within the Saccharomyces clade, the S. cerevisiae recombination hotspots are not conserved over a broader evolutionary distance. Finally and strikingly, we found evidence of frequent reversal of commitment to meiosis, resulting in return to mitotic growth after allele shuffling. Identification of this major but underestimated evolutionary phenomenon illustrates the relevance of exploring non-model species.

Evolution of intraspecific transcriptomic landscapes in yeasts.

Variations in gene expression have been widely explored in order to obtain an accurate overview of the changes in regulatory networks that underlie phenotypic diversity. Numerous studies have characterized differences in genomic expression between large numbers of individuals of model organisms such as Saccharomyces cerevisiae. To more broadly survey the evolution of the transcriptomic landscape across species, we measured whole-genome expression in a large collection of another yeast species: Lachancea kluyveri (formerly Saccharomyces kluyveri), using RNAseq. Interestingly, this species diverged from the S. cerevisiae lineage prior to its ancestral whole genome duplication. Moreover, L. kluyveri harbors a chromosome-scale compositional heterogeneity due to a 1-Mb ancestral introgressed region as well as a large set of unique unannotated genes. In this context, our comparative transcriptomic analysis clearly showed a link between gene evolutionary history and expression behavior. Indeed, genes that have been recently acquired or under function relaxation tend to be less transcribed show a higher intraspecific variation (plasticity) and are less involved in network (connectivity). Moreover, utilizing this approach in L. kluyveri also highlighted specific regulatory network signatures in aerobic respiration, amino-acid biosynthesis and glycosylation, presumably due to its different lifestyle. Our data set sheds an important light on the evolution of intraspecific transcriptomic variation across distant species.

Evolutionary metabolomics of specialized metabolism diversification in the genus Nicotiana highlights N-acylnornicotine innovations.

Specialized metabolite (SM) diversification is a core process to plants' adaptation to diverse ecological niches. Here, we implemented a computational mass spectrometry-based metabolomics approach to exploring SM diversification in tissues of 20 species covering Nicotiana phylogenetics sections. To markedly increase metabolite annotation, we created a large in silico fragmentation database, comprising >1 million structures, and scripts for connecting class prediction to consensus substructures. Together, the approach provides an unprecedented cartography of SM diversity and section-specific innovations in this genus. As a case study and in combination with nuclear magnetic resonance and mass spectrometry imaging, we explored the distribution of N-acylnornicotines, alkaloids predicted to be specific to Repandae allopolyploids, and revealed their prevalence in the genus, albeit at much lower magnitude, as well as a greater structural diversity than previously thought. Together, the data integration approaches provided here should act as a resource for future research in plant SM evolution.

The TUTase URT1 connects decapping activators and prevents the accumulation of excessively deadenylated mRNAs to avoid siRNA biogenesis.

Uridylation is a widespread modification destabilizing eukaryotic mRNAs. Yet, molecular mechanisms underlying TUTase-mediated mRNA degradation remain mostly unresolved. Here, we report that the Arabidopsis TUTase URT1 participates in a molecular network connecting several translational repressors/decapping activators. URT1 directly interacts with DECAPPING 5 (DCP5), the Arabidopsis ortholog of human LSM14 and yeast Scd6, and this interaction connects URT1 to additional decay factors like DDX6/Dhh1-like RNA helicases. Nanopore direct RNA sequencing reveals a global role of URT1 in shaping poly(A) tail length, notably by preventing the accumulation of excessively deadenylated mRNAs. Based on in vitro and in planta data, we propose a model that explains how URT1 could reduce the accumulation of oligo(A)-tailed mRNAs both by favoring their degradation and because 3' terminal uridines intrinsically hinder deadenylation. Importantly, preventing the accumulation of excessively deadenylated mRNAs avoids the biogenesis of illegitimate siRNAs that silence endogenous mRNAs and perturb Arabidopsis growth and development.

Integrated Genome-Scale Analysis and Northern Blot Detection of Retrotransposon siRNAs Across Plant Species.

Cells have sophisticated RNA-directed mechanisms to regulate genes, destroy viruses, or silence transposable elements (TEs). In terrestrial plants, a specialized non-coding RNA machinery involving RNA polymerase IV (Pol IV) and small interfering RNAs (siRNAs) targets DNA methylation and silencing to TEs. Here, we present a bioinformatics protocol for annotating and quantifying siRNAs that derive from long terminal repeat (LTR) retrotransposons. The approach was validated using small RNA northern blot analyses, comparing the species Arabidopsis thaliana and Brachypodium distachyon. To assist hybridization probe design, we configured a genome browser to show small RNA-seq mappings in distinct colors and shades according to their nucleotide lengths and abundances, respectively. Samples from wild-type and pol IV mutant plants, cross-species negative controls, and a conserved microRNA control validated the detected siRNA signals, confirming their origin from specific TEs and their Pol IV-dependent biogenesis. Moreover, an optimized labeling method yielded probes that could detect low-abundance siRNAs from B. distachyon TEs. The integration of de novo TE annotation, small RNA-seq profiling, and northern blotting, as outlined here, will facilitate the comparative genomic analysis of RNA silencing in crop plants and non-model species.

Turnip mosaic virus in oilseed rape activates networks of sRNA-mediated interactions between viral and host genomes.

Virus-induced plant diseases in cultivated plants cause important damages in yield. Although the mechanisms of virus infection are intensely studied at the cell biology level, only little is known about the molecular dialog between the invading virus and the host genome. Here we describe a combinatorial genome-wide approach to identify networks of sRNAs-guided post-transcriptional regulation within local Turnip mosaic virus (TuMV) infection sites in Brassica napus leaves. We show that the induction of host-encoded, virus-activated small interfering RNAs (vasiRNAs) observed in virus-infected tissues is accompanied by site-specific cleavage events on both viral and host RNAs that recalls the activity of small RNA-induced silencing complexes (RISC). Cleavage events also involve virus-derived siRNA (vsiRNA)-directed cleavage of target host transcripts as well as cleavage of viral RNA by both host vasiRNAs and vsiRNAs. Furthermore, certain coding genes act as virus-activated regulatory hubs to produce vasiRNAs for the targeting of other host genes. The observations draw an advanced model of plant-virus interactions and provide insights into the complex regulatory networking at the plant-virus interface within cells undergoing early stages of infection.

RST1 and RIPR connect the cytosolic RNA exosome to the Ski complex in Arabidopsis.

The RNA exosome is a key 3'-5' exoribonuclease with an evolutionarily conserved structure and function. Its cytosolic functions require the co-factors SKI7 and the Ski complex. Here we demonstrate by co-purification experiments that the ARM-repeat protein RESURRECTION1 (RST1) and RST1 INTERACTING PROTEIN (RIPR) connect the cytosolic Arabidopsis RNA exosome to the Ski complex. rst1 and ripr mutants accumulate RNA quality control siRNAs (rqc-siRNAs) produced by the post-transcriptional gene silencing (PTGS) machinery when mRNA degradation is compromised. The small RNA populations observed in rst1 and ripr mutants are also detected in mutants lacking the RRP45B/CER7 core exosome subunit. Thus, molecular and genetic evidence supports a physical and functional link between RST1, RIPR and the RNA exosome. Our data reveal the existence of additional cytosolic exosome co-factors besides the known Ski subunits. RST1 is not restricted to plants, as homologues with a similar domain architecture but unknown function exist in animals, including humans.

Differences in environmental stress response among yeasts is consistent with species-specific lifestyles.

Defining how organisms respond to environmental change has always been an important step toward understanding their adaptive capacity and physiology. Variation in transcription during stress has been widely described in model species, especially in the yeast Saccharomyces cerevisiae, which helped to shape general rules regarding how cells cope with environmental constraints, as well as to decipher the functions of many genes. Comparison of the environmental stress response (ESR) across species is essential to obtaining better insight into the common and species-specific features of stress defense. In this context, we explored the transcriptional landscape of the yeast Lachancea kluyveri (formerly Saccharomyces kluyveri) in response to diverse stresses, using RNA sequencing. We investigated variation in gene expression and observed a link between genetic plasticity and environmental sensitivity. We identified the ESR genes in this species and compared them to those already found in S. cerevisiae We observed common features between the two species, as well as divergence in the regulatory networks involved. Of interest, some changes were related to differences in species lifestyle. Thus we were able to decipher how adaptation to stress has evolved among different yeast species. Finally, by analyzing patterns of coexpression, we were able to propose potential biological functions for 42% of genes and also annotate 301 genes for which no function could be assigned by homology. This large data set allowed for the characterization of the evolution of gene regulation and provides an efficient tool for assessing gene function.

The Hidden Complexity of Mendelian Traits across Natural Yeast Populations.

Mendelian traits are considered to be at the lower end of the complexity spectrum of heritable phenotypes. However, more than a century after the rediscovery of Mendel's law, the global landscape of monogenic variants, as well as their effects and inheritance patterns within natural populations, is still not well understood. Using the yeast Saccharomyces cerevisiae, we performed a species-wide survey of Mendelian traits across a large population of isolates. We generated offspring from 41 unique parental pairs and analyzed 1,105 cross/trait combinations. We found that 8.9% of the cases were Mendelian. Further tracing of causal variants revealed background-specific expressivity and modified inheritances, gradually transitioning from Mendelian to complex traits in 30% of the cases. In fact, when taking into account the natural population diversity, the hidden complexity of traits could be substantial, confounding phenotypic predictability even for simple Mendelian traits.