The increasing implication of lymphocytes in general physiology and immune surveillance outside of infection poses the question of how their antigen receptors might be involved. Here, we show that macromolecular aggregates of intraepidermal γδ T cell antigen receptors (TCRs) in the mouse skin aligned with and depended on Skint1, a butyrophilin-like (BTNL) protein expressed by differentiated keratinocytes (KCs) at steady state. Interruption of TCR-mediated ‘normality sensing’ had no impact on γδ T cell numbers but altered their signature phenotype, while the epidermal barrier function was compromised. In addition to the regulation of steady-state physiology, normality sensing licensed intraepidermal T cells to respond rapidly to subsequent tissue perturbation by using innate tumor necrosis factor (TNF) superfamily receptors. Thus, interfering with Skint1-dependent interactions between local γδ T cells and KCs at steady state increased the susceptibility to ultraviolet B radiation (UVR)-induced DNA damage and inflammation, two cancer-disposing factors.
The steady-state immune system was long thought to be resting, awakened only by infection. However, this perspective has been transformed by examples of immune responses to non-microbial challenges and evidence that the steady-state immune system helps maintain physiological processes from metabolism to memory and barrier integrity1,2,3,4. Many such interactions happen locally, involving immune cells permanently associated with non-lymphoid tissues5,6,7. While macrophages are well-established regulators of tissue status, an increasing implication of lymphocytes in this process raises the question of whether antigen receptors are used to sense tissue perturbation, akin to recognition of foreign antigens.
To address this, we have focused on intraepithelial lymphocytes (IELs), which compose one of the largest T cell compartments and are variably conserved in multiple tissues from agnathans to humans8,9. IELs are commonly implicated in tissue stress surveillance but with little mechanistic understanding10,11,12. IELs are disproportionately enriched in γδ T cells, which express site-specific TCRs, for example, Vγ5Vδ1 in mouse epidermis, Vγ7Vδn in mouse intestine and Vγ4Vδ1 in human intestine9,13,14,15. Such selective repertoire focusing reflects the developmental dependence of IELs on cognate tissue-specific, epithelial-expressed BTNL heteromers. Thus, Skint1 and Skint2, two PD-L1-like members of the BTNL family, select murine Vγ5Vδ1+ IELs, Btnl1 and Btnl6 select murine Vγ7+ intestinal IELs, and BTNL3 and BTNL8 selectively regulate human colonic Vγ4+ IELs15,16,17,18,19. Indeed, BTNL3 is a bona fide TCRγδ ligand20,21.
Whereas Skint1 and Skint2 mediate Vγ5Vδ1+ precursor development in the thymus, they are also expressed by suprabasal KCs, among which mature Vγ5Vδ1+ IELs permanently reside17. How they might continue to influence mature IEL function is unknown. In that regard, focal clusters of Vγ5Vδ1+ TCRs constitutively engage KCs at steady state and may engage in signaling22, but neither the epithelial determinant(s) nor the biological implications of these engagements is understood.
Here, we focused on Vγ5Vδ1+ IELs, also termed dendritic epidermal T cells (DETCs), to investigate how tissue T cells might use TCRs to distinguish normal self from stressed self. We show that constitutive DETC TCR-containing foci22 align with and depend upon sustained epithelial expression of Skint1. Although these interactions did not maintain DETC numbers, they preserved DETC and KC gene expression patterns and maintained steady-state epithelial barrier function. Additionally, these interactions licensed DETCs to utilize costimulatory receptors to make rapid responses to UVR and chemical irritation. Thus, interruption of steady-state, TCR-mediated DETC–KC interactions increased epidermal susceptibility to sustained DNA damage and inflammation, two factors that jointly dispose to cancer23. As Skint1-dependent DETC–KC interactions occur at steady state, this mode of tissue immunosurveillance can be viewed as ‘normality sensing’. This offers a revised perspective of tissue immunosurveillance and of antigen receptor biology, which conventionally focuses on the recognition of molecules induced by change.
DETCs sense steady-state KCs via Skint1
At steady state, Vγ5Vδ1+ TCRs focus at apical dendrite tips in supramolecular clusters that contact neighboring suprabasal KCs, seemingly inducing constitutive TCR signaling22. TCRs within these clusters colocalized with actin (Fig. 1a), as reported22, and most DETCs in FVB mice displayed one to four clusters (Fig. 1a). For clarity, we displayed F-actin staining masked by DETC bodies, although examination of raw staining confirmed these observations (Extended Data Fig. 1a). We next examined FVB.Tac mice, which express a truncated non-functional Skint1, and Trgv5–/–Trdv4–/– mice. These two strains lack canonical Vγ5Vδ1+ DETCs12,16. Relative to wild-type mice, both strains showed increased frequencies of DETCs that lacked TCR–actin colocalization (Fig. 1a and Extended Data Fig. 1a). Transgenic expression of wild-type Skint1 in FVB.Tac mice returned TCR–actin colocalization to wild-type levels (Fig. 1a). Because DETCs in FVB.Tac and Trgv5–/–Trdv4–/– mice use non-Vγ5Vδ1+ TCRs, these results confirm reports that juxtaepithelial TCR foci require the Vγ5Vδ1 TCR22.
To identify the epithelial determinant of Vγ5Vδ1-dependent KC–DETC interactions, we focused on Skint1, given its critical and selective contribution to fetal thymic Vγ5Vδ1+ DETC development and its expression by differentiated suprabasal KCs17,24. To limit Skint1 deficiency to the epidermis, we adoptively transferred wild-type embryonic day 16 (E16)–E18 fetal thymocytes, which include Skint1-selected CD45RBhi DETC progenitors16, into neonatal Tcrd–/– or FVB.Tac Tcrd–/– mice. Epidermal Vγ5+ DETCs formed in both hosts (Extended Data Fig. 1b), suggesting that epidermal Skint1 was not strictly necessary for the epidermal seeding of DETCs. While DETCs in Tcrd–/– hosts displayed bright colocalized TCR–actin foci resembling wild-type mice, markedly fewer foci were detected in DETCs in FVB.Tac Tcrd–/– mice (Extended Data Fig. 1b), strongly suggesting that Skint1 was required to establish and/or maintain juxtaepithelial TCR focusing.
To investigate whether Skint1 protein acutely maintained TCR–KC interactions, we performed an i.d. administration of 2G2, a monoclonal antibody specific for Skint1 (ref. 19), to wild-type FVB mice. Within 6 h after injection, the frequency of TCR–actin foci was significantly reduced compared to mice administered an isotype control antibody (Fig. 1b). A similar impact was seen in C57BL/6 wild-type mice (Extended Data Fig. 1c). We next applied iSIM to steady-state epidermal sheets from wild-type mice using Skint1 antibody staining. Despite Skint1 expression being restricted to KCs17, we only detected Skint1 protein in clusters nucleated by the TCR (Fig. 1c and Extended Data Fig. 1d). Compared to wild-type mice, TCR-proximal Skint1 antibody staining was significantly reduced in FVB.Tac and Tcrd–/– mice (Extended Data Fig. 1e). While truncated Skint1 in FVB.Tac mice could affect antibody staining, Tcrd–/– mice expressed Skint1 mRNA at relatively normal levels in the epidermis (Extended Data Fig. 1f), indicating that wild-type Skint1 protein was not nucleated by TCRαβ+ cells in the Tcrd–/– epidermis25. Thus, Skint1 maintained steady-state TCR-mediated interactions of DETCs with KCs, which in turn shaped the epidermal localization of Skint1.
Skint1 promotes DETC maintenance of epidermal integrity
To analyze the impact of sustained Skint1 disruption in situ, we performed an i.d. injection of Skint1 antibody or isotype control weekly into the ears of FVB mice for 10 weeks. At this time point, mice displayed comparable numbers and frequencies of Vγ5Vδ1+ DETCs and major histocompatibility complex class II+ (MHCII+) Langerhans cells (LCs) (Fig. 2a and Extended Data Fig. 2a–c), indicating that long-term DETC population maintenance was Skint1 independent. However, relative to isotype-treated or contralateral ears, Skint1 antibody-treated ears became moderately inflamed (Fig. 2b,c), and CD5+ αβ T cells infiltrated the epidermis (Fig. 2d). Similar increased ear thickness, albeit milder, was observed in control-treated Tcrd–/– mice (Extended Data Fig. 2d), as previously reported26. Sustained Skint1 antibody treatment did not further increase ear thickness in Tcrd–/– mice (Extended Data Fig. 2d,e), indicating that Skint1 antibody treatment primarily interfered with the Skint1–TCRγδ interaction.
We additionally observed increased transepidermal water loss (TEWL) at week 5 of Skint1 antibody treatment relative to isotype treatment (Fig. 2e) or contralateral ears (Extended Data Fig. 2f). Increased baseline TEWL measurements were also detected in isotype-treated Tcrd–/– mice, as reported in ref. 4, and those were also unaffected by Skint1 antibody treatment (Fig. 2e). C57BL/6 mice lacking Skint1 or all Skint genes also had elevated baseline TEWL compared to littermate controls (Extended Data Fig. 2g). Compared to isotype-treated mice, Skint1 antibody-treated mice consistently displayed decreased epidermal expression of the DETC-derived barrier maintenance cytokine Il13 (ref. 27) and of KC genes involved in epidermal barrier formation (Tgm1 and Lor; Fig. 2f,g), as did baseline Tcrd–/– mice relative to wild-type controls (Extended Data Fig. 2h,i). KCs also displayed increased apoptosis following Skint1 antibody treatment (Fig. 2h). While keratin gene expression was unchanged, indicating normal KC differentiation, epidermal Skint1 and Skint2 expression was upregulated following Skint1 antibody treatment (Extended Data Fig. 2j), perhaps reflecting a compensatory response to Skint1–TCR disruption. Although LC abundance was unchanged (Fig. 2a and Extended Data Fig. 2a,c), Skint1 antibody treatment marginally increased the frequency of activated MHCIIhiCD86+ LCs (Fig. 2i), perhaps reflecting ongoing epidermal inflammation. Most DETCs contact LCs at steady state, and LCs may also contribute to epidermal integrity28,29. However, we did not observe differences in baseline TEWL or barrier gene expression in LC-deficient human Langerin (huLangerin)-DTA mice relative to wild-type mice (Extended Data Fig. 2k,l). Thus, the functional contributions of DETCs to effective epidermal integrity were compromised by Skint1 blockade and were LC independent.
Skint1 sensing maintains DETC immunosurveillance signatures
To investigate whether the contribution of Skint1 to skin physiology reflected downstream consequences of a primary impact on DETCs, we conducted single-cell RNA sequencing of CD45+MHCII−EpCAM– epidermal lymphocytes from wild-type FVB mice 7 d after Skint1 antibody or isotype treatment. Non-DETC populations, including αβ T cells and contaminating dermal γδ T cells, were excluded. Four DETC clusters (cluster 0–cluster 3) were identified by uniform manifold projection (UMAP) of the integrated dataset (Fig. 3a), indicating that DETCs were heterogeneous, despite their near uniform anatomical localization and TCR expression. Cluster 3 appeared segregated in two dimensions (Fig. 3a), an artifact of dimensionality reduction that could not be resolved by adjusting algorithm parameters.
Cluster 0 was enriched in transcripts related to AP1-mediated cell regulation (Fos, Junb and Jund) and the epithelia-tropic cytokines Areg (amphiregulin) and Il13 (Fig. 3b and Extended Data Fig. 3a–c), implying a potential to regulate KCs and other immunocytes27,30,31,32. Relative to other clusters, cluster 1 appeared more quiescent, with reduced expression of gene sets associated with T cell activation, immune effector processes and regulation of cell adhesion and apoptosis (Extended Data Fig. 3b,c). Nonetheless, cluster 1 was enriched in transcripts encoding galectins, S100 proteins and the aryl hydrocarbon receptor (Ahr) (Fig. 3b), of which the latter is essential for DETC population self-renewal33. Germane to that, cluster 2 was enriched in transcripts associated with proliferation (Fig. 3b and Extended Data Fig. 3c,d). Cluster 3, while small (Fig. 3a), was enriched in mRNAs ascribed to lymphocyte activation and costimulation, such as Cd82, Ramp1 and Tnfrsf9 (4-1BB), as well as effector molecules Xcl1, Ctsw, Gzmc (granzyme C) and Il13 (Fig. 3b,c and Extended Data Fig. 3b,c,e). 4-1BB is conventionally viewed as a costimulatory receptor expressed by recently activated T cells, but it is also induced by Skint1- and Btnl1-mediated selection of DETCs and gut Vγ7+ IEL progenitors, respectively15,34,35. Indeed, flow cytometry sorting indicated that 4-1BBhi DETCs expressed more Gzmc and Xcl1 than 4-1BBlo DETCs (Extended Data Fig. 3f). In Nur77–GFP reporter mice, which report calcium-dependent signaling following T cell activation36, GFP expression was significantly higher in 4-1BBhi DETCs at steady state (Fig. 3d). These results suggested that cluster 3 might represent a subset of DETCs actively responding to minor skin perturbations, such as grooming or scratching, which occur episodically in mice at steady state.
Comparison of DETCs from isotype- and Skint1 antibody-treated mice indicated that cluster 3 was selectively diminished by Skint1 antibody treatment (Fig. 3e,f), which was reflected by a loss of overall Xcl1 and Gzmc expression (Fig. 3g,h) and validated by qPCR of DETCs sorted from Skint1 antibody- and isotype-treated mice (Extended Data Fig. 3g). Consistently, 4-1BBhi DETCs were significantly reduced by Skint1 antibody treatment (Fig. 3i), and the remaining 4-1BBhi DETCs expressed less Nur77–GFP (Fig. 3d). Thus, sustained activation of a small subset of DETCs (cluster 3) requires Skint1 sensing at steady state.
We also observed broader changes in gene expression in DETCs following Skint1 antibody treatment, regardless of clusters. Ptpn6 (SHP-1 phosphatase) and overall Xcl1 expression decreased, and Rps27, Rps29 and Egr1 expression increased following Skint1 antibody treatment (Fig. 3h and Extended Data Fig. 3h), as validated by qPCR (Extended Data Fig. 3g). This indicated that the intrinsic activation state of most DETCs depended on sustained Skint1-dependent DETC–KC interactions. Peripheral Skint1-dependent regulation of Ptpn6 expression in most DETCs (Extended Data Fig. 3i) was reminiscent of the upregulation of Ptpn6 during Skint1-dependent DETC precursor selection (Supplementary Table 1). Similarly, Skint1 antibody treatment moderately reduced DETC expression of CD122 and CD45RB proteins (Extended Data Fig. 3j), which are also upregulated during BTNL-dependent IEL selection15,16. Skint1 antibody-treated wild-type DETCs resembled non-Vγ5Vδ1+ DETCs in FVB.Tac mice regarding CD45RB expression (Extended Data Fig. 3k). GSEA indicated that the collective gene signature of thymic Skint1-driven DETC progenitor selection35 was downregulated in mature DETC by Skint1 antibody treatment (Fig. 3j and Supplementary Table 1). 4-1BBhi DETC frequency and CD122 expression were unaltered in huLangerin-DTA mice (Extended Data Fig. 3l), suggesting that Skint1 regulation of DETCs was direct and LC independent. Thus, steady-state epidermal Skint1 promoted immunosurveillance by maintaining the predetermined effector phenotype in most DETCs and by underpinning the activation of a small subset of DETCs.
Skint1 sensing drives DETC responses to UVR
To further address the contribution of Skint1 to immunosurveillance, we irradiated the ears of C57BL/6 mice with UVR to induce DNA damage. Within 1 d, many KCs exhibited signs of DNA damage and repair, as measured by γH2A.X foci (Fig. 4a). Simultaneously, DETCs became rounded (Fig. 4a), consistent with reports of their response to diverse skin perturbations, including UVR37,38,39. We also observed a substantial diminution in juxtaepithelial TCR focusing (Fig. 4a) and Skint1 and Skint2 downregulation in KCs relative to unchallenged mice (Extended Data Fig. 4a). Concurrently, LCs and DETCs displayed upregulation of surface activation markers (Extended Data Fig. 4b,c).
RNA-sequencing analysis of flow cytometry-sorted DETCs 24 h after UVR revealed rapid and pleiotropic changes in gene expression compared to non-irradiated controls (Fig. 4b). This included further upregulation of Tnfrsf9 (4-1BB) and Il2rb (CD122), upregulation of Icos and Kit (c-Kit), which encode activating/renewal receptors, and upregulation of inhibitory receptors Tigit, Havcr2 (TIM-3) and Cd274 (PD-L1) (Fig. 4b). Moreover, UVR resulted in upregulation of many genes encoding potential effector molecules, including Cxcl2, Vegfa, Gzmb, Gzmf, Ctsc and Ctsh (Fig. 4b). GSEA revealed an enrichment in a cluster 0 signature in DETCs following UVR (Fig. 4c), indicating that DETCs further upregulate cluster 0-defining genes following UVR activation and/or that DETCs in other clusters assume a cluster 0 phenotype. Notably, DETCs from irradiated FVB.Tac and Trgv5–/–Trdv4–/– mice did not display upregulation of Gzmb, Gzmf or CD69 compared with wild-type FVB and Skint1-transgenic FVB.Tac mice (Fig. 4d and Extended Data Fig. 4d), strongly suggesting that the canonical Vγ5Vδ1 TCR was essential for UVR responsiveness.
To test whether Skint1-dependent TCR–KC interactions established and/or maintained the preparedness of DETCs to respond to UVR, Skint1 antibody was injected i.d. into wild-type C57BL/6 mice 1 d before UVR, and DETCs were sorted and analyzed by qPCR 24 h after UVR. The UVR-induced upregulation of Gzmb, Gzmf and Cxcl2 was significantly reduced in the Skint1 antibody-treated mice, while UVR-induced inhibitory receptor genes, such as Havcr2 and Tigit, remained unaffected, and Il13 was atypically upregulated compared to isotype-treated mice (Fig. 4e and Extended Data Fig. 4e). In addition, the UVR-induced increase in 4-1BB and CD122 expression was significantly inhibited by prior Skint1 antibody administration (Fig. 4f), while LC activation was unperturbed (Extended Data Fig. 4f). In Nur77–GFP reporter mice, DETCs displayed a basal level of GFP, which increased incrementally following UVR (Fig. 4g). Skint1 antibody treatment not only diminished the basal Nur77–GFP signal in DETCs but substantially suppressed the UVR-induced increase in Nur77–GFP expression (Fig. 4h). These observations were reproduced ex vivo in an explant model, in which DETCs in Nur77–GFP ear sheets preincubated with Skint1 antibody showed impaired upregulation of Nur77–GFP expression following in vitro exposure to UVR relative to isotype control (Extended Data Fig. 4g), indicating that Skint1 antibody acted locally. Pretreatment of ex vivo ear sheets with the Src-family kinase inhibitor PP2 produced a similar outcome (Extended Data Fig. 4h). Collectively, these observations indicate that Skint1-dependent interactions licensed DETC stress responsiveness and effector functions.
Skint1 licenses DETC responsiveness to UVR via TNF receptors
TCR signaling amplification is a defining property of costimulatory molecules, including those of the TNF receptor superfamily. To investigate whether steady-state Skint1-dependent interactions licensed DETCs to respond rapidly through co-receptors such as 4-1BB, we performed i.d. injections of an agonistic 4-1BB antibody or isotype control in the ears of Nur77–GFP mice. Akin to UVR, 4-1BB antibody drove a coordinated upregulation of CD69 and Nur77–GFP by a subset of DETCs relative to isotype control (Fig. 5a). This effect was likely direct, as epidermal 4-1BB expression was exclusive to DETCs (Extended Data Fig. 5a,b); however, i.d. injection with Skint1 antibody 24 h before 4-1BB stimulation ablated CD69 and Nur77–GFP upregulation in DETCs (Fig. 5a). Additionally, 4-1BB antibody injection induced expression of DETC effector genes related to cluster 3 (Xcl1 and Gzmc) and UVR responses (Gzmb and Gzmf), which was similarly inhibited by prior Skint1 antibody treatment (Fig. 5b). These observations indicated that steady-state TCR–Skint1 interactions with KCs licensed DETCs to respond to 4-1BB, inducing qualitative changes that partially phenocopied the UVR response.
Therefore, to test whether 4-1BB was an obligate contributor to UVR responses, we performed an i.d. injection of an antibody that blocks 4-1BB ligand (4-1BBL) binding to 4-1BB into the ears of wild-type C57BL/6 mice before UVR treatment. This antibody partially, but significantly, blocked the UVR-induced increase in epidermal expression of Gzmb and Gzmf relative to isotype control (Fig. 5c). Consistent with this, UVR-induced upregulation of Gzmb and Gzmf was also inhibited in 4-1BBL-deficient Tnfsf9–/– mice relative to wild-type controls (Extended Data Fig. 5c).
Similar to 4-1BB, expression of Tnfrsf18 (which encodes GITR) was also induced by Skint1 during thymic DETC progenitor selection (Supplementary Table 2). Unlike 4-1BB, Tnfrsf18 was not enriched in cluster 3 but was expressed by most DETCs (Extended Data Fig. 5d). Additionally, GITR mRNA and protein were expressed by KCs and a small fraction of LCs (Extended Data Fig. 5a,b). Regardless, i.d. injection of an agonistic GITR antibody also upregulated Nur77–GFP and CD69 expression by DETCs, which was almost entirely Skint1 dependent (Fig. 5d). Furthermore, antagonist antibodies that blocked GITR ligand (GITRL) interactions with GITR partially but significantly inhibited the UVR-induced upregulation of epidermal Gzmb and Gzmf (Fig. 5e). Thus, signaling through TNF receptors was one non-redundant component of the rapid DETC response to skin irradiation and was dependent on prior licensing by Skint1–TCR interactions.
Next, we asked whether other conduits of innate T cell activation required licensing via Skint1–TCR interactions. Extracellular ATP is an alarmin implicated in DETC responses to UVR together with interleukin-1 (IL-1)37. Administration (i.d.) of ATP alone induced Nur77–GFP upregulation in DETCs, which was Skint1 dependent (Fig. 5f). Thus, steady-state Skint1 sensing regulated the ability of DETCs to respond to diverse stimuli.
Skint1-licensed immunosurveillance promotes recovery from epithelial damage
To investigate the consequences of Skint1 blockade for host physiology, we performed an i.d. injection of Skint1 antibody or isotype control 24 h before and 48 h after UVR. Mutagenic cyclobutane pyrimidine dimers (CPDs) were rapidly induced by UVR in all mice but were removed less efficiently in Skint1 antibody-treated mice than in isotype control-treated mice 72 h after UVR (Fig. 6a). Likewise, higher numbers of γH2A.X foci, more ear swelling and greater epidermal thickening were detected in Skint1 antibody-treated mice (Fig. 6b–d and Extended Data Fig. 6a). Exacerbated ear swelling was similarly observed after UVR in C57BL/6 mice deficient in γδ T cells, Skint1 and the Skint locus (Extended Data Fig. 6b–d), providing genetic evidence that the pathophysiologic impacts of acute Skint1 blockade reflected consequences of disrupted Skint1-dependent TCR–KC interactions. Finally, UVR-associated defects in barrier function, as measured by toluidine blue permeation40, were also significantly exacerbated in Skint1 antibody-treated mice (Fig. 6e).
Skin is commonly exposed to non-mutagenic chemical irritants. To investigate whether Skint1 licensing applied in contexts other than DNA damage, we examined a model of irritant contact dermatitis by using a single topical dose of the hapten dinitrofluorobenzene (DNFB). DETCs upregulated Gzmb and Gzmf within 48 h after DNFB application, but this was largely inhibited by prior Skint1 antibody treatment (Fig. 7a). Similarly, Skint1 blockade suppressed the DNFB-dependent induction of Nur77–GFP in DETCs (Fig. 7b). Moreover, acute DNFB-induced dermatitis was amplified in mice treated with Skint1 antibody (Fig. 7c). Thus, Skint1 licensed prompt DETC effector responses in settings of mutagenic and non-mutagenic epidermal stress, naturally limiting inflammatory lesions and impairments in functional barrier integrity.
Here we showed that steady-state Vγ5Vδ1+ TCR-mediated interactions between DETCs and KCs involve and depend on Skint1. This crosstalk did not obviously regulate the tissue residency of DETCs nor was it required for DETCs to populate the epidermis. However, it was essential for DETC function, as disruption of Skint1 sensing impeded DETC maintenance of epidermal homeostasis, interfered with steady-state immunosurveillance signatures, and compromised efficient DETC responsiveness to epithelial stresses. Indeed, TCR-dependent, Skint1-dependent interactions with KCs licensed DETCs for innate-like responsiveness through TNF receptors, which limited tissue pathology following insults that included UVR and a chemical irritant.
Because Skint1 is expressed by differentiated suprabasal epithelial cells and is stress-labile, DETC antigen receptors sense normality. This contrasts with the conventional consideration that TCRs engage products of change, such as microbial- or cancer-associated neoantigens and stress antigens upregulated by cellular dysregulation. It is possible that the DETC TCR also transduces signals from an unelucidated stress-induced ligand distinct from Skint1 (ref. 41). Nonetheless, stress sensing by DETCs depended largely on innate receptors 4-1BB and GITR, drawing parallels to NKG2D-driven DETC surveillance12. Critically, these receptors were acquired during thymic Skint1-dependent DETC progenitor selection35 and were sustained in the periphery by ongoing Skint1–TCR-mediated interactions. In this way, among others, steady-state TCR-mediated interactions licensed DETCs to respond rapidly to gross insults.
There are several major consequences of TCR-mediated normality sensing. First, Skint1-dependent, TCR-mediated interactions regulated barrier function and tissue homeostasis, which may occur by regulation of KC-tropic growth factor expression by DETCs. While IL-13 contributes to inflammatory pathologies42, steady-state IL-13 production by DETCs reportedly benefits epidermal integrity27. Although mice lacking Vγ5Vδ1+ DETCs do not suffer life-threatening barrier deficiencies, they do display cutaneous inflammation that could reflect compromised barrier function4,16, as we observed following sustained Skint1 blockade.
Second, Skint1 sustained the signature phenotype of mature DETCs, largely recapitulating its central role in DETC progenitor selection16,17,24,35. This suggests that while the effector potential of DETCs, including predetermined expression of TNFRs and other receptors, is programmed by thymic Skint1 selection, maintaining this profile requires continued peripheral Skint1 sensing. Cluster 3 DETCs, displaying hallmarks of recent activation in unchallenged mice, could represent Skint1-licensed natural immunosurveillance responses to spontaneous epithelial perturbations from scratching/grooming or mutagenesis. Intriguingly, DETC-derived XCL1 reportedly recruits cross-presenting dendritic cells (DCs) to the epidermis43. Thus, XCL1 release from cluster 3 DETCs responding to episodic, low-level KC dysregulation might, in cases of infection or genotoxicity, facilitate neoantigen uptake by DCs, thereby initiating cytotoxic CD8+ T cell immunosurveillance.
Third, Skint1 blockade reduced the capacity of DETCs to sense tissue perturbation through innate receptors. In the context of UVR, the prime etiologic agent of human skin cancer, we observed sustained genotoxicity and inflammatory pathology, which may in part explain the regulation of cutaneous carcinogenesis by DETCs and analogous roles fulfilled by other mouse and human γδ T cell compartments11,44. Additionally, normality sensing primed DETC responses to an irritant and to ATP, and it is conceivable that it might prime for wound healing and tissue regeneration. Indeed, disruption of TCR signaling in mature DETCs impaired wound healing45, a process reportedly modulated by Skint family members other than Skint1 and Skint2 (ref. 46).
Normality sensing may license invariant DETCs to respond through a panoply of innate receptors that collectively offer versatility to detect myriad manifestations of tissue perturbations with which IELs have been associated. Tissue-wide DNA damage, like that induced by UVR, resulted in overt granzyme upregulation, which may facilitate DETC killing of potentially malignant KCs47. Granzyme B also cleaves extracellular targets, including junctional molecules, surface receptors and extracellular matrix, whereas granzyme F is cytotoxic in vitro48,49. Other aspects of the DETC UVR response, including rounding up and TIM-3 upregulation, did not obviously require Skint1 licensing. Possibly, licensing regulates those traits that are inherently more dangerous if precociously expressed.
Overall, our results provide experimental evidence supporting the lymphoid stress surveillance hypothesis, according to which lymphocytes display myelomonocytic traits, such as responding rapidly through innate receptors10. However, the results also establish the need for prior licensing through the TCR. Whereas there is no direct DETC equivalent in humans, this concept may pertain to the selection and regulation of human intestinal innate-like Vγ4+ IELs by BTNL3 and BTNL8, a bona fide TCR-Vγ ligand that is similarly expressed at steady state and is stress-labile15,20. This normality sensing mode of TCR stimulation may be linked to the unusually high TCR activation threshold acquired during Skint1 selection50, which could conceivably limit overt IEL activation in the absence of stress. In contrast to natural killer cells, which sense stress through missing self, DETC Skint1 sensing is a positive regulator of many components of tissue immunosurveillance, a role that may be phenocopied by BTNL sensing by murine and human gut IELs. The clinical implication of our findings is that epithelial pathologies may not always reflect failures to detect infectious, allergic or carcinogenic challenges but rather failings in the maintenance of normal homeostatic interactions before exposure to those challenges occur.
FVB/NJ strains (wild-type, FVB.Tac16, Trgv5–/–Trdv4–/– (ref. 12), Skint1-transgenic FVB.Tac24, huLangerin-DTA51, Tcrd–/– (ref. 52) and FVB.Tac.Tcrd–/–), C57BL/6J strains (wild-type, Nur77–GFP36, Tcrd–/– and Tnfsf9–/– (ref. 53)) and C57BL/6N strains (Skint1–/– and Skint locus–/–)54 were bred and maintained under specific pathogen-free conditions in the Biological Research Facility at the Francis Crick Institute with a 12-h light/12-h dark cycle and access to food and water ad libitum at a temperature of 19–23 °C and humidity of 55 ± 10%. Procedures were performed on mice of both sexes aged 6–15 weeks (with the exception of reconstitution experiments below) and were authorized by a UK Home Office License to A.C.H. and the Animal Welfare and Ethical Review Board of The Francis Crick Institute.
Neonatal DETC reconstitution
FVB Tcrd–/– and FVB.Tac Tcrd–/– pups at 1–3 d of age were transferred with a thymocyte suspension prepared from E16–E18 FVB wild-type fetuses by intraperitoneal injection. The equivalent of four to five fetal thymi were transferred into each pup. Reconstituted mice were analyzed at 8 weeks of age.
All antibodies were administered by a single i.d. injection into the ear pinna using a 30-gauge needle under inhalation anesthesia. For inhibition of Skint1 for confocal imaging, 2 μg of anti-Skint1 (clone 2G2) or rat IgG2a isotype control (BioLegend) were administered in 25 μl of PBS. For all other inhibition experiments, mice were administered 10 μg of anti-Skint1 (BioLegend), anti-4-1BBL (BioXCell), anti-GITRL (BioLegend) or rat IgG2a isotype control (BioLegend) in 10 μl of PBS. For antibody-mediated stimulation of receptors, mice were administered 2 μg of agonistic anti-4-1BB (R&D Systems) and anti-GITR (BioXCell) or respective isotype controls rat IgG2a (BioLegend) and rat IgG2b (BioXCell) in 10 μl of PBS. Anti-Skint1 hybridoma supernatant was produced by the Monoclonal Antibody Core Facility at the Helmholtz Centre Munich, and the antibody was purified by Cell Services at The Francis Crick Institute.
Under inhalation anesthesia, TEWL of the ear pinnae was measured using a Tewameter TM300 probe (Ck Electronic) under constant temperature and humidity conditions. Values were recorded as the average of ten consecutive measurements once readings were stabilized using MPA software (Ck Electronic).
Under inhalation anesthesia, mice were positioned under a cardboard box, and one or both ears were exposed and fixed flat using Blu-Tac. UVB radiation (302 nm) was administered to dorsal ear skin for a dose of 45 mJ cm–2 (C57BL/6) or 30 mJ cm–2 (FVB) using an XX-15 Series UV bench lamp (UVP) mounted on a Perspex and plastic housing.
Ex vivo ear skin culture
Ears were split and floated on complete RPMI (10% fetal calf serum, penicillin/streptomycin, 1× non-essential amino acids, 1 mM sodium pyruvate, 10 mM HEPES and 50 μM β-mercaptoethanol) containing 10 μg ml–1 anti-Skint1 or isotype control antibodies or 20 μM PP2 (Sigma) or DMSO vehicle alone. After 4 h of culture at 37 °C, ear skin was irradiated with 45 mJ cm–2 UVB as described above and returned to 37 °C for 24 h.
Mice were injected (i.d.) in the ear with 400 nmol ATP (Sigma) in 10 μl of PBS.
Irritant contact dermatitis model
Under inhalation anesthesia, 8 μl of 0.5% DNFB (Sigma) in 4:1 acetone:olive oil vehicle was applied to the dorsal ear skin and allowed to dry.
Epidermis was separated from dermis by floating split ear sheets on trypsin-GNK (2.5 mg ml–1 trypsin (Sigma-Aldrich), 7.35 mg ml–1 NaCl, 0.35 mg ml–1 KCl and 0.85 mg ml–1 glucose) for 1 h at 37 °C. Single-cell epidermal suspensions were prepared by subsequently digesting epidermal sheets in trypsin-GNK with 1× DNase buffer (1.21 mg ml–1 Tris base, 0.5 mg ml–1 MgCl2 and 73 μg ml–1 CaCl2) and 100 μg ml–1 DNase I (Roche) for 15 min at 37 °C with shaking at 700 r.p.m. Trypsin was neutralized with complete DMEM (Gibco), and the suspension was washed. Cells were stained with Live/Dead Fixable Near-IR dye (Molecular Probes) in PBS for 10 min at room temperature and washed in FACS buffer (PBS, 2% fetal calf serum, 2 mM EDTA and 0.04% sodium azide). Cells were blocked with mouse Fc Block (BD) and stained with combinations of the following conjugated antibodies, which are also listed in Supplementary Table 3, for 20 min at 4 °C: anti-CD3-BV510, anti-CD45-APC, anti-CD45-FITC, anti-CD45-PerCPCy5.5, anti-CD45.1-BV510, anti-CD69-PECy7, anti-CD86-APC, anti-CD86-FITC, anti-I-A/I-E-BV605, anti-TCRδ-BV421, anti-TCRδ-FITC, anti-GITR-APC (all BioLegend), anti-CD69-PE (eBioscience), anti-CD5-BV711 and anti-TCRβ-BV605 (all BD). Cells were then washed before immediate acquisition or alternatively fixed in CellFix (BD) for 10 min at 4 °C. Absolute cell counts were calculated using CountBright beads (Invitrogen). For apoptosis measurement, cells were washed in Annexin-V-binding buffer (BD) following surface staining, incubated with Annexin-V–APC (BioLegend) for 20 min at 4 °C, washed and immediately acquired. For staining of specific γδ TCRs, epidermal cells were rested overnight in complete RPMI at 37 °C before live/dead dye staining and were stained with 17D1 hybridoma supernatant55 at 4 °C for 20 min. Cells were then washed and stained with surface antibodies as above, including anti-Vγ5-APC (BioLegend) and anti-rat IgM-PE (eBioscience). For analysis of other trypsin-sensitive epitopes, the epidermis was separated and digested as above but with TrypLE Express (Gibco) in place of trypsin-GNK, and the epidermis was separated for 2 h. Cells were stained as above but including combinations of the following antibodies: anti-CD122-PE, anti-4-1BB-APC and anti-4-1BB-PE (all BioLegend). Cells were acquired with an LSR II or Fortessa (BD) and analyzed with FlowJo (TreeStar). Gating strategies are detailed in Supplementary Fig. 1.
For analysis of DETC dendrites and Vγ5 staining, ears were depilated with Nair and split, and the epidermis was separated by incubation in 0.5 M NH4SCN for 30 min at 37 °C and fixed in ice-cold acetone for 20 min at –20 °C.
For analysis of γH2A.X staining and for KC F-actin visualization in Extended Data Fig. 1, whole depilated split ear sheets were flattened in a histology cassette and fixed in neutral buffered formalin for 2 h at room temperature and permeabilized in 0.5% Triton X-100 for at least 1 h at room temperature.
All ear tissue was then washed in PBS with 1% bovine serum albumin (with 0.5% Triton X-100 for whole sheets), blocked with 5% normal goat serum for 1 h at room temperature and stained with combinations of the following antibodies and dyes (see also Supplementary Table 3) overnight at 4 °C: anti-CD3-Pacific Blue, anti-I-A/I-E-AF647, anti-Vγ5-APC (all BioLegend), anti-Vγ5-FITC (BD), anti-γH2A.X-FITC (Sigma-Aldrich), AF546-phalloidin or rhodamine-phalloidin (all Molecular Probes). Tissue was then washed and mounted in Vectashield Antifade mounting medium. Microscopy was performed using an LSM 710 (Zeiss; DETC dendrites) or TCS SP5 (Leica; γH2A.X) laser-scanning confocal microscope using a ×40 oil-immersion lens. All images were analyzed in Imaris (Bitplane).
For dendrite analysis, CD3+ cells (or Vγ5+ cells in the reconstitution experiments) were identified by the surfaces tool and used to mask F-actin staining. TCR- and F-actin-rich dendrites were identified by using the spots tool to detect 3-μm-diameter volumes containing an MFI in the top 25% of their respective channels. Colocalized CD3 and F-actin spots were selected for quantification on a per cell basis.
For γH2A.X analysis, the cell function was used to identify γ-H2A.X+ nuclei by filtering for the top 20% of cells by quality. Nuclei were quantified per image and averaged over three images per mouse.
Epidermal sheets were prepared, separated and fixed in acetone as above and blocked in 5% normal goat serum for 1 h at room temperature. Sheets were stained with anti-Skint1 for 1 h at 37 °C, washed and stained with anti-rat IgG AF568 (Invitrogen), anti-CD3-FITC, anti-Vγ5-FITC (all BD) and anti-TCRδ-AF647 (BioLegend). iSIM was performed using a VT-iSIM microscope (VisiTech International) with a ×100/1.45-NA objective. Raw images were processed with the Richardson–Lucy deconvolution algorithm and analyzed with NIS-Elements (Nikon) or Fiji (NIH) software. Skint1 staining relative to T cell bodies was quantified using Definiens Developer XD software (Definiens). Individual cells and Skint1 particles were detected using Otsu’s automatic multilevel thresholding and multiresolution segmentation algorithms. The MFI of Skint1 particles within and outside of T cell body masks was quantified.
Whole ears were flattened in a histology cassette and fixed in neutral buffered formalin for 24 h and dehydrated in 70% ethanol for 24 h. Ears were embedded, sectioned and stained with H&E by the Experimental Histopathology Facility at The Francis Crick Institute. Sections were photographed using a SlideScanner (Zeiss) at 10× magnification. Epidermal thickness was determined as the average of 30 manual measurements using QuPath software.
Ear swelling measurement
Total ear thickness was calculated as the average of three measurements of the ear pinna with a micrometer (Mitutoyo) after discarding the first measurement.
Measurement of CPDs
CPDs were measured by ELISA. Epidermal sheets were separated by NH4SCN as above. Epidermal DNA was isolated using a DNeasy Blood and Tissue kit (Qiagen) and denatured at 100 °C for 10 min. DNA (10 ng) was dried onto protamine sulfate-coated 96-well plates (Nunc), and ELISAs were performed using anti-CPDs (Cosmo Bio), biotinylated goat anti-mouse IgG (Invitrogen), streptavidin–horseradish peroxidase (R&D) and a TMB Substrate Reagent set (BD), as per the anti-CPDs manufacturer’s protocol. Absorbance was measured using an Infinite F200 Pro spectrophotometer (Tecan).
Toluidine blue permeability
Whole ears were dehydrated and rehydrated by sequential 2-min incubations in ice-cold 25%, 50%, 75%, 100%, 75%, 50% and 25% methanol in PBS and incubated in cold PBS for 5 min. Ears were stained in 0.1% (wt/vol) toluidine blue (Sigma-Aldrich) for 2 min and then destained in PBS for 1 min. Ears were photographed with an Stemi SV 11 stereomicroscope (Zeiss). Toluidine blue staining was quantified in Fiji by color deconvolution (H&E setting) and thresholding to measure total ear area (green channel) and toluidine blue area (indigo channel).
For analysis of total epidermal gene expression, epidermal sheets were separated by NH4SCN as above and homogenized in buffer RLT (Qiagen) using stainless steel beads (Qiagen) and a TissueLyser II (Qiagen). RNA was isolated using the RNeasy Mini kit (Qiagen) and an RNase-free DNase set (Qiagen).
For analysis of gene expression from DETCs, LCs and KCs, single-cell epidermal suspensions were prepared from ear sheets using trypsin-GNK and stained as above with Live/Dead Fixable Near-IR dye (Molecular Probes) and combinations of anti-CD3-BV510, anti-CD45-APC, anti-CD45-FITC, anti-CD45-Pacific Blue, anti-I-A/I-E-BV605, anti-I-A/I-E-FITC, anti-TCRδ-APC, anti-TCRδ-BV421 (all BioLegend) or anti-TCRδ-PerCPeFluor710 (eBioscience). CD45+TCRδ+ DETCs, CD45+MHCII+ LCs and CD45– KCs were sorted directly into buffer RLT Plus (Qiagen) using an Aria Fusion (BD). RNA was isolated using an RNeasy Micro Plus kit (Qiagen).
For analysis of gene expression from 4-1BBhi and 4-1BBlo DETCs, mice were shaved, and epidermal suspensions from ears and trunk skin were prepared using TrypLE Express (Gibco) as described above but with a 4-h incubation to separate the epidermis. Cells were stained with Live/Dead Fixable Near-IR dye (Molecular Probes) and anti-CD45-FITC, anti-TCRδ-BV421 and anti-4-1BB-PE (all BioLegend) and sorted. RNA was isolated as described above.
cDNA was synthesized using Superscript II Reverse Transcriptase, dNTP set and Anchored Oligo(dT)20 primer (Invitrogen). qPCRs were run with PowerUp SYBR Green (Applied Biosystems) and 0.4 μM primer sets listed in Supplementary Table 5 on a ViiA 7 Real-Time PCR System (Thermo Fisher). Threshold cycle (Ct) values were processed using QuantStudio software, and ΔCt values were calculated as 2(reference – gene), where the reference was CypA (cyclophilin).
Single-cell RNA sequencing
Wil-type FVB mice were treated with 10 μg of anti-Skint1 or rat IgG2a isotype control by i.d. injection into the ear pinna, as described above. Seven days later, epidermal suspensions were made from ear skin and were stained as described above using Live/Dead Fixable Near-IR dye (Molecular Probes), anti-EpCAM-PE (eBioscience) and anti-CD45-APC and anti-I-A/I-E-FITC (all BioLegend), and live CD45+MHCII–EpCAM– lymphocytes were sorted using an Aria Fusion (BD). Samples were prepared using the Chromium Single Cell 3′ GEM, Library and Gel Bead kit v3 (10x Genomics). Cell viability was assessed using an EVE cell counter and trypan blue viability stain. Approximately 10,000 cells were loaded into the 10x Chromium, which was operated according to manufacturer’s instructions. Sequencing was performed on an Illumina HiSeq 4000 with a read configuration recommended by 10x Genomics for the sequencing of these libraries (28-base pair (bp) read 1, 98-bp read 2, 8-bp index 1).
Bioinformatic analysis of single-cell RNA sequencing
Data were processed using CellRanger v3.0.2 using the prebuilt mm10 index v3.0.0, and individual Seurat objects were created. All of the following functions belong to the Seurat v3.1.1 package56, and their parameters are mentioned only if differing from the default. The following steps were applied to the data:
Cells with less than 500 features, more than 2,500 features and more than 12% mitochondrial reads were excluded.
Feature counts were normalized per cell using the NormalizeData function.
The FindVariableFeatures function was used to select the 2,000 most variable genes.
These genes were used to identify integration anchors (FindIntegrationAnchors, dims = 1:20, anchor.features = 2050), and the two samples (isotype control and anti-Skint1) were integrated using the IntegrateData function.
The expression of each gene was centered using the ScaleData function parametrized with all the available features.
Principal component analysis was performed with the RunPCA function (npcs = 100).
Clustering was performed using the FindNeighbors function (dims = 1:40) and the FindClusters function (resolution = 1.18).
Dimension reduction was performed using RunUMAP (dims = 1:40).
Clusters were manually annotated using the following marker genes: Ptprc, Cd3d, Cd3e, Trdc, Trgv5, Trdv4, Il17a, Scart2, Trac, Cd8a, Cd8b1, Cd4, Foxp3, Il10, Ncr1, Rora, Il7r, Csf3r, Cd207, Itgam, Krt5 and Krt10.
DETCs were identified by the expression of Ptprc, Cd3d, Cd3e, Trdc, Trgv5 and Trdv4, and their cell IDs were extracted.
The original Seurat objects were filtered using DETC IDs to remove other cell types from the analysis. Steps 1 to 8 were reapplied to the resulting data with the following differences: FindNeighbors (dims = 1:23), FindClusters (resolution = 0.28) and RunUMAP (dims = 1:23).
Marker genes for each cluster were identified using the FindAllMarkers function (only.pos = F, logfc.threshold = 0, min.pct = 0, min.cells.feature = 0 and min.cells.group = 0 and return.thresh = 2). For each cluster, marker genes were sorted by avg_logFC (increasing order), and GSEA v3.0 (ref. 57) was run for the C5 v6.2 ontology gene sets with the following parameters: -nperm 1000, -scoring_scheme classic and -set_min 3 and -set_max 50000. Likewise, for each cluster, marker genes were filtered (avg_logFC>0, p_val_adj<0.05) and used to perform a gene ontology analysis using the enrichGO function of the clusterProfiler v3.12.0 package58 (pAdjustMethod = BH, pvalueCutoff = 0.01 and qvalueCutoff = 0.05) with the org.Mm.eg.db v3.8.2 annotation.
Differential expression analysis between the two experimental conditions was performed using the FindMarkers function of Seurat (logfc.threshold = 0, min.pct = 0, min.cells.feature = 0 and min.cells.group = 0). The resulting gene list was sorted by avg_logFC (decreasing order), and GSEA v3.0 was run for the Skint1 selection signature with the following parameters: -nperm 1000, -scoring_scheme classic, -set_min 3 and -set_max 50000. The Skint1 selection signature was generated from the top 200 genes overexpressed in E15 Vγ5+ thymocytes from wild-type FVB versus Tac embryos (Supplementary Table 1), as determined by Illumina microarray analysis and described previously35.
Expression plots and violin plots were produced using the normalized counts.
Bulk RNA sequencing
Wild-type C57BL/6 mice were irradiated with 45 mJ cm–2 UVB as described above or were left unchallenged. Dorsal ear skin was collected at 24 h after UVR, and single-cell suspensions were made and stained as described above with Live/Dead Fixable Near-IR dye (Molecular Probes) and anti-CD3-PerCPCy5.5, anti-CD45-APC, anti-I-A/I-E-BV605 and anti-TCRδ-BV421 (all BioLegend). Live CD45+CD3+TCRδ+ DETCs were sorted directly into buffer RLT Plus (Qiagen) using an Aria Fusion (BD), and RNA was isolated using an RNeasy Micro Plus kit (Qiagen). cDNA was prepared from 10 ng of input RNA using the Ovation RNA-Seq System V2 (NuGEN). Final libraries were then constructed using the Ultralow Library System V2 (NuGEN). The resulting libraries were pooled for sequencing on an Illumina HiSeq 2500 platform with single-ended 75-bp reads.
Bioinformatic analysis of bulk RNA sequencing
Raw reads were quality and adapter trimmed using cutadapt v1.9.1 (ref. 59) before alignment. Reads were then aligned and quantified using RSEM v1.3.0 and STAR v2.5.2 (refs. 60,61) against the mouse genome GRCm38 and annotation release 86, both from Ensembl. Differential gene expression analysis was performed in R using the DESeq2 v1.24.0 package62. Differential genes were selected using a 0.05 false discovery rate threshold. Normalization and VST was applied on raw counts before performing principal component analysis and Euclidean distance-based clustering. For the heat maps, Z scores were computed on the VST of the raw counts. Heat maps were made in R using the ComplexHeatmap v2.0.0 package63. GSEA was run as described above using gene signatures created from the top 200 overexpressed genes in each DETC cluster from the single-cell RNA-sequencing experiment.
Data were analyzed in Prism 9.3.0 (GraphPad) and R v3.6.0. Statistical tests were performed as indicated in the figure legends, and n values are also provided. All error bars represent mean ± s.d. Unless explicitly stated in the figure legends, all dots and n values refer to individual mice. Mice were randomly assigned to treatment groups where applicable. No data were excluded. Data were presumed to be normally distributed, and all tests were two tailed where applicable. Statistical significance was defined as P < 0.05, and all significant P values are shown in the figures. Welch’s correction was applied to t-tests where variance was significantly different between sample groups as determined by an F-test.
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
RNA-sequencing datasets generated in this paper are available under accession codes GSE160476 (single cell) and GSE160477 (bulk). Previously published datasets used in this paper are provided in the Supplementary Tables. The public gene set used in Extended Data Fig. 3e is available at GO:0050778. Source data are provided with this paper.
Code used for single-cell RNA-sequencing processing is available at https://github.com/bahnk/scRNASeq-DETCs.
Kipnis, J., Gadani, S. & Derecki, N. C. Pro-cognitive properties of T cells. Nat. Rev. Immunol. 12, 663–669 (2012).
Schreiber, R. D., Old, L. J. & Smyth, M. J. Cancer immunoediting: integrating immunity’s roles in cancer suppression and promotion. Science 331, 1565–1570 (2011).
Zmora, N., Bashiardes, S., Levy, M. & Elinav, E. The role of the immune system in metabolic health and disease. Cell Metab. 25, 506–521 (2017).
Girardi, M., Lewis, J. M., Filler, R. B., Hayday, A. C. & Tigelaar, R. E. Environmentally responsive and reversible regulation of epidermal barrier function by γδ T cells. J. Invest. Dermatol. 126, 808–814 (2006).
Varol, C., Mildner, A. & Jung, S. Macrophages: development and tissue specialization. Annu. Rev. Immunol. 33, 643–675 (2015).
Vivier, E. et al. Innate lymphoid cells: 10 years on. Cell 174, 1054–1066 (2018).
Panduro, M., Benoist, C. & Mathis, D. Tissue Tregs. Annu. Rev. Immunol. 34, 609–633 (2016).
Hirano, M. et al. Evolutionary implications of a third lymphocyte lineage in lampreys. Nature 501, 435–438 (2013).
Landau, S. B., Aziz, W. I., Woodcock-Mitchell, J. & Melamede, R. Vγ (I) expression in human intestinal lymphocytes is restricted. Immunol. Invest. 24, 947–955 (1995).
Hayday, A. C. γδ T cells and the lymphoid stress-surveillance response. Immunity 31, 184–196 (2009).
Girardi, M. et al. Regulation of cutaneous malignancy by γδ T cells. Science 294, 605–609 (2001).
Strid, J. et al. Acute upregulation of an NKG2D ligand promotes rapid reorganization of a local immune compartment with pleiotropic effects on carcinogenesis. Nat. Immunol. 9, 146–154 (2008).
Asarnow, D. M. et al. Limited diversity of γδ antigen receptor genes of Thy-1+ dendritic epidermal cells. Cell 55, 837–847 (1988).
Kyes, S., Carew, E., Carding, S. R., Janeway, C. A. & Hayday, A. Diversity in T-cell receptor γ gene usage in intestinal epithelium. Proc. Natl Acad. Sci. USA 86, 5527–5531 (1989).
Di Marco Barros, R. et al. Epithelia use butyrophilin-like molecules to shape organ-specific γδ T cell compartments. Cell 167, 203–218 (2016).
Lewis, J. M. et al. Selection of the cutaneous intraepithelial γδ+ T cell repertoire by a thymic stromal determinant. Nat. Immunol. 7, 843–850 (2006).
Boyden, L. M. et al. Skint1, the prototype of a newly identified immunoglobulin superfamily gene cluster, positively selects epidermal γδ T cells. Nat. Genet. 40, 656–662 (2008).
Jandke, A. et al. Butyrophilin-like proteins display combinatorial diversity in selecting and maintaining signature intraepithelial γδ T cell compartments. Nat. Commun. 11, 3769 (2020).
Salim, M. et al. Characterization of a putative receptor binding surface on Skint-1, a critical determinant of dendritic epidermal T cell selection. J. Biol. Chem. 291, 9310–9321 (2016).
Willcox, C. R. et al. Butyrophilin-like 3 directly binds a human Vγ4+ T cell receptor using a modality distinct from clonally-restricted antigen. Immunity 51, 813–825 (2019).
Melandri, D. et al. The γδTCR combines innate immunity with adaptive immunity by utilizing spatially distinct regions for agonist selection and antigen responsiveness. Nat. Immunol. 19, 1352–1365 (2018).
Chodaczek, G., Papanna, V., Zal, M. A. & Zal, T. Body-barrier surveillance by epidermal γδ TCRs. Nat. Immunol. 13, 272–282 (2012).
Hanahan, D. & Weinberg, R. A. Hallmarks of cancer: the next generation. Cell 144, 646–674 (2011).
Barbee, S. D. et al. Skint-1 is a highly specific, unique selecting component for epidermal T cells. Proc. Natl Acad. Sci. USA 108, 3330–3335 (2011).
Jameson, J. M., Cauvi, G., Witherden, D. A. & Havran, W. L. A keratinocyte-responsive γδ TCR is necessary for dendritic epidermal T cell activation by damaged keratinocytes and maintenance in the epidermis. J. Immunol. 172, 3573–3579 (2004).
Girardi, M. et al. Resident skin-specific γδ T cells provide local, nonredundant regulation of cutaneous inflammation. J. Exp. Med. 195, 855–867 (2002).
Dalessandri, T., Crawford, G., Hayes, M., Castro Seoane, R. & Strid, J. IL-13 from intraepithelial lymphocytes regulates tissue homeostasis and protects against carcinogenesis in the skin. Nat. Commun. 7, 12080 (2016).
Lee, H.-J. et al. Epidermal barrier function is impaired in Langerhans cell-depleted mice. J. Invest. Dermatol. 139, 1182–1185 (2019).
Abeler-Dörner, L. et al. High-throughput phenotyping reveals expansive genetic and structural underpinnings of immune variation. Nat. Immunol. 21, 86–100 (2020).
Zaiss, D. M. W. et al. Amphiregulin enhances regulatory T cell-suppressive function via the epidermal growth factor receptor. Immunity 38, 275–284 (2013).
Shoyab, M., Plowman, G. D., McDonald, V. L., Bradley, J. G. & Todaro, G. J. Structure and function of human amphiregulin: a member of the epidermal growth factor family. Science 243, 1074–1076 (1989).
Defrance, T. et al. Interleukin 13 is a B cell stimulating factor. J. Exp. Med. 179, 135–143 (1994).
Kadow, S. et al. Aryl hydrocarbon receptor is critical for homeostasis of invariant γδ T cells in the murine epidermis. J. Immunol. 187, 3104–3110 (2011).
DeBenedette, M. A. et al. Role of 4-1BB ligand in costimulation of T lymphocyte growth and its upregulation on M12 B lymphomas by cAMP. J. Exp. Med. 181, 985–992 (1995).
Turchinovich, G. & Hayday, A. C. Skint-1 identifies a common molecular mechanism for the development of interferon-γ-secreting versus interleukin-17-secreting γδ T cells. Immunity 35, 59–68 (2011).
Moran, A. E. et al. T cell receptor signal strength in Treg and iNKT cell development demonstrated by a novel fluorescent reporter mouse. J. Exp. Med. 208, 1279–1289 (2011).
MacLeod, A. S. et al. Skin-resident T cells sense ultraviolet radiation-induced injury and contribute to DNA repair. J. Immunol. 192, 5695–5702 (2014).
Strid, J., Sobolev, O., Zafirova, B., Polic, B. & Hayday, A. The intraepithelial T cell response to NKG2D-ligands links lymphoid stress surveillance to atopy. Science 334, 1293–1297 (2011).
Jameson, J. et al. A role for skin γδ T cells in wound repair. Science 296, 747–749 (2002).
Shipman, W. D. et al. A protective Langerhans cell–keratinocyte axis that is dysfunctional in photosensitivity. Sci. Transl. Med. 10, eaap9527 (2018).
Komori, H. K. et al. Cutting edge: dendritic epidermal γδ T cell ligands are rapidly and locally expressed by keratinocytes following cutaneous wounding. J. Immunol. 188, 2972–2976 (2012).
Wollenberg, A. et al. Treatment of atopic dermatitis with tralokinumab, an anti-IL-13 mAb. J. Allergy Clin. Immunol. 143, 135–141 (2019).
Ghilas, S. et al. Natural killer cells and dendritic epidermal γδ T cells orchestrate type 1 conventional DC spatiotemporal repositioning toward CD8+ T cells. iScience 24, 103059 (2021).
Silva-Santos, B., Serre, K. & Norell, H. γδ T cells in cancer. Nat. Rev. Immunol. 15, 683–691 (2015).
Zhang, B. et al. Differential requirements of TCR signaling in homeostatic maintenance and function of dendritic epidermal T cells. J. Immunol. 195, 4282–4291 (2015).
Keyes, B. E. et al. Impaired epidermal to dendritic T cell signaling slows wound repair in aged skin. Cell 167, 1323–1338 (2016).
Mohamadzadeh, M. et al. Functional roles for granzymes in murine epidermal γδ T-cell-mediated killing of tumor targets. J. Invest. Dermatol. 107, 738–742 (1996).
Hiebert, P. R. & Granville, D. J. Granzyme B in injury, inflammation, and repair. Trends Mol. Med. 18, 732–741 (2012).
Shi, L., Wu, L., Wang, S. & Fan, Z. Granzyme F induces a novel death pathway characterized by Bid-independent cytochrome c release without caspase activation. Cell Death Differ. 16, 1694–1706 (2009).
Wencker, M. et al. Innate-like T cells straddle innate and adaptive immunity by altering antigen-receptor responsiveness. Nat. Immunol. 15, 80–87 (2014).
Kaplan, D. H., Jenison, M. C., Saeland, S., Shlomchik, W. D. & Shlomchik, M. J. Epidermal Langerhans cell-deficient mice develop enhanced contact hypersensitivity. Immunity 23, 611–620 (2005).
Itohara, S. et al. T cell receptor δ gene mutant mice: independent generation of αβ T cells and programmed rearrangements of γδ TCR genes. Cell 72, 337–348 (1993).
Skarnes, W. C. et al. A conditional knockout resource for the genome-wide study of mouse gene function. Nature 474, 337–342 (2011).
Narita, T., Nitta, T., Nitta, S., Okamura, T. & Takayanagi, H. Mice lacking all of the Skint family genes. Int. Immunol. 30, 301–309 (2018).
Mallick-Wood, C. A. et al. Conservation of T cell receptor conformation in epidermal γδ cells with disrupted primary Vγ gene usage. Science 279, 1729–1733 (1998).
Butler, A., Hoffman, P., Smibert, P., Papalexi, E. & Satija, R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat. Biotechnol. 36, 411–420 (2018).
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).
Yu, G., Wang, L.-G., Han, Y. & He, Q.-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284–287 (2012).
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet j. 17, 10 (2011).
Li, B. & Dewey, C. N. RSEM: accurate transcript quantification from RNA-seq data with or without a reference genome. BMC Bioinformatics 12, 323 (2011).
Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21 (2013).
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Gu, Z., Eils, R. & Schlesner, M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 32, 2847–2849 (2016).
We thank J. Strid (Imperial College London) and T. Zal (UT Southwestern) for advice; M. Battilocci and F. Watt (KCL) for advice and equipment and D. Calado (Francis Crick Institute) and T. Nitta, T. Narita and H. Takayanagi (University of Tokyo) for mice. We thank present and past members of our laboratory and the laboratory of D. Gibbons for discussions and advice, in particular O. Sobolev, A. Jandke, L. Monin, D. Moreno-Vicencio, V. Sofra, F. Cano and N. Roberts. We thank staff from the following Francis Crick Institute science technology platforms: biological research facility, flow cytometry, experimental histopathology, light microscopy, advanced sequencing, genomics equipment park (particularly O. O’Neill) and cell services and staff at the Nikon Imaging Centre at KCL Guy’s Hospital campus. This work was supported by The Francis Crick Institute, which receives core funding from Cancer Research UK (FC001093), the MRC (FC001093) and the Wellcome Trust (FC001093) (A.C.H.); by the Cancer Research UK, King’s and City of London Cancer Centre; by the Wellcome Trust (106292/Z/14/Z and 100156/Z/12/Z; A.C.H.); by a Cancer Research Institute Irvington fellowship to D.R.M. and by a Cancer Research UK postdoctoral fellowship to R.H.
A.C.H. is an equity holder in and consultant to GammaDelta Therapeutics, Adaptate Biotherapeutics and ImmunoQure AG. The remaining authors declare no competing interests.
Peer review information
Nature Immunology thanks Kenji Kabashima and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Ioana Visan was the primary editor on this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended Data Fig. 1 Steady state DETC-keratinocyte interactions depend on the Vγ5Vδ1 TCR and Skint1.
a) Confocal microscopy of CD3+ DETCs (cyan) and total F-actin or that within CD3 mask (yellow) in epidermal sheets from wild-type and FVB.Tac mice. Representative of n=4/group. Scale bar 20μm. b) Neonatal Tcrd-/- (n=13) and FVB.Tac.Tcrd-/- (n=11) pups were reconstituted with wild-type embryonic thymocytes containing mature Skint1-selected DETC precursors. Confocal microscopy of CD3+ DETCs (cyan) and F-actin within CD3 mask (yellow) in epidermal sheets at 8 weeks of age. Scale bar 20μm. Frequency distribution of number of interacting dendrites within Vγ5+ DETCs (CD3- and F-actin-rich aggregates) per cell using mean values per mouse. c) Confocal microscopy of CD3+ DETCs (cyan) and total F-actin or that within CD3 mask (yellow) in epidermal sheets from wild-type C57BL/6 mice pre-treated with Skint1 Ab or isotype control (n=5/group). Scale bar 20μm. Frequency distribution of number of interacting dendrites (CD3- and F-actin-rich aggregates) per cell using mean values per mouse. d) Instant structured illumination microscopy (iSIM) of Skint1 (magenta) or isotype control in whole mount epidermal sheets from wild-type FVB mice. Scale bar 10μm. e) iSIM of CD3 (cyan), TCRδ (white) and Skint1 (magenta) staining in whole mount epidermal sheets from wild-type FVB(n=3), FVB.Tac (n=5) and Tcrd-/- (n=2) mice. Scale bar 10μm. Mean Skint1 MFI quantified within TCRδ (wild-type, FVB.Tac) or CD3 (Tcrd-/-) masks relative to mean image background value. Each dot represents a single cell. f) qPCR for Skint gene expression from epidermis of wild-type FVB (n=11) and Tcrd-/- mice (n=6). Representative images of and cumulative data from 4 (b) or 3 (e) independent experiments. Pooled data from 2 independent experiments (f). Matched two-way (b-c) or ordinary one-way (e) ANOVA with Sidak’s multiple comparisons test, unpaired two-tailed t-test (f). Values in pie charts and error bars on graph represent mean ± SD. ns: not significant.
a) Flow cytometry of epidermal CD45+ cellular composition and b) canonical TCR expression by DETCs from wild-type FVB mice administered with Skint1 Ab or isotype control by intradermal injection into the ear weekly for 10 weeks (n=8/group). c) Confocal microscopy of Vγ5+ DETCs (green) and MHC-II+ LCs (white) from mice treated as in a. Representative of n=8/group, scale bar 50 μm. d) Total ear thickness of wild-type FVB (n=5/group) or Tcrd-/- mice (n=4/group) administered with Skint1 Ab or isotype control by intradermal injection into the ear weekly for 5 weeks. e) Epidermal thickness measured from H&E sections of ear skin of mice treated as in d. f) Transepidermal water loss (TEWL) of injected and contralateral ear skin of wild-type mice treated as in d. g) Baseline TEWL in C57BL/6 Skint1+/- (n=7), Skint1-/- (n=7) and Skint locus-/- mice (n=5). h) qPCR of Il13 expression in epidermis of untreated wild-type FVB (n=14) and Tcrd-/- mice (n=9). i) qPCR of keratinocyte gene expression in epidermis of untreated wild-type FVB (n=15) and Tcrd-/- mice (n=9). j) qPCR of keratinocyte gene expression in epidermis of mice treated as in a (n=8/group). k) Baseline TEWL in wild-type FVB (n=18), Langerin-DTA (n=11) or Tcrd-/- mice (n=6). l) qPCR of keratinocyte gene expression in epidermis of wild-type FVB (n=7) and Langerin-DTA mice (n=5). Data representative of at least 2 independent experiments (a-b, d, f, k-l) or pooled from 3 (h-i) 2 (j) independent experiments. Ordinary two-way ANOVA with Tukey’s multiple comparisons test (d-e), ordinary one-way ANOVA with Dunnett’s multiple comparisons test relative to Skint1+/- (g), paired two-way ANOVA with Sidak’s multiple comparisons test (f), unpaired two-tailed t-test (h-j, l) or ordinary one-way ANOVA with Tukey’s multiple comparisons test (k). Error bars represent mean ± SD.
a) Expression of Areg in DETCs from the combined dataset of single cell RNA sequencing of DETCs from wild-type FVB mice treated 7d prior with Skint1 Ab or isotype control by intradermal injection (n=1/group). b) Expression of selected Cluster 0 marker genes from the combined dataset as in a. c) Selected enriched gene ontology terms in each DETC cluster from the combined dataset as in a. d) Seurat cell cycle analysis of DETCs from the combined dataset as in a. e) Gene set enrichment analysis (GSEA) comparing Cluster 3 with other clusters in the combined dataset as in a, using publicly available gene set “GO positive regulation of immune response”. f) qPCR for cluster 3 marker gene expression within sorted 4-1BBhi and 4-1BBlo DETCs from untreated wild-type FVB mice (n=9). g) qPCR validation of differentially expressed genes in total DETC sorted from mice treated as in a (isotype n=9, Skint1 Ab n=8). h) Differentially expressed genes in Clusters 0 and 1 comparing Skint1 Ab and isotype treatments as in a. i) Expression of Ptpn6 in DETCs comparing Skint1 Ab and isotype treatments as in a. j) Flow cytometry of CD122 and CD45RB expression by DETCs from mice treated as in a (n=5/group). k) Flow cytometry of CD45RB expression by TCRγδ+ DETCs from wild-type FVB and FVB.Tac mice treated as in a (n=4/group). l) Flow cytometry of 4-1BBhi DETCs and CD122 expression by DETCs from wild-type FVB (n=7) and Langerin-DTA mice (n=5). Pooled data from two experiments (f-g) or representative of at least two experiments (j, l). Hypergeometric test with Benjamini-Hochberg correction (c), paired two-tailed t-test (f), unpaired two-tailed t-test (g, j, l), Wilcoxon test with Bonferroni’s correction (h), ordinary two-way ANOVA with Sidak’s multiple comparisons test (k). Error bars represent mean ± SD.
a) qPCR expression of Skint1 and Skint2 in KCs sorted from wild-type C57BL/6 mice at indicated timepoints following UVR (n=3/timepoint). b) Flow cytometry frequency of activated LCs following UVR of ear skin in wild-type C57BL/6 mice (n=3/timepoint). c) Flow cytometry of CD69 expression by DETCs following UVR of ear skin in wild-type C57BL/6 mice (n=3/timepoint). d) Flow cytometry analysis of fold change in CD69 gMFI on DETCs 24hr post-UVR relative to unchallenged skin from wild-type FVB (n=25), FVB.Tac (n=13), Trgv5-/-Trdv4-/- (n=13) and Skint1-transgenic FVB.Tac mice (n=6). e) qPCR for effector genes in DETCs sorted from wild-type C57BL/6 mice at 24hr post-UVR, following pre-treatment with intradermal injection of Skint1 Ab or isotype control (isotype normal n=8, Skint1 Ab normal n=8, isotype UVR n=10, Skint1 Ab UVR n=10). f) Flow cytometry frequency of activated LCs from mice treated as in e) (n=7/group). Flow cytometry of GFP expression by DETCs from ear sheets of Nur77-GFP mice which were cultured for 4hr in medium containing g) anti-Skint1 or isotype control antibodies and h) PP2 or DMSO vehicle prior to UVR, and analyzed 24hr later (n=3/group). Pooled data from 4 (d) or 2 (e) independent experiments, or representative of 3 experiments (g-h). Ordinary one-way ANOVA with Dunnett’s multiple comparisons test relative to 0hr (a) or wild-type (d), paired two-way ANOVA with Sidak’s multiple comparisons test (e, f-h). Error bars represent mean ± SD.
a) qPCR for gene expression (DETCs/LCs n=4, KCs n=3) and b) flow cytometry for surface expression of 4-1BB and GITR in epidermal cell populations from unchallenged wild-type C57BL/6 mice. Representative of n=3 mice (b). c) qPCR for DETC gene expression in epidermis at 24hr post-UVR from wild-type C57BL/6 (n=14) and Tnfsf9-/- mice (n=15). d) Tnfrsf18 expression in steady state DETC clusters identified from scRNAseq as in Fig. 3. Pooled data from 3 independent experiments (c). Paired two-way ANOVA with Sidak’s multiple comparisons test (c). Error bars represent mean ± SD.
a) Wild-type C57BL/6 mice were administered Skint1 Ab or isotype control by intradermal injection 24hr prior to and 48hr post-ultraviolet radiation (UVR) of dorsal ear skin. Representative confocal microscopy images of γH2A.X nuclei (magenta) and CD3+ DETCs (cyan) in ear sheets, as quantified in Fig. 6b at d3 post-UVR. Righthand images overlay F-actin staining (yellow). Scale bar 50μm (n=10/group). Relative increase in ear thickness at 5d post-UVR relative to unchallenged ear in b) wild-type C57BL/6 (n=6) vs Tcrd-/- mice (n=5), c) littermate control (n=10) vs Skint1-/- mice (n=10) and d) littermate control (n=4) vs Skint locus-/- mice (n=5). Representative of (a) or pooled data from (c) 2 independent experiments. Unpaired two-tailed t-test (b-d) with Welch’s correction where necessary (b-c). Error bars represent mean ± SD.
Statistical source data.
Statistical source data.
Statistical source data.
Statistical source data.
Statistical source data.
Statistical source data.
Statistical source data.
Statistical source data.
Statistical source data.
Statistical source data.
Statistical source data.
Statistical source data.
Statistical source data.
About this article
Cite this article
McKenzie, D.R., Hart, R., Bah, N. et al. Normality sensing licenses local T cells for innate-like tissue surveillance. Nat Immunol 23, 411–422 (2022). https://doi.org/10.1038/s41590-021-01124-8
Nature Communications (2022)
Nature Reviews Immunology (2022)
Nature Immunology (2022)