Protein Abundance Control by Non-coding Antisense Transcription

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.

dsRNA formation leads to preferential nuclear export and gene expression

When mRNAs have been transcribed and processed in the nucleus, they are exported to the cytoplasm for translation. This export is mediated by the export receptor heterodimer Mex67-Mtr2 in the yeast Saccharomyces cerevisiae (TAP-p15 in humans)1,2. Interestingly, many long non-coding RNAs (lncRNAs) also leave the nucleus but it is currently unclear why they move to the cytoplasm3. Here we show that antisense RNAs (asRNAs) accelerate mRNA export by annealing with their sense counterparts through the helicase Dbp2. These double-stranded RNAs (dsRNAs) dominate export compared with single-stranded RNAs (ssRNAs) because they have a higher capacity and affinity for the export receptor Mex67. In this way, asRNAs boost gene expression, which is beneficial for cells. This is particularly important when the expression program changes. Consequently, the degradation of dsRNA, or the prevention of its formation, is toxic for cells. This mechanism illuminates the general cellular occurrence of asRNAs and explains their nuclear export.

Post-transcriptional regulation shapes the transcriptome of quiescent budding yeast

To facilitate long-term survival, cells must exit the cell cycle and enter quiescence, a reversible non-replicative state. Budding yeast cells reprogram their gene expression during quiescence entry to silence transcription, but how the nascent transcriptome changes in quiescence is unknown. By investigating the nascent transcriptome, we identified over a thousand noncoding RNAs in quiescent and G1 yeast cells, and found noncoding transcription represented a larger portion of the quiescent transcriptome than in G1. Additionally, both mRNA and ncRNA are subject to increased post-transcriptional regulation in quiescence compared to G1. We found that, in quiescence, the nuclear exosome-NNS pathway suppresses over one thousand mRNAs, in addition to canonical noncoding RNAs. RNA sequencing through quiescent entry revealed two distinct time points at which the nuclear exosome controls the abundance of mRNAs involved in protein production, cellular organization, and metabolism, thereby facilitating efficient quiescence entry. Our work identified a previously unknown key biological role for the nuclear exosome-NNS pathway in mRNA regulation and uncovered a novel layer of gene-expression control in quiescence.

Differential expression patterns of long noncoding RNAs in a pleiomorphic diatom and relation to hyposalinity

Long non-coding (lnc)RNAs have been shown to have central roles in stress responses, cell identity and developmental processes in multicellular organisms as well as in unicellular fungi. Previous works have shown the occurrence of lncRNAs in diatoms, namely in Phaeodactylum tricornutum, many of which being expressed under specific stress conditions. Interestingly, P. tricornutum is the only known diatom that has a demonstrated morphological plasticity, occurring in three distinct morphotypes: fusiform, triradiate and oval. Although the morphotypes are interchangeable, the fusiform is the dominant one while both the triradiate and the oval forms are less common, the latter often being associated with stress conditions such as low salinity and solid culture media, amongst others. Nonetheless, the molecular basis underpinning morphotype identity in P. tricornutum remains elusive. Using twelve previously published transcriptomic datasets originating from the three morphotypes of P. tricornutum, we sought to investigate the expression patterns of lncRNAs (lincRNAs and NATs) in these distinct morphotypes, using pairwise comparisons, in order to explore the putative involvement of these noncoding molecules in morphotype identity. We found that differentially expressed lncRNAs cluster according to morphotype, indicating that lncRNAs are not randomly expressed, but rather seem to provide a specific (noncoding) transcriptomic signature of the morphotype. We also present evidence to suggest that the major differences in DE genes (both noncoding and coding) between the stress related oval morphotype and the most common fusiform morphotype could be due, to a large extent, to the hyposaline culture conditions rather than to the morphotype itself. However, several lncRNAs associated to each one of the three morphotypes were identified, which could have a potential role in morphotype (or cell) identity in P. tricornutum, similar to what has been found in both animals and plant development.

The RNA-binding protein RBP33 dampens non-productive transcription in trypanosomes

In-depth analysis of the transcriptomes of several model organisms has revealed that genomes are pervasively transcribed, giving rise to an abundance of non-canonical and mainly antisense RNA polymerase II-derived transcripts that are produced from almost any genomic context. Pervasive RNAs are degraded by surveillance mechanisms, but the repertoire of proteins that control the fate of these non-productive transcripts is still incomplete. Trypanosomes are single-celled eukaryotes that show constitutive RNA polymerase II transcription and in which initiation and termination of transcription occur at a limited number of sites per chromosome. It is not known whether pervasive transcription exists in organisms with unregulated RNA polymerase II activity, and which factors could be involved in the process. We show here that depletion of RBP33 results in overexpression of ∼40% of all annotated genes in the genome, with a marked accumulation of sense and antisense transcripts derived from silenced regions. RBP33 loss does not result in a significant increase in chromatin accessibility. Finally, we have found that transcripts that increase in abundance upon RBP33 knockdown are significantly more stable in RBP33-depleted trypanosomes, and that the exosome complex is responsible for their degradation. Our results provide strong evidence that RBP33 dampens non-productive transcription in trypanosomes.

Antisense non-coding transcription represses the PHO5 model gene at the level of promoter chromatin structure

Pervasive transcription of eukaryotic genomes generates non-coding transcripts with regulatory potential. We examined the effects of non-coding antisense transcription on the regulation of expression of the yeast PHO5 gene, a paradigmatic case for gene regulation through promoter chromatin remodeling. A negative role for antisense transcription at the PHO5 gene locus was demonstrated by leveraging the level of overlapping antisense transcription through specific mutant backgrounds, expression from a strong promoter in cis, and use of the CRISPRi system. Furthermore, we showed that enhanced elongation of PHO5 antisense leads to a more repressive chromatin conformation at the PHO5 gene promoter, which is more slowly remodeled upon gene induction. The negative effect of antisense transcription on PHO5 gene transcription is mitigated upon inactivation of the histone deacetylase Rpd3, showing that PHO5 antisense RNA acts via histone deacetylation. This regulatory pathway leads to Rpd3-dependent decreased recruitment of the RSC chromatin remodeling complex to the PHO5 gene promoter upon induction of antisense transcription. Overall, the data in this work reveal an additional level in the complex regulatory mechanism of PHO5 gene expression by showing antisense transcription-mediated repression at the level of promoter chromatin structure remodeling.

Pervasive transcription: a controlled risk

Transcriptome-wide interrogation of eukaryotic genomes has unveiled the pervasive nature of RNA polymerase II transcription. Virtually, any DNA region with an accessible chromatin structure can be transcribed, resulting in a mass production of noncoding RNAs (ncRNAs) with the potential of interfering with gene expression programs. Budding yeast has proved to be a powerful model organism to understand the mechanisms at play to control pervasive transcription and overcome the risks of hazardous disruption of cellular functions. In this review, we focus on the actors and strategies yeasts employ to govern ncRNA production, and we discuss recent findings highlighting the dangers of losing control over pervasive transcription.

New data and collaborations at the Saccharomyces Genome Database: updated reference genome, alleles, and the Alliance of Genome Resources

Saccharomyces cerevisiae is used to provide fundamental understanding of eukaryotic genetics, gene product function, and cellular biological processes. Saccharomyces Genome Database (SGD) has been supporting the yeast research community since 1993, serving as its de facto hub. Over the years, SGD has maintained the genetic nomenclature, chromosome maps, and functional annotation, and developed various tools and methods for analysis and curation of a variety of emerging data types. More recently, SGD and six other model organism focused knowledgebases have come together to create the Alliance of Genome Resources to develop sustainable genome information resources that promote and support the use of various model organisms to understand the genetic and genomic bases of human biology and disease. Here we describe recent activities at SGD, including the latest reference genome annotation update, the development of a curation system for mutant alleles, and new pages addressing homology across model organisms as well as the use of yeast to study human disease.

Modulated termination of non-coding transcription partakes in the regulation of gene expression

Pervasive transcription is a universal phenomenon leading to the production of a plethora of non-coding RNAs. If left uncontrolled, pervasive transcription can be harmful for genome expression and stability. However, non-coding transcription can also play important regulatory roles, for instance by promoting the repression of specific genes by a mechanism of transcriptional interference. The efficiency of transcription termination can strongly influence the regulatory capacity of non-coding transcription events, yet very little is known about the mechanisms modulating the termination of non-coding transcription in response to environmental cues. Here, we address this question by investigating the mechanisms that regulate the activity of the main actor in termination of non-coding transcription in budding yeast, the helicase Sen1. We identify a phosphorylation at a conserved threonine of the catalytic domain of Sen1 and we provide evidence that phosphorylation at this site reduces the efficiency of Sen1-mediated termination. Interestingly, we find that this phosphorylation impairs termination at an unannotated non-coding gene, thus repressing the expression of a downstream gene encoding the master regulator of Zn homeostasis, Zap1. Consequently, many additional genes exhibit an expression pattern mimicking conditions of Zn excess, where ZAP1 is naturally repressed. Our findings provide a novel paradigm of gene regulatory mechanism relying on the direct modulation of non-coding transcription termination.

Functional profiling of long intergenic non-coding RNAs in fission yeast

Eukaryotic genomes express numerous long intergenic non-coding RNAs (lincRNAs) that do not overlap any coding genes. Some lincRNAs function in various aspects of gene regulation, but it is not clear in general to what extent lincRNAs contribute to the information flow from genotype to phenotype. To explore this question, we systematically analysed cellular roles of lincRNAs in Schizosaccharomyces pombe. Using seamless CRISPR/Cas9-based genome editing, we deleted 141 lincRNA genes to broadly phenotype these mutants, together with 238 diverse coding-gene mutants for functional context. We applied high-throughput colony-based assays to determine mutant growth and viability in benign conditions and in response to 145 different nutrient, drug, and stress conditions. These analyses uncovered phenotypes for 47.5% of the lincRNAs and 96% of the protein-coding genes. For 110 lincRNA mutants, we also performed high-throughput microscopy and flow cytometry assays, linking 37% of these lincRNAs with cell-size and/or cell-cycle control. With all assays combined, we detected phenotypes for 84 (59.6%) of all lincRNA deletion mutants tested. For complementary functional inference, we analysed colony growth of strains ectopically overexpressing 113 lincRNA genes under 47 different conditions. Of these overexpression strains, 102 (90.3%) showed altered growth under certain conditions. Clustering analyses provided further functional clues and relationships for some of the lincRNAs. These rich phenomics datasets associate lincRNA mutants with hundreds of phenotypes, indicating that most of the lincRNAs analysed exert cellular functions in specific environmental or physiological contexts. This study provides groundwork to further dissect the roles of these lincRNAs in the relevant conditions.

Uncovering de novo gene birth in yeast using deep transcriptomics

De novo gene origination has been recently established as an important mechanism for the formation of new genes. In organisms with a large genome, intergenic and intronic regions provide plenty of raw material for new transcriptional events to occur, but little is know about how de novo transcripts originate in more densely-packed genomes. Here, we identify 213 de novo originated transcripts in Saccharomyces cerevisiae using deep transcriptomics and genomic synteny information from multiple yeast species grown in two different conditions. We find that about half of the de novo transcripts are expressed from regions which already harbor other genes in the opposite orientation; these transcripts show similar expression changes in response to stress as their overlapping counterparts, and some appear to translate small proteins. Thus, a large fraction of de novo genes in yeast are likely to co-evolve with already existing genes.

Epigenomic Reprogramming as a Driver of Malignant Glioma

Malignant gliomas are central nervous system tumors and remain among the most treatment-resistant cancers. Exome sequencing has revealed significant heterogeneity and important insights into the molecular pathogenesis of gliomas. Mutations in chromatin modifiers-proteins that shape the epigenomic landscape through remodeling and regulation of post-translational modifications on chromatin-are very frequent and often define specific glioma subtypes. This suggests that epigenomic reprogramming may be a fundamental driver of glioma. Here, we describe the key chromatin regulatory pathways disrupted in gliomas, delineating their physiological function and our current understanding of how their dysregulation may contribute to gliomagenesis.

Advances in RNAi-Assisted Strain Engineering in Saccharomyces cerevisiae

Saccharomyces cerevisiae is a widely used eukaryotic model and microbial cell factory. RNA interference (RNAi) is a conserved regulatory mechanism among eukaryotes but absent from S. cerevisiae. Recent reconstitution of RNAi machinery in S. cerevisiae enables the use of this powerful tool for strain engineering. Here we first discuss the introduction of heterologous RNAi pathways in S. cerevisiae, and the design of various expression cassettes of RNAi precursor reagents for tunable, dynamic, and genome-wide regulation. We then summarize notable examples of RNAi-assisted functional genomics and metabolic engineering studies in S. cerevisiae. We conclude with the future challenges and opportunities of RNAi-based approaches, as well as the potential of other regulatory RNAs in advancing yeast engineering.

Unique and Shared Roles for Histone H3K36 Methylation States in Transcription Regulation Functions

Set2 co-transcriptionally methylates lysine 36 of histone H3 (H3K36), producing mono-, di-, and trimethylation (H3K36me1/2/3). These modifications recruit or repel chromatin effector proteins important for transcriptional fidelity, mRNA splicing, and DNA repair. However, it was not known whether the different methylation states of H3K36 have distinct biological functions. Here, we use engineered forms of Set2 that produce different lysine methylation states to identify unique and shared functions for H3K36 modifications. Although H3K36me1/2 and H3K36me3 are functionally redundant in many SET2 deletion phenotypes, we found that H3K36me3 has a unique function related to Bur1 kinase activity and FACT (facilitates chromatin transcription) complex function. Further, during nutrient stress, either H3K36me1/2 or H3K36me3 represses high levels of histone acetylation and cryptic transcription that arises from within genes. Our findings uncover the potential for the regulation of diverse chromatin functions by different H3K36 methylation states.

Non-coding RNAs as cell wall regulators in Saccharomyces cerevisiae

The cell wall of Saccharomyces cerevisiae is an extracellular organelle crucial for preserving its cellular integrity and detecting environmental cues. The cell wall is composed of mannoproteins attached to a polysaccharide network and is continuously remodelled as cells undergo cell division, mating, gametogenesis or adapt to stressors. This makes yeast an excellent model to study the regulation of genes important for cell wall formation and maintenance. Given that certain yeast strains are pathogenic, a better understanding of their life cycle is of clinical relevance. This is why transcriptional regulatory mechanisms governing genes involved in cell wall biogenesis or maintenance have been the focus of numerous studies. However, little is known about the roles of long non-coding RNAs (lncRNAs), a class of transcripts that are thought to possess little or no protein coding potential, in controlling the expression of cell wall-related genes. This review outlines currently known mechanisms of lncRNA-mediated regulation of gene expression in S. cerevisiae and describes examples of lncRNA-regulated genes encoding cell wall proteins. We suggest that the association of currently annotated lncRNAs with the coding sequences and/or promoters of cell wall-related genes highlights a potential role for lncRNAs as important regulators of the yeast cell wall structure.

Histone H2B Ubiquitylation Regulates Histone Gene Expression by Suppressing Antisense Transcription in Fission Yeast

Histone H2B monoubiquitylation (H2Bub1) is tightly linked to RNA polymerase II transcription elongation, and is also directly implicated in DNA replication and repair. Loss of H2Bub1 is associated with defects in cell cycle progression, but how these are related to its various functions, and the underlying mechanisms involved, is not understood. Here we describe a role for H2Bub1 in the regulation of replication-dependent histone genes in the fission yeast Schizosaccharomyces pombe H2Bub1 activates histone genes indirectly by suppressing antisense transcription of ams2+ -a gene encoding a GATA-type transcription factor that activates histone genes and is required for assembly of centromeric chromatin. Mutants lacking the ubiquitylation site in H2B or the H2B-specific E3 ubiquitin ligase Brl2 had elevated levels of ams2+ antisense transcripts and reduced Ams2 protein levels. These defects were reversed upon inhibition of Cdk9-an ortholog of the kinase component of positive transcription elongation factor b (P-TEFb)-indicating that they likely resulted from aberrant transcription elongation. Reduced Cdk9 activity also partially rescued chromosome segregation phenotypes of H2Bub1 mutants. In a genome-wide analysis, loss of H2Bub1 led to increased antisense transcripts at over 500 protein-coding genes in H2Bub1 mutants; for a subset of these, including several genes involved in chromosome segregation and chromatin assembly, antisense derepression was Cdk9-dependent. Our results highlight antisense suppression as a key feature of cell cycle-dependent gene regulation by H2Bub1, and suggest that aberrant transcription elongation may underlie the effects of H2Bub1 loss on cell cycle progression.

Antisense transcriptional interference mediates condition-specific gene repression in budding yeast

Pervasive transcription generates many unstable non-coding transcripts in budding yeast. The transcription of such noncoding RNAs, in particular antisense RNAs (asRNAs), has been shown in a few examples to repress the expression of the associated mRNAs. Yet, such mechanism is not known to commonly contribute to the regulation of a given class of genes. Using a mutant context that stabilized pervasive transcripts, we observed that the least expressed mRNAs during the exponential phase were associated with high levels of asRNAs. These asRNAs also overlapped their corresponding gene promoters with a much higher frequency than average. Interrupting antisense transcription of a subset of genes corresponding to quiescence-enriched mRNAs restored their expression. The underlying mechanism acts in cis and involves several chromatin modifiers. Our results convey that transcription interference represses up to 30% of the 590 least expressed genes, which includes 163 genes with quiescence-enriched mRNAs. We also found that pervasive transcripts constitute a higher fraction of the transcriptome in quiescence relative to the exponential phase, consistent with gene expression itself playing an important role to suppress pervasive transcription. Accordingly, the HIS1 asRNA, normally only present in quiescence, is expressed in exponential phase upon HIS1 mRNA transcription interruption.

Set2 methyltransferase facilitates cell cycle progression by maintaining transcriptional fidelity

Methylation of histone H3 lysine 36 (H3K36me) by yeast Set2 is critical for the maintenance of chromatin structure and transcriptional fidelity. However, we do not know the full range of Set2/H3K36me functions or the scope of mechanisms that regulate Set2-dependent H3K36 methylation. Here, we show that the APC/C^CDC20 complex regulates Set2 protein abundance during the cell cycle. Significantly, absence of Set2-mediated H3K36me causes a loss of cell cycle control and pronounced defects in the transcriptional fidelity of cell cycle regulatory genes, a class of genes that are generally long, hence highly dependent on Set2/H3K36me for their transcriptional fidelity. Because APC/C also controls human SETD2, and SETD2 likewise regulates cell cycle progression, our data imply an evolutionarily conserved cell cycle function for Set2/SETD2 that may explain why recurrent mutations of SETD2 contribute to human disease.

Characterizing Transcriptional Interference between Converging Genes in Bacteria

Antisense transcription is common in naturally occurring genomes and is increasingly being used in synthetic genetic circuitry as a tool for gene expression control. Mutual influence on the expression of convergent genes can be mediated by antisense RNA effects and by transcriptional interference (TI). We aimed to quantitatively characterize long-range TI between convergent genes with untranslated intergenic spacers of increasing length. After controlling for antisense RNA-mediated effects, which contributed about half of the observed total expression inhibition, the TI effect was modeled. To achieve model convergence, RNA polymerase processivity and collision resistance were assumed to be modulated by ribosome trailing. The spontaneous transcription termination rate in regions of untranslated DNA was experimentally determined. Our modeling suggests that an elongating RNA polymerase with a trailing ribosome is about 13 times more likely to resume transcription than an opposing RNA polymerase without a trailing ribosome, upon head-on collision of the two.

Sensitive high-throughput single-cell RNA-seq reveals within-clonal transcript correlations in yeast populations

Single-cell RNA-seq has revealed extensive cellular heterogeneity within many organisms, but few methods have been developed for microbial clonal populations. The yeast genome displays unusually dense transcript spacing, with interleaved and overlapping transcription from both strands, resulting in a minuscule but complex pool of RNA protected by a resilient cell wall. Here, we have developed a sensitive, scalable, and inexpensive yeast single-cell RNA-seq (yscRNA-seq) method that digitally counts transcript start sites in a strand- and isoform-specific manner. YscRNA-seq detects the expression of low-abundant, non-coding RNAs, and at least half of the protein-coding genome in each cell. Within clonal cells, we observed a negative correlation for the expression of sense/antisense pairs, while paralogs and divergent transcripts co-express. Combining yscRNA-seq with index sorting, we uncovered a linear relationship between cell size and RNA content. Although we detected an average of ~3.5 molecules/gene, the number of expressed isoforms are restricted at the single-cell level. Remarkably, the expression of metabolic genes is highly variable, while their stochastic expression primes cells for increased fitness towards the corresponding environmental challenge. These findings suggest that functional transcript diversity acts as a mechanism for providing a selective advantage to individual cells within otherwise transcriptionally heterogeneous populations.

Bases of antisense lncRNA-associated regulation of gene expression in fission yeast

Antisense (as)lncRNAs can regulate gene expression but the underlying mechanisms and the different cofactors involved remain unclear. Using Native Elongating Transcript sequencing, here we show that stabilization of antisense Exo2-sensitivite lncRNAs (XUTs) results in the attenuation, at the nascent transcription level, of a subset of highly expressed genes displaying prominent promoter-proximal nucleosome depletion and histone acetylation. Mechanistic investigations on the catalase gene ctt1 revealed that its induction following oxidative stress is impaired in Exo2-deficient cells, correlating with the accumulation of an asXUT. Interestingly, expression of this asXUT was also activated in wild-type cells upon oxidative stress, concomitant to ctt1 induction, indicating a potential attenuation feedback. This attenuation correlates with asXUT abundance, it is transcriptional, characterized by low RNAPII-ser5 phosphorylation, and it requires an histone deacetylase activity and the conserved Set2 histone methyltransferase. Finally, we identified Dicer as another RNA processing factor acting on ctt1 induction, but independently of Exo2. We propose that asXUTs could modulate the expression of their paired-sense genes when it exceeds a critical threshold, using a conserved mechanism independent of RNAi.

A current view on long noncoding RNAs in yeast and filamentous fungi

Long noncoding RNAs (lncRNAs) are crucial players in epigenetic regulation. They were initially discovered in human, yet they emerged as common factors involved in a number of central cellular processes in several eukaryotes. For example, in the past decade, research on lncRNAs in yeast has steadily increased. Several examples of lncRNAs were described in Saccharomyces cerevisiae and Schizosaccharomyces pombe. Also, screenings for lncRNAs in ascomycetes were performed and, just recently, the first full characterization of a lncRNA was performed in the filamentous fungus Trichoderma reesei. In this review, we provide a broad overview about currently known fugal lncRNAs. We make an attempt to categorize them according to their functional context, regulatory strategies or special properties. Moreover, the potential of lncRNAs as a biotechnological tool is discussed.

GDV1 induces sexual commitment of malaria parasites by antagonizing HP1-dependent gene silencing

Malaria is caused by Plasmodium parasites that proliferate in the bloodstream. During each replication cycle, some parasites differentiate into gametocytes, the only forms able to infect the mosquito vector and transmit malaria. Sexual commitment is triggered by activation of AP2-G, the master transcriptional regulator of gametocytogenesis. Heterochromatin protein 1 (HP1)-dependent silencing of ap2-g prevents sexual conversion in proliferating parasites. In this study, we identified Plasmodium falciparum gametocyte development 1 (GDV1) as an upstream activator of sexual commitment. We found that GDV1 targeted heterochromatin and triggered HP1 eviction, thus derepressing ap2-g Expression of GDV1 was responsive to environmental triggers of sexual conversion and controlled via a gdv1 antisense RNA. Hence, GDV1 appears to act as an effector protein that induces sexual differentiation by antagonizing HP1-dependent gene silencing.

H3K36 Methylation Regulates Nutrient Stress Response in Saccharomyces cerevisiae by Enforcing Transcriptional Fidelity

Set2-mediated histone methylation at H3K36 regulates diverse activities, including DNA repair, mRNA splicing, and suppression of inappropriate (cryptic) transcription. Although failure of Set2 to suppress cryptic transcription has been linked to decreased lifespan, the extent to which cryptic transcription influences other cellular functions is poorly understood. Here, we uncover a role for H3K36 methylation in the regulation of the nutrient stress response pathway. We found that the transcriptional response to nutrient stress was dysregulated in SET2-deleted (set2Δ) cells and was correlated with genome-wide bi-directional cryptic transcription that originated from within gene bodies. Antisense transcripts arising from these cryptic events extended into the promoters of the genes from which they arose and were associated with decreased sense transcription under nutrient stress conditions. These results suggest that Set2-enforced transcriptional fidelity is critical to the proper regulation of inducible and highly regulated transcription programs.

Large-scale profiling of noncoding RNA function in yeast

Noncoding RNAs (ncRNAs) are emerging as key regulators of cellular function. We have exploited the recently developed barcoded ncRNA gene deletion strain collections in the yeast Saccharomyces cerevisiae to investigate the numerous ncRNAs in yeast with no known function. The ncRNA deletion collection contains deletions of tRNAs, snoRNAs, snRNAs, stable unannotated transcripts (SUTs), cryptic unstable transcripts (CUTs) and other annotated ncRNAs encompassing 532 different individual ncRNA deletions. We have profiled the fitness of the diploid heterozygous ncRNA deletion strain collection in six conditions using batch and continuous liquid culture, as well as the haploid ncRNA deletion strain collections arrayed individually onto solid rich media. These analyses revealed many novel environmental-specific haplo-insufficient and haplo-proficient phenotypes providing key information on the importance of each specific ncRNA in every condition. Co-fitness analysis using fitness data from the heterozygous ncRNA deletion strain collection identified two ncRNA groups required for growth during heat stress and nutrient deprivation. The extensive fitness data for each ncRNA deletion strain has been compiled into an easy to navigate database called Yeast ncRNA Analysis (YNCA). By expanding the original ncRNA deletion strain collection we identified four novel essential ncRNAs; SUT527, SUT075, SUT367 and SUT259/691. We defined the effects of each new essential ncRNA on adjacent gene expression in the heterozygote background identifying both repression and induction of nearby genes. Additionally, we discovered a function for SUT527 in the expression, 3' end formation and localization of SEC4, an essential protein coding mRNA. Finally, using plasmid complementation we rescued the SUT075 lethal phenotype revealing that this ncRNA acts in trans. Overall, our findings provide important new insights into the function of ncRNAs.

Antisense transcription-dependent chromatin signature modulates sense transcript dynamics

Antisense transcription is widespread in genomes. Despite large differences in gene size and architecture, we find that yeast and human genes share a unique, antisense transcription-associated chromatin signature. We asked whether this signature is related to a biological function for antisense transcription. Using quantitative RNA-FISH, we observed changes in sense transcript distributions in nuclei and cytoplasm as antisense transcript levels were altered. To determine the mechanistic differences underlying these distributions, we developed a mathematical framework describing transcription from initiation to transcript degradation. At GAL1, high levels of antisense transcription alter sense transcription dynamics, reducing rates of transcript production and processing, while increasing transcript stability. This relationship with transcript stability is also observed as a genome-wide association. Establishing the antisense transcription-associated chromatin signature through disruption of the Set3C histone deacetylase activity is sufficient to similarly change these rates even in the absence of antisense transcription. Thus, antisense transcription alters sense transcription dynamics in a chromatin-dependent manner.

A cautionary tale of sense-antisense gene pairs: independent regulation despite inverse correlation of expression

Long non-coding RNAs (lncRNAs) have been proven to play important roles in diverse cellular processes including the DNA damage response. Nearly 40% of annotated lncRNAs are transcribed in antisense direction to other genes and have often been implicated in their regulation via transcript- or transcription-dependent mechanisms. However, it remains unclear whether inverse correlation of gene expression would generally point toward a regulatory interaction between the genes. Here, we profiled lncRNA and mRNA expression in lung and liver cancer cells after exposure to DNA damage. Our analysis revealed two pairs of mRNA-lncRNA sense-antisense transcripts being inversely expressed upon DNA damage. The lncRNA NOP14-AS1 was strongly upregulated upon DNA damage, while the mRNA for NOP14 was downregulated, both in a p53-dependent manner. For another pair, the lncRNA LIPE-AS1 was downregulated, while its antisense mRNA CEACAM1 was upregulated. To test whether as expected the antisense genes would regulate each other resulting in this highly significant inverse correlation, we employed antisense oligonucleotides and RNAi to study transcript-dependent effects as well as dCas9-based transcriptional modulation by CRISPRi/CRISPRa for transcription-dependent effects. Surprisingly, despite the strong stimulus-dependent inverse correlation, our data indicate that neither transcript- nor transcription-dependent mechanisms explain the inverse regulation of NOP14-AS1:NOP14 or LIPE-AS1:CEACAM1 expression. Hence, sense-antisense pairs whose expression is strongly—positively or negatively—correlated can be nonetheless regulated independently. This highlights the requirement of individual experimental studies for each antisense pair and prohibits drawing conclusions on regulatory mechanisms from expression correlations.

Identification and characterization of L1-specific endo-siRNAs essential for early embryonic development in pig

Small noncoding RNAs (sncRNAs) play important roles in RNA interference (RNAi). In addition to microRNA (miRNA) and Piwi-interacting RNA (piRNA), one key member of sncRNAs group is endogenous small interfering RNA (endo-siRNA). Some studies do show the role of endo-siRNAs in Dicer and/or Ago mutants, however, the biological functions of specific endo-siRNAs remains mostly unanswered. In the study, we have performed a comparative analysis of endo-siRNAs present in porcine sperms, oocytes and zygotes, identified by deep sequencing and bioinformatics analysis. Further, we observe a large amount of endo-siRNAs specific binding on ORF2 and 3′ UTR of porcine L1 (L1-siRNAs). And, 9 L1-siRNAs generated from a dsRNA formed between L1 transcript and a newly identified an antisense noncoding RNA was characterized. We show the L1-siRNAs regulate early embryonic development by inhibiting the activity of L1 retrotransposition. This work can contribute to understanding the functional role of abundant endo-siRNAs in embryonic development.

Upregulation of SPS100 gene expression by an antisense RNA via a switch of mRNA isoforms with different stabilities

Pervasive transcription of genomes generates multiple classes of non-coding RNAs. One of these classes are stable long non-coding RNAs which overlap coding genes in antisense direction (asRNAs). The function of such asRNAs is not fully understood but several cases of antisense-dependent gene expression regulation affecting the overlapping genes have been demonstrated. Using high-throughput yeast genetics and a limited set of four growth conditions we previously reported a regulatory function for ∼25% of asRNAs, most of which repress the expression of the sense gene. To further explore the roles of asRNAs we tested more conditions and identified 15 conditionally antisense-regulated genes, 6 of which exhibited antisense-dependent enhancement of gene expression. We focused on the sporulation-specific gene SPS100, which becomes upregulated upon entry into starvation or sporulation as a function of the antisense transcript SUT169. We demonstrate that the antisense effect is mediated by its 3' intergenic region (3'-IGR) and that this regulation can be transferred to other genes. Genetic analysis revealed that SUT169 functions by changing the relative expression of SPS100 mRNA isoforms from a short and unstable transcript to a long and stable species. These results suggest a novel mechanism of antisense-dependent gene regulation via mRNA isoform switching.

Deep learning of the regulatory grammar of yeast 5' untranslated regions from 500,000 random sequences

Our ability to predict protein expression from DNA sequence alone remains poor, reflecting our limited understanding of cis-regulatory grammar and hampering the design of engineered genes for synthetic biology applications. Here, we generate a model that predicts the protein expression of the 5′ untranslated region (UTR) of mRNAs in the yeast Saccharomyces cerevisiae. We constructed a library of half a million 50-nucleotide-long random 5′ UTRs and assayed their activity in a massively parallel growth selection experiment. The resulting data allow us to quantify the impact on protein expression of Kozak sequence composition, upstream open reading frames (uORFs), and secondary structure. We trained a convolutional neural network (CNN) on the random library and showed that it performs well at predicting the protein expression of both a held-out set of the random 5′ UTRs as well as native S. cerevisiae 5′ UTRs. The model additionally was used to computationally evolve highly active 5′ UTRs. We confirmed experimentally that the great majority of the evolved sequences led to higher protein expression rates than the starting sequences, demonstrating the predictive power of this model.

Stop that Noise and Turn Up the Antisense Transcription

Many genes are not only transcribed in the sense direction but also yield antisense transcripts. In this issue of Cell Reports, Huber et al. (2016) report that some of these transcripts may serve to suppress sense transcription and noise.

A resource for functional profiling of noncoding RNA in the yeast Saccharomyces cerevisiae

Eukaryotic genomes are extensively transcribed, generating many different RNAs with no known function. We have constructed 1502 molecular barcoded ncRNA gene deletion strains encompassing 443 ncRNAs in the yeast Saccharomyces cerevisiae as tools for ncRNA functional analysis. This resource includes deletions of small nuclear RNAs (snRNAs), transfer RNAs (tRNAs), small nucleolar RNAs (snoRNAs), and other annotated ncRNAs as well as the more recently identified stable unannotated transcripts (SUTs) and cryptic unstable transcripts (CUTs) whose functions are largely unknown. Specifically, deletions have been constructed for ncRNAs found in the intergenic regions, not overlapping genes or their promoters (i.e., at least 200 bp minimum distance from the closest gene start codon). The deletion strains carry molecular barcodes designed to be complementary with the protein gene deletion collection enabling parallel analysis experiments. These strains will be useful for the numerous genomic and molecular techniques that utilize deletion strains, including genome-wide phenotypic screens under different growth conditions, pooled chemogenomic screens with drugs or chemicals, synthetic genetic array analysis to uncover novel genetic interactions, and synthetic dosage lethality screens to analyze gene dosage. Overall, we created a valuable resource for the RNA community and for future ncRNA research.

Mot1, Ino80C, and NC2 Function Coordinately to Regulate Pervasive Transcription in Yeast and Mammals

Pervasive transcription initiates from cryptic promoters and is observed in eukaryotes ranging from yeast to mammals. The Set2-Rpd3 regulatory system prevents cryptic promoter function within expressed genes. However, conserved systems that control pervasive transcription within intergenic regions have not been well established. Here we show that Mot1, Ino80 chromatin remodeling complex (Ino80C), and NC2 co-localize on chromatin and coordinately suppress pervasive transcription in S. cerevisiae and murine embryonic stem cells (mESCs). In yeast, all three proteins bind subtelomeric heterochromatin through a Sir3-stimulated mechanism and to euchromatin via a TBP-stimulated mechanism. In mESCs, the proteins bind to active and poised TBP-bound promoters along with promoters of polycomb-silenced genes apparently lacking TBP. Depletion of Mot1, Ino80C, or NC2 by anchor away in yeast or RNAi in mESCs leads to near-identical transcriptome phenotypes, with new subtelomeric transcription in yeast, and greatly increased pervasive transcription in both yeast and mESCs.

Shaping the cellular landscape with Set2/SETD2 methylation

Chromatin structure is a major barrier to gene transcription that must be disrupted and re-set during each round of transcription. Central to this process is the Set2/SETD2 methyltransferase that mediates co-transcriptional methylation to histone H3 at lysine 36 (H3K36me). Studies reveal that H3K36me not only prevents inappropriate transcriptional initiation from arising within gene bodies, but that it has other conserved functions that include the repair of damaged DNA and regulation of pre-mRNA splicing. Consistent with the importance of Set2/SETD2 in chromatin biology, mutations of SETD2, or mutations at or near H3K36 in H3.3, have recently been found to underlie cancer development. This review will summarize the latest insights into the functions of Set2/SETD2 in genome regulation and cancer development.

Transcriptional interference by RNA polymerase III affects expression of the Polr3e gene

A mammalian interspersed repeat (MIR) nested in antisense orientation within the first intron of the Polr3e gene, encoding an RNA polymerase (Pol) III subunit, is conserved in mammals and highly occupied by Pol III. Here, Yeganeh et al. show that the MIR affects Polr3e expression through transcriptional interference.

Overlapping gene arrangements can potentially contribute to gene expression regulation. A mammalian interspersed repeat (MIR) nested in antisense orientation within the first intron of the Polr3e gene, encoding an RNA polymerase III (Pol III) subunit, is conserved in mammals and highly occupied by Pol III. Using a fluorescence assay, CRISPR/Cas9-mediated deletion of the MIR in mouse embryonic stem cells, and chromatin immunoprecipitation assays, we show that the MIR affects Polr3e expression through transcriptional interference. Our study reveals a mechanism by which a Pol II gene can be regulated at the transcription elongation level by transcription of an embedded antisense Pol III gene.

Long Noncoding RNAs in the Yeast S. cerevisiae

Long noncoding RNAs have recently been discovered to comprise a sizeable fraction of the RNA World. The scope of their functions, physical organization, and disease relevance remain in the early stages of characterization. Although many thousands of lncRNA transcripts recently have been found to emanate from the expansive DNA between protein-coding genes in animals, there are also hundreds that have been found in simple eukaryotes. Furthermore, lncRNAs have been found in the bacterial and archaeal branches of the tree of life, suggesting they are ubiquitous. In this chapter, we focus primarily on what has been learned so far about lncRNAs from the greatly studied single-celled eukaryote, the yeast Saccharomyces cerevisiae. Most lncRNAs examined in yeast have been implicated in transcriptional regulation of protein-coding genes-often in response to forms of stress-whereas a select few have been ascribed yet other functions. Of those known to be involved in transcriptional regulation of protein-coding genes, the vast majority function in cis. There are also some yeast lncRNAs identified that are not directly involved in regulation of transcription. Examples of these include the telomerase RNA and telomere-encoded transcripts. In addition to its role as a template-encoding telomeric DNA synthesis, telomerase RNA has been shown to function as a flexible scaffold for protein subunits of the RNP holoenzyme. The flexible scaffold model provides a specific mechanistic paradigm that is likely to apply to many other lncRNAs that assemble and orchestrate large RNP complexes, even in humans. Looking to the future, it is clear that considerable fundamental knowledge remains to be obtained about the architecture and functions of lncRNAs. Using genetically tractable unicellular model organisms should facilitate lncRNA characterization. The acquired basic knowledge will ultimately translate to better understanding of the growing list of lncRNAs linked to human maladies.

Functional characterisation of long intergenic non-coding RNAs through genetic interaction profiling in Saccharomyces cerevisiae

BACKGROUND

Transcriptome studies have revealed that many eukaryotic genomes are pervasively transcribed producing numerous long non-coding RNAs (lncRNAs). However, only a few lncRNAs have been ascribed a cellular role thus far, with most regulating the expression of adjacent genes. Even less lncRNAs have been annotated as essential hence implying that the majority may be functionally redundant. Therefore, the function of lncRNAs could be illuminated through systematic analysis of their synthetic genetic interactions (GIs).

RESULTS

Here, we employ synthetic genetic array (SGA) in Saccharomyces cerevisiae to identify GIs between long intergenic non-coding RNAs (lincRNAs) and protein-coding genes. We first validate this approach by demonstrating that the telomerase RNA TLC1 displays a GI network that corresponds to its well-described function in telomere length maintenance. We subsequently performed SGA screens on a set of uncharacterised lincRNAs and uncover their connection to diverse cellular processes. One of these lincRNAs, SUT457, exhibits a GI profile associating it to telomere organisation and we consistently demonstrate that SUT457 is required for telomeric overhang homeostasis through an Exo1-dependent pathway. Furthermore, the GI profile of SUT457 is distinct from that of its neighbouring genes suggesting a function independent to its genomic location. Accordingly, we show that ectopic expression of this lincRNA suppresses telomeric overhang accumulation in sut457Δ cells assigning a trans-acting role for SUT457 in telomere biology.

CONCLUSIONS

Overall, our work proposes that systematic application of this genetic approach could determine the functional significance of individual lncRNAs in yeast and other complex organisms.