Cd Library Dc101 Software Engineer
It selects the indicated CD or CDs. It's implemented as a small console application. When you start it the first time it will ask you to enter some configuration information.
Cd Library Ii Software
Abstract
Glioblastoma (GBM) is a hypervascular and malignant form of brain tumors. Anti-angiogenic therapies (AAT) were used as an adjuvant against VEGF-VEGFR pathway to normalize blood vessels in clinical and preclinical studies, which resulted into marked hypoxia and recruited bone marrow derived cells (BMDCs) to the tumor microenvironment (TME). In vivo animal models to track BMDCs and investigate molecular mechanisms in AAT resistance are rare. We exploited recently established chimeric mouse to develop orthotopic U251 tumor, which uses as low as 5×106 GFP+ BM cells in athymic nude mice and engrafted >70% GFP+ cells within 14 days. Our unpublished data and published studies have indicated the involvement of immunosuppressive myeloid cells in therapeutic resistance in glioma. Similarly, in the present study, vatalanib significantly increased CD68+ myeloid cells, and CD133+, CD34+ and Tie2+ endothelial cell signatures. Therefore, we tested inhibition of CSF1R+ myeloid cells using GW2580 that reduced tumor growth by decreasing myeloid (Gr1+ CD11b+ and F4/80+) and angiogenic (CD202b+ and VEGFR2+) cell signatures in TME. CSF1R blockade significantly decreased inflammatory, proangiogenic and immunosuppressive molecular signatures compared to vehicle, vatalanib or combination. TCK1 or CXCL7, a potent chemoattractant and activator of neutrophils, was observed as most significantly decreased cytokine in CSF1R blockade. ERK MAPK pathway was involved in cytokine network regulation. In conclusion, present study confirmed the contribution of myeloid cells in GBM development and therapeutic resistance using chimeric mouse model. We identified novel molecular networks including CXCL7 chemokine as a promising target for future studies. Nonetheless, survival studies are required to assess the beneficial effect of CSF1R blockade.
1. INTRODUCTION
Glioblastoma (GBM), a grade IV glioma classified by World Health Organization (WHO), is considered highly malignant, vascular and invasive subtype[55]. GBM is most lethal during first year after initial diagnosis despite surgical resection, radiotherapy and/or chemotherapy [55, 74]. Median survival of patients diagnosed with GBM is only 12 to 15 months [55, 74]. Hypoxia and neovascularization are histopathologic features of GBM [12]. Anti-angiogenic therapies (AAT) were used as an adjuvant mainly against VEGF-VEGFR pathway to normalize tumor vasculatures. However, all of them provided minimal to none effect with no change in overall survival [38, 82]. Therefore, current challenge is to investigate mechanisms of undesirable outcomes in GBM clinical trials. Preclinical studies involving AATs have shown marked hypoxia, increased homing of bone marrow derived cells (BMDCs) to the tumor and activation of alternative pathways of neovascularization [9, 41].
BMDCs play a pivotal role in tumor development[66] and endothelial progenitor cells (EPCs) from BM pool is recruited to tumor microenvironment (TME) [3, 15, 19, 24]. BMDCs associated therapeutic resistance falls under evasive or adaptive resistance, where tumor itself after an initial response phase acquires evasive properties against therapeutic blockade by inducing alternate mechanisms that enable neovascularization, leading to renewed tumor growth and progression[10]. Distinct potential mechanism of resistance might be at cellular level [33] mediated through up-regulation of HIF1-α followed by induction of SDF1α, secretion of proangiogenic factors and recruitment of CXCR4+ BMDCs to the tumor [3, 4, 15, 19]. These recruited cells were characterized as pro-angiogenic CD45+VEGFR2+ EPCs, or CD45+Tie2+ monocytes [5, 20]. Interestingly, lin-ckit+Sca-1+ and their derived cells demonstrated recruitment to tumor but do not functionally contribute to tumor neovascularization [69]. BMDCs derived MMP9 modulates neovessel remodeling, thereby playing role in tumor growth [4, 67].
Studies have indicated that resistance to AAT has profound involvement of immune system [1, 16, 23, 45, 56, 68, 70]. Role of myeloid cells in tumor angiogenesis is an established phenomenon as shown by other authors [2, 46, 71, 72, 78] and supported by our previous study (unpublished), which showed that majority of bone marrow derived GFP+ cells acquire both myeloid signature (CD68) and endothelial signatures (CD202b and CD34) in the TME under vatalanib (VEGFR tyrosine kinase inhibitor) treatment. Tissue-resident macrophages originate from yolk-sac-derived erythro-myeloid progenitors [28], whereas, SDF-1 played an important role in the invasiveness of brain tumor and infiltration of macrophages from bone marrow [75]. AATs in glioma were associated with increased myeloid cell infiltration and stem cell accumulation [57]. However, investigations whether those phenotypes have bone marrow component, were lacking. In other study, authors noticed infiltration in myeloid populations in the tumor bulk and in the infiltrative regions after AAT [45]. Together, studies suggest that immune suppressive myeloid cells, especially myeloid derived suppressor cells (MDSCs) and tumor associated macrophages (TAM) [27, 35, 73], may participate in escape from AATs and represent a potential biomarker of resistance with potential therapeutic target in GBM [45, 76].
Several chemokines such as macrophage colony-stimulating factor-1 (MCSF/CSF1) and monocyte chemotactic protein-1 (MCP1/CCL2) are known to contribute in the recruitment of TAMs to the tumor [43, 65]. CSF1R expression has been reported on immunosuppressive myeloid cells and dendritic cells [47, 51, 77]. CSF1-CSF1R signaling regulates survival, differentiation, and proliferation of monocytes and macrophages [31, 34], and have critical role in angiogenesis and tumor progression [58, 83]. Therefore, goals of present study are to (1) investigate the effect of CSF1R blockade on orthotopic glioma development in a recently established preclinical chimeric mouse model, (2) to evaluate whether CSF1R blockade alone or in combination with VEGFR2 blockade could inhibit the homing of myeloid BMDCs to the glioma, (3) to identify signature immune cell populations following CSF1R inhibition that could have profound role in glioma growth and (4) to investigate key secreted molecular signatures in GBM TME following CSF1R inhibition.
2. MATERIALS AND METHODS
All Animal related experimental procedures were approved by the Institutional Animal Care and Use Committee and Institutional Review Board of Georgia Regents University (animal protocol #2014-0625). All efforts were made to ameliorate suffering of animals. CO2 with secondary method was used to euthanize animals for tissue collection at the end of the study.
2.1 Establishing chimeric mouse
Chimeric mouse for orthotopic U251 glioma was established with IACUC approved protocol and published method [6]. Transgenic mice with universally expressing GFP under the human ubiquitin C promoter (C57BL/6-tg(UBC-GFP)30Scha) were used as donors (Jackson Laboratory, Main, USA). Athymic nude mice (NCr-nu/nu) were used as recipients (Charles River, Frederick, MD, USA), and were whole body irradiated with sub-lethal dose of 6Gy (using Cs-137 source). After 24 hours, recipient mice were injected intravenous (n=3) and intraperitoneal (n=3) routes with BM cells (5×106 cells) collected from donor transgenic mice. All mononuclear cells were separated from red blood cells using lymphocyte cell separation media (Corning, Cellgro, USA), counted and 5×106 cells/100μl were injected into each mouse. Ten microliter of blood (from orbital sinus) were collected from each mouse on days 7 to 56 (n=3 each time point) following transplantation of BM to determine BM engraftment efficiency (GFP positivity) in peripheral blood using flowcytometer. Cells from athymic mice without irradiation and GFP+ cell transplantation were used as control for flow cytometry.
2.2 Animal model of human glioma
Precisely, animals were anesthetized with 100 mg/kg ketamine and 15 mg/kg xylazine i.p. The surgical zone was swabbed with betadine solution, the eyes coated with Lacri-lube and the animals were immobilized in a small animal stereotactic device (Kopf, Cayunga, CA). After draping, a 1-cm incision was made 2 mm to the right of the midline 1 mm retro-orbitally; the skull exposed with cotton-tip applicators and a 23G needle tip was used to drill a hole 2 mm to the right of the bregma, taking care not to penetrate the dura. A 10 μL Hamilton syringe with a 26G-needle containing tumor cells (2.4×105) in 3 μl was lowered to a depth of 2.5 mm, and then raised to a depth of 2 mm. During and after the injection, careful note was made of any reflux from the injection site. After completing the injection, we waited 2-3 minutes before withdrawing in a stepwise manner. The surgical hole was sealed with bone wax. Finally, the skull was swabbed with betadine before suturing the skin over the injection site.
2.3 In vivo multispectral optical imaging
Multispectral optical images were acquired using excitation profiles of 460-480 nm range and emission of 535 nm to monitor the GFP positivity on days 7, 14 and 21 after tumor cell implantation. RFP positivity was determined at excitation 587 nm and emission 610 nm. All optical imaging data was acquired by Spectral AMI (Spectral Instruments Imaging, LLC) machines and analyzed by AMI view software.
2.4 Drug treatments
Orthotopically implanted chimeric mice with U251 tumor cells were allowed to grow for 7 days and then started oral treatments of either vehicle or receptor tyrosine kinase inhibitors (vatalanib (50mg/kg/day) [81], GW2580 (160 mg/kg/day) [77] and combination of both GW2580+ vatalanib, daily for two weeks. Seven days waiting period was followed after tumor implantation to mimic clinical scenario, where treatment is being done following detection of tumor.
2.5 In vivo magnetic resonance imaging (MRI)
All MRI experiments were conducted using a 7 Tesla 12 cm (clear bore) magnet interfaced to a varian console with actively shielded gradients of 49 gauss/cm and 100μs rise times or a horizontal 7 Tesla BioSpec MRI spectrometer (Brucker Instruments, Bellerica, MA) equipped with a 12 cm self-shielded gradient set (45 gauss/cm max). Detailed MRI procedure was adopted from our several previous publications [18, 22, 25, 53, 54]. An appropriate state of anesthesia was obtained with isoflurane (2.5% for induction, 0.7% to 1.5% for maintenance in a 2:1 mixture of N2:O2). After positioning using a triplanar FLASH sequence, MR studies were performed using pre-contrast T1, T2-weighted and post contrast T1-weighted MRI scans with following parameters (1) Standard T1-weighted multislice sequence (TR/TE=500/10 ms, 256×256 matrix, 13-15 slices, 1 mm thick slice, 32 mm field of view (FOV), # of averages=4). (2) T2-mapping sequence (2D multi-slice, multi-echo (MSME) sequence, TE=10, 20, 30, 40, 50, 60 msec, TR=3000 msec, 256×256 matrix, 13-15 slices, 1 mm thick slice, 32 mm field of view (FOV), # of averages=2). Post contrast T1WI was used to determine volume of tumors in vehicle and drug treated mice by drawing irregular ROI to encircle whole tumor in each image section containing tumor using ImageJ software, and area was then multiplied by thickness of image slice to determine volume (cm3). Two investigators blinded to the animal groups determined tumor volume.
2.6 Collection of GFP+ cells and determination of immune cell populations
Freshly isolated brain samples were separated into left and right (tumor bearing) hemispheres from each group and were homogenized to pass through 40μ cell strainer to make single cell. Similarly, cells were collected from spleen and BM. Cells were labeled with antibodies (BioLegend) such as CD45, CD4, CD8, Gr1, CD11b, F4/80, CD68, CD133, CD31, CD34, CD202b (Tie2), and CD309 (VEGFR2) (other than FITC) to identify BM recruited cell types (GFP+) in the tumor as well as phenotypes of GFP+ cells in spleen and bone marrow. Flow cytometry data was acquired using Accuri C6 machine (BD Biosciences) and analyzed by BD Accuri C6 software.
2.7 Immunofluorescence study
After MRI at day 22, animals were euthanized and brains were collected for frozen tissue sections and later stained for immunofluorescence study to determine expression of myeloid cell signature markers CD11b (Abcam), F4/80 (Santa Cruz Biotechnology), p-ERK (Cell Signaling), Ki67 (DAKO), CXCL7 (Millipore) and CD31 (Abcam) at the site of tumor. Migration and incorporation pattern of GFP+ BMDCs was determined in different regions of the tumor.
2.8 Western Blot Analysis
Tumors were collected and processed for protein isolation using TPER, tissue protein extraction reagent for tissue and Pierce RIPA buffer for tumor cells (Thermo Scientific, USA). Protein concentrations were estimated with Pierce, BCA protein assay kit (Thermo Scientific, USA) and separated by standard Tris/Glysine/SDS gel electrophoresis. Membranes were incubated with primary antibodies against IDO (1:1000, Santa Cruz Biotechnology), p-ERK and ERK (1:1000, Cell Signaling), and α-Tubulin (1:5000, Abcam) followed by horseradish peroxidase-conjugated secondary antibody (1:5000, Biorad). The blots were developed using a Pierce SuperSignal West Pico Chemiluminescent substrate kit (Thermo Scientific, USA). Western blot images were acquired by Las-3000 imaging machine (Fuji Film, Japan).
2.9 Protein Array
Differentially treated tumors were collected and isolated proteins were processed for mouse cytokine array (AAM-CYT-1000-8, Ray Biotech). Membranes were imaged using Las-3000 imaging machine (Fuji Film, Japan). All signals (expression intensity) emitted from the membrane were normalized to the positive control spots of corresponding membrane using Image J software. Later, expression intensities were imported as an excel file in Partek Genomic Suit 6.6 for heat map and box-plot analysis, and web based Ingenuity Pathway Analysis (IPA) interface for network analysis.
2.10 Statistical analysis
Quantitative data was expressed as mean ± SD and analyzed through Student t-test using GraphPad Prism. Differences were considered statistically significant at p value <0.05.
3. RESULTS
3.1 Establishing chimeric mouse model
As reported before athymic nude mice (n=3, each group) were transplanted with 5×106 BM cells from GFP+ mouse after whole body irradiation of sub-lethal dose (6Gy). Engraftment efficiency (GFP+) was determined at each week (Days 7 to 56). Surprisingly, i.v. injection (n=3) of 5×106 mononuclear cells achieved >70% engraftment by day 14 and >80 % engraftment by day 56 (unpublished data). However, i.p. injection (n=3) could achieve highest engraftment efficiency of 22.3 % on day 56 (unpublished data). We used athymic nude mice in Balb/c background to establish chimera. Surprisingly, we did not observe any GVHD associated symptoms or mortality. Flow cytometry analysis using day 22 peripheral blood showed that 10-13% CD4 and 1-2% CD8 cells were transfused while establishing GFP+ chimeric mice (Supplementary fig 1). Further, we decided to use i.v. injection of 5×106 BM cells for 14 days to establish chimera for tumor studies.
3.2 Effect of GW2580 on tumor growth and GFP+ cell infiltration to the tumor
Chimeric animals were implanted with orthotopic glioma (U251) and treated with vehicle (n=17), vatalanib (n=10), GW2580 (n=5) and GW2580+ vatalanib (n=5) from day 8-21. All animals underwent MRI on day 22 following implantation of tumors. We observed significant reduction in tumor growth in GW2580 and GW2580+ vatalanib groups (Fig. 1A, upper left and right panels). In vivo optical images obtained on day 21 also showed increasing trend of cherry red+ tumor in vatalanib group and decreasing trend of cherry red+ tumor after treatments of GW2580 and GW2580+ vatalanib groups compared to vehicle (Fig. 1A, Lower panel). Vatalanib treatment resulted in increasing trend of GFP+ cells accumulation in tumor. Decreasing trend of accumulation of GFP+ cells were observed in GW2580 and GW2580+ vatalanib groups compared to vehicle (Fig. 1B). Optical imaging data was validated by flow cytometry analysis, which proved increased infiltration of GFP+ cells in vatalanib group and decreased infiltration of GFP+ cells in tumors that were treated with GW2580 and GW2580+ vatalanib (Fig. 1C, upper and lower panels).
Effect of vatalanib and CSF1R inhibitor on tumor growth and GFP+ cell infiltration to the tumor(A) MRI images showing significant reduction in tumor growth in GW2580 and GW2580+ vatalanib groups (upper left and right panels) compared to vehicle. Vatalanib did not decrease tumor growth. In vivo optical images showing increasing trend of cherry red+ tumor in vatalanib group and decreasing trend of cherry red+ tumor after treatments of GW2580 and GW2580+ vatalanib groups compared to vehicle (Lower panels). (B)In vivo optical images showing increasing trend of GFP+ cells accumulation in tumor following vatalanib treatment and decreasing trend of accumulation of GFP+ cells in GW2580 and GW2580+ vatalanib groups compared to vehicle (C) Flow cytometry analysis showing increased infiltration of GFP+ cells in vatalanib group and decreased infiltration of GFP+ cells in tumors that were treated with GW2580 and GW2580+ vatalanib (upper and lower panels). Shown is one of the two experiments performed. Quantitative data is expressed in mean ±SD. * P<0.05 and **P<0.01.
3.3 Effect of GW2580 on myeloid and endothelial cell signatures in tumor
Chimeric animals were implanted with orthotopic glioma (U251) and treated with vehicle, vatalanib, GW2580 and GW2580+vatalanib from day 8-21. Brain, spleen and BM (n=3 each) were collected after MRI for flow cytometry analysis (Fig. 2A, upper and lower panels, spleen and BM data not shown). Infiltration of GFP+ cells was mostly seen at the periphery of the vatalanib treated tumor (Supplementary fig. 2). Vatalanib treatment significantly increased CD68+ myeloid cells, and CD133+, CD34+ and CD202b+ or Tie2+ endothelial cell signatures. Significant decrease in myeloid (Gr1+ CD11b+ and F4/80+) and endothelial (CD202b+ and CD309+ or VEGFR2+) cell signatures were observed in GW2580 group compared to vehicle. However, CD133+ population was significantly increased in GW2580+ vatalanib group (upper and lower panels). Further, we observed that vatalanib induced endothelial marker expression was associated with the increased vasculature compared to vehicle group. This was characterized by increased CD31 expression on GFP+ cells. GW2580 treated tumors showed significantly decreased CD31 expression compared to vehicle and vatalanib treated groups (Supplementary fig. 3). In immunofluorescence study, vatalanib treatment showed increased myeloid cell signatures (CD11b and F4/80) compared to vehicle (Fig. 2B). Decreased expression of myeloid cell signatures were observed following GW2580 and GW2580+ vatalanib treatments, compared to vehicle and vatalanib groups (Fig. 2B). Increasing trend of expression of immunosuppressive molecule, IDO, was seen with vatalanib, which was significantly decreased in GW2580 and GW2580+ vatalanib treated tumors (Fig. 2C).
Effect of vatalanib and CSF1R inhibitor on myeloid and endothelial cell signatures in GBM tumor(A) Flowcytometry analysis showing significantly increased CD68+ myeloid cells, and CD133+, CD34+ and CD202b+ or Tie2+ endothelial cell signatures following vatalanib treatment. Significant decrease in myeloid (Gr1+ CD11b+ and F4/80+) and endothelial (CD202b+ or Tie2+ and CD309+ or VEGFR2+) cell signatures was seen following GW2580 treatment compared to vehicle. CD133+ population was significantly increased in GW2580+ vatalanib treated group (upper and lower panels). (B) Immunofluorescence study showing increased myeloid cell signatures (CD11b and F4/80) in vatalanib and decreased myeloid (CD11b and F4/80) markers following GW2580 and GW2580+ vatalanib treatments compared to vehicle (yellow arrows). (C) Western blot data showing increased IDO (immune suppressive) expression in vatalanib treated tumors, which was decreased in GW2580 and GW2580+ vatalanib groups (lower panels). Shown is one of the two experiments performed. Quantitative data is expressed in mean ±SD. * P<0.05, **P<0.01 and *** P<0.001.
3.4 Treatment induced molecular networks in GBM TME
U251 tumors were collected following vehicle, vatalanib, GW2580 and GW2580+ vatalanib treatment groups (n=3 each). Tumor isolated proteins were processed for 96-plex protein array (Ray Biotech) followed by analysis of blots for expression intensities. Expression intensities were used to generate heatmap, which showed elevated expression of cytokines in following order of treatments; vatalanib > GW2580+ vatalanib > GW2580. Cytokines were mostly downregulated in GW2580 group and upregulated in vatalanib groups. Surprisingly, combination of vatalanib with GW2580 increased cytokine expression, which were decreased in GW2580 group (Fig. 3A and supplementary fig. 4). As represented by box-plot data, overall expression intensities were changed among groups. GW2580 group showed significantly decreased expression intensities of cytokines (Fig. 3B). By grouping cytokines based on their function, we found that GW2580 significantly decreased cytokines involved in inflammation (TCA3 or CCL11, TNFRI, CXCL9, GITR IL15, IL17, IL17RB, lungkine or CXCL15, MDC or CCL22, OPG, Resistin, Shh-N, TCK1 or CXCL7, TROY and TSLP), angiogenesis (bFGF, ICAM1, ITAC or CXCL11, VEGFR1,2,3 and VEGFD), proliferation and invasion (IL2, IL3, IL3RB, TIMP1,TIMP2, HGFR, IGFBP2, IGF1, IGF2, MMP2, MMP3 and Pro-MMP2). Notably, TCK1/CXCL7/macrophage derived growth factor/neutrophil activating peptide2 was observed as most significantly decreased cytokine in CSF1R inhibitor treatment compared to vehicle and other groups. In addition to TCK1, other decreased factors in GW2580 group, were elevated in vatalanib group (Fig. 3A and supplementary fig. 4). Since, we used mouse protein array; therefore, it is expected that tumor infiltrated BMDCs were sole source of all cytokines. Network analysis showed that ERK map kinase or P44/42 is a central and critical molecule for most of cytokines evaluated in protein array (Fig. 3C). Further, we tested if drug treated tumor tissues have dysregulated ERK map kinase pathway. Surprisingly, phosphorylated or active form (p-ERK) of ERK was increased after vatalanib treatment and decreased in GW2580 and GW2580+ vatalanib groups compared to vehicle (Fig. 3D).
Molecular networks in GBM TME(A) Heatmap showing expression intensity of cytokine molecules following 96-plex protein array (Ray Biotech. These molecules were elevated in following order of treatment groups; vatalanib > GW2580+ vatalanib > GW2580. White star showing decreased expression of TCK1 or CCL7 following CSF1R blockade (B) Box-plot showing overall expression intensity changes among groups. GW2580 group showed significantly decreased expression intensity of molecules compared to other groups. (C) Network analysis showing ERK map kinase or P44/42 is a central molecule for cytokines (left panel). (D) Western blot using glioma tumor protein lysate showing increased expression of phosphorylated or active form (p-ERK) of ERK in vatalanib treatment and decreased in GW2580 and GW2580+ vatalanib groups (upper and lower panels). Shown is one of the two experiments performed. Quantitative data is expressed in mean ±SD. * P<0.05 and **P<0.01.
3.5 Association of p-ERK expression with CXCL7 expression and proliferation
Immunofluorescence staining was performed to analyze the p-ERK expression on tumor tissues in all treated groups. We noticed expression of p-ERK in both tumor cells and tumor recruited GFP+ cells (more in tumor cells) in vehicle and vatalanib groups. Phospho-ERK expression in GFP+ cells is shown by white arrow heads. Please note significant decrease of GFP+ cells in both anti-CSF1R and combination with vatalanib and thus, decreased p-ERK expression in GFP+ cells (Fig. 4A). Intracellular staining of CXCL7 on tumor tissues was analyzed with all treated groups. We noticed expression of CXCL7 in both tumor cells and tumor recruited GFP+ cells (more in tumor cells) in vehicle and vatalanib groups. CXCL7 expression in GFP+ cells is shown by white arrow heads. Significant decrease of GFP+ cells was seen in both anti-CSF1R and combination with vatalanib and thus, decreased CXCL7 expression in GFP+ cells (Fig. 4B). Phospho-ERK has been known to involve in mitogen associated protein kinase (MAPK) pathway regulated proliferation hallmarks. We evaluated if p-ERK expression correlates with Ki67 proliferation data. Immunofluorescence staining of Ki67 on tumor tissues was analyzed in all treated groups. We noticed expression of Ki67 in both tumor cells and tumor recruited GFP+ cells (more in tumor cells) in vehicle and vatalanib groups. Ki67 expression in GFP+ cells is shown by white arrow heads. We found decreased GFP+ cells in both anti-CSF1R and combination with vatalanib and thus, decreased Ki67 expression in GFP+ cells (Fig. 4C).
Phospho-ERK expression is correlated with CXCL7 expression and proliferation(A) Immunofluorescence analysis showing p-ERK expression in treated tumor tissues. Phospho-ERK expression was seen in both tumor cells and tumor recruited GFP+ cells in vehicle and vatalanib groups. Phospho-ERK expression in GFP+ cells is shown by white arrow heads. Anti-CSF1R and combination with vatalanib significant decreased of GFP+ cells in tumor and thus, decreased p-ERK expression in GFP+ cells. (B) Intracellular staining of CXCL7 on tumor tissues showing expression of CXCL7 in both tumor cells and tumor recruited GFP+ cells in vehicle and vatalanib groups. CXCL7 expression in GFP+ cells is shown by white arrow heads. GW2580 and its combination with vatalanib decreased GFP+ cells and thus, decreased CXCL7 expression in GFP+ cells. (C) Similarly, immunofluorescence staining of Ki67 on tumor tissues noticed expression of Ki67 in both tumor cells and tumor recruited GFP+ cells in vehicle and vatalanib groups. Ki67 expression in GFP+ cells is shown by white arrow heads. We found decreased GFP+ cells in both anti-CSF1R and combination with vatalanib and thus, decreased Ki67 expression in GFP+ cells.
4 DISCUSSION
Since, GBM is hypervascular in nature, different drugs e.g. vatalanib, cediranib, sunitinib, etc. have been used against VEGF-VEGFR pathway to control abnormal angiogenesis in clinical trials [14, 21, 50, 62-64]. Regrettably, benefits of AATs are at best transitory, and this period of clinical benefit is followed by restoration of tumor growth [49]. Evidence of relapse to progressive tumor growth reflects development of resistance to AATs [6, 11]. In our previous study (under review in journal), we established a novel chimeric mouse model to study the contribution of BMDCs in tumor development and therapeutic resistance in GBM. We noticed increased infiltration of GFP+ BMDCs in TME and increased tumor growth as shown by MRI data following vatalanib, nintedanib and CXCR4 antagonist (AMD 3100). Treatments induced polarization of GFP+ BMDCs to the immunosuppressive myeloid cells. Therefore, in the present study, we intended to deplete the myeloid cells in the GBM TME. Similar to our report, other study reported that AATs increased infiltration in myeloid populations in the tumor bulk and in the infiltrative regions of the recurrent GBM (Supplementary fig 2) [45], which indicated requirement of myeloid cell inhibition in GBM.
So far, many myeloid inhibitors have been developed and have shown minimal to broad effect on myeloid cells [7, 13, 32, 48, 52, 60, 76]. GW2580, a selective small molecule kinase inhibitor of CSF1R, is an orally bioavailable competitive inhibitor of adenosine triphosphate binding [17]. GW2580 completely prevented CSF1R dependent growth of macrophages in vitro and in vivo at therapeutically relevant doses [17]. In vitro kinase assays with GW2580 demonstrated at least 100-fold selectivity for CSF1R compared with approximately 300 other structurally related and unrelated kinases.[39] Therefore, GW2580 have great therapeutic potential due to high selectivity and specificity for investigating the contributions of CSF1R signaling in the recruitment and polarization of different myeloid subsets in vivo.
In the present study, GW2580 was tested to inhibit CSF1R+ myeloid cells, which decreased tumor growth by decreasing GFP+ BMDCs acquiring myeloid cells(Gr1+ CD11b+ and F4/80+), angiogenic cells (CD202b+ and CD309+) and immune suppressive signature (IDO) [80]. Interestingly, expression if angiogenic markers on GFP+ cells was associated with the CD31+ vasculature in tumor tissues. In addition, mouse specific protein array data showed decreased expression of vasculature associated markers such as basic fibroblast growth factor (bFGF), vascular cell adhesion molecule 1(VCAM1), E-selectin or CD62 or endothelialleukocyte adhesion molecule 1 (ELAM-1), intercellular Adhesion Molecule 1 (ICAM1) or CD54, vascular endothelial growth factor receptor 2 (VEGFR2) and VEGFR3 in GW2580 group. Altogether, data suggests that anti-CSF1R has profound effect on bone marrow cell contributed (GFP+) vasculature development (vasculogenesis). Previously, CSF1R inhibition using BLZ945 significantly increased survival and regressed established tumors. CSF1R blockade additionally slowed intracranial growth of patient-derived glioma xenografts. Surprisingly, TAMs were not depleted in treated mice. Instead, glioma-secreted factors, including GMCSF and IFN-γ, facilitated TAM survival in the context of CSF1R inhibition. [59]. Our study suggested that GW2580 based CSF1R targeting has better myeloid inhibition compared to BLZ945 drug, which is further confirmed by study done by Xu et al. [77]. Authors showed that CSF1R signaling blockade stanched tumor-infiltrating myeloid cells and improved the efficacy of radiotherapy in prostate cancer [77].
Additionally, to evaluate any changes in cytotoxic T cell functions, we estimated cytokines involved in the T cell function using mouse specific cytokine array, for example, T cell activation (IL1), potent T cell growth factor (IL2), immuno-suppressive cytokine involved in T cell regulatory functions (IL10), helper T cell secreted cytokines that promote proliferation and differentiation of activated B cells (IL4 and IL5), cytokine involved in Th17 cell development (IL6) and Th17 secreted signature cytokine (IL17). We observed significant decreased IL2 cytokine in anti-CSF1R group and trend of increased IL6 and IL10 in anti-CSF1R+ vatalanib group compared to vehicle. This suggests that anti-CSF1R is better that combination treatment with vatalanib. Together, we believe that anti-CSF1R works profoundly through myeloid lineage with modest effect on T cell function.
We observed that combined treatments of GW2580 and vatalanib decreased tumor growth but promoted GFP+ cells to acquire CD133 phenotype in TME. CD133 phenotype has been shown to be associated with cancer stem cell and invasiveness properties [8]. This data clearly confirmed disadvantages associated with VEGFR2 targeting as shown by previous studies [21, 58]. It is possible that decreased tumor growth by combined treatment is transient and tumor could regress due to increased CD133+ cells. However, survival studies are required to test this possibility. Our lab has previously observed similar data with sunitinib treatment against rat glioma, where sunitinib targeted both VEGFR2 and CSFR1. Sunitinib-treated animals showed significantly higher migration of the invasive cells [6]. Similar to our findings, other authors reported that angiogenesis inhibitor anti VEGFR2 antibody (DC101) induces morbidity, when combined with irradiation. In addition, combination therapy using GW2580 with an anti VEGFR2 antibody (DC101) synergistically suppressed tumor growth and severely impaired tumor angiogenesis along with reverting TAM mediated antiangiogenic compensatory mechanism involving MMP-9[58]. We did not observed any decrease in myeloid and endothelial signature cells. There could be two reasons for this differential effect, firstly, our combination treatment started after the day 8 of the tumor implantation that allowed the tumor growth with microenvironment establishment. Our treatment schedule reflected clinical condition where glioma is diagnosed at later stage. On the other hand, study done by Priceman et al. included same day treatment and tumor implantation with no sign of tumor establishment. Secondly, we cannot ignore the differential effects of antibody based and small molecule based protein kinase inhibitor targeting[37].
Several mechanisms have been known to regulate mobilization and recruitment immature myeloid cells into the TME, e.g. IL17 induced expression of GCSF through NF-κB and ERK signaling helped homing of myeloid cells to the tumor [16]. Bv8 modulated mobilization of MDSCs from BM to the tumor and promoted angiogenesis [71]. MDSC can be produced in BM in response to tumor derived factors i.e. GCSF, IL6, GMCSF, IL1β, PGE2 and TNFα, and were recruited to tumor site by CXCL12 and CXCL5 [65]. TGFβ signaling in BMDCs is important and recruits MDSCs via CCL2 in TME [26]. CEACAM1 is identified as negative regulator of myeloid cell expansion and recruitment by inhibiting GCSF-Bv8 axis [46]. Similarly, TIMP2 was shown to down regulate expression of immunosuppressive genes controlling MDSC growth such as IL10, IL13, IL11 and chemokine ligand (CCL5/RANTES), and increased IFN-γ and decreased CD40L [30]. In the present study, we performed 96-plex cytokine array to understand detailed molecular mechanisms and found aforementioned molecules and several other novel molecules involved in inflammation, angiogenesis, invasion and proliferation were downregulated after CSF1R inhibition. TCK1 or CXCL7, a potent chemoattractant and activator of neutrophils, was observed as most significantly decreased cytokine in CSF1R inhibition compared to vehicle and other groups (Fig 3A, shown by star). Since, infiltration of GFP+ cells were significantly decreased after CSF1R blockade (Figure 1C), we believe that decreased CXCL7 is associated with number of GFP+ cells in the TME. Notably, expression of CXCL7 reappeared in GW2580 and vatalanib combined treated tumors, which was associated with the increased GFP+ infiltration in tumor compared with GW2580 only group (Figure 1C).
Biological significance of CXCL7 has been studied before, for example, CXCL7 co-localized with blood vessels, platelets and CD68+ macrophages. CXCL7 expression in synovial macrophages was associated with the inflamed rheumatoid arthritis [79]. In breast tumor, CXCL7, accompanied with IL6 has been shown to modulate human bone marrow mesenchymal cells traffic from the bone marrow to accelerate metastasis at distant sites by expanding the cancer stem cell population [44]. Authors showed that granulocyte recruitment in tumor was dependent on the secretion of CXCL5 and CXCL7 chemokines by platelets upon contact with tumor cells. Blockade of the CXCL5/7 receptor CXCR2, or transient depletion of either platelets or granulocytes prevents the formation of early metastatic niches and significantly reduces metastatic seeding and progression [42]. Proangiogenic and proinflammatory cytokine CXCL7 was an independent prognostic factor for overall survival of clear cell renal cell carcinoma (ccRCC). Inhibition of CXCL7 receptors CXCR1 and CXCR2 was sufficient to inhibit endothelial cell proliferation and ccRCC growth [29]. CXCR2 axis plays a key role in the recruitment of tumor-promoting neutrophils [61] and CXCR2-expressing MDSCs are essential to promote tumorigenesis [40]. In addition, we showed ERK MAPK pathway was involved in cytokine network regulation, which is associated with the CSF1R pathway in macrophages [36]. We believe that decreased p-ERK expression in GFP+ cells is associated with decreased cytokine expression (mouse specific) in GW2580 treated tumors. Moreover, Phospho-ERK has been known to involve in mitogen associated protein kinase (MAPK) pathway regulated proliferation hallmarks. We observed that p-ERK expression was correlated with Ki67 proliferation data. In addition, please note that mouse specific protein array data showed decreased expression of proliferation associated proteins such as IL2, IL3, IL3RB, TIMP1,TIMP2, HGFR, IGFBP2, IGF1, IGF2, MMP2, MMP3 and Pro-MMP2 in GW2580 group. Altogether, data suggests that anti-CSF1R has profound effect on bone marrow cell contributed (GFP+) proliferation events in GBM development.
In conclusion, current study suggested following, (1) myeloid cells are key players in GBM development and AAT resistance (2) since CSF1R inhibition decreased cytokines (e.g.CXCL7) involved in inflammation, angiogenesis, invasion and proliferation; therefore, myeloid cells are major targets for therapeutic targeting in GBM. Clinical trials involving AATs were failed so far; therefore, our study for the first time provided a detailed investigation of altered myeloid cells and cytokine networks, sole cause of therapeutic resistance in GBM (Fig. 5). Nonetheless, survival studies are required to assess the beneficial effect of CSF1R blockade.
Schematic representation of GBM TMEOrthotropic implantation of U251 cells in chimeric mouse initiated migration and homing of GFP+ BMDCs from BM to the tumor. GFP+ cells were found more concentrated at periphery or invasive front of tumor. These infiltrated GFP+ cells acquired myeloid phenotype under influence of tumor cell secreted factors and promoted tumor growth. Blocking CSF1R+ myeloid cells significantly reduced infiltration of GFP+ cells to the tumor and inhibited tumor growth by decreasing CXCL7 expression inflammation, angiogenesis, invasiveness and proliferation via inhibited myeloid cells.
- Established novel chimeric mouse, which uses as low as 5×106 GFP+ BM cells in athymic nude mice and engrafted >70% GFP+ cells within 14 days.
- Orthotopic U251 tumor growth in chimeric mice recapitulated BM mediated vasculogenic in vivo model.
- VEGFR2 blockade increased myeloid and angiogenic cell signatures and thus increased vasculature.
- CSF1R blockade decreased myeloid and angiogenic cell signatures, and thus decreased vasculature.
- CSF1R blockade identified CXCL7 as a mediator of antiangiogenic resistance.
- ERK pathway is critical in CXCL7 mediated myeloid cell regulation.
Supplementary Material
1
9
10
2
3
4
5
6
7
8
ACKNOWLEDGEMENT
Authors thank GRU cancer center small animal imaging core facility for finishing project on timely manner. Authors thank Dr. Hasan Korkaya for providing 96-plex cytokine array (Ray Biotech).
FUNDING
Supported by the National Institutes of Health (NIH) grants R01CA160216 and R01CA172048.
Cd Library Dc101 Software Engineering
Footnotes
Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
REFERENCES
You may think that software media players like iTunes, Windows Media Player, etc., can find and download all the album art you need for your digital music library. However, there are times when you'll need to look further afield in order to successfully populate your music collection with the right CD covers.
You may, for example, have a digital music collection that is mainly made up of a lot of old analog recordings that you have—digitized vinyl records and cassette tapes, for example. Then there's rare compilations, bootleg recordings, and promotional material—album art for these types of audio collections are almost impossible to find using common methods that automatically add metadata tags; MP3 tagging software and music management programs for instance that have built-in ID3 tools.
To help you with this task, take a look at the following list (in no particular order) which showcases some of the best resources on the Internet for finding cover art for your digital music library.
Discogs
What We Like
Large selection.
Sort by genre, style, and format.
Sign up quickly via Facebook or Google.
Rare music selections.
What We Don't Like
Registration required.
Somewhat cluttered album pages.
Discogs is one of the largest online databases for audio. This rich audio catalog resource can be particularly useful for non-mainstream recordings where software media players such as iTunes or Windows Media Player might not be able to find the correct artwork. If you've got hard-to-find commercial releases, bootlegs, white label (promo) material, etc., then you might be able to source the correct album art using Discogs.
The website is easy to use for finding album covers not only for digital music releases but for older mediums, too, like vinyl records, CDs, etc. For digital music, you can also fine-tune your search with a handy filtering option that can be used to only display certain audio formats like AAC, MP3, etc.
Musicbrainz
What We Like
Large mix of search options.
Can download entire database for free.
Helpful documentation.
What We Don't Like
Outdated website design.
Confusing site to navigate.
Musicbrainz is another online audio database that has a huge catalog of music information with included artwork. It was originally conceived as an alternative to CDDB (short for Compact Disc Database) but has now been developed into an online encyclopedia of music that sports a lot more information on artists and albums than simple CD metadata does. For instance, searching for your favorite artist will usually yield information such as all albums released by them (including compilations), audio formats, music labels, background information (relationships to others), and the all-important cover art!