Long-Read Nanopore Sequencing of RPGR ORF15 is Enhanced Following DNase I Treatment of MinION Flow Cells.

INTRODUCTION: RPGR ORF15 is an exon present almost exclusively in the retinal transcript of RPGR. It is purine-rich, repetitive and notoriously hard to sequence, but is a hotspot for mutations causing X-linked retinitis pigmentosa. METHODS: Long-read nanopore sequencing on MinION and Flongle flow cells was used to sequence RPGR ORF15 in genomic DNA from patients with inherited retinal dystrophy. A flow cell wash kit was used on a MinION flow cell to increase yield. Findings were confirmed by PacBio SMRT long-read sequencing. RESULTS: We showed that long-read nanopore sequencing successfully reads through a 2 kb PCR-amplified fragment containing ORF15. We generated reads of sufficient quality and cumulative read-depth to detect pathogenic RP-causing variants. However, we observed that this G-rich, repetitive DNA segment rapidly blocks the available pores, resulting in sequence yields less than 5% of the expected output. This limited the extent to which samples could be pooled, increasing cost. We tested the utility of a MinION wash kit containing DNase I to digest DNA fragments remaining on the flow cell, regenerating the pores. Use of the DNase I treatment allowed repeated re-loading, increasing the sequence reads obtained. Our customised workflow was used to screen pooled amplification products from previously unsolved inherited retinal disease (IRD) in patients, identifying two new cases with pathogenic ORF15 variants. DISCUSSION: We report the novel finding that long-read nanopore sequencing can read through RPGR-ORF15, a DNA sequence not captured by short-read next-generation sequencing (NGS), but with a more reduced yield. Use of a flow cell wash kit containing DNase I unblocks the pores, allowing reloading of further library aliquots over a 72-h period, increasing yield. The workflow we describe provides a novel solution to the need for a rapid, robust, scalable, cost-effective ORF15 screening protocol.

Haplotyping Using Long-Range PCR and Nanopore Sequencing to Phase Variants: Lessons Learned From the ABCA4 Locus.

Short-read next-generation sequencing has revolutionized our ability to identify variants underlying inherited diseases; however, it does not allow the phasing of variants to clarify their diagnostic interpretation. The advent of widespread, increasingly accurate long-read sequencing has opened up new applications not currently available through short-read next-generation sequencing. One such use is the ability to phase variants to clarify their diagnostic interpretation and to investigate the increasingly prevalent role of cis-acting variants in the pathogenesis of the inherited disease, so-called complex alleles. Complex alleles are becoming an increasingly prevalent part of the study of genes associated with inherited diseases, for example, in ABCA4-related diseases. We sought to establish a cost-effective method to phase contiguous segments of the 130-kb ABCA4 locus by long-read sequencing of overlapping amplification products. Using the comprehensively characterized CEPH sample, NA12878, we verified the accuracy and robustness of our assay. However, in-field assessment of its utility using clinical test cases was hampered by the paucity and distribution of identified variants and by PCR chimerism, particularly where the number of PCR cycles was high. Despite this, we were able to construct robust phase blocks of up to 94.9 kb, representing 73% of the ABCA4 locus. We conclude that, although haplotype analysis of variants located within discrete amplification products was robust and informative, the stitching together of larger phase blocks using overlapping single-molecule reads remained practically challenging.

Targeted nanopore sequencing enables complete characterisation of structural deletions initially identified using exon-based short-read sequencing strategies.

BACKGROUND: The widespread adoption of exome sequencing has greatly increased the rate of genetic diagnosis for inherited conditions. However, the detection and validation of large deletions remains challenging. While numerous bioinformatics approaches have been developed to detect deletions from whole - exome sequencing and targeted panels, further work is typically required to define the physical breakpoints or integration sites. Accurate characterisation requires either expensive follow - up whole - genome sequencing or the time - consuming, laborious process of PCR walking, both of which are challenging when dealing with the repeat sequences which frequently intersect deletion breakpoints. The aim of this study was to develop a cost-effective, long-range sequencing method to characterise deletions. METHODS: Genomic DNA was amplified with primers spanning the deletion using long-range PCR and the products purified. Sequencing was performed on MinION flongle flowcells. The resulting fast5 files were basecalled using Guppy, trimmed using Porechop and aligned using Minimap2. Filtering was performed using NanoFilt. Nanopore sequencing results were verified by Sanger sequencing. RESULTS: Four cases with deletions detected following comparative read-depth analysis of targeted short-read sequencing were analysed. Nanopore sequencing defined breakpoints at the molecular level in all cases including homozygous breakpoints in EYS, CNGA1 and CNGB1 and a heterozygous deletion in PRPF31. All breakpoints were verified by Sanger sequencing. CONCLUSIONS: In this study, a quick, accurate and cost - effective method is described to characterise deletions identified from exome, and similar data, using nanopore sequencing.

Cost-effective sequence analysis of 113 genes in 1,192 probands with retinitis pigmentosa and Leber congenital amaurosis.

Introduction: Retinitis pigmentosa (RP) and Leber congenital amaurosis (LCA) are two groups of inherited retinal diseases (IRDs) where the rod photoreceptors degenerate followed by the cone photoreceptors of the retina. A genetic diagnosis for IRDs is challenging since >280 genes are associated with these conditions. While whole exome sequencing (WES) is commonly used by diagnostic facilities, the costs and required infrastructure prevent its global applicability. Previous studies have shown the cost-effectiveness of sequence analysis using single molecule Molecular Inversion Probes (smMIPs) in a cohort of patients diagnosed with Stargardt disease and other maculopathies. Methods: Here, we introduce a smMIPs panel that targets the exons and splice sites of all currently known genes associated with RP and LCA, the entire RPE65 gene, known causative deep-intronic variants leading to pseudo-exons, and part of the RP17 region associated with autosomal dominant RP, by using a total of 16,812 smMIPs. The RP-LCA smMIPs panel was used to screen 1,192 probands from an international cohort of predominantly RP and LCA cases. Results and discussion: After genetic analysis, a diagnostic yield of 56% was obtained which is on par with results from WES analysis. The effectiveness and the reduced costs compared to WES renders the RP-LCA smMIPs panel a competitive approach to provide IRD patients with a genetic diagnosis, especially in countries with restricted access to genetic testing.

Effective smMIPs-Based Sequencing of Maculopathy-Associated Genes in Stargardt Disease Cases and Allied Maculopathies from the UK.

Macular dystrophies are a group of individually rare but collectively common inherited retinal dystrophies characterised by central vision loss and loss of visual acuity. Single molecule Molecular Inversion Probes (smMIPs) have proved effective in identifying genetic variants causing macular dystrophy. Here, a previously established smMIPs panel tailored for genes associated with macular diseases has been used to examine 57 UK macular dystrophy cases, achieving a high solve rate of 63.2% (36/57). Among 27 bi-allelic STGD1 cases, only three novel ABCA4 variants were identified, illustrating that the majority of ABCA4 variants in Caucasian STGD1 cases are currently known. We examined cases with ABCA4-associated disease in detail, comparing our results with a previously reported variant grading system, and found this model to be accurate and clinically useful. In this study, we showed that ABCA4-associated disease could be distinguished from other forms of macular dystrophy based on clinical evaluation in the majority of cases (34/36).

Late-Onset Autosomal Dominant Macular Degeneration Caused by Deletion of the CRX Gene.

PURPOSE: To characterize the phenotype observed in a case series with macular disease and determine the cause. DESIGN: Multicenter case series. PARTICIPANTS: Six families (7 patients) with sporadic or multiplex macular disease with onset at 20 to 78 years, and 1 patient with age-related macular degeneration. METHODS: Patients underwent ophthalmic examination; exome, genome, or targeted sequencing; and/or polymerase chain reaction (PCR) amplification of the breakpoint, followed by cloning and Sanger sequencing or direct Sanger sequencing. MAIN OUTCOME MEASURES: Clinical phenotypes, genomic findings, and a hypothesis explaining the mechanism underlying disease in these patients. RESULTS: All 8 cases carried the same deletion encompassing the genes TPRX1, CRX, and SULT2A1, which was absent from 382 control individuals screened by breakpoint PCR and 13 096 Clinical Genetics patients with a range of other inherited conditions screened by array comparative genomic hybridization. Microsatellite genotypes showed that these 7 families are not closely related, but genotypes immediately adjacent to the deletion breakpoints suggest they may share a distant common ancestor. CONCLUSIONS: Previous studies had found that carriers for a single defective CRX allele that was predicted to produce no functional CRX protein had a normal ocular phenotype. Here, we show that CRX whole-gene deletion in fact does cause a dominant late-onset macular disease.

Novel SIX6 mutations cause recessively inherited congenital cataract, microcornea, and corneal opacification with or without coloboma and microphthalmia.

Purpose: To investigate the molecular basis of recessively inherited congenital cataract, microcornea, and corneal opacification with or without coloboma and microphthalmia in two consanguineous families. Methods: Conventional autozygosity mapping was performed using single nucleotide polymorphism (SNP) microarrays. Whole-exome sequencing was completed on genomic DNA from one affected member of each family. Exome sequence data were also used for homozygosity mapping and copy number variation analysis. PCR and Sanger sequencing were used to confirm the identification of mutations and to screen further patients. Evolutionary conservation of protein sequences was assessed using CLUSTALW, and protein structures were modeled using PyMol. Results: In family MEP68, a novel homozygous nucleotide substitution in SIX6 was found, c.547G>C, that converts the evolutionarily conserved aspartic acid residue at the 183rd amino acid in the protein to a histidine, p.(Asp183His). This residue mapped to the third helix of the DNA-binding homeobox domain in SIX6, which interacts with the major groove of double-stranded DNA. This interaction is likely to be disrupted by the mutation. In family F1332, a novel homozygous 1034 bp deletion that encompasses the first exon of SIX6 was identified, chr14:g.60975890_60976923del. Both mutations segregated with the disease phenotype as expected for a recessive condition and were absent from publicly available variant databases. Conclusions: Our findings expand the mutation spectrum in this form of inherited eye disease and confirm that homozygous human SIX6 mutations cause a developmental spectrum of ocular phenotypes that includes not only the previously described features of microphthalmia, coloboma, and congenital cataract but also corneal abnormalities.

Novel homozygous mutations in the transcription factor NRL cause non-syndromic retinitis pigmentosa.

Purpose: To describe the clinical phenotype and genetic basis of non-syndromic retinitis pigmentosa (RP) in one family and two sporadic cases with biallelic mutations in the transcription factor neural retina leucine zipper (NRL). Methods: Exome sequencing was performed in one affected family member. Microsatellite genotyping was used for haplotype analysis. PCR and Sanger sequencing were used to confirm mutations in and screen other family members where they were available. The SMART tool for domain prediction helped us build the protein schematic diagram. Results: For family MM1 of Pakistani origin, whole-exome sequencing and microsatellite genotyping revealed homozygosity on chromosome 14 and identified a homozygous stop-loss mutation in NRL, NM_006177.5: c.713G>T, p.*238Lext57, which is predicted to add an extra 57 amino acids to the normal protein chain. The variant segregated with disease symptoms in the family. For case RP-3051 of Spanish ancestry, clinical exome sequencing focusing on the morbid genome highlighted a homozygous nonsense mutation in NRL, c.238C>T, p.Gln80*, as the most likely disease candidate. For case RP-1553 of Romanian ethnicity, targeted-exome sequencing of 73 RP/LCA genes identified a homozygous nonsense mutation in NRL, c.544C>T, p.Gln182*. The variants were either rare or absent in the gnomAD database. Conclusions: NRL mutations predominantly cause dominant retinal disease, but there have been five published reports of mutations causing recessive disease. Here, we present three further examples of recessive RP due to NRL mutations. The phenotypes observed are consistent with those in the previous reports, and the observed mutation types and distribution further confirm distinct patterns for variants in NRL causing recessive and dominant diseases.

Uncovering the burden of hidden ciliopathies in the 100 000 Genomes Project: a reverse phenotyping approach.

BACKGROUND: The 100 000 Genomes Project (100K) recruited National Health Service patients with eligible rare diseases and cancer between 2016 and 2018. PanelApp virtual gene panels were applied to whole genome sequencing data according to Human Phenotyping Ontology (HPO) terms entered by recruiting clinicians to guide focused analysis. METHODS: We developed a reverse phenotyping strategy to identify 100K participants with pathogenic variants in nine prioritised disease genes (BBS1, BBS10, ALMS1, OFD1, DYNC2H1, WDR34, NPHP1, TMEM67, CEP290), representative of the full phenotypic spectrum of multisystemic primary ciliopathies. We mapped genotype data 'backwards' onto available clinical data to assess potential matches against phenotypes. Participants with novel molecular diagnoses and key clinical features compatible with the identified disease gene were reported to recruiting clinicians. RESULTS: We identified 62 reportable molecular diagnoses with variants in these nine ciliopathy genes. Forty-four have been reported by 100K, 5 were previously unreported and 13 are new diagnoses. We identified 11 participants with unreportable, novel molecular diagnoses, who lacked key clinical features to justify reporting to recruiting clinicians. Two participants had likely pathogenic structural variants and one a deep intronic predicted splice variant. These variants would not be prioritised for review by standard 100K diagnostic pipelines. CONCLUSION: Reverse phenotyping improves the rate of successful molecular diagnosis for unsolved 100K participants with primary ciliopathies. Previous analyses likely missed these diagnoses because incomplete HPO term entry led to incorrect gene panel choice, meaning that pathogenic variants were not prioritised. Better phenotyping data are therefore essential for accurate variant interpretation and improved patient benefit.

Bi-allelic mutation of CTNNB1 causes a severe form of syndromic microphthalmia, persistent foetal vasculature and vitreoretinal dysplasia.

BACKGROUND: Inherited vitreoretinopathies arise as a consequence of congenital retinal vascularisation abnormalities. They represent a phenotypically and genetically heterogeneous group of disorders that can have a major impact on vision. Several genes encoding proteins and effectors of the canonical Wnt/ß-catenin pathway have been associated and precise diagnosis, although difficult, is essential for proper clinical management including syndrome specific management where appropriate. This work aimed to investigate the molecular basis of disease in a single proband born to consanguineous parents, who presented with microphthalmia, persistent foetal vasculature, posterior lens vacuoles, vitreoretinal dysplasia, microcephaly, hypotelorism and global developmental delay, and was registered severely visually impaired by 5 months of age. METHODS: Extensive genomic pre-screening, including microarray comparative genomic hybridisation and sequencing of a 114 gene panel associated with cataract and congenital ophthalmic disorders was conducted by an accredited clinical laboratory. Whole exome sequencing (WES) was undertaken on a research basis and in vitro TOPflash transcriptional reporter assay was utilised to assess the impact of the putative causal variant. RESULTS: In the proband, WES revealed a novel, likely pathogenic homozygous mutation in the cadherin-associated protein beta-1 gene (CTNNB1), c.884C>G; p.(Ala295Gly), which encodes a co-effector molecule of the Wnt/ß-catenin pathway. The proband's parents were shown to be heterozygous carriers but ophthalmic examination did not detect any abnormalities. Functional assessment of the missense variant demonstrated significant reduction of ß-catenin activity. CONCLUSIONS: This is the first report of a biallelic disease-causing variation in CTNNB1. We conclude that this biallelic, transcriptional inactivating mutation of CTNNB1 causes a severe, syndromic form of microphthalmia, persistent foetal vasculature and vitreoretinal dysplasia that results in serious visual loss in infancy.

Molecular diagnoses in the congenital malformations caused by ciliopathies cohort of the 100,000 Genomes Project.

BACKGROUND: Primary ciliopathies represent a group of inherited disorders due to defects in the primary cilium, the 'cell's antenna'. The 100,000 Genomes Project was launched in 2012 by Genomics England (GEL), recruiting National Health Service (NHS) patients with eligible rare diseases and cancer. Sequence data were linked to Human Phenotype Ontology (HPO) terms entered by recruiting clinicians. METHODS: Eighty-three prescreened probands were recruited to the 100,000 Genomes Project suspected to have congenital malformations caused by ciliopathies in the following disease categories: Bardet-Biedl syndrome (n=45), Joubert syndrome (n=14) and 'Rare Multisystem Ciliopathy Disorders' (n=24). We implemented a bespoke variant filtering and analysis strategy to improve molecular diagnostic rates for these participants. RESULTS: We determined a research molecular diagnosis for n=43/83 (51.8%) probands. This is 19.3% higher than previously reported by GEL (n=27/83 (32.5%)). A high proportion of diagnoses are due to variants in non-ciliopathy disease genes (n=19/43, 44.2%) which may reflect difficulties in clinical recognition of ciliopathies. n=11/83 probands (13.3%) had at least one causative variant outside the tiers 1 and 2 variant prioritisation categories (GEL's automated triaging procedure), which would not be reviewed in standard 100,000 Genomes Project diagnostic strategies. These include four structural variants and three predicted to cause non-canonical splicing defects. Two unrelated participants have biallelic likely pathogenic variants in LRRC45, a putative novel ciliopathy disease gene. CONCLUSION: These data illustrate the power of linking large-scale genome sequence to phenotype information. They demonstrate the value of research collaborations in order to maximise interpretation of genomic data.

A multi-ethnic genome-wide association study implicates collagen matrix integrity and cell differentiation pathways in keratoconus.

Keratoconus is characterised by reduced rigidity of the cornea with distortion and focal thinning that causes blurred vision, however, the pathogenetic mechanisms are unknown. It can lead to severe visual morbidity in children and young adults and is a common indication for corneal transplantation worldwide. Here we report the first large scale genome-wide association study of keratoconus including 4,669 cases and 116,547 controls. We have identified significant association with 36 genomic loci that, for the first time, implicate both dysregulation of corneal collagen matrix integrity and cell differentiation pathways as primary disease-causing mechanisms. The results also suggest pleiotropy, with some disease mechanisms shared with other corneal diseases, such as Fuchs endothelial corneal dystrophy. The common variants associated with keratoconus explain 12.5% of the genetic variance, which shows potential for the future development of a diagnostic test to detect susceptibility to disease.

The Role of Csmd1 during Mammary Gland Development.

The Cub Sushi Multiple Domains-1 (CSMD1) protein is a tumour suppressor which has been shown to play a role in regulating human mammary duct development in vitro. CSMD1 knockdown in vitro demonstrated increased cell proliferation, invasion and motility. However, the role of Csmd1 in vivo is poorly characterised when it comes to ductal development and is therefore an area which warrants further exploration. In this study a Csmd1 knockout (KO) mouse model was used to identify the role of Csmd1 in regulating mammary gland development during puberty. Changes in duct development and protein expression patterns were analysed by immunohistochemistry. This study identified increased ductal development during the early stages of puberty in the KO mice, characterised by increased ductal area and terminal end bud number at 6 weeks. Furthermore, increased expression of various proteins (Stat1, Fak, Akt, Slug/Snail and Progesterone receptor) was shown at 4 weeks in the KO mice, followed by lower expression levels from 6 weeks in the KO mice compared to the wild type mice. This study identifies a novel role for Csmd1 in mammary gland development, with Csmd1 KO causing significantly more rapid mammary gland development, suggesting an earlier adult mammary gland formation.

New variants and in silico analyses in GRK1 associated Oguchi disease.

Biallelic mutations in G-Protein coupled receptor kinase 1 (GRK1) cause Oguchi disease, a rare subtype of congenital stationary night blindness (CSNB). The purpose of this study was to identify disease causing GRK1 variants and use in-depth bioinformatic analyses to evaluate how their impact on protein structure could lead to pathogenicity. Patients' genomic DNA was sequenced by whole genome, whole exome or focused exome sequencing. Disease associated variants, published and novel, were compared to nondisease associated missense variants. The impact of GRK1 missense variants at the protein level were then predicted using a series of computational tools. We identified twelve previously unpublished cases with biallelic disease associated GRK1 variants, including eight novel variants, and reviewed all GRK1 disease associated variants. Further structure-based scoring revealed a hotspot for missense variants in the kinase domain. In addition, to aid future clinical interpretation, we identified the bioinformatics tools best able to differentiate disease associated from nondisease associated variants. We identified GRK1 variants in Oguchi disease patients and investigated how disease-causing variants may impede protein function in-silico.

DYNC2H1 hypomorphic or retina-predominant variants cause nonsyndromic retinal degeneration.

PURPOSE: Determining the role of DYNC2H1 variants in nonsyndromic inherited retinal disease (IRD). METHODS: Genome and exome sequencing were performed for five unrelated cases of IRD with no identified variant. In vitro assays were developed to validate the variants identified (fibroblast assay, induced pluripotent stem cell [iPSC] derived retinal organoids, and a dynein motility assay). RESULTS: Four novel DYNC2H1 variants (V1, g.103327020_103327021dup; V2, g.103055779A>T; V3, g.103112272C>G; V4, g.103070104A>C) and one previously reported variant (V5, g.103339363T>G) were identified. In proband 1 (V1/V2), V1 was predicted to introduce a premature termination codon (PTC), whereas V2 disrupted the exon 41 splice donor site causing incomplete skipping of exon 41. V1 and V2 impaired dynein-2 motility in vitro and perturbed IFT88 distribution within cilia. V3, homozygous in probands 2-4, is predicted to cause a PTC in a retina-predominant transcript. Analysis of retinal organoids showed that this new transcript expression increased with organoid differentiation. V4, a novel missense variant, was in trans with V5, previously associated with Jeune asphyxiating thoracic dystrophy (JATD). CONCLUSION: The DYNC2H1 variants discussed herein were either hypomorphic or affecting a retina-predominant transcript and caused nonsyndromic IRD. Dynein variants, specifically DYNC2H1 variants are reported as a cause of non syndromic IRD.

The ciliary Frizzled-like receptor Tmem67 regulates canonical Wnt/ß-catenin signalling in the developing cerebellum via Hoxb5.

Primary cilia defects result in a group of related pleiotropic malformation syndromes known as ciliopathies, often characterised by cerebellar developmental and foliation defects. Here, we describe the cerebellar anatomical and signalling defects in the Tmem67tm1(Dgen)/H knockout mouse. At mid-gestation, Tmem67 mutant cerebella were hypoplastic and had aberrantly high canonical Wnt/ß-catenin signalling, proliferation and apoptosis. Later in development, mutant cerebellar hemispheres had severe foliation defects and inferior lobe malformation, characterized by immature Purkinje cells (PCs). Early postnatal Tmem67 mutant cerebellum had disrupted ciliogenesis and reduced responsiveness to Shh signalling. Transcriptome profiling of Tmem67 mutant cerebella identified ectopic increased expression of homeobox-type transcription factors (Hoxa5, Hoxa4, Hoxb5 and Hoxd3), normally required for early rostral hindbrain patterning. HOXB5 protein levels were increased in the inferior lobe, and increased canonical Wnt signalling, following loss of TMEM67, was dependent on HOXB5. HOXB5 occupancy at the ß-catenin promoter was significantly increased by activation of canonical Wnt signalling in Tmem67-/- mutant cerebellar neurones, suggesting that increased canonical Wnt signalling following mutation or loss of TMEM67 was directly dependent on HOXB5. Our results link dysregulated expression of Hox group genes with ciliary Wnt signalling defects in the developing cerebellum, providing new mechanistic insights into ciliopathy cerebellar hypoplasia phenotypes.

Correction: Biallelic sequence and structural variants in RAX2 are a novel cause for autosomal recessive inherited retinal disease.

The original version of this Article contained an incorrect version of Fig. 3, which included two variants initially shown in black text in Fig. 3a that the authors removed from the final manuscript. The correct version of Fig. 3 without the two variants now appears in the PDF and HTML versions of the Article.

Biallelic sequence and structural variants in RAX2 are a novel cause for autosomal recessive inherited retinal disease.

PURPOSE: RAX2 encodes a homeobox-containing transcription factor, in which four monoallelic pathogenic variants have been described in autosomal dominant cone-dominated retinal disease. METHODS: Exome sequencing in a European cohort with inherited retinal disease (IRD) (n = 2086) was combined with protein structure modeling of RAX2 missense variants, bioinformatics analysis of deletion breakpoints, haplotyping of RAX2 variant c.335dup, and clinical assessment of biallelic RAX2-positive cases and carrier family members. RESULTS: Biallelic RAX2 sequence and structural variants were found in five unrelated European index cases, displaying nonsyndromic autosomal recessive retinitis pigmentosa (ARRP) with an age of onset ranging from childhood to the mid-40s (average mid-30s). Protein structure modeling points to loss of function of the novel recessive missense variants and to a dominant-negative effect of the reported dominant RAX2 alleles. Structural variants were fine-mapped to disentangle their underlying mechanisms. Haplotyping of c.335dup in two cases suggests a common ancestry. CONCLUSION: This study supports a role for RAX2 as a novel disease gene for recessive IRD, broadening the mutation spectrum from sequence to structural variants and revealing a founder effect. The identification of biallelic RAX2 pathogenic variants in five unrelated families shows that RAX2 loss of function may be a nonnegligible cause of IRD in unsolved ARRP cases.