Multiomics analysis of serial PARP inhibitor treated metastatic TNBC inform on rational combination therapies

In a pilot study, we evaluated the feasibility of real-time deep analysis of serial tumor samples from triple negative breast cancer patients to identify mechanisms of resistance and treatment opportunities as they emerge under therapeutic stress engendered by poly-ADP-ribose polymerase (PARP) inhibitors (PARPi). In a BRCA-mutant basal breast cancer exceptional long-term survivor, a striking tumor destruction was accompanied by a marked infiltration of immune cells containing CD8 effector cells, consistent with pre-clinical evidence for association between STING mediated immune activation and benefit from PARPi and immunotherapy. Tumor cells in the exceptional responder underwent extensive protein network rewiring in response to PARP inhibition. In contrast, there were minimal changes in the ecosystem of a luminal androgen receptor rapid progressor, likely due to indifference to the effects of PARP inhibition. Together, identification of PARPi-induced emergent changes could be used to select patient specific combination therapies, based on tumor and immune state changes.

Analyzing mRNAsi-Related Genes Identifies Novel Prognostic Markers and Potential Drug Combination for Patients with Basal Breast Cancer

Basal breast cancer subtype is the worst prognosis subtypes among all breast cancer subtypes. Recently, a new tumor stemness index-mRNAsi is found to be able to measure the degree of oncogenic differentiation of tissues. The mRNAsi involved in a variety of cancer processes is derived from the innovative application of one-class logistic regression (OCLR) machine learning algorithm to the whole genome expression of various stem cells and tumor cells. However, it is largely unknown about mRNAsi in basal breast cancer. Here, we find that basal breast cancer carries the highest mRNAsi among all four subtypes of breast cancer, especially 385 mRNAsi-related genes are positively related to the high mRNAsi value in basal breast cancer. This high mRNAsi is also closely related to active cell cycle, DNA replication, and metabolic reprogramming in basal breast cancer. Intriguingly, in the 385 genes, TRIM59, SEPT3, RAD51AP1, and EXO1 can act as independent protective prognostic factors, but CTSF and ABHD4B can serve as independent bad prognostic factors in patients with basal breast cancer. Remarkably, we establish a robust prognostic model containing the 6 mRNAsi-related genes that can effectively predict the survival rate of patients with the basal breast cancer subtype. Finally, the drug sensitivity analysis reveals that some drug combinations may be effectively against basal breast cancer via targeting the mRNAsi-related genes. Taken together, our study not only identifies novel prognostic biomarkers for basal breast cancers but also provides the drug sensitivity data by establishing an mRNAsi-related prognostic model.

ISGylation drives basal breast tumour progression by promoting EGFR recycling and Akt signalling

ISG15 is an ubiquitin-like modifier that is associated with reduced survival rates in breast cancer patients. The mechanism by which ISG15 achieves this however remains elusive. We demonstrate that modification of Rab GDP-Dissociation Inhibitor Beta (GDI2) by ISG15 (ISGylation) alters endocytic recycling of the EGF receptor (EGFR) in non-interferon stimulated cells using CRISPR-knock out models for ISGylation. By regulating EGFR trafficking, ISGylation enhances EGFR recycling and sustains Akt-signalling. We further show that Akt signalling positively correlates with levels of ISG15 and its E2-ligase in basal breast cancer cohorts, confirming the link between ISGylation and Akt signalling in human tumours. Persistent and enhanced Akt activation explains the more aggressive tumour behaviour observed in human breast cancers. We show that ISGylation can act as a driver of tumour progression rather than merely being a bystander.

A systems-based framework to computationally describe putative transcription factors and signaling pathways regulating glycan biosynthesis

Glycosylation is a common posttranslational modification, and glycan biosynthesis is regulated by a set of glycogenes. The role of transcription factors (TFs) in regulating the glycogenes and related glycosylation pathways is largely unknown. In this work, we performed data mining of TF–glycogene relationships from the Cistrome Cancer database (DB), which integrates chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA-Seq data to constitute regulatory relationships. In total, we observed 22,654 potentially significant TF–glycogene relationships, which include interactions involving 526 unique TFs and 341 glycogenes that span 29 the Cancer Genome Atlas (TCGA) cancer types. Here, TF–glycogene interactions appeared in clusters or so-called communities, suggesting that changes in single TF expression during both health and disease may affect multiple carbohydrate structures. Upon applying the Fisher’s exact test along with glycogene pathway classification, we identified TFs that may specifically regulate the biosynthesis of individual glycan types. Integration with Reactome DB knowledge provided an avenue to relate cell-signaling pathways to TFs and cellular glycosylation state. Whereas analysis results are presented for all 29 cancer types, specific focus is placed on human luminal and basal breast cancer disease progression. Overall, the article presents a computational approach to describe TF–glycogene relationships, the starting point for experimental system-wide validation.

Syndecan-1 Promotes Angiogenesis in Triple-Negative Breast Cancer through the Prognostically Relevant Tissue Factor Pathway and Additional Angiogenic Routes

Triple-negative breast cancer (TNBC) is characterized by increased angiogenesis, metastasis, and poor survival. Dysregulation of the cell surface heparan sulfate proteoglycan and signaling co-receptor Syndecan-1 is linked to poor prognosis. To study its role in angiogenesis, we silenced Syndecan-1 in TNBC cell lines using a 3D human umbilical vein endothelial cell (HUVEC) co-culture system. Syndecan-1 siRNA depletion in SUM-149, MDA-MB-468, and MDA-MB-231 cells decreased HUVEC tubule network formation. Angiogenesis array revealed reduced VEGF-A and tissue factor (TF) in the Syndecan-1-silenced secretome. qPCR independently confirmed altered expression of F3, F7, F2R/PAR1, F2RL1/PAR2, VEGF-A, EDN1, IGFBP1, and IGFBP2 in SUM-149, MDA-MB-231, and MDA-MB-468 cells. ELISA revealed reduced secreted endothelin-1 (SUM-149, MDA-MB-468) and TF (all cell lines) upon Syndecan-1 depletion, while TF pathway inhibitor treatment impaired angiogenesis. Survival analysis of 3951 patients demonstrated that high expression of F3 and F7 are associated with better relapse-free survival, whereas poor survival was observed in TNBC and p53 mutant basal breast cancer (F3) and in ER-negative and HER2-positive breast cancer (F2R, F2RL1). STRING protein network analysis revealed associations of Syndecan-1 with VEGF-A and IGFBP1, further associated with the TF and ET-1 pathways. Our study suggests that TNBC Syndecan-1 regulates angiogenesis via the TF and additional angiogenic pathways and marks its constituents as novel prognostic markers and therapeutic targets.

Triple negative aggressive phenotype controlled by miR-135b and miR-365: new theranostics candidates

Triple negative breast cancer (TNBC) accounts for about a fifth of all breast cancers and includes a diverse group of cancers. The heterogeneity of TNBC and the lack of target receptors on the cell surface make it difficult to develop specific therapeutic treatments. These aspects cause the high negative prognosis of patients with this type of tumor. The analysis of the molecular profiles of TNBC samples has allowed a better characterization of this tumor, supporting the search for new reliable diagnostic markers. To this end, we have developed a bioinformatic approach to integrate networks of genes differentially expressed in basal breast cancer compared to healthy tissues, with miRNAs able to regulate their expression. We studied the role of these miRNAs in TNBC subtype cell lines. We therefore identified two miRNAs, namely miR-135b and miR-365, with a central role in regulating the altered functional pathways in basal breast cancer. These two miRNAs are differentially expressed in human TNBC immunohistochemistry-selected tissues, and their modulation has been shown to play a role in the proliferation of tumor control and its migratory and invasive capacity in TNBC subtype cell lines. From the perspective of personalized medicine, we managed to modulate the expression of the two miRNAs in organotypic cultures, suggesting their possible use as diagnostic and therapeutic molecules. miR-135b and miR-365 have a key role in TNBC, controlling proliferation and invasion. Their detection could be helpful in TNBC diagnosis, while their modulation could become a new therapeutic tool for TNBC.

Patient-Specific Network for Personalized Breast Cancer Therapy with Multi-Omics Data

The development of new computational approaches that are able to design the correct personalized drugs is the crucial therapeutic issue in cancer research. However, tumor heterogeneity is the main obstacle to developing patient-specific single drugs or combinations of drugs that already exist in clinics. In this study, we developed a computational approach that integrates copy number alteration, gene expression, and a protein interaction network of 73 basal breast cancer samples. 2509 prognostic genes harboring a copy number alteration were identified using survival analysis, and a protein–protein interaction network considering the direct interactions was created. Each patient was described by a specific combination of seven altered hub proteins that fully characterize the 73 basal breast cancer patients. We suggested the optimal combination therapy for each patient considering drug–protein interactions. Our approach is able to confirm well-known cancer related genes and suggest novel potential drug target genes. In conclusion, we presented a new computational approach in breast cancer to deal with the intra-tumor heterogeneity towards personalized cancer therapy.