Identification of Deep-Intronic Splice Mutations in a Large Cohort of Patients With Inherited Retinal Diseases

Clinically relevant pseudoexons of the GALNS gene and their antisense-based correction

BACKGROUND

Biallelic pathogenic variants in the GALNS gene lead to Mucopolysaccharidosis Type IVA (MPS IVA), a rare lysosomal storage disorder. GALNS encodes the enzyme N-acetylgalactosamine-6-sulfatase, whose deficiency causes accumulation of glycosaminoglycans and leads to a broad spectrum of clinical manifestations primarily affecting the osteoarticular system. Several studies have shown that, in 10%-15% of patients with the biochemical phenotype of MPS IVA, standard molecular genetic testing fails to identify one or both causative variants in the GALNS gene.

METHODS

We performed an in-depth investigation of GALNS' splicing, with a special focus on deep-intronic mutations that lead to activation of pseudoexons (PEs). Using bioinformatic tools, we analyzed all deep-intronic variants in GALNS available in public databases and subjected the most relevant ones to in vitro analyses using minigenes.

RESULTS

We characterized eight PE-activating variants, one of which (c.121-210C > T) represents a recurrent pathogenic variant which has long been hidden behind the mask of a polymorphic variant. In addition, we demonstrate that GALNS' splicing can produce a diverse range of mRNA isoforms containing so-called wild-type PEs, which are present at low levels as part of non-productive splicing, and weak canonical exons which are prone to skipping. We show that PE-activating variants cluster within wild-type PEs, highlighting the need for closer scrutiny of these regions during genetic testing. Finally, we applied modified U7 small nuclear RNAs and circular RNAs to efficiently block the identified PEs and pave the way for personalized antisense-based therapy for MPS IVA patients.

CONCLUSION

The results of this study expand the understanding of GALNS gene splicing, indicating hotspots for splicing mutations. The presented data not only help to increase the diagnostic yield for MPS IVA but also unveil new therapeutic approaches for a number of MPS IVA patients.

Targeted long-read cDNA sequencing reveals novel splice-altering pathogenic variants causing retinal dystrophies

Splice-altering variants are suggested to be responsible for part of the missing heritability of inherited retinal diseases (IRDs). The interpretation of these variants is challenging as functional evidence is required to validate pathogenicity. We explored the diagnostic value of a targeted long-read cDNA sequencing (lrcDNA-seq) approach to investigate IRD-associated splicing defects. For each affected individual, RNA was isolated from blood, and for each candidate gene, cDNA amplicons, spanning the complete open reading frame or multiple exons, were generated and subjected to long-read sequencing. We validated our approach by assessing previously described pathogenic splice-altering variants in IRD-associated genes. Next, we investigated six genetically unexplained affected individuals, each carrying pathogenic variant(s) in NMNAT1. In two probands, we provided functional validation for previously identified variants of uncertain significance present on the second allele. In four other subjects, lrcDNA-seq revealed the partial inclusion of an SVA_F retrotransposon in the NMNAT1 mRNA, predicted to introduce a premature stop codon. We showed that targeted lrcDNA-seq is effective in characterizing splice defects and in identifying novel splice-altering variants and uncovered the IRD genetic basis for six previously unexplained subjects. We believe that the implementation of this technique has the potential to contribute to an increased diagnostic rate of IRDs.

Non-coding single-nucleotide and structural variants affecting the EYS putative promoter cause autosomal recessive retinitis pigmentosa

PURPOSE

Variants in untranslated genomic regions are difficult to identify as pathogenic but are capable of causing disease by interfering with gene expression. This study aimed to characterize the effect of variants identified in the 5'-untranslated region of EYS in patients with autosomal recessive retinitis pigmentosa (RP).

METHODS

Variant screening included gene panels, Sanger, exome, and genome sequencing. Functional validation included an electrophoretic mobility shift assay and various luciferase assays.

RESULTS

Patients with RP from 6 EYS biallelic Arab-Muslim families harbored a 5' noncoding EYS variant, c.-453G>T, and 4 harbored a structural variant affecting the 5' noncoding exons. Electrophoretic mobility shift assay analysis revealed an effect on binding of transcription factors for c.-453G>T and a neighboring variant c.-454G>T. Dual luciferase assays using overexpression of various transcription factors showed distinct effects on expression. c.-453G>T was associated with higher luciferase expression with CRX overexpression and c.-454G>C with OTX2 overexpression. In addition, the 2 variants were found to influence translation by affecting upstream initiation codons. Interestingly, visual function of EYS RP patients who harbor c.-453G>T are better than those with biallelic null EYS variants.

CONCLUSION

Our analysis revealed both single-nucleotide and structural variants in the EYS promoter as the cause of autosomal recessive RP. These variants may affect EYS expression via a dual mechanism by altering transcription factor binding affinity at the EYS promoter and by affecting upstream open reading frames.

OpenSpliceAI: An efficient, modular implementation of SpliceAI enabling easy retraining on non-human species

The SpliceAI deep learning system is currently one of the most accurate methods for identifying splicing signals directly from DNA sequences. However, its utility is limited by its reliance on older software frameworks and human-centric training data. Here we introduce OpenSpliceAI, a trainable, open-source version of SpliceAI implemented in PyTorch to address these challenges. OpenSpliceAI supports both training from scratch and transfer learning, enabling seamless re-training on species-specific datasets and mitigating human-centric biases. Our experiments show that it achieves faster processing speeds and lower memory usage than the original SpliceAI code, allowing large-scale analyses of extensive genomic regions on a single GPU. Additionally, OpenSpliceAI's flexible architecture makes for easier integration with established machine learning ecosystems, simplifying the development of custom splicing models for different species and applications. We demonstrate that OpenSpliceAI's output is highly concordant with SpliceAI. In silico mutagenesis (ISM) analyses confirm that both models rely on similar sequence features, and calibration experiments demonstrate similar score probability estimates.

Developing splice-switching oligonucleotides for urea cycle disorder using an integrated diagnostic and therapeutic platform

BACKGROUNDS & AIMS

Citrin deficiency (CD) is an autosomal recessive urea cycle disorder caused by biallelic loss-of-function variants in the SLC25A13 gene, leading to life-threatening hyperammonemia and hypoglycemia. Variants in deep introns can cause genetic diseases by altering splicing and are often missed by current diagnostic tools. Splice-switching oligonucleotides (SSOs) can resolve certain intronic variants, but patients harboring such variants need to be identified. We present a lean workflow from molecular diagnostics to SSO development to resolve splice-altering variants in deep introns that is applicable to other genetic disorders.

METHODS

A deep intronic-gene panel was designed to identify deep intronic variants. SSOs were then developed and validated in vitro using a minigene assay and induced hepatocytes, and target engagement was verified in vivo by hydrodynamic tail vein injection of minigenes and SSOs.

RESULTS

With the deep intronic-gene panel and RNA analysis, we identified a novel SLC25A13 c.469-2922G>T variant that promotes the inclusion of a premature stop codon-containing pseudo-exon, SLC25A13-PE5, thereby causing CD. By a stepwise rational SSO design approach, we identified potent candidates inhibiting SLC25A13-PE5 at EC50 <2 nM in vitro. Upon conjugating the SSOs with GalNAc (N-acetylgalactosamine), they were validated to rescue normal protein expression and restore ureagenesis and ammonia clearance, key urea cycle functions, in patient-derived induced hepatocytes. In vivo on-target efficacy of the clinical GalNAc-SSO candidate, in the absence of acute toxicity and inflammation, was observed in a mouse model with exogenous hepatic minigene expression.

CONCLUSIONS

Our data validates a platform to redefine the molecular diagnosis of urea cycle disorders and provides proof-of-concept for a precision therapy for patients with CD, for whom the only effective treatment is liver transplantation.

IMPACT AND IMPLICATIONS

Deep intronic variants are common causes of genetic diseases that are commonly neglected. In this study, we demonstrate an integrated precision diagnostic and therapeutic approach for urea cycle disorders. Specifically, we focus on citrin deficiency, going from the discovery of a novel splice variant in the SLC25A13 gene with our novel deep intronic-gene panel for urea cycle disorders, to the development and in vivo validation of an efficacious splice-switching oligonucleotide candidate for the pathogenic splice variant. We envision the possibility of extrapolating this pipeline to the diagnosis and development of treatments for other rare genetic diseases.

Loss of MEF2C function by enhancer mutation leads to neuronal mitochondria dysfunction and motor deficits in mice

BACKGROUND

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder characterized by the loss of both upper and lower motor neurons, leading to progressive paralysis. Both genetic alterations and epigenetic modifications contribute to neuronal dysfunction in the pathogenesis of ALS. However, the mechanism behind genetic mutations in the non-coding region of genes that affect epigenetic modifications remains unclear.

METHODS

Convolutional neural network was used to identify an ALS-associated SNP located in the intronic region of MEF2C (rs304152), residing in a putative enhancer element. To examine the alteration of MEF2C transcription by the SNP, we generated HEK293T cells carrying the major or minor allele by CRISPR-Cas9. To verify the role of MEF2C-knockdown (MEF2C-KD) in mice, we developed AAV expressing shRNA for MEF2C based on AAV-U6 promoter vector. Neuropathological alterations of MEF2C-KD mice with mitochondrial dysfunction and motor neuronal damage were observed by confocal microscopy and transmission electron microscope (TEM). Behavioral changes of mice were examined through longitudinal study by tail suspension, inverted grid test and automated gait analysis.

RESULTS

Here, we show that enhancer mutation of MEF2C reduces own gene expression and consequently impairs mitochondrial function in motor neurons. MEF2C localizes and binds to the mitochondria DNA, and directly modulates mitochondria-encoded gene expression. CRISPR/Cas-9-induced mutation of the MEF2C enhancer decreases expression of mitochondria-encoded genes. Moreover, MEF2C mutant cells show reduction of mitochondrial membrane potential, ATP level but elevation of oxidative stress. MEF2C deficiency in the upper and lower motor neurons of mice impairs mitochondria-encoded genes, and leads to mitochondrial metabolic disruption and progressive motor behavioral deficits.

CONCLUSIONS

Together, MEF2C dysregulation by the enhancer mutation leads to mitochondrial dysfunction and oxidative stress, which are prevalent features in motor neuronal damage and ALS pathogenesis. This genetic and epigenetic crosstalk mechanism provides insights for advancing our understanding of motor neuron disease and developing effective treatments.

Critical Considerations in Calling Disease-Causing EDAR Mutations in Nonsyndromic Oligodontia

Background/Objectives: Oligodontia, the absence of six or more teeth excluding third molars, is a rare genetic condition, unlike hypodontia (missing one or more teeth), which is a relatively common human disease. Methods: To identify the genetic etiology of nonsyndromic oligodontia (NSO) families, we performed mutational analysis and investigated the functional effects of identified EDAR mutations. Whole-exome sequencing was conducted on recruited families with NSO. Bioinformatic analysis identified mutations in oligodontia-causing candidate genes, which were confirmed by Sanger sequencing and segregation within families. The impact of EDAR mutations on the EDA signaling pathway was assessed using luciferase activity analysis. Results:EDAR mutations were identified in three NSO families. A homozygous missense EDAR mutation (NM_022336.4: c.319A>G p.(Met107Val)) was found in the singleton proband of family 1. The proband of family 2 carried compound heterozygous EDAR mutations: a maternal missense mutation (c.319A>G p.(Met107Val)) and a paternal missense variant (c.1270G>A p.(Val424Met)). The proband of family 3 had heterozygous EDAR mutations: a maternal missense mutation (c.389T>A p.(Ile130Asn)) and paternal intronic variants in cis (c.[357-4G>A;440+50C>T]). Luciferase assays confirmed reduced transcriptional activity for all identified missense mutations, while splicing assays revealed altered splicing patterns. Conclusions: In conclusion, recessive EDAR mutations, including novel ones, were identified in NSO families, and their pathological mechanism was explored through transcriptional activity measurements.

Novel compound heterozygous variants in LTBP2 associated with relative anterior microphthalmos

PURPOSE

Relative anterior microphthalmos (RAM) is a rare congenital defect associated with severe vision impairment that is primarily caused by genetic alterations. The purpose of this study was to identify the causative genetic variants in two Chinese families with RAM with an autosomal recessive inheritance pattern.

METHODS

DNA samples were obtained from two probands and their family members. Targeted next-generation sequencing (NGS) was used to screen 425 genes associated with inherited eye diseases to identify possible disease-causing variants in the two patients. Sanger sequencing was subsequently used to validate the results in both families.

RESULTS

The targeted NGS panel identified potentially causative novel variants of the latent transforming growth factor beta binding protein 2 (LTBP2) gene in the two RAM families: a missense variant (c.2771C > T; p.Ala924Val) and an intronic variant (c.4582 + 9A > G) in Family A and a different missense variant (c.5239C > A; p.Arg1747Ser) and a synonymous variant (c.951G > A; p.Pro317Pro) in Family B. These four novel variants all cosegregated with the disease phenotype.

CONCLUSION

To our knowledge, this is the first study to report novel LTBP2 gene variants related to RAM. Considering the importance of LTBP2 in ocular development, we provide initial insights into the potential pathogenic mechanisms of LTBP2 in RAM.

18-Years of single-centre DNA testing in over 7000 index cases with inherited retinal dystrophies and optic neuropathies

Inherited retinal dystrophies (IRDs) and inherited optic neuropathies (IONs) are characterized by distinct genetic causes and molecular mechanisms that can lead to varying degrees of visual impairment. The discovery of pathogenic variants in numerous genes associated with these conditions has deepened our understanding of the molecular pathways that influence both vision and disease manifestation and may ultimately lead to novel therapeutic approaches. Over the past 18 years, our DNA diagnostics unit has been performing genetic testing on patients suspected of having IRD or ION, using state-of-the-art mutation detection technologies that are continuously updated. This report presents a retrospective analysis of genetic data from 6237 IRD and 780 ION patients. Out of these, 3054 IRD patients (49.0%) and 211 ION patients (27.1%) received a definitive molecular diagnosis, with disease-causing variants identified in 139 different genes. The genes most implicated in disease pathologies are ABCA4, accounting for 23.8% of all IRD/ION index cases, followed by BEST1 (7.8%), USH2A (6.2%), PRPH2 (5.7%), RPGR (5.6%), RS1 (5.5%), OPA1 (4.3%), and RHO (3.1%). Our study has compiled the most extensive dataset in combined IRD/ION diagnostics to date and offers valuable insights into the frequencies of mutant alleles and the efficiency of mutation detection in various inherited retinal conditions.

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.

Investigating Splice Defects in USH2A Using Targeted Long-Read Sequencing

Biallelic variants in USH2A are associated with retinitis pigmentosa (RP) and Type 2 Usher Syndrome (USH2), leading to impaired vision and, additionally, hearing loss in the latter. Although the introduction of next-generation sequencing into clinical diagnostics has led to a significant uplift in molecular diagnostic rates, many patients remain molecularly unsolved. It is thought that non-coding variants or variants of uncertain significance contribute significantly to this diagnostic gap. This study aims to demonstrate the clinical utility of the reverse transcription-polymerase chain reaction (RT-PCR)-Oxford Nanopore Technology (ONT) sequencing of USH2A mRNA transcripts from nasal epithelial cells to determine the splice-altering effect of candidate variants. Five affected individuals with USH2 or non-syndromic RP who had undergone whole genome sequencing were recruited for further investigation. All individuals had uncertain genotypes in USH2A, including deep intronic rare variants, c.8682-654C>G, c.9055+389G>A, and c.9959-2971C>T; a synonymous variant of uncertain significance, c.2139C>T; p.(Gly713=); and a predicted loss of function duplication spanning an intron/exon boundary, c.3812-3_3837dup p.(Met1280Ter). In silico assessment using SpliceAI provided splice-altering predictions for all candidate variants which were investigated using ONT sequencing. All predictions were found to be accurate; however, in the case of c.3812-3_3837dup, the outcome was a complex cryptic splicing pattern with predominant in-frame exon 18 skipping and a low level of exon 18 inclusion leading to the predicted stop gain. This study detected and functionally characterised simple and complex mis-splicing patterns in USH2A arising from previously unknown deep intronic variants and previously reported variants of uncertain significance, confirming the pathogenicity of the variants.

NGS data analysis for molecular diagnosis of Inborn Errors of Immunity

Inborn errors of immunity (IEI) encompass a group of disorders with a strong genetic component. Prompt and accurate diagnosis of these disorders is essential for effective clinical management. Next-generation sequencing (NGS) has significantly enhanced the diagnostic process by offering a comprehensive and scalable approach for identifying genomic variations causal for these disorders. Nevertheless, the bioinformatics analysis of NGS data poses several challenges. In this review, we explore these challenges and share our insights on addressing them, aiming to improve the overall diagnostic yield.

A deep intronic variant in MME causes autosomal recessive Charcot-Marie-Tooth neuropathy through aberrant splicing

BACKGROUND

Loss-of-function variants in MME (membrane metalloendopeptidase) are a known cause of recessive Charcot-Marie-Tooth Neuropathy (CMT). A deep intronic variant, MME c.1188+428A>G (NM_000902.5), was identified through whole genome sequencing (WGS) of two Australian families with recessive inheritance of axonal CMT using the seqr platform. MME c.1188+428A>G was detected in a homozygous state in Family 1, and in a compound heterozygous state with a known pathogenic MME variant (c.467del; p.Pro156Leufs*14) in Family 2.

AIMS

We aimed to determine the pathogenicity of the MME c.1188+428A>G variant through segregation and splicing analysis.

METHODS

The splicing impact of the deep intronic MME variant c.1188+428A>G was assessed using an in vitro exon-trapping assay.

RESULTS

The exon-trapping assay demonstrated that the MME c.1188+428A>G variant created a novel splice donor site resulting in the inclusion of an 83 bp pseudoexon between MME exons 12 and 13. The incorporation of the pseudoexon into MME transcript is predicted to lead to a coding frameshift and premature termination codon (PTC) in MME exon 14 (p.Ala397ProfsTer47). This PTC is likely to result in nonsense mediated decay (NMD) of MME transcript leading to a pathogenic loss-of-function.

INTERPRETATION

To our knowledge, this is the first report of a pathogenic deep intronic MME variant causing CMT. This is of significance as deep intronic variants are missed using whole exome sequencing screening methods. Individuals with CMT should be reassessed for deep intronic variants, with splicing impacts being considered in relation to the potential pathogenicity of variants.

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.

Reduced cone photoreceptor function and subtle systemic manifestations in two siblings with loss of SCLT1

BACKGROUND

The sodium channel and clathrin linker 1 gene (SCLT1) has been involved in the pathogenesis of various ciliopathy disorders such as Bardet-Biedl syndrome, orofaciodigital syndrome type IX, and Senior-Løken syndrome. Detailed exams are warranted to outline all clinical features. Here, we present a family with a milder phenotype of SCLT1-related disease.

MATERIAL AND METHODS

Comprehensive eye examination including fundus images, OCT, color vision, visual fields and electroretinography were performed. Affected individuals were assessed by a pediatrician and a medical geneticist for systemic features of ciliopathy. Investigations included echocardiography, abdominal ultrasonography, blood work-up for diabetes, liver and kidney function. Genetic testing included NGS retinal dystrophy panel, segregation analysis and transcriptome sequencing.

RESULTS

Two male children, age 10 and 8 years, were affected with attention deficit hyperactivity disorder (ADHD), obesity and mild photophobia. The ophthalmic exam revealed reduced best-corrected visual acuity (BCVA), strabismus, hyperopia, astigmatism and moderate red-green defects. Milder changes suggesting photoreceptors disease were found on retinal imaging. Electroretinogram confirmed cone photoreceptors dysfunction. Genetic testing revealed a homozygous likely pathogenic, splice-site variant in SCLT1 gene NM_144643.3: c.1439 + 1del in the proband and in the affected brother. The unaffected parents were heterozygous for the SCLT1 variant. Transcriptome sequencing showed retention of intron 16 in the proband.

CONCLUSIONS

In this report, we highlight the importance of further extensive diagnostics in patients with unexplained reduced vision, strabismus, refractive errors and ADHD spectrum disorders. SCLT1-related retinal degeneration is very rare and isolated reduced function of cone photoreceptors has not previously been observed.

Combining a prioritization strategy and functional studies nominates 5'UTR variants underlying inherited retinal disease

BACKGROUND

5' untranslated regions (5'UTRs) are essential modulators of protein translation. Predicting the impact of 5'UTR variants is challenging and rarely performed in routine diagnostics. Here, we present a combined approach of a comprehensive prioritization strategy and functional assays to evaluate 5'UTR variation in two large cohorts of patients with inherited retinal diseases (IRDs).

METHODS

We performed an isoform-level re-analysis of retinal RNA-seq data to identify the protein-coding transcripts of 378 IRD genes with highest expression in retina. We evaluated the coverage of their 5'UTRs by different whole exome sequencing (WES) kits. The selected 5'UTRs were analyzed in whole genome sequencing (WGS) and WES data from IRD sub-cohorts from the 100,000 Genomes Project (n = 2397 WGS) and an in-house database (n = 1682 WES), respectively. Identified variants were annotated for 5'UTR-relevant features and classified into seven categories based on their predicted functional consequence. We developed a variant prioritization strategy by integrating population frequency, specific criteria for each category, and family and phenotypic data. A selection of candidate variants underwent functional validation using diverse approaches.

RESULTS

Isoform-level re-quantification of retinal gene expression revealed 76 IRD genes with a non-canonical retina-enriched isoform, of which 20 display a fully distinct 5'UTR compared to that of their canonical isoform. Depending on the probe design, 3-20% of IRD genes have 5'UTRs fully captured by WES. After analyzing these regions in both cohorts, we prioritized 11 (likely) pathogenic variants in 10 genes (ARL3, MERTK, NDP, NMNAT1, NPHP4, PAX6, PRPF31, PRPF4, RDH12, RD3), of which 7 were novel. Functional analyses further supported the pathogenicity of three variants. Mis-splicing was demonstrated for the PRPF31:c.-9+1G>T variant. The MERTK:c.-125G>A variant, overlapping a transcriptional start site, was shown to significantly reduce both luciferase mRNA levels and activity. The RDH12:c.-123C>T variant was found in cis with the hypomorphic RDH12:c.701G>A (p.Arg234His) variant in 11 patients. This 5'UTR variant, predicted to introduce an upstream open reading frame, was shown to result in reduced RDH12 protein but unaltered mRNA levels.

CONCLUSIONS

This study demonstrates the importance of 5'UTR variants implicated in IRDs and provides a systematic approach for 5'UTR annotation and validation that is applicable to other inherited diseases.

The Value of Electroretinography in Identifying Candidate Genes for Inherited Retinal Dystrophies: A Diagnostic Guide

Inherited retinal dystrophies (IRDs) are a group of heterogeneous diseases caused by genetic mutations that specifically affect the function of the rod, cone, or bipolar cells in the retina. Electroretinography (ERG) is a diagnostic tool that measures the electrical activity of the retina in response to light stimuli, and it can help to determine the function of these cells. A normal ERG response consists of two waves, the a-wave and the b-wave, which reflect the activity of the photoreceptor cells and the bipolar and Muller cells, respectively. Despite the growing availability of next-generation sequencing (NGS) technology, identifying the precise genetic mutation causing an IRD can be challenging and costly. However, certain types of IRDs present with unique ERG features that can help guide genetic testing. By combining these ERG findings with other clinical information, such as on family history and retinal imaging, physicians can effectively narrow down the list of candidate genes to be sequenced, thereby reducing the cost of genetic testing. This review article focuses on certain types of IRDs with unique ERG features. We will discuss the pathophysiology and clinical presentation of, and ERG findings on, these disorders, emphasizing the unique role ERG plays in their diagnosis and genetic testing.

EyeG2P: an automated variant filtering approach improves efficiency of diagnostic genomic testing for inherited ophthalmic disorders

BACKGROUND

Genomic variant prioritisation is one of the most significant bottlenecks to mainstream genomic testing in healthcare. Tools to improve precision while ensuring high recall are critical to successful mainstream clinical genomic testing, in particular for whole genome sequencing where millions of variants must be considered for each patient.

METHODS

We developed EyeG2P, a publicly available database and web application using the Ensembl Variant Effect Predictor. EyeG2P is tailored for efficient variant prioritisation for individuals with inherited ophthalmic conditions. We assessed the sensitivity of EyeG2P in 1234 individuals with a broad range of eye conditions who had previously received a confirmed molecular diagnosis through routine genomic diagnostic approaches. For a prospective cohort of 83 individuals, we assessed the precision of EyeG2P in comparison with routine diagnostic approaches. For 10 additional individuals, we assessed the utility of EyeG2P for whole genome analysis.

RESULTS

EyeG2P had 99.5% sensitivity for genomic variants previously identified as clinically relevant through routine diagnostic analysis (n=1234 individuals). Prospectively, EyeG2P enabled a significant increase in precision (35% on average) in comparison with routine testing strategies (p<0.001). We demonstrate that incorporation of EyeG2P into whole genome sequencing analysis strategies can reduce the number of variants for analysis to six variants, on average, while maintaining high diagnostic yield.

CONCLUSION

Automated filtering of genomic variants through EyeG2P can increase the efficiency of diagnostic testing for individuals with a broad range of inherited ophthalmic disorders.

Intronic position +9 and -9 are potentially splicing sites boundary from intronic variants analysis of whole exome sequencing data

Whole exome sequencing (WES) can also detect some intronic variants, which may affect splicing and gene expression, but how to use these intronic variants, and the characteristics about them has not been reported. This study aims to reveal the characteristics of intronic variant in WES data, to further improve the clinical diagnostic value of WES. A total of 269 WES data was analyzed, 688,778 raw variants were called, among these 367,469 intronic variants were in intronic regions flanking exons which was upstream/downstream region of the exon (default is 200 bps). Contrary to expectation, the number of intronic variants with quality control (QC) passed was the lowest at the +2 and -2 positions but not at the +1 and -1 positions. The plausible explanation was that the former had the worst effect on trans-splicing, whereas the latter did not completely abolish splicing. And surprisingly, the number of intronic variants that passed QC was the highest at the +9 and -9 positions, indicating a potential splicing site boundary. The proportion of variants which could not pass QC filtering (false variants) in the intronic regions flanking exons generally accord with "S"-shaped curve. At +5 and -5 positions, the number of variants predicted damaging by software was most. This was also the position at which many pathogenic variants had been reported in recent years. Our study revealed the characteristics of intronic variant in WES data for the first time, we found the +9 and -9 positions might be a potentially splicing sites boundary and +5 and -5 positions were potentially important sites affecting splicing or gene expression, the +2 and -2 positions seem more important splicing site than +1 and -1 positions, and we found variants in intronic regions flanking exons over ± 50 bps may be unreliable. This result can help researchers find more useful variants and demonstrate that WES data is valuable for intronic variants analysis.

Systematic assessment of the contribution of structural variants to inherited retinal diseases

Despite increasing success in determining genetic diagnosis for patients with inherited retinal diseases (IRDs), mutations in about 30% of the IRD cases remain unclear or unsettled after targeted gene panel or whole exome sequencing. In this study, we aimed to investigate the contributions of structural variants (SVs) to settling the molecular diagnosis of IRD with whole-genome sequencing (WGS). A cohort of 755 IRD patients whose pathogenic mutations remain undefined were subjected to WGS. Four SV calling algorithms including include MANTA, DELLY, LUMPY and CNVnator were used to detect SVs throughout the genome. All SVs identified by any one of these four algorithms were included for further analysis. AnnotSV was used to annotate these SVs. SVs that overlap with known IRD-associated genes were examined with sequencing coverage, junction reads and discordant read pairs. Polymerase Chain Reaction (PCR) followed by Sanger sequencing was used to further confirm the SVs and identify the breakpoints. Segregation of the candidate pathogenic alleles with the disease was performed when possible. A total of 16 candidate pathogenic SVs were identified in 16 families, including deletions and inversions, representing 2.1% of patients with previously unsolved IRDs. Autosomal dominant, autosomal recessive and X-linked inheritance of disease-causing SVs were observed in 12 different genes. Among these, SVs in CLN3, EYS and PRPF31 were found in multiple families. Our study suggests that the contribution of SVs detected by short-read WGS is about 0.25% of our IRD patient cohort and is significantly lower than that of single nucleotide changes and small insertions and deletions.

Introns: the "dark matter" of the eukaryotic genome

The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the "dark matter" of the eukaryotic genome and discuss some of the outstanding questions in the field.

Novel Pathogenic Mutations Identified from Whole-Genome Sequencing in Unsolved Cases of Patients Affected with Inherited Retinal Diseases

Inherited retinal diseases (IRDs) are a diverse set of visual disorders that collectively represent a major cause of early-onset blindness. With the reduction in sequencing costs in recent years, whole-genome sequencing (WGS) is being used more frequently, particularly when targeted gene panels and whole-exome sequencing (WES) fail to detect pathogenic mutations in patients. In this study, we performed mutation screens using WGS for a cohort of 311 IRD patients whose mutations were undetermined. A total of nine putative pathogenic mutations in six IRD patients were identified, including six novel mutations. Among them, four were deep intronic mutations that affected mRNA splicing, while the other five affected protein-coding sequences. Our results suggested that the rate of resolution of unsolved cases via targeted gene panels and WES can be further enhanced with WGS; however, the overall improvement may be limited.

Long-Read Sequencing Identifies Novel Pathogenic Intronic Variants in Gitelman Syndrome

BACKGROUND

Gitelman syndrome is a salt-losing tubulopathy characterized by hypokalemic alkalosis and hypomagnesemia. It is caused by homozygous recessive or compound heterozygous pathogenic variants in SLC12A3 , which encodes the Na + -Cl - cotransporter (NCC). In up to 10% of patients with Gitelman syndrome, current genetic techniques detect only one specific pathogenic variant. This study aimed to identify a second pathogenic variant in introns, splice sites, or promoters to increase the diagnostic yield.

METHODS

Long-read sequencing of SLC12A3 was performed in 67 DNA samples from individuals with suspected Gitelman syndrome in whom a single likely pathogenic or pathogenic variant was previously detected. In addition, we sequenced DNA samples from 28 individuals with one variant of uncertain significance or no candidate variant. Midigene splice assays assessed the pathogenicity of novel intronic variants.

RESULTS

A second likely pathogenic/pathogenic variant was identified in 45 (67%) patients. Those with two likely pathogenic/pathogenic variants had a more severe electrolyte phenotype than other patients. Of the 45 patients, 16 had intronic variants outside of canonic splice sites (nine variants, mostly deep intronic, six novel), whereas 29 patients had an exonic variant or canonic splice site variant. Midigene splice assays of the previously known c.1670-191C>T variant and intronic candidate variants demonstrated aberrant splicing patterns.

CONCLUSION

Intronic pathogenic variants explain an important part of the missing heritability in Gitelman syndrome. Long-read sequencing should be considered in diagnostic workflows for Gitelman syndrome.

Systematic assessment of the contribution of structural variants to inherited retinal diseases

Despite increasing success in determining genetic diagnosis for patients with inherited retinal diseases (IRDs), mutations in about 30% of the IRD cases remain unclear or unsettled after targeted gene panel or whole exome sequencing. In this study, we aimed to investigate the contributions of structural variants (SVs) to settling the molecular diagnosis of IRD with whole-genome sequencing (WGS). A cohort of 755 IRD patients whose pathogenic mutations remain undefined was subjected to WGS. Four SV calling algorithms including include MANTA, DELLY, LUMPY, and CNVnator were used to detect SVs throughout the genome. All SVs identified by any one of these four algorithms were included for further analysis. AnnotSV was used to annotate these SVs. SVs that overlap with known IRD-associated genes were examined with sequencing coverage, junction reads, and discordant read pairs. PCR followed by Sanger sequencing was used to further confirm the SVs and identify the breakpoints. Segregation of the candidate pathogenic alleles with the disease was performed when possible. In total, sixteen candidate pathogenic SVs were identified in sixteen families, including deletions and inversions, representing 2.1% of patients with previously unsolved IRDs. Autosomal dominant, autosomal recessive, and X-linked inheritance of disease-causing SVs were observed in 12 different genes. Among these, SVs in CLN3, EYS, PRPF31 were found in multiple families. Our study suggests that the contribution of SVs detected by short-read WGS is about 0.25% of our IRD patient cohort and is significantly lower than that of single nucleotide changes and small insertions and deletions.

Multi-omics approach dissects cis-regulatory mechanisms underlying North Carolina macular dystrophy, a retinal enhanceropathy

North Carolina macular dystrophy (NCMD) is a rare autosomal-dominant disease affecting macular development. The disease is caused by non-coding single-nucleotide variants (SNVs) in two hotspot regions near PRDM13 and by duplications in two distinct chromosomal loci, overlapping DNase I hypersensitive sites near either PRDM13 or IRX1. To unravel the mechanisms by which these variants cause disease, we first established a genome-wide multi-omics retinal database, RegRet. Integration of UMI-4C profiles we generated on adult human retina then allowed fine-mapping of the interactions of the PRDM13 and IRX1 promoters and the identification of eighteen candidate cis-regulatory elements (cCREs), the activity of which was investigated by luciferase and Xenopus enhancer assays. Next, luciferase assays showed that the non-coding SNVs located in the two hotspot regions of PRDM13 affect cCRE activity, including two NCMD-associated non-coding SNVs that we identified herein. Interestingly, the cCRE containing one of these SNVs was shown to interact with the PRDM13 promoter, demonstrated in vivo activity in Xenopus, and is active at the developmental stage when progenitor cells of the central retina exit mitosis, suggesting that this region is a PRDM13 enhancer. Finally, mining of single-cell transcriptional data of embryonic and adult retina revealed the highest expression of PRDM13 and IRX1 when amacrine cells start to synapse with retinal ganglion cells, supporting the hypothesis that altered PRDM13 or IRX1 expression impairs interactions between these cells during retinogenesis. Overall, this study provides insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy.

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.

CRISPR-Cas9-mediated functional dissection of the foxc1 genomic region in zebrafish identifies critical conserved cis-regulatory elements

FOXC1 encodes a forkhead-domain transcription factor associated with several ocular disorders. Correct FOXC1 dosage is critical to normal development, yet the mechanisms controlling its expression remain unknown. Together with FOXQ1 and FOXF2, FOXC1 is part of a cluster of FOX genes conserved in vertebrates. CRISPR-Cas9-mediated dissection of genomic sequences surrounding two zebrafish orthologs of FOXC1 was performed. This included five zebrafish-human conserved regions, three downstream of foxc1a and two remotely upstream of foxf2a/foxc1a or foxf2b/foxc1b clusters, as well as two intergenic regions between foxc1a/b and foxf2a/b lacking sequence conservation but positionally corresponding to the area encompassing a previously reported glaucoma-associated SNP in humans. Removal of downstream sequences altered foxc1a expression; moreover, zebrafish carrying deletions of two or three downstream elements demonstrated abnormal phenotypes including enlargement of the anterior chamber of the eye reminiscent of human congenital glaucoma. Deletions of distant upstream conserved elements influenced the expression of foxf2a/b or foxq1a/b but not foxc1a/b within each cluster. Removal of either intergenic sequence reduced foxc1a or foxc1b expression during late development, suggesting a role in transcriptional regulation despite the lack of conservation at the nucleotide level. Further studies of the identified regions in human patients may explain additional individuals with developmental ocular disorders.

Noncanonical Splice Site and Deep Intronic FRMD7 Variants Activate Cryptic Exons in X-linked Infantile Nystagmus

Purpose

We aim to report noncoding pathogenic variants in patients with FRMD7-related infantile nystagmus (FIN).

Methods

Genome sequencing (n = 2 families) and reanalysis of targeted panel next generation sequencing (n = 2 families) was performed in genetically unsolved cases of suspected FIN. Previous sequence analysis showed no pathogenic coding variants in genes associated with infantile nystagmus. SpliceAI, SpliceRover, and Alamut consensus programs were used to annotate noncoding variants. Minigene splicing assay was performed to confirm aberrant splicing. In silico analysis of exonic splicing enhancer and silencer was also performed.

Results

FRMD7 intronic variants were identified based on genome sequencing and targeted next-generation sequencing analysis. These included c.285-12A>G (pedigree 1), c.284+63T>A (pedigrees 2 and 3), and c. 383-1368A>G (pedigree 4). All variants were absent in gnomAD, and the both c.285-12A>G and c.284+63T>A variants were predicted to enhance new splicing acceptor gains with SpliceAI, SpliceRover, and Alamut consensus approaches. However, the c.383-1368 A>G variant only had a significant impact score on the SpliceRover program. The c.383-1368A>G variant was predicted to promote pseudoexon inclusion by binding of exonic splicing enhancer. Aberrant exonizations were validated through minigene constructs, and all variants were segregated in the families.

Conclusions

Deep learning–based annotation of noncoding variants facilitates the discovery of hidden genetic variations in patients with FIN. This study provides evidence of effectiveness of combined deep learning–based splicing tools to identify hidden pathogenic variants in previously unsolved patients with infantile nystagmus.

Translational Relevance

These results demonstrate robust analysis using two deep learning splicing predictions and in vitro functional study can lead to finding hidden genetic variations in unsolved patients.

Case Review: Whole-Exome Sequencing Analyses Identify Carriers of a Known Likely Pathogenic Intronic BRCA1 Variant in Ovarian Cancer Cases Clinically Negative for Pathogenic BRCA1 and BRCA2 Variants

Background: Detecting pathogenic intronic variants resulting in aberrant splicing remains a challenge in routine genetic testing. We describe germline whole-exome sequencing (WES) analyses and apply in silico predictive tools of familial ovarian cancer (OC) cases reported clinically negative for pathogenic BRCA1 and BRCA2 variants. Methods: WES data from 27 familial OC cases reported clinically negative for pathogenic BRCA1 and BRCA2 variants and 53 sporadic early-onset OC cases were analyzed for pathogenic variants in BRCA1 or BRCA2. WES data from carriers of pathogenic BRCA1 or BRCA2 variants were analyzed for pathogenic variants in 10 other OC predisposing genes. Loss of heterozygosity analysis was performed on tumor DNA from variant carriers. Results: BRCA1 c.5407-25T>A intronic variant, identified in two affected sisters and one sporadic OC case, is predicted to create a new splice effecting transcription of BRCA1. WES data from BRCA1 c.5407-25T>A carriers showed no evidence of pathogenic variants in other OC predisposing genes. Sequencing the tumor DNA from the variant carrier showed complete loss of the wild-type allele. Conclusions: The findings support BRCA1 c.5407-25T>A as a likely pathogenic variant and highlight the importance of investigating intronic sequences as causal variants in OC families where the involvement of BRCA1 is highly suggestive.

A Homozygous Deep Intronic Variant Causes Von Willebrand Factor Deficiency and Lack of Endothelial-Specific Secretory Organelles, Weibel-Palade Bodies

A type 3 von Willebrand disease (VWD) index patient (IP) remains mutation-negative after completion of the conventional diagnostic analysis, including multiplex ligation-dependent probe amplification and sequencing of the promoter, exons, and flanking intronic regions of the VWF gene (VWF). In this study, we intended to elucidate causative mutation through next-generation sequencing (NGS) of the whole VWF (including complete intronic region), mRNA analysis, and study of the patient-derived endothelial colony-forming cells (ECFCs). The NGS revealed a variant in the intronic region of VWF (997 + 118 T > G in intron 8), for the first time. The bioinformatics assessments (e.g., SpliceAl) predicted this variant creates a new donor splice site (ss), which could outcompete the consensus 5′ donor ss at exon/intron 8. This would lead to an aberrant mRNA that contains a premature stop codon, targeting it to nonsense-mediated mRNA decay. The subsequent quantitative real-time PCR confirmed the virtual absence of VWF mRNA in IP ECFCs. Additionally, the IP ECFCs demonstrated a considerable reduction in VWF secretion (~6% of healthy donors), and they were devoid of endothelial-specific secretory organelles, Weibel−Palade bodies. Our findings underline the potential of NGS in conjunction with RNA analysis and patient-derived cell studies for genetic diagnosis of mutation-negative type 3 VWD patients.

Gene-Based Therapeutics for Inherited Retinal Diseases

Inherited retinal diseases (IRDs) are a heterogenous group of orphan eye diseases that typically result from monogenic mutations and are considered attractive targets for gene-based therapeutics. Following the approval of an IRD gene replacement therapy for Leber's congenital amaurosis due to RPE65 mutations, there has been an intensive international research effort to identify the optimal gene therapy approaches for a range of IRDs and many are now undergoing clinical trials. In this review we explore therapeutic challenges posed by IRDs and review current and future approaches that may be applicable to different subsets of IRD mutations. Emphasis is placed on five distinct approaches to gene-based therapy that have potential to treat the full spectrum of IRDs: 1) gene replacement using adeno-associated virus (AAV) and nonviral delivery vectors, 2) genome editing via the CRISPR/Cas9 system, 3) RNA editing by endogenous and exogenous ADAR, 4) mRNA targeting with antisense oligonucleotides for gene knockdown and splicing modification, and 5) optogenetic approaches that aim to replace the function of native retinal photoreceptors by engineering other retinal cell types to become capable of phototransduction.

Clinical and Genetic Re-Evaluation of Inherited Retinal Degeneration Pedigrees following Initial Negative Findings on Panel-Based Next Generation Sequencing

Although rare, inherited retinal degenerations (IRDs) are the most common reason for blind registration in the working age population. They are highly genetically heterogeneous (>300 known genetic loci), and confirmation of a molecular diagnosis is a prerequisite for many therapeutic clinical trials and approved treatments. First-tier genetic testing of IRDs with panel-based next-generation sequencing (pNGS) has a diagnostic yield of ≈70-80%, leaving the remaining more challenging cases to be resolved by second-tier testing methods. This study describes the phenotypic reassessment of patients with a negative result from first-tier pNGS and the rationale, outcomes, and cost of second-tier genetic testing approaches. Removing non-IRD cases from consideration and utilizing case-appropriate second-tier genetic testing techniques, we genetically resolved 56% of previously unresolved pedigrees, bringing the overall resolve rate to 92% (388/423). At present, pNGS remains the most cost-effective first-tier approach for the molecular assessment of diverse IRD populations Second-tier genetic testing should be guided by clinical (i.e., reassessment, multimodal imaging, electrophysiology), and genetic (i.e., single alleles in autosomal recessive disease) indications to achieve a genetic diagnosis in the most cost-effective manner.

RNA Analysis and Clinical Characterization of a Novel Splice Variant in the NSD1 Gene Causing Familial Sotos Syndrome

BACKGROUND

Sotos syndrome is an autosomal dominant disorder characterized by overgrowth, macrocephaly, distinctive facial features and learning disabilities. Haploinsufficiency of the nuclear receptor SET domain-containing protein 1 (NSD1) gene located on chromosome 5q35 is the major cause of the syndrome. This syndrome shares characteristics with other overgrowth syndromes, which can complicate the differential diagnosis.

METHODS

Genomic DNA was extracted from peripheral blood samples of members of the same family and targeted exome analysis was performed. In silico study of the variant found by next-generation sequencing was used to predict disruption/creation of splice sites and the identification of potential cryptic splice sites. RNA was extracted from peripheral blood samples of patients and functional analyses were performed to confirm the pathogenicity.

RESULTS

We found a novel c.6463 + 5G>A heterozygous NSD1 gene pathogenic variant in a son and his father. Molecular analyses revealed that part of the intron 22 of NSD1 is retained due to the destruction of the splicing donor site, causing the appearance of a premature stop codon in the NSD1 protein.

CONCLUSIONS

Our findings underline the importance of performing RNA functional assays in order to determine the clinical significance of intronic variants, and contribute to the genetic counseling and clinical management of patients and their relatives. Our work also highlights the relevance of using in silico prediction tools to detect a potential alteration in the splicing process.

Inherited retinal diseases: Linking genes, disease-causing variants, and relevant therapeutic modalities

Inherited retinal diseases (IRDs) are a clinically complex and heterogenous group of visual impairment phenotypes caused by pathogenic variants in at least 277 nuclear and mitochondrial genes, affecting different retinal regions, and depleting the vision of affected individuals. Genes that cause IRDs when mutated are unique by possessing differing genotype-phenotype correlations, varying inheritance patterns, hypomorphic alleles, and modifier genes thus complicating genetic interpretation. Next-generation sequencing has greatly advanced the identification of novel IRD-related genes and pathogenic variants in the last decade. For this review, we performed an in-depth literature search which allowed for compilation of the Global Retinal Inherited Disease (GRID) dataset containing 4,798 discrete variants and 17,299 alleles published in 31 papers, showing a wide range of frequencies and complexities among the 194 genes reported in GRID, with 65% of pathogenic variants being unique to a single individual. A better understanding of IRD-related gene distribution, gene complexity, and variant types allow for improved genetic testing and therapies. Current genetic therapeutic methods are also quite diverse and rely on variant identification, and range from whole gene replacement to single nucleotide editing at the DNA or RNA levels. IRDs and their suitable therapies thus require a range of effective disease modelling in human cells, granting insight into disease mechanisms and testing of possible treatments. This review summarizes genetic and therapeutic modalities of IRDs, provides new analyses of IRD-related genes (GRID and complexity scores), and provides information to match genetic-based therapies such as gene-specific and variant-specific therapies to the appropriate individuals.

Performance Evaluation of SpliceAI for the Prediction of Splicing of NF1 Variants

Neurofibromatosis type 1, characterized by neurofibromas and café-au-lait macules, is one of the most common genetic disorders caused by pathogenic NF1 variants. Because of the high proportion of splicing mutations in NF1, identifying variants that alter splicing may be an essential issue for laboratories. Here, we investigated the sensitivity and specificity of SpliceAI, a recently introduced in silico splicing prediction algorithm in conjunction with other in silico tools. We evaluated 285 NF1 variants identified from 653 patients. The effect on variants on splicing alteration was confirmed by complementary DNA sequencing followed by genomic DNA sequencing. For in silico prediction of splicing effects, we used SpliceAI, MaxEntScan (MES), and Splice Site Finder-like (SSF). The sensitivity and specificity of SpliceAI were 94.5% and 94.3%, respectively, with a cut-off value of Δ Score > 0.22. The area under the curve of SpliceAI was 0.975 (p < 0.0001). Combined analysis of MES/SSF showed a sensitivity of 83.6% and specificity of 82.5%. The concordance rate between SpliceAI and MES/SSF was 84.2%. SpliceAI showed better performance for the prediction of splicing alteration for NF1 variants compared with MES/SSF. As a convenient web-based tool, SpliceAI may be helpful in clinical laboratories conducting DNA-based NF1 sequencing.

Splicing Outcomes of 5' Splice Site GT>GC Variants That Generate Wild-Type Transcripts Differ Significantly Between Full-Length and Minigene Splicing Assays

Combining data derived from a meta-analysis of human disease-associated 5' splice site GT>GC (i.e., +2T>C) variants and a cell culture-based full-length gene splicing assay (FLGSA) of forward engineered +2T>C substitutions, we recently estimated that ∼15-18% of +2T>C variants can generate up to 84% wild-type transcripts relative to their wild-type counterparts. Herein, we analyzed the splicing outcomes of 20 +2T>C variants that generate some wild-type transcripts in two minigene assays. We found a high discordance rate in terms of the generation of wild-type transcripts, not only between FLGSA and the minigene assays but also between the different minigene assays. In the pET01 context, all 20 wild-type minigene constructs generated the expected wild-type transcripts; of the 20 corresponding variant minigene constructs, 14 (70%) generated wild-type transcripts. In the pSPL3 context, only 18 of the 20 wild-type minigene constructs generated the expected wild-type transcripts whereas 8 of the 18 (44%) corresponding variant minigene constructs generated wild-type transcripts. Thus, in the context of a particular type of variant, we raise awareness of the limitations of minigene splicing assays and emphasize the importance of sequence context in regulating splicing. Whether or not our findings apply to other types of splice-altering variant remains to be investigated.

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.