Molecular Vision 2011; 17:2292-2296 <http://www.molvis.org/molvis/v17/a249>
Received 30 May 2011 | Accepted 19 August 2011 | Published 25 August 2011

Major intrinsic protein (MIP) polymorphism is associated with age-related cataract in Chinese

Zhou Zhou,1 Binbin Wang,2 Yongfeng Luo,1 Guangkai Zhou,1 Shanshan Hu,1 Han Zhang,3 Xu Ma,2 Yanhua Qi1

The first two authors contributed equally to this work

1Department of Ophthalmology, the Second Affiliated Hospital of Harbin Medical University, Harbin, China; 2Department of Genetics, National Research Institute for Family Planning, Beijing, China; 32007 Clinical Medicine, Harbin Medical University, Harbin, China

Correspondence to: Dr. Yanhua Qi, Department of Ophthalmology, the Second Affiliated Hospital of Harbin Medical University, 246 Xuefu Road, Nan Gang District, Harbin, Heilongjiang, 150086, China; Phone: 86-451-86605851; FAX: 86-451-86605116; email: qi_yanhua@yahoo.com

Abstract

Purpose: Age-related cataract (ARC) is a complex multi-factorial disorder involving several genetic and environmental factors. The major intrinsic protein of lens fiber gene (MIP) encodes the most abundant junctional membrane protein in the mature lens and plays a critical role in maintainace of lens normal structure and internal circulation. To determine the relationship between single nucleotide polymorphisms (SNPs) in MIP and the susceptibility to ARC in a Chinese population, we conducted this case-control study.

Methods: A total of 164 unrelated ARC patients and 132 normal controls were involved in the study. All participants completed full physical and ophthalmic examinations and provided a blood sample for DNA extraction. Seven SNPs (rs2269348, rs61759527, c.-4T>C, rs77163805, rs74641138, rs35033450, and rs36032520) in MIP were amplified by polymerase chain reaction (PCR) and then sequenced. Statistical analysis was performed using SNPstats.

Results: Polymorphisms rs61759527, rs77163805, rs35033450, and rs36032520 were not detected in all 296 subjects. There were no statistical differences in genotype or allele frequency of rs2269348 and rs74641138 between ARC cases and controls. But in c.-4C>T, cataract patients had a higher TC genotype and C allele frequencies (p=0.0018 and p=0.017, respectively) compared to healthy controls. The haplotype CCG of rs2269348, c.-4T>C and rs74641138 also exhibited a significantly higher distribution in cases than controls (OR=8.83, p=0.0024).

Conclusions: Our findings indicate that the genotype TC in polymorphism c.-4T>C and haplotype CCG of rs2269348, c.-4T>C, and rs74641138 in MIP may attach an additional genetic risk factor for ARC in Chinese. This is the first association study about SNPs in MIP and susceptibility to ARC in Chinese population.

Introduction

Age-related cataract (ARC) is the leading cause of low vision and blindness all over the world [1]. As the world’s population ages, cataract-induced visual dysfunction and blindness is on the increase [2]. The development of ARC is complex and multifactorial; the risk factors for cataract including being of the female sex, having lower socioeconomic status, having diabetes mellitus, smoking, and lower body mass index [3-5]. Recently, the contribution of genetic factors in the pathogenesis of ARC was confirmed by some twin studies. In 2000 and 2001, Hammond et al. recognized the heritability for age-related cortical cataract was 53%–58% [6] and approximately 48% for nuclear cataract [7]. Furthermore, several genes were proven to link to ARC, including EPH receptor A2 (EPHA2, chromosome 1p) [8], gap junction protein α8 (GJA8, chromosome 1q) [9], galactose-1-phosphate uridylyltransferase (GALT, chromosome 9p) [10],solute carrier family 16, member 12 (SLC16A12, chromosome 10q) [11], heat shock transcription factor 4 (HSF4, chromosome 16q) [12], galactokinase 1 (GALK1, chromosome 17q) [13], ferritin light polypeptide (FTL, chromosome 19q) [14], xeroderma pigmentosum complementation group (XPD, chromosome 19q), X-ray complementing group1 (XRCC1, chromosome 19q) [15], crystallin αA (CRYAA, chromosome 21q) [16], crystalline βB2 (CRYBB2, chromosome 22q) [17], and the glutathione-S-transferases family (GSTM1, chromosome 1p, GSTP1, chromosome 11q, and GSTT1, chromosome 22q) [18-20].

The lens is a major component of the dynamic vision system. It develops and maintains a tissue structure predicated on the formation and precise organization of specialized cells, the lens fibers. As an avascular tissue, the maintenance of homeostasis and normal structure in lens fiber cells is highly dependent on the channels and adhesion molecules on the cells membrane [21]. The major intrinsic protein of lens fiber gene (MIP), which is also called aquaporin0 (AQP0), encods the most abundant junctional membrane protein in the mature lens. MIP functions not only as a lens-specific water channel but also a cell-to-cell adhesion molecule [22]. Mutations in MIP have been linked to mouse and human congenital cataracts, but there was no report about the association between the single nucleotide polymorphisms (SNPs) in MIP and ARC until now. In this study, we tested seven SNPs in MIP for association with ARC in a Chinese case-control cohort. We substantiated that cataract patients had a higher TC genotype and C allele frequencies in c.-4T>C (p=0.0018 and p=0.017, respectively) compared to healthy controls. The haplotype CCG of rs2269348, c.-4T>C and rs74641138 exhibited a significantly higher distribution (OR=8.83, p=0.0024) in cases than controls as well. Our results suggested there may be an association between the SNP c.-4T>C in MIP and the susceptibility of ARC in the Chinese population.

Methods

Patients and controls

Patients with age-related cataract and normal controls were collected from the northeast of China during our clinical work. Cataract diagnosis was determined according to lens opacities classification system III (LOCSIII) [23]. Patients with secondary cataracts due to trauma, toxins, inflammation, and degenerative ocular diseases were excluded from the study. In addition, patients with diabetes, hypertension, high myopia, glaucoma, any syndrome, and those who were a smoker, an alcoholic, exposed to UVB radiation, or under medication like steroids were also not considered. The control subjects from the same ethnic background were selected during routine medical fitness examination which included ophthalmic examination. From all the patients included in the study, information pertaining to sex, age, age at onset, and family history were collected using a specified form.

Informed consent was obtained from all subjects after explaining the nature and possible consequences of the study. All experiments were approved by the Institutional Review Board of Harbin Medical University (Harbin, China) and conducted according to the principles in the Declaration of Helsinki. To date, all of the patients have had cataract surgery.

SNP selection

We used the SNP database at NCBI and selected SNPs in the 5′UTR, exon 1, and exon 2 of MIP for this study. Six SNPs were chosen for genotyping, including rs2269348, rs61759527 (c.-10C>A), rs77163805 (c.199G>A, p.67V>I), rs74641138 (c.319G>A, p. 107V>I), rs35033450 (c.378C>T) and rs36032520 (c.516C>T). In the process of this study, we found a polymorphism c.-4C>T which had no record in the NCBISNP database but was in 1000 Genomes and we added it into our project.

DNA extraction and genotyping

Five milliliter of venous blood samples were collected in EDTA vaccutainers from ARC patients and control subjects. Genome DNA was extracted from peripheral blood leukocytes with QIAamp DNA Blood Mini Kits (QIAGEN Science, Germantown, MD). The first exon in MIP and its flanking sequences were amplified by polymerase chain reaction (PCR). The primers used in the reaction were 5′-TCT CGG CTC ATC TCC CAG TT-3′ (forward), 5′-GGC AAT AGA GAG ACA GGA CAC-3′ (reverse) and 5′-TGA AGG AGC ACT GTT AGG AGA TG-3′ (forward), 5′- AGA GGG ATA GGG CAG AGT TGA TT −3′ (reverse). The primers were designed with Primer3 according to the reference sequences in the NCBI. We sequenced the PCR products using an ABI3730 Automated Sequencer (PE Biosystems, Foster City, CA), and analyzed the sequencing results using Lasergene SeqMan (DNASTAR, Madison, WI).

Statistical analysis

Allele and genotype frequencies of MIPrs2269348, rs61759527, c.-4T>C, rs77163805, rs74641138, rs35033450, and rs36032520 were compared between patients and controls. We evaluated the frequency of genotypes and alleles in this study using the χ2 test and logistic regression. Paired SNP linkage disequilibrium analysis and haplotype analysis were done as well. The risk for patients was estimated using the odds ratio (OR) and 95% confidence interval (CI). Statistical analysis was performed using SNPstats [24]. A p<0.05 was considered statistically significant.

Results

Characteristics of patients and controls

A total of 296 individuals were included in this study (Table 1). Of these, 164 unrelated individuals were ARC patients with a mean onset age of 66.9±7.3 years. In addition, 132 unaffected individuals with a mean age of 54.2±9.0 years were also recruited. There were significantly more females than males in ARC patients (106:58), and a higher percentage of females overall compared with males (194:102). The patients and controls were matched on gender (χ2 test, p=0.70). The difference of age was adjusted by logistic regression.

Genotype and allele frequencies

The distribution of all the SNPs were consistent with the Hardy–Weinberg equilibrium (HWE) in both case and control groups. Polymorphisms rs61759527, rs77163805, rs35033450, and rs36032520 were not detected in all 296 subjects.

Table 2 showed the distribution of the genotype and allele frequencies of rs2269348, c.-4T>C, and rs74641138 in the study population. Compared with controls, a significant higher distribution of ARC cases was observed in either the TC genotype frequency [OR (95%CI)=6.18 (1.84–20.78), p=0.0018] or C allele frequency (χ2=5.731, p=0.017) in c.-4T>C. There were no statistical differences in the rs2269348 and rs74641138 genotype and allele frequencies between ARC cases and controls.

Paired SNP linkage disequilibrium analysis and haplotype analysis

Paired SNP linkage disequilibrium analysis revealed that rs2269348 was in linkage disequilibrium with c.-4T>C and rs74641138 (D’>0.75; Figure 1). Haplotype case–control analysis denoted that the haplotype CCG with the C allele of rs2269348, the C allele of c.-4T>C, and the G allele of rs74641138 exhibited a significantly higher distribution in cases than controls (OR=8.83,p=0.0024; Table 3).

Discussion

In this study, we investigated the association of MIP polymorphisms with ARC in the Chinese Han population and found a statistically significant difference between cases and controls. There was an association between the gene polymorphisms of MIP and the susceptibility to ARC.

The adult-onset and progressive pattern of inherited cataracts suggested the majority of age-related cataract-associated genes are hereditary cataract candidate genes [25]. So far, several congenital cataract-causing genes are reported to contribute to ARC. Variants in EPHA2 were shown to associate with age-related cortical cataract in both human and mice [8]. Targeted knockout of the mouse Crybb2 can induce age-related cataract in mice [17]. A mutation (F71L) causing defective chaperone-like function in αA-crystallin was reported to be associated with ARC [16]. Shi et al. [12] identified five new HSF4 mutations in 150 ARC patients, finding that these mutations might account for a small fraction of age-related cataracts. Recently, polymorphisms in GJA8 were found to associate with ARC in Chinese people [9]. In SLC16A12, a heterozygous sequence alteration in the 5′UTR (c.-17A>G), which was found in an ARC patient, is in cis with the minor G allele in SNP rs3740030 (c.-42T>G) [11].

MIP (OMIM 154050) was initially a candidate gene for hereditary cataract. Its product, which is also called AQP0, is expressed in lens fiber cell membranes. It is critical for lens transparency and homeostasis as MIP mutations in human and mice have resulted in dominant lens cataract [26,27]. The function of MIP, both as a water channel and an adhesive molecule in the lens fiber, contribute to the narrow intercellular space of the lens fibers that is required for lens transparency and accommodation [8]. Moreover, it is involved in the maintenance of fiber-fiber interactions that regulated the morphology and arrangment of lens fiber cells [28] and could interact with many other lens components including crystallins, lipids, and cytoskeletal proteins [29-31]. With age, significant modifications of AQP0 have already been characterized in human lenses [32]. A continuous age-dependent process of AQP0 truncation was observed and it may contribute to the formation of the lens barrier [33]. As suggested above, AQP0 should play an important role in the process of age-related cataract.

In this study, we found the association of SNP c.-4T>C and haplotype CCG of rs2269348, c.-4T>C and rs74641138 in MIP with ARC, which is presumed to be responsible for the molecular changes leading to cataract formation in aged people. The TC heterozygotes is obviously higher in patients with an odds ratio of 6.18 (95%CI=1.84–20.78, p=0.0018). The haplotype CCG of rs2269348, c.-4T>C, and rs74641138 showed a significant higher frequency in patients than controls with an OR of 8.83 (95%CI=2.19–35.64, p=0.0024). These indicated that the TC genotype and CCG haplotype possibly increased the diease risk. The functional consequences of the c.-4T>C polymorphism is unknown yet. Since it lays in the 5′UTR, 4 bases upstream from the start codon, it may influence the efficiency of protein translation. Previous researches showed that disturbance of the regulation of the translational machinery leads to perturbed cellular metabolism and may tilt the physiologic balance from healthy to diseased states [34]. Linkage disequilibrium analysis revealed that rs2269348, c.-4T>C, and rs74641138 are in linkage disequilibrium. More experiments are needed to explore whether the C allele work alone or with the other alleles in haplotype CCG. We need further studies to find out the precise mechanisms by which the polymorphisms influence the natural history of age-related cataract development.

This is the first association study about MIP gene polymorphisms and age-related cataract in Chinese people. However, the etiology and pathogenesis of ARC is complex and has multiple factors. More research should be done to investigate the relationship between MIP and ARC.

Acknowledgments

The authors are grateful to all patients and normal volunteers for their participation in this investigation. This study was supported by the National Natural Science Foundation of China (Grant No. 30801274) and the Science and Technology Department of Heilongjiang Province Projects of International Cooperation (WB09C104). Professors Xu Ma (genetic@263.net.cn) and Yanhua Qi contributed equally to the research project and can be considered co-corresponding authors.

References

  1. Foster A, Resnikoff S. The impact of Vision 2020 on global blindness. Eye (Lond). 2005; 19:1133-5. [PMID: 16304595]
  2. Brian G, Taylor H. Cataract blindness—challenges for the 21st century. Bull World Health Organ. 2001; 79:249-56. [PMID: 11285671]
  3. Tsai SY, Hsu WM, Cheng CY, Liu JH, Chou P. Epidemiologic study of age-related cataracts among an elderly Chinese population in Shih-Pai, Taiwan. Ophthalmology. 2003; 110:1089-95. [PMID: 12799231]
  4. Nirmalan PK, Krishnadas R, Ramakrishnan R, Thulasiraj RD, Katz J, Tielsch JM, Robin AL. Lens opacities in arural population of southern India: the Aravind Comprehensive Eye Study. Invest Ophthalmol Vis Sci. 2003; 44:4639-43. [PMID: 14578379]
  5. Krishnaiah S, Vilas K, Shamanna BR, Rao GN, Thomas R, Balasubramanian D. Smoking and its association with cataract: results of the Andhra Pradesh eye disease study from India. Invest Ophthalmol Vis Sci. 2005; 46:58-65. [PMID: 15623755]
  6. Hammond CJ, Duncan DD, Snieder H, de Lange M, West SK, Spector TD, Gilbert CE. The heritability of age-related cortical cataract: the twin eye study. Invest Ophthalmol Vis Sci. 2001; 42:601-5. [PMID: 11222516]
  7. Hammond CJ, Sneider H, Spector TD, Gilbert CE. Genetic and environmental factors in age-related nuclear cataracts in monozygotic and dizygotic twins. N Engl J Med. 2000; 342:1786-90. [PMID: 10853001]
  8. Jun G, Guo H, Klein BE, Klein R, Wang JJ, Mitchell P, Miao H, Lee KE, Joshi T, Buck M, Chugha P, Bardenstein D, Klein AP, Bailey-Wilson JE, Gong X, Spector TD, Andrew T, Hammond CJ, Elston RC, Iyengar SK, Wang B. EPHA2 Is Associated with Age-Related Cortical Cataract in Mice and Humans. PLoS Genet. 2009; 5:e1000584 [PMID: 19649315]
  9. Liu Y, Ke M, Yan M, Guo S, Mothobi ME, Chen Q, Zheng F. Association between gap junction protein-alpha 8 polymorphisms and age-related cataract. Mol Biol Rep. 2011; 38:1301-7. [PMID: 20582632]
  10. Karas N, Gobec L, Pfeifer V, Mlinar B, Battelino T, Lukac-Bajalo J. Mutations in galactose-1-phosphate uridyltransferase gene in patients with idiopathic presenile cataract. J Inherit Metab Dis. 2003; 26:699-704. [PMID: 14707519]
  11. Zuercher J, Neidhardt J, Magyar I, Labs S, Moore AT, Tanner F, Waseem NH, Schorderet D, Munier FL, Bhattacharya SS, Berger W, Kloeckener-Gruissem B. Alterations of the 5′ untranslated leader region of SLC16A12 lead to age-related cataract. Invest Ophthalmol Vis Sci. 2010; 51:3354-61.. [PMID: 14707519]
  12. Shi Y, Shi X, Jin Y, Miao A, Bu L, He J, Jiang H, Lu Y, Kong X, Hu L. Mutation screening of HSF4 in 150 age-related cataract patients. Mol Vis. 2008; 14:1850-5. [PMID: 18941546]
  13. Okano Y, Asada M, Fujimoto A, Ohtake A, Murayama K, Hsiao KJ, Choeh K, Yang Y, Cao Q, Reichardt JK, Niihira S, Imamura T, Yamano T. A genetic factor for age-related cataract: identification and characterization of a novel galactokinase variant, “Osaka,” in Asians. Am J Hum Genet. 2001; 68:1036-42. [PMID: 11231902]
  14. Faniello MC, Di Sanzo M, Quaresima B, Nistico A, Fregola A, Grosso M, Cuda G, Costanzo F. Bilateral cataract in a subject carrying a C to A transition in the L ferritin promoter region. Clin Biochem. 2009; 42:911-4. [PMID: 19254706]
  15. Padma G, Mamata M, Reddy KR, Padma T. Polymorphisms in two DNA repair genes (XPD and XRCC1)–association with age related cataracts. Mol Vis. 2011; 17:127-33. [PMID: 21245954]
  16. Bhagyalaxmi SG, Srinivas P, Barton KA, Kumar KR, Vidyavathi M, Petrash JM, Bhanuprakash Reddy G, Padma T. A novel mutation (F71L) in alphaA-crystallin with defective chaperone-like function associated with age-related cataract. Biochim Biophys Acta. 2009; 1792:974-81. [PMID: 19595763]
  17. Zhang J, Li J, Huang C, Xue L, Peng Y, Fu Q, Gao L, Zhang J, Li W. Targeted knockout of the mouse betaB2-crystallin gene (Crybb2) induces age-related cataract. Invest Ophthalmol Vis Sci. 2008; 49:5476-83. [PMID: 18719080]
  18. Juronen E, Tasa G, Veromann S, Parts L, Tiidla A, Pulges R, Panov A, Soovere L, Koka K, Mikelsaar AV. Polymorphic glutathione S-transferases as genetic risk factors for senile cortical cataract in Estonians. Invest Ophthalmol Vis Sci. 2000; 41:2262-7. [PMID: 10892871]
  19. Güven M, Unal M, Sarici A, Ozaydin A, Batar B, Devranoglu K. Glutathione-S-transferase M1 and T1 genetic polymorphisms and the risk of cataract development: a study in the Turkish population. Curr Eye Res. 2007; 32:447-54. [PMID: 17514530]
  20. Zhou J, Hu J, Guan H. The Association between Copy Number Variations in Glutathione S-transferase M1 and T1 and Age-Related Cataract in a Han Chinese Population. Invest Ophthalmol Vis Sci. 2010; 51:3924-8. [PMID: 20335620]
  21. Al-Ghoul KJ, Kirk T, Kuszak AJ, Zoltoski RK, Shiels A, Kuszak JR. Lens structure in MIP-deficient mice. Anat Rec A Discov Mol Cell Evol Biol. 2003; 273:714-30. [PMID: 12845708]
  22. Chepelinsky AB. Structural Function of MIP/Aquaporin0 in the Eye Lens; Genetic Defects Lead to Congenital Inherited Cataracts. Handb Exp Pharmacol. 2009; 190:265-97. [PMID: 19096783]
  23. Chylack LT, , Jr Wolfe JK, Singer DM, Leske MC, Bullimore MA, Bailey IL, Friend J, McCarthy D, Wu SY. The Lens Opacities Classification System III. Arch Ophthalmol. 1993; 111:831-6. [PMID: 8512486]
  24. Solé X, Guino E, Valls J, Iniesta R, Moreno V. SNPStats: a web tool for the analysis of association studies. Bioinformatics. 2006; 22:1928-9. [PMID: 16720584]
  25. Hejtmancik JF, Kantorow M. Molecular genetics of age-related cataract. Exp Eye Res. 2004; 79:3-9. [PMID: 15183095]
  26. 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]
  27. Shiels A, Bassnett S. Mutations in the founder of the MIP gene family underlie cataract development in the mouse. Nat Genet. 1996; 12:212-5. [PMID: 8563764]
  28. Al-Ghoul KJ, Kirk T, Kuszak AJ, Zoltoski RK, Shiels A, Kuszak JR. Lens structure in MIP-deficient mice. Anat Rec A Discov Mol Cell Evol Biol. 2003; 273:714-30. [PMID: 12845708]
  29. Yu XS, Jiang JX. Interaction of major intrinsic protein (aquaporin-0) with fiber connexins in lens development. J Cell Sci. 2004; 117:871-80. [PMID: 14762116]
  30. Fan J, Fariss RN, Purkiss AG, Slingsby C, Sandilands A, Quinlan R, Wistow G, Chepelinsky AB. Specific interaction between lens MIP/Aquaporin-0 and two members of the gamma-crystallin family. Mol Vis. 2005; 11:76-87. [PMID: 15692460]
  31. Gonen T, Cheng Y, Sliz P, Hiroaki Y, Fujiyoshi Y, Harrison SC, Walz T. Lipid-protein interactions in double-layered two-dimensional AQP0 crystals. Nature. 2005; 438:633-8. [PMID: 16319884]
  32. Takemoto L, Takehana M, Horwitz J. Covalent changes in MIP26K during aging of the human lens membrane. Invest Ophthalmol Vis Sci. 1986; 27:443-6. [PMID: 3949474]
  33. Korlimbinis A, Berry Y, Thibault D, Schey K, Truscott R. Protein aging: Truncation of aquaporin 0 in human lens regions is a continuous age-dependent process. Exp Eye Res. 2009; 88:966-73. [PMID: 19135052]
  34. Chatterjee S, Pal JK. Role of 5′- and 3′-untranslated regions of mRNAs in human diseases. Biol Cell. 2009; 101:251-62. [PMID: 19275763]