Skip to main content

An immune-centric exploration of BRCA1 and BRCA2 germline mutation related breast and ovarian cancers

Abstract

Background

BRCA1/2 germline mutation related cancers are candidates for new immune therapeutic interventions. This study was a hypothesis generating exploration of genomic data collected at diagnosis for 19 patients. The prominent tumor mutation burden (TMB) in hereditary breast and ovarian cancers in this cohort was not correlated with high global immune activity in their microenvironments. More information is needed about the relationship between genomic instability, phenotypes and immune microenvironments of these hereditary tumors in order to find appropriate markers of immune activity and the most effective anticancer immune strategies.

Methods

Mining and statistical analyses of the original DNA and RNA sequencing data and The Cancer Genome Atlas data were performed. To interpret the data, we have used published literature and web available resources such as Gene Ontology, The Cancer immunome Atlas and the Cancer Research Institute iAtlas.

Results

We found that BRCA1/2 germline related breast and ovarian cancers do not represent a unique phenotypic identity, but they express a range of phenotypes similar to sporadic cancers. All breast and ovarian BRCA1/2 related tumors are characterized by high homologous recombination deficiency (HRD) and low aneuploidy. Interestingly, all sporadic high grade serous ovarian cancers (HGSOC) and most of the subtypes of triple negative breast cancers (TNBC) also express a high degree of HRD.

Conclusions

TMB is not associated with the magnitude of the immune response in hereditary BRCA1/2 related breast and ovarian cancers or in sporadic TNBC and sporadic HGSOC. Hereditary tumors express phenotypes as heterogenous as sporadic tumors with various degree of “BRCAness” and various characteristics of the immune microenvironments. The subtyping criteria developed for sporadic tumors can be applied for the classification of hereditary tumors and possibly also characterization of their immune microenvironment. A high HRD score may be a good candidate biomarker for response to platinum, and potentially PARP-inhibition.

Trial registration

Phase I Study of the Oral PI3kinase Inhibitor BKM120 or BYL719 and the Oral PARP Inhibitor Olaparib in Patients With Recurrent TNBC or HGSOC (NCT01623349), first posted on June 20, 2012. The design and the outcome of the clinical trial is not in the scope of this study.

Peer Review reports

Background

The concept of cancer immunosurveillance, which claimed that the immune system can protect the host against the development of cancer, was proposed over 50 years ago by Burnet and Thomas [1, 2]. Recently, the evidence in favor of cancer immunosurveillance has been translated into new therapeutic approaches. DNA damage and genomic instability are closely linked to immunity. The production of tumor specific neoantigens is believed to be triggered by various mutations in the unstable cancer genome. Thus, immunosurveillance should be particularly relevant to BRCA1/2 germline mutation carriers, whose tumors have dysfunctional homologous recombination (HR), the main pathway for DNA double strand break repair [3].

The HR deficiency of hereditary breast and ovarian cancers makes them vulnerable to the inhibition of alternative pathways of DNA repair with inhibitors of Poly (ADP-Ribose) Polymerase (PARP) [4]. There are interests in expanding the use of PARP inhibitors to sporadic breast and ovarian cancers, some of which express phenotypes similar to hereditary tumors. For example, many sporadic TNBCs show deficiency in HR and demonstrate “BRCA-like” clinicopathological features, often referred to as “BRCAness” [5, 6]. “BRCAness” phenotype is also attributed to many hereditary and sporadic HGSOCs. However, the “BRCAness” phenotype is still poorly defined [7].

Due to having high TMB, BRCA1/2 germline mutation related tumors are considered to be candidates for immune checkpoint inhibition strategies, which were successful in highly mutated melanoma and lung cancers [8]. However, it has been shown that the BRCA1 gene product is a versatile regulator involved in many cellular functions in addition to its role in the DNA repair [9]. Moreover, the BRCA1 and BRCA2 gene products contribute in different ways to the tumorigenesis [10]. To find effective immune therapeutic strategies against hereditary breast and ovarian cancers, more information is needed about the relationship between genomic instability, phenotypes and immune microenvironments of those tumors.

The goal of this study was to explore genomic instability and phenotypes of hereditary and sporadic breast and ovarian cancers in relation to their immune microenvironments. Our results may help to find appropriate ways to stratify those tumors for testing various immune interventions. They will also help clarify the differences and similarities between BRCA1/2 germline mutation related phenotypes versus sporadic phenotypes of TNBC and HGSOC, and will help to define more precisely the elusive “BRCAness” phenotype.

Methods

Patients

The patients contributed to this study were selected for a clinical trial (#NCT01623349).

The genomic data was acquired from 19 patients out of total of 118 enrolled in the trial. Genetic material was extracted from Formalin-Fixed Paraffin-Embedded (FFPE) blocks prepared from tumors at diagnosis, before any treatment was administered to the patients. Eventually, all the 19 patients were treated heavily with conventional chemotherapy and fail the treatments. Details about the line of treatments are shown in Additional file 1: Table S1. This information may be relevant since it suggests that all the patients in this cohort could be considered resistant to conventional therapy. Design of a subsequent trial and the outcome of the trial are not in the scope of this hypothesis generating study and are available elsewhere (https://clinicaltrials.gov/ct2/show/NCT01623349).

The cohort was enriched in BRCA1/2 germline mutation carriers. The BRCA1/2 germline status was determined by a clinical test: MKS IMPACT™ tumor-profiling multiplex panel [11]. BRCA1/2 proteins were expressed in all samples, as determined at RNA level (data not shown).

RNA sequencing

RNA was extracted from Formalin-Fixed Paraffin-Embedded (FFPE) samples.

Qiagen RNeasy FFPE kit was used to extract RNA. TruSeq RNA and Access library prep kit was then used for preparing library for IIlumina RNA sequencing.

Illumina Sequencing: Illumina NextSeq 500 High Output v2 sequencer has been used to generate sequences in the FASTQ format. The 150-cycle kit for paired end 2 × 75 bp sequencing has been used with estimated 60 million total paired end raw reads per sample.

Sample extraction, library preparation and sequencing were done at the Center for Genome Innovation (CGI), Institute for System Genomics, University of Connecticut.

RNA-seq data analysis

Quality Check: FASTQ file quality was checked using FASTQC v0.11.2. The summary reports showed no potential errors or warnings.

Alignment and Pre-processing: Reads were mapped using STAR Aligner tool v020201 to the human genome reference (hg19) downloaded from UCSC genome browser.

Transcripts quantification: Gene expression levels were obtained from the RNA-seq dataset using RSEM v1.2.31 with Ensembl gene annotation database.

Differential expression analysis: we have used EBSeq v1.21.0 for differential gene expression analysis of the RNA-seq data.

Whole exome sequencing

FFPE samples were used for extracting DNA. Whole exome sequencing has been done at Memorial Sloan Kettering Cancer Center using Illumina sequencers. FASTQC v0.11.2 was used to check the quality of the paired end raw sequencing data in FASTQ format. The summary reports showed no potential errors or warnings.

Reads were aligned to hg19 genome reference using BWA v0.7.12-r1039 mem software tools.

Subtyping breast and ovarian tumors

TNBC clinical trial samples were subtyped according to Lehmann et al. [12] into 6 subtypes, using their TNBC type tool run on genome-wide gene expression matrices for each sample [13], (http://cbc.mc.vanderbilt.edu/tnbc). Ovarian clinical trial samples were subtyped using the Classification of Ovarian Cancer (CLOVAR) scheme proposed by Verhaak et al. [14]. They defined a gene signature- set of 100 genes, used for classifying ovarian cancer into four subtypes. Single sample gene set enrichment analysis (SSGSEA) [15] was performed on each sample using these CLOVAR gene set. For every sample, SSGSEA outputs a score for each of the four subtypes. The highest score defines the classification for that sample. TNBC and CLOVAR subtypes for the The Cancer Genome Atlas (TCGA) dataset were downloaded from Lehmann et al. and Verhaak at al., respectively [14, 16]. Immune Subtyping on the clinical trial samples was performed using the Immune Subtype Classifier available from The Cancer Research Institutes iAtlas (https://www.cri-iatlas.org/about/). Immune Subtypes for TCGA data were download from iAtlas.

Mutation burden analysis

Each patient’s tumor and normal BAM files were input into samtools v1.7 mpileup. Varscan somatic was called on each mpileup file yielding unfiltered vcf files. Varscan processSomatic was used to isolate high confidence SNV and indel calls, which were then false positive filtered using bam-readcount v0.8.0 and Varscan FPfilter. These high confidence, false positive filtered vcf files were used for analysis.

TCGA Mutation Annotation files for breast and ovarian cancer were downloaded from FireBrowse data version 2016_01_28 (firebrowse.org/).

Leukocyte fraction and homologous recombination

Breast and Ovarian Leukocyte fraction and Homologous Recombination data was downloaded from iAtlas data portal (https://www.cri-iatlas.org/about/).

Statistics

All statistical analysis was carrier out in R. Statistical significance was defined at a p-value < 0.05: **** < 0.0001, *** < 0.001, ** < 0.01, * < 0.05, measured by nonparametric Wilcoxon test, unless otherwise specified.

Results

Breast and ovarian cancers in BRCA1/2 germline mutation carriers show relatively low overall immune activity at diagnosis, compared to very immune active non-carriers

In our clinical trial samples, we observed a striking difference in the gene expression profiles between germline mutation carriers and non-carriers. There were 1308 genes differentially expressed between carriers and non-carriers (Posterior Probability of equal expression < 0.05). Of these, 813 showed significantly higher expression in non-carriers (log fold change > 1.5) The biological processes most highly enriched in non-carriers identified with Gene ontology tool (Panther Classification System: http://www.pantherdb.org) were all related to immune functions (Fig. 1). Other biological processes which were also enriched in non-carriers include calcium ion transport and signaling, regulation of cell adhesion, motility and chemotaxis, protein secretion, cell signaling (MAPK, ERK1/2 and JNK), cell proliferation, differentiation and cell death (Additional file 2: Table S2). Many of these processes are related to biology of immune cells. Genes overexpressed in carriers, on the other hand, were not enriched for any particular biological process (data not shown).

Fig. 1
figure1

Biological processes enriched in breast and ovarian non-carriers from the clinical trial. The list of 813 genes was analyzed with Panther classification system (http://www.pantherdb.org). The table shows the top most significantly enriched biological process. The complete list of enriched processes is shown in Additional file 2: Table S2

We have focused on the 500 biological processes, highly enriched in non-carriers, which were related to immune functions such as T cells differentiation and selection, B cells activation and regulation, production of various Interleukins and signaling via TNF alpha and interferon gamma. This data was highly significant suggesting that the immune environment of sporadic breast and ovarian cancers in our cohort was much more active relative to that of carriers of germline mutations in BRCA1/2 genes. This was independent from the type of germline mutation, BRCA1 or BRCA2 (Additional file 3: Table S3 and Additional file 4: Table S4) and was true for both types of cancers when analyzed independently (Additional file 5: Table S5 and Additional file 6: Table S6). Many genes overexpressed in breast non-carriers overlapped with those overexpressed in ovarian non-carriers (60 genes). The commonly upregulated genes in breast and ovarian non-carriers were all involved in immune functions (Fig. 2 and Additional file 7: Table S7).

Fig. 2
figure2

The common genes upregulated in breast and ovarian non-carriers from the clinical trial are involved in immune functions. 60 genes overexpressed in breast non-carriers overlapped with those overexpressed in ovarian non-carriers. The list of 60 genes was analyzed with Panther classification system (http://www.pantherdb.org). The table shows the top most significantly enriched biological process. The complete list of processes is shown in Additional file 7: Table S7

Recently there have been attempts to characterize the immune components of the tumor microenvironment from high-throughput expression data [17,18,19,20,21]. The most complete analysis of immune infiltrates in tumor microenvironment was performed by the group of Trojanoski [21]. They developed a comprehensive and interactive database for immunogenomic studies: The Cancer Immunome Atlas (TCIA) (https://tcia.at/home), which allows exploration of specific immune related gene sets and assessment of cellular composition of infiltrates from 20 solid cancers. We have used their list of 782 genes, which characterize 28 different cell types present in tumor infiltrates [22] to analyze the global immune landscapes of individual carriers and non-carriers in our cohort (Fig. 3). The gene list is shown in Additional file 8: Table S8. All four breast carriers of germline BRCA1/2 mutation showed overall low expression of genes associated with various immune cell types, while three non-carriers showed relatively high expression of most of those genes. The picture was different for ovarian cancers, where some carriers and some non-carriers showed various expression of immune genes consistent with less robust differential expression results. Thus, the expression of 28 meta-gene sets validated our results obtained from differential expression analysis. Expression of these meta-gene sets can be a convenient way of representing global immune activity of tumors.

Fig. 3
figure3

Patterns of expression of 782 genes representing 28 immune cell types, in samples from the clinical trial. Heat-maps represent expression of 782 genes in breast (a) and ovarian (b) samples from our cohort. The 782 gene list is shown in Additional file 8: Table S8. Breast carriers (BC), breast non-carriers (BN), ovarian carriers (OC), ovarian non-carriers (ON). The numbers correspond to the patient number (Fig. 4a)

BRCA1/2 germline mutation related breast and ovarian cancers show a range of phenotypes similar to that of sporadic cancers

There is still controversy if hereditary BRCA1/2 mutation related tumors represent a separate phenotypic identity. Both TNBC and HGSOC represent heterogenous groups of cancers and recently both tumor types were subdivided into several subtypes [12, 14, 16, 23,24,25,26,27]. Six subtypes of TNBC (IM, BL1, BL2, LAR, M and MSL) were identified from clustering of gene expression data [12]. The Immunomodulatory (IM) subtype is enriched in immune cell signaling. Two other subtypes (basal-like 1 and basal-like 2 (BL1 and BL2)) express high levels of the genes involved in cell proliferation and DNA damage response (DDR), however BL2 is of basal myoepithelial origin and can be distinguished by activated signaling pathways (EGF, NGF, MET, Wnt/β catenin and IGFR1) and glycolysis. Luminal androgen receptor (LAR) subtype is the most distinct of all subtypes, characterized by luminal features and expression of androgen receptor. Mesenchymal (M) and mesenchymal-stem like (MSL) subtypes are characterized by expression of genes involved in epithelial/mesenchymal transition. Patients with BL1 tumors show relatively good prognosis, while patients with BL2 tumors have very poor outcome [28].

Four subtypes of HGSOC (IMR, DIF, MES and PRO) were identified by gene expression profiling. The immunoreactive subtype (IMR) is enriched in immune cell signature, the differentiated subtype (DIF) expresses differentiation markers, the mesenchymal subtype (MES) is characterized by stromal expression signature indicating activated stroma, while the proliferative (PRO) subtype is characterized by low expression of ovarian cancer markers, but overexpression of proliferation and extracellular matrix (ECM) related genes. Importantly, the expression clusters distinguishing the subtypes strongly correlate with histological types of HGSOC [25]. Among all subtypes, the IMR shows the best prognosis and MES subtype has relatively poor outcome [14].

Only one of six subtypes of TNBC (IM) and one of four subtypes of HGSOC (IMR) are characterized by a highly immune active microenvironment. We used a publicly available tool for TNBC classification developed by Lehmann to classify breast tumors from the clinical trial samples [13], ( http://cbc.mc.vanderbilt.edu/tnbc). The classification of HGSOC was obtained using the CLOVAR signature (see Methods section for details). Indeed, one of the three sporadic TNBC in this cohort was immunomodulatory, while two others belonged to different categories (MSL and BL2) (Fig. 4a). Breast tumors from BRCA1/2 germline mutation carriers expressed M and LAR subtypes and none were immunomodulatory. Interestingly, two BRCA2 germline mutation related breast tumors were classified as not TNBC. Most of the HGSOC from carriers and non-carriers of germline mutations belonged to MES subtype and none were immunomodulatory. Thus, none of the patients in this cohort, who carried germline mutation in BRCA1/2, developed highly immune-active tumors at diagnosis (Fig. 4a). In addition, none of the TNBC were classified as BL1, which is associated with good prognosis and the majority of HGSOC (70%) expressed MES subtype associated with the poor prognosis. This is consistent with the history of the patients in this cohort (lack of response to conventional therapies and progression to metastasis).

Fig. 4
figure4

Subtypes of hereditary and sporadic breast and ovarian cancers in the clinical trial and in TCGA database. a) List of clinical trial samples. Subtyping of tumors from this cohort was obtained using TNBCtype tool (http://cbc.mc.vanderbilt.edu/tnbc) for breast cancers and CLOVAR scheme [14] for ovarian cancers. b) List of hereditary breast tumors from TCGA. Subtyping of these tumors was acquired from Lehmann at al [16]. c) Distribution of TNBC subtypes within TCGA breast cancers (sporadic TNBC and hereditary BRCA1 and BRCA2 related breast tumors). d) Distribution of HGSOC subtypes within TCGA ovarian cancers (sporadic HGSOC and hereditary BRCA1 and BRCA2 related ovarian tumors). The list of breast germline mutation carriers was established according the information acquired from CBioPortal (http://www.cbioportal.org) and iAtlas https://www.cri-iatlas.org/about/. The list of ovarian germline mutation carriers was established from CBioPortal (http://www.cbioportal.org) and it is shown in Additional file 9: Table S9. Immune Subtypes for our cohort were identified using tool available in iAtlas interactive platform and for TCGA samples were download from the site

To put this data into perspective we examined the classification of all BRCA1/2 germline mutation related breast and ovarian tumors from The Cancer Genome Atlas (TCGA) datasets (Fig. 4b and Additional file 9: Table S9). Consistent with the subtyping in our clinical trial samples, few BRCA1 germline mutation related breast tumors in TCGA database are immunomodulatory (7% versus 21% of TNBC from non-carriers) and most BRCA2 germline mutation related breast cancers do not classify as TNBC (12 out of 15, 80%) (Fig. 4b and c). The results for ovarian cancers show a similar pattern. However, it is important to emphasize that HGSOC often express multiple signatures. Therefore, classification into mutually exclusive subtypes may be less specific than in other cancers [14]. Nevertheless, BRCA1 /2 germline mutation related HGSOC are not enriched in immunoreactive phenotype (Fig. 4d).

Thus, indeed BRCA1/2 germline mutation related tumors do not belong to the most immune active category of breast and ovarian cancers. The data also suggest that BRCA1/2 germline mutation related breast and ovarian cancers express range of phenotypes similar to sporadic cancers and therefore it is unlikely that they represent unique phenotypic identity within TNBC or HGSOC.

However, BRCA1/2 hereditary tumors have unique mutational signature [29] and BRCA 1 tumors have characteristic genomic copy number alterations [30]. Thus, it seems that mostly genotypes, but not phenotypes, make tumors related to BRCA1/2 germline mutation carriers unique.

BRCA1/2 germline mutation related breast and ovarian cancers show relatively low overall immune activity in their microenvironment despite having elevated mutation burden

The relatively low immune activity in cancers (breast and ovarian) from BRCA1/2 germline mutation carriers is counterintuitive. Tumors with compromised DNA repair usually have a high mutational load and would be expected to generate a high number of neo-antigens [31]. In addition, hypermutated cancers such as melanoma or lung cancer, as well as colon cancer deficient in mismatch repair show positive response to immunotherapy [32,33,34].

As expected, germline mutation carriers from our clinical trial samples show a higher tumor mutational burden (TMB) compared to non-carriers (Fig. 5a) and this is in contrast to global immune activity, which is lower in mutation carriers (Fig. 5b). Thus, we asked if there is a correlation between TMB and global immune activity in TCGA.

Fig. 5
figure5

Hereditary breast and ovarian cancers from the clinical trial and from TCGA database show high TMB and low overall immune activity relative to the sporadic tumors. Data obtained for our cohort (a, b), data acquired for TCGA breast (c-e) and ovarian (f-h) cancers. c and f) Somatic mutation count acquired from CBioPortal (http://www.cbioportal.org), d and g) Global immune gene expression representing averaged expression of genes from 28 meta-gene sets. Expression data was downloaded from FireBrowse data version 2016_01_28 (this link http://firebrowse.org/) e and h) Leukocyte fraction acquired from Cancer Research Institute iAtlas https://www.cri-iatlas.org/about/. The dotted lines indicate the average value for all the samples in each panel

Within breast cancers, TMB was higher for BRCA1 and BRCA2 germline mutation carriers relative to non-carriers and was also elevated in BL1 subtype. Within HGSOC, TMB was higher only for germline mutation carriers and did not vary among other subtypes (Fig. 5c and f). Remarkably however, the global immune activity of tumor microenvironments, calculated as averaged expression of genes from 28 meta-gene sets, varied widely between subtypes (Fig. 5d and g).

Another measure of global immune activity is the leukocyte fraction of tumors. The leukocyte fraction for samples from TCGA is available on the web-based interactive platform: the Cancer Research Institute iAtlas https://www.cri-iatlas.org/about/. iAtlas was designed from extensive immunogenomic analysis and integration of the data for 33 cancer types [35]. The leukocyte fractions in subtypes of hereditary and sporadic TNBC and HGSOC from the TCGA database showed a very similar pattern to the expression of 28 meta-gene sets (Fig. 5e and h) and also did not correlate with TMB. Thus, BRCA 1/2 germline mutation related hereditary breast and ovarian tumors, have low overall immune activity within their tumor microenvironments despite their elevated TMB. The data suggest that diversity of immune responses in the microenvironments of hereditary and sporadic TNBC and HGSOC is likely determined by factors other than TMB.

Pattern of genomic instability is different in BRCA1 versus BRCA2 germline related tumors

TNBC and HGSOC are characterized by frequent mutations in TP53 gene and a high degree of genomic instability. Considering that elevated TMB in hereditary breast and ovarian cancers was not associated with high immune activity in the tumor microenvironment, we looked at other measures of instability that potentially could influence immune response in breast and ovarian cancers. Recently, the extensive Pan-Cancer analysis of DNA damage repair (DDR) deficiencies in cancer was published [36] and the results were made available in iAtlas (https://www.cri-iatlas.org/about/). Using this resource, we explored several measures of genomic instability including: mutation load (expressed as non-silent mutation rate and SNV neoantigen count), CNV load (expressed as number of segments and fraction genome altered), aneuploidy and HR deficiency. Genomic instability varies widely between the subtypes of breast and ovarian cancers. As expected, all tumors from germline mutation carriers display high HR deficiency but also relatively low aneuploidy. Consistent with the results shown in Fig. 5, breast and ovarian cancers from germline mutation carriers have a relatively high mutation load compared to non-carriers. (Fig. 6a and b). BRCA2 related tumors reveal a very different pattern of instability compared to BRCA1 germline related tumors with a low CNV load. This confirms that the characteristic copy number pattern published earlier for hereditary breast cancers applies only to BRCA1-related tumors [30, 37]. The relationship between measures of genomic instability and the immune activity in tumors may be complex and require further investigation.

Fig. 6
figure6

Pattern of genomic instability vary widely within hereditary and sporadic breast and ovarian cancers and it is different in BRCA1 versus BRCA2 germline related tumors. Heat-maps represent genomic instability measures in breast (a) and ovarian (b) cancers from TCGA. The data was acquired from Cancer Research Institute iAtlas (https://www.cri-iatlas.org/about/)

High HR deficiency score characterize all BRCA1/2 germline mutation carriers and is predictive of response to platinum in HGSOC

HR deficiency is particularly relevant for hereditary TNBC and HGSOC. Ovarian cancer has the highest HR deficiency score of all 33 cancers included in TCGA (average value > 40) while breast cancers show much lower HR deficiency score (average value > 20) (Fig. 7a) [36]. However, TNBCs show a HR deficiency score as high as ovarian cancers (average value > 40), with the only exception of the LAR subtype (Fig. 7b). As expected, breast and ovarian tumors from BRCA1/2 germline mutation carriers have even higher HR deficiency scores (average value > 50 for BRCA2 and > 60 for BRCA1 mutation carriers) (Fig. 7b and c). Similar to TMB, HR deficiency did not correlate with immune activity. However, HR deficiency in ovarian cancers did correlate with platinum sensitivity (Fig. 7d). The sensitive and resistant ovarian cancers were selected from TCGA database. Tumors were defined as sensitive if there was no evidence of progression or recurrence at least 6 months from the date of primary platinum treatment. Tumors that recurred within 6 months of primary treatment were considered resistant [27]. The ovarian cancers sensitive to platinum had average HR deficiency score of 46.5 and resistant tumors had the score of 36.4. The difference was statistically significant.

Fig. 7
figure7

High HR deficiency score characterize most of TNBC and predicts platinum sensitivity in HGSOC. a) Distribution of HR deficiency score across 33 TCGA cancer types, b) across the breast cancer subtypes and c) across the HGSOC subtypes. The data was acquired from Cancer Research Institute iAtlas (https://www.cri-iatlas.org/about/). d) HR deficiency score in HGSOC, which are resistant or sensitive to platinum-based therapy. The sensitivity/resistance criteria were established according to Integrated genomic analysis of ovarian carcinoma [27] and applied to TCGA data (Additional file 10: Table S10). The dotted lines indicate mean HR deficiency score for all HGSOC (top line) and all breast cancers (bottom line)

Distribution of “BRCAness” in subtypes of breast and ovarian cancers

The term “BRCAness” phenotype was coined to describe sporadic breast and ovarian cancers that behave like hereditary BRCA1/2-related tumors [5, 7].

The “BRCAness” characteristics of the subtypes of breast and HGSOC including BRCA1/2 germline mutation carriers from TCGA database are presented in Tables 1 and 2. The most important aspects of “BRCAness” phenotype chosen from literature were as follows: deficiency in HR, high genomic instability, frequent P53 mutations, but infrequent PI3K mutations in breast and ovarian cancers, in addition to basal like classification and high probability of pathological complete response (pCR) in breast cancers [6, 7, 28, 38, 39, 41, 42]. “BRCAness” is most often found in the BL1 and M subtypes of TNBC. Consistent with these results, most of the BRCA1 germline mutation carriers belong to BL1 or M subtype (Fig. 4c) and the “BRCA1-like” tumors selected according to copy number criteria also belong mostly to the BL1 and M category [30]. The LAR subtype, on the other hand, has frequent PIK3CA mutations and a low HR deficiency score. The IM subtype does not meet genomic instability criteria, MSL is not basal type and BL2 subtype is characterized by very low pCR. Importantly, BRCA2 germline related tumors do not express any attributes of “BRCAness” except high genomic instability.

Table 1 The “BRCAness” characteristics of breast tumors from TCGA database
Table 2 The “BRCAness” characteristic of ovarian tumors from TCGA database

Similar analysis was performed for HGSOC subtypes (Table 2). According to our criteria, all subtypes of HGSOC score high on “BRCAness”.

PD-L1 expression reflects overall magnitude of the immune response in breast and ovarian cancers

PD-L1 is the target for anti-PD-L1 antibodies, which are currently being examined in a phase II clinical trial (NCT02849496). PD-L1 RNA expression was significantly higher in samples from non-carriers of germline mutations compared to the carriers in our clinical trial samples (#NCT01623349) (Fig. 8a). Thus, higher overall immune activity corresponded with higher expression of this marker. We verified the expression of this marker in all subtypes of TNBC and HGSOC from TCGA database. All tumors expressed the protein, and the pattern of expression followed the pattern of overall immune activity in all samples including those from BRCA1/2 germline mutation carriers (see also Fig. 5c-h). The IMR subtype of HGSOC had the highest expression of PD-L1.

Fig. 8
figure8

The immune response pattern in hereditary and sporadic breast and ovarian cancers from the clinical trial and from TCGA database. Expression of PD-L1 (CD274) in carriers and non-carriers from our cohort (a) in breast tumors from TCGA (b) and in HGSOC from TCGA (c). Heat-map showing relative contribution (percentage) of six universal intratumor immune states (C1-C6) within microenvironments of tumors from our cohort and from hereditary and sporadic breast and ovarian cancers form TCGA (d). PD-L1 for TCGA samples was downloaded from The Cancer Immunome Atlas (TCIA) (https://tcia.at/home). Immune Subtyping on the clinical trial samples was performed using the Immune Subtype Classifier available from The Cancer Research Institutes iAtlas (https://www.cri-iatlas.org/about/) and is also shown in Fig. 4a. Immune Subtypes for TCGA data were download from iAtlas

The immune response patterns in TNBC and HGSOC

The immune landscape of 33 cancer types was recently published and made available on the web-based interactive platform [35], Cancer Research Institute iAtlas https://www.cri-iatlas.org/about/). They identified six universal intratumor immune states or response patterns. Briefly, C1, wound healing subtype, have elevated expression of angiogenic genes and high proliferative rate, C2, INF-γ subtype, have the highest M1/M2 macrophage polarization, C3 is an inflammatory subtype, C4 is lymphocytes depleted type displaying a more prominent macrophage signature, C5 is an immunologically quiet type and exhibit the lowest lymphocyte and highest macrophage response dominated by M2 and finally C6 is a TGF-β dominant type. When we applied the signatures for intratumor immune types (C1-C6) to our clinical trial samples, we found that the majority of non-carriers expressed C3 (inflammatory subtype), while majority of carriers expressed C1 (wound healing) subtype (Fig. 4a and Fig. 8d). The composition of the immune microenvironments within TNBC and HGSOC from TCGA varied widely, but almost universally the predominant subtypes were C2 (INF-γ and macrophage-enriched) and C1(wound healing). Some HGSOC expressed also C4 (lymphocytes depleted) subtype. Interestingly, two the most “BRCAness” expressing TNBC showed very different immune environments. BL1 tumors with higher overall immune activity relative to M tumors are predominantly (82.8%) associated with macrophage-enriched (C2) immune signature, while M tumors, which have overall very low immunoactivity, are predominantly (77.8%) associated with wound healing (C1) signature (Fig. 8d).

Discussion

Cells that carry BRCA1/2 germline mutations have a high degree of genomic instability due to dysfunctional HR repair mechanisms and consequently a high TMB. High TMB has been associated with immunogenicity and response to immune checkpoint inhibitors such PD-1/PD-L1 antibodies in melanoma and lung cancer [8]. We found that breast and ovarian cancers from the clinical trial #NCT01623349, which carry BRCA1/2 germline mutation, as well as germline mutation carriers from TCGA database had a high TMB as reported previously [29, 43, 44] but did not score high on overall immune activity and PD-L1 expression relative to non-carriers. These data are consistent with recently published larger cohorts [45] who also reported no increase in PD-L1 expression or TILs in BRCA1-like tumors. However, it has to be kept in mind that the TMB in breast and ovarian cancers, even that from BRCA1 germline mutation carriers, is an order of magnitude lower compared to hypermutated cancers such as melanoma and lung cancer (The Cancer Immunome Atlas (TCIA) https://tcia.at/home). A lack of a positive correlation between TMB and immune infiltration in various types of tumors was recently described by others [46]. There are at least two possible caveats for the evaluation of TMB from genomic data. For one, the sequencing data is prone to inconsistencies due to various ways of processing NGS data and diverse criteria for relevant mutations [47]. Secondly, DNA sequencing data for highly immune infiltrated tumors can be affected by the immune cell “contamination”, reducing the readings for genomic instability.

The breast and ovarian cancers from BRCA1/2 germline mutation carriers show a unique responsiveness to PARP-inhibition and it was suggested that they express distinctive phenotype, which they share with some sporadic breast and ovarian cancers [7, 48]. Hence, the idea of “BRCAness”. We have addressed the extent of “BRCAness” phenotype in hereditary and sporadic TNBC and HGSOC from TCGA database. According to our results, “BRCAness” is not a unique phenotype of BRCA1/2 mutation carriers, but rather an attribute of majority of HGSOC and of a few subtypes of sporadic TNBC, which also happened to be the most frequent subtypes found within hereditary BRCA1 germline mutation carriers (BL1 and M). BRCA2 germline mutation related breast cancers on the other hand, do not express any “BRCAness” features except high genomic instability, and express the range of PAM50 phenotypes, similar to all sporadic breast cancers [40]. Consequently, the subtyping criteria developed for sporadic tumors can be applied for identifying “BRCAness” in sporadic and germline mutation associated tumors.

While BRCA1/2 related hereditary tumors may not have a unique phenotype, the breast cancers that carry BRCA1/2 deficiency have a unique genotype characterized by distinctive mutation profile [29, 30]. Distinct copy number signature (“BRCA1-like”) is also shared between cancers related to germline mutation in BRCA1 gene and sporadic cancers whose BRCA1 protein was inactivated through other mechanisms [49, 50]. Remarkably, the “BRCA1-like” subgroup distinguished with the copy number criteria had down regulated expression of proteins related to immune functions and was associated mostly with BL1 and M subtypes of Lehmann [30].

The global immune activity varies widely between breast and ovarian cancer subtypes and the immune microenvironments are heterogenous. Especially interesting finding is that two subtypes of TNBC (BL1 and M), which score high on “BRCAness” and are the most prevalent in hereditary BRCA1 germline mutation carriers have fundamentally different immune profiles, predominantly C2 and C1, respectively. The data suggest that the diversity of immune responses in the microenvironments of hereditary and sporadic TNBC and HGSOC may be associated with their particular phenotypes. In this respect, the subtyping of breast and ovarian tumors according to criteria developed for sporadic tumors may also be useful for testing various immune interventions. In our clinical trial samples, the majority of non-carriers expressed at diagnosis the C3 inflammatory immune profile, defined by elevated Th17 and Th1 genes, while majority of carriers expressed C1 wound healing profile. Further studies, with bigger cohorts, are needed to confirm this finding and to explore its significance.

In iAtlas, germline mutation carriers had the highest degree of HR deficiency of all breast and ovarian cancers. However, the surprising finding was that all ovarian cancers and almost all TNBC (with exception of LAR subtype), but not ER+ or HER2+ breast cancers, had a high degree of HR deficiency as well. This could explain the sensitivity of those tumors (TNBC and HGSOC) to DNA damaging chemotherapy and PARP inhibition. Indeed, we have shown that the platinum sensitive HGSOC from TCGA had significantly higher HR deficiency compared to resistant tumors. Similar results were obtained by Telli et al. Their study included both TNBC and ovarian cancers and used combined HR deficiency score, defined in a similar way as the HR score provided in iAtlas [51]. They suggested a clinical application of the score to identify TNBC (not deficient in BRCA1/2), which likely respond to platinum. If we apply their criteria for response, (HRD score ≥ 42), only two subtypes (BL1 and M) of sporadic TNBC and most of the subtypes of HGSOC will qualified as possibly sensitive to the therapy. Interestingly, the ovarian cancers from TCGA database, sensitive to platinum had HRD score above 42 (46.5) and resistant tumors had the score below 42 (36.4). Thus, our results validated the HR deficiency score as candidate biomarker for resistance to platinum.

BRCA1/2 germline related breast and ovarian cancers have the highest HR deficiency score of 33 tumor types analyzed in iAtlas and rather low mutational burden relative to hypermutated tumors (Cancer Research Institute iAtlas https://www.cri-iatlas.org/about/). We have found that they also show a low aneuploidy score relative to sporadic tumors. The low CNV load was also evident, but only in BRCA2 related breast and ovarian cancers. Genomic instability is the outcome of few processes going on simultaneously: DNA damage, DNA repair and immunoediting that is the clearance of cells with unrepaired damage by the immune system [52]. In BRCA1/2 germline mutation related tumors relatively low TMB (compared to hypermutated tumors) may be explained by effective immunoediting, rather than by effective DNA repair. It is tempting to speculate that immunoediting can compensate for the lack of adequate repair. If so, the immunoediting will have particular impact on hereditary tumors or sporadic tumors with BRCA1/2 dysfunction [53]. In support of this, low aneuploidy in BRCA1/2 germline related breast and ovarian cancers may also suggest more active immunosurveillence against cancer associated hyperploidy [54,55,56]. Active immunoediting in tumors from BRCA1/2 germline mutation carries may be consistent with a positive response to immune therapies in TNBC tumors despite the fact that they have relatively low TMB [57].

Looking from immune-centric point of view on patients with BRCA1/2 germline mutations, one wonders how BRCA1/2 deficiency influences the systemic immunity independent of its role in breast and ovarian tumorigenesis? Indeed, B-cell differentiation and maturation requires DDR [58] and BRCA1 protein have a direct role in B cells lymphomagenesis [59] suggesting the possible alterations in the systemic immunity of germline mutation carriers. In keeping with this, BRCA1/2 germline mutation carriers have higher risk of developing certain leukemias and lymphomas [60, 61].

The possibility of systemic immunity playing a role in tumorigenesis in the carriers of BRCA1/2 germline mutations was already suggested by others. It has been shown that DNA damage in cells that carry BRCA1/2 germline mutation was often overestimated and can hardly account for tumorigenesis [62]. Therefore, it was proposed that other factors, such as local inflammation and/or viral infections may put stress on immune system, which is already compromised by the germline mutation and promote the tumor formation in specific tissues such as breast and ovaries [63].

Cancer immunosuiveillence has been studied extensively for the last decade leading to the successful immune therapies [64]. Even though the immunotherapies are aimed at tumor microenvironment, systemic immunity is required for the process of tumor rejection after the therapy [65]. Thus, better understanding of the effects of cancer promoting hereditary mutations on the function of the systemic immunity may be very important. It can help to develop better immunotherapeutic strategies and new approaches for preventing and/or delaying the hereditary cancer, giving hope to many affected families.

Conclusions

We have shown that TMB and other genomic instability measures are not associated with the magnitude of the immune response in hereditary BRCA1/2 related breast and ovarian cancers, as well as in sporadic TNBC and sporadic HGSOC. However, high HR deficiency score, characteristic to all ovarian cancers and most subtypes of TNBC may be associated with the sensitivity to platinum and potentially also PARP inhibition. Hereditary tumors express phenotypes as heterogenous as sporadic tumors with various degree of “BRCAness” and various characteristics of the immune microenvironment. The subtyping criteria developed for sporadic tumors can be applied for the classification of hereditary tumors and possibly also characterization of their immune microenvironment. Further studies are needed to clarify if immunoediting plays a particular role in protecting against accumulation of genetic damage in BRCA1/2 germline mutation carriers, and what is the impact of germline mutations on the systemic immunity.

Availability of data and materials

The datasets analysed during the current study are available in CBioPortal (http://www.cbioportal.org), FireBrowse data version 2016_01_28 (http://firebrowse.org/), Cancer Research Institute iAtlas (https://www.cri-iatlas.org/about/), and The Cancer Immunome Atlas (TCIA) https://tcia.at/home).

The raw expression data for the cohort used in this study in TPM are given in Additional file 11: Table S11.

The RNA-seq raw data are available at the Gene Expression Omnibus under accession number GSE141142.

Abbreviations

BL1:

Basal like 1 subtype of TNBC

BL2:

Basal-like 2 subtype of TNBC

BRCA1/2:

Breast cancer type 1/2 susceptibility protein

CLOVAR:

Classification of Ovarian Cancer

CNV:

Copy number variation

DDR:

DNA damage response

DIF:

Differentiated subtype of ovarian cancer

ECM:

Extracellular matrix

EGF:

Epidermal growth factor

ER + :

Estrogen-receptor-positive

ERK1/2:

Extracellular signal–regulated kinases

FFPE:

Formalin-Fixed Paraffin-Embedded

HER2 + :

Human epidermal growth factor receptor 2 positive

HGSOC:

High grade serous ovarian cancers

HR:

Homologous recombination

HRD:

Homologous recombination deficiency

IGFR1:

Insulin-like growth factor 1

IMR:

Immunoreactive subtype of ovarian cancer

INF-γ:

Interferon gamma

JNK:

The c-Jun N-terminal kinase

LAR:

Luminal androgen receptor subtype of TNBC

M:

Mesenchymal subtype of TNBC

MAPK:

Mitogen-activated protein kinase

MES:

Mesenchymal subtype of ovarian cancer

MET:

A single pass tyrosine kinase receptor or hepatocyte growth factor receptor (HGFR)

MSL:

Mesenchymal-stem like subtype of TNBC

NGF:

Nerve growth factor

NGS:

Next generation sequencing

PARP:

Poly (ADP-ribose) polymerase

pCR:

Pathological complete response

PD-1:

Programmed cell death protein 1

PD-L1:

Programmed death-ligand 1

PI3K:

Phosphoinositide 3-kinases

PI3kinase:

Phosphoinositide 3-kinase

PRO:

Proliferative subtype of ovarian cancer

SNV:

Single nucleotide variation

SSGSEA:

Single sample gene set enrichment analysis

TCGA:

The Cancer Genome Atlas

TCIA:

The Cancer immunome Atlas

TGF-β:

Transforming growth factor beta 1

Th1:

T-helper cell 1

Th17:

T-helper cell 17

Tils:

Tumor infiltrating lymphocytes

TMB:

Tumor mutation burden

TNBC:

Triple negative breast cancer

TNF:

Tumor necrosis factor

Wnt/β catenin:

The canonical Wnt pathway

References

  1. 1.

    Burnet M. Cancer: a biological approach. III. Viruses associated with neoplastic conditions. IV. Practical applications. Br Med J. 1957;1:841–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Thomas L. Cellular and Humoral aspects of the hypersensitive states. New York: Hoeber-Harper; 1959. p. 529–32.

    Google Scholar 

  3. 3.

    Mouw KW, Goldberg MS, Konstantinopoulos PA, D’Andrea AD. DNA damage and repair biomarkers of immunotherapy response. Cancer Discov. 2017;7:675–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Ohmoto A, Yachida S. Current status of poly (ADP-ribose) polymerase inhibitors and future directions. OTT. 2017;10:5195–208.

    Article  Google Scholar 

  5. 5.

    Turner N, Tutt A, Ashworth A. Hallmarks of “BRCAness” in sporadic cancers. Nat Rev Cancer. 2004;4:814–9.

    Article  CAS  PubMed  Google Scholar 

  6. 6.

    Wang Y, Ung MH, Cantor S, Cheng C. Computational investigation of homologous recombination DNA repair deficiency in sporadic breast Cancer. Sci Rep. 2017;7:15742.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Lord CJ, Ashworth A. BRCAness revisited. Nat Rev Cancer. 2016;16:110–20.

    Article  CAS  PubMed  Google Scholar 

  8. 8.

    Goodman AM, Kato S, Bazhenova L, Patel SP, Frampton GM, Miller V, et al. Tumor mutational burden as an independent predictor of response to immunotherapy in diverse cancers. Mol Cancer Ther. 2017;16:2598–608.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Sedic M, Kuperwasser C. BRCA1-hapoinsufficiency: unraveling the molecular and cellular basis for tissue-specific cancer. Cell Cycle. 2016;15:621–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Joosse SA. BRCA1 and BRCA2: a common pathway of genome protection but different breast cancer subtypes. Nat Rev Cancer. 2012;12:372.

    Article  CAS  PubMed  Google Scholar 

  11. 11.

    Hyman DM, Solit DB, Arcila ME, Cheng DT, Sabbatini P, Baselga J, et al. Precision medicine at memorial Sloan Kettering Cancer center: clinical next-generation sequencing enabling next-generation targeted therapy trials. Drug Discov Today. 2015;20:1422–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Lehmann BD, Bauer JA, Chen X, Sanders ME, Chakravarthy AB, Shyr Y, et al. Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. J Clin Invest. 2011;121:2750–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Chen X, Li J, Gray WH, Lehmann BD, Bauer JA, Shyr Y, et al. TNBCtype: a subtyping tool for triple-negative breast Cancer. Cancer Informat. 2012;11:147–56.

    Article  Google Scholar 

  14. 14.

    Verhaak RGW, Tamayo P, Yang J-Y, Hubbard D, Zhang H, Creighton CJ, et al. Prognostically relevant gene signatures of high-grade serous ovarian carcinoma. J Clin Invest. 2013;123:517–25.

    CAS  PubMed  Google Scholar 

  15. 15.

    Barbie DA, Tamayo P, Boehm JS, Kim SY, Moody SE, Dunn IF, et al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. 2009;462:108–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Lehmann BD, Jovanović B, Chen X, Estrada MV, Johnson KN, Shyr Y, et al. Refinement of triple-negative breast Cancer molecular subtypes: implications for Neoadjuvant chemotherapy selection. PLoS One. 2016;11:e0157368.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Rooney MS, Shukla SA, Wu CJ, Getz G, Hacohen N. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell. 2015;160:48–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Iglesia MD, Parker JS, Hoadley KA, Serody JS, Perou CM, Vincent BG. Genomic analysis of immune cell infiltrates across 11 tumor types. J Natl Cancer Inst. 2016;108:1.

    Article  CAS  Google Scholar 

  19. 19.

    Iglesia MD, Vincent BG, Parker JS, Hoadley KA, Carey LA, Perou CM, et al. Prognostic B-cell signatures using mRNA-seq in patients with subtype-specific breast and ovarian cancer. Clin Cancer Res. 2014;20:3818–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Hackl H, Charoentong P, Finotello F, Trajanoski Z. Computational genomics tools for dissecting tumour-immune cell interactions. Nat Rev Genet. 2016;17:441–58.

    Article  CAS  PubMed  Google Scholar 

  21. 21.

    Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, et al. Pan-cancer Immunogenomic analyses reveal genotype-Immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. 2017;18:248–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Bindea G, Mlecnik B, Tosolini M, Kirilovsky A, Waldner M, Obenauf AC, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity. 2013;39:782–95.

    Article  CAS  PubMed  Google Scholar 

  23. 23.

    Hubalek M, Czech T, Müller H. Biological subtypes of triple-negative breast Cancer. Breast Care (Basel). 2017;12:8–14.

    Article  Google Scholar 

  24. 24.

    Burstein MD, Tsimelzon A, Poage GM, Covington KR, Contreras A, Fuqua SAW, et al. Comprehensive genomic analysis identifies novel subtypes and targets of triple-negative breast cancer. Clin Cancer Res. 2015;21:1688–98.

    Article  CAS  PubMed  Google Scholar 

  25. 25.

    Murakami R, Matsumura N, Mandai M, Yoshihara K, Tanabe H, Nakai H, et al. Establishment of a novel Histopathological classification of high-grade serous ovarian carcinoma correlated with Prognostically distinct gene expression subtypes. Am J Pathol. 2016;186:1103–13.

    Article  CAS  PubMed  Google Scholar 

  26. 26.

    Tothill RW, Tinker AV, George J, Brown R, Fox SB, Lade S, et al. Novel molecular subtypes of serous and endometrioid ovarian cancer linked to clinical outcome. Clin Cancer Res. 2008;14:5198–208.

    Article  CAS  PubMed  Google Scholar 

  27. 27.

    Cancer Genome Atlas Research Network. Integrated genomic analyses of ovarian carcinoma. Nature. 2011;474:609–15.

    Article  CAS  Google Scholar 

  28. 28.

    Omarini C, Guaitoli G, Pipitone S, Moscetti L, Cortesi L, Cascinu S, et al. Neoadjuvant treatments in triple-negative breast cancer patients: where we are now and where we are going. Cancer Manag Res. 2018;10:91–103.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Nik-Zainal S, Alexandrov LB, Wedge DC, Van Loo P, Greenman CD, Raine K, et al. Mutational processes molding the genomes of 21 breast cancers. Cell. 2012;149:979–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Severson TM, Peeters J, Majewski I, Michaut M, Bosma A, Schouten PC, et al. BRCA1-like signature in triple negative breast cancer: molecular and clinical characterization reveals subgroups with therapeutic potential. Mol Oncol. 2015;9:1528–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Schumacher TN, Schreiber RD. Neoantigens in cancer immunotherapy. Science. 2015;348:69–74.

    Article  CAS  Google Scholar 

  32. 32.

    Daud AI, Wolchok JD, Robert C, Hwu W-J, Weber JS, Ribas A, et al. Programmed death-ligand 1 expression and response to the anti-programmed death 1 antibody Pembrolizumab in melanoma. J Clin Oncol. 2016;34:4102–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Spigel DR, Schrock AB, Fabrizio D, Frampton GM, Sun J, He J, et al. Total mutation burden (TMB) in lung cancer (LC) and relationship with response to PD-1/PD-L1 targeted therapies. J Clin Oncol. 2016;34:9017.

    Article  Google Scholar 

  34. 34.

    Le DT, Uram JN, Wang H, Bartlett BR, Kemberling H, Eyring AD, et al. PD-1 blockade in tumors with mismatch-repair deficiency. N Engl J Med. 2015;372:2509–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang T-H, et al. The Immune Landscape of Cancer. Immunity. 2018;48:812–830.e14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Knijnenburg TA, Wang L, Zimmermann MT, Chambwe N, Gao GF, Cherniack AD, et al. Genomic and Molecular Landscape of DNA Damage Repair Deficiency across The Cancer Genome Atlas. Cell Rep. 2018;23:239–254.e6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. 37.

    van Beers EH, van Welsem T, Wessels LFA, Li Y, Oldenburg RA, Devilee P, et al. Comparative genomic hybridization profiles in human BRCA1 and BRCA2 breast tumors highlight differential sets of genomic aberrations. Cancer Res. 2005;65:822–7.

    PubMed  Google Scholar 

  38. 38.

    Wang C, Zhang J, Wang Y, Ouyang T, Li J, Wang T, et al. Prevalence of BRCA1 mutations and responses to neoadjuvant chemotherapy among BRCA1 carriers and non-carriers with triple-negative breast cancer. Ann Oncol. 2015;26:523–8.

    Article  CAS  PubMed  Google Scholar 

  39. 39.

    Wunderle M, Gass P, Häberle L, Flesch VM, Rauh C, Bani MR, et al. BRCA mutations and their influence on pathological complete response and prognosis in a clinical cohort of neoadjuvantly treated breast cancer patients. Breast Cancer Res Treat. 2018;171:85–94.

    Article  CAS  PubMed  Google Scholar 

  40. 40.

    Nicolas E, Bertucci F, Sabatier R, Gonçalves A. Targeting BRCA Deficiency in Breast Cancer: What are the Clinical Evidences and the Next Perspectives? Cancers (Basel). 2018;10:1.

    Article  CAS  Google Scholar 

  41. 41.

    Masuda H, Baggerly KA, Wang Y, Zhang Y, Gonzalez-Angulo AM, Meric-Bernstam F, et al. Differential response to neoadjuvant chemotherapy among 7 triple-negative breast cancer molecular subtypes. Clin Cancer Res. 2013;19:5533–40.

    Article  CAS  PubMed  Google Scholar 

  42. 42.

    Yee D, DeMichele A, Isaacs C, Symmans F, Yau C, Albain K, et al. Abstract GS3–08: Pathological complete response predicts event-free and distant disease-free survival in the I-SPY2 TRIAL. General Session Abstracts. Am Assoc Cancer Res. 2018:GS3–08 [Cited 2019 May 22]. Available from: http://cancerres.aacrjournals.org/lookup/doi/10.1158/1538-7445.SABCS17-GS3-08.

  43. 43.

    Strickland KC, Howitt BE, Shukla SA, Rodig S, Ritterhouse LL, Liu JF, et al. Association and prognostic significance of BRCA1/2-mutation status with neoantigen load, number of tumor-infiltrating lymphocytes and expression of PD-1/PD-L1 in high grade serous ovarian cancer. Oncotarget. 2016;7:13587–98.

    Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Nolan E, Savas P, Policheni AN, Darcy PK, Vaillant F, Mintoff CP, et al. Combined immune checkpoint blockade as a therapeutic strategy for BRCA1-mutated breast cancer. Sci Transl Med. 2017;9:eaal4922.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Sobral-Leite M, Van de Vijver K, Michaut M, van der Linden R, Hooijer GKJ, Horlings HM, et al. Assessment of PD-L1 expression across breast cancer molecular subtypes, in relation to mutation rate, BRCA1-like status, tumor-infiltrating immune cells and survival. Oncoimmunology. 2018;7:e1509820.

    Article  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Wang X, Li M. Correlate tumor mutation burden with immune signatures in human cancers. BMC Immunol. 2019;20:4.

    Article  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Chan TA, Yarchoan M, Jaffee E, Swanton C, Quezada SA, Stenzinger A, et al. Development of tumor mutation burden as an immunotherapy biomarker: utility for the oncology clinic. Ann Oncol. 2019;30:44–56.

    Article  CAS  PubMed  Google Scholar 

  48. 48.

    Burgess M, Puhalla S. BRCA 1/2-mutation related and sporadic breast and ovarian cancers: more alike than different. Front Oncol. 2014;4:19.

    PubMed  PubMed Central  Google Scholar 

  49. 49.

    Joosse SA, van Beers EH, Tielen IHG, Horlings H, Peterse JL, Hoogerbrugge N, et al. Prediction of BRCA1-association in hereditary non-BRCA1/2 breast carcinomas with array-CGH. Breast Cancer Res Treat. 2009;116:479–89.

    Article  CAS  PubMed  Google Scholar 

  50. 50.

    Schouten PC, Grigoriadis A, Kuilman T, Mirza H, Watkins JA, Cooke SA, et al. Robust BRCA1-like classification of copy number profiles of samples repeated across different datasets and platforms. Mol Oncol. 2015;9:1274–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Telli ML, Timms KM, Reid J, Hennessy B, Mills GB, Jensen KC, et al. Homologous recombination deficiency (HRD) score predicts response to platinum-containing Neoadjuvant chemotherapy in patients with triple-negative breast Cancer. Clin Cancer Res. 2016;22:3764–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Dunn GP, Koebel CM, Schreiber RD. Interferons, immunity and cancer immunoediting. Nat Rev Immunol. 2006;6:836–48.

    Article  CAS  PubMed  Google Scholar 

  53. 53.

    Safonov A, Jiang T, Bianchini G, Győrffy B, Karn T, Hatzis C, et al. Immune gene expression is associated with genomic aberrations in breast Cancer. Cancer Res. 2017;77:3317–24.

    Article  CAS  PubMed  Google Scholar 

  54. 54.

    Senovilla L, Vitale I, Martins I, Kepp O, Galluzzi L, Zitvogel L, et al. An anticancer therapy-elicited immunosurveillance system that eliminates tetraploid cells. Oncoimmunology. 2013;2:e22409.

    Article  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Senovilla L, Galluzzi L, Marino G, Vitale I, Castedo M, Kroemer G. Immunosurveillance against cancer-associated hyperploidy. Oncotarget. 2012;3:1270–1.

    Article  PubMed  PubMed Central  Google Scholar 

  56. 56.

    Davoli T, Uno H, Wooten EC, Elledge SJ. Tumor aneuploidy correlates with markers of immune evasion and with reduced response to immunotherapy. Science. 2017;355:eaaf8399.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Pusztai L, Karn T, Safonov A, Abu-Khalaf MM, Bianchini G. New strategies in breast Cancer: immunotherapy. Clin Cancer Res. 2016;22:2105–10.

    Article  CAS  PubMed  Google Scholar 

  58. 58.

    Sherman MH, Kuraishy AI, Deshpande C, Hong JS, Cacalano NA, Gatti RA, et al. AID-induced genotoxic stress promotes B cell differentiation in the germinal center via ATM and LKB1 signaling. Mol Cell. 2010;39:873–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. 59.

    Björkman A, Qvist P, Du L, Bartish M, Zaravinos A, Georgiou K, et al. Aberrant recombination and repair during immunoglobulin class switching in BRCA1-deficient human B cells. Proc Natl Acad Sci U S A. 2015;112:2157–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. 60.

    Friedenson B. The BRCA1/2 pathway prevents hematologic cancers in addition to breast and ovarian cancers. BMC Cancer. 2007;7:152.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. 61.

    Kim HS, Lee SW, Choi YJ, Shin SW, Kim YH, Cho MS, et al. Novel Germline mutation of BRCA1 gene in a 56-year-old woman with breast Cancer, ovarian Cancer, and diffuse large B-cell lymphoma. Cancer Res Treat. 2015;47:534–8.

    Article  CAS  PubMed  Google Scholar 

  62. 62.

    Levin B, Lech D, Friedenson B. Evidence that BRCA1- or BRCA2-associated cancers are not inevitable. Mol Med. 2012;18:1327–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. 63.

    Friedenson B. A theory that explains the tissue specificity of BRCA1/2 related and other hereditary cancers. J Med Med Sci. 2010;1:372–84.

    Google Scholar 

  64. 64.

    Finn OJ. A Believer’s overview of Cancer Immunosurveillance and immunotherapy. J Immunol. 2018;200:385–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Spitzer MH, Carmi Y, Reticker-Flynn NE, Kwek SS, Madhireddy D, Martins MM, et al. Systemic Immunity Is Required for Effective Cancer Immunotherapy. Cell. 2017;168:487–502.e15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors acknowledge the University of Connecticut Computational Biology Core (CBC) (https://bioinformatics.uconn.edu/) and the high performance computing facilities (https://health.uconn.edu/high-performance-computing) made available for conducting the research reported in this paper.

Funding

Generating RNA sequencing data and conducting this study was supported by a grant from the National Institutes of Health (NIH, R00LM011595, PI: Nabavi). The funding body had no influence in the design of the study, collection, analysis and interpretation of data and in writing the manuscript.

Author information

Affiliations

Authors

Contributions

EP designed the study, implemented the analysis, interpreted the results, and wrote the manuscript. TD implemented the analysis, and wrote the manuscript. AH processed the raw sequencing data. JE processed clinical data. UM conducted the clinical trial. GW designed the study. SN designed the study and wrote the manuscript. All authors read and approved the final version of the manuscript.

Corresponding author

Correspondence to Sheida Nabavi.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Board of Dana-Farber Harvard Cancer Center in Boston, Massachusetts, USA, protocol legacy 12–159. Clinical trials registration is with Clinicaltrials.gov (https://clinicaltrials.gov), with registration number NCT01623349. All procedures involving human participants were carried out in accordance with the Declaration of Helsinki. Written, informed consent was obtained from patients or guardians before enrolment in the study.

Consent for publication

Not applicable.

Competing interests

GMW reports grants from SU2C-AACR-DT0209, grants from Mary Kay Ash Foundation, grants from Ovarian Cancer Research Foundation, grants from Breast Cancer Alliance, grants from Breast Cancer Research Foundation, grants from NIH RO1 1R01CA226776–01, grants from Merck&Co, during the conduct of the study; In addition, Dr. Wulf has a patent Application 14/348810, Compositions and Methods for the Treatment of proliferative diseases pending, and a patent US 20090258352 A1, Pin1 as a marker for abnormal cell growth licensed to Cell Signaling; R&D Systems.

UAM reports personal fees from Astrazeneca, personal fees from Myriad Genetics, personal fees from Clovis, personal fees from Merck, personal fees from Eli Lilly, personal fees from Mersana, personal fees from Geneos, personal fees from Fuji Film, from 2X Oncology, personal fees from Cerulean, personal fees from Immunogen, personal fees from Novartis, outside the submitted work .

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Lines of treatment given to the patients, who were enroled into the trial (# NCT01623349) and whose expression profiles at diagnosis were analyzed in this study.

Additional file 2: Table S2.

Panther analysis of genes overexpressed in non-carriers. All carriers vs. non-car (PPEE 0.05 and 1.5 LogFC).

Additional file 3: Table S3.

Panther analysis of genes overexpressed in non-carriers. BRCA1 carriers versus non-carriers.

Additional file 4: Table S4.

Panther analysis of genes overexpressed in non-carriers. BRCA2 carriers versus non-carriers.

Additional file 5: Table S5.

Panther analysis of genes overexpressed in brest non-carriers. Breast carriers versus breast non-carriers.

Additional file 6: Table S6

. Panther analysis of genes overexpressed in ovarian non-carriers. Ovarian carriers versus ovarian non-carriers.

Additional file 7: Table S7.

Panther analysis of genes commonly overexpressed in breast and ovarian non-carriers.

Additional file 8: Table S8.

 Genes used for the analysis of the global immune activity of tumors.

Additional file 9: Table S9.

 Subtypes of hereditary and sporadic HGSOC from TCGA database.

Additional file 10: Table S10.

 HR deficiencies of platinum sensitive and resistant HGSOC from TCGA database.

Additional file 11: Table S11.

 The raw expression data for the cohort used in this study in TPM.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Przybytkowski, E., Davis, T., Hosny, A. et al. An immune-centric exploration of BRCA1 and BRCA2 germline mutation related breast and ovarian cancers. BMC Cancer 20, 197 (2020). https://doi.org/10.1186/s12885-020-6605-1

Download citation

Keywords

  • BRCA1
  • BRCA2
  • Breast cancer
  • Ovarian cancer
  • Tumor mutation burden
  • Homologous recombination deficiency
  • Immunotherapy
  • Biomarkers
  • BRCAness
  • Platinum resistance
  • PARP