Article Text

PDF

Gamma-D crystallin gene (CRYGD) mutation causes autosomal dominant congenital cerulean cataracts
  1. E Nandrot1,
  2. C Slingsby2,
  3. A Basak2,
  4. M Cherif-Chefchaouni3,
  5. B Benazzouz4,
  6. Y Hajaji4,
  7. S Boutayeb4,
  8. O Gribouval5,
  9. L Arbogast1,
  10. A Berraho3,
  11. M Abitbol1,*,
  12. L Hilal4,*
  1. 1Centre de Recherches Thérapeutiques en Ophtalmologie de la Faculté de Médecine Necker, EA No 2502 du Ministére de la Recherche et de l’Enseignement Supérieur, Université René Descartes, Paris, France
  2. 2Birkbeck College, Department of Crystallography, Malet Street, London WC1E 7HX, UK
  3. 3Département d’Oto-Neuro-Ophtalmologie, Service d’Ophtalmologie B, Hôpital des Spécialités, Rabat, Morocco
  4. 4Laboratoire de Génétique et Biologie Moléculaire de la Faculté des Sciences, Université Ibn Tofail, Kénitra, Morocco
  5. 5Unité INSERM 423, CHU Necker-Enfants Malades, Paris, France
  1. Correspondence to:
 Dr M Abitbol, Centre de Recherches Thérapeutiques en Ophtalmologie, Faculté de Médecine Necker- Enfant Malades, 156 Rue de Vaugirard, 75015 Paris cedex, France; 
 abitbol{at}necker.fr or Pr L Hilal, Laboratoire de Génétique et Biologie Moléculaire, Département de Biologie, Faculté des Sciences, Université Ibn Tofail, Kénitra, Morocco; lhilal{at}yahoo.fr

Abstract

Congenital cataracts are a major cause of bilateral visual impairment in childhood. We mapped the gene responsible for autosomal congenital cerulean cataracts to chromosome 2q33–35 in a four generation family of Moroccan descent. The maximum lod score (7.19 at recombination fraction θ=0) was obtained for marker D2S2208 near the γ-crystallin gene (CRYG) cluster. Sequencing of the coding regions of the CRYGA, B, C, and D genes showed the presence of a heterozygous C>A transversion in exon 2 of CRYGD that is associated with cataracts in this family. This mutation resulted in a proline to threonine substitution at amino acid 23 of the protein in the first of the four Greek key motifs that characterise this protein. We show that although the x ray crystallography modelling does not indicate any change of the backbone conformation, the mutation affects a region of the Greek key motif that is important for determining the topology of this protein fold. Our data suggest strongly that the proline to threonine substitution may alter the protein folding or decrease the thermodynamic stability or solubility of the protein. Furthermore, this is the first report of a mutation in this gene resulting in autosomal dominant congenital cerulean cataracts.

Statistics from Altmetric.com

Congenital cataracts (ADCC) are a significant cause of visual impairment in childhood. The estimated prevalence is 1–6/10 000 live births.1,2 In non-consanguineous populations, most cases of inherited non-syndromic cataracts show an autosomal dominant (AD) mode of transmission, but X linked and autosomal recessive forms have also been observed. Inherited cataracts are clinically highly heterogeneous and show considerable inter- and intrafamilial variability.3–5 To date, more than 15 independent loci and genes on different chromosomes have been shown to be associated with autosomal dominant congenital cataracts (ADCC),2,6–9 and causative mutations have been identified in several distinct genes, including genes encoding crystallins,10–18 lens specific connexins (gap junction proteins),19–22 aquaporin (MIP),23,24 and cytoskeletal structural proteins,25,26 as well as in the homeobox gene PITX327 and the heat shock transcription factor 4 gene (HSF4).28 Congenital cataracts are genetically heterogeneous as clinically identical cataracts have been mapped to different loci.2,29 Conversely, an identical mutation can result in different phenotypes.5,10,14

Cerulean cataracts, which were first described by Vogt in 1922,30 are characterised by early onset, progressive development of bilateral lens opacities and complete penetrance. Affected newborns appear to be asymptomatic until the age of 18 to 24 months. At this age they can be clinically diagnosed by slit lamp examination, which shows the presence of tiny blue or white opacities that first develop in the superficial layers of the fetal lens nucleus. Armitage et al31 suggested that cerulean cataracts should be classified as developmental cataracts rather than as congenital cataracts. The opacities spread throughout the adult lens nucleus and cortex and form concentric layers with radially orientated central lesions.

Two cerulean ADCCs (ADCCC) have been mapped in distinct, large families: one to 17q24 (CCA1, MIM 115660) and the other to 22q11.2-q12.2 (CCA2, MIM 601547), which is close to the β-crystallin gene (CRYBB2, MIM 123620).31,32 Litt et al10 identified a chain termination mutation in the CRYBB2 gene and showed that this mutation is associated with the CAA2 type. However, no mutated genes have yet been found to be associated with the CCA1 type. We recently reported a large Moroccan pedigree with an unusual form of cerulean blue dot cataract, characterised by a very early onset of lens opacities and a rapid progression towards blindness in the absence of rapid surgical intervention. In this family, we excluded any linkage to previously mapped loci for ADCCCs (17q24 and 22q).33 We report here the linkage of the ADCCC affecting this Moroccan family to a region of chromosome 2q33-q35, spanning the γ-crystallin gene cluster. We also identified a missense mutation in the CRYGD gene, causing this precocious and rapidly progressing form of ADCCC.

MATERIALS AND METHODS

Family data

We have studied a large Moroccan family with ADCCC. The family consists of 19 affected and 24 unaffected subjects spanning four generations. Appropriate informed consent was obtained from all patients participating in the study. Clinical and ophthalmological examination of patients affected by cerulean cataract as well as unaffected family members were carried out as described previously. The diagnosis of cerulean cataract was confirmed in each affected patient by four independent ophthalmologists. There was no history of other ocular or systemic abnormalities in this family. The clinical aspects of the present study have been published recently.33 Blood samples were collected from 43 family members. Genomic DNA was extracted by standard techniques.

Genotyping and linkage analysis

The genotyping was performed as described previously33 using firstly 30 microsatellite markers, corresponding to 13 known candidate loci for autosomal dominant congenital cataracts, and then 21 markers localised to 2q33–35. The oligonucleotide primer sequences were taken from Généthon (http://www.genethon.fr). Two point disease to marker linkage analysis was conducted by MLINK from the FASTLINK (version 3.0P) software package.34 The mode of inheritance was considered to be autosomal dominant with full penetrance. The gene frequency was set at 1/20 000. As the allelic frequencies of the polymorphic markers were unknown in the Moroccan population, they were considered to be equal. Order and genetic distances were taken from the Marshfield database (http://research.marshfieldclinic.org) and ensembl genome data resources (http://www.ensembl.org). Multipoint analysis was computed using Genehunter software.35

DNA sequencing

Genomic DNA samples from all affected and unaffected family members of the pedigree and from control subjects were screened for mutations in CRYGA, B, C, and D by direct cycle sequencing. Gene specific PCR primers were used to amplify the three exons and flanking introns sequences of CRYGA, CRYGB (sequences available upon request), CRYGC, and CRYGD.12,15 PCR products were purified by means of “invitrogen” rapid PCR purification systems and were sequenced on ABIprism A310 and ABI A377 automated sequencer (PE Biosystems, USA) using the original and additional internal primers.

Molecular modelling

A model of the P23T mutant structure was built based on the human γD crystal coordinates using the program “O”36 by placing the mutant side chain in the most favourable conformation using the side chain database within the program and using the lego-side chain option. The model was energy minimised using 200 cycles of Powell energy minimisation using the CNS program37 and the quality of the model was checked with Procheck.38

RESULTS

We analysed, clinically and genetically, a large Moroccan family affected by ADCCC. The congenital cerulean cataract phenotype was clearly distinct from the phenotypes of the lamellar, coralliform, aculeiform, Coppock-like, and polymorphic congenital cataracts phenotypes.11,15,39–41

Forty-three members of the pedigree, including 20 affected subjects, 17 unaffected family members, and six unaffected spouses (fig 1A), were genotyped. The preliminary linkage analysis performed in this family using 30 microsatellite markers allowed us to exclude 12 candidate loci for ADCC (1pter-p36.13, 1q21.2, 1q21.2, 3p21.2–22.3, 10q23.3–25, 12q12–14.1, 13q11–12, 16q22.1, 17p12–13, 17q11.2-q12.2, 17q24, 21q22.3, and 22q11.2-q12.2) (data not shown). Positive two point lod scores were obtained for markers D2S72 (Zmax=3.47 at θ=0.01) and CRYGA (Zmax=1.75 at θ=0.00) (table 1). D2S72 and CRYGA flank the 4.5 cM 2q33–35 interval, where the loci corresponding to polymorphic congenital cataracts (MIM 601286),39 Coppock-like cataracts (MIM 123660, MIM 116200),40 juvenile onset punctate cataracts (MIM 123690),11 aculeiform congenital cataracts (MIM 115700),41 and variable zonular pulverulent cataracts15 have been previously mapped. Twenty-one additional markers spanning this critical interval were subsequently used for further genotyping of all family members. The two point lod scores for 10 of these 21 markers are summarised in table 1. Significant positive lod scores (Zmax>3 at θ=0) were found for 11 markers. The maximum two point lod score was obtained with marker D2S2208 (Zmax=7.19 at θ=0.00) (table 1). Multipoint analysis with the most informative markers confirmed that this locus mapped to chromosome 2q33-q35 (data not shown).

Table 1

Two point lod scores (Z) for linkage between autosomal dominant congenital cerulean cataract CCA3 and markers in the chromosomal region 2q33–35

Figure 1

(A) Pedigree of the large Moroccan family affected by autosomal dominant congenital cerulean cataracts (CCA3). (B) a Markers of the chromosomal region 2q33-q35 used for linkage analysis. b Common haplotype present in all affected family members. c Haplotypes observed in patients III.6, III.7, III.11, IV.5 and IV.21, displaying critical recombination events. Asterisks indicate subjects who underwent linkage analysis.

Haplotype analysis showed that the affected patients of this family shared a common haplotype involving 12 markers (D2S2237, D2S155, D2S2358, D2S355, D2S2192, D2S325, D2S2242, D2S2208, D2S157, CRYGA, D2S22322, D2S128). Critical recombination events were detected in affected subjects III.6, III.7, III.11, IV.21, and IV.5. This allowed us to define a disease gene containing an interval of about 11.5 cM between markers D2S72 and D2S2361 (fig 1B). All affected subjects had an affected parent, and no unaffected subjects carried the disease haplotype. Thus, penetrance appears to be virtually complete in this family. Several candidate genes have been mapped in this interval, the obvious one being the γ-crystallin gene cluster, CRYG, which has been mapped to 2q33–35.42 Another crystallin gene, CRYBA2, has been mapped to the 2q34–36 region.43CRYBA2 is 10.5 cM telomeric to CRYGA and 1 cM telomeric to D2S173. PAX3, a major developmental gene, is localised 1 cM telomeric to D2S126 and is thus outside the critical interval. The observation of recombination events, which are centromeric to these two genes, in our family excluded the CRYBA2 and PAX3 genes as candidate genes for the cerulean cataracts affecting this family. Thus, this linkage analysis identified a third locus (CCA3) associated with the cerulean blue dot cataract type. This type differs from those previously described and is characterised by early onset (diagnosed at birth) and by a faster progression.

Given our linkage analysis results and the fact that γ-crystallin genes are expressed early in development44,45 and are associated with hereditary cataracts in mice and humans,2,6,46–49 we focused on these genes.. The analysis of the sequences of all three exons of CRYGA, CRYGB, CRYGC, and CRYGD did not show any disease causing mutations in the CRYGA, CRYGB, or CRYGC gene in any of the members of this large family. Only the single nucleotide polymorphisms previously described were identified, like 2437C/T (M19364) (P64P) in exon 2 of CRYGB,12 5391C>A (M193364) (L111I) in exon 3 of CRYGB,17 18229A>G (M19364) (IVS-70 intron A) in CRYGC,12 and 286T>C (K03005) (Y16Y) in exon 2 of CRYGD.12 A unique heterozygous C>A transversion was identified at nucleotide 305 in exon 2 (K03005) of the CRYGD gene in all affected family members. This transversion led to the replacement of a proline residue at amino acid 23 by a threonine (P23T) in the first “Greek key motif” (motif 1) of the γD-crystallin protein. This substitution was not found in any of the unaffected members of the pedigree. This mutation cosegregated strictly with the trait in the kindred studied with complete penetrance. This nucleotide sequence alteration was not found in a panel of 100 unrelated control DNA samples collected from normal Moroccan subjects or in a panel of 100 unrelated control DNA samples collected from people of various ethnic backgrounds. An alignment of γ-crystallin protein sequences from different species showed that Pro 23 is conserved in all major human γ-crystallins as well as in rat and mouse γA, γB, and γC-crystallins and in bovine γB and γC-crystallins (fig 2C). Furthermore, a threonine is never found in any crystallin, whatever the species analysed, at this position. The identified CRYGD proline to a threonine (P23T) substitution at amino acid 23 is a non-conservative amino acid change as proline has a hydrophobic side chain.

Figure 2

Mutation analysis of CRYGD. (A) Sequence chromatograms of the wild type allele showing a proline (CCC) at amino acid 23. (B) Sequence chromatograms of the mutant allele showing the C to A transversion that changed proline 23 to threonine. (C) Sequence alignment of the γ-crystallin proteins in different species. Shown below are the amino acid residues 12–33 and 45–53 of human γD-crystallin. Proline 23 of human CRYGD protein is underlined and appears in bold type. The arrow indicates the P23T mutated position. Residues in the other γ-crystallins identical to those at homologous positions of human CRYGD are indicated by dots.

DISCUSSION

In this report, after excluding most known loci corresponding to ADCC, we identified a third locus (CCA3) on 2q33–35, associated with the cerulean blue dot cataract type in a large Moroccan family. We then found a C>A transversion in exon 2 of GRYGD only in all affected members of the large family studied. Crystallin genes encode a superfamily of major soluble structural proteins in the lens. There are three major classes of crystallins in humans, the α-, β-, and γ-crystallins.50,51 All three types of crystallins are β pleated sheets. Protein sequence analysis has shown homology between the β- and γ-crystallins. The γ-crystallin gene cluster comprises six genes; γA, γB, γC, γD, γE, and γF. Different levels of expression of these genes are associated with the lens refractive index characteristics of the vertebrate enabling adaptation to varying optical requirements.52 In mammals, these genes each consist of three exons.53,54 Only γC- and γD encode abundant lens γ-crystallins in humans.44,45 γE and γF are pseudogenes with in frame stop codons. γD is one of the only two γ-crystallins to be expressed at high concentrations in the fibre cells of the embryonic human lens. These cells subsequently form the lens nucleus fibres.

The identified C>A transversion found exclusively in all affected members of the pedigree results in a P23T substitution. Proline has a hydrophobic side chain often associated with turns in polypeptide chain conformation, whereas the threonine side chain has both hydrophilic and hydrophobic functions. Furthermore, a threonine is never found in any crystallin, whatever the species analysed, at this position. Proline is known to have a strong ability to break β-strands in soluble proteins.55,56 These data strongly suggest that the 305C>A transversion is indeed the CCA3 causing mutation rather than a rare benign polymorphism. It is noteworthy that the same mutation has recently been suggested to cause an autosomal dominant congenital form of lamellar cataracts in a proband and her affected father in a small nuclear Indian family.17 The difference between the phenotype reported in this Indian family and that observed in our Moroccan family may be related to the effect of an unknown modifier gene or to sequence variations within regulatory regions that could affect the expression of the γD-crystallin gene. This phenotypic heterogeneity caused by an identical genetic crystallin mutation is strongly reminiscent of that previously reported for the βB2-crystallin gene. The same chain termination mutation in the βB2-crystallin gene caused autosomal dominant congenital Coppock-like cataracts in a four generation Swiss family14 and autosomal dominant congenital cerulean cataracts (CCA2) in a very large American family.10,32,57 These data further strengthen the hypothesis that the γD-crystallin P23T mutation is pathogenic and prompted us to try to obtain further insights into the molecular mechanisms underlying cataractogenesis in the large Moroccan kindred studied.

The 3D structures of several members of the monomeric γ-crystallin family are known, showing that the polypeptide folds into four similar “Greek key motifs”. These motifs then form two similar domains joined by a linker (recently reviewed in Jaenicke et al58). The 3D structure of human γD-crystallin has now been determined at a very high resolution (1.25 Å) by x ray crystallography.59 We may try to gain some insight into the pathogenicity of the P23T γD-crystallin mutation from the crystallographic structure of the native human protein. Its conformation is very similar to that of γ-crystallins from other species60 and, as expected, very similar to that of the orthologous bovine γD.61 Each γ-crystallin domain is formed from two consecutive Greek key motifs and each motif comprises four consecutive β-strands. Thus, the N-terminal domain can be described as comprising a1 b1 c1 d1 and a2 b2 c2 d2 β-strands. These eight β-strands form two β-sheets (b1 a1 d1 c2 and b2 a2 d2 c1) that pack together to form a β-sandwich domain. Although β-sandwich domains are common in proteins, in βγ-crystallins they are characterised by their high internal conformational symmetry and by a conserved folded hairpin structure between a and b β-strands. The mutation affects a region of the first Greek key motif that is located immediately after the c1 strand, a region that is important for determining the topology of this protein fold, as after crossing from one β-sheet to the next the polypeptide chain has to get back to the original β-sheet to place the d1 strand within it. Obviously, a single Greek key cannot fold on its own; it is dependent on the concomitant folding of its partner motif. The degree of sequence conservation of residue 23 in motif 1 and the corresponding region in the other three motifs gives an idea of the importance of this region. Although not absolutely conserved, this residue tends to be a proline or a serine in motifs 1, 2, and 3, whereas it is an arginine in motif 4. The mutation site thus substitutes a threonine residue into a strategic position in terms of the fold topology, a site that has been occupied, however, by a very similar residue, serine. Fig 3 shows the N-terminal domain along the second β-sheet from the edge of the short c1 strand and centred around residue 23. In wild type human γD-crystallin, the carbonyl oxygen of Pro 23 forms a hydrogen bond with the backbone amide nitrogen of Asn 49 (fig 3B). This backbone interaction helps to stabilise the folded ab hairpin onto the sheet and is thus an important tertiary structure determinant. However, as it involves only backbone interactions, it is not affected by side chain alterations, assuming that in the mutant structure the backbone is in the same conformation as in the native structure. When threonine is modelled into this position, the new side chain is readily accommodated into the space without changing the backbone conformation (fig 3C). Threonine is more polar than proline as it has a polar backbone amide nitrogen and a side chain hydroxyl function. In this predicted model the mutant side chain appears to strengthen the interaction with tyrosine 50. The aromatic side chain of this residue can play an important role in stabilising the folded ab hairpin. Threonine 23 can make a stronger hydrogen bond with the side chain hydroxyl function of Tyr 50 using its backbone amide as well as an additional hydrogen bond with its own side chain hydroxyl. It is difficult to imagine that the mutation has such a catastrophic effect on the 3D structure of the human γD-crystallin based on this model. There are three possible scenarios: the protein does not fold properly because a side chain with a branched Cβ atom has a propensity to extend the β-strand rather than make the necessary turn; the substitution of threonine decreases the thermodynamic stability of the N-terminal domain in a subtle way so that the protein is more prone to denaturation; or thirdly the threonine side chain somehow decreases the solubility of the protein in water. The discovery of this mutation site, with its major clinical consequences, challenges our understanding of crystallin structure. This is clearly shown also by studies on other cataractogenic mutations of human γD-crystallin; recent biophysical studies suggest that the R14C γD-crystallin mutant increases the phase separation temperature62 and that the two γD-crystallin mutant proteins, R58H and R36S, crystallise far more readily than the native protein.63 However, none of these mutant proteins has been shown to change the conformation of the protein significantly. In contrast, the T5P CRYGC mutant protein induces obvious folding changes and decreases conformational stability.64 Thus, further biophysical characterisation of the P23T mutant protein is important to provide more information about the nature of the disease and improve our understanding of the molecular phenotype of lens crystallins.

Figure 3

(A) Cartoon topology diagram of the characteristic domain structure found in all βγ-crystallins showing how two “Greek key motifs”, each comprising four β-strands (a, b, c, and d), associate to form two β-sheets. The red circle shows the site of the P23T mutation in the N-terminal domain of γD-crystallin. (B) Schematic view of the protein fold determined by x ray crystallography centred on residue P23 and looking along β-sheet 2 of the normal CRYGD protein. (C) The same view centred on the mutation site but based on a predicted structure modelled on the crystallographic coordinates of the normal CRYGD protein. Light blue = β-sheet strands, green = loops connecting β-sheet strands, dashed green lines = hydrogen bonds between the region adjacent to the β-sheet c strand and residues in the folded ab β-hairpin, dashed yellow lines = local hydrogen bond in the region of the β-sheet c strand, red = oxygen atoms and blue = nitrogen atoms.

Acknowledgments

We are grateful to the members of the family for their participation in this study. We warmly thank Dr Jayanti Pande for critical reading of the manuscript. We thank Université René Descartes and Professor Patrick Berche, Dean of Faculté de Médecine Necker for his permanent support for the CERTO Laboratory. We thank Ibn Tofail University, Faculty of Sciences of Kénitra and Centre National pour la Recherche Scientifique et Technique of Morocco. We thank Professor Jean-Louis Dufier for his constant encouragement and efforts for the development of the CERTO Laboratory. We also thank Dr Cécile Marsac and Dr Maurice Menasche for their scientific support. We thank all members of the medical staff of the Department of Ophthalmology B (Hôpital des Spécialités de Rabat, Morocco). This work was mainly supported by a grant to the CERTO Laboratory from RETINA France and also by a “Programme d’Appui à la Recherche Scientifique” (Pars Médecine No 151), by a research grant from INSERM (Programme PROGRES), and by several grants from Faculté de Médecine Necker and Université René Descartes. Emeline Nandrot has been supported by a PhD grant from RETINA-France for three years and by a PhD grant for one year from Fondation pour la Recherche Médicale. We gratefully acknowledge the finacial support of the Medical Research Council, London.

REFERENCES

View Abstract

Footnotes

  • * The last two authors contributed equally to this work

Request permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.