ABCA4 c.6480-35A>G, a novel branchpoint variant associated with Stargardt disease

Introduction: Inherited retinal dystrophies (IRDs) can be caused by variants in more than 280 genes. The ATP-binding cassette transporter type A4 (ABCA4) gene is one of these genes and has been linked to Stargardt disease type 1 (STGD1), fundus flavimaculatus, cone-rod dystrophy (CRD), and pan-retinal CRD. Approximately 25% of the reported ABCA4 variants affect RNA splicing. In most cases, it is necessary to perform a functional assay to determine the effect of these variants. Methods: Whole genome sequencing (WGS) was performed in one Spanish proband with Stargardt disease. The putative pathogenicity of c.6480-35A>G on splicing was investigated both in silico and in vitro. The in silico approach was based on the deep-learning tool SpliceAI. For the in vitro approach we used a midigene splice assay in HEK293T cells, based on a previously established wild-type midigene (BA29) containing ABCA4 exons 46 to 48. Results: Through the analysis of WGS data, we identified two candidate variants in ABCA4 in one proband: a previously described deletion, c.699_768+342del (p.(Gln234Phefs*5)), and a novel branchpoint variant, c.6480-35A>G. Segregation analysis confirmed that the variants were in trans. For the branchpoint variant, SpliceAI predicted an acceptor gain with a high score (0.47) at position c.6480-47. A midigene splice assay in HEK293T cells revealed the inclusion of the last 47 nucleotides of intron 47 creating a premature stop codon and allowed to categorize the variant as moderately severe. Subsequent analysis revealed the presence of this variant as a second allele besides c.1958G>A p.(Arg653His) in an additional Spanish proband in a large cohort of IRD cases. Conclusion: A splice-altering effect of the branchpoint variant, confirmed by the midigene splice assay, along with the identification of this variant in a second unrelated individual affected with STGD, provides sufficient evidence to classify the variant as likely pathogenic. In addition, this research highlights the importance of studying non-coding regions and performing functional assays to provide a conclusive molecular diagnosis.

Combined Single Gene Testing and Genome Sequencing as an Effective Diagnostic Approach for Anophthalmia and Microphthalmia Patients

Anophthalmia and microphthalmia (A/M) are among the most severe congenital developmental eye disorders. Despite the advancements in genome screening technologies, more than half of A/M patients do not receive a molecular diagnosis. We included seven consanguineous families affected with A/M from Pakistani cohort and an unknown molecular basis. Single gene testing of FOXE3 was performed, followed by genome sequencing for unsolved probands in order to establish a genetic diagnosis for these families. All seven families were provided with a genetic diagnosis. The identified variants were all homozygous, classified as (likely) pathogenic and present in an A/M-associated gene. Targeted FOXE3 sequencing revealed two previously reported pathogenic FOXE3 variants in four families. In the remaining families, genome sequencing revealed a known pathogenic PXDN variant, a novel 13bp deletion in VSX2, and one novel deep intronic splice variant in PXDN. An in vitro splice assay was performed for the PXDN splice variant which revealed a severe splicing defect. Our study confirmed the utility of genome sequencing as a diagnostic tool for A/M-affected individuals. Furthermore, the identification of a novel deep intronic pathogenic variant in PXDN highlights the role of non-coding variants in A/M-disorders and the value of genome sequencing for the identification of this type of variants.

Optical genome mapping and revisiting short-read genome sequencing data reveal previously overlooked structural variants disrupting retinal disease-associated genes

PURPOSE

Structural variants (SVs) play an important role in inherited retinal diseases (IRD). Although the identification of SVs significantly improved upon the availability of genome sequencing, it is expected that involvement of SVs in IRDs is higher than anticipated. We revisited short-read genome sequencing data to enhance the identification of gene-disruptive SVs.

METHODS

Optical genome mapping was performed to improve SV detection in short-read genome sequencing-negative cases. In addition, reanalysis of short-read genome sequencing data was performed to improve the interpretation of SVs and to re-establish SV prioritization criteria.

RESULTS

In a monoallelic USH2A case, optical genome mapping identified a pericentric inversion (173 megabase), with 1 breakpoint disrupting USH2A. Retrospectively, the variant could be observed in genome sequencing data but was previously deemed false positive. Reanalysis of short-read genome sequencing data (427 IRD cases) was performed which yielded 30 pathogenic SVs affecting, among other genes, USH2A (n = 15), PRPF31 (n = 3), and EYS (n = 2). Eight of these (>25%) were overlooked during previous analyses.

CONCLUSION

Critical evaluation of our findings allowed us to re-establish and improve our SV prioritization and interpretation guidelines, which will prevent missing pathogenic events in future analyses. Our data suggest that more attention should be paid to SV interpretation and the current contribution of SVs in IRDs is still underestimated.

Cost-effective sequence analysis of 113 genes in 1,192 probands with retinitis pigmentosa and Leber congenital amaurosis

Introduction: Retinitis pigmentosa (RP) and Leber congenital amaurosis (LCA) are two groups of inherited retinal diseases (IRDs) where the rod photoreceptors degenerate followed by the cone photoreceptors of the retina. A genetic diagnosis for IRDs is challenging since >280 genes are associated with these conditions. While whole exome sequencing (WES) is commonly used by diagnostic facilities, the costs and required infrastructure prevent its global applicability. Previous studies have shown the cost-effectiveness of sequence analysis using single molecule Molecular Inversion Probes (smMIPs) in a cohort of patients diagnosed with Stargardt disease and other maculopathies. Methods: Here, we introduce a smMIPs panel that targets the exons and splice sites of all currently known genes associated with RP and LCA, the entire RPE65 gene, known causative deep-intronic variants leading to pseudo-exons, and part of the RP17 region associated with autosomal dominant RP, by using a total of 16,812 smMIPs. The RP-LCA smMIPs panel was used to screen 1,192 probands from an international cohort of predominantly RP and LCA cases. Results and discussion: After genetic analysis, a diagnostic yield of 56% was obtained which is on par with results from WES analysis. The effectiveness and the reduced costs compared to WES renders the RP-LCA smMIPs panel a competitive approach to provide IRD patients with a genetic diagnosis, especially in countries with restricted access to genetic testing.

Efficient Generation of Knock-In Zebrafish Models for Inherited Disorders Using CRISPR-Cas9 Ribonucleoprotein Complexes

CRISPR-Cas9-based genome-editing is a highly efficient and cost-effective method to generate zebrafish loss-of-function alleles. However, introducing patient-specific variants into the zebrafish genome with CRISPR-Cas9 remains challenging. Targeting options can be limited by the predetermined genetic context, and the efficiency of the homology-directed DNA repair pathway is relatively low. Here, we illustrate our efficient approach to develop knock-in zebrafish models using two previously variants associated with hereditary sensory deficits. We employ sgRNA-Cas9 ribonucleoprotein (RNP) complexes that are micro-injected into the first cell of fertilized zebrafish eggs together with an asymmetric, single-stranded DNA template containing the variant of interest. The introduction of knock-in events was confirmed by massive parallel sequencing of genomic DNA extracted from a pool of injected embryos. Simultaneous morpholino-induced blocking of a key component of the non-homologous end joining DNA repair pathway, Ku70, improved the knock-in efficiency for one of the targets. Our use of RNP complexes provides an improved knock-in efficiency as compared to previously published studies. Correct knock-in events were identified in 3-8% of alleles, and 30-45% of injected animals had the target variant in their germline. The detailed technical and procedural insights described here provide a valuable framework for the efficient development of knock-in zebrafish models.

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

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

The Impact of Modern Technologies on Molecular Diagnostic Success Rates, with a Focus on Inherited Retinal Dystrophy and Hearing Loss

The identification of pathogenic variants in monogenic diseases has been of interest to researchers and clinicians for several decades. However, for inherited diseases with extremely high genetic heterogeneity, such as hearing loss and retinal dystrophies, establishing a molecular diagnosis requires an enormous effort. In this review, we use these two genetic conditions as examples to describe the initial molecular genetic identification approaches, as performed since the early 90s, and subsequent improvements and refinements introduced over the years. Next, the history of DNA sequencing from conventional Sanger sequencing to high-throughput massive parallel sequencing, a.k.a. next-generation sequencing, is outlined, including their advantages and limitations and their impact on identifying the remaining genetic defects. Moreover, the development of recent technologies, also coined "third-generation" sequencing, is reviewed, which holds the promise to overcome these limitations. Furthermore, we outline the importance and complexity of variant interpretation in clinical diagnostic settings concerning the massive number of different variants identified by these methods. Finally, we briefly mention the development of novel approaches such as optical mapping and multiomics, which can help to further identify genetic defects in the near future.

Structural Variants Create New Topological-Associated Domains and Ectopic Retinal Enhancer-Gene Contact in Dominant Retinitis Pigmentosa

The cause of autosomal-dominant retinitis pigmentosa (adRP), which leads to loss of vision and blindness, was investigated in families lacking a molecular diagnosis. A refined locus for adRP on Chr17q22 (RP17) was delineated through genotyping and genome sequencing, leading to the identification of structural variants (SVs) that segregate with disease. Eight different complex SVs were characterized in 22 adRP-affected families with >300 affected individuals. All RP17 SVs had breakpoints within a genomic region spanning YPEL2 to LINC01476. To investigate the mechanism of disease, we reprogrammed fibroblasts from affected individuals and controls into induced pluripotent stem cells (iPSCs) and differentiated them into photoreceptor precursor cells (PPCs) or retinal organoids (ROs). Hi-C was performed on ROs, and differential expression of regional genes and a retinal enhancer RNA at this locus was assessed by qPCR. The epigenetic landscape of the region, and Hi-C RO data, showed that YPEL2 sits within its own topologically associating domain (TAD), rich in enhancers with binding sites for retinal transcription factors. The Hi-C map of RP17 ROs revealed creation of a neo-TAD with ectopic contacts between GDPD1 and retinal enhancers, and modeling of all RP17 SVs was consistent with neo-TADs leading to ectopic retinal-specific enhancer-GDPD1 accessibility. qPCR confirmed increased expression of GDPD1 and increased expression of the retinal enhancer that enters the neo-TAD. Altered TAD structure resulting in increased retinal expression of GDPD1 is the likely convergent mechanism of disease, consistent with a dominant gain of function. Our study highlights the importance of SVs as a genomic mechanism in unsolved Mendelian diseases.

A RIPOR2 in-frame deletion is a frequent and highly penetrant cause of adult-onset hearing loss

BACKGROUND

Hearing loss is one of the most prevalent disabilities worldwide, and has a significant impact on quality of life. The adult-onset type of the condition is highly heritable but the genetic causes are largely unknown, which is in contrast to childhood-onset hearing loss.

METHODS

Family and cohort studies included exome sequencing and characterisation of the hearing phenotype. Ex vivo protein expression addressed the functional effect of a DNA variant.

RESULTS

An in-frame deletion of 12 nucleotides in RIPOR2 was identified as a highly penetrant cause of adult-onset progressive hearing loss that segregated as an autosomal dominant trait in 12 families from the Netherlands. Hearing loss associated with the deletion in 63 subjects displayed variable audiometric characteristics and an average (SD) age of onset of 30.6 (14.9) years (range 0-70 years). A functional effect of the RIPOR2 variant was demonstrated by aberrant localisation of the mutant RIPOR2 in the stereocilia of cochlear hair cells and failure to rescue morphological defects in RIPOR2-deficient hair cells, in contrast to the wild-type protein. Strikingly, the RIPOR2 variant is present in 18 of 22 952 individuals not selected for hearing loss in the Southeast Netherlands.

CONCLUSION

Collectively, the presented data demonstrate that an inherited form of adult-onset hearing loss is relatively common, with potentially thousands of individuals at risk in the Netherlands and beyond, which makes it an attractive target for developing a (genetic) therapy.

Homozygous variants in KIAA1549, encoding a ciliary protein, are associated with autosomal recessive retinitis pigmentosa

BACKGROUND

Retinitis pigmentosa (RP) shows substantial genetic heterogeneity. It has been estimated that in approximately 60%-80% of RP cases, the genetic diagnosis can be found using whole exome sequencing (WES). In this study, the purpose was to identify causative variants in individuals with genetically unexplained retinal disease, which included one consanguineous family with two affected siblings and one case with RP.

METHODS

To identify the genetic defect, WES was performed in both probands, and clinical analysis was performed. To obtain insight into the function of KIAA1549 in photoreceptors, mRNA expression, knockdown and protein localisation studies were performed.

RESULTS

Through analysis of WES data, based on population allele frequencies, and in silico prediction tools, we identified a homozygous missense variant and a homozygous frameshift variant in KIAA1549 that segregate in two unrelated families. Kiaa1549 was found to localise at the connecting cilium of the photoreceptor cells and the synapses of the mouse retina. Both variants affect the long transcript of KIAA1549, which encodes a 1950 amino acid protein and shows prominent brain expression. The shorter transcript encodes a 734 amino acid protein with a high retinal expression and is affected by the identified missense variant. Strikingly, knockdown of the long transcript also leads to decreased expression of the short transcript likely explaining the non-syndromic retinal phenotype caused by the two variants targeting different transcripts.

CONCLUSION

In conclusion, our results underscore the causality of segregating variants in KIAA1549 for autosomal recessive RP. Moreover, our data indicate that KIAA1549 plays a role in photoreceptor function.

MPZL2, Encoding the Epithelial Junctional Protein Myelin Protein Zero-like 2, Is Essential for Hearing in Man and Mouse

In a Dutch consanguineous family with recessively inherited nonsyndromic hearing impairment (HI), homozygosity mapping combined with whole-exome sequencing revealed a MPZL2 homozygous truncating variant, c.72del (p.Ile24Metfs^∗22). By screening a cohort of phenotype-matched subjects and a cohort of HI subjects in whom WES had been performed previously, we identified two additional families with biallelic truncating variants of MPZL2. Affected individuals demonstrated symmetric, progressive, mild to moderate sensorineural HI. Onset of HI was in the first decade, and high-frequency hearing was more severely affected. There was no vestibular involvement. MPZL2 encodes myelin protein zero-like 2, an adhesion molecule that mediates epithelial cell-cell interactions in several (developing) tissues. Involvement of MPZL2 in hearing was confirmed by audiometric evaluation of Mpzl2-mutant mice. These displayed early-onset progressive sensorineural HI that was more pronounced in the high frequencies. Histological analysis of adult mutant mice demonstrated an altered organization of outer hair cells and supporting cells and degeneration of the organ of Corti. In addition, we observed mild degeneration of spiral ganglion neurons, and this degeneration was most pronounced at the cochlear base. Although MPZL2 is known to function in cell adhesion in several tissues, no phenotypes other than HI were found to be associated with MPZL2 defects. This indicates that MPZL2 has a unique function in the inner ear. The present study suggests that deleterious variants of Mplz2/MPZL2 affect adhesion of the inner-ear epithelium and result in loss of structural integrity of the organ of Corti and progressive degeneration of hair cells, supporting cells, and spiral ganglion neurons.