The interferon-rich skin environment regulates Langerhans cell ADAM17 to promote photosensitivity in lupus
Abstract
The autoimmune disease lupus erythematosus (lupus) is characterized by photosensitivity, where even ambient ultraviolet radiation (UVR) exposure can lead to development of inflammatory skin lesions. We have previously shown that Langerhans cells (LCs) limit keratinocyte apoptosis and photosensitivity via a disintegrin and metalloprotease 17 (ADAM17)-mediated release of epidermal growth factor receptor (EGFR) ligands and that LC ADAM17 sheddase activity is reduced in lupus. Here, we sought to understand how the lupus skin environment contributes to LC ADAM17 dysfunction and, in the process, differentiate between effects on LC ADAM17 sheddase function, LC ADAM17 expression, and LC numbers. We show through transcriptomic analysis a shared IFN-rich environment in non-lesional skin across human lupus and three murine models: MRL/lpr, B6.Sle1yaa, and imiquimod (IMQ) mice. IFN-I inhibits LC ADAM17 sheddase activity in murine and human LCs, and IFNAR blockade in lupus model mice restores LC ADAM17 sheddase activity, all without consistent effects on LC ADAM17 protein expression or LC numbers. Anti-IFNAR-mediated LC ADAM17 sheddase function restoration is associated with reduced photosensitive responses that are dependent on EGFR signaling and LC ADAM17. Reactive oxygen species (ROS) is a known mediator of ADAM17 activity; we show that UVR-induced LC ROS production is reduced in lupus model mice, restored by anti-IFNAR, and is cytoplasmic in origin. Our findings suggest that IFN-I promotes photosensitivity at least in part by inhibiting UVR-induced LC ADAM17 sheddase function and raise the possibility that anifrolumab ameliorates lupus skin disease in part by restoring this function. This work provides insight into IFN-I-mediated disease mechanisms, LC regulation, and a potential mechanism of action for anifrolumab in lupus.
Editor's evaluation
This study presents a useful assessment of the possible role of type I interferons in inhibiting Adam17 protease/sheddase activity and their correlation with decreased Langerhans Cells signature in lesional and nonlesional CLE and murine models as cause of photosensitive lupus. The data were collected and analyzed using solid methodology. This work will be of interest to scientists interested in photosensitivity in the setting of lupus.
https://doi.org/10.7554/eLife.85914.sa0Introduction
The autoimmune disease lupus erythematosus (lupus) manifests in part by photosensitivity, a sensitivity to ultraviolet radiation (UVR) whereby even ambient sunlight can trigger the development of inflammatory skin lesions (Foering et al., 2013; Kim and Chong, 2013). In patients with the systemic form of the disease (SLE), with or without diagnosed cutaneous lupus (CLE), photosensitive skin responses can also be associated with worsening systemic autoimmunity and end organ damage. Photosensitive skin inflammation, the accompanying risk of systemic disease flares, and the lifestyle modifications such as sun avoidance needed to prevent UVR-induced skin inflammation all contribute to reduced quality of life (Bachen et al., 2009; Foering et al., 2012). Insight into mechanisms that underlie photosensitivity remain limited (Estadt et al., 2022; Sim et al., 2021), and better understanding of pathophysiology will help to advance therapeutic options.
To study the pathophysiology that drives photosensitivity, we have focused on studying non-lesional skin to delineate mechanisms that have already gone awry that lead to photosensitivity. We recently reported that Langerhans cells (LCs) can limit UVR-induced skin inflammation by expressing a disintegrin and metalloprotease 17 (ADAM17), which releases membrane-bound epidermal growth factor receptor (EGFR) ligands that then preserve epidermal integrity via EGFR stimulation of keratinocytes (Shipman et al., 2018). Two distinct photosensitive SLE mouse models showed a reduction in LC Adam17 mRNA in non-lesional skin which corresponded to reduced LC ADAM17 sheddase activity, and we could reduce photosensitivity by circumventing the reduced EGFR ligand shedding with topical EGFR ligand application. Together, these data suggested that reduction of LC ADAM17-dependent EGFR-ligand sheddase activity contributes to photosensitivity. However, no murine ADAM17 antibody was available to differentiate between diminished LC ADAM17 protein level versus sheddase function. Non-lesional skin of SLE patients showed reduced epidermal EGFR phosphorylation and a subset had reduced LC numbers, suggesting that reduced LC ADAM17 sheddase function, LC ADAM17 expression and/or LC numbers could potentially contribute to pathology. Thus, our previous study put forth a new model for LC dysfunction in lupus photosensitivity but the mechanisms that lead to LC dysfunction and distinguishing among the regulation of LC ADAM17 sheddase function, protein expression, and LC numbers remained to be fully elucidated.
A prominent interferon (IFN) signature indicative of exposure to type I interferon (IFN-I) is found in tissues and circulating cells in lupus (Baechler et al., 2003; Bennett et al., 2003; Catalina et al., 2019; Der et al., 2019; Jabbari et al., 2014; Kirou et al., 2004; Martínez et al., 2022; Psarras et al., 2020; Psarras et al., 2022), and the recent FDA approval of anifrolumab (anti-IFNAR1) for SLE highlights the importance of IFN-I in disease pathogenesis. Numerous studies have pointed to an IFN-I-rich environment in even non-lesional skin based on transcriptomic signatures of skin from SLE and CLE patients (Billi et al., 2022; Der et al., 2017; Der et al., 2019; Psarras et al., 2020), upregulation of keratinocyte IFN-κ expression in keratinocytes of CLE, SLE, and ANA +patients (Psarras et al., 2020; Stannard et al., 2017), and upregulation of IFN-stimulated genes (ISGs) such as MX1 on tissue sections in SLE, incomplete SLE, and ANA +patients (Lambers et al., 2019; Psarras et al., 2020; Reefman et al., 2008). Remarkably, the IFN signature in both lesional and non-lesional skin is greatly enriched compared to the signature in blood (Psarras et al., 2020), further pointing to a potential pathogenic role for IFN-I in lupus skin disease, and the IFN-I-rich environment in grossly unaffected skin points to a role for IFN-I in potentially predisposing to lesion development. Functionally, the reduced skin lesions in IFNAR-deficient SLE model mice (Nickerson et al., 2013) is consistent with this idea, and, indeed, anifrolmab is especially efficacious for skin disease in SLE patients (Furie et al., 2017; Merrill et al., 2018; Morand et al., 2020). Mechanisms by which IFN-I contributes to and anifrolumab exerts this ameliorative effect on lupus skin disease, however, remains to be better understood.
Here, we consider LC dysfunction in the context of the high IFN levels in lupus, and test the hypothesis that the IFN-I-rich environment present in even non-lesional lupus skin promotes LC dysfunction to contribute to photosensitivity. We extend accumulating human data showing an IFN-I signature in non-lesional CLE skin and establish that non-lesional skin from photosensitive MRL/lpr, B6.Sle1yaa, and imiquimod (IMQ) SLE models share IFN signatures with human skin. IFN-I was sufficient to inhibit LC ADAM17 sheddase activity without consistent effects on surface expression levels or altering LC numbers in healthy murine and human LCs. Conversely, anti-IFNAR treatment of SLE mouse models restored LC ADAM17 sheddase function without affecting expression levels or LC numbers. Anti-IFNAR restoration of LC ADAM17 sheddase function was associated with reduced photosensitive response, and this anti-IFNAR effect was dependent on EGFR signaling and LC ADAM17. Lastly, we show that UVR-induced reactive oxygen species (ROS) known to promote ADAM17 sheddase activity (Singh et al., 2009), was reduced in lupus model LCs, that anti-IFNAR restored ROS expression, and that ROS was of cytoplasmic origin. Together, our results establish that multiple photosensitive murine SLE models are similar to human lupus patients in demonstrating IFN signatures in non-lesional skin and put forth a mechanism whereby IFN-I contributes to photosensitivity at least in part by inhibiting LC sheddase ADAM17 activity. These data provide insight into IFN-I-mediated contributions to photosensitivity, delineate a driver of LC dysfunction, and suggest the possibility that restoration of LC ADAM17 function is a mechanism of action for anifrolumab in lupus skin disease.
Results
Non-lesional skin from human lupus and murine models share an IFN-rich environment
To understand mechanisms that lead to the LC defects seen in non-lesional skin of SLE mouse models and human SLE (Shipman et al., 2018), we assessed the non-lesional skin transcriptomic profiles from human lupus patients and multiple murine models. We examined a non-lesional skin microarray dataset from a CLE cohort with the discoid form of CLE (DLE) that had not been previously analyzed and was part of a cohort whose lesional skin was shown to have a Th1 signature when compared to healthy controls and psoriatic skin (Jabbari et al., 2014). The gene expression profile of the non-lesional DLE skin was significantly different from that of healthy skin (Supplementary file 1), lesional skin, and psoriatic skin (Figure 1A). Differential pathway analysis via Quantitative Set Analysis for Gene Expression (QuSAGE) showed that IFN, IFN-I and IFN-γ pathways were among the upregulated pathways in non-lesional DLE skin compared to normal controls (Figure 1B and Supplementary file 2). The IFN-I pathway was more prominent than the IFNg pathway, with upregulation of a number of IFN-I-associated transcription factors such as IRF1, 6, 7, 8 and IFN-stimulated genes including MX1, MX2, XAF1, IFI27, and ISG15 in non-lesional skin (Figure 1C–D, Figure 1—figure supplement 1). Comparison of non-lesional skin with lesional skin showed that this IFN response was overall less dramatic in non-lesional skin (Figure 1D, Figure 1—figure supplement 1), consistent with findings of previous studies (Billi et al., 2022; Martínez et al., 2022; Psarras et al., 2020). The more muted response did not apply to every gene in the pathway, however, with ADAR, IRF6, and PTPN1 expressed at higher levels in non-lesional skin (Figure 1D, Figure 1—figure supplement 1) and suggesting a distinct biology in non-lesional and lesional skin. These data are consistent with that from other cohorts (Billi et al., 2022; Der et al., 2017; Der et al., 2019; Psarras et al., 2020) showing that non-lesional skin, similar to lesional skin, has an IFN-rich environment.
We also analyzed the DLE data using Gene Set Variation Analysis (GSVA) to examine expression of defined lupus-relevant gene sets. We had previously used this approach to analyze publicly available bulk transcriptomic data from non-lesional and lesional CLE skin that included the DLE lesional skin data (Martínez et al., 2022). Reanalysis of the lesional skin here reiterated the GSVA results of multiple datasets (Martínez et al., 2022), showing enrichment of immune cell gene sets including plasma cells, T cells, monocyte/myeloid cells, neutrophils, and plasmacytoid dendritic cells (pDCs) when compared to controls (Figure 1E; Supplementary file 3). The LC gene set was de-enriched in lesional skin, echoing the LC loss seen in other studies (Billi et al., 2022; Bos et al., 1986; Sontheimer and Bergstresser, 1982), and this effect was more dramatic than the LC de-enrichment in lesional psoriatic skin (Figure 1E). Lesional DLE samples were also enriched for inflammatory cytokine gene sets, including IL12, IL21, IL23, TNF, and a highly upregulated IFN gene signature as previously seen (Figure 1E; Martínez et al., 2022). The upregulation of IFN response genes in lesional DLE skin extended to gene signatures indicative of response to specific IFNs, including IFN-Is (IFN-α, β, κ, and ω) and type II IFNs (IFN-γ; Figure 1F). The IFN-γ signature is consistent with the Th1 signature identified in the original study examining these gene expression data along with cellular phenotyping (Jabbari et al., 2014).
Non-lesional skin also showed upregulation of the IFN gene set in GSVA when compared to healthy controls, although this change was more modest than in lesional skin (Figure 1E–G), similar to that seen in the QuSAGE analysis (Figure 1D, Figure 1—figure supplement 1). IFN gene set changes in non-lesional skin were driven by IFN-α, IFN-κ, and IFN-ω gene set upregulation (Figure 1H). Non-lesional skin showed de-enrichment of the LC gene set and this effect was more modest than in lesional skin (Figure 1E and G), consistent with the noticeable loss of LCs in lesional skin compared to adjacent non-lesional skin by immunostaining of tissue sections (Sontheimer and Bergstresser, 1982). In contrast, non-lesional skin showed de-enrichment of the keratinocyte gene set, an effect not seen in lesional skin (Figure 1E and G), suggesting a distinct biology in non-lesional and lesional skin. GSVA, then, showed that while lesional skin was characterized by a proinflammatory IFN-rich environment and reduction of LC signals, non-lesional skin had a less inflammatory environment but one that was notable for an IFN response and LC reduction.
We next examined non-lesional ear skin from MRL/lpr, B6.Sle1yaa, and imiquimod (IMQ) SLE mouse models by RNAseq. The photosensitive MRL/lpr model has a Fas gene mutation and develops a lupus-like phenotype especially well on the MRL genetic background (Menke et al., 2008; Theofilopoulos and Dixon, 1985). MRL/lpr mice have reduced LC ADAM17 sheddase function and LC Adam17 mRNA, and topical EGFR ligand supplementation ameliorated photosensitivity (Shipman et al., 2018). Although the differences in gene expression between MRL/lpr mice and control MRL-MpJ (MRL/+) mice by RNAseq involved fewer genes than those in the human cohort (Figures 1A–D ,–2A–D; Supplementary file 4), QuSAGE pathway analysis showed that the IFN-α/β and IFN-γ pathways were among the most highly expressed in MRL/lpr mice compared to control MRL/+mice (Figure 2B and Supplementary file 5). Consistent with the activation of IFN pathways, IRF 1, 2, 7, and 9 transcription factors and a number of their targets were expressed at higher levels in MRL/lpr mice (Figure 2C–D and Figure 2—figure supplement 1A).
The B6.Sle1yaa lupus model is driven by the Sle1 lupus susceptibility locus derived from lupus-prone NZB2410 mice in combination with the Y chromosome autoimmune accelerator locus whose effects are attributed to Tlr7 duplication (Pisitkun et al., 2006; Subramanian et al., 2006). B6.Sle1yaa mice are photosensitive and have reduced LC Adam17 mRNA expression (Shipman et al., 2018). As expected, compared to B6 controls, B6.Sle1yaa mice expressed elevated levels of several lupus related genes from the Yaa locus including Tlr7 (Figure 2—figure supplement 1C). Although differences between the B6 and B6.Sle1yaa transcriptomes were modest (Figure 2E; Supplementary file 6), several genes involved in the regulation of and response to IFN were expressed at higher levels, including Irf7, Irf9, Isg15, Xaf1, Selp, and Ube2l6 (Figure 2—figure supplement 1D). LC numbers as a proportion of total cells were not different at 10 months of age when this model is fully diseased (Shipman et al., 2018), but Cd207 expression was downregulated in these younger (5 month old) diseased mice when compared to B6 controls (Figure 2—figure supplement 1D), potentially hinting at either reduced LC numbers at this time point or altered CD207 expression on LCs.
The IMQ model is induced by 4–5 weeks of topical IMQ, a TLR 7/8 agonist, resulting in autoantibody production, enlarged spleen and lymph nodes, nephritis, and photosensitivity characteristic of SLE (Ambler et al., 2022; Goel et al., 2020; Yokogawa et al., 2014). The model can be induced in B6 mice, allowing for the use of transgenic models to assess for effects on lupus features without crossing onto a different genetic background. We applied IMQ to only one ear of B6 mice and assessed gene expression by RNAseq in both the painted ipsilateral ear and the unpainted contralateral ‘non-lesional’ ear. Similar to non-lesional skin in the MRL/lpr mice, the contralateral ear in IMQ mice demonstrated an IFN-I signature, with upregulation of a number of ISGs and IRF 7 and 9 (Figure 2F–I, Figure 2—figure supplement 1B, Supplementary file 7, Supplementary file 8). In contrast to MRL/lpr mice, IFN gamma pathways were not upregulated in IMQ mice. By QuSAGE analysis, then, the non-lesional skin of MLR/lpr, B6.Sle1yaa, and IMQ SLE models showed similarities to human non-lesional skin in demonstrating an IFN-rich environment.
By GSVA using lupus-relevant murine gene sets (Supplementary file 9), MRL/lpr mice showed increased enrichment of a number of inflammatory cell and cytokine gene sets including T cells, myeloid cells, neutrophils, pDCs, IFN, IL1, IL21, and TNF (Figure 3A). The keratinocyte signature was enriched in MRL/lpr mice, in contrast to the de-enrichment of this gene set in human DLE (Figure 1G, Figure 3A). The LC signature did not show any quantitative changes; assuming that this gene set could in part reflect LC numbers, this result was consistent with unchanged LC numbers in the skin of MRL/lpr mice (Shipman et al., 2018). Similar to human DLE, the IFN gene set upregulation in MRL/lpr non-lesional skin was broad, with significant enrichment of gene signatures induced by stimulation from multiple IFNs (Figure 3B). B6.Sle1yaa mice showed less dramatic differences than the MRL/lpr model, but also showed upregulation of the IFN gene set without a particular contribution from any one IFN (Figure 3C–D). The LC gene set was de-enriched in B6.Sle1yaa mice (Figure 3C), driven in part by the reduction in Cd207 (Figure 2—figure supplement 1C), potentially reflecting decreased LC numbers or altered LC phenotypes. The contralateral ‘non-lesional’ ears from IMQ-treated mice also exhibited significant upregulation of the IFN gene signature and were similar to MRL/lpr mice in showing enriched gene signatures of myeloid cells, IL21, and TNF (Figure 3E). T cell genes were uniquely de-enriched in IMQ skin (Figure 3E). The LC gene set was not altered (Figure 3E), although there was de-enrichment in IMQ-painted ipsilateral skin (Figure 3—figure supplement 1). Similar to human DLE and MRL/lpr mice, the IFN signature in non-lesional IMQ skin reflected enrichment of multiple IFN-I gene sets (Figure 3F). Interestingly, the overall IFN gene set was not enriched in the IMQ-painted ipsilateral skin (Figure 3—figure supplement 1), although there was enrichment in two of the specific IFN signatures (IFNB1 and IFNW1; Figure 3—figure supplement 1). This suggested a distinct biology such as tissue damage in the skin that received repeated direct exposure to IMQ. Together, the gene expression analysis revealed unique transcriptomic profiles across multiple murine SLE models and pointed to a shared IFN signature across non-lesional skin of all models and human DLE.
We assessed for common upregulated genes in human DLE, MRL/lpr, and IMQ datasets. While a direct compassion was complicated by the difference in data acquisition methods (microarray versus RNAseq), species and the magnitude of response between mouse models and human patients, utilizing a common size effect (logFC = 0.5) and FDR cutoffs (<0.05) yielded a small group of upregulated genes common to all three datasets (Figure 4A). All of these genes (Sprr1b, Isg15, Ddx60, Bst2, Xaf1, Rsad2, Slc7a8, Ifit3) were interferon-inducible genes. Similarly, in pairwise set overlaps (DLE – MRL/lpr, MRL/lpr – IMQ or IMQ – MRL/lpr) large proportion of overlapping genes were IFN-inducible targets including Ifi genes (Ifi44, 203, 204), IRF transcription factors (Irf 1,2,7, 8, and 9), Oas 2,3 and L2, Mx1 (Supplementary file 10).
Comparisons of the GSVA analyses of non-lesional skin among the datasets (Figure 4B) also emphasized the prominent enrichment of IFN genes in human DLE and MRL/lpr, B6.Sle1yaa, and B6 IMQ SLE models. Notably, IFN upregulation correlated with loss of LCs in human DLE and B6.Sle1yaa mice but not in MRL/lpr and IMQ mice (Figure 4B). The lack of correlation in IMQ mice also extended to the painted ipsilateral ear, where there was a reduced LC signature in the absence of IFN upregulation (Figure 3—figure supplement 1). Overall, these data are consistent with the idea that, similar to human lupus skin, non-lesional skin in multiple SLE models is characterized by an IFN-rich environment that potentially primes the skin for photosensitive responses. These data also suggest that LCs sit within an IFN-rich environment in both human lupus and multiple murine models that may cause their dysfunction.
IFN-I inhibits LC ADAM17 sheddase function
We examined the effects of IFN-I on UVR-induced LC ADAM17-mediated shedding and also sought to assess the extent to which altered sheddase activity reflected altered LC ADAM17 protein expression. We had previously established a facs-based LC ADAM17-mediated shedding assay by quantifying UVR-induced cell surface TNFR receptor 1 (TNFR1) loss at 45 min after UVR exposure, extensively validating that this loss reflected shedding activity, paralleled shedding of EGFR ligands into the supernatant, and was dependent on LC ADAM17 (Shipman et al., 2018). As we had studied LC ADAM17 sheddase activity using isolated LCs, we confirmed here that UVR treatment of a mixture of digested epidermal cells also induced LC cell surface TNFR1 loss (Figure 5A). This effect was LC ADAM17-dependent (Figure 5A), validating UVR-induced LC cell surface TNFR1 loss as a readout of LC ADAM17 activity in the context of a mixture of epidermal cells.
Anti-murine ADAM17 antibodies to quantify cell surface ADAM17 have not been available until one recently described by our group (Lora et al., 2021). This antibody was able to detect cell surface ADAM17 on LCs and other cells digested from skin, and this signal was specific, as it was lost from LCs in cells from mice lacking LC ADAM17 (Figure 5B; Shipman et al., 2018). Using the validated LC ADAM17 sheddase assay and new anti-murine ADAM17, we proceeded to examine the role of IFN-I in regulating LC ADAM17 activity and expression.
IFN-κ, an IFN-I that is highly expressed by lupus keratinocytes (Tsoi et al., 2019), reduced LC ADAM17 sheddase activity by nearly 40% when used to treat LCs ex vivo (Figure 5C). In contrast, IFN-κ reduced LC ADAM17 protein levels only modestly by 8% and had no effect on cell viability (Figure 5D–E). The IFN-κ-driven inhibition of ADAM17 activity was abrogated by the JAK1/3 inhibitor tofacitinib (Figure 5F) without altering ADAM17 protein levels (Figure 5G) or cell viability (Figure 5H), supporting the idea that IFN-κ mainly affected LC ADAM17 by inhibiting sheddase activity. Another IFN-I, IFN-β, similarly reduced LC ADAM17 activity without effects on cell viability ex vivo (Figure 5I–J). In vivo, application of IFN-κ to the skin of wildtype B6 mice reduced LC ADAM17-mediated shedding by 57% (Figure 5K) while cell surface ADAM17 was reduced by 9% (Figure 5L). There was no effect on LC numbers (Figure 5M). The ex vivo and in vivo murine data together suggested that IFN-I downregulates UVR-induced LC ADAM17-mediated shedding mainly by inhibiting ADAM17 sheddase function rather than by reducing ADAM17 protein expression.
We also asked how IFN-I affects human LCs. While LCs can be obtained from skin biopsies or discarded tissues associated with surgical procedures, these tissues often require lengthy enzymatic digestion to dissociate the cells (Shipman et al., 2018), which has the potential to disrupt cell function. We found that epidermal roofs yielded from suction blistering could be digested for a short amount of time to obtain an epidermal cell mixture containing LCs (Figure 5—figure supplement 1; Strassner et al., 2017). Addition of IFN-β reduced LC ADAM17 activity by 62% without altering ADAM17 expression or LC viability (Figure 5N–P). Together, these results suggested that IFN-I can inhibit LC ADAM17 sheddase function in both murine and human LCs and supported the idea that high levels of IFN-I could potentially contribute to LC ADAM17 dysfunction in disease.
IFNAR is important for LC ADAM17 dysfunction in multiple lupus models
We asked whether IFN-I was an important regulator of LC ADAM17 function in lupus models. Consistent with previous findings (Shipman et al., 2018), MRL/lpr mice showed a 42% reduction in LC ADAM17 sheddase activity when compared to control MRL/+mice (Figure 6A). LC ADAM17 protein levels, however, were less affected, showing a 21% reduction (Figure 6B). Intraperitoneal (i.p.) anti-IFNAR treatment of MRL/lpr mice at 500 μg/dose restored UVR-induced LC TNFR1 loss (Figure 6A), suggesting restoration of LC ADAM17 sheddase activity. Anti-IFNAR did not alter cell surface ADAM17 protein expression (Figure 6B), suggesting that IFN-I regulated LC ADAM17 sheddase function rather than expression. A lower 125 μg/dose of anti-IFNAR had similar effects in restoring LC ADAM17 sheddase activity without affecting LC ADAM17 protein levels (Figure 6C–D). These results together suggested that IFN-I is important for inhibiting LC ADAM17 sheddase function in the MRL/lpr model.
LCs from the B6.Sle1yaa model showed a 63% reduction in UVR-induced TNFR1 loss when compared to control B6 mice without a corresponding change in cell surface ADAM17 protein levels (Figure 6E–F), suggesting reduced LC ADAM17 sheddase function in this SLE model. Similar to MRL/lpr mice, anti-IFNAR restored LC ADAM17 activity in B6.Sle1yaa mice without increasing ADAM17 protein expression (Figure 6E–F). These results suggest that IFN-I is an important inhibitor of LC ADAM17 sheddase function in the B6.Sle1yaa model.
Similar to the B6.Sle1yaa model, IMQ model mice showed reduced UVR-induced TNFR1 loss without alterations in cell surface ADAM17 protein levels when compared to controls (Figure 6G–H), suggesting reduced LC ADAM17 sheddase function. Similar to the MRL/lpr and B6.Sle1yaa models (Shipman et al., 2018) and consistent with the unchanged LC gene set by GSVA (Figure 3E), LC numbers were unchanged from controls in the IMQ model (Figure 6I). Anti-IFNAR restored LC ADAM17 sheddase activity in IMQ mice and did not alter LC numbers (Figure 6G,I). Together, these data suggesting that anti-IFNAR restored LC ADAM17 sheddase function in non-lesional skin across three distinct SLE models pointed toward a scenario whereby the IFN-I-rich environment in non-lesional lupus skin inhibits LC ADAM17 sheddase function.
We summarized the alterations in LC numbers, ADAM17 protein levels, and ADAM17 sheddase activity and effects of IFN-I modulation shown in Figures 5–6 and Shipman et al., 2018 to assess for consistent changes. While LC ADAM17 protein levels were not consistently affected and LC numbers were unchanged in the three lupus models examined, LC ADAM17 sheddase activity was consistently reduced (Table 1). Furthermore, manipulation of IFNAR signaling in these models, similar to treating human and murine LCs with IFN-I, consistently regulated only LC ADAM17 sheddase activity (Table 1). The uncoupling of IFN-I-mediated regulation of LC ADAM17 sheddase activity from LC numbers and LC ADAM17 expression suggests that IFN-I regulates LC ADAM17 sheddase function and points to an IFN-LC ADAM17 sheddase function axis as a potential mediator of IFN-I effects in skin.
Anti-IFNAR reduces photosensitivity in a manner dependent on EGFR signaling and LC ADAM17
We next asked whether the anti-IFNAR-mediated restoration of LC ADAM17 sheddase activity impacted photosensitive skin responses. Anti-IFNAR treatment at 500 μg/dose reduced ear swelling in MRL/lpr mice (Figure 7A–B) but did not reduce epidermal permeability, a readout of skin barrier function (Figure 7C). Surprisingly, in control MRL/+mice, anti-IFNAR actually increased UVR-induced ear swelling and epidermal permeability (Figure 7B–C), suggesting that anti-IFNAR at this dose-induced photosensitivity in non-lupus mice. These results echoed findings by other groups that IFNAR deficiency in non-lupus B6 mice induced photosensitivity (Sontheimer et al., 2017) and reduced skin wound healing (Gregorio et al., 2010; Wolf et al., 2022), and we considered the possibility that the 500 μg anti-IFNAR dose perhaps mimicked complete IFNAR deficiency and disrupted normal skin functions, thus obscuring the effects of LC ADAM17 activity restoration. Reducing anti-IFNAR to 125 μg/dose still increased UVR-induced MRL/+ear swelling but no longer increased epidermal permeability (Figure 7D–E), suggesting a decreased pathologic effect on physiologic skin functions. This lower dose, which was sufficient to restore LC ADAM17 activity in MRL/lpr mice (Figure 6C), was sufficient to reduce both UVR-induced skin swelling and epidermal permeability (Figure 7D–E) and limit neutrophil and monocyte infiltration after UVR exposure (Figure 7F). Thus, anti-IFNAR at 125 μg/dose that restored LC ADAM17 sheddase function also reduced UVR-induced skin swelling, epidermal permeability, and inflammatory infiltrates in MRL/lpr mice.
We sought to understand the extent to which the anti-IFNAR-driven reduction in photosensitivity reflected the rescue of LC ADAM17 activity. As LC ADAM17 limits photosensitivity by stimulating EGFR (Shipman et al., 2018), we blocked EGFR signaling upon anti-IFNAR treatment of MRL/lpr mice (Figure 7G). The small molecule EGFR tyrosine kinase inhibitor PD168393 (PD16) abrogated the anti-IFNAR-mediated reduction in ear swelling and epidermal permeability (Figure 7H–I), suggesting that anti-IFNAR requires EGFR signaling to ameliorate photosensitivity in MRL/lpr mice. These data are consistent with the idea that anti-IFNAR reduces photosensitivity by restoring LC ADAM17 sheddase function.
IMQ mice were characterized by monocyte accumulation while neutrophils were lower in number (Figure 7J–K). Anti-IFNAR reduced the monocyte cell infiltrate in UVR-treated skin and had no effect on neutrophil numbers (Figure 7K). LC numbers were unaffected (Figure 7L). In IMQ-treated LCAd17 mice lacking ADAM17 in LCs (Shipman et al., 2018), however, anti-IFNAR failed to reduce monocyte accumulation (Figure 7K), suggesting that the ameliorative effects of anti-IFNAR on UVR-induced inflammation are, in part, dependent on LC ADAM17. That anti-IFNAR reduction of photosensitive responses across two lupus models depended on EGFR signaling and LC ADAM17 provided strong support for the idea that IFN-I contributes to photosensitivity at least in part by promoting LC ADAM17 sheddase dysfunction.
IFN-I inhibits UVR-induced LC ROS expression
UVR has been shown to trigger ADAM17 sheddase activity in fibroblasts via ROS generation (Singh et al., 2009), and we asked if IFN-I modulation of ROS was associated with changes in LC ADAM17 activity. We observed the expected UVR-induced ROS generation in LCs from WT skin (Figure 8A; de Jager et al., 2017). Skin from IMQ-treated mice showed a reduction in UVR-induced LC ROS generation when compared to control skin, and anti-IFNAR treatment of IMQ mice restored UVR-driven ROS generation (Figure 8A). Assessment of cytoplasmic versus mitochondrial LC ROS showed that UVR exposure upregulated mainly cytoplasmic ROS (Figure 8B–C). These results suggested that IFN-I could potentially limit UVR-induced LC cytoplasmic ROS to inhibit LC ADAM17 sheddase function.
Discussion
In this study, we demonstrate that non-lesional skin in three different murine SLE models is similar to human lupus skin in expressing an IFN signature and show that IFN-I contributes to the inhibition of LC ADAM17 sheddase dysfunction, which leads to photosensitivity. These results suggest that the IFN-rich environment in even non-lesional lupus skin promotes photosensitivity by at least in part by causing LC ADAM17 dysfunction and raises the possibility that the beneficial effects of anifrolumab on skin disease in SLE patients reflects the restoration of LC ADAM17 sheddase function. Our results also suggest that photosensitivity in the MRL/lpr, B6.Sle1yaa, and IMQ lupus models shares pathogenic mechanisms with human lupus, underscoring the utility of these models for the study of lupus skin pathophysiology.
Our transcriptomic data both expand on the characterization of non-lesional DLE skin and provides new characterization of non-lesional skin in MRL/lpr, B6.Sle1yaa, and IMQ lupus models. Comprehensive transcriptomic analyses of non-lesional lupus model skin has not been available previously and our data serve as a resource for understanding shared or distinct pathways between human lupus and murine models or among the distinct lupus models. Lupus patients are heterogeneous, which in part contributes to the difficulties in clinical trials (Allen et al., 2021). Understanding how specific lupus models can be used for studying different aspects of human lupus will help to advance disease insights and precision medicine approaches, especially as we consider how different tissues communicate with and regulate one another in disease.
Our data here, and especially the experiments showing the importance of EGFR signaling and LC ADAM17 in the ameliorative effects of anti-IFNAR on photosensitive responses, suggest that LC ADAM17 sheddase inhibition is a mechanism by which the IFN-rich environment primes non-lesional skin for photosensitive responses. This mechanism does not preclude parallel mechanisms such as potentiation of IL-6 expression and apoptosis of keratinocytes and activation of multiple other cell tyoes shown in the human system (Sarkar et al., 2018; Stannard et al., 2017; Billi et al., 2022). Given that LCs are positioned among the suprabasal layer of the epidermis near basal and spinous keratinocytes that demonstrate high IFN expression and signaling (Billi et al., 2022; Psarras et al., 2020; Sarkar et al., 2018), keratinocytes are likely to be an important source of IFN-I that causes LC ADAM17 sheddase dysfunction. This possibility coupled with our data here point to a potential pathogenic feed-forward loop between LCs and keratinocytes in lupus skin whereby keratinocyte-derived IFN-I causes LC ADAM17 sheddase dysfunction and this dysfunction contributes to greater UVR-induced keratinocyte damage, skin inflammation, and, over time, lesion formation.
Our findings suggest the possibility that IFNs have distinct functions in non-lesional and lesional skin. Physiologically, skin upregulates IFN-I with injury and infection, and the IFN is important in part for tissue repair and wound healing (Gregorio et al., 2010; Wolf et al., 2022; Zhang et al., 2016). This may explain why anti-IFNAR treatment of control MRL/+mice with physiologic IFN-I levels increased UVR-induced skin inflammation and suggests that the IFN-I in non-lesional skin in mic and humans reflects a repair response to a mild insult such as immune complex deposition that forms the ‘lupus band’ found in non-lesional human lupus skin. The upregulated IFN in non-lesional skin, by causing LC ADAM17 dysfunction or heightened keratinocyte responses, serves to ‘prime’ the skin for further injury and a second insult such as UVR manifests as a photosensitive response (Billi et al., 2022). Because the repair response is compromised, IFN signaling continues to be upregulated, further inhibiting repair responses and promoting inflammation, tissue damage processes, and lesion formation. In this model, then, IFN in non-lesional skin primes for further injury, and the (higher) IFN in lesional skin reflecting frustrated tissue repair, continues to prime but also contributes to tissue damage. IFNAR blockade, then, should both prevent new lesion formation and treat existing lesions as long as physiologic IFN levels for wound healing are preserved. Understanding the regulation and sources of the excess IFN-I in non-lesional and lesion skin will help to further develop therapies that ameliorate disease.
Our results extend the understanding of ADAM17 dysregulation as a contributor to disease. Increased ADAM17 activity can promote inflammation, such as by release of TNF or EGFR ligands within inflamed joints or lupus kidneys (Issuree et al., 2013; Qing et al., 2018); yet loss of ADAM17 can compromise barrier surfaces (Blaydon et al., 2011; Franzke et al., 2012). Here, we identify IFN-I as a negative regulator of LC ADAM17 sheddase function, and our data suggest a scenario whereby IFN-I reduces LC ROS generation to inhibit ADAM17 function. Further investigation will better delineate the importance of LC ROS generation as well as contributions of other molecules such as the inactive Rhomboids (iRhoms) that regulate ADAM17 activity (Geesala et al., 2019).
This report importantly distinguishes among IFN-I regulation of LC ADAM17 sheddase function, LC ADAM17 surface molecule expression, and LC numbers. Our data here showing that the regulation of LC ADAM17 sheddase activity by IFN-I is uncoupled from LC ADAM17 protein expression and LC numbers (summarized in Table 1) supports the idea that LC ADAM17 sheddase function dysregulation is a key contributor to the propensity for photosensitivity and is thus a therapeutic target. The uncoupling of ADAM17 sheddase activity from protein expression is consistent with the known biology of ADAM17 on other cell types, where sheddase function is triggered by many stimuli on other cell types including UVR and G-protein-coupled receptors (GPCRs) and controlled by mediators such as ROS and iRhoms, as we and others have shown (Li et al., 2015; Zunke and Rose-John, 2017). We speculate that the ability to modulate ADAM17 enzymatic function offers the ability to regulate UVR responses in a more timely manner and with a much greater and tunable dynamic range than control of Adam17 transcription and translation or LC numbers. Targeting the effects of IFN-I on LC ADAM17 sheddase function, then, offers the opportunity to have a large impact in protecting the skin and systemic complications.
A future direction will be to further understand LC functional phenotypes in human lupus and murine models. Our data along with recent reports indicates that LC loss is a feature of human lupus skin (Billi et al., 2022), with the effect being variable in non-lesional skin (Martínez et al., 2022; Shipman et al., 2018; Sontheimer and Bergstresser, 1982; this report) and a greater more consistent effect in lesional skin (Bos et al., 1986; Martínez et al., 2022; Sontheimer and Bergstresser, 1982). On the other hand, LC numbers appeared less affected in lupus mouse models, with reduction of the LC signature in the 5-month-old B6.Sle1yaa non-lesional skin and in the IMQ painted ear of IMQ mice, but otherwise with normal LC numbers in non-lesional skin of MRL/lpr, 10-month-old B6.Sle1yaa (Shipman et al., 2018) and IMQ mice (this report). Notably, the LC numbers and LC signature in non-lesional IMQ skin, despite being statistically unchanged from that of controls, is variable, suggesting LC loss in some mice and bearing similarities to the variability in human non-lesional skin. These data suggest similar trends between human and mouse models in LC loss, with mouse models perhaps being less severely affected. The extent to which LC ADAM17 sheddase dysfunction represents a milder or just a different dysfunction than loss of LCs, the regulation of LC numbers via regulation of survival, migration, replenishment, the consequences to skin and immune function, and similarities and differences between LCs in human disease and murine systems awaits further study.
High IFN-I states are also seen in viral infections, genetic interferonopathies, and other autoimmune diseases including rheumatoid arthritis and dermatomyositis (Muskardin and Niewold, 2018) in which photosensitivity is found in the former and is a disease hallmark in the latter (Werth et al., 2004; Wysenbeek et al., 1989). Autoimmune disease-like manifestations that include photosensitivity are among the immune-related adverse events in cancer patients treated with checkpoint inhibitors (Brahmer et al., 2018), a situation where distinct IFN-I-mediated processes may have specific roles in promoting anti-tumor immunity and treatment resistance (Benci et al., 2019). Our findings and the work that will follow will benefit patients with lupus and a wide spectrum of diseases.
Materials and methods
Study design
Request a detailed protocolThe purpose of the study was to understand how LC function is regulated by the environment in non-lesional skin to contribute to photosensitivity. The subjects are both humans and mice. We used gene expression analyses to understand the gene expression profiles in human lupus and murine lupus models. We used ex vivo flow cytometry-based assays, assessment of epidermal permeability, and skin swelling assays to understand the effects of IFN on LC function and implications for photosensitivity. Sample sizes and replicate numbers are provided in each figure legend. Sample sizes were determined based on previous similar types of experiments. Mice were randomly assigned to treatment groups. No criteria were set for excluding certain data points and no data points were excluded. No specific confounder variables such as order of treatments or cage location were controlled for except that our experiments used multiple independent cohorts of mice over replicate experiments. Assessments and analyses were not blinded. ARRIVE reporting guidelines Percie du Sert et al., 2020 have been used.
Mice
Eight- to 16-week-old mice were used unless otherwise specified. Mice were sex and age-matched. For C57Bl/6 mice, LCAd17 (Langerin-cre; Adam17 f/f) mice (Shipman et al., 2018), and wildtype littermate controls, males and females were used. Unless otherwise indicated, for MRL/+and MRL/lpr mice, females were used at 8–12 weeks, while males were used at 12–16 weeks as females present with a disease phenotype earlier than males. For B6.Sle1yaa mice and their B6 controls, only males were used because the model depends on the autoimmune accelerator locus located on the Y chromosome. B6 mice were either bred in-house or purchased from Jackson Laboratory (Bar Harbor, ME). MRL/+, MRL/lpr mice, B6.Sle1yaa mice, and their B6 controls were either bred in-house or purchased from Jackson Laboratory. LCAd17 mice were bred in-house. Mice were kept in a specific pathogen-free barrier facility and all animal procedures including protocols to reduce pain, suffering, and distress were performed in accordance with the regulations of the Institutional Animal Care and Use Committee at Weill Cornell Medicine (New York, NY; protocol number 2015–0067).
Mouse treatments
Request a detailed protocolFor IFN-κ painting, 37.5 ng IFN-κ (R&D Systems, Minneapolis, MN) was dissolved in DMSO (Sigma-Aldrich, St. Louis, MO) was applied topically onto each ear (18.75 ng on dorsal and ventral sides, each). For anti-IFNAR1 treatments, mice were injected IP with 500 µg/dose or 125 µg/dose of anti-IFNAR1 (MAR1-5A3) or isotype control (MOPC-21; Bio X Cell, West Lebanon, NH). For UVR treatments, mice were exposed to 1000 J UVB/m2 (100 mJ/cm2) daily for 4 consecutive days using a bank of four FS40T12 sunlamps that emitted a combination of UVA and UVB radiation, similar to described (Shipman et al., 2018). Timing of UVR exposure was made based on lamp emission of 1.31 mW/cm2 UVB and distance of 9 cm away from lamp. To measure ear swelling after UVR exposure, a caliper (Mitutoyo, Aurora, IL) was used to take three measurements throughout each ear, and the reported value is the average of both ears. Each measurement was taken 22–24 hr after the previous measurement. For the IMQ-induced lupus mouse model, mice were painted on the dorsal and ventral sides of the right (ipsilateral) ear with 5% imiquimod cream (42 mg/ear/mouse) 3 x/week for 4–6 weeks (Yokogawa et al., 2014) prior to, and continued through experiments. The unpainted (contralateral) left ear was taken as non-lesional skin for flow cytometry experiments. For RNAseq, the painted ipsilateral ear and the unpainted contralateral ‘non-lesional’ left ear were taken and analyzed separately. Indicated MRL/lpr mice received 2.95 mg of the irreversible EGFR inhibitor PD168393 (0.74 mg on dorsal and ventral sides, each; Selleck Chem, Pittsburgh, PA).
Human subjects
Request a detailed protocolFor microarray analysis, seven discoid CLE patients were examined. Both lesional and non-lesional skin were collected and the gene expression data from lesional tissue of these patients have been published, and patient characteristics are described (Jabbari et al., 2014). Healthy controls include the three that were previously described (Jabbari et al., 2014). All human tissue collection and research use adhered to protocols approved by the Institutional Review and Privacy Boards at the Rockefeller University (IRB# AJA-00740) and New York University (IRB# 10–02117), where participants signed written informed consents.
For suction blistering, one male and three female healthy participants between the ages of 23 and 55 participated. All human tissue collection and research use adhered to protocols approved by the Institutional Review and Privacy Board at the Hospital for Special Surgery (IRB# 2019–1998), where participants signed written informed consents.
Human skin cell collection and preparation
Request a detailed protocolSuction blistering was performed as previously described (27). Eight 5 mm suction blisters were generated with a negative pressure instrument (NP-4 model, Electronic Diversities, Finksburg, MD) on the arm of healthy donors over 30–60 min. The blister fluid was then collected by aspiration, and epidermal roofs were collected. After collection, epidermal roofs from suction blisters were digested with dispase (2.42 U/mL; Thermo Fisher Scientific, Waltham, MA), collagenase type II (616 U/mL; Worthington Biochemical Corporation, Lakewood, NJ), and DNAseI (80 μg/mL; Sigma Aldrich, St. Louis, MO) for 40 min to generate single cell suspensions that were used for the ADAM17 sheddase assay.
Skin cell collection and flow cytometry staining and quantification
Request a detailed protocolFor murine epidermal single cell suspensions, skin on the trunk was excised and subcutaneous fat scraped off. Skin was incubated in dispase (2.42 U/mL) at 37 °C. The epidermis was then gently peeled, finely minced, digested in type II collagenase (616 U/mL) and DNAse I (80 μg/mL).
For collection of ear skin, ears were cut along the cartilage ridge, and the dorsal and ventral sides were manually peeled. Ear skin was then finely minced and digested with a buffer containing dispase, type II collagenase, and DNAse I as described above. For murine LC TNFR1 staining for the ADAM17 sheddase activity assay (see below), epidermal cells were preincubated with Fc block and then stained so that we could gate on CD45+, CD3-, CD11c+, I-Ab + for B6 mice and I-Ak + for MRL mice, and EpCAM + LCs and assess for TNFR1. LCs were sometimes stained with rabbit-anti-ADAM17 (Lora et al., 2021) followed by donkey anti-rabbit IgG (Jackson Immunoresearch, West Grove, PA). Monocytes and T cells were identified by gating for CD45+, CD11b+, Ly6C+, and CD45+, CD11c-, I-Ak/I-Ab-, and CD3 + cells, respectively. 4’,6-diamidino-2-phenylindole dihydrochloride (DAPI) was used to exclude dead cells and debris (Sigma Aldrich, St. Louis, MO). All antibodies are from Biolegend, San Diego, CA unless otherwise indicated.
For human LC TNFR1 staining for the ADAM17 sheddase activity assay, single cell suspensions of the epidermal roofs of suction blisters and blister fluid were preincubated with Fc block and subsequently stained for CD45+, CD11c+, CD14-, CD3- CD19-, CD56-, CD66b-, EpCAM+, and HLA-DR +to identify LCs. DAPI was used to exclude dead cells and debris, and TNFR1 was also stained to assess ADAM17 activity. All antibodies are from Biolegend, San Diego, CA unless otherwise indicated. Cells were analyzed using a FACS Canto or FACS Symphony A3 (BD Biosciences, San Jose, CA) and data were analyzed with FlowJo V10 Software (TreeStar, Ashland, OR).
For antibodies used, please see Supplementary file 11.
Ex vivo ADAM17 sheddase activity assay
Request a detailed protocolFor murine skin, cell suspensions of epidermis were suspended at 200,000–500,000 cells/mL medium. For MRL/LPR and B6.Sle1yaa skin, medium was RPMI with 2% FBS, HEPES buffer, L-glutamine, penicillin-streptomycin (all Thermo Fisher Scientific, Waltham, MA). For IMQ mice skin, medium was phenol red-free RPMI (Corning, Edison, New Jersey), no serum, with HEPES buffer, L-glutamine, penicillin-streptomycin. Cells were exposed to 2000 (LPR and B6.Sle1yaa model mice) or 1000 (IMQ model mice) J/m2 UVR incubated for 45 min, and then stained for cell surface markers described above. For IFN-κ treatment, 3.125 ng/mL murine IFN-κ (R&D Systems, Minneapolis, MN) was added to the epidermal cell suspension that were incubated for 16–20 hr prior to UVR exposure. Indicated cells were treated with 1 nM tofacitinib (Selleck Chem, Pittsburgh, PA) during IFN-κ incubation. For IFN-β treatment, 30 ng/ml murine IFN-β (R&D Systems) was added to the epidermal cell suspensions and assessed as above.
For human skin, cells were incubated for 1 hr with recombinant human IFN-β (35.7 ng/mL; activity = 2.8 × 108 UI/mg) and then exposed to 1000 J/m2 UVR as above. Cells were further incubated for 45 min and collected for preincubation with Fc block. Cells were then stained for cell surface markers described above.
Calculating the relative change in cell surface TNFR1 levels upon UVR exposure was as previous (Shipman et al., 2018). The geometric MFI of TNFR1 was divided by the MFI for the isotype control to generate a relative MFI. The percent TNFR1 MFI change was calculated by setting the relative MFI of the control sample to 1 and normalizing the value from the experimental samples to that control. In experiments with more than one control sample, the control values were averaged, and the experimental samples were each normalized to the averaged value.
Epidermal permeability assay
Request a detailed protocolEpidermal barrier function was assessed as described (Shipman et al., 2018). Ear skin was dehydrated and rehydrated in graded methanol and then incubated in 0.1% toluidine blue dye (Sigma Aldrich, St. Louis, MO). The dye was extracted with a buffer consisting of 95% methanol, 2.5% sulfuric acid, and 2.5% water and quantified with a colorimeter at 620 nm. Dye concentration was calculated using a standard curve.
Ex vivo ROS assays
Request a detailed protocolEpidermal cell suspensions in serum-free, phenol red-free DMEM/F-12 media were exposed to 1000 J/m2 UVR and the general ROS indicator 5-(and-6)-chloromethyl-2’,7’-dichlorodihydrofluorescein diacetate, acetyl ester (CM-H2DCFDA) (Invitrogen, Waltham, MA) at 10 uM was added immediately after UVR exposure. CM-H2DCFDA oxidation by ROS inside the cells leads to fluorescence. Cells were further incubated for 45 min prior to cell surface marker staining to identify LCs by flow cytometry.
To assess cytoplasmic vs mitochondrial sources of ROS, epidermal cell suspensions were incubated with MitoSOX Green superoxide indicator (Thermo Fisher) and CellROX Deep Red Reagent (Thermo Fisher; final concentration of 5 µM for each) or vehicle (DMF and DMSO, respectively). Fifteen minutes after addition of reagents, cells were exposed to 750 J/m2 UVR, incubated another 15 min, and then washed three times with warm PBS prior to staining for extracellular surface markers for flow cytometry.
RNA extraction
Request a detailed protocolFemale MRL/+and MRL/lpr 12–14 weeks old, 12-week-old B6 and B6.Sle1yaa male mice, and 12- to 14-week-old IMQ mice were used for RNA sequencing. RNA was extracted from ear skin using an RNEasy Kit (Qiagen, Germantown, MD) and quality confirmed on a BioAnalyzer 2100 (Agilent Technologies, Santa Clara, CA).
RNA sequencing
Request a detailed protocolLibraries were prepared with non-stranded and poly-A selection TruSeq RNA Sample Preparation kits. Sequencing libraries were sequenced by the Genomics Core Facility at Weill Cornell Medicine using a HiSeq 2500 (MRL/lpr, MRL/+samples) and NovaSeq 6000 (samples from B6.Sle1yaa and IMQ models and their controls) at a depth of 15–30 million mappable paired-end reads per sample.
Computational analyses of gene expression data
Request a detailed protocolMicroarray analysis was performed using healthy controls and DLE lesional skin samples from an existing GEO dataset (GSE52471) obtained from GEO using GEOquery and combined with previously unpublished, unsubmitted non-lesional skin sample data obtained as part of the original study (Jabbari et al., 2014). Microarray data were normalized with gcRMA normalization using the gcrma R package (Wu et al., 2021). Differential gene expression analysis was performed with limma using the empirical bayes method (Ritchie et al., 2015) and controlling for microarray hybridization kit and sex. Batch correction was performed using ComBat (Kammers et al., 2021). Visualizations were generated using ComBat corrected data with plotly and ggplot2 in R.
For RNA sequencing analysis, read quality filtering and adapter trimming were performed using fastp (Chen et al., 2018). Filtered reads were mapped to the mouse genome (mm10), and exonic reads were counted against GENCODE release 27 annotation with the STAR aligner (Dobin et al., 2013) using a customized pipeline available at https://gitlab.com/hssgenomics/pipelines (Oliver, 2023). Differential gene expression analysis was performed with edgeR 3.30.2 and 3.38.4 using quasi-likelihood framework. Genes with low expression levels (<2 counts per million in at least one group) were filtered from all downstream analyses. p-Values were adjusted to correct for multiple testing using the Benhamini-Hochberg FDR procedure. Genes with adjusted p-values <0.05 were considered differentially expressed. Downstream analyses were performed in R and visualized using a Shiny-driven visualization platform (RNAseq DRaMA) developed at the HSS David Z. Rosensweig Genomics Research Center.
Differentially regulated pathways were identified using R QuSAGE 2.22.0 package (Yaari et al., 2013) with MSigDB C2 set (curated gene sets). All gene sets with less than 10 genes were excluded from the analysis. Pathways with Benhamini-Hochberg adjusted p-values <0.005 (non-lesional skin versus healthy control), <0.05 (MRL/lpr versus MRL/+, IMQ contralateral versus healthy B6) or unadjusted p-values <0.05 (Sle.yaa vs B6) were used for further analyses.
The R/Bioconductor package, GSVA (Hänzelmann et al., 2013) (v1.36.3), was utilized as a non-parametric, unsupervised method to estimate enrichment of pre-defined gene sets in gene expression data from human DLE skin and non-lesional skin from MRL/lpr, B6.Sle1yaa, and IMQ mice. The inputs for GSVA were a matrix of expression values for all samples and curated gene sets describing select immune/tissue cell types and inflammatory cytokines. For the analysis of human DLE microarray data, low-intensity probes were filtered out if the interquartile range (IQR) of their expression values across all samples was not greater than 0 prior to the analysis. Enrichment scores were calculated using a Kolgomorov Smirnoff (KS)-like random walk statistic and represented the greatest deviation from zero for a particular sample and gene set. Scores across all samples were normalized to values between –1 (indicating no enrichment) and +1 (indicating enrichment). Gene sets used as input for GSVA are listed inSupplementary files 3 and 9. The human gene sets were previously published (Martínez et al., 2022). The mouse plasma cell, T cell, myeloid cell, neutrophil, pDC, dendritic cell, and keratinocyte gene sets were derived from Mouse CellScan, a tool for the identification of cellular origin of mouse gene expression datasets. The mouse IFN gene sets have been previously described (Kingsmore et al., 2021). The mouse cytokine signatures were generated by an iterative process involving derivation through literature mining and GO term definitions provided by the Mouse Genome Informatics (MGI) GO Browser (Bult et al., 2019).
Statistical analyses
Request a detailed protocolWe determined the normality of data distribution using the Shapiro-Wilk test. For normally distributed data, we used two-tailed unpaired and paired Student’s t-tests for comparisons between two conditions to evaluate p-values as indicated. For data that were not normally distributed, the Mann-Whitney test was used for unpaired and Wilcoxon matched-pairs signed rank test for paired comparisons. Significance was defined as p<0.05. GraphPad Prism software was used. For figures showing normalized values, the control sample was set to one, and the experimental samples were normalized relative to the control for that experiment. For experiments that contained more than one control sample, the mean was obtained for the control samples, and the individual control and experimental samples were calculated relative to this mean.
Data availability
The murine RNAseq data have been deposited in GEO (GSE222573 for MRL/lpr and B6.Sle1yaa mice; GSE255519 for IMQ mice). Non-lesional human DLE microarray data were deposited with accession number GSE227329 with reanalysis of healthy control samples from GSE52471. All other data supporting the findings of this study are available within the paper and its Supplementary Materials.
-
NCBI Gene Expression OmnibusID GSE222573. The interferon-rich skin environment regulates Langerhans cell ADAM17 to promote photosensitivity in lupus.
-
NCBI Gene Expression OmnibusID GSE255519. The interferon-rich skin environment regulates Langerhans cell ADAM17 to promote photosensitivity in lupus.
-
NCBI Gene Expression OmnibusID GSE227329. The interferon-rich skin environment regulates Langerhans cell ADAM17 to promote photosensitivity in lupus [array].
-
NCBI Gene Expression OmnibusID GSE52471. Dominant Th1 and Minimal Th17 Skewing in Discoid Lupus Revealed by Transcriptomic Comparison with Psoriasis.
References
-
Leveraging heterogeneity in systemic lupus erythematosus for new therapiesTrends in Molecular Medicine 27:152–171.https://doi.org/10.1016/j.molmed.2020.09.009
-
Prevalence of mood and anxiety disorders in women with systemic lupus erythematosusArthritis and Rheumatism 61:822–829.https://doi.org/10.1002/art.24519
-
Interferon and granulopoiesis signatures in systemic lupus erythematosus bloodThe Journal of Experimental Medicine 197:711–723.https://doi.org/10.1084/jem.20021553
-
Nonlesional lupus skin contributes to inflammatory education of myeloid cells and primes for cutaneous inflammationScience Translational Medicine 14:eabn2263.https://doi.org/10.1126/scitranslmed.abn2263
-
Inflammatory skin and bowel disease linked to ADAM17 deletionThe New England Journal of Medicine 365:1502–1508.https://doi.org/10.1056/NEJMoa1100721
-
Mouse genome database (MGD) 2019Nucleic Acids Research 47:D801–D806.https://doi.org/10.1093/nar/gky1056
-
fastp: an ultra-fast all-in-one FASTQ preprocessorBioinformatics 34:i884–i890.https://doi.org/10.1093/bioinformatics/bty560
-
Ultraviolet light induced generation of reactive oxygen speciesAdvances in Experimental Medicine and Biology 996:15–23.https://doi.org/10.1007/978-3-319-56017-5_2
-
STAR: ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
Mechanisms of photosensitivity in autoimmunityThe Journal of Investigative Dermatology 142:849–856.https://doi.org/10.1016/j.jid.2021.05.007
-
Prevalence of self-report photosensitivity in cutaneous lupus erythematosusJournal of the American Academy of Dermatology 66:220–228.https://doi.org/10.1016/j.jaad.2010.12.006
-
Characterization of clinical photosensitivity in cutaneous lupus erythematosusJournal of the American Academy of Dermatology 69:205–213.https://doi.org/10.1016/j.jaad.2013.03.015
-
Epidermal ADAM17 maintains the skin barrier by regulating EGFR ligand-dependent terminal keratinocyte differentiationThe Journal of Experimental Medicine 209:1105–1119.https://doi.org/10.1084/jem.20112258
-
Novel functions of inactive rhomboid proteins in immunity and diseaseJournal of Leukocyte Biology 106:823–835.https://doi.org/10.1002/JLB.3VMR0219-069R
-
Plasmacytoid dendritic cells sense skin injury and promote wound healing through type I interferonsThe Journal of Experimental Medicine 207:2921–2930.https://doi.org/10.1084/jem.20101102
-
iRHOM2 is a critical pathogenic mediator of inflammatory arthritisThe Journal of Clinical Investigation 123:928–932.https://doi.org/10.1172/JCI66168
-
Dominant Th1 and minimal Th17 skewing in discoid lupus revealed by transcriptomic comparison with psoriasisThe Journal of Investigative Dermatology 134:87–95.https://doi.org/10.1038/jid.2013.269
-
Photosensitivity in cutaneous lupus erythematosusPhotodermatology, Photoimmunology & Photomedicine 29:4–11.https://doi.org/10.1111/phpp.12018
-
Coordinate overexpression of interferon-alpha-induced genes in systemic lupus erythematosusArthritis and Rheumatism 50:3958–3967.https://doi.org/10.1002/art.20798
-
Targeted truncation of the ADAM17 cytoplasmic domain in mice results in protein destabilization and a hypomorphic phenotypeThe Journal of Biological Chemistry 296:100733.https://doi.org/10.1016/j.jbc.2021.100733
-
Sunlight triggers cutaneous lupus through a CSF-1-dependent mechanism in MRL-Fas(lpr) miceJournal of Immunology 181:7367–7379.https://doi.org/10.4049/jimmunol.181.10.7367
-
Trial of anifrolumab in active systemic lupus erythematosusThe New England Journal of Medicine 382:211–221.https://doi.org/10.1056/NEJMoa1912196
-
Type I interferon in rheumatic diseasesNature Reviews. Rheumatology 14:214–228.https://doi.org/10.1038/nrrheum.2018.31
-
Exacerbated autoimmunity in the absence of TLR9 in MRLThe Journal of Immunology 190:3889–3894.https://doi.org/10.4049/jimmunol.1203525
-
Emerging concepts of type I interferons in SLE pathogenesis and therapyNature Reviews. Rheumatology 18:575–590.https://doi.org/10.1038/s41584-022-00826-z
-
iRhom2 promotes lupus nephritis through TNF-α and EGFR signalingThe Journal of Clinical Investigation 128:1397–1412.https://doi.org/10.1172/JCI97650
-
Photosensitivity and type I IFN responses in cutaneous lupus are driven by epidermal-derived interferon kappaAnnals of the Rheumatic Diseases 77:1653–1664.https://doi.org/10.1136/annrheumdis-2018-213197
-
A protective Langerhans cell-keratinocyte axis that is dysfunctional in photosensitivityScience Translational Medicine 10:eaap9527.https://doi.org/10.1126/scitranslmed.aap9527
-
Immune cell-stromal circuitry in lupus photosensitivityJournal of Immunology 206:302–309.https://doi.org/10.4049/jimmunol.2000905
-
UV-induced EGFR signal transactivation is dependent on proligand shedding by activated metalloproteases in skin cancer cell linesInternational Journal of Cancer 124:531–539.https://doi.org/10.1002/ijc.23974
-
Epidermal Langerhans cell involvement in cutaneous lupus erythematosusThe Journal of Investigative Dermatology 79:237–243.https://doi.org/10.1111/1523-1747.ep12500069
-
Lupus Skin Is Primed for IL-6 inflammatory responses through a keratinocyte-mediated autocrine typeiinterferon loopThe Journal of Investigative Dermatology 137:115–122.https://doi.org/10.1016/j.jid.2016.09.008
-
Suction blistering the lesional skin of vitiligo patients reveals useful biomarkers of disease activityJournal of the American Academy of Dermatology 76:847–855.https://doi.org/10.1016/j.jaad.2016.12.021
-
Murine models of systemic lupus erythematosusAdvances in Immunology 37:269–390.https://doi.org/10.1016/s0065-2776(08)60342-9
-
Photosensitivity in rheumatic diseasesThe Journal of Investigative Dermatology 9:57–63.https://doi.org/10.1111/j.1087-0024.2004.00839.x
-
SoftwareGcrma: background adjustment using sequence information, version 2.64.0R Package.
-
Prevalence and expression of photosensitivity in systemic lupus erythematosusAnnals of the Rheumatic Diseases 48:461–463.https://doi.org/10.1136/ard.48.6.461
-
The shedding protease ADAM17: Physiology and pathophysiologyBiochimica et Biophysica Acta 1864:2059–2070.https://doi.org/10.1016/j.bbamcr.2017.07.001
Article and author information
Author details
Funding
Alpha Omega Alpha Honor Medical Society Carolyn L. Kuckein fellowship
- Thomas Morgan Li
HSS Medical Student Summer Research Fellowhship
- Thomas Morgan Li
Erwin Schrodinger Fellowship (J 4638-B FWF)
- Victoria Zyulina
NIH (T32AR071302)
- Noa Schwartz
- William D Shipman
NIH MSTP grant (T32GM007739)
- William D Shipman
Tow Foundation
- Yurii Chinenov
- David J Oliver
NIH (K08 AR069111)
- Ali Jabbari
Veterans Administration VA Merit (I01 BX004907)
- Ali Jabbari
Dermatology Foundation Physician Scientist Career Development Award
- Ali Jabbari
NIH (R01AR077194)
- Ali Jabbari
NIH (DK099087)
- Inez Rogatsky
NIH (R35GM134907)
- Carl P Blobel
NIH (R01AI079178)
- Theresa T Lu
NIH (R21 AR081493)
- Theresa T Lu
DOD (W81XWH-21-LRP-IPA)
- Theresa T Lu
Lupus Research Alliance
- Theresa T Lu
St. Giles Foundation
- Theresa T Lu
Barbara Volcker Center for Women and Rheumatic Diseases
- Theresa T Lu
A Lasting Mark Foundation
- Theresa T Lu
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors thank the members of the Lu Lab for helpful discussions. Alpha Omega Alpha Carolyn L Kuckein Student Research Fellowship (TML). Erwin Schrodinger Fellowship J 4638-B FWF (VZ). National Institutes of Health grant T32AR071302-01 to the Hospital for Special Surgery Research Institute Rheumatology Training Program (NS and WDS). National Institutes of Health grant MSTP T32GM007739 to the Weill Cornell/Rockefeller/Sloan-Kettering Tri-Institutional MD-PhD Program (WDS). Tow Foundation (YC and DJO). National Institutes of Health grant K08 AR069111 to the University of Iowa Department of Dermatology (AJ). Veterans Administration VA Merit I01 BX004907 (AJ), Physician Scientist Career Development Award from the Dermatology Foundation (AJ). National Institutes of Health grant R01AR077194 (AJ). National Institutes of Health grant R01 DK099087 (IR). National Institutes of Health grant R35GM134907 (CPB). National Institutes of Health grant R01AI079178 (TTL). National Institutes of Health grant R21 AR081493 (TTL). Department of Defense grant W81XWH-21-LRP-IPA (TTL). Lupus Research Alliance (TTL). St. Giles Foundation (TTL). Barbara Volcker Center for Women and Rheumatic Diseases (TTL). A Lasting Mark Foundation (TTL). National Institutes of Health Office of the Director grant S10OD019986 to Hospital for Special Surgery.
Ethics
For samples used in microarray analysis - human tissue collection and research use adhered to protocols approved by the Institutional Review and Privacy Boards at the Rockefeller University (IRB# AJA-00740) and New York University (IRB# 10-02117), where participants signed written informed consents. For samples used in suction blistering - human tissue collection and research use adhered to protocols approved by the Institutional Review and Privacy Board at the Hospital for Special Surgery (IRB# 2019-1998), where participants signed written informed consents.
All animal procedures were performed in accordance with the regulations of the Institutional Animal Care and Use Committee at Weill Cornell Medicine (New York, NY) (Protocol number 2015-0067).
Copyright
© 2024, Li et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 788
- views
-
- 147
- downloads
-
- 1
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Immunology and Inflammation
During thymic development, most γδ T cells acquire innate-like characteristics that are critical for their function in tumor surveillance, infectious disease, and tissue repair. The mechanisms, however, that regulate γδ T cell developmental programming remain unclear. Recently, we demonstrated that the SLAM/SAP signaling pathway regulates the development and function of multiple innate-like γδ T cell subsets. Here, we used a single-cell proteogenomics approach to identify SAP-dependent developmental checkpoints and to define the SAP-dependent γδ TCR repertoire in mice. SAP deficiency resulted in both a significant loss of an immature Gzma+Blk+Etv5+Tox2+ γδT17 precursor population and a significant increase in Cd4+Cd8+Rorc+Ptcra+Rag1+ thymic γδ T cells. SAP-dependent diversion of embryonic day 17 thymic γδ T cell clonotypes into the αβ T cell developmental pathway was associated with a decreased frequency of mature clonotypes in neonatal thymus, and an altered γδ TCR repertoire in the periphery. Finally, we identify TRGV4/TRAV13-4(DV7)-expressing T cells as a novel, SAP-dependent Vγ4 γδT1 subset. Together, the data support a model in which SAP-dependent γδ/αβ T cell lineage commitment regulates γδ T cell developmental programming and shapes the γδ TCR repertoire.
-
- Immunology and Inflammation
Here, we sequenced rearranged TCRβ and TCRα chain sequences in CD4+CD8+ double positive (DP), CD4+CD8- single positive (SP4) and CD4-CD8+ (SP8) thymocyte populations from the foetus and young adult mouse. We found that life-stage had a greater impact on TCRβ and TCRα gene segment usage than cell-type. Foetal repertoires showed bias towards 3’TRAV and 5’TRAJ rearrangements in all populations, whereas adult repertoires used more 5’TRAV gene segments, suggesting that progressive TCRα rearrangements occur less frequently in foetal DP cells. When we synchronised young adult DP thymocyte differentiation by hydrocortisone treatment the new recovering DP thymocyte population showed more foetal-like 3’TRAV and 5’TRAJ gene segment usage. In foetus we identified less influence of MHC-restriction on α-chain and β-chain combinatorial VxJ usage and CDR1xCDR2 (V region) usage in SP compared to adult, indicating weaker impact of MHC-restriction on the foetal TCR repertoire. The foetal TCRβ repertoire was less diverse, less evenly distributed, with fewer non-template insertions, and all foetal populations contained more clonotypic expansions than adult. The differences between the foetal and adult thymus TCR repertoires are consistent with the foetal thymus producing αβT-cells with properties and functions that are distinct from adult T-cells: their repertoire is less governed by MHC-restriction, with preference for particular gene segment usage, less diverse with more clonotypic expansions, and more closely encoded by genomic sequence.