Molecular evolution of eukaryotic genomes: hemiascomycetous yeast spliceosomal introns

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.

Missplicing suppressor alleles of Arabidopsis PRE-MRNA PROCESSING FACTOR 8 increase splicing fidelity by reducing the use of novel splice sites

Efficient splicing requires a balance between high-fidelity splice-site (SS) selection and speed. In Saccharomyces cerevisiae, Pre-mRNA processing factor 8 (Prp8) helps to balance precise SS selection and rapid, efficient intron excision and exon joining. argonaute1-52 (ago1-52) and incurvata13 (icu13) are hypomorphic alleles of the Arabidopsis thaliana genes ARGONAUTE1 (AGO1) and AUXIN RESISTANT6 (AXR6) that harbor point mutations creating a novel 3'SS and 5'SS, respectively. The spliceosome recognizes these novel SSs, as well as the intact genuine SSs, producing a mixture of wild-type and aberrant mature mRNAs. Here, we characterized five novel mutant alleles of PRP8 (one of the two Arabidopsis co-orthologs of yeast Prp8), naming these alleles morphology of ago1-52 suppressed5 (mas5). In the mas5-1 background, the spliceosome preferentially recognizes the intact genuine 3'SS of ago1-52 and 5'SS of icu13. Since point mutations that damage genuine SSs make the spliceosome prone to recognizing cryptic SSs, we also tested alleles of four genes carrying damaged genuine SSs, finding that mas5-1 did not suppress their missplicing. The mas5-1 and mas5-3 mutations represent a novel class of missplicing suppressors that increase splicing fidelity by hampering the use of novel SSs, but do not alter general pre-mRNA splicing.

The Arabidopsis cyclophilin CYP18-1 facilitates PRP18 dephosphorylation and the splicing of introns retained under heat stress

In plants, heat stress induces changes in alternative splicing, including intron retention; these events can rapidly alter proteins or downregulate protein activity, producing nonfunctional isoforms or inducing nonsense-mediated decay of messenger RNA (mRNA). Nuclear cyclophilins (CYPs) are accessory proteins in the spliceosome complexes of multicellular eukaryotes. However, whether plant CYPs are involved in pre-mRNA splicing remain unknown. Here, we found that Arabidopsis thaliana CYP18-1 is necessary for the efficient removal of introns that are retained in response to heat stress during germination. CYP18-1 interacts with Step II splicing factors (PRP18a, PRP22, and SWELLMAP1) and associates with the U2 and U5 small nuclear RNAs in response to heat stress. CYP18-1 binds to phospho-PRP18a, and increasing concentrations of CYP18-1 are associated with increasing dephosphorylation of PRP18a. Furthermore, interaction and protoplast transfection assays revealed that CYP18-1 and the PP2A-type phosphatase PP2A B'η co-regulate PRP18a dephosphorylation. RNA-seq and RT-qPCR analysis confirmed that CYP18-1 is essential for splicing introns that are retained under heat stress. Overall, we reveal the mechanism of action by which CYP18-1 activates the dephosphorylation of PRP18 and show that CYP18-1 is crucial for the efficient splicing of retained introns and rapid responses to heat stress in plants.

Analysis of fungal genomes reveals commonalities of intron gain or loss and functions in intron-poor species

Previous evolutionary reconstructions have concluded that early eukaryotic ancestors including both the last common ancestor of eukaryotes and of all fungi had intron-rich genomes. By contrast, some extant eukaryotes have few introns, underscoring the complex histories of intron–exon structures, and raising the question as to why these few introns are retained. Here, we have used recently available fungal genomes to address a variety of questions related to intron evolution. Evolutionary reconstruction of intron presence and absence using 263 diverse fungal species supports the idea that massive intron reduction through intron loss has occurred in multiple clades. The intron densities estimated in various fungal ancestors differ from zero to 7.6 introns per 1 kb of protein-coding sequence. Massive intron loss has occurred not only in microsporidian parasites and saccharomycetous yeasts, but also in diverse smuts and allies. To investigate the roles of the remaining introns in highly-reduced species, we have searched for their special characteristics in eight intron-poor fungi. Notably, the introns of ribosome-associated genes RPL7 and NOG2 have conserved positions; both intron-containing genes encoding snoRNAs. Furthermore, both the proteins and snoRNAs are involved in ribosome biogenesis, suggesting that the expression of the protein-coding genes and noncoding snoRNAs may be functionally coordinated. Indeed, these introns are also conserved in three-quarters of fungi species. Our study shows that fungal introns have a complex evolutionary history and underappreciated roles in gene expression.

Modular engineering of Shiraia bambusicola for hypocrellin production through an efficient CRISPR system

Shiraia bambusicola exhibits an excellent capability to produce high-value pharmacological drugs, such as hypocrellin. However, less effective molecular tools hamper the processes to discover or exploit these metabolites. To address this issue, the more effective CRISPR/Cas9 system was constructed by optimizing the sgRNA transcription elements and disrupting the endogenous non-homologous end-joining pathway. These tactics prompted the gene-targeting frequency of 100% and simultaneously multiplex genome editing in S. bambusicola. This optimal CRISPR system encouraged us to rewire the entire hypocrellin flux and improve the yield by orchestrating the substrate pool supply, the central hypocrellin pathway, and the antioxidant system. Thus, 8632 mg/L hypocrellin was obtained, resulting in a 12-fold increase than that of the wild-type strain. This engineered S. bambusicola can still endure oxidative stresses from higher target metabolites and sustain an excellent biological activity. This study provides a whole conception to establish the more efficient genome-editing system. Higher conserved transcription elements for sgRNA expressions inspire us to adopt this system for gene modifications of other filamentous fungi. The rational and global biosystems outline will offer guidance to modulate metabolite productivity in other filamentous fungi.

Animal, Fungi, and Plant Genome Sequences Harbor Different Non-Canonical Splice Sites

Most protein-encoding genes in eukaryotes contain introns, which are interwoven with exons. Introns need to be removed from initial transcripts in order to generate the final messenger RNA (mRNA), which can be translated into an amino acid sequence. Precise excision of introns by the spliceosome requires conserved dinucleotides, which mark the splice sites. However, there are variations of the highly conserved combination of GT at the 5' end and AG at the 3' end of an intron in the genome. GC-AG and AT-AC are two major non-canonical splice site combinations, which have been known for years. Recently, various minor non-canonical splice site combinations were detected with numerous dinucleotide permutations. Here, we expand systematic investigations of non-canonical splice site combinations in plants across eukaryotes by analyzing fungal and animal genome sequences. Comparisons of splice site combinations between these three kingdoms revealed several differences, such as an apparently increased CT-AC frequency in fungal genome sequences. Canonical GT-AG splice site combinations in antisense transcripts are a likely explanation for this observation, thus indicating annotation errors. In addition, high numbers of GA-AG splice site combinations were observed in Eurytemoraaffinis and Oikopleuradioica. A variant in one U1 small nuclear RNA (snRNA) isoform might allow the recognition of GA as a 5' splice site. In depth investigation of splice site usage based on RNA-Seq read mappings indicates a generally higher flexibility of the 3' splice site compared to the 5' splice site across animals, fungi, and plants.

Patterns of conservation of spliceosomal intron structures and spliceosome divergence in representatives of the diplomonad and parabasalid lineages

Background

Two spliceosomal intron types co-exist in eukaryotic precursor mRNAs and are excised by distinct U2-dependent and U12-dependent spliceosomes. In the diplomonad Giardia lamblia, small nuclear (sn) RNAs show hybrid characteristics of U2- and U12-dependent spliceosomal snRNAs and 5 of 11 identified remaining spliceosomal introns are trans-spliced. It is unknown whether unusual intron and spliceosome features are conserved in other diplomonads.

Results

We have identified spliceosomal introns, snRNAs and proteins from two additional diplomonads for which genome information is currently available, Spironucleus vortens and Spironucleus salmonicida, as well as relatives, including 6 verified cis-spliceosomal introns in S. vortens. Intron splicing signals are mostly conserved between the Spironucleus species and G. lamblia. Similar to ‘long’ G. lamblia introns, RNA secondary structural potential is evident for ‘long’ (> 50 nt) Spironucleus introns as well as introns identified in the parabasalid Trichomonas vaginalis. Base pairing within these introns is predicted to constrain spatial distances between splice junctions to similar distances seen in the shorter and uniformly-sized introns in these organisms. We find that several remaining Spironucleus spliceosomal introns are ancient. We identified a candidate U2 snRNA from S. vortens, and U2 and U5 snRNAs in S. salmonicida; cumulatively, illustrating significant snRNA differences within some diplomonads. Finally, we studied spliceosomal protein complements and find protein sets in Giardia, Spironucleus and Trepomonas sp. PC1 highly- reduced but well conserved across the clade, with between 44 and 62 out of 174 studied spliceosomal proteins detectable. Comparison with more distant relatives revealed a highly nested pattern, with the more intron-rich fornicate Kipferlia bialata retaining 87 total proteins including nearly all those observed in the diplomonad representatives, and the oxymonad Monocercomonoides retaining 115 total proteins including nearly all those observed in K. bialata.

Conclusions

Comparisons in diplomonad representatives and species of other closely-related metamonad groups indicates similar patterns of intron structural conservation and spliceosomal protein composition but significant divergence of snRNA structure in genomically-reduced species. Relative to other eukaryotes, loss of evolutionarily-conserved snRNA domains and common sets of spliceosomal proteins point to a more streamlined splicing mechanism, where intron sequences and structures may be functionally compensating for the minimalization of spliceosome components.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1488-y) contains supplementary material, which is available to authorized users.

Evolution of Alternative Splicing in Eudicots

Alternative pre-mRNA splicing (AS) is prevalent in plants and is involved in many interactions between plants and environmental stresses. However, the patterns and underlying mechanisms of AS evolution in plants remain unclear. By analyzing the transcriptomes of four eudicot species, we revealed that the divergence of AS is largely due to the gains and losses of AS events among orthologous genes. Furthermore, based on a subset of AS, in which AS can be directly associated with specific transcripts, we found that AS that generates transcripts containing premature termination codons (PTC), are likely more conserved than those that generate non-PTC containing transcripts. This suggests that AS coupled with nonsense-mediated decay (NMD) might play an important role in affecting mRNA levels post-transcriptionally. To understand the mechanisms underlying the divergence of AS, we analyzed the key determinants of AS using a machine learning approach. We found that the presence/absence of alternative splice site (SS) within the junction, the distance between the authentic SS and the nearest alternative SS, the size of exon-exon junctions were the major determinants for both alternative 5' donor site and 3' acceptor site among the studied species, suggesting a relatively conserved AS mechanism. The comparative analysis further demonstrated that variations of the identified AS determinants significantly contributed to the AS divergence among closely related species in both Solanaceae and Brassicaceae taxa. Together, these results provide detailed insights into the evolution of AS in plants.

Molecular Mechanisms of pre-mRNA Splicing through Structural Biology of the Spliceosome

Precursor messenger RNA (pre-mRNA) splicing is executed by the spliceosome. In the past 3 years, cryoelectron microscopy (cryo-EM) structures have been elucidated for a majority of the yeast spliceosomal complexes and for a few human spliceosomes. During the splicing reaction, the dynamic spliceosome has an immobile core of about 20 protein and RNA components, which are organized around a conserved splicing active site. The divalent metal ions, coordinated by U6 small nuclear RNA (snRNA), catalyze the branching reaction and exon ligation. The spliceosome also contains a mobile but compositionally stable group of about 13 proteins and a portion of U2 snRNA, which facilitate substrate delivery into the splicing active site. The spliceosomal transitions are driven by the RNA-dependent ATPase/helicases, resulting in the recruitment and dissociation of specific splicing factors that enable the reaction. In summary, the spliceosome is a protein-directed metalloribozyme.

Genome-wide identification and expression profile analysis of the HOG gene family in Aspergillus oryzae

The High osmolarity glycerol (HOG) gene family plays crucial roles in various developmental and physiological processes in fungi, such as the permeability of cell membrane, chlamydospore formation and stress signaling. Although the function of HOG genes has been investigated in Saccharomyces cerevisiae and some filamentous fungi, a comprehensive analysis of HOG gene family has not been performed in Aspergillus oryzae, a fungi mainly used for the production of soy sauce. In this study, we identified and corrected a total of 90 HOG genes from the A. oryzae genome. According to the phylogenetic relationship, they were divided into four discrete groups (Group A-D) comprising of 16, 24, 30 and 20 proteins, respectively. Six conserved motifs and exon-intron structures were examined among all HOG proteins to reveal the diversity of AoHOG genes. Based on transcriptome technology, the expression patterns of AoHOG genes across all developmental stages was identified, suggesting that the AoHOG gene family mainly functions in the logarithmic phase of development. The expression profiles of AoHOG genes under different concentrations of salt stress indicated that AoHOG genes are extensively involved in salt stress response, with possibly different mechanisms. The genome-wide identification, evolutionary, gene structures and expression analyses of AoHOG genes provide a comprehensive overview of this gene family as well as their potential involvements in development and stress responses. Our results will facilitate further research on HOG gene family regarding their physiological and biochemical functions.

Structure of the human activated spliceosome in three conformational states

During each cycle of pre-mRNA splicing, the pre-catalytic spliceosome (B complex) is converted into the activated spliceosome (B^act complex), which has a well-formed active site but cannot proceed to the branching reaction. Here, we present the cryo-EM structure of the human B^act complex in three distinct conformational states. The EM map allows atomic modeling of nearly all protein components of the U2 small nuclear ribonucleoprotein (snRNP), including three of the SF3a complex and seven of the SF3b complex. The structure of the human B^act complex contains 52 proteins, U2, U5, and U6 small nuclear RNA (snRNA), and a pre-mRNA. Three distinct conformations have been captured, representing the early, mature, and late states of the human B^act complex. These complexes differ in the orientation of the Switch loop of Prp8, the splicing factors RNF113A and NY-CO-10, and most components of the NineTeen complex (NTC) and the NTC-related complex. Analysis of these three complexes and comparison with the B and C complexes reveal an ordered flux of components in the B-to-B^act and the B^act-to-B^* transitions, which ultimately prime the active site for the branching reaction.

Structure of an Intron Lariat Spliceosome from Saccharomyces cerevisiae

The disassembly of the intron lariat spliceosome (ILS) marks the end of a splicing cycle. Here we report a cryoelectron microscopy structure of the ILS complex from Saccharomyces cerevisiae at an average resolution of 3.5 Å. The intron lariat remains bound in the spliceosome whereas the ligated exon is already dissociated. The step II splicing factors Prp17 and Prp18, along with Cwc21 and Cwc22 that stabilize the 5' exon binding to loop I of U5 small nuclear RNA (snRNA), have been released from the active site assembly. The DEAH family ATPase/helicase Prp43 binds Syf1 at the periphery of the spliceosome, with its RNA-binding site close to the 3' end of U6 snRNA. The C-terminal domain of Ntr1/Spp382 associates with the GTPase Snu114, and Ntr2 is anchored to Prp8 while interacting with the superhelical domain of Ntr1. These structural features suggest a plausible mechanism for the disassembly of the ILS complex.

Posttranscriptional Regulation of Gcr1 Expression and Activity Is Crucial for Metabolic Adjustment in Response to Glucose Availability

The transcription factor Gcr1 controls expression of over 75% of the genes in actively growing yeast. Yet despite its widespread effects, regulation of Gcr1 itself remains poorly understood. Here, we show that posttranscriptional Gcr1 regulation is nutrient dependent. Moreover, GCR1 RNA contains a long, highly conserved intron, which allows the cell to generate multiple RNA and protein isoforms whose levels change upon glucose depletion. Intriguingly, an isoform generated by intron retention is exported from the nucleus, and its translation is initiated from a conserved, intronic translation start site. Expression of gene products from both the spliced and unspliced RNAs is essential, as cells expressing only one isoform cannot adjust their metabolic program in response to environmental changes. Finally, we show that the Gcr1 proteins form dimers, providing an elegant mechanism by which this one gene, through its regulation, can perform the repertoire of transcriptional activities necessary for fine-tuned environmental response.

Genome Diversity and Evolution in the Budding Yeasts (Saccharomycotina)

Considerable progress in our understanding of yeast genomes and their evolution has been made over the last decade with the sequencing, analysis, and comparisons of numerous species, strains, or isolates of diverse origins. The role played by yeasts in natural environments as well as in artificial manufactures, combined with the importance of some species as model experimental systems sustained this effort. At the same time, their enormous evolutionary diversity (there are yeast species in every subphylum of Dikarya) sparked curiosity but necessitated further efforts to obtain appropriate reference genomes. Today, yeast genomes have been very informative about basic mechanisms of evolution, speciation, hybridization, domestication, as well as about the molecular machineries underlying them. They are also irreplaceable to investigate in detail the complex relationship between genotypes and phenotypes with both theoretical and practical implications. This review examines these questions at two distinct levels offered by the broad evolutionary range of yeasts: inside the best-studied Saccharomyces species complex, and across the entire and diversified subphylum of Saccharomycotina. While obviously revealing evolutionary histories at different scales, data converge to a remarkably coherent picture in which one can estimate the relative importance of intrinsic genome dynamics, including gene birth and loss, vs. horizontal genetic accidents in the making of populations. The facility with which novel yeast genomes can now be studied, combined with the already numerous available reference genomes, offer privileged perspectives to further examine these fundamental biological questions using yeasts both as eukaryotic models and as fungi of practical importance.

New kids on the block: emerging oleaginous yeast of biotechnological importance

There is growing interest in using oleaginous yeast for the production of a variety of fatty acids and fatty acid-derived oleochemicals. This is motivated by natural propensity for high flux through lipid biosynthesis that has naturally evolved, making them a logical starting point for additional genetic engineering to improve titers and productivities. Much of the academic and industrial focus has centered on yeast that have significant genetic engineering tool capabilities, such as Yarrowia lipolytica, and those that have naturally high lipid accumulation, such as Rhodosporidium toruloides and Lipomyces starkeyi; however, there are oleaginous yeast with phenotypes better aligned with typically inhibitory process conditions, such as high salt concentrations and lignocellulosic derived inhibitors. This review addresses the foundational work in characterizing two emerging oleaginous yeast of interest: Debaryomyces hansenii and Trichosporon oleaginosus. We focus on the physiological and metabolic properties of these yeast that make each attractive for bioprocessing of lignocellulose to fuels and chemicals, discuss their respective genetic engineering tools and highlight the critical barriers facing the broader implementation of these oleaginous yeast.

Identification of a second PAD1 in Brettanomyces bruxellensis LAMAP2480

Volatile phenols are aromatic compounds produced by some yeasts of the genus Brettanomyces as defense against the toxicity of hydroxycinnamic acids (p-coumaric acid, ferulic acid and caffeic acid). The origin of these compounds in winemaking involves the sequential action of two enzymes: coumarate decarboxylase and vinylphenol reductase. The first one converts hydroxycinnamic acids into hydroxystyrenes, which are then reduced to ethyl derivatives by vinylphenol reductase. Volatile phenols derived from p-coumaric acid (4-vinylphenol and 4-ethylphenol) have been described as the major contributors to self-defeating aromas associated with stable, gouache, wet mouse, etc., which generates large economic losses in the wine industry. The gene responsible for the production of 4-vinylphenol from p-coumaric acid has been identified as PAD1, which encodes a phenylacrylic acid decarboxylase. PAD1 has been described for many species, among them Candida albicans, Candida dubliniensis, Debaryomyces hansenii and Pichia anomala. In Brettanomyces bruxellensis LAMAP2480, a 666 bp reading frame (DbPAD) encodes a coumarate decarboxylase. Recent studies have reported the existence of a new reading frame belonging to DbPAD called DbPAD2 of 531 bp, which could encode a protein with similar enzymatic activity to PAD1. The present study confirmed that the transformation of Saccharomyces cerevisiae strain BY4722 with reading frame DbPAD2 under the control of the B. bruxellensis ACT1 promoter, encodes an enzyme with coumarate decarboxylase activity. This work has provided deeper insight into the origin of aroma defects in wine due to contamination by Brettanomyces spp.

A native promoter and inclusion of an intron is necessary for efficient expression of GFP or mRFP in Armillaria mellea

Armillaria mellea is a significant pathogen that causes Armillaria root disease on numerous hosts in forests, gardens and agricultural environments worldwide. Using a yeast-adapted pCAMBIA0380 Agrobacterium vector, we have constructed a series of vectors for transformation of A. mellea, assembled using yeast-based recombination methods. These have been designed to allow easy exchange of promoters and inclusion of introns. The vectors were first tested by transformation into basidiomycete Clitopilus passeckerianus to ascertain vector functionality then used to transform A. mellea. We show that heterologous promoters from the basidiomycetes Agaricus bisporus and Phanerochaete chrysosporium that were used successfully to control the hygromycin resistance cassette were not able to support expression of mRFP or GFP in A. mellea. The endogenous A. mellea gpd promoter delivered efficient expression, and we show that inclusion of an intron was also required for transgene expression. GFP and mRFP expression was stable in mycelia and fluorescence was visible in transgenic fruiting bodies and GFP was detectable in planta. Use of these vectors has been successful in giving expression of the fluorescent proteins GFP and mRFP in A. mellea, providing an additional molecular tool for this pathogen.

Tyrosine Recombinase Retrotransposons and Transposons

Retrotransposons carrying tyrosine recombinases (YR) are widespread in eukaryotes. The first described tyrosine recombinase mobile element, DIRS1, is a retroelement from the slime mold Dictyostelium discoideum. The YR elements are bordered by terminal repeats related to their replication via free circular dsDNA intermediates. Site-specific recombination is believed to integrate the circle without creating duplications of the target sites. Recently a large number of YR retrotransposons have been described, including elements from fungi (mucorales and basidiomycetes), plants (green algae) and a wide range of animals including nematodes, insects, sea urchins, fish, amphibia and reptiles. YR retrotransposons can be divided into three major groups: the DIRS elements, PAT-like and the Ngaro elements. The three groups form distinct clades on phylogenetic trees based on alignments of reverse transcriptase/ribonuclease H (RT/RH) and YR sequences, and also having some structural distinctions. A group of eukaryote DNA transposons, cryptons, also carry tyrosine recombinases. These DNA transposons do not encode a reverse transcriptase. They have been detected in several pathogenic fungi and oomycetes. Sequence comparisons suggest that the crypton YRs are related to those of the YR retrotransposons. We suggest that the YR retrotransposons arose from the combination of a crypton-like YR DNA transposon and the RT/RH encoding sequence of a retrotransposon. This acquisition must have occurred at a very early point in the evolution of eukaryotes.

Evolutionarily conserved autoregulation of alternative pre-mRNA splicing by ribosomal protein L10a

Alternative splicing of pre-mRNAs can regulate expression of protein-coding genes by generating unproductive mRNAs rapidly degraded by nonsense-mediated mRNA decay (NMD). Many of the genes directly regulated by alternative splicing coupled with NMD (AS-NMD) are related to RNA metabolism, but the repertoire of genes regulated by AS-NMD in vivo is to be determined. Here, we analyzed transcriptome data of wild-type and NMD-defective mutant strains of the nematode worm Caenorhabditis elegans and demonstrate that eight of the 82 cytoplasmic ribosomal protein (rp) genes generate unproductively spliced mRNAs. Knockdown of any of the eight rp genes exerted a dynamic and compensatory effect on alternative splicing of its own transcript and inverse effects on that of the other rp genes. A large subunit protein L10a, termed RPL-1 in nematodes, directly and specifically binds to an evolutionarily conserved 39-nt stretch termed L10ARE between the two alternative 5′ splice sites in its own pre-mRNA to switch the splice site choice. Furthermore, L10ARE-mediated splicing autoregulation of the L10a-coding gene is conserved in vertebrates. These results indicate that L10a is an evolutionarily conserved splicing regulator and that homeostasis of a subset of the rp genes are regulated at the level of pre-mRNA splicing in vivo.

The 3.8 Å structure of the U4/U6.U5 tri-snRNP: Insights into spliceosome assembly and catalysis

Splicing of precursor messenger RNA is accomplished by a dynamic megacomplex known as the spliceosome. Assembly of a functional spliceosome requires a preassembled U4/U6.U5 tri-snRNP complex, which comprises the U5 small nuclear ribonucleoprotein (snRNP), the U4 and U6 small nuclear RNA (snRNA) duplex, and a number of protein factors. Here we report the three-dimensional structure of a Saccharomyces cerevisiae U4/U6.U5 tri-snRNP at an overall resolution of 3.8 angstroms by single-particle electron cryomicroscopy. The local resolution for the core regions of the tri-snRNP reaches 3.0 to 3.5 angstroms, allowing construction of a refined atomic model. Our structure contains U5 snRNA, the extensively base-paired U4/U6 snRNA, and 30 proteins including Prp8 and Snu114, which amount to 8495 amino acids and 263 nucleotides with a combined molecular mass of ~1 megadalton. The catalytic nucleotide U80 from U6 snRNA exists in an inactive conformation, stabilized by its base-pairing interactions with U4 snRNA and protected by Prp3. Pre-messenger RNA is bound in the tri-snRNP through base-pairing interactions with U6 snRNA and loop I of U5 snRNA. This structure, together with that of the spliceosome, reveals the molecular choreography of the snRNAs in the activation process of the spliceosomal ribozyme.

Molecular characterization and heterologous expression of a Xanthophyllomyces dendrorhous α-glucosidase with potential for prebiotics production

Basidiomycetous yeast Xanthophyllomyces dendrorhous expresses an α-glucosidase with strong transglycosylation activity producing prebiotic sugars such as panose and an unusual tetrasaccharides mixture including α-(1-6) bonds as major products, which makes it of biotechnological interest. Initial analysis pointed to a homodimeric protein of 60 kDa subunit as responsible for this activity. In this study, the gene Xd-AlphaGlu was characterized. The 4131-bp-long gene is interrupted by 13 short introns and encodes a protein of 990 amino acids (Xd-AlphaGlu). The N-terminal sequence of the previously detected 60 kDa protein resides in this larger protein at residues 583-602. Functionality of the gene was proved in Saccharomyces cerevisiae, which produced a protein of about 130 kDa containing Xd-AlphaGlu sequences. All properties of the heterologously expressed protein, including thermal and pH profiles, activity on different substrates, and ability to produce prebiotic sugars were similar to that of the α-glucosidase produced in X. dendrorhous. No activity was detected in S. cerevisiae containing exclusively the 1256-bp from gene Xd-AlphaGlu that would encode synthesis of the 60 kDa protein previously detected. Data were compatible with an active monomeric α-glucosidase of 990 amino acids and an inactive hydrolysis product of 60 kDa. Protein Xd-AlphaGlu contained most of the elements characteristic of α-glucosidases included in the glycoside hydrolases family GH31 and its structural model based on the homologous human maltase-glucoamylase was obtained. Remarkably, the Xd-AlphaGlu C-terminal domain presents an unusually long 115-residue insertion that could be involved in this enzyme's activity against long-size substrates such as maltoheptaose and soluble starch.

Genome editing in Ustilago maydis using the CRISPR-Cas system

This communication describes the establishment of the type II bacterial CRISPR-Cas9 system to efficiently disrupt target genes in the fungal maize pathogen Ustilago maydis. A single step transformation of a self-replicating plasmid constitutively expressing the U. maydis codon-optimized cas9 gene and a suitable sgRNA under control of the U. maydis U6 snRNA promoter was sufficient to induce genome editing. On average 70% of the progeny of a single transformant were disrupted within the respective b gene. Without selection the self-replicating plasmid was lost rapidly allowing transient expression of the CRISPR-Cas9 system to minimize potential long-term negative effects of Cas9. This technology will be an important advance for the simultaneous disruption of functionally redundant genes and gene families to investigate their contribution to virulence of U. maydis.

Intron evolution in Saccharomycetaceae

Introns in protein-coding genes are very rare in hemiascomycetous yeast genomes. It has been suggested that these species have experienced extensive intron loss during their evolution from the postulated intron-rich fungal ancestor. However, no intron-devoid yeast species have been identified and some of the introns remaining within the genomes of intron-poor species, such as Saccharomyces cerevisiae, appear to be beneficial during growth under stress conditions. In order to reveal the pattern of intron retention within intron-poor yeast species and better understand the mechanisms of intron evolution, we generated a comprehensive set of 250 orthologous introns in the 20 species that comprise the Saccharomycetaceae, by analyzing RNA deep-sequencing data and alignments of intron-containing genes. Analysis of these intron sets shows that intron loss is at least two orders of magnitude more frequent than intron gain. Fine mapping of intron positions shows that intron sliding is rare, and that introns are almost always removed without changing the primary sequence of the encoded protein. The latter finding is consistent with the prevailing view that homologous recombination between reverse-transcribed mature mRNAs and the corresponding genomic locus is the primary mechanism of intron loss. However, we also find evidence that loss of a small number of introns is mediated by micro-homology, and that the number of intron losses is diminished in yeast species that have lost the microhomology end joining and nonhomologous end joining machinery.

Structure and Organization of the Engraulidae Family U2 snRNA: An Evolutionary Model Gene?

The U2 snRNA multigene family has been analyzed in four species of the Engraulidae family--Engraulis encrasicolus, Engraulis mordax, Engraulis ringens, and Engraulis japonicas--with the object of understanding more about the structure of this multigene family in these pelagic species and studying their phylogenetic relationships. The results showed that the cluster of this gene family in the Engraulis genus is formed by the U2-U5 snRNA with highly conserved sequences of mini- and micro-satellites, such as (CTGT)n, embedded downstream of the transcription unit; findings indicate that this gene family evolved following the concerted model. The phylogenetic analysis of the non-transcribed spacer of cluster U2-U5 snDNA in the 4 species showed that the sequences of the species E. encrasicolus and E. japonicus are closely related; these two are genetically close to E. mordax and slightly more distant from E. ringens. The data obtained by molecular analysis of U2-U5 snDNA and their secondary structure, with the presence of the micro-satellite (CTGT)n and mini-satellites, show clearly that the species E. encrasicolus and E. japonicus are closely related and would be older than E. mordax and E. ringens.

The pre-mRNA splicing reaction

In eukaryotic organisms, nascent transcripts of protein-coding genes contain intronic sequences that are not present in mature mRNAs. Pre-mRNA splicing removes introns and joins exons to form mature mRNAs. It is catalyzed by a large RNP complex called the spliceosome. Sequences within the pre-mRNA determine intron recognition and excision. This process occurs with a high degree of accuracy to generate the functional transcriptome of a cell.

Phylogenetics and Gene Structure Dynamics of Polygalacturonase Genes in Aspergillus and Neurospora crassa

Polygalacturonase (PG) gene is a typical gene family present in eukaryotes. Forty-nine PGs were mined from the genomes of Neurospora crassa and five Aspergillus species. The PGs were classified into 3 clades such as clade 1 for rhamno-PGs, clade 2 for exo-PGs and clade 3 for exo- and endo-PGs, which were further grouped into 13 sub-clades based on the polypeptide sequence similarity. In gene structure analysis, a total of 124 introns were present in 44 genes and five genes lacked introns to give an average of 2.5 introns per gene. Intron phase distribution was 64.5% for phase 0, 21.8% for phase 1, and 13.7% for phase 2, respectively. The introns varied in their sequences and their lengths ranged from 20 bp to 424 bp with an average of 65.9 bp, which is approximately half the size of introns in other fungal genes. There were 29 homologous intron blocks and 26 of those were sub-clade specific. Intron losses were counted in 18 introns in which no obvious phase preference for intron loss was observed. Eighteen introns were placed at novel positions, which is considerably higher than those of plant PGs. In an evolutionary sense both intron loss and gain must have taken place for shaping the current PGs in these fungi. Together with the small intron size, low conservation of homologous intron blocks and higher number of novel introns, PGs of fungal species seem to have recently undergone highly dynamic evolution.

Whence genes in pieces: reconstruction of the exon-intron gene structures of the last eukaryotic common ancestor and other ancestral eukaryotes

In eukaryotes, protein-coding sequences are interrupted by non-coding sequences known as introns. During mRNA maturation, introns are excised by the spliceosome and the coding regions, exons, are spliced to form the mature coding region. The intron densities widely differ between eukaryotic lineages, from 6 to 7 introns per kb of coding sequence in vertebrates, some invertebrates and green plants, to only a few introns across the entire genome in many unicellular eukaryotes. Evolutionary reconstructions using maximum likelihood methods suggest intron-rich ancestors for each major group of eukaryotes. For the last common ancestor of animals, the highest intron density of all extant and extinct eukaryotes was inferred, at 120-130% of the human intron density. Furthermore, an intron density within 53-74% of the human values was inferred for the last eukaryotic common ancestor. Accordingly, evolution of eukaryotic genes in all lines of descent involved primarily intron loss, with substantial gain only at the bases of several branches including plants and animals. These conclusions have substantial biological implications indicating that the common ancestor of all modern eukaryotes was a complex organism with a gene architecture resembling those in multicellular organisms. Alternative splicing most likely initially appeared as an inevitable result of splicing errors and only later was employed to generate structural and functional diversification of proteins.

Origin and evolution of spliceosomal introns

Evolution of exon-intron structure of eukaryotic genes has been a matter of long-standing, intensive debate. The introns-early concept, later rebranded ‘introns first’ held that protein-coding genes were interrupted by numerous introns even at the earliest stages of life's evolution and that introns played a major role in the origin of proteins by facilitating recombination of sequences coding for small protein/peptide modules. The introns-late concept held that introns emerged only in eukaryotes and new introns have been accumulating continuously throughout eukaryotic evolution. Analysis of orthologous genes from completely sequenced eukaryotic genomes revealed numerous shared intron positions in orthologous genes from animals and plants and even between animals, plants and protists, suggesting that many ancestral introns have persisted since the last eukaryotic common ancestor (LECA). Reconstructions of intron gain and loss using the growing collection of genomes of diverse eukaryotes and increasingly advanced probabilistic models convincingly show that the LECA and the ancestors of each eukaryotic supergroup had intron-rich genes, with intron densities comparable to those in the most intron-rich modern genomes such as those of vertebrates. The subsequent evolution in most lineages of eukaryotes involved primarily loss of introns, with only a few episodes of substantial intron gain that might have accompanied major evolutionary innovations such as the origin of metazoa. The original invasion of self-splicing Group II introns, presumably originating from the mitochondrial endosymbiont, into the genome of the emerging eukaryote might have been a key factor of eukaryogenesis that in particular triggered the origin of endomembranes and the nucleus. Conversely, splicing errors gave rise to alternative splicing, a major contribution to the biological complexity of multicellular eukaryotes. There is no indication that any prokaryote has ever possessed a spliceosome or introns in protein-coding genes, other than relatively rare mobile self-splicing introns. Thus, the introns-first scenario is not supported by any evidence but exon-intron structure of protein-coding genes appears to have evolved concomitantly with the eukaryotic cell, and introns were a major factor of evolution throughout the history of eukaryotes. This article was reviewed by I. King Jordan, Manuel Irimia (nominated by Anthony Poole), Tobias Mourier (nominated by Anthony Poole), and Fyodor Kondrashov. For the complete reports, see the Reviewers’ Reports section.

Diverse forms of RPS9 splicing are part of an evolving autoregulatory circuit

Ribosomal proteins are essential to life. While the functions of ribosomal protein-encoding genes (RPGs) are highly conserved, the evolution of their regulatory mechanisms is remarkably dynamic. In Saccharomyces cerevisiae, RPGs are unusual in that they are commonly present as two highly similar gene copies and in that they are over-represented among intron-containing genes. To investigate the role of introns in the regulation of RPG expression, we constructed 16 S. cerevisiae strains with precise deletions of RPG introns. We found that several yeast introns function to repress rather than to increase steady-state mRNA levels. Among these, the RPS9A and RPS9B introns were required for cross-regulation of the two paralogous gene copies, which is consistent with the duplication of an autoregulatory circuit. To test for similar intron function in animals, we performed an experimental test and comparative analyses for autoregulation among distantly related animal RPS9 orthologs. Overexpression of an exogenous RpS9 copy in Drosophila melanogaster S2 cells induced alternative splicing and degradation of the endogenous copy by nonsense-mediated decay (NMD). Also, analysis of expressed sequence tag data from distantly related animals, including Homo sapiens and Ciona intestinalis, revealed diverse alternatively-spliced RPS9 isoforms predicted to elicit NMD. We propose that multiple forms of splicing regulation among RPS9 orthologs from various eukaryotes operate analogously to translational repression of the alpha operon by S4, the distant prokaryotic ortholog. Thus, RPS9 orthologs appear to have independently evolved variations on a fundamental autoregulatory circuit.

Eukaryotic genes are littered with non-coding intervening sequences, or introns, that must be precisely excised from a messenger RNA before it can be properly translated into protein. Despite their ubiquity, the evolution and function of introns remain poorly understood. Consequently, we cannot accurately predict the functions of individual introns in any organism. In this manuscript, we used a combination of comparative genomics and experimental tests to identify functional introns. First, we looked for signatures of selection to identify important introns in the model yeast Saccharomyces cerevisiae, which focused our attention on the introns of ribosomal protein genes. We then genetically deleted these introns to assess their function. Unlike mammalian introns, we found that yeast introns were not required for high levels of gene expression. Instead, particular introns (we focus on those within genes encoding ribosomal protein S9) were required to fine-tune gene expression through autoregulation. Surprisingly, animal orthologs of these genes also use introns to autoregulate through multiple forms of alternative splicing. We speculate that the introns of ribosomal protein genes, in particular, readily evolve means for autoregulation to meet the demanding requirements of ribosomal protein genes to maintain tight control of gene expression.

Intraspecific variations of Dekkera/Brettanomyces bruxellensis genome studied by capillary electrophoresis separation of the intron splice site profiles

In enology, "Brett" character refers to the wine spoilage caused by the yeast Dekkera/Brettanomyces bruxellensis and its production of volatile phenolic off-flavours. However, the spoilage potential of this yeast is strain-dependent. Therefore, a rapid and reliable recognition at the strain level is a key point to avoid serious economic losses. The present work provides an operative tool to assess the genetic intraspecific variation in this species through the use of introns as molecular targets. Firstly, the available partial D./B. bruxellensis genome sequence was investigated in order to build primers annealing to introns 5' splice site sequence (ISS). This analysis allowed the detection of a non-random vocabulary flanking the site and, exploiting this feature, the creation of specific probes for strain discrimination. Secondly, the separation of the intron splice site PCR fragments was obtained throughout the set up of a capillary electrophoresis protocol, giving a 94% repeatability threshold in our experimental conditions. The comparison of results obtained with ISS-PCR/CE versus the ones performed by mtDNA RFLP revealed that the former protocol is more discriminating and allowed a reliable identification at strain level. Actually sixty D./B. bruxellensis isolates were recognised as unique strains, showing a level of similarity below 79% and confirming the high genetic polymorphism existing within the species. Two main clusters were grouped at similarity levels of about 46% and 47%, respectively, showing a poor correlation with the geographic area of isolation. Moreover, from the evolutionary point of view, the proposed technique could determine the frequency of the genome rearrangements that can occur in D./B. bruxellesis populations.

TEMPERATURE-DEPENDENT STRUCTURAL VARIABILITY OF RNAs: SPLICED LEADER RNAs AND THEIR EVOLUTIONARY HISTORY

The structures attained by RNA molecules depend not only on their sequence but also on environmental parameters such as their temperature. So far, this effect has been largely neglected in bioinformatics studies. Here, we show that structural comparisons can be facilitated and more coherent structural models can be obtained when differences in environmental parameters are taken into account. We re-evaluate the secondary structures of the spliced leader (SL) RNAs from the seven eukaryotic phyla in which SL RNA trans-splicing has been described. Adjusting structure prediction to the natural growth temperatures and considering energetically similar secondary structures, we observe striking similarities among Euglenida, Kinetoplastida, Dinophyceae, Cnidaria, Rotifera, Nematoda, Platyhelminthes, and Tunicata that cannot be explained easily by the independent innovation of SL RNAs in each of these phyla. Supplementary Table is available at .

SUS1 introns are required for efficient mRNA nuclear export in yeast

Efficient coupling between mRNA synthesis and export is essential for gene expression. Sus1/ENY2, a component of the SAGA and TREX-2 complexes, is involved in both transcription and mRNA export. While most yeast genes lack introns, we previously reported that yeast SUS1 bears two. Here we show that this feature is evolutionarily conserved and critical for Sus1 function. We determine that while SUS1 splicing is inefficient, it responds to cellular conditions, and intronic mutations either promoting or blocking splicing lead to defects in mRNA export and cell growth. Consistent with this, we find that an intron-less SUS1 only partially rescues sus1Δ phenotypes. Remarkably, splicing of each SUS1 intron is also affected by the presence of the other and by SUS1 exonic sequences. Moreover, by following SUS1 RNA and protein levels we establish that nonsense-mediated decay (NMD) pathway and the splicing factor Mud2 both play a role in SUS1 expression. Our data (and those of the accompanying work by Hossain et al.) provide evidence of the involvement of splicing, translation, and decay in the regulation of early events in mRNP biogenesis; and imply the additional requirement for a balance in splicing isoforms from a single gene.

The intronome of budding yeasts

Whatever their abundance in genomes, spliceosomal introns are the signature of eukaryotic genes. The sequence of Saccharomyces cerevisiae, achieved fifteen years ago, revealed that this yeast has very few introns, but conserved intron boundaries typical for an intron definition mechanism. With the improvement and the development of new sequencing technologies, yeast genomes have been extensively sequenced during the last decade. We took advantage of this plethora of data to compile and assess the intron content of the protein-coding genes of 13 genomes representative of the evolution of hemiascomycetous yeasts. We first observed that intron paucity is a general rule and that the fastest evolving genomes tend to lose their introns more rapidly (e.g. S. cerevisiae versus Yarrowia lipolytica). Noticeable differences were also confirmed for 5' splice sites and branch point sites (BP) as well as for the relative position of the BP. These changes seemed to be correlated with the lineage specific evolution of splicing factors.

Genome-wide association between branch point properties and alternative splicing

The branch point (BP) is one of the three obligatory signals required for pre-mRNA splicing. In mammals, the degeneracy of the motif combined with the lack of a large set of experimentally verified BPs complicates the task of modeling it in silico, and therefore of predicting the location of natural BPs. Consequently, BPs have been disregarded in a considerable fraction of the genome-wide studies on the regulation of splicing in mammals. We present a new computational approach for mammalian BP prediction. Using sequence conservation and positional bias we obtained a set of motifs with good agreement with U2 snRNA binding stability. Using a Support Vector Machine algorithm, we created a model complemented with polypyrimidine tract features, which considerably improves the prediction accuracy over previously published methods. Applying our algorithm to human introns, we show that BP position is highly dependent on the presence of AG dinucleotides in the 3' end of introns, with distance to the 3' splice site and BP strength strongly correlating with alternative splicing. Furthermore, experimental BP mapping for five exons preceded by long AG-dinucleotide exclusion zones revealed that, for a given intron, more than one BP can be chosen throughout the course of splicing. Finally, the comparison between exons of different evolutionary ages and pseudo exons suggests a key role of the BP in the pathway of exon creation in human. Our computational and experimental analyses suggest that BP recognition is more flexible than previously assumed, and it appears highly dependent on the presence of downstream polypyrimidine tracts. The reported association between BP features and the splicing outcome suggests that this, so far disregarded but yet crucial, element buries information that can complement current acceptor site models.

Detection and analysis of alternative splicing in Yarrowia lipolytica reveal structural constraints facilitating nonsense-mediated decay of intron-retaining transcripts

Background

Hemiascomycetous yeasts have intron-poor genomes with very few cases of alternative splicing. Most of the reported examples result from intron retention in Saccharomyces cerevisiae and some have been shown to be functionally significant. Here we used transcriptome-wide approaches to evaluate the mechanisms underlying the generation of alternative transcripts in Yarrowia lipolytica, a yeast highly divergent from S. cerevisiae.

Results

Experimental investigation of Y. lipolytica gene models identified several cases of alternative splicing, mostly generated by intron retention, principally affecting the first intron of the gene. The retention of introns almost invariably creates a premature termination codon, as a direct consequence of the structure of intron boundaries. An analysis of Y. lipolytica introns revealed that introns of multiples of three nucleotides in length, particularly those without stop codons, were underrepresented. In other organisms, premature termination codon-containing transcripts are targeted for degradation by the nonsense-mediated mRNA decay (NMD) machinery. In Y. lipolytica, homologs of S. cerevisiae UPF1 and UPF2 genes were identified, but not UPF3. The inactivation of Y. lipolytica UPF1 and UPF2 resulted in the accumulation of unspliced transcripts of a test set of genes.

Conclusions

Y. lipolytica is the hemiascomycete with the most intron-rich genome sequenced to date, and it has several unusual genes with large introns or alternative transcription start sites, or introns in the 5' UTR. Our results suggest Y. lipolytica intron structure is subject to significant constraints, leading to the under-representation of stop-free introns. Consequently, intron-containing transcripts are degraded by a functional NMD pathway.

Analysis of gene evolution and metabolic pathways using the Candida Gene Order Browser

BACKGROUND

Candida species are the most common cause of opportunistic fungal infection worldwide. Recent sequencing efforts have provided a wealth of Candida genomic data. We have developed the Candida Gene Order Browser (CGOB), an online tool that aids comparative syntenic analyses of Candida species. CGOB incorporates all available Candida clade genome sequences including two Candida albicans isolates (SC5314 and WO-1) and 8 closely related species (Candida dubliniensis, Candida tropicalis, Candida parapsilosis, Lodderomyces elongisporus, Debaryomyces hansenii, Pichia stipitis, Candida guilliermondii and Candida lusitaniae). Saccharomyces cerevisiae is also included as a reference genome.

RESULTS

CGOB assignments of homology were manually curated based on sequence similarity and synteny. In total CGOB includes 65617 genes arranged into 13625 homology columns. We have also generated improved Candida gene sets by merging/removing partial genes in each genome. Interrogation of CGOB revealed that the majority of tandemly duplicated genes are under strong purifying selection in all Candida species. We identified clusters of adjacent genes involved in the same metabolic pathways (such as catabolism of biotin, galactose and N-acetyl glucosamine) and we showed that some clusters are species or lineage-specific. We also identified one example of intron gain in C. albicans.

CONCLUSIONS

Our analysis provides an important resource that is now available for the Candida community. CGOB is available at http://cgob.ucd.ie.

Constrained intron structures in a microsporidian

The 2.9-Mbp genome of the microsporidian Encephalitozoon cuniculi is severely reduced and compacted, possessing only 16 known tiny spliceosomal introns. Based on motif and expression data, intron profiles were constructed to screen the genome. Twenty additional introns were predicted and verified, doubling the previous estimate. We further predict that accurate 3' splice site (3'SS) selection is accomplished via a scanning mechanism with specificity achieved by maintaining a constrained variable length between the branch point motif and 3'SS. Only introns in ribosomal protein genes exhibit positional bias, and we hypothesize that splicing could be regulating expression of these genes. The large set of new introns in non-ribosomal protein genes suggests that current models of intron loss are unlikely sufficient to explain the distribution of introns. Together, these results extend our understanding of the role of intron loss in genome evolution and contribute to a novel model for splice site selection.

The N-terminus of Prp1 (Prp6/U5-102 K) is essential for spliceosome activation in vivo

The spliceosomal protein Prp1 (Prp6/U5-102 K) is necessary for the integrity of pre-catalytic spliceosomal complexes. We have identified a novel regulatory function for Prp1. Expression of mutations in the N-terminus of Prp1 leads to the accumulation of pre-catalytic spliceosomal complexes containing the five snRNAs U1, U2, U5 and U4/U6 and pre-mRNAs. The mutations in the N-terminus, which prevent splicing to occur, include in vitro and in vivo identified phosphorylation sites of Prp4 kinase. These sites are highly conserved in the human ortholog U5-102 K. The results presented here demonstrate that structural integrity of the N-terminus is required to mediate a splicing event, but is not necessary for the assembly of spliceosomes.

Population genomics of intron splicing in 38 Saccharomyces cerevisiae genome sequences

Introns are a ubiquitous feature of eukaryotic genomes, and the dynamics of intron evolution between species has been extensively studied. However, comparatively few analyses have focused on the evolutionary forces shaping patterns of intron variation within species. To better understand the population genetic characteristics of introns, we performed an extensive population genetics analysis on key intron splice sequences obtained from 38 strains of Saccharomyces cerevisiae. As expected, we found that purifying selection is the dominant force governing intron splice sequence evolution in yeast, formally confirming that intron-containing alleles are a mutational liability. In addition, through extensive coalescent simulations, we obtain quantitative estimates of the strength of purifying selection (2N_es ≈ 19) and use diffusion approximations to provide insights into the evolutionary dynamics and sojourn times of newly arising splice sequence mutations in natural yeast populations. In contrast to previous functional studies, evolutionary analyses comparing the prevalence of introns in essential and nonessential genes suggest that introns in nonribosomal protein genes are functionally important and tend to be actively maintained in natural populations of S. cerevisiae. Finally, we demonstrate that heritable variation in splicing efficiency is common in intron-containing genes with splice sequence polymorphisms. More generally, our study highlights the advantages of population genomics analyses for exploring the forces that have generated extant patterns of genome variation and for illuminating basic biological processes.

Molecular cloning and stress-dependent expression of a gene encoding Delta(12)-fatty acid desaturase in the Antarctic microalga Chlorella vulgaris NJ-7

The psychrotrophic Antarctic alga, Chlorella vulgaris NJ-7, grows under an extreme environment of low temperature and high salinity. In an effort to better understand the correlation between fatty acid metabolism and acclimation to Antarctic environment, we analyzed its fatty acid compositions. An extremely high amount of Delta(12) unsaturated fatty acids was identified which prompted us to speculate about the involvement of Delta(12) fatty acid desaturase in the process of acclimation. A full-length cDNA sequence, designated CvFAD2, was isolated from C. vulgaris NJ-7 via reverse transcription polymerase chain reaction (RT-PCR) and RACE methods. Sequence alignment and phylogenetic analysis showed that the gene was homologous to known microsomal Delta(12)-FADs with the conserved histidine motifs. Heterologous expression in yeast was used to confirm the regioselectivity and the function of CvFAD2. Linoleic acid (18:2), normally not present in wild-type yeast cells, was detected in transformants of CvFAD2. The induction of CvFAD2 at an mRNA level under cold stress and high salinity is detected by real-time PCR. The results showed that both temperature and salinity motivated the upregulation of CvFAD2 expression. The accumulation of CvFAD2 increased 2.2-fold at 15 degrees C and 3.9-fold at 4 degrees C compared to the alga at 25 degrees C. Meanwhile a 1.7- and 8.5-fold increase at 3 and 6% NaCl was detected. These data suggest that CvFAD2 is the enzyme responsible for the Delta(12) fatty acids desaturation involved in the adaption to cold and high salinity for Antarctic C. vugaris NJ-7.

Alternative splicing of PTC7 in Saccharomyces cerevisiae determines protein localization

It is well established that higher eukaryotes use alternative splicing to increase proteome complexity. In contrast, Saccharomyces cerevisiae, a single-cell eukaryote, conducts predominantly regulated splicing through retention of nonfunctional introns. In this article we describe our discovery of a functional intron in the PTC7 (YHR076W) gene that can be alternatively spliced to create two mRNAs that code for distinct proteins. These two proteins localize to different cellular compartments and have distinct cellular roles. The protein translated from the spliced mRNA localizes to the mitochondria and its expression is carbon-source dependent. In comparison, the protein translated from the unspliced mRNA contains a transmembrane domain, localizes to the nuclear envelope, and mediates the toxic effects of Latrunculin A exposure. In conclusion, we identified a definitive example of functional alternative splicing in S. cerevisiae that confers a measurable fitness benefit.

Comparative genomics of protoploid Saccharomycetaceae

Our knowledge of yeast genomes remains largely dominated by the extensive studies on Saccharomyces cerevisiae and the consequences of its ancestral duplication, leaving the evolution of the entire class of hemiascomycetes only partly explored. We concentrate here on five species of Saccharomycetaceae, a large subdivision of hemiascomycetes, that we call "protoploid" because they diverged from the S. cerevisiae lineage prior to its genome duplication. We determined the complete genome sequences of three of these species: Kluyveromyces (Lachancea) thermotolerans and Saccharomyces (Lachancea) kluyveri (two members of the newly described Lachancea clade), and Zygosaccharomyces rouxii. We included in our comparisons the previously available sequences of Kluyveromyces lactis and Ashbya (Eremothecium) gossypii. Despite their broad evolutionary range and significant individual variations in each lineage, the five protoploid Saccharomycetaceae share a core repertoire of approximately 3300 protein families and a high degree of conserved synteny. Synteny blocks were used to define gene orthology and to infer ancestors. Far from representing minimal genomes without redundancy, the five protoploid yeasts contain numerous copies of paralogous genes, either dispersed or in tandem arrays, that, altogether, constitute a third of each genome. Ancient, conserved paralogs as well as novel, lineage-specific paralogs were identified.

Evolution of ultrasmall spliceosomal introns in highly reduced nuclear genomes

Intron reduction and loss is a significant component of genome compaction in many eukaryotic lineages, including yeasts, microsporidia, and some nucleomorphs. Nucleomorphs are the extremely reduced relicts of algal endosymbiont nuclei found in two lineages, cryptomonads and chlorarachniophytes. In cryptomonads, introns are rare or even lost altogether. In contrast, the nucleomorph of the chlorarachniophyte Bigelowiella natans contains the smallest nuclear genome known but paradoxically also retained over 800 tiny spliceosomal introns, ranging from 18 to 21 nt in length. Because introns have not been described in any other chlorarachniophyte nucleomorph, we do not know when these introns were reduced or whether they have been lost in other lineages. To gain insight into the evolution of these unique introns, we sequenced more than 150 spliceosomal introns in the nucleomorph of the chlorarachniophyte Gymnochlora stellata and compared size distribution, sequence features, and patterns of gain/loss. To clarify the possible relationship between intron size and splicing efficiency, we also analyzed the outcome of 580 splicing events. Overall, these data indicate that the radical intron size reduction took place in the ancestor of all extant chlorarachniophytes and that although most introns have been retained through this reductive process, intron loss has also occurred. We also show that intron size is not static, and splicing is not determined strictly by size, but that size does play a strong role in splicing efficiency, likely as part of a combination of sequence features and size.