Alternative splicing and subfunctionalization generates functional diversity in fungal proteomes

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human transcription factor (TF) genes encode multiple protein isoforms differing in DNA-binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators," both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Resolving discrepancies between chimeric and multiplicative measures of higher-order epistasis

Epistasis - the interaction between alleles at different genetic loci - plays a fundamental role in biology. However, several recent approaches quantify epistasis using a chimeric formula that measures deviations from a multiplicative fitness model on an additive scale, thus mixing two scales. Here, we show that for pairwise interactions, the chimeric formula yields a different magnitude but the same sign of epistasis compared to the multiplicative formula that measures both fitness and deviations on a multiplicative scale. However, for higher-order interactions, we show that the chimeric formula can have both different magnitude and sign compared to the multiplicative formula. We resolve these inconsistencies by deriving mathematical relationships between the different epistasis formulae and different parametrizations of the multivariate Bernoulli distribution. We argue that the chimeric formula does not appropriately model interactions between the Bernoulli random variables. In simulations, we show that the chimeric formula is less accurate than the classical multiplicative/additive epistasis formulae and may falsely detect higher-order epistasis. Analyzing multi-gene knockouts in yeast, multi-way drug interactions in E. coli, and deep mutational scanning of several proteins, we find that approximately 10% to 60% of inferred higher-order interactions change sign using the multiplicative/additive formula compared to the chimeric formula.

Current Technologies in Snake Venom Analysis and Applications

This comprehensive review explores the cutting-edge advancements in snake venom research, focusing on the integration of proteomics, genomics, transcriptomics, and bioinformatics. Highlighting the transformative impact of these technologies, the review delves into the genetic and ecological factors driving venom evolution, the complex molecular composition of venoms, and the regulatory mechanisms underlying toxin production. The application of synthetic biology and multi-omics approaches, collectively known as venomics, has revolutionized the field, providing deeper insights into venom function and its therapeutic potential. Despite significant progress, challenges such as the functional characterization of toxins and the development of cost-effective antivenoms remain. This review also discusses the future directions of venom research, emphasizing the need for interdisciplinary collaborations and new technologies (mRNAs, cryo-electron microscopy for structural determinations of toxin complexes, synthetic biology, and other technologies) to fully harness the biomedical potential of venoms and toxins from snakes and other animals.

Quantifying higher-order epistasis: beware the chimera

Epistasis, or interactions in which alleles at one locus modify the fitness effects of alleles at other loci, plays a fundamental role in genetics, protein evolution, and many other areas of biology. Epistasis is typically quantified by computing the deviation from the expected fitness under an additive or multiplicative model using one of several formulae. However, these formulae are not all equivalent. Importantly, one widely used formula - which we call the chimeric formula - measures deviations from a multiplicative fitness model on an additive scale, thus mixing two measurement scales. We show that for pairwise interactions, the chimeric formula yields a different magnitude, but the same sign (synergistic vs. antagonistic) of epistasis compared to the multiplicative formula that measures both fitness and deviations on a multiplicative scale. However, for higher-order interactions, we show that the chimeric formula can have both different magnitude and sign compared to the multiplicative formula - thus confusing negative epistatic interactions with positive interactions, and vice versa. We resolve these inconsistencies by deriving fundamental connections between the different epistasis formulae and the parameters of the multivariate Bernoulli distribution . Our results demonstrate that the additive and multiplicative epistasis formulae are more mathematically sound than the chimeric formula. Moreover, we demonstrate that the mathematical issues with the chimeric epistasis formula lead to markedly different biological interpretations of real data. Analyzing multi-gene knockout data in yeast, multi-way drug interactions in E. coli , and deep mutational scanning (DMS) of several proteins, we find that 10 - 60% of higher-order interactions have a change in sign with the multiplicative or additive epistasis formula. These sign changes result in qualitatively different findings on functional divergence in the yeast genome, synergistic vs. antagonistic drug interactions, and and epistasis between protein mutations. In particular, in the yeast data, the more appropriate multiplicative formula identifies nearly 500 additional negative three-way interactions, thus extending the trigenic interaction network by 25%.

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Functional Diversification and the Plant Secondary Cell Wall

Much evidence exists suggesting the presence of genetic functional diversification in plants, though literature associated with the role of functional diversification in the evolution of the plant secondary cell wall (SCW) has sparsely been compiled and reviewed in a recent context. This review aims to elucidate, through the examination of gene phylogenies associated with its biosynthesis and maintenance, the role of functional diversification in shaping the critical, dynamic, and characteristic organelle, the secondary cell wall. It will be asserted that gene families resulting from gene duplication and subsequent functional divergence are present and are heavily involved in SCW biosynthesis and maintenance. Furthermore, diversification will be presented as a significant driver behind the evolution of the many functional characteristics of the SCW. The structure and function of the plant cell wall and its constituents will first be explored, followed by a discussion on the phenomenon of gene duplication and the resulting genetic functional divergence that can emerge. Finally, the major constituents of the SCW and their individual relationships with duplication and divergence will be reviewed to the extent of current knowledge on the subject.

Concerted structural rearrangements enable RNA channeling into the cytoplasmic Ski238-Ski7-exosome assembly

The Ski2-Ski3-Ski8 (Ski238) helicase complex directs cytoplasmic mRNAs toward the nucleolytic exosome complex for degradation. In yeast, the interaction between Ski238 and exosome requires the adaptor protein Ski7. We determined different cryo-EM structures of the Ski238 complex depicting the transition from a rigid autoinhibited closed conformation to a flexible active open conformation in which the Ski2 helicase module has detached from the rest of Ski238. The open conformation favors the interaction of the Ski3 subunit with exosome-bound Ski7, leading to the recruitment of the exosome. In the Ski238-Ski7-exosome holocomplex, the Ski2 helicase module binds the exosome cap, enabling the RNA to traverse from the helicase through the internal exosome channel to the Rrp44 exoribonuclease. Our study pinpoints how conformational changes within the Ski238 complex regulate exosome recruitment for RNA degradation. We also reveal the remarkable conservation of helicase-exosome RNA channeling mechanisms throughout eukaryotic nuclear and cytoplasmic exosome complexes.

SKI complex: A multifaceted cytoplasmic RNA exosome cofactor in mRNA metabolism with links to disease, developmental processes, and antiviral responses

RNA stability and quality control are integral parts of gene expression regulation. A key factor shaping eukaryotic transcriptomes, mainly via 3'-5' exoribonucleolytic trimming or degradation of diverse transcripts in nuclear and cytoplasmic compartments, is the RNA exosome. Precise exosome targeting to various RNA molecules requires strict collaboration with specialized auxiliary factors, which facilitate interactions with its substrates. The predominant class of cytoplasmic RNA targeted by the exosome are protein-coding transcripts, which are carefully scrutinized for errors during translation. Normal, functional mRNAs are turned over following protein synthesis by the exosome or by Xrn1 5'-3'-exonuclease, acting in concert with Dcp1/2 decapping complex. In turn, aberrant transcripts are eliminated by dedicated surveillance pathways, triggered whenever ribosome translocation is impaired. Cytoplasmic 3'-5' mRNA decay and surveillance are dependent on the tight cooperation between the exosome and its evolutionary conserved co-factor-the SKI (superkiller) complex (SKIc). Here, we summarize recent findings from structural, biochemical, and functional studies of SKIc roles in controlling cytoplasmic RNA metabolism, including links to various cellular processes. Mechanism of SKIc action is illuminated by presentation of its spatial structure and details of its interactions with exosome and ribosome. Furthermore, contribution of SKIc and exosome to various mRNA decay pathways, usually converging on recycling of ribosomal subunits, is delineated. A crucial physiological role of SKIc is emphasized by describing association between its dysfunction and devastating human disease-a trichohepatoenteric syndrome (THES). Eventually, we discuss SKIc functions in the regulation of antiviral defense systems, cell signaling and developmental transitions, emerging from interdisciplinary investigations. This article is categorized under: RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms RNA Turnover and Surveillance > Regulation of RNA Stability RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

Alternative splicing diversifies the transcriptome and proteome of the rice blast fungus during host infection

Alternative splicing (AS) contributes to diversifying and regulating cellular responses to environmental conditions and developmental cues by differentially producing multiple mRNA and protein isoforms from a single gene. Previous studies on AS in pathogenic fungi focused on profiling AS isoforms under a limited number of conditions. We analysed AS profiles in the rice blast fungus Magnaporthe oryzae, a global threat to rice production, using high-quality transcriptome data representing its vegetative growth (mycelia) and multiple host infection stages. We identified 4,270 AS isoforms derived from 2,413 genes, including 499 genes presumably regulated by infection-specific AS. AS appears to increase during infection, with 32.7% of the AS isoforms being produced during infection but absent in mycelia. Analysis of the isoforms observed at each infection stage showed that 636 AS isoforms were more abundant than corresponding annotated mRNAs, especially after initial hyphal penetration into host cell. Many such dominant isoforms were predicted to encode regulatory proteins such as transcription factors and phospho-transferases. We also identified the genes encoding distinct proteins via AS and confirmed the translation of some isoforms via a proteomic analysis, suggesting potential AS-mediated neo-functionalization of some genes during infection. Comprehensive profiling of the pattern of genome-wide AS during multiple stages of rice-M. oryzae interaction established a foundational resource that will help investigate the role and regulation of AS during rice infection.

Retention of duplicated genes in evolution

Gene duplication is a prevalent phenomenon across the tree of life. The processes that lead to the retention of duplicated genes are not well understood. Functional genomics approaches in model organisms, such as yeast, provide useful tools to test the mechanisms underlying retention with functional redundancy and divergence of duplicated genes, including fates associated with neofunctionalization, subfunctionalization, back-up compensation, and dosage amplification. Duplicated genes may also be retained as a consequence of structural and functional entanglement. Advances in human gene editing have enabled the interrogation of duplicated genes in the human genome, providing new tools to evaluate the relative contributions of each of these factors to duplicate gene retention and the evolution of genome structure.

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.

Transcriptional Landscapes of Long Non-coding RNAs and Alternative Splicing in Pyricularia oryzae Revealed by RNA-Seq

Pyricularia oryzae causes the rice blast, which is one of the most devastating crop diseases worldwide, and is a model fungal pathogen widely used for dissecting the molecular mechanisms underlying fungal virulence/pathogenicity. Although the whole genome sequence of P. oryzae is publicly available, its current transcriptomes remain incomplete, lacking the information on non-protein coding genes and alternative splicing. Here, we performed and analyzed RNA-Seq of conidia and hyphae, resulting in the identification of 3,374 novel genes. Interestingly, the vast majority of these novel genes likely transcribed long non-coding RNAs (lncRNAs), and most of them were localized in the intergenic regions. Notably, their expressions were concomitant with the transcription of neighboring genes thereof in conidia and hyphae. In addition, 2,358 genes were found to undergo alternative splicing events. Furthermore, we exemplified that a lncRNA was important for hyphal growth likely by regulating the neighboring protein-coding gene and that alternative splicing of the transcription factor gene CON7 was required for appressorium formation. In summary, results from this study indicate that lncRNA transcripts and alternative splicing events are two important mechanisms for regulating the expression of genes important for conidiation, hyphal growth, and pathogenesis, and provide new insights into transcriptomes and gene regulation in the rice blast fungus.

Alternatively Splicing Interactomes Identify Novel Isoform-Specific Partners for NSD2

Nuclear receptor SET domain protein (NSD2) plays a fundamental role in the pathogenesis of Wolf-Hirschhorn Syndrome (WHS) and is overexpressed in multiple human myelomas, but its protein-protein interaction (PPI) patterns, particularly at the isoform/exon levels, are poorly understood. We explored the subcellular localizations of four representative NSD2 transcripts with immunofluorescence microscopy. Next, we used label-free quantification to perform immunoprecipitation mass spectrometry (IP-MS) analyses of the transcripts. Using the interaction partners for each transcript detected in the IP-MS results, we identified 890 isoform-specific PPI partners (83% are novel). These PPI networks were further divided into four categories of the exon-specific interactome. In these exon-specific PPI partners, two genes, RPL10 and HSPA8, were successfully confirmed by co-immunoprecipitation and Western blotting. RPL10 primarily interacted with Isoforms 1, 3, and 5, and HSPA8 interacted with all four isoforms, respectively. Using our extended NSD2 protein interactions, we constructed an isoform-level PPI landscape for NSD2 to serve as reference interactome data for NSD2 spliceosome-level studies. Furthermore, the RNA splicing processes supported by these isoform partners shed light on the diverse roles NSD2 plays in WHS and myeloma development. We also validated the interactions using Western blotting, RPL10, and the three NSD2 (Isoform 1, 3, and 5). Our results expand gene-level NSD2 PPI networks and provide a basis for the treatment of NSD2-related developmental diseases.

Zebrafish Ski7 tunes RNA levels during the oocyte-to-embryo transition

Post-transcriptional regulation of gene expression is crucial during the oocyte-to-embryo transition, a highly dynamic process characterized by the absence of nuclear transcription. Thus, changes to the RNA content are solely dependent on RNA degradation. Although several mechanisms that promote RNA decay during embryogenesis have been identified, it remains unclear which machineries contribute to remodeling the maternal transcriptome. Here, we focused on the degradation factor Ski7 in zebrafish. Homozygous ski7 mutant fish had higher proportions of both poor quality eggs and eggs that were unable to develop beyond the one-cell stage. Consistent with the idea that Ski7 participates in remodeling the maternal RNA content, transcriptome profiling identified hundreds of misregulated mRNAs in the absence of Ski7. Furthermore, upregulated genes were generally lowly expressed in wild type, suggesting that Ski7 maintains low transcript levels for this subset of genes. Finally, GO enrichment and proteomic analyses of misregulated factors implicated Ski7 in the regulation of redox processes. This was confirmed experimentally by an increased resistance of ski7 mutant embryos to reductive stress. Our results provide first insights into the physiological role of vertebrate Ski7 as a post-transcriptional regulator during the oocyte-to-embryo transition.

Mitochondrial protein phosphorylation in yeast revisited

Protein phosphorylation is one of the best-known post-translational modifications occurring in all domains of life. In eukaryotes, protein phosphorylation affects all cellular compartments including mitochondria. High-throughput techniques of mass spectrometry combined with cell fractionation and biochemical methods yielded thousands of phospho-sites on hundreds of mitochondrial proteins. We have compiled the information on mitochondrial protein kinases and phosphatases and their substrates in Saccharomyces cerevisiae and provide the current state-of-the-art overview of mitochondrial protein phosphorylation in this model eukaryote. Using several examples, we describe emerging features of the yeast mitochondrial phosphoproteome and present challenges lying ahead in this exciting field.

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.

Recurrent sequence evolution after independent gene duplication

Background

Convergent and parallel evolution provide unique insights into the mechanisms of natural selection. Some of the most striking convergent and parallel (collectively recurrent) amino acid substitutions in proteins are adaptive, but there are also many that are selectively neutral. Accordingly, genome-wide assessment has shown that recurrent sequence evolution in orthologs is chiefly explained by nearly neutral evolution. For paralogs, more frequent functional change is expected because additional copies are generally not retained if they do not acquire their own niche. Yet, it is unknown to what extent recurrent sequence differentiation is discernible after independent gene duplications in different eukaryotic taxa.

Results

We develop a framework that detects patterns of recurrent sequence evolution in duplicated genes. This is used to analyze the genomes of 90 diverse eukaryotes. We find a remarkable number of families with a potentially predictable functional differentiation following gene duplication. In some protein families, more than ten independent duplications show a similar sequence-level differentiation between paralogs. Based on further analysis, the sequence divergence is found to be generally asymmetric. Moreover, about 6% of the recurrent sequence evolution between paralog pairs can be attributed to recurrent differentiation of subcellular localization. Finally, we reveal the specific recurrent patterns for the gene families Hint1/Hint2, Sco1/Sco2 and vma11/vma3.

Conclusions

The presented methodology provides a means to study the biochemical underpinning of functional differentiation between paralogs. For instance, two abundantly repeated substitutions are identified between independently derived Sco1 and Sco2 paralogs. Such identified substitutions allow direct experimental testing of the biological role of these residues for the repeated functional differentiation. We also uncover a diverse set of families with recurrent sequence evolution and reveal trends in the functional and evolutionary trajectories of this hitherto understudied phenomenon.

Exploring whole-genome duplicate gene retention with complex genetic interaction analysis

Whole-genome duplication has played a central role in the genome evolution of many organisms, including the human genome. Most duplicated genes are eliminated, and factors that influence the retention of persisting duplicates remain poorly understood. We describe a systematic complex genetic interaction analysis with yeast paralogs derived from the whole-genome duplication event. Mapping of digenic interactions for a deletion mutant of each paralog, and of trigenic interactions for the double mutant, provides insight into their roles and a quantitative measure of their functional redundancy. Trigenic interaction analysis distinguishes two classes of paralogs: a more functionally divergent subset and another that retained more functional overlap. Gene feature analysis and modeling suggest that evolutionary trajectories of duplicated genes are dictated by combined functional and structural entanglement factors.

RNA Exosomes and Their Cofactors

The RNA exosome is a ribonucleolytic multiprotein complex that is conserved and essential in all eukaryotes. Although we tend to speak of "the" exosome complex, it should be more correctly viewed as several different subtypes that share a common core. Subtypes of the exosome complex are present in the cytoplasm, the nucleus and the nucleolus of all eukaryotic cells, and carry out the 3'-5' processing and/or degradation of a wide range of RNA substrates.Because the substrate specificity of the exosome complex is determined by cofactors, the system is highly adaptable, and different organisms have adjusted the machinery to their specific needs. Here, we present an overview of exosome complexes and their cofactors that have been described in different eukaryotes.

Genome-wide evolution of wobble base-pairing nucleotides of branchpoint motifs with increasing organismal complexity

How have the branchpoint motifs evolved in organisms of different complexity? Here we identified and examined the consensus motifs (R₁C₂T₃R₄A₅Y₆, R: A or G, Y: C or T) of 898 fungal genomes. In Ascomycota unicellular yeasts, the G₄/A₄ ratio is mostly (98%) below 0.125 but increases sharply in multicellular species by about 40 times on average, and in the more complex Basidiomycota, it increases further by about 7 times. The global G₄ increase is consistent with A₄ to G₄ transitions in evolution. Of the G₄/A₄-interacting amino acids of the branchpoint binding protein MSL5 (SF1) and the HSH155 (SF3B1), as well as the 5' splice sites (SS) and U2 snRNA genes, the 5' SS G³/A³ co-vary with the G₄ to some extent. However, corresponding increase of the G₄-complementary GCAGTA-U2 gene is rare, suggesting wobble-base pairing between the G₄-containing branchpoint motif and GTAGTA-U2 in most of these species. Interestingly, the G₄/A₄ ratio correlates well with the abundance of alternative splicing in the two phyla, and G₄ enriched significantly at the alternative 3' SS of genes in RNA metabolism, kinases and membrane proteins. Similar wobble nucleotides also enriched at the 3' SS of multicellular fungi with only thousands of protein-coding genes. Thus, branchpoint motifs have evolved U2-complementarity in unicellular Ascomycota yeasts, but have gradually gained more wobble base-pairing nucleotides in fungi of higher complexity, likely to destabilize branchpoint motif-U2 interaction and/or branchpoint A protrusion for alternative splicing. This implies an important role of relaxing the branchpoint signals in the multicellularity and further complexity of fungi.

Molecular evolution and lineage-specific expansion of the PP2C family in Zea mays

97 ZmPP2Cs were clustered into 10 subfamilies with biased subfamily evolution and lineage-specific expansion. Segmental duplication after the divergence of maize and sorghum might have led to primary expansion of ZmPP2Cs. The protein phosphatase 2C (PP2C) enzymes control many stress responses and developmental processes in plants. In Zea mays, a comprehensive understanding of the evolution and expansion of the PP2C family is still lacking. In the current study, 97 ZmPP2Cs were identified and clustered into 10 subfamilies. Through the analysis of the PP2C family in monocots, the ZmPP2C subfamilies displayed biased subfamily molecular evolution and lineage-specific expansion, as evidenced by their differing numbers of member genes, expansion and evolutionary rates, conserved subdomains, chromosomal distributions, expression levels, responsive-regulatory elements and regulatory networks. Moreover, while segmental duplication events have caused the primary expansion of the ZmPP2Cs, the majority of their diversification occurred following the additional whole-genome duplication that took place after the divergence of maize and sorghum (Sorghum bicolor). After this event, the PP2C subfamilies showed asymmetric evolutionary rates, with the D, F₂ and H subfamily likely the most closely to resemble its ancestral subfamily's genes. These findings could provide novel insights into the molecular evolution and expansion of the PP2C family in maize, and lay the foundation for the functional analysis of these enzymes in maize and related monocots.

Regulation of the Hsf1-dependent transcriptome via conserved bipartite contacts with Hsp70 promotes survival in yeast

Protein homeostasis and cellular fitness in the presence of proteotoxic stress is promoted by heat shock factor 1 (Hsf1), which controls basal and stress-induced expression of molecular chaperones and other targets. The major heat shock proteins and molecular chaperones Hsp70 and Hsp90, in turn, participate in a negative feedback loop that ensures appropriate coordination of the heat shock response with environmental conditions. Features of this regulatory circuit in the budding yeast Saccharomyces cerevisiae have been recently defined, most notably regarding direct interaction between Hsf1 and the constitutively expressed Hsp70 protein Ssa1. Here, we sought to further examine the Ssa1/Hsf1 regulation. We found that Ssa1 interacts independently with both the previously defined CE2 site in the Hsf1 C-terminal transcriptional activation domain and with an additional site that we identified within the N-terminal activation domain. Consistent with both sites bearing a recognition signature for Hsp70, we demonstrate that Ssa1 contacts Hsf1 via its substrate-binding domain and that abolishing either regulatory site results in loss of Ssa1 interaction. Removing Hsp70 regulation of Hsf1 globally dysregulated Hsf1 transcriptional activity, with synergistic effects on both gene expression and cellular fitness when both sites are disrupted together. Finally, we report that Hsp70 interacts with both transcriptional activation domains of Hsf1 in the related yeast Lachancea kluyveri Our findings indicate that Hsf1 transcriptional activity is tightly regulated to ensure cellular fitness and that a general and conserved Hsp70-HSF1 feedback loop regulates cellular proteostasis in yeast.

Ser/Thr protein phosphatases in fungi: structure, regulation and function

Reversible phospho-dephosphorylation of proteins is a major mechanism for the control of cellular functions. By large, Ser and Thr are the most frequently residues phosphorylated in eukar-yotes. Removal of phosphate from these amino acids is catalyzed by a large family of well-conserved enzymes, collectively called Ser/Thr protein phosphatases. The activity of these enzymes has an enormous impact on cellular functioning. In this work we pre-sent the members of this family in S. cerevisiae and other fungal species, and review the most recent findings concerning their regu-lation and the roles they play in the most diverse aspects of cell biology.

A short splicing isoform of HBS1L links the cytoplasmic exosome and SKI complexes in humans

The exosome complex is a major eukaryotic exoribonuclease that requires the SKI complex for its activity in the cytoplasm. In yeast, the Ski7 protein links both complexes, whereas a functional equivalent of the Ski7 has remained unknown in the human genome.

Proteomic analysis revealed that a previously uncharacterized short splicing isoform of HBS1L (HBS1LV3) is the long-sought factor linking the exosome and SKI complexes in humans. In contrast, the canonical HBS1L variant, HBS1LV1, which acts as a ribosome dissociation factor, does not associate with the exosome and instead interacts with the mRNA surveillance factor PELOTA. Interestingly, both HBS1LV1 and HBS1LV3 interact with the SKI complex and HBS1LV1 seems to antagonize SKI/exosome supercomplex formation. HBS1LV3 contains a unique C-terminal region of unknown structure, with a conserved RxxxFxxxL motif responsible for exosome binding and may interact with the exosome core subunit RRP43 in a way that resembles the association between Rrp6 RNase and Rrp43 in yeast. HBS1LV3 or the SKI complex helicase (SKI2W) depletion similarly affected the transcriptome, deregulating multiple genes. Furthermore, half-lives of representative upregulated mRNAs were increased, supporting the involvement of HBS1LV3 and SKI2W in the same mRNA degradation pathway, essential for transcriptome homeostasis in the cytoplasm.

A Two-Headed Monster to Avert Disaster: HBS1/SKI7 Is Alternatively Spliced to Build Eukaryotic RNA Surveillance Complexes

The cytosolic RNA exosome, a 3'→5' exoribonuclease complex, contributes to mRNA degradation in eukaryotes, limiting the accumulation of poorly-translated, improperly translated, or aberrant mRNA species. Disruption of cytosolic RNA exosome activity allows aberrant RNA species to accumulate, which can then be detected by host antiviral immune systems as a signature of pathogen infection, activating antiviral defenses. SKI7 is a critical component of the cytosolic RNA exosome in yeast, bridging the catalytic exoribonuclease core with the SKI2/SKI3/SKI8 adaptor complex that guides aberrant RNA substrates into the exosome. The ortholog of SKI7 was only recently identified in humans as an alternative splice form of the HBS1 gene, which encodes a decoding factor translational GTPase that rescues stalled ribosomes. Here, we identify the plant orthologs of HBS1/SKI7. We found that HBS1 and SKI7 are typically encoded by alternative splice forms of a single locus, although some plant lineages have evolved subfunctionalized genes that apparently encode only HBS1 or only SKI7. In all plant lineages examined, the SKI7 gene is subject to regulation by alternative splicing that can yield unproductive transcripts, either by removing deeply conserved SKI7 coding sequences, or by introducing premature stop codons that render SKI7 susceptible to nonsense-mediated decay. Taking a comparative, evolutionary approach, we define crucial features of the SKI7 protein shared by all eukaryotes, and use these deeply conserved features to identify SKI7 proteins in invertebrate lineages. We conclude that SKI7 is a conserved cytosolic RNA exosome subunit across eukaryotic lineages, and that SKI7 is consistently regulated by alternative splicing, suggesting broad coordination of nuclear and cytosolic RNA metabolism.

Conservation of mRNA quality control factor Ski7 and its diversification through changes in alternative splicing and gene duplication

Eukaryotes maintain fidelity of gene expression by preferential degradation of aberrant mRNAs that arise by errors in RNA processing reactions. In Saccharomyces cerevisiae, Ski7 plays an important role in this mRNA quality control by mediating mRNA degradation by the RNA exosome. Ski7 was initially thought to be restricted to Saccharomyces cerevisiae and close relatives because the SKI7 gene and its paralog HBS1 arose by whole genome duplication (WGD) in a recent ancestor. We have recently shown that the preduplication gene was alternatively spliced and that Ski7 function predates WGD. Here, we use transcriptome analysis of diverse eukaryotes to show that diverse eukaryotes use alternative splicing of SKI7/HBS1 to encode two proteins. Although alternative splicing affects the same intrinsically disordered region of the protein, the pattern of splice site usage varies. This alternative splicing event arose in an early eukaryote that is a common ancestor of plants, animals, and fungi. Remarkably, through changes in alternative splicing and gene duplication, the Ski7 protein has diversified such that different species express one of four distinct Ski7-like proteins. We also show experimentally that the Saccharomyces cerevisiae SKI7 gene has undergone multiple changes that are incompatible with the Hbs1 function and may also have undergone additional changes to optimize mRNA quality control. The combination of transcriptome analysis in diverse eukaryotes and genetic analysis in yeast clarifies the mechanism by which a Ski7-like protein is expressed across eukaryotes and provides a unique view of changes in alternative splicing patterns of one gene over long evolutionary time.

The nonstop decay and the RNA silencing systems operate cooperatively in plants

Translation-dependent mRNA quality control systems protect the protein homeostasis of eukaryotic cells by eliminating aberrant transcripts and stimulating the decay of their protein products. Although these systems are intensively studied in animals, little is known about the translation-dependent quality control systems in plants. Here, we characterize the mechanism of nonstop decay (NSD) system in Nicotiana benthamiana model plant. We show that plant NSD efficiently degrades nonstop mRNAs, which can be generated by premature polyadenylation, and stop codon-less transcripts, which are produced by endonucleolytic cleavage. We demonstrate that in plants, like in animals, Pelota, Hbs1 and SKI2 proteins are required for NSD, supporting that NSD is an ancient and conserved eukaryotic quality control system. Relevantly, we found that NSD and RNA silencing systems cooperate in plants. Plant silencing predominantly represses target mRNAs through endonucleolytic cleavage in the coding region. Here we show that NSD is required for the elimination of 5' cleavage product of mi- or siRNA-guided silencing complex when the cleavage occurs in the coding region. We also show that NSD and nonsense-mediated decay (NMD) quality control systems operate independently in plants.

Alternative splice variants of rhomboid proteins: Comparative analysis of database entries for select model organisms and validation of functional potential

Background: Rhomboid serine proteases are present across many species and are often encoded in each species by more than one predicted gene. Based on protein sequence comparisons, rhomboids can be differentiated into groups - secretases, presenilin-like associated rhomboid-like (PARL) proteases, iRhoms, and "inactive" rhomboid proteins. Although these rhomboid groups are distinct, the different types can operate simultaneously. Studies in Arabidopsis showed that the number of rhomboid proteins working simultaneously can be further diversified by alternative splicing. This phenomenon was confirmed for the Arabidopsis plastid rhomboid proteins At1g25290 and At1g74130. Although alternative splicing was determined to be a significant mechanism for diversifying these two Arabidopsis plastid rhomboids, there has yet to be an assessment as to whether this mechanism extends to other rhomboids and to other species.  Methods: We thus conducted a comparative analysis of select databases to determine if the alternative splicing mechanism observed for the two Arabidopsis plastid rhomboids was utilized in other species to expand the repertoire of rhomboid proteins. To help verify the in silico observations, select splice variants from different groups were tested for activity using transgenic- and additive-based assays. These assays aimed to uncover evidence that the selected splice variants display capacities to influence processes like antimicrobial sensitivity. Results: A comparison of database entries of six widely used eukaryotic experimental models  (human, mouse, Arabidopsis, Drosophila, nematode, and yeast) revealed robust usage of alternative splicing to diversify rhomboid protein structure across the various motifs or regions, especially in human, mouse and Arabidopsis. Subsequent validation studies uncover evidence that the splice variants selected for testing displayed functionality in the different activity assays. Conclusions: The combined results support the hypothesis that alternative splicing is likely used to diversify and expand rhomboid protein functionality, and this potentially occurred across the various motifs or regions of the protein.

Alternative splicing and the evolution of phenotypic novelty

Alternative splicing, a mechanism of post-transcriptional RNA processing whereby a single gene can encode multiple distinct transcripts, has been proposed to underlie morphological innovations in multicellular organisms. Genes with developmental functions are enriched for alternative splicing events, suggestive of a contribution of alternative splicing to developmental programmes. The role of alternative splicing as a source of transcript diversification has previously been compared to that of gene duplication, with the relationship between the two extensively explored. Alternative splicing is reduced following gene duplication with the retention of duplicate copies higher for genes which were alternatively spliced prior to duplication. Furthermore, and unlike the case for overall gene number, the proportion of alternatively spliced genes has also increased in line with the evolutionary diversification of cell types, suggesting alternative splicing may contribute to the complexity of developmental programmes. Together these observations suggest a prominent role for alternative splicing as a source of functional innovation. However, it is unknown whether the proliferation of alternative splicing events indeed reflects a functional expansion of the transcriptome or instead results from weaker selection acting on larger species, which tend to have a higher number of cell types and lower population sizes.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

The non-stop decay mRNA surveillance pathway is required for oxidative stress tolerance

Reactive oxygen species (ROS) are toxic by-products of normal aerobic metabolism. ROS can damage mRNAs and the translational apparatus resulting in translational defects and aberrant protein production. Three mRNA quality control systems monitor mRNAs for translational errors: nonsense-mediated decay, non-stop decay (NSD) and no-go decay (NGD) pathways. Here, we show that factors required for the recognition of NSD substrates and components of the SKI complex are required for oxidant tolerance. We found an overlapping requirement for Ski7, which bridges the interaction between the SKI complex and the exosome, and NGD components (Dom34/Hbs1) which have been shown to function in both NSD and NGD. We show that ski7 dom34 and ski7 hbs1 mutants are sensitive to hydrogen peroxide stress and accumulate an NSD substrate. We further show that NSD substrates are generated during ROS exposure as a result of aggregation of the Sup35 translation termination factor, which increases stop codon read-through allowing ribosomes to translate into the 3΄-end of mRNAs. Overexpression of Sup35 decreases stop codon read-through and rescues oxidant tolerance consistent with this model. Our data reveal an unanticipated requirement for the NSD pathway during oxidative stress conditions which prevents the production of aberrant proteins from NSD mRNAs.

Chromatin-remodeling SWI/SNF complex regulates coenzyme Q6 synthesis and a metabolic shift to respiration in yeast

Despite its relatively streamlined genome, there are many important examples of regulated RNA splicing in Saccharomyces cerevisiae Here, we report a role for the chromatin remodeler SWI/SNF in respiration, partially via the regulation of splicing. We find that a nutrient-dependent decrease in Snf2 leads to an increase in splicing of the PTC7 transcript. The spliced PTC7 transcript encodes a mitochondrial phosphatase regulator of biosynthesis of coenzyme Q₆ (ubiquinone or CoQ₆) and a mitochondrial redox-active lipid essential for electron and proton transport in respiration. Increased splicing of PTC7 increases CoQ₆ levels. The increase in PTC7 splicing occurs at least in part due to down-regulation of ribosomal protein gene expression, leading to the redistribution of spliceosomes from this abundant class of intron-containing RNAs to otherwise poorly spliced transcripts. In contrast, a protein encoded by the nonspliced isoform of PTC7 represses CoQ₆ biosynthesis. Taken together, these findings uncover a link between Snf2 expression and the splicing of PTC7 and establish a previously unknown role for the SWI/SNF complex in the transition of yeast cells from fermentative to respiratory modes of metabolism.

Translation of poly(A) tails leads to precise mRNA cleavage

Translation of poly(A) tails leads to mRNA cleavage but the mechanism and global pervasiveness of this "nonstop/no-go" decay process is not understood. Here we performed ribosome profiling (in a yeast strain lacking exosome function) of short 15-18 nucleotides mRNA footprints to identify ribosomes stalled at 3' ends of mRNA decay intermediates. In this background, we found mRNA cleavage extending hundreds of nucleotides upstream of ribosome stalling in poly(A) and predominantly in one reading frame. These observations suggest that decay-triggering endonucleolytic cleavage is closely associated with the ribosome. Surprisingly, ribosomes appeared to accumulate (i.e., stall) in the transcriptome when as few as three consecutive ORF-internal lysine codons were positioned in the A, P, and E sites though significant mRNA degradation was not observed. Endonucleolytic cleavage was found, however, at sites of premature polyadenylation (encoding polylysine) and rescue of the ribosomes stalled at these sites was dependent on Dom34. These results suggest this process may be critical when changes in the polyadenylation site occur during development, tumorigenesis, or when translation termination/recycling is impaired.

The Evolutionary Relationship between Alternative Splicing and Gene Duplication

The protein diversity that exists today has resulted from various evolutionary processes. It is well known that gene duplication (GD) along with the accumulation of mutations are responsible, among other factors, for an increase in the number of different proteins. The gene structure in eukaryotes requires the removal of non-coding sequences, introns, to produce mature mRNAs. This process, known as cis-splicing, referred to here as splicing, is regulated by several factors which can lead to numerous splicing arrangements, commonly designated as alternative splicing (AS). AS, producing several transcripts isoforms form a single gene, also increases the protein diversity. However, the evolution and manner for increasing protein variation differs between AS and GD. An important question is how are patterns of AS affected after a GD event. Here, we review the current knowledge of AS and GD, focusing on their evolutionary relationship. These two processes are now considered the main contributors to the increasing protein diversity and therefore their relationship is a relevant, yet understudied, area of evolutionary study.

Plant Actin-Depolymerizing Factors Possess Opposing Biochemical Properties Arising from Key Amino Acid Changes throughout Evolution

Functional divergence in paralogs is an important genetic source of evolutionary innovation. Actin-depolymerizing factors (ADFs) are among the most important actin binding proteins and are involved in generating and remodeling actin cytoskeletal architecture via their conserved F-actin severing or depolymerizing activity. In plants, ADFs coevolved with actin, but their biochemical properties are diverse. Unfortunately, the biochemical function of most plant ADFs and the potential mechanisms of their functional divergence remain unclear. Here, in vitro biochemical analyses demonstrated that all 11 ADF genes in Arabidopsis thaliana exhibit opposing biochemical properties. Subclass III ADFs evolved F-actin bundling (B-type) function from conserved F-actin depolymerizing (D-type) function, and subclass I ADFs have enhanced D-type function. By tracking historical mutation sites on ancestral proteins, several fundamental amino acid residues affecting the biochemical functions of these proteins were identified in Arabidopsis and various plants, suggesting that the biochemical divergence of ADFs has been conserved during the evolution of angiosperm plants. Importantly, N-terminal extensions on subclass III ADFs that arose from intron-sliding events are indispensable for the alteration of D-type to B-type function. We conclude that the evolution of these N-terminal extensions and several conserved mutations produced the diverse biochemical functions of plant ADFs from a putative ancestor.

Cytoplasmic RNA decay pathways - Enzymes and mechanisms

RNA decay plays a crucial role in post-transcriptional regulation of gene expression. Work conducted over the last decades has defined the major mRNA decay pathways, as well as enzymes and their cofactors responsible for these processes. In contrast, our knowledge of the mechanisms of degradation of non-protein coding RNA species is more fragmentary. This review is focused on the cytoplasmic pathways of mRNA and ncRNA degradation in eukaryotes. The major 3' to 5' and 5' to 3' mRNA decay pathways are described with emphasis on the mechanisms of their activation by the deprotection of RNA ends. More recently discovered 3'-end modifications such as uridylation, and their relevance to cytoplasmic mRNA decay in various model organisms, are also discussed. Finally, we provide up-to-date findings concerning various pathways of non-coding RNA decay in the cytoplasm.

The power of fission: yeast as a tool for understanding complex splicing

Pre-mRNA splicing is an essential component of eukaryotic gene expression. Many metazoans, including humans, regulate alternative splicing patterns to generate expansions of their proteome from a limited number of genes. Importantly, a considerable fraction of human disease causing mutations manifest themselves through altering the sequences that shape the splicing patterns of genes. Thus, understanding the mechanistic bases of this complex pathway will be an essential component of combating these diseases. Dating almost to the initial discovery of splicing, researchers have taken advantage of the genetic tractability of budding yeast to identify the components and decipher the mechanisms of splicing. However, budding yeast lacks the complex splicing machinery and alternative splicing patterns most relevant to humans. More recently, many researchers have turned their efforts to study the fission yeast, Schizosaccharomyces pombe, which has retained many features of complex splicing, including degenerate splice site sequences, the usage of exonic splicing enhancers, and SR proteins. Here, we review recent work using fission yeast genetics to examine pre-mRNA splicing, highlighting its promise for modeling the complex splicing seen in higher eukaryotes.

Coregulation of tandem duplicate genes slows evolution of subfunctionalization in mammals

Gene duplication is a fundamental process in genome evolution. However, most young duplicates are degraded by loss-of-function mutations, and the factors that allow some duplicate pairs to survive long-term remain controversial. One class of models to explain duplicate retention invokes sub- or neofunctionalization, whereas others focus on sharing of gene dosage. RNA-sequencing data from 46 human and 26 mouse tissues indicate that subfunctionalization of expression evolves slowly and is rare among duplicates that arose within the placental mammals, possibly because tandem duplicates are coregulated by shared genomic elements. Instead, consistent with the dosage-sharing hypothesis, most young duplicates are down-regulated to match expression levels of single-copy genes. Thus, dosage sharing of expression allows for the initial survival of mammalian duplicates, followed by slower functional adaptation enabling long-term preservation.

Mutations in the G-domain of Ski7 cause specific dysfunction in non-stop decay

Ski7 functions as a cofactor in both normal mRNA turnover and non-stop mRNA decay (NSD) mRNA surveillance in budding yeast. The N-terminal region of Ski7 (Ski7_N) interacts with the ski-complex and the exosome. The C-terminal region of Ski7 (Ski7_C) binds guanine nucleotides and shares overall sequence and structural homology with the proteins of the translational GTPase superfamily, especially the tRNA/tRNA-mimic carrier protein subfamilies such as EF1α, eRF3, and Hbs1. Previous reports showed that Ski7_N polypeptide functions adequately in vivo, while Ski7_C, if any, only slightly. Furthermore, Ski7_C does not exhibit GTP-hydrolysing activities under normal conditions. Therefore, the physiological and functional significance of the conserved Ski7_C is unclear. Here, we report strong genetic evidence suggesting differential roles for Ski7_N and Ski7_C in normal and specific mRNA turnover pathways by creating/isolating mutations in both Ski7_N and Ski7_C conserved motifs using indicator yeast strains. We concluded that Ski7_C participates in mRNA surveillance as a regulatory module competitively with the Hbs1/Dom34 complex. Our results provide insights into the molecular regulatory mechanisms underlying mRNA surveillance.

Structure of a Cytoplasmic 11-Subunit RNA Exosome Complex

The RNA exosome complex associates with nuclear and cytoplasmic cofactors to mediate the decay, surveillance, or processing of a wide variety of transcripts. In the cytoplasm, the conserved core of the exosome (Exo10) functions together with the conserved Ski complex. The interaction of S. cerevisiae Exo10 and Ski is not direct but requires a bridging cofactor, Ski7. Here, we report the 2.65 Å resolution structure of S. cerevisiae Exo10 bound to the interacting domain of Ski7. Extensive hydrophobic interactions rationalize the high affinity and stability of this complex, pointing to Ski7 as a constitutive component of the cytosolic exosome. Despite the absence of sequence homology, cytoplasmic Ski7 and nuclear Rrp6 bind Exo10 using similar surfaces and recognition motifs. Knowledge of the interacting residues in the yeast complexes allowed us to identify a splice variant of human HBS1-Like as a Ski7-like exosome-binding protein, revealing the evolutionary conservation of this cytoplasmic cofactor.

Selection for reduced translation costs at the intronic 5' end in fungi

It is generally believed that introns are not translated; therefore, the potential intronic features that may be related to the translation step (occurring after splicing) have yet to be thoroughly studied. Here, focusing on four fungi, we performed for the first time a comprehensive study aimed at characterizing how translation efficiency is encoded in introns and affects their evolution. By analysing their intronome we provide evidence of selection for STOP codons close to the intronic 5′ end, and show that the beginning of introns are selected for significantly high translation, presumably to reduce translation and metabolic costs in cases of non-spliced introns. Ribosomal profiling data analysis in Saccharomyces cerevisiae supports the conjecture that in this organism intron retention frequently occurs, introns are partially translated, and their translation efficiency affects organismal fitness. We show that the reported results are more significant in highly translated and highly spliced genes, but are not associated only with genes with a specific function. We also discuss the potential relation of the reported signals to efficient nonsense-mediated decay due to splicing errors. These new discoveries are supported by population-genetics considerations. In addition, they are contributory steps towards a broader understanding of intron evolution and the effect of silent mutations on gene expression and organismal fitness.

Cytosolic splice isoform of Hsp70 nucleotide exchange factor Fes1 is required for the degradation of misfolded proteins in yeast

Yeast Hsp70 nucleotide exchange factor Fes1 is expressed by rare alternative splicing as two isoforms. Fes1L is targeted to the nucleus, and Fes1S localizes to the cytosol and is required for the efficient proteasomal degradation of cytosolic misfolded proteins, as well as of species that are imported into the nucleus for degradation.

Cells maintain proteostasis by selectively recognizing and targeting misfolded proteins for degradation. In Saccharomyces cerevisiae, the Hsp70 nucleotide exchange factor Fes1 is essential for the degradation of chaperone-associated misfolded proteins by the ubiquitin-proteasome system. Here we show that the FES1 transcript undergoes unique 3′ alternative splicing that results in two equally active isoforms with alternative C-termini, Fes1L and Fes1S. Fes1L is actively targeted to the nucleus and represents the first identified nuclear Hsp70 nucleotide exchange factor. In contrast, Fes1S localizes to the cytosol and is essential to maintain proteostasis. In the absence of Fes1S, the heat-shock response is constitutively induced at normally nonstressful conditions. Moreover, cells display severe growth defects when elevated temperatures, amino acid analogues, or the ectopic expression of misfolded proteins, induce protein misfolding. Importantly, misfolded proteins are not targeted for degradation by the ubiquitin-proteasome system. These observations support the notion that cytosolic Fes1S maintains proteostasis by supporting the removal of toxic misfolded proteins by proteasomal degradation. This study provides key findings for the understanding of the organization of protein quality control mechanisms in the cytosol and nucleus.

Nucleotide sequence composition adjacent to intronic splice sites improves splicing efficiency via its effect on pre-mRNA local folding in fungi

RNA splicing is the central process of intron removal in eukaryotes known to regulate various cellular functions such as growth, development, and response to external signals. The canonical sequences indicating the splicing sites needed for intronic boundary recognition are well known. However, the roles and evolution of the local folding of intronic and exonic sequence features adjacent to splice sites has yet to be thoroughly studied. Here, focusing on four fungi (Saccharomyces cerevisiae, Schizosaccharomyces pombe, Aspergillus nidulans, and Candida albicans), we performed for the first time a comprehensive high-resolution study aimed at characterizing the encoding of intronic splicing efficiency in pre-mRNA transcripts and its effect on intron evolution. Our analysis supports the conjecture that pre-mRNA local folding strength at intronic boundaries is under selective pressure, as it significantly affects splicing efficiency. Specifically, we show that in the immediate region of 12-30 nucleotides (nt) surrounding the intronic donor site there is a preference for weak pre-mRNA folding; similarly, in the region of 15-33 nt surrounding the acceptor and branch sites there is a preference for weak pre-mRNA folding. We also show that in most cases there is a preference for strong pre-mRNA folding further away from intronic splice sites. In addition, we demonstrate that these signals are not associated with gene-specific functions, and they correlate with splicing efficiency measurements (r = 0.77, P = 2.98 × 10(-21)) and with expression levels of the corresponding genes (P = 1.24 × 10(-19)). We suggest that pre-mRNA folding strength in the above-mentioned regions has a direct effect on splicing efficiency by improving the recognition of intronic boundaries. These new discoveries are contributory steps toward a broader understanding of splicing regulation and intronic/transcript evolution.

Proteins involved in the degradation of cytoplasmic mRNA in the major eukaryotic model systems

The process of mRNA decay and surveillance is considered to be one of the main posttranscriptional gene expression regulation platforms in eukaryotes. The degradation of stable, protein-coding transcripts is normally initiated by removal of the poly(A) tail followed by 5'-cap hydrolysis and degradation of the remaining mRNA body by Xrn1. Alternatively, the exosome complex degrades mRNA in the 3'>5'direction. The newly discovered uridinylation-dependent pathway, which is present in many different organisms, also seems to play a role in bulk mRNA degradation. Simultaneously, to avoid the synthesis of incorrect proteins, special cellular machinery is responsible for the removal of faulty transcripts via nonsense-mediated, no-go, non-stop or non-functional 18S rRNA decay. This review is focused on the major eukaryotic cytoplasmic mRNA degradation pathways showing many similarities and pointing out main differences between the main model-species: yeast, Drosophila, plants and mammals.

Origin of the Yeast Whole-Genome Duplication

Whole-genome duplications (WGDs) are rare evolutionary events with profound consequences. They double an organism’s genetic content, immediately creating a reproductive barrier between it and its ancestors and providing raw material for the divergence of gene functions between paralogs. Almost all eukaryotic genome sequences bear evidence of ancient WGDs, but the causes of these events and the timing of intermediate steps have been difficult to discern. One of the best-characterized WGDs occurred in the lineage leading to the baker’s yeast Saccharomyces cerevisiae. Marcet-Houben and Gabaldón now show that, rather than simply doubling the DNA of a single ancestor, the yeast WGD likely involved mating between two different ancestral species followed by a doubling of the genome to restore fertility.

The ancient whole-genome duplication in the baker's yeast lineage was described almost 20 years ago; in this Primer, one of its discoverers assesses the implications of evidence that it was preceded by hybridization between two species.

Saccharomyces cerevisiae Ski7 Is a GTP-Binding Protein Adopting the Characteristic Conformation of Active Translational GTPases

Ski7 is a cofactor of the cytoplasmic exosome in budding yeast, functioning in both mRNA turnover and non-stop decay (NSD), a surveillance pathway that degrades faulty mRNAs lacking a stop codon. The C-terminal region of Ski7 (Ski7C) shares overall sequence similarity with the translational GTPase (trGTPase) Hbs1, but whether Ski7 has retained the properties of a trGTPase is unclear. Here, we report the high-resolution structures of Ski7C bound to either intact guanosine triphosphate (GTP) or guanosine diphosphate-Pi. The individual domains of Ski7C adopt the conformation characteristic of active trGTPases. Furthermore, the nucleotide-binding site of Ski7C shares similar features compared with active trGTPases, notably the presence of a characteristic monovalent cation. However, a suboptimal polar residue at the putative catalytic site and an unusual polar residue that interacts with the γ-phosphate of GTP distinguish Ski7 from other trGTPases, suggesting it might function rather as a GTP-binding protein than as a GTP-hydrolyzing enzyme.

The evolutionary fate of alternatively spliced homologous exons after gene duplication

Alternative splicing and gene duplication are the two main processes responsible for expanding protein functional diversity. Although gene duplication can generate new genes and alternative splicing can introduce variation through alternative gene products, the interplay between the two processes is complex and poorly understood. Here, we have carried out a study of the evolution of alternatively spliced exons after gene duplication to better understand the interaction between the two processes. We created a manually curated set of 97 human genes with mutually exclusively spliced homologous exons and analyzed the evolution of these exons across five distantly related vertebrates (lamprey, spotted gar, zebrafish, fugu, and coelacanth). Most of these exons had an ancient origin (more than 400 Ma). We found examples supporting two extreme evolutionary models for the behaviour of homologous axons after gene duplication. We observed 11 events in which gene duplication was accompanied by splice isoform separation, that is, each paralog specifically conserved just one distinct ancestral homologous exon. At other extreme, we identified genes in which the homologous exons were always conserved within paralogs, suggesting that the alternative splicing event cannot easily be separated from the function in these genes. That many homologous exons fall in between these two extremes highlights the diversity of biological systems and suggests that the subtle balance between alternative splicing and gene duplication is adjusted to the specific cellular context of each gene.

Evidence for widespread subfunctionalization of splice forms in vertebrate genomes

Gene duplication and alternative splicing are important sources of proteomic diversity. Despite research indicating that gene duplication and alternative splicing are negatively correlated, the evolutionary relationship between the two remains unclear. One manner in which alternative splicing and gene duplication may be related is through the process of subfunctionalization, in which an alternatively spliced gene upon duplication divides distinct splice isoforms among the newly generated daughter genes, in this way reducing the number of alternatively spliced transcripts duplicate genes produce. Previously, it has been shown that splice form subfunctionalization will result in duplicate pairs with divergent exon structure when distinct isoforms become fixed in each paralog. However, the effects of exon structure divergence between paralogs have never before been studied on a genome-wide scale. Here, using genomic data from human, mouse, and zebrafish, we demonstrate that gene duplication followed by exon structure divergence between paralogs results in a significant reduction in levels of alternative splicing. In addition, by comparing the exon structure of zebrafish duplicates to the co-orthologous human gene, we have demonstrated that a considerable fraction of exon divergent duplicates maintain the structural signature of splice form subfunctionalization. Furthermore, we find that paralogs with divergent exon structure demonstrate reduced breadth of expression in a variety of tissues when compared to paralogs with identical exon structures and singletons. Taken together, our results are consistent with subfunctionalization partitioning alternatively spliced isoforms among duplicate genes and as such highlight the relationship between gene duplication and alternative splicing.