Structure-based network analysis predicts pathogenic variants in human proteins associated with inherited retinal disease.

Advances in gene sequencing technologies have accelerated the identification of genetic variants, but better tools are needed to understand which are causal of disease. This would be particularly useful in fields where gene therapy is a potential therapeutic modality for a disease-causing variant such as inherited retinal disease (IRD). Here, we apply structure-based network analysis (SBNA), which has been successfully utilized to identify variant-constrained amino acid residues in viral proteins, to identify residues that may cause IRD if subject to missense mutation. SBNA is based entirely on structural first principles and is not fit to specific outcome data, which makes it distinct from other contemporary missense prediction tools. In 4 well-studied human disease-associated proteins (BRCA1, HRAS, PTEN, and ERK2) with high-quality structural data, we find that SBNA scores correlate strongly with deep mutagenesis data. When applied to 47 IRD genes with available high-quality crystal structure data, SBNA scores reliably identified disease-causing variants according to phenotype definitions from the ClinVar database. Finally, we applied this approach to 63 patients at Massachusetts Eye and Ear (MEE) with IRD but for whom no genetic cause had been identified. Untrained models built using SBNA scores and BLOSUM62 scores for IRD-associated genes successfully predicted the pathogenicity of novel variants (AUC = 0.851), allowing us to identify likely causative disease variants in 40 IRD patients. Model performance was further augmented by incorporating orthogonal data from EVE scores (AUC = 0.927), which are based on evolutionary multiple sequence alignments. In conclusion, SBNA can used to successfully identify variants as causal of disease in human proteins and may help predict variants causative of IRD in an unbiased fashion.

Exome copy number variant detection, analysis, and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and, with new innovative methods, can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the Genomics Research to Elucidate the Genetics of Rare Diseases consortium and analyzed using the seqr platform. The addition of CNV detection to exome analysis identified causal CNVs for 171 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb. The causal CNVs consisted of 140 deletions, 15 duplications, 3 suspected complex structural variants (SVs), 3 insertions, and 10 complex SVs, the latter two groups being identified by orthogonal confirmation methods. To classify CNV variant pathogenicity, we used the 2020 American College of Medical Genetics and Genomics/ClinGen CNV interpretation standards and developed additional criteria to evaluate allelic and functional data as well as variants on the X chromosome to further advance the framework. We interpreted 151 CNVs as likely pathogenic/pathogenic and 20 CNVs as high-interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher-resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

Coding and non-coding variants in the ciliopathy gene CFAP410 cause early-onset non-syndromic retinal degeneration.

Inherited retinal degenerations are blinding genetic disorders characterized by high genetic and phenotypic heterogeneity. The implementation of next-generation sequencing in routine diagnostics, together with advanced clinical phenotyping including multimodal retinal imaging, have contributed to the increase of reports describing novel genotype-phenotype associations and phenotypic expansions. In this study, we describe sixteen families with early-onset non-syndromic retinal degenerations in which affected probands carried rare bi-allelic variants in CFAP410, a ciliary gene previously associated with syndromic recessive Jeune syndrome. The most common retinal phenotypes were cone-rod and rod-cone dystrophies, but the clinical presentations were unified by their early onset as well as the severe impact on central visual function. Twelve variants were detected (three pathogenic, seven likely pathogenic, two of uncertain significance), eight of which were novel. One deep intronic change, c.373+91A>G, led to the creation of a cryptic splice acceptor site in intron four, followed by the inclusion of a 200- base pair pseudoexon and subsequent premature stop codon formation. To our knowledge this is the first likely pathogenic deep-intronic variant identified in this gene. Meta-analysis of all published and novel CFAP410 variants revealed no clear correlation between the severity of the CFAP410-associated phenotypes and the identified causal variants. This is supported by the fact that the frequently encountered missense variant p.(Arg73Pro), often found in syndromic cases, was also associated with non-syndromic retinal degeneration. This study expands the current knowledge of CFAP410-associated ciliopathy by enriching its mutational landscape and supports its association with non-syndromic retinal degeneration.

Heterozygous loss-of-function SMC3 variants are associated with variable growth and developmental features.

Heterozygous missense variants and in-frame indels in SMC3 are a cause of Cornelia de Lange syndrome (CdLS), marked by intellectual disability, growth deficiency, and dysmorphism, via an apparent dominant-negative mechanism. However, the spectrum of manifestations associated with SMC3 loss-of-function variants has not been reported, leading to hypotheses of alternative phenotypes or even developmental lethality. We used matchmaking servers, patient registries, and other resources to identify individuals with heterozygous, predicted loss-of-function (pLoF) variants in SMC3, and analyzed population databases to characterize mutational intolerance in this gene. Here, we show that SMC3 behaves as an archetypal haploinsufficient gene: it is highly constrained against pLoF variants, strongly depleted for missense variants, and pLoF variants are associated with a range of developmental phenotypes. Among 14 individuals with SMC3 pLoF variants, phenotypes were variable but coalesced on low growth parameters, developmental delay/intellectual disability, and dysmorphism, reminiscent of atypical CdLS. Comparisons to individuals with SMC3 missense/in-frame indel variants demonstrated an overall milder presentation in pLoF carriers. Furthermore, several individuals harboring pLoF variants in SMC3 were nonpenetrant for growth, developmental, and/or dysmorphic features, and some had alternative symptomatologies with rational biological links to SMC3. Analyses of tumor and model system transcriptomic data and epigenetic data in a subset of cases suggest that SMC3 pLoF variants reduce SMC3 expression but do not strongly support clustering with functional genomic signatures of typical CdLS. Our finding of substantial population-scale LoF intolerance in concert with variable growth and developmental features in subjects with SMC3 pLoF variants expands the scope of cohesinopathies, informs on their allelic architecture, and suggests the existence of additional clearly LoF-constrained genes whose disease links will be confirmed only by multilayered genomic data paired with careful phenotyping.

Generation of a human induced pluripotent stem cell line (OGIi001) from peripheral blood mononuclear cells of a healthy male donor.

We have successfully derived a novel human induced pluripotent stem cell (hiPSC) line using non-integrative Sendai virus. This hiPSC line was generated from a healthy male adult donor, aged 55, and subjected to thorough characterization and extensive quality control. The analysis confirmed the expression of undifferentiated stem cell markers, demonstrated the ability to differentiate into the three germ layers, and revealed the absence of any chromosomal abnormalities.

WIDEFIELD SWEPT-SOURCE OPTICAL COHERENCE TOMOGRAPHY ANGIOGRAPHY FINDINGS IN WAGNER SYNDROME.

PURPOSE: To describe novel clinical and angiographic findings in Wagner syndrome. METHODS: A retrospective case series of three related patients with Wagner syndrome. Patients underwent standard optical coherence tomography (OCT), B-scan ultrasonography, and fluorescein angiography in addition to wide field swept-source OCT angiography (WF SS-OCTA) (PLEX Elite 9000, Carl Zeiss Meditec Inc). Patients underwent genetic testing for a panel of hereditary vitreoretinopathies. RESULTS: Three related patients with Wagner syndrome were identified. All were found to have prominent vitreous strands, abnormal vitreoretinal adhesions, peripheral retinal holes, and varying degrees of myopia. A mid-peripheral tractional ridge was identified in all six eyes. All patients were positive for a known pathologic intron variant in the VCAN gene (4004-5T-A). Wide field swept-source OCT angiography (12 mm × 12 mm) was performed in two patients and demonstrated perivascular capillary loss in the superficial capillary plexus along the arcades bilaterally. One patient demonstrated associated retinal atrophy within the area of capillary loss. The capillary loss extended beyond the margin of retinal atrophy. CONCLUSION: The unusual finding of a mid-peripheral tractional ridge of the retina associated with myopia led to a genetic diagnosis of Wagner syndrome. Widefield swept-source OCT angiography demonstrated a novel feature of perivascular loss of the superficial retinal capillary plexus. This result suggests that vitreous traction may cause localized microvasculature dysfunction and subsequent retinal atrophy in Wagner syndrome. This is the first known evaluation of Wagner syndrome using OCT angiography.

Heterozygous loss-of-function SMC3 variants are associated with variable and incompletely penetrant growth and developmental features.

Heterozygous missense variants and in-frame indels in SMC3 are a cause of Cornelia de Lange syndrome (CdLS), marked by intellectual disability, growth deficiency, and dysmorphism, via an apparent dominant-negative mechanism. However, the spectrum of manifestations associated with SMC3 loss-of-function variants has not been reported, leading to hypotheses of alternative phenotypes or even developmental lethality. We used matchmaking servers, patient registries, and other resources to identify individuals with heterozygous, predicted loss-of-function (pLoF) variants in SMC3, and analyzed population databases to characterize mutational intolerance in this gene. Here, we show that SMC3 behaves as an archetypal haploinsufficient gene: it is highly constrained against pLoF variants, strongly depleted for missense variants, and pLoF variants are associated with a range of developmental phenotypes. Among 13 individuals with SMC3 pLoF variants, phenotypes were variable but coalesced on low growth parameters, developmental delay/intellectual disability, and dysmorphism reminiscent of atypical CdLS. Comparisons to individuals with SMC3 missense/in-frame indel variants demonstrated a milder presentation in pLoF carriers. Furthermore, several individuals harboring pLoF variants in SMC3 were nonpenetrant for growth, developmental, and/or dysmorphic features, some instead having intriguing symptomatologies with rational biological links to SMC3 including bone marrow failure, acute myeloid leukemia, and Coats retinal vasculopathy. Analyses of transcriptomic and epigenetic data suggest that SMC3 pLoF variants reduce SMC3 expression but do not result in a blood DNA methylation signature clustering with that of CdLS, and that the global transcriptional signature of SMC3 loss is model-dependent. Our finding of substantial population-scale LoF intolerance in concert with variable penetrance in subjects with SMC3 pLoF variants expands the scope of cohesinopathies, informs on their allelic architecture, and suggests the existence of additional clearly LoF-constrained genes whose disease links will be confirmed only by multi-layered genomic data paired with careful phenotyping.

Exome copy number variant detection, analysis and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and with new innovative methods can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the GREGoR consortium. Each family's CNV data was analyzed using the seqr platform and candidate CNVs classified using the 2020 ACMG/ClinGen CNV interpretation standards. We developed additional evidence criteria to address situations not covered by the current standards. The addition of CNV calling to exome analysis identified causal CNVs for 173 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb with estimates that 44% would not have been detected by standard chromosomal microarrays. The causal CNVs consisted of 141 deletions, 15 duplications, 4 suspected complex structural variants (SVs), 3 insertions and 10 complex SVs, the latter two groups being identified by orthogonal validation methods. We interpreted 153 CNVs as likely pathogenic/pathogenic and 20 CNVs as high interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

RP2 X-LINKED RETINITIS PIGMENTOSA CARRIER STATE PRESENTING WITH VASCULAR LEAKAGE AND UNILATERAL MACULAR ATROPHY.

PURPOSE: We describe the unusual clinical presentation of a 33-year-old woman subsequently identified as a carrier of RP2-associated X-linked retinitis pigmentosa. METHODS: Case report. RESULTS: A 33-year-old woman without a known family history of retinal disease presented with unilateral reduced visual acuity and central scotoma in the left eye. Examination showed underlying macular atrophy in the left eye and a bilateral tapetal-like reflex. Full-field electroretinogram was abnormal in the left eye but normal in the right eye. Notable findings on wide-field imaging included bilateral peripheral vascular leakage on fluorescein angiography and a bilaterally symmetric radial pattern of hyperfluorescence on fundus autofluorescence. Genetic testing demonstrated a pathogenic variant in the gene RP2 confirming that she was a carrier of X-linked retinitis pigmentosa. CONCLUSION: We describe clinical features of the carrier state of RP2-XLRP and expand potential findings to include peripheral vascular leakage. This case highlights the importance of awareness of the carrier state, particularly if a family history cannot be provided.

Structure-based network analysis predicts mutations associated with inherited retinal disease.

With continued advances in gene sequencing technologies comes the need to develop better tools to understand which mutations cause disease. Here we validate structure-based network analysis (SBNA)1,2 in well-studied human proteins and report results of using SBNA to identify critical amino acids that may cause retinal disease if subject to missense mutation. We computed SBNA scores for genes with high-quality structural data, starting with validating the method using 4 well-studied human disease-associated proteins. We then analyzed 47 inherited retinal disease (IRD) genes. We compared SBNA scores to phenotype data from the ClinVar database and found a significant difference between benign and pathogenic mutations with respect to network score. Finally, we applied this approach to 65 patients at Massachusetts Eye and Ear (MEE) who were diagnosed with IRD but for whom no genetic cause was found. Multivariable logistic regression models built using SBNA scores for IRD-associated genes successfully predicted pathogenicity of novel mutations, allowing us to identify likely causative disease variants in 37 patients with IRD from our clinic. In conclusion, SBNA can be meaningfully applied to human proteins and may help predict mutations causative of IRD.

Natural history of retinitis pigmentosa based on genotype, vitamin A/E supplementation, and an electroretinogram biomarker.

BACKGROUNDA randomized clinical trial from 1984 to 1992 indicated that vitamin A supplementation had a beneficial effect on the progression of retinitis pigmentosa (RP), while vitamin E had an adverse effect.METHODSSequencing of banked DNA samples from that trial provided the opportunity to determine whether certain genotypes responded preferentially to vitamin supplementation.RESULTSThe genetic solution rate was 587 out of 765 (77%) of sequenced samples. Combining genetic solutions with electroretinogram outcomes showed that there were systematic differences in severity and progression seen among different genetic subtypes of RP, extending findings made for USH2A, RHO, RPGR, PRPF31, and EYS. Baseline electroretinogram 30-Hz flicker implicit time was an independent, strong predictor of progression rate. Using additional data and baseline implicit time as a predictor, the deleterious effect of vitamin E was still present. Surprisingly, the effect of vitamin A progression in the cohort as a whole was not detectable, with or without data from subsequent trials. Subgroup analyses are also discussed.CONCLUSIONOverall, genetic subtype and implicit time have significant predictive power for a patient's rate of progression, which is useful prognostically. While vitamin E supplementation should still be avoided, these data do not support a generalized neuroprotective effect of vitamin A for all types of RP.TRIAL REGISTRATIONClinicalTrials.gov NCT00000114, NCT00000116, and NCT00346333.FUNDINGFoundation Fighting Blindness and the National Eye Institute: R01 EY012910, R01 EY031036, R01 EY026904, and P30 EY014104.

Perspectives of Rare Disease Experts on Newborn Genome Sequencing.

Importance: Newborn genome sequencing (NBSeq) can detect infants at risk for treatable disorders currently undetected by conventional newborn screening. Despite broad stakeholder support for NBSeq, the perspectives of rare disease experts regarding which diseases should be screened have not been ascertained. Objective: To query rare disease experts about their perspectives on NBSeq and which gene-disease pairs they consider appropriate to evaluate in apparently healthy newborns. Design, Setting, and Participants: This survey study, designed between November 2, 2021, and February 11, 2022, assessed experts' perspectives on 6 statements related to NBSeq. Experts were also asked to indicate whether they would recommend including each of 649 gene-disease pairs associated with potentially treatable conditions in NBSeq. The survey was administered between February 11 and September 23, 2022, to 386 experts, including all 144 directors of accredited medical and laboratory genetics training programs in the US. Exposures: Expert perspectives on newborn screening using genome sequencing. Main Outcomes and Measures: The proportion of experts indicating agreement or disagreement with each survey statement and those who selected inclusion of each gene-disease pair were tabulated. Exploratory analyses of responses by gender and age were conducted using t and χ2 tests. Results: Of 386 experts invited, 238 (61.7%) responded (mean [SD] age, 52.6 [12.8] years [range 27-93 years]; 126 [52.9%] women and 112 [47.1%] men). Among the experts who responded, 161 (87.9%) agreed that NBSeq for monogenic treatable disorders should be made available to all newborns; 107 (58.5%) agreed that NBSeq should include genes associated with treatable disorders, even if those conditions were low penetrance; 68 (37.2%) agreed that actionable adult-onset conditions should be sequenced in newborns to facilitate cascade testing in parents, and 51 (27.9%) agreed that NBSeq should include screening for conditions with no established therapies or management guidelines. The following 25 genes were recommended by 85% or more of the experts: OTC, G6PC, SLC37A4, CYP11B1, ARSB, F8, F9, SLC2A1, CYP17A1, RB1, IDS, GUSB, DMD, GLUD1, CYP11A1, GALNS, CPS1, PLPBP, ALDH7A1, SLC26A3, SLC25A15, SMPD1, GATM, SLC7A7, and NAGS. Including these, 42 gene-disease pairs were endorsed by at least 80% of experts, and 432 genes were endorsed by at least 50% of experts. Conclusions and Relevance: In this survey study, rare disease experts broadly supported NBSeq for treatable conditions and demonstrated substantial concordance regarding the inclusion of a specific subset of genes in NBSeq.

Combined X-linked familial exudative vitreoretinopathy and retinopathy of prematurity phenotype in an infant with mosaic turner syndrome with ring X chromosome.

BACKGROUND: Retinopathy of prematurity (ROP) and familial exudative vitreoretinopathy (FEVR) are two distinct pathologies of retinal angiogenesis with overlapping clinical features. METHODS: Examination, multimodal imaging, and genetic testing were used to guide diagnosis and treatment. RESULTS: We report a combined phenotype of X-linked FEVR and ROP in a 4-month-old girl with mosaic Turner syndrome with ring X chromosome born at 26 weeks gestational age. She was initially diagnosed with atypical ROP with a vitreous band causing a localized traction retinal detachment, inferotemporal to the macula in the right eye, vessels to posterior zone 2 with no clear ridge temporally in the left eye, and fluorescein leakage in both eyes. Due to the suspicion of concurrent FEVR, genetic testing using a vitreoretinopathy panel was performed which revealed a mosaic Turner syndrome associated with 45,X/46,X,r(X), subsequently confirmed by chromosome analysis. The deleted region in the ring X chromosome included the NDP and RS1 genes. The patient was treated with laser photocoagulation of the peripheral avascular retina and sub-Tenon's triamcinolone injection in both eyes, intravitreal injection of bevacizumab in the left eye, and pars plicata vitrectomy in the right eye. CONCLUSIONS: In premature neonates with atypical ROP, a clinical suspicion of concurrent FEVR or similar vasculopathy is important and genetic testing may elucidate a genetic etiology, which could influence management and prognosis. Turner syndrome can be connected with co-occurring Mendelian gene disorders, particularly in individuals with mosaicism. The concurrence of FEVR and ROP appears to result in atypical and possibly more severe phenotypes.

Disparities in Inherited Retinal Degenerations.

To review disparities in the field of inherited retinal degenerations to establish foundations for future discussions oriented toward finding possible solutions. A narrative overview of the literature. Despite collective efforts towards democratization of genetic testing and investigation, genetic databases containing primarily European populations are heavily relied upon. Access to specialized care and other resources is also still not available to all. Recognizing and addressing disparities and inequities within the field of inherited retinal degenerations will improve our care of these patients and our knowledge of their conditions.

Unique Capabilities of Genome Sequencing for Rare Disease Diagnosis.

Background: Causal variants underlying rare disorders may remain elusive even after expansive gene panels or exome sequencing (ES). Clinicians and researchers may then turn to genome sequencing (GS), though the added value of this technique and its optimal use remain poorly defined. We therefore investigated the advantages of GS within a phenotypically diverse cohort. Methods: GS was performed for 744 individuals with rare disease who were genetically undiagnosed. Analysis included review of single nucleotide, indel, structural, and mitochondrial variants. Results: We successfully solved 218/744 (29.3%) cases using GS, with most solves involving established disease genes (157/218, 72.0%). Of all solved cases, 148 (67.9%) had previously had non-diagnostic ES. We systematically evaluated the 218 causal variants for features requiring GS to identify and 61/218 (28.0%) met these criteria, representing 8.2% of the entire cohort. These included small structural variants (13), copy neutral inversions and complex rearrangements (8), tandem repeat expansions (6), deep intronic variants (15), and coding variants that may be more easily found using GS related to uniformity of coverage (19). Conclusion: We describe the diagnostic yield of GS in a large and diverse cohort, illustrating several types of pathogenic variation eluding ES or other techniques. Our results reveal a higher diagnostic yield of GS, supporting the utility of a genome-first approach, with consideration of GS as a secondary or tertiary test when higher-resolution structural variant analysis is needed or there is a strong clinical suspicion for a condition and prior targeted genetic testing has been negative.

Identification of a novel large multigene deletion and a frameshift indel in PDE6B as the underlying cause of early-onset recessive rod-cone degeneration.

A family, with two affected identical twins with early-onset recessive inherited retinal degeneration, was analyzed to determine the underlying genetic cause of pathology. Exome sequencing revealed a rare and previously reported causative variant (c.1923_1969delinsTCTGGG; p.Asn643Glyfs*29) in the PDE6B gene in the affected twins and their unaffected father. Further investigation, using genome sequencing, identified a novel ∼7.5-kb deletion (Chr 4:670,405-677,862del) encompassing the ATP5ME gene, part of the 5' UTR of MYL5, and a 378-bp (Chr 4:670,405-670,782) region from the 3' UTR of PDE6B in the affected twins and their unaffected mother. Both variants segregated with disease in the family. Analysis of the relative expression of PDE6B, in peripheral blood cells, also revealed a significantly lower level of PDE6B transcript in affected siblings compared to a normal control. PDE6B is associated with recessive rod-cone degeneration and autosomal dominant congenital stationary night blindness. Ophthalmic evaluation of these patients showed night blindness, fundus abnormalities, and peripheral vision loss, which are consistent with PDE6B-associated recessive retinal degeneration. These findings suggest that the loss of PDE6B transcript resulting from the compound heterozygous pathogenic variants is the underlying cause of recessive rod-cone degeneration in the study family.

Novel RCBTB1 variants causing later-onset non-syndromic retinal dystrophy with macular chorioretinal atrophy.

BACKGROUND: Variants in RCBTB1 were recently described to cause a retinal dystrophy with only eight families described to date and a predominant phenotype of macular atrophy and peripheral reticular degeneration. Here, we further evaluate the genotypic and phenotypic characteristics of biallelic RCBTB1-associated retinal dystrophy in a North American clinic population. METHODS: A retrospective analysis of genetic and clinical features was performed in individuals with biallelic variants in RCBTB1. RESULTS: Three unrelated individuals of French-Canadian descent with rare biallelic RCBTB1 variants were identified. All individuals shared a novel p.(Ser342Leu) missense variant; one patient was homozygous whereas the other two each possessed a second unique novel variant p.(Gln120*) and p.(Pro224Leu). All three had macula-predominant disease with symptom onset in the fifth decade of life. CONCLUSION: This report adds to the genetic diversity of RCBTB1-associated disease. These cases confirm the later-onset, relative to many other retinal dystrophies, and macular focus of disease described in most cases to-date. They are thus a reminder of considering hereditary disease in the differential for later-onset macular atrophy.

A hidden structural variation in a known IRD gene: a cautionary tale of two new disease candidate genes.

Rod-cone dystrophy (RCD), also known as retinitis pigmentosa, is an inherited condition leading to vision loss, affecting 1 in 3500 people. More than 270 genes are known to be implicated in the inherited retinal degenerations (IRDs), yet genetic diagnosis for ∼30% of IRD of patients remains elusive despite advances in sequencing technologies. The goal of this study was to determine the genetic causality in a family with RCD. Family members were given a full ophthalmic exam at the Retinal Service at Massachusetts Eye and Ear and consented to genetic testing. Whole-exome sequencing (WES) was performed and variants of interest were Sanger-validated. Functional assays were conducted in zebrafish along with splicing assays in relevant cell lines to determine the impact on retinal function. WES identified variants in two potential candidate genes that segregated with disease: GNL3 (G Protein Nucleolar 3) c.1187 + 3A > C and c.1568-8C > A; and PDE4DIP (Phosphodiester 4D Interacting Protein) c.3868G > A (p.Glu1290Lys) and c.4603G > A (p.Ala1535Thr). Both genes were promising candidates based on their retinal involvement (development and interactions with IRD-associated proteins); however, the functional assays did not validate either gene. Subsequent WES reanalysis with an updated bioinformatics pipeline and widened search parameters led to the detection of a 94-bp duplication in PRPF31 (pre-mRNA Processing Factor 31) c.73_266dup (p.Asp56GlyfsTer33) as the causal variant. Our study demonstrates the importance of thorough functional characterization of new disease candidate genes and the value of reanalyzing next-generation sequencing sequence data, which in our case led to identification of a hidden pathogenic variant in a known IRD gene.

The importance of automation in genetic diagnosis: Lessons from analyzing an inherited retinal degeneration cohort with the Mendelian Analysis Toolkit (MATK).

PURPOSE: In Mendelian disease diagnosis, variant analysis is a repetitive, error-prone, and time consuming process. To address this, we have developed the Mendelian Analysis Toolkit (MATK), a configurable, automated variant ranking program. METHODS: MATK aggregates variant information from multiple annotation sources and uses expert-designed rules with parameterized weights to produce a ranked list of potentially causal solutions. MATK performance was measured by a comparison between MATK-aided and human-domain expert analyses of 1060 families with inherited retinal degeneration (IRD), analyzed using an IRD-specific gene panel (589 individuals) and exome sequencing (471 families). RESULTS: When comparing MATK-assisted analysis with expert curation in both the IRD-specific gene panel and exome sequencing (1060 subjects), 97.3% of potential solutions found by experts were also identified by the MATK-assisted analysis (541 solutions identified with MATK of 556 solutions found by conventional analysis). Furthermore, MATK-assisted analysis identified 114 additional potential solutions from the 504 cases unsolved by conventional analysis. CONCLUSION: MATK expedites the process of identification of likely solving variants in Mendelian traits, and reduces variability stemming from human error and researcher bias. MATK facilitates data reanalysis to keep up with the constantly improving annotation sources and next-generation sequencing processing pipelines. The software is open source and available at https://gitlab.com/matthew_maher/mendelanalysis.