Skip to main content
Download PDF

Senescence-associated secretory phenotypes in mesenchymal cells contribute to cytotoxic immune response in oral lichen planus

Immunity & Ageing volume 20, Article number: 72 (2023) Cite this article

Abstract

Oral lichen planus is a chronic inflammatory condition that adversely affects the oral mucosa; however, its etiology remains elusive. Consequently, therapeutic interventions for oral lichen planus are limited to symptomatic management. This study provides evidence of the accumulation of senescent mesenchymal cells, CD8 + T cells, and natural killer cells in patients with oral lichen planus. We profiled the patients’ tissues using the National Center for Biotechnology Information Gene Expression Omnibus database and found that senescence-related genes were upregulated in these tissues by gene set enrichment analysis. Immunohistochemical analysis showed increased senescent mesenchymal cells in the subepithelial layer of patients with oral lichen planus. Single-cell RNA-seq data retrieved from the Gene Expression Omnibus database of patients with oral lichen planus revealed that mesenchymal cells were marked by the upregulation of senescence-related genes. Cell-cell communication analysis using CellChat showed that senescent mesenchymal cells significantly influenced CD8 + T cells and natural killer cells via CXCL12-CXCR4 signaling, which is known to activate and recruit CD8 + T cells and NK cells. Finally, in vitro assays demonstrated that the secretion of senescence-associated factors from mesenchymal cells stimulated the activation of T cells and natural killer cells and promoted epithelial cell senescence and cytotoxicity. These findings suggest that the accumulation of mesenchymal cells with senescence-associated secretory phenotype may be a key driver of oral lichen planus pathogenesis.

Introduction

Oral lichen planus (OLP) is a persistent inflammatory condition that affects the oral mucosa [1]. It manifests as a painful erythematous center with erosions, ulcerations, and a painless or less painful white variant with papular or reticular patches on the cheeks, lips, tongue, and palate [1]. OLP has a prevalence rate of 0.38–2.05% in the general population [2], with a higher prevalence in middle-aged females [1,2,3]. The etiology of OLP remains poorly understood, although various factors have been postulated, including genetics, dental materials, medications, infectious agents, autoimmune disorders, immunodeficiency, dietary factors, allergies, stress, habits, trauma, diabetes, hypertension, cancer, and gastrointestinal diseases [4, 5]. Therefore, the current treatment for OLP primarily focuses on symptom relief rather than a definitive cure. A better understanding of the underlying pathogenesis of OLP is crucial to develop novel therapeutic targets [6].

Cellular senescence is characterized by irreversible cell cycle arrest, which can be triggered by many factors, including DNA damage, telomere dysfunction, oncogene activation, and organelle stress [7, 8]. The main physiological purpose of cellular senescence is to prevent the proliferation of damaged cells and trigger tissue repair by secreting various soluble molecules, a phenotype termed senescence-associated secretory phenotype (SASP) [7, 8]. The list of these molecules is not comprehensive, and the molecules can vary based on cell type and triggering factors but usually include interleukins, chemokines, growth factors, metalloproteinases, and other insoluble proteins and extracellular matrix components [7,8,9]. However, aging or persistent damage causes the accumulation of senescent cells and impairs cell removal by the immune system, leading to the accumulation of chronic senescent cells and the promotion of chronic inflammatory pathologies [10, 11].

Histopathologically, OLP tissues exhibit a band of T lymphocytes that infiltrate the lamina propria and cytoid bodies in the epithelial layer [3, 12]. Both CD4 + helper T cells and CD8 + cytotoxic T cells are involved in OLP. The activated helper T cells can secrete interleukin (IL)-2 and interferon (IFN)-γ, which activate the cytotoxic T cells and promote their proliferation [4]. Activated cytotoxic T cells can trigger the apoptosis of basal keratinocytes and result in the liquefaction degeneration of basal epithelial cells, which are typically found in OLP lesions and contribute to the chronic inflammation observed in OLP [13]. Natural killer (NK) cells and T cells originate from lymphoid progenitor cells and share developmental proximity. NK cells have cytolytic functions that resemble those of CD8 + cytotoxic T lymphocytes14. The CD56highCD16− NK cell subset found in cutaneous lichen planus (LP) lesions are highly positive for perforin and natural cytotoxic receptors NKG2D and NKp44 [14], suggesting that NK cells are involved in LP pathology [14]. Although cutaneous LP and OLP are distinguished by a marked heterogeneity of clinical course, both share similar histological features. Therefore, both cytotoxic CD8 + T cells and NK cells are speculated to contribute to OLP pathology [12, 15].

Although SASP attracts various immune cells, including NK and CD8 + T cells [16], to sites of chronic inflammation, studies on senescence and the SASP in OLP are limited. Previous studies have reported an increase in the expression levels of NF-kappa B-associated cytokines, such as IL-1 A, IL-6, IL-8, and TNF-alpha, in the saliva and oral fluid of patients with OLP [17, 18]. These inflammatory cytokines are secreted by senescent cells as significant components of the SASP in chronic inflammatory pathology in various organs [19,20,21,22]. In addition, several studies have reported that the number of p21CIP1/WAF1, p16INK4A, and p53-positive cells, which are cell senescence biomarkers, was increased in the epithelial layer of patients with OLP [23,24,25]. Based on these findings, a hypothesis was proposed regarding an association between senescence and OLP.

Here, we show that senescence-related genes are enriched in the tissues of patients with OLP based on gene set enrichment analysis (GSEA) results. Single-cell RNA sequencing (scRNA-seq) and immunohistochemical (IHC) analysis showed that mesenchymal cells upregulate senescence-related genes in patients with OLP, which accumulate in the subepithelial layer that is heavily infiltrated by CD45 + immune cells. Cell-cell communication analysis using scRNA-Seq revealed that senescent mesenchymal cells significantly influenced CD8 + T cells and NK cells via CXCL12-CXCR4 signaling. Lastly, in vitro experiments demonstrated that SASP from senescent fibroblasts enhances the activation and proliferation of NK cells and T cells and induces senescence and cytotoxicity in epithelial cells.

Results

Expression of cellular senescence-related gene sets was enhanced in patients with OLP

To explore whether cellular senescence contributes to OLP, we performed GSEA using the GSE38616 and GSE52130 datasets available in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database. We used several gene sets related to cellular senescence, including the Global Senescence Literature Curated 2020, SASP Literature Curated gene set, SenMayo, Fridman_Up, Purcell, Casella_Up, and Hernandez [26,27,28,29,30,31]. GSEA showed that OLP patient-derived oral mucosa (GSE38616) and oral epithelium (GSE52130) were significantly enriched in the Global Senescence Literature Curated 2020, SASP Literature Curated gene set, SenMayo, Fridman_Up, Purcell, and Casella_Up (Fig. 1a).

Fig. 1
figure 1

Senescence-related gene set enrichment analysis in patients with oral lichen planus and control subjects. Gene set enrichment analysis (GSEA) was performed for six gene sets related to senescence and senescence-associated secretory phenotype, including the Global Senescence Literature Curated 2020, SASP Literature Curated gene set, SenMayo, Fridman_Up, Purcell, Cassela_Up, and Hernandez, in patients with oral lichen planus and control subjects from the Gene Enrichment Omnibus database ((a) GSE38616 and (b) GSE52130). The normalized enrichment scores and false discovery rate q-values are listed in each GSEA plot

Full size image

Senescent mesenchymal cells were increased in patients with OLP

We performed scRNA-seq data analysis using the GSE211630 dataset in the NCBI GEO database to identify the cell types that showed senescence. The datasets included buccal mucosal tissue harvested from one healthy control, three patients with non-erosive OLP (NEOLP), and two patients with erosive OLP (EOLP). Ten distinct cell clusters were annotated using uniform manifold approximation and projection (UMAP) (Fig. 2a, Supplementary Fig. 1). Highly expressed genes are shown in the heatmap of each cell cluster (Fig. 2b). Single-sample GSEA (ssGSEA) was performed to identify senescence-related gene expression. The Global Senescence Literature Curated 2020, SASP Literature Curated gene set, SenMayo, Fridman_Up, Purcell, Casella_Up, and Hernadez senescence-related gene sets were used for ssGSEA [26,27,28,29,30,31]. We found that mesenchymal cell clusters highly expressed senescence-related gene sets in EOLP and NEOLP in the heatmap and violin plot, which showed a mean enrichment score of seven senescence-related gene sets (Fig. 2c, d). To further analyze mesenchymal cell clusters, we performed a sub-cluster analysis. We confirmed that the mesenchymal cell clusters expressed PDGFRA, a marker of mesenchymal cells. UMAP showed eight distinct clusters, and ssGSEA revealed that clusters 0, 1, 2, and 3 in EOLP and NEOLP indicated higher expression of senescence-related gene sets (Supplementary Fig. 2a, b). The cluster with higher expression of senescence-related gene sets was annotated as cluster MC1, and the cluster with lower expression of senescence-related gene sets was annotated as cluster MC2 (Fig. 3b). The heatmap and violin plot of ssGSEA showed that cluster MC1 of EOLP and NEOLP indicated higher expression of senescence-related gene sets, and cluster MC1 of control and cluster MC2 of control, EOLP, and NEOLP indicated lower expression of senescence-related gene sets (Fig. 3c, d). Based on these results, we defined cluster MC1 of EOLP and NEOLP as senescent cell (SEN) clusters and cluster MC1 of control and cluster MC2 of control, EOLP, and NEOLP as non-senescent cell (NSEN) clusters. The SEN cluster also increased the expression of NFKB1 and NFKB2, which are triggers of the SASPs (Fig. 3e) [32]. We also confirmed the presence of senescent mesenchymal cells by IHC analysis, and p21CIP1/WAF1 and PDGFRα double-positive cells were increased in the subepithelial layer of patients with OLP (Fig. 3f, g). p16INK4A, another senescence marker, was also increased in the OLP subepithelial layer (Supplementary Fig. 3a, b).

Fig. 2
figure 2

Mesenchymal cells had increased senescence-related gene expression in oral lichen planus. (a) UMAP plot of buccal mucosa cells from control, EOLP, and NEOLP patients, colored by cell clusters. (b) Heatmap of top differentially expressed genes in each cell cluster). (c) Heatmap of the ssGSEA enrichment scores for seven senescence-related gene sets. Higher ssGSEA enrichment scores are depicted in red, and genes with lower expression are depicted in blue. (d) Violin plot of mean ssGSEA enrichment scores for seven senescence-related gene sets in each cell cluster and disease. UMAP, uniform manifold approximation and projection; EOLP, erosive oral lichen planus; NEOLP, non-erosive oral lichen planus; ssGSEA, single-sample gene set enrichment analysis

Full size image
Fig. 3
figure 3

Identification of senescent mesenchymal cells in oral lichen planus (OLP). (a) UMAP plot of mesenchymal cells from control, EOLP, and NEOLP patients colored according to PDGFRA expression. (b) UMAP plot of mesenchymal cells from control, EOLP, and NEOLP patients., colored by cluster. (c) Heatmap of the ssGSEA enrichment scores for seven senescence-related gene sets. Higher ssGSEA enrichment scores are depicted in red, and genes with lower expression are depicted in blue. (d) Violin plot of mean ssGSEA enrichment scores for seven senescence-related gene sets in each cell cluster and disease. (e) UMAP plot of mesenchymal cells from EOLP and NEOLP patients colored by senescent cluster and feature plot of NFKB1 and NFKB2. Quantification of the gene expression of NFKB1 and NFKB2 in each cluster. (f) Representative images of immunohistochemical staining for p21, PDGFRα, and CD45 in the subepithelial layer of oral mucosa sections from OLP and control groups. (g) Quantification of p21 + PDGFRα + cells in the subepithelial layer. Data are presented as mean and standard error (SE) with dot plots. The circular dot plots represent females, and the triangles represent males. P-values were determined using “FindMarkers” function of Seurat for single-cell RNA-seq data and two-tailed Student’s t-test for immunohistochemical staining

UMAP, uniform manifold approximation and projection; EOLP, erosive oral lichen planus; NEOLP, non-erosive oral lichen planus; ssGSEA, single-sample gene set enrichment analysis

Full size image

Ligand-receptor analysis of senescent and non-senescent mesenchymal cells in patients with OLP

Using newly defined cell clusters, including SEN_Mesenchymal, NSEN_Mesenchymal, Monocyte, SMC, CD8, CD4, DC, Endothelial, NK, Treg, Epithelial, Neutrophil, Plasma, and B cell clusters, we performed ligand-receptor analysis using CellChat [33]. Ligand-receptor analysis showed that the SEN_Mesenchymal cluster had higher cell-cell communication (Fig. 4a). In-depth analysis showed the top 20 ligand-receptor signaling pathways and the strength of outgoing/incoming signaling in each cell cluster (Fig. 4b). The data showed that the SEN_Mesenchymal cluster was the most expressed outgoing signal, and the CD8 cluster was the most expressed incoming signal (Fig. 4b). NK and Neutrophil clusters also showed high expression of incoming signals (Fig. 4b). We next investigated the dominant signals in the SEN_Mesenchymal cluster instead of the NSEN_Mesenchymal cluster and found that COLLAGEN, LAMININ, and CXCL signals were dominant (Fig. 4b, c). The bubble plot showed significantly increased cell-cell communication in NSEN_Mesenchymal and SEN_Mesenchymal cell clusters, and the data revealed that COL1A1-CD44, COL4A1-CD44, COL4A2-CD44, COL6A1-CD44, COL6A2-CD44, LAMB1-CD44, and CXCL12-CXCR4 were significantly higher in the communication between SEN_Mesenchymal cells, epithelial cells, and immune cells, including the CD4, CD8, NK, Monocyte, Treg, Neutrophil, Plasma, and B cell clusters (Fig. 4d). These results indicate that senescent mesenchymal cells may strongly affect other cell populations.

Fig. 4
figure 4

Senescent mesenchymal cells have higher cell − cell communications with immune cells. (a) Circles plot significant communication in each cell cluster. The thicker line indicates strong communication. (b) Heatmap of the top 20 outgoing/incoming signaling patterns in each cell cluster. (c) Chord diagram of the significant network of COLLAGEN, LAMININ, and CXCL signaling pathways in each cell cluster. (d) Bubble plot of ligand-receptor pairs in mesenchymal cell clusters. The colors in the bubble plot are proportional to the communication probability, where blue and red correspond to the smallest and largest values, respectively

Full size image

Immune cells were increased in the senescent mesenchymal cells accumulated region

Previous studies have shown that accumulated senescent cells contribute to chronic inflammation via excessive immune cell infiltration through SASP expression [34]. Excessive immune cell infiltration into the subepithelial layer is also a histopathological characteristic of OLP [3, 12]. We performed an IHC analysis to investigate whether senescent mesenchymal cells increased in regions with high immune cell infiltration. We found that the extent of highly CD45 + immune cell–infiltrated regions increased in p21-positive mesenchymal cells compared to that in both of Ctrl tissue and weak to moderate immune cell–infiltrated regions (Fig. 5a, b).

Fig. 5
figure 5

Senescent mesenchymal cell count increased in the strongly immune cell–infiltrated regions. (a) Representative images of immunohistochemical staining for p21, PDGFRα, and CD45 in the subepithelial layer of oral mucosa sections from patients with oral lichen planus (OLP) and control groups. CD45-positive immune cell infiltration was graded as weak-moderate (≤ 40% area of CD45-positivity) and strong (> 40% area of CD45-positivity). (b) Quantification of p21 positive cells in control and patients with OLP. Data are presented as the mean and standard error (SE) with dot plots. The circular dot plots represent females, and the triangles represent males. P-values were determined using Tukey’s method for one-way analysis of variance (**p < 0.001)

Full size image

Senescent mesenchymal cell-derived SASP induces epithelial cell senescence and NK and T cell proliferation and activation in vitro

To investigate whether senescent mesenchymal cells contribute to OLP pathology, we performed an in vitro assay using senescent mesenchymal-derived SASP-containing conditioned media. First, we confirmed that TIG-118 expressed PDGFRα, a mesenchymal cell marker (Supplementary Fig. 4). Doxorubicin (DOX) was used to induce senescence. The DOXO-treated TIG-118 increased the percentage of senescence-associated beta-galactosidase (SA-β-Gal)-positive cells and the expression levels of CDKN1A, CDKN2A, and NFKB2 (Fig. 6a-c), which triggers CXCL12 [35]. Next, to investigate the effect of the SASP on epithelial and immune cells, we collected conditioned media from control TIG-118 and DOX-induced senescent TIG-118 cells. HaCaT, a human skin-derived epithelial cell line; Jurkat, a human T cell line; and KHYG-1, a human NK cell line, were exposed to the conditioned media (Fig. 6d). From scRNA-seq data analysis and GSEA analysis for KEGG NATURAL KILLER CELL-MEDIATED CYTOTOXICITY and GOBP NATURAL KILLER CELL ACTIVATION using GSE38616 and GSE52130 datasets (Supplementary Fig. 5), we focused on T cells and NK cells. Senescent TIG-118-derived conditioned media (SEN-CM)-exposed HaCaT cells showed a significant increase in lactate dehydrogenase (LDH) release and expression of CDKN2A, IL6, and SERPINE1, which are senescence and SASP markers, compared to control TIG-118-derived conditioned media (Ctrl-CM) (Fig. 6e). SEN-CM-exposed Jurkat cells showed increased cell proliferation by assessing WST-8 assay and MKI67 mRNA expression and increased the mRNA expression of CD25, whose expression level was increased in proliferating Jurkat cells (Fig. 6f) [36]. In addition, SEN-CM-exposed KHYG-1 cells also showed increased cell proliferation and mRNA expression of PRF1 and IFNA13, which encode perforin and type I IFN associated with NK cell-mediated cytotoxicity (Fig. 6g). These results suggest that senescent mesenchymal cells have the potential to contribute to OLP pathology by inducing epithelial cell damage and senescence, as well as T cell and NK cell proliferation and activation.

Fig. 6
figure 6

Senescent mesenchymal cell-derived SASP induces epithelial cell senescence and NK and T cell proliferation and activation in vitro. (a) Representative images of SA-β-Gal expression in TIG-118 cells with or without doxorubicin (DOXO)-treatment. (b) Quantification of SA-β-Gal-positive cells. (c) Relative mRNA expression of CDKN1A, CDKN2A, and NFKB2 in control and DOXO-induced senescent cells. (d) Protocols for in vitro study. (e) Quantification of LDH cytotoxicity and relative mRNA expression of CDKN2A, IL6, and SERPINE1 in HaCaT cells treated with control cell-derived conditioned media (Ctrl-CM) and senescent cell-derived conditioned media (SEN-CM). (f) Quantification of WST-8 proliferation and relative mRNA expression of MKI67, CD25, and CD69 in Ctrl-CM- and SEN-CM-treated Jurkat cells. (g) Quantification of WST-8 proliferation and relative mRNA expression of PRF1, GPR12, IFNA10, IFNA13, and IFNG in Ctrl-CM or SEN-CM-treated KHYG-1 cells. Data are presented as mean and standard error (SE) with dot plots. P-values were determined using a two-tailed Student’s t-test. (*P < 0.05 and **P < 0.01)

NK, natural killer; SASP, senescence-associated secretory phenotype

Full size image

Discussion

This study revealed that senescent mesenchymal cells are associated with OLP pathology; they activate immune cells, including T cells and NK cells, and induce senescence and cytotoxicity in epithelial cells.

Several studies have reported increased p21CIP1/WAF1, p16INK4A, and p53-positive cells in patients with OLP by IHC analysis [23,24,25]. Other studies showed that SMP-30, a marker of senescence and anti-apoptosis, was upregulated in OLP [37, 38]. Although these findings indicate a relationship between senescence and OLP pathogenesis, using a single senescence marker is insufficient to detect senescent cells, and applying a validated senescence gene set is considered a viable method for detecting senescence in vivo [27]. Our study analyzed seven validated senescence gene sets and two transcriptome data cohorts of patients with OLP from GEO to explore the association between senescence and OLP. The results of GSEA indicated significant enrichment of senescence gene sets in patients with OLP, confirming the possibility of a connection between senescence and OLP.

scRNA-seq and IHC analyses of patients with OLP revealed that senescence is exhibited by mesenchymal cells located in the subepithelial layer36. We utilized PDGFRα as a marker for mesenchymal cells or fibroblasts. There is significant documentation of the features of senescent PDGFRα-positive cells in various organs, such as the skin, skeletal muscle, lung, and adipose tissue [39,40,41,42,43]. These studies reported that senescent PDGFRα-positive cells exhibited increased SASP expression, which can promote the recruitment of immune cells. In the scRNA-seq analysis, we found that the mesenchymal cell cluster, which expressed PDGFRA, increased senescence-related gene expression in EOLP and NEOLP but not in the control tissue. Using CellChat, a ligand-receptor analysis tool, we found that senescent mesenchymal cell clusters strongly interact with immune cells, especially CD8 + T cells and NK cells. CellChat analysis also showed that senescent mesenchymal cell clusters have significant connections with immune cells via CXCL12-CXCR4. CXCL12 serves as ligands for CXCR4, a receptor that plays an important role in activating and homing CD8 + T cells and NK cells towards CXCL12-expressing cells [44,45,46]. Our findings suggest that augmented CXCL12-CXCR4 signaling in senescent mesenchymal cells may cause the accumulation of CD8 + T cells and NK cells in the vicinity of senescent mesenchymal cells in patients with OLP.

In the present study, the scRNA-seq analysis revealed that senescent mesenchymal cells are connected to immune cells via collagen IV, VI-CD44, and laminin-CD44. CD44 is a cellular adhesion molecule expressed in many different cell types and a receptor for hyaluronic acid, laminin, collagen, or fibronectin [47, 48]. CD44 promotes the migration of many cell types, including T cells and NK cells, by regulating their adhesion and extravasation [49]. CD44-positive cell-matrix interactions promote cell proliferation, adhesion, migration, and lymphocyte activation [50]. Although the functional relationship between senescent mesenchymal cells generating collagen, laminin, and CD44 remains unclear, the cell-matrix interaction with CD44 may be related to the OLP pathomechanism.

In our in vitro experiments, we explored the effects of senescent cells on immune and epithelial cells using a SASP-containing conditioned media (SEN-CM) isolated from senescent fibroblasts. Our findings indicate that SEN-CM stimulates the proliferation of Jurkat cells and augments the expression of the cellular proliferation markers MKI67 and CD25 [36]. SEN-CM also activates KHYG-1 cells, which are cytotoxic to NK cells [51]. SEN-CM treatment increased the expression of PRF1, which encodes perforin, and GPR12, which regulates the cytolytic activity of NK cells [52]. Moreover, SEN-CM increased the cytotoxicity of HaCaT cells, an epithelial cell line, and upregulated the expression of CDKN2A, IL6, and SERPINE1, which are typical markers of senescence. These findings suggest that senescent mesenchymal cell-derived SASP may contribute to OLP pathology by activating T cells and NK cells and inducing epithelial damage and senescence.

Although OLP is a T cell-mediated inflammatory disease of the oral mucosa, this study identified an NK cell population that may be associated with OLP pathology based on scRNA-seq analysis and GSEA. GSEA revealed significant enrichment of KEGG NATURAL KILLER CELL-MEDIATED CYTOTOXICITY and GOBP NATURAL KILLER CELL ACTIVATION in patients with OLP. The role of NK cells in OLP pathogenesis remains unclear; however, previous studies have reported that NK cells constituted 2.3–15% of cellular infiltrates in cutaneous lichen planus [14, 53]. Therefore, NK cells may provide an early stimulatory signal for the mobilization of CD4 + and CD8 + T cells upon their arrival in lichen planus-inflamed areas [54]. Our data from the bioinformatics analysis suggest the possibility of NK cell involvement in OLP pathogenesis, which may support the hypothesis mentioned above.

Senescent cells are implicated in the pathogenesis of periodontal diseases and oral submucous fibrosis [55, 56]. In periodontal disease, persistent infection with Gram-negative bacteria leads to DNA damage and subsequent senescence of host cells, resulting in chronic inflammation [55]. In oral submucous fibrosis, senescence occurs in keratinocytes, fibroblasts, and endothelial cells induced by exposure to areca nut alkaloids, causing DNA damage and increased ROS levels [56, 57]. While senescence is a known suppressor of oral squamous cell carcinoma (OSCC) development owing to permanent cell cycle arrest, senescence of the non-epithelial components, especially cancer-associated fibroblasts (CAFs), contributes to OSCC development [58]. CAFs are critical components of the tumor microenvironment with diverse functions, including matrix deposition and remodeling, and extensive reciprocal signaling between cancer cells to promote cancer progression [59]. CAFs originate from normal fibroblasts, pericytes, smooth muscle cells, or mesenchymal cells [60]. Further investigation is necessary to determine whether senescent mesenchymal cells in OLP contribute to tumor suppression or development, which may help in the development of senescent cell-targeting therapies for OLP.

There are several limitations in this study. First, human skin fibroblast (TIG-118) was used in our in vitro experiments. Fibroblasts are thought to play an important role in wound healing and scar formation and oral mucosal fibroblasts proliferate slightly more than skin fibroblasts [61]. We also used HaCaT cells which are a spontaneously immortalized keratinocyte cell line derived from adult human skin. HaCaT cells were previously used in oral cell biology as an oral keratinocyte model [62]. However, keratinocytes from the skin and oral mucosa have differences in morphology, gene expression, and proliferative capacity [63]. These differences may have affected the results of our in vitro experiments. Second, this study showed more CD45 + cells with senescent mesenchymal cells and CD8+/NK proliferation and secretory function. Although cell–cell communication analysis using scRNA-Seq revealed that senescent mesenchymal cells significantly influenced CD8 + T and NK cells, further experiments are needed to elucidate the details of the induction and migration of CD8 + T and NK cells by senescent mesenchymal cells in OLP pathology. Third, we used healthy regions from patients with leukoplakia as controls, acquired as part of the usual diagnostic procedure. Previous research used healthy tissue and healthy regions from patients with leukoplakia as controls [64]. The study showed that healthy tissue and healthy regions from patients with leukoplakia have significant differences from leukoplakia and OSCC tissues, as determined using histological and mechanical evaluations by measuring tissue stiffness and the expression levels of p16 and BCL-2, which showed that healthy regions from patients with leukoplakia could be used as a control [64]. Although healthy regions from patients with leukoplakia can be used as a control, healthy tissue could also be analyzed as a control.

In conclusion, this study indicates that senescent mesenchymal cells are associated with OLP pathology by activating immune cells, including T cells and NK cells, and inducing senescence and cytotoxicity in epithelial cells (Fig. 7). These findings reveal a new pathological mechanism of OLP, which is mediated by senescence, and provide potential therapies targeting senescent cells.

Fig. 7
figure 7

Overview of key findings

Full size image

Materials and methods

GEO data

Gene expression data for the OLP and control oral tissue samples were obtained from the NCBI GEO using the GEOquery R package [65]. The accession numbers of the datasets are GSE38616 and GSE52130. Both datasets included seven OLP and seven control samples. scRNA-seq dataset GSE211630 was obtained from the NCBI GEO. The dataset included the buccal mucosa from two EOLP patients, three NEOLP patients, and one healthy control tissue [15].

GSEA

GSEA was performed using the GSEA software (https://www.gsea-msigdb.org/gsea/index.jsp), a joint project of the UC San Diego and Broad Institute [66, 67]. GSEA analyses were performed using seven gene sets related to senescence and SASP, including the Global Senescence Literature Curated 2020, SASP Literature Curated gene set, SenMayo, Fridman Senescence Up, Purcell, Casella Up, and Hernandez gene sets [26,27,28,29,30,31], and two gene sets related to NK cells, including KEGG NATURAL KILLER CELL-MEDIATED CYTOTOXICITY and GOBP NATURAL KILLER CELL ACTIVATION.

ScRNA-seq analysis

The scRNA-seq dataset GSE211630 was imported into R (version 4.2.2), and Seurat (version 4.3.0) was used [68, 69]. The data were filtered to include genes detected in > 3 cells, cells with 100–5,000 detected genes, and a mitochondrial gene expression ratio of < 5%. The scRNA-seq datasets were normalized using the “NormalizeData” function, and 2000 highly variable genes were identified using the “FindVariableFeatures” function. Then, the canonical correction analysis method with the “FindIntegrationAnchors” and “IntegrateData” functions was used to remove batch effects. Next, we performed data scaling, used the PCA method to reduce the dimensions, and used the UMAP to visualize the data. Cell clusters were identified using the “FindNeighbors” function and the “FindClusters” function with a resolution of 0.5. We utilized already known marker genes to annotate cell clusters: mesenchymal cell (PDGFRA), smooth muscle cell (SMC) (ACTA2), endothelial cell (VWF), epithelial cell (KRT5), neuron (NRXN1), skeletal muscle (DES), immune cells (PTPRC), T cell (CD3E), CD4 + T cell (CD4), regulatory T cell (Treg) (FOXP3), CD8 + T cell (CD8A), NK cell (GNLY), monocyte (AIF1), Dendritic cell (DC) (FCER1A), plasma cell (MZB1), neutrophil (S100A8), and B cell (MSA4A1).

ssGSEA

ssGSEA was performed using the Easy single-cell analysis platform for enrichment (version 1.8.0) [70]. Enrichment scores were calculated for each of the seven senescence-related gene sets, including the Global Senescence Literature Curated 2020, SASP Literature Curated gene set, SenMayo, Fridman_Up, Purcell, Casella_Up, and Hernadez senescence-related gene sets [26,27,28,29,30,31], and visualized in a heatmap and violin plot using dittoSeq (version 1.10.0) [71].

Cell-cell communication analysis based on ligand-receptor interaction

CellChat, which quantitatively infers and analyzes intercellular communication networks from scRNA-seq data, was used for cell-cell communication analysis based on ligand-receptor interaction (version 1.6.1) [33]. CellChat assigns a communication probability value based on the law of mass action and the average expression values of a ligand by one cell cluster and that of a receptor by another cell cluster. The statistical significance of communication probability values was assessed using a permutation test. P < 0.05 was considered significant.

Patients and tissue samples

The study protocol was approved by the Institutional Review Board of Sapporo Medical University (approval number: 332 − 153).

Twenty patients diagnosed with OLP at the Department of Oral Surgery, Sapporo Medical University, were enrolled in this study, including 10 males and 10 females aged 41–89. Healthy regions from patients with leukoplakia, including five males and five females aged 44–71 years, were used as controls. Hematoxylin and eosin staining images are shown in Supplementary Fig. 6.

Immunohistochemistry

Paraffin-embedded tissues were sectioned (3-µm thickness). Sections were deparaffinized and rehydrated for immunostaining. Antigen retrieval was performed in a microwave oven (95–98°C for 10 min) using a citrate buffer (10 mM sodium citrate, pH 6.0). After cooling, the slides were washed twice with deionized water and quenched with 3% hydrogen peroxide for 10 min. Then sections were blocked with 1% bovine serum albumin in Tris-buffered saline containing Tween-20 (TBST) for 15 min at room temperature (RT: 20–25°C) and then incubated with primary antibodies overnight at 4°C or for 1 h at RT. After washing three times with TBST for 5 min each, the sections were incubated with SignalStain Boost IHC Detection Reagent (Cell Signaling Technology, MA, USA) for 30 min at RT in the dark. The sections were then washed in TBST three times for 5 min each and treated with TSA Plus Working Solution (Fluorescein, Cyanine 3, and Cyanine 5; Akoya Biosciences, MA, USA) for 10 min at RT in the dark. For multiplex staining, stripping was performed in a microwave oven (95–98°C for 10 min) using citrate buffer. After cooling, the cells were stained with different tyramide fluorescent labels, as described above. The following antibodies were used: anti-CD45 (ab8216, mouse monoclonal, 1:1,000; Abcam, Cambridge, UK), anti-p21 (ZBR1141, rabbit monoclonal, 1:1,000; Sigma-Aldrich, MO, USA), anti–PDGFRα (#3174, 1:1,000; Cell Signaling Technology), and anti–p16INK4A (ZBR1437, rabbit monoclonal, 1:200; Sigma-Aldrich). Nuclei were stained with 4’,6-diamidino-2-phenylindole (Dojindo Laboratories, Kumamoto, Japan). The sections were observed under a fluorescence microscope (Axio Observer 7; ZEISS, Baden-Wurttemberg, Germany). Images of randomly selected epithelial and subepithelial regions were acquired, and the average value of three fields per patient was used for statistical analysis. The fluorescence intensities of p21 and p16INK4A were analyzed using the image analysis function of ZEN (ZEISS). CD45-positive immune cell infiltration was graded as weak-moderate (area of CD45-positivity is 40% or less) and strong (area of CD45-positivity is more than 40%) [34].

Cell culture

TIG-118, a human skin fibroblast (JCRB cell bank: JCRB0535); HaCaT, a human skin epithelial cell line (Sapporo Medical University Cell Bank); Jurkat, a human T cell line (RIKEN cell bank: RCB3052); and KHYG-1, a human NK cell line (JCRB cell bank: JCRB0156) were used for the in vitro study. TIG-118 and HaCaT cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) containing 10% fetal bovine serum (FBS). Jurkat and KHYG-1 cells were cultured in RPMI1640 medium containing 10% FBS, and IL-2 was added to the KHYG-1 cells to promote cell proliferation.

Western blotting

Cell pellets were lysed in a cell lysis buffer (Nacalai Tesque Inc., Kyoto, Japan) with a protease and phosphatase inhibitor cocktail (Nacalai). All lysates were centrifuged, and proteins were quantified using the BCA protein assay (Thermo Fisher Scientific, MA, USA). Then, 8 µg of protein from each sample was electrophoresed on 4–12% Bris-Tris protein gels and transferred to low-fluorescence PVDF membranes (Bio-Rad Laboratories, CA, USA). Blots were probed with the following primary antibodies: anti-PDGFRα (1:1000, Cell Signaling Technology), anti-β actin (1:1000, MEDICAL&BIOLOGICAL LABORATORIES (MBL), Tokyo, Japan), and anti-α tubulin (1:1000, MBL). The blots were probed with anti-rabbit IgG (IRDye® 800CW, 1:20000, LI-COR Biosciences, NE, USA) and anti-mouse IgG (IRDye® 680RD, 1:20000, LI-COR) secondary antibodies. Fluorescence-based detection was performed using LI-COR Odyssey XF (LI-COR).

SA-β-Gal staining

To detect cellular senescence, we performed SA-β-Gal staining using the Senescence b-Galactosidase Staining Kit (Cell Signaling Technology). Cells were observed using an inverted microscope (Primovert, ZEISS), and the percentage of SA-β-Gal-positive cells was calculated by dividing the number of SA-β-Gal-positive cells by the total number of cells observed.

Senescence induction and collection of conditioned media

TIG-118 cells were seeded in culture flasks and cultured overnight in DMEM supplemented with 10% FBS. The next day, the complete medium was replaced with 250 nM DOX-containing complete medium and incubated for 24 h to induce senescence. After DOX treatment, TIG-118 cells were washed with phosphate-buffered saline and cultured in a fresh complete medium. Eight days after DOX treatment, the complete medium was replaced with a low-serum medium containing 0.2% FBS and incubated for 48 h. After 48 h, the medium was replaced with a serum-free and phenol-red-free medium, and the cells were incubated for 24 h. After 24 h, the SASP-containing conditioned media was collected. The collected conditioned media was centrifuged, and the supernatants were stored in a -80 °C refrigerator for subsequent experiments.

Cell proliferation and cytotoxicity lactate dehydrogenase assay

HaCaT, Jurkat, and KHYG-1 cells were cultured in conditioned media supplemented with 2% FBS for 48 h. After 48 h, cell proliferation or cytotoxicity lactate dehydrogenase assay was performed using the Cell Counting Kit-8 (Dojindo) or Cytotoxicity LDH Assay kit (Dojindo).

RNA extraction and quantitative real-time PCR

Total RNA was isolated from HaCaT, Jurkat, and KHYG-1 cells cultured in a conditioned media containing 2% FBS for 48 h using TRI Reagent (Molecular Research Center Inc., OH, USA). The isolated RNA was reverse-transcribed into cDNA using the iScript Advanced cDNA Synthesis Kit (Bio-Rad). Quantitative PCR was performed using the SsoAdvanced Universal SYBR Green Supermix (Bio-Rad) in a CFX Connect Real-Time PCR Detection System (Bio-Rad) under the following cycling conditions: 95 °C for 30 s, followed by 40 cycles of amplification (95 °C for 10 s and 60 °C for 30 s). The primer sequences used for the PCR are listed in Supplementary Table 1. The samples were compared using the ΔΔCt method.

Statistical analysis

Quantitative data are reported as means and standard errors (SE). Each data point was plotted using dot and violin plots with ggplot2, a plotting system for R based on the Grammar of Graphics. Normality was assessed using the Shapiro–Wilk test. Pairwise t-tests were used to compare two groups. One-way analysis of variance was conducted to assess differences between three or more groups. P-values for multiple comparisons were adjusted using Tukey’s method. Statistical analyses were performed using EZR, a graphical user interface for R [72]. Two-sided P-values less than 0.05 were considered statistically significant. For the statistical analysis of GSEA, the Singal2Noise method, a ranking method based on t-statistics, was used to calculate the P-value, and the false discovery rate method was used to calculate the q-value.

Data Availability

Data and materials not indicated in this manuscript are available from the corresponding author.

References

  1. Alrashdan MS, Cirillo N, McCullough M. Oral lichen planus: a literature review and update. Arch Dermatol Res. 2016;308:539–51.

    Article  PubMed  Google Scholar 

  2. Li C, Tang X, Zheng X, Ge S, Wen H, Lin X, et al. Global prevalence and incidence estimates of oral Lichen Planus: a systematic review and Meta-analysis. JAMA Dermatol. 2020;156:172–81.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Olson MA, Rogers RS, Bruce AJ. Oral lichen planus. Clin Dermatol. 2016;34:495–504.

    Article  PubMed  Google Scholar 

  4. Roopashree MR, Gondhalekar RV, Shashikanth MC, George J, Thippeswamy SH, Shukla A. Pathogenesis of oral lichen planus–a review. J Oral Pathol Med. 2010;39:729–34.

    Article  PubMed  CAS  Google Scholar 

  5. Aly EAH, Burgess P. Use of laser in the relief of malignant dysphagia: a district hospital experience. Dig Surg. 2002;19:3–8.

    Article  PubMed  Google Scholar 

  6. Cheng Y-SL, Gould A, Kurago Z, Fantasia J, Muller S. Diagnosis of oral lichen planus: a position paper of the American Academy of Oral and Maxillofacial Pathology. Oral Surg Oral Med Oral Pathol Oral Radiol. 2016;122:332–54.

    Article  PubMed  Google Scholar 

  7. Herranz N, Gil J. Mechanisms and functions of cellular senescence. J Clin Invest. 2018;128:1238–46.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Gorgoulis V, Adams PD, Alimonti A, Bennett DC, Bischof O, Bishop C, et al. Cellular Senescence: defining a path Forward. Cell. 2019;179:813–27.

    Article  PubMed  CAS  Google Scholar 

  9. Ferrucci L, Fabbri E. Inflammageing: chronic inflammation in ageing, Cardiovascular Disease, and frailty. Nat Rev Cardiol. 2018;15:505–22.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Childs BG, Gluscevic M, Baker DJ, Laberge R-M, Marquess D, Dananberg J, et al. Senescent cells: an emerging target for Diseases of ageing. Nat Rev Drug Discov. 2017;16:718–35.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. van Deursen JM. The role of senescent cells in ageing. Nature. 2014;509:439–46.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Dissemond J. Oral lichen planus: an overview. J Dermatolog Treat. 2004;15:136–40.

    Article  PubMed  Google Scholar 

  13. Xue N, Wang Y, Cheng H, Liang H, Fan X, Zuo F, et al. Regulatory T cell therapy suppresses inflammation of oral mucosa. Front Immunol. 2022;13:1009742.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Carbone T, Nasorri F, Pennino D, Donnarumma M, Garcovich S, Eyerich K, et al. CD56 highCD16 - NK cell involvement in cutaneous lichen planus. Eur J Dermatol. 2010;20:724–30.

    PubMed  CAS  Google Scholar 

  15. Qing M, Shang Q, Yang D, Peng J, Deng J, Jiang L et al. CD8 + tissue-resident memory T cells triggered the erosion of oral lichen planus by the cytokine network. medRxiv. 2022;:2022.10.18.22281149.

  16. Kale A, Sharma A, Stolzing A, Desprez P-Y, Campisi J. Role of immune cells in the removal of deleterious senescent cells. Immun Ageing. 2020;17:16.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Rhodus NL, Cheng B, Myers S, Bowles W, Ho V, Ondrey F. A comparison of the pro-inflammatory, NF-κB-dependent cytokines: TNF-alpha, IL-1-alpha, IL-6, and IL-8 in different oral fluids from oral lichen planus patients. Clin Immunol. 2005;114:278–83.

    Article  PubMed  CAS  Google Scholar 

  18. Rhodus NL, Cheng B, Myers S, Miller L, Ho V, Ondrey F. The feasibility of monitoring NF-κB associated cytokines: TNF-α, IL-1α, IL-6, and IL-8 in whole saliva for the malignant transformation of oral lichen planus. Mol Carcinog. 2005;44:77–82.

    Article  PubMed  CAS  Google Scholar 

  19. Wang B, Varela-Eirin M, Brandenburg SM, Hernandez-Segura A, van Vliet T, Jongbloed EM, et al. Pharmacological CDK4/6 inhibition reveals a p53-dependent senescent state with restricted toxicity. EMBO J. 2022;41:e108946.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Zhang Y, Lin M, Zhang S, Wang Z, Jiang L, Shen J, et al. NF-kappaB-dependent cytokines in saliva and serum from patients with oral lichen planus: a study in an ethnic Chinese population. Cytokine. 2008;41:144–9.

    Article  PubMed  CAS  Google Scholar 

  21. Rhodus NL, Cheng B, Myers S, Bowles W, Ho V, Ondrey F. A comparison of the pro-inflammatory, NF-kappaB-dependent cytokines: TNF-alpha, IL-1-alpha, IL-6, and IL-8 in different oral fluids from oral lichen planus patients. Clin Immunol. 2005;114:278–83.

    Article  PubMed  CAS  Google Scholar 

  22. Rhodus NL, Cheng B, Myers S, Miller L, Ho V, Ondrey F. The feasibility of monitoring NF-kappaB associated cytokines: TNF-alpha, IL-1alpha, IL-6, and IL-8 in whole saliva for the malignant transformation of oral lichen planus. Mol Carcinog. 2005;44:77–82.

    Article  PubMed  CAS  Google Scholar 

  23. Montebugnoli L, Venturi M, Gissi DB, Leonardi E, Farnedi A, Foschini MP. Immunohistochemical expression of p16(INK4A) protein in oral lichen planus. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2011;112:222–7.

    Article  PubMed  CAS  Google Scholar 

  24. Poomsawat S, Buajeeb W, Khovidhunkit S-OP, Punyasingh J. Overexpression of cdk4 and p16 in oral lichen planus supports the concept of premalignancy. J Oral Pathol Med. 2011;40:294–9.

    Article  PubMed  Google Scholar 

  25. González-Moles MA, Bascones-Ilundain C, Gil Montoya JA, Ruiz-Avila I, Delgado-Rodríguez M, Bascones-Martínez A. Cell cycle regulating mechanisms in oral lichen planus: molecular bases in epithelium predisposed to malignant transformation. Arch Oral Biol. 2006;51:1093–103.

    Article  PubMed  Google Scholar 

  26. Binet F, Cagnone G, Crespo-Garcia S, Hata M, Neault M, Dejda A et al. Neutrophil extracellular traps target senescent vasculature for tissue remodeling in retinopathy. Science. 2020;369.

  27. Saul D, Kosinsky RL, Atkinson EJ, Doolittle ML, Zhang X, LeBrasseur NK, et al. A new gene set identifies senescent cells and predicts senescence-associated pathways across tissues. Nat Commun. 2022;13:4827.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Fridman AL, Tainsky MA. Critical pathways in cellular senescence and immortalization revealed by gene expression profiling. Oncogene. 2008;27:5975–87.

    Article  PubMed  CAS  Google Scholar 

  29. Purcell M, Kruger A, Tainsky MA. Gene expression profiling of replicative and induced senescence. Cell Cycle. 2014;13:3927–37.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Casella G, Munk R, Kim KM, Piao Y, De S, Abdelmohsen K, et al. Transcriptome signature of cellular senescence. Nucleic Acids Res. 2019;47:7294–305.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Hernandez-Segura A, de Jong TV, Melov S, Guryev V, Campisi J, Demaria M. Unmasking transcriptional heterogeneity in senescent cells. Curr Biol. 2017;27:2652–2660e4.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Salminen A, Kauppinen A, Kaarniranta K. Emerging role of NF-κB signaling in the induction of senescence-associated secretory phenotype (SASP). Cell Signal. 2012;24:835–45.

    Article  PubMed  CAS  Google Scholar 

  33. Jin S, Guerrero-Juarez CF, Zhang L, Chang I, Ramos R, Kuan C-H, et al. Inference and analysis of cell-cell communication using CellChat. Nat Commun. 2021;12:1088.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Choi YW, Kim YH, Oh SY, Suh KW, Kim Y-S, Lee G-Y, et al. Senescent Tumor cells build a Cytokine Shield in Colorectal Cancer. Adv Sci (Weinh). 2021;8:2002497.

    Article  PubMed  CAS  Google Scholar 

  35. Mukherjee T, Chatterjee B, Dhar A, Bais SS, Chawla M, Roy P, et al. A TNF-p100 pathway subverts noncanonical NF-κB signaling in inflamed secondary lymphoid organs. EMBO J. 2017;36:3501–16.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Shatrova AN, Mityushova EV, Aksenov NA, Marakhova II. CD25 expression on the surface of jurkat cells. Cell Tissue biol. 2015;9:364–70.

    Article  Google Scholar 

  37. Peter CD, Shashidara R, Jain V, Haragannavar VC, Samuel P, Nayak SR. Senescence in oral lichen planus as assessed by the immunohistochemical evaluation of senescence marker protein-30 (regucalcin). Indian J Pathol Microbiol. 2023;66:9–13.

    Article  PubMed  Google Scholar 

  38. Fujita T, Inoue H, Kitamura T, Sato N, Shimosawa T, Maruyama N. Senescence marker protein-30 (SMP30) rescues cell death by enhancing plasma membrane ca(2+)-pumping activity in Hep G2 cells. Biochem Biophys Res Commun. 1998;250:374–80.

    Article  PubMed  CAS  Google Scholar 

  39. Kita A, Saito Y, Miura N, Miyajima M, Yamamoto S, Sato T, et al. Altered regulation of mesenchymal cell senescence in adipose tissue promotes pathological changes associated with diabetic wound healing. Commun Biol. 2022;5:310.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Salzer MC, Lafzi A, Berenguer-Llergo A, Youssif C, Castellanos A, Solanas G, et al. Identity noise and adipogenic traits characterize dermal fibroblast aging. Cell. 2018;175:1575–1590e22.

    Article  PubMed  CAS  Google Scholar 

  41. Saito Y, Chikenji TS, Matsumura T, Nakano M, Fujimiya M. Exercise enhances skeletal muscle regeneration by promoting senescence in fibro-adipogenic progenitors. Nat Commun. 2020;11.

  42. Schafer MJ, White TA, Iijima K, Haak AJ, Ligresti G, Atkinson EJ, et al. Cellular senescence mediates fibrotic pulmonary Disease. Nat Commun. 2017;8:14532.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Reyes NS, Krasilnikov M, Allen NC, Lee JY, Hyams B, Zhou M, et al. Sentinel p16INK4a + cells in the basement membrane form a reparative niche in the lung. Science. 2022;378:192–201.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Goedhart M, Gessel S, van der Voort R, Slot E, Lucas B, Gielen E, et al. CXCR4, but not CXCR3, drives CD8 + T-cell entry into and migration through the murine bone marrow. Eur J Immunol. 2019;49:576–89.

    Article  PubMed  CAS  Google Scholar 

  45. Levy E, Reger R, Segerberg F, Lambert M, Leijonhufvud C, Baumer Y et al. Enhanced Bone Marrow Homing of Natural Killer Cells Following mRNA Transfection With Gain-of-Function Variant CXCR4R334X. Front Immunol. 2019;10 JUN:1262.

  46. Kohli K, Pillarisetty VG, Kim TS. Key chemokines direct migration of immune cells in solid tumors. Cancer Gene Ther. 2022;29:10–21.

    Article  PubMed  CAS  Google Scholar 

  47. Ishii S, Ford R, Thomas P, Nachman A, Steele G, Jessup JM. CD44 participates in the adhesion of human colorectal carcinoma cells to laminin and type IV collagen. Surg Oncol. 1993;2:255–64.

    Article  PubMed  CAS  Google Scholar 

  48. Jalkanen S, Jalkanen M. Lymphocyte CD44 binds the COOH-terminal heparin-binding domain of fibronectin. J Cell Biol. 1992;116:817–25.

    Article  PubMed  CAS  Google Scholar 

  49. Baaten BJG, Tinoco R, Chen AT, Bradley LM. Regulation of Antigen-experienced T cells: lessons from the quintessential memory marker CD44. Front Immunol. 2012;3.

  50. Jordan AR, Racine RR, Hennig MJP, Lokeshwar VB. The role of CD44 in Disease Pathophysiology and targeted treatment. Front Immunol. 2015;6.

  51. Suck G, Branch DR, Smyth MJ, Miller RG, Vergidis J, Fahim S, et al. KHYG-1, a model for the study of enhanced natural killer cell cytotoxicity. Exp Hematol. 2005;33:1160–71.

    Article  PubMed  CAS  Google Scholar 

  52. Hossain MS, Mawatari S, Fujino T. Plasmalogen-mediated activation of GPCR21 regulates cytolytic activity of NK cells against the Target cells. J Immunol. 2022;209:310–25.

    Article  PubMed  CAS  Google Scholar 

  53. Żychowska M, Woźniak Z, Baran W. Immunohistochemical analysis of the expression of selected cell lineage markers (CD4, CD8, CD68, c-Kit, Foxp3, CD56, CD20) in cutaneous variant of lichen planus. Int J Dermatol. 2021;60:1097–101.

    Article  PubMed  Google Scholar 

  54. Vičić M, Hlača N, Kaštelan M, Brajac I, Sotošek V. Prpić Massari L. Comprehensive Insight into Lichen Planus Immunopathogenesis. Int J Mol Sci. 2023;24:3038.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Aquino-Martinez R, Khosla S, Farr JN, Monroe DG. Periodontal Disease and senescent cells: New players for an old oral Health Problem? Int J Mol Sci. 2020;21:1–22.

    Article  Google Scholar 

  56. Parkinson EK, Prime SS. Oral senescence: from Molecular Biology to Clinical Research. Front Dent Med. 2022;3:6.

    Article  Google Scholar 

  57. Sharma M, Fonseca FP, Hunter KD, Radhakrishnan R. Loss of oral mucosal stem cell markers in oral submucous fibrosis and their reactivation in malignant transformation. Int J Oral Sci. 2020;12:23.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  58. Parkinson EK. Senescence as a modulator of oral squamous cell carcinoma development. Oral Oncol. 2010;46:840–53.

    Article  PubMed  CAS  Google Scholar 

  59. Sahai E, Astsaturov I, Cukierman E, DeNardo DG, Egeblad M, Evans RM, et al. A framework for advancing our understanding of cancer-associated fibroblasts. Nat Rev Cancer. 2020;20:174–86.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  60. Liu T, Han C, Wang S, Fang P, Ma Z, Xu L, et al. Cancer-associated fibroblasts: an emerging target of anti-cancer immunotherapy. J Hematol Oncol. 2019;12:86.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  61. Lee HG, Eun HC. Differences between fibroblasts cultured from oral mucosa and normal skin: implication to wound healing. J Dermatol Sci. 1999;21:176–82.

    Article  PubMed  CAS  Google Scholar 

  62. Chen P, Lingen M, Sonis ST, Walsh-Reitz MM, Toback FG. Role of AMP-18 in oral mucositis. Oral Oncol. 2011;47:831–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  63. Turabelidze A, Guo S, Chung AY, Chen L, Dai Y, Marucha PT, et al. Intrinsic differences between oral and skin keratinocytes. PLoS ONE. 2014;9:e101480.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Pogoda K, Cieśluk M, Deptuła P, Tokajuk G, Piktel E, Król G, et al. Inhomogeneity of stiffness and density of the extracellular matrix within the leukoplakia of human oral mucosa as potential physicochemical factors leading to carcinogenesis. Transl Oncol. 2021;14:101105.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  65. Davis S, Meltzer PS. GEOquery: a bridge between the Gene expression Omnibus (GEO) and BioConductor. Bioinformatics. 2007;23:1846–7.

    Article  PubMed  Google Scholar 

  66. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  67. Mootha VK, Lindgren CM, Eriksson K-F, Subramanian A, Sihag S, Lehar J, et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human Diabetes. Nat Genet. 2003;34:267–73.

    Article  PubMed  CAS  Google Scholar 

  68. Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM, et al. Comprehensive Integration of Single-Cell Data. Cell. 2019;177:1888–1902e21.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  69. Hao Y, Hao S, Andersen-Nissen E, Mauck WM, Zheng S, Butler A, et al. Integrated analysis of multimodal single-cell data. Cell. 2021;184:3573–3587e29.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  70. Borcherding N, Vishwakarma A, Voigt AP, Bellizzi A, Kaplan J, Nepple K et al. Mapping the immune environment in clear cell renal carcinoma by single-cell genomics. Communications Biology 2021 4:1. 2021;4:1–11.

  71. Bunis DG, Andrews J, Fragiadakis GK, Burt TD, Sirota M. dittoSeq: Universal user-friendly single-cell and bulk RNA sequencing visualization Toolkit. Bioinformatics. 2020;36:5535–6.

    Article  PubMed Central  CAS  Google Scholar 

  72. Kanda Y. Investigation of the freely available easy-to-use software EZR for medical statistics. Bone Marrow Transplant. 2013;48:452–8.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

We thank Yumiko Takagi for technical support and Editage (www.editage.com) for the English language editing. This study was supported by JSPS KAKENHI (Grant Numbers 23K18443, 22K19757, 21H03049, and 21H03293), JST SPRING (Grant Number JPMJSP2119), and Takeda Science Foundation.

Author information

Authors and Affiliations

  1. Department of Oral Surgery, Sapporo Medical University School of Medicine, Sapporo, 060-8556, Japan

    Shogo Ijima & Akihiro Miyazaki

  2. Department of Anatomy, Sapporo Medical University School of Medicine, Sapporo, 060-8556, Japan

    Yuki Saito, Kentaro Nagaoka, Taiki Iwamoto & Arisa Kita

  3. Graduate School of Health Sciences, Hokkaido University, Sapporo, 060-0812, Japan

    Sena Yamamoto, Maki Miyajima, Tsukasa Sato & Takako S. Chikenji

Authors
  1. Shogo Ijima
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yuki Saito
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Sena Yamamoto
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Kentaro Nagaoka
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Taiki Iwamoto
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Arisa Kita
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Maki Miyajima
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Tsukasa Sato
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Akihiro Miyazaki
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Takako S. Chikenji
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

SI contributed to the investigation, data collection, and analysis; YS and TSC contributed to the design of the study, conceptualization, investigation, visualization, writing–original draft, review and editing, and supervision; SY, KN, TI, AK, MM, TS, and AM contributed to data collection and analysis. All authors have read and approved the final version of the manuscript and agree with the order of presentation.

Corresponding authors

Correspondence to Yuki Saito or Takako S. Chikenji.

Ethics declarations

Conflict of interest

The authors declare no conflicts of interest.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

: Supplementary Figure 1. UMAP plot of each cell marker gene. Supplementary Figure 2. UMAP plot of mesenchymal cell clusters and ssGSEA. Supplementary Figure 3. P16INK4A expression in the subepithelial layer of oral mucosa sections from patients with OLP and controls. Supplementary Figure 4. Detection of PDGFRα in TIG-118 cells by western blotting. Supplementary Figure 5. NK cells and T cells were activated in OLP. Supplementary Figure 6. Hematoxylin and eosin–stained images of samples from patients with oral lichen planus (OLP) and healthy regions from patients with leukoplakia as controls. Supplementary Table 1. The primer sequences used for the q-PCR

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ijima, S., Saito, Y., Yamamoto, S. et al. Senescence-associated secretory phenotypes in mesenchymal cells contribute to cytotoxic immune response in oral lichen planus. Immun Ageing 20, 72 (2023). https://doi.org/10.1186/s12979-023-00400-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12979-023-00400-5

Keywords

Immunity & Ageing

ISSN: 1742-4933

Contact us