Mutations in SAMD7 cause autosomal-recessive macular dystrophy with or without cone dysfunction.

Sterile alpha motif domain containing 7 (SAMD7) is a component of the Polycomb repressive complex 1, which inhibits transcription of many genes, including those activated by the transcription factor Cone-Rod Homeobox (CRX). Here we report bi-allelic mutations in SAMD7 as a cause of autosomal-recessive macular dystrophy with or without cone dysfunction. Four of these mutations affect splicing, while another mutation is a missense variant that alters the repressive effect of SAMD7 on CRX-dependent promoter activity, as shown by in vitro assays. Immunostaining of human retinal sections revealed that SAMD7 is localized in the nuclei of both rods and cones, as well as in those of cells belonging to the inner nuclear layer. These results place SAMD7 as a gene crucial for human retinal function and demonstrate a significant difference in the role of SAMD7 between the human and the mouse retina.

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.

Autozygome-guided exome-first study in a consanguineous cohort with early-onset retinal disease uncovers an isolated RIMS2 phenotype and a retina-enriched RIMS2 isoform.

Leber congenital amaurosis (LCA) and early-onset retinal degeneration (EORD) are inherited retinal diseases (IRD) characterized by early-onset vision impairment. Herein, we studied 15 Saudi families by whole exome sequencing (WES) and run-of-homozygosity (ROH) detection via AutoMap in 12/15 consanguineous families. This revealed (likely) pathogenic variants in 11/15 families (73%). A potential founder variant was found in RPGRIP1. Homozygous pathogenic variants were identified in known IRD genes (ATF6, CRB1, CABP4, RDH12, RIMS2, RPGRIP1, SPATA7). We established genotype-driven clinical reclassifications for ATF6, CABP4, and RIMS2. Specifically, we observed isolated IRD in the individual with the novel RIMS2 variant, and we found a retina-enriched RIMS2 isoform conserved but not annotated in mouse. The latter illustrates potential different phenotypic consequences of pathogenic variants depending on the particular tissue/cell-type specific isoforms they affect. Lastly, a compound heterozygous genotype in GUCY2D in one non-consanguineous family was demonstrated, and homozygous variants in novel candidate genes ATG2B and RUFY3 were found in the two remaining consanguineous families. Reporting these genes will allow to validate them in other IRD cohorts. Finally, the missing heritability of the two unsolved IRD cases may be attributed to variants in non-coding regions or structural variants that remained undetected, warranting future WGS studies.

Expanding the genetic landscape of Usher syndrome type IV caused by pathogenic ARSG variants.

Usher syndrome (USH) is the most common cause of deafblindness. USH is autosomal recessively inherited and characterized by rod-cone dystrophy or retinitis pigmentosa (RP), often accompanied by sensorineural hearing loss. Variants in >15 genes have been identified as causative for clinically and genetically distinct subtypes. Among the ultra-rare and recently discovered genes is ARSG, coding for the lysosomal sulfatase Arylsulfatase G. This subtype was assigned as "USH IV" with a late onset of RP and usually late-onset progressive SNHL without vestibular involvement. Here, we describe nine new subjects and the clinical description of four cases with the USH IV phenotype bearing seven novel and two known pathogenic variants. Functional experiments indicated the complete loss of sulfatase enzymatic activity upon ectopic expression of mutated ARSG cDNA. Interestingly, we identified a homozygous missense variant, p.(Arg99His), previously described in dogs with neuronal ceroid lipofuscinosis. Our study expands the genetic landscape of ARSG-USH IV and the number of known subjects by more than 30%. These findings highlight that USH IV likely has been underdiagnosed and emphasize the need to test molecularly unresolved subjects with deafblindness syndrome. Finally, testing of ARSG should be considered for the genetic work-up of apparent isolated inherited retinal diseases.

ISOLATED MACULOPATHY AND MODERATE ROD-CONE DYSTROPHY REPRESENT THE MILDER END OF THE RDH12-RELATED RETINAL DYSTROPHY SPECTRUM.

PURPOSE: To describe an isolated maculopathy and an intermediate rod-cone dystrophy phenotype as the milder end of the RDH12-related retinal dystrophy spectrum. METHODS: Seven patients (17-34 years of age) underwent an extensive ophthalmic workup including psychophysical and electrophysiological testing and multimodal imaging. RESULTS: Three patients have isolated macular disease. Best-corrected visual acuity (BCVA) ranges from 20/125 to 20/40 with normal visual fields or only limited central, relative scotomata, and normal full-field ERGs. Both optical coherence tomography scans and autofluorescent imaging hint at relatively better-preserved foveal quality initially. An intermediate rod-cone phenotype in four patients is characterized by a central retinal dystrophy extending just beyond the vascular arcades, characteristic peripapillary sparing, and additional scattered atrophic patches. Again, foveal quality is initially better on optical coherence tomography scans. Best-corrected visual acuity ranges from counting fingers to 20/32. Goldmann visual fields vary from central scotomata to severe generalized abnormalities. ERGs range between mild and severe rod-cone dysfunction. Nine distinct RDH12 pathogenic variants, two of which are novel, are identified. CONCLUSION: The classic phenotype of RDH12-related early-onset retinal dystrophy is expanded to include an isolated maculopathy and intermediate dystrophy phenotype, characterized by its later onset and milder course with a fair visual potential until much later in life, emphasizing the phenotypic heterogeneity of RDH12-related retinopathy.

Correction: Mapping the genomic landscape of inherited retinal disease genes prioritizes genes prone to coding and noncoding copy-number variations.

The original version of this Article contained an error in the spelling of the author Anja K. Mayer, which was incorrectly given as Anja Kathrin Mayer. This has now been corrected in both the PDF and HTML versions of the Article.

Mapping the genomic landscape of inherited retinal disease genes prioritizes genes prone to coding and noncoding copy-number variations.

PurposePart of the hidden genetic variation in heterogeneous genetic conditions such as inherited retinal diseases (IRDs) can be explained by copy-number variations (CNVs). Here, we explored the genomic landscape of IRD genes listed in RetNet to identify and prioritize those genes susceptible to CNV formation.MethodsRetNet genes underwent an assessment of genomic features and of CNV occurrence in the Database of Genomic Variants and literature. CNVs identified in an IRD cohort were characterized using targeted locus amplification (TLA) on extracted genomic DNA.ResultsExhaustive literature mining revealed 1,345 reported CNVs in 81 different IRD genes. Correlation analysis between rankings of genomic features and CNV occurrence demonstrated the strongest correlation between gene size and CNV occurrence of IRD genes. Moreover, we identified and delineated 30 new CNVs in IRD cases, 13 of which are novel and three of which affect noncoding, putative cis-regulatory regions. Finally, the breakpoints of six complex CNVs were determined using TLA in a hypothesis-neutral manner.ConclusionWe propose a ranking of CNV-prone IRD genes and demonstrate the efficacy of TLA for the characterization of CNVs on extracted DNA. Finally, this IRD-oriented CNV study can serve as a paradigm for other genetically heterogeneous Mendelian diseases with hidden genetic variation.

Flexible, scalable, and efficient targeted resequencing on a benchtop sequencer for variant detection in clinical practice.

The release of benchtop next-generation sequencing (NGS) instruments has paved the way to implement the technology in clinical setting. The need for flexible, qualitative, and cost-efficient workflows is high. We used singleplex-PCR for highly efficient target enrichment, allowing us to reach the quality standards set in Sanger sequencing-based diagnostics. For the library preparation, a modified NexteraXT protocol was used, followed by sequencing on a MiSeq instrument. With an innovative pooling strategy, high flexibility, scalability, and cost-efficiency were obtained, independent of the availability of commercial kits. The approach was validated for ∼250 genes associated with monogenic disorders. An overall sensitivity (>99%) similar to Sanger sequencing was observed in combination with a positive predictive value of >98%. The distribution of coverage was highly uniform, guaranteeing a minimal number of gaps to be filled with alternative methods. ISO15189-accreditation was obtained for the workflow. A major asset of the singleplex PCR-based enrichment is that new targets can be easily implemented. Diagnostic laboratories have validated assays available ensuring that the proposed workflow can easily be adopted. Although our platform was optimized for constitutional variant detection of monogenic disease genes, it is now also used as a model for somatic mutation detection in acquired diseases.

saseR: Juggling offsets unlocks RNA-seq tools for fast and Scalable differential usage, Aberrant Splicing and Expression Retrieval.

RNA-seq data analysis relies on many different tools, each tailored to specific applications and coming with unique assumptions and restrictions. Indeed, tools for differential transcript usage, or diagnosing patients with rare diseases through splicing and expression outliers, either lack in performance, discard information, or do not scale to massive data compendia. Here, we show that replacing the normalisation offsets unlocks bulk RNA-seq workflows for scalable differential usage, aberrant splicing and expression analyses. Our method, saseR, is much faster than state-of-the-art methods, dramatically outperforms these to detect aberrant splicing, and provides a single workflow for various short- and long-read RNA-seq applications.

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.

CEP162 deficiency causes human retinal degeneration and reveals a dual role in ciliogenesis and neurogenesis.

Defects in primary or motile cilia result in a variety of human pathologies, and retinal degeneration is frequently associated with these so-called ciliopathies. We found that homozygosity for a truncating variant in CEP162, a centrosome and microtubule-associated protein required for transition zone assembly during ciliogenesis and neuronal differentiation in the retina, caused late-onset retinitis pigmentosa in 2 unrelated families. The mutant CEP162-E646R*5 protein was expressed and properly localized to the mitotic spindle, but it was missing from the basal body in primary and photoreceptor cilia. This impaired recruitment of transition zone components to the basal body and corresponded to complete loss of CEP162 function at the ciliary compartment, reflected by delayed formation of dysmorphic cilia. In contrast, shRNA knockdown of Cep162 in the developing mouse retina increased cell death, which was rescued by expression of CEP162-E646R*5, indicating that the mutant retains its role for retinal neurogenesis. Human retinal degeneration thus resulted from specific loss of the ciliary function of CEP162.

Mild Leber hereditary optic neuropathy (LHON) in a Western European family due to the rare Asian m.14502T>C variant in the MT-ND6 gene.

BACKGROUND: Leber hereditary optic neuropathy (LHON) is a mitochondrial neurodegenerative disease. The majority (>90%) is related to three primary mitochondrial DNA (mtDNA) variants: ND1 m.3460G>A, ND4 m.11778G>A and ND6 m.14484T>C. The remaining 10% is associated with >40 secondary variants with variable penetrance and incidence between different ethnic backgrounds. MATERIALS AND METHODS: Five sisters underwent an extensive ophthalmic workup including psychophysical, electrophysiological, multimodal brain imaging, biochemical testing and molecular screening. MT-ND6 protein modelling was performed. RESULTS: A 23-year-old woman presented with acute central visual loss to counting fingers in the right eye. She developed a central visual field scotoma, severe color vision deficiencies and impaired pattern visual evoked responses. Progressive optic atrophy ensued. The left eye was unremarkable, except for borderline thinning of the temporal retinal nerve fiber layer. Alcohol use and passive smoking were noted. MtDNA analysis revealed a rare variant, m.14502T>C in MT-ND6, exclusively known to cause optic neuropathy in an Asian population. Three sisters of the proband, two of whom reported tobacco and alcohol abuse, had bilateral temporal optic disc pallor without functional impact. A fourth non-smoker sister had a completely normal eye exam. CONCLUSIONS: The rare Asian m.14502T>C variant in the MT-ND6 gene was linked to a mild LHON phenotype in a Western European family. Penetrance in this family was likely triggered by alcohol and tobacco abuse. A full mtDNA sequencing is warranted in the case of high clinical suspicion of LHON when mutation analysis for the three common pathogenic variants is negative.

Expanding the clinical spectrum and management of Traboulsi syndrome: report on two siblings homozygous for a novel pathogenic variant in ASPH.

BACKGROUND: Traboulsi syndrome is a very rare, syndromic form of ectopia lentis that is potentially sight-threatening at a young age. It is characterized by typical facial, skeletal and ocular signs. MATERIALS AND METHODS: Two siblings, born to consanguineous parents, with a clinical phenotype consistent with Traboulsi syndrome, underwent extensive ophthalmic imaging and exome-based genetic testing. Both were treated with unilateral clear lens extraction via a limbal approach. RESULTS: Two siblings, one male and one female, presented with systemic and ocular features consistent with Traboulsi syndrome. Lens subluxation was present in all 4fouraffected eyes, and spontaneous subconjunctival bleb formation was detected in one eye. This eye also showed evidence of keratoconus-related corneal thinning. The clinical diagnosis of Traboulsi syndrome was confirmed molecularly. A homozygous, novel, pathogenic nonsense variant was identified in exon 25 of the ASPH gene: c.2181_2183dup, p.(Val727_Trp728insTer). Excellent visual outcomes following clear lens extraction and postoperative rigid gas-permeable contact lens fitting were obtained. CONCLUSIONS: We expanded the genetic spectrum of Traboulsi syndrome with a novel frameshift variant in the ASPH gene. We showed that lensectomy followed by gas-permeable contact lenses is an efficient therapeutic approach to treat lens subluxation in Traboulsi syndrome. However, lifelong follow-up is crucial to avoid (late) postoperative complications.

Long-Read Sequencing to Unravel Complex Structural Variants of CEP78 Leading to Cone-Rod Dystrophy and Hearing Loss.

Inactivating variants as well as a missense variant in the centrosomal CEP78 gene have been identified in autosomal recessive cone-rod dystrophy with hearing loss (CRDHL), a rare syndromic inherited retinal disease distinct from Usher syndrome. Apart from this, a complex structural variant (SV) implicating CEP78 has been reported in CRDHL. Here we aimed to expand the genetic architecture of typical CRDHL by the identification of complex SVs of the CEP78 region and characterization of their underlying mechanisms. Approaches used for the identification of the SVs are shallow whole-genome sequencing (sWGS) combined with quantitative polymerase chain reaction (PCR) and long-range PCR, or ExomeDepth analysis on whole-exome sequencing (WES) data. Targeted or whole-genome nanopore long-read sequencing (LRS) was used to delineate breakpoint junctions at the nucleotide level. For all SVs cases, the effect of the SVs on CEP78 expression was assessed using quantitative PCR on patient-derived RNA. Apart from two novel canonical CEP78 splice variants and a frameshifting single-nucleotide variant (SNV), two SVs affecting CEP78 were identified in three unrelated individuals with CRDHL: a heterozygous total gene deletion of 235 kb and a partial gene deletion of 15 kb in a heterozygous and homozygous state, respectively. Assessment of the molecular consequences of the SVs on patient's materials displayed a loss-of-function effect. Delineation and characterization of the 15-kb deletion using targeted LRS revealed the previously described complex CEP78 SV, suggestive of a recurrent genomic rearrangement. A founder haplotype was demonstrated for the latter SV in cases of Belgian and British origin, respectively. The novel 235-kb deletion was delineated using whole-genome LRS. Breakpoint analysis showed microhomology and pointed to a replication-based underlying mechanism. Moreover, data mining of bulk and single-cell human and mouse transcriptional datasets, together with CEP78 immunostaining on human retina, linked the CEP78 expression domain with its phenotypic manifestations. Overall, this study supports that the CEP78 locus is prone to distinct SVs and that SV analysis should be considered in a genetic workup of CRDHL. Finally, it demonstrated the power of sWGS and both targeted and whole-genome LRS in identifying and characterizing complex SVs in patients with ocular diseases.

Prevalence of Germline Pathogenic Variants in Cancer Predisposing Genes in Czech and Belgian Pancreatic Cancer Patients.

(1) Background: The proportion and spectrum of germline pathogenic variants (PV) associated with an increased risk for pancreatic ductal adenocarcinoma (PDAC) varies among populations. (2) Methods: We analyzed 72 Belgian and 226 Czech PDAC patients by multigene panel testing. The prevalence of pathogenic variants (PV) in relation to personal/family cancer history were evaluated. PDAC risks were calculated using both gnomAD-NFE and population-matched controls. (3) Results: In 35/298 (11.7%) patients a PV in an established PDAC-predisposition gene was found. BRCA1/2 PV conferred a high risk in both populations, ATM and Lynch genes only in the Belgian subgroup. PV in other known PDAC-predisposition genes were rarer. Interestingly, a high frequency of CHEK2 PV was observed in both patient populations. PV in PDAC-predisposition genes were more frequent in patients with (i) multiple primary cancers (12/38; 32%), (ii) relatives with PDAC (15/56; 27%), (iii) relatives with breast/ovarian/colorectal cancer or melanoma (15/86; 17%) but more rare in sporadic PDAC (5/149; 3.4%). PV in homologous recombination genes were associated with improved overall survival (HR = 0.51; 95% CI 0.34-0.77). (4) Conclusions: Our analysis emphasizes the value of multigene panel testing in PDAC patients, especially in individuals with a positive family cancer history, and underlines the importance of population-matched controls for risk assessment.

Accurate detection and quantification of epigenetic and genetic second hits in BRCA1 and BRCA2-associated hereditary breast and ovarian cancer reveals multiple co-acting second hits.

BACKGROUND: This study characterizes the second hit spectrum in BRCA1 and BRCA2-associated breast and ovarian cancers at both gene loci to investigate if second hit mechanisms are mutually exclusive or able to coincide within the same tumor. METHODS: Loss of heterozygosity, somatic point mutations and copy number alterations along with promoter methylation were studied in 56 breast and 15 ovarian cancers from BRCA1 and BRCA2 germline mutation carriers. A mathematical methodology was introduced to quantify the tumor cell population carrying a second hit. RESULTS: Copy neutral LOH was the most prevalent LOH mechanism in this cohort (BC 69%, OC 67%). However, only 36% of BC and 47% of OC showed LOH in all cancerous cells. Somatic intragenic deletions and methylated subclones were also found in combination with (partial) loss of heterozygosity. Unequivocal deleterious somatic point mutations were not identified in this cohort. CONCLUSION: Different mechanisms inactivating the wild type allele are present within the same tumor sample at various extents. Results indicate that BRCA1/2-linked breast and ovarian cancer cells are predominantly characterized by LOH, but harbor a complex combination of second hits at various frequencies.

Evaluation of relative quantification of alternatively spliced transcripts using droplet digital PCR.

INTRODUCTION: For the relative quantification of isoform expression, RT-qPCR has been the gold standard for over a decade. More recently, digital PCR is becoming widely implemented, as it is promised to be more accurate, sensitive and less affected by inhibitors, without the need for standard curves. In this study we evaluated RT-qPCR versus RT-droplet digital PCR (ddPCR) for the relative quantification of isoforms in controls and carriers of the splice site mutation BRCA1 c.212+3A>G, associated with increased expression of several isoforms. MATERIALS AND METHODS: RNA was extracted from EBV cell lines of controls and heterozygous BRCA1 c.212+3A>G carriers. Transcript-specific plasmids were available to determine the efficiency, precision, reproducibility and accuracy of each method. RESULTS: Both ddPCR and RT-qPCR were able to accurately quantify all targets and showed the same LOB, LOD and LOQ; also precision and reproducibility were similar. Both techniques have the same dynamic range and linearity at biologically relevant template concentrations. However, a significantly higher cost and workload was required for ddPCR experiments. CONCLUSIONS: Our study recognizes the potential and validity of digital PCR but shows the value of a highly optimized qPCR for the relative quantification of isoforms. Cost efficiency and simplicity turned out to be better for RT-qPCR.