Digital histologic analysis reveals morphometric patterns of age-related involution in breast epithelium and stroma

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.

Toward a digital analysis of environmental impacts on rodent mammary gland density during critical developmental windows

While mammographic breast density is associated with breast cancer risk in humans, there is no comparable surrogate risk measure in mouse and rat mammary glands following various environmental exposures. In the current study, mammary glands from mice and rats subjected to reproductive factors and exposures to environmental chemicals that have been shown to influence mammary gland development and/or susceptibility to mammary tumors were evaluated for histologic density by manual and automated digital methods. Digital histological density detected changes due to hormonal stimuli/reproductive factors (parity), dietary fat, and exposure to environmental chemicals, such as benzophenone-3 and a combination of perfluorooctanoic acid and zeranol. Thus, digital analysis of mammary gland density offers a high throughput method that can provide a highly reproducible means of comparing a measure of histological density across independent experiments, experimental systems, and laboratories. This methodology holds promise for the detection of environmental impacts on mammary gland structure in mice and rats that may be comparable to human breast density, thus potentially allowing comparisons between rodent models and human breast cancer studies.

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.

Molecular and Clinical Characterization of Postpartum-Associated Breast Cancer in the Carolina Breast Cancer Study Phase I-III, 1993-2013

BACKGROUND

Breast cancers in recently postpartum women may have worse outcomes, but studies examining tumor molecular features by pregnancy recency have shown conflicting results.

METHODS

This analysis used Carolina Breast Cancer Study data to examine clinical and molecular tumor features among women less than 50 years of age who were recently (≤10 years prior) or remotely (>10 years prior) postpartum, or nulliparous. Prevalence odds ratios (POR) and 95% confidence intervals (CI) were estimated using multivariable models.

RESULTS

Recently postpartum women (N = 618) were more frequently lymph node-positive [POR (95% CI): 1.66 (1.26-2.19)], estrogen receptor (ER)-negative [1.37 (1.02-1.83)], and IHC-based triple negative [1.57 (1.00-2.47)] compared with nulliparous (N = 360) women. Some differences were identified between recent versus remotely postpartum; smaller tumor size [0.67 (0.52-0.86)], p53 wildtype [0.53 (0.36-0.77)], and non-basal-like phenotype [0.53 (0.33-0.84)] were more common among recently postpartum. Recently postpartum (vs. nulliparous) had significant enrichment for adaptive immunity, T cells, B cells, CD8 T cells, activated CD8 T cells/natural killer (NK) cells, and T follicular helper (Tfh) cells and higher overall immune cell scores. These differences were attenuated in remotely (compared with recently) postpartum women.

CONCLUSIONS

These results suggest a dominant effect of parity (vs. nulliparity) and a lesser effect of pregnancy recency on tumor molecular features, although tumor immune microenvironments were altered in association with pregnancy recency.

IMPACT

Our study is unique in examining tumor immune microenvironment and RNA-based markers according to time since last childbirth. Future studies should examine the interplay between tumor features, postdiagnostic treatment, and outcomes among recently postpartum women. See related commentary by McDonald et al., p. 518.

Immune checkpoint blockade reprograms systemic immune landscape and tumor microenvironment in obesity-associated breast cancer

Immune checkpoint blockade (ICB) has improved outcomes in some cancers. A major limitation of ICB is that most patients fail to respond, which is partly attributable to immunosuppression. Obesity appears to improve immune checkpoint therapies in some cancers, but impacts on breast cancer (BC) remain unknown. In lean and obese mice, tumor progression and immune reprogramming were quantified in BC tumors treated with anti-programmed death-1 (PD-1) or control. Obesity augments tumor incidence and progression. Anti-PD-1 induces regression in lean mice and potently abrogates progression in obese mice. BC primes systemic immunity to be highly responsive to obesity, leading to greater immunosuppression, which may explain greater anti-PD-1 efficacy. Anti-PD-1 significantly reinvigorates antitumor immunity despite persistent obesity. Laminin subunit beta-2 (Lamb2), downregulated by anti-PD-1, significantly predicts patient survival. Lastly, a microbial signature associated with anti-PD-1 efficacy is identified. Thus, anti-PD-1 is highly efficacious in obese mice by reinvigorating durable antitumor immunity. VIDEO ABSTRACT.

DeCompress: tissue compartment deconvolution of targeted mRNA expression panels using compressed sensing

Targeted mRNA expression panels, measuring up to 800 genes, are used in academic and clinical settings due to low cost and high sensitivity for archived samples. Most samples assayed on targeted panels originate from bulk tissue comprised of many cell types, and cell-type heterogeneity confounds biological signals. Reference-free methods are used when cell-type-specific expression references are unavailable, but limited feature spaces render implementation challenging in targeted panels. Here, we present DeCompress, a semi-reference-free deconvolution method for targeted panels. DeCompress leverages a reference RNA-seq or microarray dataset from similar tissue to expand the feature space of targeted panels using compressed sensing. Ensemble reference-free deconvolution is performed on this artificially expanded dataset to estimate cell-type proportions and gene signatures. In simulated mixtures, four public cell line mixtures, and a targeted panel (1199 samples; 406 genes) from the Carolina Breast Cancer Study, DeCompress recapitulates cell-type proportions with less error than reference-free methods and finds biologically relevant compartments. We integrate compartment estimates into cis-eQTL mapping in breast cancer, identifying a tumor-specific cis-eQTL for CCR3 (C-C Motif Chemokine Receptor 3) at a risk locus. DeCompress improves upon reference-free methods without requiring expression profiles from pure cell populations, with applications in genomic analyses and clinical settings.

Local Transdermal Delivery of Telapristone Acetate Through Breast Skin, Compared With Oral Treatment: A Randomized Double-Blind, Placebo-Controlled Phase II Trial

Oral breast cancer prevention medications entail systemic exposure, limiting acceptance by high-risk women. Delivery through the breast skin, although an attractive alternative, requires demonstration of drug distribution throughout the breast. We conducted a randomized double-blind, placebo-controlled phase II clinical trial comparing telapristone acetate, a progesterone receptor antagonist, administered orally (12 mg/day) or transdermally (12 mg/breast) for 4 ± 1 weeks to women planning mastectomy. Plasma and tissue concentrations, measured at five locations in the mastectomy specimen using liquid chromatography tandem mass spectrometry were compared. In 60 evaluable subjects, median drug concentration (ng/g tissue) was 103 (interquartile range (IQR): 46.3-336) in the oral vs. 2.82 (IQR: 1.4-5.5) in the transdermal group. Despite poor dermal permeation, within-breast drug distribution pattern was identical in both groups (R2  = 0.88, P = 0.006), demonstrating that transdermally and orally delivered drug is distributed similarly through the breast, and is strongly influenced by tissue adiposity (P < 0.0001). Other skin-penetrant drugs should be tested for breast cancer prevention.

Using Digital Pathology to Understand Epithelial Characteristics of Benign Breast Disease among Women Undergoing Diagnostic Image-Guided Breast Biopsy

Delayed terminal duct lobular unit (TDLU) involution is associated with elevated mammographic breast density (MD). Both are independent breast cancer risk factors among women with benign breast disease (BBD). Prior digital analyses of normal breast tissues revealed that epithelial nuclear density (END) and TDLU involution are inversely correlated. Accordingly, we examined associations of END, TDLU involution, and MD in BBD clinical biopsies. This study included digitized images of 262 representative image-guided hematoxylin and eosin-stained biopsies from 224 women diagnosed with BBD, enrolled within the cross-sectional BREAST-Stamp project that were visually assessed for TDLU involution (TDLU count/100 mm², median TDLU span and median acini count per TDLU). A digital algorithm estimated nuclei count per unit epithelial area, or END. Single X-ray absorptiometry of prebiopsy ipsilateral craniocaudal digital mammograms measured global and localized MD surrounding the biopsy region. Adjusted ordinal logistic regression models assessed relationships between tertiles of TDLU and END measures. Analysis of covariance examined mean differences in MD across END tertiles. TDLU measures were positively associated with increasing END tertiles [TDLU count/100 mm², OR_T3vsT1: 3.42, 95% confidence interval (CI), 1.87-6.28; acini count/TDLU_T3vsT1, OR: 2.40, 95% CI, 1.39-4.15]. END was significantly associated with localized, but not, global MD. Relationships were most apparent among patients with nonproliferative BBD. These findings suggest that quantitative END reflects different but complementary information to the histologic information captured by visual TDLU and radiologic MD measures and merits continued evaluation in assessing cellularity of breast parenchyma to understand the etiology of BBD.

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.

Hormone therapy use and breast tissue DNA methylation: analysis of epigenome wide data from the normal breast study

Hormone therapy (HT) is associated with increased risk of breast cancer, strongly dependent on type, duration, and recency of use. HT use could affect cancer risk by changing breast tissue transcriptional programs. We hypothesize that these changes are preceded by changes in DNA methylation. To explore this hypothesis we used histologically normal-appearing breast tissue from the Normal Breast Study (NBS). DNA methylation β-values were obtained using the Illumina HumanMethylation 450 BeadChips for 90 samples including all NBS-participants who used HT within 5 y before surgery. Data were analyzed using the reference-free cell mixture method. Cancer Genome Atlas (TCGA) mRNA-Seq data were used to assess correlation between DNA methylation and gene expression. We identified 527 CpG sites in 403 genes that were associated with ever using HT at genome wide significance (FDR q < 0.05), of these, 68 sites were also significantly associated with duration of use or recency of use. Twelve sites reached significance in all analyses one of which was cg01382688 in ARHGEF4 (p < 1.2x10⁻⁷). Mutations in ARHGEF4 have been reported in breast tumors, but this is the first report of possible breast cancer-related DNA methylation changes. In addition, 22 genes included more than one significant CpG site and a majority of these sites were significantly correlated with gene expression. Although based on small numbers, these findings support the hypothesis that HT is associated with epigenetic alterations in breast tissue, and identifies genes with altered DNA methylation states which could be linked to breast cancer development.

Stroma modifies relationships between risk factor exposure and age-related epithelial involution in benign breast

Delayed age-related lobular involution has been previously associated with elevated breast cancer risk. However, intraindividual variability in epithelial involution status within a woman is undefined. We developed a novel measure of age-related epithelial involution, density of epithelial nuclei in epithelial areas using digital image analysis in combination with stromal characteristics (percentage of section area comprising stroma). Approximately 1800 hematoxylin and eosin stained sections of benign breast tissue were evaluated from 416 participants having breast surgery for cancer or benign conditions. Two to sixteen slides per woman from different regions of the breast were studied. Epithelial involution status varied within a woman and as a function of stromal area. Percentage stromal area varied between samples from the same woman (median difference between highest and lowest stromal area within a woman was 7.5%, but ranged from 0.01 to 86.7%). Restricting to women with at least 10% stromal area (N = 317), epithelial nuclear density decreased with age (-637.1 cells/mm² per decade of life after age 40, p < 0.0001), increased with mammographic density (457.8 cells/mm² per increasing BI-RADs density category p = 0.002), and increased non-significantly with recent parity, later age at first pregnancy, and longer and more recent oral contraceptive use. These associations were attenuated in women with mostly fat samples (<10% stroma (N = 99)). Thirty-one percent of women evaluated had both adequate stroma (≥10%) and mostly fat (<10% stroma) regions of breast tissue, with the probability of having both types increasing with the number breast tissue samplings. Several breast cancer risk factors are associated with elevated age-related epithelial content, but associations depend upon stromal context. Stromal characteristics appear to modify relationships between risk factor exposures and breast epithelial involution.

Intra-individual Gene Expression Variability of Histologically Normal Breast Tissue

Several studies have sought to identify novel transcriptional biomarkers in normal breast or breast microenvironment to predict tumor risk and prognosis. However, systematic efforts to evaluate intra-individual variability of gene expression within normal breast have not been reported. This study analyzed the microarray gene expression data of 288 samples from 170 women in the Normal Breast Study (NBS), wherein multiple histologically normal breast samples were collected from different block regions and different sections at a given region. Intra-individual differences in global gene expression and selected gene expression signatures were quantified and evaluated in association with other patient-level factors. We found that intra-individual reliability was relatively high in global gene expression, but differed by signatures, with composition-related signatures (i.e., stroma) having higher intra-individual variability and tumorigenesis-related signatures (i.e., proliferation) having lower intra-individual variability. Histological stroma composition was the only factor significantly associated with heterogeneous breast tissue (defined as > median intra-individual variation; high nuclear density, odds ratio [OR] = 3.42, 95% confidence interval [CI] = 1.15–10.15; low area, OR = 0.29, 95% CI = 0.10–0.86). Other factors suggestively influencing the variability included age, BMI, and adipose nuclear density. Our results underscore the importance of considering intra-individual variability in tissue-based biomarker development, and have important implications for normal breast research.

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.

Mammary Gland Evaluation in Juvenile Toxicity Studies: Temporal Developmental Patterns in the Male and Female Harlan Sprague-Dawley Rat

There are currently no reports describing mammary gland development in the Harlan Sprague-Dawley (HSD) rat, the current strain of choice for National Toxicology Program (NTP) testing. Our goals were to empower the NTP, contract labs, and other researchers in understanding and interpreting chemical effects in this rat strain. To delineate similarities/differences between the female and male mammary gland, data were compiled starting on embryonic day 15.5 through postnatal day 70. Mammary gland whole mounts, histology sections, and immunohistochemically stained tissues for estrogen, progesterone, and androgen receptors were evaluated in both sexes; qualitative and quantitative differences are highlighted using a comprehensive visual timeline. Research on endocrine disrupting chemicals in animal models has highlighted chemically induced mammary gland anomalies that may potentially impact human health. In order to investigate these effects within the HSD strain, 2,3,7,8-tetrachlorodibenzo-p-dioxin, diethylstilbestrol, or vehicle control was gavage dosed on gestation day 15 and 18 to demonstrate delayed, accelerated, and control mammary gland growth in offspring, respectively. We provide illustrations of normal and chemically altered mammary gland development in HSD male and female rats to help inform researchers unfamiliar with the tissue and may facilitate enhanced evaluation of both male and female mammary glands in juvenile toxicity studies.

Standardized measures of lobular involution and subsequent breast cancer risk among women with benign breast disease: a nested case-control study

Lesser degrees of terminal duct-lobular unit (TDLU) involution predict higher breast cancer risk; however, standardized measures to quantitate levels of TDLU involution have only recently been developed. We assessed whether three standardized measures of TDLU involution, with high intra/inter pathologist reproducibility in normal breast tissue, predict subsequent breast cancer risk among women in the Mayo benign breast disease (BBD) cohort. We performed a masked evaluation of biopsies from 99 women with BBD who subsequently developed breast cancer (cases) after a median of 16.9 years and 145 age-matched controls. We assessed three metrics inversely related to TDLU involution: TDLU count/mm(2), median TDLU span (microns, which approximates acini content), and median category of acini counts/TDLU (0-10; 11-20; 21-30; 31-50; >50). Associations with subsequent breast cancer risk for quartiles (or categories of acini counts) of each of these measures were assessed with multivariable conditional logistic regression to estimate odds ratios (ORs) and 95 % confidence intervals (CI). In multivariable models, women in the highest quartile compared to the lowest quartiles of TDLU counts and TDLU span measures were significantly associated with subsequent breast cancer diagnoses; TDLU counts quartile4 versus quartile1, OR = 2.44, 95 %CI 0.96-6.19, p-trend = 0.02; and TDLU spans, quartile4 versus quartile1, OR = 2.83, 95 %CI = 1.13-7.06, p-trend = 0.03. Significant associations with categorical measures of acini counts/TDLU were also observed: compared to women with median category of <10 acini/TDLU, women with >25 acini counts/TDLU were at significantly higher risk, OR = 3.40, 95 %CI 1.03-11.17, p-trend = 0.032. Women with TDLU spans and TDLU count measures above the median were at further increased risk, OR = 3.75 (95 %CI 1.40-10.00, p-trend = 0.008), compared with women below the median for both of these metrics. Similar results were observed for combinatorial metrics of TDLU acini counts/TDLU, and TDLU count. Standardized quantitative measures of TDLU counts and acini counts approximated by TDLU span measures or visually assessed in categories are independently associated with breast cancer risk. Visual assessment of TDLU numbers and acini content, which are highly reproducible between pathologists, could help identify women at high risk for subsequent breast cancer among the million women diagnosed annually with BBD in the US.