Identifying biomarkers of differential chemotherapy response in TNBC patient-derived xenografts with a CTD/WGCNA approach.

Although systemic chemotherapy remains the standard of care for TNBC, even combination chemotherapy is often ineffective. The identification of biomarkers for differential chemotherapy response would allow for the selection of responsive patients, thus maximizing efficacy and minimizing toxicities. Here, we leverage TNBC PDXs to identify biomarkers of response. To demonstrate their ability to function as a preclinical cohort, PDXs were characterized using DNA sequencing, transcriptomics, and proteomics to show consistency with clinical samples. We then developed a network-based approach (CTD/WGCNA) to identify biomarkers of response to carboplatin (MSI1, TMSB15A, ARHGDIB, GGT1, SV2A, SEC14L2, SERPINI1, ADAMTS20, DGKQ) and docetaxel (c, MAGED4, CERS1, ST8SIA2, KIF24, PARPBP). CTD/WGCNA multigene biomarkers are predictive in PDX datasets (RNAseq and Affymetrix) for both taxane- (docetaxel or paclitaxel) and platinum-based (carboplatin or cisplatin) response, thereby demonstrating cross-expression platform and cross-drug class robustness. These biomarkers were also predictive in clinical datasets, thus demonstrating translational potential.

Enrichment of Cancer Stem Cells in a Tumorsphere Assay.

Cancer stem cells (CSCs) are a small subpopulation of self-renewing cancer cells that are present within tumors. CSCs possess tumor initiation potential as well as the ability to resist toxic compounds and chemotherapeutic agents through the upregulation of drug efflux transporters, DNA repair pathways, and survival cascades. Accumulating evidence suggests that CSCs are responsible for tumor relapse and resistance to chemotherapeutic agents and that targeting CSCs is critical to inhibition of cancer progression. Therefore, isolation and characterization of CSCs is important in studying tumor initiation and progression. In this chapter, we provide a detailed method for the identification and isolation of CSCs.

In Vitro Quantification of Cancer Stem Cells Using a Mammosphere Formation Assay.

Cancer stem cells (CSCs) are a small subpopulation of self-renewing cancer cells that are present within tumors. In this chapter, we provide a detailed method for the quantification of CSCs in vitro through mammosphere formation.

Limiting Dilution Tumor Initiation Assay: An In Vivo Approach for the Study of Cancer Stem Cells.

Cancer stem cells (CSCs) are a small subpopulation of self-renewing cancer cells that are present within tumors. Calculating the frequency of tumor-initiating cells is important in the assessment of the number of CSCs present in a cell population. In this chapter, we present a protocol developed for quantification of CSCs from breast cancer tumors that can be adapted to CSCs from other types of tumors.

NPY1R exerts inhibitory action on estradiol-stimulated growth and predicts endocrine sensitivity and better survival in ER-positive breast cancer.

G Protein-Coupled Receptors (GPCRs) represent the largest superfamily of cell-surface proteins. However, the expression and function of majority of GPCRs remain unexplored in breast cancer (BC). We interrogated the expression and phosphorylation status of 398 non-sensory GPCRs using the landmark BC proteogenomics and phosphoproteomic dataset from The Cancer Genome Atlas. Neuropeptide Y Receptor Y1 (NPY1R) gene and protein expression were significantly higher in Luminal A tumors versus other BC subtypes. The trend of NPY1R gene, protein, and phosphosite (NPY1R-S368s) expression was decreasing in the order of Luminal A, Luminal B, Basal, and human epidermal growth factor receptor 2 (HER2) subtypes. NPY1R gene expression increased in response to estrogen and reduced with endocrine therapy in estrogen receptor-positive (ER+) BC cells and xenograft models. Conversely, NPY1R expression decreased in ER+ BC cells resistant to endocrine therapies (estrogen deprivation, tamoxifen, and fulvestrant) in vitro and in vivo. NPY treatment reduced estradiol-stimulated cell growth, which was reversed by NPY1R antagonist (BIBP-3226) in ER+ BC cells. Higher NPY1R gene expression predicted better relapse-free survival and overall survival in ER+ BC. Our study demonstrates that NPY1R mediates the inhibitory action of NPY on estradiol-stimulated growth of ER+ BC cells, and its expression serves as a biomarker to predict endocrine sensitivity and survival in ER+ BC patients.

A proteogenomic portrait of lung squamous cell carcinoma.

Lung squamous cell carcinoma (LSCC) remains a leading cause of cancer death with few therapeutic options. We characterized the proteogenomic landscape of LSCC, providing a deeper exposition of LSCC biology with potential therapeutic implications. We identify NSD3 as an alternative driver in FGFR1-amplified tumors and low-p63 tumors overexpressing the therapeutic target survivin. SOX2 is considered undruggable, but our analyses provide rationale for exploring chromatin modifiers such as LSD1 and EZH2 to target SOX2-overexpressing tumors. Our data support complex regulation of metabolic pathways by crosstalk between post-translational modifications including ubiquitylation. Numerous immune-related proteogenomic observations suggest directions for further investigation. Proteogenomic dissection of CDKN2A mutations argue for more nuanced assessment of RB1 protein expression and phosphorylation before declaring CDK4/6 inhibition unsuccessful. Finally, triangulation between LSCC, LUAD, and HNSCC identified both unique and common therapeutic vulnerabilities. These observations and proteogenomics data resources may guide research into the biology and treatment of LSCC.

A novel role of ADGRF1 (GPR110) in promoting cellular quiescence and chemoresistance in human epidermal growth factor receptor 2-positive breast cancer.

While G protein-coupled receptors (GPCRs) are known to be excellent drug targets, the second largest family of adhesion-GPCRs is less explored for their role in health and disease. ADGRF1 (GPR110) is an adhesion-GPCR and has an important function in neurodevelopment and cancer. Despite serving as a poor predictor of survival, ADGRF1's coupling to G proteins and downstream pathways remain unknown in cancer. We evaluated the effects of ADGRF1 overexpression on tumorigenesis and signaling pathways using two human epidermal growth factor receptor-2-positive (HER2+) breast cancer (BC) cell-line models. We also interrogated publicly available clinical datasets to determine the expression of ADGRF1 in various BC subtypes and its impact on BC-specific survival (BCSS) and overall survival (OS) in patients. ADGRF1 overexpression in HER2+ BC cells increased secondary mammosphere formation, soft agar colony formation, and % of Aldefluor-positive tumorigenic population in vitro and promoted tumor growth in vivo. ADGRF1 co-immunoprecipitated with both Gαs and Gαq proteins and increased cAMP and IP1 when overexpressed. However, inhibition of only the Gαs pathway by SQ22536 reversed the pro-tumorigenic effects of ADGRF1 overexpression. RNA-sequencing and RPPA analysis revealed inhibition of cell cycle pathways with ADGRF1 overexpression, suggesting cellular quiescence, as also evidenced by cell cycle arrest at the G0/1 phase and resistance to chemotherapy in HER2+ BC. ADGRF1 was significantly overexpressed in the HER2-enriched BC compared to luminal A and B subtypes and predicted worse BCSS and OS in these patients. Therefore, ADGRF1 represents a novel drug target in HER2+ BC, warranting discovery of novel ADGRF1 antagonists.

Identification of EMT signaling cross-talk and gene regulatory networks by single-cell RNA sequencing.

The epithelial-to-mesenchymal transition (EMT) plays a critical role during normal development and in cancer progression. EMT is induced by various signaling pathways, including TGF-ß, BMP, Wnt-ß-catenin, NOTCH, Shh, and receptor tyrosine kinases. In this study, we performed single-cell RNA sequencing on MCF10A cells undergoing EMT by TGF-ß1 stimulation. Our comprehensive analysis revealed that cells progress through EMT at different paces. Using pseudotime clustering reconstruction of gene-expression profiles during EMT, we found sequential and parallel activation of EMT signaling pathways. We also observed various transitional cellular states during EMT. We identified regulatory signaling nodes that drive EMT with the expression of important microRNAs and transcription factors. Using a random circuit perturbation methodology, we demonstrate that the NOTCH signaling pathway acts as a key driver of TGF-ß-induced EMT. Furthermore, we demonstrate that the gene signatures of pseudotime clusters corresponding to the intermediate hybrid EMT state are associated with poor patient outcome. Overall, this study provides insight into context-specific drivers of cancer progression and highlights the complexities of the EMT process.

CD8+ T cells inhibit metastasis and CXCL4 regulates its function.

BACKGROUND: The mechanism by which immune cells regulate metastasis is unclear. Understanding the role of immune cells in metastasis will guide the development of treatments improving patient survival. METHODS: We used syngeneic orthotopic mouse tumour models (wild-type, NOD/scid and Nude), employed knockout (CD8 and CD4) models and administered CXCL4. Tumours and lungs were analysed for cancer cells by bioluminescence, and circulating tumour cells were isolated from blood. Immunohistochemistry on the mouse tumours was performed to confirm cell type, and on a tissue microarray with 180 TNBCs for human relevance. TCGA data from over 10,000 patients were analysed as well. RESULTS: We reveal that intratumoral immune infiltration differs between metastatic and non-metastatic tumours. The non-metastatic tumours harbour high levels of CD8+ T cells and low levels of platelets, which is reverse in metastatic tumours. During tumour progression, platelets and CXCL4 induce differentiation of monocytes into myeloid-derived suppressor cells (MDSCs), which inhibit CD8+ T-cell function. TCGA pan-cancer data confirmed that CD8lowPlatelethigh patients have a significantly lower survival probability compared to CD8highPlateletlow. CONCLUSIONS: CD8+ T cells inhibit metastasis. When the balance between CD8+ T cells and platelets is disrupted, platelets produce CXCL4, which induces MDSCs thereby inhibiting the CD8+ T-cell function.

EMTome: a resource for pan-cancer analysis of epithelial-mesenchymal transition genes and signatures.

BACKGROUND: The epithelial-mesenchymal transition (EMT) enables dissociation of tumour cells from the primary tumour mass, invasion through the extracellular matrix, intravasation into blood vessels and colonisation of distant organs. Cells that revert to the epithelial state via the mesenchymal-epithelial transition cause metastases, the primary cause of death in cancer patients. EMT also empowers cancer cells with stem-cell properties and induces resistance to chemotherapeutic drugs. Understanding the driving factors of EMT is critical for the development of effective therapeutic interventions. METHODS: This manuscript describes the generation of a database containing EMT gene signatures derived from cell lines, patient-derived xenografts and patient studies across cancer types and multiomics data and the creation of a web-based portal to provide a comprehensive analysis resource. RESULTS: EMTome incorporates (i) EMT gene signatures; (ii) EMT-related genes with multiomics features across different cancer types; (iii) interactomes of EMT-related genes (miRNAs, transcription factors, and proteins); (iv) immune profiles identified from The Cancer Genome Atlas (TCGA) cohorts by exploring transcriptomics, epigenomics, and proteomics, and drug sensitivity and (iv) clinical outcomes of cancer cohorts linked to EMT gene signatures. CONCLUSION: The web-based EMTome portal is a resource for primary and metastatic tumour research publicly available at www.emtome.org .

Proteogenomic Characterization Reveals Therapeutic Vulnerabilities in Lung Adenocarcinoma.

To explore the biology of lung adenocarcinoma (LUAD) and identify new therapeutic opportunities, we performed comprehensive proteogenomic characterization of 110 tumors and 101 matched normal adjacent tissues (NATs) incorporating genomics, epigenomics, deep-scale proteomics, phosphoproteomics, and acetylproteomics. Multi-omics clustering revealed four subgroups defined by key driver mutations, country, and gender. Proteomic and phosphoproteomic data illuminated biology downstream of copy number aberrations, somatic mutations, and fusions and identified therapeutic vulnerabilities associated with driver events involving KRAS, EGFR, and ALK. Immune subtyping revealed a complex landscape, reinforced the association of STK11 with immune-cold behavior, and underscored a potential immunosuppressive role of neutrophil degranulation. Smoking-associated LUADs showed correlation with other environmental exposure signatures and a field effect in NATs. Matched NATs allowed identification of differentially expressed proteins with potential diagnostic and therapeutic utility. This proteogenomics dataset represents a unique public resource for researchers and clinicians seeking to better understand and treat lung adenocarcinomas.

Neurofibromin Is an Estrogen Receptor-α Transcriptional Co-repressor in Breast Cancer.

We report that neurofibromin, a tumor suppressor and Ras-GAP (GTPase-activating protein), is also an estrogen receptor-α (ER) transcriptional co-repressor through leucine/isoleucine-rich motifs that are functionally independent of GAP activity. GAP activity, in turn, does not affect ER binding. Consequently, neurofibromin depletion causes estradiol hypersensitivity and tamoxifen agonism, explaining the poor prognosis associated with neurofibromin loss in endocrine therapy-treated ER+ breast cancer. Neurofibromin-deficient ER+ breast cancer cells initially retain sensitivity to selective ER degraders (SERDs). However, Ras activation does play a role in acquired SERD resistance, which can be reversed upon MEK inhibitor addition, and SERD/MEK inhibitor combinations induce tumor regression. Thus, neurofibromin is a dual repressor for both Ras and ER signaling, and co-targeting may treat neurofibromin-deficient ER+ breast tumors.

Proteogenomic Characterization of Endometrial Carcinoma.

We undertook a comprehensive proteogenomic characterization of 95 prospectively collected endometrial carcinomas, comprising 83 endometrioid and 12 serous tumors. This analysis revealed possible new consequences of perturbations to the p53 and Wnt/ß-catenin pathways, identified a potential role for circRNAs in the epithelial-mesenchymal transition, and provided new information about proteomic markers of clinical and genomic tumor subgroups, including relationships to known druggable pathways. An extensive genome-wide acetylation survey yielded insights into regulatory mechanisms linking Wnt signaling and histone acetylation. We also characterized aspects of the tumor immune landscape, including immunogenic alterations, neoantigens, common cancer/testis antigens, and the immune microenvironment, all of which can inform immunotherapy decisions. Collectively, our multi-omic analyses provide a valuable resource for researchers and clinicians, identify new molecular associations of potential mechanistic significance in the development of endometrial cancers, and suggest novel approaches for identifying potential therapeutic targets.

Microscaled proteogenomic methods for precision oncology.

Cancer proteogenomics promises new insights into cancer biology and treatment efficacy by integrating genomics, transcriptomics and protein profiling including modifications by mass spectrometry (MS). A critical limitation is sample input requirements that exceed many sources of clinically important material. Here we report a proteogenomics approach for core biopsies using tissue-sparing specimen processing and microscaled proteomics. As a demonstration, we analyze core needle biopsies from ERBB2 positive breast cancers before and 48-72 h after initiating neoadjuvant trastuzumab-based chemotherapy. We show greater suppression of ERBB2 protein and both ERBB2 and mTOR target phosphosite levels in cases associated with pathological complete response, and identify potential causes of treatment resistance including the absence of ERBB2 amplification, insufficient ERBB2 activity for therapeutic sensitivity despite ERBB2 amplification, and candidate resistance mechanisms including androgen receptor signaling, mucin overexpression and an inactive immune microenvironment. The clinical utility and discovery potential of proteogenomics at biopsy-scale warrants further investigation.

Immune Checkpoint Profiles in Luminal B Breast Cancer (Alliance).

BACKGROUND: Unlike estrogen receptor (ER)-negative breast cancer, ER-positive breast cancer outcome is less influenced by lymphocyte content, indicating the presence of immune tolerance mechanisms that may be specific to this disease subset. METHODS: A supervised analysis of microarray data from the ACOSOG Z1031 (Alliance) neoadjuvant aromatase inhibitor (AI) trial identified upregulated genes in Luminal (Lum) B breast cancers that correlated with AI-resistant tumor proliferation (percentage of Ki67-positive cancer nuclei, Pearson r > 0.4) (33 cases Ki67 > 10% on AI) vs LumB breast cancers that were more AI sensitive (33 cases Ki67 < 10% on AI). Overrepresentation analysis was performed using WebGestalt. All statistical tests were two-sided. RESULTS: Thirty candidate genes positively correlated (r ≥ 0.4) with AI-resistant proliferation in LumB and were upregulated greater than twofold. Gene ontologies identified that the targetable immune checkpoint (IC) components IDO1, LAG3, and PD1 were overrepresented resistance candidates (P ≤ .001). High IDO1 mRNA was associated with poor prognosis in LumB disease (Molecular Taxonomy of Breast Cancer International Consortium, hazard ratio = 1.43, 95% confidence interval = 1.04 to 1.98, P = .03). IDO1 also statistically significantly correlated with STAT1 at protein level in LumB disease (Pearson r = 0.74). As a composite immune tolerance signature, expression of IFN-γ/STAT1 pathway components was associated with higher baseline Ki67, lower estrogen, and progesterone receptor mRNA levels and worse disease-specific survival (P = .002). In a tissue microarray analysis, IDO1 was observed in stromal cells and tumor-associated macrophages, with a higher incidence in LumB cases. Furthermore, IDO1 expression was associated with a macrophage mRNA signature (M1 by CIBERSORT Pearson r = 0.62 ) and by tissue microarray analysis. CONCLUSIONS: Targetable IC components are upregulated in the majority of endocrine therapy-resistant LumB cases. Our findings provide rationale for IC inhibition in poor-outcome ER-positive breast cancer.

Integrated Proteogenomic Characterization of Clear Cell Renal Cell Carcinoma.

To elucidate the deregulated functional modules that drive clear cell renal cell carcinoma (ccRCC), we performed comprehensive genomic, epigenomic, transcriptomic, proteomic, and phosphoproteomic characterization of treatment-naive ccRCC and paired normal adjacent tissue samples. Genomic analyses identified a distinct molecular subgroup associated with genomic instability. Integration of proteogenomic measurements uniquely identified protein dysregulation of cellular mechanisms impacted by genomic alterations, including oxidative phosphorylation-related metabolism, protein translation processes, and phospho-signaling modules. To assess the degree of immune infiltration in individual tumors, we identified microenvironment cell signatures that delineated four immune-based ccRCC subtypes characterized by distinct cellular pathways. This study reports a large-scale proteogenomic analysis of ccRCC to discern the functional impact of genomic alterations and provides evidence for rational treatment selection stemming from ccRCC pathobiology.

EPHB2 carried on small extracellular vesicles induces tumor angiogenesis via activation of ephrin reverse signaling.

Angiogenesis is a key process that allows nutrient uptake and cellular trafficking and is coopted in cancer to enable tumor growth and metastasis. Recently, extracellular vesicles (EVs) have been shown to promote angiogenesis; however, it is unclear what unique features EVs contribute to the process. Here, we studied the role of EVs derived from head and neck squamous cell carcinoma (HNSCC) in driving tumor angiogenesis. Small EVs (SEVs), in the size range of exosomes (50-150 nm), induced angiogenesis both in vitro and in vivo. Proteomic analysis of HNSCC SEVs revealed the cell-to-cell signaling receptor ephrin type B receptor 2 (EPHB2) as a promising candidate cargo to promote angiogenesis. Analysis of patient data further identified EPHB2 overexpression in HNSCC tumors to be associated with poor patient prognosis and tumor angiogenesis, especially in the context of overexpression of the exosome secretion regulator cortactin. Functional experiments revealed that EPHB2 expression in SEVs regulated angiogenesis both in vitro and in vivo and that EPHB2 carried by SEVs stimulates ephrin-B reverse signaling, inducing STAT3 phosphorylation. A STAT3 inhibitor greatly reduced SEV-induced angiogenesis. These data suggest a model in which EVs uniquely promote angiogenesis by transporting Eph transmembrane receptors to nonadjacent endothelial cells to induce ephrin reverse signaling.

The Epithelial to Mesenchymal Transition Promotes Glutamine Independence by Suppressing GLS2 Expression.

Identifying bioenergetics that facilitate the epithelial to mesenchymal transition (EMT) in breast cancer cells may uncover targets to treat incurable metastatic disease. Metastasis is the number one cause of cancer-related deaths; therefore, it is urgent to identify new treatment strategies to prevent the initiation of metastasis. To characterize the bioenergetics of EMT, we compared metabolic activities and gene expression in cells induced to differentiate into the mesenchymal state with their epithelial counterparts. We found that levels of GLS2, which encodes a glutaminase, are inversely associated with EMT. GLS2 down-regulation was correlated with reduced mitochondrial activity and glutamine independence even in low-glucose conditions. Restoration of GLS2 expression in GLS2-negative breast cancer cells rescued mitochondrial activity, enhanced glutamine utilization, and inhibited stem-cell properties. Additionally, inhibition of expression of the transcription factor FOXC2, a critical regulator of EMT in GLS2-negative cells, restored GLS2 expression and glutamine utilization. Furthermore, in breast cancer patients, high GLS2 expression is associated with improved survival. These findings suggest that epithelial cancer cells rely on glutamine and that cells induced to undergo EMT become glutamine independent. Moreover, the inhibition of EMT leads to a GLS2-directed metabolic shift in mesenchymal cancer cells, which may make these cells susceptible to chemotherapies.

A CTC-Cluster-Specific Signature Derived from OMICS Analysis of Patient-Derived Xenograft Tumors Predicts Outcomes in Basal-Like Breast Cancer.

Circulating tumor cell clusters (CTCcl) have a higher metastatic potential compared to single CTCs and predict long-term outcomes in breast cancer (BC) patients. Because of the rarity of CTCcls, molecular characterization of primary tumors that give rise to CTCcl hold significant promise for better diagnosis and target discovery to combat metastatic BC. In our study, we utilized the reverse-phase protein array (RPPA) and transcriptomic (RNA-Seq) data of 10 triple-negative BC patient-derived xenograft (TNBC PDX) transplantable models with CTCs and evaluated expression of upregulated candidate protein Bcl2 (B-cell lymphoma 2) by immunohistochemistry (IHC). The sample-set consisted of six CTCcl-negative (CTCcl-) and four CTCcl-positive (CTCcl+) models. We analyzed the RPPA and transcriptomic profiles of CTCcl- and CTCcl+ TNBC PDX models. In addition, we derived a CTCcl-specific gene signature for testing if it predicted outcomes using a publicly available dataset from 360 patients with basal-like BC. The RPPA analysis of CTCcl+ vs. CTCcl- TNBC PDX tumors revealed elevated expression of Bcl2 (false discovery rate (FDR) < 0.0001, fold change (FC) = 3.5) and reduced acetyl coenzyme A carboxylase-1 (ACC1) (FDR = 0.0005, FC = 0.3) in CTCcl+ compared to CTCcl- tumors. Genome-wide transcriptomic analysis of CTCcl+ vs. CTCcl- tumors revealed 549 differentially expressed genes associated with the presence of CTCcls. Apoptosis was one of the significantly downregulated pathways (normalized enrichment score (NES) = -1.69; FDR < 0.05) in TNBC PDX tumors associated with CTCcl positivity. Two out of four CTCcl+ TNBC PDX primary tumors had high Bcl2 expression by IHC (H-score > 34); whereas, only one of six CTCcl- TNBC PDX primary tumors met this criterion. Evaluation of epithelial-mesenchymal transition (EMT)-specific signature did not show significant differences between CTCcl+ and CTCcl- tumors. However, a gene signature associated with the presence of CTCcls in TNBC PDX models was associated with worse relapse-free survival in the publicly available dataset from 360 patients with basal-like BC. In summary, we identified the multigene signature of primary PDX tumors associated with the presence of CTCcls. Evaluation of additional TNBC PDX models and patients can further illuminate cellular and molecular pathways facilitating CTCcl formation.