|Home | About | Journals | Submit | Contact Us | Français|
Continuous cell renewal in mouse epidermis is at the expense of a pool of pluripotent cells that lie in a well defined niche in the hair follicle known as the bulge. To identify mechanisms controlling hair follicle stem cell homeostasis, we developed a strategy to isolate adult bulge stem cells in mice and to define their transcriptional profile. We observed that a large number of transcripts are underexpressed in hair follicle stem cells when compared to non-stem cells. Importantly, the majority of these downregulated genes are involved in cell cycle. Using bioinformatics tools, we identified the E2F transcription factor family as a potential element involved in the regulation of these transcripts. To determine their functional role, we used engineered mice lacking Rb gene in epidermis, which showed increased expression of most E2F family members and increased E2F transcriptional activity. Experiments designed to analyze epidermal stem cell functionality (i.e.: hair regrowth and wound healing) imply a role of the Rb-E2F axis in the control of stem cell quiescence in epidermis.
The online version of this article (doi:10.1007/s12015-010-9139-0) contains supplementary material, which is available to authorized users.
Continually renewing tissues maintain their homeostasis at the expense of stem cells (SC) with the properties to self-renew and differentiate into the cell lineages that constitute their tissue of origin. Such characteristics are sustained by the niches where the SC reside, which provide specialized environments. SC functionality is critical for regenerating tissues damaged during injury . In addition, the relevance of SC in pathology is highlighted by multiple observations pointing to them as a possible origin of aggressive metastatic tumors . These particular properties, as well as their potential use in regenerative medicine, have led many groups to analyze and characterize different adult SC populations. In consequence, multiple approaches, including functional genomics, have been carried out to define specific biochemical and genetic determinants that define the actual SC attributes.
The control of stem cell properties in multiple tissues is linked to cell cycle, according to the cell quiescence demonstrated by most adult stem cells. The Rb family members (Rb1, Rbl1 and Rbl2) are negative regulators of cell cycle progression and have also been involved in the regulation of stem cell biology, including lineage specification of several types of stem cells . The Rb family functions in cell cycle control are mostly mediated by their ability to interact with E2F transcription factors, while their modulation of stem cell lineage specification may involve the regulation of tissue-specific transcription factors (reviewed in (3)). The E2F family of transcription factors comprises eight members (E2F1-8) and two related proteins (DP1 and DP2), essential for E2Fs binding to DNA (reviewed in (4)). Different functions have been assigned to E2Fs, E2F1 to −3 are considered activators and E2F4 and −5 repressors of target genes, in an Rb family (or pocket protein family)-dependent manner. E2F6 to −8 lack the Rb family binding domain, being considered pocket protein-independent repressors . Importantly, increased E2F activity due to inactivation of the pRb pathway has been found in most tumors and is associated to the increased proliferation in cancer cells .
The epidermis is the primary barrier that protects our bodies from the environment thus facing permanent wearing. To maintain its integrity this epithelium undergoes continuous cell renewal throughout the life of the individual at the expense of a small population of adult SC. In mouse epidermis, SC reside in a well established niche termed the bulge, a specialized portion of the hair follicle . Hair follicles periodically undergo cycles of growth (anagen), regression (catagen) and rest (telogen) . During the resting phase of the hair cycle (telogen), hair follicle SC (HFSC) are located at the bottom of the follicle, close to a specialized group of mesenchymal cells termed the dermal papilla (DP) and, at the anagen onset, DP cells activate HFSC forming a highly proliferative hair germ (HG). As the germ grows a proliferative compartment of transient amplifying cells (matrix cells) surround the DP at the base. The matrix cells terminally differentiate to produce the inner/outer root sheath and the hair shaft. When anagen ends, matrix cells apoptose and as a consequence the lower follicle degenerates and regresses (catagen), by this moment the DP is again located near the bulge and the sebaceous gland, resting until a new cycle starts. The source of cells that feeds these cycles comes from the bulge [9, 10]; however, bulge SC can also contribute to regenerate interfollicular epidermis after injury [9, 11] and sebaceous gland . Therefore, HFSC are multipotent progenitors capable of differentiating in all the cell lineages of the epidermis.
To determine which genes and signaling pathways operate within the HFSC, several groups have performed transcriptional profiling on isolated bulge cells [13–16]. Isolation procedures are mostly based on the specific expression of SC markers and/or the slow-cycling properties of these cells. Initial microarray experiments were designed to characterize cells that divide infrequently in the adult epidermis of transgenic mice expressing a repressible GFP marked histone ; however, hair follicle stem cells may not be as quiescent as initially thought [17, 18]. Other approaches took advantage of the restricted expression of cytokeratin 15 (K15) to mouse and human HFSC [19, 20]. Again, due to the intracellular nature of K15, it was necessary to develop transgenic mice that express GFP under the control of the K15 promoter to isolate viable cell populations with promoter activity . Direct isolation of HFSC circumventing the use of transgenic mice was possible using fluorescence-activated cells sorting (FACS) of CD34 positive cells. CD34 is a cell surface glycoprotein that has been shown to be distinctively expressed in the hair follicle bulge in murine epidermis [21, 22] where it may serve as an adhesion molecule, but its precise function is yet unknown. In mouse skin, CD34 expression is needed for tumor development , and CD34 positive populations isolated from epidermal tumors contained cells with tumor-initiating properties , supporting the notion that the HFSC population could represent an important target in the process of tumor development.
Gene expression studies in HFSC have led to the identification of genes and molecular pathways involved in SC maintenance, which in some cases have been validated through genetic approaches in transgenic mice [18, 25–27]. Similarly, the factors regulating HFSC gene expression are beginning to be elucidated and some transcription factors, identified in transcriptional profiling studies, have shown to be involved in the function and maintenance of this population, as it is the case of NFATc1 , Tcf3 [27, 28], Lhx2  and Sox9 .
We have used magnetic-activated cell sorting (MACS), as a novel way to isolate HFSC based on their CD34 cell surface expression , to develop a CD34+ HFSC signature that is robust and contains more genes than in previously reported datasets. This has allowed us to undertake the study of the regulation of gene expression in HFSC, in addition to the classical gene expression analyses. Our data suggest that the Rb/E2F families could be important regulators of HFSC gene expression. In agreement, functional data using mice deficient for Rb1 gene in epidermis, suggest that pRb ablation, which led to increased E2F activity, results in the dysregulation of stem cell quiescence.
Mice expressing GFP from the cytokeratin 15 (Krt1-15) gene promoter (Krt1-15-EGFP) were obtained from Jackson Laboratory (Jax Mice, Bar Harbour, MN). Mice lacking Rb1 gene in epidermis, through tissue specific ablation using K14Cre, have been previously described . All the animal experiments were approved by the Animal Ethical Committee (CEEA) and conducted in compliance with Centro de Investigaciones Energéticas, Medioambietales y Tecnológicas (CIEMAT) guidelines. For histological analysis, dorsal skin samples from wild type B6D2F1 mice were fixed in 70% ethanol and embedded in paraffin prior to sectioning. Sections were stained and processed as described . Antibodies used were anti-K15 (Neomarkers, Lab Vision, Fremont, CA), anti-CD34 (clone RAM 34, eBioscience, San Diego, CA), anti cytokeratin 5 (Covance, Emeryville, CA) and anti-BrdU (Roche, Indianapolis, IN). FITC, TexasRed or peroxidase conjugated secondary antibodies were purchased from Jackson ImmunoResearch (West Grove, PA). Control slides were obtained by replacing primary antibodies with PBS or preimmune sera (data not shown).
Keratinocyte isolation from mice back skin, immunofluorescent and magnetic labeling followed by MACS isolation of the stem population have been previously described . Briefly, adult mice in telogen (resting) phase of the hair cycle were used. Dorsal hair was removed and skin was placed on trypsin overnight at 4°C. Keratinocytes were collected scraping the epidermis with a scalpel and cell viability was determined by trypan blue exclusion. We confirmed that in the isolation procedure contamination with dermal cells is negligible . Cells were labeled (typically 1 μg antibody per 106 cells) with anti-mouse CD34 (eBioscience), anti-rat FITC (Jackson ImmunoResearch) and anti-FITC MicroBeads (Miltenyi Biotech, Auburn, CA). CD34-magnetic labeled (CD34+) and CD34 unlabelled fractions (CD34−) were separated using LS columns and LS column adapters (Miltenyi Biotec) following the manufacturer’s guidelines. Purity of the isolated populations was assessed by FACS and by direct visualization under the microscope (FITC fluorescence) and typically exceeded 90%.
For analysis of CD34 and cytokeratin 15 expression we used transgenic mice that express EGFP under the control of the Krt1-15 promoter. Keratinocyte populations isolated from Krt1-15-EGFP mice were labeled with anti-CD34 (eBioscience), washed in FACS Buffer (PBS containing 5% fetal bovine serum), stained with anti-rat PE (Jackson ImmunoResearch) and assessed for GFP and PE fluorescence. Cell preparations from wild type, RbF19/F19 and RbF19/F19;K14Cre mice were used to analyze CD34 and Itgα6 surface expression. CD34 was stained as described and anti-Itga6 antibody directly conjugated to FITC (BD Pharmingen, Franklin Lakes, NJ). Cells were analyzed in an EPICS XL flow cytometer (Coulter Electronics, Hialeah, FL). Dead cells and debris were excluded from the analysis.
Total RNA was prepared from CD34+ and CD34—keratinocytes obtained from pools of at least six mice per MACS isolation procedure . A total of three isolation procedures were performed. RNA extracted (RNAEasy, Qiagen) from each fraction (CD34+ / CD34−) was assessed for its quality using Agilent 2100 Bioanalyzer (Agilent Technologies) and hybridized to two Affimetrix Mouse 430 2.0 arrays, to ensure reproducibility of the process. This resulted in a total of 12 chips (6 CD34+ / 6 CD34−). Microarray data are accessible at NCBI Gene Expression Omnibus, http:www.ncbi.nlm.nih.gov/geo/, accession number GSE19448. Microarray normalization was done using RMA (http://rmaexpress.bmbolstad.com) . Differential expression analysis was done using Ttest (p-val<0.01, fold change>2).
Transcription factor consensus binding motifs search in genes specifically deregulated in CD34+ versus CD34−cells was done using Gene Set Enrichment Analysis software (GSEA, http://www.broadinstitute.org/gsea/) , using a collection of 615 motif gene sets, including i) gene sets that contain genes that share a cis-regulatory motif that is conserved across the human, mouse, rat, and dog genomes, and ii) gene sets that contain genes that share a transcription factor binding site defined in the TRANSFAC (version 7.4, http://www.gene-regulation.com/) database.
In order to find similar expression patterns between our MACS-purified CD34+ HFSC cells and published studies, we performed Gene Set Enrichment Analysis using our HFSC signature as the gene set, and the full microarray datasets of the other analyses. Log2 expression raw data of microarray analyses of skin stem cells of either mouse or human species was obtained online (see Table 1).
Ten-day old RbF19/F19 and RbF19/F19;K14Cre pups were injected with bromodeoxyuridine (BrdU) (20 μl of a 12,5 mg/ml dilution in NaCl 0,9%) every 12 h for a total of 4 injections. Skin sections were collected 75 days after the last injection and BrdU incorporation was measured as the percentage of hair follicles containing positive cells. Four different animals of RbF19/F19 and RbF19/F19;K14Cre were used to count at least 100 follicles.
Primary keratinocytes were obtained from newborn mice skin and cultured as previously described . Cells pooled from 5 to 8 different mice of each genotype were used for luciferase experiments or real time RT-PCR. E2F activity was measured in primary mouse keratinocytes transfected using Superfect (Qiagen) with pSV-Renilla (0,1 μg Promega) and pE2F-Luc (2.4 μg, a generous gift of Dr. X. Lu) as previously described . Luciferase activity, measured with the Dual Luciferase Reporter Assay system (Promega) was determined as the ratio of firefly to Renilla luciferase activity. Transfections were performed in triplicate, and the mean and standard error were calculated for each condition. Two independent transfection experiments were performed. Gene expression levels for each E2F family member was determined upon reverse transcription of 1 μg of total RNA using Transcriptor first strand cDNA synthesis kit (Roche, Basel, Switzerland). cDNA PCR amplification and analysis was done using Platinum® SYBR® Green qPCR SuperMix-UDG (Invitrogen, Carlsbad, CA) and Rotor Gene 3000 systems (Corbett Research, Cambridge UK). PCR primer sequences for E2F genes were previously reported . The expression values were first normalized using the 18 S ribosomal RNA (Universal 18 S rRNA kit, Ambion) and subsequently normalized to RbF19/F19 primary keratinocytes.
Two-month old RbF19/F19 and RbF19/F19;K14Cre mice (n=4 each) were anesthetized, shaven and 4 mm skin punch (4 per mouse) were made on the back. Wound healing was monitored by external caliper measurements every day. Data (mean±SE) are shown as the percentage of initial maximal wound diameter as a function of time after wounding. To monitor the histological process of wound healing, another cohort of RbF19/F19 and RbF19/F19;K14Cre mice (n=8 each) was used. In this case animals were sacrificed at different time points after wounding and the wound excised and fixed for analysis of the histology, BrdU incorporation, keratin 5 expression and collagen remodeling using sirius red staining. Mice were injected i.p. with 100 μg/g of body weight BrdU (Sigma Aldrich, St. Louis, MO) 1 h before they were sacrificed to analyze BrdU incorporation on the healing fronts.
To induce hair regrowth, two-month old RbF19/F19 and RbF19/F19;K14Cre mice (n=10 each) were shaven and depilatory cream was applied to the back skin. Animals were sacrificed at different time points after depilation and monitored for hair growth. Depilated skin was processed for histology and at least 5 sections (2 cm length) of each mouse were analyzed for the presence of anagen hairs.
To date, a variety of molecular markers and strategies have been employed to isolate hair follicle stem cells; particularly, the finding that the cell surface marker CD34 is enriched in hair follicle bulge stem cells  has eased direct studies on live cells avoiding the use of transgenic mice. Immunofluorescence and flow cytometry analysis showed partial overlap between CD34 and K15 expression in resting hair follicles (Fig. 1a–c), where K15 is present not only in bulge SC but also in cells of the hair germ, a small bulge-derived cell cluster that contributes to the formation of the new hair but shows less proliferative potential in vitro [16, 38]. Detailed FACS analysis of this CD34 population revealed the existence of two subsets of CD34+ keratinocytes with differential expression of the cell surface integrin alpha 6 (Itgα6) (Fig. 1d). Itgα6 mediates attachment to the basal lamina, thus these two CD34+ subpopulations have been dubbed basal (Itgα6 high) and suprabasal bulge cells (Itgα6 low) . Despite differences in Itgα6 surface expression both populations behave like true follicle stem cells: they display similar self-renewal capacity, proliferation potential and ability to terminally differentiate in vitro; and are fairly quiescent in vivo . At present, CD34 is the best characterized HFSC marker.
Based on the expression levels of the aforementioned markers, other groups have isolated follicle stem cells for transcriptional profiling using fluorescence-activated cell sorting (FACS) (Table 1). We decided to isolate CD34+ HFSC from non-genetically manipulated mice in the resting phase of the hair cycle using a magnetic-activated cell sorting (MACS) based method  to perform the gene expression analysis. This isolation procedure is efficient, fast and cell friendly, facilitating the recovery of high quality RNA. Purity of MACS isolated cells can be assessed by FACS analysis and is fairly similar to that obtained using FACS-sorting (Fig. 2). To our knowledge, this is the first time that MACS has been used to isolate CD34+ stem cells from the epidermis to perform genomic analyses.
We performed comparative gene expression analyses on CD34+ HFSC versus CD34−cells and identified 2,297 differentially expressed probes with statistical significance p<0.01 and a expression change (up or down) above 2-fold (Supplementary Table S1). Among them, there are many genes, such as Tcf3, Tcf4, Nfatc1, Lgr5 and Lhx2, that are well characterized as HFSC markers and are involved in the functionality of these cells (Fig. 3a). Interestingly, our HFSC signature contains more underexpressed (1,634 probe sets out of 2,297) than overexpressed genes, this feature is in agreement with previous studies [15, 39], and may indicate the status of a transcriptionally quiescent cell. Also, this HFSC signature contains more differentially expressed genes than previous datasets [13–16, 39], this is likely due to the combination of i) MACS cell isolation procedure that yields good-quality RNA, and ii) the use of the most comprehensive expression array available for our analysis.
Currently existing HFSC gene signatures have been developed using different species (human, mouse), selection markers, platforms, expression arrays and statistical software. To compare our HFSC signature with others we applied Gene Set Enrichment Analysis (GSEA). This analysis searches whether genes present in our 2,297-probe HFSC signature are also specifically regulated in bulge stem cells from previously reported datasets, including human (Table 1). The results from the GSEA showed that genes in our HFSC signature are significantly present in stem cells of all datasets (FDR q-val <0.0002), thus validating our experimental approach to isolate bulge stem cells based on the expression of the CD34 marker. Selecting for those genes enriched in HFSC from all datasets, we obtained a group of 156 genes (69 up / 87 down) characteristic of murine stem cells, termed moHFSC signature, and a subgroup of 59 genes (25 up / 34 down) that are also common to human hair follicle stem cells (Supplementary Table S2), termed hu-moHFSC gene signature. This gene list could represent bona fide molecular markers of skin stem cells in mice and humans. Interestingly, many of the genes common to all datasets are involved in cell cycle regulation.
We performed functional analyses of the HFSC-signature genes, and we searched for significant enrichment in Gene Ontology Biological Processes (GOBP) terms and KEGG pathways which are relevant for HFSC identity (Fig. 3b, Tables 2 and S3). We found that the most significantly enriched GOBP terms and pathways contain genes involved in cell cycle (Fig. S1). In particular, there is an extensive underexpression of genes implicated in cell cycle progression in the HFSC population, consistent with the essentially non-proliferative status of these cells. The identification of p53 and apoptosis as important pathways in HFSC is also a novel finding. Other pathways identified in our analysis, such as Wnt and Hedgehog/Bmp, have been previously reported to be regulated in HFSC [40–42], and our data add new genes regulated within these pathways which functionality may be also relevant.
Differential expression analyses have demonstrated a strong deregulation between CD34+ HFSC and CD34−mouse skin cells. We wanted to analyze possible transcription factors responsible for the expression differences between these two cell populations. For this purpose we searched for statistically significant enrichment in specific transcription factor sequence motifs within the promoters of genes overexpressed and underexpressed within the HFSC compartment using GSEA. This kind of analyses requires a critical number of genes to work with, and most of the previously reported datasets [13, 14, 39] were not large enough to undertake this study. Tables 3 and and44 show the analyses results, with the most significantly enriched transcription factor gene set motifs. Candidate regulators of the genes overexpressed in HFSC are homeobox, paired box and forkhead transcription factors (HFSC ON, Table 3); in fact it has been shown for some members of these families their role as regulators of stem cell homeostasis [26, 43, 44]. Using MatInspector , we studied the presence of binding sites in genes that are known HFSC markers and also regulate the biological properties of these cells, and found that the 5′ ends proximal to the transcriptional initiation sites of Lhx2, Lgr5 and Nfatc1 are loaded with potential binding sites for the above mentioned transcription factors (Fig. 4a). Similar observations were obtained using CD34, Sox9, Tcf3 and Tcf4 genes (not shown). We also found potential involvement of C/EBPs in regulating HFSC gene expression. We have recently reported that transcription factors of this family modulate epidermal terminal differentiation and loss of C/EBPα/β induces stem cell gene expression signatures in the epidermis, thus suggesting that they may act to restrict the epidermal stem cell compartment .
On the other hand, our results demonstrate a strong underexpression of E2F-dependent genes within the HFSC compartment, confirmed by the presence of various E2F-related gene set motifs (HFSC OFF, Table 4). More specifically, from 165 genes containing E2F binding sites within the expression array, 75 of them are specifically underexpressed in the CD34+ HFSC compartment (Fig. 4b), while the expression of most of the rest do not significantly distinguish between CD34+ and CD34−cells. We observed that, among the potentially E2F-regulated genes, E2F1, E2F3, and E2F7 are underexpressed in CD34+ HFSC. This suggests that the low expression of these E2F-dependent genes is not due to an active E2F-mediated repression, but rather to the relatively low expression level of these transcription factors.
The above reported findings led us to investigate the possible relevance of E2Fs in HFSC maintenance. Different E2Fs are present in epidermis where they exert specific roles in morphogenesis, differentiation and pathogenesis [46, 47]; however, the inherent complexity of this system, composed by at least ten different members with activator and repressor activities (see Introduction), makes it difficult to analyze their specific functions due to potential redundancy and functional compensation. In this regard, we have observed no evident phenotype in the skin of mice lacking E2F1 or E2F4 (CC, MS and JMP, unpublished data) in the absence of external stimuli, whereas overexpression of these proteins in epidermis leads to hyperplasic phenotype and different susceptibility to skin carcinogenesis protocols [48, 49].
Therefore, to analyze the relevance of E2Fs in HFSC homeostasis, we turned to a different approach using mice bearing epidermal specific ablation of the Rb gene using Cre-loxP technology. These mice display hyperplasia and hyperkeratosis in the epidermis with increased proliferation and aberrant expression of differentiation markers . We found that in Rb-deficient keratinocytes there is a significant overexpression of all E2F family members compared to wild-type cells, with the exception of E2F4, which appeared repressed (Fig. 5a). Such altered expression is accompanied by increased overall E2F activity as measured by luciferase reporter studies (Fig. 5b). Next, we studied how Rb deficiency and subsequent E2F activation affected HFSC homeostasis. The overall size of the HFSC population, measured as the percentage of cells showing CD34/Itgα6 positive phenotype in FACS analyses, is similar in Rb-deficient and control mice (Fig. 5c). In LRC experiments, Rb-deficient mice show a significant reduction in the number of hair follicles labeled 75 days after BrdU injection (Fig. 5d). BrdU-positive bulge cells appear faintly labeled in Rb-deficient mice (Fig. 5e´) as compared to the strong label displayed in control mice (Fig. 5e); also, we observed that, in Rb-deficient hair follicles, BrdU-labeled cells appeared usually in clusters rather than isolated cells (Fig. 5e and e´). Altogether, our data suggest an increased proliferation/division of the HFSC population in Rb-deficient mice leading to a dilution of the BrdU label. However, since Rb deficiency also induces disengagement between proliferation and differentiation in epidermal keratinocytes , we cannot rule out the possibility that a similar imbalance may occur during the differentiation of Rb-deficient HFSC, thus perturbing the functionality of these cells.
To explore the consequences of pRb ablation in HFSC behavior, we conducted functional studies. External stimuli, such as wounding or hair depilation, activate bulge stem cells to migrate upwards to contribute to epidermal regeneration, or downwards to regrow new hair. Time-course experiments analyzing hair regrowth showed accelerated anagen induction in Rb-deficient mice when compared to control littermates (Fig. 6a and andc),c), and only by day 8 after depilation a similar appearance was observed between Rb-deficient and control adult mice. This quickening in anagen induction resulted in evident faster hair regrowth in mutant mice even up to 12 days after depilation (Fig. 6b).
When the process of wound healing was analyzed, we observed a faster wound closure in Rb-deficient as compared to control mice (Fig. 7a and andb).b). Histological analyses demonstrate that this is accompanied by a faster re-epithelialization in mutant mice (Fig. 7c and andd).d). We also observed an increased cell proliferation at wound margins of mutant mice as compared to control (Fig. 7e). Finally, in the remodeling phase of wound healing, sirius red staining revealed a mild increase in the maturation of collagen fibers in Rb-deficient mice compared to control (Fig. 7f). These data are in agreement with the reported delayed in wound healing in E2F1 null mice, associated with substantially reduced rates of re-epithelialization without alterations in epidermal stratification and remodeling .
Overall, these data suggest that Rb loss results in increased stem cell proliferation, without an apparent loss of the HFSC pool or overt failure of stem cell function. Increased expression and transcriptional activity of E2Fs due to pRb depletion might indicate a functional role of this family of transcription factors in the maintenance of epidermal stem cell functionality; however, we cannot rule out the possibility that some of the observed effects can be driven by the interaction of pRb with other proteins . Ongoing studies, in which the epidermal specific ablation of pRb is combined with defined E2F deficiencies, will help to discern such possibilities.
The essential role of SC in tissue homeostasis and their possible involvement in multiple pathologies, have posed numerous efforts aimed to understand the specific and differential characteristics of this cell population. We have combined magnetic-activated cell sorting, as an isolation procedure, with functional genomic approaches to study the regulation of adult HFSC gene expression. Our findings point to a possible role of E2F transcription factors in the regulation of epidermal stem cell exit from quiescence.
In epidermis, several groups have reported gene expression profiling studies using FACS and/or sophisticated transgenic models [13–16, 39]. We have used MACS as an approachable and cell-friendly system to isolate CD34+ hair follicle bulge stem cells. In the resulting molecular signature we find upregulation of multiple genes that have been previously characterized as HFSC markers, such as Tcf3, Tcf4, Nfatc1, Lgr5 and Lhx2. We have validated this HFSC-signature through comparative analysis with existing ones. Using GSEA, we have found that a relevant number of genes identified by others are significantly present in our signature. An overall comparison of the gene expression profile revealed that ~70% of the differentially expressed genes are downregulated in the stem cell population, indicating that HFSC are less transcriptionally active than other epidermal cell populations; this feature is shared with other HFSC genomic analyses [15, 39]. Comparative analyses have allowed us also to identify genes with common expression patterns in all the published genomic datasets of HFSC, including mouse and human species. The mo-huHFSC signature represents a statistically robust list of genes over- or under-expressed in HFSC that could be useful for future analysis of HFSC markers, and highlights the importance of cell cycle regulation in this cell type.
Analysis of the molecular pathways in HFSC revealed that cell cycle control is the most relevant category enriched in this cell population, and may partly reflect the relative quiescence of these cells. We observed that genes implicated in cell cycle progression (cyclins and cyclin-dependent kinases) or mechanistically involved in DNA replication (Mcm genes, PCNA), as well as potential inhibitors of cell cycle progression, such as Cdkn1a, Rbl1 and Rbl2, are downregulated in HFSC. We have previously reported that, in the absence of p107 and p130 (Rbl1 and Rbl2 gene products, respectively), epidermal differentiation, including hair formation, is severely compromised  and Wnt signaling is altered . Furthermore, the epidermal ablation of Rbl1 dramatically aggravates the phenotype observed in Rb-deficient mice  and confers susceptibility to spontaneous  and chemically induced tumorigenesis . Altogether, these data suggest that the members of the Rb family could be playing a central role in epidermal stem cell activation.
Also, we have found downregulation of Trp53 and p53-dependent genes, including the aforementioned Cdkn1a. This is an apparent controversy with the widely accepted view that HFSC must possess robust mechanisms to detect possible genetic abnormalities and preclude appearance of pathologies such as tumors. However, since the CD34+ HFSC population is a fairly quiescent cell population, it cannot be excluded that transient activation of the p53-dependent pathway may occur during the process of HFSC activation to prevent expansion of genomic alterations, or during the process of hair regression as previously reported . In this regard, we have observed that mice lacking p53 in epidermis develop spontaneous tumors arising from hair follicles  associated with altered centrosome number and localization . Importantly, upon a chemical carcinogenesis challenge, the absence of Rb results in decreased tumor formation due to p53 activation  and the simultaneous absence of Rb and Trp53 results in accelerated tumor formation compared with mouse bearing only the epidermal ablation of Trp53 . Also, array data suggest that HFSC are relatively protected from apoptosis, as they show upregulation of anti-apoptotic or survival genes such as Bcl2 and Akt3, and downregulation of pro-apoptotic genes such as Trp53, Casp8, Dffa, Tradd or Ripk1. These observations may indicate that both Rb-and p53-dependent pathways could be involved in HFSC homeostasis during physiological or induced activation processes and merit future investigation.
The use of MACS and a large expression array has led to the development of a HFSC signature larger than previously reported ones [13, 14, 39]. This signature has allowed us perform a comprehensive bioinformatics analysis to characterize potential transcription factors regulating gene expression in HFSC. This kind of studies is significantly relevant only when the number of genes within a geneset or molecular signature is large enough. We observed that genes upregulated in HFSC contain potential binding sites for homeobox, paired-box and forkhead transcription factors. Some of them have been previously involved in the homeostasis of embryonic and adult stem cell populations, including follicle stem cells [26, 43, 44, 59–61], suggesting that a molecular circuitry common to multiple stem cells may exist. We also found potential involvement of C/EBPs in regulating HFSC gene expression. In this regard, we have recently reported that transcription factors of this family modulate epidermal terminal differentiation and they may act to restrict the epidermal stem cell compartment, since loss of C/EBPα/β induces stem cell gene expression signatures in the epidermis .
Our studies show the presence of E2F sites among a considerable number of genes underexpressed in HFSC. In this regard, it is important to note that the actual number of E2F-regulated genes is likely to be even higher than revealed in our study, taking in account the fact that E2F, in addition to its consensus binding sites, is capable of binding many gene targets that do not have typical E2F-consensus binding sites . This led us to investigate the possible involvement of E2Fs in the regulation of epidermal stem cell homeostasis; however, in vivo analysis is not straightforward due to the complexity of this transcription factor family, which comprises ten identified members with different activities (activation or repression) on gene transcription [4–6]. Therefore, we took advantage on the fact that the Rb ablation leads to increased expression of most E2F genes with a subsequent increase in E2F activity. We observed that the absence of pRb did not seem to alter the distribution of CD34+/Itgα6+ population in the skin. On the other hand, there was a significant reduction of the LRC population in the hair follicles of Rb-deficient mice, which frequently displayed cell clusters with decreased BrdU immunoreactivity. This label dilution points to an increased rate of HFSC proliferation and/or division in Rb-deficient mice. Another non-excluding explanation is that Rb gene ablation may affect the equilibrium between proliferation and differentiation in this cell population. In this regard, we have previously shown that epidermal keratinocytes lacking pRb display defective differentiation, and once the cells are committed into the differentiation program they are still able to proliferate [31, 37]. Future studies using CD34 purified cells from these mutant mice, or in vitro proliferation experiments, will help to clarify such possibilities.
In addition, our functional analysis did not reveal defective stem cell functionality or activation, but rather a more efficient wound healing and hair regrowth. Both are highly dynamic processes that require activation of epidermal stem cells and their interaction with dermal cells. Our data may suggest that Rb deficiency reduces the level of stem cell quiescence while maintaining the normal homeostatic functions and stemness of the epidermal SC population. This could be explained because these processes require that, upon HFSC activation, cell division proceeds; a mechanism that is facilitated by the absence of pRb and the subsequent increase in E2F activity . However, we cannot exclude the possibility that E2F-independent mechanisms triggered by Rb deficiency are partially responsible for these effects, in this sense it has been shown that pRb binds and regulates transcription factors, inhibitors of differentiation and chromatin modifiers to promote differentiation [64–67]. Nonetheless, the report on impaired wound healing in E2F1-null mice, due to defective cell proliferation and thus re-epithelialization , would reinforce the relevance of the Rb-E2F axis in maintaining appropriate stem cell proliferation in epidermis. In addition, our findings showing that Rb-deficient mice display increased E2F expression and activity and our expression data showing that many HFSC genes are E2F-targets, suggest that the consequences of pRb loss on stem cell quiescence are likely to be E2F-mediated.
In summary, we have applied a magnetic-based method to isolate CD34+ hair follicle stem cells to perform gene expression studies. We obtained a robust molecular signature that reveals the importance of the regulation of cell cycle related genes in this cell population. Analysis of the regulation of gene expression suggests that homebox, paired-box, forkhead and E2F transcription factors are important regulators of HFSC gene expression, and in vivo experiments using Rb-deficient mice point to a possible role of Rb/E2F in HFSC activation. The analysis of the different available E2F deficient mice, or the different combinations among them and with Rb family-deficient mice, will undoubtedly help to discern their possible specific functions in HFSC.
Grant support: Ministerio de Ciencia e Innovación (MICINN) grants SAF2008-0121, SAF2007-66227 and PS-010000-2008-7, Comunidad Autónoma de Madrid Oncocycle Program grant S2006/BIO-0232, and Ministerio de Sanidad y Consumo grant ISCIII-RETIC RD06/0020. JMA is a MICINN fellow, FPI-CIEMAT programme.
Conflict of interest The authors declare no potential conflicts of interest.
Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
Ramón García-Escudero and Carmen Segrelles contributed equally to this work.