Clinician-Driven Reanalysis of Exome Sequencing Data From Patients With Inherited Retinal Diseases.

Importance: Despite advances in next-generation sequencing (NGS), a significant proportion of patients with inherited retinal disease (IRD) remain undiagnosed after initial genetic testing. Exome sequencing (ES) reanalysis in the clinical setting has been suggested as one method for improving diagnosis of IRD. Objective: To investigate the association of clinician-led reanalysis of ES data, which incorporates updated clinical information and comprehensive bioinformatic analysis, with the diagnostic yield in a cohort of patients with IRDs in Korea. Design, Setting, and Participants: This was a multicenter prospective cohort study involving 264 unrelated patients with IRDs, conducted in Korea between March 2018 and February 2020. Comprehensive ophthalmologic examinations and ES analyses were performed, and ES data were reanalyzed by an IRD specialist for single nucleotide variants, copy number variants, mobile element insertions, and mitochondrial variants. Data were analyzed from March to July 2023. Main Outcomes and Measures: Diagnostic rate of conventional bioinformatic analysis and clinician-driven ES reanalysis. Results: A total of 264 participants (151 [57.2%] male; mean [SD] age at genetic testing, 33.6 [18.9] years) were enrolled, including 129 patients (48.9%) with retinitis pigmentosa and 26 patients (9.8%) with Stargardt disease or macular dystrophy. Initial bioinformatic analysis diagnosed 166 patients (62.9%). Clinician-driven reanalysis identified the molecular cause of diseases in an additional 22 patients, corresponding to an 8.3-percentage point increase in diagnostic rate. Key factors associated with new molecular diagnoses included clinical phenotype updates (4 patients) and detection of previously overlooked variation, such as structural variants (9 patients), mitochondrial variants (3 patients), filtered or not captured variants (4 patients), and noncanonical splicing variants (2 patients). Among the 22 patients, variants in 7 patients (31.8%) were observed in the initial analysis but not reported to patients, while those in the remaining 15 patients (68.2%) were newly detected by the ES reanalysis. Conclusions and Relevance: In this cohort study, clinician-centered reanalysis of ES data was associated with improved molecular diagnostic yields in patients with IRD. This approach is important for uncovering missed genetic causes of retinal disease.

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.

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.

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.

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 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.

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.

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.

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.

Broadening INPP5E phenotypic spectrum: detection of rare variants in syndromic and non-syndromic IRD.

Pathogenic variants in INPP5E cause Joubert syndrome (JBTS), a ciliopathy with retinal involvement. However, despite sporadic cases in large cohort sequencing studies, a clear association with non-syndromic inherited retinal degenerations (IRDs) has not been made. We validate this association by reporting 16 non-syndromic IRD patients from ten families with bi-allelic mutations in INPP5E. Additional two patients showed early onset IRD with limited JBTS features. Detailed phenotypic description for all probands is presented. We report 14 rare INPP5E variants, 12 of which have not been reported in previous studies. We present tertiary protein modeling and analyze all INPP5E variants for deleteriousness and phenotypic correlation. We observe that the combined impact of INPP5E variants in JBTS and non-syndromic IRD patients does not reveal a clear genotype-phenotype correlation, suggesting the involvement of genetic modifiers. Our study cements the wide phenotypic spectrum of INPP5E disease, adding proof that sequence defects in this gene can lead to early-onset non-syndromic IRD.

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.