Search

Abstract Triple negative breast cancer (TNBC) accounts for 15–20% of breast cancer cases in the United States. Systemic neoadjuvant chemotherapy (NACT), with or without immunotherapy, is the current standard of care for patients with early-stage TNBC. However, up to 70% of TNBC patients have significant residual disease once NACT is completed, which is associated with a high risk of developing recurrence within two to three years of surgical resection. To identify targetable vulnerabilities in chemoresistant TNBC, we generated longitudinal patient-derived xenograft (PDX) models from TNBC tumors before and after patients received NACT. We then compiled transcriptomes and drug response profiles for all models. Transcriptomic analysis identified the enrichment of aberrant protein homeostasis pathways in models from post-NACT tumors relative to pre-NACT tumors. This observation correlated with increased sensitivity in vitro to inhibitors targeting the proteasome, heat shock proteins, and neddylation pathways. Pevonedistat, a drug annotated as a NEDD8-activating enzyme (NAE) inhibitor, was prioritized for validation in vivo and demonstrated efficacy as a single agent in multiple PDX models of TNBC. Pharmacotranscriptomic analysis identified a pathway-level correlation between pevonedistat activity and post-translational modification (PTM) machinery, particularly involving neddylation and sumoylation targets. Elevated levels of both NEDD8 and SUMO1 were observed in models exhibiting a favorable response to pevonedistat compared to those with a less favorable response in vivo. Moreover, a correlation emerged between the expression of neddylation-regulated pathways and tumor response to pevonedistat, indicating that targeting these PTM pathways may prove effective in combating chemoresistant TNBC.

Abstract Tumors are comprised of subpopulations of cancer cells that harbor distinct genetic profiles and phenotypes that evolve over time and during treatment. By reconstructing the course of cancer evolution, we can understand the acquisition of the malignant properties that drive tumor progression. Unfortunately, recovering the evolutionary relationships of individual cancer cells linked to their phenotypes remains a difficult challenge. To address this need, we have developed PhylinSic, a method that reconstructs the phylogenetic relationships among cells linked to their gene expression profiles from single cell RNA-sequencing (scRNA-Seq) data. This method calls nucleotide bases using a probabilistic smoothing approach and then estimates a phylogenetic tree using a Bayesian modeling algorithm. We showed that PhylinSic identified evolutionary relationships underpinning drug selection and metastasis and was sensitive enough to identify subclones from genetic drift. We found that breast cancer tumors resistant to chemotherapies harbored multiple genetic lineages that independently acquired high K-Ras and β-catenin, suggesting that therapeutic strategies may need to control multiple lineages to be durable. These results demonstrated that PhylinSic can reconstruct evolution and link the genotypes and phenotypes of cells across monophyletic tumors using scRNA-Seq.

Abstract Publications comparing variant caller algorithms present discordant results with contradictory rankings. Caller performances are inconsistent and wide ranging, and dependent upon input data, application, parameter settings, and evaluation metric. With no single variant caller emerging as a superior standard, combinations or ensembles of variant callers have appeared in the literature. In this study, a whole genome somatic reference standard was used to derive principles to guide strategies for combining variant calls. Then, manually annotated variants called from the whole exome sequencing of a tumor were used to corroborate these general principles. Finally, we examined the ability of these principles to reduce noise in targeted sequencing.

KRAS is a small GTPase that regulates cell proliferation. Here, the authors show that a subset of cell surface glycosphingolipids regulate KRAS plasma membrane localization by modulating inner leaflet lipid composition, uncovering a requirement for KRAS oncogenesis that may have therapeutic potential.

Abstract Patient-derived xenograft (PDX) models of breast cancer are an effective discovery platform and tool for preclinical pharmacologic testing and biomarker identification. We established orthotopic PDX models of triple negative breast cancer (TNBC) from the primary breast tumors of patients prior to and following neoadjuvant chemotherapy (NACT) while they were enrolled in the ARTEMIS trial (NCT02276443). Serial biopsies were obtained from patients prior to treatment (pre-NACT), from poorly responsive disease after four cycles of Adriamycin and cyclophosphamide (AC, mid-NACT), and in cases of AC-resistance, after a 3-month course of different experimental therapies and/or additional chemotherapy (post-NACT). Our study cohort includes a total of 269 fine needle aspirates (FNAs) from 217 women, generating a total of 62 PDX models (overall success-rate = 23%). Success of PDX engraftment was generally higher from those cancers that proved to be treatment-resistant, whether poorly responsive to AC as determined by ultrasound measurements mid-NACT (p = 0.063), RCB II/III status after NACT (p = 0.046), or metastatic relapse within 2 years of surgery (p = 0.008). TNBC molecular subtype determined from gene expression microarrays of pre-NACT tumors revealed no significant association with PDX engraftment rate (p = 0.877). Finally, we developed a statistical model predictive of PDX engraftment using percent Ki67 positive cells in the patient’s diagnostic biopsy, positive lymph node status at diagnosis, and low volumetric reduction of the patient’s tumor following AC treatment. This novel bank of 62 PDX models of TNBC provides a valuable resource for biomarker discovery and preclinical therapeutic trials aimed at improving neoadjuvant response rates for patients with TNBC.

Background: Recent advances have been made in targeting the phosphoinositide 3-kinase pathway in breast cancer. Phosphatase and tensin homolog (PTEN) is a key component of that pathway. Objective: To understand the changes in PTEN expression over the course of the disease in patients with triple-negative breast cancer (TNBC) and whether PTEN copy number variation (CNV) by next-generation sequencing (NGS) can serve as an alternative to immunohistochemistry (IHC) to identify PTEN loss. Methods: We compared PTEN expression by IHC between pretreatment tumors and residual tumors in the breast and lymph nodes after neoadjuvant chemotherapy in 96 patients enrolled in a TNBC clinical trial. A correlative analysis between PTEN protein expression and PTEN CNV by NGS was also performed. Results: With a stringent cutoff for PTEN IHC scoring, PTEN expression was discordant between pretreatment and posttreatment primary tumors in 5% of patients ( n = 96) and between posttreatment primary tumors and lymph node metastases in 9% ( n = 33). A less stringent cutoff yielded similar discordance rates. Intratumoral heterogeneity for PTEN loss was observed in 7% of the patients. Among pretreatment tumors, PTEN copy numbers by whole exome sequencing ( n = 72) were significantly higher in the PTEN-positive tumors by IHC compared with the IHC PTEN-loss tumors ( p

Abstract Germline mutations in BRCA1 or BRCA2 exist in ~2–7% of breast cancer patients, which has led to the approval of PARP inhibitors in the advanced setting. We have previously reported a phase II neoadjuvant trial of single agent talazoparib for patients with germline BRCA pathogenic variants with a pathologic complete response (pCR) rate of 53%. As nearly half of the patients treated did not have pCR, better strategies are needed to overcome treatment resistance. To this end, we conducted multi-omic analysis of 13 treatment naïve breast cancer tumors from patients that went on to receive single-agent neoadjuvant talazoparib. We looked for biomarkers that were predictive of response (assessed by residual cancer burden) after 6 months of therapy. We found that all resistant tumors exhibited either the loss of SHLD2, expression of a hypoxia signature, or expression of a stem cell signature. These results indicate that the deep analysis of pre-treatment tumors can identify biomarkers that are predictive of response to talazoparib and potentially other PARP inhibitors, and provides a framework that will allow for better selection of patients for treatment, as well as a roadmap for the development of novel combination therapies to prevent emergence of resistance.

Epithelial-mesenchymal transition (EMT) contributes to tumor invasion, metastasis and drug resistance. AKT activation is key in a number of cellular processes. While many positive regulators for either EMT or AKT activation have been reported, few negative regulators are established. Through kinase cDNA screen, we identified brain-type creatine kinase (CKB or BCK) as a potent suppressor for both. As a ubiquitously expressed kinase in normal tissues, CKB is significantly downregulated in several solid cancer types. Lower CKB expression is significantly associated with worse prognosis. Phenotypically, CKB overexpression suppresses, while its silencing promotes, EMT and cell migration, xenograft tumor growth and metastasis of prostate cancer cells. AKT activation is one of the most prominent signaling events upon CKB silencing in prostate cancer cells, which is in line with prostate cancer TCGA data. EMT enhanced by CKB silencing is abolished by AKT inhibition. Mechanistically, CKB interacts with AKT and sequestrates it from activation by mTOR. We further elucidated that an 84aa fragment at C-terminus of CKB protein interacts with AKT's PH domain. Ectopic expression of the 84aa CKB fragment inhibits AKT activation, EMT and cell proliferation. Interestingly, molecular dynamics simulation on crystal structures of AKT and CKB independently demonstrates that AKT's PH domain and CKB's 84aa fragment establish their major interaction interface. In summary, we have discovered CKB as a negative regulator of EMT and AKT activation, revealing a new mode of their regulation . We have also demonstrated that CKB downregulation is a poor prognosticator, which is sufficient to promote prostate cancer progression.

Human metaplastic breast cancers (MpBC) are a rare, aggressive subclass of triple-negative breast cancers. Here, the authors show over-expression of histone reader TRIM24 is sufficient to generate tumors with a molecular signature of metabolic dysfunction and EMT in a mouse model of human MpBC.

Endocrine therapy remains the primary treatment choice for ER+ breast cancers. However, most advanced ER+ breast cancers ultimately develop resistance to endocrine. This acquired resistance to endocrine therapy is often driven by the activation of the PI3K/AKT/mTOR signaling pathway. Everolimus, a drug that targets and inhibits the mTOR complex has been shown to improve clinical outcomes in metastatic ER+ breast cancers. However, there are no biomarkers currently available to guide the use of everolimus in the clinic for progressive patients, where multiple therapeutic options are available. Here, we utilized gene expression signatures from 9 ER+ breast cancer cell lines and 23 patients treated with everolimus to develop and validate an integrative machine learning biomarker of mTOR inhibitor response. Our results show that the machine learning biomarker can successfully distinguish responders from non-responders and can be applied to identify patients that will most likely benefit from everolimus treatment.

High-grade serous ovarian cancer (HGSOC) is prone to developing resistance to treatment. Here, the authors use single-cell RNA-seq and an analysis of archetypes, and find that shifts in metabolism and proliferation are associated with the response to treatment and clonal heterogeneity in HGSOC.

Effective cancer chemotherapy treatment requires both therapy delivery and retention by malignant cells. Cancer cell overexpression of the multidrug transmembrane transporter gene ABCB1 (MDR1, multi-drug resistance protein 1) thwarts therapy retention, leading to a drug-resistant phenotype. We explored the phenotypic impact of ABCB1 overexpression in normal human mammary epithelial cells (HMECs) via acute adenoviral delivery and in breast cancer cell lines with stable integration of inducible ABCB1 expression. One hundred sixty-two genes were differentially expressed between ABCB1-expressing and GFP-expressing HMECs, including the gene encoding the cyclooxygenase-2 protein, PTGS2. Several breast cancer cell lines with inducible ABCB1 expression demonstrated sensitivity to the 5-lipoxygenase, cyclooxygenase-1/2 inhibitor tepoxalin in two-dimensional drug response assays, and combination treatment of tepoxalin either with chemotherapies or with histone deacetylase (HDAC) inhibitors improved therapeutic efficacy in these lines. Moreover, selection for the ABCB1-expressing cell population was reduced in three-dimensional co-cultures of ABCB1-expressing and GFP-expressing cells when chemotherapy was given in combination with tepoxalin. Further study is warranted to ascertain the clinical potential of tepoxalin, an FDA-approved therapeutic for use in domesticated mammals, to restore chemosensitivity and improve drug response in patients with ABCB1-overexpressing drug-resistant breast cancers.

Abstract Background Triple-negative breast cancers (TNBCs), which lack receptors for estrogen, progesterone, and amplification of epidermal growth factor receptor 2, are highly aggressive. Consequently, patients diagnosed with TNBCs have reduced overall and disease-free survival rates compared to patients with other subtypes of breast cancer. TNBCs are characterized by the presence of cancer cells with mesenchymal properties, indicating that the epithelial to mesenchymal transition (EMT) plays a major role in the progression of this disease. The EMT program has also been implicated in chemoresistance, tumor recurrence, and induction of cancer stem cell (CSC) properties. Currently, there are no targeted therapies for TNBC, and hence, it is critical to identify the novel targets to treat TNBC. Methods A library of compounds was screened for their ability to inhibit EMT in cells with mesenchymal phenotype as assessed using the previously described Z-cad reporters. Of the several drugs tested, GSK3β inhibitors were identified as EMT inhibitors. The effects of GSK3β inhibitors on the properties of TNBC cells with a mesenchymal phenotype were assessed using qRT-PCR, flow cytometry, western blot, mammosphere, and migration and cell viability assays. Publicly available datasets also were analyzed to examine if the expression of GSK3β correlates with the overall survival of breast cancer patients. Results We identified a GSK3β inhibitor, BIO, in a drug screen as one of the most potent inhibitors of EMT. BIO and two other GSK3β inhibitors, TWS119 and LiCl, also decreased the expression of mesenchymal markers in several different cell lines with a mesenchymal phenotype. Further, inhibition of GSK3β reduced EMT-related migratory properties of cells with mesenchymal properties. To determine if GSK3β inhibitors target mesenchymal-like cells by affecting the CSC population, we employed mammosphere assays and profiled the stem cell-related cell surface marker CD44+/24− in cells after exposure to GSK3β inhibitors. We found that GSK3β inhibitors indeed decreased the CSC properties of cell types with mesenchymal properties. We treated cells with epithelial and mesenchymal properties with GSK3β inhibitors and found that GSK3β inhibitors selectively kill cells with mesenchymal attributes while sparing cells with epithelial properties. We analyzed patient data to identify genes predictive of poor clinical outcome that could serve as novel therapeutic targets for TNBC. The Wnt signaling pathway is critical to EMT, but among the various factors known to be involved in Wnt signaling, only the higher expression of GSK3β correlated with poorer overall patient survival. Conclusions Taken together, our data demonstrate that GSK3β is a potential target for TNBCs and suggest that GSK3β inhibitors could serve as selective inhibitors of EMT and CSC properties for the treatment of a subset of aggressive TNBC. GSK3β inhibitors should be tested for use in combination with standard-of-care drugs in preclinical TNBC models.

Objective Polycystic kidney disease (PKD) is the major cause of kidney failure and mortality in humans. It has always been suspected that the development of cystic kidney disease shares features with tumorigenesis, although the evidence is unclear. Methods We crossed p53 mutant mice (p53N236S, p53S) with Werner syndrome mice and analyzed the pathological phenotypes. The RNA-seq, ssGSEA analysis, and real-time PCR were performed to dissect the gene signatures involved in the development of disease phenotypes.Results We found enlarged kidneys with fluid-filled cysts in offspring mice with a genotype of G3mTerc-/-WRN-/-p53S/S (G3TM). Pathology analysis confirmed the occurrence of PKD, and it was highly correlated with the incidence of tumorigenesis. RNA-seq data revealed the gene signatures involved in PKD development, and demonstrated that PKD and tumorigenesis shared common pathways, including complement pathways, lipid metabolism, mitochondria energy homeostasis and others. Interestingly, this G3TM PKD and the classical PKD1/2 deficient PKD shared common pathways, possibly because the mutant p53S could regulate the expression levels of PKD1/2, Pkhd1, and Hnf1b.Conclusions We established a dual mouse model for PKD and tumorigenesis derived from abnormal cellular proliferation and telomere dysfunction. The innovative point of our study is to report PKD occurring in conjunction with tumorigenesis. The gene signatures revealed might shed new light on the pathogenesis of PKD, and provide new molecular biomarkers for clinical diagnosis and prognosis.

Cancer cell phenotypes evolve during a tumor's treatment. In some cases, tumor cells acquire cancer stem cell-like (CSL) traits such as resistance to chemotherapy and diminished differentiation; therefore, targeting these cells may be therapeutically beneficial. In this study we show that in progressive estrogen receptor positive (ER+) metastatic breast cancer tumors, resistant subclones that emerge following chemotherapy have increased CSL abundance. Further, in vitro organoid growth of ER+ patient cancer cells also shows that chemotherapy treatment leads to increased abundance of ALDH+/CD44+ CSL cells. Chemotherapy induced CSL abundance is blocked by treatment with a pan-HDAC inhibitor, belinostat. Belinostat treatment diminished both mammosphere formation and size following chemotherapy, indicating a decrease in progenitor CSL traits. HDAC inhibitors specific to class IIa (HDAC4, HDAC5) and IIb (HDAC6) were shown to primarily reverse the chemo-resistant CSL state. Single-cell RNA sequencing analysis with patient samples showed that HDAC targets and MYC signaling were promoted by chemotherapy and inhibited upon HDAC inhibitor treatment. In summary, HDAC inhibition can block chemotherapy-induced drug resistant phenotypes with ‘one-two punch’ strategy in refractory breast cancer cells.

The ATP-Binding Cassette transporter A1 (ABCA1) reverse cholesterol transport channel has been associated with a number of phenotypes in breast cancer, including reduced proliferation and increased metastatic capacity. It is induced in an epithelial–mesenchymal transition (EMT), but little is known about how this occurs, and whether it is sufficient to promote metastatic phenotypes. To address these questions, we have deciphered the transcriptional regulation of ABCA1 across EMT states and found that it is repressed by MYC via an E-box element in its P1 alternative promoter. De-repression of the promoter by MYC knockdown leads to induction of ABCA1 expression. This indicates that ABCA1 expression is regulated in an EMT, revealing another link between ABCA1 and malignant phenotypes.

It is unclear how intra-tumoral heterogeneity contributes to metastasis. Here the authors study the clonal dynamics of triple negative breast cancer metastasis using patient derived xenografts and demonstrate that primary tumor clones harbor properties that support seeding and colonization of multiple organs.

Abstract Background The RhoA activating protein Net1 contributes to breast cancer cell proliferation, motility, and invasion in vitro, yet little is known about its roles in mammary gland tumorigenesis and metastasis. Methods Net1 knockout (KO) mice were bred to mice with mammary gland specific expression of the polyoma middle T antigen (PyMT) oncogene. Mammary gland tumorigenesis and lung metastasis were monitored. Individual tumors were assessed for proliferation, apoptosis, angiogenesis, RhoA activation, and activation of PyMT-dependent signaling pathways. Primary tumor cells from wild-type and Net1 KO mice were transplanted into the mammary glands of wild-type, nontumor-bearing mice, and tumor growth and metastasis were assessed. Gene expression in wild-type and Net1 KO tumors was analyzed by gene ontology enrichment and for relative activation of gene expression signatures indicative of signaling pathways important for breast cancer initiation and progression. A gene expression signature indicative of Net1 function was identified. Human breast cancer gene expression profiles were screened for the presence of a Net1 gene expression signature. Results We show that Net1 makes fundamental contributions to mammary gland tumorigenesis and metastasis. Net1 deletion delays tumorigenesis and strongly suppresses metastasis in PyMT-expressing mice. Moreover, we observe that loss of Net1 reduces cancer cell proliferation, inhibits tumor angiogenesis, and promotes tumor cell apoptosis. Net1 is required for maximal RhoA activation within tumors and for primary tumor cell motility. Furthermore, the ability of PyMT to initiate oncogenic signaling to ERK1/2 and PI3K/Akt1 is inhibited by Net1 deletion. Primary tumor cell transplantation indicates that the reduction in tumor angiogenesis and lung metastasis observed upon Net1 deletion are tumor cell autonomous effects. Using a gene expression signature indicative of Net1 activity, we show that Net1 signaling is activated in 10% of human breast cancers, and that this correlates with elevated proliferation and PI3K pathway activity. We also demonstrate that human breast cancer patients with a high Net1 gene expression signature experience shorter distant metastasis-free survival. Conclusions These data indicate that Net1 is required for tumor progression in the PyMT mouse model and suggest that Net1 may contribute to breast cancer progression in humans.

Proprotein convertase subtilisin/kexin type 9 (PCSK9) interacts directly with cytoplasmic apoB and prevents its degradation via the autophagosome/lysosome pathway. This process affects VLDL and LDL production and influences atherogenesis. Here, we investigated the molecular machinery by which PCSK9 modulates autophagy and affects atherogenesis. We backcrossed Pcsk9−/− mice with atherosclerosis-prone Ldlr−/−Apobec1−/− (LDb) mice to generate Ldlr−/−Apobec1−/−Pcsk9−/− (LTp) mice. Deletion of PCSK9 resulted in decreased hepatic apoB secretion, increased autophagic flux, and decreased plasma levels of IDL and LDL particles. The LDLs from LTp mice (LTp-LDLs) were less atherogenic and contained less cholesteryl ester and phospholipids than LDb-LDLs. Moreover LTp-LDLs induced lower endothelial expression of the genes encoding TLR2, Lox-1, ICAM-1, CCL2, CCL7, IL-6, IL-1β, Beclin-1, p62, and TRAF6. Collectively, these effects were associated with substantially less atherosclerosis development (>4-fold) in LTp mice. The absence of PCSK9 in LDb mice results in decreased lipid and apoB levels, fewer atherogenic LDLs, and marked reduction of atherosclerosis. The effect on atherogenesis may be mediated in part by the effects of modified LDLs on endothelial cell receptors and proinflammatory and autophagy molecules. These findings suggest that there may be clinical benefits of PCSK9 inhibition due to mechanisms unrelated to increased LDL receptor activity.

Abstract Epithelial ovarian cancer (EOC) is hallmarked by a high degree of heterogeneity. To address this heterogeneity, a classification scheme was developed based on gene expression patterns of 1538 tumours. Five, biologically distinct subgroups — Epi‐A, Epi‐B, Mes, Stem‐A and Stem‐B — exhibited significantly distinct clinicopathological characteristics, deregulated pathways and patient prognoses, and were validated using independent datasets. To identify subtype‐specific molecular targets, ovarian cancer cell lines representing these molecular subtypes were screened against a genome‐wide shRNA library. Focusing on the poor‐prognosis Stem‐A subtype, we found that two genes involved in tubulin processing, TUBGCP4 and NAT10, were essential for cell growth, an observation supported by a pathway analysis that also predicted involvement of microtubule‐related processes. Furthermore, we observed that Stem‐A cell lines were indeed more sensitive to inhibitors of tubulin polymerization, vincristine and vinorelbine, than the other subtypes. This subtyping offers new insights into the development of novel diagnostic and personalized treatment for EOC patients.

Command-line software plays a critical role in biology research. However, processes for installing and executing software differ widely. The Common Workflow Language (CWL) is a community standard that addresses this problem. Using CWL, tool developers can formally describe a tool’s inputs, outputs, and other execution details. CWL documents can include instructions for executing tools inside software containers. Accordingly, CWL tools are portable—they can be executed on diverse computers—including personal workstations, high-performance clusters, or the cloud. CWL also supports workflows, which describe dependencies among tools and using outputs from one tool as inputs to others. To date, CWL has been used primarily for batch processing of large datasets, especially in genomics. But it can also be used for analytical steps of a study. This article explains key concepts about CWL and software containers and provides examples for using CWL in biology research. CWL documents are text-based, so they can be created manually, without computer programming. However, ensuring that these documents conform to the CWL specification may prevent some users from adopting it. To address this gap, we created ToolJig, a Web application that enables researchers to create CWL documents interactively. ToolJig validates information provided by the user to ensure it is complete and valid. After creating a CWL tool or workflow, the user can create ‘input-object’ files, which store values for a particular invocation of a tool or workflow. In addition, ToolJig provides examples of how to execute the tool or workflow via a workflow engine. ToolJig and our examples are available at https://github.com/srp33/ToolJig.