NAD+ biosynthesis metabolism predicts prognosis and indicates immune microenvironment for breast cancer
Yuting Yang
Ze Wang2†, - 1Department of Immunology, Medical College of Qingdao University, Qingdao, Shandong, China
- 2State Key Laboratory of Proteomics, Beijing Proteome Research Center, National Center for Protein Sciences, Beijing Institute of Lifeomics, Beijing, China
- 3School of Traditional Chinese Medicine, Beijing University of Chinese Medicine, Beijing, China
- 4School of Life Sciences, Hebei University, Baoding, Hebei, China
The growing evidence implies that tumor cells need to increase NAD+ levels by upregulating NAD+ biosynthesis to satisfy their growth demand. NAD+ biosynthesis metabolism is implicated in tumor progression. Breast cancer (BC) is the most common malignant malignancy in the world. Nevertheless, the prognostic significance of NAD+ biosynthesis and its relationship with the tumor immune microenvironment in breast cancer still need further investigation. In this study, we obtained the mRNA expression data and clinical information of BC samples from public databases and calculated the level of NAD+ biosynthesis activity by single-sample gene set enrichment analysis (ssGSEA). We then explored the relationship between the NAD+ biosynthesis score, infiltrating immune cells, prognosis significance, immunogenicity and immune checkpoint molecules. The results demonstrated that patients with high NAD+ biosynthetic score displayed poor prognosis, high immune infiltration, high immunogenicity, elevated PD-L1 expression, and might more benefit from immunotherapy. Taken together, our studies not only deepened the understanding of NAD+ biosynthesis metabolism of breast cancer but also provided new insights into personalized treatment strategies and immunological therapies to improve the outcomes of breast cancer patients.
Introduction
Breast cancer (BC) is the most common malignancy in the world and the fifth leading cause of cancer death worldwide (1). With advance in the diagnosis and treatment of breast cancer, the survival rate of breast cancer patients has significantly improved (2). However, as a malignancy with high heterogeneity at the molecular level, BC patients with similar clinical features may have different prognosis (3). Therefore, it is essential to consider other important factors in guiding clinical practice.
The tumor microenvironment (TME) is a highly complex cellular network, including tumor cells, stromal cells, fibroblasts, immune cells, soluble factors, signaling molecules and extracellular matrix components (4). Recent studies suggest that the immune components in the TME are closely associated with tumor development, recurrence, and metastasis. These components are called the tumor immune microenvironment (TIME) (5). Immune cells or immune-related genes in TIME can predict the prognosis and treatment efficacy of cancer patients (6–8). For instance, tumor-infiltrating lymphocytes (TILs) have been shown to have a strong prognostic effect in patients with early stage TNBC and HER2-positive breast cancer (9, 10). Exploring TIME in BC would conductive to guide and optimize immunotherapy and improve the prognosis of BC patients.
NAD+ (nicotinamide adenine dinucleotide) is a coenzyme for redox reactions and a substrate for different signaling enzymes that can directly or indirectly affect many important cellular functions, such as metabolic pathways, DNA repair, and immune cell function (11, 12). Many of these processes are associated with cancer development. Given that NAD+ -dependent signaling responses involve the degradation of molecules, sustained NAD+ production via different biosynthetic pathways is a hallmark of many types of tumor (13, 14). Enzymes involved in NAD+ biosynthesis have been reported to be aberrantly expressed or dysregulated in many cancer types, including breast cancer. NAMPT, one of the rate-limiting enzymes of NAD+ biosynthesis, has been shown to be overexpressed in breast cancer and associated with breast cancer proliferation and invasiveness (15, 16). However, the relationship between NAD+ biosynthesis and the immune microenvironment of breast cancer has not been systematically investigated.
In the current study, a series of bioinformatic methods were used to analyze the features of NAD+ biosynthesis in BC based on transcriptional profiling data from public databases. BC patients were classified into the high and low NAD+ biosynthetic subtypes based on the NAD+ biosynthesis score and survival analysis was performed. We further explored the hallmark pathways and TME immune cell infiltration characteristics of the two subtypes. Moreover, we analyzed the relationship between the NAD+ biosynthesis score and several immunotherapy biomarkers. These findings may provide a new perspective for exploring the metabolic mechanism and treatment of breast cancer.
Materials and methods
Retrieval of NAD+ biosynthesis-related genes
NAD+ biosynthesis-related genes were obtained from the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database (Pathway: hsa00760) and Reactome database (R-HSA-196807) (17). We also reviewed the literature and added the previously reported genes (11). After selecting genes associated with NAD+ biosynthesis from all gene sets, a total of 19 genes were retrieved (Supplementary Table S1).
Datasets and data preprocessing
The mRNA sequencing, gene somatic mutations, and clinical data of breast invasive carcinoma (BRCA) patients were downloaded from The Cancer Genome Atlas (TCGA) database and included data form 99 normal samples and 1069 tumor samples. RNA-sequencing data (FPKM values) were transformed into transcripts per million (TPM) values and normalized into log2 (TPM +1) for the following analysis. Datasets of GSE20711, GSE48390 and GSE88770 from the Gene Expression Omnibus (GEO) were selected to validate the results of TCGA data analysis. The GSE20711 dataset was used by Dedeurwaerder et al. to study epigenetic variation (methylation) related to gene expression in breast cancer (18). GSE20711 consists of 88 breast cancer and 2 normal breast tissue samples, with only the gene expression profiles of 88 breast cancer patients retained for subsequent analysis. The GSE48390 dataset was used by Huang et al. to identify genes with coherent patterns of both copy number variation (CNV) and differential gene expression, and to use these genes to derive signatures related to clinical ER and HER2 status and disease-free survival (19). GSE48390 consists of 81 breast cancer samples. The GSE88770 dataset was used by Metzger-Filho et al. to analyze the prognostic value of histological grading (HG) in breast cancer (20). GSE88770 contains 117 breast cancer samples. The three datasets were based on GPL570 platforms (Affymetrix Human Genome U133 Plus 2.0 Array). We downloaded the original expression profile and used the robust multi-array average (RMA) algorithm to perform background correction and quantile normalization (21). We merged GEO datasets and used the combat algorithm to eliminate the batch effects by R package “sva” (22).
Evaluation of tumor NAD+ biosynthesis score
We applied the single-sample gene-set enrichment analysis (ssGSEA) for the NAD+ biosynthesis gene set to quantify the NAD+ biosynthesis activity (NAD+ biosynthesis score) (23). The optimal cutoff threshold for classifying breast cancer (BC) patients into high and low NAD+ biosynthetic subtypes was determined using the “surv_cutpoint” function from the R package “survminer”.
Differential gene expression analysis and enrichment analysis
We identified differentially expressed genes (DEGs) between the high and low NAD+ biosynthetic subtypes using the R package “limma,” with a significance threshold of adjusted p-value <0.05 and |log2FoldChange| > 1. To assess the differences in biological signaling pathway between the high and low NAD+ biosynthetic subtypes, GSEA analysis was performed using the gseKEGG function of the R package “clusterProfiler,” with p < 0.05 considered statistically significant (24).
Immune infiltration analysis
Immune cell gene signatures and immune-related signatures were collected from previously published works (25, 26), and ssGSEA was used to calculate the enrichment scores of infiltrating immune cells and immune functions. In addition, the “ESTIMATE” package was used to assess the composition of the immune stroma in the tumor microenvironment of breast cancer patients, and the immune score, stromal score, and estimated score (ESTIMATE Score) were calculated (27).
Construction of TF-IRG networks
We extracted differentially expressed immune-related genes (DEIRGs) and transcription factors (DETFs) from DEGs based on the lists obtained from ImmPort and Cistrome Cancer database (Supplementary Tables S2, S3). We then calculated the correlation between DETFs and DEIRGs. Correlation coefficients >0.5 (or less than −0.5) and adjusted p-value <0.01 were considered to be significantly correlated (Supplementary Table S4).
Mutation and evaluation of the therapeutic efficacy
The R package “maftools” was used to estimate the tumour mutation burden (TMB) for each patient between the two subtypes (28). Additionally, the neoantigen load and homologous recombination deficiency (HRD) in BRCA patients were collected from published studies (29, 30). For the immunotherapy sensitivity prediction analysis, the immunophenoscore (IPS) of breast cancer patients was downloaded from the TCIA database (25). For drug sensitivity prediction analysis, the “pRRophetic” R package was used to calculate the IC50 of common chemotherapeutic agents (31). The difference in the IC50 of agents between the high and low NAD+ biosynthetic subtypes was evaluated with the Wilcoxon test, and p-value <0.05 was considered significant.
Statistical and computational analysis
All statistical analyses were done in R (version 4.1). The significance of differences in continuous variables between the two subtypes was calculated using the Wilcoxon test. The Kruskal-Wallis test and Dunn’s post hoc multiple comparison test were used for more than two subtypes. The optimal cutoff values for each cohort were evaluated using the “surv_cutpoint” function in the “survminer” package. For prognostic analysis, the survival curves were constructed using the Kaplan-Meier method and log-rank test was used to judge differences between the subtypes. Spearman correlation analysis was carried out to determine the correlation coefficient. Statistical significance was considered as p-value <0.05.
Results
Identification of NAD+ biosynthetic subtypes and their correlation with biological functions in BRCA
NAD+ is essential for the proper function and metabolism of all living cells, including cancer cells. We used ssGSEA to calculate the NAD+ biosynthesis score of patients from TCGA database and found that NAD+ biosynthesis score of breast cancer tissue was significantly higher than that of normal breast tissue, which is consistent with previous findings (Supplementary Figure S1). To better understand the effect of NAD+ biosynthesis on breast cancer prognosis, the BC patients were divided into the high and low NAD+ biosynthetic subtypes based on the optimal cut-off value (−0.0743) in TCGA cohort, and subjected to survival comparison. Figure 1A showed that patients in the low NAD+ biosynthetic subtype had significantly better overall survival and disease specific survival than patients in the high NAD+ biosynthetic subtype. The prognostic significance of the NAD+ biosynthesis score was validated in GEO cohort (Supplementary Figure S2A). Compared to low NAD+ biosynthetic subtype, twelve NAD+ biosynthesis-related genes were significantly upregulated in the high NAD+ biosynthetic subtype, which indicated that NAD+ biosynthesis score may reflect the activity of NAD+ biosynthesis to some extent (Figure 1B). Next, we analyzed the clinical features differences between the two subtypes, including pathological stage, molecular subtypes and TNM stage (Supplementary Figure S3). We found that only molecular subtypes and NAD+ biosynthesis score exhibited statistical significance (Figures 1C, D). In detail, the HER2-enriched subtype had significantly higher NAD+ biosynthesis score than other four molecular subtypes, whereas the LumA subtype had the lowest NAD+ biosynthesis score.
FIGURE 1
FIGURE 1. Subtypes of NAD+ biosynthesis and biological characteristics in BRCA. (A) Survival analysis of the high and low NAD+ biosynthetic subtypes in TCGA cohort. (B) Differential expression of NAD+ biosynthesis-related genes between the high and low NAD+ biosynthetic subtypes in TCGA cohort. (C) Violin plots showed the correlation between the NAD+ biosynthesis score and molecular subtypes in TCGA cohort. (D) Alluvial plot for the NAD+ biosynthetic subtypes versus different molecular subtypes in TCGA cohort. (E) Differential pathway activities scored by GSEA between the high and low biosynthetic subtypes in TCGA cohort. The red bars indicated the upregulated pathways, while the blue bars indicated the downregulated pathways. (F) Correlation between the NAD+ biosynthesis score and ssGSEA enrichment scores of cancer hallmark pathways in TCGA cohort. *p < 0.05; **p < 0.01; ***p < 0.001, and ****p < 0.0001.
GSEA was used to investigate the biological differences between the two subtypes. As shown in Figure 1E, the results showed that many immune-related pathways were significantly upregulated in the high NAD+ biosynthetic subtype, including antigen processing and presentation, IL-17 signaling pathway and chemokine signaling pathway. The low NAD+ biosynthetic subtype was enriched in pathways related to ECM-receptor interaction and regulation of lipolysis in adipocytes. Besides, to explore the potential pathways which NAD+ biosynthesis was involved in, the correlation between the NAD+ biosynthesis score and cancer hallmarks was analyzed. We found a significant positive association between the NAD+ biosynthesis score and some hallmark pathways, including reactive oxygen species, glycolysis, and DNA repair, which may explain the inconsistent clinical outcome between the two subtypes. We also found that the NAD+ biosynthesis score was positively correlated with immune-related pathways such as allograft rejection, IFN-γsignaling, inflammatory response. While pathways such as estrogen response, TGF-beta and hedgehog signaling were negatively correlated with the NAD+ biosynthesis score (Figure 1F). Similar biological differences between the two subtypes were observed in GEO cohort (Supplementary Figures S2C, D). Collectively, these findings suggested that NAD+ biosynthesis metabolism may play an important role not only in the development of breast cancer, but also in tumor immune environment.
Immune landscape of NAD+ biosynthetic subtypes
The state of TIME determines the fate of cancer cells. As indicated by the above results, NAD+ biosynthesis had a certain connection with immune response. As such, we assessed the differences in the immune microenvironment of patients in the two subtypes. The immune score, stromal score, and ESTIMATE score for the two subtypes were calculated using the ESTIMATE algorithm. The findings suggested that high NAD+ biosynthetic subtype had a higher ESTIMATE score, immune score, but a lower stromal score compared to the low NAD+ biosynthetic subtype (Figure 2A). We also found a significant positive correlation between the NAD+ biosynthesis score and immune score in breast cancer (Figure 2B), indicating that increased NAD+ biosynthesis was associated with immune activation and high immune infiltration in breast cancer. By comparing the ssGSEA score of 28 types of immune cells using the Wilcoxon test, we found that most immune cells, such as activated B cells, activated CD8+ T cells, activated dendritic cells, myeloid-derived suppressor cells (MDSCs) and regulatory T cells, were significantly enriched in the high NAD+ biosynthetic subtype, while few immune cells, such as central memory CD8 T cells and mast cells, were significantly enriched in the low NAD+ biosynthetic subtype (Figure 2C). The NAD+ biosynthesis score was observed closely associated with various categories of infiltrating immune cell, especially MDSCs (Figure 2D). The findings were similar in GEO cohort (Supplementary Figure S4). MDSCs are important components of TIME and can suppress the anti-tumor functions of T cells and natural killer cells (32). We believed that this immunosuppressive cell type may play an important role in the adverse clinical outcome of patients with high NAD+ biosynthesis score. The correlation analysis showed that most NAD+ biosynthesis-related genes, including NMRK1, PNP, IDO1 and NAMPT, were positively correlated with most immune cells, while negative correlations were demonstrated in molecules such as NMNAT2 and AFMID (Figure 2E). Moreover, the scores of most immune-related pathways were significantly higher in the high NAD+ biosynthetic subtype, including major histocompatibility complex (MHC) class I, and interferon response (Figure 2F). These results suggested that patients with high NAD+ biosynthetic activity tend to present a high immune infiltration and suppressive immune microenvironment.
FIGURE 2
FIGURE 2. Correlation of the NAD+ biosynthesis score with the immune landscape of patients with breast cancer. (A) Comparison of the stromal score, immune score, and ESTIMATE score between the high and low NAD+ biosynthetic subtypes. (B) The correlation between the NAD+ biosynthesis score and immune score. (C) Relative infiltration of 28 types of immune cells in the high and low NAD+ biosynthetic subtypes. (D) The correlation between the NAD+ biosynthesis score and the ssGSEA enrichment scores of immune cells. (E) The heatmap showed the associations between the 19 NAD+ biosynthesis-related genes and immune cells. Red indicated positive correlations, and blue indicated negative correlations. Asterisks denoted p-value. Blank cells represented no statistical significance of the correlation. (F) Relative enrichment score of 12 immune-related signatures in the high and low NAD+ biosynthetic subtypes. *p < 0.05; **p < 0.01; ***p < 0.001, and ****p < 0.0001.
Potential TF-IRG regulatory network
To further explore the potential regulatory mechanism of NAD+ biosynthesis affecting TIME, we identified 59 differentially expressed immune-related genes (DEIRGs), including 45 upregulated DEIRGs and 14 downregulated DEIRGs, which may play a major role in regulating TIME. DEIRGs are strongly associated with chemokine-mediated signaling (Figure 3A). Enhanced expression levels of T-cells recruiting chemokines such as CXCL9, CXCL10, and CXCL11 were found in the high NAD+ biosynthetic subtype, consistent with their abundant T cells infiltration. Additionally, the overexpression of CXCL13 (which can recruit MDSC) and other tumor suppressive chemokines (e.g., CCL18, CCL8 and CCL20) was also observed in the high NAD+ biosynthetic subtype (Supplementary Table S2). Transcriptional factors (TFs) can control the expression of critical genes and thus play an important role in the regulation of tumor inflammation and immunity (33, 34). To explore the potential upstream regulatory mechanism of DEIRGs, 11 differentially expressed TFs (DETFs) were also found from DEGs (Supplementary Table S3). A TF-IRG network was constructed based on the correlation of gene expression between DETFs and DEIRGs. Only gene pairs with correlation coefficients >0.5 (or less than −0.5) and adjusted p-value <0.01 were incorporated into the network (Supplementary Table S4). The TF-IRG network was visualized in Cytoscape (Figure 3B). We found that two TF genes, BATF2 and SPIB were significantly elevated in the high NAD+ biosynthetic subtype. Notably, SPIB was significantly associated with PDCD1 and CTLA4 (Figures 3C, D). The relationships between SPIB, PDCD1 and CTLA4 were also verified in GEO cohort (Supplementary Figure S5B). SPIB has been reported to play a negative role in immune cell regulation (35). Therefore, the SPIB-centered regulatory network may play an important role in immune infiltration and immune escape in the high NAD+ biosynthetic subtypes. We further confirmed that the expression level of SPIB was significantly higher in tumor tissues than in adjacent normal tissues (Figure 3E).
FIGURE 3
FIGURE 3. Construction of the TF-IRG network in BRCA. (A) GO enrichment analysis of 59 differentially expressed immune-related genes (DEIRGs). (B) A TF-IRG network in TCGA cohort. The pink line is the positive correlation between differentially expressed transcription factors (DETFs) and differentially expressed immune-related genes (DEIRGs). The green line is the negative correlation between DETFs and DEIRGs. (C) The correlation between the expression of SPIB and PDCD1 in TCGA cohort. (D) The correlation between the expression of SPIB and CTLA4 in TCGA cohort. (E) Differential expression of SPIB between breast normal and tumor tissues.
Mutation and immunotherapy sensitivity prediction of NAD+ biosynthetic subtypes
Immunotherapy has become an emerging clinical strategy for treating cancer (36). Considering that enhanced tumor immunogenicity predicts long-term clinical benefits for patients from immune checkpoint inhibitors (ICIs), we wanted to explore the relationship between NAD+ biosynthesis scores and tumor immunogenicity. Therefore, we took the perspective of the tumor mutational burden (TMB), neoantigen load and homologous recombination deficiency (HRD). Results showed that patients in the high NAD+ biosynthetic subtype had higher TMB, neoantigen and HRD (Figures 4A–C). The subsequent correlation analysis showed that TMB, neoantigen load and HRD were positively correlated with the NAD+ biosynthesis score (Supplementary Figure S6). These results suggested that patients in the high NAD+ biosynthetic subtype had relatively high immunogenicity. Somatic mutation analysis revealed higher TP53 mutation rates in the high NAD+ biosynthetic subtype and higher PIK3CA mutation rates in the low NAD+ biosynthetic subtype. (Figures 4D, E). Additionally, the expression of immune checkpoint molecules such as PD-L1, LAG3, PD-1, and CTLA-4 was higher in the high NAD+ biosynthetic subtype (Figure 4F) We also assessed the association between IPS and NAD+ biosynthesis score. The results illustrated that the IPS for PD1/PD-L1/PD-L2 blocker, CTLA-4 blocker and PD1/PD-L1/PD-L2 plus CTLA-4 blocker in high NAD+ biosynthetic subtype were significantly higher than those in low NAD+ biosynthetic subtype (Figure 4G). These findings suggested that patients with high NAD+ biosynthetic score may be appropriate candidates for immunotherapy.
FIGURE 4
FIGURE 4. Mutation and Immunotherapy Sensitivity Prediction of the NAD+ biosynthesis signature. (A–C) TMB, Neoantigen and HRD between the high and low NAD+ biosynthetic subtypes. (D,E) Waterfall plot of tumor somatic mutation in the high and low NAD+ biosynthetic subtypes. (F) Comparison of immune checkpoints expression levels between the high and low NAD+ biosynthetic subtypes. (G) Comparison of the immunophenoscore (IPS) between the high and low NAD+ biosynthetic subtypes. *p < 0.05; **p < 0.01; ***p < 0.001, and ****p < 0.0001.
Drug sensitivity prediction
Considering that drug chemosensitivity can influence the clinical outcomes of breast cancer treatment, we predicted the IC50 of usual chemotherapeutic drugs and compared them between the two subtypes. Results showed that the estimated IC50 for doxorubicin was lower in the low NAD+ biosynthetic subtype, while the estimated IC50s for docetaxel, paclitaxel and gefitinib were lower in the high NAD+ biosynthetic subtype. However, there was no significant difference in the estimated IC50s for cisplatin and etoposide (Figures 5A–F). These results suggested that the NAD+ biosynthesis score may differentiate more patients for appropriate therapy.
FIGURE 5
FIGURE 5. Drug Sensitivity Prediction. (A–F) Box plots of the estimated IC50 for Cisplatin (A), Doxorubicin (B), Docetaxel (C), Etoposide(D), Paclitaxel (E) and Gefitinib (F). *p < 0.05; **p < 0.01; ***p < 0.001.
Discussion
The role of NAD+ in cellular metabolism and signal transduction is crucial as it cannot be obtained directly from food, making NAD+ biosynthesis essential (12). Previous works have shown that NAD+ is critical for the development and progression of breast cancer. For example, QPRT, one of the rate-limiting enzymes of NAD+ biosynthesis, was shown to potentially enhance breast cancer invasiveness through purinergic signaling (37). However, little work has explored the relationship between NAD+ biosynthesis metabolism and TIME in breast cancer. Thus, in this study, we quantified NAD+ biosynthesis metabolism based on ssGSEA and demonstrated that the NAD+ biosynthesis score was significantly associated with the overall survival. Further bioinformatics analyses revealed that the NAD+ biosynthesis score was significantly related to different tumor-infiltrating immune cells and other immune-related signatures. Interestingly, we also found that the NAD+ biosynthesis score was positively associated with several immunotherapy biomarkers. This indicates that patients in the high NAD+ biosynthesis subtype may be more suitable for immunotherapy. As far as we know, this study is the first study to systematically explore the relationship between NAD+ biosynthesis metabolism and TIME in breast cancer.
In this study, we divided patients with breast cancer into the high and low biosynthetic subtypes based on the NAD+ biosynthesis score. We found that high NAD+ biosynthetic subtype was associated with poor prognosis. This correlation was validated in GEO cohort. GSEA analysis revealed that the NAD+ biosynthesis score was positively correlated with immune-related biological processes and pathways, suggesting that the difference in survival between the two subtypes may be due to immune heterogeneity. Then, analysis of the immune microenvironment showed that high NAD+ biosynthetic subtype was characterized by immune activation and accompanied with immune suppression. The ssGSEA scores of most immune cells in high NAD+ biosynthetic subtypes were higher than those in low NAD+ biosynthetic subtypes, which indicated that the immune microenvironment in BC with high biosynthesis score tended to be “hot.” Our results also showed that HER2-enriched and Basal subtypes had higher NAD+ biosynthesis score, which were consistent with higher TILs infiltration in the two molecular subtypes (38). However, high NAD+ biosynthetic subtype did not display a corresponding survival advantage and we found MDSC might play a key role. MDSC could promote immune evasion, tumor angiogenesis and metastasis (39). Besides, there is a positive association between IFN-γ signaling and NAD+ biosynthesis score. IFN-γ has a dual role. On the one hand, IFN-γ can enhance anticancer activity by increasing MHC class I and cytotoxic proteins associated with the CTL response. On the other hand, it also promotes the expression of PD-L1, CTLA-4, IDO1, CXCL12 and other molecules to promote tumor immune escape(40, 41). In fact, it has been reported that NAD+ could enhance IFN-γ-induced PD-L1 expression and promote tumor immune evasion (42). Therefore, we believe that IFN-γ signaling and MDSC could contribute to the immunosuppressive microenvironment in NAD+ biosynthesis-activated breast cancer.
Besides, the TF-IRG regulatory network implied that two transcription factors, BATF2 and SPIB, may play an active role in transcriptional events associated with the infiltration of immune cells in the high NAD+ biosynthetic subtype. BATF2, a member of the basic leucine zipper transcription factor family, can be involved in the gene regulation of IFN-γ-activated classical macrophages and induce pro-inflammatory responses (43, 44). SPIB, belongs to the ETS family, plays an important role in the differentiation of B cells and plasmacytoid dendritic cells (45, 46). SPIB has been shown to promote tumor aerobic glycolysis (47) and participate in the recruitment of tumor-associated macrophages (TAMs) (35), which promotes cancer progression. In conclusion, the construction of the TF-IRG network provided some insights into the interaction between NAD+ biosynthesis and TIME in breast cancer.
Immunotherapy, represented by immune checkpoint inhibitors (ICIs), is changing the treatment of cancer. The higher immunogenicity in the high NAD+ biosynthetic subtype may lead to a better prognosis after receiving immunotherapy (48). TMB is a more reliable biomarker for predicting immunotherapy efficacy. Patients with high TMB will achieve higher objective response rates to immunotherapy (49, 50). The positive association between HRD and TMB has been reported in BC (51). The underlying mechanism for this correlation is that HRD may drive tumorigenesis, increasing the number of tumor mutations and the neoantigen rate in BC (52, 53). Neoantigens can promote CD8+ T cell infiltration and cytolytic activity, which are closely associated with response to immunogenicity. Mutation analysis showed high NAD+ biosynthesis was associated with high TP53 mutation rate. One study reported that patients with mutant TP53 showed stronger tumor antigenicity and tumor antigen presentation and were more likely to benefit from immunotherapy (54), which was consistent with our findings. Moreover, we also found that patients with different NAD+ biosynthesis scores might exhibit distinct sensitivity to doxorubicin, docetaxel, paclitaxel, and gefitinib. These results suggested that NAD+ biosynthesis score can assess the efficacy of immunotherapy and chemotherapy, facilitating personalized treatment for breast cancer.
This study has some limitations. Firstly, in our study, the clinical cohorts were derived from public databases and need further validated by prospective studies. Secondly, the association between NAD+ biosynthesis and TIME requires additional experimental validation.
In conclusion, our study suggests that NAD+ biosynthetic activity was correlated with the prognosis of breast cancer and the immune microenvironment. NAD+ biosynthetic activity may also serve as a potential prognostic biomarker for ICI immunotherapy, which may help to screen patients eligible for immunotherapy and guide future individualized precision therapy.
Data availability statement
The public datasets to support the results of this research can be gained from TCGA (https://portal.gdc.cancer.gov/) and GEO (https://www.ncbi.nlm.nih.gov/geo/).
Author contributions
YY, ZW and MH drafted the manuscript. YY, LD and BY contributed to the data collection and analysis. DL designed and supervised this study. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by National key Research and Development Program of China (2020YFE0202200, 2021YFA1301603), National Natural Science Foundation of China (32088101, 32271518).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors express sincere gratitude for the availability of data by public databases such as ImmPort, GEO and TCGA.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.por-journal.com/articles/10.3389/pore.2023.1610956/full#supplementary-material
References
1. Sung, H, Ferlay, J, Siegel, RL, Laversanne, M, Soerjomataram, I, Jemal, A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin (2021) 71:209–49. doi:10.3322/caac.21660
2. McDonald, ES, Clark, AS, Tchou, J, Zhang, P, and Freedman, GM. Clinical diagnosis and management of breast cancer. J Nucl Med Off Publ Soc Nucl Med (2016) 57(1):9S–16S. doi:10.2967/jnumed.115.157834
3. Januškevičienė, I, and Petrikaitė, V. Heterogeneity of breast cancer: The importance of interaction between different tumor cell populations. Life Sci (2019) 239:117009. doi:10.1016/j.lfs.2019.117009
4. Catalano, V, Turdo, A, Di Franco, S, Dieli, F, Todaro, M, and Stassi, G. Tumor and its microenvironment: A synergistic interplay. Semin Cancer Biol (2013) 23:522–32. doi:10.1016/j.semcancer.2013.08.007
5. Fu, T, Dai, L-J, Wu, S-Y, Xiao, Y, Ma, D, Jiang, YZ, et al. Spatial architecture of the immune microenvironment orchestrates tumor immunity and therapeutic response. J Hematol Oncolj Hematol Oncol (2021) 14:98. doi:10.1186/s13045-021-01103-4
6. Wang, S, Zhang, Q, Yu, C, Cao, Y, Zuo, Y, and Yang, L. Immune cell infiltration-based signature for prognosis and immunogenomic analysis in breast cancer. Brief Bioinform (2021) 22:2020–31. doi:10.1093/bib/bbaa026
7. Oshi, M, Asaoka, M, Tokumaru, Y, Yan, L, Matsuyama, R, Ishikawa, T, et al. CD8 T cell score as a prognostic biomarker for triple negative breast cancer. Int J Mol Sci (2020) 21:E6968. doi:10.3390/ijms21186968
8. Gao, W, Wen, H, Liang, L, Dong, X, Du, R, Zhou, W, et al. IL20RA signaling enhances stemness and promotes the formation of an immunosuppressive microenvironment in breast cancer. Theranostics (2021) 11:2564–80. doi:10.7150/thno.45280
9. Savas, P, Salgado, R, Denkert, C, Sotiriou, C, Darcy, PK, Smyth, MJ, et al. Clinical relevance of host immunity in breast cancer: From TILs to the clinic. Nat Rev Clin Oncol (2016) 13:228–41. doi:10.1038/nrclinonc.2015.215
10. Dieci, MV, Miglietta, F, and Guarneri, V. Immune infiltrates in breast cancer: Recent updates and clinical implications. Cells (2021) 10:223. doi:10.3390/cells10020223
11. Covarrubias, AJ, Perrone, R, Grozio, A, and Verdin, E. NAD+ metabolism and its roles in cellular processes during ageing. Nat Rev Mol Cel Biol (2021) 22:119–41. doi:10.1038/s41580-020-00313-x
12. Navas, LE, and Carnero, A. NAD+ metabolism, stemness, the immune response, and cancer. Signal Transduct Target Ther (2021) 6:2. doi:10.1038/s41392-020-00354-w
13. Chiarugi, A, Dölle, C, Felici, R, and Ziegler, M. The NAD metabolome--a key determinant of cancer cell biology. Nat Rev Cancer (2012) 12:741–52. doi:10.1038/nrc3340
14. Ghanem, MS, Monacelli, F, and Nencioni, A. Advances in NAD-lowering agents for cancer treatment. Nutrients (2021) 13:1665. doi:10.3390/nu13051665
15. Zhang, H, Zhang, N, Liu, Y, Su, P, Liang, Y, Li, Y, et al. Epigenetic regulation of NAMPT by NAMPT-AS drives metastatic progression in triple-negative breast cancer. Cancer Res (2019) 79:3347–59. doi:10.1158/0008-5472.CAN-18-3418
16. Gholinejad, Z, Kheiripour, N, Nourbakhsh, M, Ilbeigi, D, Behroozfar, K, Hesari, Z, et al. Extracellular NAMPT/Visfatin induces proliferation through ERK1/2 and AKT and inhibits apoptosis in breast cancer cells. Peptides (2017) 92:9–15. doi:10.1016/j.peptides.2017.04.007
17. Li, C, Zhu, Y, Chen, W, Li, M, Yang, M, Shen, Z, et al. Circulating NAD+ metabolism-derived genes unveils prognostic and peripheral immune infiltration in amyotrophic lateral sclerosis. Front Cel Dev Biol (2022) 10:831273. doi:10.3389/fcell.2022.831273
18. Dedeurwaerder, S, Desmedt, C, Calonne, E, Singhal, SK, Haibe-Kains, B, Defrance, M, et al. DNA methylation profiling reveals a predominant immune component in breast cancers. EMBO Mol Med (2011) 3:726–41. doi:10.1002/emmm.201100801
19. Huang, C-C, Tu, S-H, Lien, H-H, Jeng, JY, Huang, CS, Huang, CJ, et al. Concurrent gene signatures for han Chinese breast cancers. PLoS One (2013) 8:e76421. doi:10.1371/journal.pone.0076421
20. Metzger-Filho, O, Michiels, S, Bertucci, F, CatteAu, A, Salgado, R, Galant, C, et al. Genomic grade adds prognostic value in invasive lobular carcinoma. Ann Oncol Off J Eur Soc Med Oncol (2013) 24:377–84. doi:10.1093/annonc/mds280
21. Gautier, L, Cope, L, Bolstad, BM, and Irizarry, RA. affy--analysis of Affymetrix GeneChip data at the probe level. Bioinforma Oxf Engl (2004) 20:307–15. doi:10.1093/bioinformatics/btg405
22. Johnson, WE, Li, C, and Rabinovic, A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostat Oxf Engl (2007) 8:118–27. doi:10.1093/biostatistics/kxj037
23. Hänzelmann, S, Castelo, R, and Guinney, J. Gsva: Gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics (2013) 14:7. doi:10.1186/1471-2105-14-7
24. Yu, G, Wang, L-G, Han, Y, and He, Q-Y. clusterProfiler: an R Package for comparing biological themes among gene clusters. OMICS J Integr Biol (2012) 16:284–7. doi:10.1089/omi.2011.0118
25. Charoentong, P, Finotello, F, Angelova, M, Mayer, C, Efremova, M, Rieder, D, et al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep (2017) 18:248–62. doi:10.1016/j.celrep.2016.12.019
26. Bindea, G, Mlecnik, B, Tosolini, M, Kirilovsky, A, Waldner, M, Obenauf, AC, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity (2013) 39:782–95. doi:10.1016/j.immuni.2013.10.003
27. Yoshihara, K, Shahmoradgoli, M, Martínez, E, Vegesna, R, Kim, H, Torres-Garcia, W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun (2013) 4:2612. doi:10.1038/ncomms3612
28. Mayakonda, A, Lin, D-C, Assenov, Y, Plass, C, and Koeffler, HP. Maftools: Efficient and comprehensive analysis of somatic variants in cancer. Genome Res (2018) 28:1747–56. doi:10.1101/gr.239244.118
29. Shi, Z, Shen, J, Qiu, J, Zhao, Q, Hua, K, and Wang, H. CXCL10 potentiates immune checkpoint blockade therapy in homologous recombination-deficient tumors. Theranostics (2021) 11:7175–87. doi:10.7150/thno.59056
30. Thorsson, V, Gibbs, DL, Brown, SD, Wolf, D, Bortone, DS, Ou Yang, TH, et al. The immune landscape of cancer. Immunity (2018) 48:812–30.e14. doi:10.1016/j.immuni.2018.03.023
31. Geeleher, P, Cox, N, and Huang, RS. pRRophetic: an R package for prediction of clinical chemotherapeutic response from tumor gene expression levels. PLoS One (2014) 9:e107468. doi:10.1371/journal.pone.0107468
32. Hinshaw, DC, and Shevde, LA. The tumor microenvironment innately modulates cancer progression. Cancer Res (2019) 79:4557–66. doi:10.1158/0008-5472.CAN-18-3962
33. Yu, H, Pardoll, D, and Jove, R. STATs in cancer inflammation and immunity: A leading role for STAT3. Nat Rev Cancer (2009) 9:798–809. doi:10.1038/nrc2734
34. Marigo, I, Bosio, E, Solito, S, Mesa, C, Fernandez, A, Dolcetti, L, et al. Tumor-induced tolerance and immune suppression depend on the C/EBPbeta transcription factor. Immunity (2010) 32:790–802. doi:10.1016/j.immuni.2010.05.010
35. Huang, Q, Liu, J, Wu, S, Zhang, X, Xiao, Z, Liu, Z, et al. Spi-B promotes the recruitment of tumor-associated macrophages via enhancing CCL4 expression in lung cancer. Front Oncol (2021) 11:659131. doi:10.3389/fonc.2021.659131
36. Zhao, Z, Zheng, L, Chen, W, Weng, W, Song, J, and Ji, J. Delivery strategies of cancer immunotherapy: Recent advances and future perspectives. J Hematol Oncolj Hematol Oncol (2019) 12:126. doi:10.1186/s13045-019-0817-3
37. Liu, C-L, Cheng, S-P, Chen, M-J, Lin, CH, Chen, SN, Kuo, YH, et al. Quinolinate phosphoribosyltransferase promotes invasiveness of breast cancer through myosin light chain phosphorylation. Front Endocrinol (2020) 11:621944. doi:10.3389/fendo.2020.621944
38. Hammerl, D, Smid, M, Timmermans, AM, Sleijfer, S, Martens, JWM, and Debets, R. Breast cancer genomics and immuno-oncological markers to guide immune therapies. Semin Cancer Biol (2018) 52:178–88. doi:10.1016/j.semcancer.2017.11.003
39. Gabrilovich, DI. Myeloid-derived suppressor cells. Cancer Immunol Res (2017) 5:3–8. doi:10.1158/2326-6066.CIR-16-0297
40. Lin, C-F, Lin, C-M, Lee, K-Y, Wu, SY, Feng, PH, Chen, KY, et al. Escape from IFN-γ-dependent immunosurveillance in tumorigenesis. J Biomed Sci (2017) 24:10. doi:10.1186/s12929-017-0317-0
41. Mojic, M, Takeda, K, and Hayakawa, Y. The dark side of IFN-γ: Its role in promoting cancer immunoevasion. Int J Mol Sci (2017) 19:89. doi:10.3390/ijms19010089
42. Lv, H, Lv, G, Chen, C, Zong, Q, Jiang, G, Ye, D, et al. NAD+ metabolism maintains inducible PD-L1 expression to drive tumor immune evasion. Cell Metab (2021) 33:110–27.e5. doi:10.1016/j.cmet.2020.10.021
43. Guler, R, Mpotje, T, Ozturk, M, Nono, JK, Parihar, SP, Chia, JE, et al. Batf2 differentially regulates tissue immunopathology in Type 1 and Type 2 diseases. Mucosal Immunol (2019) 12:390–402. doi:10.1038/s41385-018-0108-2
44. Roy, S, Guler, R, Parihar, SP, Schmeier, S, Kaczkowski, B, Nishimura, H, et al. Batf2/Irf1 induces inflammatory responses in classically activated macrophages, lipopolysaccharides, and mycobacterial infection. J Immunol Baltim (2015) 194:6035–44. doi:10.4049/jimmunol.1402521
45. Garrett-Sinha, LA, Su, GH, Rao, S, Kabak, S, Hao, Z, Clark, MR, et al. PU.1 and Spi-B are required for normal B cell receptor-mediated signal transduction. Immunity (1999) 10:399–408. doi:10.1016/s1074-7613(00)80040-0
46. Schotte, R, Nagasawa, M, Weijer, K, Spits, H, and Blom, B. The ETS transcription factor Spi-B is required for human plasmacytoid dendritic cell development. J Exp Med (2004) 200:1503–9. doi:10.1084/jem.20041231
47. Wang, J, Wang, X, Guo, Y, Ye, L, Li, D, Hu, A, et al. Therapeutic targeting of SPIB/SPI1-facilitated interplay of cancer cells and neutrophils inhibits aerobic glycolysis and cancer progression. Clin Transl Med (2021) 11:e588. doi:10.1002/ctm2.588
48. Wang, S, He, Z, Wang, X, Li, H, and Liu, XS. Antigen presentation and tumor immunogenicity in cancer immunotherapy response prediction. eLife (2019) 8:e49020. doi:10.7554/eLife.49020
49. Chan, TA, Yarchoan, M, Jaffee, E, Swanton, C, Quezada, SA, Stenzinger, A, et al. Development of tumor mutation burden as an immunotherapy biomarker: Utility for the oncology clinic. Ann Oncol Off J Eur Soc Med Oncol (2019) 30:44–56. doi:10.1093/annonc/mdy495
50. Barroso-Sousa, R, Jain, E, Cohen, O, Kim, D, Buendia-Buendia, J, WinEr, E, et al. Prevalence and mutational determinants of high tumor mutation burden in breast cancer. Ann Oncol Off J Eur Soc Med Oncol (2020) 31:387–94. doi:10.1016/j.annonc.2019.11.010
51. Budczies, J, Kluck, K, Beck, S, Ourailidis, I, Allgauer, M, Menzel, M, et al. Homologous recombination deficiency is inversely correlated with microsatellite instability and identifies immunologically cold tumors in most cancer types. J Pathol Clin Res (2022) 8:371–82. doi:10.1002/cjp2.271
52. Pellegrino, B, Musolino, A, Llop-Guevara, A, Serra, V, De Silva, P, Hlavata, Z, et al. Homologous recombination repair deficiency and the immune response in breast cancer: A literature review. Transl Oncol (2020) 13:410–22. doi:10.1016/j.tranon.2019.10.010
53. Solinas, C, Marcoux, D, Garaud, S, Vitoria, JR, Van den Eynden, G, de Wind, A, et al. BRCA gene mutations do not shape the extent and organization of tumor infiltrating lymphocytes in triple negative breast cancer. Cancer Lett (2019) 450:88–97. doi:10.1016/j.canlet.2019.02.027
Keywords: immunotherapy, breast cancer, prognosis, immune microenvironment, NAD+ biosynthesis
Citation: Yang Y, Wang Z, He M, Diao L, Yu B and Li D (2023) NAD+ biosynthesis metabolism predicts prognosis and indicates immune microenvironment for breast cancer. Pathol. Oncol. Res. 29:1610956. doi: 10.3389/pore.2023.1610956
Received: 18 November 2022; Accepted: 28 February 2023;
Published: 17 March 2023.
Edited by:
Andrea Ladányi, National Institute of Oncology (NIO), HungaryCopyright © 2023 Yang, Wang, He, Diao, Yu and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dong Li, lidong.bprc@foxmail.com
†These authors have contributed equally to this work
