Investigators, kConFab

Abstract Background Breast cancer is etiologically heterogeneous, but which risk factors differ in their associations across tumor subtypes remains unclear. We conducted a large, pooled analysis to evaluate independent, dose-response associations between breast cancer risk factors and quantitative tumor features. Methods Analyses of 15,731 invasive breast cancers from 24 studies evaluated associations (p-trend) between reproductive and hormonal factors, body mass index (BMI), alcohol, smoking, and family history in relation to quantitative immunohistochemistry measures on tissue microarrays (ER, PR, HER2, KI67, TP53) and tumor grade. Analyses in a subset of 10 population-based studies estimated subtype-specific odds ratios (ORs) comparing cases to controls. A Bayesian False Discovery Probability (BFDP) <0.2 was used to identify associations with strong statistical evidence. Results Nulliparity and later age at menopause were associated with higher ER-positivity (p-trend=0.021 and 0.001, respectively), with corresponding OR[ER+] (95% CI) = 1.49 (1.16-1.90) for nulliparous vs. parous and 1.07 (1.03-1.11) per 5 years. Current combined menopausal hormone therapy (MHT) use was associated with lower grade (p-trend<0.001), with OR [grade1] = 3.37 (2.69-4.21) for current vs. never users. Higher BMI was associated with lower ER-positivity and higher grade in premenopausal women (p-trend<0.001 and <0.001), with OR[ER+] = 0.80 (0.74-0.87) and OR[grade1] = 0.75 (0.63-0.88) per 5 units, and with higher PR-positivity and higher grade in postmenopausal women (p-trend<0.001 and <0.001), with OR[PR+] = 1.08 (1.03-1.14) and OR[grade 3] = 1.10 (1.04-1.17) per 5 units. Conclusion This pooled analysis of 15,731 cases showed that nulliparity, age at menopause, MHT, and BMI have independent, dose-response associations with ER, PR, and grade, clarifying patterns of etiologic heterogeneity. Associations with HER2, KI67 and TP53, or other risk factors did not meet our threshold for strong evidence.

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)

Abstract BackgroundCopy number variants (CNVs) are pervasive in the human genome but potential disease associations with rare CNVs have not been comprehensively assessed in large datasets. We analysed rare CNVs in genes and non-coding regions for 86,788 breast cancer cases and 76,122 controls of European ancestry with genome-wide array data.ResultsGene burden tests detected the strongest association for deletions in BRCA1 (P= 3.7E-18). Nine other genes were associated with a p-value < 0.01 including known susceptibility genes CHEK2 (P= 0.0008), ATM (P= 0.002) and BRCA2 (P= 0.008). Outside the known genes we detected associations with p-values < 0.001 for either overall or subtype-specific breast cancer at nine deletion regions and four duplication regions. Three of the deletion regions were in established common susceptibility loci.ConclusionsThis is the first genome-wide analysis of rare CNVs in a large breast cancer case-control dataset. We detected associations with exonic deletions in established breast cancer susceptibility genes. We also detected suggestive associations with non-coding CNVs in known and novel loci with large effects sizes. Larger sample sizes will be required to reach robust levels of statistical significance.

Abstract Background The purpose of this study was to estimate precise age-specific tubo-ovarian carcinoma (TOC) and breast cancer (BC) risks for carriers of pathogenic variants in RAD51C and RAD51D. Methods We analyzed data from 6178 families, 125 with pathogenic variants in RAD51C, and 6690 families, 60 with pathogenic variants in RAD51D. TOC and BC relative and cumulative risks were estimated using complex segregation analysis to model the cancer inheritance patterns in families while adjusting for the mode of ascertainment of each family. All statistical tests were two-sided. Results Pathogenic variants in both RAD51C and RAD51D were associated with TOC (RAD51C: relative risk [RR] = 7.55, 95% confidence interval [CI] = 5.60 to 10.19; P = 5 × 10-40; RAD51D: RR = 7.60, 95% CI = 5.61 to 10.30; P = 5 × 10-39) and BC (RAD51C: RR = 1.99, 95% CI = 1.39 to 2.85; P = 1.55 × 10-4; RAD51D: RR = 1.83, 95% CI = 1.24 to 2.72; P = .002). For both RAD51C and RAD51D, there was a suggestion that the TOC relative risks increased with age until around age 60 years and decreased thereafter. The estimated cumulative risks of developing TOC to age 80 years were 11% (95% CI = 6% to 21%) for RAD51C and 13% (95% CI = 7% to 23%) for RAD51D pathogenic variant carriers. The estimated cumulative risks of developing BC to 80 years were 21% (95% CI = 15% to 29%) for RAD51C and 20% (95% CI = 14% to 28%) for RAD51D pathogenic variant carriers. Both TOC and BC risks for RAD51C and RAD51D pathogenic variant carriers varied by cancer family history and could be as high as 32–36% for TOC, for carriers with two first-degree relatives diagnosed with TOC, or 44–46% for BC, for carriers with two first-degree relatives diagnosed with BC. Conclusions These estimates will facilitate the genetic counseling of RAD51C and RAD51D pathogenic variant carriers and justify the incorporation of RAD51C and RAD51D into cancer risk prediction models.

Abstract Pathogenic sequence variants (PSV) in BRCA1 or BRCA2 (BRCA1/2) are associated with increased risk and severity of prostate cancer. We evaluated whether PSVs in BRCA1/2 were associated with risk of overall prostate cancer or high grade (Gleason 8+) prostate cancer using an international sample of 65 BRCA1 and 171 BRCA2 male PSV carriers with prostate cancer, and 3,388 BRCA1 and 2,880 BRCA2 male PSV carriers without prostate cancer. PSVs in the 3′ region of BRCA2 (c.7914+) were significantly associated with elevated risk of prostate cancer compared with reference bin c.1001-c.7913 [HR = 1.78; 95% confidence interval (CI), 1.25–2.52; P = 0.001], as well as elevated risk of Gleason 8+ prostate cancer (HR = 3.11; 95% CI, 1.63–5.95; P = 0.001). c.756-c.1000 was also associated with elevated prostate cancer risk (HR = 2.83; 95% CI, 1.71–4.68; P = 0.00004) and elevated risk of Gleason 8+ prostate cancer (HR = 4.95; 95% CI, 2.12–11.54; P = 0.0002). No genotype–phenotype associations were detected for PSVs in BRCA1. These results demonstrate that specific BRCA2 PSVs may be associated with elevated risk of developing aggressive prostate cancer. Significance: Aggressive prostate cancer risk in BRCA2 mutation carriers may vary according to the specific BRCA2 mutation inherited by the at-risk individual.

Limited epidemiological evidence suggests that the etiology of hormone receptor positive (HR+) breast cancer may differ by levels of histologic grade and proliferation. We pooled risk factor and pathology data on 5,905 HR+ breast cancer cases and 26,281 controls from 11 epidemiological studies. Proliferation was determined by centralized automated measures of KI67 in tissue microarrays. Odds ratios (OR), 95% confidence intervals (CI) and p‐values for case–case and case–control comparisons for risk factors in relation to levels of grade and quartiles (Q1–Q4) of KI67 were estimated using polytomous logistic regression models. Case–case comparisons showed associations between nulliparity and high KI67 [OR (95% CI) for Q4 vs. Q1 = 1.54 (1.22, 1.95)]; obesity and high grade [grade 3 vs. 1 = 1.68 (1.31, 2.16)] and current use of combined hormone therapy (HT) and low grade [grade 3 vs. 1 = 0.27 (0.16, 0.44)] tumors. In case–control comparisons, nulliparity was associated with elevated risk of tumors with high but not low levels of proliferation [1.43 (1.14, 1.81) for KI67 Q4 vs. 0.83 (0.60, 1.14) for KI67 Q1]; obesity among women ≥50 years with high but not low grade tumors [1.55 (1.17, 2.06) for grade 3 vs. 0.88 (0.66, 1.16) for grade 1] and HT with low but not high grade tumors [3.07 (2.22, 4.23) for grade 1 vs. 0.85 (0.55, 1.30) for grade 3]. Menarcheal age and family history were similarly associated with HR+ tumors of different grade or KI67 levels. These findings provide insights into the etiologic heterogeneity of HR+ tumors.

Background The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study. Methods We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant. Results For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10−5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10−8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants. Conclusions This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.

AbstractAutomated methods are needed to facilitate high‐throughput and reproducible scoring of Ki67 and other markers in breast cancer tissue microarrays (TMAs) in large‐scale studies. To address this need, we developed an automated protocol for Ki67 scoring and evaluated its performance in studies from the Breast Cancer Association Consortium. We utilized 166 TMAs containing 16,953 tumour cores representing 9,059 breast cancer cases, from 13 studies, with information on other clinical and pathological characteristics. TMAs were stained for Ki67 using standard immunohistochemical procedures, and scanned and digitized using the Ariol system. An automated algorithm was developed for the scoring of Ki67, and scores were compared to computer assisted visual (CAV) scores in a subset of 15 TMAs in a training set. We also assessed the correlation between automated Ki67 scores and other clinical and pathological characteristics. Overall, we observed good discriminatory accuracy (AUC = 85%) and good agreement (kappa = 0.64) between the automated and CAV scoring methods in the training set. The performance of the automated method varied by TMA (kappa range= 0.37–0.87) and study (kappa range = 0.39–0.69). The automated method performed better in satisfactory cores (kappa = 0.68) than suboptimal (kappa = 0.51) cores (p‐value for comparison = 0.005); and among cores with higher total nuclei counted by the machine (4,000–4,500 cells: kappa = 0.78) than those with lower counts (50–500 cells: kappa = 0.41; p‐value = 0.010). Among the 9,059 cases in this study, the correlations between automated Ki67 and clinical and pathological characteristics were found to be in the expected directions. Our findings indicate that automated scoring of Ki67 can be an efficient method to obtain good quality data across large numbers of TMAs from multicentre studies. However, robust algorithm development and rigorous pre‐ and post‐analytical quality control procedures are necessary in order to ensure satisfactory performance.