- 1State Key Laboratory of Bioactive Substances and Functions of Natural Medicines, Institute of Materia Medica, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
- 2Beijing Key Laboratory of Non-Clinical Drug Metabolism and PK/PD Study, Institute of Materia Medica, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
- 3Beijing Key Laboratory of New Drug Mechanisms and Pharmacological Evaluation Study, Institute of Materia Medica, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
Glioblastoma multiforme (GBM) is the most common malignant tumor of the central nervous system. Temozolomide (TMZ)–based adjuvant treatment has improved overall survival, but clinical outcomes remain poor; TMZ resistance is one of the main reasons for this. Here, we report a new phosphatidylinositide 3-kinase inhibitor, XH30; this study aimed to assess the antitumor activity of this compound against TMZ-resistant GBM. XH30 inhibited cell proliferation in TMZ-resistant GBM cells (U251/TMZ and T98G) and induced cell cycle arrest in the G1 phase. In an orthotopic mouse model, XH30 suppressed TMZ-resistant tumor growth. XH30 was also shown to enhance TMZ cytotoxicity both in vitro and in vivo. Mechanistically, the synergistic effect of XH30 may be attributed to its repression of the key transcription factor GLI1 via the noncanonical hedgehog signaling pathway. XH30 reversed sonic hedgehog–triggered GLI1 activation and decreased GLI1 activation by insulin-like growth factor 1 via the noncanonical hedgehog signaling pathway. These results indicate that XH30 may represent a novel therapeutic option for TMZ-resistant GBM.
Introduction
Glioblastoma multiforme (GBM) is the most aggressive tumor of the central nervous system in adults (DeAngelis, 2001). Despite recent research focusing on various glioblastoma therapies, the clinical benefit from treatment for people with GBM remains unsatisfactory (Touat et al., 2017; Castresana and Melendez, 2021; Liau, 2021). Mean overall survival is currently estimated to be only 9.7 months after primary treatment.
Temozolomide (TMZ) chemotherapy is still recommended as the standard care for GBM. However, GBM recurs in most people with this tumor type after primary treatment. The development of resistance to TMZ is often the primary limiting factor for treatment success (Bocangel et al., 2002; Lee, 2016; Dymova et al., 2021). Many factors contribute to TMZ resistance, such as overexpression of O6-methylguanine methyltransferase (MGMT), lack of a DNA repair pathway, activation of the hedgehog signaling pathway, the presence of glioma stem cells, and metabolic dysfunction (Happold and Weller, 2015; Perazzoli et al., 2015; Yun et al., 2020; Chien et al., 2021; Zheng et al., 2021). The hedgehog signaling pathway has been a focus of research in this area recently (Melamed et al., 2018; Avery et al., 2021). The glioma-associated oncogene GLI, a zinc finger protein, is a key component in this pathway. Pharmacological inhibition of GLI-1 has been shown to enhance the cytotoxicity of TMZ in GBM and overcome TMZ resistance (Li J. et al., 2016; Ji et al., 2018).
The phosphatidylinositide 3-kinase (PI3K) pathway is frequently overactivated in glioblastoma due to PIK3CA mutations, loss of phosphotase and tensin homolog (PTEN) gene function, and amplification of epidermal growth factor receptor (EGFR) gene expression (Sami and Karsy, 2013; Langhans et al., 2017; Colardo et al., 2021). PI3K is an attractive therapeutic target for glioblastoma (Li X. et al., 2016; Colardo et al., 2021). In animal models, PI3K inhibitors have been shown to have substantial antitumor activity against GBM (Zhao et al., 2017). The effects of PI3K inhibitors, such as paxalisib, have been tested in patients with newly diagnosed GBM (Wen et al., 2019). Preliminary data have shown encouraging survival outcomes in these patients. In addition, PI3K inhibitors have also been shown to enhance TMZ cytotoxicity in GBM via distinct mechanisms, such as downregulation of ATP-binding cassette subfamily E member 1, inhibition of autophagy, promotion of apoptosis, and inhibition of DNA double-strand break repair (Gil del Alcazar et al., 2014; Radoul et al., 2016; Zhang et al., 2018; Zajac et al., 2021). However, the antitumor activity of PI3K inhibitors in patients with TMZ-resistant GBM currently remains unclear (Hainsworth et al., 2019).
Many studies have reported crosstalk between the PI3K and hedgehog signaling pathways (Ranjan and Srivastava, 2017). The PI3K signaling pathway is a crucial non-canonical activator of GLI1 (Riobo et al., 2006; Zhou et al., 2016; Liang et al., 2017; Ranjan and Srivastava, 2017). Activation of this pathway has been found to enhance GLI1 protein stability because the serine/threonine kinase in this pathway, AKT, can extend the half-life of GLI proteins in the cells by alleviating the inhibitory effect of protein kinase A, thus facilitating nuclear translocation (Singh et al., 2017). Meanwhile, PI3K signaling activates GLI1 via its downstream effector, ribosomal S6 kinase (p70S6K) (Wang et al., 2012). Activated p70S6K promotes GLI1 disassociation from suppressor of fused homolog (SUFU) by phosphorylating GLI1 at Ser84 and enhancing GLI1 transcriptional activity. In addition, p70S6K2 has been shown to inhibit glycogen synthase kinase (GSK3) by phosphorylating GLI1 at Ser9, leading to a decrease in GSK3β-mediated GLI1 degradation (Mizuarai et al., 2009).
XH30 is a PI3K inhibitor that can cross the blood–brain barrier (Lin et al., 2018). XH30 has previously been demonstrated to have robust antitumor activity in GBM and brain metastases of lung cancer in vivo (Lin et al., 2018; Ji et al., 2021). In this study, we aimed to assess the capacity of XH30 to inhibit the growth of GBM cells with natural or TMZ-induced drug resistance both in vitro and in vivo, in an orthotopic mouse model. We also explored the underlying mechanisms of the antitumor effects of XH30.
Materials and Methods
Cell Lines
Two cell lines were used in this study. The T98G human glioma tumor cell line was purchased from ATCC (Manassas, VA, United States); this is a naturally TMZ-resistant cell line with elevated levels of MGMT. The U251/TMZ cell line was a gift from Dr. Yuhui Zou of General Hospital of Guangzhou Military Command of People’s Liberation Army, as previously reported (Ji et al., 2018); this line has acquired TMZ resistance, with overactivated hedgehog signaling. Both cell lines were cultured in Dulbecco’s modified eagle medium (Gibco, TX, United States) with 10% (v/v) fetal bovine serum (Gibco), penicillin (100 units/ml), and streptomycin (100 units/ml) in a humidified atmosphere of 5% CO2 at 37°C.
Antibodies and Reagents
XH30 was synthesized in house as described previously (Lin et al., 2018). PI3K inhibitor PF-04691502 was purchased from Selleck Chemicals (Houston, TX, United States). TMZ was purchased from J&K Scientific (Beijing, China). All antibodies [AKT, phosphor-AKT (S473), phosphor-AKT (T308), mammalian target of rapamycin (mTOR), phosphor-mTOR (S2448), GSK3β, phosphor-GSK3β (S9), proline-rich AKT substrate of 40 kDa (PRAS40), phosphor-PRAS40 (T246), p70S6K, phosphor-p70S6K (T389), S6 ribosomal protein (S6RP), phosphor-S6RP (S240/244), smoothened (SMO), cyclin D1, and anti-cyclin–dependent kinase (CDK2)] were purchased from Cell Signaling Technology (Danvers, MA, United States). Anti-GLI1 and anti–β-actin antibodies were purchased from Abcam (Cambridge, United Kingdom) and Santa Cruz Biotechnology (Dallas, TX, United States), respectively.
Cell Viability Assay
Cell viability was assessed using the Cell Counting Kit-8 (CCK-8; Solarbio, Beijing, China). Briefly, 2,000 cells per well were seeded into a 96-well plate. After incubation overnight, the cells were treated with different concentrations (0.0512, 0.256, 1.28, 6.4, 32, 160, 800, 4,000, and 20,000 nM) of XH30 or PF-04691502 with three replicates for 72 h. Then, CCK-8 solution was added and incubated for 3 h, after which, absorbance values were measured at 450 nM using a microplate reader (BioTek Instruments, Inc., United States). The half-maximal inhibitory concentration (IC50) was calculated using GraphPad Prism v8.0.1 (La Jolla, CA, United States).
Colony Formation Assay
U251/TMZ and T98G cells were seeded at a density of 200 cells per well into six-well plates. After 24 h, the cells were treated with indicated concentrations (4, 20, 100, and 500 nM) of XH30 with three replicates. The culture medium with test compound was replaced every 3 days. After cell colonies formed, cells were washed with phosphate-buffered saline (PBS) and fixed with 4% paraformaldehyde for 15 min, then stained with crystal violet for 30 min, and washed with PBS. Finally, colonies were recorded with a photograph and measured by a microplate reader (BioTek).
Cell Cycle Analysis
Flow cytometry assays were used to analyze the cell cycle distribution as previously reported (Ji et al., 2018). In brief, U251/TMZ and T98G cells were dispensed into six-well plates at a density of 50,000 cells per well. After growing overnight in a humidified atmosphere of 5% CO2 at 37°C, cells were treated with indicated concentrations (4, 20, 100, and 500 nM) of XH30 for 24 h. Then, cells were harvested and fixed ice cold 70% ethanol overnight at –20°C, washed with PBS, and stained with propidium iodide (PI) solution containing PI (20 mg/ml) and RNase A (20 mg/ml) in PBS for 30 min. DNA contents were measured using the BD fluorescence-activated cell sorting verse flow cytometer (BD Biosciences, NJ, United States), and the cell cycle distribution was analyzed.
Immunoblotting Analysis
Cells or mice tumor tissues were collected and lysed in RIPA lysate buffer supplemented with 1% protease inhibitor cocktail and 1% phosphatase inhibitor cocktail (Solarbio, Beijing, China). Lysates were then centrifuged at 12,000 g for 30 min. Proteins were quantified using a bicinchoninic acid assay kit (Solarbio, Beijing, China). Resultant samples containing equal amounts of proteins were subjected to sodium dodecylsulfate–polyacrylamide gel electrophoresis and transferred to a polyvinylidene fluoride membrane (Millipore, Darmstadt, Germany). The membrane was blocked with TBST buffer containing 5% non-fat milk for 30 min and incubated with appropriate primary antibodies (1:1,000 dilution) in TBST at 4°C overnight. After washing with TBST, the membrane was incubated with horseradish peroxidase–conjugated secondary antibodies (1:2,000 dilution; Cell Signaling Technologies, Boston, MA) for 1 h at room temperature. Bound proteins were visualized using enhanced chemiluminescence and detected using ImageQuant LAS 4000 software.
Quantitative Real-Time Polymerase Chain Reaction Analysis
Total RNA from XH30-treated U251/TMZ or T98G cells was isolated using TRIzol reagent (Bioteke Corporation, China) according to the recommended procedures of the manufacturer. First-strand cDNA was synthesized from 1 μg of total RNA using the ReverTra Ace® qPCR RT Master Mix with gDNA Remover (Toyobo, Japan). Real-time polymerase chain reaction (PCR) was performed using the Analytikjena qTOWER detection system and the SYBR® Green RT-PCR master mix (Toyobo). Target sequences were amplified at 95°C for 1 min, followed by 40 cycles at 95°C for 15 s, 60°C for 15 s, and 72°C for 45 s. Fold changes in GLI1, paired box protein 6 (PAX6), and O6-methylguanine-DNA-MGMT (MGMT) gene expression were calculated according to the 2−ΔΔCt method. The primers sequences used to amplify specific regions of the indicated genes were as follows: GLI1 forward, ATGTTCAACTCGATGACCCCAC; GLI1 reverse, CAACTTGACTTCTGTCCCCACA; MGMT forward, ATGGAT GTTTGAGCGACACA; MGMT reverse, ATAGAGCAAGGGCAGCGTTA; PAX6 forward, AACGATAACATACCAAGCGTGT; PAX6 reverse, GGTCTGCCCGTTCAACATC.
Orthotopic Mouse Tumor Model and Subcutaneous Mouse Tumor Model
Eight- to 10-week-old female Balb/c athymic nude mice (SPF Biotechnology, Beijing, China) were housed in standard facilities. Human U251/TMZ cells in PBS were injected intracranially, 2.0 mm below the skull surface, according to a previously published protocol (Ji et al., 2018). Three days after surgery, mice were randomized to receive one of the four following treatments: a drug vehicle (delivered orally once per day for 9 days). TMZ (50 mg/kg, delivered orally once per day for 5 days), or XH30 (5 mg/kg, delivered orally once per day for 9 days). TMZ and XH30 were dissolved in 0.5% carboxymethylcellulose solution. Tumor volumes were monitored using an animal magnetic resonance imaging (MRI) scanner (PharmaScan 70/16 US, Bruker, Germany). The parameters for the MRI scans were as follows: a T2_TurboRARE, with TR/TE = 5,000/40, 6 averages, 20 × 20 field of view, and 0.5-mm slice thickness. The tumor volume on the basis of MRI was calculated as V = L × W × T, where L is the maximum length of tumor, W is the maximum width perpendicular to L, and T is the thickness of the tumor slice (set at 0.5 mm).
For the subcutaneous mice tumor model, female Balb/c athymic nude mice (eight to 10 weeks of age) were subcutaneously implanted with 1 × 107 U251/TMZ cells in 0.2 ml matrigel solution in the right flank. After 2 weeks, tumor issue was harvested sterilely, and tumor cells were extracted from the tissue homogenate. Then, the mice were implanted with 2 × 106 tumor cells each in the right flank. Seven days later, when the average tumor volumes reached to 100–300 mm3, the mice were randomized into four groups, in which either alone treatment or combination was administered, respectively, using the same dose regime as in the orthotopic model. Tumor volume and body weight were monitored twice a week. Tumor volume was calculated as V = 1/2 × L × W2, where L is the maximum length of tumor and W is the maximum width of tumor. The mice were euthanized at day 14, and tumor tissues were collected for immunoblotting.
All procedures were approved by the Ethics Committee for Animal Experiments of the Institute of Materia Medica, Chinese Academy of Medical Sciences and Peking Union Medical College and conducted under the Guidelines for Animal Experiments of Peking Union Medical College.
Statistical Analysis
Most statistical analyses were performed utilizing GraphPad Prism v8.0.1 (La Jolla, CA, United States), and significance levels were evaluated using analysis of variance (ANOVA) or T-tests, as appropriate. In our experiments, we distinguish between three of significance (***p < 0.001, **p < 0.01, and *p < 0.05, respectively).
On the basis of the cell viability assay, the combination index (CI) was calculated using the Chou–Talalay method, with CI = 1, CI < 1 and CI > 1 denoting an additive effect, synergism, and antagonism, respectively. CI = (D)1/(Dx)1 + (D)2/(Dx)2, where (Dx)1 and (Dx)2 represented concentrations of each drug alone to exert x% effect, while (D)1 and (D)2 were concentrations of drugs in combination to elicit the same effect.
Results
The PI3K Inhibitor XH30 Inhibited TMZ-Resistant GBM Cell Growth in vitro
Here, we assessed the anti-tumor activity of XH30 in two TMZ-resistant GBM cell lines: one with acquired TMZ resistance and overactivated hedgehog signaling (U251/TMZ) and one with natural TMZ resistance and elevated MGMT levels (T98G). As shown in Figure 1A, TMZ did not exhibit cytotoxicity at the highest concentration of 1,000 μM. XH30 suppressed cell proliferation of both U251/TMZ and T98G cells, with 72 h IC50 values of 191 and 183 nM, respectively. The cytotoxic effect was stronger at 72 h than at 24 or 48 h (Supplementary Figure S1). The positive control, PF-04691502, also inhibited the proliferation of both cell types. In a colony formation assay, the formation of cell colonies dose-dependently decreased after exposure to XH30 in both U251/TMZ and T98G cells (Figures 1B–E). These results indicated that XH30 exhibited strong inhibitory effects on the proliferation of TMZ-resistant GBM cells.
FIGURE 1. The PI3K inhibitor XH30 inhibited cell proliferation in temozolomide (TMZ)–resistant glioblastoma multiforme (GBM) cells. (A) The IC50 values of XH30, PF-04691502 (control), and TMZ in TMZ-resistant GBM cells over 72 h. Data are presented as means ± SD, n = 3. (B,D) Colony formation assay of XH30 in TMZ-resistant U251/TMZ and T98G cells, respectively. Representative images are shown for each group. (C,E) The inhibition ratio of colony formation assay of XH30 in U251/TMZ and T98G cells, respectively. Data are presented as means ± SD, n = 3.
XH30 Reduced Downstream Molecules in the PI3K Signaling Pathway in TMZ-Resistant GBM Cells
We assessed the inhibitory activity of XH30 on the PI3K signaling pathway in both U251/TMZ and T98G cells. XH30 dose-dependently blocked downstream molecules in the PI3K pathway including p-AKT, p-GSK3β, p-PRAS40, p-p70S6K, and p-S6RP (Figures 2A–D). At a concentration of 100 nM, XH30 strongly suppressed the phosphorylation of these signaling molecules. This inhibitory activity was more potent than the inhibitory effect of the positive control, PF-04691502, at the same concentrations.
FIGURE 2. The PI3K inhibitor XH30 suppressed the PI3K signaling pathway in TMZ-resistant GBM cells. (A) XH30 dose-dependently inhibited PI3K pathway signaling in TMZ-resistant U251/TMZ cells. The cells were incubated with XH30 at indicated concentrations (0.8, 4, 20, and 100 nmol/L) or PF-04691502 (100 nmol/L) as a control for 3 h. The experiment was repeated three times. (B) Relative expression levels of the major proteins in (A). Data are presented as means ± SD, n = 3. t-test, *p < 0.05, **p < 0.01, and ***p < 0.001 compared with control. (C) XH30 dose-dependently inhibited PI3K pathway signaling in TMZ-resistant T98G cells. The cells were incubated with XH30 at indicated concentrations (0.8, 4, 20, and 100 nmol/L) or PF-04691502 (100 nmol/L) as a control for 3 h. The experiment was repeated three times. (D) Relative expression levels of major proteins in (C). Data are presented as means ± SD, n = 3. t-test, *p < 0.05, **p < 0.01, and ***p < 0.001, compared with control.
XH30 Induced Cell Cycle Arrest in TMZ-Resistant GBM Cells
Next, we investigated whether XH30 could induce cell cycle arrest or apoptosis in both TMZ-resistant GBM cell lines. In U251/TMZ cells, both PF-04691502 and XH30 induced cell cycle arrest in the G0/G1 phases. The percentage of G1 phase increased from 48.62% to 76.30% after a concentration titration of XH30 (Figure 3A). In T98MG cells, cell cycle arrest was observed during the G1 phase after exposure to XH30. The percentage of G1 phase at a XH30 concentration of 500 nM was increased compared to the control group (77.78% vs. 50.07%) (Figure 3B). Moreover, XH30 downregulated the expression of cyclin D1 and CDK2, which are markers of the G1 phase (Figures 3C,D). However, apoptosis was not observed in either U251/TMZ or T98G cells after exposure to XH30 for 48 h. These data demonstrated that XH30 induces cell cycle arrest in TMZ-resistant GBM cells (Figures 3E,F).
FIGURE 3. XH30 induced G1 cell-cycle arrest in TMZ-resistant GBM cells. (A, B) XH30 induced cell cycle arrest in TMZ-resistant U251/TMZ and T98G cells in the G1 phase. Exposure to various concentrations of XH30 at indicated concentrations (4, 20, 100, and 500 nM) or the control, PF-04691502 (500 nM), for 24 h. The cells were stained with propidium iodide (PI) for flow cytometry analysis, n = 3. (C, D) XH30 downregulated markers of the G1 phase of the cell cycle. The protein levels of cyclin D1 and CDK2 were detected via immunoblotting in both U251/TMZ and T98G cells exposed to XH30 for 24 h. (E, F) The effect of XH30 on apoptosis in both U251/TMZ and T98G cells exposed to XH30 for 48 h.
XH30 Exhibited Antitumor Activity in TMZ-Resistant GBM in vivo
To further explore the antitumor activity of XH30 in TMZ-resistant GBM, we used a U251/TMZ orthotopic mouse model. Mice were given TMZ orally at doses of 50 mg/kg/day for 5 days or XH30 at either 5 mg/kg or 10 mg/kg daily for 9 days. In our MRI analysis, the images showed that XH30 suppressed tumor growth in brain (Figures 4A,B). At a dose of 10 mg/kg/day, tumor volume was significantly lower compared to the control group that received the drug vehicle (p < 0.01; Figure 4C; Supplementary Table S1). In this model, TMZ did not reduce tumor volume at a dose of 50 mg/kg/day with 6.2% tumor growth inhibition (TGI), whereas the TGI of TMZ at the same dose was 98.4% in a U251 (TMZ sensitive) orthotopic mice model (Supplementary Figures S3A, B). During the experiments, the body weight of mice in the XH30 group did not significantly decrease compared to that of the control group (Supplementary Figure S2A). Together, these data indicate that XH30 suppressed TMZ-resistant GBM growth in vivo.
FIGURE 4. XH30 repressed TMZ-resistant GBM growth in a mouse orthotopic xenograft model. (A) magnetic resonance imaging (MRI) T2-weighted image of intracranial tumors from the various groups from the U251/TMZ orthotopic model at day 9. Images in each lane show all tumor slices from one representative mouse in each group. The red curves indicate the tumors. (B) Representative MRI images from the U251/TMZ orthotopic model at day 9. The red arrows indicate the tumor. (C) Tumor volumes in the U251/TMZ orthotopic model at day 9. Data are presented as means ± SD, n = 5. ANOVA, **p < 0.01 compared with the vehicle control group.
XH30 Enhanced TMZ Cytotoxicity in TMZ-Resistant GBM
Because XH30 showed excellent antitumor activity against TMZ-resistant GBM and TMZ resistance is a current barrier to effective treatment in GBM, we investigated whether XH30 could sensitize GBM to TMZ. To do this, XH30 and TMZ were administered together in GBM cells. XH30 at concentrations of 20 and 100 nM enhanced TMZ cytotoxicity in both U251/TMZ and T98G cells (Figure 5A). Calculation of the CI gave values of <1 for combined XH30 and TMZ in both cell types, indicating a synergic effect of these drugs (Supplementary Table S2). Immunoblotting results showed that the combination of XH30 and TMZ increased the level of the DNA damage marker λH2AX to a greater extent compared with TMZ or XH30 alone, indicating that XH30 increases TMZ cytotoxicity (Figure 5B).
FIGURE 5. XH30 enhanced TMZ response in TMZ-resistant GBM. (A) The curve for combined XH30 and TMZ treatment in GBM cells for 72 h. Data are presented as means ± SD, n = 3. (B) Combination of XH30 with TMZ increased the level of λH2AX. Cells were incubated with XH30 (100 nM) or TMZ (250 μM) for 12 h. (C) Antitumor activity of XH30 in an U251/TMZ subcutaneous xenograft model. t-test, *p < 0.05 and ***p < 0.001, compared with the TMZ or XH30 group. Data are presented as means ± SD, n = 6. (D) Tumor weight in an U251/TMZ xenograft model. t-test, *p < 0.05 and ***p < 0.001. (E) Tumor issues in a U251/TMZ xenograft model. (F) Levels of proteins in the PI3K signaling pathway in tumor tissue from the U251/TMZ xenograft mice model.
We then evaluated whether combining XH30 and TMZ had synergistic antitumor effects in vivo. We employed a subcutaneous mouse model implanted with U251/TMZ cells. The groups that received XH30 or a combination of XH30 and TMZ both had significantly delayed tumor growth (Figure 5C). Treatment with either TMZ at 50 mg/kg/day or XH30 at 5 mg/kg/day suppressed tumor growth, with TGI values of 58.9% and 69.9%, respectively (Figures 5D,E). Combined TMZ and XH30 significantly suppressed tumor growth compared with the groups that received either TMZ or XH30 alone, with a TGI of 90.1%. Although TMZ exhibited antitumor activity, the effect was weaker than in the parent cell line U251 (TMZ sensitive) in the subcutaneous mouse model (58.9% vs. 94.9% TGI, Figure S3C and D). The body weight change with combined TMZ and XH30 treatment was within acceptable limits (Supplementary Figure S4). The immunoblotting data also showed that the phosphorylation proteins downstream to PI3K including mTOR, AKT, and S6RP decreased in tumor tissues with XH30 alone and in combination with TMZ (Figure 5F).
XH30 Repressed GLI1 via the Noncanonical Hedgehog Signaling Pathway to Increase TMZ Cytotoxicity
The hedgehog signaling pathway has been demonstrated to have a role in TMZ resistance (Li J. et al., 2016; Lee, 2016), and crosstalk between noncanonical hedgehog and PI3K signaling pathways has been documented (Riobo et al., 2006). Previously, we showed that the hedgehog pathway is overactive in both U251/TMZ and T98G cells (Ji et al., 2018). Therefore, we hypothesized that inhibition of PI3K may also inhibit the noncanonical hedgehog pathway to increase the response to TMZ. As shown in Figure 6A, the protein levels of GLI1, a key factor in the hedgehog signaling pathway, was dose-dependently decreased in U251/TMZ cells exposed to XH30 in various concentrations, along with decreased levels of phosphorylated AKT. SMO protein levels did not decrease after treatment with XH30. Similar results were also observed in T98G cells. GLI1 target genes, such as PAX6 and GLI1 itself, were downregulated in the presence of XH30 at 100 nM (Figure 6B), likely because XH30 suppresses noncanonical hedgehog signaling pathway via the blockade of PI3K. In T98G cells, another GLI1 target gene, MGMT, was also downregulated by XH30. Next, we investigated whether XH30 was able to block the hedgehog pathway in the presence of the hedgehog ligand SHH (sonic hedgehog). As shown in Figure 6C, SHH activated the hedgehog pathway in U251/TMZ and T98G, reflecting upregulated GLI1 expression, whereas XH30 reduced GLI1 expression. In the presence of SHH, XH30 partially reversed SHH-mediated GLI1 activation. The phosphorylation of downstream proteins in the PI3K pathway, including AKT and S6RP, was consistently decreased in cells exposed to XH30. Moreover, insulin-like growth factor 1 (IGF-1), which activates PI3K, also upregulated GLI1 expression, which can be attributed to the crosstalk between the PI3K signaling pathway and GLI1 (Figure 6D). In the presence of XH30, IGF-mediated elevation of GLI1 was partially attenuated. There was no obvious change in SMO protein levels. We also observed that XH30 and TMZ in combination maintained lower level of GLI1 expression in both mRNA and protein level in TMZ-resistant GBM cells (Figures 6E,F). Immunoblotting results of tumor tissues also showed that the protein levels of GLI1 decreased in XH30 and TMZ combination group compared to vehicle group (Figure 6G). Together, our results demonstrate that XH30 suppressed GLI1 via blockade of the noncanonical hedgehog signaling pathway.
FIGURE 6. XH30 suppressed GLI1 via blockade of the noncanonical hedgehog signal pathway. (A) XH30 dose-dependently reduced GLI1 protein levels in TMZ-resistant U251/TMZ and T98G cells. Cells were incubated with XH30 at indicated concentrations (0.8, 4, 20, and 100 nmol/L) for 12 h (B) XH30 downregulated the GLI1 target gene expression, as indicated by mRNA levels, in both U251/TMZ and T98G cells. Data are presented as means ± SD, n = 3. t-test, *p < 0.05, **p < 0.01, and ***p < 0.001, compared with control. (C) XH30 attenuated SHH-triggered GLI1 activation in U251/TMZ and T98G cells. The cells were incubated with XH30 at a concentration of 100 nM in the absence or presence of recombinant SHH (500 ng/ml) for 24 h. (D) XH30 reversed GLI1 activation by IGF-1 in U251/TMZ and T98G cells. Cells were incubated with XH30 at a concentration of 100 nM in the absence or presence of recombinant IGF-1 (50 ng/ml) for 12 h. (E) XH30 and TMZ in combination maintained lower mRNA levels of GLI1 in GBM cells. The cells were incubated with XH30 (100 nM) or TMZ (250 μM) for 12 h. t-test, **p < 0.01 and ***p < 0.001, compared with control. (F) XH30 and TMZ in combination maintained lower protein levels of GLI1 in GBM cells. The cells were incubated with XH30 (100 nM) or TMZ (250 μM) for 12 h. (G) Protein levels of GLI1 in tumor tissue from the U251/TMZ xenograft mice model.
Discussion
Glioblastoma is the most aggressive cancer of the brain in adults, and its prevalence is increasing, especially in China (Yang et al., 2013). Standard treatment includes surgery combined with radiotherapy and chemotherapy (Weller et al., 2014). So far, TMZ is the only chemotherapeutic option with confirmed efficacy and an acceptable safety profile in this cancer type (Rajaratnam et al., 2020). In addition, vascular endothelial growth factor inhibitors and tumor-treating fields have also been approved for the treatment of glioblastoma (Lassman et al., 2020; Seystahl et al., 2020; Zhang et al., 2020). However, patient outcomes from these treatments remain contentious. Clinical trials of immunotherapy with anti-programmed cell death 1 antibodies also did not meet clinical survival endpoints in patients with GBM (Khasraw et al., 2020; Maghrouni et al., 2021). Overactivation of the PI3K pathway occurs frequently in GBM, including in those with TMZ resistance (Sami and Karsy, 2013). In our previous study, XH30 exerted robust antitumor activity against TMZ-sensitive glioblastoma (Lin et al., 2018). Therefore, in this study, we tested this compound in GBM with either natural or acquired TMZ resistance. As expected, XH30 was shown to have acceptable antitumor activity against TMZ-resistant GBM both in vitro and in vivo, via inhibition of PI3K and downstream proteins and induction of cell cycle arrest. Similarly to other pan-PI3K inhibitors, we observed reductions in total white blood cells, neutrophils, and lymphocytes after XH30 treatment in our in vivo mouse model (Supplementary Figure S2B).
After standard treatment, GMB recurs in most people with this tumor type; TMZ resistance is a primary factor contributing to this process (Bocangel et al., 2002). For recurring tumors, repeated treatment with low doses of TMZ or treatment with lomustine (CCNU) has been recommended (Birzu et al., 2020; Di Nunno et al., 2020; Weller and Le Rhun, 2020). However, the clinical benefits of this approach are limited. The mechanism of TMZ resistance is extremely complex and includes overexpression of MGMT, aberrant activation of the hedgehog signaling pathway, overexpression of P-glycoprotein, and even metabolic reprogramming (Lee, 2016). Many approaches for overcoming TMZ resistance have been assessed (Happold and Weller, 2015), but none have been particularly successful. The hedgehog signal pathway represents an attractive target for glioblastoma treatment; inhibition of this pathway could overcome TMZ resistance (Shahi et al., 2008; Braun et al., 2012). In this context, a particularly interesting finding in our study is that XH30 dose-dependently decreased GLI-1 protein levels and downregulated its target genes including PAX6 and GLI1 itself in both U251/TMZ and T98G. In T98G cells, the mRNA level of MGMT, which may be regulated by GLI1, was reduced after XH30 treatment as well. This observation triggered us to explore the possible mechanism of XH30 as part of the hedgehog signaling pathway. There have been reports that PI3K signaling pathway is an important non-canonical activator of GLI1 and that targeting the PI3K/AKT pathway via GLI inhibition enhances drug sensitivity (Liang et al., 2017). Conversely, GLI1 reduces drug sensitivity via direct activation of the PI3K pathway in acute myeloid leukemia (Zhou et al., 2021). In our study, in the presence of SHH, XH30 suppressed the hedgehog pathway and partially reversed GLI1 activation. In addition, when we added IGF-1 to activate the PI3K signaling pathway, we found that GLI1 protein levels increased after IGF-1 stimulation. In the presence of XH30, this GLI1 level increase was partially attenuated. These results provided us with an additional clue that XH30 may play other roles in GBM treatment, which warrants further research.
Previous research has reported that inhibition of hedgehog signal pathway can enhance TMZ cytotoxicity and overcome TMZ resistance (Li J. et al., 2016; Ji et al., 2018). Therefore, we predicted that XH30 may also increase the response to TMZ via blockade of the non-canonical hedgehog pathway, allowing direct antitumor activity against TMZ-resistant GBM. In our experiments, XH30 enhanced the cytotoxicity of TMZ in both TMZ-resistant cell types in vitro. Treatment with combined XH30 and TMZ increased the level of λH2AX, a marker of DNA damage. Our in vivo studies in orthotopic mice also confirmed this synergistic effect. The combination of XH30 with TMZ yielded an improved antitumor activity compared with XH30 or TMZ treatment alone. This suggests that PI3K inhibitors could be tested as adjuvant treatment along with TMZ in patients with recurrent GMB.
In conclusion, the PI3K inhibitor XH30 exhibited robust antitumor activity in TMZ-resistant GBM; this compound is therefore a novel potential therapeutic option for TMZ-resistant GBM.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding authors.
Ethics Statement
The animal study was reviewed and approved by the Ethics Committee for Animal Experiments of the Institute of Materia Medica, Chinese Academy of Medical Sciences and Peking Union Medical College.
Author Contributions
MJ designed the experiments. MJ, ZZ, and CW conducted all the experiments. SL synthesized the compound. MJ drafted the manuscript. NX, HX, and XC reviewed and edited the manuscript and supervised the entire study. All authors reviewed and approved the final version of the manuscript.
Funding
This work was supported by the Non-Profit Central Research Institute Fund of Chinese Academy of Medical Sciences (Grant 2018PT35003).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank Charlesworth (https://www.cwauthors.com/) for its linguistic assistance during the preparation of this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.749242/full#supplementary-material
References
Avery, J. T., Zhang, R., and Boohaker, R. J. (2021). GLI1: A Therapeutic Target for Cancer. Front. Oncol. 11, 673154. doi:10.3389/fonc.2021.673154
Birzu, C., French, P., Caccese, M., Cerretti, G., Idbaih, A., Zagonel, V., et al. (2020). Recurrent Glioblastoma: From Molecular Landscape to New Treatment Perspectives. Cancers (Basel) 13. doi:10.3390/cancers13010047
Bocangel, D. B., Finkelstein, S., Schold, S. C., Bhakat, K. K., Mitra, S., and Kokkinakis, D. M. (2002). Multifaceted Resistance of Gliomas to Temozolomide. Clin. Cancer Res. 8, 2725–2734.
Braun, S., Oppermann, H., Mueller, A., Renner, C., Hovhannisyan, A., Baran-Schmidt, R., et al. (2012). Hedgehog Signaling in Glioblastoma Multiforme. Cancer Biol. Ther. 13, 487–495. doi:10.4161/cbt.19591
Castresana, J. S., and Melendez, B. (2021). Molecular and Cellular Mechanisms of Glioblastoma. Cells 10. doi:10.3390/cells10061456
Chien, C. H., Hsueh, W. T., Chuang, J. Y., and Chang, K. Y. (2021). Dissecting the Mechanism of Temozolomide Resistance and its Association with the Regulatory Roles of Intracellular Reactive Oxygen Species in Glioblastoma. J. Biomed. Sci. 28, 18. doi:10.1186/s12929-021-00717-7
Colardo, M., Segatto, M., and Di Bartolomeo, S. (2021). Targeting RTK-Pi3k-mTOR Axis in Gliomas: An Update. Int. J. Mol. Sci. 22. doi:10.3390/ijms22094899
Deangelis, L. M. (2001). Brain Tumors. N. Engl. J. Med. 344, 114–123. doi:10.1056/NEJM200101113440207
Di Nunno, V., Franceschi, E., Tosoni, A., Di Battista, M., Gatto, L., Lamperini, C., et al. (2020). Treatment of Recurrent Glioblastoma: State-Of-The-Art and Future Perspectives. Expert Rev. Anticancer Ther. 20, 785–795. doi:10.1080/14737140.2020.1807949
Dymova, M. A., Kuligina, E. V., and Richter, V. A. (2021). Molecular Mechanisms of Drug Resistance in Glioblastoma. Int. J. Mol. Sci. 22. doi:10.3390/ijms22126385
Gil Del Alcazar, C. R., Hardebeck, M. C., Mukherjee, B., Tomimatsu, N., Gao, X., Yan, J., et al. (2014). Inhibition of DNA Double-Strand Break Repair by the Dual PI3K/mTOR Inhibitor NVP-Bez235 as a Strategy for Radiosensitization of Glioblastoma. Clin. Cancer Res. 20, 1235–1248. doi:10.1158/1078-0432.CCR-13-1607
Hainsworth, J. D., Becker, K. P., Mekhail, T., Chowdhary, S. A., Eakle, J. F., Wright, D., et al. (2019). Phase I/II Study of Bevacizumab with BKM120, an Oral PI3K Inhibitor, in Patients with Refractory Solid Tumors (Phase I) and Relapsed/refractory Glioblastoma (Phase II). J. Neurooncol. 144, 303–311. doi:10.1007/s11060-019-03227-7
Happold, C., and Weller, M. (2015). New Insights into Acquired Temozolomide Resistance in Glioblastoma? Brain 138, 3468–3470. doi:10.1093/brain/awv301
Ji, M., Wang, L., Chen, J., Xue, N., Wang, C., Lai, F., et al. (2018). CAT3, a Prodrug of 13A(S)-3-hydroxyl-6,7-dimethoxyphenanthro[9,10-b]-indolizidine, Circumvents Temozolomide-Resistant Glioblastoma via the Hedgehog Signaling Pathway, Independently of O6-Methylguanine DNA Methyltransferase Expression. Onco. Targets Ther. 11, 3671–3684. doi:10.2147/OTT.S163535
Ji, M., Wang, D., Lin, S., Wang, C., Li, L., Zhang, Z., et al. (2021). A Novel PI3K Inhibitor XH30 Suppresses Orthotopic Glioblastoma and Brain Metastasis in Mice Models. Acta Pharm. Sin. B.
Khasraw, M., Reardon, D. A., Weller, M., and Sampson, J. H. (2020). PD-1 Inhibitors: Do They Have a Future in the Treatment of Glioblastoma? Clin. Cancer Res. 26, 5287–5296. doi:10.1158/1078-0432.CCR-20-1135
Langhans, J., Schneele, L., Trenkler, N., Von Bandemer, H., Nonnenmacher, L., Karpel-Massler, G., et al. (2017). The Effects of PI3K-Mediated Signalling on Glioblastoma Cell Behaviour. Oncogenesis 6, 398. doi:10.1038/s41389-017-0004-8
Lassman, A. B., Joanta-Gomez, A. E., Pan, P. C., and Wick, W. (2020). Current Usage of Tumor Treating fields for Glioblastoma. Neurooncol. Adv. 2, vdaa069. doi:10.1093/noajnl/vdaa069
Lee, S. Y. (2016). Temozolomide Resistance in Glioblastoma Multiforme. Genes Dis. 3, 198–210. doi:10.1016/j.gendis.2016.04.007
Li, J., Cai, J., Zhao, S., Yao, K., Sun, Y., Li, Y., et al. (2016a). GANT61, a GLI Inhibitor, Sensitizes Glioma Cells to the Temozolomide Treatment. J. Exp. Clin. Cancer Res. 35, 184. doi:10.1186/s13046-016-0463-3
Li, X., Wu, C., Chen, N., Gu, H., Yen, A., Cao, L., et al. (2016b). PI3K/Akt/mTOR Signaling Pathway and Targeted Therapy for Glioblastoma. Oncotarget 7, 33440–33450. doi:10.18632/oncotarget.7961
Liang, H., Zheng, Q. L., Fang, P., Zhang, J., Zhang, T., Liu, W., et al. (2017). Targeting the PI3K/AKT Pathway via GLI1 Inhibition Enhanced the Drug Sensitivity of Acute Myeloid Leukemia Cells. Sci. Rep. 7, 40361. doi:10.1038/srep40361
Liau, L. M. (2021). Glioblastoma: Molecular Mechanisms and Clinical Trials. Neurosurg. Clin. N. Am. 32, xiii. doi:10.1016/j.nec.2021.03.001
Lin, S., Wang, C., Ji, M., Wu, D., Lv, Y., Zhang, K., et al. (2018). Discovery and Optimization of 2-Amino-4-Methylquinazoline Derivatives as Highly Potent Phosphatidylinositol 3-Kinase Inhibitors for Cancer Treatment. J. Med. Chem. 61, 6087–6109. doi:10.1021/acs.jmedchem.8b00416
Maghrouni, A., Givari, M., Jalili-Nik, M., Mollazadeh, H., Bibak, B., Sadeghi, M. M., et al. (2021). Targeting the PD-1/pd-L1 Pathway in Glioblastoma Multiforme: Preclinical Evidence and Clinical Interventions. Int. Immunopharmacol. 93, 107403. doi:10.1016/j.intimp.2021.107403
Melamed, J. R., Morgan, J. T., Ioele, S. A., Gleghorn, J. P., Sims-Mourtada, J., and Day, E. S. (2018). Investigating the Role of Hedgehog/GLI1 Signaling in Glioblastoma Cell Response to Temozolomide. Oncotarget 9, 27000–27015. doi:10.18632/oncotarget.25467
Mizuarai, S., Kawagishi, A., and Kotani, H. (2009). Inhibition of p70S6K2 Down-Regulates Hedgehog/GLI Pathway in Non-small Cell Lung Cancer Cell Lines. Mol. Cancer 8, 44. doi:10.1186/1476-4598-8-44
Perazzoli, G., Prados, J., Ortiz, R., Caba, O., Cabeza, L., Berdasco, M., et al. (2015). Temozolomide Resistance in Glioblastoma Cell Lines: Implication of MGMT, MMR, P-Glycoprotein and CD133 Expression. PLoS One 10, e0140131. doi:10.1371/journal.pone.0140131
Radoul, M., Chaumeil, M. M., Eriksson, P., Wang, A. S., Phillips, J. J., and Ronen, S. M. (2016). MR Studies of Glioblastoma Models Treated with Dual PI3K/mTOR Inhibitor and Temozolomide:Metabolic Changes Are Associated with Enhanced Survival. Mol. Cancer Ther. 15, 1113–1122. doi:10.1158/1535-7163.MCT-15-0769
Rajaratnam, V., Islam, M. M., Yang, M., Slaby, R., Ramirez, H. M., and Mirza, S. P. (2020). Glioblastoma: Pathogenesis and Current Status of Chemotherapy and Other Novel Treatments. Cancers (Basel) 12. doi:10.3390/cancers12040937
Ranjan, A., and Srivastava, S. K. (2017). Penfluridol Suppresses Glioblastoma Tumor Growth by Akt-Mediated Inhibition of GLI1. Oncotarget 8, 32960–32976. doi:10.18632/oncotarget.16515
Riobó, N. A., Lu, K., Ai, X., Haines, G. M., and Emerson, C. P. (2006). Phosphoinositide 3-kinase and Akt Are Essential for Sonic Hedgehog Signaling. Proc. Natl. Acad. Sci. U S A. 103, 4505–4510. doi:10.1073/pnas.0504337103
Sami, A., and Karsy, M. (2013). Targeting the PI3K/AKT/mTOR Signaling Pathway in Glioblastoma: Novel Therapeutic Agents and Advances in Understanding. Tumour Biol. 34, 1991–2002. doi:10.1007/s13277-013-0800-5
Seystahl, K., Hentschel, B., Loew, S., Gramatzki, D., Felsberg, J., Herrlinger, U., et al. (2020). Bevacizumab versus Alkylating Chemotherapy in Recurrent Glioblastoma. J. Cancer Res. Clin. Oncol. 146, 659–670. doi:10.1007/s00432-019-03086-9
Shahi, M. H., Lorente, A., and Castresana, J. S. (2008). Hedgehog Signalling in Medulloblastoma, Glioblastoma and Neuroblastoma. Oncol. Rep. 19, 681–688. doi:10.3892/or.19.3.681
Singh, R., Dhanyamraju, P. K., and Lauth, M. (2017). DYRK1B Blocks Canonical and Promotes Non-canonical Hedgehog Signaling through Activation of the mTOR/AKT Pathway. Oncotarget 8, 833–845. doi:10.18632/oncotarget.13662
Touat, M., Idbaih, A., Sanson, M., and Ligon, K. L. (2017). Glioblastoma Targeted Therapy: Updated Approaches from Recent Biological Insights. Ann. Oncol. 28, 1457–1472. doi:10.1093/annonc/mdx106
Wang, Y., Ding, Q., Yen, C. J., Xia, W., Izzo, J. G., Lang, J. Y., et al. (2012). The Crosstalk of mTOR/S6K1 and Hedgehog Pathways. Cancer Cell 21, 374–387. doi:10.1016/j.ccr.2011.12.028
Weller, M., and Le Rhun, E. (2020). How Did Lomustine Become Standard of Care in Recurrent Glioblastoma? Cancer Treat. Rev. 87, 102029. doi:10.1016/j.ctrv.2020.102029
Weller, M., Van Den Bent, M., Hopkins, K., Tonn, J. C., Stupp, R., Falini, A., et al. (2014). European Association for Neuro-Oncology Task Force on Malignant, GEANO Guideline for the Diagnosis and Treatment of Anaplastic Gliomas and Glioblastoma. Lancet Oncol. 15, e395–e403. doi:10.1016/s1470-2045(14)70011-7
Wen, P. Y., Touat, M., Alexander, B. M., Mellinghoff, I. K., Ramkissoon, S., Mccluskey, C. S., et al. (2019). Buparlisib in Patients with Recurrent Glioblastoma Harboring Phosphatidylinositol 3-Kinase Pathway Activation: An Open-Label, Multicenter, Multi-Arm, Phase II Trial. J. Clin. Oncol. 37, 741–750. doi:10.1200/JCO.18.01207
Yang, P., Wang, Y., Peng, X., You, G., Zhang, W., Yan, W., et al. (2013). Management and Survival Rates in Patients with Glioma in China (2004-2010): a Retrospective Study from a Single-Institution. J. Neurooncol. 113, 259–266. doi:10.1007/s11060-013-1103-9
Yun, E. J., Kim, S., Hsieh, J. T., and Baek, S. T. (2020). Wnt/β-catenin Signaling Pathway Induces Autophagy-Mediated Temozolomide-Resistance in Human Glioblastoma. Cell Death Dis 11, 771. doi:10.1038/s41419-020-02988-8
Zając, A., Sumorek-Wiadro, J., Langner, E., Wertel, I., Maciejczyk, A., Pawlikowska-Pawlęga, B., et al. (2021). Involvement of PI3K Pathway in Glioma Cell Resistance to Temozolomide Treatment. Int. J. Mol. Sci. 22. doi:10.3390/ijms22105155
Zhang, C., Du, J., Xu, W., Huang, H., and Gao, L. (2020). The Value of Tumor Treating Fields in Glioblastoma. J. Korean Neurosurg. Soc. 63, 681–688. doi:10.3340/jkns.2019.0224
Zhang, P., Chen, X. B., Ding, B. Q., Liu, H. L., and He, T. (2018). Down-regulation of ABCE1 Inhibits Temozolomide Resistance in Glioma through the PI3K/Akt/NF-Κb Signaling Pathway. Biosci. Rep. 38. doi:10.1042/BSR20181711
Zhao, H. F., Wang, J., Shao, W., Wu, C. P., Chen, Z. P., To, S. T., et al. (2017). Recent Advances in the Use of PI3K Inhibitors for Glioblastoma Multiforme: Current Preclinical and Clinical Development. Mol. Cancer 16, 100. doi:10.1186/s12943-017-0670-3
Zheng, Y., Liu, L., Wang, Y., Xiao, S., Mai, R., Zhu, Z., et al. (2021). Glioblastoma Stem Cell (GSC)-derived PD-L1-Containing Exosomes Activates AMPK/ULK1 Pathway Mediated Autophagy to Increase Temozolomide-Resistance in Glioblastoma. Cell. Biosci. 11, 63. doi:10.1186/s13578-021-00575-8
Zhou, C., Du, J., Zhao, L., Liu, W., Zhao, T., Liang, H., et al. (2021). GLI1 Reduces Drug Sensitivity by Regulating Cell Cycle through PI3K/AKT/GSK3/CDK Pathway in Acute Myeloid Leukemia. Cell. Death Dis. 12, 231. doi:10.1038/s41419-021-03504-2
Keywords: glioblastoma, TMZ, PI3K, hedgehog, GLI1
Citation: Ji M, Zhang Z, Lin S, Wang C, Jin J, Xue N, Xu H and Chen X (2021) The PI3K Inhibitor XH30 Enhances Response to Temozolomide in Drug-Resistant Glioblastoma via the Noncanonical Hedgehog Signaling Pathway. Front. Pharmacol. 12:749242. doi: 10.3389/fphar.2021.749242
Received: 29 July 2021; Accepted: 03 November 2021;
Published: 26 November 2021.
Edited by:
Qingbin Cui, University of Toledo, United StatesReviewed by:
Roberto Würth, German Cancer Research Center (DKFZ), GermanyMounir Tilaoui, Waterford Institute of Technology, Ireland
Copyright © 2021 Ji, Zhang, Lin, Wang, Jin, Xue, Xu and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nina Xue, angelnina@imm.ac.cn; Heng Xu, xuheng@imm.ac.cn; Xiaoguang Chen, chxg@imm.ac.cn
†These authors have contributed equally to this work