Alternative Splicing in Next Generation Sequencing Data of Saccharomyces cerevisiae

Psr1 phosphatase regulates pre-mRNA splicing through spliceosomal B complex factor Snu66

Regulated precursor messenger RNA (pre-mRNA) splicing modulates gene expression and promotes alternative splicing. The process is regulated by modifications of spliceosomal proteins and small nuclear RNAs (snRNAs). Here, we show that the protein phosphatase Psr1, known for its plasma membrane localisation and function in general stress response in Saccharomyces cerevisiae, also plays a regulatory role in pre-mRNA splicing. Independently of its presence at the plasma membrane, Psr1 binds and dephosphorylates the core splicing factor Snu66. The enzyme is not an integral component of the spliceosome. Psr1 deletion in yeast, or tethering of its catalytic mutant to Snu66, results in splicing defects of introns with non-canonical 5' splice sites (ss). While the Psr1 binding site on Snu66 is distinct from the Hub1 interaction domains (HIND), Hub1 displaces Psr1 from Snu66. Thus, Psr1 phosphatase plays a regulatory role in pre-mRNA splicing by modulating Snu66 functions.

Long-read direct RNA sequencing of the mitochondrial transcriptome of Saccharomyces cerevisiae reveals condition-dependent intron abundance

Mitochondria fulfil many essential roles and have their own genome, which is expressed as polycistronic transcripts that undergo co- or posttranscriptional processing and splicing. Due to the inherent complexity and limited technical accessibility of the mitochondrial transcriptome, fundamental questions regarding mitochondrial gene expression and splicing remain unresolved, even in the model eukaryote Saccharomyces cerevisiae. Long-read sequencing could address these fundamental questions. Therefore, a method for the enrichment of mitochondrial RNA and sequencing using Nanopore technology was developed, enabling the resolution of splicing of polycistronic genes and the quantification of spliced RNA. This method successfully captured the full mitochondrial transcriptome and resolved RNA splicing patterns with single-base resolution and was applied to explore the transcriptome of S. cerevisiae grown with glucose or ethanol as the sole carbon source, revealing the impact of growth conditions on mitochondrial RNA expression and splicing. This study uncovered a remarkable difference in the turnover of Group II introns between yeast grown in either mostly fermentative or fully respiratory conditions. Whether this accumulation of introns in glucose medium has an impact on mitochondrial functions remains to be explored. Combined with the high tractability of the model yeast S. cerevisiae, the developed method enables to monitor mitochondrial transcriptome responses in a broad range of relevant contexts, including oxidative stress, apoptosis and mitochondrial diseases.

Replicative aging in yeast involves dynamic intron retention patterns associated with mRNA processing/export and protein ubiquitination

Saccharomyces cerevisiae (baker's yeast) has yielded relevant insights into some of the basic mechanisms of organismal aging. Among these are genomic instability, oxidative stress, caloric restriction and mitochondrial dysfunction. Several genes are known to have an impact on the aging process, with corresponding mutants exhibiting short- or long-lived phenotypes. Research dedicated to unraveling the underlying cellular mechanisms can support the identification of conserved mechanisms of aging in other species. One of the hitherto less studied fields in yeast aging is how the organism regulates its gene expression at the transcriptional level. To our knowledge, we present the first investigation into alternative splicing, particularly intron retention, during replicative aging of S. cerevisiae. This was achieved by utilizing the IRFinder algorithm on a previously published RNA-seq data set by Janssens et al. (2015). In the present work, 44 differentially retained introns in 43 genes were identified during replicative aging. We found that genes with altered intron retention do not display significant changes in overall transcript levels. It was possible to functionally assign distinct groups of these genes to the cellular processes of mRNA processing and export (e.g., YRA1) in early and middle-aged yeast, and protein ubiquitination (e.g., UBC5) in older cells. In summary, our work uncovers a previously unexplored layer of the transcriptional program of yeast aging and, more generally, expands the knowledge on the occurrence of alternative splicing in baker's yeast.

Splicing factor Prp18p promotes genome-wide fidelity of consensus 3'-splice sites

The fidelity of splice site selection is critical for proper gene expression. In particular, proper recognition of 3'-splice site (3'SS) sequences by the spliceosome is challenging considering the low complexity of the 3'SS consensus sequence YAG. Here, we show that absence of the Prp18p splicing factor results in genome-wide activation of alternative 3'SS in S. cerevisiae, including highly unusual non-YAG sequences. Usage of these non-canonical 3'SS in the absence of Prp18p is enhanced by upstream poly(U) tracts and by their potential to interact with the first intronic nucleoside, allowing them to dock in the spliceosome active site instead of the normal 3'SS. The role of Prp18p in 3'SS fidelity is facilitated by interactions with Slu7p and Prp8p, but cannot be fulfilled by Slu7p, identifying a unique role for Prp18p in 3'SS fidelity. This fidelity function is synergized by the downstream proofreading activity of the Prp22p helicase, but is independent from another late splicing helicase, Prp43p. Our results show that spliceosomes exhibit remarkably relaxed 3'SS sequence usage in the absence of Prp18p and identify a network of spliceosomal interactions centered on Prp18p which are required to promote the fidelity of the recognition of consensus 3'SS sequences.

RNA recognition by Npl3p reveals U2 snRNA-binding compatible with a chaperone role during splicing

The conserved SR-like protein Npl3 promotes splicing of diverse pre-mRNAs. However, the RNA sequence(s) recognized by the RNA Recognition Motifs (RRM1 & RRM2) of Npl3 during the splicing reaction remain elusive. Here, we developed a split-iCRAC approach in yeast to uncover the consensus sequence bound to each RRM. High-resolution NMR structures show that RRM2 recognizes a 5´-GNGG-3´ motif leading to an unusual mille-feuille topology. These structures also reveal how RRM1 preferentially interacts with a CC-dinucleotide upstream of this motif, and how the inter-RRM linker and the region C-terminal to RRM2 contribute to cooperative RNA-binding. Structure-guided functional studies show that Npl3 genetically interacts with U2 snRNP specific factors and we provide evidence that Npl3 melts U2 snRNA stem-loop I, a prerequisite for U2/U6 duplex formation within the catalytic center of the Bact spliceosomal complex. Thus, our findings suggest an unanticipated RNA chaperoning role for Npl3 during spliceosome active site formation.

Promoters and introns as key drivers for enhanced gene expression in Saccharomyces cerevisiae

The transcription of genes in the yeast Saccharomyces cerevisiae is governed by multiple layers of regulatory elements and proteins, cooperating to ensure optimum expression of the final protein product based on the cellular requirements. Promoters have always been regarded as the most important determinant of gene transcription, but introns also play a key role in the expression of intron-encoding genes. Some introns can enhance transcription when introduced either promoter-proximal or embedded in the open reading frame of genes. However, the outcome is seldom predictable, with some introns increasing or decreasing transcription depending on the promoter and reporter gene employed. This chapter provides an overview of the general structure and function of promoters and introns and how they may cooperate during transcription to allow intron-mediated enhancement of gene expression. Since S. cerevisiae is a suitable host for recombinant protein production on a commercial level, stronger and more controllable promoters are in high demand. Enhanced gene expression can be achieved via promoter engineering, which may include introns that increase the efficacy of recombinant expression cassettes. Different models for the role of introns in transcription are briefly discussed to show how these intervening sequences can actively interact with the transcription machinery. Furthermore, recent examples of improved protein production via the introduction of promoter-proximal introns are highlighted to showcase the potential value of intron-mediated enhancement of gene expression.

Regulation of 3' splice site selection after step 1 of splicing by spliceosomal C* proteins

Alternative precursor messenger RNA splicing is instrumental in expanding the proteome of higher eukaryotes, and changes in 3' splice site (3'ss) usage contribute to human disease. We demonstrate by small interfering RNA-mediated knockdowns, followed by RNA sequencing, that many proteins first recruited to human C* spliceosomes, which catalyze step 2 of splicing, regulate alternative splicing, including the selection of alternatively spliced NAGNAG 3'ss. Cryo-electron microscopy and protein cross-linking reveal the molecular architecture of these proteins in C* spliceosomes, providing mechanistic and structural insights into how they influence 3'ss usage. They further elucidate the path of the 3' region of the intron, allowing a structure-based model for how the C* spliceosome potentially scans for the proximal 3'ss. By combining biochemical and structural approaches with genome-wide functional analyses, our studies reveal widespread regulation of alternative 3'ss usage after step 1 of splicing and the likely mechanisms whereby C* proteins influence NAGNAG 3'ss choices.

A dynamical stochastic model of yeast translation across the cell cycle

Translation is a central step in gene expression, however its quantitative and time-resolved regulation is poorly understood. We developed a discrete, stochastic model for protein translation in S. cerevisiae in a whole-transcriptome, single-cell context. A "base case" scenario representing an average cell highlights translation initiation rates as the main co-translational regulatory parameters. Codon usage bias emerges as a secondary regulatory mechanism through ribosome stalling. Demand for anticodons with low abundancy is shown to cause above-average ribosome dwelling times. Codon usage bias correlates strongly both with protein synthesis rates and elongation rates. Applying the model to a time-resolved transcriptome estimated by combining data from FISH and RNA-Seq experiments, it could be shown that increased total transcript abundance during the cell cycle decreases translation efficiency at single transcript level. Translation efficiency grouped by gene function shows highest values for ribosomal and glycolytic genes. Ribosomal proteins peak in S phase while glycolytic proteins rank highest in later cell cycle phases.

Spliceosomal mutations decouple 3' splice site fidelity from cellular fitness

The fidelity of splice site selection is thought to be critical for proper gene expression and cellular fitness. In particular, proper recognition of 3'-splice site (3'SS) sequences by the spliceosome is a daunting task considering the low complexity of the 3'SS consensus sequence YAG. Here we show that inactivating the near-essential splicing factor Prp18p results in a global activation of alternative 3'SS, many of which harbor sequences that highly diverge from the YAG consensus, including some highly unusual non-AG 3'SS. We show that the role of Prp18p in 3'SS fidelity is promoted by physical interactions with the essential splicing factors Slu7p and Prp8p and synergized by the proofreading activity of the Prp22p helicase. Strikingly, structure-guided point mutations that disrupt Prp18p-Slu7p and Prp18p-Prp8p interactions mimic the loss of 3'SS fidelity without any impact on cellular growth, suggesting that accumulation of incorrectly spliced transcripts does not have a major deleterious effect on cellular viability. These results show that spliceosomes exhibit remarkably relaxed fidelity in the absence of Prp18p, and that new 3'SS sampling can be achieved genome-wide without a major negative impact on cellular fitness, a feature that could be used during evolution to explore new productive alternative splice sites.

Intron distribution and emerging role of alternative splicing in fungi

Spliceosomal introns are noncoding sequences that are spliced from pre-mRNA. They are ubiquitous in eukaryotic genomes, although the average number of introns per gene varies considerably between different eukaryotic species. Fungi are diverse in terms of intron numbers ranging from 4% to 99% genes with introns. Alternative splicing is one of the most common modes of posttranscriptional regulation in eukaryotes, giving rise to multiple transcripts from a single pre-mRNA and is widespread in metazoans and drives extensive proteome diversity. Earlier, alternative splicing was considered to be rare in fungi, but recently, increasing numbers of studies have revealed that alternative splicing is also widespread in fungi and has been implicated in the regulation of fungal growth and development, protein localization and the improvement of survivability, likely underlying their unique capacity to adapt to changing environmental conditions. However, the role of alternative splicing in pathogenicity and development of drug resistance is only recently gaining attention. In this review, we describe the intronic landscape in fungi. We also present in detail the newly discovered functions of alternative splicing in various cellular processes and outline areas particularly in pathogenesis and clinical drug resistance for future studies that could lead to the development of much needed new therapeutics.

Post-transcriptional and translational control of the morphology and virulence in human fungal pathogens

Host-pathogen interactions at the molecular level are the key to fungal pathogenesis. Fungal pathogens utilize several mechanisms such as adhesion, invasion, phenotype switching and metabolic adaptations, to survive in the host environment and respond. Post-transcriptional and translational regulations have emerged as key regulatory mechanisms ensuring the virulence and survival of fungal pathogens. Through these regulations, fungal pathogens effectively alter their protein pool, respond to various stress, and undergo morphogenesis, leading to efficient and comprehensive changes in fungal physiology. The regulation of virulence through post-transcriptional and translational regulatory mechanisms is mediated through mRNA elements (cis factors) or effector molecules (trans factors). The untranslated regions upstream and downstream of the mRNA, as well as various RNA-binding proteins involved in translation initiation or circularization of the mRNA, play pivotal roles in the regulation of morphology and virulence by influencing protein synthesis, protein isoforms, and mRNA stability. Therefore, post-transcriptional and translational mechanisms regulating the morphology, virulence and drug-resistance processes in fungal pathogens can be the target for new therapeutics. With improved "omics" technologies, these regulatory mechanisms are increasingly coming to the forefront of basic biology and drug discovery. This review aims to discuss various modes of post-transcriptional and translation regulations, and how these mechanisms exert influence in the virulence and morphogenesis of fungal pathogens.

Origin, conservation, and loss of alternative splicing events that diversify the proteome in Saccharomycotina budding yeasts

Many eukaryotes use RNA processing, including alternative splicing, to express multiple gene products from the same gene. The budding yeast Saccharomyces cerevisiae has been successfully used to study the mechanism of splicing and the splicing machinery, but alternative splicing in yeast is relatively rare and has not been extensively studied. Alternative splicing of SKI7/HBS1 is widely conserved, but yeast and a few other eukaryotes have replaced this one alternatively spliced gene with a pair of duplicated, unspliced genes as part of a whole genome doubling (WGD). We show that other examples of alternative splicing known to have functional consequences are widely conserved within Saccharomycotina. A common mechanism by which alternative splicing has disappeared is by replacement of an alternatively spliced gene with duplicate unspliced genes. This loss of alternative splicing does not always take place soon after duplication, but can take place after sufficient time has elapsed for speciation. Saccharomycetaceae that diverged before WGD use alternative splicing more frequently than S. cerevisiae, suggesting that WGD is a major reason for infrequent alternative splicing in yeast. We anticipate that WGDs in other lineages may have had the same effect. Having observed that two functionally distinct splice-isoforms are often replaced by duplicated genes allowed us to reverse the reasoning. We thereby identify several splice isoforms that are likely to produce two functionally distinct proteins because we find them replaced by duplicated genes in related species. We also identify some alternative splicing events that are not conserved in closely related species and unlikely to produce functionally distinct proteins.

Identification of Genomewide Alternative Splicing Events in Sequential, Isogenic Clinical Isolates of Candida albicans Reveals a Novel Mechanism of Drug Resistance and Tolerance to Cellular Stresses

The emergence of resistance in Candida albicans, an opportunistic pathogen, against the commonly used antifungals is becoming a major obstacle in its treatment. The necessity to identify new drug targets demands fundamental insights into the mechanisms used by this organism to develop drug resistance. C. albicans has introns in 4 to 6% of its genes, the functions of which remain largely unknown. Using the RNA-sequencing data from isogenic pairs of azole-sensitive and -resistant isolates of C. albicans, here, we show how C. albicans uses modulations in mRNA splicing to overcome antifungal drug stress.

Alternative splicing (AS)—a process by which a single gene gives rise to different protein isoforms in eukaryotes—has been implicated in many basic cellular processes, but little is known about its role in drug resistance and fungal pathogenesis. The most common human fungal pathogen, Candida albicans, has introns in 4 to 6% of its genes, the functions of which remain largely unknown. Here, we report AS regulating drug resistance in C. albicans. Comparative RNA-sequencing of two different sets of sequential, isogenic azole-sensitive and -resistant isolates of C. albicans revealed differential expression of splice isoforms of 14 genes. One of these was the superoxide dismutase gene SOD3, which contains a single intron. The sod3Δ/Δ mutant was susceptible to the antifungals amphotericin B (AMB) and menadione (MND). While AMB susceptibility was rescued by overexpression of both the spliced and unspliced SOD3 isoforms, only the spliced isoform could overcome MND susceptibility, demonstrating the functional relevance of this splicing in developing drug resistance. Furthermore, unlike AMB, MND inhibits SOD3 splicing and acts as a splicing inhibitor. Consistent with these observations, MND exposure resulted in increased levels of unspliced SOD3 isoform that are unable to scavenge reactive oxygen species (ROS), resulting in increased drug susceptibility. Collectively, these observations suggest that AS is a novel mechanism for stress adaptation and overcoming drug susceptibility in C. albicans.

IMPORTANCE The emergence of resistance in Candida albicans, an opportunistic pathogen, against the commonly used antifungals is becoming a major obstacle in its treatment. The necessity to identify new drug targets demands fundamental insights into the mechanisms used by this organism to develop drug resistance. C. albicans has introns in 4 to 6% of its genes, the functions of which remain largely unknown. Using the RNA-sequencing data from isogenic pairs of azole-sensitive and -resistant isolates of C. albicans, here, we show how C. albicans uses modulations in mRNA splicing to overcome antifungal drug stress.

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

QTL analysis of natural Saccharomyces cerevisiae isolates reveals unique alleles involved in lignocellulosic inhibitor tolerance

Decoding the genetic basis of lignocellulosic inhibitor tolerance in Saccharomyces cerevisiae is crucial for rational engineering of bioethanol strains with enhanced robustness. The genetic diversity of natural strains present an invaluable resource for the exploration of complex traits of industrial importance from a pan-genomic perspective to complement the limited range of specialised, tolerant industrial strains. Natural S. cerevisiae isolates have lately garnered interest as a promising toolbox for engineering novel, genetically encoded tolerance phenotypes into commercial strains. To this end, we investigated the genetic basis for lignocellulosic inhibitor tolerance of natural S. cerevisiae isolates. A total of 12 quantitative trait loci underpinning tolerance were identified by next-generation sequencing linked bulk-segregant analysis of superior interbred pools. Our findings corroborate the current perspective of lignocellulosic inhibitor tolerance as a multigenic, complex trait. Apart from a core set of genetic variants required for inhibitor tolerance, an additional genetic background-specific response was observed. Functional analyses of the identified genetic loci revealed the uncharacterised ORF, YGL176C and the bud-site selection XRN1/BUD13 as potentially beneficial alleles contributing to tolerance to a complex lignocellulosic inhibitor mixture. We present evidence for the consideration of both regulatory and coding sequence variants for strain improvement.

5' untranslated regions: the next regulatory sequence in yeast synthetic biology

When developing industrial biotechnology processes, Saccharomyces cerevisiae (baker's yeast or brewer's yeast) is a popular choice as a microbial host. Many tools have been developed in the fields of synthetic biology and metabolic engineering to introduce heterologous pathways and tune their expression in yeast. Such tools mainly focus on controlling transcription, whereas post-transcriptional regulation is often overlooked. Herein we discuss regulatory elements found in the 5' untranslated region (UTR) and their influence on protein synthesis. We provide not only an overall picture, but also a set of design rules on how to engineer a 5' UTR. The reader is also referred to currently available models that allow gene expression to be tuned predictably using different 5' UTRs.

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.

Comparative Study on Alternative Splicing in Human Fungal Pathogens Suggests Its Involvement During Host Invasion

Alternative splicing (AS) is an important regulatory mechanism in eukaryotes but only little is known about its impact in fungi. Human fungal pathogens are of high clinical interest causing recurrent or life-threatening infections. AS can be well-investigated genome-wide and quantitatively with the powerful technology of RNA-Seq. Here, we systematically studied AS in human fungal pathogens based on RNA-Seq data. To do so, we investigated its effect in seven fungi during conditions simulating ex vivo infection processes and during in vitro stress. Genes undergoing AS are species-specific and act independently from differentially expressed genes pointing to an independent mechanism to change abundance and functionality. Candida species stand out with a low number of introns with higher and more varying lengths and more alternative splice sites. Moreover, we identified a functional difference between response to host and other stress conditions: During stress, AS affects more genes and is involved in diverse regulatory functions. In contrast, during response-to-host conditions, genes undergoing AS have membrane functionalities and might be involved in the interaction with the host. We assume that AS plays a crucial regulatory role in pathogenic fungi and is important in both response to host and stress conditions.

Proteome-scale relationships between local amino acid composition and protein fates and functions

Proteins with low-complexity domains continue to emerge as key players in both normal and pathological cellular processes. Although low-complexity domains are often grouped into a single class, individual low-complexity domains can differ substantially with respect to amino acid composition. These differences may strongly influence the physical properties, cellular regulation, and molecular functions of low-complexity domains. Therefore, we developed a bioinformatic approach to explore relationships between amino acid composition, protein metabolism, and protein function. We find that local compositional enrichment within protein sequences is associated with differences in translation efficiency, abundance, half-life, protein-protein interaction promiscuity, subcellular localization, and molecular functions of proteins on a proteome-wide scale. However, local enrichment of related amino acids is sometimes associated with opposite effects on protein regulation and function, highlighting the importance of distinguishing between different types of low-complexity domains. Furthermore, many of these effects are discernible at amino acid compositions below those required for classification as low-complexity or statistically-biased by traditional methods and in the absence of homopolymeric amino acid repeats, indicating that thresholds employed by classical methods may not reflect biologically relevant criteria. Application of our analyses to composition-driven processes, such as the formation of membraneless organelles, reveals distinct composition profiles even for closely related organelles. Collectively, these results provide a unique perspective and detailed insights into relationships between amino acid composition, protein metabolism, and protein functions.

Low-complexity domains in protein sequences are regions that are composed of only a few amino acids in the protein “alphabet”. These domains often have unique chemical properties and play important biological roles in both normal and disease-related processes. While a number of approaches have been developed to define low-complexity domains, these methods each possess conceptual limitations. Therefore, we developed a complementary approach that focuses on local amino acid composition (i.e. the amino acid composition within small regions of proteins). We find that high local composition of individual amino acids is associated with pervasive effects on protein metabolism, subcellular localization, and molecular function on a proteome-wide scale. Importantly, the nature of the effects depend on the type of amino acid enriched within the examined domains, and are observable in the absence of classically-defined low-complexity (and related) domains. Furthermore, we define the compositions of proteins involved in the formation of membraneless, protein-rich organelles such as stress granules and P-bodies. Our results provide a coherent view and unprecedented resolution of the effects of local amino acid enrichment on protein biology.

Regulation of lipid droplets by metabolically controlled Ldo isoforms

Storage and consumption of neutral lipids in lipid droplets (LDs) are essential for energy homeostasis and tightly coupled to cellular metabolism. However, how metabolic cues are integrated in the life cycle of LDs is unclear. In this study, we characterize the function of Ldo16 and Ldo45, two splicing isoforms of the same protein in budding yeast. We show that Ldo proteins interact with the seipin complex, which regulates contacts between LDs and the endoplasmic reticulum (ER). Moreover, we show that the levels of Ldo16 and Ldo45 depend on the growth stage of cells and that deregulation of their relative abundance alters LD morphology, protein localization, and triglyceride content. Finally, we show that absence of Ldo proteins results in defects in LD morphology and consumption by lipophagy. Our findings support a model in which Ldo proteins modulate the activity of the seipin complex, thereby affecting LD properties. Moreover, we identify ER-LD contacts as regulatory targets coupling energy storage to cellular metabolism.

Identification of seipin-linked factors that act as determinants of a lipid droplet subpopulation

Functional heterogeneity within the lipid droplet (LD) pool of a single cell has been observed, yet the underlying mechanisms remain enigmatic. Here, we report on identification of a specialized LD subpopulation characterized by a unique proteome and a defined geographical location at the nucleus-vacuole junction contact site. In search for factors determining identity of these LDs, we screened ∼6,000 yeast mutants for loss of targeting of the subpopulation marker Pdr16 and identified Ldo45 (LD organization protein of 45 kD) as a crucial targeting determinant. Ldo45 is the product of a splicing event connecting two adjacent genes (YMR147W and YMR148W/OSW5/LDO16). We show that Ldo proteins cooperate with the LD biogenesis component seipin and establish LD identity by defining positioning and surface-protein composition. Our studies suggest a mechanism to establish functional differentiation of organelles, opening the door to better understanding of metabolic decisions in cells.

Transcription rate strongly affects splicing fidelity and cotranscriptionality in budding yeast

The functional consequences of alternative splicing on altering the transcription rate have been the subject of intensive study in mammalian cells but less is known about effects of splicing on changing the transcription rate in yeast. We present several lines of evidence showing that slow RNA polymerase II elongation increases both cotranscriptional splicing and splicing efficiency and that faster elongation reduces cotranscriptional splicing and splicing efficiency in budding yeast, suggesting that splicing is more efficient when cotranscriptional. Moreover, we demonstrate that altering the RNA polymerase II elongation rate in either direction compromises splicing fidelity, and we reveal that splicing fidelity depends largely on intron length together with secondary structure and splice site score. These effects are notably stronger for the highly expressed ribosomal protein coding transcripts. We propose that transcription by RNA polymerase II is tuned to optimize the efficiency and accuracy of ribosomal protein gene expression, while allowing flexibility in splice site choice with the nonribosomal protein transcripts.

From transcriptional complexity to cellular phenotypes: Lessons from yeast

Pervasive transcription has been reported in many eukaryotic organisms, revealing a highly interleaved transcriptome organization that involves thousands of coding and non-coding RNAs. However, to date, the biological impact of transcriptome complexity is still poorly understood. Here I will review how subtle variations of the transcriptome can lead to divergent cellular phenotypes by fine-tuning both its coding potential and regulation. I will discuss strategies that can be used to link molecular variations with divergent biological outcomes. Finally, I will explore the implication of transcriptional complexity for our understanding of gene expression in the context of cell-to-cell phenotypic variability. Copyright © 2017 John Wiley & Sons, Ltd.

Nineteen complex-related factor Prp45 is required for the early stages of cotranscriptional spliceosome assembly

Splicing in S. cerevisiae has been shown to proceed cotranscriptionally, but the nature of the coupling remains a subject of debate. Here, we examine the effect of nineteen complex-related splicing factor Prp45 (a homolog of SNW1/SKIP) on cotranscriptional splicing. RNA-sequencing and RT-qPCR showed elevated pre-mRNA levels but only limited reduction of spliced mRNAs in cells expressing C-terminally truncated Prp45, Prp45(1-169). Assays with a series of reporters containing the AMA1 intron with regulatable splicing confirmed decreased splicing efficiency and showed the leakage of unspliced RNAs in prp45(1-169) cells. We also measured pre-mRNA accumulation of the meiotic MER2 gene, which depends on the expression of Mer1 factor for splicing. prp45(1-169) cells accumulated approximately threefold higher levels of MER2 pre-mRNA than WT cells only when splicing was induced. To monitor cotranscriptional splicing, we determined the presence of early spliceosome assembly factors and snRNP complexes along the ECM33 and ACT1 genes. We found that prp45(1-169) hampered the cotranscriptional recruitment of U2 and, to a larger extent, U5 and NTC, while the U1 profile was unaffected. The recruitment of Prp45(1-169) was impaired similarly to U5 snRNP and NTC. Our results imply that Prp45 is required for timely formation of complex A, prior to stable physical association of U5/NTC with the emerging pre-mRNA substrate. We suggest that Prp45 facilitates conformational rearrangements and/or contacts that couple U1 snRNP-recognition to downstream assembly events.

Intron retention-dependent gene regulation in Cryptococcus neoformans

The biological impact of alternative splicing is poorly understood in fungi, although recent studies have shown that these microorganisms are usually intron-rich. In this study, we re-annotated the genome of C. neoformans var. neoformans using RNA-Seq data. Comparison with C. neoformans var. grubii revealed that more than 99% of ORF-introns are in the same exact position in the two varieties whereas UTR-introns are much less evolutionary conserved. We also confirmed that alternative splicing is very common in C. neoformans, affecting nearly all expressed genes. We also observed specific regulation of alternative splicing by environmental cues in this yeast. However, alternative splicing does not appear to be an efficient method to diversify the C. neoformans proteome. Instead, our data suggest the existence of an intron retention-dependent mechanism of gene expression regulation that is not dependent on NMD. This regulatory process represents an additional layer of gene expression regulation in fungi and provides a mechanism to tune gene expression levels in response to any environmental modification.