Received 19 April 2010 | Accepted 3 December 2010 | Published 8 December 2010
Namerah Sabir,1 S. Amer Riazuddin,1,2 Tariq Butt,1 Farheena Iqbal,1 Idrees A. Nasir,1 Ahmad U. Zafar,1 Zaheeruddin A. Qazi,3 Nadeem H. Butt,4 Shaheen N. Khan,1 Tayyab Husnain,1 J. Fielding Hejtmancik,5 Sheikh Riazuddin1,4
Namerah Sabir, S. Amer Riazuddin, J. Fielding Hejtmancik, and Sheikh Riazuddin contributed equally to this work
1National Centre of Excellence in Molecular Biology, University of the Punjab, Lahore, Pakistan; 2The Wilmer Eye Institute, Johns Hopkins University School of Medicine, Baltimore, MD; 3Layton Rahmatulla Benevolent Trust Hospital Lahore, Pakistan; 4Allama Iqbal Medical College, University of Health Sciences, Lahore, Pakistan; 5Ophthalmic Genetics and Visual Function Branch, National Eye Institute, National Institutes of Health, Bethesda, MD
Correspondence to: S. Riazuddin, Ph.D., National Centre of Excellence in Molecular Biology, 87 West Canal Bank Road, Lahore 53700, Pakistan; Phone: 92-42-542-1235, FAX: 92-42-542-1316; email: firstname.lastname@example.org
Purpose: To localize the disease interval for autosomal recessive congenital cataracts in a consanguineous Pakistani family.
Methods: All affected individuals underwent detailed ophthalmologic examination. Blood samples were collected and genomic DNA was extracted. A genome-wide scan was performed with fluorescently-labeled microsatellite markers on genomic DNA from affected and unaffected family members and logarithm of odds (LOD) scores were calculated.
Results: Clinical records and ophthalmological examinations suggested that affected individuals have bilateral congenital cataracts. Genome-wide linkage analysis localized the critical interval to chromosome 3q with a maximum LOD score of 3.87 at θ=0; with marker D3S3609. Haplotype analyses refined the critical interval to a 23.39 cM (18.01 Mb) interval on chromosome 3q, flanked by D3S1614 proximally and D3S1262, distally.
Conclusions: Here, we report a new locus for autosomal recessive congenital cataract localized to chromosome 3q in a consanguineous Pakistani family.
Congenital cataracts are one of the major causes of vision loss in children worldwide [1,2]. They can occur in an isolated fashion or as one component of a syndrome affecting multiple tissues . Congenital cataracts can lead to permanent blindness by interfering with the sharp focus of light on the retina, especially during the early developmental periods. Congenital cataracts vary markedly in severity and morphology, affecting the nuclear, cortical, polar, or subcapsular parts of the lens, or in severe cases the entire lens .
To-date fourteen loci have been associated with autosomal recessive congenital cataracts [5-18]. Of these loci, mutations in nine genes, Eph-receptor type-A2 (EPHA2), connexin50 (GJA8), glucosaminyl (N-acetyl) transferase 2 (GCNT2), heat-shock transcription factor 4 (HSF4), lens intrinsic membrane protein (LIM2), beaded filament structural protein 1 (BFSP1), alphaA-crystallin (CRYαA), betaB1-crystallin (CRYβB1), and betaB3-crystallin (CRYβB3) have been found [5,7,9,12,14-18].
Here, we report a consanguineous Pakistani family with congenital cataracts. Linkage analysis localized the critical interval to chromosome 3q with a maximum LOD score of 3.87 and haplotype analyses refined the critical interval to a 23.39 cM (18.01 Mb). This is the first report localizing autosomal recessive congenital cataracts to chromosome 3q.
A total of 125 consanguineous Pakistani families with nonsyndromic cataract were recruited to participate in a collaborative study between the National Centre of Excellence in Molecular Biology, Lahore, Pakistan, and the National Eye Institute, Bethesda, MD, to identify novel loci associated congenital cataracts. Institutional Review Board (IRB) approval was obtained from both Institutes. The participating subjects gave informed consent consistent with the tenets of the Declaration of Helsinki. A detailed medical history was obtained by interviewing family members. Ophthalmic examinations were conducted with slit-lamp microscopy. Approximately 10 ml of blood samples were drawn from affected and unaffected members of the family and stored in 50 ml Sterilin® falcon tubes (BD Biosciences, San Jose, CA) containing 400 μl of 0.5 M EDTA. Blood samples were kept at −20 °C for long- term storage.
DNA was extracted by a nonorganic method as described previously [5,10]. Briefly, aliquots of 10 ml blood samples were mixed with 35 ml of TE buffer (10 mM Tris-HCl, 2 mM EDTA, pH 8.0) and the TE-blood mixture was centrifuged at 1,800× g for 20 min. The supernatant was discarded and the pellet was re-suspended in 35 ml of TE buffer and centrifuged at 1,800× g for 20 min. The TE washing was repeated for 2–3 times and the washed pellet was re-suspended in 2 ml of TE. A total of 6.25 ml of protein digestion cocktail (50 μl [10 mg/ml] of proteinase K, 6 ml TNE buffer [10 mM Tris HCl, 2 mM EDTA, 400 mM NaCl], and 200 μl of 10% sodium dodecyl sulfate) was added to the re-suspended pellets and incubated overnight in a shaker (250 rpm) at 37 °C. The digested proteins were precipitated by adding 1 ml of 5 M NaCl, followed by vigorous shaking and chilling on ice for 15 min. The precipitated proteins were pelleted by centrifugation at 1,800× g for 20 min and removed. The supernatant was mixed with equal volumes of phenol/chloroform/isoamyl alcohol (25:24:1) and the aqueous layer containing the genomic DNA was carefully collected. The DNA was precipitated with isopropanol and pelleted by centrifugation at 2,400× g for 15 min. The DNA pellets were washed with 70% ethanol and dissolved in TE buffer. The DNA concentration was determined with a SmartSpec plus Bio-Rad Spectrophotometer (Bio-Rad, Hercules, CA).
A genome-wide scan was performed with 382 highly polymorphic fluorescent markers from the ABI PRISM Linkage Mapping Set MD-10 (Applied Biosystems, Foster City, CA) having an average spacing of 10 cM. Multiplex polymerase chain reaction (PCR) was completed in a GeneAmp PCR System 9700 thermocycler (Applied Biosystems). Briefly, each reaction was performed in a 5 μl mixture containing 40 ng genomic DNA, various combinations of 10 mM dye-labeled primer pairs, 0.5 ml 10× GeneAmp PCR Buffer (Applied Biosystems), 1 mM dNTP mix, 2.5 mM MgCl2, and 0.2 U Taq DNA polymerase (Applied Biosystems). Initial denaturation was performed for 5 min at 95 °C, followed by 10 cycles of 15 s at 94 °C, 15 s at 55 °C, and 30 s at 72 °C and then 20 cycles of 15 s at 89 °C, 15 s at 55 °C, and 30 s at 72 °C. The final extension was performed for 10 min at 72 °C. PCR products from each DNA sample were pooled and mixed with a loading cocktail containing HD-400 size standards (Applied Biosystems). The resulting PCR products were separated in an ABI 3100 DNA Analyzer (Applied Biosystems) and genotypes were assigned with GeneMapper software (Applied Biosystems).
Two-point linkage analyses were performed using the FASTLINK version of MLINK from the LINKAGE Program Package, whereas the Maximum LOD scores were calculated with ILINK from the LINKAGE Program Package [19,20]. Autosomal recessive cataract was analyzed as a fully penetrant trait with an affected allele frequency of 0.001. The marker order and distances between the markers were obtained from the Marshfield database and the National Center for Biotechnology Information (NCBI) chromosome 3 sequence maps. For the initial genome scan, equal allele frequencies were assumed, while for fine mapping allele frequencies were estimated from 96 unrelated and unaffected individuals from the Punjab province of Pakistan.
A large consanguineous family, PKCC144 consisting of six affected and three unaffected individuals was recruited from the Punjab province of Pakistan (Figure 1). A detailed medical and family history was obtained from all affected and unaffected members of family. The clinical records available to us were suggestive of bi-lateral congenital cataracts with no other ocular or systemic abnormalities present in the family (Table 1). As all the affected individuals have had cataract surgery prior to the enrollment, therefore no photographic evidence of cataract was available to classify the cataracts.
Initially, linkage to all reported loci associated with autosomal recessive cataract loci were excluded by haplotype analyses with closely spaced fluorescently-labeled microsatellite markers (data not shown). Subsequently, a genome-wide linkage scan was completed and a maximum two-point LOD score of 3.87 at θ=0 was obtained with D3S3609 during the genome-wide scan (Table 2). Additional STR markers were selected to further analyze the critical interval, which further provided evidence for linkage to chromosome 3q with LOD scores of 3.42, 2.02, 2.33, and 2.04, at θ=0 with markers D3S1565, D3S3699, D3S2314, and D3S1618, respectively (Table 2).
Visual inspection of the haplotypes supported the results of linkage analysis (Figure 1). There is a recombination in affected individual 14 and D3S1614 that defines the proximal boundary (Figure 1). Similarly, recombination events in affected individuals 13 and 14 at D3S1262 define the distal boundary (Figure 1). This places the pathogenic locus in a 23.39 cM (18.01 Mb) interval. Additionally, a lack of homozygosity in affected individuals 10, 11, 12, 13, 14, and 15 at D3S3592 further suggests that the pathogenic mutations lies proximal to D3S3592, in a 20.93 cM (16.23 Mb) interval flanked by D3S1614 proximally and D3S3592 distally.
Here, we report a new locus for autosomal recessive congenital cataracts localized to chromosome 3q, in a consanguineous Pakistani family. Genome-wide linkage analyses localized the critical interval to chromosome 3q with significant LOD scores and these results were further supported by haplotype analyses, which places the critical interval 23.39 cM (18.01 Mb) flanked by markers, D3S1614 proximally and D3S1262 distally. Additionally, lack of any significant Lod scores obtained during the genome-wide scan other than with marker D3S3609, further supports linkage to chromosome 3q. This is the first report associating chromosome 3q with autosomal recessive congenital cataracts.
The critical interval is a gene rich region that harbors more than 100 genes according to the UCSC database. Subsequent to prioritizing the genes present in the critical interval based on their known or proposed function and expression especially in the lens, we are sequencing coding exons of these candidate genes in one affected and one unaffected individual. Identification of the gene responsible for recessive congenital cataract is expected not only to provide insight and clues to the allied clinical findings but it will also increase our understanding of molecular events involved in the lens biology.
The authors are thankful to all the family members for their participation in this study. This study was supported, in part by Higher Education Commission and Ministry of Science and Technology Islamabad, Pakistan, and the National Academy of Sciences, and the U. S. Department of State, Washington DC USA (Note: All findings and conclusions are those of the authors and do not necessarily reflect the views of the National Academy of Sciences and the U.S. Department of State).