p63 Regulates Human Keratinocyte Proliferation via MYC-regulated Gene Network and Differentiation Commitment through Cell Adhesion-related Gene Network

Background: Although the role of p63 in skin development is well known, its action in regulating human adult epidermal homeostasis is not completely elucidated. Results: p63 regulates keratinocyte proliferation via a MYC-regulated gene network and differentiation through a cell adhesion-related gene network. Conclusion: The balance between these two networks may control epidermal homeostasis. Significance: This is important to understand epithelial carcinogenesis. Although p63 and MYC are important in the control of epidermal homeostasis, the underlying molecular mechanisms governing keratinocyte proliferation or differentiation downstream of these two genes are not completely understood. By analyzing the transcriptional changes and phenotypic consequences of the loss of either p63 or MYC in human developmentally mature keratinocytes, we have characterized the networks acting downstream of these two genes to control epidermal homeostasis. We show that p63 is required to maintain growth and to commit to differentiation by two distinct mechanisms. Knockdown of p63 led to down-regulation of MYC via the Wnt/β-catenin and Notch signaling pathways and in turn reduced keratinocyte proliferation. We demonstrate that a p63-controlled keratinocyte cell fate network is essential to induce the onset of keratinocyte differentiation. This network contains several secreted proteins involved in cell migration/adhesion, including fibronectin 1 (FN1), interleukin-1β (IL1B), cysteine-rich protein 61 (CYR61), and jagged-1 (JAG1), that act downstream of p63 as key effectors to trigger differentiation. Our results characterized for the first time a connection between p63 and MYC and a cell adhesion-related network that controls differentiation. Furthermore, we show that the balance between the MYC-controlled cell cycle progression network and the p63-controlled cell adhesion-related network could dictate skin cell fate.

The outer layer of human skin, the epidermis, is a self-renewing, stratified squamous epithelial tissue. In adult tissue, the maintenance of epidermal homeostasis depends on an exquisite regulation of the balance between keratinocyte proliferation and differentiation (1,2). However, the specific molecular mechanisms governing each of these processes are not completely understood.
p63, a member of the p53 tumor suppressor gene family, is known as a key regulator of epidermis development and keratinocyte differentiation. Striking developmental defects have been discovered during embryonic development in p63-knockout mice (3,4). In addition, p63 is required for the maintenance of proliferative potential in epithelial stem cells (3,5,6), epithelial lineage commitment (4), differentiation of keratinocytes (7), and epithelial cell adhesion and survival (8). Dual roles of p63 in the initiation of epithelial stratification and maintenance of stem cell proliferative potential have been established (6). Because simultaneous p63 and p53 knockdown rescued the cell proliferation defect of p63 knockdown alone but failed to restore differentiation, the roles of p63 in proliferation and differentiation of developmentally mature keratinocytes appear to be distinct (7).
In vivo, as keratinocytes commit to differentiation, they detach from the basal layer and migrate outward into a spinous layer, accompanied by the expression of early differentiation markers keratin 1 (K1) and keratin 10 (K10) (9). The effects of p63 in keratinocyte differentiation have been investigated by gain and loss of function experiments. Ectopic expression of p63 induces the early marker K1 but suppresses the expression of late differentiation markers, such as loricrin and filagrin (10). Loss of p63 inhibited both stratification and differentiation in skin organotypic culture (7). However, little is known about key genes acting downstream of p63 to regulate the dynamic equilibrium of keratinocytes differentiation and proliferation as well as epidermal homeostasis (11).
Similar to p63, the MYC oncogene is predominantly expressed in the basal cell layers of epidermis and is absent in suprabasal layers (12,13). MYC plays a vital role not only in keratinocyte proliferation (14 -16) but also in accelerating the differentiation of epidermal stem cells (17)(18)(19)(20). Specifically, MYC activation leads to epidermal stem cell proliferation, and the sustained elevated expression of MYC in turn stimulates the proliferating stem cells to enter the transit-amplifying compartment, thereby initiating terminal differentiation (21). Although p63 and MYC have been independently demonstrated to be key regulators in the dynamic equilibrium between proliferation and differentiation in epidermal stem cells, no connection between p63 and MYC has been reported so far. The genetic basis underlying the processes regulated by p63 and MYC in mature keratinocytes, as well as the nature of their downstream effectors, remains mostly unknown.
We hypothesized that the regulation of the equilibrium between proliferation and differentiation of keratinocytes could rely on gene networks acting downstream of p63 and MYC rather than on a single gene. To this end, we compared phenotypic outcomes and global gene expression profiles in differentiating human keratinocytes transiently depleted of either p63 or MYC. We demonstrate that p63 regulates the proliferation and commitment to differentiation of human mature keratinocytes by two independent mechanisms. p63 is necessary for the proper expression of MYC and in turn controls keratinocyte proliferation by regulating the expression of MYC via the Wnt/␤-catenin and Notch signaling pathways. Moreover, p63 regulates human keratinocyte differentiation by controlling a network of genes implicated in cell migration and adhesion.

EXPERIMENTAL PROCEDURES
Cell Culture-Primary human keratinocytes were isolated from human skin biopsy and cultured in KGM2 medium (Clonetics) on flasks coated with collagen type I (Falcon Biocoat) at 37°C and 5% CO 2 . This study was approved by the ethical research committee "Comité de Protection des Personnes Sud-Est II" (CODECOH Number DC-2008-262).
Cell Proliferation and Cell Death Assays-12 h post-transfection, HaCaT cells were cultured in serum-free medium for 24 h. After a 40-h culture in complete medium, cell growth was performed with the ViaLight Plus kit (Cambrex) according to the manufacturer's protocol. For EdU staining, HaCaT cells were pulse-labeled for 24 h with 10 M EdU (Invitrogen) 24 h after transfection by siRNA. Cells were subsequently fixed on glass coverslips with paraformaldehyde immediately after the EdU pulse labeling. The Click-it EdU reaction was performed according to the manufacturer's instructions (Invitrogen) to label and detect EdU with the Alexa Fluor azide 555 dye. EdU staining was analyzed by BD LSRII (BD Biosciences). Student's t test was computed for statistical analysis. We used the fluorescence-based method LIVE/DEAD viability/cytotoxicity kit (Invitrogen) to evaluate cell death. HaCaT cells transfected by siRNA at 48 h after transfection were used for staining. Live cells (i.e. cells with esterase activity and intact membranes) were labeled with calcein AM, a green fluorescent product, and dead cells were labeled with ethidium homodimer-1, a red fluorescent nucleic acid stain. Cell staining was visualized with an Olympus BX61 straight microscope controlled with Metamorph software (Molecular Devices, Downington, PA).
Microarray Experiments and Genetic Network Analysis-For transcriptome analysis, microarrays with 25,342 spotted human oligonucleotides were used (Human-25K, Réseau National des Génopoles, France) (25). RNA was isolated, amplified, labeled, and hybridized following a published protocol (25). For each condition, three independent biological replicates were obtained, and the knockdown of corresponding genes was verified before RNA extraction (supplemental Fig.  S2). Each biological replicate was hybridized to the arrays in four replicates using a dye swap strategy. Slides were scanned with an Agilent G2565AA Microarray Scanner (Agilent Technologies). Data analyses, including intensity-dependent Lowess normalization of raw data and differential analysis, were performed with GeneSpring 7.0 software (Silicon Genetics). Differentially expressed genes were identified using analysis of variance (p value corrected with Benjamini and Hoehberg false discovery rate), and only genes with a -fold change of Ն1.2 and p Ͻ 0.01 were used for further analysis.
Gene pathway and genetic network analyses were performed by Ingenuity Pathway Analysis (IPA) (Ingenuity Systems). Ontologies attached to each gene were used to classify altered genes according to main biological themes. The following criteria enter into generating the networks in the Ingenuity Pathways Analysis application. 1) We have designated the molecules of interest in the analysis parameters before running the analysis. Molecules that (a) meet the cut-off and/or filter criteria and (b) interact with other molecules in Ingenuity's Knowledge Base are identified as focus molecules (also called "network eli-gible molecules"). Focus molecules serve as the "seeds," or focal points, for generating networks. 2) Networks are preferentially enriched for focus molecules with the most extensive interactions and for which interactions are specific with the other molecules in the network (rather than molecules that are promiscuous, those that interact with a broad selection of molecules throughout the Ingenuity knowledge base). 3) Additional nonfocus molecules from the data set and from the Ingenuity knowledge base are then recruited and added to the growing networks. 4) Networks are scored for the likelihood of finding the focus molecule(s) in that given network. The higher the score, the lower the probability that you would find the focus molecules(s) you see in a given network by random chance. 5) In the current version of the application, there is a cut-off of 35 molecules/network to keep networks to a usable size.
mRNA Expression Analysis-RNA was extracted with an RNeasy minikit (Qiagen). For real-time quantitative PCR, 2 g of RNA was reverse-transcribed in a total volume of 20 l with the SuperScript II RNase H reverse transcriptase system (Invitrogen) and random primers according the manufacturer's instructions. Reverse transcription reactions were diluted to 500 l of water, and 5 l of the diluted cDNA was used for each quantitative PCR. Quantitative PCR was carried out with a Platinum quantitative PCR SuperMix-UDG Kit (Invitrogen) using an ABI 7500 fast real-time PCR system (Applied Biosystems). All experiments were run in triplicate, and the results were normalized to 18 S rRNA expression. Primer sequences are listed in supplemental Table S1.
Cell Cycle Analysis-Cell cycle phase was determined by analyzing total DNA content using flow cytometry. HaCaT cells were synchronized by culturing in serum-free DMEM 24 h before transfection. Cells were transfected using INTERFERin with 10 nM siRNA. 48 h post-transfection, cells were collected, fixed, and stained with a cell cycle phase determination kit (Cayman Chemical). Flow cytometry analysis was performed in a MoFlo cell sorter (DakoCytomation).
Luciferase Reporter Assay-HaCaT cells were cotransfected with plasmid and siRNA in 24-well plates using an Amaxa Nucleofector II (Amaxa Biosystems) according to the manufacturer's recommendations. The human MYC promoter cloned into a luciferase plasmid was described previously (26). A TK-Renilla reporter was used as an internal normalization control. The ratio of firefly luciferase plasmid to Renilla luciferase was 10:1 in each nucleofection (2 g:200 ng). At 48 h post-transfection, luciferase activity was measured using a Dual-Luciferase reporter assay system (Promega). The luciferase reporter constructions were obtained as follows: pDel-1ϳpDel-4 reporter (Addgene plasmid 16601, 16602, 16603, 16604) (26), TCF-4/ LEF reporter plasmid from SuperArray Bioscience, YY-1 reporter plasmid from Panomics, and JUN and FOS reporter plasmid from H. van Dam (Leiden University, Netherlands).

siRNA-mediated Depletion of either p63 or MYC Decreased
Mature Keratinocyte Proliferation, but Only Loss of p63 Inhibited Differentiation-To investigate the specific roles of MYC and p63 in human keratinocyte differentiation, we knocked down each gene in HaCaT cells using specific siRNA and analyzed the phenotypic consequences. Because of the existence of six different isoforms of p63, we used an siRNA targeting the conserved DNA binding domain in all genes to achieve ablation of all p63 isoforms. We first verified that the expression of either MYC or p63 in siRNA-transfected keratinocytes was specifically knocked down at both the transcript (Fig. 1A) and protein levels ( Fig. 1B) 48 h post-transfection. Compared with cells treated with control siRNA, we observed a 70 and 90% downregulation of MYC and p63 expression, respectively. It is noteworthy that a significant loss of both genes was still observable after 10 days of culture (supplemental Fig. S1). Two days after siRNA transfection, we observed a reduced ATP content (Fig.  1C), whereas both EdU (Fig. 1D) and Ki67 staining (Fig. 1E) significantly decreased, both in p63-and MYC-knockdown HaCaT cells. Together, these data demonstrate a defect in keratinocyte proliferation upon ablation of either p63 or MYC. We next monitored the capacity of these cells to differentiate in vitro during 10 days of culture in the appropriate medium. As we focused on the commitment to differentiation, we chose to monitor the expression of two early differentiation markers, K1 and K10, rather than later ones, such as involucrin or filagrin. Indeed, K1 and K10 are markers of the basal-spinous layer transition in epidermis. Reduction of p63 levels in keratinocytes significantly inhibited the expression of K1 and K10 at the transcript level (Fig. 1F). In contrast, cells lacking MYC still expressed high levels of both markers, suggesting that these cells were still able to differentiate and did so even faster than the control (Fig. 1F). On the protein level, K1 was continuously induced except in cells lacking p63, which exhibited up-regulation of K1 until day 8, followed by a weak decrease in expression (Fig. 1G). Expression of the K10 protein started late in all conditions but became significant in MYC-depleted cells after 8 days of culture (Fig. 1G). In addition, we observed that loss of p63 in cultures of primary human keratinocytes inhibited expression of both K1 and K10 compared with control or MYCdepleted cells (Fig. 1H). These results demonstrate that the ablation of either MYC or p63 significantly reduced human mature keratinocyte proliferation, whereas only the siRNAmediated loss of p63 inhibited differentiation.

Knockdown of either MYC or p63 Leads to Cell Cycle Arrest-
To determine the potential origin of the proliferation defect in keratinocytes lacking either p63 or MYC, we analyzed cell death and the cell cycle. Knockdown of either p63 or MYC triggered a significant arrest in G 0 /G 1 ( Fig. 2A) without affecting cell death (Fig. 2B). To identify the molecular mechanisms controlling this cell cycle arrest, we analyzed the gene expression profiles of keratinocytes depleted of either MYC or p63 after transfection (supplemental Table S2). For these experiments, we have particularly insisted on robustness of the data. We performed three independent siRNA transfections and three independent RNA extractions (supplemental Fig. S2), along with several technical replicates, to generate three independent expression profiles that we then averaged for either MYC-or p63-depleted keratinocytes. This transcriptome analysis showed that a small net-work of genes strongly associated with cell cycle regulation (p Ͻ 6.2 ϫ 10 Ϫ35 ) was significantly down-regulated in keratinocytes lacking either MYC or p63 (supplemental Fig. S3A and Fig. 2C). The majority of genes in this network were similarly regulated in both p63-depleted and MYC-depleted cells, except CSK2, GADD45A, and CCND2 (supplemental Fig. S3A). Analysis by qRT-PCR of the major cell cycle inhibitors confirmed that p15 and p21 were both significantly up-regulated in response to the knockdown of either MYC or p63 in HaCaT cells (Fig. 2, D and E). Interestingly, the same trend was observed in primary keratinocytes (Fig. 2, F and G). We also observed that MYC, a major hub (highly connected node) in this cell cycle-controlling network, was downregulated in both HaCaT and primary keratinocyte cells lacking p63 (Fig. 2, C, D, and F). This result suggests that p63 is necessary for the proper expression of MYC in human keratinocytes.  Fig. S4) (27). Using specific siRNA targeting either ⌬Np63 or TAp63, we were able to demonstrate that only the knockdown of ⌬Np63 isoforms triggered downregulation of MYC expression (Fig. 3C). This was consistent with the fact that ⌬Np63 is the isoform predominantly expressed in adult keratinocytes. Conversely, overexpression of either ⌬Np63 or TAp63 in HaCaT cells had no effect on MYC expression (supplemental Fig. S5). Together, these results demonstrated that p63 is necessary for the proper expression of MYC.
We next investigated the molecular mechanisms underlying the p63 knockdown-triggered down-regulation of MYC expression. Using luciferase reporter constructs fused to truncated MYC promoter regions (Fig. 3D), we characterized the sequence upstream of MYC and identified a putative p63-controlled region. As demonstrated in Fig. 3E, this region extended from Ϫ349 to Ϫ607 bp upstream of the transcription start site.  FEBRUARY 17, 2012 • VOLUME 287 • NUMBER 8 A close-up view of this region revealed the absence of the p63 consensus binding site, in agreement with our previous published results, which showed that MYC was not a direct target of p63 in human keratinocytes (28,29). These results suggested that the expression of MYC is under the indirect control of p63. The region upstream of the MYC gene also contains several other binding sites for transcription factors (TFs), including TCF-4 (TF responding to the Wnt/␤-catenin signaling pathway), YY-1 (TF responding to the Notch signaling pathway), c-FOS, and c-JUN (AP1 TFs) (Fig. 3F). Because these TFs are known to play important roles in the control of keratinocyte proliferation and differentiation (30 -32), we investigated whether or not they are involved in the p63 knockdown-triggered down-regulation of MYC promoter activity. In keratinocytes lacking p63, we observed a moderate inhibition of both TCF4 or YY1-dependent luciferase activities but no change or a slight activation in JUN-and FOS-driven luciferase expression (Fig. 3G). These results suggest that TCF4 and YY1, the TFs responding to the Wnt/␤-catenin and Notch signaling path-ways, respectively, were partially responsible for the p63-dependent down-regulation of MYC expression.

p63 and Keratinocyte Cell Fate
p63

Knockdown-triggered Down-regulation of MYC Expression Is Mediated by Notch and Wnt/␤-Catenin Signaling
Pathways-To determine whether the signaling cascades leading to the activation of YY1 and TCF4 are affected by p63 levels, we scrutinized the expression profiles of keratinocytes lacking either p63 or MYC. Interestingly, we found that the Wnt/␤catenin and Notch signaling cascades were indeed potently inhibited in keratinocytes lacking p63. JAG1 and DLL1, ligands of the Notch signaling pathway, were both down-regulated in p63-knockdown cells (Fig. 4A). In addition, several inhibitors of the Wnt/␤-catenin pathway were up-regulated, whereas several activators were down-regulated, probably resulting in the inhibition of this signaling pathway (Fig. 4A). At the protein level, the ablation of p63 resulted in the down-regulation of both cleaved Notch1 (NICD) and total ␤-catenin (Fig. 4B). Taken together, these results suggest that both pathways were turned down in p63-knockdown cells. Inhibition of the Notch

p63 and Keratinocyte Cell Fate
signaling pathway in the presence of p63-targeting siRNA was also indicated by the reduced expression of two reporter genes of the Notch-signaling cascade, HES1 and HEY1 (Fig. 4C). To further confirm the potential involvement of the Wnt and Notch pathways in down-regulation of MYC expression, we restored downstream signaling by transfecting HaCaT cells with either WNT3 or NICD expression vectors, together with p63-targeting siRNA or control siRNA. We observed that the overexpression of either WNT3 or NCID induced the up-regulation of MYC expression both at the transcript (Fig. 4D) and protein levels (Fig. 4E). This would partially restore keratinocyte proliferation. However, overexpression of those two genes was not sufficient to restore proper differentiation of keratinocytes lacking p63 (data not shown). These results suggest that the molecular role of p63 in differentiation of human mature keratinocytes is probably different from the ones that dictate their proliferation.
Cell Migration/Adhesion-related Network Acts Downstream of p63 to Induce Onset of Keratinocyte Differentiation-Our results demonstrated that MYC is down-regulated in cells lacking p63 (Figs. 2 and 3), thus functionally corresponding, at least partially, to a siRNA-mediated knockdown of MYC. However, these cells exhibited completely opposite differentiation outcomes (Fig. 1E). To investigate the molecular mechanisms enabling keratinocyte differentiation downstream of p63, we compared the expression profiles of p63-depleted and MYCdepleted cells.
As demonstrated in Fig. 5A, 546 genes were common to both expression profiles. It is noteworthy that there were more genes common to both profiles than specific to the p63-depleted cells. This again suggests that part of the transcriptional response to p63 ablation in human keratinocytes was also due to the downregulation of MYC. The down-regulated genes in keratinocytes lacking either MYC (supplemental Fig. S6A) or p63 (supplemental Fig. S6B) shared similar gene ontology terms (e.g. cell cycle, DNA replication, and DNA repair). Genes up-regulated in either MYC-depleted (supplemental Fig. S6C) or p63-depleted keratinocytes (supplemental Fig. S6D) also shared some gene ontology terms, such as cellular movement or cell death.
Among the 546 genes common to both expression profiles, we found 71 genes that were antagonistically regulated (supplemental Table S3). We hypothesized that these antagonistically regulated genes could mechanistically explain, at least partially, the opposite differentiation outcomes between p63-and MYClacking keratinocytes. We used the Ingenuity knowledge base using IPA software to analyze the networks and functions associated with these 71 genes. Strikingly, a network of 41 nodes was extracted and significantly associated with a single function, cell migration/adhesion (p Ͻ 3 ϫ 10 Ϫ14 ). In cells lacking p63, this network was strongly down-regulated (Fig. 5B), whereas in MYC-depleted keratinocytes, this same network was up-regulated (supplemental Fig. S3B). We further validated the expression of some genes in the network by qRT-PCR in HaCaT cells (Fig. 5C) and in primary keratinocytes (Fig. 5D). As expected, the expression of these genes was up-regulated in MYC-depleted cells and down-regulated in keratinocytes lacking p63. These results suggest that this migration/adhesion-related gene network could contain potential effectors acting down- To validate our hypothesis, we used several functional approaches. First, we searched for known phenotypes associated with this 41-gene network in the Mouse Genome Informatics database. There are 19 knock-out mice with abnormal skin phenotypes reported in that database. Strikingly, 15 genes of these 19 KO mice were present in this network (Table 1). These data suggest that nearly 80% of all known skin dysfunction-related genes belong to the network we have characterized and functionally validate it. We also monitored in vitro differentiation of keratinocytes lacking different genes belonging to this migration/adhesion network: FN1, MMP13, JAG1, IL1B, and CYR61 (Fig. 6A). Except for MMP13, we observed a strong inhibition of differentiation in cells lacking any of these genes, as demonstrated by the delayed expression of K1 (Fig. 6B) and K10 (Fig. 6C) transcripts. It is noteworthy that although the ablation of IL1B was only partial (40%), it significantly inhibited differentiation. Furthermore, ablation of all of these genes together (siCocktail) strongly inhibited keratinocyte differenti-ation. Finally, if the genes belonging to this network promote commitment to differentiation, we postulated that their expression should be down-regulated in non-differentiated and/or pluripotent cells. We data-mined the NCBI Gene Expression Omnibus database, and interestingly, we found that seven hubs in this network, PLAU, FN1, IL1B, ADM, DUSP10, GADD45A, and RAC2, were significantly down-regulated in induced pluripotent stem cells (iPS) and are even part of the iPS transcriptomic signature (Fig. 6D) (33). Together, these results show that a p63-controlled migration/adhesion-related network plays a key role in the onset of human mature keratinocyte differentiation. As a consequence, we named this network the keratinocyte cell fate (KCF) network.

DISCUSSION
Our findings confirm that sustained expression of both p63 and MYC, two major regulators of epidermal homeostasis (5-7, 18, 21, 34, 35), is required to maintain growth and differentiation of human developmentally mature keratinocytes. However, we demonstrate that their respective roles are very differ- p63 and Keratinocyte Cell Fate ent, as already suggested by some authors (36). We propose a model to illustrate the distinct mechanisms of action of p63 on human developmentally mature keratinocyte proliferation and differentiation (Fig. 7). p63 is required for the proper expression of MYC through the combined regulation of the Wnt/␤catenin and Notch signaling pathways, leading in turn to cell cycle regulation and cell proliferation. p63 also regulates a KCF network that contains several potential "differentiation effectors," some located in the extracellular space. The up-regulation of the KCF network would promote the onset of terminal differentiation of keratinocytes (Fig. 7). In this study, we show that the siRNA-mediated loss of MYC triggers down-regulation of the "proliferation network" and up-regulation of the KCF network, promoting human keratinocyte differentiation; in

List of all knockout mice exhibiting abnormal skin phenotypes reported in the Mouse Genome Informatics database
Genes that are also present in the KCF network (Fig. 5B) and are antagonistically regulated in MYC-and p63-deficient keratinocytes are shown in boldface type. JNK and PI3K families figure in the network, and therefore no specific accession number could be given.   contrast, p63 knockdown down-regulates both cell proliferation and mobility/adhesion-related networks, thus inhibiting differentiation. These results were observed both in the HaCaT cell line (mutated p53) and in normal human primary keratinocytes (wild-type p53), suggesting that the model we propose is independent of the p53 status of skin cells. Because MYC is a transcriptional repressor of p15 and p21, two cyclin-dependent kinase inhibitors (37)(38)(39), and p63 was reported to repress p21 and p16 (40 -42), it was not a surprise to observe reduced proliferation in cells lacking either gene. More surprising was the p63 knockdown-triggered down-regulation of MYC expression that we report here for the first time. This result suggests that although MYC is not a direct target of p63 (28,29), the expression of this transcription factor is necessary for the proper expression of MYC both in HaCaT cell line and primary human keratinocyte culture. Our results establish that MYC expression is mediated, at least in part, by the Wnt/␤-catenin and Notch signaling pathways. Interestingly, both pathways are important regulators of proliferation and differentiation in epidermal stem cell maintenance and wound healing (1,(43)(44)(45)(46). The interplay between the Wnt/␤catenin and Notch pathways has been reported in epidermal homeostasis and differentiation as well (47,48). We show that both pathways act in concert downstream of p63 to control the proper expression of MYC and, in turn, regulate the keratinocyte cell cycle.

Gene ID GenBank TM accession number Phenotype
Although the inhibition of proliferation of cells lacking p63 was partially MYC-dependent, the differentiation defects appeared to be independent of MYC. Indeed, the siRNAmediated loss of MYC did not impair keratinocyte commitment to terminal differentiation. The differentiation was even slightly accelerated in MYC-depleted cells. By comparing expression profiles from mature keratinocytes lacking either MYC or p63, we found 546 common genes, of which 71 were antagonistically regulated. Strikingly, from that list of 71 genes, we extracted a gene network containing 41 genes that were significantly associated (p Ͻ3 ϫ 10 Ϫ14 ) with a single function, cell migration/adhesion. These results are consistent with recent reports showing that p63 functions as an inhibitor of cell migration (27) and that a p53 mutant forms a complex with p63 to antagonize its cell migration-inhibitory function, leading to TGF␤-dependent metastasis (49). Similarly, it was shown that p63 regulates a cell adhesion program, including integrins, in epithelial cells (8). Although we cannot exclude the possibility that other genes that are among the 546 common genes but are not found in our network might regulate early differentiation, we have clearly demonstrated that this p63-controlled cell migration/adhesion network contains several effectors acting downstream of p63 to trigger differentiation of mature keratinocytes.
Some of the differentiation effectors we have identified (such as integrins, FN1, PLAU, JAG1, IL1, and CYR61) are also involved in cancer progression in inducible human tissue neoplasia (50). Integrin signaling plays an important role in epidermal adhesion, growth, and differentiation (51,52). JAG1 is a ligand of the Notch signaling pathway and acts on keratinocyte differentiation (32,53). JAG1 is also a transcriptional target of p63 (54). Interleukin-1 (IL1) is implicated in human epidermal keratinocyte proliferation (55) and even in the regeneration of epidermal tissue in vitro (56). IL1␣ is the active form of interleukin-1 in human epidermis. IL1␤ was considered non-functional in keratinocytes, but our data suggest that IL1␤ is necessary to induce keratinocyte differentiation.
Although most of these "differentiation effectors" do not seem to be direct transcriptional targets of p63, the control exerted by p63 on this network was dominant. Indeed, we were unable to rescue differentiation of p63-depleted keratinocytes by ectopic expression of NICD or JAG1 with the use of JAG1 or IL1␤ recombinant proteins or even with an acellular matrix obtained from normal fully differentiated keratinocytes (data not shown).
We report for the first time the role of this p63-regulated cell migration/adhesion network in the commitment of developmentally mature keratinocytes to differentiation. A normal expression of this network seems to be required to trigger differentiation, whereas its down-regulation prevents it. Furthermore, we believe that misregulation of this gene network may play a major role in tumorigenesis. Indeed, Khavari's group recently reported a core tumor progression signature network in keratinocytes, which contained 282 nodes and was involved in carcinogenesis (50). This core tumor progression signature network contained several oncogene hubs, and eight of the top 10 nodes were extracellular or cell surface proteins. It is noteworthy that four of these eight extracellular oncogene hubs, PLAU, CYR61, FN1, and IL1, also belong to the p63-regulated cell migration/adhesion network we describe in this study. Other oncogene hubs reported in the core tumor progression signature network (50), such as SERPINE1 and ITGA6, were also down-regulated in human keratinocytes depleted of p63 (supplemental Table S2).
In conclusion, the siRNA-mediated loss of MYC triggers down-regulation of the "proliferation network" and up-regulation of the "KCF migration/adhesion-related network," promoting human keratinocyte differentiation; in contrast, p63 knockdown down-regulates both cell proliferation and KCF FIGURE 7. Schematic representation of the mechanism of commitment of keratinocyte differentiation controlled by p63. p63 controls keratinocyte proliferation and differentiation independently via different genetic networks. p63 regulates the cell cycle in part by regulating MYC expression through both the Wnt/␤-catenin and Notch pathways. To trigger the onset of differentiation, p63 controls the KCF network composed of keratinocyte differentiation effectors. Most of these effectors, such as FN1, IL1B, JAG1, and CYR61, are located in the extracellular matrix and are implicated in cell migration.

p63 and Keratinocyte Cell Fate
networks, thus inhibiting differentiation. We believe that the balance between levels of expression of both cell proliferation and KCF networks could dictate keratinocyte cell fate. Furthermore, we think that this network approach would reconcile much of the existing data on the regulation of the balance between proliferation and differentiation in skin.