Improved detection of aberrant splicing with FRASER 2.0 and the intron Jaccard index.

Detection of aberrantly spliced genes is an important step in RNA-seq-based rare-disease diagnostics. We recently developed FRASER, a denoising autoencoder-based method that outperformed alternative methods of detecting aberrant splicing. However, because FRASER's three splice metrics are partially redundant and tend to be sensitive to sequencing depth, we introduce here a more robust intron-excision metric, the intron Jaccard index, that combines the alternative donor, alternative acceptor, and intron-retention signal into a single value. Moreover, we optimized model parameters and filter cutoffs by using candidate rare-splice-disrupting variants as independent evidence. On 16,213 GTEx samples, our improved algorithm, FRASER 2.0, called typically 10 times fewer splicing outliers while increasing the proportion of candidate rare-splice-disrupting variants by 10-fold and substantially decreasing the effect of sequencing depth on the number of reported outliers. To lower the multiple-testing correction burden, we introduce an option to select the genes to be tested for each sample instead of a transcriptome-wide approach. This option can be particularly useful when prior information, such as candidate variants or genes, is available. Application on 303 rare-disease samples confirmed the relative reduction in the number of outlier calls for a slight loss of sensitivity; FRASER 2.0 recovered 22 out of 26 previously identified pathogenic splicing cases with default cutoffs and 24 when multiple-testing correction was limited to OMIM genes containing rare variants. Altogether, these methodological improvements contribute to more effective RNA-seq-based rare diagnostics by drastically reducing the amount of splicing outlier calls per sample at minimal loss of sensitivity.

Computational approaches for detecting disease-associated alternative splicing events.

Alternative splicing (AS) is a key transcriptional regulation pathway. Recent studies have shown that AS events are associated with the occurrence of complex diseases. Various computational approaches have been developed for the detection of disease-associated AS events. In this review, we first describe the metrics used for quantitative characterization of AS events. Second, we review and discuss the three types of methods for detecting disease-associated splicing events, which are differential splicing analysis, aberrant splicing detection and splicing-related network analysis. Third, to further exploit the genetic mechanism of disease-associated AS events, we describe the methods for detecting genetic variants that potentially regulate splicing. For each type of methods, we conducted experimental comparison to illustrate their performance. Finally, we discuss the limitations of these methods and point out potential ways to address them. We anticipate that this review provides a systematic understanding of computational approaches for the analysis of disease-associated splicing.

Comprehensive splicing analysis of the alternatively spliced CHEK2 exons 8 and 10 reveals three enhancer/silencer-rich regions and 38 spliceogenic variants.

Splicing is controlled by a large set of regulatory elements (SREs) including splicing enhancers and silencers, which are involved in exon recognition. Variants at these motifs may dysregulate splicing and trigger loss-of-function transcripts associated with disease. Our goal here was to study the alternatively spliced exons 8 and 10 of the breast cancer susceptibility gene CHEK2. For this purpose, we used a previously published minigene with exons 6-10 that produced the expected minigene full-length transcript and replicated the naturally occurring events of exon 8 [Δ(E8)] and exon 10 [Δ(E10)] skipping. We then introduced 12 internal microdeletions of exons 8 and 10 by mutagenesis in order to map SRE-rich intervals by splicing assays in MCF-7 cells. We identified three minimal (10-, 11-, 15-nt) regions essential for exon recognition: c.863_877del [ex8, Δ(E8): 75%] and c.1073_1083del and c.1083_1092del [ex10, Δ(E10): 97% and 62%, respectively]. Then 87 variants found within these intervals were introduced into the wild-type minigene and tested functionally. Thirty-eight of them (44%) impaired splicing, four of which (c.883G>A, c.883G>T, c.884A>T, and c.1080G>T) induced negligible amounts (<5%) of the minigene full-length transcript. Another six variants (c.886G>A, c.886G>T, c.1075G>A, c.1075G>T, c.1076A>T, and c.1078G>T) showed significantly strong impacts (20-50% of the minigene full-length transcript). Thirty-three of the 38 spliceogenic variants were annotated as missense, three as nonsense, and two as synonymous, underlying the fact that any exonic change is capable of disrupting splicing. Moreover, c.883G>A, c.883G>T, and c.884A>T were classified as pathogenic/likely pathogenic variants according to ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based criteria. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Exon 1-targeting miRNA reduces the pathogenic exon 1 HTT protein in Huntington disease models.

Huntington disease (HD) is a fatal neurodegenerative disease caused by a trinucleotide repeat expansion in exon 1 of the huntingtin gene (HTT) resulting in toxic gain-of-function and cell death. Despite its monogenic cause, the pathogenesis of HD is highly complex and increasing evidence indicates that, in addition to the full-length (FL) mutant HTT protein, the expanded exon 1 HTT (HTTexon1) protein that is translated from the HTT1a transcript generated by aberrant splicing is prone to aggregate and may contribute to HD pathology. This finding suggests that reducing the expression of HTT1a may achieve a greater therapeutic benefit than targeting only FL mutant HTT. Conversely, strategies that exclusively target FL HTT may not fully prevent the pathogenesis of HD. We have developed an engineered microRNA targeting the HTT exon 1 sequence (miHTT), delivered via adeno-associated virus serotype 5 (AAV5). The target sequence of miHTT is present in both FL HTT and HTT1a transcripts. Preclinical studies with AAV5-miHTT have demonstrated efficacy in several rodent and large animal models by reducing FL HTT mRNA and protein and rescuing HD-like phenotypes, and have been the rationale for phase I/II clinical studies now ongoing in the US and Europe. In the present study, we evaluated the ability of AAV5-miHTT to reduce the levels of aberrantly spliced HTT1a mRNA and the HTTexon1 protein in the brain of two mouse models of HD (heterozygous zQ175 knock-in mice and humanized Hu128/21 mice). Polyadenylated HTT1a mRNA and HTTexon1 protein were detected in the striatum and cortex of heterozygous zQ175 knock-in mice, but not in wild-type, littermate control mice. Intrastriatal administration of AAV5-miHTT resulted in dose-dependent expression of mature miHTT microRNA in cortical brain regions, accompanied by significant lowering of both FL HTT and HTT1a mRNA expression at two months post-injection. Mutant HTT and HTTexon1 protein levels were also significantly reduced in the striatum and cortex of heterozygous zQ175 knock-in at 2 months after AAV5-miHTT treatment and in humanized Hu128/21 mice 7 months post-treatment. The effects were confirmed in primary Hu128/21 neuronal cultures. These results demonstrate that AAV5-miHTT gene therapy is an effective approach to lower both FL HTT and the pathogenic HTTexon1 levels, which could potentially have an additive therapeutic benefit compared to other HTT-targeting modalities.

Loss-of-function variants affecting the STAGA complex component SUPT7L cause a developmental disorder with generalized lipodystrophy.

Generalized lipodystrophy is a feature of various hereditary disorders, often leading to a progeroid appearance. In the present study we identified a missense and a frameshift variant in a compound heterozygous state in SUPT7L in a boy with intrauterine growth retardation, generalized lipodystrophy, and additional progeroid features. SUPT7L encodes a component of the transcriptional coactivator complex STAGA. By transcriptome sequencing, we showed the predicted missense variant to cause aberrant splicing, leading to exon truncation and thereby to a complete absence of SUPT7L in dermal fibroblasts. In addition, we found altered expression of genes encoding DNA repair pathway components. This pathway was further investigated and an increased rate of DNA damage was detected in proband-derived fibroblasts and genome-edited HeLa cells. Finally, we performed transient overexpression of wildtype SUPT7L in both cellular systems, which normalizes the number of DNA damage events. Our findings suggest SUPT7L as a novel disease gene and underline the link between genome instability and progeroid phenotypes.

Therapeutic strategies for aberrant splicing in cancer and genetic disorders.

Accurate pre-mRNA splicing is essential for proper protein translation; however, aberrant splicing is commonly observed in the context of cancer and genetic disorders. Notably, in genetic diseases, these splicing abnormalities often play a pivotal role. Substantial challenges persist in accurately identifying and classifying disease-induced aberrant splicing, as well as in development of targeted therapeutic strategies. In this review, we examine prevalent forms of aberrant splicing and explore potential therapeutic approaches aimed at addressing these splicing-related diseases. This summary contributes to a deeper understanding of the complexities about aberrant splicing and provide a foundation for the development of effective therapeutic interventions in the field of genetic disorders and cancer.

RNA structure in alternative splicing regulation: from mechanism to therapy.

Alternative splicing is a highly intricate process that plays a crucial role in post-transcriptional regulation and significantly expands the functional proteome of a limited number of coding genes in eukaryotes. Its regulation is multifactorial, with RNA structure exerting a significant impact. Aberrant RNA conformations lead to dysregulation of splicing patterns, which directly affects the manifestation of disease symptoms. In this review, the molecular mechanisms of RNA secondary structure-mediated splicing regulation are summarized, with a focus on the complex interplay between aberrant RNA conformations and disease phenotypes resulted from splicing defects. This study also explores additional factors that reshape structural conformations, enriching our understanding of the mechanistic network underlying structure-mediated splicing regulation. In addition, an emphasis has been placed on the clinical role of targeting aberrant splicing corrections in human diseases. The principal mechanisms of action behind this phenomenon are described, followed by a discussion of prospective development strategies and pertinent challenges.

Parallel DNA/RNA NGS Using an Identical Target Enrichment Panel in the Analysis of Hereditary Cancer Predisposition.

Germline DNA testing using the next-gene-ration sequencing (NGS) technology has become the analytical standard for the diagnostics of hereditary diseases, including cancer. Its increasing use places high demands on correct sample identification, independent confirmation of prioritized variants, and their functional and clinical interpretation. To streamline these processes, we introduced parallel DNA and RNA capture-based NGS using identical capture panel CZECANCA, which is routinely used for DNA analysis of hereditary cancer predisposition. Here, we present the analytical workflow for RNA sample processing and its analytical and diagnostic performance. Parallel DNA/RNA analysis allowed credible sample identification by calculating the kinship coefficient. The RNA capture-based approach enriched transcriptional targets for the majority of clinically relevant cancer predisposition genes to a degree that allowed analysis of the effect of identified DNA variants on mRNA processing. By comparing the panel and whole-exome RNA enrichment, we demonstrated that the tissue-specific gene expression pattern is independent of the capture panel. Moreover, technical replicates confirmed high reproducibility of the tested RNA analysis. We concluded that parallel DNA/RNA NGS using the identical gene panel is a robust and cost-effective diagnostic strategy. In our setting, it allows routine analysis of 48 DNA/RNA pairs using NextSeq 500/550 Mid Output Kit v2.5 (150 cycles) in a single run with sufficient coverage to analyse 226 cancer predisposition and candidate ge-nes. This approach can replace laborious Sanger confirmatory sequencing, increase testing turnaround, reduce analysis costs, and improve interpretation of the impact of variants by analysing their effect on mRNA processing.

Aberrant splicing of mutant huntingtin in Huntington's disease knock-in pigs.

Huntington's disease (HD) is an autosomal-dominant inherited neurodegenerative disease caused by a CAG repeat expansion in exon1 of the huntingtin gene (HTT). This expansion leads to the production of N-terminal mutant huntingtin protein (mHtt) that contains an expanded polyglutamine tract, which is toxic to neurons and causes neurodegeneration. While the production of N-terminal mHtt can be mediated by proteolytic cleavage of full-length mHtt, abnormal splicing of exon1-intron1 of mHtt has also been identified in the brains of HD mice and patients. However, the proportion of aberrantly spliced exon1 mHTT in relation to normal mHTT exon remains to be defined. In this study, HTT exon1 production was examined in the HD knock-in (KI) pig model, which more closely recapitulates neuropathology seen in HD patient brains than HD mouse models. The study revealed that aberrant spliced HTT exon1 is also present in the brains of HD pigs, but it is expressed at a much lower level than the normally spliced HTT exon products. These findings suggest that careful consideration is needed when assessing the contribution of aberrantly spliced mHTT exon1 to HD pathogenesis, and further rigorous investigation is required.

Splicing predictions, minigene analyses, and ACMG-AMP clinical classification of 42 germline PALB2 splice-site variants.

PALB2 loss-of-function variants confer high risk of developing breast cancer. Here we present a systematic functional analysis of PALB2 splice-site variants detected in approximately 113,000 women in the large-scale sequencing project Breast Cancer After Diagnostic Gene Sequencing (BRIDGES; https://bridges-research.eu/). Eighty-two PALB2 variants at the intron-exon boundaries were analyzed with MaxEntScan. Forty-two variants were selected for the subsequent splicing functional assays. For this purpose, three splicing reporter minigenes comprising exons 1-12 were constructed. The 42 potential spliceogenic variants were introduced into the minigenes by site-directed mutagenesis and assayed in MCF-7/MDA-MB-231 cells. Splicing anomalies were observed in 35 variants, 23 of which showed no traces or minimal amounts of the expected full-length transcripts of each minigene. More than 30 different variant-induced transcripts were characterized, 23 of which were predicted to truncate the PALB2 protein. The pathogenicity of all variants was interpreted according to an in-house adaptation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG-AMP) variant classification scheme. Up to 23 variants were classified as pathogenic/likely pathogenic. Remarkably, three ±1,2 variants (c.49-2A>T, c.108+2T>C, and c.211+1G>A) were classified as variants of unknown significance, as they produced significant amounts of either in-frame transcripts of unknown impact on the PALB2 protein function or the minigene full-length transcripts. In conclusion, we have significantly contributed to the ongoing effort of identifying spliceogenic variants in the clinically relevant PALB2 cancer susceptibility gene. Moreover, we suggest some approaches to classify the findings in accordance with the ACMG-AMP rationale. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Minigene-based splicing analysis and ACMG/AMP-based tentative classification of 56 ATM variants.

The ataxia telangiectasia-mutated (ATM) protein is a major coordinator of the DNA damage response pathway. ATM loss-of-function variants are associated with 2-fold increased breast cancer risk. We aimed at identifying and classifying spliceogenic ATM variants detected in subjects of the large-scale sequencing project BRIDGES. A total of 381 variants at the intron-exon boundaries were identified, 128 of which were predicted to be spliceogenic. After further filtering, we ended up selecting 56 variants for splicing analysis. Four functional minigenes (mgATM) spanning exons 4-9, 11-17, 25-29, and 49-52 were constructed in the splicing plasmid pSAD. Selected variants were genetically engineered into the four constructs and assayed in MCF-7/HeLa cells. Forty-eight variants (85.7%) impaired splicing, 32 of which did not show any trace of the full-length (FL) transcript. A total of 43 transcripts were identified where the most prevalent event was exon/multi-exon skipping. Twenty-seven transcripts were predicted to truncate the ATM protein. A tentative ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme that integrates mgATM data allowed us to classify 29 ATM variants as pathogenic/likely pathogenic and seven variants as likely benign. Interestingly, the likely pathogenic variant c.1898+2T>G generated 13% of the minigene FL-transcript due to the use of a noncanonical GG-5'-splice-site (0.014% of human donor sites). Circumstantial evidence in three ATM variants (leakiness uncovered by our mgATM analysis together with clinical data) provides some support for a dosage-sensitive expression model in which variants producing ≥30% of FL-transcripts would be predicted benign, while variants producing ≤13% of FL-transcripts might be pathogenic. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

All reported non-canonical splice site variants in GLA cause aberrant splicing.

BACKGROUND: Fabry disease is an X-linked lysosomal storage disorder caused by insufficient α-galactosidase A (GLA) activity resulting from variants in the GLA gene, which leads to glycosphingolipid accumulation and life-threatening, multi-organ complications. Approximately 50 variants have been reported that cause splicing abnormalities in GLA. Most were found within canonical splice sites, which are highly conserved GT and AG splice acceptor and donor dinucleotides, whereas one-third were located outside canonical splice sites, making it difficult to interpret their pathogenicity. In this study, we aimed to investigate the genetic pathogenicity of variants located in non-canonical splice sites within the GLA gene. METHODS: 13 variants, including four deep intronic variants, were selected from the Human Gene Variant Database Professional. We performed an in vitro splicing assay to identify splicing abnormalities in the variants. RESULTS: All candidate non-canonical splice site variants in GLA caused aberrant splicing. Additionally, all but one variant was protein-truncating. The four deep intronic variants generated abnormal transcripts, including a cryptic exon, as well as normal transcripts, with the proportion of each differing in a cell-specific manner. CONCLUSIONS: Validation of splicing effects using an in vitro splicing assay is useful for confirming pathogenicity and determining associations with clinical phenotypes.

Aberrant splicing caused by exonic single nucleotide variants positioned 2nd or 3rd to the last nucleotide in the COL4A5 gene.

BACKGROUND AND OBJECTIVES: The evident genotype-phenotype correlation shown by the X-linked Alport syndrome warrants the assessment of the impact of identified gene variants on aberrant splicing. We previously reported that single nucleotide variants (SNVs) in the last nucleotide of exons in COL4A5 cause aberrant splicing. It is known that the nucleotides located 2nd and 3rd to the last nucleotides of exons can also play an essential role in the first step of the splicing process. In this study, we aimed to investigate whether SNVs positioned 2nd or 3rd to the last nucleotide of exons in COL4A5 resulted in aberrant splicing. METHODS: We selected eight candidate variants: six from the Human Gene Variant Database Professional and two from our cohort. We performed an in-vitro splicing assay and reverse transcription-polymerase chain reaction (RT-PCR) for messenger RNA obtained from patients, if available. RESULTS: The candidate variants were initially classified into the following groups: three nonsense, two missense, and three synonymous variants. Splicing assays and RT-PCR for messenger RNA revealed that six of the eight variants caused aberrant splicing. Four variants, initially classified as non-truncating variants, were found to be truncating ones, which usually show relatively more severe phenotypes. CONCLUSION: We revealed that exonic SNVs positioned 2nd or 3rd to the last nucleotide of exons in the COL4A5 were responsible for aberrant splicing. The results of our study suggest that attention should be paid when interpreting the pathogenicity of exonic SNVs near the 5' splice site.

Frameshift coding sequence variants in the LPL gene: identification of two novel events and exploration of the genotype-phenotype relationship for variants reported to date.

BACKGROUND: Lipoprotein lipase (LPL) is the rate-limiting enzyme for triglyceride hydrolysis. Homozygous or compound heterozygous LPL variants cause autosomal recessive familial chylomicronemia syndrome (FCS), whereas simple heterozygous LPL variants are associated with hypertriglyceridemia (HTG) and HTG-related disorders. LPL frameshift coding sequence variants usually cause complete functional loss of the affected allele, thereby allowing exploration of the impact of different levels of LPL function in human disease. METHODS: All exons and flanking intronic regions of LPL were Sanger sequenced in patients with HTG-related acute pancreatitis (HTG-AP) or HTG-AP in pregnancy. Previously reported LPL frameshift coding sequence variants were collated from the Human Gene Mutation Database and through PubMed keyword searching. Original reports were manually evaluated for the following information: zygosity status of the variant, plasma LPL activity of the variant carrier, disease referred for genetic analysis, patient's age at genetic analysis, and patient's disease history. SpliceAI was employed to predict the potential impact of collated variants on splicing. RESULTS: Two novel rare variants were identified, and 53 known LPL frameshift coding sequence variants were collated. Of the 51 variants informative for zygosity, 30 were simple heterozygotes, 12 were homozygotes, and 9 were compound heterozygotes. Careful evaluation of the 55 variants with respect to their clinical and genetic data generated several interesting findings. First, we conclude that 6-7% residual LPL function could significantly delay the age of onset of FCS and reduce the prevalence of FCS-associated syndromes. Second, whereas a large majority of LPL frameshift coding sequence variants completely disrupt gene function through their "frameshift" nature, a small fraction of these variants may act wholly or partly as "in-frame" variants, leading to the generation of protein products with some residual LPL function. Third, we identified two candidate LPL frameshift coding sequence variants that may retain residual function based on genotype-phenotype correlation or SpliceAI-predicted data. CONCLUSIONS: This study reported two novel LPL variants and yielded new insights into the genotype-phenotype relationship as it pertains to LPL frameshift coding sequence variants.

A Novel Intronic Deletion in PDE6B Causes Autosomal Recessive Retinitis Pigmentosa by Interfering with RNA Splicing.

INTRODUCTION: Retinitis pigmentosa (RP) is a rare degenerative retinal disease caused by mutations in approximately seventy genes. Currently, despite the availability of large-scale DNA sequencing technologies, ∼30-40% of patients still cannot be diagnosed at the molecular level. In this study, we investigated a novel intronic deletion of PDE6B, encoding the beta subunit of phosphodiesterase 6 in association with recessive RP. METHODS: Three unrelated consanguineous families were recruited from the northwestern part of Pakistan. Whole exome sequencing was performed for the proband of each family, and the data were analyzed according to an in-house computer pipeline. Relevant DNA variants in all available members of these families were assessed through Sanger sequencing. A minigene-based splicing assay was also performed. RESULTS: The clinical phenotype for all patients was compatible with rod cone degeneration, with the onset during childhood. Whole exome sequencing revealed a homozygous 18 bp intronic deletion (NM_000283.3:c.1921-20_1921-3del) in PDE6B, which co-segregated with disease in 10 affected individuals. In vitro splicing tests showed that this deletion causes aberrant RNA splicing of the gene, leading to the in-frame deletion of 6 codons and, likely, to disease. CONCLUSION: Our findings further expand the mutational spectrum of the PDE6B gene.

Pathogenic variants detected by RNA sequencing in Cornelia de Lange syndrome.

Recent studies suggest that transcript isoforms significantly overlap (approximately 60%) between brain tissue and Epstein-Barr virus-transformed lymphoblastoid cell lines (LCLs). Interestingly, 14 cohesion-related genes with variants that cause Cornelia de Lange Syndrome (CdLS) are highly expressed in the brain and LCLs. In this context, we first performed RNA sequencing of LCLs from 22 solved (with pathogenic variants) and 19 unsolved (with no confirmed variants) CdLS cases. Next, an RNA sequencing pipeline was developed using solved cases with two different methods: short variant analysis (for single-nucleotide and indel variants) and aberrant splicing detection analysis. Then, 19 unsolved cases were subsequently applied to our pipeline, and four pathogenic variants in NIPBL (one inframe deletion and three intronic variants) were newly identified. Two of three intronic variants were located at Alu elements in deep-intronic regions, creating cryptic exons. RNA sequencing with LCLs was useful for identifying hidden variants in exome-negative cases.

Short amplicon reverse transcription-polymerase chain reaction detects aberrant splicing in genes with low expression in blood missed by ribonucleic acid sequencing analysis for clinical diagnosis.

Use of blood RNA sequencing (RNA-seq) as a splicing analysis tool for clinical interpretation of variants of uncertain significance (VUSs) found via whole-genome and exome sequencing can be difficult for genes that have low expression in the blood due to insufficient read count coverage aligned to specific genes of interest. Here, we present a short amplicon reverse transcription-polymerase chain reaction(RT-PCR) for the detection of genes with low blood expression. Short amplicon RT-PCR, is designed to span three exons where an exon harboring a variant is flanked by one upstream and one downstream exon. We tested short amplicon RT-PCRs for genes that have median transcripts per million (TPM) values less than one according to the genotype-tissue expression database. Median TPM values of genes analyzed in this study are SYN1 = 0.8549, COL1A1 = 0.6275, TCF4 = 0.4009, DSP = .2894, TTN = 0.2851, COL5A2 = 0.1036, TERT = 0.04452, NTRK2 = 0.0344, ABCA4 = 0.00744, PRPH = 0, and WT1 = 0. All these genes show insufficient exon-spanning read coverage in our RNA-seq data to allow splicing analysis. We successfully detected all genes tested except PRPH and WT1. Aberrant splicing was detected in SYN1, TCF4, NTRK2, TTN, and TERT VUSs. Therefore, our results show short amplicon RT-PCR is a useful alternative for the analysis of splicing events in genes with low TPM in blood RNA for clinical diagnostics.

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.

Vulnerable exons, like ACADM exon 5, are highly dependent on maintaining a correct balance between splicing enhancers and silencers.

It is now widely accepted that aberrant splicing of constitutive exons is often caused by mutations affecting cis-acting splicing regulatory elements, but there is a misconception that all exons have an equal dependency on splicing regulatory elements and thus a similar susceptibility to aberrant splicing. We investigated exonic mutations in ACADM exon 5 to experimentally examine their effect on splicing and found that 7 out of 11 tested mutations affected exon inclusion, demonstrating that this constitutive exon is particularly vulnerable to exonic splicing mutations. Employing ACADM exon 5 and 6 as models, we demonstrate that the balance between splicing enhancers and silencers, flanking intron length, and flanking splice site strength are important factors that determine exon definition and splicing efficiency of the exon in question. Our study shows that two constitutive exons in ACADM have different inherent vulnerabilities to exonic splicing mutations. This suggests that in silico prediction of potential pathogenic effects on splicing from exonic mutations may be improved by also considering the inherent vulnerability of the exon. Moreover, we show that single nucleotide polymorphism that affect either of two different exonic splicing silencers, located far apart in exon 5, all protect against both immediately flanking and more distant exonic splicing mutations.

Identification of a novel TBX5 c.755 + 1 G > A variant and related pathogenesis in a family with Holt-Oram syndrome.

The proband with congenital heart disease and abnormal thumb was clinically diagnosed as Holt-Oram syndrome (HOS). A novel variant, T-box transcription factor 5 (TBX5) c.755 + 1 G > A, was identified in the proband via whole exome sequencing and validated using Sanger sequencing. Pedigree analysis and clinical examinations revealed three/seven individuals over three generations within the family, with features suggestive of HOS. Deep amplicon sequencing confirmed that the allele frequencies of the novel variant in the proband (III-1), her brother (III-2), and her mother (II-2) were 50%, 48.3%, and 38.1%, respectively, indicating that III-1 and III-2 harbored heterozygous variants, while II-2 harbored mosaic heterozygous variants. The minigene splicing assay showed that the novel variant affected the normal splicing of exon 7, resulting in the production of abnormal TBX5 transcripts. Reverse transcription-quantitative polymerase chain reaction and western blot analyses revealed that the novel variant upregulated TBX5 expression at the transcriptional and translational levels. Nuclear localization assay demonstrated impaired nuclear localization of the mutant TBX5. Cell viability assay revealed the inhibition of cell activity by the mutant TBX5. Our findings indicate that the novel variant was potentially induced HOS, probably by causing aberrant splicing, reducing the enrichment of nuclear TBX5 protein, and inhibiting cellular proliferation.