Deep Intronic Variant Research Articles

Reclassifying inconclusive diagnosis after newborn screening for cystic fibrosis. Moving forward

BackgroundNewborn screening for Cystic Fibrosis (CF) is associated with situations where the diagnosis of CF or CFTR related disorders (CFTR-RD) cannot be clearly ruled out. Materials/patients and methodsWe report a case series of 23 children with unconclusive diagnosis after newborn screening for CF and a mean follow-up of 7.7 years (4–13). Comprehensive investigations including whole CFTR gene sequencing, in vivo intestinal current measurement (ICM), nasal potential difference (NPD), and in vitro functional studies of variants of unknown significance, helped to reclassify the patients. ResultsExtensive genetic testing identified, in trans with a CF causing mutation, variants with varying clinical consequences and 3 variants of unknown significance (VUS). Eighteen deep intronic variants were identified by deep resequencing of the whole CFTR gene in 13 patients and were finally considered as non-pathogenic. All patients had normal CFTR dependent chloride transport in ICM. NPD differentiated 3 different profiles: CF-like tracings qualifying the patients as CF, such as F508del/D1152H patients; normal responses, suggesting an extremely low likelihood of developing a CFTR-RD such as F508del/TG11T5 patients; partial CFTR dysfunction above 20% of the normal, highlighting a remaining risk of developing CFTR-RD such as F508del/F1052V patients. The 3 VUS were reclassified as variant with defective maturation (D537N), defective expression (T582I) or with no clinical consequence (M952T). ConclusionThis study demonstrates the usefulness of combining genetic and functional investigations to assess the possibility of evolving to CF or CFTR-RD in babies with inconclusive diagnosis at neonatal screening.

Hereditary Multiple Exostoses-A Review of the Molecular Background, Diagnostics, and Potential Therapeutic Strategies.

Hereditary multiple exostoses (HMEs) syndrome, also known as multiple osteochondromas, represents a rare and severe human skeletal disorder. The disease is characterized by multiple benign cartilage-capped bony outgrowths, termed exostoses or osteochondromas, that locate most commonly in the juxta-epiphyseal portions of long bones. Affected individuals usually complain of persistent pain caused by the pressure on neighboring tissues, disturbance of blood circulation, or rarely by spinal cord compression. However, the most severe complication of this condition is malignant transformation into chondrosarcoma, occurring in up to 3.9% of HMEs patients. The disease results mainly from heterozygous loss-of-function alterations in the EXT1 or EXT2 genes, encoding Golgi-associated glycosyltransferases, responsible for heparan sulfate biosynthesis. Some of the patients with HMEs do not carry pathogenic variants in those genes, hence the presence of somatic mutations, deep intronic variants, or another genes/loci is suggested. This review presents the systematic analysis of current cellular and molecular concepts of HMEs along with clinical characteristics, clinical and molecular diagnostic methods, differential diagnosis, and potential treatment options.

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. Aspseudoexon 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.

Whole genome sequencing and in vitro splice assays reveal genetic causes for inherited retinal diseases

Inherited retinal diseases (IRDs) are a major cause of visual impairment. These clinically heterogeneous disorders are caused by pathogenic variants in more than 270 genes. As 30–40% of cases remain genetically unexplained following conventional genetic testing, we aimed to obtain a genetic diagnosis in an IRD cohort in which the genetic cause was not found using whole-exome sequencing or targeted capture sequencing. We performed whole-genome sequencing (WGS) to identify causative variants in 100 unresolved cases. After initial prioritization, we performed an in-depth interrogation of all noncoding and structural variants in genes when one candidate variant was detected. In addition, functional analysis of putative splice-altering variants was performed using in vitro splice assays. We identified the genetic cause of the disease in 24 patients. Causative coding variants were observed in genes such as ATXN7, CEP78, EYS, FAM161A, and HGSNAT. Gene disrupting structural variants were also detected in ATXN7, PRPF31, and RPGRIP1. In 14 monoallelic cases, we prioritized candidate noncanonical splice sites or deep-intronic variants that were predicted to disrupt the splicing process based on in silico analyses. Of these, seven cases were resolved as they carried pathogenic splice defects. WGS is a powerful tool to identify causative variants residing outside coding regions or heterozygous structural variants. This approach was most efficient in cases with a distinct clinical diagnosis. In addition, in vitro splice assays provide important evidence of the pathogenicity of rare variants.

Abstract 14385: Genome Sequencing and Re-Engineered Heart Cell-Disease Modeling Identifies a Novel Deep Intronic Variant in a Multi-Generational Long QT Syndrome Pedigree

Introduction: Most of long QT syndrome (LQTS) is explained by pathogenic variants in KCNQ1 , KCNH2 , or SCN5A . However, ~10-20% of LQTS index cases remain genetically elusive following commercial genetic testing. Here, we identified and functionally characterized a novel LQTS genetic substrate in a multi-generational, “genotype negative” LQTS pedigree. Methods: The index case was a 40-year-old female with a history of syncope, seizures, ventricular fibrillation, and a family history of LQTS and sudden death. Although anticipated to have KCNH2- mediated LQT2, commercial genetic testing of all LQTS-causative genes was negative. Genome sequencing and variant analysis was performed on 6 affected family members. Patient-specific and CRISPR/Cas9 “gene-corrected” isogenic control induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) were generated. RT-PCR was used to determine the effects of the identified variant on splicing. Action potential duration (APD) and field potential duration (FPD) were measured using voltage-sensing dye (Fluovolt) and microelectrode array, respectively. Results: No ultra-rare, nonsynonymous heterozygous variants co-segregated among the 6 LQTS phenotype-positive individuals. Instead, a deep intronic KCNH2 variant (c.3331-316 G>T) was present in all affected individuals. RT-PCR analysis of patient-specific iPSC-CM-derived RNA revealed that c.3331-316 G>T creates a novel 89 base-pair exon residing between canonical exons 14 and 15 of KCNH2 that results in a frame-shift variant annotated as p.S1112Pfs*171. The APD 90 was significantly longer in p.S1112Pfs*171-iPSC-CMs (647 ± 12 ms, n=50) compared to isogenic control iPSC-CMs (412 ± 10 ms, n=50, p < 0.0001). Further, the FPD was significantly longer in p.S1112Pfs*171-iPSC-CMs (349 ± 8 ms, n=48) compared to isogenic control iPSC-CMs (209 ± 10 ms, n=20, p < 0.0001). Conclusions: A novel deep intronic KCNH2 variant was identified in a multi-generational genetically elusive LQTS pedigree. The patient-derived and isogenic control iPSC-CMs establish that the p.S1112Pfs*171-KCNH2 variant is the monogenetic cause for this family’s LQTS. As suggested by the phenotype, this family indeed has KCNH2 -mediated LQT2 albeit stemming from a deep intronic variant .

Identification of novel deep intronic PAH gene variants in patients diagnosed with phenylketonuria.

Phenylketonuria (PKU) is caused by phenylalanine hydroxylase (PAH) gene variants. Previously, 94.21% of variants were identified using Sanger sequencing and multiplex ligation-dependent probe amplification. To investigate the remaining variants, we performed whole-genome sequencing for four patients with PKU and unknown genotypes to identify deep intronic or structural variants. We identified three novel heterozygous variants (c.706+368T>C, c.1065+241C>A, and c.1199+502A>T) in a deep PAH gene intron. We detected a c.1199+502A>T variant in 60% (6/10) of PKU patients with genetically undetermined PKU. In silico predictions indicated that the three deep variants may impact splice site selection and result in the inclusion of a pseudo-exon. A c.1199+502A>T PAH minigene and reverse transcription PCR (RT-PCR) on blood RNA from a PKU patient with biallelic variants c.1199+502A>T and c.1199G>A confirmed that the c.1199+502A>T variant may strengthen the predicted branch point and leads to the inclusion of a 25-nt pseudo-exon in the PAH mRNA. Reverse transcription polymerase chain reaction (RT-PCR) on the minigene revealed that c.706+368T>C may create an SRSF2 (SC35) binding site via a 313-nt pseudo-exon, whereas c.1065+241C>A may produce an 81-nt pseudo-exon that strengthens the predicted SRSF1 (SF2/ASF) binding site. These results augment current knowledge of PAH genotypes and show that deep intronic analysis of PAH can genetically diagnose PKU.

Detection of single nucleotide and copy number variants in the Fabry disease-associated GLA gene using nanopore sequencing

More than 900 variants have been described in the GLA gene. Some intronic variants and copy number variants in GLA can cause Fabry disease but will not be detected by classical Sanger sequence. We aimed to design and validate a method for sequencing the GLA gene using long-read Oxford Nanopore sequencing technology. Twelve Fabry patients were blindly analyzed, both by conventional Sanger sequence and by long-read sequencing of a 13 kb PCR amplicon. We used minimap2 to align the long-read data and Nanopolish and Sniffles to call variants. All the variants detected by Sanger (including a deep intronic variant) were also detected by long-read sequencing. One patient had a deletion that was not detected by Sanger sequencing but was detected by the new technology. Our long-read sequencing-based method was able to detect missense variants and an exonic deletion, with the added advantage of intronic analysis. It can be used as an efficient and cost-effective tool for screening and diagnosing Fabry disease.

Targeted massively parallel sequencing of candidate regions on chromosome 22q predisposing to multiple schwannomas: An analysis of 51 individuals in a single-center experience.

Constitutional LZTR1 or SMARCB1 pathogenic variants (PVs) have been found in ∼86% of familial and ∼40% of sporadic schwannomatosis cases. Hence, we performed massively parallel sequencing of the entire LZTR1, SMARCB1, and NF2 genomic loci in 35 individuals with schwannomas negative for constitutional first-hit PVs in the LZTR1/SMARCB1/NF2 coding sequences; however, with 22q deletion and/or a different NF2 PV in each tumor, including six cases with only one tumor available. Furthermore, we verified whether any other LZTR1/SMARCB1/NF2 (likely) PVscould be found in 16 cases carrying a SMARCB1 constitutional variant in the 3'-untranslated region (3'-UTR) c.*17C>T, c.*70C>T, or c.*82C>T. As no additional variants were found, functional studies were performed to clarify the effect of these 3'-UTR variants on the transcript. The 3'-UTR variants c.*17C>T and c.*82C>T showed pathogenicity by negatively affecting theSMARCB1 transcript level. Two novel deep intronic SMARCB1 variants, c.500+883T>G and c.500+887G>A, resulting in out-of-frame missplicing of intron 4, were identified in two unrelated individuals. Further resequencing of the entire repeat-masked genomics sequences of chromosome 22q in individuals negative for PVs in the SMARCB1/LZTR1/NF2 coding- and noncoding regions revealed five potential schwannomatosis-predisposing candidate genes, that is,MYO18B, NEFH, SGSM1, SGSM3, and SBF1, pending further verification.

Pathogenic neurofibromatosis type 1 (NF1) RNA splicing resolved by targeted RNAseq

Neurofibromatosis type 1 (NF1) is caused by loss-of-function variants in the NF1 gene. Approximately 10% of these variants affect RNA splicing and are either missed by conventional DNA diagnostics or are misinterpreted by in silico splicing predictions. Therefore, a targeted RNAseq-based approach was designed to detect pathogenic RNA splicing and associated pathogenic DNA variants. For this method RNA was extracted from lymphocytes, followed by targeted RNAseq. Next, an in-house developed tool (QURNAs) was used to calculate the enrichment score (ERS) for each splicing event. This method was thoroughly tested using two different patient cohorts with known pathogenic splice-variants in NF1. In both cohorts all 56 normal reference transcript exon splice junctions, 24 previously described and 45 novel non-reference splicing events were detected. Additionally, all expected pathogenic splice-variants were detected. Eleven patients with NF1 symptoms were subsequently tested, three of which have a known NF1 DNA variant with a putative effect on RNA splicing. This effect could be confirmed for all 3. The other eight patients were previously without any molecular confirmation of their NF1-diagnosis. A deep-intronic pathogenic splice variant could now be identified for two of them (25%). These results suggest that targeted RNAseq can be successfully used to detect pathogenic RNA splicing variants in NF1.

A novel deep intronic variant strongly associates with Alkaptonuria

Alkaptonuria is a rare autosomal recessive inherited disorder of tyrosine metabolism, which causes ochronosis, arthropathy, cardiac valvular calcification, and urolithiasis. The epidemiology of alkaptonuria in East Asia is not clear. In this study, patients diagnosed with alkaptonuria from January 2010 to June 2020 were reviewed. Their clinical and molecular features were further compared with those of patients from other countries. Three patients were found to have alkaptonuria. Mutation analyses of the homogentisate 1,2-dioxygenase gene (HGD) showed four novel variants c.16-2063 A > C, p.(Thr196Ile), p.(Gly344AspfsTer25), and p.(Gly362Arg) in six mutated alleles (83.3%). RNA sequencing revealed that c.16-2063 A > C activates a cryptic exon, causing protein truncation p.(Tyr5_Ile6insValTer17). A literature search identified another 6 patients with alkaptonuria in East Asia; including our cases, 13 of the 18 mutated alleles have not been reported elsewhere in the world. Alkaptonuria is rare in Taiwan and East Asia, with HGD variants being mostly novel and private.

657: Importance of screening for splice site and deep intronic variants—Insights from cystic fibrosis genetic testing on patients from the Indian subcontinent

657: Importance of screening for splice site and deep intronic variants—Insights from cystic fibrosis genetic testing on patients from the Indian subcontinent

COLLAGEN RELATED MUSCLE DISEASES: EP.57 Identifying preclinical outcome measures of a novel humanized knock-in mouse modeling a collagen VI deep intronic pathogenic variant

COLLAGEN RELATED MUSCLE DISEASES: EP.57 Identifying preclinical outcome measures of a novel humanized knock-in mouse modeling a collagen VI deep intronic pathogenic variant

MRNA analysis identifies deep intronic variants causing Alport syndrome and overcomes the problem of negative results of exome sequencing

Mutations in COL4A3, COL4A4 and COL4A5 genes lead to Alport syndrome (AS). However, pathogenic variants in some AS patients are not detected by exome sequencing. The aim of this study was to identify the underlying genetic causes of five unrelated AS probands with negative next-generation sequencing (NGS) test results. Urine COL4A3–5 mRNAs were analyzed in the probands with an uncertain inherited mode of AS, and COL4A5 mRNA of skin fibroblasts was analyzed in the probands with X-linked AS. RT-PCR and direct sequencing were performed to detect mRNA abnormalities. PCR and direct sequencing were used to analyze the exons with flanking intronic sequences corresponding to mRNA abnormalities. Six novel deep intronic splicing variants in COL4A4 and COL4A5 genes that cannot be captured by exome sequencing were identified in the four AS probands. Skipping of an exon was caused by an intronic variant, and retention of an intron fragment caused by five variants. In the remaining AS proband, COL4A5 variants c.2677 + 646 C > T and r.2678_r.2767del were detected at the DNA and RNA level, respectively, whereas it is unclear whether c.2677 + 646 C > T may not lead to r.2678_r.2767del. Our results reveal that mRNA analysis for AS genes from either urine or skin fibroblasts can resolve genetic diagnosis in AS patients with negative NGS results. We recommend analyzing COL4A3–5 mRNA from urine as the first choice for these patients because it is feasible and non-invasive.

Deep intronic variant in the ARSB gene as the genetic cause for Maroteaux-Lamy syndrome (MPS VI).

Maroteaux-Lamy syndrome (MPS-VI) is a rare autosomal-recessive disorder with a wide spectrum of clinical manifestations, ranging from an attenuated to a rapidly progressive disease. It is caused by variants in ARSB, which encodes the lysosomal arylsulfatase B (ARSB) enzyme, part of the degradation process of glycosaminoglycans in lysosomes. Over 220 variants have been reported so far, with a majority of missense variants. We hereby report two siblings of Bedouin origin with a diagnosis of MPS-VI. Western blots in patient fibroblasts revealed total absence of ARSB protein production. Complete sequencing of the coding region of ARSB did not identify a candidate disease-associated variant. However, deep sequencing of the noncoding region of ARSB by whole genome sequencing (WGS) revealed a c.1142+581A to G variant. The variant is located within intron 5 and fully segregated with the disease in the family. Determination of the genetic cause for these patients enabled targeted treatment by enzyme replacement therapy, along with appropriate genetic counseling and prenatal diagnosis for the family. These results highlight the advantage of WGS as a powerful tool, for improving the diagnostic rate of rare disease-causing variants, and emphasize the importance of studying deep intronic sequence variation as a cause of monogenic disorders.

High-Throughput Sequencing to Identify Mutations Associated with Retinal Dystrophies.

Retinal dystrophies (RD) are clinically and genetically heterogenous disorders showing mutations in over 270 disease-associated genes. Several millions of people worldwide are affected with different types of RD. Studying the relevance of disease-associated sequence alterations will assist in understanding disorders and may lead to the development of therapeutic approaches. Here, we established a whole exome sequencing (WES) pipeline to rapidly identify disease-associated mutations in patients. Sanger sequencing was applied to identify deep-intronic variants and to verify the co-segregation of WES results within families. We analyzed 26 unrelated patients with different syndromic and non-syndromic clinical manifestations of RD. All patients underwent ophthalmic examinations. We identified nine novel disease-associated sequence variants among 37 variants identified in total. The sequence variants located to 17 different genes. Interestingly, two cases presenting with Stargardt disease carried deep-intronic variants in ABCA4. We have classified 21 variants as pathogenic variants, 4 as benign/likely benign variants, and 12 as variants of uncertain significance. This study highlights the importance of WES-based mutation analyses in RD patients supporting clinical decisions, broadly based genetic diagnosis and support genetic counselling. It is essential for any genetic therapy to expand the mutation spectrum, understand the genes’ function, and correlate phenotypes with genotypes.

Exploring the missing heritability in subjects with hearing loss, enlarged vestibular aqueducts, and a single or no pathogenic SLC26A4 variant

Pathogenic variants in SLC26A4 have been associated with autosomal recessive hearing loss (arHL) and a unilateral or bilateral enlarged vestibular aqueduct (EVA). SLC26A4 is the second most frequently mutated gene in arHL. Despite the strong genotype–phenotype correlation, a significant part of cases remains genetically unresolved. In this study, we investigated a cohort of 28 Dutch index cases diagnosed with HL in combination with an EVA but without (M0) or with a single (M1) pathogenic variant in SLC26A4. To explore the missing heritability, we first determined the presence of the previously described EVA-associated haplotype (Caucasian EVA (CEVA)), characterized by 12 single nucleotide variants located upstream of SLC26A4. We found this haplotype and a delimited V1-CEVA haplotype to be significantly enriched in our M1 patient cohort (10/16 cases). The CEVA haplotype was also present in two M0 cases (2/12). Short- and long-read whole genome sequencing and optical genome mapping could not prioritize any of the variants present within the CEVA haplotype as the likely pathogenic defect. Short-read whole-genome sequencing of the six M1 cases without this haplotype and the two M0/CEVA cases only revealed previously overlooked or misinterpreted splice-altering SLC26A4 variants in two cases, who are now genetically explained. No deep-intronic or structural variants were identified in any of the M1 subjects. With this study, we have provided important insights that will pave the way for elucidating the missing heritability in M0 and M1 SLC26A4 cases. For pinpointing the pathogenic effect of the CEVA haplotype, additional analyses are required addressing defect(s) at the RNA, protein, or epigenetic level.

An Overview of the Genetics of ABCA4 Retinopathies, an Evolving Story.

Stargardt disease (STGD1) and ABCA4 retinopathies (ABCA4R) are caused by pathogenic variants in the ABCA4 gene inherited in an autosomal recessive manner. The gene encodes an importer flippase protein that prevents the build-up of vitamin A derivatives that are toxic to the RPE. Diagnosing ABCA4R is complex due to its phenotypic variability and the presence of other inherited retinal dystrophy phenocopies. ABCA4 is a large gene, comprising 50 exons; to date > 2000 variants have been described. These include missense, nonsense, splicing, structural, and deep intronic variants. Missense variants account for the majority of variants in ABCA4. However, in a significant proportion of patients with an ABCA4R phenotype, a second variant in ABCA4 is not identified. This could be due to the presence of yet unknown variants, or hypomorphic alleles being incorrectly classified as benign, or the possibility that the disease is caused by a variant in another gene. This underlines the importance of accurate genetic testing. The pathogenicity of novel variants can be predicted using in silico programs, but these rely on databases that are not ethnically diverse, thus highlighting the need for studies in differing populations. Functional studies in vitro are useful towards assessing protein function but do not directly measure the flippase activity. Obtaining an accurate molecular diagnosis is becoming increasingly more important as targeted therapeutic options become available; these include pharmacological, gene-based, and cell replacement-based therapies. The aim of this review is to provide an update on the current status of genotyping in ABCA4 and the status of the therapeutic approaches being investigated.

Growth Hormone Receptor (Ghr) 6ω Pseudoexon Activation: A Novel Cause Of Severe Growth Hormone Insensitivity (Ghi).

ContextSevere forms of growth hormone insensitivity (GHI) are characterized by extreme short stature, dysmorphism, and metabolic anomalies.ObjectiveThis work aims to identify the genetic cause of growth failure in 3 “classical” GHI individuals.MethodsA novel intronic growth hormone receptor gene (GHR) variant was identified, and in vitro splicing assays confirmed aberrant splicing. A 6Ω pseudoexon GHR vector and patient fibroblast analysis assessed the consequences of the novel pseudoexon inclusion and the impact on GHR function.ResultsWe identified a novel homozygous intronic GHR variant (g.5:42700940T > G, c.618+836T > G), 44 bp downstream of the previously recognized intronic 6Ψ GHR pseudoexon mutation in the index patient. Two siblings also harbored the novel intronic 6Ω pseudoexon GHR variant in compound heterozygosity with the known GHR c.181C > T (R43X) mutation. In vitro splicing analysis confirmed inclusion of a 151-bp mutant 6Ω pseudoexon not identified in wild-type constructs. Inclusion of the 6Ω pseudoexon causes a frameshift resulting in a nonfunctional truncated GHR lacking the transmembrane and intracellular domains. The truncated 6Ω pseudoexon protein demonstrated extracellular accumulation and diminished activation of STAT5B signaling following GH stimulation.ConclusionNovel GHR 6Ω pseudoexon inclusion results in loss of GHR function consistent with a severe GHI phenotype. This represents a novel mechanism of Laron syndrome and is the first deep intronic variant identified causing severe postnatal growth failure. The 2 kindreds originate from the same town in Campania, Southern Italy, implying common ancestry. Our findings highlight the importance of studying variation in deep intronic regions as a cause of monogenic disorders.

Updating the Genetic Landscape of Inherited Retinal Dystrophies.

Inherited retinal dystrophies (IRD) are a group of diseases characterized by the loss or dysfunction of photoreceptors and a high genetic and clinical heterogeneity. Currently, over 270 genes have been associated with IRD which makes genetic diagnosis very difficult. The recent advent of next generation sequencing has greatly facilitated the diagnostic process, enabling to provide the patients with accurate genetic counseling in some cases. We studied 92 patients who were clinically diagnosed with IRD with two different custom panels. In total, we resolved 53 patients (57.6%); in 12 patients (13%), we found only one mutation in a gene with a known autosomal recessive pattern of inheritance; and 27 patients (29.3%) remained unsolved. We identified 120 pathogenic or likely pathogenic variants; 30 of them were novel. Among the cone-rod dystrophy patients, ABCA4 was the most common mutated gene, meanwhile, USH2A was the most prevalent among the retinitis pigmentosa patients. Interestingly, 10 families carried pathogenic variants in more than one IRD gene, and we identified two deep-intronic variants previously described as pathogenic in ABCA4 and CEP290. In conclusion, the IRD study through custom panel sequencing demonstrates its efficacy for genetic diagnosis, as well as the importance of including deep-intronic regions in their design. This genetic diagnosis will allow patients to make accurate reproductive decisions, enroll in gene-based clinical trials, and benefit from future gene-based treatments.

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 SummaryThere 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.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.