Analysis of more than 400,000 women provides case-control evidence for BRCA1 and BRCA2 variant classification.

Clinical genetic testing identifies variants causal for hereditary cancer, information that is used for risk assessment and clinical management. Unfortunately, some variants identified are of uncertain clinical significance (VUS), complicating patient management. Case-control data is one evidence type used to classify VUS, and previous findings indicate that case-control likelihood ratios (LRs) outperform odds ratios for variant classification. As an initiative of the Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) Analytical Working Group we analyzed germline sequencing data of BRCA1 and BRCA2 from 96,691 female breast cancer cases and 303,925 unaffected controls from three studies: the BRIDGES study of the Breast Cancer Association Consortium, the Cancer Risk Estimates Related to Susceptibility consortium, and the UK Biobank. We observed 11,227 BRCA1 and BRCA2 variants, with 6,921 being coding, covering 23.4% of BRCA1 and BRCA2 VUS in ClinVar and 19.2% of ClinVar curated (likely) benign or pathogenic variants. Case-control LR evidence was highly consistent with ClinVar assertions for (likely) benign or pathogenic variants; exhibiting 99.1% sensitivity and 95.4% specificity for BRCA1 and 92.2% sensitivity and 86.6% specificity for BRCA2. This approach provides case-control evidence for 785 unclassified variants, that can serve as a valuable element for clinical classification.

Disentangling the relationships of body mass index and circulating sex hormone concentrations in mammographic density using Mendelian randomization.

PURPOSE: Mammographic density phenotypes, adjusted for age and body mass index (BMI), are strong predictors of breast cancer risk. BMI is associated with mammographic density measures, but the role of circulating sex hormone concentrations is less clear. We investigated the relationship between BMI, circulating sex hormone concentrations, and mammographic density phenotypes using Mendelian randomization (MR). METHODS: We applied two-sample MR approaches to assess the association between genetically predicted circulating concentrations of sex hormones [estradiol, testosterone, sex hormone-binding globulin (SHBG)], BMI, and mammographic density phenotypes (dense and non-dense area). We created instrumental variables from large European ancestry-based genome-wide association studies and applied estimates to mammographic density phenotypes in up to 14,000 women of European ancestry. We performed analyses overall and by menopausal status. RESULTS: Genetically predicted BMI was positively associated with non-dense area (IVW: ß = 1.79; 95% CI = 1.58, 2.00; p = 9.57 × 10-63) and inversely associated with dense area (IVW: ß = - 0.37; 95% CI = - 0.51,- 0.23; p = 4.7 × 10-7). We observed weak evidence for an association of circulating sex hormone concentrations with mammographic density phenotypes, specifically inverse associations between genetically predicted testosterone concentration and dense area (ß = - 0.22; 95% CI = - 0.38, - 0.053; p = 0.009) and between genetically predicted estradiol concentration and non-dense area (ß = - 3.32; 95% CI = - 5.83, - 0.82; p = 0.009), although results were not consistent across a range of MR approaches. CONCLUSION: Our findings support a positive causal association between BMI and mammographic non-dense area and an inverse association between BMI and dense area. Evidence was weaker and inconsistent for a causal effect of circulating sex hormone concentrations on mammographic density phenotypes. Based on our findings, associations between circulating sex hormone concentrations and mammographic density phenotypes are weak at best.

Differences in polygenic score distributions in European ancestry populations: implications for breast cancer risk prediction.

The 313-variant polygenic risk score (PRS313) provides a promising tool for breast cancer risk prediction. However, evaluation of the PRS313 across different European populations which could influence risk estimation has not been performed. Here, we explored the distribution of PRS313 across European populations using genotype data from 94,072 females without breast cancer, of European-ancestry from 21 countries participating in the Breast Cancer Association Consortium (BCAC) and 225,105 female participants from the UK Biobank. The mean PRS313 differed markedly across European countries, being highest in south-eastern Europe and lowest in north-western Europe. Using the overall European PRS313 distribution to categorise individuals leads to overestimation and underestimation of risk in some individuals from south-eastern and north-western countries, respectively. Adjustment for principal components explained most of the observed heterogeneity in mean PRS. Country-specific PRS distributions may be used to calibrate risk categories in individuals from different countries.

Age-specific breast and ovarian cancer risks associated with germline BRCA1 or BRCA2 pathogenic variants - an Asian study of 572 families.

Background: Clinical management of Asian BRCA1 and BRCA2 pathogenic variants (PV) carriers remains challenging due to imprecise age-specific breast (BC) and ovarian cancer (OC) risks estimates. We aimed to refine these estimates using six multi-ethnic studies in Asia. Methods: Data were collected on 271 BRCA1 and 301 BRCA2 families from Malaysia and Singapore, ascertained through population/hospital-based case-series (88%) and genetic clinics (12%). Age-specific cancer risks were estimated using a modified segregation analysis method, adjusted for ascertainment. Findings: BC and OC relative risks (RRs) varied across age groups for both BRCA1 and BRCA2. The age-specific RR estimates were similar across ethnicities and country of residence. For BRCA1 carriers of Malay, Indian and Chinese ancestry born between 1950 and 1959 in Malaysia, the cumulative risk (95% CI) of BC by age 80 was 40% (36%-44%), 49% (44%-53%) and 55% (51%-60%), respectively. The corresponding estimates for BRCA2 were 29% (26-32%), 36% (33%-40%) and 42% (38%-45%). The corresponding cumulative BC risks for Singapore residents from the same birth cohort, where the underlying population cancer incidences are higher compared to Malaysia, were higher, varying by ancestry group between 57 and 61% for BRCA1, and between 43 and 47% for BRCA2 carriers. The cumulative risk of OC by age 80 was 31% (27-36%) for BRCA1 and 12% (10%-15%) for BRCA2 carriers in Malaysia born between 1950 and 1959; and 42% (34-50%) for BRCA1 and 20% (14-27%) for BRCA2 carriers of the same birth cohort in Singapore. There was evidence of increased BC and OC risks for women from >1960 birth cohorts (p-value = 3.6 × 10-5 for BRCA1 and 0.018 for BRCA2). Interpretation: The absolute age-specific cancer risks of Asian carriers vary depending on the underlying population-specific cancer incidences, and hence should be customised to allow for more accurate cancer risk management. Funding: Wellcome Trust [grant no: v203477/Z/16/Z]; CRUK (PPRPGM-Nov20∖100002).

Novel breast cancer susceptibility loci under linkage peaks identified in African ancestry consortia.

BACKGROUND: Expansion of genome-wide association studies across population groups is needed to improve our understanding of shared and unique genetic contributions to breast cancer. We performed association and replication studies guided by a priori linkage findings from African ancestry (AA) relative pairs. METHODS: We performed fixed-effect inverse-variance weighted meta-analysis under three significant AA breast cancer linkage peaks (3q26-27, 12q22-23, and 16q21-22) in 9241 AA cases and 10 193 AA controls. We examined associations with overall breast cancer as well as estrogen receptor (ER)-positive and negative subtypes (193,132 SNPs). We replicated associations in the African-ancestry Breast Cancer Genetic Consortium (AABCG). RESULTS: In AA women, we identified two associations on chr12q for overall breast cancer (rs1420647, OR = 1.15, p = 2.50×10-6; rs12322371, OR = 1.14, p = 3.15×10-6), and one for ER-negative breast cancer (rs77006600, OR = 1.67, p = 3.51×10-6). On chr3, we identified two associations with ER-negative disease (rs184090918, OR = 3.70, p = 1.23×10-5; rs76959804, OR = 3.57, p = 1.77×10-5) and on chr16q we identified an association with ER-negative disease (rs34147411, OR = 1.62, p = 8.82×10-6). In the replication study, the chr3 associations were significant and effect sizes were larger (rs184090918, OR: 6.66, 95% CI: 1.43, 31.01; rs76959804, OR: 5.24, 95% CI: 1.70, 16.16). CONCLUSION: The two chr3 SNPs are upstream to open chromatin ENSR00000710716, a regulatory feature that is actively regulated in mammary tissues, providing evidence that variants in this chr3 region may have a regulatory role in our target organ. Our study provides support for breast cancer variant discovery using prioritization based on linkage evidence.

Polymorphisms in genes of melatonin biosynthesis and signaling support the light-at-night hypothesis for breast cancer.

Light-at-night triggers the decline of pineal gland melatonin biosynthesis and secretion and is an IARC-classified probable breast-cancer risk factor. We applied a large-scale molecular epidemiology approach to shed light on the putative role of melatonin in breast cancer. We investigated associations between breast-cancer risk and polymorphisms at genes of melatonin biosynthesis/signaling using a study population of 44,405 women from the Breast Cancer Association Consortium (22,992 cases, 21,413 population-based controls). Genotype data of 97 candidate single nucleotide polymorphisms (SNPs) at 18 defined gene regions were investigated for breast-cancer risk effects. We calculated adjusted odds ratios (ORs) and 95% confidence intervals (CI) by logistic regression for the main-effect analysis as well as stratified analyses by estrogen- and progesterone-receptor (ER, PR) status. SNP-SNP interactions were analyzed via a two-step procedure based on logic regression. The Bayesian false-discovery probability (BFDP) was used for all analyses to account for multiple testing. Noteworthy associations (BFDP < 0.8) included 10 linked SNPs in tryptophan hydroxylase 2 (TPH2) (e.g. rs1386492: OR = 1.07, 95% CI 1.02-1.12), and a SNP in the mitogen-activated protein kinase 8 (MAPK8) (rs10857561: OR = 1.11, 95% CI 1.04-1.18). The SNP-SNP interaction analysis revealed noteworthy interaction terms with TPH2- and MAPK-related SNPs (e.g. rs1386483R ∧ rs1473473D ∧ rs3729931D: OR = 1.20, 95% CI 1.09-1.32). In line with the light-at-night hypothesis that links shift work with elevated breast-cancer risks our results point to SNPs in TPH2 and MAPK-genes that may impact the intricate network of circadian regulation.

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival.

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.

Genome-wide association study of treatment-related toxicity two years following radiotherapy for breast cancer.

BACKGROUND AND PURPOSE: Up to a quarter of breast cancer patients treated by surgery and radiotherapy experience clinically significant toxicity. If patients at high risk of adverse effects could be identified at diagnosis, their treatment could be tailored accordingly. This study was designed to identify common single nucleotide polymorphisms (SNPs) associated with toxicity two years following whole breast radiotherapy. MATERIALS AND METHODS: A genome-wide association study (GWAS) was performed in 1,640 breast cancer patients with complete SNP, clinical, treatment and toxicity data, recruited across 18 European and US centres into the prospective REQUITE cohort study. Toxicity data (CTCAE v4.0) were collected at baseline, end of radiotherapy, and annual follow-up. A total of 7,097,340 SNPs were tested for association with the residuals of toxicity endpoints, adjusted for clinical, treatment co-variates and population substructure. RESULTS: Quantile-quantile plots showed more associations with toxicity above the p < 5 × 10-5 level than expected by chance. Eight SNPs reached genome-wide significance. Nipple retraction grade ≥ 2 was associated with the rs188287402 variant (p = 2.80 × 10-8), breast oedema grade ≥ 2 with rs12657177 (p = 1.12 × 10-10), rs75912034 (p = 1.12 × 10-10), rs145328458 (p = 1.06 × 10-9) and rs61966612 (p = 1.23 × 10-9), induration grade ≥ 2 with rs77311050 (p = 2.54 × 10-8) and rs34063419 (p = 1.21 × 10-8), and arm lymphoedema grade ≥ 1 with rs643644 (p = 3.54 × 10-8). Heritability estimates across significant endpoints ranged from 25% to 39%. Our study did not replicate previously reported SNPs associated with breast radiation toxicity at the pre-specified significance level. CONCLUSIONS: This GWAS for long-term breast radiation toxicity provides further evidence for significant association of common SNPs with distinct toxicity endpoints.

PREDICT validity for prognosis of breast cancer patients with pathogenic BRCA1/2 variants.

We assessed the PREDICT v 2.2 for prognosis of breast cancer patients with pathogenic germline BRCA1 and BRCA2 variants, using follow-up data from 5453 BRCA1/2 carriers from the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) and the Breast Cancer Association Consortium (BCAC). PREDICT for estrogen receptor (ER)-negative breast cancer had modest discrimination for BRCA1 carrier patients overall (Gönen & Heller unbiased concordance 0.65 in CIMBA, 0.64 in BCAC), but it distinguished clearly the high-mortality group from lower risk categories. In an analysis of low to high risk categories by PREDICT score percentiles, the observed mortality was consistently lower than the expected mortality, but the confidence intervals always included the calibration slope. Altogether, our results encourage the use of the PREDICT ER-negative model in management of breast cancer patients with germline BRCA1 variants. For the PREDICT ER-positive model, the discrimination was slightly lower in BRCA2 variant carriers (concordance 0.60 in CIMBA, 0.65 in BCAC). Especially, inclusion of the tumor grade distorted the prognostic estimates. The breast cancer mortality of BRCA2 carriers was underestimated at the low end of the PREDICT score distribution, whereas at the high end, the mortality was overestimated. These data suggest that BRCA2 status should also be taken into consideration with tumor characteristics, when estimating the prognosis of ER-positive breast cancer patients.

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival.

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

The impact of coding germline variants on contralateral breast cancer risk and survival.

Evidence linking coding germline variants in breast cancer (BC)-susceptibility genes other than BRCA1, BRCA2, and CHEK2 with contralateral breast cancer (CBC) risk and breast cancer-specific survival (BCSS) is scarce. The aim of this study was to assess the association of protein-truncating variants (PTVs) and rare missense variants (MSVs) in nine known (ATM, BARD1, BRCA1, BRCA2, CHEK2, PALB2, RAD51C, RAD51D, and TP53) and 25 suspected BC-susceptibility genes with CBC risk and BCSS. Hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated with Cox regression models. Analyses included 34,401 women of European ancestry diagnosed with BC, including 676 CBCs and 3,449 BC deaths; the median follow-up was 10.9 years. Subtype analyses were based on estrogen receptor (ER) status of the first BC. Combined PTVs and pathogenic/likely pathogenic MSVs in BRCA1, BRCA2, and TP53 and PTVs in CHEK2 and PALB2 were associated with increased CBC risk [HRs (95% CIs): 2.88 (1.70-4.87), 2.31 (1.39-3.85), 8.29 (2.53-27.21), 2.25 (1.55-3.27), and 2.67 (1.33-5.35), respectively]. The strongest evidence of association with BCSS was for PTVs and pathogenic/likely pathogenic MSVs in BRCA2 (ER-positive BC) and TP53 and PTVs in CHEK2 [HRs (95% CIs): 1.53 (1.13-2.07), 2.08 (0.95-4.57), and 1.39 (1.13-1.72), respectively, after adjusting for tumor characteristics and treatment]. HRs were essentially unchanged when censoring for CBC, suggesting that these associations are not completely explained by increased CBC risk, tumor characteristics, or treatment. There was limited evidence of associations of PTVs and/or rare MSVs with CBC risk or BCSS for the 25 suspected BC genes. The CBC findings are relevant to treatment decisions, follow-up, and screening after BC diagnosis.

A likelihood ratio approach for utilizing case-control data in the clinical classification of rare sequence variants: application to BRCA1 and BRCA2.

A large number of variants identified through clinical genetic testing in disease susceptibility genes, are of uncertain significance (VUS). Following the recommendations of the American College of Medical Genetics and Genomics (ACMG) and Association for Molecular Pathology (AMP), the frequency in case-control datasets (PS4 criterion), can inform their interpretation. We present a novel case-control likelihood ratio-based method that incorporates gene-specific age-related penetrance. We demonstrate the utility of this method in the analysis of simulated and real datasets. In the analyses of simulated data, the likelihood ratio method was more powerful compared to other methods. Likelihood ratios were calculated for a case-control dataset of BRCA1 and BRCA2 variants from the Breast Cancer Association Consortium (BCAC), and compared with logistic regression results. A larger number of variants reached evidence in favor of pathogenicity, and a substantial number of variants had evidence against pathogenicity - findings that would not have been reached using other case-control analysis methods. Our novel method provides greater power to classify rare variants compared to classical case-control methods. As an initiative from the ENIGMA Analytical Working Group, we provide user-friendly scripts and pre-formatted excel calculators for implementation of the method for rare variants in BRCA1, BRCA2 and other high-risk genes with known penetrance.

Genome- and transcriptome-wide association studies of 386,000 Asian and European-ancestry women provide new insights into breast cancer genetics.

By combining data from 160,500 individuals with breast cancer and 226,196 controls of Asian and European ancestry, we conducted genome- and transcriptome-wide association studies of breast cancer. We identified 222 genetic risk loci and 137 genes that were associated with breast cancer risk at a p < 5.0 × 10-8 and a Bonferroni-corrected p < 4.6 × 10-6, respectively. Of them, 32 loci and 15 genes showed a significantly different association between ER-positive and ER-negative breast cancer after Bonferroni correction. Significant ancestral differences in risk variant allele frequencies and their association strengths with breast cancer risk were identified. Of the significant associations identified in this study, 17 loci and 14 genes are located 1Mb away from any of the previously reported breast cancer risk variants. Pathways analyses including 221 putative risk genes identified multiple signaling pathways that may play a significant role in the development of breast cancer. Our study provides a comprehensive understanding of and new biological insights into the genetics of this common malignancy.

Segregation analysis of 17,425 population-based breast cancer families: Evidence for genetic susceptibility and risk prediction.

Rare pathogenic variants in known breast cancer-susceptibility genes and known common susceptibility variants do not fully explain the familial aggregation of breast cancer. To investigate plausible genetic models for the residual familial aggregation, we studied 17,425 families ascertained through population-based probands, 86% of whom were screened for pathogenic variants in BRCA1, BRCA2, PALB2, CHEK2, ATM, and TP53 via gene-panel sequencing. We conducted complex segregation analyses and fitted genetic models in which breast cancer incidence depended on the effects of known susceptibility genes and other unidentified major genes and a normally distributed polygenic component. The proportion of familial variance explained by the six genes was 46% at age 20-29 years and decreased steadily with age thereafter. After allowing for these genes, the best fitting model for the residual familial variance included a recessive risk component with a combined genotype frequency of 1.7% (95% CI: 0.3%-5.4%) and a penetrance to age 80 years of 69% (95% CI: 38%-95%) for homozygotes, which may reflect the combined effects of multiple variants acting in a recessive manner, and a polygenic variance of 1.27 (95% CI: 0.94%-1.65), which did not vary with age. The proportion of the residual familial variance explained by the recessive risk component was 40% at age 20-29 years and decreased with age thereafter. The model predicted age-specific familial relative risks consistent with those observed by large epidemiological studies. The findings have implications for strategies to identify new breast cancer-susceptibility genes and improve disease-risk prediction, especially at a young age.

A genome-wide gene-based gene-environment interaction study of breast cancer in more than 90,000 women.

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 (PGXE=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.

Treatment time and circadian genotype interact to influence radiotherapy side-effects. A prospective European validation study using the REQUITE cohort.

BACKGROUND: Circadian rhythm impacts broad biological processes, including response to cancer treatment. Evidence conflicts on whether treatment time affects risk of radiotherapy side-effects, likely because of differing time analyses and target tissues. We previously showed interactive effects of time and genotypes of circadian genes on late toxicity after breast radiotherapy and aimed to validate those results in a multi-centre cohort. METHODS: Clinical and genotype data from 1690 REQUITE breast cancer patients were used with erythema (acute; n=340) and breast atrophy (two years post-radiotherapy; n=514) as primary endpoints. Local datetimes per fraction were converted into solar times as predictors. Genetic chronotype markers were included in logistic regressions to identify primary endpoint predictors. FINDINGS: Significant predictors for erythema included BMI, radiation dose and PER3 genotype (OR 1.27(95%CI 1.03-1.56); P < 0.03). Effect of treatment time effect on acute toxicity was inconclusive, with no interaction between time and genotype. For late toxicity (breast atrophy), predictors included BMI, radiation dose, surgery type, treatment time and SNPs in CLOCK (OR 0.62 (95%CI 0.4-0.9); P < 0.01), PER3 (OR 0.65 (95%CI 0.44-0.97); P < 0.04) and RASD1 (OR 0.56 (95%CI 0.35-0.89); P < 0.02). There was a statistically significant interaction between time and genotypes of circadian rhythm genes (CLOCK OR 1.13 (95%CI 1.03-1.23), P < 0.01; PER3 OR 1.1 (95%CI 1.01-1.2), P < 0.04; RASD1 OR 1.15 (95%CI 1.04-1.28), P < 0.008), with peak time for toxicity determined by genotype. INTERPRETATION: Late atrophy can be mitigated by selecting optimal treatment time according to circadian genotypes (e.g. treat PER3 rs2087947C/C genotypes in mornings; T/T in afternoons). We predict triple-homozygous patients (14%) reduce chance of atrophy from 70% to 33% by treating in mornings as opposed to mid-afternoon. Future clinical trials could stratify patients treated at optimal times compared to those scheduled normally. FUNDING: EU-FP7.

Incorporating progesterone receptor expression into the PREDICT breast prognostic model.

BACKGROUND: Predict Breast (www.predict.nhs.uk) is an online prognostication and treatment benefit tool for early invasive breast cancer. The aim of this study was to incorporate the prognostic effect of progesterone receptor (PR) status into a new version of PREDICT and to compare its performance to the current version (2.2). METHOD: The prognostic effect of PR status was based on the analysis of data from 45,088 European patients with breast cancer from 49 studies in the Breast Cancer Association Consortium. Cox proportional hazard models were used to estimate the hazard ratio for PR status. Data from a New Zealand study of 11,365 patients with early invasive breast cancer were used for external validation. Model calibration and discrimination were used to test the model performance. RESULTS: Having a PR-positive tumour was associated with a 23% and 28% lower risk of dying from breast cancer for women with oestrogen receptor (ER)-negative and ER-positive breast cancer, respectively. The area under the ROC curve increased with the addition of PR status from 0.807 to 0.809 for patients with ER-negative tumours (p = 0.023) and from 0.898 to 0.902 for patients with ER-positive tumours (p = 2.3 × 10-6) in the New Zealand cohort. Model calibration was modest with 940 observed deaths compared to 1151 predicted. CONCLUSION: The inclusion of the prognostic effect of PR status to PREDICT Breast has led to an improvement of model performance and more accurate absolute treatment benefit predictions for individual patients. Further studies should determine whether the baseline hazard function requires recalibration.

Uncovering the Contribution of Moderate-Penetrance Susceptibility Genes to Breast Cancer by Whole-Exome Sequencing and Targeted Enrichment Sequencing of Candidate Genes in Women of European Ancestry.

Rare variants in at least 10 genes, including BRCA1, BRCA2, PALB2, ATM, and CHEK2, are associated with increased risk of breast cancer; however, these variants, in combination with common variants identified through genome-wide association studies, explain only a fraction of the familial aggregation of the disease. To identify further susceptibility genes, we performed a two-stage whole-exome sequencing study. In the discovery stage, samples from 1528 breast cancer cases enriched for breast cancer susceptibility and 3733 geographically matched unaffected controls were sequenced. Using five different filtering and gene prioritization strategies, 198 genes were selected for further validation. These genes, and a panel of 32 known or suspected breast cancer susceptibility genes, were assessed in a validation set of 6211 cases and 6019 controls for their association with risk of breast cancer overall, and by estrogen receptor (ER) disease subtypes, using gene burden tests applied to loss-of-function and rare missense variants. Twenty genes showed nominal evidence of association (p-value < 0.05) with either overall or subtype-specific breast cancer. Our study had the statistical power to detect susceptibility genes with effect sizes similar to ATM, CHEK2, and PALB2, however, it was underpowered to identify genes in which susceptibility variants are rarer or confer smaller effect sizes. Larger sample sizes would be required in order to identify such genes.

No Association Between Polygenic Risk Scores for Cancer and Development of Radiation Therapy Toxicity.

PURPOSE: Our aim was to test whether updated polygenic risk scores (PRS) for susceptibility to cancer affect risk of radiation therapy toxicity. METHODS AND MATERIALS: Analyses included 9,717 patients with breast (n=3,078), prostate (n=5,748) or lung (n=891) cancer from Radiogenomics and REQUITE Consortia cohorts. Patients underwent potentially curative radiation therapy and were assessed prospectively for toxicity. Germline genotyping involved genome-wide single nucleotide polymorphism (SNP) arrays with nontyped SNPs imputed. PRS for each cancer were generated by summing literature-identified cancer susceptibility risk alleles: 352 breast, 136 prostate, and 24 lung. Weighted PRS were generated using log odds ratio (ORs) for cancer susceptibility. Standardized total average toxicity (STAT) scores at 2 and 5 years (breast, prostate) or 6 to 12 months (lung) quantified toxicity. Primary analysis tested late STAT, secondary analyses investigated acute STAT, and individual endpoints and SNPs using multivariable regression. RESULTS: Increasing PRS did not increase risk of late toxicity in patients with breast (OR, 1.000; 95% confidence interval [CI], 0.997-1.002), prostate (OR, 0.99; 95% CI, 0.98-1.00; weighted PRS OR, 0.93; 95% CI, 0.83-1.03), or lung (OR, 0.93; 95% CI, 0.87-1.00; weighted PRS OR, 0.68; 95% CI, 0.45-1.03) cancer. Similar results were seen for acute toxicity. Secondary analyses identified rs138944387 associated with breast pain (OR, 3.05; 95% CI, 1.86-5.01; P = 1.09 × 10-5) and rs17513613 with breast edema (OR, 0.94; 95% CI, 0.92-0.97; P = 1.08 × 10-5). CONCLUSIONS: Patients with increased polygenic predisposition to breast, prostate, or lung cancer can safely undergo radiation therapy with no anticipated excess toxicity risk. Some individual SNPs increase the likelihood of a specific toxicity endpoint, warranting validation in independent cohorts and functional studies to elucidate biologic mechanisms.

Development and Optimization of a Machine-Learning Prediction Model for Acute Desquamation After Breast Radiation Therapy in the Multicenter REQUITE Cohort.

Purpose: Some patients with breast cancer treated by surgery and radiation therapy experience clinically significant toxicity, which may adversely affect cosmesis and quality of life. There is a paucity of validated clinical prediction models for radiation toxicity. We used machine learning (ML) algorithms to develop and optimise a clinical prediction model for acute breast desquamation after whole breast external beam radiation therapy in the prospective multicenter REQUITE cohort study. Methods and Materials: Using demographic and treatment-related features (m = 122) from patients (n = 2058) at 26 centers, we trained 8 ML algorithms with 10-fold cross-validation in a 50:50 random-split data set with class stratification to predict acute breast desquamation. Based on performance in the validation data set, the logistic model tree, random forest, and naïve Bayes models were taken forward to cost-sensitive learning optimisation. Results: One hundred and ninety-two patients experienced acute desquamation. Resampling and cost-sensitive learning optimisation facilitated an improvement in classification performance. Based on maximising sensitivity (true positives), the "hero" model was the cost-sensitive random forest algorithm with a false-negative: false-positive misclassification penalty of 90:1 containing m = 114 predictive features. Model sensitivity and specificity were 0.77 and 0.66, respectively, with an area under the curve of 0.77 in the validation cohort. Conclusions: ML algorithms with resampling and cost-sensitive learning generated clinically valid prediction models for acute desquamation using patient demographic and treatment features. Further external validation and inclusion of genomic markers in ML prediction models are worthwhile, to identify patients at increased risk of toxicity who may benefit from supportive intervention or even a change in treatment plan.