Functional Homologous Recombination (HR) Screening Shows the Majority of BRCA1/2-Mutant Breast and Ovarian Cancer Cell Lines Are HR-Proficient.

Tumors with a pathogenic BRCA1/2 mutation are homologous recombination (HR)-deficient (HRD) and consequently sensitive to platinum-based chemotherapy and Poly-[ADP-Ribose]-Polymerase inhibitors (PARPi). We hypothesized that functional HR status better reflects real-time HR status than BRCA1/2 mutation status. Therefore, we determined the functional HR status of 53 breast cancer (BC) and 38 ovarian cancer (OC) cell lines by measuring the formation of RAD51 foci after irradiation. Discrepancies between functional HR and BRCA1/2 mutation status were investigated using exome sequencing, methylation and gene expression data from 50 HR-related genes. A pathogenic BRCA1/2 mutation was found in 10/53 (18.9%) of BC and 7/38 (18.4%) of OC cell lines. Among BRCA1/2-mutant cell lines, 14/17 (82.4%) were HR-proficient (HRP), while 1/74 (1.4%) wild-type cell lines was HRD. For most (80%) cell lines, we explained the discrepancy between functional HR and BRCA1/2 mutation status. Importantly, 12/14 (85.7%) BRCA1/2-mutant HRP cell lines were explained by mechanisms directly acting on BRCA1/2. Finally, functional HR status was strongly associated with COSMIC single base substitution signature 3, but not BRCA1/2 mutation status. Thus, the majority of BRCA1/2-mutant cell lines do not represent a suitable model for HRD. Moreover, exclusively determining BRCA1/2 mutation status may not suffice for platinum-based chemotherapy or PARPi patient selection.

The benefit of adding polygenic risk scores, lifestyle factors, and breast density to family history and genetic status for breast cancer risk and surveillance classification of unaffected women from germline CHEK2 c.1100delC families.

To determine the changes in surveillance category by adding a polygenic risk score based on 311 breast cancer (BC)-associated variants (PRS311), questionnaire-based risk factors and breast density on personalized BC risk in unaffected women from Dutch CHEK2 c.1100delC families. In total, 117 unaffected women (58 heterozygotes and 59 non-carriers) from CHEK2 families were included. Blood-derived DNA samples were genotyped with the GSAMDv3-array to determine PRS311. Lifetime BC risk was calculated in CanRisk, which uses data from the Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm (BOADICEA). Women, were categorized into three surveillance groups. The surveillance advice was reclassified in 37.9 % of heterozygotes and 32.2 % of non-carriers after adding PRS311. Including questionnaire-based risk factors resulted in an additional change in 20.0 % of heterozygotes and 13.2 % of non-carriers; and a subanalysis showed that adding breast density on top shifted another 17.9 % of heterozygotes and 33.3 % of non-carriers. Overall, the majority of heterozygotes were reclassified to a less intensive surveillance, while non-carriers would require intensified surveillance. The addition of PRS311, questionnaire-based risk factors and breast density to family history resulted in a more personalized BC surveillance advice in CHEK2-families, which may lead to more efficient use of surveillance.

Breast cancer genomes from CHEK2 c.1100delC mutation carriers lack somatic TP53 mutations and display a unique structural variant size distribution profile.

BACKGROUND: CHEK2 c.1100delC was the first moderate-risk breast cancer (BC) susceptibility allele discovered. Despite several genomic, transcriptomic and functional studies, however, it is still unclear how exactly CHEK2 c.1100delC promotes tumorigenesis. Since the mutational landscape of a tumor reflects the processes that have operated on its development, the aim of this study was to uncover the somatic genomic landscape of CHEK2-associated BC. METHODS: We sequenced primary BC (pBC) and normal genomes of 20 CHEK2 c.1100delC mutation carriers as well as their pBC transcriptomes. Including pre-existing cohorts, we exhaustively compared CHEK2 pBC genomes to those from BRCA1/2 mutation carriers, those that displayed homologous recombination deficiency (HRD) and ER- and ER+ pBCs, totaling to 574 pBC genomes. Findings were validated in 517 metastatic BC genomes subdivided into the same subgroups. Transcriptome data from 168 ER+ pBCs were used to derive a TP53-mutant gene expression signature and perform cluster analysis with CHEK2 BC transcriptomes. Finally, clinical outcome of CHEK2 c.1100delC carriers was compared with BC patients displaying somatic TP53 mutations in two well-described retrospective cohorts totaling to 942 independent pBC cases. RESULTS: BC genomes from CHEK2 mutation carriers were most similar to ER+ BC genomes and least similar to those of BRCA1/2 mutation carriers in terms of tumor mutational burden as well as mutational signatures. Moreover, CHEK2 BC genomes did not show any evidence of HRD. Somatic TP53 mutation frequency and the size distribution of structural variants (SVs), however, were different compared to ER+ BC. Interestingly, BC genomes with bi-allelic CHEK2 inactivation lacked somatic TP53 mutations and transcriptomic analysis indicated a shared biology with TP53 mutant BC. Moreover, CHEK2 BC genomes had an increased frequency of > 1 Mb deletions, inversions and tandem duplications with peaks at specific sizes. The high chromothripsis frequency among CHEK2 BC genomes appeared, however, not associated with this unique SV size distribution profile. CONCLUSIONS: CHEK2 BC genomes are most similar to ER+ BC genomes, but display unique features that may further unravel CHEK2-driven tumorigenesis. Increased insight into this mechanism could explain the shorter survival of CHEK2 mutation carriers that is likely driven by intrinsic tumor aggressiveness rather than endocrine resistance.

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.

Clinical applicability of the Polygenic Risk Score for breast cancer risk prediction in familial cases.

BACKGROUND: Common low-risk variants are presently not used to guide clinical management of familial breast cancer (BC). We explored the additive impact of a 313-variant-based Polygenic Risk Score (PRS313) relative to standard gene testing in non-BRCA1/2 Dutch BC families. METHODS: We included 3918 BC cases from 3492 Dutch non-BRCA1/2 BC families and 3474 Dutch population controls. The association of the standardised PRS313 with BC was estimated using a logistic regression model, adjusted for pedigree-based family history. Family history of the controls was imputed for this analysis. SEs were corrected to account for relatedness of individuals. Using the BOADICEA (Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm) V.5 model, lifetime risks were retrospectively calculated with and without individual PRS313. For 2586 cases and 2584 controls, the carrier status of pathogenic variants (PVs) in ATM, CHEK2 and PALB2 was known. RESULTS: The family history-adjusted PRS313 was significantly associated with BC (per SD OR=1.97, 95% CI 1.84 to 2.11). Including the PRS313 in BOADICEA family-based risk prediction would have changed screening recommendations in up to 27%, 36% and 34% of cases according to BC screening guidelines from the USA, UK and the Netherlands (National Comprehensive Cancer Network, National Institute for Health and Care Excellence, and Netherlands Comprehensive Cancer Organisation), respectively. For the population controls, without information on family history, this was up to 39%, 44% and 58%, respectively. Among carriers of PVs in known moderate BC susceptibility genes, the PRS313 had the largest impact for CHEK2 and ATM. CONCLUSIONS: Our results support the application of the PRS313 in risk prediction for genetically uninformative BC families and families with a PV in moderate BC risk genes.

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.

Expression and Localization of Ferritin-Heavy Chain Predicts Recurrence for Breast Cancer Patients with a BRCA1/2 Mutation.

The ferritin-heavy chain (FTH1) is the catalytic subunit of the ferroxidase ferritin, which prevents oxidative DNA damage via intracellular iron storage. FTH1 was shown to be a prognostic marker for triple-negative breast cancer (BC) patients and associated with an enrichment of CD8+ effector T cells. However, whether the expression and localization of FTH1 are also associated with clinical outcome in other BC subtypes is unknown. Here, we investigated the association of FTH1 with time to survival in BCs from 222 BRCA1/2 mutation carriers by immunohistochemistry on tissue microarrays. In addition, for 51 of these patients, the association between FTH1 and specific subsets of T cells was evaluated on whole slides using automatic scoring algorithms. We revealed that nuclear FTH1 (nFTH1) expression, in multivariable analyses, was associated with a shorter disease-free (HR = 2.71, 95% CI = 1.49-4.92, p = 0.001) and metastasis-free survival (HR = 3.54, 95% CI = 1.45-8.66, p = 0.006) in patients carrying a BRCA1/2 mutation. However, we found no relation between cytoplasmic FTH1 expression and survival of BRCA1/2 mutation carriers. Moreover, we did not detect an association between FTH1 expression and the amount of CD45+ (p = 0.13), CD8+ (p = 0.18), CD4+ (p = 0.20) or FOXP3+ cells (p = 0.17). Consequently, the mechanism underlying the worse recurrence-free survival of nFTH1 expression in BRCA1/2 mutation carriers needs further investigation.

Progression-free survival and overall survival after BRCA1/2-associated epithelial ovarian cancer: A matched cohort study.

INTRODUCTION: Germline BRCA1/2-associated epithelial ovarian cancer has been associated with better progression-free survival and overall survival than sporadic epithelial ovarian cancer, but conclusive data are lacking. METHODS: We matched 389 BRCA1-associated and 123 BRCA2-associated epithelial ovarian cancer patients 1:1 to sporadic epithelial ovarian cancer patients on year of birth, year of diagnosis, and FIGO stage (< = IIA/> = IIB). Germline DNA test was performed before or after epithelial ovarian cancer diagnosis. All patients received chemotherapy. We used Cox proportional hazards models to estimate the associations between mutation status (BRCA1 or BRCA2 versus sporadic) and progression-free survival and overall survival. To investigate whether DNA testing after epithelial ovarian cancer diagnosis resulted in survival bias, we performed additional analyses limited to BRCA1/2-associated epithelial ovarian cancer patients with a DNA test result before cancer diagnosis (n = 73 BRCA1; n = 9 BRCA2) and their matched sporadic controls. RESULTS: The median follow-up was 4.4 years (range 0.1-30.1). During the first three years after epithelial ovarian cancer diagnosis, progression-free survival was better for BRCA1 (HR 0.88, 95% CI 0.74-1.04) and BRCA2 (HR 0.58, 95% CI 0.41-0.81) patients than for sporadic patients. Overall survival was better during the first six years after epithelial ovarian cancer for BRCA1 (HR 0.7, 95% CI 0.58-0.84) and BRCA2 (HR 0.41, 95% CI 0.29-0.59) patients. After surviving these years, survival benefits disappeared or were in favor of the sporadic patients. CONCLUSION: For epithelial ovarian cancer patients who received chemotherapy, we confirmed survival benefit for BRCA1 and BRCA2 germline pathogenic variant carriers. This may indicate higher sensitivity to chemotherapy, both in first line treatment and in the recurrent setting. The observed benefit appears to be limited to a relatively short period after epithelial ovarian cancer diagnosis.

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.

Breast cancer risks associated with missense variants in breast cancer susceptibility genes.

BACKGROUND: Protein truncating variants in ATM, BRCA1, BRCA2, CHEK2, and PALB2 are associated with increased breast cancer risk, but risks associated with missense variants in these genes are uncertain. METHODS: We analyzed data on 59,639 breast cancer cases and 53,165 controls from studies participating in the Breast Cancer Association Consortium BRIDGES project. We sampled training (80%) and validation (20%) sets to analyze rare missense variants in ATM (1146 training variants), BRCA1 (644), BRCA2 (1425), CHEK2 (325), and PALB2 (472). We evaluated breast cancer risks according to five in silico prediction-of-deleteriousness algorithms, functional protein domain, and frequency, using logistic regression models and also mixture models in which a subset of variants was assumed to be risk-associated. RESULTS: The most predictive in silico algorithms were Helix (BRCA1, BRCA2 and CHEK2) and CADD (ATM). Increased risks appeared restricted to functional protein domains for ATM (FAT and PIK domains) and BRCA1 (RING and BRCT domains). For ATM, BRCA1, and BRCA2, data were compatible with small subsets (approximately 7%, 2%, and 0.6%, respectively) of rare missense variants giving similar risk to those of protein truncating variants in the same gene. For CHEK2, data were more consistent with a large fraction (approximately 60%) of rare missense variants giving a lower risk (OR 1.75, 95% CI (1.47-2.08)) than CHEK2 protein truncating variants. There was little evidence for an association with risk for missense variants in PALB2. The best fitting models were well calibrated in the validation set. CONCLUSIONS: These results will inform risk prediction models and the selection of candidate variants for functional assays and could contribute to the clinical reporting of gene panel testing for breast cancer susceptibility.

Identifying Transcripts with Tandem Duplications from RNA-Sequencing Data to Predict BRCA1-Type Primary Breast Cancer.

Patients with cancers that are deficient for homologous recombination repair (HRD) may benefit from PARP inhibitor treatment. Therefore, methods that identify such cancers are crucial. Using whole genome sequencing data, specific genomic scars derived from somatic mutations and genomic rearrangements can identify HRD tumors, with only BRCA1-like HRD cancers profoundly displaying small (<10 kb) tandem duplications (TDs). In this manuscript we describe a method of detecting BRCA1-type HRD in breast cancer (BC) solely from RNA sequencing data by identifying TDs surfacing in transcribed genes. We find that the number of identified TDs (TD-score) is significantly higher in BRCA1-type vs. BRCA2-type BCs, or vs. HR-proficient BCs (p = 2.4 × 10-6 and p = 2.7 × 10-12, respectively). A TD-score ≥2 shows an 88.2% sensitivity (30 out of 34) to detect a BRCA1-type BC, with a specificity of 64.7% (143 out of 221). Pathway enrichment analyses showed genes implicated in cancer to be affected by TDs of which PTEN was found significantly more frequently affected by a TD in BRCA1-type BC. In conclusion, we here describe a novel method to identify TDs in transcripts and classify BRCA1-type BCs with high sensitivity.

A pipeline for copy number profiling of single circulating tumour cells to assess intrapatient tumour heterogeneity.

Intrapatient tumour heterogeneity is likely a major determinant of clinical outcome in cancer patients. To assess heterogeneity in a minimally invasive manner, methods to perform single circulating tumour cell (CTC) genomics at high resolution are necessary. However, due to the rarity of CTCs, development of such methods is challenging. Here, we developed a modular single CTC analysis pipeline to assess intrapatient heterogeneity by copy number (CN) profiling. To optimize this pipeline, spike-in experiments using MCF-7 breast cancer cells were performed. The VyCAP puncher system was used to isolate single cells. The quality of whole genome amplification (WGA) products generated by REPLI-g and Ampli1™ methods, as well as the results from the Illumina Truseq and the Ampli1™ LowPass library preparation techniques, was compared. Moreover, a bioinformatic pipeline was designed to generate CN profiles from single CTCs. The optimal combination of Ampli1™ WGA and Illumina Truseq library preparation was successfully validated on patient-derived CTCs. In conclusion, we developed a novel modular pipeline to isolate single CTCs and subsequently generate detailed patient-derived CN profiles that allow assessment of intrapatient heterogeneity in future studies.

Germline variants and breast cancer survival in patients with distant metastases at primary breast cancer diagnosis.

Breast cancer metastasis accounts for most of the deaths from breast cancer. Identification of germline variants associated with survival in aggressive types of breast cancer may inform understanding of breast cancer progression and assist treatment. In this analysis, we studied the associations between germline variants and breast cancer survival for patients with distant metastases at primary breast cancer diagnosis. We used data from the Breast Cancer Association Consortium (BCAC) including 1062 women of European ancestry with metastatic breast cancer, 606 of whom died of breast cancer. We identified two germline variants on chromosome 1, rs138569520 and rs146023652, significantly associated with breast cancer-specific survival (P = 3.19 × 10-8 and 4.42 × 10-8). In silico analysis suggested a potential regulatory effect of the variants on the nearby target genes SDE2 and H3F3A. However, the variants showed no evidence of association in a smaller replication dataset. The validation dataset was obtained from the SNPs to Risk of Metastasis (StoRM) study and included 293 patients with metastatic primary breast cancer at diagnosis. Ultimately, larger replication studies are needed to confirm the identified associations.

Common Susceptibility Loci for Male Breast Cancer.

BACKGROUND: The etiology of male breast cancer (MBC) is poorly understood. In particular, the extent to which the genetic basis of MBC differs from female breast cancer (FBC) is unknown. A previous genome-wide association study of MBC identified 2 predisposition loci for the disease, both of which were also associated with risk of FBC. METHODS: We performed genome-wide single nucleotide polymorphism genotyping of European ancestry MBC case subjects and controls in 3 stages. Associations between directly genotyped and imputed single nucleotide polymorphisms with MBC were assessed using fixed-effects meta-analysis of 1380 cases and 3620 controls. Replication genotyping of 810 cases and 1026 controls was used to validate variants with P values less than 1 × 10-06. Genetic correlation with FBC was evaluated using linkage disequilibrium score regression, by comprehensively examining the associations of published FBC risk loci with risk of MBC and by assessing associations between a FBC polygenic risk score and MBC. All statistical tests were 2-sided. RESULTS: The genome-wide association study identified 3 novel MBC susceptibility loci that attained genome-wide statistical significance (P < 5 × 10-08). Genetic correlation analysis revealed a strong shared genetic basis with estrogen receptor-positive FBC. Men in the top quintile of genetic risk had a fourfold increased risk of breast cancer relative to those in the bottom quintile (odds ratio = 3.86, 95% confidence interval = 3.07 to 4.87, P = 2.08 × 10-30). CONCLUSIONS: These findings advance our understanding of the genetic basis of MBC, providing support for an overlapping genetic etiology with FBC and identifying a fourfold high-risk group of susceptible men.

Genetically Predicted Levels of DNA Methylation Biomarkers and Breast Cancer Risk: Data From 228 951 Women of European Descent.

BACKGROUND: DNA methylation plays a critical role in breast cancer development. Previous studies have identified DNA methylation marks in white blood cells as promising biomarkers for breast cancer. However, these studies were limited by low statistical power and potential biases. Using a new methodology, we investigated DNA methylation marks for their associations with breast cancer risk. METHODS: Statistical models were built to predict levels of DNA methylation marks using genetic data and DNA methylation data from HumanMethylation450 BeadChip from the Framingham Heart Study (n = 1595). The prediction models were validated using data from the Women's Health Initiative (n = 883). We applied these models to genomewide association study (GWAS) data of 122 977 breast cancer patients and 105 974 controls to evaluate if the genetically predicted DNA methylation levels at CpG sites (CpGs) are associated with breast cancer risk. All statistical tests were two-sided. RESULTS: Of the 62 938 CpG sites CpGs investigated, statistically significant associations with breast cancer risk were observed for 450 CpGs at a Bonferroni-corrected threshold of P less than 7.94 × 10-7, including 45 CpGs residing in 18 genomic regions, that have not previously been associated with breast cancer risk. Of the remaining 405 CpGs located within 500 kilobase flaking regions of 70 GWAS-identified breast cancer risk variants, the associations for 11 CpGs were independent of GWAS-identified variants. Integrative analyses of genetic, DNA methylation, and gene expression data found that 38 CpGs may affect breast cancer risk through regulating expression of 21 genes. CONCLUSION: Our new methodology can identify novel DNA methylation biomarkers for breast cancer risk and can be applied to other diseases.

Addition of a 161-SNP polygenic risk score to family history-based risk prediction: impact on clinical management in non-BRCA1/2 breast cancer families.

BACKGROUND: The currently known breast cancer-associated single nucleotide polymorphisms (SNPs) are presently not used to guide clinical management. We explored whether a genetic test that incorporates a SNP-based polygenic risk score (PRS) is clinically meaningful in non-BRCA1/2 high-risk breast cancer families. METHODS: 101 non-BRCA1/2 high-risk breast cancer families were included; 323 cases and 262 unaffected female relatives were genotyped. The 161-SNP PRS was calculated and standardised to 327 population controls (sPRS). Association analysis was performed using a Cox-type random effect regression model adjusted by family history. Updated individualised breast cancer lifetime risk scores were derived by combining the Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm breast cancer lifetime risk with the effect of the sPRS. RESULTS: The mean sPRS for cases and their unaffected relatives was 0.70 (SD=0.9) and 0.53 (SD=0.9), respectively. A significant association was found between sPRS and breast cancer, HR=1.16, 95% CI 1.03 to 1.28, p=0.026. Addition of the sPRS to risk prediction based on family history alone changed screening recommendations in 11.5%, 14.7% and 19.8 % of the women according to breast screening guidelines from the USA (National Comprehensive Cancer Network), UK (National Institute for Health and Care Excellence and the Netherlands (Netherlands Comprehensive Cancer Organisation), respectively. CONCLUSION: Our results support the application of the PRS in risk prediction and clinical management of women from genetically unexplained breast cancer families.

Elucidating the Underlying Functional Mechanisms of Breast Cancer Susceptibility Through Post-GWAS Analyses.

Genome-wide association studies (GWAS) have identified more than 170 single nucleotide polymorphisms (SNPs) associated with the susceptibility to breast cancer. Together, these SNPs explain 18% of the familial relative risk, which is estimated to be nearly half of the total familial breast cancer risk that is collectively explained by low-risk susceptibility alleles. An important aspect of this success has been the access to large sample sizes through collaborative efforts within the Breast Cancer Association Consortium (BCAC), but also collaborations between cancer association consortia. Despite these achievements, however, understanding of each variant's underlying mechanism and how these SNPs predispose women to breast cancer remains limited and represents a major challenge in the field, particularly since the vast majority of the GWAS-identified SNPs are located in non-coding regions of the genome and are merely tags for the causal variants. In recent years, fine-scale mapping studies followed by functional evaluation of putative causal variants have begun to elucidate the biological function of several GWAS-identified variants. In this review, we discuss the findings and lessons learned from these post-GWAS analyses of 22 risk loci. Identifying the true causal variants underlying breast cancer susceptibility and their function not only provides better estimates of the explained familial relative risk thereby improving polygenetic risk scores (PRSs), it also increases our understanding of the biological mechanisms responsible for causing susceptibility to breast cancer. This will facilitate the identification of further breast cancer risk alleles and the development of preventive medicine for those women at increased risk for developing the disease.

The Prevalence of CD146 Expression in Breast Cancer Subtypes and Its Relation to Outcome.

CD146, involved in epithelial-to-mesenchymal transition (EMT), might affect cancer aggressiveness. We here investigated the prevalence of CD146 expression in breast cancer subtypes, its relation to prognosis, the relation between CD146 and EMT and the outcome to tamoxifen. Primary breast cancer tissues from 1342 patients were available for this retrospective study and immunohistochemically stained for CD146. For survival analyses, pure prognosis was studied by only including lymph-node negative patients who did not receive (neo)adjuvant systemic treatment (n = 551). 11% of the tumors showed CD146 expression. CD146 expression was most prevalent in triple-negative cases (64%, p < 0.001). In univariable analysis, CD146 expression was a prognostic factor for both metastasis-free survival (MFS) (p = 0.020) and overall survival (OS) (p = 0.037), but not in multivariable analysis (including age, tumor size, grade, estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and Ki-67). No correlation between CD146 and EMT nor difference in outcome to first-line tamoxifen was seen. In this large series, our data showed that CD146 is present in primary breast cancer and is a pure prognostic factor for MFS and OS in breast cancer patients. We did not see an association between CD146 expression and EMT nor on outcome to tamoxifen.

Etiology of hormone receptor positive breast cancer differs by levels of histologic grade and proliferation.

Limited epidemiological evidence suggests that the etiology of hormone receptor positive (HR+) breast cancer may differ by levels of histologic grade and proliferation. We pooled risk factor and pathology data on 5,905 HR+ breast cancer cases and 26,281 controls from 11 epidemiological studies. Proliferation was determined by centralized automated measures of KI67 in tissue microarrays. Odds ratios (OR), 95% confidence intervals (CI) and p-values for case-case and case-control comparisons for risk factors in relation to levels of grade and quartiles (Q1-Q4) of KI67 were estimated using polytomous logistic regression models. Case-case comparisons showed associations between nulliparity and high KI67 [OR (95% CI) for Q4 vs. Q1 = 1.54 (1.22, 1.95)]; obesity and high grade [grade 3 vs. 1 = 1.68 (1.31, 2.16)] and current use of combined hormone therapy (HT) and low grade [grade 3 vs. 1 = 0.27 (0.16, 0.44)] tumors. In case-control comparisons, nulliparity was associated with elevated risk of tumors with high but not low levels of proliferation [1.43 (1.14, 1.81) for KI67 Q4 vs. 0.83 (0.60, 1.14) for KI67 Q1]; obesity among women ≥50 years with high but not low grade tumors [1.55 (1.17, 2.06) for grade 3 vs. 0.88 (0.66, 1.16) for grade 1] and HT with low but not high grade tumors [3.07 (2.22, 4.23) for grade 1 vs. 0.85 (0.55, 1.30) for grade 3]. Menarcheal age and family history were similarly associated with HR+ tumors of different grade or KI67 levels. These findings provide insights into the etiologic heterogeneity of HR+ tumors.

Reproductive profiles and risk of breast cancer subtypes: a multi-center case-only study.

BACKGROUND: Previous studies have shown that reproductive factors are differentially associated with breast cancer (BC) risk by subtypes. The aim of this study was to investigate associations between reproductive factors and BC subtypes, and whether these vary by age at diagnosis. METHODS: We used pooled data on tumor markers (estrogen and progesterone receptor, human epidermal growth factor receptor-2 (HER2)) and reproductive risk factors (parity, age at first full-time pregnancy (FFTP) and age at menarche) from 28,095 patients with invasive BC from 34 studies participating in the Breast Cancer Association Consortium (BCAC). In a case-only analysis, we used logistic regression to assess associations between reproductive factors and BC subtype compared to luminal A tumors as a reference. The interaction between age and parity in BC subtype risk was also tested, across all ages and, because age was modeled non-linearly, specifically at ages 35, 55 and 75 years. RESULTS: Parous women were more likely to be diagnosed with triple negative BC (TNBC) than with luminal A BC, irrespective of age (OR for parity = 1.38, 95% CI 1.16-1.65, p = 0.0004; p for interaction with age = 0.076). Parous women were also more likely to be diagnosed with luminal and non-luminal HER2-like BCs and this effect was slightly more pronounced at an early age (p for interaction with age = 0.037 and 0.030, respectively). For instance, women diagnosed at age 35 were 1.48 (CI 1.01-2.16) more likely to have luminal HER2-like BC than luminal A BC, while this association was not significant at age 75 (OR = 0.72, CI 0.45-1.14). While age at menarche was not significantly associated with BC subtype, increasing age at FFTP was non-linearly associated with TNBC relative to luminal A BC. An age at FFTP of 25 versus 20 years lowered the risk for TNBC (OR = 0.78, CI 0.70-0.88, p < 0.0001), but this effect was not apparent at a later FFTP. CONCLUSIONS: Our main findings suggest that parity is associated with TNBC across all ages at BC diagnosis, whereas the association with luminal HER2-like BC was present only for early onset BC.