Olshan, Andrew F.

Background Polygenic risk score (PRS), calculated based on genome-wide association studies (GWASs), can improve breast cancer (BC) risk assessment. To date, most BC GWASs have been performed in individuals of European (EUR) ancestry, and the generalisation of EUR-based PRS to other populations is a major challenge. In this study, we examined the performance of EUR-based BC PRS models in Ashkenazi Jewish (AJ) women. Methods We generated PRSs based on data on EUR women from the Breast Cancer Association Consortium (BCAC). We tested the performance of the PRSs in a cohort of 2161 AJ women from Israel (1437 cases and 724 controls) from BCAC (BCAC cohort from Israel (BCAC-IL)). In addition, we tested the performance of these EUR-based BC PRSs, as well as the established 313-SNP EUR BC PRS, in an independent cohort of 181 AJ women from Hadassah Medical Center (HMC) in Israel. Results In the BCAC-IL cohort, the highest OR per 1 SD was 1.56 (±0.09). The OR for AJ women at the top 10% of the PRS distribution compared with the middle quintile was 2.10 (±0.24). In the HMC cohort, the OR per 1 SD of the EUR-based PRS that performed best in the BCAC-IL cohort was 1.58±0.27. The OR per 1 SD of the commonly used 313-SNP BC PRS was 1.64 (±0.28). Conclusions Extant EUR GWAS data can be used for generating PRSs that identify AJ women with markedly elevated risk of BC and therefore hold promise for improving BC risk assessment in AJ women.

Abstract Background Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC‐specific survival (BCSS) compared to non‐carriers. Aim To assessed the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. Methods Analyses were based on 82,701 women diagnosed with a first primary invasive BC including 963 CHEK2 c.1100delC carriers; median follow‐up was 9.1 years. Differential associations with treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi‐state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. Results There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status. The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR (95% CI): 0.66 (0.55–0.78)]. No association was observed with radiotherapy. Results from the multi‐state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non‐carriers also after accounting for CBC occurrence [HR (95% CI): 1.30 (1.09–1.56)]. Conclusion Systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk.

Abstract Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers. We aimed to assess the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. Analyses were based on 82,701 women diagnosed with invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations of treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR(95%CI): 0.66 (0.55–0.78)]. No association was observed with radiotherapy. Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR(95%CI) :1.30 (1.09–1.56)]. In conclusion, systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk. (Main MS: 3201 words)

AbstractBackgroundFertility preservation (FP) may be underused after cancer diagnosis because of uncertainty around delays to cancer treatment and subsequent reproductive success.MethodsWomen aged 15 to 39 years diagnosed with cancer between 2004 and 2015 were identified from the North Carolina Central Cancer Registry. Use of assisted reproductive technology (ART) after cancer diagnosis between 2004 and 2018 (including FP) was assessed through linkage to the Society for Assisted Reproductive Technology. Linear regression was used to examine time to cancer treatment among women who did (n = 95) or did not (n = 469) use FP. Modified Poisson regression was used to estimate risk ratios (RRs) and 95% CIs for pregnancy and birth based on timing of ART initiation relative to cancer treatment (n = 18 initiated before treatment for FP vs n = 26 initiated after treatment without FP).ResultsThe median time to cancer treatment was 9 to 33 days longer among women who used FP compared with women who did not, matched on clinical factors. Women who initiated ART before cancer treatment may be more likely to have a live birth given pregnancy compared with women who initiated ART after cancer treatment (age‐adjusted RR, 1.47; 95% CI, 0.98‐2.23), though this may be affected by the more frequent use of gestational carriers in the former group (47% vs 20% of transfer cycles, respectively).ConclusionsFP delayed gonadotoxic cancer treatment by up to 4.5 weeks, a delay that would not be expected to alter prognosis for many women. Further study of the use of gestational carriers in cancer populations is warranted to better understand its effect on reproductive outcomes.

AbstractBackgroundReproductive factors have been shown to be differentially associated with risk of estrogen receptor (ER)-positive and ER-negative breast cancer. However, their associations with intrinsic-like subtypes are less clear.MethodsAnalyses included up to 23 353 cases and 71 072 controls pooled from 31 population-based case-control or cohort studies in the Breast Cancer Association Consortium across 16 countries on 4 continents. Polytomous logistic regression was used to estimate the association between reproductive factors and risk of breast cancer by intrinsic-like subtypes (luminal A-like, luminal B-like, luminal B-HER2–like, HER2-enriched–like, and triple-negative breast cancer) and by invasiveness. All statistical tests were 2-sided.ResultsCompared with nulliparous women, parous women had a lower risk of luminal A-like, luminal B-like, luminal B-HER2–like, and HER2-enriched–like disease. This association was apparent only after approximately 10 years since last birth and became stronger with increasing time (odds ratio [OR] = 0.59, 95% confidence interval [CI] = 0.49 to 0.71; and OR = 0.36, 95% CI = 0.28 to 0.46 for multiparous women with luminal A-like tumors 20 to less than 25 years after last birth and 45 to less than 50 years after last birth, respectively). In contrast, parous women had a higher risk of triple-negative breast cancer right after their last birth (for multiparous women: OR = 3.12, 95% CI = 2.02 to 4.83) that was attenuated with time but persisted for decades (OR = 1.03, 95% CI = 0.79 to 1.34, for multiparous women 25 to less than 30 years after last birth). Older age at first birth (Pheterogeneity < .001 for triple-negative compared with luminal A-like breast cancer) and breastfeeding (Pheterogeneity < .001 for triple-negative compared with luminal A-like breast cancer) were associated with lower risk of triple-negative breast cancer but not with other disease subtypes. Younger age at menarche was associated with higher risk of all subtypes; older age at menopause was associated with higher risk of luminal A-like but not triple-negative breast cancer. Associations for in situ tumors were similar to luminal A-like.ConclusionsThis large and comprehensive study demonstrates a distinct reproductive risk factor profile for triple-negative breast cancer compared with other subtypes, with implications for the understanding of disease etiology and risk prediction.

Abstract 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

Abstract Background: Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. Methods: We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values<5x10-8 as genome-wide significant, and p-values<1x10-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants.Results: None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values<1x105. The strongest evidence was found for rs4674019 (p-value=2.27x10-7), which showed genome-wide significant interaction (p-value=3.8x10-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. Conclusions: In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT–breast cancer risk association.

Abstract Background: Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. Methods: We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values<5x10-8 as genome-wide significant, and p-values<1x10-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants.Results: None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values<1x105. The strongest evidence was found for rs4674019 (p-value=2.27x10-7), which showed genome-wide significant interaction (p-value=3.8x10-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. Conclusions: In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT–breast cancer risk association.

Abstract Background: It is not known whether modifiable lifestyle factors that predict survival after invasive breast cancer differ by subtype. Methods: We analyzed data for 121,435 women diagnosed with breast cancer from 67 studies in the Breast Cancer Association Consortium with 16,890 deaths (8,554 breast cancer specific) over 10 years. Cox regression was used to estimate associations between risk factors and 10-year all-cause mortality and breast cancer–specific mortality overall, by estrogen receptor (ER) status, and by intrinsic-like subtype. Results: There was no evidence of heterogeneous associations between risk factors and mortality by subtype (Padj &gt; 0.30). The strongest associations were between all-cause mortality and BMI ≥30 versus 18.5–25 kg/m2 [HR (95% confidence interval (CI), 1.19 (1.06–1.34)]; current versus never smoking [1.37 (1.27–1.47)], high versus low physical activity [0.43 (0.21–0.86)], age ≥30 years versus &lt;20 years at first pregnancy [0.79 (0.72–0.86)]; &gt;0–&lt;5 years versus ≥10 years since last full-term birth [1.31 (1.11–1.55)]; ever versus never use of oral contraceptives [0.91 (0.87–0.96)]; ever versus never use of menopausal hormone therapy, including current estrogen–progestin therapy [0.61 (0.54–0.69)]. Similar associations with breast cancer mortality were weaker; for example, 1.11 (1.02–1.21) for current versus never smoking. Conclusions: We confirm associations between modifiable lifestyle factors and 10-year all-cause mortality. There was no strong evidence that associations differed by ER status or intrinsic-like subtype. Impact: Given the large dataset and lack of evidence that associations between modifiable risk factors and 10-year mortality differed by subtype, these associations could be cautiously used in prognostication models to inform patient-centered care.

Abstract Background: Persons of African ancestry (AA) experience a 1.5-2-fold risk of multiple myeloma (MM) compared to persons of European ancestry (EA). We assembled a set of MM patients with self-reported AA in order to evaluate the contribution of genetics to etiology in this high-risk group. Methods: Here we present the results of a meta-analysis of two GWAS in 1,813 cases and 8,871 controls of AA. We also conducted an admixture mapping scan to identify risk alleles associated with local ancestry, fine-mapped the 23 known susceptibility loci to find markers that could better capture MM risk in individuals of AA, and constructed a polygenic risk score (PRS) to assess the aggregated effect of known MM risk alleles. Finally, we conducted an eQTL analysis measuring gene expression in those genes harboring a risk variant in malignant plasma cells from 292 of the patients from a single site. Results: In GWAS analysis, we identified two suggestive novel loci located at 9p24.3 and 9p13.1 at P&lt;1 × 10-6, but no genome-wide significant association was noted. In admixture mapping, we observed a genome-wide significant inverse association between local AA at 2p24.1-23.1 and MM risk in AA individuals. 20 of the 23 known EA risk variants showed directional consistency and 9 replicated at P&lt;0.05 in AA individuals. In eight regions, we identified markers that better capture MM risk in persons of AA. AA individuals with a PRS in the top 10% had a 1.82-fold (95%CI: 1.56, 2.11) increased MM risk compared to those with average risk (25-75%). The strongest functional association was between the risk allele for variant rs56219066 at 5q15 and lower ELL2 expression (P= 5.1 × 10–12). Conclusion: Our study shows that common genetic variation contributes to MM risk individuals of AA. This abstract is also being presented as Poster C040. Citation Format: Zhaohui Du, Niels Weinhold, Gregory Chi Song, Kristen A. Rand, David J. Van Den Berg, Amie E. Hwang, Xin Sheng, Victor Hom, Sikander Ailawadhi, Ajay K. Nooka, Seema Singhal, Karen Pawlish, Edward Peters, Cathryn Bock, Ann Mohrbacher, Alexander Stram, Sonja I. Berndt, William J. Blot, Graham Casey, Victoria L. Stevens, Rick Kittles, Phyllis J. Goodman, W. Ryan Diver, Anselm Hennis, Barbara Nemesure, Eric A. Klein, Benjamin A. Rybicki, Janet L. Stanford, John S. Witte, Lisa Signorello, Esther M. John, Leslie Bernstein, Antoinette Stroup, Owen W. Stephens, Maurizio Zangari, Frits Van Rhee, Andrew Olshan, Wei Zheng, Jennifer J. Hu, Regina Ziegler, Sarah J. Nyante, Sue Ann Ingles, Michael Press, John David Carpten, Stephen Chanock, Jayesh Mehta, Graham A Colditz, Jeffrey Wolf, Thomas G. Martin, Michael Tomasson, Mark A. Fiala, Howard Terebelo, Nalini Janakiraman, Laurence Kolonel, Kenneth C. Anderson, Loic Le Marchand, Daniel Auclair, Brian C.-H. Chiu, Elad Ziv, Daniel Stram, Ravi Vij, Leon Bernal-Mizrachi, Gareth J. Morgan, Jeffrey A. Zonder, Carol Ann Huff, Sagar Lonial, Robert Z. Orlowski, David V. Conti, Christopher A. Haiman, Wendy Cozen. A meta-analysis of genome-wide association study and eQTL analysis of multiple myeloma among African Americans [abstract]. In: Proceedings of the Twelfth AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2019 Sep 20-23; San Francisco, CA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2020;29(6 Suppl_2):Abstract nr PR05.