Identification of core genes associated with prostate cancer progression and outcome via bioinformatics analysis in multiple databases

An analysis of the prognostic role of reactive oxygen species-associated genes in breast cancer

BACKGROUND

This study aimed to type breast cancer in relation to reactive oxygen species (ROS), clinical indicators, single nucleotide variant (SNV) mutations, functional differences, immune infiltration, and predictive responses to immunotherapy or chemotherapy, and constructing a prognostic model.

METHODS

We used uniCox analysis, ConsensusClusterPlus, and the proportion of ambiguous clustering (PAC) to analyze The Cancer Genome Atlas (TCGA) data to determine optimal groupings and obtain differentially expressed ROS-related genes. Clinical indicators were then combined with the classification results and the Chi-square test was used to assess differences. We further examined SNV mutations, and functional differences using gene set enrichment analysis (GSEA) analysis, the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, immune cell infiltration, and response to immunotherapy and chemotherapy. A prognostic model for breast cancer was constructed using these differentially expressed genes, immunotherapy or chemotherapy responses, and survival curves. RT-qPCR was used to detect the differences in the expression of LCE3D, CA1, PIRT and SMR3A in breast cancer cell lines and normal breast epithelial cell line.

RESULTS

We identified two distinct tumor types with significant differences in ROS-related gene expression, clinical indicators, SNV mutations, functional pathways, and immune infiltration. The response to specific chemotherapy drugs and immunotherapy treatments also documented significant differences. The prognostic model constructed with 16 genes linked to survival could efficiently divide patients into high- and low-risk groups. The high-risk group showed a poorer prognosis, higher tumor purity, distinct immune microenvironment, and lower immunotherapy response. RT-qPCR results showed that LCE3D, CA1, PIRT and SMR3A are highly expressed in breast cancer.

CONCLUSION

Our methodical examination presented an enhanced insight into the molecular and immunological heterogeneity of breast cancer. It can contribute to the understanding of prognosis and offer valuable insights for personalized treatment strategies. Further, the prognostic model can potentially serve as a powerful tool for risk stratification and therapeutic decision-making in clinical settings.

N6-Methyladenosine reader HNRNPC-mediated downregulation of circITCH prevents miR-224-3p sequestering and contributes to tumorigenesis in nasopharyngeal carcinoma

BACKGROUND

N6-Methyladenosine (m6A) RNA methylation modulators are implicated in nasopharyngeal carcinoma (NPC). Circular RNAs (circRNAs) stimulate/inhibit the development of NPC by sponging microRNAs (miRNAs). Herein, m6A modifications affecting the circRNA/miRNA axis in NPC were explored.

METHODS

Twenty prognostic m6A RNA methylation regulators were identified from 504 head/neck squamous cell carcinoma and 44 normal samples from The Cancer Genome Atlas (TCGA). Differentially expressed miRNAs were screened from the TCGA and Gene Expression Omnibus (GEO) databases. RNA-binding protein (RBP)-circRNA and circRNA-miRNA interactive pairs were verified using RBPmap and RNAhybrid, respectively. The RBP/circRNA/miRNA network was constructed using Cytoscape. Furthermore, CircITCH (hsa_circ_00059948), HNRNPC, and miR-224-3p expressions were detected by western blotting and quantitative polymerase chain reaction. The role of circITCH in NPC was examined using apoptosis, scratch wound healing, transwell invasion, and cell counting kit-8 assays. Finally, CircITCH-miR-224-3p and circITCH-HNRNPC interactions were assessed by dual-luciferase reporter and RNA-immunoprecipitation (RIP) assays, respectively.

RESULTS

Bioinformatics analysis revealed that high pathological grade, late-stage tumors, and low survival were associated with increased HNRNPC expression. MiR-224-3p was upregulated in NPC and sequestered by circITCH. Construction of the RBP/circRNA/miRNA network highlighted the HNRNPC/circITCH/miR-224-3p axis. In vitro experiments demonstrated decreased circITCH expression and increased HNRNPC and miR-224-3p expressions in NPC. In NPC cells overexpressing circITCH, HNRNPC and miR-224-3p expressions were significantly decreased. Dual-luciferase assays demonstrated a targeting relationship between circITCH and miR-224-3p, and RIP assays demonstrated interaction of HNRNPC targets with circITCH.

CONCLUSION

CircITCH overexpression inhibited NPC progression by sequestering miR-224-3p, and HNRNPC reduced circITCH expression through direct interaction.

Marker Ki-67 is a potential biomarker for the diagnosis and prognosis of prostate cancer based on two cohorts

BACKGROUND

Prostate cancer (PCa) is a widespread malignancy, predominantly affecting elderly males, and current methods for diagnosis and treatment of this disease continue to fall short. The marker Ki-67 (MKI67) has been previously demonstrated to correlate with the proliferation and metastasis of various cancer cells, including those of PCa. Hence, verifying the association between MKI67 and the diagnosis and prognosis of PCa, using bioinformatics databases and clinical data analysis, carries significant clinical implications.

AIM

To explore the diagnostic and prognostic efficacy of antigens identified by MKI67 expression in PCa.

METHODS

For cohort 1, the efficacy of MKI67 diagnosis was evaluated using data from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases. For cohort 2, the diagnostic and prognostic power of MKI67 expression was further validated using data from 271 patients with clinical PCa.

RESULTS

In cohort 1, MKI67 expression was correlated with prostate-specific antigen (PSA), Gleason Score, T stage, and N stage. The receiver operating characteristic (ROC) curve showed a strong diagnostic ability, and the Kaplan-Meier method demonstrated that MKI67 expression was negatively associated with the progression-free interval (PFI). The time-ROC curve displayed a weak prognostic capability for MKI67 expression in PCa. In cohort 2, MKI67 expression was significantly related to the Gleason Score, T stage, and N stage; however, it was negatively associated with the PFI. The time-ROC curve revealed the stronger prognostic capability of MKI67 in patients with PCa. Multivariate COX regression analysis was performed to select risk factors, including PSA level, N stage, and MKI67 expression. A nomogram was established to predict the 3-year PFI.

CONCLUSION

MKI67 expression was positively associated with the Gleason Score, T stage, and N stage and showed a strong diagnostic and prognostic ability in PCa.

Metabolism-related gene vaccines and immune infiltration in ovarian cancer: A novel risk score model of machine learning

BACKGROUND

The present study aims to develop a metabolic gene signature to evaluate the survival rate of ovarian cancer (OC) patients and analyze the potential mechanisms of metabolic genes in OC because the difficulty in early detection of OC often leads to poor treatment outcomes.

METHODS

A non-negative matrix factorization algorithm was applied to determine molecular subtypes according to metabolism genes. To build a risk prognosis model, least absolute shrinkage and selection operator multivariate Cox analysis was carried out with weighted correlation network analysis (WCGNA). Glycolytic flux and mitochondrial function were evaluated by conducting seahorse analysis.

RESULTS

On the basis of metabolism-related genes, the two subtypes of OC samples present in The Cancer Genome Atlas database were distinguished. An analysis of WGCNA identified 1056 genes. Lastly, a 10-gene signature (CMAS, ADH1B, PLA2G2D, BHMT, CACNA1C, AADAC, ALOX12, CYP2R1, SCN1B and ME1) was constructed that demonstrated promising performance in predicting outcome in patients with OC. The RiskScore of the gene signature was linked to microenvironment cell infiltration and immune checkpoint. Higher RiskScores were associated with poorer results for OC patients. Seahorse analysis shows the influence of CMAS in cell energy metabolism.

CONCLUSIONS

In the present study, a novel marker for evaluating the survival of OC patients was developed through the creation of a gene signature incorporating metabolism-related genes. Our knowledge of immunotherapy and microenvironment cell infiltration may be enriched by evaluating metabolism-related gene modification patterns.

Integrative analysis of cuproptosis-associated genes for predicting immunotherapy response in single-cell and multi-cohort studies

BACKGROUND

The role of genes associated with the cuproptosis cell signaling pathway in prognosis and immunotherapy in ovarian cancer (OC) has been extensively investigated. In this study, we aimed to explore these mechanisms and establish a prognostic model for patients with OC using bioinformatics techniques.

METHODS

We obtained the single cell sequencing data of ovarian cancer from the Gene Expression Omnibus (GEO) database and preprocessed the data. We analyzed a variety of factors including cuproptosis cell signal score, transcription factors, tumorigenesis and progression signals, gene set variation analysis (GSVA) and intercellular communication. Differential gene analysis was performed between groups with high and low cuproptosis cell signal scores, as well as Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses. Using bulk RNA sequencing data from The Cancer Genome Atlas, we used the least absolute shrinkage and selection operator (LASSO)-Cox algorithm to develop cuproptosis cell signaling pathword-related gene signatures and validated them with GEO ovarian cancer datasets. In addition, we analyzed the inherent rules of the genes involved in building the model using a variety of bioinformatics methods, including immune-related analyses and single nucleotide polymorphisms. Molecular docking is used to screen potential therapeutic drugs. To confirm the analysis results, we performed various wet experiments such as western blot, cell counting kit 8 (CCK8) and clonogenesis tests to verify the role of the Von Willebrand Factor (VWF) gene in two ovarian cancer cell lines.

RESULTS

Based on single-cell data analysis, we found that endothelial cells and fibroblasts showed active substance synthesis and signaling pathway activation in OC, which further promoted immune cell suppression, cancer cell proliferation and metastasis. Ovarian cancer has a high tendency to metastasize, and cancer cells cooperate with other cells to promote disease progression. We developed a signature consisting of eight cuproptosis-related genes (CRGs) (MAGEF1, DNPH1, RARRES1, NBL1, IFI27, VWF, OLFML3 and IGFBP4) that predicted overall survival in patients with ovarian cancer. The validity of this model is verified in an external GEO validation set. We observed active infiltrating states of immune cells in both the high- and low-risk groups, although the specific cells, genes and pathways of activation differed. Gene mutation analysis revealed that TP53 is the most frequently mutated gene in ovarian cancer. We also predict small molecule drugs associated with CRGs and identify several potential candidates. VWF was identified as an oncogene in ovarian cancer, and the protein was expressed at significantly higher levels in tumor samples than in normal samples. The high-score model of the cuproptosis cell signaling pathway was associated with the sensitivity of OC patients to immunotherapy.

CONCLUSIONS

Our study provides greater insight into the mechanisms of action of genes associated with the cuproptosis cell signaling pathway in ovarian cancer, highlighting potential targets for future therapeutic interventions.

Prognostic role of long non-coding RNA USP30-AS1 in ovarian cancer: insights into immune cell infiltration in the tumor microenvironment

Ovarian cancer represents a formidable gynecologic malignancy bearing a dismal prognosis owing to the dearth of reliable early detection approaches and a high recurrence rate. Long non-coding RNAs (lncRNAs) have garnered immense attention as key orchestrators involved in diverse biological processes and take part in cancer initiation and progression. The present study investigated the potential significance of lncRNA USP30-AS1 in ovarian cancer prognosis, as well as its putative association with immune cell infiltration in tumor immune microenvironment (TIME). By analyzing publicly available datasets, we identified six lncRNAs with prognostic prediction ability, including USP30-AS1. The results revealed a significant positive correlation of USP30-AS1 expression with the infiltration of immune cells such as Th1 cells, TFH, CD8 T cells, B cells, antigen-presenting dendritic cells (aDC), and plasmacytoid dendritic cells (pDC) in ovarian cancer specimens. These findings provide compelling evidence of the potential involvement of lncRNA in the regulation of the TME in ovarian carcinoma. The outcomes from this study underscore the potential of USP30-AS1 as a promising prognostic biomarker for ovarian cancer. Additionally, the findings offer significant insights into the plausible role of lncRNAs in modulating immune activities, thus adding to our understanding of the disease biology. Additional investigations are necessary to unravel the molecular mechanisms underpinning these connections and validate the results seen in independent cohorts and experimental models.

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.

Construction of a cuproptosis-associated lncRNA prognostic signature for bladder cancer and experimental validation of cuproptosis-related lncRNA UBE2Q1-AS1

Introduction

Bladder cancer (BLCA) is the ninth most common malignancy worldwide and the fourth most common cancer in men. Copper levels are significantly altered in patients with thyroid, breast, lung, cervical, ovarian, pancreatic, oral, gastric, bladder, and prostate cancers. Outcomes can be predicted by constructing signatures using lncRNA-related genes associated with outcomes.

Methods

We identified lncRNAs related to outcomes, those differentially expressed in bladder cancer, and cuproptosis-related lncRNAs from TCGA. We identified the intersection to obtain 12 genes and established a prognostic risk signature consisting of eight genes using LASSO-penalized multivariate Cox analysis. We constructed a training set, performed survival analysis on the high-and low-risk groups, and performed validation in the test and full sets. There existed a substantial contrast in the likelihood of survival among the cohorts of high and low risk. An in-depth analysis of the gene mutations associated with tumors was conducted to evaluate the risk of developing cancer. We also performed gene analysis on neoadjuvant chemotherapy. We conducted experimental validation on the key gene UBE2Q1-AS1 in our prognostic signature.

Results

The risk signature we constructed shows significant differences between the high-risk group and the low-risk group. Univariate survival analysis of the eight genes in our signature showed that each gene distinguished between high- and low-risk groups. Sub-group analysis revealed that our risk score differed significantly in tumor stage, age, and gender. The analysis results of the tumor mutation burden (TMB) showed a significant difference in the TMB between the low- and high-risk groups, which had a direct impact on the outcomes. These findings highlight the importance of TMB as a potential prognostic marker in cancer detection and prevention. We analyzed the immune microenvironment and found significant differences in immune function, validation responses, immunotherapy-related positive markers, and critical steps in the tumor immunity cycle between the high- and low-risk groups. We found that the effect of anti-CTLA4 and PD-1 was higher in the high-risk group than in the low-risk group.Gene analysis of neoadjuvant chemotherapy revealed that the treatment effect in the high-risk group was better than in the low-risk group. The key gene UBE2Q1-AS1 in our prognostic signature can significantly influence the cell viability, migration, and proliferation of cancer cells.

Discussion

We established a signature consisting of eight genes constructed from cuproptosis-related lncRNAs that have potential clinical applications for outcomes prediction, diagnosis, and treatment.

A novel machine learning-based programmed cell death-related clinical diagnostic and prognostic model associated with immune infiltration in endometrial cancer

Background

To explore the underlying mechanism of programmed cell death (PCD)-related genes in patients with endometrial cancer (EC) and establish a prognostic model.

Methods

The RNA sequencing data (RNAseq), single nucleotide variation (SNV) data, and corresponding clinical data were downloaded from TCGA. The prognostic PCD-related genes were screened and subjected to consensus clustering analysis. The two clusters were compared by weighted correlation network analysis (WGCNA), immune infiltration analysis, and other analyses. The least absolute shrinkage and selection operator (LASSO) algorithm was used to construct the PCD-related prognostic model. The biological significance of the PCD-related gene signature was evaluated through various bioinformatics methods.

Results

We identified 43 PCD-related genes that were significantly related to prognoses of EC patients, and classified them into two clusters via consistent clustering analysis. Patients in cluster B had higher tumor purity, higher T stage, and worse prognoses compared to those in cluster A. The latter generally showed higher immune infiltration. A prognostic model was constructed using 11 genes (GZMA, ASNS, GLS, PRKAA2, VLDLR, PRDX6, PSAT1, CDKN2A, SIRT3, TNFRSF1A, LRPPRC), and exhibited good diagnostic performance. Patients with high-risk scores were older, and had higher stage and grade tumors, along with worse prognoses. The frequency of mutations in PCD-related genes was correlated with the risk score. LRPPRC, an adverse prognostic gene in EC, was strongly correlated with proliferation-related genes and multiple PCD-related genes. LRPPRC expression was higher in patients with higher clinical staging and in the deceased patients. In addition, a positive correlation was observed between LRPPRC and infiltration of multiple immune cell types.

Conclusion

We identified a PCD-related gene signature that can predict the prognosis of EC patients and offer potential targets for therapeutic interventions.

Novel methylation-related long non-coding RNA clinical outcome prediction method: the clinical phenotype and immune infiltration research in low-grade gliomas

Background

Recent studies have suggested that long non-coding RNAs (lncRNAs) may play crucial role in low-grade glioma; however, the underlying mechanisms linking them to epigenetic methylation remain unclear.

Methods

We downloaded expression level data for regulators associated with N1 methyladenosine (m1A), 5-methyladenine (m5C), and N6 methyladenosine (m6A) (M1A/M5C/M6A) methylation from the Cancer Genome Atlas-low-grade glioma (TCGA-LGG) database. We identified the expression patterns of lncRNAs, and selected methylation-related lncRNAs using Pearson correlation coefficient>0.4. Non-negative matrix dimensionality reduction was then used to determine the expression patterns of the methylation-associated lncRNAs. We constructed a weighted gene co-expression network analysis (WGCNA) network to explore the co-expression networks between the two expression patterns. Functional enrichment of the co-expression network was performed to identify biological differences between the expression patterns of different lncRNAs. We also constructed prognostic networks based on the methylation presence in lncRNAs in low-grade gliomas.

Results

We identified 44 regulators by literature review. Using a correlation coefficient greater than 0.4, we identified 2330 lncRNAs, among which 108 lncRNAs with independent prognostic values were further screened using univariate Cox regression at P< 0.05. Functional enrichment of the co-expression networks revealed that regulation of trans-synaptic signaling, modulation of chemical synaptic transmission, calmodulin binding, and SNARE binding were mostly enriched in the blue module. The calcium and CA2 signaling pathways were associated with different methylation-related long non-coding chains. Using the Least Absolute Shrinkage Selector Operator (LASSO) regression analysis, we analyzed a prognostic model containing four lncRNAs. The model's risk score was 1.12 *AC012063 + 0.74 * AC022382 + 0.32 * AL049712 + 0.16 * GSEC. Gene set variation analysis (GSVA) revealed significant differences in mismatch repair, cell cycle, WNT signaling pathway, NOTCH signaling pathway, Complement and Cascades, and cancer pathways at different GSEC expression levels. Thus, these results suggest that GSEC may be involved in the proliferation and invasion of low-grade glioma, making it a prognostic risk factor for low-grade glioma.

Conclusion

Our analysis identified methylation-related lncRNAs in low-grade gliomas, providing a foundation for further research on lncRNA methylation. We found that GSEC could serve as a candidate methylation marker and a prognostic risk factor for overall survival in low-grade glioma patients. These findings shed light on the underlying mechanisms of low-grade glioma development and may facilitate the development of new treatment strategies.

Endocrine Disrupting Chemicals Influence Hub Genes Associated with Aggressive Prostate Cancer

Prostate cancer (PCa) is one of the most frequently diagnosed cancers among men in the world. Its prevention has been limited because of an incomplete understanding of how environmental exposures to chemicals contribute to the molecular pathogenesis of aggressive PCa. Environmental exposures to endocrine-disrupting chemicals (EDCs) may mimic hormones involved in PCa development. This research aims to identify EDCs associated with PCa hub genes and/or transcription factors (TF) of these hub genes in addition to their protein-protein interaction (PPI) network. We are expanding upon the scope of our previous work, using six PCa microarray datasets, namely, GSE46602, GSE38241, GSE69223, GSE32571, GSE55945, and GSE26126, from the NCBI/GEO, to select differentially expressed genes based on |log2FC| (fold change) ≥ 1 and an adjusted p-value < 0.05. An integrated bioinformatics analysis was used for enrichment analysis (using DAVID.6.8, GO, KEGG, STRING, MCODE, CytoHubba, and GeneMANIA). Next, we validated the association of these PCa hub genes in RNA-seq PCa cases and controls from TCGA. The influence of environmental chemical exposures, including EDCs, was extrapolated using the chemical toxicogenomic database (CTD). A total of 369 overlapping DEGs were identified associated with biological processes, such as cancer pathways, cell division, response to estradiol, peptide hormone processing, and the p53 signaling pathway. Enrichment analysis revealed five up-regulated (NCAPG, MKI67, TPX2, CCNA2, CCNB1) and seven down-regulated (CDK1, CCNB2, AURKA, UBE2C, BUB1B, CENPF, RRM2) hub gene expressions. Expression levels of these hub genes were significant in PCa tissues with high Gleason scores ≥ 7. These identified hub genes influenced disease-free survival and overall survival of patients 60-80 years of age. The CTD studies showed 17 recognized EDCs that affect TFs (NFY, CETS1P54, OLF1, SRF, COMP1) that are known to bind to our PCa hub genes, namely, NCAPG, MKI67, CCNA2, CDK1, UBE2C, and CENPF. These validated differentially expressed hub genes can be potentially developed as molecular biomarkers with a systems perspective for risk assessment of a wide-ranging list of EDCs that may play overlapping and important role(s) in the prognosis of aggressive PCa.

Assessing the Potential Prognostic and Immunological Role of TK1 in Prostate Cancer

Background: It has been reported that thymidine kinase 1 (TK1) was up-regulated in multiple malignancies and participated in the regulation of tumor malignant behavior. However, its specific role in prostate cancer (PCa) remains unclear. Methods: TK1 expression in PCa patients and cell lines was identified via crossover analysis of the public datasets. A series of in vitro experiments and in vivo models was applied to investigate the function of TK1 in PCa. Functional enrichment analyses were further conducted to explore the underlying mechanism. Additionally, TISIDB was applied to explore the correlation between TK1 expression and tumor-infiltrating lymphocytes, immune subtypes, and immune regulatory factors. Results: TK1 expression was significantly up-regulated in PCa patients and cell lines. TK1 ablation inhibited tumor cell proliferation and migration potential, and in vivo experiments showed that TK1 inactivation can significantly restrain tumor growth. Functional enrichment analysis revealed TK1-related hub genes (AURKB, CCNB2, CDC20, CDCA5, CDK1, CENPA, CENPM, KIF2C, NDC80, NUF2, PLK1, SKA1, SPC25, ZWINT), and found that TK1 was closely involved in the regulation of cell cycle. Moreover, elevated mRNA expression of TK1 was related with higher Gleason score, higher clinical stage, higher pathological stage, higher lymph node stage, shorter overall survival, and DFS in PCa patients. Particularly, TK1 represented attenuated expression in C3 PCa and was related with infiltration of CD4+, CD8+ T cells, and dendritic cells as well as immunomodulator expression. Conclusion: Our study indicates that TK1 is a prognostic predictor correlated with poor outcomes of PCa patients, and for the first time represented that TK1 can promote the progression of PCa. Therefore, TK1 may be a potential diagnostic and prognostic biomarker, as well as a therapeutic target for PCa.

Environmental Phenol and Paraben Exposure Risks and Their Potential Influence on the Gene Expression Involved in the Prognosis of Prostate Cancer

Prostate cancer (PCa) is one of the leading malignant tumors in US men. The lack of understanding of the molecular pathology on the risk of food supply chain exposures of environmental phenol (EP) and paraben (PB) chemicals limits the prevention, diagnosis, and treatment options. This research aims to utilize a risk assessment approach to demonstrate the association of EP and PB exposures detected in the urine samples along with PCa in US men (NHANES data 2005−2015). Further, we employ integrated bioinformatics to examine how EP and PB exposure influences the molecular pathways associated with the progression of PCa. The odds ratio, multiple regression model, and Pearson coefficients were used to evaluate goodness-of-fit analyses. The results demonstrated associations of EPs, PBs, and their metabolites, qualitative and quantitative variables, with PCa. The genes responsive to EP and PB exposures were identified using the Comparative Toxicogenomic Database (CTD). DAVID.6.8, GO, and KEGG enrichment analyses were used to delineate their roles in prostate carcinogenesis. The plug-in CytoHubba and MCODE completed identification of the hub genes in Cytoscape software for their roles in the PCa prognosis. It was then validated by using the UALCAN database by evaluating the expression levels and predictive values of the identified hub genes in prostate cancer prognosis using TCGA data. We demonstrate a significant association of higher levels of EPs and PBs in the urine samples, categorical and numerical confounders, with self-reported PCa cases. The higher expression levels of the hub genes (BUB1B, TOP2A, UBE2C, RRM2, and CENPF) in the aggressive stages (Gleason score > 8) of PCa tissues indicate their potential role(s) in the carcinogenic pathways. Our results present an innovative approach to extrapolate and validate hub genes responsive to the EPs and PBs, which may contribute to the severity of the disease prognosis, especially in the older population of US men.

Integration of quantitative diffusion kurtosis imaging and prostate specific antigen in differential diagnostic of prostate cancer

This study aimed to evaluate the diagnostic performance of diffusion kurtosis imaging (DKI) and prostate-specific antigen (PSA) biomarkers in differentiating prostate cancer (PCa) and benign prostatic hyperplasia (BPH).A total of 43 cases of prostate diseases verified by pathology were enrolled in the present study. These cases were assigned to the BPH group (n = 20, 68.85±10.81 years old) and PCa group (n = 23, 74.13 ± 7.37 years old). All patients underwent routine prostate magnetic resonance imaging and DKI examinations, and the mean diffusivity (MD), mean kurtosis (MK), and fractional anisotropy (FA) values were calculated. Three serum indicators (PSA, free PSA [fPSA], and f/t PSA) were collected. We used univariate logistic regression to analyze the above quantitative parameters between the 2 groups, and the independent factors were further incorporated into the multivariate logistic regression model. The area under the receiver operating characteristic curve (AUC) was used to evaluate the diagnostic efficacy of the single indicator and combined model.The difference in PSA, f/t PSA, MK, and FA between PCa and BPH was statistically significant (P < .05). The AUC for the combined model (f/t PSA, MK, and FA) of 0.972 (95% confidence interval [CI]: 0.928, 1.000) was higher than the AUC of 0.902 (95% CI: 0.801, 1.000) for f/t PSA, 0.833 (95% CI: 0.707, 0.958) for MK, and 0.807 (95% CI: 0.679, 0.934) for FA.The MK and FA values for DKI and f/t PSA effectively identify PCa and BPH, compared to the PSA indicators. Combining DKI and PSA derivatives can further improve the diagnosis efficiency and might help in the clinical setting.

Identification of crucial genes and pathways associated with prostate cancer in multiple databases

OBJECTIVE

Prostate cancer (PCa) is a malignant neoplasm of the urinary system. This study aimed to use bioinformatics to screen for core genes and biological pathways related to PCa.

METHODS

The GSE5957 gene expression profiles were obtained from the Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs). Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the DEGs were constructed by R language. Furthermore, protein-protein interaction (PPI) networks were generated to predict core genes. The expression levels of core genes were examined in the Tumor Immune Estimation Resource (TIMER) and Oncomine databases. The cBioPortal tool was used to study the co-expression and prognostic factors of the core genes. Finally, the core genes of signaling pathways were determined using gene set enrichment analysis (GSEA).

RESULTS

Overall, 874 DEGs were identified. Hierarchical clustering analysis revealed that these 24 core genes have significant association with carcinogenesis and development. LONRF1, CDK1, RPS18, GNB2L1 (RACK1), RPL30, and SEC61A1 directly related to the recurrence and prognosis of PCa.

CONCLUSIONS

This study identified the core genes and pathways in PCa and provides candidate targets for diagnosis, prognosis, and treatment.

Screening of differentially expressed genes and identification of AMACR as a prognostic marker in prostate cancer

Prostate cancer, the second most common cancer found in male over the world, was estimated to have 191,930 new cases and 33,330 deaths in 2020 in the United States. Prostate cancer is very common in male, about 12.1% of men will acquire this cancer in their lifetime, and a higher risk was reported in older men and African American men. Gene deregulations have been found to be extensively associated with cancer development. To gain further insight into how gene deregulation affects prostate cancer, we analysed three gene profiling datasets of prostate cancer from Gene Expression Omnibus (GEO) applying bioinformatic tools in our study. Firstly, we identified common differently expressed genes (DEGs) shared by the three gene profiling datasets, constructed protein-protein interaction network and determined top 10 hub genes. Further DEGs validation in TCGA and Human Protein Atlas Database identified AMACR as the core gene. We then analysed the role of AMACR in prostate cancer cell lines and found that AMACR-knockdown resulted in the decreased cell proliferation and increased apoptosis. These results suggest an oncogenic role of AMACR in prostate cancer, and it could be a potential biomarker for the diagnosis of prostate cancer.

Identification of UBE2C as hub gene in driving prostate cancer by integrated bioinformatics analysis

BACKGROUND

The aim of this study was to identify novel genes in promoting primary prostate cancer (PCa) progression and to explore its role in the prognosis of prostate cancer.

METHODS

Four microarray datasets containing primary prostate cancer samples and benign prostate samples were downloaded from Gene Expression Omnibus (GEO), then differentially expressed genes (DEGs) were identified by R software (version 3.6.2). Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed to identify the function of DEGs. Using STRING and Cytoscape (version 3.7.1), we constructed a protein-protein interaction (PPI) network and identified the hub gene of prostate cancer. Clinical data on GSE70770 and TCGA was collected to show the role of hub gene in prostate cancer progression. The correlations between hub gene and clinical parameters were also indicated by cox regression analysis. Gene Set Enrichment Analysis (GSEA) was performed to highlight the function of Ubiquitin-conjugating enzyme complex (UBE2C) in prostate cancer.

RESULTS

243 upregulated genes and 298 downregulated genes that changed in at least two microarrays have been identified. GO and KEGG analysis indicated significant changes in the oxidation-reduction process, angiogenesis, TGF-beta signaling pathway. UBE2C, PDZ-binding kinase (PBK), cyclin B1 (CCNB1), Cyclin-dependent kinase inhibitor 3 (CDKN3), topoisomerase II alpha (TOP2A), Aurora kinase A (AURKA) and MKI67 were identified as the candidate hub genes, which were all correlated with prostate cancer patient' disease-free survival in TCGA. In fact, only UBE2C was highly expressed in prostate cancer when compared with benign prostate tissue in TCGA and the expression of UBE2C was also in parallel with the Gleason score of prostate cancer. Cox regression analysis has indicated UBE2C could function as the independent prognostic factor of prostate cancer. GSEA showed UBE2C had played an important role in the pathway of prostate cancer, such as NOTCH signaling pathway, WNT-β-catenin signaling pathway.

CONCLUSIONS

UBE2C was pivotal for the progression of prostate cancer and the level of UBE2C was important to predict the prognosis of patients.