AL355338 acts as an oncogenic lncRNA by interacting with protein ENO1 to regulate EGFR/AKT pathway in NSCLC

doi:10.21203/rs.3.rs-815448/v1

Background

Non-small cell lung cancer (NSCLC) is a malignancy with considerable morbidity and mortality. Abnormal metabolism is a hallmark of cancer; however, the mechanism of glycolysis regulation in NSCLC progression is not completely understood. Recent studies suggest that some dysregulated long non-coding RNAs (lncRNAs) play important roles in tumor metabolic reprogramming.

Methods

To identify glycolysis-associated-lncRNAs in NSCLC, we compared RNA-sequencing results between high ¹⁸F-fluorodeoxyglucose (FDG)-uptake NSCLC tissues and paired paratumor tissues. The transcript abundance of AL355338 in 80 pairs of clinical samples was evaluated by quantitative real-time PCR assay and fluorescence in situ hybridization. The biological role of AL355338 on NSCLC cells were evaluated by functional experiments in vitro and in vivo. Moreover, RNA pull-down, mass spectrometry and RNA immunoprecipitation (RIP) assays were used to identify the protein interacted with AL355338. Co-immunoprecipitation, in situ proximity ligation assays and western blotting were applied to define the potential downstream pathways of AL355338.

Results

AL355338 was an upregulated glycolysis-associated lncRNA in NSCLC. Functional assays revealed that AL355338 was critical for promoting aerobic glycolysis and NSCLC progression. Mechanistic investigations showed that AL355338 directly bound with alpha-enolase (ENO1) and enhanced the protein’s stability by modulating its degradation and ubiquitination. A positive correlation was observed between AL355338 and ENO1 in NSCLC, and ENO1 was subsequently confirmed to be responsible for the oncogenic role of AL355338. Furthermore, AL355338 was capable of modulating ENO1/EGFR complex interaction and further activating EGFR-AKT signaling.

Conclusions

This study indicates that AL355338 confers an aggressive phenotype to NSCLC, and targeting it might be an effective therapeutic strategy.

Cancer Biology

Long noncoding RNAs

aerobic glycolysis

AL355338

ENO1

EGFR

Among malignancies, lung cancer is one of the most lethal with high global mortality and incidence rates[1]. Non-small cell lung cancer (NSCLC) is the most commonly diagnosed type, accounting for ~ 85% of cases[2]. Early diagnosis of NSCLC is difficult, and most patients are diagnosed at advanced stages and have poor clinical outcomes[2]. In-depth investigation into the molecular alterations and their mechanisms that underlie NSCLC progression is extremely essential to develop and optimize therapeutic strategies.

Cancer cell growth depends on sufficient energy. Previous studies suggested that glucose metabolism reprogramming, also known as the Warburg effect, provided fuel for cancer’s malignant behavior[3]. ¹⁸F-fluorodeoxyglucose (FDG) is an analogue of glucose and radioactive tracer that enables tracking of glucose uptake in solid tumors in positron emission tomography/computed tomography (PET/CT) imaging, which is the gold standard for quantifying glycolytic metabolism[4]. Enhanced glucose uptake ability has been observed in patients with various types of cancer, and especially in NSCLC[5]. Cancer biology studies using genomic techniques revealed that the orchestrated performance of oncogenes and downstream pathway activation exacerbate metabolic disorder[6]. However, the underlying mechanistic details underlying glycolysis regulation in NSCLC progression remain to be determined.

Long non-coding RNAs (lncRNAs) are RNA transcripts (> 200 nucleotides) with limited coding potential[7] that are abnormally expressed in all stages of cancer development and play vital regulatory roles in oncogenic processes[8]. A recent study reported that lncRNA dysregulation is closely related to tumor metabolic reprogramming[9]. Further work is needed to identify glycolysis-related-lncRNAs in NSCLC and clarify their potential roles and specific mechanisms.

Alpha-enolase (ENO1) is an enolase isoform originally characterized as a crucial glycolytic enzyme catalyzing the dehydration of 2-phosphoglycerate into phosphoenolpyruvate; it also generates ATP during glycolysis[10]. Besides its well-known enzymatic functions, decades of research have shown that ENO1 is aberrantly expressed in multiple tumor types and participates in diverse cellular processes through non-metabolic mechanisms[11]. Interestingly, some studies reported that lncRNAs interact with ENO1 and affect its expression or activity[12, 13], highlighting a novel approach to elucidate the non-metabolic roles of metabolic genes in cancer.

Recently, we aimed to identify glycolysis-associated-lncRNAs through RNA-sequencing (RNA-seq) of high ¹⁸F-FDG-uptake NSCLC tissues and paired paratumor tissues. AL355338 was the top-scoring upregulated lncRNA. In this present study, we observed that high AL355338 expression positively correlated with NSCLC clinical stage and worse patient prognosis. Functional assays demonstrated the influences of AL355338 silencing on NSCLC proliferation, metastasis, and glycolysis in vitro and in vivo. Mechanistically, our evidence indicated that AL355338 could directly interact with ENO1 and stabilize the protein, preventing its ubiquitination and degradation. AL355338 exerted its biological functions partially by upregulating ENO1 expression. Further experiments showed that AL355338 positively modulated the interaction between ENO1 and epidermal growth factor receptor (EGFR), subsequently activating EGFR/AKT signaling. Collectively, these results strongly suggest that the AL355338-ENO1-EGFR/AKT regulatory axis may harbor NSCLC diagnostic biomarkers and therapeutic targets.

Patient and clinical information

Eighty paired NSCLC tissues and paratumor tissues were collected from the surgery specimen archives of the Second Affiliated Hospital of Soochow University. Patient follow-up continued to December 2020. No recruited patients received preoperative chemotherapy or radiotherapy. This research was approved by the institutional clinical research ethics committee of the Second Affiliated Hospital of Soochow University. Patient stratification was based on the new TNM staging criteria in lung cancer (8th edition)[14].

Cell culture and treatment

The human normal lung epithelial cell line HBE and four human NSCLC cell lines (A549, PC9, H1650, H1299) were obtained from ATCC (Manassas, VA, UA) and cultured in Dulbecco's modified Eagle's medium (Gibco, Grand Island, NY, USA) supplemented with 10% fetal bovine serum (Gibco) at 37°C, in a humid 5% CO₂ incubator. Two small interfering RNA (siRNA) constructs (GenePharma, Shanghai, China) were designed to knock down AL355338 expression, and non-specific siRNA served as the negative control; all siRNAs were transfected with Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA). For ectopic AL355338, ENO1, and EGFR expression, the appropriate cDNA was amplified with reverse transcription polymerase chain reaction (RT-PCR) and subcloned into a pcDNA-4/TO vector. Neofect transfection reagent (Neofect biotech, Beijing, China) was used for plasmid transfection. The sh-AL355338, sh-ENO1 and sh-NC lentiviruses were purchased from GenePharma. The siRNA and sh-RNA sequences are listed in Table S1.

Total RNA isolation and quantitative RT-PCR (qRT-PCR) assays

Total RNA was extracted from cells or tissues with TRIzol reagent (Omega, Norcross, GA, USA) according to the manufacturer’s instructions. RNA quality and quantity were assessed with the NanoDrop2000 platform (Thermo Fisher Scientific, Waltham, MA, USA). cDNA Synthesis kits (TaKaRa, Otsu, Japan) were used for mRNA reverse transcription. Real-time PCR analysis was performed with SYBR Green PCR Kits (TaKaRa) on the 7500 Sequence Detection System (Applied Biosystems, Foster City, CA, USA). The amplified transcript levels of genes were normalized to β-actin with the optimized comparative 2^-ΔΔCt value method. The primer sequences are in Table S1.

Nuclear-cytoplasmic fractionation

The nuclear and cytoplasmic fractions were generated with a PARIS™ Kit (Ambion, Austin, TX, USA) according to the manufacturer protocol. Nuclear and cytoplasmic RNA levels were normalized to U6 and β-actin, respectively using real-time PCR.

Fluorescence in situ hybridization (FISH)

FISH analyses were performed with RNA FISH kits (RiboBio, Guangzhou, China) using a previously described method[15]. The AL35538 probe sequence was 5’-DIG-ACTTCGAGACCAGCATGGCCAACATGGTGAAGCCC-DIG-3’. Nuclei were stained with DAPI, and high-resolution images were taken with a confocal microscope (ZEISS, Oberkochen, Germany). Percentages of positive stained cells were quantified as: 1 = 0-25%, 2 = 26-50%, 3 = 51-75%, and 4 = 76-100%. Staining intensity scores were: 0 = none, 1 = weak, 2 = moderate, and 3 = strong. The final scores were calculated by multiplying the percentage and intensity scores of positive cells. Samples were divided into low (score 0-6) and high (score 7-12) expression groups.

Flow cytometry assay

Transfected cells were harvested and fixed in pre-cooled 70% ethanol for 1 h on ice. Then cells were washed with cold phosphate-buffered saline (PBS) three times and incubated with RNase for 30 min at 37°C. After staining the cells with propidium Iodide for 15 minutes at room temperature, they were analyzed with a flow cytometer (BD Biosciences, Franklin Lakes, NJ, USA). The percentages of cells in the G0/G1, S, and G2/M phases were calculated using ModFit software (Verity Software, Topsham, ME, USA).

Cell proliferation assay

Cell proliferation rates were measured with Cell Counting Kit-8 (bimake, Shanghai, China) according to the manufacturer’s protocol. Briefly, 5 × 10³ cells per well were seeded into 96-well plates. At 0, 24, 48, and 72 h after cell transfection, 10 μL CCK-8 solution was added to each well. After incubating he plates for 2 h at 37°C, the optical density at 450 nm was measured for each well. For colony formation assays, transfected or control A549 and PC9 cells were harvested 24 h after transfection, and 500 cells per well were seeded into 6-well plates. After 2 weeks of incubation, the colonies were fixed in 4% paraformaldehyde and stained with crystal violet solution. Cell proliferation was also assessed using ethynyldeoxyuridine (EdU) analyses (RiboBio) according to the manufacturer’s instructions as previously reported[16].

Cell migration and invasion assay

Wound healing assays were carried out as previously described [17]. Transwell assays were performed using Boyden chambers containing transwell membrane filters (Corning, Corning, NY, USA). The detailed protocol was published elsewhere[17]. Images were collected with an inverted microscope (Olympus, Tokyo, Japan) and analyzed with ImageJ software (National Institutes of Health, Bethesda, MD, USA). At least five random fields of view were analyzed in each chamber.

Measurement of glucose uptake, lactate production, and ATP

As a glucose analogue, ¹⁸F-FDG uptake reflects the intracellular glucose uptake level. In vitro experiments were performed as previously reported[16]. Lactate production was measured in cells that were seeded into 12-well plates in triplicate for 24 h and then refreshed with 1 mM glucose culture medium overnight. Culture supernatants were collected to measure lactate concentrations (Nanjing Jiancheng Bioengineering Institute, Nanjing, China), and lysed cell pellets were used to measure ATP levels (Nanjing Jiancheng) using the manufacturer’s instructions. ¹⁸F-FDG uptake and lactate production were normalized to number of cells in each sample; ATP levels were normalized to cell protein mass.

Western blot

Western blotting was performed following a standard protocol[18]. The antibodies used are listed in Table S2. Immunoreactive bands were detected using ECL western blot kits (Amersham Biosciences, Little Chalfont, UK).

Immunofluorescence (IF)

IF was performed as described elsewhere[18]. Protein expression and localization were observed under a confocal microscope (ZEISS). Antibodies used for IF are listed in Table S2.

RNA pull-down assay

Biotin-RNA pull-down assays were performed as described elsewhere[15]. Briefly, full-length sense, antisense, and serial deletion sequences of AL355338 were amplified with a T7-containing primer and reverse transcribed using MAXIscript™ T7 Transcription Kits (Thermo Fisher Scientific). Pierce™ RNA 3’ End Desthiobiotinylation Kits (Thermo Fisher Scientific) were used to label targeted RNA. The biotin-labeled AL355338 probe was incubated with total cell lysates of A549 cells using Pierce™ Magnetic RNA-Protein Pull-Down Kits (Thermo Fisher Scientific). Eluted proteins were purified by silver staining and subjected to mass spectrometry (MS) or western blotting.

RNA immunoprecipitation (RIP) assay

Magna RNA-binding protein immunoprecipitation kits (Millipore, Billerica, MA, USA) were used to determine specific protein binding with targeted RNA according to the manufacturer’s protocol and as described previously[15]. Brieﬂy, cells were lysed with RIP lysis buffer, and lysates were incubated overnight at 4°C with magnetic beads conjugated with human anti-ENO1 antibody or anti-mouse immunoglobulin G (IgG). The next day, co-immunoprecipitated RNA was extracted for qRT-PCR analyses. Total RNA (input controls) and IgG controls were simultaneously assayed to confirm that the RNA specifically bound to ENO1.

Co-immunoprecipitation (Co-IP)

Reciprocal exogenous and endogenous co-IP assays were performed as described elsewhere repor[18]. IP samples and input controls were analyzed by western blot. The co-IP antibodies are listed in Table S2.

Proximity ligation assay (PLA)

PLA (Duolink®, Sigma-Aldrich, St. Louis, MO, USA) was performed to evaluate interactions between ENO1 and EGFR as previously reported[19]. Briefly, A549 and PC9 cells were grown on glass coverslips in 24-well plates at 37 °C overnight. The next day, cells were fixed in 4% paraformaldehyde for 20 min and permeabilized with 0.1% Triton X-100 for 5 min. After blocking for 30 min, the cells were incubated with anti-ENO1 and anti-EGFR antibodies diluted in blocking buffer at 4°C overnight. The next day, cells were incubated with species-specific PLA probes at 37°C for 1 h, then 1× ligase for 30 min, and finally 1× polymerase diluted with amplification stock solution at 37°C for 100 min to amplify the fluorescence signal. Finally, coverslips were mounted on slides with anti-fluorescence quenching sealing liquid containing DAPI. Images were acquired under confocal microscopy (ZEISS).

In vivo xenograft assays

Four-week-old male BALB/c nude mice were purchased from the Animal Center of Soochow University. Animal experiments were carried out according to the protocol approved by the Animal Care and Ethics Committee of the Second Affiliated Hospital of Soochow University. A total of 20 mice were subcutaneously injected with 5 × 10⁶ A549 cells in the right flank with sh-AL355339, and sh-NC was given in the left flank to establish an NSCLC xenograft model. Tumor volume was measured every fifth day and calculated according to the equation V = length × width²× 1/2. Micro-PET/CT (Super Nova® PET/CT, Pingseng, China) was performed 3 weeks after cell injection to measure ¹⁸F-FDG uptake by xenograft tumors. All tumor-bearing mice were fasted for 8 h prior to PET/CT. ¹⁸F-FDG (0.1 mL, 7.4 MBq) was injected via tail vein. After 30 min, they were anesthetized with 2% isoflurane in oxygen and immobilized for 15-min scans. ¹⁸F-FDG uptake was measured as maximum standard uptake (SUVmax) in the regions of interest covering the tumors. After 3 weeks, mice were humanely sacrificed by euthanasia. The mice were inhaled with CO₂ for 5 min. The CO₂ was gradually filled into an airtight chamber at 30% of the chamber volume/min rate. The xenograft tumors were harvested and measured.

In vivo metastasis assays

Stably transfected A549 cells (1 × 10⁶/0.2 ml PBS) were injected into the tail vein of nude mice. Micro-CT scans were performed 5 weeks later to observe lung metastasis. After 5 weeks, mice were humanely sacrificed by euthanasia. Following sacrifice, the lungs were surgically dissected and embedded in paraffin for hematoxylin and eosin (HE) staining. The numbers of pulmonary metastatic nodules were compared between groups.

Immunohistochemistry (IHC)

IHC was performed and quantified as previously described[20].

Statistical analysis

All data are expressed as mean ± standard deviation. Statistical analysis was performed with SPSS (version 22.0, IBM Corp., Armonk, NY, USA) and graphs were generated with GraphPad Prism software (version 8.0, GraphPad Inc., San Diego, CA, USA). Student’s t-tests or one-way analyses of variance were used to assess significant differences between groups. Categorical data were compared with Chi-square tests. Survival analyses were performed using Kaplan-Meier analyses.

AL355338 was upregulated in aggressive NSCLC and associated with poor patient prognosis

To identify glycolysis-related lncRNAs in NSCLC, we performed RNA-seq assay with ¹⁸F-FDG high-uptake NSCLC tissues and pair-matched normal tissues from our previous study [16]. We detected 222 upregulated genes (fold change >1.5, false discovery rate <0.005). The top 20 significantly upregulated lncRNAs (Figure S1A) were submitted for Gene Ontology analysis, and the results showed that they were most enriched in metabolic pathways (Figure S1B). To screen more functional lncRNAs, our next investigation used more stringent filtering criteria (>400 bp, not fully overlapped with other coding transcripts, efficient PCR amplification). In the current study, we focused on the uncharacterized AL355338, which was the top-scoring upregulated lncRNA in the RNA-seq data, and explored its function in NSCLC development. AL355338 is a 913-bp transcript on human chromosome 13 with 2 exons (Figure S1C). The Coding-Potential Assessment Tool (CPAT), PhyloCSF tool, and ORFfinder all predicted that AL355338 had no coding potential (Figure S1D-S1E). To explore the clinical significance of AL355338, we analyzed RNA-seq data from the publicly available TCGA dataset (StarBase V3.0 platform). AL355338 expression was significantly upregulated in lung adenocarcinoma and lung squamous carcinoma tissues (Figure S2A-S2B). Kaplan–Meier survival analysis revealed that patients with high AL355338 expression has shorter overall survival than patients with low AL355338 expression (Figure S2C-S2D).

QRT-PCR analysis confirmed that AL355338 was increased in tumor tissues compared with the corresponding paratumor tissues in a panel of 80 specimens from NSCLC patients (Figure 1A). Furthermore, AL355338 expression was remarkably higher in metastatic samples compared with non-metastatic samples, and elevated AL355338 expression correlated with advanced TNM stage (Figure 1B-1C). Kaplan–Meier survival analysis of overall survival showed that high AL355338 expression was associated with shorter survival (Figure 1D; P = 0.032). We next detected AL355338 with FISH and evaluated the relationship between its expression and clinicopathological characteristics. Consistent with the qRT-PCR results, more positive AL355338 expression signals were observed in tumor tissues compared to adjacent normal tissues (Figure 1E). Expression score analyses demonstrated that expression percentages of AL355338 were higher in stage III + IV than in stage I + II NSCLC samples (Figure 1F). High AL355338 expression was positively associated with T stage, N stage, distant metastasis, and TNM stage (Table 1). However, it was not related to patient age, sex, or tumor location. Uni- and multivariate regression analyses revealed that AL355338 expression was an independent predictor of NSCLC aggressiveness and could predict clinical outcome (Figure 1G-1H). Collectively, these results suggested that AL355338 was upregulated in NSCLC and might be a useful prognostic biomarker.

AL355338 promoted NSCLC cell proliferation in vitro

To explore the potential function of AL355338 on cellular behavior, we first detected the expression levels of AL355338 in NSCLC cell lines. Its expression in NSCLC cell lines was significantly higher than that in the normal lung epithelial cell line HBE (Figure S3A). We then examined its subcellular localization in NSCLC cells. Cellular FISH experiments and analyses of RNA isolated from the cytoplasmic and nuclear compartments of NSCLC cells showed that AL355338 was mainly present in the cytoplasm (Figure 2A-2B). Then, we downregulated AL355338 by transfecting two different siRNAs in A549 and PC9 cells. qRT-PCR was performed to verify the knockdown efficiency of AL355338; both si-AL355338#1 and si-AL355338#2 successfully downregulated AL355338 expression to 20-30% in A549 and PC9 cells (Figure 2C).

Flow cytometry analyses of cell cycle distribution indicated that AL355338 downregulation significantly decreased the number of cells in S phase, with an increase of cells in G0/G1 phase for both A549 and PC9 cells (Figure 2D-2F). CCK-8 and cell colony formation assays demonstrated that the proliferation and colony formation ability of A549 and PC9 cells were impaired after AL355338 knockdown (Figure 2E, 2G, 2H). Subsequent EdU assays indicated that DNA synthesis by A549 and PC9 cells was markedly reduced after silencing AL355338 (Figure 2I). These results indicated that AL355338 was effective in facilitating NSCLC cell proliferation.

AL355338 facilitated aggressive phenotypes and glycolytic metabolism of NSCLC cells

To further assess the migratory and invasive ability conferred by AL355338, we performed transwell and wound healing assays. The transwell results in Figure 3A-3C show that downregulation of AL355338 suppressed the migratory and invasive activities of A549 and PC9 cells compared to the negative control. Wound healing assays showed the same pattern (Figure 3D-3F). The initiation of epithelial-mesenchymal transition is a driving force of tumor invasion and metastasis. Analyses of the protein expression of epithelial and mesenchymal markers revealed that AL355338 silencing reduced N-cadherin and vimentin expression but enhanced E-cadherin expression in A549 and PC9 cells (Figure 3G). Consistently, IF staining also showed that the fluorescence intensity of Vimentin was attenuated in AL355338 downregulated A549 cells (Figure S3B).

Considering AL355338 expression was upregulated in association with high ¹⁸F-FDG uptake, we investigated the effect of AL355338 in glycolysis regulation. Glucose enzyme expression levels were examined, and the results showed that alteration of AL355338 expression could positively affect hexokinase II, pyruvate kinase isozymes M2, and lactate dehydrogenase A protein levels (Figure 3H). To confirm the effect of AL355338 on glycolysis, we measured ¹⁸F-FDG uptake, lactate release, and ATP production to reflect the flux distribution of glucose. AL355338 knockdown led to a decreases in all three parameters (Figure 3I-3K). These results demonstrated that AL355338 could regulate glycolytic metabolism in NSCLC.

Silencing AL355338 inhibited tumorigenesis, aerobic glycolysis, and metastasis in vivo

To determine whether AL355338 expression could affect NSCLC cell growth in vivo, A549 cells were stably transfected with sh-AL355338 or negative control sh-NC lentivirus. In vivo xenograft experiments showed that tumor volumes in the sh-AL355338 group were obviously smaller than in the sh-NC group (Figure 4A). Notably, tumor weights were markedly reduced in the sh-AL355338 group compared with controls (Figure 4B). In addition, xenograft glycolysis uptake ability was observed by ¹⁸F-FDG micro-PET/CT imaging. Tumor glycolysis was significantly blocked in sh-AL355338 xenografts, with lower SUVmax (mean SUVmax: 1.15) compared with sh-NC xenografts (mean SUVmax: 2.08; Figure 4C-4D). IHC revealed that the tumor tissue grown from sh-AL355338 cells displayed lower Ki-67, ENO1, and EGFR staining than those formed from sh-NC cells. FISH analyses confirmed the inhibitory effect of knocking down AL355338 expression in A549 cells (Figure 4E).

We then investigated the function of AL355338 on tumor metastasis in vivo. Micro-CT images showed more metastatic lesions in the lungs of mice in the sh-NC group, while fewer and smaller metastatic lung nodules were seen in the AL355338 knockdown group (Figure 4F). Summarized data on HE staining analysis demonstrated that more metastatic lung foci were formed in the sh-NC group compared to the AL355338 knockdown group (Figure 4G). These results indicated that inhibiting AL355338 expression in NSCLC cells could efficiently suppress tumor metastasis in vivo.

AL355338 specifically interacted with ENO1 and enhanced ENO1 protein stability

To investigate the underlying mechanism of AL355338-mediated glycolysis and carcinogenesis in NSCLC, we set out to identify intracellular AL355338-binding factors in an unbiased manner. Biotinylated AL355338 and negative control were incubated with total protein extracts from A549 cells and pulled down with streptavidin (Figure S4A). MS was conducted following RNA pull-down to screen endogenous AL355338 binding proteins (Table S3). From the MS data files, we selected ENO1 for further binding validation as it is a key enzyme responsible for the glycolytic pathway. Subsequent western blotting of proteins retrieved from AL355338 pull-down assays confirmed that the ENO1 protein directly bound to the sense AL355338 sequence but did not bind the antisense control (Figure 5A, upper). In addition, RIP assays with an ENO1 antibody and IgG antibody control were performed using A549 and PC9 cellular extracts. Correspondingly, we observed significantly higher AL355338 enrichment with the ENO1 antibody compared to IgG (Figure 5A, lower). Both assays demonstrated that AL355338 directly interacted with ENO1. To better identify which region of ENO1 interacted with AL355338, we performed RIP assays using FLAG-tagged truncated forms of ENO1. The results demonstrated that the ENO1-C domain (237–434 aa) was essential for AL355338 binding (Figure 5B). To map the ENO1-binding regions of AL355338, we constructed four AL355338 deletion fragments. Subsequent RNA pull-down assays followed by western blotting revealed that the ENO1 specific binding sequence was located between 528 and 913 nucleotides (Figure 5C). RNA FISH of AL355338 combined with IF of ENO1 confirmed cytoplasmic co-localization of AL355338 and ENO1 in A549 cells (Figure 5D).

A recent study demonstrated that lncRNAs interact with proteins through several modes of action, such as regulating protein expression and function, modulating protein–protein interactions, or guiding protein subcellular localization[21]. To investigate the molecular effect of AL355338 on ENO1, western blotting was performed to measure ENO1 expression after A549 and PC9 cells were subjected to different treatments with AL355338. The results showed that AL355338 positively regulated ENO1 protein expression (Figure 5E). Similar trends were observed in IF assays (Figure S4B). Interestingly, loss or gain of function of AL355338 had no influence on ENO1 mRNA expression (Figure S4C-S4D), suggesting that AL355338 may regulate ENO1 protein stability. Next, cycloheximide (CHX) chase assays were performed to observe the effect of AL355338 on ENO1 protein degradation. Western blotting showed that ENO1 protein was more stable in AL355338-transfected cells, while it had a shorter half-life in AL355338-silenced cells (Figure 5F-5G). Moreover, the proteasome inhibitor MG132 rescued the reduction of ENO1 caused by AL355338 repression (Figure 5H), suggesting that AL355338 elevated ENO1 levels by reducing its degradation. As reported, cytoplasmic protein degradation involves the ubiquitin-proteasome pathway[22]. Thus, whether AL355338 was involved in the ubiquitin-mediated ENO1 degradation was tested with in vitro ubiquitination assays. The results illustrated that upregulation of AL355338 decreased ENO1 ubiquitination, whereas silencing of AL355338 enhanced ENO1 ubiquitination (Figure 5I). Collectively, our data indicated that AL355338 directly bound the ENO1 protein and enhanced its stability in NSCLC cells.

ENO1 participated in the oncogenic function of AL355338

To clarify whether the RNA–protein interaction between AL355338 and ENO1 mediated the key oncogenic function of AL355338 in NSCLC, complementary rescue experiments were performed. Stable ENO1-knockdown A549 and PC9 cells were constructed and then transfected with AL355338 overexpression plasmid. Transfection efficiency was confirmed by qRT-PCR and western blot (Figure 6A-6B). As was expected, CCK-8 and cell colony formation assays showed that ectopic AL355338 expression increased the proliferation activity and colony formation ability of A549 and PC9 cells, while ENO1-knockdown partially abrogated the effects induced by AL355338 upregulation (Figure 6C-6D). EdU assays revealed that increased DNA synthesis capability of A549 and PC9 cells induced with AL355338 was reversed by ENO1 silencing (Figure 6E). Cell migration and invasion abilities were also assessed, and results demonstrated that AL355338 exerted minimal effects on cell migration and invasion in the absence of ENO1 (Figure 6F). Importantly, the stimulatory effect of AL355338 on epithelial and mesenchymal marker expression was recovered by repressing ENO1 (Figure 6G). In glycolysis metabolism assays, ENO1 downregulation dramatically blocked the AL355338-induced increase in ¹⁸F-FDG uptake, as well as lactate production by A549 and PC9 cells (Figure 6H-6I). The above results suggested that AL355338 regulation of glycolytic metabolism and cell proliferation were dependent on ENO1 stability in NSCLC cells.

AL355338 activated the EGFR/AKT pathway by targeting the ENO1/EGFR interaction

The EGFR/AKT signaling pathway inﬂuences tumor cell survival in many cancers. As previous reports have demonstrated that ENO1 could regulate the AMPK-mediated or FAK-mediated AKT pathway[23, 24], we sought to determine whether ENO1 could interact with EGFR and subsequently activate its downstream AKT kinase signaling pathway. To verify the relationship between ENO1 and EGFR in NSCLC cells, we first examined their TCGA database expression levels. We extracted ENO1 and EGFR data from 34 tumors and found that both were notably upregulated in diverse cancer types, especially lung adenocarcinoma and lung squamous carcinoma tissues (Figure S5A-S5B). As shown in Figure S5C, the correlation coefficient of ENO1 and EGFR reached 0.32 in the NSCLC cohort (P < 0.001, Spearman correlation analysis). Besides, ENO1 shared a similar mutation status with EGFR in NSCLC samples (P<0.001, Figure S5D).

Next, we aimed to investigate whether ENO1 directly interacted with EGFR. Reciprocal co-IP assays were performed using both exogenous labeled antibodies and endogenous antibodies to precipitate ENO1 and EGFR, respectively. The co-IP results confirmed physical binding between ENO1 and EGFR proteins in NSCLC cells (Figure 7A-7B). IF staining confirmed ENO1 and EGFR co-localization in A549 and PC9 cells (Figure 7C). Furthermore, in situ PLA was used to visualize native protein complexes. High densities of PLA-positive protein complexes that included ENO1 and EGFR were observed in A549 and PC9 cells, which was in accordance with the co-IP results (Figure 7D). Intriguingly, AL355338 knockdown noticeably impaired the ENO1 and EGFR interaction, while AL355338 overexpression had the opposite effect (Figure 7E-7F).

To understand the potential relationship between AL355338 and the EGFR/AKT pathway, we investigated EGFR and AKT expression. Western blot analyses showed that AL355338 suppression decreased the levels of phosphorylated EGFR and AKT, but total protein levels were not affected. Overexpression of AL355338 enhanced EGFR and AKT phosphorylation, indicating that the EGFR/AKT pathway was a downstream effector of AL355338. We then assessed the ability of ENO1 to overcome AL355338-mediated activation of the EGFR/AKT pathway in NSCLC cells. Notably, ENO1 knockdown abrogated AL355338-mediated increases of EGFR and AKT phosphorylation, whereas ENO1 overexpression reversed sh-AL355338-mediated inactivation of the EGFR/AKT pathway (Figure 7G-7H). Based on these results, we inferred that AL355338 could regulate EGFR/AKT pathway activity by targeting ENO1/EGFR interaction.

ENO1 was upregulated in NSCLC and positively correlated with AL355338 expression

To further establish the relationship between AL355338 and ENO1 in NSCLC tissues, we first evaluated ENO1 protein expression with IHC in 80 NSCLC and pair-matched normal tissue samples (Figure 8A). The results revealed that 76.3% (61/80) NSCLC samples showed increased ENO1 expression compared with their pair-matched normal tissues (Figure 8B). ENO1 protein levels positively correlated with AL355338 levels (r = 0.382, P < 0.001, Figure 8C), suggesting potential regulation of ENO1 by AL355338 in clinical NSCLC tissues. Meanwhile, TCGA datasets revealed that ENO1 expression was significantly upregulated in lung adenocarcinoma and lung squamous carcinoma tissues compared to normal tissues (Figure 8D-8E). Kaplan-Meier analysis of the 80-patient NSCLC cohort showed that high ENO1 expression was associated with poor prognosis (P = 0.009, Figure 8F). This was validated in the TCGA cohort; survival analyses showed that patients with low ENO1 expression had longer overall survival (P = 0.025) and disease-free survival (P = 0.044) compared to those with high ENO1 expression (Figure 8G-8H).

Despite developments in therapeutic approaches, the prognosis for NSCLC patients diagnosed in an advanced stage remains poor[25]. the enhanced Warburg effect is a leading cause of malignant features in NSCLC. Multiple glycolysis-associated genes have been reported as contributors to cancer progression[26]. Earlier studies indicated that numerous lncRNAs were associated with NSCLC initiation and development of NSCLC[7, 27]. Since ¹⁸F-FDG uptake can reflect the Warburg effect, we previously subjected high ¹⁸F-FDG-uptake NSCLC tissues to RNA-seq analysis to identify glycolysis-related lncRNAs in our previous study[16]. AL355338 was a newly identified oncogenic lncRNA among the top-scoring upregulated genes. Clinical analysis showed that AL355338 expression was notably higher in tumor tissues than in corresponding adjacent tissues. Importantly, high AL355338 expression in patients with NSCLC was related to advanced TNM stage and poor clinical outcome. Functionally, the positive regulation of AL355338 in promoting NSCLC proliferation and metastasis was identified in a series of in vitro and in vivo experiments. The observation of high AL355338 expression in glycolytically active tissues led us to investigate whether AL355338 could regulate cancer metabolism. In support of our hypothesis, downregulation of AL355338 dramatically weakened the Warburg effect in NSCLC cells. The results indicate an oncogenic role of AL355338 in regulating malignant phenotype, since abnormal cancer cell metabolism is likely to induce cancer aggressiveness.

In addition to its biological importance, we dissected the mechanisms by which AL355338 mediated glycolysis pathway activation. A cross-linking RNA pull-down/MS strategy was used to detect proteins that directly interacted with AL355338. Among the MS data, the glycolytic enzyme ENO1 piqued our interest. As a well-documented glycolytic enzyme responsible for ATP generation, ENO1 is related to mediate energy metabolism and cancer development[11]. Recent studies reported that ENO1 was correlated with cancer differentiation and progression in several tumor types by interacting with lncRNAs. Hong et al[12] concluded that lncRNA ENO1-IT1, which was activated by Fusobacterium nucleatum abundance, could promote glycolysis and oncogenesis by targeting ENO1 and guiding its histone modification pattern in colorectal cancer. Another group confirmed that lncRNA SNHG18 promoted glioma cell motility by disrupting the nucleocytoplasmic shuttling of ENO1[28]. Consistent with the previous report, our results suggested that ENO1 was responsible for the oncogenic role of AL355338 in mediating NSCLC progression. ENO1 silencing weakened the exogenous AL355338-induced metabolic shift in NSCLC cells. AL355338 also maintained ENO1 protein stabilization through decreasing ubiquitin-mediated degradation. Notably, we observed simultaneous upregulation and correlation of AL355338 and ENO1 in NSCLC tissues. ENO1 and AL355338 overexpression were strongly associated with poor prognosis in NSCLC patients, suggesting that AL355338 and ENO1 levels may predict patient outcome.

In addition to its canonical roles in glycolysis, accumulating evidence revealed ENO1 as a multifunctional protein involved in regulating several signaling pathways[29]. A recent study reported that ENO1 promoted the self-renewal and malignant phenotype of lung cancer stem cells by affecting the AMPK/mTOR pathway[23]. Fu et al[24] showed that ENO1 was involved in NSCLC glycolysis, proliferation, and metastasis through the FAK-mediated PI3K/AKT pathway. It is rather remarkable that EGFR is the best-characterized receptor tyrosine kinase and EGFR signaling is among the most dysregulated molecular pathways in NSCLC[30]. Upon ligand binding to EGFR, the AKT signaling pathway will be activated, leading to enhanced glycolysis and oncogenic phenotypes in tumor cells[31]. Given the positive correlation of ENO1 and EGFR expression in NSCLC and their similar functions in AKT pathway regulation, these reports inspired us to explore whether ENO1 might be an interaction partner of EGFR. Notably, a direct interaction between ENO1 and EGFR was observed by reverse co-IP, and native ENO1/EGFR protein complexes were visualized with PLA assays. Moreover, our data clearly showed that ENO1 positively regulated EGFR phosphorylation and thus activated downstream AKT kinase signaling.

Previous studies revealed that lncRNAs might act as a tumorigenic regulator by modulating signaling networks[32]. Apart from directly modulating ENO1 expression, we hypothesized that AL355338 functions through the ENO1-mediated EGFR/AKT pathway in NSCLC cells. As expected, ectopic AL355338 expression could enhance ENO1/EGFR interaction, which subsequently led to higher levels of phosphorylated EGFR and AKT. Opposite results were observed in AL355338-depleted cells. However, ENO1 silencing attenuated the increases in phosphorylation modulated by AL355338 induction. ENO1 overexpression restored the reduction in phosphorylated EGFR and AKT levels induced by AL355338 knockdown. Collectively, these results indicated that AL355338 might regulate the ENO1-mediated EGFR/AKT pathway.

Our study mainly focused on the glycolysis-associated lncRNA AL355338 in NSCLC. Remarkably, AL355338 was increased in NSCLC tissues, and its levels correlated with clinical outcomes, suggesting it could be an independent biomarker in patients with NSCLC. Further functional assays suggested that AL355338 functions as an oncogene by mediating glycolysis activation and promoting NSCLC progression through selective targeting of ENO1 protein. Additionally, we demonstrated that AL355338 stabilized ENO1 protein to prevent ubiquitin-mediated degradation. AL355338 was also capable of modulating interaction between ENO1 and EGFR, further activating EGFR-AKT signaling. In conclusion, our results suggested that upregulated lncRNA AL355338 might serve as a potential biomarker and therapeutic target in the clinical treatment of NSCLC.

lncRNA: long noncoding RNA; NSCLC: non-small cell lung cancer; EMT: epithelial-mesenchymal transition; FDG: fluorodeoxyglucose; TNM: tumor-node-metastasis staging system; IP, immunoprecipitation; RIP: RNA immunoprecipitation; FISH: fluorescence in situ hybridization; IHC: immunohistochemistry; IF: Immunofluorescence; PLA: Proximity Ligation Assay; ENO1: enolase 1; HK2: hexokinase II; PKM2: Pyruvate kinase isozymes M2; LDHA: lactate dehydrogenase A; EGFR: epidermal growth factor receptor; AKT: serine/threonine kinase; TCGA: the cancer genome atlas.

Acknowledgments

Not applicable.

Authors’ contributions

Baoming Mi and Zengli Liu conceived and designed the study. Qian Hua and Dongliang Wang performed the experiments and collected data; Lin Zhao provided technical support and performed the data analysis; Zhihui Hong, Kairu Ni, Yizhen Shi provided and collected the clinical data. Qian Hua drafted the manuscript. Baoming Mi supervised the overall research and revision of the manuscript. All authors read and approved the final manuscript.

Funding

Our study was supported by research grants from the National Natural Science Foundation of China (Grant No. 81703159, 82073482), the Natural Science Foundation of Jiangsu Province (Grant No. 81830052), the Postdoctoral research funding plan in Jiangsu province (1601121C)，the Suzhou Science and Technology Program (Grant No. KJXW2020016), the doctoral research funding project of the Second Affiliated Hospital of Soochow University (Grant No. SDFEYBS2012, SDFEYBS2013, XKTJ-HRC2021003).

Availability of data and materials

All data generated or analyzed during this study are included in this published article.

Declarations

Ethics approval and consent to participate

This study was reviewed and approved by the Ethics Committee of the Second Affiliated Hospital of Soochow University. Written informed consent was obtained from each patient. And this study was performed in accordance with the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflict of interest.

Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer, Statistics. 2021. CA Cancer J Clin. 2021;71:7–33.
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71:209 – 49.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144:646 – 74.
Ntziachristos V, Pleitez MA, Aime S, Brindle KM. Emerging Technologies to Image Tissue Metabolism. Cell Metab. 2019;29:518 – 38.
Torresano L, Nuevo-Tapioles C, Santacatterina F, Cuezva JM. Metabolic reprogramming and disease progression in cancer patients. Biochim Biophys Acta Mol Basis Dis. 2020;1866:165721.
Martinez-Outschoorn UE, Peiris-Pagés M, Pestell RG, Sotgia F, Lisanti MP. Cancer metabolism: a therapeutic perspective. Nat Rev Clin Oncol. 2017;14:113.
Anastasiadou E, Jacob LS, Slack FJ. Non-coding RNA networks in cancer. Nat Rev Cancer. 2018;18:5–18.
Tsagakis I, Douka K, Birds I, Aspden JL. Long non-coding RNAs in development and disease: conservation to mechanisms. J Pathol. 2020;250:480 – 95.
Xu Y, Qiu M, Shen M, Dong S, Ye G, Shi X, et al. The emerging regulatory roles of long non-coding RNAs implicated in cancer metabolism. Mol Ther. 2021;29:2209-18.
Kim JW, Dang CV. Multifaceted roles of glycolytic enzymes. Trends Biochem Sci. 2005;30:142 – 50.
Almaguel FA, Sanchez TW, Ortiz-Hernandez GL, Casiano CA. Alpha-Enolase: Emerging Tumor-Associated Antigen, Cancer Biomarker, and Oncotherapeutic Target. Front Genet. 2020;11:614726.
Hong J, Guo F, Lu SY, Shen C, Ma D, Zhang X, et al. F. nucleatum targets lncRNA ENO1-IT1 to promote glycolysis and oncogenesis in colorectal cancer. Gut. 2020.
Yu S, Li N, Huang Z, Chen R, Yi P, Kang R, et al. A novel lncRNA, TCONS_00006195, represses hepatocellular carcinoma progression by inhibiting enzymatic activity of ENO1. Cell Death Dis. 2018;9:1184.
Detterbeck FC, Boffa DJ, Kim AW, Tanoue LT. The Eighth Edition Lung Cancer Stage Classification. Chest. 2017;151:193–203.
Tang J, Yan T, Bao Y, Shen C, Yu C, Zhu X, et al. LncRNA GLCC1 promotes colorectal carcinogenesis and glucose metabolism by stabilizing c-Myc. Nat Commun. 2019;10:3499.
Hua Q, Jin M, Mi B, Xu F, Li T, Zhao L, et al. LINC01123, a c-Myc-activated long non-coding RNA, promotes proliferation and aerobic glycolysis of non-small cell lung cancer through miR-199a-5p/c-Myc axis. J Hematol Oncol. 2019;12:91.
Shen M, Zhao X, Zhao L, Shi L, An S, Huang G, et al. Met is involved in TIGAR-regulated metastasis of non-small-cell lung cancer. Mol Cancer. 2018;17:88.
Jia W, Zhao X, Zhao L, Yan H, Li J, Yang H, et al. Non-canonical roles of PFKFB3 in regulation of cell cycle through binding to CDK4. Oncogene. 2018;37:1685-98.
Hua Q, Mi B, Xu F, Wen J, Zhao L, Liu J, et al. Hypoxia-induced lncRNA-AC020978 promotes proliferation and glycolytic metabolism of non-small cell lung cancer by regulating PKM2/HIF-1α axis. Theranostics. 2020;10:4762-78.
Li L, Liang Y, Kang L, Liu Y, Gao S, Chen S, et al. Transcriptional Regulation of the Warburg Effect in Cancer by SIX1. Cancer Cell. 2018;33:368 – 85.e7.
Díaz-Ramos A, Roig-Borrellas A, García-Melero A, López-Alemany R α-Enolase, a multifunctional protein: its role on pathophysiological situations. J Biomed Biotechnol. 2012;2012:156795.
Wang J, Maldonado MA. The ubiquitin-proteasome system and its role in inflammatory and autoimmune diseases. Cell Mol Immunol. 2006;3:255 – 61.
Dai J, Zhou Q, Chen J, Rexius-Hall ML, Rehman J, Zhou G. Alpha-enolase regulates the malignant phenotype of pulmonary artery smooth muscle cells via the AMPK-Akt pathway. Nat Commun. 2018;9:3850.
Fu QF, Liu Y, Fan Y, Hua SN, Qu HY, Dong SW, et al. Alpha-enolase promotes cell glycolysis, growth, migration, and invasion in non-small cell lung cancer through FAK-mediated PI3K/AKT pathway. J Hematol Oncol. 2015;8:22.
Ramalingam SS, Owonikoko TK, Khuri FR. Lung cancer: New biological insights and recent therapeutic advances. CA Cancer J Clin. 2011;61:91–112.
Abdel-Wahab AF, Mahmoud W, Al-Harizy RM. Targeting glucose metabolism to suppress cancer progression: prospective of anti-glycolytic cancer therapy. Pharmacol Res. 2019;150:104511.
Vencken SF, Greene CM, McKiernan PJ. Non-coding RNA as lung disease biomarkers. Thorax. 2015;70:501-3.
Zheng R, Yao Q, Li X, Xu B. Long Noncoding Ribonucleic Acid SNHG18 Promotes Glioma Cell Motility via Disruption of α-Enolase Nucleocytoplasmic Transport. Front Genet. 2019;10:1140.
Cappello P, Principe M, Bulfamante S, Novelli F. Alpha-Enolase (ENO1), a potential target in novel immunotherapies. Front Biosci (Landmark Ed). 2017;22:944 – 59.
Stencel K, Chmielewska I, Milanowski J, Ramlau R Non-Small-Cell Lung Cancer: New Rare Targets-New Targeted Therapies-State of The Art and Future Directions. Cancers (Basel). 2021;13.
Zaryouh H, De Pauw I, Baysal H, Peeters M, Vermorken JB, Lardon F, et al. Recent insights in the PI3K/Akt pathway as a promising therapeutic target in combination with EGFR-targeting agents to treat head and neck squamous cell carcinoma. Med Res Rev. 2021.
Beermann J, Piccoli MT, Viereck J, Thum T. Non-coding RNAs in Development and Disease: Background, Mechanisms, and Therapeutic Approaches. Physiol Rev. 2016;96:1297 – 325.

Table 1. Correlation of the expression of AL355338 in NSCLC with clinicopathologic features

Characteristics		AL355338		P value
Characteristics		Low (n)	High (n)
Age (years)
<60	34 (42.5%)	16	18	0.54
≥60	46 (57.5%)	21	25
Gender
Male	41 (51.2%)	17	24	0.26
Female	39 (48.8%)	20	19
TNM stage
I	23 (28.7%)	19	4	<0.01**
II	19 (23.8%)	10	9
III	22 (27.5%)	4	18
IV	16 (20.0%)	0	16
T stage
T1	21 (26.3%)	12	9	0.04*
T2	37 (46.3%)	19	18
T3	16 (20.0%)	5	11
T4	6 (7.5%)	1	5
N stage
N0	32 (40.0%)	24	8	0.03*
N1	17 (21.3%)	7	10
N2	23 (28.7%)	6	17
N3	8 (10.0%)	1	7
M stage
M0	65 (81.3%)	37	28	<0.01**
M1	15 (18.8%)	0	15
Mortality
Survive	30 (37.5%)	17	13	0.04*
die	50 (62.5%)	20	30
Histological grade
I	5 (6.2%)	2	3	0.38
II	42 (52.5%)	17	25
III	33 (41.3%)	11	22
Tumor location
left	35 (43.8%)	16	19	0.55
right	45 (56.2%)	22	23

*P<0.05; **P<0.01

Supplementarymaterials.pdf

AL355338 acts as an oncogenic lncRNA by interacting with protein ENO1 to regulate EGFR/AKT pathway in NSCLC

Status:

Version 1

Abstract

Background

Methods

Results

Conclusions

Figures

Background

Materials And Methods

Results

Discussion

Conclusions

Abbreviations

Declarations

References

Tables

Supplementary Files

Status:

Version 1