- Research article
- Open Access
Screening of potential target genes for cataract by analyzing mRNA expression profile of mouse Hsf4-null lens
BMC Ophthalmology volume 15, Article number: 76 (2015)
Hsf4 is closely related to the development of cataract. However, the molecular mechanisms remain unknown. This study aimed to explore the molecular mechanisms that how Hsf4 mutations influence development of lens and thus lead to cataract in mouse.
The mRNA expression profile of mouse tissue samples from Hsf4-null and wile-type lenses was downloaded from Gene Expression Omnibus database. Then the LIMMA package was used to screen differentially expressed genes (DEGs) and DAVID was applied to identify the significantly enriched Gene Ontology (GO) categories for DEGs. Furthermore, the protein-protein interaction (PPI) network of DEGs was constructed using Cytoscape and the key modules were selected from the PPI network based on the MCODE analysis.
A total of 216 DEGs were screened, including 51 up- and 165 down-regulated genes. Meanwhile, nine GO terms were obtained, and DEGs such as SGK1, CRY2 and REV1 were enriched in response to DNA damage stimulus. Furthermore, 89 DEGs and 99 gene pairs were mapped into the PPI network and Ubc was the hob node. Two key modules, which contained the genes (e.g. Ubc, Egr1, Ptgs2, Hmox1, Cd44, Btg2, Cyr61 and Fos) were related to response to DNA damage stimulus.
The deletion of Hsf4 affects the expression of many genes, such as Ubc, Ptgs2, Egr1 and Fos. These genes may be involved in the development of cataract and could be used as therapeutic targets for cataract.
Cataract is a visible opacity in the lens substance which leads to a decrease in vision. The lens is a critical refractive element of the eye which, with the cornea, focuses images of the visual world onto the retina . Previous study has indicated that both the structure and stability of lens crystallins and maintenance of strong cellular homeostatic systems are required for sustaining normal function of lens . Age and genetic component are main factors to influence the development of lens and hence cause cataract . Nowadays, cataract remains the leading cause of blindness in the world, especially in developing countries .
Heat shock transcription factor 4 (Hsf4), a member of Hsf family, is the common gene linked to cataractogenesis and it has been regarded as a causative gene for congenital cataract . Hsf4 is expressed exclusively in the ocular lens and acts a key role in the lens formation and differentiation . Besides, Hsf4 regulates DLAD expression and promotes lens de-nucleation , and it is involved in the negative regulation of DNA binding activity . Furthermore, Cui et al. have found that Hsf4 promotes DNA damage repair through the regulation of Rad51 expression . It has been also reported that Hsf4 mutations may also be associated with age-related cataract  and mutations in the DNA binding domain (A20D, I87V, L115P, R120C and R74H) of Hsf4 cause autosomal dominant cataract [7, 10, 11]. Meanwhile, several studies have reported that mutations located in the other domain of Hsf4 contribute to the autosomal recessive cataract [12–14]. Also, removal of Hsf4 has been found to lead to cataract development in the Hsf4-null (Hsf4-/-) mice through reduction of γS-crystallin and Bfsp expression .
Gene microarray analysis provides a powerful method for rapid, comprehensive, and quantitative analysis of gene expression profiles of normal/disease states and developmental processes . Thus, the expression levels of thousands of genes can be quantified simultaneously with this technology . Using gene microarray analysis, He et al.  have found that Brg1, Hsf4 and Pax6 exert their functions through commonly regulating other genes. However, the molecular mechanisms of Hsf4 are still not fully understood. To further investigate the molecular mechanisms that how mouse Hsf4 (mHsf4) mutations influence lens development and lead to cataract, the mRNA expression profile of mHsf4-null mutation and wide-type lens was downloaded from Gene Expression Omnibus (GEO) database deposited by He et al. . Then the differentially expressed genes (DEGs) were identified and used to construct the protein-protein interaction (PPI) network. Furthermore, the significantly enriched functions and important modules were screened and analyzed.
Microarray data and data preprocessing
The mRNA expression profile of GSE22362  was obtained from GEO (available at http://www.ncbi.nlm.nih.gov/geo/) database . The total microarray contains six chips of mouse tissue samples from Hsf4-null and wild-type lenses, which were described as a previous study . The expression profile was analyzed by the platform of GPL8321 [Mouse430A_2] Affymetrix Mouse Genome 430A 2.0 Array (Affymetrix, Inc., Santa Clara, CA, USA). Raw data were preprocessed via background correction, quantile normalization and probe summarization using Affy software package  of R. Then the probe-level data in CEL files were converted into the mRNA expression values. In the case, if there was more than one probe in a single gene, the average expression values of all probes for a given gene were defined as the mRNA expression value. Meanwhile, when several mRNAs were mapped by one probe, this probe was thought to lack specificity, and was removed from the analysis.
Screening of DEGs
The wild-type samples were classed as the controls ant the normalized data were analyzed using LIMMA (Linear Models for Microarray Data, available at http://www.bioconductor.org/packages/release/bioc/html/limma.html, V 3.22.1) package . Then the p-value was adjusted into FDR (false discovery rate)  by Bonferroni method  in multtest package. The mRNAs with the cutoff criteria of |log2fold change (FC)| >1 and FDR <0.05 were considered to be DEGs. Furthermore, to explore whether the mRNAs were samples-specific, Pheatmap package (available at http://cran.r-project.org/web/packages/pheatmap/index.html, V 0.7.7)  in R was used to perform hierarchical clustering by comparing the value of each mRNA in six samples.
Functional enrichment analysis
The Database for Annotation, Visualization and Integrated Discovery (DAVID)  is the most common tool to analysis the functional enrichment of genes. To identify the functions of DEGs, the DAVID was used to identify the significantly enriched GO (Gene Ontology) categories. The p-value <0.05 was selected as cutoff criterion.
Construction of PPI network
The Search Tool for the Retrieval of Interacting Genes (STRING, available at http://www.string-db.org/) database is a useful tool that provides lots of experimental and predicted information of proteins [26, 27]. In order to research the relationship between genes, the DEGs were scanned by the STRING and the PPI pairs were selected with the cutoff criterion of combined score >0.4. Then the PPI network was visualized using Cytoscape (available at http://www.cytoscape.org/) .
Screening of modules
Proteins encoded by genes in the same module may perform the same or similar functions. To further explore the functions of proteins, the network modules were obtained from the PPI network based on the MCODE analysis . Default parameters (Degree Cutoff: 2, K-Core: 2) were set as the threshold for modules screening.
A set of 216 DEGs were identified in the Hsf4-null samples compared with wide-type samples, including 51 up- and 165 down-regulated DEGs. Besides, the hierarchical clustering analysis indicated that the DEGs in Hsf4-null samples were distinguished from that in wide-type controls (Fig. 1).
Functional enrichment analysis
In total, nine GO biological processes were obtained (Fig. 2). Among these functions, DEGs such as BTG2, HMOX1 and REV1 were significantly enriched in response to DNA damage stimulus (p = 5.01E-05); DEGs such as FOS, EGR1 and MSX1 were distinctly enriched in response to protein stimulus (p = 1.16E-04); DEGs such as PTGS2, ACVR2A and ALOX15 were markedly enriched in skeletal system development (p = 9.74E-04) (Table 1).
PPI network construction
Among the 261 DGEs, 99 gene pairs were identified with the combined score >0.4 and 89 DEGs (25 up- and 64 down-regulated) were mapped into the PPI network (Fig. 3). Several nodes had higher connectivity degrees in the PPI network, such as Ubc (degree = 18), Ptgs2 (degree = 10), Fos (degree = 9), Cd44 (degree =8), Hsph1 (degree = 6), and Gnajb1 (degree = 6) (Table 2).
Screening of modules
The GO enrichment analysis indicated that response to DNA damage stimulus was the most significant function. To further understand the relationships between DNA damage and cataract, MCODE was used to identify the functional modules of the genes that were related to DNA damage. As a result, two modules were screened. Module one contained Ubc, Egr1, Ptgs2, Hmox1 and Cd44 and module two contained Btg2, Cyr61 and Fos (Fig. 4).
Cataract is the opacification of the eye lens, and is the leading cause of blindness worldwide . Cataractogenesis has multiple causes and is often associated with an abnormality of the lens microarchitecture . Hsf4 is prominently expressed in the lens compared with in other tissues and closely related to the development of cataract [5, 32]. In the present study, we aimed to extend our understanding of the influence of lens development caused by Hsf4 mutations. Results revealed that expression levels of 216 genes were altered in mHsf4-null lens compared with wide-type controls. Functional enrichment results showed that response to DNA damage stimulus was the most significant function in mHsf4-null lens. In addition, we identified two key modules correlated with response to DNA damage stimulus from the PPI network.
Cataracts may be caused primarily by the DNA damage, such as oxidized purines  and DNA single strand breaks . Study has found that oxidative DNA damage is significantly high in the lens epithelial cells (LECs) of cataract patients [35, 36]. Besides, loss-of-function mutations in TBC1D20 cause cataracts in blind sterile mice . The DNA damage response (DDR) is a signal transduction pathway that senses DNA damage and sets a response to protect the cell and moderate the threat to the organism [38, 39]. Mice with the knock-down of HSF4 have cataract because of an increased proliferation of LECs in the lens as well as an abnormal lens fiber cell differentiation . Two modules which were closely related to DDR, were isolated from the PPI network. Ubc (ubiquitin C/polyubiquitin-C), the module-related gene, was also showed to be the hub node in the PPI network. Ubc is one of the sources of ubiquitin during cell proliferation and stress that cannot be compensated by other ubiquitin . Ubiquitin is a normal component in the lens and a ubiquitin-dependent proteolytic system exists in lens , which consists of fiber cells that differentiate from epithelial cells and undergo programmed organelle degradation during terminal differentiation . Furthermore, ubiquitination is a reversible post-translational modification of cellular proteins and is considered to play key roles in the regulation of varieties of cellular processes, such as protein degradation, cell-cycle regulation, DNA repair, apoptosis and signal transduction . Moreover, the ubiquitin proteasome system is found to be essential to cell proliferation of the lens epithelium and required for differentiation of lens fiber cells in zebrafish . The overexpression of ubiquitin affects ubiquitin proteasome system and thus disorders cell proliferation and differentiation of lens. Thus our results suggested that the lack of Hsf4 up-regulated the expression of Ubc, which might be involved in the development of cataract by regulating the cellular processes of lens.
Ptgs2 (prostaglandin-endoperoxide synthase 2) is also known as cyclooxygenase-2 or COX-2, which is involved in the conversion of arachidonic acid prostaglandin H2. In human fibroblasts, Ptgs2 has been showed to interact with Cav1 (Caveolin 1) , which is the main component of the caveolae plasma membranes. Caveolae is cholesterol-rich lipid rafts that are likely to play important roles in lens . What is more, Cav1 was found to participate in repair of DNA damage through regulating the important molecules involved in maintaining genomic integrity . Besides, redundant Cav1 has been reported to play a role in age-dependent hyporesponsiveness to growth factors in vitro and may act as an indicator of wound-healing capacity in aged human corneal epithelium . Therefore, Ptgs2 plays a role in the response to DNA damage and may be related to the repair of DNA damage in lens through the interaction with Cav1.
Egr1 (early growth response 1) belongs to EGR family of zinc finger proteins and functions as a transcriptional regulator. It has been reported that the mRNA expression of Egr1 can be used as a marker for the direction of mammalian ocular growth . In addition, Fos (FBJ murine osteosarcoma viral oncogene homolog), also named c-FOS, can be induced by a variety of extracellular stimuli  and interact with Jun (jun proto-oncogene, c-JUN) to form the transcription factor AP-1 (activating protein 1) , which regulates cell adaptation to environmental changes . Furthermore, Fos and Jun are differentially regulated during terminal differentiation of lens fiber cells . Thus, Egr1 and Fos may be involved in the cell cycle and apoptosis of lens.
However, there were some limitations in this study. For example, there were no experiments to confirm our predictions. The number of samples were also small. Considering these issues, the experimental studies will be subsequently conducted later using more samples.
In conclusion, the deficiency of Hsf4 affect the expression of a set of genes, especially Ubc, Ptgs2, Egr1 and Fos, which are closely related to the response to DNA damage stimulus. These genes may be participated in the development of cataract by influencing the cellular activities of lens and could be used as therapeutic targets for cataract if they were validated by the further experiments which would be conducted later.
Differentially expressed genes
Heat shock transcription factor 4
Early growth response 1
Michael R, Bron A. The ageing lens and cataract: a model of normal and pathological ageing. Philos Trans R Soc Lond B Biol Sci. 2011;366:1278–92.
Shiels A, Hejtmancik J. Genetics of human cataract. Clin Genet. 2013;84:120–7.
Courtney P. The National Cataract Surgery Survey: I. Method and descriptive features. Eye. 1992;6:487–92.
Qureshi N, Ahmed T, Ahmed T. Opacities in optical media to cause diminished vision. Pak J Surg. 2014;30:63–6.
Cui X, Wang L, Zhang J, Du R, Liao S, Li D, et al. HSF4 regulates DLAD expression and promotes lens de-nucleation. Biochim Biophys Acta. 1832;2013:1167–72.
Fujimoto M, Izu H, Seki K, Fukuda K, Nishida T, Yamada S, et al. HSF4 is required for normal cell growth and differentiation during mouse lens development. EMBO J. 2004;23:4297–306.
Nakai A, Tanabe M, Kawazoe Y, Inazawa J, Morimoto RI, Nagata K. HSF4, a new member of the human heat shock factor family which lacks properties of a transcriptional activator. Mol Cell Biol. 1997;17:469–81.
Cui X, Zhang J, Du R, Wang L, Archacki S, Zhang Y, et al. HSF4 is involved in DNA damage repair through regulation of Rad51. Biochim Biophys Acta. 1822;2012:1308–15.
Shi Y, Shi X, Jin Y, Miao A, Bu L, He J, et al. Mutation screening of HSF4 in 150 age-related cataract patients. Mol Vis. 2008;14:1850–5.
Bu L, Jin Y, Shi Y, Chu R, Ban A, Eiberg H, et al. Mutant DNA-binding domain of HSF4 is associated with autosomal dominant lamellar and Marner cataract. Nat Genet. 2002;31:276–8.
Ke T, Wang QK, Ji B, Wang X, Liu P, Zhang X, et al. Novel HSF4 mutation causes congenital total white cataract in a Chinese family. Am J Ophthalmol. 2006;142:298–303.
Smaoui N, Beltaief O, BenHamed S, M’Rad R, Maazoul F, Ouertani A, et al. A homozygous splice mutation in the HSF4 gene is associated with an autosomal recessive congenital cataract. Invest Ophthalmol Vis Sci. 2004;45:2716–21.
Forshew T, Johnson CA, Khaliq S, Pasha S, Willis C, Abbasi R, et al. Locus heterogeneity in autosomal recessive congenital cataracts: linkage to 9q and germline HSF4 mutations. Hum Genet. 2005;117:452–9.
Sajjad N, Goebel I, Kakar N, Cheema AM, Kubisch C, Ahmad J. A novel HSF4 gene mutation (p.R405X) causing autosomal recessive congenital cataracts in a large consanguineous family from Pakistan. BMC Med Genet. 2008;9:1471–2350.
Shi X, Cui B, Wang Z, Weng L, Xu Z, Ma J, et al. Removal of Hsf4 leads to cataract development in mice through down-regulation of γS-crystallin and Bfsp expression. BMC Mol Biol. 2009;10:10.
Lovén J, Orlando DA, Sigova AA, Lin CY, Rahl PB, Burge CB, et al. Revisiting global gene expression analysis. Cell. 2012;151:476–82.
Ueda M, Ota J, Yamashita Y, Choi YL, Ohki R, Wada T, et al. DNA microarray analysis of stage progression mechanism in myelodysplastic syndrome. Br J Haematol. 2003;123:288–96.
He S, Pirity MK, Wang W-L, Wolf L, Chauhan BK, Cveklova K, et al. Chromatin remodeling enzyme Brg1 is required for mouse lens fiber cell terminal differentiation and its denucleation. Epigenetics Chromatin. 2010;3:21.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2013;41:D991–D5.
Gautier L, Cope L, Bolstad BM, Irizarry RA. affy—analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004;20:307–15.
Smyth GK. Limma: linear models for microarray data. Bioinformatics and computational biology solutions using R and Bioconductor. New York: Springer; 2005. p. 397–420.
Benjamini Y. Discovering the false discovery rate. J Roy Stat Soc B. 2010;72:405–16.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B. 1995;57:289–300.
Wang L, Cao C, Ma Q, Zeng Q, Wang H, Cheng Z, et al. RNA-seq analyses of multiple meristems of soybean: novel and alternative transcripts, evolutionary and functional implications. BMC Plant Biol. 2014;14:169.
da Huang W, Sherman BT, Tan Q, Collins JR, Alvord WG, Roayaei J, et al. The DAVID Gene Functional Classification Tool: a novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol. 2007;8:R183.
Von Mering C, Huynen M, Jaeggi D, Schmidt S, Bork P, Snel B. STRING: a database of predicted functional associations between proteins. Nucleic Acids Res. 2003;31:258–61.
Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, et al. STRING v9. 1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013;41:D808–D15.
Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011;27:431–2.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4:2.
Watkins R. Foundations of a solution to cataract blindness. Clin Exp Optom. 2002;85:59–60.
Enoki Y, Mukoda Y, Furutani C, Sakurai H. DNA-binding and transcriptional activities of human HSF4 containing mutations that associate with congenital and age-related cataracts. Biochim Biophys Acta. 2010;1802:749–53.
Lv H, Huang C, Zhang J, Liu Z, Zhang Z, Xu H, et al. A Novel HSF4 Gene Mutation Causes Autosomal-Dominant Cataracts in a Chinese Family. G3 (Bethesda). 2014;4:823–8.
Øsnes-Ringen O, Azqueta AO, Moe MC, Zetterström C, Røger M, Nicolaissen B, et al. DNA damage in lens epithelium of cataract patients in vivo and ex vivo. Acta Ophthalmol. 2013;91:652–6.
Kleiman NJ, Spector A. DNA single strand breaks in human lens epithelial cells from patients with cataract. Curr Eye Res. 1993;12:423–31.
Sorte K, Sune P, Bhake A, Shivkumar V, Gangane N, Basak A. Quantitative assessment of DNA damage directly in lens epithelial cells from senile cataract patients. Mol Vis. 2011;17:1.
Zhang J, Wu J, Yang L, Zhu R, Yang M, Qin B, et al. DNA damage in lens epithelial cells and peripheral lymphocytes from age-related cataract patients. Ophthalmic Res. 2013;51:124–8.
Liegel RP, Handley MT, Ronchetti A, Brown S, Langemeyer L, Linford A, et al. Loss-of-function mutations in TBC1D20 cause cataracts and male infertility in blind sterile mice and Warburg micro syndrome in humans. Am J Hum Genet. 2013;93:1001–14.
Harper JW, Elledge SJ. The DNA damage response: ten years after. Mol Cell. 2007;28:739–45.
Jackson SP, Bartek J. The DNA-damage response in human biology and disease. Nature. 2009;461:1071–8.
Ryu KY, Maehr R, Gilchrist CA, Long MA, Bouley DM, Mueller B, et al. The mouse polyubiquitin gene UbC is essential for fetal liver development, cell-cycle progression and stress tolerance. EMBO J. 2007;26:2693–706.
Jahngen-Hodge J, Cyr D, Laxman E, Taylor A. Ubiquitin and ubiquitin conjugates in human lens. Exp Eye Res. 1992;55:897–902.
Morishita H, Eguchi S, Kimura H, Sasaki J, Sakamaki Y, Robinson ML, et al. Deletion of autophagy-related 5 (Atg5) and Pik3c3 genes in the lens causes cataract independent of programmed organelle degradation. J Biol Chem. 2013;288:11436–47.
Kimura Y, Tanaka K. Regulatory mechanisms involved in the control of ubiquitin homeostasis. J Biochem. 2010;147:793–8.
Imai F, Yoshizawa A, Fujimori-Tonou N, Kawakami K, Masai I. The ubiquitin proteasome system is required for cell proliferation of the lens epithelium and for differentiation of lens fiber cells in zebrafish. Development. 2010;137:3257–68.
Liou J-Y, Deng W-G, Gilroy DW, Shyue S-K, Wu KK. Colocalization and interaction of cyclooxygenase-2 with caveolin-1 in human fibroblasts. J Biol Chem. 2001;276:34975–82.
Lo W-K, C-j Z, Reddan J. Identification of caveolae and their signature proteins caveolin 1 and 2 in the lens. Exp Eye Res. 2004;79:487–98.
Zhu H, Yue J, Pan Z, Wu H, Cheng Y, Lu H, et al. Involvement of Caveolin-1 in repair of DNA damage through both homologous recombination and non-homologous end joining. PLoS One. 2010;5:e12055.
Rhim JH, Kim JH, Yeo E-J, Kim JC, Park SC. Caveolin-1 as a novel indicator of wound-healing capacity in aged human corneal epithelium. Mol Med. 2010;16:527.
Ashby RS, Zeng G, Leotta AJ, Dennis YT, McFadden SA. Egr-1 mRNA Expression Is a Marker for the Direction of Mammalian Ocular Growth. Invest Ophthalmol Vis Sci. 2014;55:5911–21.
Abate C, Luk D, Gagne E, Roeder RG, Curran T. Fos and jun cooperate in transcriptional regulation via heterologous activation domains. Mol Cell Biol. 1990;10:5532–5.
Mahner S, Baasch C, Schwarz J, Hein S, Wolber L, Janicke F, et al. C-Fos expression is a molecular predictor of progression and survival in epithelial ovarian carcinoma. Br J Cancer. 2008;99:1269–75.
Bossis G, Malnou CE, Farras R, Andermarcher E, Hipskind R, Rodriguez M, et al. Down-regulation of c-Fos/c-Jun AP-1 dimer activity by sumoylation. Mol Cell Biol. 2005;25:6964–79.
Rinaudo JAS, Zelenka PS. Expression of c-fos and c-jun mRNA in the developing chicken lens: relationship to cell proliferation, quiescence, and differentiation. Exp Cell Res. 1992;199:147–53.
The authors declare that they have no competing interests.
JL design and conceive the experiments, approved the manuscript; WjZ and WqZ acquired and analyzed the data; JZ acquired the data; DW drafted the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Zhao, W., Zhao, W., Zhao, J. et al. Screening of potential target genes for cataract by analyzing mRNA expression profile of mouse Hsf4-null lens. BMC Ophthalmol 15, 76 (2015). https://doi.org/10.1186/s12886-015-0066-3
- Differentially expressed gene
- Heat shock transcription factor 4