frontiers

Check for updates

OPEN ACCESS

EDITED BY Lei Huang, University of Massachusetts Medical School, United States

REVIEWED BY Yuejia Deng, Texas A and M University, United States Yue Liu, The University of Texas at Austin, United States

*CORRESPONDENCE Dejun Cui ☒ cuidejun@gz5055.com

Chen Li ☒ chen.li.scholar@outlook.com

These authors have contributed equally to this work and share first authorship

RECEIVED 12 January 2025 ACCEPTED 17 February 2025 PUBLISHED 03 March 2025

CITATION

Chen L, He Y, Duan M, Yang T, Chen Y, Wang B, Cui D and Li C (2025) Exploring NUP62’s role in cancer progression, tumor immunity, and treatment response: insights from multi-omics analysis. Front. Immunol. 16:1559396. doi: 10.3389/fimmu.2025.1559396

COPYRIGHT

@ 2025 Chen, He, Duan, Yang, Chen, Wang, Cui and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Exploring NUP62’s role in cancer progression, tumor immunity, and treatment response: insights from multi-omics analysis

Lihong Chen1t, Youfu He21, Menghui Duan3, Tianwen Yang4, Yin Chen1, Bo Wang5, Dejun Cui1* and Chen Li6*

1Department of Gastroenterology, Guizhou Provincial People’s Hospital, Guiyang, Guizhou, China, 2Department of Cardiology, Guizhou Provincial People’s Hospital, Guiyang, Guizhou, China, 3Department of Critical Care Medicine, The First Hospital of China Medical University, Shenyang, Liaoning, China, 4Department of Orthopaedics, Guizhou Second People’s Hospital, Guiyang, Guizhou, China, 5Department of General Surgery, Guizhou Provincial People’s Hospital, Guiyang, Guizhou, China, 6Department of Pharmacy, the First Affiliated Hospital of Guangxi Medical University, Nanning, Guangxi, China

Background: NUP62, a key component of the nuclear pore complex, is closely associated with cellular functions and cancer progression. However, its expression patterns, prognostic value, and relationship with tumour immunity and drug sensitivity across multiple cancers have not been systematically studied. This study used multi-omics analyses combined with experimental validation in gastric cancer to investigate the expression, functional characteristics, and clinical relevance of NUP62 in cancer.

Methods: Data from TCGA, GTEx, and CPTAC databases were used to analyse the expression, mutation characteristics, and clinical associations of NUP62. Tools such as SangerBox, TIMER 2.0, and GSEA were employed to evaluate the relationship between NUP62 and the tumour immune microenvironment, as well as its involvement in signalling pathways. Immunohistochemistry and RT-PCR were used to validate the expression of NUP62 in gastric cancer tissues. PRISM and CTRP databases were utilised to assess the correlation between NUP62 expression and drug sensitivity.

Results: NUP62 was significantly upregulated in multiple cancers and was associated with poor prognosis in cancers such as clear cell renal carcinoma (KIRC), lower-grade glioma (LGG), and adrenocortical carcinoma (ACC), while playing a protective role in others, such as bladder cancer (BLCA) and stomach cancer (STAD). Functional analyses showed that NUP62 is involved in cell cycle regulation, DNA damage repair, and tumour immunity. High NUP62 expression was significantly correlated with increased infiltration of immune cells, such as macrophages and T cells, and a higher response rate to immunotherapy. Drug sensitivity analysis identified NUP62 as a marker of sensitivity to various chemotherapeutic agents. Validation experiments demonstrated that NUP62 mRNA and protein levels were significantly higher in gastric cancer tissues than in adjacent normal tissues.

Conclusions: NUP62 plays a critical role in multiple cancers and shows potential as a biomarker for cancer diagnosis, prognosis, and therapeutic response prediction. Its role in tumour immunity and signalling pathways highlights its potential as a target for immunotherapy and precision medicine.

KEYWORDS

Nup62, immune, pan cancer, gastric cancer, immune cell infiltration

1 Background

Cancer represents a major public health issue worldwide and is the second leading cause of death in the United States. In 2022, there will be nearly 20 million new cases of cancer, while 9.7 million people will die of cancer. About one in five men or women will develop cancer in their lifetime, while about one in nine men and one in 12 women will die of cancer (1). The rate of decline in cancer mortality has risen from about 1 percent per year in the 1990s, to 1.5 percent per year in the 2000s, to 2 percent per year between 2015 and 2020 (2).This trend reflects the increasing depth of human research into cancer, which is closely related to advancements in early cancer diagnosis and targeted therapy (3-5). Pan-cancer analysis, which integrates multiple databases, can aid in the identification of cancer biomarkers and therapeutic targets (6, 7).

Metabolic reprogramming is a hallmark of tumors, whereby tumors undergo reprogramming of nutrient acquisition and metabolic pathways to meet the bioenergetic, biosynthetic, and redox demands of malignant cells (8, 9). There exists a close interaction between metabolism and signaling pathways in cancer cells, and several signaling pathways associated with cell proliferation also regulate metabolic pathways that integrate nutrients into biomacromolecules (10). Consequently, certain cancer-related mutations enable cancer cells to acquire and metabolize nutrients in a manner conducive to proliferation, rather than efficiently producing ATP (11, 12). To fulfill the demands of proliferation, cancer cells utilize processes such as glycolysis, glutaminolysis, and fatty acid oxidation to meet their energy requirements and metabolic synthesis processes (13, 14). Studies have indicated that oxidative phosphorylation represents an emerging target in cancer therapy (15). Furthermore, immune cells within the tumor microenvironment (TME) exhibit metabolic shifts similar to the glycolytic metabolic profile, leading to competition for nutrients between cancer cells and tumor-infiltrating cells (16). On the other hand, metabolic disturbances characterized by hypoxia and elevated metabolite levels (especially lactate) in the TME can result in immunosuppression (17, 18). Simultaneously, metabolic dysregulation and imbalances in immune cells within the TME can drive immune evasion and compromise treatment outcomes (19).

The primary function of the nuclear pore complex (NPC) is to mediate transport between the nucleus and cytoplasm. Studies have

shown that nucleoporins exert effects independently of the NPC in both nuclear and cytoplasmic compartments, directly regulating gene expression and participating in the regulation of development and the cell cycle (20). Nucleoporin 62 (NUP62), a structural component of the NPC, may mediate the localization of glucocorticoid receptors to the nucleus after binding to steroids (21). Additionally, NUP62 can interact with OSBP-related protein 8 (ORP8), which has the ability to regulate hepatic lipogenesis and plasma lipid levels (22). Furthermore, numerous studies have demonstrated that NUP62 mediates selective nucleocytoplasmic transport and is associated with various cancers through chromosomal translocations that generate fusion proteins, alterations in protein expression levels, and single-point mutations (23- 27). Specifically, high expression of NUP62 contributes to preventing epidermal differentiation in squamous cell carcinomas originating from stratified epithelia (27, 28).

2 Materials and methods

2.1 Data collection and processing

We obtained expression levels and related clinical characteristics of NUP62 from the Cancer Genome Atlas (TCGA, http:// cancergenome.nih.gov/) and the Genotype-Tissue Expression (GTEx, https://gtexportal.org/) databases through the Xena platform at the University of California, San Diego. The mutation frequency of NUP62 in the TCGA cohort was calculated using the cBioPortal database (https://www.cbioportal.org/). Pan-cancer analysis of TCGA samples was conducted using the online bioinformatics tool SangerBox 3.0 (http://sangerbox.com/) (29). Various immune infiltration algorithms in the TIMER 2.0 database (http://timer.cistrome.org) were utilized to characterize the correlation between NUP62 expression and the tumor immune microenvironment.

Expression profiles of human normal tissues and cancer cell lines were obtained from the Human Protein Atlas (HPA, https:// www.proteinatlas.org/). Relevant chemotherapy data were retrieved from the Genomics of Drug Sensitivity in Cancer (GDSC, https:// www.cancerrxgene.org/), Cancer Therapeutics Response Portal (CTRP, http://portals.broadinstitute.org/ctrp/), and PRISM databases to illustrate the relationship between NUP62 expression and drug sensitivity. Cancer immune cycle data originated from the

Tracking Immune Phenotypes in Cancer (TIP, https:// biocc.hrbmu.edu.cn/TIP/), while Immune Phenotype Scores (IPS) data were obtained from the TCIA website (https://tcia.at/home).

2.2 Expression and variation analysis

To investigate whether there were differences in NUP62 expression between tumor and normal tissues, we first compared the expression of NUP62 mRNA between tumor and normal tissues using the “wilcox” test. Subsequently, a paired “wilcox” test was conducted on matched samples to validate protein expression at the external gene transcriptome level in the GEO and CPTAC databases. Additionally, we used the “gganatogram” R package to visualize the expression of NUP62 in different human organs. The cBioPortal website (http://www.cbioportal.org) served as a powerful tool to search for NUP62 mutation frequency, types, copy number alteration (CNA) data, and gene alteration traits. Furthermore, we employed the “pROC” R package to calculate the Area Under the Curve (AUC) value, illustrating the importance of NUP62 in pan- cancer diagnosis.

2.3 Survival and clinical outcome analysis

Survival data were sourced from the TCGA database. We analyzed the relationship between NUP62 expression and these prognostic indicators, including Overall Survival (OS), Disease- Specific Survival (DSS), Progression-Free Interval (PFI), and Disease-Free Interval (DFI), using the “survival” and “survminer” R packages. Kaplan-Meier (KM) analysis and univariate COX analysis were combined to assess whether NUP62 was a protective or risk factor, resulting in the creation of high- confidence survival landscapes for NUP62. Additionally, the “forestplot” R package was utilized for visual analysis of COX survival data.

2.4 Subtyping of NUP62 and immunotherapy analysis

Researchers conducted extensive immunogenomic analysis on over 10,000 tumor samples encompassing 33 different cancer types from the TCGA database. Across tumor types, they successfully identified six immune subtypes by evaluating macrophage or lymphocyte markers, Th1 to Th2 cell ratios, the range of genetic heterogeneity among tumors, aneuploidy, neoantigen burden range, total cell landscape, immune-regulatory gene expression, and prognosis. The six subtypes are introduced as follows:

C1 (Tissue Healing) subtype exhibits elevated angiogenic gene expression, a high proliferation rate, and significant acquired immune infiltration with a Th2 bias. Primary associated cancers include colorectal adenocarcinoma (COAD), rectal adenocarcinoma (READ), lung squamous cell carcinoma (LUSC), luminal A subtype of breast invasive ductal carcinoma (BRCA), classical subtype of head and neck

squamous cell carcinoma (HNSC), and chromosomally unstable gastrointestinal subtype.

C2 (IFN-y Dominant) subtype demonstrates the highest M1/ M2 macrophage polarization, strong CD8 signaling, and, like the C6 subtype, the highest TCR diversity. Additionally, the C2 subtype has a high proliferation rate, with primary associated cancers including hypermutated BRCA, gastric cancer, ovarian cancer, HNSC, and cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC).

C3 (Inflammatory) subtype shows elevated Th17 and Th1 gene expression but fails to effectively inhibit tumor cell proliferation. Like the C5 subtype, C3 has fewer aneuploidies and overall cellular copy number alterations compared to other subtypes. Primary associated cancers include most kidney cancers, prostate adenocarcinoma (PRAD), pancreatic adenocarcinoma (PAAD), and thyroid papillary carcinoma (THCA).

C4 (Lymphocyte-Depleted) subtype is characterized by prominent macrophage features, Th1 suppression, and high M2 responsiveness. Primary associated cancers include specific subtypes of adrenocortical carcinoma (ACC), pheochromocytoma and paraganglioma (PCPG), hepatocellular carcinoma (LIHC), and gliomas.

C5 (Immunologically Silent) subtype has the lowest lymphocyte count, the highest macrophage response, and is dominated by M2 macrophages. This subtype is primarily a subtype of low-grade gliomas (LGG) of the brain, including glioma CpG island methylator phenotype-high (CIMP-H), 1p/19q codeleted subtype, and fibrillary astrocytoma-like (PA-like) type. Additionally, the remaining types are mainly within the C4 subtype, while isocitrate dehydrogenase mutant types are more prevalent in the C5 subtype than in the C4 subtype.

C6 (TGF-ß Dominant) subtype is a smaller group composed of a mix of cancers and does not dominate any TCGA subtype. This subtype has the highest TGF-B signature and high lymphatic infiltration, with equal distribution of Type I and Type II T cells.

This study further explored the relationship between NUP62 and patient prognosis in external datasets using the BEST database and investigated the correlation between NUP62 expression and the immunotherapy response in cancer patients (30).

2.5 Pathway and mechanism of action analysis

To delve into the mechanism of action of NUP62-related pathways, we classified various tumor samples based on NUP62 expression levels (the top 30% as the high-expression group and the bottom 30% as the low-expression group). Subsequently, we employed Gene Set Enrichment Analysis (GSEA) to investigate the differential activation or inhibition states of 50 characteristic genomic signatures and 83 metabolic genomic signatures between the high and low NUP62 expression groups across different malignancies. Utilizing the “GSVA” R package, we quantitatively analyzed 14 functional status genomic signatures using the “z- score” algorithm. Based on the z-score values of each genomic

signature, we conducted subsequent data processing and analysis. Furthermore, we applied “Pearson” correlation analysis to assess the statistical correlation between the z-score of each genomic signature and NUP62 expression levels. Additionally, we identified genes that underwent significant changes in both high and low NUP62 expression groups.

2.6 Identification of chemicals interacting with NUP62

To explore the potential association between NUP62 expression and drug sensitivity, we leveraged the GSCA database (http:// bioinfo.life.edu.cn/web/GSCALite/). Through the GSCA platform, we obtained information on small molecule drugs from the GDSC and CTRP databases. Furthermore, we identified genes differentially expressed between the high and low NUP62 expression groups across various malignancies. To screen for biomarkers closely related to NUP62, we selected the top 150 significantly upregulated and downregulated genes. Simultaneously, we downloaded the CMAP_gene_signatures RData file, which contains 1,288 attributes associated with compound functional characteristics.

3 Results

3.1 Expression and mutation of NUP62 in humans

Using unpaired and paired sample data from the TCGA database, we conducted an in-depth analysis at the mRNA level and observed an upward trend in NUP62 expression across multiple cancer types, including Bladder Urothelial Carcinoma (BLCA), Breast Invasive Carcinoma (BRCA), COAD, Esophageal Carcinoma (ESCA), HNSC, Kidney Renal Clear Cell Carcinoma (KIRC), Kidney Renal Papillary Cell Carcinoma (KIRP), LIHC, Lung Adenocarcinoma (LUAD), LUSC, and Stomach Adenocarcinoma (STAD) (Figures 1A, C). Additionally, we further validated these findings in the CPTAC database (Figure 1B). Figure 1D depicts the expression profile of NUP62 in different human organs. Notably, NUP62 expression levels were significantly upregulated in most cancer types, while they exhibited a downward trend in testicular tissue. NUP62 displayed variable expression patterns in most cancers, with mutation sites distributed as shown in Figure 1E. Leveraging the cBioPortal (TCGA, Pan- cancer Atlas) database, we comprehensively assessed the pan-cancer mutational characteristics of the NUP62 gene. The results indicated that the most common types of variations in NUP62 included mutations, amplifications, and deep deletions, with Endometrial Cancer, Cervical Cancer, and Bladder Cancer having the highest mutation rates (Figure 1F). Subsequently, we conducted an in-depth analysis of the distribution characteristics of NUP62 mutations across different cancers and explored the types of single nucleotide variants (SNVs). The results showed that missense mutations dominated among various types of variations (Figure 1G). Finally,

we conducted a comprehensive analysis of the distribution of NUP62 and other genes in cancer samples. The results revealed that among all cancer types, NUP62 was most widely distributed in UCEC, while the TP53 signal was significantly distributed in most cancers (Figure 1H).

3.2 Expression of NUP62 in cells

Based on fluorescence images provided by the Human Protein Atlas (HPA), we observed that NUP62 is primarily localized within the cell nucleus (Supplementary Figure 1A). Researchers systematically analyzed normal tissues such as the cerebral cortex, liver, colon, kidney, and pancreas, as well as tumor tissues such as lung cancer, liver cancer, colorectal cancer, pancreatic cancer, and breast cancer, using immunohistochemical techniques. The results indicated that compared to normal tissues, the expression level of NUP62 in tumor tissues exhibited a significant upward trend (Supplementary Figure 1B).

3.3 Correlation between NUP62 and clinical characteristics across various cancer types

As shown in Figure 2A, gender differences in NUP62 expression levels are observed in sarcoma (SARC), KIRP, and the Pan-kidney cohort (KIPAN), with significantly higher expression in female patients compared to male patients. Additionally, the expression level of NUP62 is closely related to age. In KIRP, PCPG, and ESCA, there is a negative correlation between age and NUP62 expression (Figure 2G). Furthermore, the expression of NUP62 is associated with various cancer stages. In KIPAN and LIHC, differences in NUP62 expression are observed among tumor patients with different T stages (Figure 2B). In KIRP and ACC, NUP62 expression exhibits certain differences across different N stages (Figure 2C). Only in ACC is there an association between NUP62 expression and its M stage (Figure 2D). In glioma (GBMLGG), LGG, HNSC, and LIHC, significant differences in NUP62 expression are found among tumor patients with different G stages (Figure 2E). In LIHC, NUP62 expression is correlated with its stage classification (Figure 2F).

3.4 Diagnostic value of NUP62 in pan-cancer

We evaluated the diagnostic capacity of NUP62 for various cancers in both the TCGA dataset and the combined TCGA-GTEx dataset (Figure 3A). The results indicated that the AUC value of NUP62 in ESCA was >0.9, suggesting a high diagnostic value for ESCA. Additionally, the AUC values in BRCA, LUSC, COAD, STAD, READ, KIRC, and HNSC were all between 0.8 and 0.9, indicating that NUP62 also has diagnostic value for these tumors. Figures 3B, C depict the ROC curves of NUP62 in ESCA and BRCA, respectively. Therefore, NUP62 is an effective diagnostic biomarker across various cancer types.

A

D

-

-

Normal

Tumor

MUPER (TOGA)

0

DE

C

Zscore

1

0

-

8

9-

6

2

1

0

AC

Ccing

1

4

AR

PA

CARE

100

*

POHO

READ

Ag

0

B

Cancer Tyow

-1

NÃO

DE

a

as

2

U

L

V

=

14

Normal

Tumor

5

£

A

4

4

-

5

.

1

1

5

BEIM

A

4

PRAD_BARNIN

2

Cance Type

C

-

Zscore

-

-

,

1

0

-1

L

-

-

-

-

THÍCH

E

5-

A279_T2800H

# patients

+

1-

Nspt_C

0

100

200

300

400

527mm

PTM (dbP’TM)

Phosphorylation

Acetylation

Ubiquitination

Tapology

F

9

G

Altered in 6051 (100%) of 6061 samples.

A

4

!

N

ni-

NA

KRAS

13%

For

NOTÇAT

B

DOFR

7%

%

294

HA

8%

%

HANTI

%

Cancer

OT . TRA

G

Missionsa Muzarion 1

Barwation_Start_840

CHOL

PAD

UVM

H

BLCA CESC

1 Mười_HỆ

ESCA

PRADO READ

STAD TGOT

THICA THƠM

KIRG

LING

UCEC

NUP62

5

7

8

8

1

1

4

1

1

1

2

3

7

2

1

1

10

9

22

1

NOTCH-Signaling: NOTCH1

13

10

20

1

20

15

2

89

1

3

5

39

6

25

39

3

2

3

2

5

38

30

2

2

1

54

2

TP53-Signaling: TP53

17

201

350

24

4

223

5

169

127

358

21

12

7

9

268

110

290

406

14

373

109

1

59

108

87

72

214

1

3

4

204

52

WNT-Signaling: WNT1

1

2

1

10

2

2

1

3

7

1

1

1

1

1

8

4

RTK-RAS-Signaling: MAPK1

1

2

2

11

3

1

3

9

1

1

2

3

3

4

1

1

1

5

5

9

Mutation Frequency (%)

RTK-RAS-Signaling: EGFR

2

7

15

7

11

5

95

14

2

1

36

6

72

12

6

1

3

3

4

34

21

1

36

75

RTK-RAS-Signaling: BRAF

1

13

8

3

1

58

1

2

8

9

1

4

3

42

18

1

2

2

1

7

5

1

245

14

296

25

1

50

RTK-RAS-Signaling: KRAS

16

6

16

2

174

2

2

2

1

2

5

6

2

5

168

8

1

6

119

2

60

2

13

40

13

4

1

102

7

25

TGF-Beta: TGFBR1

8

5

2

14

2

1

2

1

1

1

7

4

1

4

2

4

2

6

12

25

a

PI3K-Signaling: PIK3CA

1

90

352

84

2

129

17

39

92

6

3

45

13

28

56

2

8

5

11

22

7

22

71

3

4

266

20

MYC-Signaling: MYC-

3

3

3

9

3

1

6

2

1

4

4

1

2

1

1

1

11

7

20

2

Hippo-Signaling: YAP1

6

1

3

4

3

1

2

1

1

3

5

2

1

1

10

3

17

1

CellCycle-Signaling: CDKN1A

39

2

4

2

1

1

2

1

1

7

2

3

1

2

1

2

10

S

1

10

ACC (n=92)

BLCA (n=411)-

BRCA (n=1026)

CESC (n=291)

CHOL (n=36)

COAD (n=407)

DLBC (n=37)-

ESCA (n=185)

GBM (n=403)

HNSC (n=509)

KICH (n=66)-

KIRC (n=370)

KIRP (n=282)

LAML (n=140)

LGG (n=526)

LIHC (n=365)

LUAD (n=567)

LUSC (n=485)-

MESO (n=82)

OV (n=412)

PAAD (n=178)-

PCPG (n=184)

PRAD (n=498)

READ (n=149)

SARC (n=239)

SKCM (n=488)

STAD (n=439)

TGCT (n=151)-

THCA (n=500)

THYM (n=123)-

UCEC (n=531)

UCS (n=57)-

UVM (n=80)

FIGURE 1 Expression and Mutation of NUP62. (A) Expression of NUP62 mRNA in tumor and normal tissues. (B) Expression of NUP62 in tumor and normal tissues from the CPTAC database. (C) Expression of NUP62 in tumor and paired adjacent normal tissues. (D) Expression and distribution of NUP62 in organs of tumor and normal tissues. (E) Mutation sites of NUP62. (F) Mutation frequencies and corresponding mutations of NUP62 in different cancers. (G) Gene alteration sites and numbers of NUP62 in different cancers. (H) Distribution of NUP62 and other signals in cancers. * P<0.05. ** P<0.01. *** P<0.001. **** P<0.0001.

3.5 Correlation between NUP62 expression and pan-cancer prognosis

To gain deeper insights into the clinical significance of NUP62 in the field of cancer, we analyzed its prognostic value in multiple malignant tumors. The results in Figure 4A clearly show that NUP62 expression levels are significantly associated with poor prognosis in various cancers, serving as a risk factor for their

prognosis. Figures 4B-M present the KM survival curves for overall survival (OS) based on NUP62 expression levels in some cancers. The results indicate that high expression of NUP62 is significantly associated with shorter OS in KIRP, KIRC, ACC, LGG, LIHC, mesothelioma (MESO), skin cutaneous melanoma (SKCM), and SARC (p < 0.05). However, in uveal melanoma (UVM), BLCA, STAD, and CESC, high expression of NUP62 is associated with longer survival, potentially playing a protective role.

FIGURE 2 Relationship between NUP62 Expression and Clinical Characteristics Across Various Cancer Types. (A) Correlation between gender and NUP62 expression in pan-cancer. (B) Correlation between NUP62 expression and T stage in pan-cancer. (C) Correlation between NUP62 expression and N stage in pan-cancer. (D) Correlation between NUP62 expression and M stage in pan-cancer. (E) Correlation between NUP62 expression and G stage in pan-cancer. (F) Correlation between NUP62 expression and tumor stage in pan-cancer. (G) Correlation between age and NUP62 expression. * P<0.05. ** P<0.01. *** P<0.001. **** P<0.0001.

A

B

MALE FEMALE

1 8

n

F

T4

GEM(FEMALE-54,MALE-98 GBMICHXFEMALE-2N MALE-3813

LOFAFEMALE-225 MALE-288

LUAIXFEMALE-27%,MALE-237)

COAIXFEMALE-170MALE-156

COADREAIXFEMALE-172MALE-298

LAMLEFEMALE-KO,MALE-93 BRCA(FEMALE-182) MALE-12)

ESCA(FEMALE-26,MALE -155

STESIFEMALE-17ZABALE-423

BARCEFEMALE-14 AALE-117

KERRFEMALE-75,ALALE-213)

KIPANEFEMALE-2X,MALE-SHE

HENSCEFEMALE-186,MALE-3KP STAD(FEMALE. HA.MALE-36K

KERCEFEMALE-186,MALE-346)

LUSCIFEMALE-129,MALE-36)

THYM(FEMALE-STMALE-63

LAICIFEMALE-SMO MALE-20

THOAIFEMALE-SOR MALE-1MG)

MESO(FEMALE=16,ALALE-71) READ(FEMALENE) MALE-49

PAADUJEMALE-SOMALE IREAFEMALE.100 MALE-77 SKCM(FEMALE-42 MALE-60

UVMEFEMALE-35 MALE-46

BECAIFEMALE-106,MALE-381)

ACCEFEMALE-40,MALE-31)

KITEFEMALE JIMALE-39

CHOLEFEMALE-20,MALE-16)

DLECEFEMALE-23,MALE-22)

CESCYTE-140.72-72.73-36,76-35) LUANTI-10AT2-2MTI),14-18

CORDREAIXYI=1&T2-57,73-259,74-56) COATXTI-6.72-44,73-196,74048)

BRCA(TP-280,72-630,73-198,74-48)

ESCA(TI-31.12-43,23-97,76-6)

STES(TI=62,72-132.73-299,74-121)

KENNT1-199.77-11,13-60)

RIPAN(Y)=484,12-127,13-257,76=16)

STAINTI=21,T2-89,13=180,74-115)

PRADET2-188,73-290,74-10)

HNSC(11=11,72-189,73-120.74-489Y

KERCYT1-272.12-0417-179,76-0)

LUISCETI-112,12-292,T3-71,74-253

LICYTI-179.12-94.73-86,74-13)

THCAKTI-142,72-166,73-171,74-25)

MESKTI-14.T2-26.79-32.74-13

READY(T]=4,22-11,73-6),76-56)

PAAD(T)=7,12-24,13=142,Ta-))

TOCHT-X17-KT3-6)

SKCM(T2-5,13-10,74-84)

UVART2-14.37-32,74-33)

BECA(TI-XT2-150,73-195,74-58

ACCITHRTAGIT3KY4.17Y

KICH(TP-21,73-25,13-1)

CHOLTHIRTY-12.79-5

C

D

5%

10

*

MI

Expression

2

0

CESCENE-133.N16

LUAD(NO 329,N1-96,N2-24)

CORINND-166,N1-72.83-49; COADREADY(NO-204,N=100,N2-76

BRCAPR-SIENI-MALNI-HIRNI-TH

ESCAING-76,31-78,N2-13,N3-7

STES(-200,NI-199,N2-91,N3-89)

KKAMO~144Nl-25.83-4

KPPANCO-GUN-3N24)

STARING-12ANI-111,22-78.83-82)

PRAIXNO-343.NI-19;

IENSCOMO-228.NI-TAN2-205,N3-133

KCXNORNt.19

LUSKYNO-317N1-136X2-10N3-

LEKING-25ON-6

TECAONG-179,31-225)

NEESKKNU-44301-802-2683-8

READ(N-38.NI-23N2-21)

PARDENO-49,NI-120)

TGE TỪNG-19NI-22NHAN3-6

SKEMING-58281-8.N2-10,N9-10)

BECA(NG-31721-46,22-74,N3-8

ACLINO MANG

KKHENGER.NJ-3)

CHOLENO-2ANI-S

2

CESKYMO-116,MI-16

LUJADEMO-344,M1-25%

COAX380-193.M:+46

COADREADOR-256.Mi-533

ESCAOP-14MI-19 BRCA(MO-907341-22Y

STESOMO 507,385-43)

KERRMP-286.1410123

KIPAN(MD-6HK,A1] -43)

STAD(M9-366,MI -27Y

PRAD(MO-453385-3)

BNESKIND-SE2.MI-S

KERCINO-4UMI-NY

LUISCI380-408.33-7

LSKEM9-264M6-8

THECA(MO-281,M3-9)

MESOME-SEMI-3

READ(MO-63,381-12)

PAAD(MO-80 M3-4)

TOCT-IMMIG

BKCMEM9-07M3-3)

UVMEM9-123-4

BECAEMO-196,383-11)

ACCOR seMi-19

CHOL30-28345-55

E

G

UVM(N=79)

pValue

KIRP(N=285)

-0.00

PCPG/N=17%

1

BOA(N=181)

0.02

BRCA(N=1090)

DLBC(N=47

-0.04

PAAD(N=178)

LIHC(N=368)

-0.06

READ(N=91)

0.08

1

BLCACAO

GIMLOGHI2-247,G3-2663

LOGKG2-347,G3-2667

CESCHA-18,G3-125,03-1186

ESCAFGI-1863-74,03-49;

STES(GE-30,07-2223-24)

KIPAN(3 -14/2-228/63-296,04-74)

SYAD(GI-12.62-188,03-245)

UCH101-14022-21,63-1413

HENSC(G1-61.822-308323-124,66-7/

KIRDIGI-14,202-228,01-200,06-76

LZ一公分×328-1204-1

PAD(C1-31/32-95,63-65)

09(3-47,63-30)

CHOLARS-15.03-18

THCA(N=504)

0.10

HA~ 517)-

THYM(N=118)

STES(N-590)

SARC(N-258)

COADREAD(N-377)

OV(N-419)

F

COADO-286

UCECON-177

Group

PRADIN-4OS

Stage II

LAML(N-173)-

Stige 1

KIRC(N-530)

T

Stige EN Stage TV

MESO(N-87)

LUAD(N-494)

CESCYN-304)

KIPAN(N-881)

KEEP(Stage |-177,3-25331-32.39-16) KEYNESough Two4.8-1993-18570653

THYM( Stage 1:3631-61311-1439-6) LUISCEStage1-242.8-161.31-8419-7)

READ(Stage 1-1231-2411-3339-13) PAIXStage 1-21.3-14738-3-4

GEM(N-152)

CESCEStage 1:162,31-6RJ-45JV-21)

LEAD(Stage 1-274 3-122.18-83.TV-26)

COAD(Stage1-448-118311-8239-40

CONIREADISIng toNell=134.81-11839-33)

BRCA(Stage 1-182.38-617.311-248.19-20)

BC入Sage下一共用的,并一次8-6)

STES(Stage 1-76. 8-201.18-230.19-$7)

STADEStegn 1-88,31-121.311-169./V-41)

UCDCESTAK PK.#-NE-V-19)

HNSC(Stage 1-2131-32.81-93.3-346)

KERKCeStage 1-266 3-57,JE-123.1V-81)

LIC(Stage 1-169.1-SEJR-8SIV-SY

THCA(Stage1-2x3.#~$2.38-112TY-55)

ME.SOEStage PpHetH isH 46Y ES)

OV(Stage #~24.38-328.19-63)

TOACT(Stage 1-164.8-13.18-14)

SKCSH Stage II-66.311-26.39-3)

USMrttage #~39.38-36.TY-4)

UCS(Stage 7-22.8-S.18-20.JY-19

Bi.CA:Stage E=130.11-100.19/-133)

XCC(Stage 1-9.31-363-18JV-15

KICHStage 1-21,8-25.18-14.19-6

CHOL/Sugs 1-19.3-91-1)

ĐLBC(Sage P-831-4638-53V-12)

LGG(N-508)

STAD(N-409)

LUSC(N-489)

TGCT(N-132)

CHOL(N-36)

GBMLGG(N-660)

SKCM(N-102)

KICH(N-66)

UCS(N-57)

-0.2

-0.1

0.0

0.1

0.2

Correlation coefficient( pearson)

TCGA

A

DataBase

TCGA-GTEx

B

ESCA

ESCA

0

0

UCEC

0.62

0.831

VA

0.8

THCA

0.547

0.685

STAD

0.947

0.861

Sensitivity

Sensitivity

0.6

READ

0.866

0.943

AUC: 0.978

95% CI: 0.942-0.998

AUC: 0.934

95% CI: 0.910-0.956

PRAD

0.4

0.617

0.684

PCPG

0.522

0.701

8

8

PAAD

0.968

0.673

8

LUSC

3

0.894

0.941

Cancer Type

0.0

0.2

0.4

0.6

0.8

1.0

LUAD

0.778

0.0

0.6

0.894

0.2

0.4

0.8

1.0

1 - Specificity

1 - Specificity

LIHC

0.667

0.788

C

BRCA

BRCA

KIRP

0.677

0.787

0

2

KIRC

0.848

0.965

0

0.8

KICH -

0.605

0.633

HNSC

0.876

0.861

Sensitivity

0.6

Sensitivity

0.6

ESCA

0.935

0.981

AUC: 0.919

95% CI: 0.895-0.940

AUC: 0.895

0,4

0.4

95% CI: 0.877-0.911

COAD

0.895

0.881

CHOL

1

à

0.2

BRCA-

0.895

0.9

919

0.0

BLCA

0.764

0.858

0.0

0.2

0.4

0.6

0.8

1.0

0.0

0.2

0.4

0.6

0.8

1.0

-1.0

-0.5

0.0

0.5

1.0

1 - Specificity

1 - Specificity

AUC

FIGURE 3 Correlation between NUP62 Expression and Diagnosis of Various Cancer Types. (A) Diagnostic value of NUP62 across various cancer types; (B, C) ROC curves of NUP62 in ESCA and BRCA.

FIGURE 4 Survival Analysis Spectrum of NUP62 Across Multiple Cancer Types. (A) Correlation between NUP62 expression levels and OS, disease-specific survival (DSS), disease-free interval (DFI), and progression-free interval (PFI). (B-M) Relationship between NUP62 expression levels and OS in patients with different cancer types.

A

Method

B

NUP62 Survival Curve in KIRP(p<0.001)

C NUP62 Survival Curve in KIRC(p<0.001)

Survival Type

ACC

1.0

0

BLCA

BRCA

0.8

0.8

CESC

CHOL

Overall Survival

0.6

Overall Survival

COAD

0.6

DLBC

ESCA

0.4

0.4

GBM

HNSC

0.2

0.2

KICH

High (n= 95)

High (n= 174)

KIRC

0.0

Low (n= 193)

0.0

Low (n= 358)

KIRP

LGG

Prognostic role Method

0

5

10

15

0

2

4

6

8

10

12

LAML

Risky

Cox

LIHC

Protective N/A

Log-rank

Time (year)

Time (year)

LUAD

Nonsense

Survival Type

LUSC

DFI

D

NUP62 Survival Curve in ACC(p<0.001)

E NUP62 Survival Curve in LGG(p<0.001)

MESO

DSS

OV

OS

9

PFI

1.0

PAAD

PCPG

0.8

0.8

PRAD

READ

Overall Survival

0.6

Overall Survival

0.6

SARC

SKCM

0.4

STAD

0.4

TGCT

THCA

0.2

0.2

THYM

High (n= 26)

High (n= 357)

UCEC

0.0

Low (n= 53)

0.0

Low (n= 168)

UCS

0

2

4

6

8

10

12

0

5

10

15

UVM

Time (year)

Time (year)

F

NUP62 Survival Curve in LIHC(p<0.001)

G

NUP62 Survival Curve in MESO(p<0.001)

H

NUP62 Survival Curve in SKCM(p=0.004)

I NUP62 Survival Curve in SARC(p=0.004)

0

0

1.0

0

0.8

0.8

0.8

0.8

Overall Survival

0.6

Overall Survival

0.6

Overall Survival

0.6

Overall Survival

0.6

0.4

0.4

0.4

0.4

0.2

0.2

0.2

0.2

High (n= 111)

High (n= 59)

0.0

Low (n= 256)

0.0

Low (n= 26)

High (n= 177)

High (n= 156)

0.0

Low (n= 279)

0.0

Low (n= 107)

0

2

4

6

8

10

0

2

4

6

0

5

10

15

20

25

30

0

5

10

15

Time (year)

Time (year)

Time (year)

Time (year)

J NUP62 Survival Curve in UVM(p=0.006)

K NUP62 Survival Curve in BLCA(p=0.022)

L NUP62 Survival Curve in STAD(p=0.002)

M NUP62 Survival Curve in CESC(p=0.026)

0

1.0

1.0

0

0.8

0.8

0.8

0.8

Overall Survival

Overall Survival

0.6

Overall Survival

0.6

Overall Survival

6

0.6

0.4

0.4

0.4

0.4

0.2

0.2

0.2

0.2

High (n= 33)

High (n= 257)

High (n= 180)

0.0

Low (n= 47)

0.0

Low (n= 149)

High (n= 233)

0.0

Low (n= 156)

0.0

Low (n= 113)

0

1

2

3

4

5

6

7

0

2

4

6

8

10

12

14

0

2

4

6

8

10

0

5

10

15

Time (year)

Time (year)

Time (year)

Time (year)

3.6 The role of NUP62 in cancer pathways

To delve deeper into the potential roles and specific functions of NUP62 in cancer, we employed an integrated analysis of signature gene expression to assess the activity status of NUP62 pathways. We utilized the z-score parameter from Gene Set Variation Analysis (GSVA) to quantify and evaluate 14 functional status gene sets, covering angiogenesis, apoptosis, cell cycle regulation, cell differentiation, DNA damage response, DNA repair mechanisms,

epithelial-mesenchymal transition (EMT), hypoxic adaptation, inflammatory response, invasive behavior, metastatic potential, cell proliferation, cellular quiescence, and cellular stemness features. From this, composite z-scores were obtained. Subsequently, we calculated the Pearson correlation coefficients between NUP62 and each functional status gene set score (Figure 5A). The results revealed significant positive correlations between NUP62 expression levels and cell cycle regulation, DNA damage response, and DNA repair mechanisms. Additionally, we

A

Angiogenesis

Apoptosis

Cell Cycle

Differentiation

DNA damage

DNA repair

EMT

4

2

Oncogenic combined z-scores (zscore)

0-

-2

-4

R = 0.064; p = 6.5e-11

R=0.21.p<2.26-16

R=0.37.p<2.2e-16

R=0.16.p <2.2e-16

R=0.4. p<2:28-16

R=0.32, p<2.2e-16

R=0.12, p <2.2e-16

Hypoxia

Inflammation

Invasion

Metastasis

Proliferation

Quiescence

Stemness

4-

2

0-

-2-

-4

R=0.015, p = 0.13

R = 0.11, p <2.2e-16

R=0.11, p <2.2e-16

R=0.12, p<2.2e-16

R=0.31,p <2.2e-16

R =0.029, p = 0.0031

R=0.1,p < 2.2e-16

-2

0

2

-2

0

2

-2

0

2

-2

0

2

-2

0

2

-2

0

2

-2

0

2

B

NUP62 (zscore)

KEGG: Vitamin B6 Metabolism

Leucine And Isoleucine Degradation

-KEGG: Valine Leucine And Isoleucine Biosynthesis

KEGG: Ubiquinone And Other Terpenoid-Quinone Biosynthesis

KEES: Tyros ne Metabo

IREGG: Themine Metabolism

İbase: Tamne Backbone Biosynthesis

[BECS aPerdu Busymines

-KEGG: Synthesis And Degradation Of Ketone Bodies

-KEGG: Sulfur Metabolism

KEGG: Steroid Hormone Biosynthesis

-KEGG: Biosynthesis

-KEGG: Starch And Sucrose Metabolism

KEGG: Sphingoapie werd bonm

BER: Bochumh Metcolem RECO: Patate Metabolism

-KEGG: Pyrimidine Metabolism

-KEGG: Purine Metabolism

KEGG: Propanoate Metabolism

-KEGG: Primary Bile Acid Biosynthesis

KEGG: Porphyrin And Chlorophyll Metabolism

Phosphonsle And Phosphinale Metabolism

KEGG: Phenylalanine Tyrosine And Tryptophan Biosynthesis

IREGG: Pentose Phosphate Pathway

-KEGG: Pantose And Glucuronale

KEGG: Pantothenate And Coa Biosynthesis

REGGIOwat Phosphorylation

-KEGG: Other Types Of O-Glycan Biosynthesis

-KEGG: Other Glycan Degradation

KEGG: One Carbon Pool By Folate

-KEGG: BEGG: Nicotinate And Nicotinamide Metabolism

KEGG: Neomycin Kanamycin And Gentamicin Biosynthesis

Mucin Type O-Glycan Biosynthesis

-KEGG: Metabolism Of Xenobiotics By Cytochrome P450

KEGG: Mannose Type O-Glycan Biosynthesis Lysine Degradation

KEGG: Lipoic And Metabolism

-KEGG: Linoleic Acid Metabolism

-KEGG: Glyoxylate And Dicarboxylate Metabolism

BEEG: Glycosylphosphatidylinositol Gplanchor Biosynthesis

KEGG: Glycosphingolipid Biosynthesis Globo And Isoglobo Series

-KEGG: Glycosphingolipid Biosynthesis Ganglio Series

-KEGG: Glycosaminoglycan Degradation

Keratan

-KEGG: Glycosaminoglycan Biosynthesis Heparan Sulfate Heparin

-KEGG: Glycosaminoglycan Biosynthesis Chondroitin Sulfate Dermatan Sulfate

KEGG Glycolysis Gluceneogenesis RECO Che Serine And Threonine Metabolism

KEGG: Glycerollipid Metabolism

NES

-KEGG: Glutathione KEGG: Galactose Metabolism

-KEGG: Fructose And Mannose Metabolism KEGG: Folate Biosynthesis

2

KEGG: Fatty And E ongation

0

KEGG: Fatty And Degradation

BEGGE

LKEGG: Drug Metabolism Other Enzymes

-2

three: ong Cytochrome Paso

IKEGG: D-Glutamine And D-Glutamate Metabolism

-KEGG: D-Arginine And D-Ornithine Metabolism

-log10[p.adjust)

KEGG: Cysteine And Methionine Metabolism

TREE’S: Corale Cycle Tea Cycle

· 0.0

these: gumane maraton

TRERA: ON Machen Biotin Metabolism

25 5.0

KEGG: Biosynthesis Of Unsaturated Fatly Acids

Rua Alanine Metabolism

7.5

-KEGG: Ascorbate And Aldarate Metabolism

Biosynthesis

-KEGG: Arginine And Proline Metabolism

-KEGG: Arachidonic Acid Metabolism

KEGG: Amino Sugar And Nucleotide Sugar Metabolism

-KEGG: Alanine Aspartate And Glutamate Metabolism

habeYanoheng Metabolism

alles Wat Bela Casin Sanalin Hallmark vit Beta Canenin Signaling

Hallmark Uv Response Dn

- Hallmark Unfolded Protein Response

-Hallmarkc: Tnfa Signaling Via Nfkb

-Hallmark Tgf Beta Signaling

Hallmark Spermatogenesis

Hallmark Reactive Oxygen Species Pathway

Hallmark Piak Sign

Hallmaric Pisk Akt Mtor Signaling

Hallmark Pancreas Beta Cells

Hallmark P53 Pathway

Hallmark Oxidative Phosphorylation

-Hallmark Notch Signaling

Hallmaric Myogenesis

Hallmark Myc Targets V2

Hallmark Myc Targets Vi

-Hallmarkc Kras Signaling Up

Hallmarkc De

Hallmark Interferon Gamma Response

Hallmark Interferon Alpha Response

Tallene Inflammatory Response

-Hallmark 115 Jak Stat3 Signaling

-Hallmarkc 112 Stat5 Signaling

-Hallmark Hypoxi8

Hate

Hallmark Herre Metabolism

Hallmaric Hedgehog Signaling

Hallmark Gam Checkpoint

-Hallmark Fatty Acid Metabolism

Hallmark Estrogen Response Late

-Hallmark Estrogen Response Early

☒ ☒

Hallmark Epithelial Mesenchymal Transition

-Hallmanic E2f Targets

Hallmanc Dna Repair

Coagulation

Hallmark Cholesterol Homeostasis

-Hallmark Bile Acid Metabolism

Hallmarkc Apoptosis

Hallmark Apical Surface

-Hallmark Apical Junction

Hallmarkc Angiogenesis

Hallmanc Androgen Response

-Hallmark Allograft Rejection

Hallmark Adipogenesis

ACC

BLCA

BRCA

CESC

CHOL

COAD

DLBC

ESCA

GBM

HNSC

KICH

KIRC

KIRP

LAML

LGG

LIHC

LUAD

LUSC

MESO

OV

PAAD

PCPG

PRAD

READ

SARC

SKCM

STAD

TGCT

THCA

THYM

UCEC

UCS

UVM

FIGURE 5 (A) Relationships between NUP62 and 14 malignant features of cancer; (B) Enrichment differences of NUP62 in 50 HALLMARK and 83 metabolic gene sets.

collected datasets from cancer patients and conducted Gene Set Enrichment Analysis (GSEA). The analysis results indicated that NUP62 may be involved in various tumor-related pathways and metabolic regulatory processes, such as the G2/M checkpoint, E2F transcription factor target genes, allograft rejection, and mitotic spindle formation (Figure 5B).

3.7 The relationship between NUP62 and cellular pathways

Utilizing the TCGA database, we conducted an in-depth exploration of the interaction between NUP62 and functional proteins within the TCGA database, and found that NUP62

exhibits significant correlations with multiple key functional proteins across various cancer types. Specifically, in UVM, the expression level of NUP62 is significantly positively correlated with CASPASE8, CYCLIN D1, CRAF (pS338), AKT, and RB proteins, while it is significantly negatively correlated with the expression levels of ACC- pS79, ACC1, CKIT, ATM, and PKCơ. In ACC, the expression of NUP62 is positively correlated with the expression of CABL protein and negatively correlated with the expression of FOXO3A protein (Figure 6A). Figure 6B illustrates the relationships between cancer and 10 cancer-related pathways, where the activation of apoptosis, activation, cell cycle, and DNA damage promotes the occurrence of cancer. Subsequently, we detailed the most significantly correlated functional proteins of NUP62 in ACC and UVM from the TCPA database (Figures 6C, D).

3.8 The relationship between NUP62 and immune subtypes and immunotherapy

Based on nearly 10,000 cancer samples from the TCGA database, we subclassified them into six distinct immune subtypes. The analysis revealed that C1, C2, C3, and C4 subtypes dominate among all immune subtypes. Additionally, in cancer samples with high NUP62

expression, the proportion of the C2 subtype was significantly higher than that in samples with low NUP62 expression; conversely, in patients with low NUP62 expression, the proportions of C3 and C5 subtypes were significantly higher than those in patients with high NUP62 expression (Figure 7A). Subsequently, we delved into the correlation between NUP62 expression levels and immunotherapy response in clinical trials. Among the two immunotherapy cohorts included in the study, patients with high NUP62 expression had a significantly higher response rate to immunotherapy than those with low expression (Figure 7B). Notably, we evaluated the potential efficacy of NUP62 as a predictor of immunotherapy response using ROC curve analysis (Figure 7C). By comparing the impact of two different immunotherapy regimens on patients’ overall survival, we found that patients with elevated NUP62 expression had a significantly longer overall survival compared to those with decreased NUP62 expression (Figure 7D).

3.9 The relationship between NUP62 and immunity

Immunotherapy has emerged as a primary modality in cancer treatment, prompting our investigation into the potential association

FIGURE 6 (A) Top five functionally related proteins of NUP62 across various cancers in the TCGA database; (B) Relationships between cancer and 10 cancer- related pathways; (C, D) NUP62-related functional proteins in ACC and UVM from the TCPA database.

A

5

Y

2

ACC

class

BRCA

CHOL

ACC

AM

BRCA

DLBC

M

ESCA

HIER

GBM

Ferd

HNSC

ESCA

KICH

GBM

HNSC

KIRC

KICH

KIRP

KIRC KIRP

GG

LIHIG

LGG

m

LUAD

LIHC

MESO

LUAD

PCPG

MESO

READ

PCPG

.

READ

SARC

W

STAD

SARC

STAD

1997

THCA THYM

POST

THCA THYM

UCEC

UCEC

UCS UVM

UCS

UVM

class

CASPASER

CYCLIND

CRAF_PS359

AKT

RB

ACC_p579

ACC1

CKIT

ATM

PKCALPHA

STATSALPHA

CMET

FOXO3A

BC

CD31

PpGa

ACETYLATUBULINLYS40

PKCDELTA_PS664

GEPD

P38MAPK

DUSPA

PASIN

P62LCKLIGAND

CHK1

SHC_PY317

PCADHERIN

YB1_PS102

E

PERO

EPPK1

PI3KP85

FOXO3A_pS31BS321

CIUN PS73

SRC_PY416

INPP4B

SCD1

BOK

PKCPANBETAll_pS660

FIBRONECTIN

EIF4E

MEK

CYCLINE1

FOXM

X4EBPI

EEF2K

YAP_p$127

ARIDIA

CASPASES

YB

IGF1R_pY1135Y1138

RAB11

CRAP

HER3

NAPSINA

CLAUDIN

XRCC1

SRC_PY527

ERCCT

PCNA

NF2

TeZ

PKCALPHA_PS657

ANNEXINVI

DIRASS

PREX1

RAPTOR

CHK2

PI3KP110ALPHA

HER3_PY 1268

CASPASE7CLEAVEDD198

VEGFR2

ERK

BRD

GSK3ALPHABETA

BAX

SYK

YAP

ECADHERIN

BETACATENIN

Tomate

WE

LYON

CABL

RBM15

ERALPHA

RAK

ARAF_PS291

EGFR_PT1065

P38_PT1801162

STAT3_PY705

MYH11

RICTOR

CYCLING!

COK1PM

CDN

LCK

RB_058078811

Correlation

0.6

0.2

-02

-0.6 1

B

C

Pathway

0.50

0.8

R =- 0.6, p= 1.8e-05

ACC

BLCA

R=0.56, p =6.6e-06

BRCA

0.25

0.6

CESC

CHOL

COAD

CABL (ACC)

FOXOBA (ACC)

0.00

DLBC

0.3

ESCA

GBM

-0.25

HNSC

0.0

KICH

KIRC

Pathway

Pathway

-0.50

KIRP

LGG

Activate

Apoptosis

Inhibit

CellCycle

-0.3

LIHC

Nonsense

DNADamage

D

10

11

12

10

11

12

LUAD

EMT

0.8

LUSC

Hormone_AR

MESO

Hormone_ER

PI3KAKT

0.5

R=

-0.92, p < 2.20-16

OV

RASMAPK RTK

PAAD

0.4

PCPG

TSCmTOR

R=0.8. p=0.0027

PRAD

CYCLIND1 (UVM)

ACC_p879 (UVM)

0.0

READ

SARC

0.0

SKCM

05

STAD

TGCT

THCA

-0.4

-10-

THYM

UCEC

UCS

UVM

10.4

10.6

10.8

NUP62

11.0

11.2

10.4

10.6

10.8

11.0

NUP62

11.2

9079 TCGA patients

A

C1

C2

C3

C4

C5

C6

Expression groupssubtype group (n=9079)
C1 (n=2403, 26%)C2 (n=2580, 28%)C3 (n=2384, 26%)C4 (n=1147, 13%)C5 (n=385, 4%)C6 (n=180, 2%)P-value
Expression-Low (n=4540)1181(26%)1138(25%)1272(28%)624(14%)235(5%)90(2%)
Expression-High (n=4539)1222(27%)1442(32%)1112(24%)523(12%)150(3%)90(2%)p<0.001
FIGURE 7 (A) Differences in the number of pan-cancer immune subtypes between the high and low NUP62 expression groups. (B) Expression of NUP62 in responders and non-responders within different immunotherapy cohorts. (C) Diagnostic value of NUP62 in different immunotherapy settings. (D) Correlation between NUP62 expression and survival rates in different immunotherapy cohorts.

Expression-low Expression-high

B

Wolf cohort 2021 Anti-PD-L1

C

Ascierto cohort 2016 (Anti-PD-1)

D

Wilcoxon, p = 0.0015

4 responders and 7 nonresponders

Nathanson cohort 2017 (Anti-CTLA-4)

1.00

2

1.00

NUP62 Expression (z-score)

0.75

0.75

Sensitivity

Overall survival

0.50

-

0

0.50

0.25

Log-rank

p = 0.023

0.00

0.25

0

2

4

6

-2

Number at risk

AUC = 0.964

NUP62

High

9

2

1

0.00

Low

0.25

0.50

0.75

15

1.00

-7

0

0.00

2

NR

R

1-Specificity

0

Time in years

2

4

6

Response

Ascierto cohort 2016 Anti-PD-1

Kim cohort 2019 (Anti-PD-1/PD-L1)

T-test, p = 0.0052

8 responders and 19 nonresponders

Hugo cohort 2016 (Anti-PD-1)

1.00

1.00

1

0.75

NUP62 Expression (z-score)

Sensitivity

Overall survival

0.75

0.50

·

0

0.50

0.25

Log-rank

p = 0.038

0.00

0

1

2

3

0.25

Number at risk

-1

AUC = 0.704

NUP62

High

21

15

6

0

0.00

Low

6

2

7

0

0.00

0.25

0.50

1-Specificity

0.75

1.00

NR

Response

R

0

Time in years

1

2

3

between NUP62 and cancer immunity. Our results indicate that across various cancer types, the expression level of NUP62 significantly correlates with the expression of immune-related genes, such as those in THCA, KIRC, KIRP, and PCPG (Figure 8A). Furthermore, by examining the levels of immune cell infiltration in the TME, we found that the expression of NUP62 is significantly correlated with the infiltration levels of multiple immune cell types, including Th1 cells, Th2 cells, dendritic cells (DC-TIMER),

M0 and M1 macrophages (Macrophages-M0-CIBERSORT), neutrophils (neutrophil-TIMER), T cells (T-cells-MCPcounter), and regulatory T cells (T-cells-regulatory-CIBERSORT), among others, across several cancer types (Figure 8B). These findings suggest that NUP62 may play a crucial role in the tumor immune microenvironment, potentially affecting tumor growth and the efficacy of immunotherapy by regulating the expression of immune cells or immune-modulating genes.

FIGURE 8

A

-

B

B_col_TIMER

-B_cells_memory_CIBERSORT

B_bffs_memory_CIBERSORT_ABS

B_cols_name CIBERSORT

_B_cells_naive_CIBERSORT_ABS

LB cels_quaniseq

_ B_lineage_MCPcounter

PR_EPIC Claus-switched_memory_B-cells_xCel Memory_U-cells_xCell naive_B-cells_xCel

-CD4_Tcells_EPIC

CD4T-T-cells_xcel

_CD4+_TO_xCell

CD4+ Tem xCell

T_cel_CD4_TIMER

T_cells_CD4_memory_activated_CIBERSORT

cats CD4 memory activated_CIBERSORT ABS

T cells CD4 memory resting CIBERSORT

T cells_CD4 memory_resting CIBERSORT ABS

[_c_COM MINA_CIBERSORT_ABS

T_cells_Cl_quantiseq

T_cels_folicular_helper_CIBERSORT

cas_folicular_helper_CIBERSORT_ABS

Thi_cells_xCel Th2_cells_xCel

.CD8_T_cells_MCPcourir

CDS_Tcells_EPIC

“LCDB+_naive_T-cells_xCel

L CDS+ T-cells xcel

.CDB+_Tom_xCell

COBe Tem SCell

Trel COB TIMER

POR CIBERSORT ABS T_cat_OD_MBERSORT_ABS T_cels_CD8_quantiseq

aDC_xCel DC TIMER

.DC_xCel

.Dendrise_cols_Activated_CIBERSORT

Dendrido_cells_activaled_CIBERSORT_ABS

-Deromac teis Nuantiseq

.Dendride calls festing_CIBERSORT

Dendrido_cels_resting_CIBERSORT_ABS

-Myeloid_cencritic_cells_MCPcounter .PDC_xCel

.Endotheist_cells_MCPcounter Endothelal_cells_xCell Endothelal_EPIC

Eosinophils_CIBERSORT

Eosinophils CIBERSORT ABS

Eosinophils_xCell

+44L CAFS EPIC

Fibroblasts_MCPcounter Fibroblasts xCell

.Macrophages_MO_CIBERSORT Macrophages_MO_CIBERSORT_ABS

Macrophages_MI_CIBERSORT Macrophages_MI_CIBERSORT_ADS.
Macrophages_M_quantiseq
Macrophages_M1_xCell
Macrophages_M2_CIBERSORT
Macrophages_M2_CIBERSORT_ABS-
Macrophages_Na_quantisis-
Macrophages_M2_xCell_
Macrophage_TIMER
Macrophages EPIC. Macrophages_xCell,
Mast_cels_activaled_CIBERSORT
Mast cells activaled_CIBERSORT_ABS_
Mast_cells_resting_CIBERSORT
Want sulle pasting CHEBSORT ARE Mast_cells_resting_CdEnt-d
Mast_cells_xCell
Monocytes_CIBERSORT
Monocytes_CIBERSORT_ABS.
Monocytes quantisco Monocytes xCell Monocytes_Auch
Monocytie_lineage_MCPcounter_
Neutrophil_TIMER
Neutrophils CIBERSORT
Neutrophils_CIBERSORT ABS
Neutrophils MCPcounter Neutrophils_quantiseq. Neutrophils_xCell
NK_cols_activmed_CIBERSORT,
NK_cels_activaled CIBERBORT-AOS
NK_CBB_MCPcounter NK_cells_quantiseq.
NK_cells_resting_CIBERSORT.
NK_cells_resting_CIBERSORT_ABS
NK_cells_xCell NKcells EPIC NKT_xCell
Plasma_cols_CIBERSORT
Plasma_cels_CIBERSORT_ABS.
Plasma_cells_xCell,
CLP_xCell
CMP_xCell.
GMP_xCell
T_cells_MCPcounter
T_cells_gamma_delta_CIBERSORT T_cells_gamma_delta_CIBERSORT_ABS
Tgd_cells_xCell.
Cytotoxic_lymphocytes_MCPcounter
ESTIMATEScore_estimate
ImmuneScore_estimate-
StromalScore_estimate_
T_cells_regulatory_Tregs_CIBERSORT
T_cells_regulatory_Tregs_CIBERSORT_ABS
Tregs_quantiseq
Tregs_xCell

The relationship between NUP62 expression and immune infiltration. (A) Heatmap showing the correlation between NUP62 mRNA expression and

the expression of chemokines, chemokine receptors, immune inhibitors, immune stimulators, and MHC genes. (B) Correlation between NUP62

expression and tumor immune infiltration.

3.10 NUP62 may influence chemotherapy response

We further delved into the potential association between NUP62 expression levels and drug sensitivity. Results from both PRISM and CTRP drug sensitivity experiments revealed a significant correlation between NUP62 mRNA expression levels and drug sensitivity. The PRISM drug sensitivity experiment showed that the top three drugs positively correlated with NUP62 expression levels were trametinib, Ro-4987655, and AS-703026, while the top three negatively correlated drugs were idronoxil, sirolimus, and LY3023414 (Figure 9A). Similarly, the CTRP drug sensitivity experiment also demonstrated that the top three drugs positively correlated with NUP62 expression levels were selumetinib, saracatinib, and PD318088, while the top three negatively correlated drugs were axitinib, indisulam, and PRIMA- 1 (Figure 9B). The discovery of compounds that can significantly modulate NUP62 activity holds important potential value for developing novel and effective tumor treatment regimens. We found that arachidonyltrifluoromethane significantly impacts NUP62 expression across multiple tumor types (Figure 9C). Subsequently, we detailed the specific effects of various compounds on NUP62 in different tumor types and highlighted the significant diagnostic and prognostic predictive value of NUP62 in these tumor types (Figures 9D-F).

3.11 Validation of NUP62 expression levels in gastric cancer

To validate our initial hypothesis, we employed immunohistochemical techniques to conduct an in-depth analysis of NUP62 protein expression levels in gastric cancer tissues and their adjacent non-cancerous tissues. Figure 10A visually demonstrates the distribution and intensity of NUP62 expression in both non-cancerous and gastric cancer tissues. The research

findings indicate that the expression level of NUP62 in gastric

cancer tissues is significantly higher than that in their adjacent non- cancerous tissues (see Figure 10B). Furthermore, through quantitative PCR analysis, we found that the mRNA expression level of NUP62 in gastric cancer tissues is also significantly elevated compared to their adjacent non-cancerous tissues.

4 Discussion

Firstly, our study determined the expression and mutation status of NUP62 in humans. The results indicated that NUP62 expression is generally higher in cancer cells compared to normal cells. Mutations in NUP62 within tumors may lead to functional abnormalities of NUP62, affecting normal cellular functions and playing a role in tumorigenesis and progression. HPA fluorescence

A

C

The Connectivity Map

trametinib

<0.001

Ro-4987655

<0.001

AS-703026

<0.001

ACC

TAK-733

<0.001

BLCA

cobimetinib

<0.001

BRCA

AV-412

<0.001

ibrutinib

50.001

CESC

AZD8330

0.001

CHOL

neratinib

<0.001

COAD

XL-647

<0.001

<0.001

DLBC

afatinib

AZDB931

<0.001

ESCA

PD-153035

<0.001

GBM

finafloxacin

<0.001

MEK162

<0.001

HNSC

ZLN005

<0.001

KICH

AZ-628

<0.001

KIRC

erlotinib

<0.001

PD-198306

0.001

KIRP

bosutinib

<0.001

LAML

PD-318088

<0.001

LGG

XSum CMAP

pelitinib

<0.001

LIHC

0

saracatinib

<0.001

BIX-01294

<0.001

LUAD

BMS-690514

<0.001

LUSC

-0.5

ciclesonide

<0.001

MESO

gefitinib

<0.001

-1

LY3023414

<0.001

OV

sirolimus

<0.001

PAAD

idronoxil

<0.001

PCPG

-0.3

-0.2

-0.1

0.0

0.1

0.2

0.3

PRAD

B

NUP62 in PRISM (Relevance to AUC)

READ

SARC

selumetinib

<0.001

saracatinib

<0.001

SKCM

PD318088

<0.001

STAD

paclitaxel

<0.001

TGCT

phloretin

<0.001

ML050

<0.001

THCA

gemcitabine

<0.001

THYM

PF-184

<0.001

UCEC

topotecan

<0.001

niclosamide

<0.001

UCS

KX2-391

<0.001

UVM

ABT-737

<0.001

CD-437

<0.001

15

triazolothiadiazine

<0.001

ML239

<0.001

10

Number

KW-2449

<0.001

CCT036477

=0.001

5

LY-2183240

<0.001

etoposide

<0.001

fasudil

H.7

imatinib

STOCK1N.35696

trichostatin.A

X0225151.0000

clofibrate

STOCK1N.35874

TTNPB

NU.1025

AH.6809

exisulind

butein

mercaptopurine

PHA.00816795

W.13

MK.886

MS.275

X4.5.dianilinophthalimide

arachidonyltrifluoromethane

0

SB-743921

<0.001

mitomycin

<0.001

BRD-K34222889

<0.001

ceranib-2

<0.001

cucurbitacin

<0.001

AT7867

<0.001

BRD-K13999467

0.001

cerulenin

<0.001

PRIMA-1

<0.001

indisulam

<0.001

axitinib

<0.001

-0.3

-0.2

-0.1

0.0

0.1

0.2

0.3

NUP62 in CTRP (Relevance to AUC)

D

E

F

1.0

1.0

1.0-

0.5

0.5

0.5

score (ACC)

score (BLCA)

score (COAD)

0.0

0.0

0.0

MS.275

-0.5-

-0.5-

-0.5-

MK.886

mercaptopurine

PHA.00816795

MS.275

-1.0

arachidonyltrifluoromethane

-1.0

MK.88G

arachidonyItrifluoromethane

-1.0-

arachidonyltrifluoromethane

0

500

1000

0

500

1000

0

500

1000

number

number

number

FIGURE 9 Chemotherapy resistance analysis. Two different databases were utilized to investigate the association between NUP62 expression and drug sensitivity. (A) PRISM. (B) CTRP. (C) Identification of NUP62-targeting compounds through cMap analysis. (D-F) NUP62-targeting compounds with clinical significance.

images and staining data demonstrated that NUP62 is primarily localized in the nucleus, consistent with its role as a component of the NPC (20, 31). Further stratification according to patients’ clinical characteristics revealed that NUP62 is associated with the stages of multiple cancers, suggesting its role in tumor progression, invasion, and metastasis (32-34). The expression of NUP62 is also closely related to gender and age, which may be associated with

hormones, cellular aging, decreased DNA damage repair capacity, and increased cancer risk (35-37). In the TCGA dataset and the combined TCGA-GTEx dataset, NUP62 exhibited high AUC values across various cancers. Notably, in ESCA, the AUC value was >0.9, indicating that NUP62 has significant value for the screening and diagnosis of multiple cancers, particularly ESCA (38-40). Next, our research found that NUP62 expression is associated with the

A

E

- 200 um

200 μm

Paraneoplastic tissue

:

200 um

200 um

Gastric cancer tissue

B

Immunohistochemical Scoring Component

C

3-

P=0.0053

1.5-

P=0.0018

Relative mRNA expression of NUP62

1

2-

1.0-

-

1-

0.5.

0

0.0

paraneoplastic tissue

Gastric cancer tissue

paraneoplastic tissue

Gastric cancer tissue

FIGURE 10 (A, B) Represent the protein expression levels of NUP62 in gastric cancer tissues and their adjacent non-cancerous tissues. (C) Illustrates the mRNA expression levels of NUP62 in gastric cancer tissues versus their adjacent non-cancerous tissues.

prognosis of cancer patients, with differential expression observed in different types of cancer patients, as confirmed by KM curves. This suggests that NUP62 can serve as a prognostic marker for some cancers to guide the selection of treatment regimens for cancer patients, taking into account individual differences and formulating personalized treatment plans (41-43). Further GSEA analysis revealed that NUP62 is significantly associated with the cell

cycle, DNA damage, and DNA repair, which may be related to NUP62’s involvement in numerous tumor-related pathways and metabolic processes, including the G2/M checkpoint, E2F target genes, allograft rejection, and mitotic spindle. These pathways and processes play crucial roles in cancer cell proliferation, invasion, metastasis, and drug resistance, emphasizing the complexity and diversity of NUP62 in tumors (44-48). Subsequently, we

investigated the relationship between NUP62 and TCPA functional proteins, further demonstrating the important role of NUP62 in regulating the cell cycle of cancer patients.

We noted a decreasing trend in the expression level of NUP62 in testicular cancer tissues, which contrasts with the upregulation of NUP62 expression in most cancers. This difference in expression patterns may reflect the unique biology of testicular cancer. Therefore, we believe that testicular cancer is a validated experimental subject worthy of further study. By deeply exploring the expression regulation mechanism of NUP62 in testicular cancer, it is expected to provide new ideas and methods for the diagnosis and treatment of testicular cancer. We have not validated this for the time being, but it will be an important direction for future research.

Although our study found that high expression of NUP62 was associated with longer survival in these tumors, this does not mean that NUP62 is necessarily protective in these tumors. It is possible that the expression level of NUP62 is affected by a variety of factors, which may have different roles in different types of tumors. Therefore, the specific mechanism of action of NUP62 in these tumors and its interrelationship with other biomarkers can be further explored in subsequent studies, with the aim of providing more accurate targets for the treatment of these tumors.

Nuclear pore complex core proteins were shown to be extensively involved in metabolic reprogramming, such as NUP37, Nup210 (49-51). Therefore we speculated that NUP62 might be involved in metabolic pathways in tumors. In the subsequent analysis, we comprehensively analyzed the metabolic pathways in which NUP62 might be involved in pan-cancer by performing GSEA analysis on the KEGG metabolic gene set. The results indicated that NUP62 may inhibit multiple metabolic pathways in most tumors, such as steroid hormone synthesis, fatty acid synthesis, and glucose metabolism. This finding reveals a novel function of NUP62 as a nuclear pore complex core protein in tumor metabolism. Traditionally, NUP62 is mainly thought to be involved in the regulation of nuclear-plasmic transport and cellular signaling, whereas this study links it to metabolic pathways in tumors, providing a new perspective for understanding the role of NUP62 in tumorigenesis and progression.

Simultaneously, our research also revealed that the expression of NUP62 is significantly correlated with the infiltration levels of various immune cells, indicating that NUP62 can regulate the tumor immune microenvironment by influencing immune cells (52-55).

We delved into the correlation between NUP62 expression levels and immunotherapy response. The results showed that in clinical trials of immunotherapy, the response rate to immunotherapy was significantly higher in patients with high NUP62 expression than in patients with low expression. This suggests that NUP62 may be an important biomarker for predicting response to immunotherapy (56). By detecting patients’ NUP62 expression levels, physicians may be able to more accurately predict which patients are more likely to benefit from immunotherapy, thereby optimizing the treatment regimen and improving the therapeutic efficacy. The results of the ROC curve analysis showed the potential efficacy of NUP62 in predicting immunotherapy response, which provides a new way of thinking for the development of immunotherapy evaluation methods (57).

Traditional immunotherapy assessment methods mainly rely on indicators such as changes in tumor size, but these indicators often suffer from lag and inaccuracy. In contrast, biomarker-based prediction methods may be able to assess immunotherapy effects earlier and more accurately, thus guiding the adjustment and optimization of treatment regimens (58, 59). The study also found that the overall survival of patients with high NUP62 expression in immunotherapy was significantly longer than that of patients with low expression. This suggests that NUP62 may be associated with immunotherapy resistance (60). Through in-depth study of the mechanism of NUP62’s role in immunotherapy resistance, it may provide new strategies and methods for overcoming immunotherapy resistance, thereby improving the effectiveness of immunotherapy and patient survival (61, 62).

The significant negative correlation of NUP62 with CTLA-4 may imply that the low expression of NUP62 in these two cancers is associated with an attenuated immunosuppressive state (63). This suggests that NUP62 expression in these cancers may not be the primary factor promoting immune evasion or that it interacts in a complex manner with other immune regulatory mechanisms. The significant correlation of NUP62 with the expression of immune- related genes (e.g., THCA, KIRC, KIRP, and PCPG) in multiple cancer types suggests that it may play a role in a wide range of immune regulatory processes (64). This further supports the idea that NUP62 serves as an important node in the cancer immunoregulatory network. Given the significant correlation of NUP62 with immune-related genes, it could be a potential target for immunotherapy (65). By regulating the expression or function of NUP62, it may be able to influence the activity of immune cells and the immune status in the tumor microenvironment, thus providing new strategies for immunotherapy. In the future, further in-depth studies can be conducted to investigate how NUP62 interacts with immune-related genes and how these interactions affect the tumor immune microenvironment and immunotherapeutic effects (66).

The findings suggest that NUP62 expression level may become a marker for predicting the efficacy of specific drugs (67). This means that in future clinical practice, physicians can select the most likely effective drugs based on a patient’s NUP62 expression level, enabling more precise and personalized treatment (68). By understanding the relationship between NUP62 and drug sensitivity, drug developers can design new drugs in a more targeted manner, especially for those drugs that are positively or negatively correlated with NUP62 expression levels (69). This biomarker-based drug development strategy is expected to improve the success rate and clinical application of new drugs. In addition, by monitoring changes in NUP62 expression levels, physicians can make timely adjustments to treatment regimens to avoid or delay the onset of drug resistance. In addition, the efficacy of certain drugs, such as those negatively correlated with NUP62 expression levels, may be affected by NUP62 expression levels, and therefore, other alternative drugs or combination strategies may need to be explored in specific patient populations (70, 71). The present study reveals a novel role of NUP62 in drug sensitivity, which provides important clues for an in-depth understanding of its biological functions and regulatory mechanisms. This can help to further expand the application areas of biomarker research and improve the diagnosis and treatment of diseases such as cancer (72-74).

Additionally, we validated the differential expression of NUP62 between gastric cancer tissues and their adjacent non- cancerous tissues.

In summary, this study revealed the expression characteristics of NUP62 in multiple malignancies and its associations with tumor- related pathways, clinical features, prognosis, and tumor immunity. These findings provide crucial insights into our understanding of the role of NUP62 in tumors and offer new targets for future tumor diagnosis and immunotherapy. However, to fully uncover the mechanisms of NUP62 in tumors, further in-depth research and exploration are required.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Files, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by the Ethics Committee of the Guizhou Provincial People’s Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

LC: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing - original draft. YH: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing - original draft. MD: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Supervision, Validation, Visualization, Writing - original draft. TY: Data curation, Investigation, Methodology, Software, Writing - original draft. YC: Data curation, Investigation, Methodology, Writing - original draft. BW: Data curation, Investigation, Methodology, Software, Writing - original draft. DC: Investigation, Software, Writing - review & editing. CL: Investigation, Methodology, Software, Supervision, Writing - review & editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the [Guizhou Provincial Health and Wellness Commission] under Grant [gzwkj 2021-102]; The National Natural Science Foundation of China (No. 82260084); Guizhou Provincial Science and Technology Agency Project (Qian Ke He Foundation ZK [2023] General 217); and Support by Key Advantageous Discipline Construction Project of Guizhou Provincial Health Commission in 2023.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2025.1559396/ full#supplementary-material

SUPPLEMENTARY FIGURE 1

Expression Levels of NUP62 in Cell Line Models, Cancerous, and Normal Tissues. (A) Immunofluorescence images and merged images of the nuclei, NUP62 protein, microtubules, and endoplasmic reticulum (ER) in A-431, U2OS, and U-251MG cell lines obtained from the HPA database. (B) IHC images of NUP62 protein expression in normal and tumor tissues obtained from the HPA database.

References

1. Sun S, Wang YH, Gao X, Wang HY, Zhang L, Wang N, et al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA. Cancer J Clin. (2024) 74:229-63. doi: 10.3322/ caac.21834

2. Siegel RL, Miller KD, Wagle NS, Jemal A. Cancer statistics, 2023. CA. Cancer J Clin. (2023) 73:17-48. doi: 10.3322/caac.21763

3. Wu F, Yang J, Liu J, Wang Y, Mu J, Zeng Q, et al. Signaling pathways in cancer- associated fibroblasts and targeted therapy for cancer. Signal Transduction Targeting Ther. (2021) 6:1-35. doi: 10.1038/s41392-021-00641-0

4. Punekar SR, Velcheti V, Neel BG, Wong K-K. The current state of the art and future trends in RAS-targeted cancer therapies. Nat Rev Clin Oncol. (2022) 19:637-55. doi: 10.1038/s41571-022-00671-9

5. Pérez-Herrero E, Fernández-Medarde A. Advanced targeted therapies in cancer: Drug nanocarriers, the future of chemotherapy. Eur J Pharm Biopharm. (2015) 93:52- 79. doi: 10.1016/j.ejpb.2015.03.018

6. Zhang Y, Yan L, Liu A, Li F, Li Y, Zhang Y. Cell polarity-related gene PTK7, a potential diagnostic biomarker in pan-cancer. Curr Med Chem. (2024) 31:1-18. doi: 10.2174/0109298673313999240816103054

7. Zhang Y, Xue Y, Gao Y, Zhang Y. Prognostic and predictive value of pathohistological features in gastric cancer and identification of SLITRK4 as a potential biomarker for gastric cancer. Sci Rep. (2024) 14:29241. doi: 10.1038/ s41598-024-80292-7

8. Cordani M, Dando I, Ambrosini G, González-Menéndez P. Signaling, cancer cell plasticity, and intratumor heterogeneity. Cell Commun Signal. (2024) 22:255. doi: 10.1186/s12964-024-01643-5

9. DeBerardinis RJ, Chandel NS. Fundamentals of cancer metabolism. Sci Advances. 2(5):e1600200. doi: 10.1126/sciadv.1600200

10. Vaghari-Tabari M, Ferns GA, Qujeq D, Andevari AN, Sabahi Z, Moein S. Signaling, metabolism, and cancer: An important relationship for therapeutic intervention. J Cell Physiol. (2021) 236:5512-32. doi: 10.1002/jcp.v236.8

11. Vander Heiden MG, Cantley LC, Thompson CB. Understanding the warburg effect: the metabolic requirements of cell proliferation. Science. (2009) 324:1029-33. doi: 10.1126/science.1160809

12. Gyamfi J, Kim J, Choi J. Cancer as a metabolic disorder. Int J Mol Sci. (2022) 23:1155. doi: 10.3390/ijms23031155

13. Dey P, Kimmelman AC, DePinho RA. Metabolic codependencies in the tumor microenvironment. Cancer Discovery. (2021) 11:1067-81. doi: 10.1158/2159-8290.CD- 20-1211

14. Koppenol WH, Bounds PL, Dang CV. Otto Warburg’s contributions to current concepts of cancer metabolism. Nat Rev Cancer. (2011) 11:325-37. doi: 10.1038/ nrc3038

15. Ashton TM, McKenna WG, Kunz-Schughart LA, Higgins GS. Oxidative phosphorylation as an emerging target in cancer therapy. Clin Cancer Res. (2018) 24:2482-90. doi: 10.1158/1078-0432.CCR-17-3070

16. Chelakkot C, Chelakkot VS, Shin Y, Song K. Modulating glycolysis to improve cancer therapy. Int J Mol Sci. (2023) 24:2606. doi: 10.3390/ijms24032606

17. Wang Z-H, Peng W-B, Zhang P, Yang X-P, Zhou Q. Lactate in the tumour microenvironment: From immune modulation to therapy. eBioMedicine. (2021) 73. doi: 10.1016/j.ebiom.2021.103627

18. Ngwa VM, Edwards DN, Philip M, Chen J. Microenvironmental metabolism regulates antitumor immunity. Cancer Res. (2019) 79:4003-8. doi: 10.1158/0008- 5472.CAN-19-0617

19. Chuang Y-M, Tzeng S-F, Ho P-C, Tsai C-H. Immunosurveillance encounters cancer metabolism. EMBO Rep. (2024) 25:471-88. doi: 10.1038/s44319-023-00038-w

20. Kalverda B, Pickersgill H, Shloma VV, Fornerod M. Nucleoporins directly stimulate expression of developmental and cell-cycle genes inside the nucleoplasm. Cell. (2010) 140:360-71. doi: 10.1016/j.cell.2010.01.011

21. Echeverría PC, Mazaira G, Erlejman A, Gomez-Sanchez C, Piwien Pilipuk G, Galigniana MD. Nuclear import of the glucocorticoid receptor-hsp90 complex through the nuclear pore complex is mediated by its interaction with nup62 and importin B. Mol Cell Biol. (2009) 29:4788-97. doi: 10.1128/MCB.00649-09

22. Zhou T, Li S, Zhong W, Vihervaara T, Beaslas O, Perttila J, et al. OSBP-related protein 8 (ORP8) regulates plasma and liver tissue lipid levels and interacts with the nucleoporin nup62. PloS One. (2011) 6:e21078. doi: 10.1371/journal.pone.0021078

23. Simon DN, Rout MP. Cancer and the nuclear pore complex. In: Schirmer EC, de las Heras JI, editors. Cancer biology and the nuclear envelope: recent advances may elucidate past paradoxes. Springer, New York, NY (2014). p. 285-307. doi: 10.1007/978- 1-4899-8032-8_13

24. Singh U, Bindra D, Samaiya A, Mishra RK. Overexpressed Nup88 stabilized through interaction with Nup62 promotes NF-KB dependent pathways in cancer. Front Oncol. (2023) 13. doi: 10.3389/fonc.2023.1095046

25. Betancor YZ, Ferreiro-Pantin M, Anido-Herranz U, Fuentes-Losada M, Leon- Mateos L, Garcia-Acuna SM, et al. A three-gene expression score for predicting clinical benefit to anti-PD-1 blockade in advanced renal cell carcinoma. Front Immunol. (2024) 15. doi: 10.3389/fimmu.2024.1374728

26. Liu H, Yan Y, Chen R, Zhu M, Lin J, He C, et al. Integrated nomogram based on five stage-related genes and TNM stage to predict 1-year recurrence in hepatocellular carcinoma. Cancer Cell Int. (2020) 20:140. doi: 10.1186/s12935-020-01216-9

27. Borlido J, D’Angelo MA. Nup62-mediated nuclear import of p63 in squamous cell carcinoma. EMBO Rep. (2018) 19:3-4. doi: 10.15252/embr.201745497

28. Hazawa M, Lin DC, Kobayashi A, Jiang YY, Xu L, Dewi FRP, et al. ROCK-dependent phosphorylation of NUP62 regulates p63 nuclear transport and squamous cell carcinoma proliferation. EMBO Rep. (2018) 19:73-88. doi: 10.15252/embr.201744523

29. Shen W, Song Z, Zhong X, Huang M, Shen D, Gao P, et al. Sangerbox: A comprehensive, interaction-friendly clinical bioinformatics analysis platform. Imeta (2022) 1(3):e36. doi: 10.1002/imt2.36

30. Liu Z, Liu L, Weng S, Xu H, Xing Z, Ren Y, et al. BEST: a web application for comprehensive biomarker exploration on large-scale data in solid tumors. J Big Data. (2023) 10:165. doi: 10.1186/s40537-023-00844-y

31. Madheshiya PK, Shukla E, Singh J, Bawaria S, Ansari MY, Chauhan RJMBotC. Insights into the role of Nup62 and Nup93 in assembling cytoplasmic ring and central transport channel of the nuclear pore complex. Mol Biol Cell. (2022) 33(14):ar139. doi: 10.1091/mbc.E22-01-0027

32. Klein CA. Cancer progression and the invisible phase of metastatic colonization. Nat Rev Cancer. (2020) 20:681-94. doi: 10.1038/s41568-020-00300-6

33. Kumar R, Paul AM, Amjesh R, George B, Pillai MR. Coordinated dysregulation of cancer progression by the HER family and p21-activated kinases. Cancer Metastasis Rev. (2020) 39:583-601. doi: 10.1007/s10555-020-09922-6

34. Wang L, Guo X-L. Molecular regulation of galectin-3 expression and therapeutic implication in cancer progression. Biomed Pharmacother. (2016) 78:165-71. doi: 10.1016/j.biopha.2016.01.014

35. Iwatsubo T, Ishihara R, Morishima T, Maekawa A, Nakagawa K, Arao M, et al. Impact of age at diagnosis of head and neck cancer on incidence of metachronous cancer. BMC Cancer. (2019) 19:3. doi: 10.1186/s12885-018-5231-7

36. Pawelec G, Derhovanessian E, Larbi A. Immunosenescence and cancer. Crit Rev Oncol Hematol. (2010) 75:165-72. doi: 10.1016/j.critrevonc.2010.06.012

37. Muermann MM, Wassersug RJ. Prostate cancer from a sex and gender perspective: A review. Sex Med Rev. (2022) 10:142-54. doi: 10.1016/j.sxmr.2021.03.001

38. Menz A, Weitbrecht T, Gorbokon N, Buscheck F, Luebke AM, Kluth M, et al. Diagnostic and prognostic impact of cytokeratin 18 expression in human tumors: a tissue microarray study on 11,952 tumors. Mol Med. (2021) 27:16. doi: 10.1186/s10020- 021-00274-7

39. Krings G, Nystrom M, Mehdi I, Vohra P, Chen Y-Y. Diagnostic utility and sensitivities of GATA3 antibodies in triple-negative breast cancer. Hum Pathol. (2014) 45:2225-32. doi: 10.1016/j.humpath.2014.06.022

40. Ding N, Li M, Zhao X. PHF5A is a potential diagnostic, prognostic, and immunological biomarker in pan-cancer. Sci Rep. (2023) 13:17521. doi: 10.1038/ s41598-023-44899-6

41. Zhang L, Barritt GJ. TRPM8 in prostate cancer cells: a potential diagnostic and prognostic marker with a secretory function? Endocr. Relat Cancer. (2006) 13:27-38. doi: 10.1677/erc.1.01093

42. Mesci A, Lucien F, Huang X, Wang EH, Shin D, Meringer M, et al. RSPO3 is a prognostic biomarker and mediator of invasiveness in prostate cancer. J Transl Med. (2019) 17:125. doi: 10.1186/s12967-019-1878-3

43. Schinke H, Shi E, Lin Z, Quadt T, Kranz G, Zhou J, et al. A transcriptomic map of EGFR-induced epithelial-to-mesenchymal transition identifies prognostic and therapeutic targets for head and neck cancer. Mol Cancer. (2022) 21:178. doi: 10.1186/s12943-022-01646-1

44. Icard P, Fournel L, Wu Z, Alifano M, Lincet H. Interconnection between metabolism and cell cycle in cancer. Trends Biochem Sci. (2019) 44:490-501. doi: 10.1016/j.tibs.2018.12.007

45. Evan GI, Vousden KH. Proliferation, cell cycle and apoptosis in cancer. Nature. (2001) 411:342-8. doi: 10.1038/35077213

46. Lord CJ, Ashworth A. The DNA damage response and cancer therapy. Nature. (2012) 481:287-94. doi: 10.1038/nature10760

47. Basu AK. DNA damage, mutagenesis and cancer. Int J Mol Sci. (2018) 19:970. doi: 10.3390/ijms19040970

48. Kang T-H. DNA damage, repair, and cancer metabolism. Int J Mol Sci. (2023) 24:16430. doi: 10.3390/ijms242216430

49. Liu Z, Hu Q, Luo Q, Zhang G, Yang W, Cao K, et al. NUP37 accumulation mediated by TRIM28 enhances lipid synthesis to accelerate HCC progression. Oncogene. (2024) 43:3255-67. doi: 10.1038/s41388-024-03167-1

50. Zaitsava H, Gachowska M, Bartoszewska E, Kmiecik A, Kulbacka J. The potential of nuclear pore complexes in cancer therapy. Molecules. (2024) 29:4832. doi: 10.3390/molecules29204832

51. Han F, Fan X, Hu M, Wen J, Wang J, Zhang D, et al. Nup210 promotes colorectal cancer progression by regulating nuclear plasma transport. Lab Invest. (2024) 104. doi: 10.1016/j.labinv.2024.102149

52. Gajewski TF, Schreiber H, Fu Y-X. Innate and adaptive immune cells in the tumor microenvironment. Nat Immunol. (2013) 14:1014-22. doi: 10.1038/ni.2703

53. Fu T, Dai L-J, Wu S-Y, Xiao Y, Ma D, Jiang Y-Z, et al. Spatial architecture of the immune microenvironment orchestrates tumor immunity and therapeutic response. J Hematol Oncol.J Hematol Oncol. (2021) 14:98. doi: 10.1186/s13045-021-01103-4

54. Lv B, Wang Y, Ma D, Cheng W, Liu J, Yong T, et al. Immunotherapy: reshape the tumor immune microenvironment. Front Immunol. (2022) 13. doi: 10.3389/ fimmu.2022.844142

55. Mei K, Chen Z, Huang L, Wang J, Wei Y. Correlation between the immune microenvironment and bladder cancer based on a prognostic miRNA risk model. Cancer Insight. (2024) 3:14-25. doi: 10.58567/ci03020002

56. Liu Y, Zhang B, Wu X, Wang F, Yang Z, Li M, et al. A facile liquid biopsy assay for highly efficient CTCs capture and reagent-less monitoring of immune checkpoint PD-L1 expression on CTCs with non-small cell lung cancer patients. Biosens. Bioelectron. (2025) 275:117236. doi: 10.1016/j.bios.2025.117236

57. Zhang N, Li C, Zhao Z, Jiang B, Wang W, Sun F, et al. Immune microenvironment features underlying the superior efficacy of neoadjuvant

immunochemotherapy over chemotherapy in local advanced gastric cancer. Front Immunol. (2025) 16. doi: 10.3389/fimmu.2025.1497004

58. Wang X, Li J, Zhu Y, Shen H, Ding J, Zeng T, et al. Targeting ADAR1 with a small molecule for the treatment of prostate cancer. Nat Cancer. (2025) 1-19. doi: 10.1038/s43018-025-00907-4

59. Naing A, Mckean M, Tolcher A, Victor A, Hu P, Gao W, et al. TIGIT inhibitor M6223 as monotherapy or in combination with bintrafusp alfa in patients with advanced solid tumors: a first-in-human, phase 1, dose-escalation trial. J Immunother. Cancer. (2025) 13:e010584. doi: 10.1136/jitc-2024-010584

60. Ghorbaninezhad F, Nour MA, Farzam OR, Saeedi H, Vanan AG, Bakhshivand M, et al. The tumor microenvironment and dendritic cells: Developers of pioneering strategies in colorectal cancer immunotherapy? Biochim Biophys Acta BBA - Rev Cancer. (2025) 1880:189281. doi: 10.1016/j.bbcan.2025.189281

61. Sardar P, Beresford-Jones BS, Xia W, Shabana O, Suyama S, Ramos RJF, et al. Gut microbiota-derived hexa-acylated lipopolysaccharides enhance cancer immunotherapy responses. Nat Microbiol 1-13. (2025). doi: 10.1038/s41564-025-01930-y

62. Pan S, Wang Z. Antiviral therapy can effectively suppress irAEs in HBV positive hepatocellular carcinoma treated with ICIs: validation based on multi machine learning. Front Immunol. (2025) 15. doi: 10.3389/fimmu.2024.1516524

63. Wu Y, Yin M, Xia W, Dou B, Liu X, Sun R. Enhancing NK cell antitumor activity with natural compounds: research advances and molecular mechanisms. Phytother. Res. (2025). doi: 10.1002/ptr.8456

64. Yang J, Tang S, Saba NF, Shay C, Teng Y. Tumor secretome shapes the immune landscape during cancer progression. J Exp Clin Cancer Res. (2025) 44:47. doi: 10.1186/ s13046-025-03302-0

65. Huang M, Huang Z, Miao S, Chen X, Tan Y, Zhou Y, et al. Bioinformatics Analysis of coagulation-related genes in lung adenocarcinoma: unveiling prognostic indicators and treatment pathways. Sci Rep. (2025) 15:4972. doi: 10.1038/s41598-025- 87669-2

66. Liu R, Ji Z, Wang X, Zhu L, Xin J, Ma L, et al. Regorafenib plus sintilimab as a salvage treatment for microsatellite stable metastatic colorectal cancer: a single-arm, open-label, phase II clinical trial. Nat Commun. (2025) 16:1481. doi: 10.1038/s41467- 025-56748-3

67. Sharma D, Arumugam S. Pharmacophore-based identification and in Silico characterization of microbial metabolites as potential modulators of Wnt signaling pathway in colorectal cancer therapy. Mol Divers. (2025). doi: 10.1007/s11030-024- 11103-4

68. Wan X, Jiang M, Madan S. Research progress of nanomedicine for tumor immunotherapy. Cancer Insight. (2023) 3:83-98. doi: 10.2174/ 0109298673313999240816103054

69. Zhang C, Chen YY, Chen S, Wang Y, Yuan Y, Yang X, et al. Characterization of an enhancer RNA signature reveals treatment strategies for improving immunotherapy efficacy in cancer. Cancer Res. (2025). doi: 10.1158/0008-5472.CAN-24-2289

70. Long H, Zhou J, Zhou C, Xie S, Wang J, Tan M, et al. Proteomic characterization of liver cancer cells treated with clinical targeted drugs for hepatocellular carcinoma. Biomedicines. (2025) 13:152. doi: 10.3390/biomedicines13010152

71. Zhao L, Shang J, Meng X, He X, Zhang Y, Liu JX. Adaptive multi-kernel graph neural network for drug-drug interaction prediction. Interdiscip. Sci Comput Life Sci. (2025). doi: 10.1007/s12539-024-00684-1

72. Kuai X, Wei C, He X, Wang F, Wang C, Ji J. The potential value of RPS27A in prognosis and immunotherapy: from pan-cancer analysis to hepatocellular carcinoma validation. ImmunoTargets Ther. (2024) 13:673-90. doi: 10.2147/ITT.S493217

73. Liu S, He M, Sun H, Wu Y, Jin W. 5-hydroxytryptamine G-protein-coupled receptor family genes: key players in cancer prognosis, immune regulation, and therapeutic response. Genes. (2024) 15:1541. doi: 10.3390/genes15121541

74. Zhou T, Zhang DD, Jin J, Xie J, Yu J, Zhu C, et al. Multiomic characterization, immunological and prognostic potential of SMAD3 in pan-cancer and validation in LIHC. Sci Rep. (2025) 15:657. doi: 10.1038/s41598-024-84553-3