Identification of genes associated with castration‑resistant prostate cancer by gene expression profile analysis

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.

Identification of Critical Molecular Factors and Side Effects Underlying the Response to Thalicthuberine in Prostate Cancer: A Systems Biology Approach

Background

Uncontrolled mitosis of cancer cells and resistance cells to chemotherapy drugs are the challenges of prostate cancer. Thalicthuberine causes a mitotic arrest and a reduction of the effects of drug resistance, resulting in cell death. In this study, we applied bioinformatics and computational biology methods to identify functional pathways and side effects in response to Thalicthuberine in prostate cancer patients.

Methods

Microarray data were retrieved from Gene Expression Omnibus (GEO), and protein-protein interactions and gene regulatory networks were constructed, using the Cytoscape software. The critical genes and molecular mechanisms in response to Thalicthuberine and its side effects were identified, using the Cytoscape software and WebGestalt server, respectively. Finally, GEPIA2 was used to predict the relationship between critical genes and prostate cancer.

Results

The POLQ, EGR1, CDKN1A, FOS, MDM2, CDC20, CCNB1, and CCNB2 were identified as critical genes in response to this drug. The functional mechanisms of Thalicthuberine include a response to oxygen levels, toxic substances and immobilization stress, cell cycle regulation, regeneration, the p53 signaling pathway, the action of the parathyroid hormone, and the FoxO signaling pathway. Besides, the drug has side effects including muscle cramping, abdominal pains, paresthesia, and metabolic diseases.

Conclusion

Our model suggested newly predicted crucial genes, molecular mechanisms, and possible side effects of this drug. However, further studies are required.

Identification of a 12-Gene Signature and Hub Genes Involved in Kidney Wilms Tumor via Integrated Bioinformatics Analysis

Wilms tumor (WT), also known as nephroblastoma, is a rare primary malignancy in all kinds of tumor. With the development of second-generation sequencing, the discovery of new tumor markers and potential therapeutic targets has become easier. This study aimed to explore new WT prognostic biomarkers. In this study, WT-miRNA datasets GSE57370 and GSE73209 were selected for expression profiling to identify differentially expressed genes. The key gene miRNA, namely hsa-miR-30c-5p, was identified by overlapping, and the target gene of candidate hsa-miR-30c-5p was predicted using an online database. Furthermore, 384 genes were obtained by intersecting them with differentially expressed genes in the TARGET-WT database, and the genes were analyzed for pathway and functional enrichment. Kaplan–Meier survival analysis of the 384 genes yielded a total of 25 key genes associated with WT prognosis. Subsequently, a prediction model with 12 gene signatures (BCL6, CCNA1, CTHRC1, DGKD, EPB41L4B, ERRFI1, LRRC40, NCEH1, NEBL, PDSS1, ROR1, and RTKN2) was developed. The model had good predictive power for the WT prognosis at 1, 3, and 5 years (AUC: 0.684, 0.762, and 0.774). Finally, ERRFI1 (hazard ratios [HR] = 1.858, 95% confidence intervals [CI]: 1.298–2.660) and ROR1 (HR = 0.780, 95% CI: 0.609–0.998) were obtained as independent predictors of prognosis in WT patients by single, multifactorial Cox analysis.

Blood-based gene expression signature associated with metastatic castrate-resistant prostate cancer patient response to abiraterone plus prednisone or enzalutamide

BACKGROUND

Precision medicine approaches for managing patients with metastatic castrate-resistant prostate cancer (mCRPC) are lacking. Non-invasive approaches for molecular monitoring of disease are urgently needed, especially for patients suffering from bone metastases for whom tissue biopsy is challenging. Here we utilized baseline blood samples to identify mCRPC patients most likely to benefit from abiraterone plus prednisone (AAP) or enzalutamide.

METHODS

Baseline blood samples were collected for circulating tumor cell (CTC) enumeration and qPCR-based gene expression analysis from 51 men with mCRPC beginning treatment with abiraterone or enzalutamide.

RESULTS

Of 51 patients (median age 68 years [51-82]), 22 received AAP (abiraterone 1000 mg/day plus prednisone 10 mg/day) and 29 received enzalutamide (160 mg/day). The cohort was randomly divided into training (n = 37) and test (n = 14) sets. Baseline clinical variables (Gleason score, PSA, testosterone, and hemoglobin), CTC count, and qPCR-based gene expression data for 141 genes/isoforms in CTC-enriched blood were analyzed with respect to overall survival (OS). Genes with expression most associated with OS included MSLN, ARG2, FGF8, KLK3, ESRP2, NPR3, CCND1, and WNT5A. Using a Cox-elastic net model for our test set, the 8-gene expression signature had a c-index of 0.87 (95% CI [0.80, 0.94]) and was more strongly associated with OS than clinical variables or CTC count alone, or a combination of the three variables. For patients with a low-risk vs. high-risk gene expression signature, median OS was not reached vs. 18 months, respectively (HR 5.32 [1.91-14.80], p = 0.001). For the subset of 41 patients for whom progression-free survival (PFS) data was available, the median PFS for patients with a low-risk vs high-risk gene expression signature was 20 vs. 5 months, respectively (HR 2.95 [1.46-5.98], p = 0.003).

CONCLUSIONS

If validated in a larger prospective study, this test may predict patients most likely to benefit from second-generation antiandrogen therapy.

Stemness Related Genes Revealed by Network Analysis Associated With Tumor Immune Microenvironment and the Clinical Outcome in Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) is one of the leading fatal malignancy with high morbidity and mortality worldwide. However, due to its complicated mechanism and lack of effective clinical therapeutics, early diagnosis and prognosis are still unsatisfactory. Most of the previous studies focused on cancer stem cells (CSCs), the relationship between cancer stemness (stem-like characteristics) and anti-tumor immunity has not been clearly revealed. Therefore, this study aimed to comprehensively analyze the role of cancer stemness and tumor microenvironment (TME) in LUAD using weighted gene co-expression network analysis (WGCNA). We constructed a gene co-expression network, identified key modules, and hub genes, and further explored the relationship between hub gene expression and cancer immunological characteristics through a variety of algorithms, including Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) and Gene Set Enrichment Analysis (GSEA). The hub genes were renamed stemness related genes (SRGs), whose functions were examined at the transcription and protein levels through survival analysis with additional samples, Oncomine database, immunohistochemistry, single cell RNA sequencing (scRNA-seq) and single-sample Gene Set Enrichment Analysis (ssGSEA). Subsequently, Tumor Immune Dysfunction and Exclusion (TIDE) and Connectivity Map (CMap) were implemented for treatment and prognosis analyses. As a result, 15 co-expressed SRGs (CCNA2, CCNB1, CDC20, CDCA5, CDCA8, FEN1, KIF2C, KPNA2, MCM6, NUSAP1, RACGAP1, RRM2, SPAG5, TOP2A, and TPX2) were identified. The overexpression of which was discovered to be associated with reduced immune infiltration in LUAD. It was discovered that there was a general negative correlation between cancer stemness and immunity. The expression of SRGs could probably affect our tumor occurrence, progression, the efficacy of chemotherapy and immunotherapy, and clinical outcomes. In conclusion, the 15 SRGs reported in our study may be used as potential candidate biomarkers for prognostic indicators and therapeutic targets after further validation.

Comprehensive analysis of competing endogenous RNA network in Wilms tumor based on the TARGET database

Background

Wilms tumor (WT) was the most common malignant tumor of urinary system in children. With the advances in gene sequencing, research of molecular mechanism of WT tumor was gradually increasing. However, few studies have explored the influence of competing endogenous RNA (ceRNA) in WT. Accordingly, we aimed to explore the mechanisms of ceRNA co-expression network in WT.

Methods

A total of 6 non-tumor controls and 127 WT patients’ RNA-seq data combined with clinical data was acquired from Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database. Differentially expressed lncRNA, miRNA and mRNA between WT tissues and normal tissues were analyzed using “edgeR” package in R software. Weighted gene co-expression network analysis (WGCNA) was utilized to construct the ceRNA co-expression network while Molecular Complex Detection (MCODE) algorithm was used to extract the pivotal sub-network. Function annotation of mRNA was performed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Survival analysis was then conducted based on long-rank test and Kaplan-Meier curves using the survival package.

Results

By applying the “edgeR” package in R, the transcriptome expression data of 127 WT tissues with 6 normal tissues were normalized. Moreover, 146 DElncRNAs, 62 DEmiRNAs, 287 DEmRNAs of them were involved in ceRNA network after applying WGCNA. According to MCODE, we identified that the interactions between LINC002253 (lncRNA) and TRIM71 (mRNA) was mediated by hsa-mir-301a and hsa-mir-301b (miRNA). Furthermore, we detected 13 DElncRNAs which were significantly associated with the progression of WT.

Conclusions

We used WGCNA method to construct the WT ceRNA network for the first time. TRIM71 was identified to be the most closely related genes involved in hub sub-network by MCDOE, suggesting it might act as a new drug target and prognostic factor based on our comprehensive results.

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

Abstract

The morbidity and mortality of prostate carcinoma has increased in recent years and has become the second most common ale malignant carcinoma worldwide. The interaction mechanisms between different genes and signaling pathways, however, are still unclear.

Methods

Variation analysis of GSE38241, GSE69223, GSE46602 and GSE104749 were realized by GEO2R in Gene Expression Omnibus database. Function enrichment was analyzed by DAVID.6.8. Furthermore, the PPI network and the significant module were analyzed by Cytoscape, STRING and MCODE.GO. Pathway analysis showed that the 20 candidate genes were closely related to mitosis, cell division, cell cycle phases and the p53 signaling pathway. A total of six independent prognostic factors were identified in GSE21032 and TCGA PRAD. Oncomine database and The Human Protein Atlas were applied to explicit that six core genes were over expression in prostate cancer compared to normal prostate tissue in the process of transcriptional and translational. Finally, gene set enrichment were performed to identified the related pathway of core genes involved in prostate cancer.

Result

Hierarchical clustering analysis revealed that these 20 core genes were mostly related to carcinogenesis and development. CKS2, TK1, MKI67, TOP2A, CCNB1 and RRM2 directly related to the recurrence and prognosis of prostate cancer. This result was verified by TCGA database and GSE21032.

Conclusion

These core genes play a crucial role in tumor carcinogenesis, development, recurrence, metastasis and progression. Identifying these genes could help us to understand the molecular mechanisms and provide potential biomarkers for the diagnosis and treatment of prostate cancer.

Weighted gene co‑expression network analysis for identifying hub genes in association with prognosis in Wilms tumor

Wilms tumor (WT) is the most common type of renal malignancy in children. Survival rates are low and high-risk WT generally still carries a poor prognosis. To better elucidate the pathogenesis and tumorigenic pathways of high-risk WT, the present study presents an integrated analysis of RNA expression profiles of high-risk WT to identify predictive molecular biomarkers, for the improvement of therapeutic decision-making. mRNA sequence data from high-risk WT and adjacent normal samples were downloaded from The Cancer Genome Atlas to screen for differentially expressed genes (DEGs) using R software. From 132 Wilms tumor samples and six normal samples, 2,089 downregulated and 941 upregulated DEGs were identified. In order to identify hub DEGs that regulate target genes, weighted gene co-expression network analysis (WGCNA) was used to identify 11 free-scale gene co-expressed clusters. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were annotated using KEGG Orthology Based Annotation System annotation of different module genes. The Search Tool for the Retrieval of Interacting Genes was used to construct a protein-protein interaction network for the identified DEGs, and the hub genes of WGCNA modules were identified using the Cytohubb plugin with Cytoscape software. Survival analysis was subsequently performed to highlight hub genes with a clinical signature. The present results suggest that epidermal growth factor, cyclin dependent kinase 1, endothelin receptor type A, nerve growth factor receptor, opa-interacting protein 5, NDC80 kinetochore complex component and cell division cycle associated 8 are essential to high-risk WT pathogenesis, and they are closely associated with clinical prognosis.

Identification of the Transcription Factor Relationships Associated with Androgen Deprivation Therapy Response and Metastatic Progression in Prostate Cancer

Background: Patients with locally advanced or recurrent prostate cancer typically undergo androgen deprivation therapy (ADT), but the benefits are often short-lived and the responses variable. ADT failure results in castration-resistant prostate cancer (CRPC), which inevitably leads to metastasis. We hypothesized that differences in tumor transcriptional programs may reflect differential responses to ADT and subsequent metastasis. Results: We performed whole transcriptome analysis of 20 patient-matched Pre-ADT biopsies and 20 Post-ADT prostatectomy specimens, and identified two subgroups of patients (high impact and low impact groups) that exhibited distinct transcriptional changes in response to ADT. We found that all patients lost the AR-dependent subtype (PCS2) transcriptional signatures. The high impact group maintained the more aggressive subtype (PCS1) signal, while the low impact group more resembled an AR-suppressed (PCS3) subtype. Computational analyses identified transcription factor coordinated groups (TFCGs) enriched in the high impact group network. Leveraging a large public dataset of over 800 metastatic and primary samples, we identified 33 TFCGs in common between the high impact group and metastatic lesions, including SOX4/FOXA2/GATA4, and a TFCG containing JUN, JUNB, JUND, FOS, FOSB, and FOSL1. The majority of metastatic TFCGs were subsets of larger TFCGs in the high impact group network, suggesting a refinement of critical TFCGs in prostate cancer progression. Conclusions: We have identified TFCGs associated with pronounced initial transcriptional response to ADT, aggressive signatures, and metastasis. Our findings suggest multiple new hypotheses that could lead to novel combination therapies to prevent the development of CRPC following ADT.

Prognostic Genes of Breast Cancer Identified by Gene Co-expression Network Analysis

Breast cancer is one of the most common malignancies. The molecular mechanisms of its pathogenesis are still to be investigated. The aim of this study was to identify the potential genes associated with the progression of breast cancer. Weighted gene co-expression network analysis (WGCNA) was used to construct free-scale gene co-expression networks to explore the associations between gene sets and clinical features, and to identify candidate biomarkers. The gene expression profiles of GSE1561 were selected from the Gene Expression Omnibus (GEO) database. RNA-seq data and clinical information of breast cancer from TCGA were used for validation. A total of 18 modules were identified via the average linkage hierarchical clustering. In the significant module (R² = 0.48), 42 network hub genes were identified. Based on the Cancer Genome Atlas (TCGA) data, 5 hub genes (CCNB2, FBXO5, KIF4A, MCM10, and TPX2) were correlated with poor prognosis. Receiver operating characteristic (ROC) curve validated that the mRNA levels of these 5 genes exhibited excellent diagnostic efficiency for normal and tumor tissues. In addition, the protein levels of these 5 genes were also significantly higher in tumor tissues compared with normal tissues. Among them, CCNB2, KIF4A, and TPX2 were further upregulated in advanced tumor stage. In conclusion, 5 candidate biomarkers were identified for further basic and clinical research on breast cancer with co-expression network analysis.