A Biobank of Breast Cancer Explants with Preserved Intra-tumor Heterogeneity to Screen Anticancer Compounds.

The inter- and intra-tumor heterogeneity of breast cancer needs to be adequately captured in pre-clinical models. We have created a large collection of breast cancer patient-derived tumor xenografts (PDTXs), in which the morphological and molecular characteristics of the originating tumor are preserved through passaging in the mouse. An integrated platform combining in vivo maintenance of these PDTXs along with short-term cultures of PDTX-derived tumor cells (PDTCs) was optimized. Remarkably, the intra-tumor genomic clonal architecture present in the originating breast cancers was mostly preserved upon serial passaging in xenografts and in short-term cultured PDTCs. We assessed drug responses in PDTCs on a high-throughput platform and validated several ex vivo responses in vivo. The biobank represents a powerful resource for pre-clinical breast cancer pharmacogenomic studies (http://caldaslab.cruk.cam.ac.uk/bcape), including identification of biomarkers of response or resistance.

Spatial predictors of immunotherapy response in triple-negative breast cancer.

Immune checkpoint blockade (ICB) benefits some patients with triple-negative breast cancer, but what distinguishes responders from non-responders is unclear1. Because ICB targets cell-cell interactions2, we investigated the impact of multicellular spatial organization on response, and explored how ICB remodels the tumour microenvironment. We show that cell phenotype, activation state and spatial location are intimately linked, influence ICB effect and differ in sensitive versus resistant tumours early on-treatment. We used imaging mass cytometry3 to profile the in situ expression of 43 proteins in tumours from patients in a randomized trial of neoadjuvant ICB, sampled at three timepoints (baseline, n = 243; early on-treatment, n = 207; post-treatment, n = 210). Multivariate modelling showed that the fractions of proliferating CD8+TCF1+T cells and MHCII+ cancer cells were dominant predictors of response, followed by cancer-immune interactions with B cells and granzyme B+ T cells. On-treatment, responsive tumours contained abundant granzyme B+ T cells, whereas resistant tumours were characterized by CD15+ cancer cells. Response was best predicted by combining tissue features before and on-treatment, pointing to a role for early biopsies in guiding adaptive therapy. Our findings show that multicellular spatial organization is a major determinant of ICB effect and suggest that its systematic enumeration in situ could help realize precision immuno-oncology.

Dynamics of breast-cancer relapse reveal late-recurring ER-positive genomic subgroups.

The rates and routes of lethal systemic spread in breast cancer are poorly understood owing to a lack of molecularly characterized patient cohorts with long-term, detailed follow-up data. Long-term follow-up is especially important for those with oestrogen-receptor (ER)-positive breast cancers, which can recur up to two decades after initial diagnosis1-6. It is therefore essential to identify patients who have a high risk of late relapse7-9. Here we present a statistical framework that models distinct disease stages (locoregional recurrence, distant recurrence, breast-cancer-related death and death from other causes) and competing risks of mortality from breast cancer, while yielding individual risk-of-recurrence predictions. We apply this model to 3,240 patients with breast cancer, including 1,980 for whom molecular data are available, and delineate spatiotemporal patterns of relapse across different categories of molecular information (namely immunohistochemical subtypes; PAM50 subtypes, which are based on gene-expression patterns10,11; and integrative or IntClust subtypes, which are based on patterns of genomic copy-number alterations and gene expression12,13). We identify four late-recurring integrative subtypes, comprising about one quarter (26%) of tumours that are both positive for ER and negative for human epidermal growth factor receptor 2, each with characteristic tumour-driving alterations in genomic copy number and a high risk of recurrence (mean 47-62%) up to 20 years after diagnosis. We also define a subgroup of triple-negative breast cancers in which cancer rarely recurs after five years, and a separate subgroup in which patients remain at risk. Use of the integrative subtypes improves the prediction of late, distant relapse beyond what is possible with clinical covariates (nodal status, tumour size, tumour grade and immunohistochemical subtype). These findings highlight opportunities for improved patient stratification and biomarker-driven clinical trials.

Breast cancers with high proliferation and low ER-related signalling have poor prognosis and unique molecular features with implications for therapy.

BACKGROUND: Luminal breast cancers with high proliferation (MKShi) and low ER-related signalling (ERSlo) have a poor prognosis. We investigated treatment responses and molecular features of MKShi/ERSlo tumours to inform potential therapies. METHODS: Gene expression data from patients who received neoadjuvant chemotherapy (NAC) without (MDACC, N = 199) or with pembrolizumab (I-SPY2, N = 40), or endocrine therapy (NET) without (POETIC, N = 172) or with palbociclib (NeoPalAna, N = 32) were analyzed to assess treatment response by MKS/ERS-subgroups. TCGA was used to assess the mutational landscape and biomarkers associated with palbociclib-resistance (Cyclin-E, RBsig, IRPR) and immunotherapy-response (TMB, TILs, T-cell inflamed) by MKS/ERS-subgroups. RESULTS: Compared to MKShi/ERShi tumours, MKShi/ERSlo tumours had higher pathological response rates to NAC (22% vs 8%, p = 0.06) but a higher recurrence risk (4-year metastasis-free survival 70% vs 94%, p = 0.01). MKShi/ERSlo tumours frequently harboured TP53 (34%) and PIK3CA (33%) mutations, and showed high expression of Cyclin-E, RBsig and IRPR, high TMB and elevated TIL and T-cell inflamed metagene expression. MKShi/ERSlo tumours retained high proliferation after NET with or without palbociclib but had higher pathological complete response rates when pembrolizumab was added to NAC (42% vs 21%, p = 0.07). CONCLUSIONS: MKShi/ERSlo tumours have dismal outcomes and are enriched in chemotherapy-sensitive but ET- and palbociclib-resistant tumours. Biomarker analysis and clinical data suggest a potential role for immunotherapy in this group.

Mechanisms of responsiveness to and resistance against trabectedin in murine models of human myxoid liposarcoma.

Myxoid liposarcoma (MLPS) is a rare soft-tissue sarcoma characterised by the expression of FUS-DDIT3 chimera. Trabectedin has shown significant clinical anti-tumour activity against MLPS. To characterise the molecular mechanism of trabectedin sensitivity and of resistance against it, we integrated genomic and transcriptomic data from treated mice bearing ML017 or ML017/ET, two patient-derived MLPS xenograft models, sensitive to and resistant against trabectedin, respectively. Longitudinal RNA-Seq analysis of ML017 showed that trabectedin acts mainly as a transcriptional regulator: 15 days after the third dose trabectedin modulates the transcription of 4883 genes involved in processes that sustain adipocyte differentiation. No such differences were observed in ML017/ET. Genomic analysis showed that prolonged treatment causes losses in 4p15.2, 4p16.3 and 17q21.3 cytobands leading to acquired-resistance against the drug. The results dissect the complex mechanism of action of trabectedin and provide the basis for novel combinatorial approaches for the treatment of MLPS that could overcome drug-resistance.

Epithelioid Pleural Mesothelioma Is Characterized by Tertiary Lymphoid Structures in Long Survivors: Results from the MATCH Study.

Pleural mesothelioma (PM) is an aggressive tumor with few therapeutic options. Although patients with epithelioid PM (ePM) survive longer than non-epithelioid PM (non-ePM), heterogeneity of tumor response in ePM is observed. The role of the tumor immune microenvironment (TIME) in the development and progression of PM is currently considered a promising biomarker. A few studies have used high-throughput technologies correlated with TIME evaluation and morphologic and clinical data. This study aimed to identify different morphological, immunohistochemical, and transcriptional profiles that could potentially predict the outcome. A retrospective multicenter cohort of 129 chemonaive PM patients was recruited. Tissue slides were reviewed by dedicated pathologists for histotype classification and immunophenotype of tumor-infiltrating lymphocytes (TILs) and lymphoid aggregates or tertiary lymphoid structures (TLS). ePM (n = 99) survivors were further classified into long (>36 months) or short (<12 months) survivors. RNAseq was performed on a subset of 69 samples. Distinct transcriptional profiling in long and short ePM survivors was found. An inflammatory background with a higher number of B lymphocytes and a prevalence of TLS formations were detected in long compared to short ePM survivors. These results suggest that B cell infiltration could be important in modulating disease aggressiveness, opening a pathway for novel immunotherapeutic approaches.

Computational approach to discriminate human and mouse sequences in patient-derived tumour xenografts.

BACKGROUND: Patient-Derived Tumour Xenografts (PDTXs) have emerged as the pre-clinical models that best represent clinical tumour diversity and intra-tumour heterogeneity. The molecular characterization of PDTXs using High-Throughput Sequencing (HTS) is essential; however, the presence of mouse stroma is challenging for HTS data analysis. Indeed, the high homology between the two genomes results in a proportion of mouse reads being mapped as human. RESULTS: In this study we generated Whole Exome Sequencing (WES), Reduced Representation Bisulfite Sequencing (RRBS) and RNA sequencing (RNA-seq) data from samples with known mixtures of mouse and human DNA or RNA and from a cohort of human breast cancers and their derived PDTXs. We show that using an In silico Combined human-mouse Reference Genome (ICRG) for alignment discriminates between human and mouse reads with up to 99.9% accuracy and decreases the number of false positive somatic mutations caused by misalignment by >99.9%. We also derived a model to estimate the human DNA content in independent PDTX samples. For RNA-seq and RRBS data analysis, the use of the ICRG allows dissecting computationally the transcriptome and methylome of human tumour cells and mouse stroma. In a direct comparison with previously reported approaches, our method showed similar or higher accuracy while requiring significantly less computing time. CONCLUSIONS: The computational pipeline we describe here is a valuable tool for the molecular analysis of PDTXs as well as any other mixture of DNA or RNA species.

Diagnostic role of circulating extracellular matrix-related proteins in non-small cell lung cancer.

BACKGROUND: Interactions between cancer cells and the surrounding microenvironment are crucial determinants of cancer progression. During this process, bi-directional communication among tumor cells and cancer associated fibroblasts (CAF) regulate extracellular matrix (ECM) deposition and remodeling. As a result of this dynamic process, soluble ECM proteins can be released into the bloodstream and may represent novel circulating biomarkers useful for cancer diagnosis. The aim of the present study was to measure the levels of three circulating ECM related proteins (COL11A1, COL10A1 and SPARC) in plasma samples of lung cancer patients and in healthy heavy-smokers controls and test whether such measurements have diagnostic or prognostic value. METHODS: Gene expression profiling of lung fibroblasts isolated from paired normal and cancer tissue of NSCLC patients was performed by gene expression microarrays. The prioritization of the candidates for the study of circulating proteins in plasma was based on the most differentially expressed genes in cancer associated fibroblasts. Soluble ECM proteins were assessed by western blot in the conditioned medium of lung fibroblasts and by ELISA assays in plasma samples. RESULTS: Plasma samples from lung cancer patients and healthy heavy-smokers controls were tested for levels of COL11A1 and COL10A1 (n = 57 each) and SPARC (n = 90 each). Higher plasma levels of COL10A1 were detected in patients (p ≤ 0.001), a difference that was driven specifically by females (p < 0.001). No difference in COL11A1 levels between patients and controls was found. SPARC levels were also higher in plasma patients than controls (p < 0.001) with good performance in discriminating the two groups (AUC = 0.744). No significant association was observed between plasma proteins levels and clinicopathological features or survival. CONCLUSION: Soluble factors related to proficient tumor-stroma cross-talk are detectable in plasma of primary lung cancer patients and may represent a valuable complementary diagnostic tool to discriminate lung cancer patients from healthy heavy-smokers individuals as shown for the SPARC protein.

Dissecting Time- from Tumor-Related Gene Expression Variability in Bilateral Breast Cancer.

Metachronous (MBC) and synchronous bilateral breast tumors (SBC) are mostly distinct primaries, whereas paired primaries and their local recurrences (LRC) share a common origin. Intra-pair gene expression variability in MBC, SBC, and LRC derives from time/tumor microenvironment-related and tumor genetic background-related factors and pairs represents an ideal model for trying to dissect tumor-related from microenvironment-related variability. Pairs of tumors derived from women with SBC (n = 18), MBC (n = 11), and LRC (n = 10) undergoing local-regional treatment were profiled for gene expression; similarity between pairs was measured using an intraclass correlation coefficient (ICC) computed for each gene and compared using analysis of variance (ANOVA). When considering biologically unselected genes, the highest correlations were found for primaries and paired LRC, and the lowest for MBC pairs. By instead limiting the analysis to the breast cancer intrinsic genes, correlations between primaries and paired LRC were enhanced, while lower similarities were observed for SBC and MBC. Focusing on stromal-related genes, the ICC values decreased for MBC and were significantly different from SBC. These findings indicate that it is possible to dissect intra-pair gene expression variability into components that are associated with genetic origin or with time and microenvironment by using specific gene subsets.

Complexity in the tumour microenvironment: Cancer associated fibroblast gene expression patterns identify both common and unique features of tumour-stroma crosstalk across cancer types.

Cancer is a complex disease, driven by the accumulation of several somatic aberrations but fostered by a two-way interaction between tumour cells and the surrounding microenvironment. Cancer associated fibroblasts (CAFs) represent one of the major players in tumour-stroma crosstalk. Recent in vitro and in vivo studies, often conducted by employing high throughput approaches, have started unravelling the key pathways involved in their functional effects. This review focus on open challenges in the study of CAF properties and function, highlighting at the same time the existence of common mechanisms as well as peculiarities in different cancer types (breast, prostate and lung cancer). Although still limited by current experimental models, which are unable to deal with the full level of complexity of the tumour microenvironment, a better understanding of these mechanisms may enable the identification of new biomarkers and therapeutic targets, to improve current strategies for cancer diagnosis and treatment.

Proposal of supervised data analysis strategy of plasma miRNAs from hybridisation array data with an application to assess hemolysis-related deregulation.

BACKGROUND: Plasma miRNAs have the potential as cancer biomarkers but no consolidated guidelines for data mining in this field are available. The purpose of the study was to apply a supervised data analysis strategy in a context where prior knowledge is available, i.e., that of hemolysis-related miRNAs deregulation, so as to compare our results with existing evidence. RESULTS: We developed a structured strategy with innovative applications of existing bioinformatics methods for supervised analyses including: 1) the combination of two statistical (t- and Anderson-Darling) test results to detect miRNAs with significant fold change or general distributional differences in class comparison, which could reveal hidden differential biological processes worth to be considered for building predictive tools; 2) a bootstrap selection procedure together with machine learning techniques in class prediction to guarantee the transferability of results and explore the interconnections among the selected miRNAs, which is important for highlighting their inherent biological dependences. The strategy was applied to develop a classifier for discriminating between hemolyzed and not hemolyzed plasma samples, defined according to a recently published hemolysis score. We identified five miRNAs with increased expression in hemolyzed plasma samples (miR-486-5p, miR-92a, miR-451, miR-16, miR-22). CONCLUSIONS: We identified four miRNAs previously reported in the literature as hemolysis related together with a new one (miR-22).which needs further investigations. Our findings confirm the validity of the proposed strategy and, in parallel, the hemolysis score capability to be used as pre-analytic hemolysis detector. R codes for implementing the approaches are provided.

Gene expression profiling of circulating tumor cells in breast cancer.

BACKGROUND: Determining the transcriptional profile of circulating tumor cells (CTCs) may allow the acquisition of clinically relevant information while overcoming tumor heterogeneity-related biases associated with use of tissue samples for biomarker assessment. However, such molecular characterization is challenging because CTCs are rare and outnumbered by blood cells. METHODS: Here, we describe a technical protocol to measure the expression of >29 000 genes in CTCs captured from whole blood with magnetic beads linked with antibodies against epithelial cell adhesion molecule (EpCAM) and the carcinoma-associated mucin, MUC1, designed to be used for CTC characterization in clinical samples. Low numbers of cells (5-200) from the MCF7 and MDA-MB-468 breast cancer cell lines were spiked in healthy donor blood samples and isolated with the AdnaTest EMT-1/Stem CellSelect kit. Gene expression profiles (GEPs) were obtained with the WG-DASL HT assay and compared with GEPs obtained from RNA isolated from cultured cell lines and unspiked samples. RESULTS: GEPs from samples containing 25 or more spiked cells correlated (r = 0.95) with cognate 100-ng RNA input samples, clustered separately from blood control samples, and allowed MCF7 and MDA-MB-468 cells to be distinguished. GEPs with comparable technical quality were also obtained in a preliminary series of clinical samples. CONCLUSIONS: Our approach allows technically reliable GEPs to be obtained from isolated CTCs for the acquisition of biologically useful information. It is reproducible and suitable for application in prospective studies to assess the clinical utility of CTC GEPs, provided that >25 CTCs can be isolated.

Bioinformatics tools for secretome analysis.

Over recent years, analyses of secretomes (complete sets of secreted proteins) have been reported in various organisms, cell types, and pathologies and such studies are quickly gaining popularity. Fungi secrete enzymes can break down potential food sources; plant secreted proteins are primarily parts of the cell wall proteome; and human secreted proteins are involved in cellular immunity and communication, and provide useful information for the discovery of novel biomarkers, such as for cancer diagnosis. Continuous development of methodologies supports the wide identification and quantification of secreted proteins in a given cellular state. The role of secreted factors is also investigated in the context of the regulation of major signaling events, and connectivity maps are built to describe the differential expression and dynamic changes of secretomes. Bioinformatics has become the bridge between secretome data and computational tasks for managing, mining, and retrieving information. Predictions can be made based on this information, contributing to the elucidation of a given organism's physiological state and the determination of the specific malfunction in disease states. Here we provide an overview of the available bioinformatics databases and software that are used to analyze the biological meaning of secretome data, including descriptions of the main functions and limitations of these tools. The important challenges of data analysis are mainly related to the integration of biological information from dissimilar sources. Improvements in databases and developments in software will likely substantially contribute to the usefulness and reliability of secretome studies. This article is part of a Special Issue entitled: An Updated Secretome.

Mesenchymal to amoeboid transition is associated with stem-like features of melanoma cells.

BACKGROUND: Cellular plasticity confers cancer cells the ability to adapt to microenvironmental changes, a fundamental requirement for tumour progression and metastasis. The epithelial to mesenchymal transition (EMT) is a transcriptional programme associated with increased cell motility and stemness. Besides EMT, the mesenchymal to amoeboid transition (MAT) has been described during tumour progression but to date, little is known about its transcriptional control and involvement in stemness. The aim of this manuscript is to investigate (i) the transcriptional profile associated with the MAT programme and (ii) to study whether MAT acquisition in melanoma cancer cells correlates with clonogenic potential to promote tumour growth. RESULTS: By using a multidisciplinary approach, we identified four different treatments able to induce MAT in melanoma cells: EphA2 overexpression, Rac1 functional inhibition using its RacN17 dominant negative mutant, stimulation with Ilomastat or treatment with the RhoA activator Calpeptin. First, gene expression profiling identified the transcriptional pathways associated with MAT, independently of the stimulus that induces the MAT programme. Notably, gene sets associated with the repression of mesenchymal traits, decrease in the secretion of extracellular matrix components as well as increase of cellular stemness positively correlate with MAT. Second, the link between MAT and stemness has been investigated in vitro by analysing stemness markers and clonogenic potential of melanoma cells undergoing MAT. Finally, the link between MAT inducing treatments and tumour initiating capability has been validated in vivo. CONCLUSION: Taken together, our results demonstrate that MAT programme in melanoma is characterised by increased stemness and clonogenic features of cancer cells, thus sustaining tumour progression. Furthermore, these data suggest that stemness is not an exclusive feature of cells undergoing EMT, but more generally is associated with an increase in cellular plasticity of cancer cells.

Unlocking the potential of Molecular Tumor Boards: from cutting-edge data interpretation to innovative clinical pathways.

The emerging era of precision medicine is characterized by an increasing availability of targeted anticancer therapies and by the parallel development of techniques to obtain more refined molecular data, whose interpretation may not always be straightforward. Molecular tumor boards gather various professional figures, in order to leverage the analysis of molecular data and provide prognostic and predictive insights for clinicians. In addition to healthcare development, they could also become a tool to promote knowledge and research spreading. A growing body of evidence on the application of molecular tumor boards to clinical practice is forming and positive signals are emerging, although a certain degree of heterogeneity exists. This work analyzes molecular tumor boards' potential workflows, figures involved, data sources, sample matrices and eligible patients, as well as available evidence and learning examples. The emerging concept of multi-institutional, disease-specific molecular tumor boards is also considered by presenting two ongoing nationwide experiences.

Converging and evolving immuno-genomic routes toward immune escape in breast cancer.

The interactions between tumor and immune cells along the course of breast cancer progression remain largely unknown. Here, we extensively characterize multiple sequential and parallel multiregion tumor and blood specimens of an index patient and a cohort of metastatic triple-negative breast cancers. We demonstrate that a continuous increase in tumor genomic heterogeneity and distinct molecular clocks correlated with resistance to treatment, eventually allowing tumors to escape from immune control. TCR repertoire loses diversity over time, leading to convergent evolution as breast cancer progresses. Although mixed populations of effector memory and cytotoxic single T cells coexist in the peripheral blood, defects in the antigen presentation machinery coupled with subdued T cell recruitment into metastases are observed, indicating a potent immune avoidance microenvironment not compatible with an effective antitumor response in lethal metastatic disease. Our results demonstrate that the immune responses against cancer are not static, but rather follow dynamic processes that match cancer genomic progression, illustrating the complex nature of tumor and immune cell interactions.

Feasibility of circulating miRNA microarray analysis from archival plasma samples.

MicroRNAs have been found to be deregulated in several diseases and, due to their high stability in body fluids, represent promising noninvasively detectable biomarkers. However, numerous technical variables can affect accurate measurement of circulating miRNAs. Using a microarray-based method we assessed the: (i) adequate intra- and inter-array reproducibility of miRNA profiling; (ii) feasibility of using archival plasma samples stored for an extended period of time and available in limited amounts; (iii) good correlation between different batches; and (iv) time-dependent increase of background signals close to the chip expiration date.

Cancer-specific association between Tau (MAPT) and cellular pathways, clinical outcome, and drug response.

Tau (MAPT) is a microtubule-associated protein causing common neurodegenerative diseases or rare inherited frontotemporal lobar degenerations. Emerging evidence for non-canonical functions of Tau in DNA repair and P53 regulation suggests its involvement in cancer. To bring new evidence for a relevant role of Tau in cancer, we carried out an in-silico pan-cancer analysis of MAPT transcriptomic profile in over 10000 clinical samples from 32 cancer types and over 1300 pre-clinical samples from 28 cancer types provided by the TCGA and the DEPMAP datasets respectively. MAPT expression associated with key cancer hallmarks including inflammation, proliferation, and epithelial to mesenchymal transition, showing cancer-specific patterns. In some cancer types, MAPT functional networks were affected by P53 mutational status. We identified new associations of MAPT with clinical outcomes and drug response in a context-specific manner. Overall, our findings indicate that the MAPT gene is a potential major player in multiple types of cancer. Importantly, the impact of Tau on cancer seems to be heavily influenced by the specific cellular environment.

Pan-cancer analysis of antibody-drug conjugate targets and putative predictors of treatment response.

BACKGROUND: Antibody-drug conjugates (ADCs) are a rapidly expanding class of compounds in oncology. Our goal was to assess the expression of ADC targets and potential downstream determining factors of activity across pan-cancer and normal tissues. MATERIALS AND METHODS: ADCs in clinical trials (n = 121) were identified through ClinicalTrials.gov, corresponding to 54 targets. Genes potentially implicated in treatment response were identified in the literature. Gene expression from The Cancer Genome Atlas (9000+ cancers of 31 cancer types), the Genotype-Tissue Expression database (n = 19,000 samples from 31 normal tissue types), and the TNMplot.com (n = 12,494 unmatched primary and metastatic samples) were used in this analysis. To compare relative expression across and within tumour types we used pooled normal tissues as reference. RESULTS: For most ADC targets, mRNA levels correlated with protein expression. Pan-cancer target expression distributions identified appealing cancer types for each ADC development. Co-expression of multiple targets was common and suggested opportunities for ADC combinations. Expression levels of genes potentially implicated in ADC response downstream of the target might provide additional information (e.g. TOP1 was highly expressed in many tumour types, including breast and lung cancers). Metastatic compared to primary tissues overexpressed some ADCs targets. Single sample "targetgram" plots were generated to visualise the expression of potentially competing ADC targets and resistance/sensitivity markers highlighting high inter-patient heterogeneity. Off-cancer target expression only partially explains adverse events, while expression of determinants of payload activity explained more of the observed toxicities. CONCLUSION: Our findings draw attention to new therapeutic opportunities for ADCs that can be tested in the clinic and our web platform (https://tnmplot.com) can assist in prioritising upcoming ADC targets for clinical development.