Literature DB >> 22184289

Cord blood stem cells revert glioma stem cell EMT by down regulating transcriptional activation of Sox2 and Twist1.

Kiran Kumar Velpula1, Venkata Ramesh Dasari, Andrew J Tsung, Dzung H Dinh, Jasti S Rao.   

Abstract

The dynamic nature of cancer stem cells that underlie metastasis or their ability to switch between different cellular identities, as in EMT and MET, has profound implications for cancer therapy. The functional relationship between molecules involved in cancer cell stemness and metastasis is not clear. In this regard, our studies on hGBM tissue grade IV specimens showed significant expression of Twist1 and Sox2, known mesenchymal and stemness related markers, respectively, indicating their association with glial tumor genesis and metastasis. The glioma stem cells obtained from CD133+ cells demonstrated increased expression of Twist1 and Sox2 accompanied by significant increase in the mesenchymal markers such as N-cadherin, vimentin and β-catenin. Our studies on glioma stem cells treatment with human umbilical cord blood derived- mesenchymal stem cells, showed down regulation of Twist1 and Sox2 proteins, apart from other mesenchymal stem cell markers. Based on the in vitro experiments and in vivo intracranial xenograft mouse model studies, we elucidated the potential therapeutic role of hUCBSC in suppressing glioma cancer stemness by the induction of MET.

Entities:  

Mesh:

Substances:

Year:  2011        PMID: 22184289      PMCID: PMC3282065          DOI: 10.18632/oncotarget.367

Source DB:  PubMed          Journal:  Oncotarget        ISSN: 1949-2553


INTRODUCTION

Glioblastoma multiforme represents one of the most frequent brain tumors with high cellular complexity. Cancer stem cells (CSC) are cells within tumor cells that possess self-renewal capacity, causing heterogeneous lineages of cancer cells with distinct characteristics of tumor-initiating, invasion, neurosphere formation and differentiation [1-5] contributing tumor maintenance and recurrence [6, 7]. Tumor recurrence and metastasis are identified as the primary cause of treatment failure in several glioblastoma patients. It is thus essential to understand the molecular mechanisms that restore the GSCs balance towards differentiation. Several investigators believed induction of EMT (epithelial to mesenchymal transition) is crucial for carcinogenesis and cancer progression [8, 9]. EMT of cancer cells leads to malignant carcinomas with an increased invasive or metastatic phenotype leading to tumor progression [10]. At the molecular level, EMT is interpreted by the down-regulation of epithelial markers and transcriptional induction of mesenchymal markers [11]. Transcription factors like Twist1 and SNAIL are known to be associated with cellular EMT, promoting cadherin switching. Twist1 is reported to be a master regulator of EMT and metastasis in breast, gastric, hepatocellular, prostate, and brain cancers including glioblastoma [11-21]. Further, growing evidences suggest that pathways regulating self-renewal are deregulated in cancer stem cells resulting in the continuous expansion of self-renewing cancer cells and tumor formation [22]. Sox2, commonly expressed by glioma cells and stem cells of the embryonic and adult brain [23] maintains proliferative potential of neural precursor cells, glioma-initiating cells and drives oncogenesis- associated aggressive tumor phenotype. Downregulation of Sox2 blocks proliferation and induces neuronal differentiation [24]. Variable percentage expression of Sox2 has been observed in gastric cancer, breast cancer, pancreatic cancer and glioblastoma [25], suggesting its relevance to cancer cell aberrant growth. Various experiments conducted on Sox2 have been shown to regulate the expression of other essential genes as well as its own transcription by positive feedback loops [26, 27]. Interestingly, several transcriptional networks with vital functions in neural stem cell populations appear to be co-expressed in gliomas. Studies are required to characterize the expression of different stem cell regulatory components in gliomas and to better understand their function in cancer progression. In exploring new therapeutic strategies to counter glioblastoma progression, our previous studies have shown that human umbilical cord blood stem cells (hUCBSC) provides a promising platform to treat glioblastoma [28-30]. We investigated the co-regulatory role of Sox2 and Twist1 in GSC and their combinatorial effect on stemness maintenance and tumor progression. We also showed that restoration of glioma cancer cell EMT to MET pathway by hUCBSC treatment was accompanied by de-differentiation and dissemination of the tumor cells resulting in a less aggressive tumor phenotype. Our studies are focused on the role of hUCBSC in retarding GSC invasion and metastasis by suppressing the transcriptional activity and association of Sox2 and Twist1. To the best of our knowledge, this is the first study to discover that hUCBSC treatment per se might be significant enough to promote MET in vivo. However, further studies are needed to elucidate the underlying mechanisms.

RESULTS

Sox2 and Twist1 are expressed in human glioblastoma (hGBM) patient-derived specimens

Although earlier studies have indicated the independent expression of Sox2 [31-33] and Twist1 [12-14, 16, 18, 19] in different cancers including glioma [17, 25], the functional relationship of Sox2 and Twist1 together in GSC remains unknown. To evaluate this relationship, we examined seven tissue samples of grade IV gliomas (from surgical biopsy specimens) labeled with antibodies against Sox2 and Twist1. We detected high Sox2 and Twist1 expression levels in these samples with differential expression in non-tumor areas, suggesting their expression in cancerous tissues and not in normal brain tissue (Figs. 1A-1B). (Five hGBMs 2,4,5,6 and 7 of seven glioma specimens expressed variable levels of both Sox2 and Twist1). We next studied the global expression of these two regulatory molecules in different cancers. Twist1 expression was detected in medulloblastoma, meningioma, lung cancer, sarcoid lung cancer, and pancreatic cancer but not in prostate cancer (Supp. Fig. S1A) while Sox2 expression was seen in lung cancer, sarcoid lung cancer and meningioma but not in medulloblastoma, pancreatic and prostrate cancer (Supp. Fig. S1B). We also observed expression of Sox2 and Twist1 in U251 and 5310 parental cell lines (Figs. 1C-1D). Over expression of Sox2 and Twist1 in glioblastoma enhances tumor cell invasion and metastasis, thus supporting roles for Twist1 and Sox2 in glial tumorigenesis and progression [34, 35]. Our present data indicate that both Sox2 and Twist1 expression escorts glial tumorigenesis and that their higher expression levels may be associated with grade IV glial tumors.
Figure 1

Primary glioblastoma tumors express mesenchymal makers Twist1 and Sox2 phenotype

(A-B) Representative immunohistochemical staining of Twist1 and Sox2 in human GBM patient specimens and normal brain. Bar = 200 μm. Strong positivity of Twist1 and Sox2 was detected in five specimens (hGBMs 2,4,5,6 and 7). Negative staining was seen in normal human brain samples. (C-D) Specific expression of Twist1 and Sox2 in U251 and 5310 mouse xenografts models. Strong positivity of Twist1 and Sox2 was detected by H&E staining. Bar = 200 μm;. (E) GBM tissue specimen lysates for western blot analysis were prepared as described previously [30]. In vivo expression was studied by loading equal amounts of protein (40 μg) from tissue lysates onto 10-12% SDS-PAGE gels. (F) Semi-quantitative RT-PCR analysis was done to determine the expression of various stem cell markers in human GBM-derived biopsies. All determinations were done in duplicate.

Primary glioblastoma tumors express mesenchymal makers Twist1 and Sox2 phenotype

(A-B) Representative immunohistochemical staining of Twist1 and Sox2 in human GBM patient specimens and normal brain. Bar = 200 μm. Strong positivity of Twist1 and Sox2 was detected in five specimens (hGBMs 2,4,5,6 and 7). Negative staining was seen in normal human brain samples. (C-D) Specific expression of Twist1 and Sox2 in U251 and 5310 mouse xenografts models. Strong positivity of Twist1 and Sox2 was detected by H&E staining. Bar = 200 μm;. (E) GBM tissue specimen lysates for western blot analysis were prepared as described previously [30]. In vivo expression was studied by loading equal amounts of protein (40 μg) from tissue lysates onto 10-12% SDS-PAGE gels. (F) Semi-quantitative RT-PCR analysis was done to determine the expression of various stem cell markers in human GBM-derived biopsies. All determinations were done in duplicate. To further characterize Sox2 and Twist1 association with malignant progression, we examined changes in protein and mRNA levels in hGBM samples. Western blotting along with RT-PCR analysis of these hGBM specimens demonstrated correlative expression of Sox2 and Twist1 (Figs. 1E-1F). Nonetheless, both Twist1 and Sox2 levels varied considerably among the seven hGBM samples studied but synergistic expression was observed between Twist1 and Sox2. The degree of malignancy of these GBM tumors may in part be reflected by the expression levels of both Twist1 and Sox2. Epithelial malignancy triggers the transition from an epithelial to a mesenchymal cell and is characterized by the loss of E-cadherin and gain of N-cadherin or vimentin expression. Immunoblot analysis showed increased expression of mesenchymal markers such as N-cadherin, Vimentin and decrease of E-cadherin. GSC markers, including CD133, CD44, Msi-1 and Stro-1, showed differential expression in the hGBM samples when compared to human control brain samples (Fig. 1E). These findings from western blot analysis were confirmed using RT-PCR to identify the expression of Sox2 and Twist1 along with other glioma stem cell markers in a panel of seven hGBM specimens (Fig. 1F).

GSC characterization

As shown by immunocytochemistry, GSCs obtained from U251, U87, 4910 and 5310 cells were positive for the GSC markers of proliferation, glial differentiation, and neuronal differentiation including CD133, CD44, GFAP, Ki-67, Msi-1, Nestin, Sox2, Stro-1 and Tuj1 (Fig. 2A), while the non-GSCs evaluated for the expression of the afore-mentioned markers showed reduced expression levels (Supp. Fig. S3D). To further assess the in vitro expression of Twist1 and Sox2 with regard to other GSC markers, we cultured CD133-sorted GSCs from U251, U87, 4910 and 5310 cell lines as neurospheres. Approximately 87.57% (U251 GSC), 69.58% (U87 GSC), 83.75% (4910 GSC) and 28.49% (5310 GSC) of CD133-labeled/sorted cells were positive for CD44; 80.5% (U251 GSC), 76.36% (U87 GSC), 87.8% (4910 GSC) and 91% (5310 GSC) of cells were positive for Stro1; 71.2% (U251 GSC), 64% (U87 GSC), 93.02% (4910 GSC) and 90% (5310 GSC) were positive for Sox2; and 73% (U251 GSC), 70.58% (U87 GSC), 73.75% (4910 GSC) and 82.49% (5310 GSC) were positive for Msi1 (Fig. 2B and Supp. Fig. S2). The GSCs obtained from U251, U87, 4910 and 5310 showed about 70-90% increase in the invasive potential while the non-GSCs showed 40-60% invasive potential by matrigel assay (Supp. Fig. S3A). In another experiment, hUCBSC co-culture challenged GSC- and their respective non GSC-induced microvessel formation (Supp. Fig. S3B). Alternatively, RT-PCR and western blot analysis confirmed abundant expression of aforementioned stem cell markers compared to their parent counterparts (Figs. 2C-2D). Earlier studies have shown that CD133 expression in tumor cells drives tumor progression [36], but the expression of other genes related to stemness on CD133-positive cells is not clearly understood. The increase in the expression of stemness related markers in GSCs (Figs. 2C-2D) was correlated with the increase in the expression of Sox2, Twist1, N-cadherin and Vimentin (Figs. 2C & 2E).
Figure 2

Evaluation of stem cell-like characteristics of U251, U87, 4910 and 5310 neurospheres

(A) Immunocytochemical fluorescence analysis of U251, U87, 4910 and 5310 neurospheres (passage 10). The majority of neurospheres are immunopositive for CD133, CD44, GFAP, Ki-67, Msi-1, Nestin, Sox2, Stro-1 and Tuj-1. Inset pictures show DAPI (blue) staining. Bar = 100 μm. (B) Double positive co-expression of CD133+ GSCs isolated from U251, U87, 4910 and 5310 neurospheres were studied by flow cytometry. CD133-sorted cells from all four GSCs were labeled individually with CD44, Stro-1, Sox2 and Msi-1 with their respective antibodies. In all the experiments, CD133 was conjugated with Alexa Fluor-594 (red), and CD44, Stro-1, Sox2 and Msi-1 was conjugated with Alexa Fluor-488 (green). The experiment was performed in duplicate, and error bars represent standard deviation. (C) Expression levels of various glioma stem cell markers and other molecules were analyzed in U251, U87, 4910 and 5310 neurospheres in comparison to their respective non-GSC cells using semi-quantitative RT-PCR. β-actin served as an internal control for equal loading of the PCR products. Results presented are representative images of three independent experiments (n=3). (D-E) Western blot analysis of various neuronal and EMT makers in U251, U87, 4910 and 5310 neurospheres. Data presented here are a representation of three individual experiments (n=3). All cell lysates were obtained from cells at passage 10.

Evaluation of stem cell-like characteristics of U251, U87, 4910 and 5310 neurospheres

(A) Immunocytochemical fluorescence analysis of U251, U87, 4910 and 5310 neurospheres (passage 10). The majority of neurospheres are immunopositive for CD133, CD44, GFAP, Ki-67, Msi-1, Nestin, Sox2, Stro-1 and Tuj-1. Inset pictures show DAPI (blue) staining. Bar = 100 μm. (B) Double positive co-expression of CD133+ GSCs isolated from U251, U87, 4910 and 5310 neurospheres were studied by flow cytometry. CD133-sorted cells from all four GSCs were labeled individually with CD44, Stro-1, Sox2 and Msi-1 with their respective antibodies. In all the experiments, CD133 was conjugated with Alexa Fluor-594 (red), and CD44, Stro-1, Sox2 and Msi-1 was conjugated with Alexa Fluor-488 (green). The experiment was performed in duplicate, and error bars represent standard deviation. (C) Expression levels of various glioma stem cell markers and other molecules were analyzed in U251, U87, 4910 and 5310 neurospheres in comparison to their respective non-GSC cells using semi-quantitative RT-PCR. β-actin served as an internal control for equal loading of the PCR products. Results presented are representative images of three independent experiments (n=3). (D-E) Western blot analysis of various neuronal and EMT makers in U251, U87, 4910 and 5310 neurospheres. Data presented here are a representation of three individual experiments (n=3). All cell lysates were obtained from cells at passage 10.

Synergistic expression of Twist1 and Sox2 promotes EMT

Previous reports indicate the role of Twist1 in increased tumor metastatic ability leading to the epithelial-to-mesenchymal transition [11]. In the present study, we found that down-regulation of either Twist1 or Sox2 in U251, U87, 4910 and 5310 GSC resulted in phenotypic change from neurospheres to differentiated colonies (Fig. 3A).
Figure 3

Analysis of U251, U87, 4910 and 5310 GSCs treated with Twist1- and Sox2- siRNA

DAB immunocytochemistry was done to study the phenotypical changes of the U251, U87, 4910 and 5310 GSC when treated with (A) Twist1 siRNA and Sox2 siRNA. Cell lystates (40 μg protein for each sample) obtained by treating U251, U87, 4910 and 5310 GSCs with (B) siTwist and (C) siSox2 were loaded onto 8-12% SDS-PAGE gels, transferred onto nitrocellulose membranes, and probed with respective antibodies showing MET transition. (D-E) Nuclear extracts of the same samples were subjected to Western blotting (n=3). GSC (+) indicates the lysate obtained from scrambled vector in the immunoblots.

Analysis of U251, U87, 4910 and 5310 GSCs treated with Twist1- and Sox2- siRNA

DAB immunocytochemistry was done to study the phenotypical changes of the U251, U87, 4910 and 5310 GSC when treated with (A) Twist1 siRNA and Sox2 siRNA. Cell lystates (40 μg protein for each sample) obtained by treating U251, U87, 4910 and 5310 GSCs with (B) siTwist and (C) siSox2 were loaded onto 8-12% SDS-PAGE gels, transferred onto nitrocellulose membranes, and probed with respective antibodies showing MET transition. (D-E) Nuclear extracts of the same samples were subjected to Western blotting (n=3). GSC (+) indicates the lysate obtained from scrambled vector in the immunoblots. The increased expression of N-cadherin, vimentin, β-catenin and Sox2 in U251, U87, 4910 and 5310 control GSC when compared to the si-Twist-transfected samples is suggestive of the molecular alterations during mesenchymal transition (Figs. 3A & 3B). Apart from the EMT markers in Twist1 knockdown cells, we also observed a significant decrease in Sox2 levels, which is the first finding of this sort. To study in detail, our hypothesis that Twist1 regulates Sox2 expression, we transfected U251, U87, 4910 and 5310 GSCs with si-Sox2. We observed increased expression levels of E-cadherin but reduced N-cadherin, vimentin, β-catenin and Twist1 expression in the si-Sox2 transfectants (Fig. 3C). Nuclear extracts of U251, U87, 4910 and 5310 GSCs from both si-Twist1 and si-Sox2 transfected samples demonstrated reduced expression levels when compared to their respective control GSCs (Figs. 3D-3E). Collectively, the above data suggests that the expression of both Twist1 and Sox2 in accord may be necessary to maintain the stemness and induce EMT in GSC.

Expression studies of Twist1 and Sox2 in si-Twist, si-Sox2 transfected GSC

To further confirm expression levels of Twist1 and Sox2 in both si-Twist1 and si-Sox2 transfected GSCs, we did FACS analysis. Contour plots were plotted to study frequency distribution of the expression levels of Twist1 and Sox2. si-Twist1 reduced the expression levels of both Twist1 and Sox2, while si-Sox2 reduced both Sox2 and Twist1 expression levels in all four different GSCs studied (Figs. 4A-4D). Overlay of the dot plots, shows reduced expression levels of Sox2 and Twist1 as observed in contour plots to provide clear understanding of the synergistic phenomenon. The percent expression of Twist1 after cells were treated with si-Twist was observed to be reduced to 43%, 40%, 28% and 35% in U251, U87, 4910 and 5310 GSCs, respectively. The percent expression of Sox2 after cells were treated with si-Twist was observed to be reduced to 25%, 16%, 18% and 14% in U251, U87, 4910 and 5310 GSCs, respectively (Figs. 4E-4F). Similarly, the percent expression of Sox2 after cells were treated with si-Sox2 was reduced to 29%, 27.8%, 29% and 31% in U251, U87, 4910 and 5310 GSCs, respectively. The percent expression of Twist after cells were treated with si-Sox2 was observed to be reduced to 18.6%, 19%, 16% and 16% in U251, U87, 4910 and 5310 GSCs, respectively (Figs. 4G-4H). These results are in agreement with our earlier findings that knockdown of Twist1 led to reduced expression of Sox2 and similarly, knockdown of Sox2 decreased Twist1 expression. These results suggest that their combinatorial effect may contribute to the maintenance of stemness and metastatic ability of the GSC.
Figure 4

Effect of siTwist and siSox2 on Twist1 and Sox2 expression using frequency distribution contour plots

Approximately 1×105 U251 (A), U87 (B), 4910 (C), and 5310 (D) cells were transfected with either siTwist or si-Sox2. Samples were harvested after 72 hrs. siTwist1 and siSox2 expression was determined by staining with the mAb to human Twist1, Sox2 or isotype control mouse or goat IgG1, and then followed by a secondary antibody goat anti-mouse IgG or donkey anti-goat IgG conjugated to Alexa Fluor red or green dye. The background was subtracted using different isotypic controls for respective antibodies. In overlays under Figures 4A-4D, purple color represents the negative control, green color indicates the siTwist or siSox2 treatment, orange color represents Sox2 expression, and red color indicates Twist1 expression. Bar graph (E) represents the expression level of Twist1 in siTwist1 knock down (F) represents the expression level of Sox2 in siSox2 knockdown (G) represents the expression level of Sox2 in siTwist1 knock down while (H) represents the percent expression levels of Twist1 in siSox2 knock down as quantified from the above contour plots. Cells were analyzed with a FACSCalibur flow cytometer and CellQuest software.

Effect of siTwist and siSox2 on Twist1 and Sox2 expression using frequency distribution contour plots

Approximately 1×105 U251 (A), U87 (B), 4910 (C), and 5310 (D) cells were transfected with either siTwist or si-Sox2. Samples were harvested after 72 hrs. siTwist1 and siSox2 expression was determined by staining with the mAb to human Twist1, Sox2 or isotype control mouse or goat IgG1, and then followed by a secondary antibody goat anti-mouse IgG or donkey anti-goat IgG conjugated to Alexa Fluor red or green dye. The background was subtracted using different isotypic controls for respective antibodies. In overlays under Figures 4A-4D, purple color represents the negative control, green color indicates the siTwist or siSox2 treatment, orange color represents Sox2 expression, and red color indicates Twist1 expression. Bar graph (E) represents the expression level of Twist1 in siTwist1 knock down (F) represents the expression level of Sox2 in siSox2 knockdown (G) represents the expression level of Sox2 in siTwist1 knock down while (H) represents the percent expression levels of Twist1 in siSox2 knock down as quantified from the above contour plots. Cells were analyzed with a FACSCalibur flow cytometer and CellQuest software.

Treatment of GSC with hUCBSC induces MET

Previously, we have shown that hUCBSC regulates glioma tumor progression [28, 29]. Taking cue from our previous findings, we tried to evaluate the efficacy of hUCBSC against U251, U87, 4910 and 5310 GSCs. To start with, 1×106 GSCs were co-cultured with 1×106 hUCBSC (1:1 ratio) for 72 hrs. We observed that hUCBSC co-culture retrogress the aggressive mesenchymal glioma stem cells to their epithelial mode via differentiation reverting back to the parent phenotype (Fig. 5A). To further analyze the kinetics of hUCBSC-induced differentiation and proliferation, different phases of the cell cycle were determined by FACS analysis after 72 hrs. U251, U87, 4910 and 5310 control GSCs showed 64%, 79.01%, 75.75%, and 75.69% of the G0-G1 phase; 4.3%, 12% 1.89%, and 6.1% of the S phase; and 6.38%, 1.69%, 7.09%, and 6.82% of the G2-M phase, respectively. In contrast, hUCBSC-treated cells showed 34.6%, 29.75%, 43.81%, and 43.82% of the G0-G1 phase; 13.3%, 18% 6.89%, and 8.1% of the S phase; and 41.22%, 40.59%, 21.17%, and 17% of the G2-M phase, respectively (Figs. 5B-5C). To determine if the treatment of hUCBSC resulted in a decrease in the invasive ability of GSCs, we monitored the effect of serum-induced invasion through Matrigel. GSC/ non-GSCs migrated aggressively through Matrigel pores while hUCBSC treatment reduced the invasive potential by around 40-70% (Supp. Fig. S3A). Similarly, angiogenic tube formation was reduced both in GSC/non-GSCs when co-cultured with hUCBSC (Supp. Fig. S3B).
Figure 5

Effect of hUCBSC on GSC

(A) Co-culture experiments showing induction of MET transition by hUCBSC in GSCs. (B) Various stages of cell cycle of GSCs influenced by hUCBSC. In all of the GSCs, cell cycle was arrested at the G2-M stage. (C) Bar graph showing the results of (B). Western blot analysis showing MET pathway molecules in U251, U87, 4910 and 5310 GSCs when co-cultured with hUCBSC (D) in vitro (72 hrs) and (E) in vivo. (F) Semi-quantitative RT-PCR analysis showing MET transition of GSCs.

Effect of hUCBSC on GSC

(A) Co-culture experiments showing induction of MET transition by hUCBSC in GSCs. (B) Various stages of cell cycle of GSCs influenced by hUCBSC. In all of the GSCs, cell cycle was arrested at the G2-M stage. (C) Bar graph showing the results of (B). Western blot analysis showing MET pathway molecules in U251, U87, 4910 and 5310 GSCs when co-cultured with hUCBSC (D) in vitro (72 hrs) and (E) in vivo. (F) Semi-quantitative RT-PCR analysis showing MET transition of GSCs. Since it is understood that Twist1 and Sox2 regulate the metastatic ability of tumor cells and their involvement in stemness preservation, we analyzed the hUCBSC co-culture treated GSC for Sox2 and Twist1 expression. Immunoblotting analysis demonstrate that treatment of U251, U87, 4910 and 5310 GSC with hUCBSC reduced N-cadherin, β-catenin, vimentin, Sox2 and Twist1 expression, with increased E-cadherin levels indicative of MET course (Fig. 5D). To further explore the effect of hUCBSC on GSCs, we looked for changes in cellular functions frequently associated with EMT by hUCBSC treatment. We examined the expression levels of the above mentioned proteins in U251 and 5310 GSC xenografts by western blot analysis; the expression levels of mesenchymal markers were decreased upon hUCBSC treatment (Fig. 5E). RT-PCR results confirmed our above results obtained from western blot analysis (Fig. 5F).

hUCBSC treatment inhibits Sox2 and Twist1 expression in the intracranial GSC induced tumors of mice

To determine the analeptic of hUCBSC on GSC tumor formation in vivo, athymic nude mice were injected with U251 and 5310 GSCs (1×105) in the left frontal lobe. Mice injected with U251 and 5310 GSCs showed detectable tumor characteristics, such as hunch back appearance, paralyzed limbs, severe weight loss and neurologic pain within 12-18 days. After three weeks, an equal number of hUCBSC (1×105) were injected into the right lobe of the mice, already having GSC induced tumors. The hUCBSC treatment showed regression in GSC induced tumors and tumor development. Mice treated with only hUCBSC did not form tumors for 120 days (Data not shown). Pathological examination of GSC induced tumors in the nude mice showed a vast necrotic area associated with the tumor tissues; hUCBSC treatment significantly reduced GSC induced tumor volumes (Supp. Fig. S3C). The tumor volumes were reduced to 65-75% in U251 and 75-90% in 5310 GSC- treated mice upon hUCBSC treatment as evidenced by H&E staining (Fig. 6C). To corroborate our findings and to extend our results from hGBM tissue specimens, we analyzed the expression of both Twist1 and Sox2 by immunohistochemistry. Adjacent sections were used to study the expression of both Twist1 and Sox2. High expression of Twist1 and Sox2 was observed in the necrotic regions of U251 and 5310 GSC induced tumor, while hUCBSC-treated U251 and 5310 xenografts showed reduced or no expression of either Twist1 or Sox2 in individual experiments (Fig 6C).
Figure 6

Twist1 co-expresses with Sox2 in U251, U87, 4910 and 5310 GSCs

(A) Immunostaining analysis from neurospheres generated from GSC-enriched populations (U251, U87, 4910 and 5310 cells) shows Twist1 staining that appears to correlate with high Sox2 expression. Photographs depicted show Twist1 (red) is co-expressed/associated with Sox2 (green). All nuclei are counterstained with DAPI in blue. Scale bars represent 100 μm. (B) Immunostaining of GBM mouse xenografts (U251 and 5310) demonstrate that Twist1 (red) is expressed in the necrotic area and is co-expressed with Sox2 (green). Scale bars represent 50 μm. DAB immunohistochemical analysis of Twist1 and Sox2 expression in U251 GSC and 5310 GSC xenografts in nude mice brains. (C) Representative microphotographs of U251 and 5310 GSCs injected intracranially, which show infiltrative/migratory nature of xenograft-derived tumors. Insets in figures are representative H&E stained sections. Immunohistochemical staining was done using specific antibodies for Twist1 and Sox2 (dark brown in color). Scale bar represents 100 μm.

Twist1 co-expresses with Sox2 in U251, U87, 4910 and 5310 GSCs

(A) Immunostaining analysis from neurospheres generated from GSC-enriched populations (U251, U87, 4910 and 5310 cells) shows Twist1 staining that appears to correlate with high Sox2 expression. Photographs depicted show Twist1 (red) is co-expressed/associated with Sox2 (green). All nuclei are counterstained with DAPI in blue. Scale bars represent 100 μm. (B) Immunostaining of GBM mouse xenografts (U251 and 5310) demonstrate that Twist1 (red) is expressed in the necrotic area and is co-expressed with Sox2 (green). Scale bars represent 50 μm. DAB immunohistochemical analysis of Twist1 and Sox2 expression in U251 GSC and 5310 GSC xenografts in nude mice brains. (C) Representative microphotographs of U251 and 5310 GSCs injected intracranially, which show infiltrative/migratory nature of xenograft-derived tumors. Insets in figures are representative H&E stained sections. Immunohistochemical staining was done using specific antibodies for Twist1 and Sox2 (dark brown in color). Scale bar represents 100 μm. As evidenced by the significant increase in the expression of both Twist1 and Sox2 in the necrotic tumor areas, and their co-regulatory effect in siRNA transfectants, we next asked whether this expression is the co-regulatory mechanism necessitated by the interaction between Twist1 and Sox2. Initial experiments in this direction were targeted at neurospheres at the in vitro level. Control GSCs obtained from U251, U87, 4910 and 5310 showed clear interaction and co-localization between Twist1 and Sox2 while hUCBSC treatment showed decreased expression and co-localization of both these proteins (Fig. 6A). These results prompted us to study co-localization of these proteins in an in vivo xenograft model. U251 and 5310 GSC control tissue samples showed co-localization of Sox2 with Twist1 which was significantly reduced upon hUCBSC treatment (Fig. 6B). The association of transcription factors, Sox2 (responsible for stemness) and Twist1 (responsible for EMT) thus provides a platform to correlate the stem cell pathway with EMT pathway, which can be further studied to address potential targets for glioblastoma cure.

DISCUSSION

Cancer stem cells render abnormal growth, frequency and radio-resistance in malignant glioblastoma. Therapeutic strategies in treating glioblastoma require a better understanding of the mechanisms that promote stem cell maintenance that sustain tumor growth. Elias, et al. [34] proposed neural precursor cells as the likely targets for malignant transformation and showed that deregulation of the self-renewal capacity of normal stem cells resulted in stable augmentation of self-renewing cancer cells and tumor formation. Investigations to study this tumorigenic process will provide a powerful tool to uncover new targets for glioma therapeutic intervention. On par with the available evidence, it is understood that cancer stem cells (CSC) play an important role in tumor initiation, growth, invasion and recurrence. Therefore, identification and study of CSC will provide insights into better understanding of brain tumors and their cellular targets. In this study, we have obtained GSCs from U251, U87, 4910 and 5310 cell lines using the surface protein CD133 and FACS sorting as reported previously [37-39]. In the present study, we show that malignant human grade IV gliomas display significant expression of Twist1 and Sox2. Twist1 expression was observed in hGBMs as well as in medulloblastoma, meningioma, lung and pancreatic cancer specimens, while Sox2 expression was observed in hGBM specimens along with that of meningioma, lung and pancreatic cancers. These results indicate the participation of both Twist1 and Sox2 in the genesis of a wide range of neuroectodermal tumors. Twist1 encourages cancer cell survival and tumor progression by increasing resistance to cytotoxic therapies resulting in increased cellular proliferation [40]. Twist1 is considered to be a key player of EMT, a metastatic cascade in which epithelial cells take over a mesenchymal phenotype [10]. EMT phenotype expresses cancer stem-like cells signatures associated with tumor recurrence and a drug-resistant phenotype [41]. Immunoblotting analysis conducted on the hGBM specimens have shown reduced epithelial E-cadherin expression and upregulated mesenchymal N-cadherin, vimentin, β-catenin and Twist1, suggesting the EMT mode. High expression levels of other stem cell markers such as CD44, Msi-1, Stro-1, and EMT markers, such as Sox2, Twist1, vimentin suggested an invasive phenotype. The synergistic up-regulation of Twist1 and Sox2 expression in hGBM specimens may instigate E-cadherin replacement by N-cadherin expression triggering cadherin switching [42, 43]. Numerous lines of evidences suggest that the Twist1 regulates programmed cell death and induces apoptosis [44, 45]. Moreover, it has been observed that transfection of Bel7402 (highly expressing Twist1) cells with Twist1 shRNA plasmid resulted in suppression of Twist1 expression and subsequent inhibition of cell invasion and migration in vitro, indicating its possible role in tumor cell plasticity [46]. The widespread expression of Twist1 in various cancer phenotypes suggests its action upon a diverse set of downstream target genes, depending upon the tissue, to elicit a variety of cellular responses. The chimeric nature of beta-Helix loop Helix (bHLH) proteins, such as Twist1 are implicated to possibly combine heterodimerically with other bHLH proteins to activate or repress diverse downstream targets generating different types of biological responses [47]. Cimadamore et al., [48] reported that knockdown of Sox2, the SRY (sex-determining region)-box 2 gene, resulted in downregulation of proneural bHLH transcription factors. Twist1-mediated molecular changes in our study provided important insight into its role in understanding mesenchymal change with genes related to EMT being upregulated in U251, U87, 4910 and 5310 GSCs. Since we observed the co-expression of Sox2 along with Twist1, further studies are warranted to examine their role in acquiring an invasive malignant phenotype. Knockdown of Twist1 reduced the expression levels of Sox2, while knockdown of Sox2 reduced the expression levels of Twist1. In both experiments, increased expression of E-cadherin was observed, which suggests that both Twist1 and Sox2 function in harmony and that their concurrence is essential for the maintenance of tumor plasticity. Inhibition of Twist1 or Sox2 expression by their specific siRNAs resulted in differentiation of GSC formation and growth and is suggestive of the above discussed mechanism. Targeted migration of mesenchymal stem cells (MSCs) to tumors makes them a very promising strategy for anti-tumor therapy [49]. Previously, we have demonstrated the potential therapeutic applications of hUCBSC in the treatment of gliomas [28, 29]. In addition, with increased understanding of the mechanistic underpinnings of MSC action, the potential for enhancing MSC-targeted therapies appears promising. This report provides a strong evidence that hUCBSC inhibit the metastatic and invasive potential of U251, U87, 4910 and 5310 GSCs. Co-culture of hUCBSC with these GSCs changed their morphology and clearly reduced their migration and invasion abilities. We also showed that hUCBSC induced MET in U251, U87, 4910 and 5310 GSCs, which might be responsible for inhibition of metastasis-related abilities. It appears that hUCBSC counteract and repress Twist1 expression by increasing E-cadherin levels, a phenomenon generally observed at sites of EMT. Our present data demonstrated the ability of hUCBSC to inhibit GSC phenotypes in co-culture experiments both in vitro and in vivo. Immunoblotting experiments showed a decrease in the expression levels of both Twist1 and Sox2 along with other EMT markers, suggesting that hUCBSC treatment reversed the EMT to MET. Based on all our findings, we propose that hUCBSC inhibit GSC invasion, migration by suppressing key molecules involved in EMT. These actions resulted in induction of a MET state, which is a key step in reversal and treatment of the cellular metastatic process. In conclusion, for the first time, we showed that interactive and co-regulatory role of Twist1 and Sox2 necessary to maintain glioma stemness was inhibited by hUCBSC. Here, hUCBSC was administered to study key aspects in treating the glioma stem cell machinery. We engineered hUCBSC near the tumor site and observed that hUCBSC are effective in regressing tumor growth in the intracranial xenograft mouse model. Based on our findings, we conclude that hUCBSC can be used as a means of therapy and will be of great interest for the clinical application of stem cell-based cancer therapy.

MATERIALS AND METHODS

Tumor samples

Glioblastoma specimens were obtained during autopsies of glioblastoma patients within 24 hrs of death or from patients who underwent surgery at Saint Francis Medical Center (Peoria, IL). All samples were collected under protocols approved by the UICOMP (Peoria, IL) Institutional Review Board. All tumors without prior evidence of progression from a lower-grade tumor were clinically classified as primary glioblastoma.

Culture of hUCBSC

Human umbilical cord blood stem cells were isolated as described previously (See supplemental methods) [29].

Enrichment for glioma stem cells

Human glioma cell lines parental and non-GSC (U251 and U87) and xenografts (4910 and 5310) were cultured using conventional tissue culture media: DMEM/RPMI-1640 (Hyclone, Logan, UT) supplemented with 10% fetal bovine serum (Gibco BRL, Grand Island, NY), 100 units/mL penicillin/streptomycin (Lonza, Walkersville, MD). GBM neurospheres of U87, U251, and xenografts 4910, 5310 were cultured in neurobasal medium supplemented with N2 (1%), B27 (1%), 20 ng/mL of bFGF and EGF, 10 ng/mL of LIF (Millipore, Billerica, MA)[50]. Experiments for the present study were conducted using the neurospheres obtained between the 5-14th passages.

Immuno- cyto- histochemical staining

Experiments were performed to detect the expressi1 on of various stem cell markers CD133, CD44, Msi-1, Stro1, Nestin, Sox2, Ki-67, GFAP and Tuj-1 in both GSC and non-GSCs (see Supplemental methods).

Semi-Quantitative RT-PCR

Total RNA (5 μg) from cells, brains or tumor samples extracted using Qiagen RNeasy mini kit (Valencia, CA) was reverse-transcribed using the Transcriptor First Strand cDNA Synthesis Kit (Roche, Indianapolis, IN). (See Supplementary methods)

CD133 sorting, cell cycle analysis

GSCs obtained from U251, U87, 4910 and 5310 cells dissociated by trituration into single cell suspension were incubated for 60 min at 4°C with anti-CD133 phycoerythrin (PE)-conjugated antibody (clone AC133-PE, mouse IgG1; Miltenyi Biotec, Bergisch Gladbach, Germany) according to the manufacturer's instructions. PE-conjugated mouse IgG1 isotype antibody was used as control. The cells re-suspended in 0.5 mL PBS-BSA were analyzed and sorted by flow cytometry. Acquisition was performed on a FACSCalibur flow cytometer (Becton Dickinson, San Jose, CA) and viable cells were analyzed with CellQuest software. In a separate experiment, the CD133-positive cells were labeled with anti-CD44, anti-Stro-1, anti-Msi-1 and anti-Sox2 in independent experiments using a similar staining protocol. The control for each sample was prepared identically except that an isotype-specific antibody was used. Progression through different cell cycle phases in U251, U87, 4910 and 5310 cells alone and in co-culture with hUCBSC for 72 hrs was monitored by flow cytometric analysis of the DNA content of cell populations stained with propidium iodide (BioSure, Grass Valley, CA).

Matrigel invasion and In vitro angiogenesis assay

The Matrigel invasion assay and tumor conditioned medium-induced microtubule network formation were studied in U251, U87, 4910 and 5310 GSC, non-GSC alone or in co-culture with hUCBSC (See supplemental methods and materials).

siRNA-based knockdown of Twist1 and Sox2

The siRNA targeting human Twist1 (sc-38604), Sox2 (sc-38408) and scrambled negative control RNA (Santa Cruz, CA). were used to transfect U251, U87, 4910 and 5310 GSCs (2.5×105 cells) using 6 μL of Fugene HD reagent (Roche, Indianapolis, IN). The siRNA/Fugene HD complex was added to the GSCs and the transfection was carried out for 48-72 hrs. Nuclear extracts were prepared from these siRNA-transfected samples using a nuclear extraction kit from Panomics, Inc. (Fremont, CA) according to the manufacturer's instructions.

Contour plot FACS analysis

Twist1 and Sox2 expression was measured by intracellular protein staining using antibodies specific for human Twist1 and Sox2 (Santa Cruz Technologies, Santa Cruz, CA). The Twist1 siRNA- or Sox2 siRNA-transfected cells were incubated for 1 hr with gene-specific Twist1 or Sox2. Twist1 expression in si-Sox2-treated samples and Sox2 expression in si-Twist1 treated samples were studied by measuring cellular fluorescence using the FACSCalibur flow cytometer (Becton Dickinson, San Jose, CA).

Western blotting

GBM parental or stem cells alone or in co-culture with hUCBSC for 72 hrs were collected, lysed in RIPA buffer and then resolved via SDS-PAGE (See supplemental methods).

Intracranial administration of GSC and hUCBSC

Glioma stem cells were injected intracerebrally into the left side of the brain of nude mice with 10 μL aliquots of U251 and 5310 GSC (1×105) under isofluorane anesthesia with the aid of a stereotactic frame. Three weeks after tumor implantation, hUCBSC were injected towards the right side of the brain in 1:1 ratio. Six mice from each group were sacrificed by cardiac perfusion with 4% formaldehyde in PBS, and paraffin sections were prepared. Sections were stained with H&E to visualize tumor cells and to examine tumor volume. The sections were blindly reviewed and scored semi-quantitatively for tumor size. The average visible tumor area per section integrated to the number of sections was calculated using the formula 1/6 II (Rmax) × (Rmin)2, where Rmax and Rmin are the maximum and minimum tumor radii, respectively. The calculated tumor volume was compared between controls and treated groups.
Table 1

PCR primers

GeneForward PrimerReverse Primer
CD1335'CCAAGTTCTACCTCATGTTTGG 3'5'ACCAACAGGGAGATTGCAAAGC3'
CD445' AATCCCTGCTACCAATATGGACT 3'5' AGCCTTCAGAATGATTTGGGTC 3'
N-cadherin5'CAACTTGCCAGAAAACTCCAGG 3'5' ATGAAACCGGGCTATCTGCTC 3'
E-cadherin5' CGAGAGCTACACGTTCACGG 3'5' GGCCTTTTGACTGTAATCACACC 3'
Sox25' CACAACTCGGAGATCAGCAA 3'5' CTCCGGGAAGCGTGTACTTA 3'
Twist15'CGGACAAGCTGAGCAAGATT 3'5' CCTTCTCTGGAAACAATGAC 3'
Nestin5' GGCATCCTCACCCTGAAGTA 3'5' CTTGGGGTCCTGAAAGCTG 3'
b-catenin5'GAAACGGCTTTCAGTTGAG 3'5' CTGGCCATATCCACCAGAGT 3'
Vimentin5'GAACGCCAGATGCGTGAAATG3'5'CCAGAGGGAGTGAATCCAGATTA3'
b-actin5'GGCATCCTCACCCTGAAGTA3'5'GGGGTGTTGAAGGTCTCAAA3'
  50 in total

1.  Sox2 is translationally activated by eukaryotic initiation factor 4E in human glioma-initiating cells.

Authors:  Yuqing Ge; Fengbiao Zhou; Hong Chen; Chunhong Cui; Dan Liu; Qiuping Li; Zhiyuan Yang; Guoqiang Wu; Shuhui Sun; Jianxin Gu; Yuanyan Wei; Jianhai Jiang
Journal:  Biochem Biophys Res Commun       Date:  2010-06-09       Impact factor: 3.575

2.  Transcriptional repression of Mad-Max complex by human umbilical cord blood stem cells downregulates extracellular signal-regulated kinase in glioblastoma.

Authors:  Kiran Kumar Velpula; Venkata Ramesh Dasari; Andrew J Tsung; Dzung H Dinh; Jasti S Rao
Journal:  Stem Cells Dev       Date:  2011-11-21       Impact factor: 3.272

Review 3.  Epithelial-mesenchymal transitions in tumour progression.

Authors:  Jean Paul Thiery
Journal:  Nat Rev Cancer       Date:  2002-06       Impact factor: 60.716

4.  DNA methylation markers for prostate cancer with a stem cell twist.

Authors:  Carmen Jerónimo; Manel Esteller
Journal:  Cancer Prev Res (Phila)       Date:  2010-08-10

5.  Snail2 is an essential mediator of Twist1-induced epithelial mesenchymal transition and metastasis.

Authors:  Esmeralda Casas; Jihoon Kim; Andrés Bendesky; Lucila Ohno-Machado; Cecily J Wolfe; Jing Yang
Journal:  Cancer Res       Date:  2011-01-01       Impact factor: 12.701

6.  Twist haploinsufficiency in Saethre-Chotzen syndrome induces calvarial osteoblast apoptosis due to increased TNFalpha expression and caspase-2 activation.

Authors:  Malika Yousfi; Francoise Lasmoles; Vincent El Ghouzzi; Pierre J Marie
Journal:  Hum Mol Genet       Date:  2002-02-15       Impact factor: 6.150

7.  TWIST1 promotes invasion through mesenchymal change in human glioblastoma.

Authors:  Svetlana A Mikheeva; Andrei M Mikheev; Audrey Petit; Richard Beyer; Robert G Oxford; Leila Khorasani; John-Patrick Maxwell; Carlotta A Glackin; Hiroaki Wakimoto; Inés González-Herrero; Isidro Sánchez-García; John R Silber; Philip J Horner; Robert C Rostomily
Journal:  Mol Cancer       Date:  2010-07-20       Impact factor: 27.401

8.  Regulation of glioblastoma progression by cord blood stem cells is mediated by downregulation of cyclin D1.

Authors:  Kiran Kumar Velpula; Venkata Ramesh Dasari; Andrew J Tsung; Christopher S Gondi; Jeffrey D Klopfenstein; Sanjeeva Mohanam; Jasti S Rao
Journal:  PLoS One       Date:  2011-03-24       Impact factor: 3.240

9.  The SOX2 response program in glioblastoma multiforme: an integrated ChIP-seq, expression microarray, and microRNA analysis.

Authors:  Xuefeng Fang; Jae-Geun Yoon; Lisha Li; Wei Yu; Jiaofang Shao; Dasong Hua; Shu Zheng; Leroy Hood; David R Goodlett; Gregory Foltz; Biaoyang Lin
Journal:  BMC Genomics       Date:  2011-01-06       Impact factor: 3.969

Review 10.  Redundant or separate entities?--roles of Twist1 and Twist2 as molecular switches during gene transcription.

Authors:  Hector L Franco; José Casasnovas; José R Rodríguez-Medina; Carmen L Cadilla
Journal:  Nucleic Acids Res       Date:  2010-10-08       Impact factor: 16.971

View more
  36 in total

1.  Autocrine IL-8 promotes F-actin polymerization and mediate mesenchymal transition via ELMO1-NF-κB-Snail signaling in glioma.

Authors:  Baogang Zhang; Lihong Shi; Shijun Lu; Xiuning Sun; Yuqing Liu; Hongli Li; Xuejian Wang; Chunzhen Zhao; Heng Zhang; Ying Wang
Journal:  Cancer Biol Ther       Date:  2015-04-14       Impact factor: 4.742

2.  miR-429 inhibits glioma invasion through BMK1 suppression.

Authors:  Weiyi Chen; Baogang Zhang; Wenjun Guo; Linlin Gao; Lihong Shi; Hongli Li; Shijun Lu; Yuqing Liu; Xiaolong Li
Journal:  J Neurooncol       Date:  2015-08-14       Impact factor: 4.130

3.  Conditioned Medium from Adipose-Derived Stem Cells (ADSCs) Promotes Epithelial-to-Mesenchymal-Like Transition (EMT-Like) in Glioma Cells In vitro.

Authors:  Isabele C Iser; Stefanie M Ceschini; Giovana R Onzi; Ana Paula S Bertoni; Guido Lenz; Márcia R Wink
Journal:  Mol Neurobiol       Date:  2015-12-19       Impact factor: 5.590

4.  CDA-2 induces cell differentiation through suppressing Twist/SLUG signaling via miR-124 in glioma.

Authors:  You-ke Xie; Shao-fen Huo; Gong Zhang; Fu Zhang; Zu-ping Lian; Xiong-lin Tang; Chuan Jin
Journal:  J Neurooncol       Date:  2012-08-24       Impact factor: 4.130

5.  The homing of human cord blood stem cells to sites of inflammation: unfolding mysteries of a novel therapeutic paradigm for glioblastoma multiforme.

Authors:  Kiran Kumar Velpula; Venkata Ramesh Dasari; Jasti S Rao
Journal:  Cell Cycle       Date:  2012-06-15       Impact factor: 4.534

6.  Urokinase-type plasminogen activator receptor (uPAR)-mediated regulation of WNT/β-catenin signaling is enhanced in irradiated medulloblastoma cells.

Authors:  Swapna Asuthkar; Christopher S Gondi; Arun Kumar Nalla; Kiran Kumar Velpula; Bharathi Gorantla; Jasti S Rao
Journal:  J Biol Chem       Date:  2012-04-16       Impact factor: 5.157

7.  Serine Protease Activation Essential for Endothelial-Mesenchymal Transition in Vascular Calcification.

Authors:  Jiayi Yao; Pierre J Guihard; Ana M Blazquez-Medela; Yina Guo; Jeremiah H Moon; Medet Jumabay; Kristina I Boström; Yucheng Yao
Journal:  Circ Res       Date:  2015-08-11       Impact factor: 17.367

8.  The enforced expression of c-Myc in pig fibroblasts triggers mesenchymal-epithelial transition (MET) via F-actin reorganization and RhoA/Rock pathway inactivation.

Authors:  Jun-Wen Shi; Wei Liu; Ting-Ting Zhang; Sheng-Chun Wang; Xiao-Lin Lin; Jing Li; Jun-Shuang Jia; Hong-Fen Sheng; Zhi-Fang Yao; Wen-Tao Zhao; Zun-Lan Zhao; Rao-Ying Xie; Sheng Yang; Fei Gao; Quan-Rong Fan; Meng-Ya Zhang; Min Yue; Jin Yuan; Wei-Wang Gu; Kai-Tai Yao; Dong Xiao
Journal:  Cell Cycle       Date:  2013-03-06       Impact factor: 4.534

9.  Tumor-derived exosomal lnc-Sox2ot promotes EMT and stemness by acting as a ceRNA in pancreatic ductal adenocarcinoma.

Authors:  Zhonghu Li; Peng Jiang; Jie Li; Minjie Peng; Xin Zhao; Xi Zhang; Kai Chen; Yujun Zhang; Hui Liu; Lang Gan; Huaqiang Bi; Ping Zhen; Jin Zhu; Xiaowu Li
Journal:  Oncogene       Date:  2018-04-12       Impact factor: 9.867

10.  Krüppel like factor 6 splice variant 1 (KLF6-SV1) overexpression recruits macrophages to participate in lung cancer metastasis by up-regulating TWIST1.

Authors:  Jian Wang; Xiao Wang; Yawei Wang; Shuguang Li; Xiuwen Wang
Journal:  Cancer Biol Ther       Date:  2018-12-27       Impact factor: 4.742

View more

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