The intramolecular stem-loop structure of U6 snRNA can functionally replace the U6atac snRNA stem-loop

Clinical interpretation of variants identified in RNU4ATAC, a non-coding spliceosomal gene

Biallelic variants in RNU4ATAC, a non-coding gene transcribed into the minor spliceosome component U4atac snRNA, are responsible for three rare recessive developmental diseases, namely Taybi-Linder/MOPD1, Roifman and Lowry-Wood syndromes. Next-generation sequencing of clinically heterogeneous cohorts (children with either a suspected genetic disorder or a congenital microcephaly) recently identified mutations in this gene, illustrating how profoundly these technologies are modifying genetic testing and assessment. As RNU4ATAC has a single non-coding exon, the bioinformatic prediction algorithms assessing the effect of sequence variants on splicing or protein function are irrelevant, which makes variant interpretation challenging to molecular diagnostic laboratories. In order to facilitate and improve clinical diagnostic assessment and genetic counseling, we present i) an update of the previously reported RNU4ATAC mutations and an analysis of the genetic variations affecting this gene using the Genome Aggregation Database (gnomAD) resource; ii) the pathogenicity prediction performances of scores computed based on an RNA structure prediction tool and of those produced by the Combined Annotation Dependent Depletion tool for the 285 RNU4ATAC variants identified in patients or in large-scale sequencing projects; iii) a method, based on a cellular assay, that allows to measure the effect of RNU4ATAC variants on splicing efficiency of a minor (U12-type) reporter intron. Lastly, the concordance of bioinformatic predictions and cellular assay results was investigated.

Rearrangements within the U6 snRNA Core during the Transition between the Two Catalytic Steps of Splicing

The RNA catalytic core of spliceosomes as visualized by cryoelectron microscopy (cryo-EM) remains unchanged at different stages of splicing. However, we demonstrate that mutations within the core of yeast U6 snRNA modulate conformational changes between the two catalytic steps. We propose that the intramolecular stem-loop (ISL) of U6 exists in two competing states, changing between a default, non-catalytic conformation and a transient, catalytic conformation. Whereas stable interactions in the catalytic triplex promote catalysis and their disruptions favor exit from the catalytic conformation, destabilization of the lower ISL stem promotes catalysis and its stabilization supports exit from the catalytic conformation. Thus, in addition to the catalytic triplex, U6-ISL acts as an important dynamic component of the catalytic center. The relative flexibility of the lower U6-ISL stem is conserved across eukaryotes. Similar features are found in U6atac and domain V of group II introns, arguing for the generality of the proposed mechanism.

The life of U6 small nuclear RNA, from cradle to grave

Removal of introns from precursor messenger RNA (pre-mRNA) and some noncoding transcripts is an essential step in eukaryotic gene expression. In the nucleus, this process of RNA splicing is carried out by the spliceosome, a multi-megaDalton macromolecular machine whose core components are conserved from yeast to humans. In addition to many proteins, the spliceosome contains five uridine-rich small nuclear RNAs (snRNAs) that undergo an elaborate series of conformational changes to correctly recognize the splice sites and catalyze intron removal. Decades of biochemical and genetic data, along with recent cryo-EM structures, unequivocally demonstrate that U6 snRNA forms much of the catalytic core of the spliceosome and is highly dynamic, interacting with three snRNAs, the pre-mRNA substrate, and >25 protein partners throughout the splicing cycle. This review summarizes the current state of knowledge on how U6 snRNA is synthesized, modified, incorporated into snRNPs and spliceosomes, recycled, and degraded.

Splicing diversity revealed by reduced spliceosomes in C. merolae and other organisms

Pre-mRNA splicing has been considered one of the hallmarks of eukaryotes, yet its diversity is astonishing: the number of substrate introns for splicing ranges from hundreds of thousands in humans to a mere handful in certain parasites. The catalytic machinery that carries out splicing, the spliceosome, is similarly diverse, with over 300 associated proteins in humans to a few tens in other organisms. In this Point of View, we discuss recent work characterizing the reduced spliceosome of the acidophilic red alga Cyanidioschyzon merolae, which further highlights the diversity of splicing in that it does not possess the U1 snRNP that is characteristically responsible for 5′ splice site recognition. Comparisons to other organisms with reduced spliceosomes, such as microsporidia, trypanosomes, and Giardia, help to identify the most highly conserved splicing factors, pointing to the essential core of this complex machine. These observations argue for increased exploration of important biochemical processes through study of a wider ranger of organisms.

Biochemical defects in minor spliceosome function in the developmental disorder MOPD I

Biallelic mutations of the human RNU4ATAC gene, which codes for the minor spliceosomal U4atac snRNA, cause the developmental disorder, MOPD I/TALS. To date, nine separate mutations in RNU4ATAC have been identified in MOPD I patients. Evidence suggests that all of these mutations lead to abrogation of U4atac snRNA function and impaired minor intron splicing. However, the molecular basis of these effects is unknown. Here, we use a variety of in vitro and in vivo assays to address this question. We find that only one mutation, 124G>A, leads to significantly reduced expression of U4atac snRNA, whereas four mutations, 30G>A, 50G>A, 50G>C and 51G>A, show impaired binding of essential protein components of the U4atac/U6atac di-snRNP in vitro and in vivo. Analysis of MOPD I patient fibroblasts and iPS cells homozygous for the most common mutation, 51G>A, shows reduced levels of the U4atac/U6atac.U5 tri-snRNP complex as determined by glycerol gradient sedimentation and immunoprecipitation. In this report, we establish a mechanistic basis for MOPD I disease and show that the inefficient splicing of genes containing U12-dependent introns in patient cells is due to defects in minor tri-snRNP formation, and the MOPD I-associated RNU4ATAC mutations can affect multiple facets of minor snRNA function.

The significant other: splicing by the minor spliceosome

The removal of non-coding sequences, introns, from the mRNA precursors is an essential step in eukaryotic gene expression. U12-type introns are a minor subgroup of introns, distinct from the major or U2-type introns. U12-type introns are present in most eukaryotes but only account for less than 0.5% of all introns in any given genome. They are processed by a specific U12-dependent spliceosome, which is similar to, but distinct from, the major spliceosome. U12-type introns are spliced somewhat less efficiently than the major introns, and it is believed that this limits the expression of the genes containing such introns. Recent findings on the role of U12-dependent splicing in development and human disease have shown that it can also affect multiple cellular processes not directly related to the functions of the host genes of U12-type introns. At the same time, advances in understanding the regulation and phylogenetic distribution of the minor spliceosome are starting to shed light on how the U12-type introns and the minor spliceosome may have evolved. © 2012 John Wiley & Sons, Ltd.

Functionally important structural elements of U12 snRNA

U12 snRNA is analogous to U2 snRNA of the U2-dependent spliceosome and is essential for the splicing of U12-dependent introns in metazoan cells. The essential region of U12 snRNA, which base pairs to the branch site of minor class introns is well characterized. However, other regions which are outside of the branch site base pairing region are not yet characterized and the requirement of these structures in U12-dependent splicing is not clear. U12 snRNA is predicted to form an intricate secondary structure containing several stem-loops and single-stranded regions. Using a previously characterized branch site genetic suppression assay, we generated second-site mutations in the suppressor U12 snRNA to investigate the in vivo requirement of structural elements in U12-dependent splicing. Our results show that stem-loop IIa is essential and required for in vivo splicing. Interestingly, an evolutionarily conserved stem-loop IIb is dispensable for splicing. We also show that stem-loop III, which binds to a p65 RNA binding protein of the U11-U12 di.snRNP complex, is essential for in vivo splicing. The data validate the existence of proposed stem-loops of U12 snRNA and provide experimental support for individual secondary structures.

The conserved 3' end domain of U6atac snRNA can direct U6 snRNA to the minor spliceosome

U6 and U6atac snRNAs play analogous critical roles in the major U2-dependent and minor U12-dependent spliceosomes, respectively. Previous results have shown that most of the functional cores of these two snRNAs are either highly similar in sequence or functionally interchangeable. Thus, a mechanism must exist to restrict each snRNA to its own spliceosome. Here we show that a chimeric U6 snRNA containing the unique and highly conserved 3' end domain of U6atac snRNA is able to function in vivo in U12-dependent spliceosomal splicing. Function of this chimera required the coexpression of a modified U4atac snRNA; U4 snRNA could not substitute. Partial deletions of this element in vivo, as well as in vitro antisense experiments, showed that the 3' end domain of U6atac snRNA is necessary to direct the U4atac/U6atac.U5 tri-snRNP to the forming U12-dependent spliceosome. In vitro experiments also uncovered a role for U4atac snRNA in this targeting.

U12-dependent intron splicing in plants

U12-dependent (U12) introns have persisted in the genomes of plants since the ancestral divergence between plants and metazoans. These introns, which are rare, are found in a range of genes that include essential functions in DNA replication and RNA metabolism and are implicated in regulating the expression of their host genes. U12 introns are removed from pre-mRNAs by a U12 intron-specific spliceosome. Although this spliceosome shares many properties with the more abundant U2-dependent (U2) intron spliceosome, four of the five small nuclear RNAs (snRNAs) required for splicing are different and specific for the unique splicing of U12 introns. Evidence in plants so far indicates that splicing signals of plant U12 introns and their splicing machinery are similar to U12 intron splicing in other eukaryotes. In addition to the high conservation of splicing signals, plant U12 introns also retain unique characteristic features of plant U2 introns, such as UA-richness, which suggests a requirement for plant-specific components for both the U2 and U12 splicing reaction. This chapter compares U12 and U2 splicing and reviews what is known about plant U12 introns and their possible role in gene expression.

Spliceosomal snRNAs in the unicellular eukaryote Trichomonas vaginalis are structurally conserved but lack a 5'-cap structure

Few genes in the divergent eukaryote Trichomonas vaginalis have introns, despite the unusually large gene repertoire of this human-infective parasite. These introns are characterized by extended conserved regulatory motifs at the 5' and 3' boundaries, a feature shared with another divergent eukaryote, Giardia lamblia, but not with metazoan introns. This unusual characteristic of T. vaginalis introns led us to examine spliceosomal small nuclear RNAs (snRNAs) predicted to mediate splicing reactions via interaction with intron motifs. Here we identify T. vaginalis U1, U2, U4, U5, and U6 snRNAs, present predictions of their secondary structures, and provide evidence for interaction between the U2/U6 snRNA complex and a T. vaginalis intron. Structural models predict that T. vaginalis snRNAs contain conserved sequences and motifs similar to those found in other examined eukaryotes. These data indicate that mechanisms of intron recognition as well as coordination of the two catalytic steps of splicing have been conserved throughout eukaryotic evolution. Unexpectedly, we found that T. vaginalis spliceosomal snRNAs lack the 5' trimethylguanosine cap typical of snRNAs and appear to possess unmodified 5' ends. Despite the lack of a cap structure, U1, U2, U4, and U5 genes are transcribed by RNA polymerase II, whereas the U6 gene is transcribed by RNA polymerase III.

Free energy landscapes of RNA/RNA complexes: with applications to snRNA complexes in spliceosomes

We develop a statistical mechanical model for RNA/RNA complexes with both intramolecular and intermolecular interactions. As an application of the model, we compute the free energy landscapes, which give the full distribution for all the possible conformations, for U4/U6 and U2/U6 in major spliceosome and U4atac/U6atac and U12/U6atac in minor spliceosome. Different snRNA experiments found contrasting structures, our free energy landscape theory shows why these structures emerge and how they compete with each other. For yeast U2/U6, the model predicts that the two distinct experimental structures, the four-helix junction structure and the helix Ib-containing structure, can actually coexist and specifically compete with each other. In addition, the energy landscapes suggest possible mechanisms for the conformational switches in splicing. For instance, our calculation shows that coaxial stacking is essential for stabilizing the four-helix junction in yeast U2/U6. Therefore, inhibition of the coaxial stacking possibly by protein-binding may activate the conformational switch from the four-helix junction to the helix Ib-containing structure. Moreover, the change of the energy landscape shape gives information about the conformational changes. We find multiple (native-like and misfolded) intermediates formed through base-pairing rearrangements in snRNA complexes. For example, the unfolding of the U2/U6 undergoes a transition to a misfolded state which is functional, while in the unfolding of U12/U6atac, the functional helix Ib is found to be the last one to unfold and is thus the most stable structural component. Furthermore, the energy landscape gives the stabilities of all the possible (functional) intermediates and such information is directly related to splicing efficiency.

RNA structure and RNA-protein interactions in purified yeast U6 snRNPs

The U6 small nuclear RNA (snRNA) undergoes major conformational changes during the assembly of the spliceosome and catalysis of splicing. It associates with the specific protein Prp24p, and a set of seven LSm2p-8p proteins, to form the U6 small nuclear ribonucleoprotein (snRNP). These proteins have been proposed to act as RNA chaperones that stimulate pairing of U6 with U4 snRNA to form the intermolecular stem I and stem II of the U4/U6 duplex, whose formation is essential for spliceosomal function. However, the mechanism whereby Prp24p and the LSm complex facilitate U4/U6 base-pairing, as well as the exact binding site(s) of Prp24p in the native U6 snRNP, are not well understood. Here, we have investigated the secondary structure of the U6 snRNA in purified U6 snRNPs and compared it with its naked form. Using RNA structure-probing techniques, we demonstrate that within the U6 snRNP a large internal region of the U6 snRNA is unpaired and protected from chemical modification by bound Prp24p. Several of these U6 nucleotides are available for base-pairing interaction, as only their sugar backbone is contacted by Prp24p. Thus, Prp24p can present them to the U4 snRNA and facilitate formation of U4/U6 stem I. We show that the 3' stem-loop is not bound strongly by U6 proteins in native particles. However, when compared to the 3' stem-loop in the naked U6 snRNA, it has a more open conformation, which would facilitate formation of stem II with the U4 snRNA. Our data suggest that the combined association of Prp24p and the LSm complex confers upon U6 nucleotides a conformation favourable for U4/U6 base-pairing. Interestingly, we find that the open structure of the yeast U6 snRNA in native snRNPs can also be adopted by human U6 and U6atac snRNAs.

Heterogeneous nuclear ribonucleoprotein H is required for optimal U11 small nuclear ribonucleoprotein binding to a retroviral RNA-processing control element: implications for U12-dependent RNA splicing

An RNA-processing element from Rous sarcoma virus, the negative regulator of splicing (NRS), represses splicing to generate unspliced RNA that serves as mRNA and as genomic RNA for progeny virions and also promotes polyadenylation of the unspliced RNA. Integral to NRS function is the binding of U1 small nuclear ribonucleoprotein (snRNP), but its binding is controlled by U11 snRNP that binds to an overlapping site. U11 snRNP, the U1 counterpart for splicing of U12-dependent introns, binds the NRS remarkably well and requires G-rich elements just downstream of the consensus U11 binding site. We present evidence that heterogeneous nuclear ribonucleoprotein (hnRNP) H binds to the NRS G-rich elements and that hnRNP H is required for optimal U11 binding in vitro. It is further shown that hnRNP H (but not hnRNP F) can promote U11 binding and splicing from the NRS in vivo when tethered to the RNA as an MS2 fusion protein. Interestingly, 17% of the naturally occurring U12-dependent introns have at least two potential hnRNP H binding sites positioned similarly to the NRS. For two such introns from the SCN4A and P120 genes, we show that hnRNP H binds to each in a G-tract-dependent manner, that G-tract mutations strongly reduce splicing of minigene RNA, and that tethered hnRNP H restores splicing to mutant RNA. In support of a role for hnRNP H in both splicing pathways, hnRNP H antibodies co-precipitate U1 and U11 small nuclear ribonucleoproteins. These results indicate that hnRNP H is an auxiliary factor for U11 binding to the NRS and that, more generally, hnRNP H is a splicing factor for a subset of U12-dependent introns that harbor G-rich elements.

Splicing of a rare class of introns by the U12-dependent spliceosome

Pre-mRNA splicing is catalyzed by two unique spliceosomes, designated U2- or U12-dependent. In contrast to the well-characterized U2-dependent spliceosome, much remains to be learned about the less abundant U12-type spliceosome. This review focuses on recent advances in elucidating the structure and function of the minor U12-dependent spliceosome. Interesting similarities and differences between the U12- and U2-dependent spliceosomes are also highlighted.

Determinants of plant U12-dependent intron splicing efficiency

Factors affecting splicing of plant U12-dependent introns have been examined by extensive mutational analyses in an in vivo tobacco (Nicotiana tabacum) protoplast system using introns from three different Arabidopsis thaliana genes: CBP20, GSH2, and LD. The results provide evidence that splicing efficiency of plant U12 introns depends on a combination of factors, including UA content, exon bridging interactions between the U12 intron and flanking U2-dependent introns, and exon splicing enhancer sequences (ESEs). Unexpectedly, all three plant U12 introns required an adenosine at the upstream purine position in the branchpoint consensus UCCUURAUY. The exon upstream of the LD U12 intron is a major determinant of its higher level of splicing efficiency and potentially contains two ESE regions. These results suggest that in plants, U12 introns represent a level at which expression of their host genes can be regulated.

U4 small nuclear RNA can function in both the major and minor spliceosomes

U4 small nuclear RNA (snRNA) and U6 snRNA form a base-paired di-snRNP complex that is essential for pre-mRNA splicing of the major class of metazoan nuclear introns. The functionally analogous but highly diverged U4atac and U6atac snRNAs form a similar complex that is involved in splicing of the minor class of introns. Previous results with mutants of U6atac in which a substructure was replaced by the analogous structure from U6 snRNA suggested that wild-type U4 snRNA might be able to interact productively with the mutant U6atac snRNA. Here we show that a mutant U4 snRNA designed to base pair with a mutant U6atac snRNA can activate U12-dependent splicing when coexpressed in an in vivo genetic suppression assay. This genetic interaction could also be demonstrated in an in vitro crosslinking assay. These results show that a U4/U6atac di-snRNP can correctly splice a U12-dependent intron and suggest that the specificity for spliceosome type resides in the U6 and U6atac small nuclear ribonucleoproteins. Further experiments suggest that expression of a mutant U4 snRNA that can bind to wild-type U6atac snRNA alters the specificity of some splice sites, providing an evolutionary rationale for maintaining two U4-like snRNAs.

Identification, characterization and molecular phylogeny of U12-dependent introns in the Arabidopsis thaliana genome

U12-dependent introns are spliced by the minor U12-type spliceosome and occur in a variety of eukaryotic organisms, including Arabidopsis. In this study, a set of putative U12-dependent introns was compiled from a large collection of cDNA/EST- confirmed introns in the Arabidopsis thaliana genome by means of high-throughput bioinformatic analysis combined with manual scrutiny. A total of 165 U12-type introns were identified based upon stringent criteria. This number of sequences well exceeds the total number of U12-type introns previously reported for plants and allows a more thorough statistical analysis of U12-type signals. Of particular note is the discovery that the distance between the branch site adenosine and the acceptor site ranges from 10 to 39 nt, significantly longer than the previously postulated limit of 21 bp. Further analysis indicates that, in addition to the spacing constraint, the sequence context of the potential acceptor site may have an important role in 3' splice site selection. Several alternative splicing events involving U12-type introns were also captured in this study, providing evidence that U12-dependent acceptor sites can also be recognized by the U2-type spliceosome. Furthermore, phylogenetic analysis suggests that both U12-type AT-AC and U12-type GT-AG introns occurred in Na+/H+ antiporters in a progenitor of animals and plants.

Domains of human U4atac snRNA required for U12-dependent splicing in vivo

U4atac snRNA forms a base-paired complex with U6atac snRNA. Both snRNAs are required for the splicing of the minor U12-dependent class of eukaryotic nuclear introns. We have developed a new genetic suppression assay to investigate the in vivo roles of several regions of U4atac snRNA in U12-dependent splicing. We show that both the stem I and stem II regions, which have been proposed to pair with U6atac snRNA, are required for in vivo splicing. Splicing activity also requires U4atac sequences in the 5' stem-loop element that bind a 15.5 kDa protein that also binds to a similar region of U4 snRNA. In contrast, mutations in the region immediately following the stem I interaction region, as well as a deletion of the distal portion of the 3' stem-loop element, were active for splicing. Complete deletion of the 3' stem-loop element abolished in vivo splicing function as did a mutation of the Sm protein binding site. These results show that the in vivo sequence requirements of U4atac snRNA are similar to those described previously for U4 snRNA using in vitro assays and provide experimental support for models of the U4atac/U6atac snRNA interaction.

The conserved central domain of yeast U6 snRNA: importance of U2-U6 helix Ia in spliceosome assembly

In the pre-mRNA processing machinery of eukaryotic cells, U6 snRNA is located at or near the active site for pre-mRNA splicing catalysis, and U6 is involved in catalyzing the first chemical step of splicing. We have further defined the roles of key features of yeast U6 snRNA in the splicing process. By assaying spliceosome assembly and splicing in yeast extracts, we found that mutations of yeast U6 nt 56 and 57 are similar to previously reported deletions of U2 nt 27 or 28, all within yeast U2-U6 helix Ia. These mutations lead to the accumulation of yeast A1 spliceosomes, which form just prior to the Prp2 ATPase step and the first chemical step of splicing. These results strongly suggest that, at a late stage of spliceosome assembly, the presence of U2-U6 helix Ia is important for promoting the first chemical step of splicing, presumably by bringing together the 5' splice site region of pre-mRNA, which is base paired to U6 snRNA, and the branchsite region of the intron, which is base paired to U2 snRNA, for activation of the first chemical step of splicing, as previously proposed by Madhani and Guthrie [Cell, 1992, 71: 803-817]. In the 3' intramolecular stem-loop of U6, mutation G81C causes an allele-specific accumulation of U6 snRNP. Base pairing of the U6 3' stem-loop in yeast spliceosomes does not extend as far as to include the U6 sequence of U2-U6 helix Ib, in contrast to the human U6 3' stem-loop structure.

A catalytically active group II intron domain 5 can function in the U12-dependent spliceosome

Both spliceosomal and self-splicing group II introns require the function of similar small, metal binding RNA stem-loop elements located in U6 or U6atac snRNAs of the spliceosome or domain 5 (D5) of group II introns. Here we report that two different D5 elements can functionally replace the U6atac snRNA stem-loop in an in vivo splicing assay. For efficient function in vivo, a single base pair from the upper helical section of the D5 sequence had to be removed. Introducing the equivalent base pair deletion into the D5 element of a group II intron reduced but did not eliminate self-splicing activity. Our results strengthen the case that these RNA elements play similar roles in the catalytic centers of both the spliceosome and a self-splicing ribozyme.

Human U4/U6.U5 and U4atac/U6atac.U5 tri-snRNPs exhibit similar protein compositions

In the U12-dependent spliceosome, the U4atac/U6atac snRNP represents the functional analogue of the major U4/U6 snRNP. Little information is available presently regarding the protein composition of the former snRNP and its association with other snRNPs. In this report we show that human U4atac/U6atac di-snRNPs associate with U5 snRNPs to form a 25S U4atac/U6atac.U5 trimeric particle. Comparative analysis of minor and major tri-snRNPs by using immunoprecipitation experiments revealed that their protein compositions are very similar, if not identical. Not only U5-specific proteins but, surprisingly, all tested U4/U6- and major tri-snRNP-specific proteins were detected in the minor tri-snRNP complex. Significantly, the major tri-snRNP-specific proteins 65K and 110K, which are required for integration of the major tri-snRNP into the U2-dependent spliceosome, were among those proteins detected in the minor tri-snRNP, raising an interesting question as to how the specificity of addition of tri-snRNP to the corresponding spliceosome is maintained. Moreover, immunodepletion studies demonstrated that the U4/U6-specific 61K protein, which is involved in the formation of major tri-snRNPs, is essential for the association of the U4atac/U6atac di-snRNP with U5 to form the U4atac/U6atac.U5 tri-snRNP. Subsequent immunoprecipitation studies demonstrated that those proteins detected in the minor tri-snRNP complex are also incorporated into U12-dependent spliceosomes. This remarkable conservation of polypeptides between minor and major spliceosomes, coupled with the absence of significant sequence similarity between the functionally analogous snRNAs, supports an evolutionary model in which most major and minor spliceosomal proteins, but not snRNAs, are derived from a common ancestor.

A computational scan for U12-dependent introns in the human genome sequence

U12-dependent introns are found in small numbers in most eukaryotic genomes, but their scarcity makes accurate characterisation of their properties challenging. A computational search for U12-dependent introns was performed using the draft version of the human genome sequence. Human expressed sequences confirmed 404 U12-dependent introns within the human genome, a 6-fold increase over the total number of non-redundant U12-dependent introns previously identified in all genomes. Although most of these introns had AT-AC or GT-AG terminal dinucleotides, small numbers of introns with a surprising diversity of termini were found, suggesting that many of the non-canonical introns found in the human genome may be variants of U12-dependent introns and, thus, spliced by the minor spliceosome. Comparisons with U2-dependent introns revealed that the U12-dependent intron set lacks the 'short intron' peak characteristic of U2-dependent introns. Analysis of this U12-dependent intron set confirmed reports of a biased distribution of U12-dependent introns in the genome and allowed the identification of several alternative splicing events as well as a surprising number of apparent splicing errors. This new larger reference set of U12-dependent introns will serve as a resource for future studies of both the properties and evolution of the U12 spliceosome.