Identification of Candidate Klf4 Target Genes Reveals the Molecular Basis of the Diverse Regulatory Roles of Klf4 in the Mouse Cornea

Abstract

purpose. Krüppel-like factor 4 (Klf4) plays a crucial role in the development and maintenance of the mouse cornea. In the current study, wild-type (WT) and Klf4-conditional null (Klf4CN) corneal gene expression patterns were examined, to gain understanding of the molecular basis of the Klf4CN corneal phenotype.

methods. Expression of more than 22,000 genes in 10 WT and Klf4CN corneas was compared by microarrays, analyzed using BRB ArrayTools (National Cancer Institute, Bethesda, MD) and validated by Q-RT-PCR. Transient cotransfections were used to test whether Klf4 activates the aquaporin-3, Aldh3a1, and TKT promoters.

results. Scatterplot analysis identified 740 and 529 genes up- and downregulated by more than twofold, respectively, in the Klf4CN corneas. Cell cycle activators were upregulated, whereas the inhibitors were downregulated, consistent with the increased Klf4CN corneal epithelial cell proliferation. Desmosomal components were downregulated, consistent with the Klf4CN corneal epithelial fragility. Downregulation of aquaporin-3, detected by microarray, was confirmed by immunoblot and immunohistochemistry. Aquaporin-3 promoter activity was stimulated 7- to 10-fold by cotransfection with pCI-Klf4. The corneal crystallins Aldh3A1 and TKT were downregulated in the Klf4CN cornea, and their respective promoter activities were upregulated 16- and 9-fold by pCI-Klf4 in cotransfections. The expression of epidermal keratinocyte differentiation markers was affected in the Klf4CN cornea. Although the cornea-specific keratin-12 was downregulated, most other keratins were upregulated, suggesting hyperkeratosis.

conclusions. Functionally diverse candidate Klf4 target genes were identified, revealing the molecular basis of the diverse aspects of the Klf4CN corneal phenotype. These results establish Klf4 as an important node in the genetic network of transcription factors regulating the corneal homeostasis.

Clear vision requires proper development and maintenance of the cornea, a multilayered tissue comprising an outer stratified squamous epithelium, an inner monolayered endothelium and a central stroma with regularly arranged collagen lamellae sparsely populated by keratocytes. The molecular and cellular mechanisms involved in the development and maintenance of the transparence, refractive, and barrier functions of the cornea are exquisitely regulated. ^{1 2 3 4 5 6 7 8} Defective development and maintenance of the cornea result in severe defects in vision. ^{9 10} Different transcription factors influencing corneal morphogenesis and their target genes are known. ^{11 12 13 14 15 16 17 18 19 20 21 22 23 24 25} Despite these advances, our knowledge of the genetic network of transcription factors regulating embryonic morphogenesis, postnatal maturation, and maintenance of cornea remains incomplete. 

Klf4, a member of the Krüppel-like transcription factors (Klf) subfamily of Cys2-His2 zinc finger proteins capable of binding the “GT box” or “CACCC” element, is one of the most highly expressed transcription factors in both 9-day- and 6-week-old mouse corneas. ^{20 26 27 28 29} During development, the expression of Klf4 begins in a stripe of mesenchymal cells extending from the forelimb bud to the developing eye at approximately embryonic day (E)10. ³⁰ In the adult mouse, Klf4 is widely expressed in postmitotic epithelia of diverse tissues, including skin, gastrointestinal tract, and cornea. ^{27 31 32} Klf4-null mice die within 15 hours after birth because of late-stage defects in skin barrier formation. ³³ Klf4-conditional null (Klf4CN) corneas, derived by mating Klf4-LoxP mice ³⁴ with Le-Cre mice, ^{17 35} have multiple ocular defects, including corneal epithelial fragility, stromal edema, smaller, vacuolated lens, and loss of conjunctival goblet cells. ³⁶ To investigate the changes in gene expression underlying the Klf4CN corneal phenotype, we have compared the gene expression patterns of the wild-type (WT) and Klf4CN corneas by microarray hybridization in the present report. Our results show that Klf4 plays a significant role in the maintenance of corneal homeostasis by regulating a wide array of genes encompassing a diverse spectrum of functional subgroups, such as regulators of cell proliferation, cell adhesion molecules, corneal crystallins, components of epithelial barrier function, and regulators of stromal hydration. 

Materials and Methods

Conditional Deletion of Klf4

Klf4CN mice were generated on a mixed background by mating Klf4 ^loxP/loxP , Le-Cre/⁻ mice with Klf4 ^loxP/loxP mice as described before. ^{34 35 36} This mating scheme generated equal numbers of Klf4 ^loxP/loxP , Le-Cre/⁻ (Klf4CN), and Klf4 ^loxP/loxP (control) offspring. Genomic DNA isolated from tail clippings was assayed for the presence of the Klf4-LoxP and Le-Cre transgenes by PCR using specific primers. The mice in the present study were maintained in accordance with the guidelines set forth by the Animal Care and Use Committee of the National Eye Institute, National Institutes of Health, and the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research. 

Isolation of Total RNA, Quality Control, Labeling, and Microarray Analysis

In the present analysis, we used the whole cornea, comprising epithelial cells, stromal keratocytes, and endothelial cells, as well as a small number of infiltrating leukocytes. Similar microarray analyses of whole corneas have been useful in identifying the corneal responses to Aspergillus fumigatus ³⁷ or Pseudomonas aeruginosa ³⁸ infections and diabetic conditions ³⁹ and in characterizing the healing process after laser ablation ⁴⁰ or keratectomy. ⁴¹ Five age-matched, 8-week-old WT and Klf4CN mice each were used for comparison of corneal gene expression by microarray analysis. Two dissected corneas from each mouse were combined for isolation of total RNA (RNeasy Mini kit; Qiagen, Valencia, CA). The RNA was sent to Expression Analysis (Durham, NC) which confirmed the integrity of the isolated total RNA by using a bioanalyzer (Agilent Technologies, Palo Alto, CA) (Supplementary Fig. S1) and subsequently performed the microarray hybridizations. Labeled samples were hybridized with mouse gene arrays (Mouse 430 2.0; Affymetrix, Santa Clara, CA) containing 45,101 panels, each targeting a specific nucleic acid sequence. In these arrays, approximately 22,000 transcripts identified by Entrez Gene numbers (Entrez Database; National Center for Biotechnology Information, Bethesda, MD, available at www.ncbi.nlm.nih.gov/sites/entrez) are redundantly targeted by 41,400 panels: The remaining panels target relatively less-characterized sequences. The raw data obtained from Expression Analysis was analyzed using BRB-ArrayTools software (developed by Richard Simon and Amy Peng Lamthe, Biometric Research Branch, National Cancer Institute, Bethesda, MD, available at linus.nci.nih.gov/brb-arraytools.html/). Microarray data were normalized by using median over entire array and filtered according to the following criteria: The genes were excluded from analysis if (1) the probability was greater than 0.005, (2) the percentage of data missing or filtered out were greater than 50%, (3) greater than 20% of expression data showed more than a 1.5-fold change in either direction from the median value, or (4) the detection call was “absent” in more than 50% in both WT and Klf4CN. The microarray results provided in this report are log transformed. 

Validation of Microarray Results Using RT-PCR and Real-Time Quantitative RT-PCR

Total RNA isolated from the WT or Klf4CN corneas was quantified, the concentration adjusted with RNase-free water to 100 ng/μL, and one-step RT-PCR performed with 100 ng total RNA and RT-PCR beads (Ready-To-Go; GE Healthcare, Piscataway, NJ). To distinguish the products amplified from the contaminating genomic DNA if any, from those originating from the mRNA, the forward and reverse primers used in RT-PCR were chosen from adjacent exons. The sequence of primers used for RT-PCR is provided in Supplementary Table S1. The RT-PCR products were separated on a 1.5% agarose gel with TBE buffer. 

Applied Biosystems, Inc. ([ABI], Foster City, CA) was the source of the reagents, equipment, and software for gene expression quantitative real-time RT-PCR assays (Q-RT-PCR; TaqMan; ABI). cDNA was generated (High Capacity cDNA Archive Kit; ABI), and total RNA was isolated from pooled corneas of 10 WT or Klf4CN mice. Q-RT-PCR assays with prestandardized gene-specific probes for different transcripts were performed in a thermocycler (model 7900HT; ABI), with 18S rRNA as the endogenous control; the results were analyzed with commercial software (SDS software ver. 2.1; ABI). 

Immunoblots and Immunohistochemistry

Total protein was extracted by homogenizing dissected corneas in 8.0 M urea, 0.08% Triton X-100, 0.2% sodium dodecyl sulfate, 3% β-mercaptoethanol, and proteinase inhibitors and was quantified by the bicinchoninic acid method (Pierce, Rockford, IL). Equal amounts of protein were electrophoresed in sodium dodecyl sulfate-polyacrylamide gels, transferred to polyvinylidene difluoride membranes and subjected to immunoblot analysis. Rabbit anti-aquaporin-3 (Calbiochem, La Jolla, CA), and anti-actin antibody (Sigma-Aldrich, St. Louis, MO) were used as primary antibodies at a 1:1000 dilution in PBST. Horseradish peroxidase-coupled anti-rabbit immunoglobulin G (GE Healthcare, Piscataway, NJ) was used as a secondary antibody at a 1:5000 dilution. Immunoreactive bands were visualized by chemiluminescence (Super Signal West Pico solutions; Pierce). 

For immunohistochemistry, 10-μm-thick cryosections from OCT-embedded eyeballs were fixed in freshly prepared, buffered 4% paraformaldehyde for 30 minutes, blocked with 10% sheep serum in PBST for 1 hour at room temperature in a humidified chamber, washed twice with PBST for 5 minutes each, incubated with a 1:100 dilution of the primary antibody for 1 hour at room temperature, washed thrice with PBST for 10 minutes each, incubated with secondary antibody (AlexaFluor 555-coupled goat anti-rabbit IgG antibody; Invitrogen-Molecular Probes, Carlsbad, CA) at a 1:300 dilution for 1 hour at room temperature, washed thrice with PBST for 10 minutes each, mounted with antifade reagent with DAPI (Prolong Gold; Molecular Probes), and observed with a fluorescence microscope (Axioplan 2; Carl Zeiss Meditec, Dublin, CA). 

Reporter Vectors, Cell Culture, and Promoter Activities

Mouse genomic DNA was used to amplify different promoter fragments used in the cotransfection assays reported herein. Aquaporin3 (Aqp3) −502/+42- and −262/+42-bp promoter fragments were amplified by using the downstream Aqp3 +42/+22C HindIII (+42 ATGCAAGCTTGTCCGGCGGCGTACGAGTGC +22C) and upstream Aqp3 −502/−482 XhoI (−502 ATGCCTCGAGCACGAAGCGCTGGTGAATTC −482) or Aqp3 −262/−245 XhoI (−262 ATGCCTCGAGGGAGACCGCTTGCTCTTC −245) primers. Transketolase (TKT) −518/+104-bp promoter fragment was amplified by using upstream TKT −518/−491 KpnI (−518 GGCCGGTACCGGCAAACCCAGTAATCTC −491) and downstream TKT +104/+87-bp HindIII (+104 GGCCAAGCTTCCTTCCATGGCGTGGTAGG +87) primers. Aldh3a1 −1050/+3486-bp promoter fragment was isolated as described previously. ⁴² These promoter fragments were cloned upstream of the luciferase reporter gene in a vector (pGL3Basic; Promega, Madison, WI) to generate the reporter vectors. Full-length Klf4 was transiently expressed by using the CMV promoter in pCI-Klf4. Monkey kidney Cos7 cells were grown in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum, penicillin, and streptomycin at 37°C in a humidified chamber containing 5% CO₂ and 95% air. Simian virus SV40-transformed human corneal epithelial (HCE) cells ⁴³ were grown at 37°C in Ham’s F-12 supplemented with 10% fetal bovine serum, 0.5% (vol/vol) dimethyl sulfoxide, cholera toxin (0.1 μg/mL), epidermal growth factor (10 ng/mL), insulin (5 μg/mL), gentamicin (40 μg/mL), and glutamine (20 mM) in a humidified chamber containing 5% CO₂ and 95% air. Mid-log phase cells in six-well plates were transfected with 0.5 μg of pAldh3A1-Luc, pTKT-Luc, or pAqp3-Luc, along with 10 ng pRL-SV40, for normalization of transfection efficiency (Promega), and 0.5 μg of pCI or pCI-Klf4, using 3 μL of transfection reagent (Fugene 6; Roche Molecular Biochemicals, Indianapolis, IN). After 2 days, the cells were washed with cold PBS and lysed with 500 μL passive lysis buffer (Promega). The lysate was clarified by centrifugation and 50 μg protein in the lysate was analyzed with a dual-luciferase assay kit (Promega) and a microplate luminometer (Victor; Perkin-Elmer, Wellesley, MA). The measurement was integrated over 10 seconds with a delay of 2 seconds. Results from at least three independent experiments, normalized for transfection efficiency using the SV40 promoter-driven Renilla luciferase activity, were used to obtain mean promoter activities and standard deviation. Activation (x-fold multiples of change) was determined by dividing mean promoter activity by the promoter activity without added pCI/pCI-Klf4. 

Results

Microarray Analysis and Validation of Results

To obtain mechanistic insight into the diverse ocular surface phenotype observed in the Klf4CN mice,³⁶we compared the gene expression patterns between the WT and Klf4CN corneas by microarray hybridization. Scatterplot analysis of the 6333 genes that passed the filtering criteria described in Materials and Methods showed that 529 genes were downregulated and 740 genes were upregulated by more than twofold in the Klf4CN compared with the WT corneas (Fig. 1 , Supplementary Tables S2A, S2B). Microarray results were validated by quantitative real-time RT-PCR comparison of expression levels of 19 genes (Table 1) . There was a general conformity between the microarray and Q-RT-PCR results, indicating that the microarray results reflect the extent of changes in gene expression in the Klf4CN corneas (Table 1) . The candidate Klf4 target genes, whose expression was significantly affected in the Klf4CN corneas, fall into different functional subgroups regulating diverse functions such as cell proliferation, apoptosis, development, immune response, and barrier function. Below, we have analyzed the microarray results further, to correlate different aspects of the Klf4CN corneal phenotype with specific changes in gene expression, thus revealing the diverse contributions of Klf4 to corneal physiology. 

Role of Klf4 in Regulating the Progression of Corneal Epithelial Cell Cycle

In view of the increased rate of Klf4CN corneal epithelial cell proliferation, ³⁶ we examined the expression pattern of a subset of genes known to participate in or regulate cell cycle. Several cell proliferation activators were upregulated while cell cycle suppressors were downregulated in the Klf4CN cornea, compared to the WT (Table 2) . Important activators of the cell cycle that were upregulated include cyclin D2, cyclin-dependent kinase 6, receptor type protein tyrosine phosphatase-b and -c, and FMS-like tyrosine kinase-1. Significantly downregulated inhibitors of cell cycle include retinoblastoma-1 (Rb) and cyclin-dependent kinase inhibitor 1a (p21). The upregulation of activators of the cell cycle and the downregulation of inhibitors of the cell cycle are consistent with the increased cell proliferation in the Klf4CN corneal epithelium. ³⁶  

Role of Klf4 in Regulation of Epithelial Cell–Cell Adhesion

Despite the increased rate of cell proliferation, the Klf4CN cornea possessed fewer epithelial cell layers than the WT, suggesting reduced cell–cell adhesion at the Klf4CN ocular surface. ³⁶ To explore this possibility further, we examined the expression of different components of desmosomes, the intercellular adhesion complexes prevalent in the stratified epithelia. ⁴⁴ The expression of most of the desmosomal components was significantly downregulated in the Klf4CN cornea, strongly supporting the notion that compromised cell–cell adhesion is responsible for the reduced numbers of Klf4CN corneal epithelial cell layers (Fig. 2A) . Comparison of expression of these desmosomal components in the WT and Klf4CN corneas by RT-PCR validated the microarray results (Fig. 2B) . 

Role of Klf4 in Regulation of Corneal Hydration

Aquaporin-1, -3, and -5 are responsible for keeping the hydrophilic corneal stroma from abnormal swelling. ^{45 46} We have reported that the downregulation of Aqp5 may be responsible for the Klf4CN corneal stromal edema. ³⁶ In addition to confirming the downregulation of Aqp5, the present microarray analysis revealed that the expression of Aqp3 also is reduced in the Klf4CN cornea to approximately half that in the WT, whereas the expression of Aqp1 remains relatively unaffected. Downregulation of Aqp3 transcripts and the 38-kDa Aqp3 protein was validated by Q-RT-PCR (Table 1 ; Fig. 3A ) and immunoblot analysis (Fig. 3B) , respectively. Equal loading ofproteins for immunoblot analysis was confirmed by stripping the membrane of antibodies and reprobing with an anti-actin antibody, which did not show any difference between WT and Klf4CN corneas (Fig. 3B) . Immunohistochemistry confirmed the reduction in the amount of Aqp3 in the Klf4CN corneal epithelium (Fig. 3C)consistent with the reduced amount of Aqp3 in the Klf4CN corneal extracts compared with the WT extracts (Fig. 3B) . In view of the reduced stratification of the Klf4CN corneal epithelium (Fig. 3C) , it was uncertain to what extent the reduced Aqp3 protein band in the immunoblot (Fig. 3B)was due to reduced expression of Aqp3 in the epithelium or to the relatively large stromal cell representation in the Klf4CN corneal extracts inasmuch as stromal cells do not express Aqp3 (present study and Refs. ^{47 48} ), yet contain actin which was used for normalization. To test more directly whether Klf4 stimulates Aqp3 promoter activity, we used cotransfection assays in cultured cells. Both −502/+42- and −262/+42-bp Aqp3 promoter activities were upregulated 8- to 9-fold by cotransfection with pCI-Klf4, indicating that the activating sites reside within the −262/+42-bp region of the Aqp3 proximal promoter (Fig. 3D) . Taken together, these results show that the Aqp3 promoter is regulated by Klf4 and suggest that the cumulative effect of the reduced expression of Aqp3 and Aqp5 is responsible for the observed Klf4CN stromal edema. ³⁶  

Role of Klf4 in the Regulation of Corneal Crystallins

Corneas, like lenses, accumulate unusually high proportions of a few water-soluble proteins termed corneal crystallins, in a taxon-specific manner. ⁴⁹ In the mouse corneal epithelium, aldehyde dehydrogenase-3a1 (Aldh3a1) and transketolase (TKT) comprise 20% to 50% and 10% to 15%, respectively, of the water-soluble protein. ^{50 51} Microarray comparisons indicated that the expression of Aldh3a1 and TKT genes was downregulated in the Klf4CN cornea to 90% and 30% of the WT, respectively (Fig. 4A) . This downregulation was estimated more quantitatively by Q-RT-PCR, which showed that the expression of Aldh3a1 and TKT genes in the Klf4CN cornea was downregulated to 41% and 56% of the WT, respectively (Fig. 4B) . To test directly the influence of Klf4 on Aldh3A1 and TKT promoter activities, we used cotransfection assays in cultured Cos7 cells, in which the −1050/+3486-bp Aldh3A1 promoter and −518/+104-bp TKT promoter activities were upregulated by 16- and 9-fold, respectively, in cells cotransfected with pCI-Klf4, compared to those cotransfected with pCI (Fig. 4C) . Taken together, these results suggest that Klf4 contributes to the corneal transparence and refractive properties, by activating the expression of corneal crystallins Aldh3A1 and TKT. 

Role of Klf4 in Regulating the Barrier Function of Corneal Epithelium

Several genes encoding proteins related directly or indirectly to epithelial barrier function were downregulated in the Klf4CN cornea. Uroplakin 3B, an integral membrane protein that contributes to the apical membrane permeability barrier function, ⁵² as well as gastrokine-1 (also known as antrum mucosal protein-18), a mitogenic protein abundantly expressed in the superficial gastric epithelium and downregulated in gastric carcinoma, ⁵³ were significantly downregulated in the Klf4CN cornea (Supplementary Table S2B). The expression of a large number of immune response-related genes was affected in the Klf4CN cornea (Supplementary Tables S2A, S2B). The upregulated immune response–related genes belonged to diverse groups, such as complement components, chemokines, and chemokine receptors. Infiltration of the immune cells into the Klf4CN cornea in response to the inflammatory signals generated by the fragile corneal epithelium and edematous stroma may account for the observed upregulation of the immune response–related genes in the Klf4CN cornea. 

Keratins, the fibrous intermediate-filament protein polymers, play critical roles in maintenance of the epithelial cell structure, protecting them from mechanical trauma. We have shown that Krt12 promoter is bound and upregulated by Klf4 and that the downregulation of Krt12 may be responsible for the Klf4CN corneal epithelial fragility. ³⁶ Microarray analysis showed that while Krt12 was indeed downregulated, most of the other keratins were upregulated, indicating hyperkeratosis in the Klf4CN cornea (Fig. 5A) . RT-PCR assays validated these microarray results (Fig. 5B) . 

Role of Klf4 in Regulation of Genes Involved in Mouse Epidermal Keratinocyte Differentiation

In view of the fact that Klf4 plays a crucial role in the mouse epidermal development, ³³ we examined the expression pattern of different genes involved in the mouse epidermal keratinocyte differentiation. Oncogenes Ets1 and Jun, as well as tumor necrosis factor receptor Tnfrsf1b, were upregulated, whereas the expression of epidermal growth factor (egf), protein kinase-Cβ1, and mitogen-activated protein kinase kinase-6 (Mapkk6) were downregulated by more than twofold in the Klf4CN cornea, compared with the WT (Supplementary Table S2A; Table 3 ). Expression of epidermal keratinocyte differentiation markers, such as involucrin and transglutaminase was affected in the Klf4CN cornea. Expression of several members of the Sprr family of proteins, constituents of the cornified envelope in the skin, was upregulated in the Klf4CN cornea (Supplementary Table S2A; Table 3 ). 

Role of Klf4 in the Regulation of Gene Regulatory Networks in the Cornea

To understand the impact of conditional deletion of Klf4 on gene regulatory networks in the cornea, we examined the expression levels of different transcription factors in the Klf4CN cornea. Downregulation of the paired box–homeobox transcription factor Pax6, central to eye development, ⁵⁴ by more than twofold, as detected by microarray analysis (Table 4) , was validated by Q-RT-PCR (Table 1) . In addition, the special AT-rich sequence-binding protein-1 (SATB1), which organizes nuclear architecture and recruits chromatin remodeling factors to specific promoters, ^{55 56 57} was downregulated by approximately threefold (Table 4) . Another significantly downregulated transcription factor, inhibitor of DNA binding-2 (Id2), is known to contribute to epithelial cell differentiation and suppress tumor formation. ⁵⁸ In contrast, several stress response–related transcription factors, such as Ahr, ⁵⁹ E2F5 and EGR1, ⁶⁰ NFκB2, ⁶¹ and ATF3 ⁶² were significantly upregulated in the Klf4CN cornea, compared with the WT cornea (Table 4) . 

Discussion

Previously, we demonstrated that the structural integrity of the mouse corneal epithelium, stroma, and endothelium was affected when the Klf4 gene was deleted, indicating that Klf4 plays a crucial role in the development and maintenance of the mouse cornea. ³⁶ We also demonstrated that the expression levels of Krt12 and Aqp5 are reduced in the Klf4CN cornea, consistent with the Klf4CN corneal epithelial fragility and stromal edema, respectively. Here, we have used microarray analysis to obtain a comprehensive view of the changes in corneal gene expression on deletion of Klf4. Our findings reveal the molecular basis of the wide-ranging influence of Klf4 on corneal homeostasis and identify candidate target genes of Klf4 in the adult mouse cornea. 

Previous attempts at Klf4 target gene profiling used a human cell line containing inducible Klf4, ^{63 64} or transgenic mice overexpressing Klf4 in the epidermis under the control of keratin-5 promoter (Krt5-Klf4) and the Klf4 ^−/− mice. ⁶⁵ In addition to being largely consistent with the earlier reports, our results identify novel Klf4 target genes with specialized functions in the mouse cornea, such as aquaporins-3 and -5, the corneal crystallins Aldh3a1 and transketolase, and desmosomal components. 

Consistent with the established mechanism of suppression of cell proliferation by Klf4, ^{66 67} cyclin-D2 was upregulated and cdkn1a (p21) was downregulated in the Klf4CN cornea. Furthermore, our results are consistent with those in a previous report that showed several other inhibitors of cell division to be upregulated, and activators of cell cycle to be downregulated by Klf4, ⁶⁴ suggesting additional mechanisms by which Klf4 can suppress the progression of cell cycle (Tables 2A 2B) . 

Many members of the large family of keratins, the major components of the intermediate filaments protecting the epithelial cells from mechanical and nonmechanical stresses, ⁶⁸ are upregulated in the Klf4CN corneas. In contrast to our results, most keratins were upregulated by Klf4 in a Klf4-inducible cell culture system. ⁶³ This difference could be due to the different experimental systems used and/or to an indirect hyperkeratinizing response in the Klf4CN cornea, as a consequence of the loss of corneal barrier function. 

A striking feature of our results is the large number of genes whose expression is affected by the absence of Klf4 in the mouse cornea. Although some of these genes are likely to be direct targets of Klf4, the remaining genes may be indirect targets through other transcription factors regulated by Klf4. The expression of Pax6 is reduced to about half in the Klf4CN cornea compared to the WT cornea. Thus, Pax6 target genes are likely to be included in the list of genes affected by the absence of Klf4. In this regard, it is noteworthy that the Pax6^+/− corneal epithelium appears remarkably similar to the Klf4CN corneal epithelium, with fewer cell layers, roughened ocular surface, and epithelial vacuolation. ^{15 69} Other transcription factors with reduced (Id2, SATB1, and Elk4) or increased (NFκB, Sox4, Atf3, C/EBPδ. GATA2, Ahr and Egr1) expression in the Klf4CN cornea may indirectly contribute to the list of genes with changed expression in the Klf4CN cornea. 

In addition to these changes, the expression of many genes with no established functions in the cornea and/or no apparent relevance to the Klf4CN corneal phenotype was significantly affected in the Klf4CN cornea. Three members of the Ly6/Plaur domain containing the proteins Slurp1, one of the most abundant transcripts in the mouse cornea, ²⁷ a ligand for nicotinic acetylcholine receptors and a late marker of epidermal differentiation associated with the inflammatory palmoplantar keratoderma disease Mal de Meleda ^{70 71 72} ; Lynx1 (also a ligand for nicotinic acetylcholine receptors ⁷³ ); and Lypd2 were significantly downregulated in the Klf4CN cornea (Supplementary Table S2B). Similarly, the expression of 15 and 9 different members of the solute carrier family of proteins was up- and downregulated respectively, in the Klf4CN compared with the WT cornea (Supplementary Tables S2A, S2B). Whether these changes contribute to any aspect of the Klf4CN corneal phenotype remains to be established. 

The results presented in this report show that Klf4 coordinately regulates functionally related subsets of genes, such as those contributing to the control of corneal epithelial cell cycle progression, intercellular adhesion, corneal crystallins, the Ly6/Plaur domain containing proteins Slurp1, Lypd2, and Lynx1, ^{70 71 72 73} and the small proline-rich proteins (SPRR), the primary constituents of the cornified cell envelope and integral components of the surface barrier. ^{74 75} We have also shown that Klf4 stimulates the promoter activities of aquaporin-3 and -5, ³⁶ and corneal crystallins Aldh3A1 and TKT in cultured cells. It remains to be established whether Klf4 plays a direct role in the coordinate regulation of the remaining groups of genes whose expression is affected in the Klf4CN cornea. A fraction of the observed changes in gene expression could be indirect, such as a response to the inflammatory conditions caused by the fragile Klf4CN corneal epithelium. The loss of epithelial barrier function may be responsible for the overexpression of several stress-related genes in the Klf4CN cornea, such as the antioxidant enzyme ceruloplasmin, which is upregulated in different neurodegenerative disorders including glaucoma ^{76 77} ; arachidonate lipoxygenase-12 and -15, which promote epithelial wound healing and host defense ⁷⁸ ; and carbonic anhydrase-2, -12, and -13, regulators of corneal ion transport, that are overexpressed in human glaucoma ^{79 80} (Supplementary Tables S2A, S2B). 

In summary, the changes in gene expression patterns detected by the present microarray analysis are consistent with the phenotypic changes in the Klf4CN cornea. Our results show that Klf4 contributes to corneal homeostasis by coordinately regulating the expression of subsets of genes involved in specific functions such as progression of the cell cycle, cell–cell adhesion, epithelial barrier formation, expression of corneal crystallins, and maintenance of corneal hydration. Taken together with the findings in our earlier report, ³⁶ the present studies establish Klf4 as an important component in the genetic network of transcription factors necessary for proper development and maintenance of the ocular surface. 

 Supported by the intramural research program of the National Eye Institute, NEI Career Development Award 1 K22 EY016875-01 (SKS), startup funds from the Department of Ophthalmology, Core Grant for Vision Research 5P30 EY08098-19, Research to Prevent Blindness, and the Eye and Ear Foundation, Pittsburgh.

 Submitted for publication January 30, 2008; revised March 27, 2008; accepted June 16, 2008.

 Disclosure: S.K. Swamynathan, None; J. Davis, None; J. Piatigorsky, None

 The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked “advertisement” in accordance with 18 U.S.C. §1734 solely to indicate this fact.

 

Figure 1.

 

Scatterplot analysis of 6333 genes passing the filter described in Materials and Methods. Expression of 529 genes was downregulated and 740 genes upregulated by more than twofold in Klf4CN compared with the WT corneas. The expression of approximately 80% of the genes (5064 genes) passing the filter remained relatively stable with less than twofold change.

Table 1.

 

Validation of Microarray Analysis Results by Real-Time Quantitative RT-PCR of Selected Genes

Gene Symbol	Difference
		Microarray	Q-RT-PCR
Alox12e	19.64	61.34
Alox15	86.21	137.26
Aqp3	0.44	0.45
Aqp5	0.09	0.24
Krt1-12	0.91	0.33
Krt1-17	57.79	80.59
Krt2-4	5.98	10.7
Lumican	1.08	1.36
CcnD2	5.86	5.71
Lamb1-1	0.29	0.34
Muc1	3.94	7.04
Sprr2A	77.43	20.01
ALDH3A1	0.92	0.41
ELF3	0.76	0.65
IRF1	0.6	0.56
Pax6	0.38	0.55
SLURP1	0.02	0.1
GSTO1	1.33	0.84
SPARC	1.41	1.7

 Data are expressed as x-fold difference.

Table 2.

 

Expression Levels of Genes Involved in Cell Cycle Regulation in Klf4CN Compared with WT Corneas

Description	Gene Symbol	WT	Klf4CN	Difference (x-Fold)
Upregulated genes
Cyclin D2 (Validated by Real time Q-RT-PCR)	Ccnd2	6.86	9.41	5.86
Cyclin-dependent kinase 6	Cdk6	9.34	10.71	2.58
Protein kinase, cGMP-dependent, type II	Prkg2	5.80	8.29	5.64
Protein kinase, cAMP-dependent regulatory, type I, alpha	Prkar1a	5.68	6.86	2.27
FMS-like tyrosine kinase 1	Flt1	5.64	6.95	2.48
Protein kinase C, mu	Prkcm	5.71	7.20	2.80
Protein tyrosine phosphatase, receptor type, B	Ptprb	5.84	7.25	2.66
Protein tyrosine phosphatase, receptor type, C	Ptprc	6.04	8.41	5.19
Early growth response 1	Egr1	9.70	11.38	3.20
Jun oncogene	Jun	9.69	10.86	2.26
Expressed in non-metastatic cells 1, protein	Nme1	11.70	12.77	2.11
Platelet derived growth factor receptor, beta polypeptide	Pdgfrb	8.08	9.25	2.26
Eph receptor A4	Epha4	9.38	10.72	2.52
Downregulated genes
Eph receptor B6	Ephb6	8.07	6.95	0.46
Ephrin A1	Efna1	10.04	8.9	0.45
Ephrin B1	Efnb1	9.12	7.96	0.45
Cyclin-dependent kinase inhibitor 1a (p21)	Cdkn1a	11.76	10.92	0.56
Adenylate kinase 1	Ak1	9.41	8.11	0.41
NIMA (never in mitosis gene a)-related expressed kinase 3	Nek3	9.11	7.4	0.31
Macrophage stimulating 1 receptor (c-met–related tyrosine kinase)	Mst1r	8.9	7.89	0.5
Protein kinase C, beta 1	Prkcb1	7.56	6.53	0.49
Protein kinase C, iota	Prkci	10.67	9.28	0.38
Noncatalytic region of tyrosine kinase adaptor protein 1	Nck1	11.46	10.44	0.49
Rous sarcoma oncogene	Src	9.01	7.99	0.49
Fibroblast growth factor receptor 2	Fgfr2	10.28	9.08	0.43
RAN binding protein 2	Ranbp2	5.5	4.19	0.4
Retinoblastoma 1	Rb1	6.48	4.24	0.21
Growth arrest and DNA-damage-inducible 45 alpha	Gadd45a	10.95	9.82	0.46
Growth arrest specific 2	Gas2	8.55	7.32	0.43

 The values shown for the WT and Klf4CN corneas are log transformed.

Figure 2.

 

(A) Expression levels of different desmosomal components in the WT and Klf4CN corneas as detected by microarray. The values shown are log transformed. (B) Validation of downregulation of genes encoding desmosomal components in the Klf4CN (CN) compared with the WT corneas, by RT-PCR. M, molecular weight markers.

Figure 3.

 

Aquaporin-3 expression was reduced in the Klf4CN relative to the WT cornea. (A) Q-RT-PCR comparison of expression levels of Aqp3 in the WT and Klf4CN corneas. (B) Immunoblot analysis of Aqp3 levels in total proteins extracted from WT and Klf4CN corneas. The blot probed with anti-Aqp3 antibody (left) was stripped of the antibody and reprobed with anti-actin antibody (right) to ensure equal loading of proteins. (C) Immunohistochemistry with anti-Aqp3 antibody. Images show (top row) nuclei stained with DAPI (bottom row) fluorescence coming from the secondary antibody bound to the primary anti-Aqp3 antibody. The expression of Aqp3 localized to the epithelial and endothelial cell membranes in the WT (bottom, left) was reduced in the Klf4CN cornea (bottom, center). The sections processed in a similar manner without the primary antibody served as controls (bottom, right). (D) Aquaporin-3 promoter was stimulated by KLF4 in cultured human corneal epithelial cells cotransfected with Aqp3-Luc and pCI-Klf4.

Figure 4.

 

Expression of corneal crystallins Aldh3a1 and TKT was regulated by Klf4. (A) Comparison of expression of Aldh3a1 and TKT in the WT and Klf4CN cornea, as detected by microarray. (B) Q-RT-PCR comparison of expression of Aldh3a1 and TKT in WT and Klf4CN corneas. (C) Aldh3a1 and TKT promoter activities are stimulated by Klf4 in Cos-7 cells on cotransfection with pCI-Klf4.

Figure 5.

 

(A) Expression levels of different keratins in the WT and Klf4CN corneas as detected by microarray. The values shown are log transformed. (B) Validation of upregulation of genes encoding different keratin genes in the Klf4CN (CN) compared with the WT corneas, by RT-PCR. M, molecular weight markers.

Table 3.

 

Expression Levels of Mouse Epidermal Keratinocyte Differentiation Markers in WT and Klf4CN Corneas

Description	Gene Symbol	WT	Klf4CN	Difference (x-Fold)
Tumor necrosis factor receptor superfamily, member 1b	Tnfrsf1b	5.91	7.30	2.62
E26 avian leukemia oncogene 1, 5′ domain	Ets1	6.84	8.32	2.77
Jun oncogene	Jun	9.69	10.86	2.26
Small proline-rich protein 2a	Sprr2a	5.96	12.24	77.43
Small proline-rich protein 2f	Sprr2f	5.64	10.41	27.25
Small proline-rich protein 2d	Sprr2d	2.26	5.49	25.06
Cystatin B	Cstb	9.91	10.30	1.47
Transglutaminase 1	Tgm1	5.56	6.02	1.58
Transglutaminase 2	Tgm2	9.08	9.57	1.63
Epidermal growth factor	Egf	8.31	7.02	0.41
Protein kinase C, beta 1	Prkcb1	2.02	1.88	0.49
Mitogen activated protein kinase kinase 6	Map2k6	8.39	7.37	0.49
Involucrin	Ivl	5.78	5.18	0.55
Transglutaminase 4	Tgm4	5.43	5.15	0.75
Cystatin A	Csta	9.59	9.10	0.61
Cystatin E/M	Cst6	9.98	8.03	0.26

 The values shown for the WT and Klf4CN corneas are log transformed.

Table 4.

 

Comparison of Expression Levels of Selected Transcription Factors in WT and Klf4CN Corneas

Description	Gene Symbol	WT	Klf4CN	Difference (x-Fold)
Upregulated genes
SRY-box containing gene 4	Sox4	4.49	7.52	8.19
Nuclear Factor Kappa b2	Nfkb2	5.31	7.49	4.52
Activating transcription factor 3	Atf3	5.96	7.73	3.41
E2F transcription factor 5	E2f5	6.08	7.92	3.57
Aryl-hydrocarbon receptor	Ahr	9.87	11.36	2.82
Ribosomal protein S9	Rps9	10.01	12.97	7.81
Early growth response 1	Egr1	9.70	11.38	3.20
CCAAT/enhancer binding protein (C/EBP), delta	Cebpd	8.86	11.41	5.86
Eukaryotic translation initiation factor 4E member 3	Eif4e3	6.13	8.48	5.10
GATA binding protein 2	Gata2	6.11	7.41	2.47
Downregulated genes
Paired box gene 6	Pax6	11.02	9.62	0.38
Special AT-rich sequence binding protein 1	Satb1	9.17	7.14	0.24
Inhibitor of DNA binding 2	Id2	12.35	10.56	0.29
Hepatic leukemia factor	Hlf	13.11	11.88	0.43
ELK4, member of ETS oncogene family	Elk4	5.2	3.5	0.31
Eukaryotic translation initiation factor 2 alpha kinase 4	Eif2ak4	9.30	8.12	0.44

 The values shown for the WT and Klf4CN corneas are log transformed.

Supplementary Materials

Supplementary Figure S1 - (PDF) 

Supplementary Table S1 - (.xls) 

Supplementary Table S2a - (.xls) 

Supplementary Table S2b - (.xls) 

The authors thank Stephen Harvey, University of Pittsburgh, for insightful comments on the manuscript.