Transcriptome analysis of differentially expressed genes and pathways associated with mitoxantrone treatment prostate cancer

Prognostic signatures associated with high infiltration of Tregs in bone metastatic prostate cancer

Metastatic cancer especially bone metastasis (BM) is the lethal end-stage of castration-resistant prostate cancer (CRPC). To understand the possible molecular mechanisms underlying the development of the distant metastasis is of potential clinical value. We sought to identify differentially expressed genes between patient-matched primary and bone metastatic CRPC tumors. Functional enrichment, protein-protein interaction networks, and survival analysis of DEGs were performed. DEGs with a prognostic value considered as candidate genes were evaluated, followed by genetic analysis of tumor infiltrating immune cells based on Wilcoxon test and immunofluorescence identification. Expression profiles analysis showed that 381 overlapping genes were screened as differentially expressed genes (DEGs), of which 16 DEGs were randomly selected to be validated and revealed that most of these genes showed a transcriptional profile similar to that seen in the datasets (Pearson’s r = 0.76). Six core genes were found to be involved in regulation of extracellular matrix receptor interaction and chemotactic activity, and four of them were significantly correlated with the survival of PCa patients with bone metastases. Immune infiltration analysis showed that the expressions levels of COL3A1, RAC1, FN1, and SDC2 in CD4+T cells were significantly higher than those in tumor cells, especially regulatory T cell infiltration was significantly increased in BM tumors. We analyzed gene expression signatures specifically associated with the development of bone metastases of CRPC patients. Characterization of genes associated with BM of mCRPC is critical for identification of predictive biomarkers and potential therapeutic targets.

Selective killing of homologous recombination-deficient cancer cell lines by inhibitors of the RPA:RAD52 protein-protein interaction

Synthetic lethality is a successful strategy employed to develop selective chemotherapeutics against cancer cells. Inactivation of RAD52 is synthetically lethal to homologous recombination (HR) deficient cancer cell lines. Replication protein A (RPA) recruits RAD52 to repair sites, and the formation of this protein-protein complex is critical for RAD52 activity. To discover small molecules that inhibit the RPA:RAD52 protein-protein interaction (PPI), we screened chemical libraries with our newly developed Fluorescence-based protein-protein Interaction Assay (FluorIA). Eleven compounds were identified, including FDA-approved drugs (quinacrine, mitoxantrone, and doxorubicin). The FluorIA was used to rank the compounds by their ability to inhibit the RPA:RAD52 PPI and showed mitoxantrone and doxorubicin to be the most effective. Initial studies using the three FDA-approved drugs showed selective killing of BRCA1-mutated breast cancer cells (HCC1937), BRCA2-mutated ovarian cancer cells (PE01), and BRCA1-mutated ovarian cancer cells (UWB1.289). It was noteworthy that selective killing was seen in cells known to be resistant to PARP inhibitors (HCC1937 and UWB1 SYr13). A cell-based double-strand break (DSB) repair assay indicated that mitoxantrone significantly suppressed RAD52-dependent single-strand annealing (SSA) and mitoxantrone treatment disrupted the RPA:RAD52 PPI in cells. Furthermore, mitoxantrone reduced radiation-induced foci-formation of RAD52 with no significant activity against RAD51 foci formation. The results indicate that the RPA:RAD52 PPI could be a therapeutic target for HR-deficient cancers. These data also suggest that RAD52 is one of the targets of mitoxantrone and related compounds.

In silico identification of MAPK14-related lncRNAs and assessment of their expression in breast cancer samples

Mitogen-activated protein kinase (MAP kinase) pathways participate in regulation of several cellular processes involved in breast carcinogenesis. A number of non-coding RNAs including both microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) regulate or being regulated by MAPKs. We performed an in-silico method for identification of MAPKs with high number of interactions with miRNAs and lncRNAs. Bioinformatics approaches revealed that MAPK14 ranked first among MAPKs. Subsequently, we identified miRNAs and lncRNAs that were predicted to be associated with MAPK14. Finally, we selected four lncRNAs with higher predicted scores (NORAD, HCG11, ZNRD1ASP and TTN-AS1) and assessed their expression in 80 breast cancer tissues and their adjacent non-cancerous tissues (ANCTs). Expressions of HCG11 and ZNRD1ASP were lower in tumoral tissues compared with ANCTs (P values < 0.0001). However, expression levels of MAPK14 and NORAD were not significantly different between breast cancer tissues and ANCTs. A significant association was detected between expression of HCG11 and estrogen receptor (ER) status in a way that tumors with up-regulation of this lncRNA were mostly ER negative (P value = 0.04). Expressions of ZNRD1ASP and HCG11 were associated with menopause age and breast feeding duration respectively (P values = 0.02 and 0.04 respectively). There was a trend towards association between ZNRD1ASP expression and patients' age of cancer diagnosis. Finally, we detected a trend toward association between expression of NORAD and history of hormone replacement therapy (P value = 0.06). Expression of MAPK14 was significantly higher in grade 1 tumors compared with grade 2 tumors (P value = 0.02). Consequently, the current study provides evidences for association between lncRNA expressions and reproductive factors or tumor features.

Expression profile-based screening for critical genes reveals S100A4, ACKR3 and CDH1 in docetaxel-resistant prostate cancer cells

Docetaxel is a first-line anticancer drug widely used in the treatment of advanced prostate cancer. However, its therapeutic efficacy is limited by its side effects and the development of chemoresistance by the tumor. Using a gene differential expression microarray, we identified 449 genes differentially expressed in docetaxel-resistant DU145 and PC3 cell lines as compared to docetaxel-sensitive controls. Moreover, western blotting and immunohistochemistry revealed altered expression of S100A4, ACKR3 and CDH1in clinical tumor samples. Cytoscape software was used to investigate the relationship between critical proteins and their signaling transduction networks. Functional and pathway enrichment analyses revealed that these signaling pathways were closely related to cellular proliferation, cell adhesion, cell migration and metastasis. In addition, ACKR3 knockout using the crispr/cas9 method andS100A4knockdownusing targeted shRNA exerted additive effects suppressing cancer cell proliferation and migration. This exploratory analysis provides information about potential candidate genes. It also provides new insight into the molecular mechanism underlying docetaxel-resistance in androgen-independent prostate cancer and highlights potential targets to improve therapeutic outcomes.

ZeOncoTest: Refining and Automating the Zebrafish Xenograft Model for Drug Discovery in Cancer

The xenograft of human cancer cells in model animals is a powerful tool for understanding tumor progression and metastatic potential. Mice represent a validated host, but their use is limited by the elevated experimental costs and low throughput. To overcome these restrictions, zebrafish larvae might represent a valuable alternative. Their small size and transparency allow the tracking of transplanted cells. Therefore, tumor growth and early steps of metastasis, which are difficult to evaluate in mice, can be addressed. In spite of its advantages, the use of this model has been hindered by lack of experimental homogeneity and validation. Considering these facts, the aim of our work was to standardize, automate, and validate a zebrafish larvae xenograft assay with increased translatability and higher drug screening throughput. The ZeOncoTest reliability is based on the optimization of different experimental parameters, such as cell labeling, injection site, automated individual sample image acquisition, and analysis. This workflow implementation finally allows a higher precision and experimental throughput increase, when compared to previous reports. The approach was validated with the breast cancer cell line MDA-MB-231, the colorectal cancer cells HCT116, and the prostate cancer cells PC3; and known drugs, respectively RKI-1447, Docetaxel, and Mitoxantrone. The results recapitulate growth and invasion for all tested tumor cells, along with expected efficacy of the compounds. Finally, the methodology has proven useful for understanding specific drugs mode of action. The insights gained bring a step further for zebrafish larvae xenografts to enter the regulated preclinical drug discovery path.

Transcriptome analysis of differentially expressed genes and pathways associated with mitoxantrone treatment prostate cancer

The global physiological function of specifically expressed genes of mitoxantrone (MTX)‐resistant prostate cancer (PCa) is unclear. In this study, gene expression pattern from microarray data was investigated for identifying differentially expressed genes (DEGs) in MTX‐resistant PCa xenografts. Human PCa cell lines DU145 and PC3 were cultured in vitro and xenografted into severe combined immunodeficiency (SCID) mice, treated with MTX intragastrically, three times a week until all mice relapsed. Gene expression profiles of the xenografts from castrated mice were performed with Affymetrix human whole genomic oligonucleotide microarray. The Cytoscape software was used to investigate the relationship between proteins and the signalling transduction network. A total of 355 overlapping genes were differentially expressed in MTX‐resistant DU145R and PC3R xenografts. Of these, 16 genes were selected to be validated by quantitative real‐time PCR (qRT‐PCR) in these xenografts, and further tested in a set of formalin‐fixed, paraffin‐embedded and optimal cutting temperature (OCT) clinical tumour samples. Functional and pathway enrichment analyses revealed that these DEGs were closely related to cellular activity, androgen synthesis, DNA damage and repair, also involved in the ERK/MAPK, PI3K/serine‐threonine protein kinase, also known as protein kinase B, PKB (AKT) and apoptosis signalling pathways. This exploratory analysis provides information about potential candidate genes and may bring new insights into the molecular cascade involvement in MTX‐resistant PCa.