The life history of 21 breast cancers.

Cancer evolves dynamically as clonal expansions supersede one another driven by shifting selective pressures, mutational processes, and disrupted cancer genes. These processes mark the genome, such that a cancer's life history is encrypted in the somatic mutations present. We developed algorithms to decipher this narrative and applied them to 21 breast cancers. Mutational processes evolve across a cancer's lifespan, with many emerging late but contributing extensive genetic variation. Subclonal diversification is prominent, and most mutations are found in just a fraction of tumor cells. Every tumor has a dominant subclonal lineage, representing more than 50% of tumor cells. Minimal expansion of these subclones occurs until many hundreds to thousands of mutations have accumulated, implying the existence of long-lived, quiescent cell lineages capable of substantial proliferation upon acquisition of enabling genomic changes. Expansion of the dominant subclone to an appreciable mass may therefore represent the final rate-limiting step in a breast cancer's development, triggering diagnosis.

Mutational processes molding the genomes of 21 breast cancers.

All cancers carry somatic mutations. The patterns of mutation in cancer genomes reflect the DNA damage and repair processes to which cancer cells and their precursors have been exposed. To explore these mechanisms further, we generated catalogs of somatic mutation from 21 breast cancers and applied mathematical methods to extract mutational signatures of the underlying processes. Multiple distinct single- and double-nucleotide substitution signatures were discernible. Cancers with BRCA1 or BRCA2 mutations exhibited a characteristic combination of substitution mutation signatures and a distinctive profile of deletions. Complex relationships between somatic mutation prevalence and transcription were detected. A remarkable phenomenon of localized hypermutation, termed "kataegis," was observed. Regions of kataegis differed between cancers but usually colocalized with somatic rearrangements. Base substitutions in these regions were almost exclusively of cytosine at TpC dinucleotides. The mechanisms underlying most of these mutational signatures are unknown. However, a role for the APOBEC family of cytidine deaminases is proposed.

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.

Matched analysis of circulating selenium with the breast cancer selenotranscriptome: a multicentre prospective study.

INTRODUCTION: Low serum selenium and altered tumour RNA expression of certain selenoproteins are associated with a poor breast cancer prognosis. Selenoprotein expression stringently depends on selenium availability, hence circulating selenium may interact with tumour selenoprotein expression. However, there is no matched analysis to date. METHODS: This study included 1453 patients with newly diagnosed breast cancer from the multicentric prospective Sweden Cancerome Analysis Network - Breast study. Total serum selenium, selenoprotein P and glutathione peroxidase 3 were analysed at time of diagnosis. Bulk RNA-sequencing was conducted in matched tumour tissues. Fully adjusted Cox regression models with an interaction term were employed to detect dose-dependent interactions of circulating selenium with the associations of tumour selenoprotein mRNA expression and mortality. RESULTS: 237 deaths were recorded within ~ 9 years follow-up. All three serum selenium biomarkers correlated positively (p < 0.001). All selenoproteins except for GPX6 were expressed in tumour tissues. Single cell RNA-sequencing revealed a heterogeneous expression pattern in the tumour microenvironment. Circulating selenium correlated positively with tumour SELENOW and SELENON expression (p < 0.001). In fully adjusted models, the associations of DIO1, DIO3 and SELENOM with mortality were dose-dependently modified by serum selenium (p < 0.001, p = 0.020, p = 0.038, respectively). With increasing selenium, DIO1 and SELENOM associated with lower, whereas DIO3 expression associated with higher mortality. Association of DIO1 with lower mortality was only apparent in patients with high selenium [above median (70.36 µg/L)], and the HR (95%CI) for one-unit increase in log(FPKM + 1) was 0.70 (0.50-0.98). CONCLUSIONS: This first unbiased analysis of serum selenium with the breast cancer selenotranscriptome identified an effect-modification of selenium on the associations of DIO1, SELENOM, and DIO3 with prognosis. Selenium substitution in patients with DIO1-expressing tumours merits consideration to improve survival.

Extended genetic analysis and tumor characteristics in over 4600 women with suspected hereditary breast and ovarian cancer.

BACKGROUND: Genetic screening for pathogenic variants (PVs) in cancer predisposition genes can affect treatment strategies, risk prediction and preventive measures for patients and families. For decades, hereditary breast and ovarian cancer (HBOC) has been attributed to PVs in the genes BRCA1 and BRCA2, and more recently other rare alleles have been firmly established as associated with a high or moderate increased risk of developing breast and/or ovarian cancer. Here, we assess the genetic variation and tumor characteristics in a large cohort of women with suspected HBOC in a clinical oncogenetic setting. METHODS: Women with suspected HBOC referred from all oncogenetic clinics in Sweden over a six-year inclusion period were screened for PVs in 13 clinically relevant genes. The genetic outcome was compared with tumor characteristics and other clinical data collected from national cancer registries and hospital records. RESULTS: In 4622 women with breast and/or ovarian cancer the overall diagnostic yield (the proportion of women carrying at least one PV) was 16.6%. BRCA1/2 PVs were found in 8.9% of women (BRCA1 5.95% and BRCA2 2.94%) and PVs in the other breast and ovarian cancer predisposition genes in 8.2%: ATM (1.58%), BARD1 (0.45%), BRIP1 (0.43%), CDH1 (0.11%), CHEK2 (3.46%), PALB2 (0.84%), PTEN (0.02%), RAD51C (0.54%), RAD51D (0.15%), STK11 (0) and TP53 (0.56%). Thus, inclusion of the 11 genes in addition to BRCA1/2 increased diagnostic yield by 7.7%. The yield was, as expected, significantly higher in certain subgroups such as younger patients, medullary breast cancer, higher Nottingham Histologic Grade, ER-negative breast cancer, triple-negative breast cancer and high grade serous ovarian cancer. Age and tumor subtype distributions differed substantially depending on genetic finding. CONCLUSIONS: This study contributes to understanding the clinical and genetic landscape of breast and ovarian cancer susceptibility. Extending clinical genetic screening from BRCA1 and BRCA2 to 13 established cancer predisposition genes almost doubles the diagnostic yield, which has implications for genetic counseling and clinical guidelines. The very low yield in the syndrome genes CDH1, PTEN and STK11 questions the usefulness of including these genes on routine gene panels.

Somatic mutations reveal asymmetric cellular dynamics in the early human embryo.

Somatic cells acquire mutations throughout the course of an individual's life. Mutations occurring early in embryogenesis are often present in a substantial proportion of, but not all, cells in postnatal humans and thus have particular characteristics and effects. Depending on their location in the genome and the proportion of cells they are present in, these mosaic mutations can cause a wide range of genetic disease syndromes and predispose carriers to cancer. They have a high chance of being transmitted to offspring as de novo germline mutations and, in principle, can provide insights into early human embryonic cell lineages and their contributions to adult tissues. Although it is known that gross chromosomal abnormalities are remarkably common in early human embryos, our understanding of early embryonic somatic mutations is very limited. Here we use whole-genome sequences of normal blood from 241 adults to identify 163 early embryonic mutations. We estimate that approximately three base substitution mutations occur per cell per cell-doubling event in early human embryogenesis and these are mainly attributable to two known mutational signatures. We used the mutations to reconstruct developmental lineages of adult cells and demonstrate that the two daughter cells of many early embryonic cell-doubling events contribute asymmetrically to adult blood at an approximately 2:1 ratio. This study therefore provides insights into the mutation rates, mutational processes and developmental outcomes of cell dynamics that operate during early human embryogenesis.

Merged testing for colorectal cancer syndromes and re-evaluation of genetic variants improve diagnostic yield: Results from a nationwide prospective cohort.

Approximately 5% of patients with colorectal cancer (CRC) have a Mendelian predisposition for the disease. Identification of the disease-causing genetic variant enables carrier testing and tailored cancer prevention within affected families. To determine the panorama and genetic variation of Mendelian CRC syndromes among referrals at the cancer genetics clinics in Sweden, 850 patients clinically selected for CRC genetic investigation were included in a prospective study that tested for all major hereditary polyposis and nonpolyposis CRC conditions. Genetically defined syndromes were diagnosed in 11% of the patients. Lynch syndrome was predominant (n = 73) followed by familial adenomatous polyposis (n = 12) and MUTYH-associated polyposis (n = 8); the latter of which two patients presented with CRC before polyposis was evident. One patient with a history of adolescent-onset CRC and polyposis had biallelic disease-causing variants diagnostic for constitutional mismatch repair deficiency syndrome. Post-study review of detected variants of unknown clinical significance (n = 129) resulted in the reclassification of variants as likely benign (n = 59) or as diagnostic for Lynch syndrome (n = 2). Our results reveal the panorama of Mendelian CRC syndromes at the cancer genetics clinics in Sweden and show that unified testing for polyposis and nonpolyposis CRC conditions as well as regular reexamination of sequence data improve the diagnostic yield.

Clinical, splicing, and functional analysis to classify BRCA2 exon 3 variants: Application of a points-based ACMG/AMP approach.

Skipping of BRCA2 exon 3 (∆E3) is a naturally occurring splicing event, complicating clinical classification of variants that may alter ∆E3 expression. This study used multiple evidence types to assess pathogenicity of 85 variants in/near BRCA2 exon 3. Bioinformatically predicted spliceogenic variants underwent mRNA splicing analysis using minigenes and/or patient samples. ∆E3 was measured using quantitative analysis. A mouse embryonic stem cell (mESC) based assay was used to determine the impact of 18 variants on mRNA splicing and protein function. For each variant, population frequency, bioinformatic predictions, clinical data, and existing mRNA splicing and functional results were collated. Variant class was assigned using a gene-specific adaptation of ACMG/AMP guidelines, following a recently proposed points-based system. mRNA and mESC analysis combined identified six variants with transcript and/or functional profiles interpreted as loss of function. Cryptic splice site use for acceptor site variants generated a transcript encoding a shorter protein that retains activity. Overall, 69/85 (81%) variants were classified using the points-based approach. Our analysis shows the value of applying gene-specific ACMG/AMP guidelines using a points-based approach and highlights the consideration of cryptic splice site usage to appropriately assign PVS1 code strength.

High-definition spatial transcriptomics for in situ tissue profiling.

Spatial and molecular characteristics determine tissue function, yet high-resolution methods to capture both concurrently are lacking. Here, we developed high-definition spatial transcriptomics, which captures RNA from histological tissue sections on a dense, spatially barcoded bead array. Each experiment recovers several hundred thousand transcript-coupled spatial barcodes at 2-µm resolution, as demonstrated in mouse brain and primary breast cancer. This opens the way to high-resolution spatial analysis of cells and tissues.

Risks of breast and ovarian cancer for women harboring pathogenic missense variants in BRCA1 and BRCA2 compared with those harboring protein truncating variants.

PURPOSE: Germline genetic testing for BRCA1 and BRCA2 variants has been a part of clinical practice for >2 decades. However, no studies have compared the cancer risks associated with missense pathogenic variants (PVs) with those associated with protein truncating (PTC) variants. METHODS: We collected 582 informative pedigrees segregating 1 of 28 missense PVs in BRCA1 and 153 pedigrees segregating 1 of 12 missense PVs in BRCA2. We analyzed 324 pedigrees with PTC variants in BRCA1 and 214 pedigrees with PTC variants in BRCA2. Cancer risks were estimated using modified segregation analysis. RESULTS: Estimated breast cancer risks were markedly lower for women aged >50 years carrying BRCA1 missense PVs than for the women carrying BRCA1 PTC variants (hazard ratio [HR] = 3.9 [2.4-6.2] for PVs vs 12.8 [5.7-28.7] for PTC variants; P = .01), particularly for missense PVs in the BRCA1 C-terminal domain (HR = 2.8 [1.4-5.6]; P = .005). In case of BRCA2, for women aged >50 years, the HR was 3.9 (2.0-7.2) for those heterozygous for missense PVs compared with 7.0 (3.3-14.7) for those harboring PTC variants. BRCA1 p.[Cys64Arg] and BRCA2 p.[Trp2626Cys] were associated with particularly low risks of breast cancer compared with other PVs. CONCLUSION: These results have important implications for the counseling of at-risk women who harbor missense PVs in the BRCA1/2 genes.

Transcriptome-wide association study of breast cancer risk by estrogen-receptor status.

Previous transcriptome-wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome-wide association studies (GWAS), but analyses of breast cancer subtype-specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta-analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER-). We further compared associations with ER+ and ER- subtypes, using a case-only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER- breast cancer. We further identified 30 TWAS-significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast-cancer gene in three of six regions harboring multiple TWAS-significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER- breast cancer.

Molecular analyses of triple-negative breast cancer in the young and elderly.

BACKGROUND: Breast cancer in young adults has been implicated with a worse outcome. Analyses of genomic traits associated with age have been heterogenous, likely because of an incomplete accounting for underlying molecular subtypes. We aimed to resolve whether triple-negative breast cancer (TNBC) in younger versus older patients represent similar or different molecular diseases in the context of genetic and transcriptional subtypes and immune cell infiltration. PATIENTS AND METHODS: In total, 237 patients from a reported population-based south Swedish TNBC cohort profiled by RNA sequencing and whole-genome sequencing (WGS) were included. Patients were binned in 10-year intervals. Complimentary PD-L1 and CD20 immunohistochemistry and estimation of tumor-infiltrating lymphocytes (TILs) were performed. Cases were analyzed for differences in patient outcome, genomic, transcriptional, and immune landscape features versus age at diagnosis. Additionally, 560 public WGS breast cancer profiles were used for validation. RESULTS: Median age at diagnosis was 62 years (range 26-91). Age was not associated with invasive disease-free survival or overall survival after adjuvant chemotherapy. Among the BRCA1-deficient cases (82/237), 90% were diagnosed before the age of 70 and were predominantly of the basal-like subtype. In the full TNBC cohort, reported associations of patient age with changes in Ki67 expression, PIK3CA mutations, and a luminal androgen receptor subtype were confirmed. Within DNA repair deficiency or gene expression defined molecular subgroups, age-related alterations in, e.g., overall gene expression, immune cell marker gene expression, genetic mutational and rearrangement signatures, amount of copy number alterations, and tumor mutational burden did, however, not appear distinct. Similar non-significant associations for genetic alterations with age were obtained for other breast cancer subgroups in public WGS data. Consistent with age-related immunosenescence, TIL counts decreased linearly with patient age across different genetic TNBC subtypes. CONCLUSIONS: Age-related alterations in TNBC, as well as breast cancer in general, need to be viewed in the context of underlying genomic phenotypes. Based on this notion, age at diagnosis alone does not appear to provide an additional layer of biological complexity above that of proposed genetic and transcriptional phenotypes of TNBC. Consequently, treatment decisions should be less influenced by age and more driven by tumor biology.

Distinct mechanisms of resistance to fulvestrant treatment dictate level of ER independence and selective response to CDK inhibitors in metastatic breast cancer.

BACKGROUND: Resistance to endocrine treatment in metastatic breast cancer is a major clinical challenge. Clinical tools to predict both drug resistance and possible treatment combination approaches to overcome it are lacking. This unmet need is mainly due to the heterogeneity underlying both the mechanisms involved in resistance development and breast cancer itself. METHODS: To study the complexity of the mechanisms involved in the resistance to the selective estrogen receptor degrader (SERD) fulvestrant, we performed comprehensive biomarker analyses using several in vitro models that recapitulate the heterogeneity of developed resistance. We further corroborated our findings in tissue samples from patients treated with fulvestrant. RESULTS: We found that different in vitro models of fulvestrant resistance show variable stability in their phenotypes, which corresponded with distinct genomic alterations. Notably, the studied models presented adaptation at different cell cycle nodes to facilitate progression through the cell cycle and responded differently to CDK inhibitors. Cyclin E2 overexpression was identified as a biomarker of a persistent fulvestrant-resistant phenotype. Comparison of pre- and post-treatment paired tumor biopsies from patients treated with fulvestrant revealed an upregulation of cyclin E2 upon development of resistance. Moreover, overexpression of this cyclin was found to be a prognostic factor determining resistance to fulvestrant and shorter progression-free survival. CONCLUSIONS: These data highlight the complexity of estrogen receptor positive breast cancer and suggest that the development of diverse resistance mechanisms dictate levels of ER independence and potentially cross-resistance to CDK inhibitors.

Association between breast cancer risk and disease aggressiveness: Characterizing underlying gene expression patterns.

The association between breast cancer risk defined by the Tyrer-Cuzick score (TC) and disease prognosis is not well established. Here, we investigated the relationship between 5-year TC and disease aggressiveness and then characterized underlying molecular processes. In a case-only study (n = 2474), we studied the association of TC with molecular subtypes and tumor characteristics. In a subset of patients (n = 672), we correlated gene expression to TC and computed a low-risk TC gene expression (TC-Gx) profile, that is, a profile expected to be negatively associated with risk, which we used to test for association with disease aggressiveness. We performed enrichment analysis to pinpoint molecular processes likely to be altered in low-risk tumors. A higher TC was found to be inversely associated with more aggressive surrogate molecular subtypes and tumor characteristics (P < .05) including Ki-67 proliferation status (P < 5 × 10-07 ). Our low-risk TC-Gx, based on the weighted sum of 37 expression values of genes strongly correlated with TC, was associated with basal-like (P < 5 × 10-13 ), HER2-enriched subtype (P < 5 × 10-07 ) and worse 10-year breast cancer-specific survival (log-rank P < 5 × 10-04 ). Associations between low-risk TC-Gx and more aggressive molecular subtypes were replicated in an independent cohort from The Cancer Genome Atlas database (n = 975). Gene expression that correlated with low TC was enriched in proliferation and oncogenic signaling pathways (FDR < 0.05). Moreover, higher proliferation was a key factor explaining the association with worse survival. Women who developed breast cancer despite having a low risk were diagnosed with more aggressive tumors and had a worse prognosis, most likely driven by increased proliferation. Our findings imply the need to establish risk factors associated with more aggressive breast cancer subtypes.

The predictive ability of the 313 variant-based polygenic risk score for contralateral breast cancer risk prediction in women of European ancestry with a heterozygous BRCA1 or BRCA2 pathogenic variant.

PURPOSE: To evaluate the association between a previously published 313 variant-based breast cancer (BC) polygenic risk score (PRS313) and contralateral breast cancer (CBC) risk, in BRCA1 and BRCA2 pathogenic variant heterozygotes. METHODS: We included women of European ancestry with a prevalent first primary invasive BC (BRCA1 = 6,591 with 1,402 prevalent CBC cases; BRCA2 = 4,208 with 647 prevalent CBC cases) from the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA), a large international retrospective series. Cox regression analysis was performed to assess the association between overall and ER-specific PRS313 and CBC risk. RESULTS: For BRCA1 heterozygotes the estrogen receptor (ER)-negative PRS313 showed the largest association with CBC risk, hazard ratio (HR) per SD = 1.12, 95% confidence interval (CI) (1.06-1.18), C-index = 0.53; for BRCA2 heterozygotes, this was the ER-positive PRS313, HR = 1.15, 95% CI (1.07-1.25), C-index = 0.57. Adjusting for family history, age at diagnosis, treatment, or pathological characteristics for the first BC did not change association effect sizes. For women developing first BC < age 40 years, the cumulative PRS313 5th and 95th percentile 10-year CBC risks were 22% and 32% for BRCA1 and 13% and 23% for BRCA2 heterozygotes, respectively. CONCLUSION: The PRS313 can be used to refine individual CBC risks for BRCA1/2 heterozygotes of European ancestry, however the PRS313 needs to be considered in the context of a multifactorial risk model to evaluate whether it might influence clinical decision-making.

CDKN2A genetic testing in melanoma-prone families in Sweden in the years 2015-2020: implications for novel national recommendations.

Background: Inherited pathogenic variants (PVs) in the CDKN2A gene are among the strongest known risk factors for cutaneous melanoma. Carriers are at high risks to develop multiple primary melanomas and other cancers, in particular pancreatic cancer. In this study, the CDKN2A testing, carried out in Sweden in the years 2015-2020, was evaluated.Materials and methods: Included families had (1) three or more cases of melanoma and/or pancreatic cancer, (2) two melanomas in first-degree relatives, the youngest case <55 years or (3) individuals with three or more multiple primary melanomas, the first before the age of 55 years, and no other affected family members. The included families had at least one affected member that had been tested for CDKN2A PVs.Results: In total, 403 families were included, whereof 913 family members had been diagnosed with cutaneous melanoma and 129 with pancreatic cancer, 33 (8.2%) were found to have PVs in CDKN2A. Frequencies ranged from 0.9% in families with only two melanomas to 43.2% in families with three or more melanoma cases and pancreatic cancer (p < 0.001). The frequency of PVs ranged from 2.1% to 16.5% in families where the youngest case was ≥55 years or <35 years (p = 0.040). In families with or without CDKN2A PVs, 37.6% and 10.0% had melanoma cases that had died from melanoma, respectively (p < 0.001).Discussion: Significant differences were seen in the frequencies of CDKN2A PVs, dependent on numbers or age at diagnosis of melanomas and diagnoses of pancreatic cancers in the family. Further, melanoma cases belonging to families that tested positive for CDKN2A PVs had a significantly higher mortality. To summarize, the current evaluation shows that, with adequately selected criteria to guide genetic testing, CDKN2A PVs are identified at significant frequencies. Identification of carrier families is of importance to ensure that members are enrolled in a preventive surveillance program.

Analysis of fusion transcripts indicates widespread deregulation of snoRNAs and their host genes in breast cancer.

Genomic rearrangements in cancer can join the sequences of two separate genes. Studies of such gene fusion events have mainly focused on identification of fusion proteins from the chimeric transcripts. We have previously investigated how fusions instead can affect the expression of intronic microRNA (miRNA) genes that are encoded within fusion gene partners. Here, we extend our analysis to small nucleolar RNAs (snoRNAs) that also are embedded within protein-coding or noncoding host genes. We found that snoRNA hosts are selectively enriched in fusion transcripts, like miRNA host genes, and that this enrichment is associated with all snoRNA classes. These structural changes may have functional consequences for the cell; proteins involved in the protein translation machinery are overrepresented among snoRNA host genes, a gene architecture assumed to be needed for closely coordinated expression of snoRNAs and host proteins. Our data indicate that this structure is frequently disrupted in cancer. We furthermore observed that snoRNA genes involved in fusions tend to associate with stronger promoters than the natural host, suggesting a mechanism that selects for snoRNA overexpression. In summary, we highlight a previously unexplored frequent structural change in cancer that affects important components of cellular physiology.

Breast cancer survival in Nordic BRCA2 mutation carriers-unconventional association with oestrogen receptor status.

BACKGROUND: The natural history of breast cancer among BRCA2 carriers has not been clearly established. In a previous study from Iceland, positive ER status was a negative prognostic factor. We sought to identify factors that predicted survival after invasive breast cancer in an expanded cohort of BRCA2 carriers. METHODS: We studied 608 women with invasive breast cancer and a pathogenic BRCA2 mutation (variant) from four Nordic countries. Information on prognostic factors and treatment was retrieved from health records and by analysis of archived tissue specimens. Hazard ratios (HR) were estimated for breast cancer-specific survival using Cox regression. RESULTS: About 77% of cancers were ER-positive, with the highest proportion (83%) in patients under 40 years. ER-positive breast cancers were more likely to be node-positive (59%) than ER-negative cancers (34%) (P < 0.001). The survival analysis included 584 patients. Positive ER status was protective in the first 5 years from diagnosis (multivariate HR = 0.49; 95% CI 0.26-0.93, P = 0.03); thereafter, the effect was adverse (HR = 1.91; 95% CI 1.07-3.39, P = 0.03). The adverse effect of positive ER status was limited to women who did not undergo endocrine treatment (HR = 2.36; 95% CI 1.26-4.44, P = 0.01) and patients with intact ovaries (HR = 1.99; 95% CI 1.11-3.59, P = 0.02). CONCLUSIONS: The adverse effect of a positive ER status in BRCA2 carriers with breast cancer may be contingent on exposure to ovarian hormones.

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.

Alternative splicing and ACMG-AMP-2015-based classification of PALB2 genetic variants: an ENIGMA report.

BACKGROUND: PALB2 monoallelic loss-of-function germ-line variants confer a breast cancer risk comparable to the average BRCA2 pathogenic variant. Recommendations for risk reduction strategies in carriers are similar. Elaborating robust criteria to identify loss-of-function variants in PALB2-without incurring overprediction-is thus of paramount clinical relevance. Towards this aim, we have performed a comprehensive characterisation of alternative splicing in PALB2, analysing its relevance for the classification of truncating and splice site variants according to the 2015 American College of Medical Genetics and Genomics-Association for Molecular Pathology guidelines. METHODS: Alternative splicing was characterised in RNAs extracted from blood, breast and fimbriae/ovary-related human specimens (n=112). RNAseq, RT-PCR/CE and CloneSeq experiments were performed by five contributing laboratories. Centralised revision/curation was performed to assure high-quality annotations. Additional splicing analyses were performed in PALB2 c.212-1G>A, c.1684+1G>A, c.2748+2T>G, c.3113+5G>A, c.3350+1G>A, c.3350+4A>C and c.3350+5G>A carriers. The impact of the findings on PVS1 status was evaluated for truncating and splice site variant. RESULTS: We identified 88 naturally occurring alternative splicing events (81 newly described), including 4 in-frame events predicted relevant to evaluate PVS1 status of splice site variants. We did not identify tissue-specific alternate gene transcripts in breast or ovarian-related samples, supporting the clinical relevance of blood-based splicing studies. CONCLUSIONS: PVS1 is not necessarily warranted for splice site variants targeting four PALB2 acceptor sites (exons 2, 5, 7 and 10). As a result, rare variants at these splice sites cannot be assumed pathogenic/likely pathogenic without further evidences. Our study puts a warning in up to five PALB2 genetic variants that are currently reported as pathogenic/likely pathogenic in ClinVar.