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.

The single-cell pathology landscape of breast cancer.

Single-cell analyses have revealed extensive heterogeneity between and within human tumours1-4, but complex single-cell phenotypes and their spatial context are not at present reflected in the histological stratification that is the foundation of many clinical decisions. Here we use imaging mass cytometry5 to simultaneously quantify 35 biomarkers, resulting in 720 high-dimensional pathology images of tumour tissue from 352 patients with breast cancer, with long-term survival data available for 281 patients. Spatially resolved, single-cell analysis identified the phenotypes of tumour and stromal single cells, their organization and their heterogeneity, and enabled the cellular architecture of breast cancer tissue to be characterized on the basis of cellular composition and tissue organization. Our analysis reveals multicellular features of the tumour microenvironment and novel subgroups of breast cancer that are associated with distinct clinical outcomes. Thus, spatially resolved, single-cell analysis can characterize intratumour phenotypic heterogeneity in a disease-relevant manner, with the potential to inform patient-specific diagnosis.

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.

Combined quantitative measures of ER, PR, HER2, and KI67 provide more prognostic information than categorical combinations in luminal breast cancer.

Although most women with luminal breast cancer do well on endocrine therapy alone, some will develop fatal recurrence thereby necessitating the need to prospectively determine those for whom additional cytotoxic therapy will be beneficial. Categorical combinations of immunohistochemical measures of ER, PR, HER2, and KI67 are traditionally used to classify patients into luminal A-like and B-like subtypes for chemotherapeutic reasons, but this may lead to the loss of prognostically relevant information. Here, we compared the prognostic value of quantitative measures of these markers, combined in the IHC4-score, to categorical combinations in subtypes. Using image analysis-based scores for all four markers, we computed the IHC4-score for 2498 patients with luminal breast cancer from two European study populations. We defined subtypes (A-like (ER + and PR + : and HER2- and low KI67) and B-like (ER + and/or PR + : and HER2 + or high KI67)) by combining binary categories of these markers. Hazard ratios and 95% confidence intervals for associations with 10-year breast cancer-specific survival were estimated in Cox proportional-hazard models. We accounted for clinical prognostic factors, including grade, tumor size, lymph-nodal involvement, and age, by using the PREDICT-score. Overall, Subtypes [hazard ratio (95% confidence interval) B-like vs. A-like = 1.64 (1.25-2.14); P-value < 0.001] and IHC4-score [hazard ratio (95% confidence interval)/1 standard deviation = 1.32 (1.20-1.44); P-value < 0.001] were prognostic in univariable models. However, IHC4-score [hazard ratio (95% confidence interval)/1 standard deviation = 1.24 (1.11-1.37); P-value < 0.001; likelihood ratio chi-square (LRχ2) = 12.5] provided more prognostic information than Subtype [hazard ratio (95% confidence interval) B-like vs. A-like = 1.38 (1.02-1.88); P-value = 0.04; LRχ2 = 4.3] in multivariable models. Further, higher values of the IHC4-score were associated with worse prognosis, regardless of subtype (P-heterogeneity = 0.97). These findings enhance the value of the IHC4-score as an adjunct to clinical prognostication tools for aiding chemotherapy decision-making in luminal breast cancer patients, irrespective of subtype.

Etiology of hormone receptor positive breast cancer differs by levels of histologic grade and proliferation.

Limited epidemiological evidence suggests that the etiology of hormone receptor positive (HR+) breast cancer may differ by levels of histologic grade and proliferation. We pooled risk factor and pathology data on 5,905 HR+ breast cancer cases and 26,281 controls from 11 epidemiological studies. Proliferation was determined by centralized automated measures of KI67 in tissue microarrays. Odds ratios (OR), 95% confidence intervals (CI) and p-values for case-case and case-control comparisons for risk factors in relation to levels of grade and quartiles (Q1-Q4) of KI67 were estimated using polytomous logistic regression models. Case-case comparisons showed associations between nulliparity and high KI67 [OR (95% CI) for Q4 vs. Q1 = 1.54 (1.22, 1.95)]; obesity and high grade [grade 3 vs. 1 = 1.68 (1.31, 2.16)] and current use of combined hormone therapy (HT) and low grade [grade 3 vs. 1 = 0.27 (0.16, 0.44)] tumors. In case-control comparisons, nulliparity was associated with elevated risk of tumors with high but not low levels of proliferation [1.43 (1.14, 1.81) for KI67 Q4 vs. 0.83 (0.60, 1.14) for KI67 Q1]; obesity among women ≥50 years with high but not low grade tumors [1.55 (1.17, 2.06) for grade 3 vs. 0.88 (0.66, 1.16) for grade 1] and HT with low but not high grade tumors [3.07 (2.22, 4.23) for grade 1 vs. 0.85 (0.55, 1.30) for grade 3]. Menarcheal age and family history were similarly associated with HR+ tumors of different grade or KI67 levels. These findings provide insights into the etiologic heterogeneity of HR+ tumors.

Differential oestrogen receptor binding is associated with clinical outcome in breast cancer.

Oestrogen receptor-α (ER) is the defining and driving transcription factor in the majority of breast cancers and its target genes dictate cell growth and endocrine response, yet genomic understanding of ER function has been restricted to model systems. Here we map genome-wide ER-binding events, by chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq), in primary breast cancers from patients with different clinical outcomes and in distant ER-positive metastases. We find that drug-resistant cancers still recruit ER to the chromatin, but that ER binding is a dynamic process, with the acquisition of unique ER-binding regions in tumours from patients that are likely to relapse. The acquired ER regulatory regions associated with poor clinical outcome observed in primary tumours reveal gene signatures that predict clinical outcome in ER-positive disease exclusively. We find that the differential ER-binding programme observed in tumours from patients with poor outcome is not due to the selection of a rare subpopulation of cells, but is due to the FOXA1-mediated reprogramming of ER binding on a rapid timescale. The parallel redistribution of ER and FOXA1 binding events in drug-resistant cellular contexts is supported by histological co-expression of ER and FOXA1 in metastatic samples. By establishing transcription-factor mapping in primary tumour material, we show that there is plasticity in ER-binding capacity, with distinct combinations of cis-regulatory elements linked with the different clinical outcomes.

Association between tumour infiltrating lymphocytes, histotype and clinical outcome in epithelial ovarian cancer.

BACKGROUND: There is evidence that some ovarian tumours evoke an immune response, which can be assessed by tumour infiltrating lymphocytes (TILs). To facilitate adoption of TILs as a clinical biomarker, a standardised method for their H&E visual evaluation has been validated in breast cancer. METHODS: We sought to investigate the prognostic significance of TILs in a study of 953 invasive epithelial ovarian cancer tumour samples, both primary and metastatic, from 707 patients from the prospective population-based SEARCH study. TILs were analysed using a standardised method based on H&E staining producing a percentage score for stromal and intratumoral compartments. We used Cox regression to estimate hazard ratios of the association between TILs and survival. RESULTS: The extent of stromal and intra-tumoral TILs were correlated in the primary tumours (n = 679, Spearman's rank correlation = 0.60, P < 0.001) with a similar correlation in secondary tumours (n = 224, Spearman's rank correlation = 0.62, P < 0.001). There was a weak correlation between stromal TIL levels in primary and secondary tumour samples (Spearman's rank correlation = 0.29, P < 0.001) and intra-tumoral TIL levels in primary and secondary tumour samples (Spearman's rank correlation = 0.19, P = 0.0094). The extent of stromal TILs differed between histotypes (Pearson chi2 (12d.f.) 54.1, P < 0.0001) with higher levels of stromal infiltration in the high-grade serous and endometriod cases. A significant association was observed for higher intratumoral TIL levels and a favourable prognosis (HR 0.74 95% CI 0.55-1.00 p = 0.047). CONCLUSION: This study is the largest collection of epithelial ovarian tumour samples evaluated for TILs. We have shown that stromal and intratumoral TIL levels are correlated and that their levels correlate with clinical variables such as tumour histological subtype. We have also shown that increased levels of both intratumoral and stromal TILs are associated with a better prognosis; however, this is only statistically significant for intratumoral TILs. This study suggests that a clinically useful immune prognostic indicator in epithelial ovarian cancer could be developed using this technique.

Computational pathology of pre-treatment biopsies identifies lymphocyte density as a predictor of response to neoadjuvant chemotherapy in breast cancer.

BACKGROUND: There is a need to improve prediction of response to chemotherapy in breast cancer in order to improve clinical management and this may be achieved by harnessing computational metrics of tissue pathology. We investigated the association between quantitative image metrics derived from computational analysis of digital pathology slides and response to chemotherapy in women with breast cancer who received neoadjuvant chemotherapy. METHODS: We digitised tissue sections of both diagnostic and surgical samples of breast tumours from 768 patients enrolled in the Neo-tAnGo randomized controlled trial. We subjected digital images to systematic analysis optimised for detection of single cells. Machine-learning methods were used to classify cells as cancer, stromal or lymphocyte and we computed estimates of absolute numbers, relative fractions and cell densities using these data. Pathological complete response (pCR), a histological indicator of chemotherapy response, was the primary endpoint. Fifteen image metrics were tested for their association with pCR using univariate and multivariate logistic regression. RESULTS: Median lymphocyte density proved most strongly associated with pCR on univariate analysis (OR 4.46, 95 % CI 2.34-8.50, p < 0.0001; observations = 614) and on multivariate analysis (OR 2.42, 95 % CI 1.08-5.40, p = 0.03; observations = 406) after adjustment for clinical factors. Further exploratory analyses revealed that in approximately one quarter of cases there was an increase in lymphocyte density in the tumour removed at surgery compared to diagnostic biopsies. A reduction in lymphocyte density at surgery was strongly associated with pCR (OR 0.28, 95 % CI 0.17-0.47, p < 0.0001; observations = 553). CONCLUSIONS: A data-driven analysis of computational pathology reveals lymphocyte density as an independent predictor of pCR. Paradoxically an increase in lymphocyte density, following exposure to chemotherapy, is associated with a lack of pCR. Computational pathology can provide objective, quantitative and reproducible tissue metrics and represents a viable means of outcome prediction in breast cancer. TRIAL REGISTRATION: ClinicalTrials.gov NCT00070278 ; 03/10/2003.

Prognostic value of automated KI67 scoring in breast cancer: a centralised evaluation of 8088 patients from 10 study groups.

BACKGROUND: The value of KI67 in breast cancer prognostication has been questioned due to concerns on the analytical validity of visual KI67 assessment and methodological limitations of published studies. Here, we investigate the prognostic value of automated KI67 scoring in a large, multicentre study, and compare this with pathologists' visual scores available in a subset of patients. METHODS: We utilised 143 tissue microarrays containing 15,313 tumour tissue cores from 8088 breast cancer patients in 10 collaborating studies. A total of 1401 deaths occurred during a median follow-up of 7.5 years. Centralised KI67 assessment was performed using an automated scoring protocol. The relationship of KI67 levels with 10-year breast cancer specific survival (BCSS) was investigated using Kaplan-Meier survival curves and Cox proportional hazard regression models adjusted for known prognostic factors. RESULTS: Patients in the highest quartile of KI67 (>12 % positive KI67 cells) had a worse 10-year BCSS than patients in the lower three quartiles. This association was statistically significant for ER-positive patients (hazard ratio (HR) (95 % CI) at baseline = 1.96 (1.31-2.93); P = 0.001) but not for ER-negative patients (1.23 (0.86-1.77); P = 0.248) (P-heterogeneity = 0.064). In spite of differences in characteristics of the study populations, the estimates of HR were consistent across all studies (P-heterogeneity = 0.941 for ER-positive and P-heterogeneity = 0.866 for ER-negative). Among ER-positive cancers, KI67 was associated with worse prognosis in both node-negative (2.47 (1.16-5.27)) and node-positive (1.74 (1.05-2.86)) tumours (P-heterogeneity = 0.671). Further classification according to ER, PR and HER2 showed statistically significant associations with prognosis among hormone receptor-positive patients regardless of HER2 status (P-heterogeneity = 0.270) and among triple-negative patients (1.70 (1.02-2.84)). Model fit parameters were similar for visual and automated measures of KI67 in a subset of 2440 patients with information from both sources. CONCLUSIONS: Findings from this large-scale multicentre analysis with centrally generated automated KI67 scores show strong evidence in support of a prognostic value for automated KI67 scoring in breast cancer. Given the advantages of automated scoring in terms of its potential for standardisation, reproducibility and throughput, automated methods appear to be promising alternatives to visual scoring for KI67 assessment.

Patterns of Immune Infiltration in Breast Cancer and Their Clinical Implications: A Gene-Expression-Based Retrospective Study.

BACKGROUND: Immune infiltration of breast tumours is associated with clinical outcome. However, past work has not accounted for the diversity of functionally distinct cell types that make up the immune response. The aim of this study was to determine whether differences in the cellular composition of the immune infiltrate in breast tumours influence survival and treatment response, and whether these effects differ by molecular subtype. METHODS AND FINDINGS: We applied an established computational approach (CIBERSORT) to bulk gene expression profiles of almost 11,000 tumours to infer the proportions of 22 subsets of immune cells. We investigated associations between each cell type and survival and response to chemotherapy, modelling cellular proportions as quartiles. We found that tumours with little or no immune infiltration were associated with different survival patterns according to oestrogen receptor (ER) status. In ER-negative disease, tumours lacking immune infiltration were associated with the poorest prognosis, whereas in ER-positive disease, they were associated with intermediate prognosis. Of the cell subsets investigated, T regulatory cells and M0 and M2 macrophages emerged as the most strongly associated with poor outcome, regardless of ER status. Among ER-negative tumours, CD8+ T cells (hazard ratio [HR] = 0.89, 95% CI 0.80-0.98; p = 0.02) and activated memory T cells (HR 0.88, 95% CI 0.80-0.97; p = 0.01) were associated with favourable outcome. T follicular helper cells (odds ratio [OR] = 1.34, 95% CI 1.14-1.57; p < 0.001) and memory B cells (OR = 1.18, 95% CI 1.0-1.39; p = 0.04) were associated with pathological complete response to neoadjuvant chemotherapy in ER-negative disease, suggesting a role for humoral immunity in mediating response to cytotoxic therapy. Unsupervised clustering analysis using immune cell proportions revealed eight subgroups of tumours, largely defined by the balance between M0, M1, and M2 macrophages, with distinct survival patterns by ER status and associations with patient age at diagnosis. The main limitations of this study are the use of diverse platforms for measuring gene expression, including some not previously used with CIBERSORT, and the combined analysis of different forms of follow-up across studies. CONCLUSIONS: Large differences in the cellular composition of the immune infiltrate in breast tumours appear to exist, and these differences are likely to be important determinants of both prognosis and response to treatment. In particular, macrophages emerge as a possible target for novel therapies. Detailed analysis of the cellular immune response in tumours has the potential to enhance clinical prediction and to identify candidates for immunotherapy.

Heterogeneity of luminal breast cancer characterised by immunohistochemical expression of basal markers.

BACKGROUND: Luminal A breast cancer defined as hormone receptor positive and human epidermal growth factor receptor 2 (HER2) negative is known to be heterogeneous. Previous study showed that luminal A tumours with the expression of basal markers ((cytokeratin (CK) 5 or CK5/6) or epidermal growth factor receptor (EGFR)) were associated with poorer prognosis compared with those that stained negative for basal markers. Prompted by this study, we assessed whether tumour characteristics and risk factors differed by basal marker status within luminal A tumours. METHODS: We pooled 5040 luminal A cases defined by immunohistochemistry (4490 basal-negative ((CK5 (or CK5/6))- and EGFR-) and 550 basal-positive ((CK5 (or CK5/6+)) or EGFR+)) from eight studies participating in the Breast Cancer Association Consortium. Case-case comparison was performed using unconditional logistic regression. RESULTS: Tumour characteristics and risk factors did not vary significantly by the expression of basal markers, although results suggested that basal-positive luminal tumours tended to be smaller and node negative, and were more common in women with a positive family history and lower body mass index. CONCLUSIONS: Most established breast cancer risk factors were similar in basal-positive and basal-negative luminal A tumours. The non-significant but suggestive differences in tumour features and family history warrant further investigations.

Spatial Biology of Breast Cancer.

Spatial findings have shaped on our understanding of breast cancer. In this review, we discuss how spatial methods, including spatial transcriptomics and proteomics and the resultant understanding of spatial relationships, have contributed to concepts regarding cancer progression and treatment. In addition to discussing traditional approaches, we examine how emerging multiplex imaging technologies have contributed to the field and how they might influence future research.

Biological and prognostic associations of miR-205 and let-7b in breast cancer revealed by in situ hybridization analysis of micro-RNA expression in arrays of archival tumour tissue.

Micro-RNAs (miRNAs) are frequently dysregulated in a range of human malignancies, many have been shown to act either as tumour supressors or oncogenes and several have been implicated in breast cancer. However, breast cancer is a diverse disease and little is known about the relationships between miRNA expression, clinical outcome and tumour subtype. We used locked nucleic acid probe in situ hybridization (LNA-ISH) to visualize, in tissue micro-arrays (TMAs) of 2919 formalin-fixed paraffin-embedded (FFPE) archival breast tumours, the expression of two key miRNAs that are frequently lost in a range of solid malignancies, let-7b and miR-205. These miRNAs were also quantified by quantitative reverse transcription PCR in cores of FFPE tissue from 40 of these cases, demonstrating that LNA-ISH is semi-quantitative. The tumours in the TMAs were assigned to subtypes based on their immunohistochemical (IHC) staining with ER, PR, HER2, CK5/6 and EGFR. let-7b expression was shown to be associated with luminal tumours and to have an independent significant positive prognostic value in this group. miR-205 is associated with tumours of ductal morphology and is of significant positive prognostic value within these tumours. We propose that the expression of miR-205 may contribute to ductal tumour morphology.

A Ki67/BCL2 index based on immunohistochemistry is highly prognostic in ER-positive breast cancer.

There is an urgent need to improve prognostic classifiers in breast cancer. Ki67 and B-cell lymphoma 2 protein (BCL2) are established prognostic markers which have traditionally been assessed separately, in a dichotomous manner. This study was conducted to test the hypothesis that combinatorial assessment of these markers would provide superior prognostic information and improve their clinical utility. Tissue microarrays were used to assess the expression of Ki67 and BCL2 in 2749 cases of invasive breast cancer. We devised a Ki67/BCL2 index representing the relative expression of each protein and assessed its association with breast cancer-specific survival (BCSS) using a Cox proportional-hazards model. Based on our findings, an independent cohort of 3992 cases was used to validate the prognostic significance of the Ki67/BCL2 index. All survival analyses were conducted on complete data as well as data where missing values were resolved using multiple imputation. This study complied with reporting recommendations for tumour marker prognostic studies (REMARK) criteria. The Ki67/BCL2 index showed a significant association with BCSS at 10 years in estrogen receptor (ER)-positive disease. In multivariate analysis, adjusting for major clinical and molecular markers, the Ki67/BCL2 index retained prognostic significance, robustly classifying cases into three risk groups [intermediate- versus low-risk hazard ratio (HR), 1.5; 95% confidence interval (95% CI), 1.0-2.0; p = 0.031; high- versus low-risk HR, 2.6; 95% CI, 1.3-5.0; p = 0.005]. This finding was validated in an independent cohort of 3992 tumours containing 2761 ER-positive tumours (intermediate- versus low-risk HR, 1.7; 95% CI, 1.3-2.1; p < 0.001; high- versus low-risk HR, 2.0; 95% CI, 1.4-2.9; p < 0.001). Ki67 and BCL2 can be effectively combined to produce an index which is an independent predictor of BCSS in ER-positive breast cancer, enhancing their potential prognostic utility.

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.

Breast tumor microenvironment structures are associated with genomic features and clinical outcome.

The functions of the tumor microenvironment (TME) are orchestrated by precise spatial organization of specialized cells, yet little is known about the multicellular structures that form within the TME. Here we systematically mapped TME structures in situ using imaging mass cytometry and multitiered spatial analysis of 693 breast tumors linked to genomic and clinical data. We identified ten recurrent TME structures that varied by vascular content, stromal quiescence versus activation, and leukocyte composition. These TME structures had distinct enrichment patterns among breast cancer subtypes, and some were associated with genomic profiles indicative of immune escape. Regulatory and dysfunctional T cells co-occurred in large 'suppressed expansion' structures. These structures were characterized by high cellular diversity, proliferating cells and enrichment for BRCA1 and CASP8 mutations and predicted poor outcome in estrogen-receptor-positive disease. The multicellular structures revealed here link conserved spatial organization to local TME function and could improve patient stratification.

Cancer stem cell markers in breast cancer: pathological, clinical and prognostic significance.

INTRODUCTION: The cancer stem cell (CSC) hypothesis states that tumours consist of a cellular hierarchy with CSCs at the apex driving tumour recurrence and metastasis. Hence, CSCs are potentially of profound clinical importance. We set out to establish the clinical relevance of breast CSC markers by profiling a large cohort of breast tumours in tissue microarrays (TMAs) using immunohistochemistry (IHC). METHODS: We included 4, 125 patients enrolled in the SEARCH population-based study with tumours represented in TMAs and classified into molecular subtype according to a validated IHC-based five-marker scheme. IHC was used to detect CD44/CD24, ALDH1A1, aldehyde dehydrogenase family 1 member A3 (ALDH1A3) and integrin alpha-6 (ITGA6). A 'Total CSC' score representing expression of all four CSC markers was also investigated. Association with breast cancer specific survival (BCSS) at 10 years was assessed using a Cox proportional-hazards model. This study was complied with REMARK criteria. RESULTS: In ER negative cases, multivariate analysis showed that ITGA6 was an independent prognostic factor with a time-dependent effect restricted to the first two years of follow-up (hazard ratio (HR) for 0 to 2 years follow-up, 2.4; 95% confidence interval (95% CI), 1.2 to 4.8; P = 0.009). The composite 'Total CSC' score carried independent prognostic significance in ER negative cases for the first four years of follow-up (HR for 0 to 4 years follow-up, 1.3; 95% CI, 1.1 to 1.6; P = 0.006). CONCLUSIONS: Breast CSC markers do not identify identical subpopulations in primary tumours. Both ITGA6 and a composite Total CSC score show independent prognostic significance in ER negative disease. The use of multiple markers to identify tumours enriched for CSCs has the greatest prognostic value. In the absence of more specific markers, we propose that the effective translation of the CSC hypothesis into patient benefit will necessitate the use of a panel of markers to robustly identify tumours enriched for CSCs.

Landscapes of cellular phenotypic diversity in breast cancer xenografts and their impact on drug response.

The heterogeneity of breast cancer plays a major role in drug response and resistance and has been extensively characterized at the genomic level. Here, a single-cell breast cancer mass cytometry (BCMC) panel is optimized to identify cell phenotypes and their oncogenic signalling states in a biobank of patient-derived tumour xenograft (PDTX) models representing the diversity of human breast cancer. The BCMC panel identifies 13 cellular phenotypes (11 human and 2 murine), associated with both breast cancer subtypes and specific genomic features. Pre-treatment cellular phenotypic composition is a determinant of response to anticancer therapies. Single-cell profiling also reveals drug-induced cellular phenotypic dynamics, unravelling previously unnoticed intra-tumour response diversity. The comprehensive view of the landscapes of cellular phenotypic heterogeneity in PDTXs uncovered by the BCMC panel, which is mirrored in primary human tumours, has profound implications for understanding and predicting therapy response and resistance.

DNA methylation landscapes of 1538 breast cancers reveal a replication-linked clock, epigenomic instability and cis-regulation.

DNA methylation is aberrant in cancer, but the dynamics, regulatory role and clinical implications of such epigenetic changes are still poorly understood. Here, reduced representation bisulfite sequencing (RRBS) profiles of 1538 breast tumors and 244 normal breast tissues from the METABRIC cohort are reported, facilitating detailed analysis of DNA methylation within a rich context of genomic, transcriptional, and clinical data. Tumor methylation from immune and stromal signatures are deconvoluted leading to the discovery of a tumor replication-linked clock with genome-wide methylation loss in non-CpG island sites. Unexpectedly, methylation in most tumor CpG islands follows two replication-independent processes of gain (MG) or loss (ML) that we term epigenomic instability. Epigenomic instability is correlated with tumor grade and stage, TP53 mutations and poorer prognosis. After controlling for these global trans-acting trends, as well as for X-linked dosage compensation effects, cis-specific methylation and expression correlations are uncovered at hundreds of promoters and over a thousand distal elements. Some of these targeted known tumor suppressors and oncogenes. In conclusion, this study demonstrates that global epigenetic instability can erode cancer methylomes and expose them to localized methylation aberrations in-cis resulting in transcriptional changes seen in tumors.

Imaging mass cytometry and multiplatform genomics define the phenogenomic landscape of breast cancer.

Genomic alterations shape cell phenotypes and the structure of tumor ecosystems in poorly defined ways. To investigate these relationships, we used imaging mass cytometry to quantify the expression of 37 proteins with subcellular spatial resolution in 483 tumors from the METABRIC cohort. Single-cell analysis revealed cell phenotypes spanning epithelial, stromal and immune types. Distinct combinations of cell phenotypes and cell-cell interactions were associated with genomic subtypes of breast cancer. Epithelial luminal cell phenotypes separated into those predominantly impacted by mutations and those affected by copy number aberrations. Several features of tumor ecosystems, including cellular neighborhoods, were linked to prognosis, illustrating their clinical relevance. In summary, systematic analysis of single-cell phenotypic and spatial correlates of genomic alterations in cancer revealed how genomes shape both the composition and architecture of breast tumor ecosystems and will enable greater understanding of the phenotypic impact of genomic alterations.