A new prognostic model for RHOV, ABCC2, and CYP4B1 to predict the prognosis and association with immune infiltration of lung adenocarcinoma
Original Article

A new prognostic model for RHOV, ABCC2, and CYP4B1 to predict the prognosis and association with immune infiltration of lung adenocarcinoma

Qiao Li1#, Xiao-Li Liu2#, Ni Jiang3#, Qian-Yun Li4, Yong-Xiang Song1, Xi-Xian Ke1, Hao Han1, Qian Luo1, Qiang Guo1,5, Xiang-Yu Luo5, Cheng Chen1

1Department of Thoracic Surgery, Affiliated Hospital of Zunyi Medical University, Zunyi, China; 2Department of Ultrasound, The People’s Hospital of Jianyang City, Jianyang, China; 3Department of Obstetrics and Gynecology, Women and Children’s Hospital of Chongqing Medical University, Chongqing, China; 4The Fourth Affiliated Hospital, Zhejiang University School of Medicine, Yiwu, China; 5Department of Cardiothoracic Surgery, Taihe Hospital, Hubei University of Medicine, Shiyan, China

Contributions: (I) Conception and design: C Chen, XY Luo; (II) Administrative support: C Chen, YX Song, XX Ke; (III) Provision of study materials or patients: QY Li, Q Guo; (IV) Collection and assembly of data: Q Li, H Han, Q Luo, Q Guo; (V) Data analysis and interpretation: Q Li, XL Liu, N Jiang, Q Guo; (VI) Manuscript writing: All authors; (VII) Final approval of manuscript: All authors.

#These authors contributed equally to this work.

Correspondence to: Cheng Chen. Department of Thoracic Surgery, Affiliated Hospital of Zunyi Medical University, 149 Dalian Road, Zunyi 563000, Guizhou, China. Email: 29217036@qq.com; Xiang-Yu Luo. Department of Cardiothoracic Surgery, Taihe Hospital, Hubei University of Medicine, 32 Renmin South Road, Shiyan 442012, Hubei, China. Email: 3282174937@qq.com.

Background: Lymph node metastasis is one of the important factors affecting the prognosis of lung adenocarcinoma (LUAD) patients. The key molecules in lymph node metastasis have not yet been fully revealed. Therefore, we aimed to construct a prognostic model based on lymph node metastasis-related genes to evaluate the prognosis of LUAD patients.

Methods: The differentially expressed genes (DEGs) in the process of LUAD metastasis were identified in The Cancer Genome Atlas (TCGA) database, and the biological roles of the DEGs were depicted using Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and a protein-protein interaction (PPI) network. Survival analysis and Cox regression analysis were used to identify the genes related to the prognosis of patients with LUAD, and a nomogram and a prognostic model were constructed. The potential prognostic value, immune escape, and regulatory mechanisms of the prognostic model in LUAD progression were explored through survival analysis and gene set enrichment analysis (GSEA).

Results: A total of 75 genes were upregulated, and 138 genes were downregulated in tissues of lymph node metastasis. The expression levels of STC1, CYP17A1, RHOV, GUCA2B, TM4SF20, DEFB1, CRHR2, ABCC2, CYP4B1, KRT16, and NTS were revealed as risk factors for a poor prognosis in LUAD patients. High-risk LUAD patients had a poor prognosis in the prognostic model based on RHOV, ABCC2, and CYP4B1. The clinical stage and the risk score were found to be independent risk factors for a poor prognosis in LUAD patients, and the risk score was associated with the tumor purity, T cell, natural killer (NK) cell, and other immune cells. The prognostic model might affect the progression of LUAD using DNA replication, the cell cycle, P53, and other signaling pathways.

Conclusions: Lymph node metastasis-related genes RHOV, ABCC2, and CYP4B1 are associated with a poor prognosis in LUAD. A prognostic model based on RHOV, ABCC2, and CYP4B1 might predict the prognosis of LUAD patients and be associated with immune infiltration.

Keywords: RHOV; ABCC2; CYP4B1; lymph node metastasis; lung adenocarcinoma (LUAD)


Submitted Feb 07, 2023. Accepted for publication Mar 17, 2023. Published online Apr 10, 2023.

doi: 10.21037/jtd-23-265


Highlight box

Key findings

• Our new prognostic model based on lymph node metastasis-related genes RHOV, ABCC2, and CYP4B1 might predict the prognosis of LUAD patients and be associated with immune infiltration.

What is known and what is new?

• Lymph node metastasis is related to the poor prognosis of patients with LUAD.

• The molecular mechanisms of lymph node metastasis have not been fully understood in LUAD. Lymph node metastasis-related genes RHOV, ABCC2, and CYP4B1 were found to be associated with the prognosis of LUAD patients.

What is the implication, and what should change now?

• A prognostic model based on RHOV, ABCC2, and CYP4B1 might predict the prognosis of cancer patients, which might become a tool to predict the prognosis of LUAD patients.


Introduction

Lung cancer is one of the most common malignant tumors worldwide, as well as one of the main causes of cancer-related death (1,2). In China, the morbidity and mortality of lung cancer remain high throughout the entire year (3). Lung adenocarcinoma (LUAD) is one of the common subtypes of lung cancer (4). Although there are many effective targeted drugs for treating cancer patients, LUAD remains one of the most common and fatal cancers worldwide.

Early metastasis of LUAD is one of the key factors leading to cancer progression to the middle and advanced stages (5-7). Lymph node metastasis is one of the common modes of LUAD metastasis and is one of the risk factors affecting the long-term survival of patients with LUAD (6,7). Luo et al. reported that the metastasis rates of N1 and N2 stages increased with the tumor diameter (6). Dai et al. reported that 15% of LUAD patients had lymph node metastasis. Compared with node-negative patients, recurrence-free survival (RFS) and overall survival (OS) have been shown to be significantly decreased in patients with metastasis (7). However, the molecules and the signaling mechanisms of lymph node metastasis in LUAD remain unelucidated. The Cancer Genome Atlas (TCGA) project aims to improve the ability to prevent, diagnose, and treat cancer using high-throughput genome analysis technology. Multiple cancer types and data of genes, microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and others are displayed in the TCGA database (8,9). Lymph node metastasis was a risk factor for poor prognosis in LUAD patients. However, the roles of lymph node metastasis-related genes in the progression of LUAD has not been fully revealed. Therefore, important molecular markers in the process of lymph node metastasis were explored based on the data from the TCGA database. The roles of the established lymph node metastasis-related the prognostic model were investigated in LUAD metastasis to improve the treatment value for LUAD patients. The relationship between the prognostic model and the prognosis, as well as immune cell infiltration of LUAD, were explored to understand the key molecules in LUAD metastasis. We present the following article in accordance with the TRIPOD reporting checklist (available at https://jtd.amegroups.com/article/view/10.21037/jtd-23-265/rc).


Methods

TCGA data download and visualization analysis

The gene expression data of 594 cases of LUAD high-throughput sequence-fragments per kilobase million (HTSeq-FPKM) and the clinical data of 522 cancer patients were downloaded from the TCGA database. The gene expression data of 535 LUAD tissues were included, and LUAD gene expression data with the N0-3 stage were extracted in our study. A total of 330 LUAD patients were lymph node-negative, whereas 171 LUAD patients were lymph node-positive. Differentially expressed genes (DEGs) of LUAD metastasis were screened using the limma package of R software (version 4.0.2; https://www.r-project.org/) with the criteria of a false discovery rate (FDR) <0.05 and |logFC| >1, and the visualization results were displayed using a volcano map. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).

Biological functions and protein-protein interaction (PPI) network of lymph node metastasis-related genes

Biological processes, cellular components, molecular functions, and signaling mechanisms involved in lymph node metastasis-related genes were investigated using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis in the Database for Annotation, Visualization, and Integrated Discovery (DAVID) with the criterion of FDR <0.05. A PPI network of lymph node metastasis-related genes was visualized in the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database and enriched for analysis using the MCODE method in Cytoscape software (10).

Kaplan-Meier (K-M) survival analysis

Grouped by the median value of DEGs, the relationship between high- and low-expressed DEGs and the OS of LUAD patients were explored via K-M survival analysis and with P<0.05 for filtering and screening. According to the median score, the patients were divided into high-risk and low-risk groups. The prognostic value of LUAD patients in the high-risk and low-risk groups was presented using K-M survival analysis.

Construction of nomogram and prognostic model of lymph node metastasis-related genes

The roles of DEGs in lymph node metastasis were identified in the prognosis of LUAD patients using univariate Cox regression analysis with the criterion of P<0.001. On this basis, multivariate Cox regression analysis and the Akaike information Criterion (AIC) method were performed to screen the factors affecting a poor prognosis of LUAD patients, and a nomogram and a prognostic model were constructed (11).

Construction of prognostic model-related nomogram

The risk score data and the clinicopathological characteristic data of LUAD patients were matched. The relationship between clinicopathological characteristics and the lymph node metastasis prognostic model, and the prognosis of LUAD patients were investigated using univariate and multivariate Cox regression analysis, and a nomogram was constructed based on the multivariate Cox analysis results.

Gene set enrichment analysis (GSEA)

The regulatory mechanisms in which genes might be involved were explored using GSEA (12,13). Grouping of the gene expression data of LUAD was performed according to the median prognostic model score, and the signaling pathways of the lymph node metastasis-related prognostic model were explored. The screening criterion of GSEA: NOM P<0.05.

Immune analysis of prognostic model

Tissue samples from patients with LUAD were scored using Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT), microenvironment cell populations (MCP)-counter, and estimation methods (14,15). Pearson correlation analysis was used to explore the relationship between the immune score and the levels of tumor purity, immune cells, and immune cell markers, and P<0.05 was considered a significant screening criterion.

Statistical analysis

Perl (https://www.perl.org/) and R were used for data processing and statistical analysis. Cox regression and K-M survival analyses were conducted to filter the risk factors of OS in patients with LUAD, and ROC analysis was carried out to assess the role of a gene-associated nomogram in lymph node metastasis. The relationship between the immune score and the levels of tumor purity, immune cells, and immune cell markers was explored using Pearson correlation analysis. The expression of DEGs in high- and low-risk groups was detected using a t-test, and P<0.05 was considered statistically significant.


Results

Identification of DEGs related to metastasis in LUAD

Comparing the tissues with lymph node-negative LUAD patients, there were 213 DEGs in the tissues of lymph node-positive LUAD patients (Table S1 and Table S2). Among them, 75 DEGs were upregulated (Table S1), and 138 DEGs were downregulated (Table S2). The top 15 DEGs in LUAD tissues were shown by fold changes (Figure 1). In detail, the expression levels of LRRC38, TAC1, CALB1, CGB5, KRT20, TRIM48, TM4SF20, NNAT, GCG, PI3, CYP2B6, SPAG11B, GUCA2B, RHCG, and MUC2 were increased in the tissues of lymph node-positive LUAD patients (Figure 1A), whereas the expression of ALB, NPY, AFP, SPINK4, WFDC12, DKK4, WFDC5, DLK1, MSTN, FABP7, HIST1H4C, HIST1H1B, GC, VTN, and FTHL17 was decreased in the tissues of lymph node-positive LUAD patients (Figure 1B).

Figure 1 Fifteen DEGs of lymph node metastasis in LUAD shown using heatmap and histogram. (A) Overexpressed genes; (B) lowly expressed genes. DEGs, differentially expressed genes; LUAD, lung adenocarcinoma; FC, fold change.

Biological functions and PPI network of lymph node metastasis-related genes

The DEGs of lymph node metastasis were found to be involved in the DNA replication-dependent nucleosome assembly, negative regulation of gene expression, cellular protein metabolic process, extracellular exosome, DNA-templated transcription and initiation, positive regulation of cytokine secretion, drug metabolic process, WNT signaling pathway, chemokine production, receptor binding, toll-like receptor 4 binding, and other functions (Figure 2A-2C and Table S3). The signaling pathways in which the DEGs of lymph node metastasis were involved were viral carcinogenesis, steroid hormone biosynthesis, drug metabolism-cytochrome P450, chemical carcinogenesis, the PPAR signaling pathway, and others (Figure 2D and Table 1). The PPI network was presented and was enriched for analysis using the MCODE method (Figure 3 and Figure S1).

Figure 2 Functions and mechanisms of LUAD metastasis-related DEGs using GO and KEGG analysis. BP, biological processes; GO, Gene Ontology; CC, cellular components; KEGG, Kyoto Encyclopedia of Genes and Genomes; LUAD, lung adenocarcinoma; DEGs, differentially expressed genes.

Table 1

Signaling pathways in which DEGs of lymph node metastasis were involved

Term Content Count P value
Hsa05322 Systemic lupus erythematosus 23 9.88E-20
Hsa05034 Alcoholism 24 2.99E-18
Hsa05203 Viral carcinogenesis 13 5.36E-06
Hsa00140 Steroid hormone biosynthesis 7 6.31E-05
Hsa00980 Metabolism of xenobiotics by cytochrome P450 7 2.48E-04
Hsa00830 Retinol metabolism 6 9.89E-04
Hsa04913 Ovarian steroidogenesis 5 0.00276348
Hsa00982 Drug metabolism-cytochrome P450 5 0.008937351
Hsa04975 Fat digestion and absorption 4 0.011324964
Hsa05204 Chemical carcinogenesis 5 0.015572222
Hsa03320 PPAR signaling pathway 4 0.046733944
Hsa05202 Transcriptional misregulation in cancer 6 0.050691025

DEGs, differentially expressed genes; PPAR, peroxisome proliferator activated receptor.

Figure 3 PPI network of LUAD metastasis-related DEGs. PPI, protein-protein interaction; LUAD, lung adenocarcinoma; DEGs, differentially expressed genes.

Metastasis genes related to prognosis of LUAD

The results of K-M survival analysis revealed that the expression levels of PI3, CALB1, STC1, STAR, HIST1H4B, CYP17A1, HIST2H2AB, RHOV, GUCA2B, TM4SF20, KRT20, HIST1H4A, PI15, GLP2R, KRT78, DEFB1, CRHR2, ABCC2, HIST1H2BO, CYP4B1, LIPF, S100G, CPB1, OTX2, KRT16, CYP2A6, and NTS were related to the OS of LUAD patients (Table 2).

Table 2

Screening of prognostic metastasis-related genes using K-M survival analysis

Gene P value
ABCC2 2.504e-02
CALB1 4.097e-02
CPB1 1.504e-02
CRHR2 2.322e-02
CYP2A6 7.923e-03
CYP4B1 2.655e-03
CYP17A1 9.490e-05
DEFB1 1.287e-02
GLP2R 2.748e-02
GUCA2B 2.236e-03
HIST1H2BO 2.663e-04
HIST1H4A 3.322e-02
HIST1H4B 4.416e-02
HIST2H2AB 1.243e-02
KRT16 6.119e-04
KRT20 4.358e-02
KRT78 1.235e-02
LIPF 3.390e-02
NTS 4.799e-02
OTX2 1.315e-02
PI15 2.613e-02
PI3 3.250e-02
RHOV 7.444e-03
STAR 4.710e-03
STC1 3.702e-03
TM4SF20 4.258e-02
S100G 1.731e-03

LUAD, lung adenocarcinoma; K-M, Kaplan-Meier.

Construction of the prognostic model

Univariate Cox regression analysis revealed that the expression levels of STC1, RHOV, GUCA2B, ABCC2, CYP4B1, KRT16, and NTS might be risk factors for the OS of LUAD patients (Figure S2). Multivariate Cox regression analysis and the AIC method revealed that RHOV, ABCC2, and CYP4B1 were independent risk factors for a poor prognosis in patients with LUAD, and a prognostic model based on the expression levels of RHOV, ABCC2, and CYP4B1 was constructed. In addition, a nomogram based on the expression levels of RHOV, ABCC2, and CYP4B1 was constructed (Figure 4).

Figure 4 The nomograms of prognostic genes on overall survival in LUAD. LUAD, lung adenocarcinoma.

Risk score of lymph node metastasis-related genes was associated with poor prognosis in LUAD patients

Figure 5A-5C depict the relationship between the risk score and the prognosis of cancer patients. K-M survival analysis illustrated that OS in the low-risk group was significantly higher than in the high-risk group in LUAD patients (Figure 5D). Univariate Cox regression analysis demonstrated that the clinical stage, T stage, lymph node metastasis, and risk score were the influencing factors of a poor prognosis in LUAD patients (Figure S3A). Multivariate Cox regression analysis depicted that the clinical stage and prognostic model score were independent risk factors for a poor prognosis in LUAD patients (Figure S3B). A nomogram based on the results of the Cox regression analysis was constructed to assess the prognosis of cancer patients (Figure 6).

Figure 5 Evaluation of survival time of patients with LUAD in prognostic model. (A) Prognostic model-related genes showed using heatmap; (B,C) relationship between risk score and prognosis of cancer patients; (D) Kaplan-Meier survival analysis showing OS of LUAD patients in low- and high-risk groups. LUAD, lung adenocarcinoma; OS, overall survival.
Figure 6 Nomogram related to prognostic model.

Prognostic model participation in the signaling mechanisms of LUAD

In the signaling mechanism module, DNA replication, cell cycle, homologous recombination, mismatch repair, proteasome, pyrimidine metabolism, base excision repair, pentose phosphate pathway, spliceosome, P53 signaling pathway, oocyte meiosis, ubiquitin-mediated proteolysis, basal transcription factors, and other pathways were significantly enriched in the high-risk group (Table 3).

Table 3

Signaling pathways enriched in high-risk score based on genes of lymph node metastasis

Name Size NES NOM P value
DNA replication 36 2.188415 0
Cell cycle 124 2.101134 0
Homologous recombination 28 2.0641458 0
Mismatch repair 23 2.0362506 0
Proteasome 44 2.0210252 0
Pyrimidine metabolism 97 2.0000293 0
Base excision repair 33 1.98535 0.002083333
Pentose phosphate pathway 27 1.9286441 0.004024145
Nucleotide excision repair 44 1.8649865 0.005725191
P53 signaling pathway 68 1.7573924 0.007905139
Riboflavin metabolism 15 1.6678965 0.011494253
Spliceosome 126 1.8515968 0.014084507
Pathogenic escherichia coli infection 55 1.7213273 0.015717093
Fructose and mannose metabolism 32 1.6948974 0.015904572
N glycan biosynthesis 46 1.7191821 0.018181818
Protein export 23 1.7690526 0.019646365
Oocyte meiosis 112 1.5964646 0.023715414
RNA degradation 57 1.628982 0.028957529
Glycolysis gluconeogenesis 61 1.6014928 0.029166667
Ubiquitin mediated proteolysis 133 1.5566719 0.037328094
Purine metabolism 156 1.4549121 0.042769857
Glyoxylate and dicarboxylate metabolism 16 1.6099515 0.042857144
Basal transcription factors 35 1.4993141 0.04828974
Galactose metabolism 25 1.5636381 0.05

NES, normalized enrichment score; NOM, nominal.

Prognostic model in LUAD immune microenvironment

Based on the results of estimation methods concerning LUAD tissues, the risk score was significantly correlated with tumor purity, stromal score, immune score, and estimated score expression levels using correlation analysis (Figure 7A-7D). In the high-risk group, the expression of tumor purity was increased, and the expression levels of the stromal score, immune score, and estimate score were decreased (Figure 7E-7H). Based on the results of MCP-counter methods concerning LUAD tissues, the risk score was significantly associated with T cells, endothelial cells, neutrophils, and other immune cells (Table 4). Based on the results of CIBERSORT methods concerning LUAD tissues, the risk score was significantly associated with the levels of B cell memory, T cell CD8, T cell follicular helper, and other immune cells (Figure 8 and Table 5). In the LUAD tissues, the risk score was significantly associated with the levels of immune cell markers (Table 6). More specifically, the risk score was significantly correlated with the expression levels of BCL6, CCR7, CCR8, and other markers of immune cells.

Figure 7 prognostic model in LUAD immune microenvironment. (A) Tumor purity; (B) immune score; (C) stromal score; (D) estimate score. *, P<0.05; **, P<0.01. LUAD, lung adenocarcinoma.

Table 4

Risk score was significantly associated with immune cells based on data of MCP-counter analysis

Immune cells Correlation coefficient P value
T cells −0.149 0.001
CD8 T cells 0.024 0.597
Cytotoxic lymphocytes 0.03 0.505
B lineage −0.126 0.006
NK cells 0.007 0.886
Monocytic lineage −0.068 0.137
Myeloid dendritic cells −0.290 <0.001
Neutrophils −0.109 0.017
Endothelial cells −0.321 <0.001
Fibroblasts 0.057 0.214

MCP, microenvironment cell populations; NK, natural killer.

Figure 8 prognostic model score is significantly associated with immune cell infiltration. (A) Plasma cells. (B) T cells follicular helper. (C) T cells CD8. (D) T cells CD4 memory resting. (E) T cells CD4 memory activated. (F) B cells memory. (G) Monocytes. (H) NK cells resting. (I) NK cells activated. (J) Macrophages M0. (K) Mast cells resting. (L) Mast cells activated. NK, natural killer.

Table 5

Significant associations of risk score with immune cells based on data of CIBERSORT analysis

Immune cells Correlation coefficient P value
B cells memory −0.123 0.007
B cells naive 0.055 0.227
Plasma cells 0.108 0.018
T cells CD8 0.097 0.034
T cells CD4 memory resting −0.278 <0.001
T cells CD4 memory activated 0.160 <0.001
T cells follicular helper 0.220 <0.001
T cells regulatory −0.021 0.652
T cells gamma delta 0.012 0.786
NK cells resting 0.113 0.013
NK cells activated 0.132 0.004
Monocytes −0.184 <0.001
Macrophages M0 0.220 <0.001
Macrophages M1 0.077 0.091
Macrophages M2 −0.164 <0.001
Dendritic cells resting −0.220 <0.001
Dendritic cells activated 0.093 0.043
Mast cells resting −0.252 <0.001
Mast cells activated 0.134 0.003
Eosinophils 0.019 0.672
Neutrophils 0.095 0.037

NK, natural killer.

Table 6

Significant associations of risk score with markers of immune cells

Cell markers Correlation coefficient P value
BCL6 −0.102 0.025
CCR8 −0.106 0.020
CD2 −0.100 0.028
CD3E −0.101 0.027
CD8B 0.024 0.597
CD79A −0.081 0.076
CEACAM8 −0.156 <0.001
FOXP3 −0.091 0.046
GZMB 0.170 <0.001
HLA-DPA1 −0.307 <0.001
HLA-DQB1 −0.256 <0.001
IFNG 0.094 0.040
IL17A −0.007 0.877
IRF5 −0.046 0.311
ITGAX −0.154 <0.001
MS4A4A −0.165 <0.001
NRP1 −0.167 <0.001
PTGS2 0.039 0.392
STAT3 −0.072 0.114
STAT5B −0.239 <0.001
TBX21 0.014 0.755
TNF −0.002 0.971
CCR7 −0.180 <0.001
CD1C −0.268 <0.001
CD3D −0.058 0.208
CD8A 0.022 0.626
CD19 −0.113 0.013
CD163 −0.092 0.044
CTLA4 −0.052 0.251
GATA3 0.185 <0.001
HAVCR2 −0.114 0.012
HLA-DPB1 −0.326 <0.001
HLA-DRA −0.273 <0.001
IL13 −0.100 0.029
IL21 −0.012 0.785
ITGAM −0.201 <0.001
LAG3 0.023 0.613
NOS2 −0.009 0.845
PDCD1 0.018 0.701
STAT1 0.185 <0.001
STAT5A −0.182 <0.001
STAT6 −0.180 <0.001
TGFB1 −0.112 0.014
VSIG4 −0.142 0.002

Discussion

Lung cancer is the most common malignant tumor worldwide and has the highest mortality rate. LUAD accounts for about 40% of lung cancer (1,3,4). Currently, the annual mortality rate of LUAD patients remains high. In the past few decades, the application and mining of big data have constituted one of the important means by which to diagnose, treat, and evaluate the prognosis of patients with cancer. The TCGA database contains high-quality tumor genome data and clinical information of patients. In our study, we found that PI3, CALB1, STC1, STAR, HIST1H4B, CYP17A1, HIST2H2AB, RHOV, GUCA2B, TM4SF20, KRT20, HIST1H4A, PI15, GLP2R, KRT78, DEFB1, CRHR2, ABCC2, HIST1H2BO, CYP4B1, LIPF, S100G, CPB1, OTX2, KRT16, CYP2A6, and NTS were unusually expressed, and related to the OS of LUAD using the TCGA database. RHOV, ABCC2, and CYP4B1 were independent risk factors for a poor prognosis in LUAD patients and were correlated with the prognostic model. Currently, RHOV, ABCC2, and CYP4B1 have important biological roles in cancer (16-21). For example, Shepelev et al. found that the expression of RHOV was increased in lung cancer tissues and cells, and the increased expression of RHOV was associated with a poor prognosis of patients (16). Chen et al. reported that ABCC2 was overexpressed in various human cancers (18). The expression of ABCC2 was upregulated in cisplatin-resistant A549 cells (A549/DDP). Interfering with the expression of ABCC2 could reverse the resistance of A549/DDP cells to cisplatin in vitro, promote G1 phase arrest, and activate the expression of PARP and caspase-3 proteins. The knockout of ABCC2 expression in vivo could enhance the cytotoxicity of cisplatin to subcutaneous transplanted tumors (18). The expression of CYP4B1 in LUAD decreased, which was related to the history of drug treatment, radiotherapy, and the survival status of cancer patients (21). In addition, we established prognostic model for RHOV, ABCC2, and CYP4B1 and demonstrated that RHOV and ABCC2 were overexpressed and CYP4B1 was underexpressed in the high-risk group, and the OS of LUAD patients in the low-risk group was significantly higher than that in the high-risk group. Univariate Cox regression analysis showed that clinical stage, T stage, lymph node metastasis, and risk score were influencing factors of a poor prognosis in patients with LUAD. Multivariate Cox regression analysis demonstrated that clinical stage and prognostic model score were independent risk factors for a poor prognosis in LUAD patients. This demonstrated that RHOV, ABCC2, and CYP4B1 played an important role in lung cancer and indicated that the prognostic model and the nomogram based on the genes RHOV, ABCC2, and CYP4B1 of LUAD metastasis have important predictive value.

The process of lung cancer metastasis involves a variety of biological processes and changes of molecular markers (22-27). The expression of long-chain noncoding RNA NSCLCAT1 was upregulated in the NSCLC tissues. NSCLCAT1 could increase the viability, migration, and invasion of NSCLC cells and reduce apoptosis by inhibiting the expression of CDH1 and mediating the hippo signaling pathway (22). In NSCLC cells, the inhibition of PLK1 expression could change the expression of genes related to DNA damage, replication, and repair (23). Rig-G was frequently downregulated in lung cancer tissues and cell lines and associated with a poor prognosis in lung cancer patients. The overexpression of Rig-G has been shown to result in a significant reduction in cell growth and migration inhibition in A549 and NCI-H1944 cells, along with a reduced epithelial-to-mesenchymal transition. Rig-G acted as a tumor suppressor through the p53 signaling mechanism (27). The metastasis genes and the prognostic model had important biological value in the cell cycle, DNA replication, p53 signaling pathway, and other mechanisms using GO, KEGG, and GSEA, whereby proving that our prognostic model had good predictive value in the progression of LUAD.

Recently, immunotherapy has become a dominant therapeutic theme. Immunotherapy could improve long-term survival and the chances of surgery in patients with LUAD (28,29). For example, the CDK5 inhibitor resulted in decreased PD-L1 protein expression in human lung adenocarcinoma (LLC) cells. PD-L1 protein degradation was mediated using the E3 ligase TRIM21 ubiquitination-proteasome pathway (29). In vitro, the deletion of CDK5 in the LLC of mice has not been shown to affect cell proliferation. However, the attenuation of CDK5 or binding to anti-PD-L1 was shown to greatly inhibit tumor growth in vivo mouse model of LLC implantation. CDK5 disruption caused higher levels of CD3, CD4, and CD8 T cells in the spleen and decreased PD-1 expression in CD4 and CD8 T cells, which provided a potential therapeutic target for LUAD combination immunotherapy (29). In our research, the prognostic model score based on the RHOV, ABCC2, and CYP4B1 was significantly correlated with tumor purity, stromal score, immune score, and estimate score expression levels. The risk score was significantly associated with T cells, endothelial cells, neutrophils, CD8, T cell follicular helper cells, and other immune cells. In LUAD tissues, the risk score was significantly associated with the levels of immune cell markers. In addition, study has confirmed that RHOV is related to the regulation of immune cells (30). Specifically, RHOV expression increased during the differentiation of macrophages into osteoclasts, while a large number of macrophages showed apoptosis. When osteoprotegerin (OPG) inhibits the differentiation of macrophages into osteoclasts, and then OPG can inhibit apoptosis, which is related to the down-regulation of RHOV expression level (30). However, the relationship between ABCC2, and CYP4B1 and immune cell regulation has not been reported in the literature, which will be our research direction in the future.

The molecular mechanisms of lymph node metastasis have not been fully understood in LUAD. Lymph node metastasis-related genes RHOV, ABCC2, and CYP4B1 were found to be associated with the prognosis of LUAD patients. In addition, this study used big data samples to provide new candidate biomarkers related to LUAD metastasis for the prognosis of LUAD patients, which has the advantage of high reliability. However, our study also had some limitations. First, the expression levels of RHOV, ABCC2, and CYP4B1 in clinical LUAD samples need to be verified. Moreover, the expression levels of RHOV, ABCC2, and CYP4B1 and the values of their constructed prognostic model in the prognosis of LUAD patients need to be explored. The roles and the underlying signaling mechanisms of our constructed prognostic model were explored in LUAD using basic research in the future.


Conclusions

There are many molecular DEGs in the process of LUAD metastasis. RHOV, ABCC2, and CYP4B1 are influencing factors of a poor prognosis, and clinical stage and risk score are independent risk factors for a poor prognosis in patients with LUAD. The prognostic model might be involved in the progression of LUAD through the cell cycle, DNA replication, p53 signaling pathway, and others. A prognostic model based on RHOV, ABCC2, and CYP4B1 might predict the prognosis of LUAD patients and be associated with immune infiltration.


Acknowledgments

Thanks to the TCGA database for providing the data used for analysis.

Funding: None.


Footnote

Reporting Checklist: The authors have completed the TRIPOD reporting checklist. Available at https://jtd.amegroups.com/article/view/10.21037/jtd-23-265/rc

Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://jtd.amegroups.com/article/view/10.21037/jtd-23-265/coif). The authors have no conflicts of interest to declare.

Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).

Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.


References

  1. Chen PW, Huang SK, Chou WC, et al. Severinia buxifolia-isolated acridones inhibit lung cancer invasion and decrease HIFα protein synthesis involving 5'UTR-mediated translation inhibition. Phytomedicine 2023;109:154570. [Crossref] [PubMed]
  2. Ding DX, Li Q, Shi K, et al. LncRNA NEAT1-miR-101-3p/miR-335-5p/miR-374a-3p/miR-628-5p-TRIM6 axis identified as the prognostic biomarker for lung adenocarcinoma via bioinformatics and meta-analysis. Transl Cancer Res 2021;10:4870-83. [Crossref] [PubMed]
  3. Feng RM, Zong YN, Cao SM, et al. Current cancer situation in China: good or bad news from the 2018 Global Cancer Statistics? Cancer Commun (Lond) 2019;39:22. [Crossref] [PubMed]
  4. Kleczko EK, Kwak JW, Schenk EL, et al. Targeting the Complement Pathway as a Therapeutic Strategy in Lung Cancer. Front Immunol 2019;10:954. [Crossref] [PubMed]
  5. Yang J, Peng A, Wang B, et al. The prognostic impact of lymph node metastasis in patients with non-small cell lung cancer and distant organ metastasis. Clin Exp Metastasis 2019;36:457-66. [Crossref] [PubMed]
  6. Luo T, Chen Q, Zeng J. Analysis of lymph node metastasis in 200 patients with non-small cell lung cancer. Transl Cancer Res 2020;9:1577-83. [Crossref] [PubMed]
  7. Dai C, Xie H, Kadeer X, et al. Relationship of Lymph Node Micrometastasis and Micropapillary Component and Their Joint Influence on Prognosis of Patients With Stage I Lung Adenocarcinoma. Am J Surg Pathol 2017;41:1212-20. [Crossref] [PubMed]
  8. Guo Q, Wu CY, Jiang N, et al. Downregulation of T-cell cytotoxic marker IL18R1 promotes cancer proliferation and migration and is associated with dismal prognosis and immunity in lung squamous cell carcinoma. Front Immunol 2022;13:986447. [Crossref] [PubMed]
  9. Gong WJ, Zhou T, Wu SL, et al. A novel immune-related ceRNA network that predicts prognosis and immunotherapy response in lung adenocarcinoma. Ann Transl Med 2021;9:1484. [Crossref] [PubMed]
  10. Du J, Zheng L, Chen S, et al. NFIL3 and its immunoregulatory role in rheumatoid arthritis patients. Front Immunol 2022;13:950144. [Crossref] [PubMed]
  11. Guo Q, Liu XL, Liu HS, et al. The Risk Model Based on the Three Oxidative Stress-Related Genes Evaluates the Prognosis of LAC Patients. Oxid Med Cell Longev 2022;2022:4022896. [Crossref] [PubMed]
  12. Huang H, Wu W, Lu Y, et al. The development and validation of a m6A-lncRNAs based prognostic model for overall survival in lung squamous cell carcinoma. J Thorac Dis 2022;14:4055-72. [Crossref] [PubMed]
  13. Zhang YQ, Li K, Guo Q, et al. A New Risk Model Based on 7 Quercetin-Related Target Genes for Predicting the Prognosis of Patients With Lung Adenocarcinoma. Front Genet 2022;13:890079. [Crossref] [PubMed]
  14. Guo Q, Xiao XY, Wu CY, et al. Clinical Roles of Risk Model Based on Differentially Expressed Genes in Mesenchymal Stem Cells in Prognosis and Immunity of Non-small Cell Lung Cancer. Front Genet 2022;13:823075. [Crossref] [PubMed]
  15. Shen A, Ye Y, Chen F, et al. Integrated multi-omics analysis identifies CD73 as a prognostic biomarker and immunotherapy response predictor in head and neck squamous cell carcinoma. Front Immunol 2022;13:969034. [Crossref] [PubMed]
  16. Shepelev MV, Korobko IV. The RHOV gene is overexpressed in human non-small cell lung cancer. Cancer Genet 2013;206:393-7. [Crossref] [PubMed]
  17. Zhang D, Jiang Q, Ge X, et al. RHOV promotes lung adenocarcinoma cell growth and metastasis through JNK/c-Jun pathway. Int J Biol Sci 2021;17:2622-32. [Crossref] [PubMed]
  18. Chen Y, Zhou H, Yang S, et al. Increased ABCC2 expression predicts cisplatin resistance in non-small cell lung cancer. Cell Biochem Funct 2021;39:277-86. [Crossref] [PubMed]
  19. Qian CY, Zheng Y, Wang Y, et al. Associations of genetic polymorphisms of the transporters organic cation transporter 2 (OCT2), multidrug and toxin extrusion 1 (MATE1), and ATP-binding cassette subfamily C member 2 (ABCC2) with platinum-based chemotherapy response and toxicity in non-small cell lung cancer patients. Chin J Cancer 2016;35:85. [Crossref] [PubMed]
  20. Lim S, Alshagga M, Ong CE, et al. Cytochrome P450 4B1 (CYP4B1) as a target in cancer treatment. Hum Exp Toxicol 2020;39:785-96. [Crossref] [PubMed]
  21. Liu X, Jia Y, Shi C, et al. CYP4B1 is a prognostic biomarker and potential therapeutic target in lung adenocarcinoma. PLoS One 2021;16:e0247020. [Crossref] [PubMed]
  22. Zhao W, Zhang LN, Wang XL, et al. Long noncoding RNA NSCLCAT1 increases non-small cell lung cancer cell invasion and migration through the Hippo signaling pathway by interacting with CDH1. FASEB J 2019;33:1151-66. [Crossref] [PubMed]
  23. Yao D, Gu P, Wang Y, et al. Inhibiting polo-like kinase 1 enhances radiosensitization via modulating DNA repair proteins in non-small-cell lung cancer. Biochem Cell Biol 2018;96:317-25. [Crossref] [PubMed]
  24. Wang Y, Zhang J, Li YJ, et al. MEST promotes lung cancer invasion and metastasis by interacting with VCP to activate NF-κB signaling. J Exp Clin Cancer Res 2021;40:301. [Crossref] [PubMed]
  25. Zhang Q, Liu H, Zhang J, et al. MiR-142-5p Suppresses Lung Cancer Cell Metastasis by Targeting Yin Yang 1 to Regulate Epithelial-Mesenchymal Transition. Cell Reprogram 2020;22:328-36. [Crossref] [PubMed]
  26. Chao CC, Lee WF, Wang SW, et al. CXC chemokine ligand-13 promotes metastasis via CXCR5-dependent signaling pathway in non-small cell lung cancer. J Cell Mol Med 2021;25:9128-40. [Crossref] [PubMed]
  27. Sun J, Wang X, Liu W, et al. Novel evidence for retinoic acid-induced G (Rig-G) as a tumor suppressor by activating p53 signaling pathway in lung cancer. FASEB J 2020;34:11900-12. [Crossref] [PubMed]
  28. Zhang YQ, Yuan Y, Zhang J, et al. Evaluation of the roles and regulatory mechanisms of PD-1 target molecules in NSCLC progression. Ann Transl Med 2021;9:1168. [Crossref] [PubMed]
  29. Gao L, Xia L, Ji W, et al. Knockdown of CDK5 down-regulates PD-L1 via the ubiquitination-proteasome pathway and improves antitumor immunity in lung adenocarcinoma. Transl Oncol 2021;14:101148. [Crossref] [PubMed]
  30. Song R, Liu X, Zhu J, et al. RhoV mediates apoptosis of RAW264.7 macrophages caused by osteoclast differentiation. Mol Med Rep 2015;11:1153-9. [Crossref] [PubMed]

(English Language Editor: J. Jones)

Cite this article as: Li Q, Liu XL, Jiang N, Li QY, Song YX, Ke XX, Han H, Luo Q, Guo Q, Luo XY, Chen C. A new prognostic model for RHOV, ABCC2, and CYP4B1 to predict the prognosis and association with immune infiltration of lung adenocarcinoma. J Thorac Dis 2023;15(4):1919-1934. doi: 10.21037/jtd-23-265

Download Citation