E-cadherin expression and gene expression profiles in corticotroph pituitary neuroendocrine tumor subtypes

Corticotroph adenomas/pituitary neuroendocrine tumors (PitNETs) are associated with significant morbidity and mortality. Predictors of tumor behavior have not shown high prognostic accuracy. For somatotroph adenomas/PitNETs, E-cadherin expression correlates strongly with prognosis. E-cadherin expression has not been investigated in other PitNETs. A retrospective chart review of adults with corticotroph adenomas/PitNETs was conducted to assess correlation between E-cadherin expression and tumor characteristics. In addition, gene expression microarray was performed in subset of tumors (n = 16). Seventy-seven patients were identified; 71% were female, with median age of cohort 45.2 years. Seventy-five percent had macroadenomas, of which 22% were hormonally active. Ninety-five percent of microadenomas were hormonally active. Adrenocorticotropic hormone granulation pattern by IHC identified 63% as densely granulated (DG) and 34% as sparsely granulated (SG). All microadenomas were DG (p < .001); 50% of macroadenomas were DG associated with increased tumor invasion compared to SG. E-cadherin IHC was positive in 80%, diminished in 17%, and absent in 20% and did not correlate with corticotroph PitNETs subtype, size, or prognosis. In contrast to the distinct transcriptomes of corticotroph PitNETs and normal pituitaries, a comparison of clinically active and silent corticotroph PitNETs demonstrated similar molecular signatures indicating their common origin, but with unique differences related to their secretory status.

NPEPPS Is a Druggable Driver of Platinum Resistance

There is an unmet need to improve the efficacy of platinum-based cancer chemotherapy, which is used in primary and metastatic settings in many cancer types. In bladder cancer, platinum-based chemotherapy leads to better outcomes in a subset of patients when used in the neoadjuvant setting or in combination with immunotherapy for advanced disease. Despite such promising results, extending the benefits of platinum drugs to a greater number of patients is highly desirable. Using the multiomic assessment of cisplatin-responsive and -resistant human bladder cancer cell lines and whole-genome CRISPR screens, we identified puromycin-sensitive aminopeptidase (NPEPPS) as a driver of cisplatin resistance. NPEPPS depletion sensitized resistant bladder cancer cells to cisplatin in vitro and in vivo. Conversely, overexpression of NPEPPS in sensitive cells increased cisplatin resistance. NPEPPS affected treatment response by regulating intracellular cisplatin concentrations. Patient-derived organoids (PDO) generated from bladder cancer samples before and after cisplatin-based treatment, and from patients who did not receive cisplatin, were evaluated for sensitivity to cisplatin, which was concordant with clinical response. In the PDOs, depletion or pharmacologic inhibition of NPEPPS increased cisplatin sensitivity, while NPEPPS overexpression conferred resistance. Our data present NPEPPS as a druggable driver of cisplatin resistance by regulating intracellular cisplatin concentrations.

SIGNIFICANCE

Targeting NPEPPS, which induces cisplatin resistance by controlling intracellular drug concentrations, is a potential strategy to improve patient responses to platinum-based therapies and lower treatment-associated toxicities.

ProstaMine: a bioinformatics tool for identifying subtype-specific co-alterations associated with aggressiveness in prostate cancer

Background

Prostate cancer is a leading cause of cancer-related deaths among men, marked by heterogeneous clinical and molecular characteristics. The complexity of the molecular landscape necessitates tools for identifying multi-gene co-alteration patterns that are associated with aggressive disease. The identification of such gene sets will allow for deeper characterization of the processes underlying prostate cancer progression and potentially lead to novel strategies for treatment.

Methods

We developed ProstaMine to systematically identify co-alterations associated with aggressiveness in prostate cancer molecular subtypes defined by high-fidelity alterations in primary prostate cancer. ProstaMine integrates genomic, transcriptomic, and clinical data from five primary and one metastatic prostate cancer cohorts to prioritize co-alterations enriched in metastatic disease and associated with disease progression.

Results

Integrated analysis of primary tumors defined a set of 17 prostate cancer alterations associated with aggressive characteristics. We applied ProstaMine to NKX3-1-loss and RB1-loss tumors and identified subtype-specific co-alterations associated with metastasis and biochemical relapse in these molecular subtypes. In NKX3-1-loss prostate cancer, ProstaMine identified novel subtype-specific co-alterations known to regulate prostate cancer signaling pathways including MAPK, NF-kB, p53, PI3K, and Sonic hedgehog. In RB1-loss prostate cancer, ProstaMine identified novel subtype-specific co-alterations involved in p53, STAT6, and MHC class I antigen presentation. Co-alterations impacting autophagy were noted in both molecular subtypes.

Conclusion

ProstaMine is a method to systematically identify novel subtype-specific co-alterations associated with aggressive characteristics in prostate cancer. The results from ProstaMine provide insights into potential subtype-specific mechanisms of prostate cancer progression which can be formed into testable experimental hypotheses. ProstaMine is publicly available at: https://bioinformatics.cuanschutz.edu/prostamine.

Microbiome preterm birth DREAM challenge: Crowdsourcing machine learning approaches to advance preterm birth research

Every year, 11% of infants are born preterm with significant health consequences, with the vaginal microbiome a risk factor for preterm birth. We crowdsource models to predict (1) preterm birth (PTB; <37 weeks) or (2) early preterm birth (ePTB; <32 weeks) from 9 vaginal microbiome studies representing 3,578 samples from 1,268 pregnant individuals, aggregated from public raw data via phylogenetic harmonization. The predictive models are validated on two independent unpublished datasets representing 331 samples from 148 pregnant individuals. The top-performing models (among 148 and 121 submissions from 318 teams) achieve area under the receiver operator characteristic (AUROC) curve scores of 0.69 and 0.87 predicting PTB and ePTB, respectively. Alpha diversity, VALENCIA community state types, and composition are important features in the top-performing models, most of which are tree-based methods. This work is a model for translation of microbiome data into clinically relevant predictive models and to better understand preterm birth.

Down syndrome is associated with altered frequency and functioning of tracheal multiciliated cells, and response to influenza virus infection

Individuals with Down syndrome (DS) clinically manifest severe respiratory illnesses; however, there is a paucity of data on how DS influences homeostatic physiology of lung airway, and its reactive responses to pulmonary pathogens. We generated well-differentiated ciliated airway epithelia using tracheas from wild-type and Dp(16)1/Yey mice in vitro, and discovered that Dp(16)1/Yey epithelia have significantly lower abundance of ciliated cells, an altered ciliary beating profile, and reduced mucociliary transport. Interestingly, both sets of differentiated epithelia released similar quantities of viral particles after infection with influenza A virus (IAV). However, RNA-sequencing and proteomic analyses revealed an immune hyperreactive phenotype particularly for monocyte-recruiting chemokines in Dp(16)1/Yey epithelia. Importantly, when we challenged mice in vivo with IAV, we observed immune hyper-responsiveness in Dp(16)1/Yey mice, evidenced by higher quantities of lung airway infiltrated monocytes, and elevated levels of pro-inflammatory cytokines in bronchoalveolar lavage fluid. Our findings illuminate mechanisms underlying DS-mediated pathophysiological changes in airway epithelium.

SIX1 and EWS/FLI1 co-regulate an anti-metastatic gene network in Ewing Sarcoma

Ewing sarcoma (ES), which is characterized by the presence of oncogenic fusion proteins such as EWS/FLI1, is an aggressive pediatric malignancy with a high rate of early dissemination and poor outcome after distant spread. Here we demonstrate that the SIX1 homeoprotein, which enhances metastasis in most tumor types, suppresses ES metastasis by co-regulating EWS/FLI1 target genes. Like EWS/FLI1, SIX1 promotes cell growth/transformation, yet dramatically inhibits migration and invasion, as well as metastasis in vivo. We show that EWS/FLI1 promotes SIX1 protein expression, and that the two proteins share genome-wide binding profiles and transcriptional regulatory targets, including many metastasis-associated genes such as integrins, which they co-regulate. We further show that SIX1 downregulation of integrins is critical to its ability to inhibit invasion, a key characteristic of metastatic cells. These data demonstrate an unexpected anti-metastatic function for SIX1, through coordinate gene regulation with the key oncoprotein in ES, EWS/FLI1.

A harmonized resource of integrated prostate cancer clinical, -omic, and signature features

Genomic and transcriptomic data have been generated across a wide range of prostate cancer (PCa) study cohorts. These data can be used to better characterize the molecular features associated with clinical outcomes and to test hypotheses across multiple, independent patient cohorts. In addition, derived features, such as estimates of cell composition, risk scores, and androgen receptor (AR) scores, can be used to develop novel hypotheses leveraging existing multi-omic datasets. The full potential of such data is yet to be realized as independent datasets exist in different repositories, have been processed using different pipelines, and derived and clinical features are often not provided or  not standardized. Here, we present the curatedPCaData R package, a harmonized data resource representing >2900 primary tumor, >200 normal tissue, and >500 metastatic PCa samples across 19 datasets processed using standardized pipelines with updated gene annotations. We show that meta-analysis across harmonized studies has great potential for robust and clinically meaningful insights. curatedPCaData is an open and accessible community resource with code made available for reproducibility.

Non-muscle-invasive bladder cancer molecular subtypes predict differential response to intravesical Bacillus Calmette-Guérin

The recommended treatment for patients with high-risk non-muscle-invasive bladder cancer (HR-NMIBC) is tumor resection followed by adjuvant Bacillus Calmette-Guérin (BCG) bladder instillations. However, only 50% of patients benefit from this therapy. If progression to advanced disease occurs, then patients must undergo a radical cystectomy with risks of substantial morbidity and poor clinical outcome. Identifying tumors unlikely to respond to BCG can translate into alternative treatments, such as early radical cystectomy, targeted therapies, or immunotherapies. Here, we conducted molecular profiling of 132 patients with BCG-naive HR-NMIBC and 44 patients with recurrences after BCG (34 matched), which uncovered three distinct BCG response subtypes (BRS1, 2 and BRS3). Patients with BRS3 tumors had a reduced recurrence-free and progression-free survival compared with BRS1/2. BRS3 tumors expressed high epithelial-to-mesenchymal transition and basal markers and had an immunosuppressive profile, which was confirmed with spatial proteomics. Tumors that recurred after BCG were enriched for BRS3. BRS stratification was validated in a second cohort of 151 BCG-naive patients with HR-NMIBC, and the molecular subtypes outperformed guideline-recommended risk stratification based on clinicopathological variables. For clinical application, we confirmed that a commercially approved assay was able to predict BRS3 tumors with an area under the curve of 0.87. These BCG response subtypes will allow for improved identification of patients with HR-NMIBC at the highest risk of progression and have the potential to be used to select more appropriate treatments for patients unlikely to respond to BCG.

Microparticle-Delivered Cxcl9 Prolongs Braf Inhibitor Efficacy in Melanoma

Patients with BRAF-mutant melanoma show substantial responses to combined BRAF and MEK inhibition, but most relapse within 2 years. A major reservoir for drug resistance is minimal residual disease (MRD), comprised of drug-tolerant tumor cells laying in a dormant state. Towards exploiting potential therapeutic vulnerabilities of MRD, we established a genetically engineered mouse model of BrafV600E-driven melanoma MRD wherein genetic BrafV600E extinction leads to strong but incomplete tumor regression. Transcriptional time-course analysis after BrafV600E extinction revealed that after an initial surge of immune activation, tumors later became immunologically "cold" after MRD establishment. Computational analysis identified candidate T-cell recruiting chemokines as strongly upregulated initially and steeply decreasing as the immune response faded. Therefore, we hypothesized that sustaining chemokine signaling could impair MRD maintenance through increased recruitment of effector T cells. We found that intratumoral administration of recombinant Cxcl9 (rCxcl9), either naked or loaded in microparticles, significantly impaired MRD relapse in BRAF-inhibited tumors, including several complete pathologic responses after microparticle-delivered rCxcl9 combined with BRAF and MEK inhibition. Our experiments constitute proof of concept that chemokine-based microparticle delivery systems are a potential strategy to forestall tumor relapse and thus improve the clinical success of first-line treatment methods.

Microbiome Preterm Birth DREAM Challenge: Crowdsourcing Machine Learning Approaches to Advance Preterm Birth Research

Globally, every year about 11% of infants are born preterm, defined as a birth prior to 37 weeks of gestation, with significant and lingering health consequences. Multiple studies have related the vaginal microbiome to preterm birth. We present a crowdsourcing approach to predict: (a) preterm or (b) early preterm birth from 9 publicly available vaginal microbiome studies representing 3,578 samples from 1,268 pregnant individuals, aggregated from raw sequences via an open-source tool, MaLiAmPi. We validated the crowdsourced models on novel datasets representing 331 samples from 148 pregnant individuals. From 318 DREAM challenge participants we received 148 and 121 submissions for our two separate prediction sub-challenges with top-ranking submissions achieving bootstrapped AUROC scores of 0.69 and 0.87, respectively. Alpha diversity, VALENCIA community state types, and composition (via phylotype relative abundance) were important features in the top performing models, most of which were tree based methods. This work serves as the foundation for subsequent efforts to translate predictive tests into clinical practice, and to better understand and prevent preterm birth.

curatedPCaData: Integration of clinical, genomic, and signature features in a curated and harmonized prostate cancer data resource

Genomic and transcriptomic data have been generated across a wide range of prostate cancer (PCa) study cohorts. These data can be used to better characterize the molecular features associated with clinical outcomes and to test hypotheses across multiple, independent patient cohorts. In addition, derived features, such as estimates of cell composition, risk scores, and androgen receptor (AR) scores, can be used to develop novel hypotheses leveraging existing multi-omic datasets. The full potential of such data is yet to be realized as independent datasets exist in different repositories, have been processed using different pipelines, and derived and clinical features are often not provided or unstandardized. Here, we present the curatedPCaData R package, a harmonized data resource representing >2900 primary tumor, >200 normal tissue, and >500 metastatic PCa samples across 19 datasets processed using standardized pipelines with updated gene annotations. We show that meta-analysis across harmonized studies has great potential for robust and clinically meaningful insights. curatedPCaData is an open and accessible community resource with code made available for reproducibility.

SOPHIE: Generative Neural Networks Separate Common and Specific Transcriptional Responses

Genome-wide transcriptome profiling identifies genes that are prone to differential expression (DE) across contexts, as well as genes with changes specific to the experimental manipulation. Distinguishing genes that are specifically changed in a context of interest from common differentially expressed genes (DEGs) allows more efficient prediction of which genes are specific to a given biological process under scrutiny. Currently, common DEGs or pathways can only be identified through the laborious manual curation of experiments, an inordinately time-consuming endeavor. Here we pioneer an approach, Specific cOntext Pattern Highlighting In Expression data (SOPHIE), for distinguishing between common and specific transcriptional patterns using a generative neural network to create a background set of experiments from which a null distribution of gene and pathway changes can be generated. We apply SOPHIE to diverse datasets including those from human, human cancer, and bacterial pathogen Pseudomonas aeruginosa. SOPHIE identifies common DEGs in concordance with previously described, manually and systematically determined common DEGs. Further molecular validation indicates that SOPHIE detects highly specific but low-magnitude biologically relevant transcriptional changes. SOPHIE's measure of specificity can complement log2 fold change values generated from traditional DE analyses. For example, by filtering the set of DEGs, one can identify genes that are specifically relevant to the experimental condition of interest. Consequently, these results can inform future research directions. All scripts used in these analyses are available at https://github.com/greenelab/generic-expression-patterns. Users can access https://github.com/greenelab/sophie to run SOPHIE on their own data.

Abstract IA020: Reprogramming of myogenic transcription factors in rhabdomyosarcoma

Rhabdomyosarcoma, a pediatric malignancy with partial resemblance to undifferentiated skeletal muscle, is characterized by high expression of myogenic-lineage transcription factors such as MYOD1 and MYOG. Despite high expression of these transcription factors, which in normal muscle result in differentiation, RMS cells fail to differentiate, suggesting the presence of factors that inhibit their normal differentiation-promoting functions. In this talk, I will present data that the key muscle transcriptional regulator, SIX1, which in development activates the myogenic regulatory factors (MRFs) and promotes muscle differentiation, in fact inhibits differentiation in fusion-negative (FN) RMS. SIX1 holds FN-RMS cells in a progenitor-like state by altering the chromatin landscape and causing MYOD1, a key MRF, to preferentially bind to regulatory regions of genes permissive to growth rather than differentiation. Loss of SIX1 results in re-localization of MYOD1 to promoters/enhancers of genes associated with differentiation, and further results in increased binding of MYOG at such loci. Altered binding of MYOD1 and MYOG in response to SIX1 loss results in marked inhibition of RMS growth in vivo, via induction of differentiation. These data suggest that SIX1 acts as a master regulatory factor controlling the fate of RMS cells. Data will be presented that suggest dynamic actions of SIX1 and its co-factors throughout normal muscle differentiation, whereby high levels of SIX1 expressed in early muscle differentiation may mimic a transcriptional state seen in RMS. We hypothesize that the specific levels of SIX1, combined with a unique combination of transcriptional co-factors, reprogram genome-wide binding of MRFs to different promoter/enhancer sites at specific developmental time points, and that RMS is trapped in an early developmental state where SIX1 represses differentiation via genome-wide alterations in MRF binding that favor growth. Understanding the co-factors that work with SIX1 to alter chromatin state and MRF binding may enable the discovery of novel targets whose inhibition could serve as a relatively non-toxic treatment to restore normal developmental processes and inhibit RMS progression.

Citation Format: Heide L. Ford, Jessica Y. Hsu, Etienne P. Danis, Stephanie Nance, Jenean H. O’Brien, Annika L. Gustafson, Veronica M. Wessells, Andrew E. Goodspeed, Jared C. Talbot, Sharon L. Amacher, Paul Jedlicka, Joshua C. Black, James C. Costello, Adam D. Durbin, Kristin B. Artinger. Reprogramming of myogenic transcription factors in rhabdomyosarcoma [abstract]. In: Proceedings of the AACR Special Conference: Sarcomas; 2022 May 9-12; Montreal, QC, Canada. Philadelphia (PA): AACR; Clin Cancer Res 2022;28(18_Suppl):Abstract nr IA020.

Antigens Expressed by Breast Cancer Cells Undergoing EMT Stimulate Cytotoxic CD8+ T Cell Immunity

Simple Summary

The transition of cells with epithelial characteristics to those with mesenchymal characteristics (termed EMT) facilitates breast cancer invasive capacity. The EMT program can also contribute to immunosuppressive and immunoevasive properties, altering susceptibility to immune cell recognition and killing. The goal of our study was to manipulate EMT to reveal potential neoantigens that might affect the ability of tumor cells to circumvent immune escape and/or be utilized as an anticancer vaccine to kill cancer cells exhibiting the cellular plasticity that permits therapy resistance and metastatic progression. We identified potential neoantigens resulting from EMT-associated altered gene expression and alternative splicing events and observed increased immunogenicity and susceptibility to killing of the more epithelial-like cancer cells. Although the tested peptides did not protect from tumor growth, a limited number of predicted neoantigens derived from intron retention events were tested. In the future, refined prediction programs may facilitate exciting antigen discoveries.

Abstract

Antigenic differences formed by alterations in gene expression and alternative splicing are predicted in breast cancer cells undergoing epithelial to mesenchymal transition (EMT) and the reverse plasticity known as MET. How these antigenic differences impact immune interactions and the degree to which they can be exploited to enhance immune responses against mesenchymal cells is not fully understood. We utilized a master microRNA regulator of EMT to alter mesenchymal-like EO771 mammary carcinoma cells to a more epithelial phenotype. A computational approach was used to identify neoantigens derived from the resultant differentially expressed somatic variants (SNV) and alternative splicing events (neojunctions). Using whole cell vaccines and peptide-based vaccines, we find superior cytotoxicity against the more-epithelial cells and explore the potential of neojunction-derived antigens to elicit T cell responses through experiments designed to validate the computationally predicted neoantigens. Overall, results identify EMT-associated splicing factors common to both mouse and human breast cancer cells as well as immunogenic SNV- and neojunction-derived neoantigens in mammary carcinoma cells.

A Crowdsourcing Approach to Develop Machine Learning Models to Quantify Radiographic Joint Damage in Rheumatoid Arthritis

Importance

An automated, accurate method is needed for unbiased assessment quantifying accrual of joint space narrowing and erosions on radiographic images of the hands and wrists, and feet for clinical trials, monitoring of joint damage over time, assisting rheumatologists with treatment decisions. Such a method has the potential to be directly integrated into electronic health records.

Objectives

To design and implement an international crowdsourcing competition to catalyze the development of machine learning methods to quantify radiographic damage in rheumatoid arthritis (RA).

Design, Setting, and Participants

This diagnostic/prognostic study describes the Rheumatoid Arthritis 2-Dialogue for Reverse Engineering Assessment and Methods (RA2-DREAM Challenge), which used existing radiographic images and expert-curated Sharp-van der Heijde (SvH) scores from 2 clinical studies (674 radiographic sets from 562 patients) for training (367 sets), leaderboard (119 sets), and final evaluation (188 sets). Challenge participants were tasked with developing methods to automatically quantify overall damage (subchallenge 1), joint space narrowing (subchallenge 2), and erosions (subchallenge 3). The challenge was finished on June 30, 2020.

Main Outcomes and Measures

Scores derived from submitted algorithms were compared with the expert-curated SvH scores, and a baseline model was created for benchmark comparison. Performances were ranked using weighted root mean square error (RMSE). The performance and reproductivity of each algorithm was assessed using Bayes factor from bootstrapped data, and further evaluated with a postchallenge independent validation data set.

Results

The RA2-DREAM Challenge received a total of 173 submissions from 26 participants or teams in 7 countries for the leaderboard round, and 13 submissions were included in the final evaluation. The weighted RMSEs metric showed that the winning algorithms produced scores that were very close to the expert-curated SvH scores. Top teams included Team Shirin for subchallenge 1 (weighted RMSE, 0.44), HYL-YFG (Hongyang Li and Yuanfang Guan) subchallenge 2 (weighted RMSE, 0.38), and Gold Therapy for subchallenge 3 (weighted RMSE, 0.43). Bootstrapping/Bayes factor approach and the postchallenge independent validation confirmed the reproducibility and the estimation concordance indices between final evaluation and postchallenge independent validation data set were 0.71 for subchallenge 1, 0.78 for subchallenge 2, and 0.82 for subchallenge 3.

Conclusions and Relevance

The RA2-DREAM Challenge resulted in the development of algorithms that provide feasible, quick, and accurate methods to quantify joint damage in RA. Ultimately, these methods could help research studies on RA joint damage and may be integrated into electronic health records to help clinicians serve patients better by providing timely, reliable, and quantitative information for making treatment decisions to prevent further damage.

PRDM paralogs antagonistically balance Wnt/β-catenin activity during craniofacial chondrocyte differentiation

Cranial neural crest cell (NCC)-derived chondrocyte precursors undergo a dynamic differentiation and maturation process to establish a scaffold for subsequent bone formation, alterations in which contribute to congenital birth defects. Here, we demonstrate that transcription factor and histone methyltransferase proteins Prdm3 and Prdm16 control the differentiation switch of cranial NCCs to craniofacial cartilage. Loss of either paralog results in hypoplastic and disorganized chondrocytes due to impaired cellular orientation and polarity. We show that these proteins regulate cartilage differentiation by controlling the timing of Wnt/β-catenin activity in strikingly different ways: Prdm3 represses whereas Prdm16 activates global gene expression, although both act by regulating Wnt enhanceosome activity and chromatin accessibility. Finally, we show that manipulating Wnt/β-catenin signaling pharmacologically or generating prdm3-/-;prdm16-/- double mutants rescues craniofacial cartilage defects. Our findings reveal upstream regulatory roles for Prdm3 and Prdm16 in cranial NCCs to control Wnt/β-catenin transcriptional activity during chondrocyte differentiation to ensure proper development of the craniofacial skeleton.

SIX1 reprograms myogenic transcription factors to maintain the rhabdomyosarcoma undifferentiated state

Rhabdomyosarcoma (RMS) is a pediatric muscle sarcoma characterized by expression of the myogenic lineage transcription factors (TFs) MYOD1 and MYOG. Despite high expression of these TFs, RMS cells fail to terminally differentiate, suggesting the presence of factors that alter their functions. Here, we demonstrate that the developmental TF SIX1 is highly expressed in RMS and critical for maintaining a muscle progenitor-like state. SIX1 loss induces differentiation of RMS cells into myotube-like cells and impedes tumor growth in vivo. We show that SIX1 maintains the RMS undifferentiated state by controlling enhancer activity and MYOD1 occupancy at loci more permissive to tumor growth over muscle differentiation. Finally, we demonstrate that a gene signature derived from SIX1 loss correlates with differentiation status and predicts RMS progression in human disease. Our findings demonstrate a master regulatory role of SIX1 in repression of RMS differentiation via genome-wide alterations in MYOD1 and MYOG-mediated transcription.

Abstract P177: NPEPPS regulates cisplatin-resistance and can be targeted to overcome treatment resistance in patient-derived bladder cancer tumoroids

Gemcitabine/cisplatin combination is the recommended neoadjuvant chemotherapy (NAC) for muscle-invasive bladder cancer (MIBC) patients undergoing surgical removal of the bladder (cystectomy). Patients without residual tumor at cystectomy ,complete pathological response (pCR), have a good 5-yr OS of 80%. However, due to NAC resistance, only 25% of the patients achieve a pCR. This study aims to increase the number of patients reaching pCR by discovery and functional interrogation of novel and pharmacologically targetable mechanisms causal to NAC resistance. We have identified an M1 aminopeptidase, NPEPPS, to regulate intracellular cisplatin import through volume regulated anion channels (VRACs). We show that genetic inhibition of NPEPPS increases the intracellular concentrations of cisplatin, hereby re-sensitizing gem/cis-resistant BC cell lines. On the other side, NPEPPS over-expression enhanced resistance to cisplatin. Pharmacological inhibition of NPEPPS with tosedostat results re-sensitization of gem/cis resistant BC cells in vitro and in vivo. Additional therapeutically validation was performed in gem/cis resistant patient-derived tumoroids. For this purpose, tumoroid cultures were generated from bladder cancer patients whom did not respond to pre-operative gem/cis, molecularly characterized and compared to the tumor of origin. Comparison of bladder cancer specific hotspot mutations, copy-number aberrations and H&E staining confirmed that patient-specific tumor traits were maintained in ex vivo tumoroid cultures. Moreover, tumoroid cultures generated from NAC-resistant patients were also resistant to cisplatin concentrations that exceeded physiological serum concentrations. Interestingly, addition of NPEPPS-inhibitor tosedostat sensitized these NAC resistant tumoroids to serum concentrations of cisplatin. These findings have potential for rapid translation into the clinic and invite trials investigating tosedostat to overcome chemoresistance.

Citation Format: Mathijs P. Scholtes, Maryam Akbarzadeh, Dan Theodorescu, James C. Costello, Tokameh Mahmoudi, Tahlita C.M. Zuiverloon. NPEPPS regulates cisplatin-resistance and can be targeted to overcome treatment resistance in patient-derived bladder cancer tumoroids [abstract]. In: Proceedings of the AACR-NCI-EORTC Virtual International Conference on Molecular Targets and Cancer Therapeutics; 2021 Oct 7-10. Philadelphia (PA): AACR; Mol Cancer Ther 2021;20(12 Suppl):Abstract nr P177.

IL13Rα2 Promotes Proliferation and Outgrowth of Breast Cancer Brain Metastases

PURPOSE

The survival of women with brain metastases (BM) from breast cancer remains very poor, with over 80% dying within a year of their diagnosis. Here, we define the function of IL13Rα2 in outgrowth of breast cancer brain metastases (BCBM) in vitro and in vivo, and postulate IL13Rα2 as a suitable therapeutic target for BM.

EXPERIMENTAL DESIGN

We performed IHC staining of IL13Rα2 in BCBM to define its prognostic value. Using inducible shRNAs in TNBC and HER2+ breast-brain metastatic models, we assessed IL13Rα2 function in vitro and in vivo. We performed RNAseq and functional studies to define the molecular mechanisms underlying IL13Rα2 function in BCBM.

RESULTS

High IL13Rα2 expression in BCBM predicted worse survival after BM diagnoses. IL13Rα2 was essential for cancer-cell survival, promoting proliferation while repressing invasion. IL13Rα2 KD resulted in FAK downregulation, repression of cell cycle and proliferation mediators, and upregulation of Ephrin B1 signaling. Ephrin-B1 (i) promoted invasion of BC cells in vitro, (ii) marked micrometastasis and invasive fronts in BCBM, and (iii) predicted shorter disease-free survival and BM-free survival (BMFS) in breast primary tumors known to metastasize to the brain. In experimental metastases models, which bypass early tumor invasion, downregulation of IL13Rα2 before or after tumor seeding and brain intravasation decreased BMs, suggesting that IL13Rα2 and the promotion of a proliferative phenotype is critical to BM progression.

CONCLUSIONS

Non-genomic phenotypic adaptations at metastatic sites are critical to BM progression and patients' prognosis. This study opens the road to use IL13Rα2 targeting as a therapeutic strategy for BM.

Mediator of DNA Damage Checkpoint 1 (MDC1) Is a Novel Estrogen Receptor Coregulator in Invasive Lobular Carcinoma of the Breast

Invasive lobular carcinoma (ILC) is the most common special histologic subtype of breast cancer, and nearly all ILC tumors express estrogen receptor alpha (ER). However, clinical and laboratory data suggest ILC are strongly estrogen-driven but not equally antiestrogen-sensitive. We hypothesized ILC-specific ER coregulators mediate ER functions and antiestrogen resistance in ILC, and profiled ER-associated proteins by mass spectrometry. Three ER+ ILC cell lines (MDA MB 134VI, SUM44PE, and BCK4) were compared with ER+ invasive ductal carcinoma (IDC) line data, and we examined whether siRNA of identified proteins suppressed ER-driven proliferation in ILC cells. This identified mediator of DNA damage checkpoint 1 (MDC1), a tumor suppressor in DNA damage response (DDR), as a novel ER coregulator in ILC. We confirmed ER:MDC1 interaction was specific to ILC versus IDC cells, and found MDC1 knockdown suppressed ILC cell proliferation and tamoxifen resistance. Using RNA-sequencing, we found in ILC cells MDC1 knockdown broadly dysregulates the ER transcriptome, with ER:MDC1 target genes enriched for promoter hormone response elements. Importantly, our data are inconsistent with MDC1 tumor suppressor functions in DDR, but suggest a novel oncogenic role for MDC1 as an ER coregulator. Supporting this, in breast tumor tissue microarrays, MDC1 protein was frequently low or absent in IDC, but MDC1 loss was rare in ER+ ILC. ER:MDC1 interaction and MDC1 coregulator functions may underlie ER function in ILC and serve as targets to overcome antiestrogen resistance in ILC. IMPLICATIONS: MDC1 has novel ER coregulator activity in ILC, which may underlie ILC-specific ER functions, estrogen response, and antiestrogen resistance.

Abstract 615: Differential pathway analyses of BCG-treated T1HG bladder cancer using Philips OncoSignal: A pilot study

Intravesical instillations with Bacillus Calmette-Guérin (BCG) is the recommended treatment for T1 high-grade (HG) bladder cancer (BC) patients. Unfortunately, risk stratification is insufficient to identify patients at risk of BCG treatment failure. Another challenge is the worldwide BCG-shortage, which emphasizes the need for alternative therapeutic strategies. Thus, we aimed to investigate molecular signal transduction pathway activities using Philips' OncoSignal consumer-ready test in BCG treated T1HG BC in order to identify differentially expressed pathways that may predict clinical outcome or provide a rationale for alternative therapeutic strategies. Primary T1HG BC patients treated with ≥5/6 BCG induction instillations were retrospectively included. RNA-sequencing (50 million paired-end reads, containing ≥80% tumor) was performed on: 1) BCG-naïve tumors from patients with complete response to BCG at study inclusion (BCG-responders) and 2) BCG-naïve tumors with matching muscle-invasive bladder cancer (MIBC) recurrences during BCG therapy (BCG-failures). Philips OncoSignal analysis was performed to quantify functional signal transduction pathway activities. Molecular subtyping was performed using non-muscle invasive and muscle-invasive classifiers. Clinical and molecular subtyping results were combined with Philips OncoSignal to explore differences between BCG responders and BCG failures and between molecular subtypes. OncoSignal was performed in 20 tumors: 14 BCG-naïve (7 BCG-responders vs 7 BCG-failures) and 6 MIBC recurrences. Median follow-up for BCG-responders was 99 months (IQR 74-120); median follow-up for BCG-failures was 49 months (IQR 25-79), with a median time to progression of 14 months (IQR 10-41). Based on pilot data, MAPK and TGFβ pathway activity was significantly higher in tumors from BCG failures vs BCG-responders (both p&lt;0.05). Using TCGAs MIBC classifier, 10/14 pre-BCG samples were luminal-papillary (LumP), 1/14 was luminal-infiltrated (lum-inf) and 3/14 basal-squamous (bas-sq). In post-BCG MIBC, 2/6 tumors were lum-inf, 3/6 bas-sq and 1/6 neuronal. Using Erasmus MCs T1-BCG classifier (unpublished), 7/14 pre-BCG tumors were luminal-like, the other 7/14 were lum-inf/basal-like. All 6 post-BCG MIBCs were lum-inf/basal-like. Importantly, TCGAs bas-sq/lum-inf subtype and Erasmus MCs lum-inf/basal-like subtype had the highest MAPK, TGFβ and JAK-STAT1/2 activity (all p&lt;0.05). Based on our pilot study, results indicate that Philips OncoSignal is able to stratify BCG failures from BCG responders. Furthermore, OncoSignal results discern LumP tumors from lum-inf/basal-like tumors. Philips OncoSignal might be useful to stratify patients for alternative treatments. Further investigation on a large cohort of BC patients is required to determine whether Philips OncoSignal can be of value during clinical-decision making.

Citation Format: Florus C. De Jong, Teemu D. Laajala, Robert F. Hoedemaeker, Sébastien Rinaldetti, Jolien T. Mensink, Angelique C. Van Der Made, Deric K. Van der Schoot, Egbert R. Boevé, Ellen C. Zwarthoff, Joost L. Boormans, Dan Theodorescu, James C. Costello, Tahlita C. Zuiverloon. Differential pathway analyses of BCG-treated T1HG bladder cancer using Philips OncoSignal: A pilot study [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 615.

Abstract 614: Transcriptomic analysis of BCG-treated T1HG bladder cancer patients identifies an EMT-basal subgroup with immune suppressive characteristics at high risk of BCG-failure

Intravesical instillations with Bacillus Calmette-Guérin (BCG) are the recommended treatment in T1 high-grade (HG) bladder cancer (BC) patients. However, even with BCG treatment, 50% of patients develop a HG recurrence or progression to advanced disease. Current risk stratification is insufficient to identify patients at risk of BCG failure. We aimed to identify molecular predictors of BCG failure in T1HG BC patients with the objective to improve clinical decision-making. Primary T1HG BC patients treated with ≥5/6 BCG induction instillations were retrospectively included. RNA-sequencing (50 million paired-end reads, ≥80% tumor) was performed on centrally reviewed tumors from: 1) BCG-naïve tumors and recurrent tumors from patients that failed BCG therapy (BCG non-responders), matched to an equal number of 2) patients who had an ongoing complete response until study inclusion (BCG responders). Endpoints were HG recurrence-free survival defined as BCG-failure (RFS) and progression-free survival (PFS), assessed by survival analyses. BCG-failure was defined as development of a biopsy-proven HG recurrence or progression to muscle-muscle invasion. Using RNA-seq, we sought to identify genes, pathways and signatures that stratified patients with respect to these endpoints. 132 BCG-naïve T1HG tumors were sequenced and 45 post-BCG tumors. Median follow-up for all patients was 75 months (IQR 48-106), 99 months for n=63 BCG-responders (IQR 74-120) and 49 months for n=69 BCG-non-responders (IQR 25-79). Three molecular subtypes were identified: BCG1 (32%), BCG2 (40%) and BCG3 (28%). BCG3 patients had a significantly worse RFS (HR 2.4 p&lt;0.01) and PFS (HR 2.7; p&lt;0.01) compared to BCG1 or BCG2 patients. BCG3 tumors expressed high EMT and basal markers and had an immune suppressive tumor microenvironment. Immune deconvolution revealed significantly increased immune cell infiltration (i.e. B cells, T regulatory cells and tumor associated macrophages), while CD14 and FOXP3 were identified as significant regulons in BCG3 patients. BCG1 patients showed a favorable PFS and overexpression of genes associated with BCG processing and antigen presentation. BCG2 patients overexpressed luminal BC markers, such as FGFR3, showed MYC pathway activity and were mostly luminal papillary (consensus muscle-invasive BC classifier). Interestingly, post-BCG tumor recurrences were strongly enriched for the BCG3 subtype; 30/45 (67%) of recurring tumors were classified as BCG3. Gene expression profiling of BCG-naïve primary T1HG BC patients identified three molecular subtypes that corresponded to clinical outcome and response to BCG therapy. The EMT-basal/immune suppressive BCG3 patients are not candidates for treatment with BCG given the poor RFS and PFS and might be candidates for early radical cystectomy or immune-modulating therapy.

Citation Format: Florus C. De Jong, Teemu D. Laajala, Robert F. Hoedemaeker, Sébastien Rinaldetti, Jolien T. Mensink, Angelique C. Van der Made, Deric K. Van der Schoot, Egbert R. Boevé, Ellen C. Zwarthoff, Joost L. Boormans, Dan Theodorescu, James C. Costello, Tahlita C. Zuiverloon. Transcriptomic analysis of BCG-treated T1HG bladder cancer patients identifies an EMT-basal subgroup with immune suppressive characteristics at high risk of BCG-failure [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 614.

MAP3K7 Loss Drives Enhanced Androgen Signaling and Independently Confers Risk of Recurrence in Prostate Cancer with Joint Loss of CHD1

Prostate cancer genomic subtypes that stratify aggressive disease and inform treatment decisions at the primary stage are currently limited. Previously, we functionally validated an aggressive subtype present in 15% of prostate cancer characterized by dual deletion of MAP3K7 and CHD1. Recent studies in the field have focused on deletion of CHD1 and its role in androgen receptor (AR) chromatin distribution and resistance to AR-targeted therapy; however, CHD1 is rarely lost without codeletion of MAP3K7. Here, we show that in the clinically relevant context of co-loss of MAP3K7 and CHD1 there are significant, collective changes to aspects of AR signaling. Although CHD1 loss mainly impacts the expansion of the AR cistrome, loss of MAP3K7 drives increased AR target gene expression. Prostate cancer cell line models engineered to cosuppress MAP3K7 and CHD1 also demonstrated increased AR-v7 expression and resistance to the AR-targeting drug enzalutamide. Furthermore, we determined that low protein expression of both genes is significantly associated with biochemical recurrence (BCR) in a clinical cohort of radical prostatectomy specimens. Low MAP3K7 expression, however, was the strongest independent predictor for risk of BCR over all other tested clinicopathologic factors including CHD1 expression. Collectively, these findings illustrate the importance of MAP3K7 loss in a molecular subtype of prostate cancer that poses challenges to conventional therapeutic approaches. IMPLICATIONS: These findings strongly implicate MAP3K7 loss as a biomarker for aggressive prostate cancer with significant risk for recurrence that poses challenges for conventional androgen receptor-targeted therapies.

Androgen Receptor Signaling in Prostate Cancer Genomic Subtypes

While many prostate cancer (PCa) cases remain indolent and treatable, others are aggressive and progress to the metastatic stage where there are limited curative therapies. Androgen receptor (AR) signaling remains an important pathway for proliferative and survival programs in PCa, making disruption of AR signaling a viable therapy option. However, most patients develop resistance to AR-targeted therapies or inherently never respond. The field has turned to PCa genomics to aid in stratifying high risk patients, and to better understand the mechanisms driving aggressive PCa and therapy resistance. While alterations to the AR gene itself occur at later stages, genomic changes at the primary stage can affect the AR axis and impact response to AR-directed therapies. Here, we review common genomic alterations in primary PCa and their influence on AR function and activity. Through a meta-analysis of multiple independent primary PCa databases, we also identified subtypes of significantly co-occurring alterations and examined their combinatorial effects on the AR axis. Further, we discussed the subsequent implications for response to AR-targeted therapies and other treatments. We identified multiple primary PCa genomic subtypes, and given their differing effects on AR activity, patient tumor genetics may be an important stratifying factor for AR therapy resistance.

Crowdsourcing assessment of maternal blood multi-omics for predicting gestational age and preterm birth

Identification of pregnancies at risk of preterm birth (PTB), the leading cause of newborn deaths, remains challenging given the syndromic nature of the disease. We report a longitudinal multi-omics study coupled with a DREAM challenge to develop predictive models of PTB. The findings indicate that whole-blood gene expression predicts ultrasound-based gestational ages in normal and complicated pregnancies (r = 0.83) and, using data collected before 37 weeks of gestation, also predicts the delivery date in both normal pregnancies (r = 0.86) and those with spontaneous preterm birth (r = 0.75). Based on samples collected before 33 weeks in asymptomatic women, our analysis suggests that expression changes preceding preterm prelabor rupture of the membranes are consistent across time points and cohorts and involve leukocyte-mediated immunity. Models built from plasma proteomic data predict spontaneous preterm delivery with intact membranes with higher accuracy and earlier in pregnancy than transcriptomic models (AUROC = 0.76 versus AUROC = 0.6 at 27-33 weeks of gestation).

Induction of ADAM10 by Radiation Therapy Drives Fibrosis, Resistance, and Epithelial-to-Mesenchyal Transition in Pancreatic Cancer

Stromal fibrosis activates prosurvival and proepithelial-to-mesenchymal transition (EMT) pathways in pancreatic ductal adenocarcinoma (PDAC). In patient tumors treated with neoadjuvant stereotactic body radiation therapy (SBRT), we found upregulation of fibrosis, extracellular matrix (ECM), and EMT gene signatures, which can drive therapeutic resistance and tumor invasion. Molecular, functional, and translational analysis identified two cell-surface proteins, a disintegrin and metalloprotease 10 (ADAM10) and ephrinB2, as drivers of fibrosis and tumor progression after radiation therapy (RT). RT resulted in increased ADAM10 expression in tumor cells, leading to cleavage of ephrinB2, which was also detected in plasma. Pharmacologic or genetic targeting of ADAM10 decreased RT-induced fibrosis and tissue tension, tumor cell migration, and invasion, sensitizing orthotopic tumors to radiation killing and prolonging mouse survival. Inhibition of ADAM10 and genetic ablation of ephrinB2 in fibroblasts reduced the metastatic potential of tumor cells after RT. Stimulation of tumor cells with ephrinB2 FC protein reversed the reduction in tumor cell invasion with ADAM10 ablation. These findings represent a model of PDAC adaptation that explains resistance and metastasis after RT and identifies a targetable pathway to enhance RT efficacy. SIGNIFICANCE: Targeting a previously unidentified adaptive resistance mechanism to radiation therapy in PDAC tumors in combination with radiation therapy could increase survival of the 40% of PDAC patients with locally advanced disease.See related commentary by Garcia Garcia et al., p. 3158 GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/81/12/3255/F1.large.jpg.

Inhibition of the CCL2 receptor, CCR2, enhances tumor response to immune checkpoint therapy

Immunotherapies targeting the PD-1/PD-L1 axis are now a mainstay in the clinical management of multiple cancer types, however, many tumors still fail to respond. CCL2 is highly expressed in various cancer types and has been shown to be associated with poor prognosis. Inhibition or blockade of the CCL2/CCR2 signaling axis has thus been an area of interest for cancer therapy. Here we show across multiple murine tumor and metastasis models that CCR2 antagonism in combination with anti-PD-1 therapy leads to sensitization and enhanced tumor response over anti-PD-1 monotherapy. We show that enhanced treatment response correlates with enhanced CD8+ T cell recruitment and activation and a concomitant decrease in CD4+ regulatory T cell. These results provide strong preclinical rationale for further clinical exploration of combining CCR2 antagonism with PD-1/PD-L1-directed immunotherapies across multiple tumor types especially given the availability of small molecule CCR2 inhibitors and antibodies.

Biologically Informed Neural Networks Predict Drug Responses

Deep neural networks often achieve high predictive accuracy on biological problems, but it can be hard to contextualize how and explain why predictions are made. In this issue, Kuenzi et al. model the sensitivity of cancers to drugs using deep neural networks with a hierarchical structure derived from the Gene Ontology.

KDM5A and PHF2 positively control expression of pro-metastatic genes repressed by EWS/Fli1, and promote growth and metastatic properties in Ewing sarcoma

Ewing sarcoma is an aggressive malignant neoplasm with high propensity for metastasis and poor clinical outcomes. The EWS/Fli1 oncofusion is the disease driver in > 90% of cases, but presents a difficult therapeutic target. Moreover, EWS/Fli1 plays a complex role in disease progression, with inhibitory effects on critical steps of metastasis. Like many other pediatric cancers, Ewing sarcoma is a disease marked by epigenetic dysregulation. Epigenetic mechanisms present alternative targeting opportunities, but their contributions to Ewing sarcoma metastasis and disease progression remain poorly understood. Here, we show that the epigenetic regulators KDM5A and PHF2 promote growth and metastatic properties in Ewing sarcoma, and, strikingly, activate expression many pro-metastatic genes repressed by EWS/Fli1. These genes include L1CAM, which is associated with adverse outcomes in Ewing sarcoma, and promotes migratory and invasive properties. KDM5A and PHF2 retain their growth promoting effects in more metastatically potent EWS/Fli1^low cells, and PHF2 promotes both invasion and L1CAM expression in this cell population. Furthermore, KDM5A and PHF2 each contribute to the increased metastatic potency of EWS/Fli1^low cells in vivo. Together, these studies identify KDM5A and PHF2 as novel disease-promoting factors, and potential new targets, in Ewing sarcoma, including the more metastatically potent EWS/Fli1^low cell population.

The Capacity of the Ovarian Cancer Tumor Microenvironment to Integrate Inflammation Signaling Conveys a Shorter Disease-free Interval

PURPOSE

Ovarian cancer has one of the highest deaths to incidence ratios across all cancers. Initial chemotherapy is effective, but most patients develop chemoresistant disease. Mechanisms driving clinical chemo-response or -resistance are not well-understood. However, achieving optimal surgical cytoreduction improves survival, and cytoreduction is improved by neoadjuvant chemotherapy (NACT). NACT offers a window to profile pre- versus post-NACT tumors, which we used to identify chemotherapy-induced changes to the tumor microenvironment.

EXPERIMENTAL DESIGN

We obtained matched pre- and post-NACT archival tumor tissues from patients with high-grade serous ovarian cancer (patient, n = 6). We measured mRNA levels of 770 genes (756 genes/14 housekeeping genes, NanoString Technologies), and performed reverse phase protein array (RPPA) on a subset of matched tumors. We examined cytokine levels in pre-NACT ascites samples (n = 39) by ELISAs. A tissue microarray with 128 annotated ovarian tumors expanded the transcriptional, RPPA, and cytokine data by multispectral IHC.

RESULTS

The most upregulated gene post-NACT was IL6 (16.79-fold). RPPA data were concordant with mRNA, consistent with elevated immune infiltration. Elevated IL6 in pre-NACT ascites specimens correlated with a shorter time to recurrence. Integrating NanoString (n = 12), RPPA (n = 4), and cytokine (n = 39) studies identified an activated inflammatory signaling network and induced IL6 and IER3 (immediate early response 3) post-NACT, associated with poor chemo-response and time to recurrence.

CONCLUSIONS

Multiomics profiling of ovarian tumor samples pre- and post-NACT provides unique insight into chemo-induced changes to the tumor microenvironment. We identified a novel IL6/IER3 signaling axis that may drive chemoresistance and disease recurrence.

Abstract 5305: Characterization of the ovarian tumor transcriptome and microenvironment in pre and post-chemotherapy treated patients

Purpose: Ovarian cancer has one of the highest deaths to incidence ratios across all cancers. Initial chemotherapy is typically effective, but most patients will experience recurrence and enter a cycle of response and recurrence, ultimately developing resistant disease. Importantly, the understanding of mechanisms contributing to clinical chemoresistance in ovarian cancer is limited. Achieving optimal cytoreduction improves survival. Gynecologic oncologists therefore perform laparoscopy to evaluate resectability and collect tumor biopsies. Patients identified with un-resectable tumor receive three cycles of neoadjuvant platinum/taxane-based chemotherapy (NACT) to improve the likelihood of optimal cytoreductive surgery. NACT offers a previously unavailable window for identifying therapy-induced remodeling of the tumor, potentially providing insight into drivers of chemo-resistance.

Experimental Design: We accessed formalin fixed (FFPE) tumor tissue from high-grade serous ovarian cancer NACT patients (n=6) pre- and post-chemotherapy. RNA was extracted and the mRNA transcript levels were examined for 791 genes via the Immunoncology 360 NanoString panel. Further, using the pre-treated ascites fluid we examined the cytokine environment through MesoScale Discovery multiplex ELISA. Transcriptional and cytokine profiles were correlated to clinical outcomes, such as days to recurrence.

Results: Component analysis found that transcriptional profiles separated based on pre- and post-chemotherapy status. Comparing the pre- (n=6) and post-chemotherapy (n=6), the most upregulated gene following therapy was IL6 (4.1 log2 fold change, adj. p = 0.045) and the most downregulated gene was UBE2C (-3.9 log2 fold change, adj. p = 0.001). Elevated IL-6 within ascites correlated to a shorter time to recurrence. On the patient-by-patient basis we compared transcriptional changes between pre- and post therapy. We correlated these changes to IL6 expression and time to recurrence uncovered that an increase in Immediate Early Response 3 (IER3) is a poor prognostic indicator. IER3/IEX-1 is a NF-κB target gene involved in apoptosis and potentiates MAPK activation.

Conclusion: Using ovarian cancer tumors from NACT patients provides a unique insight into chemo-induced transcriptional changes. To-date little is known of IER3's function in ovarian cancer thus future work will focus on a potential novel IL-6/IER3 signaling axis that protects tumors from chemotherapy-induced apoptosis. Also, computer modeling and integration of transcriptional and cytokine date could also provide an approach to predict recurrence.

Citation Format: Kimberly R. Jordan, Jill Slansky, Angela Minic, Jennifer K. Richer, James C. Costello, Aaron J. Clauset, Rajendra T. Kumar, Kian Behbakht, Matthew J. Sikora, Benjamin G. Bitler. Characterization of the ovarian tumor transcriptome and microenvironment in pre and post-chemotherapy treated patients [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5305.

KDM3A/Ets1 epigenetic axis contributes to PAX3/FOXO1-driven and independent disease-promoting gene expression in fusion-positive Rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is the most common soft tissue sarcoma in children and young adults. RMS exists as two major disease subtypes, oncofusion-negative RMS (FN-RMS) and oncofusion-positive RMS (FP-RMS). FP-RMS is characterized by recurrent PAX3/7-FOXO1 driver oncofusions and is a biologically and clinically aggressive disease. Recent studies have revealed FP-RMS to have a strong epigenetic basis. Epigenetic mechanisms represent potential new therapeutic vulnerabilities in FP-RMS, but their complex details remain to be defined. We previously identified a new disease-promoting epigenetic axis in RMS, involving the chromatin factor KDM3A and the Ets1 transcription factor. In the present study, we define the KDM3A and Ets1 FP-RMS transcriptomes and show that these interface with the recently characterized PAX3/FOXO1-driven gene expression program. KDM3A and Ets1 positively control numerous known and candidate novel PAX3/FOXO1-induced RMS-promoting genes, including subsets under control of PAX3/FOXO1-associated superenhancers (SE), such as MEST. Interestingly, KDM3A and Ets1 also positively control a number of known and candidate novel FP-RMS-promoting, but not PAX3/FOXO1-dependent, genes. Epistatically, Ets1 is downstream of, and exerts disease-promoting effects similar to, both KDM3A and PAX3/FOXO1. MEST also manifests disease-promoting properties in FP-RMS, and KDM3A and Ets1 each impacts activation of the PAX3/FOXO1-associated MEST SE. Taken together, our studies show that the KDM3A/Ets1 epigenetic axis plays an important role in disease promotion in FP-RMS, and provide insight into potential new ways to target aggressive phenotypes in this disease.

Abstract IA10: Developing rational combination therapy with checkpoint inhibitors

Targeting antibodies to programmed cell death protein-1 (PD-1) is an effective treatment across multiple cancer types. While a subset of patients receiving these therapies experience favorable responses, many still show disease progression, highlighting the importance of other mechanisms influencing immune responsiveness in these tumors. Therefore, combining therapies that enhance antitumor immunity has been an area of great interest to the entire cancer community. We have recently tackled this challenge in the rapidly evolving field of cancer immunotherapy by using in vivo functional genomics to identify genes whose inhibition potentiates the response to anti-PD-1 immunotherapy. Using an in vivo screening approach with a customized shRNA pooled library, we identified a number of candidates including DDR2 as promising targets for the enhancement of response to anti-PD-1 immunotherapy. In the case of DDR2, using isogenic in vivo murine models across five different tumor histologies—bladder, breast, colon, sarcoma, and melanoma—we show that DDR2 depletion increases sensitivity to anti-PD-1 treatment compared to monotherapy. Combination treatment of tumor-bearing mice with anti-PD-1 and dasatinib, a tyrosine kinase inhibitor of DDR2, also led to tumor load reduction and in some cases, complete clearance. RNAseq and CyTOF analysis revealed higher CD8+ T-cell populations in tumors with DDR2 depletion and those treated with dasatinib when either was combined with anti-PD-1 treatment. Our work provides strong scientific rationale for targeting DDR2 in combination with PD-1 inhibitors. In addition, a number of other potential druggable targets have been identified in our screen that we are currently pursuing.

Citation Format: Megan M. Tu, Francis Y. F. Lee, Robert T. Jones, Abigail K. Kimball, Elizabeth Saravia, Robert F. Graziano, Brianne Coleman, Krista Menard, Jun Yan, Erin Michaud, Han Chang, Hany A. Abdel-Hafiz, Andrii I. Rozhok, Jason E. Duex, Neeraj Agarwal, Ana Chauca-Diaz, Linda K. Johnson, Terry L. Ng, John C. Cambier, Eric T. Clambey, James C. Costello, Alan J. Korman, Dan Theodorescu. Developing rational combination therapy with checkpoint inhibitors [abstract]. In: Proceedings of the AACR Special Conference on Bladder Cancer: Transforming the Field; 2019 May 18-21; Denver, CO. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(15_Suppl):Abstract nr IA10.

Abstract B11: Multi-omic interrogation of gemcitabine and cisplatin-resistant bladder cancer cell lines identifies unique and shared mediators of chemosensitivity and resistance

Neoadjuvant platinum-based chemotherapy (NAC) followed by radical cystectomy is the preferred first-line option for treating patients diagnosed with muscle-invasive bladder cancer (MIBC). One of the two most commonly utilized chemotherapy regimens for MIBC consists of combination therapy with gemcitabine and cisplatin (GC). Importantly, cisplatin-based chemotherapy leads to complete response (pT0) in approximately 25% of patients, while the remaining 75% are left with residual disease. For MIBC patients, the presence of residual disease at cystectomy is associated with significantly poorer prognosis. Previous efforts to screen tumor-derived specimens have identified a number of mutational biomarkers associated with NAC response; however, to date these efforts have been limited by the number of patients profiled and the molecular and mutational heterogeneity intrinsic to bladder cancers. Moreover, these efforts have primarily relied on a single technology (e.g., whole-exome sequencing or transcriptome profiling), and were not designed to look at GC-resistance patterns in a way that integrates these data types. Therefore, we sought to systematically define the molecular phenotypes associated with chemoresistant bladder cancer cells, and to functionally interrogate these features to better define essential genes and pathways that govern chemoresponsiveness. To accomplish this, we have extensively characterized a panel of six bladder cancer cell lines (KU-19-19, T24, TCCSUP, 253J, 5637, and RT112), each of which has three derivative cell lines with acquired resistance to gemcitabine, cisplatin, and GC combination. For all of these cell lines we have performed whole-transcriptome mRNA sequencing, whole-exome sequencing, and reverse phase protein microarray. Additionally, we have performed functional genomic screening using whole-genome CRISPR knockout libraries, in the combination GC-resistant derivative cell lines. Using these data we identified multiple important chemoresistance genes and pathways, both previously known and novel. Interestingly, even though we found heterogeneous expression patterns seen across different –omic platforms, the functional genomic screens in GC-resistant cell lines converged on many genes and pathways that are essential to DNA damage response. These data suggest that functional genomic screens are perhaps more likely to identify genes and pathways that are directly linked to each of these drugs' specific mechanisms of action, as opposed to context-specific molecular differences seen across cell lines. Importantly, these results could enable improved prediction and personalization of therapy through identifying biomarkers that predict chemosensitivity and resistance. This study also suggests that identifying more generalizable strategies to overcome chemoresistance could increase the number of bladder cancer patients experiencing deeper and more complete responses to these common chemotherapeutic drugs.

Citation Format: Robert T. Jones, Tahlita C. M. Zuiverloon, Hedvig Vekony, Andrew Goodspeed, Teemu D. Laajala, Molishree Joshi, Colin Sempeck, Annie Jean, Megan Tu, Julia D. Wulfkuhle, Emanuel F. Petricoin, James C. Costello, Dan Theodorescu. Multi-omic interrogation of gemcitabine and cisplatin-resistant bladder cancer cell lines identifies unique and shared mediators of chemosensitivity and resistance [abstract]. In: Proceedings of the AACR Special Conference on Bladder Cancer: Transforming the Field; 2019 May 18-21; Denver, CO. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(15_Suppl):Abstract nr B11.

Abstract B08: MSH2 downregulation induces chemoresistance in bladder cancer

Treatment of muscle-invasive bladder cancer has seen little improvement in three decades. The objective of this study is to identify genes and pathways that associate with cisplatin resistance in bladder cancer. A whole-genome, CRISPR-based screen was performed in the MGHU4 bladder cancer cell line treated with cisplatin to identify genes that mediate resistance to cisplatin. Targeted validation was performed in vitro in another bladder cancer cell line, 253J. MSH2 and MLH1 are prominent components of the mismatch repair. These are the top two genes that have the largest effect on cisplatin resistance according to our synthetic lethal screen. Bladder cancer cells with an MSH2 knockdown are resistant to cisplatin in vitro, in part due to a reduction in apoptosis. MSH2 was also investigated in a publicly available bladder cancer dataset containing 340 bladder cancer patients with treatment, protein, and survival information. Bladder cancer patients with low protein levels of MSH2 have poorer overall survival when treated with cisplatin- or carboplatin-based therapy. If confirmed in follow-up studies, MSH2 protein level may serve as a prospective biomarker of chemotherapy response in this cancer type.

Note: This abstract was not presented at the conference.

Citation Format: Andrew E. Goodspeed, Annie Jean, James C. Costello. MSH2 downregulation induces chemoresistance in bladder cancer [abstract]. In: Proceedings of the AACR Special Conference on Bladder Cancer: Transforming the Field; 2019 May 18-21; Denver, CO. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(15_Suppl):Abstract nr B08.

Mesenchymal and MAPK Expression Signatures Associate with Telomerase Promoter Mutations in Multiple Cancers

In a substantial fraction of cancers TERT promoter (TERTp) mutations drive expression of the catalytic subunit of telomerase, contributing to their proliferative immortality. We conducted a pan-cancer analysis of cell lines and find a TERTp mutation expression signature dominated by epithelial-to-mesenchymal transition and MAPK signaling. These data indicate that TERTp mutants are likely to generate distinctive tumor microenvironments and intercellular interactions. Analysis of high-throughput screening tests of 546 small molecules on cell line growth indicated that TERTp mutants displayed heightened sensitivity to specific drugs, including RAS pathway inhibitors, and we found that inhibition of MEK1 and 2, key RAS/MAPK pathway effectors, inhibited TERT mRNA expression. Consistent with an enrichment of mesenchymal states in TERTp mutants, cell lines and some patient tumors displayed low expression of the central adherens junction protein E-cadherin, and we provide evidence that its expression in these cells is regulated by MEK1/2. Several mesenchymal transcription factors displayed elevated expression in TERTp mutants including ZEB1 and 2, TWIST1 and 2, and SNAI1. Of note, the developmental transcription factor SNAI2/SLUG was conspicuously elevated in a significant majority of TERTp-mutant cell lines, and knock-down experiments suggest that it promotes TERT expression. IMPLICATIONS: Cancers harboring TERT promoter mutations are often more lethal, but the basis for this higher mortality remains unknown. Our study identifies that TERTp mutants, as a class, associate with a distinct gene and protein expression signature likely to impact their biological and clinical behavior and provide new directions for investigating treatment approaches for these cancers.

Identification of a Small-Molecule Inhibitor That Disrupts the SIX1/EYA2 Complex, EMT, and Metastasis

Metastasis is the major cause of mortality for patients with cancer, and dysregulation of developmental signaling pathways can significantly contribute to the metastatic process. The Sine oculis homeobox homolog 1 (SIX1)/eyes absent (EYA) transcriptional complex plays a critical role in the development of multiple organs and is typically downregulated after development is complete. In breast cancer, aberrant expression of SIX1 has been demonstrated to stimulate metastasis through activation of TGFβ signaling and subsequent induction of epithelial-mesenchymal transition (EMT). In addition, SIX1 can induce metastasis via non-cell autonomous means, including activation of GLI-signaling in neighboring tumor cells and activation of VEGFC-induced lymphangiogenesis. Thus, targeting SIX1 would be expected to inhibit metastasis while conferring limited side effects. However, transcription factors are notoriously difficult to target, and thus novel approaches to inhibit their action must be taken. Here we identified a novel small molecule compound, NCGC00378430 (abbreviated as 8430), that reduces the SIX1/EYA2 interaction. 8430 partially reversed transcriptional and metabolic profiles mediated by SIX1 overexpression and reversed SIX1-induced TGFβ signaling and EMT. 8430 was well tolerated when delivered to mice and significantly suppressed breast cancer-associated metastasis in vivo without significantly altering primary tumor growth. Thus, we have demonstrated for the first time that pharmacologic inhibition of the SIX1/EYA2 complex and associated phenotypes is sufficient to suppress breast cancer metastasis. SIGNIFICANCE: These findings identify and characterize a novel inhibitor of the SIX1/EYA2 complex that reverses EMT phenotypes suppressing breast cancer metastasis.

OR12-05 MDC1 Is a Novel Estrogen Receptor Co-Regulator in Invasive Lobular Carcinoma of the Breast

Invasive Lobular Carcinoma (ILC) is the 2^nd most common histotype of breast cancer, but is critically understudied. ~95% of ILC are estrogen receptor (ER) positive, and previous studies demonstrate the importance of estrogen in ILC etiology. However, retrospective studies show that anti-estrogens are substantially less effective in ILC than in ER+ Invasive Ductal Carcinoma (IDC). This strongly suggests that regulation of ER function is unique in ILC, and we hypothesize that this is due to an ILC-specific cohort of ER co-regulators. We performed Rapid Immunoprecipitation Mass Spectrometry of Endogenous Proteins (RIME) to determine ILC-specific ER-interacting proteins, and identified Mediator of DNA Damage Checkpoint 1 (MDC1) as a novel ER co-regulator in ILC cells. We confirmed ER:MDC1 interaction by co-immunoprecipitation and proximity ligation assays (PLA); interaction was specifically observed in ILC cell lines but not IDC cell lines. Consistent with co-regulator function, we found MDC1 is essential for ER-driven proliferation of ILC cells. MDC1 knockdown dysregulates transcription of ER target genes in ILC cells (e.g. IGFBP4, WNT4). Moreover, RNA-seq analysis showed that in ILC cell line MDA MB 134VI, >50% of ER target genes require MDC1 for their regulation. To understand how MDC1 controls ER transcriptional activity, we performed ChIP-qPCR and found that MDC1 controls ER binding to DNA in ILC cells. Further, MDC1 controls binding of the pioneer factor FOXA1 to DNA, and Dual PLA studies of ER:MDC1 and ER:FOXA1 interaction revealed that MDC1 knockdown decreased ER:FOXA1 interaction. MDC1 canonically functions in DNA damage response, but our preliminary data suggest MDC1 is decoupled from its canonical role in DDR in the context of ER co-regulator activity in ILC cells. Together, these data suggest MDC1, independent of its role in DDR, acts as a novel ER co-regulator in ILC and regulates ER:DNA binding and ER transcriptional function to drive ILC cell proliferation and survival.

Monocytic Subclones Confer Resistance to Venetoclax-Based Therapy in Patients with Acute Myeloid Leukemia

Venetoclax-based therapy can induce responses in approximately 70% of older previously untreated patients with acute myeloid leukemia (AML). However, up-front resistance as well as relapse following initial response demonstrates the need for a deeper understanding of resistance mechanisms. In the present study, we report that responses to venetoclax +azacitidine in patients with AML correlate closely with developmental stage, where phenotypically primitive AML is sensitive, but monocytic AML is more resistant. Mechanistically, resistant monocytic AML has a distinct transcriptomic profile, loses expression of venetoclax target BCL2, and relies on MCL1 to mediate oxidative phosphorylation and survival. This differential sensitivity drives a selective process in patients which favors the outgrowth of monocytic subpopulations at relapse. Based on these findings, we conclude that resistance to venetoclax + azacitidine can arise due to biological properties intrinsic to monocytic differentiation. We propose that optimal AML therapies should be designed so as to independently target AML subclones that may arise at differing stages of pathogenesis. SIGNIFICANCE: Identifying characteristics of patients who respond poorly to venetoclax-based therapy and devising alternative therapeutic strategies for such patients are important topics in AML. We show that venetoclax resistance can arise due to intrinsic molecular/metabolic properties of monocytic AML cells and that such properties can potentially be targeted with alternative strategies.

Trisomy 21 activates the kynurenine pathway via increased dosage of interferon receptors

Trisomy 21 (T21) causes Down syndrome (DS), affecting immune and neurological function by ill-defined mechanisms. Here we report a large metabolomics study of plasma and cerebrospinal fluid, showing in independent cohorts that people with DS produce elevated levels of kynurenine and quinolinic acid, two tryptophan catabolites with potent immunosuppressive and neurotoxic properties, respectively. Immune cells of people with DS overexpress IDO1, the rate-limiting enzyme in the kynurenine pathway (KP) and a known interferon (IFN)-stimulated gene. Furthermore, the levels of IFN-inducible cytokines positively correlate with KP dysregulation. Using metabolic tracing assays, we show that overexpression of IFN receptors encoded on chromosome 21 contribute to enhanced IFN stimulation, thereby causing IDO1 overexpression and kynurenine overproduction in cells with T21. Finally, a mouse model of DS carrying triplication of IFN receptors exhibits KP dysregulation. Together, our results reveal a mechanism by which T21 could drive immunosuppression and neurotoxicity in DS.

Modeling genetic heterogeneity of drug response and resistance in cancer

Heterogeneity in tumors is recognized as a key contributor to drug resistance and spread of advanced disease, but deep characterization of genetic variation within tumors has only recently been quantifiable with the advancement of next generation sequencing and single cell technologies. These data have been essential in developing molecular models of how tumors develop, evolve, and respond to environmental changes, such as therapeutic intervention. A deeper understanding of tumor evolution has subsequently opened up new research efforts to develop mathematical models that account for evolutionary dynamics with the goal of predicting drug response and resistance in cancer. Here, we describe recent advances and limitations of how models of tumor evolution can impact treatment strategies for cancer patients.

Improved inference of chromosome conformation from images of labeled loci

We previously published a method that infers chromosome conformation from images of fluorescently-tagged genomic loci, for the case when there are many loci labeled with each distinguishable color.  Here we build on our previous work and improve the reconstruction algorithm to address previous limitations.  We show that these improvements 1) increase the reconstruction accuracy and 2) allow the method to be used on large-scale problems involving several hundred labeled loci.  Simulations indicate that full-chromosome reconstructions at 1/2 Mb resolution are possible using existing labeling and imaging technologies.  The updated reconstruction code and the script files used for this paper are available at:  https://github.com/heltilda/align3d.

GSEA-InContext: identifying novel and common patterns in expression experiments

Motivation

Gene Set Enrichment Analysis (GSEA) is routinely used to analyze and interpret coordinate pathway-level changes in transcriptomics experiments. For an experiment where less than seven samples per condition are compared, GSEA employs a competitive null hypothesis to test significance. A gene set enrichment score is tested against a null distribution of enrichment scores generated from permuted gene sets, where genes are randomly selected from the input experiment. Looking across a variety of biological conditions, however, genes are not randomly distributed with many showing consistent patterns of up- or down-regulation. As a result, common patterns of positively and negatively enriched gene sets are observed across experiments. Placing a single experiment into the context of a relevant set of background experiments allows us to identify both the common and experiment-specific patterns of gene set enrichment.

Results

We compiled a compendium of 442 small molecule transcriptomic experiments and used GSEA to characterize common patterns of positively and negatively enriched gene sets. To identify experiment-specific gene set enrichment, we developed the GSEA-InContext method that accounts for gene expression patterns within a background set of experiments to identify statistically significantly enriched gene sets. We evaluated GSEA-InContext on experiments using small molecules with known targets to show that it successfully prioritizes gene sets that are specific to each experiment, thus providing valuable insights that complement standard GSEA analysis.

Availability and implementation

GSEA-InContext implemented in Python, Supplementary results and the background expression compendium are available at: https://github.com/CostelloLab/GSEA-InContext.

Loss of MAP3K7 Sensitizes Prostate Cancer Cells to CDK1/2 Inhibition and DNA Damage by Disrupting Homologous Recombination

The combined loss of CHD1 and MAP3K7 promotes aggressive prostate cancer by unknown mechanisms. Because both of these genes are lost genetically in prostate cancer, they cannot be directly targeted. We applied an established computational systems pharmacology approach (TRAP) to identify altered signaling pathways and associated druggable targets. We compared gene expression profiles of prostate cancer with coloss of CHD1 and MAP3K7 with prostate cancer diploid for these genes using The Cancer Genome Atlas patient samples. This analysis prioritized druggable target genes that included CDK1 and CDK2. We validated that inhibitors of these druggable target genes, including the CDK1/CDK2 inhibitor dinaciclib, had antiproliferative and cytotoxic effects selectively on mouse prostate cells with knockdown of Chd1 and Map3k7. Dinaciclib had stronger effects on prostate cells with suppression of Map3k7 independent of Chd1 and also compared with cells without loss of Map3k7. Dinaciclib treatment reduced expression of homologous recombination (HR) repair genes such as ATM, ATR, BRCA2, and RAD51, blocked BRCA1 phosphorylation, reduced RAD51 foci formation, and increased γH2AX foci selectively in prostate cells with suppression of Map3k7, thus inhibiting HR repair of chromosomal double-strand breaks. Dinaciclib-induced HR disruption was also observed in human prostate cells with knockdown of MAP3K7. Cotreatment of dinaciclib with DNA-damaging agents or PARP inhibitor resulted in a stronger cytotoxic effect on prostate cells with suppression of MAP3K7 compared with those without loss of MAP3K7, or to each single agent. IMPLICATIONS: These findings demonstrate that loss of MAP3K7 is a main contributing factor to drug response through disruption of HR in prostate cancer.

A DREAM Challenge to Build Prediction Models for Short-Term Discontinuation of Docetaxel in Metastatic Castration-Resistant Prostate Cancer

Purpose

Docetaxel has a demonstrated survival benefit for patients with metastatic castration-resistant prostate cancer (mCRPC); however, 10% to 20% of patients discontinue docetaxel prematurely because of toxicity-induced adverse events, and the management of risk factors for toxicity remains a challenge.

Patients and Methods

The comparator arms of four phase III clinical trials in first-line mCRPC were collected, annotated, and compiled, with a total of 2,070 patients. Early discontinuation was defined as treatment stoppage within 3 months as a result of adverse treatment effects; 10% of patients discontinued treatment. We designed an open-data, crowd-sourced DREAM Challenge for developing models with which to predict early discontinuation of docetaxel treatment. Clinical features for all four trials and outcomes for three of the four trials were made publicly available, with the outcomes of the fourth trial held back for unbiased model evaluation. Challenge participants from around the world trained models and submitted their predictions. Area under the precision-recall curve was the primary metric used for performance assessment.

Results

In total, 34 separate teams submitted predictions. Seven models with statistically similar area under precision-recall curves (Bayes factor ≤ 3) outperformed all other models. A postchallenge analysis of risk prediction using these seven models revealed three patient subgroups: high risk, low risk, or discordant risk. Early discontinuation events were two times higher in the high-risk subgroup compared with the low-risk subgroup. Simulation studies demonstrated that use of patient discontinuation prediction models could reduce patient enrollment in clinical trials without the loss of statistical power.

Conclusion

This work represents a successful collaboration between 34 international teams that leveraged open clinical trial data. Our results demonstrate that routinely collected clinical features can be used to identify patients with mCRPC who are likely to discontinue treatment because of adverse events and establishes a robust benchmark with implications for clinical trial design.

Abstract 2107: Integrative molecular and functional genomic analysis of chemotherapy resistant bladder cancer cell lines identifies novel mediators of therapeutic response

Gemcitabine and cisplatin (GC) combination neoadjuvant chemotherapy (NAC) followed by radical cystectomy (RC) is the standard of care for treating muscle-invasive bladder cancer (MIBC). Approximately 25% of patients experience complete pathologic response (pT0), while the remaining ~75% have residual disease and poorer prognosis. Therefore, defining mechanisms that enable bladder cancer cells to survive NAC will allow us to I) identify predictive biomarkers of response to NAC and improve MIBC patient stratification for NAC, and II) define novel genetic dependencies in GC-resistant bladder cancers for future drug development. We have performed RNA sequencing on a panel of six genetically diverse urothelial cancer cell lines, each with matched derivatives possessing acquired resistance to single agent gemcitabine, cisplatin, and combination GC. Initial analyses of RNA sequencing data has revealed hundreds of differentially expressed genes and associated pathways as potential drivers of resistance to these agents. Examples of significantly dysregulated pathways seen in gene set enrichment analyses implicate a range of pathways known to be associated with many diverse cellular processes such as: extracellular matrix organization, invasion and migration, interferon signaling, lipid and cholesterol metabolism and cellular differentiation status, among others. Importantly, these candidates include both shared and context-specific resistance phenotypes, which vary based on the genetic and molecular alterations found in the parental cell line. To complement these transcriptomic profiling experiments and to better understand how these changes in gene expression relate to mutation status and pathway activation status on the protein level, we have performed whole-exome sequencing, and are currently conducting phospho-proteomic profiling with reverse-phase protein microarray to enable more integrative analyses. Together the complement of transcriptomic, genomic and proteomic analyses will facilitate a more comprehensive assessment of which molecular patterns are drivers of resistance in our resistant cell lines. Finally, we have taken a more systematic approach to define functional genetic dependencies acquired in these NAC-resistant cells. To accomplish this we have performed whole-genome CRISPR knockout screening on the combined GC-resistant cells to identify vulnerabilities associated with GC-resistance. Taken together, our data provide a deep look into the context-specific and shared molecular mechanisms driving resistance to GC chemotherapy and provide a rational basis for the future development of molecularly guided therapies for NAC resistant bladder cancers.

Citation Format: Robert T. Jones, Tahlita C. Zuiverloon, Hedvig Vekony, Andrew Goodspeed, Teemu D. Laajala, Molishree Joshi, Colin Sempeck, Megan Tu, James C. Costello, Dan Theodorescu. Integrative molecular and functional genomic analysis of chemotherapy resistant bladder cancer cell lines identifies novel mediators of therapeutic response [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 2107.

Abstract 4085: DDR2 inhibition enhances response to anti-PD-1 immunotherapy

Therapies targeting PD-1 are used in multiple cancer types. While a fraction of patients show durable therapeutic responses, most remain unresponsive, highlighting an urgent need to better understand and improve these therapies. Using an in vivo screening approach with a customized shRNA pooled library, we identified DDR2 as a leading target for the enhancement of response to anti-PD-1 immunotherapy. Using isogenic in vivo murine models across five different tumor histologies, bladder, breast, colon, sarcoma and melanoma, we show that DDR2 depletion increases sensitivity to anti-PD-1 treatment compared to monotherapy. Combination treatment of tumor-bearing mice with anti-PD-1 and dasatinib, a tyrosine kinase inhibitor of DDR2, also led to tumor load reduction and in some cases, complete clearance. RNAseq and CyTOF analysis revealed higher CD8+ T cell populations in tumors with DDR2 depletion and those treated with dasatinib when either was combined with anti-PD-1 treatment. Our work provides strong scientific rationale for targeting DDR2 in combination with PD-1 inhibitors.

Citation Format: Megan M. Tu, Francis Y. Lee, Robert T. Jones, Abigail K. Kimball, Elizabeth Saravia, Robert F. Graziano, Brianne Coleman, Krista Menard, Jun Yan, Erin Michaud, Han Chang, Hany A. Abdel-Hafiz, Andrii I. Rozhok, Jason E. Duex, Neeraj Agarwal, Ana Chauca-Diaz, Linda K. Johnson, Terry L. Ng, John C. Cambier, Eric T. Clambey, James C. Costello, Alan J. Korman, Dan Theodorescu. DDR2 inhibition enhances response to anti-PD-1 immunotherapy [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4085.

Targeting DDR2 enhances tumor response to anti-PD-1 immunotherapy

While a fraction of cancer patients treated with anti-PD-1 show durable therapeutic responses, most remain unresponsive, highlighting the need to better understand and improve these therapies. Using an in vivo screening approach with a customized shRNA pooled library, we identified DDR2 as a leading target for the enhancement of response to anti-PD-1 immunotherapy. Using isogenic in vivo murine models across five different tumor histologies-bladder, breast, colon, sarcoma, and melanoma-we show that DDR2 depletion increases sensitivity to anti-PD-1 treatment compared to monotherapy. Combination treatment of tumor-bearing mice with anti-PD-1 and dasatinib, a tyrosine kinase inhibitor of DDR2, led to tumor load reduction. RNA-seq and CyTOF analysis revealed higher CD8⁺ T cell populations in tumors with DDR2 depletion and those treated with dasatinib when either was combined with anti-PD-1 treatment. Our work provides strong scientific rationale for targeting DDR2 in combination with PD-1 inhibitors.

SIX2 Mediates Late-Stage Metastasis via Direct Regulation of SOX2 and Induction of a Cancer Stem Cell Program

The capacity for tumor cells to metastasize efficiently is directly linked to their ability to colonize secondary sites. Here we identify Six2, a developmental transcription factor, as a critical regulator of a breast cancer stem cell program that enables metastatic colonization. In several triple-negative breast cancer (TNBC) models, Six2 enhanced the expression of genes associated with embryonic stem cell programs. Six2 directly bound the Sox2 Srr2 enhancer, promoting Sox2 expression and downstream expression of Nanog, which are both key pluripotency factors. Regulation of Sox2 by Six2 enhanced cancer stem cell properties and increased metastatic colonization. Six2 and Sox2 expression correlated highly in breast cancers including TNBC, where a Six2 expression signature was predictive of metastatic burden and poor clinical outcome. Our findings demonstrate that a SIX2/SOX2 axis is required for efficient metastatic colonization, underscoring a key role for stemness factors in outgrowth at secondary sites. SIGNIFICANCE: These findings provide novel mechanistic insight into stemness and the metastatic outgrowth of triple-negative breast cancer cells.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/79/4/720/F1.large.jpg.