Background Due to the abundant usage of chemotherapy in young triple-negative breast cancer (TNBC) patients, the unbiased prognostic value of BRCA1-related biomarkers in this population remains... Show moreBackground Due to the abundant usage of chemotherapy in young triple-negative breast cancer (TNBC) patients, the unbiased prognostic value of BRCA1-related biomarkers in this population remains unclear. In addition, whether BRCA1-related biomarkers modify the well-established prognostic value of stromal tumor-infiltrating lymphocytes (sTILs) is unknown. This study aimed to compare the outcomes of young, node-negative, chemotherapy-na & iuml;ve TNBC patients according to BRCA1 status, taking sTILs into account. Methods We included 485 Dutch women diagnosed with node-negative TNBC under age 40 between 1989 and 2000. During this period, these women were considered low-risk and did not receive chemotherapy. BRCA1 status, including pathogenic germline BRCA1 mutation (gBRCA1m), somatic BRCA1 mutation (sBRCA1m), and tumor BRCA1 promoter methylation (BRCA1-PM), was assessed using DNA from formalin-fixed paraffin-embedded tissue. sTILs were assessed according to the international guideline. Patients' outcomes were compared using Cox regression and competing risk models. Results Among the 399 patients with BRCA1 status, 26.3% had a gBRCA1m, 5.3% had a sBRCA1m, 36.6% had tumor BRCA1-PM, and 31.8% had BRCA1-non-altered tumors. Compared to BRCA1-non-alteration, gBRCA1m was associated with worse overall survival (OS) from the fourth year after diagnosis (adjusted HR, 2.11; 95% CI, 1.18-3.75), and this association attenuated after adjustment for second primary tumors. Every 10% sTIL increment was associated with 16% higher OS (adjusted HR, 0.84; 95% CI, 0.78-0.90) in gBRCA1m, sBRCA1m, or BRCA1-non-altered patients and 31% higher OS in tumor BRCA1-PM patients. Among the 66 patients with tumor BRCA1-PM and >= 50% sTILs, we observed excellent 15-year OS (97.0%; 95% CI, 92.9-100%). Conversely, among the 61 patients with gBRCA1m and < 50% sTILs, we observed poor 15-year OS (50.8%; 95% CI, 39.7-65.0%). Furthermore, gBRCA1m was associated with higher (adjusted subdistribution HR, 4.04; 95% CI, 2.29-7.13) and tumor BRCA1-PM with lower (adjusted subdistribution HR, 0.42; 95% CI, 0.19-0.95) incidence of second primary tumors, compared to BRCA1-non-alteration. Conclusions Although both gBRCA1m and tumor BRCA1-PM alter BRCA1 gene transcription, they are associated with different outcomes in young, node-negative, chemotherapy-na & iuml;ve TNBC patients. By combining sTILs and BRCA1 status for risk classification, we were able to identify potential subgroups in this population to intensify and optimize adjuvant treatment. Show less
Wilcox, N.; Dumont, M.; González-Neira, A.; Carvalho, S.; Beauparlant, C.J.; Crotti, M.; ... ; Simard, J. 2023
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... Show moreLinkage 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 x 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 x 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.Meta-analysis of three large whole-exome sequencing datasets highlights protein-truncating and rare missense variants associated with breast cancer susceptibility. Show less
Smid, M.; Schmidt, M.K.; Smissen, W.J.C. van der; Ruigrok-Ritstier, K.; Schreurs, M.A.C.; Cornelissen, S.; ... ; Hollestelle, A. 2023
BackgroundCHEK2 c.1100delC was the first moderate-risk breast cancer (BC) susceptibility allele discovered. Despite several genomic, transcriptomic and functional studies, however, it is still... Show moreBackgroundCHEK2 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.MethodsWe 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.ResultsBC 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.ConclusionsCHEK2 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. Show less
Neoadjuvant ipilimumab + nivolumab has demonstrated high pathologic response rates in stage III melanoma. Patients with low intra-tumoral interferon-γ (IFN-γ) signatures are less likely to benefit... Show moreNeoadjuvant ipilimumab + nivolumab has demonstrated high pathologic response rates in stage III melanoma. Patients with low intra-tumoral interferon-γ (IFN-γ) signatures are less likely to benefit. We show that domatinostat (a class I histone deacetylase inhibitor) addition to anti-PD-1 + anti-CTLA-4 increased the IFN-γ response and reduced tumor growth in our murine melanoma model, rationalizing evaluation in patients. To stratify patients into IFN-γ high and low cohorts, we developed a baseline IFN-γ signature expression algorithm, which was prospectively tested in the DONIMI trial. Patients with stage III melanoma and high intra-tumoral IFN-γ scores were randomized to neoadjuvant nivolumab or nivolumab + domatinostat, while patients with low IFN-γ scores received nivolumab + domatinostat or ipilimumab + nivolumab + domatinostat. Domatinostat addition to neoadjuvant nivolumab ± ipilimumab did not delay surgery but induced unexpected severe skin toxicity, hampering domatinostat dose escalation. At studied dose levels, domatinostat addition did not increase treatment efficacy. The baseline IFN-γ score adequately differentiated patients who were likely to benefit from nivolumab alone versus patients who require other therapies. Show less
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... Show moreBackground 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 x 10(-6)) and AC058822.1 (P = 1.47 x 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 x 10(-5)), demonstrating the importance of diversifying study cohorts. Show less
Jong, V.M.T. de; Wang, Y.W.; Hoeve, N.D. ter; Opdam, M.; Stathonikos, N.; Józwiak, K.; ... ; Linn, S.C. 2022
PURPOSE Triple-negative breast cancer (TNBC) is considered aggressive, and therefore, virtually all young patients with TNBC receive (neo)adjuvant chemotherapy. Increased stromal tumor-infiltrating... Show morePURPOSE Triple-negative breast cancer (TNBC) is considered aggressive, and therefore, virtually all young patients with TNBC receive (neo)adjuvant chemotherapy. Increased stromal tumor-infiltrating lymphocytes (sTILs) have been associated with a favorable prognosis in TNBC. However, whether this association holds for patients who are node-negative (N0), young (< 40 years), and chemotherapy-naive, and thus can be used for chemotherapy de-escalation strategies, is unknown.METHODS We selected all patients with N0 TNBC diagnosed between 1989 and 2000 from a Dutch population-based registry. Patients were age < 40 years at diagnosis and had not received (neo)adjuvant systemic therapy, as was standard practice at the time. Formalin-fixed paraffin-embedded blocks were retrieved (PALGA: Dutch Pathology Registry), and a pathology review including sTILs was performed. Patients were categorized according to sTILs (< 30%, 30%-75%, and >= 75%). Multivariable Cox regression was performed for overall survival, with or without sTILs as a covariate. Cumulative incidence of distant metastasis or death was analyzed in a competing risk model, with second primary tumors as competing risk.RESULTS sTILs were scored for 441 patients. High sTILs (>= 75%; 21%) translated into an excellent prognosis with a 15-year cumulative incidence of a distant metastasis or death of only 2.1% (95% CI, 0 to 5.0), whereas low sTILs (< 30%; 52%) had an unfavorable prognosis with a 15-year cumulative incidence of a distant metastasis or death of 38.4% (32.1 to 44.6). In addition, every 10% increment of sTILs decreased the risk of death by 19% (adjusted hazard ratio: 0.81; 95% CI, 0.76 to 0.87), which are an independent predictor adding prognostic information to standard clinicopathologic variables (chi(2) = 46.7, P < .001).CONCLUSION Chemotherapy-naive, young patients with N0 TNBC with high sTILs (>= 75%) have an excellent long-term prognosis. Therefore, sTILs should be considered for prospective clinical trials investigating (neo)adjuvant chemotherapy de-escalation strategies. (c) 2022 by American Society of Clinical Oncology Show less
Dumont, M.; Weber-Lassalle, N.; Joly-Beauparlant, C.; Ernst, C.; Droit, A.; Feng, B.J.; ... ; Simard, J. 2022
Simple Summary Genetic variants explaining approximately 40% of familial breast cancer risk have been identified, thus leaving a significant fraction of the heritability of this disease still... Show moreSimple Summary Genetic variants explaining approximately 40% of familial breast cancer risk have been identified, thus leaving a significant fraction of the heritability of this disease still unexplained. The exact nature of this missing fraction is unknown; more extensive sequencing efforts could potentially identify new moderate-penetrance breast cancer risk alleles. The aim of this study was to perform a large-scale whole-exome sequencing study, followed by a targeted validation, in breast cancer patients and healthy women of European descent. We identified 20 novel genes with modest evidence of association (p-value < 0.05) for either overall or subtype-specific breast cancer; however, much larger studies are needed to confirm the exact role of these genes in susceptibility to breast cancer. 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. Show less
Neoadjuvant ipilimumab plus nivolumab showed high pathologic response rates (pRRs) in patients with macroscopic stage III melanoma in the phase 1b OpACIN () and phase 2 OpACIN-neo () studies(1,2).... Show moreNeoadjuvant ipilimumab plus nivolumab showed high pathologic response rates (pRRs) in patients with macroscopic stage III melanoma in the phase 1b OpACIN () and phase 2 OpACIN-neo () studies(1,2). While the results are promising, data on the durability of these pathologic responses and baseline biomarkers for response and survival were lacking. After a median follow-up of 4 years, none of the patients with a pathologic response (n = 7/9 patients) in the OpACIN study had relapsed. In OpACIN-neo (n = 86), the 2-year estimated relapse-free survival was 84% for all patients, 97% for patients achieving a pathologic response and 36% for nonresponders (P < 0.001). High tumor mutational burden (TMB) and high interferon-gamma-related gene expression signature score (IFN-gamma score) were associated with pathologic response and low risk of relapse; pRR was 100% in patients with high IFN-gamma score/high TMB; patients with high IFN-gamma score/low TMB or low IFN-gamma score/high TMB had pRRs of 91% and 88%; while patients with low IFN-gamma score/low TMB had a pRR of only 39%. These data demonstrate long-term benefit in patients with a pathologic response and show the predictive potential of TMB and IFN-gamma score. Our findings provide a strong rationale for a randomized phase 3 study comparing neoadjuvant ipilimumab plus nivolumab versus standard adjuvant therapy with antibodies against the programmed cell death protein-1 (anti-PD-1) in macroscopic stage III melanoma. Show less
BACKGROUNDGenetic testing for breast cancer susceptibility is widely used, but for many genes, evidence of an association with breast cancer is weak, underlying risk estimates are imprecise, and... Show moreBACKGROUNDGenetic testing for breast cancer susceptibility is widely used, but for many genes, evidence of an association with breast cancer is weak, underlying risk estimates are imprecise, and reliable subtype-specific risk estimates are lacking.METHODSWe used a panel of 34 putative susceptibility genes to perform sequencing on samples from 60,466 women with breast cancer and 53,461 controls. In separate analyses for protein-truncating variants and rare missense variants in these genes, we estimated odds ratios for breast cancer overall and tumor subtypes. We evaluated missense-variant associations according to domain and classification of pathogenicity.RESULTSProtein-truncating variants in 5 genes (ATM, BRCA1, BRCA2, CHEK2, and PALB2) were associated with a risk of breast cancer overall with a P value of less than 0.0001. Protein-truncating variants in 4 other genes (BARD1, RAD51C, RAD51D, and TP53) were associated with a risk of breast cancer overall with a P value of less than 0.05 and a Bayesian false-discovery probability of less than 0.05. For protein-truncating variants in 19 of the remaining 25 genes, the upper limit of the 95% confidence interval of the odds ratio for breast cancer overall was less than 2.0. For protein-truncating variants in ATM and CHEK2, odds ratios were higher for estrogen receptor (ER)-positive disease than for ER-negative disease; for protein-truncating variants in BARD1, BRCA1, BRCA2, PALB2, RAD51C, and RAD51D, odds ratios were higher for ER-negative disease than for ER-positive disease. Rare missense variants (in aggregate) in ATM, CHEK2, and TP53 were associated with a risk of breast cancer overall with a P value of less than 0.001. For BRCA1, BRCA2, and TP53, missense variants (in aggregate) that would be classified as pathogenic according to standard criteria were associated with a risk of breast cancer overall, with the risk being similar to that of protein-truncating variants.CONCLUSIONSThe results of this study define the genes that are most clinically useful for inclusion on panels for the prediction of breast cancer risk, as well as provide estimates of the risks associated with protein-truncating variants, to guide genetic counseling. (Funded by European Union Horizon 2020 programs and others.) Show less
Genome-wide analysis identifies 32 loci associated with breast cancer susceptibility, accounting for estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and... Show moreGenome-wide analysis identifies 32 loci associated with breast cancer susceptibility, accounting for estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade.Breast cancer susceptibility variants frequently show heterogeneity in associations by tumor subtype(1-3). To identify novel loci, we performed a genome-wide association study including 133,384 breast cancer cases and 113,789 controls, plus 18,908 BRCA1 mutation carriers (9,414 with breast cancer) of European ancestry, using both standard and novel methodologies that account for underlying tumor heterogeneity by estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade. We identified 32 novel susceptibility loci (P < 5.0 x 10(-8)), 15 of which showed evidence for associations with at least one tumor feature (false discovery rate < 0.05). Five loci showed associations (P < 0.05) in opposite directions between luminal and non-luminal subtypes. In silico analyses showed that these five loci contained cell-specific enhancers that differed between normal luminal and basal mammary cells. The genetic correlations between five intrinsic-like subtypes ranged from 0.35 to 0.80. The proportion of genome-wide chip heritability explained by all known susceptibility loci was 54.2% for luminal A-like disease and 37.6% for triple-negative disease. The odds ratios of polygenic risk scores, which included 330 variants, for the highest 1% of quantiles compared with middle quantiles were 5.63 and 3.02 for luminal A-like and triple-negative disease, respectively. These findings provide an improved understanding of genetic predisposition to breast cancer subtypes and will inform the development of subtype-specific polygenic risk scores. Show less
Fine-mapping of causal variants and integration of epigenetic and chromatin conformation data identify likely target genes for 150 breast cancer risk regions.Genome-wide association studies have... Show moreFine-mapping of causal variants and integration of epigenetic and chromatin conformation data identify likely target genes for 150 breast cancer risk regions.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. Show less
Yeo, L.L.L.; Andersson, T.; Holmberg, A.; Mpotsaris, A.; Soderman, M.; Holmin, S.; ... ; Cornelissen, S. 2019
Stratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for... Show moreStratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for prediction of estrogen receptor (ER)-specific disease, from the largest available genome-wide association dataset and to empirically validate the PRSs in prospective studies. The development dataset comprised 94,075 case subjects and 75,017 control subjects of European ancestry from 69 studies, divided into training and validation sets. Samples were genotyped using genome-wide arrays, and single-nucleotide polymorphisms (SNPs) were selected by stepwise regression or lasso penalized regression. The best performing PRSs were validated in an independent test set comprising 11,428 case subjects and 18,323 control subjects from 10 prospective studies and 190,040 women from UK Biobank (3,215 incident breast cancers). For the best PRSs (313 SNPs), the odds ratio for overall disease per 1 standard deviation in ten prospective studies was 1.61 (95%CI: 1.57-1.65) with area under receiver-operator curve (AUC) = 0.630 (95%CI: 0.628-0.651). The lifetime risk of overall breast cancer in the top centile of the PRSs was 32.6%. Compared with women in the middle quintile, those in the highest 1% of risk had 4.37- and 2.78-fold risks, and those in the lowest 1% of risk had 0.16- and 0.27-fold risks, of developing ER-positive and ER-negative disease, respectively. Goodness-of-fit tests indicated that this PRS was well calibrated and predicts disease risk accurately in the tails of the distribution. This PRS is a powerful and reliable predictor of breast cancer risk that may improve breast cancer prevention programs. Show less