Validation of the BOADICEA model in a prospective cohort of BRCA1/2 pathogenic variant carriers.

BACKGROUND: No validation has been conducted for the BOADICEA multifactorial breast cancer risk prediction model specifically in BRCA1/2 pathogenic variant (PV) carriers to date. Here, we evaluated the performance of BOADICEA in predicting 5-year breast cancer risks in a prospective cohort of BRCA1/2 PV carriers ascertained through clinical genetic centres. METHODS: We evaluated the model calibration and discriminatory ability in the prospective TRANsIBCCS cohort study comprising 1614 BRCA1 and 1365 BRCA2 PV carriers (209 incident cases). Study participants had lifestyle, reproductive, hormonal, anthropometric risk factor information, a polygenic risk score based on 313 SNPs and family history information. RESULTS: The full multifactorial model considering family history together with all other risk factors was well calibrated overall (E/O=1.07, 95% CI: 0.92 to 1.24) and in quintiles of predicted risk. Discrimination was maximised when all risk factors were considered (Harrell's C-index=0.70, 95% CI: 0.67 to 0.74; area under the curve=0.79, 95% CI: 0.76 to 0.82). The model performance was similar when evaluated separately in BRCA1 or BRCA2 PV carriers. The full model identified 5.8%, 12.9% and 24.0% of BRCA1/2 PV carriers with 5-year breast cancer risks of <1.65%, <3% and <5%, respectively, risk thresholds commonly used for different management and risk-reduction options. CONCLUSION: BOADICEA may be used to aid personalised cancer risk management and decision-making for BRCA1 and BRCA2 PV carriers. It is implemented in the free-access CanRisk tool (https://www.canrisk.org/).

Associations of height, body mass index, and weight gain with breast cancer risk in carriers of a pathogenic variant in BRCA1 or BRCA2: the BRCA1 and BRCA2 Cohort Consortium.

INTRODUCTION: Height, body mass index (BMI), and weight gain are associated with breast cancer risk in the general population. It is unclear whether these associations also exist for carriers of pathogenic variants in the BRCA1 or BRCA2 genes. PATIENTS AND METHODS: An international pooled cohort of 8091 BRCA1/2 variant carriers was used for retrospective and prospective analyses separately for premenopausal and postmenopausal women. Cox regression was used to estimate breast cancer risk associations with height, BMI, and weight change. RESULTS: In the retrospective analysis, taller height was associated with risk of premenopausal breast cancer for BRCA2 variant carriers (HR 1.20 per 10 cm increase, 95% CI 1.04-1.38). Higher young-adult BMI was associated with lower premenopausal breast cancer risk for both BRCA1 (HR 0.75 per 5 kg/m2, 95% CI 0.66-0.84) and BRCA2 (HR 0.76, 95% CI 0.65-0.89) variant carriers in the retrospective analysis, with consistent, though not statistically significant, findings from the prospective analysis. In the prospective analysis, higher BMI and adult weight gain were associated with higher postmenopausal breast cancer risk for BRCA1 carriers (HR 1.20 per 5 kg/m2, 95% CI 1.02-1.42; and HR 1.10 per 5 kg weight gain, 95% CI 1.01-1.19, respectively). CONCLUSION: Anthropometric measures are associated with breast cancer risk for BRCA1 and BRCA2 variant carriers, with relative risk estimates that are generally consistent with those for women from the general population.

Boosting GWAS using biological networks: A study on susceptibility to familial breast cancer.

Genome-wide association studies (GWAS) explore the genetic causes of complex diseases. However, classical approaches ignore the biological context of the genetic variants and genes under study. To address this shortcoming, one can use biological networks, which model functional relationships, to search for functionally related susceptibility loci. Many such network methods exist, each arising from different mathematical frameworks, pre-processing steps, and assumptions about the network properties of the susceptibility mechanism. Unsurprisingly, this results in disparate solutions. To explore how to exploit these heterogeneous approaches, we selected six network methods and applied them to GENESIS, a nationwide French study on familial breast cancer. First, we verified that network methods recovered more interpretable results than a standard GWAS. We addressed the heterogeneity of their solutions by studying their overlap, computing what we called the consensus. The key gene in this consensus solution was COPS5, a gene related to multiple cancer hallmarks. Another issue we observed was that network methods were unstable, selecting very different genes on different subsamples of GENESIS. Therefore, we proposed a stable consensus solution formed by the 68 genes most consistently selected across multiple subsamples. This solution was also enriched in genes known to be associated with breast cancer susceptibility (BLM, CASP8, CASP10, DNAJC1, FGFR2, MRPS30, and SLC4A7, P-value = 3 × 10-4). The most connected gene was CUL3, a regulator of several genes linked to cancer progression. Lastly, we evaluated the biases of each method and the impact of their parameters on the outcome. In general, network methods preferred highly connected genes, even after random rewirings that stripped the connections of any biological meaning. In conclusion, we present the advantages of network-guided GWAS, characterize their shortcomings, and provide strategies to address them. To compute the consensus networks, implementations of all six methods are available at https://github.com/hclimente/gwas-tools.

Gene- and pathway-level analyses of iCOGS variants highlight novel signaling pathways underlying familial breast cancer susceptibility.

Single-nucleotide polymorphisms (SNPs) in over 180 loci have been associated with breast cancer (BC) through genome-wide association studies involving mostly unselected population-based case-control series. Some of them modify BC risk of women carrying a BRCA1 or BRCA2 (BRCA1/2) mutation and may also explain BC risk variability in BC-prone families with no BRCA1/2 mutation. Here, we assessed the contribution of SNPs of the iCOGS array in GENESIS consisting of BC cases with no BRCA1/2 mutation and a sister with BC, and population controls. Genotyping data were available for 1281 index cases, 731 sisters with BC, 457 unaffected sisters and 1272 controls. In addition to the standard SNP-level analysis using index cases and controls, we performed pedigree-based association tests to capture transmission information in the sibships. We also performed gene- and pathway-level analyses to maximize the power to detect associations with lower-frequency SNPs or those with modest effect sizes. While SNP-level analyses identified 18 loci, gene-level analyses identified 112 genes. Furthermore, 31 Kyoto Encyclopedia of Genes and Genomes and 7 Atlas of Cancer Signaling Network pathways were highlighted (false discovery rate of 5%). Using results from the "index case-control" analysis, we built pathway-derived polygenic risk scores (PRS) and assessed their performance in the population-based CECILE study and in a data set composed of GENESIS-affected sisters and CECILE controls. Although these PRS had poor predictive value in the general population, they performed better than a PRS built using our SNP-level findings, and we found that the joint effect of family history and PRS needs to be considered in risk prediction models.

Oral contraceptive use and ovarian cancer risk for BRCA1/2 mutation carriers: an international cohort study.

BACKGROUND: Ovarian cancer risk in BRCA1 and BRCA2 mutation carriers has been shown to decrease with longer duration of oral contraceptive use. Although the effects of using oral contraceptives in the general population are well established (approximately 50% risk reduction in ovarian cancer), the estimated risk reduction in mutation carriers is much less precise because of potential bias and small sample sizes. In addition, only a few studies on oral contraceptive use have examined the associations of duration of use, time since last use, starting age, and calendar year of start with risk of ovarian cancer. OBJECTIVE: This study aimed to investigate in more detail the associations of various characteristics of oral contraceptive use and risk of ovarian cancer, to provide healthcare providers and carriers with better risk estimates. STUDY DESIGN: In this international retrospective study, ovarian cancer risk associations were assessed using oral contraceptives data on 3989 BRCA1 and 2445 BRCA2 mutation carriers. Age-dependent-weighted Cox regression analyses were stratified by study and birth cohort and included breast cancer diagnosis as a covariate. To minimize survival bias, analyses were left truncated at 5 years before baseline questionnaire. Separate analyses were conducted for each aspect of oral contraceptive use and in a multivariate analysis, including all these aspects. In addition, the analysis of duration of oral contraceptive use was stratified by recency of use. RESULTS: Oral contraceptives were less often used by mutation carriers who were diagnosed with ovarian cancer (ever use: 58.6% for BRCA1 and 53.5% BRCA2) than by unaffected carriers (ever use: 88.9% for BRCA1 and 80.7% for BRCA2). The median duration of use was 7 years for both BRCA1 and BRCA2 carriers who developed ovarian cancer and 9 and 8 years for unaffected BRCA1 and BRCA2 carriers with ovarian cancer, respectively. For BRCA1 mutation carriers, univariate analyses have shown that both a longer duration of oral contraceptive use and more recent oral contraceptive use were associated with a reduction in the risk of ovarian cancer. However, in multivariate analyses, including duration of use, age at first use, and time since last use, duration of oral contraceptive use proved to be the prominent protective factor (compared with <5 years: 5-9 years [hazard ratio, 0.67; 95% confidence interval, 0.40-1.12]; >10 years [hazard ratio, 0.37; 95% confidence interval, 0.19-0.73]; Ptrend=.008). The inverse association between duration of use and ovarian cancer risk persisted for more than 15 years (duration of ≥10 years; BRCA1 <15 years since last use [hazard ratio, 0.24; 95% confidence interval, 0.14-0.43]; BRCA1 >15 years since last use [hazard ratio, 0.56; 95% confidence interval, 0.18-0.59]). Univariate results for BRCA2 mutation carriers were similar but were inconclusive because of limited sample size. CONCLUSION: For BRCA1 mutation carriers, longer duration of oral contraceptive use is associated with a greater reduction in ovarian cancer risk, and the protection is long term.

A new hybrid record linkage process to make epidemiological databases interoperable: application to the GEMO and GENEPSO studies involving BRCA1 and BRCA2 mutation carriers.

BACKGROUND: Linking independent sources of data describing the same individuals enable innovative epidemiological and health studies but require a robust record linkage approach. We describe a hybrid record linkage process to link databases from two independent ongoing French national studies, GEMO (Genetic Modifiers of BRCA1 and BRCA2), which focuses on the identification of genetic factors modifying cancer risk of BRCA1 and BRCA2 mutation carriers, and GENEPSO (prospective cohort of BRCAx mutation carriers), which focuses on environmental and lifestyle risk factors. METHODS: To identify as many as possible of the individuals participating in the two studies but not registered by a shared identifier, we combined probabilistic record linkage (PRL) and supervised machine learning (ML). This approach (named "PRL + ML") combined together the candidate matches identified by both approaches. We built the ML model using the gold standard on a first version of the two databases as a training dataset. This gold standard was obtained from PRL-derived matches verified by an exhaustive manual review. Results The Random Forest (RF) algorithm showed a highest recall (0.985) among six widely used ML algorithms: RF, Bagged trees, AdaBoost, Support Vector Machine, Neural Network. Therefore, RF was selected to build the ML model since our goal was to identify the maximum number of true matches. Our combined linkage PRL + ML showed a higher recall (range 0.988-0.992) than either PRL (range 0.916-0.991) or ML (0.981) alone. It identified 1995 individuals participating in both GEMO (6375 participants) and GENEPSO (4925 participants). CONCLUSIONS: Our hybrid linkage process represents an efficient tool for linking GEMO and GENEPSO. It may be generalizable to other epidemiological studies involving other databases and registries.

Polygenic risk scores and breast and epithelial ovarian cancer risks for carriers of BRCA1 and BRCA2 pathogenic variants.

PURPOSE: We assessed the associations between population-based polygenic risk scores (PRS) for breast (BC) or epithelial ovarian cancer (EOC) with cancer risks for BRCA1 and BRCA2 pathogenic variant carriers. METHODS: Retrospective cohort data on 18,935 BRCA1 and 12,339 BRCA2 female pathogenic variant carriers of European ancestry were available. Three versions of a 313 single-nucleotide polymorphism (SNP) BC PRS were evaluated based on whether they predict overall, estrogen receptor (ER)-negative, or ER-positive BC, and two PRS for overall or high-grade serous EOC. Associations were validated in a prospective cohort. RESULTS: The ER-negative PRS showed the strongest association with BC risk for BRCA1 carriers (hazard ratio [HR] per standard deviation = 1.29 [95% CI 1.25-1.33], P = 3×10-72). For BRCA2, the strongest association was with overall BC PRS (HR = 1.31 [95% CI 1.27-1.36], P = 7×10-50). HR estimates decreased significantly with age and there was evidence for differences in associations by predicted variant effects on protein expression. The HR estimates were smaller than general population estimates. The high-grade serous PRS yielded the strongest associations with EOC risk for BRCA1 (HR = 1.32 [95% CI 1.25-1.40], P = 3×10-22) and BRCA2 (HR = 1.44 [95% CI 1.30-1.60], P = 4×10-12) carriers. The associations in the prospective cohort were similar. CONCLUSION: Population-based PRS are strongly associated with BC and EOC risks for BRCA1/2 carriers and predict substantial absolute risk differences for women at PRS distribution extremes.

Morphology and genomic hallmarks of breast tumours developed by ATM deleterious variant carriers.

BACKGROUND: The ataxia telangiectasia mutated (ATM) gene is a moderate-risk breast cancer susceptibility gene; germline loss-of-function variants are found in up to 3% of hereditary breast and ovarian cancer (HBOC) families who undergo genetic testing. So far, no clear histopathological and molecular features of breast tumours occurring in ATM deleterious variant carriers have been described, but identification of an ATM-associated tumour signature may help in patient management. METHODS: To characterise hallmarks of ATM-associated tumours, we performed systematic pathology review of tumours from 21 participants from ataxia-telangiectasia families and 18 participants from HBOC families, as well as copy number profiling on a subset of 23 tumours. Morphology of ATM-associated tumours was compared with that of 599 patients with no BRCA1 and BRCA2 mutations from a hospital-based series, as well as with data from The Cancer Genome Atlas. Absolute copy number and loss of heterozygosity (LOH) profiles were obtained from the OncoScan SNP array. In addition, we performed whole-genome sequencing on four tumours from ATM loss-of-function variant carriers with available frozen material. RESULTS: We found that ATM-associated tumours belong mostly to the luminal B subtype, are tetraploid and show LOH at the ATM locus at 11q22-23. Unlike tumours in which BRCA1 or BRCA2 is inactivated, tumours arising in ATM deleterious variant carriers are not associated with increased large-scale genomic instability as measured by the large-scale state transitions signature. Losses at 13q14.11-q14.3, 17p13.2-p12, 21p11.2-p11.1 and 22q11.23 were observed. Somatic alterations at these loci may therefore represent biomarkers for ATM testing and harbour driver mutations in potentially 'druggable' genes that would allow patients to be directed towards tailored therapeutic strategies. CONCLUSIONS: Although ATM is involved in the DNA damage response, ATM-associated tumours are distinct from BRCA1-associated tumours in terms of morphological characteristics and genomic alterations, and they are also distinguishable from sporadic breast tumours, thus opening up the possibility to identify ATM variant carriers outside the ataxia-telangiectasia disorder and direct them towards effective cancer risk management and therapeutic strategies.

Telomere length, ATM mutation status and cancer risk in Ataxia-Telangiectasia families.

Recent studies have linked constitutive telomere length (TL) to aging-related diseases including cancer at different sites. ATM participates in the signaling of telomere erosion, and inherited mutations in ATM have been associated with increased risk of cancer, particularly breast cancer. The goal of this study was to investigate whether carriage of an ATM mutation and TL interplay to modify cancer risk in ataxia-telangiectasia (A-T) families.The study population consisted of 284 heterozygous ATM mutation carriers (HetAT) and 174 non-carriers (non-HetAT) from 103 A-T families. Forty-eight HetAT and 14 non-HetAT individuals had cancer, among them 25 HetAT and 6 non-HetAT were diagnosed after blood sample collection. We measured mean TL using a quantitative PCR assay and genotyped seven single-nucleotide polymorphisms (SNPs) recurrently associated with TL in large population-based studies.HetAT individuals were at increased risk of cancer (OR = 2.3, 95%CI = 1.2-4.4, P = 0.01), and particularly of breast cancer for women (OR = 2.9, 95%CI = 1.2-7.1, P = 0.02), in comparison to their non-HetAT relatives. HetAT individuals had longer telomeres than non-HetAT individuals (P = 0.0008) but TL was not associated with cancer risk, and no significant interaction was observed between ATM mutation status and TL. Furthermore, rs9257445 (ZNF311) was associated with TL in HetAT subjects and rs6060627 (BCL2L1) modified cancer risk in HetAT and non-HetAT women.Our findings suggest that carriage of an ATM mutation impacts on the age-related TL shortening and that TL per se is not related to cancer risk in ATM carriers. TL measurement alone is not a good marker for predicting cancer risk in A-T families.

FANCM c.5791C>T nonsense mutation (rs144567652) induces exon skipping, affects DNA repair activity and is a familial breast cancer risk factor.

Numerous genetic factors that influence breast cancer risk are known. However, approximately two-thirds of the overall familial risk remain unexplained. To determine whether some of the missing heritability is due to rare variants conferring high to moderate risk, we tested for an association between the c.5791C>T nonsense mutation (p.Arg1931*; rs144567652) in exon 22 of FANCM gene and breast cancer. An analysis of genotyping data from 8635 familial breast cancer cases and 6625 controls from different countries yielded an association between the c.5791C>T mutation and breast cancer risk [odds ratio (OR) = 3.93 (95% confidence interval (CI) = 1.28-12.11; P = 0.017)]. Moreover, we performed two meta-analyses of studies from countries with carriers in both cases and controls and of all available data. These analyses showed breast cancer associations with OR = 3.67 (95% CI = 1.04-12.87; P = 0.043) and OR = 3.33 (95% CI = 1.09-13.62; P = 0.032), respectively. Based on information theory-based prediction, we established that the mutation caused an out-of-frame deletion of exon 22, due to the creation of a binding site for the pre-mRNA processing protein hnRNP A1. Furthermore, genetic complementation analyses showed that the mutation influenced the DNA repair activity of the FANCM protein. In summary, we provide evidence for the first time showing that the common p.Arg1931* loss-of-function variant in FANCM is a risk factor for familial breast cancer.