Integrative multi-omics analyses to identify the genetic and functional mechanisms underlying ovarian cancer risk regions.

To identify credible causal risk variants (CCVs) associated with different histotypes of epithelial ovarian cancer (EOC), we performed genome-wide association analysis for 470,825 genotyped and 10,163,797 imputed SNPs in 25,981 EOC cases and 105,724 controls of European origin. We identified five histotype-specific EOC risk regions (p value <5 × 10-8) and confirmed previously reported associations for 27 risk regions. Conditional analyses identified an additional 11 signals independent of the primary signal at six risk regions (p value <10-5). Fine mapping identified 4,008 CCVs in these regions, of which 1,452 CCVs were located in ovarian cancer-related chromatin marks with significant enrichment in active enhancers, active promoters, and active regions for CCVs from each EOC histotype. Transcriptome-wide association and colocalization analyses across histotypes using tissue-specific and cross-tissue datasets identified 86 candidate susceptibility genes in known EOC risk regions and 32 genes in 23 additional genomic regions that may represent novel EOC risk loci (false discovery rate <0.05). Finally, by integrating genome-wide HiChIP interactome analysis with transcriptome-wide association study (TWAS), variant effect predictor, transcription factor ChIP-seq, and motifbreakR data, we identified candidate gene-CCV interactions at each locus. This included risk loci where TWAS identified one or more candidate susceptibility genes (e.g., HOXD-AS2, HOXD8, and HOXD3 at 2q31) and other loci where no candidate gene was identified (e.g., MYC and PVT1 at 8q24) by TWAS. In summary, this study describes a functional framework and provides a greater understanding of the biological significance of risk alleles and candidate gene targets at EOC susceptibility loci identified by a genome-wide association study.

Exome sequencing identifies HELB as a novel susceptibility gene for non-mucinous, non-high-grade-serous epithelial ovarian cancer.

Rare, germline loss-of-function variants in a handful of genes that encode DNA repair proteins have been shown to be associated with epithelial ovarian cancer with a stronger association for the high-grade serous hiostotype. The aim of this study was to collate exome sequencing data from multiple epithelial ovarian cancer case cohorts and controls in order to systematically evaluate the role of coding, loss-of-function variants across the genome in epithelial ovarian cancer risk. We assembled exome data for a total of 2,573 non-mucinous cases (1,876 high-grade serous and 697 non-high grade serous) and 13,925 controls. Harmonised variant calling and quality control filtering was applied across the different data sets. We carried out a gene-by-gene simple burden test for association of rare loss-of-function variants (minor allele frequency < 0.1%) with all non-mucinous ovarian cancer, high grade serous ovarian cancer and non-high grade serous ovarian cancer using logistic regression adjusted for the top four principal components to account for cryptic population structure and genetic ancestry. Seven of the top 10 associated genes were associations of the known ovarian cancer susceptibility genes BRCA1, BRCA2, BRIP1, RAD51C, RAD51D, MSH6 and PALB2 (false discovery probability < 0.1). A further four genes (HELB, OR2T35, NBN and MYO1A) had a false discovery rate of less than 0.1. Of these, HELB was most strongly associated with the non-high grade serous histotype (P = 1.3×10-6, FDR = 9.1×10-4). Further support for this association comes from the observation that loss of function variants in this gene are also associated with age at natural menopause and Mendelian randomisation analysis shows an association between genetically predicted age at natural menopause and endometrioid ovarian cancer, but not high-grade serous ovarian cancer.

Large-scale genome-wide association study of 398,238 women unveils seven novel loci associated with high-grade serous epithelial ovarian cancer risk.

Background: Nineteen genomic regions have been associated with high-grade serous ovarian cancer (HGSOC). We used data from the Ovarian Cancer Association Consortium (OCAC), Consortium of Investigators of Modifiers of BRCA1/BRCA2 (CIMBA), UK Biobank (UKBB), and FinnGen to identify novel HGSOC susceptibility loci and develop polygenic scores (PGS). Methods: We analyzed >22 million variants for 398,238 women. Associations were assessed separately by consortium and meta-analysed. OCAC and CIMBA data were used to develop PGS which were trained on FinnGen data and validated in UKBB and BioBank Japan. Results: Eight novel variants were associated with HGSOC risk. An interesting discovery biologically was finding that TP53 3'-UTR SNP rs78378222 was associated with HGSOC (per T allele relative risk (RR)=1.44, 95%CI:1.28-1.62, P=1.76×10-9). The optimal PGS included 64,518 variants and was associated with an odds ratio of 1.46 (95%CI:1.37-1.54) per standard deviation in the UKBB validation (AUROC curve=0.61, 95%CI:0.59-0.62). Conclusions: This study represents the largest GWAS for HGSOC to date. The results highlight that improvements in imputation reference panels and increased sample sizes can identify HGSOC associated variants that previously went undetected, resulting in improved PGS. The use of updated PGS in cancer risk prediction algorithms will then improve personalized risk prediction for HGSOC.

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.

Revisiting the association of sudden infant death syndrome (SIDS) with polymorphisms of NHE3 and IL13.

OBJECTIVES: Disturbances of the central nervous system and immune system are thought to play a role in sudden infant death syndrome (SIDS). Dysregulated expression of sodium (Na+)/hydrogen (H+) exchanger 3 (NHE3) in the brainstem and of interleukin 13 (IL13) in the lungs has been observed in SIDS. An association of single-nucleotide polymorphisms (SNPs) in NHE3 and IL13 with SIDS has been proposed, but controversial results were reported. Therefore, there is a need to revisit the association of SNPs in NHE3 and IL13 with SIDS. METHODS: Genotyping of rs71597645 (G1131A) and rs2247114 (C2405T) in NHE3 and rs20541 (+ 4464A/G) in IL13 was performed in 201 SIDS cases and 338 controls. A meta-analysis was performed after merging our data with previously published data (all from European populations). RESULTS: Polymorphisms rs2247114 (NHE3) and rs20541 (IL13) were significantly associated with SIDS overall and in multiple subgroups, but no association was found for rs71597645 (NHE3). After combining our data with previously published data, a fixed-effect meta-analysis showed that rs2247114 in NHE3 retained a significant association with SIDS under a recessive model (OR 2.78, 95%CI 1.53 to 5.06; p = 0.0008). CONCLUSION: Our findings suggest an association of NHE3 variant rs2247114 (C2405T), though not rs71597645 (NHE3), with SIDS. A potential role of rs20541 (IL13) still has to be elucidated. Especially NHE3 seems to be an interesting topic for future SIDS research.

Bone Marrow Failure and Immunodeficiency Associated with Human RAD50 Variants.

PURPOSE: The MRE11-RAD50-NBN (MRN) complex plays a key role in recognizing and signaling DNA double-strand breaks. Pathogenic variants in NBN and MRE11 give rise to the autosomal-recessive diseases, Nijmegen breakage syndrome (NBS) and ataxia telangiectasia-like disorder, respectively. The clinical consequences of pathogenic variants in RAD50 are incompletely understood. We aimed to characterize a newly identified RAD50 deficiency/NBS-like disorder (NBSLD) patient with bone marrow failure and immunodeficiency. METHODS: We report on a girl with microcephaly, mental retardation, bird-like face, short stature, bone marrow failure and B-cell immunodeficiency. We searched for candidate gene by whole-exome sequencing and analyzed the cellular phenotype of patient-derived fibroblasts using immunoblotting, radiation sensitivity assays and lentiviral complementation experiments. RESULTS: Compound heterozygosity for two variants in the RAD50 gene (p.Arg83His and p.Glu485Ter) was identified in this patient. The expression of RAD50 protein and MRN complex formation was maintained in the cells derived from this patient. DNA damage-induced activation of the ATM kinase was markedly decreased, which was restored by the expression of wild-type (WT) RAD50. Radiosensitivity appeared inconspicuous in the patient-derived cell line as assessed by colony formation assay. The RAD50R83H missense substitution did not rescue the mitotic defect in complementation experiments using RAD50-deficient fibroblasts, whereas RAD50WT did. The RAD50E485X nonsense variant was associated with in-frame skipping of exon 10 (p.Glu485_545del). CONCLUSION: These findings indicate important roles of RAD50 in human bone marrow and immune cells. RAD50 deficiency/NBSLD can manifest as a distinct inborn error of immunity characterized by bone marrow failure and B-cell immunodeficiency.

Exome sequencing identifies breast cancer susceptibility genes and defines the contribution of coding variants to breast cancer risk.

Linkage and candidate gene studies have identified several breast cancer susceptibility genes, but the overall contribution of coding variation to breast cancer is unclear. To evaluate the role of rare coding variants more comprehensively, we performed a meta-analysis across three large whole-exome sequencing datasets, containing 26,368 female cases and 217,673 female controls. Burden tests were performed for protein-truncating and rare missense variants in 15,616 and 18,601 genes, respectively. Associations between protein-truncating variants and breast cancer were identified for the following six genes at exome-wide significance (P < 2.5 × 10-6): the five known susceptibility genes ATM, BRCA1, BRCA2, CHEK2 and PALB2, together with MAP3K1. Associations were also observed for LZTR1, ATR and BARD1 with P < 1 × 10-4. Associations between predicted deleterious rare missense or protein-truncating variants and breast cancer were additionally identified for CDKN2A at exome-wide significance. The overall contribution of coding variants in genes beyond the previously known genes is estimated to be small.

Spectrum and Frequency of Germline FANCM Protein-Truncating Variants in 44,803 European Female Breast Cancer Cases.

FANCM germline protein truncating variants (PTVs) are moderate-risk factors for ER-negative breast cancer. We previously described the spectrum of FANCM PTVs in 114 European breast cancer cases. In the present, larger cohort, we report the spectrum and frequency of four common and 62 rare FANCM PTVs found in 274 carriers detected among 44,803 breast cancer cases. We confirmed that p.Gln1701* was the most common PTV in Northern Europe with lower frequencies in Southern Europe. In contrast, p.Gly1906Alafs*12 was the most common PTV in Southern Europe with decreasing frequencies in Central and Northern Europe. We verified that p.Arg658* was prevalent in Central Europe and had highest frequencies in Eastern Europe. We also confirmed that the fourth most common PTV, p.Gln498Thrfs*7, might be a founder variant from Lithuania. Based on the frequency distribution of the carriers of rare PTVs, we showed that the FANCM PTVs spectra in Southwestern and Central Europe were much more heterogeneous than those from Northeastern Europe. These findings will inform the development of more efficient FANCM genetic testing strategies for breast cancer cases from specific European populations.

Evaluation of European-based polygenic risk score for breast cancer in Ashkenazi Jewish women in Israel.

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.

Association of two genomic variants with HPV type-specific risk of cervical cancer.

PROBLEM: Human papillomavirus infection is integral to developing invasive cervical cancer in the majority of patients. In a recent genome-wide association study, rs9357152 and rs4243652 have been associated with seropositivity for HPV16 or HPV18, respectively. It is unknown whether these variants also associate with cervical cancer triggered by either HPV16 or HPV18. METHODS: We investigate whether the two HPV susceptibility variants show association with type-specific cervical cancer in a genetic case-control study with cases stratified by HPV16 or HPV18, respectively. We further tested whether rs9357152 modulates gene expression of any of 36 genes at the human leukocyte antigen locus in 256 cervical tissues. RESULTS: rs9357152 was associated with invasive HPV16-positive cervical cancer (OR 1.33, 95%CI 1.03-1.70, p = 0.03), and rs4243652 was associated with HPV18-positive adenocarcinomas (OR 2.96, 95%CI 1.18-7.41, p = 0.02). These associations remained borderline significant after testing against different sets of controls. rs9357152 was found to be an eQTL for HLA-DRB1 in HPV-positive cervical tissues (pANOVA = 0.0009), with the risk allele lowering mRNA levels. CONCLUSIONS: We find evidence that HPV seropositivity variants at chromosome 6 and 14 may modulate type-specific cervical cancer risk. rs9357152 may exert its effect through regulating HLA-DRB1 induction in the presence of HPV. In regard of multiple testing, these results need to be confirmed in larger studies.

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.

Clinical parameters affecting survival outcomes in patients with low-grade serous ovarian carcinoma: an international multicentre analysis.

BACKGROUND: Women with low-grade ovarian serous carcinoma (LGSC) benefit from surgical treatment; however, the role of chemotherapy is controversial. We examined an international database through the Ovarian Cancer Association Consortium to identify factors that affect survival in LGSC. METHODS: We performed a retrospective cohort analysis of patients with LGSC who had had primary surgery and had overall survival data available. We performed univariate and multivariate analyses of progression-free survival and overall survival, and generated Kaplan-Meier survival curves. RESULTS: Of the 707 patients with LGSC, 680 (96.2%) had available overall survival data. The patients' median age overall was 54 years. Of the 659 patients with International Federation of Obstetrics and Gynecology stage data, 156 (23.7%) had stage I disease, 64 (9.7%) had stage II, 395 (59.9%) had stage III, and 44 (6.7%) had stage IV. Of the 377 patients with surgical data, 200 (53.0%) had no visible residual disease. Of the 361 patients with chemotherapy data, 330 (91.4%) received first-line platinum-based chemotherapy. The median follow-up duration was 5.0 years. The median progression-free survival and overall survival were 43.2 months and 110.4 months, respectively. Multivariate analysis indicated a statistically significant impact of stage and residual disease on progression-free survival and overall survival. Platinum-based chemotherapy was not associated with a survival advantage. CONCLUSION: This multicentre analysis indicates that complete surgical cytoreduction to no visible residual disease has the most impact on improved survival in LGSC. This finding could immediately inform and change practice.

Multi-trait genome-wide association study identifies a novel endometrial cancer risk locus that associates with testosterone levels.

To detect novel endometrial cancer risk variants, we leveraged information from endometrial cancer risk factors in a multi-trait GWAS analysis. We first assessed causal relationships between established and suspected endometrial cancer risk factors, and endometrial cancer using Mendelian randomization. Following multivariable analysis, five independent risk factors (waist circumference, testosterone levels, sex hormone binding globulin levels, age at menarche, and age at natural menopause) were included in a multi-trait Bayesian GWAS analysis. We identified three potentially novel loci that associate with endometrial cancer risk, one of which (7q22.1) replicated in an independent endometrial cancer GWAS dataset and was genome-wide significant in a meta-analysis. This locus may affect endometrial cancer risk through altered testosterone levels. Consistent with this, we observed colocalization between the signals for endometrial cancer risk and expression of CYP3A7, a gene involved in testosterone metabolism. Thus, our findings suggest opportunities for hormone therapy to prevent or treat endometrial cancer.

Polymorphisms of the hypothalamic-pituitary-adrenal axis may lead to an inadequate response to stress and contribute to sudden infant death syndrome.

AIM: Impaired resilience to stress may be a factor in sudden infant death syndrome (SIDS). However, no comprehensive studies have been performed on polymorphisms that are relevant to the hypothalamic-pituitary-adrenal (HPA) axis, which regulates the stress hormone cortisol. METHODS: We analysed 22 relevant single nucleotide polymorphisms (SNPs) in 206 anonymised SIDS cases who died at a mean of 131 days (range: 5-343) and 256 adult controls who were recruited from paternity testing cases. Additional stratified analyses were performed for sex, age and season of death. Both the cases and the controls were Caucasian. RESULTS: Variants for rs2235543 (HSD11B1) and rs3779250 (CRHR2) were associated with SIDS in the overall analysis, and borderline for rs2446432 (CRH), at least before corrections for multiple testing. A combination of these three variants was observed in 52.9% of SIDS cases but only 43.0% of controls (p = 0.039). Five or more variants showed an association in the subgroups. CONCLUSION: Our findings suggest that the HPA axis influences SIDS and supports the hypothesis that an inadequate stress response may add to the risk. The associated variants for rs2235543, rs3779250 and rs2446432 appeared to decrease the cortisol concentration and impair an appropriate stress response.

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.

Aggregation tests identify new gene associations with breast cancer in populations with diverse ancestry.

BACKGROUND: Low-frequency variants play an important role in breast cancer (BC) susceptibility. Gene-based methods can increase power by combining multiple variants in the same gene and help identify target genes. METHODS: We evaluated the potential of gene-based aggregation in the Breast Cancer Association Consortium cohorts including 83,471 cases and 59,199 controls. Low-frequency variants were aggregated for individual genes' coding and regulatory regions. Association results in European ancestry samples were compared to single-marker association results in the same cohort. Gene-based associations were also combined in meta-analysis across individuals with European, Asian, African, and Latin American and Hispanic ancestry. RESULTS: In European ancestry samples, 14 genes were significantly associated (q < 0.05) with BC. Of those, two genes, FMNL3 (P = 6.11 × 10-6) and AC058822.1 (P = 1.47 × 10-4), represent new associations. High FMNL3 expression has previously been linked to poor prognosis in several other cancers. Meta-analysis of samples with diverse ancestry discovered further associations including established candidate genes ESR1 and CBLB. Furthermore, literature review and database query found further support for a biologically plausible link with cancer for genes CBLB, FMNL3, FGFR2, LSP1, MAP3K1, and SRGAP2C. CONCLUSIONS: Using extended gene-based aggregation tests including coding and regulatory variation, we report identification of plausible target genes for previously identified single-marker associations with BC as well as the discovery of novel genes implicated in BC development. Including multi ancestral cohorts in this study enabled the identification of otherwise missed disease associations as ESR1 (P = 1.31 × 10-5), demonstrating the importance of diversifying study cohorts.

FANCM missense variants and breast cancer risk: a case-control association study of 75,156 European women.

Evidence from literature, including the BRIDGES study, indicates that germline protein truncating variants (PTVs) in FANCM confer moderately increased risk of ER-negative and triple-negative breast cancer (TNBC), especially for women with a family history of the disease. Association between FANCM missense variants (MVs) and breast cancer risk has been postulated. In this study, we further used the BRIDGES study to test 689 FANCM MVs for association with breast cancer risk, overall and in ER-negative and TNBC subtypes, in 39,885 cases (7566 selected for family history) and 35,271 controls of European ancestry. Sixteen common MVs were tested individually; the remaining rare 673 MVs were tested by burden analyses considering their position and pathogenicity score. We also conducted a meta-analysis of our results and those from published studies. We did not find evidence for association for any of the 16 variants individually tested. The rare MVs were significantly associated with increased risk of ER-negative breast cancer by burden analysis comparing familial cases to controls (OR = 1.48; 95% CI 1.07-2.04; P = 0.017). Higher ORs were found for the subgroup of MVs located in functional domains or predicted to be pathogenic. The meta-analysis indicated that FANCM MVs overall are associated with breast cancer risk (OR = 1.22; 95% CI 1.08-1.38; P = 0.002). Our results support the definition from previous analyses of FANCM as a moderate-risk breast cancer gene and provide evidence that FANCM MVs could be low/moderate risk factors for ER-negative and TNBC subtypes. Further genetic and functional analyses are necessary to clarify better the increased risks due to FANCM MVs.

Associations of a Breast Cancer Polygenic Risk Score With Tumor Characteristics and Survival.

PURPOSE: A polygenic risk score (PRS) consisting of 313 common genetic variants (PRS313) is associated with risk of breast cancer and contralateral breast cancer. This study aimed to evaluate the association of the PRS313 with clinicopathologic characteristics of, and survival following, breast cancer. METHODS: Women with invasive breast cancer were included, 98,397 of European ancestry and 12,920 of Asian ancestry, from the Breast Cancer Association Consortium (BCAC), and 683 women from the European MINDACT trial. Associations between PRS313 and clinicopathologic characteristics, including the 70-gene signature for MINDACT, were evaluated using logistic regression analyses. Associations of PRS313 (continuous, per standard deviation) with overall survival (OS) and breast cancer-specific survival (BCSS) were evaluated with Cox regression, adjusted for clinicopathologic characteristics and treatment. RESULTS: The PRS313 was associated with more favorable tumor characteristics. In BCAC, increasing PRS313 was associated with lower grade, hormone receptor-positive status, and smaller tumor size. In MINDACT, PRS313 was associated with a low risk 70-gene signature. In European women from BCAC, higher PRS313 was associated with better OS and BCSS: hazard ratio (HR) 0.96 (95% CI, 0.94 to 0.97) and 0.96 (95% CI, 0.94 to 0.98), but the association disappeared after adjustment for clinicopathologic characteristics (and treatment): OS HR, 1.01 (95% CI, 0.98 to 1.05) and BCSS HR, 1.02 (95% CI, 0.98 to 1.07). The results in MINDACT and Asian women from BCAC were consistent. CONCLUSION: An increased PRS313 is associated with favorable tumor characteristics, but is not independently associated with prognosis. Thus, PRS313 has no role in the clinical management of primary breast cancer at the time of diagnosis. Nevertheless, breast cancer mortality rates will be higher for women with higher PRS313 as increasing PRS313 is associated with an increased risk of disease. This information is crucial for modeling effective stratified screening programs.