Noncanonical splice-site mutations are an important cause of inherited diseases. Based on in vitro and stem-cell-based studies, some splice-site variants show a stronger splice defect than expected based on their predicted effects, suggesting that other sequence motifs influence the outcome. We investigated whether splice defects due to human-inherited-disease-associated variants in noncanonical splice-site sequences in ABCA4, DMD, and TMC1 could be rescued by strengthening the splice site on the other side of the exon. Noncanonical 5′- and 3′-splice-site variants were selected. Rescue variants were introduced based on an increase in predicted splice-site strength, and the effects of these variants were analyzed using in vitro splice assays in HEK293T cells. Exon skipping due to five variants in noncanonical splice sites of exons in ABCA4, DMD, and TMC1 could be partially or completely rescued by increasing the predicted strengths of the other splice site of the same exon. We named this mechanism “splicing interdependency”, and it is likely based on exon recognition by splicing machinery. Awareness of this interdependency is of importance in the classification of noncanonical splice-site variants associated with disease and may open new opportunities for treatments.

Missing heritability in human diseases represents a major challenge. Although whole-genome sequencing enables the analysis of coding and non-coding sequences, substantial costs and data storage requirements hamper its large-scale use to (re)sequence genes in genetically unsolved cases. The ABCA4 gene implicated in Stargardt disease (STGD1) has been studied extensively for 22 years, but thousands of cases remained unsolved. Therefore, single molecule molecular inversion probes were designed that enabled an automated and cost-effective sequence analysis of the complete 128-kb ABCA4 gene. Analysis of 1,054 unsolved STGD and STGD-like probands resulted in bi-allelic variations in 448 probands. Twenty-seven different causal deep-intronic variants were identified in 117 alleles. Based on in vitro splice assays, the 13 novel causal deep-intronic variants were found to result in pseudo-exon (PE) insertions (n=10) or exon elongations (n=3). 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 that were accompanied by flanking exon deletions. Structural variant analysis revealed 11 distinct deletions, two of which contained small inverted segments. Uniparental isodisomy of chromosome 1 was identified in one proband. Integrated complete gene sequencing combined with transcript analysis, identified pathogenic deep-intronic and structural variants in 26% of bi-allelic cases not solved previously by sequencing of coding regions. This strategy serves as a model study that can be applied to other inherited diseases in which only one or a few genes are involved in the majority of cases.