Single-cell RNA sequencing reveals rebalancing of immunological response in patients with periodontitis after non-surgical periodontal therapy
Journal of Translational Medicine volume 20, Article number: 504 (2022)
Periodontitis is a major inflammatory disease of the oral mucosa that is not limited to the oral cavity but also has systemic consequences. Although the importance of chronic periodontitis has been emphasized, the systemic immune response induced by periodontitis and its therapeutic effects remain elusive. Here, we report the transcriptomes of peripheral blood mononuclear cells (PBMCs) from patients with periodontitis.
Using single-cell RNA sequencing, we profiled PBMCs from healthy controls and paired pre- and post-treatment patients with periodontitis. We extracted differentially expressed genes and biological pathways for each cell type and calculated activity scores reflecting cellular characteristics. Intercellular crosstalk was classified into therapy-responsive and -nonresponsive pathways.
We analyzed pan-cellular differentially expressed genes caused by periodontitis and found that most cell types showed a significant increase in CRIP1, which was further supported by the increased levels of plasma CRIP1 observed in patients with periodontitis. In addition, activated cell type-specific ligand-receptor interactions, including the BTLA, IFN-γ, and RESISTIN pathways, were prominent in patients with periodontitis. Both the BTLA and IFN-γ pathways returned to similar levels in healthy controls after periodontal therapy, whereas the RESISTIN pathway was still activated even after therapy.
These data collectively provide insights into the transcriptome changes and molecular interactions that are responsive to periodontal treatment. We identified periodontitis-specific systemic inflammatory indicators and suggest unresolved signals of non-surgical therapy as future therapeutic targets.
Periodontitis is a highly prevalent disease in humans and is the sixth most common human disease . As reclassified at the World Workshop in 2017, periodontitis has a clear distinction from gingivitis in that periodontitis manifests as a loss of alveolar bone and periodontal tissue and has permanent bone defects even after periodontal therapy [1, 2]. In the past, periodontitis has been suggested to be a simple inflammatory disease in the periodontium that is induced by providing a habitat for microbes. However, it is now defined as a complex disease induced by repeated interactions between host modulation of inflammatory and immune activities, periodontal pathogens, and environmental factors . Thus, the susceptibilities to periodontitis are determined by the host response, specifically the magnitude of the inflammatory response and the differential activation of immune pathways rather than the amount and/or the pathogenicity of periodontal bacteria [4, 5].
Although periodontitis is an inflammatory disease localized to the oral cavity, it has been known to be associated with various systemic diseases, such as cardiovascular diseases, diabetes mellitus, and rheumatoid arthritis [6, 7]. Several possibilities have been proposed for the mechanism by which local tissues exert systemic effects. The surface of ulcerated periodontal pockets may allow microbial products, such as lipopolysaccharide (LPS) or proteases, to enter the circulation, resulting in bacteremia . In addition, patients with severe periodontitis have elevated levels of inflammatory mediators, such as IL-1, IL-6, C-reactive protein (CRP), and fibrinogen, in their blood [8,9,10]. Thus, periodontitis can affect the whole body and induce a systemic inflammatory response through leakage of microbial products into the circulatory system [7, 11]. Therefore, treatment for periodontitis is fundamental not only to relieve immediate pain, but also to attenuate the subsequent systemic response.
Recent advances in single-cell RNA sequencing (scRNA-seq) technology have enabled us to understand and elucidate cellular complexities . Since each cell type, especially immune cells, has unique functions, pathways, and regulatory mechanisms, it is essential to characterize cell type-specific genes to identify disease biomarkers and therapeutic targets, which is feasible using scRNA-seq technology. To examine the main factors for chronic inflammation induced by periodontitis, we compared the expression profiles of peripheral blood mononuclear cells (PBMCs) from patients with periodontitis and healthy controls. In addition, to determine whether therapeutic intervention can redirect the immunologic status of periodontitis, we delineated therapy-responsive and -nonresponsive intercellular pathways.
Sample preparation and inclusion criteria
Participants for scRNA-seq analysis and enzyme-linked immunosorbent assay (ELISA) were recruited from the Department of Periodontology, Pusan National University Dental Hospital, under protocols approved by the Institutional Review Board (IRB no. PNUDH-2020-001, and IRB no. PNUDH-2020-032 for scRNA-seq analysis and ELISA, respectively). All participants provided written informed consent. All patients diagnosed with periodontitis stage III were included in this study, and the severity of periodontitis was defined according to the 2017 World Workshop criteria, a new classification scheme for periodontal and peri-implant diseases and conditions . Healthy periodontal status was defined as a pocket depth ≤ 3 mm and no bleeding on probing, with no signs of clinical inflammation, including redness and swelling. The inclusion criteria were subjects with 10 or more teeth and no systemic diseases. In addition, we excluded subjects who had received periodontal treatment within the last 6 months or anti-inflammatory drugs or antibiotics within the last 6 weeks. Those who were previously or currently heavy smokers (which we defined as people who smoked more than one pack per day) were also excluded. All participants with periodontitis received non-surgical periodontal treatment from the same periodontist. Non-surgical periodontal treatment consisted of full-mouth scaling and root planning (SRP) performed in a single appointment under local anesthesia using an ultrasonic device and hand instruments. We attempted to minimize confounding factors that could arise between treatments. Oral hygiene instruction was provided and consisted of the Bass brushing technique and use of an interdental brush and dental floss. The participants were instructed not to use any antimicrobial mouth-rinsing solutions for the duration of the study. The periodontal clinical parameters before and after one month of treatment, as well as the demographic features, are presented in Additional file 1: Table S1.
Blood collection and PBMC isolation
Peripheral blood samples were collected in plastic blood collection tubes containing EDTA. PBMCs for scRNA-seq were isolated using SepMate (Stemcell Technologies Inc.) within 30 min of collection according to the manufacturer’s instructions. Briefly, density gradient medium and diluted blood samples were added to a SepMate tube. After carefully mixing the medium and samples, the tubes were centrifuged at 1200 × g for 10 min. The top layers were poured into a new tube and washed twice with phosphate-buffered saline containing 2% fetal bovine serum. The tubes were then centrifuged at 300 × g for 8 min at room temperature. The collected PBMCs were frozen and stored at − 80 °C until sequencing. Plasma for ELISA was collected after centrifugation at 2,000 × g for 10 min at room temperature and stored at − 80 °C until use.
Enzyme-linked immunosorbent assay (ELISA)
The blood samples in EDTA were centrifuged at 2,000 × g for 10 min at 4 °C, and the supernatant was stored at − 70 °C until analysis. For the measurement of TNF-α, IFITM1 and CRIP1 concentrations, Human TNF-α Uncoated ELISA (Invitrogen, Vienna, Austria), Human Interferon Induced Transmembrane Protein 1 (IFITM1) ELISA (MyBioSource, Vancouver, Canada) and Cysteine-rich protein 1 ELISA (Mybiosoruce, San Diego, CA, USA) kits were used, respectively. Briefly, plasma samples were added to assay plates pre-coated with anti-TNF-α, IFITM1 or CRIP1 antibodies and incubated for 2 h. The plates were then incubated with diluted detection antibody for 1 h, reacted with the substrate solution for 30 min, followed by the addition of the stop solution. All reactions for TNF-α, IFITM1 and CRIP1 were performed at RT and 37 °C, respectively. Standard curves were plotted as control versus the mean optical density at 450 nm. TNFα, IFITM1 and CRIP1 concentrations in each sample were quantified based on a standard curve. Clinical information about the participants whose plasma was used for ELISA is shown in Additional files 2, 3: Tables S2, S3.
Library preparation and sequencing of scRNA
Libraries were prepared using the chromium controller according to the 10× chromium Next GEM Single Cell 3ʹ v3.1 protocol. The cell suspension was mixed with the master mix and loaded with Single Cell 3ʹ v3.1 Gel Beads and Partitioning Oil into a chromium Next GEM chip G. RNA transcripts from single cells were uniquely barcoded and reverse-transcribed within droplets. cDNA molecules were pooled and then subjected to end repair, addition of a single ‘A’ base, and ligation of the adapters. Next, the products were purified and enriched using PCR to create a final cDNA library. Finally, the libraries were sequenced using the Illumina HiSeq platform according to the read length provided in the user guide.
scRNA-seq data pre-processing
Single-cell gene expression data were processed using 10× Genomics Cell Ranger v3.1.0. Raw BCL files from the Illumina sequencing platform were demultiplexed to generate FASTQ files using the ‘cellranger mkfastq’ pipeline. Then, raw FASTQ files were analyzed using the ‘cellranger count’ pipeline. This step includes alignment to the human reference genome (GRCh38, v3.0.0) and measurement of gene expression with a unique molecular identifier (UMI) and cell barcode. Consequently, a cell-by-gene count matrix was generated. To remove low-quality cells, cells with less than 500 UMIs or more than 20,000 UMIs and > 20% mitochondrial genes were filtered out. In addition, we removed cells with fewer than 250 genes or more than 5000 genes, as well as cells with less than 80% complexity (number of genes detected per UMI with log transformation), which could be interpreted as specific cell types, artifacts, or contaminants. In addition, we included genes expressed in more than 0.1% of the cells, not only to eliminate zero counts, but also to prevent genes expressed in a few cells from lowering the average of all other cells. As some samples had a large number of cells (maximum 12,177 cells), possible doublets were estimated using Scrublet, and 3.5% of cells were eliminated (maximum 11,852 cells) .
Analysis of scRNA-seq data
The Seurat (v3.2.2) R package was used to integrate, scale, cluster, and visualize data. The remaining count data were normalized using the SCTransform function based on regularized negative binomial regression on total cellular read counts [14, 15]. After normalization of each sample, the FindIntegrationAnchors and IntegrateData functions were used for integration with the largest dataset among the 12 samples, which was used as the reference dataset. Then, scaling and principal component analysis (PCA) were performed using the ScaleData and RunPCA functions, and the first 30 principal components that were selected depending on the elbow plot were utilized to construct the UMAP dimension reduction and shared nearest-neighbor graph (SNN) using RunUMAP and the FindNeighbors function. Then, cell clusters were distinguished using the graph-based modularity optimization algorithm of the Louvain method for detecting communities with a setting resolution of 0.6, which is able to clearly specify the cell type as well as detect tiny molecular signals ; ultimately, 26 clusters were differentiated. Cell identity markers were identified using the FindAllMarkers function, and genes with a log-fold change threshold > 0.5 and false discovery rate (FDR) < 0.01 were regarded as significant differentially expressed genes (DEGs).
To annotate each cluster into immune cell types, these cell identity markers were supplied into the SingleR (v3.12) R package . SingleR first assigned a single cell into highly matched cell types and then annotated and combined clusters into the most likely cell types. Most of the cells in the same cluster were suggested to be of consistent cell type, usually over 70%, whereas clusters with low consistency were annotated using multiple genes representing cell-specific signatures.
DEG analysis was conducted using edgeR and flexible zero-inflated negative binomial-based wanted variation extraction (ZINB-WaVE), which accounts for zero inflation and over-dispersion [18, 19]. To identify periodontitis-related genes, we compared the expression differences between the three groups as follows: (1) healthy control vs. pre-treatment and (2) pre-treatment vs. post-treatment. In case (2), we compared gene expression with a paired test. DEGs with log2FC > 0.3 and FDR value < 0.05, were extracted and log2FC was recalculated to compare chronic periodontitis status with that of the other two groups.
Gene ontology (GO) analysis
To investigate the biological function of the DEGs detected in the pre-treatment group, we used the R package ‘clusterProfiler’ (ver. 4.0.5) . The upregulated or downregulated genes in each cell type were used to determine the biological function, and the gene set with p < 0.05 and more than one enriched gene was considered to be significant. For genes with more than five biological functions, the top five, in the order of the lowest p-value, were extracted.
We inferred cell‒cell communication based on the expression of ligand-receptor pairs from the manually curated human database CellChatDB. We followed the official workflow for the R package ‘CellChat’ . We created normalized count data from the Seurat object to the CellChat object and pre-processed with the identifyOverExpressedGenes, identifyOverExpressedInteractions, and projectData functions. We generated a ligand‒receptor interaction database comprising secreted signaling, extracellular matrix-receptor, and cell‒cell contact interactions. Then, we calculated the communication probability and inferred cellular communication networks with the computeCommunProb and computeCommunProbPathway functions. Finally, we aggregated the inferred signaling pathways, compared their strength and cellular components, and extracted statistically significant pathways with greater than 1.5-fold signal intensity between the two groups.
Statistical analyses were performed using R version 4.0.3. Kruskal‒Wallis and Wilcoxon rank-sum tests were used to compare cellular composition by setting the significance level at 0.05. PCA plots for clinical information were obtained using the prcomp() function. Receiver operating characteristic (ROC) curves were plotted using the R package ‘ROCit’ .
Single-cell transcriptional landscape of the PBMCs from patients with periodontitis
We isolated PBMCs from four healthy donors and four paired pre-and post-treatment patients for scRNA-seq and obtained 111,213 cells, with an average of 9,268 cells per participant (Fig. 1 A, B). Distinct cell populations were observed after dimensionality reduction and graph-based clustering.
We assigned clusters to the cell type with the highest probability and confirmed the identities using known markers (Fig. 1 C, D). We found 13 major cell compartments that are important for the immune response, and each donor’s PBMCs contained all of these clusters (Fig. 1E). The PBMC populations comprised T cells (60–73%), B cells (4–15%), natural killer (NK) cells (5–23%), monocytes (5–24%), and dendritic cells (1–2%) with normal distributions, regardless of the presence of disease or chronic inflammation-resolved status via periodontal therapy. Although most of the samples showed no significant differences in cellular composition among the three groups, in patients with periodontitis the proportion of CD8+ T cells was slightly lower, while naïve CD8T cells were slightly higher compared to those in healthy controls, but these changes were partially rescued after therapeutic intervention (p = 0.062 and p =, respectively), (Additional file 4: Fig. S1A, B). To verify the clinical inflammatory status of the patients with periodontitis and resolution by periodontal therapy, we performed PCA using six variables: erythrocyte sedimentation rate (ESR) and CRP, which are indicators of inflammation; probing pocket depth (PPD), clinical attachment level (CAL), plaque index (PI), and gingival index (GI), which are indicators of the severity of periodontitis. The healthy controls clustered together but were separated from the pre- and post-treatment groups (Fig. 1 F and Additional file 5: Fig. S2). The pre-treatment group was distinct from the healthy group, while the post-treatment group was more similar to the controls, although the clinical variables of post-treated patients were not comparable to those of the healthy controls. This shows that therapeutic intervention improved the inflammatory status and altered the clinical characteristics of periodontitis.
Gene expression alterations in innate immune cells of patients with chronic periodontitis
Next, we investigated the transcriptome changes in innate immune cells, including monocytes, dendritic cells (DCs), and NK cells. A subtype of DCs, called plasmacytoid dendritic cells (pDCs), which are responsible for antiviral immunity, were excluded as they had no detectable transcriptional changes. The monocytes obtained from Fig. 1 C were re-clustered and classified into three subtypes: classical, non-classical, and intermediate monocytes (Fig. 2 A). The identities of the sub-clusters were determined based on the expression of CD14 and FCGR3A (Additional file 6: Fig. S3A). The DEGs in innate immune cells are shown in Fig. 2B‒D. An increase in cysteine-rich protein 1 (CRIP1) in patients with periodontitis is observed in all types of innate immune cells, and CRIP1 levels after therapy are similar to that in healthy controls, suggesting that CRIP1 may be an important factor in maintaining chronic inflammation. When GO terms were investigated using pre-treatment group-specific DEGs, classical monocytes had low levels of MHC protein complex assembly genes and reduced antigen-presenting activity. In addition, the lymphocyte-activating activity was reduced in non-classical monocytes (Additional file 7: Fig. S4). However, genes responding to LPS stimulation were increased in myeloid DCs (mDCs) of the pre-treatment groups and approached intensity levels similar to the controls following treatment (Fig. 2E, Additional file 8: Fig. S5). These results imply that although monocytes had reduced antigen-presenting function, mDCs actively responded after detection of LPS, a virulence factor of the periodontitis bacteria Porphyromonas gingivalis and could be rescued by periodontal therapy .
Systemic properties of adaptive immune cells found in patients with chronic periodontitis
Next, we investigated the transcriptomic changes in adaptive immune cells. The B, CD4T, and CD8T cell groups from Fig. 1 C were extracted and separated into more specific subtypes. B cells were separated into non-switched memory B cells (non-switched MB), switched memory B cells (switched MB), exhausted B cells, and plasmablasts (Fig. 3 A). CD4+ T cells were divided into four subtypes: T helper 1 (Th1), T helper 17 (Th17), follicular helper T (Tfh), and regulatory T (Treg) cells (Fig. 3B). CD8+ T cells were classified into three subtypes: central memory T (TCM), effector memory T (TEM), and terminal effector T (TTE) cells (Fig. 3 C). The markers for each compartment are described in Additional file 6: Fig. S3B‒D.
The genes that showed altered expression levels in the pre-treatment samples compared to those in both the healthy and post-treatment samples are shown in Fig. 3D‒F. The GO terms of these genes are described in Additional files 9‒11: Fig. S6‒S8. The CD4+ and CD8+ T cell subsets showed decreased cell migration and motility in pre-treatment patients. However, there were no notable changes in the expression of exhausted B cells, plasmablasts, and TEM cells. The lack of differences in TEM cells between the groups was expected since TEM cells were previously reported to be a part of the acute (but not chronic) immune response . This suggests that these nonresponsive cells are not affected by the periodontal bacteria-induced systemic inflammation and are not the main drivers of chronic inflammation.
Common DEGs in immune cells are responsive to systemic chronic inflammation
Across all immune cells, we found common DEGs, such as CRIP1, macrophage migration inhibitory factor (MIF), interferon-induced transmembrane protein 1 (IFITM1), and ribosomal protein S17 (RPS17). DEGs that exhibited contrasting expression patterns in the pre-treatment group compared to those in the healthy and post-treatment groups and their GO terms in the remaining cell types that were not divided into subtypes are shown in Additional files 12, 13: Fig. S9, S10. Notably, in more than half of the cells, all four genes were highly expressed only in pre-treatment patients (Additional File 14: Table S4). The relative expression of these four genes in the healthy group was visualized along with their localization according to clinical inflammatory and periodontal measurements in Additional file 15: Fig. S11. Across all cell types where differential expression was observed, the expression of the four genes increased. CRIP1 was specifically upregulated in the B cell lineage, and IFITM1 was upregulated in the T cell lineage. However, they were all reduced following therapeutic intervention, thus demonstrating a relationship between clinical variables and the expression patterns of these genes.
Of the four genes, CRIP1 and IFITM1 were overexpressed in approximately 97% and 71% of the total immune cells, respectively (Additional file 14: Table S4). An ELISA was conducted on the plasma of the healthy and periodontitis groups to measure changes in CRIP1 and IFITM1 protein concentrations, and compared with those of the inflammatory markers TNF-α, CRP, and ESR. Although there was a similar concentration of the three inflammatory indicators between the healthy and periodontitis groups, we found a difference in CRIP1 and IFITM1 levels (Fig. 3G). In addition, CRIP1 levels positively correlated with TNF-α levels (Fig. 3 H, Additional file 16: Fig. S12). The ROC curves showed that the area under the ROC curve (AUC) for CRIP1 and IFITM1 in the healthy and periodontitis groups was 0.94 and 0.66, respectively (Fig. 3I, both p-values < 0.001). Thus, these results indicate that the levels of these four genes are sensitive to the systemic immune response and that plasma CRIP1 could be a marker for the systemic inflammation caused by periodontitis.
Cell‒cell communication in immune cells of patients with chronic periodontitis
To investigate the intercellular communication between immune cells, which could induce or suppress inflammatory signals, we inferred the cell‒cell interactions within each group. We then compared the signal length of the interactions between the groups (Additional file 17: Fig. S13) and found several chronic inflammation-specific crosstalk proteins: B and T lymphocyte attenuator (BTLA), CCL, CD99, SEMA4, RESISTIN, interferon-γ (IFNG), THBS, and CD48.
Changes in the BTLA and IFNG pathways by periodontal infection are reversible with treatment
The BTLA-mediated pathway was significantly upregulated in pre-treatment patients. The BTLA pathway consists of BTLA‒TNFRSF14 as a ligand‒receptor pair in every group, and naïve B cells predominantly acted as the signal senders (Fig. 4 A‒C). However, unlike the other groups, pDCs also exerted a sending function. In addition, the proportion of pDCs expressing BTLA was higher in the pre-treatment group (Fig. 4D). Based on a study where BTLA‒TNFRSF14 crosstalk promoted Foxp3 expression in T cells via upregulation of CD5, we examined the proportional change in CD5+-expressing T cells and FOXP3+ Tregs. As a result, the ratio of CD5+ T cells and FOXP3+ Tregs increased in the periodontitis pre-treatment group (Additional file 18: Fig. S14). Moreover, in Treg cells, there was a significantly higher propensity for immune tolerance in the pre-treatment group (Fig. 4E). This is consistent with a previous study showing that BTLA in DCs modulates the immunotolerance of T cells . After non-surgical treatment, the BTLA+ pDC ratio and the involvement of pDCs in BTLA interactions, and the tolerance of Tregs, were similar to those of the healthy group.
Conversely, IFNG-mediated interactions only appeared in the periodontal inflammatory group (Fig. 5 A). IFNG secreted by NK cells binds to IFNG receptors (IFNGR) on mDCs and monocytes. The expression level of IFNG and the ratio of IFNG+ NK cells were the highest in the pre-treatment group (Fig. 5B, C). In addition, there were significant increases in pro-inflammatory cytokines released from DCs only in the IFNGR+ mDCs of pre-treatment individuals (Fig. 5D) . This suggests that NK cell-secreted IFNG stimulates mDCs to induce and maintain an inflammatory response in patients with periodontitis but can be ameliorated by non-surgical therapy [27, 28].
Persistent RESISTIN pathway activation despite periodontitis treatment
The RESISTIN pathway was observed in all groups, with the weakest signal intensity in healthy subjects and increased signal intensity in periodontitis patients (Fig. 6 A‒C). In all groups, monocytes interact with all cell types as the signal sender, but in the disease state, mDCs are also able to initiate signal transduction. The RESISTIN pathway consists of interactions with CAP1 and TLR4 receptors, and there was no compositional difference in the ligand-receptor pairs between the groups (Fig. 6D). Examining the expression pattern of the resistin-coding gene RETN in signal sender cells, RETN was most highly expressed in monocytes of the periodontitis group, but levels were significantly decreased after treatment, although it was not completely rescued (Fig. 6E). In addition, the proportion of RETN+ monocytes in the periodontitis group was twice that of the healthy controls, but decreased to levels similar to the healthy group after treatment (Fig. 6 F). In contrast, in mDCs, the expression level of RETN and the proportion of RETN+ cells were only slightly increased (Fig. 6G, H). Therefore, it is noteworthy that, unlike the preceding signaling pathways, treatment intervention did not rescue the changes in RESISTIN pathway activity due to periodontitis. In conclusion, although there was partial rescue in some cells that activated the RESISTIN signaling pathway, mDCs involved in periodontitis were not affected by non-surgical treatment and continued to activate RESISTIN signaling.
In addition, there are other therapy-responsive and nonresponsive intercellular pathways. The CCL and SEMA4 signals, which were elevated in pretreatment patients, were weakened after periodontal treatment, while CD99, THBS, and CD48 signals were maintained. The ligand-receptor pair and cellular networks of these pathways are illustrated in Additional file 19: Fig. S15.
Periodontitis is caused by periodontal dysbiosis of the microbiota, which promotes loss of gingival tissue and bone destruction . The risk of periodontitis affects not only the local site, but also induces a systemic immune response. For example, periodontitis-induced chronic inflammation can disturb typical immunological mechanisms, including bystander activation, amplification of cytokines, epitope spreading, autoantigen complementarity, and activation or inhibition of receptors related to microorganism regulation . Consequently, pathogenesis models leading to extra-oral diseases, such as cardiovascular disease, rheumatoid arthritis, adverse pregnancy outcomes, and respiratory disease, have been widely investigated . However, beyond the pathophysiology of periodontitis, the systemic immune response to periodontal treatment remains unclear. To expand our knowledge, this study shows the periodontal therapy-induced process and identifies genes that are differentially altered according to the disease condition. Moreover, we focused on the intercellular signals between immune cells and further uncovered the interactions that did not respond to periodontal treatment.
In this study, we examined the transcriptional changes and related functions in different immune cells. Each immune cell had pro-inflammatory cytokines that were elevated only in pre-treated periodontitis patients. For instance, CCL3L1 (in NK cells), which enhances the inflammatory response and increases the risk of autoimmune disease , S100A12 (in mDCs), which is a valuable serum inflammatory marker , TNFRSF1B (in mDCs), which regulates TNF-α levels , and CCR7 (in Th17 cells), which promotes joint inflammation , were found to be highly expressed in patients with periodontitis. Unexpectedly, we observed that only a few pro-inflammatory cytokines were altered in periodontitis. In addition, as shown by the ELISA results, the prevalent inflammatory markers only showed slightly increased levels in periodontitis. Altogether, our results indicate that periodontitis induces mild systemic inflammation, consistent with previous studies [7, 35].
Remarkably, the expression of the CRIP1 gene was differentially increased in almost every cell type in the periodontitis samples and the concentration of plasma CRIP1 protein also increased in the periodontitis group. The functional role of CRIP1 has been reported in a previous study, where transgenic mice overexpressing Crip1 produced high concentrations of IL-6 and IL-10 after LPS treatment . In addition, its role has been studied in diverse cancer types, suggesting that CRIP1 might be a risk factor for chronic inflammation and systemically affects numerous tissues and organs [37,38,39,40]. Moreover, CRIP1 levels reflect periodontitis-induced inflammation more sensitively than typical indicators. Taken together, the upregulation of CRIP1 may be a predictive marker for systemic inflammation induced by periodontitis.
In the intercellular network, BTLA and IFNG signals are strongly activated in the circulatory system during periodontitis. The BTLA‒TNFRSF14 interaction has been reported to regulate T-cell responses and increase disease susceptibility. Here, we confirmed that pDCs formed distinct signaling networks, and the immunological tolerance of T cells was enhanced in periodontal disease . It is speculated that this is because BTLA-expressing DCs induce and modulate subsequent immune responses of Treg cells by promoting Foxp3 and CD5 expression . Additionally, mDCs and monocytes are involved in IFNG network. In particular, mDCs expressing IFNGR secrete high levels of inflammatory cytokines, which may lead to abnormal localization or activation of the immune response . Similarly, the CCL pathway, which is known to augment the production of IFNG from the CCL5‒CCR1 interaction, was strongly upregulated in periodontitis patients, which supports the higher activity of INFG . From these results, we can infer that patients with periodontitis exhibit antigen-specific tolerance and an imbalanced immune response. However, these intercellular networks can be resolved using nonsurgical periodontal treatment.
In contrast, the RESISTIN pathway was consistently active regardless of non-surgical treatment. Unlike rodent adipocytes, human RESISTIN is predominantly expressed in macrophages and is produced in response to inflammatory stimuli [43,44,45,46]. As in previous research, most RESISTIN transcripts were secreted by monocytes from healthy donors and patients with periodontitis. However, interaction with mDCs appeared in patients with periodontitis and was sustained even after periodontal therapy. In addition, RESISTIN binds primarily to the adenylyl cyclase-associated protein 1 (CAP1) receptor, which is known to upregulate circulating AMP (cAMP) concentrations, protein kinase A (PKA) activity, and NF-κB-associated transcription of inflammatory cytokines . Elevated RESISTIN levels have also been demonstrated in various inflammatory diseases, such as rheumatoid arthritis (RA), inflammatory bowel disease (IBD), type II diabetes, and sepsis. [48,49,50]. Thus, RESISTIN signaling may systemically increase susceptibility to diverse diseases. Moreover, the THBS network, whose ligand is highly expressed after LPS stimulation, increases the risk of inflammation [51, 52], and the CD48 network, which is known to activate NK cells and contribute to diverse autoimmune diseases, is maintained even after periodontal therapy [53, 54]. All this evidence indicates that typical periodontal treatment cannot be the sole strategy for complete immune recovery and necessitates consideration of other therapeutic approaches.
Nevertheless, this study had some limitations. First, there is a gap between the ages of healthy individuals and patients with periodontitis. Age is one of the major factors that enhance the prevalence of periodontitis and chronic inflammation [55,56,57]. Although all subjects in this study were middle-aged, there was a 6-year difference between the groups on average. This could indicate that the patients had higher exposure to a known risk factor for chronic inflammation; thus, further research under strict control of external factors is required. Next, cryopreserved PBMCs were used to generate RNA sequence profiles. Despite the general use of frozen samples in scRNA-seq, the identification of immune cell subpopulations that have not been fully characterized may be challenging, considering the comparatively lower detection of UMIs in frozen PBMCs than in fresh PBMCs [58, 59]. Therefore, fresh samples can assist with producing more precise results.
In summary, our current study reports transcriptome changes in periodontitis at single-cell resolution. In this study, we revealed the systemic immunological effects of periodontitis and identified periodontitis-specific predictors of inflammation. In addition, based on the immune pathways that were responsive to therapy, we found nonresponsive pathways to treatment that can increase the risk of comorbidity. Therefore, we suggest target pathways that can resolve chronic systemic inflammation but these therapeutics targets will require experimental validation.
The datasets generated during this study are available at GEO (https://www.ncbi.nlm.nih.gov/geo/) under accession number GSE174609.
area under the ROC curve
B and T lymphocyte attenuator
clinical attachment level
cysteine-rich protein 1
differentially expressed gene
Enzyme-linked immunosorbent assay
erythrocyte sedimentation rate
interferon-induced transmembrane protein 1
myeloid dendritic cell
macrophage migration inhibitory factor
- non-switched MB:
non-switched memory B cells
peripheral blood mononuclear cells
plasmacytoid dendritic cell
probing pocket depth
- ROC curve:
receiver operating characteristic curve
ribosomal protein S17
single-cell RNA sequencing
- switched MB:
switched memory B cells
central memory T
effector memory T
follicular helper T
T helper 1
T helper 17
terminal effector T
Genco RJ, Sanz M. Clinical and public health implications of periodontal and systemic diseases: An overview. Periodontol 2000 2020, 83:7–13.
Caton JG, Armitage G, Berglundh T, Chapple ILC, Jepsen S, Kornman KS, Mealey BL, Papapanou PN, Sanz M, Tonetti MS. A new classification scheme for periodontal and peri-implant diseases and conditions - Introduction and key changes from the 1999 classification. J Clin Periodontol. 2018;45(Suppl 20):1–8.
Balta MG, Papathanasiou E, Blix IJ, Van Dyke TE. Host Modulation and Treatment of Periodontal Disease. J Dent Res. 2021;100:798–809.
Preshaw PM. Host modulation therapy with anti-inflammatory agents. Periodontol 2000 2018, 76:131–149.
Van Dyke TE. Shifting the paradigm from inhibitors of inflammation to resolvers of inflammation in periodontitis. J Periodontol. 2020;91(Suppl 1:):19–25.
Seymour G, Ford P, Cullinan M, Leishman S, Yamazaki K. Relationship between periodontal infections and systemic disease. Clin Microbiol Infect. 2007;13:3–10.
Hajishengallis G, Chavakis T. Local and systemic mechanisms linking periodontal disease and inflammatory comorbidities. Nat Rev Immunol. 2021;21:426–40.
D’Aiuto F, Orlandi M, Gunsolley JC. Evidence that periodontal treatment improves biomarkers and CVD outcomes. J Clin Periodontol. 2013;40:85–105.
Genco RJ, Van Dyke TE. Reducing the risk of CVD in patients with periodontitis. Nat Reviews Cardiol. 2010;7:479–80.
Schenkein HA, Papapanou PN, Genco R, Sanz M. Mechanisms underlying the association between periodontitis and atherosclerotic disease. Periodontol 2000. 2020;83:90–106.
Bui FQ, Almeida-da-Silva CLC, Huynh B, Trinh A, Liu J, Woodward J, Asadi H, Ojcius DM. Association between periodontal pathogens and systemic disease. Biomedical J. 2019;42:27–35.
Papalexi E, Satija R. Single-cell RNA sequencing to explore immune cell heterogeneity. Nat Rev Immunol. 2018;18:35–45.
Wolock SL, Lopez R, Klein AM. Scrublet: Computational Identification of Cell Doublets in Single-Cell Transcriptomic Data. Cell Syst. 2019;8:281–91 e289.
Butler A, Hoffman P, Smibert P, Papalexi E, Satija R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol. 2018;36:411–20.
Hafemeister C, Satija R. Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression. Genome Biol. 2019;20:1–15.
Waltman L, Van Eck NJ. A smart local moving algorithm for large-scale modularity-based community detection. Eur Phys J B. 2013;86:1–14.
Aran D, Looney AP, Liu L, Wu E, Fong V, Hsu A, Chak S, Naikawadi RP, Wolters PJ, Abate AR, et al. Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage. Nat Immunol. 2019;20:163–72.
McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012;40:4288–97.
Risso D, Perraudeau F, Gribkova S, Dudoit S, Vert JP. A general and flexible method for signal extraction from single-cell RNA-seq data. Nat Commun. 2018;9:284.
Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7.
Jin S, Guerrero-Juarez CF, Zhang L, Chang I, Ramos R, Kuan C-H, Myung P, Plikus MV, Nie Q. Inference and analysis of cell-cell communication using CellChat. Nat Commun. 2021;12:1–20.
Khan MRAA: ROCit-An R Package for Performance Assessment of Binary Classifier with Visualization. 2019.
Wang PL, Ohura K. Porphyromonas gingivalis lipopolysaccharide signaling in gingival fibroblasts-CD14 and Toll-like receptors. Crit Rev Oral Biol Med. 2002;13:132–42.
Harris NL, Watt V, Ronchese F, Le Gros G. Differential T cell function and fate in lymph node and nonlymphoid tissues. J Exp Med. 2002;195:317–26.
Jones A, Bourque J, Kuehm L, Opejin A, Teague RM, Gross C, Hawiger D. Immunomodulatory Functions of BTLA and HVEM Govern Induction of Extrathymic Regulatory T Cells and Tolerance by Dendritic Cells. Immunity. 2016;45:1066–77.
Blanco P, Palucka AK, Pascual V, Banchereau J. Dendritic cells and cytokines in human inflammatory and autoimmune diseases. Cytokine Growth Factor Rev. 2008;19:41–52.
Clavijo-Salomon MA, Salcedo R, Roy S, das Neves RX, Dzutsev A, Sales-Campos H, Borbely KS, Silla L, Orange JS, Mace EM, et al. Human NK cells prime inflammatory DC precursors to induce Tc17 differentiation. Blood Adv. 2020;4:3990–4006.
Goldszmid RS, Caspar P, Rivollier A, White S, Dzutsev A, Hieny S, Kelsall B, Trinchieri G, Sher A. NK cell-derived interferon-gamma orchestrates cellular dynamics and the differentiation of monocytes into dendritic cells at the site of infection. Immunity. 2012;36:1047–59.
Suarez LJ, Garzon H, Arboleda S, Rodriguez A. Oral Dysbiosis and Autoimmunity: From Local Periodontal Responses to an Imbalanced Systemic Immunity. A Review. Front Immunol. 2020;11:591255.
Hajishengallis G. Periodontitis: from microbial immune subversion to systemic inflammation. Nat Rev Immunol. 2015;15:30–44.
McKinney C, Merriman ME, Chapman PT, Gow PJ, Harrison AA, Highton J, Jones PB, McLean L, O’Donnell JL, Pokorny V, et al. Evidence for an influence of chemokine ligand 3-like 1 (CCL3L1) gene copy number on susceptibility to rheumatoid arthritis. Ann Rheum Dis. 2008;67:409–13.
Meijer B, Gearry RB, Day AS. The role of S100A12 as a systemic marker of inflammation. Int J Inflam. 2012;2012:907078.
Martin EM, Remke A, Pfeifer E, Polz J, Pietryga-Krieger A, Steffens-Weber D, Freudenberg MA, Mostbock S, Mannel DN. TNFR2 maintains adequate IL-12 production by dendritic cells in inflammatory responses by regulating endogenous TNF levels. Innate Immun. 2014;20:712–20.
Van Raemdonck K, Umar S, Palasiewicz K, Volkov S, Volin MV, Arami S, Chang HJ, Zanotti B, Sweiss N, Shahrara S. CCL21/CCR7 signaling in macrophages promotes joint inflammation and Th17-mediated osteoclast formation in rheumatoid arthritis. Cell Mol Life Sci. 2020;77:1387–99.
Chen Y, Liu S, Leng SX. Chronic Low-grade Inflammatory Phenotype (CLIP) and Senescent Immune Dysregulation. Clin Ther. 2019;41:400–9.
Lanningham-Foster L, Green CL, Langkamp-Henken B, Davis BA, Nguyen KT, Bender BS, Cousins RJ. Overexpression of CRIP in transgenic mice alters cytokine patterns and the immune response. Am J Physiology-Endocrinology Metabolism. 2002;282:E1197–203.
Groene J, Mansmann U, Meister R, Staub E, Roepcke S, Heinze M, Klaman I, Brümmendorf T, Hermann K, Loddenkemper C. Transcriptional census of 36 microdissected colorectal cancers yields a gene signature to distinguish UICC II and III. Int J Cancer. 2006;119:1829–36.
Ludyga N, Englert S, Pflieger K, Rauser S, Braselmann H, Walch A, Auer G, Höfler H, Aubele M. The impact of cysteine-rich intestinal protein 1 (CRIP1) in human breast cancer. Mol Cancer. 2013;12:1–10.
Baumhoer D, Elsner M, Smida J, Zillmer S, Rauser S, Schoene C, Balluff B, Bielack S, Jundt G, Walch A, Nathrath M. CRIP1 expression is correlated with a favorable outcome and less metastases in osteosarcoma patients. Oncotarget. 2011;2:970–5.
Balluff B, Rauser S, Meding S, Elsner M, Schöne C, Feuchtinger A, Schuhmacher C, Novotny A, Jütting U, Maccarrone G. MALDI imaging identifies prognostic seven-protein signature of novel tissue markers in intestinal-type gastric cancer. Am J Pathol. 2011;179:2720–9.
Gavrieli M, Sedy J, Nelson CA, Murphy KM. BTLA and HVEM cross talk regulates inhibition and costimulation. Adv Immunol. 2006;92:157–85.
Ness TL, Carpenter KJ, Ewing JL, Gerard CJ, Hogaboam CM, Kunkel SL. CCR1 and CC chemokine ligand 5 interactions exacerbate innate immune responses during sepsis. J Immunol. 2004;173:6938–48.
Kunnari AM, Savolainen ER, Ukkola OH, Kesaniemi YA, Jokela MA. The expression of human resistin in different leucocyte lineages is modulated by LPS and TNFalpha. Regul Pept. 2009;157:57–63.
Kaser S, Kaser A, Sandhofer A, Ebenbichler CF, Tilg H, Patsch JR. Resistin messenger-RNA expression is increased by proinflammatory cytokines in vitro. Biochem Biophys Res Commun. 2003;309:286–90.
Lu SC, Shieh WY, Chen CY, Hsu SC, Chen HL. Lipopolysaccharide increases resistin gene expression in vivo and in vitro. FEBS Lett. 2002;530:158–62.
Schwartz DR, Lazar MA. Human resistin: found in translation from mouse to man. Trends Endocrinol Metab. 2011;22:259–65.
Lee S, Lee HC, Kwon YW, Lee SE, Cho Y, Kim J, Lee S, Kim JY, Lee J, Yang HM, et al. Adenylyl cyclase-associated protein 1 is a receptor for human resistin and mediates inflammatory actions of human monocytes. Cell Metab. 2014;19:484–97.
Sunden-Cullberg J, Nystrom T, Lee ML, Mullins GE, Tokics L, Andersson J, Norrby-Teglund A, Treutiger CJ. Pronounced elevation of resistin correlates with severity of disease in severe sepsis and septic shock. Crit Care Med. 2007;35:1536–42.
Filkova M, Haluzik M, Gay S, Senolt L. The role of resistin as a regulator of inflammation: Implications for various human pathologies. Clin Immunol. 2009;133:157–70.
Chen BH, Song Y, Ding EL, Manson JE, Roberts CK, Rifai N, Buring JE, Gaziano JM, Liu S. Association of resistin promoter polymorphisms with plasma resistin levels and type 2 diabetes in women and men. Int J Mol Epidemiol Genet. 2010;1:167–74.
Fordham JB, Hua J, Morwood SR, Schewitz-Bowers LP, Copland DA, Dick AD, Nicholson LB. Environmental conditioning in the control of macrophage thrombospondin-1 production. Sci Rep. 2012;2:1–9.
Gao L, Yang TT, Zhang JS, Liu HX, Cai DC, Wang LT, Wang J, Li XW, Gao K, Zhang SY, et al. THBS1/CD47 Modulates the Interaction of gamma-Catenin With E-Cadherin and Participates in Epithelial-Mesenchymal Transformation in Lipid Nephrotoxicity. Front Cell Dev Biol. 2020;8:601521.
Messmer B, Eissmann P, Stark S, Watzl C. CD48 stimulation by 2B4 (CD244)-expressing targets activates human NK cells. J Immunol. 2006;176:4646–50.
Sun L, Gang X, Li Z, Zhao X, Zhou T, Zhang S, Wang G. Advances in understanding the roles of CD244 (SLAMF4) in immune regulation and associated diseases. Front Immunol. 2021;12:648182.
Chung HY, Kim DH, Lee EK, Chung KW, Chung S, Lee B, Seo AY, Chung JH, Jung YS, Im E, et al. Redefining Chronic Inflammation in Aging and Age-Related Diseases: Proposal of the Senoinflammation Concept. Aging Dis. 2019;10:367–82.
Clark D, Kotronia E, Ramsay SE. Frailty, aging, and periodontal disease: Basic biologic considerations. Periodontol 2000. 2021;87:143–56.
Van Dyke TE, Sheilesh D. Risk factors for periodontitis. J Int Acad Periodontol. 2005;7:3–7.
Chen J, Cheung F, Shi R, Zhou H, Lu W, Consortium CHI. PBMC fixation and processing for Chromium single-cell RNA sequencing. J Transl Med. 2018;16:198.
Yao C, Bora SA, Chen P, Goodridge HS, Gharib SA. Sample processing and single cell RNA-sequencing of peripheral blood immune cells from COVID-19 patients. STAR Protoc. 2021;2:100582.
The authors thank all members of ROKIT GENOMICS for their helpful comments and discussions.
This work was supported by the National Research Foundation of Korea (NRF-2020R1A2C1005203, NRF-2020R1C1C1003741, and NRF-2018R1A5A2023879) and the Korea Medical Institute (KMI).
Ethics approval and consent to participate
Eight study participants were recruited from the Department of Periodontology, Pusan National University Dental Hospital, under protocols approved by the Institutional Review Board (IRB no. PNUDH-2020-001). All participants provided written informed consent.
Consent for publication
Patents pending by the authors and their institutions.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Lee, H., Joo, JY., Sohn, D.H. et al. Single-cell RNA sequencing reveals rebalancing of immunological response in patients with periodontitis after non-surgical periodontal therapy. J Transl Med 20, 504 (2022). https://doi.org/10.1186/s12967-022-03702-2