Search

Abstract Despite our understanding of the genetic basis of intra-tumoral heterogeneity, the role of stromal heterogeneity arising from an altered tumor microenvironment in affecting tumorigenesis is poorly understood. In particular, extensive study on the peri-tumoral stroma in the morphologically normal tissues surrounding the tumor is lacking. Here, we examine the heterogeneity in tumors and peri-tumoral stroma from 8 ER+/PR+/HER2− invasive breast carcinomas, through multi-region transcriptomic profiling by microarray. We describe the regional heterogeneity observed at the intrinsic molecular subtype, pathway enrichment, and cell type composition levels within each tumor and its peri-tumoral region, up to 7 cm from the tumor margins. Moreover, we identify a pro-inflammatory adipose-enriched peri-tumoral subtype which was significantly associated with poorer overall survival in breast cancer patients, in contrast to an adaptive immune cell- and myofibroblast-enriched subtype. These data together suggest that peri-tumoral heterogeneity may be an important determinant of the evolution and treatment of breast cancers.

Abstract B cells are known to contribute to the anti-tumor immune response, especially in immunogenic tumors such as melanoma, yet humoral immunity has not been characterized in these cancers to detail. Here we show comprehensive phenotyping in samples of circulating and tumor-resident B cells as well as serum antibodies in melanoma patients. Memory B cells are enriched in tumors compared to blood in paired samples and feature distinct antibody repertoires, linked to specific isotypes. Tumor-associated B cells undergo clonal expansion, class switch recombination, somatic hypermutation and receptor revision. Compared with blood, tumor-associated B cells produce antibodies with proportionally higher levels of unproductive sequences and distinct complementarity determining region 3 properties. The observed features are signs of affinity maturation and polyreactivity and suggest an active and aberrant autoimmune-like reaction in the tumor microenvironment. Consistent with this, tumor-derived antibodies are polyreactive and characterized by autoantigen recognition. Serum antibodies show reactivity to antigens attributed to autoimmune diseases and cancer, and their levels are higher in patients with active disease compared to post-resection state. Our findings thus reveal B cell lineage dysregulation with distinct antibody repertoire and specificity, alongside clonally-expanded tumor-infiltrating B cells with autoimmune-like features, shaping the humoral immune response in melanoma.

Abstract Background Chimeric antigen receptor (CAR)-T cells can induce powerful immune responses in patients with hematological malignancies but have had limited success against solid tumors. This is in part due to the immunosuppressive tumor microenvironment (TME) which limits the activity of tumor-infiltrating lymphocytes (TILs) including CAR-T cells. We have developed a next-generation armored CAR (F i-CAR) targeting receptor tyrosine kinase-like orphan receptor 1 (ROR1), which is expressed at high levels in a range of aggressive tumors including poorly prognostic triple-negative breast cancer (TNBC). The F i-CAR-T is designed to release an anti-PD-1 checkpoint inhibitor upon CAR-T cell activation within the TME, facilitating activation of CAR-T cells and TILs while limiting toxicity. Methods To bolster potency, we developed a F i-CAR construct capable of IL-2-mediated, NFAT-induced secretion of anti-PD-1 single-chain variable fragments (scFv) within the tumor microenvironment, following ROR1-mediated activation. Cytotoxic responses against TNBC cell lines as well as levels and binding functionality of released payload were analyzed in vitro by ELISA and flow cytometry. In vivo assessment of potency of F i-CAR-T cells was performed in a TNBC NSG mouse model. Results F i-CAR-T cells released measurable levels of anti-PD-1 payload with 5 h of binding to ROR1 on tumor and enhanced the cytotoxic effects at challenging 1:10 E:T ratios. Treatment of established PDL1 + TNBC xenograft model with F i-CAR-T cells resulted in significant abrogation in tumor growth and improved survival of mice (71 days), compared to non-armored CAR cells targeting ROR1 (F CAR-T) alone (49 days) or in combination with systemically administered anti-PD-1 antibody (57 days). Crucially, a threefold increase in tumor-infiltrating T cells was observed with F i-CAR-T cells and was associated with increased expression of genes related to cytotoxicity, migration and proliferation. Conclusions Our next-generation of ROR1-targeting inducible armored CAR platform enables the release of an immune stimulating payload only in the presence of target tumor cells, enhancing the therapeutic activity of the CAR-T cells. This technology provided a significant survival advantage in TNBC xenograft models. This coupled with its potential safety attributes merits further clinical evaluation of this approach in TNBC patients.

Abstract The level of stromal tumor-infiltrating lymphocytes (sTILs) in triple-negative (TNBC) and HER2-positive breast cancers convey prognostic information. The importance of systemic immunity to local immunity is unknown in breast cancer. We previously demonstrated that histological alterations in axillary lymph nodes (LNs) carry clinical relevance. Here, we capture local immune responses by scoring TILs at the primary tumor and systemic immune responses by recording the formation of secondary follicles, also known as germinal centers, in 2,857 cancer-free and involved axillary LNs on haematoxylin and eosin (H&E) stained sections from a retrospective cohort of 161 LN-positive triple-negative and HER2-positive breast cancer patients. Our data demonstrate that the number of germinal center formations across all cancer-free LNs, similar to high levels of TILs, is associated with a good prognosis in low TILs TNBC. This highlights the importance of assessing both primary and LN immune responses for prognostication and for future breast cancer research.

Background: Major advances in breast cancer treatment have led to a reducuction in mortality. However, there are still women who are not cured. We hypothesize there is a sub-group of women with treatment-resistant cancers causing early death. Methods: Between 1975 and 2006, 5392 women with invasive breast cancer underwent surgery at Guy’s Hospital, London. Data on patient demographics, tumour characteristics, treatment regimens, local recurrence, secondary metastasis, and death were prospectively recorded. We considered four time periods (1975–1982, 1983–1990, 1991–1998, 1999–2006). Risks and time to event analysis were performed with Cox proportional hazards model and Kaplan-Meier estimation. Results: Unadjusted hazard ratios for developing metastasis and overall mortality relative to the 1975–1982 cohort decreased steadily to 0.23 and 0.63, respectively in 1999–2006. However, metastasis-free interval shortened, with the proportion of women developing metastasis ≤5 years increasing from 73.9% to 83.0%. Furthermore, median post-metastatic survival decreased from 1.49 years to 0.94 years. Applying our risk criteria identified the presence of ±200 patients in each cohort who developed metastasis early and died within a much shorter time frame. Conclusions: Advances in treatment have decreased the risk of metastasis and improved survival in women with invasive breast cancer over the last 40 years. Despite this, a subpopulation with shorter metastasis-free and post-metastatic survival who are unresponsive to available treatment remains. This may be due to the ATRESS phenomenon (adjuvant therapy-related shortening of survival) secondary to preselection inherent in adjuvant therapy, successful treatment of less malignant tumour cells and treatment-induced resistance in the remaining tumour clones.

Abstract Background Metabolites are genetically and environmentally determined. Consequently, they can be used to characterize environmental exposures and reveal biochemical mechanisms that link exposure to disease. To explore disease susceptibility and improve population risk stratification, we aimed to identify metabolic profiles linked to carcinogenesis and mortality and their intrinsic associations by characterizing subgroups of individuals based on serum biomarker measurements. We included 13,615 participants from the Swedish Apolipoprotein MOrtality RISk Study who had measurements for 19 biomarkers representative of central metabolic pathways. Latent Class Analysis (LCA) was applied to characterise individuals based on their biomarker values (according to medical cut-offs), which were then examined as predictors of cancer and death using multivariable Cox proportional hazards models. Results LCA identified four metabolic profiles within the population: (1) normal values for all markers (63% of population); (2) abnormal values for lipids (22%); (3) abnormal values for liver functioning (9%); (4) abnormal values for iron and inflammation metabolism (6%). All metabolic profiles (classes 2–4) increased risk of cancer and mortality, compared to class 1 (e.g. HR for overall death was 1.26 (95% CI: 1.16–1.37), 1.67 (95% CI: 1.47–1.90), and 1.21 (95% CI: 1.05–1.41) for class 2, 3, and 4, respectively). Conclusion We present an innovative approach to risk stratify a well-defined population based on LCA metabolic-defined subgroups for cancer and mortality. Our results indicate that standard of care baseline serum markers, when assembled into meaningful metabolic profiles, could help assess long term risk of disease and provide insight in disease susceptibility and etiology.

Lymph nodes (LNs) are highly organized secondary lymphoid organs, and reflective of immune responses to infection, injuries, or the presence of cancer. Extensive molecular and morphological analyses of immune and stromal features in tumors and LNs of breast cancer patients have revealed novel patterns indicative of disease progression. Within LNs, there are dynamic structures called germinal centers (GCs), that act as the immunological hubs for B cell development and generation of affinity matured memory B and antibody-producing plasma cells. Acting as a bridge between systemic and local immunity, associations are observed between the frequency of GCs within cancer-free LNs, the levels of stromal tumor infiltrating lymphocytes, and cancer progression. Scattered throughout the tumor microenvironment (TME) or aggregated in clusters forming tertiary lymphoid structures (TLS), the occurrence of tumor infiltrating B cells (TIL-Bs) has been linked mostly to superior disease trajectories in solid cancers. Recent TIL-Bs profiling studies have revealed a plethora of different TIL-B populations, their functional roles, and whether they are derived from GC reactions in the LN, and/or locally from GC-like structures within the TME remains to be investigated. However, parallels between the immunogenic nature of LNs as a pre-metastatic niche, TIL-B populations within the TME, and the presence of TLS will help to decipher local and widespread TIL-Bs responses and their influence on cancer progression to the lymphatics. Therapies that enhance TIL-Bs responses in the LN GC and/or in GC-like structures in the TME are thus emerging management strategies for breast and other cancer patients.

Abstract Lymph node (LN) metastasis is an important prognostic parameter in breast carcinoma, a crucial site for tumour–immune cell interaction and a gateway for further dissemination of tumour cells to other metastatic sites. To gain insight into the underlying molecular changes from the pre-metastatic, via initial colonisation to the fully involved LN, we reviewed transcriptional research along the evolving microenvironment of LNs in human breast cancers patients. Gene expression studies were compiled and subjected to pathway-based analyses, with an emphasis on immune cell-related genes. Of 366 studies, 14 performed genome-wide gene expression comparisons and were divided into six clinical-biological scenarios capturing different stages of the metastatic pathway in the LN, as follows: metastatically involved LNs are compared to their patient-matched primary breast carcinomas (scenario 1) or the normal breast tissue (scenario 2). In scenario 3, uninvolved LNs were compared between LN-positive patients and LN-negative patients. Scenario 4 homed in on the residual uninvolved portion of involved LNs and compared it to the patient-matched uninvolved LNs. Scenario 5 contrasted uninvolved and involved LNs, whilst in scenario 6 involved (sentinel) LNs were assessed between patients with other either positive or negative LNs (non-sentinel). Gene lists from these chronological steps of LN metastasis indicated that gene patterns reflecting deficiencies in dendritic cells and hyper-proliferation of B cells parallel to tumour promoting pathways, including cell adhesion, extracellular matrix remodelling, cell motility and DNA repair, play key roles in the changing microenvironment of a pro-metastatic to a metastatically involved LN. Similarities between uninvolved LNs and the residual uninvolved portion of involved LNs hinted that LN alterations expose systemic tumour-related immune responses in breast cancer patients. Despite the diverse settings, gene expression patterns at different stages of metastatic colonisation in LNs were recognised and may provide potential avenues for clinical interventions to counteract disease progression for breast cancer patients.

Triple negative breast cancers harbor multiple copy number aberrations driving gene expression changes thought to underpin their malignant phenotypes. Here the authors integrate these features, finding and functionally validating 37 gene addictions among which they identify the mechanism of addiction to KIFC1, a potential selective drug target.

Abstract Background Cell migration is essential for development and tissue repair, but it also contributes to disease. Rho GTPases regulate cell migration, but a comprehensive analysis of how each Rho signalling component affects migration has not been carried out. Results Through an RNA interference screen, and using a prostate cancer cell line, we find that approximately 25% of Rho network components alter migration. Some genes enhance migration while others decrease basal and/or hepatocyte growth factor-stimulated migration. Surprisingly, we identify RhoH as a screen hit. RhoH expression is normally restricted to haematopoietic cells, but we find it is expressed in multiple epithelial cancer cell lines. High RhoH expression in samples from prostate cancer patients correlates with earlier relapse. RhoH depletion reduces cell speed and persistence and decreases migratory polarity. Rac1 activity normally localizes to the front of migrating cells at areas of dynamic membrane movement, but in RhoH-depleted cells active Rac1 is localised around the whole cell periphery and associated with membrane regions that are not extending or retracting. RhoH interacts with Rac1 and with several p21-activated kinases (PAKs), which are Rac effectors. Similar to RhoH depletion, PAK2 depletion increases cell spread area and reduces cell migration. In addition, RhoH depletion reduces lamellipodium extension induced by PAK2 overexpression. Conclusions We describe a novel role for RhoH in prostate cancer cell migration. We propose that RhoH promotes cell migration by coupling Rac1 activity and PAK2 to membrane protrusion. Our results also suggest that RhoH expression levels correlate with prostate cancer progression.

Abstract Background Metastases from primary breast cancers can involve single or multiple organs at metastatic disease diagnosis. Molecular risk factors for particular patterns of metastastic spread in a clinical population are limited. Methods A case-control design including 1357 primary breast cancers was used to study three distinct clinical patterns of metastasis, which occur within the first six months of metastatic disease: bone and visceral metasynchronous spread, bone-only, and visceral-only metastasis. Whole-genome expression profiles were obtained using whole genome (WG)-DASL assays from formalin-fixed paraffin-embedded (FFPE) samples. A systematic protocol was developed for handling FFPE samples together with stringent data quality controls to identify robust expression profiling data. A panel of published and novel gene sets were tested for association with these specific patterns of metastatic spread and odds ratios (ORs) were calculated. Results Metasynchronous metastasis to bone and viscera was found in all intrinsic breast cancer subtypes, while immunohistochemically (IHC)-defined receptor status and specific IntClust subgroups were risk factors for visceral-only or bone-only first metastases. Among gene modules, those related to proliferation increased the risk of metasynchronous metastasis (OR (95% CI) = 2.3 (1.1–4.8)) and visceral-only first metastasis (OR (95% CI) = 2.5 (1.2–5.1)) but not bone-only metastasis (OR (95% CI) = 0.97 (0.56–1.7)). A 21-gene module (BV) was identified in estrogen-receptor-positive breast cancers with metasynchronous metastasis to bone and viscera (area under the curve = 0.77), and its expression increased the risk of bone and visceral metasynchronous spread in this population. BV was further orthogonally validated with NanoString nCounter in primary breast cancers, and was reproducible in their matched lymph nodes metastases and an external cohort. Conclusion This case-control study of WG-DASL global expression profiles from FFPE tumour samples, after careful quality control and RNA selection, revealed that gene modules in the primary tumour have differing risks for clinical patterns of metasynchronous first metastases. Moreover, a novel gene module was identified as a putative risk factor for metasynchronous bone and visceral first metastatic spread, with potential implications for disease monitoring and treatment planning.

Summary: Lymphatic vasculature is crucial for metastasis in triple-negative breast cancer (TNBC); however, cellular and molecular drivers controlling lymphovascular metastasis are poorly understood. We define a macrophage-dependent signaling cascade that facilitates metastasis through lymphovascular remodeling. TNBC cells instigate mRNA changes in macrophages, resulting in β4 integrin-dependent adhesion to the lymphovasculature. β4 integrin retains macrophages proximal to lymphatic endothelial cells (LECs), where release of TGF-β1 drives LEC contraction via RhoA activation. Macrophages promote gross architectural changes to lymphovasculature by increasing dilation, hyperpermeability, and disorganization. TGF-β1 drives β4 integrin clustering at the macrophage plasma membrane, further promoting macrophage adhesion and demonstrating the dual functionality of TGF-β1 signaling in this context. β4 integrin-expressing macrophages were identified in human breast tumors, and a combination of vascular-remodeling macrophage gene signature and TGF-β signaling scores correlates with metastasis. We postulate that future clinical strategies for patients with TNBC should target crosstalk between β4 integrin and TGF-β1. : Breast cancer metastasis through lymphatic vessels is associated with poor prognosis. Evans et al. describe β4 integrin-expressing macrophages that regulate lymphatic vessel structure in breast cancer. Macrophage-released TGF-β1 drives lymphatic cell contraction via RhoA activation, culminating in lymphatic hyperpermeability. This study defines a signaling cascade that could be targeted therapeutically. Keywords: lymphovasculature, macrophages, cancer, remodeling, adhesion, contraction, β4 integrin, TGF-β1, RhoA

Summary: The SRC-family kinase LYN is highly expressed in triple-negative/basal-like breast cancer (TNBC) and in the cell of origin of these tumors, c-KIT-positive luminal progenitors. Here, we demonstrate LYN is a downstream effector of c-KIT in normal mammary cells and protective of apoptosis upon genotoxic stress. LYN activity is modulated by PIN1, a prolyl isomerase, and in BRCA1 mutant TNBC PIN1 upregulation activates LYN independently of c-KIT. Furthermore, the full-length LYN splice isoform (as opposed to the Δaa25–45 variant) drives migration and invasion of aggressive TNBC cells, while the ratio of splice variants is informative for breast cancer-specific survival across all breast cancers. Thus, dual mechanisms—uncoupling from upstream signals and splice isoform ratios—drive the activity of LYN in aggressive breast cancers. : Tornillo et al. show that in aggressive breast cancers, LYN activity is deregulated by a change in patterns of splice isoform expression. In BRCA1-dysfunctional breast cancers, LYN activity is upregulated by a prolyl isomerase (PIN1) that is normally repressed by BRCA1. Keywords: BRCA1, LYN kinase, triple-negative/basal-like breast cancer, PIN1, ESRP1, c-KIT

Breast cancers with BRCA1 germline mutation have a characteristic DNA copy number (CN) pattern. We developed a test that assigns CN profiles to be ‘BRCA1‐like’ or ‘non‐BRCA1‐like’, which refers to resembling a BRCA1‐mutated tumor or resembling a tumor without a BRCA1 mutation, respectively. Approximately one third of the BRCA1‐like breast cancers have a BRCA1 mutation, one third has hypermethylation of the BRCA1 promoter and one third has an unknown reason for being BRCA1‐like. This classification is indicative of patients' response to high dose alkylating and platinum containing chemotherapy regimens, which targets the inability of BRCA1 deficient cells to repair DNA double strand breaks. We investigated whether this classification can be reliably obtained with next generation sequencing and copy number platforms other than the bacterial artificial chromosome (BAC) array Comparative Genomic Hybridization (aCGH) on which it was originally developed. We investigated samples from 230 breast cancer patients for which a CN profile had been generated on two to five platforms, comprising low coverage CN sequencing, CN extraction from targeted sequencing panels (CopywriteR), Affymetrix SNP6.0, 135K/720K oligonucleotide aCGH, Affymetrix Oncoscan FFPE (MIP) technology, 3K BAC and 32K BAC aCGH. Pairwise comparison of genomic position‐mapped profiles from the original aCGH platform and other platforms revealed concordance. For most cases, biological differences between samples exceeded the differences between platforms within one sample. We observed the same classification across different platforms in over 80% of the patients and kappa values of at least 0.36. Differential classification could be attributed to CN profiles that were not strongly associated to one class. In conclusion, we have shown that the genomic regions that define our BRCA1‐like classifier are robustly measured by different CN profiling technologies, providing the possibility to retro‐ and prospectively investigate BRCA1‐like classification across a wide range of CN platforms.

Molecular, cellular, and clinical studies have combined to demonstrate a contribution from the DNA cytosine deaminase APOBEC3B (A3B) to the overall mutation load in breast, head/neck, lung, bladder, cervical, ovarian, and other cancer types. However, the complete landscape of mutations attributable to this enzyme has yet to be determined in a controlled human cell system. We report a conditional and isogenic system for A3B induction, genomic DNA deamination, and mutagenesis. Human 293-derived cells were engineered to express doxycycline-inducible A3B-eGFP or eGFP constructs. Cells were subjected to 10 rounds of A3B-eGFP exposure that each caused 80-90% cell death. Control pools were subjected to parallel rounds of non-toxic eGFP exposure, and dilutions were done each round to mimic A3B-eGFP induced population fluctuations. Targeted sequencing of portions of TP53 and MYC demonstrated greater mutation accumulation in the A3B-eGFP exposed pools. Clones were generated and microarray analyses were used to identify those with the greatest number of SNP alterations for whole genome sequencing. A3B-eGFP exposed clones showed global increases in C-to-T transition mutations, enrichments for cytosine mutations within A3B-preferred trinucleotide motifs, and more copy number aberrations. Surprisingly, both control and A3B-eGFP clones also elicited strong mutator phenotypes characteristic of defective mismatch repair. Despite this additional mutational process, the 293-based system characterized here still yielded a genome-wide view of A3B-catalyzed mutagenesis in human cells and a system for additional studies on the compounded effects of simultaneous mutation mechanisms in cancer cells.

Despite the wide use of cell lines in cancer research, the extent to which their surface properties correspond to those of primary tumors is poorly characterized. The present study addresses this problem from a transcriptional standpoint, analyzing the expression of membrane protein genes--the Membranome--in primary tumors and immortalized in-vitro cultured tumor cells. 409 human samples, deriving from ten independent studies, were analyzed. These comprise normal tissues, primary tumors and tumor derived cell lines deriving from eight different tissues: brain, breast, colon, kidney, leukemia, lung, melanoma, and ovary. We demonstrated that the Membranome has greater power than the remainder of the transcriptome when used as input for the automatic classification of tumor samples. This feature is maintained in tumor derived cell lines. In most cases primary tumors show maximal similarity in Membranome expression with cell lines of same tissue origin. Differences in Membranome expression between tumors and cell lines were analyzed also at the pathway level and biological themes were identified that were differentially regulated in the two settings. Moreover, by including normal samples in the analysis, we quantified the degree to which cell lines retain the Membranome up- and down-regulations observed in primary tumors with respect to their normal counterparts. We showed that most of the Membranome up-regulations observed in primary tumors are lost in the in-vitro cultured cells. Conversely, the majority of Membranome genes down-regulated upon tumor transformation maintain lower expression levels also in the cell lines. This study points towards a central role of Membranome genes in the definition of the tumor phenotype. The comparative analysis of primary tumors and cell lines identifies the limits of cell lines as a model for the study of cancer-related processes mediated by the cell surface. Results presented allow for a more rational use of the cell lines as a model of cancer.

The majority of new drug approvals for cancer are based on existing therapeutic targets. One approach to the identification of novel targets is to perform high-throughput RNA interference (RNAi) cellular viability screens. We describe a novel approach combining RNAi screening in multiple cell lines with gene expression and genomic profiling to identify novel cancer targets. We performed parallel RNAi screens in multiple cancer cell lines to identify genes that are essential for viability in some cell lines but not others, suggesting that these genes constitute key drivers of cellular survival in specific cancer cells. This approach was verified by the identification of PIK3CA, silencing of which was selectively lethal to the MCF7 cell line, which harbours an activating oncogenic PIK3CA mutation. We combined our functional RNAi approach with gene expression and genomic analysis, allowing the identification of several novel kinases, including WEE1, that are essential for viability only in cell lines that have an elevated level of expression of this kinase. Furthermore, we identified a subset of breast tumours that highly express WEE1 suggesting that WEE1 could be a novel therapeutic target in breast cancer. In conclusion, this strategy represents a novel and effective strategy for the identification of functionally important therapeutic targets in cancer.

Although paediatric high grade gliomas resemble their adult counterparts in many ways, there appear to be distinct clinical and biological differences. One important factor hampering the development of new targeted therapies is the relative lack of cell lines derived from childhood glioma patients, as it is unclear whether the well-established adult lines commonly used are representative of the underlying molecular genetics of childhood tumours. We have carried out a detailed molecular and phenotypic characterisation of a series of paediatric high grade glioma cell lines in comparison to routinely used adult lines.All lines proliferate as adherent monolayers and express glial markers. Copy number profiling revealed complex genomes including amplification and deletions of genes known to be pivotal in core glioblastoma signalling pathways. Expression profiling identified 93 differentially expressed genes which were able to distinguish between the adult and paediatric high grade cell lines, including a number of kinases and co-ordinated sets of genes associated with DNA integrity and the immune response.These data demonstrate that glioma cell lines derived from paediatric patients show key molecular differences to those from adults, some of which are well known, whilst others may provide novel targets for evaluation in primary tumours. We thus provide the rationale and demonstrate the practicability of using paediatric glioma cell lines for preclinical and mechanistic studies.