Molecular Vision 2016; 22:610-625
Received 16 October 2015 | Accepted 08 June 2016 | Published 10 June 2016
Firoz Kabir,1 Inayat Ullah,2 Shahbaz Ali,2 Alexander D.H. Gottsch,1 Muhammad Asif Naeem,2 Muhammad Zaman Assir,3,4 Shaheen N. Khan,2 Javed Akram,3,4 Sheikh Riazuddin,2, 3, 4 Radha Ayyagari,5 J. Fielding Hejtmancik,6 S. Amer Riazuddin1
The first two and last two authors contributed equally to this work.
1The Wilmer Eye Institute, Johns Hopkins University School of Medicine, Baltimore, MD; 2National Centre of Excellence in Molecular Biology, University of the Punjab, Lahore, Pakistan; 3Allama Iqbal Medical College, University of Health Sciences, Lahore, Pakistan; 4National Centre for Genetic Diseases, Shaheed Zulfiqar Ali Bhutto Medical University, Islamabad, Pakistan; 5Shiley Eye Institute, University of California, San Diego, CA; 6Ophthalmic Genetics and Visual Function Branch, National Eye Institute, National Institutes of Health, Bethesda, MD
Correspondence to: S. Amer Riazuddin, Johns Hopkins University School of Medicine, 600 N. Wolfe Street; Maumenee 840, Baltimore, MD 21287; Phone: (410) 955-3656; FAX: (410) 955-3656; email: email@example.com
Purpose: This study was undertaken to identify causal mutations responsible for autosomal recessive retinitis pigmentosa (arRP) in consanguineous families.
Methods: Large consanguineous families were ascertained from the Punjab province of Pakistan. An ophthalmic examination consisting of a fundus evaluation and electroretinography (ERG) was completed, and small aliquots of blood were collected from all participating individuals. Genomic DNA was extracted from white blood cells, and a genome-wide linkage or a locus-specific exclusion analysis was completed with polymorphic short tandem repeats (STRs). Two-point logarithm of odds (LOD) scores were calculated, and all coding exons and exon–intron boundaries of RP1 were sequenced to identify the causal mutation.
Results: The ophthalmic examination showed that affected individuals in all families manifest cardinal symptoms of RP. Genome-wide scans localized the disease phenotype to chromosome 8q, a region harboring RP1, a gene previously implicated in the pathogenesis of RP. Sanger sequencing identified a homozygous single base deletion in exon 4: c.3697delT (p.S1233Pfs22*), a single base substitution in intron 3: c.787+1G>A (p.I263Nfs8*), a 2 bp duplication in exon 2: c.551_552dupTA (p.Q185Yfs4*) and an 11,117 bp deletion that removes all three coding exons of RP1. These variations segregated with the disease phenotype within the respective families and were not present in ethnically matched control samples.
Conclusions: These results strongly suggest that these mutations in RP1 are responsible for the retinal phenotype in affected individuals of all four consanguineous families.
Retinitis pigmentosa (RP) is the most common inherited retinal dystrophy, affecting approximately 1 in 5,000 individuals worldwide [1,2]. RP primarily affects the rod photoreceptors, while the cone cells are compromised as the disease progresses . Affected individuals exhibit night blindness in the initial stages of the disease followed by a progressive reduction in the visual field . Ocular findings include atrophic changes in the photoreceptors and the RPE followed by the appearance of melanin-containing structures in the retinal vascular layer . The fundus changes include a pale optic nerve, attenuation of the retinal vessels, and bone spicule-like pigmentation in the mid-peripheral retina . Electroretinography (ERG) recordings show severely diminished or completely extinguished rod response while the cone response is somewhat normal in early stages but is undetectable as the disease progresses .
RP is a genetically heterogeneous disorder that manifests as an autosomal dominant, autosomal recessive, and X-linked trait. To date, 73 genes have been implicated in the pathogenesis of RP. Of these genes, 27 have been associated with autosomal dominant RP (adRP) [4-30] while mutations in 50 genes have been identified in patients with autosomal recessive RP (arRP; RetNet) [31-77]. Interestingly, mutations in RHO (Gene ID: 6010; OMIM: 180380), RP1 (Gene ID: 6101; OMIM: 603937), NRL (Gene ID: 4901; OMIM: 162080), RPE65 (Gene ID: 6121; OMIM: 180069), BEST1 (Gene ID: 7439; OMIM: 607854), NR2E3 (Gene ID: 10002; OMIM: 604485), IMPDH1 (Gene ID: 3614; OMIM: 146690) have been identified in familial cases of both adRP and arRP. Likewise, causal mutations in OFD1 (Gene ID: 8481; OMIM: 300170), RP2 (Gene ID: 6102; OMIM: 300757), and RPGR (Gene ID: 6103; OMIM: 312610) have been identified in RP cases with an X-linked inheritance pattern [78-80].
RP1 was localized to chromosome 8q and consists of four exons that encode for a 2,156 amino acid protein . Pierce and colleagues first identified mutations in RP1 responsible for adRP, and subsequently, they estimated that the nonsense mutation at codon 677 (p.R677*) is present in approximately 3% of the dominant RP cases in North America . The RP1 protein localizes to the connecting cilia of the rod and cone cells in the ocular retina and is required for correct stacking of the outer segment disc [81,82].
Here, we report four consanguineous familial cases with multiple members who manifest cardinal symptoms of RP. Genome-wide linkage analyses localized the disease phenotype to chromosome 8q, harboring RP1, while bidirectional Sanger sequencing identified causal mutations in RP1 that segregated with the disease phenotype in their respective families and were absent in the ethnically matched controls and the genome-variant databases.
More than 300 consanguineous Pakistani families with non-syndromic retinal dystrophies were recruited to identify new disease loci responsible for inherited visual diseases. The institutional review boards (IRBs) of the National Centre of Excellence in Molecular Biology (Lahore, Pakistan), the National Eye Institute (Bethesda, MD) and the Johns Hopkins University (Baltimore, MD) approved the study. All participating family members provided informed written consent that was endorsed by the respective IRBs and is consistent with the tenets of the Declaration of Helsinki.
A detailed clinical and medical history was obtained from the individual families. Funduscopy was performed at Layton Rehmatulla Benevolent Trust (LRBT) Hospital (Lahore, Pakistan). ERG measurements were recorded by using equipment manufactured by LKC (Gaithersburg, MD). Dark-adapted rod responses were determined through incident ﬂash attenuated by −25 dB, whereas rod–cone responses were measured at 0 dB. The 30 Hz flicker responses were recorded at 0 dB to a background illumination of 17 to 34 cd/m2.
All participating members voluntarily provided an approximately 10 ml blood sample that was stored in 50 ml Sterilin® Falcon (Sarstedt, Inc. Newton, NC) tubes containing 400 μl of 0.5 M EDTA. Blood samples were stored at −20 °C for long-term storage.
Genomic DNA was extracted from white blood cells using a non-organic modified procedure as described previously . The concentration of the extracted genomic DNA was estimated with a SmartSpec plus Bio-Rad Spectrophotometer (Bio-Rad, Hercules, CA).
The Applied Biosystems MD-10 linkage mapping panels (Applied Biosystems, Foster City, CA) were used to complete a genome-wide scan for family PKRP117. PCR was completed in a 5 μl reaction volume containing 40 ng of genomic DNA, various combinations of 10 μM fluorescently labeled primer pairs, 10X PCR buffer (100 mM Tris HCl pH 8.4, 400 mM NaCl, 15 mM MgCl2, 2.5 mM Spermidine), 2 mM deoxynucleotide triphosphate (dNTP) mix, and 0.2 U OneTaq DNA polymerase (New England BioLabs Inc, Ipswich, MA). Initial denaturation was performed for 5 minutes (min) at 95 °C, followed by 10 cycles consisting of denaturation at 94 °C for 15 s, annealing at 55 °C for 15 s and elongating at 72 °C for 30 s, and then 20 cycles consisting of denaturation at 89 °C for 15 s, annealing at 55 °C for 15 s, and elongation at 72 °C for 30 s. The final extension was performed for 10 min at 72 °C, followed by a final hold at 16 °C. PCR products were mixed with a loading cocktail containing HD-400 size standards (Applied Biosystems) and resolved in an Applied Biosystems 3100 DNA Analyzer. Genotypes were assigned using the Gene Mapper software from Applied Biosystems. Exclusion analyses were completed for PKRP262, PKRP344, and PKRP358 using closely spaced short tandem repeat (STR) markers.
Linkage analysis was performed with alleles of PKRP117 obtained through the genome-wide scan and the alleles of PKRP262, PKRP344, and PKRP358 obtained through exclusion analysis using the FASTLINK version of MLINK from the LINKAGE Program Package [84,85]. Maximum LOD scores were calculated using ILINK from the LINKAGE Program Package. Autosomal recessive retinitis pigmentosa was investigated as a fully penetrant disorder that has an affected allele frequency of 0.001.
The sequences of primer pairs used to amplify RP1 exons were designed using the Primer3 software. The sequences of the primer pairs used for sequencing are shown in Appendix 1. PCR reactions were completed in 10 μl volume containing 20 ng of genomic DNA. PCR amplification consisted of a denaturation step at 95 °C for 5 min followed by a two-step touchdown procedure. The first step of ten cycles consisted of denaturation at 95 °C for 30 s, followed by a primer set-specific annealing for 30 s (annealing temperature decreases by 1 °C per cycle) and elongation at 72 °C for 45 s. The second step of 30 cycles consisted of denaturation at 95 °C for 30 s followed by annealing (annealing temperature −10 °C) for 30 s and elongation at 72 °C for 45 s, followed by a final elongation at 72 °C for 5 min.
The PCR primers for each exon were used for bidirectional sequencing using the BigDye Terminator Ready reaction mix (Applied Biosystems, Foster City, CA), according to the manufacturer’s instructions. The sequencing products were resolved on an ABI PRISM 3100 DNA analyzer (Applied Biosystems), and results were analyzed with Applied Biosystems SeqScape software.
The degree of evolutionary conservation of the splice donor site (c.787+1G) in other RP1 orthologs was examined using the UCSC Genome Browser (Genome). The effect of the c.787+1G>A mutation on RP1 mRNA splicing was predicted with an online bioinformatics tool, the Human Splicing Finder 2.4.1 (HSF).
In an ongoing effort to investigate the genetic load of retinal dystrophies in the Pakistani population, we recruited a large cohort of familial cases with multiple members in these families manifesting cardinal symptoms of early onset RP. Among these families, PKRP117, PKRP262, PKRP344, and PKRP358 were recruited from the Punjab province of Pakistan. In PKRP117, we were able to enroll a total of 19 family members, including nine affected individuals (Figure 1A). A detailed medical history, including the onset and progression of the ocular disease, was obtained by interviewing family elders, especially the parents of the affected individuals, which revealed that all affected individuals complained of night blindness during the early years of their life. Exclusion analysis with closely spaced STR markers spanning known RP loci suggested linkage to chromosome 8q harboring RP1. Maximum two-point LOD scores of 3.21, 6.95, and 5.95 at θ = 0 were obtained with markers D8S1737, D8S509, and D8S2332, respectively (Table 1).
In parallel, we enrolled three affected and nine unaffected members of PKRP262, four affected and three unaffected members of PKRP344, and four affected and three unaffected members of PKRP358 (Figure 1B-D). Fundus photographs of the affected individuals revealed typical symptoms of RP, including attenuated retinal arteries, waxy, pale optic disc, and bone spicule-like pigment deposits in the lateral and mid-periphery of the retina (Figure 2). Likewise, scotopic ERG recordings measured at −25 dB and photopic responses at 0 dB (30 Hz flicker) were undetectable in the affected individuals, suggestive of compromised rod photoreceptor and cone cells, while unaffected individuals exhibited rod and cone responses in the normal range (Figure 3).
Exclusion analysis localized all three familial cases (PKRP262, PKRP344, and PKRP358) to a region of chromosome 8q harboring RP1. Alleles for markers D8S1737, D8S509, and D8S2332 yielded 2-point LOD scores of 2.54, 2.51, and 2.51 at θ = 0 and 1.42, 2.13, and 2.52 at θ = 0 for families PKRP262 and PKRP358, respectively, while the alleles for markers D8S1737 and D8S509 yielded 2-point LOD scores of 1.62 and 1.15 at θ = 0 for PKRP344 (Table 1). Although these LOD scores are less than a 2-point LOD score of 3.0, which is traditionally considered sufficient for linkage, they are the maximum two-point LOD scores attainable by PKRP262, PKRP344, and PKRP358 and were considered worthy of further investigation because the known RP locus, RP1, was included in the region.
To identify the causal mutation responsible for the RP phenotype in these four families, we sequenced all coding exons and the exon–intron boundaries of RP1. In PKRP117, we identified a 1-bp homozygous deletion in exon 4, c.3697delT, that is predicted to result in a frameshift p.S1233Pfs22* (Figure 4A-C). Likewise, we identified a homozygous variation in intron 3, c.787+1G>A (p.I263Nfs8*), in PKRP262 that affects the conserved splice donor site (Figure 4D-F) and a 2 bp duplication in exon 2, c.551_552dupTA (p.Q185Yfs4*), in PKRP344 (Figure 4G-I).
In PKRP358, PCR of all three coding exons of RP1 did not yield any amplification products for the affected individuals while the genomic DNA of the unaffected individuals produced amplified products of the appropriate size. One plausible explanation is that the affected individuals of PKRP358 harbor a homozygous deletion that removes the coding region of RP1. We designed six primer pairs between exons 1 and 2 of RP1 and another six primer pairs downstream of RP1 using Primer3 software. The sequences are available upon request. Briefly, each primer pair was PCR-amplified in a 10 μl reaction volume containing 20 ng of genomic DNA, 1 μl of 10X PCR buffer, 2 mM dNTP mix, 1 μl of 10 μM forward and reverse primer, 500 mM Betaine, 700 mM Dimethyl sulfoxide (DMSO), and 0.2 U OneTaq DNA polymerase. The initial denaturation step was at 95 °C for 5 min followed by a two-step touchdown procedure. The first step of 10 cycles consisted of denaturation at 95 °C for 30 s, followed by a 66 °C annealing for 30 s (annealing temperature decreased by 1 °C after every cycle) and elongation at 72 °C for 45 s. The second step of 25 cycles consisted of denaturation at 95 °C for 30 s followed by 56 °C annealing for 30 s and elongation at 72 °C for 45 s, followed by a final elongation at 72 °C for 5 min. The amplification pattern of these 12 primer pairs helped us identify a 11,117 bp deletion (chr8:55,531,690-55, 542, 807) that removes all three coding exons of RP1 (Figure 5A-C).
All four mutations segregated with the disease phenotype in their respective families. All affected individuals were homozygous for the mutant allele while unaffected individuals were either heterozygous carriers or homozygous for the wild-type allele. These mutations were absent in ethnically matched control chromosomes and were not found in the 1000 Genomes, the NHLBI Exome Sequencing Project, and the dbSNP databases.
The mutation identified in PKRP117, c.3697delT (p.S1233Pfs22*), is expected to produce a truncated protein lacking 903 amino acids of the C-terminal, while the transcript harboring the mutation identified in PKRP344, c.551_552dupTA (p.Q185Yfs4*), is expected to degrade through nonsense-mediated decay. Moreover, affected individuals in PKRP358 are believed to have no expression of RP1 due to the large deletion that removes all coding exons of RP1. Therefore, the respective mutations in PKRP117, PKRP344, and PKRP358 are likely to have caused the RP phenotype in these families. However, we sought additional evidence to strengthen the candidacy of the splice donor variation identified in PKRP262.
First, we examined the evolutionary conservation of c.787+1G and found that the +1G of the splice donor site is completely conserved in RP1 orthologs in general and mammals in particular (Appendix 2). Second, we evaluated the effect of the c.787+1G>A variation on RP1 mRNA splicing using Human Splice Finder 2 (HSF2). The HSF2 generated consensus values of 80.96 and 48.29 for the wild-type (c.787+1G) and mutant (c.787+1A) nucleotides, respectively. The predicted consensus value deviation of −32.67 for c.787+1G>A suggests that the loss of the wild-type splice site will result in the retention of intron 3 of RP1 (Figure 5D), resulting in a frame shift and eventually a premature stop codon (p.I263Nfs8*).
We previously reported homozygous mutations responsible for arRP in three consanguineous familial cases that implicated RP1 in the pathogenesis of arRP . Here, we report four consanguineous families recruited from the Punjab province of Pakistan with multiple members who manifest cardinal symptoms of RP. Linkage analyses with closely spaced STR markers localized the linkage interval in all four families to chromosome 8q12.1, harboring RP1, while Sanger sequencing of RP1 identified mutations that segregate with the disease phenotype in their respective families and are predicted to produce truncated RP1 proteins.
A total of 55 mutations have been identified in RP1 associated with RP in multiple ethnic populations [7,48,86-112]. 37 result in adRP while 17 mutations are responsible for arRP and one mutation responsible for sporadic case. The majority of the mutations identified in RP1 are either nonsense codons or lead to premature termination of RP1. Nonsense mutations in mammalian genes generally lead to unstable mRNAs that are degraded by nonsense-mediated decay . However, nonsense-mediated decay does not occur if the mutation is present in the last exon, and therefore, most of the reported mutations are expected to produce stable transcripts translated into truncated proteins.
The precise mechanism that dictates whether a single allele would be sufficient or whether a homozygous mutation would be required for manifestation of the disease phenotype is not yet completely understood and probably varies among different examples. Nonetheless, it is conceivable that variations resulting in a mutant protein with a deleterious effect would manifest as an autosomal dominant trait while functionally null alleles would manifest as an autosomal recessive disease. We previously speculated that disruption of RP1 within or immediately after the bifocal gene product (BIF) domain might result in a protein with a deleterious effect whereas truncation of RP1 before the BIF motif or within the terminal portion of the protein would result in a loss of RP1 function .
More than 50 causal mutations have been identified in patients with adRP and arRP within the past decade, including our earlier observation (Table 2). As shown in Figure 6, pathogenic mutations in the heterozygous state in RP1 responsible for adRP seem to reside between amino acid residues 617–1551. In contrast, the homozygous mutations responsible for arRP cluster in two regions: first, between amino acids 193–736 and, second, between amino acids 1243–1890 (Figure 6). Taken together, these data support our earlier speculation : Truncation mutations especially in the central region (800–1,200) of RP1 that is structurally not well defined, produce mutant proteins with a deleterious effect while mutations close to the N- and C-terminals of RP1 (>600 or <1,600) produce loss of function of the mutant proteins.
Moreover, the regions that include amino acid residues 600–750 and 1,250–1,550 include mutations that have been associated with both adRP and arRP, although no single mutation has been associated with both adRP and arRP. Thus, it is tempting to speculate that “zones” exist in RP1 in which the nature of the mutation instead of its location in the polypeptide dictates the inheritance pattern. Additional functional investigations are required to define these zones and elucidate mechanistic details of mutations leading to the particular genetic trait.
Appendix 1. The primer sequences for the amplification of RP1.
Appendix 2. Sequence alignment of the exon-intron junction illustrating conservations of splice-donor site of intron 3 of RP1. The c.787+1G shown in red are fully conserved in RP1 orthologs.
We are thankful to all family members for their participation in this study. This study was supported in part by Higher Education Commission, Islamabad, Pakistan, and by the National Eye Institute Grant R01EY021237 (RA and SAR).