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.
Androgen receptor–positive prostate cancer (PCa) and estrogen receptor–positive luminal breast cancer (BCa) are generally less responsive to immunotherapy compared with certain tumor types such as melanoma. However, the underlying mechanisms are not fully elucidated. In this study, we found that FOXA1 overexpression inversely correlated with interferon (IFN) signature and antigen presentation gene expression in PCa and BCa patients. FOXA1 bound the STAT2 DNA-binding domain and suppressed STAT2 DNA-binding activity, IFN signaling gene expression, and cancer immune response independently of the transactivation activity of FOXA1 and its mutations detected in PCa and BCa. Increased FOXA1 expression promoted cancer immuno- and chemotherapy resistance in mice and PCa and BCa patients. These findings were also validated in bladder cancer expressing high levels of FOXA1. FOXA1 overexpression could be a prognostic factor to predict therapy resistance and a viable target to sensitize luminal PCa, BCa, and bladder cancer to immuno- and chemotherapy.
Postmenopausal obesity is associated with increased circulating levels of androgens and estrogens and elevated breast cancer risk. Crown-like structures (CLS; microscopic foci of dying adipocytes surrounded by macrophages) are proposed to represent sites of increased aromatization of androgens to estrogens. Accordingly, we examined relationships between CLS and sex-steroid hormones in breast adipose tissue and serum from postmenopausal breast cancer patients. Formalin-fixed paraffin embedded benign breast tissues collected for research from postmenopausal women (n = 83) diagnosed with invasive breast cancer in the Polish Breast Cancer Study (PBCS) were evaluated. Tissues were immunohistochemically stained for CD68 to determine the presence of CLS per unit area of adipose tissue. Relationships were assessed between CD68 density and CLS and previously reported sex-steroid hormones quantified using radioimmunoassays in serum taken at the time of diagnosis and in fresh frozen adipose tissue taken at the time of surgery. Logistic regression analysis was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs) for the relationships between hormones (in tertiles) and CLS. CLS were observed in 36% of benign breast tissues, with a higher frequency among obese versus lean women (54% versus 17%, p = 0.03). Detection of CLS was not related to individual hormone levels or breast tumor pathology characteristics. However, detection of CLS was associated with hormone ratios. Compared with women in the highest tertile of estrone:androstenedione ratio in fat, those in the lowest tertile were less likely to have CLS (OR 0.12, 95% CI 0.03–0.59). A similar pattern was observed with estradiol:testosterone ratio in serum and CLS (lowest versus highest tertile, OR 0.18, 95% CI 0.04–0.72). CLS were more frequently identified in the breast fat of obese women and were associated with increased ratios of select estrogens:androgens in the blood and tissues, but not with individual hormones. Additional studies on CLS, tissue and blood hormone levels, and breast cancer risk are needed to understand and confirm these findings.
Archived formalin fixed paraffin embedded (FFPE) samples are valuable clinical resources to examine clinically relevant morphology features and also to study genetic changes. However, DNA quality and quantity of FFPE samples are often sub-optimal, and resulting NGS-based genetics variant detections are prone to false positives. Evaluations of wet-lab and bioinformatics approaches are needed to optimize variant detection from FFPE samples. As a pilot study, we designed within-subject triplicate samples of DNA derived from paired FFPE and fresh frozen breast tissues to highlight FFPE-specific artifacts. For FFPE samples, we tested two FFPE DNA extraction methods to determine impact of wet-lab procedures on variant calling: QIAGEN QIAamp DNA Mini Kit (“QA”), and QIAGEN GeneRead DNA FFPE Kit (“QGR”). We also used negative-control (NA12891) and positive control samples (Horizon Discovery Reference Standard FFPE). All DNA sample libraries were prepared for NGS according to the QIAseq Human Breast Cancer Targeted DNA Panel protocol and sequenced on the HiSeq 4000. Variant calling and filtering were performed using QIAGEN Gene Globe Data Portal. Detailed variant concordance comparisons and mutational signature analysis were performed to investigate effects of FFPE samples compared to paired fresh frozen samples, along with different DNA extraction methods. In this study, we found that five times or more variants were called with FFPE samples, compared to their paired fresh-frozen tissue samples even after applying molecular barcoding error-correction and default bioinformatics filtering recommended by the vendor. We also found that QGR as an optimized FFPE-DNA extraction approach leads to much fewer discordant variants between paired fresh frozen and FFPE samples. Approximately 92% of the uniquely called FFPE variants were of low allelic frequency range (< 5%), and collectively shared a “C > T|G > A” mutational signature known to be representative of FFPE artifacts resulting from cytosine deamination. Based on control samples and FFPE-frozen replicates, we derived an effective filtering strategy with associated empirical false-discovery estimates. Through this study, we demonstrated feasibility of calling and filtering genetic variants from FFPE tissue samples using a combined strategy with molecular barcodes, optimized DNA extraction, and bioinformatics methods incorporating genomics context such as mutational signature and variant allelic frequency.
Purpose: Programmed death ligand 1 [PD-(L)1]-targeted therapies have shown modest survival benefit in triple-negative breast cancer (TNBC). PD-L1+ microenvironments in TNBC are not well characterized and may inform combinatorial immune therapies. Herein, we characterized clinicopathologic features, RNA-based immune signatures, and spatially defined protein-based tumor–immune microenvironments (TIME) in early-stage PD-L1+ and PD-L1− TNBC.Experimental Design: From a large cohort of chemotherapy-naïve TNBC, clinicopathologic features, deconvoluted RNA immune signatures, and intraepithelial and stromal TIME (Nanostring GeoMX) were identified in subsets of PD-L1+ and PD-L1− TNBC, as defined by FDA-approved PD-L1 companion assays.Results: 228 of 499 (46%) TNBC were PD-L1+ (SP142: ≥1% immune cells-positive). Using PD-L1 22C3, 46% had combined positive score (CPS) ≥ 1 and 16% had CPS ≥10. PD-L1+ TNBC were higher grade with higher tumor-infiltrating lymphocytes (TIL; P < 0.05). PD-L1 was not associated with improved survival following adjustment for TILs and other variables. RNA profiles of PD-L1+ TNBC had increased dendritic cell, macrophage, and T/B cell subset features; and decreased myeloid-derived suppressor cells. PD-L1+ stromal and intraepithelial TIMEs were highly enriched in IDO-1, HLA-DR, CD40, and CD163 compared with PD-L1-TIME, with spatially specific alterations in CTLA-4, Stimulator of Interferon Genes (STING), and fibronectin. Macrophage- and antigen presentation–related proteins correlated most strongly with PD-L1 protein.Conclusions: In this early-stage TNBC cohort, nearly 50% were PD-L1+ (SP142 companion assay) while 16% were PD-L1+ with the 22C3 companion assay. PD-L1+ TNBC had specific myeloid-derived and lymphoid features. Spatially defined PD-L1+ TIME were enriched in several clinically actionable immune proteins. These data may inform future studies on combinatorial immunotherapies for patients with PD-L1+ TNBC.See related commentary by Symmans, p. 5446