Search

Background: Somatic mutational signatures elucidate molecular vulnerabilities to therapy, and therefore detecting signatures and classifying tumors with respect to signatures has clinical value. However, identifying the etiology of the mutational signatures remains a statistical challenge, with both small sample sizes and high variability in classification algorithms posing barriers. As a result, few signatures have been strongly linked to particular risk factors., Methods: Here, we develop a statistical model, Diffsig, for estimating the association of one or more continuous or categorical risk factors with DNA mutational signatures. Diffsig takes into account the uncertainty associated with assigning signatures to samples as well as multiple risk factors' simultaneous effect on observed DNA mutations., Results: We applied Diffsig to breast cancer data to assess relationships between five established breast-relevant mutational signatures and etiologic variables, confirming known mechanisms of cancer development. In simulation, our model was capable of accurately estimating expected associations in a variety of contexts., Conclusions: Diffsig allows researchers to quantify and perform inference on the associations of risk factors with mutational signatures., Impact: We expect Diffsig to provide more robust associations of risk factors with signatures to lead to better understanding of the tumor development process and improved models of tumorigenesis., (©2024 American Association for Cancer Research.)

Markers of genomic instability, including TP53 status and homologous recombination deficiency (HRD), are candidate biomarkers of immunogenicity and immune-mediated survival, but little is known about the distribution of these markers in large, population-based cohorts of racially diverse patients with breast cancer. In prior clinical trials, DNA-based approaches have been emphasized, but recent data suggest that RNA-based assessment can capture pathway differences conveniently and may be streamlined with other RNA-based genomic risk scores. Thus, we used RNA expression to study genomic instability (HRD and TP53 pathways) in context of the breast cancer immune microenvironment in three datasets (total n = 4,892), including 1,942 samples from the Carolina Breast Cancer Study, a population-based study that oversampled Black ( n = 1,026) and younger women ( n = 1,032). Across all studies, 36.9% of estrogen receptor (ER)-positive and 92.6% of ER-negative breast cancer had presence of at least one genomic instability signature. TP53 and HRD status were significantly associated with immune expression in both ER-positive and ER-negative breast cancer. RNA-based genomic instability signatures were associated with higher PD-L1, CD8 T-cell marker, and global and multimarker immune cell expression. Among tumors with genomic instability signatures, adaptive immune response was associated with improved recurrence-free survival regardless of ER status, highlighting genomic instability as a candidate marker for predicting immunotherapy response. Leveraging a convenient, integrated RNA-based approach, this analysis shows that genomic instability interacts with immune response, an important target in breast cancer overall and in Black women who experience higher frequency of TP53 and HR deficiency., Significance: Despite promising advances in breast cancer immunotherapy, predictive biomarkers that are valid across diverse populations and breast cancer subtypes are needed. Genomic instability signatures can be coordinated with other RNA-based scores to define immunogenic breast cancers and may have value in stratifying immunotherapy trial participants., Competing Interests: S.C. Van Alsten reports grants from NCI during the conduct of the study. B.C. Calhoun reports personal fees from Luminex Corp outside the submitted work. No disclosures were reported by the other authors., (© 2023 The Authors; Published by the American Association for Cancer Research.)

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., (© 2022. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

Background: Genome-wide association studies (GWAS) have identified more than 200 susceptibility loci for breast cancer, but these variants explain less than a fifth of the disease risk. Although gene-environment interactions have been proposed to account for some of the remaining heritability, few studies have empirically assessed this., Methods: We obtained genotype and risk factor data from 46,060 cases and 47,929 controls of European ancestry from population-based studies within the Breast Cancer Association Consortium (BCAC). We built gene expression prediction models for 4,864 genes with a significant (P<0.01) heritable component using the transcriptome and genotype data from the Genotype-Tissue Expression (GTEx) project. We leveraged predicted gene expression information to investigate the interactions between gene-centric genetic variation and 14 established risk factors in association with breast cancer risk, using a mixed-effects score test., Results: After adjusting for number of tests using Bonferroni correction, no interaction remained statistically significant. The strongest interaction observed was between the predicted expression of the C13orf45 gene and age at first full-term pregnancy (P GXE =4.44×10 -6 )., Conclusion: In this transcriptome-informed genome-wide gene-environment interaction study of breast cancer, we found no strong support for the role of gene expression in modifying the associations between established risk factors and breast cancer risk., Impact: Our study suggests a limited role of gene-environment interactions in breast cancer risk., Competing Interests: Conflict of interest disclosure statement: The authors declare no potential conflicts of interest.

Purpose: To describe breast cancer treatment patterns among premenopausal women by age and time since last pregnancy., Methods: Data were analyzed from 1179 women diagnosed with premenopausal breast cancer in the Carolina Breast Cancer Study. Of these, 160 had a recent pregnancy (within 5 years of cancer diagnosis). Relative frequency differences (RFDs) and 95% confidence intervals (CIs) were used to compare cancer stage, treatment modality received, treatment initiation delay (> 30 days), and prolonged treatment duration (> 2 to > 8 months depending on the treatment received) by age and recency of pregnancy., Results: Recently postpartum women were significantly more likely to have stage III disease [RFD (95% CI) 12.2% (3.6%, 20.8%)] and to receive more aggressive treatment compared to nulliparous women. After adjustment for age, race and standard clinical tumor characteristics, recently postpartum women were significantly less likely to have delayed treatment initiation [RFD (95% CI) - 11.2% (- 21.4%, - 1.0%)] and prolonged treatment duration [RFD (95% CI) - 17.5% (- 28.0%, - 7.1%)] and were more likely to have mastectomy [RFD (95% CI) 14.9% (4.8%, 25.0%)] compared to nulliparous. Similarly, younger women (< 40 years of age) were significantly less likely to experience prolonged treatment duration [RFD (95% CI) - 5.6% (- 11.1%, - 0.0%)] and more likely to undergo mastectomy [RFD (95% CI) 10.6% (5.2%, 16.0%)] compared to the study population as a whole., Conclusion: These results suggest that recently postpartum and younger women often received prompt and aggressive breast cancer treatment. Higher mortality and recurrence among recently pregnant women are unlikely to be related to undertreatment., (© 2022. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

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., (©2021 American Association for Cancer Research.)

Background: TP53 and estrogen receptor (ER) both play essential roles in breast cancer development and progression, with recent research revealing cross-talk between TP53 and ER signaling pathways. Although many studies have demonstrated heterogeneity of risk factor associations across ER subtypes, associations by TP53 status have been inconsistent., Methods: This case-case analysis included incident breast cancer cases (47% Black) from the Carolina Breast Cancer Study (1993-2013). Formalin-fixed paraffin-embedded tumor samples were classified for TP53 functional status (mutant-like/wild-type-like) using a validated RNA signature. For IHC-based TP53 status, mutant-like was classified as at least 10% positivity. We used two-stage polytomous logistic regression to evaluate risk factor heterogeneity due to RNA-based TP53 and/or ER, adjusting for each other and for PR, HER2, and grade. We then compared this with the results when using IHC-based TP53 classification., Results: The RNA-based classifier identified 55% of tumors as TP53 wild-type-like and 45% as mutant-like. Several hormone-related factors (oral contraceptive use, menopausal status, age at menopause, and pre- and postmenopausal body mass index) were associated with TP53 mutant-like status, whereas reproductive factors (age at first birth and parity) and smoking were associated with ER status. Multiparity was associated with both TP53 and ER. When classifying TP53 status using IHC methods, no associations were observed with TP53. Associations observed with RNA-based TP53 remained after accounting for basal-like subtype., Conclusions: This case-case study found breast cancer risk factors associated with RNA-based TP53 and ER., Impact: RNA-based TP53 and ER represent an emerging etiologic schema of interest in breast cancer prevention research., (©2021 The Authors; Published by the American Association for Cancer Research.)

Continuous risk of recurrence scores (CRS) based on tumor gene expression are vital prognostic tools for breast cancer. Studies have shown that Black women (BW) have higher CRS than White women (WW). Although systemic injustices contribute substantially to breast cancer disparities, evidence of biological and germline contributions is emerging. In this study, we investigated germline genetic associations with CRS and CRS disparity using approaches modeled after transcriptome-wide association studies (TWAS). In the Carolina Breast Cancer Study, using race-specific predictive models of tumor expression from germline genetics, we performed race-stratified ( N = 1,043 WW, 1,083 BW) linear regressions of three CRS (ROR-S: PAM50 subtype score; proliferation score; ROR-P: ROR-S plus proliferation score) on imputed tumor genetically regulated tumor expression (GReX). Bayesian multivariate regression and adaptive shrinkage tested GReX-prioritized genes for associations with tumor PAM50 expression and subtype to elucidate patterns of germline regulation underlying GReX-CRS associations. At FDR-adjusted P < 0.10, 7 and 1 GReX prioritized genes among WW and BW, respectively. Among WW, CRS were positively associated with MCM10 , FAM64A , CCNB2 , and MMP1 GReX and negatively associated with VAV3 , PCSK6 , and GNG11 GReX. Among BW, higher MMP1 GReX predicted lower proliferation score and ROR-P. GReX-prioritized gene and PAM50 tumor expression associations highlighted potential mechanisms for GReX-prioritized gene to CRS associations. Among patients with breast cancer, differential germline associations with CRS were found by race, underscoring the need for larger, diverse datasets in molecular studies of breast cancer. These findings also suggest possible germline trans -regulation of PAM50 tumor expression, with potential implications for CRS interpretation in clinical settings. SIGNIFICANCE: This study identifies race-specific genetic associations with breast cancer risk of recurrence scores and suggests mediation of these associations by PAM50 subtype and expression, with implications for clinical interpretation of these scores., (©2021 The Authors; Published by the American Association for Cancer Research.)

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., (© 2021. The Author(s).)

The NanoString RNA counting assay for formalin-fixed paraffin embedded samples is unique in its sensitivity, technical reproducibility and robustness for analysis of clinical and archival samples. While commercial normalization methods are provided by NanoString, they are not optimal for all settings, particularly when samples exhibit strong technical or biological variation or where housekeeping genes have variable performance across the cohort. Here, we develop and evaluate a more comprehensive normalization procedure for NanoString data with steps for quality control, selection of housekeeping targets, normalization and iterative data visualization and biological validation. The approach was evaluated using a large cohort ($N=\kern0.5em 1649$) from the Carolina Breast Cancer Study, two cohorts of moderate sample size ($N=359$ and$130$) and a small published dataset ($N=12$). The iterative process developed here eliminates technical variation (e.g. from different study phases or sites) more reliably than the three other methods, including NanoString's commercial package, without diminishing biological variation, especially in long-term longitudinal multiphase or multisite cohorts. We also find that probe sets validated for nCounter, such as the PAM50 gene signature, are impervious to batch issues. This work emphasizes that systematic quality control, normalization and visualization of NanoString nCounter data are an imperative component of study design that influences results in downstream analyses., (© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

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., (© The Author(s) 2020. Published by Oxford University Press.)

Background: The relationship between germline genetic variation and breast cancer survival is largely unknown, especially in understudied minority populations who often have poorer survival. Genome-wide association studies (GWAS) have interrogated breast cancer survival but often are underpowered due to subtype heterogeneity and clinical covariates and detect loci in non-coding regions that are difficult to interpret. Transcriptome-wide association studies (TWAS) show increased power in detecting functionally relevant loci by leveraging expression quantitative trait loci (eQTLs) from external reference panels in relevant tissues. However, ancestry- or race-specific reference panels may be needed to draw correct inference in ancestrally diverse cohorts. Such panels for breast cancer are lacking., Results: We provide a framework for TWAS for breast cancer in diverse populations, using data from the Carolina Breast Cancer Study (CBCS), a population-based cohort that oversampled black women. We perform eQTL analysis for 406 breast cancer-related genes to train race-stratified predictive models of tumor expression from germline genotypes. Using these models, we impute expression in independent data from CBCS and TCGA, accounting for sampling variability in assessing performance. These models are not applicable across race, and their predictive performance varies across tumor subtype. Within CBCS (N = 3,828), at a false discovery-adjusted significance of 0.10 and stratifying for race, we identify associations in black women near AURKA, CAPN13, PIK3CA, and SERPINB5 via TWAS that are underpowered in GWAS., Conclusions: We show that carefully implemented and thoroughly validated TWAS is an efficient approach for understanding the genetics underpinning breast cancer outcomes in diverse populations.

Background: The dominant invasive breast cancer histologic subtype, ductal carcinoma, shows intrinsic subtype diversity. However, lobular breast cancers are predominantly Luminal A. Both histologic subtypes show distinct relationships with patient and tumor characteristics, but it is unclear if these associations remain after accounting for intrinsic subtype., Methods: Generalized linear models were used to estimate relative frequency differences (RFDs) and 95% confidence intervals (95% CIs) for the associations between age, race, tumor characteristics, immunohistochemistry (IHC) and RNA-based intrinsic subtype, TP53 status, and histologic subtype in the Carolina Breast Cancer Study (CBCS, n = 3,182) and The Cancer Genome Atlas (TCGA, n = 808)., Results: Relative to ductal tumors, lobular tumors were significantly more likely to be Luminal A [CBCS RNA RFD: 44.9%, 95% CI (39.6, 50.1); TCGA: RFD: 50.5%, 95% CI (43.9, 57.1)], were less frequent among young (≤ 50 years) and black women, were larger in size, low grade, less frequently had TP53 pathway defects, and were diagnosed at later stages. These associations persisted among Luminal A tumors (n = 242)., Conclusions: While histology is strongly associated with molecular characteristics, histologic associations with age, race, size, grade, and stage persisted after restricting to Luminal A subtype. Histology may continue to be clinically relevant among Luminal A breast cancers.

Background: Invasive lobular breast tumors display unique reproductive risk factor profiles. Lobular tumors are predominantly Luminal A subtype, and it is unclear whether reported risk factor associations are independent of molecular subtype., Methods: Polytomous logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (95% CIs) for the associations between risk factors and histologic subtype [ductal (n = 2,856), lobular (n = 326), and mixed ductal-lobular (n = 473)] in the Carolina Breast Cancer Study (1993-2013). Three-marker immunohistochemical clinical subtypes were defined as Luminal A (ER+ or PR+/HER2-), Luminal B (ER+ or PR+/HER2+), Triple Negative (ER-/PR-/HER2-), and HER2+ (ER-/PR-/HER2+)., Results: In case-case analyses compared to ductal, lobular tumors were significantly associated with lactation duration > 12 months [OR 1.86, 95% CI (1.33-2.60)], age at first birth ≥ 26 years [OR: 1.35, 95% CI: (1.03-1.78)], and current oral contraceptive use [OR: 1.86, 95% CI: (1.08-3.20)]. Differences in risk factor associations between ductal and lobular tumors persisted after restricting to Luminal A subtype., Conclusions: Lobular tumors were associated with older age at first birth, increased lactation duration, and current oral contraceptive use. Etiologic heterogeneity by histology persisted after restricting to Luminal A subtype, suggesting both tumor histology and intrinsic subtype play integral parts in breast cancer risk.