The study protocol was approved by the Institutional Review Board and adhered to the tenets of the Declaration of Helsinki. Probands with a family history suggestive of recessive RP were included in the study and available family members were enrolled. All subjects, both affected and unaffected, were clinically evaluated and informed consent was obtained. A total of 34 families with 2 or more affected individuals were recruited; 24 families were consanguineous and 10 were nonconsanguineous. Essential diagnostic criteria for inclusion included bilateral, diffuse, and widespread retinal pigment epithelial degeneration, arterial narrowing, commensurate visual field loss, and reduced amplitudes on electroretinogram (ERG) reduced to less than 25% of the maximum retinal response in normal individuals (normal amplitude of b-wave >350 μV and a-wave >110 μV) with evidence of rod and cone involvement. Other clinical signs that were supportive but not essential for diagnosis of RP included pigment migration including bone–corpuscular pigmentation, vitreous opacities and vitreous pigments, associated retinal pigment epithelium atrophic changes in the macular area, and diffuse disc pallor. Excluded were patients who had unilateral disease, nystagmus, and eye-poking behavior in childhood; exudative retinal detachment; retinal vasculitis; chorioretinitis; or any other secondary cause of pigmentary retinal changes. Clinical features of patients were reviewed and confirmed independently by two ophthalmologists.
Blood samples were collected from the affected and unaffected members of the families by venipuncture. The DNA was extracted from this blood by the phenol-chloroform method.
Twenty-three candidate genes were selected for screening, which included 14 known genes for ARRP: phosphodiesterase 6A (PDE6A), phosphodiesterase 6B (PDE6B), rhodopsin (RHO), cyclic nucleotide gated channel alpha 1 (CNGA1), cyclic nucleotide gated channel beta 1 (CNGB1), crumbs homolog 1 (CRB1), retinitis pigmentosa 1 (RP1), neural retina leucine zipper (NRL), ATP-binding cassette subfamily A member 4 (ABCA4), cellular retinaldehyde binding protein 1(RLBP1), retinal pigment epithelium protein 65 kDa (RPE65), retinal G-protein coupled receptor (RGR), tubby-like protein 1 (TULP1), and prominin 1 (PROM1); 7 genes for related disorders such as Leber congenital amaurosis (LCA), cone-rod dystrophy (CRD), and dominant/digenic RP: guanylate cyclase 2D, membrane (retina-specific) (GUCY2D), guanylate cyclase activator 1A (GUCA1A), rod outer segment membrane protein 1 (ROM1), retinal degeneration slow (RDS), cone-rod homeobox (CRX), aryl hydrocarbon receptor interacting protein-like 1 (AIPL1), and RPGR interacting protein 1 (RPGRIP1); and 2 genes that are candidates for retinal dystrophy but have not yet been shown to have mutations in humans: phosphodiesterase 6G (PDE6G) and cellular retinol binding protein 1 (RBP1). Each locus was screened for homozygosity by genotyping 2 or more microsatellite markers (total of 57 markers). Microsatellite markers were selected based on reported high heterozygosity (0.7 or more) and were generally located within an interval of ∼5 to 10 Mb of the candidate gene. Information on the primers for amplification of microsatellite markers, marker heterozygosity, and location was obtained from the UniSTS (http://www.ncbi.nlm.nih.gov/sites/entrez?db=unists) Human Genome Database and NCBI Mapview, (http://www.ncbi.nlm.nih.gov/mapview/) databases (both provided in the public domain by the National Center for Biotechnology Information, Bethesda, MD). The detection of homozygosity at a given locus shared only by affected members but not by unaffected family members was investigated further by typing additional markers at the locus for confirming homozygosity and subsequent screening of the relevant gene for mutations. Uninformative loci in which affected as well as one or more unaffected members were homozygous, were genotyped with additional markers. Genotyping was performed for 76 affected and 88 unaffected individuals from 34 families. Genotyping was performed in multiplex PCR reactions followed by electrophoresis (model 310 genetic analyzer; ABI). Alleles were determined by using genotyping software (GeneScan; Applied Biosystems Inc, Foster City, CA). Screening of coding regions of genes was performed by PCR amplification of exons and adjacent intronic regions, followed by direct automated sequencing. The sequence changes observed were checked for cosegregation in the family and for presence or absence in at least 100 healthy control individuals by RFLP or direct sequencing. For RFLP, restriction enzyme-digested products were resolved on 8% or 10% acrylamide gels and visualized after ethidium bromide staining.
Multiple sequence alignment of protein sequences was performed using ClustalW2 (http://www.ebi.ac.uk/Tools/clustalw2/index.html; European Bioinformatics Institute, Cambridge, UK). Sorting intolerant from tolerant (SIFT) analysis (http://blocks.fhcrc.org/sift/ provided in the public domain by the Fred Hutchinson Cancer Research Center, Seattle, WA) was used to predict the potential impact of a missense substitution on protein function. A SIFT score below the cutoff of 0.05 for a given substitution is classified as not tolerated while those with scores higher than this value are considered tolerated.