Strand-specific RNA sequencing reveals extensive regulated long antisense transcripts that are conserved across yeast species

Efficient production of RNA in Saccharomyces cerevisiae through inducing high level transcription of functional ncRNA-SRG1

RNA (Ribonucleic Acid) is an essential component of organisms and is widely used in the food and pharmaceutical industries. Saccharomyces cerevisiae, recognized as a safe strain, is widely used for RNA production. In this study, the S. cerevisiae W303-1a was used as a starting strain and molecular modifications were made to the functional ncRNA-SRG1 to evaluate the effect on RNA production. At the same time, its transcriptionally associated helper genes (Spt2, Spt6 and Cha4) were overexpressed and the culture medium was supplemented with serine to induce SRG1 transcription, to increase SRG1 transcription levels and investigate its effect on intracellular RNA levels. The results showed that the intracellular RNA content of the recombinant strain W303-1a-SRG1 was 10.27 %, an increase of 11.15 % compared to the starting strain (W303-1a, with an intracellular RNA content of 9.24 %). On this basis, a gene co-overexpression strain-W303-1a-SRG1-Spt6 was constructed. Simultaneously, the addition of 2 % serine strategy was used to increase the transcription level of SRG1 and RNA content of the recombinant strain. The intracellular RNA of the recombinant strain reached 11.41 %, an increase of 23.38 % compared to the starting strain (W303-1a, without serine supplementation). In addition, the growth performance of the strain was assessed by measuring the SRG1 transcription level in the strain and plotting the growth curve. Therefore, we found that improving the transcription level of ncRNA can be used as a new idea to construct S. cerevisiae with high RNA content, which provides a strong help for subsequent research in related fields. This work provides a new strategy for increasing the nucleic acid content of S. cerevisiae.

Adapting cytoskeleton-mitochondria patterning with myocyte differentiation by promyogenic PRR33

Coordinated cytoskeleton-mitochondria organization during myogenesis is crucial for muscle development and function. Our understanding of the underlying regulatory mechanisms remains inadequate. Here, we identified a novel muscle-enriched protein, PRR33, which is upregulated during myogenesis and acts as a promyogenic factor. Depletion of Prr33 in C2C12 represses myoblast differentiation. Genetic deletion of Prr33 in mice reduces myofiber size and decreases muscle strength. The Prr33 mutant mice also exhibit impaired myogenesis and defects in muscle regeneration in response to injury. Interactome and transcriptome analyses reveal that PRR33 regulates cytoskeleton and mitochondrial function. Remarkably, PRR33 interacts with DESMIN, a key regulator of cytoskeleton-mitochondria organization in muscle cells. Abrogation of PRR33 in myocytes substantially abolishes the interaction of DESMIN filaments with mitochondria, leading to abnormal intracellular accumulation of DESMIN and mitochondrial disorganization/dysfunction in myofibers. Together, our findings demonstrate that PRR33 and DESMIN constitute an important regulatory module coordinating mitochondrial organization with muscle differentiation.

Natural antisense transcripts as versatile regulators of gene expression

Long non-coding RNAs (lncRNAs) are emerging as a major class of gene products that have central roles in cell and developmental biology. Natural antisense transcripts (NATs) are an important subset of lncRNAs that are expressed from the opposite strand of protein-coding and non-coding genes and are a genome-wide phenomenon in both eukaryotes and prokaryotes. In eukaryotes, a myriad of NATs participate in regulatory pathways that affect expression of their cognate sense genes. Recent developments in the study of NATs and lncRNAs and large-scale sequencing and bioinformatics projects suggest that whether NATs regulate expression, splicing, stability or translation of the sense transcript is influenced by the pattern and degrees of overlap between the sense-antisense pair. Moreover, epigenetic gene regulatory mechanisms prevail in somatic cells whereas mechanisms dependent on the formation of double-stranded RNA intermediates are prevalent in germ cells. The modulating effects of NATs on sense transcript expression make NATs rational targets for therapeutic interventions.

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.

Exogenous and endogenous dsRNAs perceived by plant Dicer-like 4 protein in the RNAi-depleted cellular context

BACKGROUND

In plants, RNase III Dicer-like proteins (DCLs) act as sensors of dsRNAs and process them into short 21- to 24-nucleotide (nt) (s)RNAs. Plant DCL4 is involved in the biogenesis of either functional endogenous or exogenous (i.e. viral) short interfering (si)RNAs, thus playing crucial antiviral roles.

METHODS

In this study we expressed plant DCL4 in Saccharomyces cerevisiae, an RNAi-depleted organism, in which we could highlight the role of dicing as neither Argonautes nor RNA-dependent RNA polymerase is present. We have therefore tested the DCL4 functionality in processing exogenous dsRNA-like substrates, such as a replicase-assisted viral replicon defective-interfering RNA and RNA hairpin substrates, or endogenous antisense transcripts.

RESULTS

DCL4 was shown to be functional in processing dsRNA-like molecules in vitro and in vivo into 21- and 22-nt sRNAs. Conversely, DCL4 did not efficiently process a replicase-assisted viral replicon in vivo, providing evidence that viral RNAs are not accessible to DCL4 in membranes associated in active replication. Worthy of note, in yeast cells expressing DCL4, 21- and 22-nt sRNAs are associated with endogenous loci.

CONCLUSIONS

We provide new keys to interpret what was studied so far on antiviral DCL4 in the host system. The results all together confirm the role of sense/antisense RNA-based regulation of gene expression, expanding the sense/antisense atlas of S. cerevisiae. The results described herein show that S. cerevisiae can provide insights into the functionality of plant dicers and extend the S. cerevisiae tool to new biotechnological applications.

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.

The evolution, evolvability and engineering of gene regulatory DNA

Mutations in non-coding regulatory DNA sequences can alter gene expression, organismal phenotype and fitness^1-3. Constructing complete fitness landscapes, in which DNA sequences are mapped to fitness, is a long-standing goal in biology, but has remained elusive because it is challenging to generalize reliably to vast sequence spaces^4-6. Here we build sequence-to-expression models that capture fitness landscapes and use them to decipher principles of regulatory evolution. Using millions of randomly sampled promoter DNA sequences and their measured expression levels in the yeast Saccharomyces cerevisiae, we learn deep neural network models that generalize with excellent prediction performance, and enable sequence design for expression engineering. Using our models, we study expression divergence under genetic drift and strong-selection weak-mutation regimes to find that regulatory evolution is rapid and subject to diminishing returns epistasis; that conflicting expression objectives in different environments constrain expression adaptation; and that stabilizing selection on gene expression leads to the moderation of regulatory complexity. We present an approach for using such models to detect signatures of selection on expression from natural variation in regulatory sequences and use it to discover an instance of convergent regulatory evolution. We assess mutational robustness, finding that regulatory mutation effect sizes follow a power law, characterize regulatory evolvability, visualize promoter fitness landscapes, discover evolvability archetypes and illustrate the mutational robustness of natural regulatory sequence populations. Our work provides a general framework for designing regulatory sequences and addressing fundamental questions in regulatory evolution.

Advances in Non-Coding RNA Sequencing

Non-coding RNAs (ncRNAs) comprise a set of abundant and functionally diverse RNA molecules. Since the discovery of the first ncRNA in the 1960s, ncRNAs have been shown to be involved in nearly all steps of the central dogma of molecular biology. In recent years, the pace of discovery of novel ncRNAs and their cellular roles has been greatly accelerated by high-throughput sequencing. Advances in sequencing technology, library preparation protocols as well as computational biology helped to greatly expand our knowledge of which ncRNAs exist throughout the kingdoms of life. Moreover, RNA sequencing revealed crucial roles of many ncRNAs in human health and disease. In this review, we discuss the most recent methodological advancements in the rapidly evolving field of high-throughput sequencing and how it has greatly expanded our understanding of ncRNA biology across a large number of different organisms.

Extracting novel hypotheses and findings from RNA-seq data

Over the past decade, improvements in technology and methods have enabled rapid and relatively inexpensive generation of high-quality RNA-seq datasets. These datasets have been used to characterize gene expression for several yeast species and have provided systems-level insights for basic biology, biotechnology and medicine. Herein, we discuss new techniques that have emerged and existing techniques that enable analysts to extract information from multifactorial yeast RNA-seq datasets. Ultimately, this minireview seeks to inspire readers to query datasets, whether previously published or freshly obtained, with creative and diverse methods to discover and support novel hypotheses.

Identification of long noncoding natural antisense transcripts (lncNATs) correlated with drought stress response in wild rice (Oryza nivara)

Background

Wild rice, including Oryza nivara and Oryza rufipogon, which are considered as the ancestors of Asian cultivated rice (Oryza sativa), possess high genetic diversity and serve as a crucial resource for breeding novel cultivars of cultivated rice. Although rice domestication related traits, such as seed shattering and plant architecture, have been intensively studied at the phenotypic and genomic levels, further investigation is needed to understand the molecular basis of phenotypic differences between cultivated and wild rice. Drought stress is one of the most severe abiotic stresses affecting rice growth and production. Adaptation to drought stress involves a cascade of genes and regulatory factors that form complex networks. O. nivara inhabits swampy areas with a seasonally dry climate, which is an ideal material to discover drought tolerance alleles. Long noncoding natural antisense transcripts (lncNATs), a class of long noncoding RNAs (lncRNAs), regulate the corresponding sense transcripts and play an important role in plant growth and development. However, the contribution of lncNATs to drought stress response in wild rice remains largely unknown.

Results

Here, we conducted strand-specific RNA sequencing (ssRNA-seq) analysis of Nipponbare (O. sativa) and two O. nivara accessions (BJ89 and BJ278) to determine the role of lncNATs in drought stress response in wild rice. A total of 1246 lncRNAs were identified, including 1091 coding–noncoding NAT pairs, of which 50 were expressed only in Nipponbare, and 77 were expressed only in BJ89 and/or BJ278. Of the 1091 coding–noncoding NAT pairs, 240 were differentially expressed between control and drought stress conditions. Among these 240 NAT pairs, 12 were detected only in Nipponbare, and 187 were detected uniquely in O. nivara. Furthermore, 10 of the 240 coding–noncoding NAT pairs were correlated with genes enriched in stress responsive GO terms; among these, nine pairs were uniquely found in O. nivara, and one pair was shared between O. nivara and Nipponbare.

Conclusion

We identified lncNATs associated with drought stress response in cultivated rice and O. nivara. These results will improve our understanding of the function of lncNATs in drought tolerance and accelerate rice breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07754-4.

Natural antisense transcript of Period2, Per2AS, regulates the amplitude of the mouse circadian clock

In mammals, a set of core clock genes form transcription-translation feedback loops to generate circadian oscillations. We and others recently identified a novel transcript at the Period2 (Per2) locus that is transcribed from the antisense strand of Per2 This transcript, Per2AS, is expressed rhythmically and antiphasic to Per2 mRNA, leading to our hypothesis that Per2AS and Per2 mutually inhibit each other's expression and form a double negative feedback loop. By perturbing the expression of Per2AS, we found that Per2AS transcription, but not transcript, represses Per2 However, Per2 does not repress Per2AS, as Per2 knockdown led to a decrease in the Per2AS level, indicating that Per2AS forms a single negative feedback loop with Per2 and maintains the level of Per2 within the oscillatory range. Per2AS also regulates the amplitude of the circadian clock, and this function cannot be solely explained through its interaction with Per2, as Per2 knockdown does not recapitulate the phenotypes of Per2AS perturbation. Overall, our data indicate that Per2AS is an important regulatory molecule in the mammalian circadian clock machinery. Our work also supports the idea that antisense transcripts of core clock genes constitute a common feature of circadian clocks, as they are found in other organisms.

HIV-1 Natural Antisense Transcription and Its Role in Viral Persistence

Natural antisense transcripts (NATs) represent a class of RNA molecules that are transcribed from the opposite strand of a protein-coding gene, and that have the ability to regulate the expression of their cognate protein-coding gene via multiple mechanisms. NATs have been described in many prokaryotic and eukaryotic systems, as well as in the viruses that infect them. The human immunodeficiency virus (HIV-1) is no exception, and produces one or more NAT from a promoter within the 3’ long terminal repeat. HIV-1 antisense transcripts have been the focus of several studies spanning over 30 years. However, a complete appreciation of the role that these transcripts play in the virus lifecycle is still lacking. In this review, we cover the current knowledge about HIV-1 NATs, discuss some of the questions that are still open and identify possible areas of future research.

UNAGI: an automated pipeline for nanopore full-length cDNA sequencing uncovers novel transcripts and isoforms in yeast

Sequencing the entire RNA molecule leads to a better understanding of the transcriptome architecture. SMARTer (Switching Mechanism at 5'-End of RNA Template) is a technology aimed at generating full-length cDNA from low amounts of mRNA for sequencing by short-read sequencers such as those from Illumina. However, short read sequencing such as Illumina technology includes fragmentation that results in bias and information loss. Here, we built a pipeline, UNAGI or UNAnnotated Gene Identifier, to process long reads obtained with nanopore sequencing and compared this pipeline with the standard Illumina pipeline by studying the Saccharomyces cerevisiae transcriptome in full-length cDNA samples generated from two different biological samples: haploid and diploid cells. Additionally, we processed the long reads with another long read tool, FLAIR. Our strand-aware method revealed significant differential gene expression that was masked in Illumina data by antisense transcripts. Our pipeline, UNAGI, outperformed the Illumina pipeline and FLAIR in transcript reconstruction (sensitivity and specificity of 80% and 40% vs. 18% and 34% and 79% and 32%, respectively). Moreover, UNAGI discovered 3877 unannotated transcripts including 1282 intergenic transcripts while the Illumina pipeline discovered only 238 unannotated transcripts. For isoforms profiling, UNAGI also outperformed the Illumina pipeline and FLAIR in terms of sensitivity (91% vs. 82% and 63%, respectively). But the low accuracy of nanopore sequencing led to a closer gap in terms of specificity with Illumina pipeline (70% vs. 63%) and to a huge gap with FLAIR (70% vs 0.02%).

Transcriptome visualization and data availability at the Saccharomyces Genome Database

The Saccharomyces Genome Database (SGD; www.yeastgenome.org) maintains the official annotation of all genes in the Saccharomyces cerevisiae reference genome and aims to elucidate the function of these genes and their products by integrating manually curated experimental data. Technological advances have allowed researchers to profile RNA expression and identify transcripts at high resolution. These data can be configured in web-based genome browser applications for display to the general public. Accordingly, SGD has incorporated published transcript isoform data in our instance of JBrowse, a genome visualization platform. This resource will help clarify S. cerevisiae biological processes by furthering studies of transcriptional regulation, untranslated regions, genome engineering, and expression quantification in S. cerevisiae.

Accurate and Sensitive Quantitation of the Dynamic Heat Shock Proteome Using Tandem Mass Tags

Cells respond to environmental perturbations and insults through modulating protein abundance and function. However, the majority of studies have focused on changes in RNA abundance because quantitative transcriptomics has historically been more facile than quantitative proteomics. Modern Orbitrap mass spectrometers now provide sensitive and deep proteome coverage, allowing direct, global quantification of not only protein abundance but also post-translational modifications (PTMs) that regulate protein activity. We implemented and validated using the well-characterized heat shock response of budding yeast, a tandem mass tagging (TMT), triple-stage mass spectrometry (MS³) strategy to measure global changes in the proteome during the yeast heat shock response over nine time points. We report that basic-pH, ultra-high performance liquid chromatography (UPLC) fractionation of tryptic peptides yields superfractions of minimal redundancy, a crucial requirement for deep coverage and quantification by subsequent LC-MS³. We quantified 2275 proteins across three biological replicates and found that differential expression peaked near 90 min following heat shock (with 868 differentially expressed proteins at 5% false discovery rate). The sensitivity of the approach also allowed us to detect changes in the relative abundance of ubiquitination and phosphorylation PTMs over time. Remarkably, relative quantification of post-translationally modified peptides revealed striking evidence of regulation of the heat shock response by protein PTMs. These data demonstrate that the high precision of TMT-MS³ enables peptide-level quantification of samples, which can reveal important regulation of protein abundance and regulatory PTMs under various experimental conditions.

Rapidly evolving protointrons in Saccharomyces genomes revealed by a hungry spliceosome

Introns are a prevalent feature of eukaryotic genomes, yet their origins and contributions to genome function and evolution remain mysterious. In budding yeast, repression of the highly transcribed intron-containing ribosomal protein genes (RPGs) globally increases splicing of non-RPG transcripts through reduced competition for the spliceosome. We show that under these “hungry spliceosome” conditions, splicing occurs at more than 150 previously unannotated locations we call protointrons that do not overlap known introns. Protointrons use a less constrained set of splice sites and branchpoints than standard introns, including in one case AT-AC in place of GT-AG. Protointrons are not conserved in all closely related species, suggesting that most are not under positive selection and are fated to disappear. Some are found in non-coding RNAs (e. g. CUTs and SUTs), where they may contribute to the creation of new genes. Others are found across boundaries between noncoding and coding sequences, or within coding sequences, where they offer pathways to the creation of new protein variants, or new regulatory controls for existing genes. We define protointrons as (1) nonconserved intron-like sequences that are (2) infrequently spliced, and importantly (3) are not currently understood to contribute to gene expression or regulation in the way that standard introns function. A very few protointrons in S. cerevisiae challenge this classification by their increased splicing frequency and potential function, consistent with the proposed evolutionary process of “intronization”, whereby new standard introns are created. This snapshot of intron evolution highlights the important role of the spliceosome in the expansion of transcribed genomic sequence space, providing a pathway for the rare events that may lead to the birth of new eukaryotic genes and the refinement of existing gene function.

The protein coding information in eukaryotic genes is broken by intervening sequences called introns that are removed from RNA during transcription by a large protein-RNA complex called the spliceosome. Where introns come from and how the spliceosome contributes to genome evolution are open questions. In this study, we find more than 150 new places in the yeast genome that are recognized by the spliceosome and spliced out as introns. Since they appear to have arisen very recently in evolution by sequence drift and do not appear to contribute to gene expression or its regulation, we call these protointrons. Protointrons are found in both protein-coding and non-coding RNAs and are not efficiently removed by the splicing machinery. Although most protointrons are not conserved and will likely disappear as evolution proceeds, a few are spliced more efficiently, and are located where they might begin to play functional roles in gene expression, as predicted by the proposed process of intronization. The challenge now is to understand how spontaneously appearing splicing events like protointrons might contribute to the creation of new genes, new genetic controls, and new protein isoforms as genomes evolve.

Rapidly evolving protointrons in Saccharomyces genomes revealed by a hungry spliceosome

Introns are a prevalent feature of eukaryotic genomes, yet their origins and contributions to genome function and evolution remain mysterious. In budding yeast, repression of the highly transcribed intron-containing ribosomal protein genes (RPGs) globally increases splicing of non-RPG transcripts through reduced competition for the spliceosome. We show that under these "hungry spliceosome" conditions, splicing occurs at more than 150 previously unannotated locations we call protointrons that do not overlap known introns. Protointrons use a less constrained set of splice sites and branchpoints than standard introns, including in one case AT-AC in place of GT-AG. Protointrons are not conserved in all closely related species, suggesting that most are not under positive selection and are fated to disappear. Some are found in non-coding RNAs (e. g. CUTs and SUTs), where they may contribute to the creation of new genes. Others are found across boundaries between noncoding and coding sequences, or within coding sequences, where they offer pathways to the creation of new protein variants, or new regulatory controls for existing genes. We define protointrons as (1) nonconserved intron-like sequences that are (2) infrequently spliced, and importantly (3) are not currently understood to contribute to gene expression or regulation in the way that standard introns function. A very few protointrons in S. cerevisiae challenge this classification by their increased splicing frequency and potential function, consistent with the proposed evolutionary process of "intronization", whereby new standard introns are created. This snapshot of intron evolution highlights the important role of the spliceosome in the expansion of transcribed genomic sequence space, providing a pathway for the rare events that may lead to the birth of new eukaryotic genes and the refinement of existing gene function.

Robustness by intrinsically disordered C-termini and translational readthrough

During protein synthesis genetic instructions are passed from DNA via mRNA to the ribosome to assemble a protein chain. Occasionally, stop codons in the mRNA are bypassed and translation continues into the untranslated region (3′-UTR). This process, called translational readthrough (TR), yields a protein chain that becomes longer than would be predicted from the DNA sequence alone. Protein sequences vary in propensity for translational errors, which may yield evolutionary constraints by limiting evolutionary paths. Here we investigated TR in Saccharomyces cerevisiae by analysing ribosome profiling data. We clustered proteins as either prone or non-prone to TR, and conducted comparative analyses. We find that a relatively high frequency (5%) of genes undergo TR, including ribosomal subunit proteins. Our main finding is that proteins undergoing TR are highly expressed and have a higher proportion of intrinsically disordered C-termini. We suggest that highly expressed proteins may compensate for the deleterious effects of TR by having intrinsically disordered C-termini, which may provide conformational flexibility but without distorting native function. Moreover, we discuss whether minimizing deleterious effects of TR is also enabling exploration of the phenotypic landscape of protein isoforms.

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.

The Paf1 Complex Broadly Impacts the Transcriptome of Saccharomyces cerevisiae

The Polymerase Associated Factor 1 complex (Paf1C) is a multifunctional regulator of eukaryotic gene expression important for the coordination of transcription with chromatin modification and post-transcriptional processes. In this study, we investigated the extent to which the functions of Paf1C combine to regulate the Saccharomyces cerevisiae transcriptome. While previous studies focused on the roles of Paf1C in controlling mRNA levels, here, we took advantage of a genetic background that enriches for unstable transcripts, and demonstrate that deletion of PAF1 affects all classes of Pol II transcripts including multiple classes of noncoding RNAs (ncRNAs). By conducting a de novo differential expression analysis independent of gene annotations, we found that Paf1 positively and negatively regulates antisense transcription at multiple loci. Comparisons with nascent transcript data revealed that many, but not all, changes in RNA levels detected by our analysis are due to changes in transcription instead of post-transcriptional events. To investigate the mechanisms by which Paf1 regulates protein-coding genes, we focused on genes involved in iron and phosphate homeostasis, which were differentially affected by PAF1 deletion. Our results indicate that Paf1 stimulates phosphate gene expression through a mechanism that is independent of any individual Paf1C-dependent histone modification. In contrast, the inhibition of iron gene expression by Paf1 correlates with a defect in H3 K36 trimethylation. Finally, we showed that one iron regulon gene, FET4, is coordinately controlled by Paf1 and transcription of upstream noncoding DNA. Together, these data identify roles for Paf1C in controlling both coding and noncoding regions of the yeast genome.

A new experimental platform facilitates assessment of the transcriptional and chromatin landscapes of aging yeast

Replicative aging of Saccharomyces cerevisiae is an established model system for eukaryotic cellular aging. A limitation in yeast lifespan studies has been the difficulty of separating old cells from young cells in large quantities. We engineered a new platform, the Miniature-chemostat Aging Device (MAD), that enables purification of aged cells at sufficient quantities for genomic and biochemical characterization of aging yeast populations. Using MAD, we measured DNA accessibility and gene expression changes in aging cells. Our data highlight an intimate connection between aging, growth rate, and stress. Stress-independent genes that change with age are highly enriched for targets of the signal recognition particle (SRP). Combining MAD with an improved ATAC-seq method, we find that increasing proteasome activity reduces rDNA instability usually observed in aging cells and, contrary to published findings, provide evidence that global nucleosome occupancy does not change significantly with age.

Catabolic repression in early-diverging anaerobic fungi is partially mediated by natural antisense transcripts

Early-diverging anaerobic fungi (order: Neocallimastigomycota), lignocelluolytic chytrid-like fungi central to fiber degradation in the digestive tracts of large herbivores, are attractive sources of cellulases and hemicellulases for biotechnology. Enzyme expression is tightly regulated and coordinated through mechanisms that remain unelucidated to optimize hydrolytic efficiency. Our analysis of anaerobic fungal transcriptomes reveals hundreds of cis-natural antisense transcripts (cis-NATs), which we hypothesize play an integral role in this regulation. Through integrated genomic and transcriptomic sequencing on a range of catabolic substrates, we validate these NATs in three species (Anaeromyces robustus, Neocallimasix californiae, and Piromyces finnis), and analyze their expression patterns and prevalence to gain insight into their function. NAT function was diverse and conserved across the three fungal genomes studied, with 10% of all metabolic process NATs associated with lignocellulose hydrolysis. Despite these similarities, however, only eleven gene targets were conserved orthologs. Several NATs were dynamically regulated by lignocellulosic substrates while their gene targets were unregulated. This observation is consistent with a hypothesized, but untested, regulatory mechanism where selected genes are exclusively regulated at the transcriptional/post-transcriptional level by NATs. However, only genes with high NAT relative expression levels displayed this phenomenon, suggesting a selection mechanism that favors larger dynamic ranges for more precise control of gene expression. In addition to this mode, we observed two other possible regulatory fates: canonical transcriptional regulation with no NAT response, and positive co-regulation of target mRNA and cognate NAT, which we hypothesize is a fine-tuning strategy to locally negate control outputs from global regulators. Our work reveals the complex contributions of antisense RNA to the catabolic response in anaerobic fungi, highlighting its importance in understanding lignocellulolytic activity for bioenergy applications. More importantly, the relative expression of NAT to target may form a critical determinant of transcriptional vs post-transcriptional (NAT) control of gene expression in primitive anaerobic fungi.

Widespread antisense transcription of Populus genome under drought

Antisense transcription is widespread in many genomes and plays important regulatory roles in gene expression. The objective of our study was to investigate the extent and functional relevance of antisense transcription in forest trees. We employed Populus, a model tree species, to probe the antisense transcriptional response of tree genome under drought, through stranded RNA-seq analysis. We detected nearly 48% of annotated Populus gene loci with antisense transcripts and 44% of them with co-transcription from both DNA strands. Global distribution of reads pattern across annotated gene regions uncovered that antisense transcription was enriched in untranslated regions while sense reads were predominantly mapped in coding exons. We further detected 1185 drought-responsive sense and antisense gene loci and identified a strong positive correlation between the expression of antisense and sense transcripts. Additionally, we assessed the antisense expression in introns and found a strong correlation between intronic expression and exonic expression, confirming antisense transcription of introns contributes to transcriptional activity of Populus genome under drought. Finally, we functionally characterized drought-responsive sense-antisense transcript pairs through gene ontology analysis and discovered that functional groups including transcription factors and histones were concordantly regulated at both sense and antisense transcriptional level. Overall, our study demonstrated the extensive occurrence of antisense transcripts of Populus genes under drought and provided insights into genome structure, regulation pattern and functional significance of drought-responsive antisense genes in forest trees. Datasets generated in this study serve as a foundation for future genetic analysis to improve our understanding of gene regulation by antisense transcription.

DMS-Seq for In Vivo Genome-wide Mapping of Protein-DNA Interactions and Nucleosome Centers

Protein-DNA interactions provide the basis for chromatin structure and gene regulation. Comprehensive identification of protein-occupied sites is thus vital to an in-depth understanding of genome function. Dimethyl sulfate (DMS) is a chemical probe that has long been used to detect footprints of DNA-bound proteins in vitro and in vivo. Here, we describe a genomic footprinting method, dimethyl sulfate sequencing (DMS-seq), which exploits the cell-permeable nature of DMS to obviate the need for nuclear isolation. This feature makes DMS-seq simple in practice and removes the potential risk of protein re-localization during nuclear isolation. DMS-seq successfully detects transcription factors bound to cis-regulatory elements and non-canonical chromatin particles in nucleosome-free regions. Furthermore, an unexpected preference of DMS confers on DMS-seq a unique potential to directly detect nucleosome centers without using genetic manipulation. We expect that DMS-seq will serve as a characteristic method for genome-wide interrogation of in vivo protein-DNA interactions.

Long Noncoding RNAs in Yeast Cells and Differentiated Subpopulations of Yeast Colonies and Biofilms

We summarize current knowledge regarding regulatory functions of long noncoding RNAs (lncRNAs) in yeast, with emphasis on lncRNAs identified recently in yeast colonies and biofilms. Potential regulatory functions of these lncRNAs in differentiated cells of domesticated colonies adapted to plentiful conditions versus yeast colony biofilms are discussed. We show that specific cell types differ in their complements of lncRNA, that this complement changes over time in differentiating upper cells, and that these lncRNAs target diverse functional categories of genes in different cell subpopulations and specific colony types.

Transcriptome Remodeling of Differentiated Cells during Chronological Ageing of Yeast Colonies: New Insights into Metabolic Differentiation

We present the spatiotemporal metabolic differentiation of yeast cell subpopulations from upper, lower, and margin regions of colonies of different ages, based on comprehensive transcriptomic analysis. Furthermore, the analysis was extended to include smaller cell subpopulations identified previously by microscopy within fully differentiated U and L cells of aged colonies. New data from RNA-seq provides both spatial and temporal information on cell metabolic reprogramming during colony ageing and shows that cells at marginal positions are similar to upper cells, but both these cell types are metabolically distinct from cells localized to lower colony regions. As colonies age, dramatic metabolic reprogramming occurs in cells of upper regions, while changes in margin and lower cells are less prominent. Interestingly, whereas clear expression differences were identified between two L cell subpopulations, U cells (which adopt metabolic profiles, similar to those of tumor cells) form a more homogeneous cell population. The data identified crucial metabolic reprogramming events that arise de novo during colony ageing and are linked to U and L cell colony differentiation and support a role for mitochondria in this differentiation process.

Quantitating translational control: mRNA abundance-dependent and independent contributions and the mRNA sequences that specify them

Translation rate per mRNA molecule correlates positively with mRNA abundance. As a result, protein levels do not scale linearly with mRNA levels, but instead scale with the abundance of mRNA raised to the power of an ‘amplification exponent’. Here we show that to quantitate translational control, the translation rate must be decomposed into two components. One, TR_mD, depends on the mRNA level and defines the amplification exponent. The other, TR_mIND, is independent of mRNA amount and impacts the correlation coefficient between protein and mRNA levels. We show that in Saccharomyces cerevisiae TR_mD represents ∼20% of the variance in translation and directs an amplification exponent of 1.20 with a 95% confidence interval [1.14, 1.26]. TR_mIND constitutes the remaining ∼80% of the variance in translation and explains ∼5% of the variance in protein expression. We also find that TR_mD and TR_mIND are preferentially determined by different mRNA sequence features: TR_mIND by the length of the open reading frame and TR_mD both by a ∼60 nucleotide element that spans the initiating AUG and by codon and amino acid frequency. Our work provides more appropriate estimates of translational control and implies that TR_mIND is under different evolutionary selective pressures than TR_mD.

Metabolic differentiation of surface and invasive cells of yeast colony biofilms revealed by gene expression profiling

Background

Yeast infections are often connected with formation of biofilms that are extremely difficult to eradicate. An excellent model system for deciphering multifactorial determinants of yeast biofilm development is the colony biofilm, composed of surface (“aerial”) and invasive (“root”) cells. While surface cells have been partially analyzed before, we know little about invasive root cells. In particular, information on the metabolic, chemical and morphogenetic properties of invasive versus surface cells is lacking. In this study, we used a new strategy to isolate invasive cells from agar and extracellular matrix, and employed it to perform genome wide expression profiling and biochemical analyses of surface and invasive cells.

Results

RNA sequencing revealed expression differences in 1245 genes with high statistical significance, indicating large genetically regulated metabolic differences between surface and invasive cells. Functional annotation analyses implicated genes involved in stress defense, peroxisomal fatty acid β-oxidation, autophagy, protein degradation, storage compound metabolism and meiosis as being important in surface cells. In contrast, numerous genes with functions in nutrient transport and diverse synthetic metabolic reactions, including genes involved in ribosome biogenesis, biosynthesis and translation, were found to be important in invasive cells. Variation in gene expression correlated significantly with cell-type specific processes such as autophagy and storage compound accumulation as identified by microscopic and biochemical analyses. Expression profiling also provided indications of cell-specific regulations. Subsequent knockout strain analyses identified Gip2p, a regulatory subunit of type 1 protein phosphatase Glc7p, to be essential for glycogen accumulation in surface cells.

Conclusions

This is the first study reporting genome wide differences between surface and invasive cells of yeast colony biofilms. New findings show that surface and invasive cells display very different physiology, adapting to different conditions in different colony areas and contributing to development and survival of the colony biofilm as a whole. Notably, surface and invasive cells of colony biofilms differ significantly from upper and lower cells of smooth colonies adapted to plentiful laboratory conditions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4214-4) contains supplementary material, which is available to authorized users.

Similarity-Based Segmentation of Multi-Dimensional Signals

The segmentation of time series and genomic data is a common problem in computational biology. With increasingly complex measurement procedures individual data points are often not just numbers or simple vectors in which all components are of the same kind. Analysis methods that capitalize on slopes in a single real-valued data track or that make explicit use of the vectorial nature of the data are not applicable in such scenaria. We develop here a framework for segmentation in arbitrary data domains that only requires a minimal notion of similarity. Using unsupervised clustering of (a sample of) the input yields an approximate segmentation algorithm that is efficient enough for genome-wide applications. As a showcase application we segment a time-series of transcriptome sequencing data from budding yeast, in high temporal resolution over ca. 2.5 cycles of the short-period respiratory oscillation. The algorithm is used with a similarity measure focussing on periodic expression profiles across the metabolic cycle rather than coverage per time point.

Novel Stress-Inducible Antisense RNAs of Protein-Coding Loci Are Synthesized by RNA-Dependent RNA Polymerase

Our previous study identified approximately 6,000 abiotic stress-responsive noncoding transcripts existing on the antisense strand of protein-coding genes and implied that a type of antisense RNA was synthesized from a sense RNA template by RNA-dependent RNA polymerase (RDR). Expression analyses revealed that the expression of novel abiotic stress-induced antisense RNA on 1,136 gene loci was reduced in the rdr1/2/6 mutants. RNase protection indicated that the RD29A antisense RNA and other RDR1/2/6-dependent antisense RNAs are involved in the formation of dsRNA. The accumulation of stress-inducible antisense RNA was decreased and increased in dcp5 and xrn4, respectively, but not changed in dcl2/3/4, nrpd1a and nrpd1b RNA-seq analyses revealed that the majority of the RDR1/2/6-dependent antisense RNA loci did not overlap with RDR1/2/6-dependent 20-30 nt RNA loci. Additionally, rdr1/2/6 mutants decreased the degradation rate of the sense RNA and exhibited arrested root growth during the recovery stage following a drought stress, whereas dcl2/3/4 mutants did not. Collectively, these results indicate that RDRs have stress-inducible antisense RNA synthesis activity and a novel biological function that is different from the known endogenous small RNA pathways from protein-coding genes. These data reveal a novel mechanism of RNA regulation during abiotic stress response that involves complex RNA degradation pathways.

Natural antisense transcripts are significantly involved in regulation of drought stress in maize

Natural antisense transcripts (NATs) are a prominent and complex class of regulatory RNAs. Using strand-specific RNA sequencing, we identified 1769 sense and antisense transcript pairs (NAT pairs) in two maize inbreds with different sensitivity to drought, as well as in two derivative recombination inbred lines (RILs). A significantly higher proportion of NATs relative to non-NATs are specifically expressed under water stress (WS). Surprisingly, expression of sense and antisense transcripts produced by NAT pairs is significantly correlated, particularly under WS. We found an unexpected large proportion of NATs with protein coding potential, as estimated by ribosome release scores. Small RNAs significantly accumulate within NAT pairs, with 21 nt smRNA particularly enriched in overlapping regions of these pairs of genes. The abundance of these smRNAs is significantly altered in the leafbladeless1 mutant, suggesting that these genes may be regulated by the tasiRNA pathway. Further, NATs are significantly hypomethylated and include fewer transposable element sequences relative to non-NAT genes. NAT gene regions also exhibit higher levels of H3K36me3, H3K9ac, and H3K4me3, but lower levels of H3K27me3, indicating that NAT gene pairs generally exhibit an open chromatin configuration. Finally, NAT pairs in 368 diverse maize inbreds and 19 segregating populations were specifically enriched for polymorphisms associated with drought tolerance. Taken together, the data highlight the potential impact of that small RNAs and histone modifications have in regulation of NAT expression, and the significance of NATs in response to WS.

Yeast RNA-Binding Protein Nab3 Regulates Genes Involved in Nitrogen Metabolism

Termination of Saccharomyces cerevisiae RNA polymerase II (Pol II) transcripts occurs through two alternative pathways. Termination of mRNAs is coupled to cleavage and polyadenylation while noncoding transcripts are terminated through the Nrd1-Nab3-Sen1 (NNS) pathway in a process that is linked to RNA degradation by the nuclear exosome. Some mRNA transcripts are also attenuated through premature termination directed by the NNS complex. In this paper we present the results of nuclear depletion of the NNS component Nab3. As expected, many noncoding RNAs fail to terminate properly. In addition, we observe that nitrogen catabolite-repressed genes are upregulated by Nab3 depletion.

Genome-wide Identification and Characterization of Natural Antisense Transcripts by Strand-specific RNA Sequencing in Ganoderma lucidum

Ganoderma lucidum is a white-rot fungus best-known for its medicinal and ligninolytic activities. To discover the underlying genes responsible for these activities, we identified and characterized the natural antisense transcripts (NATs) using strand-specific (ss) RNA-seq data obtained from the mycelia, primordia and fruiting bodies. NATs were identified using a custom pipeline and then subjected to functional enrichment and differential expression analyses. A total of 1613 cis- and 244 trans- sense and antisense transcripts were identified. Mapping to GO terms and KEGG pathways revealed that NATs were frequently associated with genes of particular functional categories in particular stages. ssRT-qPCR experiments showed that the expression profiles of 30 of 50 (60%) transcripts were highly correlated with those of the RNA-seq results (r ≥ 0.9). Expression profiles of 22 of 25 (88%) pairs of NATs and STs were highly correlated (p ≤ 0.01), with 15 having r ≥ 0.8 and 4 having r ≤ -0.8. Six lignin-modifying genes and their NATs were analyzed in detail. Diverse patterns of differential expression among different stages and positive and negative correlations were observed. These results suggested that NATs were implicated in gene expression regulation in a function-group and developmental-stage specific manner through complex mechanisms.

Developmental transitions in Arabidopsis are regulated by antisense RNAs resulting from bidirectionally transcribed genes

Transcription terminators are DNA elements located at the 3' end of genes that ensure efficient cleavage of nascent RNA generating the 3' end of mRNA, as well as facilitating disengagement of elongating DNA-dependent RNA polymerase II. Surprisingly, terminators are also a potent source of antisense transcription. We have recently described an Arabidopsis antisense transcript originating from the 3' end of a master regulator of Arabidopsis thaliana seed dormancy DOG1. In this review, we discuss the broader implications of our discovery in light of recent developments in yeast and Arabidopsis. We show that, surprisingly, the key features of terminators that give rise to antisense transcription are preserved between Arabidopsis and yeast, suggesting a conserved mechanism. We also compare our discovery to known antisense-based regulatory mechanisms, highlighting the link between antisense-based gene expression regulation and major developmental transitions in plants.

Kinetic CRAC uncovers a role for Nab3 in determining gene expression profiles during stress

RNA-binding proteins play a key role in shaping gene expression profiles during stress, however, little is known about the dynamic nature of these interactions and how this influences the kinetics of gene expression. To address this, we developed kinetic cross-linking and analysis of cDNAs (χCRAC), an ultraviolet cross-linking method that enabled us to quantitatively measure the dynamics of protein-RNA interactions in vivo on a minute time-scale. Here, using χCRAC we measure the global RNA-binding dynamics of the yeast transcription termination factor Nab3 in response to glucose starvation. These measurements reveal rapid changes in protein-RNA interactions within 1 min following stress imposition. Changes in Nab3 binding are largely independent of alterations in transcription rate during the early stages of stress response, indicating orthogonal transcriptional control mechanisms. We also uncover a function for Nab3 in dampening expression of stress-responsive genes. χCRAC has the potential to greatly enhance our understanding of in vivo dynamics of protein-RNA interactions.Protein RNA interactions are dynamic and regulated in response to environmental changes. Here the authors describe 'kinetic CRAC', an approach that allows time resolved analyses of protein RNA interactions with minute time point resolution and apply it to gain insight into the function of the RNA-binding protein Nab3.

Backmasking in the yeast genome: encoding overlapping information for protein-coding and RNA degradation

Backmasking is a recording technique used to hide a sound or message in a music track in reverse, meaning that it is only audible when the record is played backwards. Analogously, the compact yeast genome encodes for diverse sources of information such as overlapping coding and non-coding transcripts, and protein-binding sites on the two complementary DNA strands. Examples are the consensus binding site sequences of the RNA-binding proteins Nrd1 and Nab3 that target non-coding transcripts for degradation. Here, by examining the overlap of stable (SUTs, stable unannotated transcripts) and unstable (CUTs, cryptic unstable transcripts) transcripts with protein-coding genes, we show that the predicted Nrd1 and Nab3-binding site sequences occur at differing frequencies. They are always depleted in the sense direction of protein-coding genes, thus avoiding degradation of the transcript. However in the antisense direction, predicted binding sites occur at high frequencies in genes with overlapping unstable ncRNAs (CUTs), so limiting the availability of non-functional transcripts. In contrast they are depleted in genes with overlapping stable ncRNAs (SUTs), presumably to avoid degrading the non-coding transcript. The protein-coding genes maintain similar amino-acid contents, but they display distinct codon usages so that Nrd1 and Nab3-binding sites can arise at differing frequencies in antisense depending on the overlapping transcript type. Our study demonstrates how yeast has evolved to encode multiple layers of information-protein-coding genes in one strand and the relative chance of degrading antisense RNA in the other strand-in the same regions of a compact genome.

Transcription-coupled changes to chromatin underpin gene silencing by transcriptional interference

Long non-coding RNA (lncRNA) transcription into a downstream promoter frequently results in transcriptional interference. However, the mechanism of this repression is not fully understood. We recently showed that drug tolerance in fission yeast Schizosaccharomyces pombe is controlled by lncRNA transcription upstream of the tgp1⁺ permease gene. Here we demonstrate that transcriptional interference of tgp1⁺ involves several transcription-coupled chromatin changes mediated by conserved elongation factors Set2, Clr6CII, Spt6 and FACT. These factors are known to travel with RNAPII and establish repressive chromatin in order to limit aberrant transcription initiation from cryptic promoters present in gene bodies. We therefore conclude that conserved RNAPII-associated mechanisms exist to both suppress intragenic cryptic promoters during genic transcription and to repress gene promoters by transcriptional interference. Our analyses also demonstrate that key mechanistic features of transcriptional interference are shared between S. pombe and the highly divergent budding yeast Saccharomyces cerevisiae. Thus, transcriptional interference is an ancient, conserved mechanism for tightly controlling gene expression. Our mechanistic insights allowed us to predict and validate a second example of transcriptional interference involving the S. pombe pho1⁺ gene. Given that eukaryotic genomes are pervasively transcribed, transcriptional interference likely represents a more general feature of gene regulation than is currently appreciated.

Expansion of antisense lncRNA transcriptomes in budding yeast species since the loss of RNAi

Antisense long noncoding RNAs (ASlncRNAs) have been implicated in regulating gene expression in response to physiological cues. However, little is known about the evolutionary dynamics of ASlncRNA and what underlies the evolution of their expression. Here, using budding yeast Saccharomyces spp. and Naumovozyma castellii as models, we show that ASlncRNA repertoires have expanded since the loss of RNA interference (RNAi), in terms of their expression levels, their lengths and their degree of overlap with coding genes. Furthermore, we show that RNAi is inhibitory to ASlncRNA transcriptomes and that increased expression of ASlncRNAs in the presence of RNAi is deleterious to N. castellii, which has retained RNAi. Together, our results suggest that the loss of RNAi had substantial effects on the genome-wide increase in expression of ASlncRNAs during the evolution of budding yeasts.

Survey of cryptic unstable transcripts in yeast

Background

Cryptic unstable transcripts (CUTs) are a largely unexplored class of nuclear exosome degraded, non-coding RNAs in budding yeast. It is highly debated whether CUT transcription has a functional role in the cell or whether CUTs represent noise in the yeast transcriptome.

We sought to ascertain the extent of conserved CUT expression across a variety of Saccharomyces yeast strains to further understand and characterize the nature of CUT expression.

Results

We sequenced the WT and rrp6Δ transcriptomes of three S.cerevisiae strains: S288c, Σ1278b, JAY291 and the S.paradoxus strain N17 and utilized a hidden Markov model to annotate CUTs in these four strains. Utilizing a four-way genomic alignment we identified a large population of CUTs with conserved syntenic expression across all four strains. By identifying configurations of gene-CUT pairs, where CUT expression originates from the gene 5’ or 3′ nucleosome free region, we observed distinct gene expression trends specific to these configurations which were most prevalent in the presence of conserved CUT expression. Divergent pairs correlate with higher expression of genes, and convergent pairs correlate with reduced gene expression.

Conclusions

Our RNA-seq based method has greatly expanded upon previous CUT annotations in S.cerevisiae underscoring the extensive and pervasive nature of unstable transcription. Furthermore we provide the first assessment of conserved CUT expression in yeast and globally demonstrate possible modes of CUT-based regulation of gene expression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2622-5) contains supplementary material, which is available to authorized users.

Early-branching gut fungi possess a large, comprehensive array of biomass-degrading enzymes

The fungal kingdom is the source of almost all industrial enzymes in use for lignocellulose bioprocessing. We developed a systems-level approach that integrates transcriptomic sequencing, proteomics, phenotype, and biochemical studies of relatively unexplored basal fungi. Anaerobic gut fungi isolated from herbivores produce a large array of biomass-degrading enzymes that synergistically degrade crude, untreated plant biomass and are competitive with optimized commercial preparations from Aspergillus and Trichoderma. Compared to these model platforms, gut fungal enzymes are unbiased in substrate preference due to a wealth of xylan-degrading enzymes. These enzymes are universally catabolite-repressed and are further regulated by a rich landscape of noncoding regulatory RNAs. Additionally, we identified several promising sequence-divergent enzyme candidates for lignocellulosic bioprocessing.

Functional elucidation of the non-coding RNAs of Kluyveromyces marxianus in the exponential growth phase

BACKGROUND

Non-coding RNAs (ncRNAs), which perform diverse regulatory roles, have been found in organisms from all superkingdoms of life. However, there have been limited numbers of studies on the functions of ncRNAs, especially in nonmodel organisms such as Kluyveromyces marxianus that is widely used in the field of industrial biotechnology.

RESULTS

In this study, we measured changes in transcriptome at three time points during the exponential growth phase of K. marxianus by using strand-specific RNA-seq. We found that approximately 60% of the transcriptome consists of ncRNAs transcribed from antisense and intergenic regions of the genome that were transcribed at lower levels than mRNA. In the transcriptome, a substantial number of long antisense ncRNAs (lancRNAs) are differentially expressed and enriched in carbohydrate and energy metabolism pathways. Furthermore, this enrichment is evolutionarily conserved, at least in yeast. Particularly, the mode of regulation of mRNA/lancRNA pairs is associated with mRNA transcription levels; the correlation between the pairs is positive at high mRNA transcriptional levels and negative at low levels. In addition, significant induction of mRNA and coverage of more than half of the mRNA sequence by a lancRNA strengthens the positive correlation between mRNA/lancRNA pairs.

CONCLUSIONS

Transcriptome sequencing of K. marxianus in the exponential growth phase reveals pervasive transcription of ncRNAs with evolutionarily conserved functions. Studies of the mode of regulation of mRNA/lancRNA pairs suggest that induction of lancRNA may be associated with switch-like behavior of mRNA/lancRNA pairs and efficient regulation of the carbohydrate and energy metabolism pathways in the exponential growth phase of K. marxianus being used in industrial applications.

Loss of the Yeast SR Protein Npl3 Alters Gene Expression Due to Transcription Readthrough

Yeast Npl3 is a highly abundant, nuclear-cytoplasmic shuttling, RNA-binding protein, related to metazoan SR proteins. Reported functions of Npl3 include transcription elongation, splicing and RNA 3’ end processing. We used UV crosslinking and analysis of cDNA (CRAC) to map precise RNA binding sites, and strand-specific tiling arrays to look at the effects of loss of Npl3 on all transcripts across the genome. We found that Npl3 binds diverse RNA species, both coding and non-coding, at sites indicative of roles in both early pre-mRNA processing and 3’ end formation. Tiling arrays and RNAPII mapping data revealed 3’ extended RNAPII-transcribed RNAs in the absence of Npl3, suggesting that defects in pre-mRNA packaging events result in termination readthrough. Transcription readthrough was widespread and frequently resulted in down-regulation of neighboring genes. We conclude that the absence of Npl3 results in widespread 3' extension of transcripts with pervasive effects on gene expression.

Npl3 is a yeast mRNA binding protein with many reported functions in RNA processing. We wanted to identify direct targets and therefore combined analyses of the transcriptome-wide effects of the loss of Npl3 on gene expression with UV crosslinking and bioinformatics to identify RNA-binding sites for Npl3. We found that Npl3 binds diverse sites on large numbers of transcripts, and that the loss of Npl3 results in transcriptional readthrough on many genes. One effect of this transcription readthrough is that the expression of numerous flanking genes is strongly down regulated. This underlines the importance of faithful termination for the correct regulation of gene expression. The effects of the loss of Npl3 are seen on both mRNAs and non-protein coding RNAs. These have distinct but overlapping termination mechanisms, with both classes requiring Npl3 for correct RNA packaging.

Evolution of regulatory networks in Candida glabrata: learning to live with the human host

The opportunistic human fungal pathogen Candida glabrata is second only to C. albicans as the cause of Candida infections and yet is more closely related to Saccharomyces cerevisiae. Recent advances in functional genomics technologies and computational approaches to decipher regulatory networks, and the comparison of these networks among these and other Ascomycete species, have revealed both unique and shared strategies in adaptation to a human commensal/opportunistic pathogen lifestyle and antifungal drug resistance in C. glabrata. Recently, several C. glabrata sister species in the Nakeseomyces clade representing both human associated (commensal) and environmental isolates have had their genomes sequenced and analyzed. This has paved the way for comparative functional genomics studies to characterize the regulatory networks in these species to identify informative patterns of conservation and divergence linked to phenotypic evolution in the Nakaseomyces lineage.

An Ancient Yeast for Young Geneticists: A Primer on the Schizosaccharomyces pombe Model System

The fission yeast Schizosaccharomyces pombe is an important model organism for the study of eukaryotic molecular and cellular biology. Studies of S. pombe, together with studies of its distant cousin, Saccharomyces cerevisiae, have led to the discovery of genes involved in fundamental mechanisms of transcription, translation, DNA replication, cell cycle control, and signal transduction, to name but a few processes. However, since the divergence of the two species approximately 350 million years ago, S. pombe appears to have evolved less rapidly than S. cerevisiae so that it retains more characteristics of the common ancient yeast ancestor, causing it to share more features with metazoan cells. This Primer introduces S. pombe by describing the yeast itself, providing a brief description of the origins of fission yeast research, and illustrating some genetic and bioinformatics tools used to study protein function in fission yeast. In addition, a section on some key differences between S. pombe and S. cerevisiae is included for readers with some familiarity with budding yeast research but who may have an interest in developing research projects using S. pombe.

The Set3 Complex Antagonizes the MYST Acetyltransferase Esa1 in the DNA Damage Response

Acetylation is a dynamic posttranslational modification that contributes to chromatin-regulated processes, including DNA replication, repair, recombination, and gene expression. Acetylation is controlled by complexes containing opposing lysine and histone acetyltransferase (KAT and HAT) and deacetylase (KDAC and HDAC) activities. The essential MYST family Esa1 KAT acetylates core histones and many nonhistone substrates. Phenotypes of esa1 mutants include transcriptional silencing and activation defects, impaired growth at high temperatures, and sensitivity to DNA damage. The KDAC Rpd3 was previously identified as an activity opposing Esa1, as its deletion suppresses growth and silencing defects of esa1 mutants. However, loss of Rpd3 does not suppress esa1 DNA damage sensitivity. In this work, we identified Hos2 as a KDAC counteracting ESA1 in the damage response. Deletion of HOS2 resulted in changes of esa1's transcriptional response upon damage. Further, loss of HOS2 or components of the Set3 complex (Set3C) in which it acts specifically suppressed damage sensitivity and restored esa1 histone H4 acetylation. This rescue was mediated via loss of either Set3C integrity or of its binding to dimethylated histone H3K4. Our results thus add new insight into the interactions of an essential MYST acetyltransferase with diverse deacetylases to respond specifically to environmental and physiological challenges.

Genome analysis of rice-blast fungus Magnaporthe oryzae field isolates from southern India

The Indian subcontinent is the center of origin and diversity for rice (Oryza sativa L.). The O. sativa ssp. indica is a major food crop grown in India, which occupies the first and second position in area and production, respectively. Blast disease caused by Magnaporthe oryzae is a major constraint to rice production. Here, we report the analysis of genome architecture and sequence variation of two field isolates, B157 and MG01, of the blast fungus from southern India. The 40 Mb genome of B157 and 43 Mb genome of MG01 contained 11,344 and 11,733 predicted genes, respectively. Genomic comparisons unveiled a large set of SNPs and several isolate specific genes in the Indian blast isolates. Avr genes were analyzed in several sequenced Magnaporthe strains; this analysis revealed the presence of Avr-Pizt and Avr-Ace1 genes in all the sequenced isolates. Availability of whole genomes of field isolates from India will contribute to global efforts to understand genetic diversity of M. oryzae population and to track the emergence of virulent pathotypes.

•

The first genomic study of Magnaporthe from Indian subcontinent

•

Provided information about genomic variations in terms of SNPs, InDels and ICVs due to transposable elements

•

Identified novel genes specific to Indian isolates

•

Genome wide antisense transcripts identified from this study

•

Identified Magnaporthe specific pathogenicity genes that are absent in non-pathogenic Ascomycetes fungi

The landscape of antisense gene expression in human cancers

High-throughput RNA sequencing has revealed more pervasive transcription of the human genome than previously anticipated. However, the extent of natural antisense transcripts’ (NATs) expression, their regulation of cognate sense genes, and the role of NATs in cancer remain poorly understood. Here, we use strand-specific paired-end RNA sequencing (ssRNA-seq) data from 376 cancer samples covering nine tissue types to comprehensively characterize the landscape of antisense expression. We found consistent antisense expression in at least 38% of annotated transcripts, which in general is positively correlated with sense gene expression. Investigation of sense/antisense pair expressions across tissue types revealed lineage-specific, ubiquitous and cancer-specific antisense loci transcription. Comparisons between tumor and normal samples identified both concordant (same direction) and discordant (opposite direction) sense/antisense expression patterns. Finally, we provide OncoNAT, a catalog of cancer-related genes with significant antisense transcription, which will enable future investigations of sense/antisense regulation in cancer. Using OncoNAT we identified several functional NATs, including NKX2-1-AS1 that regulates the NKX2-1 oncogene and cell proliferation in lung cancer cells. Overall, this study provides a comprehensive account of NATs and supports a role for NATs' regulation of tumor suppressors and oncogenes in cancer biology.