Next Article in Journal
Oestrogenic Regulation of Mitochondrial Dynamics
Previous Article in Journal
Air Plasma-Activated Medium Evokes a Death-Associated Perinuclear Mitochondrial Clustering
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 3
10.3390/ijms23031123
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

EN1 Regulates Cell Growth and Proliferation in Human Glioma Cells via Hedgehog Signaling

by
Jinchun Chang
1,2,
Chenjia Guo
3,
Jianyu Li
3,
Zhangqian Liang
3,
Yankai Wang
2,
Anliang Yu
2,
Runze Liu
2,
Yuting Guo
2,
Jian Chen
3,4,* and
Song Huang
2,*
1
College of Life Sciences, Beijing Normal University, Beijing 100875, China
2
National Institute of Biological Sciences, Beijing 102206, China
3
Chinese Institute for Brain Research, Beijing 102206, China
4
Institute of Functional Nano and Soft Materials (FUNSOM) & Collaborative Innovation Center of Suzhou Nano Science and Technology, Soochow University, Suzhou 215123, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(3), 1123; https://doi.org/10.3390/ijms23031123
Submission received: 10 December 2021 / Revised: 9 January 2022 / Accepted: 17 January 2022 / Published: 20 January 2022
(This article belongs to the Section Molecular Neurobiology)
Browse Figures
Versions Notes

Abstract

:
Glioblastoma is an aggressive cancer of the nervous system that accounts for the majority of brain cancer-related deaths. Through cross-species transcriptome studies, we found that Engrailed 1 (EN1) is highly expressed in serum-free cultured glioma cells as well as glioma tissues, and increased expression level predicts a worse prognosis. EN1 controls glioma cell proliferation, colony formation, migration, and tumorigenic capacity in vivo. It also influences sensitivity of glioma cells to γ-ray irradiation by regulating intracellular ROS levels. Mechanistically, EN1 influences Hedgehog signaling by regulating the level of Gli1 as well as primary cilia length and the primary cilia transport-related protein TULP3. In conclusion, we demonstrate that EN1 acts as an oncogenic regulator that contributes to glioblastoma pathogenesis and could serve as a diagnostic/prognostic marker and therapeutic target for glioblastoma.

1. Introduction

Gliomas are the most common types of brain cancer, and glioblastoma is the most aggressive form of cancer that begins within the brain [1]. It accounts for more than 60% of all brain tumors in adults, and has been designated as a Grade IV tumor of the central nervous system by the World Health Organization [2]. The median survival of glioblastoma patients is generally 14–16 months, with fewer than 3–5% of patients surviving >5 years [3]. Currently, the standard treatment for glioblastoma is maximal safe resection followed by radiotherapy and chemotherapy with temozolomide [4]. In spite of intense efforts to optimize the treatment of gliomas over the past decades, no significant advancement in patient survival has been achieved. Notably, glioblastoma is also prone to chemotherapy and radiation resistance [5].
The gene Engrailed was first discovered by Ecker in a spontaneous mutant of Drosophila melanogaster in 1929, and it has been identified as a conserved gene cassette throughout holometabolous insects [6]. Vertebrates have two Engrailed genes (EN1 and EN2), and these homeobox-containing transcription factors function in regulating early embryonic development [7]. During the development of mouse, EN1 and EN2 have overlapping functions in regulating formation of additional fissures and for extensive cerebellar growth [8,9], but mutations in EN1 and EN2 produce different phenotypes [10]. In humans, aberrant expression of these genes has been associated with the pathogenesis of many types of tumors, including breast cancer [11,12], prostate cancer [13,14], colorectal cancer [15], and ovarian cancer [16,17]. Elevated EN2 expression in particular has been associated with poor survival, and this gene has been proposed as a biomarker for breast cancer screening [18]. The molecular basis of Engrailed-mediated pathogenic impacts in glioma remains unknown.
The Hedgehog signaling pathway was first identified in the Drosophila melanogaster; it is an evolutionarily conserved pathway of signal transmission from the cell membrane to the nucleus [19,20]. Its plays a significant role in the normal embryonic development and organogenesis of almost all organs in mammals [21]. Aberrant activation of Hedgehog signaling is one of the key oncogenic events in pediatric brain tumor medulloblastomas [22] and is also involved in glioma pathogenesis and glioma stem-cell maintenance [23].
Here, seeking to better understand the aggressive nature of glioblastomas, we investigated glioma stem cells obtained from isolated stem-cell cultures of mouse spontaneous glioma tumor samples [24]. Comparative transcriptomics of mouse glioma and normal stem cells, glioblastoma tumors, and normal tissues revealed that the mouse ortholog of Engrailed is highly expressed in gliomas. Our follow-up analyses ultimately indicate that EN1 is tumorigenic for gliomas, regulating the proliferation and growth of cancerous cells by modulating Hedgehog signaling activity.

2. Results

2.1. EN1 Is Upregulated in Glioblastoma and Is Related to Poor Patient Prognosis

Glioma stem cells (GSCs) are believed to contribute to tumor maintenance, therapeutic resistance, and tumor recurrence [25,26]. In order to identify key molecules regulating GSC, we examined glioma stem cells obtained from mouse spontaneous glioma tumor samples [24]. Our transcriptomics analysis of glioma stem cells and normal neural stem cells (NSCs) indicated that Engrailed was the top-ranking differentially expressed gene between GSCs and NSCs (Figure 1A). We also found that the Engrailed gene was expressed at higher levels in mouse glioblastoma tumors than in normal tissues (Figure 1A), and it was mostly expressed by tumor cells (Figure S7).
To examine the potential clinical relevance of Engrailed expression in gliomas, we analyzed clinical sample data from TCGA database and found that Engrailed 1 expression was significantly increased in glioblastoma tissue samples as compared to normal tissue specimens and to other low-grade gliomas (p < 0.001, Figure 1B). Similar expression trends were detected in the datasets studied by Henry Ford Hospital [27] (Figure 1C), further indicating that increased EN1 expression is associated with higher glioma grade.
To further assess the correlation between expression level of EN1 and patient survival with glioblastoma, we performed Kaplan–Meier survival analysis on the basis of the EN1 expression levels of glioblastoma patients from TCGA databases and the Chinese genome project database [28]. The overall survival of glioblastoma patients with high levels of EN1 expression was significantly poorer than that of patients with lower EN1 expression (p < 0.01) (Figure 1D and Figure S1). These findings indicated that EN1 may be associated with progression of glioblastoma.

2.2. EN1 Knockdown Restrains the Proliferation and Migration of Glioblastoma Cells In Vitro and In Vivo

To evaluate the functions of EN1 in glioma cells, we first assessed EN1 expression in four glioma cell lines, and we found that U-118 MG, U251, and A172 had higher EN1 expression than did the U-87 MG cell line (Figure 1E). We next used EN1-shRNA to knockdown (KD) the expression of EN1 in the three glioblastoma cell lines with relatively high EN1 expression, which led to a clear inhibition of proliferation (Figure 2A–C and Figure S2A,B). These findings were later supported by our observation of inhibited proliferation of EN1 knockout lines by CRISPR/Cas9 technology (Table S1 and Figure S3). We also observed that deficiency of EN1 reduced the colony formation abilities of U251 cells (Figure 2D). Conversely, we found that overexpression of EN1 did not increase cellular proliferation, but enhanced the clonogenicity of U-87 MG glioma cells (Figure 2H).
We performed wound-healing assays to evaluate the effect of EN1 on cell migration and discovered that EN1 KD significantly repressed glioblastoma cell migration (Figure 2E,F), whereas increased migration was observed for glioblastoma cells with EN1 overexpression (Figure S2C).
Finally, to study the roles of EN1 in tumor formation in vivo, we transplanted EN1 KD and control U251 cells intracranially into nude mice and found that reducing the EN1 level significantly prolonged the survival of tumor-bearing mice (Figure 2G). Collectively, these results indicate that EN1 functions directly in glioblastoma cell growth in vitro and in vivo.

2.3. EN1 Knockdown Elevated Intracellular ROS Levels and Susceptibility to γ-ray Irradiation

A basal level of reactive oxygen species (ROS) is necessary for cellular proliferation and differentiation [29]; to determine the cellular ROS change after EN1 knockdown, fluorescence sensitive probe H2DCF-DA was used for ROS characterization by flow cytometry measurement. The analyses demonstrated that EN1 knockdown caused ROS increase in glioma cells compared with scramble controls (Figure 3A). In contrast, overexpression of EN1 reduced ROS generation (27.15% vs. 51.6% in the control group, p < 0.05), which may allow for cell survival during various stressors. These data demonstrated that EN1 is an important regulator of glioma cell ROS generation.
Radiation therapy usually follows surgery in treatment of glioma, especially for high-grade gliomas [30,31]. Ionizing radiation produces ROS, which induce metabolic oxidative stress and prolonged cell injury [32], and the increased ROS level also inhibits the survival and self-renewal of GSCs [33]. Since EN1 knockdown caused ROS increase, we assessed the impact of modulating EN1 on responsivity to radiotherapy using colony formation assays. U-118 MG, U251, and A172 cells were treated with graded doses of radiation (0, 0.5, 1, 1.5, and 2 Gy). Colony-forming efficiency was assessed 2 weeks later, and the surviving fraction was calculated. The colony number of EN1 KD cell lines after radiotherapy was significantly lower than the scramble control cell lines (Figure 3B), suggesting that EN1 might modulate the cellular radiation sensitivity of glioma cells.

2.4. EN1 Increases Hedgehog Signaling Activity in Glioma Cells

To investigate the potential mechanisms of EN1 in the glioma cell line, we used RNA sequencing to assess glioma cells after knocking down EN1. RNA-seq analysis of these EN1 KD cell lines revealed that EN1 disruption induces differential expression of genes in multiple signaling pathways, including the smoothened signaling pathway (Hedgehog signaling pathway), glutathione metabolic process, protein localization to cilium, regulation of cellular senescence, and DNA replication (Figure 4A).
As the Hedgehog signaling pathway is a major regulator of many fundamental processes including stem-cell maintenance and oncogenesis [34], we next examined the expression of given glioma-associated oncogene homolog 1 (Gli1), which is known to function as a nuclear mediator of the Hedgehog pathway [35]. Gli1 expression was found to be significantly decreased in the EN1 KD cell lines as compared to scramble control cells (Figure 4B and Figure S4). We next carried out functional recovery experiments to determine whether EN1 exerts its functions in glioma cells via the Hedgehog pathway. We evaluated the colony-forming capacity of U251 and U251 overexpressing Gli1 (U251-OE-Gli1) with EN1 knockout mediated by CRISPR/Cas9. Briefly, we found that the restrained colony-forming abilities caused by the knockout of EN1 was reversed in the Gli1-OE cell background (Figure 4C,D). These data suggested that EN1 promotes glioma proliferation by regulating the Hedgehog pathway.

2.5. EN1 Knockdown Upregulates the Expression of TULP3

RNA-seq results also revealed that alteration of TULP3 was common to all examined EN1 KD glioma cell lines (Figure 5A and Figure S5A,B). The TULP3 gene product is an adaptor protein that traffics membrane proteins to primary cilia [36], which is the central organelle for the transduction of the Hedgehog signaling pathway [37]. A previous study showed that TULP3 is necessary for the trafficking of Gpr161 to cilia, while Gpr161 functions as a negative regulator of the Sonic hedgehog (SHH) pathway [38]. Furthermore, in mouse embryo studies, TULP3 was similarly described as a negative regulator of Hedgehog signaling [39]. Our qPCR and immunoblotting verified that TULP3 level was significantly increased in EN1 KD cell lines (Figure 5B,C). In addition, we also detected a significant decrease in the number of primary cilia in EN1 KD cells (Figure 5D,E and Figure S5C).

3. Discussion

Through cross-species transcriptome studies, we found that EN1 is highly expressed in gliomas, that high expression of EN1 in gliomas affects cell proliferation, colony formation, and resistance to radiotherapy, as well as promoting tumor growth in vivo, and that EN1 may function as an oncogene in gliomas. Therefore, EN1 expression alone can be a strong predictor of patient prognosis; even in patients with glioblastoma, where survival time is often very short, patients with high EN1 expression have significantly worse survival outcomes (Figure 1D).
EN1 has been reported as a key regulator that contributes to the progression of a variety of human cancers, including breast cancer [11,12] and salivary gland adenoid cystic carcinoma [40]. Previous studies found that EN1 could negatively regulate Wnt/β-catenin transcriptional activity in human colon cancer cells [41]. In addition, cytoplasmic or secreted EN1 has been previously reported as a tumor progression marker, while most of its function in tumor development is exercised by repressing or activating the expression of genes in the nucleus. Here, through RNA-seq analysis and further validation, we found that part of the function of EN1 in glioma cells was achieved through regulating SHH pathway activity.
The SHH pathway is also activated in multiple types of tumors, where it affects cancer cell proliferation, growth, and survival [42,43]. Both EN1 and Hedgehog signaling play crucial roles in developmental processes. EN1 is essential in the early establishment and maintenance of the midbrain/hindbrain [44]; interestingly, SHH signaling together with Wnt signaling also plays essential roles in the patterning and development of the midbrain/hindbrain [7]. In Drosophila, interactions between Engrailed (en) and Hedgehog signaling pathways revealed that the Hedgehog pathway is transcriptionally activated by the Engrailed protein [45]. Our current study adds further layers to understanding of the regulation of SHH by EN1, through modulating the number of primary cilia and protein trafficking in primary cilia by the TULP3 protein.
In conclusion, our study indicates that EN1 is frequently elevated in glioma tissues and associated with poor prognosis, and it may, therefore, be a promising diagnostic and prognostic marker for glioma patients. Mechanistically, we show that EN1 promotes glioma cell proliferation and migration, doing so at least partially through regulating Hedgehog pathway activity. Taken together, these findings suggest that EN1 is an oncogene in glioma that can be a promising target for the development of innovative therapies against glioma, a deadly disease largely lacking effective treatments.

4. Materials and Methods

4.1. Cell Lines and Culture

The glioblastoma cell lines U-118 MG, U251, A172, U-87 MG, and T98G were purchased from the Procell Life Science & Technology Co., Ltd. (Wuhan, China). The HEK 293T cell lines were provided by Dr. Ting Han of NIBS (Beijing, China). All cell lines were maintained in Dulbecco’s modified Eagle’s medium (DMEM) (Gibco, Waltham, MA, USA) and were supplemented with 10% fetal bovine serum (FBS) (Gibco, USA) and 1% penicillin/streptomycin antibiotic mixture (Gibco, Waltham, MA, USA). All cells were maintained at 37 °C with 5% CO2.
The glioma stem cells were obtained from isolated stem cell cultures of mouse spontaneous glioma tumor samples. The genotypes of the spontaneously tumor-bearing mice were hGFAP-Cre; p53fl/fl; ptenfl/+; nf1fl/+ [24]. The tumor samples were cut into small pieces, digested and blown, and cultured in vitro using glioma stem-cell medium. The cell line obtained after adherent growth was the glioma stem-cell line.

4.2. Plasmids and Lentiviral Transduction

The PLKO.1 vector was used to clone the shRNAs targeting EN1 to obtain EN1 knockdown cell lines, and the EN1 gene-deficiency monoclonal cell line was obtained by CRISPR/Cas9-mediated genome editing. For overexpression of EN1 or Gli1 in glioma cells, the full-length EN1 or Gli1 cDNA was amplified and then inserted into the FCTP vector (Table S2 and Figure S6). EN1 loss- and gain-of-function stable cell lines were generated via retroviral infection. Briefly, the HEK 293T cells were transfected with the target plasmids, pSPAX2 and pMD2G, at a ratio of 5:3:2 with Lipofectamine™ 3000 Transfection Reagent (Thermo Scientific, Waltham, MA, USA); virus particles were collected 48 h after transfection. The U-118 MG, U251, U-87 MG, A172, and T98G cells were infected with recombinant lentivirus transducing units using 1 μg/mL polybrene (Sigma-Aldrich, Burlington, MA, USA), and then 1 or 2 μg/mL puromycin (Thermo Scientific, Waltham, MA, USA) was used to screen positive cells for 3 days to acquire stable cell lines.

4.3. Quantitative PCR (qPCR)

Total RNA from shRNA-treated and EN1-overexpression cells was isolated from cells using a HP Total RNA Kit (Omega, Boca Raton, FL, USA). cDNAs were synthesized using KOD SYBR qPCR Mix (TOYOBO, Osaka, Japan) according to the manufacturer’s instructions. The 2−ΔΔCq method was used to calculate the expression level; GAPDH was used as a reference gene for quantification assays.

4.4. Immunoblotting

Cells were lysed with RIPA buffer, and protein concentrations were determined using the BCA Protein Assay Reagent (Thermo Scientific, Waltham, MA, USA). Equal amounts of protein were subjected to 10% SDS-PAGE before being transferred to a nitrocellulose membrane. The membrane was incubated (overnight, 4 °C) with the following primary antibodies: (EN1-2377, NIBS, Beijing, China), (EN1-2378, NIBS, Beijing, China), (EN1 Polyclonal Antibody, Invitrogen, Waltham, MA, USA), (EN1-4D9, University of California, USA), (EN1-4G11, Columbia University, USA), (EN1, Abcam, Cambridge, MA, USA), (EN1, Sigma-Aldrich, USA), (EN1, Atlas Antibodies, Stockholm, Sweden), (EN1, Cusabio, Wuhan, China) (only EN1-2377 recognizes the correct size of EN1 band), (GAPDH, Origene, USA), and (TULP3, Abcam, Cambridge, MA, USA), (Gli1, Cell Signaling Technology, Danvers, MA, USA), before washing three times with TBST. The membrane was then incubated (1 h) with horseradish peroxidase-conjugated secondary goat anti-rabbit or anti-mouse antibodies (Cell Signaling Technology, Danvers, MA, USA), before being washed again with TBST (Table S3). Signals were detected using an ECL solution (Thermo Fisher, Waltham, MA, USA).

4.5. Cell Proliferation Assays

The effect of EN1 on cell proliferation was measured using a CellTiter-Glo assay kit (Promega, Madison, WI, USA). The glioblastoma cells were seeded at a density of 2500 cells in 96-well plates and incubated. CellTiter-Glo® reagents were mixed and added to the wells after 24/72/120 h. After 2 min of mixing and 10 min of incubation at room temperature, luminescence was recorded on an EnSpire® Multimode Plate Reader (PerkinElmer, Waltham, MA, USA). The luminescence was normalized to the value assessed on day 0, and it is presented as the mean to show the relative proliferation rates of cells after different incubation times. Each experiment was performed in triplicate.
For the EdU incorporation assay, the 5 × 104 scramble control and EN1-KD cells were cultured in six-well plates for 16 h at 37 °C and then exposed to 20 µM 5-ethynyl-2′-deoxyuridine (EdU, Beyotime, Beijing, China) for an additional 3 h. The cells were fixed with 4% formaldehyde for 20 min and permeabilized with 0.5% Triton X-100 at room temperature for 10 min. After washing three times with PBS, Click Reaction Buffer (100 µL) was added to each well and incubated for 30 min, while being protected from light. The cells were subsequently stained with Hoechst 33342 for 10 min, washed three times with PBS, and visualized under a Nikon fluorescent microscope (Nikon, Tokyo, Japan). The EdU-positive cells (green cells) were counted using NIS-Elements BR 3.0 software (Nikon, Tokyo, Japan). The EdU incorporation rate was expressed as the ratio of the number of EdU-positive cells (green cells) to the total number of Hoechst 33342-positive cells (blue cells).

4.6. Wound-Healing Assay

The effect of gene silencing on cell proliferation was measured using a wound-healing assay [46]. Glioblastoma cells were seeded in six-well plates and incubated to reach 90% confluence. Monolayers were scratched with a sterile 200 µL pipette tip to create a wound, and cells were then washed twice with PBS to remove floating cells. Subsequently, cells were cultured in 2% serum DMEM medium for an additional 48 h. Images were captured using a microscope at 0 and 24/48 h, and they were analyzed using the wound-healing macro of ImageJ. The rate of wound closure was calculated as a function of the area of the gap.

4.7. Colony Formation Assays

To evaluate the effect of EN1 in the response to radiotherapy, glioblastoma cells were treated using a γ-ray generator Gammacell 1000 (MDS Nordion, Ottawa, ON, Canada) with doses ranging 0–2 Gy. After different doses of γ-ray irradiation, the glioblastoma cells were suspended and plated in 6 cm plates at 1000 cells/plate for colony formation. Two weeks later, the colonies were fixed in paraformaldehyde (PFA) 4%, then stained with 0.5% crystal violet for 30 min, and photographed. The colonies with more than 50 cells were identified and enumerated using an inverted microscope. Survival was calculated as the average number of colonies counted divided by the number of cells plated multiplicated by plating efficiency (PE), where plating efficiency (PE) is defined as the number of colonies observed/the number of cells plated without radiation [47].

4.8. ROS Level

The scramble control and EN1-KD cell intracellular ROS levels were measured and analyzed by the DCFDA/H2DCFDA assay. Cells were harvested and trypsinized as a single cell, washed with DPBS to remove traces of the original medium, resuspended in 10 μM H2DCFDA (Thermo Fisher Scientific, Waltham, MA, USA) in DMEM medium without serum, and stained for 25 min in the dark at 37 °C. After staining, cells were washed twice with PBS and then resuspended in PBS. H2DCFDA fluoresce signaling was detected by flow cytometer (BD Biosciences, San Jose, CA, USA) at the 488 nm channel, and a total of 10,000 cells were analyzed per sample. Data analysis was performed with FlowJo 7.6 software. Relative fluorescence was calculated by normalization to the control group.

4.9. Immunofluorescence

For immunostaining, the glioma cells were seeded in six-well plates with coverslips (105 cell/well) overnight and fixed in 4% paraformaldehyde. After permeabilizing with 0.25% Triton X-100 for 10 min, cells were blocked in 10% goat serum and incubated with primary antibodies (acetylated tubulin (Lys40), Mouse mAb; ARL13B, Rabbit pAb) (overnight, 4 °C) and washed (30 min) with 0.1% Triton X-100 in PBS. The cells were then incubated for 2 h at room temperature with a secondary Alexa Fluor™ 594 goat anti-mouse IgG (H + L) antibody or 488 goat anti-ribbit IgG (H + L) antibody (Invitrogen, Waltham, MA, USA) and washed again (30 min) with 0.1% Triton X-100 in PBS. Slides were mounted using Prolong Gold antifade with DAPI. Images were acquired using a Nikon fluorescent microscope (NIKON, Tokyo, Japan). The cilia positive cells were expressed as the ratio of the number of cilia-positive cells (green + red cells) to the total number of Hoechst 33342-positive cells (blue cells), counting at least 100 cells
Due to antibody limitations, we could not perform the experiments of detect EN1 expression by IHC, Meanwhile, we checked the single-cell RNA-seq data of glioblastoma [48] and confirmed that EN1 and EN2 was mostly expressed by the tumor cells, and that the EN1 level was higher than EN2 (Figure S7).

4.10. RNA-Seq

Total RNA was isolated from cell samples with Trizol reagent according to the manufacturer’s protocol and sequenced by GENEWIZ (China). Reads were mapped to GRCh37.p13 genome and annotated with STAR, using the ENCODES recommended arguments. Gene-level read counts generated by STAR were used for differential expression analysis with the DESeq2 package. Log2 fold change and adjusted p-values of genes generated by DESeq2 were used for visualization using R.

4.11. Bioinformatics Analysis

TCGA glioma dataset was retrieved from Broad GDAC Firehose (https://gdac.broadinstitute.org/, data version 28 January 2016), and the Henry Ford Hospital dataset was retrieved from the Gene Expression Omnibus (GEO) with accession number GSE4290 following by quantile normalization and log2 transformation with R software. The dataset with 525 glioma cancer samples and 148 normal glioma samples was retrieved from GEO (GSE16011), which were used to detect the expression of EN1 and EN2 in glioma samples. Next, we verified the expression of EN1 in glioma tissues in CGGA (CGGA, Chinese Glioma Genome Atlas; http://cgga.org.cn, accessed on 8 December 2021, DataSet ID: mRNAseq_325) (325 samples) and the Henry Ford Hospital datasets. The expression of EN1 was evaluated with means and standard errors of the mean (SEM) with Graphpad Prism Software 8.0 (San Diego, CA, USA). The power of EN1 to differentiate glioma cancer from normal tissues was evaluated according to ROC curves. We performed correlation analysis to screen co-expression genes of EN1 in glioma patients, with screening criteria as follows: p < 0.05 and |Pearson correlation coefficient| ≥ 0.3. To determine how EN1 affected the prognosis of glioma patients, we performed Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of EN1 co-expressed genes.

4.12. Animal Studies

All animal experiments were conducted following the National Guidelines for Housing and Care of Laboratory Animals in China and performed under the approved Institutional Animal Care and Use Committee protocols.
For mouse tumor models, U251-shcon, U251-shEN1-3898, or U251-shEN1-3899 cells were administered intracranially into the 6–8 week old female BALB/c nude mice (n = 5). The cells were suspended in phosphate-buffered saline (Thermo Scientific, Waltham, MA, USA): Matrigel (Corning, NY, USA) = 1:1 culture medium (5 × 105 cells/5 µL) and injected with a glass micropipette at a rate of 1 µL/min with a Hamilton syringe (10 µL) and a stereotaxic microinjector (RWD, Shenzhen, China). The animals were then monitored daily, and mice were killed when they had a weight decrease >30%. The survival of each mouse was recorded and overall survival was calculated.

4.13. Statistical Analyses

All statistical analyses in this study were conducted using GraphPad Prism 8.0 software (GraphPad, USA). For experimental work, all assays were performed in triplicate, and the measurement data are presented as the mean ± standard deviation (SD). Statistical comparisons between two groups were based on Student’s t-tests, with multiple t-test used for multigroup comparisons. A p-value < 0.05 was considered statistically significant.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms23031123/s1.

Author Contributions

J.C. (Jinchun Chang) and J.C. (Jian Chen) conceptualized and designed the experiments; J.C. (Jinchun Chang), C.G., J.L., Z.L. and R.L. performed the experiments; J.C. (Jinchun Chang) and Y.W. analyzed the data; A.Y. contributed reagents/materials; J.C. (Jinchun Chang) wrote the original draft; Y.G. reviewed and edited the manuscript; J.C. (Jian Chen) and S.H. supervised the research and the preparation of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key R&D Program of China (2016YFA0503100), the Major Project for Natural Science research in Jiangsu Province (17KJA320004), the Jiangsu Basic Research Funds for Young Scientists BK2015035, and institutional grants from the Chinese Ministry of Science and Technology and Beijing Municipal Commission of Science and Technology to the National Institute of Biological Sciences (NIBS).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of Chinese Institute for Brain Research.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are contained within this manuscript.

Acknowledgments

We thank Yun Zhao at the Shanghai Institute of Biochemistry and Cell Biology for providing the hGli1 plasmid, as well as John Snyder for editing our manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Medikonda, R.; Dunn, G.; Rahman, M.; Fecci, P.; Lim, M. A review of glioblastoma immunotherapy. J. Neuro-Oncol. 2021, 151, 41–53. [Google Scholar] [CrossRef] [PubMed]
  2. Hanif, F.; Muzaffar, K.; Perveen, K.; Malhi, S.M.; Simjee, S.U. Glioblastoma Multiforme: A Review of its Epidemiology and Pathogenesis through Clinical Presentation and Treatment. Asian Pac. J. Cancer Prev. 2017, 18, 3–9. [Google Scholar] [PubMed]
  3. Van Meir, E.G.; Hadjipanayis, C.G.; Norden, A.D.; Shu, H.K.; Wen, P.Y.; Olson, J.J. Exciting New Advances in Neuro-Oncology: The Avenue to a Cure for Malignant Glioma. CA A Cancer J. Clin. 2010, 60, 166–193. [Google Scholar] [CrossRef] [PubMed]
  4. Davis, M. Glioblastoma: Overview of Disease and Treatment. Clin. J. Oncol. Nurs. 2016, 20, S2–S8. [Google Scholar] [CrossRef] [Green Version]
  5. Osuka, S.; Van Meir, E.G. Overcoming therapeutic resistance in glioblastoma: The way forward. J. Clin. Investig. 2017, 127, 415–426. [Google Scholar] [CrossRef] [Green Version]
  6. García-Bellido, A.; Santamaría, P. Developmental analysis of the wing disc in the mutant engrailed of Drosophila melanogaster. Genetics 1972, 72, 87–104. [Google Scholar] [CrossRef] [PubMed]
  7. dos Santos, M.T.A.; Smidt, M.P. En1 and Wnt signaling in midbrain dopaminergic neuronal development. Neural Dev. 2011, 6, 23. [Google Scholar] [CrossRef] [Green Version]
  8. Orvis, G.D.; Hartzell, A.L.; Smith, J.B.; Barraza, L.H.; Wilson, S.L.; Szulc, K.U.; Turnbull, D.H.; Joyner, A.L. The engrailed homeobox genes are required in multiple cell lineages to coordinate sequential formation of fissures and growth of the cerebellum. Dev. Biol. 2012, 367, 25–39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Willett, R.T.; Bayin, N.S.; Lee, A.S.; Krishnamurthy, A.; Wojcinski, A.; Lao, Z.; Stephen, D.; Rosello-Diez, A.; Dauber-Decker, K.L.; Orvis, G.D.; et al. Cerebellar nuclei excitatory neurons regulate developmental scaling of presynaptic Purkinje cell number and organ growth. eLife 2019, 8, e50617. [Google Scholar] [CrossRef]
  10. Hanks, M.; Wurst, W.; Anson-Cartwright, L.; Auerbach, A.B.; Joyner, A.L. Rescue of the En-1 mutant phenotype by replacement of En-1 with En-2. Science 1995, 269, 679–682. [Google Scholar] [CrossRef]
  11. Kim, Y.J.; Sung, M.; Oh, E.; Vrancken, M.V.; Song, J.Y.; Jung, K.; Choi, Y.L. Engrailed 1 overexpression as a potential prognostic marker in quintuple-negative breast cancer. Cancer Biol. Ther. 2018, 19, 335–345. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Beltran, A.S.; Graves, L.M.; Blancafort, P. Novel role of Engrailed 1 as a prosurvival transcription factor in basal-like breast cancer and engineering of interference peptides block its oncogenic function. Oncogene 2013, 33, 4767–4777. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Do Carmo Silva, J.; Vesely, S.; Novak, V.; Luksanova, H.; Prusa, R.; Babjuk, M. Engrailed 2 in urine fails as a marker in prostate cancer detection. Eur. Urol. Suppl. 2017, 16, e2646. [Google Scholar] [CrossRef]
  14. Punia, N.; Primon, M.; Simpson, G.R.; Pandha, H.S.; Morgan, R. Membrane insertion and secretion of the Engrailed-2 (EN2) transcription factor by prostate cancer cells may induce antiviral activity in the stroma. Sci. Rep. 2019, 9, 5138. [Google Scholar] [CrossRef] [PubMed]
  15. Li, Y.; Liu, J.; Xiao, Q.; Tian, R.; Zhou, Z.; Gan, Y.; Li, Y.; Shu, G.; Yin, G. EN2 as an oncogene promotes tumor progression via regulating CCL20 in colorectal cancer. Cell Death Dis. 2020, 11, 604. [Google Scholar] [CrossRef] [PubMed]
  16. Michael, A.; McGrath, S.; Annels, N.; Denyer, M.; Morgan, R.; Pandha, H. 438 Engrailed 2 protein (EN2) as a novel biomarker in epithelial ovarian cancer. Eur. J. Cancer 2015, 51, S91. [Google Scholar] [CrossRef]
  17. McGrath, S.E.; Annels, N.; Madhuri, T.K.; Tailor, A.; Butler-Manuel, S.A.; Morgan, R.; Pandha, H.; Michael, A. Engrailed-2 (EN2)—A novel biomarker in epithelial ovarian cancer. BMC Cancer 2018, 18, 943. [Google Scholar] [CrossRef] [Green Version]
  18. Martin, N.L.; Saba-El-Leil, M.K.; Sadekova, S.; Meloche, S.; Sauvageau, G. EN2 is a candidate oncogene in human breast cancer. Oncogene 2005, 24, 6890–6901. [Google Scholar] [CrossRef] [Green Version]
  19. Skoda, A.M.; Simovic, D.; Karin, V.; Kardum, V.; Vranic, S.; Serman, L. The role of the Hedgehog signaling pathway in cancer: A comprehensive review. Bosn. J. Basic Med. Sci. 2018, 18, 8–20. [Google Scholar] [CrossRef]
  20. Chen, M.H.; Wilson, C.W.; Li, Y.J.; Law, K.K.; Lu, C.S.; Gacayan, R.; Zhang, X.; Hui, C.C.; Chuang, P.T. Cilium-independent regulation of Gli protein function by Sufu in Hedgehog signaling is evolutionarily conserved. Genes Dev. 2009, 23, 1910–1928. [Google Scholar] [CrossRef] [Green Version]
  21. Jeng, K.S.; Chang, C.F.; Lin, S.S. Sonic Hedgehog Signaling in Organogenesis, Tumors, and Tumor Microenvironments. Int. J. Mol. Sci. 2020, 21, 758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Yao, M.; Ventura, P.B.; Jiang, Y.; Rodriguez, F.J.; Wang, L.; Perry, J.S.A.; Yang, Y.; Wahl, K.; Crittenden, R.B.; Bennett, M.L.; et al. Astrocytic trans-Differentiation Completes a Multicellular Paracrine Feedback Loop Required for Medulloblastoma Tumor Growth. Cell 2020, 180, 502–520.e19. [Google Scholar] [CrossRef] [PubMed]
  23. Santoni, M.; Burattini, L.; Nabissi, M.; Beatrice Morelli, M.; Berardi, R.; Santoni, G.; Cascinu, S. Essential Role of Gli Proteins in Glioblastoma Multiforme. Curr. Protein Pept. Sci. 2013, 14, 133–140. [Google Scholar] [CrossRef] [PubMed]
  24. Alcantara Llaguno, S.; Chen, J.; Kwon, C.H.; Jackson, E.L.; Li, Y.; Burns, D.K.; Alvarez-Buylla, A.; Parada, L.F. Malignant astrocytomas originate from neural stem/progenitor cells in a somatic tumor suppressor mouse model. Cancer Cell 2009, 15, 45–56. [Google Scholar] [CrossRef] [Green Version]
  25. Chen, J.; Li, Y.; Yu, T.S.; McKay, R.M.; Burns, D.K.; Kernie, S.G.; Parada, L.F. A restricted cell population propagates glioblastoma growth after chemotherapy. Nature 2012, 488, 522–526. [Google Scholar] [CrossRef] [Green Version]
  26. Auffinger, B.; Spencer, D.; Pytel, P.; Ahmed, A.U.; Lesniak, M.S. The role of glioma stem cells in chemotherapy resistance and glioblastoma multiforme recurrence. Expert Rev. Neurother. 2015, 15, 741–752. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Gusev, Y.; Bhuvaneshwar, K.; Song, L.; Zenklusen, J.C.; Fine, H.; Madhavan, S. The REMBRANDT study, a large collection of genomic data from brain cancer patients. Sci. Data 2018, 5, 180158. [Google Scholar] [CrossRef] [PubMed]
  28. Zhao, Z.; Zhang, K.N.; Wang, Q.; Li, G.; Zeng, F.; Zhang, Y.; Wu, F.; Chai, R.; Wang, Z.; Zhang, C.; et al. Chinese Glioma Genome Atlas (CGGA): A Comprehensive Resource with Functional Genomic Data from Chinese Glioma Patients. Genom. Proteom. Bioinform. 2021, 19, 1–12. [Google Scholar] [CrossRef] [PubMed]
  29. Mittler, R. ROS Are Good. Trends. Plant Sci. 2017, 22, 11–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Lindell Jonsson, E.; Erngren, I.; Engskog, M.; Haglof, J.; Arvidsson, T.; Hedeland, M.; Petterson, C.; Laurell, G.; Nestor, M. Exploring Radiation Response in Two Head and Neck Squamous Carcinoma Cell Lines Through Metabolic Profiling. Front. Oncol. 2019, 9, 825. [Google Scholar] [CrossRef] [Green Version]
  31. Azzam, E.I.; Jay-Gerin, J.P.; Pain, D. Ionizing radiation-induced metabolic oxidative stress and prolonged cell injury. Cancer Lett. 2012, 327, 48–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Yamamori, T.; Yasui, H.; Yamazumi, M.; Wada, Y.; Nakamura, Y.; Nakamura, H.; Inanami, O. Ionizing radiation induces mitochondrial reactive oxygen species production accompanied by upregulation of mitochondrial electron transport chain function and mitochondrial content under control of the cell cycle checkpoint. Free Radic. Biol. Med. 2012, 53, 260–270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Singer, E.; Judkins, J.; Salomonis, N.; Matlaf, L.; Soteropoulos, P.; McAllister, S.; Soroceanu, L. Reactive oxygen species-mediated therapeutic response and resistance in glioblastoma. Cell Death Dis. 2015, 6, e1601–e1611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Gupta, S.; Takebe, N.; Lorusso, P. Targeting the Hedgehog pathway in cancer. Ther. Adv. Med. Oncol. 2010, 2, 237–250. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Carpenter, R.L.; Lo, H.W. Hedgehog pathway and GLI1 isoforms in human cancer. Discov. Med. 2012, 13, 105. [Google Scholar] [PubMed]
  36. Badgandi, H.B.; Hwang, S.H.; Shimada, I.S.; Loriot, E.; Mukhopadhyay, S. Tubby family proteins are adapters for ciliary trafficking of integral membrane proteins. J. Cell Biol. 2017, 216, 743–760. [Google Scholar] [CrossRef] [Green Version]
  37. Bangs, F.; Anderson, K.V. Primary cilia and mammalian hedgehog signaling. Cold Spring Harb. Perspect. Biol. 2017, 9, a028175. [Google Scholar] [CrossRef]
  38. Mukhopadhyay, S.; Wen, X.; Ratti, N.; Loktev, A.; Rangell, L.; Scales, S.J.; Jackson, P.K. The ciliary G-protein-coupled receptor Gpr161 negatively regulates the Sonic hedgehog pathway via cAMP signaling. Cell 2013, 152, 210–223. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Norman, R.X.; Ko, H.W.; Huang, V.; Eun, C.M.; Abler, L.L.; Zhang, Z.; Sun, X.; Eggenschwiler, J.T. Tubby-like protein 3 (TULP3) regulates patterning in the mouse embryo through inhibition of Hedgehog signaling. Hum. Mol. Genet. 2009, 18, 1740–1754. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Bell, D.; Bell, A.; Roberts, D.; Weber, R.S.; El-Naggar, A.K. Developmental transcription factor EN1- novel biomarker in human salivary gland adenoid cystic carcinoma. Cancer 2012, 118, 1288–1292. [Google Scholar] [CrossRef] [Green Version]
  41. Bachar-Dahan, L.; Goltzmann, J.; Yaniv, A.; Gazit, A. Engrailed-1 negatively regulates beta-catenin transcriptional activity by destabilizing beta-catenin via a glycogen synthase kinase-3beta-independent pathway. Mol. Biol. Cell 2006, 17, 2572–2580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Clement, V.; Sanchez, P.; de Tribolet, N.; Radovanovic, I.; Ruiz i Altaba, A. HEDGEHOG-GLI1 Signaling Regulates Human Glioma Growth, Cancer Stem Cell Self-Renewal, and Tumorigenicity. Curr. Biol. 2007, 17, 165–172. [Google Scholar] [CrossRef] [PubMed]
  43. Lee, D.H.; Lee, S.Y.; Oh, S.C. Hedgehog signaling pathway as a potential target in the treatment of advanced gastric cancer. Tumor Biol. 2017, 39, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Bilovocky, N.A.; Romito-DiGiacomo, R.R.; Murcia, C.L.; Maricich, S.M.; Herrup, K. Factors in the genetic background suppress the engrailed-1 cerebellar phenotype. J. Neurosci. 2003, 23, 5105–5112. [Google Scholar] [CrossRef]
  45. Tabata, T.; Eaton, S.; Kornberg, T.B. The Drosophila hedgehog gene is expressed specifically in posterior compartment cells and is a target of engrailed regulation. Genes Dev. 1992, 6, 2635–2645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Jonkman, J.E.; Cathcart, J.A.; Xu, F.; Bartolini, M.E.; Amon, J.E.; Stevens, K.M.; Colarusso, P. An introduction to the wound healing assay using live-cell microscopy. Cell Adh. Migr. 2014, 8, 440–451. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Franken, N.A.; Rodermond, H.M.; Stap, J.; Haveman, J.; van Bree, C. Clonogenic assay of cells in vitro. Nat. Protoc. 2006, 1, 2315–2319. [Google Scholar] [CrossRef]
  48. Darmanis, S.; Sloan, S.A.; Croote, D.; Mignardi, M.; Chernikova, S.; Samghababi, P.; Zhang, Y.; Neff, N.; Kowarsky, M.; Caneda, C. Single-cell RNA-seq analysis of infiltrating neoplastic cells at the migrating front of human glioblastoma. Cell Rep. 2017, 21, 1399–1410. [Google Scholar] [CrossRef] [Green Version]
Ijms 23 01123 g001 550
Figure 1. Elevated EN expression is associated with glioma cancer progression. (A) RNA-seq results for mouse glioma stem cells lines vs. neural stem cells (left) and for mouse brain tumor tissues vs. normal brain tissues (right). (B) The expression of EN1 (left) and EN2 (right) in glioma samples from TCGA database. Black, blue and red color is for Grade 2, 3 and 4 respectively. (C) Analysis using the Henry Ford Hospital dataset, showing correlation between the EN1 (left) and EN2 (right) gene expression levels and glioma patient tumor grade (○ for Non-tumor, □ and ▲ for Oligodendroglioma-2 and 3, ◊ and ∆ for Asreocytoma-2 and 3, red ● for GBM). (D) Correlation between glioma patient survival and the EN1 expression level from TCGA glioblastoma database (left) and the Chinese Glioma genome Atlas gliomas database (right). (E) For the indicated glioma cell lines, Engrailed 1 expression assessed by qPCR (left) and Engrailed 1 accumulation assessed by immunoblotting (right). * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001 and ns (no significant) according to Student’s t-test.
Figure 1. Elevated EN expression is associated with glioma cancer progression. (A) RNA-seq results for mouse glioma stem cells lines vs. neural stem cells (left) and for mouse brain tumor tissues vs. normal brain tissues (right). (B) The expression of EN1 (left) and EN2 (right) in glioma samples from TCGA database. Black, blue and red color is for Grade 2, 3 and 4 respectively. (C) Analysis using the Henry Ford Hospital dataset, showing correlation between the EN1 (left) and EN2 (right) gene expression levels and glioma patient tumor grade (○ for Non-tumor, □ and ▲ for Oligodendroglioma-2 and 3, ◊ and ∆ for Asreocytoma-2 and 3, red ● for GBM). (D) Correlation between glioma patient survival and the EN1 expression level from TCGA glioblastoma database (left) and the Chinese Glioma genome Atlas gliomas database (right). (E) For the indicated glioma cell lines, Engrailed 1 expression assessed by qPCR (left) and Engrailed 1 accumulation assessed by immunoblotting (right). * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001 and ns (no significant) according to Student’s t-test.
Ijms 23 01123 g001
Ijms 23 01123 g002 550
Figure 2. EN1 knockdown restrains the proliferation and migration of glioma cells. (A) Assessment of knockdown efficiency for EN1 with immunoblotting, and confirmation of EN1 knockout in U251 cells with immunoblotting and genome sequencing. Cell proliferation was assessed using CellTiter-Glo assays (B), 5-ethynyl-2′-deoxyuridine (Edu) assays (C), and colony formation assays (D) (● for shcon, ■ for shEN1-3898, ▲ for shEN1-3899, ◆ for shEN1-3901). (E) Cell migration was assessed with a wound-healing assay. (F) Quantification of wound-healing assay. (G) Survival curves for mice bearing xenograft tumors derived from shRNA-EN1 cells or scramble control cells. (H) Cell proliferation was assessed by CellTiter-Glo and colony formation assays in U-87 MG cells after EN1 overexpression. Scale bar: 100 μm. Data are presented as the mean ± SEM from three independent experiments. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001 according to Student’s t-test.
Figure 2. EN1 knockdown restrains the proliferation and migration of glioma cells. (A) Assessment of knockdown efficiency for EN1 with immunoblotting, and confirmation of EN1 knockout in U251 cells with immunoblotting and genome sequencing. Cell proliferation was assessed using CellTiter-Glo assays (B), 5-ethynyl-2′-deoxyuridine (Edu) assays (C), and colony formation assays (D) (● for shcon, ■ for shEN1-3898, ▲ for shEN1-3899, ◆ for shEN1-3901). (E) Cell migration was assessed with a wound-healing assay. (F) Quantification of wound-healing assay. (G) Survival curves for mice bearing xenograft tumors derived from shRNA-EN1 cells or scramble control cells. (H) Cell proliferation was assessed by CellTiter-Glo and colony formation assays in U-87 MG cells after EN1 overexpression. Scale bar: 100 μm. Data are presented as the mean ± SEM from three independent experiments. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001 according to Student’s t-test.
Ijms 23 01123 g002
Ijms 23 01123 g003 550
Figure 3. EN1 knockdown increased intracellular ROS levels and elevated susceptibility to γ-ray irradiation. (A) EN1 knockdown elevated the intracellular ROS levels of glioma cells, assessed with H2DCFDA assays (● for shcon, ■ for shEN1-3898, ▲ for shEN1-3899, ◆ for shEN1-3901). (B) Colony forming assays demonstrating that EN1 knockdown elevated the susceptibility of glioma cells to γ-ray irradiation, assessed with plating efficiency (PE), where PE is the fraction of colonies counted divided by cells plated without radiation. Data are presented as the mean from three independent experiments. * p < 0.05, ** p < 0.01, *** p < 0.001 according to Student’s t-test.
Figure 3. EN1 knockdown increased intracellular ROS levels and elevated susceptibility to γ-ray irradiation. (A) EN1 knockdown elevated the intracellular ROS levels of glioma cells, assessed with H2DCFDA assays (● for shcon, ■ for shEN1-3898, ▲ for shEN1-3899, ◆ for shEN1-3901). (B) Colony forming assays demonstrating that EN1 knockdown elevated the susceptibility of glioma cells to γ-ray irradiation, assessed with plating efficiency (PE), where PE is the fraction of colonies counted divided by cells plated without radiation. Data are presented as the mean from three independent experiments. * p < 0.05, ** p < 0.01, *** p < 0.001 according to Student’s t-test.
Ijms 23 01123 g003
Ijms 23 01123 g004 550
Figure 4. EN1 increases the activity of proliferation related Hedgehog signaling pathways. (A) Overrepresented Gene Ontology (GO) terms from RNA−seq analysis of upregulated gene sets and downregulated gene sets in EN1 knockdown group compared to control cells. (B) Immunoblotting assays showing that EN1 knockdown significantly reduces the Gli1 level in glioma cells. (C) The EN1 and Gli1 transcript levels in U251 cells transduced with scramble control, shEN1, shcon + Gli1, or shEN1 + Gli1. (D) Overexpression of Gli1 rescues the restrained colony-forming abilities caused by the knockout of EN1 (□ for U251-sgLacZ, ● for U251-sgEN1, ■ for U251+Gli1-sgEN1). Data are presented as the mean from three independent experiments. * p < 0.05, ### p < 0.001 according to Student’s t-test.
Figure 4. EN1 increases the activity of proliferation related Hedgehog signaling pathways. (A) Overrepresented Gene Ontology (GO) terms from RNA−seq analysis of upregulated gene sets and downregulated gene sets in EN1 knockdown group compared to control cells. (B) Immunoblotting assays showing that EN1 knockdown significantly reduces the Gli1 level in glioma cells. (C) The EN1 and Gli1 transcript levels in U251 cells transduced with scramble control, shEN1, shcon + Gli1, or shEN1 + Gli1. (D) Overexpression of Gli1 rescues the restrained colony-forming abilities caused by the knockout of EN1 (□ for U251-sgLacZ, ● for U251-sgEN1, ■ for U251+Gli1-sgEN1). Data are presented as the mean from three independent experiments. * p < 0.05, ### p < 0.001 according to Student’s t-test.
Ijms 23 01123 g004
Ijms 23 01123 g005 550
Figure 5. Silencing EN1 upregulates the expression of TULP3. (A) Volcano plot of differentially expressed genes (DEGs) in EN1 knockdown glioma cells (vs. scramble control cells). qPCR (B) and immunoblotting (C) confirmed that EN1 knockdown glioma cells have elevated TULP3. (D) Staining of glioma cells with antibodies against ARL13B and acetylated tubulin. (E) EN1 knockdown significantly decreased the number of primary cilia. Scale bar: 50 µm. Data are presented as the mean from three independent experiments (● for shcon, ■ for shEN1-3898, ▲ for shEN1-3899, ◆ for shEN1-3901). * p < 0.05 according to Student’s t-test.
Figure 5. Silencing EN1 upregulates the expression of TULP3. (A) Volcano plot of differentially expressed genes (DEGs) in EN1 knockdown glioma cells (vs. scramble control cells). qPCR (B) and immunoblotting (C) confirmed that EN1 knockdown glioma cells have elevated TULP3. (D) Staining of glioma cells with antibodies against ARL13B and acetylated tubulin. (E) EN1 knockdown significantly decreased the number of primary cilia. Scale bar: 50 µm. Data are presented as the mean from three independent experiments (● for shcon, ■ for shEN1-3898, ▲ for shEN1-3899, ◆ for shEN1-3901). * p < 0.05 according to Student’s t-test.
Ijms 23 01123 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chang, J.; Guo, C.; Li, J.; Liang, Z.; Wang, Y.; Yu, A.; Liu, R.; Guo, Y.; Chen, J.; Huang, S. EN1 Regulates Cell Growth and Proliferation in Human Glioma Cells via Hedgehog Signaling. Int. J. Mol. Sci. 2022, 23, 1123. https://doi.org/10.3390/ijms23031123

AMA Style

Chang J, Guo C, Li J, Liang Z, Wang Y, Yu A, Liu R, Guo Y, Chen J, Huang S. EN1 Regulates Cell Growth and Proliferation in Human Glioma Cells via Hedgehog Signaling. International Journal of Molecular Sciences. 2022; 23(3):1123. https://doi.org/10.3390/ijms23031123

Chicago/Turabian Style

Chang, Jinchun, Chenjia Guo, Jianyu Li, Zhangqian Liang, Yankai Wang, Anliang Yu, Runze Liu, Yuting Guo, Jian Chen, and Song Huang. 2022. "EN1 Regulates Cell Growth and Proliferation in Human Glioma Cells via Hedgehog Signaling" International Journal of Molecular Sciences 23, no. 3: 1123. https://doi.org/10.3390/ijms23031123

APA Style

Chang, J., Guo, C., Li, J., Liang, Z., Wang, Y., Yu, A., Liu, R., Guo, Y., Chen, J., & Huang, S. (2022). EN1 Regulates Cell Growth and Proliferation in Human Glioma Cells via Hedgehog Signaling. International Journal of Molecular Sciences, 23(3), 1123. https://doi.org/10.3390/ijms23031123

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop