Comparative analysis of in-silico tools in identifying pathogenic variants in dominant inherited retinal diseases.

Inherited retinal diseases (IRDs) are a group of rare genetic eye conditions that cause blindness. Despite progress in identifying genes associated with IRDs, improvements are necessary for classifying rare autosomal dominant (AD) disorders. AD diseases are highly heterogenous, with causal variants being restricted to specific amino acid changes within certain protein domains, making AD conditions difficult to classify. Here, we aim to determine the top-performing in-silico tools for predicting the pathogenicity of AD IRD variants. We annotated variants from ClinVar and benchmarked 39 variant classifier tools on IRD genes, split by inheritance pattern. Using area-under-the-curve (AUC) analysis, we determined the top-performing tools and defined thresholds for variant pathogenicity. Top-performing tools were assessed using genome sequencing on a cohort of participants with IRDs of unknown etiology. MutScore achieved the highest accuracy within AD genes, yielding an AUC of 0.969. When filtering for AD gain-of-function and dominant negative variants, BayesDel had the highest accuracy with an AUC of 0.997. Five participants with variants in NR2E3, RHO, GUCA1A, and GUCY2D were confirmed to have dominantly inherited disease based on pedigree, phenotype, and segregation analysis. We identified two uncharacterized variants in GUCA1A (c.428T>A, p.Ile143Thr) and RHO (c.631C>G, p.His211Asp) in three participants. Our findings support using a multi-classifier approach comprised of new missense classifier tools to identify pathogenic variants in participants with AD IRDs. Our results provide a foundation for improved genetic diagnosis for people with IRDs.

Loss-of-function variants in UBAP1L cause autosomal recessive retinal degeneration.

PURPOSE: Inherited retinal diseases (IRDs) are a group of monogenic conditions that can lead to progressive blindness. Their missing heritability is still considerable, due in part to the presence of disease genes that await molecular identification. The purpose of this work was to identify novel genetic associations with IRDs. METHODS: Patients underwent a comprehensive ophthalmological evaluation using standard-of-care tests, such as detailed retinal imaging (macular optical coherence tomography and short-wavelength fundus autofluorescence) and electrophysiological testing. Exome and genome sequencing, as well as computer-assisted data analysis were used for genotyping and detection of DNA variants. A minigene-driven splicing assay was performed to validate the deleterious effects of 1 of such variants. RESULTS: We identified 8 unrelated families from Hungary, the United States, Israel, and The Netherlands with members presenting with a form of autosomal recessive and nonsyndromic retinal degeneration, predominantly described as rod-cone dystrophy but also including cases of cone/cone-rod dystrophy. Age of disease onset was very variable, with some patients experiencing first symptoms during their fourth decade of life or later. Myopia greater than 5 diopters was present in 5 of 7 cases with available refractive data, and retinal detachment was reported in 2 cases. All ascertained patients carried biallelic loss-of-function variants in UBAP1L (HGNC: 40028), a gene with unknown function and with homologies to UBAP1, encoding a protein involved in ubiquitin metabolism. One of these pathogenic variants, the intronic NM_001163692.2:c.910-7G>A substitution, was identified in 5 unrelated families. Minigene-driven splicing assays in HEK293T cells confirmed that this DNA change is responsible for the creation of a new acceptor splice site, resulting in aberrant splicing. CONCLUSION: We identified UBAP1L as a novel IRD gene. Although its function is currently unknown, UBAP1L is almost exclusively expressed in photoreceptors and the retinal pigment epithelium, hence possibly explaining the link between pathogenic variants in this gene and an ocular phenotype.

Expanding the phenotypic and genotypic spectrum of patients with HGSNAT-related retinopathy.

BACKGROUND: Variants in HGSNAT have historically been associated with syndromic mucopolysaccharidosis type IIIC (MPSIIIC) but more recent studies demonstrate cases of HGSNAT-related non-syndromic retinitis pigmentosa. We describe and expand the genotypic and phenotypic spectrum of this disease. MATERIALS AND METHODS: This is a retrospective, observational, case series of 11 patients with pericentral retinitis pigmentosa due to variants in HGSNAT gene without a syndromic diagnosis of MPSIIIC. We reviewed ophthalmologic data extracted from medical records, genetic testing, color fundus photos, fundus autofluorescence (FAF), and optical coherence tomography (OCT). RESULTS: Of the 11 patients, the mean age was 52 years (range: 26-78). The mean age of ophthalmologic symptoms onset was 45 years (range: 15-72). The visual acuity varied from 20/20 to 20/80 (mean 20/30 median 20/20). We described five novel variants in HGSNAT: c.715del (p.Arg239Alafs *37), c.118 G>A (p.Asp40Asn), c.1218_1220delinsTAT, c.1297A>G (p.Asn433Asp), and c.1726 G>T (p.Gly576*). CONCLUSIONS: HGSNAT has high phenotypic heterogeneity. Data from our cohort showed that all patients who had at least one variant of c.1843 G>A (p.Ala615Thr) presented with the onset of ocular symptoms after the fourth decade of life. The two patients with onset of ocular symptoms before the fourth decade did not carry this variant. This may suggest that c.1843 G>A variant is associated with a later onset of retinopathy.

The genetic counselor workforce in inherited retinal disease clinics: a descriptive assessment.

BACKGROUND: Genetic counselors (GCs) have practiced in Inherited Retinal Disease (IRD) clinics for several decades. In this small subspecialty of genetic counseling, GCs are critical for patient understanding of genetic information, which can have prognostic, systemic, family planning and therapeutic implications. Recently, both access to genetic testing for IRDs and the number of genes associated with IRDs (>350) has increased dramatically. However, the practice models and roles of IRD GCs have not been previously described. MATERIALS AND METHODS: GCs working in academic IRD clinics were surveyed to assess their experience, clinical practices, and roles performed. The collected data was compared to the broader genetic counseling profession and to other specialties using publicly available data on GC professional practices. RESULTS: While roles of IRD GCs were overlapping with those of the overall genetic counseling profession, all survey respondents reported diverse roles that included both clinical and non-clinical duties, spending up to half their time on research and educational responsibilities. Most respondents (89%) felt that their clinic's MD to GC ratio was too high, while clinical load varied. IRD GCs report varying degrees of prior genetic counseling and ophthalmology-specific experience but unanimously desire additional subspecialty-specific training. CONCLUSIONS: This descriptive assessment of a small subspecialty suggests a need for growth in the number of GCs practicing in IRD clinics and could help to inform development of new GC positions in IRD centers. It also highlights the desire for additional GC-specific education and may be relevant to curriculum development within GC programs.

Vitelliform maculopathy in MELAS syndrome.

Purpose: We present a unique case of foveomacular vitelliform lesions in a patient with metabolic encephalomyopathy, lactic acidosis and stroke-like episodes (MELAS). Observations: After performing large panel next generation sequencing genetic testing, there was no likely alternative genetic etiology for vitelliform maculopathy in this patient. Conclusions and Importance: We present a rare case of a visually asymptomatic pediatric patient with MELAS and vitelliform maculopathy, which may be part of the spectrum of retinal manifestations in MELAS. Pediatric-onset vitelliform maculopathy in MELAS may be under-diagnosed due to its asymptomatic nature. Given the known risk of choroidal neovascularization in vitelliform maculopathy, it is important to identify these patients for proper surveillance.

Novel Pathogenic Mutations Identified from Whole-Genome Sequencing in Unsolved Cases of Patients Affected with Inherited Retinal Diseases.

Inherited retinal diseases (IRDs) are a diverse set of visual disorders that collectively represent a major cause of early-onset blindness. With the reduction in sequencing costs in recent years, whole-genome sequencing (WGS) is being used more frequently, particularly when targeted gene panels and whole-exome sequencing (WES) fail to detect pathogenic mutations in patients. In this study, we performed mutation screens using WGS for a cohort of 311 IRD patients whose mutations were undetermined. A total of nine putative pathogenic mutations in six IRD patients were identified, including six novel mutations. Among them, four were deep intronic mutations that affected mRNA splicing, while the other five affected protein-coding sequences. Our results suggested that the rate of resolution of unsolved cases via targeted gene panels and WES can be further enhanced with WGS; however, the overall improvement may be limited.

Systematic assessment of the contribution of structural variants to inherited retinal diseases.

Despite increasing success in determining genetic diagnosis for patients with inherited retinal diseases (IRDs), mutations in about 30% of the IRD cases remain unclear or unsettled after targeted gene panel or whole exome sequencing. In this study, we aimed to investigate the contributions of structural variants (SVs) to settling the molecular diagnosis of IRD with whole-genome sequencing (WGS). A cohort of 755 IRD patients whose pathogenic mutations remain undefined was subjected to WGS. Four SV calling algorithms including include MANTA, DELLY, LUMPY, and CNVnator were used to detect SVs throughout the genome. All SVs identified by any one of these four algorithms were included for further analysis. AnnotSV was used to annotate these SVs. SVs that overlap with known IRD-associated genes were examined with sequencing coverage, junction reads, and discordant read pairs. PCR followed by Sanger sequencing was used to further confirm the SVs and identify the breakpoints. Segregation of the candidate pathogenic alleles with the disease was performed when possible. In total, sixteen candidate pathogenic SVs were identified in sixteen families, including deletions and inversions, representing 2.1% of patients with previously unsolved IRDs. Autosomal dominant, autosomal recessive, and X-linked inheritance of disease-causing SVs were observed in 12 different genes. Among these, SVs in CLN3, EYS, PRPF31 were found in multiple families. Our study suggests that the contribution of SVs detected by short-read WGS is about 0.25% of our IRD patient cohort and is significantly lower than that of single nucleotide changes and small insertions and deletions.

Systematic assessment of the contribution of structural variants to inherited retinal diseases.

Despite increasing success in determining genetic diagnosis for patients with inherited retinal diseases (IRDs), mutations in about 30% of the IRD cases remain unclear or unsettled after targeted gene panel or whole exome sequencing. In this study, we aimed to investigate the contributions of structural variants (SVs) to settling the molecular diagnosis of IRD with whole-genome sequencing (WGS). A cohort of 755 IRD patients whose pathogenic mutations remain undefined were subjected to WGS. Four SV calling algorithms including include MANTA, DELLY, LUMPY and CNVnator were used to detect SVs throughout the genome. All SVs identified by any one of these four algorithms were included for further analysis. AnnotSV was used to annotate these SVs. SVs that overlap with known IRD-associated genes were examined with sequencing coverage, junction reads and discordant read pairs. Polymerase Chain Reaction (PCR) followed by Sanger sequencing was used to further confirm the SVs and identify the breakpoints. Segregation of the candidate pathogenic alleles with the disease was performed when possible. A total of 16 candidate pathogenic SVs were identified in 16 families, including deletions and inversions, representing 2.1% of patients with previously unsolved IRDs. Autosomal dominant, autosomal recessive and X-linked inheritance of disease-causing SVs were observed in 12 different genes. Among these, SVs in CLN3, EYS and PRPF31 were found in multiple families. Our study suggests that the contribution of SVs detected by short-read WGS is about 0.25% of our IRD patient cohort and is significantly lower than that of single nucleotide changes and small insertions and deletions.

A DOUBLE HYPERAUTOFLUORESCENT RING IN A 33-YEAR-OLD-FEMALE PATIENT.

PURPOSE: To describe the clinical phenotype and molecular diagnosis of a patient with atypical presentation of enhanced S-cone syndrome. METHODS: This is a case report of a patient who underwent best-corrected visual acuity, slit-lamp exam, fundus examination, autofluorescence, optical coherence tomography, kinetic perimetry, and full-field electroretinography. Genetic testing was performed via next-generation sequencing. RESULTS: A 33-year-old female patient presented with mild nyctalopia, but normal rod function measured by electroretinogram and foveoschisis on optical coherence tomography. She also presented a double hyperautofluorescent ring on autofluorescence. Genetic testing found a pathogenic variant c.925C>G (p.Arg309Gly) and a likely pathogenic variant c.299C>T (p.Arg77Trp) in NR2E3 gene. CONCLUSION: Enhanced S-cone syndrome may present without the pathognomonic findings of decreased rod function on electroretinogram, suggesting the importance of genetic testing in retinal diseases for diagnosis.

Signal Peptide Variants in Inherited Retinal Diseases: A Multi-Institutional Case Series.

Signal peptide (SP) mutations are an infrequent cause of inherited retinal diseases (IRDs). We report the genes currently associated with an IRD that possess an SP sequence and assess the prevalence of these variants in a multi-institutional retrospective review of clinical genetic testing records. The online databases, RetNet and UniProt, were used to determine which IRD genes possess a SP. A multicenter retrospective review was performed to retrieve cases of patients with a confirmed diagnosis of an IRD and a concurrent SP variant. In silico evaluations were performed with MutPred, MutationTaster, and the signal peptide prediction tool, SignalP 6.0. SignalP 6.0 was further used to determine the locations of the three SP regions in each gene: the N-terminal region, hydrophobic core, and C-terminal region. Fifty-six (56) genes currently associated with an IRD possess a SP sequence. Based on the records review, a total of 505 variants were present in the 56 SP-possessing genes. Six (1.18%) of these variants were within the SP sequence and likely associated with the patients' disease based on in silico predictions and clinical correlation. These six SP variants were in the CRB1 (early-onset retinal dystrophy), NDP (familial exudative vitreoretinopathy) (FEVR), FZD4 (FEVR), EYS (retinitis pigmentosa), and RS1 (X-linked juvenile retinoschisis) genes. It is important to be aware of SP mutations as an exceedingly rare cause of IRDs. Future studies will help refine our understanding of their role in each disease process and assess therapeutic approaches.

Retrograde Ret signaling controls sensory pioneer axon outgrowth.

The trafficking mechanisms and transcriptional targets downstream of long-range neurotrophic factor ligand/receptor signaling that promote axon growth are incompletely understood. Zebrafish carrying a null mutation in a neurotrophic factor receptor, Ret, displayed defects in peripheral sensory axon growth cone morphology and dynamics. Ret receptor was highly enriched in sensory pioneer neurons and Ret51 isoform was required for pioneer axon outgrowth. Loss-of-function of a cargo adaptor, Jip3, partially phenocopied Ret axonal defects, led to accumulation of activated Ret in pioneer growth cones, and reduced retrograde Ret51 transport. Jip3 and Ret51 were also retrogradely co-transported, ultimately suggesting Jip3 is a retrograde adapter of active Ret51. Finally, loss of Ret reduced transcription and growth cone localization of Myosin-X, an initiator of filopodial formation. These results show a specific role for Ret51 in pioneer axon growth, and suggest a critical role for long-range retrograde Ret signaling in regulating growth cone dynamics through downstream transcriptional changes.

In Vivo Time-Lapse Imaging in the Zebrafish Lateral Line: A Flexible, Open-Ended Research Project for an Undergraduate Neurobiology Laboratory Course.

The lateral line sensory system in fish detects movements in the water and allows fish to respond to predators, prey, and other stimuli. As the lateral line forms in the first two days of zebrafish development, axons extend caudally along the lateral surface of the fish, eventually forming synapses with hair cells of neuromasts. Growing lateral line axons are located superficially under the skin and can be labeled in living zebrafish using fluorescent protein expression. This system provides a relatively straightforward approach for in vivo time-lapse imaging of neuronal development in an undergraduate setting. Here we describe an upper-level neurobiology laboratory module in which students investigate aspects of axonal development in the zebrafish lateral line system. Students learn to handle and image living fish, collect time-lapse videos of moving mitochondria, and quantitatively measure mitochondrial dynamics by generating and analyzing kymographs of their movements. Energy demands may differ between axons with extending growth cones versus axons that have already reached their targets and are forming synapses. Since relatively little is known about this process in developing lateral line axons, students generate and test their own hypotheses regarding how mitochondrial dynamics may differ at two different time points in axonal development. Students also learn to incorporate into their analysis a powerful yet accessible quantitative tool, the kymograph, which is used to graph movement over time. After students measure and quantify dynamics in living fish at 1 and 2 days post fertilization, this module extends into independent projects, in which students can expand their studies in a number of different, inquiry-driven directions. The project can also be pared down for courses that wish to focus solely on the quantitative analysis (without fish handling), or vice versa. This research module provides a useful approach for the design of open-ended laboratory research projects that integrate the scientific process into undergraduate Biology courses, as encouraged by the AAAS and NSF Vision and Change Initiative.