Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 04 October 2024
Sec. Cancer Immunity and Immunotherapy
This article is part of the Research Topic Unveiling Biomarkers and Mechanisms in the Tumor-Immune Nexus View all 28 articles

Integrated transcriptome analysis of CSE1L regarding poor prognosis and immune infiltration in bladder urothelial carcinoma and experimental verification

  • 1Department of Urology, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
  • 2Beijing National Day School, Beijing, China

Background: Bladder urothelial carcinoma (BLCA) is one of the most prevalent tumors globally, with its incidence rising notably in developed countries, significantly affecting human health. CSE1L encodes a protein that is involved in various cellular processes and plays a critical role in cancer initiation and progression. However, its role in BLCA remains underexplored.

Methods: CSE1L expression in BLCA was analyzed using TCGA data and validated by qRT-PCR and Western blot in clinical samples. Survival analysis and Cox regression models were used to evaluate its prognostic value. Functional enrichment and protein interaction analyses were performed, and immune cell infiltration was assessed using CIBERSORT. Drug sensitivity was analyzed using GDSC data. In vitro assays evaluated the effects of CSE1L knockdown on cell proliferation, migration, and invasion.

Results: CSE1L was found to be significantly overexpressed in BLCA tissues compared to normal tissues. High CSE1L expression was associated with poor overall survival and unfavorable clinicopathological features. Functional enrichment analysis revealed that DEGs related to CSE1L were involved in cell cycle regulation and immune-related pathways. Immune infiltration analysis indicated a significant correlation between CSE1L expression and various immune cell types, particularly T cells and macrophages. Drug sensitivity analysis identified several chemotherapeutic agents, including MG-132, Palbociclib, and Nutlin-3a, which were more effective in the low-CSE1L expression group, while the high-CSE1L expression group showed sensitivity to drugs like S-Trityl-L-cysteine, Bleomycin, and Cisplatin. In vitro knockdown of CSE1L in BLCA cell lines inhibited cell proliferation, migration, and invasion.

Conclusions: The overexpression of CSE1L is associated with the progression and poor prognosis of bladder cancer, suggesting it could be a promising target for bladder cancer in the future.

Introduction

Bladder urothelial carcinoma (BLCA) remains a significant global health concern, ranking 9th in incidence and 13th in cancer-related deaths worldwide (1). Treatment options for BLCA depend on various factors, including tumor stage, grade, and patient characteristics. Transurethral resection of bladder tumor (TURBT) is the standard treatment for non-muscle-invasive bladder cancer (NMIBC). For muscle-invasive bladder cancer (MIBC), radical cystectomy, with or without neoadjuvant chemotherapy, is often recommended (2, 3). Despite advances in understanding BLCA, challenges persist, including tumor heterogeneity, risk stratification, and treatment response prediction. Additionally, the identification of novel biomarkers, genetic alterations, and therapeutic targets holds promise for personalized medicine approaches.

CSE1L (chromosome segregation 1-like protein), also known as CAS (cellular apoptosis susceptibility protein), encodes a protein involved in various cellular processes and is found in humans and many other organisms. CSE1L is primarily recognized for its role in nucleocytoplasmic transport, the process by which molecules such as proteins and RNA move between the nucleus and the cytoplasm of a cell (4). CSE1L functions as part of the nuclear pore complex (NPC), serving as a gateway for the transport of molecules in and out of the nucleus (5). Additionally, CSE1L has been implicated in other cellular functions, including cell cycle regulation, cell division, and apoptosis (6, 7). Research on the CSE1L gene and protein has highlighted its importance in maintaining cellular homeostasis and its potential involvement in disease. Dysregulation of regulatory mechanisms can lead to aberrant CSE1L expression, contributing to developmental defects or disease progression (8). Further research is needed to deepen our understanding of the precise mechanisms by which CSE1L exerts its functions in development and disease. Unraveling the upstream regulators and downstream targets of CSE1L will provide insights into its molecular networks and potential therapeutic targets. Additionally, investigating the utility of CSE1L as a diagnostic or prognostic marker in BLCA may open avenues for personalized medicine approaches.

The present study aimed to explore the correlation between CSE1L and BLCA and identify potential biomarkers for this disease. The prognostic value of CSE1L expression was analyzed using survival analysis, and its correlation between CSE1L expression and clinicopathological characteristics in BLCA patients was examined.

Materials and methods

Data acquisition

Gene expression profile data and clinical information were collected from the TCGA-BLCA project, specifically utilizing the HTSeq-FPKM workflow. The data were sourced from the TCGA database (https://portal.gdc.cancer.gov/). To ensure a comprehensive analysis, RNA-seq data lacking clinical information were excluded, resulting in a total of 431 cases. These cases included both BLCA specimens and adjacent non-tumor specimens. Any unavailable or unknown clinical features were treated as missing values in the analysis.

Survival analysis and construction of nomogram

Survival analyses of CSE1L were conducted using the R package survival. Patients were divided into low and high CSE1L expression groups. Univariate and multivariate Cox regression analyses of the clinical characteristics were used to establish a risk score as an independent prognostic predictor. To accurately estimate each patient’s prognosis, a nomogram was constructed based on the risk score and conventional clinical characteristics, including age, gender, stage, and TNM classification. The nomogram was then evaluated using calibration curves.

GO and KEGG enrichment analyses

The clusterProfiler R package was used to identify significantly differentially expressed genes (DEGs) between BLCA samples and normal samples using an unpaired t-test. The Benjamini-Hochberg method was applied to adjust the threshold value, with an adjusted P <0.05 and |log FC| >1. Gene Ontology (GO) analysis revealed that these genes were associated with a diverse range of functional categories, including biological process and molecular function. Additionally, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs were conducted using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) online tools (https://david.ncifcrf.gov/), with a P-value cut-off criterion of <0.05.

Composition of invasive immune cells in BLCA

The marker genes for 24 immune cell types were identified based on a literature review. The infiltration of these immune cells in BLCA was analyzed using the CIBERSORT algorithm and single-sample Gene Set Enrichment Analysis(ssGSEA) (9). The correlation between CSE1L expression and immune cell infiltration was assessed using the Wilcoxon method.

Protein-protein interaction networks

The Pearson method was applied to calculate the correlation coefficient, with an adjusted |Cor| ≥ 0.3 and P ≤ 0.05. A protein-protein interaction (PPI) network was constructed using the STRING database (https://string-db.org) based on the intersection of genes between the screened DEGs and CSE1L-correlation genes. Subsequently, Gene Ontology Biological Process (GO-BP) enrichment analysis was performed on intersection genes, and the PPI network was annotated.

Drug sensitivity analysis

Based on the Genomics of Drug Sensitivity in Cancer (GDSC, https://www.cancerrxgene.org/), we used the R package prophytic to predict the drug sensitivity of each tumor sample. The estimated half-maximal inhibitory concentration (IC50) of each compound was calculated. Representative drugs were then selected for Spearman correlation analysis.

Patient samples

We collected BLCA tissues and matched normal bladder tissues from BLCA patients who underwent TURBT in the Urology Department of Beijing Tiantan Hospital, Capital Medical University, between February 2024 and May 2024. During the operation, the surgical specimens were promptly immersed in liquid nitrogen and stored at -80°C for subsequent use. All participants provided informed consent by signing the required form, and the project was approved by the Ethics Committee of Beijing Tiantan Hospital, Capital Medical University [KY2024-005-01].

Cell lines and cell culture

The human normal urothelial cell line (SV-HUC-1) and BLCA cell lines (5637 and T24) were purchased from the Cell Bank of the Chinese Academy of Sciences. SV-HUC-1 cells were cultured in F-12K medium (Gibco), while T24 and 5637 cells were cultured in RPMI-1640 medium (Gibco). The medium was supplemented with 10% fetal bovine serum (FBS) and 1% penicillin and streptomycin. All cells were incubated at 37°C in a humidified environment with 5% carbon dioxide.

Transfection

Small interfering CSE1L RNA (siR-CSE1L) and its corresponding negative control siRNAs (NC) were synthesized by General Biol Company (China). NC and siR-CSE1L were transfected into BLCA cell lines using Lipofectamine 3000 (Invitrogen, USA). All primer sequences are listed in Table 1.

Table 1
www.frontiersin.org

Table 1. Primer sequences used in this research.

RNA extraction, reverse transcription, and quantitative real-time PCR

Total RNA was extracted from freshly frozen tissues or bladder cancer cell lines using QIAGEN/74136/RNeazy Plus Mini Kit. Reverse transcription was then performed using the Vazyme/R312-02/HiScript III 1st Strand cDNA Synthesis Kit (+gDNA wiper). Subsequently, qRT-PCR was carried out using Vazyme/Q712-03/Taq Pro Universal SYBR qPCR Master Mix. The housekeeping gene GAPDH was used as an internal reference, and the relative expression level of the target gene was calculated using the 2−ΔΔCT method. Each experiment was performed in three replicates, and the primer sequences used are provided in Table 1.

Western blot

Protein extraction and Western blotting were performed as follows. Tumor tissue was cut into small pieces, each weighing approximately 30mg, and placed in 500 μL RIPA lysis buffer (pre-added with protease and phosphatase inhibitors). A sample of 20 μg of protein was then loaded onto an SDS- PAGE gel and run under standard conditions, with the same volume of culture medium used as a control. The proteins were transferred to a PVDF membrane using iBlot 2 system following the standard protocol. The membrane was blocked with 5% nonfat milk in TBST and incubated for 1 hour at room temperature, then incubated overnight at 4°C with agitation in a solution containing the primary antibody diluted to the recommended concentration in TBST. After washing the blot in TBST three to four times for 5 minutes each at room temperature, the membrane was incubated with a horseradish peroxidase (HRP)-conjugated secondary antibody, diluted to the recommended concentration, for 1 hour at room temperature. Protein bands were visualized using a Clinx Gel Documentation and Analysis instrument, while ImageJ was used for the analysis of the gray values of the bands.

Cell counting kit-8 assay

5637 and T24 cells transfected with NC and siR-CSE1L were plated in 96-well plates at a density of 5×10^3 cells per well. The cells were pre-cultured in a CO2 incubator at 37°C for 24 hours. CCK-8 reagent was added at 24, 48, 72 and 96 hours. After a 2-hour incubation at 37°C, absorbance at 450nm was measured using a spectrophotometer.

Cell migration assay

Cells were seeded at a density of 2×10^5 cells per well in a 6-well plate and cultured for 24 hours until they adhered. A P200 pipette tip was used to create a vertical scratch. The cells were then washed three times with PBS to remove the dislodged cells, and images were captured immediately using a microscope (0h point). Afterward, serum-free medium was added, and the cells were incubated for an additional 12 hours before capturing images at the same locations. The results were analyzed using ImageJ software, with the wound healing area as an indicator of cell migration.

Transwell invasion assay

For the cell invasion assay, the transwell plate with an 8-mm diameter polycarbonate membrane was used. The membrane in the upper chamber was coated with a layer of Matrigel diluted in medium, and the cells were resuspended in serum-free medium. The lower chamber was filled with complete medium. After incubation at 37°C in a humidified atmosphere with 5% carbon dioxide for 24-72 hours, the cells on the lower surface of the membrane were fixed with 100% methanol, stained with 0.1% crystal violet for 15-30 minutes, and counted under the microscope.

Statistical analysis

The associations between clinical factors and CSE1L were analyzed using various statistical tests, including the Wilcoxon rank-sum test and logistic regression, and the t-test for numerical variables. The chi-square test was also employed for categorical data. Prognostic data were obtained from TCGA, and Cox regression analyses, along with the Kaplan-Meier method, were used to evaluate prognostic factors. In all tests, P value <0.05 was considered statistically significant.

Results

CSE1L expression in pan−cancer and BLCA

First, pan-cancer analyses were conducted to compare the CSE1L expression between tumor samples and their corresponding normal samples, using the Wilcoxon rank-sum test. Most human cancers, including BLCA, BRCA, CESC, CHOL, COAD, ESCA, HNSC, KICH, KIRC, KIRP, LIHC, LUAD, LUSC, READ, STAD and UCEC exhibit overexpression of CSE1L. In BLCA, a significant increase in CSE1L expression was observed in tumors compared to normal tissues (p < 0.05) (Figure 1A). Next, we compared the expression of CSE1L in 19 normal samples and 412 BLCA samples from the TCGA-BLCA dataset. CSE1L expression was significantly high in cancer samples, with a median level of 6.924 and 6.211 in tumor and normal tissues, respectively (p < 0.001) (Figure 1B). Subsequently, we grouped the samples based on high- and low-CSE1L expression groups, analyzed the DEGs, and generated heatmaps (Figures 1B, C).

Figure 1
www.frontiersin.org

Figure 1. CSE1L expression in TCGA database and DEGs. (A) The mRNA expression levels of CSE1L in pan-cancer; (B) Density heatmap of DEGs; (C) Heatmap of top 100 DEGs. **p < 0.01; ***p < 0.001.

CSE1L could be an independent prognostic factor for BLCA

Clinical data and gene expression data for 412 BLCA samples were downloaded from TCGA database. We investigated the clinicopathological characteristics of BLCA patients with differential CSE1L expression. As shown in Table 2, compared with the low-CSE1L group, patients in the high-CSE1L group exhibited a higher proportion of histologic grade (P < 0.001), and worse primary therapy outcome (P = 0.003). There was also a significant correlation between CSE1L levels and smoking status (P = 0.013) as well as overall survival (OS) (P < 0.001). To explore the correlation between CSE1L expression and BLCA prognosis, we compared the survival rates of high-and low-CSE1L expression groups using Kaplan-Meier analysis. The results showed that high CSE1L expression was significantly associated with poor OS in BLCA patients (P=0.002) (Figure 2C).

Table 2
www.frontiersin.org

Table 2. Clinical pathological characteristics of the patients.

Figure 2
www.frontiersin.org

Figure 2. High expression of CSE1L indicates worse prognosis. (A) The nomogram to estimate the 1-, 3-, and 5-year OS of bladder cancer patients. (B) The calibration curves of actual values and predicted values provided by the nomogram. (C) Kaplan-Meier survival curves indicated that bladder cancer patients with high CASE1L mRNA expression had a shorter OS.

A nomogram model was constructed using multivariate Cox regression based on the relative expression of the key gene CSE1L and its main clinical observation indicators. The results indicate that CSE1L is an independent prognostic factor in patients with BLCA, with the prognosis worsening as the risk score increased (Figure 2A). The calibration curve demonstrated that the actual OS values were consistent with those predicted by the nomogram, with an overall concordance index of 0.704 (Figure 2B).

Functional enrichment and analyses

We performed DEGs analysis using data from the TCGA cohort. The GO enrichment for DEGs showed that they were primarily associated with chromosome activity, the cell cycle, and mitosis (Figures 3A, C). The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated significant enrichment in pathways related to cell cycle, systemic lupus erythematosus, retinol metabolism, motor proteins, neutrophil extracellular trap formation, alcoholism, and amoebiasis (Figure 3B). Figure 3D illustrates the DEGs enriched in corresponding pathways. Overall, the functions of DEGs are significantly related to the cell cycle.

Figure 3
www.frontiersin.org

Figure 3. Enrichment plots from GO-KEGG analysis. (A) GO enrichment semantic similarity matrix. (B) Stacked chart of GO terms; (C) KEGG enrichment significance and the descriptions of each corresponding pathway. (D) DEGs enriched under the corresponding pathways.

PPI network analysis

PPI network analysis is essential for studying the pathway involved in tumor development. To better understand the potential biological functions of CSE1L in BLCA, gene correlations were assessed using Pearson correlation analysis. Using the STRING database, we constructed a PPI network. Based on the intersecting genes, we performed Gene Ontology Biological Process (GO-BP) enrichment analysis and classified their functions into seven categories for PPI network annotation. The inner circle in the PPI network represents the regulation of gene expression, while the outer circle indicates the different functions of the gene/protein (Figure 4).

Figure 4
www.frontiersin.org

Figure 4. Protein-protein interaction (PPI) network of the intersecting genes.

Correlation between CSE1L expression and immune infiltration

First, the CIBERSORT method was used to assess the infiltration of 22 different immune cell types in BLCA. The Wilcoxon method was employed to explore differences in immune cell infiltration between the high- and low-CSE1L expression groups. Our analysis revealed that CSE1L expression primarily affected T cells, macrophages, and mast cells (Figure 5A). We further analyzed the correlation between CSE1L expression and immune infiltration. As illustrated in Figure 5B, CSE1L expression was positively correlated with the CD4 memory-activated T cells, M0 macrophages, M1 macrophages, and resting NK cells, while it was negatively correlated with the infiltration level of regulatory Tregs and plasma cells. Next, we used ssGSEA to determine the infiltration of different immune cell types in BLCA. Spearman analysis was conducted to investigate the relationship between CSE1L expression and immune cell infiltration. Our analysis revealed that CSE1L expression was negatively correlated with the pDCs and CD56bright NK cells and positively correlated with the infiltration of Th2 cells (Figure 5C).

Figure 5
www.frontiersin.org

Figure 5. Association analysis of CSE1L gene expression and immune infiltration (A) Comparison of immune cells between high- and low-CSE1L expression groups (*P<0.05, **P<0.01, ***P<0.001.). (B) Immuno-infiltration correlation heat map and association with CSE1L gene expression (C) Correlation between CSE1L expression and Th2, pDC based on ssGSEA.

Drug sensitivity analysis

We conducted a chemosensitivity analysis of the TCGA dataset using the GDSC database through the R software package prophytic. The relationship between CSE1L expression levels and drug sensitivity was explored by calculating the correlation between the half-maximal inhibitory concentration (IC50) of chemotherapeutic drugs and CSE1L expression. A total of 70 drugs with IC50 values significantly correlated with CSE1L expression (p < 0.05) were identified. The results showed that the low CSE1L expression group was more sensitive to MG-132, Palbociclib (PD-0332991), Nutlin-3a, and Selumetinib (AZD6244), while the high CSE1L expression group was more sensitive to S-Trityl-L-cysteine (NSC 83265), Bleomycin, BAY 61-3606, and Cisplatin. These findings have implications for selecting specific medications based on anti-tumor drug sensitivity (Spearman │Cor│> 0.2 and p < 0.05, Figures 6A–H).

Figure 6
www.frontiersin.org

Figure 6. Drug sensitivity analysis. (A–H) There were significant differences in IC50 between the high and low CSE1L expression groups.

CSE1L expression in BLCA tissues and cell lines

We collected tumor tissues and adjacent normal tissues from six pairs of BLCA patients for qRT-PCR and Western blot validation. The results revealed a significantly elevated expression level of CSE1L mRNA in the tumor tissues (Figures 7A, B). This finding was similarly confirmed in BLCA cell lines. Compared to the normal urothelial cell line SV-HUC1, CSE1L mRNA expression was significantly elevated in all BLCA cell lines tested (Figure 7C: 5637, P < 0.01; T24, P < 0.001).

Figure 7
www.frontiersin.org

Figure 7. mRNA expression levels of CSE1L in BLCA tissues and cell lines. (A, B) mRNA expression levels of CSE1L in BLCA tissues; (C) mRNA expression levels of CSE1L in BLCA cell lines. (**p < 0.01, ***P < 0.001).

CSE1L exerted a promoting effect on BLCA in vitro

To explore the biological functions of CSE1L in BLCA, we transfected 5637 and T24 cell lines with siR-CSE1L and a negative control (NC) separately. Knockdown efficiency was validated by qRT-PCR, which demonstrated a significant decrease in CSE1L expression at the mRNA level (Figure 8D). We then assessed the impact of CSE1L on the migration and invasion of BLCA cells using wound healing and transwell invasion assay. The wound healing assay indicated a significant reduction in the migration ability of 5637 and T24 cells after CSE1L knockdown (Figures 8A, E). Similarly, the transwell invasion assay showed a decrease in the invasive ability of BLCA cells as CSE1L expression was downregulated (Figures 8B, F). These results indicate that CSE1L promotes the migration and invasion of BLCA cells in vitro. Additionally, we investigated the effect of CSE1L on the proliferation of BLCA cells using CCK-8 assay. The CCK-8 assay showed that silencing CSE1L significantly attenuated the proliferation of 5637 and T24 cells (Figure 8C). In summary, these findings suggest that CSE1L promotes the development of BLCA in vitro.

Figure 8
www.frontiersin.org

Figure 8. CSE1L promotes bladder cancer cells proliferation, migration and invasion in vitro. (A, E) Wound healing assay indicated that knockdown of CSE1L reduced the migratory ability of bladder cancer cells; (B, F) Transwell invasion assay indicated that knockdown of CSE1L reduced the invasive activity of bladder cancer cells; (C) CCK-8 assay indicated that knockdown of CSE1L reduced the proliferation ability of bladder cancer cells; (D) The mRNA expression level of CSE1L measured by qRT-PCR after cell transfection. *p < 0.05; **p < 0.01; ***p < 0.001.

Discussion

Bladder cancer is the fifth most common cancer worldwide, characterized by a high incidence, poor prognosis, and a significant likelihood of recurrence. Currently, surgical resection is the most effective treatment (10). However, due to the high recurrence rate and substantial heterogeneity of bladder tumors, there is an urgent need for new molecular markers to predict patient survival in this disease.

In 1995, Brinkmann et al. first isolated a CSE1L-DNA fragment in breast cancer cells (11). CSE1L is the human homolog of the yeast gene CSE1, encoding a protein distributed in the cell nucleus and cytoplasm (12). The gene is located on 20q13, a site frequently amplified in various cancers and associated with genetic instability (13). CSE1L has been found to be involved in apoptosis, proliferation, survival, nuclear-cytoplasmic transport, and cancer metastasis (6, 14, 15). Additionally, CSE1L is a novel micro-vesicle membrane protein that may serve as a potential target for the development of efficient antibody-drug conjugates (ADCs) for cancer therapy (16).

As a key factor in the nuclear transport pathway, CSE1L participates in nucleocytoplasmic transport, a crucial process in tumor growth and development. Consequently, it may have applications in clinical diagnosis and treatment. Moreover, CSE1L expression is associated with tumor progression in various types of cancers. However, the clinical significance and biological functions of CSE1L in bladder cancer remain unclear. In this study, CSE1L was found to promote the initiation and progression of tumors. Evidence from the Kaplan-Meier plotter and nomogram analyses indicated that CSE1L possesses both diagnostic and prognostic value in distinguishing bladder cancer patients from healthy individuals. Furthermore, it was revealed that high CSE1L expression is associated with poor prognosis.

Based on the differentially expressed genes we screened and CSE1L-related genes, a PPI network was constructed. GO-BP enrichment analysis was performed on the intersecting genes, categorizing their functions into seven group: cell composition, cell cycle, enzymatic activity, metabolize, protein modifications, signaling and others. Genes that primarily interact with CSE1L include AURKA, NUP155, BIRC5, and TPX2. Interestingly, many of these genes are associated with unfavorable patient prognosis. Our study also revealed that CSE1L expression in bladder cancer is correlated with the type and density of infiltrating immune cells. The degree of infiltration of various immune cells, including T cells and macrophages, was significantly correlated with the CSE1L expression. CSE1L is believed to play a role in nucleus-to-cytoplasmic transport, chromosome separation during mitosis, proliferation, and apoptosis. In our study, the decrease in immune infiltration of regulatory T cells may be due to the high expression of CSE1L in these cells, leading to their apoptosis. The relationship between CSE1L expression and immune regulatory cells warrants further investigation.

Previous studies have shown that the expression of CSE1L is positively correlated with the immune checkpoint molecules PD-L1 (CD274) and PDL2 (PDCD1LG2). Progression-free survival (PFS) was significantly shorter in patients with high CSE1L expression compared to those with low expression, suggesting that CSE1L expression may affect the therapeutic efficacy of PD-1 monoclonal antibodies (17). Studies have also indicated that PD-L1 protein levels are reduced following CSE1L silencing (18). Additionally, Lee et al. found that serum phospho-CSE1L could be used for early detection of tumor-targeted therapy efficacy and to monitor secondary resistance in mouse tumor xenograft models. The CSE1L high expression group was more sensitive to S-Trityl-L-cysteine (NSC 83265), Bleomycin, BAY 61-3606, and Cisplatin, while the CSE1L low expression group was more sensitive to MG-132, Palbociclib (PD-0332991), Nutlin-3a, and Selumetinib (AZD6244). Notably, Nutlin-3a, which inhibits the MDM2-p53 interaction and stabilizes p53 protein, induces autophagy and apoptosis. In future studies, we will verify these findings using a larger cancer cell panel and in vivo animal models.

In this study, a series of experiments was conducted to elucidate the biological functions of CSE1L in bladder tumor cells. We found that silencing CSE1L expression in bladder cancer cells inhibited cell proliferation and reduced cell invasion and migration. CSE1L promotes tumor progression by enhancing cell proliferation and contributing to drug resistance. These findings align with our understanding of CSE1L’s impact on cancer cell growth and are consistent with previous studies.

Multiple potential underlying mechanisms of CSE1L action have been documented in various types of cancer. Studies have proposed that CSE1L interacts with the cAMP/PKA and RAS/ERK signaling pathways in melanoma. According to Jiang et al., the phosphorylation of CSE1L is regulated by extracellular signal-regulated kinase (ERK) (19, 20). Notably, CSE1L promotes melanin formation and melanoma progression by influencing the expression and phosphorylation of CREB (cAMP response element binding protein) and MITF (microphthalmia-associated transcription factor) (20). The PKA signaling pathway induces CREB phosphorylation, thereby prompting MITF expression. Furthermore, the PKA pathway also interacts with the MAPK/ERK pathway, leading to ERK activation. ERK, in turn, interacts with CSE1L, facilitating the phosphorylation of MITF (21). In gastric cancer, CSE1L inhibition has been shown to activate the PI3K/AKT/mTOR and MEK/ERK pathways by downregulation of MITF and GPNMB (15). In pancreatic cancer, CSE1L may regulate proliferation through the AKT signaling pathway (22). In lung cancer, CSE1L interacts with p65 and regulates MAPK signaling (23, 24). CSE1L could also promotes the proliferative and migratory abilities in oral cancer by positively regulating MITF and activating the Akt/mTOR pathway (25). In osteosarcoma cells, CSE1L, a positive regulator of MSH6 protein, is associated with poor patient prognosis. A rescue study suggested that the inhibitory effect of CSE1L knockdown on osteosarcoma cell growth could be reversed by overexpression of MSH6 (26). Current research has investigated whether CSE1L can be used to predict the prognosis of bladder cancer (BLCA) and provide treatment recommendations. However, this study had several limitations. First, due to the lack of additional data sources, the study may be subject to bias. Second, prospective studies are needed to validate our findings. Third, the stability and accuracy of CSE1L as a prognostic marker must be verified using data from other sources. Fourth, the molecular mechanisms by which CSE1L influences the development and progression of BLCA require further exploration through basic experiments.

Conclusions

The expression level of CSE1L is upregulated in BLCA cells and is correlated with tumor progression, poor prognosis, and immune infiltration in bladder cancer. In vitro experiments demonstrated that the proliferation, migration, and invasion of BLCA cells can be inhibited by knocking down CSE1L expression.

Data availability statement

The original contributions presented in the study are included in the article/supplementary materials. Further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by IRB of Beijing Tiantan Hospital, Capital Medical University {KY2024-005-01}. The studies were conducted in accordance with the local legislation and institutional requirements. The human samples used in this study were acquired from primarily isolated as part of your previous study for which ethical approval was obtained. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

RL: Data curation, Formal analysis, Investigation, Resources, Writing – original draft, Writing – review & editing. ZZ: Investigation, Methodology, Validation, Writing – review & editing. JM: Resources, Software, Writing – review & editing. YZ: Conceptualization, Investigation, Project administration, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the National Nature Science Foundation of China (Grant No.82270782).

Acknowledgments

Thanks to Dr. Ning and Dr. Zhai for their technical support and guidance for this study.

Conflict of interest

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

The reviewer WD declared a shared parent affiliation, with no collaboration, with the authors to the handling editor at the time of the review.

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.

References

1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2021) 71:209–49. doi: 10.3322/caac.21660

PubMed Abstract | Crossref Full Text | Google Scholar

2. Babjuk M, Burger M, Capoun O, Cohen D, Comperat EM, Dominguez Escrig JL, et al. European association of urology guidelines on non-muscle-invasive bladder cancer (Ta, T1, and carcinoma in situ). Eur Urol. (2022) 81:75–94. doi: 10.1016/j.eururo.2021.08.010

PubMed Abstract | Crossref Full Text | Google Scholar

3. Witjes JA, Bruins HM, Cathomas R, Comperat EM, Cowan NC, Gakis G, et al. European association of urology guidelines on muscle-invasive and metastatic bladder cancer: summary of the 2020 guidelines. Eur Urol. (2021) 79:82–104. doi: 10.1016/j.eururo.2020.03.055

PubMed Abstract | Crossref Full Text | Google Scholar

4. Cole CN, Scarcelli JJ. Transport of messenger RNA from the nucleus to the cytoplasm. Curr Opin Cell Biol. (2006) 18:299–306. doi: 10.1016/j.ceb.2006.04.006

PubMed Abstract | Crossref Full Text | Google Scholar

5. Beck M, Hurt E. The nuclear pore complex: understanding its function through structural insight. Nat Rev Mol Cell Biol. (2017) 18:73–89. doi: 10.1038/nrm.2016.147

PubMed Abstract | Crossref Full Text | Google Scholar

6. Behrens P, Brinkmann U, Wellmann A. CSE1L/CAS: its role in proliferation and apoptosis. Apoptosis. (2003) 8:39–44. doi: 10.1023/A:1021644918117

PubMed Abstract | Crossref Full Text | Google Scholar

7. Jiang MC. CAS (CSE1L) signaling pathway in tumor progression and its potential as a biomarker and target for targeted therapy. Tumour Biol. (2016) 37:13077–90. doi: 10.1007/s13277-016-5301-x

PubMed Abstract | Crossref Full Text | Google Scholar

8. Nagashima S, Maruyama J, Honda K, Kondoh Y, Osada H, Nawa M, et al. CSE1L promotes nuclear accumulation of transcriptional coactivator TAZ and enhances invasiveness of human cancer cells. J Biol Chem. (2021) 297:100803. doi: 10.1016/j.jbc.2021.100803

PubMed Abstract | Crossref Full Text | Google Scholar

9. Chen B, Khodadoust MS, Liu CL, Newman AM, Alizadeh AA. Profiling tumor infiltrating immune cells with CIBERSORT. Methods Mol Biol. (2018) 1711:243–59. doi: 10.1007/978-1-4939-7493-1_12

PubMed Abstract | Crossref Full Text | Google Scholar

10. Tran L, Xiao JF, Agarwal N, Duex JE, Theodorescu D. Advances in bladder cancer biology and therapy. Nat Rev Cancer. (2021) 21:104–21. doi: 10.1038/s41568-020-00313-1

PubMed Abstract | Crossref Full Text | Google Scholar

11. Brinkmann U, Brinkmann E, Gallo M, Pastan I. Cloning and characterization of a cellular apoptosis susceptibility gene, the human homologue to the yeast chromosome segregation gene CSE1. Proc Natl Acad Sci USA. (1995) 92:10427–31. doi: 10.1073/pnas.92.22.10427

PubMed Abstract | Crossref Full Text | Google Scholar

12. Scherf U, Pastan I, Willingham MC, Brinkmann U. The human CAS protein which is homologous to the CSE1 yeast chromosome segregation gene product is associated with microtubules and mitotic spindle. Proc Natl Acad Sci USA. (1996) 93:2670–4. doi: 10.1073/pnas.93.7.2670

PubMed Abstract | Crossref Full Text | Google Scholar

13. Tanaka T, Ohkubo S, Tatsuno I, Prives C. hCAS/CSE1L associates with chromatin and regulates expression of select p53 target genes. Cell. (2007) 130:638–50. doi: 10.1016/j.cell.2007.08.001

PubMed Abstract | Crossref Full Text | Google Scholar

14. Dong Q, Li X, Wang CZ, Xu S, Yuan G, Shao W, et al. Roles of the CSE1L-mediated nuclear import pathway in epigenetic silencing. Proc Natl Acad Sci USA. (2018) 115:E4013–E22. doi: 10.1073/pnas.1800505115

PubMed Abstract | Crossref Full Text | Google Scholar

15. Li Y, Yuan S, Liu J, Wang Y, Zhang Y, Chen X, et al. CSE1L silence inhibits the growth and metastasis in gastric cancer by repressing GPNMB via positively regulating transcription factor MITF. J Cell Physiol. (2020) 235:2071–9. doi: 10.1002/jcp.v235.3

PubMed Abstract | Crossref Full Text | Google Scholar

16. Liao CF, Lin SH, Chen HC, Tai CJ, Chang CC, Li LT, et al. CSE1L, a novel microvesicle membrane protein, mediates Ras-triggered microvesicle generation and metastasis of tumor cells. Mol Med. (2012) 18:1269–80. doi: 10.2119/molmed.2012.00205

PubMed Abstract | Crossref Full Text | Google Scholar

17. Li H, Wang L, Ruan Z, Li X, Yang Y, Fang J, et al. CSE1L as a prognostic biomarker associated with pan cancer immune infiltration and drug sensitivity. Expert Rev Clin Immunol. (2024) 20(9):1113–25. doi: 10.1080/1744666X.2024.2356747

PubMed Abstract | Crossref Full Text | Google Scholar

18. Liu W, Zhou Z, Li Y, Xu J, Shen Y, Luo S, et al. CSE1L silencing impairs tumor progression via MET/STAT3/PD-L1 signaling in lung cancer. Am J Cancer Res. (2021) 11:4380–93.

PubMed Abstract | Google Scholar

19. Lee WR, Shen SC, Shih YH, Chou CL, Tseng JT, Chin SY, et al. Early decline in serum phospho-CSE1L levels in vemurafenib/sunitinib-treated melanoma and sorafenib/lapatinib-treated colorectal tumor xenografts. J Transl Med. (2015) 13:191. doi: 10.1186/s12967-015-0553-6

PubMed Abstract | Crossref Full Text | Google Scholar

20. Lee WR, Shen SC, Wu PR, Chou CL, Shih YH, Yeh CM, et al. CSE1L Links cAMP/PKA and Ras/ERK pathways and regulates the expressions and phosphorylations of ERK1/2, CREB, and MITF in melanoma cells. Mol Carcinog. (2016) 55:1542–52. doi: 10.1002/mc.v55.11

PubMed Abstract | Crossref Full Text | Google Scholar

21. Levy C, Khaled M, Fisher DE. MITF: master regulator of melanocyte development and melanoma oncogene. Trends Mol Med. (2006) 12:406–14. doi: 10.1016/j.molmed.2006.07.008

PubMed Abstract | Crossref Full Text | Google Scholar

22. Zhang X, Zhang X, Mao T, Xu H, Cui J, Lin H, et al. CSE1L, as a novel prognostic marker, promotes pancreatic cancer proliferation by regulating the AKT/mTOR signaling pathway. J Cancer. (2021) 12:2797–806. doi: 10.7150/jca.54482

PubMed Abstract | Crossref Full Text | Google Scholar

23. Lin HC, Li J, Cheng DD, Zhang X, Yu T, Zhao FY, et al. Nuclear export protein CSE1L interacts with P65 and promotes NSCLC growth via NF-kappaB/MAPK pathway. Mol Ther Oncolytics. (2021) 21:23–36. doi: 10.1016/j.omto.2021.02.015

PubMed Abstract | Crossref Full Text | Google Scholar

24. Duan L, Tadi MJ, Maki CG. CSE1L is a negative regulator of the RB-DREAM pathway in p53 wild-type NSCLC and can be targeted using an HDAC1/2 inhibitor. Sci Rep. (2023) 13:16271. doi: 10.1038/s41598-023-43218-3

PubMed Abstract | Crossref Full Text | Google Scholar

25. Wang YS, Peng C, Guo Y, Li Y. CSE1L promotes proliferation and migration in oral cancer through positively regulating MITF. Eur Rev Med Pharmacol Sci. (2020) 24:5429–35. doi: 10.26355/eurrev_202005_21327

PubMed Abstract | Crossref Full Text | Google Scholar

26. Cheng DD, Lin HC, Li SJ, Yao M, Yang QC, Fan CY. CSE1L interaction with MSH6 promotes osteosarcoma progression and predicts poor patient survival. Sci Rep. (2017) 7:46238. doi: 10.1038/srep46238

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: CSE1L, bladder urothelial carcinoma, prognosis, therapeutic target, immune infiltration

Citation: Liu R, Ma J, Zhang Y and Zhou Z (2024) Integrated transcriptome analysis of CSE1L regarding poor prognosis and immune infiltration in bladder urothelial carcinoma and experimental verification. Front. Immunol. 15:1449251. doi: 10.3389/fimmu.2024.1449251

Received: 14 June 2024; Accepted: 19 September 2024;
Published: 04 October 2024.

Edited by:

Pengpeng Zhang, Nanjing Medical University, China

Reviewed by:

Baosheng Li, Shandong Cancer Hospital, China
Shuanghu Yuan, Anhui Provincial Hospital, China
Wei Deng, Capital Medical University, China

Copyright © 2024 Liu, Ma, Zhang and Zhou. 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.

*Correspondence: Yong Zhang, doctorzzb@163.com; Zhongbao Zhou, 346894383@qq.com

Disclaimer: 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.