Quantitative analysis of breast cancer tissue composition and associations with tumor subtype

The tumor microenvironment is an important determinant of breast cancer progression, but standard methods for describing the tumor microenvironment are lacking. Measures of microenvironment composition such as stromal area and immune infiltrate are labor-intensive and few large studies have systematically collected this data. However, digital histologic approaches are becoming more widely available, allowing high-throughput, quantitative estimation. We applied such methods to tissue microarrays of tumors from 1687 women (mean 4 cores per case) in the Carolina Breast Cancer Study Phase 3. Tumor composition was quantified as percentage of epithelium, stroma, adipose, and lymphocytic infiltrate (with the latter as presence/absence using a ≥1% cutoff). Composition proportions and presence/absence were evaluated in association with clinical and molecular features of breast cancer (intrinsic subtype and RNA-based risk of recurrence [ROR] scores) using multivariable linear and logistic regression. Lower stromal content was associated with aggressive tumor phenotypes, including triple-negative (31.1% vs. 41.6% in HR+/HER2-; RFD [95% CI]: −10.5%, [-13.1, −7.9]), Basal-like subtypes (29.0% vs. 44.0% in Luminal A; RFD [95% CI]: −14.9%, [-17.8, −12.0]), and high RNA-based PAM50 ROR scores (27.6% vs. 48.1% in ROR low; RFD [95% CI]: −20.5%, [24.3, 16.7]), after adjusting for age and race. HER2+ tumors also had lower stromal content, particularly among RNA-based HER2-enriched (35.2% vs. 44.0% in Luminal A; RFD [95% CI]: −8.8%, [-13.8, −3.8]). Similar associations were observed between immune infiltrate and tumor phenotypes. Quantitative digital image analysis of the breast cancer microenvironment showed significant associations with demographic characteristics and biological indicators of aggressive behavior

Vascular density of histologically benign breast tissue from women with breast cancer: associations with tissue composition and tumor characteristics

In breast tumors, it is well established that intratumoral angiogenesis is crucial for malignant progression, but little is known about the vascular characteristics of extratumoral, cancer-adjacent breast. Genome-wide transcriptional data suggest that extratumoral microenvironments may influence breast cancer phenotypes; thus, histologic features of cancer-adjacent tissue may also have clinical implications. To this end, we developed a digital algorithm to quantitate vascular density in approximately 300 histologically benign tissue specimens from breast cancer patients enrolled in the UNC Normal Breast Study (NBS). Specimens were stained for CD31, and vascular content was compared to demographic variables, tissue composition metrics, and tumor molecular features. We observed that the vascular density of cancer-adjacent breast was significantly higher in older and obese women, and was strongly associated with breast adipose tissue content. Consistent with observations that older and heavier women experience higher frequencies of ER+ disease, higher extratumoral vessel density was also significantly associated with positive prognostic tumor features such as lower stage, negative nodal status, and smaller size (<2 cm). These results reveal biological relationships between extratumoral vascular content and body size, breast tissue composition, and tumor characteristics, and suggest biological plausibility for the relationship between weight gain (and corresponding breast tissue changes) and breast cancer progression

Reproducibility and intratumoral heterogeneity of the PAM50 breast cancer assay

Background: The PAM50 assay is used routinely in clinical practice to determine breast cancer prognosis and management; however, research assessing how technical variation and intratumoral heterogeneity contribute to misclassification and reproducibility of these tests is limited. Methods: We evaluated the impact of intratumoral heterogeneity on the reproducibility of results for the PAM50 assay by testing RNA extracted from formalin-fixed paraffin embedded breast cancer blocks sampled at distinct spatial locations. Samples were classified according to intrinsic subtype (Luminal A, Luminal B, HER2-enriched, Basal-like, or Normal-like) and risk of recurrence with proliferation score (ROR-P, high, medium, or low). Intratumoral heterogeneity and technical reproducibility (replicate assays on the same RNA) were assessed as percent categorical agreement between paired intratumoral and replicate samples. Euclidean distances between samples, calculated across the PAM50 genes and the ROR-P score, were compared for concordant vs. discordant samples. Results: Technical replicates (N = 144) achieved 93% agreement for ROR-P group and 90% agreement on PAM50 subtype. For spatially distinct biological replicates (N = 40 intratumoral replicates), agreement was lower (81% for ROR-P and 76% for PAM50 subtype). The Euclidean distances between discordant technical replicates were bimodal, with discordant samples showing higher Euclidian distance and biologic heterogeneity. Conclusion: The PAM50 assay achieved very high technical reproducibility for breast cancer subtyping and ROR-P, but intratumoral heterogeneity is revealed by the assay in a small proportion of cases

Correlated metabolomic, genomic, and histologic phenotypes in histologically normal breast tissue

Breast carcinogenesis is a multistep process accompanied by widespread molecular and genomic alterations, both in tumor and in surrounding microenvironment. It is known that tumors have altered metabolism, but the metabolic changes in normal or cancer-adjacent, nonmalignant normal tissues and how these changes relate to alterations in gene expression and histological composition are not well understood. Normal or cancer-adjacent normal breast tissues from 99 women of the Normal Breast Study (NBS) were evaluated. Data of metabolomics, gene expression and histological composition was collected by mass spectrometry, whole genome microarray, and digital image, respectively. Unsupervised clustering analysis determined metabolomics-derived subtypes. Their association with genomic and histological features, as well as other breast cancer risk factors, genomic and histological features were evaluated using logistic regression. Unsupervised clustering of metabolites resulted in two main clusters. The metabolite differences between the two clusters suggested enrichment of pathways involved in lipid metabolism, cell growth and proliferation, and migration. Compared with Cluster 1, subjects in Cluster 2 were more likely to be obese (body mass index 30 kg/m2, p<0.05), have increased adipose proportion (p<0.01) and associated with a previously defined Active genomic subtype (p<0.01). By the integrated analyses of histological, metabolomics and transcriptional data, we characterized two distinct subtypes of non-malignant breast tissue. Further research is needed to validate our findings, and understand the potential role of these alternations in breast cancer initiation, progression and recurrence

Metabolic reprogramming of macrophages: Glucose transporter 1 (GLUT1)-mediated glucose metabolism drives a proinflammatory phenotype

Glucose is a critical component in the proinflammatory response of macrophages (MΦs). However, the contribution of glucose transporters (GLUTs) and the mechanisms regulating subsequent glucose metabolism in the inflammatory response are not well understood. Because MΦs contribute to obesity-induced inflammation, it is important to understand how substrate metabolism may alter inflammatory function. We report that GLUT1 (SLC2A1) is the primary rate-limiting glucose transporter on proinflammatory-polarized MΦs. Furthermore, in high fat diet-fed rodents, MΦs in crown-like structures and inflammatory loci in adipose and liver, respectively, stain positively for GLUT1. We hypothesized that metabolic reprogramming via increased glucose availability could modulate the MΦ inflammatory response. To increase glucose uptake, we stably overexpressed the GLUT1 transporter in RAW264.7 MΦs (GLUT1-OE MΦs). Cellular bioenergetics analysis, metabolomics, and radiotracer studies demonstrated that GLUT1 overexpression resulted in elevated glucose uptake and metabolism, increased pentose phosphate pathway intermediates, with a complimentary reduction in cellular oxygen consumption rates. Gene expression and proteome profiling analysis revealed that GLUT1-OE MΦs demonstrated a hyperinflammatory state characterized by elevated secretion of inflammatory mediators and that this effect could be blunted by pharmacologic inhibition of glycolysis. Finally, reactive oxygen species production and evidence of oxidative stress were significantly enhanced in GLUT1-OEMΦs; antioxidant treatment blunted the expression of inflammatory mediators such as PAI-1 (plasminogen activator inhibitor 1), suggesting that glucose-mediated oxidative stress was driving the proinflammatory response. Our results indicate that increased utilization of glucose induced a ROS-driven proinflammatory phenotype in MΦs, which may play an integral role in the promotion of obesity-associated insulin resistance

Hepatocyte growth factor pathway expression in breast cancer by race and subtype

Background: African American women have the highest risk of breast cancer mortality compared to other racial groups. Differences in tumor characteristics have been implicated as a possible cause; however, the tumor microenvironment may also contribute to this disparity in mortality. Hepatocyte growth factor (HGF) is a stroma-derived marker of the tumor microenvironment that may affect tumor progression differentially by race. Objective: To examine whether an HGF gene expression signature is differentially expressed by race and tumor characteristics. Methods: Invasive breast tumors from 1957 patients were assessed for a 38-gene RNA-based HGF gene expression signature. Participants were black (n = 1033) and non-black (n = 924) women from the population-based Carolina Breast Cancer Study (1993–2013). Generalized linear models were used to estimate the relative frequency differences (RFD) in HGF status by race, clinical, and demographic factors. Results: Thirty-two percent of tumors were positive for the HGF signature. Black women were more likely [42% vs. 21%; RFD = + 19.93% (95% CI 16.00, 23.87)] to have HGF-positive tumors compared to non-black women. Triple-negative patients had a higher frequency of HGF positivity [82% vs. 13% in non-triple-negative; RFD = + 65.85% (95% CI 61.71, 69.98)], and HGF positivity was a defining feature of basal-like subtype [92% vs. 8% in non-basal; RFD = + 81.84% (95% CI 78.84, 84.83)]. HGF positivity was associated with younger age, stage, higher grade, and high genomic risk of recurrence (ROR-PT) score. Conclusion: HGF expression is a defining feature of basal-like tumors, and its association with black race and young women suggests it may be a candidate pathway for understanding breast cancer disparities

The Landscape of Immune Microenvironments in Racially Diverse Breast Cancer Patients

Background: Immunotherapy is a rapidly evolving treatment option in breast cancer; However, the breast cancer immune microenvironment is understudied in Black and younger (<50 years) patients. Methods: We used histologic and RNA-based immunoprofiling methods to characterize the breast cancer immune landscape in 1,952 tumors from the Carolina Breast Cancer Study (CBCS), a population-based study that oversampled Black (n ¼ 1,030) and young women (n ¼ 1,039). We evaluated immune response leveraging markers for 10 immune cell populations, compared profiles to those in The Cancer Genome Atlas (TCGA) Project [n ¼ 1,095 tumors, Black (n ¼ 183), and young women (n ¼ 295)], and evaluated in association with clinical and demographic variables, including recurrence. Results: Consensus clustering identified three immune clusters in CBCS (adaptive-enriched, innate-enriched, or immune-quiet) that varied in frequency by race, age, tumor grade and subtype; however, only two clusters were identified in TCGA, which were predominantly comprised of adaptive-enriched and innate-enriched tumors. In CBCS, the strongest adaptive immune response was observed for basal-like, HER2-positive (HER2þ), triple-negative breast cancer (TNBC), and high-grade tumors. Younger patients had higher proportions of adaptive-enriched tumors, particularly among estrogen receptor (ER)-negative (ER-) cases. Black patients had higher frequencies of both adaptive-enriched and innate-enriched tumors. Immune clusters were associated with recurrence among ER- tumors, with adaptive-enriched showing the best and innate-enriched showing the poorest 5-year recurrence-free survival. Conclusions: These data suggest that immune microenvironments are intricately related to race, age, tumor subtype, and grade. Impact: Given higher mortality among Black and young women, more defined immune classification using cell-type–specific panels could help explain higher recurrence and ultimately lead to target-able interventions

Racial differences in breast cancer outcomes by hepatocyte growth factor pathway expression

Purpose: Black women have a 40% increased risk of breast cancer-related mortality. These outcome disparities may reflect differences in tumor pathways and a lack of targetable therapies for specific subtypes that are more common in Black women. Hepatocyte growth factor (HGF) is a targetable pathway that promotes breast cancer tumorigenesis, is associated with basal-like breast cancer, and is differentially expressed by race. This study assessed whether a 38-gene HGF expression signature is associated with recurrence and survival in Black and non-Black women. Methods: Study participants included 1957 invasive breast cancer cases from the Carolina Breast Cancer Study. The HGF signature was evaluated in association with recurrence (n = 1251, 171 recurrences), overall, and breast cancer-specific mortality (n = 706, 190/328 breast cancer/overall deaths) using Cox proportional hazard models. Results: Women with HGF-positive tumors had higher recurrence rates [HR 1.88, 95% CI (1.19, 2.98)], breast cancer-specific mortality [HR 1.90, 95% CI (1.26, 2.85)], and overall mortality [HR 1.69; 95% CI (1.17, 2.43)]. Among Black women, HGF positivity was significantly associated with higher 5-year rate of recurrence [HR 1.73; 95% CI (1.01, 2.99)], but this association was not significant in non-Black women [HR 1.68; 95% CI (0.72, 3.90)]. Among Black women, HGF-positive tumors had elevated breast cancer-specific mortality [HR 1.80, 95% CI (1.05, 3.09)], which was not significant in non-Black women [HR 1.52; 95% CI (0.78, 2.99)]. Conclusion: This multi-gene HGF signature is a poor-prognosis feature for breast cancer and may identify patients who could benefit from HGF-targeted treatments, an unmet need for Black and triple-negative patients

Bimodal age distribution at diagnosis in breast cancer persists across molecular and genomic classifications

Purpose: Female breast cancer demonstrates bimodal age frequency distribution patterns at diagnosis, interpretable as two main etiologic subtypes or groupings of tumors with shared risk factors. While RNA-based methods including PAM50 have identified well-established clinical subtypes, age distribution patterns at diagnosis as a proxy for etiologic subtype are not established for molecular and genomic tumor classifications. Methods: We evaluated smoothed age frequency distributions at diagnosis for Carolina Breast Cancer Study cases within immunohistochemistry-based and RNA-based expression categories. Akaike information criterion (AIC) values compared the fit of single density versus two-component mixture models. Two-component mixture models estimated the proportion of early-onset and late-onset categories by immunohistochemistry-based ER (n = 2860), and by RNA-based ESR1 and PAM50 subtype (n = 1965). PAM50 findings were validated using pooled publicly available data (n = 8103). Results: Breast cancers were best characterized by bimodal age distribution at diagnosis with incidence peaks near 45 and 65 years, regardless of molecular characteristics. However, proportional composition of early-onset and late-onset age distributions varied by molecular and genomic characteristics. Higher ER-protein and ESR1-RNA categories showed a greater proportion of late age-at-onset. Similarly, PAM50 subtypes showed a shifting age-at-onset distribution, with most pronounced early-onset and late-onset peaks found in Basal-like and Luminal A, respectively. Conclusions: Bimodal age distribution at diagnosis was detected in the Carolina Breast Cancer Study, similar to national cancer registry data. Our data support two fundamental age-defined etiologic breast cancer subtypes that persist across molecular and genomic characteristics. Better criteria to distinguish etiologic subtypes could improve understanding of breast cancer etiology and contribute to prevention efforts

Outcomes of Hormone-Receptor Positive, HER2-Negative Breast Cancers by Race and Tumor Biological Features

Background: Black women have higher hormone receptor positive (HR+) breast cancer mortality than White women. Early recurrence rates differ by race, but little is known about genomic predictors of early recurrence among HR+ women. Methods: Using data from the Carolina Breast Cancer Study (phase III, 2008-2013), we estimated associations between race and recurrence among nonmetastatic HR+/HER2-negative tumors, overall and by PAM50 Risk of Recurrence score, PAM50 intrinsic subtype, and tumor grade using survival curves and Cox models standardized for age and stage. Relative frequency differences (RFD) were estimated using multivariable linear regression. To assess intervention opportunities, we evaluated treatment patterns by race among patients with high-risk disease. Results: Black women had higher recurrence risk relative to White women (crude hazard ratio = 1.81, 95% confidence interval [CI] = 1.34 to 2.46), which remained elevated after standardizing for clinical covariates (hazard ratio = 1.42, 95% CI = 1.05 to 1.93). Racial disparities were most pronounced among those with high PAM50 Risk of Recurrence score (5-year standardized recurrence risk = 18.9%, 95% CI = 8.6% to 29.1% in Black women vs 12.5%, 95% CI = 2.0% to 23.0% in White women) and high grade (5-year standardized recurrence risk = 16.6%, 95% CI = 11.7% to 21.5% in Black women vs 12.0%, 95% CI = 7.3% to 16.7% in White women). However, Black women with high-grade tumors were statistically significantly less likely to initiate endocrine therapy (RFD = -8.3%, 95% CI = -15.9% to -0.6%) and experienced treatment delay more often than White women (RFD = +9.0%, 95% CI = 0.3% to 17.8%). Conclusions: Differences in recurrence by race appear greatest among women with aggressive tumors and may be influenced by treatment differences. Efforts to identify causes of variation in cancer treatment are critical to reducing outcome disparities