Molecular Vision 2011; 17:448-455 <http://www.molvis.org/molvis/v17/a51>
Received 10 December 2010 | Accepted 5 February 2011 | Published 10 February 2011

A 556 kb deletion in the downstream region of the PAX6 gene causes familial aniridia and other eye anomalies in a Chinese family

Fang Cheng,1 Wulian Song,1 Yang Kang,2 Shihui Yu,3 Huiping Yuan1

The first two authors contributed equally to this work

1Department of Ophthalmology, the 2nd Affiliated Hospital of Harbin Medical University, Harbin, China; 2Department of Ophthalmology, the 1st Affiliated Hospital of Harbin Medical University, Harbin, China; 3Department of Pathology, Children’s Mercy Hospitals and Clinics and University of Missouri-Kansas City School of Medicine, Kansas City, MO

Correspondence to: Huiping Yuan, Department of Ophthalmology, the 2nd Affiliated Hospital of Harbin Medical University, Harbin, China; Phone: 86-451-86605757; FAX: 86-451-86605757; email: yuanhp@yahoo.com

Abstract

Purpose: The paired box gene 6 (PAX6) on human chromosome 11p13 is an essential transcription factor for eye formation in animals. Mutations in PAX6 can lead to varieties of autosomal-dominant ocular malformations with aniridia as the major clinical signs. Known genetic alterations causing haplo-insufficiency of PAX6 include nonsense mutations, frame-shift mutations, splicing errors, or genomic deletions. The purpose of this study was to identify genetic defects as the underlying cause of familial aniridia in a large Chinese family.

Methods: All exons of PAX6 in the proband were sequenced by the Sanger sequencing technique. The genome of the proband was evaluated by a microarray-based comparative genomic hybridization (aCGH). Quantitative real-time PCR was applied to verify the abnormal aCGH findings in the proband and to test five other family members.

Results: There were no detectable pathogenic mutations in the exons of PAX6 in the proband. The aCGH analysis showed two copies of PAX6 but revealed a 566 kb hemizygous deletion of chromosome 11p13, including four annotated genes doublecortin domain containing 1 (DCDC1), DnaJ homolog subfamily C member 24 (DNAJC24), IMP1 inner mitochondrial membrane(IMMP1L), andelongation factor protein 4 (ELP4) downstream of PAX6. Quantitative real-time PCR verified the deletion in the proband and further identified the deletion in a blind fashion in four affected family members but not in the one with a normal phenotype.

Conclusions: The 566 kb hemizygous deletion of chromosome 11p13 downstream of PAX6 should be the cause of the familial aniridia in this Chinese family, although two copies of PAX6 are intact. aCGH evaluation should be applied if there is a negative result for the mutation detection of PAX6 in patients with aniridia.

Introduction

Aniridia (OMIM 106210) is a congenital eye disorder characterized by the complete or partial absence of the iris and iris hypoplasia. Eight-five percent of individuals with aniridia inherit this disorder as an autosomal-dominant trait, 13% occur as part of the autosomal-dominant WAGR syndrome (Wilms’ tumor, aniridia, genitourinary abnormalities, and mental retardation), and the remaining 2% occur as part of other disorders, including Peters’ anomaly and Gillespie’s syndrome, in either autosomal-dominant or autosomal-recessive inheritances [1-3].

The paired box gene 6 (PAX6) is located on chromosome 11p13 and contains 14 exons encoding a protein (PAX6) with 422 amino acids. PAX6 is a transcriptional factor controlling development of a diversity of organs and tissues (forebrain, pancreas, and ocular tissues, including corneal epithelium, lens, and retina) by recognizing specific DNA sequences of its downstream target genes [4]. Nonsense mutations or deletions of the PAX6 gene primarily cause aniridia due to its haplo-insufficiency [5], while missense mutations of this gene are associated with a diversity of eye abnormalities through gain-of-function of the mutated protein, such as Peters’ anomaly, corneal dystrophy and opacification, congenital cataracts, glaucoma, and foveal hypoplasia [2,6,7]. Whereas mutations or intragenic deletions of PAX6 represent the major causes of aniridia, genomic rearrangements involving the downstream region of PAX6 were identified in patients with aniridia although both copies of PAX6 are intact in these patients [8-13].

In the present study, we report a genomic microdeletion in the downstream region of PAX6 in a large Chinese family with familial aniridia and other eye anomalies, using microarray-based comparative genomic hybridization (aCGH) techniques. To our knowledge this is the first case found in an Asian population and is one of few similar cases with this kind of genetic mechanism [9,11,12].

Methods

Patients

A five-generation Chinese family with familial aniridia was recruited from Heilongjiang province, northeastern China. There were 15 affected individuals in this family from which five affected members (II:2, III:4, III:6, III:13, and IV:11) and one unaffected individual (IV:4) participated in this study (Figure 1). Ocular tests for these six members included visual acuity of naked eyes and corrected visual acuity, which were recorded using the logarithm of the minimum angle of resolution E chart (Precision Vision, Villa Park, IL), slit-lamp biomicroscopy, measurement of intraocular pressure by applanation tonometry, and gonioscopic evaluation of the anterior chamber angle. Systemic evaluation was performed in the six subjects in this study to exclude WAGR syndrome (Wilms’ tumor, aniridia, genitourinary abnormalities, and mental retardation), Axenfeld–Rieger syndrome, iridocorneal endothelial syndromes, sclerocornea, and Peter’s anomaly. Informed consents were obtained from the six individuals in this study. Five affected patients and one family member with normal phenotype were further investigated using molecular techniques. The research protocol for this study was approved by the Harbin Medical University Ethics Committee (Harbin, China).

Experiments

All exons of PAX6 in the proband (III:4) were amplified and sequenced with the ABI BigDye Terminator Cycle Sequencing kit v3.1, (ABI Applied Biosystems, Foster City, CA) according to the described protocol [14]. The superimposed mutant PCR products were subcloned into pGEM-T vector (Promega, Madison, WI) and sequenced to identify the mutation. The genomic DNA from the proband was analyzed by aCGH using the Agilent Human Genome Microarray Kit 244K (Agilent Technologies, Santa Clara, CA) based on the published procedures [15]. The digested test DNA and reference DNA were labeled with cyanine (Cy)3-deoxyuridine triphosphate (dUTP) or Cy5-dUTP, respectively. Following purification, the appropriate Cy3-labeled test DNA and Cy5-labeled reference DNA were mixed together and combined with 2× Hybridization buffer (Agilent), 10× blocking agent (Agilent), and Human Cot-1 DNA (Invitrogen, Carlsbad, CA). The hybridization mixture was slowly dispensed to a microarray chip and assembled with an Agilent SureHyb chamber. After washing, all microarray slides were scanned on an Agilent Microarray Scanner G2565BA with 5-μm resolution. Captured images were transformed to data with Feature Extraction Software, version 9.5 (Agilent), and then imported into Agilent CGH Analytics 3.2.5 software for analysis. Quantitative real-time PCR (qPCR) was performed to verify the abnormal aCGH findings in the proband and to test other family members, according to published guidelines [16]. Individual optimization qPCRs were performed in a 20 μl volume including 25 ng template DNA, 200 nM of each primer, and 1× Premix-Choice with ROX reference dye in an initial denaturation of 93 °C for 10min, followed by 35 cycles of 93 °C for 1 min and 61 °C for 1 min, and 72 °C for 1 min. The three sets of primers for qPCR analysis (Table 1) include primers targeted to the elongation factor protein 4 (ELP4) gene within the deleted genomic region, primers within the PAX6 gene, and reference primers within the Oral-facial-digital syndrome 1 protein (OFD1) gene on the human X chromosome. Both ELP4 and PAX6 were assayed as test genes compared with OFD1, and both ELP4 and PAX6 were assayed as test/controls.

Results

Clinical findings

The clinical findings are summarized in Table 2. The proband (III:4), was diagnosed with bilateral iris coloboma, optic atrophy of the left eye, right corneal opacity, and strabismus. Because of glaucoma, the proband had undergone cyclocryotherapy on the left eye in 2002. Other recurrent symptoms in the other affected participants include optic atrophy, retinitis pigmentosa, cataract, subluxation of lens, and dysplasia of trabecular meshwork, fovea, and optic nerve (Table 2). Poor vision is due to foveal or optic nerve hypoplasia, cataract, glaucoma and amblyopia. No visible autistic problems or intellectual disabilities were identified in these affected individuals.

Sequencing results for PAX6

All exons of PAX6 from the six participants of the family were amplified and sequenced using standard techniques. No intragenic point mutation or deletion could be identified (data not shown).

Microarray-based comparative genomic hybridization findings in the proband

The aCGH analysis detected 35 copy number variations (CNVs) in the proband’s genome (Table 3). One of them is a 566 kb hemizygous deletion of chromosome 11p13 (chr:31,074,403–31,640,263) with approximately a 123 kb distance from the 3′ end of PAX6 (chr11:31,762,916–31,796,085) according to HG18 (NCBI 36, March 2006; Figure 2. This deletion contains four annotated genes: doublecortin domain containing 1 (DCDC1), DnaJ homolog subfamily C member 24 (DNAJC24), IMP1 inner mitochondrial membrane (IMMP1L), and elongation factor protein 4 (ELP4). The remaining 34 CNVs were considered to be benign because these CNVs were identified in healthy individuals documented in the Database of Genomic Variation or in recent publications [17,18].

Quantitative real-time PCR assays in the proband and other family members

The 566 kb hemizygous deletion of chromosome 11p13 in the proband was verified by qPCR using primers targeted to ELP within the deleted region, which showed one threshold cycle difference between the patient and reference DNA samples for the amplification of the test gene (Figure 3). Two copies of PAX6 were confirmed by qPCR, consistent with the aCGH results in this individual (data not shown). We used the qPCR method for testing the other five family members in a blind fashion, four affected family members with eye anomalies, and one healthy individual. All four patients were confirmed to carry the deletion, including ELP, but not PAX6. The phenotypic normal individual showed two copies for both ELP and PAX6 (data not shown).

Discussion

The major finding in this Chinese family with familial aniridia is the presence of a 566 kb heterozygous deletion containing four annotated genes: DCDC1, DNAJC24, IMMP1L, and ELP4. The proximal breakpoint of this deletion is approximately 123 kb from the 3′ end of PAX6. We postulate that this 566 kb heterozygous deletion is the underlying cause of the familial aniridia and acts by disrupting the transcription in one of the two PAX6 alleles, even though the two copies of PAX6 were intact in all individuals investigated in this study.

Our finding provides further evidence of the existence of the remote 3′ regulatory elements in the downstream region of PAX6 controlling the expression of this gene, if disrupted, leading to aniridia and other eye anomalies. To our knowledge, this is the first case found in Asian patients with aniridia and one of few similar cases with this kind of genetic mechanism.

Our postulation is based on following evidence: 1) The heterozygous deletion segregated with aniridia in the five affected individuals but not in the unaffected individual, while the exons and splicing regions of PAX6 are apparently free of mutations. 2) Several publications reported similar observations in patients with aniridia, but the chromosomal breakpoint from the 3′ end of PAX6 and the fragment of deletion were different. For example, two aniridia pedigrees have been characterized in which the disease segregates with chromosomal re-arrangements that involve 11p13 but do not disrupt the PAX6 gene since the chromosomal breakpoint is at least 85 kb away from the 3′ end of PAX6 [8]. Two aniridia pedigrees have also been described in which deletion in the ELP4 gene region, not involving PAX6, was present in all subjects with aniridia but not in the investigated normal relatives [12]. A 1.3 Mb deletion (including seven annotated genes: metallophosphoesterase domain containing 2 (MPPED2), doublecortin domain containing 5 (DCDC5), DCDC1, DNAJC24, IMMP1L, zinc finger CLS domain containing 3 (DPH4), and ELP4) has been characterized that starts 35 kb from the 3′ end of PAX6 in a patient with aniridia, autism, and mental retardation [11]. Recently a ~406 kb heterozygous genomic deletion containing four annotated genes (DCDC1, DNAJC24, IMMP1L, and ELP4) was found in patients with aniridia [13]; apparently the gene contents in this deletion are the same as the 566 kb heterozygous deletion in the family we report here. 3) Functional studies in both human cells and animal models confirmed the existence of remote 3′ regulatory elements in the downstream region of PAX6. Deletions in this region have been shown to abolish PAX6 expression and cause aniridia and other eye anomalies due to loss of enhancers and a downstream regulatory region [9,19,20].

Little information is known about the four annotated genes within the 566 kb deletion in this study. It is also unknown whether these genes are involved in any of the phenotypic features found in these individuals carrying the deletion. The DCDC1 gene encodes a member of the doublecortin family, which is highly expressed in testis and fetal brain [21]. DNAJC24 is one of several enzymes involved in synthesis of diphthamide, which is a unique posttranslationally modified histidine found only in translation elongation factor 2 [22]. IMMP1L encodes a peptidase similar to mitochondrial inner membrane peptidase (IMP1), which is one of the catalytic subunits of the IMP complex proteolytically removing the mitochondrial targeting presequence of nuclear-encoded proteins [23]. ELP4 encodes a component of the six subunit elongator complex, a histone acetyltransferase complex that associates directly with RNA polymerase II during transcriptional elongation. Two recent reports indicate that ELP4 is possibly associated with the centrotemporal sharp wave electroencephalogram (EEG) trait in rolandic epilepsy and speech sound disorder [24,25].

Submicroscopic copy number variations may play a role in human diseases either by loss of gene expression regulatory elements or by disrupting coding sequences. As well as the point mutations in PAX6 exons, copy number variation should be investigated in the flanking regions of PAX6. We suggest that patients, such as the subjects reported here, should be investigated using high resolution aCGH techniques in a clinical setting if sequencing analyses for PAX6 in patients with aniridia is negative.

Acknowledgments

This study was supported by the innovation funding of the graduated students of Harbin Medical University (No. Hcxb2008017). We are grateful to all the patients, their families and all the people who ever provided help to this study.

References

  1. Kokotas H, Petersen MB. Clinical and molecular aspects of aniridia. Clin Genet. 2010; 77:409-20. [PMID: 20132240]
  2. Hanson IM, Fletcher JM, Jordan T, Brown A, Taylor D, Adams RJ, Punnett HH, van Heyningen V. Mutations at the PAX6 locus are found in heterogeneous anterior segment malformations including Peters' anomaly. Nat Genet. 1994; 6:168-73. [PMID: 8162071]
  3. François J, Lentini F, de Rouck F. Gillespie's syndrome (incomplete aniridia, cerebellar ataxia and oligophrenia). Ophthalmic Paediatr Genet. 1984; 4:29-32. [PMID: 6544390]
  4. Azuma N, Yamaguchi Y, Handa H, Hayakawa M, Kanai A, Yamada M. Missense mutation in the alternative splice region of the PAX6 gene in eye anomalies. Am J Hum Genet. 1999; 65:656-63. [PMID: 10441571]
  5. Martha A, Ferrell RE, Mintz-Hittner H, Lyons LA, Saunders GF. Paired box mutations in familial and sporadic aniridia predicts truncated aniridia proteins. Am J Hum Genet. 1994; 54:801-11. [PMID: 7909985]
  6. Azuma N, Hotta Y, Tanaka H, Yamada M. Missense mutations in the PAX6 gene in aniridia. Invest Ophthalmol Vis Sci. 1998; 39:2524-8. [PMID: 9856761]
  7. Xu HE, Rould MA, Xu W, Epstein JA, Maas RL, Pabo CO. Crystal structure of the human Pax6 paired domain-DNA complex reveals specific roles for the linker region and C-terminal subdomain in DNA binding. Genes Dev. 1999; 13:1263-75. [PMID: 10346815]
  8. Fantes J, Redeker B, Breen M, Boyle S, Brown J, Fletcher J, Jones S, Bickmore W, Fukushima Y, Mannens M, Danes S, Heyningen V, Hanson I. Aniridia-associated cytogenetic rearrangements suggest that a position effect may cause the mutant phenotype. Hum Mol Genet. 1995; 4:415-22. [PMID: 7795596]
  9. Lauderdale JD, Wilensky JS, Oliver ER, Walton DS, Glaser T. 3′ deletions cause aniridia by preventing PAX6 gene expression. Proc Natl Acad Sci USA. 2000; 97:13755-9. [PMID: 11087823]
  10. Crolla JA, van Heyningen V. Frequent chromosome aberrations revealed by molecular cytogenetic studies in patients with aniridia. Am J Hum Genet. 2002; 71:1138-49. [PMID: 12386836]
  11. Davis LK, Meyer KJ, Rudd DS, Librant AL, Epping EA, Sheffield VC, Wassink TH. Pax6 3′ deletion results in aniridia, autism and mental retardation. Hum Genet. 2008; 123:371-8. [PMID: 18322702]
  12. D'Elia AV, Pellizzari L, Fabbro D, Pianta A, Divizia MT, Rinaldi R, Grammatico B, Grammatico P, Arduino C, Damante G. A deletion 3′ to the PAX6 gene in familial aniridia cases. Mol Vis. 2007; 13:1245-50. [PMID: 17679951]
  13. Bayrakli F, Guney I, Bayri Y, Ercan-Sencicek AG, Ceyhan D, Cankaya T, Mason C, Bilguvar K, Bayrakli S, Mane SM, State MW, Gunel M. A novel heterozygous deletion within the 3′ region of the PAX6 gene causing isolated aniridia in a large family group. J Clin Neurosci. 2009; 16:1610-4. [PMID: 19793656]
  14. Yuan H, Kang Y, Shao Z, Li Y, Yang G, Xu N. Two novel PAX6 mutations identified in northeastern Chinese patients with aniridia. Mol Vis. 2007; 13:1555-61. [PMID: 17893655]
  15. Yu S, Bittel DC, Kibiryeva N, Zwick DL, Cooley LD. Validation of the Agilent 244K oligonucleotide array-based comparative genomic hybridization platform for clinical cytogenetic diagnosis. Am J Clin Pathol. 2009; 132:349-60. [PMID: 19687311]
  16. Yu S, Kielt M, Stegner AL, Kibiryeva N, Bittel DC, Cooley LD. Quantitative real-time polymerase chain reaction for the verification of genomic imbalances detected by microarray-based comparative genomic hybridization. Genet Test Mol Biomarkers. 2009; 13:751-60. [PMID: 20001581]
  17. Shaikh TH, Gai X, Perin JC, Glessner JT, Xie H, Murphy K, O'Hara R, Casalunovo T, Conlin LK, D'Arcy M, Frackelton EC, Geiger EA, Haldeman-Englert C, Imielinski M, Kim CE, Medne L, Annaiah K, Bradfield JP, Dabaghyan E, Eckert A, Onyiah CC, Ostapenko S, Otieno FG, Santa E, Shaner JL, Skraban R, Smith RM, Elia J, Goldmuntz E, Spinner NB, Zackai EH, Chiavacci RM, Grundmeier R, Rappaport EF, Grant SF, White PS, Hakonarson H. High-resolution mapping and analysis of copy number variations in the human genome: a data resource for clinical and research applications. Genome Res. 2009; 19:1682-90. [PMID: 19592680]
  18. Itsara A, Cooper GM, Baker C, Girirajan S, Li J, Absher D, Krauss RM, Myers RM, Ridker PM, Chasman DI, Mefford H, Ying P, Nickerson DA, Eichler EE. Population analysis of large copy number variants and hotspots of human genetic disease. Am J Hum Genet. 2009; 84:148-61. [PMID: 19166990]
  19. Griffin C, Kleinjan DA, Doe B, van Heyningen V. New 3′ elements control Pax6 expression in the developing pretectum, neural retina and olfactory region. Mech Dev. 2002; 112:89-100. [PMID: 11850181]
  20. Kleinjan DA, Seawright A, Mella S, Carr CB, Tyas DA, Simpson TI, Mason JO, Price DJ, van Heyningen V. Long-range downstream enhancers are essential for Pax6 expression. Dev Biol. 2006; 299:563-81. [PMID: 17014839]
  21. Zeng L, Gu S, Li Y, Zhao E, Xu J, Ye X, Wu Q, Wang L, Xie Y, Mao Y. Identification of a novel human doublecortin-domain-containing gene (DCDC1) expressed mainly in testis. J Hum Genet. 2003; 48:393-6. [PMID: 12820024]
  22. Liu S, Milne GT, Kuremsky JG, Fink GR, Leppla SH. Identification of the proteins required for biosynthesis of diphthamide, the target of bacterial ADP-ribosylating toxins on translation elongation factor 2. Mol Cell Biol. 2004; 24:9487-97. [PMID: 15485916]
  23. Burri L, Strahm Y, Hawkins CJ, Gentle IE, Puryer MA, Verhagen A, Callus B, Vaux D, Lithgow T. Mature DIABLO/Smac is produced by the IMP protease complex on the mitochondrial inner membrane. Mol Biol Cell. 2005; 16:2926-33. [PMID: 15814844]
  24. Strug LJ, Clarke T, Chiang T, Chien M, Baskurt Z, Li W, Dorfman R, Bali B, Wirrell E, Kugler SL, Mandelbaum DE, Wolf SM, McGoldrick P, Hardison H, Novotny EJ, Ju J, Greenberg DA, Russo JJ, Pal DK. Centrotemporal sharp wave EEG trait in rolandic epilepsy maps to Elongator Protein Complex 4 (ELP4). Eur J Hum Genet. 2009; 17:1171-81. [PMID: 19172991]
  25. Pal DK, Li W, Clarke T, Lieberman P, Strug LJ. Pleiotropic effects of the 11p13 locus on developmental verbal dyspraxia and EEG centrotemporal sharp waves. Genes Brain Behav. 2010; 9:1004-12. [PMID: 20825490]