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.

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.

Two novel non-coding single nucleotide variants in the DNase1 hypersensitivity site of PRDM13 causing North Carolina macular dystrophy in Korea

Purpose

Pathogenic variants in North Carolina macular dystrophy (NCMD) have rarely been reported in the East Asian population. Herein, we reported novel variants of NCMD in 2 Korean families.

Methods

The regions associated with NCMD were analyzed with genome sequencing, and variants were filtered based on the minor allele frequency (0.5%) and heterozygosity. Non-coding variants were functionally annotated using multiple computational tools.

Results

We identified two rare novel variants, chr6:g.99,598,914T>C (hg38; V17) and chr6:g.99,598,926G>A (hg38; V18) upstream of PRDM13 in families A and B, respectively. In Family 1, Grade 2 NCMD and a best-corrected visual acuity of 20/25 and 20/200 in the right and left eyes, respectively, were observed. In Family B, all affected individuals had Grade 1 NCMD with characteristic confluent drusen at the fovea and a best-corrected visual acuity of 20/20 in both eyes. These two variants are 10-22 bp downstream of the reported V10 variant within the DNase1 hypersensitivity site. This site is associated with progressive bifocal chorioretinal atrophy and congenital posterior polar chorioretinal hypertrophy and lies in the putative enhancer site of PRDM13.

Conclusion

We identified two novel NCMD variants in the Korean population and further validated the regulatory role of the DNase1 hypersensitivity site upstream of PRDM13.

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.

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.

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.

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.

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.

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.

Non-syndromic Retinal Degeneration Caused by Pathogenic Variants in Joubert Syndrome Genes

Inherited retinal degenerations (IRDs) are a group of genetic disorders characterized by progressive dysfunction and loss of photoreceptors. IRDs are classified as non-syndromic or syndromic, depending on whether retinal degeneration manifests alone or in combination with other associated symptoms. Joubert syndrome (JBTS) is a genetically and clinically heterogeneous disorder affecting the central nervous system and other organs and tissues, including the neuroretina. To date, 39 genes have been associated with JBTS, a majority of which encode structural or functional components of the primary cilium, a specialized sensory organelle present in most post-mitotic cells, including photoreceptors. The use of whole exome and IRD panel next-generation sequencing in routine diagnostics of non-syndromic IRD cases led to the discovery of pathogenic variants in JBTS genes that cause photoreceptor loss without other syndromic features. Here, we recapitulate these findings, describing the JBTS gene defects leading to non-syndromic IRDs.

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.

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.

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/3500 people. Over 270 genes are known to be implicated in the inherited retinal degenerations (IRDs), yet genetic diagnosis for ~30% 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 Rod-cone dystrophy (RCD). Family members were given a full ophthalmic exam at the Retinal Service at MEE 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 94bp 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 NGS sequence data, which in our case led to identification of a hidden pathogenic variant in a known IRD gene.

WDR34, a candidate gene for non-syndromic rod-cone dystrophy

Rod-cone dystrophy (RCD), also called retinitis pigmentosa, is characterized by rod followed by cone photoreceptor degeneration, leading to gradual visual loss. Mutations in over 65 genes have been associated with non-syndromic RCD explaining 60% to 70% of cases, with novel gene defects possibly accounting for the unsolved cases. Homozygosity mapping and whole-exome sequencing applied to a case of autosomal recessive non-syndromic RCD from a consanguineous union identified a homozygous variant in WDR34. Mutations in WDR34 have been previously associated with severe ciliopathy syndromes possibly associated with a retinal dystrophy. This is the first report of a homozygous mutation in WDR34 associated with non-syndromic RCD.

Mutations in the Kinesin-2 Motor KIF3B Cause an Autosomal-Dominant Ciliopathy

Kinesin-2 enables ciliary assembly and maintenance as an anterograde intraflagellar transport (IFT) motor. Molecular motor activity is driven by a heterotrimeric complex comprised of KIF3A and KIF3B or KIF3C plus one non-motor subunit, KIFAP3. Using exome sequencing, we identified heterozygous KIF3B variants in two unrelated families with hallmark ciliopathy phenotypes. In the first family, the proband presents with hepatic fibrosis, retinitis pigmentosa, and postaxial polydactyly; he harbors a de novo c.748G>C (p.Glu250Gln) variant affecting the kinesin motor domain encoded by KIF3B. The second family is a six-generation pedigree affected predominantly by retinitis pigmentosa. Affected individuals carry a heterozygous c.1568T>C (p.Leu523Pro) KIF3B variant segregating in an autosomal-dominant pattern. We observed a significant increase in primary cilia length in vitro in the context of either of the two mutations while variant KIF3B proteins retained stability indistinguishable from wild type. Furthermore, we tested the effects of KIF3B mutant mRNA expression in the developing zebrafish retina. In the presence of either missense variant, rhodopsin was sequestered to the photoreceptor rod inner segment layer with a concomitant increase in photoreceptor cilia length. Notably, impaired rhodopsin trafficking is also characteristic of recessive KIF3B models as exemplified by an early-onset, autosomal-recessive, progressive retinal degeneration in Bengal cats; we identified a c.1000G>A (p.Ala334Thr) KIF3B variant by genome-wide association study and whole-genome sequencing. Together, our genetic, cell-based, and in vivo modeling data delineate an autosomal-dominant syndromic retinal ciliopathy in humans and suggest that multiple KIF3B pathomechanisms can impair kinesin-driven ciliary transport in the photoreceptor.

A combined RNA-seq and whole genome sequencing approach for identification of non-coding pathogenic variants in single families

Inherited retinal degenerations (IRDs) are at the focus of current genetic therapeutic advancements. For a genetic treatment such as gene therapy to be successful, an accurate genetic diagnostic is required. Genetic diagnostics relies on the assessment of the probability that a given DNA variant is pathogenic. Non-coding variants present a unique challenge for such assessments as compared to coding variants. For one, non-coding variants are present at much higher number in the genome than coding variants. In addition, our understanding of the rules that govern the non-coding regions of the genome is less complete than our understanding of the coding regions. Methods that allow for both the identification of candidate non-coding pathogenic variants and their functional validation may help overcome these caveats allowing for a greater number of patients to benefit from advancements in genetic therapeutics. We present here an unbiased approach combining whole genome sequencing (WGS) with patient-induced pluripotent stem cell (iPSC)-derived retinal organoids (ROs) transcriptome analysis. With this approach, we identified and functionally validated a novel pathogenic non-coding variant in a small family with a previously unresolved genetic diagnosis.

Longitudinal Clinical Follow-up and Genetic Spectrum of Patients With Rod-Cone Dystrophy Associated With Mutations in PDE6A and PDE6B

Importance

A precise phenotypic characterization of retinal dystrophies is needed for disease modeling as a basis for future therapeutic interventions.

Objective

To compare genotype, phenotype, and structural changes in patients with rod-cone dystrophy (RCD) associated with mutations in PDE6A or PDE6B.

Design, Setting, and Participants

In a retrospective cohort study conducted in Paris, France, from January 2007 to September 2017, 54 patients from a cohort of 1095 index patients with RCD underwent clinical examination, including personal and familial history, best-corrected visual acuity (BCVA), color vision, slitlamp examination, full-field electroretinography, kinetic visual fields (VFs), retinophotography, optical coherence tomography, near-infrared fundus autofluorescence, and short-wavelength fundus autofluorescence imaging. Genotyping was performed using microarray analysis, targeted next-generation sequencing, and Sanger sequencing validation with familial segregation when possible. Data were analyzed from September 1, 2017, to February 1, 2018. Clinical variables were subsequently analyzed in 2018.

Main Outcomes and Measures

Phenotype and genotype comparison of patients carrying mutations in PDE6A or PDE6B.

Results

Of the 54 patients included in the study, 19 patients of 17 families (11 women [58%]; mean [SD] age at diagnosis, 14.83 [10.63] years) carried pathogenic mutations in PDE6A, and 35 patients of 26 families (17 women [49%]; mean [SD] age at diagnosis, 21.10 [11.56] years) had mutations in PDE6B, accounting for prevalences of 1.6% and 2.4%, respectively. Among 49 identified genetic variants, 14 in PDE6A and 15 in PDE6B were novel. Overall, phenotypic analysis revealed no substantial differences between the 2 groups except for night blindness as a presenting symptom that was noted to be more prevalent in the PDE6A than PDE6B group (80% vs 37%, respectively; P = .005). The mean binocular BCVA and VF decrease over time (measured as mean individual slopes coefficients) was comparable between patients with PDE6A and PDE6B mutations: 0.04 (0.12) vs 0.02 (0.05) for BCVA (P = .89) and 14.33 (7.12) vs 13.27 (6.77) for VF (P = .48).

Conclusions and Relevance

Mutations in PDE6A and PDE6B accounted for 1.6% and 2.4%, respectively, in a cohort of French patients with RCD. The functional and structural findings reported may constitute the basis of disease modeling that might be used for better prognostic estimation and candidate selection for photoreceptor therapeutic rescue.

Expanding the phenotypic spectrum in RDH12-associated retinal disease

Retinol dehydrogenase 12, RDH12, plays a pivotal role in the visual cycle to ensure the maintenance of normal vision. Alterations in activity of this protein result in photoreceptor death and decreased vision beginning at an early age and progressing to substantial vision loss later in life. Here we describe 11 patients with retinal degeneration that underwent next-generation sequencing (NGS) with a targeted panel of all currently known inherited retinal degeneration (IRD) genes and whole-exome sequencing to identify the genetic causality of their retinal disease. These patients display a range of phenotypic severity prompting clinical diagnoses of macular dystrophy, cone-rod dystrophy, retinitis pigmentosa, and early-onset severe retinal dystrophy all attributed to biallelic recessive mutations in RDH12. We report 15 causal alleles and expand the repertoire of known RDH12 mutations with four novel variants: c.215A > G (p.Asp72Gly); c.362T > C (p.Ile121Thr); c.440A > C (p.Asn147Thr); and c.697G > A (p.Val233Ille). The broad phenotypic spectrum observed with biallelic RDH12 mutations has been observed in other genetic forms of IRDs, but the diversity is particularly notable here given the prior association of RDH12 primarily with severe early-onset disease. This breadth emphasizes the importance of broad genetic testing for inherited retinal disorders and extends the pool of individuals who may benefit from imminent gene-targeted therapies.

Clinical Features of a Retinopathy Associated With a Dominant Allele of the RGR Gene

Purpose

We describe the clinical features in two pedigrees with dominantly inherited retinopathy segregating the previously reported frameshifting mutation, c.836dupG (p.Ile280Asn*78) in the terminal exon of the RGR gene, and compare their haplotypes to that of the previously reported pedigree.

Methods

The probands were ascertained at West Virginia University Eye Institute (WVU) and Moorfields Eye Hospital (MEH) through next generation sequencing (NGS) and whole genome sequencing (WGS) respectively. Clinical data included visual acuity (VA), visual fields, fundus autofluorescence (FAF), optical coherence tomography (OCT), and electroretinography (ERG). Haplotype analysis was performed using Sanger sequencing of the DNA from the molecularly ascertained individuals from the three pedigrees.

Results

Nine heterozygous mutation carriers were identified in two families. Four carriers were asymptomatic; five carriers had variable VA reduction, visual field constriction, and experienced difficulty under dim illumination. Fundus examination of the asymptomatic carriers showed diffuse or reticular pigmentation of the retina; the symptomatic carriers had chorioretinal atrophy. FAF imaging showed widespread signal loss in advanced retinopathy, and reticular hyperautofluorescence in mild cases. OCT showed loss of outer retinal lamina in advanced disease. ERG showed moderate-to-severe rod–cone dysfunction in two symptomatic carriers; and was normal in three asymptomatic carriers. A shared haplotype flanking the mutation of up to 6.67 Mb was identified in both families. Within this region, 1.27 Mb were shared with the first family reported with this retinopathy.

Conclusions

The clinical data suggest a variable and slow degeneration of the RPE. A shared chromosomal segment surrounding the RGR gene suggests a single ancestral mutational event underlying all three families.

Changes in extracellular matrix cause RPE cells to make basal deposits and activate the alternative complement pathway

The design of efficient therapies for age-related macular degeneration (AMD) is limited by our understanding of the pathogenesis of basal deposits, which form between retinal pigment epithelium (RPE) and Bruch's membrane (BrM) early in disease, and involve activation of the complement system. To investigate the roles of BrM, RPE and complement in an AMD, we generated abnormal extracellular matrix (ECM) using CRISPR-edited ARPE-19 cells. We introduced to these cells the p.R345W mutation in EFEMP1, which causes early-onset macular degeneration. The abnormal ECM binds active complement C3 and causes the formation of basal deposits by normal human fetal (hf)RPE cells. Human fetal RPE (hfRPE) cells grown on abnormal ECM or BrM explants from AMD donors show chronic activation of the alternative complement pathway by excessive deposition of C3b. This process is exacerbated by impaired ECM turnover via increased matrix metalloproteinase-2 activity. The local cleavage of C3 via convertase-independent mechanisms can be a new therapeutic target for early AMD.

Contribution of noncoding pathogenic variants to RPGRIP1-mediated inherited retinal degeneration

Purpose:

With the advent of gene therapies for inherited retinal degenerations (IRDs), genetic diagnostics will have an increasing role in clinical decision-making. Yet the genetic cause of disease cannot be identified using exon-based sequencing for a significant portion of patients. We hypothesized that non-coding mutations contribute significantly to the genetic causality of IRDs and evaluated patients with single coding mutations in RPGRIP1 to test this hypothesis.

Methods:

IRD families underwent targeted panel sequencing. Unsolved cases were explored by whole exome and genome sequencing looking for additional mutations. Candidate mutations were then validated by Sanger sequencing, quantitative PCR, and in vitro splicing assays in two cell lines analyzed through amplicon sequencing.

Results:

Among 1722 families, three had biallelic loss of function mutations in RPGRIP1 while seven had a single disruptive coding mutation. Whole exome and genome sequencing revealed potential non-coding mutations in these seven families. In six, the non-coding mutations were shown to lead to loss of function in vitro.

Conclusion:

Non-coding mutations were identified in 6 of 7 families with single coding mutations in RPGRIP1. The results suggest that non-coding mutations contribute significantly to the genetic causality of IRDs and RPGRIP1–mediated IRDs are more common than previously thought.

Ift172 conditional knock-out mice exhibit rapid retinal degeneration and protein trafficking defects

Intraflagellar transport (IFT) is a bidirectional transport process that occurs along primary cilia and specialized sensory cilia, such as photoreceptor outersegments. Genes coding for various IFT components are associated with ciliopathies. Mutations in IFT172 lead to diseases ranging from isolated retinal degeneration to severe syndromic ciliopathies. In this study, we created a mouse model of IFT172-associated retinal degeneration to investigate the ocular disease mechanism. We found that depletion of IFT172 in rod photoreceptors leads to a rapid degeneration of the retina, with severely reduced electroretinography (ERG) responses by 1 month and complete outer-nuclear layer (ONL) degeneration by 2 months. We investigated molecular mechanisms of degeneration and show that IFT172 protein reduction leads to mislocalization of specific photoreceptor outersegment (OS) proteins (RHO, RP1, IFT139), aberrant light-driven translocation of alpha transducin and altered localization of glioma-associated oncogene family member 1 (GLI1). This mouse model exhibits key features of the retinal phenotype observed in patients with IFT172-associated blindness and can be used for in vivo testing of ciliopathy therapies.

The importance of genetic testing as demonstrated by two cases of CACNA1F-associated retinal generation misdiagnosed as LCA

PURPOSE

To describe in detail cases with an initial diagnosis of Leber congenital amaurosis that were later found to have a hemizygous mutation in the CACNA1F gene.

METHODS

The patients underwent a detailed ophthalmological evaluation and full-field electroretinography (ERG). Selective targeted capture and whole-exome next-generation sequencing (NGS) were used to find the disease-causing mutations.

RESULTS

Patient 1 presented at age 3 months with nystagmus, normal visual attention, and a normal fundus exam. ERG responses were severely decreased. Patient 2 presented with nystagmus, severe hyperopia, esotropia, and visual acuity of 20/360 oculus dexter (OD) and 20/270 oculus sinister (OS) at age 5 months. His fundus exam showed slightly increased pigmentation around the foveae. The scotopic ERG responses were severely decreased and photopic responses mildly decreased. Based on the initial presentation, both patients received the clinical diagnosis of Leber congenital amaurosis (LCA). However, genetic testing showed no mutations in known LCA genes. Instead, broader genetic testing using NGS showed point mutations in the CACNA1F gene, which is reported to be associated with type 2 congenital stationary night blindness (CSNB2).

CONCLUSIONS

These two cases demonstrate the clinical overlap between LCA and CSNB in infants and young children. Genetic testing is an essential tool in these cases and provides a more accurate diagnosis and prognosis for patients with inherited retinal degenerative disorders.

Photoreceptor Cilia and Retinal Ciliopathies

Photoreceptors are sensory neurons designed to convert light stimuli into neurological responses. This process, called phototransduction, takes place in the outer segments (OS) of rod and cone photoreceptors. OS are specialized sensory cilia, with analogous structures to those present in other nonmotile cilia. Deficient morphogenesis and/or dysfunction of photoreceptor sensory cilia (PSC) caused by mutations in a variety of photoreceptor-specific and common cilia genes can lead to inherited retinal degenerations (IRDs). IRDs can manifest as isolated retinal diseases or syndromic diseases. In this review, we describe the structure and composition of PSC and different forms of ciliopathies with retinal involvement. We review the genetics of the IRDs, which are monogenic disorders but genetically diverse with regard to causality.

Copy-number variation is an important contributor to the genetic causality of inherited retinal degenerations

PURPOSE

Despite substantial progress in sequencing, current strategies can genetically solve only approximately 55-60% of inherited retinal degeneration (IRD) cases. This can be partially attributed to elusive mutations in the known IRD genes, which are not easily identified by the targeted next-generation sequencing (NGS) or Sanger sequencing approaches. We hypothesized that copy-number variations (CNVs) are a major contributor to the elusive genetic causality of IRDs.

METHODS

Twenty-eight cases previously unsolved with a targeted NGS were investigated with whole-genome single-nucleotide polymorphism (SNP) and comparative genomic hybridization (CGH) arrays.

RESULTS

Deletions in the IRD genes were detected in 5 of 28 families, including a de novo deletion. We suggest that the de novo deletion occurred through nonallelic homologous recombination (NAHR) and we constructed a genomic map of NAHR-prone regions with overlapping IRD genes. In this article, we also report an unusual case of recessive retinitis pigmentosa due to compound heterozygous mutations in SNRNP200, a gene that is typically associated with the dominant form of this disease.

CONCLUSIONS

CNV mapping substantially increased the genetic diagnostic rate of IRDs, detecting genetic causality in 18% of previously unsolved cases. Extending the search to other structural variations will probably demonstrate an even higher contribution to genetic causality of IRDs.Genet Med advance online publication 13 October 2016.

Efficient In Silico Identification of a Common Insertion in the MAK Gene which Causes Retinitis Pigmentosa

Background

Next generation sequencing (NGS) offers a rapid and comprehensive method of screening for mutations associated with retinitis pigmentosa and related disorders. However, certain sequence alterations such as large insertions or deletions may remain undetected using standard NGS pipelines. One such mutation is a recently-identified Alu insertion into the Male Germ Cell-Associated Kinase (MAK) gene, which is missed by standard NGS-based variant callers. Here, we developed an in silico method of searching NGS raw sequence reads to detect this mutation, without the need to recalculate sequence alignments or to screen every sample by PCR.

Methods

The Linux program grep was used to search for a 23 bp “probe” sequence containing the known junction sequence of the insert. A corresponding search was performed with the wildtype sequence. The matching reads were counted and further compared to the known sequences of the full wildtype and mutant genomic loci. (See https://github.com/MEEIBioinformaticsCenter/grepsearch.)

Results

In a test sample set consisting of eleven previously published homozygous mutants, detection of the MAK-Alu insertion was validated with 100% sensitivity and specificity. As a discovery cohort, raw NGS reads from 1,847 samples (including custom and whole exome selective capture) were searched in ~1 hour on a local computer cluster, yielding an additional five samples with MAK-Alu insertions and solving two previously unsolved pedigrees. Of these, one patient was homozygous for the insertion, one compound heterozygous with a missense change on the other allele (c. 46G>A; p.Gly16Arg), and three were heterozygous carriers.

Conclusions

Using the MAK-Alu grep program proved to be a rapid and effective method of finding a known, disease-causing Alu insertion in a large cohort of patients with NGS data. This simple approach avoids wet-lab assays or computationally expensive algorithms, and could also be used for other known disease-causing insertions and deletions.

Targeted exon sequencing in Usher syndrome type I

PURPOSE

Patients with Usher syndrome type I (USH1) have retinitis pigmentosa, profound congenital hearing loss, and vestibular ataxia. This syndrome is currently thought to be associated with at least six genes, which are encoded by over 180 exons. Here, we present the use of state-of-the-art techniques in the molecular diagnosis of a cohort of 47 USH1 probands.

METHODS

The cohort was studied with selective exon capture and next-generation sequencing of currently known inherited retinal degeneration genes, comparative genomic hybridization, and Sanger sequencing of new USH1 exons identified by human retinal transcriptome analysis.

RESULTS

With this approach, we were able to genetically solve 14 of the 47 probands by confirming the biallelic inheritance of mutations. We detected two likely pathogenic variants in an additional 19 patients, for whom family members were not available for cosegregation analysis to confirm biallelic inheritance. Ten patients, in addition to primary disease-causing mutations, carried rare likely pathogenic USH1 alleles or variants in other genes associated with deaf-blindness, which may influence disease phenotype. Twenty-one of the identified mutations were novel among the 33 definite or likely solved patients. Here, we also present a clinical description of the studied cohort at their initial visits.

CONCLUSIONS

We found a remarkable genetic heterogeneity in the studied USH1 cohort with multiplicity of mutations, of which many were novel. No obvious influence of genotype on phenotype was found, possibly due to small sample sizes of the genotypes under study.

Mutations in IFT172 cause isolated retinal degeneration and Bardet-Biedl syndrome

Primary cilia are sensory organelles present on most mammalian cells. The assembly and maintenance of primary cilia are facilitated by intraflagellar transport (IFT), a bidirectional protein trafficking along the cilium. Mutations in genes coding for IFT components have been associated with a group of diseases called ciliopathies. These genetic disorders can affect a variety of organs including the retina. Using whole exome sequencing in three families, we identified mutations in Intraflagellar Transport 172 Homolog [IFT172 (Chlamydomonas)] that underlie an isolated retinal degeneration and Bardet-Biedl syndrome. Extensive functional analyses of the identified mutations in cell culture, rat retina and in zebrafish demonstrated their hypomorphic or null nature. It has recently been reported that mutations in IFT172 cause a severe ciliopathy syndrome involving skeletal, renal, hepatic and retinal abnormalities (Jeune and Mainzer-Saldino syndromes). Here, we report for the first time that mutations in this gene can also lead to an isolated form of retinal degeneration. The functional data for the mutations can partially explain milder phenotypes; however, the involvement of modifying alleles in the IFT172-associated phenotypes cannot be excluded. These findings expand the spectrum of disease associated with mutations in IFT172 and suggest that mutations in genes originally reported to be associated with syndromic ciliopathies should also be considered in subjects with non-syndromic retinal dystrophy.

Further insights into GPR179: expression, localization, and associated pathogenic mechanisms leading to complete congenital stationary night blindness

PURPOSE

Mutations in GPR179, which encodes the G protein-coupled receptor 179, lead to autosomal recessive complete (c) congenital stationary night blindness (CSNB), which is characterized by an ON-bipolar retinal cell dysfunction. This study further defined the exact site of Gpr179 expression and its protein localization in human retina and elucidated the pathogenic mechanism of the reported missense and splice site mutations.

METHODS

RNA in situ hybridization was performed with mouse retinal sections. A commercially available antibody was validated with GPR179-overexpressing COS-1 cells and applied to human retinal sections. Live-cell extracellular staining along with subsequent intracellular immunolocalization and ELISA studies were performed using mammalian cells overexpressing wild-type or missense mutated GPR179. Wild-type and splice site-mutated mini-gene constructs were transiently transfected, and RNA was extracted. RT-PCR-amplified products were cloned, and Sanger sequenced.

RESULTS

Mouse Gpr179 transcript was expressed in the upper part of the inner nuclear layer, and the respective human protein localized at the dendritic tips of bipolar cells in human retina. The missense mutations p.Tyr220Cys, p.Gly455Asp, and p.His603Tyr led to severely reduced cell surface localization, whereas p.Asp126His did not. The mutated splice donor site altered GPR179 splicing.

CONCLUSIONS

Our findings indicate that the site of expression and protein localization of human and mouse GPR179 is similar to that of other proteins implicated in cCSNB. For most of the mutations identified so far, loss of the GPR179 protein function seems to be the underlying pathogenic mechanism leading to this form of cCSNB.

A Study into the Evolutionary Divergence of the Core Promoter Elements of PRPF31 and TFPT

Mutations in PRPF31 have been implicated in retinitis pigmentosa, a blinding disease caused by degeneration of rod photoreceptors. The disease mechanism in the majority of cases is haploinsufficiency. Crucially, attempts at generation of animal models of disease have proved unsuccessful, yielding animals with a visual phenotype that does not mirror human disease. This suggests that, in these animals, the transcriptional regulation of PRPF31 is different to humans and compared to other species. Study of the evolution of the PRPF31 core promoter has important implications for our understanding of human disease, as disease phenotype is modified by differentially expressed alleles in the population.

PRPF31 lies in a head-to-head arrangement with TFPT, a gene involved in cellular apoptosis. The two genes were shown to share common regulatory elements in the human genome. In this study, the core promoters of PRPF31 and TFPT were characterised by dual-luciferase reporter assay using genomic DNA from the green monkey, domestic dog and house mouse. It was found that the core promoters were conserved between human and monkey.

In dog, the TFPT core promoter was conserved, but different PRPF31 gene architecture meant the gene was controlled by a long-range promoter lying some 2000bp from the transcription start site.

There was very low level of conservation (<20%) of the PRPF31 5′ region between mouse and human. It was shown that mouse populations did not show variable Prpf31 expression levels, revealing a potential explanation for the lack of phenotype observed in the Prpf31 knock-out mouse model.

Development and application of a next-generation-sequencing (NGS) approach to detect known and novel gene defects underlying retinal diseases

Background

Inherited retinal disorders are clinically and genetically heterogeneous with more than 150 gene defects accounting for the diversity of disease phenotypes. So far, mutation detection was mainly performed by APEX technology and direct Sanger sequencing of known genes. However, these methods are time consuming, expensive and unable to provide a result if the patient carries a new gene mutation. In addition, multiplicity of phenotypes associated with the same gene defect may be overlooked.

Methods

To overcome these challenges, we designed an exon sequencing array to target 254 known and candidate genes using Agilent capture. Subsequently, 20 DNA samples from 17 different families, including four patients with known mutations were sequenced using Illumina Genome Analyzer IIx next-generation-sequencing (NGS) platform. Different filtering approaches were applied to identify the genetic defect. The most likely disease causing variants were analyzed by Sanger sequencing. Co-segregation and sequencing analysis of control samples validated the pathogenicity of the observed variants.

Results

The phenotype of the patients included retinitis pigmentosa, congenital stationary night blindness, Best disease, early-onset cone dystrophy and Stargardt disease. In three of four control samples with known genotypes NGS detected the expected mutations. Three known and five novel mutations were identified in NR2E3, PRPF3, EYS, PRPF8, CRB1, TRPM1 and CACNA1F. One of the control samples with a known genotype belongs to a family with two clinical phenotypes (Best and CSNB), where a novel mutation was identified for CSNB. In six families the disease associated mutations were not found, indicating that novel gene defects remain to be identified.

Conclusions

In summary, this unbiased and time-efficient NGS approach allowed mutation detection in 75% of control cases and in 57% of test cases. Furthermore, it has the possibility of associating known gene defects with novel phenotypes and mode of inheritance.