==== Front Front Immunol Front Immunol Front. Immunol. Frontiers in Immunology 1664-3224 Frontiers Media S.A. 10.3389/fimmu.2020.558757 Immunology Original Research Molecular and Clinical Characterization of PD-1 in Breast Cancer Using Large-Scale Transcriptome Data Liu Qiang † Cheng Ran † Kong Xiangyi † Wang Zhongzhao * Fang Yi * Wang Jing * Department of Breast Surgical Oncology, National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China Edited by: Ali A. Zarrin, TRex Bio, United States Reviewed by: Ye Li, University of Texas MD Anderson Cancer Center, United States; Bipulendu Jena, Independent Researcher, San Diego, CA, United States *Correspondence: Jing Wang, wangjing@cicams.ac.cn; Yi Fang, fangyi@cicams.ac.cn; Zhongzhao Wang, wangzhongzhao@cicams.ac.cnThis article was submitted to Cancer Immunity and Immunotherapy, a section of the journal Frontiers in Immunology †These authors have contributed equally to this work 17 11 2020 2020 11 55875704 5 2020 16 10 2020 Copyright © 2020 Liu, Cheng, Kong, Wang, Fang and Wang2020Liu, Cheng, Kong, Wang, Fang and WangThis 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.Despite the impressive impact of PD-1 (programmed cell death protein 1)-targeted cancer immunotherapy, a great part of patients with cancer fail to respond. PD-1 impact on immune cells in addition to T cells, and the synergistic role of PD-1 with other immune modulators remain largely unknown. To fill this gap, we systematically investigated PD-1-related transcriptome data and relevant clinical information derived from TCGA (The Cancer Genome Atlas) and METABRIC (Molecular Taxonomy of Breast Cancer International Consortium), which involved a total of 2,994 breast cancer patients. Our results revealed the relationship among PD-1 and major molecular and clinical characteristics in breast cancer. More importantly, we depicted the association landscape between PD-1 and other immune cell populations. Gene ontology analyses and gene set variation analysis (GSVA) of genes correlated with PD-1 revealed that PD-1 was mainly involved in immune responses and inflammatory activities. We also elucidated the association of PD-1 with other immune modulators in pan-cancer level, especially the potential synergistic relationship between PD-1 and other immune checkpoints members in breast cancer. In short, the expression level of PD-1 was bound up with breast cancer malignancy, which could be used as a potential biomarker; PD-1 might manipulate the anti-tumor immune response by impacting not just T cells, and this might vary among different tumor types. Furthermore, PD-1 might synergize with other immune checkpoint members to modulate the immune microenvironment in breast cancer. breast cancercancer immunotherapyPD-1immune responseinflammatory activityNational Natural Science Foundation of China10.13039/50110000180981872160 ==== Body Introduction Breast cancer is the most frequently diagnosed women cancer all over the world (1). As development of comprehensive treatments, including resection by surgery, chemotherapy, radiotherapy, endocrine therapy, and targeted therapy, many early stage breast cancers can be effectively controlled (2). However, for advanced breast cancers, the existing standard treatments have limited effectiveness because of the aggressiveness of tumors, resistance to treatments, recurrence, and metastasis during or after treatments (3). In the past decade, researches suggested that tumor cells’ immunologic escape of and aberrant human immune surveillance play essential roles in the carcinogenesis, progression, and metastasis of cancers (4), and studies focused on anticancer immune responses have achieved marked success of many malignant tumors in preclinical and clinical trials (5–9), the PD-1 (programmed cell death protein 1) and PD-L1 (programmed cell death-ligand 1) axis has been identified as one of the most encouraging findings in cancer immunotherapy (10). PD-1 is a receptor that is expressed on the surface of activated T cells, and the PD-L1 and PD-L2 are ligands of PD-1 that are expressed on the surface of antigen-presenting cells (11). The PD-1 and PD-L1 interaction can result in T cells inactivation and ensure that the immune system can be activated at the appropriate time, thus minimize the possibility of chronic autoimmune inflammation. However, tumor cells or non-transformed cells in the tumor microenvironment overexpressing PD-L1, leading to generate an adaptive immune resistance in response to endogenous immune anti-tumor activity (12). The antibody of PD-1/PD-L1 can block the immune escape mediated by PD-1 and PD-L1 interaction and has been approved by the FDA in a fast speed (13). In the first clinical trials of breast cancer, the inhibitors of PD-1/PD-L1 showed promising activity (14, 15). However, many problems still need to be resolved, including the lack of clearly available data in breast cancer, notably regarding PD-1 expression and its prognostic value, and the application of PD-1/PD-L1 inhibitors in combination with other immune checkpoint inhibitors. Despite the impressive impact of PD-1/PD-L1-targeted cancer immunotherapy, a large proportion of cancer patients fail to respond (16). Furthermore, although the combination of PD-1/PD-L1 blockade with complementary checkpoint inhibitors has achieved some success for some malignant tumors in the preclinical and clinical trials, the impact of PD-1 on immune cells in addition to T cells and the synergistic role of PD-1 with other immune modulators remain mostly unknown (17). In the present study, we systematically investigated the PD-1-related transcriptome profile and revealed its potential role in inducing immune responses and inflammatory activities as well as its potential relationship with immune modulators. Methods Data Collection Transcriptome data from TCGA (The Cancer Genome Atlas) were downloaded by GDCRNA tools (access date: February 01, 2020) (18). The edgeR (19) and limma packages (20) available from the Bioconductor project (21) offer a well-developed suite of statistical methods for dealing with this question for RNA-seq data. Raw count data were normalized using TMM implemented in edgeR, and then were transformed by voom in limma, only genes with cpm > 1 in more than half of the samples were kept. Sieved TCGA breast cancer clinical data were kindly provided by Dr. Hai Hu and Dr. Jianfang Liu (Chan Soon-Shiong Institute of Molecular Medicine, Windber). Human epidermal growth factor receptor 2 (HER2) status was recalled using DNA copy number for cases without immunohistochemistry or fluorescence in situ hybridization (FISH) status. Standardized survival data of TCGA cohort were downloaded from TCGA-CDR (TCGA Pan-Cancer Clinical Data Resource) (22). The METABRIC (Molecular Taxonomy of Breast Cancer International Consortium) dataset (23) containing 1904 tumor cases was downloaded from the cBioPortal database (http://www.cbioportal.org/) (access date: Feb 01, 2019). A total of 2,994 samples with full clinical characteristics and transcriptome data were used to perform the following data exploration. The detailed clinical characteristics of breast cancer patients from TCGA and METABRIC are listed in Tables S1 and S2, respectively). Bioinformatics Analysis Gene ontology analyses of the genes that correlated with PD-1 were performed using clusterProfiler package (24). Immunologically related genes were collected from The ImmPort (Immunology Database and Analysis Portal) database (https://www.immport.org/home) (25). The absolute abundance of immune cell populations was estimated using Microenvironment Cell Populations-counter algorithm (26). GSVA (Gene Set Variation Analysis) (27) was used to calculate scores of gene sets that correlated with immune functions and inflammatory activities (28). Association of PD-1 and other immune modulators in pan-cancer were depicted through the database of TISIDB (29), which is an integrated repository portal for tumor-immune system interactions. The study summary diagram is shown in Figure 1. Figure 1 Summary diagram of the present study. Statistical Analysis Spearman correlation method was used to estimate the correlations between continuous variables. Student t-test, one-way ANOVA, or Pearson’s Chi-squared test were used to determine any differences in variables between groups. R language was used to perform all statistical tests. The prognostic value of PD-1 was evaluated through Cox proportional hazards model analysis. Several packages including ggplot2 pheatmap, pROC (30), circlize (31), and corrgram (32) were used to perform other statistical calculations and graphical work (33), and P < 0.05 was considered to have the statistically significant difference. Results PD-1 Expression Pattern in Breast Cancer To characterize the relationship between PD-1 expression and molecular and clinical features in breast cancer, individuals were dichotomized into high and low groups based on the expression of PD-1 using median cut. We found that PD-1 was asscociated with patient age, American Joint Committee on Cancer (AJCC) stage, tumor grade, estrogen receptor (ER) status, progesterone receptor (PR) status and HER2 status (Tables 1 and 2). Next, we further detected that the PD-1 expression was upregulated in tumor tissues, the ER-negative group (ER-) and PR-negative group (PR-) in both TCGA and METABRIC datasets, while upregulation of PD-1 in the HER2-positive group was only observed in the METABRIC dataset (Figures 2A–F). Meanwhile, we observed that PD-1 expression was upregulated in the molecular subtypes such as basal-like and HER2-enriched when compared with luminal A, while no difference was found between luminal A and luminal B subtypes. While no claudin-low subtype was found in the TCGA cohort, the other four subtypes were consistently in both the TCGA and METABRIC cohorts (Figures 2G, H). Subsequently, we also found that PD-1 expression was enriched in higher grade tumors in the METABRIC dataset (Figure 2I). PD-1 overexpression was also oberved in the triple-negative breast cancer (TNBC) subtype and could be as a predictor for TNBC subtype in both TCGA (AUC = 0.671) and METABRIC (AUC = 0.672) databases (Figures 3A–D). One potential limitation of this result was that this was a primary result based on only one gene. Future studies focusing on investigating robust biomarkers for TNBC might consider combining PD-1 and other biomarkers together, and comparing the predictive capacity between each other. In summary, these findings suggest that PD-1 expression is enriched in higher malignant breast cancer and might be a potential biomaker in TNBC. Table 1 Association between PD-1 mRNA expression and clinicopathologic characteristics in TCGA cohort. Expression Total (n = 1090) PD-1 high (n = 545) PD-1 low (n = 545) P-value Age (years)  >=55 517 (47.4%) 277 (50.8%) 240 (44.0%) 0.025  <55 573 (52.6%) 268 (49.2%) 305 (56.0%) T stage  T1 279 (25.6%) 130 (23.9%) 149 (27.3%) 0.028  T2 631 (57.9%) 330 (60.6%) 301 (55.2%)  T3 137 (12.6%) 72 (13.2%) 65 (11.9%)  T4 40 (3.7%) 13 (2.4%) 27 (5.0%) Unknown 3 (0.3%) 0 (0%) 3 (0.6%) N stage  N0 514 (47.2%) 257 (47.2%) 257 (47.2%) 0.065  N1 360 (33.0%) 177 (32.5%) 183 (33.6%)  N2 120 (11.0%) 57 (10.5%) 63 (11.6%)  N3 76 (7.0%) 48 (8.8%) 28 (5.1%)  Unknown 20 (1.8%) 6 (1.1%) 14 (2.6%) M stage  M0 907 (83.2%) 452 (82.9%) 455 (83.5%) 0.061  M1 22 (2.0%) 6 (1.1%) 16 (2.9%)  Unknown 161 (14.8%) 87 (16.0%) 74 (13.6%) AJCC stage  I 181 (16.6%) 82 (15.0%) 99 (18.2%) 0.097  II 621 (57.0%) 319 (58.5%) 302 (55.4%)  III 250 (22.9%) 131 (24.0%) 119 (21.8%)  IV 20 (1.8%) 5 (0.9%) 15 (2.8%)  Unknown 18 (1.7%) 8 (1.5%) 10 (1.8%) ER status  Negative 236 (21.7%) 155 (28.4%) 81 (14.9%) <0.001  Positive 803 (73.7%) 372 (68.3%) 431 (79.1%)  Unknown 51 (4.7%) 18 (3.3%) 33 (6.1%) PR status  Negative 343 (31.5%) 203 (37.2%) 140 (25.7%) <0.001  Positive 694 (63.7%) 323 (59.3%) 371 (68.1%)  Unknown 53 (4.9%) 19 (3.5%) 34 (6.2%) HER2 status  Negative 895 (82.1%) 449 (82.4%) 446 (81.8%) 0.032  Positive 168 (15.4%) 89 (16.3%) 79 (14.5%)  Unknown 27 (2.5%) 7 (1.3%) 20 (3.7%) Table 2 Association between PD-1 mRNA expression and clinicopathologic characteristics in METABRIC cohort. Expression Total (n = 1904) PD-1 high (n = 952) PD-1 low (n = 952) P-value Age (years)  >=55 952 (50.0%) 522 (54.8%) 430 (45.2%) <0.001  <55 952 (50.0%) 430 (45.2%) 522 (54.8%) Tumor size  >=2cm 592 (31.1%) 293 (30.8%) 299 (31.4%) 0.079  <2cm 1292 (67.9%) 644 (67.6%) 648 (68.1%)  Unknown 20 (1.1%) 15 (1.6%) 5 (0.5%) AJCC stage  0 4 (0.2%) 3 (0.3%) 1 (0.1%) 0.039  I 475 (24.9%) 218 (22.9%) 257 (27.0%)  II 800 (42.0%) 413 (43.4%) 387 (40.7%)  III 115 (6.0%) 68 (7.1%) 47 (4.9%)  IV 9 (0.5%) 2 (0.2%) 7 (0.7%)  Unknown 501 (26.3%) 248 (26.1%) 253 (26.6%) Tumor grade  I 165 (8.7%) 50 (5.3%) 115 (12.1%) <0.001  II 740 (38.9%) 311 (32.7%) 429 (45.1%)  III 927 (48.7%) 551 (57.9%) 376 (39.5%)  Unknown 72 (3.8%) 40 (4.2%) 32 (3.4%) ER status  Negative 445 (23.4%) 309 (32.5%) 136 (14.3%) <0.001  Positive 1459 (76.6%) 643 (67.5%) 816 (85.7%) PR status  Negative 895 (47.0%) 526 (55.3%) 369 (38.8%) <0.001  Positive 1009 (53.0%) 426 (44.7%) 583 (61.2%) HER2 status  Negative 1668 (87.6%) 800 (84.0%) 868 (91.2%) <0.001  Positive 236 (12.4%) 152 (16.0%) 84 (8.8%) Figure 2 The expression of PD-1 in different ER, PR, and HER2 status (A–F), molecular subtypes (G, H), and grades (I) in TCGA or METABRIC cohort. (*: P < 0.05, ****: P < 0.0001, ns: no significant difference). Figure 3 PD-1 serves as a potential biomarker. PD-1 expression pattern between TNBC and none-TNBC tissues in TCGA and METABRIC (A, B). ROC curves predict PD-1 as a biomarker of TNBC (C, D). (****: P < 0.0001). PD-1 Was Bound Up With Immune Functions in Breast Cancer To further explore PD-1 related biological processes in breast cancer, a total of 1008 genes and 449 genes, seived from TCGA and METABRIC datasets, respectively, have strongly correlation with PD-1 according to Spearman correlation analysis (|R| > 0.4 and P < 0.05). Subsequently, GO (gene ontology) enrichment analyses were performed to investigate PD-1’s potential biological functions. We found that PD-1-related genes were mainly involved in immune-related pathways and inflamatory pathways (Figure 4A), including T cell regulation-related biological processes and leukocyte regulation-related pathways. Importantly, these results were also validated in METABRIC datasets (Figure 4B). Figure 4 PD-1 is closely related to immune functions in breast cancer. Gene ontology analysis shows that PD-1 is mainly involved in immune response and inflammatory response in TCGA and METABRIC databases (A, B). PD-1 Related Immune Response To further investigate PD-1-related immune functions in breast cancer, 4723 immunologically related genes were retrieved from the ImmPort (25). Genes that were most correlated with PD-1 (Spearman |R| > 0.4, P < 0.05) were seived to depict the expression pattern of these immunologically related genes in breast cancer. Subsequently, we found that 508 and 248 immune-related genes in TCGA and METABRIC datasets, respectively, were correlated with PD-1 postively, while only 6 and no immunologically related genes, respectively, were correlated with PD-1 negatively (Figures S1A, B). The Relationship Between PD-1 Expression and Immune Cell Populations Previously, Tao Jiang et al. reported PD-1 is a positively correlated with T cells, monocytic lineage cells, and myeloid dendritic cells, but not with cytotoxic lymphocytes, natural killer (NK) cells, or B lineage cells in diffuse glioma (34). In breast cancer, to further clarify the immune manipulative functions of PD-1, the Microenvironment Cell Populations-counter algorithm was used to calculate the absolute abundance of immune cell populations (26). The abundance pattern of these cell populations in breast cancer is depicted in Figures 5A, B. We observed that PD-1 expression was strongly correlated with immune cell population scores of T cells, CD8+ T cells, cytotoxic lymphocytes, NK cells, B lineage cells, monocytic lineage cells, and myeloid dendritic cells, but not neutrophils, endothelial cells, or fibroblasts (Figures 5C, D). These findings suggest that PD-1 may not just be involved in regualting T cell immunity, other immune cell immunity might also be involved. Furthermore, the immune regulatory pattern of PD-1 may be varied in different tumors. Figure 5 The relationship between PD-1 expression and immune cell populations in TCGA and METABRIC databases (A–D). Subtypes denotes breast cancer molecular subtypes including Basal, basal-like; Her2, Her2-enriched; LumA, luminal A; LumB, luminal B. Relationship of the Expression of PD-1 and Immune Modulators in Pan-Cancer To explore the synergistic role of PD-1 and other immune modulators in pan-cancer, we systematically analyzed the correlations between PD-1 expression and three types of immune modulators described in the previous study conducted by Charoentong et al. (35). Interestingly, we found a similar correlation pattern between immune modulators and PD-1 across 30 tumor types in which the majority of immunoinhibitors and immunostimulators were correlated with PD-1 positively (Figures S2 and S3), while only a small number of immunoinhibitors and immunostimulators were negatively correlated with PD-1. More interestingly, we observed that PD-1 was positively correlated with almost all MHC molecules across 30 tumor types (Figure S4). These results provide a landscape perspective regarding the correlation of PD-1 with immune modulators, and we could compare this correlation pattern of PD-1 with immune modulators in various tumors. The general PD-1 correlation pattern trends were similar in various tumors. These observations suggest that PD-1 might synergize with other immune modulators in manipulating the anti-tumor immune response. PD-1 Is Synergistic With Other Immune Checkpoint Members in Breast Cancer-Induced Immune Response We estimated the association of PD-1 with other immune checkpoint members to further explore the synergistic role of PD-1 in breast cancer-induced immune responses (Figures 6A–D). The detailed R and p-values of correlations between PD-1 and other checkpoint members are listed in Tables S3. Interestingly, we found that PD-1 was strongly correlated with similar checkpoint members including PD-L2, CD274 (PD-L1) CTLA4, IDO1, and LAG3, as well as other checkpoint members including BTLA, ICOS, CD27, CD40, and CD48. These findings revealed that PD-1 might manipulate anti-tumor immune responses through co-regulation with the above mentioned immune checkpoint molecules, thereby lending support to using combination cancer immunotherapy targeting these molecules in future studies. Figure 6 PD-1 expression is correlated with immune checkpoint members in TCGA and METABRIC databases (A–D). Asscociation Between PD-1 and Specific Cell Immune Responses The specific immune regulatory role of PD-1 in breast cancer remains largely unknown. To further explore the association between PD-1 and specific immune responses, GSVA of gene ontology biological pathways was performed. Consistent with the above results, we found that PD-1 was strongly correlated with both T and B cell immunity (Figures 7A, B). PD-1 was positively correlated with T cell proliferation, T cell differentiation, T cell activation and T cell receptor signaling pathways. Meanwhile, we also observed that PD-1 was positively correlated with B cell activation. Figure 7 PD-1-related cell immunity and inflammatory activities in breast cancer. The relationship between PD-1 and cell immunity in TCGA and METABRIC datasets (A, B). The relationship between PD-1 and inflammatory activities in TCGA and METABRIC datasets (C, D). GO:0030217: T cell differentiation; GO:0042098:T cell proliferation; GO:0042102: positive regulation of T cell proliferation; GO:0042113: B cell activation; GO:0042129: regulation of T cell proliferation; GO:0050852: T cell receptor signaling pathway; GO:0050870: positive regulation of T cell activation. The pie denotes the correlation coefficient of PD-1 and GO term. Relationship Between PD-1 and Inflammatory Activities To further understand the role of PD-1 in mediating inflammatory activities, 104 genes were derived from seven clusters and defined as metagenes using the GSVA algorithm (28), implicating different types of inflammation and immune functions. Detailed information of these metagenes is shown in Table S4. We found that LAG3 was positively correlated with MHC-I, MHC-II, LCK, STAT1, HCK, and interferon, but not IgG (Figures 7C, D). Among these seven clusters, PD-1 showed the strongest correlation with MHC-II and LCK in both TCGA and METABRIC databases. These results further suggested that PD-1 not only correlated with T cell immunity but also with other immune cells. In summary, these findings indicated that PD-1 has important immune and inflammatory functions in breast cancer. Prognostic Value of PD-1 in Breast Cancer To explore the influence on breast cancer survival, we evaluated the prognostic value of PD-1 in both TCGA (n = 1090) and METABRIC (n = 1994) cohorts. Interestingly, our data indicated that PD-1 was an independent prognostic factor in breast cancer on the basis of TCGA cohort multivariate analysis after adjusting for patient age, AJCC stage, ER, PR, and HER2 status (Figure 8A). In the METABRIC cohort, the PD-1 expression level was also an independent prognostic indicator for breast cancer after adjusting for tumor grade, AJCC stage, ER, PR, and HER2 status (Figure 8B). Despite the fact that PD-1 was upregulated in higher malignant tumors, our results suggest that PD-1 is predictive of good prognosis in breast cancer patients. Figure 8 Prognostic value of PD-1 in breast cancer using TCGA and METABRIC databases (A and B). Discussion The development of immune therapies for solid tumors have promoted clinical advances (36), PD-1/PD-L1 are critical biological suppressors of cytotoxic immune reactions, and the PD-L1 expression is one of the major immunologic escape mechanisms in tumors (37). However, despite the impressive impact of PD-1/PD-L1-targeted cancer immunotherapy, a large proportion of cancer patients fail to respond. To understand why this occurs, we need to have a better understanding of the molecular regulatory mechanisms of PD-1/PD-L1. Whether PD-1/PL-L1 predicts prognosis in breast cancer patients is currently being debated, and large-scale investigations are still required to further confirm the specific relationship between PD-1/PD-L1 and prognosis of breast cancer (38). PD-1 was found to be significantly associated with better DFS and OS in Ren et al.’s study of only 195 TNBC patients (39). The present study used large-scale transcriptome data to provide additional strong evidence to support the correlation between high PD-1 expression and good prognosis in breast cancer patients. We estimated the association of PD-1 with other checkpoint proteins and found that PD-1 was strongly correlated with similar checkpoint members. Some of these results were consistent with previous studies, while some results have never been reported before. For instance, previous studies reported that LAG-3 and PD-1 were co-expressed on tumor-infiltrating lymphocytes, and blockade of both pathways had synergistic effects on anti-tumor CD8+ T cell stimulation and response (40, 41). In the present study, one interesting finding was that the immune inhibitors and immune stimulators were positively and concomitantly correlated with PD-1 expression. For primary and secondary resistance to immunotherapy, the etiologies are multifaceted, tumor intrinsic factors, and the complex interplay between cancer and its microenvironment all should to be taken into consideration (42). One possible explanation for immunotherpy resistance might be that there are complex interactions between PD-1 expression, immune inhibitors, and immune stimulators. Hence, the result that immune inhibitors and immune stimulators are positively correlated at the same time in the same sets of patients is reflective of a complex tumor microenvironment. More importantly, these results also suggest that when using immunotherapy, both activation of stimulatory pathways and blockade of inhibitory checkpoints can occur and should therefore be taken into consideration. It is generally agreed that PD-1 mainly inhibits the activation and immunologic function of T cells (43). We found that PD-1 was not only strongly correlated with T cell immunity, but also with B cell immunity. Tumor-infiltrating B cells with distinct phenotypes and functions, which might play specific roles in the anti-tumor responses (44). However, there is still a lack of direct evidence to support B cells have the immunosuppressive role in human cancers (45). Further study is needed to explore the possible manipulation by B cells-mediated immune suppression through the overexpression of PD-1. Moreover, analysis of the relationship between PD-1 and inflammatory activity, further suggested that PD-1 was not only correlated with T cell immunity but also with other immune cells. These findings indicated that PD-1 plays essential immune and inflammatory functions in breast cancer. In our study, we focused on characterizing the role of PD-1; however, future studies should address whether there is any concomitant over-expression of PD-L1 on tumors, on immune cell populations, or on any of these immune cells in parallel with overexpression of PD-1. This would be helpful in gaining a deeper understanding of the association of PD1/PDL1 with related immune cell populations. In summary, the expression of PD-1 was closely associated with the malignancy and might be as a potential biomaker in breast cancer, especially for TNBC. PD-1 might manipulate the anti-tumor immune response by impacting multiple immune cells, and this could vary with different tumors. Furthermore, PD-1 might synergize with other immune checkpoint members to modulate the immune microenvironment in breast cancer, which could be applied in the development of new targeted drugs for immunotherapy. This study had limitations. Because detailed treatment information was not available in the TCGA database, treatment effects were potential confounders that should be considered when available and adjusted for appropriately. In certain cases, a proxy for standard of care, such as age, treating hospital, and year of diagnosis can alleviate some of the bias when treatment is unknown. A future study using an alterative source of data should address this problem. Data Availability Statement Publicly available datasets were analyzed in this study. This data can be found here—TCGA: https://portal.gdc.cancer.gov/; METABRIC: http://www.cbioportal.org. Author Contributions JW: Conception and design and study supervision. YF and ZW: Acquisition of data and development of methodology. QL and RC: Analysis and interpretation of data and writing the original manuscript. XK: Verification and revision of the manuscript. All authors contributed to the article and approved the submitted version. Funding This work was supported by the Natural Science Foundation of China (No. 81872160); the Beijing Municipal Natural Science Foundation (Key Project) (No. 7191009); the Beijing Municipal Natural Science Foundation (No. 7204293); the Special Research Fund for Central Universities, Peking Union Medical College (No. 3332019053); the Beijing Hope Run Special Fund of Cancer Foundation of China (No. LC2019B03); the Beijing Hope Run Special Fund of Cancer Foundation of China (No. LC2019L07); the Golden Bridge Project Seed Fund of Beijing Association for Science and Technology; the PhD. Innovation Fund of Cancer Hospital, Chinese Academy of Medical Sciences (No. C2019-1051-09). 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 would like to thank Professor Hai Hu and Doctor Jianfang Liu for their kindly help providing us high-quality TCGA breast cancer clinical data. We wish to thank TCGA project organizers as well as all study participants. The first author is particularly indebted to Ping He for help during the study. Supplementary Material The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2020.558757/full#supplementary-material Supplementary Figure 1 PD-1-related immune responses. Most of immune-related genes are positively correlated with PD-1 expression in TCGA and METABRIC databases, while a small number of genes are negatively associated (A and B), Subtypes denotes breast cancer molecular subtypes including Basal, basal-like; Her2, Her2-enriched; LumA, luminal A; LumB, luminal B. Click here for additional data file. Supplementary Figure 2 PD-1 expression is correlated with immunoinhibitors in pan-cancer. Click here for additional data file. Supplementary Figure 3 PD-1 expression is correlated with immunostimulators in pan-cancer. Click here for additional data file. Supplementary Figure 4 PD-1 expression is correlated with MHC molecules in pan-cancer. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. ==== Refs References 1 DeSantis CE Ma J Gaudet MM Newman LA Miller KD Goding Sauer A Breast cancer statistics, 2019 . CA Cancer J Clin (2019 ) 69 :438–51. 10.3322/caac.21583 2 Moo TA Sanford R Dang C Morrow M Overview of Breast Cancer Therapy . PET Clin (2018 ) 13 :339–54. 10.1016/j.cpet.2018.02.006 3 Gerber B Freund M Reimer T Recurrent breast cancer: treatment strategies for maintaining and prolonging good quality of life . Dtsch Arztebl Int (2010 ) 107 :85 –91 . 10.3238/arztebl.2010.0085 20204119 4 Ghahremanloo A Soltani A Modaresi SMS Hashemy SI Recent advances in the clinical development of immune checkpoint blockade therapy . Cell Oncol (Dordr) (2019 ) 42 :609–26. 10.1007/s13402-019-00456-w 5 Wolchok JD Chiarion-Sileni V Gonzalez R Rutkowski P Grob JJ Cowey CL Overall Survival with Combined Nivolumab and Ipilimumab in Advanced Melanoma . N Engl J Med (2017 ) 377 :1345–56. 10.1056/NEJMoa1709684 6 Reck M Rodriguez-Abreu D Robinson AG Hui R Csoszi T Fulop A Pembrolizumab versus Chemotherapy for PD-L1-Positive Non-Small-Cell Lung Cancer . N Engl J Med (2016 ) 375 :1823–33. 10.1056/NEJMoa1606774 7 Ferris RL Blumenschein G Jr.Fayette J Guigay J Colevas AD Licitra L Nivolumab for Recurrent Squamous-Cell Carcinoma of the Head and Neck . N Engl J Med (2016 ) 375 :1856–67. 10.1056/NEJMoa1602252 8 Motzer RJ Escudier B McDermott DF George S Hammers HJ Srinivas S Nivolumab versus Everolimus in Advanced Renal-Cell Carcinoma . N Engl J Med (2015 ) 373 :1803–13. 10.1056/NEJMoa1510665 9 Brahmer JR Tykodi SS Chow LQ Hwu WJ Topalian SL Hwu P Safety and activity of anti-PD-L1 antibody in patients with advanced cancer . N Engl J Med (2012 ) 366 :2455–65. 10.1056/NEJMoa1200694 10 Yasunaga M Antibody therapeutics and immunoregulation in cancer and autoimmune disease . Semin Cancer Biol (2020 ) 64 :1 –12 . 10.1016/j.semcancer.2019.06.001 31181267 11 Zhao Y Harrison DL Song Y Ji J Huang J Hui E Antigen-Presenting Cell-Intrinsic PD-1 Neutralizes PD-L1 in cis to Attenuate PD-1 Signaling in T Cells . Cell Rep (2018 ) 24 :379 –390 e6 . 10.1016/j.celrep.2018.06.054 29996099 12 Pardoll DM The blockade of immune checkpoints in cancer immunotherapy . Nat Rev Cancer (2012 ) 12 :252–64. 10.1038/nrc3239 13 O’Donnell JS Long GV Scolyer RA Teng MW Smyth MJ Resistance to PD1/PDL1 checkpoint inhibition . Cancer Treat Rev (2017 ) 52 :71 –81 . 10.1016/j.ctrv.2016.11.007 27951441 14 Bedognetti D Maccalli C Al Bader SB Marincola FM Seliger B Checkpoint inhibitors and their application in breast cancer . Breast Care (2016 ) 11 :108–15. 10.1159/000445335 15 Hartkopf AD Taran F-A Wallwiener M Walter CB Krämer B Grischke E-M PD-1 and PD-L1 immune checkpoint blockade to treat breast cancer . Breast Care (2016 ) 11 :385–90. 10.1159/000453569 16 Bersanelli M Buti S From targeting the tumor to targeting the immune system: Transversal challenges in oncology with the inhibition of the PD-1/PD-L1 axis . World J Clin Oncol (2017 ) 8 :37 –53 . 10.5306/wjco.v8.i1.37 28246584 17 Ott PA Hodi FS Kaufman HL Wigginton JM Wolchok JD Combination immunotherapy: a road map . J Immunother Cancer (2017 ) 5 :16 . 10.1186/s40425-017-0218-5 28239469 18 Li R Qu H Wang S Wei J Zhang L Ma R GDCRNATools: an R/Bioconductor package for integrative analysis of lncRNA, miRNA and mRNA data in GDC . Bioinf (Oxford England) (2018 ) 34 :2515–7. 10.1093/bioinformatics/bty124 19 Robinson MD McCarthy DJ Smyth GK edgeR: a Bioconductor package for differential expression analysis of digital gene expression data . Bioinf (Oxford England) (2010 ) 26 :139–40. 10.1093/bioinformatics/btp616 20 Ritchie ME Phipson B Wu D Hu Y Law CW Shi W limma powers differential expression analyses for RNA-sequencing and microarray studies . Nucleic Acids Res (2015 ) 43 :e47 . 10.1093/nar/gkv007 25605792 21 Huber W Carey VJ Gentleman R Anders S Carlson M Carvalho BS Orchestrating high-throughput genomic analysis with Bioconductor . Nat Methods (2015 ) 12 :115–21. 10.1038/nmeth.3252 22 Liu J Lichtenberg T Hoadley KA Poisson LM Lazar AJ Cherniack AD An integrated TCGA pan-cancer clinical data resource to drive high-quality survival outcome analytics . Cell (2018 ) 173 :400 –16.e11 . 10.1158/1538-7445.AM2018-3287 29625055 23 Curtis C Shah SP Chin S-F Turashvili G Rueda OM Dunning MJ The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups . Nature (2012 ) 486 :346–52. 10.1038/nature10983 24 Yu G Wang L-G Han Y He Q-Y clusterProfiler: an R package for comparing biological themes among gene clusters . Omics: J Integr Biol (2012 ) 16 :284–7. 10.1089/omi.2011.0118 25 Bhattacharya S Andorf S Gomes L Dunn P Schaefer H Pontius J ImmPort: disseminating data to the public for the future of immunology . Immunol Res (2014 ) 58 :234–9. 10.1007/s12026-014-8516-1 26 Becht E Giraldo NA Lacroix L Buttard B Elarouci N Petitprez F Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression . Genome Biol (2016 ) 17 :218 . 10.1186/s13059-016-1113-y 27765066 27 Hänzelmann S Castelo R Guinney J GSVA: gene set variation analysis for microarray and RNA-seq data . BMC Bioinf (2013 ) 14 :7 . 10.1186/1471-2105-14-7 28 Rody A Holtrich U Pusztai L Liedtke C Gaetje R Ruckhaeberle E T-cell metagene predicts a favorable prognosis in estrogen receptor-negative and HER2-positive breast cancers . Breast Cancer Res: BCR (2009 ) 11 :R15 . 10.1186/bcr2234 19272155 29 Ru B Wong CN Tong Y Zhong JY Zhong SSW Wu WC TISIDB: an integrated repository portal for tumor–immune system interactions . Bioinf (Oxford England) (2019 ) 35 :4200–2. 10.1093/bioinformatics/btz210 30 Robin X Turck N Hainard A Tiberti N Lisacek F Sanchez J-C pROC: an open-source package for R and S+ to analyze and compare ROC curves . BMC Bioinf (2011 ) 12 :77 . 10.1186/1471-2105-12-77 31 Gu Z Gu L Eils R Schlesner M Brors B circlize implements and enhances circular visualization in R . Bioinf (Oxford England) (2014 ) 30 :2811–2. 10.1093/bioinformatics/btu393 32 Friendly M Corrgrams: Exploratory displays for correlation matrices . Am Statistician (2002 ) 56 :316–24. 10.1198/000313002533 33 Wickham H Wiley Interdiscip Rev Comput Stat. (2011 ) 3 :180–5. 10.1002/wics.147 34 Liu S Wang Z Wang Y Fan X Zhang C Ma W PD-1 related transcriptome profile and clinical outcome in diffuse gliomas . Oncoimmunology (2018 ) 7 :e1382792. 10.1080/2162402X.2017.1382792 29308304 35 Charoentong P Finotello F Angelova M Mayer C Efremova M Rieder D Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade . Cell Rep (2017 ) 18 :248–62. 10.1016/j.celrep.2016.12.019 36 Zhang H Chen J Current status and future directions of cancer immunotherapy . J Cancer (2018 ) 9 :1773–81. 10.7150/jca.24577 37 Schutz F Stefanovic S Mayer L von Au A Domschke C Sohn C PD-1/PD-L1 Pathway in Breast Cancer . Oncol Res Treat (2017 ) 40 :294–7. 10.1159/000464353 38 Jiang C Cao S Li N Jiang L Sun T PD-1 and PD-L1 correlated gene expression profiles and their association with clinical outcomes of breast cancer . Cancer Cell Int (2019 ) 19 :233 . 10.1186/s12935-019-0955-2 31516390 39 Ren X Wu H Lu J Zhang Y Luo Y Xu Q PD1 protein expression in tumor infiltrated lymphocytes rather than PDL1 in tumor cells predicts survival in triple-negative breast cancer . Cancer Biol Ther (2018 ) 19 :373–80. 10.1080/15384047.2018.1423919 40 Woo SR Turnis ME Goldberg MV Bankoti J Selby M Nirschl CJ Immune inhibitory molecules LAG-3 and PD-1 synergistically regulate T-cell function to promote tumoral immune escape . Cancer Res (2012 ) 72 :917–27. 10.1158/0008-5472.CAN-11-1620 41 Lichtenegger FS Rothe M Schnorfeil FM Deiser K Krupka C Augsberger C Targeting LAG-3 and PD-1 to Enhance T Cell Activation by Antigen-Presenting Cells . Front Immunol (2018 ) 9 :385. 10.3389/fimmu.2018.00385 29535740 42 Murciano-Goroff YR Warner AB Wolchok JD The future of cancer immunotherapy: microenvironment-targeting combinations . Cell Res (2020 ) 30 (6 ):507–19. 10.1038/s41422-020-0337-2 43 Mizuno R Sugiura D Shimizu K Maruhashi T Watada M Okazaki IM PD-1 Primarily Targets TCR Signal in the Inhibition of Functional T Cell Activation . Front Immunol (2019 ) 10 :630. 10.3389/fimmu.2019.00630 31001256 44 Zhang Y Gallastegui N Rosenblatt JD Regulatory B cells in anti-tumor immunity . Int Immunol (2015 ) 27 :521–30. 10.1093/intimm/dxv034 45 Ren Z Peng H Fu YX PD-1 Shapes B Cells as Evildoers in the Tumor Microenvironment . Cancer Discovery (2016 ) 6 :477–8. 10.1158/2159-8290.CD-16-0307