Comprehensive characterisation of intronic mis-splicing mutations in human cancers

Alternative Splicing Regulation in Metabolic Disorders

Alternative splicing (AS) is a fundamental mechanism for enhancing transcriptome diversity and regulating gene expression, crucial for various cellular processes and the development of complex traits. This review examines the role of AS in metabolic disorders, including obesity, weight loss, dyslipidemias, and metabolic syndrome. We explore the molecular mechanisms underlying AS regulation, focusing on the interplay between cis-acting elements and trans-acting factors, and the influence of RNA-binding proteins (RBPs). Advances in high-throughput sequencing and bioinformatics have unveiled the extensive landscape of AS events across different tissues and conditions, highlighting the importance of tissue-specific splicing in metabolic regulation. We discuss the impact of genetic variants on AS, with a particular emphasis on splicing quantitative trait loci (sQTLs) and their association with cardiometabolic traits. The review also covers the regulation of spliceosome components by phosphorylation, the role of m6A modification in AS, and the interaction between transcription and splicing. Additionally, we address the clinical relevance of AS, illustrating how splicing misregulation contributes to metabolic diseases and the potential for therapeutic interventions targeting splicing mechanisms. This comprehensive overview underscores the significance of AS in metabolic health and disease, advocating for further research to harness its therapeutic potential.

Clinicopathologic and molecular characteristics of neuroendocrine carcinomas of the gallbladder

Gallbladder neuroendocrine carcinomas (GB-NECs) are a rare subtype of malignant gallbladder cancer (GBC). The genetic and molecular characteristics of GB-NECs are rarely reported. This study aims to assess the frequency of microsatellite instability (MSI) in GB-NECs and characterize their clinicopathologic and molecular features in comparison with gallbladder adenocarcinomas (GB-ADCs). Data from six patients with primary GB-NECs and 13 with GB-ADCs were collected and reevaluated. MSI assay, immunohistochemistry for mismatch repair proteins (MLH1, MSH2, MSH6, and PMS2), comprehensive genomic profiling (CGP) via next-generation sequencing (NGS), and evaluation of tumor mutation burden (TMB) were conducted on these samples. The six GB-NEC cases were all female, with a mean age of 62.0±9.2 years. Of these, two cases were diagnosed as large cell neuroendocrine carcinomas (LCNECs), while the remaining four were small cell neuroendocrine carcinomas (SCNECs). Microsatellite states observed in both GB-NECs and GB-ADCs were consistently microsatellite stable (MSS). Notably, TP53 (100%, 6/6) and RB1 (100%, 6/6) exhibited the highest mutation frequency in GB-NECs, followed by SMAD4 (50%, 3/6), GNAS (50%, 3/6), and RICTOR (33%, 2/6), with RB1, GNAS, and RICTOR specifically present in GB-NECs. Immunohistochemical (IHC) assays of p53 and Rb in the six GB-NECs were highly consistent with genetic mutations detected by targeted NGS. Moreover, no statistical difference was observed in TMB between GB-NECs and GB-ADCs (p=0.864). Although overall survival in GB-NEC patients tended to be worse than in GB-ADC patients, this difference did not reach statistical significance (p=0.119). This study has identified the microsatellite states and molecular mutation features of GB-NECs, suggesting that co-mutations in TP53 and RB1 may signify a neuroendocrine inclination in GB-NECs. The IHC assay provides an effective complement to targeted NGS for determining the functional status of p53 and Rb in clinical practice.

iSoMAs: Finding isoform expression and somatic mutation associations in human cancers

Aberrant alternative splicing, prevalent in cancer, impacts various cancer hallmarks involving proliferation, angiogenesis, and invasion. Splicing disruption often results from somatic point mutations rewiring functional pathways to support cancer cell survival. We introduce iSoMAs (iSoform expression and somatic Mutation Association), an efficient computational pipeline leveraging principal component analysis technique, to explore how somatic mutations influence transcriptome-wide gene expression at the isoform level. Applying iSoMAs to 33 cancer types comprising 9,738 tumor samples in The Cancer Genome Atlas, we identified 908 somatically mutated genes significantly associated with altered isoform expression across three or more cancer types. Mutations linked to differential isoform expression occurred through both cis- and trans-acting mechanisms, involving well-known oncogenes/suppressor genes, RNA binding protein and splicing factor genes. With wet-lab experiments, we verified direct association between TP53 mutations and differential isoform expression in cell cycle genes. Additional iSoMAs genes have been validated in the literature with independent cohorts and/or methods. Despite the complexity of cancer, iSoMAs attains computational efficiency via dimension reduction strategy and reveals critical associations between regulatory factors and transcriptional landscapes.

Activation of Wnt/β-catenin signaling is critical for the tumorigenesis of choroid plexus

BACKGROUND

The choroid plexus (ChP) is the secretory epithelial structure located in the brain ventricles. Choroid plexus tumors (CPTs) are rare neoplasms predominantly occurring in young patients with intensified malignancy in children. CPT treatment is hindered by insufficient knowledge of tumor pathology and the limited availability of valid models.

METHODS

Genomic and transcriptomic data from CPT patients were analyzed to identify the putative pathological pathway. Cellular and molecular techniques were employed to validate bioinformatic results in CPT patient samples. Pharmacologic inhibition of Wnt/β-catenin signaling was assessed in CPT cells. Cell-based assays of ChP cell lines were performed following CRISPR-Cas9-derived knockout and overexpression of Wnt/β-catenin pathway genes. A 3D CPT model was generated through CRISPR-Cas9-derived knockout of APC.

RESULTS

We discovered that Wnt/β-catenin signaling is activated in human CPTs, likely as a consequence of large-scale chromosomal instability events of the CPT genomes. We demonstrated that CPT-derived cells depend on autocrine Wnt/β-catenin signaling for survival. Constitutive Wnt/β-catenin pathway activation, either through knockout of the negative regulator APC or overexpression of the ligand WNT3A, induced tumorigenic properties in ChP 2D in vitro models. Increased activation of the Wnt/β-catenin pathway in ChP organoids, through treatment with a potent GSK3β inhibitor, reduced the differentiation of mature ChP epithelial cells. Remarkably, the depletion of APC was sufficient to induce the oncogenic transformation of ChP organoids.

CONCLUSIONS

Our research identifies Wnt/β-catenin signaling as a critical driver of CPT tumorigenesis and provides the first 3D in vitro model for future pathological and therapeutic studies of CPT.

A novel RUNX2 splice site mutation in Chinese associated with cleidocranial dysplasia

Pathogenic genes in most patients with cleidocranial dysplasia have been confirmed to be runt-related transcription factor 2 (RUNX2), which controls mutations in specific osteoblast transcription factors and affects skull ossification and suture adhesion. This study aimed to explore the role of RUNX2 mutations. Here, we report a rare case of a splice site mutation in a Chinese population with typical cleidocranial dysplasia symptoms, cranial suture insufficiency, clavicle dysplasia, and dental anomalies. Peripheral blood samples from the proband and her mother were subjected to Sanger sequencing. The expression levels of RUNX2 before and after mutation were verified using digital PCR (dPCR). The results revealed a classic mutation at the fifth base of the intron 5 initiation splicing sequence (NM001024630.4: C.685+5G > A). The mutation rate in the proband was 53 %, while the mother did not have any mutations. The secondary RNA structure of the RUNX2 gene in the progenitor was predicted to change, and the structural free energy was low in the wild-type, with the stem folded first and the structure being relatively stable. After the mutation, the free energy increased. This finding enriches the RUNX2 mutation library of CD-related genes in Chinese individuals.

SPLICER: a highly efficient base editing toolbox that enables in vivo therapeutic exon skipping

Exon skipping technologies enable exclusion of targeted exons from mature mRNA transcripts, which have broad applications in medicine and biotechnology. Existing techniques including antisense oligonucleotides, targetable nucleases, and base editors, while effective for specific applications, remain hindered by transient effects, genotoxicity, and inconsistent exon skipping. To overcome these limitations, here we develop SPLICER, a toolbox of next-generation base editors containing near-PAMless Cas9 nickase variants fused to adenosine or cytosine deaminases for the simultaneous editing of splice acceptor (SA) and splice donor (SD) sequences. Synchronized SA and SD editing improves exon skipping, reduces aberrant splicing, and enables skipping of exons refractory to single splice site editing. To demonstrate the therapeutic potential of SPLICER, we target APP exon 17, which encodes amino acids that are cleaved to form Aβ plaques in Alzheimer's disease. SPLICER reduces the formation of Aβ42 peptides in vitro and enables efficient exon skipping in a mouse model of Alzheimer's disease. Overall, SPLICER is a widely applicable and efficient exon skipping toolbox.

DOCK8 deficiency due to a deep intronic variant in two kindreds with hyper-IgE syndrome

Dedicator of cytokinesis 8 (DOCK8) deficiency underlies the majority of cases of patients with autosomal recessive form of the hyper-immunoglobulin E syndrome (HIES). Most DOCK8 mutations involve deletions and splice junction mutations that abrogate protein expression. However, a few patients whose presentation is reminiscent of DOCK8 deficiency have no identifiable mutations. Using Whole Exome Sequencing (WES), we identified a deep intronic homozygous DOCK8 variant located in intron 36 (c.4626 + 76 A > G) in two unrelated patients with features of HIES that resulted in an in-frame 75 base pair intronic sequence insertion in DOCK8 cDNA, resulting in a premature stop codon (p.S1542ins6Ter). This variant resulted in variable decrease in DOCK8 expression that was associated with impaired T cell receptor-triggered actin polymerization, decreased IL-6-induced STAT3 phosphorylation, reduced expression of the Th17 cell markers CCR6 and IL-17, and higher frequencies of GATA3+ T cells indicative of Th2 skewing. Our approach extends the reach of WES in identifying disease-related intronic variants. It highlights the role of non-coding mutations in immunodeficiency disorders, including DOCK8 deficiency, and emphasizes the need to explore these mutations in unexplained inborn errors of immunity.

Autophagy Gene BECN1 Intronic Variant rs9890617 Predisposes Individuals to Hepatitis B Virus Infection

Beclin 1 protein encoded by the BECN1 gene plays a critical role in the autophagy pathway which is utilized by the Hepatitis B virus (HBV) for its replication. HBV is known for the subversion of the host's autophagy process for its multiplication. The aim of this study was to determine the role of BECN1 intronic variants in HBV susceptibility. Intronic region variant rs9890617 was analyzed using Human splicing finder v3.1 and was found to alter splicing signals. A total of 712 individuals (494 HBV infected and 218 healthy controls) were recruited in the study and genotyped by applying Polymerase Chain Reaction-Restriction Fragment Length Polymorphism (PCR-RFLP). Statistical analysis revealed that the mutant allele T of rs9890617 was significantly associated with the overall disease risk in the allelic model (OR 1.41; 95%CI 1.00-1.99, p = 0.04). On stratifying the data based on the different stages of HBV infection, the mutant genotype showed a significant association with the chronic group in allelic (OR 1.62; 95%CI 1.11-2.39, p = 0.01), dominant (OR 1.64; 95%CI 1.07-2.52, p = 0.02), and co-dominant (OR 1.55; 95%CI 1.00-2.40, p = 0.04) models. Overall, this is the first study regarding beclin 1 variant rs9890617 and we found a significant association of the mutant T allele with the genetic predisposition to HBV infection.

Systems modeling of oncogenic G-protein and GPCR signaling reveals unexpected differences in downstream pathway activation

Mathematical models of biochemical reaction networks are an important and emerging tool for the study of cell signaling networks involved in disease processes. One promising potential application of such mathematical models is the study of how disease-causing mutations promote the signaling phenotype that contributes to the disease. It is commonly assumed that one must have a thorough characterization of the network readily available for mathematical modeling to be useful, but we hypothesized that mathematical modeling could be useful when there is incomplete knowledge and that it could be a tool for discovery that opens new areas for further exploration. In the present study, we first develop a mechanistic mathematical model of a G-protein coupled receptor signaling network that is mutated in almost all cases of uveal melanoma and use model-driven explorations to uncover and explore multiple new areas for investigating this disease. Modeling the two major, mutually-exclusive, oncogenic mutations (Gαq/11 and CysLT2R) revealed the potential for previously unknown qualitative differences between seemingly interchangeable disease-promoting mutations, and our experiments confirmed oncogenic CysLT2R was impaired at activating the FAK/YAP/TAZ pathway relative to Gαq/11. This led us to hypothesize that CYSLTR2 mutations in UM must co-occur with other mutations to activate FAK/YAP/TAZ signaling, and our bioinformatic analysis uncovers a role for co-occurring mutations involving the plexin/semaphorin pathway, which has been shown capable of activating this pathway. Overall, this work highlights the power of mechanism-based computational systems biology as a discovery tool that can leverage available information to open new research areas.

The Neo-Open Reading Frame Peptides That Comprise the Tumor Framome Are a Rich Source of Neoantigens for Cancer Immunotherapy

Identification of immunogenic cancer neoantigens as targets for therapy is challenging. Here, we integrate the whole-genome and long-read transcript sequencing of cancers to identify the collection of neo-open reading frame peptides (NOP) expressed in tumors. We termed this collection of NOPs the tumor framome. NOPs represent tumor-specific peptides that are different from wild-type proteins and may be strongly immunogenic. We describe a class of hidden NOPs that derive from structural genomic variants involving an upstream protein coding gene driving expression and translation of noncoding regions of the genome downstream of a rearrangement breakpoint, i.e., where no gene annotation or evidence for transcription exists. The entire collection of NOPs represents a vast number of possible neoantigens particularly in tumors with many structural genomic variants and a low number of missense mutations. We show that NOPs are immunogenic and epitopes derived from NOPs can bind to MHC class I molecules. Finally, we provide evidence for the presence of memory T cells specific for hidden NOPs in peripheral blood from a patient with lung cancer. This work highlights NOPs as a major source of possible neoantigens for personalized cancer immunotherapy and provides a rationale for analyzing the complete cancer genome and transcriptome as a basis for the detection of NOPs.

A single-cell strategy for the identification of intronic variants related to mis-splicing in pancreatic cancer

Most clinical diagnostic and genomic research setups focus almost exclusively on coding regions and essential splice sites, thereby overlooking other non-coding variants. As a result, intronic variants that can promote mis-splicing events across a range of diseases, including cancer, are yet to be systematically investigated. Such investigations would require both genomic and transcriptomic data, but there currently exist very few datasets that satisfy these requirements. We address this by developing a single-nucleus full-length RNA-sequencing approach that allows for the detection of potentially pathogenic intronic variants. We exemplify the potency of our approach by applying pancreatic cancer tumor and tumor-derived specimens and linking intronic variants to splicing dysregulation. We specifically find that prominent intron retention and pseudo-exon activation events are shared by the tumors and affect genes encoding key transcriptional regulators. Our work paves the way for the assessment and exploitation of intronic mutations as powerful prognostic markers and potential therapeutic targets in cancer.

Characterization and Engineered U1 snRNA Rescue of Splicing Variants in a Turkish Neurodevelopmental Disease Cohort

Although they are rare in the population, rare neurodevelopmental disorders (RNDDs) constitute a significant portion of all rare diseases. While advancements in sequencing technologies led to improvements in diagnosing and managing rare neurodevelopmental diseases, accurate pathogenicity classification of the identified variants is still challenging. Sequence variants altering pre-mRNA splicing make up a significant part of pathogenic variants. Despite advances in the in silico prediction tools, noncanonical splice site variants are one of the groups of variants that pose a challenge in their clinical interpretation. In this study, we analyzed the effects of seven splicing variants we had previously proposed as disease-causing and demonstrated that all but one of the seven variants had a strong or moderate effect on splicing, as assessed by a minigene assay. Next, applying U1 snRNAs engineered for different splicing variants in the corresponding genes and expressed with minigene plasmids in HeLa cells provided a partial correction in four of the studied genes to varying degrees. Findings from our study highlight the importance of in vitro minigene-based assays for the reclassification of putative splice-altering variants of uncertain significance and the therapeutic potential of modified U1 snRNAs in RNDDs.

SPLICER: A Highly Efficient Base Editing Toolbox That Enables In Vivo Therapeutic Exon Skipping

Exon skipping technologies enable exclusion of targeted exons from mature mRNA transcripts, which has broad applications in molecular biology, medicine, and biotechnology. Existing exon skipping techniques include antisense oligonucleotides, targetable nucleases, and base editors, which, while effective for specific applications at some target exons, remain hindered by shortcomings, including transient effects for oligonucleotides, genotoxicity for nucleases and inconsistent exon skipping for base editors. To overcome these limitations, we created SPLICER, a toolbox of next-generation base editors consisting of near-PAMless Cas9 nickase variants fused to adenosine or cytosine deaminases for the simultaneous editing of splice acceptor (SA) and splice donor (SD) sequences. Synchronized SA and SD editing with SPLICER improves exon skipping, reduces aberrant outcomes, including cryptic splicing and intron retention, and enables skipping of exons refractory to single splice-site editing. To demonstrate the therapeutic potential of SPLICER, we targeted APP exon 17, which encodes the amino acid residues that are cleaved to form the Aβ plaques in Alzheimer's disease. SPLICER reduced the formation of Aβ42 peptides in vitro and enabled efficient exon skipping in a mouse model of Alzheimer's disease. Overall, SPLICER is a widely applicable and efficient toolbox for exon skipping with broad therapeutic applications.

Functional requirement of alternative splicing in epithelial-mesenchymal transition of pancreatic circulating tumor

Circulating tumor cells (CTCs) that undergo epithelial-to-mesenchymal transition (EMT) can provide valuable information regarding metastasis and potential therapies. However, current studies on the EMT overlook alternative splicing. Here, we used single-cell full-length transcriptome data and mRNA sequencing of CTCs to identify stage-specific alternative splicing of partial EMT and mesenchymal states during pancreatic cancer metastasis. We classified definitive tumor and normal epithelial cells via genetic aberrations and demonstrated dynamic changes in the epithelial-mesenchymal continuum in both epithelial cancer cells and CTCs. We provide the landscape of alternative splicing in CTCs at different stages of EMT, uncovering cell-type-specific splicing patterns and splicing events in cell surface proteins suitable for therapies. We show that MBNL1 governs cell fate through alternative splicing independently of changes in gene expression and affects the splicing pattern during EMT. We found a high frequency of events that contained multiple premature termination codons and were enriched with C and G nucleotides in close proximity, which influence the likelihood of stop codon readthrough and expand the range of potential therapeutic targets. Our study provides insights into the EMT transcriptome's dynamic changes and identifies potential diagnostic and therapeutic targets in pancreatic cancer.

Detecting and understanding meaningful cancerous mutations based on computational models of mRNA splicing

Cancer research has long relied on non-silent mutations. Yet, it has become overwhelmingly clear that silent mutations can affect gene expression and cancer cell fitness. One fundamental mechanism that apparently silent mutations can severely disrupt is alternative splicing. Here we introduce Oncosplice, a tool that scores mutations based on models of proteomes generated using aberrant splicing predictions. Oncosplice leverages a highly accurate neural network that predicts splice sites within arbitrary mRNA sequences, a greedy transcript constructor that considers alternate arrangements of splicing blueprints, and an algorithm that grades the functional divergence between proteins based on evolutionary conservation. By applying this tool to 12M somatic mutations we identify 8K deleterious variants that are significantly depleted within the healthy population; we demonstrate the tool's ability to identify clinically validated pathogenic variants with a positive predictive value of 94%; we show strong enrichment of predicted deleterious mutations across pan-cancer drivers. We also achieve improved patient survival estimation using a proposed set of novel cancer-involved genes. Ultimately, this pipeline enables accelerated insight-gathering of sequence-specific consequences for a class of understudied mutations and provides an efficient way of filtering through massive variant datasets - functionalities with immediate experimental and clinical applications.

Activation-induced cytidine deaminase causes recurrent splicing mutations in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma. A major mutagenic process in DLBCL is aberrant somatic hypermutation (aSHM) by activation-induced cytidine deaminase (AID), which occurs preferentially at RCH/TW sequence motifs proximal to transcription start sites. Splice sequences are highly conserved, rich in RCH/TW motifs, and recurrently mutated in DLBCL. Therefore, we hypothesized that aSHM may cause recurrent splicing mutations in DLBCL. In a meta-cohort of > 1,800 DLBCLs, we found that 77.5% of splicing mutations in 29 recurrently mutated genes followed aSHM patterns. In addition, in whole-genome sequencing (WGS) data from 153 DLBCLs, proximal mutations in splice sequences, especially in donors, were significantly enriched in RCH/TW motifs (p < 0.01). We validated this enrichment in two additional DLBCL cohorts (N > 2,000; p < 0.0001) and confirmed its absence in 12 cancer types without aSHM (N > 6,300). Comparing sequencing data from mouse models with and without AID activity showed that the splice donor sequences were the top genomic feature enriched in AID-induced mutations (p < 0.0001). Finally, we observed that most AID-related splice site mutations are clonal within a sample, indicating that aSHM may cause early loss-of-function events in lymphomagenesis. Overall, these findings support that AID causes an overrepresentation of clonal splicing mutations in DLBCL.

Cryptic splice mutation in the fumarate hydratase gene in patients with clinical manifestations of Hereditary Leiomyomatosis and Renal Cell Cancer

Hereditary leiomyomatosis and renal cell carcinoma (HLRCC) is an autosomal dominant condition characterized by the development of cutaneous and uterine leiomyomas and risk for development of an aggressive form of papillary renal cell cancer. HLRCC is caused by germline inactivating pathogenic variants in the fumarate hydratase (FH) gene, which encodes the enzyme that catalyzes the interconversion of fumarate and L-malate. We utilized enzyme and protein mobility assays to evaluate the FH enzyme in a cohort of patients who showed clinical manifestations of HLRCC but were negative for known pathogenic FH gene variants. FH enzyme activity and protein levels were decreased by 50% or greater in three family members, despite normal FH mRNA expression levels as measured by quantitative PCR. Direct Nanopore RNA sequencing demonstrated 57 base pairs of retained intron sequence between exons 9 and 10 of polyadenylated FH mRNA in these patients, resulting in a truncated FH protein. Genomic sequencing revealed a heterozygous intronic alteration of the FH gene (chr1: 241498239 T/C) resulting in formation of a splice acceptor site near a polypyrimidine tract, and a uterine fibroid obtained from a patient showed loss of heterozygosity at this site. The same intronic FH variant was identified in an unrelated patient who also showed a clinical phenotype of HLRCC. These data demonstrate that careful clinical assessment as well as biochemical characterization of FH enzyme activity, protein expression, direct RNA sequencing, and genomic DNA sequencing of patient-derived cells can identify pathogenic variants outside of the protein coding regions of the FH gene.

Genomic variants exclusively identified in children with birth defects and concurrent malignant tumors predispose to cancer development

Children with birth defects (BD) express distinct clinical features that often have various medical consequences, one of which is predisposition to the development of cancers. Identification of the underlying genetic mechanisms related to the development of cancer in BD patients would allow for preventive measures. We performed a whole genome sequencing (WGS) study on blood-derived DNA samples from 1566 individuals without chromosomal anomalies, including 454 BD probands with at least one type of malignant tumors, 767 cancer-free BD probands, and 345 healthy individuals. Exclusive recurrent variants were identified in BD-cancer and BD-only patients and mapped to their corresponding genomic regions. We observed statistically significant overlaps for protein-coding/ncRNA with exclusive variants in exons, introns, ncRNAs, and 3'UTR regions. Exclusive exonic variants, especially synonymous variants, tend to occur in prior exons locus in BD-cancer children. Intronic variants close to splicing site (< 500 bp from exon) have little overlaps in BD-cancer and BD-only patients. Exonic variants in non-coding RNA (ncRNA) tend to occur in different ncRNAs exons regardless of the overlaps. Notably, genes with 5' UTR variants are almost mutually exclusive between the two phenotypes. In conclusion, we conducted the first genomic study to explore the impact of recurrent variants exclusive to the two distinguished clinical phenotypes under study, BD with or without cancer, demonstrating enrichment of selective protein-coding/ncRNAs differentially expressed between these two phenotypes, suggesting that selective genetic factors may underlie the molecular processes of pediatric cancer development in BD children.

Role of RNA Splicing Mutations in Diffuse Large B Cell Lymphoma

Ribonucleic acid splicing is a crucial process to create a mature mRNA molecule by removing introns and ligating exons. This is a highly regulated process, but any alteration in splicing factors, splicing sites, or auxiliary components affects the final products of the gene. In diffuse large B-cell lymphoma, splicing mutations such as mutant splice sites, aberrant alternative splicing, exon skipping, and intron retention are detected. The alteration affects tumor suppression, DNA repair, cell cycle, cell differentiation, cell proliferation, and apoptosis. As a result, malignant transformation, cancer progression, and metastasis occurred in B cells at the germinal center. B-cell lymphoma 7 protein family member A (BCL7A), cluster of differentiation 79B (CD79B), myeloid differentiation primary response gene 88 (MYD88), tumor protein P53 (TP53), signal transducer and activator of transcription (STAT), serum- and glucose-regulated kinase 1 (SGK1), Pou class 2 associating factor 1 (POU2AF1), and neurogenic locus notch homolog protein 1 (NOTCH) are the most common genes affected by splicing mutations in diffuse large B cell lymphoma.

Detection of germline variants with pathogenic potential in 48 patients with familial colorectal cancer by using whole exome sequencing

BACKGROUND

Hereditary genetic mutations causing predisposition to colorectal cancer are accountable for approximately 30% of all colorectal cancer cases. However, only a small fraction of these are high penetrant mutations occurring in DNA mismatch repair genes, causing one of several types of familial colorectal cancer (CRC) syndromes. Most of the mutations are low-penetrant variants, contributing to an increased risk of familial colorectal cancer, and they are often found in additional genes and pathways not previously associated with CRC. The aim of this study was to identify such variants, both high-penetrant and low-penetrant ones.

METHODS

We performed whole exome sequencing on constitutional DNA extracted from blood of 48 patients suspected of familial colorectal cancer and used multiple in silico prediction tools and available literature-based evidence to detect and investigate genetic variants.

RESULTS

We identified several causative and some potentially causative germline variants in genes known for their association with colorectal cancer. In addition, we identified several variants in genes not typically included in relevant gene panels for colorectal cancer, including CFTR, PABPC1 and TYRO3, which may be associated with an increased risk for cancer.

CONCLUSIONS

Identification of variants in additional genes that potentially can be associated with familial colorectal cancer indicates a larger genetic spectrum of this disease, not limited only to mismatch repair genes. Usage of multiple in silico tools based on different methods and combined through a consensus approach increases the sensitivity of predictions and narrows down a large list of variants to the ones that are most likely to be significant.

The role of noncoding genetic variants in cardiomyopathy

Cardiomyopathies remain one of the leading causes of morbidity and mortality worldwide. Environmental risk factors and genetic predisposition account for most cardiomyopathy cases. As with all complex diseases, there are significant challenges in the interpretation of the molecular mechanisms underlying cardiomyopathy-associated genetic variants. Given the technical improvements and reduced costs of DNA sequence technologies, an increasing number of patients are now undergoing genetic testing, resulting in a continuously expanding list of novel mutations. However, many patients carry noncoding genetic variants, and although emerging evidence supports their contribution to cardiac disease, their role in cardiomyopathies remains largely understudied. In this review, we summarize published studies reporting on the association of different types of noncoding variants with various types of cardiomyopathies. We focus on variants within transcriptional enhancers, promoters, intronic sites, and untranslated regions that are likely associated with cardiac disease. Given the broad nature of this topic, we provide an overview of studies that are relatively recent and have sufficient evidence to support a significant degree of causality. We believe that more research with additional validation of noncoding genetic variants will provide further mechanistic insights on the development of cardiac disease, and noncoding variants will be increasingly incorporated in future genetic screening tests.

Splicing accuracy varies across human introns, tissues and age

Alternative splicing impacts most multi-exonic human genes. Inaccuracies during this process may have an important role in ageing and disease. Here, we investigated mis-splicing using RNA-sequencing data from ~14K control samples and 42 human body sites, focusing on split reads partially mapping to known transcripts in annotation. We show that mis-splicing occurs at different rates across introns and tissues and that these splicing inaccuracies are primarily affected by the abundance of core components of the spliceosome assembly and its regulators. Using publicly available data on short-hairpin RNA-knockdowns of numerous spliceosomal components and related regulators, we found support for the importance of RNA-binding proteins in mis-splicing. We also demonstrated that age is positively correlated with mis-splicing, and it affects genes implicated in neurodegenerative diseases. This in-depth characterisation of mis-splicing can have important implications for our understanding of the role of splicing inaccuracies in human disease and the interpretation of long-read RNA-sequencing data.

Pathogenicity of Intronic and Synonymous Variants of ATP7B in Wilson Disease

Wilson disease (WD) is a hereditary disorder of copper metabolism, resulting from mutations within ATP7B. Early diagnosis is essential for affected individuals. However, there are still patients with clinically suspected WD who do not have detectable pathogenic variants, which makes diagnosis difficult and delays treatment. This study included such patients from the authors' center and screened for the full-length sequence of ATP7B by next-generation sequencing. Newly identified synonymous and intronic variants were then analyzed with in silico tools. A minigene system was constructed to determine the pathogenicity of these variants in terms of splicing and blood RNA extraction, and RT-PCR experiments were performed on several patients to verify the splicing alterations. The phenotypes of the patients were also analyzed. Fourteen suspected pathogenic variants, including nine synonymous and five intronic variants, were detected in 12 patients with clinically suspected WD. Among them, four synonymous variants (c.1050G>A, c.1122C>G, c.3243G>A, and c.4014T>A) and four intronic variants (c.1543 +40G>A, c.1707+6_1707+16del, c.1870-49A>G, and c.2731-67A>G) resulted in splicing changes in ATP7B. After the above analysis, the diagnosis of WD could be confirmed in eight clinically suspected patients with WD who showed a late age of onset.

Aberrant splicing in human cancer: An RNA structural code point of view

Alternative splicing represents an essential process that occurs widely in eukaryotes. In humans, most genes undergo alternative splicing to ensure transcriptome and proteome diversity reflecting their functional complexity. Over the last decade, aberrantly spliced transcripts due to mutations in cis- or trans-acting splicing regulators have been tightly associated with cancer development, largely drawing scientific attention. Although a plethora of single proteins, ribonucleoproteins, complexed RNAs, and short RNA sequences have emerged as nodal contributors to the splicing cascade, the role of RNA secondary structures in warranting splicing fidelity has been underestimated. Recent studies have leveraged the establishment of novel high-throughput methodologies and bioinformatic tools to shed light on an additional layer of splicing regulation in the context of RNA structural elements. This short review focuses on the most recent available data on splicing mechanism regulation on the basis of RNA secondary structure, emphasizing the importance of the complex RNA G-quadruplex structures (rG4s), and other specific RNA motifs identified as splicing silencers or enhancers. Moreover, it intends to provide knowledge on newly established techniques that allow the identification of RNA structural elements and highlight the potential to develop new RNA-oriented therapeutic strategies against cancer.

Computational prediction of human deep intronic variation

BACKGROUND

The adoption of whole-genome sequencing in genetic screens has facilitated the detection of genetic variation in the intronic regions of genes, far from annotated splice sites. However, selecting an appropriate computational tool to discriminate functionally relevant genetic variants from those with no effect is challenging, particularly for deep intronic regions where independent benchmarks are scarce.

RESULTS

In this study, we have provided an overview of the computational methods available and the extent to which they can be used to analyze deep intronic variation. We leveraged diverse datasets to extensively evaluate tool performance across different intronic regions, distinguishing between variants that are expected to disrupt splicing through different molecular mechanisms. Notably, we compared the performance of SpliceAI, a widely used sequence-based deep learning model, with that of more recent methods that extend its original implementation. We observed considerable differences in tool performance depending on the region considered, with variants generating cryptic splice sites being better predicted than those that potentially affect splicing regulatory elements. Finally, we devised a novel quantitative assessment of tool interpretability and found that tools providing mechanistic explanations of their predictions are often correct with respect to the ground - information, but the use of these tools results in decreased predictive power when compared to black box methods.

CONCLUSIONS

Our findings translate into practical recommendations for tool usage and provide a reference framework for applying prediction tools in deep intronic regions, enabling more informed decision-making by practitioners.

POU1F1/Pou1f1 c.143-83A > G Variant Disrupts the Branch Site in Pre-mRNA and Leads to Dwarfism

POU Class 1 Homeobox1 (POU1F1/Pou1f1) is a well-established pituitary-specific transcription factor, and causes, when mutated, combined pituitary hormone deficiency in humans and mice. POU1F1/Pou1f1 has 2 isoforms: the alpha and beta isoforms. Recently, pathogenic variants in the unique coding region of the beta isoform (beta domain) and the intron near the exon-intron boundary for the beta domain were reported, although their functional consequences remain obscure. In this study, we generated mice carrying the Pou1f1 c.143-83A>G substitution that recapitulates the human intronic variant near the exon-intron boundary for the beta domain. Homozygous mice showed postnatal growth failure, with an average body weight that was 35% of wild-type littermates at 12 weeks, which was accompanied by anterior pituitary hypoplasia and deficiency of circulating insulin-like growth factor 1 and thyroxine. The results of RNA-seq analysis of the pituitary gland were consistent with reduction of somatotrophs, and this was confirmed immunohistochemically. Reverse transcription polymerase chain reaction of pituitary Pou1f1 mRNA showed abnormal splicing in homozygous mice, with a decrease in the alpha isoform, an increase in the beta isoform, and the emergence of the exon-skipped transcript. We further characterized artificial variants in or near the beta domain, which were candidate positions of the branch site in pre-mRNA, using cultured cell-basis analysis and found that only c.143-83A>G produced transcripts similar to the mice model. Our report is the first to show that the c.143-83A>G variant leads to splicing disruption and causes morphological and functional abnormalities in the pituitary gland. Furthermore, our mice will contribute understanding the role of POU1F1/Pou1f1 transcripts in pituitary development.

Transcriptome analysis provides critical answers to the "variants of uncertain significance" conundrum

While whole-genome and exome sequencing have transformed our collective understanding of genetics' role in disease pathogenesis, there are certain conditions and populations for whom DNA-level data fails to identify the underlying genetic etiology. Specifically, patients of non-White race and non-European ancestry are disproportionately affected by "variants of unknown/uncertain significance" (VUS), limiting the scope of precision medicine for minority patients and perpetuating health disparities. VUS often include deep intronic and splicing variants which are difficult to interpret from DNA data alone. RNA analysis can illuminate the consequences of VUS, thereby allowing for their reclassification as pathogenic versus benign. Here we review the critical role transcriptome analysis plays in clarifying VUS in both neoplastic and non-neoplastic diseases.

Rb Tumor Suppressor in Small Cell Lung Cancer: Combined Genomic and IHC Analysis with a Description of a Distinct Rb-Proficient Subset

PURPOSE

RB1 mutations and loss of retinoblastoma (Rb) expression represent consistent but not entirely invariable hallmarks of small cell lung cancer (SCLC). The prevalence and characteristics of SCLC retaining wild-type Rb are not well-established. Furthermore, the performance of targeted next-generation sequencing (NGS) versus immunohistochemistry for Rb assessment is not well-defined.

EXPERIMENTAL DESIGN

A total of 208 clinical SCLC samples were analyzed by comprehensive targeted NGS, covering all exons of RB1, and Rb IHC. On the basis of established coordination of Rb/p16/cyclinD1 expression, p16-high/cyclinD1-low profile was used as a marker of constitutive Rb deficiency.

RESULTS

Fourteen of 208 (6%) SCLC expressed wild-type Rb, accompanied by a unique p16-low/cyclinD1-high profile supporting Rb proficiency. Rb-proficient SCLC was associated with neuroendocrine-low phenotype, combined SCLC with non-SCLC (NSCLC) histology and aggressive behavior. These tumors exclusively harbored CCND1 amplification (29%), and were markedly enriched in CDKN2A mutations (50%) and NSCLC-type alterations (KEAP1, STK11, FGFR1). The remaining 194 of 208 SCLC were Rb-deficient (p16-high/cyclinD1-low), including 184 cases with Rb loss (of which 29% lacked detectable RB1 alterations by clinical NGS pipeline), and 10 cases with mutated but expressed Rb.

CONCLUSIONS

This is the largest study to date to concurrently analyze Rb by NGS and IHC in SCLC, identifying a 6% rate of Rb proficiency. Pathologic-genomic data implicate NSCLC-related progenitors as a putative source of Rb-proficient SCLC. Consistent upstream Rb inactivation via CDKN2A/p16↓ and CCND1/cyclinD1↑ suggests the potential utility of CDK4/6 inhibitors in this aggressive SCLC subset. The study also clarifies technical aspects of Rb status determination in clinical practice, highlighting the limitations of exon-only sequencing for RB1 interrogation. See related commentary by Mahadevan and Sholl, p. 4603.

Genetic risk variants for multiple sclerosis are linked to differences in alternative pre-mRNA splicing

BACKGROUND

Multiple sclerosis (MS) is a chronic immune-mediated disease of the central nervous system to which a genetic predisposition contributes. Over 200 genetic regions have been associated with increased disease risk, but the disease-causing variants and their functional impact at the molecular level are mostly poorly defined. We hypothesized that single-nucleotide polymorphisms (SNPs) have an impact on pre-mRNA splicing in MS.

METHODS

Our study focused on 10 bioinformatically prioritized SNP-gene pairs, in which the SNP has a high potential to alter alternative splicing events (ASEs). We tested for differential gene expression and differential alternative splicing in B cells from MS patients and healthy controls. We further examined the impact of the SNP genotypes on ASEs and on splice isoform expression levels. Novel genotype-dependent effects on splicing were verified with splicing reporter minigene assays.

RESULTS

We were able to confirm previously described findings regarding the relation of MS-associated SNPs with the ASEs of the pre-mRNAs from GSDMB and SP140. We also observed an increased IL7R exon 6 skipping when comparing relapsing and progressive MS patients to healthy subjects. Moreover, we found evidence that the MS risk alleles of the SNPs rs3851808 (EFCAB13), rs1131123 (HLA-C), rs10783847 (TSFM), and rs2014886 (TSFM) may contribute to a differential splicing pattern. Of particular interest is the genotype-dependent exon skipping of TSFM due to the SNP rs2014886. The minor allele T creates a donor splice site, resulting in the expression of the exon 3 and 4 of a short TSFM transcript isoform, whereas in the presence of the MS risk allele C, this donor site is absent, and thus the short transcript isoform is not expressed.

CONCLUSION

In summary, we found that genetic variants from MS risk loci affect pre-mRNA splicing. Our findings substantiate the role of ASEs with respect to the genetics of MS. Further studies on how disease-causing genetic variants may modify the interactions between splicing regulatory sequence elements and RNA-binding proteins can help to deepen our understanding of the genetic susceptibility to MS.

In Silico Study to Predict the Structural and Functional Consequences of SNPs on Biomarkers of Ovarian Cancer (OC) and BPA Exposure-Associated OC

Background: Recently, we have shown that seven genes, namely GBP5, IRS2, KRT4, LINCOO707, MRPL55, RRS1 and SLC4A11, have prognostic power for the overall survival in ovarian cancer (OC). Methods: We present an analysis on the association of these genes with any phenotypes and mutations indicative of involvement in female cancers and predict the structural and functional consequences of those SNPS using in silico tools. Results: These seven genes present with 976 SNPs/mutations that are associated with human cancers, out of which 284 related to female cancers. We have then analysed the mutation impact on amino acid polarity, charge and water affinity, leading to the identification of 30 mutations in gynaecological cancers where amino acid (aa) changes lead to opposite polarity, charges and water affinity. Out of these 30 mutations identified, only a missense mutation (i.e., R831C/R804C in uterine corpus endometrial carcinomas, UCEC) was suggestive of structural damage on the SLC4A11 protein. Conclusions: We demonstrate that the R831C/R804C mutation is deleterious and the predicted ΔΔG values suggest that the mutation reduces the stability of the protein. Future in vitro studies should provide further insight into the role of this transporter protein in UCEC.

Pseudoexon activation in disease by non-splice site deep intronic sequence variation - wild type pseudoexons constitute high-risk sites in the human genome

Accuracy of pre-messenger RNA (pre-mRNA) splicing is crucial for normal gene expression. Complex regulation supports the spliceosomal distinction between authentic exons and the many seemingly functional splice sites delimiting pseudoexons. Pseudoexons are nonfunctional intronic sequences that can be activated for aberrant inclusion in mRNA, which may cause disease. Pseudoexon activation is very challenging to predict, in particular when activation occurs by sequence variants that alter the splicing regulatory environment without directly affecting splice sites. As pseudoexon inclusion often evades detection due to activation of nonsense-mediated mRNA decay, and because conventional diagnostic procedures miss deep intronic sequence variation, pseudoexon activation is a heavily underreported disease mechanism. Pseudoexon characteristics have mainly been studied based on in silico predicted sequences. Moreover, because recognition of sequence variants that create or strengthen splice sites is possible by comparison with well-established consensus sequences, this type of pseudoexon activation is by far the most frequently reported. Here we review all known human disease-associated pseudoexons that carry functional splice sites and are activated by deep intronic sequence variants located outside splice site sequences. We delineate common characteristics that make this type of wild type pseudoexons distinct high-risk sites in the human genome.

Role of Splicing Regulatory Elements and In Silico Tools Usage in the Identification of Deep Intronic Splicing Variants in Hereditary Breast/Ovarian Cancer Genes

Simple Summary

There is a significant percentage of hereditary breast and ovarian cancer (HBOC) cases that remain undiagnosed, because no pathogenic variant is detected through massively parallel sequencing of coding exons and exon-intron boundaries of high-moderate susceptibility risk genes. Deep intronic regions may contain variants affecting RNA splicing, leading ultimately to disease, and hence they may explain several cases where the genetic cause of HBOC is unknown. This study aims to characterize intronic regions to identify a landscape of “exonizable” zones and test the efficiency of two in silico tools to detect deep intronic variants affecting the mRNA splicing process.

Abstract

The contribution of deep intronic splice-altering variants to hereditary breast and ovarian cancer (HBOC) is unknown. Current computational in silico tools to predict spliceogenic variants leading to pseudoexons have limited efficiency. We assessed the performance of the SpliceAI tool combined with ESRseq scores to identify spliceogenic deep intronic variants by affecting cryptic sites or splicing regulatory elements (SREs) using literature and experimental datasets. Our results with 233 published deep intronic variants showed that SpliceAI, with a 0.05 threshold, predicts spliceogenic deep intronic variants affecting cryptic splice sites, but is less effective in detecting those affecting SREs. Next, we characterized the SRE profiles using ESRseq, showing that pseudoexons are significantly enriched in SRE-enhancers compared to adjacent intronic regions. Although the combination of SpliceAI with ESRseq scores (considering ∆ESRseq and SRE landscape) showed higher sensitivity, the global performance did not improve because of the higher number of false positives. The combination of both tools was tested in a tumor RNA dataset with 207 intronic variants disrupting splicing, showing a sensitivity of 86%. Following the pipeline, five spliceogenic deep intronic variants were experimentally identified from 33 variants in HBOC genes. Overall, our results provide a framework to detect deep intronic variants disrupting splicing.

Splicing in the Diagnosis of Rare Disease: Advances and Challenges

Mutations which affect splicing are significant contributors to rare disease, but are frequently overlooked by diagnostic sequencing pipelines. Greater ascertainment of pathogenic splicing variants will increase diagnostic yields, ending the diagnostic odyssey for patients and families affected by rare disorders, and improving treatment and care strategies. Advances in sequencing technologies, predictive modeling, and understanding of the mechanisms of splicing in recent years pave the way for improved detection and interpretation of splice affecting variants, yet several limitations still prohibit their routine ascertainment in diagnostic testing. This review explores some of these advances in the context of clinical application and discusses challenges to be overcome before these variants are comprehensively and routinely recognized in diagnostics.