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.

Ovarian cancer pathology characteristics as predictors of variant pathogenicity in BRCA1 and BRCA2.

BACKGROUND: The distribution of ovarian tumour characteristics differs between germline BRCA1 and BRCA2 pathogenic variant carriers and non-carriers. In this study, we assessed the utility of ovarian tumour characteristics as predictors of BRCA1 and BRCA2 variant pathogenicity, for application using the American College of Medical Genetics and the Association for Molecular Pathology (ACMG/AMP) variant classification system. METHODS: Data for 10,373 ovarian cancer cases, including carriers and non-carriers of BRCA1 or BRCA2 pathogenic variants, were collected from unpublished international cohorts and consortia and published studies. Likelihood ratios (LR) were calculated for the association of ovarian cancer histology and other characteristics, with BRCA1 and BRCA2 variant pathogenicity. Estimates were aligned to ACMG/AMP code strengths (supporting, moderate, strong). RESULTS: No histological subtype provided informative ACMG/AMP evidence in favour of BRCA1 and BRCA2 variant pathogenicity. Evidence against variant pathogenicity was estimated for the mucinous and clear cell histologies (supporting) and borderline cases (moderate). Refined associations are provided according to tumour grade, invasion and age at diagnosis. CONCLUSIONS: We provide detailed estimates for predicting BRCA1 and BRCA2 variant pathogenicity based on ovarian tumour characteristics. This evidence can be combined with other variant information under the ACMG/AMP classification system, to improve classification and carrier clinical management.

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.

CDK5RAP3, a New BRCA2 Partner That Regulates DNA Repair, Is Associated with Breast Cancer Survival.

BRCA2 is essential for homologous recombination DNA repair. BRCA2 mutations lead to genome instability and increased risk of breast and ovarian cancer. Similarly, mutations in BRCA2-interacting proteins are also known to modulate sensitivity to DNA damage agents and are established cancer risk factors. Here we identify the tumor suppressor CDK5RAP3 as a novel BRCA2 helical domain-interacting protein. CDK5RAP3 depletion induced DNA damage resistance, homologous recombination and single-strand annealing upregulation, and reduced spontaneous and DNA damage-induced genomic instability, suggesting that CDK5RAP3 negatively regulates double-strand break repair in the S-phase. Consistent with this cellular phenotype, analysis of transcriptomic data revealed an association between low CDK5RAP3 tumor expression and poor survival of breast cancer patients. Finally, we identified common genetic variations in the CDK5RAP3 locus as potentially associated with breast and ovarian cancer risk in BRCA1 and BRCA2 mutation carriers. Our results uncover CDK5RAP3 as a critical player in DNA repair and breast cancer outcomes.

Evaluation of the association of heterozygous germline variants in NTHL1 with breast cancer predisposition: an international multi-center study of 47,180 subjects.

Bi-allelic loss-of-function (LoF) variants in the base excision repair (BER) gene NTHL1 cause a high-risk hereditary multi-tumor syndrome that includes breast cancer, but the contribution of heterozygous variants to hereditary breast cancer is unknown. An analysis of 4985 women with breast cancer, enriched for familial features, and 4786 cancer-free women revealed significant enrichment for NTHL1 LoF variants. Immunohistochemistry confirmed reduced NTHL1 expression in tumors from heterozygous carriers but the NTHL1 bi-allelic loss characteristic mutational signature (SBS 30) was not present. The analysis was extended to 27,421 breast cancer cases and 19,759 controls from 10 international studies revealing 138 cases and 93 controls with a heterozygous LoF variant (OR 1.06, 95% CI: 0.82-1.39) and 316 cases and 179 controls with a missense variant (OR 1.31, 95% CI: 1.09-1.57). Missense variants selected for deleterious features by a number of in silico bioinformatic prediction tools or located within the endonuclease III functional domain showed a stronger association with breast cancer. Somatic sequencing of breast cancers from carriers indicated that the risk associated with NTHL1 appears to operate through haploinsufficiency, consistent with other described low-penetrance breast cancer genes. Data from this very large international multicenter study suggests that heterozygous pathogenic germline coding variants in NTHL1 may be associated with low- to moderate- increased risk of breast cancer.

Association of germline variation with the survival of women with BRCA1/2 pathogenic variants and breast cancer.

Germline genetic variation has been suggested to influence the survival of breast cancer patients independently of tumor pathology. We have studied survival associations of genetic variants in two etiologically unique groups of breast cancer patients, the carriers of germline pathogenic variants in BRCA1 or BRCA2 genes. We found that rs57025206 was significantly associated with the overall survival, predicting higher mortality of BRCA1 carrier patients with estrogen receptor-negative breast cancer, with a hazard ratio 4.37 (95% confidence interval 3.03-6.30, P = 3.1 × 10-9). Multivariable analysis adjusted for tumor characteristics suggested that rs57025206 was an independent survival marker. In addition, our exploratory analyses suggest that the associations between genetic variants and breast cancer patient survival may depend on tumor biological subgroup and clinical patient characteristics.

Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes.

Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.

Functional Analysis of Promoter Variants in Genes Involved in Sex Steroid Action, DNA Repair and Cell Cycle Control.

Genetic variants affecting the regulation of gene expression are among the main causes of human diversity. The potential importance of regulatory polymorphisms is underscored by results from Genome Wide Association Studies, which have already implicated such polymorphisms in the susceptibility to complex diseases such as breast cancer. In this study, we re-sequenced the promoter regions of 24 genes involved in pathways related to breast cancer including sex steroid action, DNA repair, and cell cycle control in 60 unrelated Caucasian individuals. We constructed haplotypes and assessed the functional impact of promoter variants using gene reporter assays and electrophoretic mobility shift assays. We identified putative functional variants within the promoter regions of estrogen receptor 1 (ESR1), ESR2, forkhead box A1 (FOXA1), ubiquitin interaction motif containing 1 (UIMC1) and cell division cycle 7 (CDC7). The functional polymorphism on CDC7, rs13447455, influences CDC7 transcriptional activity in an allele-specific manner and alters DNA⁻protein complex formation in breast cancer cell lines. Moreover, results from the Breast Cancer Association Consortium show a marginal association between rs13447455 and breast cancer risk (p=9.3x10-5), thus warranting further investigation. Furthermore, our study has helped provide methodological solutions to some technical difficulties that were encountered with gene reporter assays, particularly regarding inter-clone variability and statistical consistency.

DSNetwork: An Integrative Approach to Visualize Predictions of Variants' Deleteriousness.

One of the most challenging tasks of the post-genome-wide association studies (GWAS) research era is the identification of functional variants among those associated with a trait for an observed GWAS signal. Several methods have been developed to evaluate the potential functional implications of genetic variants. Each of these tools has its own scoring system, which forces users to become acquainted with each approach to interpret their results. From an awareness of the amount of work needed to analyze and integrate results for a single locus, we proposed a flexible and versatile approach designed to help the prioritization of variants by aggregating the predictions of their potential functional implications. This approach has been made available through a graphical user interface called DSNetwork, which acts as a single point of entry to almost 60 reference predictors for both coding and non-coding variants and displays predictions in an easy-to-interpret visualization. We confirmed the usefulness of our methodology by successfully identifying functional variants in four breast cancer and nine schizophrenia susceptibility loci.

Association analysis identifies 65 new breast cancer risk loci.

Breast cancer risk is influenced by rare coding variants in susceptibility genes, such as BRCA1, and many common, mostly non-coding variants. However, much of the genetic contribution to breast cancer risk remains unknown. Here we report the results of a genome-wide association study of breast cancer in 122,977 cases and 105,974 controls of European ancestry and 14,068 cases and 13,104 controls of East Asian ancestry. We identified 65 new loci that are associated with overall breast cancer risk at P < 5 × 10-8. The majority of credible risk single-nucleotide polymorphisms in these loci fall in distal regulatory elements, and by integrating in silico data to predict target genes in breast cells at each locus, we demonstrate a strong overlap between candidate target genes and somatic driver genes in breast tumours. We also find that heritability of breast cancer due to all single-nucleotide polymorphisms in regulatory features was 2-5-fold enriched relative to the genome-wide average, with strong enrichment for particular transcription factor binding sites. These results provide further insight into genetic susceptibility to breast cancer and will improve the use of genetic risk scores for individualized screening and prevention.

Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer.

Most common breast cancer susceptibility variants have been identified through genome-wide association studies (GWAS) of predominantly estrogen receptor (ER)-positive disease. We conducted a GWAS using 21,468 ER-negative cases and 100,594 controls combined with 18,908 BRCA1 mutation carriers (9,414 with breast cancer), all of European origin. We identified independent associations at P < 5 × 10-8 with ten variants at nine new loci. At P < 0.05, we replicated associations with 10 of 11 variants previously reported in ER-negative disease or BRCA1 mutation carrier GWAS and observed consistent associations with ER-negative disease for 105 susceptibility variants identified by other studies. These 125 variants explain approximately 16% of the familial risk of this breast cancer subtype. There was high genetic correlation (0.72) between risk of ER-negative breast cancer and breast cancer risk for BRCA1 mutation carriers. These findings may lead to improved risk prediction and inform further fine-mapping and functional work to better understand the biological basis of ER-negative breast cancer.

Prediction of Breast and Prostate Cancer Risks in Male BRCA1 and BRCA2 Mutation Carriers Using Polygenic Risk Scores.

Purpose BRCA1/2 mutations increase the risk of breast and prostate cancer in men. Common genetic variants modify cancer risks for female carriers of BRCA1/2 mutations. We investigated-for the first time to our knowledge-associations of common genetic variants with breast and prostate cancer risks for male carriers of BRCA1/ 2 mutations and implications for cancer risk prediction. Materials and Methods We genotyped 1,802 male carriers of BRCA1/2 mutations from the Consortium of Investigators of Modifiers of BRCA1/2 by using the custom Illumina OncoArray. We investigated the combined effects of established breast and prostate cancer susceptibility variants on cancer risks for male carriers of BRCA1/2 mutations by constructing weighted polygenic risk scores (PRSs) using published effect estimates as weights. Results In male carriers of BRCA1/2 mutations, PRS that was based on 88 female breast cancer susceptibility variants was associated with breast cancer risk (odds ratio per standard deviation of PRS, 1.36; 95% CI, 1.19 to 1.56; P = 8.6 × 10-6). Similarly, PRS that was based on 103 prostate cancer susceptibility variants was associated with prostate cancer risk (odds ratio per SD of PRS, 1.56; 95% CI, 1.35 to 1.81; P = 3.2 × 10-9). Large differences in absolute cancer risks were observed at the extremes of the PRS distribution. For example, prostate cancer risk by age 80 years at the 5th and 95th percentiles of the PRS varies from 7% to 26% for carriers of BRCA1 mutations and from 19% to 61% for carriers of BRCA2 mutations, respectively. Conclusion PRSs may provide informative cancer risk stratification for male carriers of BRCA1/2 mutations that might enable these men and their physicians to make informed decisions on the type and timing of breast and prostate cancer risk management.

VEXOR: an integrative environment for prioritization of functional variants in fine-mapping analysis.

Motivation: The identification of the functional variants responsible for observed genome-wide association studies (GWAS) signals is one of the most challenging tasks of the post-GWAS research era. Several tools have been developed to annotate genetic variants by their genomic location and potential functional implications. Each of these tools has its own requirements and internal logic, which forces the user to become acquainted with each interface. Results: From an awareness of the amount of work needed to analyze a single locus, we have built a flexible, versatile and easy-to-use web interface designed to help in prioritizing variants and predicting their potential functional implications. This interface acts as a single-point of entry linking association results with reference tools and relevant experiments. Availability and Implementation: VEXOR is an integrative web application implemented through the Shiny framework and available at: http://romix.genome.ulaval.ca/vexor. Contact: arnaud.droit@crchuq.ulaval.ca. Supplementary information: Supplementary data are available at Bioinformatics online.

Evaluation of Polygenic Risk Scores for Breast and Ovarian Cancer Risk Prediction in BRCA1 and BRCA2 Mutation Carriers.

Background: Genome-wide association studies (GWAS) have identified 94 common single-nucleotide polymorphisms (SNPs) associated with breast cancer (BC) risk and 18 associated with ovarian cancer (OC) risk. Several of these are also associated with risk of BC or OC for women who carry a pathogenic mutation in the high-risk BC and OC genes BRCA1 or BRCA2. The combined effects of these variants on BC or OC risk for BRCA1 and BRCA2 mutation carriers have not yet been assessed while their clinical management could benefit from improved personalized risk estimates. Methods: We constructed polygenic risk scores (PRS) using BC and OC susceptibility SNPs identified through population-based GWAS: for BC (overall, estrogen receptor [ER]-positive, and ER-negative) and for OC. Using data from 15 252 female BRCA1 and 8211 BRCA2 carriers, the association of each PRS with BC or OC risk was evaluated using a weighted cohort approach, with time to diagnosis as the outcome and estimation of the hazard ratios (HRs) per standard deviation increase in the PRS. Results: The PRS for ER-negative BC displayed the strongest association with BC risk in BRCA1 carriers (HR = 1.27, 95% confidence interval [CI] = 1.23 to 1.31, P = 8.2×10 -53 ). In BRCA2 carriers, the strongest association with BC risk was seen for the overall BC PRS (HR = 1.22, 95% CI = 1.17 to 1.28, P = 7.2×10 -20 ). The OC PRS was strongly associated with OC risk for both BRCA1 and BRCA2 carriers. These translate to differences in absolute risks (more than 10% in each case) between the top and bottom deciles of the PRS distribution; for example, the OC risk was 6% by age 80 years for BRCA2 carriers at the 10th percentile of the OC PRS compared with 19% risk for those at the 90th percentile of PRS. Conclusions: BC and OC PRS are predictive of cancer risk in BRCA1 and BRCA2 carriers. Incorporation of the PRS into risk prediction models has promise to better inform decisions on cancer risk management.

Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer.

To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC.

The OncoArray Consortium: A Network for Understanding the Genetic Architecture of Common Cancers.

BACKGROUND: Common cancers develop through a multistep process often including inherited susceptibility. Collaboration among multiple institutions, and funding from multiple sources, has allowed the development of an inexpensive genotyping microarray, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits. METHODS: The OncoArray can be genotyped using a novel technology developed by Illumina to facilitate efficient genotyping. The consortium developed standard approaches for selecting SNPs for study, for quality control of markers, and for ancestry analysis. The array was genotyped at selected sites and with prespecified replicate samples to permit evaluation of genotyping accuracy among centers and by ethnic background. RESULTS: The OncoArray consortium genotyped 447,705 samples. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis. CONCLUSIONS: Results from these analyses will enable researchers to identify new susceptibility loci, perform fine-mapping of new or known loci associated with either single or multiple cancers, assess the degree of overlap in cancer causation and pleiotropic effects of loci that have been identified for disease-specific risk, and jointly model genetic, environmental, and lifestyle-related exposures. IMPACT: Ongoing analyses will shed light on etiology and risk assessment for many types of cancer. Cancer Epidemiol Biomarkers Prev; 26(1); 126-35. ©2016 AACR.

Association of breast cancer risk in BRCA1 and BRCA2 mutation carriers with genetic variants showing differential allelic expression: identification of a modifier of breast cancer risk at locus 11q22.3.

PURPOSE: Cis-acting regulatory SNPs resulting in differential allelic expression (DAE) may, in part, explain the underlying phenotypic variation associated with many complex diseases. To investigate whether common variants associated with DAE were involved in breast cancer susceptibility among BRCA1 and BRCA2 mutation carriers, a list of 175 genes was developed based of their involvement in cancer-related pathways. METHODS: Using data from a genome-wide map of SNPs associated with allelic expression, we assessed the association of ~320 SNPs located in the vicinity of these genes with breast and ovarian cancer risks in 15,252 BRCA1 and 8211 BRCA2 mutation carriers ascertained from 54 studies participating in the Consortium of Investigators of Modifiers of BRCA1/2. RESULTS: We identified a region on 11q22.3 that is significantly associated with breast cancer risk in BRCA1 mutation carriers (most significant SNP rs228595 p = 7 × 10-6). This association was absent in BRCA2 carriers (p = 0.57). The 11q22.3 region notably encompasses genes such as ACAT1, NPAT, and ATM. Expression quantitative trait loci associations were observed in both normal breast and tumors across this region, namely for ACAT1, ATM, and other genes. In silico analysis revealed some overlap between top risk-associated SNPs and relevant biological features in mammary cell data, which suggests potential functional significance. CONCLUSION: We identified 11q22.3 as a new modifier locus in BRCA1 carriers. Replication in larger studies using estrogen receptor (ER)-negative or triple-negative (i.e., ER-, progesterone receptor-, and HER2-negative) cases could therefore be helpful to confirm the association of this locus with breast cancer risk.

Association of breast cancer risk with genetic variants showing differential allelic expression: Identification of a novel breast cancer susceptibility locus at 4q21.

There are significant inter-individual differences in the levels of gene expression. Through modulation of gene expression, cis-acting variants represent an important source of phenotypic variation. Consequently, cis-regulatory SNPs associated with differential allelic expression are functional candidates for further investigation as disease-causing variants. To investigate whether common variants associated with differential allelic expression were involved in breast cancer susceptibility, a list of genes was established on the basis of their involvement in cancer related pathways and/or mechanisms. Thereafter, using data from a genome-wide map of allelic expression associated SNPs, 313 genetic variants were selected and their association with breast cancer risk was then evaluated in 46,451 breast cancer cases and 42,599 controls of European ancestry ascertained from 41 studies participating in the Breast Cancer Association Consortium. The associations were evaluated with overall breast cancer risk and with estrogen receptor negative and positive disease. One novel breast cancer susceptibility locus on 4q21 (rs11099601) was identified (OR = 1.05, P = 5.6x10-6). rs11099601 lies in a 135 kb linkage disequilibrium block containing several genes, including, HELQ, encoding the protein HEL308 a DNA dependant ATPase and DNA Helicase involved in DNA repair, MRPS18C encoding the Mitochondrial Ribosomal Protein S18C and FAM175A (ABRAXAS), encoding a BRCA1 BRCT domain-interacting protein involved in DNA damage response and double-strand break (DSB) repair. Expression QTL analysis in breast cancer tissue showed rs11099601 to be associated with HELQ (P = 8.28x10-14), MRPS18C (P = 1.94x10-27) and FAM175A (P = 3.83x10-3), explaining about 20%, 14% and 1%, respectively of the variance inexpression of these genes in breast carcinomas.

Fine scale mapping of the 17q22 breast cancer locus using dense SNPs, genotyped within the Collaborative Oncological Gene-Environment Study (COGs).

Genome-wide association studies have found SNPs at 17q22 to be associated with breast cancer risk. To identify potential causal variants related to breast cancer risk, we performed a high resolution fine-mapping analysis that involved genotyping 517 SNPs using a custom Illumina iSelect array (iCOGS) followed by imputation of genotypes for 3,134 SNPs in more than 89,000 participants of European ancestry from the Breast Cancer Association Consortium (BCAC). We identified 28 highly correlated common variants, in a 53 Kb region spanning two introns of the STXBP4 gene, that are strong candidates for driving breast cancer risk (lead SNP rs2787486 (OR = 0.92; CI 0.90-0.94; P = 8.96 × 10(-15))) and are correlated with two previously reported risk-associated variants at this locus, SNPs rs6504950 (OR = 0.94, P = 2.04 × 10(-09), r(2) = 0.73 with lead SNP) and rs1156287 (OR = 0.93, P = 3.41 × 10(-11), r(2) = 0.83 with lead SNP). Analyses indicate only one causal SNP in the region and several enhancer elements targeting STXBP4 are located within the 53 kb association signal. Expression studies in breast tumor tissues found SNP rs2787486 to be associated with increased STXBP4 expression, suggesting this may be a target gene of this locus.