Comparative analysis of information contents relevant to recognition of introns in many species

DRANetSplicer: A Splice Site Prediction Model Based on Deep Residual Attention Networks

The precise identification of splice sites is essential for unraveling the structure and function of genes, constituting a pivotal step in the gene annotation process. In this study, we developed a novel deep learning model, DRANetSplicer, that integrates residual learning and attention mechanisms for enhanced accuracy in capturing the intricate features of splice sites. We constructed multiple datasets using the most recent versions of genomic data from three different organisms, Oryza sativa japonica, Arabidopsis thaliana and Homo sapiens. This approach allows us to train models with a richer set of high-quality data. DRANetSplicer outperformed benchmark methods on donor and acceptor splice site datasets, achieving an average accuracy of (96.57%, 95.82%) across the three organisms. Comparative analyses with benchmark methods, including SpliceFinder, Splice2Deep, Deep Splicer, EnsembleSplice, and DNABERT, revealed DRANetSplicer's superior predictive performance, resulting in at least a (4.2%, 11.6%) relative reduction in average error rate. We utilized the DRANetSplicer model trained on O. sativa japonica data to predict splice sites in A. thaliana, achieving accuracies for donor and acceptor sites of (94.89%, 94.25%). These results indicate that DRANetSplicer possesses excellent cross-organism predictive capabilities, with its performance in cross-organism predictions even surpassing that of benchmark methods in non-cross-organism predictions. Cross-organism validation showcased DRANetSplicer's excellence in predicting splice sites across similar organisms, supporting its applicability in gene annotation for understudied organisms. We employed multiple methods to visualize the decision-making process of the model. The visualization results indicate that DRANetSplicer can learn and interpret well-known biological features, further validating its overall performance. Our study systematically examined and confirmed the predictive ability of DRANetSplicer from various levels and perspectives, indicating that its practical application in gene annotation is justified.

Conserved and divergent signals in 5' splice site sequences across fungi, metazoa and plants

In eukaryotic organisms the ensemble of 5' splice site sequences reflects the balance between natural nucleotide variability and minimal molecular constraints necessary to ensure splicing fidelity. This compromise shapes the underlying statistical patterns in the composition of donor splice site sequences. The scope of this study was to mine conserved and divergent signals in the composition of 5' splice site sequences. Because 5' donor sequences are a major cue for proper recognition of splice sites, we reasoned that statistical regularities in their composition could reflect the biological functionality and evolutionary history associated with splicing mechanisms. Results: We considered a regularized maximum entropy modeling framework to mine for non-trivial two-site correlations in donor sequence datasets corresponding to 30 different eukaryotes. For each analyzed species, we identified minimal sets of two-site coupling patterns that were able to replicate, at a given regularization level, the observed one-site and two-site frequencies in donor sequences. By performing a systematic and comparative analysis of 5'splice sites we showed that lineage information could be traced from joint di-nucleotide probabilities. We were able to identify characteristic two-site coupling patterns for plants and animals, and propose that they may echo differences in splicing regulation previously reported between these groups.

Alternative splicing signature of alveolar type II epithelial cells of Tibetan pigs under hypoxia-induced

Alternative splicing (AS) allows the generation of multiple transcript variants from a single gene and affects biological processes by generating protein diversity in organisms. In total, 41,642 AS events corresponding to 9,924 genes were identified, and SE is the most abundant alternatively spliced type. The analysis of functional categories demonstrates that alternatively spliced differentially expressed genes (DEGs) were enriched in the MAPK signaling pathway and hypoxia-inducible factor 1 (HIF-1) signaling pathway. Proteoglycans in cancer between the normoxic (21% O2, TN and LN) and hypoxic (2% O2, TL and LL) groups, such as SLC2A1, HK1, HK2, ENO3, and PFKFB3, have the potential to rapidly proliferate alveolar type II epithelial (ATII) cells by increasing the intracellular levels of glucose and quickly divert to anabolic pathways by glycolysis intermediates under hypoxia. ACADL, EHHADH, and CPT1A undergo one or two AS types with different frequencies in ATII cells between TN and TL groups (excluding alternatively spliced DEGs shared between normoxic and hypoxic groups), and a constant supply of lipids might be obtained either from the circulation or de novo synthesis for better growth of ATII cells under hypoxia condition. MCM7 and MCM3 undergo different AS types between LN and LL groups (excluding alternatively spliced DEGs shared between normoxic and hypoxic groups), which may bind to the amino-terminal PER-SIM-ARNT domain and the carboxyl terminus of HIF-1α to maintain their stability. Overall, AS and expression levels of candidate mRNAs between Tibetan pigs and Landrace pigs revealed by RNA-seq suggest their potential involvement in the ATII cells grown under hypoxia conditions.

Regulation of alternative splicing of PaFT and PaFDL1, the FT and FD homologs in Platanus acerifolia

Alternative splicing (AS) selects different alternative splice sites and produces a variety of transcripts with different exon/intron combinations, which may result in multiple protein isoforms. The splicing signals include cis-elements and RNA structures; however, the mechanisms of AS regulation in plants have yet to be elucidated. Previous studies have shown that in Platanus acerifolia, the FLOWERING LOCUS T (FT) homolog PaFT has a unique and complex AS pattern, in which most of the splice forms of PaFT involve the first and/or second intron, and the FD homolog PaFDL1 produces two transcripts via AS, whereas the other FT homolog PaFTL is not regulated by AS. In this study, the regulatory mechanism of the AS of PaFT was demonstrated to be conserved in different plant species. To define the distribution of the AS regulatory signals, the intron-swap, site-directed mutagenesis of alternative splice sites, and deletion experiment were performed. For the PaFT gene, all the signals that regulate the AS of the first intron were located within this intron, while the usage of the first alternative splice site in the second intron was determined by the first intron. Meanwhile, the AS of PaFDL1 might be co-regulated by exons and the first intron. Additionally, the first alternative splice site and adjacent region in PaFT intron 1 might contain cis-elements and/or RNA structures that affect the use of the other sites. This study had provided a deeper insight into the distribution of AS signals in plants, namely the AS signals of different splice sites might exist in the intron where the sites were present, and might also be distributed in exons or other introns.

RNA-seq Analysis Reveals Alternative Splicing Under Heat Stress in Rainbow Trout (Oncorhynchus mykiss)

Rainbow trout (Oncorhynchus mykiss) is one of the most economically important cold-water farmed species in the world, and transcriptomic studies in response to heat stress have been conducted and will be studied in depth. Alternative splicing (AS), a post-transcriptional regulatory process that regulates gene expression and increases proteomic diversity, is still poorly understood in rainbow trout under heat stress. In the present study, 18,623 alternative splicing events were identified from 9936 genes using RNA transcriptome sequencing technology (RNA-Seq) and genomic information. A total of 2731 differential alternative splicing (DAS) events were found among 2179 differentially expressed genes (DEGs). Gene ontology analysis revealed that the DEGs were mainly enriched in cellular metabolic process, cell part, and organic cyclic compound binding under heat stress. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis displayed that the DEGs were enriched for 39 pathways, and some key pathways, such as lysine degradation, are involved in the regulation of heat stress in liver tissues of rainbow trout. The results were validated by qRT-PCR, confirming reliability of our bioinformatics analysis.

Uncovering the impacts of alternative splicing on the proteome with current omics techniques

The high‐throughput sequencing of cellular RNAs has underscored a broad effect of isoform diversification through alternative splicing on the transcriptome. Moreover, the differential production of transcript isoforms from gene loci has been recognized as a critical mechanism in cell differentiation, organismal development, and disease. Yet, the extent of the impact of alternative splicing on protein production and cellular function remains a matter of debate. Multiple experimental and computational approaches have been developed in recent years to address this question. These studies have unveiled how molecular changes at different steps in the RNA processing pathway can lead to differences in protein production and have functional effects. New and emerging experimental technologies open exciting new opportunities to develop new methods to fully establish the connection between messenger RNA expression and protein production and to further investigate how RNA variation impacts the proteome and cell function.

This article is categorized under:

RNA Processing > Splicing Regulation/Alternative Splicing

Translation > Regulation

RNA Evolution and Genomics > Computational Analyses of RNA

Arabidopsis exoribonuclease USB1 interacts with the PPR-domain protein SOAR1 to negatively regulate abscisic acid signaling

Signaling by the phytohormone abscisic acid (ABA) involves pre-mRNA splicing, a key process of post-transcriptional regulation of gene expression. However, the regulatory mechanism of alternative pre-mRNA splicing in ABA signaling remains largely unknown. We previously identified a pentatricopeptide repeat protein SOAR1 (suppressor of the ABAR-overexpressor 1) as a crucial player downstream of ABAR (putative ABA receptor) in ABA signaling. In this study, we identified a SOAR1 interaction partner USB1, which is an exoribonuclease catalyzing U6 production for spliceosome assembly. We reveal that together USB1 and SOAR1 negatively regulate ABA signaling in early seedling development. USB1 and SOAR1 are both required for the splicing of transcripts of numerous genes, including those involved in ABA signaling pathways, suggesting that USB1 and SOAR1 collaborate to regulate ABA signaling by affecting spliceosome assembly. These findings provide important new insights into the mechanistic control of alternative pre-mRNA splicing in the regulation of ABA-mediated plant responses to environmental cues.

Reconstruction of the Carbohydrate 6-O Sulfotransferase Gene Family Evolution in Vertebrates Reveals Novel Member, CHST16, Lost in Amniotes

Glycosaminoglycans are sulfated polysaccharide molecules, essential for many biological processes. The 6-O sulfation of glycosaminoglycans is carried out by carbohydrate 6-O sulfotransferases (C6OSTs), previously named Gal/GalNAc/GlcNAc 6-O sulfotransferases. Here, for the first time, we present a detailed phylogenetic reconstruction, analysis of gene synteny conservation and propose an evolutionary scenario for the C6OST family in major vertebrate groups, including mammals, birds, nonavian reptiles, amphibians, lobe-finned fishes, ray-finned fishes, cartilaginous fishes, and jawless vertebrates. The C6OST gene expansion likely started early in the chordate lineage, giving rise to four ancestral genes after the divergence of tunicates and before the emergence of extant vertebrates. The two rounds of whole-genome duplication in early vertebrate evolution (1R/2R) only contributed two additional C6OST subtype genes, increasing the vertebrate repertoire from four genes to six, divided into two branches. The first branch includes CHST1 and CHST3 as well as a previously unrecognized subtype, CHST16 that was lost in amniotes. The second branch includes CHST2, CHST7, and CHST5. Subsequently, local duplications of CHST5 gave rise to CHST4 in the ancestor of tetrapods, and to CHST6 in the ancestor of primates. The teleost-specific gene duplicates were identified for CHST1, CHST2, and CHST3 and are result of whole-genome duplication (3R) in the teleost lineage. We could also detect multiple, more recent lineage-specific duplicates. Thus, the vertebrate repertoire of C6OST genes has been shaped by gene duplications and gene losses at several stages of vertebrate evolution, with implications for the evolution of skeleton, nervous system, and cell-cell interactions.

Modeling one thousand intron length distributions with fitild

Motivation

Intron length distribution (ILD) is a specific feature of a genome that exhibits extensive species-specific variation. Whereas ILD contributes to up to 30% of the total information content for intron recognition in some species, rendering it an important component of computational gene prediction, very few studies have been conducted to quantitatively characterize ILDs of various species.

Results

We developed a set of computer programs (fitild, compild, etc.) to build statistical models of ILDs and compare them with one another. Each ILD of more than 1000 genomes was fitted with fitild to a statistical model consisting of one, two, or three components of Frechet distributions. Several measures of distances between ILDs were calculated by compild. A theoretical model was presented to better understand the origin of the observed shape of an ILD.

Availability and implementation

The C++ source codes are available at https://github.com/ogotoh/fitild.git/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Architecture and Distribution of Introns in Core Genes of Four Fusarium Species

Removal of introns from transcribed RNA represents a crucial step during the production of mRNA in eukaryotes. Available whole-genome sequences and expressed sequence tags (ESTs) have increased our knowledge of this process and revealed various commonalities among eukaryotes. However, certain aspects of intron structure and diversity are taxon-specific, which can complicate the accuracy of in silico gene prediction methods. Using core genes, we evaluated the distribution and architecture of Fusarium circinatum spliceosomal introns, and linked these characteristics to the accuracy of the predicted gene models of the genome of this fungus. We also evaluated intron distribution and architecture in F. verticillioides, F. oxysporum, and F. graminearum, and made comparisons with F. circinatum. Results indicated that F. circinatum and the three other Fusarium species have canonical 5′ and 3′ splice sites, but with subtle differences that are apparently not shared with those of other fungal genera. The polypyrimidine tract of Fusarium introns was also found to be highly divergent among species and genes. Furthermore, the conserved adenosine nucleoside required during the first step of splicing is contained within unique branch site motifs in certain Fusarium introns. Data generated here show that introns of F. circinatum, as well as F. verticillioides, F. oxysporum, and F. graminearum, are characterized by a number of unique features such as the CTHAH and ACCAT motifs of the branch site. Incorporation of such information into genome annotation software will undoubtedly improve the accuracy of gene prediction methods used for Fusarium species and related fungi.

A Novel Dual Expression Platform for High Throughput Functional Screening of Phage Libraries in Product like Format

High throughput screenings of single chain Fv (scFv) antibody phage display libraries are currently done as soluble scFvs produced in E.coli. Due to endotoxin contaminations from bacterial cells these preparations cannot be reliably used in mammalian cell based assays. The monovalent nature and lack of Fc in soluble scFvs prevent functional assays that are dependent on target cross linking and/or Fc functions. A convenient approach is to convert scFvs into scFv.Fc fusion proteins and express them in mammalian cell lines for screening. This approach is low throughput and is only taken after primary screening of monovalent scFvs that are expressed in bacteria. There is no platform at present that combines the benefits of both bacterial and mammalian expression system for screening phage library output. We have, therefore, developed a novel dual expression vector, called pSplice, which can be used to express scFv.Fc fusion proteins both in E.coli and mammalian cell lines. The hallmark of the vector is an engineered intron which houses the bacterial promoter and signal peptide for expression and secretion of scFv.Fc in E.coli. When the vector is transfected into a mammalian cell line, the intron is efficiently spliced out resulting in a functional operon for expression and secretion of the scFv.Fc fusion protein into the culture medium. By applying basic knowledge of mammalian introns and splisosome, we designed this vector to enable screening of phage libraries in a product like format. Like IgG, the scFv.Fc fusion protein is bi-valent for the antigen and possesses Fc effector functions. Expression in E.coli maintains the speed of the bacterial expression platform and is used to triage clones based on binding and other assays that are not sensitive to endotoxin. Triaged clones are then expressed in a mammalian cell line without the need for any additional cloning steps. Conditioned media from the mammalian cell line containing the fusion proteins are then used for different types of cell based assays. Thus this system retains the speed of the current screening system for phage libraries and adds additional functionality to it.

Pre-mRNA Splicing in Plants: In Vivo Functions of RNA-Binding Proteins Implicated in the Splicing Process

Alternative pre-messenger RNA splicing in higher plants emerges as an important layer of regulation upon exposure to exogenous and endogenous cues. Accordingly, mutants defective in RNA-binding proteins predicted to function in the splicing process show severe phenotypic alterations. Among those are developmental defects, impaired responses to pathogen threat or abiotic stress factors, and misregulation of the circadian timing system. A suite of splicing factors has been identified in the model plant Arabidopsis thaliana. Here we summarize recent insights on how defects in these splicing factors impair plant performance.

Identification and Validation of Evolutionarily Conserved Unusually Short Pre-mRNA Introns in the Human Genome

According to the length distribution of human introns, there is a large population of short introns with a threshold of 65 nucleotides (nt) and a peak at 85 nt. Using human genome and transcriptome databases, we investigated the introns shorter than 66 nt, termed ultra-short introns, the identities of which are scarcely known. Here, we provide for the first time a list of bona fide human ultra-short introns, which have never been characterized elsewhere. By conducting BLAST searches of the databases, we screened 22 introns (37–65 nt) with conserved lengths and sequences among closely related species. We then provide experimental and bioinformatic evidence for the splicing of 15 introns, of which 12 introns were remarkably G-rich and 9 introns contained completely inefficient splice sites and/or branch sites. These unorthodox characteristics of ultra-short introns suggest that there are unknown splicing mechanisms that differ from the well-established mechanism.

The complex intron landscape and massive intron invasion in a picoeukaryote provides insights into intron evolution

Genes in pieces and spliceosomal introns are a landmark of eukaryotes, with intron invasion usually assumed to have happened early on in evolution. Here, we analyze the intron landscape of Micromonas, a unicellular green alga in the Mamiellophyceae lineage, demonstrating the coexistence of several classes of introns and the occurrence of recent massive intron invasion. This study focuses on two strains, CCMP1545 and RCC299, and their related individuals from ocean samplings, showing that they not only harbor different classes of introns depending on their location in the genome, as for other Mamiellophyceae, but also uniquely carry several classes of repeat introns. These introns, dubbed introner elements (IEs), are found at novel positions in genes and have conserved sequences, contrary to canonical introns. This IE invasion has a huge impact on the genome, doubling the number of introns in the CCMP1545 strain. We hypothesize that each IE class originated from a single ancestral IE that has been colonizing the genome after strain divergence by inserting copies of itself into genes by intron transposition, likely involving reverse splicing. Along with similar cases recently observed in other organisms, our observations in Micromonas strains shed a new light on the evolution of introns, suggesting that intron gain is more widespread than previously thought.

Assessment and refinement of eukaryotic gene structure prediction with gene-structure-aware multiple protein sequence alignment

Background

Accurate computational identification of eukaryotic gene organization is a long-standing problem. Despite the fundamental importance of precise annotation of genes encoded in newly sequenced genomes, the accuracy of predicted gene structures has not been critically evaluated, mostly due to the scarcity of proper assessment methods.

Results

We present a gene-structure-aware multiple sequence alignment method for gene prediction using amino acid sequences translated from homologous genes from many genomes. The approach provides rich information concerning the reliability of each predicted gene structure. We have also devised an iterative method that attempts to improve the structures of suspiciously predicted genes based on a spliced alignment algorithm using consensus sequences or reliable homologs as templates. Application of our methods to cytochrome P450 and ribosomal proteins from 47 plant genomes indicated that 50 ~ 60 % of the annotated gene structures are likely to contain some defects. Whereas more than half of the defect-containing genes may be intrinsically broken, i.e. they are pseudogenes or gene fragments, located in unfinished sequencing areas, or corresponding to non-productive isoforms, the defects found in a majority of the remaining gene candidates can be remedied by our iterative refinement method.

Conclusions

Refinement of eukaryotic gene structures mediated by gene-structure-aware multiple protein sequence alignment is a useful strategy to dramatically improve the overall prediction quality of a set of homologous genes. Our method will be applicable to various families of protein-coding genes if their domain structures are evolutionarily stable. It is also feasible to apply our method to gene families from all kingdoms of life, not just plants.

Global dissection of alternative splicing in paleopolyploid soybean

Alternative splicing (AS) is common in higher eukaryotes and plays an important role in gene posttranscriptional regulation. It has been suggested that AS varies dramatically among species, tissues, and duplicated gene families of different sizes. However, the genomic forces that govern AS variation remain poorly understood. Here, through genome-wide identification of AS events in the soybean (Glycine max) genome using high-throughput RNA sequencing of 28 samples from different developmental stages, we found that more than 63% of multiexonic genes underwent AS. More AS events occurred in the younger developmental stages than in the older developmental stages for the same type of tissue, and the four main AS types, exon skipping, intron retention, alternative donor sites, and alternative acceptor sites, exhibited different characteristics. Global computational analysis demonstrated that the variations of AS frequency and AS types were significantly correlated with the changes of gene features and gene transcriptional level. Further investigation suggested that the decrease of AS within the genome-wide duplicated genes were due to the diminution of intron length, exon number, and transcriptional level. Altogether, our study revealed that a large number of genes were alternatively spliced in the soybean genome and that variations in gene structure and transcriptional level may play important roles in regulating AS.

Systematic evaluation of spliced alignment programs for RNA-seq data

Authors compare RNA-seq aligners on mouse and human data sets using benchmarks such as alignment yield, splice junction accuracy and suitability for transcript reconstruction. The work highlights the strength of each program and discusses outstanding needs in RNA-seq analysis.

High-throughput RNA sequencing is an increasingly accessible method for studying gene structure and activity on a genome-wide scale. A critical step in RNA-seq data analysis is the alignment of partial transcript reads to a reference genome sequence. To assess the performance of current mapping software, we invited developers of RNA-seq aligners to process four large human and mouse RNA-seq data sets. In total, we compared 26 mapping protocols based on 11 programs and pipelines and found major performance differences between methods on numerous benchmarks, including alignment yield, basewise accuracy, mismatch and gap placement, exon junction discovery and suitability of alignments for transcript reconstruction. We observed concordant results on real and simulated RNA-seq data, confirming the relevance of the metrics employed. Future developments in RNA-seq alignment methods would benefit from improved placement of multimapped reads, balanced utilization of existing gene annotation and a reduced false discovery rate for splice junctions.

Complexity of the alternative splicing landscape in plants

Alternative splicing (AS) of precursor mRNAs (pre-mRNAs) from multiexon genes allows organisms to increase their coding potential and regulate gene expression through multiple mechanisms. Recent transcriptome-wide analysis of AS using RNA sequencing has revealed that AS is highly pervasive in plants. Pre-mRNAs from over 60% of intron-containing genes undergo AS to produce a vast repertoire of mRNA isoforms. The functions of most splice variants are unknown. However, emerging evidence indicates that splice variants increase the functional diversity of proteins. Furthermore, AS is coupled to transcript stability and translation through nonsense-mediated decay and microRNA-mediated gene regulation. Widespread changes in AS in response to developmental cues and stresses suggest a role for regulated splicing in plant development and stress responses. Here, we review recent progress in uncovering the extent and complexity of the AS landscape in plants, its regulation, and the roles of AS in gene regulation. The prevalence of AS in plants has raised many new questions that require additional studies. New tools based on recent technological advances are allowing genome-wide analysis of RNA elements in transcripts and of chromatin modifications that regulate AS. Application of these tools in plants will provide significant new insights into AS regulation and crosstalk between AS and other layers of gene regulation.

Gene make-up: rapid and massive intron gains after horizontal transfer of a bacterial α-amylase gene to Basidiomycetes

BACKGROUND

Increasing genome data show that introns, a hallmark of eukaryotes, already existed at a high density in the last common ancestor of extant eukaryotes. However, intron content is highly variable among species. The tempo of intron gains and losses has been irregular and several factors may explain why some genomes are intron-poor whereas other are intron-rich.

RESULTS

We studied the dynamics of intron gains and losses in an α-amylase gene, whose product breaks down starch and other polysaccharides. It was transferred from an Actinobacterium to an ancestor of Agaricomycotina. This gene underwent further duplications in several species. The results indicate a high rate of intron insertions soon after the gene settled in the fungal genome. A number of these oldest introns, regularly scattered along the gene, remained conserved. Subsequent gains and losses were lineage dependent, with a majority of losses. Moreover, a few species exhibited a high number of both specific intron gains and losses in recent periods. There was little sequence conservation around insertion sites, then probably little information for splicing, whereas splicing sites, inside introns, showed typical and conserved patterns. There was little variation of intron size.

CONCLUSIONS

Since most Basidiomycetes have intron-rich genomes and this richness was ancestral in Fungi, long before the transfer event, we suggest that the new gene was shaped to comply with requirements of the splicing machinery, such as short exon and intron sizes, in order to be correctly processed.

Benchmarking spliced alignment programs including Spaln2, an extended version of Spaln that incorporates additional species-specific features

Spliced alignment plays a central role in the precise identification of eukaryotic gene structures. Even though many spliced alignment programs have been developed, recent rapid progress in DNA sequencing technologies demands further improvements in software tools. Benchmarking algorithms under various conditions is an indispensable task for the development of better software; however, there is a dire lack of appropriate datasets usable for benchmarking spliced alignment programs. In this study, we have constructed two types of datasets: simulated sequence datasets and actual cross-species datasets. The datasets are designed to correspond to various real situations, i.e. divergent eukaryotic species, different types of reference sequences, and the wide divergence between query and target sequences. In addition, we have developed an extended version of our program Spaln, which incorporates two additional features to the scoring scheme of the original version, and examined this extended version, Spaln2, together with the original Spaln and other representative aligners based on our benchmark datasets. Although the effects of the modifications are not individually striking, Spaln2 is consistently most accurate and reasonably fast in most practical cases, especially for plants and fungi and for increasingly divergent pairs of target and query sequences.

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.

Compensatory signals associated with the activation of human GC 5' splice sites

GC 5′ splice sites (5′ss) are present in ∼1% of human introns, but factors promoting their efficient selection are poorly understood. Here, we describe a case of X-linked agammaglobulinemia resulting from a GC 5′ss activated by a mutation in BTK intron 3. This GC 5′ss was intrinsically weak, yet it was selected in >90% primary transcripts in the presence of a strong and intact natural GT counterpart. We show that efficient selection of this GC 5′ss required a high density of GAA/CAA-containing splicing enhancers in the exonized segment and was promoted by SR proteins 9G8, Tra2β and SC35. The GC 5′ss was efficiently inhibited by splice-switching oligonucleotides targeting either the GC 5′ss itself or the enhancer. Comprehensive analysis of natural GC-AG introns and previously reported pathogenic GC 5′ss showed that their efficient activation was facilitated by higher densities of splicing enhancers and lower densities of silencers than their GT 5′ss equivalents. Removal of the GC-AG introns was promoted to a minor extent by the splice-site strength of adjacent exons and inhibited by flanking Alu repeats, with the first downstream Alus located on average at a longer distance from the GC 5′ss than other transposable elements. These results provide new insights into the splicing code that governs selection of noncanonical splice sites.