Gene expression profiling assigns CHEK2 1100delC breast cancers to the luminal intrinsic subtypes

Clinicopathologic and genetic analysis of invasive breast carcinomas in women with germline CHEK2 variants

PURPOSE

Germline pathogenic variants in checkpoint kinase 2 (CHEK2) are associated with a moderately increased risk of breast cancer (BC). The spectrum of clinicopathologic features and genetics of these tumors has not been fully established.

METHODS

We characterized the histopathologic and clinicopathologic features of 44 CHEK2-associated BCs from 35 women, and assessed responses to neoadjuvant chemotherapy. A subset of cases (n = 23) was additionally analyzed using targeted next-generation DNA sequencing (NGS).

RESULTS

Most (94%, 33/35) patients were heterozygous carriers for germline CHEK2 variants, and 40% had the c.1100delC allele. Two patients were homozygous, and five had additional germline pathogenic variants in ATM (2), PALB2 (1), RAD50 (1), or MUTYH (1). CHEK2-associated BCs occurred in younger women (median age 45 years, range 25-75) and were often multifocal (20%) or bilateral (11%). Most (86%, 38/44) were invasive ductal carcinomas of no special type (IDC-NST). Almost all (95%, 41/43) BCs were ER + (79% ER + HER2-, 16% ER + HER2 + , 5% ER-HER2 +), and most (69%) were luminal B. Nottingham grade, proliferation index, and results of multiparametric molecular testing were heterogeneous. Biallelic CHEK2 alteration with loss of heterozygosity was identified in most BCs (57%, 13/23) by NGS. Additional recurrent alterations included GATA3 (26%), PIK3CA (226%), CCND1 (22%), FGFR1 (22%), ERBB2 (17%), ZNF703 (17%), TP53 (9%), and PPM1D (9%), among others. Responses to neoadjuvant chemotherapy were variable, but few patients (21%, 3/14) achieved pathologic complete response. Most patients (85%) were without evidence of disease at time of study (n = 34). Five patients (15%) developed distant metastasis, and one (3%) died (mean follow-up 50 months).

CONCLUSION

Almost all CHEK2-associated BCs were ER + IDC-NST, with most classified as luminal B with or without HER2 overexpression. NGS supported the luminal-like phenotype and confirmed CHEK2 as an oncogenic driver in the majority of cases. Responses to neoadjuvant chemotherapy were variable but mostly incomplete.

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.

HER2-low breast cancer shows a lower immune response compared to HER2-negative cases

Currently, the human epidermal growth factor receptor 2 (HER2) status of breast cancer is classified dichotomously as negative or positive to select patients for HER2-targeted therapy. However, with the introduction of novel treatment options, it is important to get more insight in the biology of cancers with low HER2 expression. Therefore, we studied several clinicopathologic characteristics in relation to the level of HER2 expression (HER2- versus HER2low). We used a well-documented cohort of breast cancer patients (n = 529), with available tissue microarrays and Affymetrix mRNA expression data. HER2 status was scored as negative (immunohistochemistry 0) or low (immunohistochemistry 1 + or 2 + without amplification). We associated HER2 status with several clinicopathologic characteristics, gene-expression data and survival, stratified for estrogen receptor (ER) status. Overall, breast cancers were scored as HER2- (n = 429) or HER2low (n = 100). Within the ER+ cohort (n = 305), no significant associations were found between the HER2 groups and clinicopathologic features. However, HER2low tumors showed several differentially expressed genes compared to HER2- cases, including genes that are associated with worse outcome and depletion of immunity. In ER- cases (n = 224), HER2low status was significantly associated with increased regional nodal positivity, lower density of tumor infiltrating lymphocyte and a lower protein expression of Ki-67 and EGFR compared to HER2- cases. After multivariate analysis, only density of tumor infiltrating lymphocytes remained significantly associated with HER2low status (P = 0.035). No difference in survival was observed between HER2low and HER2- patients, neither in the ER+ nor ER- cohort. In conclusion, our data suggests that HER2low breast cancer is associated with a lower immune response compared to HER2- breast cancer.

Somatic loss of the remaining allele occurs approximately in half of CHEK2-driven breast cancers and is accompanied by a border-line increase of chromosomal instability

PURPOSE

Germline mutations in CHEK2 gene represent the second most frequent cause of hereditary breast cancer (BC) after BRCA1/2 lesions. This study aimed to identify the molecular characteristics of CHEK2-driven BCs.

METHODS

Loss of heterozygosity (LOH) for the remaining CHEK2 allele was examined in 50 CHEK2-driven BCs using allele-specific PCR assays for the germline mutations and analysis of surrounding single-nucleotide polymorphisms (SNPs). Paired tumor and normal DNA samples from 25 cases were subjected to next-generation sequencing analysis.

RESULTS

CHEK2 LOH was detected in 28/50 (56%) BCs. LOH involved the wild-type allele in 24 BCs, mutant CHEK2 copy was deleted in 3 carcinomas, while in one case the origin of the deleted allele could not be identified. Somatic PIK3CA and TP53 mutations were present in 13/25 (52%) and 4/25 (16%) tumors, respectively. Genomic features of homologous recombination deficiency (HRD), including the HRD score ≥ 42, the predominance of BRCA-related mutational signature 3, and the high proportion of long (≥ 5 bp) indels, were observed only in 1/20 (5%) BC analyzed for chromosomal instability. Tumors with the deleted wild-type CHEK2 allele differed from LOH-negative cases by elevated HRD scores (median 23 vs. 7, p = 0.010) and higher numbers of chromosomal segments affected by copy number aberrations (p = 0.008).

CONCLUSION

Somatic loss of the wild-type CHEK2 allele is observed in approximately half of CHEK2-driven BCs. Tumors without CHEK2 LOH are chromosomally stable. BCs with LOH demonstrate some signs of chromosomal instability; however, its degree is significantly lower as compared to BRCA1/2-associated cancers.

Population-based estimates of breast cancer risk for carriers of pathogenic variants identified by gene-panel testing

Population-based estimates of breast cancer risk for carriers of pathogenic variants identified by gene-panel testing are urgently required. Most prior research has been based on women selected for high-risk features and more data is needed to make inference about breast cancer risk for women unselected for family history, an important consideration of population screening. We tested 1464 women diagnosed with breast cancer and 862 age-matched controls participating in the Australian Breast Cancer Family Study (ABCFS), and 6549 healthy, older Australian women enroled in the ASPirin in Reducing Events in the Elderly (ASPREE) study for rare germline variants using a 24-gene-panel. Odds ratios (ORs) were estimated using unconditional logistic regression adjusted for age and other potential confounders. We identified pathogenic variants in 11.1% of the ABCFS cases, 3.7% of the ABCFS controls and 2.2% of the ASPREE (control) participants. The estimated breast cancer OR [95% confidence interval] was 5.3 [2.1–16.2] for BRCA1, 4.0 [1.9–9.1] for BRCA2, 3.4 [1.4–8.4] for ATM and 4.3 [1.0–17.0] for PALB2. Our findings provide a population-based perspective to gene-panel testing for breast cancer predisposition and opportunities to improve predictors for identifying women who carry pathogenic variants in breast cancer predisposition genes.

Increased gene expression variability in BRCA1-associated and basal-like breast tumours

Purpose

Inherited variants in the cancer susceptibility genes, BRCA1 and BRCA2 account for up to 5% of breast cancers. Multiple gene expression studies have analysed gene expression patterns that maybe associated with BRCA12 pathogenic variant status; however, results from these studies lack consensus. These studies have focused on the differences in population means to identified genes associated with BRCA1/2-carriers with little consideration for gene expression variability, which is also under genetic control and is a feature of cellular function.

Methods

We measured differential gene expression variability in three of the largest familial breast cancer datasets and a 2116 breast cancer meta-cohort. Additionally, we used RNA in situ hybridisation to confirm expression variability of EN1 in an independent cohort of more than 500 breast tumours.

Results

BRCA1-associated breast tumours exhibited a 22.8% (95% CI 22.3–23.2) increase in transcriptome-wide gene expression variability compared to BRCAx tumours. Additionally, 40 genes were associated with BRCA1-related breast cancers that had ChIP-seq data suggestive of enriched EZH2 binding. Of these, two genes (EN1 and IGF2BP3) were significantly variable in both BRCA1-associated and basal-like breast tumours. RNA in situ analysis of EN1 supported a significant (p = 6.3 × 10⁻⁰⁴) increase in expression variability in BRCA1-associated breast tumours.

Conclusion

Our novel results describe a state of increased gene expression variability in BRCA1-related and basal-like breast tumours. Furthermore, genes with increased variability may be driven by changes in DNA occupancy of epigenetic effectors. The variation in gene expression is replicable and led to the identification of novel associations between genes and disease phenotypes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10549-021-06328-y.

Investigation of monogenic causes of familial breast cancer: data from the BEACCON case-control study

Breast cancer (BC) has a significant heritable component but the genetic contribution remains unresolved in the majority of high-risk BC families. This study aims to investigate the monogenic causes underlying the familial aggregation of BC beyond BRCA1 and BRCA2, including the identification of new predisposing genes. A total of 11,511 non-BRCA familial BC cases and population-matched cancer-free female controls in the BEACCON study were investigated in two sequencing phases: 1303 candidate genes in up to 3892 cases and controls, followed by validation of 145 shortlisted genes in an additional 7619 subjects. The coding regions and exon–intron boundaries of all candidate genes and 14 previously proposed BC genes were sequenced using custom designed sequencing panels. Pedigree and pathology data were analysed to identify genotype-specific associations. The contribution of ATM, PALB2 and CHEK2 to BC predisposition was confirmed, but not RAD50 and NBN. An overall excess of loss-of-function (LoF) (OR 1.27, p = 9.05 × 10⁻⁹) and missense (OR 1.27, p = 3.96 × 10⁻⁷³) variants was observed in the cases for the 145 candidate genes. Leading candidates harbored LoF variants with observed ORs of 2–4 and individually accounted for no more than 0.79% of the cases. New genes proposed by this study include NTHL1, WRN, PARP2, CTH and CDK9. The new candidate BC predisposition genes identified in BEACCON indicate that much of the remaining genetic causes of high-risk BC families are due to genes in which pathogenic variants are both very rare and convey only low to moderate risk.

Population-Based Estimates of the Age-Specific Cumulative Risk of Breast Cancer for Pathogenic Variants in CHEK2: Findings from the Australian Breast Cancer Family Registry

Simple Summary

It is well established that women who carry pathogenic CHEK2 variants have about a 3-fold increased risk of developing breast cancer. CHEK2 is now commonly included in genetic tests for breast cancer predisposition and increasingly used to inform the clinical management of women who are identified to carry pathogenic variants. Important information for counselling these women includes knowing how breast cancer risk, due to having a pathogenic variant in CHEK2, changes over a woman’s lifetime. This information is currently not well established. By conducting a population-based case-control-family study of pathogenic CHEK2 variants we aimed to provide this information and estimated the penetrance (age-specific cumulative risk) of breast cancer to be 18% (95% CI 11–30%) to age 60 years and 33% (95% CI 21–48%) to age 80 years. These findings provide new and important information for the clinical management of breast cancer risk for women carrying pathogenic variants in CHEK2.

Abstract

Case-control studies of breast cancer have consistently shown that pathogenic variants in CHEK2 are associated with about a 3-fold increased risk of breast cancer. Information about the recurrent protein-truncating variant CHEK2 c.1100delC dominates this estimate. There have been no formal estimates of age-specific cumulative risk of breast cancer for all CHEK2 pathogenic (including likely pathogenic) variants combined. We conducted a population-based case-control-family study of pathogenic CHEK2 variants (26 families, 1071 relatives) and estimated the age-specific cumulative risk of breast cancer using segregation analysis. The estimated hazard ratio for carriers of pathogenic CHEK2 variants (combined) was 4.9 (95% CI 2.5–9.5) relative to non-carriers. The HR for carriers of the CHEK2 c.1100delC variant was estimated to be 3.5 (95% CI 1.02–11.6) and the HR for carriers of all other CHEK2 variants combined was estimated to be 5.7 (95% CI 2.5–12.9). The age-specific cumulative risk of breast cancer was estimated to be 18% (95% CI 11–30%) and 33% (95% CI 21–48%) to age 60 and 80 years, respectively. These findings provide important information for the clinical management of breast cancer risk for women carrying pathogenic variants in CHEK2.

CHEK2 Germline Variants in Cancer Predisposition: Stalemate Rather than Checkmate

Germline alterations in many genes coding for proteins regulating DNA repair and DNA damage response (DDR) to DNA double-strand breaks (DDSB) have been recognized as pathogenic factors in hereditary cancer predisposition. The ATM-CHEK2-p53 axis has been documented as a backbone for DDR and hypothesized as a barrier against cancer initiation. However, although CHK2 kinase coded by the CHEK2 gene expedites the DDR signal, its function in activation of p53-dependent cell cycle arrest is dispensable. CHEK2 mutations rank among the most frequent germline alterations revealed by germline genetic testing for various hereditary cancer predispositions, but their interpretation is not trivial. From the perspective of interpretation of germline CHEK2 variants, we review the current knowledge related to the structure of the CHEK2 gene, the function of CHK2 kinase, and the clinical significance of CHEK2 germline mutations in patients with hereditary breast, prostate, kidney, thyroid, and colon cancers.

Genetics and Genomics of Breast Cancer: update and translational perspectives

In the recent years the rapid scientific innovation in the evaluation of the individual's genome have allowed the identification of variants associated with the onset, treatment and prognosis of various pathologies including cancer, and with a potential impact in the assessment of therapy responses. Despite the analysis and interpretation of genomic information is considered incomplete, in many cases the identification of specific genomic profile has allowed the stratification of subgroups of patients characterized by a better response to drug therapies. Individual genome analysis has changed profoundly the diagnostic and therapeutic approach of breast cancer in the last 15 years by identifying selective molecular lesions that drive the development of neoplasms, showing that each tumor has its own genomic signature, with some specific features and some features common to several sub-types. Several personalized therapies have been (and still are being) developed showing a remarkable efficacy in the treatment of breast cancer.

TEffectR: an R package for studying the potential effects of transposable elements on gene expression with linear regression model

Introduction

Recent studies highlight the crucial regulatory roles of transposable elements (TEs) on proximal gene expression in distinct biological contexts such as disease and development. However, computational tools extracting potential TE -proximal gene expression associations from RNA-sequencing data are still missing.

Implementation

Herein, we developed a novel R package, using a linear regression model, for studying the potential influence of TE species on proximal gene expression from a given RNA-sequencing data set. Our R package, namely TEffectR, makes use of publicly available RepeatMasker TE and Ensembl gene annotations as well as several functions of other R-packages. It calculates total read counts of TEs from sorted and indexed genome aligned BAM files provided by the user, and determines statistically significant relations between TE expression and the transcription of nearby genes under diverse biological conditions.

Availability

TEffectR is freely available at https://github.com/karakulahg/TEffectR along with a handy tutorial as exemplified by the analysis of RNA-sequencing data including normal and tumour tissue specimens obtained from breast cancer patients.

Clonality, Antigen Recognition, and Suppression of CD8+ T Cells Differentially Affect Prognosis of Breast Cancer Subtypes

PURPOSE

In breast cancer, response rates to immune therapies are generally low and differ significantly across molecular subtypes, urging a better understanding of immunogenicity and immune evasion.

EXPERIMENTAL DESIGN

We interrogated large gene-expression data sets including 867 node-negative, treatment-naïve breast cancer patients (microarray data) and 347 breast cancer patients (whole-genome sequencing and transcriptome data) according to parameters of T cells as well as immune microenvironment in relation to patient survival.

RESULTS

We developed a 109-immune gene signature that captures abundance of CD8 tumor-infiltrating lymphocytes (TIL) and is prognostic in basal-like, her2, and luminal B breast cancer, but not in luminal A or normal-like breast cancer. Basal-like and her2 are characterized by highest CD8 TIL abundance, highest T-cell clonality, highest frequencies of memory T cells, and highest antigenicity, yet only the former shows highest expression level of immune and metabolic checkpoints and highest frequency of myeloid suppressor cells. Also, luminal B shows a high antigenicity and T-cell clonality, yet a low abundance of CD8 TILs. In contrast, luminal A and normal-like both show a low antigenicity, and notably, a low and high abundance of CD8 TILs, respectively, which associates with T-cell influx parameters, such as expression of adhesion molecules.

CONCLUSIONS

Collectively, our data argue that not only CD8 T-cell presence itself, but rather T-cell clonality, T-cell subset distribution, coinhibition, and antigen presentation reflect occurrence of a CD8 T-cell response in breast cancer subtypes, which have been aborted by distinct T-cell-suppressive mechanisms, providing a rationale for subtype-specific combination immune therapies.

Germline Genetic Testing for Breast Cancer Risk: The Past, Present, and Future

The field of germline genetic testing for breast cancer (BC) risk has evolved substantially in the last decade. The introduction of multigene panel testing (MGPT) led to an urgent need to understand the cancer risk associated with several genes included in the panels. Although the research on understanding the cancer risk associated with mutations in several genes continues, there is also a need to understand the modifying effects of race and ethnicity, family history, and BC pathology on the prevalence of germline mutations and associated BC risk. Furthermore, polygenic risk scores (PRSs) to predict BC risk in patients with or without germline mutations in cancer-predisposition genes are now available for clinical use, although data on the clinical utility of PRSs are lacking. In patients with advanced BC associated with BRCA1/2 mutation, olaparib and talazoparib are now approved for treatment. In addition, molecular profiling studies are being used to clarify the BC tumor biology in mutation carriers to identify potential therapeutic options. In this article, we discuss these advances in the field of germline genetic testing and highlight current limitations and implications for clinical care.

Transcriptome mining of non-BRCA1/A2 and BRCA1/A2 familial breast cancer

About 10% of all breast cancer cases are the familial type. Mutations in two highly penetrance breast cancer susceptibility genes, BRCA1 and BRCA2, can only explain 20% to 25% of genetic susceptibility to breast cancer, and most familial breast cancer cases have intact BRCA1 and BRCA2 genes that refer to non-BRCA1/A2 or BRCAX familial breast cancer. Despite extensive studies, more than 50% of genetic susceptibility to breast cancer remained to be disclosed. Finding the differences between these two types of breast cancer (non-BRCA1/A2 and BRCA1/A2) at genomic, transcriptomic, and proteomic levels can help us to elucidate fundamental molecular processes and develope more promising therapeutic targets. Here, we used expression data of 391 patients with familial breast cancer including 195 non-BRCA1/A2 and 196 BRCA1 and/or BRCA2 cases from four independent studies by means of meta-analysis to find differences in gene expression signature between these two types of familial breast cancer. As well as, we applied comprehensive network analysis to find crucial protein complexes and regulators for each condition. Our results revealed significant overexpression of cell cycle processes in BRCA1/A2 patients and significant overexpression of estrogen axis in non-BRCA1/A2 patients. Moreover, we found FOXM1 as the central regulator of cell cycle processes and GATA3, FOXA1, and ESR1 as the main regulators of estrogen axis.

Identification and analysis of CHEK2 germline mutations in Chinese BRCA1/2-negative breast cancer patients

PURPOSE

Cell-cycle-checkpoint kinase 2 (CHEK2) is an important moderate-penetrance breast cancer predisposition gene; however, recurrent CHEK2 mutations found in Caucasian women are very rare in Chinese population. We investigated the mutation spectrum and clinical relevance of CHEK2 germline mutations in Chinese breast cancer patients.

METHODS

The entire coding regions and splicing sites of CHEK2 were screened in 7657 Chinese BRCA1/2-negative breast cancer patients, using 62-gene panel-based sequencing.

RESULTS

Out of 7657 BRCA1/2-negative breast cancer patients, 26 (0.34%) carried CHEK2 pathogenic germline mutations. Most of these mutations (92.3%, 24/26) were nonsense or frameshift mutations; 84.6% (22/26) of them were in forkhead-associated (FHA) or kinase domains. Of the 18 types of CHEK2 mutations we found, 61.1% (11/18) of were novel mutations and two recurrent mutations (Y139X and R137X) were found in this cohort. Patients with CHEK2 mutations were significantly more likely to have family histories of breast and/or ovarian cancer (23.1% vs. 8.6%, p = 0.022) and family histories of any cancer (50.0% vs. 31.6%, p = 0.044); and were significantly more likely to have lymph node-positive (53.8% vs. 27.3%, p = 0.002) and progesterone receptor (PR)-positive (88.5% vs. 64.5%, p = 0.011) breast cancers.

CONCLUSIONS

Among Chinese breast cancer patients, the CHEK2 germline mutation rate is approximately 0.34% and two specific mutations (Y139X and R137X) are recurrent. Patients with CHEK2 mutations are significantly more likely to have family histories of cancer, and to develop lymph node-positive and/or PR-positive breast cancers.

Computational Investigation of Homologous Recombination DNA Repair Deficiency in Sporadic Breast Cancer

BRCAness has important implications in the management and treatment of patients with breast and ovarian cancer. In this study, we propose a computational framework to measure the BRCAness of breast and ovarian tumor samples based on their gene expression profiles. We define a characteristic profile for BRCAness by comparing gene expression differences between BRCA1/2 mutant familial tumors and sporadic breast cancer tumors while adjusting for relevant clinical factors. With this BRCAness profile, our framework calculates sample-specific BRCA scores, which indicates homologous recombination (HR)-mediated DNA repair pathway activity of samples. We found that in sporadic breast cancer high BRCAness score is associated with aberrant copy number of HR genes rather than somatic mutation and other genomic features. Moreover, we observed significant correlations of BRCA score with genome instability and neoadjuvant chemotherapy. More importantly, BRCA score provides significant prognostic value in both breast and ovarian cancers after considering established clinical variables. In summary, the inferred BRCAness from our framework can be used as a robust biomarker for the prediction of prognosis and treatment response in breast and ovarian cancers.

Genomic profiling of pelvic genital type leiomyosarcoma in a woman with a germline CHEK2:c.1100delC mutation and a concomitant diagnosis of metastatic invasive ductal breast carcinoma

We describe a woman with the known pathogenic germline variant CHEK2:c.1100delC and synchronous diagnoses of both pelvic genital type leiomyosarcoma (LMS) and metastatic invasive ductal breast carcinoma. CHEK2 (checkpoint kinase 2) is a tumor-suppressor gene encoding a serine/threonine-protein kinase (CHEK2) involved in double-strand DNA break repair and cell cycle arrest. The CHEK2:c.1100delC variant is a moderate penetrance allele resulting in an approximately twofold increase in breast cancer risk. Whole-genome and whole-transcriptome sequencing were performed on the leiomyosarcoma and matched blood-derived DNA. Despite the presence of several genomic hits within the double-strand DNA damage pathway (CHEK2 germline variant and multiple RAD51B somatic structural variants), tumor profiling did not show an obvious DNA repair deficiency signature. However, even though the LMS displayed clear malignant features, its genomic profiling revealed several characteristics classically associated with leiomyomas including a translocation, t(12;14), with one breakpoint disrupting RAD51B and the other breakpoint upstream of HMGA2 with very high expression of HMGA2 and PLAG1. This is the first report of LMS genomic profiling in a patient with the germline CHEK2:c.1100delC variant and an additional diagnosis of metastatic invasive ductal breast carcinoma. We also describe a possible mechanistic relationship between leiomyoma and LMS based on genomic and transcriptome data. Our findings suggest that RAD51B translocation and HMGA2 overexpression may play an important role in LMS oncogenesis.

Transcriptional signature of lymphoblastoid cell lines of BRCA1, BRCA2 and non-BRCA1/2 high risk breast cancer families

Approximately 25% of hereditary breast cancer cases are associated with a strong familial history which can be explained by mutations in BRCA1 or BRCA2 and other lower penetrance genes. The remaining high-risk families could be classified as BRCAX (non-BRCA1/2) families.

Gene expression involving alternative splicing represents a well-known mechanism regulating the expression of multiple transcripts, which could be involved in cancer development. Thus using RNA-seq methodology, the analysis of transcriptome was undertaken to potentially reveal transcripts implicated in breast cancer susceptibility and development.

RNA was extracted from immortalized lymphoblastoid cell lines of 117 women (affected and unaffected) coming from BRCA1, BRCA2 and BRCAX families. Anova analysis revealed a total of 95 transcripts corresponding to 85 different genes differentially expressed (Bonferroni corrected p-value <0.01) between those groups. Hierarchical clustering allowed distinctive subgrouping of BRCA1/2 subgroups from BRCAX individuals. We found 67 transcripts, which could discriminate BRCAX from BRCA1/BRCA2 individuals while 28 transcripts discriminate affected from unaffected BRCAX individuals.

To our knowledge, this represents the first study identifying transcripts differentially expressed in lymphoblastoid cell lines from major classes of mutation-related breast cancer subgroups, namely BRCA1, BRCA2 and BRCAX. Moreover, some transcripts could discriminate affected from unaffected BRCAX individuals, which could represent potential therapeutic targets for breast cancer treatment.

Patient survival and tumor characteristics associated with CHEK2:p.I157T - findings from the Breast Cancer Association Consortium

BACKGROUND

P.I157T is a CHEK2 missense mutation associated with a modest increase in breast cancer risk. Previously, another CHEK2 mutation, the protein truncating c.1100delC has been associated with poor prognosis of breast cancer patients. Here, we have investigated patient survival and characteristics of breast tumors of germ line p.I157T carriers.

METHODS

We included in the analyses 26,801 European female breast cancer patients from 15 studies participating in the Breast Cancer Association Consortium. We analyzed the association between p.I157T and the clinico-pathological breast cancer characteristics by comparing the p.I157T carrier tumors to non-carrier and c.1100delC carrier tumors. Similarly, we investigated the p.I157T associated risk of early death, breast cancer-associated death, distant metastasis, locoregional relapse and second breast cancer using Cox proportional hazards models. Additionally, we explored the p.I157T-associated genomic gene expression profile using data from breast tumors of 183 Finnish female breast cancer patients (ten p.I157T carriers) (GEO: GSE24450). Differential gene expression analysis was performed using a moderated t test. Functional enrichment was investigated using the DAVID functional annotation tool and gene set enrichment analysis (GSEA). The tumors were classified into molecular subtypes according to the St Gallen 2013 criteria and the PAM50 gene expression signature.

RESULTS

P.I157T was not associated with increased risk of early death, breast cancer-associated death or distant metastasis relapse, and there was a significant difference in prognosis associated with the two CHEK2 mutations, p.I157T and c.1100delC. Furthermore, p.I157T was associated with lobular histological type and clinico-pathological markers of good prognosis, such as ER and PR expression, low TP53 expression and low grade. Gene expression analysis suggested luminal A to be the most common subtype for p.I157T carriers and CDH1 (cadherin 1) target genes to be significantly enriched among genes, whose expression differed between p.I157T and non-carrier tumors.

CONCLUSIONS

Our analyses suggest that there are fundamental differences in breast tumors of CHEK2:p.I157T and c.1100delC carriers. The poor prognosis associated with c.1100delC cannot be generalized to other CHEK2 mutations.

Age- and Tumor Subtype-Specific Breast Cancer Risk Estimates for CHEK2*1100delC Carriers

PURPOSE

CHEK2*1100delC is a well-established breast cancer risk variant that is most prevalent in European populations; however, there are limited data on risk of breast cancer by age and tumor subtype, which limits its usefulness in breast cancer risk prediction. We aimed to generate tumor subtype- and age-specific risk estimates by using data from the Breast Cancer Association Consortium, including 44,777 patients with breast cancer and 42,997 controls from 33 studies genotyped for CHEK2*1100delC.

PATIENTS AND METHODS

CHEK2*1100delC genotyping was mostly done by a custom Taqman assay. Breast cancer odds ratios (ORs) for CHEK2*1100delC carriers versus noncarriers were estimated by using logistic regression and adjusted for study (categorical) and age. Main analyses included patients with invasive breast cancer from population- and hospital-based studies.

RESULTS

Proportions of heterozygous CHEK2*1100delC carriers in controls, in patients with breast cancer from population- and hospital-based studies, and in patients with breast cancer from familial- and clinical genetics center-based studies were 0.5%, 1.3%, and 3.0%, respectively. The estimated OR for invasive breast cancer was 2.26 (95%CI, 1.90 to 2.69; P = 2.3 × 10(-20)). The OR was higher for estrogen receptor (ER)-positive disease (2.55 [95%CI, 2.10 to 3.10; P = 4.9 × 10(-21)]) than it was for ER-negative disease (1.32 [95%CI, 0.93 to 1.88; P = .12]; P interaction = 9.9 × 10(-4)). The OR significantly declined with attained age for breast cancer overall (P = .001) and for ER-positive tumors (P = .001). Estimated cumulative risks for development of ER-positive and ER-negative tumors by age 80 in CHEK2*1100delC carriers were 20% and 3%, respectively, compared with 9% and 2%, respectively, in the general population of the United Kingdom.

CONCLUSION

These CHEK2*1100delC breast cancer risk estimates provide a basis for incorporating CHEK2*1100delC into breast cancer risk prediction models and into guidelines for intensified screening and follow-up.

Women at high risk of breast cancer: Molecular characteristics, clinical presentation and management

The presence of breast cancer in any first-degree female relative in general nearly doubles the risk for a proband and the risk gradually increases with the number of affected relatives. Current advances in molecular oncology and oncogenetics may enable the identification of high-risk individuals with breast-cancer predisposition. The best-known forms of hereditary breast cancer (HBC) are caused by mutations in the high-penetrance genes BRCA1 and BRCA2. Other genes, including PTEN, TP53, STK11/LKB1, CDH1, PALB2, CHEK2, ATM, MRE11, RAD50, NBS1, BRIP1, FANCA, FANCC, FANCM, RAD51, RAD51B, RAD51C, RAD51D, and XRCC2 have been described as high- or moderate-penetrance breast cancer-susceptibility genes. The majority of breast cancer-susceptibility genes code for tumor suppressor proteins that are involved in critical processes of DNA repair pathways. This is of particular importance for those women who, due to their increased risk of breast cancer, may be subjected to more frequent screening but due to their repair deficiency might be at the risk of developing radiation-induced malignancies. It has been proven that cancers arising from the most frequent BRCA1 gene mutation carriers differ significantly from the sporadic disease of age-matched controls in their histopathological appearances and molecular characteristics. The increased depth of mutation detection brought by next-generation sequencing and a better understanding of the mechanisms through which these mutations cause the disease will bring novel insights in terms of oncological prevention, diagnostics, and therapeutic options for HBC patients.

Genomic profiling of CHEK2*1100delC-mutated breast carcinomas

BACKGROUND

CHEK2*1100delC is a moderate-risk breast cancer susceptibility allele with a high prevalence in the Netherlands. We performed copy number and gene expression profiling to investigate whether CHEK2*1100delC breast cancers harbor characteristic genomic aberrations, as seen for BRCA1 mutated breast cancers.

METHODS

We performed high-resolution SNP array and gene expression profiling of 120 familial breast carcinomas selected from a larger cohort of 155 familial breast tumors, including BRCA1, BRCA2, and CHEK2 mutant tumors. Gene expression analyses based on a mRNA immune signature was used to identify samples with relative low amounts of tumor infiltrating lymphocytes (TILs), which were previously found to disturb tumor copy number and LOH (loss of heterozygosity) profiling. We specifically compared the genomic and gene expression profiles of CHEK2*1100delC breast cancers (n = 14) with BRCAX (familial non-BRCA1/BRCA2/CHEK2*1100delC mutated) breast cancers (n = 34) of the luminal intrinsic subtypes for which both SNP-array and gene expression data is available.

RESULTS

High amounts of TILs were found in a relatively small number of luminal breast cancers as compared to breast cancers of the basal-like subtype. As expected, these samples mostly have very few copy number aberrations and no detectable regions of LOH. By unsupervised hierarchical clustering of copy number data we observed a great degree of heterogeneity amongst the CHEK2*1100delC breast cancers, comparable to the BRCAX breast cancers. Furthermore, copy number aberrations were mostly seen at low frequencies in both the CHEK2*1100delC and BRCAX group of breast cancers. However, supervised class comparison identified copy number loss of chromosomal arm 1p to be associated with CHEK2*1100delC status.

CONCLUSIONS

In conclusion, in contrast to basal-like BRCA1 mutated breast cancers, no apparent specific somatic copy number aberration (CNA) profile for CHEK2*1100delC breast cancers was found. With the possible exception of copy number loss of chromosomal arm 1p in a subset of tumors, which might be involved in CHEK2 tumorigenesis. This difference in CNAs profiles might be explained by the need for BRCA1-deficient tumor cells to acquire survival factors, by for example specific copy number aberrations, to expand. Such factors may not be needed for breast tumors with a defect in a non-essential gene such as CHEK2.

Genotype/Phenotype correlations in patients with hereditary breast cancer

Of all breast cancer cases, 5-10% can be attributed to germline mutations, and the high-susceptibility genes BRCA1 and BRCA2 account for about 25-28% of these cases. For the remainder, several genes of moderate and low penetrance have been discovered. Histopathologic characteristics have been studied in small cohorts, but for most of the known non-BRCA1/2-associated hereditary breast cancers, the histologic and immunohistochemical phenotypes are not yet identified. Particularly BRCA1 tumors are associated with a distinct morphology and immunohistochemical characteristics that differ from sporadic breast cancer of age-matched controls. The recognition of features characteristic of these mutations can be helpful to identify patients likely to carry a germline mutation and to assess which gene should be screened for first, in families with a high occurrence of breast and ovarian cancer.

Sensitivity to systemic therapy for metastatic breast cancer in CHEK2 1100delC mutation carriers

PURPOSE

The role of CHEK2 in DNA repair by homologous recombination suggests that CHEK2-associated breast cancer (BC) patients might be more sensitive to chemotherapy inducing double-strand DNA breaks, but results hereon are lacking. We compared the sensitivity to first-line chemotherapy and endocrine therapy between CHEK2 1100delC and non-CHEK2 metastatic breast cancer (MBC) patients.

METHODS

Sixty-two CHEK2 1100delC MBC patients were selected from three cohorts genotyped for CHEK2 1100delC (one non-BRCA1/2 cohort and two sporadic cohorts). Controls were 62 non-CHEK2 MBC patients, matched for age at and year of primary BC diagnosis, and year of metastatic disease. Objective response rate (complete and partial response) to, and progression-free survival (PFS) and overall survival (OS) after start of first-line chemotherapy and endocrine therapy were compared between CHEK2 and non-CHEK2 patients.

RESULTS

Median age at BC diagnosis was 46 and 51 years at MBC diagnosis. First-line chemotherapy consisted of anthracycline-based chemotherapy (n = 73), taxanes (n = 16), CMF(-like) chemotherapy (n = 33) and taxane/anthracycline regimens (n = 2). CHEK2 and non-CHEK2 patients had a comparable objective response rate (44 vs. 52 %). Also, PFS and OS after start of chemotherapy were comparable between both patient groups (hazard ratio 0.91; 95 % confidence interval 0.63-1.30 and 1.03; 95 % CI 0.71-1.49, respectively). Thirty-six CHEK2 and 32 non-CHEK2 patients received first-line endocrine therapy (mainly tamoxifen) for MBC. No significant differences were observed in objective response rate to, and PFS and OS after start of endocrine therapy.

CONCLUSION

No differential efficacy of chemotherapy and endocrine therapy given for MBC was observed in CHEK2 versus non-CHEK2 patients.

Proper genomic profiling of (BRCA1-mutated) basal-like breast carcinomas requires prior removal of tumor infiltrating lymphocytes

INTRODUCTION

BRCA1-mutated breast carcinomas may have distinct biological features, suggesting the involvement of specific oncogenic pathways in tumor development. The identification of genomic aberrations characteristic for BRCA1-mutated breast carcinomas could lead to a better understanding of BRCA1-associated oncogenic events and could prove valuable in clinical testing for BRCA1-involvement in patients.

METHODS

For this purpose, genomic and gene expression profiles of basal-like BRCA1-mutated breast tumors (n = 27) were compared with basal-like familial BRCAX (non-BRCA1/2/CHEK2*1100delC) tumors (n = 14) in a familial cohort of 120 breast carcinomas.

RESULTS

Genome wide copy number profiles of the BRCA1-mutated breast carcinomas in our data appeared heterogeneous. Gene expression analyses identified varying amounts of tumor infiltrating lymphocytes (TILs) as a major cause for this heterogeneity. Indeed, selecting tumors with relative low amounts of TILs, resulted in the identification of three known but also five previously unrecognized BRCA1-associated copy number aberrations. Moreover, these aberrations occurred with high frequencies in the BRCA1-mutated tumor samples. Using these regions it was possible to discriminate BRCA1-mutated from BRCAX breast carcinomas, and they were validated in two independent cohorts. To further substantiate our findings, we used flow cytometry to isolate cancer cells from formalin-fixed, paraffin-embedded, BRCA1-mutated triple negative breast carcinomas with estimated TIL percentages of 40% and higher. Genomic profiles of sorted and unsorted fractions were compared by shallow whole genome sequencing and confirm our findings.

CONCLUSION

This study shows that genomic profiling of in particular basal-like, and thus BRCA1-mutated, breast carcinomas is severely affected by the presence of high numbers of TILs. Previous reports on genomic profiling of BRCA1-mutated breast carcinomas have largely neglected this. Therefore, our findings have direct consequences on the interpretation of published genomic data. Also, these findings could prove valuable in light of currently used genomic tools for assessing BRCA1-involvement in breast cancer patients and pathogenicity assessment of BRCA1 variants of unknown significance. The BRCA1-associated genomic aberrations identified in this study provide possible leads to a better understanding of BRCA1-associated oncogenesis.

CHEK2 c.1100delC allele is rarely identified in Greek breast cancer cases

The CHEK2 gene encodes a protein kinase that plays a crucial role in maintenance of genomic integrity and the DNA repair mechanism. CHEK2 germline mutations are associated with increased risk of breast cancer and other malignancies. From a clinical perspective, the most significant mutation identified is the c.1100delC mutation, which is associated with an approximately 25% lifetime breast cancer risk. The distribution of this mutation shows wide geographical variation; it is more prevalent in the Northern European countries and less common, or even absent, in Southern Europe. In order to estimate the frequency of the CHEK2 c.1100delC mutation in Greek breast cancer patients, we genotyped 2,449 patients (2,408 females and 41 males), which was the largest series ever tested for c.1100delC. The mean age of female and male breast cancer diagnosis was 49 and 59 years, respectively. All patients had previously tested negative for the Greek BRCA1 founder and recurrent mutations. The CHEK2 c.1100delC mutation was detected in 0.16% (4 of 2,408) of females, all of whom were diagnosed with breast cancer before the age of 50 years. Only one c.1100delC carrier was reported with breast cancer family history. The present study indicates that the CHEK2 c.1100delC mutation does not contribute substantially to hereditary breast cancer in patients of Greek descent.

Hereditary breast cancer: clinical, pathological and molecular characteristics

Pathogenic mutations in BRCA1 or BRCA2 are only detected in 25% of families with a strong history of breast cancer, though hereditary factors are expected to be involved in the remaining families with no recognized mutation. Molecular characterization is expected to provide new insight into the tumor biology to guide the search of new high-risk alleles and provide better classification of the growing number of BRCA1/2 variants of unknown significance (VUS). In this review, we provide an overview of hereditary breast cancer, its genetic background, and clinical implications, before focusing on the pathologically and molecular features associated with the disease. Recent transcriptome and genome profiling studies of tumor series from BRCA1/2 mutation carriers as well as familial non-BRCA1/2 will be discussed. Special attention is paid to its association with molecular breast cancer subtypes as well as the latest advances in predicting BRCA1/2 involvement (BRCAness) using molecular signatures, for improved diagnostics and selection of patients sensitive to targeted therapeutics.

Targeted sequencing by proximity ligation for comprehensive variant detection and local haplotyping

Despite developments in targeted gene sequencing and whole-genome analysis techniques, the robust detection of all genetic variation, including structural variants, in and around genes of interest and in an allele-specific manner remains a challenge. Here we present targeted locus amplification (TLA), a strategy to selectively amplify and sequence entire genes on the basis of the crosslinking of physically proximal sequences. We show that, unlike other targeted re-sequencing methods, TLA works without detailed prior locus information, as one or a few primer pairs are sufficient for sequencing tens to hundreds of kilobases of surrounding DNA. This enables robust detection of single nucleotide variants, structural variants and gene fusions in clinically relevant genes, including BRCA1 and BRCA2, and enables haplotyping. We show that TLA can also be used to uncover insertion sites and sequences of integrated transgenes and viruses. TLA therefore promises to be a useful method in genetic research and diagnostics when comprehensive or allele-specific genetic information is needed.

Survival and contralateral breast cancer in CHEK2 1100delC breast cancer patients: impact of adjuvant chemotherapy

Background:

We assessed the sensitivity to adjuvant chemotherapy in cell cycle checkpoint kinase 2 (CHEK2) vs non-CHEK2 breast cancer patients by comparing the contralateral breast cancer incidence and distant disease-free and breast cancer-specific survival between both groups, stratified for adjuvant chemotherapy.

Methods:

One Dutch hereditary non-BRCA1/2 breast cancer patient cohort (n=1220) and two Dutch cohorts unselected for family history (n=1014 and n=2488, respectively) were genotyped for CHEK2 1100delC. Hazard ratios for contralateral breast cancer, distant disease-free and breast cancer-specific death for mutation carriers vs noncarriers were calculated using the Cox proportional hazard method, stratified for adjuvant chemotherapy.

Results:

The CHEK2 mutation carriers (n=193) had an increased incidence of contralateral breast cancer (multivariate hazard ratio 3.97, 95% confidence interval 2.59–6.07). Distant disease-free and breast cancer-specific survival were similar in the first 6 years in mutation carriers compared with noncarriers, but diverted as of 6 years after breast cancer diagnosis (multivariate hazard ratios and 95% confidence intervals 2.65 (1.79–3.93) and 2.05 (1.41–2.99), respectively). No significant interaction between CHEK2 and adjuvant chemotherapy was observed.

Conclusions:

The CHEK2 1100delC-associated breast cancer is associated with a higher contralateral breast cancer rate as well as worse survival measures beyond 6 years after diagnosis. No differential sensitivity to adjuvant chemotherapy was observed in CHEK2 patients.

RNA profiling reveals familial aggregation of molecular subtypes in non-BRCA1/2 breast cancer families

BACKGROUND

In more than 70% of families with a strong history of breast and ovarian cancers, pathogenic mutation in BRCA1 or BRCA2 cannot be identified, even though hereditary factors are expected to be involved. It has been proposed that tumors with similar molecular phenotypes also share similar underlying pathophysiological mechanisms. In the current study, the aim was to investigate if global RNA profiling can be used to identify functional subgroups within breast tumors from families tested negative for BRCA1/2 germline mutations and how these subgroupings relate to different breast cancer patients within the same family.

METHODS

In the current study we analyzed a collection of 70 frozen breast tumor biopsies from a total of 58 families by global RNA profiling and promoter methylation analysis.

RESULTS

We show that distinct functional subgroupings, similar to the intrinsic molecular breast cancer subtypes, exist among non-BRCA1/2 breast cancers. The distribution of subtypes was markedly different from the distribution found among BRCA1/2 mutation carriers. From 11 breast cancer families, breast tumor biopsies from more than one affected family member were included in the study. Notably, in 8 of these families we found that patients from the same family shared the same tumor subtype, showing a tendency of familial aggregation of tumor subtypes (p-value = 1.7e-3). Using our previously developed BRCA1/2-signatures, we identified 7 non-BRCA1/2 tumors with a BRCA1-like molecular phenotype and provide evidence for epigenetic inactivation of BRCA1 in three of the tumors. In addition, 7 BRCA2-like tumors were found.

CONCLUSIONS

Our finding indicates involvement of hereditary factors in non-BRCA1/2 breast cancer families in which family members may carry genetic susceptibility not just to breast cancer but to a particular subtype of breast cancer. This is the first study to provide a biological link between breast cancers from family members of high-risk non-BRCA1/2 families in a systematic manner, suggesting that future genetic analysis may benefit from subgrouping families into molecularly homogeneous subtypes in order to search for new high penetrance susceptibility genes.

Classifications within molecular subtypes enables identification of BRCA1/BRCA2 mutation carriers by RNA tumor profiling

Pathogenic germline mutations in BRCA1 or BRCA2 are detected in less than one third of families with a strong history of breast cancer. It is therefore expected that mutations still remain undetected by currently used screening methods. In addition, a growing number of BRCA1/2 sequence variants of unclear pathogen significance are found in the families, constituting an increasing clinical challenge. New methods are therefore needed to improve the detection rate and aid the interpretation of the clinically uncertain variants. In this study we analyzed a series of 33 BRCA1, 22 BRCA2, and 128 sporadic tumors by RNA profiling to investigate the classification potential of RNA profiles to predict BRCA1/2 mutation status. We found that breast tumors from BRCA1 and BRCA2 mutation carriers display characteristic RNA expression patterns, allowing them to be distinguished from sporadic tumors. The majority of BRCA1 tumors were basal-like while BRCA2 tumors were mainly luminal B. Using RNA profiles, we were able to distinguish BRCA1 tumors from sporadic tumors among basal-like tumors with 83% accuracy and BRCA2 from sporadic tumors among luminal B tumors with 89% accuracy. Furthermore, subtype-specific BRCA1/2 gene signatures were successfully validated in two independent data sets with high accuracies. Although additional validation studies are required, indication of BRCA1/2 involvement (“BRCAness”) by RNA profiling could potentially be valuable as a tool for distinguishing pathogenic mutations from benign variants, for identification of undetected mutation carriers, and for selecting patients sensitive to new therapeutics such as PARP inhibitors.

CHEK2*1100delC homozygosity in the Netherlands--prevalence and risk of breast and lung cancer

The 1100delC mutation in the CHEK2 gene has a carrier frequency of up to 1.5% in individuals from North-West Europe. Women heterozygous for 1100delC have an increased breast cancer risk (odds ratio 2.7). To explore the prevalence and clinical consequences of 1100delC homozygosity in the Netherlands, we genotyped a sporadic breast cancer hospital-based cohort, a group of non-BRCA1/2 breast cancer families, and breast tumors from a tumor tissue bank. Three 1100delC homozygous patients were found in the cohort of 1434 sporadic breast cancer patients, suggesting an increased breast cancer risk for 1100delC homozygotes (odds ratio 3.4, 95% confidence interval 0.4-32.6, P=0.3). Another 1100delC homozygote was found in 592 individuals from 108 non-BRCA1/2 breast cancer families, and two more were found after testing 1706 breast tumors and confirming homozygosity on their wild-type DNA. Follow-up data was available for five homozygous patients, and remarkably, three of them had developed contralateral breast cancer. A possible relationship between 1100delC and lung cancer risk was investigated in 457 unrelated lung cancer patients but could not be confirmed. Due to the small number of 1100delC homozygotes identified, the breast cancer risk estimate associated with this genotype had limited accuracy but is probably higher than the risk in heterozygous females. Screening for CHEK2 1100delC could be beneficial in countries with a relatively high allele frequency.