Functional Characterization of the Human BRCA1 ∆11 Splicing Isoforms in Yeast.

BRCA1, a crucial tumor suppressor gene, has several splicing isoforms, including Δ9-11, Δ11, and Δ11q, which lack exon 11, coding for significant portions of the protein. These isoforms are naturally present in both normal and cancerous cells, exhibiting altered activity compared to the full-length BRCA1. Despite this, the impact on cancer risk of the germline intronic variants promoting the exclusive expression of these Δ11 isoforms remains uncertain. Consequently, they are classified as variants of uncertain significance (VUS), posing challenges for traditional genetic classification methods due to their rarity and complexity. Our research utilizes a yeast-based functional assay, previously validated for assessing missense BRCA1 variants, to compare the activity of the Δ11 splicing isoforms with known pathogenic missense variants. This approach allows us to elucidate the functional implications of these isoforms and determine whether their exclusive expression could contribute to increased cancer risk. By doing so, we aim to provide insights into the pathogenic potential of intronic VUS-generating BRCA1 splicing isoforms and improve the classification of BRCA1 variants.

Limited Added Diagnostic Value of Whole Genome Sequencing in Genetic Testing of Inherited Retinal Diseases in a Swiss Patient Cohort.

The purpose of this study was to assess the added diagnostic value of whole genome sequencing (WGS) for patients with inherited retinal diseases (IRDs) who remained undiagnosed after whole exome sequencing (WES). WGS was performed for index patients in 66 families. The datasets were analyzed according to GATK's guidelines. Additionally, DeepVariant was complemented by GATK's workflow, and a novel structural variant pipeline was developed. Overall, a molecular diagnosis was established in 19/66 (28.8%) index patients. Pathogenic deletions and one deep-intronic variant contributed to the diagnostic yield in 4/19 and 1/19 index patients, respectively. The remaining diagnoses (14/19) were attributed to exonic variants that were missed during WES analysis due to bioinformatic limitations, newly described loci, or unclear pathogenicity. The added diagnostic value of WGS equals 5/66 (9.6%) for our cohort, which is comparable to previous studies. This figure would decrease further to 1/66 (1.5%) with a standardized and reliable copy number variant workflow during WES analysis. Given the higher costs and limited added value, the implementation of WGS as a first-tier assay for inherited eye disorders in a diagnostic laboratory remains untimely. Instead, progress in bioinformatic tools and communication between diagnostic and clinical teams have the potential to ameliorate diagnostic yields.

Molecular Challenges in the Diagnosis of X-Linked Chronic Granulomatous Disease: CNVs, Intronic Variants, Skewed X-Chromosome Inactivation, and Gonosomal Mosaicism.

Chronic granulomatous disease (CGD) is a prototypical inborn error of immunity affecting phagocytes, in which these cells are unable to produce reactive oxygen species. CGD is caused by defects in genes encoding subunits of the NADPH oxidase enzyme complex (CYBA, CYBB, CYBC1, NCF1, NCF2, NCF4); inflammatory responses are dysregulated, and patients are highly susceptible to recurrent severe bacterial and fungal infections. X-linked CGD (XL-CGD), caused by mutations in the CYBB gene, is the most common and severe form of CGD. In this study, we describe the analytical processes undertaken in 3 families affected with XL-CGD to illustrate several molecular challenges in the genetic diagnosis of this condition: in family 1, a girl with a heterozygous deletion of CYBB exon 13 and skewed X-chromosome inactivation (XCI); in family 2, a boy with a hemizygous deletion of CYBB exon 7, defining its consequences at the mRNA level; and in family 3, 2 boys with the same novel intronic variant in CYBB (c.1151 + 6 T > A). The variant affected the splicing process, although a small fraction of wild-type mRNA was produced. Their mother was a heterozygous carrier, while their maternal grandmother was a carrier in form of gonosomal mosaicism. In summary, using a variety of techniques, including an NGS-based targeted gene panel and deep amplicon sequencing, copy number variation calling strategies, microarray-based comparative genomic hybridization, and cDNA analysis to define splicing defects and skewed XCI, we show how to face and solve some uncommon genetic mechanisms in the diagnosis of XL-CGD.

Overcoming the challenges associated with identification of deep intronic variants by whole genome sequencing.

Whole-genome sequencing (WGS) now allows identification of multiple variants in non-coding regions. The large number of variants identified by WGS however complicates their interpretation. Through identification of the first deep intronic variant in NPHS2, which encodes podocin, a protein implicated in autosomal recessive steroid resistant nephrotic syndrome (SRNS), we compare herein three different tools including a newly developed targeted NGS-based RNA-sequencing to explore the splicing effect of intronic variations. WGS identified two different variants in NPHS2 eventually involved in the disease. Through RT-PCR, exon-trapping Minigene assay and targeted RNA sequencing, we were able to identify the splicing defect in NPHS2 mRNA from patient kidney tissue. Only targeted RNA-seq simultaneously analyzed the effect of multiple variants and offered the opportunity to quantify consequences on splicing. Identifying deep intronic variants and their role in disease is of utmost importance. Alternative splicing can be predicted by in silico tools but always requires confirmation through functional testing with RNA analysis from the implicated tissue remaining the gold standard. When several variants with potential effects on splicing are identified by WGS, a targeted RNA sequencing panel could be of great value.

Minigene Splicing Assays and Long-Read Sequencing to Unravel Pathogenic Deep-Intronic Variants in PAX6 in Congenital Aniridia.

PAX6 haploinsufficiency causes aniridia, a congenital eye disorder that involves the iris, and foveal hypoplasia. Comprehensive screening of the PAX6 locus, including the non-coding regions, by next-generation sequencing revealed four deep-intronic variants with potential effects on pre-RNA splicing. Nevertheless, without a functional analysis, their pathogenicity could not be established. We aimed to decipher their impact on the canonical PAX6 splicing using in vitro minigene splicing assays and nanopore-based long-read sequencing. Two multi-exonic PAX6 constructs were generated, and minigene assays were carried out. An aberrant splicing pattern was observed for two variants in intron 6, c.357+136G>A and c.357+334G>A. In both cases, several exonization events, such as pseudoexon inclusions and partial intronic retention, were observed due to the creation or activation of new/cryptic non-canonical splicing sites, including a shared intronic donor site. In contrast, two variants identified in intron 11, c.1032+170A>T and c.1033-275A>C, seemed not to affect splicing processes. We confirmed the high complexity of alternative splicing of PAX6 exon 6, which also involves unreported cryptic intronic sites. Our study highlights the importance of integrating functional studies into diagnostic algorithms to decipher the potential implication of non-coding variants, usually classified as variants of unknown significance, thus allowing variant reclassification to achieve a conclusive genetic diagnosis.

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.

Unpredicted Aberrant Splicing Products Identified in Postmortem Sudden Cardiac Death Samples.

Molecular screening for pathogenic mutations in sudden cardiac death (SCD)-related genes is common practice for SCD cases. However, test results may lead to uncertainty because of the identification of variants of unknown significance (VUS) occurring in up to 70% of total identified variants due to a lack of experimental studies. Genetic variants affecting potential splice site variants are among the most difficult to interpret. The aim of this study was to examine rare intronic variants identified in the exonic flanking sequence to meet two main objectives: first, to validate that canonical intronic variants produce aberrant splicing; second, to determine whether rare intronic variants predicted as VUS may affect the splicing product. To achieve these objectives, 28 heart samples of cases of SCD carrying rare intronic variants were studied. Samples were analyzed using 85 SCD genes in custom panel sequencing. Our results showed that rare intronic variants affecting the most canonical splice sites displayed in 100% of cases that they would affect the splicing product, possibly causing aberrant isoforms. However, 25% of these cases (1/4) showed normal splicing, contradicting the in silico results. On the contrary, in silico results predicted an effect in 0% of cases, and experimental results showed >20% (3/14) unpredicted aberrant splicing. Thus, deep intron variants are likely predicted to not have an effect, which, based on our results, might be an underestimation of their effect and, therefore, of their pathogenicity classification and family members' follow-up.

Identification of Spliceogenic Variants beyond Canonical GT-AG Splice Sites in Hereditary Cancer Genes.

Pathogenic/likely pathogenic variants in susceptibility genes that interrupt RNA splicing are a well-documented mechanism of hereditary cancer syndromes development. However, if RNA studies are not performed, most of the variants beyond the canonical GT-AG splice site are characterized as variants of uncertain significance (VUS). To decrease the VUS burden, we have bioinformatically evaluated all novel VUS detected in 732 consecutive patients tested in the routine genetic counseling process. Twelve VUS that were predicted to cause splicing defects were selected for mRNA analysis. Here, we report a functional characterization of 12 variants located beyond the first two intronic nucleotides using RNAseq in APC, ATM, FH, LZTR1, MSH6, PALB2, RAD51C, and TP53 genes. Based on the analysis of mRNA, we have successfully reclassified 50% of investigated variants. 25% of variants were downgraded to likely benign, whereas 25% were upgraded to likely pathogenic leading to improved clinical management of the patient and the family members.

Intronic Haplotypes in the GBA Gene Do Not Predict Age at Diagnosis of Parkinson's Disease.

BACKGROUND: GBA mutations are a common risk factor for Parkinson's disease (PD). A recent study has suggested that GBA haplotypes, identified by intronic variants, can affect age at diagnosis of PD. OBJECTIVES: In this study, we assess this hypothesis using long reads across a large cohort and the publicly available Accelerating Medicines Partnership-Parkinson's Disease (AMP-PD) cohort. METHODS: We recruited a PD cohort through the Remote Assessment of Parkinsonism Supporting Ongoing Development of Interventions in Gaucher Disease study (RAPSODI) and sequenced GBA using Oxford Nanopore technology. Genetic and clinical data on the full AMP-PD cohort were obtained from the online portal of the consortium. RESULTS: A total of 1417 participants were analyzed. There was no significant difference in age at PD diagnosis between the two main haplotypes of the GBA gene. CONCLUSIONS: GBA haplotypes do not affect age at diagnosis of PD in the two independent cohorts studied. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

The Alter Retina: Alternative Splicing of Retinal Genes in Health and Disease.

Alternative splicing of mRNA is an essential mechanism to regulate and increase the diversity of the transcriptome and proteome. Alternative splicing frequently occurs in a tissue- or time-specific manner, contributing to differential gene expression between cell types during development. Neural tissues present extremely complex splicing programs and display the highest number of alternative splicing events. As an extension of the central nervous system, the retina constitutes an excellent system to illustrate the high diversity of neural transcripts. The retina expresses retinal specific splicing factors and produces a large number of alternative transcripts, including exclusive tissue-specific exons, which require an exquisite regulation. In fact, a current challenge in the genetic diagnosis of inherited retinal diseases stems from the lack of information regarding alternative splicing of retinal genes, as a considerable percentage of mutations alter splicing or the relative production of alternative transcripts. Modulation of alternative splicing in the retina is also instrumental in the design of novel therapeutic approaches for retinal dystrophies, since it enables precision medicine for specific mutations.

RT-PCR analysis of mRNA revealed the splice-altering effect of rare intronic variants in monogenic disorders.

BACKGROUND: Variants perturbing the normal splicing of pre-mRNA can lead to human diseases. The splice-altering effect and eventual consequence on gene function was sometimes uncertain and hinders a definitive molecular diagnosis. METHODS: The impact of four rare intronic variants on splicing was analyzed through reverse transcription - polymerase chain reaction (RT-PCR) analysis of mRNA derived from the peripheral blood of patients. The results were compared with in-silico prediction. Potential implication on molecular diagnosis was discussed. RESULTS: Four rare intronic variants of SLC9A6, DLG3, GAA, and OCRL were identified in patients with suspected disorders, respectively. Although these four variants were all predicted to alter splicing by in-silico tools, RT-PCR analysis of mRNA derived from peripheral blood showed these variants affected splicing in different ways: c.899+3_899+6del of SLC9A6 resulted in one-exon skipping and an out-of-frame transcript; c.905-2A > G of DLG3 resulted in a mix of in-frame transcripts; c.1195-11T > A of GAA resulted in the in-frame insertion of nine nucleotides; c.723-2A > C of OCRL resulted in one-exon skipping and in-frame deletion of 102 nucleotides. The consequence revealed by mRNA analysis is essential for accurate interpretation of pathogenicity. CONCLUSION: Four intronic variants all caused aberrant mRNA splicing. For intronic variants with uncertain impact on splicing, mRNA analysis is helpful for ascertainment of alternative splicing and accurate interpretation of pathogenicity.

Resolving the dark matter of ABCA4 for 1054 Stargardt disease probands through integrated genomics and transcriptomics.

PURPOSE: Missing heritability in human diseases represents a major challenge, and this is particularly true for ABCA4-associated Stargardt disease (STGD1). We aimed to elucidate the genomic and transcriptomic variation in 1054 unsolved STGD and STGD-like probands. METHODS: Sequencing of the complete 128-kb ABCA4 gene was performed using single-molecule molecular inversion probes (smMIPs), based on a semiautomated and cost-effective method. Structural variants (SVs) were identified using relative read coverage analyses and putative splice defects were studied using in vitro assays. RESULTS: In 448 biallelic probands 14 known and 13 novel deep-intronic variants were found, resulting in pseudoexon (PE) insertions or exon elongations in 105 alleles. Intriguingly, intron 13 variants c.1938-621G>A and c.1938-514G>A resulted in dual PE insertions consisting of the same upstream, but different downstream PEs. The intron 44 variant c.6148-84A>T resulted in two PE insertions and flanking exon deletions. Eleven distinct large deletions were found, two of which contained small inverted segments. Uniparental isodisomy of chromosome 1 was identified in one proband. CONCLUSION: Deep sequencing of ABCA4 and midigene-based splice assays allowed the identification of SVs and causal deep-intronic variants in 25% of biallelic STGD1 cases, which represents a model study that can be applied to other inherited diseases.

Toward a clinical diagnostic pipeline for SPINK1 intronic variants.

BACKGROUND: The clinical significance of SPINK1 intronic variants in chronic pancreatitis has been previously assessed by various approaches including a cell culture-based full-length gene assay. A close correlation between the results of this assay and in silico splicing prediction was apparent. However, until now, a clinical diagnostic pipeline specifically designed to classify SPINK1 intronic variants accurately and efficiently has been lacking. Herein, we present just such a pipeline and explore its efficacy and potential utility in potentiating the classification of newly described SPINK1 intronic variants. RESULTS: We confirm a close correlation between in silico splicing prediction and results from the cell culture-based full-length gene assay in the context of three recently reported pathogenic SPINK1 intronic variants. We then integrated in silico splicing prediction and the full-length gene assay into a stepwise approach and tested its utility in the classification of two novel datasets of SPINK1 intronic variants. The first dataset comprised 16 deep intronic variants identified in 52 genetically unexplained Chinese chronic pancreatitis patients by sequencing the entire intronic sequence of the SPINK1 gene. The second dataset comprised five novel rare proximal intronic variants identified through the routine analysis of the SPINK1 gene in French pancreatitis patients. Employing a minor allele frequency of > 5% as a population frequency filter, 6 of the 16 deep intronic variants were immediately classified as benign. In silico prediction of the remaining ten deep intronic variants and the five rare proximal intronic variants with respect to their likely impact on splice site selection suggested that only one proximal intronic variant, c.194 + 5G > A, was likely to be of functional significance. Employing the cell culture-based full-length gene assay, we functionally analyzed c.194 + 5G > A, together with seven predicted non-functional variants, thereby validating their predicted effects on splicing in all cases. CONCLUSIONS: We demonstrated the accuracy and efficiency of in silico prediction in combination with the cell culture-based full-length gene assay for the classification of SPINK1 intronic variants. Based upon these findings, we propose an operational pipeline for classifying SPINK1 intronic variants in the clinical diagnostic setting.

Cost-effective molecular inversion probe-based ABCA4 sequencing reveals deep-intronic variants in Stargardt disease.

PURPOSE: Stargardt disease (STGD1) is caused by biallelic mutations in ABCA4, but many patients are genetically unsolved due to insensitive mutation-scanning methods. We aimed to develop a cost-effective sequencing method for ABCA4 exons and regions carrying known causal deep-intronic variants. METHODS: Fifty exons and 12 regions containing 14 deep-intronic variants of ABCA4 were sequenced using double-tiled single molecule Molecular Inversion Probe (smMIP)-based next-generation sequencing. DNAs of 16 STGD1 cases carrying 29 ABCA4 alleles and of four healthy persons were sequenced using 483 smMIPs. Thereafter, DNAs of 411 STGD1 cases with one or no ABCA4 variant were sequenced. The effect of novel noncoding variants on splicing was analyzed using in vitro splice assays. RESULTS: Thirty-four ABCA4 variants previously identified in 16 STGD1 cases were reliably identified. In 155/411 probands (38%), two causal variants were identified. We identified 11 deep-intronic variants present in 62 alleles. Two known and two new noncanonical splice site variants showed splice defects, and one novel deep-intronic variant (c.4539+2065C>G) resulted in a 170-nt mRNA pseudoexon insertion (p.[Arg1514Lysfs*35,=]). CONCLUSIONS: smMIPs-based sequence analysis of coding and selected noncoding regions of ABCA4 enabled cost-effective mutation detection in STGD1 cases in previously unsolved cases.

Identification of splice defects due to noncanonical splice site or deep-intronic variants in ABCA4.

Pathogenic variants in the ATP-binding cassette transporter A4 (ABCA4) gene cause a continuum of retinal disease phenotypes, including Stargardt disease. Noncanonical splice site (NCSS) and deep-intronic variants constitute a large fraction of disease-causing alleles, defining the functional consequences of which remains a challenge. We aimed to determine the effect on splicing of nine previously reported or unpublished NCSS variants, one near exon splice variant and nine deep-intronic variants in ABCA4, using in vitro splice assays in human embryonic kidney 293T cells. Reverse transcription-polymerase chain reaction and Sanger sequence analysis revealed splicing defects for 12 out of 19 variants. Four deep-intronic variants create pseudoexons or elongate the upstream exon. Furthermore, eight NCSS variants cause a partial deletion or skipping of one or more exons in messenger RNAs. Among the 12 variants, nine lead to premature stop codons and predicted truncated ABCA4 proteins. At least two deep-intronic variants affect splice enhancer and silencer motifs and, therefore, these conserved sequences should be carefully evaluated when predicting the outcome of NCSS and deep-intronic variants.

Where are the missing gene defects in inherited retinal disorders? Intronic and synonymous variants contribute at least to 4% of CACNA1F-mediated inherited retinal disorders.

Inherited retinal disorders (IRD) represent clinically and genetically heterogeneous diseases. To date, pathogenic variants have been identified in ~260 genes. Albeit that many genes are implicated in IRD, for 30-50% of the cases, the gene defect is unknown. These cases may be explained by novel gene defects, by overlooked structural variants, by variants in intronic, promoter or more distant regulatory regions, and represent synonymous variants of known genes contributing to the dysfunction of the respective proteins. Patients with one subgroup of IRD, namely incomplete congenital stationary night blindness (icCSNB), show a very specific phenotype. The major cause of this condition is the presence of a hemizygous pathogenic variant in CACNA1F. A comprehensive study applying direct Sanger sequencing of the gene-coding regions, exome and genome sequencing applied to a large cohort of patients with a clinical diagnosis of icCSNB revealed indeed that seven of the 189 CACNA1F-related cases have intronic and synonymous disease-causing variants leading to missplicing as validated by minigene approaches. These findings highlight that gene-locus sequencing may be a very efficient method in detecting disease-causing variants in clinically well-characterized patients with a diagnosis of IRD, like icCSNB.

Prevalence of ABCA4 Deep-Intronic Variants and Related Phenotype in An Unsolved "One-Hit" Cohort with Stargardt Disease.

We investigated the prevalence of reported deep-intronic variants in a French cohort of 70 patients with Stargardt disease harboring a monoallelic pathogenic variant on the exonic regions of ABCA4. Direct Sanger sequencing of selected intronic regions of ABCA4 was conducted. Complete phenotypic analysis and correlation with the genotype was performed in case a known intronic pathogenic variant was identified. All other variants found on the analyzed sequences were queried for minor allele frequency and possible pathogenicity by in silico predictions. The second mutated allele was found in 14 (20%) subjects. The three known deep-intronic variants found were c.5196+1137G>A in intron 36 (6 subjects), c.4539+2064C>T in intron 30 (4 subjects) and c.4253+43G>A in intron 28 (4 subjects). Even though the phenotype depends on the compound effect of the biallelic variants, a genotype-phenotype correlation suggests that the c.5196+1137G>A was mostly associated with a mild phenotype and the c.4539+2064C>T with a more severe one. A variable effect was instead associated with the variant c.4253+43G>A. In addition, two novel variants, c.768+508A>G and c.859-245_859-243delinsTGA never associated with Stargardt disease before, were identified and a possible splice defect was predicted in silico. Our study calls for a larger cohort analysis including targeted locus sequencing and 3D protein modeling to better understand phenotype-genotype correlations associated with deep-intronic changes and patients' selection for clinical trials.

In silico prioritization and further functional characterization of SPINK1 intronic variants.

BACKGROUND: SPINK1 (serine protease inhibitor, kazal-type, 1), which encodes human pancreatic secretory trypsin inhibitor, is one of the most extensively studied genes underlying chronic pancreatitis. Recently, based upon data from qualitative reverse transcription-PCR (RT-PCR) analyses of transfected HEK293T cells, we concluded that 24 studied SPINK1 intronic variants were not of pathological significance, the sole exceptions being two canonical splice site variants (i.e., c.87 + 1G > A and c.194 + 2T > C). Herein, we employed the splicing prediction tools included within the Alamut software suite to prioritize the 'non-pathological' SPINK1 intronic variants for further quantitative RT-PCR analysis. RESULTS: Although our results demonstrated the utility of in silico prediction in classifying and prioritizing intronic variants, we made two observations worth noting. First, we established that most of the prediction tools employed ignored the general rule that GC is a weaker donor splice site than the canonical GT site. This finding is potentially important because for a given disease gene, a GC variant donor splice site may be associated with a milder clinical manifestation. Second, the non-pathological c.194 + 13T > G variant was consistently predicted by different programs to generate a new and viable donor splice site, the prediction scores being comparable to those for the physiological c.194 + 2T donor splice site and even higher than those for the physiological c.87 + 1G donor splice site. We do however provide convincing in vitro evidence that the predicted donor splice site was not entirely spurious. CONCLUSIONS: Our findings, taken together, serve to emphasize the importance of functional analysis in helping to establish or refute the pathogenicity of specific intronic variants.

Characterization of Two Deep Intronic Variants on the ß-Globin Gene with Inconsistent Interpretations of Clinical Significance.

Sequence variants located in the introns of the ß-globin gene may affect the mRNA processing and cause ß-thalassemia (ß-thal). Sequence variants that change one of the invariant dinucleotides at the exon-intron boundaries may have fatal consequences for normal mRNA splicing. Intronic variants located far from obvious regulatory sequences can be more difficult to evaluate. There is a potential for misinterpretation of such sequence variants. Hence, thorough evaluation of patient data together with critical use of databases and in silico prediction tools are important. Here, we describe two rare sequence variants in the second intron of the ß-globin gene, HBB: c.316-70C>G and HBB: c.316-125A>G (NM_000518.4), both previously reported as variants causing ß-thal, and later as benign sequence variants. Due to the limited number of published cases and inconsistent interpretations, the significance of these sequence variants has been unclear. We have identified these two sequence variants in multiple individuals, alone and in a variety of combinations with other δ- and ß-globin defects, and we find no influence of the sequence variants on the phenotype.