Back to Journals » OncoTargets and Therapy » Volume 14

CBX7, a Potential Prognostic Biomarker in Lung Adenocarcinoma

Authors Yang Y, Hu Z, Sun H, Yu Q, Yang L, Yin F, Sun Y, Pu L, Zhu X, Li S, Chen X, Zhao Y

Received 22 June 2021

Accepted for publication 28 September 2021

Published 24 December 2021 Volume 2021:14 Pages 5477—5492

DOI https://doi.org/10.2147/OTT.S325203

Checked for plagiarism Yes

Review by Single anonymous peer review

Peer reviewer comments 4

Editor who approved publication: Dr Arseniy Yuzhalin



Yanlong Yang,1,* Zaoxiu Hu,2,* Hongwen Sun,1 Qinghe Yu,1 Linzhu Yang,1 Fang Yin,1 Yongmen Sun,1 Lisha Pu,3 Xingming Zhu,1 Shen Li,4 Xiaobo Chen,1 Yunping Zhao1

1Department of Thoracic Surgery, The First Affiliated Hospital of Kunming Medical University, Kunming, 650032, People’s Republic of China; 2Department of Pathology, The Third Affiliated Hospital of Kunming Medical University (Yunnan Cancer Hospital, Yunnan Cancer Center), Kunming, 650118, People’s Republic of China; 3Department of Nephrology, Kunming Yanan Hospital, Kunming, 650051, People’s Republic of China; 4Department of Thoracic Surgery, Baoshan People’s Hospital, Baoshan, 678000, People’s Republic of China

*These authors contributed equally to this work

Correspondence: Xiaobo Chen; Yunping Zhao
Department of Thoracic Surgery, The First Affiliated Hospital of Kunming Medical University, No. 295 Xichang Road, Wuhua District, Kunming, 650032, Yunnan Province, People’s Republic of China
Tel/Fax +86 87165324888
Email [email protected]; [email protected]

Background: Lung adenocarcinoma (LUAD) is a major type of NSCLC and has high morbidity and mortality. The identification of useful prognostic biomarkers for LUAD is important. CBX7 has been reported in various cancers yet its expression level and potential roles have not been fully understood.
Methods: GEPIA, Oncomine, TCGA, KM plotter and OSluca databases were used to explore the expression profile and prognostic effects of CBX7 mRNA expression in patients with LUAD. TIMER was used to explore the relationship between CBX7 and immune infiltrating cells. GSEA was used to further explore the potential biological process and pathways regulated by CBX7 in LUAD. Lastly, IHC detection of CBX7 in 95 samples was used to validate the result.
Results: We found CBX7 was downregulated in LUAD in GEPIA, Oncomine and TCGA databases. TCGA, KM plotter and OSluca databases suggested that CBX7 was associated with poor clinical outcomes and low survival rate. Using TIMER, we found that CBX7 might be associated with immune infiltration. Via gene set enrichment analysis, we found that tumor-associated biological processes and signaling pathways were enriched in the CBX7 downregulated group. Using clinical samples, we found that CBX7 protein has low expression in LUAD and was associated with poor survival.
Conclusion: CBX7 might serve as a promising biomarker and potential molecular target in LUAD.

Keywords: CBX7, LUAD, prognosis, biomarker, immune infiltration

Introduction

Lung cancer has the highest morbidity and mortality not only in China but also worldwide. About 80% of lung cancers are non-small-cell lung cancers.1 NSCLC consisted of 3 subtypes, including lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), and large cell carcinoma.2 Among these, LUAD has been widely accepted as the dominant subtype, especially among Asian women.3 Though various efforts have been made, the prognosis of lung cancer remains poor. In addition, clinical outcomes may vary even when similar treatments are applied to patients with the same histology and TNM stages. Such instability could be caused by genetic intra-tumor heterogeneity, which has been a major obstacle for lung cancer treatment.4 As considerable progress is being made regarding molecular biological techniques, such as next-generation sequencing, it is now possible to further understand the molecular features in lung cancer.

Polycomb group (PcG) complexes are epigenetic regulatory complexes that conduct transcriptional repression of target genes via chromatin modification. Polycomb repressive complexes 1 and 2 (PRC1 and PRC2) are the two major forms of PcG, which are required for maintaining the stemness of embryonic stem cells and numerous types of adult stem cells. Abnormally increased expression of either PRC1 or PRC2 triggers aberrant silencing of corresponding target genes, which leads to the appearance of a broad spectrum of carcinomas.5,6 Chromobox7 (CBX7) belongs to the CBX protein family, which is involved in the organization of polycomb repressive complex 1 (PRC1). All of these proteins are involved in the regulation of heterochromatin, gene expression, and developmental progress.5 Recent studies have reported that CBX7 expression is decreased in many cancers, such as cervical cancer,7 pancreatic cancer,8 thyroid, colorectal, and breast cancer.5 However, the expression profile and prognostic value of CBX7 in LUAD samples remains unclear at present. This study aims to explore the differential expression and potential role of CBX7 in LUAD.

Therefore, in the current research, we investigated the expression profile and prognostic effects of CBX7 mRNA expression in patients with LUAD using various databases, including TCGA, Oncomine, GEPIA, OSluca and Kaplan–Meier plotter. TIMER was used to explore the relationship between CBX7 and immune infiltrating cells. GSEA was used to further explore the potential biological processes and pathways regulated by CBX7 in LUAD. Lastly, IHC detection of CBX7 in 95 samples was used to validate the result. Overall, our research preliminarily yet systematically characterizes the expression profiles of CBX7 in LUAD and reveals that the detection of the CBX7 expression status in LUAD patients is instructive and may be applied as potential biomarkers for prognostic assessment.

Materials and Methods

Data Mining in TCGA

The UCSC Xena browser (https://xenabrowser.net/) was used to download the Cancer Genome Atlas (TCGA) dataset of lung adenocarcinoma (LUAD). The expression level of CBX7 mRNA expression level in tumor and normal lung tissues was collected. In addition, the clinicopathological characters and survival status were collected.

Oncomine and GEPIA Database Analysis

538 samples from 5 Oncomine datasets (www.oncomine.org)9 were used to analyze the CBX7 expression in LUAD vs normal tissues. Gene Expression Profiling Interactive Analysis (GEPIA, http://gepia.cancer-pku.cn/)10 was used to explore the expression profile in pan-cancers.

Meta-Analyses of CBX7 Expression in GEO Dataset in OSluca Online Database

The OSluca online database (http://bioinfo.henu.edu.cn/LUCA/LUCAList.jsp)11 contains multiple lung cancer datasets, which can perform hazard ratio analysis on the impact of CBX7 in LUAD. The hazard ratio (HR) and its 95% confidence interval (CI) in each eligible dataset were extracted for meta-analysis. The combined HRs were estimated graphically by Forest plots. Heterogeneity between studies was detected by the Q test and the I2 metric (no heterogeneity: I2 = 0–25%; moderate heterogeneity: 25–50%; large heterogeneity: 50–75%; and extreme heterogeneity: 75–100%). If P≥0.10 in the Q test or I2 < 50%, the fixed-effect model (the Mantel Haenszel method)12 was used. Otherwise, a random effect model analysis was conducted.13 Subgroup analyses based on sex, age, smoking status and stage were conducted. In addition, publication bias was assessed by the method reported by Begg and Mazumdar.14 The meta-analysis was conducted by STATA v11.2 (StataCorp, College Station, Texas, USA).

Data Mining in Kaplan-Meier Plotter

Kaplan-Meier plotter (www.kmplot.com)15 was used to analyze the relationship between CBX7 mRNA expression and overall survival (OS) as well as recurrence-free survival (RFS). The hazard ratio (HR) with 95% confidence intervals (CI) and log-rank P value were calculated

The Relationship Between CBX7 and Immune Microenvironment by TIMER

TIMER(https://cistrome.shinyapps.io/timer/)16 is a web resource for systematical evaluations of the clinical impact of different immune cells among diverse cancer types. The level of abundance regarding six immune cell types: B cell, CD4 T cell, CD8 T cell, neutrophil, macrophage and dendritic cell in the tumor microenvironment can be estimated. Gene markers of tumor-infiltrating immune cells included markers of CD8+T cells, T cells (general), B cells, monocytes, TAMs, M1macrophages, M2 macrophages, neutrophils, natural killer (NK)cells, dendritic cells (DCs), T-helper 1 (Th1) cells, T-helper 2(Th2) cells, follicular helper T (Tfh) cells, T-helper 17 (Th17) cells, Tregs, and exhausted T cells. These gene markers have been referred in prior studies. In current study, we used TIMER to explore the relationship between CBX7 and gene markers of immune infiltrating cells in LUAD.

Gene Set Enrichment Analysis (GSEA)

GSEA was used to enrich the mRNAs predicted to be correlated with pathways. P < 0.05 was considered to indicate statistical significance, with reference to gene sets from the Molecular Signatures Database (MSigDB) of GO gene sets and KEGG gene sets, respectively. The number of permutations was set at 1000. Enrichment results satisfying a numerical P-value <0.05 and a false discovery rate (FDR) <0.25 were considered to be statistically significant.

Immunohistochemistry (IHC) Staining

95 LUAD tissues and corresponding adjacent normal lung tissues were collected from patients receiving surgery at the Department of Thoracic Surgery, the First Affiliated Hospital of Kunming Medical University from Jan 2015 and Dec 2016. Anti-CBX7 (1:100 dilution, Cat. 14876-1-AP, ProteinTech, Wuhan, China) were used for IHC staining. Antibody staining was visualized with DAB and hematoxylin counterstain (ZSGB-Bio). The percentage of positively stained cells was scored on a scale of 0 to 4 as follows: 0: < 1%, 1: 1–25%, 2: 25–50%, 3: 50–75% and 4: > 75%. The staining intensity was scored from 0 to 3 as follows: 0: negative, 1: weak, 2: moderate, and 3: strong. The immunoreactivity score (IRS) for each case was generated by multiplying percentages of positive cells and staining intensities. Immuno-stained sections were photographed using a microscope (Olympus Corporation, Tokyo, Japan).

Statistical Analysis

Kaplan-Meier survival analysis using log rank test was used to calculate the median survival time (MST). Univariate and multivariate Cox proportional hazards models were used to calculate the HR and 95% CI. Statistical analyses were performed using R software v3.5.0 (R Foundation for Statistical Computing, Vienna, Austria) and GraphPad Prism v7.00 (GraphPad Software Inc., USA). P < 0.05 was considered statistically significant. Receiver operating characteristic (ROC) analysis for death and recurrence detection was performed to identify the Youden index for CBX7 expression (the cutoff to separate patients) in Kaplan Meier curves of OS and RFS.

Result

Expression Profile and Clinicopathological Characteristics of CBX7 in LUAD in Oncomine and TCGA Datasets

Five datasets17–21 from Oncomine were used to evaluate the CBX7 expression in LUAD and normal lung tissue, all of which demonstrated CBX7 was downregulated in LUAD compared with normal lung tissues. Fold change (FC) of LUAD exceeded 2 in 4 studies (Table 1). Furthermore, data of 59 normal lung tissues and 526 LUAD tissues were obtained from TCGA dataset to evaluate the expression and clinicopathological characteristics of CBX7 in LUAD. Expression level of CBX7 was found to be significantly downregulated in LUAD compared with normal lung tissues (P<0.001) (Figure 1A).

Table 1 Comparison of CBX7 in LUAD and Normal Lung Tissues Across 5 Datasets in Oncomine

In addition, the expression profile of CBX7 in pan-cancer was also investigated by Oncomine and GEPIA databases. We found CBX7 was not only downregulated in LUAD but also downregulated in most cancers (Supplementary Figure S1), indicated CBX7 may play a tumor suppressor gene in many cancers.

495 LUAD patients with full clinicopathological characteristics were further analyzed. Low CBX7 expression was found to be more common in male, age <65 years, as well as smokers. Low CBX7 expression seemed to be associated with increased tumor size, lymph node metastasis, distant metastasis and advanced TNM stage (P<0.05) (Figure 1). ROC analysis yielded the cutoff value for CBX7 as 2.0119. According to this cut-off value, 495 LUAD patients were divided into two groups, low CBX7 group (n=246) and high CBX7 group (n=249). Similar results were obtained when continuous data was transformed into binary data. CBX7 expression level was found to be associated with age, sex, tumor size (T), lymph node metastasis (N), distant metastasis (M) and TNM stage. In addition, CBX7 expression was associated with recurrence and living status (Table 2).

Table 2 The Relationship Between CBX7 and the Clinicopathological Characteristics of Patients with LUAD in TCGA Dataset

Figure 1 The expression level of CBX7 in LUAD and its relationship with clinicopathological characteristics in TCGA. (A). CBX7 was downregulated in LUAD when compared with normal tissues. (BH) the relationship between CBX7 expression and sex (B), age (C), smoking status (D), T (E), N (F), M (G) and stage (H). *, ** and *** represented P <0.05, 0.01 and 0.001 respectively. NS represented not significance.

CBX7 Expression and OS in LUAD by Meta-Analysis of OSluca Online Database

A total of 20 GEO datasets including 2864 LUAD patients from OSluca were included for analysis. The pooled analysis suggested reduced CBX7 was associated with worse OS (HR=0.612, 95% CI = 1.39–1.80, p<0.001), with no heterogeneity between studies (I2=18.7%, P=0.222) (Figure 2A and Table 3). The shape of the funnel plot was symmetrical for all comparisons suggesting no publication bias detected (Figure 2B). In subgroup analysis, most subgroups such as sex, age, smoking status and stage all suggested a significant association between low CBX7 and poor OS with no heterogeneity (Table 3). Eight representative KM plots from OSluca online database (GSE13213, GSE30219, GSE31210, GSE3141, GSE41271, GSE50081, GSE67639 and GSE68465) are showed in Figure 2C.

Table 3 Summarized Meta-Analysis Results of the Impact of CBX7 Expression on Overall Survival in Patients with LUAD

Figure 2 Meta-analyses of CBX7 expression in GEO dataset in OSluca online database. (A). Forrest plot of hazard ratio (HR) for the association of CBX7 expression with overall survival (OS) in LUAD. (B).The funnel plot of the meta-analysis of the impact of CBX7 expression on overall survival (OS) in LUAD. (C). Kaplan-Meier curves of OS for CBX7 expression in LUAD in GSE13213, GSE30219, GSE31210, GSE3141, GSE41271, GSE50081, GSE67639 and GSE68465.

Figure 3 Kaplan-Meier curves of overall survival (OS) and recurrence-free survival (RFS) for CBX7 expression in LUAD. (A and B) Kaplan-Meier curves of OS (A) and RFS (B) in KM plotter. (C and D) Kaplan-Meier curves of OS (C) and RFS (D) in TCGA.

Prognostic Role of CBX7 in LUAD in KM Plotter and TCGA

Firstly, KM plotter was used to assess the prognostic role of CBX7 in LUAD. Low CBX7 expression was found to be associated with poor OS (HR=0.50, 95% CI=0.39–0.64, P=1.1×10−8) (Figure 3A) and RFS (HR=0.34, 95% CI=0.24–0.47, P=4×10−8) (Figure 3B). TCGA dataset analysis showed that patients with low CBX7 expression tend to have worse OS (Figure 3C) and RFS (Figure 3D). Moreover, univariate and multivariate Cox regression analyses were conducted. In univariate analysis, T, N, M pathologic stage and low CBX7 protein level were correlated with worse OS, while T, N pathologic stage and low CBX7 protein level were correlated with worse RFS (Table 4). Further multivariate showed that only T, N and CBX7 were associated with poor OS and RFS (Table 4).

Table 4 Univariate and Multivariate Analysis of OS and RFS in LUAD Patients in TCGA Dataset

CBX7 Expression Was Associated with Tumor-Infiltrating Immune Cells in LUAD

The relationship between CBX7 and immune infiltration in LUAD was evaluated by TIMER. CBX7 was found to be negatively correlated with tumor purity (r=−0.246, P=3.2e-8). There was a positive correlation between CBX7 expression and B cell (r=0.492, P=7.52e-31), CD8+ T cell (r=0.16, P=2.97e-4), CD4+ T cell (r=0.49, P=1.33e-33), macrophage (r=0.265, P=3e-9), neutrophil (r=0.215, P=1.81e-6), and dendritic cell (r=0.383, P=1.75e-18) (Figure 4). In addition, various immune marker genes corresponding to different tumor-infiltrating immune cells (CD8+ T cells, T cells (general), B cells, monocytes, TAMs, M1 and M2 macrophages, neutrophils, NK cells, Th1, Th2, Tfh, Th17, Treg and T cell exhaustion) were also investigated by TIMER and CBX7 were found to be associated with most of these marker genes. Generally, CBX7 expression were positively correlated with the majority of gene markers of different functional T cells (CD8+ T, Th1, Th2, Tfh, Th17, Treg and T cell exhaustion), B cells, dendritic cells, neutrophils and monocyte after purifying adjustments. CBX7 was related to tumor-associated macrophage (TAM, M1 and M2 macrophage) infiltration (Table 5), which suggested that CBX7 play critical roles in immune microenvironment in LUAD.

Table 5 The Relationship Between CBX7 and Relate Genes and Markers of Immune Cells in TIMER

Figure 4 CBX7 expression level has significant negative correlations with infiltrating levels of B cell, CD4+ T cells, Macrophages, Neutrophils and DCs in LUAD.

GSEA Result of CBX7

Potential mechanism of CBX7 was explored and analyzed by GSEA. Pathways that were observed to be significantly enriched between CBX7 high 10% group and low 10% groups were analyzed. The result suggested that CBX7 may participate in Aminoacyl tRNA biosynthesis, base excision repair, cell cycle, citrate cycle TCA cycle, DNA replication, homologous recombination, mismatch repair, and nucleotide excision repair pathways (Figure 5 and Table 6).

Table 6 Gene Set Enrichment Analysis Result of CBX7 Low Expression Group

Figure 5 Gene Set Enrichment Analysis (GESA) results for tumors with low CBX7 expression in LUAD.

IHC Validation of CBX7 in LUAD

95 paraffin-embedded LUAD and corresponding adjacent normal lung tissues were included for analysis. The IHC scores of paired LUAD tissues and their adjacent normal tissues showed that protein level of CBX7 is down-regulated in 95 LUAD patients (P<0.001) (Figure 6 and Figure 7A). Low CBX7 protein level was found to be correlated with differentiation degree, tumor size and lymph node metastasis (N). No significant association was found in region, sex, smoking status and TNM stage (Figure 7). The cutoff value yielded by ROC analysis for CBX7 IRS was 7.5, according to which patients were divided into low CBX7 (n=46) and high CBX7 groups (n=49). CBX7 expression was associated with living status, no significant association was found in age, region, sex, smoking history, differentiation degree, T, N, M, TNM stage and recurrence status (Table 7). The prognostic role of CBX7 in LUAD was evaluated by univariate and multivariate Cox regression analyses. In univariate analysis, T, N, M, pathologic stage and low CBX7 protein level were correlated with the worst OS, while differentiation degree, N, M, pathologic stage and low CBX7 protein level were correlated with unfavorable RFS (Table 8 and Figure 8). In multivariate analysis, only pathologic M and CBX7 protein expression were associated with poor OS and RFS (Table 8).

Table 7 The Relationship Between CBX7 and the Clinicopathological Characteristics in 95 patients with LUAD

Table 8 Univariate and Multivariate Analysis of OS and RFS in 95 Patients with LUAD

Figure 6 Immunostaining of CBX7 in normal lung tissues and LUAD.

Figure 7 The expression level of CBX7 and its relationship with clinicopathological characteristics in 95 patients with LUAD. (A) CBX7 was downregulated in LUAD when compared with normal tissues. (BH) The relationship between CBX7 expression and region (B), sex (C), smoking status (D), differentiation (E), T (F), N (G) and stage (H). *, ** and *** represented P <0.05, 0.01 and 0.001 respectively. NS represented not significance.

Figure 8 Kaplan-Meier curves of overall survival (OS)(A) and recurrence-free survival(RFS) (B) for CBX7 expression in 95 patients with LUAD.

Discussion

Our current study comprehensively investigated the prognostic role of CBX7 in LUAD. Through online database and our individual dataset, we demonstrated that CBX7 was down-regulated in LUAD tissues and might serve as an independent prognostic biomarker in LUAD.

LUAD and LUSC were the major types for NSCLC. With increasing prevalence and its dominant histology, LUAD has been selected as the focus of this study. Another important reason stands as to avoid bias that could possibly arise due to different biological processes and molecular profiles between LUAD and LUSC. Further studies are necessary to investigate the roles of CBX7 in LUSC and other subtypes in order to confirm whether CBX7 could serve as a potential biomarker in lung cancer treatment. As previous studies reported that CBX7 was downregulated in many cancers including hepatocellular carcinoma, glioma, pancreatic cancer, colon, urothelial and thyroid cancer, such characteristic was also associated with poor prognosis in hepatocellular, glioma, pancreatic cancer and colon cancer. As far as human lung carcinoma is concerned, no detectable level of CBX7 has been detected in the samples mentioned, whereas CBX7 is abundantly expressed in normal lung tissue. Interestingly, LOH was reported in 50% of informative cases. Additionally, the expression of CBX7 was lost in the “normal” area surrounding carcinomas (non-pathological at the morphological evaluation) in 50% of the analyzed cases. Intriguingly, these cases also involved LOH, suggesting that the decreased expression of CBX7 could be involved in the transition of epithelial lung cells toward a fully transformed and malignant phenotype. In this study, the expression profile and the prognostic role of CBX7 in LUAD have been comprehensively investigated using both online database and private database. Experimental results suggested that low CBX7 was associated with aggressive phenotypes, such as stage and definition, which is consistent with TCGA and private database analysis. Hence, this study is presented with rounded and scientifically sound evidence.

Reduced CBX7 expression was associated with aggressive phenotypes and poor survival status among lung cancer and other cancers, suggesting that CBX7 serves as a novel tumor-suppressive gene. Pallante et al stated that the loss of CBX7 promoted cancer development in mice, and CBX7-knockout (CBX7-KO) mice were more likely to develop liver and lung neoplasia compared with wide-type. Mouse embryonic fibroblasts (MEFs) derived from CBX7-KO mice was used to explore the effect of CBX7 regarding cellular proliferation and CBX7+/– and CBX7–/– MEFs were found to grow faster than CBX7+/+ MEF. In addition, the percentage of S-phase cell in CBX7-KO MEF was higher than wild-type MEF.22 Several studies have shed light on the mechanism of reduced CBX7 promoting cancer development. It has been reported that CBX7 was able to sustain the expression of the E-cadherin gene (CDH1) through binding to its promoter, as well as counteracting with anti-transcriptional effect of histone deacetylase 2 (HDAC2) and Protein Arginine Methyltransferase 1(PRMT1).23,24 As a result, the reduced expression of CBX7 decreases E-cadherin expression, which then leads to the critical event epithelial-to-mesenchymal transition (EMT).23,24 CBX7 binds the E-box to inhibit TWIST-1 function, tumorigenicity, and metastatic potential. Deletion of CBX7 in sEOC can sufficiently reactivate TWIST-1-induced transcription, prompt mesenchymal transformation, and enhanced tumorigenicity in vivo.25 CBX7 interacts with HMGA1b, which binds to and regulates the acetylation of CCNE1 promoter. As a coding gene for cyclin E, CCNE1 regulates cell cycle that is crucial for G1/S transition. Another mechanism through which CBX7 involves with cancer progression is that CBX7 negatively modulates the expression of SPP1 gene that encodes the osteopontin protein, a chemokine promoting cell migration.23,26

This study also investigated the potential role of CBX7 in LUAD via GSEA analysis. Several cancer-related pathways were identified including aminoacyl tRNA biosynthesis, base excision repair, cell cycle, citrate cycle TCA cycle, DNA replication, homologous recombination, mismatch repair, and nucleotide excision repair pathways. All of these pathways have been proved to be closely related to cancer development, and some, for example, cell cycle pathways, have been reported and investigated by previous studies.22,23 How CBX7 participate in DNA damage repair pathways (base excision repair, DNA replication, homologous recombination, mismatch repair, and nucleotide excision repair pathways)? PcG proteins participate in DNA damage response (DDR). The majority of PcG proteins, including EZH2 and the CBXs, localize to regions of damaged chromatin in order to recruit DNA damage repair machinery.27,28 Inhibiting CBXs may enhance chemotherapeutic response by preventing DNA damage repair, allowing massive accumulation of DNA damage.29 However, the role of CBX7 in the process is still unknown. Further studies need to be conducted for verification.

Tumor microenvironments, especially tumor-infiltrating immune cells, play an important role in the prognosis of lung cancer. This study demonstrated that CBX7 affects and mediates immune infiltration level in LUAD. CBX7 expression was positively related to the level of B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils and DCs infiltration. CBX7 expression was also positively related to most gene markers of T cells, B cells, tumor-associated macrophages (TAM), M1, M2, DC cell, and T helper cells (Th1, Th2, Tfh, and Th17). High level of these tumor-infiltrating immune was associated with better survival in lung cancer. CBX7 could regulate immune infiltration level by mediating gene expression in these cells. Also, in the pathway enriched analysis, we found that CBX may participate in DNA damage repair pathways (base excision repair, DNA replication, homologous recombination, mismatch repair, and nucleotide excision repair pathways). The deficiency of DNA damage repair pathways may result in the high tumor mutation burden (TMB), which may be associated with immune infiltration. However, the precise role of CBX7 in tumor immune microenvironment and tumor progression still needs exploration through future studies.

To validate the potential role of CBX7, 95 paraffin-embedded LUAD and corresponding adjacent normal lung tissue samples were included for analysis. CBX7 was found to be downregulated in LUAD tissues and was associated with malignant phenotypes such as T/N/M. In addition, multiple cox analysis showed that CBX7 is an independent prognostic factor, indicating that CBX7 might serve as a useful prognostic biomarker in LUAD.

In conclusion, this study demonstrated that CBX7 is downregulated in LUAD, which is associated with unappealing clinicopathological features and poor prognosis proved by both online datasets and our clinical samples. CBX7 stands as a crucial mediating factor regarding LUAD tumorigenesis due to its correlation with several cancer-related pathways and its regulating effects on immune infiltration level. Hence, CBX7 might serve as a promising biomarker and potential molecular target regarding clinical treatments of LUAD. Future large cohort and mechanism studies are necessary to validate and further reveal the molecular mechanisms and pathways of CBX7 in LUAD.

Ethical Approval and Consent to Participate

95 surgical LUAD samples used in this manuscript have been conducted in accordance with the principles of the 1964 Helsinki declaration and its later amendments or comparable ethical standards. Our study was also approved by the Ethics Committee of the First Affiliated Hospital of Kunming Medical University. All patients were obtained with informed consent.

Author Contributions

Yunping Zhao and Xiaobo Chen: conception and design of the work. Yanlong Yang and Zaoxiu Hu: conception, design of the work and have drafted the work. Linzhu Yang: the acquisition, analysis. Hongwen Sun and Qinghe Yu: substantively revised. Fang Yin, Yongmeng Sun, Xinming Zhu, Shen Li, Zaoxiu Hu and Lisha Pu collected data and analyzed. All authors contributed to data analysis, drafting or revising the article, have agreed on the journal to which the article will be submitted, gave final approval of the version to be published, and agree to be accountable for all aspects of the work.

Funding

This study was supported by Scientific Research Foundation of Yunnan Education Department (2019J1234 and 2019J1282), Yunnan Provincial Department of Science and Technology-Kunming Medical University Joint Foundation (202001AY070001-144 and 202101AY070001-172), Doctoral Research Foundation of the First Affiliated Hospital of Kunming Medical University (2018BS012) and Science and Technology Personnel Training Program of High-level Health Technicians in Yunnan Province (H-2017014).

Disclosure

The authors have no conflicts of interest to declare.

References

1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424. doi:10.3322/caac.21492

2. Herbst RS, Heymach JV, Lippman SM. Lung cancer. N Engl J Med. 2008;359(13):1367–1380. doi:10.1056/NEJMra0802714

3. Zappa C, Mousa SA. Non-small cell lung cancer: current treatment and future advances. Transl Lung Cancer Res. 2016;5(3):288–300. doi:10.21037/tlcr.2016.06.07

4. Denisenko TV, Budkevich IN, Zhivotovsky B. Cell death-based treatment of lung adenocarcinoma. Cell Death Dis. 2018;9(2):117. doi:10.1038/s41419-017-0063-y

5. Pallante P, Forzati F, Federico A, Arra C, Fusco A. Polycomb protein family member CBX7 plays a critical role in cancer progression. Am J Cancer Res. 2015;5(5):1594–1601.

6. Ma RG, Zhang Y, Sun TT, Cheng B. Epigenetic regulation by polycomb group complexes: focus on roles of CBX proteins. J Zhejiang Univ Sci B. 2014;15(5):412–428. doi:10.1631/jzus.B1400077

7. Li R, Yan Q, Tian P, et al. CBX7 inhibits cell growth and motility and induces apoptosis in cervical cancer cells. Mol Ther Oncolytics. 2019;15:108–116. doi:10.1016/j.omto.2019.09.002

8. Ni S, Wang H, Zhu X, et al. CBX7 suppresses cell proliferation, migration, and invasion through the inhibition of PTEN/Akt signaling in pancreatic cancer. Oncotarget. 2017;8(5):8010–8021. doi:10.18632/oncotarget.14037

9. Rhodes DR, Yu J, Shanker K, et al. ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia. 2004;6(1):1–6. doi:10.1016/S1476-5586(04)80047-2

10. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–W102. doi:10.1093/nar/gkx247

11. Yan Z, Wang Q, Lu Z, et al. OSluca: an interactive web server to evaluate prognostic biomarkers for lung cancer. Front Genet. 2020;11:420. doi:10.3389/fgene.2020.00420

12. Mantel N, Haenszel W. Statistical aspects of the analysis of data from retrospective studies of disease. J Natl Cancer Inst. 1959;22(4):719–748.

13. DerSimonian R, Laird N. Meta-analysis in clinical trials. Control Clin Trials. 1986;7(3):177–188. doi:10.1016/0197-2456(86)90046-2

14. Begg CB, Mazumdar M. Operating characteristics of a rank correlation test for publication bias. Biometrics. 1994;50(4):1088–1101. doi:10.2307/2533446

15. Gyorffy B, Surowiak P, Budczies J, Lanczky A. Online survival analysis software to assess the prognostic value of biomarkers using transcriptomic data in non-small-cell lung cancer. PLoS One. 2013;8(12):e82241. doi:10.1371/journal.pone.0082241

16. Li T, Fan J, Wang B, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–e110. doi:10.1158/0008-5472.CAN-17-0307

17. Hou J, Aerts J, den Hamer B, et al. Gene expression-based classification of non-small cell lung carcinomas and survival prediction. PLoS One. 2010;5(4):e10312. doi:10.1371/journal.pone.0010312

18. Su LJ, Chang CW, Wu YC, et al. Selection of DDX5 as a novel internal control for Q-RT-PCR from microarray data using a block bootstrap re-sampling scheme. BMC Genomics. 2007;8:140. doi:10.1186/1471-2164-8-140

19. Selamat SA, Chung BS, Girard L, et al. Genome-scale analysis of DNA methylation in lung adenocarcinoma and integration with mRNA expression. Genome Res. 2012;22(7):1197–1211. doi:10.1101/gr.132662.111

20. Landi MT, Dracheva T, Rotunno M, et al. Gene expression signature of cigarette smoking and its role in lung adenocarcinoma development and survival. PLoS One. 2008;3(2):e1651. doi:10.1371/journal.pone.0001651

21. Okayama H, Kohno T, Ishii Y, et al. Identification of genes upregulated in ALK-positive and EGFR/KRAS/ALK-negative lung adenocarcinomas. Cancer Res. 2012;72(1):100–111. doi:10.1158/0008-5472.CAN-11-1403

22. Forzati F, Federico A, Pallante P, et al. CBX7 is a tumor suppressor in mice and humans. J Clin Invest. 2012;122(2):612–623. doi:10.1172/JCI58620

23. Federico A, Pallante P, Bianco M, et al. Chromobox protein homologue 7 protein, with decreased expression in human carcinomas, positively regulates E-cadherin expression by interacting with the histone deacetylase 2 protein. Cancer Res. 2009;69(17):7079–7087. doi:10.1158/0008-5472.CAN-09-1542

24. Federico A, Sepe R, Cozzolino F, et al. The complex CBX7-PRMT1 has a critical role in regulating E-cadherin gene expression and cell migration. Biochim Biophys Acta Gene Regul Mech. 2019;1862(4):509–521. doi:10.1016/j.bbagrm.2019.02.006

25. Li J, Alvero AB, Nuti S, et al. CBX7 binds the E-box to inhibit TWIST-1 function and inhibit tumorigenicity and metastatic potential. Oncogene. 2020;39(20):3965–3979. doi:10.1038/s41388-020-1269-5

26. Sepe R, Formisano U, Federico A, et al. CBX7 and HMGA1b proteins act in opposite way on the regulation of the SPP1 gene expression. Oncotarget. 2015;6(5):2680–2692. doi:10.18632/oncotarget.2777

27. Chou DM, Adamson B, Dephoure NE, et al. A chromatin localization screen reveals poly (ADP ribose)-regulated recruitment of the repressive polycomb and NuRD complexes to sites of DNA damage. Proc Natl Acad Sci U S A. 2010;107(43):18475–18480. doi:10.1073/pnas.1012946107

28. Gieni RS, Ismail IH, Campbell S, Hendzel MJ. Polycomb group proteins in the DNA damage response: a link between radiation resistance and “stemness”. Cell Cycle. 2011;10(6):883–894. doi:10.4161/cc.10.6.14907

29. Connelly KE, Martin EC, Dykhuizen EC. CBX chromodomain inhibition enhances chemotherapy response in glioblastoma multiforme. Yale J Biol Med. 2016;89(4):431–440.

Creative Commons License © 2021 The Author(s). This work is published and licensed by Dove Medical Press Limited. The full terms of this license are available at https://www.dovepress.com/terms.php and incorporate the Creative Commons Attribution - Non Commercial (unported, v3.0) License. By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted without any further permission from Dove Medical Press Limited, provided the work is properly attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.