Abstract
Objective: To investigate the prognostic value of arachidonate lipoxygenases 5 (ALOX5) expression and methylation, and explore the immune functions of arachidonate lipoxygenases 5 expression in low-grade glioma (LGG).
Materials and Methods: Using efficient bioinformatics approaches, the differential expression of arachidonate lipoxygenases 5 and the association of its expression with clinicopathological characteristics were evaluated. Then, we analyzed the prognostic significance of arachidonate lipoxygenases 5 expression and its methylation level followed by immune cell infiltration analysis. The functional enrichment analysis was conducted to determine the possible regulatory pathways of arachidonate lipoxygenases 5 in low-grade glioma. Finally, the drug sensitivity analysis was performed to explore the correlation between arachidonate lipoxygenases 5 expression and chemotherapeutic drugs.
Results: arachidonate lipoxygenases 5 mRNA expression was increased in low-grade glioma and its expression had a notable relation with age and subtype (p < 0.05). The elevated mRNA level of arachidonate lipoxygenases 5 could independently predict the disease-specific survival (DSS), overall survival (OS), and progression-free interval (PFI) (p < 0.05). Besides, arachidonate lipoxygenases 5 expression was negatively correlated with its methylation level and the arachidonate lipoxygenases 5 hypomethylation led to a worse prognosis (p < 0.05). The arachidonate lipoxygenases 5 expression also showed a positive connection with immune cells, while low-grade glioma patients with higher immune cell infiltration had poor survival probability (p < 0.05). Further, arachidonate lipoxygenases 5 might be involved in immune- and inflammation-related pathways. Importantly, arachidonate lipoxygenases 5 expression was negatively related to drug sensitivity.
Conclusion: arachidonate lipoxygenases 5 might be a promising biomarker, and it probably occupies a vital role in immune cell infiltration in low-grade glioma.
Introduction
Glioma is the most common primary malignant brain tumor characterized by a high disability rate, high recurrence rate, and high mortality rate (). The molecular parameters and histology were incorporated into the 2016 version of the WHO classification of central nervous system tumors (). As of 2016, gliomas are divided into circumscribed gliomas (WHO grade I) and diffusely infiltrating gliomas (WHO grades II-IV) based on their pattern of growth and the IDH mutation status (). Surgery, neuroimaging, chemotherapy, radiation therapy, and neuropathology are available approaches for the management of gliomas (). Circumscribed gliomas are benign and curable by complete surgical resection, while diffuse gliomas are almost never cured by resection alone (). Histologically, grade II (low-grade) diffuse gliomas display nuclear atypia, grade III (anaplastic) show increased mitotic activity, and grade IV (glioblastomas) display additional microvascular proliferation, necrosis, or both (). Diffuse low-grade gliomas (LGG) include grade II astrocytomas, oligodendrogliomas, and oligoastrocytomas (). After comprehensive treatments including surgical resection, chemotherapy, and radiotherapy, LGG patients have a better prognosis with a median survival of more than 10 years (). However, they have inevitably suffered from recurrence and malignant progression due to the highly invasive nature of LGG (). Currently, biomarkers such as O6-methylguanine DNA methyltransferase (MGMT) and isocitrate dehydrogenase 1 (IDH1) have become vital markers of LGG clinical behavior and are closely connected with patient prognosis (). Further development of prognostic indicators for LGG is urgently needed to gain additional insights to provide additional potential therapeutic targets.
Lipid peroxidation, which preferentially oxidizes polyunsaturated fatty acids, has been involved in the etiology of diverse pathological conditions such as diabetes and cancers (; ; ). Arachidonic acid (AA) is one of the most important polyunsaturated fatty acids in mammalian cells, which is required for normal cellular membrane fluidity and is the direct precursor of various bioactive mediators including endocannabinoids and leukotrienes (LTs) associated with inflammation (). AA is mainly metabolized by arachidonate lipoxygenases (ALOXs) which have six functional subtypes (). ALOX5 is a non-heme iron-containing dioxygenase that encodes lipoxygenase and metabolizes AA into hydroperoxyderivatives (5-HPETE) and further into 5-hydroxyeicosatetraenoic acid (5-HETE), conferring growth, invasion, and chemopreventive advantage in cancer cells (; ). Kumar et al. demonstrated that the serum level of ALOX5 was elevated in breast cancer, and serves as a promising therapeutic target (). In addition, ALOX5 is upregulated in gastrointestinal cancer and its overexpression is related to the cancer progression (). However, the comprehensive role of ALOX5 in LGG remains unclarified.
In this study, we analyzed the expression level of ALOX5 in LGG and normal tissues followed by evaluating its prognostic value in LGG. Besides, the survival analysis of ALOX5 methylation and immune cell infiltration was conducted. We also performed functional enrichment analyses and revealed the molecular characteristic of ALOX5 in LGG. Finally, the drug sensitivity analysis was performed to explore the correlation between ALOX5 expression and sensitivity of chemotherapeutic drugs.
Materials and methods
Arachidonate lipoxygenases 5 gene expression analysis
The GEPIA database (http://gepia.cancer-pku.cn/) integrates TCGA big data for cancerous tissue and GTEx big data for normal tissue. This database was used to analyze the ALOX5 gene expression across all tumor samples and normal tissues, and to examine the expression of the ALOX5 gene in LGG and normal groups. Next, the immunohistochemistry (IHC) analysis was conducted to assess the ALOX5 protein level in LGG and normal tissues following the manufacturer’s protocol.
Then, normalized RNA-seq data and clinical information for LGG samples from the TCGA cohort were obtained from the UCSC Xena database (https://xenabrowser.net/), which has recomputed all raw expression data from TCGA. Inclusion criteria: 1) patients in WHO grades II gliomas; 2) with matched survival and expression data. Finally, 260 samples were enrolled in the study. The patients were divided into high- and low- ALOX5 expression groups with the median expression of ALOX5, and the association of ALOX5 with clinicopathological characteristics was analyzed using the chi-square test. Besides, we assessed the correlation between ALOX5 as a continuous variable and clinical parameters by t-test or ANOVA.
Arachidonate lipoxygenases 5 gene survival analysis
Kaplan-Meier plotter curves were drawn to examine the correlation between ALOX5 mRNA expression and DFI, DSS, OS, and PFI of patients using the TCGA-LGG data. The log-rank test was used to assess the survival differences between high-and low- ALOX5 expression groups. Next, ROC analysis was conducted to evaluate the value of ALOX5 in predicting the survival status of LGG patients by using R package “timeROC” (), and the area under the curve (AUC) was calculated using the R package “pROC” (). For validation, the glioma mRNA-seq dataset “mRNAseq_693” was acquired from the CGGA database (http://www.cgga.org.cn/). The Kaplan-Meier plotter method and ROC analysis were employed to verify the prognostic role of ALOX5 in LGG patients.
Following this, Cox regression analysis was performed to assess the independent prognostic factors for LGG patients using the TCGA and CGGA data. Nomogram models were constructed based on the multivariate Cox regression analysis results and calibration curves were drawn.
Arachidonate lipoxygenases 5 methylation analysis
Gene Set Cancer Analysis (GSCA) (http://bioinfo.life.hust.edu.cn/GSCA/#/) is an integrated platform for genomic, pharmacogenomic, and immunogenomic gene set cancer analysis. We adopted this database to analyze the relationship between ALOX5 mRNA expression and its methylation level.
MEXPRESS (https://mexpress.be/) is a data visualization tool designed for the easy visualization of TCGA expression, DNA methylation, and clinical data, as well as the relationships between them. We downloaded the DNA methylation data of ALOX5 in LGG for analyzing its prognostic value.
Immune cell infiltration analysis
To evaluate the correlation between ALOX5 expression and tumor immune microenvironment in LGG, the ESTIMATE algorithm was adopted to analyze the Immune score, Stromal score, and Estimate score. The effect of Immune score, Stromal score, and Estimate score on patient prognosis was explored using the Kaplan-Meier plotter method. Additionally, the association of ALOX5 expression with B cell, CD8 + T cell, CD4 + T cell, macrophage, neutrophil, and dendritic cell was evaluated using the TIMER algorithm. Similarly, the survival analyses of these immune cells were conducted by the Kaplan-Meier plotter method.
Enrichment analysis
We used the R package “limma” to obtain differentially expressed genes (DEGs) between the high and low ALOX5 expression groups according to the median value of ALOX5 using the TCGA data (). The significant DEGs were screened using the threshold of log 2 |foldchange (FC)| >1 and p-value < 0.05, and these DEGs were considered as LGG-related genes. Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of LGG-related genes were performed using the “clusterProfiler” in R ().
Further, the gene set enrichment analysis (GSEA) was performed to illustrate the significant survival difference between the two expression groups. The gene set was permutated 1,000 times and the expression level of F12 was used as a phenotypic label. A nominal p-value <0.05 and an FDR q-value <0.25 were considered to be statistically significant.
Subsequently, the association of ALOX5 expression with immune activation-related genes and immune checkpoint-related genes in LGG was analyzed using the Pearson correlation test.
Drug sensitivity analysis
We investigated the correlation between ALOX5 expression and drug response based on GDSC and CTRP in the GSCA database.
Statistical analysis
All statistical analyses were performed in SPSS software (version 23.0) and packages of R (version 3.6.3). The t-test was used to analyze differences in each two-group comparison, and one-way ANOVA was employed to assess differences among at least three groups. Survival curves were drawn by the Kaplan-Meier method and differences in survival were compared by log-rank tests. p < 0.05 was considered statistically significant.
Results
Arachidonate lipoxygenases 5 gene expression analysis
Firstly, the ALOX5 gene expression in pan-cancer was explored through the GEPIA database as shown in Figure 1A. The ALOX5 mRNA level was significantly higher in the LGG group than that in the normal group (p < 0.01) (Figure 1B). Figure 1C showed the ALOX5 protein level in LGG and normal tissues. Besides, ALOX5 was significantly related to age and subtype (all p < 0.01), but had no significant relation with gender, IDH status, race, and laterality (Table 1) (Figure 1D). Of note, ALOX5 expression was upregulated in LGG patients aged ≤40 years compared with those aged >40 years. Astrocytoma patients seems to exhibit the highest ALOX5 mRNA expression, followed by the mixed type, and oligodendroglioma patients (all p < 0.01) (Figure 1D). These results suggested that high ALOX5 expression might be involved in the occurrence of LGG.
FIGURE 1
TABLE 1
| Characteristics | Low ALOX5 | High ALOX5 | p-value |
|---|---|---|---|
| Age | 0.001 | ||
| ≤40 | 63 (48.5%) | 89 (68.5%) | |
| >40 | 67 (51.5%) | 41 (31.5%) | |
| Gender | 0.619 | ||
| Male | 63 (48.5%) | 59 (45.4%) | |
| Female | 67 (51.5%) | 71 (54.6%) | |
| Subtype | <0.001 | ||
| Astrocytoma | 18 (13.8%) | 47 (36.2%) | |
| Oligodendroglioma | 78 (60.0%) | 39 (30.0) | |
| Mixed | 34 (26.2) | 44 (33.8%) | |
| IDH status | 0.385 | ||
| Wild type | 72 (69.9%) | 82 (75.2%) | |
| Mutation | 31 (30.1%) | 27 (24.8%) | |
| Race | 0.571 | ||
| Asian | 3 (2.4%) | 1 (0.8%) | |
| African-American | 6 (4.7%) | 7 (5.6%) | |
| Caucasian | 118 (92.9%) | 118 (93.7%) | |
| Laterality | 0.380 | ||
| Left | 61 (48.0%) | 68 (53.5%) | |
| Right | 66 (52.0%) | 59 (46.5%) |
Association of ALOX5 expression with clinicopathological characteristics in TCGA-LGG dataset.
High arachidonate lipoxygenases 5 gene expression predicted poor prognosis
To investigate the prognostic significance of ALOX5 gene in LGG, the survival curves were drawn using Kaplan-Meier plotter method based on TCGA-LGG data. ALOX5 expression was not significantly related to the DFI (p > 0.05) (Figure 2A), but patients with higher ALOX5 expression had shorter DSS, OS, and PFI with statistical significance (all p < 0.01) (Figures 2B–D). Although ALOX5 could not predict the DFI status well (AUC: 0.664; p > 0.05), it exhibited acceptable performance in predicting the statuses of DSS, OS, and PFI with AUCs of 0.757, 0.785, and 0.644, respectively (all p < 0.001) (Figures 3A–D).
FIGURE 2
FIGURE 3
For verification, the survival analysis was conducted using CGGA-LGG data. Expectedly, the increased ALOX5 expression contributed to unfavorable OS (p < 0.01) (Figure 4A), and its expression had the ability in distinguishing the OS status of LGG patients (AUC: 0.714; p < 0.001) (Figure 4B).
FIGURE 4
To determine the independent prognostic role of ALOX5 in LGG, univariate and multivariate Cox regression analyses were performed using the TCGA-LGG data. The age, subtype, and ALOX5 were independent predictors for DSS and OS (all p < 0.05) (Table 2; Table 3; Table 4); while literality and ALOX5 were significantly related to PFI in the multivariate analysis (all p < 0.05) (Table 4). These independent factors were enrolled into the construction of nomogram models to predict the probability of DSS, OS, and PFI in LGG. The total number of points was calculated from the sum of the points assigned to each variable in the line graph. The calibration curves demonstrated an agreement between actual observations and predictions of the nomogram in the TCGA-LGG samples (Figures 5A–C). For validation, the Cox regression analysis was conducted using the CGGA-LGG data. The age, IDH status, subtype, and ALOX5 were independent predictors for OS (all p < 0.05) (Table 5). Among these independent factors, ALOX5 contributed the most to predicting the probability of OS in patients with LGG. The calibration curves also showed an agreement between actual observations and predictions of the nomogram in the CGGA-LGG samples (Figure 6). The above findings indicated that ALOX5 gene was significantly related to the LGG progression.
TABLE 2
| Characteristics (References) | Univariate analysis | Multivariate analysis | ||
|---|---|---|---|---|
| HR (95% CI) | p-Value | HR (95% CI) | p-Value | |
| Age | 1.04 (1.01–1.07) | 0.005 | 1.05 (1.02–1.08) | 0.001 |
| Gender (female) | 0.93 (0.49–1.76) | 0.812 | 0.82 (0.38–1.78) | 0.617 |
| IDH status (wild type) | 1.92 (0.84–4.38) | 0.122 | 1.39 (0.55–3.52) | 0.486 |
| Race (Asian) | 1.08 (0.30–3.83) | 0.909 | 1.03 (0.26–4.16) | 0.964 |
| Subtype (astrocytoma) | 1.37 (0.85–2.22) | 0.196 | 1.85 (1.06–3.23) | 0.032 |
| Laterality (left) | 0.56 (0.29–1.01) | 0.087 | 0.49 (0.22–1.10) | 0.085 |
| ALOX5 | 1.53 (1.03–2.28) | 0.034 | 1.59 (1.02–2.47) | 0.039 |
Cox regression analysis of the association of ALOX5 and clinical variables with DSS in the TCGA-LGG dataset.
Abbreviations: DSS, disease-specific survival; HR, hazard ratio; 95% CI, 95% confidence interval.
TABLE 3
| Characteristics (References) | Univariate analysis | Multivariate analysis | ||
|---|---|---|---|---|
| HR (95% CI) | p-Value | HR (95% CI) | p-Value | |
| Age | 1.04 (1.01–1.06) | 0.005 | 1.05 (1.02–1.08) | 0.001 |
| Gender (female) | 0.96 (0.52–1.76) | 0.891 | 0.83 (0.39–1.74) | 0.615 |
| IDH status (wild type) | 1.80 (0.80–4.07) | 0.157 | 1.32 (0.53–3.28) | 0.557 |
| Race (Asian) | 1.11 (0.31–3.95) | 0.871 | 1.08 (0.27–4.26) | 0.914 |
| Subtype (astrocytoma) | 1.37 (0.87–2.15) | 0.174 | 1.90 (1.10–3.25) | 0.021 |
| Laterality (left) | 0.58 (0.31–1.10) | 0.096 | 0.51 (0.23–1.10) | 0.807 |
| ALOX5 | 1.62 (1.12–2.35) | 0.011 | 1.66 (1.09–2.53) | 0.018 |
Cox regression analysis of the association of ALOX5 and clinical variables with OS in the TCGA-LGG dataset.
Abbreviations: OS, overall survival; HR, hazard ratio; 95% CI, 95% confidence interval.
TABLE 4
| Characteristics (References) | Univariate analysis | Multivariate analysis | ||
|---|---|---|---|---|
| HR (95% CI) | p-Value | HR (95% CI) | p-Value | |
| Age | 1.01 (0.99–1.02) | 0.699 | 1.01 (0.98–1.02) | 0.820 |
| Gender (female) | 0.88 (0.58–1.33) | 0.538 | 0.71 (0.46–1.12) | 0.142 |
| IDH status (wild type) | 1.05 (0.60–1.83) | 0.874 | 0.52 (0.32–0.84) | 0.731 |
| Race (Asian) | 0.89 (0.44–1.81) | 0.747 | 0.72 (0.32–1.63) | 0.431 |
| Subtype (astrocytoma) | 0.94 (0.71–1.25) | 0.674 | 0.89 (0.66–1.19) | 0.416 |
| Laterality (left) | 0.50 (0.32–0.78) | 0.002 | 0.52 (0.32–0.84) | 0.007 |
| ALOX5 | 1.55 (1.20–1.99) | 0.001 | 1.55 (1.18–2.05) | 0.002 |
Cox regression analysis of the association of ALOX5 and clinical variables with PFI in the TCGA-LGG dataset.
Abbreviations: PFI, progression-free interval; HR, hazard ratio; 95% CI, 95% confidence interval.
FIGURE 5
TABLE 5
| Characteristics (References) | Univariate analysis | Multivariate analysis | ||
|---|---|---|---|---|
| HR (95% CI) | p-value | HR (95% CI) | p-value | |
| Age | 1.03 (0.99–1.06) | 0.096 | 1.04 (1.01–1.07) | 0.018 |
| Gender (female) | 1.03 (0.57–1.87) | 0.929 | 1.15 (0.61–2.17) | 0.665 |
| IDH status (wild type) | 0.44 (0.23–0.85) | 0.015 | 0.44 (0.21–0.89) | 0.023 |
| Subtype (astrocytoma) | 1.46 (1.02–2.09) | 0.040 | 1.56 (1.10–2.22) | 0.013 |
| ALOX5 | 1.06 (1.03–1.09) | <0.001 | 1.06 (1.03–1.09) | <0.001 |
Cox regression analysis of the association of ALOX5 and clinical variables with OS in the CGGA-LGG dataset.
Abbreviations: OS, overall survival; HR, hazard ratio; 95% CI, 95% confidence interval.
FIGURE 6
High methylation level of arachidonate lipoxygenases 5 was related to favorable prognosis
DNA methylation is an epigenetic alteration that plays an essential role in the development of several cancers (). Pearson correlation test unveiled the negative correlation between ALOX5 mRNA expression and its methylation level (Figure 7A). Using the DNA methylation data in the TCGA-LGG cohort, we found that low methylation level groups of cg06127294, cg06935264, cg09832911, cg10069493, cg10909790, cg12063947, cg14514237, cg15590007, cg19021,328, cg19517653, cg22770815, cg23054840, cg23084016, and cg24302529 led to worse OS than their related high methylation level groups (all p < 0.05) (Figure 7B) (Table 6). The result suggested that the low methylation level of ALOX5 might affect its gene expression, and lead to unfavorable clinical outcomes for LGG patients.
FIGURE 7
TABLE 6
| Methylated sites | Hazard ratio | 95% CI-low | 95% CI-high | p-Value |
|---|---|---|---|---|
| cg06127294 | 0.43 | 0.21 | 0.87 | 0.020 |
| cg06935264 | 0.36 | 0.18 | 0.72 | 0.003 |
| cg09832911 | 0.44 | 0.23 | 0.87 | 0.020 |
| cg10069493 | 0.48 | 0.24 | 0.96 | 0.030 |
| cg10909790 | 0.23 | 0.11 | 0.49 | <0.001 |
| cg12063947 | 0.35 | 0.15 | 0.81 | 0.010 |
| cg14514237 | 0.38 | 0.19 | 0.75 | 0.004 |
| cg15590007 | 0.34 | 0.17 | 0.66 | 0.001 |
| cg19021328 | 0.37 | 0.18 | 0.77 | 0.006 |
| cg19517653 | 0.40 | 0.20 | 0.80 | 0.007 |
| cg22770815 | 0.37 | 0.18 | 0.75 | 0.004 |
| cg23054840 | 0.37 | 0.19 | 0.74 | 0.003 |
| cg23084016 | 0.35 | 0.17 | 0.71 | 0.002 |
| cg24302529 | 0.32 | 0.16 | 0.64 | <0.001 |
Summary of the Kaplan-Meier curve data for the 14 methylated sites of ALOX5 in the TCGA-LGG dataset.
Abbreviations: 95% CI, 95% confidence interval.
High immune cell infiltration level was correlated with poor prognosis
The tumor microenvironment (TME) plays an essential role in the progression and pathological features of malignancies (). To explore the relationship between ALOX5 expression and TME in LGG, we assessed the association of ALOX5 expression with Immune score, Stromal score, and Estimate score using the ESTIMATE algorithm. The significant positive associations were illustrated in Figures 8A–C. To exhibit the significance of the TME in LGG, Kaplan-Meier analysis was performed to evaluate the correlation between the Immune score-, Stromal score-, and Estimate score-based groups, and sample survival. The Immune score-high, Stromal score-high, and Estimate score-high groups had decreased survival probability (all p < 0.001) (Figures 8D–F). Consistent with this, the ALOX5 mRNA expression was positively related to B cell, CD8 + T cell, CD4 + T cell, macrophage, neutrophil, and dendritic cell (all p < 0.001) (Figure 9A). In addition, patients with high immune cell infiltration levels had unfavorable clinical outcomes (Figure 9B). These results suggested that the elevated level of ALOX5 might affect the prognosis of LGG patients partially due to immune infiltration.
FIGURE 8
FIGURE 9
Arachidonate lipoxygenases 5 is involved in immune-related pathways
To reveal the pathological role of ALOX5 in LGG, the DEGs between the high- and low- ALOX5 expression groups were identified using the TCGA-LGG data and were shown in the volcano plot (Figure 10A) and the heat map (Figure 10B). As for biological processes, these DEGs were mainly involved in regulation of immune system process, cell activation, and immune effector process (Figure 10C). As for cellular component, they mainly participated in cell surface, side of membrane, and vesicle membrane (Figure 10D). In terms of molecular function, they were mainly involved in molecular transduce activity, signaling receptor binding, and immune receptor activity (Figure 10E). The major KEGG pathways were cytokine-cytokine receptor interaction, cell adhesion molecules, Toll-like receptor signaling pathway, NOD-like receptor signaling pathway, and B cell receptor signaling pathway (Figure 10F).
FIGURE 10
To decipher the potential mechanisms of ALOX5, we performed GSEA and found that several immune-related pathways such as B cell receptor signaling pathway, Toll like receptor signaling pathway, Nod like receptor signaling pathway, and T cell receptor signaling pathway were enriched in high ALOX5 expression group (Figure 11A). Therefore, we further explored the relationship between the ALOX5 and immune-related gene sets. ALOX5 was significantly positively correlated not only with immune activation-related genes, but also with immune checkpoint-related genes (Figures 11B,C). These results revealed that ALOX5 might affect the progression of LGG through activating the immune- and inflammation-related pathways.
FIGURE 11
Drug sensitivity analysis
Genes associated with the sensitivity and resistance of cancer drugs have been extensively studied (; ). We predicted the relationship between ALOX5 expression level and drug sensitivity using the GSCA database. ALOX5 was negatively correlated with the sensitivity to various chemotherapeutic drugs such as GSK690693, I-BET-762, PHA-793887, and afatinib (Figure 12).
FIGURE 12
Discussion
Herein, we found that the expression of ALOX5 was higher in LGG than that in normal tissue, and its overexpression was closely associated with poor prognosis of patients with LGG. Patients with low methylation levels of ALOX5 had unfavorable clinical outcomes, while those with high immune cell infiltration levels presented shorter survival. The enrichment analysis revealed the possible immune-related pathways regulated by ALOX5. Finally, our result demonstrated that ALOX5 mRNA expression was negatively correlated with drug sensitivity.
The specific function of ALOX5 is to catalyze AA into 5-HPETE, then into 5-HETE, and further into 5-oxo-eicosatetraenoic acid, which was further catalyzed into LTA4 generating proinflammatory factors including LTB4 and LTC4 by LTA4 hydrolase and LTC4 synthase, respectively (). In addition, ALOX5 is an essential enzyme that mediates lipid peroxidation by producing lipid peroxides (). Excessive lipid peroxidation is easy to occur in phospholipids, contributing to membrane rupture and evoking cell death such as pyroptosis, ferroptosis, and apoptosis (). Pyroptosis is a proinflammatory form of programmed cell death, which is dependent on the activity of caspase acid-specific proteases. In the coupling of the amino-terminal and carboxy-terminal linkers of gasdermin D (GSDMD) by caspases, GSDMD is transferred to the membrane and perforated, causing moisture penetration, cell swelling, and release of inflammatory factors followed by pyroptosis (). Ferroptosis is executed by phospholipid peroxidation which is a process relying on the transition metal iron, phospholipids containing polyunsaturated fatty acid chains (PUFA-PLs), and reactive oxygen species (ROS). Proinflammatory stimuli could alter the cellular metabolism, influencing levels of PUFA-PLs and ROS. Ferroptosis may also have physiological functions in tumor suppression and immune surveillance (; ). The activation of the cell death pathway can also trigger different immune and inflammatory reactions through the release of damage-associated molecular patterns (DAMPs) such as HMGB1, IL-33, IL-18, IL-1α, and IL-1β, further evoking inflammation and cell death (; ). Inflammation and cell death can affect each other, regulating organ homeostasis, and once dysregulated, it leads to pathological conditions including cancer and therapeutic resistance (). Therefore, the authors speculated that high ALOX5 expression might affect the poor prognosis of patients with LGG by triggering immune and inflammatory reactions and activating the cell death pathway. On the other hand, LGG patients with high ALOX5 expression might exhibit strong therapeutic resistance, thus having shorter survival. The possible mechanism of ALOX5 involvement in LGG was shown in Figure 13.
FIGURE 13
Abnormal methylation occupies a crucial role in the induction and progression of cancer including LGG (; ). Hypermethylation of the CpG sites in promoters frequently results in transcriptional silencing, while the hypomethylation of CpG sites in a gene body typically leads to upregulated gene expression (). This may explain that ALOX5 methylation had a strong negative connection with its mRNA expression. High expression of ALOX5 in LGG tissues might be affected by its low methylation level. Subsequently, we investigated the prognostic significance of ALOX5 DNA methylation and found that hypomethylation of 14 methylated sites in the ALOX5 gene correlated with a worse prognosis. Similarly, the methylation level of GALNT2/14 was negatively correlated with its expression, and patients with GALNT2/14 hypomethylation had shorter OS in lung adenocarcinoma (). Therefore, DNA methylation of ALOX5 might provide additional insight into the treatment and prognosis of LGG.
Increasing evidence proves that tumor-infiltrating immune cells partially influence the prognosis of LGG (; ; ). Tumor-infiltrating immune cells are part of the complex microenvironment closely related to LGG biological behavior and patient survival (; ). In this analysis, we observed that ALOX5 expression was positively associated with an immune score and several immune cells. Furthermore, high immune infiltration levels lead to a poor prognosis. Therefore, the authors inferred that ALOX5 expression might be implicated in the tumor immune microenvironment, leading to unfavorable clinical outcomes.
Following this, we performed the enrichment analysis to explore the underlying mechanism of ALOX5 in LGG. The results showed that ALOX5 was involved in immune and inflammation-related pathways such as the Toll-like receptor signaling pathway, Nod-like receptor signaling pathway, and antigen processing and presentation. Toll-like receptors are involved in the initiation of the innate and adaptive immune responses and the activation of these receptors triggers an inflammatory response that activates the regulatory pathways of innate and adaptive immunity (). Nod-like receptors are also involved in inflammatory responses and tumorigenesis (). Hence, ALOX5 expression might participate in the LGG progression by activation of these pathways. Finally, patients with high ALOX5 expression have higher immune checkpoint molecule expression, indicating that treatment with immune checkpoint inhibitors in the high ALOX5 expression group could attenuate immunosuppression to enhance existing antitumor immunity.
Our study only conducted IHC analysis, but the exploration of potential pathways provided insight into the future experiment. Second, the prognostic role of ALOX5 needs to be validated in prospective cohorts yet may require decades of follow-up.
In conclusion, the high ALOX5 expression serves as an independent prognostic predictor of poor prognosis in LGG. Additionally, ALOX5 expression might be regulated by its methylation level, and ALOX5 expression is also involved in the tumor immune microenvironment. Moreover, ALOX5 expression might be a robust biomarker for improving the prognosis of LGG patients.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://www.jianguoyun.com/p/DXtAeYQQp4jxChiLtNMEIAA.
Ethics statement
The studies involving human participants were reviewed and approved by The First Affiliated Hospital of Zhejiang Chinese Medical University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
R-HP and XZ contributed to the conception and design. Z-PC analyzed and interpreted the data. Y-JL contributed to the collection and assembly of data. All authors wrote and approved the final manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AnthonymuthuT. S.KennyE. M.BayirH. (2016). Therapies targeting lipid peroxidation in traumatic brain injury. Brain Res.1640, 57–76. 10.1016/j.brainres.2016.02.006
2
BlancheP.DartiguesJ. F.Jacqmin-GaddaH. (2013). Estimating and comparing time-dependent areas under receiver operating characteristic curves for censored event times with competing risks. Stat. Med.32 (30), 5381–5397. 10.1002/sim.5958
3
BurgosR. A.AlarconP.QuirogaJ.ManosalvaC.HanckeJ. (2020). Andrographolide, an anti-inflammatory multitarget drug: All roads lead to cellular metabolism. Molecules26 (1), 5. 10.3390/molecules26010005
4
ChenP.LiY.LiN.ShenL.LiZ. (2021). Comprehensive analysis of pyroptosis-associated in molecular classification, immunity and prognostic of glioma. Front. Genet.12, 781538. 10.3389/fgene.2021.781538
5
de AlmeidaB. P.ApolonioJ. D.BinnieA.Castelo-BrancoP. (2019). Roadmap of DNA methylation in breast cancer identifies novel prognostic biomarkers. BMC Cancer19 (1), 219. 10.1186/s12885-019-5403-0
6
DingM.LiY. A.LuZ.HouG. (2022). Identification of potential immune checkpoint inhibitor targets in gliomas via bioinformatic analyses. Biomed. Res. Int.2022, 1734847. 10.1155/2022/1734847
7
JiangT.MaoY.MaW.MaoQ.YouY.YangX.et al (2016). CGCG clinical practice guidelines for the management of adult diffuse gliomas. Cancer Lett.375 (2), 263–273. 10.1016/j.canlet.2016.01.024
8
JiangT.NamD. H.RamZ.PoonW. S.WangJ.BoldbaatarD.et al (2021a). Clinical practice guidelines for the management of adult diffuse gliomas. Cancer Lett.499, 60–72. 10.1016/j.canlet.2020.10.050
9
JiangX.StockwellB. R.ConradM. (2021b). Ferroptosis: Mechanisms, biology and role in disease. Nat. Rev. Mol. Cell Biol.22 (4), 266–282. 10.1038/s41580-020-00324-8
10
KhanA. Q.MohamedE. A. N.HakeemI.NazeerA.KuttikrishnanS.PrabhuK. S.et al (2020). Sanguinarine induces apoptosis in papillary thyroid cancer cells via generation of reactive oxygen species. Molecules25 (5), 1229. 10.3390/molecules25051229
11
KumarR.SinghA. K.KumarM.ShekharS.RaiN.KaurP.et al (2016). Serum 5-LOX: A progressive protein marker for breast cancer and new approach for therapeutic target. Carcinogenesis37 (9), 912–917. 10.1093/carcin/bgw075
12
KuoY. R.ChienC. M.KuoM. J.WangF. S.HuangE. Y.WangC. J. (2016). Endothelin-1 expression associated with lipid peroxidation and nuclear factor-κb activation in type 2 diabetes mellitus patients with angiopathy and limb amputation. Plast. Reconstr. Surg.137 (1), 187e–195e. 10.1097/PRS.0000000000001886
13
LapointeS.PerryA.ButowskiN. A. (2018). Primary brain tumours in adults. Lancet392 (10145), 432–446. 10.1016/S0140-6736(18)30990-5
14
LiuC.WeiD.XiangJ.RenF.HuangL.LangJ.et al (2020). An improved anticancer drug-response prediction based on an ensemble method integrating matrix completion and ridge regression. Mol. Ther. Nucleic Acids21, 676–686. 10.1016/j.omtn.2020.07.003
15
LiuP.LuZ.LiuL.LiR.LiangZ.ShenM.et al (2019). NOD-like receptor signaling in inflammation-associated cancers: From functions to targeted therapies. Phytomedicine64, 152925. 10.1016/j.phymed.2019.152925
16
LiuX.YangJ.ZhangY.FangY.WangF.WangJ.et al (2016). A systematic study on drug-response associated genes using baseline gene expressions of the Cancer Cell Line Encyclopedia. Sci. Rep.6, 22811. 10.1038/srep22811
17
LouisD. N.PerryA.WesselingP.BratD. J.CreeI. A.Figarella-BrangerD.et al (2021). The 2021 WHO classification of tumors of the central nervous system: A summary. Neuro Oncol.23 (8), 1231–1251. 10.1093/neuonc/noab106
18
MathurR.ZhangY.GrimmerM. R.HongC.ZhangM.BollamS.et al (2020). MGMT promoter methylation level in newly diagnosed low-grade glioma is a predictor of hypermutation at recurrence. Neuro Oncol.22 (11), 1580–1590. 10.1093/neuonc/noaa059
19
MerchantN.BhaskarL.MominS.SujathaP.ReddyA. B. M.NagarajuG. P. (2018). 5-Lipoxygenase: Its involvement in gastrointestinal malignancies. Crit. Rev. Oncol. Hematol.127, 50–55. 10.1016/j.critrevonc.2018.05.012
20
MolloyA. R.NajacC.ViswanathP.LakhaniA.SubramaniE.BatsiosG.et al (2020). MR-detectable metabolic biomarkers of response to mutant IDH inhibition in low-grade glioma. Theranostics10 (19), 8757–8770. 10.7150/thno.47317
21
Moradi-MarjanehR.HassanianS. M.FiujiH.SoleimanpourS.FernsG. A.AvanA.et al (2018). Toll like receptor signaling pathway as a potential therapeutic target in colorectal cancer. J. Cell Physiol.233 (8), 5613–5622. 10.1002/jcp.26273
22
MorganA. E.DaviesT. J.Mc AuleyM. T. (2018). The role of DNA methylation in ageing and cancer. Proc. Nutr. Soc.77 (4), 412–422. 10.1017/S0029665118000150
23
MouchlisV. D.DennisE. A. (2019). Phospholipase A2 catalysis and lipid mediator lipidomics. Biochim. Biophys. Acta Mol. Cell Biol. Lipids1864 (6), 766–771. 10.1016/j.bbalip.2018.08.010
24
PengH.WangY.WangP.HuangC.LiuZ.WuC. (2022). A risk model developed based on homologous recombination deficiency predicts overall survival in patients with lower grade glioma. Front. Genet.13, 919391. 10.3389/fgene.2022.919391
25
PengQ.LiuY.KongX.XianJ.YeL.YangL.et al (2021). The novel methylation biomarker SCARA5 sensitizes cancer cells to DNA damage chemotherapy drugs in NSCLC. Front. Oncol.11, 666589. 10.3389/fonc.2021.666589
26
RobinX.TurckN.HainardA.TibertiN.LisacekF.SanchezJ. C.et al (2011). pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinforma.12, 77. 10.1186/1471-2105-12-77
27
SinghN. K.RaoG. N. (2019). Emerging role of 12/15-Lipoxygenase (ALOX15) in human pathologies. Prog. Lipid Res.73, 28–45. 10.1016/j.plipres.2018.11.001
28
TangD.ChenX.KangR.KroemerG. (2021a). Ferroptosis: Molecular mechanisms and health implications. Cell Res.31 (2), 107–125. 10.1038/s41422-020-00441-1
29
TangJ.ZhangC.LinJ.DuanP.LongJ.ZhuH. (2021b). ALOX5-5-HETE promotes gastric cancer growth and alleviates chemotherapy toxicity via MEK/ERK activation. Cancer Med.10 (15), 5246–5255. 10.1002/cam4.4066
30
WangF.TaoZ.TianZ.JinJ.DongJ.DaiY.et al (2022). CCR5 as a prognostic factor in lower-grade glioma is involved in the remodeling of the tumor microenvironment. Front. Genet.13, 874896. 10.3389/fgene.2022.874896
31
WangH.RenS.XuY.MiaoW.HuangX.QuZ.et al (2019). MicroRNA-195 reverses the resistance to temozolomide through targeting cyclin E1 in glioma cells. Anticancer Drugs30 (1), 81–88. 10.1097/CAD.0000000000000700
32
WeeninkB.DraaismaK.OoiH. Z.KrosJ. M.Sillevis SmittP. A. E.DebetsR.et al (2019a). Low-grade glioma harbors few CD8 T cells, which is accompanied by decreased expression of chemo-attractants, not immunogenic antigens. Sci. Rep.9 (1), 14643. 10.1038/s41598-019-51063-6
33
WeeninkB.van BrakelM.WijersR.Sillevis SmittP. A. E.FrenchP. J.DebetsR. (2019b). Lack of B and T cell reactivity towards IDH1(R132H) in blood and tumor tissue from LGG patients. J. Neurooncol144 (1), 79–87. 10.1007/s11060-019-03228-6
34
WiernickiB.DuboisH.TyurinaY. Y.HassanniaB.BayirH.KaganV. E.et al (2020). Excessive phospholipid peroxidation distinguishes ferroptosis from other cell death modes including pyroptosis. Cell Death Dis.11 (10), 922. 10.1038/s41419-020-03118-0
35
WojtonJ.MeisenW. H.KaurB. (2016). How to train glioma cells to die: Molecular challenges in cell death. J. Neurooncol126 (3), 377–384. 10.1007/s11060-015-1980-1
36
WuF.LiG. Z.LiuH. J.ZhaoZ.ChaiR. C.LiuY. Q.et al (2020). Molecular subtyping reveals immune alterations in IDH wild-type lower-grade diffuse glioma. J. Pathol.251 (3), 272–283. 10.1002/path.5468
37
YangH. B.JiangJ.LiL. L.YangH. Q.ZhangX. Y. (2018). Biomarker identification of thyroid associated ophthalmopathy using microarray data. Int. J. Ophthalmol.11 (9), 1482–1488. 10.18240/ijo.2018.09.09
38
YangJ.WeiX.HuF.DongW.SunL. (2022). Development and validation of a novel 3-gene prognostic model for pancreatic adenocarcinoma based on ferroptosis-related genes. Cancer Cell Int.22 (1), 21. 10.1186/s12935-021-02431-8
39
YuG.WangL. G.HanY.HeQ. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS16 (5), 284–287. 10.1089/omi.2011.0118
40
YuY.WangZ.ZhengQ.LiJ. (2021). GALNT2/14 overexpression correlate with prognosis and methylation: Potential therapeutic targets for lung adenocarcinoma. Gene790, 145689. 10.1016/j.gene.2021.145689
41
ZhangH. C.DengS. H.PiY. N.GuoJ. N.XiH.ShiX.et al (2022). Identification and validation in a novel quantification system of ferroptosis patterns for the prediction of prognosis and immunotherapy response in left- and right-sided colon cancer. Front. Immunol.13, 855849. 10.3389/fimmu.2022.855849
42
ZhaoX.WangY.LiJ.QuF.FuX.LiuS.et al (2022). RFC2: A prognosis biomarker correlated with the immune signature in diffuse lower-grade gliomas. Sci. Rep.12 (1), 3122. 10.1038/s41598-022-06197-5
43
ZivkovicN. P.PetroveckiM.LoncaricC. T.NikolicI.WaegG.JaganjacM.et al (2017). Positron emission tomography-computed tomography and 4-hydroxynonenal-histidine immunohistochemistry reveal differential onset of lipid peroxidation in primary lung cancer and in pulmonary metastasis of remote malignancies. Redox Biol.11, 600–605. 10.1016/j.redox.2017.01.005
Summary
Keywords
arachidonate lipoxygenases 5 expression, methylation, immune cell infiltration, low-grade glioma, prognosis
Citation
Pan R, Zhang X, Chen Z and Liu Y (2023) Arachidonate lipoxygenases 5 is a novel prognostic biomarker and correlates with high tumor immune infiltration in low-grade glioma. Front. Genet. 14:1027690. doi: 10.3389/fgene.2023.1027690
Received
25 August 2022
Accepted
16 January 2023
Published
27 January 2023
Volume
14 - 2023
Edited by
Fang Wang, The Second Affiliated Hospital of Harbin Medical University, China
Reviewed by
Liangfang Shen, Central South University, China
Zihao Liu, Shandong Provincial Hospital, China
Xing Fan, Capital Medical University, China
Updates
Copyright
© 2023 Pan, Zhang, Chen and Liu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ya-jun Liu, yajunl475@163.com
This article was submitted to RNA, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.