An axonemal intron splicing program sustains Plasmodium male development.

Differentiation of male gametocytes into flagellated fertile male gametes relies on the assembly of axoneme, a major component of male development for mosquito transmission of the malaria parasite. RNA-binding protein (RBP)-mediated post-transcriptional regulation of mRNA plays important roles in eukaryotic sexual development, including the development of female Plasmodium. However, the role of RBP in defining the Plasmodium male transcriptome and its function in male gametogenesis remains incompletely understood. Here, we performed genome-wide screening for gender-specific RBPs and identified an undescribed male-specific RBP gene Rbpm1 in the Plasmodium. RBPm1 is localized in the nucleus of male gametocytes. RBPm1-deficient parasites fail to assemble the axoneme for male gametogenesis and thus mosquito transmission. RBPm1 interacts with the spliceosome E complex and regulates the splicing initiation of certain introns in a group of 26 axonemal genes. RBPm1 deficiency results in intron retention and protein loss of these axonemal genes. Intron deletion restores axonemal protein expression and partially rectifies axonemal defects in RBPm1-null gametocytes. Further splicing assays in both reporter and endogenous genes exhibit stringent recognition of the axonemal introns by RBPm1. The splicing activator RBPm1 and its target introns constitute an axonemal intron splicing program in the post-transcriptional regulation essential for Plasmodium male development.

Age-dependent regulation of ELP1 exon 20 splicing in Familial Dysautonomia by RNA Polymerase II kinetics and chromatin structure.

Familial Dysautonomia (FD) is a rare disease caused by ELP1 exon 20 skipping. Here we clarify the role of RNA Polymerase II (RNAPII) and chromatin on this splicing event. A slow RNAPII mutant and chromatin-modifying chemicals that reduce the rate of RNAPII elongation induce exon skipping whereas chemicals that create a more relaxed chromatin exon inclusion. In the brain of a mouse transgenic for the human FD-ELP1 we observed on this gene an age-dependent decrease in the RNAPII density profile that was most pronounced on the alternative exon, a robust increase in the repressive marks H3K27me3 and H3K9me3 and a decrease of H3K27Ac, together with a progressive reduction in ELP1 exon 20 inclusion level. In HEK 293T cells, selective drug-induced demethylation of H3K27 increased RNAPII elongation on ELP1 and SMN2, promoted the inclusion of the corresponding alternative exons, and, by RNA-sequencing analysis, induced changes in several alternative splicing events. These data suggest a co-transcriptional model of splicing regulation in which age-dependent changes in H3K27me3/Ac modify the rate of RNAPII elongation and affect processing of ELP1 alternative exon 20.

Elevated nuclear TDP-43 induces constitutive exon skipping.

BACKGROUND: Cytoplasmic inclusions and loss of nuclear TDP-43 are key pathological features found in several neurodegenerative disorders, suggesting both gain- and loss-of-function mechanisms of disease. To study gain-of-function, TDP-43 overexpression has been used to generate in vitro and in vivo model systems. METHODS: We analyzed RNA-seq datasets from mouse and human neurons overexpressing TDP-43 to explore species specific splicing patterns. We explored the dynamics between TDP-43 levels and exon repression in vitro. Furthermore we analyzed human brain samples and publicly available RNA datasets to explore the relationship between exon repression and disease. RESULTS: Our study shows that excessive levels of nuclear TDP-43 protein lead to constitutive exon skipping that is largely species-specific. Furthermore, while aberrant exon skipping is detected in some human brains, it is not correlated with disease, unlike the incorporation of cryptic exons that occurs after loss of TDP-43. CONCLUSIONS: Our findings emphasize the need for caution in interpreting TDP-43 overexpression data and stress the importance of controlling for exon skipping when generating models of TDP-43 proteinopathy.

LINC00982-encoded protein PRDM16-DT regulates CHEK2 splicing to suppress colorectal cancer metastasis and chemoresistance.

Metastasis is one of the key factors of treatment failure in late-stage colorectal cancer (CRC). Metastatic CRC frequently develops resistance to chemotherapeutic agents. This study aimed to identify the novel regulators from "hidden" proteins encoded by long noncoding RNAs (lncRNAs) involved in tumor metastasis and chemoresistance. Methods: CRISPR/Cas9 library functional screening was employed to identify the critical suppressor of cancer metastasis in highly invasive CRC models. Western blotting, immunofluorescence staining, invasion, migration, wound healing, WST-1, colony formation, gain- and loss-of-function experiments, in vivo experimental metastasis models, multiplex immunohistochemical staining, immunohistochemistry, qRT-PCR, and RT-PCR were used to assess the functional and clinical significance of FOXP3, PRDM16-DT, HNRNPA2B1, and L-CHEK2. RNA-sequencing, co-immunoprecipitation, qRT-PCR, RT-PCR, RNA affinity purification, RNA immunoprecipitation, MeRIP-quantitative PCR, fluorescence in situ hybridization, chromatin immunoprecipitation and luciferase reporter assay were performed to gain mechanistic insights into the role of PRDM16-DT in cancer metastasis and chemoresistance. An oxaliplatin-resistant CRC cell line was established by in vivo selection. WST-1, colony formation, invasion, migration, Biacore technology, gain- and loss-of-function experiments and an in vivo experimental metastasis model were used to determine the function and mechanism of cimicifugoside H-1 in CRC. Results: The novel protein PRDM16-DT, encoded by LINC00982, was identified as a cancer metastasis and chemoresistance suppressor. The down-regulated level of PRDM16-DT was positively associated with malignant phenotypes and poor prognosis of CRC patients. Transcriptionally regulated by FOXP3, PRDM16-DT directly interacted with HNRNPA2B1 and competitively decreased HNRNPA2B1 binding to exon 9 of CHEK2, resulting in the formation of long CHEK2 (L-CHEK2), subsequently promoting E-cadherin secretion. PRDM16-DT-induced E-cadherin secretion inhibited fibroblast activation, which in turn suppressed CRC metastasis by decreasing MMP9 secretion. Cimicifugoside H-1, a natural compound, can bind to LEU89, HIS91, and LEU92 of FOXP3 and significantly upregulated PRDM16-DT expression to repress CRC metastasis and reverse oxaliplatin resistance. Conclusions: lncRNA LINC00982 can express a new protein PRDM16-DT to function as a novel regulator in cancer metastasis and drug resistance of CRC. Cimicifugoside H-1 can act on the upstream of the PRDM16-DT signaling pathway to alleviate cancer chemoresistance.

Exon-junction complex association with stalled ribosomes and slow translation-independent disassembly.

Exon junction complexes are deposited at exon-exon junctions during splicing. They are primarily known to activate non-sense mediated degradation of transcripts harbouring premature stop codons before the last intron. According to a popular model, exon-junction complexes accompany mRNAs to the cytoplasm where the first translating ribosome pushes them out. However, they are also removed by uncharacterized, translation-independent mechanisms. Little is known about kinetic and transcript specificity of these processes. Here we tag core subunits of exon-junction complexes with complementary split nanoluciferase fragments to obtain sensitive and quantitative assays for complex formation. Unexpectedly, exon-junction complexes form large stable mRNPs containing stalled ribosomes. Complex assembly and disassembly rates are determined after an arrest in transcription and/or translation. 85% of newly deposited exon-junction complexes are disassembled by a translation-dependent mechanism. However as this process is much faster than the translation-independent one, only 30% of the exon-junction complexes present in cells at steady state require translation for disassembly. Deep RNA sequencing shows a bias of exon-junction complex bound transcripts towards microtubule and centrosome coding ones and demonstrate that the lifetimes of exon-junction complexes are transcript-specific. This study provides a dynamic vision of exon-junction complexes and uncovers their unexpected stable association with ribosomes.

An interpretable model of pre-mRNA splicing for animal and plant genes.

Pre-mRNA splicing is a fundamental step in gene expression, conserved across eukaryotes, in which the spliceosome recognizes motifs at the 3' and 5' splice sites (SSs), excises introns, and ligates exons. SS recognition and pairing is often influenced by protein splicing factors (SFs) that bind to splicing regulatory elements (SREs). Here, we describe SMsplice, a fully interpretable model of pre-mRNA splicing that combines models of core SS motifs, SREs, and exonic and intronic length preferences. We learn models that predict SS locations with 83 to 86% accuracy in fish, insects, and plants and about 70% in mammals. Learned SRE motifs include both known SF binding motifs and unfamiliar motifs, and both motif classes are supported by genetic analyses. Our comparisons across species highlight similarities between non-mammals, increased reliance on intronic SREs in plant splicing, and a greater reliance on SREs in mammalian splicing.

Analysis of 3760 hematologic malignancies reveals rare transcriptomic aberrations of driver genes.

BACKGROUND: Rare oncogenic driver events, particularly affecting the expression or splicing of driver genes, are suspected to substantially contribute to the large heterogeneity of hematologic malignancies. However, their identification remains challenging. METHODS: To address this issue, we generated the largest dataset to date of matched whole genome sequencing and total RNA sequencing of hematologic malignancies from 3760 patients spanning 24 disease entities. Taking advantage of our dataset size, we focused on discovering rare regulatory aberrations. Therefore, we called expression and splicing outliers using an extension of the workflow DROP (Detection of RNA Outliers Pipeline) and AbSplice, a variant effect predictor that identifies genetic variants causing aberrant splicing. We next trained a machine learning model integrating these results to prioritize new candidate disease-specific driver genes. RESULTS: We found a median of seven expression outlier genes, two splicing outlier genes, and two rare splice-affecting variants per sample. Each category showed significant enrichment for already well-characterized driver genes, with odds ratios exceeding three among genes called in more than five samples. On held-out data, our integrative modeling significantly outperformed modeling based solely on genomic data and revealed promising novel candidate driver genes. Remarkably, we found a truncated form of the low density lipoprotein receptor LRP1B transcript to be aberrantly overexpressed in about half of hairy cell leukemia variant (HCL-V) samples and, to a lesser extent, in closely related B-cell neoplasms. This observation, which was confirmed in an independent cohort, suggests LRP1B as a novel marker for a HCL-V subclass and a yet unreported functional role of LRP1B within these rare entities. CONCLUSIONS: Altogether, our census of expression and splicing outliers for 24 hematologic malignancy entities and the companion computational workflow constitute unique resources to deepen our understanding of rare oncogenic events in hematologic cancers.

Identification and characterization of whole blood gene expression and splicing quantitative trait loci during early to mid-lactation of dairy cattle.

BACKGROUND: Characterization of regulatory variants (e.g., gene expression quantitative trait loci, eQTL; gene splicing QTL, sQTL) is crucial for biologically interpreting molecular mechanisms underlying loci associated with complex traits. However, regulatory variants in dairy cattle, particularly in specific biological contexts (e.g., distinct lactation stages), remain largely unknown. In this study, we explored regulatory variants in whole blood samples collected during early to mid-lactation (22-150 days after calving) of 101 Holstein cows and analyzed them to decipher the regulatory mechanisms underlying complex traits in dairy cattle. RESULTS: We identified 14,303 genes and 227,705 intron clusters expressed in the white blood cells of 101 cattle. The average heritability of gene expression and intron excision ratio explained by cis-SNPs is 0.28 ± 0.13 and 0.25 ± 0.13, respectively. We identified 23,485 SNP-gene expression pairs and 18,166 SNP-intron cluster pairs in dairy cattle during early to mid-lactation. Compared with the 2,380,457 cis-eQTLs reported to be present in blood in the Cattle Genotype-Tissue Expression atlas (CattleGTEx), only 6,114 cis-eQTLs (P < 0.05) were detected in the present study. By conducting colocalization analysis between cis-e/sQTL and the results of genome-wide association studies (GWAS) from four traits, we identified a cis-e/sQTL (rs109421300) of the DGAT1 gene that might be a key marker in early to mid-lactation for milk yield, fat yield, protein yield, and somatic cell score (PP4 > 0.6). Finally, transcriptome-wide association studies (TWAS) revealed certain genes (e.g., FAM83H and TBC1D17) whose expression in white blood cells was significantly (P < 0.05) associated with complex traits. CONCLUSIONS: This study investigated the genetic regulation of gene expression and alternative splicing in dairy cows during early to mid-lactation and provided new insights into the regulatory mechanisms underlying complex traits of economic importance.

Identification of two novel heterozygous variants of SMC3 with Cornelia de Lange syndrome.

BACKGROUND: Cornelia de Lange syndrome (CdLS) is a multisystem genetic disorder, and cases caused by variants in the structural maintenance of chromosomes protein 3 (SMC3) gene are uncommon. Here, we report two cases of CdLS associated with novel pathogenic variants in SMC3 from two Chinese families. METHODS: Clinical presentations of two patients with CdLS were evaluated, and specimens from the patients and other family members were collected for Trio-based whole-exome sequencing. Pyrosequencing, chip-based digital PCR, minigene splicing assay, and in silico analysis were carried out to elucidate the impact of novel variants. RESULTS: Novel heterozygous variants in SMC3 were identified in each proband. One harbored a novel splicing and mosaic variant (c.2535+1G>A) in SMC3. The mutated allele G>A conversion was approximately 23.1% by digital PCR, which indicated that 46.2% of peripheral blood cells had this variant. Additionally, in vitro minigene splicing analysis validated that the c.2535+1G>A variant led to an exon skipping in messenger RNA splicing. The other carried a heterozygous variant (c.435C>A), which was predicted to be pathogenic as well as significantly altered in local electrical potential. The former showed multiple abnormalities and marked clinical severity, and the latter mainly exhibited a speech developmental disorder and slightly facial anomalies. CONCLUSION: Both patients were clinically diagnosed with Cornelia de Lange syndrome 3 (CdLS3). The newly identified SMC3 gene variants can expand the understanding of CdLS3 and provide reliable evidence for genetic counseling to the affected family.

Accelerating spliced alignment of long RNA sequencing reads using parallel maximal exact match retrieval.

The genomics landscape has undergone a revolutionary transformation with the emergence of third-generation sequencing technologies. Fueled by the exponential surge in sequencing data, there is an urgent demand for accurate and rapid algorithms to effectively handle this burgeoning influx. Under such circumstances, we developed a parallelized, yet accuracy-lossless algorithm for maximal exact match (MEM) retrieval to strategically address the computational bottleneck of uLTRA, a leading spliced alignment algorithm known for its precision in handling long RNA sequencing (RNA-seq) reads. The design of the algorithm incorporates a multi-threaded strategy, enabling the concurrent processing of multiple reads simultaneously. Additionally, we implemented the serialization of index required for MEM retrieval to facilitate its reuse, resulting in accelerated startup for practical tasks. Extensive experiments demonstrate that our parallel algorithm achieves significant improvements in runtime, speedup, throughput, and memory usage. When applied to the largest human dataset, the algorithm achieves an impressive speedup of 10.78 × , significantly improving throughput on a large scale. Moreover, the integration of the parallel MEM retrieval algorithm into the uLTRA pipeline introduces a dual-layered parallel capability, consistently yielding a speedup of 4.99 × compared to the multi-process and single-threaded execution of uLTRA. The thorough analysis of experimental results underscores the adept utilization of parallel processing capabilities and its advantageous performance in handling large datasets. This study provides a showcase of parallelized strategies for MEM retrieval within the context of spliced alignment algorithm, effectively facilitating the process of RNA-seq data analysis. The code is available at https://github.com/RongxingWong/AcceleratingSplicedAlignment.

Intramolecular autoinhibition regulates the selectivity of PRPF40A tandem WW domains for proline-rich motifs.

PRPF40A plays an important role in the regulation of pre-mRNA splicing by mediating protein-protein interactions in the early steps of spliceosome assembly. By binding to proteins at the 5´ and 3´ splice sites, PRPF40A promotes spliceosome assembly by bridging the recognition of the splices. The PRPF40A WW domains are expected to recognize proline-rich sequences in SF1 and SF3A1 in the early spliceosome complexes E and A, respectively. Here, we combine NMR, SAXS and ITC to determine the structure of the PRPF40A tandem WW domains in solution and characterize the binding specificity and mechanism for proline-rich motifs recognition. Our structure of the PRPF40A WW tandem in complex with a high-affinity SF1 peptide reveals contributions of both WW domains, which also enables tryptophan sandwiching by two proline residues in the ligand. Unexpectedly, a proline-rich motif in the N-terminal region of PRPF40A mediates intramolecular interactions with the WW tandem. Using NMR, ITC, mutational analysis in vitro, and immunoprecipitation experiments in cells, we show that the intramolecular interaction acts as an autoinhibitory filter for proof-reading of high-affinity proline-rich motifs in bona fide PRPF40A binding partners. We propose that similar autoinhibitory mechanisms are present in most WW tandem-containing proteins to enhance binding selectivity and regulation of WW/proline-rich peptide interaction networks.

Genome organization around nuclear speckles drives mRNA splicing efficiency.

The nucleus is highly organized, such that factors involved in the transcription and processing of distinct classes of RNA are confined within specific nuclear bodies1,2. One example is the nuclear speckle, which is defined by high concentrations of protein and noncoding RNA regulators of pre-mRNA splicing3. What functional role, if any, speckles might play in the process of mRNA splicing is unclear4,5. Here we show that genes localized near nuclear speckles display higher spliceosome concentrations, increased spliceosome binding to their pre-mRNAs and higher co-transcriptional splicing levels than genes that are located farther from nuclear speckles. Gene organization around nuclear speckles is dynamic between cell types, and changes in speckle proximity lead to differences in splicing efficiency. Finally, directed recruitment of a pre-mRNA to nuclear speckles is sufficient to increase mRNA splicing levels. Together, our results integrate the long-standing observations of nuclear speckles with the biochemistry of mRNA splicing and demonstrate a crucial role for dynamic three-dimensional spatial organization of genomic DNA in driving spliceosome concentrations and controlling the efficiency of mRNA splicing.

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.

The debranching enzyme Dbr1 regulates lariat turnover and intron splicing.

The majority of genic transcription is intronic. Introns are removed by splicing as branched lariat RNAs which require rapid recycling. The branch site is recognized during splicing catalysis and later debranched by Dbr1 in the rate-limiting step of lariat turnover. Through generation of a viable DBR1 knockout cell line, we find the predominantly nuclear Dbr1 enzyme to encode the sole debranching activity in human cells. Dbr1 preferentially debranches substrates that contain canonical U2 binding motifs, suggesting that branchsites discovered through sequencing do not necessarily represent those favored by the spliceosome. We find that Dbr1 also exhibits specificity for particular 5' splice site sequences. We identify Dbr1 interactors through co-immunoprecipitation mass spectrometry. We present a mechanistic model for Dbr1 recruitment to the branchpoint through the intron-binding protein AQR. In addition to a 20-fold increase in lariats, Dbr1 depletion increases exon skipping. Using ADAR fusions to timestamp lariats, we demonstrate a defect in spliceosome recycling. In the absence of Dbr1, spliceosomal components remain associated with the lariat for a longer period of time. As splicing is co-transcriptional, slower recycling increases the likelihood that downstream exons will be available for exon skipping.

DNMT3B splicing dysregulation mediated by SMCHD1 loss contributes to DUX4 overexpression and FSHD pathogenesis.

Structural maintenance of chromosomes flexible hinge domain-containing 1 (SMCHD1) is a noncanonical SMC protein and an epigenetic regulator. Mutations in SMCHD1 cause facioscapulohumeral muscular dystrophy (FSHD), by overexpressing DUX4 in muscle cells. Here, we demonstrate that SMCHD1 is a key regulator of alternative splicing in various cell types. We show how SMCHD1 loss causes splicing alterations of DNMT3B, which can lead to hypomethylation and DUX4 overexpression. Analyzing RNA sequencing data from muscle biopsies of patients with FSHD and Smchd1 knocked out cells, we found mis-splicing of hundreds of genes upon SMCHD1 loss. We conducted a high-throughput screen of splicing factors, revealing the involvement of the splicing factor RBM5 in the mis-splicing of DNMT3B. Subsequent RNA immunoprecipitation experiments confirmed that SMCHD1 is required for RBM5 recruitment. Last, we show that mis-splicing of DNMT3B leads to hypomethylation of the D4Z4 region and to DUX4 overexpression. These results suggest that DNMT3B mis-splicing due to SMCHD1 loss plays a major role in FSHD pathogenesis.

Assessing the Impact of Novel BRCA1 Exon 11 Variants on Pre-mRNA Splicing.

Our study focused on assessing the effects of three newly identified BRCA1 exon 11 variants (c.1019T>C, c.2363T>G, and c.3192T>C) on breast cancer susceptibility. Using computational predictions and experimental splicing assays, we evaluated their potential as pathogenic mutations. Our in silico analyses suggested that the c.2363T>G and c.3192T>C variants could impact both splicing and protein function, resulting in the V340A and V788G mutations, respectively. We further examined their splicing effects using minigene assays in MCF7 and SKBR3 breast cancer cell lines. Interestingly, we found that the c.2363T>G variant significantly altered splicing patterns in MCF7 cells but not in SKBR3 cells. This finding suggests a potential influence of cellular context on the variant's effects. While attempts to correlate in silico predictions with RNA binding factors were inconclusive, this observation underscores the complexity of splicing regulation. Splicing is governed by various factors, including cellular contexts and protein interactions, making it challenging to predict outcomes accurately. Further research is needed to fully understand the functional consequences of the c.2363T>G variant in breast cancer pathogenesis. Integrating computational predictions with experimental data will provide valuable insights into the role of alternative splicing regulation in different breast cancer types and stages.

MAPP unravels frequent co-regulation of splicing and polyadenylation by RNA-binding proteins and their dysregulation in cancer.

Maturation of eukaryotic pre-mRNAs via splicing and polyadenylation is modulated across cell types and conditions by a variety of RNA-binding proteins (RBPs). Although there exist over 1,500 RBPs in human cells, their binding motifs and functions still remain to be elucidated, especially in the complex environment of tissues and in the context of diseases. To overcome the lack of methods for the systematic and automated detection of sequence motif-guided pre-mRNA processing regulation from RNA sequencing (RNA-Seq) data we have developed MAPP (Motif Activity on Pre-mRNA Processing). Applying MAPP to RBP knock-down experiments reveals that many RBPs regulate both splicing and polyadenylation of nascent transcripts by acting on similar sequence motifs. MAPP not only infers these sequence motifs, but also unravels the position-dependent impact of the RBPs on pre-mRNA processing. Interestingly, all investigated RBPs that act on both splicing and 3' end processing exhibit a consistently repressive or activating effect on both processes, providing a first glimpse on the underlying mechanism. Applying MAPP to normal and malignant brain tissue samples unveils that the motifs bound by the PTBP1 and RBFOX RBPs coordinately drive the oncogenic splicing program active in glioblastomas demonstrating that MAPP paves the way for characterizing pre-mRNA processing regulators under physiological and pathological conditions.

Comprehensive profiling of cancer neoantigens from aberrant RNA splicing.

BACKGROUND: Cancer neoantigens arise from protein-altering somatic mutations in tumor and rank among the most promising next-generation immuno-oncology agents when used in combination with immune checkpoint inhibitors. We previously developed a computational framework, REAL-neo, for identification, quality control, and prioritization of both class-I and class-II human leucocyte antigen (HLA)-presented neoantigens resulting from somatic single-nucleotide mutations, small insertions and deletions, and gene fusions. In this study, we developed a new module, SPLICE-neo, to identify neoantigens from aberrant RNA transcripts from two distinct sources: (1) DNA mutations within splice sites and (2) de novo RNA aberrant splicings. METHODS: First, SPLICE-neo was used to profile all DNA splice-site mutations in 11,892 tumors from The Cancer Genome Atlas (TCGA) and identified 11 profiles of splicing donor or acceptor site gains or losses. Transcript isoforms resulting from the top seven most frequent profiles were computed using novel logic models. Second, SPLICE-neo identified de novo RNA splicing events using RNA sequencing reads mapped to novel exon junctions from either single, double, or multiple exon-skipping events. The aberrant transcripts from both sources were then ranked based on isoform expression levels and z-scores assuming that individual aberrant splicing events are rare. Finally, top-ranked novel isoforms were translated into protein, and the resulting neoepitopes were evaluated for neoantigen potential using REAL-neo. The top splicing neoantigen candidates binding to HLA-A*02:01 were validated using in vitro T2 binding assays. RESULTS: We identified abundant splicing neoantigens in four representative TCGA cancers: BRCA, LUAD, LUSC, and LIHC. In addition to their substantial contribution to neoantigen load, several splicing neoantigens were potent tumor antigens with stronger bindings to HLA compared with the positive control of antigens from influenza virus. CONCLUSIONS: SPLICE-neo is the first tool to comprehensively identify and prioritize splicing neoantigens from both DNA splice-site mutations and de novo RNA aberrant splicings. There are two major advances of SPLICE-neo. First, we developed novel logic models that assemble and prioritize full-length aberrant transcripts from DNA splice-site mutations. Second, SPLICE-neo can identify exon-skipping events involving more than two exons, which account for a quarter to one-third of all skipping events.

Copper-binding proteins and exonic splicing enhancers and silencers.

Eukaryotic DNA codes not only for proteins but contains a wealth of information required for accurate splicing of messenger RNA precursors and inclusion of constitutively or alternatively spliced exons in mature transcripts. This "auxiliary" splicing code has been characterized as exonic splicing enhancers and silencers (ESE and ESS). The exact interplay between protein and splicing codes is, however, poorly understood. Here, we show that exons encoding copper-coordinating amino acids in human cuproproteins lack ESEs and/or have an excess of ESSs, yet RNA sequencing and expressed sequence tags data show that they are more efficiently included in mature transcripts by the splicing machinery than average exons. Their largely constitutive inclusion in messenger RNA is facilitated by stronger splice sites, including polypyrimidine tracts, consistent with an important role of the surrounding intron architecture in ensuring high expression of metal-binding residues during evolution. ESE/ESS profiles of codons and entire exons that code for copper-coordinating residues were very similar to those encoding residues that coordinate zinc but markedly different from those that coordinate calcium. Together, these results reveal how the traditional and auxiliary splicing motifs responded to constraints of metal coordination in proteins.