Distinct Pigmentary and Melanocortin 1 Receptor–Dependent Components of Cutaneous Defense against Ultraviolet Radiation

Genetic variation at the melanocortin 1 receptor (MC1R) is an important risk factor for developing ultraviolet (UV) radiation–induced skin cancer, the most common form of cancer in humans. The underlying mechanisms by which the MC1R defends against UV-induced skin cancer are not known. We used neonatal mouse skin (which, like human skin, contains a mixture of melanocytes and keratinocytes) to study how pigment cells and Mc1r genotype affect the genome-level response to UV radiation. Animals without viable melanocytes (KitW-v/KitW-v) or animals lacking a functional Mc1r (Mc1re/Mc1re) were exposed to sunburn-level doses of UVB radiation, and the patterns of large-scale gene expression in the basal epidermis were compared to each other and to nonmutant animals. Our analysis revealed discrete Kit- and Mc1r-dependent UVB transcriptional responses in the basal epidermis. The Kit-dependent UVB response was characterized largely by an enrichment of oxidative and endoplasmic reticulum stress genes, highlighting a distinctive role for pigmented melanocytes in mediating antioxidant defenses against genotoxic stresses within the basal epidermal environment. By contrast, the Mc1r-dependent UVB response contained an abundance of genes associated with regulating the cell cycle and oncogenesis. To test the clinical relevance of these observations, we analyzed publicly available data sets for primary melanoma and melanoma metastases and found that the set of genes specific for the Mc1r-dependent UVB response was able to differentiate between different clinical subtypes. Our analysis also revealed that the classes of genes induced by UVB differ from those repressed by UVB with regard to their biological functions, their overall number, and their size. The findings described here offer new insights into the transcriptional nature of the UV response in the skin and provide a molecular framework for the underlying mechanisms by which melanocytes and the Mc1r independently mediate and afford protection against UV radiation.


Introduction
One of the most important functions of cutaneous pigmentation in humans is protection against the damaging effects of ultraviolet (UV) radiation. Depending on wavelength and intensity, UV radiation can have broad-ranging effects on DNA damage, cell cycle arrest, and apoptosis [1] in virtually every one of the more than 25 differentiated cell types in the skin, as well as more specialized responses such as immunosuppression, vitamin D synthesis, and sunburn/tanning [2]. Inadequate protection against UV radiation is a major contributor to melanoma and nonmelanoma skin cancer and a significant public health concern in many populations [3].
Pigmentary defenses against UV radiation depend on both quantitative variation in the number, size, and arrangement of melanosomes-pigment granules transferred from melanocytes to surrounding keratinocytes-and qualitative variation in the type of pigment, eumelanin or pheomelanin, made within those granules. In general, increasing skin darkness correlates closely with an increased number and size and greater dispersal of melanosomes [4,5] and with increasing amounts of both eumelanin and pheomelanin (although eumelanin predominates) [6,7]; this quantitative variation is controlled by polygenic inheritance. By contrast, a single major gene, melanocortin 1 receptor (MC1R), can regulate the ratio of eumelanin to pheomelanin. Complete loss-offunction for MC1R, which encodes a seven-transmembrane receptor coupled to adenylate cyclase, produces a so-called red hair color phenotype in individuals of northern European ancestry-bright red hair, fair skin, freckling, inability to tan, and increased susceptibility to sunburn-and is associated with nearly exclusive production of pheomelanin [8][9][10]. Pedigree studies are generally consistent with Mendelian expectations for recessive inheritance of the red hair color phenotype [11,12], but in many populations there is considerable MC1R diversity with a range of hypomorphic alleles [13][14][15][16], and several association studies suggest a more complex relationship between MC1R genotype and pigmentary phenotype that depends on both gene dosage [12,17,18] and biochemical activity [19][20][21] of specific alleles.
In addition to effects on pigmentation, MC1R loss of function is associated with an increased incidence of melanoma and nonmelanoma skin cancer in Australian [22,23], Mediterranean [24][25][26], and northern European [27,28] populations. Although fair skin is a well-known risk factor for skin cancer, most of the aforementioned studies conclude that increased skin cancer susceptibility due to MC1R variation cannot be explained solely by pigmentary phenotype. Potential explanations for this remarkable finding include differential effects of melanin subtypes on cancer susceptibility that have not or cannot be captured by skin spectrophotometry or categorical classification of skin color and/or a protective effect of MC1R signaling that occurs independently of pigmentation. Supporting the former idea are in vitro and/or cell culture studies demonstrating that pheomelanin is phototoxic [29,30], generating reactive oxygen species [31][32][33] and promoting cell death [34,35] in response to UV radiation. However, differences in melanin subtype seem unlikely to explain all of the protective effects against skin cancer afforded by normal MC1R genotype, since both pheomelanin and eumelanin are positively correlated with skin color darkness and the content of both increases following UV radiation [6]. Indeed, several investigators have suggested that UV and MC1R signaling have synergistic actions on both pigmentation and susceptibility to skin cancer [27,[36][37][38].
Much of what we know about how UV radiation affects skin physiology is based on cell culture studies in which particular behavioral aspects, like the extent of survival, melanogenesis, proliferation, dendricity, or DNA repair, were measured in a specific skin cell type, like melanocytes, upon exposure to UV radiation [39][40][41]. More recently, large-scale gene expression studies have been carried out on a variety of primary or established cell lines exposed to different doses and wave-lengths of UV. These approaches have begun to show how the response to DNA damage elicits changes in cell biology that ultimately translate into whole animal phenotypes. To better understand how the skin responds to UV radiation as an integrated tissue, we have examined the whole genome transcriptional response to UVB radiation in neonatal mouse skin using a combination of genomic, genetic, and bioinformatic approaches. Although adult mouse skin lacks epidermal melanocytes in hairy areas (and therefore differs substantially from the histologic architecture of human skin), neonatal mouse skin has an abundance of active epidermal melanocytes [42,43] and a histologic architecture that shares many features with human skin [44]. Most important, a mouse model provides genetic tools to directly assess the role of melanocytes and MC1R genotype. Notably, our results identify a specific signature for Mc1r-dependent UVB transcriptional responses, characterized by a predominance of genes implicated in the cell cycle and oncogenesis. The Mc1rdependent UVB transcriptional response is clinically relevant, since the genes identified in our study can partition two human melanoma data sets into clinical subtypes. This work provides new insight into the cutaneous response to UV radiation and reveals mechanisms whereby interactions between UV radiation and MC1R signaling may predispose toward the pathogenesis of skin cancer.

A UVB-Responsive Transcriptional Profile in the Basal Epidermis
Between postnatal day 1.5 (P1.5) and P3.5 when hair follicles are still developing, melanocytes are relatively abundant in the basal epidermis, which provides a useful model to study melanocyte-keratinocyte interactions. In control experiments to examine how the transcriptional response to UV radiation changes over time, we determined that 24 h captures both early and intermediate response genes and provides substantial information in terms of the underlying gene expression profiles (unpublished data). For the experiments described below, P1.5 mice were exposed to a single dose of 100 mJ/cm 2 UVB radiation, which corresponds to mild sunburn, approximately one minimal erythemal dose; control animals underwent the identical procedure but were shielded from UV radiation with a black Perspex box. At 24 h after exposure, we isolated basal epidermis and compared the patterns of gene expression among different experimental groups.
Our preparation of basal epidermis is based on a combination of enzymatic and physical dissociation steps that separates whole epidermal sheets, basal epidermal cells, and whole skins; in previous studies we have shown that this approach yields skin layer-specific gene expression signatures that correspond to the known biological functions of these compartments and in which the basal epidermal compartment is composed of mostly keratinocytes and melanocytes. All experiments were carried out as triplicate biological replicates (each replicate consisted of a separate pool containing two to four animals) with a common reference as described previously.
To establish a baseline, we first determined how UVB affects gene expression in C57BL/6J animals. From spotted cDNA arrays containing 31,756 features, we obtained a data

Author Summary
Skin cancer is the most common type of cancer in humans and annually accounts for approximately 60,000 deaths worldwide. The most important factors causally linked to skin cancer susceptibility are inadequate protection against ultraviolet (UV) B radiation, fair skin color, and variation of the melanocortin 1 receptor (MC1R) gene. We used cDNA microarrays to measure the genome-wide transcriptional responses to UVB irradiation in the epidermis of neonatal mice (which approximates the human basal epidermis in its cellular composition and general physiology). To investigate how pigment cells (melanocytes) and MC1R afford protection against UVB radiation, we compared results from normal mice to those from mutant mice that lacked either melanocytes (Kit W-v /Kit W-v ) or a functional Mc1r (Mc1r e /Mc1r e ). We identified melanocyte-and Mc1rdependent UVB gene expression profiles in the basal epidermis. Surprisingly, the melanocyte-and Mc1r-dependent UVB responses highlighted distinct functions, with the former largely mediating antioxidant defenses and the latter regulating the cell cycle and susceptibility to oncogenesis. We also demonstrated that a subset of Mc1r-dependent UVB-responsive genes could discriminate among human melanoma subtypes, thereby suggesting a mechanism by which MC1R gene variants may predispose toward skin cancer.
set of 13,508 cDNAs that passed quality control filters for signal intensity and replication. At a false discovery rate (FDR) of 0.05 [45], 1,434 of the 13,508 cDNAs were differentially regulated by UVB ( Figure 1A and Table S1). The distribution of quantitative changes in expression was similar for induction compared to repression ( Figure 1B); however, there was a striking difference in the number of cDNAs affected, with many more repressed than induced by UVB: 1,065 versus 369. Large-scale transcriptional repression has also been observed in other microarray studies of UV radiation. Indeed, the cDNA whose UVB-induced repression exhibited the most significant p-value in our data set was Polr2e, one of a dozen subunits of RNA polymerase II.
Annotation using UniGene (http://www.ncbi.nlm.nih.gov/ entrez/query.fcgi?db¼unigene) Cluster IDs (assigned for the majority of the UVB-responsive features on the array) revealed that 770 of 1,065 of the UVB-repressed (UVBr) and 332 of 369 of the UVB-induced (UVBi) cDNAs were unique. Readily apparent are several representatives of gene classes that are well known to participate in the response to UV radiation. Examples include those involved in cell cycle arrest (Cdk2 and Cyclin E1), DNA repair (Ercc1, Fanc1, and Rad51), and antiapoptosis (Aatf and Dad1), all of which were induced by UVB. The Tp53 pathway is known to play an important role in mediating many of these stress responses, and indeed, we found several Tp53-dependent genes in our UVB- The ordered columns represent replicate samples of sham-(black) and UVB-(blue) irradiated mice, and the rows represent log 2 ratios for 1,434 UVB-responsive genes. Data are displayed as a pseudocolor heatmap where red, green, and black represent induction, repression, and no change, respectively, relative to the median. Gray represents missing data. The ratio fold-change is indicated by a scale bar. Genes are ranked by their average fold change (D) value. The average fold D value was obtained by dividing the average expression level of each gene in the experimental (UVB-irradiated) group by that in the control (sham-irradiated) group. Average expression levels for the experimental and control groups are represented as log 2 ratios 6SEM. The top ten ranked genes for the repressed and induced responses are listed along with their corresponding average fold change (D) and Student's t-test p-values. In some cases, official mouse gene symbols are used, instead of the lengthier gene names. (B) Distribution of fold-change differences for UVB-responsive genes. Fold-change differences were obtained by comparing sham-and UVB-irradiated mice. Data were binned into one of five fold-change classes for each of the repressed and induced responses. The total number of genes sorted into each bin is shown above each class. doi:10.1371/journal.pgen.0030009.g001 responsive gene list, including Mcm2, Hif1a, and Bcl3, a finding that, in general, is further supported by a recent microarray study reporting a set of p53-dependent UVB-responsive genes in primary human melanocyte cultures [46].
In a recent study of UV-sensitive and p53-dependent gene expression profiles in a colorectal carcinoma cell line, McKay et al. [47] noted that genes induced by UV tended to be more compact and to have smaller introns and proposed that a gene size constraint on UV-induced mRNA expression played a key role in the evolution of UV-response pathways. Remarkably, we found a similar bias in our system, with the mean gene size for UVBi genes (approximately 30 kb) approximately 1.8 times smaller than that of the UVBr genes (approximately 54 kb) and a mean intron size for the UVBi genes (approximately 4,000 bp) approximately 1.5 times smaller than that of the UVBr genes (approximately 5,900 bp) ( Figure 2). The mean number of introns per gene was also significantly smaller in the UVBi compared to the UVBr gene class (unpublished data). There were no differences in gene size, intron size, and the number of introns per gene between the set of UVBr genes and the entire set of approximately 35,000 genes on the array.

Functional Annotation of UV-Responsive Genes by Over-Representation Analysis
To better characterize the UVB transcriptional skin response in normal C57BL/6J mice, we extracted Gene Ontology (GO) terms from all UVBr or UVBi genes and compared the relative frequencies of these GO terms to those for all of the genes on the 35,000 array. Using EASE software [48], we observed 220 and 199 GO categories that were significantly enriched (p ¼ 0.05) for the UVBr and UVBi genes, respectively. Of these, 119 UVBr and 64 UVBi categories were supported by at least three gene hits and were enriched by at least 2-fold over the 35,000 array. We focused on nonoverlapping biological process, molecular function, or cellular component terms for further study; results are presented in Figure 3.
For the 1,065 UVBr genes ( Figure 3A), over-represented GO terms were largely consistent with those which we had observed in a previous study using basal epidermis [49]. Biological processes over-represented by UVBr genes in the basal epidermis included many metabolic and signaling functions: sterol and fatty acid metabolism, ATP metabolism, purine ribonucleoside triphosphate biosynthesis, regulation of Wnt receptor signaling pathway, antigen processing, endogenous antigen via MHC class I, and ribosome biogenesis. Molecular functions over-represented by UVBr genes in the basal epidermis largely corresponded to the aforementioned processes: MHC class I receptor activity, acetylcoenzyme A C-acetyltransferase activity, clathrin binding, protein disulfide isomerase activity, selenium binding, and iron binding.
The pattern of GO terms over-represented for the 369 UVBi genes in the basal epidermis was very different from that of the UVBr genes ( Figure 3B versus Figure 3A). Biological processes over-represented by UVBi genes included DNA repair (nucleotide-excision repair) and response to DNA damage stimulus, S-phase and G2/M transition of mitotic cell cycle, DNA replication initiation, rRNA metabolism, and nuclear mRNA splicing. Molecular functions overrepresented by UVBi genes included nuclease and helicase activities, damaged DNA binding, pre-mRNA splicing factor activity, DNA-directed DNA polymerase activity, and cyclindependent protein kinase activity.
As described further below, UVB-responsive genes identified in our study and others exhibit only a partial overlap when considered from the perspective of individual genes, in part due to differences in array platforms and in part due to different patterns of gene expression between skin in vivo and cells in culture. However, large-scale transcriptional repression, downregulation of basal metabolic functions, upregulation of DNA repair functions, and bias in gene and Figure 2. UVB Response and Gene Size Gene (A) and intron (B) size data were retrieved for 305 of 332 UVBinduced, 705 of 770 UVB-repressed, and 10,316 of 14,605 35,000 unique mouse genes (by UniGene Cluster ID criteria) using the Ensembl (National Center for Biotechnology Information m34 Assembly) and National Institute on Aging Gene Index 5.0 databases. The size data were averaged for each UVB response class and were then compared, using a Student's t-test, to the gene and intron size averages for 10,316 unique genes on the 35,000 array. Significance values are shown above each UVB response class, and error bars are 6SEM. Gene and intron sizes are represented here on a linear scale, but were log 2 -transformed prior to ttest analyses, due to the non-normal distribution of the data. Nonparametric Mann-Whitney U tests on the untransformed data yielded similar results. doi:10.1371/journal.pgen.0030009.g002 intron size between these functional classes are general themes emerging from our analysis of the UVB response in vivo and are concordant with earlier work carried out on a variety of cultured cells and UV exposures.

Genetic Deconvolution of Gene Expression in the Basal Epidermis: Melanocytes (Kit-Dependent) and Mc1r Signaling
To investigate how melanocytes and Mc1r signaling contribute to epidermal gene expression, we made use of two classic pigmentation mutants, viable dominant spotting (Kit W-v /Kit W-v ) and recessive yellow (Mc1r e /Mc1r e ). The Kit W-v allele is a missense mutation that severely impairs tyrosine kinase activity [50], yielding animals that, post embryonic day 12.5, are almost completely bereft of skin and hair melanoblasts [51]; the Mc1r e allele is a complete loss-of-function that causes almost exclusive production of pheomelanin instead of eumelanin [52]. Using Mc1r e and Kit W-v alleles on a congenic C57BL/6 background, we measured gene expression profiles in basal epidermis of P2.5 mutant animals and compared the results to those obtained from nonmutant C57BL/6 animals. Kit signaling is also required for development of germ cells, erythrocytes, and mast cells. However, pigment cell development is, in general, more sensitive to Kit gene dosage than mast cell development (Kit W-v /Kit W-v mutant mice retain approximately 40% of the normal number of mast cells), and mast cells are normally found in the dermis rather than the epidermis; therefore, we expected that epidermal gene expression in Kit W-v /Kit W-v mutant mice would be influenced mostly by deficiency of pigment cells.
We first characterized patterns of gene expression in Kit W-v / Kit W-v and Mc1r e /Mc1r e animals in the absence of UV radiation. We identified 430 genes whose expression in the basal epidermis was altered in Kit W-v /Kit W-v animals: 140 downregulated and 290 upregulated ( Figure 4A). Among the downregulated group, over-represented GO terms included melanin biosynthesis (48.7-fold enriched; p ¼ 6.51EÀ04), protein and electron transport, adenyl and guanyl nucleotide binding, organelle organization, and biogenesis, as well as small GTPase-mediated signal transduction (Table S2). Among the upregulated group, over-represented GO terms included embryogenesis and morphogenesis, development, extracellular space, cell adhesion molecule, and signal transducer activity (Table S3). While some of the downregulated genes simply reflect the loss of melanocytes (and melanocytespecific genes), most genes whose expression was altered in Kit W-v /Kit W-v animals are not melanocyte-specific. These results indicate that the gene expression profile in Kit mutant mice represents a transcriptional response rather than a simple change in cellular composition and suggest that, in  [48] was used to map GO terms to (A) 1,065 UVB-repressed (green bars) and (B) 369 UVB-induced (red bars) genes. GO terms that were significantly enriched (p ¼ 0.05), by 2-fold, relative to all of the assayed genes on the 35,000 array (35,328 genes) (black bars) and that contained at least three gene hits, are displayed. GO terms are ranked, within high-level GO branches, by fold-enrichment relative to the 35,000 array. A thumbnail image provides a reference for each of the analyzed UVB response classes. doi:10.1371/journal.pgen.0030009.g003 addition to pigmentation, melanocytes influence the development and maintenance of the extracellular environment in the basal epidermis.
We also identified 284 genes whose expression was altered in Mc1r e /Mc1r e compared to nonmutant animals (in the absence of UV radiation): 184 downregulated and 100 upregulated ( Figure 4B). Approximately 10% of the genes downregulated in Mc1r e /Mc1r e animals were also downregulated in Kit W-v /Kit W-v animals ( Figure 4C), including two genes required for eumelanogenesis, Dopachrome tautomerase and Silver. However, there was virtually no overlap between Mc1rand Kit-dependent genes when conditioning on UVB exposure (described further below); thus, Mc1r signaling is an important component in how the epidermis senses and responds to pigment cells in the basal state, but Mc1r signaling plays a specific and largely independent role from Kit (and the presence of melanocytes) in the transcriptional response to UV radiation.
We compared the profiles of Mc1r-and Kit-dependent genes in the basal epidermis with the 1,102 unique genes that were either induced or repressed by UVB in nonmutant skin as depicted in Figure 1  To investigate if melanocytes and Mc1r signaling are required specifically for the cutaneous response to UV radiation, we used the same UVB exposure paradigm described earlier and compared gene expression profiles in the basal epidermis of P2.5 Kit W-v /Kit W-v and Mc1r e /Mc1r e to that obtained from nonmutant animals. The total data set consisted of 18 arrays with three biological replicates per sample (three genotypes, either exposed or not exposed to UV radiation).
We used two-way analysis of variance (ANOVA) to characterize the interaction between genotype and UVB exposure and identified 120 and 147 genes that showed a significant interaction (p ¼ 0.05) between UVB and Mc1r and between UVB and Kit, respectively. Post hoc analyses performed separately on each of the 120 and 147 gene sets yielded three patterns, depending on whether a gene responded to UVB only in the mutant, only in the nonmutant, or in both mutant and nonmutant genotypes ( Figure  5). Collectively, we refer to the 120 UVB 3 Mc1r genes as the Mc1r-dependent UVB-responsive genes (Table S4) and the 147 UVB 3 Kit genes as the Kit-dependent UVB-responsive genes (Table S5).
Surprisingly, we found very little overlap between the Mc1r-and the Kit-dependent UVB-responsive genes (only five in common of a total of 267), suggesting that the biological processes engaged by Mc1r signaling in the response to UVB are distinct from the ability of pigment cells to protect against UV radiation by providing a simple physical barrier. To further explore this idea, we extracted and compared GO terms derived from the Mc1r-and the Kitdependent UVB-responsive gene sets using the same approach described earlier (Figure 3). Significant (p ¼ 0.05) GO terms that contained at least three hits and that were enriched by at least 2-fold over the 35,000 array are depicted in Figure 6 for both sets of genes. The Mc1r-dependent UVB response was characterized largely by genes involved in the cell cycle (Cdc2l1, Fos, Hspa8, Mtcp1, Top3a, and 2810406C15Rik) and oncogenesis (Cdc2l1, Fgfr1, and Mtcp1) and included functions and components relating to GTPasemediated signal transduction (Grb2, Rab7, Rab12, Rab35, and Ralbp1), chromatin (H2faz, Scmh1, and 5430405G24Rik), ribosomes (Rpl8, Rpl13, and Rps16), and mitochondria (Atpif1, C1qbp, Cndp2, Mtcp1, Pdha1, and Phb).
To better understand how the set of Mc1r-dependent UVBresponsive genes in mice might be related to MC1R genotype as a risk factor for human skin cancer, we considered Mc1rdependent UVB-responsive genes known to have GO-anno- We also considered whether any of the Mc1r-dependent UVB-responsive genes had previously been implicated in cancer biology by virtue of their aberrant regulation. Of 20 genes for which roles in cancer susceptibility or progression had previously been suggested (Table 2), 75% exhibited the direction of change expected for a protective Mc1r-dependent UVB response in the basal epidermis. For example, whereas Arpc1b is normally induced in the Mc1r þ /Mc1r þ basal epidermis, it remained unresponsive in Mc1r e /Mc1r e animals and is frequently silenced in gastric cancers [53]. In contrast, while Mapkapk2 was normally unresponsive to UVB in the basal epidermis of Mc1r þ /Mc1r þ animals, it was induced in Mc1r e /Mc1r e animals and has been reported to be activated in some breast cancers [54].

MC1R-Dependent UVB-Responsive Genes, Melanoma Profiling, and Cancer Biology
Given the over-representation of cell cycle and oncogenesis genes in our Mc1r-dependent UVB response, we sought to determine whether the expression patterns of these genes might be informative in a clinical setting. Among a number of previous microarray studies of human melanoma or melanoma cell lines, we identified two large data sets based on tissue samples from either primary melanoma [55] or melanoma metastases [56].
Starting with our 120 Mc1r-dependent UVB-responsive genes, we were able to retrieve gene expression data for 86 and 53 orthologous human genes from the primary melanoma and the melanoma metastases data sets, respectively (see Materials and Methods for further details), of which approximately one third (47) were shared between the two data sets. Hierarchical clustering revealed that our Mc1r-dependent  UVB-responsive genes were able to partition the primary melanoma ( Figure 7A) and the melanoma metastases ( Figure  7B) samples into two subtypes each. Both sets of melanoma subclusters obtained here were nearly identical to those originally identified, with the primary melanoma subtypes [55] corresponding to TNM stage (I-II versus III-IV), and the melanoma metastases subtypes [56] corresponding to growth phase (type I radial, in situ and minimally invasive tumors, versus type II vertical, aggressive metastatic tumors). Of the 86 genes tested on the primary melanoma data set, expression levels for 19 distinguished the different stages, with all of the 19 underexpressed in the stage III-IV samples relative to the stage I-II samples ( Figure 7A). Of the 53 genes tested on the melanoma metastases data set, expression levels for 33 distinguished the different growth phases, with 10 and 23 genes overexpressed and underexpressed in the type I radial growth relative to the type II vertical growth samples, respectively ( Figure 7B). Closer inspection of the overlap between the 19 and 33 genes used to distinguish the primary melanoma and melanoma metastases data sets revealed nine genes in common: ACTN1, ARPC1B, IMP3, KIAA0553, MGC3731, MTCP1, RAC1, RALBP1, and SDFR1.

Validation, Comparison, and Context
The microarray platform and design used here are identical to those of a previous study in which we compared gene expression profiles in different skin layers of nonmutant and Mc1r e /Mc1r e mice and carried out selective validation by semiquantitative reverse-transcription polymerase chain reaction [49]; nearly all of the genes that we previously identified as Mc1r-dependent in the basal epidermis were also detected in the current set of experiments. To examine these findings in a broader context, we compared our UVBresponsive gene list to those reported in other large-scale UV expression experiments. From 28 published studies using a variety of different experimental designs and platforms (species, cell type, tissues, wavelength, dose, and time) (Table  S6), we retrieved the UV-responsive gene lists and then combined and sorted these lists to yield a nonredundant set of 5,178 UniGene IDs. From the cDNAs shown to be UVB responsive in the basal epidermis ( Figure 1 and Table S1), we identified a nonredundant set of 1,201 UniGene IDs. There are 502 genes in common between the two sets (5,178 from 28 earlier studies and the 1,201 reported here); this degree of overlap is highly significant (p , 0.0001), yet indicates that our approach yielded a large number of genes (699) ( Table  S7) not previously known to be UV responsive.

Discussion
Of 28 previous studies designed to study the large-scale transcriptional response to UV radiation, only three examined tissues or tissue compartments; the remainder focused on cultured cells. In two studies from Hochberg and colleagues [57,58] where the UV response of human suprabasal epidermis (obtained from suction blisters) was compared directly to that of cultured human keratinocytes using the same platform, only 207 of 1,931 genes found to be UV responsive in either the suprabasal epidermis or cultured keratinocytes were UV responsive in both samples.
Our findings provide additional support for the notion that substantial differences exist between the identities of UVresponsive genes in tissues and cultured cells, yet reinforce several common themes. First, transcriptional repression accounts for a large proportion (75%) of the UV response, a finding consistent with that seen, although not always reported, in other large-scale UV gene expression studies. Second, downregulation of basal metabolic genes together with upregulation of DNA damage and/or repair genes characterizes the UV transcriptional response across a broad set of cell types and several tissues. Finally, the quantitative characteristics of the UV response are characterized by small changes in expression that affect a large number of genes. All of these features may be linked to the observation made here and by McKay et al. [47] that UV-induced and UV-repressed genes differ not only in their biological function but also in their size, suggesting evolutionary pressure to ensure that the very genes required for a critical response to UV radiation are themselves less likely to sustain UV-induced DNA lesions.
In addition to providing a broader perspective on the transcriptional response to UV radiation, our work helps to explain the epidemiologic relationship between MC1R genotype and skin cancer susceptibility. The set of Mc1rdependent UVB-responsive genes we identified is greatly enriched for processes and functions involved in oncogenesis and the cell cycle; moreover, the pattern in which specific genes change suggests that a functional Mc1r contributes to the ability of UVB to induce cell cycle arrest in the basal epidermis ( Figure 8). Our work on the genome-level response to UV radiation was necessarily carried out in laboratory mice, where the availability of specific mutations on a defined genetic background allows a substantial degree of experi- mental control. However, our conclusions about the general relationship between MC1R genotype and cancer susceptibility in the mouse are supported by previous studies of two specific melanoma genes in humans: CDKN2A (a tumor suppressor) and BRAF (a proto-oncogene). In families carrying a germline CDKN2A mutation, the risk of melanoma was substantially increased by the presence of an MC1R variant [59][60][61][62][63]. In recent work from Landi et al. [64], germline MC1R variation was found to be a risk factor for developing melanomas with BRAF mutations. Because the association between MC1R variation and BRAF was most apparent in melanomas without evidence of chronic sun damage, Landi et al. hypothesized that UV radiation acted indirectly to promote BRAF-mutant melanoma, which is also consistent with our findings.
The nature of our experimental system-an integrated tissue collected from whole animals-makes it difficult to assess whether Mc1r signaling might affect UVB-induced cell cycle arrest in keratinocytes, melanocytes, or both. However, in recent work from Yamaguchi et al. [65] where the effects of UV radiation on cutaneous DNA damage and cell death were examined in individuals with different pigmentary phenotypes, light skin was found to exhibit a surprising lack of UVinduced apoptosis in many cells of the basal epidermis. Although MC1R genotypes were not determined in the subjects studied by Yamaguchi et al. [65], our work suggests that some of the effects they observed might be attributed to MC1R heterozygosity.
What sort of molecular mechanisms might link MC1R signaling to UV-induced apoptosis [39][40][41]? From a cell- Orthologous human gene expression data, derived from 120 mouse genes that were identified as Mc1r-dependent and UVB-responsive in the basal epidermis of P2.5 mice, were clustered in (A) a primary melanoma data set [55], here consisting of 86 genes and 45 samples, and (B) a melanoma metastases data set [56], here consisting of 53 genes and 19 samples. Melanoma subclusters correspond to the original subtypes, TNM stage I-II (blue) versus III-IV (orange) in (A), first identified by Talantov et al. [55], and growth phase type I radial (orange) versus type II vertical (blue) in (B), first identified by Haqq et al. [56]. Genes driving the observed melanoma clustering patterns are displayed to the right of each large cluster, along with their corresponding UVB response (", induction; #, repression; and Both, regulation in both autonomous perspective, phosphorylation of p53 Ser389 is a critical and major step in triggering apoptosis in response to UV-induced DNA damage [66], and it is possible that MC1R activation (either constitutive or in response to UV radiation) plays a permissive role. However, recent studies by Abdel-Malek and colleagues [67] indicate that MC1R deficiency causes increased UV-induced apoptosis in cultured melanocytes. Thus, a more likely scenario is that decreased MC1R signaling in melanocytes exposed to UV radiation in vivo impairs the ability of surrounding cells to undergo p53mediated cell cycle arrest and apoptosis. Indeed, previous work from our group and others indicates that while Mc1r expression is confined mainly to melanocytes, Mc1r deficiency leads to substantial changes in gene expression throughout the skin and suggests that melanocytes influence the behavior of surrounding cells via paracrine mechanisms. For example, L-type prostaglandin D2 synthase (Ptgds) is expressed at high levels in pigment cells [68] and is repressed in Mc1r e /Mc1r e skin (current data, [49]); prostaglandin D2 and its derivatives could be one of several mediators released from Mc1r-competent melanocytes that facilitate or potentiate p53 phosphorylation in surrounding cells [69]. Regardless of the molecular mechanisms that link Mc1r signaling to changes in the expression of cell cycle genes, our analysis of human melanoma samples indicates that those genes are relevant to melanoma biology, since our set of Mc1r-dependent UVBresponsive genes discriminates between the severity of primary melanoma subtypes and between the growth phases of melanoma metastasis subtypes.
Several known antioxidant genes, including Glo1, Prdx5, Sepp1, and Stk25, were all repressed in basal epidermis of Kit þ / Kit þ animals but not Kit W-v /Kit W-v animals. These results are consistent with reports in which UVB, in inducing oxidative stress, was reported to deplete the antioxidant defense system in the epidermis of mice [70,71] and humans [72]. We also identified additional genes, within the Kit-dependent UVB response, that were not annotated in the GO database as oxidative or endoplasmic reticulum (ER) stress response genes, including Cdk5, Eif1, Egln1, Ext1, Gpr175, Gpt2, Itm1, Nars, Pabpc1, Pecam1, Prdx4, Prkrip1, Pten, Rheb, Sh3glb1, Tbc1d20, and Tuba3, for which roles in oxidative stress/hypoxia and the ER stress/unfolded protein response have been reported [73][74][75][76]. Several other genes involved in the ER/Golgi trafficking or folding of proteins, but not previously reported to be regulated during oxidative stress/hypoxia or the unfolded protein response, were also identified, including Ext1, Ergic3, Grb2, Inpp5e, Man1a2, Prkcsh, Rer1, and Surf4, raising the possibility that these genes may represent new skin stress response genes. Given the absence of this oxidative/ER stress response signature in our Mc1r-dependent UVB response, we interpret our findings to suggest that melanocytes, despite constituting a small fraction of the skin cell population, play a major role during oxidative and ER stress in the basal epidermis and, significantly, that this role is independent of Mc1r signaling.
Taken together, our results provide a molecular framework for the underlying mechanisms by which melanocytes and the Mc1r independently mediate, and afford protection against, UVB signals within the normal skin environment. This latter finding is of particular significance as it suggests insights into the molecular mechanisms by which MC1R mutations predispose toward melanoma susceptibility.

Materials and Methods
Mice and UVB irradiation. The administered UVB dose corresponded to a mild sunburn dose, approximately one minimal erythemal dose, for the C57BL/6J strain [77]. Control pups from all three genotypes were shielded from UV irradiation with a black Perspex box. Immediately postirradiation, all pups were returned to their mothers for 24 h, after which they were killed. Because Kit W-v /Kit W-v mice are sterile, homozygous neonates were obtained by heterozygous matings followed by a polymerase chain reaction and NciI restriction digest-based diagnostic genotyping assay for the Kit W-v mutation [78]. Dorsal trunk skins were dissected and basal epidermal cells were prepared and stored as previously described [49].
Experimental design. The methodology described here conforms  RNA extraction and labeling. Total RNAs were extracted from the basal epidermal cells using TRIzol Reagent (Invitrogen) followed by further purification using a RNeasy cleanup protocol (Qiagen, http:// www.qiagen.com). Experimental and reference RNA samples (20 lg each) were reverse-transcribed and then directly labeled using an anchored oligo-dT 18 primer (Qiagen) and a CyScribe First-Strand cDNA Labeling Kit (Amersham Biosciences, http://www.amershambiosciences.com), as per the manufacturer's instructions. The RNA was degraded in 1N NaOH for 30 min at 70 8C followed by neutralization in 1N HCl. The Cy5-dUTP-labeled basal epidermal and Cy3-dUTP-labeled reference samples were combined, concentrated using Centricon YM-30 microconcentrators (Millipore/Amicon, http://www.millipore.com), and then competitively hybridized to a 35,000 spotted glass mouse microarray slide.
Data processing and statistical analyses. Images were acquired with a GenePix 4000A scanner and feature measurements were extracted using GenePix Pro 3.0 software (Axon Instruments, http://www.axon. com). Raw data were uploaded into the Stanford Microarray Database (publicly available at http://smd.stanford.edu) where backgroundsubtracted fluorescence ratios were linear-normalized such that the mean Cy5/Cy3 ratio on each array was 1. Good-quality spots were selected by manually excluding spots that were irregularly shaped or of small diameter or contained excessive background intensities. We next selected spots whose median signal intensity was 1.5-fold above the median background for both Cy3 and Cy5 channels. Our normalized, background-subtracted, raw data set totaled 31,756 spots. For our statistical analyses, the raw data set was filtered such that measurements were present in at least two of three replicates per experimental condition. This filter reduced our data set to 13,508 cDNAs (9,110 unique genes). Student's tand two-way ANOVA tests were performed using GeneSpring Software (Agilent, Palo Alto, California, United States) on centered data sets, in which the median log 2 ratio for each gene equaled 0. Over-representation analysis was performed using EASE [48], which reports functional gene categories as statistically significant GO terms. Significance values were adjusted for multiple hypothesis testing by the FDR, and values with an FDR of 0.05 were selected for further analyses.
Unsupervised hierarchical clustering of melanoma data sets. Using our list of 120 Mc1r-dependent UVB-responsive mouse genes, together with National Center for Biotechnology Information's HomoloGene and BLASTn algorithms, we retrieved gene expression data for 147 and 74 orthologous human genes from the primary melanoma and melanoma metastases data sets of Talantov et al. [55] and Haqq et al. [56], respectively. The log 2 ratios for both data sets were median-centered separately, and nonredundant human gene lists for each of the data sets were generated by averaging expression data for gene replicates. The resulting primary melanoma data set, now consisting of 86 genes and 45 samples, and the melanoma metastases data set, now consisting of 53 genes and 19 samples, were clustered separately using an unsupervised algorithm in Cluster [79], in which uncentered correlation coefficients (similarity metric) for both samples and genes underwent average linkage clustering. The results are displayed as heatmaps using Treeview software [79]. Table S1. UVB-Responsive Genes in the þ/þ Basal Epidermis of P2.5 Mice Found at doi:10.1371/journal.pgen.0030009.st001 (684 KB PDF).