Molecular Vision 2011; 17:693-707 <http://www.molvis.org/molvis/v17/a79>
Received 10 October 2010 | Accepted 7 March 2011 | Published 11 March 2011

Mutation screening and genotype phenotype correlation of α-crystallin, γ-crystallin and GJA8 gene in congenital cataract

Manoj Kumar,1 Tushar Agarwal,2 Sudarshan Khokhar,2 Manoj Kumar,3 Punit Kaur,3 Tara Sankar Roy,4 Rima Dada1

1Laboratory for Molecular Reproduction and Genetics, Department of Anatomy, All India Institute of Medical Sciences, New Delhi, India; 2Dr. Rajendra Prasad Centre for Ophthalmic Sciences, All India Institute of Medical Sciences, New Delhi, India; 3Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India; 4Department of Anatomy, All India Institute of Medical Sciences, New Delhi, India

Correspondence to: Dr. Rima Dada, Associate Professor, Laboratory for Molecular Reproduction and Genetics, Department of Anatomy, All India Institute of Medical Sciences, New Delhi-110029, India; Phone: +91 11 26543517; FAX: +91 11 26588663; email: rima.dadarima@gmail.com or rima_dada@rediffmail.com

Abstract

Purpose: To screen α-crystallin (CRYAB), γ-crystallin (CRYGC and CRYGD), and Connexin 50 (Cx-50 or GJA8) genes in congenital cataract patients and controls.

Methods: Thirty clinically diagnosed congenital cataract cases below 3 years of age from northern India, presenting at Dr. R. P. Centre for Ophthalmic Sciences (AIIMS, New Delhi, India) were enrolled in this study. Genomic DNA was extracted from peripheral blood, all coding and exon/intron regions were amplified using PCR and direct sequencing was performed to detect any nucleotide variation. ProtScale and Discovery Studio programs were used for insilico and structural analysis of non-synonymous mutations.

Results: DNA sequencing analysis of CRYAB, CRYGC, CRYGD, and GJA8 showed a total of six variations of which two were novel (CRYGC:p.R48H and GJA8:p.L281C) and four have been previously reported (CRYAB: rs11603779T>G, GJA8: p.L268L, CRYGD: p.R95R, and c.T564C). Both the novel changes, in CRYGC and GJA8 were found in 16.6% of the patients. Previously reported nucleotide alterations (CRYGD:p.R95R and c.T564C) were found in 90% of the patients. Insilico and structural analysis data suggested that two novel non-synonymous mutations altered the stability and solvent accessibility of γC-crystallin and Cx-50 proteins which may lead to lens opacification.

Conclusions: We observed two novel nonsynonymous variations and four reported variations in CRYAB, CRYGC, CRYGD, and GJA8. The p.R48H variation in γC-crystallin may disrupt the normal structure of lens and can cause cataract. Cx50 is responsible for joining the lens cells into a functional syncytium and a mutation (p.L281C) in GJA8 may lead to lens opacification resulting in cataract formation. This study further expands the mutation spectrum of congenital cataract and help understanding how mutant proteins lead to opacification of lens.

Introduction

Cataract is an opacification of the lens resulting from alterations in lens cellular architecture or in lens proteins, or both. Congenital cataract is a clinically and genetically heterogeneous lens disease responsible for a significant proportion of visual impairment and blindness in childhood [1,2]. The prevalence of congenital cataract is estimated to vary from 0.6 to 6 per 10,000 live births with an incidence of 2.2–2.49 per 10,000 live births [3]. It is estimated that globally, 20 million children under the age of 16 years suffer from cataract, and among these, 200,000 (15%) are severely visually impaired or blind [4,5]. Pediatric cataracts are responsible for 7.4% to 15.3% of childhood blindness in developing countries like india [6-8]. The prevalence of blindness in children ranges from approximately 0.3/1000 children in affluent regions to 1.5/1000 in the poorest communities. The population of India in 2001 was estimated to be 1.03 billion, approximately 420 million of whom are children under 16 years of age (40.9%). Overall, there are probably 280,000–320,000 blind children in India [9]. Cataract is responsible for ~12% of childhood blindness in India [6,10]. The occurrence of congenital cataract varies in different parts of India as it is 5.25% in northen India [7], 8.5% in northeast states [6], 7.25% in western India [11] and 11.4% in south India [12].

Currently there are ~60 loci implicated in non-syndromic congenital cataract, among these, over 22 have been associated with mutations in specific genes [13-15] including 10 crystallin genes: αA-crystallin (CRYAA), αB-crystallin (CRYAB), βA1-crystallin (CRYBA1), βA4-crystallin (CRYBA4), βB1-crystallin (CRYBB1), βB2-crystallin (CRYBB2), βB3-crystallin (CRYBB3), γC-crystallin (CRYGC), γD-crystallin (CRYGD), γS-crystallin (CRYGS) [16-21], two cytoskeletal protein genes: beaded filament structural protein 1, filensin (BFSP1), beaded filament structural protein 2, phakinin (BFSP2) [22], four membrane protein genes: gap junction protein (alpha 3, GJA3), gap junction protein (alpha 8, GJA8), major intrinsic protein of lens fiber (MIP), and lens intrinsic membrane protein 2 gene (LIM2) [23-25], three growth and transcription factor genes: heat shock transcription factor 4 (HSF4), paired-like homeodomain 3 (PITX3), Maf-like protein (MAF) [26,27], and chromatin modifying protein-4B (CHMP4B), Ephrin receptor EphA2 (EPHA2), and Nance-Horan syndrome (NHS) [27].

As crystallin genes account for nearly 90% of the water soluble proteins in lens and the encoded proteins account for around 30% of lens mass, these proteins play an essential roles in maintaining the lens transparency [28] and are good candidate genes for screening in congenital cataract patients. Crystallin mutations accounts for about 50% of the non syndromic cataract [14]. Since the lens is an avascular structure, the crystallins are retained in soluble form through the maintenance of ionic balance by the actions of gap junction proteins which allow the metabolically active epithelium to regulate the precise inter-cellular communication and transport between the lens periphery and its interior. Gap junction proteins (alpha 3 and alpha 8) are expressed in the lens vesicle and mutations in these genes have been reported to lead to cataract [29-31].

Congenital cataract is the most important treatable cause of pediatric blindness in developing countries like India. In this pilot study, we screened 30 cases of congenital cataract for sequence variations in CRYAB, CRYGC, CRYGD, and GJA8. Upon sequence analysis we detected six sequence variations. Out of six, two novel mutations were found in CRYGC (R48H) and GJA8 (L281C). The probable pathogenicity of the mutations found in this study as disease causing is discussed in light of earlier studies.

Methods

Clinical examination and selection of cases

After receiving ethical approval from the institutional review board (IRB#00006862; All India Institute of Medical Sciences, Delhi, India), 30 clinically diagnosed consecutive congenital cataract cases below 3 years of age from northern India, presenting at the Dr. R. P. Centre for Ophthalmic Sciences (AIIMS, New Delhi, India) were enrolled in this study. These congenital cataract cases had no other ocular or systemic abnormalities. Detailed history was taken from parents regarding high fever, TORCHES ([Toxoplasma gondii; T. gondii], rubella virus [RV], cytomegalovirus [CMV], herpes simplex virus [HSV], syphilis [caused by Treponema pallidum]) infection, tuberculosis, exposure to radiation, and drug intake during gestation period. Metabolic tests like serum biochemistry for levels of blood glucose, calcium and phosphorous evaluations, RBC transferase and galactokinase levels and urine test for reducing sugars (galactosemia) and for amino acids (Lowe syndrome) were also done. Cases with known cause of congenital cataract were excluded from the study. Affected status was determined by a history of cataract extraction or ophthalmologic examination. A total of 30 ethnically and age-matched normal individuals without any history of ocular or systemic disorders were enrolled as controls. They had no metabolic, genetic, or ocular disorder on examination by a ophthalmologist and an extensive history was taken regarding family, occupation of parents, any medical problem, and drug intake by parents. Informed consent in accordance with the Declaration of Helsinki was obtained from all participants or their parents and controls.

DNA isolation, PCR amplification and sequence analysis

Genomic DNA was extracted from whole blood samples of all cases and controls, using organic method as described by Sambrook et al. [32] with some modifications. Briefly, 5 ml blood was incubated in 15 ml of Red Cell Lysis Buffer (RCLB) at 4 °C and then centrifuged at 6,861× g for 15 min at 4 °C. Supernatant was discarded and pellet was given repeated washes with RCLB till the pellet became white. The white pellet was re-suspended in 5 ml of DNA extraction buffer, 40 µl of proteinase-K (40 µg/ml) and 300 µl of sodium dodecyl sulphate (SDS). The cocktail was incubated at 55 °C for 2-3 h. The digested proteins were precipitated by adding equal volumes of saturated phenol and chloroform:isoamylalcohol (24:1). The mixture was gently mixed on rotor-mixer for 15-20 min and then centrifuge at 6,861× g for 15 min at 4 °C. In upper viscous layer equal amount of chloroform:isoamylalcohol (24:1) solution was added and again mixed for 15 min in rotor-mixer. The aqueous layer containing the genomic DNA was carefully collected. Then DNA was precipitated with chilled ethanol and centrifugation at 1,461.6× g for 15 min. DNA pellet was washed with 70% ethanol and dissolved in TE buffer. The exon-intron regions of CRYAB, CRYGC, CRYGD, and GJA8 were amplified in congenital cataract patients and controls. PCR amplifications for all primer sets (Table 1) were performed in a 40 μl volume containing 1.0 μl of 20 mM stock solution for each primer (Eurofins Genomics India pvt Ltd, Bangalore, India), 100 ng of genomic DNA, 1 unit of Taq polymerase (Banglore Genei, Bengaluru, Karnataka, India), 0.1 mM of each deoxynucleotide triphosphate (dNTP), and 4 μl of 10× PCR buffer (with 15 mM MgCl2). Amplified PCR products were purified using a gel/PCR DNA fragments extraction kit (Geneaid Biotech Ltd., Sijhih City, Taiwan). Purified PCR products were sent for sequencing to Molecular Cloning Laboratories (South San Francisco, CA). All fragments were sequenced in both forward and reverse directions for confirmation of any nucleotide variation in congenital cataract patients and controls and compared to the Human Genome Reference Sequence (NC_000002.11 and NC_000001.10) provided by the National Center for Biotechnology Information (NCBI), using ClustalW2 (multiple sequence alignment program for DNA; European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, UK).

Computational assessment of missense mutations

We used an evolutionary model to predict the functional consequence of genetic variation in the ATP-binding cassette, sub-family A (ABC1), member 1 gene and tested these predictions through in vitro assessments of protein function [33]. We predicted the functional consequence of each variant using PANTHER. The probability that a given coding variant will cause a deleterious functional change is estimated by the substitution position-specific evolutionary conservation (sub-PSEC) score. SIFT (Sorting Intolerant From Tolerant) analysis tool was also used to predict the functional impact of missense changes identified in this study. SIFT is a sequence homology based tool that sorts intolerant from tolerant amino acid substitutions and predicts whether an amino acid substitution in a protein will have a phenotypic effect [34]. SIFT is based on the premise that protein evolution is correlated with protein function. Positions with normalized probabilities less than 0.05 are predicted to be deleterious and, those greater than or equal to 0.05 are predicted to be tolerated in case of SIFT. We have also used an improved splice site predictor tool [35] to predict whether a nucleotide change is likely to create a splice site.

Protein modeling

The normal and mutant proteins were analyzed for their structure. Prediction of structure differences between wild and mutant were performed using Discovery Studio (DS) 2.0 (Accelrys Inc., San Diego, CA) [36]. The first step in homology modeling method was to find suitable homologus structure (template). Comparative modeling for GJA8 was not possible as the homology model has only 21% sequence identity whereas homology model for human γC-crystallin has 84% sequence identity, thus the comparative modeling for human γC-crystallin was possible using homology model.

Comparative modeling of human γC-crystallin

The best available template for the modeling of 3-D structure of the human γC-crystallin was a high resolution (1.9 Å) crystal structure of mouse γC-crystallin (PDB ID=2V2U) [37]. The sequence identity and similarity between human and mouse γC-crystallin was found to be 84% and 91%, respectively. This template was used to build the homology model of human γC-crystallin using MODELER 9.2 program [38,39] available in Discovery Studio (DS) 2.0 (Accelrys Inc., San Diego, CA), a software package for molecular modeling and simulation. The model with the lowest energy among all the generated models was taken and its stereochemistry checked using the Ramachandran plot. The native model was solvated and further minimized using the available molecular dynamics (MD) simulation protocols to ensure the stability of the generated model. The 3-D model structure of human γC-crystallin mutant (Arg48His) was developed taking the model structure of wild type human γC-crystallin by using the “Build Mutant” protocol and altering the corresponding residue from Arg to His. The built model of mutant was optimized similar to the wild-type human γC-crystallin. The explicit solvent MD simulation was also performed similar to the wild type human γC-crystallin.

Statistical analysis

The correlation coefficient between mutations in crystallin and gap junction protein genes and parameters like degree of opacification, morphology of congenital cataract, and visual acuity were calculated by spearman's test. p-value less than 0.05 is considered as significant. Statistical analyses were performed using graphpad software (GraphPad Software, Inc., La Jolla, CA).

Results

Clinical findings

A total of 30 congenital cataract patients below 3 years of age were enrolled in this study. The mean age of the patients was 1.75±0.19 years (one month to 3 years). The age of onset was recorded as the age at which the disease was first noticed by the parents or first documented by a clinician. All cases were sporadic and were enrolled consecutively as they presented to Dr. R.P. Centre for Ophthalmic Sciences. In this study 20 cases were males and 10 were females. None of the cases were product of consanguineous marriage and all cases had bilateral congenital cataract. The cataract phenotype varied among patients as 66.66% (20/30) of patients had nuclear cataract (Figure 1A), 23.33% (7/30) had zonular/lamellar type cataract (Figure 1B), 6.66% (2/30) had anterior polar cataract (Figure 1C), and 3.3% (1/30) had total cataract (Figure 1D; Table 2). In this study 93% cases were detected with one or the other nucleotide alterations in CRYAB, CRYGC, CRYGD, and GJA8. Six nucleotide variations were detected in patients (Table 3). 66% nucleotide changes were found in crystallin genes (CRYAB, CYRGC, and CRYGD) and 44% were detected in connexin (GJA8). Of the six mutations identified, two were novel and four have been reported.

Summary of mutations in α-crystallin genes

The α-crystallin gene family consists of two similar genes coding for αA-crystallin (CRYAA located on chromosome 21q22.3) and αB-crystallin (CRYAB on chromosome 11q22.1) sharing 57% sequence identity. CRYAB contains 3 exons which encodes a 175 amino acid protein. Direct sequencing of the coding regions and of the flanking intronic sequences of CRYAB revealed one nucleotide change (rs11603779T>G) in the intronic region between exon 2 and 3 of CRYAB. The variation was found in 13.33% (4/30) case of congenital cataract. No nucleotide changes were found in controls. Improved splice site prediction for rs11603779T>G showed that this location is not present at splice site and may not create a splicing error in CRYAB.

Summary of mutations in γ-crystallin genes

The γ-crystallin gene family is mainly located in a cluster of six highly related genes (CRYGA-CRYGF) on human chromosome 2q33–35 and the seventh CRYG gene (CRYGS) on human chromosome 3. Mutations in CRYGC and CRYGD have been associated with congenital and hereditary cataract (Table 4). Direct sequencing of the coding region and of the flanking intronic sequences of CRYGC and CRYGD revealed three sequence variations. One heterozygous nucleotide change (c.G181A) was detected in exon 2 of CRYGC, resulting in the substitution of Arg to His at codon 48 (p.R48H; Figure 2) and was found in 13.33% (4/30) cases of congenital cataract. The multiple sequence alignments generated using FASTA3 (version 3 at the EBI) software showed that the Arg at position 48 of human CRYGC is highly conserved in Macaca mulatta, Canis lupus, Bos taurus, Rattus norvegicus, Mus musculus, and Pan troglodytes (Figure 3). Nucleotide change p.R48H was found to be non-pathogenic on insilico analysis (PANTHER and SIFT; Table 3). However as this change was in a highly conserved domain it may adversely affect protein function. None of the nucleotide changes were detected in control group. Two nucleotide changes in CRYGD; c.A313G in exon 3, resulting in synonymous change (p.R95R; rs2305430) and c.T564C (rs2305429) in the 3′UTR region, were also observed in 28 and 23 patients, respectively.

Summary of mutations in the GJA8 (Connexin-50) gene

Direct sequencing of the amplified fragments of GJA8 in congenital cataract patients identified two single base alterations (p.L268L and p.L281C). The change p.L268L was found in 3.33% (1/30) case of congenital cataract with anterior polar cataract whereas p.L281C (heterozygous) also found in 3.33% (1/30) cases. This case had lamellar/zonular form of cataract with nystagmoid movement. Both the nucleotide alterations (c.C857T and c.T905C) were in the second exon of GJA8. The nucleotide alteration c.T905C resulted in a novel amino acid substitution of leucine to cysteine at codon 281 (p.L281C; Figure 4) whereas the c.C857T nucleotide alteration leads to a synonymous amino acid substitution. GJA8 family protein sequences were obtained from NCBI website and multiple-sequence alignments of GJA8 family proteins from various species were obtained (Figure 5) using FASTA (version 3 at the EBI). This changed a phylogenetically conserved leucine to cysteine at codon 281 (p.L281C). Computational analysis (PANTHER and SIFT) of p.L281C predicted this nucleotide change as a pathogenic variant (Table 3). The remainder of the GJA8 coding sequence showed no change. In addition, these nucleotide changes were not detected in 30 normal unrelated individuals from the same ethnic background.

Comparative modeling study of γC-crystallin

Since the sequence identity between target (human γC-crystallin) and template (mouse γC-crystallin) was 84% (Figure 6), the structural reliability of the generated 3-dimensional homology model was high. The stability of the modeled wild-type and mutant (Arg48His) human γC-crystallin was checked by performing molecular dynamics simulation. These results indicated that the structures were stable. Both the wild type and mutant had very similar conformation with good geometry. The structure of wild type as well as mutant human γC-crystallin consisted of three small helices (each helix contains four residues), four anti-parallel β-strands and connecting loops (Figure 7).

In-silico analysis

PANTHER and SIFT online tools were used for potential functional prediction of mutant proteins. After input the amino acid sequences of the wild-type CRYGC and GJA8 protein and their mutants, the PANTHER scores were −2.72 and −3.97, respectively, whereas the SIFT scores were 1.00 and 0.00, respectively, which meant that the variant (CRYGC:p.R48H) was predicted as “non-pathogenic” and the variant (GJA8:p.L281C) was predicted as “pathogenic” with high confidence. In comparison with the wild-type CRYGC and GJA8 protein, the hydrophobicity of the mutants CRYGC and GJA8 were dramatically increased (Figure 8). The secondary structure of mutant and wild type amino acid sequences of GJA8 were analyzed by Antheprot 2000ver. 6.0 software (IBCP, Lyon,France) which showed that the mutation p.L281C lead to the replacement of random coil with extended loop (Figure 9). This replacement may be sufficient to change the secondary structure of the protein resulting in lens opacification.

Genotype-phenotypes correlation

The genotype-phenotype correlation with the different morphological types of congenital cataract, their severity, visual acuity and different mutations have been tabulated (Table 5). The genotype and phenotype correlation coefficient (r value) between parameters like Degree of opacification, morphology of congenital cataract, visual acuity and mutations, showed no significance and hence no association between mutations and different parameters. Therefore no particular type of cataract was found to be associated with any particular mutant phenotype.

Discussion

The transparency and high refractive index of the lens are achieved by the precise architecture of the fiber cells and the homeostasis of the lens proteins in terms of their concentration, stability, and supramolecular organization [40]. In this pilot study we identified six nucleotide variations (Table 3). Crystallin specific mutations (16.6%) were identified which is similar to the mutations detected in south Indian population [41]. We also detected 16.6% GJA8 specific variations in congenital cataract cases. We identified two non-synonymous novel mutations, p.R48H(4/30) and p.L281C(1/30) in CRYGC and GJA8, respectively. Crystallins (α-, β-, and γ-crystallin) encode the major proportion of water soluble structural proteins of the lens fiber cells and are ubiquitous lens proteins. Functional changes and alteration of crystallin molecular properties could cause the breakdown of the lens microstructure and result in changes in the refractive index and increased light scattering.

Description of mutations in CRYGC and associated phenotypes

It is reported that self aggregation or quaternary structural alteration of γ-crystallin is responsible for the phenotypic association with lens opacification as well as cataractogenesis [42,43]. To the best of our knowledge, four mutations in CRYGC have been reported in the literature (Table 4) [44-47]. The mutation p.R48H involves substitution of highly basic and polar charged Arginine with a neutral and less polar Histidine which may cause conformational changes. Arginine has well spread electron density enabling high solubility. It is a hydrophilic amino acid with a positive charge and lies within the extended strand on the surface of the molecule, interacting with water. Arginine is replaced by histidine, a hydrophobic amino acid compared to arginine. It has been reported that changing the solvation property of an amino acid residue on the surface of the γ-crystallin protein diminishes its solubility [48]. The distorted γC-crystallin may change its folding properties as shown in a study where the COOH-terminal domain folds before and nucleates the folding of the NH2-terminal domain in human γD-crystallin refolding [49]. The relatively loose or partially unfolded structure of mutant γC-crystallin may be susceptible to aggregation and insolubilization, which leads to cataract formation [50]. Another possible consequence of the R48H mutation may be related to the disturbances of the interactions between γC-crystallin and other crystallins [51].

The overall conformation, secondary structure elements and geometry of the conformers of wild-type and mutant human γC-crystallin were mainly similar. The most significant variation was observed in the conformation of the loop 3 region involving residues 47–54 which houses the mutation (Arg48His; Figure 7). The mutant has substitution of the longer and basic Arginine by a shorter Histidine possessing an imidazole ring. Thus this mutation alters the characteristic of this residue in both nature and length which is reflected in the difference in its interactions with neighboring amino acid residues The Arg48 in the wild type interacts with its adjacent acidic residue Glu47 which in turn forms a hydrogen bond with Gln54 thus stabilizing the loop 3 (Figure 10A). Thus the interaction of Glu47 with both residues Arg48 and Gln54 imparts an orientation to loop3. In the case of the mutant the shorter Histidine is no longer able to interact with Glu47 but instead interacts with Gln52 (Figure 10B). This loss of the stabilizing interaction enables Glu47 to adopt a different orientation and it in turn interacts with Arg77 while maintaining its interaction with Gln54. This alteration in the orientation of Glu47 and subsequently its interacting residues in the mutant, results in the modification of the relative orientation of the loop 3 comprising residues 47 to 54.

γ-Crystallins are long lived proteins of the lens and are generally characterized by both high stability and solubility. A key feature of γ-crystallins is that their surfaces are covered in ion pairs. The Arg head groups are most accessible to solvent water and therefore have a profound effect on the surrounding water, particularly important in lens. Since sequences of γ-crystallins are optimized for high solubility and minor changes to the surface can dramatically alter solution interaction properties. So the changes in protein conformations can decrease in solubility and stability of native protein which could result in aggregation of protein and thus opacification of lens. The R48H mutation interferes with the formation of two COOH-terminal Greek key motifs. Although the function of the Greek key motifs has not been elaborated in detail, computer-based analysis suggests that it may be responsible for particular protein–protein interactions in the lens, and it is postulated to be critical in the maintenance of lens transparency. The possible influence of the mutation on the structure as well as the function of γC-crystallin requires further investigation.

Description of mutations in GJA8 and associated phenotypes

GJA8 is located on chromosome 1q21.1 and encodes a 50 kDa protein (connexin 50; Cx50). Cx-50 is a member of the connexin family of proteins that are important in the formation of gap junction channels [52] in human lens which are responsible for direct intercellular transfer of ions and molecules between adjacent cells [53]. Since the eye lens is an avascular structure, it relies heavily on an intercellular communication system constructed of gap junctions for preservation of tissue homeostasis and transparency [3,13]. Cx50 contains four transmembrane domains (M1, M2, M3, and M4) linked by two extracellular loops (E1 and E2), as well as an intracellular loop (CL), and intra-cytoplasmic NH2- and COOH-termini [54]. Mutant connexins are unable to participate in gap junction formation [13] and inhibit channel formation. To date, 22 mutations have been detected in GJA8 in association with congenital cataract (Table 4) [54-69]. We identified a novel missense mutation, p.L281C, in GJA8 in cases with congenital cataract (zonular/lamellar cataract with nystagmus). Different mutations (Table 4) in the connexin gene are often associated pulverulent nuclear opacities [55,56,58,59,61,65,68,69]. However some studies [60,63,67] have detected mutations in GJA8 which were associated with zonular/lamellar cataract phenotype as we found in this study. Insilico analysis of p.L281C mutation showed that the hydrophobicity of the mutant protein increased while the hydrophobic moment decreased (Figure 8). The predicted new characteristics of the mutant protein, which include altered interactions with other proteins, altered regulation activities of GJA8, and altered assemblies, may be the cause of the disease.

Our findings further expand the mutation spectrum of GJA8 and CRYGC in congenital cataract. In summary, this study identified variations in 28 of 30 congenital cataract patients in north Indian population. Crystallin family (α- and γ-crystallin) accounts for 66% of the variations whereas connexins accounts for 44% of the total variations. It is notable that only two variations (CRYGC:p.R48H and GJA8:p.L281C) in CRYGC and GJA8 detected in this study were predicted to be pathogenic which may cause congenital cataract. This study further confirms that CRYGC and GJA8 play a major role in the maintenance of lens transparency and expands the mutation spectrum of both the genes in congenital cataract.

Acknowledgments

We thank all patients and their family members who participated in this study. This study was financially supported by the ICMR (Indian Council of Medical Research, New Delhi, India). Manoj Kumar is a Senior Research Fellow (SRF) of ICMR and gratefully acknowledges ICMR for its financial support.

References

  1. Scott MH, Hejtmancik JF, Wozencraft LA, Reuter LM, Parks MM, Kaiser-Kupfer MI. Autosomal dominant congenital cataract. Interocular phenotypic variability. Ophthalmology. 1994; 101:866-71.8190472. [PMID: 8190472]
  2. Hejtmancik JF. The genetics of cataract: our vision becomes clearer. Am J Hum Genet. 1998; 62:520-5. [PMID: 9497271]
  3. Reddy MA, Francis PJ, Berry V, Bhattacharya SS, Moore AT. Molecular genetic basis of inherited cataract and associated phenotypes. Surv Ophthalmol. 2004; 49:300-15. [PMID: 15110667]
  4. Johnson GJ, editor. The epidemiology of eye disease. 2nd ed. London: Arnold; 2003.
  5. Foster A, Gilbert C, Rahi J. Epidemiology of cataract in childhood: a global perspective. J Cataract Refract Surg. 1997; 23:601-4. [PMID: 9278811]
  6. Bhattacharjee H, Das K, Borah RR, Guha K, Gogate P, Purukayastha S, Gilbert C. Causes of childhood blindness in the northeastern states of India. Indian J Ophthalmol. 2008; 56:495-9. [PMID: 18974521]
  7. Titiyal JS, Pal N, Murthy GV, Gupta SK, Tandon R, Vajpayee RB, Gilbert CE. Causes and temporal trends of blindness and severe visual impairment in children in schools for the blind in North India. Br J Ophthalmol. 2003; 87:941-5. [PMID: 12881329]
  8. Dandona L, Williams JD, Williams BC, Rao GN. Population-based assessment of childhood blindness in southern India. Arch Ophthalmol. 1998; 116:545-6. [PMID: 9565065]
  9. [REMOVED IF= FIELD]9Census of India 2001. India: Registrar General and Census Commissioner, 2001.
  10. Santhiya ST, Kumar GS, Sudhakar P, Gupta N, Klopp N, Illig T, Söker T, Groth M, Platzer M, Gopinath PM, Graw J. Molecular analysis of cataract families in India: new mutations in the CRYBB2 and GJA3 genes and rare polymorphisms. Mol Vis. 2010; 16:1837-47. [PMID: 21031021]
  11. Johar SR, Savalia NK, Vasavada AR, Gupta PD. Epidemiology based etiological study of pediatric cataract in western India. Indian J Med Sci. 2004; 58:115-21. [PMID: 15051906]
  12. Gogate P, Kalua K, Courtright P. Blindness in childhood in developing countries: time for a reassessment? PLoS Med. 2009; 6:e1000177 [PMID: 19997501]
  13. Hejtmancik JF. Congenital cataract and their molecular genetics. Semin Cell Dev Biol. 2008; 19:134-49. [PMID: 18035564]
  14. Shiels A, Bennett TM, Hejtmancik JF. Cat-Map: putting cataract on the map. Mol Vis. 2010; 16:2007-15. [PMID: 21042563]
  15. Graw J. Mouse models of cataract. J Genet. 2009; 88:469-86. [PMID: 20090208]
  16. Berry V, Francis P, Reddy MA, Collyer D, Vithana E, MacKay I, Dawson G, Carey AH, Moore A, Bhattacharya SS, Quinlan RA. Alpha-B crystallin gene (CRYAB) mutation causes dominant congenital posterior polar cataract in humans. Am J Hum Genet. 2001; 69:1141-5. [PMID: 11577372]
  17. Xu J, Wong C, Tan X, Jing H, Zhou G, Song W. Decreasing the homodimer interaction: a common mechanism shared by the deltaG91 mutation and deamidation in betaA3-crystallin. Mol Vis. 2010; 16:438-44. [PMID: 20300566]
  18. Zhou G, Zhou N, Hu S, Zhao L, Zhang C, Qi Y. A missense mutation in CRYBA4 associated with congenital cataract and microcornea. Mol Vis. 2010; 16:1019-24. [PMID: 20577656]
  19. Billingsley G, Santhiya ST, Paterson AD, Ogata K, Wodak S, Hosseini SM, Manisastry SM, Vijayalakshmi P, Gopinath PM, Graw J, Heon E. CRYBA4, a novel human cataract gene, is also involved in microphthalmia. Am J Hum Genet. 2006; 79:702-9. [PMID: 16960806]
  20. Riazuddin SA, Yasmeen A, Yao W, Sergeev YV, Zhang Q, Zulfiqar F, Riaz A, Riazuddin S, Hejtmancik JF. Mutations in betaB3-crystallin associated with autosomal recessive cataract in two Pakistani families. Invest Ophthalmol Vis Sci. 2005; 46:2100-6. [PMID: 15914629]
  21. Sun H, Ma Z, Li Y, Liu B, Li Z, Ding X, Gao Y, Ma W, Tang X, Li X, Shen Y. Gamma-S crystallin gene (CRYGS) mutation causes dominant progressive cortical cataract in humans. J Med Genet. 2005; 42:706-10. [PMID: 16141006]
  22. Song S, Landsbury A, Dahm R, Liu Y, Zhang Q, Quinlan RA. Functions of the intermediate filament cytoskeleton in the eye lens. J Clin Invest. 2009; 119:1837-48. [PMID: 19587458]
  23. Zhou Z, Hu S, Wang B, Zhou N, Zhou S, Ma X, Qi Y. Mutation analysis of congenital cataract in a Chinese family identified a novel missense mutation in the connexin 46 gene (GJA3). Mol Vis. 2010; 16:713-9. [PMID: 20431721]
  24. Berry V, Francis P, Kaushal S, Moore A, Bhattacharya S. Missense mutations in MIP underlie autosomal dominant 'polymorphic' and lamellar cataracts linked to 12q. Nat Genet. 2000; 25:15-7. [PMID: 10802646]
  25. Pras E, Levy-Nissenbaum E, Bakhan T, Lahat H, Assia E, Geffen-Carmi N, Frydman M, Goldman B, Pras E. A missense mutation in the LIM2 gene is associated with autosomal recessive presenile cataract in an inbred IraqiJewish family. Am J Hum Genet. 2002; 70:1363-7. [PMID: 11917274]
  26. Shi X, Cui B, Wang Z, Weng L, Xu Z, Ma J, Xu G, Kong X, Hu L. Removal of Hsf4 leads to cataract development in mice through down-regulation of gamma S-crystallin and Bfsp expression. BMC Mol Biol. 2009; 10:10 [PMID: 19224648]
  27. Wang KJ, Zhu SQ, Cheng J. Progress in pathogenic genes and their functions of congenital cataract. Zhonghua Yan Ke Za Zhi. 2010; 46:280-4. [PMID: 20450675]
  28. Wistow GJ, Piatigorsky J. Lens crystallins: the evolution and expression of proteins for a highly specialized tissue. Annu Rev Biochem. 1988; 57:479-504. [PMID: 3052280]
  29. Li L, Cheng C, Xia CH, White TW, Fletcher DA, Gong X. Connexin mediated cataract prevention in mice. PLoS One. 2010; 5:e12624 [PMID: 20844585]
  30. Hu S, Wang B, Zhou Z, Zhou G, Wang J, Ma X, Qi Y. A novel mutation in GJA8 causing congenital cataract-microcornea syndrome in a Chinese pedigree. Mol Vis. 2010; 16:1585-92. [PMID: 20806042]
  31. Yan M, Xiong C, Ye SQ, Chen Y, Ke M, Zheng F, Zhou X. A novel connexin 50 (GJA8) mutation in a Chinese family with a dominant congenital pulverulent nuclear cataract. Mol Vis. 2008; 14:418-24. [PMID: 18334966]
  32. Sambrook J, Fritsch EF, Maniatus T. Molecular cloning: a laboratory manual, second ed., Cold Spring Harbour Labortory Press, 1989.
  33. Brunham LR, Singaraja RR, Pape TD, Kejariwal A, Thomas PD, Hayden MR. Accurate prediction of the functional significance of single nucleotide polymorphisms and mutations in the ABCA1 gene. PLoS Genet. 2005; 1:e83 [PMID: 16429166]
  34. Ng PC, Henikoff S. SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003; 31:3812-4. [PMID: 12824425]
  35. Reese MG, Eeckman FH, Kulp D, Haussler D. Improved splice site detection in Genie. J Comput Biol. 1997; 4:311-23.
    .9278062
    [PMID: 9278062]
  36. Discovery Studio 2.0.. Molecular modeling program package. Accelrys Software Inc; San Diego, CA USA 2006;
  37. Purkiss AG, Bateman OA, Wyatt K, Wilmarth PA, David LL, Wistow GJ, Slingsby C. Biophysical Properties of? C-Crystallin in Human and Mouse Eye Lens: The Role of Molecular Dipoles. J Mol Biol. 2007; 372:205-22. [PMID: 17659303]
  38. Sali A, Blundell TL. Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993; 234:779-815. [PMID: 8254673]
  39. Sali A, Pottertone L, Yuan F, van Vlijmen H, Karplus M. Evaluation of comparative protein modeling by MODELLER. Proteins. 1995; 23:318-26. [PMID: 8710825]
  40. Graw J. Cataract mutations and lens development. Prog Retin Eye Res. 1999; 18:235-67. [PMID: 9932285]
  41. Devi RR, Yao W, Vijayalakshmi P, Sergeev YV, Sundaresan P, Hejtmancik JF. Crystallin gene mutations in Indian families with inherited pediatric cataract. Mol Vis. 2008; 14:1157-70. [PMID: 18587492]
  42. Pande A, Pande J, Asherie N, Lomakin A, Ogun O, King JA, Lubsen NH, Walton D, Benedek GB. Molecular basis of a progressive juvenile-onset hereditary cataract. Proc Natl Acad Sci USA. 2000; 97:1993-8. [PMID: 10688888]
  43. Pande A, Pande J, Asherie N, Lomakin A, Ogun O, King J, Benedek GB. Crystal cataract: human genetic cataract caused by protein crystallization. Proc Natl Acad Sci USA. 2001; 98:6116-20. [PMID: 11371638]
  44. Devi RR, Vijayalakshmi P. Novel mutations in GJA8 associated with autosomal dominant congenital cataract and microcornea. Mol Vis. 2006; 12:190-5. [PMID: 16604058]
  45. Hansen L, Yao W, Eiberg H, Kjaer KW, Baggesen K, Hejtmancik JF, Rosenberg T. Genetic heterogeneity in microcornea-cataract: five novel mutations in CRYAA, CRYGD, and GJA8. Invest Ophthalmol Vis Sci. 2007; 48:3937-44. [PMID: 17724170]
  46. Héon E, Priston M, Schorderet DF, Billingsley GD, Girard PO, Lubsen N, Munier FL. The? -crystallins and human cataracts: a puzzle made clearer. Am J Hum Genet. 1999; 65:1261-7. [PMID: 10521291]
  47. Roshan M, Vijaya PH, Lavanya GR, Shama PK, Santhiya ST, Graw J, Gopinath PM, Satyamoorthy K. A novel human CRYGD mutation in a juvenile autosomal dominant cataract. Mol Vis. 2010; 16:887-96. [PMID: 20508808]
  48. Evans P, Wyatt K, Wistow GJ, Bateman OA, Wallace BA, Slingsby C. The P23T cataract mutation causes loss of solubility of folded gammaD-crystallin. J Mol Biol. 2004; 343:435-44. [PMID: 15451671]
  49. Flaugh SL, Kosinski-Collins MS, King J. Interdomain sidechain interactions in human gammaD crystallin influencing folding and stability. Protein Sci. 2005; 14:2030-43. [PMID: 16046626]
  50. Talla V, Narayanan C, Srinivasan N, Balasubramanian D. Mutation causing self-aggregation in human gammaC-crystallin leading to congenital cataract. Invest Ophthalmol Vis Sci. 2006; 47:5212-7. [PMID: 17122105]
  51. Fu L, Liang JJ. Detection of protein-protein interactions among lens crystallins in a mammalian two-hybrid system assay. J Biol Chem. 2002; 277:4255-60. [PMID: 11700327]
  52. Berthoud VM, Minogue PJ, Guo J, Williamson EK, Xu X, Ebohara L, Beyer EC. Loss of function and impaired degradation of a cataract-associated mutant connexin 50. Eur J Cell Biol. 2003; 82:209-21. [PMID: 12800976]
  53. Bakthavachalu B, Kalanke S, Galande S, Ramanamurthy B, Parab P, Kohale KN, Seshadri V. Dense cataract and microphthalmia (dcm) in BALB/c mice is caused by mutations in the GJA8 locus. J Genet. 2010; 89:147-54. [PMID: 20861565]
  54. Willoughby CE, Arab S, Gandhi R, Zeinali S, Arab S, Luk D, Billingsley G, Munier FL, Heon E. A novel GJA8 mutation in an Iranian family with progressive autosomal dominant congenital nuclear cataract. J Med Genet. 2003; 40:e124 [PMID: 14627691]
  55. Wang K, Wang B, Wang J, Zhou S, Yun B, Suo P, Cheng J, Ma X, Zhu S. A novel GJA8 mutation (p.I31T) causing autosomal dominant congenital cataract in a Chinese family. Mol Vis. 2009; 15:2813-20. [PMID: 20019893]
  56. Vanita V, Singh JR, Singh D, Varon R, Sperling K. A novel mutation in GJA8 associated with jellyfish-like cataract in a family of Indian origin. Mol Vis. 2008; 14:323-6. [PMID: 18334946]
  57. Arora A, Minogue PJ, Liu X, Addison PK, Russel-Eggitt I, Webster AR, Hunt DM, Ebihara L, Beyer EC, Berthoud VM, Moore AT. A novel connexin50 mutation associated with congenital nuclear pulverulent cataracts. J Med Genet. 2008; 45:155-60. [PMID: 18006672]
  58. Lin Y, Liu NN, Lei CT, Fan YC, Liu XQ, Yang Y, Wang JF, Liu B, Yang ZL. A novel GJA8 mutation in a Chinese family with autosomal dominant congenital cataract. Zhonghua Yi Xue Yi Chuan Xue Za Zhi. 2008; 25:59-62. [PMID: 18247306]
  59. Berry V, Mackay D, Khaliq S, Francis PJ, Hameed A, Anwar K, Mehdi SQ, Newbold RJ, Ionides A, Shiels A, Moore T, Bhattacharya SS. Connexin 50 mutation in a family with congenital “zonular nuclear” pulverulent cataract of Pakistani origin. Hum Genet. 1999; 105:168-70. [PMID: 10480374]
  60. Zheng JQ, Ma ZW, Sun HM. A heterozygous transversion of connexin 50 in a family with congential nuclear cataract in the northeast of China. Zhonghua Yi Xue Yi Chuan Xue Za Zhi. 2005; 22:76-8. [PMID: 15696487]
  61. Vanita V, Hennies HC, Singh D, Nürnberg P, Sperling K, Singh JR. A novel mutation in GJA8 associated with autosomal dominant congenital cataract in a family of Indian origin. Mol Vis. 2006; 12:1217-22. [PMID: 17110920]
  62. Shiels A, Mackay D, Ionides A, Berry V, Moore A, Bhattacharya S. A missense mutation in the human connexin50 gene (GJA8) underlies autosomal dominant “zonular pulverulent” cataract, on chromosome 1q. Am J Hum Genet. 1998; 62:526-32. [PMID: 9497259]
  63. Vanita V, Singh JR, Singh D, Varon R, Sperling K. A mutation in GJA8 (p.P88Q) is associated with “balloon-like” cataract with Y-sutural opacities in a family of Indian origin. Mol Vis. 2008; 14:1171-5. [PMID: 18587493]
  64. Ponnam SP, Ramesha K, Tejwani S, Ramamurthy B, Kannabiran C. Mutation of the gap junction protein alpha 8 (GJA8) gene causes autosomal recessive cataract. J Med Genet. 2007; 44:e85 [PMID: 17601931]
  65. Polyakov AV, Shagina IA, Khlebnikova OV, Evgrafov OV. Mutation in the connexin 50 gene (GJA8) in a Russian family with zonular pulverulent cataract. Clin Genet. 2001; 60:476-8. [PMID: 11846744]
  66. Schmidt W, Klopp N, Illig T, Graw J. A novel GJA8 mutation causing a recessive triangular cataract. Mol Vis. 2008; 14:851-6. [PMID: 18483562]
  67. Yan M, Xiong CL, Ye SQ, Chen YM, Ke M, Zheng F, Zhou X. A novel connexin 50 (GJA8) mutation in a Chinese family with a dominant congenital pulverulent nuclear cataract. Mol Vis. 2008; 14:418-24. [PMID: 18334966]
  68. Li F, Wang S, Gao C, Liu S, Zhao B, Zhang M, Huang S, Zhu S, Ma X. Mutation G61C in the CRYGD gene causing autosomal dominant congenital coralliform cataracts. Mol Vis. 2008; 14:378-86. [PMID: 18334953]
  69. Stephan DA, Gillanders E, Vanderveen D, Freas-Lutz D, Wistow G, Baxevanis AD, Robbins CM, van Auken A, Quesenberry MI, Bailey-Wilson J, Juo SH, Trent JM, Smith L, Brownstein MJ. Progressive juvenile-onset punctuate cataracts caused by mutation of the gD-crystallin gene. Proc Natl Acad Sci USA. 1999; 96:1008-12. [PMID: 9927684]
  70. Messina-Baas OM, Gonzalez-Huerta LM, Cuevas-Covarrubias SA. Two affected siblings with nuclear cataract associated with a novel missense mutation in the CRYGD gene. Mol Vis. 2006; 12:995-1000. [PMID: 16943771]
  71. Santhiya ST, Shyam Manohar M, Rawlley D, Vijayalakshmi P, Namperumalsamy P, Gopinath PM, Loster J, Graw J. Novel mutations in the gamma-crystallin genes cause autosomal dominant congenital cataracts. J Med Genet. 2002; 39:352-8. [PMID: 12011157]
  72. Kmoch S, Brynda J, Asfaw B, Bezouska K, Novák P, Rezácová P, Ondrová L, Filipec M, Sedlácek J, Elleder M. Link between a novel human? D-crystallin allele and a unique cataract phenotype explained by protein crystallography. Hum Mol Genet. 2000; 9:1779-86. [PMID: 10915766]
  73. Zhang LY, Gong B, Tong JP, Fan DS, Chiang SW, Lou D, Lam DS, Yam GH, Pang CP. A novel gammaD-crystallin mutation causes mild changes in protein properties but leads to congenital coralliform cataract. Mol Vis. 2009; 15:1521-9. [PMID: 19668596]
  74. Santana A, Waiswol M, Arcieri ES, Cabral de Vasconcellos JP, Barbosa de Melo M. Mutation analysis of CRYAA, CRYGC, and CRYGD associated with autosomal dominant congenital cataract in Brazilian families. Mol Vis. 2009; 15:793-800. [PMID: 19390652]
  75. Zhang LY, Yam GH, Fan DS, Tam PO, Lam DS, Pang CP. A novel deletion variant of gammaD-crystallin responsible for congenital nuclear cataract. Mol Vis. 2007; 13:2096-104. [PMID: 18079686]
  76. Ren Z, Li A, Shastry BS, Padma T, Ayyagari R, Scott MH, Parks MM, Kaiser-Kupfer MI, Hejtmancik JF. A 5-base insertion in the gammaC-crystallin gene is associated with autosomal dominant variable zonular pulverulent cataract. Hum Genet. 2000; 106:531-7. [PMID: 10914683]
  77. Yao K, Jin C, Zhu N, Wang W, Wu R, Jiang J, Shentu X. A nonsense mutation in CRYGC associated with autosomal dominant congenital nuclear cataract in a Chinese family. Mol Vis. 2008; 14:1272-6. [PMID: 18618005]