Abstract
The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the “dark matter” of the eukaryotic genome and discuss some of the outstanding questions in the field.
Introduction
Historically, introns were considered the non-coding, non-functional sequence elements which disrupt those that are protein-coding, called exons (Gilbert, 1978). While this protein-centric definition of introns (Figure 1, left) has served its purpose, their presence in long non-coding RNA reveals that introns are not specific to protein-coding genes but instead serve a broader role in eukaryotic gene expression (Krchňáková et al., 2019; ). Moreover, introns have been found to host other lariat-derived RNAs, including microRNAs, long noncoding RNAs, small nucleolar RNAs, small nuclear RNAs, and circular RNAs that are crucial for gene regulation (Liu and Maxwell, 1990; Hesselberth, 2013; Seal et al., 2020; Kumari et al., 2022; Vakirlis et al., 2022). Introns can also house enhancer elements that drive tissue-specific expression kinetics during complex vertebrate development and embryogenesis (Emera et al., 2016; ; Meng et al., 2021; Shiau et al., 2022). These intervening sequences necessitated co-evolution of splicing machinery to facilitate production of a contiguous transcript capable of encoding a functional unit (Grabowski et al., 1985; Nilsen, 2003). Inhibition of splicing results in retention of introns in the mature transcript, which often disrupts the open reading frame and ultimately dictates the fate of the final transcript (Kaida et al., 2007; Effenberger et al., 2017; Olthof et al., 2021). Since the discovery of splicing, introns have been extensively investigated and the significance of splicing in regulating gene expression is well documented (Singh and Padgett, 2009; Tellier et al., 2020; Zhu et al., 2020; ; Reimer et al., 2021). Taken together, the presence of introns has a significant impact on eukaryotic gene expression and underpins many of the complexities required to build higher eukaryotes. Therefore, here we present an intron-centric perspective (Figure 1, right) towards understanding regulation of eukaryotic gene expression.
FIGURE 1
Function and evolution of intronic elements
Introns date back to the last eukaryotic common ancestor, after invasion into the early eukaryotic genome (Russell et al., 2006; ; ). While an endogenous model has been proposed to explain the emergence of introns (), there is a general consensus that prokaryotic group II self-splicing introns underwent invasion and mutational degeneration during early eukaryogenesis, resulting in inert introns and trans-acting splicing machinery (Michel et al., 1989; Sharp, 1991; Sontheimer et al., 1999; Shukla and Padgett, 2002). As the origin of eukaryotic introns has been extensively described (Koonin, 2006; Rogozin et al., 2012; Vosseberg and Snel, 2017; ; Smathers and Robart, 2019), here, we focus on the continued maintenance and diversification of introns in eukaryotic genomes.
Prokaryotic group II self-splicing introns behaved largely as transposable elements, which may have facilitated their invasion of the eukaryotic genome (Figure 2) (Lambowitz and Zimmerly, 2011). Initially characterized in the maize genome, transposable elements are repetitive sequences found across eukaryotes and are critically known for their ability to relocate in the genome and alter gene expression (McClintock, 1950; SanMiguel et al., 1996; Elliott et al., 2005; Wells and Feschotte, 2020). Short and long interspersed retro-transposable elements (SINEs/LINEs) belong to the non-long terminal repeat class of elements which have retained transposable activity and are highly represented in the human genome as Alu and L1 elements, respectively (Kazazian and Moran, 1998; Lander et al., 2001; ). When carrying splice sites, these transposable elements can create novel exon/intron boundaries, which hold the potential to alter expression of that gene (Figure 2); a detailed description of exon/intron boundaries and their recognition by splicing machinery is discussed in the following sections. For example, a recent study queried pathogenic mutations that were associated with novel intron-exon boundaries in humans and identified those which aligned with transposable elements. They found that clusters of transposable elements are more liable to exonization, likely due to the combined effort of LTR and Alu elements in potentiating all necessary splice sites (). In another computational investigation of the human genome, mutagenesis of Alu elements into weak splice sites was found to be well-tolerated if not retained long-term and was often associated with exon skipping events (Sorek et al., 2002). Exon skipping is a frequently observed form of alternative splicing, which more broadly serves as an important regulatory node for gene expression in developing systems (). One can then speculate that Alu elements in this manner allow for transient sampling of novel functions of proteins encoded by these alternatively spliced transcripts. This idea is an extension of the already known role of Alu elements in tissue-specific transcription regulation (Franchini et al., 2011). Notably, weak splice sites in Alu elements can eventually become constitutively spliced exons, losing their capacity for transposition and become exons used in regulating tissue-specific gene expression, as is observed in the human NARF gene (Lev-Maor et al., 2007).
FIGURE 2
Inherent to the jumping nature of transposable elements is the impartiality of transposon landing. Transposon insertion would likely be deleterious in the protein-coding region of a gene, leading to evolutionary selection against that gene configuration. However, in a heterozygote, transposon-induced activation of a novel splice site within an intron could allow for a low-cost trial of differentially spliced isoforms, while still maintaining a functionally expressed copy. A susceptibility of spliceosomal introns to genomic recombination was demonstrated in two Saccharomyces cerevisiae genes, RPL8B and ADH2. Truncated versions of these genes were used in a splicing reporter construct, such that the second exon was expressed in frame with a fused EGFP cassette. Additionally, each construct carried an embedded S. pombe his5+ gene within the first intron, encoded for in the opposite direction as EGFP. Here, the his5+ gene contains an artificial intron lacking a catalytic branch point, and containing splice sites in such an orientation that they are only capable of splicing from the EGFP transcript. Thus, splicing of the artificial intron followed by transposition of the EGFP intron into the genomic loci was required to confer a positive result (Lee and Stevens, 2016). Meanwhile, Gozashti et al. (2022), has attributed rapid, lineage-specific intron gains to Introner elements derived from transposable elements. Through analysis of 1,700 species, these “intron-generating transposable elements families” were identified in approximately 5% of genomes and significantly overrepresented in aquatic lineages. Based on statistical association models and a consideration of likely propagation mechanisms, they concluded that Introner elements may facilitate recent intron gain, particularly through horizontal gene transfer in aquatic lineages. The activity of Introner elements is particularly interesting, as mechanisms of Introners in Micromonas pusilla and Aureococcus anophagefferens exhibit seemingly preferential insertion between pre-existing nucleosomes (Huff et al., 2016). The rationale here is such that the linker sequence between nucleosomes is often open and available for insertion events. Further support for this idea is seen in the unequal distribution and position of nucleosomes observed between protein-coding exons, pseudo exons, and introns in human and Caenorhabditis elegans (). Using transcriptomic and genomic sequencing data, Huff et al. (2016), reported that Introners are largely capable of co-opting splice sites and inserting by DNA transposition in both orientations, though with biases consistent with species-specific patterns in genome organization. Outside of splice site generation, transposons have also been implicated in regulation of splicing-competent snRNAs, such as those L1 transposons which are associated with formation of U6 pseudogenic snRNA (Doucet et al., 2015). Pseudogenes can encode variations of spliceosomal snRNAs, the implications of which are discussed further below. In all, transposable elements further expand gene structure by modifying intronic elements, thus revealing a critical role of non-coding intronic elements in eukaryotic genome evolution.
Classification and splicing of introns
After the discovery of splicing, identified introns appeared to show a pattern of conserved terminal di-nucleotides at the exon-intron and intron-exon boundaries, and this feature became a defining characteristic of spliced introns (; ; ). As sequencing techniques have progressed and data now includes more diverse eukaryotic genomes, it is increasingly clear that introns are defined by several extended consensus sequences. These include the 5′ splice site (5′SS), the branch point sequence (BPS), and the 3′ splice site (3′SS) (Dietrich et al., 1997; Mercer et al., 2015). Not long after their discovery, it was determined that most introns are processed by five Uridylyl-rich snRNAs—U1, U2, U4, U5, and U6—that are highly conserved between eukaryotes and assemble into a ribonucleoprotein complex, the spliceosome (; ; ; Nilsen, 2003; Wahl et al., 2009). Specifically, U1 snRNA has complementarity at the 5′ splice site, marking the exon-intron boundary, while U2 snRNA base pairs around a conserved adenosine toward the 3′ end, at what has become known as the branch point sequence (Yan and Ares, 1996; Malca et al., 2003). The direct base pairing of these snRNAs with splice site consensus sequences helps to recognize and remodel the intron during splicing, conferring the core function of the spliceosome.
As this mechanism was coming into focus, Jackson (1991), discovered spliced transcripts, that when mapped to the genome, showed intronic splice site sequences that were incompatible with the identified snRNAs. The fact that these introns were nonetheless spliced suggested the existence of a separate mechanism for their removal. This discordant finding led to sequence-based investigations for U snRNAs with complementary to non-consensus splice sites. This included an exploratory genomics investigation by Hall and Padgett (1994), and ultimately led to the hypothesis that newly identified U11 and U12 snRNAs serve in roles analogous to U1 and U2 during splicing (Montzka and Steitz, 1988). A role for U11 and U12 was confirmed in vitro (Tarn and Steitz, 1996a) and in vivo (Hall and Padgett, 1996; Kolossova and Padgett, 1997), and bolstered by the additional identification of snRNAs analogous to U4/U6, U4atac and U6atac (Tarn and Steitz, 1996b; Incorvaia and Padgett, 1998). Based on their relative abundance in analyzed genomes, the intron types and their respective spliceosomes were henceforward labeled major (U2-type) and minor (U12-type) in those eukaryotes that maintain them in parallel (; Lynch and Richardson, 2002; Lin et al., 2010). Of note, major introns and the major spliceosome are ubiquitous in the eukaryotic lineage, while minor introns and the minor spliceosome are reportedly absent in some lineages, such as Caenorhabditis elegans ().
Both the major and minor spliceosomes employ U5 snRNA, and each snRNA further associates with specific proteins in their splicing-competent forms (Tarn and Steitz, 1996a; Tarn and Steitz, 1997). Though the individual snRNAs have specific proteins associated with their regulation and maturation, many of the remaining proteins that comprise the spliceosome are shared between both the major and minor molecular machineries (Will et al., 1999); for a more comprehensive presentation of individual spliceosome components, see Olthof et al. (2022). Worth noting, the same protein can carry out different roles in each spliceosome, as is observed by URP (also called ZRSR2) (Tronchère et al., 1997; Shen et al., 2010). While the size and dynamic composition of the spliceosome can make it difficult to fully resolve, identifying the proteins involved in splicing regulation remains an area of active investigation. Recent biochemical and cryogenic electron microscopy investigations to this end have significantly enhanced our understanding of minor spliceosome-specific proteins. For example, the protein compositions of U4.U6/U5 and U4atac.U6atac/U5 tri-snRNP complexes were previously thought to be identical. However, co-immunoprecipitation and co-migration analyses have suggested that CENATAC may aid in 5′SS recognition for a subclass of minor introns characterized by AT-AN terminal di-nucleotides. Previously known as CCDC84, CENATAC was renamed following its mutagenic link to intron retention in human genes that contribute to chromosome stability and segregation (). Interestingly, phylogenetic profiling of CENATAC across 90 eukaryotic species showed that it co-enriched with other components of the minor spliceosome, including the newly characterized SCNM1 protein (). The U12 snRNA is flanked by the N-terminal C2H2 zinc fingers of SCNM1, which interacts with the U12/BPS duplex and the U12 Sm ring (). The N-terminus of SCNM1 also functions to stabilize U6atac and RNF113A at the 5′SS, maintenance of which is required for spliceosome activation in vivo (Incorvaia and Padgett, 1998; ). Structural insights were also important in identifying the novel minor spliceosome protein, RBM48, which is now known to bind ARMC7 and interact with terminal ends of U6atac snRNA via conserved RNA binding residues (; Siebert et al., 2022). Structural analyses of the minor spliceosome are a recent advancement and do not yet cover all phases of splicing, notably excluding the U11/U12 di-snRNP. As such, there remains the possibility for other unidentified components regulating the nuances of minor intron splicing.
A delineation between major versus minor intron splicing is often based on the quantitative analysis of splice site conservation, and thus relative splice site strength. Intron splice sites are generally scored based on the degree of similarity to the major versus minor intron consensus sequences found in Figure 3, using position weight matrices (Sheth et al., 2006; ; Olthof et al., 2019; Moyer et al., 2020). The resultant major or minor intron classification inherently dictates how we interpret its processing, such that bioinformatically classified minor introns are predicted to be spliced by the minor spliceosome, and vice versa. However, RNA sequencing data has shown that, upon inhibition of the minor spliceosome, not all bioinformatically classified minor introns show a splicing defect (Olthof et al., 2019). Thus, parallel existence of major and minor spliceosomes, combined with diverging intron consensus sequences, reveal an added complexity in the relationship between a given intron and its recruited spliceosome. Akin to how the concept of a single intron type was disrupted by the discovery of minor introns; it seems increasingly likely that the binary classification or major versus minor itself is insufficient to fully resolve all introns. Rather, evidence has begun to suggest that the stringency of the classification schema fails to consider the fluidity of exons and introns. For example, use of novel splice sites within exonic regions in the unicellular Paramecium is evidence of intronization activity in eukaryotes (Ryll et al., 2022). In essence, these findings increasingly suggest that the current approach to intron classification is too reductive to fully capture the complexities and dynamic regulation of eukaryotic introns. Towards this end, an examination of minor-type splice sites in Physarum polycephalum has suggested that minor introns may exist in divergent, if not degenerative, types (Larue et al., 2021) and this idea is currently being refined in other studies that combine principles of speciation and comparative genomics.
FIGURE 3
How gene architecture informs splice site selection
Spliceosomal introns are known to range from tens of base pairs in length to hundreds of kilobases in length, with a mean length that is smaller in lower eukaryotes and larger in higher eukaryotes (Sakharkar et al., 2004; Piovesan et al., 2015; ; Li et al., 2017; Jakt et al., 2022). The size of an intron has an inherent impact on gene expression, as it will take longer for transcription machinery to create nascent transcripts. In turn, this will impact the kinetics of co-transcriptional intron splicing; these ideas have been explored in depth (Herzel et al., 2017; Wallace and Beggs, 2017; Neugebauer, 2019). It is long since established that relative intron and exon lengths can differentially affect splicing efficiency due to a presence or absence of regulatory elements and differing requirements for catalysis (Fox-Walsh et al., 2005; Kandul and Noor, 2009; Pai et al., 2017). Splicing efficiency refers to the proportion of spliced versus un-spliced transcripts relative to the number of total transcripts. This is commonly assessed using computational strategies that characterize splice events in the transcriptome (; Jiang et al., 2023), followed by a validation of observed changes in expression using techniques such as RT-PCR. In one assessment of how splicing efficiency and gene expression patterns may be coupled, intron length was found to contribute to the temporal coordination that is required for co-expression of genes with interdependent biochemical functions (Keane and Seoighe, 2016). This idea is further reflected by distinct differences in splice site strength relative to intron length, and by differences in splicing efficiency and mRNA abundance relative to gene length (Gelfman et al., 2011; Sánchez-Escabias et al., 2022). Vertebrates are known to increase splicing efficiency around longer introns via cell-specific recursive splicing and transposable elements that form stems with intronic RNA loops to juxtapose splice sites (Shepard et al., 2009; Zhang et al., 2018). For details on recursive splicing, please see published reviews (Georgomanolis et al., 2016; Gehring and Roignant, 2021; Joseph et al., 2022; Pitolli et al., 2022).
Separate from this, longer introns may also have a propensity to contain multiple splice sites within one intronic feature, leading to alternative splicing from competing splice site use (Sun and Chasin, 2000; Roca et al., 2003; Kapustin et al., 2011). Meaning it becomes increasingly likely that multiple splice sites be present, in addition to exonic splicing enhancers and silencer elements, which themselves can act as determinants of splice site usage (; Wang et al., 2006). It thus follows that the sequence content of the intron to be excised can drive splicing progression. Splice site selection is thought to occur by competing intron- and exon-definition models, which describe how the spliceosome assembles either through cross-bridging interactions across the intron itself or across the flanking exon. Specifically, the intron-definition model refers to the mechanism whereby 3′ SS selection is informed by recognition of the upstream 5′SS, such that the spliceosome assembles across the intron. For exon-definition interactions, 3′ SS recognition depends instead on recognition of the downstream 5′SS (Robberson et al., 1990; ; Romfo et al., 2000; ; Olthof et al., 2021). For example, most genes in Saccharomyces cerevisiae, contain only one, short intron. With this gene architecture, it is not surprising that intron-definition interactions predominate. Surprisingly, cryo-electron microscopy structures of the pre-catalytic spliceosome demonstrated that the same splicing machinery can perform exon-definition interactions in multi-intronic genes (Li et al., 2019). This finding brings to bear uncertainty as to how and when an intron- versus exon-centric model is utilized. This becomes especially important in higher vertebrates which have a larger intronic burden.
Reconciliation between the intron- and exon-definition models is coupled with new insight into how proximity rules inform splice site selection. Based on the length of an intron, the intron-centric proximity rule dictates a preference for the spliceosome to assemble over a splice site pair that minimizes the distance between 5′ and 3′ end selection (Reed and Maniatis, 1986). More recently, computational analyses by refined the exon-centric proximity rule, by which splice sites are selected to minimize the exon-spanning distance. Meaning if one were to imagine an intron with two adjacent sets of 5′ and 3′ splice sites, the intron-centric proximity rule would employ the innermost splice site pair, maximizing the resultant exons. Meanwhile, the exon-centric proximity rule would, in contrast, use the exon-proximal splice sites to maximize the size of the intron being excised. In either case, commitment to the intron-centric or exon-centric proximity rule has commensurate intronization/exonization consequences as molecular machinery decides whether to select for the smaller or larger exonic sequences. In addition to intron size, studies suggest that GC content of the intron may also be a determinative factor in the mechanism employed for splice site selection. In one study, (Tammer et al. (2022), examined the nucleotide composition of exons versus introns and subsequently identified genes they refer to as “differential” or “leveled”. In “leveled” genes, GC content is found to be similarly high in exons and introns, while “differential” genes are ones wherein GC content is low in exons, and even lower in introns. Notably, Tammer et al. (2022), describe a partiality for intron-definition interactions across “leveled” genes, while exon-definition interactions predominate over “differential” genes. This finding is in line with previously reported links between differential GC content and splice site selection ().
Spliceosomal sRNAs
As described above, snRNAs confer the primary function of the spliceosome through formation of specific base pair interactions with consensus sequences in the intron. The presence and function of snRNAs is essential for recognition and restructuring of the nascent mRNA transcript in the sequential, exothermic transesterification reactions that constitute splicing.
In general, most snRNAs (U1, U2, U4, U5, U11, and U12) are transcribed by RNA polymerase II, while U6 and U6atac expression are largely dependent on RNA polymerase III (Reddy et al., 1987; Jawdekar and Henry, 2008; Younis et al., 2013). Initiation of transcription of these snRNAs is highly reliant on the proximal and distal sequence elements located upstream of the snRNA-encoding region. Specifically, because they serve as promoter and enhancer elements for recruitment of transcription machinery through interactions with the SNAPc transcription factor complex and stabilizing co-activators (Sadowski et al., 1993; Henry et al., 1998; Mittal et al., 1999; ). Structural insights by cryogenic electron microscopy of SNAPc during the transcription of U6 has revealed the importance of conserved subunits which recognize and bind the proximal sequence element (Sun et al., 2022). One unique exception to this rule is for the expression of human U4atac snRNA, which is embedded into an intron of CLASP1 (Edery et al., 2011). Therefore, U4atac expression relies on RNA polymerase II mediated transcription of this gene, as well as successful splicing of this intron.
Within the genome, spliceosomal snRNAs often exist both as gene copies and gene families, whereby divergent genes can encode for variant snRNAs with nucleotide polymorphisms (; ). There are both productive and unproductive variants of the snRNAs annotated; productive snRNAs are capable of splicing, while those that are not are termed pseudogenic (Mabin et al., 2021). For example, the U6 snRNA has many pseudogenes and fewer productive copies that are dispersed throughout the genome, whereas U1 and U2 snRNAs are encoded by many functional copies that are organized in homogenous repeats (Van Arsdell and Weiner, 1984; Theissen et al., 1985; Tichelaar et al., 1998; Domitrovich and Kunkel, 2003; O’Reilly et al., 2013; ). The presence of multiple gene copies may in part explain the splicing-independent roles of U1 and U2 in regulating transcription termination and 3′ end processing (Friend et al., 2007; ; So et al., 2019). Moreover, the idea that multiple gene copies exist for minor spliceosomal snRNAs, including U4atac and U6atac, warrants further investigation. Even if multiple gene copies do exist, it must be noted that U6atac expression is maintained at a lower level through rapid post-transcriptional turnover (Younis et al., 2013).
Perhaps counterintuitively, U5 snRNA has the smallest gene family, yet it is the only shared snRNA between the major and minor spliceosomes. Investigations by Mabin et al. (2021) into the relevance of snRNA variants in splicing led to the discovery of high sequence identity between U5 variants. In fact, they report several U5 variants with a conserved stem consensus sequence (CUUUU) that can be incorporated into catalytic spliceosomes. Based on these observations, it has been suggested that U5 may not have a canonical snRNA; rather, specific variants may be optimal for use in one spliceosome type over the other (Mabin et al., 2021). While mechanistically unvalidated, this logic is consistent with the analogous nature of the other major versus minor snRNAs. Yet, it also remains possible that these U5 variants are regulated in a context-dependent way, as is observed for U1 snRNA variants during human stem cell programming (Vazquez-Arango et al., 2016). Additionally, U5 snRNA variants have been identified in regulating development in humans, Drosophila, and Lytechinus variegatus (Sontheimer and Steitz, 1992; Morales et al., 1997; ). The expression of snRNA variants to specify a differentiating transcriptome is not unique to U5 snRNA, but more broadly detected for other snRNAs and across species (Lo and Mount, 1990; ; Sierra-Montes et al., 2005; O’Reilly et al., 2013; Lu and Matera, 2014).
Functional sequence variants of the snRNAs have the potential to contact cryptic or degenerating splice sites, make novel protein interactions, and adopt secondary structures that alter spliceosome conformation. It is thus possible, given our evolving understanding of consensus sequences, that these variant snRNAs do confer complementarity to specific intron splice sites. Accordingly, from an intron-centric perspective, we must allow for the possibility that seemingly unproductive snRNAs are leveraged to splice a specific subset of introns. A role for non-consensus intron classes was voiced by Hudson et al. (2019), whose bioinformatics analyses of diplomonad and parabasalid lineage eukaryotes revealed splice site sequences that diverged from both the major and minor consensus sequences. They similarly identify divergent snRNAs, though they maintained key functional structures including stem loops and putative Sm binding sites. Perhaps more compelling, the discovered snRNAs showed aggregate features of both the major- and minor-type snRNAs, suggesting a propensity for the spliceosome to adopt complementarity to trans-spliced introns.
It remains to be established if variant snRNAs are evolutionarily selected for use in differential splicing or if they arise stochastically. Though, one could imagine that selective use of a variant splicing component would provide an opportunity to splice novel or divergent splice site sequences. It is known that mutations in the snRNAs can have pathogenic effects, as demonstrated by RNU12 which is causal to early onset cerebellar ataxia (Elsaid et al., 2017). Additionally, snRNA secondary structure is important for splicing as it dictates the RNA-protein interactions necessary for spliceosome assembly. For example, U11/U12-65K binds the 3′ stem loop II (SLII) of U12 snRNA based on the integrity of this structure and its RNA binding motif. Further, 3’ truncation mutants that disrupt the U12 SLIII are targeted for degradation by the nuclear exosome targeting complex upon reimport to the nucleus (Norppa and Frilander, 2021). In another example, the U2/U6 and U12/U6atac complexes are remodeled and stabilized prior to the first catalytic step in splicing by intramolecular base pairing with RBM22 (). Regardless, developmentally regulated snRNA variants demonstrate that mutations outside of critical structures may maintain, albeit differential, functionality. Thus, it stands to reason that variant snRNAs without disease-causing consequences to splicing may have a context-dependent role in the regulation of introns with divergent consensus sequences.
The evolutionary advantages of introns
Introns have served a valuable evolutionary role for eukaryotes in that they are more prone to genetic drift compared to exons. Introns appear to be under weaker selection than exons in somatic cells, which may be due to a mismatch repair system employed for exons that is notably lacking for introns (Hoffman and Birney, 2006; Resch et al., 2007; Frigola et al., 2017; Rodriguez-Galindo et al., 2020). Using a combinatorial multi-omics approach, Huang et al. (2018), has attributed the selective protection and mismatch repair of actively transcribed genes to an enrichment of H3K36me markers, which help regulate molecular responses to DNA damages induced by prolonged euchromatic conformation. More broad analyses of the differentiating human methylome reveal distinct differences in methylation pattern between genomics features, such that methylation is generally more common to exons than introns (Laurent et al., 2010). This unequal distribution may explain the higher frequency of mismatch repair observed for exons versus introns. In this capacity, introns can essentially act as a sponge to harbor mutations that would be otherwise detrimental in exonic sequences. Nevertheless, many mutations in intronic elements are linked to diseases, suggesting that there are limits to the number of mutations an intron can withstand. Mutations at splice sites and within introns are known to underscore an array of genetic and developmental disorders, including muscular dystrophy (Dominov et al., 2019) and inherited retinal diseases (Qian et al., 2021). Pathogenic disorders due to mutation of the spliceosome, i.e., spliceosomopathies, include but are not limited to craniofacial defects, myelodysplastic syndromes, and retinitis pigmentosa (Griffin and Saint-Jeannet, 2020). For review of major and minor splicing-associated diseases, see (; Olthof et al., 2022).
While introns are seemingly advantageous, prokaryotes show that the absence of introns is not prohibitive to life. This begs the question, to what extent do eukaryotic cells really require introns? In one study, Parenteau et al. (2008), investigated the consequences of intron depletion in Saccharomyces cerevisiae (Figure 4A). Introns are far less abundant in S. cerevisiae compared to other species, such as vertebrates and land plants, making the yeast genome a strong model for intron depletion studies (). Indeed, S. cerevisiae could survive without introns, however, intron-depleted strains fared variably when subjected to drug-induced and carbon source stresses. However, transcription machinery was found to be capable of responding to expression deficits following intron-depletion by using alternative promoter selections, highlighting the role introns play in expanding the eukaryotic transcriptome (Parenteau et al., 2008). Should one suppose that introns can be leveraged to induce stress-related patterns of gene expression, it then follows that the splicing efficiency of an intron is responsive to stress application. This idea was recently explored by Frumkin et al. (2019), who employed YFP reporter constructs containing known introns with high and low splicing efficiencies embedded and fused to a kanamycin resistance cassette (Figure 4B). To test the capacity of introns and the spliceosome to respond to metabolic pressure, the constructs were expressed in S. cerevisiae cells under antibiotic selection and subjected to a lab-evolution paradigm. Growth and transcriptomic analyses of derived cell generations revealed independent, adaptive mutations occurring both cis- and trans-to improve splicing efficiency and thus antibiotic resistance and cell survivability. The cis-mutations were proposed to increase accessibility of splice site sequences, while trans-mutations might increase the cellular abundance of splicing machinery. Importantly, cis-fitness-inducing mutations could alleviate selection-independent splicing inefficiencies, however, mutations in trans-were particularly advantageous during periods of active selection (Frumkin et al., 2019). Though these experiments were performed in S. cerevisiae, one can imagine that similar mechanisms may be employed for evolutionary adaptation. For example, in ecotypic Cichlid fish, alternative splicing is a dominant mechanism for rapid changes in gene expression. Specifically, alternative splicing underpins the diversification of jaw morphology as it relates to the food they have evolved to consume in different ecological niches (Singh et al., 2017).
FIGURE 4
The influence of introns on gene expression
In both mammals and plants, the presence of introns is known to enhance gene expression in a phenomenon sometimes referred to as intron-mediated enhancement (; Furger et al., 2002; Samadder et al., 2008). The recent development of sequencing techniques such as GRO-seq, mNET-seq and long read sequencing have revealed that splicing of neither major nor minor introns occurs in isolation, but rather in a highly active genomic context where splicing and transcription are coupled both kinetically and physically (Nojima et al., 2015; Herzel et al., 2017; Sheridan et al., 2019; Drexler et al., 2020; Reimer et al., 2021; Zhang et al., 2021). In the context of splicing informing transcription, the position of the intron matters, as promoter-proximal introns are especially known to enhance transcription (Furger et al., 2002; Rose et al., 2008). The knowledge that introns may enhance transcriptional output was leveraged to modify the generally used CMV promoter for expression plasmids, whereby introduction of an intron significantly upregulated transcription of downstream coding sequences (Simari et al., 1998).
The mechanism by which 5’ introns regulate transcription involves, at least in part, control of the open chromatin signatures H3K4me3 and H3K9ac, which facilitate recruitment of RNA polymerase II and general transcription factors to promoters. These marks are deposited at the first exon-intron boundary of genes, explaining how the distance between transcription start site and the first intron can influence the expression level of a gene (, Lister, 2009). Interestingly, differential methylation patterns are not unique to protein-coding genes, as revealed through a bioinformatics model which considered the modified human nucleosome library and analysis of splicing efficiency. For example, high nucleosome density was observed in the internal exons of long non-coding RNAs, while high H3K4me3 signals were observed in upstream introns. Importantly, these signatures were often associated with exon skipping and intron retention, particularly around the first intron (). While a tissue-independent model likely obscures some of the nuanced features regulating splicing-dependent gene expression, a genome-wide comparative analysis by revealed that correlation between CpG methylation and gene expression is unique to the first exon and intron. As CpG markers of DNA methylation tend to decrease across exons and increase across introns, it is possible that methylation may inform gene expression by mediating intron splice site recognition (Laurent et al., 2010). In fact, removal of promoter-proximal introns altogether reduces levels of H3K4me3 and chromatin-bound RNA polymerase II, reducing transcriptional output (; Laxa et al., 2016). Similarly, reduction in chromatin accessibility was observed when formation of the active spliceosome was inhibited by spliceostatin A. This finding highlights an important role for the spliceosome in regulating transcriptional output. Notably, this effect was not intrinsic to the presence of introns, but dependent on their splicing ().
One caveat to the spliceostatin A experiment is that it inhibits the entire splicing machinery, without revealing the specific interactions between the spliceosome and intron consensus sequences that enhance transcription. In fact, it is not the entire spliceosome that needs to be activated for transcription enhancement, as the formation of stable interactions between U1 snRNA and the promoter-proximal 5′SS can enhance transcription (Engreitz et al., 2014). Recruitment of the U1 snRNP to the first intron enhances transcription initiation through recruitment of general transcription factors, such as TFIIH, and stabilization the first formed phosphodiester bond by RNA polymerase II (Kwek et al., 2002; ). Notably, this effect is independent of its role in major intron splicing, as mere introduction of a 5′SS sequence is sufficient to enhance transcription (). This splicing-independent function of U1 might help explain its constitutive association with the elongating RNA polymerase II and why it is likewise recruited to intronless genes (Spiluttini et al., 2010; Leader et al., 2021).
Beside the role of U1 in transcription initiation, U1 snRNA is also independently involved in preventing pre-mature transcription termination, which can occur if RNA polymerase II encounters a polyadenylation site within an intron. Surmounting 3′ end sequencing data has revealed that introns often contain cryptic or pre-mature polyadenylation sites that result in the destabilization of RNA polymerase II, thereby producing truncated transcripts incapable of encoding a protein (Di Giammartino et al., 2011). Remarkably, the production of these truncated transcripts can be blocked by the U1 snRNA in a process called telescripting. In this capacity, U1 is capable of complexing with 3′ processing factors to protect the mRNA from premature cleavage and termination (Kaida et al., 2010; ). This mechanism occurs alongside the elongating polymerase to allow for U1-mediated suppression of cryptic polyadenylation sites in the intron or 3′ UTR (). Proper transcription termination is important in regulating the length and structure of the 3’ UTR, which in turn promotes formation of the export-competent messenger ribonucleoprotein. Similar to U1, U11 is expressed more highly than is necessary for its function in splicing (). Given that U11 is more abundant than U12 though they present at the same stoichiometric ratio within the minor spliceosome, U11 may similarly have splicing-independent functions. We speculate that U11 may either function in a mechanism like telescripting or participate in an alternative function, such as the subnuclear clustering of expressed minor intron-containing genes.
Localization of spliceosome components
Genes, chromatin, and RNA polymerase II have a subnuclear organization around topologically-associated domains to phase-separate euchromatic regions of active transcription (Szentirmay and Sawadogo, 2000; Ulianov et al., 2016; Szabo et al., 2020). Alongside this, it would be reasonable to hypothesize that splicing machinery is also organized to support efficient gene expression. In fact, major and minor spliceosome snRNPs display similar partiality for nuclear localization, except for U6 and U6atac snRNPs (Spiller et al., 2007; Pessa et al., 2008; Steitz et al., 2008). In the nucleus, matured snRNPs of the major spliceosome accumulate in phase-separated speckles that serve to organize spliceosome components adjacent to perichromatin regions of active transcription. This was concluded following nonradioactive and fluorescence in situ hybridization analyses, as well as RNA and protein blotting of subcellular compartment extracts (Pessa et al., 2008). While this model is an enticing way to interpret speckles as a regulatory mechanism over major intron splicing, it does not necessarily extend to that of minor introns. Given that the major and minor spliceosomes are known to interact with each other in the splicing of minor intron-containing genes, the model does not encompass all mechanisms of splicing (; Olthof et al., 2021). Punctate subcellular localization of spliceosome machinery is not specific to core snRNP components, but also includes some of the auxiliary splicing factors that contribute to spliceosome stability, conformational changes, and catalytic activity during splicing. These non-snRNP factors are integral to spliceosome assembly and the coordinated action of snRNPs during splicing (). For example, a new model supposes that the unequal phasic separation of SR proteins and heterogenous nuclear ribonucleoproteins proteins (hnRNP) at nuclear speckles can contribute to splice site selection. Specifically, the positional distribution of SR proteins and hnRNPs around a splice site generally determines the positive or negative regulatory effect of their binding, and taken with their distinct subnuclear distributions, can dictate the use of splice sites (Liao and Regev, 2021).
In all, here through an intron-centric lens, we focus our attention on the myriad of regulatory and functional consequences that have emerged by the presence of introns in the genome. Thus, we hope that future studies will begin to shed light on this “dark matter” of the eukaryotic genome to uncover the secrets buried within. Importantly, the advent of next-generation sequencing and computational analysis will invariably play a critical role in uncovering some of these mysteries. Throughout this article, we have described several of these methods, and here we point readers to other reviews (Halperin et al., 2021; Lorenzi et al., 2021; Gondane and Itkonen, 2023).
Statements
Author contributions
KG was responsible for curation of literature, organizing, writing text, and generating figures. AO was responsible for editing, structural organization, and help with literature curation. RK was responsible for the vision, writing, figures, and structure of the document. All authors contributed to the article and approved the submitted version.
Funding
Funding for this study comes from the National Institute of Neurological Disorders and Stroke (R01NS102538 to RK).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AbebreseE. L.AliS. H.ArnoldZ. R.AndrewsV. M.ArmstrongK.BurnsL.et al (2017). Identification of human short introns. PLoS ONE12, e0175393. 10.1371/journal.pone.0175393
2
AbelS.KissT.SolymosyF. (1989). Molecular analysis of eight U1 RNA gene candidates from tomato that could potentially be transcribed into U1 RNA sequence variants differing from each other in similar regions of secondary structure. Nucleic Acids Res.17, 6319–6337. 10.1093/nar/17.15.6319
3
Abou AlezzM.CelliL.BelottiG.LisaA.BioneS. (2020). GC-AG introns features in long non-coding and protein-coding genes suggest their role in gene expression regulation. Front. Genet. 11. 10.3389/fgene.2020.00488
4
AgirreE.OldfieldA. J.BelloraN.SegelleA.LucoR. F. (2021). Splicing-associated chromatin signatures: A combinatorial and position-dependent role for histone marks in splicing definition. Nat. Commun.12, 682. 10.1038/s41467-021-20979-x
5
AkinyiM. V.FrilanderM. J. (2021). At the intersection of major and minor spliceosomes: Crosstalk mechanisms and their impact on gene expression. Front. Genet.12, 700744. 10.3389/fgene.2021.700744
6
AliotoT. S. (2007). U12DB: A database of orthologous U12-type spliceosomal introns. Nucleic Acids Res.35, D110–D115. 10.1093/nar/gkl796
7
AlvarezM. E. V.ChiversM.BorovskaI.MongerS.GiannoulatouE.KralovicovaJ.et al (2021). Transposon clusters as substrates for aberrant splice-site activation. RNA Biol.18, 354–367. 10.1080/15476286.2020.1805909
8
AmitM.DonyoM.HollanderD.GorenA.KimE.GelfmanS.et al (2012). Differential GC content between exons and introns establishes distinct strategies of splice-site recognition. Cell Rep.1, 543–556. 10.1016/j.celrep.2012.03.013
9
AnastasiadiD.Esteve-CodinaA.PiferrerF. (2018). Consistent inverse correlation between DNA methylation of the first intron and gene expression across tissues and species. Epigenetics Chromatin11, 37. 10.1186/s13072-018-0205-1
10
AnderssonR.EnrothS.Rada-IglesiasA.WadeliusC.KomorowskiJ. (2009). Nucleosomes are well positioned in exons and carry characteristic histone modifications. Genome Res.19, 1732–1741. 10.1101/gr.092353.109
11
AnjosA.Ruiz-RuanoF. J.CamachoJ. P. M.LoretoV.CabreroJ.de SouzaM. J.et al (2015). U1 snDNA clusters in grasshoppers: Chromosomal dynamics and genomic organization. Heredity114, 207–219. 10.1038/hdy.2014.87
12
AnnaA.MonikaG. (2018). Splicing mutations in human genetic disorders: Examples, detection, and confirmation. J. Appl. Genet.59, 253–268. 10.1007/s13353-018-0444-7
13
BaiR.WanR.WangL.XuK.ZhangQ.LeiJ.et al (2021). Structure of the activated human minor spliceosome. Science371, eabg0879. 10.1126/science.abg0879
14
BalachandranP.WalawalkarI. A.FloresJ. I.DaytonJ. N.AudanoP. A.BeckC. R. (2022). Transposable element-mediated rearrangements are prevalent in human genomes. Nat. Commun.13, 7115. 10.1038/s41467-022-34810-8
15
BaralleF. E.GiudiceJ. (2017). Alternative splicing as a regulator of development and tissue identity. Nat. Rev. Mol. Cell Biol.18, 437–451. 10.1038/nrm.2017.27
16
BaumgartnerM.DrakeK.KanadiaR. N. (2019). An integrated model of minor intron emergence and conservation. Front. Genet.10, 1113. 10.3389/fgene.2019.01113
17
BaumgartnerM.LemoineC.Al SeesiS.KarunakaranD. K. P.SturrockN.BandayA. R.et al (2015). Minor splicing snRNAs are enriched in the developing mouse CNS and are crucial for survival of differentiating retinal neurons. Dev. Neurobiol.75, 895–907. 10.1002/dneu.22257
18
BergM. G.SinghL. N.YounisI.LiuQ.PintoA. M.KaidaD.et al (2012). U1 snRNP determines mRNA length and regulates isoform expression. Cell150, 53–64. 10.1016/j.cell.2012.05.029
19
BergetS. M. (1995). Exon recognition in vertebrate splicing. J. Biol. Chem.270, 2411–2414. 10.1074/jbc.270.6.2411
20
BiebersteinN. I.Carrillo OesterreichF.StraubeK.NeugebauerK. M. (2012). First exon length controls active chromatin signatures and transcription. Cell Rep.2, 62–68. 10.1016/j.celrep.2012.05.019
21
BindereifA.GreenM. R. (1990). Identification and functional analysis of mammalian splicing factors. Genet. Eng. (N. Y.)12, 201–224. 10.1007/978-1-4613-0641-2_11
22
BlackD. L. (2003). Mechanisms of alternative pre-messenger RNA splicing. Annu. Rev. Biochem.72, 291–336. 10.1146/annurev.biochem.72.121801.161720
23
BlankvoortS.WitterM. P.NoonanJ.CotneyJ.KentrosC. (2018). Marked diversity of unique cortical enhancers enables neuron-specific tools by enhancer-driven gene expression. Curr. Biol. CB28, 2103–2114. 10.1016/j.cub.2018.05.015
24
BreathnachR.BenoistC.O’HareK.GannonF.ChambonP. (1978). Ovalbumin gene: Evidence for a leader sequence in mRNA and DNA sequences at the exon-intron boundaries. Proc. Natl. Acad. Sci. U. S. A.75, 4853–4857. 10.1073/pnas.75.10.4853
25
BreathnachR.ChambonP. (1981). Organization and expression of eucaryotic split genes coding for proteins. Annu. Rev. Biochem.50, 349–383. 10.1146/annurev.bi.50.070181.002025
26
BringmannP.AppelB.RinkeJ.ReuterR.TheissenH.LührmannR. (1984). Evidence for the existence of snRNAs U4 and U6 in a single ribonucleoprotein complex and for their association by intermolecular base pairing. EMBO J.3, 1357–1363. 10.1002/j.1460-2075.1984.tb01977.x
27
BringmannP.LührmannR. (1986). Purification of the individual snRNPs U1, U2, U5 and U4/U6 from HeLa cells and characterization of their protein constituents. EMBO J.5, 3509–3516. 10.1002/j.1460-2075.1986.tb04676.x
28
BringmannP.RinkeJ.AppelB.ReuterR.LührmannR. (1983). Purification of snRNPs U1, U2, U4, U5 and U6 with 2,2,7-trimethylguanosine-specific antibody and definition of their constituent proteins reacting with anti-Sm and anti-(U1)RNP antisera. EMBO J.2, 1129–1135. 10.1002/j.1460-2075.1983.tb01557.x
29
BrinsterR. L.AllenJ. M.BehringerR. R.GelinasR. E.PalmiterR. D. (1988). Introns increase transcriptional efficiency in transgenic mice. Proc. Natl. Acad. Sci. U. S. A.85, 836–840. 10.1073/pnas.85.3.836
30
BurgeC. B.PadgettR. A.SharpP. A. (1998). Evolutionary fates and origins of U12-type introns. Mol. Cell2, 773–785. 10.1016/s1097-2765(00)80292-0
31
CáceresJ. F.McKenzieD.ThimmapayaR.LundE.DahlbergJ. E. (1992). Control of mouse U1a and U1b snRNA gene expression by differential transcription. Nucleic Acids Res.20, 4247–4254. 10.1093/nar/20.16.4247
32
CarmelL.WolfY. I.RogozinI. B.KooninE. V. (2007). Three distinct modes of intron dynamics in the evolution of eukaryotes. Genome Res.17, 1034–1044. 10.1101/gr.6438607
33
CarranzaF.ShenasaH.HertelK. J. (2022). Splice site proximity influences alternative exon definition. RNA Biol.19, 829–840. 10.1080/15476286.2022.2089478
34
CataniaF.GaoX.ScofieldD. G. (2009). Endogenous mechanisms for the origins of spliceosomal introns. J. Hered.100, 591–596. 10.1093/jhered/esp062
35
ChenL.LuilloD. J.MaE.CelnikerS. E.RioD. C.DoudnaJ. A. (2005). Identification and analysis of U5 snRNA variants in Drosophila. RNA11, 1473–1477. 10.1261/rna.2141505
36
CiavarellaJ.PereaW.GreenbaumN. L. (2020). Topology of the U12–U6 atac snRNA complex of the minor spliceosome and binding by NTC-related protein RBM22. ACS Omega5, 23549–23558. 10.1021/acsomega.0c01674
37
CrickF. (1979). Split genes and RNA splicing. Science204, 264–271. 10.1126/science.373120
38
CsurosM.RogozinI. B.KooninE. V. (2011). A detailed history of intron-rich eukaryotic ancestors inferred from a global survey of 100 complete genomes. PLOS Comput. Biol.7, e1002150. 10.1371/journal.pcbi.1002150
39
DamgaardC. K.KahnsS.Lykke-AndersenS.NielsenA. L.JensenT. H.KjemsJ. (2008). A 5’ splice site enhances the recruitment of basal transcription initiation factors in vivo. Mol. Cell29, 271–278. 10.1016/j.molcel.2007.11.035
40
De ContiL.BaralleM.BurattiE. (2013). Exon and intron definition in pre-mRNA splicing. WIREs RNA4, 49–60. 10.1002/wrna.1140
41
de Melo CostaV. R.PfeufferJ.LouloupiA.ØromU. A. V.PiroR. M. (2021). SPLICE-Q: A Python tool for genome-wide quantification of splicing efficiency. BMC Bioinforma.22, 368. 10.1186/s12859-021-04282-6
42
de WolfB.OghabianA.AkinyiM. V.HanksS.TromerE. C.van HooffJ. J. E.et al (2021). Chromosomal instability by mutations in the novel minor spliceosome component CENATAC. EMBO J.40, e106536. 10.15252/embj.2020106536
43
DenisonR. A.Van ArsdellS. W.BernsteinL. B.WeinerA. M. (1981). Abundant pseudogenes for small nuclear RNAs are dispersed in the human genome. Proc. Natl. Acad. Sci.78, 810–814. 10.1073/pnas.78.2.810
44
DergaiO.CousinP.GougeJ.SatiaK.PrazV.KuhlmanT.et al (2018). Mechanism of selective recruitment of RNA polymerases II and III to snRNA gene promoters. Genes Dev.32, 711–722. 10.1101/gad.314245.118
45
DeyP.MattickJ. S. (2021). High frequency of intron retention and clustered H3K4me3-marked nucleosomes in short first introns of human long non-coding RNAs. Epigenetics Chromatin14, 45. 10.1186/s13072-021-00419-2
46
DiC.SoB. R.CaiZ.AraiC.DuanJ.DreyfussG. (2019). U1 snRNP telescripting roles in transcription and its mechanism. Cold Spring Harb. Symp. Quant. Biol.84, 115–122. 10.1101/sqb.2019.84.040451
47
Di GiammartinoD. C.NishidaK.ManleyJ. L. (2011). Mechanisms and consequences of alternative polyadenylation. Mol. Cell43, 853–866. 10.1016/j.molcel.2011.08.017
48
DietrichR. C.IncorvaiaR.PadgettR. A. (1997). Terminal intron dinucleotide sequences do not distinguish between U2- and U12-dependent introns. Mol. Cell1, 151–160. 10.1016/s1097-2765(00)80016-7
49
DominovJ. A.UyanÖ.McKenna‐YasekD.NallamilliB. R. R.KergourlayV.BartoliM.et al (2019). Correction of pseudoexon splicing caused by a novel intronic dysferlin mutation. Ann. Clin. Transl. Neurol.6, 642–654. 10.1002/acn3.738
50
DomitrovichA. M.KunkelG. R. (2003). Multiple, dispersed human U6 small nuclear RNA genes with varied transcriptional efficiencies. Nucleic Acids Res.31, 2344–2352. 10.1093/nar/gkg331
51
DoucetA. J.DrocG.SiolO.AudouxJ.GilbertN. (2015). U6 snRNA pseudogenes: Markers of retrotransposition dynamics in mammals. Mol. Biol. Evol.32, 1815–1832. 10.1093/molbev/msv062
52
DrexlerH. L.ChoquetK.ChurchmanL. S. (2020). Splicing kinetics and coordination revealed by direct nascent RNA sequencing through nanopores. Mol. Cell77, 985–998. 10.1016/j.molcel.2019.11.017
53
EderyP.MarcaillouC.SahbatouM.LabalmeA.ChastangJ.TouraineR.et al (2011). Association of TALS developmental disorder with defect in minor splicing component U4atac snRNA. Science332, 240–243. 10.1126/science.1202205
54
EffenbergerK. A.UrabeV. K.JuricaM. S. (2017). Modulating splicing with small molecular inhibitors of the spliceosome: Modulating splicing with small molecular inhibitors. Wiley Interdiscip. Rev. RNA8, e1381. 10.1002/wrna.1381
55
ElliottB.RichardsonC.JasinM. (2005). Chromosomal translocation mechanisms at intronic Alu elements in mammalian cells. Mol. Cell17, 885–894. 10.1016/j.molcel.2005.02.028
56
ElsaidM. F.ChalhoubN.Ben-OmranT.KumarP.KamelH.IbrahimK.et al (2017). Mutation in noncoding RNA RNU12 causes early onset cerebellar ataxia. Ann. Neurol.81, 68–78. 10.1002/ana.24826
57
EmeraD.YinJ.ReillyS. K.GockleyJ.NoonanJ. P. (2016). Origin and evolution of developmental enhancers in the mammalian neocortex. Proc. Natl. Acad. Sci.113, E2617–E2626. 10.1073/pnas.1603718113
58
EngreitzJ. M.SirokmanK.McDonelP.ShishkinA. A.SurkaC.RussellP.et al (2014). RNA-RNA interactions enable specific targeting of noncoding RNAs to nascent Pre-mRNAs and chromatin sites. Cell159, 188–199. 10.1016/j.cell.2014.08.018
59
Fox-WalshK. L.DouY.LamB. J.HungS.BaldiP. F.HertelK. J. (2005). The architecture of pre-mRNAs affects mechanisms of splice-site pairing. Proc. Natl. Acad. Sci.102, 16176–16181. 10.1073/pnas.0508489102
60
FranchiniL. F.López-LealR.NasifS.BeatiP.GelmanD. M.LowM. J.et al (2011). Convergent evolution of two mammalian neuronal enhancers by sequential exaptation of unrelated retroposons. Proc. Natl. Acad. Sci. USA108, 15270–15275. 10.1073/pnas.1104997108
61
FriendK.LovejoyA. F.SteitzJ. A. (2007). U2 snRNP binds intronless histone pre-mRNAs to facilitate U7-snRNP-Dependent 3′-end formation. Mol. Cell28, 240–252. 10.1016/j.molcel.2007.09.026
62
FrigolaJ.SabarinathanR.MularoniL.MuiñosF.Gonzalez-PerezA.López-BigasN. (2017). Reduced mutation rate in exons due to differential mismatch repair. Nat. Genet.49, 1684–1692. 10.1038/ng.3991
63
FrumkinI.YofeI.Bar-ZivR.GurvichY.LuY.-Y.VoichekY.et al (2019). Evolution of intron splicing towards optimized gene expression is based on various Cis- and Trans-molecular mechanisms. PLoS Biol.17, e3000423. 10.1371/journal.pbio.3000423
64
FurgerA.O’SullivanJ. M.BinnieA.LeeB. A.ProudfootN. J. (2002). Promoter proximal splice sites enhance transcription. Genes Dev.16, 2792–2799. 10.1101/gad.983602
65
GehringN. H.RoignantJ.-Y. (2021). Anything but ordinary – emerging splicing mechanisms in eukaryotic gene regulation. Trends Genet.37, 355–372. 10.1016/j.tig.2020.10.008
66
GelfmanS.BursteinD.PennO.SavchenkoA.AmitM.SchwartzS.et al (2011). Changes in exon–intron structure during vertebrate evolution affect the splicing pattern of exons. Genome Res.22, 35–50. 10.1101/gr.119834.110
67
GeorgomanolisT.SofiadisK.PapantonisA. (2016). Cutting a long intron short: Recursive splicing and its implications. Front. Physiol.7, 598. 10.3389/fphys.2016.00598
68
GilbertW. (1978). Why genes in pieces?Nature271, 501. 10.1038/271501a0
69
GondaneA.ItkonenH. M. (2023). Revealing the history and mystery of RNA-seq. Curr. Issues Mol. Biol.45, 1860–1874. 10.3390/cimb45030120
70
GozashtiL.RoyS. W.ThornlowB.KramerA.AresM.Corbett-DetigR. (2022). Transposable elements drive intron gain in diverse eukaryotes. Proc. Natl. Acad. Sci.119, e2209766119. 10.1073/pnas.2209766119
71
GrabowskiP. J.SeilerS. R.SharpP. A. (1985). A multicomponent complex is involved in the splicing of messenger RNA precursors. Cell42, 345–353. 10.1016/S0092-8674(85)80130-6
72
GriffinC.Saint-JeannetJ.-P. (2020). Spliceosomopathies: Diseases and mechanisms. Dev. Dyn.249, 1038–1046. 10.1002/dvdy.214
73
HallS. L.PadgettR. A. (1994). Conserved sequences in a class of rare eukaryotic nuclear introns with non-consensus splice sites. J. Mol. Biol.239, 357–365. 10.1006/jmbi.1994.1377
74
HallS. L.PadgettR. A. (1996). Requirement of U12 snRNA for in vivo splicing of a minor class of eukaryotic nuclear pre-mRNA introns. Science271, 1716–1718. 10.1126/science.271.5256.1716
75
HalperinR. F.HegdeA.LangJ. D.RaupachE. A.LegendreC.LiangW. S.et al (2021). Improved methods for RNAseq-based alternative splicing analysis. Sci. Rep.11, 10740. 10.1038/s41598-021-89938-2
76
HenryR. W.MittalV.MaB.KobayashiR.HernandezN. (1998). SNAP19 mediates the assembly of a functional core promoter complex (SNAPc) shared by RNA polymerases II and III. Genes Dev.12, 2664–2672. 10.1101/gad.12.17.2664
77
HerzelL.OttozD. S. M.AlpertT.NeugebauerK. M. (2017). Splicing and transcription touch base: Co-transcriptional spliceosome assembly and function. Nat. Rev. Mol. Cell Biol.18, 637–650. 10.1038/nrm.2017.63
78
HesselberthJ. R. (2013). Lives that introns lead after splicing. Wiley Interdiscip. Rev. RNA4, 677–691. 10.1002/wrna.1187
79
HoffmanM. M.BirneyE. (2006). Estimating the neutral rate of nucleotide substitution using introns. Mol. Biol. Evol.24, 522–531. 10.1093/molbev/msl179
80
HuangY.GuL.LiG.-M. (2018). H3K36me3-mediated mismatch repair preferentially protects actively transcribed genes from mutation. J. Biol. Chem.293, 7811–7823. 10.1074/jbc.RA118.002839
81
HudsonA. J.McWattersD. C.BowserB. A.MooreA. N.LarueG. E.RoyS. W.et al (2019). Patterns of conservation of spliceosomal intron structures and spliceosome divergence in representatives of the diplomonad and parabasalid lineages. BMC Evol. Biol.19, 162. 10.1186/s12862-019-1488-y
82
HuffJ. T.ZilbermanD.RoyS. W. (2016). Mechanism for DNA transposons to generate introns on genomic scales. Nature538, 533–536. 10.1038/nature20110
83
IncorvaiaR.PadgettR. A. (1998). Base pairing with U6atac snRNA is required for 5’ splice site activation of U12-dependent introns in vivo. RNA4, 709–718. 10.1017/s1355838298980207
84
JacksonI. J. (1991). A reappraisal of non-consensus mRNA splice sites. Nucleic Acids Res.19, 3795–3798. 10.1093/nar/19.14.3795
85
JaktL. M.DubinA.JohansenS. D. (2022). Intron size minimisation in teleosts. BMC Genomics23, 628. 10.1186/s12864-022-08760-w
86
JawdekarG.HenryR. (2008). Transcriptional regulation of human small nuclear RNA genes. Biochim. Biophys. Acta BBA - Gene Regul. Mech.1779, 295–305. 10.1016/j.bbagrm.2008.04.001
87
JiangM.ZhangS.YinH.ZhuoZ.MengG. (2023). A comprehensive benchmarking of differential splicing tools for RNA-seq analysis at the event level. Brief. Bioinform., bbad121. bbad121. 10.1093/bib/bbad121
88
JosephB.ScalaC.KondoS.LaiE. C. (2022). Molecular and genetic dissection of recursive splicing. Life Sci. Alliance5, e202101063. 10.26508/lsa.202101063
89
KaidaD.BergM. G.YounisI.KasimM.SinghL. N.WanL.et al (2010). U1 snRNP protects pre-mRNAs from premature cleavage and polyadenylation. Nature468, 664–668. 10.1038/nature09479
90
KaidaD.MotoyoshiH.TashiroE.NojimaT.HagiwaraM.IshigamiK.et al (2007). Spliceostatin A targets SF3b and inhibits both splicing and nuclear retention of pre-mRNA. Nat. Chem. Biol.3, 576–583. 10.1038/nchembio.2007.18
91
KandulN. P.NoorM. A. (2009). Large introns in relation to alternative splicing and gene evolution: A case study of Drosophila bruno-3. BMC Genet.10, 67. 10.1186/1471-2156-10-67
92
KapustinY.ChanE.SarkarR.WongF.VorechovskyI.WinstonR. M.et al (2011). Cryptic splice sites and split genes. Nucleic Acids Res.39, 5837–5844. 10.1093/nar/gkr203
93
KazazianH. H.MoranJ. V. (1998). The impact of L1 retrotransposons on the human genome. Nat. Genet.19, 19–24. 10.1038/ng0598-19
94
KeaneP. A.SeoigheC. (2016). Intron length coevolution across mammalian genomes. Mol. Biol. Evol.33, 2682–2691. 10.1093/molbev/msw151
95
KolossovaI.PadgettR. A. (1997). U11 snRNA interacts in vivo with the 5’ splice site of U12-dependent (AU-AC) pre-mRNA introns. RNA3, 227–233.
96
KooninE. V. (2006). The origin of introns and their role in eukaryogenesis: A compromise solution to the introns-early versus introns-late debate?Biol. Direct1, 22. 10.1186/1745-6150-1-22
97
KrchňákováZ.ThakurP. K.KrausováM.BiebersteinN.HabermanN.Müller-McNicollM.et al (2019). Splicing of long non-coding RNAs primarily depends on polypyrimidine tract and 5′ splice-site sequences due to weak interactions with SR proteins. Nucleic Acids Res.47, 911–928. 10.1093/nar/gky1147
98
KumariA.SedehizadehS.BrookJ. D.KozlowskiP.WojciechowskaM. (2022). Differential fates of introns in gene expression due to global alternative splicing. Hum. Genet.141, 31–47. 10.1007/s00439-021-02409-6
99
KwekK. Y.MurphyS.FurgerA.ThomasB.O’GormanW.KimuraH.et al (2002). U1 snRNA associates with TFIIH and regulates transcriptional initiation. Nat. Struct. Biol.9, 800–805. 10.1038/nsb862
100
LambowitzA. M.ZimmerlyS. (2011). Group II introns: Mobile ribozymes that invade DNA. Cold Spring Harb. Perspect. Biol.3, a003616. 10.1101/cshperspect.a003616
101
LanderE. S.LintonL. M.BirrenB.NusbaumC.ZodyM. C.BaldwinJ.et al (2001). Initial sequencing and analysis of the human genome. Nature409, 860–921. 10.1038/35057062
102
LarueG. E.EliášM.RoyS. W. (2021). Expansion and transformation of the minor spliceosomal system in the slime mold Physarum polycephalum. Curr. Biol.31, 3125–3131.e4. 10.1016/j.cub.2021.04.050
103
LaurentL.WongE.LiG.HuynhT.TsirigosA.OngC. T.et al (2010). Dynamic changes in the human methylome during differentiation. Genome Res.20, 320–331. 10.1101/gr.101907.109
104
LaxaM.MüllerK.LangeN.DoeringL.PruschaJ. T.PeterhänselC. (2016). The 5’UTR intron of arabidopsis GGT1 aminotransferase enhances promoter activity by recruiting RNA polymerase II. Plant Physiol.172, 313–327. 10.1104/pp.16.00881
105
LeaderY.Lev MaorG.SorekM.ShayevitchR.HusseinM.HameiriO.et al (2021). The upstream 5′ splice site remains associated to the transcription machinery during intron synthesis. Nat. Commun.12, 4545. 10.1038/s41467-021-24774-6
106
LeeS.StevensS. W. (2016). Spliceosomal intronogenesis. Proc. Natl. Acad. Sci.113, 6514–6519. 10.1073/pnas.1605113113
107
Lev-MaorG.SorekR.LevanonE. Y.PazN.EisenbergE.AstG. (2007). RNA-editing-mediated exon evolution. Genome Biol.8, R29. 10.1186/gb-2007-8-2-r29
108
LiX.LiuS.ZhangL.IssaianA.HillR. C.EspinosaS.et al (2019). A unified mechanism for intron and exon definition and back-splicing. Nature573, 375–380. 10.1038/s41586-019-1523-6
109
LiY.XuY.MaZ. (2017). Comparative analysis of the exon-intron structure in eukaryotic genomes. Yangtze Med.01, 50–64. 10.4236/ym.2017.11006
110
LiaoS. E.RegevO. (2021). Splicing at the phase-separated nuclear speckle interface: A model. Nucleic Acids Res.49, 636–645. 10.1093/nar/gkaa1209
111
LinC.-F.MountS. M.JarmołowskiA.MakałowskiW. (2010). Evolutionary dynamics of U12-type spliceosomal introns. BMC Evol. Biol.10, 47. 10.1186/1471-2148-10-47
112
LiuJ.MaxwellE. S. (1990). Mouse U14 snRNA is encoded in an intron of the mouse cognate hsc70 heat shock gene. Nucleic Acids Res.18, 6565–6571. 10.1093/nar/18.22.6565
113
LoP. C.MountS. M. (1990). Drosophila melanogaster genes for U1 snRNA variants and their expression during development. Nucleic Acids Res.18, 6971–6979. 10.1093/nar/18.23.6971
114
LorenziC.BarriereS.ArnoldK.LucoR. F.OldfieldA. J.RitchieW. (2021). IRFinder-S: A comprehensive suite to discover and explore intron retention. Genome Biol.22, 307. 10.1186/s13059-021-02515-8
115
LuZ.MateraA. G. (2014). Developmental analysis of spliceosomal snRNA isoform expression. G3 Bethesda Md5, 103–110. 10.1534/g3.114.015735
116
LynchM.RichardsonA. O. (2002). The evolution of spliceosomal introns. Curr. Opin. Genet. Dev.12, 701–710. 10.1016/s0959-437x(02)00360-x
117
MabinJ. W.LewisP. W.BrowD. A.DvingeH. (2021). Human spliceosomal snRNA sequence variants generate variant spliceosomes. RNA27, 1186–1203. 10.1261/rna.078768.121
118
MalcaH.ShomronN.AstG. (2003). The U1 snRNP base pairs with the 5′ splice site within a penta-snRNP complex. Am. Soc. Microbiol.3, 3442–3455. 10.1128/MCB.23.10.3442–3455.2003
119
McClintockB. (1950). The origin and behavior of mutable loci in maize. Proc. Natl. Acad. Sci. U. S. A.36, 344–355. 10.1073/pnas.36.6.344
120
MengF.ZhaoH.ZhuB.ZhangT.YangM.LiY.et al (2021). Genomic editing of intronic enhancers unveils their role in fine-tuning tissue-specific gene expression in Arabidopsis thaliana. Plant Cell33, 1997–2014. 10.1093/plcell/koab093
121
MercerT. R.ClarkM. B.AndersenS. B.BrunckM. E.HaertyW.CrawfordJ.et al (2015). Genome-wide discovery of human splicing branchpoints. Genome Res.25, 290–303. 10.1101/gr.182899.114
122
MichelF.UmesonoK.OzekiH. (1989). Comparative and functional anatomy of group II catalytic introns-a review. Gene82, 5–30. 10.1016/0378-1119(89)90026-7
123
MittalV.MaB.HernandezN. (1999). SNAP(c): A core promoter factor with a built-in DNA-binding damper that is deactivated by the oct-1 POU domain. Genes Dev.13, 1807–1821. 10.1101/gad.13.14.1807
124
MontzkaK. A.SteitzJ. A. (1988). Additional low-abundance human small nuclear ribonucleoproteins: U11, U12, etc. Proc. Natl. Acad. Sci. U A85, 8885–8889. 10.1073/pnas.85.23.8885
125
MoralesJ.BorreroM.SumerelJ.SantiagoC. (1997). Identification of developmentally regulated sea urchin U5 snRNA genes. DNA Seq. J. DNA Seq. Mapp.7, 243–259. 10.3109/10425179709034044
126
MoyerD. C.LarueG. E.HershbergerC. E.RoyS. W.PadgettR. A. (2020). Comprehensive database and evolutionary dynamics of U12-type introns. Nucleic Acids Res.48, 7066–7078. 10.1093/nar/gkaa464
127
NeugebauerK. M. (2019). Nascent RNA and the coordination of splicing with transcription. Cold Spring Harb. Perspect. Biol.11, a032227. 10.1101/cshperspect.a032227
128
NilsenT. W. (2003). The spliceosome: The most complex macromolecular machine in the cell?Bioessays25, 1147–1149. 10.1002/bies.10394
129
NojimaT.GomesT.GrossoA. R. F.KimuraH.DyeM. J.DhirS.et al (2015). Mammalian NET-seq reveals genome-wide nascent transcription coupled to RNA processing. Cell161, 526–540. 10.1016/j.cell.2015.03.027
130
NorppaA. J.FrilanderM. J. (2021). The integrity of the U12 snRNA 3′ stem–loop is necessary for its overall stability. Nucleic Acids Res.49, 2835–2847. 10.1093/nar/gkab048
131
OlthofA. M.HyattK. C.KanadiaR. N. (2019). Minor intron splicing revisited: Identification of new minor intron-containing genes and tissue-dependent retention and alternative splicing of minor introns. BMC Genomics20, 686. 10.1186/s12864-019-6046-x
132
OlthofA. M.WhiteA. K.KanadiaR. N. (2022). The emerging significance of splicing in vertebrate development. Development149, dev200373. 10.1242/dev.200373
133
OlthofA. M.WhiteA. K.MieruszynskiS.DoggettK.LeeM. F.ChakrounA.et al (2021). Disruption of exon-bridging interactions between the minor and major spliceosomes results in alternative splicing around minor introns. Nucleic Acids Res.49, 3524–3545. 10.1093/nar/gkab118
134
O’ReillyD.DienstbierM.CowleyS. A.VazquezP.DrozdzM.TaylorS.et al (2013). Differentially expressed, variant U1 snRNAs regulate gene expression in human cells. Genome Res.23, 281–291. 10.1101/gr.142968.112
135
PaiA. A.HenriquesT.McCueK.BurkholderA.AdelmanK.BurgeC. B. (2017). The kinetics of pre-mRNA splicing in the Drosophila genome and the influence of gene architecture. eLife6, e32537. 10.7554/eLife.32537
136
ParenteauJ.DurandM.VéronneauS.LacombeA.-A.MorinG.GuérinV.et al (2008). Deletion of many yeast introns reveals a minority of genes that require splicing for function. Mol. Biol. Cell19, 1932–1941. 10.1091/mbc.E07-12-1254
137
PessaH. K. J.WillC. L.MengX.SchneiderC.WatkinsN. J.PeräläN.et al (2008). Minor spliceosome components are predominantly localized in the nucleus. Proc. Natl. Acad. Sci.105, 8655–8660. 10.1073/pnas.0803646105
138
PiovesanA.CaracausiM.RicciM.StrippoliP.VitaleL.PelleriM. C. (2015). Identification of minimal eukaryotic introns through GeneBase, a user-friendly tool for parsing the NCBI Gene databank. DNA Res. Int. J. Rapid Publ. Rep. Genes Genomes22, 495–503. 10.1093/dnares/dsv028
139
PitolliC.MariniA.SetteC.PagliariniV. (2022). Non-canonical splicing and its implications in brain physiology and cancer. Int. J. Mol. Sci.23, 2811. 10.3390/ijms23052811
140
QianX.WangJ.WangM.IgelmanA. D.JonesK. D.LiY.et al (2021). Identification of deep-intronic splice mutations in a large cohort of patients with inherited retinal diseases. Front. Genet.12. 10.3389/fgene.2021.647400
141
ReddyR.HenningD.DasG.HarlessM.WrightD. (1987). The capped U6 small nuclear RNA is transcribed by RNA polymerase III. J. Biol. Chem.262, 75–81. 10.1016/s0021-9258(19)75890-6
142
ReedR.ManiatisT. (1986). A role for exon sequences and splice-site proximity in splice-site selection. Cell46, 681–690. 10.1016/0092-8674(86)90343-0
143
ReimerK. A.MimosoC. A.AdelmanK.NeugebauerK. M. (2021). Co-transcriptional splicing regulates 3′ end cleavage during mammalian erythropoiesis. Mol. Cell81, 998–1012.e7. 10.1016/j.molcel.2020.12.018
144
ReschA. M.CarmelL.Mariño-RamírezL.OgurtsovA. Y.ShabalinaS. A.RogozinI. B.et al (2007). Widespread positive selection in synonymous sites of mammalian genes. Mol. Biol. Evol.24, 1821–1831. 10.1093/molbev/msm100
145
RobbersonB. L.CoteG. J.BergetS. M. (1990). Exon definition may facilitate splice site selection in RNAs with multiple exons. Mol. Cell. Biol.10, 84–94. 10.1128/mcb.10.1.84
146
RocaX.SachidanandamR.KrainerA. R. (2003). Intrinsic differences between authentic and cryptic 5′ splice sites. Nucleic Acids Res.31, 6321–6333. 10.1093/nar/gkg830
147
Rodriguez-GalindoM.CasillasS.WeghornD.BarbadillaA. (2020). Germline de novo mutation rates on exons versus introns in humans. Nat. Commun.11, 3304. 10.1038/s41467-020-17162-z
148
RogozinI. B.CarmelL.CsurosM.KooninE. V. (2012). Origin and evolution of spliceosomal introns. Biol. Direct7, 11. 10.1186/1745-6150-7-11
149
RomfoC. M.AlvarezC. J.van HeeckerenW. J.WebbC. J.WiseJ. A. (2000). Evidence for splice site pairing via intron definition in Schizosaccharomyces pombe. Mol. Cell. Biol.20, 7955–7970. 10.1128/mcb.20.21.7955-7970.2000
150
RoseA. B.ElfersiT.ParraG.KorfI. (2008). Promoter-proximal introns in Arabidopsis thaliana are enriched in dispersed signals that elevate gene expression. Plant Cell20, 543–551. 10.1105/tpc.107.057190
151
RussellA. G.CharetteJ. M.SpencerD. F.GrayM. W. (2006). An early evolutionary origin for the minor spliceosome. Nature443, 863–866. 10.1038/nature05228
152
RyllJ.RotheringR.CataniaF. (2022). Intronization signatures in coding exons reveal the evolutionary fluidity of eukaryotic gene architecture. Microorganisms10, 1901. 10.3390/microorganisms10101901
153
SadowskiC. L.HenryR. W.LoboS. M.HernandezN. (1993). Targeting TBP to a non-TATA box cis-regulatory element: A TBP-containing complex activates transcription from snRNA promoters through the PSE. Genes Dev.7, 1535–1548. 10.1101/gad.7.8.1535
154
SakharkarM. K.ChowV. T. K.KangueaneP. (2004). Distributions of exons and introns in the human genome. Silico Biol.4, 387–393.
155
SamadderP.SivamaniE.LuJ.LiX.QuR. (2008). Transcriptional and post-transcriptional enhancement of gene expression by the 5’ UTR intron of rice rubi3 gene in transgenic rice cells. Mol. Genet. Genomics MGG279, 429–439. 10.1007/s00438-008-0323-8
156
Sánchez-EscabiasE.Guerrero-MartínezJ. A.ReyesJ. C. (2022). Co-transcriptional splicing efficiency is a gene-specific feature that can be regulated by TGFβ. Nat. Commun.5, 277. 10.1038/s42003-022-03224-z
157
SanMiguelP.TikhonovA.JinY.-K.MotchoulskaiaN.ZakharovD.Melake-BerhanA.et al (1996). Nested retrotransposons in the intergenic regions of the maize genome. Science274, 765–768. 10.1126/science.274.5288.765
158
SealR. L.ChenL.-L.Griffiths-JonesS.LoweT. M.MathewsM. B.O’ReillyD.et al (2020). A guide to naming human non-coding RNA genes. EMBO J.39, e103777. 10.15252/embj.2019103777
159
SharpP. A. (1991). Five easy pieces. Science254, 663. 10.1126/science.1948046
160
ShenH.ZhengX.LueckeS.GreenM. R. (2010). The U2AF35-related protein Urp contacts the 3′ splice site to promote U12-type intron splicing and the second step of U2-type intron splicing. Genes Dev.24, 2389–2394. 10.1101/gad.1974810
161
ShepardS.McCrearyM.FedorovA. (2009). The peculiarities of large intron splicing in animals. PLOS ONE4, e7853. 10.1371/journal.pone.0007853
162
SheridanR. M.FongN.D’AlessandroA.BentleyD. L. (2019). Widespread backtracking by RNA pol II is a major effector of gene activation, 5’ pause release, termination, and transcription elongation rate. Mol. Cell73, 107–118. 10.1016/j.molcel.2018.10.031
163
ShethN.RocaX.HastingsM. L.RoederT.KrainerA. R.SachidanandamR. (2006). Comprehensive splice-site analysis using comparative genomics. Nucleic Acids Res.34, 3955–3967. 10.1093/nar/gkl556
164
ShiauC.-K.HuangJ.-H.LiuY.-T.TsaiH.-K. (2022). Genome-wide identification of associations between enhancer and alternative splicing in human and mouse. BMC Genomics22, 919. 10.1186/s12864-022-08537-1
165
ShuklaG. C.PadgettR. A. (2002). A catalytically active group II intron domain 5 can function in the U12-dependent spliceosome. Mol. Cell9, 1145–1150. 10.1016/S1097-2765(02)00505-1
166
SiebertA. E.CorllJ.Paige GroneveltJ.LevineL.HobbsL. M.KenneyC.et al (2022). Genetic analysis of human RNA binding motif protein 48 (RBM48) reveals an essential role in U12-type intron splicing. Genetics222, iyac129. 10.1093/genetics/iyac129
167
Sierra-MontesJ. M.Pereira-SimonS.SmailS. S.HerreraR. J. (2005). The silk moth Bombyx mori U1 and U2 snRNA variants are differentially expressed. Gene352, 127–136. 10.1016/j.gene.2005.02.013
168
SimariR. D.YangZ.-Y.LingX.StephanD.PerkinsN. D.NabelG. J.et al (1998). Requirements for enhanced transgene expression by untranslated sequences from the human cytomegalovirus immediate-early gene. Mol. Med.4, 700–706. 10.1007/BF03401764
169
SinghJ.PadgettR. A. (2009). Rates of in situ transcription and splicing in large human genes. Nat. Struct. Mol. Biol.16, 1128–1133. 10.1038/nsmb.1666
170
SinghP.BörgerC.MoreH.SturmbauerC. (2017). The role of alternative splicing and differential gene expression in Cichlid adaptive radiation. Genome Biol. Evol.9, 2764–2781. 10.1093/gbe/evx204
171
SmathersC. M.RobartA. R. (2019). The mechanism of splicing as told by group II introns: Ancestors of the spliceosome. Biochim. Biophys. Acta Gene Regul. Mech.1862, 194390. 10.1016/j.bbagrm.2019.06.001
172
SoB. R.DiC.CaiZ.VentersC. C.GuoJ.OhJ.-M.et al (2019). A complex of U1 snRNP with cleavage and polyadenylation factors controls telescripting, regulating mRNA transcription in human cells. Mol. Cell76, 590–599. 10.1016/j.molcel.2019.08.007
173
SontheimerE. J.GordonP. M.PiccirilliJ. A. (1999). Metal ion catalysis during group II intron self-splicing: Parallels with the spliceosome. Genes Dev.13, 1729–1741. 10.1101/gad.13.13.1729
174
SontheimerE. J.SteitzJ. A. (1992). Three novel functional variants of human U5 small nuclear RNA. Mol. Cell Biol.12, 734–746. 10.1128/mcb.12.2.734
175
SorekR.AstG.GraurD. (2002). Alu-containing exons are alternatively spliced. Genome Res.12, 1060–1067. 10.1101/gr.229302
176
SpillerM. P.BoonK.-L.ReijnsM. A. M.BeggsJ. D. (2007). The Lsm2-8 complex determines nuclear localization of the spliceosomal U6 snRNA. Nucleic Acids Res.35, 923–929. 10.1093/nar/gkl1130
177
SpiluttiniB.GuB.BelagalP.SmirnovaA. S.NguyenV. T.HebertC.et al (2010). Splicing-independent recruitment of U1 snRNP to a transcription unit in living cells. J. Cell Sci.123, 2085–2093. 10.1242/jcs.061358
178
SteitzJ. A.DreyfussG.KrainerA. R.LamondA. I.MateraA. G.PadgettR. A. (2008). Where in the cell is the minor spliceosome?Proc. Natl. Acad. Sci. U. S. A.105, 8485–8486. 10.1073/pnas.0804024105
179
SunH.ChasinL. A. (2000). Multiple splicing defects in an intronic false exon. Mol. Cell. Biol.20, 6414–6425. 10.1128/MCB.20.17.6414-6425.2000
180
SunJ.LiX.HouX.CaoS.CaoW.ZhangY.et al (2022). Structural basis of human SNAPc recognizing proximal sequence element of snRNA promoter. Nat. Commun.13, 6871. 10.1038/s41467-022-34639-1
181
SzaboQ.DonjonA.JerkovićI.PapadopoulosG. L.CheutinT.BonevB.et al (2020). Regulation of single-cell genome organization into TADs and chromatin nanodomains. Nat. Genet.52, 1151–1157. 10.1038/s41588-020-00716-8
182
SzentirmayM. N.SawadogoM. (2000). Spatial organization of RNA polymerase II transcription in the nucleus. Nucleic Acids Res.28, 2019–2025. 10.1093/nar/28.10.2019
183
TammerL.HameiriO.KeydarI.RoyV. R.Ashkenazy-TitelmanA.CustódioN.et al (2022). Gene architecture directs splicing outcome in separate nuclear spatial regions. Mol. Cell82, 1021–1034.e8. 10.1016/j.molcel.2022.02.001
184
TarnW. Y.SteitzJ. A. (1996a). A novel spliceosome containing U11, U12, and U5 snRNPs excises a minor class (AT-AC) intron in vitro. Cell84, 801–811. 10.1016/s0092-8674(00)81057-0
185
TarnW. Y.SteitzJ. A. (1996b). Highly diverged U4 and U6 small nuclear RNAs required for splicing rare AT-AC introns. Science273, 1824–1832. 10.1126/science.273.5283.1824
186
TarnW. Y.SteitzJ. A. (1997). Pre-mRNA splicing: The discovery of a new spliceosome doubles the challenge. Trends Biochem. Sci.22, 132–137. 10.1016/s0968-0004(97)01018-9
187
TellierM.MaudlinI.MurphyS. (2020). Transcription and splicing: A two‐way street. WIREs RNA11, e1593. 10.1002/wrna.1593
188
TheissenH.RinkeJ.TraverC. N.LührmannR.AppelB. (1985). Novel structure of a human U6 snRNA pseudogene. Gene36, 195–199. 10.1016/0378-1119(85)90086-1
189
ThompsonP. J.MacfarlanT. S.LorinczM. C. (2016). Long terminal repeats: From parasitic elements to building blocks of the transcriptional regulatory repertoire. Mol. Cell62, 766–776. 10.1016/j.molcel.2016.03.029
190
TichelaarJ. W.WiebenE. D.ReddyR.VrabelA.CamachoP. (1998). In vivo expression of a variant human U6 RNA from a unique, internal promoter. Biochemistry37, 12943–12951. 10.1021/bi9811361
191
TronchèreH.WangJ.FuX. D. (1997). A protein related to splicing factor U2AF35 that interacts with U2AF65 and SR proteins in splicing of pre-mRNA. Nature388, 397–400. 10.1038/41137
192
UlianovS. V.KhrameevaE. E.GavrilovA. A.FlyamerI. M.KosP.MikhalevaE. A.et al (2016). Active chromatin and transcription play a key role in chromosome partitioning into topologically associating domains. Genome Res.26, 70–84. 10.1101/gr.196006.115
193
VakirlisN.VanceZ.DugganK. M.McLysaghtA. (2022). De novo birth of functional microproteins in the human lineage. Cell Rep.41, 111808. 10.1016/j.celrep.2022.111808
194
Van ArsdellS. W.WeinerA. M. (1984). Human genes for U2 small nuclear RNA are tandemly repeated. Mol. Cell. Biol.4, 492–499. 10.1128/mcb.4.3.492
195
Vazquez-ArangoP.VowlesJ.BrowneC.HartfieldE.FernandesH. J. R.MandefroB.et al (2016). Variant U1 snRNAs are implicated in human pluripotent stem cell maintenance and neuromuscular disease. Nucleic Acids Res.44, 10960–10973. 10.1093/nar/gkw711
196
VossebergJ.SnelB. (2017). Domestication of self-splicing introns during eukaryogenesis: The rise of the complex spliceosomal machinery. Biol. Direct12, 30. 10.1186/s13062-017-0201-6
197
WahlM. C.WillC. L.LührmannR. (2009). The spliceosome: Design principles of a dynamic RNP machine. Cell136, 701–718. 10.1016/j.cell.2009.02.009
198
WallaceE. W. J.BeggsJ. D. (2017). Extremely fast and incredibly close: Cotranscriptional splicing in budding yeast. RNA23, 601–610. 10.1261/rna.060830.117
199
WangZ.XiaoX.Van NostrandE.BurgeC. B. (2006). General and specific functions of exonic splicing silencers in splicing control. Mol. Cell23, 61–70. 10.1016/j.molcel.2006.05.018
200
WellsJ. N.FeschotteC. (2020). A field guide to eukaryotic transposable elements. Annu. Rev. Genet.54, 539–561. 10.1146/annurev-genet-040620-022145
201
WillC. L.SchneiderC.ReedR.LührmannR. (1999). Identification of both shared and distinct proteins in the major and minor spliceosomes. Science284, 2003–2005. 10.1126/science.284.5422.2003
202
YanD.AresM. (1996). Invariant U2 RNA sequences bordering the branchpoint recognition region are essential for interaction with yeast SF3a and SF3b subunits. Mol. Cell. Biol.16, 818–828. 10.1128/mcb.16.3.818
203
YounisI.DittmarK.WangW.FoleyS. W.BergM. G.HuK. Y.et al (2013). Minor introns are embedded molecular switches regulated by highly unstable U6atac snRNA. eLife2, e00780. 10.7554/eLife.00780
204
ZhangS.AibaraS.VosS. M.AgafonovD. E.LührmannR.CramerP. (2021). Structure of a transcribing RNA polymerase II-U1 snRNP complex. Science371, 305–309. 10.1126/science.abf1870
205
ZhangX.-O.FuY.MouH.XueW.WengZ. (2018). The temporal landscape of recursive splicing during Pol II transcription elongation in human cells. PLOS Genet.14, e1007579. 10.1371/journal.pgen.1007579
206
ZhuD.MaoF.TianY.LinX.GuL.GuH.et al (2020). The features and regulation of Co-transcriptional splicing in arabidopsis. Mol. Plant13, 278–294. 10.1016/j.molp.2019.11.004
Summary
Keywords
intron, evolution, splicing, snRNA, spliceosome, eukaryotes, gene expression
Citation
Girardini KN, Olthof AM and Kanadia RN (2023) Introns: the “dark matter” of the eukaryotic genome. Front. Genet. 14:1150212. doi: 10.3389/fgene.2023.1150212
Received
23 January 2023
Accepted
28 April 2023
Published
16 May 2023
Volume
14 - 2023
Edited by
Yadong Zheng, Zhejiang Agriculture and Forestry University, China
Reviewed by
Hari Krishna Yalamanchili, Baylor College of Medicine, United States
Bruce McKay, Carleton University, Canada
Updates
Copyright
© 2023 Girardini, Olthof and Kanadia.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rahul N. Kanadia, rahul.kanadia@uconn.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.