Abstract
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.