Skip to main content
Download PDF

Comprehensive characterization of B7 family members in breast cancer: B7-H5 switch reverses breast cancer from “immuno-cold” into “immuno-hot” status

Cancer Cell International volume 24, Article number: 205 (2024) Cite this article

Abstract

The members of the classic B7 family regulate the immune microenvironment of several malignant tumors. However, the potential relationship between the B7 family and the breast cancer (BrCa) tumor immune microenvironment has remained elusive. In the present study, we provide a comprehensive explanation of the expression, clinical significance, mutation, and immune cell infiltration of B7 family molecules in BrCa. First, we recruited 10 patients with BrCa surgery from the Wuxi Maternal and Child Health Hospital and performed single-cell RNA sequencing (scRNA-seq) analysis to investigate the distribution of B7 family members in multiple immune cell subsets. We focused on B7-2, B7-H3, and B7-H5 molecules of the B7 family and constructed tumor microarrays by self-recruiting patients to perform multiple immunohistochemical (mIHC) analyses and study tumor expression of B7-2, B7-H3, B7-H5 and CD8+ immune cell infiltration. B7-H5 displayed a strong correlation with CD8+ immune cell infiltration. In summary, B7-H5 provides a new perspective for the identification of immunothermal subtypes of BrCa and could function as a switch to reverse BrCa from an “immunologically cold” state to an “immunologically hot” state.

Graphical abstract

Introduction

Breast cancer (BrCa) ranks second in mortality among females worldwide, accounting for 31% [1]. Because the traditional treatments of BrCa are limited by several shortcomings, increasing attention has been paid to immunotherapy, which could result in long-term disease control in patients with metastatic BrCa [2]. Immune checkpoint inhibitor (ICI)-based immunotherapy is a cutting-edge method of treating cancer by hindering tumor recurrence [3]. In contrast to chemotherapy and radiation, which indirectly impair cell growth and survival, immunotherapy directly targets the disease by enhancing anticancer immune responses [4]. However, BrCa has not been considered a highly immunogenic malignancy since the 1980s and the overall response rate has remained disappointing when PD-1 or PD-L1-targeted treatments were used alone or in combination with traditional chemotherapy [5]. Consequently, efficient biomarkers are required to categorize patients with BrCa into the best immunotherapy groups and identify those who could benefit from immunotherapy.

Neoadjuvant treatment (NAT) has consistently improved clinical results in patients with BrCa such that a significant percentage of patients who were previously inoperable are now able to undergo surgical resection [6]. A gene expression-based evaluation is crucial for forecasting the course and prognosis of patients with BrCa although it is challenging to precisely estimate a patient’s treatment response to NAT because of the complexity and heterogeneity of BrCa malignancies [7]. Recent research has implicated TIME in controlling BrCa development and NAT effectiveness, and the interplay involves tumor cells, stromal cells, tumor-infiltrating lymphocytes (TILs), and cancer-associated fibroblasts [8]. Tumors’ immunological “immuno-hot and immuno-cold” states are separated according to the nature of TIME. In detail, “immune-hot” tumors consistently exhibit positive therapeutic effects to different treatments, distinguished by effective T cell infiltration and immune support duration. “Immune-cold” tumors are characterized by immunosuppressive duration and resistance to most treatments, including immunotherapy and cytotoxic chemotherapy [9, 10]. Thus, it is crucial to identify “immuno-hot” phenotypes in patients with BrCa because they will benefit more from immunotherapy and have a better prognosis.

The B7 family molecules, which include co-stimulatory and co-inhibitory molecules, intricately regulate immune responses to elicit protective immunity against cancer and infection as well as limit the overactivity of the immune system [11, 12]. B7-1, B7-2, B7-H1, B7-DC, B7-H2, B7-H3, B7-H4, B7-H5, B7-H6, and B7-H7 constitute the B7 family [13]. The common aliases of B7 family members and their respective ligands are shown in Table 1. The immunological responses mediated by T cells can be inhibited or stimulated via the interaction of B7 ligands with different receptors. These substances are involved in tumor evasion and significantly contribute to the activity, differentiation, effective functions, survival, and inhibition of T cell immune response [14, 15]. The signals produced by the B7 family can be manipulated to treat BrCa. However, further research is required to comprehensively understand the possible correlation between B7 family molecules and BrCa, particularly concerning the “immuno-hot” and “immuno-cold” aspects of these molecules and BrCa.

The importance of B7 molecules, particularly B7-H3, B7-2, and B7-H5, in regulating tumor immunity has been widely discussed in recent years. Recent data indicate that B7-H3 is connected to tumor metastasis, proliferation, and treatment resistance, which are associated with a bad prognosis for patients [16]. Furthermore, immunological cells expressing B7-H3 include neutrophils, macrophages, B cells, dendritic cells, myeloid suppressor cells (MDSCs), monocytes, and activated T cells. These findings suggest that B7-H3 regulates TIME [17]. Eltanbouly demonstrated the expression of B7-H5 in human cancer tissues, including those from colorectal, gastric, non-small cell lung, pancreatic, ovarian, prostate, hepatocellular, acute myeloid leukemia, and metastatic melanoma cancers [18]. However, little is known about the expression patterns of B7-H5, B7-H3, and B7-2 in BrCa. In the present study, we systematically characterized the expression patterns and immuno-biological correlations of B7 family members, particularly B7-2, B7-H3, and B7-H5 to investigate immune-related biomarkers and stratify patients into the best immunotherapy strategies. We found a substantial correlation between B7-H5 expression and CD8+ T cell infiltration. Consequently, we provide an original strategy to distinguish between the “immuno-hot” and “immune-cold” phenotypes of BrCa. Individuals exhibiting high levels of B7-H5 expression display an immuno-hot phenotype, which may benefit them from immunotherapy.

Table 1 Common aliases of B7 family members and their respective ligands
Full size table

Materials and methods

Single-cell RNA-seq data acquisition and analysis

Single-cell RNA-sequencing (scRNA-seq) was performed by Shanghai Genechem Co., Ltd. on an Illumina HiSeq XTen machine. Cell Ranger 3.0.2 was applied to perform sample demultiplexing, barcode processing, and generate gene count data for each cell. The R package “singleR” package was applied and the Human Primary Cell Atlas Data, Blueprint Encode Data, and Immune Cell Expression Data were used as reference data for auxiliary annotation. This was followed by applying the Cell Marker database and previous studies to identify marker genes for manual annotation of different clusters.

Acquisition and analysis of TCGA data

Normalized RNA-sequencing and clinical data of BrCa samples were downloaded from the TCGA dataset (https://portal.gdc.cancer.gov/) with clinical features, including Stage, T, N, M and survival outcomes. The association between B7 family and survival outcomes were analyzed via the GEPIA database (http://gepia.cancer-pku.cn/), including overall survival and disease-free survival. The genetic alteration data of B7 molecules was downloaded from the cBioPortal dataset (https://www.cbioportal.org/) including mutation rate, genomic alterations and lollipop plots. TIME analysis was performed using the “e1071”, “parallel”, and “preprocess Core” packages, including TIME scores, immune cell infiltrating and immune functions. The analysis for differential expression of B7 family was conducted via “limma” package while the differentially expressed genes with the criteria of fdr < 0.001 were furtherly analyzed for GO and GSEA functional enrichments through “cluster Profiler,” “org.Hs.eg.db,” “enrichplot,” and “ggplot2” packages. The risk classification of BrCa patients in TCGA database was conducted into two clusters according to B7-H5 and CD8 + immune cell infiltrating via the “Consensus Cluster Plus” package.

Clinical specimens

The BrCa tissue microarrays (TMA, Cat. T22-0792-149) were obtained from Outdo BioTech (Shanghai, China). The T22-0792-149 TMA contained 63 female BrCa tissues and 20 adjacent tissues, which were used for exploring B7-2, B7-H3, and B7-H5 expression and infiltration of CD8+ immune cell infiltration. All patients in the cohort were recruited from the Wuxi Maternal and Child Health Hospital and submitted for mIHC analysis. Detailed clinicopathological characteristics of the cohorts were provided by Freethinking Biotechnology Co., Ltd. (Nanjing, China). Ethical approval for the TMA study was granted by the Clinical Research Ethics Committee, Outdo BioTech (Shanghai, China). The Institutional Review Board of the Wuxi Maternal and Child Health Hospital approved the study’s ethical affairs, and each subject provided written consent (Wuxi Maternal and Child Health Hospital Ethical Review [2021-01-0927-28]).

MIHC

MIHC staining was performed on the TMAs of BrCa tissues according to the standardized procedure. Ethylenediaminetetraacetic acid (EDTA) was used for antigen retrieval, and the primary antibodies were incubated overnight at 4℃. The primary antibodies used were as follows: anti-B7-2 (1:8000 dilution, Cat. ab219648, clone: EPR20115; Abcam, Cambridge, UK), anti-B7-H5 (1:50 dilution, Cat. ab252438, clone: EPR23665–20, Abcam), anti-CD8 (ready-to-use, Cat. PA067, clone: 457F6F8, Abcarta, Suzhou, China) and anti-B7-H3 (1:200 dilution, Cat. ab134161, clone: EPNCIR122, Abcam).

Quantitative analysis of mIHC

MIHC staining was used to study the expression of each protein marker after panoramic scanning. The images were analyzed using the halo pathological image analysis software. The signal intensity of the channel was expressed by histochemical score (H-score). H-score ranges from 0 to 300 and is calculated as follows:

$$\eqalign{& H - score\, = \left[ {1 \times *\left( {\% {\rm{ }}cells{\rm{ }}1 + } \right) + {\rm{ }}2 \times *\left( {\% {\rm{ }}cells{\rm{ }}2 + } \right)} \right. \cr & \left. {\quad\quad\quad\quad\quad\quad+ {\rm{ }}3 \times *\left( {\% {\rm{ }}cells{\rm{ }}3 + } \right)} \right] \cr}$$

Each staining channel was used as a unit. The analysis results included statistics on the overall number of positive cells of each channel in each sample and the number of positive cells in the three grades of weak, moderate, and strong, and the corresponding positive cells in the three parts, namely, the nucleus, cytoplasm, and cytomembrane. Finally, we calculated the corresponding percentage in the whole cell and the H-score value of each channel.

Cell culture

BrCa cell lines SK-BR-3 (Cat. KG197), MDA-MB-231 (Cat. KG033), MCF-7 (Cat. KG031), BT-549 (Cat. KG413), and MCF-10 A (Cat. HH0220) were authenticated using short tandem repeat profiling and obtained from KeyGen (Nanjing, China) and Boster (Wuhan, China). SK-BR-3 cells were cultured in McCoy’s 5 A medium containing 10% fetal bovine serum (FBS) at 37℃ with 5% CO2. MCF-7 cells were cultured in Roswell Park Memorial Institute (RPMI)-1640 medium containing 10% FBS at 37 °C with 5% CO2, BT-549 cells were cultured in RPMI-1640 medium containing 10% FBS and 0.023 U/mL insulin at 37℃ with 5% CO2, and MDA-MB-231 and MCF-10 A cells were cultured in Dulbecco’s Modified Eagle Medium (DMEM) medium containing 10% FBS and 5% HS and 20 ng/mL EGF and 0.5 µg/mL hydrocortisone and 100 ng/mL cholera toxin and 10 µg/mL insulin and 1% P/S at 37℃ with 5% CO2. All assays were conducted with mycoplasma-free.

Real-time quantitative PCR (real-time qPCR)

Real-time quantitative polymerase chain reaction (qPCR) was performed using the Applied Biosystems QuantStudio 5 system (Applied Biosystems; CA, USA). The values were normalized to that of glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and reflected the mean of three separate trials. The 2-ΔΔCt method for quantification was used to assess the data. The primer sequences are shown below (“F” for “forward” and “R” for “reverse”): B7-2 F, 5’-CTGCTCATCTATACACGGTTACC-3’ and R, 5’-GGAAACGTCGTACAGTTCTGTG-3’; B7-H3 F, 5’-TCTGGGCATCCCAAGTTTTGAC-3’ and R, 5’-TCCGCCTTTTGATCTCCGATT-3’; B7-H5 F, 5’ -ACGCCGTATTCCCTGTATGTC-3’ and R, 5’ -TTGTAGAAGGTCACATCGTGC-3’; GAPDH F, 5’ -GGACCTGACCTGCCGTCTAG-3’ and R, 5’ -GTAGCCCAGGATGCCCTTGA-3’.

Western blotting

The cells were lysed using the radioimmunoprecipitation (RIPA) lysis buffer (Solarbio; Beijing, China), containing phosphatase and protease inhibitors. The protein concentration was measured using the bicinchoninic acid (BCA) protein detection kit (Vazyme; Nanjing, China). Equal amounts of a sample protein (25 µg) were loaded onto a sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) gel and then transferred to polyvinylidene fluoride (PVDF) membranes (Beyotime; Shanghai, China). The membranes were blocked for 2 h at normal temperature and incubated overnight at 4℃ with primary antibodies, including anti-B7-2 antibody (Abcam, ab239075, 1:1000), anti-B7-H3 antibody (Abcam, ab134161, 1:1000), anti-B7-H5 antibody (Abcam, ab300042, 1:1000) and anti-β-actin antibody (Abcam, ab8226, 1:1000). Afterward, the membranes were incubated with a suitable secondary antibody (Beyotime) at room temperature for 2 h. Finally, the protein bands were observed using an enhanced chemiluminescence (ECL) reagent (Vazyme).

Statistical analysis

All statistical analyses were performed using the SPSS V.24.0 (Chicago, Illinois, USA) or GraphPad Prism V.8 (La Jolla, California, USA) software. The correlation between clinical features, immune cell infiltration, and expression of B7-2, B7-H3, and B7-H5 in BrCa samples was studied using either χ2 or Fisher’s exact test. Three independent experiments were conducted and the experimental data were analyzed using Student’s t-test. For all statistical analyses, p ≤ 0.05 was considered significant.

Results

Expression and scRNA-seq distribution of B7 family members in BrCa

The members of the B7 family and their aliases are shown in Table 1. We initially used the TCGA database to compare the expression of B7 members in normal and BrCa samples. Compared with normal tissues, eight B7 family members demonstrated differential expression in BrCa, including B7-1, B7-2, B7-DC, B7-H1, B7-H3, B7-H4, B7-H5, and B7-H6 (Fig. 1A). In addition, we enrolled 10 BrCa patients from the Wuxi Maternal and Child Health Hospital for scRNA-seq analysis to investigate the distribution of different cell types in BrCa TIME, including B cells, basal cells, endothelial cells, fibroblasts, luminal cells, mast cells, myeloid cells, pericytes, plasma cells, and T/natural killer (NK) cells (Fig. 1B). Next, we explored the association between B7 family member expression and scRNA-seq distribution in BrCa TIME. The cell distribution of B7-H7 was almost invisible. B7-1 and B7-2 were mainly concentrated in myeloid and B cells. The distribution of B7-H1, B7-H2, and B7-H6 was similar, and they were distributed in cell subsets depicted in the diagram. B7-DC was mainly present in myeloid cells and fibroblasts. B7-H4 was mainly found in luminal cells, whereas B7-H3 was significantly enriched in several subgroups except for T/NK and B cells. B7-H5 was largely distributed in all involved cell subsets, except in luminal cells (Fig. 1C).

Fig. 1
figure 1

Expression and scRNA-seq distribution of B7 family members in BrCa (A) The expression of the B7 family in normal and BrCa tissues using the TCGA database. (B) The scRNA-seq analysis of 10 BrCa patients shows the distribution of different cell types in BrCa TIME. (C) The relationship between B7 family expression and scRNA-seq distribution in BrCa TIME.

Full size image

Association of B7 family with clinicopathological parameters and survival in BrCa

We next analyzed the correlation between B7 family members and clinicopathological parameters in BrCa, including TNM stage, T stage, N stage, and M stage (Fig. 2). According to the TNM stage, B7-H3 demonstrated significant differential expression among different TNM stages (p < 0.05). In addition, based on the T stage, B7-DC (p<0.05), B7-H4 (p<0.0001), and B7-H7 (p<0.01) demonstrated significant differential expression among different T stages. In addition, the expression of B7-H3 (p < 0.05), B7-H5 (p < 0.001), and B7-H6 (p < 0.01) was related to the N stage, whereas the expression of B7-H3 (p < 0.05) correlated with the M stage. Further, we investigated the relationship between the B7 family and survival in BrCa, including overall survival (OS) and disease-free survival (DFS). Thus, no significant association was present between B7 family members and patients with OS and DFS (See Fig. 3).

Fig. 2
figure 2

Association between B7 family and clinicopathological parameters and overall survival in BrCa. (A) Expression of 10 B7 family members among different TNM stages in BrCa. (B) Expression of 10 B7 family members among different T stages in BrCa. (C) Expression of 10 B7 family members among different N stages in BrCa. (D) Expression of 10 B7 family members among different M stages in BrCa

Full size image
Fig. 3
figure 3

Association between B7 family expression and survival in BrCa. (A) Association between B7 family expression and overall survival (OS) in BrCa. (B) Association between B7 family expression and disease-free survival (DFS) in BrCa

Full size image

Functional enrichment and signaling pathways of B7 family members in BrCa

The median of B7 family members was used to perform differential analysis between low-expressed and high-expressed groups, whereas differentially expressed genes of B7 family members were presented using the volcano plot (Fig. 4A). Next, we conducted the gene ontology (GO) enrichment analysis of B7 family members, including biological process (BP), cellular component (CC), and molecular function (MF) (Fig. 4B). The gene set enrichment analysis (GSEA) enrichment analysis of B7 family members was performed for potential signaling pathways (Fig. 4C).

B7-1 was enriched in immune response-activating cell surface receptors, external side of the plasma membrane, and protein serine or threonine kinase activity, whereas B7-1 was enriched in ribosomes, oxidative phosphorylation, JAK–STAT signaling pathway, T cell receptor signaling pathway, NK cell-mediated cytotoxicity, and Parkinson’s disease. Similarly, B7-2 was enriched in leukocyte migration, the external side of the plasma membrane, and nucleoside-triphosphatase regulator activity, whereas B7-2 was enriched in ribosomes, cytokine, asthma, viral myocarditis, glycosylphosphatidylinositol (GPI) anchor biosynthesis, and oxidative phosphorylation. B7-H1 was enriched in the regulation of lymphocyte activation, the external side of the plasma membrane, and antigen binding, whereas B7-H1 was enriched in the ribosomes, oxidative spliceosomes, Leishmania infection, viral myocarditis, and cytokine–cytokine receptor interaction. B7-DC was enriched in leukocyte migration, endosome membrane, and cell adhesion molecule (CAM) binding, whereas B7-DC was enriched in oxidative phosphorylation, ribosome, CAMs, viral myocarditis, Leishmania infection, and Parkinson’s disease. B7-H2 was enriched in phagocytosis, immunoglobulin complex, and antigen binding, whereas B7-H2 was enriched in intestinal immune network for IGA production, melanogenesis, cytokine, cytokine receptor interaction, endocytosis, neurotrophin signaling pathway, and cell cycle. B7-H3 was enriched in extracellular structure organization, collagen-containing extracellular matrix, and CAM binding, whereas B7-H3 was enriched in basal cell carcinoma, leukocyte, neuroactive ligand–receptor interaction, taste transduction, pyrimidine metabolism, and amyotrophic lateral sclerosis. B7-H4 was enriched in extracellular structure organization, cell–cell junction, and CAM binding, whereas B7-H4 was enriched in the proteasome, base excision repair, DNA replication, melanogenesis, hedgehog signaling pathway, and leukocyte transendothelial migration. B7-H5 was enriched in the regulation of lymphocyte activation, cell–substrate junction, and CAM binding, whereas B7-H5 was enriched in systemic lupus erythematosus, CAMs, cytokine–cytokine receptor interaction, GPI anchor biosynthesis, nucleotide excision repair, and mismatch repair. B7-H6 was enriched in proteasomal protein catabolic process, mitochondrial inner membrane, and small GTPase binding, whereas B7-H6 was enriched in ribosomes, Parkinson’s disease, oxidative phosphorylation, taste transduction, axon guidance, and oocyte meiosis. B7-H7 was enriched in the regulation of lymphocyte activation, the external side of the plasma membrane, and antigen binding, whereas B7-H7 was enriched in spliceosomes, RNA degradation, primary immunodeficiency, GPI anchor biosynthesis, intestinal immune network for IGA production, and autoimmune thyroid disease.

Fig. 4
figure 4

Functional enrichment and signaling pathways of B7 family members in BrCa. (A) The volcano plot depicts differentially expressed genes of B7 family members in BrCa. (B) The GO functional enrichment analysis of B7 family members in BrCa. (C) The GSEA signaling pathway analysis of B7 family members in BrCa

Full size image

Genomic alterations of B7 family members in BrCa

We used the cBioPortal database to investigate the relationship between B7 family members and genomic alterations in BrCa. With respect to the genomic alteration frequency, the mutation ratios of B7-H5 and B7-H6 were relatively higher (up to 3%), whereas the mutation rates of B7-1 and B7-2 were lower, less than 1% (Fig. 5A). Genomic alterations of B7 family members in different BrCa-related researches were presented, including structural variant data, mutation data, and CNA data (Fig. 5B). In addition, the lollipop plot of B7 family members in BrCa is depicted (Fig. 5C).

Fig. 5
figure 5

Genomic alterations of B7 family members in BrCa. (A) The mutation rate of B7 family members in BrCa. (B) The genomic alterations of B7 family members in different BrCa-related researches. (C) The lollipop plot of B7 family members in BrCa

Full size image

Correlation between B7 family members and TIME in BrCa

We next analyzed the relationship between the B7 family and TIME in BrCa, including TIME scores, immune cell infiltration, and immune functions. TIME score included estimate score, immune score, and stromal score, whereas the immune cell infiltration included 23 kinds of immune cells. The immune functions consisted of antigen-presenting cell (APC) co-inhibition, APC co-stimulation, CCR, checkpoint, cytolytic activity, HLA, promotion of inflammation, MHC class I, parainflammation, T cell co-stimulation, T cell co-inhibition, type I interferon (IFN) response, and type II IFN response. According to the TIME scores, 10 B7 family members exhibited significant differences in estimate scores and immune scores, whereas 8 B7 family members displayed significant differences in stromal scores, except for B7-H2 and B7-H6 (Fig. 6A-C). In addition, the relationship between B7 family members and immune cell infiltration, and immune functions were presented (Fig. 7A-B).

Fig. 6
figure 6

Correlation between B7 family members and TIME in BrCa. Relationship between B7 family members and TIME estimate score. (B) Relationship between B7 family members and TIME immune score. (C) Relationship between B7 family members and TIME stromal score

Full size image
Fig. 7
figure 7

Correlation between B7 family members and immune cell infiltration and immune functions. (A) Relationship between B7 family members and immune cell infiltration. (B) Relationship between B7 family members and immune function

Full size image

Expression of B7-2, B7-H3, and B7-H5 in BrCa

We recruited a clinical cohort from the Wuxi Maternal and Child Health Hospital. The preliminary analysis revealed the correlation between the B7 family and the clinical stage. B7-H3 and B7-H5, two molecules with clinical research potential, as well as B7-2, a B7 family member with little research in BrCa, were selected for immunohistochemical and mIHC analyses. Figure 8A depicts representative images of B7-2, B7-H3, and B7-H5 expression in mIHC-stained tumor tissue. Representative images stained with B7-2, B7-H3, and B7-H5 were divided into low expression group and high expression group for display. DAPI staining results were used as a reference. Statistical analysis of H-score staining results in BrCa microarray tumor (n = 63) and adjacent tumor (n = 20) tissues demonstrated that the expression of B7-2 and B7-H3 in BrCa tissues was significantly higher than in adjacent tumor tissues, whereas no statistical difference was found in B7-H5 between adjacent tumor tissues and BrCa tissues (Fig. 8B). In addition, we used multiple subtypes of BrCa cells to further verify, whereas we used normal breast cells MCF-10 A as controls. The results of western blotting and qRT-PCR were consistent, indicating that B7-2 and B7-H3 were highly expressed in BrCa, and the expression of B7-H5 was low in BrCa (Fig. 8C-D).

According to the clinical cohort in Wuxi Maternal and Child Health Hospital, the relationship between clinical features and expression of B7-2, B7-H3, and B7-H5 was analyzed using the chi-square test or Fisher’s test. The statistical results demonstrated that B7-2, B7-H3, and B7-H5 in BrCa tissues were not statistically correlated with clinical features (Table 2).

Table 2 Association between B7-2, B7-H3, B7-H5, and clinical features, infiltration of CD8+immune cells
Full size table
Fig. 8
figure 8

B7-2, B7-H3, and B7-H5 expression in BrCa. (A) Representative images of B7-2, B7-H3, and B7-H5 expression in BrCa tissues stained by mIHC. All representative images were divided into low expression group and high expression group. The DAPI staining results are listed as a reference. (B) The expression of B7-2, B7-H3, and B7-H5 in several BrCa cells was assessed by qRT-PCR. (C) The expression of B7-2, B7-H3, and B7-H5 in several BrCa cells was assessed by western blotting. Normal breast cell line MCF-10 A was used as a reference. Significance was calculated using Student’s t-test

Full size image

B7-H5 is significantly correlated with CD8+cell infiltration and has the potential to predict immunothermal tumors

CD8+ immune cells were divided into three layers with high, medium and low expression. We used the clinical cohort in Wuxi Maternal and Child Health Hospital to study the relationship between infiltration of CD8+ immune cells and expression of B7-2, B7-H3, and B7-H5 using chi-square test or Fisher’s test. We found that B7-H5 significantly correlated with CD8+ cells (p < 0.001), as shown in Table 2. We demonstrated the representative images of B7-H5 and CD8 high, middle, and low layers. As depicted in Fig. 9A, the chip diagram (40 ×) uses a white box to select typical mIHC images, and four single-channel display images of DAPI\B7-H5\CD8+ (200 ×) are conducted. Results demonstrated that CD8+ infiltration increased significantly with an increase in B7-H5 expression. We believe that B7-H5 is a potential molecule in predicting thermal tumors, that is, immune invasive tumors.

We performed a cluster analysis to divide the expression of B7H5 and CD8 into cluster 1 and cluster 2, in which cluster 2 referred to high expression of B7H5 and high CD8+ (Fig. 10A-C). The risk classification of patients with BrCa was performed based on the expression of B7-H5 and CD8+, which was further assessed in terms of TIME, degree of immune cell infiltration, and immune functions (Fig. 10D-F).

Fig. 9
figure 9

B7-H5 and CD8 + expression in BrCa. (A) Representative images uncovering B7-H5 and CD8+ expression in tumor tissues using mIHC staining. The images are divided into three levels: B7-H5 low and CD8+ low, B7-H5 mid and CD8+ mid, and B7-H5 high and CD8+ high. The DAPI staining results are used as a reference

Full size image
Fig. 10
figure 10

Association between B7-H5 expression and CD8+ immune cell infiltration in BrCa patients. (A-C) Risk classification of patients with BrCa according to B7-H5 and CD8+ expression. (D) TIME estimate score of risk classification in BrCa. (E) TIME immune score of risk classification in BrCa. (F) TIME stromal score of risk classification in BrCa. (G) Immune cell infiltration of risk classification in BrCa. (H) Immune functions of risk classification in BrCa

Full size image

Discussion

The B7 family has recently received considerable attention. Immunomodulatory B7 family proteins are significantly expressed in different malignancies and are known to regulate a disease’s advancement and unfavorable prognosis [12, 19]. The expression of B7 family members could be a useful predictive indicator for women with BrCa who are diagnosed and treated early. Increasing evidence suggests that B7 molecular expression dysregulation inhibits TIME and immune escape. In this work, we first studied the expression patterns of these B7 molecules in BrCa. Compared with normal tissues, eight B7 family members displayed differential expression in BrCa, including B7-1, B7-2, B7-DC, B7-H1, B7-H3, B7-H4, B7-H5, and B7-H6. Furthermore, the B7 family members are genetically stable and change with a low frequency in BrCa.

B7-H3 is a transmembrane immunoregulatory protein that is physically linked to B7 family proteins. It is a negative immunomodulatory protein, which is upregulated in different human malignant tumors and downregulated in normal tissues [20]. It is reported that B7-H3 knockdown increases the chemosensitivity of mantle cell lymphoma [21]. In addition, B7-H3 silencing increases drug-induced cytotoxicity and drug-induced apoptosis in acute monocytic leukemia (AML) cells via upregulating the caspase-3 (CASP-3) activity in vitro [22]. Artemether derivatives derived from Chinese herbs can treat a range of malignant tumors. Research has demonstrated that doxorubicin could suppress B7-H3 expression in neuroblastoma cell lines, whereas B7-H3 suppression further inhibited tumor development in cells treated with artemether [23]. The B7-H3-mediated treatment resistance in melanoma is associated with the activation of the P38–MAPK pathway [24]. Chemoresistance in ovarian cancer is attributed to B7-H3 expression, which stimulates the PI3K/AKT signaling pathway and upregulates Bcl-2 in protein levels [25, 26]. Thus, B7-H3 has emerged as a viable therapeutic target for cancer.

B7-2 functions as a ligand for CD28 and cytotoxic T lymphocyte antigen 4 (CTLA-4) [27]. B7-2 is expressed in dendritic cells, B/T cells, and APCs, whereas APCs have been shown to upregulate B7-2 expression under stimulatory circumstances [28]. CD28 has a high affinity for B7-2, whereas CTLA-4 has an even higher affinity for both B7-1 and B7-2 [29, 30]. A few studies have been conducted on the prognostic and clinical importance of B7-2 in malignant tumors. Soluble B7-2 (sB7-2) is present in the serum of healthy individuals and is present in resting monocytes, dendritic cells, and certain cancer cells [31, 32]. Elevated sB7-2 levels have been reported in different leukemias, including acute myeloid leukemia (AML), and B-chronic lymphocytic leukemia (B-CLL) [32,33,34]. Patients with myelodysplastic syndrome (MDS) exhibit greater levels of sB7-2. Compared to AML patients with normal B7-2 levels, those with higher levels of B7-2 had lower chances of complete response and shorter survival [35]. Furthermore, the anti-tumor response demonstrated by B7.2-IgG was CD8+ dependent.

The inclusion of B7-H5, a T-cell-activated V-domain immunoglobulin inhibitor, could be a viable target for immune oncology [36, 37]. B7-H5 is mostly expressed in hematopoietic tissue and tissues with a high concentration of white blood cells. In comparison to its decreased expression on CD4+ and CD8+ T lymphocytes, B7-H5 is overexpressed in bone marrow cells. B7-H5 expression is absent in B cells. However, it was highly expressed in plasma cells [18, 38, 39]. Because of its main homeostatic function, B7-H5 regulates the immune system. Consequently, B7-H5 and other immunomodulatory receptors are not the same [40]. B7-H5 is used as both a ligand and a receptor, having an identical inhibitory function [41]. When used either by itself or in conjunction with another immune checkpoint inhibitor, B7-H5 represents a potential novel target for cancer immunotherapy either by itself or in conjunction with another immune checkpoint inhibitor [38]. A recent study reported that the most effective cancer treatment could involve combining anti-B7-H5 and anti-PD-1 antibodies, as demonstrated by their effectiveness in a mouse tumor model [42].

Among the 10 B7 family members, we selected three critical B7 family members for further investigation in BrCa, including B7-2, B7-H3, and B7-H5. We constructed a TMA of 63 female BrCa patients and 20 paracarcinoma samples for mIHC analysis to investigate the relationship between the expression of B7-2, B7-H3, B7-H5, and infiltration of CD8+ immune cells. We found that B7-H5 displayed a positive correlation with CD8+ immune cell infiltrating in BrCa, providing a novel method for distinguishing between “immuno-hot” and “immuno-cold” subtypes in BrCa.

The major findings of this study could offer a strong scientific foundation for clinical diagnosis and therapy assessment. However, the study had several limitations. Initially, an internal cohort was not used to confirm the expression, correlation, and mutation of the B7 family in BrCa. Instead, the same findings were investigated using a public TCGA cohort. In addition, scRNA-seq data cannot discriminate between different cell types, and certain variations could be present in the expression at the mRNA and protein levels. Furthermore, although we suggested a unique BrCa subtype method, it has not yet been utilized in a clinical immunotherapy cohort.

Conclusion

In summary, we performed a thorough investigation of B7 family members and clinical features, OS, functional enrichments, genomic alterations, and TIME in BrCa. B7-H5 was found to be strongly associated with CD8+ immune cell infiltration in BrCa. Consequently, we provided an original strategy to distinguish between the “immuno-hot” and “immune-cold” phenotypes of BrCa: individuals exhibiting high levels of B7-H5 expression displayed an immuno-hot phenotype, which could benefit them from immunotherapy. Altogether, B7-H5 could function as a switch reversing BrCa from “immuno-cold” into “immuno-hot” status.

Data availability

No datasets were generated or analysed during the current study.

Abbreviations

BrCa:

Breast cancer

ScRNA-seq:

Single-cell RNA sequencing

TCGA:

The Cancer Genome Atlas

TIME:

Tumor immune microenvironment

TMA:

Tumor microarray

MIHC:

Multiple immunohistochemical

ICIs:

Immune checkpoint inhibitors

NAT:

Neoadjuvant treatment

TILs:

Tumor-infiltrating lymphocytes

MDSCs:

Myeloidal suppressor cells

OS:

Overall survival

DFS:

Disease-free survival

BP:

Biological process

CC:

Cellular component

MF:

Molecular function

H -score:

I Histochemical score

AML:

acute monocytic leukemia

CASP-3:

Caspase-3

CTLA-4:

Cytotoxic T lymphocyte antigen 4

APC:

Antigen presenting cells

B -CLL:

C B-chronic lymphocytic leukemia

MDS:

Myelodysplastic syndrome

References

  1. Siegel RL, Miller KD, Wagle NS, et al. Cancer statistics, 2023 [J]. Cancer J Clin. 2023;73(1):17–48.

    Article  Google Scholar 

  2. Harbeck N, Gnant M. Breast cancer [J]. Lancet (London England). 2017;389(10074):1134–50.

    Article  PubMed  Google Scholar 

  3. Abdul-Latif M, Townsend K, Dearman C, et al. Immunotherapy in gastrointestinal cancer: the current scenario and future perspectives [J]. Cancer Treat Rev. 2020;88:102030.

    Article  CAS  PubMed  Google Scholar 

  4. Ochsenbein AF. Principles of tumor immunosurveillance and implications for immunotherapy [J]. Cancer Gene Ther. 2002;9(12):1043–55.

    Article  CAS  PubMed  Google Scholar 

  5. Esteva FJ, Hubbard-Lucey VM, Tang J, et al. Immunotherapy and targeted therapy combinations in metastatic breast cancer [J]. Lancet Oncol. 2019;20(3):e175–86.

    Article  CAS  PubMed  Google Scholar 

  6. Sung H, Ferlay J, Siegel RL, et al. Global Cancer statistics 2020: GLOBOCAN estimates of incidence and Mortality Worldwide for 36 cancers in 185 countries [J]. Cancer J Clin. 2021;71(3):209–49.

    Article  Google Scholar 

  7. Quaresma M, Coleman MP, Rachet B. 40-year trends in an index of survival for all cancers combined and survival adjusted for age and sex for each cancer in England and Wales, 1971–2011: a population-based study [J]. Lancet (London England). 2015;385(9974):1206–18.

    Article  PubMed  Google Scholar 

  8. Singh HM, Bailey P, Hübschmann D, et al. Poly(ADP-ribose) polymerase inhibition in pancreatic cancer [J]. Genes Chromosomes Cancer. 2021;60(5):373–84.

    Article  CAS  PubMed  Google Scholar 

  9. Paz-Ares L, Dvorkin M, Chen Y, et al. Durvalumab plus platinum-etoposide versus platinum-etoposide in first-line treatment of extensive-stage small-cell lung cancer (CASPIAN): a randomised, controlled, open-label, phase 3 trial [J]. Lancet (London England). 2019;394(10212):1929–39.

    Article  CAS  PubMed  Google Scholar 

  10. Cristescu R, Mogg R, Ayers M, et al. Pan-tumor genomic biomarkers for PD-1 checkpoint blockade-based immunotherapy [J]. Volume 362. Science (New York, NY); 2018. 6411.

  11. Su M, Lin Y, Cui C, et al. ERMAP is a B7 family-related molecule that negatively regulates T cell and macrophage responses [J]. Cell Mol Immunol. 2021;18(8):1920–33.

    Article  CAS  PubMed  Google Scholar 

  12. Ni L, Dong C. New B7 Family checkpoints in Human cancers [J]. Mol Cancer Ther. 2017;16(7):1203–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Liu W, Zang X. Structures of Immune checkpoints: an overview on the CD28-B7 family [J]. Adv Exp Med Biol. 2019;1172:63–78.

    Article  CAS  PubMed  Google Scholar 

  14. Liu D, Jenkins RW, Sullivan RJ. Mechanisms of resistance to Immune Checkpoint blockade [J]. Am J Clin Dermatol. 2019;20(1):41–54.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Emens LA. Breast Cancer Immunotherapy: facts and hopes [J]. Clin cancer Research: Official J Am Association Cancer Res. 2018;24(3):511–20.

    Article  CAS  Google Scholar 

  16. Picarda E, Ohaegbulam KC, Zang X. Molecular pathways: Targeting B7-H3 (CD276) for Human Cancer immunotherapy [J]. Clin cancer Research: Official J Am Association Cancer Res. 2016;22(14):3425–31.

    Article  CAS  Google Scholar 

  17. Zhao B, Li H, Xia Y, et al. Immune checkpoint of B7-H3 in cancer: from immunology to clinical immunotherapy [J]. J Hematol Oncol. 2022;15(1):153.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Eltanbouly MA, Croteau W, Noelle RJ, et al. VISTA: a novel immunotherapy target for normalizing innate and adaptive immunity [J]. Semin Immunol. 2019;42:101308.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Flem-Karlsen K, Fodstad Ø, Nunes-Xavier CE. B7-H3 Immune checkpoint protein in Human Cancer [J]. Curr Med Chem. 2020;27(24):4062–86.

    Article  CAS  PubMed  Google Scholar 

  20. Getu AA, Tigabu A, Zhou M, et al. New frontiers in immune checkpoint B7-H3 (CD276) research and drug development [J]. Mol Cancer. 2023;22(1):43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Zhang W, Wang Y, Wang J, et al. B7-H3 silencing inhibits tumor progression of mantle cell lymphoma and enhances chemosensitivity [J]. Int J Oncol. 2015;46(6):2562–72.

    Article  CAS  PubMed  Google Scholar 

  22. Zhang W, Wang J, Wang Y, et al. B7-H3 silencing by RNAi inhibits tumor progression and enhances chemosensitivity in U937 cells [J]. Onco Targets Ther. 2015;8:1721–33.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Tan WQ, Chen G, Ye M et al. Artemether regulates Chemosensitivity to Doxorubicin via Regulation of B7-H3 in human neuroblastoma cells [J]. Medical science monitor: breastbe.

  24. Flem-Karlsen K, Tekle C, Øyjord T, et al. p38 MAPK activation through B7-H3-mediated DUSP10 repression promotes chemoresistance [J]. Sci Rep. 2019;9(1):5839.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Yang T, Xu F, Sheng Y, et al. A targeted proteomics approach to the quantitative analysis of ERK/Bcl-2-mediated anti-apoptosis and multi-drug resistance in breastbreast cancer [J]. Anal Bioanal Chem. 2016;408(26):7491–503.

    Article  CAS  PubMed  Google Scholar 

  26. Zhou L, Zhao Y. B7-H3 induces ovarian Cancer drugs Resistance through an PI3K/AKT/BCL-2 signaling pathway [J]. Cancer Manage Res. 2019;11:10205–14.

    Article  CAS  Google Scholar 

  27. Tang W, Zhao G. Small molecules targeting HIF-1α pathway for cancer therapy in recent years [J]. Volume 28. Bioorganic & medicinal chemistry; 2020. p. 115235. 2.

  28. Khan M, Arooj S, Wang H. Soluble B7-CD28 Family Inhibitory Immune Checkpoint proteins and Anti-cancer Immunotherapy [J]. Front Immunol. 2021;12:651634.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Engelhardt JJ, Sullivan TJ, Allison JP. CTLA-4 overexpression inhibits T cell responses through a CD28-B7-dependent mechanism [J]. J Immunol (Baltimore Md: 1950). 2006;177(2):1052–61.

    Article  CAS  Google Scholar 

  30. Buchbinder EI, Desai A. CTLA-4 and PD-1 pathways: similarities, differences, and implications of their inhibition [J]. Am J Clin Oncol. 2016;39(1):98–106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Jeannin P, Magistrelli G, Aubry JP, et al. Soluble CD86 is a costimulatory molecule for human T lymphocytes [J]. Immunity. 2000;13(3):303–12.

    Article  CAS  PubMed  Google Scholar 

  32. Hock BD, Patton WN, Budhia S, et al. Human plasma contains a soluble form of CD86 which is present at elevated levels in some leukaemia patients [J]. Leukemia. 2002;16(5):865–73.

    Article  CAS  PubMed  Google Scholar 

  33. Dong MP, Enomoto M, Thuy LTT, et al. Clinical significance of circulating soluble immune checkpoint proteins in sorafenib-treated patients with advanced hepatocellular carcinoma [J]. Sci Rep. 2020;10(1):3392.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Odagiri N, Hai H, Thuy LTT et al. Early change in the plasma levels of circulating Soluble Immune Checkpoint proteins in patients with Unresectable Hepatocellular Carcinoma Treated by Lenvatinib or Transcatheter arterial chemoembolization [J]. Cancers, 2020, 12(8).

  35. Hock BD, Mckenzie JL, Patton WN, et al. The clinical significance of soluble CD86 levels in patients with acute myeloid leukemia and myelodysplastic syndrome [J]. Cancer. 2003;98(8):1681–8.

    Article  CAS  PubMed  Google Scholar 

  36. Deng J, Le Mercier I, Kuta A, et al. A New VISTA on combination therapy for negative checkpoint regulator blockade [J]. J Immunother Cancer. 2016;4:86.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Huang X, Zhang X, Li E, et al. VISTA: an immune regulatory protein checking tumor and immune cells in cancer immunotherapy [J]. J Hematol Oncol. 2020;13(1):83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Le Mercier I, Chen W, Lines JL, et al. VISTA regulates the development of Protective Antitumor immunity [J]. Cancer Res. 2014;74(7):1933–44.

    Article  PubMed  Google Scholar 

  39. Borggrewe M, Grit C, Den Dunnen WF, A, et al. VISTA expression by microglia decreases during inflammation and is differentially regulated in CNS diseases [J]. Glia. 2018;66(12):2645–58.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Eltanbouly MA, Schaafsma E, Noelle RJ, et al. VISTA: coming of age as a multi-lineage immune checkpoint [J]. Clin Exp Immunol. 2020;200(2):120–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Nowak EC, Lines JL, Varn FS, et al. Immunoregulatory functions of VISTA [J]. Immunol Rev. 2017;276(1):66–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Liu J, Yuan Y, Chen W, et al. Immune-checkpoint proteins VISTA and PD-1 nonredundantly regulate murine T-cell responses [J]. Proc Natl Acad Sci USA. 2015;112(21):6682–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was founded by Major project of Wuxi Science and Technology Bureau (N20201006), Wuxi Double-Hundred Talent Fund Project (BJ2023075), Wuxi Health Commission Precision Medicine Project (J202106), Jiangsu Provincial Maternal and Child Health Research Project (F202034), Jiangsu Provincial Six Talent Peaks Project (YY-124), General Program of Wuxi Medical Center of Nanjing Medical University (WMCG202354), Doctoral Talent Fund of the Affiliated Wuxi People’s Hospital of Nanjing Medical University (BSRC202303).

Author information

Author notes

  1. Jiayu Liu, Cenzhu Wang, Ying Jiang and Yunxu Zhou contributed equally to this work.

Authors and Affiliations

  1. Department of Oncology, Wuxi Maternal and Child Health Care Hospital, Wuxi School of Medicine, Jiangnan University, No.48 Huaishu Road, Wuxi, Jiangsu, 214002, China

    Jiayu Liu, Ying Jiang, Yunxu Zhou & Yan Zhang

  2. Wuxi People’s Hospital, Wuxi Medical Center, The Affiliated Wuxi People’s Hospital of Nanjing Medical University, Nanjing Medical University, Jiangsu, 214000, China

    Cenzhu Wang

  3. Department of Oncology, The Affiliated Wuxi People’s Hospital of Nanjing Medical University, Jiangsu, 214023, China

    Cenzhu Wang

  4. Wuxi Maternal and Child Health Hospital, Nanjing Medical University, Jiangsu, 214000, China

    Lingyan Chen & Zhiwen Qian

  5. Wuxi Maternal and Child Health Hospital, Jiangnan University, Jiangsu, 214002, China

    Lu Liu, Danping Wu & Yan Zhang

Authors
  1. Jiayu Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Cenzhu Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Ying Jiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yunxu Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Lingyan Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Zhiwen Qian
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Lu Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Danping Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Yan Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

J.L.: data curation, methodology, formal analysis, and writing. C.W.: investigation, data curation, methodology, and writing. Y.J.: resources, data curation, and methodology. Y.Z.: formal analysis, and writing . L.C.: investigation and methodology. Z.Q.: methodology, and formal analysis. L.L.: investigation and methodology. D.W.: investigation and methodology. Y.Z.: conceptualization, resources, supervision and funding acquisition. All authors approved the final version of the manuscript.

Corresponding author

Correspondence to Yan Zhang.

Ethics declarations

Ethical approval

Ethical approval for the use of TMAs was granted by the Clinical Research Ethics Committee in Outdo Biotech. This study was reviewed and approved by the Wuxi Maternal and Child Health Hospital Ethical Review Committee, with the approval number 2021-01-0927-28. All experiments were performed in accordance with the Declaration of Helsinki, and informed consent was obtained from all subjects.

Consent for publication

All authors agreed to publish the manuscript.

Competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, J., Wang, C., Jiang, Y. et al. Comprehensive characterization of B7 family members in breast cancer: B7-H5 switch reverses breast cancer from “immuno-cold” into “immuno-hot” status. Cancer Cell Int 24, 205 (2024). https://doi.org/10.1186/s12935-024-03392-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12935-024-03392-4

Keywords

Cancer Cell International

ISSN: 1475-2867

Contact us