- 1Department of Medicine, MacKay Medical College, New Taipei City, Taiwan
- 2Department of General Surgery, MacKay Memorial Hospital, Taipei, Taiwan
- 3Department of Biomedical Sciences and Engineering, National Central University, Taoyuan, Taiwan
- 4Institute of Statistical Science, Academia Sinica, Taipei, Taiwan
- 5Department of Statistics, University of California Los Angeles, Los Angeles, CA, United States
Breast cancer is the most common invasive cancer in women worldwide. Next-generation sequencing (NGS) provides a high-resolution profile of cancer genome. Our study ultimately gives the insight for genetic screening to identify the minority of patients with breast cancer with a poor prognosis, who might benefit from the most intensive possible treatment. The detection of mutations can polish the traditional method to detect high-risk patients who experience poor prognosis, recurrence and death early. In total, 147 breast cancer tumors were sequenced with targeted sequencing using a RainDance Cancer Hotspot Panel. The average age of all 147 breast cancer patients in the study was 51.7 years, with a range of 21–77 years. The average sequencing depth was 5,222x (range 2,900x-8,633x), and the coverage was approximately 100%. A total of 235 variants in 43 genes were detected in 147 patients by high-depth Illumina sequencing. A total of 219 single nucleotide variations were found in 42 genes from 147 patients, and 16 indel mutations were found in 13 genes from 84 patients. After filtering with the 1000 Genomes database and for synonymous SNPs, we focused on 54 somatic functional point mutations. The functional point mutations contained 54 missense mutations in 22 genes. Additionally, mutation of genes within the RET, PTEN, CDH1, MAP2K4, NF1, ERBB2, RUNX1, PIK3CA, FGFR3, KIT, KDR, APC, SMO, NOTCH1, and FBXW7 in breast cancer patients were with poor prognosis. Moreover, TP53 and APC mutations were enriched in triple-negative breast cancer. APC mutations were associated with a poor prognosis in human breast cancer (log-rank P<0.001). Our study identified tumor mutation hotspot profiles in Taiwanese breast cancer patients, revealing new targetable gene mutations in Asian breast cancer patients.
Introduction
Breast cancer is the most common form and leading cause of cancer death for women worldwide, with nearly 45,000 deaths each year (1, 2), and it is also the highest incidence of all cancers in women in Taiwan (3). In the US, breast cancer is the most common cancer among women, and women have a 12% lifetime risk of developing breast cancer (4). Breast carcinoma is known to arise from a multifactorial process and is approximately twice as common in first-degree relatives of women with the disease as in the general population, consistent with genetic susceptibility to the disease. More frustrating from a cancer prevention viewpoint is that 40% of Taiwanese women who develop this disease fit none of the currently identified risk groups. This has prompted us to search for clues at the molecular level that may help in understanding breast tumor pathogenesis. Breast cancer is a heterogeneous disease, appearing in a variety of forms with different behaviors (5–7). Next-generation sequencing technology is based on massively parallel sequencing of millions of fragments using novel reversible terminator-based sequencing chemistry. This new technology provides high-throughput measurements of the whole genome sequence, single nucleotide polymorphisms (SNPs), copy number variations, insertions, and deletions and can be performed de novo, on the transcriptome or on targeted sequences.
In past few years, many genomes have been sequenced, and this has rapidly generated large volumes of genomic data (8, 9). Several large-scale genomic studies have revealed the heterogeneity of breast cancer (10–12), and these experiments provide rich genetic information for breast cancer research. Different targeted next-generation sequencing (NGS) panels have been used in breast cancer studies (13–15). The gene mutations have been different in different reported studies. Therefore, the genetic profiles might differ between geographic or ethnic populations (16). However, the genetic spectrum of Asian breast cancer remains limited, and the most common breast cancer is not well studied in the Taiwanese population. We used the RainDance Cancer Hotspot Panel for 147 Taiwanese breast cancer patients. The cancer panel contains 54 genes, including numerous known key driver genes (BRCA1, BRCA2 and ERBB2) in breast cancer. In our study, we focused on the relationship between somatic mutations and survival. The results will be useful for clinical prognosis and therapeutic target development. Cancer is a complex genetic disease driven by somatic mutations that accumulate in the genome during the lifetime of an individual. These somatic mutations could also additionally contribute to transcriptomic alterations, single nucleotide variants (SNVs), insertions and deletions (indels), copy number variations, and genome rearrangements (structural variants, SVs) (17). We hope our study ultimately paves the way for genetic screening to identify the minority of breast cancer patients with a poor prognosis, who might benefit from the most intensive possible treatment. The detection of mutations can improve our ability to detect high-risk patients who are likely to experience relapse and die early.
Materials and Methods
Clinical Patients
A total of 147 pathologically confirmed breast cancer samples were obtained from the tissue bank of MacKay Memorial Hospital. This investigation was performed after approval by the Institutional Review Board (MacKay Memorial Hospital: 12MMHIS121 and Academia Sinica: AS-IRB02-102086). The clinical information is provided in Supplementary Table 1.
Targeted Enrichment of Genomic DNA Regions for Next-Generation Sequencing
To identify rare and novel variants in genes associated with cancer, we used a microdroplet-based target enrichment (RainDance Technologies, Billerica, MA) hotspot cancer panel followed by high-throughput sequencing. A primer library for 54 susceptible genes related to cancer was designed using a proprietary pipeline. The primer library was combined dropwise with each sheared genomic DNA sample on a RainDance ThunderStorm™ Panel (RainDance). Microdroplet-based PCR of 54 genes targeting more than 13,000 known mutations found in the Catalogue of Somatic Mutations in Cancer (COSMIC) database was performed. All 40 samples enriched with targeted amplicons were sequenced on the MiSeq platform (Illumina, Inc., San Diego, CA), generating 250-bp paired-end reads according, to the manufacturer’s instructions. NGS raw sequencing data were submitted to the NCBI Sequence Read Archive with BioProject ID PRJNA758602.
Variant Calling
Reads were aligned against the human reference genome (UCSC Hg19) using BWA (7.5.a) (18). Variants called by GATK software (2.7-2) (19) and SAMtools software (0.1.19) (20) were annotated using the Annovar package (21).
Functional Enrchment Analysis of Genes With Somatic Mutations
We used Ingenuity Pathway Analysis (IPA) (QIAGEN company, CA), a web-based computational platform, to conduct functional enrichment analysis of genes. We input the set of genes with somatic mutations and used the core analysis enrichment tool with the default settings.
Sanger Sequencing Validation of APC Mutations
For APC mutations prioritized for Sanger validation, we designed PCR primers and validated the presence of candidate mutations by Sanger sequencing on samples using a 3500xL DX Genetic Analyzer (Life Technologies Co., Grand Island, New York).
Survival Analysis
We calculated the patients’ risk scores from the number of functional mutations and classified them into high-risk or low-risk groups with the medium risk score as the threshold. Kaplan-Meier survival curves were obtained and compared by log-rank tests.
Statistical Analysis
Samples with mutations and without mutations for each clinical phenotype were determined using 2 × 2 contingency tables, and Fisher’s exact test was used to calculate P-values in the detected mutated genes. All P-values were two sided, and statistical significance was defined as P < 0.05.
Results
The Cancer Hotspot Mutation Profile of Breast Cancer
In total, 147 breast cancer patients underwent targeted sequencing with the RainDance Cancer Hotspot Panel. The average age of all 147 breast cancer patients included in the study was 51.7 years, with a range of 21–77 years. Breast cancer samples were categorized as ER-positive, HER2-positive, PR-positive and triple-negative subtypes (Table 1). The average sequencing depth was 5,222x (range 2,900x-8,633x), and the coverage was approximately 100% (Table S2). A total of 235 variants in 43 genes were detected in 147 patients by high-depth Illumina sequencing (Figure 1). A total of 219 single nucleotide variations were found in 42 genes from 147 patients, and 16 indel mutations were found in 13 genes from 84 patients. After filtering with the 1000 Genomes database and for synonymous SNPs, we focused on 54 somatic functional point mutations. The functional point mutations contained 54 missense mutations in 22 genes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819555/fonc-12-819555-HTML/image_m/fonc-12-819555-g001.jpg)
Figure 1 Mutations and in the 147 breast cancer patients. The top panel shows a summary of the mutations in 43 cancer hotspot genes (see text and Table S2 for details). Patients are arranged from left to right by the number of mutations, shown in the top track. Colored rectangles indicate the mutation category observed in a given gene.
The Missense Somatic Variant Allele Frequency
Fifty-four missense somatic mutations were identified in 22 genes of 85 patients. The next-generation sequencing analysis provided a quantitative assessment of the mutation frequency in each sample. We found that the missense somatic mutations found in the tumor samples were heterozygous mutations present with high allelic frequencies, ranging from 17% to 82% of the alleles in tissue samples. Figure 2 shows that the most commonly mutated genes were PI3CA (5 mutations in 25 patients) and TP53 (17 mutations in 17 patients). The mean number of missense somatic mutations per tumor patient was 0.58 (range 0-4). The missense somatic variants found in the PIK3CA, TP53, MAP2K4 and APC genes were recurrent, with a frequency higher than 3% in the 147 breast cancer patients (Figure 2). These recurrent somatic mutations were validated by Sanger sequencing (Table S3). Finally, we found 37 functional mutations in 37 patients, and 7 patients had 2 functional mutations. Figure 3 showed that older breast cancer patients accumulated more mutation and MAP2K4 mutations occurred in patients aged 50 years and older. DNA mismatch repair-related gene mutations were found in MSH2 (2 mutations in 2 patients), MLH1 (2 mutations in 2 patients), and BRCA1 (1 mutation in 1 patient). The cancer hotspot genes were involved in 16 canonical pathways. The 22 genes with missense somatic mutations represented 10 canonical pathways (Table S4). PolyPhen database predicts the functional effect and estimates the score difference between variants. Based on the scores of PolyPhen, NOTCH1, ABL1, SMO, APC, KDR, FGFR3, MLH1, PIK3CA, MSH2, RUNX1, TP53, BRCA1, NF1 and CDH1 were the 14 most significantly mutated genes (PolyPhen2_HDIV score> 0.9, Table S4). In the present study, APC was mutated in 5 breast cancer specimens, including two frameshift mutations, namely, p.R414S, and p.V1125A with poor outcomes in our breast cancer cohort, all of which had not been reported previously and were deemed to be very strong evidence of pathogenicity in breast cancer. However, no significant association between other genes and cancer prognosis was found in our study. APC was a negative regulator of WNT signaling pathway. Mutations in APC causing loss of APC tumor-suppressive functions were the primary mechanism for hyperproliferation in colorectal cancer. APC is both a marker for the development of chemotherapy resistance and a potential therapeutic target (22). Figure 4 shows that most breast cancer patients did not have missense somatic mutations and that the missense somatic mutations detected in 85 patients affected a number of pathways enriched in genome integrity, apoptosis, MAPK signaling and PI3K signaling pathways. The IPA pathway analysis indicated that the 54 missense mutations were enriched in the hereditary breast cancer signaling pathway (P= 3.74E-13, Table S5), which is characterized by an inherited susceptibility to breast cancer due to mutations in high-penetrance susceptibility genes, including BRCA1, BRCA2, TP53 and PTEN.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819555/fonc-12-819555-HTML/image_m/fonc-12-819555-g002.jpg)
Figure 2 Somatic missense mutation profiles of breast cancers identified by a hotspot panel of 22 cancer-associated genes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819555/fonc-12-819555-HTML/image_m/fonc-12-819555-g003.jpg)
Figure 3 Genetic mutations identified by the RainDance™ gene panel in 147 breast cancers. The Oncoprint illustrated the distribution of somatic mutations according to age at diagnosis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819555/fonc-12-819555-HTML/image_m/fonc-12-819555-g004.jpg)
Figure 4 Ten interaction networks constructed from canonical maps including the mutations found in 147 luminal breast tumors. In the concentric circle diagram, tumors are arranged as radial spokes and categorized by their mutation status in each network (concentric ring color). Orange, genome integrity, proteolysis, and apoptosis; green, MAPK signaling and PI3K signaling; blue, RTK signaling and Hippo signaling; purple, Notch signaling, Hedgehog signaling, and Wnt signaling.
Correlating Mutations With Survival
We identified 101 somatic variants in breast cancer patients. In a Cox univariate analysis, only 54 nonsilent mutations were significantly associated with overall survival (P= 0.04, Table 2). All germline mutations were not significant in predicting overall survival in breast cancer patients (Table S6). Additionally, we identified that the 54 nonsilent somatic mutations were not significantly associated with overall survival prediction in breast cancer patients (log-rank P= 0.21, Figure S1). According to the functional pathways related to the identified genes, we distinguished two mutation signatures. One was a genomic instability signature including somatic mutations in RB1, TP53, BRCA1, BRCA2, MSH2, MLH1 and ABL1, and the other was a kinase signaling signature including somatic mutations in RET, PTEN, CDH1, MAP2K4, NF1, ERBB2, RUNX1, PIK3CA, FGFR3, KIT, KDR, APC, SMO, NOTCH1 and FBXW7. The kinase signaling signature significantly predicted overall survival in the 147 breast cancer patients (log-rank P=0.048 Figure 5), but the genomic instability signature did not significantly predict patient survival (log-rank P=0.51 Figure S2). Mutations in the kinase signaling pathway may have a negative prognostic role in survival in breast cancer. We compared the mutation frequency with TCGA database (Figure S3). We have provided the results in Figure S4. Kaplan-Meier survival analysis showed that there were significant in APC, FBXW7, FGFR3, KIT and RB1 groups. Despite there are some significant in those genes, only one or two patients have mutation in these genes except APC and RB1 gene.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819555/fonc-12-819555-HTML/image_m/fonc-12-819555-g005.jpg)
Figure 5 Kaplan–Meier survival curves for overall survival based on kinase signaling pathway in breast cancer patients.
Correlating Mutations with Triple-Negative Breast Cancer
Triple-negative breast cancer (TNBC) is defined as a breast tumor that is negative for estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor-2 (HER2). These tumors are usually marked by higher rates of metastases and recurrence and are insensitive to traditional hormone treatment and chemotherapy (23). We evaluated the relationship between the 22 genes with mutations and triple-negative breast cancer patients. Table 3 indicates that mutations in the APC (P =0.006) and TP53 (P =0.045) genes were enriched in triple-negative breast cancer patients. Other genes were not enriched in any hormone receptor positive (ER/PR+) and HER2+ cohorts (Table S7). TP53-mutated breast cancer patients had significantly higher activity than TP53-wild-type patients. Meanwhile, we identified APC genes enriched with ER-negative and TNBC subtypes with statistical significance. In this subset analysis, mutations in novel codons that directly led to amino acid changes were related to the poorest overall survival. We further conducted multivariate Cox proportional hazard regression analysis with a signature containing the 37 novel mutations and other prognostic factors as predictors. We also evaluated the relationship between single gene mutation and survival in triple-negative breast cancer patients. A univariate Cox proportional hazard model showed that APC had the most significant impact on overall survival (P = 0.003) and significantly predicted overall survival in the 147 breast cancer patients (log-rank P<0.001, Figure 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819555/fonc-12-819555-HTML/image_m/fonc-12-819555-t003.jpg)
Table 3 Evaluating associations between gene mutation status and triple-negative breast cancer patients by Fisher’s exact test.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819555/fonc-12-819555-HTML/image_m/fonc-12-819555-g006.jpg)
Figure 6 Kaplan–Meier survival curves for overall survival based on APC mutations in breast cancer patients.
Discussion
An average of 1.02 to 1.66 somatic mutations per Mb in coding regions have been found in breast cancers (2, 12, 24), which translates into a mean of 56.9 (range 5-374) somatic mutations per tumor (25). Furthermore, significantly mutation genes in breast cancers have been identified, and these genes significantly affect pathways of functional or cellular processes, including the TP53 or PI3K pathway (e.g., PIK3CA, PTEN, and AKT1), MAPK/JNK signaling (e.g., MAP3K1, MAP2K4, and NF1), transcription factors and regulators (e.g., GATA3, RUNX1, and CBFB), splicing factors (e.g., SF3B1), and chromatin remodelers (e.g., MLL3 and ARID1A). It has been reported that only TP53, PIK3CA and GATA33 are consistently mutated in more than 10% of unselected breast cancers, while the remaining genes are mutated in less than 7% of patients, with a very long list of genes mutated in less than 1% of patients (2). Therefore, we applied targeted next-generation sequencing panel (with a high-depth sequencing method) for the detection of novel functional point mutations. The whole-genome sequences of the breast tumors indicated that the 10 most frequently mutated genes were TP53, PIK3CA, MYC, CCND1, PTEN, ERBB2, the ZNF703/FGFR1 locus on chromosome 8, GATA3, RB1 and MAP3K1[11]. However, we found that the 10 most frequently mutated genes were PIK3CA (17%), TP53 (12%), MAP2K4 (5%), APC (3.4%), RUNX1 (2.7%), RB1 (2.7%), NF1 (1.4%), MSH2 (1.4%), MLH1 (1.4%), KIT (1.4%), and CDH1 (1.4%). Consistent with previous mutational studies, PIK3CA and TP53 in our study were the most frequently mutated genes (12, 13, 26, 27). In our Asian breast cancer patient data, MAP2K4 had a higher mutation frequency than that seen in Caucasian breast cancer patients. Recently, the Taiwanese breast cancer study showed the high mutation frequency in PIK3CA, and TP53 (28, 29), but they did not investigate the mutations associated with patient survivals.
We indicated that missense somatic mutations were the most associated with survival. Germline mutations have been shown to impact tumor evolution and inheritance mechanisms (30, 31). Somatic mutation data in the Catalogue of Somatic Mutations in Cancer (COSMIC) (32) and TCGA provide somatic mutational landscapes across all cancer types. Missense somatic mutations in protein-coding regions can generate neoantigens (33, 34), which, in turn, could be associated with patient survival (35). Higher somatic mutation and neoantigen burdens are correlated with poor progression-free survival in multiple myeloma (36). This is consistent with our findings.
Mutations of the APC tumor-suppressor gene are thought to regulate beta-catenin within the Wingless/Wnt signal transduction pathway (37). We found that the number of missense somatic mutations and mutations in APC (V1125A and R414S) were independent predictors for the overall survival in our breast cancer patients. Here, we first reported the identification of a novel APC mutation in Taiwanese breast cancers patients. Furthermore, we also revealed that mutations in the TP53 and APC genes were enriched in TNBC. The TP53 mutation frequency was higher in TNBC than in other types of breast cancer (38). TP53-mutated breast cancer patients had significantly higher activity than TP53-wild-type patients. Furthermore, we explored the molecular markers correlated with the differences in immune activities between TP53-mutated and TP53-wild-type patients (39). TP53 mutation regulates immune responses. This pathway is involved in host immunity, influencing both innate and adaptive immune responses (40). Wnt signaling dysfunction has been reported to mediate the progression of triple-negative breast cancer (41). A proteomics study showed that TNBC tumors with APC mutations expressed median levels of β-catenin that were 4-fold higher than those in tumors without APC mutations (42). β-Catenin was required for the tumor cell motility in triple-negative breast cancer cells (43). Wnt/β-catenin pathway activation has been implicated in stem cell self-renewal, maintenance, and differentiation (44, 45). We first reported that APC somatic mutations (V1125A and R414S) were enriched in TNBC and associated with reduced overall survival in all breast cancer patients. These mutations may identify new biomarkers for breast cancer patients in the future.
Conclusions
In conclusion, we have analyzed targeted sequencing data generated from 147 Taiwanese breast cancer patients by using RainDance Cancer Hotspot Panel. We observed 54 somatic functional point mutations distributed across 22 genes in 85 patients. Mutations of APC and TP53 were enriched in TPNBC and patients with Patients with APC mutations had significantly poorer overall survival than patients without APC mutations. However, the mutations in kinase genes had significantly poorer overall survival than patients without kinase gene mutations. The results indicated that APC variant V1125A and R414S might be the oncogenic and the kinase genes association of survival in breast cancer patients which provides novel target for the development of breast cancer and identifies APC as a novel putative target for targeted therapy.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
Written informed consent was not obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
P-SY, C-LL, K-CL and Y-CH made substantial contributions to conception and design. C-LL, P-SY, and Y-CC were responsible for clinical management and interpreted clinical data. P-SY, K-CL, Y-TC, and Y-CH wrote the manuscript. YCH, P-SY and C-FL interpreted the data and analyzed the data. All authors gave final approval for the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to thank the Pharmacogenomics Laboratory (PGL) of the National Research Program for Biopharmaceuticals (NRPB) for the support and the Cancer Database of MacKay Memorial Hospital in Taiwan for the clinical data in this study. This research was supported by Academia Sinica (Mathematics in Biology), the Ministry of Science and Technology (MOST 107-2314-B-008 -002 and MOST 110-2321-B-195-001) and MacKay Memorial Hospital (MMH-107-87) and Ministry of Health and Welfare and National Health Research Institutes, Taiwan (NHRI-109BCCO-MF-202019-02).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.819555/full#supplementary-material
References
1. Nandy A, Gangopadhyay S, Mukhopadhyay A. Individualizing Breast Cancer Treatment-The Dawn of Personalized Medicine. Exp Cell Res (2014) 320:1–11. doi: 10.1016/j.yexcr.2013.09.002
2. Cancer Genome Atlas N. Comprehensive Molecular Portraits of Human Breast Tumours. Nature (2012) 490:61–70. doi: 10.1038/nature11412
3. Kohler BA, Sherman RL, Howlader N, Jemal A, Ryerson AB, Henry KA, et al. Annual Report to the Nation on the Status of Cancer, 1975-2011, Featuring Incidence of Breast Cancer Subtypes by Race/Ethnicity, Poverty, and State. J Natl Cancer Inst (2015) 107:djv048. doi: 10.1093/jnci/djv048
4. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2019. CA Cancer J Clin (2019) 69:7–34. doi: 10.3322/caac.21551
5. Ng CK, Schultheis AM, Bidard FC, Weigelt B, Reis-Filho JS. Breast Cancer Genomics From Microarrays to Massively Parallel Sequencing: Paradigms and New Insights. J Natl Cancer Inst (2015) 107:djv015. doi: 10.1093/jnci/djv015
6. Reis-Filho JS, Pusztai L. Gene Expression Profiling in Breast Cancer: Classification, Prognostication, and Prediction. Lancet (2011) 378:1812–23. doi: 10.1016/S0140-6736(11)61539-0
7. Weigelt B, Reis-Filho JS. Histological and Molecular Types of Breast Cancer: Is There a Unifying Taxonomy? Nat Rev Clin Oncol (2009) 6:718–30. doi: 10.1038/nrclinonc.2009.166
8. Mardis ER, Ding L, Dooling DJ, Larson DE, McLellan MD, Chen K, et al. Recurring Mutations Found by Sequencing an Acute Myeloid Leukemia Genome. N Engl J Med (2009) 361:1058–66. doi: 10.1056/NEJMoa0903840
9. Wang J, Wang W, Li R, Li Y, Tian G, Goodman L, et al. The Diploid Genome Sequence of an Asian Individual. Nature (2008) 456:60–5. doi: 10.1038/nature07484
10. Angus L, Smid M, Wilting SM, van Riet J, Van Hoeck A, Nguyen L, et al. The Genomic Landscape of Metastatic Breast Cancer Highlights Changes in Mutation and Signature Frequencies. Nat Genet (2019) 51:1450–8. doi: 10.1038/s41588-019-0507-7
11. Nik-Zainal S, Davies H, Staaf J, Ramakrishna M, Glodzik D, Zou X, et al. Landscape of Somatic Mutations in 560 Breast Cancer Whole-Genome Sequences. Nature (2016) 534:47–54. doi: 10.1038/nature17676
12. Stephens PJ, Tarpey PS, Davies H, Van Loo P, Greenman C, Wedge DC, et al. The Landscape of Cancer Genes and Mutational Processes in Breast Cancer. Nature (2012) 486:400–4. doi: 10.1038/nature11017
13. Kota K, Wald AI, Nikiforov YE, Nikiforova MN, Oesterreich S, Lee AV. Targeted Mutation Detection in Breast Cancer Using MammaSeq. Breast Cancer Res (2019) 21:22. doi: 10.1186/s13058-019-1102-7
14. Liang X, Vacher S, Boulai A, Bernard V, Baulande S, Bohec M, et al. Targeted Next-Generation Sequencing Identifies Clinically Relevant Somatic Mutations in a Large Cohort of Inflammatory Breast Cancer. Breast Cancer Res (2018) 20:88. doi: 10.1186/s13058-018-1007-x
15. Di Cosimo S, Appierto V, Silvestri M, Pruneri G, Vingiani A, Perrone F, et al. Targeted-Gene Sequencing to Catch Triple Negative Breast Cancer Heterogeneity Before and After Neoadjuvant Chemotherapy. Cancers (Basel) (2019) 11:1753. doi: 10.3390/cancers11111753
16. Tan DS, Mok TS, Rebbeck TR. Cancer Genomics: Diversity and Disparity Across Ethnicity and Geography. J Clin Oncol (2016) 34:91–101. doi: 10.1200/JCO.2015.62.0096
17. Campbell PJ, Stephens PJ, Pleasance ED, O'Meara S, Li H, Santarius T, et al. Identification of Somatically Acquired Rearrangements in Cancer Using Genome-Wide Massively Parallel Paired-End Sequencing. Nat Genet (2008) 40:722–9. doi: 10.1038/ng.128
18. Liu H. (2013). Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv 1303.3997v2 [q-bio.GN]
19. DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A Framework for Variation Discovery and Genotyping Using Next-Generation DNA Sequencing Data. Nat Genet (2011) 43:491–8. doi: 10.1038/ng.806
20. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map Format and SAMtools. Bioinformatics (2009) 25:2078–9. doi: 10.1093/bioinformatics/btp352
21. Wang K, Li M, Hakonarson H. ANNOVAR: Functional Annotation of Genetic Variants From High-Throughput Sequencing Data. Nucleic Acids Res (2010) 38:e164. doi: 10.1093/nar/gkq603
22. Zhang L, Shay JW. Multiple Roles of APC and its Therapeutic Implications in Colorectal Cancer. J Natl Cancer Inst (2017) 109(8):djw332. doi: 10.1093/jnci/djw332
23. Dent R, Trudeau M, Pritchard KI, Hanna WM, Kahn HK, Sawka CA, et al. Triple-Negative Breast Cancer: Clinical Features and Patterns of Recurrence. Clin Cancer Res (2007) 13:4429–34. doi: 10.1158/1078-0432.CCR-06-3045
24. Banerji S, Cibulskis K, Rangel-Escareno C, Brown KK, Carter SL, Frederick AM, et al. Sequence Analysis of Mutations and Translocations Across Breast Cancer Subtypes. Nature (2012) 486:405–9. doi: 10.1038/nature11154
25. Kandoth C, McLellan MD, Vandin F, Ye K, Niu B, Lu C, et al. Mutational Landscape and Significance Across 12 Major Cancer Types. Nature (2013) 502:333–9. doi: 10.1038/nature12634
26. Razavi P, Chang MT, Xu G, Bandlamudi C, Ross DS, Vasan N, et al. The Genomic Landscape of Endocrine-Resistant Advanced Breast Cancers. Cancer Cell (2018) 34:427–438.e426. doi: 10.1016/j.ccell.2018.08.008
27. Mathioudaki A, Ljungstrom V, Melin M, Arendt ML, Nordin J, Karlsson Å, et al. Targeted Sequencing Reveals the Somatic Mutation Landscape in a Swedish Breast Cancer Cohort. Sci Rep (2020) 10:19304. doi: 10.1038/s41598-020-74580-1
28. Wu CH, Hsieh CS, Chang YC, Huang CC, Yeh HT, Hou MF, et al. Differential Whole-Genome Doubling and Homologous Recombination Deficiencies Across Breast Cancer Subtypes From the Taiwanese Population. Commun Biol (2021) 4:1052. doi: 10.1038/s42003-021-02597-x
29. Huang CC, Tsai YF, Liu CY, Chao TC, Lien PJ, Lin YS, et al. Comprehensive Molecular Profiling of Taiwanese Breast Cancers Revealed Potential Therapeutic Targets: Prevalence of Actionable Mutations Among 380 Targeted Sequencing Analyses. BMC Cancer (2021) 21:199. doi: 10.1186/s12885-021-07931-4
30. Liaw D, Marsh DJ, Li J, Dahia PL, Wang SI, Zheng Z, et al. Germline Mutations of the PTEN Gene in Cowden Disease, an Inherited Breast and Thyroid Cancer Syndrome. Nat Genet (1997) 16:64–7. doi: 10.1038/ng0597-64
31. Hu C, Hart SN, Polley EC, Gnanaolivu R, Shimelis H, Lee KY, et al. Association Between Inherited Germline Mutations in Cancer Predisposition Genes and Risk of Pancreatic Cancer. JAMA (2018) 319:2401–9. doi: 10.1001/jama.2018.6228
32. Forbes SA, Beare D, Boutselakis H, Bamford S, Bindal N, Tate J, et al. COSMIC: Somatic Cancer Genetics at High-Resolution. Nucleic Acids Res (2017) 45:D777–83. doi: 10.1093/nar/gkw1121
33. Gubin MM, Zhang X, Schuster H, Caron E, Ward JP, Noguchi T, et al. Checkpoint Blockade Cancer Immunotherapy Targets Tumour-Specific Mutant Antigens. Nature (2014) 515:577–81. doi: 10.1038/nature13988
34. Colli LM, Machiela MJ, Myers TA, Jessop L, Yu K, Chanock SJ. Burden of Nonsynonymous Mutations Among TCGA Cancers and Candidate Immune Checkpoint Inhibitor Responses. Cancer Res (2016) 76:3767–72. doi: 10.1158/0008-5472.CAN-16-0170
35. Hsu YC, Chang YH, Chang GC, Ho BC, Yuan SS, Li YC, et al. Tumor Mutation Burden and Recurrent Tumors in Hereditary Lung Cancer. Cancer Med (2019) 8:2179–87. doi: 10.1002/cam4.2120
36. Miller A, Asmann Y, Cattaneo L, Braggio E, Keats J, Auclair D, et al. High Somatic Mutation and Neoantigen Burden are Correlated With Decreased Progression-Free Survival in Multiple Myeloma. Blood Cancer J (2017) 7:e612. doi: 10.1038/bcj.2017.94
37. Fodde R, Smits R, Clevers H. APC, Signal Transduction and Genetic Instability in Colorectal Cancer. Nat Rev Cancer (2001) 1:55–67. doi: 10.1038/35094067
38. Won KA, Spruck C. Triple-negative breast cancer therapy: Current and future perspectives. Int J Oncol (2020) 57:1245–61. doi: 10.3892/ijo.2020.5135. Epub 2020 Oct 16.
39. Liu Z, Jiang Z, Gao Y, Wang L, Chen C, Wang X. TP53 Mutations Promote Immunogenic Activity in Breast Cancer. J Oncol (2019) 2019:5952836. doi: 10.1155/2019/5952836
40. Munoz-Fontela C, Mandinova A, Aaronson SA, Lee SW. Emerging Roles of P53 and Other Tumour-Suppressor Genes in Immune Regulation. Nat Rev Immunol (2016) 16:741–50. doi: 10.1038/nri.2016.99
41. Pohl SG, Brook N, Agostino M, Arfuso F, Kumar AP, Dharmarajan A. Wnt Signaling in Triple-Negative Breast Cancer. Oncogenesis (2017) 6:e310. doi: 10.1038/oncsis.2017.14
42. Lawrence RT, Perez EM, Hernandez D, Miller CP, Haas KM, Irie HY, et al. The Proteomic Landscape of Triple-Negative Breast Cancer. Cell Rep (2015) 11:630–44. doi: 10.1016/j.celrep.2015.03.050
43. Matsuda Y, Schlange T, Oakeley EJ, Boulay A, Hynes NE. WNT Signaling Enhances Breast Cancer Cell Motility and Blockade of the WNT Pathway by Sfrp1 Suppresses MDA-MB-231 Xenograft Growth. Breast Cancer Res (2009) 11:R32. doi: 10.1186/bcr2317
44. Reya T, Clevers H. Wnt Signalling in Stem Cells and Cancer. Nature (2005) 434:843–50. doi: 10.1038/nature03319
Keywords: cancer panel, next-generation sequencing, breast cancer, survival analysis, triple negative
Citation: Yang P-S, Chao Y-T, Lung C-F, Liu C-L, Chang Y-C, Li K-C and Hsu Y-C (2022) Association of Pathway Mutations With Survival in Taiwanese Breast Cancers. Front. Oncol. 12:819555. doi: 10.3389/fonc.2022.819555
Received: 21 November 2021; Accepted: 24 June 2022;
Published: 22 July 2022.
Edited by:
Giovana Tardin Torrezan, A.C.Camargo Cancer Center, BrazilReviewed by:
Wenjun Yi, Central South University, ChinaVidya Sethunath, Dana–Farber Cancer Institute, United States
Copyright © 2022 Yang, Chao, Lung, Liu, Chang, Li and Hsu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yi-Chiung Hsu, syic@ncu.edu.tw