1. Introduction
Central nervous system (CNS) cancers are rare and heterogeneous tumors with diverse biology and genetics and account for about 3% of all cancers in the world and are more common in men [
1]. Gliomas constitute the most common type of brain tumors, comprising around 24.7% of all primary tumors in the brain and other CNS tumors and 74.6% of malignant cases [
2]. The most prevalent malignant histology is glioblastoma making up 14.2% of all tumors and constituting 50.1% of all malignant CNS tumors [
3]. New 2021 WHO CNS 5 classification has divided diffuse gliomas into adult type and pediatric type. The basic molecules for the integrated diagnosis of adult-diffuse gliomas are
IDH (Isocitrate dehydrogenase
), p53, ATRX (alpha-thalassemia/mental retardation, X-linked) &
1p/19q co-deletion. Adult diffuse gliomas (IDH-mutant astrocytomas, IDH-mutant, 1p/19q codeleted oligodendroglioma and IDH-wildtype Glioblastoma) are diffusely infiltrating brain tumors [
4]. High grade Gliomas predominantly wild type glioblastomas are extremely lethal neoplasms with a poor prognosis. Despite maximum neurosurgical resection and adjuvant therapy; temozolamide (TMZ), median survival barely extends to approximately 12 months [
5]. Unfortunately, as with other solid tumors, chemo- resistance is one of the major challenges in this regard. Various mechanisms have been described for chemo resistance. Evasion of apoptosis is one of the mechanisms for tumor progression and chemo resistance in gliomas [
6] probably through either deficiency in BAX or BAK or gain of Bcl2 or BclX [
7]. Some studies have also revealed that the resistance of cancer cells to chemotherapeutic drugs may be due to the up regulation of autophagy, thereby avoiding apoptosis [
8,
9].
Autophagy is a cellular degradation pathway for the breakdown and removal of impaired long lived proteins, as well as the elimination of organelles and pathogens[
10]. Moreover, it has an important homeostatic role in maintaining the cell viability in stressed or nutritionally deprived states by recycling the cytoplasmic constituent [
11]. Thus, it controls the quality and quantity of proteins and organelles. Dysfunctional autophagy contributes in various diseases, with cancer being one of them. Nevertheless, in cancer, role of autophagy is somewhat complicated and controversial. Autophagy is a double edge sword, has dual effects on cancer. On one hand, it promotes tumor cell survival, by breaking down macromolecules into smaller components like amino acids, fatty acids, and metabolic substances. Studies have suggested that autophagy is activated in glioblastoma (GBM) as a reaction to pathophysiological challenges like necrosis and an acidic milieu. Hypoxia, a characteristic feature of high-grade gliomas, is responsible for inducing autophagy in these tumors. As the tumor progresses and reacts to therapeutic interventions, cells must adjust their metabolism to endure in hypoxic and nutrient-deficient surroundings, which is commonly linked to chemotherapy and radiotherapy resistance. This process facilitates tumor growth and viability, contributing to the maintenance of intracellular metabolism. Under unfavorable hypoxic conditions, molecules such as hypoxia-induced factors (HIF2a and HIF2b), BECN1 (Beclin-1), and BNIP3 (BCL2/adenovirus E1B 19KDa interacting protein 2) serve as survival mechanisms, promoting the progression of GBM and resistance to anticancer treatment in vivo [12-15]. Hypoxia causes BECN1 phosphorylation via the HIF-1a/BECN1 signaling pathway[
16,
17]. Beclin-1 is a key regulatory protein in autophagosome formation, which binds to the class III PI–3 kinase Vps34, thus facilitating the progression of autophagic process [
18,
19]. Remarkably, in the presence of hypoxia, the initiation of autophagy through BNIP3/BNIP3L serves as a survival mechanism, fostering the progression of GBM and enhancing resistance to anticancer treatments in vivo [
20]. A recent study showed that in U87 cells exposed to hypoxic conditions, the knockdown of HIF1A significantly reduced BNIP3 expression[
21]. This implies that tissues with low oxygen levels are prone to show heightened autophagy, indicating that an increased autophagy flux may serve an adaptive function [
15,
22]. On the other hand, autophagy can also suppress tumor progression in development, as impaired autophagy causes oxidative stress, triggers DNA damage responses and results in genomic instability, a known cause for tumor initiation. Such dual effects on tumor behavior are most likely context dependent [23-26]. Moreover research indicated that prolonged hypoxic stress in vitro, up regulates a pro-apoptotic Bcl-2 protein (BNIP3), leading to hypoxia-dependent autophagic cell death (ACD) in GBM cells [
27]. In terms of the mechanism, the elevated expression of BNIP3 leads to the release of BECN1 from complexes with Bcl-2 or Bcl-xL, thereby enabling BECN1 to activate autophagy[
28]. Hence, autophagy has been found to show conflicting functions in tumor onset and advancement. Thus, manipulation of autophagy pathways as means of cell death in cancer has led to the exploration of both inhibitors and inducers. Keeping in mind the contradictory results of autophagy in cancer, we investigated expression of autophagy markers and autophagy related genes (ATG) in our cohort of low and high grade gliomas.
3. Discussion
Autophagy is primarily a stress response process, and most of body tissues need autophagy to eliminate accumulated damaged organelles and unfolded proteins to maintain normal homeostatic milieu. In cancer, autophagy plays both tumor suppressing and tumor promoting effects depending on the specific context, thus could influence on the prognosis, either favorably or adversely. We investigated whether autophagy plays anti-tumoral or pro-tumoral functions in various grades of gliomas.
In this current research, we assessed the immunohistochemical expression of autophagy markers as well as mRNA levels of ATG in different grades of adult-type diffuse gliomas, while also assessing their correlation with clinicopathological parameters. Our investigation revealed a significant connection between the immunohistochemical expression of
LC3B, SQSTM1/p62, DRAM, and the overall autophagy status, particularly in relation to tumor grades. Notably, among the 25 positive cases of
LC3B, only two were observed in the LGG. Similarly, out of the 20 strongly positive cases of
SQSTM1/p62, none was found in the LGG category, and this difference achieved statistical significance (as indicated in
Table 3). These results were in accordance with other studies which reported that
LC3B and
SQSTM1/p62 was highly expressed in high grade gliomas [
29,
30,
31]. Remarkably there is statistically significant association between the autophagy markers
LC3B and
SQSTM1/p62 (p=0.000), highlighting the interaction among each other. These results may be attributed to the active participation of
LC3B and
SQSTM1/p62 in the process of autophagy, given their roles as essential structural components of the autophagosome[
32] and recent literature also revealed that that
SQSTM1/p62 overexpression is capable of promoting mitochondrial and classical macroautophagy [
33], which promote tumor progression and chemo-resistances. Furthermore increased expression of
LC3B and
SQSTM1/p62 in HGG compared with LGG may be in response to confer stress tolerance which is definitely more in HGG, serves to maintain tumor cell survival [
34]. Hypoxia, which is major hallmark of Glioblastoma, is known to induce autophagy in these tumors, which eventually produce survival mechanism facilitating the breakdown of various cellular components to generate ATP and metabolic precursors as a means to cope up with stress. This highlights the crucial role of autophagy in protecting cells against stressful conditions[
35]. In a study conducted by Deng et al; it was reported that both mRNA and protein levels of
SQSTM1/p62 were elevated in human glioma tissues, and the suppression of
SQSTM1/p62 had an anti-tumor effect on glioma cells [
34]. These findings may support that high levels of
LC3B and
SQSTM1/p62 expression or prompted autophagy correlated with advanced tumor grade and aggressiveness with the potential role of autophagy as tumor enhancer, thus supporting our results.
We also looked another important autophagy modulator
DRAM1 in our study and found its overexpression in 28.2% of the cases, although it was not significantly associated with tumor grades but higher expression of
DRAM1 was present in high grade gliomas, when compared with LGG.
DRAM 1 primarily resides in lysosomes and is frequently down regulated in various human cancers[
36].
DRAM1 enhances lysosomal acidification and facilitates the fusion of lysosomes with autophagosomes, thus triggers the autophagy. Moreover
DRAM1 plays pivotal role in governing the association of
SQSTM1/p62 with autophagosomes and its subsequent degradation through autophagy. Therefore, reduced
DRAM1 expression might be associated with decreased
SQSTM1/p62 localization within autophagosome, highlighting the role of
DRAM1 in
SQSTM1/p62 -mediated autophagy. Our findings align with the study by Geng et al; which also reported reduced
DRAM1 expression in non-small cell lung carcinoma, correlating with an unfavorable prognosis [
37]. However, another study has shown high expression of both
DRAM1 and
SQSTM1/p62 in Glioblastoma, where they regulate cell migration and invasion, and are associated with shorter or poor overall survival [
38]. These conflicting results could arise from variations in demographics and sample sizes. Thus utilizing the expression levels of
LC3B, SQSTM1/p62 and
DRAM1 to assess autophagy status in glioma patients could be potential predictive markers.
While investigating the autophagy gene expression patterns in both low (grade 2) and high (grade 3 and 4) grades of tumors, we noticed a remarkable and statistically significant up-regulation in autophagy critical kinase, vacuolar protein sorting 34 (
Vps34) gene expression thus indicating a tumor promoter role of
VPS34 gene in autophagy. Previous research has shown that
VPS34 is involved in initiating autophagy through its interaction with Vps15/Atg14/UVRAG/Beclin1[
39]. Consequently, this has gathered considerable attention as a potential target for inhibiting autophagy [
40,
41] .
We also evaluated the autophagy upstream pathway,
PI3K/AKT/mTOR and recorded higher expression of
mTOR transcripts in HGG as compared to LGG, which is a negative regulator of autophagy but it didn’t reach significant value(p = 0.07).
PI3K/Akt/mTOR signaling pathway is a frequently disrupted pathway across different cancer types [
42,
43] and abnormal activation of this pathway has been associated with tumor development, progression, invasion and metastasis [
44] and is indeed activated in glioma cells [
45]. Nonetheless the
PTEN, PI3K, AKT, mTOR, Beclin-1, UVRAG, ULK1 and
ULK2 genes were remained unchanged when compared in both grades. The significant increase in transcriptional expression of VPS34 implies an increased initiation of autophagy, as VPS34 plays a crucial role in the nucleation of autophagosomes. Conversely, the absence of substantial transcriptional alterations in other genes suggests that these components may not be profoundly influenced at the transcriptional level in high-grade gliomas. Nevertheless, it is important to consider the potential involvement of post-transcriptional, post-translational, or alternative regulatory mechanisms. The specific modulation of autophagy in high-grade gliomas appears to be contingent on the unique context of these tumors. Gliomas exhibit heterogeneity, and the signaling pathways governing autophagy may be subject to the diverse genetic and epigenetic profiles inherent to individual tumors.
The spearman’s correlation between autophagy genes was examined in grades 2, 3 and 4 glioma cases and we found significant positive correlation among different genes. It’s essential to highlight that these observed correlations were statistically significant, but whether they have a biological basis or not remains uncertain. For instance in grade 2 gliomas
PTEN showed a strong positive correlation with autophagy genes in line with a study by Errafiy Rajaa [
46]. However the absence of
PTEN correlation in grade 4 in current study highlights loss of
PTEN which is a hallmark of GBM, could be due to mutation or promoter methylation of
PTEN gene [
47,
48].
Similarly in HGG the significant positive correlation between
VPS34 and
ULK1, ULK2, UVRAG, Beclin1, PI3K genes can also be justified biologically, as literature prevailed that
VPS34 kinase forms a stable complex with Beclin1 and p150 which serves as binding partner for
ATG14L, UVRAG and AMBRA [
49] which are responsible for promotion of autophagy. Thus significant positive correlation indicates that expression of autophagy genes tends to change consistently together, providing insights into the regulatory connections within the autophagic pathways. This understanding may have implications for comprehending the underlying mechanisms influencing tumor progression.
Nevertheless, a limitation of our study is the absence of autophagy flux assessment in formalin fixed paraffin embedded tissue blocks (FFPE) samples. Elevated expression levels of LC3B and SQSTM1/p62 do not consistently correlate with an overall augmentation in autophagy and may be attributed to a potential hindrance in autophagy at the later stages of autophagosome processing. Thus need for a more comprehensive assessment of autophagy, particularly considering the dynamic nature of the process. Moreover, ATG5, ATG12, ATG 7 and ATG 4 genes are crucial for two conjugation systems (Atg8-Atg4 and Atg12-Atg5) involved in autophagy and would provide valuable information, as these processes contribute to the formation and elongation of autophagosomes, facilitating the degradation and recycling of cellular components. Unfortunately, due to financial constraints, we were unable to conduct this analysis. However, future comprehensive studies utilizing diverse techniques will be invaluable in addressing these limitations.
4. Materials and Methods
4.1. Patients Selection
Patients of adult –diffuse gliomas were enrolled from Histopathology department of Dow Diagnostic Reference and Research Lab (DDRRL), Dow University of Health Sciences (DUHS). The research protocol received approval from the Institutional Review board at Dow university of health sciences [Ref: IRB-1150/DUHS/Approval/2018]. Clinical information was recorded. The inclusion criteria required the inclusion of patients diagnosed with adult diffuse gliomas of all grades. Nevertheless, patients diagnosed with pediatric-type diffuse gliomas and astrocytomas with circumscribed morphology were excluded from the study. Initially, there were 50 prospective cases of adult diffuse gliomas that were included in the study. Nonetheless, due to the use of multiple immunohistochemical stains and the extraction of nucleic acids from tissue samples, some cases ran out of available tissue. Consequently, 39 patients were finally included in the study. The mRNA and protein expression of autophagy-related genes (ATG) were determined using Quantitative real-time PCR (qPCR) and IHC.
4.2. Tissue processing for histopathological analysis
Biopsy specimens of brain tumors were received in histopathology department from different hospitals and centers of city through its collection point service. For histopathological examination, every specimen was placed in an automated tissue processor (Thermo scientific, model EXCELSIOR AS) for 12 hours, where they were exposed to graded series of alcohol as 70%, 95% 100% followed by xylene, and then paraffin. The paraffin blocks obtained were sliced into sections measuring 4-5μm thickness for subsequent hematoxylin and eosin staining (H&E). Before staining slides were deparaffinized and sections were treated with xylene, followed by hydration through decreasing concentration of ethanol, as 100%, 70% and 50%. Slides were then cleaned with water before being submerged in H & E for staining. To eliminate excess water from the slide after staining, the sections were submerged in 60%, 80%, and 100% ethanol. Next, slides were mounted by coverslip with Enthelan® after being submerged in xylene (Merck, Darmstadt, Germany). Histopathological examination of the H &E-stained slides at a microscopic level was conducted utilizing the Nikon Eclipse E200 optical microscope (Nikon Instruments Inc. in Tokyo, Japan) and the tumor characterized as adult diffuse gliomas were enrolled in the study. The cases of diffuse gliomas were categorized according to CNS5 classification into distinct subtypes and malignancy grades using histopathological features such as cellularity, atypia, necrosis, micro-vascular proliferation and mitosis.
4.3. Tissue processing for immunohistochemical examination
IHC was carried out using Autostainer Link 48 (Dako North America Inc., USA, S no AS3006D1307), on formalin-fixed paraffin embedded (FFPE) tissue blocks for LC3B, p62 and DRAM antibodies. Sections were sliced into 3-4μm and affixed on charged slides (EnVision FLEX visualization systems) and dried for 60 minutes at 60-70-ºC. Following deparaffinization with xylene rinsed in decreasing concentration of ethanol and then rehydrated in distilled H2O. To facilitate antigen unmasking, retrieval solution (EnVision Flex Target Retrieval solution, pH 9.0, TRIS HCL) was used for 20 minutes in a preheated water bath. To quench endogenous peroxidase activity, slides were immersed in peroxidase blocking solution (EnVision Flex Peroxidase blocking reagent, RTU) for 10 minutes. TBST (Tris Buffer saline with Tween 20, EnVision Flex wash Buffer) was used for washing then sections were incubated for 30-60 minutes at room temperature with primary antibodies LC3B ((ABCAM, ab51520: 1: 800), p62 ((Invitrogen, clone: SOSTM1: 1: 25) and DRAM (Invitrogen, 1; 25). Subsequently, the sections were subjected to a 30-minute treatment at room temperature with a secondary antibody (EnVision Flex/HRP, RTU), after washing with PBS buffer. DAB solution 3, 3-diaminobenzidine tetrahydrochloride solution (EnVision Flex DAB+ chromogen) was applied for 10 minutes to the sections on the slides to reveal the color of antibody staining and counterstained with hematoxylin. Slides were washed for 10 minutes in running water. Sections were then dehydrated (in graded alcohol as 80%, 90%, and 100%), cleaned and cover-slipped using DPX mounting solution.
4.4. Immunohistochemical evaluation of autophagy markers
LC3B and
SQSTM1/p62 punctate/dotted cytoplasmic staining was considered as positive [
50,
51]. The immunopositivity was scored on the basis of strong intensity and percentage positive glioma cells. Immunopositivity of >50% of tumor cells was considered as positive for
LC3B, whereas for
SQSTM1/p62 it was considered positive when it exceeded >30%[
52]. The intensity and percentage scores were multiplied to find the
DRAM overexpression. The intensity of
DRAM1 staining was scored as follows: 0 (no staining), 1 (weak staining), 2 (moderate staining), and 3 (strong staining). Percentage scores were assigned as follows: 1 (0–25%), 2 (26–50%), 3 (51–75%), and 4 (71–100%). Every sample was assigned a score, which was multiplied to get a total value that ranged from 0 to 12. Score of ≥4 was defined as
DRAM1 overexpression, and score of < 4 was categorized as weak or negative expression [
53]. The immunohistochemical staining was independently scored by two pathologists at 40x objective magnification and then discrepancies were discussed on a multi head microscope and final scores were determined. Autophagy status was considered as positive when two out of the three autophagy-associated proteins were detected in each sample[
54].
4.5. RNA extraction and cDNA synthesis
For the RT-qPCR analysis, total RNA was isolated from FFPE blocks via Pure Link FFPE, total isolation (Invitrogen; Thermo Fisher Scientific, Inc) according to the manufacturer’s protocols [
55]. Subsequently DNase treatment was carried out for any DNA contamination. This involved combining 1μg of RNA template with 1μl of reaction buffer containing MgCl2, 1μl of DNase-I, RNase-free (Thermo Fisher Scientific, Cat. No. EN0521), and nuclease-free water in a 0.2 ml tube. The resultant mixture was incubated in the Master cycler X50a (Eppendorf, Germany) for 30 min at 37 °C. To prevent RNA degradation, following the DNase-I treatment, we introduced 1μl of 50 mM EDTA and the samples were incubated at 65 °C for 10 minutes. The RNA integrity was evaluated with a nanodrop, and subsequently, cDNA was generated according to the manufacturer's instructions using the Revert-Aid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Cat. No. K1612). The resulting cDNA was stored at -20°C for future applications.
4.6. Gene expression analysis via quantitative real time (qPCR)
qPCR was conducted to assess the genes expression levels of
ULK1, ULK2, Beclin 1, UVRAG, VPS34, PTEN, PI3K, AKT, and
mTOR, using PCR kit (PowerUp™ SYBR™ Green Master Mix (Thermo Fisher Scientific, Cat. No A25742), and primers (Eurofins, USA). β-actin, a housekeeping gene was used for result normalization in a qPCR assay, utilizing the corresponding primer sets. To perform qPCR analysis, cDNA samples were used and 10uL reaction mixture was prepared. Mixture comprised of 2 microliter cDNA, 1μl primers of both forward and reverse primers, and 5μl of PowerUp™ SYBR™ Green Master Mix and 1-2 microliter of nuclease free water. The thermal cycling conditions for the reaction were as follows: an initial 2 min hold at 50°C, another 2-minute hold at 95°C, and then 40 cycles of denaturation for 15 sec at 95°C followed by annealing at 60°C for 1 minute. 2-ΔΔCt (Livak’s method) was used to analyze the relative changes in gene expression and multivariate Anova test was performed for statistical significance. Following the differential expression of ATG genes, the log2 fold change was calculated. Grade 2 gliomas, which are considered low grade, were chosen as the baseline group for computing the log2 fold change (FC).
Table 8 shows the list of primers that were utilized.
4.7. Statistical Analysis
Descriptive statistics were used to express as means with standard deviation. Pearson’s Chi-square test was executed for the association of demographics, clinical pathologic parameters, and molecular markers with tumor type and grade. Multivariate Anova test was performed test to identify any notable difference in the relative gene expression among all the examined ATG genes in 2, 3 & 4 WHO grades of diffuse gliomas and p < 0.05 was considered as significant. Spearman's correlation test was executed to determine the correlation between autophagy genes in grade 2, 3 and 4. We performed all the analyses using IBM SPSS version 24, and used a significance threshold of p < 0.05 to determine statistical associations.