Search

Breast cancers in humans belong to one of several intrinsic molecular subtypes each with different tumor biology and different clinical impact. Mammary gland tumors in dogs are proposed as a relevant comparative model for human breast cancer; however, it is still unclear whether the intrinsic molecular subtypes have the same significance in dogs and humans. Using publicly available data, we analyzed gene expression and whole-exome sequencing data from 158 canine mammary gland tumors. We performed molecular subtyping using the PAM50 method followed by subtype-specific comparisons of gene expression characteristics, mutation patterns and copy number profiles between canine tumors and human breast tumors from The Cancer Genome Atlas (TCGA) breast cancer cohort (n = 1097). We found that luminal A canine tumors greatly resemble luminal A human tumors both in gene expression characteristics, mutations and copy number profiles. Also, the basal-like canine and human tumors were relatively similar, with low expression of luminal epithelial markers and high expression of genes involved in cell proliferation. There were, however, distinct differences in immune-related gene expression patterns in basal-like tumors between the two species. Characteristic HER2-enriched and luminal B subtypes were not present in the canine cohort, and we found no tumors with high-level ERBB2 amplifications. Benign and malignant canine tumors displayed similar PAM50 subtype characteristics. Our findings indicate that deeper understanding of the different molecular subtypes in canine mammary gland tumors will further improve the value of canines as comparative models for human breast cancer.

The claudin-low breast cancer subtype is defined by gene expression characteristics and encompasses a remarkably diverse range of breast tumors. Here, we investigate genomic, transcriptomic, and clinical features of claudin-low breast tumors. We show that claudin-low is not simply a subtype analogous to the intrinsic subtypes (basal-like, HER2-enriched, luminal A, luminal B and normal-like) as previously portrayed, but is a complex additional phenotype which may permeate breast tumors of various intrinsic subtypes. Claudin-low tumors are distinguished by low genomic instability, mutational burden and proliferation levels, and high levels of immune and stromal cell infiltration. In other aspects, claudin-low tumors reflect characteristics of their intrinsic subtype. Finally, we explore an alternative method for identifying claudin-low tumors and thereby uncover potential weaknesses in the established claudin-low classifier. In sum, these findings elucidate the heterogeneity in claudin-low breast tumors, and substantiate a re-definition of claudin-low as a cancer phenotype., In breast cancer, the claudin-low breast cancer subtype is remarkably diverse. Here, the authors propose that claudin-low is not a classical intrinsic breast cancer subtype, but rather a complex additional phenotype that can occur across intrinsic subtypes.

How mixtures of immune cells associate with cancer cell phenotype and affect pathogenesis is still unclear. In 15 breast cancer gene expression datasets, we invariably identify three clusters of patients with gradual levels of immune infiltration. The intermediate immune infiltration cluster (Cluster B) is associated with a worse prognosis independently of known clinicopathological features. Furthermore, immune clusters are associated with response to neoadjuvant chemotherapy. In silico dissection of the immune contexture of the clusters identified Cluster A as immune cold, Cluster C as immune hot while Cluster B has a pro-tumorigenic immune infiltration. Through phenotypical analysis, we find epithelial mesenchymal transition and proliferation associated with the immune clusters and mutually exclusive in breast cancers. Here, we describe immune clusters which improve the prognostic accuracy of immune contexture in breast cancer. Our discovery of a novel independent prognostic factor in breast cancer highlights a correlation between tumor phenotype and immune contexture., In breast cancer, the immune infiltration of the tumour associates with clinical outcome. Here, the authors infer immune context based on gene expression data and identify a new independent subtype linked to pro-tumorigenic immune infiltration.

The large heterogeneity of DCIS complicates the identification of markers for risk of progression to invasive breast cancer. Molecular and cellular heterogeneity is observed across DCIS, between multiple foci of DCIS within a lesion and within mixed lesions of invasive tumors with DCIS foci. The intrinsic breast cancer subtypes represent fundamental biological properties of the disease and are also found at the DCIS stage. However, their frequency and relevance as prognostic markers in pre-invasive stage is still not clear. Novel technologies for spatial characterization of tumors allow detailed mapping of DCIS and the microenvironment and may give us clues to which DCIS have the propensity to invade surrounding stroma. In the presentation, we will discuss results from our ongoing genomic characterization of DCIS and invasive breast carcinomas that illustrate the need for stratification when searching for prognostic markers. We have found disproportionate distribution and characteristics of the intrinsic subtypes in DCIS compared with invasive breast cancer and the dominance of Luminal A and HER2 subtypes. Bona-fide basal-like tumors may rarely be found at the DCIS stage and HER2-enriched DCIS are associated with specific immune cell profiles. Citation Format: Therese Sørlie. Molecular subtypes and spatial heterogeneity in DCIS [abstract]. In: Proceedings of the AACR Special Conference on Rethinking DCIS: An Opportunity for Prevention?; 2022 Sep 8-11; Philadelphia, PA. Philadelphia (PA): AACR; Can Prev Res 2022;15(12 Suppl_1): Abstract nr IA030.

Simple Summary In breast cancer, there is a high degree of variability in tumors and the surrounding tissue called the tumor microenvironment (TME). To better understand tumor biology and metastasis, as well as to predict response to cancer treatments or the course of the disease, it is important to characterize molecular diversity in the breast TME. The GeoMx Digital Spatial Profiler (DSP) enables researchers to spatially analyze proteins and RNA transcripts in tumors and surrounding tissues from patients or preclinical models. Using the GeoMx DSP, protein expression and RNA transcripts in the distinct regions of a tumor can be quantified up to and including the whole transcriptome level. Herein, the GeoMx Breast Cancer Consortium presents best practices for GeoMx spatial profiling of tumors to promote the collection of high-quality data, optimization of data analysis and integration of datasets to accelerate biomarker discovery. These best practices can also be applied to any tumor type to provide information about the tumor and the TME. Abstract Breast cancer is a heterogenous disease with variability in tumor cells and in the surrounding tumor microenvironment (TME). Understanding the molecular diversity in breast cancer is critical for improving prediction of therapeutic response and prognostication. High-plex spatial profiling of tumors enables characterization of heterogeneity in the breast TME, which can holistically illuminate the biology of tumor growth, dissemination and, ultimately, response to therapy. The GeoMx Digital Spatial Profiler (DSP) enables researchers to spatially resolve and quantify proteins and RNA transcripts from tissue sections. The platform is compatible with both formalin-fixed paraffin-embedded and frozen tissues. RNA profiling was developed at the whole transcriptome level for human and mouse samples and protein profiling of 100-plex for human samples. Tissue can be optically segmented for analysis of regions of interest or cell populations to study biology-directed tissue characterization. The GeoMx Breast Cancer Consortium (GBCC) is composed of breast cancer researchers who are developing innovative approaches for spatial profiling to accelerate biomarker discovery. Here, the GBCC presents best practices for GeoMx profiling to promote the collection of high-quality data, optimization of data analysis and integration of datasets to advance collaboration and meta-analyses. Although the capabilities of the platform are presented in the context of breast cancer research, they can be generalized to a variety of other tumor types that are characterized by high heterogeneity.

Claudin-low breast cancer is a molecular subtype associated with poor prognosis and without targeted treatment options. The claudin-low subtype is defined by certain biological characteristics, some of which may be clinically actionable, such as high immunogenicity. In mice, the medroxyprogesterone acetate (MPA) and 7,12-dimethylbenzanthracene (DMBA) induced mammary tumor model yields a heterogeneous set of tumors, a subset of which display claudin-low features. Neither the genomic characteristics of MPA/DMBA-induced claudin-low tumors, nor those of human claudin-low breast tumors, have been thoroughly explored.The transcriptomic characteristics and subtypes of MPA/DMBA-induced mouse mammary tumors were determined using gene expression microarrays. Somatic mutations and copy number aberrations in MPA/DMBA-induced tumors were identified from whole exome sequencing data. A publicly available dataset was queried to explore the genomic characteristics of human claudin-low breast cancer and to validate findings in the murine tumors.Half of MPA/DMBA-induced tumors showed a claudin-low-like subtype. All tumors carried mutations in known driver genes. While the specific genes carrying mutations varied between tumors, there was a consistent mutational signature with an overweight of T>A transversions in TG dinucleotides. Most tumors carried copy number aberrations with a potential oncogenic driver effect. Overall, several genomic events were observed recurrently, however none accurately delineated claudin-low-like tumors. Human claudin-low breast cancers carried a distinct set of genomic characteristics, in particular a relatively low burden of mutations and copy number aberrations. The gene expression characteristics of claudin-low-like MPA/DMBA-induced tumors accurately reflected those of human claudin-low tumors, including epithelial-mesenchymal transition phenotype, high level of immune activation and low degree of differentiation. There was an elevated expression of the immunosuppressive genes PTGS2 (encoding COX-2) and CD274 (encoding PD-L1) in human and murine claudin-low tumors. Our findings show that the claudin-low breast cancer subtype is not demarcated by specific genomic aberrations, but carries potentially targetable characteristics warranting further research.Author SummaryBreast cancer is comprised of several distinct disease subtypes with different etiologies, prognoses and therapeutic targets. The claudin-low breast cancer subtype is relatively poorly understood, and no specific treatment exists targeting its unique characteristics. Animal models accurately representing human disease counterparts are vital for developing novel therapeutics, but for the claudin-low breast cancer subtype, no such uniform model exists. Here, we show that exposing mice to the carcinogen DMBA and the hormone MPA causes a diverse range of mammary tumors to grow, and half of these have a gene expression pattern similar to that seen in human claudin-low breast cancer. These tumors have numerous changes in their DNA, with clear differences between each tumor, however no specific DNA aberrations clearly demarcate the claudin-low subtype. We also analyzed human breast cancers and show that human claudin-low tumors have several clear patterns in their DNA aberrations, but no specific features accurately distinguish claudin-low from non-claudin-low breast cancer. Finally, we show that both human and murine claudin-low tumors express high levels of genes associated with suppression of immune response. In sum, we highlight claudin-low breast cancer as a clinically relevant subtype with a complex etiology, and with potential unexploited therapeutic targets.

Physical inactivity, high postmenopausal body mass index, alcohol consumption and use of menopausal hormone therapy are established risk factors for breast cancer. Less is known about whether these factors influence the risk of progression of benign and premalignant breast lesions to invasive breast cancer. This registry-based cohort study was based on women with a precancerous lesion who were followed for breast cancer. The cohort consisted of 11 270 women with a benign lesion, 972 women with hyperplasia with atypia and 2379 women with carcinoma in situ diagnosed and treated after participation in BreastScreen Norway, 2006-2016. Information on breast cancer risk factors was collected by a questionnaire administered with the invitation letter. Cox regression analysis was used to estimate the association between breast cancer and physical activity, body mass index, alcohol consumption, tobacco smoking and menopausal hormone therapy, adjusted for age. During follow-up, 274 women with a benign lesion, 34 women with hyperplasia with atypia and 118 women with carcinoma in situ were diagnosed with invasive breast cancer. We observed an increased risk of breast cancer associated with use of menopausal hormone therapy for women with a benign or premalignant lesion. Alcohol consumption and tobacco smoking showed suggestive increased risk of breast cancer among women with a benign lesion. We were only to a limited degree able to identify associations between modifiable risk factors of breast cancer and the disease among women with a precancerous lesion, and a larger study is needed to confirm or refute associations. The work was funded by The Norwegian Cancer Society (grant no. 5746604).

Background Ductal carcinoma in situ (DCIS) comprises a diverse group of preinvasive lesions in the breast and poses a considerable clinical challenge due to lack of markers of progression. Genomic alterations are to a large extent similar in DCIS and invasive carcinomas, although differences in copy number aberrations, gene expression patterns, and mutations exist. In mixed tumors with synchronous invasive breast cancer (IBC) and DCIS, it is still unclear to what extent invasive tumor cells are directly derived from the DCIS cells. Aim Our aim was to compare cancer‐relevant mutation profiles of different cellular compartments in mixed DCIS/IBC and pure DCIS tumors. Methods and results We performed targeted sequencing of 50 oncogenes in microdissected tissue from three different epithelial cell compartments (in situ, invasive, and normal adjacent epithelium) from 26 mixed breast carcinomas. In total, 44 tissue samples (19 invasive, 16 in situ, 9 normal) were subjected to sequencing using the Ion Torrent platform and the AmpliSeq Cancer Hotspot Panel v2. For comparison, 10 additional, pure DCIS lesions were sequenced. Across all mixed samples, we detected 23 variants previously described in cancer. The most commonly affected genes were TP53, PIK3CA, and ERBB2. The PIK3CA:p.H1047R variant was found in nine samples from six patients. Most variants detected in invasive compartments were also found in the corresponding in situ cell compartment indicating a clonal relationship between the tumor stages. A lower frequency of variants were observed in pure DCIS lesions. Conclusion Similar mutation profiles between in situ and invasive cell compartments indicate a similar origin of the two tumor stages in mixed breast tumors. The lower number of potential driver variants found in pure DCIS compared with the in situ cell compartments of mixed tumors may imply that pure DCIS is captured earlier in the path of progression to invasive disease.

Summary The receptor tyrosine kinase AXL is associated with epithelial plasticity in several solid tumors including breast cancer and AXL-targeting agents are currently in clinical trials. We hypothesized that AXL is a driver of stemness traits in cancer by co-option of a regulatory function normally reserved for stem cells. AXL-expressing cells in human mammary epithelial ducts co-expressed markers associated with multipotency, and AXL inhibition abolished colony formation and self-maintenance activities while promoting terminal differentiation in vitro. Axl-null mice did not exhibit a strong developmental phenotype, but enrichment of Axl+ cells was required for mouse mammary gland reconstitution upon transplantation, and Axl-null mice had reduced incidence of Wnt1-driven mammary tumors. An AXL-dependent gene signature is a feature of transcriptomes in basal breast cancers and reduced patient survival irrespective of subtype. Our interpretation is that AXL regulates access to epithelial plasticity programs in MaSCs and, when co-opted, maintains acquired stemness in breast cancer cells., Graphical Abstract, Highlights • AXL + mammary epithelial cells have multipotent activity conserved in women and mice • AXL allows accesses to epithelial-to-mesenchymal transition genes and prevents differentiation into luminal cells • Deletion of Axl reduced incidence of Wnt1-driven tumors in mice • Provides a rationale explaining the advantage to cancer cells that co-opt AXL signaling, Cell Biology; Stem Cells Research; Cancer

The long non-coding RNA NEAT1 locus is transcribed into two overlapping isoforms, NEAT1_1 and NEAT1_2, of which the latter is essential for the assembly of nuclear paraspeckles. NEAT1 is abnormally expressed in a wide variety of human cancers. Emerging evidence suggests that the two isoforms have distinct functions in gene expression regulation, and recently it was shown that NEAT1_2, but not NEAT1_1, expression predicts poor clinical outcome in cancer. Here, we report that NEAT1_2 expression correlates with HER2-positive breast cancers and high-grade disease. We provide evidence that NEAT1_1 and NEAT1_2 have distinct expression pattern among different intrinsic breast cancer subtypes. Finally, we show that NEAT1_2 expression and paraspeckle formation increase upon lactation in humans, confirming what has previously been demonstrated in mice. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Ductal carcinoma in situ (DCIS) is a non-invasive type of breast cancer with highly variable potential of becoming invasive and affecting mortality. Currently, many patients with DCIS are overtreated due to the lack of specific biomarkers that distinguish low risk lesions from those with a higher risk of progression. In this study, we analyzed 57 pure DCIS and 313 invasive breast cancers (IBC) from different patients. Three levels of genomic data were obtained; gene expression, DNA methylation, and DNA copy number. We performed subtype stratified analyses and identified key differences between DCIS and IBC that suggest subtype specific progression. Prominent differences were found in tumors of the basal-like subtype: Basal-like DCIS were less proliferative and showed a higher degree of differentiation than basal-like IBC. Also, core basal tumors (characterized by high correlation to the basal-like centroid) were not identified amongst DCIS as opposed to IBC. At the copy number level, basal-like DCIS exhibited fewer copy number aberrations compared with basal-like IBC. An intriguing finding through analysis of the methylome was hypermethylation of multiple protocadherin genes in basal-like IBC compared with basal-like DCIS and normal tissue, possibly caused by long range epigenetic silencing. This points to silencing of cell adhesion-related genes specifically in IBC of the basal-like subtype. Our work confirms that subtype stratification is essential when studying progression from DCIS to IBC, and we provide evidence that basal-like DCIS show less aggressive characteristics and question the assumption that basal-like DCIS is a direct precursor of basal-like invasive breast cancer.

Simple Summary The PAM50 risk of recurrence (ROR) score is predictive of the risk of distant recurrence and the benefits of adjuvant therapy in early breast cancer. The Prosigna assay utilizes RNA from tumor cells selected via macrodissection of formalin-fixed, paraffin-embedded (FFPE) tissue sections to measure the activity of the PAM50 genes. An alternative and widely used extraction method is RNA purification from fresh-frozen (FF) bulk tissue without enriching the tumor cellularity. However, the impact of the RNA preparation approach on ROR scores and subsequent treatment selection has not been systematically evaluated. We compared the different approaches and found high correlation between risk of recurrence scores estimated from macrodissected FFPE tissue and scores estimated from bulk FF tumor tissue. However, important discrepancies were revealed for luminal tumors, which may have consequences for treatment recommendations for these patients. Abstract The PAM50 gene expression subtypes and the associated risk of recurrence (ROR) score are used to predict the risk of recurrence and the benefits of adjuvant therapy in early-stage breast cancer. The Prosigna assay includes the PAM50 subtypes along with their clinicopathological features, and is approved for treatment recommendations for adjuvant hormonal therapy and chemotherapy in hormone-receptor-positive early breast cancer. The Prosigna test utilizes RNA extracted from macrodissected tumor cells obtained from formalin-fixed, paraffin-embedded (FFPE) tissue sections. However, RNA extracted from fresh-frozen (FF) bulk tissue without macrodissection is widely used for research purposes, and yields high-quality RNA for downstream analyses. To investigate the impact of the sample preparation approach on ROR scores, we analyzed 94 breast carcinomas included in an observational study that had available gene expression data from macrodissected FFPE tissue and FF bulk tumor tissue, along with the clinically approved Prosigna scores for the node-negative, hormone-receptor-positive, HER2-negative cases (n = 54). ROR scores were calculated in R; the resulting two sets of scores from FFPE and FF samples were compared, and treatment recommendations were evaluated. Overall, ROR scores calculated based on the macrodissected FFPE tissue were consistent with the Prosigna scores. However, analyses from bulk tissue yielded a higher proportion of cases classified as normal-like; these were samples with relatively low tumor cellularity, leading to lower ROR scores. When comparing ROR scores (low, intermediate, and high), discordant cases between the two preparation approaches were revealed among the luminal tumors; the recommended treatment would have changed in a minority of cases.

To investigate immune escape during breast tumor progression, we analyzed the composition of leukocytes in normal breast tissues, ductal carcinoma in situ (DCIS), and invasive ductal carcinomas (IDC). We found significant tissue and tumor subtype-specific differences in multiple cell types including T cells and neutrophils. Gene expression profiling of CD45+CD3+ T cells demonstrated a decrease in CD8+ signatures in IDCs. Immunofluorescence analysis showed fewer activated GZMB+CD8+ T cells in IDC than in DCIS, including in matched DCIS and recurrent IDC. T-cell receptor clonotype diversity was significantly higher in DCIS than in IDCs. Immune checkpoint protein TIGIT-expressing T cells were more frequent in DCIS, whereas high PD-L1 expression and amplification of CD274 (encoding PD-L1) was only detected in triple-negative IDCs. Coamplification of a 17q12 chemokine cluster with ERBB2 subdivided HER2+ breast tumors into immunologically and clinically distinct subtypes. Our results show coevolution of cancer cells and the immune microenvironment during tumor progression. Significance: The design of effective cancer immunotherapies requires the understanding of mechanisms underlying immune escape during tumor progression. Here we demonstrate a switch to a less active tumor immune environment during the in situ to invasive breast carcinoma transition, and identify immune regulators and genomic alterations that shape tumor evolution. Cancer Discov; 7(10); 1098–115. ©2017 AACR. See related commentary by Speiser and Verdeil, p. 1062. This article is highlighted in the In This Issue feature, p. 1047

Background Molecular profiling has recently been included in recommendations for decisions on adjuvant treatment in breast cancer (BrCa). However, the use of molecular profiling has not yet been widely established in all countries. Additional studies may give important information about the clinical relevance of the tests. Aims The study aims to discover the long term prognostic impact of PAM50 and PAM50 ROR score on survival for early BrCa pts according to treatment, with comparison to the routine clinical and histopathological parameters. Patients and methods Unselected early BrCa pts (n=651) from the Oslo Micromet project (n=920) having available FFPE primary tumor tissue were included in the current study. The pts were enrolled from 1995-1998. Follow up status is available for distant disease (median FU 7 years) and BrCa death (16.0-19.7 years after study inclusion). Clinical and histopathological parameters have been collected from the hospital records. FFPE tissue sections were macrodissected, RNA isolated from the dissected tumor tissue, followed by analysis of the PAM50 gene list on the Nanostring Platform. The samples were run in research mode and the raw data was sent to Nanostring (Seattle) for determination of the PAM50 subtype and ROR score. Results Of the 651 included pts, 323 did not receive any adjuvant systemic treatment (pT1pN0 patients), 161 tamoxifen only, the rest chemotherapy+/-tamoxifen. Twelve preoperatively treated pts were excluded from the analyses. Of the 639 remaining pts, PAM50 molecular profiling defined 52.3% as LumA, 26.8% LumB, 10.6% HER2enriched and 10.3% Basal. Multivariate analysis showed that the PAM50 intrinsic subtypes yielded prognostic information in addition to the established clinicopathological variables (pT, Grade, pN, age HR/HER2 subgroups, systemic treatment)(BCSS: HazardR vs LumA: 2.7 (95% CI 1.7-4.1) for LumB, 3.5 (1.8-6.8) for HER2enriched, 1.8 (0.8-4.2) for Basal). For the HR+HER2- pts, the risk classification by ROR score was an independent prognostic factor (BCSS: HazardR vs low risk: 3.1 (1.2-8.1) for intermediate, 6.6 (2.5-17.1) for high risk). In univariate analysis, the PAM50 intrinsic subtype classification separated clinical outcome both for all pts, for no adjuvant treated pts (both p Conclusions PAM50 subtype classification and ROR score improves classification of BrCa pts into prognostic groups, allowing more precise identification of future recurrence risk and improved basis for adjuvant treatment decisions. Citation Format: Naume B, Borgen E, Falk RS, Ohnstad HO, Lien TG, Aaserud M, Sveli MAT, Kyte JA, Kristensen V, Geitvik G, Schlichting E, Wist E, Sørlie T, Russnes H. Establishment of molecular profiling for individual treatment decisions in early breast cancer – Clinical impact of PAM50 and PAM50 risk of recurrence score after more than 16 years follow up [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr P2-05-16.

Increased expression of GLI1, the main Hedgehog signalling pathway effector, is related to unfavourable prognosis and progressive disease of certain breast cancer subtypes. We used conditional transgenic mice induced to overexpress GLI1 in the mammary epithelium either alone or in combination with deletion of one Trp53 allele to address the role of elevated GLI1 expression in breast tumour initiation and progression. Induced GLI1 expression facilitates mammary gland tumour formation and this was further increased upon heterozygous deletion of Trp53. The GLI1-induced primary tumours were of different murine molecular subtypes, including Normal-likeEx , Class8Ex , Claudin-LowEx and Erbb2-likeEx . The gene expression profiles of some of the tumours correlated well with the PAM50 subtypes for human breast cancer. Whole-exome sequencing revealed somatic mutation profiles with only little overlap between the primary tumours. Orthotopically serially transplanted GLI1-induced tumours maintained the main morphological characteristics of the primary tumours for ≥10 generations. Independent of Trp53 status and molecular subtype, the serially transplanted GLI1-induced tumours were able to grow both in the absence of transgenic GLI1 expression and in the presence of the GLI1 inhibitor GANT61. These data suggest that elevated GLI1 expression has a determinant role in tumour initiation; however, additional genetic events are required for tumour progression.

High mammographic density (MD) is associated with a 4–6 times increase in breast cancer risk. For post-menopausal women,MD often decreases over time, but little is known about the underlying biological mechanisms. MD reflects breast tissuecomposition, and may be associated with microenvironment subtypes previously identified in tumor-adjacent normal tissue. Currently, these subtypes have not been explored in normal breast tissue. We obtained biopsies from breasts of healthy women at two different time points several years apart and performed microarray gene expression analysis. At time point 1, 65 samples with both MD and gene expression were available. At time point 2, gene expression and MD data were available from 17 women, of which 11 also had gene expression data available from the first time point. We validated findings from our previous study; negative correlation between RBL1 and MD in post-menopausal women, indicating involvement of the TGFβ pathway. We also found that breast tissue samples from women with a large decrease in MD sustained higher expression of genes in the histone family H4. In addition, we explored the previously defined active and inactive microenvironment subtypes and demonstrated that normal breast samples of the active subtype had characteristics similar to the claudin-low breast cancer subtype. Breast biopsies from healthy women are challenging to obtain, but despite a limited sample size, we have identified possible mechanisms relevant for changes in breast biology and MD over time that may be of importance for breast cancer risk and tumor initiation. publishedVersion

Background: HER2 is a well established prognostic and predictive factor in invasive breast cancer. The role of HER2 in ductal breast carcinoma in situ (DCIS) is much debated and recent data have suggested that HER2 is mainly related to in situ recurrences. This contrasts the proposed role of HER2 in the progression from in situ to invasive cancer. Our aim was to study HER2 as a prognostic factor in a large population based cohort of DCIS. Methods: All 458 women diagnosed with a primary DCIS 1986-2004 in two Swedish regions were included and tissue microarrays constructed. Silver-enhanced in situ hybridisation (SISH) and immunohistochemistry (IHC) were used for detection of HER2 amplification and IHC expression. HER2 status and its relation to invasive breast cancer recurrence (IBCR) (ipsilateral or contralateral invasive events and regional or distant metastasis) and ipsilateral events (IBE) were studied. Kaplan-Meier survival analyses and Cox proportional hazards regression models were used. Adjustments were made for age, size, radiotherapy and ER status. Results: Mean follow up was 184 months. DCIS was screening detected in 75.5% of cases. Breast conserving surgery (BCS) was performed in 78.6% of whom 44.0% received postoperative radiotherapy. No women had hormonal or chemotherapy. A total of 106 IBCRs and 105 IBEs were identified. 54 IBEs were in situ and 51 invasive cancer. Eighteen women died from breast cancer and another 114 had died from other causes. 420 tumours could be classified using available SISH or IHC data; 132 were HER2 positive (31.4%) and 288 HER2 negative. SISH and IHC data were concordant in 296 of 332 (89.2%) available cases. HER2 positivity was related to size, grade and ER and PR negativity. The risk of IBCR was statistically significantly lower subsequent to a HER2 positive DCIS, hazard ratio (HR) 0.53 (95% CI 0.31-0.90), log rank p=0.01. However, the curves did not separate until after almost ten years. HRs after adjustments were similar to the crude analyses. The risk of any IBE was not statistically differently changed by HER2 status. In women undergoing BCS, HR was 1.18 (0.77-1.82), log rank p=0.44. But interestingly, divided by type of IBE, HER2-positivity showed an increased risk of in situ IBEs, HR 1.64 (0.92-2.90), log rank p=0.09 and, a decreased risk of invasive IBEs, HR 0.76 (0.39-1.50), log rank p=0.43. These correlations were however not statistically significant. Risk of invasive and local recurrence by HER2 status in a population based cohort of women with a primary DCIS HER2 positiveHER2 negativeInvasive Breast Cancer RecurrenceHR (95% CI)Reference 1.0All patients (n=420)(events=101)0.53 (0.31-0.90)RefIpsilateral new Breast Events (IBE) BCS (n=324) (events=94)1.18 (0.77-1.82)RefBCS, in situ IBEs (events=48)1.64 (0.92-2.90)RefBCS, invasive IBEs (events=46)0.76 (0.39-1.50)Ref Conclusions: In this long term follow-up DCIS cohort positive HER2 status in the primary DCIS predicted a statistically lower risk of IBCR. This effect was seen from ten years after primary surgery and onwards. The risk of an in situ IBE was increased and the risk of an invasive IBE was decreased if the primary DCIS was HER2 positive. All together, this challenges the role of HER2 as a driving force of the progression from in situ to invasive cancer. Citation Format: Signe Borgquist, Wenjing Zhou, Karin Jirström, Rose-Marie Amini, Thomas Sollie, Therese Sørlie, Salma Butt, Carl Blomqvist, Fredrik Wärnberg. The prognostic role of HER2 expression in ductal breast carcinoma in situ [abstract]. In: Proceedings of the Thirty-Seventh Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2014 Dec 9-13; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2015;75(9 Suppl):Abstract nr P6-13-03.

Featured Article: Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, et al. Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci USA 2001;98:10869–74.2 In our 2001 article highlighted here, my colleagues and I demonstrated for the first time the clinical impact of the main 5 “intrinsic” molecular subtypes of breast carcinomas. Shortly before, we had published the identification of unique “molecular portraits” of human breast tumors based solely on variation in gene expression patterns (1). The pervasive order and interpretation of these patterns indicated that the tumors could be classified into distinct biological groups based on overall molecular commonalities. The basal-like, ERBB2 (erb-b2 receptor tyrosine kinase 2)-positive,3 and normal-like subtypes all showed low or no expression of the estrogen receptor (ER) and other genes characteristic of luminal epithelial cells and had additional features specific for each of the 3 subtypes. For example, basal-like tumors showed high expression of the basal cytokeratins 5 and 17 and epidermal growth factor receptor and represented the predominant breast cancer subtype among BRCA1 (breast cancer 1) carriers (2). Luminal …

SummaryBreast cancer consists of at least five main molecular “intrinsic” subtypes that are reflected in both pre-invasive and invasive disease. Although previous studies have suggested that many of the molecular features of invasive breast cancer are established early, it is unclear what mechanisms drive progression and whether the mechanisms of progression are dependent or independent of subtype. We have generated mRNA, miRNA, and DNA copy-number profiles from a total of 59 in situ lesions and 85 invasive tumors in order to comprehensively identify those genes, signaling pathways, processes, and cell types that are involved in breast cancer progression. Our work provides evidence that there are molecular features associated with disease progression that are unique to the intrinsic subtypes. We additionally establish subtype-specific signatures that are able to identify a small proportion of pre-invasive tumors with expression profiles that resemble invasive carcinoma, indicating a higher likelihood of future disease progression.

We use an integrated approach to understand breast cancer heterogeneity by modeling mRNA, copy number alterations, microRNAs, and methylation in a pathway context utilizing the pathway recognition algorithm using data integration on genomic models (PARADIGM). We demonstrate that combining mRNA expression and DNA copy number classified the patients in groups that provide the best predictive value with respect to prognosis and identified key molecular and stromal signatures. A chronic inflammatory signature, which promotes the development and/or progression of various epithelial tumors, is uniformly present in all breast cancers. We further demonstrate that within the adaptive immune lineage, the strongest predictor of good outcome is the acquisition of a gene signature that favors a high T-helper 1 (Th1)/cytotoxic T-lymphocyte response at the expense of Th2-driven humoral immunity. Patients who have breast cancer with a basal HER2-negative molecular profile (PDGM2) are characterized by high expression of protumorigenic Th2/humoral-related genes (24–38%) and a low Th1/Th2 ratio. The luminal molecular subtypes are again differentiated by low or high FOXM1 and ERBB4 signaling. We show that the interleukin signaling profiles observed in invasive cancers are absent or weakly expressed in healthy tissue but already prominent in ductal carcinoma in situ, together with ECM and cell-cell adhesion regulating pathways. The most prominent difference between low and high mammographic density in healthy breast tissue by PARADIGM was that of STAT4 signaling. In conclusion, by means of a pathway-based modeling methodology (PARADIGM) integrating different layers of molecular data from whole-tumor samples, we demonstrate that we can stratify immune signatures that predict patient survival.

Angiogenesis is important for the growth and spread of malignant tumours (Carmeliet, 2003). Although vascular density has been correlated to survival in different malignancies (Weidner et al, 1991, 1992; Hlatky et al, 2002; Uzzan et al, 2004), angiogenesis markers have not been associated with response to chemotherapy or anti-angiogenesis treatment (Paulsen et al, 1997; Tynninen et al, 2002; Jubb et al, 2006). Glomeruloid microvascular proliferations (GMP), a novel angiogenesis marker, are focal aggregates of small vessels resembling a renal glomerulus (Pettersson et al, 2000; Sundberg et al, 2001). In mice, GMP has been induced by local injection of an adenovirus vector directing VEGF-A expression (Sundberg et al, 2001). A human parallel appears to be the POEMS syndrome, where increased VEGF-A is associated with glomeruloid haemangiomas (Tsai et al, 2001). In humans, GMP is a defining histological feature of glioblastoma multiforme (Wesseling et al, 1993; Schiffer et al, 1999) and a prognostic factor in several other tumours (Straume et al, 2002; Tanaka et al, 2003; Foulkes et al, 2004). In this study of locally advanced breast cancer, GMP showed a significant association with lack of treatment response and progressive disease following chemotherapy.

Reactivation of immune responses against cancer cells—immunotherapy—is one of the few cancer therapies that can successfully eliminate even metastatic disease in a relatively nontoxic manner. However, its success has been limited to a subset of patients. For example, in breast cancer only ~20% of triple-negative breast cancer (TNBC) patients benefit from anti-PDL1 therapy. One reason for this limited success can be that different tumors evade the immune system via different mechanisms, which suggests that they may respond to different types of immunotherapies. Epithelial cancer cells in ductal carcinoma in situ (DCIS) are physically separated from the tumor-infiltrating leukocytes by the myoepithelial cell layer and the basement membrane, whereas in invasive ductal carcinoma (IDC), the epithelial cancer cells are intermingled with leukocytes. Therefore, we hypothesize that the DCIS to IDC transition is a key step in tumor progression as cancer cells are under different selection pressures, and only those that can evade the immune system can continue tumor progression, hence shaping subsequent tumor evolution. To dissect the role of leukocytes in the DCIS to IDC transition, we began by analyzing the composition and molecular profiles of leukocytes, with special emphasis on T cells, in normal breast tissues, DCIS, and IDC. We found that the relative frequency of leukocytes increases during tumor progression but the CD8/CD4 T cell ratio decreases. In addition, the gene expression profile of CD45+CD3+ T cells is different in DCIS compared to those isolated from normal breast tissue and IDCs. We found that gene set signatures corresponding to CD8+ T cells and NKT cells were enriched over regulatory T-cell signatures in DCIS compared to IDC. This result suggested that DCIS had a more activated immune environment compared to IDC. We further examined T-cell activation by immunofluorescence (IF) analysis and found a higher percentage of activated GZMB+CD8+ T cells in DCIS compared to IDC including a set of matched DCIS and locally recurrent IDC. We also found that the TCR clonotype was more diverse in DCIS than in IDCs. Interestingly, we detected a few relatively frequent clones that were shared among different DCIS patients, one of which was previously shown to recognize a protein from the Epstein-Bar virus. In order to dissect mechanisms of immune evasion in IDC, we analyzed immune checkpoint genes and proteins by FISH and IF. We found that TIGIT+ T cells were slightly more frequent in DCIS than in IDC. In triple-negative IDC, there was high expression of PD-L1 in epithelial cells and in 3/10 cases amplification of CD274 (encoding PD-L1), whereas DCIS had lower expression of PD-L1 and no amplification of CD274. To further elucidate mechanisms of immune evasion, we explored the significance of a cluster of genes encoding several chemokines that are located in close proximity of ERBB2 (encoding HER2). When analyzing the HER2+ samples from the TCGA, we found that coamplification of the 17q12 chemokine cluster (CC) with ERBB2 was enriched in HER2+ER+ luminal-like tumors, whereas there was either no gain or loss of the cluster in the HER2+ER breast tumors. Interestingly, we found higher expression of both T-cell activation and exhaustion-related genes in tumors that lack CC gain. Moreover, when assessing a cohort of HER2+ samples by multicolor FISH and IF, we found an inverse correlation between CC amplification and activation of CD8+ T cells. There was no correlation between CC amplification and recruitment of macrophages or myeloid-derived suppressor cells. Overall our results show coevolution of cancer cells and the immune microenvironment during tumor progression. Citation Format: Carlos R. Gil del Alcazar, SungJin Huh, Muhammad B. Ekram, Anne Trinh, Lin L. Liu, Francisco Beca, Zi Xiaoyuan, Misuk Kwak, Helga Bergholtz, Ying Su, Lina Ding, Hege G. Russnes, Andrea L. Richardson, Kirsten Babski, Elizabeth Min Hui Kim, Charles H. McDonnell, III, Jon Wagner, Ron Rowberry, Gordon J. Freeman, Deborah Dillon, Therese Sorlie, Lisa M. Coussens, Judy E. Garber, Rong Fan, Kristie Bobolis, D. Craig Allred, Joon Jeong, So Yeon Park, Franziska Michor, Kornelia Polyak. Characterization of the immune environment in the in situ to invasive breast carcinoma transition [abstract]. In: Proceedings of the AACR Special Conference: Advances in Breast Cancer Research; 2017 Oct 7-10; Hollywood, CA. Philadelphia (PA): AACR; Mol Cancer Res 2018;16(8_Suppl):Abstract nr A21.

This abstract was not presented at the symposium.

Glycosylation is the stepwise procedure of covalent attachment of oligosaccharide chains to proteins or lipids, and alterations in this process have been associated with malignant transformation. Simultaneous analysis of the expression of all glycan-related genes clearly gives the advantage of enabling a comprehensive view of the genetic background of the glycobiological changes in cancer cells. Studies focusing on the expression of the whole glycome have now become possible, which prompted us to review the present knowledge on glycosylation in relation to breast cancer diagnosis and progression, in the light of available expression data from tumors and breast tissue of healthy individuals. We used various data resources to select a set of 419 functionally relevant genes involved in synthesis, degradation and binding of N-linked and O-linked glycans, Lewis antigens, glycosaminoglycans (chondroitin, heparin and keratan sulfate in addition to hyaluronan) and glycosphingolipids. Such glycans are involved in a number of processes relevant to carcinogenesis, including regulation of growth factors/growth factor receptors, cell-cell adhesion and motility as well as immune system modulation. Expression analysis of these glycan-related genes revealed that mRNA levels for many of them differ significantly between normal and malignant breast tissue. An associative analysis of these genes in the context of current knowledge of their function in protein glycosylation and connection(s) to cancer indicated that synthesis, degradation and adhesion mediated by glycans may be altered drastically in mammary carcinomas. Although further analysis is needed to assess how changes in mRNA levels of glycan genes influence a cell's glycome and the precise role that such altered glycan structures play in the pathogenesis of the disease, lessons drawn from this study may help in determining directions for future research in the rapidly-developing field of glycobiology.

Genetic factors are thought to play a role in resistance towards chemotherapeutic agents such as 5-fluorouracil (5-FU). Approximately 30 genes are directly or indirectly involved in 5-FU metabolism, and genetic variation in any of these may contribute to anti-tumor response. Polymorphisms in these genes were analyzed in relation to tumoral mRNA levels of 5-FU metabolizing genes, response to chemotherapy and survival. A total of 21 genetic variants were studied in 35 breast cancer patients treated with 5-FluoroUracil, mitomycin (FUMI) and in a similar cohort of 90 doxorubicin treated breast cancer patients. Genotype distributions were compared using 109 healthy controls. No significant association was found between any polymorphisms and response to chemotherapy as measured by shrinkage of tumor. However, carriers of 3 copies of the TYMS 5'UTR repeat had shorter survival than noncarriers (p = .048) in the FUMI treatment group, but not in the doxorubicin treated group. Carriers of 3 copies of the repeat were also more frequently observed in both patients groups than in healthy controls (p = .034). Several highly significant associations were observed between genotypes and expression levels of 5-FU metabolizing genes. A SNP in codon 72 of TP53 was revealed to be a key regulator of 5-FU metabolizing genes such as DHFR and MTHFR, constituting 50% of all significant associations observed after FUMI therapy. These data suggest that 3 copies of the TYMS 5'UTR repeat may give a treatment specific reduced survival in breast cancer patients, and that TP53 may have a direct, allele specific, role in 5-FU mediated response.

Expression profiles of primary breast tumors were investigated in relation to disseminated tumor cells (DTCs) in bone marrow (BM) in order to increase our understanding of the dissemination process. Tumors were classified into five pre-defined molecular subtypes, and presence of DTC identified (at median 85 months follow-up) a subgroup of luminal A patients with particular poor outcome ( p ¼ 0.008). This was not apparent for other tumor subtypes. Gene expression profiles associated with DTC and with systemic relapse for luminal A patients were identified. This study suggests that DTC in BM differentially distinguishes clinical outcome in patients with luminal A type tumors and that DTC-associated gene expression analysis may identify genes of potential importance in tumor dissemination. a 2007 Federation of European Biochemical Societies.

Alterations in glycan structures are early signs of malignancy and have recently been proposed to be in part a driving force behind malignant transformation. Here, we explore whether differences in expression of genes related to the process of glycosylation exist between breast carcinoma subtypes – and look for their association to clinical parameters. Five expression datasets of 454 invasive breast carcinomas, 31 ductal carcinomas in situ (DCIS), and 79 non-malignant breast tissue samples were analysed. Results were validated in 1960 breast carcinomas. 419 genes encoding glycosylation-related proteins were selected. The DCIS samples appeared expression-wise similar to carcinomas, showing altered gene expression related to glycosaminoglycans (GAGs) and N-glycans when compared to non-malignant samples. In-situ lesions with different aggressiveness potentials demonstrated changes in glycosaminoglycan sulfation and adhesion proteins. Subtype-specific expression patterns revealed down-regulation of genes encoding glycan-binding proteins in the luminal A and B subtypes. Clustering basal-like samples using a consensus list of genes differentially expressed across discovery datasets produced two clusters with significantly differing prognosis in the validation dataset. Finally, our analyses suggest that glycolipids may play an important role in carcinogenesis of breast tumors – as demonstrated by association of B3GNT5 and UGCG genes to patient survival. In conclusion, most glycan-specific changes occur early in the carcinogenic process. We have identified glycan-related alterations specific to breast cancer subtypes including a prognostic signature for two basal-like subgroups. Future research in this area may potentially lead to markers for better prognostication and treatment stratification of breast cancer patients.

Neoadjuvant treatment offers an opportunity to correlate molecular variables to treatment response and to explore mechanisms of drug resistance in vivo. Here, we present a statistical analysis of large-scale gene expression patterns and their relationship to response following neoadjuvant chemotherapy in locally advanced breast cancers. We analyzed cDNA expression data from 81 tumors from two patient series, one treated with doxorubicin alone (51) and the other treated with 5-fluorouracil and mitomycin (30), and both were previously studied for correlations between TP53 status and response to therapy. We observed a low frequency of progressive disease within the luminal A subtype from both series (2 of 36 versus 13 of 45 patients; P = 0.0089) and a high frequency of progressive disease among patients with luminal B type tumors treated with doxorubicin (5 of 8 patients; P = 0.0078); however, aside from these two observations, no other consistent associations between response to chemotherapy and tumor subtype were observed. These specific associations could possibly be explained by covariance with TP53 mutation status, which also correlated with tumor subtype. Using supervised analysis, we could not uncover a gene profile that could reliably (>70% accuracy and specificity) predict response to either treatment regimen. [Mol Cancer Ther 2006;5(11):2914–8]

Breast cancer is a leading cause of cancer-death among women, where the clinicopathological features of tumors are used to prognosticate and guide therapy. DNA copy number alterations (CNAs), which occur frequently in breast cancer and define key pathogenetic events, are also potentially useful prognostic or predictive factors. Here, we report a genome-wide array-based comparative genomic hybridization (array CGH) survey of CNAs in 89 breast tumors from a patient cohort with locally advanced disease. Statistical analysis links distinct cytoband loci harboring CNAs to specific clinicopathological parameters, including tumor grade, estrogen receptor status, presence of TP53 mutation, and overall survival. Notably, distinct spectra of CNAs also underlie the different subtypes of breast cancer recently defined by expression-profiling, implying these subtypes develop along distinct genetic pathways. In addition, higher numbers of gains/losses are associated with the "basal-like" tumor subtype, while high-level DNA amplification is more frequent in "luminal-B" subtype tumors, suggesting also that distinct mechanisms of genomic instability might underlie their pathogenesis. The identified CNAs may provide a basis for improved patient prognostication, as well as a starting point to define important genes to further our understanding of the pathobiology of breast cancer. This article contains Supplementary Material available at http://www.interscience.wiley.com/jpages/1045-2257/suppmat

Interactions between luminal epithelial cells and their surrounding microenvironment govern the normal development and function of the mammary gland. Estradiol plays a key role in abnormal intracellular signaling, which contributes to the development and progression of breast tumors. The present article summarizes the results from a microarray whole genome gene expression analysis as well as a quantitative analysis of the mRNA expression of members of the estradiol metabolic and signaling pathways in the tumors of postmenopausal breast cancer patients. The analysis of the variation in whole genome gene expression resulted in a tumor classification comprising several distinct groups with distinct expression of the estrogen receptor (ER). The parallel study on the expression of only nine mRNA transcripts of members of the estradiol pathways resulted in two main clusters, representing ER− and ER tumors. The mRNA expression of the estradiol-metabolizing enzymes did not follow the expression of the ER in all cases, leading to the recognition of several further subclasses of tumors. When the tumor classes obtained by whole genome gene expression analysis were compared with those obtained by independent quantitation of the estradiol-metabolizing enzymes, a statistically significant association between both classification groups was observed. These findings point to a possible association between development of a tumor with a particular expression profile and its capacity to synthesize estradiol as measured by the expression of the transcripts for the necessary key enzymes. Further, whole genome expression patterns were studied in 12 patients treated with anastrozole. Using significance analysis of microarrays, we identified 298 genes significantly differently expressed between partial response and progressive disease groups.

Gene mutations play a critical role in cancer development and progression, and their identification offers possibilities for accurate diagnostics and therapeutic targeting. Finding genes undergoing mutations is challenging and slow, even in the post-genomic era. A new approach was recently developed by Noensie and Dietz to prioritize and focus the search, making use of nonsense-mediated mRNA decay (NMD) inhibition and microarray analysis (NMD microarrays) in the identification of transcripts containing nonsense mutations. We combined NMD microarrays with array-based CGH (comparative genomic hybridization) in order to identify inactivation of tumor suppressor genes in cancer. Such a “mutatomics” screening of prostate cancer cell lines led to the identification of inactivating mutations in the EPHB2 gene. Up to 8% of metastatic uncultured prostate cancers also showed mutations of this gene whose loss of function may confer loss of tissue architecture. NMD microarray analysis could turn out to be a powerful research method to identify novel mutated genes in cancer cell lines, providing targets that could then be further investigated for their clinical relevance and therapeutic potential.

Recent microarray studies have identified distinct subtypes of breast tumors that arise from different cell types and that show statistically significant differences in patient outcome. To gain insight into these differences, we identified in vitro and in vivo changes in gene expression induced by chemotherapeutics. We treated two cell lines derived from basal epithelium (immortalized human mammary epithelial cells) and two lines derived from luminal epithelium (MCF-7 and ZR-75–1) with chemotherapeutics used in the treatment of breast cancer and assayed for changes in gene expression using DNA microarrays. Treatment doses for doxorubicin and 5-fluorouracil were selected to cause comparable cytotoxicity across all four cell lines. The dominant expression response in each of the cell lines was a general stress response; however, distinct expression patterns were observed. Both cell types induced DNA damage-response genes such as p21waf1, but the response in the luminal cells showed higher fold changes and included more p53-regulated genes. Luminal cell lines repressed a large number of cell cycle-regulated genes and other genes involved in cellular proliferation, whereas the basal cell lines did not. Instead, the basal cell lines repressed genes that were involved in differentiation. These in vitro responses were compared with expression responses in breast tumors sampled before and after treatment with doxorubicin or 5-fluorouracil/mitomycin C. The in vivo data corroborated the cell-type-specific responses to chemotherapeutics observed in vitro, including the induction of p21waf1. Similarities between in vivo and in vitro responses help to identify important response mechanisms to chemotherapeutics.

Immunocompetent transgenic mouse models have for decades served as valuable tools to address the effect of oncogenes and tumor suppressors. Immunodeficient mice have been used to establish patient derived xenograft (PDX) models harboring human tumors. These types of models have been used to study cancer initiation and progression as well as preclinical evaluation of anticancer drugs. One major limiting factor for both these model systems is the immune system; the PDX models lack an immune system and the immune system in immunocompetent mouse models display differences compared to the human immune system. Thus, in cancer research there is a need for preclinical models to study the influence of human immune cells on tumor progression and response to cancer therapies. We have established and optimized protocols for intravenous (IV) or intraperitoneal (IP) injections of isolated human PBMCs, to generate humanized mice harboring human immune cells. Our protocols do not require irradiation or busulfan pretreatment of the animals. Flow cytometry analyses showed that mainly T helper cells, CD4+, and cytotoxic T cells, CD8+, were present in the humanized mice. Immunohistochemistry (IHC) analyses of humanized PDX cancer model mice showed that human lymphocytes were present in the tumor periphery and some very few cells within the tumor, in addition to the human leukocytes in the liver and spleen. The growth rates of the PDX tumors were not affected by the humanization. Our protocols require minimal preparations of the animals and generate humanized PDX mice harboring human lymphocytes locating in close proximity of, as well as inside the PDX tumors. The model systems are suitable for preclinical studies of human, adaptive immune responses in tumor progression and cancer therapies. Citation Format: Jens Henrik Norum, Dag Josefsen, Gunnar Kvalheim, Olav Engebråten, Therese Sørlie, Gunhild M. Mælandsmo. Humanization of patient derived xenograft (PDX) cancer model mice with peripheral blood mononuclear cells (PBMCs) [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1660. doi:10.1158/1538-7445.AM2017-1660

Purpose: To identify genes predicting benefit of radiotherapy in patients with high-risk breast cancer treated with systemic therapy and randomized to receive or not receive postmastectomy radiotherapy (PMRT). Experimental Design: The study was based on the Danish Breast Cancer Cooperative Group (DBCG82bc) cohort. Gene-expression analysis was performed in a training set of frozen tumor tissue from 191 patients. Genes were identified through the Lasso method with the endpoint being locoregional recurrence (LRR). A weighted gene-expression index (DBCG-RT profile) was calculated and transferred to quantitative real-time PCR (qRT-PCR) in corresponding formalin-fixed, paraffin-embedded (FFPE) samples, before validation in FFPE from 112 additional patients. Results: Seven genes were identified, and the derived DBCG-RT profile divided the 191 patients into “high LRR risk” and “low LRR risk” groups. PMRT significantly reduced risk of LRR in “high LRR risk” patients, whereas “low LRR risk” patients showed no additional reduction in LRR rate. Technical transfer of the DBCG-RT profile to FFPE/qRT-PCR was successful, and the predictive impact was successfully validated in another 112 patients. Conclusions: A DBCG-RT gene profile was identified and validated, identifying patients with very low risk of LRR and no benefit from PMRT. The profile may provide a method to individualize treatment with PMRT. Clin Cancer Res; 20(20); 5272–80. ©2014 AACR.

Introduction Dysregulated choline metabolism is a well-known feature of breast cancer, but the underlying mechanisms are not fully understood. In this study, the metabolomic and transcriptomic characteristics of a large panel of human breast cancer xenograft models were mapped, with focus on choline metabolism. Methods Tumor specimens from 34 patient-derived xenograft models were collected and divided in two. One part was examined using high-resolution magic angle spinning (HR-MAS) MR spectroscopy while another part was analyzed using gene expression microarrays. Expression data of genes encoding proteins in the choline metabolism pathway were analyzed and correlated to the levels of choline (Cho), phosphocholine (PCho) and glycerophosphocholine (GPC) using Pearson’s correlation analysis. For comparison purposes, metabolic and gene expression data were collected from human breast tumors belonging to corresponding molecular subgroups. Results Most of the xenograft models were classified as basal-like (N = 19) or luminal B (N = 7). These two subgroups showed significantly different choline metabolic and gene expression profiles. The luminal B xenografts were characterized by a high PCho/GPC ratio while the basal-like xenografts were characterized by highly variable PCho/GPC ratio. Also, Cho, PCho and GPC levels were correlated to expression of several genes encoding proteins in the choline metabolism pathway, including choline kinase alpha (CHKA) and glycerophosphodiester phosphodiesterase domain containing 5 (GDPD5). These characteristics were similar to those found in human tumor samples. Conclusion The higher PCho/GPC ratio found in luminal B compared with most basal-like breast cancer xenograft models and human tissue samples do not correspond to results observed from in vitro studies. It is likely that microenvironmental factors play a role in the in vivo regulation of choline metabolism. Cho, PCho and GPC were correlated to different choline pathway-encoding genes in luminal B compared with basal-like xenografts, suggesting that regulation of choline metabolism may vary between different breast cancer subgroups. The concordance between the metabolic and gene expression profiles from xenograft models with breast cancer tissue samples from patients indicates that these xenografts are representative models of human breast cancer and represent relevant models to study tumor metabolism in vivo. © 2014 Grinde et al.; licensee BioMed Central Ltd. This is an open access article distributed under the terms of the CreativeCommons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The purpose of this study was to classify breast carcinomas based on variations in gene expression patterns derived from cDNA microarrays and to correlate tumor characteristics to clinical outcome. A total of 85 cDNA microarray experiments representing 78 cancers, three fibroadenomas, and four normal breast tissues were analyzed by hierarchical clustering. As reported previously, the cancers could be classified into a basal epithelial-like group, an ERBB2 -overexpressing group and a normal breast-like group based on variations in gene expression. A novel finding was that the previously characterized luminal epithelial/estrogen receptor-positive group could be divided into at least two subgroups, each with a distinctive expression profile. These subtypes proved to be reasonably robust by clustering using two different gene sets: first, a set of 456 cDNA clones previously selected to reflect intrinsic properties of the tumors and, second, a gene set that highly correlated with patient outcome. Survival analyses on a subcohort of patients with locally advanced breast cancer uniformly treated in a prospective study showed significantly different outcomes for the patients belonging to the various groups, including a poor prognosis for the basal-like subtype and a significant difference in outcome for the two estrogen receptor-positive groups.

Malignant peripheral nerve sheath tumors (MPNSTs) are frequently associated with the disease neurofibromatosis type 1. Only few recurrent cytogenetic changes have been reported, including rearrangements of the short arm of chromosome 9. By fluorescence in situ hybridization with a centromere 9 probe, and by allelic imbalance studies with seven 9p21-23 markers in nine familial and three sporadic MPNSTs, we found interstitial deletions that supported CDKN2A as a possible target gene. Nine MPNSTs showed aberrations of CDKN2A by Southern blot analyses, and in four of these, expression of CDKN2A could not be detected by Northern blot analysis. No mutations of CDKN2A were identified by sequencing of the coding region, and gene inactivation by promoter methylation was not found. In the 9p allelic imbalance studies, a novel allele was detected at one locus in one tumor. Analyses of additional markers (n = 8) excluded mismatch repair deficiency as an important mechanism in the genesis of these tumors. The tumors were analyzed further for alterations in other candidate cell cycle-associated genes. In total, 11/12 MPNSTs showed DNA changes in one or more of the genes CDKN2A, CDKN2B, RB1, CDK4, MDM2, and CCND2. The present study suggests that disruption of the pRB pathway is common in MPNST, and that dose reduction of CDKN2A is particularly frequent and contributes to MPNST development. Genes Chromosomes Cancer 26:151-160, 1999.

Exposure to aromatic amines is considered a major risk factor for the development of bladder cancer. In this study, we have analysed the pattern of point mutations in several tumour genes in 21 cases of bladder cancer arising among western European workers exposed to aromatic amines in an attempt to determine whether this exposure may be associated with a unique spectrum of mutations. Of the four genes analysed (p53, p16MTS1, p21WAF1 and H-ras), only p53 showed a high frequency of mutations (in 8 out of 21 cases, 38%). Two mutations were found in p16, one in H-ras and none in p21 exon 3. All mutations were at G:C base pairs, mostly at non-CpG residues. This spectrum of mutations, which is highly suggestive of an involvement of exogenous carcinogens, is however identical to the spectrum of p53 mutations detected in bladder cancers of the general population. In exposed workers, p53 mutations were associated with tumour grade and with high occupational and tobacco exposure. Taken together, our data suggest that the same carcinogens may be responsible for the development of bladder cancers in workers exposed to aromatic amines and in the general population.

Germline mutations in the BRCA1 gene predispose to breast and ovarian cancer. An estimated 45% of families with multiple breast cancer cases and more than 80% of breast-ovarian cancer families are linked to BRCA1. Mutation analyses by collaborative laboratories have revealed around 460 distinct BRCA1 sequence alterations, mostly germline mutations from familial cases. The majority of these alterations were nonsense and frame-shift mutations. In the present study, breast tumors of both sporadic and familial origin were investigated for allelic imbalance (AI) at the BRCA1 locus. AI was observed in 52% of the sporadic cases and in 17% of the familial cases. Furthermore, 104 breast carcinomas from patients with sporadic disease and 77 patients with positive family histories of breast and/or ovarian cancer were examined for translation-terminating mutations in exon 11 of the BRCA1 gene using the protein truncation test (PTT). No somatic mutations were detected in any of the tumors analysed, and only one BRCA1 mutation carrier was found among the familial cases. The result of this study gives no indication that truncating somatic mutations in exon 11 of BRCA1 play a major role in the tumorigenesis of the breast. Furthermore, the frequency of such mutation carriers in breast cancer populations with weak family histories of breast and/or ovarian cancer seems to be low.

For a panel of cancer related proteins, the aim was to shed light on which molecular level the expression of each protein was mainly regulated in breast tumors, and to investigate whether differences in regulation were reflected in different molecular subtypes. DNA, mRNA and protein lysates from 251 breast tumor specimens were analyzed using appropriate microarray technologies. Data from all three levels were available for 52 proteins selected for their known involvement in cancer, primarily through the PI3K/Akt pathway. For every protein, in cis Spearman rank correlations between the three molecular levels were calculated across all samples and within each intrinsic gene expression subtype, enabling 63 comparisons altogether due to multiple gene probes matching to single proteins. Subtype‐specific relationships between the three molecular levels were studied by calculating the variance of subtype‐specific correlation and differences between overall and average subtype‐specific correlation. The findings were validated in an external dataset comprising 703 breast tumor specimens. The proteins were sorted into four groups based on the calculated rank correlation values between the three molecular levels. Group A consisted of eight proteins with significant correlation between DNA copy number levels and mRNA expression, and between mRNA expression and protein expression (Bonferroni adjusted p

BACKGROUND: Different molecular subtypes of breast cancer have been identified based on gene expression profiling. Treatment suggestions based on an approximation of these subtypes by immunohistochemical criteria have been published by the St Gallen international expert consensus panel. Ductal carcinoma in situ (DCIS) can be classified into the same molecular subtypes. Our aim was to study the relation between these newly defined subtypes and prognosis in DCIS. METHODS: TMA including 458 women from a population-based cohort with DCIS diagnosed 1986-2004 was used. Stainings for ER, PR, HER2 and Ki67 were used to classify the surrogate molecular subtypes according to the St Gallen criteria from 2011. The associations with prognosis were examined using Kaplan-Meier analyses and Cox proportional hazards regression models. RESULTS: Surrogate molecular subtyping could be done in 381 cases. Mean follow up was 164 months. Of the classified DCIS 186 were Luminal A (48.8%), 33 Luminal B/HER2- (8.7%), 74 Luminal B/HER2+ (17.4%), 61 HER2+/ER- (16.0%) and 27 Triple Negative (7.1%). One hundred and two women had a local recurrence of which 58 were invasive. Twenty-two women had generalised disease, 8 without a prior local recurrence. We could not find a prognostic significance of the molecular subtypes other than a higher risk of developing breast cancer after more than 10 years of follow-up among women with a Triple Negative DCIS (OR 3.2; 95% CI 1.1-9.8). CONCLUSIONS: The results from this large population-based cohort, with long-term follow up failed to demonstrate a prognostic value for the surrogate molecular subtyping of DCIS using the St Gallen criteria up to ten years after diagnosis. More than ten years after diagnosis Triple Negative DCIS had an elevated risk of recurrence.

Tankyrase enzymes play crucial roles in the Wnt signaling pathway, involved in the regulation of intestinal stem cells and tissue homeostasis. Aberrant Wnt signaling is known to cause intestinal cancers and the tankyrase inhibitor G007-LK has previously been shown to inhibit tumor growth in an APC-mutant colorectal cancer xenograft model. We have used in vivo models to address the effect of G007-LK on small intestine tissue homeostasis. The mice were treated for up to 3 weeks with the tankyrase inhibitor, without affecting the body mass or observed clinical condition of the mice. H&E staining of fixed small intestine tissue sections showed no significant differences in the general tissue morphology between control and G007-LK treated mice. In the presence of G007-LK, lineage tracing from the Lgr5+ intestinal stem cells was significantly reduced compared to control mice. When the G007-LK treatment was stopped two days prior to initiation of lineage tracing, there were no significant difference in the number of traced cells between the mice pretreated with G007-LK and untreated control group. Immunohistochemistry (IHC) staining for the proliferation marker Ki67 showed a slight reduction in the number of positive cells in the small intestine of the G007-LK treated compared to control mice. As expected, in the crypt base cells there was a reduction of nuclei staining positive for β-catenin in the G007-LK treated compared to control mice, suggesting repressed Wnt/β-catenin signaling. Taken together, our data show that the tankyrase inhibitor G007-LK is well tolerated by the mice. The tankyrase inhibitor reduced the proliferative activity of the small intestine Lgr5+ stem cells, without altering the general morphology or function of the tissue. Citation Format: Jens Henrik Norum, Ellen Skarpen, Andreas Brech, Raoul Kuiper, Jo Waaler, Stefan Krauss, Therese Sørlie. The tankyrase inhibitor G007-LK inhibits small intestine Lgr5+ stem cell proliferation without altering tissue morphology. [abstract]. In: Proceedings of the Fourth AACR International Conference on Frontiers in Basic Cancer Research; 2015 Oct 23-26; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2016;76(3 Suppl):Abstract nr A18.

Estrogen receptor is a key driver in breast cancer and is expressed in about 75% of breast tumors. ER positive tumors are susceptible to endocrine therapies; however, the major obstacle for curative treatment is recurrence due to resistance to anti-estrogens. Endocrine therapies may induce a selective pressure promoting growth of estrogen independent cell subclones. Our aim was to reveal molecular changes occurring in tumors in response to anti-estrogen treatment, and to identify subpopulations of cells able to withstand anti-estrogen treatment. A luminal-like estrogen-dependent orthotopically growing xenograft model was treated with fulvestrant, or exposed to estrogen deprivation. The effect of ER-signaling inhibition was analyzed using quantitative mass spectrometry (MS) -based proteomic analysis and high resolution magic angle spinning magnetic resonance spectroscopy (HR MAS MRS). Cell surface marker expression (CD24 and SSEA-4) was monitored by flow cytometry, allowing detailed comparison of protein expression between intratumor cell subpopulations. We found that both modes of anti-estrogen therapy restrained tumor growth and induced expression of enzymes involved in TCA cycle, oxidative phosphorylation and fatty acid beta-oxidation. This was accompanied by changes in levels of specific metabolites indicative of a possible reprogramming of cell metabolism and utilization of oxidative phosphorylation in preference to aerobic glycolysis (decrease in Warburg effect). Furthermore, anti-estrogen treatment seemed to have selective effects on intratumor cell subpopulations, specified by expression of the markers CD24 and SSEA-4. More specifically, highly tumorigenic CD24low/SSEA-4low (dbl. low) cells were eliminated and the seemingly more benign CD24high/SSEA-4high (dbl. high) cells were enriched in the residual tumor. When comparing the proteome in dbl. low verus dbl. high cells sorted from untreated tumors, metabolism was one of the most differentially enriched processes. Enzymes involved in glycolysis, TCA cycle, respiratory electron transport chain and fatty acid were more abundant in the dbl. high subpopulation. These results suggest that cancer cells may reprogram their metabolism in response to anti-estrogen therapy to support a less estrogen-dependent phenotype. Moreover, subpopulations of cells with different metabolism may exist within the growing tumor, and these may respond differently to anti-estrogen treatment. Citation Format: Nirma Skrbo, Ufuk Kirik, Alexandr Kristian, Paolo Cifani, Linn Antberg, Siver A. Moestue, Olav Engebraaten, Gunhild M. Mælandsmo, Kristin Andersen, Peter James, Therese Sørlie. Protein expression analysis of intratumor heterogeneity in a luminal-like breast cancer xenograft. [abstract]. In: Proceedings of the AACR Special Conference on Advances in Breast Cancer Research; Oct 17-20, 2015; Bellevue, WA. Philadelphia (PA): AACR; Mol Cancer Res 2016;14(2_Suppl):Abstract nr A36.

The majority of human breast cancers exhibit luminal epithelial differentiation. However, most aggressive behavior, including invasion and purported cancer stem cell activity, are considered characteristics of basal-like cells. We asked the following questions: Must luminal-like breast cancer cells become basal-like to initiate tumors or to invade? Could luminally differentiated cells within a basally initiated hierarchy also be tumorigenic? To answer these questions, we used rare and mutually exclusive lineage markers to isolate subsets of luminal-like and basal-like cells from human breast tumors. We enriched for populations with or without prominent basal-like traits from individual tumors or single cell cloning from cell lines and recovered cells with a luminal-like phenotype. Tumor cells with basal-like traits mimicked phenotypic and functional behavior associated with stem cells assessed by gene expression, mammosphere formation and lineage markers. Luminal-like cells without basal-like traits, surprisingly, were fully capable of initiating invasive tumors in NOD SCID gamma (NSG) mice. In fact, these phenotypically pure luminal-like cells generated larger and more invasive tumors than their basal-like counterparts. The tumorigenicity and invasive potential of the luminal-like cancer cells relied strongly on the expression of the gene GCNT1 , which encodes a key glycosyltransferase controlling O-glycan branching. These findings demonstrate that basal-like cells, as defined currently, are not a requirement for breast tumor aggressiveness, and that within a single tumor there are multiple “stem-like” cells with tumorigenic potential casting some doubt on the hypothesis of hierarchical or differentiative loss of tumorigenicity.

Source at https://doi.org/10.1016/j.molonc.2012.10.005. Accepted manuscript version, licensed CC BY-NC-ND 4.0. Antiangiogenic therapy with bevacizumab has shown varying results in breast cancer clinical trials. Identifying robust biomarkers for selecting patients who may benefit from such treatment and for monitoring response is important for the future use of bevacizumab. Two established xenograft models representing basal‐like and luminal‐like breast cancer were used to study bevacizumab treatment response on the metabolic and gene expression levels. Tumor samples were obtained from mice treated with bevacizumab, doxorubicin or a combination of the two drugs, and high resolution magic angle spinning magnetic resonance spectroscopy and gene expression microarray analysis was performed. Combination treatment with bevacizumab showed the strongest growth inhibiting effect in basal‐like tumors, and this was reflected by a significant change in the metabolomic and transcriptomic profiles. In the luminal‐like xenografts, addition of bevacizumab did not improve the effect of doxorubicin. On the global transcriptomic level, the largest gene expression changes were observed for the most efficient treatment in both models. Glycerophosphocholine showed opposite response in the treated xenografts compared with untreated controls; lower in basal‐like and higher in luminal‐like tumors. Comparing combination therapy with doxorubicin monotherapy in basal‐like xenografts, 14 genes showed significant differential expression, including very low density lipoprotein receptor (VLDLR) and hemoglobin, theta 1 (HBQ1). Bevacizumab‐treated tumors were associated with a more hypoxic phenotype, while no evidence was found for associations between bevacizumab treatment and vascular invasion or tumor grade. This study underlines the importance of characterizing biological differences between subtypes of breast cancer to identify personalized biomarkers for improved patient stratification and evaluation of response to therapy.

Time to freezing tumor tissue for RNA expression analysis will always vary to some extent. To evaluate the effect of ischemia time, tumor tissue from ten breast cancer patients was collected and aliquots of tissue were snap frozen at different time points after surgery (0, 0.5, 1, 3 and 6 h). Using miRNA and mRNA expression microarrays and statistical analysis, 56 miRNAs and 1788 mRNAs were found to be significantly altered with ischemia time up to six hours. Several of the 56 miRNAs have been reported to play a role in cancer, such as hsa-miR-663 and hsa-miR-125a-3p. Known stress response genes such as GADD45B, JUND and FOSB were among the mRNAs most significantly affected by time to freezing. A novel statistical method for identification of consistently correlated miRNA-mRNA pairs and miRNA-associated biological processes in time course data is presented. Application of this method revealed that several miRNAs, including hsa-miR-1228, hsa-miR-1225-5p and hsa-miR-574-5p, were associated through their correlation to mRNAs to biological processes such as "response to stimulus" and "stress response". These miRNAs also showed enrichment of predicted targets among either their positively or negatively correlated mRNAs. The induced miRNAs may play both direct and indirect roles in biological responses. Caution should be taken when the miRNAs and mRNAs reported to be affected by ischemia time are included in a prognostic or predictive signature.

Ductal carcinoma in situ (DCIS) is a non‐invasive form of breast cancer where cells restricted to the ducts exhibit an atypical phenotype. Some DCIS lesions are believed to rapidly transit to invasive ductal carcinomas (IDCs), while others remain unchanged. Existing classification systems for DCIS fail to identify those lesions that transit to IDC. We studied gene expression patterns of 31 pure DCIS, 36 pure invasive cancers and 42 cases of mixed diagnosis (invasive cancer with an in situ component) using Agilent Whole Human Genome Oligo Microarrays 44k. Six normal breast tissue samples were also included as controls. qRT‐PCR was used for validation. All DCIS and invasive samples could be classified into the “intrinsic” molecular subtypes defined for invasive breast cancer. Hierarchical clustering establishes that samples group by intrinsic subtype, and not by diagnosis. We observed heterogeneity in the transcriptomes among DCIS of high histological grade and identified a distinct subgroup containing seven of the 31 DCIS samples with gene expression characteristics more similar to advanced tumours. A set of genes independent of grade, ER‐status and HER2‐status was identified by logistic regression that univariately classified a sample as belonging to this distinct DCIS subgroup. qRT‐PCR of single markers clearly separated this DCIS subgroup from the other DCIS, and contains samples from several histopathological and intrinsic molecular subtypes. The genes that differentiate between these two types of DCIS suggest several processes related to the re‐organisation of the microenvironment. This raises interesting possibilities for identification of DCIS lesions both with and without invasive characteristics, which potentially could be used in clinical assessment of a woman's risk of progression, and lead to improved management that would avoid the current over‐ and under‐treatment of patients.

Triple-negative breast cancers (TNBC), characterized by absence of estrogen receptor (ER), progesterone receptor (PR) and lack of overexpression of human epidermal growth factor receptor 2 (HER2), are typically associated with poor prognosis, due to aggressive tumor phenotype(s), only partial response to chemotherapy and present lack of clinically established targeted therapies. Advances in the design of individualized strategies for treatment of TNBC patients require further elucidation, by combined ‘omics’ approaches, of the molecular mechanisms underlying TNBC phenotypic heterogeneity, and the still poorly understood association of TNBC with BRCA1 mutations. An overview is here presented on TNBC profiling in terms of expression signatures, within the functional genomic breast tumor classification, and ongoing efforts toward identification of new therapy targets and bioimaging markers. Due to the complexity of aberrant molecular patterns involved in expression, pathological progression and biological/clinical heterogeneity, the search for novel TNBC biomarkers and therapy targets requires collection of multi-dimensional data sets, use of robust multivariate data analysis techniques and development of innovative systems biology approaches.