Stem cells expand potency and alter tissue fitness by accumulating diverse epigenetic memories.

HFSCs can partake in long-term repair
HFSCs can participate in repairing damaged epidermis, but their relative long-term contributions remain unclear (9, 11, 12). We therefore challenged different skin SC niches to participate long-term in repairing damaged epidermis. To do so, we introduced shallow, intermediate, or full-thickness (deep) wounds during the prolonged resting phase (telogen) of the hair cycle (fig. S1, A and B). As judged by planar and sagittal immunofluorescence images, shallow wounds removed only epidermis, leaving behind both HF junctional zone (JZ) stem cells (Lrig1+) that rejuvenate the upper HF/sebaceous gland and bulge HFSCs (Krt24+) that fuel hair growth (fig. S1C). By contrast, intermediate wounds removed both epidermis and Lrig1+ JZ in their entirety, exposing the HF bulge with its characteristic hair shaft (green autofluorescence) to the overlying eschar (scab) of day-1 wounds (fig. S1C). With this documentation, we then used lineage-tracing reporter mice to monitor progeny of Krt24+ bulge HFSCs (Krt19CreER+Sox9CreER+) and Lrig1+ JZ HFSCs (Sox9CreER+Krt19CreERneg). Figure S2, A to C, shows the lineage tracings used and the fluorescence-activated cell sorting (FACS) strategy to purify and quantify SCs (α6+SCA1+) within re-epithelialized epidermis. These data were normalized to account for CreER targeting efficiency of Rosa26-fl-stop-fl-YFP within otherwise transcriptome-homogeneous bulge SCs (13).
In agreement with prior studies (9, 11) and underscoring the preferential mobilization of wound-edge EpdSCs in full-thickness wounds (14, 15), most deep wound re-epithelialization came from epidermis (Sox9YFPneg) and not HFs (Sox9YFP+) (Fig. 1A and fig. S2, D and E). However, in shallow and intermediate wounds, Sox9YFP lineage tracings showed >90% contribution, indicating little or no contribution from wound-edge EpdSCs (Fig. 1A and fig. S2F). By contrast, Krt19YFP tracings showed that bulge SCs were the major contributors to de novo repaired epidermis in intermediate wounds. There, the percentage of YFP+ HF-derived basal cells within re-epithelialized epidermis was comparable to that within the bulges of unwounded HFs, indicating that intermediate wounds use primarily bulge SCs in repair (fig. S2G). Notably, this percentage of Krt19YFP cells in de novo epidermis remained relatively steady even 2 months after wounding, indicating that long-term re-epithelialized epidermis in intermediate wounds was largely derived from bulge HFSCs (Fig. 1A). Together, these findings show that the type of wound affects which SC compartment will predominate in re-epithelialization: deep (EpdSCs), shallow (JZ HFSCs), and intermediate (bulge HFSCs).

HF-derived epidermis is fully functional in homeostasis
At day 60 after wounding, bulge SC–derived (hereafter, HF-derived) epidermis appeared similar to unwounded (native) epidermis in morphology (Fig. 1B), thickness (fig. S3A), and expression of terminal differentiation–specific markers (fig. S2, F and G). Although trans-epidermal water loss (TEWL) assays indicated a delay in skin barrier restoration for intermediate versus deep wounds, within 19 days after wounding, the eschar (scab) was gone, tight junctions were reestablished, and the skin barrier function of HF-derived epidermis was intact (Fig. 1C and fig. S3B). At this time, proliferation rates and orientations of EpdSC divisions were also comparable between HF-derived and native skins (Fig. 1D and fig. S3, C and D). Thus, in post-repaired homeostatic skin, the epidermis formed from wound-mobilized bulge HFSCs behaved much like native epidermis.
We next performed single-cell RNA sequencing (scRNA-seq) on FACS-purified live YFP+ and YFPneg keratinocytes from skins of age- and litter-matched control Krt19CreER; R26-YFP mice whose postnatal day 50 skins had been treated with tamoxifen during second telogen and then harvested either before wounding (day 0) or at day 80 after wounding (fig. S4A and table S1). After passing quality control metrics (fig. S4B), Leiden clustering of all epithelial cells at day 0 revealed six distinct subsets, readily identified by established markers (16, 17) and underscoring the marked differences between bulge HFSCs and EpdSCs (Fig. 1E and fig. S4, C to E).
Within the UMAP of EpdSC and suprabasal Epd clusters, day-80 bulge-derived wound-experienced Krt19CreER;R26YFP+ cells overlaid with day-80 unwounded native EpdSCs (Fig. 1E, lower inset). Among >16,000 detected genes, DESeq2 (18) revealed no significant differences in expression (Fig. 1E, right panels). These analyses showed that the transcriptomes of HF-derived and native EpdSCs were virtually indistinguishable; this was further confirmed by pseudospace analysis, where quantification of the relative position of each cell in the UMAP showed no difference (fig. S4F). Thus, despite originating from a markedly distinct niche in which they displayed a very different transcriptome and performed distinct regenerative tasks, bulge HFSCs (uniquely Krt19YFP+) responded to injury by re-epithelializing missing epidermis and establishing residence within a new niche. Given the comparable transcriptomes of HF-derived EpdSCs and native EpdSCs, it was not surprising that HF-derived EpdSCs were able to generate and maintain the skin’s barrier.

SCs adapt to their new niche after repair
Probing deeper for differences between HF-derived and native EpdSCs, we turned to the chromatin level. We performed ATAC-seq (assay for transposase-accessible chromatin using sequencing) (19) on FACS-purified HFSCs isolated during the repair and recovery process (Fig. 2A and fig. S5A). Principal components analyses revealed close clustering among replicates, the expected ATAC peaks centered around transcription start sites (TSSs) and CTCF chromatin looping sites, and peak distribution patterns were comparable across genomic regions (fig. S5, B to D). Genes associated with ATAC peaks were assigned first by proximity using the Genomic Regions Enrichment of Annotations Tool (GREAT) and then refined according to transcript status in bulge HFSCs, wound-activated SCs, and/or EpdSCs (12, 13, 17).
Before injury, >9000 chromatin peaks distinguished quiescent bulge HFSCs from homeostatic EpdSCs (Fig. 2B). Nearly 6000 peaks were lost within 3 days after wounding. Lost peaks associated with genes such as Nfix and Lhx2, which are essential for bulge identity and whose expression is known to be rapidly down-regulated upon wounding (Fig. 2, B to D, and table S2) (20, 21). We observed that 2360 HFSC peaks were maintained throughout the repair process before waning, including peaks associated with “lineage infidelity” genes such as Sox9 and the WNT receptor–signaling gene Tcf7l1 (TCF3) (22). Long after repair, however, ~1000 peaks established in bulge HFSCs remained open in HF-derived EpdSC chromatin (Fig. 2B). We address their importance later.
Although many HFSC chromatin peaks closed during and after wound repair, >11,000 chromatin peaks were gained. Associating with EpdSC-expressed genes, >4000 peaks were induced by day 3 after wounding, in agreement with the known induction of epidermal lineage markers by migrating bulge HFSCs during wound repair (“lineage infidelity”) (12) (Fig. 2, B and E). Even though the wound had closed by day 7, nearly 4000 homeostatic EpdSC-specific chromatin peaks only became accessible after the HF-derived EpdSCs had adapted to their new niche. Notably, these “niche-adaptive” peaks were associated with epidermal differentiation (Gata3, Krt10), barrier function (Scd1, Lce1l), and immune defense (Il20ra, Il34, Irf4) genes (Fig. 2C and table S3), reflecting the newly acquired tasks of HF-derived EpdSCs, namely skin barrier maintenance and sentinels to sense barrier breaches and recruit an immune response.
Some niche-adaptive peaks displayed heightened or even unique accessibility in HF-derived EpdSCs over that of native EpdSCs (Fig. 2F and fig. S6A). Because many of these peaks associated with epidermal differentiation and immune defense genes (Fig. 2C and table S4), their openness seemed to reflect “compensatory adaptation,” which suggests that HF-derived EpdSCs may remodel these genes to a more open chromatin state in order to transcribe the requisite levels needed to perform their new homeostatic tasks. In this way, HF-derived EpdSCs may be able to compensate for deficiencies in perceiving local environmental cues the way their native counterparts do. This notion was in agreement with the several days’ delay relative to native EpdSCs in establishing the barrier after re-epithelialization (Fig. 1C).

Epigenetic wound memories of inflammation and migration
The distinct post-repair configuration of a foreign stem cell’s chromatin relative to its native counterparts led us to wonder whether other kinds of epigenetic rearrangements might have occurred within these immigrants. Around 2500 peaks became accessible only after bulge HFSCs experienced the wound, but these peaks were retained in HF-derived EpdSCs after re-epithelialization (Fig. 3A and fig. S6A). Because these peaks were not a feature of native EpdSCs and their associated genes returned to baseline transcription levels after wound closure, they appeared to represent an epigenetic memory of the wound experience.
Wounds are known to induce inflammation, and this was reflected in the enrichment of immune-related genes and pathways for wound memory peaks (Fig. 3A and table S5). Topical skin treatment with imiquimod (IMQ), a psoriasis mimetic drug, induces long-lasting inflammatory-sensitive chromatin peaks within EpdSC chromatin independent of B and T lymphocytes and macrophages (23). We therefore looked at the ATAC peaks that were induced in IMQ-exposed EpdSCs and then retained at day 30 after inflammation (fig. S6B) and found them to also be more accessible in HF-derived EpdSCs (Fig. 3B). First described as “trained immunity” in short-lived macrophages, inflammatory memory peaks in EpdSC chromatin sensitize their associated genes to respond more robustly to secondary inflammatory assaults (24, 25).
In contrast, the large number of wound-induced memory ATAC peaks in our study went beyond those described for inflammatory memory. Many of these memory-associated genes encoded cytoskeletal organization and Rho signaling/actomyosin-regulated proteins essential for the extensive polarized cell migration that occurs during re-epithelialization (26–28) (Fig. 3A and table S5). This aspect of the wound response is not encountered when stem cells are simply exposed to a topical inflammatory stimulus.
To test the physiological relevance of this facet of the wound memory, we first examined the ability of wound-experienced HF-derived epidermis to respond to secondary wound closure in vivo. Relative to naïve native epidermis, wound-experienced HF-derived EpdSCs closed wounds markedly faster (Fig. 3C). As the rates of ethynyl-2′-deoxyuridine (EdU) incorporation were comparable (fig. S6C), the accelerated rate of wound closure appeared to reflect enhanced migration rather than proliferation. Indeed, at the wound front, the polarized epidermal tongue marked by migration-specific integrin α5 was significantly more pronounced in HF-derived epidermal wounds (Fig. 3D). Ex vivo assays further corroborated the enhanced migration of the HF-derived epidermis out of skin explants and revealed more robust signs of membrane ruffling, a hallmark of enhanced Rho/Rac-mediated actomyosin dynamics (fig. S6D) (29, 30).
To test whether the enhanced epithelial migration is intrinsic to HF-derived EpdSCs, we isolated day-80 HF-derived and native EpdSCs by FACS and adapted them to 2D culture conditions wherein no supporting cells were present (fig. S6E). Although no overt differences were observed during propagation (fig. S6F), when scratch-wounded in the presence of mitomycin C to block proliferation, HF-derived EpdSCs were significantly more efficient than their native EpdSC counterparts at migrating in and closing the gap (Fig. 3E and fig. S6G). Moreover, fluorescence microscopy and live imaging of wound-experienced green fluorescent protein (GFP)–actin+ HF-derived EpdSCs further corroborated features of active collective cell migration, including membrane ruffling, prominent integrin-mediated focal adhesions, and polarized actin at the leading edge of the migrating front (fig. S6, H and I, and movies S1 and S2) (26, 29, 30). These findings underscored the cell-intrinsic features of the enhanced performance of wound-experienced HF-derived EpdSCs unleashed upon exposure to a secondary wound environment.
Also noteworthy was the finding that although wound memory genes are largely transcriptionally silent during homeostasis, they became reactivated upon a secondary wound. Moreover, as shown by scRNA-seq analyses, the responsiveness of many memory genes was even more robust within the α5-integrin+ migrating tongue of secondary wounds (HF-derived EpdSCs) than in primary wounds (native EpdSCs) (Fig. 3F, fig. S7, and table S6). Further underscoring the potential importance of this finding, up-regulated genes were involved in actin cytoskeleton organization, Rho signaling, angiogenesis, tissue growth, and inflammation—all processes important for wound healing (fig. S7F and table S7).

HF-derived EpdSCs remember their niche origin
Curiously, 834 ATAC peaks were characteristic of homeostatic bulge HFSC chromatin and had maintained accessibility throughout both wound repair and epidermal adaptation of HF-derived EpdSCs (Fig. 4A). On the basis of this behavior and their relative paucity in native EpdSCs, these peaks in HF-derived EpdSC chromatin appeared to reflect an epigenetic memory of their bulge niche origin.
The top Gene Ontology term categories in this cohort were known to be important in HFSC and hair formation (table S8). These memory domains were associated with key genes regulating bulge HFSC quiescence, such as Bmp6, Sox9, Nfib, Col17a1, Id2, and Gli1, and also bulge-enriched genes essential for WNT signaling–mediated activation of HFSCs to launch the hair cycle, such as Tcf7l1, Lgr5, Lpp, Fzd1, and Lrp6 (22, 31–33) (Fig. 4A). Although these genes were not transcribed by HF-derived EpdSCs during normal skin homeostasis, we wondered whether this epigenetic memory might endow the cells with a heightened propensity to unleash an HFSC-like behavior when challenged to do so.
Through mechanisms poorly understood, HFSCs are known to possess higher WNT sensitivity than EpdSCs and display greater colony-forming efficiency and passaging potential when placed in culture under conditions that favor a lineage-activated (“hair germ”) state (4). True to their HF origins but distinct from native EpdSCs, HF-derived EpdSCs were more efficient in forming colonies (Fig. 4B). Additionally, in response to R-spondin1 and WNT7a, an established canonical WNT in regeneration and wound repair (33), HF-derived EpdSCs showed an increased sensitivity to WNT signaling, characteristic of activated HFSCs (Fig. 4C and fig. S8, A to C).
In their quiescent state, bulge HFSCs are in a WNT-restricted niche, where they express transcription factors (TFs) such as SOX9 and NFIB/X, which function in maintaining HFSCs and preventing EpdSC conversion (20). Bulge HFSCs also express TCF3, which is essential to mediate WNT-responsive HFSC activation to launch hair cycling (22). These genes were associated with memory of bulge niche origin, and when we switched to culture conditions that mimic quiescent bulge HFSCs in their native niche, HF-derived EpdSCs up-regulated the nuclear expression of SOX9, TCF3, and NFIB and down-regulated that of EpdSC TF KLF5 (Fig. 4D and fig. S8D). Thus, HF-derived EpdSCs retained their potential to respond to WNTs and bulge niche–promoting conditions even after taking up residence long term in the epidermal niche. Because the transcriptomes of HF-derived EpdSCs and EpdSCs within the homeostatic epidermal niche in vivo were indistinguishable, the thread connecting HFSCs and HF-derived EpdSCs appeared to be the >800 “memory of origin domains” associated with genes such as Tcf7l1, Nfib, and Sox9 that became reactivated upon secondary exposure to a bulge microenvironment.
We reasoned that if the collective “memory of bulge niche origin” is functionally important in skin physiology, it should confer an increased ability of HF-derived EpdSCs to regenerate hair in vivo. We tested this possibility with “chamber graft” assays in which we combined our cultured stem cell populations with neonatal dermal fibroblasts and engrafted them onto hairless (Nude) mice (4). In contrast to GFP-marked EpdSCs, which generated primarily epidermis in grafts, GFP-marked HF-derived EpdSCs formed both epidermis and hair efficiently (Fig. 4E and fig. S8E). These de novo hairs were replete with WNT-activated HF morphogenesis preceding hair production (fig. S8, E and F). Together, these data provide compelling evidence that HF-derived EpdSCs maintain a memory of their niche origin, and after engraftment in vivo, this memory endows them with a broadened tissue-regenerating capacity.

Insights into the multiplicity of memories, their diversity, and their ability to accumulate
We next turned to addressing how memory domains are established in space and time, and what maintains them in a poised open state when their environment shifts. We first analyzed the TF sequence motifs encompassed within ATAC peaks that are open exclusively in one particular state and then resolved in subsequent environments (Fig. 5A, fig. S9A, and table S9). As expected, the ATAC peaks exclusive to quiescent bulge HFSCs were enriched for motifs that bind bulge identity TFs (LHX2, NFIB/X, NFATc1, SOX9), wound-specific peaks were enriched for the binding of wound TFs (AP-1, CEBP, ATF/CRE, NFκB/REL), and EpdSC-specific peaks were enriched for the binding of EpdSC TFs (KLF5, GATA3, GRHL) (Fig. 5A). With the exception of the dual expression of EpdSC and HFSC TF genes during the lineage infidelity period of the wound response (12), the expression patterns of TFs correlated with this state specificity (Fig. 5B).
Displaying the greatest complexity in TF motifs were the memory-of-bulge-origin peaks; they showed enrichment for the binding of TFs from all three states (Fig. 5, A and B, and table S9). In this regard, they differed from lineage-infidelity ATAC peaks, which were typified by sequence motifs for HFSC-specific and wound-specific TFs but lacked those for EpdSC-specific TF motifs (fig. S9B and table S9). Although memory-of-bulge-origin ATAC peaks were already open in bulge HFSCs and remained open throughout environmental shifts, wound memory peaks opened de novo upon wounding. These domains were highly enriched for STAT and C/EBP, which have been implicated as key pioneer factors in opening inflammatory memory chromatin (34, 35). Like all of the memory domains analyzed, wound memory peaks were also typified by the presence of motifs for the chromatin remodeling factor AP-1 (36, 37) and EpdSC TF KLF5.
Finally, niche-adaptive and compensatory adaptation peaks were typified by the same binding motifs as native EpdSC peaks (Fig. 5A, fig. S9B, and table S9). Although knowledge of chromatin dynamics in the epidermis is still scant, this provides an explanation for how these peaks open only after re-epithelialization. These peaks show an enrichment for AP-1, which is also induced in a variety of stress situations. It will be interesting to probe whether adaptation entails a transient stress period during which these sites open and HF-derived EpdSCs adjust to their new environment.
Another intriguing feature of memory domains is that they tended to be within intergenic regions, likely enhancers, that are more open than their state-specific counterparts at the time when their establishment took place (fig. S9, C and D). We posit that this might extend the accessibility of these domains as stem cells transition from one environment to another, thereby facilitating the binding of new state-specific TFs as the expression of others wanes and maintaining the open chromatin state across these dynamics.
Our transcriptomic analysis of primary versus secondary wounds (Fig. 3F) suggests that the memory and niche adaptation domains identified here act similarly to inflammatory memory domains in that their associated genes are generally dormant unless these memory enhancers are triggered by reexposure to the stimulus that prompted their activation (25). To further bolster this point and establish a more direct connection between memory peaks and their environment-sensing activity, we selected some of the most accessible memory/adaptation domains and used them as enhancers to drive eGFP in mice (fig. S9E). These domains functioned faithfully in driving reporter expression in a state-specific manner that corresponded both to the cell’s environmental experience within the skin and to the natural transcriptional status of the genes associated with these chromatin peaks (Fig. 5C and fig. S9, F to H). Although we could not appraise the heightened accessibility of compensatory adaptation peaks by this assay, we did find that peaks associated with compensatory adaptation were active in only epidermis and not HF or wound states. Overall, these results underscore the physiological relevance of memory/adaptation domains in gene governance.
Finally, if these epigenetic memories are truly dependent on a stem cell’s prior experiences, then memories of wound and bulge niche origin should be uncoupled in re-epithelialized epidermis derived from wound-mobilized EpdSCs (deep wound) as compared to wound-mobilized bulge HFSCs (intermediate wound). Indeed, by ATAC-seq analyses, whereas wound memory domains were comparable, epidermis from wound-mobilized EpdSCs lacked robust HFSC memory domains (Fig. 5D and fig. S10, A to E). Chromatin changes associated with compensatory adaptation to the epidermal niche were also markedly different in HF-derived EpdSCs versus epidermal-derived EpdSCs (fig. S10F).
Consistent with their wound memory, wound-experienced epidermal-derived EpdSCs exhibited enhanced performance over their unwounded counterparts at wound closure both in vivo and in vitro (Fig. 5E and fig. S10, G and H). However, and in contrast to wound-experienced HF-derived EpdSCs, the sensitivity of wound-experienced epidermal-derived EpdSCs to WNTs and to hair regeneration in engraftments was no better than their native EpdSC counterparts (Fig. 5F and fig. S10I). These data provide functional evidence in support of the experience-dependent and cumulative nature of epigenetic memories and their physiological relevance in conferring heightened secondary responses to subsequent exposure to the environmental conditions that stimulated their establishment.

Discussion
Because long-lived adult stem cells can harbor inflammatory memory (23, 25, 34), we examined whether they might also possess memories of different kinds of experiences, and if so, whether such memories are cumulative. Wounding under conditions where bulge HFSCs were the primary responders to long-term epidermal repair offered an opportunity to monitor these tissue stem cells as they encountered a series of diverse experiences, first exiting their HF niche, then encountering inflammation as they migrated into the wound bed, then re-epithelializing the vacated epidermis, and eventually taking up long-term residence there. A remarkable facet of this choreographed process is that these SCs gained epigenetic memories of their diverse experiences at each step along the way.
These memories were both distinct and cumulative. For each environmental encounter, certain chromatin domains within key temporally regulated gene enhancers became highly accessible and then remained accessible long after the experience and after normal physiology had been restored. Moreover, each memory that had been acquired and stored along the way endowed the stem cell with distinct physiological advantages that were dormant during normal homeostasis, but then become unleashed upon a secondary challenge. Thus, HF-derived EpdSCs now residing as permanent immigrants in the epidermis behaved like their native EpdSC counterparts, right down to their equivalent transcriptomes. However, when roused to self-renew, regenerate skin tissues, or provide a secondary wound response, wound-experienced HF-derived EpdSCs outperformed normal EpdSCs. These features are classical ones of trained immunity, where some epigenetic marks of an inflammatory response are retained after exposure, priming the cell to respond more quickly to a subsequent exposure (24, 25).
The epigenetic memories revealed here are separable, varying according to a stem cell’s past experiences. Thus, when EpdSCs were mobilized to repair deep wounds, they retained a wound memory similar to HFSCs that were mobilized to repair intermediate wounds. However, only HF-derived EpdSCs possessed a memory of their HF origin. Moreover, although EpdSCs underwent some chromatin remodeling after their re-epithelialization of a deep wound, this adaptation differed markedly from the compensatory adaptation exhibited by HF-derived EpdSCs. These findings further support the notion that sustaining a transcriptome akin to native EpdSCs and performing the task of fueling the skin’s barrier remained an adaptation challenge for these former HF residents in a way not experienced by wound-mobilized EpdSCs.
Finally, we showed that these epigenetic memory domains possess characteristics that are adjusted to suit their particular past experiences. In cases such as Sox9, whose TF is required for maintaining HFSC fate as well as lineage infidelity, some ATAC peaks activated in wound-mobilized bulge SCs closed within 3 days, others closed only after re-epithelialization, whereas still others remained open at day 80 after epidermis repair. Sox9’s transcriptional profile—more robust in the bulge than in the wound—reflected these features. However, as is common for all genes associated with day 80 memory-of-origin domains, Sox9 was not expressed in HF-derived EpdSCs until spurred to do so by the right environmental stimuli. In this regard, it was notable that “memory of bulge niche origin” domains possessed motifs for the binding of HFSC, wound, and EpdSC TFs. This enabled the chromatin to adapt to different lineage fates, even after the TFs that initially drove their transcription were silenced.
Repairing injuries is a crucial role for all tissue stem cells. As such, our discovery that stem cells acquire and store functional epigenetic memories of the complex steps involved in wound repair is likely to have implications for tissue fitness and regenerative medicine that go beyond the skin. Our findings imply that epigenetic memories can be maladaptive (as inflammatory memory is likely to be for chronic inflammatory conditions) or beneficial (as the memories of origin, wound repair, and even niche compensatory adaptations that we describe here would appear to be). Although these memories are unlikely to persist forever, it is noteworthy that 2 months in a mouse’s life is the equivalent of 5 to 6 years in humans. Thus, the combination of epigenetic adaptations and cumulative acquisition of memories stockpiled within the chromatin of wound-experienced adult stem cells provides an altered view of tissue performance. That said, as desirable as it may be to be able to self-renew longer, repair wounds faster, and broaden the fates of tissue stem cells, these features are also ones associated with malignancies, and as such, they may blur the line between good and bad memories.

Methods summary

Deep, intermediate, and shallow wounding
Adult mice at telogen were anesthetized with isoflurane and treated with buprenorphine prior to wounding. For shallow and intermediate wounds, back skins were first shaved, then treated with hair removal cream (Veet). A Dremel Inc. tool with a polishing wheel was used to generate abrasions by polishing the skin laterally three or four times for shallow wounds and six to eight times for intermediate wounds (Dremel 100-series rotary tool and 520 polishing wheel). For deep wounds, back skins were shaved and wounds were created with a 6-mm punch biopsy tool.

ATAC-seq library preparation and sequencing
ATAC-seq libraries were made from freshly FACS-sorted cells, with two to four biologically independent replicates per cell population. Library preparation was performed as described (12). The samples were sequenced on Illumina Nextseq 500 using a 75-bp paired-end-reads setting.
Detailed materials and methods are available in the supplementary materials.

Supplementary Material
Figures, Tables, ReferencesMovies S1, S2