Cuproptosis-related lncRNAs modulate the prognosis of MIBC by regulating the expression pattern of immunosuppressive molecules within the tumor microenvironment

doi:10.21203/rs.3.rs-3132752/v1

Download PDF

Research Article

Cuproptosis-related lncRNAs modulate the prognosis of MIBC by regulating the expression pattern of immunosuppressive molecules within the tumor microenvironment

https://doi.org/10.21203/rs.3.rs-3132752/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Cuproptosis-related gene and long non-coding RNA (lncRNA) modulation of cancer regulation is well-established. This investigation aimed to elucidate the prognostic implications of cuproptosis-associated lncRNAs in muscle-invasive bladder cancer (MIBC). Employing The Cancer Genome Atlas (TCGA) and IMvigor210 cohorts, bioinformatics and statistical analyses probed the prognostic relevance of cuproptosis-related lncRNAs. Co-expression analysis revealed tight associations between lncRNA expression and cuproptosis-linked genes, with 13 cuproptosis-related lncRNAs found to correlate with MIBC prognosis. Lasso regression identified a six-lncRNA prognostic signature, enabling patient stratification into high- and low-risk categories. Tissue validation substantiated differential expression of FAM13A-AS1, GHRLOS, LINC00456, OPA1-AS1, RAP2C-AS1, and UBE2Q1-AS1 between MIBC tumor and normal tissues. Comparative analyses of tumor microenvironments and immune profiles between risk groups disclosed elevated immunosuppressive molecule expression, including programmed cell death-1 (PD-L1) and T-cell immunoglobulin-3 (TIM-3), in high-risk individuals. These findings suggest that cuproptosis-related lncRNAs modulate immune responses within the tumor microenvironment, thereby influencing MIBC tumorigenesis and progression. Further exploration is warranted to unveil novel therapeutic targets for MIBC based on the expression patterns of cuproptosis-related lncRNAs and their impact on immune responses in the tumor microenvironment.

cuproptosis-related lncRNAs

immunosuppressive molecules

MIBC

prognostic

risk scores

Bladder carcinoma (BLCA) ranks among the top ten most prevalent and costly cancers globally, accounting for approximately 550,000 new cases annually and significant morbidity and mortality(1, 2). Muscle-invasive bladder cancer (MIBC) is initially diagnosed in 25% of patients, while 10%-20% of non-muscle invasive bladder cancer (NMIBC) cases progress to MIBC(3). Despite advancements in surgical and treatment modalities, MIBC prognosis remains poor, with a 5-year post-therapy survival rate of 40%-60%(4, 5). Investigating MIBC development mechanisms and identifying prognostic biomarkers are vital for enhancing patient outcomes and tailoring postoperative treatment strategies.

Cuproptosis, a recently discovered form of programmed cell death, involves copper, an essential cofactor in enzymatic activity and physiological homeostasis(6–9). Previous studies have shown that c opper induces tumor cell death, inhibits tumor progression, contributes to proteotoxic stress, impedes proteasomes, and promotes anti-angiogenesis(10–12). However, few studies have explored the BLCA microenvironment or its biological functions, warranting further investigation into cuproptosis' relationship with BLCA development and progression. In the meantime, long non-coding RNA (lncRNA) dysregulation is associated with various tumor types(13, 14), but the complete role of cuproptosis-related gene regulatory factors in lncRNA regulation remains unclear. Uncovering the mechanisms by which cuproptosis-related gene modification influences MIBC development may reveal potential therapeutic targets.

Using The Cancer Genome Atlas (TCGA) and IMvigor210 cohorts, we analyzed MIBC patient data to determine the prognostic significance of cuproptosis-related lncRNA. Our findings indicated that six cuproptosis-related lncRNAs possess prognostic value in MIBC patients, and we constructed a prognostic signature based on their ability to predict overall survival. Crucially, our study suggests that cuproptosis-related lncRNA modulates immunosuppressive molecule expression, impacting MIBC prognosis.

Tissue sample collection and RT-qPCR analysis

Ten pairs of fresh tumor tissue samples from MIBC patients and paired normal tissue samples were collected at the Xinhua Hospital Affiliated Shanghai Jiao Tong University School of Medicine. This research was approved by the Medical Ethics Committee of Xinhua Hospital of Affiliated Shanghai Jiao Tong University School of Medicine and informed consent was obtained from all patients. Total RNA was isolated from tissues using TRIzol reagent (Invitrogen, CA, USA) following the manufacturer’s protocol. Then, cDNA was used as a qPCR template to quantify genes using TB Green® Premix Ex TaqTM II (Takara, Shanghai, China). The relative gene expression was normalized to β-actin according to the cycle threshold (2 − ΔCt) method. The experiment was performed in triplicate, and the primer sequences employed for analysis are listed in Table 1.

Table 1

Primer sequences of lncRNA
FAM13A-AS1	F: 5′-CAAATATGGGTAAGGAGG-3′
FAM13A-AS1	R: 5′-GTTTAGAACTATGAGGGACT-3′
GHRLOS	F: 5'-TGGAAACTCCCCTAGCCACA-3′
GHRLOS	R: 5'-GCATCTCTCCTCTGTTCCGT-3′
RAP2C-AS1	F: 5'-ACTTAGCCGTGCCTGACAAA-3′
RAP2C-AS1	R: 5'-GCTCCAAAAAGGCACCCTTG-3′
UBE2Q1-AS1	F: 5'-AGTCTCTCCACCAGCCAGAT-3′
UBE2Q1-AS1	R: 5'-TGACAGAGACAGCCATTTTGC-3′
LINC00456	F: 5'-TGAGCCAAGTCCAGTGGAAA-3′
LINC00456	R: 5'- AAGTCTTGTCGCACTGTAGCA-3′
OPA1-AS1	F: 5'-TCCACAGCCATGGAATTGAGAA-3′
OPA1-AS1	R: 5'-ACGGCGTTTAGAGCAACAGA-3′

MIBC Datasets and Samples

The RNA-seq transcriptome data of patients with BLCA and the corresponding clinical data were downloaded from the TCGA data portal (https://portal.gdc.cancer.gov/). Finally, 396 MIBC patients were included for further analysis based on the simultaneous availability of mRNA expression profile data and complete clinical information. The IMvigor210 dataset was obtained from http://research-pub.gene.com/IMvigor210CoreBiologies/. The detailed clinical patient characteristics are summarized in Table 2. Patients with locally advanced or metastatic uroepithelial carcinoma were included in the study and further analyzed.

Table 2

Clinical characteristics of BLCA patients
		TCGA		IMvigor210
	Type	Test	Train	IMvigor210
Age	<=65	77	81	/
	> 65	121	117	/
Gender	Female	52	51	35
	Male	146	147	133
Grade	High	186	188	/
	Low	11	9	/
	unknown	1	1	/
Stage	I	0	0	56
	II	66	61	41
	III	67	71	35
	IV	65	66	36

Cuproptosis-related genes detection

Expression matrixes of 19 cuproptosis-related genes were extracted from the TCGA databases based on previous publications, including NFE2L2, NLRP3, ATP7B, ATP7A, SLC31A1, FDX1, LIAS, LIPT1, LIPT2, DLD, DALAT, PDHA1, PDHB, MTF1, GLS, CDKN2A, DBT, GCSH, DLST.

Identification of cuproptosis-related lncRNAs

The long non-coding RNA annotation file of Genome Reference Consortium Human Build 38 (GRCh38) was obtained from the GENCODE website for the annotation of lncRNAs in both the TCGA dataset and the Imvigor210 dataset. Pearson correlation analysis was conducted in the TCGA and IMvigor210 datasets to identify cuproptosis-related lncRNAs. Following this, 71 survival-associated cuproptosis-related lncRNAs were selected for further investigation.

Establishment of the prognostic model of cuproptosis-related lncRNAs

The least absolute shrinkage and selection operator (LASSO) Cox regression algorithm was used to establish a prognostic model based on the normalized expression matrix of candidate prognostic 6 cuproptosis-related lncRNAs. The patient's risk score was calculated according to the expression level of each gene and its corresponding regression coefficient. The formula was established as follows: risk score = gene's expression × corresponding coefficient. MIBC samples were divided into high-risk and low-risk groups based on the median value of the risk score. Ttryyyyhhe Kaplan-Meier (KM) curve, the receiver operating characteristic (ROC) and the principal component (PCA) were applied to evaluate the difference in the survival status of patients between the two groups. Finally, a nomogram model was established to enhance the predictive efficacy of our model in the TCGA dataset.

Screening and functional enrichment of differential genes between the two groups

The differential expression genes (DEGs) were filtrated by “Limma” packages between-group (high‐risk vs low‐risk). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses based on the DEGs (|log2FC| ≥ 1, FDR < 0.05) between the high-risk and low-risk groups.

Estimation of infiltrating cells in TME in MIBC samples

Single-sample gene set enrichment analysis (ssGSEA) was used to calculate the infiltration score of 16 immune cells and the activity of 13 immune-related pathways.

Calculation of Tumor Mutation Burden Scores

Tumor mutational burden (TMB) data from BLCA samples downloaded from TCGA were analyzed using the "maftools" R package. The waterfall plot displayed the relationship between the risk signature model and TMB for BLCA patients.

Statistical analysis

Kaplan-Meier analysis and log-rank test were used to compare OS between groups. Univariate and multivariate Cox regression analyses were used to determine independent predictors of OS. The chi-square test was used for the comparison of proportional differences. All statistical analysis is carried out using R software.

Identification of cuproptosis-related lncRNA in MIBC patients

To ascertain the biological function of cuproptosis-related lncRNAs in the development of MIBC, we identified these lncRNAs from the TCGA dataset and the IMvigor210 dataset. The study flowchart is depicted in Fig. 1A. We identified 71 common cuproptosis-related lncRNAs in both the TCGA and IMvigor210 datasets. Subsequently, univariate Cox regression analysis revealed 13 cuproptosis-related lncRNAs that were significantly associated with survival prognosis (Fig. 1B). We systematically examined the expression profiles of these 13 prognostic-related cuproptosis-related lncRNAs between MIBC and adjacent normal tissue pairs, finding that their expression patterns varied between tumor and normal tissues. Some cuproptosis-related lncRNAs were highly expressed in tumor tissues, while others were more highly expressed in normal tissues (Fig. 1C). Our results support the notion that cuproptosis-related lncRNAs may constitute one of the components influencing the initiation and progression of MIBC.

Construction and validation of prognostic signatures for cuproptosis-related lncRNAs

To explore the contribution of cuproptosis-related lncRNAs in the prognosis of MIBC patients, 396 MIBC patients were randomly assigned to TCGA training cohort and validation cohort at a ratio of 5:5. Meanwhile, we performed LASSO regression analysis based on the expression values of 13 cuproptosis-related lncRNAs previously screened in the TCGA training cohort, and finally selected 6 cuproptosis-related lncRNAs, including FAM13A-AS1, GHRLOS, LINC00456, OPA1-AS1, RAP2C-AS1, and UBE2Q1-AS1 (Fig. 1D-F).

Subsequently, patients were divided into either the high-risk group or the low-risk group according to the median value of the risk score of prognostic cuproptosis-related lncRNAs. In both the training and validation cohorts, overall survival (OS) was longer in the low-risk group than in the high-risk group (Fig. 2A-C). Moreover, the survival rate was almost double in the low-risk group compared with that in the high-risk group (Fig. 2D). The PCA analyses indicated that patients in different risk groups were distributed in three directions (Fig. 2E). In addition, the predictive value of the signature for OS and clinical characteristics was evaluated using time-dependent ROC curves, showing that the predictive value of the risk signature model was higher than that of clinical traits (Fig. 2F and 2G). For further evaluating the relationship between the level of lncRNAs expression and clinical features, the heatmap was used to display the relationship between the expression profiles of 6 cuproptosis-related lncRNAs and their clinical characteristics in all MIBC patients (Fig. 2H).

Construction of the cuproptosis-related lncRNAs Nomogram

Univariate and multivariate Cox regression analysis revealed that the risk score of cuproptosis-related lncRNAs was a substantial independent prognostic factor for overall survival in the TCGA dataset (Fig. 3A and 3B). To improve the accuracy of the prediction model, we constructed a nomogram for MIBC based on the risk signature and clinical characteristics. The total score was the sum of the individual scores calculated based on the use of the nomogram, and the OS of the patient gradually diminish as the score increased. Figure 3C illustrates an example of employing the nomogram to predict the survival probability for a randomly selected patient. ROC analysis demonstrated that the nomogram surpassed other factors in prognosticating MIBC outcomes (Fig. 3D). Additionally, the nomogram exhibited high concordance between the predicted and observed survival probabilities (Fig. 3E).

Functional annotation of cuproptosis-related lncRNAs risk score

To investigate the function of cuproptosis-related lncRNAs, the “Limma” package was applied to analyze DEGs between the high-risk and low‐risk groups, and 445 DEGs were filtrated. Interestingly, GO enrichment showed that DEGs were significantly enriched in many areas related to the regulation of immune responses, including antigen binding, antibody production, humoral immune response, and lymphocyte-mediated immunity (Fig. 4A). KEGG pathway analyses also indicated that the cytokine-cytokine receptor interaction, IL-17signaling pathway, and TNF signaling pathway were enriched in these DEGs (Fig. 4B).

GO and KEGG analysis displayed that DEGs in the risk score were enriched in immune interaction pathways. We further studied the relationship between the risk score of cuproptosis-related lncRNAs and tumor-infiltrating immune cells. Intriguingly, the results showed that both pro-tumor-related immune factors and anti-tumor-related immune parameters were upregulated in the high-risk group (Fig. 4C and 4D). Subsequently, we analyzed the expression of tumor immunosuppressive checkpoints between the two groups and the majority of tumor immunosuppressive checkpoints were detected to be highly expressed in the high-risk group (Fig. 4E). This observation may explain the limited clinical efficacy of single immune checkpoint inhibitors.

The relationship between cuproptosis-related lncRNAs risk scores and the tumor mutation burden

Previously reported that the TMB can influence the infiltration of tumor immune cells(15). We explored the expression pattern of TMB in the TME between the high- and low-risk groups. The results depicted that the mutation rates of the top 20 most substantially variable genes were significantly distinct between the high- and low-risk groups, as shown in Fig. 5A and 5B. Upon further analysis, we observed that a substantial number of patients in the low-risk group exhibited low TMB, and their overall survival (OS) was significantly compromised, indicating that low TMB adversely affected patient OS. This observation may also provide insight into why some patients in the previous low-risk group did not exhibit prolonged OS (Fig. 5C and 5D). Finally, we integrated the risk signature, patient survival status, tumor stage, molecular subtype, and TMB, identifying a robust correlation among these factors (Fig. 5E).

Validation of cuproptosis-related lncRNAs and association with immunotherapy response in IMvigor210 cohort

Our prior investigations established a robust association between risk scores and immune checkpoints. To validate the constructed model, BLCA patients from the IMvigor210 cohort were selected. Employing the same formula, a corresponding risk score was calculated for each patient within the IMvigor210 cohort. The complete response (CP) and stable disease (SD) rates were higher in the low-risk score group compared to the high-risk score group (Fig. 6A). Moreover, overall survival (OS) was prolonged in the low-risk score group relative to the high-risk score group (Fig. 6B). Regrettably, the scores of TCGA subtypes did not exhibit significant differences, but the risk scores of patients with CR and SD were notably lower (Fig. 6C and 6D). Utilizing the Lund subtype, the risk scores were highest in the genomically unstable subtype and lowest in the infiltrated subtype (Fig. 6E and 6F).

Relative mRNA expression of cuproptosis-related lncRNAs

The RT-qPCR validation results from MIBC patient tissues revealed differential expression of FAM13A-AS1, GHRLOS, LINC00456, OPA1-AS1, RAP2C-AS1, and UBE2Q1-AS1 between normal and tumor samples (Fig. 7A–F). Corresponding P-values for FAM13A-AS1, GHRLOS, LINC00456, and OPA1-AS1 were 0.0098, 0.0273, 0.0137, and 0.0273, respectively. Although UBE2Q1-AS1 expression did not display significant differences, a trend was observed between tumor and normal tissues, potentially warranting sample size expansion for further investigation (Fig. 7F). Collectively, these validation results suggest that most lncRNAs in the model exhibit distinct expression patterns, and additional samples may be required for comprehensive validation.

The high prevalence of BLCA poses a significant threat to human health and imposes a substantial economic burden. As the tumor progresses to MIBC, the specific survival rate declines to 65–78%(5, 16). necessitating further exploration of the biological mechanisms underlying MIBC development, particularly at the genetic level.

Cuproptosis, a novel copper-dependent cell death modality(8), is regulated by copper ions binding to lipid-acylated components of the tricarboxylic acid cycle (TCA) in mitochondrial respiration, subsequently causing proteotoxic stress and cell death(17, 18). Numerous cuproptosis-related genes have been implicated in cancer progression(19, 20). LncRNAs have emerged as crucial regulators in tumorigenesis by modulating nuclear structure, transcription, mRNA stability, translation, and post-translational modifications(21, 22). It has been reported that cuproptosis-related genes can regulate the stability of lncRNA to promote tumor progression. Exploring the expression patterns of cuproptosis-related lncRNAs may provide us with new ideas and targets for diagnosis, treatment, and prognosis assessment in MIBC.

In this study, we investigated the expression and functions of cuproptosis-related lncRNAs in MIBC and BLCA patients using TCGA and IMvigor210 databases. LASSO Cox regression analysis identified six lncRNAs (FAM13A-AS1, GHRLOS, LINC00456, OPA1-AS1, RAP2C-AS1, and UBE2Q1-AS1) with independent prognostic values for MIBC. A risk model based on these lncRNAs classified MIBC patients into high- and low-risk groups, with the latter exhibiting significantly better survival outcomes. The risk score, an independent prognostic factor, demonstrated good predictive ability for MIBC patients' clinical characteristics.

Cuproptosis-related lncRNAs may act as oncogenes or tumor suppressors. Our study identified 6 cuproptosis-related lncRNAs which had prognostic values. FAM13A-AS1 is an important tumor suppressor factor(23, 24). Previous research found that the upregulation of FAM13A-AS1 inhibits the proliferation, migration, and invasion of breast cancer cell lines by targeting the miRNA-205-3p/DDI2 axis (25). GHRLOS expression levels in colorectal cancer tissues were significantly lower than in matched normal tissues, and decreased GHRLOS expression was significantly associated with the occurrence of lymph node metastasis and distant metastasis in colorectal cancer(26). There are relatively few reports on UBE2Q1-AS1, and studies have shown that UBE2Q1-AS1 may contribute to GC development(27). RAP2C-AS1 is widely expressed in tissues and organs. Studies have indicated that RAP2C-AS1 regulates the occurrence and development of esophageal cancer through the lncRNA-miRNA-mRNA axis(28). Our RT-qPCR results confirmed the low expression of RAP2C-AS1 in tumor tissues, which may imply that it is a tumor suppressor gene.

To elucidate the role of cuproptosis-associated lncRNAs in MIBC, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses based on differentially expressed genes (DEGs) between distinct risk groups. Unexpectedly, we discovered that numerous cytokine activities, lymphocyte-mediated immune responses, and cytokine-cytokine receptor interactions were enriched. It is plausible to hypothesize that cuproptosis-related lncRNAs are intimately linked to tumor immune regulation, which aligns with prior findings(29, 30). Delving further into the impact of cuproptosis-associated lncRNAs on immune cell infiltration, we observed that various tumor-associated immune cells, including Treg cells, mast cells, and follicular helper T cells, were highly infiltrated in the high-risk group, alongside an enrichment of numerous anti-tumor immune-related signaling molecules. Intriguingly, we identified a substantial presence of immunosuppressive molecules, such as PD-1, IDO1, CTLA4, and PD-L1, in high-risk patients. Previous research has substantiated that elevated PD-L1 expression in MIBC patient tumor tissues enhances the efficacy of anti-PD-L1 therapy(31). I Our findings indicate that high-risk patients in our cohort may exhibit improved responses to anti-PD-L1 immunotherapy. We subsequently validated our hypothesis utilizing the IMvigor210 dataset of BLCA patients undergoing anti-PD-L1 treatment. The constructed predictive model also demonstrated strong prognostic accuracy for anti-PD-L1 treatment response; however, we were surprised to find that patient survival in the high-risk group did not improve following anti-PD-L1 therapy. This may be attributable to the elevated expression of immunosuppressive molecules within the high-risk group, rendering nanotherapeutic anti-immune treatments ineffective (32). These observations suggested that future therapeutic approaches may be directed toward combination therapies(33, 34).

Nonetheless, this study is not without limitations, including the restricted sample size and the absence of methylation level analysis for the candidate lncRNAs. To comprehensively elucidate the regulatory mechanisms of candidate lncRNAs in MIBC, future investigations should encompass an expanded sample size and the acquisition of solid tumor specimens. The identified six lncRNAs may offer valuable insights into the functional mechanisms of cuproptosis-related genes in MIBC patients. Moreover, it is imperative to undertake experimental research to further probe the roles of these lncRNAs.

In summary, we identified six prognostic lncRNA markers associated with cuproptosis-related genes and validated their predictive significance in the context of MIBC patient prognosis by leveraging publicly available databases. Our findings suggest that candidate cuproptosis-related lncRNAs may modulate the expression of immunosuppressive molecules, thereby influencing tumor initiation and progression. This investigation has enhanced our comprehension of cuproptosis-related lncRNAs and their mechanistic roles within the tumor microenvironment, potentially unveiling novel therapeutic targets and prognostic biomarkers for further exploration.

Acknowledgments

Thanks to professor Shen HB and doctor Liu HL for their strict guidelines on this article, and thanks to every author who participated in this research. Thanks to the reviewers and editors for their sincere comments.

Disclosure of potential conflicts of interest

No potential conflict of interest was reported by the author(s).

Data availability statement

The dataset supporting the conclusions of this article is available upon reasonable request from Haibo Shen.

Funding

This work was supported by the Science and Technology Commission of Shanghai Municipality (No. 21511902300).

Author contributions

Shen HB, Cao QF, and Liu HL conceived the project and designed experiments. Duan HQ, Shen Y, Yu SG and Wang C collected, organized, and analyzed these data. Xia WM, Xu D, Liu HL and Zhang S helped to collect databases and organized these databases. Duan HQ and Shen HB wrote the manuscript.

Richters A, Aben KKH, Kiemeney L (2020) The global burden of urinary bladder cancer: an update. World J Urol 38(8):1895–1904
Feng RM, Zong YN, Cao SM, Xu RH (2019) Current cancer situation in China: good or bad news from the 2018 Global Cancer Statistics? Cancer Commun (Lond) 39(1):22
Álvarez-Maestro M, Guerrero-Ramos F, Rodríguez-Faba O, Domínguez-Escrig JL, Fernández-Gómez JM (2021) Current treatments for BCG failure in non-muscle invasive bladder cancer (NMIBC). Actas Urol Esp (Engl Ed) 45(2):93–102
Witjes JA, Bruins HM, Cathomas R, Compérat EM, Cowan NC, Gakis G et al (2021) European Association of Urology Guidelines on Muscle-invasive and Metastatic Bladder Cancer: Summary of the 2020 Guidelines. Eur Urol 79(1):82–104
Sanli O, Dobruch J, Knowles MA, Burger M, Alemozaffar M, Nielsen ME et al (2017) Bladder cancer. Nat Rev Dis Primers 3:17022
Duan WJ, He RR (2022) Cuproptosis: copper-induced regulated cell death. Sci China Life Sci 65(8):1680–1682
Tang D, Chen X, Kroemer G (2022) Cuproptosis: a copper-triggered modality of mitochondrial cell death. Cell Res 32(5):417–418
Wang Y, Zhang L, Zhou F (2022) Cuproptosis: a new form of programmed cell death. Cell Mol Immunol 19(8):867–868
Ge EJ, Bush AI, Casini A, Cobine PA, Cross JR, DeNicola GM et al (2022) Connecting copper and cancer: from transition metal signalling to metalloplasia. Nat Rev Cancer 22(2):102–113
Lv H, Liu X, Zeng X, Liu Y, Zhang C, Zhang Q et al (2022) Comprehensive Analysis of Cuproptosis-Related Genes in Immune Infiltration and Prognosis in Melanoma. Front Pharmacol 13:930041
Bian Z, Fan R, Xie L (2022) A Novel Cuproptosis-Related Prognostic Gene Signature and Validation of Differential Expression in Clear Cell Renal Cell Carcinoma. Genes (Basel). ;13(5)
Li SR, Bu LL, Cai L (2022) Cuproptosis: lipoylated TCA cycle proteins-mediated novel cell death pathway. Signal Transduct Target Ther 7(1):158
Bhan A, Soleimani M, Mandal SS (2017) Long Noncoding RNA and Cancer: A New Paradigm. Cancer Res 77(15):3965–3981
Peng WX, Koirala P, Mo YY (2017) LncRNA-mediated regulation of cell signaling in cancer. Oncogene 36(41):5661–5667
Zhang L, Li B, Peng Y, Wu F, Li Q, Lin Z et al (2020) The prognostic value of TMB and the relationship between TMB and immune infiltration in head and neck squamous cell carcinoma: A gene expression-based study. Oral Oncol 110:104943
Woerl AC, Eckstein M, Geiger J, Wagner DC, Daher T, Stenzel P et al (2020) Deep Learning Predicts Molecular Subtype of Muscle-invasive Bladder Cancer from Conventional Histopathological Slides. Eur Urol 78(2):256–264
Cobine PA, Brady DC, Cuproptosis (2022) Cellular and molecular mechanisms underlying copper-induced cell death. Mol Cell 82(10):1786–1787
Zhao J, Guo S, Schrodi SJ, He D (2022) Cuproptosis and cuproptosis-related genes in rheumatoid arthritis: Implication, prospects, and perspectives. Front Immunol 13:930278
Bao JH, Lu WC, Duan H, Ye YQ, Li JB, Liao WT et al (2022) Identification of a novel cuproptosis-related gene signature and integrative analyses in patients with lower-grade gliomas. Front Immunol 13:933973
Zhang Z, Zeng X, Wu Y, Liu Y, Zhang X, Song Z (2022) Cuproptosis-Related Risk Score Predicts Prognosis and Characterizes the Tumor Microenvironment in Hepatocellular Carcinoma. Front Immunol 13:925618
Ahmad S, Abbas M, Ullah M, Aziz M, Beylerli O, Alam M et al (2021) Long non-coding RNAs regulated NF-κB signaling in cancer metastasis: Micromanaging by not so small non-coding RNAs. Seminars in cancer biology.
Su S, Yeh C, Lin C, Hsieh Y, Chuang C, Tang C et al (2021) A novel melatonin-regulated lncRNA suppresses TPA-induced oral cancer cell motility through replenishing PRUNE2 expression. J Pineal Res. :e12760
Sun Z, Jing C, Xiao C, Li T (2020) An autophagy-related long non-coding RNA prognostic signature accurately predicts survival outcomes in bladder urothelial carcinoma patients. Aging 12(15):15624–15637
Rao Y, Liu H, Yan X, Wang J (2020) In Silico Analysis Identifies Differently Expressed lncRNAs as Novel Biomarkers for the Prognosis of Thyroid Cancer. Comput Math Methods Med 2020:3651051
Qiu Z, He L, Yu F, Lv H, Zhou Y (2022) LncRNA FAM13A-AS1 Regulates Proliferation and Apoptosis of Cervical Cancer Cells by Targeting miRNA-205-3p/DDI2 Axis. J Oncol 2022:8411919
Wu S, Liu J, Wang X, Li M, Chen Z, Tang Y (2017) Aberrant Expression of the Long Non-coding RNA GHRLOS and Its Prognostic Significance in Patients with Colorectal Cancer. J Cancer 8(19):4040–4047
Zhang X, Jiang Y, Xie Y, Leng X, Song F (2020) Comprehensive Analysis of lncRNAs Associated with the Pathogenesis and Prognosis of Gastric Cancer. DNA Cell Biol 39(2):299–309
Liu H, Zhang Q, Lou Q, Zhang X, Cui Y, Wang P et al (2020) Differential Analysis of lncRNA, miRNA and mRNA Expression Profiles and the Prognostic Value of lncRNA in Esophageal Cancer. Pathol Oncol Res 26(2):1029–1039
Xu S, Liu D, Chang T, Wen X, Ma S, Sun G et al (2022) Cuproptosis-Associated lncRNA Establishes New Prognostic Profile and Predicts Immunotherapy Response in Clear Cell Renal Cell Carcinoma. Front Genet 13:938259
Han J, Hu Y, Liu S, Jiang J, Wang H (2022) A Newly Established Cuproptosis-Associated Long Non-Coding RNA Signature for Predicting Prognosis and Indicating Immune Microenvironment Features in Soft Tissue Sarcoma. J Oncol 2022:8489387
Mariathasan S, Turley SJ, Nickles D, Castiglioni A, Yuen K, Wang Y et al (2018) TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells. Nature 554(7693):544–548
Hayashi H, Sugawara S, Fukuda Y, Fujimoto D, Miura S, Ota K et al (2022) A Randomized Phase II Study Comparing Nivolumab with Carboplatin-Pemetrexed for EGFR-Mutated NSCLC with Resistance to EGFR Tyrosine Kinase Inhibitors (WJOG8515L). Clin Cancer Res 28(5):893–902
Sharma P, Siefker-Radtke A, de Braud F, Basso U, Calvo E, Bono P et al (2019) Nivolumab Alone and With Ipilimumab in Previously Treated Metastatic Urothelial Carcinoma: CheckMate 032 Nivolumab 1 mg/kg Plus Ipilimumab 3 mg/kg Expansion Cohort Results. J Clin Oncol 37(19):1608–1616
Necchi A, Raggi D, Gallina A, Madison R, Colecchia M, Lucianò R et al (2020) Updated Results of PURE-01 with Preliminary Activity of Neoadjuvant Pembrolizumab in Patients with Muscle-invasive Bladder Carcinoma with Variant Histologies. Eur Urol 77(4):439–446

No competing interests reported.

Download PDF

Version 1

posted

You are reading this latest preprint version

Cuproptosis-related lncRNAs modulate the prognosis of MIBC by regulating the expression pattern of immunosuppressive molecules within the tumor microenvironment

Status:

Version 1

Abstract

Figures

Introduction

Materials and methods

Tissue sample collection and RT-qPCR analysis

MIBC Datasets and Samples

Cuproptosis-related genes detection

Identification of cuproptosis-related lncRNAs

Establishment of the prognostic model of cuproptosis-related lncRNAs

Screening and functional enrichment of differential genes between the two groups

Estimation of infiltrating cells in TME in MIBC samples

Calculation of Tumor Mutation Burden Scores

Statistical analysis

Results

Identification of cuproptosis-related lncRNA in MIBC patients

Construction and validation of prognostic signatures for cuproptosis-related lncRNAs

Construction of the cuproptosis-related lncRNAs Nomogram

Functional annotation of cuproptosis-related lncRNAs risk score

The relationship between cuproptosis-related lncRNAs risk scores and the tumor mutation burden

Validation of cuproptosis-related lncRNAs and association with immunotherapy response in IMvigor210 cohort

Relative mRNA expression of cuproptosis-related lncRNAs

Discussion

Conclusion

Declarations

References

Additional Declarations

Status:

Version 1