Abstract
Purpose: The prevalence of mutations in cancer susceptibility genes such as BRCA1 and BRCA2 and other cancer susceptibility genes and their clinical relevance are largely unknown among a large series of unselected breast cancer patients in the Chinese population.
Experimental Design: A total of 8,085 consecutive unselected Chinese breast cancer patients were enrolled. Germline mutations in 46 cancer susceptibility genes were detected using a 62-gene panel.
Results: Pathogenic mutations were identified in 9.2% of patients among the 8,085 unselected breast cancer patients. Of these, 5.3% of patients carried a BRCA1 or BRCA2 mutation (1.8% in BRCA1 and 3.5% in BRCA2), 2.9% carried other breast cancer susceptibility genes (BOCG) and 1.0% carried another cancer susceptibility genes. Triple-negative breast cancers had the highest prevalence of BRCA1/2 mutations (11.2%) and other BOCG mutations (3.8%) among the four molecular subgroups, whereas ER−/PR−HER2+ breast cancers had the lowest mutations in BRCA1/2 (1.8%) and BOCG (1.6%). In addition, BRCA1 mutation carriers had a significant worse disease-free survival [unadjusted hazard ratio (HR) 1.60; 95% confidence interval (CI) 1.10–2.34; P = 0.014] and disease-specific survival (unadjusted HR 1.96; 95% CI, 1.03–3.65; P = 0.040) than did non-carriers, whereas no significant difference in survival was found between BRCA2 mutation carriers and non-carriers.
Conclusions: 9.2% of breast cancer patients carry a pathogenic mutation in cancer susceptibility genes in this large unselected series. Triple-negative breast cancers have the highest prevalence of mutations in BRCA1 /2 and other breast cancer susceptibility genes among the four molecular subgroups, whereas ER−/PR−HER2+ breast cancers had the lowest mutations in these genes. Clin Cancer Res; 23(20); 6113–9. ©2017 AACR.
Translational Relevance
The prevalence of mutations in cancer susceptibility genes and clinical relevance are not well investigated in a large unselected breast cancer series. In this study, germline mutations in 46 cancer susceptibility genes were detected in 8,085 consecutive unselected Chinese breast cancer patients. We found that 9.2% of breast cancer patients carried a pathogenic mutation in this large unselected series, and of these, 5.3% of patients carried BRCA1/2 mutations, 2.9% carried other breast cancer susceptibility genes, and 1.0% carried another cancer susceptibility genes. Triple-negative breast cancers had the highest prevalence of mutations in BRCA1/2 genes (11.2%) and other breast cancer susceptibility genes (3.8%) among the four molecular subgroups, whereas ER−/PR−HER2+ breast cancers had the lowest mutations in BRCA1/2 (1.8%) and other breast cancer susceptibility genes (1.6%). Our finding may be useful to select the right breast cancer to receive multiple-gene panel testing.
Introduction
Next-generation sequencing provides a unique platform to rapidly screen germline mutations in multiple genes in a huge number of patients. Recently, germline mutations in BRCA1/2 and other breast/ovarian cancer susceptibility genes detected by using multiple-gene panels assay have been reported in breast cancer patients, such as triple-negative breast cancer (TNBC) patients (1), early-onset or familial breast cancer patients (2–7), and those who are referred for BRCA1/2 testing (8–10), and more recently in unselected breast cancer patients with a relatively small sample size (n = 488; ref. 11).
The prevalence of mutations in BRCA1 and BRCA2 or other cancer susceptibility genes and their phenotype and clinical relevance are not well investigated in a large unselected breast cancer series. It is largely unknown whether BRCA1 and BRCA2 mutation carriers are associated with poor survival. Some previous studies suggested that BRCA1 mutation is associated with worse survival (12–18), others do not (19–24). More importantly, no studies so far are available for investigating the association between mutations in other breast cancer susceptibility genes and survival in large unselected breast cancer series. Another interesting issue is to explore whether molecular breast cancer subtypes based on estrogen receptor, progesterone receptor, and HER2 status could predict the prevalence of mutations in BRCA1 and BRCA2 and other breast cancer susceptibility genes.
In this study, we detected germline mutations in the 46 cancer susceptibility genes (including BRCA1 and BRCA2) using a 62-gene panel in a consecutive unselected large cohort of 8,085 breast cancer patients. The 46 cancer susceptibility genes were selected based on their roles in the development of breast cancer and other cancers.
We investigated the frequency of germline mutation in BRCA1 and BRCA2 genes and other breast cancer susceptibility genes and another cancer susceptibility genes in the entire cohort and in subgroups according to breast cancer family history, age at diagnosis, or molecular subtype and further explored the associations between the mutation carriers and clinicalpathological characteristics. Finally, we investigated whether mutations in BRCA1 and BRCA2 and other cancer susceptibility genes are associated with poorer survival compared with non-carriers in terms of disease-free survival (DFS) and disease-specific survival (DSS) in the entire cohort of 8,085 patients.
Materials and Methods
Study population
A total of 10,378 consecutive breast cancer patients were treated at the Breast Center of Peking University Cancer Hospital from October 2003 to May 2015. Among these, 8,085 cases have an enough genomic DNA quantity and quality drawn from blood samples to meet multi-gene panel analyses. The cohort of 8,085 breast cancer patients (91.7% were primary operable breast cancer with stage I to III disease) were unselected for age at diagnosis and family history of breast cancer. Family history of breast cancer is defined as the breast cancer patient had one or more breast and/or ovarian cancer patients in the first-, second-, or third-degree relatives; family history of any cancer is defined as the breast cancer patient had one or more cancer patients (any kind of cancers) in the first-, second-, or third-degree relatives.
Familial breast cancer is defined as breast cancer patients who had family history of breast and/or ovarian cancer; early-onset breast cancer (EBC) is defined as breast cancer patients who did not have family history of breast and/or ovarian cancer and were diagnosed at or before the age of 40 years; sporadic breast cancer (SBC) is defined as breast cancer patients who did not have a family history of breast and/or ovarian cancer and were diagnosed over the age of 40 years.
Tumor size was defined as the maximum tumor diameter measured by ultrasound at the time of diagnosis. The tumors were graded according to the modified Bloom–Richardson system. Estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2) status were determined using the breast tumor tissue obtained from a core-needle biopsy or taken from surgery. ER or PR immunostaining was considered positive when >1% of the tumor cells showed positive nuclear staining. HER2 positivity was defined as a score of 3+ via immunohistochemical staining or HER2 gene amplification via fluorescence in situ hybridization as described elsewhere (25). Informed written consent was obtained from all participants. The study was reviewed and approved by the Ethics Committee of Peking University Cancer Hospital (No. 2011KT12) and was performed in accordance with the Declaration of Helsinki.
Panel-based sequencing assay
Genomic DNA extracted from peripheral blood was captured target sequences using the SeqCap EZ hybridization and purification kit (Roche). Target sequencing was designed to cover all coding regions and intron–exon boundaries of the 62 cancer-gene panel (including 46 cancer susceptibility genes; Supplementary Table S1). Products were sequenced on a HiSeq 2500 (Illumina) to at least an average depth of 200-fold coverage, with a minimum of 30-fold required at every targeted position. Reads were aligned to the reference human genome GRCh37. Germline variations were called with GATK, Atlas2, and PlatyPus. Annotations were defined using ANNOVAR (http://annovar.openbioinformatics.org/en/latest/). Population allele frequencies were extracted from 1000 Genomes (http://www.1000genomes.org), dbSNP (version 138; http://www.ncbi.nlm.nih.gov/projects/SNP), and the Exome Variant Server (http://evs.gs.washington.edu/EVS).
Variant classification
In this study, we analyzed the genetic variants in 46 cancer susceptibility genes from the 62 cancer-gene panel. Only novel variants or variants with <1% population frequency in the 1000 Genomes were collected. The variant data were filtered to identify all protein-truncating mutations, including nonsense and frameshift mutations. For splice-site variants that have been verified by functional analyses from the literature were included. For synonymous, nonsynonymous, in-frame, and stop-loss variants, only variants classified as pathogenic or probably pathogenic by ClinVar (http://www.ncbi.nlm.nih.gov/clinvar/) were included in the analysis. The variant data were further assessed and classified according to American College of Medical Genetics and Genomics recommendations (26), with supporting data from function prediction software, public literature, and curated databases. Variants classified to be pathogenic or likely pathogenic were considered as pathogenic in this study. All the pathogenic mutations were validated by Sanger sequencing. If a patient carried more than one pathogenic mutation, mutation status was classified using the most damaging one when analyzed.
Statistical analyses
Categorical variables were compared between mutation carriers and non-carriers using the Chi-square test or the Fisher exact test, where appropriate. Continuous variables were tested with a t test. DFS was defined as from the time of diagnosis to first recurrence (local or distant), or contralateral breast cancer, or death from breast cancer (for patients without a recorded relapse) or date of last follow-up. DSS was defined as from the time of diagnosis to death from breast cancer or date of last follow-up. Survival was estimated using the Kaplan–Meier product-limit method and differences were tested for statistical significance using the log-rank test. A Cox proportional hazards model was used to determine whether a factor was associated with survival. Two-sided P values less than 0.05 were considered to be statistically significant. All analyses were performed using SPSS 20.0 software.
Results
Prevalence of germline mutations in cancer susceptibility genes
Germline mutations in 46 cancer susceptibility genes were detected in a large cohort of 8,085 unselected breast cancer patients using the 62-gene panel. The clinicopathological characteristics of these patients are presented in Supplementary Table S2. A total of 770 pathogenic mutations were identified in 743 (9.2%) of the 8,085 unselected breast cancer patients (Fig. 1; Supplementary Table S3). Twenty-seven patients (0.3%) had more than one mutation (Supplementary Table S4). Biallelic mutations in ATM and RAD51D were found in two patients (Supplementary Table S4).
Pathogenic mutations identified in 46 cancer susceptibility genes in 8,085 unselected breast cancer patients.
Four hundred and twenty-eight (5.3%) of the patients carried a germline mutation in BRCA1 or BRCA2 gene. Of these, 146 (1.8%) were in BRCA1 and 285 were (3.5%) in BRCA2; three women carried deleterious mutations in both BRCA1 and BRCA2 genes. In addition, 237 (2.9%) patients carried a total of 256 pathogenic mutations in other breast and/or ovarian cancer susceptibility genes (BOCG) beyond BRCA1 and BRCA2 genes, with PALB2 (n = 56), TP53 (n = 38), ATM (n = 31), RAD51D (n = 31), RECQL (n = 30), and CHEK2 (n = 27) ranked the most frequently mutated genes (Fig. 1).
Eighty-three (1.0%) patients carried pathogenic mutations in another cancer susceptibility genes in this cohort of 8,085 unselected breast cancer patients. Of these, 58 mutations were identified in gastrointestinal cancer susceptibility genes, in which 42 were in DNA mismatch repair (MMR) genes, 8 in POLE, 5 in EPCAM, 2 in KIT, and 1 in CTNNA1. The remaining 25 mutations were found in cancer susceptibility genes related to renal cancer (n = 15), hematologic tumors (n = 6), melanoma (n = 3), and thyroid tumors (n = 1; Supplementary Table S5). No mutations were found in AKT1, PIK3CA, APC, CDH1, SMAD4, BMPR1A, AXIN2, MEN1, FH, CDK4, and EGFR genes.
Germline mutations in cancer susceptibility genes in subgroups
BRCA1/2 mutation rates were 18.1% in familial breast cancer patients (FBC), 6.4% in EBC patients, and 3.3% in SBC patients, respectively (Table 1). The prevalence of BRCA1/2 mutations in the FBC group was significantly higher than that of EBC and SBC subgroups. Overall, the BRCA2 mutation rate was more pronounced than the BRCA1 mutation rate (3.5% versus 1.8%) and also in each subgroup (Table 1). The mutation frequencies of other BOCG and another cancer susceptibility genes were 4.6% and 1.4% in FBC, 3.6% and 0.7% in EBC, and 2.5% and 1.0% in SBC, respectively (Table 1). Furthermore, the overall prevalence of pathogenic mutations in FBC, EBC, and SBC were 24.1%, 10.7%, and 6.8%, respectively (Table 1).
Germline mutations in subgroups according to age at diagnosis and family history of breast cancer in this cohort of 8,085 unselected breast cancer patients
Next, we analyzed the pathogenic mutations in molecular subgroups based on ER, PR, and HER2 status (Table 2). The BRCA1 mutation rate was significantly higher in ER−/PR−, HER2− (TNBC) than that of the other subgroups, representing 7.4% in the TNBC group (Table 2). BRCA2 mutation rate was higher in ER+ and/or PR+, HER2− group than in other subgroups. Both BRCA1 and BRCA2 mutation rates were lower in ER−/PR−, and HER2+ group than in other groups (Table 2). Furthermore, mutations in other BOCG were more frequent in the TNBC group (3.8%), and were less in the ER−/PR−/HER2+ group (1.6%; Table 2). Mutations in another cancer susceptibility genes were not significantly different among the molecular subgroups (Table 2). Taken together, the overall prevalence of pathogenic mutations was the highest in the TNBC group (16.0%), following by ER+ and/or PR+, HER2− group (9.2%), ER+ and/or PR+/HER2+ group (6.7%), and the ER−/PR−/HER2+ group (4.7%; Table 2).
Germline mutations in molecular subgroups based on hormone receptor and HER2 status
Tumor characteristics and mutation status
The mean age at diagnosis of breast cancer was significantly younger for BRCA1 carriers (44.8 years, P < 0.001), BRCA2 carriers (47.8 years, P < 0.001), and other BOCG carriers (48.4 years, P < 0.001) when compared with non-carriers (51.4 years; Table 3). Both BRCA1 and BRCA2 mutation carriers and other BOCG mutation carriers were more likely to have a family history of cancer, especially a family history of breast and/or ovarian cancer (Table 3). Bilateral breast cancers were more frequent in BRCA1 and BRCA2 carriers compared with non-carriers (Table 3). BRCA1 mutation carriers were significantly associated with high-grade tumor (BRCA1 carriers versus non-carriers, 39.7% versus 15.3%, P < 0.001), and a triple-negative phenotype (58.2% versus 13.9%, P < 0.001; Table 3).
Clinicopathological characteristics between mutation carriers and non-carriers in this cohort of 8,085 unselected breast cancer patients
Survival
In this cohort of 8,085 unselected breast cancer patients, 114 patients with stage IV disease when administrated for treatment and 56 patients lost to follow-up were excluded for survival analysis, the remaining 7,915 (97.9%) primary operable breast cancer patients were then analyzed for survival. The median follow-up time of these 7,915 patients was 50 months (range, 1–152 months).
BRCA1 mutation carriers exhibited a significantly worse DFS and DSS than did non-carriers [DFS, unadjusted hazard ratio (HR) 1.60, 95% confidence interval (CI), 1.10–2.34, P = 0.014; DSS, unadjusted HR 1.94, 95% CI, 1.03–3.65, P = 0.040, respectively; Fig. 2; Supplementary Table S6]. Although BRCA2 mutations showed a slightly poor DFS and DSS when compared with non-carriers, the differences did not reach significance (Fig. 2; Supplementary Table S6). Mutations in other BOCG or in another cancer susceptibility genes were not significantly associated with DFS and DSS. Furthermore, in multivariate analysis, BRCA1 mutations did not show a poor survival after adjustment for age, tumor size, lymph node, tumor grade, ER, PR, HER2 status, and treatment; it did not reach a significance (Supplementary Table S7).
Survival analyses by Kaplan–Meier according to mutation status in the 7,915 primary operable breast cancer patients. A, B, DFS and DSS by mutation status, respectively. Other BOCG, other breast, and/or ovarian cancer susceptibility genes; another mutation, mutations in another cancer susceptibility genes unrelated to breast cancer.
We then analyzed the associations between BRCA1, BRCA2, and other BOCG mutations and survival in molecular subgroups stratified by ER, PR, and HER2 status. However, no significant associations were found between mutation carriers and non-carriers in the subgroups (Supplementary Table S8).
Discussion
In this study, we found that 9.8% of patients carried at least one pathogenic mutation in the 46 cancer susceptibility genes in the 8,085 unselected breast cancer patients. To our knowledge, this is the most comprehensive study undertaken to investigate the germline mutations in multiple cancer susceptibility genes in a large consecutive unselected breast cancer series to date.
Our data showed that 5.3% of patients carried a BRCA1 or BRCA2 mutation in this large series. BRCA1 and BRCA2 mutation rates in unselected breast cancer patients detected by next-generation sequencing assay have been only reported in a few studies. The prevalence of BRCA1/2 mutations in our current study is consistent with a recent report using multi-gene panel in a small series of 488 unselected breast cancer (6.1% in the entire cohort and 5.1% in non-Ashkenazi; ref. 11) and with the report from The Cancer Genome Atlas in 507 breast cancer patients (5.5%; ref. 27). Mutations in BRCA2 were almost twice compared with mutations in BRCA1 in our current large unselected series (3.5% vs. 1.8%), similar reports are also found in other Asian populations (28–31).
In addition, 2.9% of patients had a mutation in other breast cancer susceptibility genes in our large unselected series. Previous studies have reported the mutation rates in other breast cancer susceptibility genes using multi-gene panel assay in TNBC patients (3.7%; ref. 1) and in patients who are recommended for BRCA1/2 testing based on current guidelines (mutation rate ranges from 2.9% to 3.9%; refs. 8, 9). Considering our patients are unselected, the prevalence of mutations in other breast cancer susceptibility genes in our series is remarkable. Our data suggested that the most frequently mutated genes in non-BRCA1/2 breast cancer susceptibility genes were PALB2, TP53, RAD51D, and ATM, which are in agreement with most recent studies using a multiple-gene panel assay (1, 3, 5–7, 11, 32). However, the frequency of CHEK2 mutation in our large unselected breast cancer series was lower when compared with the mutation rate in Caucasian women (9, 11, 33). The possible reason might be the frequency of CHEK2 mutations depending on ethnic populations; for example, the CHEK2 recurrent mutation, 1100delC, is frequently found in Caucasian women but not in Chinese women (33). In contrast, the frequency of RAD51D mutation in this study was higher than that of the reports in Caucasian women (7, 11); the p.K91fs mutation of the RAD51D gene that is very rare in Caucasian women was very common in our series, which may have a founder effect in Chinese women.
We also reported that 1.0% of patients carried mutations in another cancer susceptibility genes in this series, particularly in mismatch repair genes, where the mutations are closely linked to Lynch syndrome. Although the association between mutations in mismatch repair genes and breast cancer risk is not fully established, increasing evidence suggested that mutations in mismatch repair genes may not only confer to a high risk for colorectal cancer, but may also increase breast cancer risk (34–36).
Clinical factors such as early-onset age and family history of breast or ovarian cancer were significantly associated with mutations in BRCA1/2 genes in our series. We also found that patients with a family history of breast cancer had a higher mutation rate in other breast cancer susceptibility genes than did SBC patients (4.6% vs. 2.5%). More importantly, germline mutation overall is significantly associated with TNBC phenotype, representing most frequently mutation rate among the four molecular subgroups, particularly for BRCA1 mutations. These data indicated that TNBCs exhibit genome instability and a more aggressive phenotype. In contrast, HER2-positive breast cancers exhibit the most less germline mutations among the four molecular subgroups not only in BRCA1/2 genes, but also in other breast cancer susceptibility genes. Interestingly, HER2 is a stronger driver gene that could initially trigger tumorigenesis per se; this notion is well demonstrated in mouse models (37).
There is conflicting evidence (18, 19, 21, 23, 38) so far regarding survival of breast cancer patients who carry BRCA1 or BRCA2 mutations because previous studies are based on selected high-risk populations, small study populations, or patients with founder mutations. In our study, we are able to comprehensively analyze survival in patients who carried a mutation in BRCA1 or BRCA2 or other cancer susceptibility genes in a large consecutive unselected breast cancer series. We found that BRCA1 mutation carriers had an unfavorable survival compared with non-carriers; this is not case for BRCA2 mutation carriers and patients who carried mutations in other breast cancer susceptibility genes. Recent studies suggested that breast cancer patients with a mutation in BRCA1/2 or other DNA repair genes are more likely to respond to poly-(ADP) ribose polymerase inhibitors (39, 40). Therefore, mutation carriers may benefit from the targeted therapy with optimal chemotherapy regimens.
There are two limitations in our study. First, although the sample size is large, it is a hospital-based population; second, the time of follow-up is relatively short (median 50 months); therefore, long-term follow-up is needed to investigate the survival in mutation carriers.
In conclusion, we found that 9.2% of breast cancer patients carry a pathogenic mutation in cancer susceptibility genes in this large consecutive unselected series. TNBCs have the highest prevalence of mutations in BRCA1 /2 and other breast cancer susceptibility genes among the four molecular subgroups, whereas ER−/PR−HER2+ breast cancers had the lowest mutations in these genes. Our findings have potential clinical implications, because strategies for early detection, prevention intervention, adjuvant chemotherapy regimens and targeted therapy are different between mutation carriers and non-carriers. As multiple-gene panel testing becomes more routine and cost-effective, our findings may be useful for selecting the right patients to receive multi-gene panel testing.
Disclosure of Potential Conflicts of Interest
J. Zhang is RD Director at Berry Genomics Ltd. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: Y. Xie
Development of methodology: J. Sun, H. Meng, M. Lv, J. Bai, J. Zhang
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): J. Sun, H. Meng, L. Yao, M. Lv, J. Bai, T. Ouyang, J. Li, T. Wang, Z. Fan, T. Fan, B. Lin, Y. Xie
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): J. Sun, H. Meng, L. Yao, M. Lv, J. Bai, L. Wang, Y. Xie
Writing, review, and/or revision of the manuscript: J. Sun, H. Meng, Y. Xie
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): L. Wang
Study supervision: Y. Xie
Grant Support
This study was supported by the National Natural Science Foundation of China (81372832), the National Science and Technology Support Program (2014BAI09B08), and the 973 project (2013CB911004).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Footnotes
Note: Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).
- Received December 23, 2016.
- Revision received February 27, 2017.
- Accepted July 13, 2017.
- ©2017 American Association for Cancer Research.