Skip to main content
Download PDF

Altered expression of ACOX2 in non-small cell lung cancer

BMC Pulmonary Medicine volume 22, Article number: 321 (2022) Cite this article

Abstract

Peroxisomes are organelles that play essential roles in many metabolic processes, but also play roles in innate immunity, signal transduction, aging and cancer. One of the main functions of peroxisomes is the processing of very-long chain fatty acids into metabolites that can be directed to the mitochondria. One key family of enzymes in this process are the peroxisomal acyl-CoA oxidases (ACOX1, ACOX2 and ACOX3), the expression of which has been shown to be dysregulated in some cancers. Very little is however known about the expression of this family of oxidases in non-small cell lung cancer (NSCLC). ACOX2 has however been suggested to be elevated at the mRNA level in over 10% of NSCLC, and in the present study using both standard and bioinformatics approaches we show that expression of ACOX2 is significantly altered in NSCLC. ACOX2 mRNA expression is linked to a number of mutated genes, and associations between ACOX2 expression and tumour mutational burden and immune cell infiltration were explored. Links between ACOX2 expression and candidate therapies for oncogenic driver mutations such as KRAS were also identified. Furthermore, levels of acyl-CoA oxidases and other associated peroxisomal genes were explored to identify further links between the peroxisomal pathway and NSCLC. The results of this biomarker driven study suggest that ACOX2 may have potential clinical utility in the diagnosis, prognosis and stratification of patients into various therapeutically targetable options.

Peer Review reports

Introduction

Lung cancer remains the most common diagnosed malignancy in the world and the most recent Globocan analysis indicates that there will be an estimated 19.3 million new cancer cases worldwide, of which 11.4% will be lung cancer and at 18% the leading cause of cancer death (with an estimated 1.8 million deaths) [1, 3]. Lung cancer itself can be separated into two subtypes small cell lung cancer (SCLC) and Non-Small Cell Lung Cancer (NSCLC), with the vast majority (approximately 85%) falling into the latter subtype. Treatment options for NSCLC are still limited, but recent advances in immunotherapy [4], targeted therapy [5] and the discovery of new actionable mutations [6, 7] have greatly increased the treatment options available for this cancer.

Since their discovery in 1954 [8], peroxisomes have emerged as key metabolic organelles with many diverse functions ranging from cellular lipid metabolism and reactive oxygen species, to non-metabolic roles such as cellular stress responses and synthesis of cellular signalling molecules [9, 10]. These are achieved through interactions with other cellular compartments such as the endoplasmic reticulum and mitochondria [10].

Evidence linking peroxisomes to the biology of NSCLC would stem strongly from the earlier studies demonstrating the roles of peroxisome proliferator-activated receptors (PPARs) in lung cancer tumorigenesis [11]. Moreover, peroxisomes have been linked to cancer via their roles in aberrant metabolism and crosstalk with mitochondria (and associated mitochondrial dysfunction) in cancer [12]. In lung cancer aberrant metabolism is well documented [13, 14], as is mitochondrial dysfunction [15, 16].

One of the main functions of peroxisomes is the processing of very-long chain fatty acids (VLCFAs) via either α- or β- oxidation into metabolites that can be directed to the mitochondria [12]. Peroxisomal acyl-CoA oxidases (ACOs) have been described as the peroxisomal equivalent of the mitochondrial acyl-CoA dehydrogenases (ACADs) [17, 18]. Belonging to the same superfamily, ACOs contain one non-covalently bound Flavin adenine dinucleotide (FAD) per subunit [17, 18]. Similar to mitochondrial ACADs, ACOs catalyse the formation of α-, β- dehydrogenation of acyl-CoA, the initial and rate-determining step of the peroxisomal fatty acid β-oxidation pathway [17, 18]. Depending on the type of very long chain fatty acid (VLCFA) and/or its linear or branched status, this step is carried out by either ACOX1, ACOX2, or ACOX3 [12]. Oxidation of straight-chain fatty acids with different chain lengths is conducted by ACOX1 [12, 19]. ACOX2 and ACOX3 have known roles in the degradation of branched-chain fatty acids [19]. ACOX2 however is the only known ACO associated with bile acid biosynthesis in humans [19].

As such these diverse roles indicate that peroxisomes by their nature play important roles in a large number of globally important human disease including obesity, cancer and age-related disorders [9]. In the disease setting peroxisomes were originally associated with various metabolic disorders [20, 21], and more recently various cancers have been identified as having dysregulation of peroxisomal genes/proteins including prostate, breast and lung cancer [12, 22,23,24,25]. In particular alterations to the β-oxidation of peroxisomal fatty acids have been extensively studied in breast cancer [12], and whilst many of the results come from individual stand-alone studies, the upregulation of ACOX1 and other members of this process in cancer cells suggests that the entire β-oxidation metabolic pathway may be affected and linking the entire β-oxidation process to breast cancer tumorigenesis [12].

This pathway has not been studied in such detail in NSCLC, although evidence is beginning to emerge linking the β-oxidation of peroxisomal fatty acids to NSCLC. In the analysis of the peroxisomal pathway in NSCLC by Zhang et al. [25], 38 differentially expressed genes were identified from analyses of The Cancer Genome Atlas (TCGA) NSCLC datasets, but none of the identified genes belonged to the ACOX family. Surprisingly ACOX2 was described as a gene commonly upregulated across several TCGA cancer datasets including NSCLC [12].

As ACOX2 was suggested to be elevated at the mRNA level in over 10% of NSCLC [12], we undertook an analysis of this gene in NSCLC to determine if it had any potential utility as a biomarker in lung cancer at both mRNA and protein level. At the same time using in silico analyses we explored all of the peroxisomal ACOX family members in NSCLC. We show that in contrast to the initial data [12] ACOX2 is predominantly downregulated in NSCLC at both the mRNA and protein level. Decreased ACOX2 expression may be associated with hyper-methylation of individual CpG residues in the ACOX2 promoter. Expression of ACOX2 is associated with prognosis (OS) predominantly in the lung adenocarcinoma (LUAD) subset. We further link altered expression of the other members of both the ACOX family and other peroxisomal proteins that associate with these enzymes to both aberrant expression and prognostic value in NSCLC, and altered ACOX2 expression may be associated sensitivity to certain drugs [12].

Materials and methods

Primary tumour samples

Twenty-two surgically resected chemotherapy naïve fresh-frozen tumour specimens were used in this study, taken from surgical resections of at St. James’s Hospital, Dublin during the period 2011–2016. Immediately following resection all samples were evaluated by a pathologist and tumour tissue along with matched normal tissue were dissected for downstream analysis. In total, 11 adenocarcinomas and 11 squamous cell carcinomas were utilized in this study and a summary of the histopathological and clinicopathological data is presented in Table 1. Informed consent was obtained from each patient, and the study was conducted after formal approval from the SJH/AMNCH Hospital Ethics Committee—Ethics REC (No.: 041018/8804).

Table 1 Details of surgically resected fresh frozen patient samples used in this study
Full size table

Fixed formalin paraffin embedded samples

A total of 204 surgically resected NSCLC tumour specimens from the period 1999–2007 were included in this study. All surgically resected tumour specimens and control specimens were fixed with 10% formalin and embedded in paraffin (FFPE). The Union for International Cancer Control Tumour-Node-Metastasis (TNM) Classification of Malignant Tumours 8th edition was used to stage the tumours [26, 27] and histologically subtyped using the World Health Organization guidelines [28, 29]. A summary of the histopathological and clinicopathological data (including age, sex, smoking status, histology, TNM stage, surgical procedure, tumour grade, and primary site) for the cohort of patients utilized are presented in Table 2.

Table 2 Patient characteristics in the SJH NSCLC TMA
Full size table

A Beecher Manual Tissue Arrayer (Model MTA-1) was used to generate of a tissue microarray (TMA) containing quadruplicate cores (0.6 mm) of the FFPE embedded samples and a 4 µm sections were subsequently used for immunohistochemistry (IHC) of ACOX2.

Immunohistochemistry

IHC was performed on TMA sections by utilizing a standard protocol to deparaffinise, rehydrated and wash the slides. Subsequently, ULTRA cell conditioning (ULTRA CC1), pH9.1, was used for heat induced epitope retrieval (HIER). For ACOX2 antibody staining, slides were incubated with rabbit polyclonal primary antibody HPA038280 (Atlas Antibodies, Merck) diluted in ventana antibody diluent (Ref: 251–018) (1:50) for 32 min at ambient temperature and stained using the OptiView™ DAB IHC detection kit on a Ventana BenchMark XT processor.

Following IHC, staining was independently assessed by two pathologists blinded to the clinical, pathological and follow-up data. Staining intensity was designated as either 0, 1 + , 2 + or 3 + and each tumour section was given a H-score between 0 and 300 = 3(% at 3 +) + 2(% at 2 +) + 1(% at 1 +). No samples were observed with a H-Score of 300.

Tumors with high ACOX2 expression were designated as those with an average H score above the median value and low expression below the median. Kaplan–Meier analyses were constructed using Prism 5.01 (GraphPad).

RNA isolation and qPCR amplification

Total RNA was isolated and converted to cDNA using our previously described methodology [30,31,32]. Briefly, total RNA was extracted using TRI reagent (Molecular Research Center, Montgomery Road, OH, USA) according to the manufacturer’s instructions. 250 ng of this total RNA was then pre-treated to remove contaminating genomic DNA with amplification grade DNase I (Sigma-Aldrich, St. Louis, MO, USA) according to the manufacturer’s instructions [30,31,32]. All-in-One cDNA Synthesis SuperMix (Bimake, W Sylvanfield Drive, Houston, TX, USA) was then used to generate first strand cDNA according to the manufacturer’s instructions [30,31,32].

qPCRs for ACOX2 (relative quantification method) were subsequently conducted on these samples using an Illumina Eco qPCR and 2 × SYBR Green qPCR Master Mix (Bimake, W Sylvanfield Drive, Houston, TX, USA) using the manufacturers protocol in a 2-step qPCR program with histone H3 as the internal control using the following primers:

ACOX2 FWD: 5’– CCAAGTGGACATGGCAAGAA – 3’

ACOX2 REV: 5’ – GTGACTTCTGAGCCCACTGGA – 3’

Histone H3 FWD: 5’ – GGTAAAGCACCCAGGAAGCA – 3’.

Histone H3 REV: 5’ – CCTCCAGTAGAGGGCGCAC– 3’.

The following cycling parameters were used:

An initial Polymerase activation of 95 °C for 2 min followed by 35 cycles of 95 °C 15 s and annealing/amplification 61 °C for 1 min. A melting curve analysis was conducted at the end of each PCR using 95 °C 15 s, 55 °C 15 s and a final 95 °C for 15 s. The final concentration of primer in each reaction was 200 nM, and the.

Data was analysed using either the default in-built Eco software or imported into EcoStudy (Illumina) for analysis using the ΔΔCq method [33].

Identification of differential expressed genes for acyl-CoA oxidases in NSCLC

To determine if acyl-CoA oxidases show significantly altered expression in NSCLC a systematic analysis was carried out using Lung Cancer Explorer (LCE) [34] to conduct a meta-analysis of standardized mean difference (tumour vs normal) in gene expression, and subsequently a meta-analysis of survival and gene expression association based on univariate Cox Proportional-Hazards Model.

Establishment of protein–protein interaction networks associated with acyl-CoA oxidases

The STRING database [35] was interrogated to conduct a functional enrichment analysis of acyl-CoA oxidases to identify Protein–Protein Interaction Networks, exported to Cytoscape [36] and first neighbours highlighted. All identified candidates were subsequently interrogated for both differential gene expression and survival associations using LCE [34].

Validations of the expression differences of the differentially expressed genes associated with acyl-CoA oxidases

Validations of altered ACOX2 expression were examined in the TCGA NSCLC cohorts using LCE. Meta—analyses of additional gene expression datasets were also conducted using LCE [34]. Tumour—Normal Gene expression differences were calculated by LCE by meta-analysis of standardized mean difference (tumour—normal) using Hedges’ G as an effect size metric. For this all studies included for meta-analysis must have at least 10 samples in each group and meta-analysis was only performed for genes with data available from at least three qualifying studies. The analysis was performed separately for lung adenocarcinoma and squamous cell carcinoma [34]. Survival association with gene expression involved meta-analysis of survival and gene expression associations based on univariate Cox Proportional-Hazards modelling. For each included study, gene expression was normalized to zero mean and unit variance for each gene prior to model fitting. Only studies that had survival data for at least 10 samples were included in the meta-analysis which was only performed for genes if data was available from at least three qualifying studies [34].

Altered protein expression in the TCGA-LUAD was validated using UALCAN [37] to interrogate the CPTAC discovery dataset. This dataset comprises 111 tumors, (with 102 tumors paired with normal adjacent tissue samples) which have global proteome and phospho-proteome data available for interrogation [38].

Prognostic value for acyl-CoA oxidases and other associated differentially expressed genes

To examine for prognostic value, five separate analyses were conducted using either KM-Plot (Microarray and RNA-Seq) [39], LCE [34], UALCAN [37] and OncoLnc [40]. LCE, UALCAN and OncoLnc all utilise TCGA-LUAD and TCGA-LUSC datasets. All analyses used the default settings to assess prognostic survival as follows: KM-Plot (median); LCE (median); UALCAN (High vs Medium/Low); OncoLnc (50:50 percentiles).

All analyses use the Cox proportional hazards model.

Effects of copy number variations (CNV) on acyl-CoA oxidases

To examine whether the altered expressions of acyl-CoA oxidases (ACOX1-3) were correlated with CNVs, the datasets for TCGA-LUAD and TCGA-LUSC were analysed through cBioPortal [41, 42].

Associations with mutated genes

To identify mutated genes from the associated TCGA datasets that result in changes of expression of USO1 in NSCLC an analysis was conducted using the MuTarget platform with ACOX2 set as the target gene, and with mutation prevalence set at 2% [43]. All targets identified were subsequently retested using correlation analysis on TIMER2.0 [44].

Methylation analysis on acyl-CoA oxidases

Correlations between DNA methylation and altered ACOX2 expression were examined initially using UALCAN [37], and subsequently with MEXPRESS [45] and WANDERER [46] using default settings.

Correlations between ACOX2 expression and immune infiltrations in NSCLC

To interrogate how altered expression of ACOX2 expression was associated with tumour infiltrating immune cells (TIICs) in NCSLC we used TIMER [47] or TIMER2 [44]. The gene expression levels against tumour purity were calculated, and assessments were calculated based on the purity-corrected partial Spearman’s rho value with associated statistical significance. For analysis in GEPIA2 [48], a similar strategy was employed correlating ACOX2 expression with a six gene T-Cell exhaustion signature selecting for Spearman’s correlation coefficient as the default setting. GEPIA2 utilizes the non-log scale for calculation of significance. Finally, the ESTIMATE score [49] was utilized to show the immune score and tumour purity.

Correlations between ACOX2 expression and anti-cancer drug sensitivity in NSCLC cell lines

The relative expressions of ACOX2 in NSCLC cell lines and associated pharmacological profiles for 24 anti-cancer drugs were downloaded from the Cancer Cell Line Encyclopaedia (CCLE) database [50]. The half maximal inhibitory concentrations (IC50) values of 24 anticancer drugs in n = 89 NSCLC cell lines and a non-parametric correlation (Spearman) between ACOX2 gene expression and IC50 was conducted. For any drugs for which correlations were observed, the cell lines were subsequently divided into “sensitive” and “resistant” groups using a drug IC50 of 8 µM as the threshold according to Xiang et al. [51]. The expression differences for ACOX2 in these two groups were evaluated using a non-parametric Mann–Whitney test.

A further analysis for drugs potentially capable of targeting NSCLC based on ACOX2 expression was conducted using the DepMap PRISM repurposing Primary Screen [52].

Statistical analysis

All data are expressed as mean ± SEM unless stated otherwise. Statistical analysis was performed with Prism 5.01 (GraphPad, Ca, USA) using either paired two-tailed Student’s t-test or the Mann–Whitney two-tailed t-test. Correlations were conducted using the Spearman non-parametric correlation coefficient. Differences were considered significant when p < 0.05. All other analyses through web based servers used the default settings in each instance.

Results

Identification of ACOX2 as differentially expressed genes in NSCLC.

We assessed the expression ACOX2 mRNA in a panel of surgically resected fresh-frozen normal/tumour matched patient samples from Stage I and II patients by qPCR (Fig. 1A). Overall, levels of ACOX2 were shown to be significantly decreased across all samples (p < 0.0001). When stratified according to histology, downregulation of ACOX2 was found to be significant in both Adenocarcinomas (LUAD) (p = 0.0151) and Squamous Cell Carcinomas (LUSC) (p = 0.0008) respectively (Fig. 1A–C). We subsequently interrogated the TCGA LUAD and LUSC datasets for altered expression of ACOX2. In agreement with our data, significantly decreased expression of ACOX2 mRNA was observed for both LUAD (Fig. 1D; p = 1.7 × 10–4) and LUSC (Fig. 1E; p = 1.2 × 10–68) respectively. The other acyl-CoA oxidases were found to have significantly decreased expression of ACOX1 mRNA and ACOX3 mRNA in the LUSC samples only (Additional file 1: Table S1).

Fig. 1
figure 1

Altered expression of ACOX2 in NSCLC. Identification of decreased expression of ACOX2 in NSCLC. Examination of changes to ACOX2 mRNA levels in fresh-frozen surgically resected patient samples comprising A All Histologies; B Adenocarcinomas alone and C Squamous Cell Carcinomas. D Confirmatory comparative analysis of ACOX2 mRNA levels in the Cancer Genome Atlas (TCGA) Lung Adenocarcinoma (LUAD) and E the Lung Squamous Cell Carcinoma (LUSC) datasets using Lung Cancer Explorer (LCE) (19). F Altered expression of ACOX2 total protein levels in LUAD as assessed using cProSite. ***p < 0.001. G Altered expression of ACOX2 total protein levels in LUSC as assessed using cProSite. ***p < 0.001

Full size image

A further meta-analysis of a larger number of gene expression datasets was conducted, using a systematic analysis (Additional file 1: Table S2—Tumour vs Normal). From this analysis decreased expression of ACOX2 mRNA was found to be only significantly altered at the mRNA level in LUSC (Additional file 1: Table S2—p = 0.0031). A meta-analysis of the same datasets for altered ACOX2 mRNA expression between Tumour/Normal sample was conducted, and in this analysis ACOX2 mRNA is not associated with any mRNA differences overall—Additional file 6: Figure S5A. However, when separated according to histology LUAD mRNA remains non-significant (Additional file 7: Figure S6, p = 0.6), whilst LUSC remained significantly differentially expressed (Additional file 7: Figure S6, p = 0.0014).

We subsequently interrogated the CPTAC dataset for LUAD (23) and the results show that ACOX2 protein levels are significantly reduced in tumour samples compared to normal (Additional file 1: Table S1). At the time of writing the CPTAC data for LUSC is unavailable for analysis in UALCAN. However, access to the proteomic data for both LUAD and LUSC can now be found at cProSite. Using this we assessed the expression of ACOX1-3 in both LUAD and LUSC and the data is presented in Additional file 1: Table S1. From this analysis, levels of ACOX1 were not significantly altered in LUAD and LUSC. ACOX2 however was significantly downregulated at the protein level in both LUAD (p < 0.0001, Fig. 1F) and LUSC (p < 0.0001, Fig. 1G). Finally, ACOX3 protein was found to be significantly downregulated in LUAD for paired samples, whereas non-paired samples had no significant change to their overall protein levels, while ACOX3 protein was only significantly downregulated in LUSC (Additional file 1: Table S1).

Altered expression of acyl-CoA oxidases and other peroxisomal genes in NSCLC

Using STRING we identified a series of protein–protein interactions for acyl-CoA oxidases (ACOX1-3). The resulting network of first neighbours is presented in Additional file 2: Figure S1. Using systematic analysis in LCE we examined these genes for significant alterations in expression between tumour and normal lung, and the results are presented in Additional file 1: Table S2. This meta-analysis is presented as the standardized mean difference of Tumour-Normal gene expression using Hedges’ G as an effect size metric [34]. From this analysis significant alterations were observed for ACOX1 (p.adj—0.0011; LUAD); ACOX2 (p.adj—0.0031; LUSC); ACOX3 (p.adj—0.0011; LUSC). Other significant associations were observed for ACOXL, ACAA1, ACADM, ACOT4, ACSL1, ACSL5, AMACR, CAT, CRAT, DECR2, EHHADH, GNPAT, HSD17B4, PEX5, PEX14 and SCP2 (Additional file 1: Table S2). Of these, CAT, HSD17B4 and ACAA1 have recently been identified as showing either altered expression, or prognostic significance, or both in NSCLC [25].

Potential prognostic value of acyl-CoA oxidases and other peroxisomal genes in NSCLC

An initial analysis of the prognostic value of ACOX2 was conducted using KM-Plot [39]. Overall high expression of ACOX2 mRNA was associated with better overall survival (OS) in NSCLC (Fig. 2A). When the analysis was stratified according to histological subtype, this OS benefit was restricted solely to LUAD (Fig. 2B) and was not seen in LUSC (Fig. 2C). We subsequently examined the expression of ACOX2 in patient samples by IHC. In this regard no OS benefit was observed for ACOX2 overall or in any histological subtype (Fig. 2D–F; Additional file 3: Figure S2), suggesting that the observed OS benefit is restricted to mRNA expression (Additional file 4: Figure S3).

Fig. 2
figure 2

Prognostic value of ACOX2 in NSCLC. The prognostic value of ACOX2 expression was assessed for overall survival (OS) using KM-Plotter [39] and by IHC on a patient TMA. Higher expression of the mRNA for ACOX2 was associated with better OS overall (A); which when stratified by tumour histology was limited to the LUAD subtype (B); whilst no difference in OS was observed for LUSC (C). In contrast no significant OS benefit was observed by IHC for high ACOX2 protein expression overall (D); or in LUAD (E); or LUSC (F). p < 0.05 was considered to be significant

Full size image

Using LCE, we conducted a systematic analysis of ACOX2 mRNA for OS benefit. Overall, ACOX2 mRNA was found to have OS benefit in LUAD (p = 0.00049) but not in LUSC (Additional file 1: Table S3), while ACOX1 or ACOX2 mRNA had no survival benefit. However, when progression free survival (PFS) was assessed, ACOX2 mRNA was not found to show any significant PFS benefit (Additional file 5: Figure S4), whereas both ACOX1 and ACOX3 mRNA had significant associations with PFS, indicating that ACOX2 mRNA is solely linked with OS benefit. A meta-analysis of the same datasets for ACOX2 mRNA shown for all histologies for expression differences (Additional file 6: Figure S5A), and in this analysis ACOX2 mRNA is associated with OS (p < 0.01—Additional file 6: Figure S5B).

To evaluate the potential prognostic value of the acyl-CoA oxidases and their identified first neighbours, a systematic analysis was initially conducted using LCE and the results are provided in Additional file 1: Table S3. In agreement with Zhang et al., we confirm that altered expression of ACAA1 and CAT were associated with prognostic value in LUAD [25]. However, systematic analysis in LCE did not demonstrate any prognostic value for HSD17B4 in contrast to that previously observed [25]. Along with these genes, other genes which appear to show prognostic value in our analysis included ACOX2, ACOXL, ACOT4, ACSL5, AGPS, HADH and SCP2 (Additional file 1: Table S3). We subsequently compared OS values across five different database analyses (Additional file 1: Table S4). Those genes showing significant prognostic OS value were predominantly linked to the LUAD subtype (Additional file 1: Table S4), but some analyses suggest that altered expression of certain peroxisomal genes are associated with OS prognostic value in LUSC (for example ACOXL, ACSL5 and others in Additional file 1: Table S4). ACOXL and AGPS were significant for LUAD across all datasets/analyses, whilst CAT was significant across four of five (Additional file 1: Table S4). From this analysis we chose ACOX1, ACOX2, ACOX3, ACOXL, AGPS and CAT for further evaluation in LUAD using KM-Plot, and the results are shown in Additional file 4: Figure S3. Surprisingly, ACOX2 mRNA was not found to have any association with PFS (Additional file 5: Figure S4), whereas both ACOX1 and ACOX3 mRNA had significant associations with PFS alongside ACOXL, AGPS and CAT (Additional file 1: Table S5). From these analyses, high expression (stratified based on the median) of all 6 mRNAs is associated with significantly better OS in LUAD, and warrant further assessment.

Influence of Copy Number Variations on the expression of acyl-CoA oxidases in NSCLC

As acyl-CoA oxidases are significantly dysregulated in NSCLC (Additional file 1: Table S1), to further study the potential effects of this dysregulation we used cBioPortal [41] to assess for any correlations between copy number variations (CNVs) and gene expression changes in the TCGA –LUAD and –LUSC datasets. As shown in Fig. 3, positive correlations between CNV and gene expression were observed for all genes in both LUAD (Fig. 3A, C, and E) and LUSC (Fig. 3B, D, and F).

Fig. 3
figure 3

Correlations between acyl-CoA oxidase expression and copy number value in NSCLC. Correlations between acyl-CoA oxidase expression and copy number value (CNV) were examined using cBioPortal [41, 42]. Significant positive correlations between expression and CNV were observed in LUAD and LUSC for ACOX1 (A, B), ACOX2 (C, D) and ACOX3 (E, F). Spearman and Pearson correlation analyses were used and p < 10–5 was considered significant

Full size image

Correlations between ACOX2 expression and mutations in NSCLC

We subsequently assessed if the significant dysregulation of ACOX2 could be linked to the mutational status of the tumours, using TIMER2.0 we examined a number of genes commonly mutated in lung cancer (e.g. TP53) to determine whether mutations within these key genes were correlated with altered ACOX2 expression levels and the results are presented in Table 3. Of the genes examined, mutations in TP53, KRAS and CDKN2A were significantly correlated with altered ACOX2 expression for both LUAD and LUSC. Mutated ROS1 was significantly associated with altered ACOX2 expression only in LUAD, while mutated ERBB2 and ALK were associated with altered ACOX2 in LUSC.

Table 3 Correlations between ACOX2 mRNA expression and mutation of key genes in NSCLC. Gene Expression Correlations
Full size table

Finally, using muTarget [43], we then analysed whether mutations in any other genes may affect ACOX2 expression in LUAD and LUSC. From this analysis 244 genes were identified in LUAD which if mutated resulted in a significant alteration in ACOX2 expression, while 35 genes were identified in LUSC (Additional file 1: Tables S6 and S7) and the results of the top five mutated genes as defined by muTarget that affect ACOX2 expression were validated using TIMER2.0 and are presented in Fig. 4 for both LUAD (TP53; FAT2; PTPRZ1; GUCY1A3 and LRRC7—Fig. 4A) and LUSC (RASA1; PTCHD2; AMER3; GLRA2 and PCSK5—Fig. 4B).

Fig. 4
figure 4

Identification of mutated genes which affect ACOX2 expression in NSCLC. ACOX2 gene expression changes and mutation status in NSCLC were examined using muTarget [43]. The resulting analysis identified several genes which if were mutated resulted in significantly altered ACOX2 expression as follows: A TP53, FAT2, PTPRZ, GUCY1A3 and LRRC7 in LUAD and B RASA1, PTCHD2, AMER3, GLRA2 and PCSK5 in LUSC. All results were reassessed and plotted as shown using TIMER2 [44]

Full size image

Methylation analysis of ACOX2

As ACOX2 is significantly downregulated in both LUAD and LUSC using UALCAN we next examined whether methylation may be an important element in the downregulation of this gene. While significant methylation is presented for both LUAD (Fig. 5A) and LUSC (Fig. 5B), the required Beta-value cut-off for hypermethylation is (0.7–0.5) [53, 54] is only present for LUSC suggesting that DNA hypermethylation is a significant element in LUSC, but not LUAD (Fig. 5). UALCAN methylation analysis involves the Illumina 450 k Infinium chip and the associated CG probes showing significance were (CG16587010, CG13705284 and CG02259384). We then re-assessed methylation at this gene using MEXPRESS [45]. Analysis found significant negative correlations between methylation and gene expression particularly within CpG residues of the ACOX2 promoter (Additional file 1: Table S8). To assess this in more detail we re-analysed methylation using TCGA-Wanderer [46]. Selecting the Illumina 450 K settings the analysis found that the majority of probes in the LUAD dataset did not show significance with the exceptions of cg13705284, cg12075202, and cg23652987 (Fig. 5B and Additional file 1: Table S8) whereas all probes remained significant in the LUSC dataset (Fig. 5C and Additional file 1: Table S8). These results suggest that the downregulation of ACOX2 particularly in LUSC may be due in part to increased DNA CpG methylation within the ACOX2 promoter.

Fig. 5
figure 5

Methylation analysis of the ACOX2 promoter in NSCLC. The potential role of DNA CpG methylation was examined in LUAD and LUSC. A Analysis of methylation changes using UALCAN [37] and using stratification according to the default Beta-value cut-off for hypermethylation (0.7–0.5) the results suggest that hypermethylation is a significant element in LUSC, but not LUAD. Re-evaluation of methylation at the ACOX2 promoter in LUAD (B) and LUSC (C) using TCGA-Wanderer [46], showing clear differences in methylation at the promoter in LUSC. (*p < 0.05 is considered significant)

Full size image

Effects of decreased ACOX2 expression on immune infiltration

To assess the potential impact of decreased ACOX2 expression on tumour immunity, an analysis of tumour-infiltrating immune cells in NSCLC was conducted using TIMER [47]. Following purity-adjustment, spearman's rho and significance for the six immune cell types were generated and the results are presented in Table 4a. The results suggest that decreased expression of ACOX2 was negatively correlated with CD8 + T cell and Neutrophil infiltration in LUAD (Table 4a), whereas in contrast, decreased ACOX2 expression in LUSC was significantly positively associated with immune cell infiltration for all six immune cell types examined (Table 4a). Subsequently, when USO1 expression and immune cell infiltrates were examined for correlations with survival, only B-cell and dendritic cell immune infiltrates had survival benefit and this was further limited to LUAD only (Table 4b).

Table 4 Correlations between ACOX2 and immune infiltrations in NSCLC
Full size table

We re-assessed the effects of ACOX2 on immune infiltrates using TIMER2 [44], which provides a provides a more robust estimation of immune infiltration levels for The Cancer Genome Atlas (TCGA) datasets by using six state-of-the-art algorithms (TIMER, xCell, MCP-counter, CIBERSORT, EPIC and quanTIseq), and the results are presented in Additional file 1: Tables S9 and S10.

Finally, the ESTIMATE package [49] was used to generate an Estimate score (inferring tumour purity linked to gene expression). The Estimate scores (StromalScore, ImmuneScore, ESTIMATEScore, TumorPurity) are provided in Additional file 1: Table S11), and the graphed using TIMER2 [44, 55] as shown in Fig. 6.

Fig. 6
figure 6

Associations between ACOX2 expression and tumour purity. Associations between Tumour Purity and ACOX2 expression were generated using Timer2 [44] for A LUAD and B LUSC

Full size image

Correlation analysis between ACOX2 expression and immune cell exhaustion

To investigate the relationship between ACOX2 we used TIMER2.0 to assess correlations between its expression and the expression of important markers of T cell exhaustion [44, 56]. The markers chosen were PD-1 (PDCD1), CTLA4, LAG3, TIM-3 (HAVCR2), and GZMB and the results of this analysis are presented in Additional file 1: Table S12. After correlation adjustment by purity, ACOX2 expression was negatively correlated with the expression levels of PD-1, LAG3 and GZMB in LUAD (Additional file 1: Table S12), while its expression was positively correlated with PD-1, CTLA4 and TIM-3 in LUSC (Additional file 1: Table S12). A second assessment of T cell exhaustion was carried out using GEPIA2 [48] and a pre-defined set of T cell exhaustion markers (PDCD1, HAVCR2, TIGIT, LAG3, CXCL13 and LAYN). A similar pattern to that observed in TIMER was observed where a negative correlation between this 6 gene signature and ACOX2 expression with T cell exhaustion occurs in LUAD (r =  − 0.15, p = 0.001) with a positive correlation observed in LUSC (r = 0.23, p = 4.4 × 10−07).

Correlation analysis between ACOX2 expression and tumour mutational burden

Tumour mutational burden (TMB) is widely considered to be a biomarker for predicting potential patient response to immune checkpoint inhibitor therapy [57, 58]. Using the methodology described by Feng and Shen [59], we analysed the correlation between ACOX2 expression and biomarkers of tumour burden mutation, genes associated with either the DNA Damage Response (DDR) pathway or the Mismatch Excision Repair (MMR) pathway as a proxy for TMB (Table 5). The combined signatures for DDR pathway found a small significant negative correlation in both LUAD (p = 0.02), a very significant negative correlation was observed for LUSC (3.2 × 10–21). The signature for genes associated with MMR showed no correlations in LUAD, but again had a very significant negative correlation in LUSC (3.2 × 10–19).

Table 5 Correlation between ACOX2 expression and markers of Tumour Mutational Burden
Full size table

Whilst this represents a proxy for assessing the effect of ACOX2 expression on TMB, we subsequently reassessed the TCGA datasets for whether ACOX2 expression is associated with TMB using functions in UCSCXenaShiny [60] and also through by querying cBioPortal [41, 42]. As shown in Fig. 7A TMB is significantly associated with ACOX2 expression in LUAD, as reflected in the cBioPortal data (Fig. 7B). In contrast, TMB has no association with ACOX2 expression for LUSC (Fig. 7C). It must be noted however, that the cBioPortal does not have TMB data on all the patients, and the correlations observed reflect only a subset of patients. The fact that TMB may not have any influence on ACOX2 expression in LUSC may reflect either the loss of expression of ACOX2 expression by DNA CpG methylation in LUSC (Fig. 5), or a limitation in the numbers of samples available in the TCGA-LUSC dataset.

Fig. 7
figure 7

Associations between TMB and ACOX2 expression in NSCLC. The associations between TMB and ACOX2 expression were assessed in the TCGA datasets as follows: A XenaShiny analysis of the correlation between ACOX2 expression and TMB across all TCGA datasets; B cBioPortal analysis for ACOX2 expression and TMB in the TCGA-LUAD dataset; and C cBioPortal analysis for ACOX2 expression and TMB in the TCGA-LUSC dataset

Full size image

Taken together, these results suggest that TMB may be associated with ACOX2 expression particularly in LUAD.

ACOX2 and anti-cancer drug sensitivity in NSCLC

To examine any associations between ACOX2 expression and drug sensitivity, we examined the effects of ACOX2 expression against the pharmacological profiles for 24 anticancer drugs from the CCLE database [50]. Using parameters previously defined for IC50 for drug sensitive versus drug-resistant NSCLC [51], we identified three drugs for which lower ACOX2 mRNA expression was associated with resistance to 17-AAG (r = − 0.2399, p = 0.0344) a HSP90 inhibitor, PD-0325901 (r = -0.3472, p = 0.0018) a MEK inhibitor, and Crizotinib (r = − 0.2307, p = 0.0422) a c-MET/ALK inhibitor as shown in Fig. 6. A separate analysis using the DepMap PRISM repurposing Primary Screen [52] confirmed the potential to target MEK as it identified Tramatenib (GSK1120212) a MEK1/2 inhibitor FDA approved for the treatment of BRAF-mutated advanced NSCLC [61], and currently in a Phase I trial for KRAS mutated NSCLC [62].

Discussion

Metabolic pathway reprogramming is now considered to be a hallmark of cancer [63], and is an evolving area of therapeutic interest. Peroxisomes are key metabolic organelles [9], with key emerging roles in various cancers [9, 12]. In this manuscript we identified the family of acyl-CoA oxidases as being significantly dysregulated in NSCLC. In particular, we demonstrate that ACOX2 is significantly downregulated in NSCLC in both LUAD and LUSC (Fig. 1A–C). Re-analysis of existing TCGA datasets confirms this significantly decreased expression in larger well established datasets for NSCLC (Fig. 1D), and further show that the level of ACOX2 protein is significantly altered in the CTPAC dataset for NSCLC LUAD (Fig. 1E). However, whilst high mRNA expression of ACOX2 was found to be a potential prognostic biomarker for overall survival, in contrast expression of ACOX2 protein was not found to have any overall prognostic value (Fig. 2). Surprisingly, ACOX2 mRNA was not found to have any association with PFS (Additional file 5: Figure S4), whereas all other associated genes examined demonstrated associations with PFS benefit (Additional file 1: Table S5). PFS is often used as a surrogate measure of clinical benefit particularly for clinical trial drug approvals [64], but crossover and post-progression treatments may bias the relationship between surrogate endpoints such as PFS and OS [65], and discrepancies between PFS and OS often arise [66]. Nevertheless, when clinical trials are biomarker driven, strong evidence links PFS to response rate and survival benefit [67, 68].

Loss of ACOX2 expression due to mutation has been shown to have a functional role in cardiac cancer [69], with associated metabolic defects [70]. Building on from this we identified a series of additional peroxisomal genes which were also altered in NSCLC (Additional file 3: Figure S2, Additional file 1: Table S2). During the preparation of this manuscript, a publication emerged that identified that other key peroxisomal genes such as ACAA1were dysregulated in NSCLC [25, 59], and our analyses confirm/validate the observations from these papers (Additional file 1: Tables S2 and S3). Additional analyses indicate that the altered expression of ACOX2 (and other acyl-CoA oxidases) has significant prognostic value at least at the mRNA level predominantly in LUAD (Additional file 4: Figure S3). Other peroxisomal genes which showed prognostic value were also predominantly associated with the LUAD histological subtype (Additional file 4: Figure S3; Additional file 1: Tables S3 and S4).

Epigenetics plays a role in the dysregulation of many genes in cancer [71]. One of the best established mechanisms by which this occurs is via DNA CpG methylation at the promoter regions of genes [72]. Analysis of existing methylation data in the TCGA datasets demonstrates that altered methylation does occur at the promoter of ACOX2, and in particular hypermethylation occurs predominantly in the ACOX2 promoter in LUSC samples (Fig. 6). Given that the expression of ACOX2 is reduced in our patient LUSC samples (Fig. 1C) and confirmed as dramatically reduced in the TCGA LUSC dataset compared to LUAD (Fig. 1D), this may indicate that aberrant DNA methylation plays more important roles in the downregulation of ACOX2 in the squamous cell subtype. This is reflected in the DNA methylation at particular promoter specific methylation probes which are predominantly significant in the LUSC samples but not in LUAD (Additional file 1: Table S5) for example cg16587010. Increased promoter methylation may be responsible for the loss of ACOX2 expression in LUSC, but it does not appear to be a factor in the loss of ACOX2 expression in LUAD. As such it would appear that multiple elements may play a role in the dysregulation of ACOX2 in NSCLC. In this regard, mutation of certain lung cancer “driver” genes such as TP53, KRAS and CDKN2A was also associated with altered expression of ACOX2 in both LUAD and LUSC (Table 2) suggesting that cohorts of patients with such mutations could have altered peroxisomal pathways that may play a role in tumorigenesis. Other additional genes were also identified which when mutated correlated with altered ACOX2 expression. Examples of these such as FAT2 and LRRC7 have known roles in various cancers including lung cancer [73,74,75,76,77,78].

One emerging area in cancer where the peroxisomal pathway plays important roles is in the regulation of cellular immune responses [12]. Indeed, similar observations for immune infiltrating scores for LUAD and LUSC have recently been identified which have prognostic value [79,80,81]. In this manuscript we linked dysregulated expression of ACOX2 in NSCLC to altered infiltration of immune cells into the tumour environment. In this regard, in LUAD a negative correlation between ACOX2 expression and CD8 + T cell infiltration (Table 4a) occurs, whereas in LUSC a positive correlation between ACOX2 expression and CD8 + T cell infiltration is observed (Table 4a). As such it is possible that this observation may have potential implications for patient management for checkpoint inhibitors given that CD8 + tumour infiltrating lymphocytes are associated with better outcome in NSCLC [82, 83]. Moreover, in colon adenocarcinomas, reduced peroxisomal number and associated enzymatic activities occurs, and most recently an at-risk subpopulation has been identified where a low peroxisome pathway score is associated with a worse clinical outcome and high immune cell infiltration in CRC patients [84], and suggesting that peroxisomal genes such as ACOX2 could potentially be used in a similar fashion to explore for further associations with respect to utility in a predictive score for checkpoint inhibitor stratification or outcome in NSCLC.

It has also been suggested that peroxisomes play a role in the Warburg effect [85], whereby cancer cells utilize glycolysis as their predominant energy source. Thus competition for nutrients between cancer cells and immune cells is proposed to extend the function of the Warburg effect to a cell-extrinsic advantage depleting extracellular glucose in the tumour microenvironment, and thus rendering tumour-infiltrating T cells dysfunctional [59, 86]. Given the altered expression of a significant number of peroxisomal associated genes identified in this and other studies [25, 59], and the demonstration that altered cholesterol metabolism can affect the anti-tumour response of CD8( +) T cells [87], the data presented here further links dysregulation of peroxisome function with effects on T cell infiltration and activity in NSCLC.

It has been suggested that it may be possible to target peroxisomes by inducing “pexophagy” or the autophagocytic degradation of peroxisomes [12]. However, no agents have currently been developed that specifically target the peroxisome [12]. It may however be possible to utilise anti-lipolytic agents, or drugs such as small molecule inhibitors of specific peroxisome proteins could alternatively be targeted [12].

However, since the vast majority of patients with NSCLC display decreased ACOX2 at both the mRNA and protein level (Fig. 1), it may not be possible or clinically relevant to specifically target ACOX2 for the treatment of NSCLC patients. However, it may be possible to identify subsets of patients that may respond to other existing targeted therapies based on ACOX2 expression. In this regard, from our analysis increased sensitivity to established agents such as Tanespimycin (HSP90—Fig. 8A) and Crizotinib (c-MET/ALK—Fig. 8C) is associated with higher levels of ACOX2 expression. Moreover, MEK was identified as a candidate drug target associated with ACOX2 expression (Fig. 8B, D). MEK inhibitors are currently FDA approved for BRAF mutated advanced NSCLC [61], and are under investigation as a therapeutic option for KRAS mutated NSCLC [62]. As such it would be very interesting to compare levels of ACOX2 in a larger cohort of patients with KRAS mutations against a similar number of NSCLC wild-type for KRAS to determine if there may be any functional or biological link to sensitivity to MEK inhibitors. Indeed, a cohort of NSCLC patients with both RASA1 and NF-1 mutations has been shown to be uniquely sensitive to MEK inhibitors [78, 88], and as our analysis found that mutated RASA1 results in increased ACOX2 (Fig. 4) this lends further support to the possibility that cancers overexpressing ACOX2 may also be sensitive to MEK inhibitors. Despite the overall downregulation of ACOX2 protein in NSCLC (Fig. 1F, G), there is a subset of patients which demonstrate medium to high expression of this gene. Given the strong links we have identified between ACOX2 and mutated KRAS in both LUAD and LUSC (Table 2), it may be possible to expand the treatment options for MEK inhibitors to include those patients with high ACOX2 expression. As such it may therefore be possible to stratify patients into treatment or develop clinical trials using these agents based on expression levels of this gene, but additional investigations will be required to fully establish these possibilities.

Fig. 8
figure 8

Associations between ACOX2 expression in NSCLC cell lines and differential anti-cancer drug sensitivity. Analysis of the CCLE database [50] demonstrated that lower ACOX2 mRNA expression differences were found to be associated with resistance to A Tanespimycin (HSP90 inhibitor), B PD-0235901 (MEK inhibitor) and C Crizotinib (c-MET/ALK inhibitor). A non-parametric Mann–Whitney test was used to assess for significance with p < 0.05 considered to be the threshold for significant. D Separate analysis of the DepMAp PRISM repurposing Primary Screen [52] also identified trametinib (GSK1120212) a separate MEK1/2 inhibitor as a potential candidate drug linked with ACOX2 expression

Full size image

Nevertheless, it may be also possible to potentially restore levels of ACOX2 in patients where decreased ACOX2 is expressed through (PPAR)-α/γ dual activators such as SN158, which was shown to significantly upregulate ACOX2 mRNA [89], although such drugs have yet to enter the clinic due to serious side effects. Nevertheless, both in vitro and in vivo studies suggest that PPAR based agonists such as pioglitazone may have potential for treating NSCLC [90,91,92], however some degree of caution is indicated as other studies of such agonists suggest that when activated in myeloid cells of the tumour microenvironment, this may aggravate and promote lung cancer progression [2].

Overall, the results presented in this study raise interesting possibilities with respect to the potential use of ACOX2 in the diagnosis/prognosis of patients with NSCLC. Moreover, the results suggest that ACOX2 levels could potentially be useful to stratify patients into subsets that may respond to checkpoint inhibitor therapy. In addition, expression of ACOX2 can also potentially be used to stratify patients that may respond to various drugs including Tanespimycin, Crizotinib and MEK inhibitors. However, further studies will be required to elucidate the exact roles of these genes in the pathogenesis of LUAD and LUSC respectively, and also to further delineate the potential for use of ACOX2 expression to stratify or link patient responses to MEK inhibitors and potential subsets of patients with actionable mutations such as BRAF, KRAS, RASA1 and NF-1.

Availability of data and materials

The data that support the findings presented in this study are available for interrogation at the following online resources: TIMER: https://cistrome.shinyapps.io/timer/. TIMER2.0: http://timer.cistrome.org/. GEPIA2.0: http://gepia2.cancer-pku.cn/#index. LCE: http://lce.biohpc.swmed.edu/lungcancer/. UALCAN: http://ualcan.path.uab.edu/index.html. Oncolnc: http://oncolnc.org. KM-PLOT: https://kmplot.com/analysis/index.php?p=background. cBioPortal: https://www.cbioportal.org/. TCGA WANDERER: http://maplab.imppc.org/wanderer/#. MEXPRESS: https://mexpress.be/index.html. CCLE: https://portals.broadinstitute.org/ccle. MuTarget: https://www.mutarget.com/. DepMap: https://depmap.org/portal/. cProSite: https://cprosite.ccr.cancer.gov/#/. ESTIMATE: https://bioinformatics.mdanderson.org/estimate/index.html.

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  2. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

  3. Ackermann CJ, Adderley H, Ortega-Franco A, Khan A, Reck M, Califano R. First-line immune checkpoint inhibition for advanced non-small-cell lung cancer: state of the art and future directions. Drugs. 2020;80(17):1783–97.

    Article  PubMed  Google Scholar 

  4. Mustachio LM, Roszik J. Current targeted therapies for the fight against non-small cell lung cancer. Pharmaceuticals (Basel). 2020. https://doi.org/10.3390/ph13110374.

    Article  Google Scholar 

  5. Markham A. Pralsetinib: first approval. Drugs. 2020;80(17):1865–70.

    Article  CAS  PubMed  Google Scholar 

  6. Russo A, Lopes AR, McCusker MG, Garrigues SG, Ricciardi GR, Arensmeyer KE, Scilla KA, Mehra R, Rolfo C. New targets in lung cancer (excluding EGFR, ALK, ROS1). Curr Oncol Rep. 2020;22(5):48.

    Article  PubMed  Google Scholar 

  7. Rhodin J (1954) Correlation of ultrastructural organization and function in normal and experimentally changed proximal convoluted tubule cells of the mouse kidney. PhD Thesis. Karolinska Institute. Aktiebolaget Godvil (Publisher), 1–76

  8. Islinger M, Voelkl A, Fahimi HD, Schrader M. The peroxisome: an update on mysteries 2.0. Histochem Cell Biol. 2018;150(5):443–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Plett A, Charton L, Linka N. Peroxisomal cofactor transport. Biomolecules. 2020. https://doi.org/10.3390/biom10081174.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Lakshmi SP, Reddy AT, Banno A, Reddy RC. PPAR agonists for the prevention and treatment of lung cancer. PPAR Res. 2017;2017:8252796.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. Dahabieh MS, Di Pietro E, Jangal M, Goncalves C, Witcher M, Braverman NE, Del Rincón SV. Peroxisomes and cancer: the role of a metabolic specialist in a disease of aberrant metabolism. Biochim Biophys Acta Rev Cancer. 2018;1870(1):103–21.

    Article  CAS  PubMed  Google Scholar 

  12. Svensson RU, Shaw RJ. Lipid synthesis is a metabolic liability of non-small cell lung cancer. Cold Spring Harb Symp Quant Biol. 2016;81:93–103.

    Article  PubMed  Google Scholar 

  13. Dowling CM, Zhang H, Chonghaile TN, Wong KK. Shining a light on metabolic vulnerabilities in non-small cell lung cancer. Biochim Biophys Acta Rev Cancer. 2021;1875(1): 188462.

    Article  CAS  PubMed  Google Scholar 

  14. Majem B, Nadal E, Muñoz-Pinedo C. Exploiting metabolic vulnerabilities of Non small cell lung carcinoma. Semin Cell Dev Biol. 2020;98:54–62.

    Article  CAS  PubMed  Google Scholar 

  15. Ng Kee Kwong F, Nicholson AG, Harrison CL, Hansbro PM, Adcock IM, Chung KF. Is mitochondrial dysfunction a driving mechanism linking COPD to nonsmall cell lung carcinoma? Eur Respir Rev. 2017. https://doi.org/10.1183/16000617.0040-2017.

    Article  PubMed  Google Scholar 

  16. Honsho M, Okumoto K, Tamura S, Fujiki Y. Peroxisome biogenesis disorders. In: Lizard G, editor. Peroxisome biology: experimental models, peroxisomal disorders and neurological diseases. Cham: Springer International Publishing; 2020. p. 45–54.

    Chapter  Google Scholar 

  17. Kim JJ, Miura R. Acyl-CoA dehydrogenases and acyl-CoA oxidases Structural basis for mechanistic similarities and differences. Eur J Biochem. 2004;271(3):483–93.

    Article  CAS  PubMed  Google Scholar 

  18. Ferdinandusse S, Denis S, van Roermund CWT, Preece MA, Koster J, Ebberink MS, Waterham HR, Wanders RJA. A novel case of ACOX2 deficiency leads to recognition of a third human peroxisomal acyl-CoA oxidase. Biochim Biophys Acta Mol Basis Dis. 2018;1864(3):952–8.

    Article  CAS  PubMed  Google Scholar 

  19. Braverman NE, Raymond GV, Rizzo WB, Moser AB, Wilkinson ME, Stone EM, Steinberg SJ, Wangler MF, Rush ET, Hacia JG, et al. Peroxisome biogenesis disorders in the Zellweger spectrum: An overview of current diagnosis, clinical manifestations, and treatment guidelines. Mol Genet Metab. 2016;117(3):313–21.

    Article  CAS  PubMed  Google Scholar 

  20. Schutgens RB, Heymans HS, Wanders RJ, van den Bosch H, Tager JM. Peroxisomal disorders: a newly recognised group of genetic diseases. Eur J Pediatr. 1986;144(5):430–40.

    Article  CAS  PubMed  Google Scholar 

  21. Bjørklund SS, Kristensen VN, Seiler M, Kumar S, Alnæs GI, Ming Y, Kerrigan J, Naume B, Sachidanandam R, Bhanot G, et al. Expression of an estrogen-regulated variant transcript of the peroxisomal branched chain fatty acid oxidase ACOX2 in breast carcinomas. BMC Cancer. 2015;15:524.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Ko HK, Berk M, Chung YM, Willard B, Bareja R, Rubin M, Sboner A, Sharifi N. Loss of an androgen-inactivating and isoform-specific HSD17B4 splice form enables emergence of castration-resistant prostate cancer. Cell Rep. 2018;22(3):809–19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Zha S, Ferdinandusse S, Hicks JL, Denis S, Dunn TA, Wanders RJ, Luo J, De Marzo AM, Isaacs WB. Peroxisomal branched chain fatty acid beta-oxidation pathway is upregulated in prostate cancer. Prostate. 2005;63(4):316–23.

    Article  CAS  PubMed  Google Scholar 

  24. Zhang X, Yang H, Zhang J, Gao F, Dai L. HSD17B4, ACAA1, and PXMP4 in peroxisome pathway are down-regulated and have clinical significance in non-small cell lung cancer. Front Genet. 2020;11:273.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Chansky K, Detterbeck FC, Nicholson AG, Rusch VW, Vallières E, Groome P, Kennedy C, Krasnik M, Peake M, Shemanski L, et al. The IASLC lung cancer staging project: external validation of the revision of the TNM stage groupings in the eighth edition of the TNM classification of lung cancer. J Thorac Oncol. 2017;12(7):1109–21.

    Article  PubMed  Google Scholar 

  26. Detterbeck FC, Chansky K, Groome P, Bolejack V, Crowley J, Shemanski L, Kennedy C, Krasnik M, Peake M, Rami-Porta R. The IASLC lung cancer staging project: methodology and validation used in the development of proposals for revision of the stage classification of NSCLC in the forthcoming (eighth) edition of the TNM classification of lung cancer. J Thorac Oncol. 2016;11(9):1433–46.

    Article  PubMed  Google Scholar 

  27. Travis WD, Brambilla E, Burke AP, Marx A, Nicholson AG. Introduction to the 2015 world health organization classification of tumors of the lung, pleura, thymus, and heart. J Thorac Oncol. 2015;10(9):1240–2.

    Article  PubMed  Google Scholar 

  28. Travis WD, Brambilla E, Nicholson AG, Yatabe Y, Austin JHM, Beasley MB, Chirieac LR, Dacic S, Duhig E, Flieder DB, et al. The 2015 world health organization classification of lung tumors: impact of genetic, clinical and radiologic advances since the 2004 classification. J Thorac Oncol. 2015;10(9):1243–60.

    Article  PubMed  Google Scholar 

  29. Baird AM, Gray SG, O’Byrne KJ. IL-20 is epigenetically regulated in NSCLC and down regulates the expression of VEGF. Eur J Cancer. 2011;47(12):1908–18.

    Article  CAS  PubMed  Google Scholar 

  30. Cregan S, McDonagh L, Gao Y, Barr MP, O’Byrne KJ, Finn SP, Cuffe S, Gray SG. KAT5 (Tip60) is a potential therapeutic target in malignant pleural mesothelioma. Int J Oncol. 2016;48(3):1290–6.

    Article  CAS  PubMed  Google Scholar 

  31. Singh AS, Heery R, Gray SG. In silico and in vitro analyses of LncRNAs as potential regulators in the transition from the epithelioid to sarcomatoid histotype of malignant pleural mesothelioma (MPM). Int J Mol Sci. 2018. https://doi.org/10.3390/ijms19051297.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  33. Cai L, Lin S, Girard L, Zhou Y, Yang L, Ci B, Zhou Q, Luo D, Yao B, Tang H, et al. LCE: an open web portal to explore gene expression and clinical associations in lung cancer. Oncogene. 2019;38(14):2551–64.

    Article  CAS  PubMed  Google Scholar 

  34. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.

    Article  CAS  PubMed  Google Scholar 

  35. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, Varambally S. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Gillette MA, Satpathy S, Cao S, Dhanasekaran SM, Vasaikar SV, Krug K, Petralia F, Li Y, Liang WW, Reva B, et al. Proteogenomic characterization reveals therapeutic vulnerabilities in lung adenocarcinoma. Cell. 2020;182(1):200-225 e235.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Győrffy B, Surowiak P, Budczies J, Lánczky A. Online survival analysis software to assess the prognostic value of biomarkers using transcriptomic data in non-small-cell lung cancer. PLoS ONE. 2013;8(12): e82241.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Anaya J. OncoLnc: linking TCGA survival data to mRNAs, miRNAs, and lncRNAs. PeerJ Computer Science. 2016;2: e67.

    Article  Google Scholar 

  40. Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2(5):401–4.

    Article  PubMed  Google Scholar 

  41. Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):l1.

    Article  CAS  Google Scholar 

  42. Nagy Á, Győrffy B. muTarget: a platform linking gene expression changes and mutation status in solid tumors. Int J Cancer. 2021;148(2):502–11.

    Article  CAS  PubMed  Google Scholar 

  43. Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, Li B, Liu XS. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48(W1):W509–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Koch A, Jeschke J, Van Criekinge W, van Engeland M, De Meyer T. MEXPRESS update 2019. Nucleic Acids Res. 2019;47(W1):W561–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Díez-Villanueva A, Mallona I, Peinado MA. Wanderer, an interactive viewer to explore DNA methylation and gene expression data in human cancer. Epigenetics Chromatin. 2015;8:22.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  46. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Tang Z, Kang B, Li C, Chen T, Zhang Z. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019;47(W1):W556–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, Treviño V, Shen H, Laird PW, Levine DA, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.

    Article  PubMed  CAS  Google Scholar 

  49. Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, Wilson CJ, Lehár J, Kryukov GV, Sonkin D, et al. The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012;483(7391):603–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Xiang Z, Song S, Zhu Z, Sun W, Gifts JE, Sun S, Li QS, Yu Y, Li KK. LncRNAs GIHCG and SPINT1-AS1 are crucial factors for pan-cancer cells sensitivity to lapatinib. Front Genet. 2019;10:25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Corsello SM, Nagari RT, Spangler RD, Rossen J, Kocak M, Bryan JG, Humeidi R, Peck D, Wu X, Tang AA, et al. Discovering the anti-cancer potential of non-oncology drugs by systematic viability profiling. Nat Cancer. 2020;1(2):235–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Men C, Chai H, Song X, Li Y, Du H, Ren Q. Identification of DNA methylation associated gene signatures in endometrial cancer via integrated analysis of DNA methylation and gene expression systematically. J Gynecol Oncol. 2017;28(6): e83.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Sturm G, Finotello F, Petitprez F, Zhang JD, Baumbach J, Fridman WH, List M, Aneichyk T. Comprehensive evaluation of transcriptome-based cell-type quantification methods for immuno-oncology. Bioinformatics. 2019;35(14):i436–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Wherry EJ, Kurachi M. Molecular and cellular insights into T cell exhaustion. Nat Rev Immunol. 2015;15(8):486–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Goodman AM, Kato S, Bazhenova L, Patel SP, Frampton GM, Miller V, Stephens PJ, Daniels GA, Kurzrock R. Tumor mutational burden as an independent predictor of response to immunotherapy in diverse cancers. Mol Cancer Ther. 2017;16(11):2598–608.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Meng G, Liu X, Ma T, Lv D, Sun G. Predictive value of tumor mutational burden for immunotherapy in non-small cell lung cancer: A systematic review and meta-analysis. PLoS ONE. 2022;17(2): e0263629.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Feng H, Shen W. ACAA1 is a predictive factor of survival and is correlated with T cell infiltration in non-small cell lung cancer. Front Oncol. 2020;10: 564796.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Wang S, Xiong Y, Zhao L, Gu K, Li Y, Zhao F, Li J, Wang M, Wang H, Tao Z, et al. UCSCXenaShiny: an R/CRAN package for interactive analysis of UCSC xena data. Bioinformatics. 2021;38(2):527–9.

    Article  PubMed Central  CAS  Google Scholar 

  59. Odogwu L, Mathieu L, Blumenthal G, Larkins E, Goldberg KB, Griffin N, Bijwaard K, Lee EY, Philip R, Jiang X, et al. FDA approval summary: Dabrafenib and Trametinib for the treatment of metastatic non-small cell lung cancers harboring BRAF V600E mutations. Oncologist. 2018;23(6):740–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Lin SH, Lin HY, Verma V, Xu-Welliver M, Thall PF, Yao L, Kim PY, Gombos DS, Kawedia JD, Komaki R, et al. Phase I trial of definitive concurrent chemoradiotherapy and trametinib for KRAS-mutated non-small cell lung cancer. Cancer Treat Res Commun. 2022;30: 100514.

    Article  PubMed  Google Scholar 

  61. Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.

    Article  CAS  PubMed  Google Scholar 

  62. Hess LM, Brnabic A, Mason O, Lee P, Barker S. Relationship between progression-free survival and overall survival in randomized clinical trials of targeted and biologic agents in oncology. J Cancer. 2019;10(16):3717–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Hashim M, Pfeiffer BM, Bartsch R, Postma M, Heeg B. Do surrogate endpoints better correlate with overall survival in studies that did not allow for crossover or reported balanced postprogression treatments? An application in advanced non-small cell lung cancer. Value Health. 2018;21(1):9–17.

    Article  PubMed  Google Scholar 

  64. Shinawi T, Hill VK, Krex D, Schackert G, Gentle D, Morris MR, Wei W, Cruickshank G, Maher ER, Latif F. DNA methylation profiles of long- and short-term glioblastoma survivors. Epigenetics. 2013;8(2):149–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Fan Y, Yin G. Concordance index: surrogacy of progression-free survival for overall survival. Contemp Clin Trials. 2021;104: 106353.

    Article  PubMed  Google Scholar 

  66. Schwaederle M, Zhao M, Lee JJ, Eggermont AM, Schilsky RL, Mendelsohn J, Lazar V, Kurzrock R. Impact of precision medicine in diverse cancers: a meta-analysis of phase II clinical trials. J Clin Oncol. 2015;33(32):3817–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Schwaederle M, Zhao M, Lee JJ, Lazar V, Leyland-Jones B, Schilsky RL, Mendelsohn J, Kurzrock R. Association of biomarker-based treatment strategies with response rates and progression-free survival in refractory malignant neoplasms: a meta-analysis. JAMA Oncol. 2016;2(11):1452–9.

    Article  PubMed  Google Scholar 

  68. Zhou X, Wang H. ACOX2 deficiency in primary malignant cardiac tumors. Proc Natl Acad Sci U S A. 2017;114(18):E3590-e3591.

    CAS  PubMed  PubMed Central  Google Scholar 

  69. Zhou X, Xu M, Zeng W, Chen Z, Lu G, Gong Y, Finnell RH, Xiao H, Qiao B, Wang H. Combined effects of FH (E404D) and ACOX2 (R409H) cause metabolic defects in primary cardiac malignant tumor. Cell Death Discov. 2018;4:18.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  70. Portela A, Esteller M. Epigenetic modifications and human disease. Nat Biotechnol. 2010;28(10):1057–68.

    Article  CAS  PubMed  Google Scholar 

  71. Loaeza-Loaeza J, Beltran AS, Hernández-Sotelo D. DNMTs and impact of CpG content, transcription factors, consensus motifs, lncRNAs, and histone marks on DNA methylation. Genes (Basel). 2020. https://doi.org/10.3390/genes11111336.

    Article  Google Scholar 

  72. Santoni MJ, Kashyap R, Camoin L, Borg JP. The Scribble family in cancer: twentieth anniversary. Oncogene. 2020;39(47):7019–33.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  73. Hu M, Tan J, Liu Z, Li L, Zhang H, Zhao D, Li B, Gao X, Che N, Zhang T. Comprehensive comparative molecular characterization of young and old lung cancer patients. Front Oncol. 2021;11: 806845.

    Article  PubMed  Google Scholar 

  74. Yuan F, Lu L, Zou Q. Analysis of gene expression profiles of lung cancer subtypes with machine learning algorithms. Biochim Biophys Acta Mol Basis Dis. 2020;1866(8): 165822.

    Article  CAS  PubMed  Google Scholar 

  75. Mohammadoo Khorasani M, Karami Tehrani F, Parizadeh SMR, Atri M. Differential expression of alternative transcripts of soluble guanylyl cyclase, GYCY1a3 and GUCY1b3 genes, in the malignant and benign breast tumors. Nitric Oxide. 2019;83:65–71.

    Article  CAS  PubMed  Google Scholar 

  76. Xia Z, Ouyang D, Li Q, Li M, Zou Q, Li L, Yi W, Zhou E. The expression, functions, interactions and prognostic values of PTPRZ1: a review and bioinformatic analysis. J Cancer. 2019;10(7):1663–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Kitajima S, Barbie DA. RASA1/NF1-mutant lung cancer: racing to the clinic? Clin Cancer Res. 2018;24(6):1243–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Yang X, Shi Y, Li M, Lu T, Xi J, Lin Z, Jiang W, Guo W, Zhan C, Wang Q. Identification and validation of an immune cell infiltrating score predicting survival in patients with lung adenocarcinoma. J Transl Med. 2019;17(1):217.

    Article  PubMed  PubMed Central  Google Scholar 

  79. Zhang L, Chen J, Cheng T, Yang H, Li H, Pan C. Identification of the key genes and characterizations of tumor immune microenvironment in lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC). J Cancer. 2020;11(17):4965–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Theelen W, Krijgsman O, Monkhorst K, Kuilman T, Peters D, Cornelissen S, Ligtenberg MA, Willems SM, Blaauwgeers JLG, van Noesel CJM, et al. Presence of a 34-gene signature is a favorable prognostic marker in squamous non-small cell lung carcinoma. J Transl Med. 2020;18(1):271.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Khagi Y, Kurzrock R, Patel SP. Next generation predictive biomarkers for immune checkpoint inhibition. Cancer Metastasis Rev. 2017;36(1):179–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Schalper KA, Brown J, Carvajal-Hausdorf D, McLaughlin J, Velcheti V, Syrigos KN, Herbst RS, Rimm DL. Objective measurement and clinical significance of TILs in non-small cell lung cancer. J Natl Cancer Inst. 2015. https://doi.org/10.1093/jnci/dju435.

    Article  PubMed  PubMed Central  Google Scholar 

  83. Yin J, Wang H, Hong Y, Ren A, Wang H, Liu L, Zhao Q. Identification of an at-risk subpopulation with high immune infiltration based on the peroxisome pathway and TIM3 in colorectal cancer. BMC Cancer. 2022;22(1):44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Stine ZE, Schug ZT, Salvino JM, Dang CV. Targeting cancer metabolism in the era of precision oncology. Nat Rev Drug Discov. 2022;21(2):141–62.

    Article  CAS  PubMed  Google Scholar 

  85. Zeng W, Yin X, Jiang Y, Jin L, Liang W. PPARα at the crossroad of metabolic-immune regulation in cancer. Febs j. 2021. https://doi.org/10.1111/febs.16181.

    Article  PubMed  Google Scholar 

  86. Yang W, Bai Y, Xiong Y, Zhang J, Chen S, Zheng X, Meng X, Li L, Wang J, Xu C, et al. Potentiating the antitumour response of CD8(+) T cells by modulating cholesterol metabolism. Nature. 2016;531(7596):651–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Hayashi T, Desmeules P, Smith RS, Drilon A, Somwar R, Ladanyi M. RASA1 and NF1 are preferentially co-mutated and define A distinct genetic subset of smoking-associated non-small cell lung carcinomas sensitive to MEK inhibition. Clin Cancer Res. 2018;24(6):1436–47.

    Article  CAS  PubMed  Google Scholar 

  88. Jung Y, Cao Y, Paudel S, Yoon G, Cheon SH, Bae GU, Jin LT, Kim YK, Kim SN. Antidiabetic effect of SN158 through PPARα/γ dual activation in ob/ob mice. Chem Biol Interact. 2017;268:24–30.

    Article  CAS  PubMed  Google Scholar 

  89. Kiran A, Kumari GK, Krishnamurthy PT. Preliminary evaluation of anticancer efficacy of pioglitazone combined with celecoxib for the treatment of non-small cell lung cancer. Invest New Drugs. 2022;40(1):1–9.

    Article  CAS  PubMed  Google Scholar 

  90. To KKW, Wu WKK, Loong HHF. PPARgamma agonists sensitize PTEN-deficient resistant lung cancer cells to EGFR tyrosine kinase inhibitors by inducing autophagy. Eur J Pharmacol. 2018;823:19–26.

    Article  CAS  PubMed  Google Scholar 

  91. Kim TW, Hong DW, Park JW, Hong SH. CB11, a novel purine-based PPARɣ ligand, overcomes radio-resistance by regulating ATM signalling and EMT in human non-small-cell lung cancer cells. Br J Cancer. 2020;123(12):1737–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  92. Sippel TR, Johnson AM, Li HY, Hanson D, Nguyen TT, Bullock BL, Poczobutt JM, Kwak JW, Kleczko EK, Weiser-Evans MC, et al. Activation of PPARγ in myeloid cells promotes progression of epithelial lung tumors through TGFβ1. Mol Cancer Res. 2019;17(8):1748–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Data used in this publication were generated by the Clinical Proteomic Tumor Analysis Consortium (NCI/NIH).

Funding

The research in this article was supported by the St James’s Hospital Foundation (N/C: 2657).

Author information

Authors and Affiliations

  1. Thoracic Oncology Research Group, Laboratory Medicine and Molecular Pathology, Central Pathology Laboratory, St. James’s Hospital, Dublin, D08RX0X, Ireland

    Jane S. Y. Sui, Petra Martin, Anna Keogh, Stephen P. Finn & Steven G. Gray

  2. Department of Medicine, Thoracic Oncology Service, Memorial Sloan Kettering Cancer Center, New York, USA

    Jane S. Y. Sui

  3. Midland Regional Hospital Tullamore, Tullamore, Ireland

    Petra Martin

  4. Department of Histopathology and Morbid Anatomy, Trinity College Dublin, Dublin, Ireland

    Pierre Murchan & Stephen P. Finn

  5. Department of Histopathology, Labmed Directorate, St. James’s Hospital, Dublin, Ireland

    Lisa Ryan, Siobhan Nicholson & Stephen P. Finn

  6. HOPE Directorate, St James’s Hospital, Dublin, Ireland

    Sinead Cuffe

  7. School of Mathematics, Statistics, and Applied Mathematics, National University of Ireland Galway, Galway, Ireland

    Pierre Murchan & Pilib Ó Broin

  8. Cancer Molecular Diagnostics, Labmed Directorate, St. James’s Hospital, Dublin, Ireland

    Stephen P. Finn

  9. Surgery, Anaesthesia and Critical Care Directorate, St James’s Hospital, Dublin, Ireland

    Gerard J. Fitzmaurice, Ronan Ryan & Vincent Young

  10. Department of Clinical Medicine, Trinity College Dublin, Dublin, Ireland

    Steven G. Gray

  11. School of Biological Sciences, Technological University Dublin, Dublin, Ireland

    Steven G. Gray

Authors
  1. Jane S. Y. Sui
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Petra Martin
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Anna Keogh
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Pierre Murchan
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Lisa Ryan
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Siobhan Nicholson
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Sinead Cuffe
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Pilib Ó Broin
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Stephen P. Finn
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Gerard J. Fitzmaurice
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Ronan Ryan
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Vincent Young
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Steven G. Gray
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

JSYS, PM (1) and SGG designed the experiments and wrote the paper; LR, JSYS, PM (1), PM (2) and SGG performed the experiments. SN, SC, SPF, RR, GJF and VY contributed access to sample collection, clinical and histopathological support. PM (2) and PÓB assisted with some of the bioinformatic analyses. AK and SPF developed the immunoscore and scored the IHC TMA. SPF, JSYS and SGG read and approved the final manuscript for publication. All authors reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Steven G. Gray.

Ethics declarations

Ethics approval and consent to participate

The study was conducted in accordance with the Declaration of Helsinki, in-line with national and international guidelines and approved by the authors’ institutional board. Informed consent was obtained from each patient, and the study was conducted after formal approval from the SJH/AMNCH Hospital Ethics Committee—Ethics REC (No.: 041018/8804).

Consent for publication

Consent for publication is not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

. Analysis of TCGA-LUAD and –LUSC datasets for alterations in acyl-CoA peroxidases at the mRNA and protein level.

Additional file 2

: Figure S1. Protein-Protein Interaction Networks associated with acyl-CoA oxidases. Functional enrichment analysis of acyl-CoA oxidases was carried out on the STRING database [34] to identify Protein-Protein Interaction Networks. The results were imported into Cytoscape [34] and first neighbours highlighted. The resulting first-neighbours are shown.

Additional file 3

: Figure S2. IHC of ACOX2. Representative examples of ACOX2 protein expression in NSCLC. (A) Negative staining (B) Positive Staining . Images were taken at 40x.

Additional file 4

: Figure S3. Prognostic value of acyl-CoA oxidases and other peroxisomal genes in NSCLC. The prognostic effects of ACOX2, other acyl-CoA oxidases and key other peroxisomal associated genes were assessed for overall survival (OS) using KM-Plotter [38]. Analyses were conducted on the gene chip datasets for LUAD. Higher expression of the mRNA for (A) ACOX1, (B) ACOX2, (C) ACOX3, (D) ACOXL, (E) AGPS and (F) CAT was found to be associated with significantly better OS in LUAD, with p<0.05 considered to be significant.

Additional file 5

: Figure S4. Prognostic value of acyl-CoA oxidases mRNA expression in NSCLC as assessed by Progression Free Survival. Progression free survival (PFS) is defined First Progression (FP) and KM-Plot was used to analyse ACOX1-3 across the NSCLC. The results are presented as follows: (A) ACOX1 – all histologies; (B) ACOX1 – LUAD; (C) ACOX1 – LUSC; (D) ACOX2 – all histologies; (E) ACOX2 – LUAD; (F) ACOX2 – LUSC; (G) ACOX3 – all histologies; (H) ACOX3 – LUAD; (I) ACOX3 – LUSC.

Additional file 6

: Figure S5. Meta-analysis of ACOX2 expression in multiple NSCLC datasets. A meta-analysis was conducted on LCE [33] for a large number of gene expression datasets generating forest plots summarizing (A) tumour - normal standardized mean difference for tumour vs normal meta-analysis and (B) hazard ratios for OS meta-analysis for all NSCLC datasets.

Additional file 7

: Figure S6. Meta-analysis of ACOX2 expression in multiple NSCLC datasets. Meta-analysis conducted on LCE [32] generating forest plots summarizing hazard ratios for OS in (A) LUAD and (B) LUSC specific datasets.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sui, J.S.Y., Martin, P., Keogh, A. et al. Altered expression of ACOX2 in non-small cell lung cancer. BMC Pulm Med 22, 321 (2022). https://doi.org/10.1186/s12890-022-02115-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12890-022-02115-7

Keywords

BMC Pulmonary Medicine

ISSN: 1471-2466

Contact us