Endogenous RNAi pathway evolutionarily shapes the destiny of the antisense lncRNAs transcriptome

Pervasive formation of double-stranded RNAs by overlapping sense/antisense transcripts in budding yeast mitosis and meiosis

Previous RNA profiling studies revealed coexpression of overlapping sense/antisense (s/a) transcripts in pro- and eukaryotic organisms. Functional analyses in yeast have shown that certain s/a mRNA/mRNA and mRNA/lncRNA pairs form stable double-stranded RNAs (dsRNAs) that affect transcript stability. Little is known, however, about the genome-wide prevalence of dsRNA formation and its potential functional implications during growth and development in diploid budding yeast. To address this question, we monitored dsRNAs in a Saccharomyces cerevisiae strain expressing the ribonuclease DCR1 and the RNA-binding protein AGO1 from Naumovozyma castellii We identify dsRNAs at 347 s/a loci that express partially or completely overlapping transcripts during mitosis, meiosis, or both stages of the diploid life cycle. We associate dsRNAs with s/a loci previously thought to be exclusively regulated by antisense interference, and others that encode antisense RNAs, which down-regulate sense mRNA-encoded protein levels. To facilitate hypothesis building, we developed the sense/antisense double-stranded RNA (SensR) expression viewer. Users are able to retrieve different graphical displays of dsRNA and RNA expression data using genome coordinates and systematic or standard names for mRNAs and different types of stable or cryptic long noncoding RNAs (lncRNAs). Our data are a useful resource for improving yeast genome annotation and for work on RNA-based regulatory mechanisms controlling transcript and protein levels. The data are also interesting from an evolutionary perspective, since natural antisense transcripts that form stable dsRNAs have been detected in many species from bacteria to humans. The SensR viewer is freely accessible at https://sensr.genouest.org.

Northern Blotting: Protocols for Radioactive and Nonradioactive Detection of RNA

Northern blotting is a common technique in RNA biology, allowing to detect and quantify RNAs of interest following separation by gel electrophoresis, transfer to a membrane, and hybridization of specific anti-complementary labelled probes. In this chapter, we describe our protocol for efficient RNA extraction from yeast, separation on agarose gel, and capillary transfer to a membrane. We provide two different methods for strand-specific detection of several types of RNAs using oligonucleotide probes, the first using radioactive 32P-labelled probes, the second based on nonradioactive digoxigenin-labelled probes.

Transcriptome-Wide Analysis of the 5' Cap Status of RNA Using 5' Monophosphate-Dependent Exonuclease Digestion and RNA Sequencing

Eukaryotic mRNAs carry an N7-methylguanosine (m7G) cap structure at their 5' extremity, which protects them from the degradation by 5'-3' exoribonucleases and plays a pivotal role in mRNA metabolism, promoting splicing, nuclear export, and translation. Decapping, the enzymatic process that removes this structure, is a key event during cytoplasmic mRNA 5'-3' decay, leading to the degradation of the transcript body by Xrn1. In this chapter, we describe a procedure to assess the cap status of RNA at the transcriptome level. It is based on a treatment of total RNA extracts with a 5' monophosphate-dependent exonuclease, which like Xrn1 specifically degrades decapped RNAs harboring 5' monophosphate extremities, but not RNAs with intact m7G cap. The digested RNAs are then analyzed by RNA sequencing.

The 5-Fluorouracil RNA Expression Viewer (5-FUR) Facilitates Interpreting the Effects of Drug Treatment and RRP6 Deletion on the Transcriptional Landscape in Yeast

Saccharomyces cerevisiae is an excellent model to study the effect of external cues on cell division and stress response. 5-Fluorocuracil (5-FU) has been used to treat solid tumors since several decades. The drug was initially designed to interfere with DNA replication but was later found to exert its antiproliferative effect also via RNA-dependent processes. Since 5-FU inhibits the activity of the 3'-5'-exoribonuclease Rrp6 in yeast and mammals, earlier work has compared the effect of 5-FU treatment and RRP6 deletion at the transcriptome level in diploid synchronized yeast cells. To facilitate interpreting the expression data we have developed an improved 5-Fluorouracil RNA (5-FUR) expression viewer. Users can access information via genome coordinates and systematic or standard names for mRNAs and Xrn1-dependent-, stable-, cryptic-, and meiotic unannotated transcripts (XUTs, SUTs, CUTs, and MUTs). Normalized log2-transformed or linear data can be displayed as filled diagrams, line graphs or color-coded heatmaps. The expression data are useful for researchers interested in processes such as cell cycle regulation, mitotic repression of meiotic genes, the effect of 5-FU treatment and Rrp6 deficiency on the transcriptome and expression profiles of sense/antisense loci that encode overlapping transcripts. The viewer is accessible at http://5fur.genouest.org.

Pervasive translation of Xrn1-sensitive unstable long noncoding RNAs in yeast

Despite being predicted to lack coding potential, cytoplasmic long noncoding (lnc)RNAs can associate with ribosomes. However, the landscape and biological relevance of lncRNA translation remain poorly studied. In yeast, cytoplasmic Xrn1-sensitive unstable transcripts (XUTs) are targeted by nonsense-mediated mRNA decay (NMD), suggesting a translation-dependent degradation process. Here, we report that XUTs are pervasively translated, which impacts their decay. We show that XUTs globally accumulate upon translation elongation inhibition, but not when initial ribosome loading is impaired. Ribo-seq confirmed ribosomes binding to XUTs and identified ribosome-associated 5'-proximal small ORFs. Mechanistically, the NMD-sensitivity of XUTs mainly depends on the 3'-untranslated region length. Finally, we show that the peptide resulting from the translation of an NMD-sensitive XUT reporter exists in NMD-competent cells. Our work highlights the role of translation in the posttranscriptional metabolism of XUTs. We propose that XUT-derived peptides could be exposed to natural selection, while NMD restricts XUT levels.

Genome-wide transcriptome analysis reveals the diversity and function of long non-coding RNAs in dinoflagellates

Dinoflagellates are a diverse group of phytoplankton, ranging from harmful bloom-forming microalgae to photosymbionts of coral reefs. Genome-scale data from dinoflagellates reveal atypical genomic features, extensive genomic divergence, and lineage-specific innovation of gene functions. Long non-coding RNAs (lncRNAs), known to regulate gene expression in eukaryotes, are largely unexplored in dinoflagellates. Here, using high-quality genome and transcriptome data, we identified 48039 polyadenylated lncRNAs in three dinoflagellate species: the coral symbionts Cladocopium proliferum and Durusdinium trenchii, and the bloom-forming species, Prorocentrum cordatum. These lncRNAs have fewer introns and lower G+C content than protein-coding sequences; 37 768 (78.6%) are unique with respect to sequence similarity. We classified all lncRNAs based on conserved motifs (k-mers) into distinct clusters, following properties of protein-binding and/or subcellular localisation. Interestingly, 3708 (7.7%) lncRNAs are differentially expressed under heat stress, algal lifestyle, and/or growth phase, and share co-expression patterns with protein-coding genes. Based on inferred triplex interactions between lncRNA and putative promoter regions, we identified 19 460 putative gene targets for 3721 lncRNAs; 907 genes exhibit differential expression under heat stress. These results reveal, for the first time, the diversity of lncRNAs in dinoflagellates and how lncRNAs may regulate gene expression as a heat-stress response in these ecologically important microbes.

Parallel Nonfunctionalization of CK1δ/ε Kinase Ohnologs Following a Whole-Genome Duplication Event

Whole-genome duplication (WGD) followed by speciation allows us to examine the parallel evolution of ohnolog pairs. In the yeast family Saccharomycetaceae, HRR25 is a rare case of repeated ohnolog maintenance. This gene has reverted to a single copy in Saccharomyces cerevisiae where it is now essential, but has been maintained as pairs in at least 7 species post-WGD. In S. cerevisiae, HRR25 encodes the casein kinase 1δ/ε and plays a role in a variety of functions through its kinase activity and protein-protein interactions (PPIs). We hypothesized that the maintenance of duplicated HRR25 ohnologs could be a result of repeated subfunctionalization. We tested this hypothesis through a functional complementation assay in S. cerevisiae, testing all pairwise combinations of 25 orthologs (including 7 ohnolog pairs). Contrary to our expectations, we observed no cases of pair-dependent complementation, which would have supported the subfunctionalization hypothesis. Instead, most post-WGD species have one ohnolog that failed to complement, suggesting their nonfunctionalization or neofunctionalization. The ohnologs incapable of complementation have undergone more rapid protein evolution, lost most PPIs that were observed for their functional counterparts and singletons from post-WGD and non-WGD species, and have nonconserved cellular localization, consistent with their ongoing loss of function. The analysis in Naumovozyma castellii shows that the noncomplementing ohnolog is expressed at a lower level and has become nonessential. Taken together, our results indicate that HRR25 orthologs are undergoing gradual nonfunctionalization.

Parallel nonfunctionalization of CK1δ/ε kinase ohnologs following a whole-genome duplication event

Whole genome duplication (WGD) followed by speciation allows us to examine the parallel evolution of ohnolog pairs. In the yeast family Saccharomycetaceae, HRR25 is a rare case of repeated ohnolog maintenance. This gene has reverted to a single copy in S. cerevisiae where it is now essential, but has been maintained as pairs in at least 7 species post WGD. In S. cerevisiae, HRR25 encodes the casein kinase (CK) 1δ/ε and plays a role in a variety of functions through its kinase activity and protein-protein interactions (PPIs). We hypothesized that the maintenance of duplicated HRR25 ohnologs could be a result of repeated subfunctionalization. We tested this hypothesis through a functional complementation assay in S. cerevisiae, testing all pairwise combinations of 25 orthologs (including 7 ohnolog pairs). Contrary to our expectations, we observed no cases of pair-dependent complementation, which would have supported the subfunctionalization hypothesis. Instead, most post-WGD species have one ohnolog that failed to complement, suggesting their nonfunctionalization or neofunctionalization. The ohnologs incapable of complementation have undergone more rapid protein evolution, lost most PPIs that were observed for their functional counterparts and singletons from post and non-WGD species, and have non-conserved cellular localization, consistent with their ongoing loss of function. The analysis in N. castelli shows that the non-complementing ohnolog is expressed at a lower level and has become non-essential. Taken together, our results indicate that HRR25 orthologs are undergoing gradual nonfunctionalization.

The RNA helicases Dbp2 and Mtr4 regulate the expression of Xrn1-sensitive long non-coding RNAs in yeast

The expression of yeast long non-coding (lnc)RNAs is restricted by RNA surveillance machineries, including the cytoplasmic 5'-3' exonuclease Xrn1 which targets a conserved family of lncRNAs defined as XUTs, and that are mainly antisense to protein-coding genes. However, the co-factors involved in the degradation of these transcripts and the underlying molecular mechanisms remain largely unknown. Here, we show that two RNA helicases, Dbp2 and Mtr4, act as global regulators of XUTs expression. Using RNA-Seq, we found that most of them accumulate upon Dbp2 inactivation or Mtr4 depletion. Mutants of the cytoplasmic RNA helicases Ecm32, Ski2, Slh1, Dbp1, and Dhh1 did not recapitulate this global stabilization of XUTs, suggesting that XUTs decay is specifically controlled by Dbp2 and Mtr4. Notably, Dbp2 and Mtr4 affect XUTs independently of their configuration relative to their paired-sense mRNAs. Finally, we show that the effect of Dbp2 on XUTs depends on a cytoplasmic localization. Overall, our data indicate that Dbp2 and Mtr4 are global regulators of lncRNAs expression and contribute to shape the non-coding transcriptome together with RNA decay machineries.

RNA-Mediated Regulation of Meiosis in Budding Yeast

Cells change their physiological state in response to environmental cues. In the absence of nutrients, unicellular fungi such as budding yeast exit mitotic proliferation and enter the meiotic cycle, leading to the production of haploid cells that are encased within spore walls. These cell state transitions are orchestrated in a developmentally coordinated manner. Execution of the meiotic cell cycle program in budding yeast, Saccharomyces cerevisiae, is regulated by the key transcription factor, Ime1. Recent developments have uncovered the role of non-coding RNA in the regulation of Ime1 and meiosis. In this review, we summarize the role of ncRNA-mediated and RNA homeostasis-based processes in the regulation of meiosis in Saccharomyces cerevisiae.

From Yeast to Mammals, the Nonsense-Mediated mRNA Decay as a Master Regulator of Long Non-Coding RNAs Functional Trajectory

The Nonsense-Mediated mRNA Decay (NMD) has been classically viewed as a translation-dependent RNA surveillance pathway degrading aberrant mRNAs containing premature stop codons. However, it is now clear that mRNA quality control represents only one face of the multiple functions of NMD. Indeed, NMD also regulates the physiological expression of normal mRNAs, and more surprisingly, of long non-coding (lnc)RNAs. Here, we review the different mechanisms of NMD activation in yeast and mammals, and we discuss the molecular bases of the NMD sensitivity of lncRNAs, considering the functional roles of NMD and of translation in the metabolism of these transcripts. In this regard, we describe several examples of functional micropeptides produced from lncRNAs. We propose that translation and NMD provide potent means to regulate the expression of lncRNAs, which might be critical for the cell to respond to environmental changes.

Xrn1p acts at multiple steps in the budding-yeast RNAi pathway to enhance the efficiency of silencing

RNA interference (RNAi) is a gene-silencing pathway that can play roles in viral defense, transposon silencing, heterochromatin formation and post-transcriptional gene silencing. Although absent from Saccharomyces cerevisiae, RNAi is present in other budding-yeast species, including Naumovozyma castellii, which have an unusual Dicer and a conventional Argonaute that are both required for gene silencing. To identify other factors that act in the budding-yeast pathway, we performed an unbiased genetic selection. This selection identified Xrn1p, the cytoplasmic 5'-to-3' exoribonuclease, as a cofactor of RNAi in budding yeast. Deletion of XRN1 impaired gene silencing in N. castellii, and this impaired silencing was attributable to multiple functions of Xrn1p, including affecting the composition of siRNA species in the cell, influencing the efficiency of siRNA loading into Argonaute, degradation of cleaved passenger strand and degradation of sliced target RNA.

Nuclear Argonaute Piwi Gene Mutation Affects rRNA by Inducing rRNA Fragment Accumulation, Antisense Expression, and Defective Processing in Drosophila Ovaries

Drosophila key nuclear piRNA silencing pathway protein Piwi of the Argonaute family has been classically studied as a factor controlling transposable elements and fertility. Piwi has been shown to concentrate in the nucleolus for reasons largely unknown. Ribosomal RNA is the main component of the nucleolus. In this work the effect of a piwi mutation on rRNA is described. This work led to three important conclusions: A mutation in piwi induces antisense 5S rRNA expression, a processing defect of 2S rRNA orthologous to the 3′-end of eukaryotic 5.8S rRNA, and accumulation of fragments of all five rRNAs in Drosophilamelanogaster ovaries. Hypotheses to explain these phenomena are proposed, possibly involving the interaction of the components of the piRNA pathway with the RNA surveillance machinery.