Injury prevents Ras mutant cell expansion in mosaic skin

Nature volume 619, pages 167–175 (2023)Cite this article 11k Accesses 154 Altmetric Metrics details Healthy skin is a mosaic of wild-type and mutant clones1,2. Although injury can cooperate with

Nature volume 619, pages 167–175 (2023)Cite this article

11k Accesses

154 Altmetric

Metrics details

Healthy skin is a mosaic of wild-type and mutant clones1,2. Although injury can cooperate with mutated Ras family proteins to promote tumorigenesis3,4,5,6,7,8,9,10,11,12, the consequences in genetically mosaic skin are unknown. Here we show that after injury, wild-type cells suppress aberrant growth induced by oncogenic Ras. HrasG12V/+ and KrasG12D/+ cells outcompete wild-type cells in uninjured, mosaic tissue but their expansion is prevented after injury owing to an increase in the fraction of proliferating wild-type cells. Mechanistically, we show that, unlike HrasG12V/+ cells, wild-type cells respond to autocrine and paracrine secretion of EGFR ligands, and this differential activation of the EGFR pathway explains the competitive switch during injury repair. Inhibition of EGFR signalling via drug or genetic approaches diminishes the proportion of dividing wild-type cells after injury, leading to the expansion of HrasG12V/+ cells. Increased proliferation of wild-type cells via constitutive loss of the cell cycle inhibitor p21 counteracts the expansion of HrasG12V/+ cells even in the absence of injury. Thus, injury has a role in switching the competitive balance between oncogenic and wild-type cells in genetically mosaic skin.

Throughout our lifetimes, we acquire mutations in our skin, owing to its constant exposure to environmental insults. As a result, phenotypically normal skin contains a mosaic of epithelial stem cells with somatic mutations, including in genes that are associated with cancer development, such as the GTPase Ras family1,2. Constitutive activation of Ras oncogenes has been identified as the initial genetic event in 3–30% of human cutaneous squamous cell carcinomas13,14,15 (cSCCs) and in experimentally induced cSCCs in mice16,17. In mouse models with mosaic epithelial expression of the constitutive active form of Hras (HrasG12V/+), mutant cells outcompete wild-type cells and expand in the uninjured skin epidermis18,19,20. Although activated-Hras mutant cells are tolerated within otherwise wild-type and uninjured skin epithelium18,19,20, injury has been shown to cooperate with oncogenic mutations to trigger tumorigenesis in various mouse models3,4,5,6,7,8,9,10,11,12. We hypothesized that the expansion of HrasG12V/+ cells in the epidermis could represent a vulnerability upon injury; HrasG12V/+ cells could futher expand and lead to tumours. For instance, the hyperproliferative environment generated during injury repair may further stimulate the proliferative behaviour of mutant cells and break the tolerance of the tissue. Here we investigated how injury affects the oncogenic potential of HrasG12V/+ within a genetically mosaic and phenotypically relevant context.

The stratified skin epidermis is uniquely accessible for direct observation, which enables the visualization of the emergence of aberrant growth at single-cell resolution. The basal layer contains epidermal stem cells, which can self-renew to generate more basal cells or differentiate and delaminate upwards to replace outer, barrier-forming cells21,22 (Extended Data Fig. 1a). We hypothesized that injury repair would cooperate with the constitutive activation of the Hras oncogene (HrasG12V) to promote tumorigenesis in phenotypically normal, genetically mosaic skin. To test this hypothesis, we generated mice in which we could induce and follow populations of HrasG12V/+ mutant cells within wild-type epithelium (Krt14-CreER; flox and replace (FR)-HrasG12V/+; Lox-STOP-Lox (LSL)-tdTomato; Krt14-H2B–GFP; Methods). In these mice, tamoxifen treatment activates Cre in keratin 14-expressing basal cells and, in turn, induces the co-expression of HrasG12V/+ from its endogenous promoter and a cytoplasmic fluorescent tdTomato reporter that provides an approximation of mutant cells. Moreover, these mice also express histone H2B-GFP in basal cells, which persists throughout differentiation23, enabling the visualization of all basal stem cells and their progeny (Fig. 1a). We treated mice with tamoxifen at 3 weeks of age and, three days later, introduced a full-thickness injury down to the cartilage (4 mm diameter punch biopsy) in one ear. We used two doses of tamoxifen to drive HrasG12V/+ expression in either approximately 99% of basal cells (HrasG12V/+ max) to recapitulate previous studies of homogeneous models3,4,5,24,25, or in approximately 65% of basal cells (HrasG12V/+ mosaic) to mimic genetically mosaic skin (Fig. 1b). As a control, we also engineered Krt14-CreER; LSL-tdTomato; Krt14-H2B–GFP mice and treated them similarly to drive tdTomato expression in approximately 65% of wild-type basal cells (wild-type mosaic) (Fig. 1b). By longitudinally imaging the same regions of the skin epidermis and always comparing HrasG12V/+ mosaic models to wild-type mosaic models, we were able to control for any potential leakiness of the CreER system and study the consequences of HrasG12V/+ expression on the tissue overall and cell behaviours (see Methods).

a, Left, cartoon depicting the 4-mm full-thickness wound (W) on a mouse ear and the imaged area. Centre, top-down (x–y; top) and transverse (x–z; bottom) views of a two-photon image of the skin epithelium at 14 days PWI from a Krt14-CreER; LSL-tdTomato; Krt14-H2B–GFP mouse (asterisks, hair canals; dashed lines, basement membrane in x–z view and wound edge in x–y view). Centre middle, cartoon schematic of wild-type (green) and recombined cells (red) after tamoxifen injection around the injury. Right, magnification of top-down (x–y) and transverse (x–z) views of the skin epithelium show epithelial cell nuclei (Krt14-H2B–GFP) in green and recombined cells expressing tdTomato in red (dashed lines mark the transversal section in the x–y view and the basement membrane in the x–z view). b, The initial tdTomato+ area 3 days PWI (n = 4 wild-type mosaic and HrasG12V/+ mosaic mice and n = 5 HrasG12V/+ max mice). At least three independent areas of approximately 300 μm2 were analysed for each mouse (Methods). Data are mean ± s.d. c, Heat maps of the top-down (x–y) view of representative two-photon images adjacent to the injury at 14 days PWI (dashed lines highlight the wound edge). Colour represents the thickness of the epithelium and identifies the presence of aberrant growth. d, The thickness of the epithelium at 14 days PWI around the wound. Solid lines represent means and dashed lines show s.d. n = 4 wild-type mosaic and HrasG12V/+ mosaic mice and n = 5 HrasG12V/+ max mice.

Source Data

We monitored the injured epithelium over time by combining deep tissue imaging with quantitative analyses using IMARIS and MatLab software, which enabled us to evaluate tissue thickness with intensity heat maps (Fig. 1c,d, Extended Data Fig. 1f, Supplementary Videos 1–4 and Methods). Notably, at 14 days post-wound induction (PWI), our HrasG12V/+ mosaic models, which are on an outbred genetic background, did not develop the aberrant growth and thick epithelium observed in HrasG12V/+ max models (Fig. 1c,d, Extended Data Fig. 1f, Supplementary Videos 1–4 and Methods). Histopathological analysis further showed that the skin epithelium around the repaired injury in the HrasG12V/+ mosaic model was normal, despite the high burden of the HrasG12V/+ mutation (Extended Data Fig. 1b,c,e,g,i). By contrast, abnormal growth formed rapidly within two weeks after injury induction in HrasG12V/+ max, as expected (Extended Data Fig. 1d,e,h,i).

Collectively, these data show that HrasG12V/+ cells break homeostatic tissue architecture during injury repair only when nearly all cells in the basal stem cell layer have this genotype.

Having found that injury repair does not trigger aberrant growth in HrasG12V/+ mosaic tissue, we next investigated how HrasG12V/+ and wild-type cells within mosaic epithelia respond to injury. Previous studies had shown that embryonically induced HrasG12V/+ basal stem cells integrate and expand in the skin epidermis, eventually outcompeting wild-type cells18,19,20. We had hypothesized that this proliferative advantage of HrasG12V/+ cells would be amplified during injury repair, which has a higher proliferative demand than uninjured skin. To test our hypothesis, we revisited the same wild-type and HrasG12V/+ cells in the skin epidermis of live mice for one month, with or without injury repair (Fig. 2a and Extended Data Fig. 2a). Epithelial cells begin to contribute to re-epithelialization approximately three days after injury26,27. We therefore started our analysis three days PWI (six days after tamoxifen-induced mosaicism) (Fig. 2a). We drew boundaries between GFP+tdTomato+ and GFP+tdTomato− regions and represented the tdTomato+ areas as a percentage of the total area (Fig. 2b–e).

a, Schematic of the experimental design using the injury repair model. b, Representative two-photon revisit images of the basal stem cell layer of the epidermis. White lines highlight the boundaries between tdTomato+ and tdTomato− populations. Epithelial nuclei are in green (Krt14-H2B–GFP) and recombined cells are in red (LSL-tdTomato). c, The increase in tdTomato+ area over time in the uninjured condition. n = 3 mice. d, Representative two-photon revisit images of the basal stem cell layer of the epidermis during injury repair. e, The increase in tdTomato+ area over time in the injured condition. n = 4 mice. Unpaired, two-tailed t-test comparing wild-type with HrasG12V/+ mice at different time points in uninjured and injured conditions. P values are shown. At least three independent areas of approximately 300 µm2 were analysed for each mouse (Methods). Data are mean ± s.d. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001; NS, not significant.

Source Data

We found that the HrasG12V/+-tdTomato+ population expanded substantially in uninjured HrasG12V/+ mosaic epithelium, to a greater extent than the tdTomato+ population in uninjured wild-type mosaic epithelium. In uninjured HrasG12V/+ mosaic mice, HrasG12V/+ cells outcompeted wild-type cells and increased their occupancy of the basal stem cell layer by approximately 20% after one month, consistent with previous work using constitutive Cre expression18 (Fig. 2b,c and Extended Data Fig. 2b). By contrast, one month after injury, HrasG12V/+ cells did not expand over wild-type cells and the proportion of tdTomato+ cells was similar between wild-type and HrasG12V/+ mosaic models (Figs. 1b and 2d,e). Notably, these results were also observed in mosaicism with a lower HrasG12V/+ mutational burden, which more closely mimic physiological conditions (Extended Data Fig. 2c–f). This finding was unexpected and indicates that HrasG12V/+ cells did not outcompete wild-type cells and expand after injury repair in HrasG12V/+ mosaic mice.

Collectively, these data demonstrate that the injury repair process does not amplify but rather abrogates the competitive advantage that HrasG12V/+ cells have over wild-type cells in the absence of injury.

The suppressed expansion of HrasG12V/+ cells after injury prompted us to investigate how this process affects different cellular behaviours of mutant and wild-type cells. Specifically, we first examined the number of dividing cells, given that Ras family proteins are key regulators of epithelial cell proliferation in the skin epithelium. Indeed, epithelial stem cells in vitro and in vivo do not proliferate upon ablation of all Ras isoforms28. To score the number of dividing cells, we immunostained cells for the mitotic marker phosphorylated histone H3 (p-histone H3). We observed an increase in the amount of epithelial cell divisions accompanying efficient wound repair at 3 days PWI in the wild-type mosaic model (Fig. 3a,b). In sharp contrast, although we observed an increase in mitotic events in wild-type cells in the HrasG12V/+ mosaic model, the proliferation of HrasG12V/+ cells was unaltered during repair (Fig. 3b). These findings were corroborated by measuring mitotic figures (Extended Data Fig. 3a–c). Therefore, wild-type cells have an unexpected and selective competitive advantage over HrasG12V/+ cells in the acute phase of injury repair.

a, Representative two-photon images of the epidermal preparation, immunostained for p-histone H3 (highlighted by white dashed circle). b, Quantification of p-histone H3-positive (p-H3+) cells in tdTomato+ and tdTomato− populations in injured (I) and uninjured (U) skin. n = 4 mice. c,d, Quantification of mitotic cells in uninjured (c) and uninjured (d) ears. n = 4 mice. e, Representative two-photon images of the epidermal preparation, immunostained for keratin-10 (white). f, Quantification of keratin-10-positive (KRT10+) cells in tdTomato+ and tdTomato− populations in injured and uninjured skin. n = 3 mice. b–d,f, Paired or unpaired two-tailed t-test for comparison between tdTomato+ and tdTomato− populations in the same group or in different groups of mice. At least three independent areas of approximately 300 µm2 were analysed for each mouse (Methods). KRT10, keratin-10. g, Uniform manifold approximation and projection (UMAP) showing epidermal keratinocyte clusters from uninjured and injured conditions. h, Violin plots (left) showing gene expression in clusters from g, together with the gene expression superimposed on the UMAP (right). i, Cell classification into homeostatic and injury-responsive populations based on clustering and gene expression, overlaid on the UMAP. g–i, n = 12 mice. j, Quantification of wild-type and HrasG12V/+ cells as a proportion of tdTomato+ cells in basal non-committed and basal committed groups for homeostatic and injury-responsive cells. Data are averaged results for each biological replicate. One-sided t-test with Holm–Sidak correction. k, Violin plots showing Krt6a and Krt10 expression in cells grouped as in j. In h,k, internal box plots denote the 25th, 50th and 75th centiles, with whiskers depicting minima and maxima, excluding outliers that are beyond 1.5× the interquartile range. j,k, Homeostatic cells from n = 12 mice and injury-responsive cells from n = 6 mice. P values are shown. Data are mean ± s.d.

Source Data

To determine whether the increased number of wild-type cell divisions that we observed at 3 days PWI was sustained over time, we scored mitotic events over 4 weeks. In the wild-type mosaic model, the initial increase in the number of mitotic cells observed in tdTomato+ and tdTomato− wild-type cells returned to baseline by 7 days PWI (10 days after tamoxifen-induced mosaicism) and appeared similar to uninjured mice, as expected (Fig. 3c,d). The initial increase in divisions of wild-type cells observed in the injured HrasG12V/+ mosaic model also decreased over time, and eventually returned to a similar level as in neighbouring HrasG12V/+ cells, but remained higher than in wild-type cells following injury in the wild-type mosaic model (Fig. 3c,d). By contrast, the proliferative capacity of HrasG12V/+ cells in the HrasG12V/+ mosaic model was not substantially affected by injury repair at any of the time points analysed (Fig. 3c,d). Thus, the balanced mitotic events of wild-type and HrasG12V/+ cells sustained at later time points after injury would effectively continue to prevent the expansion of HrasG12V/+ cells in the HrasG12V/+ mosaic model.

To test whether injury repair leads to genotype-specific changes in other cellular behaviours, we monitored apoptosis, which is an established cell competition mechanism29, and inhibited by Ras signalling30. We examined cell death by scoring for either nuclear fragmentation events or expression of an apoptotic marker, active caspase-3. The overall frequency of apoptosis was low, and we did not observe significant differences in cell death events of wild-type or HrasG12V/+ cells in mice with or without injury at 6 days after tamoxifen-induced mosaicism and at later time points (Extended Data Fig. 3d–g). Similarly, we found no evidence of cell senescence (Extended Data Fig. 3h). Differentiation is another mechanism of cell loss that could influence competition between wild-type and HrasG12V/+ cells in the skin epidermis. We therefore comprehensively evaluated the number of differentiation events and the expression of early differentiation markers using both protein (immunostaining) and transcriptional (single-cell RNA sequencing (scRNA-seq)) analyses in uninjured and injured (3 days PWI) settings of wild-type mosaic and HrasG12V/+ mosaic models (Fig. 3e–k and Extended Data Figs. 3i and 4a–h). The scRNA-seq analysis revealed an increased fraction of epidermal basal cells committing towards differentiation upon injury, which we inferred from the ratio of basal cell states (stem cells versus differentiation-committed22) and corresponding changes in differentiation and stemness transcripts (Fig. 3j,k and Methods). The proliferation signature that emerged in Gene Ontology terms for ‘biological process’ in uninjured settings did not appear in injured contexts when comparing HrasG12V/+ mosaic to wild-type mosaic models, further corroborating the live-imaging observations (Extended Data Fig. 4i,j). To determine whether the increased differentiation depended on the genotype (wild-type versus HrasG12V/+ in HrasG12V/+ mosaic mice), we quantified the expression of the early differentiation marker keratin-10 in wild-type and HrasG12V/+ cells by immunostaining. We observed a similar increase of differentiating cells for both genotypes upon injury compared to the uninjured condition, in contrast to the selective increase in the fraction of dividing wild-type cells (Fig. 3b,f). At 7 days PWI, we showed that the number of differentiating wild-type and HrasG12V/+ cells in the HrasG12V/+ mosaic model was similar, which in addition to their equal mitotic capacities, supports the maintenance of the mutant and wild-type clone sizes over time after injury (Fig. 3d and Extended Data Fig. 3i).

Collectively, these results show that injury increases the number of divisions of wild-type cells only, and therefore suppresses HrasG12V/+ cell expansion in the HrasG12V/+ mosaic model in the acute phase of injury repair.

Next, we investigated whether the injury repair process in mosaic skin could effectively suppress KrasG12D, a more aggressive Ras family oncogene. Mice with homogeneous KrasG12D/+ activation in the skin epidermis rapidly develop oncogenic growth in areas of constant abrasion31,32 (data not shown). Kras is one of the most frequently mutated oncogenes in human cancer, and is broadly activated across epithelial cancers, including cSCCs33. We generated mice in which we could induce and follow KrasG12D/+ cells within wild-type epithelium (Krt14-CreER; LSL-KrasG12D/+; LSL-tdTomato; Krt14-H2B–GFP; Methods). Similar to HrasG12V/+ tdTomato+ cells in the HrasG12V/+ mosaic model, KrasG12D/+ tdTomato+ cells expanded in the uninjured KrasG12D/+ mosaic mouse skin but did not expand after injury (Fig. 4a,b and Extended Data Fig. 5a). We scored mitotic events before and after injury induction, and again observed a selective increase in the amount of dividing wild-type cells but not of mutant cells in the KrasG12D/+ mosaic skin, similar to our observations in the HrasG12V/+ mosaic model (Fig. 4c).

a, Representative two-photon revisit images of the epidermal basal stem cell layer. White lines mark borders between mutant and wild-type cells. b, The increase in tdTomato+ area in uninjured and injured conditions following induction with tamoxifen (uninjured KrasG12D/+ mosaic, n = 4 mice; injured KrasG12D/+ mosaic, n = 5 mice). c, Quantification of mitotic cells in tdTomato+ and tdTomato− areas. Uninjured KrasG12D/+ mosaic, n = 4 mice; injured KrasG12D/+ mosaic, n = 5 mice. d, Heat maps of the top-down (x–y) view of representative two-photon images around the injury at 14 days PWI. Colour represents the thickness of the epithelium. e, Average epithelial thickness at 14 days PWI around the wound in wild-type (n = 3 mice), KrasG12D/+ mosaic (n = 4 mice) and KrasG12D/+ max (n = 3 mice). Solid lines represent means and dashed lines show s.d. b,c, Paired two-tailed t-test comparing tdTomato+ and tdTomato− populations in the same group of mice. Unpaired two-tailed t-test comparing tdTomato+ and tdTomato− populations in different groups of mice and KrasG12D/+ mutant mice in uninjured and injured conditions at different time points. P values are shown. At least three independent areas of approximately 300 µm2 were analysed for each mouse (Methods). Data are mean ± s.d.

Source Data

To monitor phenotypes at the tissue level, we applied two-photon microscopy with quantitative analyses of epidermal thickness represented by intensity heat maps. Despite the high burden of the KrasG12D/+ mutation (approximately 65% of recombined cells), the skin epithelium of KrasG12D/+ mosaic mice remained similar to wild-type mosaic models after injury (Figs. 1c and 4d,e and Extended Data Fig. 5a,b,d). By contrast, the KrasG12D/+ max model, in which nearly all basal cells had the KrasG12D/+ genotype, displayed rapid abnormal growth within the first two weeks after injury (Fig. 4d,e and Extended Data Fig. 5a,c,d).

Overall, our work strongly suggests that the selective increase in wild-type cell divisions during injury repair of mosaic skin limits the expansion of mutant cells expressing different oncogenic variants of the Ras gene family.

We found that injury repair in mosaic skin triggers a specific competitive advantage of wild-type cells but not KrasG12D/+ and HrasG12V/+ cells. However, it remained unclear whether the increased fraction of dividing wild-type cells per se prevented the expansion of Ras mutant cells. To identify the mechanisms responsible for this competitive switch, we extended our scRNA-seq analysis to epithelial cells, fibroblasts and immune cells in our wild-type and HrasG12V/+ models (Extended Data Figs. 6 and 7), to analyse the expression of growth factors, among other soluble mediators, that influence epithelial cell proliferation. Across all injured models, we found a particularly high enrichment of EGFR ligands in both fibroblasts and epithelial cells when compared to other growth factors (Fig. 5a and Extended Data Figs. 6 and 7). Notably, EGFR is one of the best characterized upstream activators of the Ras pathway and previous studies have shown that EGFR signalling promotes epithelial cell proliferation during injury repair in wild-type skin34,35,36. In agreement with this, our western blot analyses and quantification from wild-type mice showed increased EGFR activation in injured versus uninjured conditions (Fig. 5b; for gel source data, see Supplementary Fig. 1). Additionally, HrasG12V/+ cells displayed a downregulation of the receptor in uninjured and injured skin epithelium, as previously shown in vitro37 (Fig. 5b; for gel source data, see Supplementary Fig. 1). We therefore designed a strategy to inhibit EGFR to selectively reduce the number of dividing wild-type cells after injury. We used the EGFR inhibitor Gefitinib, which we verified repressed the activity of EGFR and its downstream target ERK1/2 during injury repair in wild-type mosaic mice (Extended Data Fig. 8a,b and Methods; for gel source data, see Supplementary Fig. 1). As expected, Gefitinib treatment selectively inhibited the amount of wild-type cell divisions but not those of HrasG12V/+ cells, in HrasG12V/+ mosaic models after injury (Fig. 5c). To assess how Gefitinib treatment affects cell competition during injury repair, we tracked the percentage of surface coverage of HrasG12V/+ cells (tdTomato+) at 3, 7 and 14 days PWI (6, 10 and 17 days, respectively, after tamoxifen-induced mosaicism). We found that EGFR inhibition re-establishes the competitive advantage of HrasG12V/+ cells over their wild-type neighbours during injury repair in mosaic mice (Fig. 5d and Extended Data Fig. 8e). Because this drug treatment broadly affects multiple cell types in addition to the epithelium, we generated a genetic model that targets the expression of the dominant negative (DN) form of EGFR (lacking the cytoplasmic tyrosine kinase domain) to the basal cells of the skin epidermis (Krt14-CreER; Krt14-rtTA; FR-HrasG12V/+; LSL-tdTomato; TRE-Egfr-DN; Krt14-H2B–GFP; Extended Data Fig. 8c,d and Methods). This selective EGFR inhibition in the basal stem cell layer reduced the amount of dividing wild-type cells and re-established the competitive advantage of HrasG12V/+ cells over their wild-type neighbours during injury repair of mosaic mice (Extended Data Fig. 8c,d).

a, Violin plots showing cell scoring based on expression of EGFR ligands (Extended Data Fig. 6i), separated by cell type and genotype. Two-tailed t-test comparing the averages of biological replicates according to conditions (fibroblasts: n = 3 wild-type mosaic and n = 6 HrasG12V/+ mosaic; epithelial cells n = 6 mice per genotype), P values are shown. Internal box plots denote the 25th, 50th and 75th centiles, with whiskers depicting minima and maxima, excluding outliers that are beyond 1.5× the interquartile range. b, Western blot analysis (left) and quantification (middle) of phosphorylated EGFR (p-EGFR) and total EGFR in injured and uninjured conditions (n = 3 mice). Paired, two-tailed t-test. Pairs of coloured dots represent ratios for individual mice. Right, western blot analysis of total EGFR normalized to GAPDH (n = 3 mice). Blots were processed at the same time. Unpaired, two-tailed t-test. c, Quantification of mitotic cells in tdTomato+ and tdTomato− areas. d, Left, representative two-photon revisit images following injury. Right, the increase in tdTomato+ area after tamoxifen induction. In c,d, vehicle, n = 3 mice; Gefitinib, n = 4 mice. e, Quantification of mitotic cells in tdTomato+ and tdTomato− areas in wild-type mosaic (n = 3 mice), HrasG12V/+ mosaic (n = 3 mice) and constitutive p21-null HrasG12V/+ mosaic (n = 4 mice). f, Quantification of p-histone H3-positive cells in tdTomato+ and tdTomato− areas (n = 6 mice). g, Representative two-photon revisit images following injury in wild-type mosaic, HrasG12V/+ mosaic and constitutive p21null-HrasG12V/+ mosaic mice. h, The increase in tdTomato+ area. Unpaired, ordinary one-way ANOVA comparing wild-type mosaic, HrasG12V/+ mosaic and constitutive p21null-HrasG12V/+ mosaic at different time points in the uninjured condition. In g,h, n = 3 wild-type mosaic mice, n = 3 HrasG12V/+ mosaic mice and n = 4 constitutive p21null-HrasG12V/+ mosaic mice. c–f, Unpaired or paired two-tailed t-test for comparisons between different groups or within the same group of mice. P values are shown. At least three independent areas of approximately 300 µm2 were analysed for each mouse (Methods). Data are mean ± s.d.

Source Data

Collectively, our data demonstrate that the EGFR–Ras signalling pathway is required to increase wild-type cell divisions after injury and to prevent HrasG12V/+ cell expansion in mosaic models.

To examine whether the increase in mitotic events in wild-type cells is sufficient to counteract the expansion of Ras mutant cells, we directly manipulated cell proliferation in an uninjured setting using a p21-null genetic model (constitutive p21null mice38; the p21 gene is also known as Cdkn1a). This approach, in the absence of injury, enabled us to more directly assess the role of wild-type cell proliferation in preventing Ras mutant cell expansion, as it lacks the influence of other factors that participate in injury repair (such as recruitment of immune cells and activated fibroblasts). The G1/S phase cyclin-dependent kinase inhibitor p21 is expressed in G1 phase cells to maintain skin epithelium homeostasis39,40, and its genetic ablation leads to increase proliferation in wild-type mouse skin epithelial cells in vitro and in vivo39,40,41,42. Indeed, p-histone H3 immunostaining of a p21-null model revealed a substantially increased number of dividing basal cells, but not dermal cells, in uninjured mice (Extended Data Fig. 8f–h). We hypothesized that HrasG12V/+ cells downregulate p21 expression to promote mitotic events. p21 expression was significantly reduced at both mRNA and protein levels in HrasG12V/+ mosaic models compared to wild-type mosaic models (Extended Data Fig. 9a–d). Therefore, we reasoned that constitutive p21 deletion in the HrasG12V/+ mosaic model may selectively manipulate wild-type cells without affecting HrasG12V/+ cells, which express low levels of p21. To this end, we combined the p21-null model with the tamoxifen-inducible HrasG12V/+ model (Krt14-CreER; FR-HrasG12V/+; LSL-tdTomato; constitutive p21-null; Krt14-H2B–GFP; Methods). p21 loss increased the proliferation of wild-type cells but not that of HrasG12V/+ cells in p21-null HrasG12V/+ mosaic mice, mimicking the injured context (Fig. 5e,f). This p21-null-triggered increased proliferation was sufficient to suppress the competitive advantage of HrasG12V/+ cells over time in uninjured mice (Fig. 5g,h and Extended Data Fig. 8i), recapitulating the effects of injury repair in HrasG12V/+ mosaic mice. The fraction of differentiating wild-type cells was also accelerated by the loss of p21, which explains why wild-type cells did not outcompete HrasG12V/+ cells in the p21-null HrasG12V/+ mosaic model (Extended Data Fig. 8j).

To explore the molecular mechanism that specifically increases the division of wild-type cells but not HrasG12V/+ cells in p21-null and injury models, we probed the activation status of key regulators of cell proliferation downstream of Ras43, ERK1/2 and AKT protein kinases (Extended Data Fig. 10a,b; for gel source data, see Supplementary Fig. 1). In uninjured, wild-type skin, there was a lower level of activated ERK1/2 (p-ERK1/2) when compared with HrasG12V/+ mosaic and HrasG12V/+ max models, as expected (Extended Data Fig. 10a). We observed increased levels of p-ERK1/2 in wild-type skin that was undergoing injury repair or lacked p21 (Extended Data Fig. 10a). However, the increase in p-ERK1/2 was not significantly different in HrasG12V/+ mosaic and HrasG12V/+ max models in uninjured and injured conditions (Extended Data Fig. 10a). Moreover, p-ERK1/2 levels were similar after injury in all three models (wild-type, HrasG12V/+ mosaic and HrasG12V/+ max) (Extended Data Fig. 10a). Of note, we performed p-ERK1/2 immunostaining to explore the ERK1/2 activity in individual cells. This experiment confirmed the data from the western blot analyses, with wild-type clones selectively increasing the number of cells expressing p-ERK1/2 upon wounding (Extended Data Fig. 10c,d). Finally, the activation of AKT, another downstream target of the Ras pathway, was not significantly affected by injury or p21 loss (Extended Data Fig. 10b).

Overall, these data suggest that injury repair or loss of p21 specifically increase the activity of ERK1/2, a downstream effector of Ras, in wild-type cells to increase the fraction of dividing cells, enabling them to effectively suppress the competitive advantage of oncogenic Ras mutant cells in mosaic mice.

Healthy tissues, including skin, harbour a number of somatic mutations, some of which are in known tumour driver genes1,2,44,45. Models have shown that tumours can arise from the accumulation of multiple mutations or from a lower mutational burden cooperating with additional exogenous insults, such as injury3,4,5,6,7,8,9,10,11,24,25. Here we found that Ras mutant cells break tissue architecture during injury repair only when they represent nearly all cells in the basal stem cell layer of the skin epidermis. By contrast, when Ras mutant cells coexist with wild-type neighbours, injury selectively activates the endogenous proliferation programme in wild-type cells only, which counteracts the expansion of the Ras mutant cells (Extended Data Fig. 10e). Specifically, after an initial spike following injury induction, the number of dividing wild-type cells equalizes to the level of Ras mutant cells, higher than the level of wild-type cell division in homogeneous wild-type models. Neighbouring wild-type cells exert a powerful defensive mechanism, even in the presence of a higher mutational burden from the most aggressive isoform of the Ras family, Kras.

We found that although injury repair is coordinated by various growth factors and receptors, the EGFR–Ras signalling pathway emerges as key for increasing the number of dividing wild-type cells to suppress the expansion of Ras mutant cells. In the absence of injury, the constitutive loss of p21 also leads to a selective increase in the fraction of dividing wild-type cells in uninjured mosaic skin, recapitulating the responses to injury and highlighting the sufficiency of an increased number of wild-type cell divisions as a protective mechanism against the expansion of Ras mutant cells. Thus, genetic and environmental mechanisms promote a competitive advantage for wild-type cells, whereas Ras mutant cells are insensitive. Our data suggest that the HrasG12V mutation renders cells insensitive to pro-proliferative stimuli mediated by the EGFR–Ras pathway during injury repair, in part because HrasG12V cells already exhibit high levels of ERK1/2 activation, which promotes cell proliferation. On the basis of this observation, we propose that the fold change in Ras pathway activation before and after injury contributes to the selective capacity of wild-type cells to respond to pro-proliferative stimuli during injury repair that is already maximized in HrasG12V/+ cells (Extended Data Fig. 10e).

Our findings have broad implications, given that as well as in uninjured skin, Ras mutant cells have competitive advantages over wild-type cells in other uninjured tissues such as intestinal crypts and blood46,47. Recent studies in a single-layer epithelium in vitro and in vivo have shown that RasG12V mutant cells are apically extruded when surrounded by wild-type epithelial cells48,49,50. In the stratified epithelium of skin epidermis, our studies and others have shown that Ras mutant cells expand to outcompete wild-type neighbours and integrate in healthy tissue, suggesting a different mode of cell competition compared with the systems above18,19,20. We reconcile the above evidence by noting that different tissues preserve and maintain their specific architecture through distinct cell behaviours. Endogenous behaviours for maintaining monolayer epithelia are proliferation and extrusion, whereas multilayer epithelia rely on proliferation and delamination–differentiation. Thus, in a stratified epithelium, differentiation is analogous to extrusion in monolayer epithelia. Consistent with this reasoning, our evidence suggests that apoptosis, an ectopic behaviour in the adult skin epithelium, is not involved in the competition between wild-type and HrasG12V/+ cells.

Traditional therapeutic approaches used for cancer treatment involve suppressing the proliferation of both mutant and wild-type cells. Although these approaches restrain tumour expansion, they also impair the opportunity for the tissue to deploy natural defences, such as the selective promotion of wild-type cell proliferation. The next step towards an effective therapeutic treatment would be to determine how to promote the proliferative advantage bestowed on wild-type cells in the injury environment or in the pro-proliferative p21-null state. Our data suggest that in precancerous states and to prevent tumour relapse, EGFR activation—such as through EGF treatment—might provide a competitive advantage to wild-type cells in the presence of neighbours expressing the constitutively active form of a Ras oncogene. Future studies could investigate how the competitive advantage of wild-type cells may be leveraged in the setting of field cancerization with more complex mutant clones, which would represent a powerful treatment approach for field therapy in precancerous conditions such as actinic keratosis, the precursor of cSCC.

Collectively, this work provides a way forward for future research and clinical applications to focus on the mechanisms that empower wild-type cells in the competition with mutated neighbouring cells.

We used an outbred mouse strain (CD1) that is less susceptible to generating papilloma or squamous cell carcinomas than more inbred mouse strains used in other studies11,12. Mice were generated by interbreeding mice carrying the following alleles: Krt14-CreER51 and FR-HrasG12V/+ (ref. 52), constitutive p21 (Cdkn1a) loss of function38 (JAX stock no. 016565), LSL-tdTomato53 (JAX stock no. 007909), Krt14-H2B–GFP23, Krt14-rtTA54 (JAX stock no. 008099), TRE-EGFR-DN55 (JAX stock no. 010575) and LSL-KrasG12D/+ (ref. 56). Mice expressing Krt14-H2B–GFP were used to track epithelial cells with the two-photon microscope. The tdTomato reporter line was used to visualize CreER-driven recombination upon tamoxifen injection. CreER/LoxP lines are well documented to exhibit a certain degree of tamoxifen-independent (leaky) Cre activity over time. To account for the leakiness of the system in our experiments and to correctly interpret the data, we always compared wild-type mosaic models to Ras-mosaic models as well as tracked the same cells in the same animals over time. Mice from experimental and control groups were randomly selected for either sex. No blinding was done. All procedures involving animal subjects were performed under the approval of the Institutional Animal Care and Use Committee (IACUC) of the Yale School of Medicine. The mice were sacrificed if tumours reached 1 cm3 (not allowed by IACUC) or if they presented signs of distress or weight loss. The tumour size limit was not exceeded in any of the experiments.

To induce CreER-driven recombination, mice were administered a single dose of 100 μg (mosaic induction) or 2 mg (maximal induction) tamoxifen (Sigma T5648-5G in corn oil) at postnatal day 19 by intraperitoneal injection (this time is designated day 0 for experiments). To induce rtTA-driven induction of EGFR-DN, mice were administered 2% of doxycycline (Sigma D9891) and 2% sucrose (Sigma S9378) in drinking water. All time courses began 6 days after tamoxifen injection. Gefitinib (ZD1839-Selleckchem) was resuspended in water with 0.5% (w/v) methylcellulose and 0.2% (v/v) Tween-80 (vehicle) and was administered orally (200 mg kg−1 body weight) starting 2 days before wound induction until 14 days PWI (mice were not treated at day 7 and day 8 PWI).

At postnatal day 21, mice were anaesthetized by intraperitoneal injection of a ketamine and xylazine mix (100 mg kg−1 and 10 mg kg−1, respectively in phosphate-buffered saline). Once the anaesthetized mouse did not physically respond to a noxious stimulus, a punch biopsy was performed using a 4-mm-diameter punch biopsy tool (Integra Miltex Standard Biopsy Punches). The punch biopsy tool was used to make a circular full-thickness injury on the dorsal side of a mouse ear or back skin. The injury did not penetrate the entire ear but remained above the cartilage. The skin epithelium in the mouse ear was chosen in this study for its accessibility to two-photon imaging and revisits over time. For recovery from the wound procedure Meloxicam (Metacam Loxicom) was administered via subcutaneous injection (0.3 mg kg−1).

Large-scale production and concentration of lentivirus expressing CreER (LV-CreER) was performed as previously described57,58. Detailed descriptions of in utero-guided lentiviral transduction can be found elsewhere57,59. To induce LV-CreER-mediated recombination, a maximal dose of tamoxifen (2 mg) was intraperitoneal injected at postnatal day 19.

Mice were anaesthetized by intraperitoneal injection of ketamine and xylazine cocktail mix (100 mg kg−1 and 10 mg kg−1, respectively in phosphate-buffered saline) and then anaesthesia was maintained throughout the course of the experiment with the delivery of vaporized isoflurane by a nose cone as previously described60. Image stacks were acquired with a LaVision TriM Scope II (LaVision Biotec) laser scanning microscope equipped with a tunable Two-photon Chameleon Vision II (Coherent) Ti:Sapphire laser. To acquire serial optical sections, a laser beam 940 nm was focused through a 20× water immersion lens (NA 1.0; Zeiss) or a 25× water immersion lens (NA 1.1; Nikon) and scanned with a pixel size of 0.49 × 0.49 μm2 or 0.43 × 0.43 μm2 at 800 Hz. Z-stacks were acquired in 2–3 μm steps to image a total depth of 90–900 μm of the tissue. ImSpector v7.5.2 (LaVision Biotec) was used for 3D image acquisition. To visualize large areas, 5–12 tiles of optical fields were imaged using a motorized stage to automatically acquire sequential fields of view as previously described60. Mice were imaged at different time points after tamoxifen treatment and injury induction as indicated. To revisit the same area of the skin epidermis, organizational clusters of hair follicles and vasculature were used as landmarks.

To prepare whole mounts of mouse ear, skin was separated from connective tissue using forceps and incubated in 4% paraformaldehyde (PFA) at 37 °C for 4 h. To prepare the epidermal preparation, the skin separated from the connective tissue was incubated with Dispase (5 mg ml−1, Roche) for 10 min at 37 °C and then the intact sheet of epidermis was gently peeled away from the dermis. The epidermal preparations were fixed in 4% PFA for 45 min at room temperature. Immunostaining was performed on whole mounts and epidermal preparations blocked with 5% normal goat serum/1% BSA/2% Triton X-100/PBS at room temperature. For tissue-section analysis, mouse ears were fixed in 4% PFA for 1 h at room temperature and then embedded in optimal cutting temperature (OCT; Tissue Tek). Frozen OCT blocks were sectioned at 10 μm. The skin preparations were incubated with primary antibodies (active caspase-3 (AF835-R&D Systems) 1:300, p-histone H3 (06-570-Millipore) 1:300, keratin-6A (905701-BioLegend) 1:500, p21 (also known as Cdkn1a) (ab188224-Abcam) 1:50; p-p44/42 MAPK(Erk1/2) (4370-Cell Signaling) 1:300, Alexa Fluor 647 Phalloidin (A22287-Thermofisher) 1:200 and keratin-10 (03-GP-K10-ARP) 1:200 diluted in blocking buffer (∼66 h at 37 °C for whole mounts and 4 °C overnight for epidermal preparations and OCT sections). Secondary antibodies (A-21071-Invitrogen goat anti-rabbit IgG (H + L) secondary antibody, Alexa Fluor 633, A-21105-Invitrogen goat anti-guinea pig IgG (H + L) highly cross-adsorbed secondary antibody, Alexa Fluor 633, A-21206, donkey anti-rabbit IgG (H + L) highly cross-adsorbed, Alexa Fluor 488 and A-10042 donkey anti-rabbit IgG (H+L) highly cross-adsorbed secondary antibody, Alexa Fluor 568) were diluted 1:300 in blocking buffer and applied to skin sections (∼66 h at 37 °C for whole mounts or 1 h at room temperature for epidermal preparations and OCT sections). DAPI was added to label nuclei. Image stacks of whole-mount, epidermal preparation and OCT section immunostaining were acquired with the two-photon microscope described above with an additional tunable two-photon Chameleon Discovery (Coherent) Ti:Sapphire laser. To acquire serial optical sections, laser wavelengths of 800 nm, 880 nm (Vision II) and 1080 nm (Discovery) were focused through a 20× water immersion lens (NA 1.0; Zeiss) or a 25× water immersion lens (NA 1.1; Nikon) and scanned with a pixel size of 0.49 × 0.49 μm2 or 0.43 × 0.43 μm2 at 800 Hz. Image stacks of epidermal preparation immunolabelled for p-ERK1/2 and phalloidin were acquired with confocal microscope Zeiss LSM 980, objective Zeiss 20× (0.8 NA Dry) and scanned with a pixel size of 0.20 × 0.20 μm2 with Software ZEN (blue edition). Z-stacks were acquired in 1-μm (whole mounts) or 2-μm (epidermal preparations and OCT sections) steps to image total depth of the samples.

Senescence was measured in epidermal preparation of uninjured and injured (day 3 PWI) mouse ear skin and in mouse pancreas and kidney, as positive controls, with a Beta Galactosidase Assay Kit from Abcam (ab287846). This kit uses the fluorogenic fluorescein digalactoside galactosidase substrate, which, upon cleavage by β-galactosidase, generates a fluorescent product that can be measured with an ELISA plate reader. In brief, ~10 mg of the tissues described above was lysed with protein lysis buffer (included in the kit) and then the kit guidelines were followed to prepare samples, positive controls provided by the kit and the standard curve. The fluorescence in each sample was quantified with a GloMax Plate reader (Promega) at 475 nm excitation and 500–550 nm emission at two different time points (30 min apart). β-Gal levels in each sample were calculated using a β-galactosidase standard curve.

Raw two-photon image stacks were analysed in ImageJ (NIH Image, 1.53c) or IMARIS (v. 9.9.1, Oxford Instruments). ImageJ was used to draw the boundaries between tdTomato+ and tdTomato− areas in the basal stem cell layers of the skin epidermis and to measure the percentage of coverage of these areas. The average of the percentage of coverage of tdTomato+ cells of three areas (294 μm2 each) taken at different distances from the wound and randomly in the uninjured condition were calculated for each mouse at every time point. Then, the percentage increase of tdTomato+ area over time was represented in the graphs. To be reproducible in the measurements, the areas quantified in the injured ears were not localized close to the wound edge because of the increased thickness of the epithelium in that region that prevents an accurate isolation of the basal stem cell layer due to reduced image resolution. tdTomato+ and tdTomato− regions for each uninjured or injured ear were revisited and quantified over time.

We measured the number of events (mitotic figures, nuclear fragmentation events, or cells positive for p-histone H3, keratin-10 and active caspase-3 immunostaining) per unit of surface (1 μm2) of tdTomato+ or tdTomato− areas and then multiplied that value by the total surface (294 μm2) to compare wild-type and mutant cell populations. ~900 μm2 were analysed for each uninjured or injured mouse ear. We measured the number of cells positive for p-ERK1/2 per unit of surface (1 μm2) of tdTomato+ or tdTomato− areas and then multiplied that value by the total surface (150 μm2) to compare wild-type and mutant cell populations. An area of approximately 600 μm2 was analysed for each uninjured or injured mouse ear.

To quantify the thickness of the skin epithelium, we used IMARIS (v. 9.9.1). By utilizing the second harmonic collagen signal in the dermis and, when absent near the wound, the basal epithelial cell layer, we created a surface to approximate the basement membrane of the epithelium. When unable to visualize the basal–dermal interface near the wound due to excessive epithelia thickness, such as in the case of HrasG12V/+ max or KrasG12D/+ max, we created a surface along the bottom of the stack to ensure our thickness measurements were stringent in that and reflected only as much tissue as we could effectively image. From this surface we could extract the distance to the top of the cornified layer around the injured area of the epithelium and visually depict the tissue thickness with an intensity heat map. Using MatLab (v. R2018a) we could extract individual pixel intensity values from this heat map that directly correlate to tissue thickness and plot them based on relative distance from the wound edge.

To quantify the nuclear signal of p21 in the basal stem cell layer of the skin epidermis, we used IMARIS (v. 9.9.1). In IMARIS, surfaces were created using the Krt14-H2B–GFP signal to isolate a mask of the p21 signal within the basal epidermal stem cell nuclei. A maximum intensity projection of this mask was used to quantify the mean fluorescence intensity of the p21 signal within each individual nucleus. The mean fluorescence intensity of the p21 signal of each cell was normalized for background by subtracting the average p21 fluorescent intensity of mitotic cells within the field of view, as these cells would be negative for p21.

After the sacrifice of wild-type, HrasG12V/+ mosaic and HrasG12V/+ max models at 6 days after tamoxifen injection (3 days PWI), the uninjured and injured ears of each mouse were cut in small pieces with a punch biopsy of 8 mm in diameter (the wound was kept at the centre the 8 mm biopsy to mostly isolate cells involved in injury repair). The ear epidermis was dissociated from the dermis and incubated in 0.25 % Trypsin at 37 °C for 30 min. The epidermal preparation was placed in a 70 μm cell strainer, smashed with a piston and rinsed three times with PBS + 0.04% BSA. The flow-through was subsequently filtered through a 40 μm cell strainer, spun down and resuspended in 300–400 μl of PBS + 0.04% BSA. The viability of the cell suspension was determined using trypan blue. To prepare the single-cell library, the cellular suspensions were counted and diluted to a final concentration of 1,200 cells per μl in PBS/0.04% BSA and then loaded on a Chromium Controller to generate single-cell gel bead emulsions, targeting 3′. Single-cell 3′ RNA-seq libraries were generated according to the manufacturer’s instructions (Chromium Single Cell 3′ Reagent v3 Chemistry Kit, 10X Genomics). Libraries were sequenced to an average depth of ∼20,000 reads per cell on an Illumina Novaseq 6000 system.

Single-cell data from each sample—that is, all wild-type and HrasG12V/+ mosaic and HrasG12V/+ max, uninjured and injured conditions (24 independent samples, Extended Data Figs. 4a, 6a and 7a) were first processed with SoupX61 (https://github.com/constantAmateur/SoupX) to remove barcodes that most probably represent ambient RNA as opposed to whole cells, using the algorithm’s automated method. The resulting matrix was then processed with the Seurat package62 (v.3, https://satijalab.org/seurat/index.html), to retain genes or features that are detected in at least 3 cells and include cells for which at least 200 genes or features are detected. Additionally, cells expressing greater than 12.5% of mitochondrial transcripts were filtered out as possible dead or dying cells. According to Seurat’s normal workflow, the data was log-normalized and scaled. Linear dimensionality reduction was carried out using principal component analysis and the first 15 principal components were chosen for the downstream analysis steps. Clustering was carried out using Louvain algorithm, for resolution of 0.1. Non-linear dimensionality reduction was carried out by running UMAP. Next, the DoubletFinder63 package (https://github.com/chris-mcginnis-ucsf/DoubletFinder) was used to get rid of barcodes that may represent possible doublets. The resulting cell matrix was normalized, scaled, and re-clustered, using the same steps mentioned above with the Seurat package. Three replicates for each group were then integrated using Seurat’s canonical correlation analysis (CCA). The data were re-clustered as described previously22, and the clusters were annotated using the top 5–10 highly expressed genes in each cluster. To further remove doublets that were not identified by the DoubletFinder algorithm and other contaminating cell populations, infundibulum cells, immune cells, and red blood cells (RBCs) were removed. In brief, using a chosen set of features, each cell in the Seurat object was assigned a score using the AddModuleScore function. The features used for each type of cells are listed below. Infundibulum-specific features: “Sostdc1”, “Aqp3”, “Ptn”, “Fst”, “Aldh3a1”, “Postn”, “Krt17”,“Alcam”, “Apoe”, “Sox9”, “Vdr”, “Nfib”, “App”, “Gsn”, “Hmcn1”, “Cspg4”, “Efnb2”, “Nedd4”, “Adh7”, “Defb6”, “Mgst1”, “Krt79”. Immune cell-specific features: “H2-Aa”, “H2-Ab1”, “H2-Eb1”, “Cd74”, “Ptprc”. RBC-specific features: “Hba-a1”, “Hbb-bs”, “Hba-a2”, “Hbb-bt”, “Bpgm”, “Hebp2”. For infundibulum signature, cells with scores higher than 0.4 were removed, while for immune and RBCs, cells with scores higher than 0.5 were removed. The following cell numbers (per sample and biological replicate (R)) passed QC and constitute the final dataset: HMU: 11,437 (R1), 9,864 (R2), 9,650 (R3); HMW: 9,590 (R1), 10,281 (R2), 8,922 (R3); WTMU: 11,856 (R1), 9,916 (R2), 8,416 (R3); WTMW: 9,181 (R1), 8,747 (R2), 8,320 (R3); HFU: 13,278 (R1), 10,711 (R2), 6,997 (R3); HFW: 9,613 (R1), 8,999 (R2), 7,825 (R3); WTFU: 13,379 (R1), 11,701 (R2), 6,433 (R3); WTFW: 9,649 (R1), 9,374 (R2), 6,350 (R3).

For interfollicular epidermis (IFE) keratinocytes, all mosaic samples were integrated and annotated using Scanpy (1.6-1.9)64. In brief, raw counts for the selected IFE cells were log-normalized and cell cycle stages were scored (sc.pp.score_genes_cell_cycle) based on a gene list from65. Biological replicate batches were corrected with bbknn (1.4.1)66. Next, the selected cells were scored for stress, immune and infundibulum related gene expression signatures (see the notebooks on GitHub: https://github.com/kasperlab/Gallini_et_al_2023_Nature), classified with a Gaussian mixture model (scikit-learn, 0.24.267) and positive cells were filtered out. Similarly, classification was performed to annotate cells based on tdTomato and GFP expression. The remaining healthy IFE keratinocytes were then mapped (sc.tl.ingest) onto the characterized IFE differentiation trajectory and annotated accordingly based on the basal-suprabasal status and commitment that have been previously defined22. Finally, appropriate IFE groups were integrated with CCA as described above, using the Seurat package. Differential gene expression analysis was carried out between corresponding datasets and cell types.

For fibroblast and immune cell characterization, all mosaic and max datasets were analysed in Scanpy, with similar preprocessing and batch effect removal as for keratinocytes. Mixed cell populations were removed based on shared gene expression signatures. Wound-related cells were annotated based on Leiden clustering, sample type and wound-related gene expression signatures with additional confirmation by differential abundance testing using miloR (1.2.0)68. In brief, a Milo graph was built using the integrated dataset with the following parameters: k = 20, d = 30 and differential abundance was tested for the injury condition. Differential gene expression analysis was performed with scanpy.tl.rank_genes_groups function using Wilcoxon rank-sum test and Holm–Sidak correction for multiple comparisons. Gene set enrichment analysis was performed using the enrichr method in GSEAPY package (v0.12) with Gene Ontology biological process 2021 gene sets69,70,71,72,73.

Uninjured and injured sections of ear skin were fixed in 10% neutral formalin for 24 h and stored in 70% ethanol until paraffin embedding. Haematoxylin and eosin (H&E)-stained skin sections were used for histopathology analysis. Images were taken using an Olympus BX61 microscope equipped with a SPOT flex 15.2 64-Mp shifting pixel camera, 4×, 10× and 20× objectives, and SPOT v 5.2 software.

Uninjured and injured ear skin were lysed with ice-cold RIPA buffer (Pierce) supplemented with cOmplete Protease Inhibitor Cocktail and PhosSTOP (Sigma) and centrifuged at maximum speed for 30 min to collect lysates. Protein concentration was measured with the BCA protein assay (Pierce). An aliquot of 20–30 μg of total protein per sample was loaded into 7.5 or 10% Mini-PROTEAN TGX Precast Protein Gels (BioRad) and separated by SDS–PAGE. Proteins were transferred to PVDF membranes (BioRad). The following rabbit primary antibodies were used at the given concentrations; p-p44/42 MAPK (ERK1/2) (Thr202/Tyr204) (1:500, Cell Signaling 9101), p44/42 MAPK (ERK1/2) (1:500, Cell Signaling 4695), p-EGFR (Tyr1068) (1:100, Cell Signaling 2234), EGFR (1:100 Cell Signaling 4267; Extended Data Fig. 8b), EGFR (1:100, Cell Signaling 2232; Fig. 5b), p-AKT (Ser473) (1:200, Cell Signaling 4060), AKT (1:200, Cell Signaling 9262) and GAPDH (14C10) (1:500, Cell Signaling 2118). An anti-rabbit IgG HRP (1:500, Cell Signaling 7074) secondary antibody was used. Western blot analyses were performed on whole ear skin at 6 days after tamoxifen injection (3 days PWI).

Statistical analyses were performed using an unpaired, two-tailed Student’s t-test for comparison between different groups of mice. Paired, two-tailed Student’s t-tests were used for comparison between tdTomato+ and tdTomato− populations in the same group of mice. Unpaired, ordinary one-way ANOVA was used for comparison between mice with three different genotypes. Statistical analyses were performed using Prism (v. 9) as indicated in the figure legends. Gene expression differences between different conditions from scRNA-seq data were performed with Student’s t-test and Holm–Sidak correction for multiple comparisons. P values of less than 0.05 were considered statistically significant (*P < 0.05, **P < 0.005, ***P < 0.0005, ****P < 0.0001). n is defined for each experiment, and always indicates the number of mice used for each condition examined. Box plots within violin plots denote the 25th, 50th and 75th quartiles, with whiskers depicting the minima and maxima of the data, excluding outliers that are beyond 1.5× the interquartile range.

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

All data from this study are available from the authors on request. The raw data files of the scRNA-seq analyses reported in the manuscript have been uploaded to the Gene Expression Omnibus under accession GSE195892. Previously published scRNA-seq data that were used for reference are available under accession codes GSE152044, GSE129218 and GSE67602. Source data are provided with this paper.

Annotated and analysed sequencing data have been deposited at Zenodo (https://doi.org/10.5281/zenodo.7768108) and analysis notebooks have been uploaded to GitHub (https://github.com/kasperlab/Gallini_et_al_2023_Nature).

Martincorena, I. et al. Tumor evolution. High burden and pervasive positive selection of somatic mutations in normal human skin. Science 348, 880–886 (2015).

Article ADS CAS PubMed PubMed Central Google Scholar

Fowler, J. C. et al. Selection of oncogenic mutant clones in normal human skin varies with body site. Cancer Discov. 11, 340–361 (2021).

Article CAS PubMed Google Scholar

Bailleul, B. et al. Skin hyperkeratosis and papilloma formation in transgenic mice expressing a ras oncogene from a suprabasal keratin promoter. Cell 62, 697–708 (1990).

Article CAS PubMed Google Scholar

Cannon, R. E. et al. Kinetics of wound-induced v-Ha-ras transgene expression and papilloma development in transgenic Tg.AC mice. Mol. Carcinog. 20, 108–114 (1997).

3.0.CO;2-5" data-track-action="article reference" href="https://doi.org/10.1002%2F%28SICI%291098-2744%28199709%2920%3A1%3C108%3A%3AAID-MC12%3E3.0.CO%3B2-5" aria-label="Article reference 4" data-doi="10.1002/(SICI)1098-2744(199709)20:13.0.CO;2-5">Article CAS PubMed Google Scholar

Greenhalgh, D. A. et al. Induction of epidermal hyperplasia, hyperkeratosis, and papillomas in transgenic mice by a targeted v-Ha-ras oncogene. Mol. Carcinog. 7, 99–110 (1993).

Article CAS PubMed Google Scholar

Vassar, R., Hutton, M. E. & Fuchs, E. Transgenic overexpression of transforming growth factor α bypasses the need for c-Ha-ras mutations in mouse skin tumorigenesis. Mol. Cell. Biol. 12, 4643–4653 (1992).

CAS PubMed PubMed Central Google Scholar

Dominey, A. M. et al. Targeted overexpression of transforming growth factor alpha in the epidermis of transgenic mice elicits hyperplasia, hyperkeratosis, and spontaneous, squamous papillomas. Cell Growth Differ. 4, 1071–1082 (1993).

CAS PubMed Google Scholar

Oki-Idouchi, C. E. & Lorenzo, P. S. Transgenic overexpression of RasGRP1 in mouse epidermis results in spontaneous tumors of the skin. Cancer Res. 67, 276–280 (2007).

Article CAS PubMed PubMed Central Google Scholar

Kasper, M. et al. Wounding enhances epidermal tumorigenesis by recruiting hair follicle keratinocytes. Proc. Natl Acad. Sci. USA 108, 4099–4104 (2011).

Article ADS CAS PubMed PubMed Central Google Scholar

Wong, S. Y. & Reiter, J. F. Wounding mobilizes hair follicle stem cells to form tumors. Proc. Natl Acad. Sci. USA 108, 4093–4098 (2011).

Article ADS CAS PubMed PubMed Central Google Scholar

Page, M. E., Lombard, P., Ng, F., Gottgens, B. & Jensen, K. B. The epidermis comprises autonomous compartments maintained by distinct stem cell populations. Cell Stem Cell 13, 471–482 (2013).

Article CAS PubMed PubMed Central Google Scholar

Huang, P. Y. et al. Lgr6 is a stem cell marker in mouse skin squamous cell carcinoma. Nat. Genet. 49, 1624–1632 (2017).

Article CAS PubMed PubMed Central Google Scholar

White, A. C. et al. Defining the origins of Ras/p53-mediated squamous cell carcinoma. Proc. Natl Acad. Sci. USA 108, 7425–7430 (2011).

Article ADS CAS PubMed PubMed Central Google Scholar

Spencer, J. M., Kahn, S. M., Jiang, W., DeLeo, V. A. & Weinstein, I. B. Activated ras genes occur in human actinic keratoses, premalignant precursors to squamous cell carcinomas. Arch. Dermatol. 131, 796–800 (1995).

Article CAS PubMed Google Scholar

van der Schroeff, J. G., Evers, L. M., Boot, A. J. & Bos, J. L. Ras oncogene mutations in basal cell carcinomas and squamous cell carcinomas of human skin. J. Invest. Dermatol. 94, 423–425 (1990).

Article PubMed Google Scholar

Quintanilla, M., Brown, K., Ramsden, M. & Balmain, A. Carcinogen-specific mutation and amplification of Ha-ras during mouse skin carcinogenesis. Nature 322, 78–80 (1986).

Article ADS CAS PubMed Google Scholar

Brown, K., Strathdee, D., Bryson, S., Lambie, W. & Balmain, A. The malignant capacity of skin tumours induced by expression of a mutant H-ras transgene depends on the cell type targeted. Curr. Biol. 8, 516–524 (1998).

Article CAS PubMed Google Scholar

Pineda, C. M. et al. Hair follicle regeneration suppresses Ras-driven oncogenic growth. J. Cell Biol. 218, 3212–3222 (2019).

Article CAS PubMed PubMed Central Google Scholar

Cai, E. Y. et al. Selective translation of cell fate regulators mediates tolerance to broad oncogenic stress. Cell Stem Cell 27, 270–283.e277 (2020).

Article MathSciNet CAS PubMed PubMed Central Google Scholar

Sandoval, M., Ying, Z. & Beronja, S. Interplay of opposing fate choices stalls oncogenic growth in murine skin epithelium. eLife 10, e54618 (2021).

Article CAS PubMed PubMed Central Google Scholar

Mesa, K. R. et al. Niche-induced cell death and epithelial phagocytosis regulate hair follicle stem cell pool. Nature 522, 94–97 (2015).

Article ADS CAS PubMed PubMed Central Google Scholar

Cockburn, K. et al. Gradual differentiation uncoupled from cell cycle exit generates heterogeneity in the epidermal stem cell layer. Nat. Cell Biol. 24, 1692–1700 (2022).

Article CAS PubMed PubMed Central Google Scholar

Tumbar, T. et al. Defining the epithelial stem cell niche in skin. Science 303, 359–363 (2004).

Article ADS CAS PubMed Google Scholar

Hansen, L. A. & Tennant, R. Focal transgene expression associated with papilloma development in v-Ha-ras-transgenic TG.AC mice. Mol. Carcinog. 9, 143–154 (1994).

Article CAS PubMed Google Scholar

Leder, A., Kuo, A., Cardiff, R. D., Sinn, E. & Leder, P. v-Ha-ras transgene abrogates the initiation step in mouse skin tumorigenesis: effects of phorbol esters and retinoic acid. Proc. Natl Acad. Sci. USA 87, 9178–9182 (1990).

Article ADS CAS PubMed PubMed Central Google Scholar

Park, S. et al. Tissue-scale coordination of cellular behaviour promotes epidermal wound repair in live mice. Nat. Cell Biol. 19, 155–163 (2017).

Article CAS PubMed PubMed Central Google Scholar

Aragona, M. et al. Defining stem cell dynamics and migration during wound healing in mouse skin epidermis. Nat. Commun. 8, 14684 (2017).

Article ADS PubMed PubMed Central Google Scholar

Drosten, M., Lechuga, C. G. & Barbacid, M. Ras signaling is essential for skin development. Oncogene 33, 2857–2865 (2014).

Article CAS PubMed Google Scholar

Moreno, E., Valon, L., Levillayer, F. & Levayer, R. Competition for space induces cell elimination through compaction-driven ERK downregulation. Curr. Biol. 29, 23–34.e28 (2019).

Article CAS PubMed PubMed Central Google Scholar

Calautti, E., Li, J., Saoncella, S., Brissette, J. L. & Goetinck, P. F. Phosphoinositide 3-kinase signaling to Akt promotes keratinocyte differentiation versus death. J. Biol. Chem. 280, 32856–32865 (2005).

Article CAS PubMed Google Scholar

van der Weyden, L. et al. Acute sensitivity of the oral mucosa to oncogenic K-ras. J. Pathol. 224, 22–32 (2011).

Article PubMed PubMed Central Google Scholar

Lapouge, G. et al. Identifying the cellular origin of squamous skin tumors. Proc. Natl Acad. Sci. USA 108, 7431–7436 (2011).

Article ADS CAS PubMed PubMed Central Google Scholar

Vitale-Cross, L., Amornphimoltham, P., Fisher, G., Molinolo, A. A. & Gutkind, J. S. Conditional expression of K-ras in an epithelial compartment that includes the stem cells is sufficient to promote squamous cell carcinogenesis. Cancer Res. 64, 8804–8807 (2004).

Article CAS PubMed Google Scholar

Chen, X. et al. IL-17R–EGFR axis links wound healing to tumorigenesis in Lrig1+ stem cells. J. Exp. Med. 216, 195–214 (2019).

Article CAS PubMed PubMed Central Google Scholar

Repertinger, S. K. et al. EGFR enhances early healing after cutaneous incisional wounding. J. Invest. Dermatol. 123, 982–989 (2004).

Article CAS PubMed Google Scholar

Takazawa, Y. et al. Notch down-regulation in regenerated epidermis contributes to enhanced expression of interleukin-36α and suppression of keratinocyte differentiation during wound healing. J. Dermatol. Sci. 79, 10–19 (2015).

Article CAS PubMed Google Scholar

Derer, S. et al. Oncogenic KRAS impairs EGFR antibodies’ efficiency by C/EBPβ-dependent suppression of EGFR expression. Neoplasia 14, 190–205 (2012).

Article CAS PubMed PubMed Central Google Scholar

Deng, C., Zhang, P., Harper, J. W., Elledge, S. J. & Leder, P. Mice lacking p21CIP1/WAF1 undergo normal development, but are defective in G1 checkpoint control. Cell 82, 675–684 (1995).

Article CAS PubMed Google Scholar

Lee, J. et al. Runx1 and p21 synergistically limit the extent of hair follicle stem cell quiescence in vivo. Proc. Natl Acad. Sci. USA 110, 4634–4639 (2013).

Article ADS CAS PubMed PubMed Central Google Scholar

Haensel, D. et al. Defining epidermal basal cell states during skin homeostasis and wound healing using single-cell transcriptomics. Cell Rep. 30, 3932–3947.e3936 (2020).

Article CAS PubMed PubMed Central Google Scholar

Topley, G. I., Okuyama, R., Gonzales, J. G., Conti, C. & Dotto, G. P. p21WAF1/Cip1 functions as a suppressor of malignant skin tumor formation and a determinant of keratinocyte stem-cell potential. Proc. Natl Acad. Sci. USA 96, 9089–9094 (1999).

Article ADS CAS PubMed PubMed Central Google Scholar

Bedelbaeva, K. et al. Lack of p21 expression links cell cycle control and appendage regeneration in mice. Proc. Natl Acad. Sci. USA 107, 5845–5850 (2010).

Article ADS CAS PubMed PubMed Central Google Scholar

Dumesic, P. A., Scholl, F. A., Barragan, D. I. & Khavari, P. A. Erk1/2 MAP kinases are required for epidermal G2/M progression. J. Cell Biol. 185, 409–422 (2009).

Article CAS PubMed PubMed Central Google Scholar

Martincorena, I. et al. Somatic mutant clones colonize the human esophagus with age. Science 362, 911–917 (2018).

Article ADS CAS PubMed PubMed Central Google Scholar

Suda, K. et al. Clonal expansion and diversification of cancer-associated mutations in endometriosis and normal endometrium. Cell Rep. 24, 1777–1789 (2018).

Article CAS PubMed Google Scholar

Sabnis, A. J. et al. Oncogenic Kras initiates leukemia in hematopoietic stem cells. PLoS Biol. 7, e59 (2009).

Article PubMed Google Scholar

Vermeulen, L. et al. Defining stem cell dynamics in models of intestinal tumor initiation. Science 342, 995–998 (2013).

Article ADS CAS PubMed Google Scholar

Sato, N. et al. The COX-2/PGE2 pathway suppresses apical elimination of RasV12-transformed cells from epithelia. Commun. Biol. 3, 132 (2020).

Article PubMed PubMed Central Google Scholar

Sasaki, A. et al. Obesity suppresses cell-competition-mediated apical elimination of RasV12-transformed cells from epithelial tissues. Cell Rep. 23, 974–982 (2018).

Article CAS PubMed PubMed Central Google Scholar

Hill, W. et al. EPHA2-dependent outcompetition of KRASG12D mutant cells by wild-type neighbors in the adult pancreas. Curr. Biol. 31, 2550–2560.e2555 (2021).

Article CAS PubMed PubMed Central Google Scholar

Vasioukhin, V., Degenstein, L., Wise, B. & Fuchs, E. The magical touch: genome targeting in epidermal stem cells induced by tamoxifen application to mouse skin. Proc. Natl Acad. Sci. USA 96, 8551–8556 (1999).

Article ADS CAS PubMed PubMed Central Google Scholar

Chen, X. et al. Endogenous expression of HrasG12V induces developmental defects and neoplasms with copy number imbalances of the oncogene. Proc. Natl Acad. Sci. USA 106, 7979–7984 (2009).

Article ADS CAS PubMed PubMed Central Google Scholar

Madisen, L. et al. A robust and high-throughput Cre reporting and characterization system for the whole mouse brain. Nat. Neurosci. 13, 133–140 (2010).

Article CAS PubMed Google Scholar

Nguyen, H., Rendl, M. & Fuchs, E. Tcf3 governs stem cell features and represses cell fate determination in skin. Cell 127, 171–183 (2006).

Article CAS PubMed Google Scholar

Roh, M., Paterson, A. J., Asa, S. L., Chin, E. & Kudlow, J. E. Stage-sensitive blockade of pituitary somatomammotrope development by targeted expression of a dominant negative epidermal growth factor receptor in transgenic mice. Mol. Endocrinol. 15, 600–613 (2001).

Article CAS PubMed Google Scholar

Jackson, E. L. et al. Analysis of lung tumor initiation and progression using conditional expression of oncogenic K-ras. Genes Dev. 15, 3243–3248 (2001).

Article CAS PubMed PubMed Central Google Scholar

Beronja, S., Livshits, G., Williams, S. & Fuchs, E. Rapid functional dissection of genetic networks via tissue-specific transduction and RNAi in mouse embryos. Nat. Med. 16, 821–827 (2010).

Article CAS PubMed PubMed Central Google Scholar

Ying, Z., Sandoval, M. & Beronja, S. Oncogenic activation of PI3K induces progenitor cell differentiation to suppress epidermal growth. Nat. Cell Biol. 20, 1256–1266 (2018).

Article CAS PubMed PubMed Central Google Scholar

Beronja, S. et al. RNAi screens in mice identify physiological regulators of oncogenic growth. Nature 501, 185–190 (2013).

Article ADS CAS PubMed PubMed Central Google Scholar

Pineda, C. M. et al. Intravital imaging of hair follicle regeneration in the mouse. Nat. Protoc. 10, 1116–1130 (2015).

Article CAS PubMed PubMed Central Google Scholar

Young, M. D. & Behjati, S. SoupX removes ambient RNA contamination from droplet-based single-cell RNA sequencing data. Gigascience 9, giaa151 (2020).

Article PubMed PubMed Central Google Scholar

Satija, R., Farrell, J. A., Gennert, D., Schier, A. F. & Regev, A. Spatial reconstruction of single-cell gene expression data. Nat. Biotechnol. 33, 495–502 (2015).

Article CAS PubMed PubMed Central Google Scholar

McGinnis, C. S., Murrow, L. M. & Gartner, Z. J. DoubletFinder: doublet detection in single-cell RNA sequencing data using artificial nearest neighbors. Cell Syst. 8, 329–337.e324 (2019).

Article CAS PubMed PubMed Central Google Scholar

Wolf, F. A., Angerer, P. & Theis, F. J. SCANPY: large-scale single-cell gene expression data analysis. Genome Biol. 19, 15 (2018).

Article PubMed PubMed Central Google Scholar

Tirosh, I. et al. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science 352, 189–196 (2016).

Article ADS CAS PubMed PubMed Central Google Scholar

Polanski, K. et al. BBKNN: fast batch alignment of single cell transcriptomes. Bioinformatics 36, 964–965 (2020).

Article CAS PubMed Google Scholar

Pedregosa, F. et al. Scikit-learn: Machine Learning in Python. J. Mach. Learn. Res. 12, 2825–2830 (2011).

MathSciNet MATH Google Scholar

Dann, E., Henderson, N. C., Teichmann, S. A., Morgan, M. D. & Marioni, J. C. Differential abundance testing on single-cell data using k-nearest neighbor graphs. Nat. Biotechnol. 40, 245–253 (2022).

Article CAS PubMed Google Scholar

Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl Acad. Sci. USA 102, 15545–15550 (2005).

Article ADS CAS PubMed PubMed Central Google Scholar

Mootha, V. K. et al. PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat. Genet. 34, 267–273 (2003).

Article CAS PubMed Google Scholar

Chen, E. Y. et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics 14, 128 (2013).

Article PubMed PubMed Central Google Scholar

Kuleshov, M. V. et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 44, W90–97 (2016).

Article CAS PubMed PubMed Central Google Scholar

Xie, Z. et al. Gene set knowledge discovery with Enrichr. Curr. Protoc. 1, e90 (2021).

Article CAS PubMed PubMed Central Google Scholar

Korosec, A. et al. Lineage identity and location within the dermis determine the function of papillary and reticular fibroblasts in human skin. J. Invest. Dermatol. 139, 342–351 (2019).

Article CAS PubMed Google Scholar

Schneider, M. R. & Wolf, E. The epidermal growth factor receptor ligands at a glance. J. Cell. Physiol. 218, 460–466 (2009).

Article CAS PubMed Google Scholar

Liu, Y. et al. Single-cell profiling reveals divergent, globally patterned immune responses in murine skin inflammation. iScience 23, 101582 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Download references

The authors thank all members of the Greco and Kasper laboratories, B. Slobodan, S. Regot, K. Politi, B. Hu, L. Gonzalez and A. Andersen for critical feedback on the manuscript; M. Mandar for the LSL-KrasG12D/+ mice; B. De Kumar and G. Wang for their support in scRNA-seq experiments and data analyses; and W. Damsky for the acquisition of H&E staining slides. Research reported in this publication was supported by an 55108527 HHMI Scholar award, the National Institute Of Arthritis and Musculoskeletal and Skin Diseases of the National Institutes of Health under award numbers R01AR063663 and R01AR072668, and the National Institute on Aging of the National Institutes of Health under award number DP1AG066590. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. S.G. was supported by the Human Frontiers Science Program Long-term Postdoctoral fellowship (LT000051_2017-L). M.K. and K.A. were supported by grants from Cancerfonden (CAN 21 1821 Pj), Vetenskapsrådet (VR2018-02963), Karolinska Institutet (2-2111/2019; 2016-00206) and StratRegen (SFO at Karolinska Institutet).

Open access funding provided by Karolinska Institute.

These authors contributed equally: Karl Annusver, Nur-Taz Rahman

Department of Genetics, Yale School of Medicine, New Haven, CT, USA

Sara Gallini, David G. Gonzalez, Sangwon Yun, Catherine Matte-Martone, Tianchi Xin, Elizabeth Lathrop & Valentina Greco

Department of Cell and Molecular Biology, Karolinska Institutet, Stockholm, Sweden

Karl Annusver & Maria Kasper

Bioinformatics Support Program, Cushing/Whitney Medical Library, Yale School of Medicine, New Haven, CT, USA

Nur-Taz Rahman

Dermatologic Surgery, Yale School of Medicine, New Haven, CT, USA

Kathleen C. Suozzi

Departments of Cell Biology and Dermatology, Yale Stem Cell Center, Yale Cancer Center, Yale School of Medicine, New Haven, CT, USA

Valentina Greco

You can also search for this author in PubMed Google Scholar

S.G. and V.G. designed experiments and wrote the manuscript. S.G. performed two-photon imaging, epidermal preparation immunostaining, mouse genetics, molecular assays, image analyses, designed the scRNA-seq experiment and assisted with analyses. D.G.G. assisted with IMARIS and MatLab analyses. N.-T.R., K.A. and M.K. performed scRNA-seq analyses and interpreted the data. K.A. and M.K. provided continuous feedback on the manuscript and assisted in manuscript writing. S.Y. assisted with IMARIS data analysis, LV-CreER in utero injection and provided feedback on the manuscript. C.M.-M. assisted with IMARIS data analysis, epidermal preparation immunostaining and provided feedback on the manuscript. E.L. and T.X. assisted with whole-mount tissue and OCT immunostaining. K.C.S. assisted with clinical diagnosis of histological images and provided feedback on the manuscript.

Correspondence to Maria Kasper or Valentina Greco.

The authors declare no competing interests.

Nature thanks James DeGregori, Kim Jensen, Diether Lambrechts, Carien Niessen and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer review reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

a) Skin tissue architecture cartoon schematic. b, c, d) H&E staining of uninjured and injured ear skin within 2 weeks PWI with magnified insets. n = 5 wild-type, n = 5 HrasG12V/+-max and n = 4 HrasG12V/+-mosaic. b) Wild-type tissue exhibits normal histology (left) and injured tissues (right) demonstrates mixed inflammatory infiltrate and increased fibroblast number consistent with early scar tissue. c) HrasG12V/+-mosaic uninjured tissue (left) shows normal histology. Injured ear (right) exhibits dermal fibrosis consistent with normal scar tissue. d) HrasG12V/+-max uninjured tissue (left) with normal epithelial thickness and differentiation with mild orthokeratotic hyperkeratosis (stratum corneum thickening with normal basket-weave appearance). In the HrasG12V/+-max injured ear (right) adjacent to the wound there is notable acanthosis and hypergranulosis of the epithelium with compact orthokeratotic hyperkeratosis. Cytologic atypia is present. e) Aberrant growth in the injured-HrasG12V/+-max model in contrast to the normal epithelium in the injured-HrasG12V/+-mosaic model (day-14 PWI). f) (right) Heat maps of epithelial thickness in the top-down (x-y) view of representative two-photon images adjacent to injury in the back skin of HrasG12V/+-mosaic and HrasG12V/+-max identifying aberrant growth. Scale bar, 100 μm. (left) Quantification of average epithelial thickness 14 days PWI at different distances from wound edge in HrasG12V/+-mosaic and HrasG12V/+-max. n = 3 mice (mean±s.d.). g, h) H&E staining of injured back skin within two weeks PWI. n = 3 mice. g) HrasG12V/+-mosaic injured back shows skin dermal fibrosis with vertically oriented blood vessels consistent with scar. The overlying epidermis is normal thickness with orthokeratosis. h) HrasG12V/+-max injured back skin with dermal scar and the overlying epidermis is characterized by hyperkeratosis and papillomatosis. i) Aberrant growth around the wound of the HrasG12V/+-max model in contrast to the normal epithelium of the HrasG12V/+-mosaic model.

Source Data

a) Two-photon revisit images of the same area of the basal stem cell layer of the epidermis in wild-type-mosaic and HrasG12V/+-mosaic at 3 and 14 days PWI. Dashed lines, wound edge; white squares, three representative regions tracked and quantified to measure the percent tdTomato+ area. The hair follicle pattern was used to revisit the exact same area of the skin (*=hair canals). Scale bar, 100 μm. b) Quantification of the initial percent tdTomato+ area 6 days post-tamoxifen injection in the uninjured condition. n = 3 mice. c) Representative two-photon revisit images of the same tdTomato+ clone revisited in the basal stem cell layer of the epidermis infected with lentiviral(LV)-CreER. d) Quantification of the relative ratio of the tdTomato+ area at different time points compared to the initial clone size at 6 days post-tamoxifen injection. Relative ratio was used to consider the different initial sizes of the analyzed clones. The initial clone size was between 2 to ~50 cells. n = 3 mice. e) Representative two-photon revisit images of the same tdTomato+ clone in the basal stem cell layer of the epidermis infected with LV-CreER. f) Same quantification as in d). n = 3 mice. (d, f) Statistics: Unpaired, two-tailed t-test between wild-type and mutant mice at different time points in uninjured and injured conditions. Exact p-value reported on the figure. ns indicates not statistically significant. At least three independent tdTomato+ clones were analysed for each mouse. Data are represented as means and standard deviations. Scale bar, 20 μm.

Source Data

a, b) Representative two-photon images of mitotic figures (white circles). c) Quantification of mitotic figures in tdTomato+ and tdTomato− areas in injured (I; n = 4 mice) and uninjured (U; n = 3 mice) wild-type and HrasG12V/+-mosaic. d) Representative two-photon images of apoptotic bodies (white circles) and quantification in tdTomato+ and tdTomato− areas in injured (I) and uninjured (U) wild-type and HrasG12V/+-mosaic. n = 4 mice. e) Representative two-photon images of the epidermal preparation immunofluorescence for active-Caspase-3 (white circle) and quantification of active-Caspase3+ cells in tdTomato+ and tdTomato− areas in injured (I) and uninjured (U) wild-type and HrasG12V/+-mosaic. n = 3 mice. f) Quantification of apoptotic bodies over time in tdTomato+ and tdTomato− areas in uninjured ears. Statistic between wild-type-mosaic and HrasG12V/+-mosaic in tdTomato− or tdTomato+ areas are represented with green or red dotted lines. g) Same quantification as in f) in injured ears. (f, g) n = 4 mice. h) β-galactosidase activity assay of epidermal preparations of uninjured and injured wild-type and HrasG12V/+-mosaic and positive controls: pancreas, kidney and provided by the kit. n = 3 mice for epidermal preparations. i) Quantification of cells expressing Keratin10 in tdTomato+ and tdTomato− areas in injured (I) wild-type and HrasG12V/+-mosaic. n = 3. Statistics: Paired, two-tailed t-test, comparing tdTomato+ and tdTomato− areas in the same group of mice. Unpaired, two-tailed t-test comparing tdTomato+ and tdTomato− populations in different groups of mice. Exact p-value reported on the figure. ns indicates not statistically significant. At least three independent areas of approximately 300 μm2 were analysed for each mouse (see Methods). Data are represented as means and standard deviations. Scale bar, 20 μm.

Source Data

a) Violin plots showing quality control metrics for each sample. All the box plots within violin plots denote the 25%, 50% and 75% quartiles with whiskers depicting the minima and maxima of the data, excluding outliers that are beyond 1.5x interquartile range. b) UMAP displaying the main cell populations of the integrated dataset (left), and dot-plot showing characteristic marker-gene expression (right). c) UMAPs showing the distribution of interfollicular epidermal (IFE) keratinocytes of the different conditions, coloured according to biological replicates. Grey cells denote all keratinocytes. d) Abstracted graph of neighbourhoods superimposed on IFE UMAP, showing differential abundance testing results from MiloR68. Node sizes represent the size of the neighbourhood and edges indicate number of cells common between neighbourhoods. Neighbourhoods displaying significant differential abundance are coloured according to the log-foldchange of the differential abundance testing. e) Micrographs showing Keratin6a protein expression in uninjured and injured ear tissue. W – wound, scale bar, 500 μm. n = 3 mice. f, g, h) UMAPs showing cell classification into basal non-committed, basal committed, and delaminated populations (f), cell cycle phase (g) and cell classification based on tdTomato expression (h). i, j) Bar plots showing top GO terms (upper) for differentially up-regulated genes comparing HrasG12V/+ and wild-type-mosaic cells in uninjured (i) or injured (j) condition. Respective differential gene expression analysis is shown on volcano plots (lower) with log2 fold-changes on x-axis and -log10(adjusted p-values) on y-axis. Genes were considered differentially expressed (orange dots) when they were expressed in >25% of cells in each of the compared biological replicate, had absolute log-foldchange > 0.5 and adjusted p-value < 0.05 (Wilcoxon rank-sum test with Benjamini-Hochberg correction). a-d, f-j) n = 12 independently sequenced mice (3 mice per condition and genotype).

Source Data

a) Quantification of the initial percent tdTomato+ area in the first revisit of KrasG12D/+-mosaic in uninjured (n = 4 mice) and injured (n = 5 mice) conditions and KrasG12D/+-max in injured (n = 3 mice) condition. At least three independent areas of approximately 300 μm2 were analysed for each mouse (Methods). Data are represented as means and standard deviations. b, c) Histopathologic examination via H&E staining of uninjured and injured ear skin within 2 weeks PWI. n = 4 KrasG12D/+-mosaic, n = 5 KrasG12D/+-max. b) KrasG12D/+-mosaic (left) in uninjured ear skin with normal epithelium, dermis and cartilage and injured ear (right) showing mixed inflammatory cell infiltrate and increased number of fibroblasts consistent with early scar tissue. c) KrasG12D/+-max uninjured ear (right) showing normal epithelial thickness and architecture. This is in comparison to the injured ear (left) showing significant expansion of the epithelial layer with hypergranulosis, focal parakeratosis alternative to the compact hyperkeratin. Focal cytologic atypia is present. Scale bars indicated on the figure. Magnified insets of the epidermis in the lower left corner of each image. d) Macroscopic image of the aberrant growth around the wound of the KrasG12D/+-max model in contrast to the normal epithelium of the KrasG12V/+-mosaic model.

Source Data

a) Violin plots showing quality control metrics for each individual sample. b) Unsupervised clustering of sampled fibroblasts74. c) Marker gene expression in fibroblast subclusters overlaid on UMAP. d) Distribution of cells from different mouse models among fibroblast clusters (n = 3 mice for each condition). Grey dots denote all cells. HrasG12V/+ and wild-type models have a similar distribution of cells in four clusters, depending on the tamoxifen treatment and the presence or not of the wound. e) UMAP showing cell-cycle classification (left) and injury response status (right). f) Dot-plot showing marker gene expression for the fibroblast subclusters. g) Volcano plot of differential gene expression profiles between fibroblasts from homeostatic and injured conditions showing the magnitude on the x-axis (Log2 fold change) and significance on the y-axis (-Log10 adjusted p-value, Wilcoxon rank-sum test with Benjamini-Hochberg correction). Red dots mark the 50 highest differentially expressed genes. Red names highlight the highest differentially expressed growth factors that affect epithelial cell behaviors. Blue dots with blue names represent other growth factors that take part in injury-repair and black dots with black names indicate genes involved in cell proliferation. h) Violin plots of EGFR ligands that affect epithelial cell behaviors75. Analysis is based on all mouse models combined. i) Violin plots that compare homeostatic and injury-responsive cells for the expression of EGFR ligands with a significantly different expression in (h). n = 3 mice per group. (a, h, i) Internal box plots denote the 25%, 50% and 75% quartiles with whiskers depicting the minima and maxima of the data, excluding outliers that are beyond 1.5x interquartile range. Statistics: two-tailed t-test comparing the averages of biological replicates according to conditions. (a-h) n = 24 independently sequenced mice (3 mice per condition and genotype, note that some samples did not contain fibroblasts; ‘mosaic samples’ are the same as in Extended Data Fig. 4).

a) Violin plots showing quality control metrics for each individual sample. b) Annotation of cluster identities based on marker gene expression76. c) Unsupervised clustering of sampled immune cells (see Methods)76. d) Distribution of cells from different mouse models among the identified immune cell clusters (n = 3 mice for each condition). e) UMAP representations of cell cycle classification (top) and injury status of the dataset (bottom). f) Violin plots showing immune cell EGFR ligand expression for all datasets combined (n = 24 mice, top) and comparison between wild-type-mosaic and HrasG12V/+-mosaic datasets (n = 6 mice per group, bottom). g) Violin plots showing epidermal (i.e.: IFE) EGFR ligand expression for all datasets combined (n = 12 mice, top) and comparison between wild-type-mosaic and HrasG12V/+-mosaic datasets (n = 6 mice per group, bottom). (a, f, g) Internal box plots denote the 25%, 50% and 75% quartiles with whiskers depicting the minima and maxima of the data, excluding outliers that are beyond 1.5x interquartile range. (f, g) Statistics: two-tailed t-test comparing the means of biological replicates according to conditions. (a-f) n = 24 independently sequenced mice (3 mice per condition and genotype; ‘mosaic samples’ are the same as in Extended Data Figs. 4, 6).

a) Schematic representation of the EGFR signaling pathway. HrasG12V is constitutively active and therefore less dependent on upstream activation of EGFR. b) Western blot analysis of p-EGFR(phospho-Tyr1068) normalized on total-EGFR (n = 6 mice) and of p-ERK1/2(phospho-Thr202/Tyr204) and p-AKT(phospho-Ser473) normalized on total-ERK1/2 and total-AKT (n = 3 mice). Unpaired, two-tailed t-test. c) Mitotic figure quantification in tdTomato+ and tdTomato− areas in HrasG12V/+-mosaic without or with EGFR-Dominant Negative (DN) expression (EGFR-WT or EGFR-DN). n = 3 mice. Unpaired or Paired two-tailed t-test for comparison between different or the same groups of mice. d) Revisit images of the same area of the basal stem cell layer. (left) The increase of tdTomato+ area. n = 3 mice. Unpaired, two-tailed t-test. e) Initial percent tdTomato+ area in the first revisit of injured-HrasG12V/+-mosaic treated with vehicle (n = 3 mice) or Gefitinib (n = 4 mice). f) Epidermal preparation immunofluorescence for phospho-Histone3 in wild-type (left) and constitutive-p21null (right) at postnatal-day-25. g) Phospho-Histone3+ cell quantification in wild-type, constitutive-p21null with or without LSL-tdTomato at postnatal-day-25 (n = 3 mice). Unpaired, two-tailed t-test. h) Whole mount immunofluorescence for phospho-Histone3 in wild-type (left) and constitutive-p21null (right) at postnatal-day-25. Phospho-Histone3+ dermal cell quantification (n = 4 mice). Unpaired, two-tailed t-test. i) Initial percent tdTomato+ area in the first revisit of uninjured wild-type-mosaic (n = 3 mice), HrasG12V/+-mosaic (n = 3 mice) and constitutive-p21null-HrasG12V/+-mosaic (n = 4 mice). j) Keratin10+ cell quantification in tdTomato+ and tdTomato− areas in HrasG12V/+-mosaic and constitutive-p21null-HrasG12V/+-mosaic in uninjured(U) ears (n = 3 mice). Unpaired or Paired two-tailed t-test for comparison between different or the same groups of mice. Exact p-values reported on the figure. ns indicates not statistically significant. At least three independent areas of approximately 300 μm2 were analysed for each mouse (Methods). Data are represented as means and standard deviations. Scale bar, 20 μm.

Source Data

a) Representative two-photon images of the epidermal preparation immunofluorescence for p21 (white) in constitutive-p21null (left), wild-type-mosaic (center) and HrasG12V/+-mosaic (right). The basal cells expressing K14H2B-GFP are marked in magenta. Representative images show both the overlay (top) and separate channels for K14H2B-GFP expression (middle) and p21 immunolabeling (bottom). *=mitotic cells negative for p21 (background signal). b) Quantification of p21 immunolabeling in the nucleus of the basal cells in wild-type-mosaic and HrasG12V/+-mosaic (n = 3 mice). Average of approximately 3,000 cells were analysed for each mouse. Unpaired, two-tailed t-test for the nuclear intensity of p21. Data are represented as means and standard deviations. c) p21 level quantification in the basal cells by using scRNA-seq datasets of wild-type-mosaic and HrasG12V/+-mosaic with the selection of only tdTomato+ cells. d) Same quantification as in c) but with the selection of non-cycling and tdTomato+ cells. (c, d) n = 3 mice. Statistics on scRNA-seq analysis of p21 expression is based on the non-parametric Wilcoxon rank sum test performed by the Seurat package62. Adjusted p-value is based on Bonferroni correction using all features in the dataset.

Source Data

a) Western blot analysis of p-ERK1/2(phospho-Thr202/Tyr204) normalized on total-ERK1/2. b) Western blot analysis of p-AKT(phospho-Ser473) normalized on total-AKT. (a, b) wild-type (n = 4 mice), HrasG12V/+-mosaic (n = 3 mice), HrasG12V/+-max (n=3 mice) and constitutive-p21null (n = 3 mice). Unpaired, two-tailed t-test. Data are represented as means and standard deviations. c) Confocal representative images of the epidermal preparation immunofluorescence for p-ERK1/2 and phalloidin in HrasG12V/+-mosaic. (left) Insets of p-ERK1/2+ cells: a) cells in interphase b) mitotic cell and c) footprint of a cell that is departing from the basal layer. d) Quantification of p-ERK1/2+ cells in tdTomato+ and tdTomato− areas in HrasG12V/+-mosaic in injured/3 days PWI(I) and uninjured(U) ears. At least four independent 150 μm2 areas were analysed for each mouse. (c, d) n = 3 mice. Unpaired or Paired two-tailed t-test for comparison between different or the same groups of mice. ns indicates not statistically significant. Data are represented as means and standard deviations. Scale bar, 20 μm. (a, b, c, d) Exact p-values reported on the figure. e) Final model. In the uninjured mosaic skin epidermis, Ras cells integrate and expand, outcompeting wild-type neighbors. During injury-repair of mosaic skin, the competitive advantage of Ras cells is suppressed, and oncogenic growths do not develop. EGFR signaling pathway is crucial for this selective increase of wild-type cell divisions that prevents Ras cell expansion during injury-repair. Inducing proliferation, via constitutive p21 loss, mimics the injury condition in uninjured skin, counteracting the competitive advantage of Ras cells. Our data support a model whereby injury-repair and p21 loss increase the activity of ERK1/2, a Ras downstream pathway that controls cell proliferation. This leads to comparable ERK1/2 levels between wild-type and Ras cells, resulting in an increase of dividing wild-type cells that effectively prevent Ras mutant cell expansion.

Source Data

Uncropped and unprocessed immunoblots of Fig. 5b (n = 3 mice for each condition), Extended Data Fig. 8b (n = 6 mice p-EGFR/EGFR, n = 3 mice p-ERK1/2/ERK1/2 and p-AKT/AKT) and Extended Data Fig. 10 a,b (n = 4 mice wild-type, n = 3 mice HrasG12V/+ mosaic, n = 3 mice HrasG12V/+ max). For statistical and densitometry analyses see figure legends and data sources. Membranes were cut around the molecular weight of interest before blotting. Molecular markers were acquired with colorimetric setting. Proteins were detected by western blotting, using the indicated antibodies and the red rectangles show the cropping location.

Representative z-stack (every 1 µm) of the ear epidermis at 14 days. PWI of a wild-type mosaic mouse (wound on the left of the video). All epithelial nuclei are in green (Krt14-H2B–GFP) and recombined cells are in red (LSL-tdTomato). Optical sections start from the differentiated cornified layer of the skin epithelium and proceed downward into the dermis defined by the second harmonic collagen signal (cyan).

Representative z-stack (every 1.2 µm) of the ear epidermis at 14 days PWI of a HrasG12V/+ mosaic mouse (wound on the left of the video). All epithelial nuclei are in green (Krt14-H2B–GFP) and recombined cells are in red (LSL-tdTomato). Optical sections start from the differentiated cornified layer of the skin epithelium and proceed downward into the dermis defined by the second harmonic collagen signal (cyan).

Representative z-stack (every 1.3 µm) of the ear epidermis 14 days PWI of a HrasG12V/+ max mouse (wound on the left of the video). All epithelial nuclei are in green (Krt14-H2B–GFP) and recombined cells are in red (LSL-tdTomato). Optical sections start from the differentiated cornified layer of the skin epithelium and proceed downward into the dermis defined by the second harmonic collagen signal (cyan). The aberrant growth present close to the wound edge prevents us from resolving the dermis below due to the thickness of the tissue and hyper-keratinization.

Comparison of representative z-stacks of the ear epidermis at 14 days PWI of wild-type mosaic, HrasG12V/+ mosaic and HrasG12V/+ max mouse models (wound on the left of the videos). The epidermis is depicted in white (combined Krt14-H2B–GFP and LSL-tdTomato expression). Optical sections start from the differentiated cornified layer of the skin epithelium and proceed downward into the dermis defined by the second harmonic collagen signal (cyan). The comparison of the optical sections between the three models shows the different thickness of the epidermis around the wound and exemplifies the presence of aberrant growth in the HrasG12V/+ max mouse.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Gallini, S., Annusver, K., Rahman, NT. et al. Injury prevents Ras mutant cell expansion in mosaic skin. Nature 619, 167–175 (2023). https://doi.org/10.1038/s41586-023-06198-y

Download citation

Received: 05 March 2021

Accepted: 11 May 2023

Published: 21 June 2023

Issue Date: 06 July 2023

DOI: https://doi.org/10.1038/s41586-023-06198-y

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

A History of Fatigue, Reviewed

How To Strain Pasta When You Don't Have A Colander

News

Injury prevents Ras mutant cell expansion in mosaic skin