Differences in 21-Gene and PAM50 Recurrence Scores in Younger and Black Women With Breast Cancer.

PURPOSE: Genomic tests, such as the Oncotype Dx 21-gene and Prosigna risk of recurrence (ROR-P) assay, are commonly used for breast cancer prognostication. Emerging data suggest variability between assays, but this has not been compared in diverse populations. MATERIALS AND METHODS: RNA sequencing was performed on 647 previously untreated stage I-III estrogen receptor-positive/human epidermal growth factor receptor 2-negative tumors in the Carolina Breast Cancer Study, which oversampled Black and younger women (age <50 years at diagnosis), using research versions of two common RNA-based prognostic assays: ROR-PR and the 21-gene recurrence score (RSR). Relative frequency differences and 95% CIs were estimated for associations with race and age, and hazards of 5-year local or distant recurrence were modeled with Cox regression. Proliferation and estrogen module scores from each assay, representing broad activity of genes in those pathways, were examined to guide interpretation of differences between tests. RESULTS: Among both younger and older individuals, Black women had higher frequency of intermediate and high ROR-PR scores than non-Black women. Race was not significantly associated with RSR in either age group. High (hazard ratio [HR], 4.67 [95% CI, 1.73 to 12.70]) and intermediate (HR, 2.12 [95% CI, 0.98 to 4.62]) ROR-PR scores were associated with greater risk of recurrence, but RSR did not predict recurrence. RSR emphasized estrogen over proliferation modules, whereas ROR-PR emphasized proliferation. Higher proliferation scores were associated with younger age and Black race in both assays. Modifications to the RSR algorithm that increased emphasis on proliferation improved prognostication in this diverse population. CONCLUSION: ROR-PR and the 21-gene RSR differentially emphasize estrogen-related and proliferative biology. The emphasis of 21-gene RS on estrogen-related biology and lower endocrine therapy initiation among Black women may contribute to poorer prognostic ability in heterogeneously treated populations.

BIRC5 expression by race, age and clinical factors in breast cancer patients.

PURPOSE: Survivin/BIRC5 is a proliferation marker that is associated with poor prognosis in breast cancer and an attractive therapeutic target. However, BIRC5 has not been well studied among racially diverse populations where aggressive breast cancers are prevalent. EXPERIMENTAL DESIGN: We studied BIRC5 expression in association with clinical and demographic variables and as a predictor of recurrence in 2174 participants in the Carolina Breast Cancer Study (CBCS), a population-based study that oversampled Black (n = 1113) and younger (< 50 years; n = 1137) participants with breast cancer. For comparison, similar analyses were conducted in The Cancer Genome Atlas [TCGA N = 1094, Black (n = 183), younger (n = 295)]. BIRC5 was evaluated as a continuous and categorical variable (highest quartile vs. lower three quartiles). RESULTS: Univariate, continuous BIRC5 expression was higher in breast tumors from Black women relative to non-Black women in both estrogen receptor (ER)-positive and ER-negative tumors and in analyses stratified by stage (i.e., within Stage I, Stage II, and Stage III/IV tumors). Within CBCS and TCGA, BIRC5-high was associated with young age (< 50 years) and Black race, as well as hormone receptor-negative tumors, non-Luminal A PAM50 subtypes, advanced stage, and larger tumors (> 2 cm). Relative to BIRC5-low, BIRC5-high tumors were associated with poor 5-year recurrence-free survival (RFS) among ER-positive tumors, both in unadjusted models [HR (95% CI): 2.7 (1.6, 4.6)] and after adjustment for age and stage [Adjusted HR (95% CI): 1.87 (1.07, 3.25)]. However, this relationship was not observed among ER-negative tumors [Crude HR (95% CI): 0.7 (0.39, 1.2); Adjusted HR (95% CI): 0.67 (0.37, 1.2)]. CONCLUSION: Black and younger women with breast cancer have a higher burden of BIRC5-high tumors than older and non-Black women. Emerging anti-survivin treatment strategies may be an important future direction for equitable breast cancer outcomes.

Disparities in OncotypeDx Testing and Subsequent Chemotherapy Receipt by Geography and Socioeconomic Status.

BACKGROUND: OncotypeDx is a prognostic and predictive genomic assay used in early-stage hormone receptor-positive, HER2- (HR+/HER2-) breast cancer. It is used to inform adjuvant chemotherapy decisions, but not all eligible women receive testing. We aimed to assess variation in testing by demographics and geography, and to determine whether testing was associated with chemotherapy. METHODS: For 1,615 women in the Carolina Breast Cancer Study with HR+/HER2-, Stage I-II tumors, we estimated prevalence differences (PD) and 95% confidence intervals (CI) for receipt of OncotypeDx genomic testing in association with and sociodemographic characteristics. We assessed associations between testing and chemotherapy receipt overall and by race. Finally, we calculated the proportion of eligible women receiving OncotypeDx by county-level rurality, census tract-level socioeconomic status, and Area Health Education Center regions. RESULTS: 38% (N = 609) of potentially eligible women were tested, with lower testing prevalences in Black (31%; PD, -11%; 95% CI, -16%-6%) and low-income women (24%; PD, -20%; 95% CI, -29% to -11%) relative to non-Black and higher income women. Urban participants were less likely to be tested than rural participants, though this association varied by region. Among women with low genomic risk tumors, tested participants were 29% less likely to receive chemotherapy than untested participants (95% CI, -40% to -17%). Racial differences in chemotherapy were restricted to untested women. CONCLUSIONS: Both individual and area-level socioeconomics predict likelihood of OncotypeDx testing. IMPACT: Variable adoption of OncotypeDx by socioeconomics and across geographic settings may contribute to excess chemotherapy among patients with HR+/HER2- cancers. See related In the Spotlight, p. 635.

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.

Gene Expression Profiling Identifies Two Chordoma Subtypes Associated with Distinct Molecular Mechanisms and Clinical Outcomes.

PURPOSE: Chordoma is a rare bone tumor with a high recurrence rate and limited treatment options. The aim of this study was to identify molecular subtypes of chordoma that may improve clinical management. EXPERIMENTAL DESIGN: We conducted RNA sequencing in 48 tumors from patients with Chinese skull-base chordoma and identified two major molecular subtypes. We then replicated the classification using a NanoString panel in 48 patients with chordoma from North America. RESULTS: Tumors in one subtype were more likely to have somatic mutations and reduced expression in chromatin remodeling genes, such as PBRM1 and SETD2, whereas the other subtype was characterized by the upregulation of genes in epithelial-mesenchymal transition and Sonic Hedgehog pathways. IHC staining of top differentially expressed genes between the two subtypes in 312 patients with Chinese chordoma with long-term follow-up data showed that the expression of some markers such as PTCH1 was significantly associated with survival outcomes. CONCLUSIONS: Our findings may improve the understanding of subtype-specific tumorigenesis of chordoma and inform clinical prognostication and targeted options.

RNA-Based Classification of Homologous Recombination Deficiency in Racially Diverse Patients with Breast Cancer.

BACKGROUND: Aberrant expression of DNA repair pathways such as homologous recombination (HR) can lead to DNA repair imbalance, genomic instability, and altered chemotherapy response. DNA repair imbalance may predict prognosis, but variation in DNA repair in diverse cohorts of breast cancer patients is understudied. METHODS: To identify RNA-based patterns of DNA repair expression, we performed unsupervised clustering on 51 DNA repair-related genes in the Cancer Genome Atlas Breast Cancer [TCGA BRCA (n = 1,094)] and Carolina Breast Cancer Study [CBCS (n = 1,461)]. Using published DNA-based HR deficiency (HRD) scores (high-HRD ≥ 42) from TCGA, we trained an RNA-based supervised classifier. Unsupervised and supervised HRD classifiers were evaluated in association with demographics, tumor characteristics, and clinical outcomes. RESULTS: : Unsupervised clustering on DNA repair genes identified four clusters of breast tumors, with one group having high expression of HR genes. Approximately 39.7% of CBCS and 29.3% of TCGA breast tumors had this unsupervised high-HRD (U-HRD) profile. A supervised HRD classifier (S-HRD) trained on TCGA had 84% sensitivity and 73% specificity to detect HRD-high samples. Both U-HRD and S-HRD tumors in CBCS had higher frequency of TP53 mutant-like status (45% and 41% enrichment) and basal-like subtype (63% and 58% enrichment). S-HRD high was more common among black patients. Among chemotherapy-treated participants, recurrence was associated with S-HRD high (HR: 2.38, 95% confidence interval = 1.50-3.78). CONCLUSIONS: HRD is associated with poor prognosis and enriched in the tumors of black women. IMPACT: RNA-level indicators of HRD are predictive of breast cancer outcomes in diverse populations.

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 targetable interventions.

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.

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.

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.

Identification of a Novel Inflamed Tumor Microenvironment Signature as a Predictive Biomarker of Bacillus Calmette-Guérin Immunotherapy in Non-Muscle-Invasive Bladder Cancer.

PURPOSE: Improved risk stratification and predictive biomarkers of treatment response are needed for non-muscle-invasive bladder cancer (NMIBC). Here we assessed the clinical utility of targeted RNA and DNA molecular profiling in NMIBC. EXPERIMENTAL DESIGN: Gene expression in NMIBC samples was profiled by NanoString nCounter, an RNA quantification platform, from two independent cohorts (n = 28, n = 50); targeted panel sequencing was performed in a subgroup (n = 50). Gene signatures were externally validated using two RNA sequencing datasets of NMIBC tumors (n = 438, n = 73). Established molecular subtype classifiers and novel gene expression signatures were assessed for associations with clinicopathologic characteristics, somatic tumor mutations, and treatment outcomes. RESULTS: Molecular subtypes distinguished between low-grade Ta tumors with FGFR3 mutations and overexpression (UROMOL-class 1) and tumors with more aggressive clinicopathologic characteristics (UROMOL-classes 2 and 3), which were significantly enriched with TERT promoter mutations. However, UROMOL subclasses were not associated with recurrence after bacillus Calmette-Guérin (BCG) immunotherapy in two independent cohorts. In contrast, a novel expression signature of an inflamed tumor microenvironment (TME) was associated with improved recurrence-free survival after BCG. Expression of immune checkpoint genes (PD-L1/PD-1/CTLA-4) was associated with an inflamed TME, but not with higher recurrence rates after BCG. FGFR3 mutations and overexpression were both associated with low immune signatures. CONCLUSIONS: Assessment of the immune TME, rather than molecular subtypes, is a promising predictive biomarker of BCG response. Modulating the TME in an immunologically "cold" tumor warrants further investigation. Integrated transcriptomic and exome sequencing should improve treatment selection in NMIBC.

Differences in risk factors for molecular subtypes of clear cell renal cell carcinoma.

The ccA and ccB molecular subtypes of clear cell renal cell carcinoma (ccRCC) have well-characterized prognostic relevance. However, it is not known whether they possess distinct etiologies. We investigated the relationships between these subtypes and RCC risk factors within a case-control study conducted in Eastern Europe. We analyzed risk factor data for ccA (n = 144) and ccB (n = 106) cases and 1476 controls through case-only and case-control comparisons to assess risk factor differences across subtypes using logistic and polytomous regression models. We also performed a meta-analysis summarizing case-only results from our study and three patient cohorts. Patients with ccB tumors had poorer survival than those with ccA tumors and were more likely to be male (case-only odds ratio [OR] 2.68, 95% confidence interval [CI] 1.43-5.03). In case-control analyses, body mass index was significantly associated with ccA tumors (OR 2.45, 95% CI 1.18-5.10 for ≥35 vs <25 kg/m2 ) but not with ccB tumors (1.52, 0.56-4.12), while trichloroethylene was associated with ccB but not ccA (OR 3.09, 95% CI 1.11-8.65 and 1.25, 0.36-4.39 respectively for ≥1.58 ppm-years vs unexposed). A polygenic risk score of genetic variants identified from genome-wide association studies was associated with both ccA and, in particular, ccB (OR 1.82, 1.11-2.99 and 2.87, 95% CI 1.64-5.01 respectively for 90th vs 10th percentile). In a meta-analysis of case-only results including three patient cohorts, we still observed the ccB excess for male sex and the ccA excess for obesity. In conclusion, our findings suggest the existence of etiologic heterogeneity across ccRCC molecular subtypes for several risk factors.

Protein-based immune profiles of basal-like vs. luminal breast cancers.

Tumor-infiltrating lymphocytes play an important, but incompletely understood role in chemotherapy response and prognosis. In breast cancer, there appear to be distinct immune responses by subtype, but most studies have used limited numbers of protein markers or bulk sequencing of RNA to characterize immune response, in which spatial organization cannot be assessed. To identify immune phenotypes of Basal-like vs. Luminal breast cancer we used the GeoMx® (NanoString) platform to perform digital spatial profiling of immune-related proteins in tumor whole sections and tissue microarrays (TMA). Visualization of CD45, CD68, or pan-Cytokeratin by immunofluorescence was used to select regions of interest in formalin-fixed paraffin embedded tissue sections. Forty-four antibodies representing stromal markers and multiple immune cell types were applied to quantify the tumor microenvironment. In whole tumor slides, immune hot spots (CD45+) had increased expression of many immune markers, suggesting a diverse and robust immune response. In epithelium-enriched areas, immune signals were also detectable and varied by subtype, with regulatory T-cell (Treg) markers (CD4, CD25, and FOXP3) being higher in Basal-like vs. Luminal breast cancer. Extending these findings to TMAs with more patients (n = 75), we confirmed subtype-specific immune profiles, including enrichment of Treg markers in Basal-likes. This work demonstrated that immune responses can be detected in epithelium-rich tissue, and that TMAs are a viable approach for obtaining important immunoprofiling data. In addition, we found that immune marker expression is associated with breast cancer subtype, suggesting possible prognostic, or targetable differences.

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.

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.

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 mm2, 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 mm2, ORT3vsT1: 3.42, 95% confidence interval (CI), 1.87-6.28; acini count/TDLUT3vsT1, 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-7). 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.

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.