Transcriptional profiling of epidermal differentiation - CiteSeerX

8 downloads 82374 Views 4MB Size Report
Jul 5, 2006 - were obtained with Microarray Suite version 5.0 (Affymetrix) and scaled by .... digital images were collected with the Adobe Photoshop TWAIN_32 program. ..... These genes are, of course, not expressed in the β4+ cells.
Physiol Genomics 27: 65–78, 2006. First published July 5, 2006; doi:10.1152/physiolgenomics.00031.2006.

Transcriptional profiling of epidermal differentiation Nada Radoja,1 Alix Gazel,1 Tomohiro Banno,1,2 Shoichiro Yano,1,3 and Miroslav Blumenberg1,4,5 Departments of 1Dermatology and 4Biochemistry, and the 5NYU Cancer Institute, New York University School of Medicine, New York, New York; 2Dermatology Department at the Institute of Clinical Medicine, Tsukuba University, Tsukuba, Ibaraki; and 3Department of Dermatology, Faculty of Medicine, University of Tokyo, Tokyo, Japan Submitted 20 February 2006; accepted in final form 16 June 2006

Radoja, Nada, Alix Gazel, Tomohiro Banno, Shoichiro Yano, and Miroslav Blumenberg. Transcriptional profiling of epidermal differentiation. Physiol Genomics 27: 65–78, 2006. First published July 5, 2006; doi:10.1152/physiolgenomics.00031.2006.—In epidermal differentiation basal keratinocytes detach from the basement membrane, stop proliferating, and express a new set of structural proteins and enzymes, which results in an impermeable protein/lipid barrier that protects us. To define the transcriptional changes essential for this process, we purified large quantities of basal and suprabasal cells from human epidermis, using the expression of ␤4 integrin as the discriminating factor. The expected expression differences in cytoskeletal, cell cycle, and adhesion genes confirmed the effective separation of the cell populations. Using DNA microarray chips, we comprehensively identify the differences in genes expressed in basal and differentiating layers of the epidermis, including the ECM components produced by the basal cells, the proteases in both the basal and suprabasal cells, and the lipid and steroid metabolism enzymes in suprabasal cells responsible for the permeability barrier. We identified the signaling pathways specific for the two populations and found two previously unknown paracrine and one juxtacrine signaling pathway operating between the basal and suprabasal cells. Furthermore, using specific expression signatures, we identified a new set of late differentiation markers and mapped their chromosomal loci, as well as a new set of melanocyte-specific markers. The data represent a quantum jump in understanding the mechanisms of epidermal differentiation. basal cells; epidermal differentiation complex; integrins; keratinization; melanocytes; microarrays DIFFERENTIATION IS THE PROCESS by which progenitor cells acquire functional capabilities. This process, perhaps most studied in lymphoid cells, is ideally suited for global transcriptional profiling because it comprises large transcriptional changes, as cells convert from rapidly proliferating to functionally specialized phenotypes. Indeed, microarrays have been used to study the differentiation processes in human cell lines induced to differentiate along a certain pathway by in vitro treatments with various agents (13, 45). However, so far, no such studies of normally differentiating human cells in vivo have been reported, perhaps because of difficulties in obtaining sufficient quantities of human cell populations at specific stages of differentiation. Epidermis is an ideal target tissue for studies of differentiation because both the progenitor and the differentiating cells can be easily identified and purified in sufficient quantities. Therefore, we decided to develop comprehensive transcriptional profiles of human epidermal basal cells and

Article published online before print. See web site for date of publication (http://physiolgenomics.physiology.org). Address for reprint requests and other correspondence: M. Blumenberg, NYU Cancer Institute, New York Univ. School of Medicine, 550 1st Ave., New York, NY 10016 (e-mail: [email protected]).

their differentiated descendants, using large oligonucleotide microarrays. Four layers can be distinguished morphologically in healthy epidermis, stratum basale, stratum spinosum, stratum granulosum, and stratum corneum (9, 23). The function of stratum basale is to proliferate and to keep the epidermis firmly anchored to the basal lamina. The function of stratum spinosum is to start the differentiation process by producing new cell-cell attachments, cytoskeleton, etc. The function of stratum granulosum is to produce cornified envelope proteins, crosslinking enzymes and appropriate lipids. In stratum corneum, the proteins are crosslinked, lipids are extruded, and the epidermal barrier is formed. In healthy epidermis, only the basal layer contains mitotic cells. The keratinocytes of the basal layer are characterized by tightly packed palisaded appearance and a position immediately attached to the basal lamina. The cytoskeleton in basal layer contains keratins K5, K14, and small amounts of K15; their function is to anchor the epidermis firmly to its substratum (66). These keratins attach to hemidesmosomes through ␣6␤4 integrins, the laminin-V receptors. Additional integrins, ␣2␤1, ␣3␤1, and ␣5␤1, also bind to the matrix proteins of the basal lamina and provide the basal keratinocytes with essential information about their physical location. Detachment of integrins from their ligands may play a role in initiating differentiation (1, 2). Junction proteins in the basal keratinocytes include hemidesmosomal proteins, as well as desmogleins 2 and 3, desmocollin 3, E- and P-cadherins, and ␣1-connexin (34, 40, 44, 75). Many receptors have been demonstrated on the cell surface of the basal keratinocytes, providing responsiveness to the signals from surrounding cells (37, 41, 50, 51, 59, 64). In addition, basal keratinocytes can produce a large variety of secreted peptides, growth factors, cytokines, and chemokines (43, 47, 62). Thus the basal layer is not just a passive responder to signals from other cell types but can also be an initiator and modulator of the signals in the surrounding tissue. The keratinocytes in suprabasal layers, spinosum, granulosum, and corneum, do not divide under normal conditions, but can be induced to do so, e.g., in wound healing (29). Their cytoskeleton contains keratins K1 and K10 and in upper layers K2e, in addition to the remnants of K5 and K14 from the basal layer. A new set of junctional proteins is produced, including desmocollin 1, plakophilin 1, and ␤2-connexin, along with desmogleins 1 and 3, E-cadherin, and occludin (34, 40, 44, 75), and so are the early differentiation markers, such as involucrin (19, 77). The cytoskeletal changes probably reflect the differences between the essentially immobile basal layer and the easily peripatetic suprabasal cells. These cells can also respond to a wide range of stimuli, but their responses are fundamentally different from the basal layer responses. For example, the

1094-8341/06 $8.00 Copyright © 2006 the American Physiological Society

65

66

PROFILING EPIDERMAL LAYERS

␣5␤1 integrin is not a functional signal receiver in the suprabasal layers (2). The response to IL-1 represents another clear distinction between the basal and suprabasal keratinocytes: the former do not respond to IL-1, the latter, however, do. Suprabasal keratinocytes contain a functional IL-1R, whereas the basal layer expresses the “decoy,” inactive type II receptor (43). Interestingly, both respond to TNF-␣, another proinflammatory signal (39, 41). A comprehensive examination of the signal transducing proteins and transcription factors in the epidermal layers has not been performed. Although some of the proteins that characterize the epidermal layers have been described, a systematic analysis of layer-specific epidermal markers was not feasible because the necessary technology has not been available. An exciting technology to arise from the genome sequencing projects is DNA microarray hybridization, which provides a global view into changes of expression for a large set of genes. Using this methodology, we have presented systematic analyses of gene expression regulated by UV, interferon (IFN)-␥, and TNF-␣ in epidermal keratinocytes (4 – 6, 46). To define comprehensively the genes expressed in different layers of human epidermis using microarrays, we have developed a new and simple procedure that yields large quantities of basal and suprabasal cells from human epidermis, based on the basal-specific expression of integrin-␤4. This, coupled with the microarray technology, enabled us to establish the “transcriptome,” the comprehensive databank of genes expressed in the basal and suprabasal layers of the epidermis. We believe that these will open the door to much additional exploration of epidermal biology and pathology and will not only lead to a new, significantly deeper and more detailed understanding of the epidermis but also provide a rational base for more focused and effective therapies for skin diseases in the future. Furthermore, these studies pioneer analysis of the molecular processes of differentiation, a central problem of current molecular and cell biology. MATERIALS AND METHODS

Keratinocyte culture. Normal epidermal keratinocytes from human foreskin were obtained from Dr. M. Simon (Living Skin Bank, Burn Unit SUNY Stony Brook, Stony Brook, NY). The cultures were initiated using 3T3 feeder layers as described (68) and then frozen in liquid N2 until used. Once thawed, the keratinocytes were grown without feeder cells in defined serum-free keratinocyte growth medium (KGM) supplemented with 0.05 mg/ml bovine pituitary extract, 5 ng/ml epidermal growth factor, and 1% penicillin-streptomycin (KGM from GIBCO-BRL) at 37°C, in 5% CO2. The medium was replaced every 2 days. The cells were expanded through three passages for the experiments and trypsinized with 0.025% trypsin, which was neutralized with 0.5 mg/ml of trypsin inhibitor. We avoid using serum because it can promote partial keratinocyte differentiation. For all experiments, third-passage keratinocytes were used 1 day after reaching confluence. Isolation of ␤4⫹ and ␤⫺ keratinocytes from skin. Skin, discarded after reduction mammoplasty, was obtained according to the Institutional Review Board-approved protocol, generally within 2– 6 h after surgery. Three separate donor specimens were used, harvested, and processed separately at different times. The skin was first washed six times with PBS, and excess liquid was drained. Using scissors and a scalpel, we removed as much fat and dermis as possible. The tissue was cut into 3-mm strips and incubated with dispase (2.4 U/ml, Roche) and RNase inhibitor (4 U/ml Roche) at 4°C overnight. Next day, the epidermis was gently separated from the dermis with forceps Physiol Genomics • VOL

27 •

and incubated in 0.05% trypsin, 0.02% EDTA (GIBCO-BRL) at 37°C. After 10 min, 2 volumes of 0.5 mg/ml trypsin inhibitor (Sigma) were added and the tissue filtered through Cell Strainer (Falcon). The trypsinization of the tissue was repeated twice more. The cells were collected by centrifugation, examined using trypan blue, and counted, and, if appropriate, the isolates were combined. This represented the unfractionated, total epidermal cell population (E). Magnetic beads, M-450 rat anti-mouse IgG1, were prepared as suggested by the manufacturer (Dynal). The cells were incubated with the beads in the following ratio: 100 ␮l beads: 10 –20 ␮g ␤4⫺ antibody: 4 ⫻ 106 cells (exactly) in 1⫻ PBS, 0.1% BSA, at 4° for 1–2 h. We used M-450 rat anti-mouse IgG1 beads and the 3E1 clone ␤4⫺ antibody from GIBCO-BRL. The beads were separated on a magnetic separator for 2–3 min and washed three or four times with PBS, collecting and combining the nonadherent, ␤4⫺ cells as the suprabasal cell population (S). The beads bound to the ␤4⫹ cells, the basal cells (B), were used in RNA isolation without removing the cells from the beads. Isolation of total RNA. To obtain RNA of appropriate quality for chip analysis from in vivo epidermis, we have tested several purification methods. After extensive experimentation, we settled on the following approach. First, the epidermal cells are disrupted, and the RNA is isolated with TRIzol (GIBCO). This is followed by the use of Qiashredders to homogenize cell extracts with centrifugation at 1,800 g for 2 min. DNA is removed with on-column DNase digestion using Qiagen RNases-free DNase Set. RNeasy kits from Qiagen are used to prepare the RNA according to the manufacturer’s protocols. TRIzol gives good yields and effectively disrupts the epidermis, but the purity of the RNA is inadequate for the subsequent steps; the RNA isolation kit gives adequate purity but inefficiently disrupts the tissue, which is why the two are used in series. All solutions, starting with the one containing dispase, up to the final RNA elution, contain RNase inhibitor, because the skin is particularly rich in RNases, necessitating their inhibition. Five micrograms of total RNA are reverse transcribed, amplified, and labeled as described (30). Labeled cRNA is hybridized to oligonucleotide microarrays from Affymetrix. Arrays are washed, stained with antibiotin streptavidin/phycoerythrin-labeled antibody, and scanned using the GeneChip system (Hewlett-Packard) and GeneChip 3.0 software to determine the expression of each gene. Northern and Western blot analyses. Approximately 5 ␮g of each RNA sample were loaded onto 1.5% agarose-formaldehyde gel. The RNA was transferred to a nylon membrane (Amersham) and crosslinked in a UVC Stratalinker (Stratagene). The Notch probe was synthesized using keratinocyte RNA and the RT-PCR kit (Promega) with the following primers GGTGATCGGCTCGGTAGTAA and ACACCCAGGACGAGAGATGAC; the transglutaminase 1 probe with GCTTAGCTACCTACGCACGG and GCCACTGCTAGTCTCTTGGG; and the GAPDH probe with ACAGTCAGCCGCATCTTCTT and GACAAGCTTCCCGTTCTCAG. The 32P-labeled probes were generated with [32P]dCTP (3,000 Ci/mmol, Dupont NEN) and Random Primed DNA Labeling Kit (Roche Molecular Biochemicals) and purified using D-Salt Excellulose Plastic Desalting Columns (Pierce) from the inserts of plasmids containing cloned cDNAs of K15, K14, and hypoxanthine phosphoribosyltransferase, as previously described (38, 48, 67). The probes were hybridized using ExpressHyb solution (Clontech) at 68°C for 1 h. The membrane was washed with 2⫻ SSC, 0.05% SDS solution, with continuous shaking three times for 30 min at room temperature and with 0.1⫻ SSC, 0.1% SDS at 50°C for 40 min. The membrane was exposed to BioMax film (Kodak) at ⫺80°C. The radiographs were scanned and analyzed with MultiImagine Light Cabinet AlphaImagine 2000 documentation and analysis system (Alpha Innotech). The RNA levels in Northern blots were normalized to GAPDH mRNA. For preparation of the whole cell lysates, we washed the cells with cold PBS and broke them in radioimmunoprecipitation assay buffer containing 50 mM Tris 䡠 HCl pH 7.5, 150 mM NaCl, 2 mM EDTA, 1% Nonidet P-40, 0.1% SDS, 25 mM NaF, 1 mM PMSF, 10 ␮g/ml www.physiolgenomics.org

PROFILING EPIDERMAL LAYERS

67

was performed using algorithms available at http://rana.stanford.edu/ software (21). Data from three individual donors were collected and analyzed both individually for each donor, as well as in aggregate for all three. The statistical analysis of the data included t-test and Mann-Whitney statistical evaluations. We calculated the medians of expression values for the each gene in the ␤4⫹, ␤4⫺, and total epidermal cells, determined the expression ratios between ␤4⫹ and ␤4⫺ cells, as well as between ␤4⫹ and total epidermal cells from the same donor, and then averaged these values for all three donors. We selected for further analysis the genes with at least one statistically significant difference in the t-test or Mann-Whitney analysis, and at least a twofold average difference of expression between either the ␤4⫹ vs. the ␤4⫺, or the ␤4⫹ vs. the total epidermal cells. We developed an extensive gene annotation table describing the molecular function and biological category of the genes present on the

Fig. 1. Clean separation of the basal layer keratinocytes demonstrated on Western blots. Separated preparations of ␤4⫹ and ␤4⫺ cells were tested with AE2 and AE3 antibodies, generous gifts from H. Sun. The same blot was stripped and reprobed with AE2. On the right, keratin proteins are identified. HDJ was used for loading control after stripping the involucrin antibody. Equivalent results were obtained in Northern blots using K5, K14, and K15 probes (not shown). Two independent separations from 2 individual donors gave identical results. Note that K5, basal keratin, is a stable protein that persists at the protein level in the suprabasal, ␤4⫺ cells.

aprotinin, and 10 ␮g/ml leupeptin. The lysates were centrifuged at 15,000 g, 10 min at 4°C. The protein concentration of each sample was determined using a Bio-Rad Protein assay reagent. Protein (50 ␮g) was loaded on 10% or 20% SDS-polyacrylamide gels, transferred to polyvinylidene difluoride membrane (Immobilon-P) using a semidry transfer cell (Bio-Rad), and blocked in 5% BSA in TBST (50 mM Tris, pH 7.5, 150 mM NaCl, 0.05% Tween 20). The membrane was incubated with primary antibody overnight at 4°C in 2.0% BSA in TBST, washed extensively, incubated with anti-rabbit or anti-mouse horseradish peroxidase-conjugated secondary antibodies, and visualized with the ECL detection kit (Amersham). The equivalent loading of proteins in each well was ascertained by Ponceau staining of the membrane. The primary antibodies were a gift from T.-T. Sun (New York University School of Medicine) or were purchased from Transduction Labs, Leinco Technologies, and Research Diagnostics. Array data analysis. Generally, we used the same data analysis approach as before (5, 6). Intensity values from HUG95Av2 chips were obtained with Microarray Suite version 5.0 (Affymetrix) and scaled by calculating the overall signal for each array. To compare data from multiple arrays, the signal of each probe array was scaled to the same target intensity value. To eliminate genes exhibiting potential false positive differential expression, we select only genes that possessed relative signal intensity greater than one standard deviation above average for all genes scored as “absent” in the samples. Differential expressions of transcripts were determined by calculating the fold change between the signal intensity values from the B, S, and E cells using GeneChip data mining software (Affymetrix-DMT3.0). To improve reliability, we checked individually the absolute expression levels and P values in all samples. The hierarchical clustering Physiol Genomics • VOL

27 •

Fig. 2. Known differentially expressed functional categories: hemidesmosomal and cornified envelope components, cell cycle proteins and DNA replication/repair enzymes. The genes preferentially expressed in the basal layer (B) are marked in red and have positive values for fold-differential expression; those overexpressed in the suprabasal layers (S) are marked with negative values and green color. The ␤4 integrin, used to separate the layers, is highlighted in yellow. Note the predominantly suprabasal expression of the cell cycle inhibitors and basal expression of cell cycle and DNA replication genes. Also note the unexpected basal expression of S100A genes. Column E: the mean expression levels in total epidermis from 3 individuals. Columns t: the results of statistical evaluations using the t-test; columns M: a nonparametric, Mann-Whitney test showing differential expression between the ␤4⫹ and ␤4⫺ cells, B/S, or between the ␤4 and total epidermal preparations, B/E. Identical notations are used in Figs. 2 through 10. www.physiolgenomics.org

68

PROFILING EPIDERMAL LAYERS

chip (T. Banno and M. Blumenberg, unpublished). The table is primarily based on the data by J. M. Rouillard (70) and the Gene Ontology Consortium data http://dot.ped.med.umich.edu:2000/ ourimage/pub/shared/JMR_pub_affyannot.html and http://cgap. nci.nih.gov/Genes/GOBrowser. The regulated genes were classified according to this table. In addition, we used the L2L program to identify biological processes and cell components statistically overrepresented in our lists of differentially expressed genes (61). The transcription factor binding sites in the regulated genes were identified with the oPOSSUM program (32). We first calibrated the parameters of the program using a set of identified NF-␬B-regulated genes (5) to obtain the optimal statistical P values in the one-tailed Fisher exact probability analysis. We then used the same parameters for the differentially expressed genes. The comparison of the lists of differentially expressed genes in this study with those genes found to be regulated in epidermal keratinocytes by TNF-␣ or IFN-␥ was performed using the LOLA set of programs (11). Immunohistochemical staining of normal human skin. Pieces of normal human skin were obtained immediately after surgery, mounted in tissue Tec optimum cutting temperature compound (Sakura Finetek), and immediately frozen in liquid nitrogen. Sections, 4- to 6-␮m thick, were obtained with a cryostat (Miles Laboratories). Slides containing frozen sections were dried and fixed in 100% acetone or 4% paraformaldehyde for 20 min at ⫺20°C. After being blocked with 2% BSA in PBS, they were incubated with primary antibodies (Santa Cruz) at 4°C overnight. The sections were washed with PBS three times and incubated with secondary fluorescent-conjugated antimouse IgG antibody diluted 1:60 (Sigma-Aldrich) for 1 h at room temperature. After the final wash in PBS, they were mounted with mounting media (Fluoromount-G, Southern Biotechnology Associates) and covered with coverslips. As a negative control, untreated skin was stained omitting the incubation with primary antibody. Stained sections were examined under a Carl Zeiss microscope, and digital images were collected with the Adobe Photoshop TWAIN_32 program. RESULTS AND DISCUSSION

Separation of basal and suprabasal keratinocytes from human epidermis. We have developed a method for purification of basal layer keratinocytes in sufficient quantities and of appropriate quality for microarray analysis. The method consist of the following steps: 1) separation of the epidermis from dermis using dispase; 2) meticulous, fastidious, and gentle trypsinization of the epidermis to obtain undamaged single cells in suspension; 3) tagging the basal layer cells with ␤4 integrin antibody; and 4) separation of tagged cells using magnetic beads. We have considered and tried several approaches, with less success. For example, we found that fluorescence-activated cell sorting takes longer and damages cells more that the magnetic separation. We also found that ␣5 and ␤1 integrins are present at sufficient levels on suprabasal cells to preclude a clean separation from the basal layer (not shown). We have also painstakingly optimized the times, concentrations, and mechanical procedures of the isolation. With this optimized procedure, we can prepare easily and efficiently 5–10 ⫻ 107 basal layer cells and 2–5 ⫻ 108 nonbasal layer cells in a single experiment. The quantities are amply sufficient for our purposes. The separation and purification of the epidermal layers were performed three separate times with samples from three separate donors. Reduction mammoplasty was chosen as the ideal source of skin because of the relatively narrow age range of donors (20s and early 30s), same sex, same body area, and a Physiol Genomics • VOL

27 •

-

Fig. 3. Adhesion proteins, extracellular matrix components, proteases, and their inhibitors. Integrins and ECM genes are produced by the basal cells, whereas the desmosomal genes are specific for the suprabasal ones. Different layers make different categories of junctional proteins, proteases, and their inhibitors.

uniformly sun-protected site. The uniformity of the starting material can ensure the reliability and consistency of the results. However, the transcription profiles of skin from men, other body sites, age, and different sun exposure will probably turn out to have subtle but significant differences from those presented here. www.physiolgenomics.org

PROFILING EPIDERMAL LAYERS

To ascertain the purity of the separated cell populations, we used Western and Northern blots. As shown in Fig. 1, AE3 monoclonal antibody, specific for type II keratins (15, 60), detects K1, K2, K5, and K6 keratins in total epidermal cells and in the ␤4⫺ cells, whereas only the basal layer-specific K5 was found in the cultured and ␤4⫹ cells. AE2, which recognizes differentiation-specific keratins (16, 78), detects K1 and K10 only in total epidermis and ␤4⫺ cells, but not in cultured and ␤4⫹ cells. Similarly, the differentiation markers filaggrin and involucrin were found only in the ␤4⫺, but not in ␤4⫹, cells. These results demonstrate that the preparations of ␤4⫹ and ␤4⫺ cells contain the appropriate markers that correspond to basal layer and nonbasal layer keratinocytes. Microarrays of separated epidermal ␤4⫹ and ␤4⫺ cells. The reliability of microarray analysis is considerably enhanced by multiple, redundant experiments, and these are commonly performed in triplicate. Therefore, we performed three separations of ␤4⫹ and ␤4⫺ cells, from three different donors at three separate occasions, and tested each sample as in Fig. 1A. From all three donors, we hybridized the ␤4⫹ and ␤4⫺ cells, as well as the total epidermal cell preparations before the magnetic separation, a total of nine microarray chips. Because of the interpersonal variation of skin samples (M. Blumenberg, unpublished), the data were statistically analyzed in two ways: collectively and individually. Collectively, we calculated the medians of the ␤4⫹, ␤4⫺, and total epidermal cells and performed t-test and Mann-Whitney statistical evaluations. Individually, we calculated the expression ratios between ␤4⫹ and ␤4⫺ cells, and between ␤4⫹ and total epidermal cells from the same donor and then averaged these values. We selected the genes with at least one statistically significant difference in the cumulative analysis and at least a twofold average difference of expression between the ␤4⫹ cells and either the ␤4⫺ or the total epidermal cells in the individual analysis. This selection identified 875 differentially expressed genes. The complete data are submitted to Gene Expression Omnibus repository (in process). In these samples, the keratinocytes vastly outnumber all other cell types combined, by five- to tenfold. We do observe a low level of expression of melanocyte-specific markers, see below, but these do not affect significantly our results and conclusions. Figure 2 contains three functional categories of genes known to be differentially expressed, and it serves as an independent control of our separation of the ␤4⫹ and ␤4⫺ cells. The categories are hemidesmosomal genes, which are specific for basal cell and of which ␤4 was used for cell separation, cell cycle, and DNA replication genes, which are also specific for the basal, proliferative compartment of the epidermis and the epidermal differentiation and cornification markers, which are specific for the suprabasal, differentiating cells. As expected, the hemidesmosomal genes are exclusively expressed in the ␤4⫹ cells. We note that the “fold difference” for ␤4 integrin, the cell surface marker used to separate the ␤4⫹ and ␤4⫺ cells, is 5.0. The well-recognized epidermal differentiation markers, filaggrin, involucrin, loricrin, and SPRRs, are overexpressed in the ␤4⫺ cells. The S100A calcium-binding proteins are also considered to be epidermal differentiation markers because their genes bracket the epidermal differentiation complex (EDC) on chromosome 1, which contains the SPRR and LEP gene families and the filaggrin/involucrin/loricrin Physiol Genomics • VOL

27 •

69

genes (22, 52, 55). Interestingly, here we show that the S100A genes are expressed in the basal cells, as has been reported previously (10, 22, 58, 69, 73). This suggests that the regulation of the EDC genes is more complex that believed hitherto. Only basal cells replicate in healthy epidermis, and, accordingly, most cell cycle genes are found expressed in the ␤4⫹ cells, whereas the cell cycle inhibitors are in the ␤4⫺ group. Similarly, the DNA replication enzymes are also in the ␤4⫹ cells; the two exceptions, POLB and NP, found expressed at higher levels in the ␤4⫺, play roles primarily in DNA repair and degradation. Taken together, the data in Fig. 2 agree very well with the expression studies of individual proteins (8, 20, 27) and explicitly validate our successful separation of epidermal keratinocytes from particular layers. This gave us confidence to examine several additional important and interesting functional groups. Adhesion, ECM, and proteolysis. The basal and suprabasal cells express clearly distinct adhesion genes (Fig. 3). Desmosomal genes are expressed at higher levels in ␤4⫺ cells, whereas all integrins are expressed in the ␤4⫹ cells. Interestingly, different adhesion markers are expressed in the different layers: the cadherins are in ␤4⫹ cells, conferring selective homotypic adhesion to the basal keratinocytes, whereas the

Fig. 4. Metabolic enzymes, oncogenes, and tumor suppressors. The suprabasal cells make the lipid and steroid metabolism enzymes. Basal cells make solute carrier gradient-driven transporters, whereas the suprabasal cells use ATP for transport. Note the oncogenes expressed in the suprabasal and tumor suppressors in the basal cells. www.physiolgenomics.org

70

PROFILING EPIDERMAL LAYERS

ases, each cell type produces a characteristic set of serine proteinase inhibitors (SERPINs). The genes encoding extracellular proteins listed in Fig. 3 reflect the different milieus the two cell populations come into contact with. Although some of the markers were known to be differentially expressed in epidermal layers, we were surprised by several, e.g., by the extensive and diverse contribution to the basement membrane produced by the epidermis, by the specificity of the SERPINs produced by the different layers, and by the presence of cadherin 13 in the basal layer. These novel findings represent major advances in the understanding of epidermal structure and function and will inspire focused research in these areas. Metabolism. Basal and differentiating layers of the epidermis have clearly different metabolic needs and processes and express different metabolic enzymes. In Fig. 4, we chose four

Fig. 5. Cell surface receptors and secreted signaling proteins. Several pairs of receptor/ligand combinations transmit signals between the basal and suprabasal cells, e.g., FGF receptor (FGFR) 2/FGF2 and ERBB2/NRG2.

tight- and gap-junction genes are expressed in ␤4⫺ cells, contributing to the epidermal barrier. The ECM proteins are predominantly made by the ␤4⫹ cells, as expected because they abut the basement membrane. Laminin-V components (the ligand for ␤4 integrin), as well as collagen IV (a basement membrane component), and collagen VI (a matrix assembly scaffolding protein) are found expressed in the basal keratinocytes. These may play a role in maintenance of epidermal stem cells (57, 76) However, collagen I, which constitutes ⬃80% of the dermal collagen, is absent, which leads us to believe that the contamination of fibroblasts in our preparation is negligible. The ␤4⫹ and ␤4⫺ cells produce distinct sets of proteases and their inhibitors (Fig. 3). Matrix metalloproteases MMP9 and MMP12 are produced by the ␤4⫹ cells, whereas the kallikreins are made by the ␤4⫺ cells, reflecting different extracellular proteolytic targets of the two cells. While only the ␤4⫹ cells need and produce tissue inhibitor of metalloproteinPhysiol Genomics • VOL

27 •

Fig. 6. Immunofluorescence detection of NOTCH3 and IGFBP3 expression in human skin in vivo confirms the microarray results. Antibodies specific for NOTCH3 and IGFBP3 stain the suprabasal and basal layer, respectively, in 5-␮M sections of human epidermis. The dashed line marks the basement membrane, indicating the epidermis-dermis boundary. Control without primary antibody shows weak background fluorescence of the stratum corneum. www.physiolgenomics.org

PROFILING EPIDERMAL LAYERS

71

Fig. 7. Signal transduction. The right-most column contains the designations of the signaling pathways, if known. Note the basal prevalence of Ras genes and the suprabasal prevalence of diacylglycerol (DAG), phosphatidylinositol (IP), and prostaglandin (PG) pathways.

specific categories because of their special significance. The epidermal permeability barrier is formed by the differentiating cells in the granular layer and consists of a protein and a lipid component (23). Of the known genes involved in lipid metabolism, at least 11 are produced at higher levels in the ␤4⫺ cells. The appropriate steroid metabolism in the epidermis is also essential for the proper function of the cornified layer, evidenced by the fact that steroid sulfatase deficiency leads to X-linked ichthyosis (18). Consequently, we find six enzymes of the steroid metabolism preferentially expressed in the ␤4⫺, differentiating cells (Fig. 4). Thus, related to their barrier function, the suprabasal cells contain the enzymes for steroids and other lipids metabolism, which means that the basal cells Physiol Genomics • VOL

27 •

only passively, without biochemical processing, transport the lipid barrier building blocks. An interesting difference between the ␤4⫹ and ␤4⫺ cells is found in the transporter protein genes: the ␤4⫹ cells express several solute carriers, which exchange ions with the extracellular medium, whereas the ␤4⫺ cells express ATPase family of transporters, which involve energy metabolism. Such characteristically different approaches to transport have not been described before and their further exploration may yield profound new understanding of the function of different epidermal layers. We were particularly excited by the oncogenesis-related genes are differentially expressed in the ␤4⫹ and ␤4⫺ cells. www.physiolgenomics.org

72

PROFILING EPIDERMAL LAYERS

The differentially expressed oncogenes were found only in the ␤4⫺, whereas the tumor suppressors were in the ␤4⫹ cells. This may reflect the fact that only the basal cells can proliferate and are therefore more at risk for being transformed. We note that the list in the Fig. 4 is an abbreviated one; the oncogenes and tumor suppressors that are transcription factors or are involved in the cell cycle are presented in the respective figures. Cell surface receptors and secreted signaling proteins. One of the essential protective functions of the epidermis is to alert the underlying tissues to the damage from the environment, as well as to respond to the stimuli from these tissues. In this as well, the epidermal layers have specific functions. Importantly, we find that both the basal and the suprabasal cells express genes for many cell surface receptors, as well as secreted signaling proteins, which means that they both respond, as well as signal to the surrounding tissue. The ␤4⫹ cells specifically express the decoy IL-1R to protect themselves from the IL-1 signal, and the FGF receptor (FGFR), which can respond to the proliferative signal of the keratinocyte growth factor and other FGFs (Fig. 5) (43, 49). In contrast, the ␤4⫺ cells express the IL-1R-associated kinase, as well as receptors for thrombin, ephrin, etc. Interestingly, we find that the suprabasal cells specifically express the two relatives of the EGF receptor (EGFR), ERBB2 and ERBB3. These have distinctive ligand recognition specificities, and ERBB2 dimerized with EGFR can be activated by neuregulin 2 (14, 24). Neuregulin 2 is differentially expressed by the ␤4⫹ cells (Fig. 5), which creates a possibility for paracrine signaling from the basal to the suprabasal cells. Conversely, we find FGF2 produced by the ␤4⫺ cells, potentially signaling to the FGFR-containing ␤4⫹ cells. Neither neuregulin 2 nor FGF2 was known previously to be expressed by keratinocytes! The ␤4⫹ cells also express a relatively large number of genes encoding secreted proteins, i.e., growth factors and cytokines that target, in addition to epidermal cells, the underlying fibroblasts and endothelial cells, perhaps lymphocytes as well. Prominent among these are the members of the transforming growth factor (TGF)-␤ family. We find that the ␤4⫺ cells selectively express higher levels of apolipoproteins D and E, although the basal cells were shown to secrete ApoE (7). Our data identify two potential paracrine signaling mechanisms between the basal and suprabasal cells, an area of research that has not been intensely addressed in the past. These are the neuregulin/ERBB and FGF2/FGFR pairings. The identification of the receptors in different epidermal layers suggests a possibility for targeted differential therapy of proliferative diseases, both to enhance proliferation, e.g., in chronic wounds, and to inhibit it, e.g., in psoriasis. We were particularly intrigued by the production of neuregulin 2 and its co-receptor ERBB2 in the basal and suprabasal cells, respectively. Their roles in skin cancers have not been fully described, and the discovery of this paracrine pathway correlates with the predominantly suprabasal expression of oncogenes and the basal expression of tumor suppressors. We suspect that the abundance of cell cycle and DNA replication genes makes the basal cells particularly susceptible to oncogenic transformation. The data presented here can form a basis for the analysis of genetic changes necessary for the occurrence of basal and squamous cell carcinomas. Physiol Genomics • VOL

27 •

To confirm that some the microarray findings accurately reflect the differential expression in human epidermis in vivo, we used immunohistology on skin sections (Fig. 6). We tested the expression of IGFBP3, a marker of basal ␤4⫹ cells, and NOTCH3, a marker of suprabasal, ␤4⫺ cells. We find that the in vivo data perfectly correspond to the expectations from the microarray experiments, thus validating our results. Differentially expressed signal transduction proteins and transcription factors. The differences in the signaling pathways in the two cell types are very interesting (Fig. 7). For example, the Ras pathway genes are mainly found expressed in the ␤4⫹ cells, perhaps reflecting their proliferative capability. Playing an important role in proliferation, the Ras family may regulate the cell cycle and DNA replication. Interestingly, the diacylglycerol, inositol phosphate, and prostaglandin pathways are more prominent in the ␤4⫺ cells, perhaps regulating the

Fig. 8. Transcription factors. Basal and suprabasal cells express distinctive sets of transcription factors. The 4 homeobox genes are set apart on the bottom. www.physiolgenomics.org

PROFILING EPIDERMAL LAYERS

lipid metabolism in the process of barrier formation. Various protein kinases and phosphatases are differentially expressed in the different epidermal layers, and we note their abundance in the ␤4⫺ cells, indicating that these cells, during their passage to a cornified layer, are highly responsive to signaling. We note that the ␤4⫹ and ␤4⫺ cells express Jagged 2 and Notch 3, respectively, potentially establishing another juxtacrine signaling pathway. Notch signaling regulates epidermal differentiation (63), although these two particular proteins were not known to be expressed in epidermis. It will be interesting to ascertain whether the Jagged 2 and Notch 3 expressing cells are juxtaposed in the epidermis. Similarly instructive are the differences in the transcription factors expressed in the basal and suprabasal cells, listed in Fig. 8. We do not know enough to assign a function to each of the differentially expressed gene in this functional category, but we believe that the field will be able to fill this gap. The SMAD genes are well represented in the ␤4⫹ fraction, reflecting the specific effects of TGF-␤ on the basal cells (26, 28, 37, 79). We were intrigued by the TGF-␤-responsive machinery specific for the basal cells. This correlates well with the previous results from our and other laboratories, suggesting that TGF-␤ promotes the basal cell phenotype in skin (37). STAT5, responsive to IL-3, prolactin, etc., is also preferentially found in the ␤4⫹ cells. CCAAT enhancer-binding protein/␤ is found in the ␤4⫺ cells, confirming its role in epidermal differentiation (54). NF-␬B proteins also play a role in differentiation (17, 72); interestingly, we find NFKBIE, an inhibitor of NF-␬B, in the ␤4⫹ and TANK, an activator, in the ␤4⫺ cells. Several IFN- and TNF-␣-responsive transcription factors are present in the ␤4⫺ population, demonstrating that the suprabasal cells are clearly competent to respond specifically to proinflammatory and immunomodulatory signals.

73

The ␤4⫹ and ␤4⫺ cells contain characteristic sets of zinc finger, Ets proteins, and nuclear receptors with associated co-regulators. Interestingly, we do not find significant differences in retinoic acid and vitamin D receptors, known regulators of epidermal differentiation. In addition, we do not find basonuclin differentially expressed, presumably because significant levels of this protein are made both in the early suprabasal keratinocytes and in the basal cells (74). SREBF2, a sterol regulatory factor, is preferentially expressed in the ␤4⫺ cells, reflecting the fact that the steroid metabolism specifically occurs in these cells. We would like to draw attention to KLF4, the kruppel-like factor implicated in regulation of epidermal differentiation and barrier formation (71, 35); KLF4 is overexpressed in the suprabasal, ␤4⫺ cells. We were particularly interested in homeobox-containing transcription factors, because these are thought to regulate epidermal differentiation (56, 65). Indeed, we find DLX5, a distal-less family member specifically expressed in the ␤4⫺ cells (DLX3 is not represented on the microarrays). DLX5 may play an essential role in regulating epidermal differentiation. Parenthetically, IRX5, a member of the Iroquois family, is probably expressed in melanocytes, see below. Melanocyte-derived transcripts. Although keratinocytes are by far the predominant cell type in the epidermis, melanocytes constitute a significant minority; it is therefore not surprising that several melanogenesis enzymes have been detected by the arrays (Fig. 9). Melanocytes do not express ␤4 integrin and are found in the ␤4⫺ population. We noticed a peculiarity in the expression of melanogenesis enzymes: they were expressed at significantly higher levels in the ␤4⫺ population than in the unseparated, total epidermal cells. For a great majority of genes, the expression levels in the ␤4⫺ population and in the total epidermal cells are very similar because the ␤4⫺ cells are

Fig. 9. Melanocyte-specific markers and genes with similar expression pattern. The right-most column contains the ratio of expression in total epidermal, vs. ␤4⫺ cells. Note the universally significant B/S ratios, much less prominent B/E ratios. The well-known melanocyte markers are highlighted in yellow.

Physiol Genomics • VOL

27 •

www.physiolgenomics.org

74

PROFILING EPIDERMAL LAYERS

Fig. 10. Epidermal differentiation markers. The penultimate column contains the ratio of expression in total epidermal vs. ␤4⫺ cells, whereas the right-most column contains the chromosomal location of the genes. The well-characterized epidermal differentiation markers are highlighted in yellow.

the predominant population in the epidermis. The overrepresentation of melanocytes in the ␤4⫺ cells population derives, we believe, from the isolation procedure: being in the basal layer, the melanocytes are quantitatively released by trypsinization, whereas many of the ␤4⫺ keratinocytes, especially those at advanced stages of differentiation, are more resistant to trypsinization, remain in undigested cell aggregates, and are lost from the ␤4⫺ isolates. The skewed ␤4⫺ cells vs. total cells ratio of melanogenic genes led us to ask Physiol Genomics • VOL

27 •

whether there are additional genes with similarly skewed ratios, i.e., “expression signatures.” Using a twofold cutoff for the ␤4⫺/total epidermis expression ratio, we found 49 such genes (Fig. 9). Among these is vimentin, a known melanocyte cytoskeletal component, which validated our approach. We note that for all these genes there is a statistically significant difference in expression between the ␤4⫹ and ␤4⫺ cells, which is not necessarily true for the ␤4⫹ vs. total cells comparisons (Fig. 9, columns 4 –7). We found several addiwww.physiolgenomics.org

PROFILING EPIDERMAL LAYERS

tional genes, encoding cytoskeletal, metabolic, signaling proteins, and transcription factors, as well as nine uncharacterized expressed sequence tags (ESTs), none previously known to be expressed in melanocytes. We are particularly intrigued by the relatively high levels of c-Myc, an oncogene important in melanomas. These data will undoubtedly be useful to studies of melanocytes, which, to date have been largely restricted to melanoma cells. Parenthetically, the chromosomal locations of the melanocyte markers appear random, i.e., we did not detect clusters of melanocyte genes, similar to the epidermal differentiation gene clusters. These genes, presumably, derive from melanocytes, although we cannot at present exclude the possibility that they derive from other cell types that reside in the basal layer, such as Langerhans or Merkel cells. We find no fibroblast markers, such as Coll1, at all in our microarrays, and therefore we believe that the contributions of dermal cells, fibroblasts, endothelial cells, and smooth muscle cells, are below detection levels. Keratinization markers. Although the melanocyte markers are overrepresented in the ␤4⫺ population relative to the total epidermal cells, another group of genes was underrepresented in the ␤4⫺ population relative to the total epidermal cells. Careful examination of the functions of genes with this characteristic expression signature pointed to the known cornified envelope proteins, markers of advanced epidermal differentiation. These genes are, of course, not expressed in the ␤4⫹ cells and, relative to the total epidermis, are selectively missing from the ␤4⫺ population as well. We believe that the cells in the most advanced stages of differentiation are more resistant to trypsinization, remain in clusters, and are lost during the purification of the ␤4⫺ cells, whereas the harsh RNA purification procedure releases the mRNA equally from all cells of the total epidermal preparations, including the cells in the most advanced stages of cornification. Alternatively, the mRNA for markers of advanced differentiation may be selectively degraded during the isolation procedure, targeted by the epidermal RNases. Using a 1.7-fold cutoff for the total epidermis/␤4⫺ expression ratio, we find 83 differentially expressed genes (Fig. 10). Among these are transglutaminases, loricrin, involucrin, and filaggrin, well-known markers of cornification (8, 20, 27). We also find a significant number of lipid and steroid metabolic enzymes, kallikreins, etc. In this group we find ARS, the gene that, when mutated, causes Mal de Maleda, an inherited disorder of cornification (3, 25, 53). In addition, proteases and protease inhibitors are amply represented among these genes, reflecting the extensive remodeling of the stratum corneum on one hand and protection from the environmental microorganisms on the other. Similarly, in this group we find claudin and occludin, the tight junction proteins with important functions in permeability barrier and antimicrobial defense (4, 31, 36). Several interesting regulatory genes fall in this group, encoding protein kinases and phosphatases, and we would like to draw attention to Fyn, which may play a role in regulation of epidermal differentiation (12, 33). Using this expression signature, we identified a much larger set of potential late epidermal differentiation markers, many of which we were not aware of before. The newly identified genes include signaling proteins and transcription factors, potentially regulators of the final steps in keratinocyte terminal differenPhysiol Genomics • VOL

27 •

75

tiation. Interestingly, the location of these genes in the human genome is not random. As mentioned above, many cornified envelope proteins are encoded by the EDC on Ch1q21, and indeed, we find a cluster of seven genes mapping at that locus. When we examined the genomic positions of all genes in Fig. 10, we identified additional chromosomal bands with multiple differentiation specific genes, 6p21, 11q13, and 19q13. We hypothesize that these clusters also have functional significance, similar to the already described cluster on 1q21. We are pursuing these leads in the hope of identifying additional EDCs and associating them with mapped human inherited skin diseases. We would like to point out that, due to space limitations, we have not included in the above analysis genes encoding several important functional categories, e.g., cytoskeletal and membrane proteins, apoptosis, immune response, and RNA metabolism; all these can be found in the Supplement. The online version of this article contains supplemental data. Comparison of the lists of differentially expressed genes with their assigned ontology functions confirms the above analysis (Fig. 11). This is particularly apparent in the

Fig. 11. Ontological categories of differentially expressed genes. Biological processes (left) and cellular components (right) overrepresented among the genes with ␤4⫹, ␤4⫺, upper layer cells and melanocyte signature expression. The P value shows the probability that the specific process or component is as overrepresented in a same-size list of random genes. Only categories with P value ⬍10⫺3 are shown. Note the redundancy among similar categories. No cellular component was found to be overrepresented among the genes with the melanocyte-characteristic expression signature. www.physiolgenomics.org

76

PROFILING EPIDERMAL LAYERS

cellular component categories, where in the basal cells the ECM and basement membrane genes are overrepresented, whereas the intercellular junction components are characteristic of the ␤4⫺ cells and the cornified envelopes in the cells of the “upper layer.” The biological processes characteristic of the ␤4⫹ cells include transport of nutrients and cell motility. In contrast, the ␤4⫺ cells’ epidermis development, keratinocyte differentiation, eicosanoid, lipid, and steroid biosynthesis are overrepresented. These biological processes are even more pronounced in the terminally differentiating, upper layer cells. Melanin biosynthesis and pigmentation are, as expected, greatly overrepresented among the melanocyte signature genes. When we compared the lists of differentially expressed genes in the ␤4⫹ and ␤4⫺ cells with the keratinocyte genes regulated by IFN-␥ or TNF-␣ (3, 5), we found several interesting correlations. Both IFN-␥ and TNF-␣ induced ␤4⫹, basal layer-specific genes ANXA3, C1R, C1S, IFIT3, ITGB4, LAMA3, LAMB3, MMP9, and NFKBIE. These are associated with the basement membrane and the proinflammatory processes. Both IFN-␥ and TNF-␣ induced ␤4⫺-specific genes CST6, GPR109B, PTPN12, and QPCT, and TNF-␣ induced but IFN-␥ suppressed differentiation-specific markers FLG, IVL, and SERPINB3. This result confirms our suggestion that IFN-␥ suppresses epidermal differentiation. Furthermore, IFN-␥, but not TNF-␣, induced ␤4⫺ genes CLDN4 and OCLN, encoding tight junction proteins, presumably to build up tight junctions and prevent paracellular transepidermal conduit of viruses, as we suggested previously (3). The data presented here define the transcriptional changes that occur as the epidermal keratinocytes leave the basal layer and commence differentiation. Specifically, we find the genes encoding cell cycle and DNA synthesis proteins expressed in the basal keratinocytes, whereas the cornified envelope proteins and other differentiation markers are expressed in the suprabasal ones. This is as expected and verifies the separation of the layers using our new methodology. But more than that, the data identify and define the fundamental biological differences between the progenitor and the operational cells of the epidermis. ACKNOWLEDGMENTS Special thanks go to Marjana Tomic-Canic and the members of our laboratory for advice, reagents, and encouragement. GRANTS Our research was supported by National Institute of Arthritis and Musculoskeletal and Skin Diseases Grant AR-41850. REFERENCES 1. Adams JC and Watt FM. Fibronectin inhibits the terminal differentiation of human keratinocytes. Nature 340: 307–309, 1989. 2. Adams JC and Watt FM. Changes in keratinocyte adhesion during terminal differentiation: reduction in fibronectin binding precedes alpha 5 beta 1 integrin loss from the cell surface. Cell 63: 425– 435, 1990. 3. Bakija-Konsuo A, Basta-Juzbasic A, Rudan I, Situm M, NardelliKovacic M, Levanat S, Fischer J, Hohl D, Loncaric D, Seiwert, and Campbell H. Mal de Meleda: genetic haplotype analysis and clinicopathological findings in cases originating from the island of Mljet (Meleda), Croatia. Dermatology 205: 32–39, 2002. 4. Banno T, Adachi M, Mukkamala L, and Blumenberg M. Unique keratinocyte-specific effects of interferon-gamma that protect skin from viruses, identified using transcriptional profiling. Antivir Ther 8: 541–554, 2003. Physiol Genomics • VOL

27 •

5. Banno T, Gazel A, and Blumenberg M. Pathway-specific profiling identifies the NF-kappaB-dependent tumor necrosis factor alpha-regulated genes in epidermal keratinocytes. J Biol Chem 280: 18973–18980, 2005. 6. Banno T, Gazel A, Blumenberg M, Adachi M, and Mukkamala L. Effects of tumor necrosis factor-alpha (TNF-alpha) in epidermal keratinocytes revealed using global transcriptional profiling. J Biol Chem 279: 32633–32642, 2004. 7. Barra RM, Fenjves ES, and Taichman LB. Secretion of apolipoprotein E by basal cells in cultures of epidermal keratinocytes. J Invest Dermatol 102: 61– 66, 1994. 8. Blumenberg M. Keratinocytes: biology and differentiation. In: Cutaneous Medicine and Surgery (vol. 1), edited by Arndt K, Leboit P, Robinson J, and Wintroub B. Philadelphia, PA: WB Saunders, p. 58 –74, 1994. 9. Blumenberg M and Tomic-Canic M. Human epidermal keratinocyte: keratinization processes. In: Formation and Structure of Human Hair (vol. 78), edited by Jolles P, Zahn H, and Hocker H. Basel: Birkhauser Verlag, p. 1–29, 1997. 10. Broome AM, Ryan D, and Eckert RL. S100 protein subcellular localization during epidermal differentiation and psoriasis. J Histochem Cytochem 51: 675– 685, 2003. 11. Cahan P, Ahmad AM, Burke H, Fu S, Lai Y, Florea L, Dharker N, Kobrinski T, Kale P, and McCaffrey TA. List of lists-annotated (LOLA): a database for annotation and comparison of published microarray gene lists. Gene 360: 78 – 82, 2005. 12. Calautti E, Cabodi S, Stein PL, Hatzfeld M, Kedersha N and Paolo Dotto G. Tyrosine phosphorylation and src family kinases control keratinocyte cell-cell adhesion. J Cell Biol 141: 1449 –1465, 1998. 13. Chambers RC, Leoni P, Kaminski N, Laurent GJ, and Heller RA. Global expression profiling of fibroblast responses to transforming growth factor-beta1 reveals the induction of inhibitor of differentiation-1 and provides evidence of smooth muscle cell phenotypic switching. Am J Pathol 162: 533–546, 2003. 14. Citri A, Skaria KB, and Yarden Y. The deaf and the dumb: the biology of ErbB-2 and ErbB-3. Exp Cell Res 284: 54 – 65, 2003. 15. Cooper D, Schermer A, Pruss R, and Sun TT. The use of aIF, AE1, and AE3 monoclonal antibodies for the identification and classification of mammalian epithelial keratins. Differentiation 28: 30 –35, 1984. 16. Cooper D and Sun TT. Monoclonal antibody analysis of bovine epithelial keratins. Specific pairs as defined by coexpression. J Biol Chem 261: 4646 – 4654, 1986. 17. Dajee M, Lazarov M, Zhang JY, Cai T, Green CL, Russell AJ, Marinkovich MP, Tao S, Lin Q, Kubo Y, and Khavari PA. NF-kappaB blockade and oncogenic Ras trigger invasive human epidermal neoplasia. Nature 421: 639 – 643, 2003. 18. DiGiovanna JJ and Robinson-Bostom L. Ichthyosis: etiology, diagnosis, and management. Am J Clin Dermatol 4: 81–95, 2003. 19. Eckert RL. Structure, function, and differentiation of the keratinocyte. Physiol Rev 69: 1316 –1346, 1989. 20. Eckert RL, Crish JF, and Robinson NA. The epidermal keratinocyte as a model for the study of gene regulation and cell differentiation. Physiol Rev 77: 397– 424, 1997. 21. Eisen MB, Spellman PT, Brown PO, and Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA 95: 14863–14868, 1998. 22. Elder JT and Zhao X. Evidence for local control of gene expression in the epidermal differentiation complex. Exp Dermatol 11: 406 – 412, 2002. 23. Elias PM. Stratum corneum defensive functions: an integrated view. J Invest Dermatol 125: 183–200, 2005. 24. Falls DL. Neuregulins: functions, forms, and signaling strategies. Exp Cell Res 284: 14 –30, 2003. 25. Fischer J, Bouadjar B, Heilig R, Huber M, Lefevre C, Jobard F, Macari F, Bakija-Konsuo A, Ait-Belkacem F, Weissenbach, Lathrop M, Hohl D, and Prud’homme JF. Mutations in the gene encoding SLURP-1 in Mal de Meleda. Hum Mol Genet 10: 875– 880, 2001. 26. Freedberg IM, Tomic-Canic M, Komine M, and Blumenberg M. Keratins and the keratinocyte activation cycle. J Invest Dermatol 116: 633– 640, 2001. 27. Fuchs E. Epidermal differentiation: The bare essentials. J Cell Biol 111: 2807–2814, 1990. 28. Gailit J, Welch MP, and Clark RA. TGF-beta 1 stimulates expression of keratinocyte integrins during re-epithelialization of cutaneous wounds. J Invest Dermatol 103: 221–226, 1994. www.physiolgenomics.org

PROFILING EPIDERMAL LAYERS 29. Garlick JA and Taichman LB. Fate of human keratinocytes during reepithelialization in an organotypic culture model. Lab Invest 70: 916 – 924, 1994. 30. Gazel A, Ramphal P, Rosdy M, De Wever B, Tornier C, Hosein N, Lee B, Tomic-Canic M, and Blumenberg M. Transcriptional profiling of epidermal keratinocytes: comparison of genes expressed in skin, cultured keratinocytes, and reconstituted epidermis, using large DNA microarrays. J Invest Dermatol 121: 1459 –1468, 2003. 31. Herard AL, Zahm JM, Pierrot D, Hinnrasky J, Fuchey C, and Puchelle E. Epithelial barrier integrity during in vitro wound repair of the airway epithelium. Am J Respir Cell Mol Biol 15: 624 – 632, 1996. 32. Ho Sui SJ, Mortimer JR, Arenillas DJ, Brumm J, Walsh CJ, Kennedy BP, and Wasserman WW. oPOSSUM: identification of over-represented transcription factor binding sites in co-expressed genes. Nucleic Acids Res 33: 3154 –3164. 2005. 33. Ilic D, Kanazawa S, Nishizumi H, Aizawa S, Kuroki T, Mori S, and Yamamoto T. Skin abnormality in aged fyn-/- fak⫹/- mice. Carcinogenesis 18: 1473–1476, 1997. 34. Jamora C and Fuchs E. Intercellular adhesion, signalling and the cytoskeleton. Nat Cell Biol 4: E101–E108, 2002. 35. Jaubert J, Cheng J, and Segre JA. Ectopic expression of kruppel like factor 4 (Klf4) accelerates formation of the epidermal permeability barrier. Development 130: 2767–2777, 2003. 36. Jepson MA, Collares-Buzato CB, Clark MA, Hirst BH, and Simmons NL. Rapid disruption of epithelial barrier function by Salmonella typhimurium is associated with structural modification of intercellular junctions. Infect Immun 63: 356 –359, 1995. 37. Jiang CK, Tomic-Canic M, Lucas DJ, Simon M, and Blumenberg M. TGF beta promotes the basal phenotype of epidermal keratinocytes: transcriptional induction of K#5 and K#14 keratin genes. Growth Factors 12: 87–97, 1995. 38. Komine M, Rao LS, Freedberg IM, Simon M, Milisavljevic V, and Blumenberg M. Interleukin-1 induces transcription of keratin K6 in human epidermal keratinocytes. J Invest Dermatol 116: 330 –338, 2001. 39. Komine M, Rao LS, Kaneko T, Tomic-Canic M, Tamaki K, Freedberg IM, and Blumenberg M. Inflammatory versus proliferative processes in epidermis. Tumor necrosis factor alpha induces K6b keratin synthesis through a transcriptional complex containing NFkappa B and C/EBPbeta. J Biol Chem 275: 32077–32088, 2000. 40. Kowalczyk AP, Borgwardt JE, and Green KJ. Analysis of desmosomal cadherin-adhesive function and stoichiometry of desmosomal cadherinplakoglobin complexes. J Invest Dermatol 107: 293–300, 1996. 41. Kristensen M, Chu CQ, Eedy DJ, Feldmann M, Brennan FM, and Breathnach SM. Localization of tumour necrosis factor-alpha (TNFalpha) and its receptors in normal and psoriatic skin: epidermal cells express the 55-kD but not the 75-kD TNF receptor. Clin Exp Immunol 94: 354 –362, 1993. 42. Kupper TS. The activated keratinocyte: a model for inducible cytokine production by non-bone marrow-derived cells in cutaneous inflammatory and immune responses. J Invest Dermatol 94: 146S–150S, 1990. 43. Kupper TS and Groves RW. The interleukin-1 axis and cutaneous inflammation. J Invest Dermatol 105: 62S– 66S, 1995. 44. Langbein L, Grund C, Kuhn C, Praetzel S, Kartenbeck J, Brandner JM, Moll I, and Franke WW. Tight junctions and compositionally related junctional structures in mammalian stratified epithelia and cell cultures derived therefrom. Eur J Cell Biol 81: 419 – 435, 2002. 45. Le Naour F, Hohenkirk L, Grolleau A, Misek DE, Lescure P, Geiger JD, Hanash S, and Beretta L. Profiling changes in gene expression during differentiation and maturation of monocyte-derived dendritic cells using both oligonucleotide microarrays and proteomics. J Biol Chem 276: 17920 –17931, 2001. 46. Li D, Turi TG, Schuck A, Freedberg IM, Khitrov G, and Blumenberg M. Rays and arrays: the transcriptional program in the response of human epidermal keratinocytes to UVB illumination. FASEB J 15: 2533–2535, 2001. 47. Luger TA, Schwarz T, Krutmann J, Kock A, Urbanski A, and Kirnbauer R. Cytokines and the skin. Curr Probl Dermatol 19: 35– 49, 1990. 48. Ma S, Rao L, Freedberg IM, and Blumenberg M. Transcriptional control of K5, K6, K14, and K17 keratin genes by AP-1 and NF-kappaB family members. Gene Expr 6: 361–370, 1997. 49. Maas-Szabowski N, Shimotoyodome A, and Fusenig NE. Keratinocyte growth regulation in fibroblast cocultures via a double paracrine mechanism. J Cell Sci 112: 1843–1853, 1999. Physiol Genomics • VOL

27 •

77

50. Mansbridge JN and Hanawalt PC. Role of transforming growth factor beta in the maturation of human epidermal keratinocytes. J Invest Dermatol 90: 336 –341, 1988. 51. Marchese C, Chedid M, Dirsch OR, Csaky KG, Santanelli F, Latini C, LaRochelle WJ, Torrisi MR, and Aaronson SA. Modulation of keratinocyte growth factor and its receptor in reepithelializing human skin. J Exp Med 182: 1369 –1376, 1995. 52. Marenholz I, Volz A, Ziegler A, Davies A, Ragoussis I, Korge BP, and Mischke D. Genetic analysis of the epidermal differentiation complex (EDC) on human chromosome 1q21: chromosomal orientation, new markers, and a 6-Mb YAC contig. Genomics 37: 295–302, 1996. 53. Marrakchi S, Audebert S, Bouadjar B, Has C, Lefevre C, Munro C, Cure S, Jobard F, Morlot S, Hohl D, Prud’homme JF, Zahaf A, Turki H, and Fischer J. Novel mutations in the gene encoding secreted lymphocyte antigen-6/urokinase-type plasminogen activator receptor-related protein-1 (SLURP-1) and description of five ancestral haplotypes in patients with Mal de Meleda. J Invest Dermatol 120: 351–355, 2003. 54. Maytin EV and Habener JF. Transcription factors C/EBP alpha, C/EBP beta, and CHOP (Gadd153) expressed during the differentiation program of keratinocytes in vitro and in vivo. J Invest Dermatol 110: 238 –246, 1998. 55. Mischke D, Korge BP, Marenholz I, Volz A, and Ziegler A. Genes encoding structural proteins of epidermal cornification and S100 calciumbinding proteins form a gene complex (“epidermal differentiation complex”) on human chromosome 1q21. J Invest Dermatol 106: 989 –992, 1996. 56. Morasso MI, Markova NG, and Sargent TD. Regulation of epidermal differentiation by a Distal-less homeodomain gene. J Cell Biol 135: 1879 –1887, 1996. 57. Morris RJ, Liu Y, Marles L, Yang Z, Trempus C, Li S, Lin JS, Sawicki JA, and Cotsarelis G. Capturing and profiling adult hair follicle stem cells. Nat Biotechnol 22: 411– 417, 2004. 58. Nagy N, Hoyaux D, Gielen I, Schafer BW, Pochet R, Heizmann CW, Kiss R, Salmon I, and Decaestecker C. The Ca2⫹-binding S100A2 protein is differentially expressed in epithelial tissue of glandular or squamous origin. Histol Histopathol 17: 123–130, 2002. 59. Nanney LB, Stoscheck CM, King LE Jr, Underwood RA, and Holbrook KA. Immunolocalization of epidermal growth factor receptors in normal developing human skin. J Invest Dermatol 94: 742–748, 1990. 60. Nelson WG and Sun TT. The 50- and 58-kdalton keratin classes as molecular markers for stratified squamous epithelia: cell culture studies. J Cell Biol 97: 244 –51, 1983. 61. Newman JC and Weiner AM. L2L: a simple tool for discovering the hidden significance in microarray expression data. Genome Biol 6: R81, 2005. 62. Nickoloff BJ, Griffiths CE, and Barker JN. The role of adhesion molecules, chemotactic factors, and cytokines in inflammatory and neoplastic skin disease–1990 update. J Invest Dermatol 94: 151S– 157S, 1990. 63. Nickoloff BJ, Qin JZ, Chaturvedi V, Denning MF, Bonish B, and Miele L. Jagged-1 mediated activation of notch signaling induces complete maturation of human keratinocytes through NF-kappaB and PPARgamma. Cell Death Differ 9: 842–55, 2002. 64. Nickoloff BJ and Turka LA. Keratinocytes: key immunocytes of the integument. Am J Pathol 143: 325–331, 1993. 65. Park GT and Morasso MI. Regulation of the Dlx3 homeobox gene upon differentiation of mouse keratinocytes. J Biol Chem 274: 26599 – 608, 1999. 66. Porter RM and Lane EB. Phenotypes, genotypes and their contribution to understanding keratin function. Trends Genet 19: 278 –285, 2003. 67. Radoja N, Stojadinovic O, Waseem A, Tomic-Canic M, Milisavljevic V, Teebor S, and Blumenberg M. Thyroid hormones and gamma interferon specifically increase K15 keratin gene transcription. Mol Cell Biol 24: 3168 –3179, 2004. 68. Rheinwald JG and Green H. Serial cultivation of strains of human epidermal keratinocytes: the formation of keratinizing colonies from single cells. Cell 6: 331–344, 1975. 69. Robinson NA, Lapic S, Welter JF, and Eckert RL. S100A11, S100A10, annexin I, desmosomal proteins, small proline-rich proteins, plasminogen activator inhibitor-2, and involucrin are components of the cornified envelope of cultured human epidermal keratinocytes. J Biol Chem 272: 12035–12046, 1997. www.physiolgenomics.org

78

PROFILING EPIDERMAL LAYERS

70. Rouillard JM, Herbert CJ, and Zuker M. OligoArray: genome-scale oligonucleotide design for microarrays. Bioinformatics 18: 486 – 487, 2002. 71. Segre JA, Bauer C, and Fuchs E. Klf4 is a transcription factor required for establishing the barrier function of the skin. Nat Genet 22: 356 –360, 1999. 72. Seitz CS, Lin Q, Deng H, and Khavari PA. Alterations in NF-kappaB function in transgenic epithelial tissue demonstrate a growth inhibitory role for NF-kappaB. Proc Natl Acad Sci USA 95: 2307–2312, 1998. 73. Sun L, Sun TT, and Lavker RM. CLED: a calcium-linked protein associated with early epithelial differentiation. Exp Cell Res 259: 96 –106, 2000. 74. Tseng H and Green H. Association of basonuclin with ability of keratinocytes to multiply and with absence of terminal differentiation. J Cell Biol 126: 495–506, 1994.

Physiol Genomics • VOL

27 •

75. Tsuruta D, Green KJ, Getsios S, and Jones JC. The barrier function of skin: how to keep a tight lid on water loss. Trends Cell Biol 12: 355–357, 2002. 76. Tumbar Guasch T, Greco G, Blanpain V, Lowry C, Rendl WE, Fuchs M, E. Defining the epithelial stem cell niche in skin. Science 303: 359 – 63, 2004. 77. Watt FM. Terminal differentiation of epidermal keratinocytes. Curr Opin Cell Biol 1: 1107–1115, 1989. 78. Woodcock-Mitchell J, Eichner R, Nelson WG, and Sun TT. Immunolocalization of keratin polypeptides in human epidermis using monoclonal antibodies. J Cell Biol 95: 580 –588, 1982. 79. Zambruno G, Marconi A, Melchiori A, Marchisio PC, Giannetti A, and De Luca M. Transforming growth factor (TGF) B1 enhances ␣5B1 and ␣vB5 integrin expression and migration of normal human keratinocytes. J Invest Dermatol 102: 593, 1994.

www.physiolgenomics.org