hnRNPM protects against the dsRNA-mediated interferon response by repressing LINE-associated cryptic splicing.

RNA splicing is pivotal in post-transcriptional gene regulation, yet the exponential expansion of intron length in humans poses a challenge for accurate splicing. Here, we identify hnRNPM as an essential RNA-binding protein that suppresses cryptic splicing through binding to deep introns, maintaining human transcriptome integrity. Long interspersed nuclear elements (LINEs) in introns harbor numerous pseudo splice sites. hnRNPM preferentially binds at intronic LINEs to repress pseudo splice site usage for cryptic splicing. Remarkably, cryptic exons can generate long dsRNAs through base-pairing of inverted ALU transposable elements interspersed among LINEs and consequently trigger an interferon response, a well-known antiviral defense mechanism. Significantly, hnRNPM-deficient tumors show upregulated interferon-associated pathways and elevated immune cell infiltration. These findings unveil hnRNPM as a guardian of transcriptome integrity by repressing cryptic splicing and suggest that targeting hnRNPM in tumors may be used to trigger an inflammatory immune response, thereby boosting cancer surveillance.

Interpretable prioritization of splice variants in diagnostic next-generation sequencing.

A critical challenge in genetic diagnostics is the computational assessment of candidate splice variants, specifically the interpretation of nucleotide changes located outside of the highly conserved dinucleotide sequences at the 5' and 3' ends of introns. To address this gap, we developed the Super Quick Information-content Random-forest Learning of Splice variants (SQUIRLS) algorithm. SQUIRLS generates a small set of interpretable features for machine learning by calculating the information-content of wild-type and variant sequences of canonical and cryptic splice sites, assessing changes in candidate splicing regulatory sequences, and incorporating characteristics of the sequence such as exon length, disruptions of the AG exclusion zone, and conservation. We curated a comprehensive collection of disease-associated splice-altering variants at positions outside of the highly conserved AG/GT dinucleotides at the termini of introns. SQUIRLS trains two random-forest classifiers for the donor and for the acceptor and combines their outputs by logistic regression to yield a final score. We show that SQUIRLS transcends previous state-of-the-art accuracy in classifying splice variants as assessed by rank analysis in simulated exomes, and is significantly faster than competing methods. SQUIRLS provides tabular output files for incorporation into diagnostic pipelines for exome and genome analysis, as well as visualizations that contextualize predicted effects of variants on splicing to make it easier to interpret splice variants in diagnostic settings.

RNA aptamer reveals nuclear TDP-43 pathology is an early aggregation event that coincides with STMN-2 cryptic splicing and precedes clinical manifestation in ALS.

TDP-43 is an aggregation-prone protein which accumulates in the hallmark pathological inclusions of amyotrophic lateral sclerosis (ALS). However, the analysis of deeply phenotyped human post-mortem samples has shown that TDP-43 aggregation, revealed by standard antibody methods, correlates poorly with symptom manifestation. Recent identification of cryptic-splicing events, such as the detection of Stathmin-2 (STMN-2) cryptic exons, are providing evidence implicating TDP-43 loss-of-function as a potential driving pathomechanism but the temporal nature of TDP-43 loss and its relation to the disease process and clinical phenotype is not known. To address these outstanding questions, we used a novel RNA aptamer, TDP-43APT, to detect TDP-43 pathology and used single molecule in situ hybridization to sensitively reveal TDP-43 loss-of-function and applied these in a deeply phenotyped human post-mortem tissue cohort. We demonstrate that TDP-43APT identifies pathological TDP-43, detecting aggregation events that cannot be detected by classical antibody stains. We show that nuclear TDP-43 pathology is an early event, occurring prior to cytoplasmic accumulation and is associated with loss-of-function measured by coincident STMN-2 cryptic splicing pathology. Crucially, we show that these pathological features of TDP-43 loss-of-function precede the clinical inflection point and are not required for region specific clinical manifestation. Furthermore, we demonstrate that gain-of-function in the form of extensive cytoplasmic accumulation, but not loss-of-function, is the primary molecular correlate of clinical manifestation. Taken together, our findings demonstrate implications for early diagnostics as the presence of STMN-2 cryptic exons and early TDP-43 aggregation events could be detected prior to symptom onset, holding promise for early intervention in ALS.

Combining genetic constraint with predictions of alternative splicing to prioritize deleterious splicing in rare disease studies.

BACKGROUND: Despite numerous molecular and computational advances, roughly half of patients with a rare disease remain undiagnosed after exome or genome sequencing. A particularly challenging barrier to diagnosis is identifying variants that cause deleterious alternative splicing at intronic or exonic loci outside of canonical donor or acceptor splice sites. RESULTS: Several existing tools predict the likelihood that a genetic variant causes alternative splicing. We sought to extend such methods by developing a new metric that aids in discerning whether a genetic variant leads to deleterious alternative splicing. Our metric combines genetic variation in the Genome Aggregate Database with alternative splicing predictions from SpliceAI to compare observed and expected levels of splice-altering genetic variation. We infer genic regions with significantly less splice-altering variation than expected to be constrained. The resulting model of regional splicing constraint captures differential splicing constraint across gene and exon categories, and the most constrained genic regions are enriched for pathogenic splice-altering variants. Building from this model, we developed ConSpliceML. This ensemble machine learning approach combines regional splicing constraint with multiple per-nucleotide alternative splicing scores to guide the prediction of deleterious splicing variants in protein-coding genes. ConSpliceML more accurately distinguishes deleterious and benign splicing variants than state-of-the-art splicing prediction methods, especially in "cryptic" splicing regions beyond canonical donor or acceptor splice sites. CONCLUSION: Integrating a model of genetic constraint with annotations from existing alternative splicing tools allows ConSpliceML to prioritize potentially deleterious splice-altering variants in studies of rare human diseases.

Prp8 impacts cryptic but not alternative splicing frequency.

Pre-mRNA splicing must occur with extremely high fidelity. Spliceosomes assemble onto pre-mRNA guided by specific sequences (5' splice site, 3' splice site, and branchpoint). When splice sites are mutated, as in many hereditary diseases, the spliceosome can aberrantly select nearby pseudo- or "cryptic" splice sites, often resulting in nonfunctional protein. How the spliceosome distinguishes authentic splice sites from cryptic splice sites is poorly understood. We performed a Caenorhabditis elegans genetic screen to find cellular factors that affect the frequency with which the spliceosome uses cryptic splice sites and identified two alleles in core spliceosome component Prp8 that alter cryptic splicing frequency. Subsequent complementary genetic and structural analyses in yeast implicate these alleles in the stability of the spliceosome's catalytic core. However, despite a clear effect on cryptic splicing, high-throughput mRNA sequencing of these prp-8 mutant C. elegans reveals that overall alternative splicing patterns are relatively unchanged. Our data suggest the spliceosome evolved intrinsic mechanisms to reduce the occurrence of cryptic splicing and that these mechanisms are distinct from those that impact alternative splicing.

Systematic Computational Identification of Variants That Activate Exonic and Intronic Cryptic Splice Sites.

We developed a variant-annotation method that combines sequence-based machine-learning classification with a context-dependent algorithm for selecting splice variants. Our approach is distinctive in that it compares the splice potential of a sequence bearing a variant with the splice potential of the reference sequence. After training, classification accurately identified 168 of 180 (93.3%) canonical splice sites of five genes. The combined method, CryptSplice, identified and correctly predicted the effect of 18 of 21 (86%) known splice-altering variants in CFTR, a well-studied gene whose loss-of-function variants cause cystic fibrosis (CF). Among 1,423 unannotated CFTR disease-associated variants, the method identified 32 potential exonic cryptic splice variants, two of which were experimentally evaluated and confirmed. After complete CFTR sequencing, the method found three cryptic intronic splice variants (one known and two experimentally verified) that completed the molecular diagnosis of CF in 6 of 14 individuals. CryptSplice interrogation of sequence data from six individuals with X-linked dyskeratosis congenita caused by an unknown disease-causing variant in DKC1 identified two splice-altering variants that were experimentally verified. To assess the extent to which disease-associated variants might activate cryptic splicing, we selected 458 pathogenic variants and 348 variants of uncertain significance (VUSs) classified as high confidence from ClinVar. Splice-site activation was predicted for 129 (28%) of the pathogenic variants and 75 (22%) of the VUSs. Our findings suggest that cryptic splice-site activation is more common than previously thought and should be routinely considered for all variants within the transcribed regions of genes.

Celia's encephalopathy and c.974dupG in BSCL2 gene: a hidden change in a known variant.

Celia's encephalopathy (progressive encephalopathy with/without lipodystrophy (PELD)) is a childhood neurodegenerative disorder with a fatal prognosis before the age of 10, due to the variant c.985C>T in the BSCL2 gene that causes a cryptic splicing site leading to skipping of exon 7. For years, different authors have reported cases of congenital generalized lipodystrophy due to the variant c.974dupG in BSCL2 associated with neurological manifestations of variable severity, although some of them clearly superimposable to PELD. To identify the molecular mechanisms responsible for these neurological alterations in two patients with c.974dupG. Clinical characterization, biochemistry, and neuroimaging studies of two girls carrying this variant. In silico analysis, PCR amplification, and BSCL2 cDNA sequencing. BSCL2-201 transcript expression, which lacks exon 7, by qPCR in fibroblasts from the index case, from a healthy child as a control and from two patients with PELD, and in leukocytes from the index case and her parents. One with a severe encephalopathy including a picture of intellectual deficiency, severe language impairment, myoclonic epilepsy, and lipodystrophy as described in PELD, dying at 9 years and 9 months of age. The other 2-year-old patient showed incipient signs of neurological involvement. In silico and cDNA sequencing studies showed that variant c.974dupG gives rise to skipping of exon 7. The expression of BSCL2-201 in fibroblasts was significantly higher in the index case than in the healthy child, although less than in the case with homozygous PELD due to c.985C>T variant. The expression of this transcript was approximately half in the healthy carrier parents of this patient. The c.974dupG variant leads to the skipping of exon 7 of the BSCL2 gene and is responsible for a variant of Celia's encephalopathy, with variable phenotypic expression.

Inhibition of nonsense-mediated RNA decay by ER stress.

Nonsense-mediated RNA decay (NMD) selectively degrades mutated and aberrantly processed transcripts that contain premature termination codons (PTC). Cellular NMD activity is typically assessed using exogenous PTC-containing reporters. We overcame some inherently problematic aspects of assaying endogenous targets and developed a broadly applicable strategy to reliably and easily monitor changes in cellular NMD activity. Our new method was genetically validated for distinguishing NMD regulation from transcriptional control and alternative splicing regulation, and unexpectedly disclosed a different sensitivity of NMD targets to NMD inhibition. Applying this robust method for screening, we identified NMD-inhibiting stressors but also found that NMD inactivation was not universal to cellular stresses. The high sensitivity and broad dynamic range of our method revealed a strong correlation between NMD inhibition, endoplasmic reticulum (ER) stress, and polysome disassembly upon thapsigargin treatment in a temporal and dose-dependent manner. We found little evidence of calcium signaling mediating thapsigargin-induced NMD inhibition. Instead, we discovered that of the three unfolded protein response (UPR) pathways activated by thapsigargin, mainly protein kinase RNA-like endoplasmic reticulum kinase (PERK) was required for NMD inhibition. Finally, we showed that ER stress compounded TDP-43 depletion in the up-regulation of NMD isoforms that had been implicated in the pathogenic mechanisms of amyotrophic lateral sclerosis and frontotemporal dementia, and that the additive effect of ER stress was completely blocked by PERK deficiency.

Transcriptome protection by the expanded family of hnRNPs.

The family of heterogeneous ribonucleoproteins (hnRNPs) have multiple functions in RNA metabolism. In recent years, several hnRNPs have also been shown to be essential for the maintenance of transcriptome integrity, by preventing intronic cryptic splicing signals from mis-splicing of many endogeneous pre-mRNA transcripts. Here we discuss the possibility for a general role of this family of proteins and their expansion in transcriptome protection.

MATR3 pathogenic variants differentially impair its cryptic splicing repression function.

Matrin-3 (MATR3) is an RNA-binding protein implicated in neurodegenerative and neurodevelopmental diseases. However, little is known regarding the role of MATR3 in cryptic splicing within the context of functional genes and how disease-associated variants impact this function. We show that loss of MATR3 leads to cryptic exon inclusion in many transcripts. We reveal that ALS-linked S85C pathogenic variant reduces MATR3 solubility but does not impair RNA binding. In parallel, we report a novel neurodevelopmental disease-associated M548T variant, located in the RRM2 domain, which reduces protein solubility and impairs RNA binding and cryptic splicing repression functions of MATR3. Altogether, our research identifies cryptic events within functional genes and demonstrates how disease-associated variants impact MATR3 cryptic splicing repression function.

Regulating PCCA gene expression by modulation of pseudoexon splicing patterns to rescue enzyme activity in propionic acidemia.

Pseudoexons are nonfunctional intronic sequences that can be activated by deep-intronic sequence variation. Activation increases pseudoexon inclusion in mRNA and interferes with normal gene expression. The PCCA c.1285-1416A>G variation activates a pseudoexon and causes the severe metabolic disorder propionic acidemia by deficiency of the propionyl-CoA carboxylase enzyme encoded by PCCA and PCCB. We characterized this pathogenic pseudoexon activation event in detail and identified hnRNP A1 to be important for normal repression. The PCCA c.1285-1416A>G variation disrupts an hnRNP A1-binding splicing silencer and simultaneously creates a splicing enhancer. We demonstrate that blocking this region of regulation with splice-switching antisense oligonucleotides restores normal splicing and rescues enzyme activity in patient fibroblasts and in a cellular model created by CRISPR gene editing. Interestingly, the PCCA pseudoexon offers an unexploited potential to upregulate gene expression because healthy tissues show relatively high inclusion levels. By blocking inclusion of the nonactivated wild-type pseudoexon, we can increase both PCCA and PCCB protein levels, which increases the activity of the heterododecameric enzyme. Surprisingly, we can increase enzyme activity from residual levels in not only patient fibroblasts harboring PCCA missense variants but also those harboring PCCB missense variants. This is a potential treatment strategy for propionic acidemia.

Identification of a novel TNNI3 synonymous variant causing intron retention in autosomal recessive dilated cardiomyopathy.

BACKGROUND: Dilated cardiomyopathy type 2A (DCM2A, MIM: #611880) is a rare autosomal recessive heart disease leading to heart failure and sudden cardiac death. However, the causative role of TNNI3 in DCM2A is still questioned due to few cases reported and the conflicting molecular biological evidence. METHODS: Trio whole-exome sequencing (trio-WES) was performed in a Chinese family with dilated cardiomyopathy. Sanger sequencing and real-time quantitative PCR were used to confirm the variants identified. Expression outcome caused by the synonymous mutation was validated by minigene splicing analyses. RESULTS: The one-year-old girl presented severe left ventricular enlargement and significantly reduced left ventricular systolic function and she died of respiratory and heart failure soon after her diagnosis. Trio-WES revealed a compound heterozygous variants of TNNI3, a novel c.24G＞A (p.Ala8Ala) (NM_000363.4) in exon 2 and a deletion of entire gene. Minigene splicing analyses showed it led to an intron retention (c.24 + 1_24 + 45ins) by intron 2 cryptic splicing. CONCLUSIONS: Our study describes and characterizes a synonymous mutation in TNNI3 gene, supporting the clinical diagnosis of an autosomal recessive DCM. Our study emphasizes the importance of functional analysis to assess the potential pathogenicity of synonymous mutations, especially when the synonymous variants are not annotated as benign.

Splicing quality control mediated by DHX15 and its G-patch activator SUGP1.

Pre-mRNA splicing is surveilled at different stages by quality control (QC) mechanisms. The leukemia-associated DExH-box family helicase hDHX15/scPrp43 is known to disassemble spliceosomes after splicing. Here, using rapid protein depletion and analysis of nascent and mature RNA to enrich for direct effects, we identify a widespread splicing QC function for DHX15 in human cells, consistent with recent in vitro studies. We find that suboptimal introns with weak splice sites, multiple branch points, and cryptic introns are repressed by DHX15, suggesting a general role in promoting splicing fidelity. We identify SUGP1 as a G-patch factor that activates DHX15's splicing QC function. This interaction is dependent on both DHX15's ATPase activity and on SUGP1's U2AF ligand motif (ULM) domain. Together, our results support a model in which DHX15 plays a major role in splicing QC when recruited and activated by SUGP1.

Cryptic splicing: common pathological mechanisms involved in male infertility and neuronal diseases.

High levels of transcription and alternative splicing are recognized hallmarks of gene expression in the testis and largely driven by cells in meiosis. Because of this, the male meiosis stage of the cell cycle is often viewed as having a relatively permissive environment for gene expression. In this review, we highlight recent findings that identify the RNA binding protein RBMXL2 as essential for male meiosis. RBMXL2 functions as a "guardian of the transcriptome" that protects against the use of aberrant (or "cryptic") splice sites that would disrupt gene expression. This newly discovered protective role during meiosis links with a wider field investigating mechanisms of cryptic splicing control that protect neurons from amyotrophic lateral sclerosis and Alzheimer's disease. We discuss how the mechanism repressing cryptic splicing patterns during meiosis evolved, and why it may be essential for sperm production and male fertility.

Corrigendum: Analysis of Pathogenic Pseudoexons Reveals Novel Mechanisms Driving Cryptic Splicing.

[This corrects the article DOI: 10.3389/fgene.2021.806946.].

Analysis of Pathogenic Pseudoexons Reveals Novel Mechanisms Driving Cryptic Splicing.

Understanding pre-mRNA splicing is crucial to accurately diagnosing and treating genetic diseases. However, mutations that alter splicing can exert highly diverse effects. Of all the known types of splicing mutations, perhaps the rarest and most difficult to predict are those that activate pseudoexons, sometimes also called cryptic exons. Unlike other splicing mutations that either destroy or redirect existing splice events, pseudoexon mutations appear to create entirely new exons within introns. Since exon definition in vertebrates requires coordinated arrangements of numerous RNA motifs, one might expect that pseudoexons would only arise when rearrangements of intronic DNA create novel exons by chance. Surprisingly, although such mutations do occur, a far more common cause of pseudoexons is deep-intronic single nucleotide variants, raising the question of why these latent exon-like tracts near the mutation sites have not already been purged from the genome by the evolutionary advantage of more efficient splicing. Possible answers may lie in deep intronic splicing processes such as recursive splicing or poison exon splicing. Because these processes utilize intronic motifs that benignly engage with the spliceosome, the regions involved may be more susceptible to exonization than other intronic regions would be. We speculated that a comprehensive study of reported pseudoexons might detect alignments with known deep intronic splice sites and could also permit the characterisation of novel pseudoexon categories. In this report, we present and analyse a catalogue of over 400 published pseudoexon splice events. In addition to confirming prior observations of the most common pseudoexon mutation types, the size of this catalogue also enabled us to suggest new categories for some of the rarer types of pseudoexon mutation. By comparing our catalogue against published datasets of non-canonical splice events, we also found that 15.7% of pseudoexons exhibit some splicing activity at one or both of their splice sites in non-mutant cells. Importantly, this included seven examples of experimentally confirmed recursive splice sites, confirming for the first time a long-suspected link between these two splicing phenomena. These findings have the potential to improve the fidelity of genetic diagnostics and reveal new targets for splice-modulating therapies.

The CDH1 c.1901C>T Variant: A Founder Variant in the Portuguese Population with Severe Impact in mRNA Splicing.

Hereditary diffuse gastric cancer (HDGC) caused by CDH1 variants predisposes to early-onset diffuse gastric (DGC) and lobular breast cancer (LBC). In Northern Portugal, the unusually high number of HDGC cases in unrelated families carrying the c.1901C>T variant (formerly known as p.A634V) suggested this as a CDH1-founder variant. We aimed to demonstrate that c.1901C>T is a bona fide truncating variant inducing cryptic splicing, to calculate the timing of a potential founder effect, and to characterize tumour spectrum and age of onset in carrying families. The impact in splicing was proven by using carriers' RNA for PCR-cloning sequencing and allelic expression imbalance analysis with SNaPshot. Carriers and noncarriers were haplotyped for 12 polymorphic markers, and the decay of haplotype sharing (DHS) method was used to estimate the time to the most common ancestor of c.1901C>T. Clinical information from 58 carriers was collected and analysed. We validated the cryptic splice site within CDH1-exon 12, which was preferred over the canonical one in 100% of sequenced clones. Cryptic splicing induced an out-of-frame 37bp deletion in exon 12, premature truncation (p.Ala634ProfsTer7), and consequently RNA mediated decay. The haplotypes carrying the c.1901C>T variant were found to share a common ancestral estimated at 490 years (95% Confidence Interval 445-10,900). Among 58 carriers (27 males (M)-31 females (F); 13-83 years), DGC occurred in 11 (18.9%; 4M-7F; average age 33 ± 12) and LBC in 6 females (19.4%; average age 50 ± 8). Herein, we demonstrated that the c.1901C>T variant is a loss-of-function splice-site variant that underlies the first CDH1-founder effect in Portugal. Knowledge on this founder effect will drive genetic testing of this specific variant in HDGC families in this geographical region and allow intrafamilial penetrance analysis and better estimation of variant-associated tumour risks, disease age of onset, and spectrum.

At the Intersection of Major and Minor Spliceosomes: Crosstalk Mechanisms and Their Impact on Gene Expression.

Many eukaryotic species contain two separate molecular machineries for removing non-coding intron sequences from pre-mRNA molecules. The majority of introns (more than 99.5% in humans) are recognized and excised by the major spliceosome, which utilizes relatively poorly conserved sequence elements at the 5' and 3' ends of the intron that are used for intron recognition and in subsequent catalysis. In contrast, the minor spliceosome targets a rare group of introns (approximately 0.5% in humans) with highly conserved sequences at the 5' and 3' ends of the intron. Minor introns coexist in the same genes with major introns and while the two intron types are spliced by separate spliceosomes, the two splicing machineries can interact with one another to shape mRNA processing events in genes containing minor introns. Here, we review known cooperative and competitive interactions between the two spliceosomes and discuss the mechanistic basis of the spliceosome crosstalk, its regulatory significance, and impact on spliceosome diseases.

Alkaptonuria in Russia: mutational spectrum and novel variants.

Alkaptonuria is a rare genetic disease caused by mutations in HGD gene. Here we report the results of genetic and biochemical analysis of 49 Russian patients with alkaptonuria. One of the common variants c.481G > A; p.(Gly161Arg) comprising 72.4% of identified alleles was found in 45 of 49 patients in our cohort, which is probably the highest frequency of this variant worldwide. 9 novel variants were found: 6 missense, 2 splicing and 1 loss of start-codon. For missense variants we performed bioinformatic analysis, protein 3D-modeling and molecular dynamics simulations, which strongly suggest their pathogenic effect. For the rare synonymous variant c.753C > T; p.(Gly251Gly), which was found in 3 cases and predicted to activate cryptic splice site, we performed the detailed functional analysis on patient's cDNA and minigene assay and confirmed its pathogenicity.

Expression Changes Confirm Genomic Variants Predicted to Result in Allele-Specific, Alternative mRNA Splicing.

Splice isoform structure and abundance can be affected by either noncoding or masquerading coding variants that alter the structure or abundance of transcripts. When these variants are common in the population, these nonconstitutive transcripts are sufficiently frequent so as to resemble naturally occurring, alternative mRNA splicing. Prediction of the effects of such variants has been shown to be accurate using information theory-based methods. Single nucleotide polymorphisms (SNPs) predicted to significantly alter natural and/or cryptic splice site strength were shown to affect gene expression. Splicing changes for known SNP genotypes were confirmed in HapMap lymphoblastoid cell lines with gene expression microarrays and custom designed q-RT-PCR or TaqMan assays. The majority of these SNPs (15 of 22) as well as an independent set of 24 variants were then subjected to RNAseq analysis using the ValidSpliceMut web beacon (http://validsplicemut.cytognomix.com), which is based on data from the Cancer Genome Atlas and International Cancer Genome Consortium. SNPs from different genes analyzed with gene expression microarray and q-RT-PCR exhibited significant changes in affected splice site use. Thirteen SNPs directly affected exon inclusion and 10 altered cryptic site use. Homozygous SNP genotypes resulting in stronger splice sites exhibited higher levels of processed mRNA than alleles associated with weaker sites. Four SNPs exhibited variable expression among individuals with the same genotypes, masking statistically significant expression differences between alleles. Genome-wide information theory and expression analyses (RNAseq) in tumor exomes and genomes confirmed splicing effects for 7 of the HapMap SNP and 14 SNPs identified from tumor genomes. q-RT-PCR resolved rare splice isoforms with read abundance too low for statistical significance in ValidSpliceMut. Nevertheless, the web-beacon provides evidence of unanticipated splicing outcomes, for example, intron retention due to compromised recognition of constitutive splice sites. Thus, ValidSpliceMut and q-RT-PCR represent complementary resources for identification of allele-specific, alternative splicing.