Genome-wide association studies (GWAS) have identified more than 170 breast cancer susceptibility loci. Here we hypothesize that some risk-associated variants might act in non-breast tissues, specifically adipose tissue and immune cells from blood and spleen. Using expression quantitative trait loci (eQTL) reported in these tissues, we identify 26 previously unreported, likely target genes of overall breast cancer risk variants, and 17 for estrogen receptor (ER)-negative breast cancer, several with a known immune function. We determine the directional effect of gene expression on disease risk measured based on single and multiple eQTL. In addition, using a gene-based test of association that considers eQTL from multiple tissues, we identify seven (and four) regions with variants associated with overall (and ER-negative) breast cancer risk, which were not reported in previous GWAS. Further investigation of the function of the implicated genes in breast and immune cells may provide insights into the etiology of breast cancer.

BRCA1 BRCA2 mutational spectrum in the Middle East, North Africa, and Southern Europe is not well characterized. The unique history and cultural practices characterizing these regions, often involving consanguinity and inbreeding, plausibly led to the accumulation of population-specific founder pathogenic sequence variants (PSVs). To determine recurring BRCA PSVs in these locales, a search in PUBMED, EMBASE, BIC, and CIMBA was carried out combined with outreach to researchers from the relevant countries for unpublished data. We identified 232 PSVs in BRCA1 and 239 in BRCA2 in 25 of 33 countries surveyed. Common PSVs that were detected in four or more countries were c.5266dup (p.Gln1756Profs), c.181T>G (p.Cys61Gly), c.68_69del (p.Glu23Valfs), c.5030_5033del (p.Thr1677Ilefs), c.4327C>T (p.Arg1443Ter), c.5251C>T (p.Arg1751Ter), c.1016dup (p.Val340Glyfs), c.3700_3704del (p.Val1234Glnfs), c.4065_4068del (p.Asn1355Lysfs), c.1504_1508del (p.Leu502Alafs), c.843_846del (p.Ser282Tyrfs), c.798_799del (p.Ser267Lysfs), and c.3607C>T (p.Arg1203Ter) in BRCA1 and c.2808_2811del (p.Ala938Profs), c.5722_5723del (p.Leu1908Argfs), c.9097dup (p.Thr3033Asnfs), c.1310_1313del (p. p.Lys437Ilefs), and c.5946del (p.Ser1982Argfs) for BRCA2. Notably, some mutations (e.g., p.Asn257Lysfs (c.771_775del)) were observed in unrelated populations. Thus, seemingly genotyping recurring BRCA PSVs in specific populations may provide first pass BRCA genotyping platform.

Breast cancer is a common disease partially caused by genetic risk factors. Germline pathogenic variants in DNA repair genes BRCA1, BRCA2, PALB2, ATM, and CHEK2 are associated with breast cancer risk. FANCM, which encodes for a DNA translocase, has been proposed as a breast cancer predisposition gene, with greater effects for the ER-negative and triple-negative breast cancer (TNBC) subtypes. We tested the three recurrent protein-truncating variants FANCM:p.Arg658*, p.Gln1701*, and p.Arg1931* for association with breast cancer risk in 67,112 cases, 53,766 controls, and 26,662 carriers of pathogenic variants of BRCA1 or BRCA2. These three variants were also studied functionally by measuring survival and chromosome fragility in FANCM^−/− patient-derived immortalized fibroblasts treated with diepoxybutane or olaparib. We observed that FANCM:p.Arg658* was associated with increased risk of ER-negative disease and TNBC (OR = 2.44, P = 0.034 and OR = 3.79; P = 0.009, respectively). In a country-restricted analysis, we confirmed the associations detected for FANCM:p.Arg658* and found that also FANCM:p.Arg1931* was associated with ER-negative breast cancer risk (OR = 1.96; P = 0.006). The functional results indicated that all three variants were deleterious affecting cell survival and chromosome stability with FANCM:p.Arg658* causing more severe phenotypes. In conclusion, we confirmed that the two rare FANCM deleterious variants p.Arg658* and p.Arg1931* are risk factors for ER-negative and TNBC subtypes. Overall our data suggest that the effect of truncating variants on breast cancer risk may depend on their position in the gene. Cell sensitivity to olaparib exposure, identifies a possible therapeutic option to treat FANCM-associated tumors.

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 analysis 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 pathogenicityfindings that would not have been reached using other case-control analysis methods. Our novel method provides greater power to classify rare variants compared with classical case-control methods. As an initiative from the ENIGMA Analytical Working Group, we provide user-friendly scripts and preformatted Excel calculators for implementation of the method for rare variants in BRCA1, BRCA2, and other high-risk genes with known penetrance.