Immune Cell Infiltration Characteristics of Pigmented Villous Nodular Synovitis and Prediction of Potential Diagnostic Markers Based on Bioinformatics

Zhang, Jun; Li, Bin; Zhao, Boming; Qi, Yongjian; Chen, Liaobin; Chen, Jun; Chen, Biao

doi:https://doi.org/10.1155/2022/8708692

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 8708692 | https://doi.org/10.1155/2022/8708692

Immune Cell Infiltration Characteristics of Pigmented Villous Nodular Synovitis and Prediction of Potential Diagnostic Markers Based on Bioinformatics

Jun Zhang,¹Bin Li,¹Boming Zhao,¹Yongjian Qi,²Liaobin Chen,¹Jun Chen,¹and Biao Chen¹

Academic Editor: Arun K. Sharma

Received21 Feb 2022

Accepted21 May 2022

Published07 Jun 2022

Abstract

Background. Pigmented villous nodular synovitis (PVNS) is a tumor-like proliferative disease characterized by impairment of daily activities, decreased quality of life, and a high recurrence rate. However, the specific pathological mechanisms are still ill-defined and controversial. The purpose of this study was to define potential diagnostic markers and evaluate immune cell infiltration in the pathogenesis of PVNS. Method. The expression profile of GSE3698 was reanalyzed in the Gene Expression Omnibus (GEO) database. First, differentially expressed genes (DEGs) were identified using the R package “limma” and analyzed by Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment. Next, the DEGs were imported into the STRING database and Cytoscape to construct a protein–protein interaction (PPI) network. Then, cytoHubba and ROC curve analyses were used to determine potential diagnostic biomarkers of PVNS. Finally, we used CIBERSORT to estimate the proportions of 22 immune cell subtypes in PVNS and analyzed the correlation between diagnostic markers and infiltrating immune cells. Result. We found 139 DEGs (including 93 upregulated genes and 46 downregulated genes). TYROBP, FCER1G, LAPTM5, and HLA-DPB1 were identified as potential diagnostic biomarkers of PVNS. Immune cell infiltration analysis indicated that neutrophils and M2 macrophages might be associated with the genesis and progression of PVNS. Furthermore, our correlation analysis of diagnostic markers and infiltrating immune cells found that TYROBP, FCER1G, LAPTM5, and HLA-DPB1 were positively correlated with M2 macrophage infiltration and that neutrophils, TYROBP, FCER1G, and LAPTM5 were negatively correlated with plasma cell infiltration. Conclusions. We identified TYROBP, FCER1G, LAPTM5, and HLA-DPB1 as potential diagnostic markers for PVNS and concluded that immune cell infiltration plays an important role in the genesis and progression of PVNS.

1. Introduction

Pigmented villous nodular synovitis (PVNS) is a relatively rare, tumor-like proliferative disease in the synovium of the joint tissue [1]. PVNS favors the knee and hip, mainly affecting young adults between 20 and 40 years, with an incidence estimated to be between 1.8 per million. Clinically, PVNS manifests as a locally destructive process, leading to functional decline and decreased quality of life [2–4]. Moreover, the local recurrence rate is still high even after complete surgical resection, especially for the diffuse type [5]. If not treated in time, PVNS will lead to joint degenerative changes, moderate to severe joint deformity, and cortical bone destruction. Eventually, patients will require total joint arthroplasty or even amputation [6]. Currently, the primary therapeutic approach is surgical resection of the entire pathological synovial tissue to decrease joint pain and destruction, followed by postoperative chemotherapy to reduce the risk of recurrence [7]. However, it is still controversial whether PVNS is an inflammatory reaction or a neoplastic hyperplasia disease [8, 9]. Considering the high recurrence rate, extremely aggressive invasion, and destruction caused by PVNS, it is crucial to elucidate the molecular mechanisms, understand the underlying pathological mechanisms, and identify novel early diagnostic markers.

Previous studies have shown that the proportions of T cells, natural killer (NK) cells, and natural killer T (NKT) cells are significantly increased in the peripheral blood of PVNS patients. In addition, a significant increase in osteoclastogenesis and macrophage activation was observed in the local synovial membrane [10]. Additional studies have also suggested that macrophages and macrophage-like cells are involved in the occurrence and development of PVNS [11]. West et al. [12] showed that a minority chromosome involving the colony-stimulating factor 1 (CSF1) gene ectopic in synovial tissue resulted in the accumulation of monocytes and macrophages. These observations indicated that infiltrating immune cells might play an essential role in the pathological processes of PVNS. Therefore, it is critical to evaluate whether there is immune cell infiltration in PVNS and determine the intrinsic differences in infiltrating immune cell types for PVNS from an immune system perspective, which could help elucidate the potential mechanism of disease and suggest possible novel immunotherapy targets. With the rapid development of science and technology, innovative bioinformatic approaches have been developed that offer a powerful platform to look for potential diagnostic markers and evaluate immune cell infiltration [13]. CIBERSORT is a popular analytical tool that utilizes gene expression-profiling datasets to assess the expression of immune cells and analyze various immune cell proportions in the synovial tissue of PVNS and osteoarthritis (OA) [14]. At present, CIBERSORT has been widely used for research in myocarditis [15], breast cancer [13], atherosclerosis-related cardiovascular diseases [16], and many other conditions. However, to our knowledge, no studies on the use of CIBERSORT to analyze immune cell infiltration in PVNS have been reported thus far.

In this study, we first downloaded the original microarray dataset from the Gene Expression Omnibus (GEO) database and identified the differentially expressed genes implicated in PVNS. In addition, Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to investigate the biological function of DEGs. Subsequently, the STRING database, Cytoscape, and receiver operating characteristic (ROC) curves were further used to filter and identify potential diagnostic biomarkers in PVNS. Moreover, in this study, for the first time, we applied the CIBERSORT method to analyze the differences in 22 immune cell subsets in PVNS using gene expression profiles. Finally, we investigated the relationship between the identified diagnostic marker genes and infiltrating immune cells to help better understand the pathological immune mechanisms of PVNS.

2. Materials and Methods

2.1. Data Source and Data Preprocessing

We downloaded the GSE3698 RAW dataset from the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo), which contained 11 cases of PVNS synovial tissue and 19 cases of OA synovial tissue. The gene annotation platform of the GSE3698 dataset was GPL3050. First, the data were preprocessed using the “limma” package of R/Bioconductor [17]. Subsequently, probes were annotated using an annotation platform file. Probe sets without known annotation, mapped to multiple Gene IDs, or did not map to any Gene ID were removed. In cases where a gene corresponds to a plurality of probe sets, the maximum value was used as the gene expression value. Eventually, 30 samples and 9647 genes were processed for subsequent analysis.

2.2. Differentially Expressed Gene (DEG) Identifications

The R ackage “limma” was used to screen the differentially expressed genes, with absolute log2 fold change and as threshold. Volcano plots were produced using the R package “ggplot2,” and heat maps were produced using the R package “pheatmap” for differentially expressed genes. According to the logFC values, the top 10 genes upregulated and downregulated were selected to construct the heat map.

2.3. Pathway Enrichment Analyses

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed using Database for Annotation, Visualization, and Integrated Discovery (DAVID) tool (version 6.8) [18]. The GO terms and KEGG pathways with and were considered statistically significant and visualized through the R package “Goplot.”

2.4. Protein-Protein Interaction (PPI) Network Construction and Hub Gene Identification

The Search Tool for Retrieving Interacting Genes (STRING), which contains a database of known and predicted protein-protein interactions (PPI), was used to construct PPI networks of differentially expressed genes (confidence score ) [19]. The gene clusters were then imported into Cytoscape (version 3.9.0.) to visualize the networks [20]. Three analysis methods were used to screen out the top 13 hub genes, respectively, namely, Degree, Maximum Neighborhood Component (MNC), and Maximal Clique Centrality (MCC) in the Cytoscape plugin cytoHubba [21]. In addition, we looked for shared Hub genes of selection between three analysis methods and were visualized by Venn diagrams in R package “VennDiagram.”

2.5. ROC Curve Verification of Hub Genes

A receiver operating characteristic (ROC) curve was drawn, and an area under the curve (AUC) was calculated in R package “pROC” according to gene expression profile data of hub genes. In general, we considered AUC values of >0.9 indicated high diagnostic performance [22].

2.6. Immune Cell Infiltration Analysis

First, the GSE3698 expression matrix data and LM22 matrix transform gene expression were uploaded into R (version 4.1.2); the source code of CIBERSORT and the LM22 matrix transform gene expression were obtained from the CIBERSORTx website (https://cibersortx.stanford.edu/). Filtered out samples with and removed noninfiltrating immune cells. The differences in immune cell infiltration between PVNS and OA samples were obtained in R package “vioplot.” The proportion of infiltrating immune cells in the sample is plotted in R package “ggplot2,” and the R package “corrplot” draws a heat map of infiltrating immune cell correlations.

2.7. Correlation Analysis between Diagnostic Markers and Infiltrating Immune Cells

We performed Spearman correlation analysis between the expression matrix data of diagnostic markers and infiltrating immune cells. Data were visualized using the ggplot2 package.

3. Results

3.1. DEG Identifications

We identified a total of 139 DEGs (including 93 upregulated DEGs and 46 downregulated DEGs) in PVNS synovial tissue compared with OA synovial tissue, and a volcano plot is shown in Figure 1(a). A heat map of DEGs is shown in Figure 1(b). In addition, the heat map for the top 10 genes upregulated and downregulated according to the logFC values is shown in Figure 1(c).

(a)

(b)

(c)

3.2. Pathway Enrichment Analyses

We further performed GO and KEGG pathway enrichment analyses on DEGs using the DAVID database to explore the potential pathological process of PVNS. Functional enrichment analysis revealed that 44 biological process (BP), 30 cellular component (CC), 18 molecular function (MF) (Table S1), and 22 KEGG pathways (Table S2) were statistically significant. We screened out the top three GO-BP, GO-CC, and GO-MF pathways and drew a circle plot showing the number of genes enriched in each GO pathway (Figure 2(a)). The top three GO-BPs were antigen processing and presentation of exogenous peptide antigen via MHC class II, immune response, antigen processing, and presentation of peptide or polysaccharide antigen via MHC class II. The top three GO-CC terms were extracellular exosome, MHC class II protein complex, and lysosomal membrane; the top three GO-MF terms were MHC class II receptor activity, collagen binding, and MHC class II protein complex binding. Figure 2(b) shows the top three KEGG pathways and the related genes enriched in each pathway. The top three KEGG pathways were rheumatoid arthritis, asthma, antigen processing, and presentation.

(a)

(b)

3.3. PPI Network Construction and Hub Gene Identification

We uploaded 139 DEGs into the STRING database to construct PPI networks visualized by the Cytoscape software. The PPI network showed 90 nodes, where each node represents a DEG, including 25 upregulated genes and 65 downregulated genes (Figure 3). The top 13 genome modules were obtained from the PPI network through three methods: MCC, MNC, and degree in the cytoHubba plugin of the Cytoscape software. We screened 8 hub genes (TYROBP, FCER1G, LAPTM5, HLA-DPB1, CD74, HLA-DPA1, HLA-DQB2, PLEK) with upregulated expression by taking the intersection of genome modules (Figure 4).

(a)

(b)

(c)

(d)

Figure 4

(a) Gene set network diagram of MCC algorithm in cytoHubba, the redder the color, the higher the ranking, and the yellower the color, the lower the ranking. (b) The gene set network diagram of the MNC algorithm in cytoHubba, the redder the color, the higher the ranking, and the yellower the color, the lower the ranking. (c) The gene set network diagram of the degree algorithm in cytoHubba, the redder the color in the box, the higher the ranking, and the yellower the color, the lower the ranking. (d) Venn diagrams of gene sets obtained from the three algorithms.

3.4. ROC Curve Verification of Hub Genes

We then assessed the diagnostic sensitivity of the hub genes by ROC curve analysis. The area under the curve (AUC) for TYROBP was 0.962 (95% confidence interval (CI): 0.904-1.000) (Figure 5(a)); the area under the curve (AUC) for FCER1G was 0.957 (95% confidence interval (CI): 0.890-1.000) (Figure 5(b)); the area under the curve (AUC) for LAPTM5 was 0.947 (95% confidence interval (CI): 0.851-1.000) (Figure 5(c)); the area under the curve (AUC) for HLA-DPB1 was 0.909 (95% confidence interval (CI): 0.796-1.000) (Figure 5(d)); the area under the curve (AUC) for CD74 was 0.871 (95% confidence interval (CI): 0.729-1.000) (Figure 5(e)); the area under the curve (AUC) for HLA-DPA1 was 0.727 (95% confidence interval (CI): 0.508-0.957) (Figure 5(f)); the area under the curve (AUC) for HLA-DQB2 was 0.842 (95% confidence interval (CI): 0.694-0.990) (Figure 5(g)); the area under the curve (AUC) for PLEK was 0.866 (95% confidence interval (CI): 0.704-1.000) (Figure 5(h)). Based on the high diagnostic value of [22], we identified TYROBP, FCER1G, LAPTM5, and HLA-DPB1 as potential diagnostic markers in PVNS.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

3.5. Immune Cell Infiltration Analysis

A violin plot of differences in immune cell infiltration suggested that 19 types of immune cell infiltration were present in PVNS samples. At the same time, the infiltration of monocytes, naive CD4+ T cells, and M0 macrophages was not detected in PVNS. M2 macrophages () and neutrophils () were expressed at significantly higher levels in PVNS than in OA (Figure 6(a)). Figure 6(b) shows the individual sample’s proportion of infiltrated immune cells. A correlation heat map of the 19 types of infiltrating immune cells revealed that regulatory T cells had the greatest positive correlation with resting mast cells (). CD8+ T cells had a positive correlation with activated CD4+ memory T cells (). Resting dendritic cells had a positive correlation with activated mast cells (). M1 macrophages had the greatest negative correlation with activated dendritic cells (). Gamma-delta T cells had a negative correlation with CD8+ T cells () and follicular helper T cells (). Plasma cells had a negative correlation with neutrophils (), M2 macrophages (), and CD8+ T cells (). Naive B cells had a negative correlation () with memory B cells (Figure 6(c)).

(a)

(b)

(c)

Figure 6

Assessment and visualization of immune cell infiltration. (a) Violin plot of the proportion of 22 immune cells, with red markers representing significant infiltration differences between the two groups. (b) The proportion of 22 immune cells between-sample. (c) Correlation heat map of 19 types of immune cells. The size of the colored square numbers represents the strength of the correlation; red represents positive correlation, and blue represents negative correlation. The darker the color, the stronger the correlation.

3.6. Correlation Analysis between Diagnostic Markers and Infiltrating Immune Cells

A correlation analysis was performed among four diagnostic markers, TYROBP, FCER1G, LAPTM5, and HLA-DPB1, and 19 types of infiltrating immune cells. The results showed that LAPTM5 was positively correlated with M2 macrophages (, ) and neutrophils (, ) and negatively correlated with plasma cells (, ); FCER1G was positively correlated with M2 macrophages (, ) and neutrophils (, ) and negatively correlated with plasma cells (, ); TYROBP was positively correlated with M2 macrophages (, ) and neutrophils and negatively correlated with plasma cells (, ); HLA-DPB1 was positively correlated with M2 macrophages (, ) and neutrophils (, ) (Figure 7).

(a)

(b)

(c)

(d)

Figure 7

Correlation between LAPTM5, FCER1G, TYROBP, and HLA-DPB1 with infiltrating immune cells: (a) correlation between LAPTM5 and infiltrating immune cells; (b) correlation between FCER1G and infiltrating immune cells; (c) correlation between TYROBP and infiltrating immune cells; (d) correlation between HLA-DPB1 and infiltrating immune cells. The size of the dot represents the strength of the correlation between genes and immune cells; the larger the dot, the stronger the correlation, and the smaller the dot, the weaker the correlation. The color of the point represents the value; the greener the color, the larger the value, and the redder the color, the smaller the value. was considered statistically significant.

4. Discussion

PVNS is a relatively rare tumor-like proliferative disease of the knee joint. As the disease progresses, local manifestations include artificial knee effusions, hemosiderin deposition, synovial hyperplasia, and bone erosions. Moderate to severe joint destruction and stiffness due to recurrent hemarthrosis may also occur [2]. If the articular cartilage is severely eroded, cortical bone destruction will necessitate total synovectomy, total knee replacement, or even amputation. Moreover, patients typically relapse after discontinuation of therapy [23]. PVNS progression is usually insidious and begins many years before the onset of clinical symptoms. It will cause a substantial economic and psychological burden to the patient.

The GSE3698 dataset was downloaded from the GEO database, and 139 DEGs were identified. GO analysis revealed that these genes were mainly involved in antigen processing and presentation of exogenous peptide antigens via MHC class II and the immune response. KEGG analysis revealed that these genes were mainly involved in rheumatoid arthritis, asthma, and antigen processing and presentation. The above results suggested that inflammation and immunological reactions might play an essential role in PVNS. Zhao et al. confirmed that increased immune cell infiltration and cytokine secretion in PVNS synovial tissue affected its pathological process [10]. Cao et al. [24] found that inflammatory factors were significantly upregulated in PVNS knee synovial fluid. Therefore, our conclusions are consistent with the above reports.

Cytoscape is an open-source software project for integrating biomolecular interaction networks with high-throughput expression data and other molecular states into a unified conceptual framework [20]. The cytoHubba plugin ranks nodes with various algorithms based on network characteristics [21]. In this study, based on three of these algorithms, we identified TYROBP, CD74, PLEK, HLA-DPA1, HLA-DQB2, LAPTM5, HLA-DPB1, and FCER1G as the hub genes of PVNS. ROC analysis is a tool for evaluating model accuracy and has been commonly used for disease screening, diagnosis, treatment, and prognosis [22]. We further predicted the diagnostic value of 8 hub genes using ROC curves. The results showed that TYROBP, FCER1G, LAPTM5, and HLA-DPB1 could be used as high-value diagnostic markers for PVNS.

TYROBP is a type I transmembrane protein, also known as DAP12, consisting of a leader peptide, cysteine residues, a transmembrane fragment of aspartic acid residues, and a cytoplasmic tyrosine activation motif containing an immunoreceptor domain [25]. TYROBP is expressed in immune cells such as macrophages, monocytes, and osteoclasts [26]. Studies have shown that upregulated TYROBP might promote osteoporosis through osteoclasts [27]. Given that osteoclasts are essential in the pathological development of PVNS [28], we believe that TYROBP upregulation might be involved in the pathological progression of PVNS by regulating osteoclast differentiation. FCER1G is located on chromosome 1q23 and the gamma subunit of immunoglobulin E (IgE) encoding the crystallizable fragment (Fc) region [29]. Early studies found that various immune receptor activation signals transduced by FCER1G were particularly important in chronic responses such as autoimmunity and chronic infection [30]. FCER1G recognizes and eliminates non-self-antigens in normal immune system states but may trigger destructive inflammation, immune cell activation, phagocytosis, oxidative burst, and cytokine release in pathological conditions [31]. Furthermore, Dang et al. [32] discovered that FCER1G might be an ideal candidate for a robust, universal, and optimal marker of both macrophages and immune system players. Considering the macrophage accumulation in the synovial tissue PVNS, we believe that high expression of FCER1G might be involved in PVNS inflammation and immune pathways through macrophages. LAPTM5, also known as E3 protein, is a lysosomal membrane protein preferentially expressed in immune cells that interacts with the Nedd4 family of ubiquitin ligases, which play a role in hematopoiesis by preventing lymphocyte hyperactivation [33]. LAPTM5 is a positive regulator of the inflammatory signaling cascade in macrophages, and low expression of LAPTM5 inhibits the bladder cancer cell cycle [34, 35]. The above studies indicated that high expression of LAPTM5 might participate in the cascade regulation of the inflammatory signaling pathway, regulate the cell cycle, and stimulate cell hyperproliferation in PVNS. HLA-DPB1, located in the HLA class II region, is present on the cell surface of antigen-presenting cells. Genetic variants in HLA-DPB1 are associated with various autoimmune diseases, including rheumatoid arthritis, Graves’ disease, and multiple sclerosis [36]. However, whether HLA-DPB1 is associated with PVNS remains unclear.

To further explore the pathological process of immune cell infiltration in PVNS, we performed a comprehensive assessment using CIBERSORT. It is noteworthy that we used OA synovium as the control group. A certain degree of immune infiltration changes was also observed in the OA synovium. Some studies reported that, compared with normal synovium, there was no significant difference in M2 macrophages and neutrophils in OA synovium [37, 38]. Our results suggested that the infiltration of M2 macrophages and neutrophils was significantly increased in PVNS synovial tissue compared to OA synovial tissue. Our results further suggested that these two types of immune cell infiltration might play an essential role in the pathogenesis of PVNS. Previous studies have shown abnormal proliferation of macrophages and high expression of macrophage markers in the pathological process of PVNS [11]. Another study has shown that synovial tissue in PVNS abnormally secretes colony-stimulating factor 1, which recruits macrophages through binding to colony-stimulating factor 1 receptor to abnormally accumulate and locally form tumor-like masses, which constitute most of the tumor cell population and present as a tumor landscape effect in PVNS synovium [12]. In addition, the recruited macrophages could be polarized to the M2 phenotype [39, 40]. It has been reported that tumor-associated macrophages display an M2 macrophage phenotype, which can promote tumor growth and angiogenesis and invade local tissues. Therefore, we speculate that the highly infiltrated M2 macrophages in the PVNS synovium might be involved in mass growth and local bone erosion. Neutrophils can induce RANKL (receptor activator of NF-κB ligand) expression and induce osteoclast formation [41, 42]. Several studies have shown that a large number of osteoclast-type giant cells in PVNS synovial tissue could affect local bone remodeling, further leading to osteolysis and cortical bone destruction [10, 28]. Therefore, we believe that the highly infiltrated neutrophils in the PVNS synovium might lead to local osteolysis and cortical bone destruction through osteoclasts. Additionally, we revealed details of the correlation of immune cell infiltration and found that regulatory T cells had the greatest positive correlation with resting mast cells. CD8+ T cells had a positive correlation with activated CD4+ memory T cells. Resting dendritic cells had a positive correlation with activated mast cells. M1 macrophages had the greatest negative correlation with activated dendritic cells. Gamma-delta T cells had a negative correlation with CD8+ T cells and follicular helper T cells. Plasma cells had a negative correlation with neutrophils, M2 macrophages, and CD8+ T cells. Naive B cells had a negative correlation with memory B cells.

The correlation analysis displayed that FCER1G, LAPTM5, TYROBP, and HLA-DPB1 were positively correlated with M2 macrophages and neutrophils. FCER1G, LAPTM5, and TYROBP were negatively correlated with plasma cells. Combined with the above results, these results suggested that the high expression of FCER1G, LAPTM5, TYROBP, and HLA-DPB1 might promote the recruitment, differentiation, and proliferation of macrophages and neutrophils and further participate in the occurrence and development of PVNS. The effects of plasma cells on PVNS have not been studied thus far. It has been reported that plasma cells can produce a large number of cytokines and antibodies during tumor infiltration [43]. These antibodies could further promote antitumor immunity by driving antibody-dependent cellular cytotoxicity (ADCC), phagocytosis, and complement activation and enhancing antigen presentation by dendritic cells [44]. In addition, in this study, we found that plasma cells were significantly negatively correlated with neutrophils and M2 macrophages through correlation analysis. Therefore, we inferred that high expression of FCER1G, LAPTM5, and TYROBP might weaken the antitumor protective effect of plasma cells by reducing plasma cell infiltration. Our inferences require further experimental studies to demonstrate the interaction between the hub genes and immune cells infiltrated in PVNS.

However, our research has certain limitations. First, the samples in GSE3698 are unpaired, so more paired sample data are needed to validate our conclusions. Second, although the CIBERSORT analysis showed lower estimation bias than other methods, it was based on limited genetic data, which might still bias the final results. Third, the function and immune cell infiltration of the four biomarkers in PVNS were inferred by bioinformatics analysis. An experimental study with a larger sample size should be performed to validate our conclusions.

5. Conclusions

In this study, we found for the first time that LAPTM5, FCER1G, TYROBP, and HLA-DPB1 might be novel diagnostic markers for PVNS. In addition, for the first time, the infiltration of immune cells in the pathological process of PVNS was identified by CIBERSORT, and we found that M2 macrophages and neutrophils might be involved in the immune regulation process of PVNS. LAPTM5, FCER1G, TYROBP, and HLA-DPB1 were positively correlated with M2 macrophage and neutrophil infiltration, and LAPTM5, FCER1G, and TYROBP were negatively associated with plasma cell infiltration. These abnormally expressed genes and immune cells might play vital roles in the pathological process of PVNS. Further exploration of the regulatory interaction between these genes and immune cells will provide new clinical diagnostic markers and therapeutic targets for PVNS.

Data Availability

The data used to support the findings of this study are available from the GEO database (https://www.ncbi.nlm.nih.gov/geo/).

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Authors’ Contributions

Biao Chen and Bin Li designed the study. Boming Zhao and Yongjian Qi collected the data. Jun Zhang performed the data analysis and wrote the original draft. Bin Li and Biao Chen provided the writing guidance. Jun Chen helped revise the manuscript. All authors approved of the final version to be published. Jun Zhang and Bin Li contributed equally to this study.

Acknowledgments

This work was supported by grants from the National Natural Science Foundation of China (No. 81871779), the Zhongnan Hospital of Wuhan University Science, Technology and Innovation Seed Fund (No. znpy2019043), and the Program of Excellent Doctoral (Postdoctoral) of Zhongnan Hospital of Wuhan University (No. ZNYB2019001).

Supplementary Materials

See Tables S1-S2 in Supplementary Materials for all pathway enrichment analyses. (Supplementary Materials)

References

F. G. Verspoor, I. C. van der Geest, E. Vegt, R. P. Veth, W. T. van der Graaf, and H. W. Schreuder, “Pigmented villonodular synovitis: current concepts about diagnosis and management,” Future Oncology, vol. 9, no. 10, pp. 1515–1531, 2013.
View at: Publisher Site | Google Scholar
M. Brahmi, A. Vinceneux, and P. A. Cassier, “Current systemic treatment options for tenosynovial giant cell tumor/pigmented villonodular synovitis: targeting the CSF1/CSF1R axis,” Current Treatment Options in Oncology, vol. 17, no. 2, p. 10, 2016.
View at: Publisher Site | Google Scholar
L. van der Heijden, S. R. Piner, and M. A. van de Sande, “Pigmented villonodular synovitis: a crowdsourcing study of two hundred and seventy two patients,” International Orthopaedics, vol. 40, no. 12, pp. 2459–2468, 2016.
View at: Publisher Site | Google Scholar
G.-p. Xie, N. Jiang, C.-x. Liang et al., “Pigmented villonodular synovitis: a retrospective multicenter study of 237 cases,” PLoS One, vol. 10, no. 3, article e0121451, 2015.
View at: Publisher Site | Google Scholar
J. C. Aurégan, S. Klouche, Y. Bohu, N. Lefèvre, S. Herman, and Hardy, “Treatment of pigmented villonodular synovitis of the knee,” Arthroscopy, vol. 30, no. 10, pp. 1327–1341, 2014.
View at: Publisher Site | Google Scholar
J. S. Chang, J. P. Higgins, J. D. Kosy, and J. Theodoropoulos, “Systematic arthroscopic treatment of diffuse pigmented villonodular synovitis in the knee,” Arthroscopy Techniques, vol. 6, no. 5, pp. e1547–e1551, 2017.
View at: Publisher Site | Google Scholar
Y. Duan, J. Qian, K. Chen, and Z. Zhang, “Necessity of adjuvant postoperative radiotherapy for diffuse pigmented villonodular synovitis of the knee: a case report and literature review,” Medicine (Baltimore), vol. 97, no. 3, article e9637, 2018.
View at: Publisher Site | Google Scholar
F. Peyraud, S. Cousin, and A. Italiano, “CSF-1R inhibitor development: current clinical status,” Current Oncology Reports, vol. 19, no. 11, p. 70, 2017.
View at: Publisher Site | Google Scholar
B. P. Rubin, “Tenosynovial giant cell tumor and pigmented villonodular synovitis: a proposal for unification of these clinically distinct but histologically and genetically identical lesions,” Skeletal Radiology, vol. 36, no. 4, pp. 267-268, 2007.
View at: Publisher Site | Google Scholar
Y. Zhao, J. Lv, H. Zhang, J. Xie, H. Dai, and X. Zhang, “Gene expression profiles analyzed using integrating RNA sequencing, and microarray reveals increased inflammatory response, proliferation, and osteoclastogenesis in pigmented villonodular synovitis,” Frontiers in Immunology, vol. 12, p. 665442, 2021.
View at: Publisher Site | Google Scholar
I. Berger, H. Weckauf, B. Helmchen et al., “Rheumatoid arthritis and pigmented villonodular synovitis: comparative analysis of cell polyploidy, cell cycle phases and expression of macrophage and fibroblast markers in proliferating synovial cells,” Histopathology, vol. 46, no. 5, pp. 490–497, 2005.
View at: Publisher Site | Google Scholar
R. B. West, B. P. Rubin, M. A. Miller et al., “A landscape effect in tenosynovial giant-cell tumor from activation of CSF1 expression by a translocation in a minority of tumor cells,” Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 3, pp. 690–695, 2006.
View at: Publisher Site | Google Scholar
S. Sui, X. An, C. Xu et al., “An immune cell infiltration-based immune score model predicts prognosis and chemotherapy effects in breast cancer,” Theranostics, vol. 10, no. 26, pp. 11938–11949, 2020.
View at: Publisher Site | Google Scholar
A. M. Newman, C. L. Liu, M. R. Green et al., “Robust enumeration of cell subsets from tissue expression profiles,” Nature Methods, vol. 12, no. 5, pp. 453–457, 2015.
View at: Publisher Site | Google Scholar
J. I. Kawada, S. Takeuchi, H. Imai et al., “Immune cell infiltration landscapes in pediatric acute myocarditis analyzed by CIBERSORT,” Journal of Cardiology, vol. 77, no. 2, pp. 174–178, 2021.
View at: Publisher Site | Google Scholar
J. Gao, L. Shi, J. Gu et al., “Difference of immune cell infiltration between stable and unstable carotid artery atherosclerosis,” Journal of Cellular and Molecular Medicine, vol. 25, no. 23, pp. 10973–10979, 2021.
View at: Publisher Site | Google Scholar
I. Diboun, L. Wernisch, C. A. Orengo, and M. Koltzenburg, “Microarray analysis after RNA amplification can detect pronounced differences in gene expression using limma,” BMC Genomics, vol. 7, no. 1, p. 252, 2006.
View at: Publisher Site | Google Scholar
D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources,” Nature Protocols, vol. 4, no. 1, pp. 44–57, 2009.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. L. Gable, D. Lyon et al., “STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets,” Nucleic Acids Research, vol. 47, no. D1, pp. D607–d613, 2019.
View at: Publisher Site | Google Scholar
P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar
C. H. Chin, S. H. Chen, H. H. Wu, C. W. Ho, M. T. Ko, and C. Y. Lin, “cyto Hubba: identifying hub objects and sub-networks from complex interactome,” BMC systems biology, vol. 8, Suppl 4, 2014.
View at: Publisher Site | Google Scholar
N. A. Obuchowski and J. A. Bullen, “Receiver operating characteristic (ROC) curves: review of methods with applications in diagnostic medicine,” Physics in medicine and biology, vol. 63, no. 7, 2018.
View at: Publisher Site | Google Scholar
F. G. M. Verspoor, A. A. G. Zee, G. Hannink, I. C. M. van der Geest, R. P. H. Veth, and H. W. B. Schreuder, “Long-term follow-up results of primary and recurrent pigmented villonodular synovitis,” Rheumatology, vol. 53, no. 11, 2014.
View at: Publisher Site | Google Scholar
C. Cao, F. Wu, X. Niu et al., “Cadherin-11 cooperates with inflammatory factors to promote the migration and invasion of fibroblast-like synoviocytes in pigmented villonodular synovitis,” Theranostics, vol. 10, no. 23, pp. 10573–10588, 2020.
View at: Publisher Site | Google Scholar
E. Tomasello, L. Olcese, F. Vély et al., “Gene structure, expression pattern, and biological activity of mouse killer cell activating receptor-associated protein (KARAP)/DAP-12,” The Journal of Biological Chemistry, vol. 273, no. 51, pp. 34115–34119, 1998.
View at: Publisher Site | Google Scholar
J. A. Hamerman, M. Ni, J. R. Killebrew, C. L. Chu, and C. A. Lowell, “The expanding roles of ITAM adapters FcRgamma and DAP12 in myeloid cells,” Immunological Reviews, vol. 232, no. 1, pp. 42–58, 2009.
View at: Publisher Site | Google Scholar
B. Xia, Y. Li, J. Zhou, B. Tian, and L. Feng, “Identification of potential pathogenic genes associated with osteoporosis,” Bone & Joint Research, vol. 6, no. 12, pp. 640–648, 2017.
View at: Publisher Site | Google Scholar
R. Taylor, T. G. Kashima, H. Knowles, C. L. M. H. Gibbons, D. Whitwell, and N. A. Athanasou, “Osteoclast formation and function in pigmented villonodular synovitis,” The Journal of Pathology, vol. 225, no. 1, pp. 151–156, 2011.
View at: Publisher Site | Google Scholar
N. Novak, C. Tepel, S. Koch, K. Brix, T. Bieber, and S. Kraft, “Evidence for a differential expression of the FcepsilonRIgamma chain in dendritic cells of atopic and nonatopic donors,” The Journal of Clinical Investigation, vol. 111, no. 7, pp. 1047–1056, 2003.
View at: Publisher Site | Google Scholar
R. A. Sweet, K. M. Nickerson, J. L. Cullen, Y. Wang, and M. J. Shlomchik, “B cell-extrinsic Myd88 and Fcer1g negatively regulate autoreactive and normal B cell immune responses,” Journal of immunology, vol. 199, no. 3, pp. 885–893, 2017.
View at: Publisher Site | Google Scholar
A. M. Brandsma, P. M. Hogarth, F. Nimmerjahn, and J. H. Leusen, “Clarifying the confusion between cytokine and Fc receptor “common gamma chain”,” Immunity, vol. 45, no. 2, pp. 225-226, 2016.
View at: Publisher Site | Google Scholar
D. Dang, S. Taheri, S. Das, P. Ghosh, L. S. Prince, and D. Sahoo, “Computational approach to identifying universal macrophage biomarkers,” Frontiers in Physiology, vol. 11, p. 275, 2020.
View at: Publisher Site | Google Scholar
X. Cai, Y. Qiao, Y. Chen et al., “Expression and polymorphisms of lysosome-associated protein transmembrane 5 (LAPTM5) in patients with systemic lupus erythematosus in a Chinese population,” Biochemical Genetics, vol. 53, no. 7-8, pp. 200–210, 2015.
View at: Publisher Site | Google Scholar
L. Chen, G. Wang, Y. Luo et al., “Downregulation of LAPTM5 suppresses cell proliferation and viability inducing cell cycle arrest at G0/G1 phase of bladder cancer cells,” International Journal of Oncology, vol. 50, no. 1, pp. 263–271, 2017.
View at: Publisher Site | Google Scholar
W. K. Glowacka, P. Alberts, R. Ouchida, J. Y. Wang, and D. Rotin, “LAPTM5 protein is a positive regulator of proinflammatory signaling pathways in macrophages,” The Journal of Biological Chemistry, vol. 287, no. 33, pp. 27691–27702, 2012.
View at: Publisher Site | Google Scholar
M. Burek Kamenaric, M. Maskalan, Z. Grubic et al., “HLA-DPB1 matching in unrelated hematopoietic stem cell transplantation program contributes to a higher incidence of disease relapse,” Human Immunology, vol. 78, no. 11-12, pp. 665–671, 2017.
View at: Publisher Site | Google Scholar
Y. J. Deng, E. H. Ren, W. H. Yuan, G. Z. Zhang, Z. L. Wu, and Q. Q. Xie, “GRB10 and E2F3 as diagnostic markers of osteoarthritis and their correlation with immune infiltration,” Diagnostics, vol. 10, no. 3, 2020.
View at: Publisher Site | Google Scholar
W. Cai, H. Li, Y. Zhang, and G. Han, “Identification of key biomarkers and immune infiltration in the synovial tissue of osteoarthritis by bioinformatics analysis,” PeerJ, vol. 8, article e8390, 2020.
View at: Publisher Site | Google Scholar
T. Chanmee, P. Ontong, K. Konno, and N. Itano, “Tumor-associated macrophages as major players in the tumor microenvironment,” Cancers, vol. 6, no. 3, pp. 1670–1690, 2014.
View at: Publisher Site | Google Scholar
S. M. Pyonteck, L. Akkari, A. J. Schuhmacher et al., “CSF-1R inhibition alters macrophage polarization and blocks glioma progression,” Nature Medicine, vol. 19, no. 10, pp. 1264–1272, 2013.
View at: Publisher Site | Google Scholar
X. Hu, Y. Sun, W. Xu, T. Lin, and H. Zeng, “Expression of RANKL by peripheral neutrophils and its association with bone mineral density in COPD,” Respirology, vol. 22, no. 1, pp. 126–132, 2017.
View at: Publisher Site | Google Scholar
A. R. Kim, J. H. Kim, Y. H. Choi et al., “The presence of neutrophils causes RANKL expression in periodontal tissue, giving rise to osteoclast formation,” Journal of Periodontal Research, vol. 55, no. 6, pp. 868–876, 2020.
View at: Publisher Site | Google Scholar
V. D. Dang, E. Hilgenberg, S. Ries, P. Shen, and S. Fillatreau, “From the regulatory functions of B cells to the identification of cytokine- producing plasma cell subsets,” Current Opinion in Immunology, vol. 28, pp. 77–83, 2014.
View at: Publisher Site | Google Scholar
A. E. Gilbert, P. Karagiannis, T. Dodev et al., “Monitoring the systemic human memory B cell compartment of melanoma patients for anti-tumor IgG antibodies,” PLoS One, vol. 6, no. 4, article e19330, 2011.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Jun Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

277

Downloads

425

Citations