Key Genes Are Associated with the Prognosis of Glioma, and Melittin Can Regulate the Expression of These Genes in Glioma U87 Cells

Li, Ran; Tao, Ting; Ren, Qiuyun; Xie, Sujun; Gao, Xiaofen; Wu, Jie; Chen, Diling; Xu, Changqiong

doi:https://doi.org/10.1155/2022/7033478

BioMed Research International

On this page

Abstract Methods Results Discussion Conclusion Abbreviations Data Availability Consent Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 7033478 | https://doi.org/10.1155/2022/7033478

Key Genes Are Associated with the Prognosis of Glioma, and Melittin Can Regulate the Expression of These Genes in Glioma U87 Cells

Ran Li,^1,2,3Ting Tao,³Qiuyun Ren,⁴Sujun Xie,⁵Xiaofen Gao,²Jie Wu,⁴Diling Chen,⁶and Changqiong Xu^2,3

Academic Editor: Quintino Giorgio D'alessandris

Received07 May 2022

Revised04 Nov 2022

Accepted02 Dec 2022

Published31 Dec 2022

Abstract

Glioma is the most common primary tumor of the central nervous system. Currently, there is no effective treatment for glioma. Melittin (MT) is the main component of bee venom, which was found to have therapeutic effects on a variety of tumors. In this study, we explored the relationship between key genes regulated by MT and the prognosis of glioma. In cultured glioma U87 and U251 cells, MT inhibited cell proliferation and induces cell apoptosis in a time- and concentration-dependent manner. RNA-seq revealed that MT upregulated 11 genes and downregulated 37 genes. These genes are mainly enriched in cell membrane signaling pathways, such as surface membrane, membrane-enclosed organelles, integral component of membrane, PPAR signaling pathway, and voltage-gated potassium channel. PPI network analysis and literature analysis of 48 genes were performed, and 8 key genes were identified, and these key genes were closely associated with clinical prognosis. Overexpression of PCDH18, PPL, DEPP1, VASN, KCNE4, MYBPH, and C5AR2 genes or low expression of MARCH4 gene in glioma patients was associated with poor survival. qPCR confirmed that MT can regulate the expression of these genes in glioma U87 cells. This study indicated that MT significantly inhibited the growth and regulated the expression of PCDH18, C5AR2, VASN, DEPP1, MYBPH, KCNE4, PPL, and MARCH4 genes in glioma U87 cells in vitro. These genes are closely related to the prognosis of patients with glioma and can be used as independent prognostic factors in patients with glioma. MT is a potential drug for the treatment of glioma.

1. Introductions

Glioma is the most common primary malignant tumor of the brain. It has high morbidity and low survival rate [1–3]. In the 2021 World Health Organization (WHO) classification of central nervous system (CNS) tumors, low-grade glioma (LGG) includes CNS WHO grade 1-2, and high-grade glioma (HGG) includes CNS WHO grade 3-4. LGG accounts for 6% of primary central nervous system tumors in adults and has a good prognosis [4]. However, compared with other benign intracranial tumors, there is still a higher recurrence rate. Glioblastoma multiforme (GBM) is WHO grade 4, which is the most aggressive and malignant primary brain tumor. Surgical is an important treatment for LGG. However, due to the aggressive growth of glioma, surgical treatment alone cannot completely cure the diffuse growth of the tumor. At present, adjuvant therapies such as radiotherapy and chemotherapy are mostly used to delay the tumor recurrence time. But LGG are still difficult to cure completely. For HGG, the current standard treatment includes surgical removal of the tumor [5], tumor-treating fields (TT Fields) [6], targeted therapy [7], immunotherapy [8], radiotherapy (RT) [9], and chemotherapy [10]. Despite the continuous updates and advances in diagnostic and therapeutic techniques, the prognosis of patients is not satisfactory. In recent years, TMZ has made some progress as the main drug in chemotherapy, but the prognosis of patients is still poor [11]. It is urgent to find effective treatment methods. Therefore, searching for effective treatment is still a hot spot in recent years.

MT is the main active component of honeybee venom, honeybee venom is a kind of fragrant transparent venom secreted by worker bee venom glands and accessory glands (Figure 1(a)). Among the many components of honeybee venom, the content of MT is the highest which makes up more than 40% of honeybee venom, MT is a polypeptide composed of 26 amino acid residues [12] (Figure 1(b)), the molecular formula is C131H2 29N39O31, and the relative molecular weight is 2846.46.

(a)

(b)

In recent years, there have been many studies on its pharmacological action and mechanisms. It has been reported that MT has antioxidant [13], antifungal [14], and anti-inflammatory [15, 16] pharmacological effects. Therefore, MT was used in the treatment of a variety of diseases, especially rheumatoid arthritis [17]. The mechanism may be related to regulating Th17/Treg balance [18] and inhibiting the activation of NF-κB, STAT3, and bcl-2 expression induced by IL-6/S IL-6R complex [19]. In addition, previous studies have shown that MT destroys the cell membrane and leads to cell death. It does not need to enter the cells, but can show cytotoxicity to destroy tumor cells outside the cells [20]. MT is believed to have a proapoptotic effect and antitumor activity [21] by inhibition of the HIF-1 α/Akt pathway in liver cancer [22], NF-κB signaling pathway in lung carcinoma cells [23], TGF-β-mediated ERK signaling pathway in lung cancer [24], Her2 enrichment and growth factor receptor activation in triple negative breast cancer [25], MAPK in melanoma [26], and JAK2/STAT3 in ovarian cancer [27]. In addition, some studies [28–30] have shown that MT inhibits proliferation and induces apoptosis of malignant human glioma cells. The majority of the antineoplastic activity of honeybee venom has been attributed to melittin through inhibition of VEGF, FLT-1, and MMP-9 in glioma C6 cells [31]; metaloprotease-2 in human glioblastoma cells [32]; and STAT3 and VEGF in glioma SHG44 cell [33]. However, information regarding the functional role of MT in glioma is limited. The targets and mechanism of action of MT in the treatment of glioma need to be further studied.

Therefore, on the basis of the previous studies, this study further explored the effect of MT on glioma. In this study, we evaluated the potential of MT in inhibiting U87 and U251 cells in vitro. We generated deep sequencing RNA data from U87- and MT-treated U87 cell samples, monitoring the differentially expressed genes. TCGA database was used to further screen genes related to the survival prognosis of glioma. We found that MT can regulate the expression of these genes associated with glioma prognosis in glioma U87 cells. That could be helpful for determining the molecular mechanism of MT in the treatment of glioma.

2. Methods

The methods were followed as previously described [34] with some modifications.

2.1. MT Preparation

MT (20449-79-0, Shanghai Baishun Biological Technology Co., Ltd., China) was dissolved in saline and stocked at -20°C.

2.2. Cell Culture

U251 and U87 cells were derived from glioblastoma tissue of patients. U87 cells were cultured from a grade III astrocytoma-glioblastoma of a 44-year-old woman. U251 cells were cultured from a 75-year-old patient. U87 cells (provided by professor Yi Guan) and U251 cells (Suzhou Institute of Cell Biology, Chinese Academy of Sciences) were cultured in Dulbecco modified Eagle medium (Gibco, USA) with 10% FBS (Gibco, USA) and 1% of penicillin-streptomycin at 37°C, in humidified air containing 5% of CO₂.

2.3. Cell Viability

Cell proliferation was determined by using a Cell Counting Kit-8 (CCK-8, Meilunbio, Dalian, China). Cells were seeded into a 96-well plate () with 100 μL of culture medium. After 12 h, cells were treated with various concentrations of MT (0, 30.0, 60.0, 150.0, 300.0, 600.0, and 1500.0 nM) for 24 and 48 h. 10 μL CCK-8 solution was added to each well and cells were incubated at 37°C for 2 h. Light absorbance was measured at 450 nm using a microplate reader (Bio Tek Elx800, USA).

The cell growth inhibition and the median inhibitory concentration (IC₅₀) were calculated from the cytotoxicity curves. Cell morphological changes were observed using an inverted microscope (Nikon, ECLIPSE E600, Tokyo, JPN).

2.4. Annexin V-FITC/PI Cell Apoptosis Analysis

Apoptosis was detected by flow cytometry via the examination of altered plasma membrane phospholipid packing by lipophilic dye Annexin V. Briefly, treated cells (treated with MT, 0-150.0 nM) were harvested by trypsin, washed twice with PBS, and then, resuspended in binding buffer. Thereafter, 5 μL of Annexin V-FITC was added into cell suspension and incubated for 15 min at room temperature in the dark. PI staining was performed 5 min before evaluation by flow cytometry within 30 min (BD Accuri™ C6, BD Biosciences, San Jose, CA). Data analysis was performed with BD Accuri C6 software (Version 1.0.264.21).

2.5. Total RNA Isolation and Transcriptome Sequencing

U87 cells were seeded in a 6-well plate ( cells/well). After 12 h, cells were treated with MT (0 μM and 150 nM for 12 h). Total RNA was isolated using a RNA extraction kit (Beyotime, Shanghai, China). Nanodrop 2000 spectrophotometer (Implen, Los Angeles, USA) was used to test RNA purity and concentration. After that, mRNA was enriched, the enriched mRNA was fragmented into short fragments and reverse transcribed into cDNA. Second-strand cDNA was synthesized by DNA polymerase I, RNase H, dNTP, and buffer. Then, the cDNA fragments were purified with a QIAQuick PCR extraction kit, end-repaired, poly(A) tailed, and ligated to Illumina sequencing adapters. The ligation products were size selected by agarose gel electrophoresis, PCR amplified, and sequenced using an Illumina HiSeqTM 2500 platform by Gene Denovo Biotechnology Co. (Guangzhou, China).

2.6. PPI Network Analysis and Key Genes

PPI network analysis plays a major role in predicting the functionality of interacting genes or proteins and gives an insight into the functional relationships and evolutionary conservation of interactions among the genes. An interaction network is a graphical representation of gene/protein interactome, where each gene/protein is a node, and interaction between gene/protein is an edge [35]. PPI was constructed using the String Database (https://string-db.org/), and the key genes were obtained.

2.7. Correlation Analysis between Key Genes and Survival Prognosis

Survival analysis [36] is a statistical method used to analyze the relationship between survival status and corresponding time. It has been widely used in medicine. Survival Rate represents the probability that the survival time of the patient is greater than time , denoted by . Time represents the actual survival time of the patient. The survival rate curve can be obtained if the horizontal axis is taken as the following prevention time and the vertical axis is taken as the survival probability .

Independent prognostic analysis of genes associated with glioma survival was performed to obtain independent prognostic factors associated with glioma. These genes as independent prognostic factors are of great clinical interest which suggested that these genes could be independent of other clinical traits as independent prognostic factors [37].

To further verify the relationship between the key genes and the survival prognosis of patients with glioma, RNA-seq data and clinical data of glioma were taken from The Cancer Genome Atlas (TCGA) database (https://tcga-data.nci.nih.gov/tcga/). According to the expression of DEG mRNA in glioma, DEGs were divided into high-expression group (above average) and low-expression group (below average), and the survival rate of patients was analyzed. Kaplan-Meier curve was used to draw the survival curve to detect the relationship between the expression level of DEGs and the survival of glioma patients. We then assessed the DEGs based on other clinical factors in the univariate and multivariate Cox proportional-hazard analysis.

2.8. Quantitative Real-Time PCR Analysis

Total RNA was isolated for RNA-seq. Based on the RNA-seq results, relevant genes were selected, and primers were designed using Primer Premier 5 software (primer sequences are shown in Table 1). Reverse transcription into cDNA was performed using a cDNA Synthesis SuperMix for qPCR Kit (Transgen, Beijing, China), and the cDNA was amplified using an RT-PCR amplification kit (PerfectStart Green qPCR SuperMix) (Transgen, Beijing, China). The method was used to perform analysis with GAPDH as the reference gene.

2.9. Statistical Analysis

Data were analyzed using GraphPad Prism 7.0 (GraphPad Software, Inc. La Jolla, CA, USA) and SPSS software (v.24.0). Data are presented as . Statistical analyses of multiple group comparisons were performed by two-way ANOVA (CCK8 analysis, IC₅₀, apoptotic rate, qPCR). Clinical correlation analysis was performed using R (version 4.0.4, http://www.r-project.org). Cumulative survival time was calculated using the Kaplan-Meier method, and the differences in survival curves were analyzed using the log-rank test from “survival” package [38] (version: 2.41.3). Univariate and multivariate analyses were conducted using the Cox proportional hazard regression model. For all tests, significance was defined at , , , and .

3. Results

3.1. MT Affects Cell Morphology and Suppresses Cell Proliferation

The human U87 and U251 cells were used to evaluate the pharmacological effects of MT. Cells were treated with different concentrations of MT for 48 h, then the cell morphological changes were examined. At the concentration of 60.0 nM group, cells exhibited morphological changes, such as cell rounding and the disappearance of protruding spike. At the concentration of 150 nM, treated cells showed cell shrinking, disappearance of spikes, and detachment from the substrate (Figure 2(a)).

(a)

(b)

(c)

CCK-8 assay on cell viability showed that MT treatment remarkably inhibited U87 cell growth in a concentration-dependent manner. The IC₅₀ values of MT-induced inhibition at 24 and 48 h were and , respectively. In U251 cells, the IC₅₀ values at 24 and 48 h were and , respectively (Figures 2(b) and 2(c), Table 2).

3.2. MT Induces Cell Apoptosis

To evaluate the effects of MT on cell apoptosis, U87 and U251 cells were treated with different concentrations of MT, then examined for cell apoptosis by Annexin V-FITC/PI flow cytometry (Figure 3(a)). After treatment with 30.0, 60.0, and 150.0 nM MT on U87 cells, the early apoptotic rate increased from (untreated control group) to , then decreased to and , respectively. The decrease of early apoptosis at 60.0 and 150.0 nM MT may be due to the rapid occurrence of late apoptosis and even death of tumor cells. Late apoptotic rate increased from (control) to , , and at 30.0, 60.0, and 150.0 nM MT, respectively. Thus, the total apoptotic rate increased from (control) to , , and at 30.0, 60.0, and 150.0 nM MT, respectively. There was a significant difference between the 60.0 and 150.0 nM ( vs. control) in late and total apoptosis. In U251 cells, we observed that MT had a statistically significant effect on the late apoptosis and total apoptosis of U251 cells at 60 nM ( vs. control), while MT had a statistically significant effect on the late apoptosis and total apoptosis of U251 cells at 150 nM ( vs. control) (Figure 3(b)).

(a)

(b)

Figure 3

The effect of MT on apoptosis of U87 and U251 glioma cells. (a) Cells treated with different concentrations of MT for 48 h. Cell apoptosis detected by flow cytometry. The lower left represents normal cells, the lower right represents early apoptotic cells, the upper right represents late apoptotic cells, and the upper left represents dead or late apoptotic cells. (b) Analysis of the altered rate of early apoptosis, late apoptosis, and total apoptosis after MT treatment. Results are expressed as , , vs. the control group.

3.3. Transcriptome Sequencing

To investigate the possible mechanisms of MT in inhibiting U87 glioma cells, we analyzed the transcription levels of associated mRNAs. Total RNAs were isolated from the control and MT-treated groups, and mRNA sequencing was performed. A total of 48 DEGs were obtained with a and . Among these, 37 genes were found to be upregulated, and 11 genes were downregulated in the MT-treated group compared with the control group (Figures 4(a)–4(c)). To determine the potential function of the DEGs, Gene Ontology (GO) enrichment analyses including categories of biological process (BP), cellular component (CC), and molecular function (MF) were performed. The results revealed that the DEGs were significantly enriched in GO0009987 (cellular process), GO0044464 (cell part), and GO0005623 (cell). GO0043226 (organelle), GO0016020 (membrane), GO0032501 (multicellular organismal process), GO0044425 (membrane part), GO0044422 (organelle part), GO0005576 (extracellular region), GO: 0031974 (membrane-enclosed lumen), GO0030054 (cell junction), and GO0044421 (extracellular region part) (Figure 4(d)).

(a)

(b)

(c)

(d)

(e)

Figure 4

MT-induced differentially expressed genes. (a) Number of up- and downregulated genes. (b) Groups of the differentially expressed mRNAs: CON-1, CON-2, and CON-3 (0 nM MT); MT-1, MT-2, and MT-3 (150 nM). (c) Presentation of the volcano plot of DEGs identified from the RNA-seq. Red plots stand for upregulated genes, and blue plots indicate downregulated genes with the following criteria: and . Grey plots indicate nonsignificantly expressed genes. The ordinate shows the -log10 of the adjusted value for each gene, symbolizing the strength of the association. (d) GO analysis of the differentially expressed genes in cellular component, biological process and molecular function. (e) KEGG pathway analyzes the pathways involved in MT regulation.

The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed the top 20 significantly enriched pathways that were directly associated with MT. They involved peroxisome proliferator-activated receptor (PPAR), interleukin 17 (IL-17), adenosine 5-monophosphate- (AMP-) activated protein kinase (AMPK), hematopoietic cell lineage, and cholesterol metabolism (Figure 4(e)).

3.4. PPI Network Analysis and Identification of Key Genes

STRING was used to extract a PPI of the DEGs to evaluate potential interactions of the keg genes following 150 nM MT treatment for 12 h. A total of seven genes were identified from the PPI, including Sterol regulatory element-binding transcription factor 1 (SREBF1), membrane-associated ring-CH-type finger 4 (MARCH4), potassium voltage-gated channel subfamily E regulatory subunit 4 (KCNE4), angiopoietin-like 4 (ANGPTL4), myosin-binding protein H (MYBPH), periplakin (PPL), and complement C5a receptor 2 (C5AR2) (Figure 5). In addition, there are four more important genes protocadherin 18 (PCDH18), inhibitor of DNA-binding 3, HLH protein (ID3), DEPP1 autophagy regulator (DEPP1) and vasorin (VASN) reported in the literature. These genes mainly concentrated in the ion channels: membrane-related signaling pathways.

3.5. Relationship between the Key Genes and Survival Prognosis in Patients with Glioma

RNA-seq data and clinical data of glioma were obtained from TCGA database. Kaplan-Meier analysis was performed on the prognosis of patients with GBM and LGG in TCGA database, and survival curves were drawn. The survival curve revealed that in patients with LGG, there was a significant difference between the high-expression group and the low-expression group of PCDH18 (), PPL (), DEPP1 (), KCNE4 (), MYBPH (), VASN (), ID3 (), and C5AR2 () (Figure 6(a)). Over a 15-year time horizon after diagnosis, the survival rate of patients with low expression of these genes were significantly higher than that of patients with high expression. These results suggested that overexpression of these genes in LGG patients were associated with poor survival. In addition, the survival rate of patients in LGG was significantly different between the low-expression group and the high-expression group of MARCH4 () which suggested that low expression of the gene predicted poor survival (Figure 6(b)). In patients with GBM, B3GALT4 () and MYBPH () were significantly different between the high-expression group and the low-expression group. Overexpression of these genes were associated with poor survival rate in patients with GBM (Figure 6(c)). It is worth noting that MT administration reduced the expression levels of these genes in U87 cells compared with the control cells. These results indicated that the expression level of these genes in glioma is closely associated with the survival of patients, and may play an important role in the development of glioma. MT may promote the apoptosis of tumor cells by regulating the expression level of these genes.

(a)

(b)

(c)

Figure 6

Kaplan-Meier analysis of overall survival (OS) in glioma patients with high and low expression of the key genes (red represents high expression; blue represents low expression). (a) The survival rate of patients with low expression of these genes were significantly higher than that of patients with high expression. Overexpression of PCDH18, PPL, DEPP1, VASN, KCNE4, MYBPH, C5AR2, and ID3 genes in patients with LGG was associated with poor survival. (b) Low expression of MARCH4 predicts poor survival in LGG patients. (c) Overexpression of B3GALT4 and MYBPH genes in GBM patients was associated with poor survival.

3.6. Expression of Survival Key Genes in Grade 2 (G2) and Grade 3 (G3) Glioma Patients

Of the 508 patients with glioma, the expression of PCDH18 (), DEPP1 (), VASN (), KCNE4 (), MYBPH () and MARCH4 () genes was significantly different between G2 and G3. Except that the expression of MARCH4 was decreased in G3, the expression of other genes was increased in G3 (Figures 7(a) and 7(b)).

(a)

(b)

3.7. Univariate Independent Prognostic Analysis

Univariate and multivariate Cox regression analyses were performed on 508 glioma patients to assess the influence of genes associated with glioma survival and other clinicopathological factors on survival status.

Univariate independent prognostic analysis was conducted to compare clinical traits and risk values individually with survival time and survival status to analyze whether these factors were associated with survival. Risk score indicates that the correlation between the factor and survival is significant. Univariate independent prognostic analysis showed that age (), tumor grade (), PCDH18 (), PPL (), DEPP1 (), VASN (), KCNE4 (), MYBPH (), C5AR2 (), and MARCH4 () gene expression levels were significantly different (Table 3). It suggested that age, grade, and the level of PCDH18, PPL, DEPP1, VASN, KCNE4, MYBPH, C5AR2, and MARCH4 genes were important predictors of survival.

3.8. Multivariate Factor Analysis

Multivariate factor analysis is to compare clinical characteristics and risk values together with survival time and survival status, and it takes into account the influence of factors. Risk score indicated that the risk value can be used as an independent prognostic factor [39]. Multivariate Cox regression analysis showed that the expression of PCDH18 (), PPL (), DEPP1 (), VASN (), KCNE4 (), MYBPH (), C5AR2 (), and MARCH4 () was an independent predictor of glioma which is most closely related to the survival of patients with glioma, independent of other variables, and has significant prognostic value (Figure 8).

3.9. Quantitative Real-Time PCR Analyses of the Key Genes Related to MT Treatment

In order to confirm the regulatory effect of MT on the eight key genes related to the prognosis of patients with glioma, the eight key genes between MT-treated (150 nM treatment for 12 h) and MT-untreated groups were selected and verified using qPCR. The gene expression levels, using the method and normalization to GAPDH as a reference, were determined. The results showed that PCDH18, PPL, DEPP1, VASN, KCNE4, MYBPH, and C5AR2 genes were found to be downregulated (Figure 9(a)), and MARCH4 gene was upregulated compared to the control group (Figure 9(b)). This result suggested that MT can regulate the expression of these genes in glioma U87 cells. The expression profiles of these genes was consistent with the RNA sequencing data.

(a)

(b)

4. Discussion

In our study, we identified 48 genes that differed between MT treatment group and control group by transcriptome sequencing. GO analysis showed that these genes are mostly enriched in membrane-related signaling pathways such as membrane part, membrane, and membrane-enclosed lumen. The results suggest that MT mainly changed the cell morphology and destroyed cell membrane, thus inducing U87 cell death. KEGG enrichment analysis demonstrated that the signaling pathways associated with these DEGs were pathways in PPAR, IL-17, and AMPK signaling pathways. PPAR plays an important role in the occurrence and development of glioma, and PPAR agonists might represent novel adjuvant therapeutic agents for the treatment of gliomas [40]. In glioma samples, PPAR hyperactivation is associated with immunosuppression through increased regulatory T cell expression [41], and MT has been reported to be involved in T cell response and has immunomodulatory effects [42]. Whether MT could act as a PPAR agonist to promote tumor cell apoptosis by regulating the immune system will be explored in our future research.

Studies have shown that ion channel activity changes during tumorigenesis and development, and the abnormal expression and activity changes of ion channels are closely related to tumor cell proliferation and apoptosis [43]. In particular, the migration and invasion of glioma cells are promoted by ion channels and transporters [44]. Voltage-gated potassium channels are a large family of ion channels that are important transmembrane channels in neurons, regulating cell membrane potential and cell proliferation. Some studies have found that the occurrence and development of a variety of tumors are closely related to the dysregulation of voltage-gated potassium channel, such as glioma [45] and lung cancer [46]. The effects of voltage-gated potassium channels on proliferation and apoptosis of glioma cells have been widely studied, which provides a new therapeutic target for the prevention and treatment of glioma. In this study, KEGG enrichment analysis showed that voltage-gated potassium channels were abnormally expressed in U87 cells. It is noteworthy that MT may be involved in the apoptotic process of tumor cells by regulating voltage-gated potassium channels to change membrane function and influence membrane potential.

Through TCGA database, we identified that the eight key genes (PCDH18, PPL, DEPP1, VASN, KCNE4, MYBPH, C5AR2, and MARCH4) were associated with the survival of patients with glioma. Overexpression of PCDH18, PPL, DEPP1, VASN, KCNE4, MYBPH, and C5AR2 in glioma and low expression of MARCH4 were associated with poor survival. Interestingly, MT can downregulate the overexpression of PCDH18, PPL, DEPP1, VASN, KCNE4, MYBPH, and C5AR2 genes and upregulate the expression of MARCH4 gene in U87 cells. Independent prognostic analysis showed that eight genes are strongly associated with patient outcomes. It can be used as an independent prognostic factor in patients with glioma.

KCNE4 is expressed abundantly in the brain [47]. It can control the firing rate of neurons and synaptic transmission [48]. KCNE4 is one of the genes that modify ion channels in glioblastoma [49]. KCNE subunits are present in the immune system and may play a role in the immune system via associations with leukocyte K⁺ channels, specifically KCNE4, as novel targets for immunomodulation [50]. Our study showed that KCNE4 is closely associated with the prognosis of cancer patients. MT can inhibit the expression of KCNE4 gene on the gated potassium channel in U87 cells. MT may act as an inhibitor of potassium channel to inhibit the growth and proliferation of U87 cells and induce apoptosis of U87 cells. As a member of the ID protein family, ID3 plays an important role in cell proliferation, differentiation, and senescence and is involved in metastasis and angiogenesis of malignant tumors [51]. ID3 expression is increased in advanced glioma [52], and high ID3 expression had a shorter overall survival time in glioma patients [53]. Studies have shown that VASN protein has a cancer-promoting effect [54] and may be closely related to the occurrence and development of liver cancer [55]. Some studies found that glioma patients with high VASN expression had a shorter overall survival time. VASN stimulates tumor progression and angiogenesis in glioma and represents a novel therapeutic target for glioma [56]. Our findings were inconsistent with these studies. Importantly, we found that MT can downregulate the expression of VASN in U87 cells. MARCH4 plays an important regulatory role in the immune system, and its mechanism of action in glioma warrants further studied.

Although MT has the potential to inhibit tumor cell growth, maybe it is also cytotoxic to normal cells, limiting its clinical use [57, 58]. It has been suggested that this can be solved by targeted delivery of melittin nanoparticles [59]. Someone developed transdermal delivery preparations for bee venom, such as bee venom plastics with transdermal administration; it was found to be safe in the range of experimental doses and showed no acute or long-term toxicity [60]. In addition, the bee-sting therapy was used in clinical, and it was found that there were no abnormalities in blood, urine, stool routine, electrocardiogram, liver, and kidney function within the safe dose range [61]. Bee products and cupping are also commonly used in combination with bee-sting therapy, which effectively improve the quality of life and survival time of patients [62]. It was also reported that the toxicity of the peptide to normal cells was less [25], which may be related to the concentration of the peptide. Therefore, it is possible to reduce the toxicity of MT on normal tissues and cells through dose control or the use of drug delivery vectors.

5. Conclusion

It is worth noting that this study has certain limitations. For instance, no in vivo validation was performed, and an in-depth mechanistic exploration is needed. However, the current study provides important clues regarding that MT promotes apoptosis of U87 cells by regulating key genes that are mainly enriched in membrane potential, PPAR signaling pathway, and voltage-gated potassium channel. These genes are closely associated with the prognosis of patients with glioma. MT is a potential drug for the treatment of glioma.

Abbreviations

MT:	Melittin
WHO:	World Health Organization
CNS:	Central nervous system
LGG:	Low-grade glioma
HGG:	High-grade glioma
GBM:	Glioblastoma multiforme
TT Fields:	Tumor-treating fields
RT:	Radiotherapy
DEGs:	Differentially expressed genes
RNA-seq:	Transcriptome sequencing
PPI:	Protein-protein interaction
TCGA:	The Cancer Genome Atlas
TMZ:	Temozolomide
CCK-8:	Cell Counting Kit-8
GO:	Gene Ontology
KEGG:	Kyoto Encyclopedia of Genes and Genomes
SREBF1:	Sterol regulatory element-binding transcription factor 1
KCNE4:	Potassium voltage-gated channel subfamily E regulatory subunit 4
PCDH18:	Protocadherin 18
PPL:	Periplakin
DEPP1:	DEPP1 autophagy regulator
MYBPH:	Myosin-binding protein H
C5AR2:	C5a anaphylatoxin chemotactic receptor 2
ID3:	DNA-binding protein inhibitor ID-3
VASN:	Vasorin
ANGPTL4:	Angiopoietin-like 4
B3GALT4:	Beta-1,3-galactosyltransferase 4
MARCH4:	Membrane-associated ring-CH-type finger 4.

Data Availability

RNA-seq data and clinical data for glioma patients can be found in TCGA database (https://www.cancer.gov/tcga); the rest can be obtained from the corresponding author on reasonable request.

Informed consent was obtained from all individual participants included in this study.

Conflicts of Interest

The authors declare no conflict of interests.

Authors’ Contributions

CX conceived the concept and wrote the original draft. QR, XG, and SX conducted the QPCR experiment. CX, TT, and RL conducted the cell experiment. RL, DC, and JW reviewed and edited the manuscript. All authors read and approved the final manuscript. Jie Wu, Diling Chen, and Ran Li contributed equally to this work.

Acknowledgments

We thank Professor Guan Yi for providing U87 cells for free. This work was supported by the Research Foundation of Natural Science Foundation of Hunan Province (nos. 2022JJ40460, 2022JJ30041, and 2022JJ30579), the Hunan Health Commission Project (nos. 202102081835 and 202103101455), the Research Foundation of Hunan Administration of Traditional Chinese Medicine (nos. 2021091 and 2019101), Guangdong Administration of Traditional Chinese Medicine (nos. 20201175 and 20202090), the Clinical Innovation Leading Scientific Research Projects of Hunan Science and Technology Department (no. 2020SK52603), Hunan Provincial Clinical Medical Technology Demonstration Base Project (no. 2020SK4025), and the Scientific Research Project of Chinse Ethnic Medicine Society (no. 2019KYXM-Z1130-115).

References

C. L. Gladson, R. A. Prayson, and W. M. Liu, “The pathobiology of glioma tumors,” Annual Review of Pathology, vol. 5, no. 1, pp. 33–50, 2010.
View at: Publisher Site | Google Scholar
P. D. Delgado-Lopez and E. M. Corrales-Garcia, “Survival in glioblastoma: a review on the impact of treatment modalities,” Clinical & Translational Oncology, vol. 18, no. 11, pp. 1062–1071, 2016.
View at: Publisher Site | Google Scholar
A. D. Joshi, D. W. Parsons, V. E. Velculescu, and G. J. Riggins, “Sodium ion channel mutations in glioblastoma patients correlate with shorter survival,” Molecular Cancer, vol. 10, no. 1, 2011.
View at: Publisher Site | Google Scholar
P. Śledzińska, G. Bebyn Marek, J. Furtak, J. Kowalewski, and M. A. Lewandowska, “Prognostic and predictive biomarkers in gliomas,” International Journal of Molecular Sciences, vol. 22, no. 19, article 10373, 2021.
View at: Publisher Site | Google Scholar
N. Sanai, M. Y. Polley, M. W. Mcdermott, A. T. Parsa, and M. S. Berger, “An extent of resection threshold for newly diagnosed glioblastomas,” Journal of Neurosurgery, vol. 115, no. 1, pp. 3–8, 2011.
View at: Publisher Site | Google Scholar
R. Stupp, S. Taillibert, A. Kanner et al., “Effect of tumor-treating fields plus maintenance temozolomide vs maintenance temozolomide alone on survival in patients with glioblastoma: a randomized clinical trial,” Journal of the American Medical Association, vol. 318, no. 23, pp. 2306–2316, 2017.
View at: Publisher Site | Google Scholar
S. R. y. Cajal, M. Sesé, C. Capdevila et al., “Clinical implications of intratumor heterogeneity: challenges and opportunities,” Journal of Molecular Medicine, vol. 98, no. 2, pp. 161–177, 2020.
View at: Publisher Site | Google Scholar
C. E. Brown, D. Alizadeh, R. Starr et al., “Regression of glioblastoma after chimeric antigen receptor T-cell therapy,” The New England Journal of Medicine, vol. 375, no. 26, pp. 2561–2569, 2016.
View at: Publisher Site | Google Scholar
I. Valduvieco, E. Verger, J. Bruna et al., “Impact of radiotherapy delay on survival in glioblastoma,” Clinical & Translational Oncology, vol. 15, no. 4, pp. 278–282, 2013.
View at: Publisher Site | Google Scholar
M. Esteller, J. Garcia-Foncillas, E. Andion et al., “Inactivation of the DNA-repair gene MGMT and the clinical response of gliomas to alkylating agents,” New England Journal of Medicine, vol. 343, no. 19, pp. 1350–1354, 2000.
View at: Publisher Site | Google Scholar
Z. Y. Yu, G. F. Xie, G. T. Zhou et al., “NVP-BEZ235, a novel dual PI3K-mTOR inhibitor displays anti-glioma activity and reduces chemoresistance to temozolomide in human glioma cells,” Cancer Letters, vol. 367, no. 1, pp. 58–68, 2015.
View at: Publisher Site | Google Scholar
T. C. Terwilliger and D. Eisenberg, “The structure of melittin. II. Interpretation of the structure,” The Journal of Biological Chemistry, vol. 257, no. 11, pp. 6016–6022, 1982.
View at: Publisher Site | Google Scholar
S. Somwongin, P. Chantawannakul, and W. Chaiyana, “Antioxidant activity and irritation property of venoms from Apis species,” Toxicon, vol. 145, pp. 32–39, 2018.
View at: Publisher Site | Google Scholar
M. Hamed and M. Mojtaba, “Anti-fungal properties and mechanisms of melittin,” Applied Microbiology and Biotechnology, vol. 104, no. 15, pp. 6513–6526, 2020.
View at: Publisher Site | Google Scholar
G. Lee and H. Bae, “Anti-inflammatory applications of melittin, a major component of bee venom: detailed mechanism of action and adverse effects,” Molecules, vol. 21, no. 5, p. 616, 2016.
View at: Publisher Site | Google Scholar
J. H. Park, Y. S. Kum, T. I. Lee et al., “Melittin attenuates liver injury in thioacetamide-treated mice through modulating inflammation and fibrogenesis,” Experimental Biology and Medicine, vol. 236, no. 11, pp. 1306–1313, 2011.
View at: Publisher Site | Google Scholar
S. D. He, N. Tan, C. X. Sun et al., “Treatment with melittin induces apoptosis and autophagy of fibroblastlike synoviocytes in patients with rheumatoid arthritis,” Current Pharmaceutical Biotechnology, vol. 21, no. 8, pp. 734–740, 2020.
View at: Publisher Site | Google Scholar
N. Tan, S. D. He, G. Li, P. C. Yang, Z. G. Liu, and D. Y. Li, “Effects of melittin on Th17/Treg balance and inflammation in rat model of collagen-induced arthritis,” Chinese Journal of Comparative Medicine, vol. 29, no. 2, pp. 1–6, 2019.
View at: Google Scholar
S. K. Kim, K. Y. Park, W. C. Yoon, S. H. Park, D. H. Yoo, and J. Y. Choe, “Melittin enhances apoptosis through suppression of IL-6/sIL-6R complex-induced NF-κB and STAT3 activation and Bcl-2 expression for human fibroblast-like synoviocytes in rheumatoid arthritis,” Joint, Bone, Spine, vol. 78, no. 5, pp. 471–477, 2011.
View at: Publisher Site | Google Scholar
Z. C. Wang, Q. Chen, X. Liu, X. H. Xie, X. F. Jiang, and W. H. Zhu, “Research progress on anti-tumor mechanism of melittin,” Journal of Jilin Medical College, vol. 39, no. 6, pp. 462–464, 2018.
View at: Google Scholar
I. Rady, I. A. Siddiqui, M. Rady, and H. Mukhtar, “Melittin, a major peptide component of bee venom, and its conjugates in cancer therapy,” Cancer Letters, vol. 402, pp. 16–31, 2017.
View at: Publisher Site | Google Scholar
Q. W. Chen, W. F. Lin, Z. F. Yin et al., “Melittin Inhibits Hypoxia-Induced Vasculogenic Mimicry Formation and Epithelial-Mesenchymal Transition through Suppression of HIF-1/Akt Pathway in Liver Cancer,” Evidence-based Complementary and Alternative Medicine, vol. 2019, Article ID 9602935, 10 pages, 2019.
View at: Publisher Site | Google Scholar
K. E. Choi, C. J. Hwang, S. M. Gu et al., “Cancer cell growth inhibitory effect of bee venom via increase of death receptor 3 expression and inactivation of NF-kappa B in NSCLC cells,” Toxins, vol. 6, no. 8, pp. 2210–2228, 2014.
View at: Publisher Site | Google Scholar
R. Z. Yu, M. Wang, M. H. Wang, and L. Han, “Melittin suppresses growth and induces apoptosis of non-small-cell lung cancer cells via down-regulation of TGF-β-mediated ERK signal pathway,” Brazilian Journal of Medical and Biological Research, vol. 54, no. 2, pp. 1414–1431, 2021.
View at: Publisher Site | Google Scholar
C. Duffy, A. Sorolla, E. Wang et al., “Honeybee venom and melittin suppress growth factor receptor activation in HER2-enriched and triple-negative breast cancer,” npj Precision Oncology, vol. 4, no. 1, p. 24, 2020.
View at: Publisher Site | Google Scholar
H. N. Lim, S. B. Baek, and H. J. Jung, “Bee venom and its peptide component melittin suppress growth and migration of melanoma cells via inhibition of PI3K/AKT/mTOR and MAPK pathways,” Molecules, vol. 24, no. 5, p. 929, 2019.
View at: Publisher Site | Google Scholar
M. Jo, M. H. Park, P. S. Kollipara et al., “Anti-cancer effect of bee venom toxin and melittin in ovarian cancer cells through induction of death receptors and inhibition of JAK2/STAT3 pathway,” Toxicology and Applied Pharmacology, vol. 258, no. 1, pp. 72–81, 2012.
View at: Publisher Site | Google Scholar
Z. L. Yang, Y. Q. Ke, R. X. Xu, and P. Peng, “Melittin inhibits proliferation and induces apoptosis of malignant human glioma cells,” Journal of Southern Medical University, vol. 27, no. 11, pp. 1775–1777, 2007.
View at: Google Scholar
A. A. Lebel, M. V. Kisembo, M. N. Soucy, M. P. A. Hébert, P. J. Morin, and L. H. Boudreau, “Molecular characterization of the anticancer properties associated with bee venom and its components in glioblastoma multiforme,” Chemico-Biological Interactions, vol. 347, article 109622, 2021.
View at: Publisher Site | Google Scholar
G. Gajski, T. Čimbora-Zovko, S. Rak, M. Osmak, and V. Garaj-Vrhovac, “Antitumour action on human glioblastoma A1235 cells through cooperation of bee venom and cisplatin,” Cytotechnology, vol. 68, no. 4, pp. 1197–1205, 2016.
View at: Publisher Site | Google Scholar
X. J. Yao, Y. H. Wang, G. C. Deng, M. M. Zeng, and B. X. Guang, “Effect of recombinant melittin on inhibition of glioma C6 cell growth and the expression of VEGF, Flt-1, MMP-9,” Chinese Journal of Experimental Traditional Medical Formulae, vol. 19, no. 24, pp. 267–269, 2013.
View at: Google Scholar
M. Sisakht, B. Mashkani, A. Bazi et al., “Bee venom induces apoptosis and suppresses matrix metaloprotease-2 expression in human glioblastoma cells,” Brazilian Journal of Pharmacognosy, vol. 27, no. 3, pp. 324–328, 2017.
View at: Publisher Site | Google Scholar
X. J. Yao, Y. H. Wang, W. Z. Hou, G. C. Deng, M. M. Zeng, and B. X. Guang, “Effect of recombinant melittin on inhibition of glioma SHG44 cell growth and the expression of STAT3, VEGF,” Chinese Journal of Gerontology, vol. 33, pp. 3372–3374, 2013.
View at: Google Scholar
R. Li, C. Q. Xu, J. X. Shen et al., “4-Methoxydalbergione is a potent inhibitor of human astroglioma U87 cells in vitro and in vivo,” Acta Pharmacologica Sinica, vol. 42, no. 9, pp. 1507–1515, 2021.
View at: Publisher Site | Google Scholar
S. K. Miryala, A. Anbarasu, and S. Ramaiah, “Discerning molecular interactions: a comprehensive review on biomolecular interaction databases and network analysis tools,” Gene, vol. 642, pp. 84–94, 2018.
View at: Publisher Site | Google Scholar
B. George, S. Seals, and I. Aban, “Survival analysis and regression models,” Journal of Nuclear Cardiology, vol. 21, no. 4, pp. 686–694, 2014.
View at: Publisher Site | Google Scholar
G. M. Callagy, M. J. Webber, P. D. Pharoah, and C. Caldas, “Meta-analysis confirms BCL2 is an independent prognostic marker in breast cancer,” BMC Cancer, vol. 8, no. 1, p. 153, 2008.
View at: Publisher Site | Google Scholar
M. S. Schröder, A. C. Culhane, J. Quackenbush, and B. Haibe-Kains, “Survcomp: an R/bioconductor package for performance assessment and comparison of survival models,” Bioinformatics, vol. 27, no. 22, pp. 3206–3208, 2011.
View at: Publisher Site | Google Scholar
H. H. Wang, D. H. Guo, J. J. Li, B. L. Wei, and H. Zheng, “Increased expression of osteopontin indicates poor prognosis in hepatocellular carcinoma,” International Journal of Clinical and Experimental Pathology, vol. 11, no. 12, pp. 5916–5922, 2018.
View at: Google Scholar
N. Khoo, S. C. Hebbar, W. L. Zhao, S. A. Moore, F. E. Domann, and M. E. Robbins, “Differential activation of catalase expression and activity by PPAR agonists: implications for astrocyte protection in anti-glioma therapy,” Redox Biology, vol. 1, no. 1, pp. 70–79, 2013.
View at: Publisher Site | Google Scholar
W. H. Chang and A. G. Lai, “The pan-cancer mutational landscape of the PPAR pathway reveals universal patterns of dysregulated metabolism and interactions with tumor immunity and hypoxia,” Annals of the New York Academy of Sciences, vol. 1448, no. 1, pp. 65–82, 2019.
View at: Publisher Site | Google Scholar
W. Y. Li, Bee-Sting Therapy, People’s Medical Publishing House, Beijing, 1th edition, 2003.
V. R. Rao, M. Perez-Neut, S. Kaja, and S. Gentile, “Voltage-gated ion channels in cancer cell proliferation,” Cancers, vol. 7, no. 2, pp. 849–875, 2015.
View at: Publisher Site | Google Scholar
V. A. Cuddapah and H. Sontheimer, “Ion channels and transporters [corrected] in cancer. 2. Ion channels and the control of cancer cell migration,” American Journal of Physiology. Cell Physiology, vol. 301, no. 3, pp. C541–C549, 2011.
View at: Publisher Site | Google Scholar
K. Preußat, C. Beetz, M. Schrey et al., “Expression of voltage-gated potassium channels Kv1.3 and Kv1.5 in human gliomas,” Neuroscience Letters, vol. 346, no. 1-2, pp. 33–36, 2003.
View at: Publisher Site | Google Scholar
J. H. Lee, J. W. Park, J. K. Byun et al., “Silencing of voltage-gated potassium channel KV9.3 inhibits proliferation in human colon and lung carcinoma cells,” Oncotarget, vol. 6, no. 10, pp. 8132–8143, 2015.
View at: Publisher Site | Google Scholar
L. J. Manderfield and A. L. George, “KCNE4 can co-associate with the IKs (KCNQ1–KCNE1) channel complex,” FEBS Journal, vol. 275, no. 6, pp. 1336–1349, 2008.
View at: Publisher Site | Google Scholar
M. Grunnet, H. B. Rasmussen, A. Hay-Schmidt et al., “KCNE4 is an inhibitory subunit to Kv1.1 and Kv1.3 potassium channels,” Biophysical Journal, vol. 85, no. 3, pp. 1525–1537, 2003.
View at: Publisher Site | Google Scholar
A. Biasiotta, D. D. Arcangelo, F. Passarelli, E. M. Nicodemi, and A. Facchiano, “Ion channels expression and function are strongly modified in solid tumors and vascular malformations,” Journal of Translational Medicine, vol. 14, no. 1, p. 285, 2016.
View at: Publisher Site | Google Scholar
L. Solé, M. Roura-Ferrer, M. Pérez-Verdaguer et al., “KCNE4 suppresses Kv1.3 currents by modulating trafficking, surface expression and channel gating,” Journal of Cell Science, vol. 122, no. 20, pp. 3738–3748, 2009.
View at: Publisher Site | Google Scholar
H. M. Jeon, X. Jin, J. S. Lee et al., “Inhibitor of differentiation 4 drives brain tumor-initiating cell genesis through cyclin E and notch signaling,” Genes & Development, vol. 22, no. 15, pp. 2028–2033, 2008.
View at: Publisher Site | Google Scholar
D. A. Vandeputte, D. Troast, S. Leenstra et al., “Expression and distribution of id helix-loop-helix proteins in human astrocytic tumors,” Glia, vol. 38, no. 4, pp. 329–338, 2002.
View at: Publisher Site | Google Scholar
R. K. Mai, H. T. Yu, Z. F. Zhu, J. Q. Feng, and Y. Zheng, “The relationship among the expression level of inhibitor of differentiation/DNA binding-factor (ID) gene, pathologic grade and survival in patients with gliomas,” Anatomy Research International, vol. 41, no. 4, pp. 295–298, 2019.
View at: Google Scholar
W. Wang, “Study on VASN protein as serum marker and therapeutic target of liver cancer,” Academy of Military Medical Sciences of the People’s Liberation Army, 2012.
View at: Google Scholar
A. X. Huang, “Preliminary study on the biological function and regulatory mechanism of VASN protein,” Academy of Military Medical Sciences of the Chinese People’s Liberation Army, 2015.
View at: Google Scholar
W. Y. Liang, B. Y. Guo, J. C. Ye et al., “Vasorin stimulates malignant progression and angiogenesis in glioma,” Cancer Science, vol. 110, no. 8, pp. 2558–2572, 2019.
View at: Publisher Site | Google Scholar
G. Gajski and V. Garaj-Vrhovac, “Melittin: a lytic peptide with anticancer properties,” Environmental Toxicology and Pharmacology, vol. 36, no. 2, pp. 697–705, 2013.
View at: Publisher Site | Google Scholar
G. Gajski, A. M. Domijan, B. Žegura et al., “Melittin induced cytogenetic damage, oxidative stress and changes in gene expression in human peripheral blood lymphocytes,” Toxicon, vol. 110, pp. 56–67, 2016.
View at: Publisher Site | Google Scholar
H. Pan, N. R. Soman, P. H. Schlesinger, G. M. Lanza, and S. A. Wickline, “Cytolytic peptide nanoparticles (‘NanoBees’) for cancer therapy,” Wiley Interdisciplinary Reviews. Nanomedicine and Nanobiotechnology, vol. 3, no. 3, pp. 318–327, 2011.
View at: Publisher Site | Google Scholar
Y. Gao, F. Zhu, Z. Yang et al., “Safety evaluation of bee venom plastics with transdermal administration,” China pharmacy, vol. 30, no. 16, pp. 2181–2186, 2019.
View at: Google Scholar
S. Y. Chen, The Study the Clinical Efficacy and Safety of Bee Acupuncture for Rheumatoid Arthritis, Guangzhou University of Chinese Medicine, 2014.
X. Zhang, C. Z. Tang, W. Y. Li, and R. Y. Huang, “Discussion on the treatment of the advanced lung cancer patient with bee-sting therapy,” Asia-Pacific Traditional Medicine, vol. 14, no. 2, pp. 140–143, 2018.
View at: Google Scholar

Copyright

Copyright © 2022 Ran Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

474

Downloads

487

Citations