Next Article in Journal
Regulation of Autophagy via Carbohydrate and Lipid Metabolism in Cancer
Next Article in Special Issue
Recent Advances with Precision Medicine Treatment for Breast Cancer including Triple-Negative Sub-Type
Previous Article in Journal
Artificial Intelligence Applied to Colonoscopy: Is It Time to Take a Step Forward?
Previous Article in Special Issue
CDK4/6 Inhibitors Overcome Endocrine ESR1 Mutation-Related Resistance in Metastatic Breast Cancer Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Meta-Analysis of Triple-Negative Breast Cancer Response to Neoadjuvant Chemotherapy

1
Department of Public Health Sciences, University of Miami Miller School of Medicine, Miami, FL 33136, USA
2
California Academy of Mathematics and Science, 1000 E Victoria St, Carson, CA 90747, USA
3
Sylvester Comprehensive Cancer Center, University of Miami Miller School of Medicine, Miami, FL 33136, USA
4
Department of Medicine, Vanderbilt University Medical Center, Nashville, TN 37232, USA
5
Vanderbilt-Ingram Cancer Center, Vanderbilt University Medical Center, Nashville, TN 37232, USA
*
Authors to whom correspondence should be addressed.
These authors jointly supervised this work: Brian D. Lehmann, X. Steven Chen.
Cancers 2023, 15(8), 2194; https://doi.org/10.3390/cancers15082194
Submission received: 4 March 2023 / Revised: 1 April 2023 / Accepted: 4 April 2023 / Published: 7 April 2023
(This article belongs to the Special Issue Precision Medicine in Breast Cancer Treatment)

Abstract

:

Simple Summary

This study aimed to identify genes associated with neoadjuvant chemotherapy (NAC) response and disease-free survival (DFS) of triple-negative breast cancer (TNBC) patients. We conducted a large-scale meta-analysis of gene expression data from multiple TNBC cohorts and analyzed four gene expression quadrants from an integrated analysis of NAC response and clinical outcomes. The study also highlighted the importance of incorporating both short-term and long-term clinical outcomes in NAC response evaluation. The findings provide insights into the complex molecular mechanisms underlying TNBC and may guide the development of personalized treatment strategies for TNBC patients.

Abstract

Triple-negative breast cancer (TNBC) is a heterogeneous disease with varying responses to neoadjuvant chemotherapy (NAC). The identification of biomarkers to predict NAC response and inform personalized treatment strategies is essential. In this study, we conducted large-scale gene expression meta-analyses to identify genes associated with NAC response and survival outcomes. The results showed that immune, cell cycle/mitotic, and RNA splicing-related pathways were significantly associated with favorable clinical outcomes. Furthermore, we integrated and divided the gene association results from NAC response and survival outcomes into four quadrants, which provided more insights into potential NAC response mechanisms and biomarker discovery.

1. Introduction

Triple-negative breast cancer (TNBC) refers to a heterogeneous collection of tumors that lack expression of the estrogen receptor (ER), progesterone receptor (PR), and HER2 amplification. TNBC is more prevalent among younger patients and African American women, representing around 15% to 20% of all breast cancers [1,2]. The estimated number of newly diagnosed TNBC cases is around 60,000 each year in the United States, which is comparable to, or even larger than, patient cohorts for many other cancer types, such as pancreatic, rectal, prostate, and ovarian cancers [3]. TNBC tumors are generally larger in size, are of higher grade, have lymph node involvement at diagnosis, and are biologically more aggressive [4]. Despite having higher rates of clinical response to pre-surgical (neo-adjuvant) chemotherapy, TNBC patients have a higher rate of distant recurrence and a poorer prognosis than women with other breast cancer subtypes [1,4]. Fewer than 30% of women with metastatic TNBC survive five years, and almost all die of their disease despite adjuvant chemotherapy [1].
The molecular landscape of TNBC has been extensively studied, and molecular subtypes of TNBC based on transcriptome data have been identified and utilized by the research community [5,6,7,8,9,10,11]. However, the targeted treatment strategies based on molecular subtypes requires additional clinical trials to prove their efficacy [12]. The emerging therapeutic options, including immune checkpoint inhibitors (ICIs) [13,14], poly-ADP-ribosyl polymerase (PARP) inhibitors [15,16], and antibody–drug conjugates (ADCs) [17,18], are available for metastatic TNBC only at the current stage, and the standard neoadjuvant treatment for non-metastatic TNBC is still conventional chemotherapy [12]. Treatment of TNBC often includes neoadjuvant chemotherapy with a combination of anthracyclines, alkylating agents, platinum salts, and taxanes, which has proven to be effective in the treatment of TNBC for a subset of patients. Patients who experience a pathologic complete response (pCR) and show no evidence of disease in the surgical specimen following neoadjuvant chemotherapy (NAC) have significant improvements in both disease-free and overall survival compared with patients with the residual invasive disease [19], while those patients with the residual disease have a much poorer prognosis and are 6 times more likely to have a recurrence and 12 times more likely to die [20,21].
Neoadjuvant chemotherapy is the standard for newly diagnosed TNBC. While ~30% of patients respond and have favorable outcomes, there are a large portion of patients who must endure the toxicities associated with chemotherapy. The ability to identify those patients unlikely to respond could give clinicians guidance to enroll patients on alternative compounds in clinical trials. Currently, there is no effective way to identify these TNBC patients. Gene signatures have been made to guide chemotherapy decisions for breast cancer, such as MammaPrint, but all TNBC samples fall into the poor prognosis category and are predicted to be responsive to NAC [22]. A series of gene signatures have also been developed to predict the pathologic complete response (pCR) of NAC in TNBC [23,24,25,26,27,28,29,30,31]. Although some meaningful patterns have emerged, such as the association of immune-related genes with pCR, most of the prediction performances were moderate and difficult to be implemented for clinical usage. The gene signatures and performances were also inconsistent across the studies since most of the studies were based on a few TNBC transcriptome datasets. Thus, it is still critical to have a deeper understanding of the mechanisms behind the responses of NAC of TNBC to develop more efficient experimental design and analysis strategies for TNBC biomarker discovery.
Gene expression meta-analysis provides a valuable opportunity to integrate multiple datasets from different sources, increasing statistical power and accuracy of results. In this study, we conducted a large-scale gene expression meta-analysis of TNBC NAC responses based on 1739 samples from 16 datasets, the largest gene expression data collection for TNBC with NAC responses to our knowledge. We also performed another gene expression meta-analysis for relapse-free survival outcomes of TNBC using 7 out of 16 datasets. By integrating meta-analysis results from NAC and survival outcomes, we were able to dissect the TNBC cohort into four new NAC response subgroups and identify the genes and pathways associated with each subgroup.
Our study provides a comprehensive overview of the gene expression changes associated with TNBC NAC responses and can help identify new treatment prediction biomarkers and therapeutic targets for TNBC patients. The results of this study could potentially lead to personalized treatment approaches for TNBC patients, improving their survival outcomes and overall quality of life.

2. Materials and Methods

2.1. Study Workflow

The study design and workflow are illustrated in Figure 1. We collected TNBC transcriptome data with NAC responses, including pathologic complete response (pCR) or residual disease (RD) outcomes, from 16 microarray/RNA-seq datasets. After preprocessing the microarray and RNA-seq data, we conducted differential testing for pCR outcome and association testing for survival outcome for each dataset. We then congregated the test results from each individual dataset using meta-analysis. The significant genes identified from pCR meta-analysis and survival meta-analysis were further divided into four subgroups. Finally, we performed functional enrichment analysis and gene network analysis for each subgroup.

2.2. Gene Expression Datasets

A total of 16 publicly available neoadjuvant chemotherapy in breast cancer cohorts were selected for the study and listed in Table 1. The study analyzed a total of 1739 TNBC samples from the 16 studies. Seven cohorts contained pCR status and breast cancer disease-free survival information, while nine cohorts only contained pCR status. Since overall survival information was not available in most datasets, disease-free survival (DFS), recurrence-free survival, relapse-free survival, and distant disease-free survival were used as endpoints and treated equivalently in the survival analysis, with DFS used to represent them throughout the analysis. Among the 16 datasets, 5 were RNA-seq datasets generated by the Illumina platform, and 11 were microarray datasets generated by the Affymetrix platform. Except for the VUMC, all other datasets were downloaded from GEO using the GEOquery R package except for the VUMC dataset [32], which can be found in Supplementary Table S1. Table 1 provides detailed information on each dataset, including pCR and DFS information, references to the study, and the platform used.

2.3. Data Pre-Processing

All microarray datasets were first normalized by the frozen robust multi-array analysis (fRMA) method [33], which performed background correction, normalization, and summarization, and removed the individual batch effect. This method was implemented in the fRMA R package. Samples that were not TNBC were filtered out, and for datasets lacking TNBC information, TNBCs were determined based on the status of estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2). For datasets lacking ER, PR, and HER2 information, we used a two-component Gaussian mixture model with parameters estimated by R optim function to approximate the distribution of ER, PR, and HER2. The three Affymetrix probes 205225_at, 208305_at, and 216836_s_at were used to represent ER, PR, and HER2 gene expressions, respectively [34]. The posterior probability of negative expression status of ER, PR, and HER2 can be calculated by a mixture Gaussian model, respectively. We determined the negative status of ER, PR, and HER2 with their posterior probabilities larger than 0.5 [5,8]. Two microarray datasets were determined by the approximate distribution of ER, PR, and HER2, and other microarray datasets either included the status of ER, PR, and HER2, or provided TNBC information. The genes in the microarray datasets were represented by the probes with the largest interquartile range (IQR) statistics using findLargest function in genefilter R package.
For the RNA-seq datasets, the counts and FPKM values were first converted to TPM values. Genes with a sum of TPM value of 0 were excluded from the datasets, and all genes in the datasets were then normalized by log2 transformation.

2.4. Association Testing for Individual Dataset

To assess the association of gene expression with NAC response (pCR or RD), we performed logistic regression models with the probability of pCR as the outcome and genomic features as the main independent variable for the individual cohort. Age was adjusted in the model as the clinical variable. All logistics regression models were built using glm function in R. In the survival analysis, Cox proportional-hazard (PH) regression models implemented with coxph function in survival R package were used to test the association between gene expression and DFS outcome, including event and event-free in years. We also adjusted the clinical variable age in the Cox models.

2.5. Meta-Analysis

To integrate individual test results across different cohorts, we used the inverse-variance weighted random-effects model, which was implemented in the meta R package. The random-effects model was chosen because of the significant heterogeneities of cross-study datasets [35]. Association test results of pCR and survival analysis were used to perform meta-analysis separately. Specifically, in pCR meta-analysis, we used the coefficients and standard errors obtained from logistic regression models as estimates and standard errors of estimates as input to the meta-analysis model. The estimated effect sizes and standard errors from pCR meta-analysis were then re-scaled to compute odds ratios for every one-unit increase in expression values. In the survival meta-analysis, we used the log hazard ratios and standard errors from Cox PH regression models as input. Similarly, the hazard ratios were computed by the estimated effect sizes and standard errors from the survival meta-analysis. We included genes present in least three cohorts that appeared in the pCR and survival meta-analysis, respectively.
To select the significant genes associated with pCR and DFS outcome, we first divided the genes into four quadrants based on the results from pCR meta-analysis and survival meta-analysis illustrated in Figure 1. In gene list quadrant one, we selected genes with positive effect sizes in pCR meta-analysis and with negative effect sizes in survival meta-analysis. This quadrant represents a favorable DFS outcome with pCR status. Similarly, quadrants two, three, and four represent the genes associated with favorable DFS outcomes with RD status, unfavorable DFS outcomes with RD status, and unfavorable DFS outcomes with pCR status, respectively.

2.6. Pathway and Functional Enrichment Analysis

Gene set enrichment analysis (GSEA) implemented in Bioconductor/R package fgsea was applied for the pCR and survival meta-analysis gene lists [36]. The functional enrichment analyses were performed for gene lists in each of the four quadrants, respectively, using over-representation analysis [37], which was implemented in clusterProfiler 4.0 Bioconductor/R package [38]. Briefly, this method determined whether biological processes that were over-represented in the gene list of interest using p-values calculated by hypergeometric distribution and adjusted by the Benjamini–Hochberg (BH) method to calculate false discovery rate (FDR). The pathway collection for both GSEA and functional enrichment analysis is BioCarta subset, KEGG subset, PID subset, Reactome subset, and WikiPathways subset of Canonical pathways (C2:CP) from the Molecular Signatures Database (MSigDB) [39], which can be accessed by the msigdbr Bioconductor/R package. The minimum and maximum sizes of gene sets used for analysis were set to 10 and 500, respectively. FDR 0.25 was used for the significant threshold for both GSEA and functional enrichment analysis.

2.7. Protein–Protein Interaction (PPI) Network Analysis

To assess the interactive associations among the four-quadrant gene lists at the protein level, we performed PPI network analysis using NetworkAnalyst webtool (version 3.0) [40], which can be accessed at http://www.networkanalyst.ca (accessed on 2 February 2023). NetworkAnalyst provides degree and betweenness to help identify the highly interactive hub nodes. The degree is used to measure the number of connections of a node to other nodes, and the betweenness measures the number of shortest paths to a node. We uploaded gene lists to the webtool and performed the PPI network analysis with International Molecular Exchange (IMEx) Interactome database [41]. The networks were built either based on the original seed proteins through the first-order network or the minimum network by trimming the first-order network to keep the essential nodes that connected the seed nodes.

3. Results

3.1. Meta-Analysis Results of NAC Response

We performed a meta-analysis based on association testing results from 16 datasets to assess the association of gene expression with TNBC NAC response (pCR vs. RD). Table 2 lists the top 10 positively associated genes with pCR (upregulated in pCR group) and the top 10 negatively associated genes with pCR (upregulated in RD group), while the complete gene list of pCR meta-analysis can be found in Supplementary Table S1. We observed a significantly unbalanced association pattern, with 1656 genes positively associated with pCR and only 70 genes negatively associated with pCR out of 1726 genes passing the FDR threshold of 0.05. Figure 2 displays the top 10 pathways with positive normalized enrichment score (NES) (positively associated with pCR) and the top 10 pathways with negative NES. Our pathway analysis suggests that cell cycle/mitotic and immune-related pathways are associated with pCR, while metabolism-related pathways are associated with RD. The complete GSEA analysis results can be found in Supplementary Table S2.

3.2. Meta-Analysis Results of DFS

We also performed a meta-analysis based on association testing results with TNBC DFS survival outcomes from seven datasets. Table 3 lists the top 10 positively associated genes with survival (upregulated gene associated with higher survival risk) and the top 10 negatively associated genes with survival (upregulated gene associated with lower survival risk). Due to the limited availability of TNBC datasets with survival outcomes, we used nominal p-values instead of FDR to detect significant genes. The complete gene list of survival meta-analysis and pathway list from GSEA analysis can be found in Supplementary Table S3. The unbalanced association pattern observed in the pCR meta-analysis was also evident in the survival meta-analysis, with 888 genes having a hazard ratio less than 1 (upregulated gene expression associated with lower survival risk) and only 91 genes having a hazard ratio greater than 1 out of 979 genes passing the nominal p-value 0.05. Figure 3 displays the top 10 pathways with positive NES (associated with higher survival risk) and the top 10 pathways with negative NES. Our pathway analysis suggests that RNA splicing and immune-related pathways are associated with favorable survival outcomes, while cell junction and metabolism-related pathways are associated with unfavorable survival outcomes.

3.3. Integration of pCR and DFS Meta-Analysis Results

The pCR meta-analysis and DFS meta-analysis show substantial overlap in cell cycle and immune pathways; these findings are expected as pCR was shown to be a surrogate for survival outcomes for TNBC (Tables S2 and S3) [21,42]. However, the discrepancy between the two analysis results provides the opportunity to discover new biomarkers and therapeutic targets for TNBC treatments. We further dissected the significant genes from pCR and DFS meta-analysis results into four quadrants based on their association with pCR and DFS. The criteria to select genes for each quadrant were either FDR < 0.05 in pCR meta-analysis or p-value < 0.05 in DFS meta-analysis. Since there could potentially be a large number of genes selected in quadrant one, we used a more stringent criterion requiring both FDR < 0.05 in pCR meta-analysis and p-value < 0.05 in survival meta-analysis to select genes for quadrant I. The complete gene lists for each quadrant can be found in Supplementary Table S4. We discuss the downstream analyses, including functional enrichment analysis and PPI network analysis, in the following subsections.

3.4. Quadrant I: Genes Associated with pCR and Lower DFS Risk

The number of significant genes in quadrant I is dominant compared to the other three quadrants (Figure 4). These genes are associated with favorable clinical outcomes of TNBC, including pCR and lower DFS risk. The functional enrichment analysis results (Supplementary Table S5 and Figure 5A) show that immune and mRNA preprocessing-related pathways are enriched in the top-ranking list, in addition to numerous cell cycle and DNA repair pathways which is concordant with the GSEA analysis results for pCR and DFS meta-analysis. Figure 5B displays the PPI network using the Quadrant I significant genes as the seeds. SUMO2, CUL3, CAND1, GRB2, and HDAC4 are the top hub genes with the most interactions with other genes.

3.5. Quadrant II: Genes Associated with RD and Lower DFS Risk

The genes in quadrant II are relatively sparse and associated with favorable long-term clinical outcomes of TNBC, including lower DFS risk, despite their association with unfavorable short-term clinical outcomes (RD). These tumors likely represent slower growing chemo-resistant tumors that have favorable prognosis. Only a few pathways passed the FDR 0.25 threshold and are related to the extracellular matrix and elastic fiber (Supplementary Table S5 and Figure 6A). CCND1, ZBTB16, FLNC, and PPP2R2B are displayed as major hub genes in PPI network (Figure 6B).

3.6. Quadrant III: Genes Associated with RD and Higher DFS Risk

Quadrant III contains the second-largest number of genes among the four quadrants and is associated with unfavorable clinical outcomes, including RD and higher DFS risk. Functional enrichment analysis shows that the G-protein-coupled receptor (GPCR) pathways are the major significant functional gene sets (Supplementary Table S5 and Figure 7A). The hub genes include H4C9, H3C12, CFTR, SNCA, and PRKN (Figure 7B).

3.7. Quadrant IV: Genes Associated with pCR and Higher DFS Risk

Quadrant IV has the least number of genes among the four quadrants, which are associated with pCR and high survival risk. 3q29 copy number variation syndrome gene set is significantly enriched with the quadrant IV genes (Supplementary Table S5 and Figure 8A). RAD21, SRPK1, H3C3, and ISG15 are among the most connected genes within the quadrant IV PPI network (Figure 8B).

4. Discussion

TNBC is still being treated as a single disease, and NAC is the common treatment strategy for TNBC, which involves administering chemotherapy before surgery to shrink the tumor and, potentially, improve surgical outcomes. However, the heterogeneous responses of TNBC patients for NAC indicate the need for precision medicine, including the identification of biomarkers to predict NAC response and inform personalized treatment strategies.
The molecular mechanisms underlying NAC response in TNBC are complex and involve interactions between the tumor and the host immune system, as well as between the tumor and the chemotherapy drugs. In the past, numerous studies have been conducted to discover the NAC response mechanisms and biomarkers, but the inconsistent results and inconclusive prediction performances show that it is still an open research area [23,24,25,26,27,28,29,30,31]. The large-scale gene expression meta-analyses we conducted not only increase the statistical power to detect more signals associated with NAC response and survival outcome, also provide an innovative way to integrate both NAC response and survival outcome to gain more biological insight into potential NAC response mechanisms.
The immune and proliferation pathways are significantly enriched at the top pathway list in pCR meta-analysis, and immune pathways are also strongly associated with lower DFS risk in the survival meta-analysis. Those results verified the outcomes from previous studies about the association of immune and proliferation signatures with pCR and survival outcomes of TNBC patients [23,30,43,44,45]. RNA splicing-related pathways are also among the top pathways associated with both pCR and lower DFS risk, which was shown the important role of controlling the proliferation of TNBC cells [46].
SLAMF7 and CCND1 are the most significantly associated genes with pCR (HR 1.419, adjusted p-value 2.95 × 10−10) and RD (HR, 0.770, adjusted p-value 7.85 × 10−6), respectively, in the meta-analysis. SLAMF7 encodes surface antigen CD319, which is a stable marker of both normal and malignant plasma cells in multiple myeloma [47]. Besides its expression on myeloma cells, SLAMF7 is also present in immune cells, particularly in natural killer (NK) cells, where it serves as a costimulatory receptor to stimulate innate and adaptive immunity. SLAMF7 was included in the immune gene signatures to predict TNBC NAC response [29,48], and it was also used as a target by a nanoconjugate system to activate an immune response for breast cancer [49]. CCND1 is associated with RD and plays an important role in cell cycle regulation. It encodes a protein that forms a complex with cyclin-dependent kinase 4 (CDK4) or CDK6, promoting the progression of cells through the G1 phase of the cell cycle. Aberrant expression of CCND1 has been found in a variety of human cancers, including breast cancer [50,51,52]. In breast cancer, CCND1 is frequently overexpressed due to gene amplification or other mechanisms, and its overexpression has been associated with the development and progression of the disease [53,54,55,56]. CCND1 is associated with luminal breast cancers and amplifications have been found in chemoresistant post-treatment ER+ tumors [57].
A prominent result from NAC response and DFS meta-analysis is the extremely unbalanced significant upregulated genes associated with pCR and lower survival risk, which could be caused by the sheer number of differentially expressed transcripts associated with immune cells and large number or transcript involved in cell cycle progression. This phenomenon partially explains the challenge of achieving the satisfying specificity for TNBC NAC signatures [26,27,29] since the number and effect size of the positively associated genes overwhelms the negatively associated genes.
To maximally utilize the information from NAC response and DFS meta-analysis, we further divided the gene association results into four quadrants. The explanation of quadrant I and III is straightforward. Genes in quadrant I are associated with pCR and lower survival risk, and the prevalence of immune and cell cycle genes associated with chemotherapy. Genes in quadrant III are associated with RD and higher survival risk, and G-protein-coupled receptors (GPCRs) pathways stand out in enrichment analysis, which is not at the top-ranking list in NAC response or survival analysis alone.
Quadrant II genes are associated with RD and lower survival risk. Those genes could be linked to a group of TNBC patients who are insensitive to NAC but with less aggressive growth patterns of tumors. Luminal androgen receptor (LAR) subtype TNBC patients defined by our TNBC type are less responsive to NAC compared to other TNBC subtypes [58,59], but breast cancer patients with positive androgen receptor (AR) expression have favorable long-term clinical outcomes, including overall survival and recurrence [60,61]. In our recent multi-omics analysis of TNBC, we showed that LAR subtype has a genetic dependency on CCND1 [9], which happened to be the hub gene within quadrant II.
It seems to be challenging to explain the genes in quadrant IV. However, a small group of patients with micrometastases in sentinel lymph modes could have a higher risk of recurrence even with the pCR outcome in NAC [62,63]. 3q29 copy number variation syndrome pathway is far more significantly enriched than other pathways in quadrant IV genes. Previously, 3q29 region was shown as an amplified region in basal-like breast cancer with BRCA1 mutation [64], and it was also identified as a significantly altered region in the metastatic breast cancer cohort [65].
With this large-scale gene expression meta-analysis, several conclusions are summarized as the following. First, different and appropriate clinical endpoints need to be considered for NAC response. Survival outcomes were shown to be a more effective endpoint to evaluate the combination of durvalumab, which is an immune-checkpoint inhibitor, with NAC for early-stage TNBC in phase II clinical trial [66]. As we demonstrated in the four-quadrants analysis of NAC response and survival outcomes, integration of both short-term and long-term clinical outcomes is more meaningful in guiding clinical decisions. Residual cancer burden (RCB), which is a continuous index to measure treatment responses based on primary tumor and nodal metastases, is considered the better endpoint to replace pCR for the neoadjuvant treatment responses in clinical trials [67,68].
Second, the multi-modal prediction models need to be implemented for NAC prediction. The most available prediction models were built based on gene expression data. Recently, the prediction model using multi-modal data, such as DNA mutation, gene expression, copy number, digital pathology, and imaging, etc., achieved better prediction performance for multiple cancer prognosis and treatment prediction, including breast cancer NAC [69,70,71,72]. The research community needs to make an effort to collect multi-modal data for TNBC clinical trials. The unbalanced signals issue from gene expression alone could be solved using multi-modal data of TNBCs.
Third, racial disparities need to be included during the biomarker discovery and NAC prediction model building for TNBC. The incidence of TNBC compared to other types of breast cancer has been estimated to be approximately 28% among African American women and 20% among Hispanic women compared to 12% for non-Hispanic white patients [73,74]. Compared to Whites, African Americans and Hispanics were more likely to present with regional or metastatic disease and have comorbidity among breast cancer patients. Based on the National Cancer Data Base, among TNBC patients who received NAC, black women had significantly lower rates of pCR compared with white women [75]. Current evidence from multiple studies suggests that both socioeconomic diversity and tumor biology contribute to the racial disparities in TNBC [76]. There is an urgent need to include socioeconomic status and socioeconomic-genomic interaction effects in prediction models.

5. Conclusions

In conclusion, TNBC is a heterogeneous disease that requires precision medicine approaches, including the identification of biomarkers to predict NAC response and inform personalized treatment strategies. The molecular mechanisms underlying NAC response in TNBC are complex, and large-scale gene expression meta-analyses provide a way to integrate both NAC response and survival outcome to gain more biological insight. The results of this study highlight the importance of considering appropriate clinical endpoints and multi-modal prediction models for NAC prediction, as well as racial disparities during biomarker discovery and NAC prediction model building for TNBC. Overall, this study contributes to the ongoing efforts to improve the management of TNBC and personalize treatment strategies for patients with this disease.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15082194/s1, Table S1: RNAseq normalized data matrix and clinical information for the VUMC cohort; Table S2: Meta-analysis and GSEA results for pCR; Table S3: Meta-analysis and GSEA results for DSF; Table S4: Significant genes selected in each quadrant; Table S5: Functional enrichment analysis results for each quadrant.

Author Contributions

X.S.C. and B.D.L. conceived and designed the study. W.Z. led the data analysis. E.L. contributed to the bioinformatics and statistical analysis. L.W. contributed to the interpretation of the results. X.S.C. and B.D.L. drafted the paper. X.S.C. and B.D.L. supervised the study. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Cancer Institute grants R01CA200987, P50CA098131, P30CA240139, and the Department of Defense Breast Cancer Research Program BC201286.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

RNAseq normalized data matrix for the VUMC cohort provided in supplemental Table S1.

Conflicts of Interest

All authors declare that they have no competing interests.

References

  1. Dent, R.; Trudeau, M.; Pritchard, K.I.; Hanna, W.M.; Kahn, H.K.; Sawka, C.A.; Lickley, L.A.; Rawlinson, E.; Sun, P.; Narod, S.A. Triple-negative breast cancer: Clinical features and patterns of recurrence. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2007, 13, 4429–4434. [Google Scholar] [CrossRef] [Green Version]
  2. Bianchini, G.; Balko, J.M.; Mayer, I.A.; Sanders, M.E.; Gianni, L. Triple-negative breast cancer: Challenges and opportunities of a heterogeneous disease. Nat. Rev. Clin. Oncol. 2016, 13, 674–690. [Google Scholar] [CrossRef]
  3. Siegel, R.L.; Miller, K.D.; Wagle, N.S.; Jemal, A. Cancer statistics, 2023. CA Cancer J. Clin. 2023, 73, 17–48. [Google Scholar] [CrossRef]
  4. Haffty, B.G.; Yang, Q.; Reiss, M.; Kearney, T.; Higgins, S.A.; Weidhaas, J.; Harris, L.; Hait, W.; Toppmeyer, D. Locoregional relapse and distant metastasis in conservatively managed triple negative early-stage breast cancer. J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol. 2006, 24, 5652–5657. [Google Scholar] [CrossRef]
  5. Lehmann, B.D.; Bauer, J.A.; Chen, X.; Sanders, M.E.; Chakravarthy, A.B.; Shyr, Y.; Pietenpol, J.A. Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. J. Clin. Investig. 2011, 121, 2750–2767. [Google Scholar] [CrossRef] [Green Version]
  6. Lehmann, B.D.; Jovanovic, B.; Chen, X.; Estrada, M.V.; Johnson, K.N.; Shyr, Y.; Moses, H.L.; Sanders, M.E.; Pietenpol, J.A. Refinement of Triple-Negative Breast Cancer Molecular Subtypes: Implications for Neoadjuvant Chemotherapy Selection. PLoS ONE 2016, 11, e0157368. [Google Scholar] [CrossRef]
  7. Masuda, H.; Baggerly, K.A.; Wang, Y.; Zhang, Y.; Gonzalez-Angulo, A.M.; Meric-Bernstam, F.; Valero, V.; Lehmann, B.D.; Pietenpol, J.A.; Hortobagyi, G.N.; et al. Differential response to neoadjuvant chemotherapy among 7 triple-negative breast cancer molecular subtypes. Clin. Cancer Res. 2013, 19, 5533–5540. [Google Scholar] [CrossRef] [Green Version]
  8. Chen, X.; Li, J.; Gray, W.H.; Lehmann, B.D.; Bauer, J.A.; Shyr, Y.; Pietenpol, J.A. TNBCtype: A Subtyping Tool for Triple-Negative Breast Cancer. Cancer Inform. 2012, 11, 147–156. [Google Scholar] [CrossRef] [PubMed]
  9. Lehmann, B.D.; Colaprico, A.; Silva, T.C.; Chen, J.; An, H.; Ban, Y.; Huang, H.; Wang, L.; James, J.L.; Balko, J.M.; et al. Multi-omics analysis identifies therapeutic vulnerabilities in triple-negative breast cancer subtypes. Nat. Commun. 2021, 12, 6276. [Google Scholar] [CrossRef] [PubMed]
  10. Burstein, M.D.; Tsimelzon, A.; Poage, G.M.; Covington, K.R.; Contreras, A.; Fuqua, S.A.; Savage, M.I.; Osborne, C.K.; Hilsenbeck, S.G.; Chang, J.C.; et al. Comprehensive genomic analysis identifies novel subtypes and targets of triple-negative breast cancer. Clin. Cancer Res. 2015, 21, 1688–1698. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Bareche, Y.; Venet, D.; Ignatiadis, M.; Aftimos, P.; Piccart, M.; Rothe, F.; Sotiriou, C. Unravelling triple-negative breast cancer molecular heterogeneity using an integrative multiomic analysis. Ann. Oncol. 2018, 29, 895–902. [Google Scholar] [CrossRef]
  12. Bianchini, G.; De Angelis, C.; Licata, L.; Gianni, L. Treatment landscape of triple-negative breast cancer—Expanded options, evolving needs. Nat. Rev. Clin. Oncol. 2022, 19, 91–113. [Google Scholar] [CrossRef]
  13. Schmid, P.; Adams, S.; Rugo, H.S.; Schneeweiss, A.; Barrios, C.H.; Iwata, H.; Dieras, V.; Hegg, R.; Im, S.A.; Shaw Wright, G.; et al. Atezolizumab and Nab-Paclitaxel in Advanced Triple-Negative Breast Cancer. N. Engl. J. Med. 2018, 379, 2108–2121. [Google Scholar] [CrossRef]
  14. Cortes, J.; Cescon, D.W.; Rugo, H.S.; Nowecki, Z.; Im, S.A.; Yusof, M.M.; Gallardo, C.; Lipatov, O.; Barrios, C.H.; Holgado, E.; et al. 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. Lancet 2020, 396, 1817–1828. [Google Scholar] [CrossRef] [PubMed]
  15. Robson, M.; Im, S.A.; Senkus, E.; Xu, B.; Domchek, S.M.; Masuda, N.; Delaloge, S.; Li, W.; Tung, N.; Armstrong, A.; et al. Olaparib for Metastatic Breast Cancer in Patients with a Germline BRCA Mutation. N. Engl. J. Med. 2017, 377, 523–533. [Google Scholar] [CrossRef]
  16. Litton, J.K.; Rugo, H.S.; Ettl, J.; Hurvitz, S.A.; Goncalves, A.; Lee, K.H.; Fehrenbacher, L.; Yerushalmi, R.; Mina, L.A.; Martin, M.; et al. Talazoparib in Patients with Advanced Breast Cancer and a Germline BRCA Mutation. N. Engl. J. Med. 2018, 379, 753–763. [Google Scholar] [CrossRef] [PubMed]
  17. Modi, S.; Saura, C.; Yamashita, T.; Park, Y.H.; Kim, S.B.; Tamura, K.; Andre, F.; Iwata, H.; Ito, Y.; Tsurutani, J.; et al. Trastuzumab Deruxtecan in Previously Treated HER2-Positive Breast Cancer. N. Engl. J. Med. 2020, 382, 610–621. [Google Scholar] [CrossRef]
  18. Modi, S.; Park, H.; Murthy, R.K.; Iwata, H.; Tamura, K.; Tsurutani, J.; Moreno-Aspitia, A.; Doi, T.; Sagara, Y.; Redfern, C.; et al. Antitumor Activity and Safety of Trastuzumab Deruxtecan in Patients With HER2-Low-Expressing Advanced Breast Cancer: Results From a Phase Ib Study. J. Clin. Oncol. 2020, 38, 1887–1896. [Google Scholar] [CrossRef]
  19. Mieog, J.S.; van der Hage, J.A.; van de Velde, C.J. Preoperative chemotherapy for women with operable breast cancer. Cochrane Database Syst. Rev. 2007, 2007, CD005002. [Google Scholar] [CrossRef] [PubMed]
  20. Carey, L.A.; Dees, E.C.; Sawyer, L.; Gatti, L.; Moore, D.T.; Collichio, F.; Ollila, D.W.; Sartor, C.I.; Graham, M.L.; Perou, C.M. The triple negative paradox: Primary tumor chemosensitivity of breast cancer subtypes. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2007, 13, 2329–2334. [Google Scholar] [CrossRef] [Green Version]
  21. Liedtke, C.; Mazouni, C.; Hess, K.R.; Andre, F.; Tordai, A.; Mejia, J.A.; Symmans, W.F.; Gonzalez-Angulo, A.M.; Hennessy, B.; Green, M.; et al. Response to neoadjuvant therapy and long-term survival in patients with triple-negative breast cancer. J. Clin. Oncol. 2008, 26, 1275–1281. [Google Scholar] [CrossRef]
  22. Straver, M.E.; Glas, A.M.; Hannemann, J.; Wesseling, J.; van de Vijver, M.J.; Rutgers, E.J.T.; Peeters, M.J.T.F.D.V.; van Tinteren, H.; van’t Veer, L.J.; Rodenhuis, S. The 70-gene signature as a response predictor for neoadjuvant chemotherapy in breast cancer. Breast Cancer Res. Treat. 2010, 119, 551–558. [Google Scholar] [CrossRef] [Green Version]
  23. Filho, O.M.; Stover, D.G.; Asad, S.; Ansell, P.J.; Watson, M.; Loibl, S.; Geyer, C.E., Jr.; Bae, J.; Collier, K.; Cherian, M.; et al. Association of Immunophenotype With Pathologic Complete Response to Neoadjuvant Chemotherapy for Triple-Negative Breast Cancer: A Secondary Analysis of the BrighTNess Phase 3 Randomized Clinical Trial. JAMA Oncol. 2021, 7, 603–608. [Google Scholar] [CrossRef] [PubMed]
  24. Orozco, J.I.J.; Grumley, J.G.; Matsuba, C.; Manughian-Peter, A.O.; Chang, S.C.; Chang, G.; Gago, F.E.; Salomon, M.P.; Marzese, D.M. Clinical Implications of Transcriptomic Changes After Neoadjuvant Chemotherapy in Patients with Triple-Negative Breast Cancer. Ann. Surg. Oncol. 2019, 26, 3185–3193. [Google Scholar] [CrossRef]
  25. Zhao, Y.; Schaafsma, E.; Cheng, C. Gene signature-based prediction of triple-negative breast cancer patient response to Neoadjuvant chemotherapy. Cancer Med. 2020, 9, 6281–6295. [Google Scholar] [CrossRef] [PubMed]
  26. Perez-Pena, J.; Tibor Fekete, J.; Paez, R.; Baliu-Pique, M.; Garcia-Saenz, J.A.; Garcia-Barberan, V.; Manzano, A.; Perez-Segura, P.; Esparis-Ogando, A.; Pandiella, A.; et al. A Transcriptomic Immunologic Signature Predicts Favorable Outcome in Neoadjuvant Chemotherapy Treated Triple Negative Breast Tumors. Front. Immunol. 2019, 10, 2802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Fournier, M.V.; Goodwin, E.C.; Chen, J.; Obenauer, J.C.; Tannenbaum, S.H.; Brufsky, A.M. A Predictor of Pathological Complete Response to Neoadjuvant Chemotherapy Stratifies Triple Negative Breast Cancer Patients with High Risk of Recurrence. Sci. Rep. 2019, 9, 14863. [Google Scholar] [CrossRef] [Green Version]
  28. Kim, E.K.; Park, A.K.; Ko, E.; Park, W.Y.; Lee, K.M.; Noh, D.Y.; Han, W. Risk stratification of triple-negative breast cancer with core gene signatures associated with chemoresponse and prognosis. Breast Cancer Res. Treat. 2019, 178, 185–197. [Google Scholar] [CrossRef]
  29. Vishnubalaji, R.; Abdel-Razeq, H.; Gehani, S.; Albagha, O.M.E.; Alajez, N.M. Identification of a Gene Panel Predictive of Triple-Negative Breast Cancer Response to Neoadjuvant Chemotherapy Employing Transcriptomic and Functional Validation. Int. J. Mol. Sci. 2022, 23, 10901. [Google Scholar] [CrossRef]
  30. Prat, A.; Lluch, A.; Albanell, J.; Barry, W.T.; Fan, C.; Chacon, J.I.; Parker, J.S.; Calvo, L.; Plazaola, A.; Arcusa, A.; et al. Predicting response and survival in chemotherapy-treated triple-negative breast cancer. Br. J. Cancer 2014, 111, 1532–1541. [Google Scholar] [CrossRef] [Green Version]
  31. Sinn, B.V.; Loibl, S.; Hanusch, C.A.; Zahm, D.M.; Sinn, H.P.; Untch, M.; Weber, K.; Karn, T.; Becker, C.; Marme, F.; et al. Immune-related Gene Expression Predicts Response to Neoadjuvant Chemotherapy but not Additional Benefit from PD-L1 Inhibition in Women with Early Triple-negative Breast Cancer. Clin. Cancer Res. 2021, 27, 2584–2591. [Google Scholar] [CrossRef]
  32. Jovanovic, B.; Mayer, I.A.; Mayer, E.L.; Abramson, V.G.; Bardia, A.; Sanders, M.E.; Kuba, M.G.; Estrada, M.V.; Beeler, J.S.; Shaver, T.M.; et al. A Randomized Phase II Neoadjuvant Study of Cisplatin, Paclitaxel With or Without Everolimus in Patients with Stage II/III Triple-Negative Breast Cancer (TNBC): Responses and Long-term Outcome Correlated with Increased Frequency of DNA Damage Response Gene Mutations, TNBC Subtype, AR Status, and Ki67. Clin. Cancer Res. 2017, 23, 4035–4045. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. McCall, M.N.; Bolstad, B.M.; Irizarry, R.A. Frozen robust multiarray analysis (fRMA). Biostatistics 2010, 11, 242–253. [Google Scholar] [CrossRef] [Green Version]
  34. Karn, T.; Metzler, D.; Ruckhaberle, E.; Hanker, L.; Gatje, R.; Solbach, C.; Ahr, A.; Schmidt, M.; Holtrich, U.; Kaufmann, M.; et al. Data-driven derivation of cutoffs from a pool of 3,030 Affymetrix arrays to stratify distinct clinical types of breast cancer. Breast Cancer Res. Treat. 2010, 120, 567–579. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Riley, R.D.; Higgins, J.P.; Deeks, J.J. Interpretation of random effects meta-analyses. BMJ 2011, 342, d549. [Google Scholar] [CrossRef] [Green Version]
  36. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [Green Version]
  37. Boyle, E.I.; Weng, S.; Gollub, J.; Jin, H.; Botstein, D.; Cherry, J.M.; Sherlock, G. GO::TermFinder--open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics 2004, 20, 3710–3715. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Wu, T.; Hu, E.; Xu, S.; Chen, M.; Guo, P.; Dai, Z.; Feng, T.; Zhou, L.; Tang, W.; Zhan, L.; et al. clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation 2021, 2, 100141. [Google Scholar] [CrossRef] [PubMed]
  39. Liberzon, A.; Birger, C.; Thorvaldsdottir, H.; Ghandi, M.; Mesirov, J.P.; Tamayo, P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 2015, 1, 417–425. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Zhou, G.; Soufan, O.; Ewald, J.; Hancock, R.E.W.; Basu, N.; Xia, J. NetworkAnalyst 3.0: A visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. 2019, 47, W234–W241. [Google Scholar] [CrossRef] [Green Version]
  41. Orchard, S.; Kerrien, S.; Abbani, S.; Aranda, B.; Bhate, J.; Bidwell, S.; Bridge, A.; Briganti, L.; Brinkman, F.S.; Cesareni, G.; et al. Protein interaction data curation: The International Molecular Exchange (IMEx) consortium. Nat. Methods 2012, 9, 345–350. [Google Scholar] [CrossRef]
  42. von Minckwitz, G.; Untch, M.; Blohmer, J.U.; Costa, S.D.; Eidtmann, H.; Fasching, P.A.; Gerber, B.; Eiermann, W.; Hilfrich, J.; Huober, J.; et al. Definition and impact of pathologic complete response on prognosis after neoadjuvant chemotherapy in various intrinsic breast cancer subtypes. J. Clin. Oncol. 2012, 30, 1796–1804. [Google Scholar] [CrossRef] [Green Version]
  43. Denkert, C.; von Minckwitz, G.; Brase, J.C.; Sinn, B.V.; Gade, S.; Kronenwett, R.; Pfitzner, B.M.; Salat, C.; Loi, S.; Schmitt, W.D.; et al. Tumor-infiltrating lymphocytes and response to neoadjuvant chemotherapy with or without carboplatin in human epidermal growth factor receptor 2-positive and triple-negative primary breast cancers. J. Clin. Oncol. 2015, 33, 983–991. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, Y.; Asad, S.; Weber, Z.; Tallman, D.; Nock, W.; Wyse, M.; Bey, J.F.; Dean, K.L.; Adams, E.J.; Stockard, S.; et al. Genomic features of rapid versus late relapse in triple negative breast cancer. BMC Cancer 2021, 21, 568. [Google Scholar] [CrossRef] [PubMed]
  45. Shepherd, J.H.; Ballman, K.; Polley, M.C.; Campbell, J.D.; Fan, C.; Selitsky, S.; Fernandez-Martinez, A.; Parker, J.S.; Hoadley, K.A.; Hu, Z.; et al. CALGB 40603 (Alliance): Long-Term Outcomes and Genomic Correlates of Response and Survival After Neoadjuvant Chemotherapy With or Without Carboplatin and Bevacizumab in Triple-Negative Breast Cancer. J. Clin. Oncol. 2022, 40, 1323–1334. [Google Scholar] [CrossRef] [PubMed]
  46. Koedoot, E.; van Steijn, E.; Vermeer, M.; Gonzalez-Prieto, R.; Vertegaal, A.C.O.; Martens, J.W.M.; Le Devedec, S.E.; van de Water, B. Splicing factors control triple-negative breast cancer cell mitosis through SUN2 interaction and sororin intron retention. J. Exp. Clin. Cancer Res. 2021, 40, 82. [Google Scholar] [CrossRef] [PubMed]
  47. Frigyesi, I.; Adolfsson, J.; Ali, M.; Christophersen, M.K.; Johnsson, E.; Turesson, I.; Gullberg, U.; Hansson, M.; Nilsson, B. Robust isolation of malignant plasma cells in multiple myeloma. Blood 2014, 123, 1336–1340. [Google Scholar] [CrossRef] [Green Version]
  48. Vishnubalaji, R.; Alajez, N.M. Transcriptional landscape associated with TNBC resistance to neoadjuvant chemotherapy revealed by single-cell RNA-seq. Mol. Ther. Oncolytics 2021, 23, 151–162. [Google Scholar] [CrossRef] [PubMed]
  49. Lu, Y.; Huntoon, K.; Lee, D.; Wang, Y.; Ha, J.; Qie, Y.; Li, X.; Schrank, B.R.; Dong, S.; Gallup, T.D.; et al. Immunological conversion of solid tumours using a bispecific nanobioconjugate for cancer immunotherapy. Nat. Nanotechnol. 2022, 17, 1332–1341. [Google Scholar] [CrossRef]
  50. Roy, P.G.; Thompson, A.M. Cyclin D1 and breast cancer. Breast 2006, 15, 718–727. [Google Scholar] [CrossRef]
  51. Montalto, F.I.; De Amicis, F. Cyclin D1 in Cancer: A Molecular Connection for Cell Cycle Control, Adhesion and Invasion in Tumor and Stroma. Cells 2020, 9, 2648. [Google Scholar] [CrossRef] [PubMed]
  52. Chen, Y.; Huang, Y.; Gao, X.; Li, Y.; Lin, J.; Chen, L.; Chang, L.; Chen, G.; Guan, Y.; Pan, L.K.; et al. CCND1 Amplification Contributes to Immunosuppression and Is Associated With a Poor Prognosis to Immune Checkpoint Inhibitors in Solid Tumors. Front. Immunol. 2020, 11, 1620. [Google Scholar] [CrossRef] [PubMed]
  53. Lundberg, A.; Lindstrom, L.S.; Li, J.; Harrell, J.C.; Darai-Ramqvist, E.; Sifakis, E.G.; Foukakis, T.; Perou, C.M.; Czene, K.; Bergh, J.; et al. The long-term prognostic and predictive capacity of cyclin D1 gene amplification in 2305 breast tumours. Breast Cancer Res. 2019, 21, 34. [Google Scholar] [CrossRef] [Green Version]
  54. Reis-Filho, J.S.; Savage, K.; Lambros, M.B.; James, M.; Steele, D.; Jones, R.L.; Dowsett, M. Cyclin D1 protein overexpression and CCND1 amplification in breast carcinomas: An immunohistochemical and chromogenic in situ hybridisation analysis. Mod. Pathol. 2006, 19, 999–1009. [Google Scholar] [CrossRef] [Green Version]
  55. Jeffreys, S.A.; Becker, T.M.; Khan, S.; Soon, P.; Neubauer, H.; de Souza, P.; Powter, B. Prognostic and Predictive Value of CCND1/Cyclin D1 Amplification in Breast Cancer With a Focus on Postmenopausal Patients: A Systematic Review and Meta-Analysis. Front. Endocrinol. 2022, 13, 895729. [Google Scholar] [CrossRef]
  56. Valla, M.; Klaestad, E.; Ytterhus, B.; Bofin, A.M. CCND1 Amplification in Breast Cancer -associations With Proliferation, Histopathological Grade, Molecular Subtype and Prognosis. J. Mammary Gland. Biol. Neoplasia 2022, 27, 67–77. [Google Scholar] [CrossRef] [PubMed]
  57. Hoogstraat, M.; Lips, E.H.; Mayayo-Peralta, I.; Mulder, L.; Kristel, P.; van der Heijden, I.; Annunziato, S.; van Seijen, M.; Nederlof, P.M.; Sonke, G.S.; et al. Comprehensive characterization of pre- and post-treatment samples of breast cancer reveal potential mechanisms of chemotherapy resistance. NPJ Breast Cancer 2022, 8, 60. [Google Scholar] [CrossRef]
  58. Echavarria, I.; Lopez-Tarruella, S.; Picornell, A.; Garcia-Saenz, J.A.; Jerez, Y.; Hoadley, K.; Gomez, H.L.; Moreno, F.; Monte-Millan, M.D.; Marquez-Rodas, I.; et al. Pathological Response in a Triple-Negative Breast Cancer Cohort Treated with Neoadjuvant Carboplatin and Docetaxel According to Lehmann’s Refined Classification. Clin. Cancer Res. 2018, 24, 1845–1852. [Google Scholar] [CrossRef] [Green Version]
  59. Santonja, A.; Sanchez-Munoz, A.; Lluch, A.; Chica-Parrado, M.R.; Albanell, J.; Chacon, J.I.; Antolin, S.; Jerez, J.M.; de la Haba, J.; de Luque, V.; et al. Triple negative breast cancer subtypes and pathologic complete response rate to neoadjuvant chemotherapy. Oncotarget 2018, 9, 26406–26416. [Google Scholar] [CrossRef] [Green Version]
  60. Qu, Q.; Mao, Y.; Fei, X.C.; Shen, K.W. The impact of androgen receptor expression on breast cancer survival: A retrospective study and meta-analysis. PLoS ONE 2013, 8, e82650. [Google Scholar] [CrossRef] [Green Version]
  61. Vera-Badillo, F.E.; Templeton, A.J.; de Gouveia, P.; Diaz-Padilla, I.; Bedard, P.L.; Al-Mubarak, M.; Seruga, B.; Tannock, I.F.; Ocana, A.; Amir, E. Androgen receptor expression and outcomes in early breast cancer: A systematic review and meta-analysis. J. Natl. Cancer Inst. 2014, 106, djt319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Dutta, S.W.; Volaric, A.; Morgan, J.T.; Chinn, Z.; Atkins, K.A.; Janowski, E.M. Pathologic Evaluation and Prognostic Implications of Nodal Micrometastases in Breast Cancer. Semin. Radiat. Oncol. 2019, 29, 102–110. [Google Scholar] [CrossRef]
  63. Houvenaeghel, G.; Sabatier, R.; Reyal, F.; Classe, J.M.; Giard, S.; Charitansky, H.; Rouzier, R.; Faure, C.; Garbay, J.R.; Darai, E.; et al. Axillary lymph node micrometastases decrease triple-negative early breast cancer survival. Br. J. Cancer 2016, 115, 1024–1031. [Google Scholar] [CrossRef] [Green Version]
  64. Prat, A.; Cruz, C.; Hoadley, K.A.; Diez, O.; Perou, C.M.; Balmana, J. Molecular features of the basal-like breast cancer subtype based on BRCA1 mutation status. Breast Cancer Res. Treat. 2014, 147, 185–191. [Google Scholar] [CrossRef] [Green Version]
  65. Paul, M.R.; Pan, T.C.; Pant, D.K.; Shih, N.N.; Chen, Y.; Harvey, K.L.; Solomon, A.; Lieberman, D.; Morrissette, J.J.; Soucier-Ernst, D.; et al. Genomic landscape of metastatic breast cancer identifies preferentially dysregulated pathways and targets. J. Clin. Investig. 2020, 130, 4252–4265. [Google Scholar] [CrossRef]
  66. Loibl, S.; Untch, M.; Burchardi, N.; Huober, J.; Sinn, B.V.; Blohmer, J.U.; Grischke, E.M.; Furlanetto, J.; Tesch, H.; Hanusch, C.; et al. A randomised phase II study investigating durvalumab in addition to an anthracycline taxane-based neoadjuvant therapy in early triple-negative breast cancer: Clinical results and biomarker analysis of GeparNuevo study. Ann. Oncol. 2019, 30, 1279–1288. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Agostinetto, E.; Gligorov, J.; Piccart, M. Systemic therapy for early-stage breast cancer: Learning from the past to build the future. Nat. Rev. Clin. Oncol. 2022, 19, 763–774. [Google Scholar] [CrossRef]
  68. Symmans, W.F.; Wei, C.; Gould, R.; Yu, X.; Zhang, Y.; Liu, M.; Walls, A.; Bousamra, A.; Ramineni, M.; Sinn, B.; et al. Long-Term Prognostic Risk After Neoadjuvant Chemotherapy Associated With Residual Cancer Burden and Breast Cancer Subtype. J. Clin. Oncol. 2017, 35, 1049–1060. [Google Scholar] [CrossRef] [Green Version]
  69. Sammut, S.J.; Crispin-Ortuzar, M.; Chin, S.F.; Provenzano, E.; Bardwell, H.A.; Ma, W.; Cope, W.; Dariush, A.; Dawson, S.J.; Abraham, J.E.; et al. Multi-omic machine learning predictor of breast cancer therapy response. Nature 2022, 601, 623–629. [Google Scholar] [CrossRef] [PubMed]
  70. Vanguri, R.S.; Luo, J.; Aukerman, A.T.; Egger, J.V.; Fong, C.J.; Horvat, N.; Pagano, A.; Araujo-Filho, J.A.B.; Geneslaw, L.; Rizvi, H.; et al. Multimodal integration of radiology, pathology and genomics for prediction of response to PD-(L)1 blockade in patients with non-small cell lung cancer. Nat. Cancer 2022, 3, 1151–1164. [Google Scholar] [CrossRef]
  71. Chen, R.J.; Lu, M.Y.; Williamson, D.F.K.; Chen, T.Y.; Lipkova, J.; Noor, Z.; Shaban, M.; Shady, M.; Williams, M.; Joo, B.; et al. Pan-cancer integrative histology-genomic analysis via multimodal deep learning. Cancer Cell 2022, 40, 865–878. [Google Scholar] [CrossRef] [PubMed]
  72. Boehm, K.M.; Aherne, E.A.; Ellenson, L.; Nikolovski, I.; Alghamdi, M.; Vazquez-Garcia, I.; Zamarin, D.; Long Roche, K.; Liu, Y.; Patel, D.; et al. Multimodal data integration using machine learning improves risk stratification of high-grade serous ovarian cancer. Nat. Cancer 2022, 3, 723–733. [Google Scholar] [CrossRef]
  73. Lara-Medina, F.; Perez-Sanchez, V.; Saavedra-Perez, D.; Blake-Cerda, M.; Arce, C.; Motola-Kuba, D.; Villarreal-Garza, C.; Gonzalez-Angulo, A.M.; Bargallo, E.; Aguilar, J.L.; et al. Triple-negative breast cancer in Hispanic patients: High prevalence, poor prognosis, and association with menopausal status, body mass index, and parity. Cancer 2011, 117, 3658–3669. [Google Scholar] [CrossRef]
  74. Lin, N.U.; Vanderplas, A.; Hughes, M.E.; Theriault, R.L.; Edge, S.B.; Wong, Y.N.; Blayney, D.W.; Niland, J.C.; Winer, E.P.; Weeks, J.C. Clinicopathologic features, patterns of recurrence, and survival among women with triple-negative breast cancer in the National Comprehensive Cancer Network. Cancer 2012, 118, 5463–5472. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Killelea, B.K.; Yang, V.Q.; Wang, S.Y.; Hayse, B.; Mougalian, S.; Horowitz, N.R.; Chagpar, A.B.; Pusztai, L.; Lannin, D.R. Racial Differences in the Use and Outcome of Neoadjuvant Chemotherapy for Breast Cancer: Results From the National Cancer Data Base. J. Clin. Oncol. 2015, 33, 4267–4276. [Google Scholar] [CrossRef] [PubMed]
  76. Dietze, E.C.; Sistrunk, C.; Miranda-Carboni, G.; O’Regan, R.; Seewaldt, V.L. Triple-negative breast cancer in African-American women: Disparities versus biology. Nat. Rev. Cancer 2015, 15, 248–254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The study design and workflow of the transcriptome meta-analysis of TNBC.
Figure 1. The study design and workflow of the transcriptome meta-analysis of TNBC.
Cancers 15 02194 g001
Figure 2. Top pathways of GSEA analysis for TNBC NAC response meta-analysis.
Figure 2. Top pathways of GSEA analysis for TNBC NAC response meta-analysis.
Cancers 15 02194 g002
Figure 3. Top pathways of GSEA analysis for TNBC DFS meta-analysis.
Figure 3. Top pathways of GSEA analysis for TNBC DFS meta-analysis.
Cancers 15 02194 g003
Figure 4. Quadrant plot of significant genes from integration of TNBC NAC and DFS meta-analysis.
Figure 4. Quadrant plot of significant genes from integration of TNBC NAC and DFS meta-analysis.
Cancers 15 02194 g004
Figure 5. (A) Top pathways from functional enrichment analysis of quadrant I genes (B) PPI network of quadrant I genes.
Figure 5. (A) Top pathways from functional enrichment analysis of quadrant I genes (B) PPI network of quadrant I genes.
Cancers 15 02194 g005
Figure 6. (A) Top pathways from functional enrichment analysis of quadrant II genes. (B) PPI network of quadrant II genes.
Figure 6. (A) Top pathways from functional enrichment analysis of quadrant II genes. (B) PPI network of quadrant II genes.
Cancers 15 02194 g006
Figure 7. (A) Top pathways from functional enrichment analysis of quadrant III genes. (B) PPI network of quadrant III genes.
Figure 7. (A) Top pathways from functional enrichment analysis of quadrant III genes. (B) PPI network of quadrant III genes.
Cancers 15 02194 g007
Figure 8. (A) Top pathways from functional enrichment analysis of quadrant IV genes. (B) PPI network of quadrant IV genes.
Figure 8. (A) Top pathways from functional enrichment analysis of quadrant IV genes. (B) PPI network of quadrant IV genes.
Cancers 15 02194 g008
Table 1. Patient and characteristics of gene expression datasets.
Table 1. Patient and characteristics of gene expression datasets.
DatasetsTypeNumbers of SamplepCRDSF EventMedian Follow-Up
(sd, Range), Years
Mean Age
(sd, Range)
TotalTNBCYesNoNAYesNo
GSE16446Microarray120861073318682.9 (1.5, 0.2–5.9)/
GSE20271Microarray178591346////51 (10.8, 29–74)
GSE22093Microarray103391623////46.7 (10.2, 31–67)
GSE18864Microarray8424816////49.8 (10.3, 29–68)
GSE20194Microarray278712546////50.4 (10.8, 29–75)
GSE25055Microarray3101194078137822.1 (1.7, 0.1–7.4)49.6 (10.8, 28–75)
GSE25065Microarray198632236522412.8 (1.5, 0.4–6.1)49.1 (11.2, 24–72)
GSE18728Microarray6122616/////
GSE41998Microarray279145478315////
GSE32646Microarray115261016////54.4 (12.5, 28–70)
GSE164458RNA-Seq482482236246/////
VUMCRNA-Seq145452916/8372.6 (0.8, 0.7–5.8)50.6 (10.8, 29–75)
GSE154524RNA-Seq389389210179/1152745.3 (2, 0.1–8.1)/
GSE22226RNA-Seq15027720/11163.6 (2.2, 0.5–6.7)47 (8, 33.5–63.2)
GSE192341RNA-Seq87522129210421.8 (1.2, 0.4–4.2)44.6 (11.7, 26–73)
GSE163882RNA-Seq222903852////54.1 (11.5, 23–74)
Table 2. Top 10 positive and negative associated genes with pCR outcomes.
Table 2. Top 10 positive and negative associated genes with pCR outcomes.
GeneOdds
Ratio
Estimatesez-Valuep-ValueFDRDirection *Number of Datasets
Positive associated genes
SLAMF71.4190.3500.0457.7926.62 × 10−152.95 × 10−10+++++-++++-+-+++16
GBP51.3980.3350.0437.7361.02 × 10−142.95 × 10−10+++++-+-++10
GZMB1.3560.3050.0417.4618.59 × 10−141.65 × 10−9++++-+++++++-+++16
CD2741.5900.4640.0637.4041.32 × 10−131.91 × 10−9+++++++-++10
GBP41.4740.3880.0547.1896.54 × 10−137.56 × 10−9+++++++-++10
NUP2101.5990.4700.0676.9723.14 × 10−123.02 × 10−8++++--++++++++++16
NOL72.0170.7020.1016.9473.73 × 10−123.08 × 10−8++++-+++++++++++16
MCM51.7490.5590.0826.8437.76 × 10−125.61 × 10−8+++-+++++++++++15
NKG71.3480.2980.0446.8099.84 × 10−126.32 × 10−8++++++++++--+++15
CXCL101.2440.2180.0326.7751.25 × 10−117.20 × 10−8+++++-+++++-+++15
Negative associated genes
CCND10.770−0.2610.046−5.6951.24 × 10−87.85 × 10−6-+-+-+----------16
LINC006220.635−0.4540.086−5.2601.44 × 10−74.49 × 10−5---+---+8
CTSF0.683−0.3810.074−5.1242.98 × 10−78.13 × 10−5----+---++-++--+16
SLC4A110.814−0.2060.044−4.6383.51 × 10−65.23 × 10−4----+--+--10
ZNF6970.651−0.4290.095−4.5355.77 × 10−67.57 × 10−4--------++10
AHNAK20.825−0.1930.043−4.5006.80 × 10−68.68 × 10−4-+--+--+-------+16
OPLAH0.800−0.2230.052−4.2602.05 × 10−52.04 × 10−3----++---------+16
MAPRE30.702−0.3530.086−4.0894.34 × 10−53.69 × 10−3-+--+----------+16
MGAM0.767−0.2660.066−4.0475.19 × 10−54.26 × 10−3--------------+15
MSX10.755−0.2820.071−3.9717.15 × 10−55.29 × 10−3---+-----------+16
* + repressents postive association in the single dataset; - represents negtative association in the single dataset.
Table 3. Top 10 positive and negative associated genes with DFS outcomes.
Table 3. Top 10 positive and negative associated genes with DFS outcomes.
GeneHazard
Ratio
Estimatesez-Valuep-ValueFDRDirection *Number of Datasets
Positive associated genes
MST1L1.3670.3120.0783.9856.74 × 10−58.70 × 10−2+-++4
LINC009052.0430.7150.2123.3697.56 × 10−42.55 × 10−1+++3
MYBPH1.1650.1530.0483.1631.56 × 10−33.13 × 10−1+-+-+++7
FAM20C1.3720.3160.1003.1511.63 × 10−33.16 × 10−1+++-+5
CLDN161.1450.1350.0433.1201.81 × 10−33.35 × 10−1+-++++6
SCGB2A11.0980.0930.0303.0902.00 × 10−33.46 × 10−1+++++++7
TGM51.1320.1240.0413.0452.33 × 10−33.46 × 10−1++++++6
VPS502.2020.7900.2652.9852.84 × 10−33.46 × 10−1+++-4
LINC000321.3020.2640.0902.9333.35 × 10−33.69 × 10−1+-+3
TBC1D214.3301.4660.5052.9023.70 × 10−33.77 × 10−1+-+3
Negative associated genes
PPP1R12A0.628−0.4650.107−4.3311.48 × 10−58.70 × 10−2-------7
AKAP50.627−0.4670.110−4.2362.28 × 10−58.70 × 10−2---+--6
SLAMF70.818−0.2000.049−4.0774.57 × 10−58.70 × 10−2-------7
VEZF10.640−0.4460.110−4.0754.61 × 10−58.70 × 10−2-------7
OGG10.586−0.5340.132−4.0604.90 × 10−58.70 × 10−2-----+-7
LRRK20.651−0.4290.106−4.0485.16 × 10−58.70 × 10−2-----5
RBMXL10.596−0.5170.131−3.9438.03 × 10−58.70 × 10−2---+-5
HERPUD10.670−0.4010.104−3.8371.25 × 10−49.95 × 10−2------6
PDCD1LG20.729−0.3160.083−3.7951.48 × 10−41.12 × 10−1---+---7
TXNDC50.639−0.4480.118−3.7921.49 × 10−41.12 × 10−1---3
* + repressents postive association in the single dataset; - represents negtative association in the single dataset.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, W.; Li, E.; Wang, L.; Lehmann, B.D.; Chen, X.S. Transcriptome Meta-Analysis of Triple-Negative Breast Cancer Response to Neoadjuvant Chemotherapy. Cancers 2023, 15, 2194. https://doi.org/10.3390/cancers15082194

AMA Style

Zhang W, Li E, Wang L, Lehmann BD, Chen XS. Transcriptome Meta-Analysis of Triple-Negative Breast Cancer Response to Neoadjuvant Chemotherapy. Cancers. 2023; 15(8):2194. https://doi.org/10.3390/cancers15082194

Chicago/Turabian Style

Zhang, Wei, Emma Li, Lily Wang, Brian D. Lehmann, and X. Steven Chen. 2023. "Transcriptome Meta-Analysis of Triple-Negative Breast Cancer Response to Neoadjuvant Chemotherapy" Cancers 15, no. 8: 2194. https://doi.org/10.3390/cancers15082194

APA Style

Zhang, W., Li, E., Wang, L., Lehmann, B. D., & Chen, X. S. (2023). Transcriptome Meta-Analysis of Triple-Negative Breast Cancer Response to Neoadjuvant Chemotherapy. Cancers, 15(8), 2194. https://doi.org/10.3390/cancers15082194

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop