the mRNA expression value of LNPEP in different cancer types
Firstly, we performed pan-cancer analysis derived from TCGA by adopting Timer2.0 and Gent2 to estimate the expression pattern of LNPEP. Compared with normal tissues, LNPEP was displayed to be upregulated in the liver, bile duct, whereas down-regulated was revealed in the breast, central nervous system (CNS), kidney, lung, prostate, rectum, skin according to Timer analysis (Fig. 1A). Moreover, we estimated LNPEP expression data profiled of 72 paired cancer compared to normal tissues by utilizing HG-U133 microarray (GPL570 platform) of Gent2 database (Fig. 1B). Thus, it can be observed from the results LNPEP was up-regulated in several cancer cases concerning adipose tissue, adrenal gland, cervix, endometrium, esophagus, kidney, oral, liver, pancreas, pharynx, small intestine, and vulva. Conversely, down-regulation of LNPEP was reported in the field of blood, colon, head and neck, lung, muscle, ovary, prostate, skin, teeth, and vagina (Fig. 1C and D). To evaluate the correlation between LNPEP gene expression and the patient prognosis in 33 tumors, we used gene expression profile data and single-factor regression analysis to draw forest plots. In addition, OS showed that LNPEP was notably correlated to the prognosis of KIRC (p < 0.001), OV (p = 0.006), and READ (p = 0.042) (Fig. 1E). DSS displayed that LNPEP was markedly correlated with the prognosis of KIRC (p < 0.001) and OV (p = 0.009) (Fig. 1F). PFI reflected that LNPEP was associated with the prognosis of KIRC (p < 0.001), PCPG (p = 0.024), and SKCM (p = 0.024) (Fig. 1G). The expression of LNPEP in different cancer cell lines and normal tissues was explored via the BioGPS database. Ten cancer cell lines with the highest LNPEP expression level are displayed in Figure S1A. In normal cells, the LNPEP expression level was higher in immune cells (Figure S1B). Detailed information is shown in Supplementary Tables S3 and S4. The above results suggest that LNPEP might participate in the process of immune regulation.
Expression of LNPEP and its association with Clinico-pathological Characteristics in Patients with OV
Next, the mRNA level of LNPEP was down-regulated in OV compared to normal ovary tissues which were tested using the data sourced from the available database. Based on the GEPIA with the criterion of |Log2FC|>2 and p < 0.05, LNPEP has significantly low expression in patients with OV (Fig. 2A). Similarly, we also found that LNPEP downregulation arose from the GE-mini database (Fig. 2B). Subsequently, we analyzed the diagnostic values of LNPEP in OSC with ROC curves. Based on the Xiantao Xueshu web tool (https://www.xiantao.love/products), we performed a plot of the area under the curve (AUC) of LNPEP which was estimated to be 0.859 (Fig. 2C). From the high AUC value derived, it was speculated that LNPEP could have a greater potential to be a diagnostic biomarker for patients with OV. As described above, the mRNA levels of LNPEP were downregulated in OV, so we proceeded to test the protein level of LNPEP in OV. The protein levels for LNPEP in OV and normal ovary tissue were subjected to western blot analysis. The result showed that LNPEP expression was upregulated in OV, as compared to the mRNA levels (Fig. 2D). In the HPA dataset, the IHC staining data revealed that LNPEP was mainly localized in the cytoplasmic/membranous tissues and had it showed higher expressions in OV tissues (Figs. 2E, S3A, and S3B). To validate the above findings and investigate the clinicopathological roles and distribution of LNPEP expression in OV, immunohistochemical analysis of the 60 paraffin-embedded OV tissue blocks was performed. Representative immunohistochemical staining of LNPEP in OSC is illustrated in Fig. 2F. In addition, LNPEP expression differed in different immune subtypes of OV. Meanwhile, a significant correlation between LNPEP expression and molecular subtypes existed in OV (Figure S2). To further evaluate the prognostic value of LNPEP expression, univariate and multivariate Cox regression analysis was performed based on the OS, DSS, and PFI of patients with OV based on TCGA. As shown in Tables 1, 2, and 3, LNPEP expression, age class, race, FIGO stage, primary therapy outcome, tumor residual, and tumor status were termed as independent prognostic factors for the determination of the OS, DSS, and PFI of patients with OV, respectively.
Table 1
Univariate and multivariate Cox proportional hazards analysis for overall survival in the OV cohort.
| Univariate analysis | Multivariate analysis |
Variables | HR (95%CI) | P-value | HR (95%CI) | P-value |
LNPEP:high vs low | 1.441 (1.111–1.870) | 0.006 | 1.312 (0.949–1.814) | 0.1 |
FIGO stage: I + II vs III + IV | 1.981(1.419–2.764) | < 0.001 | 2.076 (0.490–8.806) | 0.322 |
Primary therapy outcome: PD + SD vs PR + CR | 0.301 (0.204–0.444) | < 0.001 | 0.350 (0.226–0.542) | < 0.001 |
Age: <=60 vs > 60 | 1.355 (1.046–1.754) | 0.021 | 1.286 (0.931–1.776) | 0.128 |
Tumor residual:NRD vs RD | 2.313 (1.486–3.599) | < 0.001 | 1.239 (0.722–2.124) | 0.437 |
Tumor status: Tumor free vs With tumor | 9.576 (4.476–20.486) | < 0.001 | 18.787 (4.593–76.844) | < 0.001 |
Table 2
Univariate and multivariate Cox proportional hazards analysis of disease-specific survival in OV cohort.
| Univariate analysis | Multivariate analysis |
Variables | HR (95%CI) | P-value | HR (95%CI) | P-value |
LNPEP:high vs low | 1.456 (1.099–1.928) | 0.009 | 1.300 (0.941–1.794) | 0.112 |
Race:Asian Black + African American vs White | 0.592 (0.370–0.946) | 0.028 | 0.774 (0.434–1.380) | 0.385 |
Primary therapy outcome: PD + SD vs PR + CR | 0.294 (0.198–0.436) | < 0.001 | 0.310 (0.200-0.479) | < 0.001 |
Tumor residual:NRD vs RD | 2.572 (1.580–4.187) | < 0.001 | 2.192 (1.280–3.754) | 0.004 |
Table 3
Univariate and multivariate Cox proportional hazards analysis of progression-free interval in OV cohort.
| Univariate analysis | Multivariate analysis |
Variables | HR (95%CI) | P-value | HR (95%CI) | P-value |
LNPEP:high vs low | 1.219 (0.962–1.544) | 0.101 | | |
FIGO stage:I + II vs III + IV | 1.573 (0.918–2.694) | 0.099 | | |
Primary therapy outcome: PD + SD vs PR + CR | 0.457 (0.325–0.642) | < 0.001 | 0.752 (0.521–1.086) | 0.128 |
Tumor residual:NRD vs RD | 1.695 (1.219–2.358) | 0.002 | 0.885 (0.594–1.320) | 0.550 |
Tumor status: Tumor free vs With tumor | 10.045(5.758–17.526) | < 0.001 | 12.988 (6.395–26.377) | < 0.001 |
Genetic Alterations of LNPEP in OV
Genomic alterations induce changes in gene expression. We explored the mutation frequency of LNPEP using the cBioPortal website. Five datasets (MSK, MSKCC, TCGA Pan-Cancer Atlas, TCGA Firehose legacy, and TCGA nature), which included 1,737 samples, were selected for analysis. The somatic mutation frequency of LNPEP in OV was observed to be 3%. It mainly consisted of missense mutations and synonymous substitution (Fig. 3A). We observed that OSC patients possessed a high frequency of gene alterations (Figure S4). In addition, the mutation types of LNPEP were further evaluated in the COSMIC database. Two pie charts were used to show the distribution of mutation types clearly (Figs. 3B and 3C). Missense substitutions occurred in approximately 34.47% of the samples, while synonymous substitutions occurred in 11.69% of the samples. Lastly, nonsense substitution occurred in 4.4% of the samples (Fig. 3B). The substitution mutations mainly occurred at C > T (27.08%), followed by G > A (16%), G > T (11.38%), and A > G (10.15%) (Fig. 3C). In addition, altered LNPEP had a significant correlation with the OS of OV (p < 0.05) (Fig. 3D).
The Prognostic and Diagnostic Values of LNPEP in Patients with OV
Survival time was determined from the date of primary tumor surgery to the time of death or last follow-up. To assess the effect of LNPEP expression on the OS, DSS, and PFI of patients with OV, and to explore the prognostic value of LNPEP in OV, survival analysis was performed based on information from the TCGA database. Results displayed that the upregulation of LNPEP expression was associated with poor OS (HR = 1.44 (1.11–1.87), log-rank p = 0.006), DSS (HR = 1.46 (1.10–1.93), log-rank p = 0.009), PFI (HR = 1.36 (1.07–1.73), log-rank p = 0.011) of patients with OV (Fig. 4A–C). Furthermore, high LNPEP expression was also associated with poor OS/DSS/PFI of patients with OV who were white, high grade histologic, anatomic neoplasm subdivision (bilateral), age, tumor residual (NRD and RD), and tumor status (with tumor and tumor-free) (Figures S5A–I, S6A–I, and S7A–I), but not those patients who were under the anatomic neoplasm subdivision (unilateral) (Figures S5J, S6J, and S7J). In conclusion, high LNPEP expression was associated with a poor prognosis of patients with OV.
Relationship between LNPEP Expression and Immune Cell Infiltration in Pan-Cancer Patients
Tumor-infiltrating lymphocytes are self-governed predictors of tumor sentinel lymph node status and survival rate response to therapies. We investigated the relationship between LNPEP expression and immune infiltration using TIMER2.0. The correlation coefficients between LNPEP expression and the abundances of seven immune infiltrates (CD8 + T cells, CD4 + T cells, B cells, neutrophils, macrophages, monocytes, and dendritic cells) were analyzed using Spearman tests (tumor purity adjusted). To explore the pan-cancer correlation between LNPEP expression and immune infiltration, we first evaluated the abundance of immune cell infiltration. As shown in Fig. 5A–G, the profiles of LNPEP associated with various immune infiltration displayed that it was positively correlated with immune cell infiltration levels of neutrophils. However, the data indicated that it was contradicted in CD8 + T cells, CD4 + T cells, B cells, macrophages, monocyte, and DCs. Correlation analyses using data from published work, which evaluated the 24 immune cells, showed that LNPEP was positively correlated with the infiltration levels of mast cells, T cells, T helper cells, Tcm, Th2 cells, whereas it was negatively correlated with those of NK CD56 bright cells (Fig. 5H).
LNPEP is Correlated with Immune Infiltration in OV
Afterward, we then performed a correlation analysis between LNPEP and tumor-infiltrating immune cells to evaluate the immunotherapy effect. Tumor purity acted as a critical factor, which effected the analysis of immune infiltration based on the genomic approach. We found LNPEP expression had a slightly negative correlation with tumor purity (r=-0.149, p = 1.81E-02). In addition, LNPEP expression had significant positive correlations with five immune infiltrates, including CD8 + T cells (r = 0.274, p = 1.17E-05), B cells (r = 0.161, p = 1.09E-02), macrophage (r = 0.274, p = 1.16E-05), monocytes (r = 0.227, p = 3.13E-04), dendritic cells (r = 0.128, p = 4.39E-02) (Fig. 6A). In addition, we also comprehensively clarified the correlations between LNPEP expression and related immune cell gene markers. Correlation coefficients were adjusted for tumor purity. The genetic markers of immune cells were used to analyze and identify the immune cells, including CD8 + T cells, general T cells, B cells, monocytes, M1 and M2 macrophages, neutrophils, dendritic cells, Th1, Th2, Tfh, Th17, Treg and T cell exhaustion (Table 4). Interestingly, LNPEP expression was critically correlated with the expression of markers of specific immune cells, such as CD8 + T cells marker, FCGR3B (r = 0.208, p < 0.001), SIGLEC5 (r = 0.248, p < 0.001), FPR1 (r = 0.146, p < 0.05); T cells marker, CD2 (r = 0.133, p < 0.05); B cells marker, CD19 (r = 0.145, p < 0.05); monocytes marker, CD115 (r = 0.162, p < 0.05) and CD86 (r = 0.193, p < 0.01); TAM marker, CD68 (r = 0.23, p < 0.001), IL10 (r = 0.232, p < 0.001); M1 macrophages marker, PTGS2 (r = 0.175, p < 0.01), CD163 (r = 0.267, p < 0.001), IRF5 (r = 0.146, p < 0.05); M2 Macrophages marker, VSIG4 (r = 0.163, p < 0.05), MS4A4A (r = 0.17, p < 0.01); neutrophils marker, CD11b (r = 0.198, p < 0.01), CCR7 (r = 0.185, p < 0.01); dendritic cells marker, ITGAM (r = 0.198, p < 0.01), and NRP1 (r = 0.189, p < 0.01). Moreover, the expression of LNPEP was associated with the expression of markers of specific subsets of T cells in OV, which included Th1 marker, TBX21 (r = 0.181, p < 0.01), STAT4 (r = 0.177, p < 0.01), STAT1 (r = 0.269, p < 0.001), tumor necrotic factor (TNF) (r = 0.135, p < 0.05); Th2 marker, STAT5A (r = 0.215, p < 0.001), QRSL1 (r = 0.23, p < 0.05), STAT6 (r = 0.161, p < 0.01); Tfh marker, BCL6 (r = 0.158, p < 0.05); Th17 marker, STAT3 (r = 0.395, p < 0.001); Treg marker, FOXP3 (r = 0.286, p < 0.001), STAT5B (r = 0.344, p < 0.001), and CCR8 (r = 0.195, p < 0.01); T cell exhaustion marker, CTLA4 (r = 0.152, p < 0.05), HAVCR2 (r = 0.212, p < 0.001), PDCD1 (r = 0.121, p < 0.05), LAG3 (r = 0.126, p < 0.05), GAMB (r = 0.166, p < 0.05). However, LNPEP expression did not reveal any significant correlation with CD4 + T cells. LNPEP was significantly correlated with immune stimulators, such as CD80 (rho = 0.217, p = 0.000136), IL2RA (rho = 0.214, p = 0.000163), MICB (rho = 0.207, p = 0.000265), and PVR (rho = 0.205, p = 0.00031) (Fig. 6B). The expression of LNPEP was also associated with immune inhibitors, including CD274 (rho = 0.148, p < 3.3E-05), CTLA (rho = 0.148, p = 0.00933), IDO1 (rho = 0.126, p = 0.0269), TGFB1 (rho = 0.148, p = 0.00927) (Fig. 6C). LNPEP expression was significantly correlated with CCL4 (rho = 0.132, p = 0.0212), CCL11 (rho = 0.199, p = 0.00045), CXCL13 (rho = 0.141, p = 0.0132), and CXCL14 (rho = 0.209, p = 0.00023) (Fig. 6D). Meanwhile, LNPEP expression was significantly associated with chemokine receptors, including CCR1 (rho = 0.175, p = 0.00211), CCR4 (rho = 0.19, p = 0.000825), CCR8 (rho = 0.166, p = 0.00363), and CXCR6 (rho = 0.188, p = 0.000947) (Fig. 6E). In summary, comprehensive analysis indicated that LNPEP may function as an immunoregulatory factor in OV.
Table 4
Correlation analysis between LNPEP and related genes of immune cells in TIMER database.
Description | Gene markers | OV (n = 303) |
None | Purity |
rho | P | rho | P |
CD8 + T cell | CD8A | 0.161 | ** | 0.114 | 0.0731 |
| CD8B | 0.155 | ** | 0.094 | 0.14 |
| FCGR3B | 0.244 | *** | 0.208 | *** |
| SIGLEC5 | 0.245 | *** | 0.248 | *** |
| FPR1 | 0.166 | ** | 0.146 | * |
T cell (general) | CD3D | 0.125 | * | 0.067 | 0.290 |
| CD3E | 0.172 | ** | 0.119 | 0.0617 |
| CD2 | 0.135 | * | 0.133 | * |
B cell | CD19 | 0.209 | *** | 0.145 | * |
Monocyte | CD14 | 0.147 | * | 0.104 | 0.101 |
| CD115(CSF1R) | 0.184 | ** | 0.162 | * |
| CD86 | 0.215 | *** | 0.193 | ** |
TAM | CD68 | 0.246 | *** | 0.23 | *** |
| IL10 | 0.248 | *** | 0.232 | *** |
M1 Macrophage | COX(PTGS2) | 0.206 | *** | 0.175 | ** |
| CD163 | 0.27 | *** | 0.267 | *** |
| IRF5 | 0.233 | *** | 0.146 | * |
M2 Macrophage | VSIG4 | 0.181 | ** | 0.163 | * |
| MS4A4A | 0.186 | ** | 0.17 | ** |
Neutrophils | CD11b(ITGAM) | 0.218 | *** | 0.198 | ** |
| CCR7 | 0.201 | *** | 0.185 | ** |
Dentritic cell | CD11c(ITGAX) | 0.28 | *** | 0.265 | *** |
| BDCA-4(NRP1) | 0.227 | *** | 0.189 | ** |
Th1 | T-bet (TBX21) | 0.182 | ** | 0.135 | * |
| STAT4 | 0.194 | *** | 0.177 | ** |
| STAT1 | 0.27 | *** | 0.269 | *** |
| TNF-a(TNF) | 0.188 | *** | 0.181 | ** |
Th2 | STAT5A | 0.256 | *** | 0.215 | *** |
| GATA3(QRSL1) | 0.243 | *** | 0.23 | * |
| STAT6 | 0.204 | *** | 0.191 | ** |
Tfh | BCL6 | 0.143 | * | 0.158 | * |
Th17 | STAT3 | 0.401 | *** | 0.395 | *** |
Treg | FOXP3 | 0.292 | *** | 0.286 | *** |
| CCR8 | 0.201 | *** | 0.195 | ** |
| STAT5B | 0.363 | *** | 0.344 | *** |
| TGFB(TGFB1) | 0.284 | *** | 0.297 | *** |
T cell exhaustion | CTLA4 | 0.185 | ** | 0.152 | * |
| TIM3 (HAVCR2) | 0.232 | *** | 0.212 | *** |
| PD-1(PDCD1) | 0.181 | ** | 0.121 | * |
| LAG3 | 0.126 | * | 0.126 | * |
| GZMB | 0.14 | * | 0.166 | * |
LNPEP Co-Expression Network in OV
The results of the co-expression pattern of LNPEP displayed that 1,914 genes were positively correlated with LNPEP, while 1,279 genes were negatively correlated with LNPEP (Fig. 7A). Heat maps showed the top 50 genes positively and negatively associated with LNPEP (Figs. 7B, C). KEGG pathway showed that co-expressed genes of LNPEP mainly participated in the ECM-receptor interaction, allograft rejection, Th1, and Th2 cell differentiation, Th17 cell differentiation, osteoclast differentiation, graft-versus-host disease, malaria, natural killer cell-mediated cytotoxicity, leishmaniasis, and TNF signaling pathway. (Fig. 7D). GO term annotation indicated enrichment in adaptive immune response, leukocyte cell-cell adhesion, T cell activation, cellular defense response, negative regulation of cell adhesion, cell adhesion mediated by integrin protein, lymphocyte-mediated immunity, leukocyte migration, regulation of cell-cell adhesion, tolerance induction, (Fig. 7E). In addition, we evaluated the pathway through which LNPEP may involve using GSEA in OV based on TCGA. The results indicated that LNPEP was significantly associated with cell adhesion and extracellular matrix (Figure S8A and S8B). These results suggested that the LNPEP expression network significantly affected the immune microenvironment and ligand-receptor interactions in OV.
Protein-Protein Interaction (PPI) Network Analysis for LNPEP
An aberrant expression of LNPEP may be involved in the process of diverse cancer types. Herein, we implemented two web resources, GeneMANIA, and STRING, to investigate the PPI network associated with LNPEP. It is well known that the protein alters a wide variety of biological processes and cellular signaling regarding protein interaction and signal transduction pathway[30, 31]. GeneMANIA, as an integrated network, focuses on functional prediction and performs an interaction network analysis based on a network-based gene ranking algorithm. Meanwhile, STRING is more inclined to the physical and functional interactions of the gene set. As Fig. 8A showed, GeneMANIA supplied the PPI-associated protein, including ORC6, CMTR1, KIF3B, LMNB2, PGD, NEK2, MTHFD2, PTTG1, EZH2, NCAPH, PKMYT1, RCHY1, SRPK1, RBM38, DMRTC2, ASXL1, SF1, PLAS2, PIAS4, and PIAS1. STRING analysis provided a predicted PPI network, which shares interaction with TNKS2, TNKS, STX4, VAMP2, RAB10, RAB14, RAB8A, TBC1D4, SLC2A4, RAB28 (Fig. 8B). The PPI network stats are the number of nodes: 11; the number of edges: 32; average node degree: 5.82; local clustering coefficient: 0.913 and PPI enrichment p < 0.05. The related parameters predicted PLAGL2 which is involved in the progression and prognosis of cancer.
T Cell Exhaustion Analysis
T cell exhaustion is a common feature of chronic infections and cancers in patients. Most patients with tumors possess a large number of exhausted T cells. Herein, we evaluated the relationship between LNPEP and marker genes of exhausted T cells, immune stimulator-related genes, immune inhibitor-related genes, chemokines, and receptors. The results determined that LNPEP was positively correlated with marker genes of exhausted T cells, immune stimulator-related genes, and immune inhibitor-related genes in OV (Figs. 9A and B). Additionally, LNPEP expression was also positively correlated with chemokines and chemokine receptors, such as CCL4 and CCL11, and their receptors CCR1 and CCR4 (Figs. 9C and D).