Search

PURPOSE: Patients with postpartum breast cancer diagnosed after cessation of breastfeeding (postweaning, PP-BCPW) have a particularly poor prognosis compared with patients diagnosed during lactation (PP-BCDL), or to pregnant (Pr-BC) and nulliparous (NP-BC) patients, regardless of standard prognostic characteristics. Animal studies point to a role of the involution process in stimulation of tumor growth in the mammary gland. However, in women, the molecular mechanisms that underlie this poor prognosis of patients with PP-BCPW remain vastly underexplored, due to of lack of adequate patient numbers and outcome data. EXPERIMENTAL DESIGN: We explored whether distinct prognostic features, common to all breast cancer molecular subtypes, exist in postpartum tumor tissue. Using detailed breastfeeding data, we delineated the postweaning period in PP-BC as a surrogate for mammary gland involution and performed whole transcriptome sequencing, immunohistochemical, and (multiplex) immunofluorescent analyses on tumor tissue of patients with PP-BCPW, PP-BCDL, Pr-BC, and NP-BC. RESULTS: We found that patients with PP-BCPW having a low expression level of an immunoglobulin gene signature, but high infiltration of plasma B cells, have an increased risk for metastasis and death. Although PP-BCPW tumor tissue was also characterized by an increase in CD8+ cytotoxic T cells and reduced distance among these cell types, these parameters were not associated with differential clinical outcomes among groups. CONCLUSIONS: These data point to the importance of plasma B cells in the postweaning mammary tumor microenvironment regarding the poor prognosis of PP-BCPW patients. Future prospective and in-depth research needs to further explore the role of B-cell immunobiology in this specific group of young patients with breast cancer., Pattern Recognition and Bioinformatics

Purpose: To assess whether contralateral parenchymal enhancement (CPE) on MRI is associated with gene expression pathways in ER+/HER2-breast cancer, and if so, whether such pathways are related to survival. Methods: Preoperative breast MRIs were analyzed of early ER+/HER2-breast cancer patients eligible for breast-conserving surgery included in a prospective observational cohort study (MARGINS). The contralateral parenchyma was segmented and CPE was calculated as the average of the top-10% delayed enhancement. Total tumor RNA sequencing was performed and gene set enrichment analysis was used to reveal gene expression pathways associated with CPE (N = 226) and related to overall survival (OS) and invasive disease-free survival (IDFS) in multivariable survival analysis. The latter was also done for the METABRIC cohort (N = 1355). Results: CPE was most strongly correlated with proteasome pathways (normalized enrichment statistic = 2.04, false discovery rate =.11). Patients with high CPE showed lower tumor proteasome gene expression. Proteasome gene expression had a hazard ratio (HR) of 1.40 (95% CI = 0.89, 2.16; P =.143) for OS in the MARGINS cohort and 1.53 (95% CI = 1.08, 2.14; P =.017) for IDFS, in METABRIC proteasome gene expression had an HR of 1.09 (95% CI = 1.01, 1.18; P =.020) for OS and 1.10 (95% CI = 1.02, 1.18; P =.012) for IDFS. Conclusion: CPE was negatively correlated with tumor proteasome gene expression in early ER+/HER2-breast cancer patients. Low tumor proteasome gene expression was associated with improved survival in the METABRIC data., Pattern Recognition and Bioinformatics

Purpose: To assess whether contralateral parenchymal enhancement (CPE) on MRI is associated with gene expression pathways in ER+/HER2-breast cancer, and if so, whether such pathways are related to survival. Methods: Preoperative breast MRIs were analyzed of early ER+/HER2-breast cancer patients eligible for breast-conserving surgery included in a prospective observational cohort study (MARGINS). The contralateral parenchyma was segmented and CPE was calculated as the average of the top-10% delayed enhancement. Total tumor RNA sequencing was performed and gene set enrichment analysis was used to reveal gene expression pathways associated with CPE (N = 226) and related to overall survival (OS) and invasive disease-free survival (IDFS) in multivariable survival analysis. The latter was also done for the METABRIC cohort (N = 1355). Results: CPE was most strongly correlated with proteasome pathways (normalized enrichment statistic = 2.04, false discovery rate =.11). Patients with high CPE showed lower tumor proteasome gene expression. Proteasome gene expression had a hazard ratio (HR) of 1.40 (95% CI = 0.89, 2.16; P =.143) for OS in the MARGINS cohort and 1.53 (95% CI = 1.08, 2.14; P =.017) for IDFS, in METABRIC proteasome gene expression had an HR of 1.09 (95% CI = 1.01, 1.18; P =.020) for OS and 1.10 (95% CI = 1.02, 1.18; P =.012) for IDFS. Conclusion: CPE was negatively correlated with tumor proteasome gene expression in early ER+/HER2-breast cancer patients. Low tumor proteasome gene expression was associated with improved survival in the METABRIC data.

Background: Better understanding of the molecular biology associated with MRI phenotypes may aid in the diagnosis and treatment of breast cancer. Purpose: To discover the associations between MRI phenotypes of breast cancer and their underlying molecular biology derived from gene expression data. Materials and Methods: This is a secondary analysis of the Multimodality Analysis and Radiologic Guidance in Breast-Conserving Therapy, or MARGINS, study. MARGINS included patients eligible for breast-conserving therapy between November 2000 and December 2008 for preoperative breast MRI. Tumor RNA was collected for sequencing from surgical specimen. Twenty-one computer-generated MRI features of tumors were condensed into seven MRI factors related to tumor size, shape, initial enhancement, late enhancement, smoothness of enhancement, sharpness, and sharpness variation. These factors were associated with gene expression levels from RNA sequencing by using gene set enrichment analysis. Statistical significance of these associations was evaluated by using a sample permutation test and the false discovery rate. Results: Gene expression and MRI data were obtained for 295 patients (mean age, 56 years 6 10.3 [standard deviation]). Larger and more irregular tumors showed increased expression of cell cycle and DNA damage checkpoint genes (false discovery rate ,0.25; normalized enrichment statistic [NES], 2.15). Enhancement and sharpness of the tumor margin were associated with expression of ribosomal proteins (false discovery rate ,0.25; NES, 1.95). Smoothness of enhancement, tumor size, and tumor shape were associated with expression of genes involved in the extracellular matrix (false discovery rate ,0.25; NES, 2.25). Conclusion: Breast cancer MRI phenotypes were related to their underlying molecular biology revealed by using RNA sequencing. The association between enhancements and sharpness of the tumor margin with the ribosome suggests that these MRI features may be imaging bioma

Integrative analyses that summarize and link molecular data to treatment sensitivity are crucial to capture the biological complexity which is essential to further precision medicine. We introduce Weighted Orthogonal Nonnegative parallel factor analysis (WON-PARAFAC), a data integration method that identifies sparse and interpretable factors. WON-PARAFAC summarizes the GDSC1000 cell line compendium in 130 factors. We interpret the factors based on their association with recurrent molecular alterations, pathway enrichment, cancer type, and drug-response. Crucially, the cell line derived factors capture the majority of the relevant biological variation in Patient-Derived Xenograft (PDX) models, strongly suggesting our factors capture invariant and generalizable aspects of cancer biology. Furthermore, drug response in cell lines is better and more consistently translated to PDXs using factor-based predictors as compared to raw feature-based predictors. WON-PARAFAC efficiently summarizes and integrates multiway high-dimensional genomic data and enhances translatability of drug response prediction from cell lines to patient-derived xenografts., Pattern Recognition and Bioinformatics

Purpose: To retrospectively explore the relation between parenchymal enhancement of the healthy contralateral breast on dynamic contrast-enhanced magnetic resonance imaging (MRI) and genomic tests for estrogen receptor (ER)-pathway activity in patients with ER-positive/HER2-negative cancer. Methods: A subset of 227 consecutively included patients with unilateral invasive ER-positive/HER2-negative breast cancer underwent dynamic contrast-enhanced MRI prior to breast-conserving therapy between 2000 and 2008. Perfusion of the parenchyma in the healthy breast was assessed using a previously reported measure of contralateral parenchymal enhancement (CPE), consisting of the mean of the top-10% late enhancement. ER-pathway activity was assessed from the surgical resection specimen by the previously reported sensitivity to endocrine therapy (SET)-index and ER-factor. The SET-index is a genetic test to estimate survival benefit from endocrine therapy, consisting of genes related to the ESR1 gene. The ER-factor examines other factors as well including protein expression. The relation between CPE and ER-pathway activity was modeled using linear regression. Results: Patients had a median age of 59 years. CPE was not significantly associated with the SET-index (R-squared = 0.005) nor the ER-factor (R-squared = 0.0002). The only variable significantly different between low and high CPE was age at diagnosis (P < 0.001). Conclusions: Contralateral parenchymal enhancement on dynamic contrast-enhanced MRI was not associated with tumor-derived estrogen receptor pathway activity.

Effective cancer treatment is crucially dependent on the identification of the biological processes that drive a tumor. However, multiple processes may be active simultaneously in a tumor. Clustering is inherently unsuitable to this task as it assigns a tumor to a single cluster. In addition, the wide availability of multiple data types per tumor provides the opportunity to profile the processes driving a tumor more comprehensively. Here we introduce Functional Sparse-Factor Analysis (funcSFA) to address these challenges. FuncSFA integrates multiple data types to define a lower dimensional space capturing the relevant variation. A tailor-made module associates biological processes with these factors. FuncSFA is inspired by iCluster, which we improve in several key aspects. First, we increase the convergence efficiency significantly, allowing the analysis of multiple molecular datasets that have not been pre-matched to contain only concordant features. Second, FuncSFA does not assign tumors to discrete clusters, but identifies the dominant driver processes active in each tumor. This is achieved by a regression of the factors on the RNA expression data followed by a functional enrichment analysis and manual curation step. We apply FuncSFA to the TCGA breast and lung datasets. We identify EMT and Immune processes common to both cancer types. In the breast cancer dataset we recover the known intrinsic subtypes and identify additional processes. These include immune infiltration and EMT, and processes driven by copy number gains on the 8q chromosome arm. In lung cancer we recover the major types (adenocarcinoma and squamous cell carcinoma) and processes active in both of these types. These include EMT, two immune processes, and the activity of the NFE2L2 transcription factor. We validate the breast cancer findings on the METABRIC set and demonstrate the translatability of the TCGA breast cancer factors to METABRIC. In summary, FuncSFA is a robust method to perform dis, Pattern Recognition and Bioinformatics

Background: Patients with BRCA1-like tumors correlate with improved response to DNA double-strand break-inducing therapy. A gene expression-based classifier was developed to distinguish between BRCA1-like and non-BRCA1-like tumors. We hypothesized that these tumors may also be more sensitive to PARP inhibitors than standard treatments. Methods: A diagnostic gene expression signature (BRCA1ness) was developed using a centroid model with 128 triple-negative breast cancer samples from the EU FP7 RATHER project. This BRCA1 ness signature was then tested in HER2-negative patients (n= 116) from the I-SPY 2 TRIAL who received an oral PARP inhibitor veliparib in combination with carboplatin (V-C), or standard chemotherapy alone. We assessed the association between BRCA1 ness and pathologic complete response in the V-C and control arms alone using Fisher ’s exact test, and the relative performance between arms (biomarker × treatment interaction, likelihood ratio p< 0.05) using a logistic model and adjusting for hormone receptor status (HR). Results: We developed a gene expression signature to identify BRCA1-like status. In the I-SPY 2 neoadjuvant setting the BRCA1ness signature associated signif icantly with response to V-C (p= 0.03), but not in the control arm (p = 0.45). We identified a significant interaction between BRCA1ness and V-C (p= 0.023) after correcting for HR. Conclusions: A genomic-based BRCA1-like signature was successfully translated to an expression-based signature (BRC1Aness). In the I-SPY 2 neoadjuvant setting, we determined that the BRCA1ness signature is capable of predicting benefit of V-C added to standard chemotherapy compared to standard chemotherapy alone. Trial registration: I-SPY 2 TRIAL beginning December 31, 2009: Neoadjuvant and Personalized Adaptive Novel Agents to Treat Breast Cancer (I-SPY 2), NCT01042379., Pattern Recognition and Bioinformatics

Purpose: To determine whether markers of healthy breast stroma are able to select a subgroup of patients at low risk of death or metastasis from patients considered at high risk according to routine markers of the tumor. Experimental Design: Patients with ER+/HER2− breast cancer were consecutively included for retrospective analysis. The contralateral parenchyma was segmented automatically on dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), where upon the average of the top-10% late enhancement was calculated. This contralateral parenchymal enhancement (CPE) was analyzed with respect to routine prognostic models and molecular assays (Nottingham Prognostic Index, Dutch clinical chemotherapy-selection guidelines, 70-gene signature, and 21-gene recurrence score). CPE was split in tertiles and tested for overall and distant disease-free survival. CPE was adjusted for patient and tumor characteristics, as well as systemic therapy, using inverse probability weighting (IPW). Subanalyses were performed in patients at high risk according to prognostic models and molecular assays. Results: Four-hundred-and-fifteen patients were included, constituting the same group in which the association between CPE and survival was discovered. Median follow-up was 85 months, 34/415(8%) patients succumbed. After IPW-adjustment for patient and tumor characteristics, patients with high CPE had significantly better overall survival than those with low CPE in groups at high risk according to the Nottingham Prognostic Index [HR (95% CI): 0.08 (0.00–0.40), P < 0.001]; Dutch clinical guidelines [HR (95% CI): 0.22 (0.00–0.81), P = 0.021]; and 21-gene recurrence score [HR (95% CI): 0.14 (0.00–0.84), P = 0.030]. One group showed a trend [70-gene signature: HR (95% CI): 0.25 (0.00–1.02), P = 0.054]. Conclusions: In patients at high risk based on the tumor, subgroups at relatively low risk were identified using pretreatment enhancement of the stroma on breast DCE-MRI.