- Research article
- Open Access
Applying next generation sequencing with microdroplet PCR to determine the disease-causing mutations in retinal dystrophies
BMC Ophthalmology volume 17, Article number: 157 (2017)
Inherited Retinal dystrophy (IRD) is a broad group of inherited retinal disorders with heterogeneous genotypes and phenotypes. Next generation sequencing (NGS) methods have been broadly applied for analyzing patients with IRD. Here we report a novel approach to enrich the target gene panel by microdroplet PCR.
This assay involved a primer library which targeted 3071 amplicons from 2078 exons comprised of 184 genes involved in retinal function and/or retinal development. We amplified the target regions using the RainDance target enrichment PCR method and sequenced the products using the MiSeq NGS platform.
In this study, we analyzed 82 samples from 67 families with IRD. Bioinformatics analysis indicated that this procedure was able to reach 99% coverage of target sequences with an average sequence depth of reads at 119×. The variants detected by this study were filtered, validated, and prioritized by pathogenicity analysis. Genotypes and phenotypes were correlated by determining a consistent relationship in 38 propands (56.7%). Pathogenic variants in genes related to retinal function were found in another 11 probands (16.4%), but the clinical correlations showed inconsistencies and insufficiencies in these patients.
The application of NGS in IRD clinical molecular diagnosis provides a powerful approach to exploring the etiology and pathology in patients. It is important for the clinical laboratory to interpret the molecular findings in the context of patient clinical presentations because accurate interpretation of pathogenic variants is critical for delivering solid clinical molecular diagnosis to clinicians and patients and improving the standard care of patients.
Inherited Retinal dystrophy (IRD) is a clinically and genetically heterogeneous group of disorders that affects approximately 1 in 2000 individuals worldwide . Clinically IRD can be classified according to the clinical phenotype and course of the disease, as well as the underlying inheritance. However, sporadic retinal dystrophies, such as Retinitis Pigmentosa (RP), Leber Congenital Amaurosis (LCA), Cone and Cone-Rod dystrophy (CRD), Macular Dystrophy (MD), and Congenital Stationary Night Blindness (CSNB) can present difficulties in the clinical differential diagnosis. The heterogeneity of some syndromic conditions with retinal abnormalities (i.e., Usher syndrome, Bardet-Biedl syndrome, and Nephronophthisis) could add more complexity because these diseases are caused by multiple underlying gene defects [2, 3]. In molecular pathology, IRD may result from mutations in a variety of genes and may show different inheritance patterns, including autosomal dominant, autosomal recessive, X-linked, and mitochondrial inheritance [3, 4]. Additionally, mutations within the same gene may be associated with different phenotypes [2, 5]. Although monogenic forms have been reported in most families, some digenic forms have also been identified [2, 3].
Retinitis pigmentosa is the most frequent type of IRD and it is also highly genetic heterogeneous. More than 50 genes have been identified responsible for RP with similar clinical presentations. Typically, it presents night blindness and peripheral visual field decreasing. As part of progression in RP, nigh blindness will eventually become blindness. On the other hand, the night blindness is also a typical presentation of CSDB. While CSDB can be distinguished by electroretinogram (ERG), more than 10 genes have been identified responsible for CSDB and a portion of CSDB patients were found no mutation in the known CSDB genes. LCA, CRD, MD, and Occult Macular Dystrophy (OCMD) can be differentiated clinically in general, but the molecular pathology of malfunction in cone photoreceptors cannot be established without identifying the disease-causing gene in each patient.
To date, more than 200 genes have been identified for different retinal diseases by a variety of methods (http://www.sph.uth.tmc.edu/retnet/sum-dis.htm). Significant progress has been made in determining the molecular causes of IRD, but much more work remains to be done.
Genetic testing for IRD can identify causative mutations but requires sequencing of many individual candidate genes. Next generation sequencing (NGS) has been broadly applied for analyzing patients with IRD. NGS enables rapid and cost-effective parallel sequencing of a large panel of disease genes. It offers an ideal model in a clinical diagnostic setting. Many studies investigating sequence-capture technology with selected gene targets for enrichment have been published [6,7,8,9,10]. A significantly higher rate of molecular diagnosis (as compared to Sanger sequencing potential candidate genes) of well above 50% was also achieved. Furthermore, NGS can be adapted to include other retinal diseases without significantly increasing the cost. In this study, we performed a comprehensive molecular screening of patients by NGS using Illumina MiSeq platform and target enrichment using a customized primer library based microdroplet PCR. We accomplished targeted gene panel sequencing of 184 genes in 74 patients with different inheritance patterns and clinical diagnoses including Choroidermia (CHM), CRD, CSNB, LCA, MD, OCMD, and RP. Careful clinical evaluation and follow-up led to the more precise clinical diagnoses and extensive phenotyping in these RD families. Collectively, this study underscores the importance of combining comprehensive molecular screening and clinical information to accurately diagnose diverse retinal disorders.
Subjects and clinical investigation
A set of 82 genomic DNA samples was included in this study. This cohort includes 74 patients with clinical diagnosis of Choroidermia (CHM), CRD, CSNB, LCA, MD, Occult Macular Dystrophy (OCMD), and RP from the NEI Ophthalmic Genetics Clinic or tested in the NEI DNA Diagnostic Laboratory. Clinical evaluations included measurement of best-corrected visual acuity, visual field assessment (either kinetic or static perimetry depending on patient presentation), fundus biomicroscopy and indirect ophthalmoscopy. Optical coherence tomography (Stratus OCT 3 or Cirrus, Carl Zeiss Meditec, Inc., Dublin, CA) was performed. Color and autofluorescence fundus imaging was obtained. International Society for Clinical Electrophysiology of Vision (ISCEV) standard electroretinography responses were obtained using a Burian-Allen contact lens electrode and an LKC-2000 system (LKC Technologies, Gaithersburg, MD). Systemic examinations were performed whenever necessary. Genomic DNA was isolated from peripheral leukocytes using the Gentra Puregene Blood Kit (Qiagen Inc., Valencia, CA) in accordance with the manufacturer’s protocol. Whenever available, a blood sample from affected and unaffected family members was collected for co-segregation analysis. Pedigrees were constructed based on patient interviews. This study was reviewed and approved by the Combined Neuroscience Institutional Review Board of the National Institutes of Health and informed consent was obtained from each participant as adhering to tenets of the Declarations of Helsinki.
Molecular genetic analysis
Design of the microdroplet-based PCR primer library
A panel of 184 RD genes responsible for retinal related diseases was provided to RainDance Technologies Inc. (Billerica, MA) for the design of a primer library in 2011, using their custom primer design pipeline based on the Primer3 algorithm (Additional file 1: Table S1) (http://frodo.wi.mit.edu/primer3). A total of 3071 primer pairs were designed to target the 2078 coding exons of the184 RD genes (Additional file 2: Table S2). RD genes included in this study were previously associated with RD in the literature, RetNet database (http://www.sph.uth.tmc.edu/retnet/sum-dis.htm) and the OMIM database  (http://www.ncbi.nlm.nih.gov/omim).
RainDance target enrichment and NGS sequencing
Before amplification, the samples were fragmented to 5 kb by shearing the genomic DNA with the Covaris M220 instrument (Covaris, Woburn, MA) and the preparation for amplification was following the manufacturer’s recommended protocol . The samples (50 ng) then entered the sample preparation protocol (NEXTflex PCR-Free DNA Sequencing Kit and NEXTflex PCR-Free barcode 1, Illumina®-Compatible, BIOO Scientific, Austin, Texas), followed by 150 bp pair-end sequencing on a MiSeq instrument (Illumina, San Diego, CA). Samples in each batch were indexed using 12 different index tags (Nextera, Epicentre Biotechnologies, Madison, WI).
Reference gene sequences were annotated with known single-nucleotide polymorphisms (SNPs) from the NCBI dbSNP database build 130, 1000 Genomes database build 201,105 and 201,011 and mutations from the Human Gene Mutation Database (HGMD) (http://www.hgmd.cf.ac.uk/) or Leiden Open Variation Database (LOVD) (http://www.lovd.nl/3.0) or reported in the literature. Other variant reference databases used were following: HapMap variants database (Ensembl), Clinical Variants in dbSNP database (NCBI), 1000 genomes database (Ensembl), Chromosome bands ideogram (UCSC), dbSNP common variants (UCSC), Genomic Annotations (Ensembl), Cosmic Noncoding Variants & Coding Mutations (http://cancer.sanger.ac.uk/cosmic).
The sequencing reads were aligned to the reference sequence from NCBI (hg19 Build 37) and analyzed using the CLC Bio Genomics Workbench software™ (QIAGEN). Data was analyzed by using the Probabilistic Variant Detection Tool provided by the Genomics Workbench. Mutation calls were made under robust and stringent sequencing criteria set by the CLC Genomics Workbench with a probability call of 100 in order to delete any false mutation observations. Identified sequence variants were annotated according to the guidelines published by the Human Genome Variation Society (HGVS).
Determination of pathogenic variants
The potential pathogenicity of sequence variants was assessed as follows. First, reported SNP with allele frequency above 2% were excluded. Second, potential damaging nucleotide substitution variants defined as nonsense, missense, silent with predicted splicing effects, intronic splice-site variants and any novel variants in the target regions were included. Damaging indels located in the coding regions of the gene were also included. The variants that met the above criteria were selected for downstream analysis as candidate alleles. The pathogenic effect of every candidate allele on gene expression and protein function was assessed in silico using the Alamut Visual software version 2.7 (Alamut, Rouen, France) (www.interactivebiosoftware.com). The pathogenicity assessment includes PolyPhen-2 (Polymorphism Phenotyping) (http://genetics.bwh.harvard.edu/pph2/), SIFT (http://sift.jcvi.org/), and Mutation Taster (http://www.mutationtaster.org/). Synonymous and intronic sequence variants were assessed for potential deleterious effects upon messenger RNA splicing using the Human Splicing Finder V.2.4 tool (http://www.umd.be/HSF/). The Alamut analysis generated a report per variant for the final interpretation.
Validation of mutations by Sanger sequencing and segregation
All pathogenic variants identified by NGS after variant filtering were validated by Sanger sequencing as previously described . Reference DNA sequences were obtained from the NCBI. Primer3 (http://frodo.wi.mit.edu/primer3/) was used to design new primers wherever applicable otherwise primers that were used in the RainDance Inc. designed primer library were used in the study.
A total of 82 samples including 74 patients and 4 unaffected family members from 67 unrelated RD families were included in this study. In addition, 3 mutation positive samples that were previously studied in this lab and one DNA control sample acquired from Coriell Repository were included in this study . Clinical examinations were performed in the NEI Ophthalmic Genetics Clinic and patients were clinically diagnosed with CHM, CRD, CSNB, LCA, MD, OCMD, and RP. Patient samples referred from outside clinics were previously tested in the NEI DNA Diagnostic Laboratory.
Evaluation and validation of the novel RD gene panel
NGS sequencing data were processed and analyzed through a bioinformatic pipeline by the Genomic Workbench and the Alamut Visual. An average of 1.14 million reads was generated per sample, an average of 59.5% of reads mapped to the targeted regions. Bioinformatic analysis indicated that this procedure was able to reach 99% of coverage on target sequences with averages of 97% at 1X, 95% at 10X, and 90% at 20X. The average read depth was 119X. The read depth per sample was evaluated at 1X, 10X, and 20X and summarized in the Fig. 1.
Variant filtering and annotation
On average, our bioinformatic pipeline generated 500–600 variants per sample. After filtering out common polymorphisms with a frequency > 2% in any of the variant databases queried, including 1000 Genome build 201,105 and 201,011 (Genomes Project 2010) and dbSNP130 (NCBI dbSNP database build 130), an average around 30–60 variants remained per sample. The remaining variants were further filtered based on our secondary functional selection as stopgain or stoploss variants, missense variants, synonymous variants with predicted splicing effects, and intronic variants within the 25 bp of exon flanking regions with predicted splicing effects. We had an average of 10 candidate variants per sample for further Alamut pathogenicity analysis and reporting.
Identification of known and novel mutations
Sanger sequencing was used to analyze whether a variation segregated with the disease phenotype in the corresponding family (whenever additional family members were available). Using this approach, we confirmed 63 previously determined pathogenic variants and 28 novel variants predicted as pathogenic or likely pathogenic in our analyses (Additional file 3: Table S3).
The clinical correlation was performed after mutations and variations were determined and verified by Sanger sequencing. We found that phenotype-genotype relationship could be established in 38 of the 67 probands (56.7%) and the additional affected family members were confirmed with the familial mutations (Table 1). We included 6 samples from a large Occult Macular Dystrophy (OCMD) family. In the three affected members, we confirmed the known mutation, p.R45W in the RP1L1 gene and excluded the potential of additional pathogenic variants in other retinal disease related genes. We also confirmed the mutation carrying status in the other three unaffected family members.
We found inconsistences between the genotypes and the clinical diagnoses in 4 samples (Table 2). The confirmed mutations in these 4 patients had been reported with different retinal disorders for the genotype-phenotype relationship in the literature [13,14,15], but not reported for conditions that our patients were diagnosed with. Patient RD1–12 was diagnosed with sporadic CRD. Two C2orf71 gene pathogenic variants were confirmed: the p.W505* was reported in unrelated patient  and the p.S1090Ifs*17 was novel and not found in populations (ExAC: Exome Aggregation Consortium, Cambridge, MA) (http://exac.broadinstitute.org/). Additionally, three heterozygous stop-gain mutations in unrelated genes (BBS4, TYRP1, and SLC45A2) were found (Table 2). The C2orf71 gene mutations were reported in patients with autosomal recessive RP , not CRD. Our review of the patient’s clinical conditions still concluded a diagnosis of CRD (Fig. 2a and b). Patient RD11–05 was diagnosed with sporadic RP. Exam was positive for a preserved central visual field island with 20/32 acuity and some macular cystic changes. Electroretinography was consistent with retinitis pigmentosa. We identified a known mutation in the GUCA1A gene, p.P50L, which has been reported in patients with autosomal dominant CRD (Downes et al. 2001), not RP. Her clinical manifestation was consistent with RP (Fig. 2c). Patient RD14–05 was diagnosed with sporadic RP. Best corrected visual acuity measured at 20/100 and central visual field was limited to 10 degrees with a few peripheral islands. Findings were consistent with advanced RP (Fig. 2d and e). We found two variants in the TRPM1 gene, p.W398R and p.R1305H. The p.W398R was novel and was predicted as pathogenic, but the p.R1305 was found in African population with an allele frequency of 2% in our most recent ExAC search (rs13380059). Patient RD11–06 was also diagnosed with RP at age of 74 yr. and no family history was reported. Best corrected visual acuity was at 20/800. The patient reported a 10 year history of decline in peripheral and central visual function and denies any problems with vision as a child. Visual field was limited to a paracentral small island and electroretinography responses were extinguished (Fig. 2f). The patient was found carrying one candidate allele per gene in three genes including 2 reported mutations and 1 reported variant in the dbSNP database in heterozygous state. One of the reported mutations, p.P575L in the GUCY2D gene, has a very high allele frequency in the African population in ExAC (unlikely pathogenic). The other known mutation, p.P406L in the TYR gene, has been reported in patients with autosomal recessive Oculocutaneous Albinism Type I. The reported variant, IVS9 + 1G > A in the TYRP1 gene, is predicted as pathogenic for a likely splicing error, but mutations in the TYRP1 have been reported with autosomal recessive Oculocutaneous Albinism Type III.
We also found single heterozygous mutation status as insufficient to confirm the genotype-phenotype relationship in another 5 patients. The patients RD4–05 and RD11–02 were found heterozygous for novel EYS gene mutations, p.W2046* and p.Y1893Rfs*12 respectively. The mutations in the EYS gene have been reported with autosomal recessive RP. Since possible gene or partial gene copy number variation or deep intronic mutations were not excluded in this study for these patients, it was insufficient to conclude the correlation in these 2 patients. A similar situation presented itself for patients RD14–03 and RD15–01 who carried single heterozygous mutations in the ABCA4 gene, IVS + 5G > A splicing error and p.I1562T respectively (Table 2). It was reported that a portion of patients with the Stargardt disease were bearing single heterozygous mutations in the ABCA4 gene [17, 18]. It could be considered as consistent with the clinical diagnosis of RP, even if the second mutation in their ABCA4 gene was not identified. Adding further complexity to the interpretation, patient RD14–03 was also found carrying two additional USH2A gene variants, p.S2867 L and p.V3518I. Of these USH2A gene variants, the p.S2867 L has a very low allele frequency in populations (ExAC) (http://exac.broadinstitute.org/) and is predicted as pathogenic, while the p.V3518I has an allele frequency of 2% in the African population and 6 homozygous individual have been identified in the African population(ExAC) (http://exac.broadinstitute.org/). The disease-responsible gene for patient RD14–03 remains unidentified. The patient RD15–01 was found with additional pathogenic variant in the RPGRIP1 gene, p.P585S - previously reported in patient with autosomal recessive CRD .
We have also identified additional previously reported pathogenic variants in some of the probands for whom we established genotype-phenotype relationship in this study (Table 3). These additional pathogenic alleles could be sufficient for a clinical correlation by itself (RD14–08 and RD12–02), or have been reported multiple times in literature (RD20–07, RD12–02, RD12–01, RD13–02, RD12–07, and RD13–08, or are predicted as obligate disease causing variants (RD11–03). For example, the p.C984R in the GUCY2D gene was not reported in the literature and predicted as pathogenic. Mutations in the GUCY2D have been reported in autosomal dominant CRD. The p.C984R is a valid pathogenic variant by itself.
To technically verify the sequencing results, we included 4 previously tested DNA samples. Two of them were tested in the resequencing chip by Affimatrix, in which we had previously identified and validated one missense mutation in the GUCY2D gene (RD1–01) and no mutation in a Coriell DNA sample RD1–03, ND0068*B1 . The other 2 samples were from previous Sanger sequencing: one 1 bp deletion in the CNGB3 gene in one sample (RD13–08) and two OCA gene mutations in the other sample (RD13–07). All of the mutations were correctly identified through the bioinformatics pipeline. In addition to the previously identified mutation, we also detected additional pathogenic variants in the RD13–08 (Table 3).
Microdroplet PCR has been used in many studies including the sequencing of the entire human X chromosome exome with coverage of 97% . Our results demonstrated as much as 99% coverage and average read depth at 119X in our standardized procedure. We were testing the microdroplet PCR as the enrichment method for the targeted gene panel NGS with consideration of repeatability, standardization, and potential clinical application. Comparing with other enrichment technologies at the time, the microdroplet PCR provided a reliable and manageable procedure. As we observed random low coverage of a number of target exons with the capture based enrichment technologies in our studies and published data, we could address the uncovered target exons with a simultaneously pooled supplementary PCR using individually designed alternative primer pairs and unique PCR conditions. In addition, the designed primer library could be quickly adapted to verify the variants found in NGS by Sanger sequencing or non-amplifiable targets in target samples. The disadvantage of the microdroplet PCR included high up-front costs for the equipment and annual maintenance service contract even if there was no high volume sample flow.
After comparing several public and commercial data processing softwares, we chose the Qiagen Genomics Workbench as our primary sequence aligner and annotation tool. For quality control purposes, beside the general alignment and coverage analysis, this software provided a detailed amplicon based analysis table including reads per target region (amplicon), percent sequence coverage per amplicon, and depth of reads as zero covered target bases, reads at minimally covered bases, reads at maximally covered and averages of reads at non-zero covered sequences, which easily determined the amplification resistant targets. It enabled a standardized protocol with clear information on exact coverage limitations for every sample. As a microdroplet PCR based NGS, we obtained a comprehensive list of targets with no amplification. It will enable 100% coverage with follow-up Sanger sequencing (or other approaches) on targets with poor coverage quality if necessary. Although we could not extend the analysis to the copy number variation with large deletions/insertions beyond the size of amplicons in this study, we are still looking for the potential solution by implementing other aligner and calculating Algorithms.
Many NGS studies have been published recently to report pathogenic variant found in IRD patients [6,7,8,9,10]. NGS has the power to screen hundreds of genes in one reaction thereby increasing of finding incidental pathogenic variants. It will be essential to analyze the clinical validity by clinical correlation, family mutation co-segregation analysis, or other approaches to determine if the pathogenic alleles were indeed the disease-causing mutations per patient. Reporting pathogenic variants does not necessarily provide direct answers to the question regarding patient specific disease-causing mutations. Making studies available in literature could help other clinical geneticists by providing real cases with consistent clinical correlations and more importantly, detailed discussions of inconsistent/insufficient correlations.
In this study, we had the chance to evaluate pathogenic variants in the context of clinical presentations in these patients. We determined that 38 out of 67 probands had consistent clinical correlations (56.7%, Table 1). Segregation analysis and testing unaffected family members was performed wherever possible. Patients RD13–01 and RD14–04 were CRD affected son and affected mom. They were both heterozygous for the GUCA1A p.Y99C. Patients RD2–01 and RD2–02 are affected mom and daughter with peripapillary atrophy. There are two additional affected males in this family indicating a dominant inheritance (Leys et al., personal communication, manuscript is in preparation separately). It was interesting to note that patients RD10–02, RD10–03, and RD10–06 were three affected OCMD belonging to a single large family with clear autosomal dominant inheritance (Fig. 3). Six family members were included in this study. The other three family members (RD10–01, RD10–04, and RD10–05) were determined as RP1L1 p.R45W pathogenic variant positive and no clinical presentation of OCMD at the time was examined in the past. We included them to test the potential contribution of other pathogenic variants in the retinal genes. No additional pathogenic variant was found in the three affected patients in this family and no apparent co-segregation of any variants within two groups (affected three vs unaffected three). This result is consistent with the previously reported findings of incomplete penetrance of the p.R45W mutation in OCMD [20, 21]. The potential genetic modification for an in-family variation was unlikely from other retinal function related genes included in this study (Additional file 1: Table S1). We suspected that the p.R45W mutation might rather be a gain-of-function from malfunctioned products by an interaction with unknown component, which could lead to late onset macular dystrophy. The modification could be from a restricted expression of the wild type products too.
We found USH2A p.C759F mutation in a couple of RP families. Two patients, RD20–08 and RD20–02 were from a family with confirmed compound heterozygosity of the most common USH2A gene mutation, p.E766Sfs*21 and the p.C759F. Carrier status of the p.C759F was also confirmed in one available parent (RD20–01) in this family. In the other family, patient RD5–07 was apparently homozygous of the p.C759F. The USH2A gene mutations have been reported in the Usher syndrome and non-syndromic RP [22, 23]. The p.C759F was presumed as retinal disease specific allele . It was consistent in this study that the p.C759F correlated with Nonsyndromic RP. However, the pathogenicity of homozygous p.C759F was challenged by a recent clinical report that an unaffected and one affected dizygotic twin brothers shares homozygosity of p.C759F .
Many studies report putative pathogenic variants in genes that lack established correlation with the patients’ clinical manifestations [6, 15]. The inconsistencies surely raise the question each time about the correlation and deserve a careful review. Making these analyses available to our community would be of great interest for clinical diagnostic practice because of the impact on the interpretation for the clinical diagnostic laboratories. With the capacity of NGS, we will encounter cases with similar scenarios from time to time. In this study, the pathogenic variants overrepresented in patient RD1–12 could not be explained with the clinical manifestations of CRD (Table 2 and Fig. 2a and b). Two EYS gene pathogenic variants would be interpreted as putative disease causing in patients with RP as a clinical diagnosis, but those do not correlate with the CRD diagnosis in the patient RD1–12. On the other hand, we don’t know whether the additional pathogenic variants could contribute to the etiology or pathology of the clinical presentation in this patient. The patient RD11–05 was diagnosed with RP (Fig. 2 and Table 2) and we only found the GUCA1A p.P50L, which was reported in a family with autosomal dominant CRD and marked variability in expressivity , not RP. The ExAC reported an allele frequency of 0.1%, which was too high for an autosomal dominant inheritance. It could be possible that the p.P50L was coincidentally co-segregating with the CRD in the published CRD family if we disqualified the pathogenicity. In our patient RD11–05, it again could coincidentally co-segregate with RP. We should pay attention to the correlation if more cases with the same segregation can be found. The mutations in the TRPM1 gene have been reported in autosomal recessive Congenital Stationary Night Blindness . The patient RD14–05 was reviewed for clinical presentations and we could not correlate the two TRPM1 gene variants to her clinical manifestations (Table 3). The inconclusive clinical correlation did not support the predicted pathogenicity. For the patient RD11–06, the high frequent GUCY2D gene p.P575L allele does not support the pathogenicity. Two mutations in two different OCA related genes have no direct clinical relevance to the RP diagnosis in this patient. In our designing of this gene panel, we included 7 genes related to OCA. It was interesting to point out that we had observed relative high frequencies of OCA related alleles in heterozygous in 10 samples of this cohort (total 12 alleles were all in patient samples, not in the controls or family members except the OCA control sample) (Additional file 4: Table S4). We were not sure if this was just by chance or there was indeed an overrepresentation of OCA related alleles in IRD patients with certain phenotypes.
In our review of the clinical diagnosis of the 11 patients with inconclusive pathogenic variants in Table 2, we noticed comparable genotypes in some of the patients listed in Table 1 with the additional pathogenic variants. The combined genotypes were summarized in Table 3. If we made an assumption that the patients in Table 2 could have disease responsible genes not included in this study, those pathogenic variants we found in Table 2 might represent as the additional pathogenic variants like those in Table 3. It would be difficult to distinguish the two possibilities. We would assume that the likelihood in a patient with additional pathogenic variant being found and the disease-causing mutation being not found is much less likely based on overall positive clinical correlations in all studies. Interestingly, we found that one of the control samples also carried additional pathogenic variants in ABCA4 gene and RPGRIP1 gene and this sample was from a patient who was diagnosed with Achromatopsia (Table 3, RD13–08).
In addition, it is too early to suggest any pathogenic contribution from these additional pathogenic variants. In the 18 probands that had no potential pathogenic variant identified in this study, one patient was clinically diagnosed with Gyrate Atrophy and another one with Leber’s Hereditary Optic Neuropathy in our follow-up clinical review.
In summary, we tested a 184 targeted gene panel sequencing in 74 patients with different inheritance patterns and clinical diagnoses. We demonstrated that next-generation sequencing can be an effective tool for determining the pathogenic variants in inherited disease families with highly heterogeneous causes. We highlight the importance of interpreting molecular findings in the context of patient clinical presentation.
Single nucleotide polymorphism
Sohocki MM, Daiger SP, Bowne SJ, Rodriquez JA, Northrup H, Heckenlively JR, Birch DG, Mintz-Hittner H, Ruiz RS, Lewis RA, et al. Prevalence of mutations causing retinitis pigmentosa and other inherited retinopathies. Hum Mutat. 2001;17(1):42–51.
Daiger SP. Identifying retinal disease genes: how far have we come, how far do we have to go? Novartis Found Symp. 2004;255:17–27. discussion 27-36, 177-178
Rivolta C, Sharon D, DeAngelis MM, Dryja TP. Retinitis pigmentosa and allied diseases: numerous diseases, genes, and inheritance patterns. Hum Mol Genet. 2002;11(10):1219–27.
Hamel CP. Cone rod dystrophies. Orphanet J Rare Dis. 2007;2:7.
MacDonald IM, Tran M, Musarella MA. Ocular genetics: current understanding. Surv Ophthalmol. 2004;49(2):159–96.
Ge Z, Bowles K, Goetz K, Scholl HP, Wang F, Wang X, Xu S, Wang K, Wang H, Chen R. NGS-based Molecular diagnosis of 105 eyeGENE((R)) probands with Retinitis Pigmentosa. Sci Rep. 2015;5:18287.
Glockle N, Kohl S, Mohr J, Scheurenbrand T, Sprecher A, Weisschuh N, Bernd A, Rudolph G, Schubach M, Poloschek C, et al. Panel-based next generation sequencing as a reliable and efficient technique to detect mutations in unselected patients with retinal dystrophies. Eur J Hum Genet. 2014;22(1):99–104.
Neveling K, Collin RW, Gilissen C, van Huet RA, Visser L, Kwint MP, Gijsen SJ, Zonneveld MN, Wieskamp N, de Ligt J, et al. Next-generation genetic testing for retinitis pigmentosa. Hum Mutat. 2012;33(6):963–72.
O'Sullivan J, Mullaney BG, Bhaskar SS, Dickerson JE, Hall G, O'Grady A, Webster A, Ramsden SC, Black GC. A paradigm shift in the delivery of services for diagnosis of inherited retinal disease. J Med Genet. 2012;49(5):322–6.
Shanks ME, Downes SM, Copley RR, Lise S, Broxholme J, Hudspith KA, Kwasniewska A, Davies WI, Hankins MW, Packham ER, et al. Next-generation sequencing (NGS) as a diagnostic tool for retinal degeneration reveals a much higher detection rate in early-onset disease. Eur J Hum Genet. 2013;21(3):274–80.
Mondal K, Shetty AC, Patel V, Cutler DJ, Zwick ME. Targeted sequencing of the human X chromosome exome. Genomics. 2011;98(4):260–5.
Song J, Smaoui N, Ayyagari R, Stiles D, Benhamed S, MacDonald IM, Daiger SP, Tumminia SJ, Hejtmancik F, Wang X. High-throughput retina-array for screening 93 genes involved in inherited retinal dystrophy. Invest Ophthalmol Vis Sci. 2011;52(12):9053–60.
Downes SM, Payne AM, Kelsell RE, Fitzke FW, Holder GE, Hunt DM, Moore AT, Bird AC. Autosomal dominant cone-rod dystrophy with mutations in the guanylate cyclase 2D gene encoding retinal guanylate cyclase-1. Arch Ophthalmol. 2001;119(11):1667–73.
Small KW, Silva-Garcia R, Udar N, Nguyen EV, Heckenlively JR. New mutation, P575L, in the GUCY2D gene in a family with autosomal dominant progressive cone degeneration. Arch Ophthalmol. 2008;126(3):397–403.
Yang L, Cui H, Yin X, Dou H, Zhao L, Chen N, Zhang J, Zhang H, Li G, Ma Z. Dependable and Efficient Clinical Molecular Diagnosis of Chinese RP Patient with Targeted Exon Sequencing. PLoS One. 2015;10(10):e0140684.
Morimura H, Saindelle-Ribeaudeau F, Berson EL, Dryja TP. Mutations in RGR, encoding a light-sensitive opsin homologue, in patients with retinitis pigmentosa. Nat Genet. 1999;23(4):393–4.
Briggs CE, Rucinski D, Rosenfeld PJ, Hirose T, Berson EL, Dryja TP. Mutations in ABCR (ABCA4) in patients with Stargardt macular degeneration or cone-rod degeneration. Invest Ophthalmol Vis Sci. 2001;42(10):2229–36.
Yatsenko AN, Shroyer NF, Lewis RA, Lupski JR. Late-onset Stargardt disease is associated with missense mutations that map outside known functional regions of ABCR (ABCA4). Hum Genet. 2001;108(4):346–55.
Wiszniewski W, Lewis RA, Stockton DW, Peng J, Mardon G, Chen R, Lupski JR. Potential involvement of more than one locus in trait manifestation for individuals with Leber congenital amaurosis. Hum Genet. 2011;129(3):319–27.
Akahori M, Tsunoda K, Miyake Y, Fukuda Y, Ishiura H, Tsuji S, Usui T, Hatase T, Nakamura M, Ohde H, et al. Dominant mutations in RP1L1 are responsible for occult macular dystrophy. Am J Hum Genet. 2010;87(3):424–9.
Davidson AE, Sergouniotis PI, Mackay DS, Wright GA, Waseem NH, Michaelides M, Holder GE, Robson AG, Moore AT, Plagnol V, et al. RP1L1 variants are associated with a spectrum of inherited retinal diseases including retinitis pigmentosa and occult macular dystrophy. Hum Mutat. 2013;34(3):506–14.
Eudy JD, Weston MD, Yao S, Hoover DM, Rehm HL, Ma-Edmonds M, Yan D, Ahmad I, Cheng JJ, Ayuso C, et al. Mutation of a gene encoding a protein with extracellular matrix motifs in Usher syndrome type IIa. Science. 1998;280(5370):1753–7.
Lenassi E, Vincent A, Li Z, Saihan Z, Coffey AJ, Steele-Stallard HB, Moore AT, Steel KP, Luxon LM, Heon E, et al. A detailed clinical and molecular survey of subjects with nonsyndromic USH2A retinopathy reveals an allelic hierarchy of disease-causing variants. Eur J Hum Genet. 2015;23(10):1318–27.
Pozo MG, Bravo-Gil N, Mendez-Vidal C, Montero-de-Espinosa I, Millan JM, Dopazo J, Borrego S, Antinolo G. Re-evaluation casts doubt on the pathogenicity of homozygous USH2A p.C759F. Am J Med Genet A. 2015;167(7):1597–600.
Audo I, Kohl S, Leroy BP, Munier FL, Guillonneau X, Mohand-Said S, Bujakowska K, Nandrot EF, Lorenz B, Preising M, et al. TRPM1 is mutated in patients with autosomal-recessive complete congenital stationary night blindness. Am J Hum Genet. 2009;85(5):720–9.
Payne AM, Downes SM, Bessant DA, Taylor R, Holder GE, Warren MJ, Bird AC, Bhattacharya SS. A mutation in guanylate cyclase activator 1A (GUCA1A) in an autosomal dominant cone dystrophy pedigree mapping to a new locus on chromosome 6p21.1. Hum Mol Genet. 1998;7(2):273–7.
Sullivan LS, Bowne SJ, Birch DG, Hughbanks-Wheaton D, Heckenlively JR, Lewis RA, Garcia CA, Ruiz RS, Blanton SH, Northrup H, et al. Prevalence of disease-causing mutations in families with autosomal dominant retinitis pigmentosa: a screen of known genes in 200 families. Invest Ophthalmol Vis Sci. 2006;47(7):3052–64.
Poloschek CM, Bach M, Lagreze WA, Glaus E, Lemke JR, Berger W, Neidhardt J. ABCA4 and ROM1: implications for modification of the PRPH2-associated macular dystrophy phenotype. Invest Ophthalmol Vis Sci. 2010;51(8):4253–65.
Chakarova CF, Papaioannou MG, Khanna H, Lopez I, Waseem N, Shah A, Theis T, Friedman J, Maubaret C, Bujakowska K, et al. Mutations in TOPORS cause autosomal dominant retinitis pigmentosa with perivascular retinal pigment epithelium atrophy. Am J Hum Genet. 2007;81(5):1098–103.
Garcia-Garcia G, Aparisi MJ, Jaijo T, Rodrigo R, Leon AM, Avila-Fernandez A, Blanco-Kelly F, Bernal S, Navarro R, Diaz-Llopis M, et al. Mutational screening of the USH2A gene in Spanish USH patients reveals 23 novel pathogenic mutations. Orphanet J Rare Dis. 2011;6:65.
Ramsden SC, O'Grady A, Fletcher T, O'Sullivan J, Hart-Holden N, Barton SJ, Hall G, Moore AT, Webster AR, Black GC. A clinical molecular genetic service for United Kingdom families with choroideraemia. Eur J Med Genet. 2013;56(8):432–8.
Riveiro-Alvarez R, Lopez-Martinez MA, Zernant J, Aguirre-Lamban J, Cantalapiedra D, Avila-Fernandez A, Gimenez A, Lopez-Molina MI, Garcia-Sandoval B, Blanco-Kelly F, et al. Outcome of ABCA4 disease-associated alleles in autosomal recessive retinal dystrophies: retrospective analysis in 420 Spanish families. Ophthalmology. 2013;120(11):2332–7.
Beales PL, Badano JL, Ross AJ, Ansley SJ, Hoskins BE, Kirsten B, Mein CA, Froguel P, Scambler PJ, Lewis RA, et al. Genetic interaction of BBS1 mutations with alleles at other BBS loci can result in non-Mendelian Bardet-Biedl syndrome. Am J Hum Genet. 2003;72(5):1187–99.
Bonilha VL, Rayborn ME, Bell BA, Marino MJ, Pauer GJ, Beight CD, Chiang J, Traboulsi EI, Hollyfield JG, Hagstrom SA. Histopathological comparison of eyes from patients with autosomal recessive retinitis pigmentosa caused by novel EYS mutations. Graefes Arch Clin Exp Ophthalmol. 2015;253(2):295–305.
Bowne SJ, Sullivan LS, Blanton SH, Cepko CL, Blackshaw S, Birch DG, Hughbanks-Wheaton D, Heckenlively JR, Daiger SP. Mutations in the inosine monophosphate dehydrogenase 1 gene (IMPDH1) cause the RP10 form of autosomal dominant retinitis pigmentosa. Hum Mol Genet. 2002;11(5):559–68.
Dryja TP, McGee TL, Reichel E, Hahn LB, Cowley GS, Yandell DW, Sandberg MA, Berson EL. A point mutation of the rhodopsin gene in one form of retinitis pigmentosa. Nature. 1990;343(6256):364–6.
Hernan I, Gamundi MJ, Planas E, Borras E, Maseras M, Carballo M. Cellular expression and siRNA-mediated interference of rhodopsin cis-acting splicing mutants associated with autosomal dominant retinitis pigmentosa. Invest Ophthalmol Vis Sci. 2011;52(6):3723–9.
Coppieters F, Casteels I, Meire F, De Jaegere S, Hooghe S, van Regemorter N, Van Esch H, Matuleviciene A, Nunes L, Meersschaut V, et al. Genetic screening of LCA in Belgium: predominance of CEP290 and identification of potential modifier alleles in AHI1 of CEP290-related phenotypes. Hum Mutat. 2010;31(10):E1709–66.
Baux D, Blanchet C, Hamel C, Meunier I, Larrieu L, Faugere V, Vache C, Castorina P, Puech B, Bonneau D, et al. Enrichment of LOVD-USHbases with 152 USH2A genotypes defines an extensive mutational spectrum and highlights missense hotspots. Hum Mutat. 2014;35(10):1179–86.
McGee TL, Seyedahmadi BJ, Sweeney MO, Dryja TP, Berson EL. Novel mutations in the long isoform of the USH2A gene in patients with Usher syndrome type II or non-syndromic retinitis pigmentosa. J Med Genet. 2010;47(7):499–506.
Simeonov DR, Wang X, Wang C, Sergeev Y, Dolinska M, Bower M, Fischer R, Winer D, Dubrovsky G, Balog JZ, et al. DNA variations in oculocutaneous albinism: an updated mutation list and current outstanding issues in molecular diagnostics. Hum Mutat. 2013;34(6):827–35.
Giebel LB, Tripathi RK, Strunk KM, Hanifin JM, Jackson CE, King RA, Spritz RA. Tyrosinase gene mutations associated with type IB ("yellow") oculocutaneous albinism. Am J Hum Genet. 1991;48(6):1159–67.
Cremers FP, van de Pol DJ, van Driel M, den Hollander AI, van Haren FJ, Knoers NV, Tijmes N, Bergen AA, Rohrschneider K, Blankenagel A, et al. Autosomal recessive retinitis pigmentosa and cone-rod dystrophy caused by splice site mutations in the Stargardt's disease gene ABCR. Hum Mol Genet. 1998;7(3):355–62.
Allikmets R, Shroyer NF, Singh N, Seddon JM, Lewis RA, Bernstein PS, Peiffer A, Zabriskie NA, Li Y, Hutchinson A, et al. Mutation of the Stargardt disease gene (ABCR) in age-related macular degeneration. Science. 1997;277(5333):1805–7.
Fernandez-Martinez L, Letteboer S, Mardin CY, Weisschuh N, Gramer E, Weber BH, Rautenstrauss B, Ferreira PA, Kruse FE, Reis A, et al. Evidence for RPGRIP1 gene as risk factor for primary open angle glaucoma. Eur J Hum Genet. 2011;19(4):445–51.
Guo J, Cooper LF. Influence of an LRP5 cytoplasmic SNP on Wnt signaling and osteoblastic differentiation. Bone. 2007;40(1):57–67.
Tuson M, Marfany G, Gonzalez-Duarte R. Mutation of CERKL, a novel human ceramide kinase gene, causes autosomal recessive retinitis pigmentosa (RP26). Am J Hum Genet. 2004;74(1):128–38.
Chaki M, Hoefele J, Allen SJ, Ramaswami G, Janssen S, Bergmann C, Heckenlively JR, Otto EA, Hildebrandt F. Genotype-phenotype correlation in 440 patients with NPHP-related ciliopathies. Kidney Int. 2011;80(11):1239–45.
Lin YY, Wei AH, He X, Zhou ZY, Lian S, Zhu W. A comprehensive study of oculocutaneous albinism type 1 reveals three previously unidentified alleles on the TYR gene. Eur J Dermatol. 2014;24(2):168–73.
Jaakson K, Zernant J, Kulm M, Hutchinson A, Tonisson N, Glavac D, Ravnik-Glavac M, Hawlina M, Meltzer MR, Caruso RC, et al. Genotyping microarray (gene chip) for the ABCR (ABCA4) gene. Hum Mutat. 2003;22(5):395–403.
Hutton SM, Spritz RA. Comprehensive analysis of oculocutaneous albinism among non-Hispanic caucasians shows that OCA1 is the most prevalent OCA type. J Invest Dermatol. 2008;128(10):2442–50.
Sundin OH, Yang JM, Li Y, Zhu D, Hurd JN, Mitchell TN, Silva ED, Maumenee IH. Genetic basis of total colourblindness among the Pingelapese islanders. Nat Genet. 2000;25(3):289–93.
Lewis RA, Shroyer NF, Singh N, Allikmets R, Hutchinson A, Li Y, Lupski JR, Leppert M, Dean M. Genotype/Phenotype analysis of a photoreceptor-specific ATP-binding cassette transporter gene, ABCR, in Stargardt disease. Am J Hum Genet. 1999;64(2):422–34.
The authors would like to thank Dr. Monique Leys for her clinical advice, Dr. Anand Swaroop, Matthew Brooks, Harsha K Rajasimha and Delphine Blain for their scientific advice, technical assistance and support during the course of this study. Dr. J. Fielding Hejtmancik for scientific advice for this manuscript. The authors would like to thank the Exome Aggregation Consortium and the groups that provided exome variant data for comparison.
This research was supported by the US National Institutes of Health (NIH) / National Eye Institute (NEI) intramural and clinical programs.
Availability of data and materials
The data in the study would be available from the corresponding author upon request for any research use.
Ethics approval and consent to participate
This study was reviewed and approved by the Combined Neuroscience Institutional Review Board of the National Institutes of Health and written informed consent was obtained from each participant as adhering to tenets of the Declarations of Helsinki. In case of minors, the consent for participation in the study was obtained from either of the parents.
Consent for publication
Not applicable as the study does not identify the subjects and no personal information identifying the subject was included in the manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Summary of genes and their clinical correlation. A panel of 184 RD genes responsible for retinal related diseases was provided. (PDF 45 kb)
Summary of Retinal Dystrophy primer library design. Statisitc analysis of designing for the primer sets covering entire coding and intronic flanking regions. (PDF 8 kb)
Summary of candidate variants in this cohort. A complete list of candidate variants that found in this patient group before clinical correlation was established. (PDF 301 kb)
Pathogenic variants in OCA related genes. Carrying status of pathogenic variants in the OCA genes in some of the RDS patients. (PDF 275 kb)
About this article
Cite this article
Wang, X., Zein, W.M., D’Souza, L. et al. Applying next generation sequencing with microdroplet PCR to determine the disease-causing mutations in retinal dystrophies. BMC Ophthalmol 17, 157 (2017). https://doi.org/10.1186/s12886-017-0549-5
- Microdroplet PCR
- Mutation screening