Literature DB >> 29147627

STAT1 gene deficient mice develop accelerated breast cancer growth and metastasis which is reduced by IL-17 blockade.

Sanjay Varikuti1, Steve Oghumu2, Mohamad Elbaz1,3, Greta Volpedo1, Dinesh K Ahirwar1, Pablo C Alarcon1, Rachel H Sperling1, Ellen Moretti1, Marissa S Pioso1, Jennifer Kimble1, Mohd W Nasser1,4, Ramesh K Ganju1, Cesar Terrazas1, Abhay R Satoskar1.   

Abstract

Signal transducer and activator of transcription 1 (STAT1) mediates interferon gamma signaling which activates the expression of various genes related to apoptosis, inflammation, cell cycle and angiogenesis. Several experimental and clinical studies have investigated the role of STAT1 in primary tumor growth in breast cancer; however, its role in tumor metastasis remains to be determined. To determine the role of STAT1 in breast cancer metastasis, we analyzed growth and metastasis in WT or STAT1-/- mice orthotopically implanted with metastatic 4T1.2 cells. Primary tumor development was faster in STAT1-/- mice and these mice developed significantly bigger primary tumors and displayed more lung metastasis compared with WT counterparts. STAT1-/- mice showed elevated Ly6G+CD11b+ granulocytic MDSC infiltration in their primary tumors and spleens with concomitant upregulation of Mmp9 and Cxcl1 expression in tumors compared with WT counterparts. Blockade of IL-17A in primary tumor-bearing STAT1-/- mice suppressed accumulation of Ly6G+CD11b+ cells and markedly reduced lung metastasis. These data show that STAT1 is an important suppressor of primary breast tumor growth and metastasis. Importantly, we found anti-IL-17 treatment can rescue STAT1 deficient animals from developing exacerbated metastasis to the lungs which could be important for immunotherapies for immunocompromised breast cancer patients.

Entities:  

Keywords:  IL-17; breast cancer; metastasis; neutrophil

Year:  2017        PMID: 29147627      PMCID: PMC5674966          DOI: 10.1080/2162402X.2017.1361088

Source DB:  PubMed          Journal:  Oncoimmunology        ISSN: 2162-4011            Impact factor:   8.110


Chemokine (C-X-C motif) ligand 1; Vegf Receptor tyrosine-protein kinase erbB-2 Granulocyte-Colony Stimulating Factor Myeloid Derived Suppressor Cells Matrix metallopeptidase 9 Programmed Death Ligand-1 Tumor Necrosis factor- Related Apoptosis-Inducing Ligand vesicular endothelial growth factor

Introduction

Signal transducers and activators of transcription (STATs) are a family of molecules that mediate intracellular signaling triggered mainly by cytokines. STAT1 is critical for the expression of several genes in response to interferons and has recently been involved in mediating signaling by regulatory cytokines such as IL-35 and IL-27. STAT1 is important for induction of the transcription factor T-bet and development of an efficient Th1 response. In turn, Th1 cells produce IFN-γ which activates macrophages to produce nitric oxide and mediates tumoricidal activity. In tumors, STAT1 activation by IFNs affects the expression of genes related to antiproliferative and apoptotic functions such as downregulation of c-myc, induction of caspases and surface death ligands such as TRAIL. Additionally, STAT1 has been shown to be important for immunosurvelliance by enhancing MHCI expression that is required for efficient antigen presentation to cytolytic lymphocytes. Interestingly, tumor intrinsic-STAT1 can repress angiogenesis by downregulating the expression of matrix metallo-proteinases (MMPs). Despite these findings, the role of STAT1 in cancer progression is still controversial. For example STAT1 expression was categorized as an oncogene in endometrial cancer by upregulating MYC and PDL1 expression. Conversely, absence of STAT1 abrogated the suppressor activity of tumor associated macrophages (TAMs) and myeloid derived suppressor cells (MDSCs) on T cells. Also of note, MDSCs are important suppressors of the anti-tumor immune response and deletion of this population has been shown to result in tumor reduction by enhancing the infiltration of cytotoxic T cells. Previous studies investigating the role of STAT1 in primary tumor development and growth of breast cancer have reported conflicting findings. While studies using STAT1−/− mice have found that STAT1 plays a critical role in preventing mammary tumor development and restricting tumor growth, Hix et al recently reported that upregulation of STAT1 expression in tumors correlated with disease progression to invasive carcinoma in both humans and mice. This study also found that overexpression of STAT1 in breast cancer cells promotes tumor growth in mice by mediating the accumulation of MDSCs into the tumor, indicating a role for tumor cell-intrinsic STAT1 in promoting homing of MDSCs into the tumor. These conflicting results reported in the above studies could be due to differences in the experimental models or to the fact that different roles of STAT1 may depend on its expression in particular cell types. Nevertheless, to our knowledge, the role of STAT1 in breast cancer metastasis has never been evaluated. In this study we examined the role of STAT1 in metastasis of breast cancer using STAT1−/− mice and an orthotopic mouse model of breast cancer metastasis to the lungs. Our results demonstrate that STAT1 plays a critical role in controlling primary tumor growth and inhibiting metastasis in breast cancer. Additionally, our studies reveal that the effect of host STAT1 deficiency on tumor metastasis to the lung can be reversed upon IL-17-blockade treatment which correlates with decreased expression of pro-angiogenic factors in the primary tumor.

Results

STAT1 gene deficient mice develop significantly larger primary tumors compared with WT mice

Following implantation of a highly metastatic 4T1.2 cells into mammary fat pads, STAT1−/− Balb/c female mice developed faster and significantly larger primary tumors compared with WT mice. Large primary tumors were evident at earlier time points in STAT1−/− mice compared with primary tumor masses in WT mice beginning at day 14 and thereafter (Fig. 1a). Furthermore, weights of primary tumors harvested from the STAT1−/− group at the end of the experiment were significantly higher compared with primary tumors from WT mice (Fig. 1b). These data indicate that STAT1 plays a critical role in limiting primary tumor growth in breast cancer.
Figure 1.

STAT1 deficient mice develop faster tumor growth and enhanced metastasis to the lungs. WT or STAT1−/− mice were injected with 4T1.2 cells in the mammary tissue. (a) Primary breast tumor measurements in WT and STAT1−/− mice during the course of the experiment. (b) Postmortem analysis of primary breast tumor weight in WT and STAT1−/− mice. (c) Representative images and H&E sections of lungs isolated from WT and STAT1−/− primary tumor bearing mice. (d) Counts of nodules per/lung in WT and STAT1−/− mice. Data from 2 independent experiments with an n-value of 5–6 mice per group. *p = 0.05, **p = 0.01.

STAT1 deficient mice develop faster tumor growth and enhanced metastasis to the lungs. WT or STAT1−/− mice were injected with 4T1.2 cells in the mammary tissue. (a) Primary breast tumor measurements in WT and STAT1−/− mice during the course of the experiment. (b) Postmortem analysis of primary breast tumor weight in WT and STAT1−/− mice. (c) Representative images and H&E sections of lungs isolated from WT and STAT1−/− primary tumor bearing mice. (d) Counts of nodules per/lung in WT and STAT1−/− mice. Data from 2 independent experiments with an n-value of 5–6 mice per group. *p = 0.05, **p = 0.01.

Absence of STAT1 leads to enhanced tumor metastasis to the lungs

The main cause of death among breast cancer patients is the metastasis of the primary tumor to distant sites; most commonly the bone, lung and liver. Despite the accumulated experimental data regarding the role of STAT1 in breast cancer, to our knowledge no study has shown the impact of host STAT1 deficiency on breast cancer metastasis. Thus we further analyzed the specific role of STAT1 on spontaneous metastasis of breast cancer cells. To this end we analyzed the tumor metastatic nodules in the lung. At day 25, we found lungs from STAT1−/− mice contained significantly more metastatic nodules compared with WT mice (Fig. 1c and d), indicating that STAT1 also plays an important role in restricting metastasis of breast cancer cells.

STAT1−/− mice show increased granulocytic MDSCs accumulation into tumors and spleens during breast cancer

Myeloid derived suppressor cells (MDSCs) are known to play an important role in tumor growth by inhibiting the activation of tumoricidal cytotoxic and T helper cells. MDSCs were first characterized as cells expressing Gr-1 and CD11b on their membrane. However, anti Gr-1 antibody recognizes the epitopes of both Ly6C and Ly6G molecules. More recently, MDSCs have been divided into 2 different subsets, one monocytic subset which expresses Ly6Chi Ly6GCD11b+ and another granulocytic subset which expresses Ly6G+ Ly6Cmed CD11b+ 16. Due to the importance of MDSCs in cancer progression, we evaluated the accumulation of MDSCs in the spleen and tumor by flow cytometry. The myeloid compartment in the spleens of WT and STAT1−/− mice was found to be expanded in breast cancer bearing mice compared with non-cancer bearing mice. However, the expansion of myeloid cells identified by CD11b expression was 3-fold higher in tumor-bearing STAT1−/− mice compared with WT mice (Fig. 2a and c). Next, we investigated surface markers associated with MDSCs. As shown in Fig. 2b, the CD11b+ gated cells were further subdivided by Ly6C and Ly6G expression. WT and STAT1−/− mice presented enhanced accumulation of Ly6ChiLy6GCD11b+ and Ly6G+Ly6CmedCD11b+ cells in the spleens compared with naïve mice (Fig. 2b). However, the frequency of both myeloid subsets was higher in STAT1−/− mice compared with WT mice (Fig. 2d and e). Interestingly Ly6G+Ly6CmedCD11b+ cells were the predominant subset in WT and STAT1−/− mice (Fig. 2b). Thus, absence of STAT1 results in increased accumulation of myeloid cells that resemble the monocytic-like and granulocytic-like MDSCs. The presence of MDSCs usually correlates with diminished T cell response upon re-stimulation. We confirmed the identity of Ly6G+Ly6CmedCD11b+ cells as MDSCs by their ability to suppress T cell proliferation using an in vitro suppression assay. Ly6G+Ly6CmedCD11b+ cells from tumor bearing WT and STAT1−/− mice were able to suppress responder T cell proliferation compared with control CFSE labeled T cells (Fig. 2f). Interestingly, although STAT1−/− mice showed enhanced accumulation of MDSCs in the spleen, following in vitro stimulation with anti-CD3, splenocytes from tumor bearing STAT1−/− mice showed increased proliferation compared with similarly activated spleen cells from tumor bearing WT mice (Fig. 2g). We also observed increased frequency of CD11b+Ly6C−/Ly6G− cells in spleens of tumor bearing STAT1−/− mice compared with WT counterparts (Fig. 2b). Our analyses show that these cells are a diverse population, which are predominantly F4/80+, CD11c-, NK1.1 (data not shown). Similar to the spleens, tumors of STAT1−/− mice also showed enhanced accumulation of Ly6G+Ly6CmedCD11b+ cells compared with the tumors of WT mice (Fig. 2h and i). PDL1 expression on myeloid cells has been shown to be a mechanism of immunoregulation on breast cancer. Also, PDL1 expression is dependent on STAT1 expression in some circumstances. We found that Ly6ChiLy6GCD11b+ and Ly6G+Ly6CmedCD11b+ MDSCs isolated from the spleens of WT and STAT1−/− mice expressed comparable levels of PDL1. On the other hand, Ly6ChiLy6GCD11b+ and Ly6G+Ly6CmedCD11b+ cells in the primary tumors of STAT1−/− mice expressed significantly less PDL1 compared with their WT counterparts (Supplementary Fig. 1). Together, these data show enhanced MDSCs accumulation in absence of STAT1; however the proliferation response of splenic T cells was augmented in these mice. Interestingly, PDL1 expression was dependent on STAT1 only in MDSCs which accumulated in the tumor.
Figure 2.

Enhanced recruitment of Ly6G+ cells in STAT1−/− tumor bearing mice. Flow cytometry analysis of splenocytes from WT or STAT1−/− mice. (a) Myeloid cells were identified by flow cytometry by labeling splenocytes with anti CD11b. (b) Myeloid cells gated in Fig. 2a were analyzed for Ly6G and Ly6C markers. (c) Frequencies of CD11b+ cell in the spleens of WT or STAT1−/− naïve or primary tumor bearing mice. (d) Frequencies of Ly6C+CD11b+ cells gated as in Fig. 2b. (e) Frequencies of Ly6G+CD11b+ cells in spleens of WT or STAT1−/− mice gated as in Fig. 2b. (f) Proliferation of CFSE labeled T cells incubated with either sorted WT or STAT1KO CD11b+ Ly6G+ cells. (g) Proliferation of splenocytes from WT or STAT1−/− tumor bearing mice re-stimulated with anti CD3 was evaluated by alamar blue reduction after 72h. (h) Tumors were harvested; myeloid cells were gated with CD11b and analyzed for Ly6G and Ly6C markers. (i) Frequencies of CD11b+Ly6G+Ly6Cmed cells in tumors of WT and STAT1−/− mice. *p = 0.05, **p = 0.01, ***p = 0.0001.

Enhanced recruitment of Ly6G+ cells in STAT1−/− tumor bearing mice. Flow cytometry analysis of splenocytes from WT or STAT1−/− mice. (a) Myeloid cells were identified by flow cytometry by labeling splenocytes with anti CD11b. (b) Myeloid cells gated in Fig. 2a were analyzed for Ly6G and Ly6C markers. (c) Frequencies of CD11b+ cell in the spleens of WT or STAT1−/− naïve or primary tumor bearing mice. (d) Frequencies of Ly6C+CD11b+ cells gated as in Fig. 2b. (e) Frequencies of Ly6G+CD11b+ cells in spleens of WT or STAT1−/− mice gated as in Fig. 2b. (f) Proliferation of CFSE labeled T cells incubated with either sorted WT or STAT1KO CD11b+ Ly6G+ cells. (g) Proliferation of splenocytes from WT or STAT1−/− tumor bearing mice re-stimulated with anti CD3 was evaluated by alamar blue reduction after 72h. (h) Tumors were harvested; myeloid cells were gated with CD11b and analyzed for Ly6G and Ly6C markers. (i) Frequencies of CD11b+Ly6G+Ly6Cmed cells in tumors of WT and STAT1−/− mice. *p = 0.05, **p = 0.01, ***p = 0.0001.

IL-17 neutralization reduces tumor metastasis to the lung in STAT1−/− mice

Ly6G+CD11b+ cell accumulation has been shown to be largely dependent on IL-17 production in several diseases. Recent reports have shown that IL-17 mediated accumulation of neutrophils has an important role in promoting breast cancer growth and metastasis in an experimental mode l26,27. To determine whether IL-17 was involved in the observed accumulation of Ly6G+Ly6CmedCD11b+ cells in the spleens and tumors, as well as in promoting primary tumor growth and/or metastasis in STAT1−/− mice, we first examined IL-17 levels in tumor bearing WT and STAT1−/− mice. Our results indicate that STAT1−/− mice expressed slightly higher Il-17 transcripts in their tumors compared with the tumors of WT mice, although not statistically significant (Fig. 3a). IL-17 cytokine production in stimulated splenocytes of tumor bearing STAT1−/− mice was also slightly higher than in WT mice, although this difference was not statistically significant (Fig. 3b). Interestingly, STAT1−/− mice intrinsically upregulate IL-17 in T cells stimulated with CD3 in vitro, which is enhanced in the presence of 4T1.2 breast cancer cell supernatants (Supplementary Fig. 2a, b). These studies corroborate previous findings that STAT1 deficiency promotes Th17 differentiation in vitro.
Figure 3.

Enhanced metastasis in STAT1−/− mice is controlled by anti-IL-17 treatment. WT and STAT1−/− female mice were injected with 4T1.2 cells in the mammary tissue. After 25 d groups of STAT1−/− mice received either isotype control or anti-IL-17 neutralizing antibody every other day for the rest of the experiment. (a) Gene expression analysis of IL-17 in primary tumors from WT and STAT1−/− mice as determined by RT PCR. (b) IL-17 cytokine production from splenocytes of tumor bearing WT and STAT1−/− mice stimulated with CD3. (c) Progression of primary breast tumor volume in WT, STAT1−/− and STAT1−/− treated with anti-IL-17. * and ** represent significant differences between WT and STAT1−/− mice (d) Primary breast tumor weights at the end of the experiment. (e) Number of metastatic tumor nodules in the lungs of the different groups. (f) Representative pictures of lungs from different experimental groups. (g) H&E staining of lung sections. Representative data from one of 2 independent experiments with an n-value of 5–7 mice per group. *p = 0.05, **p = 0.01, ***p = 0.0001.

Enhanced metastasis in STAT1−/− mice is controlled by anti-IL-17 treatment. WT and STAT1−/− female mice were injected with 4T1.2 cells in the mammary tissue. After 25 d groups of STAT1−/− mice received either isotype control or anti-IL-17 neutralizing antibody every other day for the rest of the experiment. (a) Gene expression analysis of IL-17 in primary tumors from WT and STAT1−/− mice as determined by RT PCR. (b) IL-17 cytokine production from splenocytes of tumor bearing WT and STAT1−/− mice stimulated with CD3. (c) Progression of primary breast tumor volume in WT, STAT1−/− and STAT1−/− treated with anti-IL-17. * and ** represent significant differences between WT and STAT1−/− mice (d) Primary breast tumor weights at the end of the experiment. (e) Number of metastatic tumor nodules in the lungs of the different groups. (f) Representative pictures of lungs from different experimental groups. (g) H&E staining of lung sections. Representative data from one of 2 independent experiments with an n-value of 5–7 mice per group. *p = 0.05, **p = 0.01, ***p = 0.0001. Our observed differences in IL-17 production in tumor bearing STAT1−/− and WT mice led us to investigate the impact of anti-IL-17 treatment on MDSC recruitment and subsequent tumor growth and metastasis. Female mice with similar primary tumor size were selected and randomized in groups that received anti-IL-17 treatment or IgG control antibody administration. Administration of IL-17 neutralizing antibody resulted in a slight reduction of primary tumor progression in STAT1−/− mice, however there were no statistical differences in any of the time points compared with STAT1−/− that received the isotype control antibody. Also, no differences were noted in the primary tumor sizes and weights between anti-IL-17 versus control antibody treated STAT1−/− mice at the termination of the experiment (Fig. 3c, d). Interestingly, STAT1−/− mice treated with anti-IL-17 contained significantly fewer metastatic nodules in their lungs compared with STAT1−/− mice that received isotype antibody (Fig. 3e, f and g).

Anti-IL-17 treatment affects immune populations in the spleen

We next evaluated the impact of IL-17 blockade on immune cell populations in the spleen. As expected we found decreased accumulation of Ly6G+CD11b+ cells in STAT1−/− mice treated with anti-IL-17 blocking antibody compared with STAT1−/− mice that received the isotype antibody (Fig. 4a, b). This is in line with previous reports showing reduced accumulation of neutrophils after IL-17 blockade in breast cancer bearing WT mice. Furthermore, we investigated whether blocking IL-17had an effect on T cell population and activation status by evaluating the expression of CD69. The proportion of both CD4+ and CD8+ cells was lower in tumor bearing STAT1−/− mice in comparison to tumor bearing WT mice (Fig. 4c and d). The frequency of CD4+ CD69+ activated T cells was similar between WT and STAT1−/− mice, but it was enhanced in STAT1−/− mice upon IL-17 neutralization (Fig. 4e). Unexpectedly, the frequency of CD8+CD69+ cells was higher in STAT1−/− mice and was unchanged upon anti-IL-17 treatment (Fig. 4f).
Figure 4.

Treatment with anti-IL-17 affects the expansion of Ly6G+ cells in the spleens of STAT1−/− mice. Mice were treated as in Fig. 3. At the end of the experiment different populations of splenocytes were analyzed by flow cytometry. (a) Representative dot plot of Ly6G+CD11b+ cells from STAT1−/− and anti-IL-17 treated STAT1−/− primary tumor bearing mice. (b) Frequencies of Ly6G+CD11b+ cells in the spleens of primary tumor bearing WT, STAT1−/− or STAT1−/− treated with anti-IL-17 neutralizing antibody. Frequencies of total (c) CD4+ cells and (d) CD8+ T cells. Frequencies of activated (e) CD4+ and (f) CD8+ in the different experimental conditions (WT CNTRL: naive WT mice, STAT1−/− CNTRL: naïve STAT1−/− mice, WT: tumor bearing WT mice, STAT1−/−: tumor bearing STAT1−/− mice, STAT1−/−αIL-17: tumor bearing STAT1−/− mice treated with anti-IL-17. *p = 0.05, **p = 0.01.

Treatment with anti-IL-17 affects the expansion of Ly6G+ cells in the spleens of STAT1−/− mice. Mice were treated as in Fig. 3. At the end of the experiment different populations of splenocytes were analyzed by flow cytometry. (a) Representative dot plot of Ly6G+CD11b+ cells from STAT1−/− and anti-IL-17 treated STAT1−/− primary tumor bearing mice. (b) Frequencies of Ly6G+CD11b+ cells in the spleens of primary tumor bearing WT, STAT1−/− or STAT1−/− treated with anti-IL-17 neutralizing antibody. Frequencies of total (c) CD4+ cells and (d) CD8+ T cells. Frequencies of activated (e) CD4+ and (f) CD8+ in the different experimental conditions (WT CNTRL: naive WT mice, STAT1−/− CNTRL: naïve STAT1−/− mice, WT: tumor bearing WT mice, STAT1−/−: tumor bearing STAT1−/− mice, STAT1−/−αIL-17: tumor bearing STAT1−/− mice treated with anti-IL-17. *p = 0.05, **p = 0.01.

STAT1−/− mice present enhanced tumor promoting microenvironment that is partially dependent on IL-17

STAT1 expression in tumor cells (fribrosarcoma) is necessary to regulate angiogenesis, tumorigenesis and metastasis, however the requirement of STAT1 expression in non-tumor cells for these functions has not been evaluated. Of note, vascularization of the tumor is largely associated with tumor metastasis. Because STAT1−/− mice showed enhanced lung metastasis, we investigated vascularization markers in the primary tumor. Interestingly, we found increased CD31 positive cells within the primary tumors of STAT1−/− mice compared with WT mice (Fig. 5a and b). STAT1−/− mice that received anti-IL-17 treatment showed slightly lower expression of CD31 compared with STAT1−/− mice that received the isotype antibody, however, the difference between the groups was statistically not significant (Fig. 5a and c). A recent report showed the critical importance of Ly6G+ cells for inducing tumor growth and metastasis during breast cancer. We therefore evaluated Ly6G+ cells in the primary tumor by immunofluorescence microscopy. As noted earlier, WT and STAT1−/− mice showed Ly6G+ cell accumulation in their primary tumors, which was higher in STAT1−/− mice compared with WT mice. Additionally, Ly6G+ cell accumulation was reduced by anti-IL-17 treatment, although this was not statistically significant (Fig. 5a, d and e). Since STAT1 is a critical mediator of IFN-γ signaling, we evaluated the expression of iNOS, a molecule known to be induced by IFN-γ/STAT1 pathway in macrophages and related to their tumoricidal functions. We found increased iNOS expression in primary breast tumors of WT and STAT1−/− compared with normal mammary gland tissue, however, the expression of iNOS in STAT1−/− mice was not affected upon IL-17 neutralization (Fig. 6a). We also evaluated the expression of genes related with tumor progression and metastasis. Interestingly, the expression of arginase, a molecule expressed in tumor promoting macrophages and canonical marker for M2 macrophages, was highly increased in STAT1−/− mice compared with WT mice. Of note, primary tumors from STAT1−/− mice treated with anti-IL-17 showed significantly reduced arginase expression (Fig. 6b). To further understand the mechanism of increased influx of Ly6G+ cells into primary tumors, we evaluated the expression of Cxcl1 a chemokine involved in neutrophil and granulocytic MDSC chemotaxis and poor breast cancer prognosis. In line with elevated numbers of granulocytic MDSCs in STAT1−/− mice, the expression of Cxcl1 was higher in STAT1−/− mice than WT mice; interestingly Cxcl1 expression in STAT1−/− mice was reduced upon treatment with anti-IL-17 (Fig. 6c). Since primary tumors of STAT1−/− showed enhanced vascularization compared with WT mice, we investigated the expression of genes related with angiogenesis such as Mmp9 and Vegf. Interestingly, Mmp9 expression was higher in STAT1−/− compared with WT tumor bearing mice (Fig. 6d). Anti-IL-17 treatment resulted in reduced the expression of Mmp9 in STAT1−/− tumor bearing mice. The expression of granzyme B and Vegf were similar in WT and STAT1−/− and unaffected upon anti IL17 treatment in primary tumors (Fig. 6e and f). It is also known that in addition to IL-17, G-CSF also plays an important role in neutrophil expansion. Hence, we measured the expression of G-CSF in the tumors. Our data indicates that tumors of STAT1−/− mice expressed significantly higher levels of G-CSF, compared with the tumors of WT mice (Fig. 6g). Anti-IL-17 treatment resulted in decreased expression of G-CSF in tumor bearing mice STAT1−/− mice. Next, we examined the effects of STAT1 deficiency on the expression of genes associated with granulocytic MDSC phenotype during breast cancer progression and metastasis. Gene expression of c-Kit, Nos2, Prok2, S100a8, S100a9, which have been shown to be associated with the granulocytic MDSC phenotype in some breast cancer models, were examined in sorted CD11b+ Ly6G+ cells in spleens of tumor bearing WT and STAT 1−/− mice. Interestingly, our data shows no significant differences in the expression of these genes between WT and STAT1−/− Ly6G+ CD11b+ cells (Supplementary Fig. 2c-g). It appears from our data that the role STAT1 plays in our breast cancer model is mostly associated with granulocytic MDSC accumulation, and does not affect the phenotype of these cells.
Figure 5.

STAT1−/− mice present enhanced vasculature in the primary tumor. Cryosections of primary tumors were processed for immunostaining of CD31 and Ly6G. (a) Representative Z stack projections of sections immunostained with anti Ly6G (red) anti CD31 (Green) and stained with DAPI (Blue). (b) Number of CD31+ cells per section. (c) Percentage of area stained by anti CD31. (d) Number of Ly6G+ cells per section. (e) Percentage of positive area stained with anti Ly6G. Data represent analysis of sections from 3 mice per group. *p = 0.05.

Figure 6.

Dysregulated gene expression in the primary tumor of STAT1−/− is partially dependent on IL-17. Gene expression analysis of primary tumors from WT, STAT1−/− or STAT1−/− treated with anti-IL-17 neutralizing antibody by real time PCR. Gene expressions of (a) Nos2, (b) Arg1, (c) Cxcl1, (d) Mmp9, (e) Granzyme B, (f) Vegf (g) G-CSF respectively. Data obtained from 3 different mice per group. *p = 0.05, **p = 0.01, ***p = 0.0001.

STAT1−/− mice present enhanced vasculature in the primary tumor. Cryosections of primary tumors were processed for immunostaining of CD31 and Ly6G. (a) Representative Z stack projections of sections immunostained with anti Ly6G (red) anti CD31 (Green) and stained with DAPI (Blue). (b) Number of CD31+ cells per section. (c) Percentage of area stained by anti CD31. (d) Number of Ly6G+ cells per section. (e) Percentage of positive area stained with anti Ly6G. Data represent analysis of sections from 3 mice per group. *p = 0.05. Dysregulated gene expression in the primary tumor of STAT1−/− is partially dependent on IL-17. Gene expression analysis of primary tumors from WT, STAT1−/− or STAT1−/− treated with anti-IL-17 neutralizing antibody by real time PCR. Gene expressions of (a) Nos2, (b) Arg1, (c) Cxcl1, (d) Mmp9, (e) Granzyme B, (f) Vegf (g) G-CSF respectively. Data obtained from 3 different mice per group. *p = 0.05, **p = 0.01, ***p = 0.0001.

Deficiency in CD8+ cell infiltration in primary tumor of STAT1−/− mice is restored upon anti-IL-17 treatment

A well-known mechanism for controlling tumor progression is the infiltration of T cells. We explored the infiltration of T cells into the tumors of WT and STAT1−/− mice by flow cytometry. STAT1−/− mice presented decreased accumulation of CD4+ T cells compared with WT mice (Fig. 7a). Additionally, CD4+CD69+ activated T cells were also less frequent in tumors from STAT1−/− mice and were unaffected by the treatment with anti-IL-17 (Fig. 7c). Similar to CD4+ T cells, the frequency of CD8+ T cells was lower in tumors of STAT1−/− mice compared with WT mice (Fig. 7b). The infiltrating CD8+ cells were found to be activated in WT and STAT1−/− mice as assessed by CD69 expression. Interestingly, the treatment with anti-IL-17 in STAT1−/− mice restored the frequencies of CD8+ cells similar to those found in tumors from WT mice (Fig. 7d). However, the expression of granzyme b, an important molecule that mediates tumor toxicity by CD8 and NK cells was similar between WT and STAT1−/− mice and unaffected upon anti-IL-17 treatment (Fig. 6e).
Figure 7.

Impaired infiltration of CD4+ and CD8+ cells in the primary tumor of STAT1−/− mice. Primary tumors from WT, STAT1−/− or STAT1−/− anti IL-17 treated mice were processed to obtain single cell suspension and analyzed by flow cytometry. Frequency of total (a) CD4+ and (b) CD8+ T cells within the primary tumor. Percentage of (c) CD4+ or (d) CD8+ activated T cells. Data from 4 mice per group. *p = 0.05, **p = 0.01.

Impaired infiltration of CD4+ and CD8+ cells in the primary tumor of STAT1−/− mice. Primary tumors from WT, STAT1−/− or STAT1−/− anti IL-17 treated mice were processed to obtain single cell suspension and analyzed by flow cytometry. Frequency of total (a) CD4+ and (b) CD8+ T cells within the primary tumor. Percentage of (c) CD4+ or (d) CD8+ activated T cells. Data from 4 mice per group. *p = 0.05, **p = 0.01.

Discussion

STAT1 has important roles in the function and activation of immune cells. Here, we show that STAT1 is necessary to control tumor growth and metastasis in an orthotopic model of breast cancer. Mechanistically, host STAT1 deficiency resulted in increased accumulation of Ly6G+Ly6CCD11b+ cells resembling granulocytic MDSCs and overexpression of molecules associated with tumor progression and metastasis such as Mmp9, Cxcl1, arginase and CD31. The accumulation of Ly6G+Ly6CCD11b+ cells in the spleen of STAT1−/− mice was reduced upon IL-17 blockade along with Mmp9, Cxcl1 and arginase expression. These data identify that STAT1 deficiency in non-tumor cells generates a favorable microenvironment for tumor metastasis. Also of note, the uncontrolled metastasis found in STAT1−/− mice can be restored to basal levels upon IL-17 blockade, suggesting a negative regulation of the STAT1 on IL-17 induced genes involved in tumor metastasis. Type I and II interferons are canonical inducers of STAT1 activation. The gene expression activated by IFNs in tumor cells has been shown to be mainly related to tumor apoptosis. Earlier studies have found beneficial effects of IFN therapy against tumors, however with high toxicity effects. Similarly, different studies using different models of breast cancer in mice have found an important role of STAT1 in suppressing tumor growth. Interestingly, STAT1−/− mice were reported to have an increased susceptibility of developing spontaneous breast tumors compared with WT mice. Also, targeted deletion of STAT1 in mammary epithelia cells resulted in accelerated tumor growth. Although these findings suggest a tumor intrinsic importance for STAT1 regulating tumor growth, our studies have found an important role for STAT1 in non-tumor cells for controlling tumor growth. Orthotopic injection of 4T1.2 cells resulted in increased tumor growth and metastasis in STAT1 deficient mice. In addition, we found a remarkable increase in lung metastasis in STAT1−/− mice. The relationship between STAT1 and metastasis suppression was shown to be important in a fibrosarcoma model, where loss of STAT1 in fibrosarcoma cells resulted in enhanced metastasis and angiogenesis. In fact, we found enhanced expression of angiogenic factors (MMP9) and enhanced micro vessel density in tumors of STAT1−/− mice which may favor the enhanced metastasis to the lungs. Increased Mmp9 expression in tumors of STAT1−/− mice is perhaps not surprising as IFN-β, a STAT1 activator, has been shown to downregulate MMP9 expression in neutrophils from tumor bearing mice to basal levels. An important cytokine for breast cancer growth and metastasis is IL-17 which drives the expansion of Ly6G+Ly6CCD11b+ cells. In a model of melanoma, IFN-γ deficient mice showed enhanced tumor growth whereas IFN-γ−/−/IL-17−/− double KO mice showed smaller tumors indicating that IFN-γ could be inhibiting tumor growth by suppressing IL-17 production. In the present study, we found that IL-17 was increased in tumor bearing STAT1−/− mice, and host STAT1 deficiency also led to CD11b+Ly6G+ accumulation in the spleens which was at least in part mediated by IL-17, as STAT1−/− mice treated with anti-IL-17 antibody showed a marked reduction in CD11b+Ly6G+ accumulation. In addition, we found that neutralization of IL-17 diminished expression of proangiogenic factors such as Mmp9, which is overexpressed in tumor infiltrating CD11b+Ly6G+ cells in breast cancer. Importantly, it has been shown that MMP9+ CD11b+Ly6G+ cells induce tumor angiogenesis and vascularization. It is possible that the reduced expression of Mmp9 in tumors from anti-IL-17-treated STAT1−/− mice in the present study could be due to the IL-17 induced Mmp9 expression in CD11b+Ly6G+ cells. Previous studies have shown that CXCL1 enhances accumulation of CD11b+Ly6G+ cells into the tumors which promotes tumor cell survival, and that IL-17 is an important cytokine for CXCL1 induction. In the present study, we found that Cxcl1 expression was increased in the tumors of STAT1−/− mice in an IL-17 dependent manner. Collectively, these data suggest a pathologic loop of mediators of tumor growth and metastasis that are negatively regulated by STAT1. The expansion of Ly6G+CD11b+ cells into the tumor and spleens has been reported to be IL-17 dependent by increasing the production of granulocytic colony stimulating-factor (G-CSF) that in turn expands and maintains an elevated number of Ly6G+CD11b+ cells. In line with these studies, we observed increased G-CSF expression in tumors of STAT1−/− mice compared with WT mice, which was decreased by anti-IL-17 treatment. Interestingly, the suppressor activity of monocytic MDSCs and TAMs has been reported to be dependent on STAT1 expression. In the present study, despite the higher accumulation of MDSCs in the spleen of tumor-bearing STAT1−/− mice, the proliferative capacity of splenic T cell in response to anti-CD3 stimulation was higher than that of WT mice, suggesting an intact capacity of T cells to proliferate even in presence of high numbers of MDSCs in STAT1−/− mice. This phenomenon could be explained by lower levels of PDL-1 expressing MDSCs observed in STAT1−/− tumor bearing mice compared with the WT tumor bearing mice (Supplementary Fig. 1c, d). These results are quite different from the spleen where similar levels of PDL1 are observed among monocytic and granulocytic MDSCs between tumor bearing WT and STAT−/− mice (Supplementary Fig. 1c, d). Indeed, granulocytic MDSCs from the spleen are able to suppress T cell proliferative responses equally between tumor bearing WT and STAT1−/− mice (Fig. 2f). Additionally, neutralization of IL-17 in STAT1−/− mice resulted in increased frequencies of activated CD4+ T cells in the spleens along with increased infiltration of activated CD8+ T cells into the tumors. The increased infiltration of CD8+ T cells is in line with a recent report showing induction of CD8+ T cell infiltration after IL-17 neutralization in WT mice bearing breast cancer and correlates with the decreased accumulation of Ly6G+CD11b+ cells after IL-17 neutralization. Interestingly, despite increased infiltration of CD8+ T cells, granzyme b expression was not restored by anti-IL-17 treatment of STAT1−/− mice. Of course, granzyme b is also expressed by NK cells as well as other non-cytotoxic cells such as basophils and mast cells. Moreover, as a protease which is present in granules of activated CD8+ T cells, increased de novo gene expression of Gzmb may not be reflective of an increase in the number of activated CD8+ T cells. We observed interesting differences in activation marker expression between CD4+ and CD8+ T cells in the tumors and spleens of tumor bearing STAT1−/− mice treated with IL-17 blocking antibody. We previously showed a differential requirement for STAT1 in CD4+ vs. CD8+ T cell activation. In that study, STAT1−/− CD4+ and CD8+ T cells were activated via separate signaling pathways. It is possible that similar effects of IL-17 blockade in STAT1−/− CD4+ and CD8+ T cells are also occurring in our breast cancer model. It has been largely thought that the control of tumor growth by STAT1 is mainly related to the ability to induce apoptosis and provide immunosurveillance of tumors. Our data provide evidence of a new mechanism of how STAT1 can repress tumor growth by limiting the accumulation of MDSCs and the expression of pro-angiogenic genes such as arginase1, Cxcl1 and Mmp9, which was reversed upon IL-17 blockade. This implies that anti-IL17 treatment may decrease the expression of pro-metastatic factors even under conditions where the host is immunologically compromised and is unable to control breast cancer metastasis. Finally, our study suggests investigating the role of inhibition of this inflammatory loop by using STAT1 agonists to control breast cancer growth and metastasis.

Methods

Mice and tumor injections

WT Balb/c mice were purchased from Envigo (previously Harlan laboratories, Indianapolis, IN). STAT1−/− mice were generated and maintained as described previously. All female mice were age matched and maintained in compliance with the guidelines of OSU ULAR. 4T1.2 tumor cells were obtained from 4T1.2, a subclone of murine cell line 4T1 cells and were kindly provided by Dr. Robin Anderson. The cells were cultured in RPMI1640 medium supplemented with 10% FBS and 1% penicillin and streptomycin (Life technologies, Carlsbad, CA) at 370C in 5%CO2. All experimental mice were injected with 105 4T1.2 cells in the mammary fat pad by following the OSU IACUC guidelines.

Tumor measurements and lung metastasis counts

Breast tumors were measured every 3 d by using electronic calipers and volume was calculated by using the following formula (volume = 0.52 X a2 X b), where “a” is the smaller superficial diameter and “b” is the larger superficial diameter (https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3271140/). Experimental mice were killed at specified time intervals, and tumor weights and lung metastasis were calculated. Lungs were inflated and stored in Bouin's solution (Sigma Aldrich, St. Louis, MO). Lung metastasis was calculated by counting the total number of visible nodules present on the lungs.

Immunohistochemistry

Tumors from WT and the various treatment groups of STAT1 mice were cryopreserved in OCT Compound (Fisher Healthcare 23–730–571). Sections 7 micron thick were obtained (Leica Microsystems CM1850 Cryostat) and stained with anti CD31 (Santa Cruz biotechnologies) and anti Ly6G-PE conjugated antibody (Bio Legend) overnight at 4°C. Anti- goat secondary antibody was added during 2.5 hours to detect anti CD31. The sections were stained with DAPI (Bio Legend). Images were captured by using a confocal microscope Zeiss LSM700 (Zeiss, Dublin CA, USA). Z stack projections were obtained with the Zeiss black software and analysis of number of positive cells and percentage of positive area were analyzed by Image J software.

Flow cytometry analysis

Mice were harvested at respective time points, single cell suspensions were prepared from spleens and tumors. Briefly, spleens were mashed and red blood cells were lysed by ACK lysis buffer. Tumors were pulverized and incubated in buffer containing 1mg/ml DNase and 1mg/ml collagenase at 450C for 15 minutes. Spleens and tumor suspensions were passed through 60micron filters and cells were stained with the respective stain mixes. CD3(clone:145–2C11), CD4(clone: GK1.5), CD8(clone:53–6.7), CD69(clone:H1.2F3), CD11b(clone:M1/70), Ly6g(clone:IA8), Ly6C(clone:HK1.4) and PDL1(clone:10F.9G2) antibodies were purchased from Bio Legend (San Diego, CA). Cells were acquired through Fluorescence Activated Cell Sorter (FACS, BD biosciences, San Jose, CA, USA). Analysis was performed with Flow Jo software (Tree Star Inc., Ashland, OR, USA).

ELISA and T cell proliferation assays

Single cell suspensions were prepared from the spleens at the concentration of 5 × 106/ml. Cells were incubated with 2µg/ml of LEAF purified anti-mouse CD3e (Purchased from Bio Legend, San Diego, CA) in RPMI 1640 supplemented with 10% FBS and 1% penicillin and streptomycin at 370C in 5%CO2 for 72hrs. Cell culture supernatants were collected, production of IL-17 cytokine was analyzed by ELISA. All the capture, detection antibodies and standards of IL-17 were purchased from Bio Legend (San Diego, CA, USA). Cell proliferation was measured by Alamar blue reduction method (Bio-Rad AbD Serotec Inc., Raleigh, NC) Briefly, 10% Alamar blue (Life technologies, Carlsbad, CA) was added at 60 hrs of incubation and reduction of Alamar blue was measured at 72 hrs by measuring the absorbance at 570nm and 670 nm by using spectramax microplate reader and Softmax pro software (Molecular Devices LLC, Sunnyvale, CA, USA).

RT-PCR analysis

Total RNA was extracted from both the spleens and tumors of the experimental mice by TRIzol extraction method (purchased from Life technologies (Carlsbad, CA)). cDNAs were prepared by iScript reverse transcriptase and RT PCR reactions were done by using I Q SYBR green super mix and CFX 96 RT-PCR cycler (purchased from Bio-Rad, Hercules, CA, USA). Primers were selected from the Primer bank website (http://pga.mgh.harvard.edu/primerbank). Data obtained was normalized by using housekeeping gene β-actin and presented as the fold induction over WT mice.

In vivo administrations of neutralizing antibodies

For the in vivo antibody administration, STAT1−/− mice with similar sized tumors were selected, divided into 2 groups. One group was treated with 200µg of InVivo MAb anti-mouse IL-17A (clone: 17F3), Bio X Cell, (Lebanon, NH, USA) and the second group was treated with a similar concentration of isotype control every alternate day for 2 weeks via intraperitoneal route.

Statistics

Statistical differences between groups were analyzed in Graph Pad Prism by using Student T test. Differences were considered statistically significant when the p value was below to 0.05.
  48 in total

1.  Endogenous IL-17 as a mediator of neutrophil recruitment caused by endotoxin exposure in mouse airways.

Authors:  Masahide Miyamoto; Olof Prause; Margareta Sjöstrand; Martti Laan; Jan Lötvall; Anders Lindén
Journal:  J Immunol       Date:  2003-05-01       Impact factor: 5.422

Review 2.  Cytotoxic and non-cytotoxic roles of the CTL/NK protease granzyme B.

Authors:  Inna S Afonina; Sean P Cullen; Seamus J Martin
Journal:  Immunol Rev       Date:  2010-05       Impact factor: 12.988

Review 3.  Antitumour actions of interferons: implications for cancer therapy.

Authors:  Belinda S Parker; Jai Rautela; Paul J Hertzog
Journal:  Nat Rev Cancer       Date:  2016-03       Impact factor: 60.716

4.  STAT1 drives tumor progression in serous papillary endometrial cancer.

Authors:  Budiman Kharma; Tsukasa Baba; Noriomi Matsumura; Hyun Sook Kang; Junzo Hamanishi; Ryusuke Murakami; Melissa M McConechy; Samuel Leung; Ken Yamaguchi; Yuko Hosoe; Yumiko Yoshioka; Susan K Murphy; Masaki Mandai; David G Hunstman; Ikuo Konishi
Journal:  Cancer Res       Date:  2014-09-29       Impact factor: 12.701

Review 5.  Interferons, immunity and cancer immunoediting.

Authors:  Gavin P Dunn; Catherine M Koebel; Robert D Schreiber
Journal:  Nat Rev Immunol       Date:  2006-11       Impact factor: 53.106

6.  Stat1 is a suppressor of ErbB2/Neu-mediated cellular transformation and mouse mammary gland tumor formation.

Authors:  Jennifer F Raven; Virginie Williams; Shuo Wang; Michel L Tremblay; William J Muller; Joan E Durbin; Antonis E Koromilas
Journal:  Cell Cycle       Date:  2011-03-01       Impact factor: 4.534

7.  Genomic analysis of a spontaneous model of breast cancer metastasis to bone reveals a role for the extracellular matrix.

Authors:  Bedrich L Eckhardt; Belinda S Parker; Ryan K van Laar; Christina M Restall; Anthony L Natoli; Michael D Tavaria; Kym L Stanley; Erica K Sloan; Jane M Moseley; Robin L Anderson
Journal:  Mol Cancer Res       Date:  2005-01       Impact factor: 5.852

8.  IL17 Promotes Mammary Tumor Progression by Changing the Behavior of Tumor Cells and Eliciting Tumorigenic Neutrophils Recruitment.

Authors:  Luciana Benevides; Denise Morais da Fonseca; Paula Barbim Donate; Daniel Guimarães Tiezzi; Daniel D De Carvalho; Jurandyr M de Andrade; Gislaine A Martins; João S Silva
Journal:  Cancer Res       Date:  2015-07-24       Impact factor: 12.701

9.  Requirement of interleukin 17 receptor signaling for lung CXC chemokine and granulocyte colony-stimulating factor expression, neutrophil recruitment, and host defense.

Authors:  P Ye; F H Rodriguez; S Kanaly; K L Stocking; J Schurr; P Schwarzenberger; P Oliver; W Huang; P Zhang; J Zhang; J E Shellito; G J Bagby; S Nelson; K Charrier; J J Peschon; J K Kolls
Journal:  J Exp Med       Date:  2001-08-20       Impact factor: 14.307

10.  Elevated CXCL1 expression in breast cancer stroma predicts poor prognosis and is inversely associated with expression of TGF-β signaling proteins.

Authors:  An Zou; Diana Lambert; Henry Yeh; Ken Yasukawa; Fariba Behbod; Fang Fan; Nikki Cheng
Journal:  BMC Cancer       Date:  2014-10-24       Impact factor: 4.430

View more
  14 in total

1.  STAT1 inhibits T-cell exhaustion and myeloid derived suppressor cell accumulation to promote antitumor immune responses in head and neck squamous cell carcinoma.

Authors:  Nathan Ryan; Kelvin Anderson; Greta Volpedo; Omar Hamza; Sanjay Varikuti; Abhay R Satoskar; Steve Oghumu
Journal:  Int J Cancer       Date:  2019-11-29       Impact factor: 7.396

2.  Long non-coding RNA AFAP1-AS1 plays an oncogenic role in promoting cell migration in non-small cell lung cancer.

Authors:  Juan He; Ke Wu; Chenglin Guo; Jian-Kang Zhou; Wenchen Pu; Yulan Deng; Yuanli Zuo; Yun Zhao; Lunxu Liu; Yu-Quan Wei; Yong Peng
Journal:  Cell Mol Life Sci       Date:  2018-10-06       Impact factor: 9.261

3.  Predicting STAT1 as a prognostic marker in patients with solid cancer.

Authors:  Jinguo Zhang; Fanchen Wang; Fangran Liu; Guoxiong Xu
Journal:  Ther Adv Med Oncol       Date:  2020-05-08       Impact factor: 8.168

4.  Deficiency in STAT1 Signaling Predisposes Gut Inflammation and Prompts Colorectal Cancer Development.

Authors:  Sonia Leon-Cabrera; Armando Vázquez-Sandoval; Emmanuel Molina-Guzman; Yael Delgado-Ramirez; Norma L Delgado-Buenrostro; Blanca E Callejas; Yolanda I Chirino; Carlos Pérez-Plasencia; Miriam Rodríguez-Sosa; Jonadab E Olguín; Citlaltepetl Salinas; Abhay R Satoskar; Luis I Terrazas
Journal:  Cancers (Basel)       Date:  2018-09-19       Impact factor: 6.639

Review 5.  Metabolic Regulation of Myeloid-Derived Suppressor Cell Function in Cancer.

Authors:  Yufei Wang; Anna Jia; Yujing Bi; Yuexin Wang; Guangwei Liu
Journal:  Cells       Date:  2020-04-18       Impact factor: 6.600

Review 6.  Targeting Myeloid-Derived Suppressor Cell, a Promising Strategy to Overcome Resistance to Immune Checkpoint Inhibitors.

Authors:  Aohan Hou; Kaiyu Hou; Qiubo Huang; Yujie Lei; Wanling Chen
Journal:  Front Immunol       Date:  2020-05-15       Impact factor: 7.561

7.  Immune Suppression Mediated by STAT4 Deficiency Promotes Lymphatic Metastasis in HNSCC.

Authors:  Kelvin Anderson; Nathan Ryan; Greta Volpedo; Sanjay Varikuti; Abhay R Satoskar; Steve Oghumu
Journal:  Front Immunol       Date:  2020-01-15       Impact factor: 7.561

8.  IL-17 Affects the Progression, Metastasis, and Recurrence of Laryngeal Cancer via the Inhibition of Apoptosis through Activation of the PI3K/AKT/FAS/FASL Pathways.

Authors:  Yang Song; Ming Yang; Hongjian Zhang; Yan Sun; Ye Tao; Huihui Li; Jing Zhang; Yuncheng Li; Jianming Yang
Journal:  J Immunol Res       Date:  2020-12-18       Impact factor: 4.818

9.  Proteomic Characterization of Cytoplasmic Lipid Droplets in Human Metastatic Breast Cancer Cells.

Authors:  Alyssa S Zembroski; Chaylen Andolino; Kimberly K Buhman; Dorothy Teegarden
Journal:  Front Oncol       Date:  2021-06-01       Impact factor: 6.244

Review 10.  CXCL1: Gene, Promoter, Regulation of Expression, mRNA Stability, Regulation of Activity in the Intercellular Space.

Authors:  Jan Korbecki; Katarzyna Barczak; Izabela Gutowska; Dariusz Chlubek; Irena Baranowska-Bosiacka
Journal:  Int J Mol Sci       Date:  2022-01-12       Impact factor: 5.923

View more

北京卡尤迪生物科技股份有限公司 © 2022-2023.