Validation of a Metastatic Assay using biopsies to improve risk stratification in patients with prostate cancer treated with radical radiation therapy

Polyploid cancer cells reveal signatures of chemotherapy resistance

Therapeutic resistance in cancer significantly contributes to mortality, with many patients eventually experiencing recurrence after initial treatment responses. Recent studies have identified therapy-resistant large polyploid cancer cells in patient tissues, particularly in late-stage prostate cancer, linking them to advanced disease and relapse. Here, we analyzed bone marrow aspirates from 44 advanced prostate cancer patients and found the presence of circulating tumor cells with increased genomic content (CTC-IGC) was significantly associated with poorer progression-free survival. Single cell copy number profiling of CTC-IGC displayed clonal origins with typical CTCs, suggesting complete polyploidization. Induced polyploid cancer cells from PC3 and MDA-MB-231 cell lines treated with docetaxel or cisplatin were examined through single cell DNA sequencing, RNA sequencing, and protein immunofluorescence. Novel RNA and protein markers, including HOMER1, TNFRSF9, and LRP1, were identified as linked to chemotherapy resistance. These markers were also present in a subset of patient CTCs and are associated with recurrence in public gene expression data. This study highlights the prognostic significance of large polyploid tumor cells, their role in chemotherapy resistance, and the expression of markers tied to cancer relapse, offering new potential avenues for therapeutic development.

Integrating Radiosensitivity Index and Radiation Resistance Related Index Improves Prostate Cancer Outcome Prediction

Purpose

This study aimed to establish a nomogram combining 31-gene signature (31-GS), radiosensitivity index (RSI), and radiation-resistance-related gene index (RRRI) to predict recurrence in prostate cancer (PCa) patients.

Methods and Materials

Transcriptome data of PCa were obtained from gene expression omnibus and the cancer genome atlas to validate the predictive potential of 3 sets of published biomarkers, namely, 31-GS, RSI, and RRRI. To adjust these markers for the characteristics of PCa, we analyzed 4 PCa-associated radiosensitivity predictive indices based on 31-GS, RSI, and RRRI by the Cox analysis and least absolute shrinkage and selection operator regression analysis. Time-dependent receiver operating characteristic curves, decision curve analyses, integrated discrimination improvement, net reclassification improvement and decision tree model construction were used to compare the radiosensitivity predictive ability of these 4 gene signatures. Key modules and associated functions were identified through a weighted gene co-expression network analysis and gene function enrichment analysis. A nomogram was built to improve the recurrence-prediction capability.

Results

We validated and compared the predictive potential of 2 published predictive indices. Based on the 31-GS, RSI, and RRRI, we analyzed 4 PCa-associated radiosensitivity predictive indices: 14Genes, RSI, RRRI, and 20Genes. Among them, 14Genes showed the most promising predictive performance and discriminative capacity. Genes in the key module defined by the 14Genes model were significantly enriched in radiation therapy-related cell death pathways. The area under receiver operating characteristic curve and decision tree variable importance for 14Genes was the highest in the cancer genome atlas and Gene Expression Omnibus Series (GSE) cohorts.

Conclusions

This study successfully established a radiosensitivity-related nomogram, which had excellent performance in predicting recurrence in patients with PCa. For patients who received radiation therapy, the 20Genes and RRRI models can be used to predict recurrence-free survival, whereas 20Genes is more radiation therapy-specific but needs further external validation.

CTHRC1 expresses in cancer-associated fibroblasts and is associated with resistance to anti-androgen therapy in prostate cancer

BACKGROUND

CTHRC1 overexpresses in prostate cancer and is associated with the proliferation, invasion and migration of prostate cancer cells. However, the roles and mechanisms of CTHRC1 expression in prostate cancer prognosis and treatment outcomes remain unknown.

OBJECTIVE

This study aimed to explore the expression and gene function of CTHRC1 in prostate cancer, investigate the prognostic value and potential effect in the treatment of prostate cancer.

METHODS

Bulk and single-cell RNA sequencing analyses were used to evaluate the expression of CTHRC1 in prostate cancer. All data used in the study were obtained from publicly available sources to ensure transparency. Study enrolled 1999 cases of prostate cancer and 531 normal controls. Single-cell RNA sequencing profile included 62,995 cells from seven prostate primary tumors. CTHRC1 expression and prognosis analyses were conducted with these samples and verified by immunohistochemical staining. CIBERSORT algorithm was used to assess the tumor immune infiltrating cells based on bulk mRNA sequencing profiles. Genomics of drug sensitivity in cancer (GDSC) database was used to predict IC50 to anti-androgen therapy (ADT) drugs of the samples.

RESULTS

CTHRC1 expressed in prostate cancer was higher than that in normal prostate tissue, and the expression increased with the progress of prostate cancer. CTHRC1 was the risk factor of progression-free interval (PFI). CTHRC1 was positively correlated with the infiltration of tumor-associated macrophages (TAMs). Myofibroblast-like cancer-associated fibroblasts (myCAFs) were the major CTHRC1 expressers in prostate cancer. TGF-β signaling activated in CTHRC1-positive myCAFs and was involved in TAMs polarization. Biological functions of myCAFs were enriched in hormone response and metabolism. CTHRC1 may regulate androgen receptor signaling through CCN2/CAV1/AR pathway. Moreover, ADT drug Bicalutamide and AZD3514 were less sensitive in the high CTHRC1 expression tumors.

CONCLUSIONS

As a potential molecular target of ADT resistance in prostate cancer, CTHRC1 provides a new promising molecular approach for the diagnosis and treatment of prostate cancer.

Integration of single-cell and bulk RNA-sequencing data reveals the prognostic potential of epithelial gene markers for prostate cancer

Prognostic transcriptomic signatures for prostate cancer (PCa) often overlook the cellular origin of expression changes, an important consideration given the heterogeneity of the disorder. Current clinicopathological factors inadequately predict biochemical recurrence, a critical indicator guiding post-treatment strategies following radical prostatectomy. To address this, we conducted a meta-analysis of four large-scale PCa datasets and found 33 previously reported PCa-associated genes to be consistently up-regulated in prostate tumours. By analysing single-cell RNA-sequencing data, we found these genes predominantly as markers in epithelial cells. Subsequently, we applied 97 advanced machine-learning algorithms across five PCa cohorts and developed an 11-gene epithelial expression signature. This signature robustly predicted biochemical recurrence-free survival (BCRFS) and stratified patients into distinct risk categories, with high-risk patients showing worse survival and altered immune cell populations. The signature outperformed traditional clinical parameters in larger cohorts and was overall superior to published PCa signatures for BCRFS. By analysing peripheral blood data, four of our signature genes showed potential as biomarkers for radiation response in patients with localised cancer and effectively stratified castration-resistant patients for overall survival. In conclusion, this study developed a novel epithelial gene-expression signature that enhanced BCRFS prediction and enabled effective risk stratification compared to existing clinical- and gene-expression-derived prognostic tools. Furthermore, a set of genes from the signature demonstrated potential utility in peripheral blood, a tissue amenable to minimally invasive sampling in a primary care setting, offering significant prognostic value for PCa patients without requiring a tumour biopsy.

Machine learning-based identification of co-expressed genes in prostate cancer and CRPC and construction of prognostic models

The objective of this study was to employ machine learning to identify shared differentially expressed genes (DEGs) in prostate cancer (PCa) initiation and castration resistance, aiming to establish a robust prognostic model and enhance understanding of patient prognosis for personalized treatment strategies. mRNA transcriptome data associated with Castration-Resistant Prostate Cancer (CRPC) were obtained from the GEO database. Differential expression analysis was conducted using the limma R package to compare normal prostate samples with PCa samples, and PCa samples with CRPC samples. Next, we applied LASSO regression, univariate, and multivariate COX regression analyses to pinpoint genes linked to prognosis and build prognostic models. Validation was performed using the TCGA_PRAD dataset to confirm expression differences of hub genes and explore their correlation with clinical variables and prognostic significance. We successfully established a prostate cancer risk prognostic model containing seven genes (KIF4A, UBE2C, FAM72D, CCDC78, HOXD9, LIX1 and SLC5A8) and verified its accuracy on an independent data set. The results of calibration curve and decision curve show that the model has potential clinical application value. The nomogram can accurately predict the prognosis of patients. Additionally, elevated expression of KIF4A, UBE2C, and FAM72D, or reduced expression of LIX1, correlated with advanced pathological T and N stages, clinical T stage, prostate-specific antigen (PSA) level, age at diagnosis, Gleason score, and shorter progression-free interval (PFI) (P < 0.05). By integrating bioinformatics analysis and clinical data, we not only established a reliable prognostic model for prostate cancer but also identified key genes pivotal in disease progression and treatment resistance. These findings provide novel insights and methodologies for assessing prognosis and tailoring treatment strategies for prostate cancer patients.

Associations of Prostate Tumor Immune Landscape with Vigorous Physical Activity and Prostate Cancer Progression

BACKGROUND

Vigorous physical activity has been associated with lower risk of fatal prostate cancer. However, mechanisms contributing to this relationship are not understood.

METHODS

We studied 117 men with prostate cancer in the University of North Carolina Cancer Survivorship Cohort (UNC CSC) who underwent radical prostatectomy and 101 radiation-treated patients with prostate cancer in FASTMAN. Structured questionnaires administered in UNC CSC assessed physical activity. In both studies, digital image analysis of hematoxylin and eosin-stained tissues was applied to quantify tumor-infiltrating lymphocytes in segmented regions. NanoString gene expression profiling in UNC CSC and microarray in FASTMAN were performed on tumor tissue, and a 50-gene signature utilized to predict immune cell types.

RESULTS

Vigorous recreational activity, reported by 34 (29.1%) UNC CSC men, was inversely associated with tumor-infiltrating lymphocyte abundance. Tumors of men reporting any vigorous activity versus none showed lower gene expression-predicted abundance of Th, exhausted CD4 T cells, and macrophages. T-cell subsets, including regulatory T cells, Th, Tfh, exhausted CD4 T cells, and macrophages, were associated with an increased risk of biochemical recurrence, only among men with ERG-positive tumors.

CONCLUSIONS

Vigorous activity was associated with lower prostate tumor inflammation and immune microenvironment differences. Macrophages and T-cell subsets, including those with immunosuppressive roles and those with lower abundance in men reporting vigorous exercise, were associated with worse outcomes in ERG-positive prostate cancer.

IMPACT

Our novel findings contribute to our understanding of the role of the tumor immune microenvironment in prostate cancer progression and may provide insights into how vigorous exercise could affect prostate tumor biology.

Molecular Hallmarks of Prostate-specific Membrane Antigen in Treatment-naïve Prostate Cancer

BACKGROUND AND OBJECTIVE

We characterized tumor prostate-specific membrane antigen (PSMA) levels as a reflection of cancer biology and treatment sensitivities for treatment-naïve prostate cancer.

METHODS

We first correlated PSMA positron emission tomography (PET) maximum standardized uptake values (SUVmax) in primary prostate cancer with tumor FOLH1 (PSMA RNA abundance) to establish RNA as a proxy (n = 55). We then discovered and validated molecular pathways associated with PSMA RNA levels in two large primary tumor cohorts. We validated those associations in independent cohorts (18 total; 5684 tumor samples) to characterize the pathways and treatment responses associated with PSMA.

KEY FINDINGS AND LIMITATIONS

PSMA RNA abundance correlates moderately with SUVmax (ρ = 0.41). In independent cohorts, androgen receptor signaling is more active in tumors with high PSMA. Accordingly, patients with high PSMA tumors experienced longer cancer-specific survival when managed with androgen deprivation therapy for biochemical recurrence (adjusted hazard ratio [AHR] 0.54 [0.34-0.87]; n = 174). PSMA low tumors possess molecular markers of resistance to radiotherapy. Consistent with this, patients with high PSMA tumors experience longer time to recurrence following primary radiotherapy (AHR 0.50 [0.28-0.90]; n = 248). In the SAKK09/10 trial (n = 224), patients with high PSMA tumors who were managed with salvage radiotherapy experienced longer time to progression in the 64-Gy arm (restricted mean survival time [RMST] +7.60 [0.05-15.16]), but this effect was mitigated in the 70-Gy arm (RMST 3.52 [-3.30 to 10.33]). Limitations include using PSMA RNA as a surrogate for PET SUVmax.

CONCLUSIONS AND CLINICAL IMPLICATIONS

PSMA levels in treatment-naïve prostate cancer differentiate tumor biology and treatment susceptibilities. These results warrant validation using PET metrics to substantiate management decisions based on imaging.

ClassifieR 2.0: expanding interactive gene expression-based stratification to prostate and high-grade serous ovarian cancer

BACKGROUND

Advances in transcriptional profiling methods have enabled the discovery of molecular subtypes within and across traditional tissue-based cancer classifications. Such molecular subgroups hold potential for improving patient outcomes by guiding treatment decisions and revealing physiological distinctions and targetable pathways. Computational methods for stratifying transcriptomic data into molecular subgroups are increasingly abundant. However, assigning samples to these subtypes and other transcriptionally inferred predictions is time-consuming and requires significant bioinformatics expertise. To address this need, we recently reported "ClassifieR," a flexible, interactive cloud application for the functional annotation of colorectal and breast cancer transcriptomes. Here, we report "ClassifieR 2.0" which introduces additional modules for the molecular subtyping of prostate and high-grade serous ovarian cancer (HGSOC).

RESULTS

ClassifieR 2.0 introduces ClassifieRp and ClassifieRov, two specialised modules specifically designed to address the challenges of prostate and HGSOC molecular classification. ClassifieRp includes sigInfer, a method we developed to infer commercial prognostic prostate gene expression signatures from publicly available gene-lists or indeed any user-uploaded gene-list. ClassifieRov utilizes consensus molecular subtyping methods for HGSOC, including tools like consensusOV, for accurate ovarian cancer stratification. Both modules include functionalities present in the original ClassifieR framework for estimating cellular composition, predicting transcription factor (TF) activity and single sample gene set enrichment analysis (ssGSEA).

CONCLUSIONS

ClassifieR 2.0 combines molecular subtyping of prostate cancer and HGSOC with commonly used sample annotation tools in a single, user-friendly platform, allowing scientists without bioinformatics training to explore prostate and HGSOC transcriptional data without the need for extensive bioinformatics knowledge or manual data handling to operate various packages. Our sigInfer method within ClassifieRp enables the inference of commercially available gene signatures for prostate cancer, while ClassifieRov incorporates consensus molecular subtyping for HGSOC. Overall, ClassifieR 2.0 aims to make molecular subtyping more accessible to the wider research community. This is crucial for increased understanding of the molecular heterogeneity of these cancers and developing personalised treatment strategies.

Unmasking Neuroendocrine Prostate Cancer with a Machine Learning-Driven Seven-Gene Stemness Signature That Predicts Progression

Prostate cancer (PCa) poses a significant global health challenge, particularly due to its progression into aggressive forms like neuroendocrine prostate cancer (NEPC). This study developed and validated a stemness-associated gene signature using advanced machine learning techniques, including Random Forest and Lasso regression, applied to large-scale transcriptomic datasets. The resulting seven-gene signature (KMT5C, DPP4, TYMS, CDC25B, IRF5, MEN1, and DNMT3B) was validated across independent cohorts and patient-derived xenograft (PDX) models. This signature demonstrated strong prognostic value for progression-free, disease-free, relapse-free, metastasis-free, and overall survival. Importantly, the signature not only identified specific NEPC subtypes, such as large-cell neuroendocrine carcinoma, which is associated with very poor outcomes, but also predicted a poor prognosis for PCa cases that exhibit this molecular signature, even when they were not histopathologically classified as NEPC. This dual prognostic and classifier capability makes the seven-gene signature a robust tool for personalized medicine, providing a valuable resource for predicting disease progression and guiding treatment strategies in PCa management.

Polyploid cancer cells reveal signatures of chemotherapy resistance

Therapeutic resistance in cancer significantly contributes to mortality, with many patients eventually experiencing recurrence after initial treatment responses. Recent studies have identified therapy-resistant large polyploid cancer cells in patient tissues, particularly in late-stage prostate cancer, linking them to advanced disease and relapse. Here, we analyzed bone marrow aspirates from 44 advanced prostate cancer patients and found the presence of CTC-IGC was significantly associated with poorer progression-free survival. Single cell copy number profiling of CTC-IGC displayed clonal origins with typical CTCs, suggesting complete polyploidization. Induced polyploid cancer cells from PC3 and MDA-MB-231 cell lines treated with docetaxel or cisplatin were examined through single cell DNA sequencing, RNA sequencing, and protein immunofluorescence. Novel RNA and protein markers, including HOMER1, TNFRSF9, and LRP1, were identified as linked to chemotherapy resistance. These markers were also present in a subset of patient CTCs and associated with recurrence in public gene expression data. This study highlights the prognostic significance of large polyploid tumor cells, their role in chemotherapy resistance, and their expression of markers tied to cancer relapse, offering new potential avenues for therapeutic development.

Precision molecular insights for prostate cancer prognosis: tumor immune microenvironment and cell death analysis of senescence-related genes by machine learning and single-cell analysis

BACKGROUND

Prostate cancer (PCa) is a prevalent malignancy among men, primarily originating from the prostate epithelium. It ranks first in global cancer incidence and second in mortality rates, with a rising trend in China. PCa's subtle initial symptoms, such as urinary issues, necessitate diagnostic measures like digital rectal examination, prostate-specific antigen (PSA) testing, and tissue biopsy. Advanced PCa management typically involves a multifaceted approach encompassing surgery, radiation, chemotherapy, and hormonal therapy. The involvement of aging genes in PCa development and progression, particularly through the mTOR pathway, has garnered increasing attention.

METHODS

This study aimed to explore the association between aging genes and biochemical PCa recurrence and construct predictive models. Utilizing public gene expression datasets (GSE70768, GSE116918, and TCGA), we conducted extensive analyses, including Cox regression, functional enrichment, immune cell infiltration estimation, and drug sensitivity assessments. The constructed risk score model, based on aging-related genes (ARGs), demonstrated superior predictive capability for PCa prognosis compared to conventional clinical features. High-risk genes positively correlated with risk, while low-risk genes displayed a negative correlation.

RESULTS

An ARGs-based risk score model was developed and validated for predicting prognosis in prostate adenocarcinoma (PRAD) patients. LASSO regression analysis and cross-validation plots were employed to select ARGs with prognostic significance. The risk score outperformed traditional clinicopathological features in predicting PRAD prognosis, as evidenced by its high AUC (0.787). The model demonstrated good sensitivity and specificity, with AUC values of 0.67, 0.675, 0.696, and 0.696 at 1, 3, 5, and 8 years, respectively, in the GEO cohort. Similar AUC values were observed in the TCGA cohort at 1, 3, and 5 years (0.67, 0.659, 0.667, and 0.743). The model included 12 genes, with high-risk genes positively correlated with risk and low-risk genes negatively correlated.

CONCLUSIONS

This study presents a robust ARGs-based risk score model for predicting biochemical recurrence in PCa patients, highlighting the potential significance of aging genes in PCa prognosis and offering enhanced predictive accuracy compared to traditional clinical parameters. These findings open new avenues for research on PCa recurrence prediction and therapeutic strategies.

Unmasking Neuroendocrine Prostate Cancer with a Machine Learning-Driven 7-Gene Stemness Signature that Predicts Progression

Prostate cancer (PCa) poses a significant global health challenge, particularly due to its progression into aggressive forms like neuroendocrine prostate cancer (NEPC). This study developed and validated a stemness-associated gene signature using advanced machine learning techniques, including Random Forest and Lasso regression, applied to large-scale transcriptomic datasets. The resulting 7-gene signature (KMT5C, MEN1, TYMS, IRF5, DNMT3B, CDC25B and DPP4) was validated across independent cohorts and patient-derived xenograft (PDX) models. The signature demonstrated strong prognostic value for progression-free, disease-free, relapse-free, metastasis-free, and overall survival. Importantly, the signature not only identified specific NEPC subtypes, such as large-cell neuroendocrine carcinoma, which is associated with very poor outcomes, but also predicted a poor prognosis for PCa cases that exhibit this molecular signature, even when they were not histopathologically classified as NEPC. This dual prognostic and classifier capability makes the 7-gene signature a robust tool for personalized medicine, providing a valuable resource for predicting disease progression and guiding treatment strategies in PCa management.

Leveraging cell death patterns to predict metastasis in prostate adenocarcinoma and targeting PTGDS for tumor suppression

Metastasis is the major cause of treatment failure in patients with prostate adenocarcinoma (PRAD). Diverse programmed cell death (PCD) patterns play an important role in tumor metastasis and hold promise as predictive indicators for PRAD metastasis. Using the LASSO Cox regression method, we developed PCD score (PCDS) based on differentially expressed genes (DEGs) associated with PCD. Clinical correlation, external validation, functional enrichment analysis, mutation landscape analysis, tumor immune environment analysis, and immunotherapy analysis were conducted. The role of Prostaglandin D2 Synthase (PTGDS) in PRAD was examined through in vitro experiments, single-cell, and Mendelian randomization (MR) analysis. PCDS is elevated in patients with higher Gleason scores, higher T stage, biochemical recurrence (BCR), and higher prostate-specific antigen (PSA) levels. Individuals with higher PCDS are prone to metastasis, metastasis after BCR, BCR, and castration resistance. Moreover, PRAD patients with low PCDS responded positively to immunotherapy. Random forest analysis and Mendelian randomization analysis identified PTGDS as the top gene associated with PRAD metastasis and in vitro experiments revealed that PTGDS was considerably downregulated in PRAD cells against normal prostate cells. Furthermore, the overexpression of PTGDS was found to suppress the migration, invasion, proliferationof DU145 and LNCaP cells. To sum up, PCDS may be a useful biomarker for forecasting the possibility of metastasis, recurrence, castration resistance, and the efficacy of immunotherapy in PRAD patients. Additionally, PTGDS was identified as a viable therapeutic target for the management of PRAD.

Machine learning based androgen receptor regulatory gene-related random forest survival model for precise treatment decision in prostate cancer

Background

It has been demonstrated that aberrant androgen receptor (AR) signaling contributes to the pathogenesis of prostate cancer (PCa). To date, the most efficacious strategy for the treatment of PCa remains to target the AR signaling axis. However, numerous PCa patients still face the issue of overtreatment or undertreatment. The establishment of a precise risk prediction model is urgently needed to distinguish patients with high-risk and select appropriate treatment modalities.

Methods

In this study, a consensus AR regulatory gene-related signature (ARS) was developed by integrating a total of 101 algorithm combinations of 10 machine learning algorithms. We evaluated the value of ARS in predicting patient prognosis and the therapeutic effects of the various treatments. Additionally, we conducted a screening of therapeutic targets and agents for high-risk patients, followed by the verification in vitro and in vivo.

Results

ARS was an independent risk factor for biochemical recurrence and distant metastasis in PCa patients. The enhanced and consistent prognostic predictive capability of ARS across various platforms was confirmed when compared with 44 previously published signatures. More importantly, PCa patients in the ARShigh group benefit more from PARP inhibitors and immunotherapy, while chemotherapy, radiotherapy, and AR-targeted therapy are more effective for ARSlow patients. The results of in silico screening suggest that AURKB could potentially serve as a promising therapeutic target for ARShigh patients.

Conclusions

Collectively, this prediction model based on AR regulatory genes holds great clinical translational potential to solve the dilemma of treatment choice and identify potential novel therapeutic targets in PCa.

Identification of circadian clock-related immunological prognostic index and molecular subtypes in prostate cancer

BACKGROUND

Evidence suggests that the circadian clock (CIC) is among the important factors for tumorigenesis. We aimed to provide new insights into CIC-mediated molecular subtypes and gene prognostic indexes for prostate cancer (PCa) patients undergoing radical prostatectomy (RP) or radical radiotherapy (RT).

METHODS

PCa data from TCGA was analyzed to identify differentially expressed genes (DEGs) with significant fold changes and p-values. A prognostic index called CIC-related gene prognostic index (CICGPI) was developed through clustering methods and survival analysis and validated on multiple data sets. The diagnostic accuracy of CICGPI for resistance to chemotherapy and radiotherapy was confirmed. Additionally, the interaction between tumor immune environment and CICGPI score was explored, along with their correlation with prognosis.

RESULTS

TOP2A, APOE, and ALDH2 were used to classify the PCa patients into two subtypes. Cluster 2 had a higher risk of biochemical recurrence (BCR) than cluster 1 for PCa patients undergoing RP or RT. A CIC-related gene prognostic index (CICGPI) was constructed using the above three genes for PCa patents in the TCGA database. The CICGPI score showed good prognostic value in the TCGA database and was externally confirmed by PCa patients in GSE116918, MSKCC2010 and GSE46602. In addition, the CICGPI score had a certain and high diagnostic accuracy for tumor chemoresistance (AUC: 0.781) and radioresistance (AUC: 0.988). For gene set variation analysis, we observed that both beta alanine metabolism and limonene and pinene degradation were upregulated in cluster 1 for PCa patients undergoing RP or RT. For PCa patients undergoing RP, cell cycle, homologous recombination, mismatch repair, and DNA replication were upregulated in cluster 2. A strongly positive relationship between cancer-related fibroblasts and CICGPI score was observed in PCa patients undergoing RP or RT. Moreover, a high density of CAFs was highly closely associated with poorer BCR-free survival of PCa patients.

CONCLUSIONS

In this study, we established CIC-related immunological prognostic index and molecular subtypes, which might be useful for the clinical practice.

Polyploid cancer cells reveal signatures of chemotherapy resistance

Therapeutic resistance in cancer significantly contributes to mortality, with many patients eventually experiencing recurrence after initial treatment responses. Recent studies have identified therapy-resistant large polyploid cancer cells in patient tissues, particularly in late-stage prostate cancer, linking them to advanced disease and relapse. Here, we analyzed bone marrow aspirates from 44 advanced prostate cancer patients and found the presence of circulating tumor cells with increased genomic content (CTC-IGC) was significantly associated with poorer progression-free survival. Single cell copy number profiling of CTC-IGC displayed clonal origins with typical CTCs, suggesting complete polyploidization. Induced polyploid cancer cells from PC3 and MDA-MB-231 cell lines treated with docetaxel or cisplatin were examined through single cell DNA sequencing, RNA sequencing, and protein immunofluorescence. Novel RNA and protein markers, including HOMER1, TNFRSF9, and LRP1, were identified as linked to chemotherapy resistance. These markers were also present in a subset of patient CTCs and associated with recurrence in public gene expression data. This study highlights the prognostic significance of large polyploid tumor cells, their role in chemotherapy resistance, and their expression of markers tied to cancer relapse, offering new potential avenues for therapeutic development.

Molecular analysis of archival diagnostic prostate cancer biopsies identifies genomic similarities in cases with progression post-radiotherapy, and those with de novo metastatic disease

BACKGROUND

It is important to identify molecular features that improve prostate cancer (PCa) risk stratification before radical treatment with curative intent. Molecular analysis of historical diagnostic formalin-fixed paraffin-embedded (FFPE) prostate biopsies from cohorts with post-radiotherapy (RT) long-term clinical follow-up has been limited. Utilizing parallel sequencing modalities, we performed a proof-of-principle sequencing analysis of historical diagnostic FFPE prostate biopsies. We compared patients with (i) stable PCa (sPCa) postprimary or salvage RT, (ii) progressing PCa (pPCa) post-RT, and (iii) de novo metastatic PCa (mPCa).

METHODS

A cohort of 19 patients with diagnostic prostate biopsies (n = 6 sPCa, n = 5 pPCa, n = 8 mPCa) and mean 4 years 10 months follow-up (diagnosed 2009-2016) underwent nucleic acid extraction from demarcated malignancy. Samples underwent 3'RNA sequencing (3'RNAseq) (n = 19), nanoString analysis (n = 12), and Illumina 850k methylation (n = 8) sequencing. Bioinformatic analysis was performed to coherently identify differentially expressed genes and methylated genomic regions (MGRs).

RESULTS

Eighteen of 19 samples provided useable 3'RNAseq data. Principal component analysis (PCA) demonstrated similar expression profiles between pPCa and mPCa cases, versus sPCa. Coherently differentially methylated probes between these groups identified ~600 differentially MGRs. The top 50 genes with increased expression in pPCa patients were associated with reduced progression-free survival post-RT (p < 0.0001) in an external cohort.

CONCLUSIONS

3'RNAseq, nanoString and 850k-methylation analyses are each achievable from historical FFPE diagnostic pretreatment prostate biopsies, unlocking the potential to utilize large cohorts of historic clinical samples. Profiling similarities between individuals with pPCa and mPCa suggests biological similarities and historical radiological staging limitations, which warrant further investigation.

Characteristics of the immune environment in prostate cancer as an adjunct to immunotherapy

Background and Aims

The tumor microenvironment (TME) exerts an important role in carcinogenesis and progression. Several investigations have suggested that immune cell infiltration (ICI) is of high prognostic importance for tumor progression and patient survival in many tumors, particularly prostate cancer. The pattern of immune infiltration of PCa, on the other hand, has not been thoroughly understood.

Methods

The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) datasets on PCa were obtained, and several datasets were merged into one data set using the "ComBat" algorithm. The ICI profiles of PCa patients were then to be uncovered by two computer techniques. The unsupervised clustering method was utilized to identify three ICI patterns in tumor samples, and Principal Component Analysis (PCA) was conducted to estimate the ICI score.

Results

Three different clusters of three ICIs were identified in 1341 PCa samples, which also correlated with different clinical features/characteristics and biological pathways. Patients with PCa are classified into high and low subtypes based on the ICI scores extracted from immune-associated signature genes. High ICI score subtypes are associated with a worse prognosis, which may intrigue the activation of cancer-related and immune-related pathways such as pathways involving Toll-like receptors, T-cell receptors, JAK-STAT, and natural killer cells. The ICI score was linked to tumor mutation load and immune/cancer-relevant signaling pathways, which explain prostate cancer's poor prognosis.

Conclusion

The findings of this study not only advanced our knowledge of the mechanism of immune response in the prostate tumor microenvironment but also provided a novel biomarker, that is, the ICI score, for disease prognosis and guiding precision immunotherapy.

Predicting prostate cancer recurrence: Introducing PCRPS, an advanced online web server

Background

Prostate cancer (PCa) is one of the leading causes of cancer death in men. About 30% of PCa will develop a biochemical recurrence (BCR) following initial treatment, which significantly contributes to prostate cancer-related deaths. In clinical practice, accurate prediction of PCa recurrence is crucial for making informed treatment decisions. However, the development of reliable models and biomarkers for predicting PCa recurrence remains a challenge. In this study, the aim is to establish an effective and reliable tool for predicting the recurrence of PCa.

Methods

We systematically screened and analyzed potential datasets to predict PCa recurrence. Through quality control analysis, low-quality datasets were removed. Using meta-analysis, differential expression analysis, and feature selection, we identified key genes associated with recurrence. We also evaluated 22 previously published signatures for PCa recurrence prediction. To assess prediction performance, we employed nine machine learning algorithms. We compared the predictive capabilities of models constructed using clinical variables, expression data, and their combinations. Subsequently, we implemented these machine learning models into a user-friendly web server freely accessible to all researchers.

Results

Based on transcriptomic data derived from eight multicenter studies consisting of 733 PCa patients, we screened 23 highly influential genes for predicting prostate cancer recurrence. These genes were used to construct the Prostate Cancer Recurrence Prediction Signature (PCRPS). By comparing with 22 published signatures and four important clinicopathological features, the PCRPS exhibited a robust and significantly improved predictive capability. Among the tested algorithms, Random Forest demonstrated the highest AUC value of 0.72 in predicting PCa recurrence in the testing dataset. To facilitate access and usage of these machine learning models by all researchers and clinicians, we also developed an online web server (https://urology1926.shinyapps.io/PCRPS/) where the PCRPS model can be freely utilized. The tool can also be used to (1) predict the PCa recurrence by clinical information or expression data with high accuracy. (2) provide the possibility of PCa recurrence by nine machine learning algorithms. Furthermore, using the PCRPS scores, we predicted the sensitivity of 22 drugs from GDSC2 and 95 drugs from CTRP2 to the samples. These predictions provide valuable insights into potential drug sensitivities related to the PCRPS score groups.

Conclusion

Overall, our study provides an attractive tool to further guide the clinical management and individualized treatment for PCa.

Transcriptome analysis of polyploid giant cancer cells and their progeny reveals a functional role for p21 in polyploidization and depolyploidization

Polyploid giant cancer cells (PGCC) are frequently detected in tumors and are increasingly recognized for their roles in chromosomal instability and associated genome evolution that leads to cancer recurrence. We previously reported that therapy stress promotes polyploidy, and that acid ceramidase plays a role in depolyploidization. In this study, we used an RNA-seq approach to gain a better understanding of the underlying transcriptomic changes that occur as cancer cells progress through polyploidization and depolyploidization. Our results revealed gene signatures that are associated with disease-free and/or overall survival in several cancers and identified the cell cycle inhibitor CDKN1A/p21 as the major hub in PGCC and early progeny. Increased expression of p21 in PGCC was limited to the cytoplasm. We previously demonstrated that the sphingolipid enzyme acid ceramidase is dispensable for polyploidization upon therapy stress but plays a crucial role in depolyploidization. The current study demonstrates that treatment of cells with ceramide is not sufficient for p53-independent induction of p21 and that knockdown of acid ceramidase, which hydrolyzes ceramide, does not interfere with upregulation of p21. In contrast, blocking the expression of p21 with UC2288 prevented the induction of acid ceramidase and inhibited both the formation of PGCC from parental cells as well as the generation of progeny from PGCC. Taken together, our data suggest that p21 functions upstream of acid ceramidase and plays an important role in polyploidization and depolyploidization.

Identification of cell differentiation trajectory-related gene signature to reveal the prognostic significance and immune landscape in prostate cancer based on multiomics analysis

Background

In the context of prostate cancer (PCa), the occurrence of biochemical recurrence (BCR) stands out as a pivotal factor significantly impacting prognosis, potentially leading to metastasis and mortality. However, the early detection of BCR poses a substantial challenge for PCa patients. There is an urgent need to pinpoint hub genes that can serve as predictive indicators for BCR in PCa patients.

Methods

Our primary goal was to identify cell differentiation trajectory-related gene signature in PCa patients by pseudo-time trajectory analysis. We further explored the functional enrichment of overlapped marker genes and probed clinically relevant modules and BCR-related genes using Weighted Gene Co-expression Network Analysis (WGCNA) in PCa patients. Key genes predicting recurrence-free survival were meticulously identified through univariate and multivariate Cox regression analyses. Subsequently, these genes were utilized to construct a prognostic gene signature, the expression, predictive efficacy, putative functions, and immunological landscape of which were thoroughly validated. Additionally, we employed immunohistochemistry (IHC) and a western blotting assay to quantify the expression of PYCR1 in clinical samples.

Results

Our single-cell RNA (scRNA) sequencing analysis unveiled three subgroups characterized by distinct differentiation trajectories, and the marker genes associated with these groups were extracted from PCa patients. These marker genes successfully classified the PCa sample into two molecular subtypes, demonstrating a robust correlation with clinical characteristics and recurrence-free survival. Through WGCNA and Lasso analysis, we identified four hub genes (KLK3, CD38, FASN, and PYCR1) to construct a risk profile of prognostic genes linked to BCR. Notably, the high-risk patient group exhibited elevated levels of B cell naive, Macrophage M0, and Macrophage M2 infiltration, while the low-risk group displayed higher levels of T cells CD4 memory activated and monocyte infiltration. Furthermore, IHC and western blotting assays confirmed the heightened expression of PYCR1 in PCa tissues.

Conclusion

This study leveraged the differentiation trajectory and genetic variability of the microenvironment to uncover crucial prognostic genes associated with BCR in PCa patients. These findings present novel perspectives for tailoring treatment strategies for PCa patients on an individualized basis.

A reliable transcriptomic risk-score applicable to formalin-fixed paraffin-embedded biopsies improves outcome prediction in localized prostate cancer

BACKGROUND

Clinical manifestation of prostate cancer (PCa) is highly variable. Aggressive tumors require radical treatment while clinically non-significant ones may be suitable for active surveillance. We previously developed the prognostic ProstaTrend RNA signature based on transcriptome-wide microarray and RNA-sequencing (RNA-Seq) analyses, primarily of prostatectomy specimens. An RNA-Seq study of formalin-fixed paraffin-embedded (FFPE) tumor biopsies has now allowed us to use this test as a basis for the development of a novel test that is applicable to FFPE biopsies as a tool for early routine PCa diagnostics.

METHODS

All patients of the FFPE biopsy cohort were treated by radical prostatectomy and median follow-up for biochemical recurrence (BCR) was 9 years. Based on the transcriptome data of 176 FFPE biopsies, we filtered ProstaTrend for genes susceptible to FFPE-associated degradation via regression analysis. ProstaTrend was additionally restricted to genes with concordant prognostic effects in the RNA-Seq TCGA prostate adenocarcinoma (PRAD) cohort to ensure robust and broad applicability. The prognostic relevance of the refined Transcriptomic Risk Score (TRS) was analyzed by Kaplan-Meier curves and Cox-regression models in our FFPE-biopsy cohort and 9 other public datasets from PCa patients with BCR as primary endpoint. In addition, we developed a prostate single-cell atlas of 41 PCa patients from 5 publicly available studies to analyze gene expression of ProstaTrend genes in different cell compartments.

RESULTS

Validation of the TRS using the original ProstaTrend signature in the cohort of FFPE biopsies revealed a relevant impact of FFPE-associated degradation on gene expression and consequently no significant association with prognosis (Cox-regression, p-value > 0.05) in FFPE tissue. However, the TRS based on the new version of the ProstaTrend-ffpe signature, which included 204 genes (of originally 1396 genes), was significantly associated with BCR in the FFPE biopsy cohort (Cox-regression p-value < 0.001) and retained prognostic relevance when adjusted for Gleason Grade Groups. We confirmed a significant association with BCR in 9 independent cohorts including 1109 patients. Comparison of the prognostic performance of the TRS with 17 other prognostically relevant PCa panels revealed that ProstaTrend-ffpe was among the best-ranked panels. We generated a PCa cell atlas to associate ProstaTrend genes with cell lineages or cell types. Tumor-specific luminal cells have a significantly higher TRS than normal luminal cells in all analyzed datasets. In addition, TRS of epithelial and luminal cells was correlated with increased Gleason score in 3 studies.

CONCLUSIONS

We developed a prognostic gene-expression signature for PCa that can be applied to FFPE biopsies and may be suitable to support clinical decision-making.

Establishment of novel ferroptosis-related prognostic subtypes correlating with immune dysfunction in prostate cancer patients

Background

We aimed to identify two new prognostic subtypes and create a predictive index for prostate cancer (PCa) patients based on ferroptosis database.

Methods

The nonnegative matrix factorization approach was used to identify molecular subtypes. We investigate the differences between cluster 1 and cluster 2 in terms of clinical features, functional pathways, tumour stemness, tumour heterogeneity, gene mutation and tumour immune microenvironment score after identifying the two molecular subtypes. Colony formation assay and flow cytometry assay were performed.

Results

The stratification of two clusters was closely connected to BCR-free survival using the nonnegative matrix factorization method, which was validated in the other three datasets. Furthermore, multivariate Cox regression analysis revealed that this classification was an independent risk factor for patients with PCa. Ribosome, aminoacyl tRNA production, oxidative phosphorylation, and Parkinson's disease-related pathways were shown to be highly enriched in cluster 1. In comparison to cluster 2, patients in cluster 1 exhibited significantly reduced CD4+ T cells, CD8+ T cells, neutrophils, dendritic cells and tumor immune microenvironment scores. Only HHLA2 was more abundant in cluster 1. Moreover, we found that P4HB downregulation could significantly inhibit the colony formation ability and contributed to cell apoptosis of C4-2B and DU145 cell lines.

Conclusions

We discovered two new prognostic subtypes associated with immunological dysfunction in PCa patients based on ferroptosis-related genes and found that P4HB downregulation could significantly inhibit the colony formation ability and contributed to cell apoptosis of PCa cell lines.

Integrating single-cell and bulk RNA sequencing data unveils antigen presentation and process-related CAFS and establishes a predictive signature in prostate cancer

BACKGROUND

Cancer-associated fibroblasts (CAFs) are heterogeneous and can influence the progression of prostate cancer in multiple ways; however, their capacity to present and process antigens in PRAD has not been investigated. In this study, antigen presentation and process-related CAFs (APPCAFs) were identified using bioinformatics, and the clinical implications of APPCAF-related signatures in PRAD were investigated.

METHODS

SMART technology was used to sequence the transcriptome of primary CAFs isolated from patients undergoing different treatments. Differential expression gene (DEG) screening was conducted. A CD4 + T-cell early activation assay was used to assess the activation degree of CD4 + T cells. The datasets of PRAD were obtained from The Cancer Genome Atlas (TCGA) database and NCBI Gene Expression Omnibus (GEO), and the list of 431 antigen presentation and process-related genes was obtained from the InnateDB database. Subsequently, APP-related CAFs were identified by nonnegative matrix factorization (NMF) based on a single-cell seq (scRNA) matrix. GSVA functional enrichment analyses were performed to depict the biological functions. A risk signature based on APPCAF-related genes (APPCAFRS) was developed by least absolute shrinkage and selection operator (LASSO) regression analysis, and the independence of the risk score as a prognostic factor was evaluated by univariate and multivariate Cox regression analyses. Furthermore, a biochemical recurrence-free survival (BCRFS)-related nomogram was established, and immune-related characteristics were assessed using the ssGSEA function. The immune treatment response in PRAD was further analyzed by the Tumor Immune Dysfunction and Exclusion (TIDE) tool. The expression levels of hub genes in APPCAFRS were verified in cell models.

RESULTS

There were 134 upregulated and 147 downregulated genes, totaling 281 differentially expressed genes among the primary CAFs. The functions and pathways of 147 downregulated DEGs were significantly enriched in antigen processing and presentation processes, MHC class II protein complex and transport vesicle, MHC class II protein complex binding, and intestinal immune network for IgA production. Androgen withdrawal diminished the activation effect of CAFs on T cells. NMF clustering of CAFs was performed by APPRGs, and pseudotime analysis yielded the antigen presentation and process-related CAF subtype CTSK + MRC2 + CAF-C1. CTSK + MRC2 + CAF-C1 cells exhibited ligand‒receptor connections with epithelial cells and T cells. Additionally, we found a strong association between CTSK + MRC2 + CAF-C1 cells and inflammatory CAFs. Through differential gene expression analysis of the CTSK + MRC2 + CAF-C1 and NoneAPP-CAF-C2 subgroups, 55 significant DEGs were identified, namely, APPCAFRGs. Based on the expression profiles of APPCAFRGs, we divided the TCGA-PRAD cohort into two clusters using NMF consistent cluster analysis, with the genetic coefficient serving as the evaluation index. Four APPCAFRGs, THBS2, DPT, COL5A1, and MARCKS, were used to develop a prognostic signature capable of predicting BCR occurrence in PRAD patients. Subsequently, a nomogram with stability and accuracy in predicting BCR was constructed based on Gleason grade (p = n.s.), PSA (p < 0.001), T stage (p < 0.05), and risk score (p < 0.01). The analysis of immune infiltration showed a positive correlation between the abundance of resting memory CD4 + T cells, M1 macrophages, resting dendritic cells, and the risk score. In addition, the mRNA expression levels of THBS2, DPT, COL5A1, and MARCKS in the cell models were consistent with the results of the bioinformatics analysis.

CONCLUSIONS

APPCAFRS based on four potential APPCAFRGs was developed, and their interaction with the immune microenvironment may play a crucial role in the progression to castration resistance of PRAD. This novel approach provides valuable insights into the pathogenesis of PRAD and offers unexplored targets for future research.

Distinct mesenchymal cell states mediate prostate cancer progression

In the complex tumor microenvironment (TME), mesenchymal cells are key players, yet their specific roles in prostate cancer (PCa) progression remain to be fully deciphered. This study employs single-cell RNA sequencing to delineate molecular changes in tumor stroma that influence PCa progression and metastasis. Analyzing mesenchymal cells from four genetically engineered mouse models (GEMMs) and correlating these findings with human tumors, we identify eight stromal cell populations with distinct transcriptional identities consistent across both species. Notably, stromal signatures in advanced mouse disease reflect those in human bone metastases, highlighting periostin's role in invasion and differentiation. From these insights, we derive a gene signature that predicts metastatic progression in localized disease beyond traditional Gleason scores. Our results illuminate the critical influence of stromal dynamics on PCa progression, suggesting new prognostic tools and therapeutic targets.

Integrated analysis of single-cell and bulk transcriptomics develops a robust neuroendocrine cell-intrinsic signature to predict prostate cancer progression

Neuroendocrine prostate cancer (NEPC) typically implies severe lethality and limited treatment options. The precise identification of NEPC cells holds paramount significance for both research and clinical applications, yet valid NEPC biomarker remains to be defined. Methods: Leveraging 11 published NE-related gene sets, 11 single-cell RNA-sequencing (scRNA-seq) cohorts, 15 bulk transcriptomic cohorts, and 13 experimental models of prostate cancer (PCa), we employed multiple advanced algorithms to construct and validate a robust NEPC risk prediction model. Results: Through the compilation of a comprehensive scRNA-seq reference atlas (comprising a total of 210,879 single cells, including 66 tumor samples) from 9 multicenter datasets of PCa, we observed inconsistent and inefficient performance among the 11 published NE gene sets. Therefore, we developed an integrative analysis pipeline, identifying 762 high-quality NE markers. Subsequently, we derived the NE cell-intrinsic gene signature, and developed an R package named NEPAL, to predict NEPC risk scores. By applying to multiple independent validation datasets, NEPAL consistently and accurately assigned NE feature and delineated PCa progression. Intriguingly, NEPAL demonstrated predictive capabilities for prognosis and therapy responsiveness, as well as the identification of potential epigenetic drivers of NEPC. Conclusion: The present study furnishes a valuable tool for the identification of NEPC and the monitoring of PCa progression through transcriptomic profiles obtained from both bulk and single-cell sources.

A novel tumor purity and immune infiltration-related model for predicting distant metastasis-free survival in prostate cancer

BACKGROUND

umor cells, immune cells and stromal cells jointly modify tumor development and progression. We aim to explore the potential effects of tumor purity on the immune microenvironment, genetic landscape and prognosis in prostate cancer (PCa).

METHODS

Tumor purity of prostate cancer patients was extracted from The cancer genome atlas (TCGA). Immune cellular proportions were calculated by the CIBERSORT. To identify critical modules related to tumor purity, we used weighted gene co-expression network analysis (WGCNA). Using STRING and Cytoscape, protein-protein interaction (PPI) networks were constructed and analyzed. A Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, Disease Ontology (DO), and Gene Set Enrichment Analysis (GSEA) enrichment analysis of identified modules was conducted. To identify the expression of key genes at protein levels, we used the Human Protein Atlas (HPA) platform.

RESULTS

A model of tumor purity score (TPS) was constructed in the gene expression omnibus series (GSE) 116,918 cohort. TCGA cohort served as a validation set and was employed to validate the TPS. TPS model, as an independent prognostic factor of distant metastasis-free survival (DMFS) in PCa. Patients had higher tumor purity and better prognosis in the low-TPS group. Tumor purity was related to the infiltration of mast cells and macrophage cells positively, whereas related to the infiltration of dendritic cells, T cells and B cells negatively in PCa. The nomogram based on TPS, Age, Gleason score and T stage had a good predictive value and could evaluate the prognosis of PCa metastasis. GO and KEGG enrichment analyses showed that hub genes mainly participate in T cell activation and T-helper lymphocytes (TH) differentiation. Hub genes were mainly enriched in primary immunodeficiency disease, according to DO analysis. SLAMF8 was identified as the most critical gene by Cytoscape and HPA analysis.

CONCLUSIONS

Dynamic changes in the immune microenvironment associated with tumor purity could correlate with a poor DMFS of low-purity PCa. The TPS can predict the DMFS of PCa. In addition, prostate cancer metastases may be related to immunosuppression caused by a disorder of the immune microenvironment.

Integrative multi-omics analysis unveils stemness-associated molecular subtypes in prostate cancer and pan-cancer: prognostic and therapeutic significance

BACKGROUND

Prostate cancer (PCA) is the fifth leading cause of cancer-related deaths worldwide, with limited treatment options in the advanced stages. The immunosuppressive tumor microenvironment (TME) of PCA results in lower sensitivity to immunotherapy. Although molecular subtyping is expected to offer important clues for precision treatment of PCA, there is currently a shortage of dependable and effective molecular typing methods available for clinical practice. Therefore, we aim to propose a novel stemness-based classification approach to guide personalized clinical treatments, including immunotherapy.

METHODS

An integrative multi-omics analysis of PCA was performed to evaluate stemness-level heterogeneities. Unsupervised hierarchical clustering was used to classify PCAs based on stemness signature genes. To make stemness-based patient classification more clinically applicable, a stemness subtype predictor was jointly developed by using four PCA datasets and 76 machine learning algorithms.

RESULTS

We identified stemness signatures of PCA comprising 18 signaling pathways, by which we classified PCA samples into three stemness subtypes via unsupervised hierarchical clustering: low stemness (LS), medium stemness (MS), and high stemness (HS) subtypes. HS patients are sensitive to androgen deprivation therapy, taxanes, and immunotherapy and have the highest stemness, malignancy, tumor mutation load (TMB) levels, worst prognosis, and immunosuppression. LS patients are sensitive to platinum-based chemotherapy but resistant to immunotherapy and have the lowest stemness, malignancy, and TMB levels, best prognosis, and the highest immune infiltration. MS patients represent an intermediate status of stemness, malignancy, and TMB levels with a moderate prognosis. We further demonstrated that these three stemness subtypes are conserved across pan-tumor. Additionally, the 9-gene stemness subtype predictor we developed has a comparable capability to 18 signaling pathways to make tumor diagnosis and to predict tumor recurrence, metastasis, progression, prognosis, and efficacy of different treatments.

CONCLUSIONS

The three stemness subtypes we identified have the potential to be a powerful tool for clinical tumor molecular classification in PCA and pan-cancer, and to guide the selection of immunotherapy or other sensitive treatments for tumor patients.

Identification and validation of cancer-associated fibroblast-related subtypes and the prognosis model of biochemical recurrence in prostate cancer based on single-cell and bulk RNA sequencing

BACKGROUND

Cancer-associated fibroblasts (CAFs) are an essential component of the tumor immune microenvironment that are involved in extracellular matrix (ECM) remodeling. We aim to investigate the characteristics of CAFs in prostate cancer and develop a biochemical recurrence (BCR)-related CAF signature for predicting the prognosis of PCa patients.

METHODS

The bulk RNA-seq and relevant clinical information were obtained from the TCGA and GEO databases, respectively. The infiltration scores of CAFs in prostate cancer patients were calculated using the MCP counter and EPIC algorithms. The single-cell RNA sequencing (scRNA-seq) was downloaded from the GEO database. Subsequently, univariate Cox regression analysis was employed to identify prognostic genes associated with CAFs. We identified two subtypes (C1 and C2) of prostate cancer that were associated with CAFs via non-negative matrix factorization (NMF) clustering. In addition, the BCR-related CAF signatures were constructed using Lasso regression analysis. Finally, a nomogram model was established based on the risk score and clinical characteristics of the patients.

RESULTS

Initially, we found that patients with high CAF infiltration scores had shorter biochemical recurrence-free survival (BCRFS) times. Subsequently, CAFs in four pairs of tumors and paracancerous tissues were identified. We discovered 253 significantly differentially expressed genes, of which 13 had prognostic significance. Using NMF clustering, we divided PCa patients into C1 and C2 subgroups, with the C1 subgroup having a worse prognosis and substantially enriched cell cycle, homologous recombination, and mismatch repair pathways. Furthermore, a BCR-related CAFs signature was established. Multivariate COX regression analysis confirmed that the BCR-related CAFs signature was an independent prognostic factor for BCR in PCa. In addition, the nomogram was based on the clinical characteristics and risk scores of the patient and demonstrated high accuracy and reliability for predicting BCR. Lastly, our findings indicate that the risk score may be a useful tool for predicting PCa patients' sensitivity to immunotherapy and drug treatment.

CONCLUSION

NMF clustering based on CAF-related genes revealed distinct TME immune characteristics between groups. The BCR-related CAF signature accurately predicted prognosis and immunotherapy response in prostate cancer patients, offering a promising new approach to cancer treatment.

Identification of senescence-related lncRNA prognostic index correlating with prognosis and radiosensitivity in prostate cancer patients

BACKGROUND

An increasing number of studies are shown how crucial a role cellular senescence plays in tumor development. In this study, we developed a senescence-related lncRNA prognostic index (SRLPI) to forecast radiosensitivity and the probability of biochemical recurrence (BCR) in patients with prostate cancer (PCa).

METHODS

PCa cohorts in TCGA and GEO databases were downloaded. Senescence-and prognosis-related lncRNA with differential expression in tumor and normal samples were identified and used to establish the SRLPI score. Mutation landscape, function pathway, tumor stemness and heterogeneity and tumor immune microenvironment were also analyzed. We performed the analysis using R 3.6.3 and the appropriate tools.

RESULTS

A SRLPI score was constructed based on SNHG1 and MIAT in the TCGA cohort. Our classification of PCa patients into high- and low-risk groups was based on the median SRLPI score. When compared to the low-SRLPI group, the high-SRLPI group was more vulnerable to BCR (HR: 3.33). In terms of BCR-free survival and metastasis-free survival, the GSE116918 showed similar findings. Surprisingly, the SRLPI score demonstrated a high level of radiosensitivity for diagnosis (AUC: 0.98). Age, Gleason score, T stage, N stage, positive lymph nodes, and residual tumor were all significantly greater in patients with high SRLPI scores. Furthermore, this score was significantly related to markers of senescence. Protein secretion and androgen response were found to be substantially enriched in the low-SRLPI group, whereas E2F targets were found to be strongly enriched in the high-SRLPI group for pathway analysis. For the tumor microenvironment assessment, B cells, CD8+ T cells, immune score and TIDE score were positively related to SRLPI score while endothelial level was negatively associated with SRLPI score with statistical significance.

CONCLUSIONS

We developed a SRLPI score that was related to prognosis and radiosensitivity and might be helpful in clinical practice.

Integrating single-cell and bulk RNA sequencing to predict prognosis and immunotherapy response in prostate cancer

Prostate cancer (PCa) stands as a prominent contributor to morbidity and mortality among males on a global scale. Cancer-associated fibroblasts (CAFs) are considered to be closely connected to tumour growth, invasion, and metastasis. We explored the role and characteristics of CAFs in PCa through bioinformatics analysis and built a CAFs-based risk model to predict prognostic treatment and treatment response in PCa patients. First, we downloaded the scRNA-seq data for PCa from the GEO. We extracted bulk RNA-seq data for PCa from the TCGA and GEO and adopted "ComBat" to remove batch effects. Then, we created a Seurat object for the scRNA-seq data using the package "Seurat" in R and identified CAF clusters based on the CAF-related genes (CAFRGs). Based on CAFRGs, a prognostic model was constructed by univariate Cox, LASSO, and multivariate Cox analyses. And the model was validated internally and externally by Kaplan-Meier analysis, respectively. We further performed GO and KEGG analyses of DEGs between risk groups. Besides, we investigated differences in somatic mutations between different risk groups. We explored differences in the immune microenvironment landscape and ICG expression levels in the different groups. Finally, we predicted the response to immunotherapy and the sensitivity of antitumour drugs between the different groups. We screened 4 CAF clusters and identified 463 CAFRGs in PCa scRNA-seq. We constructed a model containing 10 prognostic CAFRGs by univariate Cox, LASSO, and multivariate Cox analysis. Somatic mutation analysis revealed that TTN and TP53 were significantly more mutated in the high-risk group. Finally, we screened 31 chemotherapeutic drugs and targeted therapeutic drugs for PCa. In conclusion, we identified four clusters based on CAFs and constructed a new CAFs-based prognostic signature that could predict PCa patient prognosis and response to immunotherapy and might suggest meaningful clinical options for the treatment of PCa.

Development of a novel, clinically relevant anoikis-related gene signature to forecast prognosis in patients with prostate cancer

Introduction: Anoikis is a specific form of programmed cell death and is related to prostate cancer (PC) metastasis. This study aimed to develop a reliable anoikis-related gene signature to accurately forecast PC prognosis. Methods: Based on anoikis-related genes and The Cancer Genome Atlas (TCGA) data, anoikis-related molecular subtypes were identified, and their differences in disease-free survival (DFS), stemness, clinical features, and immune infiltration patterns were compared. Differential expression analysis of the two subtypes and weighted gene co-expression network analysis (WGCNA) were employed to identify clinically relevant anoikis-related differentially expressed genes (DEGs) between subtypes, which were then selected to construct a prognostic signature. The clinical utility of the signature was verified using the validation datasets GSE116918 and GSE46602. A nomogram was established to predict patient survival. Finally, differentially enriched hallmark gene sets were revealed between the different risk groups. Results: Two anoikis-related molecular subtypes were identified, and cluster 1 had poor prognosis, higher stemness, advanced clinical features, and differential immune cell infiltration. Next, 13 clinically relevant anoikis-related DEGs were identified, and five of them (CKS2, CDC20, FMOD, CD38, and MSMB) were selected to build a prognostic signature. This gene signature had a high prognostic value. A nomogram that combined Gleason score, T stage, and risk score could accurately predict patient survival. Furthermore, gene sets closely related with DNA repair were differentially expressed in the different risk groups. Conclusion: A novel, clinically relevant five-anoikis-related gene signature was a powerful prognostic biomarker for PC.

STAT3/LKB1 controls metastatic prostate cancer by regulating mTORC1/CREB pathway

Prostate cancer (PCa) is a common and fatal type of cancer in men. Metastatic PCa (mPCa) is a major factor contributing to its lethality, although the mechanisms remain poorly understood. PTEN is one of the most frequently deleted genes in mPCa. Here we show a frequent genomic co-deletion of PTEN and STAT3 in liquid biopsies of patients with mPCa. Loss of Stat3 in a Pten-null mouse prostate model leads to a reduction of LKB1/pAMPK with simultaneous activation of mTOR/CREB, resulting in metastatic disease. However, constitutive activation of Stat3 led to high LKB1/pAMPK levels and suppressed mTORC1/CREB pathway, preventing mPCa development. Metformin, one of the most widely prescribed therapeutics against type 2 diabetes, inhibits mTORC1 in liver and requires LKB1 to mediate glucose homeostasis. We find that metformin treatment of STAT3/AR-expressing PCa xenografts resulted in significantly reduced tumor growth accompanied by diminished mTORC1/CREB, AR and PSA levels. PCa xenografts with deletion of STAT3/AR nearly completely abrogated mTORC1/CREB inhibition mediated by metformin. Moreover, metformin treatment of PCa patients with high Gleason grade and type 2 diabetes resulted in undetectable mTORC1 levels and upregulated STAT3 expression. Furthermore, PCa patients with high CREB expression have worse clinical outcomes and a significantly increased risk of PCa relapse and metastatic recurrence. In summary, we have shown that STAT3 controls mPCa via LKB1/pAMPK/mTORC1/CREB signaling, which we have identified as a promising novel downstream target for the treatment of lethal mPCa.

Membrane tension-mediated stiff and soft tumor subtypes closely associated with prognosis for prostate cancer patients

BACKGROUND

Prostate cancer (PCa) is usually considered as cold tumor. Malignancy is associated with cell mechanic changes that contribute to extensive cell deformation required for metastatic dissemination. Thus, we established stiff and soft tumor subtypes for PCa patients from perspective of membrane tension.

METHODS

Nonnegative matrix factorization algorithm was used to identify molecular subtypes. We completed analyses using software R 3.6.3 and its suitable packages.

RESULTS

We constructed stiff and soft tumor subtypes using eight membrane tension-related genes through lasso regression and nonnegative matrix factorization analyses. We found that patients in stiff subtype were more prone to biochemical recurrence than those in soft subtype (HR 16.18; p < 0.001), which was externally validated in other three cohorts. The top ten mutation genes between stiff and soft subtypes were DNAH, NYNRIN, PTCHD4, WNK1, ARFGEF1, HRAS, ARHGEF2, MYOM1, ITGB6 and CPS1. E2F targets, base excision repair and notch signaling pathway were highly enriched in stiff subtype. Stiff subtype had significantly higher TMB and T cells follicular helper levels than soft subtype, as well as CTLA4, CD276, CD47 and TNFRSF25.

CONCLUSIONS

From the perspective of cell membrane tension, we found that stiff and soft tumor subtypes were closely associated with BCR-free survival for PCa patients, which might be important for the future research in the field of PCa.

Radiogenomic profiling of prostate tumors prior to external beam radiotherapy converges on a transcriptomic signature of TGF-β activity driving tumor recurrence

Background

Patients with localized prostate cancer have historically been assigned to clinical risk groups based on local disease extent, serum prostate specific antigen (PSA), and tumor grade. Clinical risk grouping is used to determine the intensity of treatment with external beam radiotherapy (EBRT) and androgen deprivation therapy (ADT), yet a substantial proportion of patients with intermediate and high risk localized prostate cancer will develop biochemical recurrence (BCR) and require salvage therapy. Prospective identification of patients destined to experience BCR would allow treatment intensification or selection of alternative therapeutic strategies.

Methods

Twenty-nine individuals with intermediate or high risk prostate cancer were prospectively recruited to a clinical trial designed to profile the molecular and imaging features of prostate cancer in patients undergoing EBRT and ADT. Whole transcriptome cDNA microarray and whole exome sequencing were performed on pretreatment targeted biopsy of prostate tumors (n=60). All patients underwent pretreatment and 6-month post EBRT multiparametric MRI (mpMRI), and were followed with serial PSA to assess presence or absence of BCR. Genes differentially expressed in the tumor of patients with and without BCR were investigated using pathways analysis tools and were similarly explored in alternative datasets. Differential gene expression and predicted pathway activation were evaluated in relation to tumor response on mpMRI and tumor genomic profile. A novel TGF-β gene signature was developed in the discovery dataset and applied to a validation dataset.

Findings

Baseline MRI lesion volume and PTEN/TP53 status in prostate tumor biopsies correlated with the activation state of TGF-β signaling measured using pathway analysis. All three measures correlated with the risk of BCR after definitive RT. A prostate cancer-specific TGF-β signature discriminated between patients that experienced BCR vs. those that did not. The signature retained prognostic utility in an independent cohort.

Interpretation

TGF-β activity is a dominant feature of intermediate-to-unfavorable risk prostate tumors prone to biochemical failure after EBRT with ADT. TGF-β activity may serve as a prognostic biomarker independent of existing risk factors and clinical decision-making criteria.

Funding

This research was supported by the Prostate Cancer Foundation, the Department of Defense Congressionally Directed Medical Research Program, National Cancer Institute, and the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research.

ATP8B1: A prognostic prostate cancer biomarker identified via genetic analysis

BACKGROUND

Controlling the asymmetric distribution of phospholipids across biological membranes plays a pivotal role in the life cycle of cells; one of the most important contributors that maintain this lipid asymmetry are phospholipid-transporting adenosine triphosphatases (ATPases). Although sufficient information regarding their association with cancer exists, there is limited evidence linking the genetic variants of phospholipid-transporting ATPase family genes to prostate cancer in humans.

METHODS

In this study, we investigated the association of 222 haplotype-tagging single-nucleotide polymorphisms (SNPs) in eight phospholipid-transporting ATPase genes with cancer-specific survival (CSS) and overall survival (OS) of 630 patients treated with androgen-deprivation therapy (ADT) for prostate cancer.

RESULTS

After multivariate Cox regression analysis and multiple testing correction, we found that ATP8B1 rs7239484 was remarkably associated with CSS and OS after ADT. A pooled analysis of multiple independent gene-expression datasets demonstrated that ATP8B1 was under-expressed in tumor tissues and that a higher ATP8B1 expression was associated with a better patient prognosis. Moreover, we established highly invasive sublines using two human prostate cancer cell lines to mimic cancer progression traits in vitro. The expression of ATP8B1 was consistently downregulated in both highly invasive sublines.

CONCLUSION

Our study indicates that rs7239484 is a prognostic factor for patients treated with ADT and that ATP8B1 can potentially attenuate prostate cancer progression.

Using biological constraints to improve prediction in precision oncology

Many gene signatures have been developed by applying machine learning (ML) on omics profiles, however, their clinical utility is often hindered by limited interpretability and unstable performance. Here, we show the importance of embedding prior biological knowledge in the decision rules yielded by ML approaches to build robust classifiers. We tested this by applying different ML algorithms on gene expression data to predict three difficult cancer phenotypes: bladder cancer progression to muscle-invasive disease, response to neoadjuvant chemotherapy in triple-negative breast cancer, and prostate cancer metastatic progression. We developed two sets of classifiers: mechanistic, by restricting the training to features capturing specific biological mechanisms; and agnostic, in which the training did not use any a priori biological information. Mechanistic models had a similar or better testing performance than their agnostic counterparts, with enhanced interpretability. Our findings support the use of biological constraints to develop robust gene signatures with high translational potential.

Immune-related gene index predicts metastasis for prostate cancer patients undergoing radical radiotherapy

In this study, we established a novel immunologic gene prognostic index (IGPI) to predict metastasis and provided new insights into tumor immune microenvironment (TIME) for PCa patients receiving radical radiotherapy. GBP2 and IGF1 were independent factors associated with metastasis-free survival. IGPI score was calculated based on GBP2 and IGF1 and this score was an independent risk factor for PCa patients undergoing radical radiotherapy. Patients with higher IGPI score were at higher risk of metastasis and biochemical recurrence, which were externally validated in the TCGA database and other GEO datasets. IGPI score had demonstrated moderate diagnostic ability of radiation resistance (AUC: 0.889). This score increased with the augment of Gleason score and T stage, as well as biochemical recurrence. Using EPIC, ESTIMATE and immunophenoscore (IPS) algorithms, cancer associated fibroblasts (CAFs), macrophages, stromal score, and estimate score were significantly higher in patients with metastasis group compared to their counterpart. Besides, for CAFs, macrophages, stromal score, and estimate score, patients with higher scores were at higher risk of metastasis, and the HRs were 3.65, 4.01, 4.27, and 3.78, respectively. IGPI score was highly positively associated with stromal score (coefficient: 0.39), immune score (coefficient: 0.43), estimate score (coefficient: 0.45), CAFs (coefficient: 0.42) and macrophages (coefficient: 0.42), while showing the opposite relationship with tumor purity (coefficient: - 0.45). In conclusion, we found that IGPI based on GBP2 and IGF1 might serve as a biomarker predicting metastasis for PCa patients. Besides, the current data further highlight the importance of CAFs in the metastatic process of PCa.

Identification of cancer-associated fibroblasts subtypes in prostate cancer

Introduction

Cancer-associated fibroblasts (CAFs) are one of the most abundant cell types in tumor microenvironment. However, the phenotypic and functional heterogeneities among CAFs have not been sufficiently investigated in prostate cancer.

Methods

We obtained and analyzed the single-cell RNA-sequencing data from 26 hormone-sensitive prostate cancer samples and 8 castration-resistant prostate cancer samples, along with the analysis of bulk-sequencing datasets. Furthermore, we performed multicolor immunofluorescence staining to verify the findings from the data analysis.

Results

We identified two major CAFs subtypes with distinct molecular characteristics and biological functions in prostate cancer microenvironment, namely αSMA+ CAV1+ CAFs-C0 and FN1+ FAP+ CAFs-C1. Another single-cell RNA-sequencing dataset containing 7 bone metastatic prostate cancer samples demonstrated that osteoblasts in the bone metastatic lesions comprised two subtypes with molecular characteristics and biological functions similar to CAFs-C0 and CAFs-C1 in the primary tumor sites. In addition, we discovered a transcriptional factor regulatory network depending on CAFs-C1. CAFs-C1, but not CAFs-C0, was associated with castration resistance and poor prognosis. We also found that CAFs-C1 signature was involved in treatment resistance to immune checkpoint inhibitors.

Discussion

In summary, our results identified the presence of heterogeneous CAFs subtypes in prostate cancer microenvironment and the potential of specific CAFs subtype as therapeutic target for castration-resistant prostate cancer.

The cGAS-STING pathway-related gene signature can predict patient prognosis and immunotherapy responses in prostate adenocarcinoma

The cyclic GMP-AMP synthase-stimulator of the interferon genes (cGAS-STING) pathway is essential in inflammation-driven tumor occurrence and progression. However, the prognostic roles and immune functions of cGAS-STING pathway-related genes in patients with prostate adenocarcinoma (PRAD) remain unclear. cGAS-STING pathway-related genes were obtained from the gene set enrichment analysis (GSEA) website. Univariate Cox regression analysis was performed to screen the prognosis-related hub genes in the cancer genome atlas (TCGA) and GSE116918 datasets. Unsupervised clustering analysis was performed to identify different clusters. The least absolute shrinkage and selection operator and multivariate Cox regression analyses were applied to develop a prognostic risk model. The prognostic values and predictive performance of risk signature were assessed by the Kaplan-Meier curve and receiver operating characteristic curve. The IMvigor210 cohort was used to investigate the potential values of the risk score in immunotherapeutic responses. Two clusters were identified based on the expression matrix of 12 prognosis-related genes. Specifically, better overall survival was observed in cluster 2 than cluster 1 in both datasets. Inflammation-related pathway enrichment and immune cell infiltration levels were altered between 2 clusters. Moreover, 6 genes (CASP8, GRK6, IL3RA, PLCB1, TBKBP1, and TNFSF10) were identified to generate a cGAS-STING pathway-related signature (CPRS). Survival analysis showed that patients in the high-risk group showed a more dismal survival than those in the low-risk group in TCGA and GSE116918 datasets. Notably, the CPRS can differentiate responsive patients from non-responsive individuals treated with PD-L1 blockades in an independent cohort. In addition, higher CPRS was associated with a more favorable prognosis. The proposed risk model was developed based on 6 cGAS-STING pathway related-genes, which can be used as a promising predictor for patient survival and immunotherapeutic responses in PRAD, contributing to treatment strategy-related decision-making.

Plasma cell subtypes analyzed using artificial intelligence algorithm for predicting biochemical recurrence, immune escape potential, and immunotherapy response of prostate cancer

Background

Plasma cells as an important component of immune microenvironment plays a crucial role in immune escape and are closely related to immune therapy response. However, its role for prostate cancer is rarely understood. In this study, we intend to investigate the value of a new plasma cell molecular subtype for predicting the biochemical recurrence, immune escape and immunotherapy response in prostate cancer.

Methods

Gene expression and clinicopathological data were collected from 481 prostate cancer patients in the Cancer Genome Atlas. Then, the immune characteristics of the patients were analyzed based on plasma cell infiltration fractions. The unsupervised clustering based machine learning algorithm was used to identify the molecular subtypes of the plasma cell. And the characteristic genes of plasma cell subtypes were screened out by three types of machine learning models to establish an artificial neural network for predicting plasma cell subtypes. Finally, the prediction artificial neural network of plasma cell infiltration subtypes was validated in an independent cohort of 449 prostate cancer patients from the Gene Expression Omnibus.

Results

The plasma cell fraction in prostate cancer was significantly decreased in tumors with high T stage, high Gleason score and lymph node metastasis. In addition, low plasma cell fraction patients had a higher risk of biochemical recurrence. Based on the differential genes of plasma cells, plasma cell infiltration status of PCa patients were divided into two independent molecular subtypes(subtype 1 and subtype 2). Subtype 1 tends to be immunosuppressive plasma cells infiltrating to the PCa region, with a higher likelihood of biochemical recurrence, more active immune microenvironment, and stronger immune escape potential, leading to a poor response to immunotherapy. Subsequently, 10 characteristic genes of plasma cell subtype were screened out by three machine learning algorithms. Finally, an artificial neural network was constructed by those 10 genes to predict the plasma cell subtype of new patients. This artificial neural network was validated in an independent validation set, and the similar results were gained.

Conclusions

Plasma cell infiltration subtypes could provide a potent prognostic predictor for prostate cancer and be an option for potential responders to prostate cancer immunotherapy.

Prognostic value of PNN in prostate cancer and its correlation with therapeutic significance

Prostate cancer (PCa) is the most common malignancy. New biomarkers are in demand to facilitate the management. The role of the pinin protein (encoded by PNN gene) in PCa has not been thoroughly explored yet. Using The Cancer Genome Atlas (TCGA-PCa) dataset validated with Gene Expression Omnibus (GEO) and protein expression data retrieved from the Human Protein Atlas, the prognostic and diagnostic values of PNN were studied. Highly co-expressed genes with PNN (HCEG) were constructed for pathway enrichment analysis and drug prediction. A prognostic signature based on methylation status using HCEG was constructed. Gene set enrichment analysis (GSEA) and the TISIDB database were utilised to analyse the associations between PNN and tumour-infiltrating immune cells. The upregulated PNN expression in PCa at both transcription and protein levels suggests its potential as an independent prognostic factor of PCa. Analyses of the PNN's co-expression network indicated that PNN plays a role in RNA splicing and spliceosomes. The prognostic methylation signature demonstrated good performance for progression-free survival. Finally, our results showed that the PNN gene was involved in splicing-related pathways in PCa and identified as a potential biomarker for PCa.

Leucine zipper protein 2 serves as a prognostic biomarker for prostate cancer correlating with immune infiltration and epigenetic regulation

Background

We sought to determine whether leucine zipper protein 2 (LUZP2) could benefit men with prostate cancer (PCa) undergoing radical radiotherapy (RT) or prostatectomy (RP).

Methods

Analysis was done on differentiating expression, clinical prognosis, co-expressed genes, immune infiltration, and epigenetic changes. All of our analyses were done using the R software (version 3.6.3) and the appropriate packages.

Results

In terms of PCa, tumor samples expressed LUZP2 more than normal samples did. In the TCGA database and GSE116918, we found that LUZP2 was the only independent risk factor for PCa. The shared enriched pathways for patients undergoing RP or RT were cell-cell adhesion, regulation of filopodium assembly, and extracellular matrix containing collagen. With the exception of TNFRSF14, we discovered that LUZP2 was negatively correlated with 21 immune checkpoints in PCa patients receiving RT. We found a significant inverse relationship between LUZP2 expression and the tumor immune environment, which included B cells, CD4+ T cells, neutrophils, macrophages, dendritic cells, stromal score, immune score, and estimate score, in patients receiving RP or RT. Additionally, tumor purity was positively correlated with LUZP2. We found that the drug bortezomib may be susceptible to the LUZP2. DNA methylation was significantly associated with the mRNA expression of LUZP2 in PCa patients from the TCGA database, and LUZP2 methylation was positively correlated with immune cells. The proliferative activity of various PCa cells, which correlated to different stages of this disease, was also found to be significantly reduced by LUZP2 reduction, according to the results of our experimental work.

Conclusions

We proposed a relatively comprehensive understanding of the roles of LUZP2 on PCa from the fresh perspective of senescence.

Mitochondrial Aldehyde Dehydrogenase 2 Represents a Potential Biomarker of Biochemical Recurrence in Prostate Cancer Patients

BACKGROUND

We aimed to explore the role of mitochondrial aldehyde dehydrogenase 2 (ALDH2) in prostate cancer (PCa) patients and provide insights into the tumor immune microenvironment (TME) for those patients undergoing radical radiotherapy.

METHODS

We performed all analyses using R version 3.6.3 and its suitable packages. Cytoscape 3.8.2 was used to establish network of competing endogenous RNAs (ceRNAs).

RESULTS

Downregulation of ADLH2 was significantly associated with higher risk of BCR-free survival (HR: 0.40, 95%CI: 0.24-0.68, p = 0.001) and metastasis-free survival (HR: 0.21, 95%CI: 0.09-0.49, p = 0.002). Additionally, ALDH2 repression contributed to significantly shorter BCR-free survival in the TCGA database (HR: 0.55, 95%CI: 0.33-0.93, p = 0.027). For immune checkpoints, patients that expressed a higher level of CD96 had a higher risk of BCR than their counterparts (HR: 1.79, 95%CI: 1.06-3.03, p = 0.032), as well as NRP1 (HR: 2.18, 95%CI: 1.29-3.69, p = 0.005). In terms of the TME parameters, the spearman analysis showed that ALDH was positively associated with B cells (r: 0.13), CD8+ T cells (r: 0.19), neutrophils (r: 0.13), and macrophages (r: 0.17). Patients with higher score of neutrophils (HR: 1.75, 95%CI: 1.03-2.95, p = 0.038), immune score (HR: 1.92, 95%CI: 1.14-3.25, p = 0.017), stromal score (HR: 2.52, 95%CI: 1.49-4.26, p = 0.001), and estimate score (HR: 1.81, 95%CI: 1.07-3.06, p = 0.028) had higher risk of BCR than their counterparts. Our ceRNA network found that PART1 might regulate the expression of ALDH via has-miR-578 and has-miR-6833-3p. Besides, PHA-793887, PI-103, and piperlongumine had better correlations with ALDH2.

CONCLUSIONS

We found that ALDH2 might serve as a potential biomarker predicting biochemical recurrence for PCa patients.

Particle Swarm Optimization Artificial Intelligence technique for gene signature discovery in transcriptomic cohorts

•

EBPSO identifies unique, accurate, and succinct gene signatures.

•

Key genes within the signatures provide biological insights its associated functions.

•

A web-based micro-framework developed for ease of use and real-time visualizations.

•

A promising alternative to traditional single gene signature generation.

•

Downstream analysis will better translate these signatures towards clinical translation.

The development of gene signatures is key for delivering personalized medicine, despite only a few signatures being available for use in the clinic for cancer patients. Gene signature discovery tends to revolve around identifying a single signature. However, it has been shown that various highly predictive signatures can be produced from the same dataset. This study assumes that the presentation of top ranked signatures will allow greater efforts in the selection of gene signatures for validation on external datasets and for their clinical translation. Particle swarm optimization (PSO) is an evolutionary algorithm often used as a search strategy and largely represented as binary PSO (BPSO) in this domain. BPSO, however, fails to produce succinct feature sets for complex optimization problems, thus affecting its overall runtime and optimization performance. Enhanced BPSO (EBPSO) was developed to overcome these shortcomings. Thus, this study will validate unique candidate gene signatures for different underlying biology from EBPSO on transcriptomics cohorts. EBPSO was consistently seen to be as accurate as BPSO with substantially smaller feature signatures and significantly faster runtimes. 100% accuracy was achieved in all but two of the selected data sets. Using clinical transcriptomics cohorts, EBPSO has demonstrated the ability to identify accurate, succinct, and significantly prognostic signatures that are unique from one another. This has been proposed as a promising alternative to overcome the issues regarding traditional single gene signature generation. Interpretation of key genes within the signatures provided biological insights into the associated functions that were well correlated to their cancer type.

A cellular senescence-related gene prognostic index for biochemical recurrence and drug resistance in patients with prostate cancer

In this study, we aimed to establish a novel cellular senescence-related gene prognostic index (CSG PI) to predict biochemical recurrence (BCR) and drug resistance in patients with prostate cancer (PCa) undergoing radical radiotherapy or prostatectomy. We performed all analyses using R version 3.6.3 and its suitable packages. Cytoscape 3.8.2 was used to establish a network of transcription factors and competing endogenous RNAs. Three cellular senescence-related genes were used to establish the CSGPI. We observed that CSGPI was an independent risk factor for BCR in PCa patients (HR: 2.62; 95% CI: 1.55-4.44), consistent with the results of external validation (HR: 1.88; 95% CI: 1.12-3.14). The CSGPI had a moderate diagnostic effect on drug resistance (AUC: 0.812, 95% CI: 0.586-1.000). The lncRNA PART1 was significantly associated with BCR (HR: 0.46; 95% CI: 0.27-0.77), and might modulate the mRNA expression of definitive genes through interactions with 57 miRNAs. Gene set enrichment analysis indicated that CSGPI was closely related to ECM receptor interaction, focal adhesion, TGF beta signaling pathway, pathway in cancer, regulation of actin cytoskeleton, and so on. Immune checkpoint analysis showed that PDCD1LG2 and CD96 were significantly higher in the BCR group compared to non-BCR group, and patients with higher expression of CD96 were more prone to BCR than their counterparts (HR: 1.79; 95% CI: 1.06-3.03). In addition, the CSGPI score was significantly associated with the mRNA expression of HAVCR2, CD96, and CD47. Analysis of mismatch repair and methyltransferase genes showed that DNMT3B was more highly expressed in the BCR group and that patients with higher expression of DNMT3B experienced a higher risk of BCR (HR: 2.08; 95% CI: 1.23-3.52). We observed that M1 macrophage, CD8+ T cells, stromal score, immune score, and ESTIMATE score were higher in the BCR group. In contrast, tumor purity was less scored in the BCR group. Spearman analysis revealed a positive relationship between CSGPI and M1 macrophages, CD4+ T cells, dendritic cells, stromal score, immune score, and ESTIMATE score. In conclusion, we found that the CSGPI might serve as a biomarker to predict BCR and drug resistance in PCa patients. Moreover, CD96 and DNMT3B might be potential treatment targets, and immune evasion might contribute to the BCR process of PCa.

Utilization of hypoxia-derived gene signatures to predict clinical outcomes and immune checkpoint blockade therapy responses in prostate cancer

Background: Increasing evidences show a clinical significance in the interaction between hypoxia and prostate cancer. However, reliable prognostic signatures based on hypoxia have not been established yet.

Methods: We screened hypoxia-related gene modules by weighted gene co-expression network analysis (WGCNA) and established a hypoxia-related prognostic risk score (HPRS) model by univariate Cox and LASSO-Cox analyses. In addition, enriched pathways, genomic mutations, and tumor-infiltrating immune cells in HPRS subgroups were analyzed and compared. HPRS was also estimated to predict immune checkpoint blockade (ICB) therapy response.

Results: A hypoxia-related 22-gene prognostic model was established. Furthermore, three independent validation cohorts showed moderate performance in predicting biochemical recurrence-free (BCR-free) survival. HPRS could be a useful tool in selecting patients who can benefit from ICB therapy. The CIBERSORT results in our study demonstrated that hypoxia might act on multiple T cells, activated NK cells, and macrophages M1 in various ways, suggesting that hypoxia might exert its anti-tumor effects by suppressing T cells and NK cells.

Conclusion: Hypoxia plays an important role in the progression of prostate cancer. The hypoxia-derived signatures are promising biomarkers to predict biochemical recurrence-free survival and ICB therapy responses in patients with prostate cancer.

SLC26A4 correlates with homologous recombination deficiency and patient prognosis in prostate cancer

Background

Homologous recombination deficiency (HRD) is closely associated with patient prognosis and treatment options in prostate cancer (PCa). However, there is a lack of quantitative indicators related to HRD to predict the prognosis of PCa accurately.

Methods

We screened HRD-related genes based on the HRD scores and constructed an HRD cluster system to explore different clinicopathological, genomic, and immunogenomic patterns among the clusters. A risk signature, HRDscore, was established and evaluated by multivariate Cox regression analysis. We noticed that SLC26A4, a model gene, demonstrated unique potential to predict prognosis and HRD in PCa. Multi-omics analysis was conducted to explore its role in PCa, and the results were validated by qRT-PCR and immunohistochemistry.

Results

Three HRD clusters were identified with significant differences in patient prognosis, clinicopathological characteristics, biological pathways, immune infiltration characteristics, and regulation of immunomodulators. Further analyses revealed that the constructed HRDscore system was an independent prognostic factor of PCa patients with good stability. Finally, we identified a single gene, SLC26A4, which significantly correlated with prognosis in three independent cohorts. Importantly, SLC26A4 was confirmed to distinguish PCa (AUC for mRNA 0.845; AUC for immunohistochemistry score 0.769) and HRD (AUC for mRNA 0.911; AUC for immunohistochemistry score 0.689) at both RNA and protein levels in our cohort.

Conclusion

This study introduces HRDscore to quantify the HRD pattern of individual PCa patients. Meanwhile, SLC26A4 is a novel biomarker and can reasonably predict the prognosis and HRD in PCa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03513-5.

LTBP2 inhibits prostate cancer progression and metastasis via the PI3K/AKT signaling pathway

Biochemical recurrence (BCR) is a cause of concern in advanced prostate cancer (PCa). Thus, novel diagnostic biomarkers are required to improve clinical care. However, research on PCa immunotherapy is also scarce. Hence, the present study aimed to explore promising BCR-related diagnostic biomarkers, and their expression pattern, prognostic value, immune response effects, biological functions, and possible molecular mechanisms were evaluated. GEO datasets (GSE46602, GSE70768, and GSE116918) were downloaded and merged as the training cohort, and differential expression analysis was performed. Lasso regression and SVM-RFE algorithm, as well as PPI analysis and MCODE algorithm, were then applied to filter BCR-related biomarker genes. The CIBERSORT and estimation of stromal and immune cells in malignant tumor tissues using expression data (ESTIMATE) methods were used to calculate the fractions of tumor-infiltrating immune cells. GO/DO enrichment analyses were used to identify the biological functions. The expression of latent transforming growth factor β-binding protein 2 (LTBP2) was determined by RT-qPCR and western blotting. The role of LTBP2 in PCa was determined by CCK-8, Transwell, and the potential mechanism was investigated by KEGG and GSEA and confirmed by western blotting. In total, 44 BCR-related differentially expressed genes (DEGs) in the training cohort were screened. LTBP2 was found to be a diagnostic biomarker of BCR in PCa and was associated with CD4⁺ T-cell infiltration and response to anti-PD-1/PD-L1 immunotherapy. Subsequently, using the ESTIMATE algorithm, it was identified that LTBP2 was associated with the tumor microenvironment and could be a predictor of the clinical benefit of immune checkpoint blockade. Finally, the expression and biological function of LTBP2 were evaluated via cellular experiments. The results showed that LTBP2 was downregulated in PCa cells and inhibited PCa proliferation and metastasis via the PI3K/AKT signaling pathway in vitro. In conclusion, LTBP2 was a promising diagnostic biomarker of BCR of PCa and had an important role in CD4⁺ T-cell recruitment. Moreover, it was associated with immunotherapy in patients with PCa who developed BCR, and it inhibited PCa proliferation and metastasis via the PI3K/AKT signaling pathway in vitro.

A Novel Angiogenesis-Related Gene Signature to Predict Biochemical Recurrence of Patients with Prostate Cancer following Radical Therapy

Background

Prostate cancer (PCa) is one of the most common malignancies in males; we aim to establish a novel angiogenesis-related gene signature for biochemical recurrence (BCR) prediction in PCa patients following radical therapy.

Methods

Gene expression profiles and corresponding clinicopathological data were acquired from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. We quantified the levels of various cancer hallmarks and identified angiogenesis as the primary risk factor for BCR. Then machine learning methods combined with Cox regression analysis were used to screen prognostic genes and construct an angiogenesis-related gene signature, which was further verified in external cohorts. Furthermore, estimation of immune cell abundance and prediction of drug responses were also conducted to detect potential mechanisms.

Results

Angiogenesis was regarded as the leading risk factor for BCR in PCa patients (HR = 1.58, 95% CI: 1.38–1.81), and a novel prognostic signature based on three genes (NRP1, JAG2, and VCAN) was developed in the training cohort and successfully validated in another three independent cohorts. In all datasets, this signature was identified as a prognostic factor with promising and robust predictive values. Moreover, it also predicted higher infiltration of regulatory T cells and M2-polarized macrophages and increased drug sensitivity of angiogenesis inhibitors in high-risk patients.

Conclusions

The angiogenesis-related three-gene signature may serve as an independent prognostic factor for BCR, which would contribute to risk stratification and personalized management of PCa patients after radical therapy in clinical practice.