The exercise cytokine interleukin 15 rescues slow wound healing in aged mice

Impaired wound healing in elderly individuals increases infection risk and prolongs surgical recovery, but current treatment options are limited. Low doses of interleukin 15 (IL-15) that mimic exercise responses in the circulation improve skin structure and increase mitochondria in uninjured aged skin, suggesting that IL-15 is an essential mitochondrial signal for healing that is lost during aging. Here we used gene microarray analysis of old and young murine epidermal stem cells and demonstrate that aging results in a gene signature characteristic of bioenergetic dysfunction. Intravenous IL-15 treatment rescued chronological aging-induced healing defects and restored youthful wound closure in old, sedentary mice. Additionally, exercise-mediated improvements in the healing of aged skin depend upon circulating IL-15. We show that IL-15 induces signal transducer and activator of transcription 3 (STAT3) signaling characteristic of young animals, reduces markers of growth arrest, and increases keratinocyte and fibroblast growth. Moreover, exercise or exercise-mimicking IL-15 treatment rescued the age-associated decrease in epidermal mitochondrial complex IV activity. Overall, these results indicate that IL-15 or its analogs represent promising therapies for treating impaired wound healing in elderly patients.

The integrity of the skin barrier is critical for blocking pathogens and other environmental hazards and the proper maintenance of skin structure requires timely, coordinated tissue repair.
Skin healing is complex and requires the coordination of clotting, the recruitment of multiple types of immune cells and significant amounts of tissue re-growth (1). Skin trauma is repaired relatively quickly in young skin, however in the elderly wounds heal more slowly (2) and this older population is particularly susceptible to nonhealing, chronic wounds (3,4). Delays in skin healing dramatically increase the risk of infection, prolong recovery from surgery and increase morbidity and mortality in older adults (5). As a result, the care and management of non-healing ulcers or poorly healing wounds costs approximately $10-25 billion per year in the U.S. alone (6). Overall, improving skin healing in the elderly is a major challenge for medicine and new treatment strategies are needed.
Aging is broadly associated with various forms of cellular dysfunction such as genomic instability, inflammation and stem cell dysfunction (7,8). While the mechanisms responsible for poor healing remain unclear, there is evidence that aging delays epidermal stem cell activation (9,10), shifts immune cell populations (11) and alters vascular repair after injury (12). In association with these changes aged skin exhibits epidermal and dermal thinning, has lower rates of post-injury keratinocyte proliferation, and experiences slowed wound reepithelialization. Past work has employed expression profiling of aged keratinocytes to better understand the mechanisms of epidermal dysfunction and these data suggest that changes in the cell cycle, disrupted cytokine signaling and altered growth cues may prevent rapid healing (9). As a result of this work, keratinocyte STAT3 was Interleukin-15 rescues impaired wound healing in aged skin 2 identified as a promising target for reversing ageinduced healing defects as it is required for reepithelialization after skin injury (13) and dysfunctional STAT3 signaling at the wound edge is thought to impede keratinocyte proliferation in old mice (9). However, we do not understand what triggers age-related changes in keratinocyte STAT3 signaling over the lifespan and we currently lack treatments that can restore this pathway and augment healing in aged skin.
Several therapeutic approaches have been identified that improve skin healing and reduce scarring in young, healthy skin, such as negative pressure therapy, cellular scaffolds or supplemental growth factors, however none of these therapies has proven effective at reversing healing defects in aged skin (14). There are several roadblocks to developing effective healing treatments for the elderly: 1) few studies specifically target intrinsic, aging-induced cellular dysfunction; 2) aged skin cells have impaired signal transduction responses (15), which may reduce the effectiveness of growth factors or other mitogenic treatments in aged skin; 3) excessive exposure to exogenous growth factors may cause the malignant transformation of cells or promote tumor growth in aged, cancer-prone cells. Thus, it may be more prudent to first focus on treating the pathophysiology of skin aging in order to curtail healing defects. In a similar vein, we reasoned that therapies known to successfully reverse skin aging would also improve poor skin healing.
Aerobic exercise is a robust lifestyle intervention that produces widespread health benefits over the lifespan including reduced mortality and a lower incidence of neurologic, cardiovascular and cancer-related disease (16). Notably, pre-and post-operative exercise is associated with improved recovery from various types of surgery in elderly adults (17). Exercise is also one of the few interventions shown to enhance wound healing in the aged skin of both mice (18) and humans (19), however the mechanism underlying these benefits is not known. Since exercise has such a well-established safety profile to protect against systemic disease, mimicking exercise signaling is an ideal way to treat agerelated disorders, particularly in the frail elderly that are unable to exercise. We previously discovered that regular exercise reverses some of the structural effects of epidermal and dermal aging in the skin and increases the abundance of skin mitochondria in humans and mice (20). By screening the blood of athletes, we identified exercise-induced IL-15 as the circulating signal that increases skin mitochondria. Additionally, we found that a low, exercise-mimicking dose of recombinant mouse IL-15 (rmIL-15; 25 ng/kg) provided daily for 5 weeks to old mice was sufficient to improve skin structure, restore tissue mitochondrial abundance and reduce circulating inflammatory cytokines akin to exercise (20). This was of interest because IL-15 is considered proinflammatory and induces T-cell proliferation (21). Conversely, regular exercise is broadly antiinflammatory (22), but stimulates far lower levels of IL-15 in circulation than the concentrations associated with innate immune responses (20). Overall, these data suggested that increasing mitochondrial abundance in the skin is a potential way to treat skin aging and that the low pulses of IL-15 stimulated by exercise are a physiologic regulator of skin mitochondria. As circulating levels of IL-15 are reduced in aged mice (23) and disrupted IL-15 signaling is observed in nonhealing human wounds (24), we reasoned that our low-dose rmIL-15 therapy may restore youthful signaling and improve healing in aged animals.
We show that daily, low-dose injections of rmIL-15 restore youthful wound closure, increase mitochondrial activity, reduce cellular senescence and normalize STAT3 signaling in aged, sedentary mice. Additionally, we find that the improvement in wound re-epithelialization observed in old, exercise-trained mice is dependent on IL-15 in circulation, demonstrating that IL-15 is an important anti-aging exercise signal.

Interleukin-15 rescues impaired wound healing in aged skin
and dendritic T-cells (9,10) in skin aging and the resulting impairment in re-epithelialization. To gain additional insight into the underlying mechanisms of aging within the epidermis, we reanalyzed a publicly available gene expression dataset comparing FACS sorted young and old murine epidermal stem cells (ESCs, CD49f high /CD34 -) collected from uninjured dorsal skin that originally examined age-related changes in circadian rhythms (GSE84511) (28). We analyzed differentially expressed genes in young versus old ESCs irrespective of the time of day (n = 24 per age group), giving us substantial power to detect aging effects, and then employed two established differential gene expression analysis tools, GAGE (generally applicable gene set enrichment for pathway analysis) (29) and GSEA (gene set enrichment analysis) (30). GAGE evaluates significantly differentially expressed genes without an arbitrary fold-change cut-off to identify sets of coordinately perturbed genes (up and/or down) and their related pathways, increasing the sensitivity and reproducibility of pathway analysis (29). Our GAGE analysis identified several significantly altered pathways previously associated with ESC aging such as JAK/STAT signaling and senescence, as shown by simultaneous aging associated elevations in canonical STAT signaling with prosurvival (Akt), cell cycle arrest (Cyclin D1), and anti-apoptotic signaling (Bcl-Xl) (Fig. 1A, S1A, S1B). Additionally, GAGE highlighted pathways that have never been connected to ESC aging, including the cellular energy sensor AMP-activated protein kinase (AMPK) and metabolic peroxisome proliferator-activated receptor (PPAR) signaling (Fig. 1A). Pathview is a graphical display package which overlays gene expression metrics on Kyoto Encyclopedia of Genes and Genome (KEGG) pathway diagrams. Detailed networks of the JAK/STAT pathway and AMPK signaling alterations using Pathview show comprehensive downstream pathway dysregulation (up and down) in response to aging (Fig. S1A-B). Subsequently, GSEA analysis was used to identify cellular functions consistently up-regulated or consistently down-regulated in response to aging. We found that metabolic pathways related to glucose and carbohydrate metabolism were upregulated and that mitochondrial synthesis and ATP production were downregulated in aged ESCs (Fig. 1B), suggesting a shift away from oxidative phosphorylation-derived cellular energy. A selection of differentially expressed genes from the identified pathways highlights these changes (Fig.  1C). Age-induced mitochondrial deterioration is well-described in other stem cells such as muscle and neural progenitors (31), but these data provide the first evidence that the epidermal stem cell compartment experiences broad metabolic dysfunction during aging. Together, these data suggest that mitochondrial defects are mirrored across multiple cell types in aging skin and may drive its aging process.

Administration of IL-15 to sedentary mice rescues defective would closure
We have previously shown that IL-15 reverses skin aging in intact skin, but the impact of low dose rmIL-15 treatment on skin healing has not been established. Interestingly, IL-15 receptor a subunit expression (IL-15Ra) was identified through our pathway analysis to be increased in old versus young epidermal stem cells (Fig. 1C), which is analogous to greater IL-15Ra expression levels found in non-healing versus healing wounds in humans (24). Moreover, Keyes et al. identified IL-15 as one of a handful of downregulated genes in old versus young mouse epidermal keratinocytes isolated from the leading wound edge (9) and circulating IL-15 levels are reduced in old mice (23). In combination this suggested that IL-15 treatment may improve wound healing in aged skin. To assess this, we used our previously described regimen of IL-15 therapy in which aged mice were pre-treated for 33 days with daily intravenous rmIL-15 (old rmIL-15, 25 ng/kg) or PBS vehicle injections (old controls) as well as a young, PBSinjected control group (20). After this pre-treatment phase, excisional wound surgery was performed and the closure of wounds was tracked over 9 days, which is the approximate time to full wound closure in young mice. Splinted, full-thickness wounds were utilized to prevent confounding wound contraction prevalent in murine wounds and to promote greater re-epithelialization akin to human wound healing as described by others (32). Due to the long wound healing period after pre-treatment and the rapid turnover rate of the murine epidermis (33), we continued daily injection treatments of rmIL-15 or vehicle control until 5 days post-surgery ( Fig. 2A). The young control group was largely fully healed after 9 days (> 95% original wound Interleukin-15 rescues impaired wound healing in aged skin 4 area healed), whereas old control mice had significantly less advanced wound closure at 9 days post-wounding ( Fig. 2B-D. Remarkably, we found that rmIL-15 treatment in old mice resulted in significantly smaller wounds at 9 days postwounding compared to old controls and their degree of healing was indistinguishable from young control mice (Fig. 2B-D). We also analyzed the rate of wound closure over time and this showed a significantly decreased wound closure rate in the old PBS group relative to young PBS, but a restoration of healing rate in rmIL-15 treated old mice (Fig. 2E). Additionally, when we extrapolated these healing rates to time until full wound closure, we found that old control mice were projected to heal approximately 2 days slower than both young control and old rmIL-15 treated groups (Fig. 2F).
It has been recently shown that a transiently thickened region of Keratin 14 positive epidermal progenitors forms inward from the original wound margin, termed the proliferative hub (34). This thickened hub sustains the extending epithelium as the wound closes but resolves and thins when reepithelialization is completed. This phenomenon has not yet been investigated in aged mouse skin, so we sought to use it to better understand the impact of our therapy on the extent of reepithelialization. Thus, we stained day 9 postwounding skin cross-sections for Keratin 14 and 5ethynyl-2'-deoxyuridine (EdU) incorporation and assessed the average thickness of the Keratin 14 layer 300 microns toward the center of the wound from the original wound margin guided in part by EdU staining. We found that old control mice had a hyper-thickened proliferative hub at day 9 postwounding compared to young controls, suggesting a less advanced stage of re-epithelialization, however, the old rmIL-15 treated group had a reduced proliferative hub thickness compared to old controls and this thickness was indistinguishable from young mice ( Fig. 2G-H). Thus, short-term IL-15 therapy is able to restore a youthful rate of skin healing in aged mice.

Recombinant IL-15 therapy rescues ageassociated changes in STAT3 and cellular senescence in vivo
Since STAT3 signaling is implicated in keratinocyte aging (9), we sought to clarify how it was impacted by aging and IL-15 therapy in vivo. Using the intact skin excised during wound surgery, we found that staining for nuclear STAT3 phosphorylated at tyrosine 705 (pSTAT3 Y705 ) was higher in interfollicular epidermal cells and dermal fibroblasts (vimentin + dermal cells) of old control mice compared to young controls, however, rmIL-15 treatment largely abrogated this elevation in both cell types (Fig. 3A-B). Conversely, STAT3 phosphorylation at its S727 site (pSTAT3 S727 ) was reduced in old versus young control mice yet restored in old rmIL-15 treated mice (Fig. 3A, C). Thus, the two STAT3 regulatory sites are altered in opposite directions by both aging and rmIL-15 treatment within the epidermis and dermal fibroblasts. Since we have previously shown our rmIL-15 regimen increases mitochondrial activity in whole mouse skin and within cultured human fibroblasts (20), we sought to assess whether these changes occurred in isolated epidermis. Thus, we measured complex IV activity in whole epidermal protein lysates from each group. We found that complex IV activity tended to be reduced in young versus old control mice (p = 0.097, Fig. S2A). Additionally, old mice treated with rmIL-15 had significantly greater complex IV activity than old control mice (Fig. S2A). These data show that rmIL-15 treatment augments epidermal mitochondrial function. Cellular senescence has also been implicated as a driver of aging in numerous tissues (35) and biomarkers of senescence such as the loss of nuclear proteins HMGB1 or Lamin B1 occur more frequently in aged skin (36)(37)(38). Thus, an increased burden of growth arrested cells is a logical explanation for at least some of the ageassociated delays in cell growth during wound healing. Given the improved cellular re-growth after injury found in old rmIL-15-treated mice, we reasoned that this treatment would be associated with a lower burden of cellular senescence. To assess senescence, we quantified the presence or absence of nuclear HMGB1 and nuclear Lamin B1 levels via immunofluorescence in intact skin samples collected during wound surgery. We found that nuclear HMGB1 staining was reduced in epidermal cells and dermal fibroblasts of old control mice relative to young control, but that this loss was not observed in the old rmIL-15 treated mice ( Fig. S2B-C). However, the aging-associated loss of nuclear Lamin B1 observed in old control mice compared to young controls was only partially rescued by rmIL-15 treatment in the epidermis and Interleukin-15 rescues impaired wound healing in aged skin was not impacted by rmIL-15 therapy in dermal fibroblasts (Fig. S2B, D). The difference in response between senescence biomarkers is likely due to HMGB1 being more related to transient growth arrest (p53-related), whereas Lamin B1 is more reflective of terminal growth arrest (Rb/p16 ink4a ). Taken together, IL-15 therapy reduces the proportion of senescent cells in aged skin.
IL-15's most well-established physiologic function is in innate immunity where it activates and expands T-cells and NK cells and this function has been leveraged as an immunotherapy to reduce the growth of some types of cancers (21). However, IL-15's effects on T-cell proliferation and its antitumor activity typically require at least 3 orders of magnitude higher dosing than the dose we utilized in our experiments (39). Nevertheless, because Tcells play a prominent role in wound healing (9), we wanted to assess the impact of our recombinant IL-15 treatment on this immune population. First, we assessed spleen weights across all of our treatment groups to ascertain any broad changes in immune system activation. Old mice had greater spleen weights (as a percent of body mass) compared to young mice, but this was similar between old control and old rmIL-15 treated mice (Fig. S2E). We then assessed T-cells in the epidermis by staining for the T-cell marker CD3 in day 33 intact skin collected during wound surgery. However, we found no differences between our treatment groups in the proportion of CD3 + cells within the epidermis (Fig. S2F). Finally, we quantified circulating, peripheral blood cytotoxic T cells (CD3 + /CD8 + ) by flow cytometry in our experimental cohort after 31 days of pre-treatment and just prior to wound surgery. Similar to epidermal T-cells, we found no changes in the abundance of circulating cytotoxic T-cells across our treatment groups (Fig. S2G). Thus, our dose of IL-15 does not appear to expand circulating or resident epidermal T-cells or overtly activate the immune system in aged mice.
To gain further insight into how IL-15 was working, we next sought to examine whether this benefit was specific to aged skin. We repeated our experiment of 33 daily injections of rmIL-15 or PBS vehicle prior to wound surgery in young mice. We presumed that the rate of re-epithelialization and healing was somewhat maximal in young skin. In agreement with this hypothesis, we found that rmIL-15 treatment in young mice did not impact day 9 wound closure or the rate of wound healing compared to young control mice ( Fig. S3A-D). However, there was a modest trend for a reduction in projected time until closure (< 1 day, p = 0.0502, Fig. S3E). This suggests that young skin is able to "catch up" at later stages of healing which somewhat negates the benefits of IL-15 therapy. Since IL-15 treatment impacted pSTAT3 signaling in aged skin, we again measured the proportion of skin cells with strong nuclear pSTAT3 Y705 and cellular pSTAT3 S727 in tissue cross sections. We found that rmIL-15 treatment in young mice only produced significant changes in interfollicular epidermal cells, but not dermal fibroblasts, reducing nuclear pSTAT3 Y705 positive cells and increasing pSTAT3 S727 positive cells compared to young vehicle controls ( Fig. S3F-I). Thus, recombinant IL-15 treatment does not appreciably impact healing in young mice and modestly alters STAT3 signaling in the epidermis.

IL-15 directly increases the cell growth of mouse fibroblasts and keratinocytes
Our data clearly show that low-dose IL-15 improves the pathophysiology of aged skin and its healing, but it was still unclear whether this was a direct, cell autonomous effect. We have previously shown that post-exercise levels of recombinant human IL-15 (10 pg/mL) stimulate human dermal fibroblast proliferation (20), so we mirrored these conditions using post-exercise levels of IL-15 found in mice (100 pg/mL). To examine whether low-dose IL-15 acts directly on mouse keratinocytes, we employed an explant culture method that recapitulates age-related growth defects in vitro and exposed young and old skin explants to vehicle or 100 pg/mL of rmIL-15 immediately after plating. After 48 hours of treatment, young mouse keratinocytes had a larger average area of outgrowth than old keratinocytes and rmIL-15 treatment was able to rescue growth defects in old explants to that of young mice ( . This demonstrates that rmIL-15 treatment rescues growth defects in aged skin explants primarily by augmenting epidermal stem cell growth. Finally, since we observe many parallel effects on both keratinocytes and fibroblasts in vivo, we sought to test the effect of 100 pg/mL IL-15 on mouse fibroblast growth. Primary mouse dermal fibroblasts showed a significant increase in proliferation when exposed to rmIL-15 for 72 hours and this effect is blocked when the cells are treated with the JAK inhibitors ruxolitinib or tofacitinib prior to rmIL-15 ( Fig. 3H-I). Thus, low dose IL-15 acts cell autonomously on mouse epidermal stem cells and fibroblasts to stimulate cell growth in a JAK-dependent manner. In a separate experiment, we measured the acute response of pSTAT3 Y705 and pSTAT3 S727 relative to total STAT3 in fibroblasts exposed rmIL-15 via western blotting. Since the cytokine IL-6 is well known to signal via STAT3 Y705 phosphorylation, we used this as a reference for comparison. IL-15 treatment only increased pSTAT3 S727 , whereas IL-6 only induced pSTAT3 Y705 in comparison to untreated cells ( Fig.  3J-L). Therefore, unlike the classical STAT3 activator IL-6, low-dose IL-15 increases pSTAT3 S727 and not pSTAT3 Y705 .

Exercise improves wound healing in aged mice via IL-15 signaling
Exercise is one of the only interventions that has been shown to improve wound healing in both aged mice (18) and humans (19), however the underlying mechanism of this effect is not known. Since our rmIL-15 therapeutic regimen mimics physiologic elevations of post-exercise IL-15 in circulation, we sought to test whether this signaling was essential for the benefits of exercise on wound healing in aged skin. Thus, we subjected 2 groups of old mice to daily treadmill exercise for 33 days, with one group receiving IL-15 neutralizing antibody treatment (Old EX-nAb) and the other receiving control IgG antibody (Old EX) prior to each exercise session. Additionally, we included young and old sedentary control groups (Young SED, Old SED) that also received daily IgG control antibody injections. After 33 days of pre-treatment, we again performed splinted wound surgery and monitored healing over 9 days. To sustain our therapeutic window and in order to be congruent with the prior work showing that exercise improves wound healing in old mice (18), we continued the exercise and injection treatments for 5 days postsurgery (Fig. 4A). As might be expected, the continuation of exercise treatment after wounding surgery produced significantly more variability in day 1 to 5 of healing than in sedentary mice. Despite this, we found that aging resulted in impaired healing by day 9 compared to young sedentary controls, yet exercise treatment reversed this effect in old exercised mice versus old sedentary mice ( Fig. 4B-D). Remarkably, old exercised mice receiving IL-15 neutralizing antibody had similar day 9 wound closure as the sedentary old mice, demonstrating that IL-15 signaling is essential for the effects of exercise to speed wound closure in aged skin ( Fig. 4B-D). There was no effect on the rate of wound healing across the 4 groups (Fig. 4E). Analyzing the projected time until full wound closure revealed a significant 1.5 day difference between young, sedentary control and old sedentary mice groups, but projected healing times for old exercised and old exercised with neutralizing antibody groups were not significantly different compared to the young group (Fig. 4F). We then measured the Keratin-14 positive epidermal proliferative hub using day 9 post-wounding skin tissue cross-sections. We again found that old sedentary mice had a significantly thickened proliferative hub compared to young sedentary mice, however this was reduced to the thickness of young mice in exercised old mice (Fig. 4G-H). Importantly, old exercised mice that received IL-15 neutralizing antibody had greater proliferative hub thickness than old exercised and young sedentary mice ( Fig. 4G-H), further supporting the necessity of IL-15 for the beneficial effects of exercise on healing. These findings corroborate the work of others showing that exercise improves wound closure and re-epithelialization in old mice and we extend these results by showing that it is dependent on IL-15 in circulation.
Since recombinant rmIL-15 therapy in sedentary mice corrected age-induced changes in STAT3 signaling, we sought to examine this and additional aspects of mitochondrial and inflammatory signaling within our exercise experiment. Using intact skin excised during wound surgery, we again stained and quantified nuclear pSTAT3 Y705 -and pSTAT3 S727 -positive epidermal and dermal fibroblast cells. We found that old sedentary mice had more pSTAT3 Y705 -Interleukin-15 rescues impaired wound healing in aged skin positive cells than young mice but that this was abrogated in old exercised mice (Fig. 5A-B). Additionally, old exercised mice receiving IL-15 neutralizing antibody had more pSTAT3 Y705positive cells that young sedentary mice or old exercised mice (Fig. 5A-B). Consistent with the inverse relationship of pSTAT3 S727 -positive cells to those with pSTAT3 Y705 , we found fewer pSTAT3 S727 -positive cells in old sedentary versus young sedentary mice (Fig. 5A, C). In addition, old exercised mice had more pSTAT3 S727 -positive cells than old sedentary animals, and this effect was lost in old exercised mice receiving IL-15 neutralizing antibody (Fig. 5A, C). Thus, exercise reverses aging-induced changes in STAT3 signaling in skin and this is dependent upon circulating IL-15.
Given the dual role of the Y705 and S727 phosphorylation sites in regulating inflammatory and mitochondrial signaling, respectively, we sought to characterize downstream epidermal gene expression of these processes in response to aging and exercise. To capture in vivo gene expression changes in response to our treatment, we used whole epidermis separated from intact skin using dispase enzymatic digestion. Consistent with the downregulation of mitochondrial signaling observed in aged ESCs (Fig. 1), we found lower mRNA expression of the mitochondrial complex I subunit Ndufs2 in old sedentary mice compared to young sedentary mice (Fig. 5D). Interestingly, epidermal Ndufs2 mRNA expression in old exercised mice was indistinguishable from young sedentary mice, but remained lower in old exercised mice receiving IL-15 neutralizing antibody (Fig.  5D). Additionally, while Cox7B mRNA was not altered in young sedentary versus old sedentary mice, it was significantly elevated in old exercised mice compared to young sedentary mice (Fig. 5D). However, old exercised mice treated with IL-15 neutralizing antibody exhibited no difference in Cox7B mRNA compared to young or old sedentary mice. There were no significant changes in the complex III subunit Uqcr10 or the mitochondrial regulator Pgc1a between treatment groups (Fig.  5D). To corroborate these mitochondrial changes, we again measured Complex IV activity in epidermal protein lysates. We found that old sedentary mice, but not old exercised mice, had significantly lower complex IV activity than young sedentary mice (Fig. 5E). Additionally, old exercised mice that received IL-15 neutralizing antibody had significantly reduced complex IV activity compared to young sedentary mice and old exercised mice (Fig. 5E).
When assessing epidermal inflammatory genes, we found that old exercised mice exhibited significantly higher Cxcl1 mRNA than young control or old control mice, but old exercised mice receiving IL-15 neutralizing antibody had comparable expression to young and old control mice (Fig. 5F). Similarly, old exercised mice, but not exercised mice receiving IL-15 neutralizing antibody, had greater epidermal Cxcl2 expression than old control mice (Fig. 5F). We also assessed the mRNA expression of the pro-inflammatory cytokine IL-6 (Il6), a classical activator of pSTAT3 Y705 signaling, where we expected an ageinduced elevation in Il6 expression as seen in ESCs (Fig. 1C). Unexpectedly, we found that epidermal Il6 was reduced in old versus young control mice, but was restored in old exercised mice to that of young control mice (Fig. 5F). Moreover, old exercising mice given neutralizing antibody had significantly lower Il6 than young controls (Fig.  5F). Finally, we measured Il1b mRNA, but found no differences between treatment groups (Fig. 5F). Collectively, these data demonstrate that exercise depends upon circulating IL-15 to reverse aginginduced changes in epidermal mitochondrial and cytokine signaling. Finally, since we have shown that exercise and IL-15 therapy each independently reverse age-induced epidermal thinning, we sought to test if circulating IL-15 was essential for this effect of exercise in old mice. As we have previously observed, old control mice had reduced epidermal thickness versus young controls, but this was restored in old, exercised mice ( Fig. 5G-H). However, old exercising mice receiving IL-15 neutralizing antibody had significantly lower epidermal thickness that old exercised mice ( Fig.  5G-H), showing that exercise-induced IL-15 signaling is essential for the restoration of youthful epidermal structure.

Exercise-mediated changes in skin cell senescence partially depend upon circulating IL-15
Since recombinant IL-15 therapy was able to reverse the pattern of epidermal and dermal senescence in aged skin, we reasoned that exercise training in old mice would produce similar effects. Thus, we again assessed indicators of skin cell Interleukin-15 rescues impaired wound healing in aged skin senescence (i.e. loss of nuclear HMGB1 or Lamin B1) in our exercise cohort. We found that the number of cells with nuclear HMGB1 in epidermal cells and dermal fibroblasts was again reduced in old versus young sedentary mice (Fig. S4A-B). This loss of nuclear HMGB1 seen with aging was rescued in old exercised mice to that of young sedentary mice, however the rescue was completely abrogated by IL-15 neutralizing antibody treatment (Fig. S4A-B). Additionally, the proportion of epidermal cells and dermal fibroblasts with positive staining for nuclear Lamin B1 levels was reduced in old versus young sedentary mice but was restored to young levels in old exercised animals (Fig. S4A,  C). Old exercised mice receiving IL-15 neutralizing antibody had fewer dermal fibroblasts with nuclear Lamin B1 compared to old exercising mice, yet the proportion of epidermal cells with nuclear Lamin B1 was similar in old exercised mice treated with or without neutralizing antibody (Fig. S4A, C). Additionally, there were no significant changes in the epidermal mRNA expression of the cell cycle regulators p21 or p16 ink4a in response to aging or exercise (Fig. S4D). Therefore some, but not all, of the benefits of exercise to mitigate skin growth arrest are derived from IL-15 signaling.

Discussion
The poor healing ability of aged skin is a major medical issue (2), but we lack therapies for this debilitating problem. Wound closure defects in the elderly are especially challenging to manage because this population frequently undergoes surgery, has a greater risk of open skin lesions (i.e. bedsores), and often suffers from disorders that impede tissue repair such as vascular disease or diabetes. The identification of new therapies to treat poorly healing wounds has stalled in recent decades and only 3 healing therapies have been approved for clinical use in the last 24 years (40). One reason for this lack of progress is that the underlying cellular defects responsible for aging are not well understood and potential therapies must overcome impairments in signal transduction in aged cells (5). To begin addressing these problems we have used a highly powered bioinformatics approach to associate several known pathways of keratinocyte aging, such as JAK/STAT signaling and cellular senescence, to an area not previously connected to healing defects: bioenergetic dysfunction. We establish that daily pulses of IL-15 via injection or exercise, which boost mitochondrial abundance (20), restores youthful keratinocyte STAT3 signaling, reduces markers of growth arrest and improves healing in aged skin. We observe that aging impedes wound healing by about 2 days, which is in line with the delay in healing observed in old versus young humans (41). Moreover, our finding that IL-15 can speed wound closure by ~2 days is similar to the improvement seen following becaplermin (PDGF) treatment in ischemic and hyperglycemic mouse models (42). However, in more severe clinical healing disorders, such as lower extremity ulcers, becaplermin treatment is able to reduce wound healing time by multiple weeks (43). Thus, IL-15 may have more substantial healing benefits if it is used to treat severe skin ulcers.
We opted to use intravenous tail vein injections to mirror the post-exercise secretions of IL-15 into the blood based around our prior work (20). We understand that a more practical delivery method is therapeutically desirable, however, drug delivery to the wound bed directly has its own confounding issues. First, healing wounds have a very proteolytic environment, so the bioavailability of rmIL-15 could be severely impacted and potentially in a different manner in young and old wounds. Second, most carrier formulations (i.e. ointments) make it difficult to visualize the leading edge of the wound for closure measurements. It is likely that a water-soluble small molecule would get around these issues, but this requires further investigation. Nevertheless, our data demonstrate that intravenous IL-15 administration or exercise reverse fundamental aspects of skin aging and suggest that targeting these downstream molecular pathways can guide future efforts to improve poor skin healing.
Tissues directly involved in exercise such as skeletal muscle or the heart positively adapt over time to improve their physiological function. This adaptation occurs in response to a wide variety of cellular cues, including mechanical stress, calcium fluctuations and other complex signals that would be difficult to recapitulate with drug injections (44). Thus, it is notable that a single factor such as IL-15 can enhance wound closure in aged skin and is the responsible exercise signaling factor for healing. We believe that IL-15 is effective for several reasons. First, healing is complex and treatments likely need to impact multiple facets of tissue Interleukin-15 rescues impaired wound healing in aged skin remodeling and regrowth to be effective. We find that the actions of IL-15 are mirrored in keratinocytes and dermal fibroblasts in vitro and in vivo, suggesting direct, multi-potent benefits of IL-15 treatment. Second, we find that IL-15 is able to reverse basic aspects of cellular aging, which likely permits enhanced tissue re-growth after injury. IL-15 restored epidermal mitochondrial (Ndufs2, Cox7b) and cytokine mRNA (Il6, Cxcl1, Cxcl2) as well as mitochondrial enzyme activity in aged epidermis. However, low dose IL-15 does not increase nuclear pSTAT3 Y705 , suggesting a resensitization of keratinocyte recruitment or activation that became lost during aging. This idea fits well with prior research showing that circulating IL-15 is reduced in aged mice (23) and IL-15 mRNA expression is lower in old vs. young keratinocytes at the leading edge of a healing wound (9). However, whether the loss of IL-15 signaling itself or other factors initially trigger skin aging and defects in tissue healing remains unknown.
Our understanding of how mitochondrial signaling and energy generation impact wound healing in the skin is limited, particularly in the context of intrinsic, chronological aging. However, there is some evidence for the involvement of mitochondrial deterioration in the other form of skin aging, UV-radiation induced photoaging. This type of aging is primarily due to exposure to UV-A and UV-B radiation in sunlight, which causes both nuclear and mitochondrial DNA damage and also increases the release of mitochondrial free radicals (8). Human photoaged skin exhibits a high incidence of large-scale mtDNA deletions including two characteristic 3895-base pair and 4977-base pair common deletions which are significantly elevated compared to non-photoaged skin from the same subject (45). These mtDNA deletions target portions of Complex I and IV, which impairs mitochondrial ATP generation (46).
The role of mitochondrial function in chronological aging-associated disease has been extensively investigated in other organs, however its impact in the skin is relatively understudied. The most direct example of mitochondrial dysfunction in driving tissue aging is the progeroid aging phenotype of polymerase gamma mutator mice (PolG mice) that have error prone replication of their mitochondrial DNA (mtDNA) (47). These animals systemically accumulate excessive mtDNA mutations, have reduced mitochondrial energy generation, and experience accelerated aging throughout the body including the skin. Notably, treatments that boost cellular metabolism such as endurance exercise (48) or the drug bezafibrate (49), rescues the epidermal, dermal and sub-dermal adipose layer thinning of the skin observed in PolG mice. The striking reversal of progeroid aging in PolG mice by exercise training (48) was the initial catalyst for our previous investigation regarding exercise-induced circulating factors that impact skin aging (20). Based upon these initial observations in uninjured skin we now show that circulating IL-15 is part of a critical signaling circuit that restores aged skin repair through increased re-epithelialization and an enhancement in epidermal stem cell growth. Our findings add to the growing body of literature detailing how physical activity improves health and combats disease such as studies showing that exercise augments the post-injury repair of cardiac muscle (50) and peripheral nerves (51).
Unlike other stem cell populations, aging does not reduce the total number of epidermal stem cells, but does appear to stall their growth potential in vitro as well as their proliferation after wounding (9, 10). While slow healing in aged skin has been associated with altered keratinocyte STAT3 signaling (9), there is little known about the role of cellular senescence. The delay in keratinocyte proliferation after wounding in healthy, aged mouse skin implies transient pausing of the cell cycle rather than extensive permanent growth arrest. Our data corroborate this idea and show that aging alters the alarmin and transient senescence marker HMGB1 but does not impact permanent senescence markers such as p16 ink4a . Additionally, we show that both HMGB1 and STAT3 dynamically respond to aging, exercise and IL-15 treatment in the skin. This is not entirely unexpected as both factors activate pro-inflammatory processes after stress (38,52,53). Notably, a similar loss of HMGB1 has been associated with poor wound healing in diabetic skin and exogenous HMGB1 treatment restored diabetic healing defects (54), however this prior study did not examine skin STAT3 signaling. Thus, changes in STAT3 may drive the loss of HMGB1 from the nucleus or vice versa to stall the cell cycle and delay re-epithelialization, but more work is needed regarding the interaction of these pathways in the context of wound re-epithelialization.

Interleukin-15 rescues impaired wound healing in aged skin
While the biology of STAT3 has been extensively studied in disease, until now, all of our knowledge regarding STAT3 in aging skin has focused on its Y705 phosphorylation site, which is well-established to increase in the aged epidermis (9,55). However, no studies have previously characterized aging-induced changes to the STAT3 S727 site. Classical pro-inflammatory cytokines such as IL-6 or GM-CSF increase STAT3 Y705 phosphorylation as well as induce senescence (56), and this class of cytokines is increased in the blood during aging as a consequence of low-grade inflammation (57). However, it has not been clear how one might physiologically reverse these cellular changes. Mechanistically, increasing pSTAT3 S727 is one of the only known ways to promote dephosphorylation of pSTAT3 Y705 via upregulation of the phosphatase TC45 (58). At the same time pSTAT3 S727 also directs the mitochondrial localization of STAT3 and the resulting enhancements in bioenergetic function (59). Since we show increased cellular pSTAT3 S727 after IL-15 treatment and exercise, this seems the most likely explanation for the corresponding reduction in pSTAT3 Y705 . Thus, there appears to be a counterbalancing role for pSTAT3 S727 in both cellular metabolism and cell growth during wound repair through moderation of pSTAT3 Y705mediated signaling.
These studies expand our understanding of the health benefits of exercise in reducing the risk of age-related disease. Our findings suggest that other skin or epithelial disorders may be improved by regular exercise or conversely worsened by a sedentary lifestyle. Some indications of this connection have been established as vigorous physical activity adherence is associated with a lower risk of the inflammatory skin disorder psoriasis (60). Additionally, our findings have high translational potential as they are rooted in the physiologic exercise signaling observed in mice and humans and exercise is already known to improve healing in elderly subjects. Finally, higher doses of recombinant human IL-15 are already FDA approved for cancer immunotherapy, potentially reducing regulatory hurdles in applying a lower dose as a treatment for poorly healing wounds. Overall, these studies demonstrate that IL-15 therapy is a safe and efficacious means of reversing wound healing defects in aging skin.
Keratinocyte explant culture was performed by shaving and depilating telogen dorsal back skin from C57Bl/6 mice (3-5 months old for young and 23-24 months old for old) using clippers and NAIR before skin excision. The excised tissue was then incubated for 3 minutes in a solution of povidone-iodine, then twice in a solution of 70% ethanol for a total of 3 minutes, and finally rinsed in sterile PBS. The skin was then pinned epidermis side down and scraped with a beveled scalpel blade until all subdermal adipose had been removed and the dermis was abraded and translucent. A 4 mm biopsy punch was used to excise individual explants, which were plated dermis side down into a 24-well tissue culture plate (BioLite, Thermo) with a 1 hour drying period to promote explant adhesion. Following the drying period, 400 µL of media (DMEM supplemented with 15% FBS, 1% glutamine, 1% penicillin, and 1% streptomycin) was added in order to submerge the explant. Explants were cultured at 37°C with 5% CO2 and Interleukin-15 rescues impaired wound healing in aged skin were cultured for 48 hours after treatments before imaging cell outgrowth.
All fibroblast and keratinocyte explant experiments were performed in triplicate technical replicates to generate data for a single biological replicate (i.e. n = 1 for each unique mouse donor). Primary mouse fibroblast identity was determined by morphology. Primary mouse keratinocyte identity was confirmed by staining for Keratin 14, Keratin 10, or Loricrin as well as negative staining for Vimentin.

Tissue Lysis and Western Blotting
Freshly dissected epidermis was separated from dermis using dispase (Roche) enzymatic digestion for approximately 1.5-2 hours at 37°C. Isolated epidermis was then homogenized in lysis buffer (50 mM HEPES, 150 mM NaCl, 100 mM NaF, 10 mM Na pyrophosphate, 5 mM EDTA, 250 mM sucrose, 1 mM DTT, 1% Triton-X, and Pierce Protease and Phosphatase Inhibitor solution (Thermo Scientific)). Primary dermal fibroblast and epidermal lysates were homogenized by freeze thaw 3 times in liquid nitrogen followed by passage through a 27-gauge needle to ensure complete lysis. The protein concentration of epidermal lysates was assessed using a BCA kit (Pierce). For the fibroblast analysis of STAT3, 20 µg of protein per sample was prepared for SDS-PAGE by boiling in reducing 4x SDS-sample buffer (Boston BioProducts) for 5 minutes at 98°C. Samples were loaded on a Mini-PROTEAN 4-15% TGX Stain-Free Tris/tricine pre-cast gel (Bio-Rad) and run with Laemmli SDS-Page running buffer (Bio-Rad) for approximately 90 minutes at 100 V. Gels were semi-dry transferred onto PVDF membranes using the Trans-Blot Turbo system (Bio-Rad) and all imaging was done with a ChemiDoc MP (Bio-Rad). Membranes were blocked and antibodies were diluted in Tris-buffered saline with 0.1% Tween-20 and 5% bovine serum albumin fraction V (Bio-Rad). Primary antibodies and their usage were as follows: b-actin (rabbit, Cell Signaling #4967, 1:1000) total STAT3 (rabbit, Cell Signaling #30835S, 1:1000), pSTAT3 Y705 (rabbit, Cell Signaling #9145, 1:1000), pSTAT3 S727 (rabbit, Abcam #AB30647, 1:1000). Membranes were stripped and re-probed using Restore Western Blot Stripping Buffer (Thermo) between pSTAT3 S727 , pSTAT3 Y705 , and total STAT3. Secondary antibody used was anti-rabbit-IgG-HRP (H+L) (goat, Bio-Rad, 1:1000) with Clarity Western ECL Substrate (Bio-Rad) solution. Band densitometry was measured in Image Lab (Bio-Rad) to quantify relative phosphorylation normalized to b-actin. Western blotting was performed in triplicate with three biologically unique primary fibroblast cultures.

Animals
All animals used for the chronological aging experiments were young, adult (3-5 months old) or old (23-24 months old) female C57BL/6 mice. Female mice were used because the thickness of the skin and the rate of wound healing differs between male and female mice, and because female mice are less prone to fighting and barbering when group housed. Mice were randomly allocated to treatment groups prior to the interventions. During the study, mice were group housed except after wounding surgery when they were individually housed to minimize disturbance to the wound site. All mice were maintained in an AAALAC approved animal facility and all procedures were approved by the Northeastern University IACUC according to NIH standards.

In vivo Treatments
Old exercising mice that received IL-15 neutralizing antibody or IgG control antibody were exercised daily on a 15-degree incline using a 6lane rodent treadmill (Columbus Instruments) at 16 m/min pace after initial 1 week ramp up from 10 m/min. Exercising mice underwent training for 33 days plus 5 days post-wound surgery but did not exercise on the day of surgery. Post-surgery the exercising mice ran at 8 m/min to minimize wound aggravation. Recombinant mouse IL-15 (rmIL-15) (R&D Systems #447-ML/CF) was administered at 25 ng/kg body weight by tail vein injection. IL-15 neutralizing antibody (R&D Systems #AF447) or normal goat IgG control antibody (R&D Systems #AB-108-C) was provided at 8,750 ng/kg body weight, via tail vein injection daily prior to exercise. Compounds were administered between 2-4 hours prior to treadmill exercise sessions, once daily for 38 consecutive days.

Wounding Surgery
Prior to surgery, all mice were anesthetized with isoflurane, provided subcutaneous carprofen for pain relief (5 mg/kg) (TCL) and hair was shaved Interleukin-15 rescues impaired wound healing in aged skin with clippers. The skin was then swabbed with povidone-iodine (APLICARE) and ethanol. Thereafter, two 6 mm diameter dorsal fullthickness wounds were excised using a dermal biopsy punch (Medline) and this sample was used for intact skin histology and mRNA analysis. Wounds were made in telogen skin in all experiments. Wounds were splinted using a 10 mm inner diameter stainless steel ring (Grainger) secured with 4-6 monofilament nylon 5-0 sutures (Oasis) and cyanoacrylate adhesive (Krazy Glue) as described (32). A clean Tegaderm (3M) wound dressing was loosely applied around the circumference of the mice over the wounds to cover the splinted wounds. Carprofen (5 mg/kg) was administered subcutaneously for the following 3 post-operative days. Starting on the day of surgery, wounds were photographed daily using a DMC-FZ70 camera (Panasonic) when the Tegaderm dressing was changed until 9 days post-surgery. Wound site photography lighting conditions and imaging parameters on the camera were kept identical between days. Proliferating cells were labeled using 100 mg/kg BrdU (Sigma) that was administered intraperitoneally 24 hours before wounding or 100 mg/kg EdU (Cayman Chemical) that was administered intraperitoneally 24 hours before sacrifice. Wound surgery occurred on day 33.
For immunofluorescence staining of skin explants, samples were grown on cover glass and surrounded by a sterile 0.5 mm press-to-seal silicone gasket (Invitrogen) to reduce staining volume. After 48 hours of skin explant outgrowth, cells were fixed for 1 hour at room temperature in 10% neutral buffered formalin followed by 10 minutes of cell permeabilization using 0.1% Triton X-100 in PBS. Cells were then blocked for 30 minutes in 5% normal goat serum in PBS followed by a 2 hour room temperature incubation with primary antibody and 30 minute room temperature incubation with secondary antibody as described above. Hard set mounting media with DAPI was applied and cured prior to imaging.
Hematoxylin and Eosin staining was performed on 4 µm de-waxed sections using standard histology protocols using Gill's Hematoxylin No.1 and Eosin Y (Sigma). Colorimetric stained slides were mounted with Permount (Fisher). All brightfield imaging was done using an inverted EVOS XL (Life Technologies) microscope. All fluorescence imaging was performed on an Revolve R4 (Echo Labs) microscope equipped with Olympus UPlanFL 10x/0.30, and UPlanFL 20x/0.50 air objectives using DAPI, FITC, and Cy5 fluorescence channels.

Image Analysis
All image processing and analysis was conducted using ImageJ FIJI (version 2.0.0-rc-69/1.52i, NIH), Photoshop CC 2019 (version 20.0.4, Adobe), and Illustrator CC 2019 (version 23.0.3, Adobe). Daily images of the wound site were used to quantify wound closure over time. Images were processed to normalize camera distance and angle using the size of the metal splint as a standardization factor. Wound measurements were obtained by tracing the area of the visible wound margin on each day relative to the splint area. Wound closure over time was normalized as a percentage of the original day 0 wound. The rate of wound healing was analyzed by fitting wound closure on days 0-9 on a linear best fit line.
Immunofluorescence images first underwent thresholding over entire image set of all samples from the experiment followed by manual counting of positive cells for at least 3 unique fields from each sample. All epidermal analyses only counted interfollicular epidermal cells and excluded hair follicles. Dermal fibroblasts were identified by vimentin staining. pSTAT3 Y705 positive cells were those that had overlap with their DAPI nuclear signal. pSTAT3 S727 positive cells were those that had DAPI-associated cytoplasmic signal. HMGB1 positive cells were those that had overlapping signal with their nuclear DAPI signal. Lamin B1 positive cells were those that exhibited perinuclear signal relative to their DAPI signal. Vimentin positive cells were those that had DAPIassociated periplasmic signal. Keratin 14, Keratin 10, and Loricrin positive cells were those that had DAPI-associated cytoplasmic signal.
Average thickness of the epidermal proliferative hub was measured over 300-micron span from wound margin, guided by EdU staining along with reference to the subdermal adipose layer and serially sectioned H&E samples (34). Composite images for thickness analysis were tiled by manually registering images of 10x fields with approximately 20% overlapping content using Photoshop (Adobe), and individual thickness measurements were taken using ImageJ FIJI (NIH).
Epidermal thickness measurements were performed on H&E stained non-wounded skin cross-sections as previously described (20) by measuring the orthogonal distance from the basement membrane to the end of the spinous layer several times over the interfollicular epidermis.
Keratinocyte outgrowth was imaged in a circular overlapping clockwise pattern from 12 o'clock 24 hours after plating and areas of outgrowth were imaged at 48 hours after plating. Outgrowth at 48 hours post-explanting were measured, and the maximal area outgrowth field for each discrete outgrowth on each explant was recorded by tracing its outline and calculating the enclosed area using Illustrator (Adobe). The average of the 3 highest of these maximal outgrowth area fields was taken and used as representative growth from the explant. For each data point, this process was done 2 times on 2 different pieces of explant skin from the same mouse.

Gene Expression
The epidermis was separated from dermis using dispase (Roche) enzymatic digestion for 2 hours at 37 °C. After digestion, the epidermis was scraped using a scalpel blade and immediately frozen at -80 °C until later analysis. RNA was extracted by homogenizing separated epidermis in TRIzol (Invitrogen) reagent using a bead mill apparatus (MPBio 5G) followed by column purification (Zymo Direct-zol RNA mini prep), on column DNase treatment and subsequent elution. 300 ng of total RNA was then reversed transcribed to cDNA (ABI HC cDNA synthesis kit), diluted 1:30 in ultrapure water and mRNA expression was assessed using qPCR with TaqMan or SYBR chemistry.

Flow Cytometry
Whole blood was collected on day 31 of the rmIL-15 or exercise experiments from the submandibular vein. Blood was collected into K2-EDTA MiniCollectTubes (Greiner Bio-One) and RBCs were lysed using PharmLyse (BD Biosciences). Lymphocytes were then pelleted and incubated with FITC-CD3 (BD Pharminigen #553061) and APC-CD8 (BD Pharminigen #561093) antibodies and analyzed in conjunction with unlabeled and compensation controls on a BD FACSAria instrument. Data analysis was done using FlowJo software for MacOSX.

Microarray Analysis
Microarray data from GSE84511 (28) was obtained from GEO using the Bioconductor suite of tools in R (version 3.5.1) on MacOSX. Differential Interleukin-15 rescues impaired wound healing in aged skin gene expression between young and old epidermal stem cells were obtained using the LIMMA package (version 3.36.2) to contrast all aged animal GSM entries against all young animal GSM entries using group-means parameterization, empirical Bayes statistics for differential expression, and Benjamini-Hochberg false discovery rate adjustment (62). Gene set enrichment analysis using the FGSEA package (version 1.6.0) (default pathways, p. adj < 0.05, 10000 permutations) and bi-directional KEGG pathway analysis using the GAGE package (version 2.30.0) was performed on differentially expressed genes using mouse (mmu) gene sets for bioprocesses, molecular factors, and cellular compartments with wiring diagrams generated using the Pathview package (version 1.20.0) for GAGE (29,30).

Cytochrome C Oxidase Activity
Cytochrome C oxidase activity was measured as the decrease of absorbance due to oxidation of cytochrome c in epidermal lysates as described previously (20). Briefly, sodium dithionite (Fisher, S310-100) was used to reduce Cytochrome C (Sigma, C-2506) in assay buffer (50 mM potassium phosphate buffer with 1 mM EDTA, pH 7.4). Reduced cytochrome c was added to 20 µL of sample lysate per well in a 96-well plate and the absorbance at 550 nm was monitored over 90 seconds as a fitted linear slope (R 2 > 0.95). Samples were measured in duplicate and absolute enzyme activity was normalized to protein concentration: "[ ] × >(?@A* 2 1000 ⁄

Statistical Analysis
Graphpad Prism (version 7.0d) was used to calculate statistics. Data were log transformed or square root transformed prior to statistical analysis. Significance was set as p < 0.05. Comparisons of two groups were performed using a student's t test.
Comparisons of more than 2 treatment groups were analyzed using a one-way or two-way analysis of variance (ANOVA) as appropriate for the number of treatment factors. If statistical significance was achieved via ANOVA, specific differences between groups were identified using Holm-Sidak or Tukey's post-hoc analysis.

Interleukin-15 rescues impaired wound healing in aged skin
Interleukin-15 rescues impaired wound healing in aged skin showing sets of genes either up or down with age associated with specific gene ontology terms. Gene ontologies are ordered by their normalized enrichment score that signifies how uniform the gene expression change is across all member genes of a set (Gene Set Rank). (C) Candidate differentially expressed genes in young versus old mouse epidermal stem cells. All analysis in A-C was generated using raw data from GSE84511. Significance was determined by a Benjamini-Hochberg False Discovery Rate adjusted p < 0.05.