Front. Genet.Frontiers in GeneticsFront. Genet.1664-8021Frontiers Media S.A.112597010.3389/fgene.2023.1125970GeneticsOriginal ResearchPotential predictive value of circulating tumor DNA (ctDNA) mutations for the efficacy of immune checkpoint inhibitors in advanced triple-negative breast cancerTan et al.10.3389/fgene.2023.1125970TanQiaorui1ChiYajing12SuMu3ZhouJinxing3ZhouDongdong1ZhengFangchao1ManXiaochu1SunShujuan1HuangJie1LiHuihui1*1Department of Medical Oncology, Shandong Cancer Hospital and Institute, Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan, Shandong, China2School of Medicine, Nankai University, Tianjin, China3Berry Oncology Corporation, Beijing, China
Edited by:Changgang Sun, Affiliated Hospital of Weifang Medical University, China
Reviewed by:Jiang Chang, Huazhong University of Science and Technology, China
Kan Zhai, Capital Medical University, China
*Correspondence: Huihui Li, huihuili82@163.com
This article was submitted to Cancer Genetics and Oncogenomics, a section of the journal Frontiers in Genetics
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.
Background: In recent years, tumor immunotherapy has become a viable treatment option for triple negative breast cancer (TNBC). Among these, immune checkpoint inhibitors (ICIs) have demonstrated good efficacy in advanced TNBC patients with programmed death-ligand 1 (PD-L1) positive expression. However, only 63% of PD-L1-positive individuals showed any benefit from ICIs. Therefore, finding new predictive biomarkers will aid in identifying patients who are likely to benefit from ICIs. In this study, we used liquid biopsies and next-generation sequencing (NGS) to dynamically detect changes in circulating tumor DNA (ctDNA) in the blood of patients with advanced TNBC treated with ICIs and focused on its potential predictive value.
Methods: From May 2018 to October 2020, patients with advanced TNBC treated with ICIs at Shandong Cancer Hospital were included prospectively. Patient blood samples were obtained at the pretreatment baseline, first response evaluation, and disease progression timepoints. Furthermore, 457 cancer-related genes were evaluated by NGS, and patients’ ctDNA mutations, gene mutation rates, and other indicators were determined and coupled with clinical data for statistical analysis.
Results: A total of 11 TNBC patients were included in this study. The overall objective response rate (ORR) was 27.3%, with a 6.1-month median progression-free survival (PFS) (95% confidence interval: 3.877–8.323 months). Of the 11 baseline blood samples, 48 mutations were found, with the most common mutation types being frame shift indels, synonymous single-nucleotide variations (SNVs), frame indel missenses, splicing, and stop gains. Additionally, univariate Cox regression analysis revealed that advanced TNBC patients with one of 12 mutant genes (CYP2D6 deletion and GNAS, BCL2L1, H3F3C, LAG3, FGF23, CCND2, SESN1, SNHG16, MYC, HLA-E, and MCL1 gain) had a shorter PFS with ICI treatment (p < 0.05). To some extent, dynamic changes of ctDNA might indicate the efficacy of ICIs.
Conclusion: Our data indicate that ICI efficacy in patients with advanced TNBC may be predicted by 12 mutant ctDNA genes. Additionally, dynamic alterations in peripheral blood ctDNA might be used to track the effectiveness of ICI therapy in those with advanced TNBC.
triple-negative breast cancerimmune checkpoint inhibitorsctDNAnext-generation sequencingefficacy1 Introduction
Breast cancer has become the leading cause of cancer-related morbidity and mortality in women worldwide (Bray et al., 2018; Feng et al., 2019; Siegel et al., 2021). Triple negative breast cancer (TNBC) accounts for 12%–17% of all breast cancer cases, is commonly recurrent and metastatic, and has a median overall survival (OS) of only 10–13 months (William D Foulkes et al., 2010; Schumacher and Schreiber, 2015; Garrido-Castro et al., 2019). TNBC is the most drug resistant subtype of breast cancer because of low or absent expression of the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2) (Adams et al., 2019). Despite chemotherapy remaining the most effective course of treatment, advanced TNBC patients usually have limited response to chemotherapy and develop disease progression. Therefore, it is critical to continuously examine therapeutic efficacy and seek out new treatment options.
TNBC differs from other subtypes of breast cancer in various ways, including having a higher mutation frequency, higher percentage of tumor-infiltrating lymphocytes, and 20% of TNBC cells expressing programmed death-ligand 1 (PD-L1) (Mittendorf et al., 2014; Budczies et al., 2015; Schumacher and Schreiber, 2015; Luen et al., 2016). Because of these characteristics, immune checkpoint inhibitors (ICIs) may have a higher efficacy in TNBC than in other breast cancer subtypes. To date, clinical trials have tested a range of ICIs, both independently and in combination, such as programmed cell death protein 1 (PD-1)/PD-L1 inhibitors plus chemotherapy. In the Impassion 130 study, patients with advanced TNBC who were PD-L1 positive and treated with atezolizumab and nab-paclitaxel had prolonged progression-free survival (PFS) and OS rates (Schmid et al., 2018; Iwata et al., 2020; Chen et al., 2022). Similar findings were reported in the KEYNOTE-355 study, which found that pembrolizumab in combination with chemotherapy significantly improved OS in PD-L1 positive patients with advanced TNBC compared with chemotherapy alone (Cortes et al., 2020; Cortes et al., 2022). However, ICIs are not appropriate for all advanced TNBC patients. The KEYNOTE-119 study found that in advanced TNBC patients, pembrolizumab did not significantly prolong OS compared with chemotherapy alone (Winer et al., 2021). Similarly, the IMpassion 131 study discovered that advanced TNBC patients who received paclitaxel with atezolizumab did not display significantly improved PFS or OS (Miles et al., 2021). These results suggest that current ICIs treatment strategies may not be applicable to all advanced TNBC patients. Therefore, it is necessary to identify specific biomarkers for ICI treatment effectiveness to help determine which advanced TNBC patients could possibly benefit from such therapies.
Currently, the need for real-time illness surveillance cannot be satisfied by traditional diagnostic methods like CT imaging (Zhang et al., 2020). Interestingly, there is occasionally a genetic divergence between the primary tumor cells and metastatic breast cancer cells (Itakura et al., 2015). Circulating tumor DNA (ctDNA) is DNA that has been released to blood by primary and metastatic cancer cells during apoptosis and necrosis (Swarup and Rajeswari, 2007). Compared with traditional biopsy methods, examining ctDNA levels offers unique advantages. For example, ctDNA-based approaches overcome tumor heterogeneity and provide real-time molecular data on driver genes, drug resistance genes, and clonal organization (Abstracts of Presentations at the, 2022; Sant et al., 2022). Additionally, lesions that are difficult to evaluate using imaging can be analyzed (Appierto et al., 2017; Hrebien et al., 2019). These ctDNA-based liquid biopsies, which elucidate genetic markers in patient plasma using next-generation sequencing (NGS), can supplement the current screening techniques and have proven to be a powerful option for monitoring cancer progression (Heitzer et al., 2019; Pantel and Alix-Panabieres, 2019). Research on accurate breast cancer diagnosis and treatment has also recently utilized ctDNA. Multiple studies in early breast cancer discovered that ctDNA clearance was related to improved survival after neoadjuvant treatment (Radovich et al., 2020; Magbanua et al., 2021; Ortolan et al., 2021). Additionally, Chen et al. found that ctDNA mutations offer useful information for evaluating the effectiveness of targeted therapies in breast cancer patients who are resistant to trastuzumab and chemotherapy (Chen et al., 2020). Furthermore, their findings demonstrated that gene amplifications and decreased ctDNA levels were related to higher PFS in patients with advanced TNBC after chemotherapy (Wongchenko et al., 2020; Collier et al., 2021). However, there is a lack of studies predicting the efficacy of ICIs in advanced TNBC through ctDNA detected by NGS.
In this study, we aim to identify novel biomarkers that predict the responses to ICI treatments and explore the patterns of ctDNA mutation changes in relation to immunotherapy efficacy by analyzing the ctDNA modifications and mutations in plasma of advanced TNBC patients during ICI therapy.
2 Materials and methods2.1 Patients and samples
All patients were recruited prospectively from Shandong Cancer Hospital and Institute, Jinan, Shandong, China, from May 2018 to October 2020. The study was approved by the Institutional Review Board and the Ethics Committee of Shandong Cancer Hospital and Institute. Eligible cases had progressed TNBC, which was characterized by ER and PR immunohistochemistry (IHC) staining of less than 1% of tumor cell nuclei, a HER2 IHC score of 0-1 or 2+, and fluorescence in situ hybridization was negative. Pathological examinations were performed by the pathologists in the Department of Pathology, Shandong Cancer Hospital and Institute.
All eligible patients provided written informed consent before participating in this study. Every two cycles of ICI treatments, patients were evaluated for a response using the Response Evaluation Criteria in Solid Tumors (RECIST) 1.1 criteria. Patients’ reactions to ICIs were classified as complete response (CR), partial response (PR), stable disease (SD), or progressive disease (PD) using CT scans and the RECIST 1.1 criteria. We collected dynamic peripheral blood specimens from the patients at baseline, first efficacy assessment, and disease progression. All eligible patients had baseline blood specimens and six patients had dynamic blood specimens (Figure 1). The last follow-up appointment was 29 January 2021. At the time of the final follow-up, all eligible patients were still alive. All relevant clinical data were obtained by reviewing patients’ medical records and contacting them over the phone. The time from registration to illness progression was defined as PFS.
The timepoints for collection of venous blood specimens. The numbers represent the cycles of treatment.
2.2 Plasma isolation and ctDNA library construction
Peripheral blood samples were collected in special cell-free DNA BCT tubes (Streck Laboratories, United States). Plasma supernatant and blood cell sediment were separated after whole blood samples were centrifuged at 4°C for 10 min (1,600 x g), then the process was repeated (16,000 x g). The ctDNA was isolated using the MagMAX Cell-Free DNA Isolation Kit (Thermo Fisher Scientific, United States) according to the manufacturer’s instructions. Gel electrophoresis and a Qubit® 4.0 Fluorometer (Life Technologies, United States) were used to quantify the purity of the isolated DNA. The pre-libraries were created using the previously described technique for targeted sequencing of ctDNA (Lv et al., 2015). To create a sequencing library, ctDNA fragments were collected using the NGS gene panel, which includes 457 cancer-associated genes (Berry Oncology, Peking, China) (Supplementary Table S1). The Illumina NovaSeq 6000 platform (Illumina, San Diego, CA, United States) was used to apply the sequencing libraries in 150PE mode.
2.3 Bioinformatics analysis of ctDNA mutations
FASTP (Chen et al., 2018) was used to trim adapters and obtain clean results by removing low-quality sequences. By using Burrows-Wheeler Alignment tool (BWA), the clean reads were aligned to the Ensemble GRCh37/hg19 reference genome (Li and Durbin, 2009). PCR duplications were processed by gencore (Chen et al., 2019), then consensus reads were generated. SAMtools (Li et al., 2009) was applied for the detection of single-nucleotide variations (SNVs), insertions, and loss. The ANNOVAR software (Wang et al., 2010) was used to annotate the HGVS variant description. After annotation, we eliminated non-synonymous SNVs with PopFreqMax >0.05 and retained those with VAF >0.5% or >0.1% in cancer hotspots obtained from patient databases for further investigation. Tumor mutation burden (TMB) was defined as non-synonymous coding mutations per Megabase. For copy number variants (CNVs), the gene was regarded as gain when the copy number >3 and regarded as loss when the copy number <1.
2.4 Statistical analysis
R version 3.4.0 was used for all analyses. Survival-related ctDNA was identified using univariate and multivariate Cox regression analyses. Kaplan-Meier (KM) survival curve differences were assessed by log-rank analysis. Fisher’s exact test was performed to compare the ctDNA variations among patients with advanced TNBC. Each test was two-sided, andP-values <0.05 were considered significant.
3 Results3.1 Characteristics of enrolled patients
Using the eligibility criteria, 11 advanced TNBC patients who received ICIs were enrolled in this study. Among these 11 patients, the median age at diagnosis was 46 years (range 34–67 years). Most patients (54.5%, 6/11) were diagnosed with stage III disease, while 27.3% (3/11) and 18.2% (2/11) of patients were stage II and stage I, respectively. According to the pathological classification, 10 individuals were grade III, while just one patient was grade II. Additionally, 81.8% (9/11) of patients were diagnosed with invasive ductal carcinoma. The median PFS for the 11 patients was 6.1 months (95% confidence interval (CI): 3.877–8.323 months). All eligible patients received a combination treatment of ICIs plus chemotherapy as first-line treatment (63.6%, 7/11), second-line treatment (18.2%, 2/11), and third-line treatment (18.2%, 2/11). During the treatment procedure, the best efficacy was CR in one case (18.2%), PR in two cases (18.2%), and SD in six cases (54.5%). The clinicopathologic characteristics of the patients are shown in Table 1.
Clinical and pathological features of enrolled patients (N = 11).
3.2 Mutant genes in ctDNA may predict the efficacy of ICIs
A total of 11 baseline blood specimens were analyzed by NGS, with 457 genes detected (Supplementary Table S1). We observed 48 genetic variants in the 11 baseline blood specimens, with the mutation types primarily including frame shift indels, synonymous SNVs, frame indel missenses, splicing, and stop gains. Within the SNVs, the most frequently mutated gene was TP53 (73%, 8/11), with the other SNV frequencies<20%.
Furthermore, we identified CNVs with mutation frequencies>20% (Figure 2) and found that 64% of patients have HLA-A and HLA-C gain. HLA-B gain and CYP2D6 loss were present in 45% and 27% of patients, respectively. In addition, 27% of patients have HLA-E, MCL1, MYC, NBN, PTPN1, RAD21, and SNHG16 gain. The predictive power of ctDNA in ICI treatment management was then assessed. By univariate Cox regression analysis, we found that CYP2D6 loss and GNAS, BCL2L1, H3F3C, LAG3, FGF23, CCND2, SESN1, SNHG16, MYC, HLA-E, and MCL1 gain were associated with short PFS (p < 0.05) (Figure 3). The results indicate that patients with these particular mutations might not fully benefit from ICI therapy. However, multivariate Cox regression analysis did not identify statistically significant mutated genes associated with PFS.
Genes with a mutation frequency greater than 20%.
Correlation of circulating tumor DNA (ctDNA) gene mutations with progression-free survival (PFS) in advanced triple negative breast cancer (TNBC) patients receiving immune checkpoint inhibitor (ICI) treatment.
3.3 Dynamic changes of mutant genes in ctDNA during ICI treatment
Considering that ctDNA genomic traits change dynamically with therapy (Garcia-Murillas et al., 2015), we analyzed ctDNA collected from six patients with dynamic blood samples. We also analyzed the changes in carcinoembryonic antigen (CEA), carbohydrate antigen 125 (CA125), and carbohydrate antigen 153 (CA153) levels (Figure 4).
The dynamic changes of circulating tumor DNA (ctDNA) features (copy number variants (CNVs) and single-nucleotide variations (SNVs)), tumor burden, and CEA, CA125, and CA153 levels in six triple negative breast cancer (TNBC) patients.
For patients with the best CR, PR, or SD efficacy evaluations, we observed two patterns of ctDNA mutations: 1) the copy number or mutation frequency declined at first response evaluation, while the copy number or mutation frequency rose again at disease progression; 2) the copy number or mutation frequency rose at first response evaluation, while the copy number or mutation frequency declined at disease progression. From these two patterns of mutations, we observed that during ICI treatment, the CNV or SNV initially rises and then falls or initially falls and then rises. When this occurs with changing ctDNA mutations, it possibly indicates that the patients are responding to ICI treatment. Additionally, Patient C6 experienced disease progression at the first response evaluation, at which point the mutation frequency or copy number rose. However, we did not observe this pattern of changes in CEA, CA125 or CA153 associated with ICI efficacy among these six patients.
In future studies, the sample size should be increased and additional blood collection points should be included to further verify the relationship between the dynamic changes of ctDNA and ICI treatment efficacy.
4 Discussion
ICIs are a promising treatment method that has significantly improved TNBC patient treatment outcomes (Cortes et al., 2020; Iwata et al., 2020; Emens et al., 2021; Chen et al., 2022). However, not all patients with advanced TNBC can benefit from ICIs (Miles et al., 2021; Winer et al., 2021), so timely monitoring of ICI efficacy and choosing appropriate patients are essential. Using NGS technology, we dynamically evaluated peripheral blood ctDNA alterations in patients with advanced TNBC during treatment. Concurrent tumor load and tumor marker changes in the included patients were also recorded. We found that 48 genetic mutations occurred in the 11 baseline samples. For the SNV, TP53 was the gene with the highest mutation frequency in advanced TNBC. These results were consistent with previous reports (Jiang et al., 2021; Chen et al., 2022). Furthermore, in advanced TNBC patients who received ICIs, CYP2D6 loss and GNAS, BCL2L1, H3F3C, LAG3, FGF23, CCND2, SESN1, SNHG16, MYC, HLA-E, and MCL1 gain were related to shorter PFS. Our results suggest that advanced TNBC patients with one of these 12 mutated genes are possibly not suitable candidates for ICI treatment. As shown in the literature, BCL2L1, LAG3, CCND2, SNHG16, MYC, HLA-E, MCL-1, and GNAS have regulatory effects on apoptosis, anti-tumor immunity, the tumor cell cycle, cell growth, or metastasis. Among them, BCL2L1 and MCL-1 had an inhibitory effect on apoptosis in TNBC (Kelly and Strasser, 2011; Goodwin et al., 2015). Other studies found that LAG3 and HLA-E were associated with immunosuppressive effects and tumor immune evasion in breast cancer, gastric cancer, and ovarian cancer, among others (Woo et al., 2012; Tuncel et al., 2013; Morandi et al., 2014; Zheng et al., 2015; Su et al., 2016; Andrews et al., 2017; Du et al., 2020). In addition, upregulated CCND2, SNHG16, MYC, and GNAS levels were associated with accelerating cell cycle progression, as well as promoting growth and metastasis in lung cancer, breast cancer, and colon cancer (Cai et al., 2017; Fallah et al., 2017; Teng et al., 2017; Batistatou et al., 2018; Hung et al., 2018; Han et al., 2019; Jin et al., 2019; SM, 2020). Additionally, to explore the linkage among these 12 mutated genes, we utilized the KEGG database to examine the signaling pathways in which these genes are located. We found that several of the mutant genes are implicated in numerous important signaling pathways, both upstream and downstream, including the PI3K-AKT, MAPK, Ras, and JAK-STAT signaling pathways (Christov et al., 2011; Hernandez-Varas et al., 2011; Clinkenbeard et al., 2016; Wang et al., 2017; Luo et al., 2018; Wang et al., 2018; Wu et al., 2019; Zhang et al., 2019; Chang et al., 2020; Schafer et al., 2020; Umehara et al., 2020; Ho and Bergwitz, 2021) (Figure 5). These signaling pathways are involved in the regulation of cancer-related phenotypes, including cell cycle progression, cell survival, cell migration, signal transduction, and more. Therefore, the regulatory interactions between CCND2, FGF23, MCL1, MYC, GNAS, and BCL2L1 may alter the efficacy of ICIs. We hypothesize that the effects of these mutant genes on drug efficacy result from their promoting malignant cellular phenotypes. Future cell experiments will be conducted to confirm this regulatory impact.
(A) FGF23 and BCL2L1 in the RAS signaling pathway. (B) FGF23 and MYC in the MAPK signaling pathway. (C) MCL1, BCL2L1, MYC, and CCND2 in the JAK-STAT signaling pathway. (D) CCND2, FGF23, MCL1, MYC, and BCL2L1 in the PI3K-AKT signaling pathway. All figures are included from the KEGG database.
Serum CEA, CA153, and CA125 are tumor markers that are frequently used in breast cancer screening and therapy efficacy monitoring (Cheung and Robertson, 2000; Seregni et al., 2004; Duffy, 2006; Duffy et al., 2010). However, traditional tumor markers still have some drawbacks. First, tumor load, physiological and pathological conditions of the body, basic diseases, and medicines can all alter serum tumor marker levels, resulting in lower sensitivity and specificity. Second, the critical value of serum tumor marker levels that insinuate disease progression is unclear (Seregni et al., 2004). Third, 30% of patients who respond to therapy may see a paradoxical increase in tumor marker levels following the start of chemotherapy. This rise could be attributed to the severe neoplastic cell necrosis caused by cytotoxic treatments, and marker levels can remain elevated for up to 3 months (Seregni et al., 2004). Overall, these reasons indicate that traditional tumor markers do not effectively reflect the efficacy of ICIs. Our results suggest that the changes in CEA, CA153, and CA125 levels may not be consistent with the tumor load trends. Because of the small sample size of this study, the connection between traditional tumor markers and ICI efficacy needs to be verified in a larger population.
Currently, previous research has shown that the dynamic changes of ctDNA levels correlate with treatment efficacy in hepatocellular carcinoma, pancreatic cancer, and breast cancer (Cai et al., 2019; van der Sijde et al., 2019; Liu et al., 2022). However, there is a lack of studies on the relationship between ctDNA dynamic changes and ICI efficacy. By monitoring the changes of ctDNA-related mutations at several timepoints during the patient treatment period, we found that some changes in ctDNA were similar to tumor load changes, such as pattern 1. By dynamically monitoring the changes of ctDNA-related mutations during the ICI treatment period of advanced TNBC patients, we found that the elimination of ctDNA mutations or a decreased mutation rate often represented a better ICI efficacy. In contrast, as the disease progressed, the ctDNA mutation pattern revealed that the copy number or mutation frequency rose. However, we also observed that the copy number or mutation frequency rose in patients treated with ICIs despite being effective, such as in pattern 2. Considering the source of ctDNA, we hypothesize that this is the result of chemotherapy effectively triggering tumor cell death. This causes more ctDNA to be released into the blood and allows mutations that were not previously found in the blood to be detected throughout the treatment process (Diaz and Bardelli, 2014). Our results indicate that plasma ctDNA-based liquid biopsy could potentially be used as an additional test to monitor disease in TNBC patients.
In conclusion, our data demonstrate that CYP2D6 loss and GNAS, BCL2L1, H3F3C, LAG3, FGF23, CCND2, SESN1, SNHG16, MYC, HLA-E, and MCL1 gains were predictive of ICI effectiveness in patients with advanced TNBC. In addition, dynamic monitoring of ctDNA in patients with advanced TNBC might provide a timely indicator of sensitivity to ICI treatment. In the future, further studies are necessary to validate our observations.
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
The studies involving human participants were reviewed and approved by Institutional Review Board and the Ethics Committee of Shandong Cancer Hospital and Institute. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
Conceptualization: HL, article collection and analysis: QT, YC, MS, JZ, DZ, FZ, XM, SS, and JH.
Funding
We would like to acknowledge the funding support of the Start-up fund of Shandong Cancer Hospital (Grant No. 2020-PYB10), Beijing Science and Technology Innovation Fund (Grant No. KC2021-ZZ-0010-1), and Chinese anti-cancer association (Grant No. 2021HZ21).
We thank J. Iacona, Ph.D., from Liwen Bianji (Edanz) (www.liwenbianji.cn) for editing the English text of a draft of this manuscript.
Conflict of interest
Authors MS and JZ were employed by Berry Oncology Corporation.
The remaining 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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1125970/full#supplementary-material
ReferencesAbstracts of Presentations at the association of clinical scientists 143rd meeting louisville, KY may 11-14,2022. Ann. Clin. Lab. Sci. 2022;52(3):511–525.AdamsS.DiamondJ. R.HamiltonE.PohlmannP. R.TolaneyS. M.ChangC. W. (2019). Atezolizumab plus nab-paclitaxel in the treatment of metastatic triple-negative breast cancer with 2-year survival follow-up: A phase 1b clinical trial. JAMA Oncol.5 (3), 334–342. 10.1001/jamaoncol.2018.5152AndrewsL. P. M. A.DrakeC. G.VignaliD. A. (2017). LAG3 (CD223) as a cancer immunotherapy target. Immunol. Rev.276 (1), 80–96. 10.1111/imr.12519AppiertoV.Di CosimoS.ReduzziC.PalaV.CappellettiV.DaidoneM. G. (2017). How to study and overcome tumor heterogeneity with circulating biomarkers: The breast cancer case. Semin. Cancer Biol.44, 106–116. 10.1016/j.semcancer.2017.04.007BatistatouA.KotoulaV.BobosM.KouvatseasG.ZagouriF.TsolakiE. (2018). Correlation of MYC gene and protein status with breast cancer subtypes and outcome of patients treated with anthracycline-based adjuvant chemotherapy. Pooled analysis of 2 hellenic cooperative group phase III trials. Clin. Breast Cancer18 (1), 53–62. 10.1016/j.clbc.2017.07.004BrayF.FerlayJ.SoerjomataramI.SiegelR. L.TorreL. A.JemalA. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin.68 (6), 394–424. 10.3322/caac.21492BudcziesJ.BockmayrM.DenkertC.KlauschenF.LennerzJ. K.GyorffyB. (2015). Classical pathology and mutational load of breast cancer - integration of two worlds. J. Pathol. Clin. Res.1 (4), 225–238. 10.1002/cjp2.25CaiC.HuoQ.WangX.ChenB.YangQ. (2017). SNHG16 contributes to breast cancer cell migration by competitively binding miR-98 with E2F5. Biochem. Biophys. Res. Commun.485 (2), 272–278. 10.1016/j.bbrc.2017.02.094CaiZ.ChenG.ZengY.DongX.LiZ.HuangY. (2019). Comprehensive liquid profiling of circulating tumor DNA and protein biomarkers in long-term follow-up patients with hepatocellular carcinoma. Clin. Cancer Res.25 (17), 5284–5294. 10.1158/1078-0432.CCR-18-3477ChangY. S.ChangC. M.LinC. Y.ChaoD. S.HuangH. Y.ChangJ. G. (2020). Pathway mutations in breast cancer using whole-exome sequencing. Oncol. Res.28 (2), 107–116. 10.3727/096504019X15698362825407ChenL.JiangY. Z.WuS. Y.WuJ.DiG. H.LiuG. Y. (2022). Famitinib with camrelizumab and nab-paclitaxel for advanced immunomodulatory triple-negative breast cancer (FUTURE-C-Plus): An open-label, single-arm, phase II trial. Clin. Cancer Res.: J. Am. Ass. Cancer Res.28 (13), 2807–2817. 10.1158/1078-0432.CCR-21-4313ChenS.ZhouY.ChenY.GuJ. (2018). fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics34 (17), i884–i890. 10.1093/bioinformatics/bty560ChenS.ZhouY.ChenY.HuangT.LiaoW.XuY. (2019). Gencore: An efficient tool to generate consensus reads for error suppressing and duplicate removing of NGS data. BMC Bioinforma.20 (23), 606. 10.1186/s12859-019-3280-9ChenZ.SunT.YangZ.ZhengY.YuR.WuX. (2020). Monitoring treatment efficacy and resistance in breast cancer patients via circulating tumor DNA genomic profiling. Mol. Genet. Genomic Med.8 (2), e1079. 10.1002/mgg3.1079CheungK. L. G. C.RobertsonJ. F. (2000). Tumour marker measurements in the diagnosis and monitoring of breast cancer. Cancer Treat. Rev.26 (2), 91–102. 10.1053/ctrv.1999.0151ChristovM.KorenS.YuanQ.BaronR.LanskeB. (2011). Genetic ablation of sfrp4 in mice does not affect serum phosphate homeostasis. Endocrinology152 (5), 2031–2036. 10.1210/en.2010-1351ClinkenbeardE. L.CassT. A.NiP.HumJ. M.BellidoT.AllenM. R. (2016). Conditional deletion of murine Fgf23: Interruption of the normal skeletal responses to phosphate challenge and rescue of genetic hypophosphatemia. J. Bone Min. Res.31 (6), 1247–1257. 10.1002/jbmr.2792CollierK. A.AsadS.TallmanD.JenisonJ.RajkovicA.MardisE. R. (2021). Association of 17q22 amplicon via cell-free DNA with platinum chemotherapy response in metastatic triple-negative breast cancer. JCO Precis. Oncol.5, PO.21.00104. 10.1200/PO.21.00104CortesJ.CesconD. W.RugoH. S.NoweckiZ.ImS-A.YusofM. M. (2020). Pembrolizumab plus chemotherapy versus placebo plus chemotherapy for previously untreated locally recurrent inoperable or metastatic triple-negative breast cancer (KEYNOTE-355): A randomised, placebo-controlled, double-blind, phase 3 clinical trial. Lancet396 (10265), 1817–1828. 10.1016/S0140-6736(20)32531-9CortesJ.RugoH. S.CesconD. W.ImS. A.YusofM. M.GallardoC. (2022). Pembrolizumab plus chemotherapy in advanced triple-negative breast cancer. N. Engl. J. Med.387 (3), 217–226. 10.1056/NEJMoa2202809DiazL. A.Jr.BardelliA. (2014). Liquid biopsies: Genotyping circulating tumor DNA. J. Clin. Oncol. official J. Am. Soc. Clin. Oncol.32 (6), 579–586. 10.1200/JCO.2012.45.2011DuH.YiZ.WangL.LiZ.NiuB.RenG. (2020). The co-expression characteristics of LAG3 and PD-1 on the T cells of patients with breast cancer reveal a new therapeutic strategy. Int. Immunopharmacol.78, 106113. 10.1016/j.intimp.2019.106113DuffyM. J.EvoyD.McDermottE. W. (2010). CA 15-3: Uses and limitation as a biomarker for breast cancer. Clin. Chim. Acta411 (23-24), 1869–1874. 10.1016/j.cca.2010.08.039DuffyM. J. (2006). Serum tumor markers in breast cancer: Are they of clinical value?Clin. Chem.52 (3), 345–351. 10.1373/clinchem.2005.059832EmensL. A.AdamsS.BarriosC. H.DierasV.IwataH.LoiS. (2021). First-line atezolizumab plus nab-paclitaxel for unresectable, locally advanced, or metastatic triple-negative breast cancer: IMpassion130 final overall survival analysis. Ann. Oncol.32 (8), 983–993. 10.1016/j.annonc.2021.05.355FallahY.BrundageJ.AllegakoenP.Shajahan-HaqA. N. (2017). MYC-driven pathways in breast cancer subtypes. Biomolecules7 (3), 53. 10.3390/biom7030053FengR. M.ZongY. N.CaoS. M.XuR. H. (2019). Current cancer situation in China: Good or bad news from the 2018 global cancer statistics?Cancer Commun. (Lond)39 (1), 22. 10.1186/s40880-019-0368-6Garcia-MurillasI.SchiavonG.WeigeltB.NgC.HrebienS.CuttsR. J. (2015). Mutation tracking in circulating tumor DNA predicts relapse in early breast cancer. Sci. Transl. Med.7 (302), 302ra133. 10.1126/scitranslmed.aab0021Garrido-CastroA. C.LinN. U.PolyakK. (2019). Insights into molecular classifications of triple-negative breast cancer: Improving patient selection for treatment. Cancer Discov.9 (2), 176–198. 10.1158/2159-8290.CD-18-1177GoodwinC. M.RossaneseO. W.OlejniczakE. T.FesikS. W. (2015). Myeloid cell leukemia-1 is an important apoptotic survival factor in triple-negative breast cancer. Cell Death Differ.22 (12), 2098–2106. 10.1038/cdd.2015.73HanW.WangL.ZhangL.WangY.LiY. (2019). Circular RNA circ-RAD23B promotes cell growth and invasion by miR-593-3p/CCND2 and miR-653-5p/TIAM1 pathways in non-small cell lung cancer. Biochem. Biophys. Res. Commun.510 (3), 462–466. 10.1016/j.bbrc.2019.01.131HeitzerE.HaqueI. S.RobertsC. E. S.SpeicherM. R. (2019). Current and future perspectives of liquid biopsies in genomics-driven oncology. Nat. Rev. Genet.20 (2), 71–88. 10.1038/s41576-018-0071-5Hernandez-VarasP.ColoG. P.BartolomeR. A.PatersonA.Medrano-FernandezI.Arellano-SanchezN. (2011). Rap1-GTP-interacting adaptor molecule (RIAM) protein controls invasion and growth of melanoma cells. J. Biol. Chem.286 (21), 18492–18504. 10.1074/jbc.M110.189811HoB. B. B. C.BergwitzC. (2021). FGF23 signalling and physiology. J. Mol. Endocrinol.66 (2), R23–R32. 10.1530/JME-20-0178HrebienS.CitiV.Garcia-MurillasI.CuttsR.FenwickK.KozarewaI. (2019). Early ctDNA dynamics as a surrogate for progression-free survival in advanced breast cancer in the BEECH trial. Ann. Oncol.:J. Euro. Soc. Med. Oncol.30 (6), 945–952. 10.1093/annonc/mdz085HungC. S.WangS. C.YenY. T.LeeT. H.WenW. C.LinR. K. (2018). Hypermethylation of CCND2 in lung and breast cancer is a potential biomarker and drug target. Int. J. Mol. Sci.19 (10), 3096. 10.3390/ijms19103096ItakuraH.AchrolA. S.MitchellL. A.LoyaJ. J.LiuT.WestbroekE. M. (2015). Magnetic resonance image features identify glioblastoma phenotypic subtypes with distinct molecular pathway activities. Sci. Transl. Med.7 (303), 303ra138. 10.1126/scitranslmed.aaa7582IwataH.EmensL.AdamsS.BarriosC. H.DiérasV.LoiS. (2020). 49MO IMpassion130: Final OS analysis from the pivotal phase III study of atezolizumab + nab-paclitaxel vs placebo + nab-paclitaxel in previously untreated locally advanced or metastatic triple-negative breast cancer. Ann. Oncol.31, S1261–S1262. 10.1016/j.annonc.2020.10.069JiangY. Z.LiuY.XiaoY.HuX.JiangL.ZuoW. J. (2021). Molecular subtyping and genomic profiling expand precision medicine in refractory metastatic triple-negative breast cancer: The FUTURE trial. Cell Res.31 (2), 178–186. 10.1038/s41422-020-0375-9JinX.ZhuL.CuiZ.TangJ.XieM.RenG. (2019). Elevated expression of GNAS promotes breast cancer cell proliferation and migration via the PI3K/AKT/Snail1/E-cadherin axis. Clin. Transl. Oncol.21 (9), 1207–1219. 10.1007/s12094-019-02042-wKellyG. L.StrasserA. (2011). The essential role of evasion from cell death in cancer. Adv. Cancer Res.111, 39–96. 10.1016/B978-0-12-385524-4.00002-7LiH.DurbinR. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics25 (14), 1754–1760. 10.1093/bioinformatics/btp324LiH.HandsakerB.WysokerA.FennellT.RuanJ.HomerN. (2009). The sequence alignment/map format and SAMtools. Bioinformatics25 (16), 2078–2079. 10.1093/bioinformatics/btp352LiuB.HuZ.RanJ.XieN.TianC.TangY. (2022). The circulating tumor DNA (ctDNA) alteration level predicts therapeutic response in metastatic breast cancer: Novel prognostic indexes based on ctDNA. Breast65, 116–123. 10.1016/j.breast.2022.07.010LuenS.VirassamyB.SavasP.SalgadoR.LoiS. (2016). The genomic landscape of breast cancer and its interaction with host immunity. Breast29, 241–250. 10.1016/j.breast.2016.07.015LuoJ.YaoJ. F.DengX. F.ZhengX. D.JiaM.WangY. Q. (2018). 14, 15-EET induces breast cancer cell EMT and cisplatin resistance by up-regulating integrin αvβ3 and activating FAK/PI3K/AKT signaling. J. Exp. Clin. Cancer Res.37 (1), 23. 10.1186/s13046-018-0694-6LvW.WeiX.GuoR.LiuQ.ZhengY.ChangJ. (2015). Noninvasive prenatal testing for Wilson disease by use of circulating single-molecule amplification and resequencing technology (cSMART). Clin. Chem.61 (1), 172–181. 10.1373/clinchem.2014.229328MagbanuaM. J. M.SwigartL. B.WuH. T.HirstG. L.YauC.WolfD. M. (2021). Circulating tumor DNA in neoadjuvant-treated breast cancer reflects response and survival. Ann. Oncol.32 (2), 229–239. 10.1016/j.annonc.2020.11.007MilesD.GligorovJ.AndreF.CameronD.SchneeweissA.BarriosC. (2021). Primary results from IMpassion131, a double-blind, placebo-controlled, randomised phase III trial of first-line paclitaxel with or without atezolizumab for unresectable locally advanced/metastatic triple-negative breast cancer. Ann. Oncol.32 (8), 994–1004. 10.1016/j.annonc.2021.05.801MittendorfE. A.PhilipsA. V.Meric-BernstamF.QiaoN.WuY.HarringtonS. (2014). PD-L1 expression in triple-negative breast cancer. Cancer Immunol. Res.2 (4), 361–370. 10.1158/2326-6066.CIR-13-0127MorandiF.FainardiE.RizzoR.Rouas-FreissN. (2014). The role of HLA-class Ib molecules in immune-related diseases, tumors, and infections. J. Immunol. Res.2014, 231618. 10.1155/2014/231618OrtolanE.AppiertoV.SilvestriM.MiceliR.VeneroniS.FolliS. (2021). Blood-based genomics of triple-negative breast cancer progression in patients treated with neoadjuvant chemotherapy. ESMO Open6 (2), 100086. 10.1016/j.esmoop.2021.100086PantelK.Alix-PanabieresC. (2019). Liquid biopsy and minimal residual disease - latest advances and implications for cure. Nat. Rev. Clin. Oncol.16 (7), 409–424. 10.1038/s41571-019-0187-3RadovichM.JiangG.HancockB. A.ChitambarC.NandaR.FalksonC. (2020). Association of circulating tumor DNA and circulating tumor cells after neoadjuvant chemotherapy with disease recurrence in patients with triple-negative breast cancer: Preplanned secondary analysis of the BRE12-158 randomized clinical trial. JAMA Oncol.6 (9), 1410–1415. 10.1001/jamaoncol.2020.2295SantM.Bernat-PegueraA.FelipE.MargeliM.BentoM. J.ArdanazE. (2022). Association of sentinel node biopsy and pathological report completeness with survival benefit for cutaneous melanoma and factors influencing their different uses in European populations. Cancers (Basel)14 (2), 4379. 10.3390/cancers14184379SchaferJ. M.LehmannB. D.Gonzalez-EricssonP. I.MarshallC. B.BeelerJ. S.RedmanL. N. (2020). Targeting MYCN-expressing triple-negative breast cancer with BET and MEK inhibitors. Sci. Transl. Med.12 (534), eaaw8275. 10.1126/scitranslmed.aaw8275SchmidP.AdamsS.RugoH. S.SchneeweissA.BarriosC. H.IwataH. (2018). Atezolizumab and nab-paclitaxel in advanced triple-negative breast cancer. N. Engl. J. Med.379 (22), 2108–2121. 10.1056/NEJMoa1809615SchumacherT. N. S. R.SchreiberR. D. (2015). Neoantigens in cancer immunotherapy. Science348 (6230), 69–74. 10.1126/science.aaa4971SeregniE.ColiA.MazzuccaN.Italian Group Ria-Irma Test IAoNM (2004). Circulating tumour markers in breast cancer. Eur. J. Nucl. Med. Mol. Imaging31 (1), S15–S22. 10.1007/s00259-004-1523-zSiegelR. L.MillerK. D.FuchsH. E.JemalA. (2021). Cancer statistics, 2021. CA Cancer J. Clin.71 (1), 7–33. 10.3322/caac.21654SmD. (2020). The SNHG16/miR-30a axis promotes breast cancer cell proliferation and invasion by regulating RRM2. Neoplasma67 (3), 567–575. 10.4149/neo_2020_190625N550SuM.HuangC. X.DaiA. P. (2016). Immune checkpoint inhibitors: Therapeutic tools for breast cancer. Asian Pac J. Cancer Prev.17 (3), 905–910. 10.7314/apjcp.2016.17.3.905SwarupV.RajeswariM. R. (2007). Circulating (cell-free) nucleic acids--a promising, non-invasive tool for early detection of several human diseases. FEBS Lett.581 (5), 795–799. 10.1016/j.febslet.2007.01.051TengY.RenY.HuX.MuJ.SamykuttyA.ZhuangX. (2017). MVP-mediated exosomal sorting of miR-193a promotes colon cancer progression. Nat. Commun.8, 14448. 10.1038/ncomms14448TuncelT.KaragozB.HaholuA.OzgunA.EmirzeogluL.BilgiO. (2013). Immunoregulatory function of HLA-G in gastric cancer. Asian Pac J. Cancer Prev.14 (12), 7681–7684. 10.7314/apjcp.2013.14.12.7681UmeharaT.KagawaS.TomidaA.MuraseT.AbeY.ShinguK. (2020). Body temperature-dependent microRNA expression analysis in rats: rno-miR-374-5p regulates apoptosis in skeletal muscle cells via Mex3B under hypothermia. Sci. Rep.10 (1), 15432. 10.1038/s41598-020-71931-wvan der SijdeF.VietschE. E.MustafaD. A. M.BesselinkM. G.Groot KoerkampB.van EijckC. H. J. (2019). Circulating biomarkers for prediction of objective response to chemotherapy in pancreatic cancer patients. Cancers (Basel)11 (1), 93. 10.3390/cancers11010093WangK.LiM.HakonarsonH. (2010). Annovar: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res.38 (16), e164. 10.1093/nar/gkq603WangR.SuC.WangX.FuQ.GaoX.ZhangC. (2018). Global gene expression analysis combined with a genomics approach for the identification of signal transduction networks involved in postnatal mouse myocardial proliferation and development. Int. J. Mol. Med.41 (1), 311–321. 10.3892/ijmm.2017.3234WangW.JiangW.HouL.DuanH.WuY.XuC. (2017). Weighted gene co-expression network analysis of expression data of monozygotic twins identifies specific modules and hub genes related to BMI. BMC Genomics18 (1), 872. 10.1186/s12864-017-4257-6William D FoulkesI. E. S.Reis-FilhoJorge S.Reis-FilhoJ. S. (2010). Triple-negative breast cancer. N. Engl. J. Med.363 (20), 1938–1948. 10.1056/NEJMra1001389WinerE. P.LipatovO.ImS-A.GoncalvesA.Muñoz-CouseloE.LeeK. S. (2021). Pembrolizumab versus investigator-choice chemotherapy for metastatic triple-negative breast cancer (KEYNOTE-119): A randomised, open-label, phase 3 trial. Lancet Oncol.22 (4), 499–511. 10.1016/S1470-2045(20)30754-3WongchenkoM. J.KimS‐.B.SauraC.OliveiraM.LipsonD.KennedyM. (2020). Circulating tumor DNA and biomarker analyses from the LOTUS randomized trial of first-line ipatasertib and paclitaxel for metastatic triple-negative breast cancer. JCO Precis. Oncol.4, 1012–1024. 10.1200/PO.19.00396WooS. R.TurnisM. E.GoldbergM. V.BankotiJ.SelbyM.NirschlC. J. (2012). Immune inhibitory molecules LAG-3 and PD-1 synergistically regulate T-cell function to promote tumoral immune escape. Cancer Res.72 (4), 917–927. 10.1158/0008-5472.CAN-11-1620WuD. M.HongX. W.WenX.HanX. R.WangS.WangY. J. (2019). MCL1 gene silencing promotes senescence and apoptosis of glioma cells via inhibition of the PI3K/Akt signaling pathway. IUBMB Life71 (1), 81–92. 10.1002/iub.1944ZhangC.ChenZ.ChongX.ChenY.WangZ.YuR. (2020). Clinical implications of plasma ctDNA features and dynamics in gastric cancer treated with HER2-targeted therapies. Clin. Transl. Med.10 (8), e254. 10.1002/ctm2.254ZhangZ.ZhaiW.LiangJ.ChenZ.MaM.ZhaoY. (2019). Mutual inhibitions between epidermal growth factor receptor signaling and miR-124a control pancreatic progenitor proliferation. J. Cell Physiol.234 (8), 12978–12988. 10.1002/jcp.27967ZhengH.LuR.XieS.WenX.WangH.GaoX. (2015). Human leukocyte antigen-E alleles and expression in patients with serous ovarian cancer. Cancer Sci.106 (5), 522–528. 10.1111/cas.12641