Inheritance of deleterious mutations at both BRCA1 and BRCA2 in an international sample of 32,295 women

Background Most BRCA1 or BRCA2 mutation carriers have inherited a single (heterozygous) mutation. Transheterozygotes (TH) who have inherited deleterious mutations in both BRCA1 and BRCA2 are rare, and the consequences of transheterozygosity are poorly understood. Methods From 32,295 female BRCA1/2 mutation carriers, we identified 93 TH (0.3 %). “Cases” were defined as TH, and “controls” were single mutations at BRCA1 (SH1) or BRCA2 (SH2). Matched SH1 “controls” carried a BRCA1 mutation found in the TH “case”. Matched SH2 “controls” carried a BRCA2 mutation found in the TH “case”. After matching the TH carriers with SH1 or SH2, 91 TH were matched to 9316 SH1, and 89 TH were matched to 3370 SH2. Results The majority of TH (45.2 %) involved the three common Jewish mutations. TH were more likely than SH1 and SH2 women to have been ever diagnosed with breast cancer (BC; p = 0.002). TH were more likely to be diagnosed with ovarian cancer (OC) than SH2 (p = 0.017), but not SH1. Age at BC diagnosis was the same in TH vs. SH1 (p = 0.231), but was on average 4.5 years younger in TH than in SH2 (p < 0.001). BC in TH was more likely to be estrogen receptor (ER) positive (p = 0.010) or progesterone receptor (PR) positive (p = 0.013) than in SH1, but less likely to be ER positive (p < 0.001) or PR positive (p = 0.012) than SH2. Among 15 tumors from TH patients, there was no clear pattern of loss of heterozygosity (LOH) for BRCA1 or BRCA2 in either BC or OC. Conclusions Our observations suggest that clinical TH phenotypes resemble SH1. However, TH breast tumor marker characteristics are phenotypically intermediate to SH1 and SH2. Electronic supplementary material The online version of this article (doi:10.1186/s13058-016-0768-3) contains supplementary material, which is available to authorized users.


Background
Women who have inherited mutations in BRCA1 or BRCA2 are at greatly increased risk of developing breast cancer (BC) and ovarian cancer (OC) [25,38]. Identification of a mutation at these loci can lead to risk or mortality reduction if optimal surveillance, risk-reducing mastectomy (RRM), and risk-reducing salpingo-oophorectomy (RRSO) are applied [8,29]. In addition, treatment of cancers in mutation carriers has advanced with the development of PARP inhibitors, which take advantage of the loss of BRCA1/2 function in tumors [37]. BRCA1 and BRCA2 are tumor suppressor genes, and tumors from the majority of mutation carriers have loss of heterozygosity (LOH), with loss of the normal allele, so there is no functioning protein [6,7,13,31]. In early studies, including a small number of tumor samples obtained from large BC and OC families, it was suggested that greater than 85 % of BRCA1or BRCA2-associated cancers exhibited LOH, and all showed loss of the normal allele.
The vast majority of BRCA1 and BRCA2 mutation carriers are single heterozygotes for BRCA1 (SH1) or BRCA2 (SH2). Homozygosity of missense alleles at BRCA2 (FANCD1) leads to Fanconi Anemia and increased cancer susceptibility, notably hematological malignancies [15,22]. At least three Fanconi Anemia cases are attributable to BRCA2/FANCD1 homozygous mutations [22]. Observations of homozygosity or compound heterozygosity at BRCA1 are very rare. Domchek et al. [9] reported a female patient with short stature, microcephaly, developmental delay, significant toxicity from chemotherapy, and epithelial ovarian carcinoma diagnosed at age 28 years. This woman was a compound heterozygote at BRCA1, with mutations c.2457delC (p.Asp821Ilefs*25) and c.5207 T > C (p.Val1736Ala). Both of these mutations are likely to be deleterious variants in BRCA1-associated cancer. The only other reported case of biallelic BRCA1 mutations was in a woman with multiple congenital anomalies consistent with a Fanconi anemia-like disorder and breast cancer at age 23 [30].
Transheterozygosity (TH) is the state of heterozygosity at two different loci. Here, we define TH to be inheritance of deleterious mutations in both BRCA1 and BRCA2. Reports on several BRCA1/2 transheterozygotes (TH) have been reported in the literature, mainly without further details on tumor or patient phenotype. Ramus et al. [27] reported on one TH who had been diagnosed with both BC and OC, and was identified as having a mutation in BRCA1 c.68_69delAG (185/ 187delAG) and BRCA2 c.5946delT (6174delT). LOH in these tumors was not found. Additional reports identified TH for BRCA1 c.2389G > T and BRCA2 c.3068dupA [21], BRCA1 c.68_69delAG and a BRCA2 c.5946delT [36], and TH with BRCA1 c.68_69delAG and BRCA2 c.5946delT [11] in four cases. In addition, a number of reports of TH with LOH in cancer samples have been published. Randall et al. [28] reported one TH identified with a BRCA1 c.3770_3771delGA and BRCA2 c.5946delT, and being affected with both BC and OC. For the BC, only LOH at the BRCA1 locus was found (not at BRCA2), and the OC sustained LOH at both BRCA1 and BRCA2. Tesoriero et al. [35] reported a TH with BRCA1 c.3770_3771delGA and BRCA2 c.5946delT. The BC of this patient lost the wild-type BRCA2 allele. Bell et al. [1] reported on a TH with c.5266dupC BRCA1 and c.5946delT BRCA2 mutation having three independent BCs. They showed that LOH occurred in two BRCA2 and one BRCA1 tumor. A large clinic-based series of 1191 carriers from Israel [20] identified 16 TH females, 14 with the c.68_69delAG BRCA1 and c.5946delT BRCA2 mutations and two with the c.5266dupC BRCA1 and c.5946delT BRCA2 mutations. A study from Germany identified eight female TH from 8162 BC/OC families and compared the clinical characteristics of the TH to their SH relatives and to SH in the family-based study [14].
To characterize the nature of TH and clinical phenotypes of TH, we used the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) dataset of 32,295 female BRCA1/2 mutation carriers ascertained in highrisk clinics and population-based studies. From this dataset, we investigated the occurrence of TH, we compared the characteristics and features of BC and OC in TH and single BRCA1 or BRCA2 mutations, and we examined LOH in as many cancer samples as possible.

Study sample
Details of CIMBA participating centers and data collection have been reported previously [5]. All the included mutation carriers participated in clinical and research studies at the host institutions after providing informed consent under IRB-approved protocols. Fifty-five centers and multicenter consortia (Additional file 1: Table S1) submitted data that met the CIMBA inclusion criteria [5]. Only female carriers with pathogenic BRCA1/2 mutations, concerning TH, SH1, and SH2 mutation carriers, were included in the current analysis. Pathogenicity of mutation was defined as follows: 1) generating a premature termination codon (PTC), except variants generating a PTC in exon 27 after codon 3010 of BRCA2; 2) large in-frame deletions that span one or more exons; and 3) deletion of transcription regulatory regions (promoter and/or first exon) expected to cause lack of expression of mutant allele. We also included missense variants considered pathogenic by using multifactorial likelihood approaches [4,12]. Mutations that did not meet the above criteria but have been classified as pathogenic by Myriad Genetics, Inc. (Salt Lake City, UT, USA) were also included.
Mutations are described using the Human Genome Variation Society (HGVS) nomenclature (http://www.HGVS.org/ varnomen) where the nucleotide numbering is from the A of the ATG translation initiator codon. For deletions or insertions, the most 3′ position possible was arbitrarily assigned as the altered nucleotide. The description of mutations of all types is given at the genomic level (using cDNA reference sequences NM_007294.3/ BRCA1 and NM_000059/BRCA2). BIC nomenclature was also presented for common variants that are familiar to many researchers and clinicians by their BIC designation (http://research.nhgri.nih.gov/bic). For BIC nomenclature, cDNA sequences were used as reference sequence (Genbank: U14680/BRCA1 and NM_000059.1/BRCA2). The nucleotide numbering is from nucleotide 1 of the cDNA gene sequence and for deletions or insertions the most 3′ position possible was arbitrarily assigned as the altered nucleotide.
In order to compare the TH with SH1 and SH2 mutation carriers on phenotypes of interest, we created a matched case-control set, in which "cases" were defined as TH, and "controls" were SH1 and SH2 mutation carriers. Matched SH1 "controls" carried a BRCA1 mutation found in the TH "case". Matched SH2 "controls" carried a BRCA2 mutation found in the TH "case". SH1 and SH2 were not matched to TH for any other characteristics. Using this approach, we identified 91 TH and 9316 matched SH1 mutation carriers, and 89 TH and 3370 matched SH2 mutation carriers.

Loss of heterozygosity
From 10 TH individuals, tumor tissue was available from twelve tumors, and blood DNA from 10 TH. From one case, tumor tissue from both BC and OC was available, and from another case affected with bilateral BC, tumor samples were available from both breast tumors. Hematoxylin and eosin (H&E) slides from each tumor were examined by a specialist pathologist. Areas of >80 % tumor cells were marked for macro-dissection. DNA from two 10-micron unstained slides was extracted using the Qiagen QIAmp DNA FFPE Tissue Kit using the standard protocol but with 500 μl deparaffinization solution.
We performed micro-satellite analysis to objectively detect LOH as described previously [16]. We amplified patient tumor and blood DNA for two markers within BRCA1 (D17S855 and D17S1322) and four markers around BRCA2 (D13S290, D13S260, D13S1698, and D13S171). The heterozygosity for these markers ranged from 0.46 to 0.82 [17,26]. Primer sequences and distance from BRCA1 or BRCA2 are given in Additional file 1 (Table S2). After polymerase chain reaction (PCR) amplification, samples were size-separated on a 96 capillary DNA analyzer (Applied Biosystems 3730xl). Data were analyzed using Genemapper Software (Applied Biosystems). For micro-satellites that were heterozygous, the ratios of allele peak heights for each tumor sample were compared to the allele peak heights for the blood DNA sample using the following formula L = (at2 X an1)/(at1 X an2), where L = the ratio; a = the height of the peak; n1 and n2 = normal allele 1 and normal allele 2; t1 and t2 = tumor allele 1 and tumor allele 2. All ten cases were informative for at least one marker in BRCA1. Where cases were informative for both markers, the LOH data were consistent for the two nearby markers. All ten cases were also informative for at least one of the four markers in BRCA2. In two cases, the data were not consistent across all markers in the 1.74 MB region and the data for the marker closest to BRCA2 was used.
To complement the information obtained from microsatellite analysis, we also undertook DNA sequence analysis. For each individual, a small region (<200 bp) around each of their two mutations was PCR-amplified from both tumor and blood DNA. DNA from peripheral blood of a healthy control individual was also amplified for each fragment as a control for no mutation. We used 10 ng of DNA in the PCR reaction, using standard protocol and primer sequences (given in Additional file 1: Table S3). All three samples for each mutation were then treated with EXO-SAP-IT (Affymetrix) and Sanger sequenced using standard methods [32]. This sequencing was used to confirm the presence of each mutation in the blood DNA from the patient and not in the control sample. We also assessed the mutation status in the tumor to determine if LOH had occurred. Since we extracted areas of >80 % tumor cells, both alleles can be present even when LOH is present, due to contaminating normal tissue. Therefore, for each tumor we determined for each mutation if the two alleles were at an equal ratio compared to the germline sample or if there was a decrease in one of the two alleles.

Statistical Analysis
For comparison of TH and SH mutation carriers, contingency table analysis using a chi-square test was used for dichotomous variables, and a t test for continuous variables. Fisher's exact tests were used if sample sizes in any contingency table cell were less than five. Analyses were done in STATA, v. 13.1. Table 1 describes the 93 female TH from 84 families identified from the CIMBA database. Among the matched TH-SH1/SH2 sets, 25 had no cancer diagnosis. The average age of these women was 39 years and the average age at diagnosis of BC was 41 years. Only 16 women were less than age 41 and 9 women were over age 41 at the time of diagnosis (mean age 49.9, range 41.4-67.9). Table 2 shows that OC age for the matched BRCA1 TH cases was 51.1 years and SH1 controls was 50.9 years (p = 0.154). For the matched BRCA2 set the average OC age for the TH cases was 54.7 years and for SH2 controls was 56.8 years (p = 0.421) (Fig. 1).

Characteristics of TH versus SH1 and SH2 mutation carriers
The most common TH involved inheritance of two of the three common Jewish mutations: 5 (5.4 %) women inherited BRCA1 c.5266dupC and BRCA2 c.5946delT; 31 (33.3 %) women inherited BRCA1 c.68_69delAG and BRCA2 c.5946delT. Six (6.5 %) women carried one of the three common Jewish mutations and another mutation. The majority of the remaining TH were observed only once. The majority of the TH self-identified as non-Hispanic Caucasian or Jewish. Of the 6907 women who carried one of the Jewish founder mutations, 2732 (39.6 %) self-identified as Jewish, 947 (13.7 %) were unknown, and 3225 (46.7 %) reported an ethnicity other than Jewish. We observed two TH in Hispanics and six TH in Asians (four of which were Korean). Of the 93 TH, 51 were diagnosed with BC only, 4 with OC only, 13 with both BC and OC, and 25 with no cancer diagnosis.
The matched datasets included 91 TH and 9316 SH1 for the BRCA1 matched analysis, and 89 TH and 3370 SH2 for the BRCA2 matched analysis. Two BRCA1 mutations were observed among the TH in our dataset that were not observed among SH1 (c.1390delA and c.3196G > T), and four BRCA2 mutations were observed in the TH dataset that were not observed among the SH2 (c.8633-?_8754 + ?amp, c.739_740delAT, c.5380delG, and c.2269A > T). These six carriers were not included in the analysis (denoted by asterisk in Table 1). TH were more likely to be born more recently (i.e., since 1961) than SH2 mutation carriers but not when compared to SH1s ( Table 2). The TH group consisted of more individuals from Asian ancestry compared to the SH1 and SH2 groups, with an excess of women having a Jewish ancestry vs. the SH1 group. TH were more likely to have ever been diagnosed with BC than SH1 or SH2 individuals (68.1 % vs. 52.0 %; p = 0.002, and 67.4 % vs. 50.4 %; p = 0.002), and TH were more likely to be diagnosed with OC than SH2 women (16.9 % vs. 9.3 %; p = 0.017), which was not observed in TH vs. SH1 women, perhaps due to the lower incidence of OC in BRCA2 vs. BRCA1. Age at BC diagnosis was significantly different for TH vs. SH2 (40.5 years vs. 45.0 years; p < 0.001), but there was no difference between TH and SH1.
There were 64 TH cases with BC. Of these, 62 TH were matched to 4846 SH1s and 60 TH were matched to 1699 SH2 (Table 3). TH were more likely to have estrogen receptor (ER)-and progesterone receptor (PR)positive BC than SH1s (ER: 42.9 % vs. 24.0 %; p = 0.010; PR: 40.6 % vs. 20.0 %; p = 0.013). In contrast, the BCs of TH were less likely ER-and PR-positive than in SH2s (ER: 42.9 % vs. 76.5 %; p < 0.0001; PR: 40.6 % vs. 62.8 %; p = 0.012). The proportion of ER-and PR-positive BCs in TH was intermediate to that of SH1 and SH2. No difference was seen regarding the HER2 status between the BCs of TH and SH1s and SH2s, respectively, although the available numbers were small. No differences in other BC characteristics (morphology, grade, stage) were observed.
Only 17 TH were diagnosed with OC, and thus we had limited data on features of OC to make inferences regarding differences in TH compared with SH1 or SH2. No statistically significant differences were observed for OC traits between TH and SHs (Table 4). Surprisingly, four borderline tumors were reported in both the SH1 and SH2 groups.

Loss of heterozygosity
Due to the frequent LOH in SH individuals, we examined the hypothesis that either BRCA1 or BRCA2 would be lost in each of the TH individuals due to LOH, and that whichever gene was lost could have an impact on their tumor characteristics. Of the 68 TH individuals with cancer, LOH analysis of three tumors from two cases had previously been published by our group using the same methods as the newly identified cases [27]. In the context of the current study, 12 additional tumor samples from 10 patients were analyzed (Table 5). We first used micro-satellite markers and an objective ratio of peak heights to determine if there was loss of one of the alleles when an individual was heterozygous [3] (Additional file 1: Tables S4 and S5). LOH analysis with micro-satellite markers normally includes linkage or segregation data to determine if the normal allele is lost. Since we did not have samples from other family members, we performed Sanger sequencing at the position of the mutations in both germline and tumor samples to determine which allele was lost. One sample failed for the sequencing so it was not possible to determine whether the normal or mutated allele was lost. Some samples showed loss of the mutant allele, which would suggest random loss. Tumors that exhibited LOH by micro-satellite analysis but did not indicate a decrease of the normal allele by sequencing were not considered to exhibit classic LOH. Following both sets of analyses and including our previously published data, one breast tumor (case 8) and one OC (case 2) showed LOH for BRCA1, two breast tumors (cases 9 and 11) showed LOH of BRCA2, and the remaining tumors provided no evidence for LOH at either BRCA1 or BRCA2 (Table 5).    [20], we report that TH were more likely than both SH1 or SH2 to be diagnosed with BC, which was also observed in our series. In addition to prior reports, we observed that TH were more likely to be diagnosed with OC compared with SH2s, but not compared with with SH1s. TH breast tumors were more likely to be ER-/PR-positive than in SH1, but less likely than in SH2 patients, without other different tumor or disease characteristics. A number of TH had not been diagnosed with cancer by the time this analysis was completed. Twenty-five TH in our cohort had no BC or OC diagnosis at the time of counseling or genotyping. The average age of these TH individuals was 39.1 years (range 20-68). Of these, 16 (64 %) were less than 41 years old at the time of study, which is the average age of BC diagnosis, and 23 (92 %) were younger than the average age of OC diagnosis (54 years) in the CIMBA data. Of these 25 unaffected TH women, 7 (28 %) reported a RRSO compared to 2751 (22.6 %) who underwent RRSO among the total set of SH controls without BC or OC (12,154). Two (8.0 %) cancer-free TH underwent bilateral risk-reducing mastectomy compared to 1076 (8.9 %) SH. In addition, we had missing data for a number of relevant variables that could have impacted some inferences. For example, of the 62 breast cancers in the TH groups, only 21 (34 %) reported stage information.

Discussion
Although this is the largest series of TH women reported to date, the study is still limited in a number of ways. TH were more likely to be born more recently (i.e., since 1961) than SH2, but not SH1. Since there is evidence that birth cohort may have an important effect on cancer risk [18], the risk associations reported here may require additional evaluation in the future. The higher incidence of BC in the TH group versus both SH1 and SH2 groups, and of OC in the TH vs. the SH2 cohort could be explained by non-random inclusion of TH in the sample, leading to potential biases in associations, and this may limit generalizability of the dataset. Our analyses also do not account for potentially important confounders and the longitudinal nature of the data to follow cancer cases from time of testing to either cancer diagnosis or censoring after risk-reducing surgery. Furthermore, the great majority of missing data on cancer features avoids that certain questions may be appropriately addressed from this type of dataset. Additional future studies are required to completely evaluate these clinically important unresolved issues, and hopefully with the ongoing multinational collaboration within consortia like CIMBA this will be possible in time.  BRCA1 and BRCA2, which would be consistent with the tumors being driven by loss of either BRCA1 or BRCA2. We attempted to determine the frequency of loss of each gene in a subset of cases where tumor material was available. Previously published data suggest a high rate of LOH with loss of the normal allele in the majority of BRCA1 and BRCA2 cases with strong family history at approximately 80 and 70 years, respectively [24]. However, we did not find loss of either BRCA1 or BRCA2 in the majority of tumors. The low frequency of LOH was consistent with the results from a previously published case (case 12) where we did not find LOH of either gene in either the breast or ovarian tumor [27]. Three other papers on TH showed LOH with loss of the normal allele [1,28,35]. One potential reason for the low frequency of LOH in this study could be that seven of the breast tumor samples were areas of ductal carcinoma in situ (DCIS) with micro-invasion rather than a region of the invasive breast tumor. However, we identified two tumors with LOH in these types of samples so this explanation is unlikely to be the major cause of the low rate of LOH.
The observed ages at diagnosis of BC in TH, SH1, and SH2, and the distributions of tumor characteristics may also reveal the interactions of BRCA1 and BRCA2 mutations, which may have implications for modeling the cancer susceptibility in TH. The observed age distributions rule out a multiplicative model for the interactions of BRCA1 and BRCA2 mutations on BC risk. Given the well-established BC risks for BRCA1 and BRCA2 mutations, a multiplicative model would imply very high cancer risk at young ages. However, the present study suggests that ages at BC diagnosis in TH are not significantly different from those in BRCA1 mutation carriers. Therefore, a multiplicative model of cancer risk for BRCA1 and BRCA2 is inconsistent with the current observations. This observation, combined with the fact that the tumor characteristics are intermediate to SH1 and SH2, suggests that an additive model for the joint effects of BRCA1 and BRCA2 mutations is more plausible. These results could be used for modeling the cancer risks for TH carriers and could be incorporated into risk prediction models.
Micro-satellite analysis alone did show a decrease in one of the two alleles in more of the tumors (6 out of 12 BRCA1 and 5 out of 12 BRCA2); however, the sequencing data suggested that the mutant allele rather than the normal allele was lost in many of the tumors. Although the early publications in high-risk families showed very high rates of LOH, exclusively with loss of the normal allele, more recently there have been many publications showing larger numbers of cases with no LOH [19,23,24] and an increasing number of tumors with loss of the mutant allele [19,24]. The second hit in these tumors could be due to a somatic mutation of the normal chromosome or due to promoter methylation, rather than LOH. Unfortunately, the amount of material from these tumors was very limited, and it was not possible to preform additional experiments to investigate alternative mechanisms. Methylation of BRCA1 has been shown to be a mechanism of decreased BRCA1 expression in sporadic BC [2,34], although this is less frequent in BRCA1 carriers [10,33]. Why the mechanism of LOH with loss of the normal allele in TH might be different compared with SH is unclear. Tumor material was only available in a small proportion of the cases with cancer. Therefore, it is difficult to interpret the results of the tumor study more broadly. Despite the small numbers, we did not find evidence to support the hypothesis that the tumors would either have LOH of BRCA1 or BRCA2. The TH breast tumor characteristics, however, do appear to be intermediate in phenotype to SH1 and SH2, suggesting some cancers are being driven by inactivation of BRCA1 and some by inactivation of BRCA2. Additional studies that explore other causes of inactivation (e.g., methylation or somatic mutation) are warranted.

Conclusions
We report evidence that the BRCA1 mutation in TH may drive these clinical TH phenotypes based on elevated OC risk in TH vs. SH2 but not SH1, and earlier age of BC diagnosis in TH vs. SH2 but not SH1. Therefore, TH may be managed more like BRCA1 mutation carriers than BRCA2 mutation carriers. In contrast, TH breast tumor characteristics (e.g., ER/PR status) are intermediate in phenotype to SH1 and SH2. Future studies are warranted to understand whether TH should be managed differently to SH1 or SH2 carriers, and, if so, to enable individualized counseling and clinical management appropriate for TH mutation carriers.

Additional files
Additional file 1: Table S1. Ethics committees that granted approval for the access and use of the data for this study. Table S2. Participant counts by center and mutation. Table S3. Primers used for PCR and Sanger sequencing. Table S4. Primers used in micro-satellite analysis for loss of heterozygosity. We wish to thank Yuan Chun Ding and Linda Steele for their work in participant enrollment and biospecimen and data management.

CBCS
We thank Bent Ejlertsen Ejlertsen and Anne-Marie Gerdes for the recruitment and genetic counseling of participants CNIO This work was partially supported by Spanish Association against Cancer (AECC08), RTICC 06/0020/1060, FISPI08/1120, Mutua Madrileña Foundation (FMMA) and SAF2010-20493 We thank Alicia Barroso, Rosario Alonso and Guillermo Pita for their assistance.