Citations

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)

Abstract Background Obesity is an established risk factor for colorectal cancer (CRC), but the evidence for the association is inconsistent across molecular subtypes of the disease. Methods We pooled data on body mass index (BMI), tumor microsatellite instability status, CpG island methylator phenotype status, BRAF and KRAS mutations, and Jass classification types for 11 872 CRC cases and 11 013 controls from 11 observational studies. We used multinomial logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CI) adjusted for covariables. Results Higher BMI was associated with increased CRC risk (OR per 5 kg/m2 = 1.18, 95% CI = 1.15 to 1.22). The positive association was stronger for men than women but similar across tumor subtypes defined by individual molecular markers. In analyses by Jass type, higher BMI was associated with elevated CRC risk for types 1-4 cases but not for type 5 CRC cases (considered familial-like/Lynch syndrome microsatellite instability-H, CpG island methylator phenotype-low or negative, BRAF-wild type, KRAS-wild type, OR = 1.04, 95% CI = 0.90 to 1.20). This pattern of associations for BMI and Jass types was consistent by sex and design of contributing studies (cohort or case-control). Conclusions In contrast to previous reports with fewer study participants, we found limited evidence of heterogeneity for the association between BMI and CRC risk according to molecular subtype, suggesting that obesity influences nearly all major pathways involved in colorectal carcinogenesis. The null association observed for the Jass type 5 suggests that BMI is not a risk factor for the development of CRC for individuals with Lynch syndrome.

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.

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. 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. 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 (PGXE = 4.44 × 10−6). 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. Our study suggests a limited role of gene–environment interactions in breast cancer risk.

Abstract Acquired PARP inhibitor (PARPi) resistance in BRCA1- or BRCA2-mutant ovarian cancer often results from secondary mutations that restore expression of functional protein. RAD51C is a less commonly studied ovarian cancer susceptibility gene whose promoter is sometimes methylated, leading to homologous recombination (HR) deficiency and PARPi sensitivity. For this study, the PARPi-sensitive patient-derived ovarian cancer xenograft PH039, which lacks HR gene mutations but harbors RAD51C promoter methylation, was selected for PARPi resistance by cyclical niraparib treatment in vivo. PH039 acquired PARPi resistance by the third treatment cycle and grew through subsequent treatment with either niraparib or rucaparib. Transcriptional profiling throughout the course of resistance development showed widespread pathway level changes along with a marked increase in RAD51C mRNA, which reflected loss of RAD51C promoter methylation. Analysis of ovarian cancer samples from the ARIEL2 Part 1 clinical trial of rucaparib monotherapy likewise indicated an association between loss of RAD51C methylation prior to on-study biopsy and limited response. Interestingly, the PARPi resistant PH039 model remained platinum sensitive. Collectively, these results not only indicate that PARPi treatment pressure can reverse RAD51C methylation and restore RAD51C expression, but also provide a model for studying the clinical observation that PARPi and platinum sensitivity are sometimes dissociated.

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 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: 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.

In this study we aim to examine gene–environment interactions (GxEs) between genes involved with estrogen metabolism and environmental factors related to estrogen exposure. GxE analyses were conducted with 1970 Korean breast cancer cases and 2052 controls in the case-control study, the Seoul Breast Cancer Study (SEBCS). A total of 11,555 SNPs from the 137 candidate genes were included in the GxE analyses with eight established environmental factors. A replication test was conducted by using an independent population from the Breast Cancer Association Consortium (BCAC), with 62,485 Europeans and 9047 Asians. The GxE tests were performed by using two-step methods in GxEScan software. Two interactions were found in the SEBCS. The first interaction was shown between rs13035764 of NCOA1 and age at menarche in the GE|2df model (p-2df = 1.2 × 10−3). The age at menarche before 14 years old was associated with the high risk of breast cancer, and the risk was higher when subjects had homozygous minor allele G. The second GxE was shown between rs851998 near ESR1 and height in the GE|2df model (p-2df = 1.1 × 10−4). Height taller than 160 cm was associated with a high risk of breast cancer, and the risk increased when the minor allele was added. The findings were not replicated in the BCAC. These results would suggest specificity in Koreans for breast cancer risk.