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.

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.

Whole genome sequencing for USH2A-associated disease reveals several pathogenic deep-intronic variants that are amenable to splice correction.

A significant number of individuals with a rare disorder such as Usher syndrome (USH) and (non-)syndromic autosomal recessive retinitis pigmentosa (arRP) remain genetically unexplained. Therefore, we assessed subjects suspected of USH2A-associated disease and no or mono-allelic USH2A variants using whole genome sequencing (WGS) followed by an improved pipeline for variant interpretation to provide a conclusive diagnosis. One hundred subjects were screened using WGS to identify causative variants in USH2A or other USH/arRP-associated genes. In addition to the existing variant interpretation pipeline, a particular focus was put on assessing splice-affecting properties of variants, both in silico and in vitro. Also structural variants were extensively addressed. For variants resulting in pseudoexon inclusion, we designed and evaluated antisense oligonucleotides (AONs) using minigene splice assays and patient-derived photoreceptor precursor cells. Biallelic variants were identified in 49 of 100 subjects, including novel splice-affecting variants and structural variants, in USH2A or arRP/USH-associated genes. Thirteen variants were shown to affect USH2A pre-mRNA splicing, including four deep-intronic USH2A variants resulting in pseudoexon inclusion, which could be corrected upon AON treatment. We have shown that WGS, combined with a thorough variant interpretation pipeline focused on assessing pre-mRNA splicing defects and structural variants, is a powerful method to provide subjects with a rare genetic condition, a (likely) conclusive genetic diagnosis. This is essential for the development of future personalized treatments and for patients to be eligible for such treatments.

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.

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.

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.

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.

PISTILLATA paralogs in Tarenaya hassleriana have diverged in interaction specificity.

BACKGROUND: Floral organs are specified by MADS-domain transcription factors that act in a combinatorial manner, as summarized in the (A)BCE model. However, this evolutionarily conserved model is in contrast to a remarkable amount of morphological diversity in flowers. One of the mechanisms suggested to contribute to this diversity is duplication of floral MADS-domain transcription factors. Although gene duplication is often followed by loss of one of the copies, sometimes both copies are retained. If both copies are retained they will initially be redundant, providing freedom for one of the paralogs to change function. Here, we examine the evolutionary fate and functional consequences of a transposition event at the base of the Brassicales that resulted in the duplication of the floral regulator PISTILLATA (PI), using Tarenaya hassleriana (Cleomaceae) as a model system. RESULTS: The transposition of a genomic region containing a PI gene led to two paralogs which are located at different positions in the genome. The original PI copy is syntenic in position with most angiosperms, whereas the transposed copy is syntenic with the PI genes in Brassicaceae. The two PI paralogs of T. hassleriana have very similar expression patterns. However, they may have diverged in function, as only one of these PI proteins was able to act heterologously in the first whorl of A. thaliana flowers. We also observed differences in protein complex formation between the two paralogs, and the two paralogs exhibit subtle differences in DNA-binding specificity. Sequence analysis indicates that most of the protein sequence divergence between the two T. hassleriana paralogs emerged in a common ancestor of the Cleomaceae and the Brassicaceae. CONCLUSIONS: We found that the PI paralogs in T. hassleriana have similar expression patterns, but may have diverged at the level of protein function. Data suggest that most protein sequence divergence occurred rapidly, prior to the origin of the Brassicaceae and Cleomaceae. It is tempting to speculate that the interaction specificities of the Brassicaceae-specific PI proteins are different compared to the PI found in other angiosperms. This could lead to PI regulating partly different genes in the Brassicaceae, and ultimately might result in change floral in morphology.

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.

Altered neural responsivity to food cues in relation to food preferences, but not appetite-related hormone concentrations after RYGB-surgery.

BACKGROUND: After Roux-en-Y gastric bypass (RYGB) surgery, patients report a shift in food preferences away from high-energy foods. OBJECTIVE: We aimed to elucidate the potential mechanisms underlying this shift in food preferences by assessing changes in neural responses to food pictures and odors before and after RYGB. Additionally, we investigated whether altered neural responsivity was associated with changes in plasma endocannabinoid and ghrelin concentrations. DESIGN: 19 RYGB patients (4 men; age 41 ± 10 years; BMI 41 ± 1 kg/m2 before; BMI 36 ± 1 kg/m2 after) participated in this study. Before and two months after RYGB surgery, they rated their food preferences using the Macronutrient and Taste Preference Ranking Task and BOLD fMRI responses towards pictures and odors of high-, and low-energy foods and non-food items were measured. Blood samples were taken to determine plasma endocannabinoid and ghrelin concentrations pre- and post-surgery. RESULTS: Patients demonstrated a shift in food preferences away from high-fat/sweet and towards low-energy/savory food products, which correlated with decreased superior parietal lobule responsivity to high-energy food odor and a reduced difference in precuneus responsivity to high-energy versus low-energy food pictures. In the anteroventral prefrontal cortex (superior frontal gyrus) the difference in deactivation towards high-energy versus non-food odors reduced. The precuneus was less deactivated in response to all cues. Plasma concentrations of anandamide were higher after surgery, while plasma concentrations of other endocannabinoids and ghrelin did not change. Alterations in appetite-related hormone concentrations did not correlate with changes in neural responsivity. CONCLUSIONS: RYGB leads to changed responsivity of the frontoparietal control network that orchestrates top-down control to high-energy food compared to low-energy food and non-food cues, rather than in reward related brain regions, in a satiated state. Together with correlations with the shift in food preference from high- to low-energy foods this indicates a possible role in new food preference formation.

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.

Comparative analysis of binding patterns of MADS-domain proteins in Arabidopsis thaliana.

BACKGROUND: Correct flower formation requires highly specific temporal and spatial regulation of gene expression. In Arabidopsis thaliana the majority of the master regulators that determine flower organ identity belong to the MADS-domain transcription factor family. The canonical DNA binding motif for this transcription factor family is the CArG-box, which has the consensus CC(A/T)6GG. However, so far, a comprehensive analysis of MADS-domain binding patterns has not yet been performed. RESULTS: Eight publicly available ChIP-seq datasets of MADS-domain proteins that regulate the floral transition and flower formation were analyzed. Surprisingly, the preferred DNA binding motif of each protein was a CArG-box with an NAA extension. Furthermore, motifs of other transcription factors were found in the vicinity of binding sites of MADS-domain transcription factors, suggesting that interaction of MADS-domain proteins with other transcription factors is important for target gene regulation. Finally, conservation of CArG-boxes between Arabidopsis ecotypes was assessed to obtain information about their evolutionary importance. CArG-boxes that fully matched the consensus were more conserved than other CArG-boxes, suggesting that the perfect CArG-box is evolutionary more important than other CArG-box variants. CONCLUSION: Our analysis provides detailed insight into MADS-domain protein binding patterns. The results underline the importance of an extended version of the CArG-box and provide a first view on evolutionary conservation of MADS-domain protein binding sites in Arabidopsis ecotypes.