Lipid biomarkers that reflect postoperative recurrence risk in lung cancer patients who smoke: a case–control study
Lipids in Health and Disease volume 22, Article number: 15 (2023)
The risk of postoperative recurrence is higher in lung cancer patients who smoke than non-smokers. However, objective evaluation of the postoperative recurrence risk is difficult using conventional pathological prognostic factors because of their lack of reproducibility. Consequently, novel objective biomarkers that reflect postoperative risk in lung cancer patients who smoke must be identified. Because cigarette smoking and oncogenesis alter lipid metabolism in lung tissue, we hypothesized that the lipid profiles in lung cancer tissues are influenced by cigarette smoking and can reflect the postoperative recurrence risk in smoking lung cancer patients. This study aimed to identify lipid biomarkers that reflect the smoking status and the postoperative recurrence risk.
Primary tumor tissues of lung adenocarcinoma (ADC) (n = 26) and squamous cell carcinoma (SQCC) (n = 18) obtained from surgery were assigned to subgroups according to the patient’s smoking status. The ADC cohort was divided into never smoker and smoker groups, while the SQCC cohort was divided into moderate smoker and heavy smoker groups. Extracted lipids from the tumor tissues were subjected to liquid chromatography-tandem mass spectrometry analysis. Lipids that were influenced by smoking status and reflected postoperative recurrence and pathological prognostic factors were screened.
Two and 12 lipid peaks in the ADC and SQCC cohorts showed a significant positive correlation with the Brinkman index, respectively. Among them, in the ADC cohort, a higher lipid level consisted of three phosphatidylcholine (PC) isomers, PC (14:0_18:2), PC (16:1_16:1), and PC (16:0_16:2), was associated with a shorter recurrence free period (RFP) and a greater likelihoods of progressed T-factor (≥ pT2) and pleural invasion. In the SQCC cohort, a lower m/z 736.5276 level was associated with shorter RFP and greater likelihood of recurrence.
From our data, we propose three PC isomers, PC (14:0_18:2), PC (16:1_16:1), and PC (16:0_16:2), and a lipid peak of m/z 736.5276 as novel candidate biomarkers for postoperative recurrence risk in lung ADC and SQCC patients who are smokers.
Cigarette smoking is a strong risk factor for developing non-small cell lung cancer (NSCLC). Around 90% of NSCLC cases are associated with active cigarette smoking . Although the standard treatment modality for resectable stage I–III NSCLC is radical surgery , the postoperative recurrence risk for current smokers is reported to be high compared with non-smokers and former smokers (relative risks of recurrence: current smoker, 1; non-smoker, 0.45; former smoker; 0.54) . Hence, predicting the postoperative recurrence risk in individual NSCLC patients who are smokers will contribute to developing a qualified postoperative therapeutic strategy.
Currently, the postoperative recurrence risk of NSCLC patients is evaluated using the tumor, lymph node, and metastasis (TNM) classification system  and several histopathological prognostic factors, such as lymph node metastasis , pleural invasion , lymphatic vessel invasion , blood vessel invasion [8, 9], and spread through air space (STAS) . However, evaluating the risk of postoperative recurrence using these conventional prognostic factors remains challenging because of subjective judgment and a lack of reproducibility . Therefore, identifying novel and objective biomarkers that reflect postoperative recurrence risk in smokers NSCLC patients is needed.
With the recent advancements in mass spectrometry technology, lipidomics has been utilized as a new field for exploring cancer biomarkers . Cancer cells exhibit a common phenotype of uncontrolled proliferation, and their lipid metabolism is altered to generate energy and biomass components efficiently [13, 14]. Specific lipids in cancer tissues have been reported as potential biomarkers reflecting patient prognosis and the risk of postoperative recurrence. For example, an elevated level of phosphatidylcholine (PC) (32:1) is suggested to be a candidate predictor for postoperative recurrence of primary triple-negative breast cancer . In clear cell renal cell carcinoma, oleic acid attenuation is associated with a shorter progression-free period . In the field of NSCLC, lung cancer tissue has been demonstrated to have distinct lipid profiles for their cancer differentiation [17,18,19,20], histological subtypes [17,18,19,20], and oncogenic driver gene mutation , reflecting their lipid metabolism alterations. Furthermore, we have reported specific sphingomyelin (SM) species in lung adenocarcinoma (ADC) and squamous cell carcinoma (SQCC) as novel candidate predictors for postoperative recurrence after radical surgery [21, 22]. Thus, these previous works led us to focus on lipidomics to explore novel biomarkers that reflect postoperative recurrence risk in smoking NSCLC patients.
Cigarette smoking alters lipid metabolism in the lung tissue. Short-term cigarette exposure induces decreased PC levels in the alveolar type II cells of the mouse lung [23, 24]. Additionally, decreased levels of surfactant PC are observed in bronchoalveolar lavage fluids from human lungs exposed to cigarette smoke . PC is a significant surfactant component [23,24,25], and specific PC species can suppress lung cancer growth and metastasis by decreasing antiapoptotic factors and matrix metallopeptidases [26, 27]. From these previous studies, we hypothesized that NSCLC tissue lipid profiles are influenced by cigarette smoking and reflect the risk of postoperative recurrence in NSCLC smoking patients. However, the influence of cigarette smoking on the lipid profiles of lung cancer tissues has not been previously investigated.
Therefore, this study aims to identify NSCLC lipid biomarkers that can reflect the patient’s smoking status and postoperative recurrence risk using liquid chromatography-tandem mass spectrometry (LC–MS/MS). Such biomarkers may contribute to evaluating postoperative recurrence risk and developing a qualified postoperative therapeutic strategy for smoking NSCLC patients.
Patients and tissue samples
Among lung ADC and SQCC patients who received surgery with complete resection from January 2013 to December 2016 at Hamamatsu University Hospital, cases with available retrospective frozen tissue samples were enrolled. Primary tumor tissue samples were frozen in liquid nitrogen immediately after resection and stored at -80 °C. Postoperative patients were followed-up with computed tomography (CT) and blood examination of carcinoembryonic antigen (CEA), squamous cell carcinoma antigen (SCC), and cytokeratin 19 fragment (CYFRA) every three months during the first two years, then every six months until five years after the surgery. Tumor marker elevation (CEA ≥ 5.0 ng/mL, SCC ≥ 2.5 ng/mL, CYFRA ≥ 3.5 ng/mL) without recurrent findings on CT was screened for brain or bone metastasis with head magnetic resonance imaging (MRI) and systemic positron emission tomography (PET).
For patient selection, we retrospectively enrolled frozen tissue samples obtained from NSCLC patients with ADC or SQCC. Patients who received induction chemotherapy or radiotherapy, which could cause artifacts in lipid analysis, were excluded. We defined recurrence as radiological imaging-based findings of distant or locoregional recurrence within five years after surgery.
Three μm-thick sections of paraffin-embedded tissue blocks were used. Hematoxylin and eosin (H&E) staining was used to evaluate the adenocarcinoma subtype, tumor size, lymph node metastasis, and STAS. D2-40 and Elastica van Gienson staining were used to assess lymphatic and blood vessel invasion, respectively. For the enrolled cases, histopathological diagnoses and pathological staging were performed by experienced pathologists according to the World Health Organization criteria  and the 8th edition of the TNM classification for lung and pleural tumors , respectively. In principle, an adenocarcinoma phenotype was confirmed by the positivity of thyroid transcription factor-1 (TTF-1) or napsin A. On the other hand, a squamous phenotype was confirmed with the positivity of one of the squamous markers (p40, p63, and cytokeratin 5/6) and negativity for TTF-1.
Lipid extraction from frozen tissue samples
As described previously, lipids were extracted from the frozen tissue samples [21, 22]. Briefly, each tissue was weighed using a Sartorius analytical lab balance CPA224S (Sartorius AG, Göttingen, Germany) (Additional file 1, Supplemental Fig. 1). Then, the modified Bligh-Dyer method was applied to the tissue samples for lipid extraction. During the extraction, 1.6 mmol of 1,2-dilauroyl-sn-glycero-3-PC (Avanti Polar Lipids, Alabaster, AL, USA), PC (12:0_12:0) per 1 mg of sample was added to standardize the lipid levels. The extracted lipids were dried using miVac Duo LV (Genevac, Ipswich, UK). We dissolved the dried lipids with methanol proportionally to the tissue weights to keep the PC (12:0_12:0) levels similar among the cases.
Lipid analysis by LC–MS/MS
As reported previously, the diluted lipids were subjected to LC–MS/MS [21, 22]. Briefly, 2–10 μL of the diluted lipids were applied to an Acclaim 120 C18 column (150 mm × 2.1 mm, three μm) (Thermo Fisher Scientific, Waltham, MA, USA) and analyzed using a Q Exactive™ Hybrid Quadrupole-Orbitrap™ Mass Spectrometer (Thermo Fisher Scientific) with an UltiMate 3000 (Thermo Fisher Scientific). Spectral data were recorded using Xcalibur v3.0 Software (Thermo Fisher Scientific). Lipid peak identification and semi-quantification were performed using LipidSearch™ software version 4.2.13 (Mitsui Knowledge Industry, Tokyo, Japan) with the same parameter settings as reported previously [21, 22]. We aligned the identified peaks with a retention time (RT) tolerance of 1.0 min. Identified redundant peaks with different RTs were regarded as independent structural isomers (shown as “Duplication” in Additional File 2, sheet 1). The relative abundance of the identified lipid peaks was measured by calculating their area value. An area value of a lipid species was divided by that of the standard PC (12:0_12:0) in the corresponding case for standardization. The complete list of identified lipid peaks with standardized area values is presented in Additional File 2, sheet 1.
To screen lipids influenced by cigarette smoking, ADC and SQCC patient cohorts were divided into two subgroups according to the smoking status. The ADC cohort was divided into never smoker (Brinkman index [BI] = 0) and smoker (BI > 0) groups to detect the influence of cigarette smoking on lipid profiles as much as possible. On the other hand, the SQCC cohort was divided into moderate smoker (BI < 1590) and heavy smoker (BI > 1590) groups bordering on a median BI value of 1590, as the cohort did not include never smokers. P-values were calculated using Welch’s t-test between the subgroups for each area value of lipid species. Additionally, the fold change of respective lipid peaks was calculated as follows: mean area value of smoker divided by that of never smoker in the ADC cohort; mean area value of heavy smoker divided by that of moderate smoker in the SQCC cohort. Then, volcano plots with -log10 (P-value) for the vertical axis and log2 (fold change) for the horizontal axis were generated for the two cohorts, and significantly abundant lipids in the respective subgroups were screened. P-values of < 0.05 and fold changes ≥ 2.0 or ≤ 0.5 were considered significant. The heat map display and Pearson correlation analysis of the relative intensities of the screened significant lipid species and BI were performed using MetaboAnalyst 5.0 (https://www.metaboanalyst.ca/). The lipid peak intensity values were divided by each lipid peak's standard deviation for standardization. Thus, lipid peaks significantly correlated with BI (Pearson correlation coefficient ≥ 0.5 ∩ P-value < 0.05) were selected as lipids influenced by cigarette smoking.
We evaluated the postoperative recurrence risk by examining endpoints, including recurrence-free period (RFP) and odds ratios (ORs) on recurrence, pathological T-factor, lymph node metastasis, pleural invasion, blood vessel invasion, lymphatic vessel invasion and STAS. To calculate the odds ratio (OR), all pathological prognostic factors as objective variables were binarized as follows: T-factor; ≤ pT1 and ≥ pT2, differentiation; ≤ G1 and ≥ G2, and all other factors into negative and positive.
Associations between patient characteristics and smoking status were evaluated by the Fisher exact test (categorical variables) or Mann–Whitney U-test (continuous variables). The Welch’s t-test used in the volcano plots was performed with “TTEST” of Excel™ (Microsoft, Redmond, WA, USA). The cutoff values discriminating the recurrent and non-recurrent cases were calculated by receiver operating characteristic (ROC) curve analysis. The area under the ROC curve (AUC) values were calculated for the discrimination ability. The RFP was determined as the time from operation until the first disease recurrence. RFP curves were generated using the Kaplan–Meier (KM) method. Correlation analyses were performed by calculating Spearman’s rank correlation coefficient (rs). Except for Welch’s t-test, all statistical analyses were performed using R (The R Foundation for Statistical Computing, Vienna, Austria, version 3.6.2). P-values < 0.05 were considered statistically significant.
Patient clinicopathological characteristics
Patient characteristics are shown in Table 1. Overall, 26 patients in the ADC cohort were divided into never smoker (n = 8) and smoker (n = 18) groups, while 18 patients in the SQCC cohort were divided into moderate smoker (n = 9) and heavy smoker (n = 9) groups. In the ADC cohort, more female patients were enrolled in the never smoker group (n = 6 [75%]) than in the smoker group (n = 1 [5.9%]). Median BI values were 0 (never smoker) and 663 (smoker) in the ADC cohort and 920 (moderate smoker) and 2100 (heavy smoker) in the SQCC cohort. No significant difference was observed in any other characteristic examined. The median follow-up period for the entire cohort was 60 months (range of 4–90 months).
Screening of lipids influenced by cigarette smoking
In the LC–MS/MS analysis of the lipids extracted from the frozen tissue samples, LipidSearch™ software identified a total of 2,453 lipid peaks. The complete list of the identified lipid peaks is presented in Additional file 2, sheet 1. Lipid intensities between the never smoker and smoker groups in the ADC cohort and between the moderate smoker and heavy smoker groups in the SQCC cohort were compared using volcano plots (Fig. 1). In the ADC cohort, four and 29 significantly high lipid peaks were identified in the never smoker and smoker groups, respectively. In the SQCC cohort, 15 lipid peaks were significantly high in the heavy smoker group, but no significant lipid peaks were identified in the moderate smoker group. The relative intensities of these screened lipid peaks and BI values are depicted using heat maps (Fig. 2). Overall, each screened lipid peak reproducibly showed different intensities between the subgroups in the ADC (Fig. 2a) and SQCC (Fig. 2b). Intriguingly, the majority of high-intensity lipid peaks in the heavy smoker group of the SQCC cohort showed a trend of positive correlation with BI (Fig. 2b). In contrast, a weak positive and negative correlation trend in the smoker group of the ADC cohort was observed (Fig. 2a).
These findings led us to perform a correlation analysis between the screened lipid peaks and BI. Only two of the 29 lipid peaks in the ADC cohort (Fig. 3a) showed a significant correlation with BI. On the other hand, 12 of the 15 lipid peaks showed significance in the SQCC cohort (Fig. 3b). The correlation tables with P-values are presented in Supplemental file 2, sheet 2. From these results, we considered these lipid peaks that showed significant correlations with BI as lipids that were influenced by cigarette smoking.
Identification of lipids that reflect postoperative recurrence and histopathological prognostic factors
The fourteen lipid peaks showing a significant correlation with BI were subjected to AUC analysis to evaluate the discrimination ability between the recurrent and non-recurrent cases (Table 2). Two and four lipid peaks in the ADC and SQCC cohorts that showed AUC ≥ 0.6 were selected, respectively. Then, we performed an RFP curve analysis on these lipid peaks using the cutoff values calculated in the AUC analysis. In the case of m/z 730.5381 (ID: 944) in the ADC cohort, the RFP of the high intensity group was significantly shorter than in the low intensity group (P = 0.044, median RFP [MRFP]: 15.5 vs 60.5 months) (Fig. 4a). Notably, with an m/z of 736.5276 (ID: 1033) in the SQCC cohort, recurrence was not seen in the high intensity group. Accordingly, the P-value of this RFP analysis did not meet the significance level (P = 0.645, MRFP: 45.5 vs 56.0 months) regardless of the observed difference in recurrent prognosis between the low and high intensity groups. Hence, we regarded m/z 736.5276 as a valid predictor for postoperative recurrence.
The tandem mass spectrometry analysis of m/z 730.5381 and m/z 736.5276 in the ADC and SQCC cohorts were attempted using tandem mass spectral data. The mass spectrum peak of m/z 730.5381 was composed of three different PC isomers of PC (14:0_18:2), PC (16:1_16:1), and PC (16:0_16:2) (Additional file 1, Supplemental Fig. 2). We defined the combined area value of these three PC species as ‘PC cluster’. In contrast, the structural identification of m/z 736.5276 was not achieved because of its weak intensities in all samples. Using LipidSearch™ software, PC (O-14:1_20:5) was assigned as a candidate molecule of m/z 736.5276.
Then, we evaluated correlations between the identified lipids and recurrence or pathological prognostic factors with ORs. The higher PC cluster (m/z 730.5381) level of the ADC cohort was significantly associated with greater likelihoods of progressed T-factor (≥ pT2) (OR: 7.33, 95% confidence interval [CI]: 1.27–42.30, rs: 0.46) and pleural invasion (OR: 7.33, 95% CI: 1.27–42.30, rs: 0.46). Additionally, recurrence (OR: 3.50, 95% CI: 0.69–17.90, rs: 0.30), differentiation (≥ G2) (OR: 6.11, 95% CI: 0.60–62.20, rs: 0.32), and STAS (OR: 6.00, 95% CI: 0.92–39.20, rs: 0.39) showed tendencies of moderate likelihoods, though they did not reach statistical significance (Fig. 5a). The lower m/z 736.5276 level of the SQCC cohort was strongly associated only with a greater likelihood of recurrence (OR: 3.14E + 08, 95% CI: 0.00-infinity, rs: 0.54) (Fig. 5b).
With this result, we identified lipid biomarkers that reflected postoperative recurrence risk and pathological prognostic factors of lung ADC and SQCC patients who are smokers.
In this study, we identified lipid biomarkers influenced by cigarette smoking and reflected postoperative recurrence risk of lung ADC and SQCC. Overall, the findings are compatible with our initial hypothesis.
All the lipid peaks significantly influenced by cigarette smoking showed a positive correlation with BI; lipid peaks showing significant negative correlations were not detected. Previous studies have demonstrated negative correlations between cigarette smoking and lipid levels, mainly PC, in normal lung tissues. The mechanism of this negative correlation has largely been explained by decreased PC production in the alveolar type II cells [23,24,25]. The positive correlation we observed in this study between cigarette smoking and lipid levels in NSCLC tissues, the mechanism of which remains unknown, showed an inverse trend compared with that reported in normal lung tissues in previous studies. This mechanism may be explained by studies investigating the correlation between lipid metabolism and oxidative stress. Cigarette smoking generates reactive oxygen species (ROS) from free radicals, which is a significant oxidative stress factor for cells [29, 30]. Cancer cells exposed to oxidative stress produce a large number of lipid droplets that promotes stress resistance and homeostasis . Actually, a higher intracellular ROS level is demonstrated to be associated with a higher lipid accumulation . ROS damage mitochondria resulting in a further increase in ROS levels. Then, the increased ROS levels activate the sterol regulatory element binding protein pathway that leads to lipid accumulation .
A greater number of lipid peaks with significant positive correlations were identified in the SQCC cohort (12 peaks) than in the ADC cohort (two peaks). Cigarette smoking has the strongest association with the oncogenesis of SQCC and small cell carcinoma among all the lung cancer histological subtypes. The relative risk of non-ADCs associated with cigarette smoking is approximately ten while doubling the risk of ADC . During cancer development, lipid metabolism is reprogrammed to satisfy the strong lipid avidity of the transformed cells. Endogenous lipid synthesis and uptake of exogenous lipids are enhanced to maintain a highly proliferative phenotype [13, 14]. Therefore, a greater number of lipids with a significant positive correlation with BI in the SQCC cohort may be attributable to an indirect causal relationship between cigarette smoking and lipid metabolism alteration mediated by oncogenesis.
Among the lipid peaks that showed a significant positive correlation with BI, the high intensity of one lipid peak (m/z 730.5381) in the ADC cohort was significantly associated with shorter RFP. The PC cluster (m/z 730.5381) in the ADC cohort was revealed to consist of three PC isomers, PC (14:0_18:2), PC (16:1_16:1), and PC (16:0_16:2), and the high-intensity level of this PC cluster was significantly associated with progressed T-factor (≥ pT2) and pleural invasion. According to the 8th edition of the TNM classification for lung cancer, the pathological T-factor is determined by invasion diameter . However, the reproducibility of invasion diameter measurements and the diagnostic concordance among pathologists are low, partially because of their subjective judgement [11, 34, 35]. The routine evaluation method for pleural invasion also varies in pathologists ; this may hinder objective validation. In contrast, the PC cluster identified in our study reflected the likelihoods, displayed as ORs, of these two pathological prognostic factors; these ORs were considered more objective than conventional pathological evaluation. If the prognostic value of the PC cluster is validated in other cohorts, these lipids could be developed as an alternative parameter for the conventional pathological prognostic factors. Although the higher PC cluster showed tendencies of a moderate likelihood for recurrence, differentiation (≥ G2), and STAS, the ORs were not statistically significant; this may be attributable to the small sample size of the cohort.
The species in the PC cluster included 16-carbon fatty acids (FA) of FA (16:0) (palmitic acid) and FA (16:1) (palmitoleic acid). De novo synthesis of FA (16:0) is enhanced by increased expression of fatty acid synthase (FASN) in various cancers . FA (16:0) is desaturated by stearoyl-CoA desaturases (SCDs) to form FA (16:1) . Notably, PC species that contain FA (16:1) are increased in some cancer tissues compared with normal tissues, including lung ADC [37, 38]. Because FASN and SCDs that produce these saturated and unsaturated FAs are suggested to play essential roles in cancer progression [14, 39, 40], our findings that higher PC species, including FA (16:0) and FA (16:1), reflected shorter RFP and disease progression are compatible with these previous reports.
The lower intensity of one lipid peak (m/z 736.5276) in the SQCC cohort was strongly associated with shorter RFP and postoperative recurrence. Nevertheless, associations between m/z 736.5276 and pathological prognostic factors remained weak. Therefore, m/z 736.5276 may reflect postoperative recurrence risk in which conventional pathological prognostic factors cannot be explanatory variables. The LipidSearch™ software suggested that the candidate molecule for m/z 736.5276 is PC (O-14:1_20:5). To the best of our knowledge, the bioactivity of PC (O-14:1_20:5) in cancer biology is unknown. We believe that m/z 736.5276 can be developed as an objective predictor for postoperative recurrence.
Lastly, our previous studies reported SM (d35:1) and SM (t34:1) in NSCLC as candidate predictors for postoperative recurrence after radical surgery [21, 22]. However, we did not identify SM species that reflect BI and postoperative recurrence risk in this study. The volcano plots described on lipid peaks for which SM species were assigned as candidate molecules by LipidSearch™ software showed only one significant lipid peak (m/z 819.6950 [ID: 2132]) in the smoker group of the ADC cohort; the assigned candidate molecule was SM (t41:0) (Additional file 1, Supplemental Fig. 3). This lipid peak did not show significant correlation with BI in the subsequent correlation analysis (Fig. 3a). The candidate lipid biomarkers identified in this study were considered distinctive from the two SM species we have reported previously in that these lipids can reflect both BI and recurrence risk.
Several limitations in this study must be acknowledged. First, our study was conducted with a retrospective design and a small sample size. Therefore, reproducibility could not be validated using other cohorts. On screening of the m/z 730.5381 (ID: 944) and m/z 736.5276 (ID: 1033) using t-test in the volcano plots, the post hoc powers provided by sample sizes of ADC and SQCC cohorts were low as 39.8% and 56.8%, respectively: the low post hoc powers were attributable to the small sample size. Based on standard deviations of m/z 730.5381 (ID: 944) and m/z 736.5276 (ID: 1033) observed in this study, the required sample sizes for a detection power of 80% with a two-sided significance level of P < 0.05 were 30 and 16 patients per group for ADC and SQCC cohorts. A future large cohort validation study is necessary to develop the lipids identified in this study into robust prognostic factors. Second, the MS peak of the PC cluster composed of three different PC isomers could not be isolated into respective isomers because their RTs were too close to one another. Therefore, we could not calculate the area value of single PC isomers, and evaluating the prognostic impact of single PC isomers was difficult. The difference in biological activities among the three PC species is currently unknown. Accordingly, there is no obvious necessity to isolate these PC isomers to evaluate the prognostic impact. Third, structural identification of m/z 736.5276 could not be achieved because of its weak intensity. Hence, it remained a candidate molecule with database exploration. The weak intensity may be attributed to the small amount of extracted lipid samples. Further experiments using specimens from another cohort are needed for molecule identification. Fourth, in the screening of the lipid peaks associated with postoperative recurrence using the KM method for RFP analysis, 10 (38%) patients in the ADC cohort and one (6%) patient in the SQC cohort received adjuvant chemotherapy that was previously shown to reduce the postoperative recurrence risk [41, 42]. Accordingly, the bias of the adjuvant chemotherapy could not be eliminated from the RFP analysis. Fifth, In the ADC cohort, more female patients were enrolled in the never smoker group (75%) than in the smoker group (5.9%); thereby, the difference demonstrated in the volcano plot of the ADC cohort may have been biased by the sex difference. However, we performed further correlation analysis for screening lipids showing a significant correlation with BI. Therefore, we regarded the screened lipids as influenced by cigarette smoking. Further validation study using another cohort with no background difference is expected. Sixth, LC–MS/MS is not standard in routine clinical examinations. However, LC–MS/MS as an optional diagnostic modality in cancer screening has advanced recently [43, 44]. Introducing LC–MS/MS for screening patients with a high risk of postoperative recurrence is strongly expected.
From our data, we propose three PC isomers, PC (14:0_18:2), PC (16:1_16:1), and PC (16:0_16:2), and a lipid peak of m/z 736.5276 as novel candidate biomarkers for postoperative recurrence risk in lung ADC and SQCC patients that smoke. Our findings may contribute to developing an objective method for evaluating postoperative recurrence risk and a qualified postoperative therapeutic strategy for smoking NSCLC patients.
Availability of data and materials
The dataset supporting the conclusions of this article is included within the Additional Files.
Area under the ROC curve
Cytokeratin 19 fragment
Liquid chromatography-tandem mass spectrometry
- m/z :
Median recurrence-free period
Non-small cell lung cancer
Receiver operating characteristic
Reactive oxygen species
Spearman’s rank correlation coefficient
Squamous cell carcinoma antigen
Squamous cell carcinoma
Spread through air space
Tumor, lymph node, and metastasis
Thyroid transcription factor-1
Alberg AJ, Samet JM. Epidemiology of lung cancer. Chest. 2003;123(1 Suppl):21s–49s.
Howington JA, Blum MG, Chang AC, Balekian AA, Murthy SC. Treatment of stage I and II non-small cell lung cancer: Diagnosis and management of lung cancer, 3rd ed: American College of Chest Physicians evidence-based clinical practice guidelines. Chest. 2013;143(5 Suppl):e278S-e313S.
Sardari Nia P, Weyler J, Colpaert C, Vermeulen P, Van Marck E, Van Schil P. Prognostic value of smoking status in operated non-small cell lung cancer. Lung Cancer. 2005;47(3):351–9.
Goldstraw P, Chansky K, Crowley J, Rami-Porta R, Asamura H, Eberhardt WE, et al. The IASLC Lung Cancer Staging Project: Proposals for Revision of the TNM Stage Groupings in the Forthcoming (Eighth) Edition of the TNM Classification for Lung Cancer. J Thorac Oncol. 2016;11(1):39–51.
Nakao M, Yoshida J, Ishii G, Kawase A, Maeda R, Aokage K, et al. Prognostic impact of node involvement pattern in pN1 non-small cell lung cancer patients. J Thorac Oncol. 2010;5(10):1576–82.
Kawase A, Yoshida J, Miyaoka E, Asamura H, Fujii Y, Nakanishi Y, et al. Visceral pleural invasion classification in non-small-cell lung cancer in the 7th edition of the tumor, node, metastasis classification for lung cancer: validation analysis based on a large-scale nationwide database. J Thorac Oncol. 2013;8(5):606–11.
Funai K, Sugimura H, Morita T, Shundo Y, Shimizu K, Shiiya N. Lymphatic vessel invasion is a significant prognostic indicator in stage IA lung adenocarcinoma. Ann Surg Oncol. 2011;18(10):2968–72.
Wang J, Chen J, Chen X, Wang B, Li K, Bi J. Blood vessel invasion as a strong independent prognostic indicator in non-small cell lung cancer: a systematic review and meta-analysis. PLoS One. 2011;6(12):e28844.
Samejima J, Yokose T, Ito H, Nakayama H, Nagashima T, Suzuki M, et al. Prognostic significance of blood and lymphatic vessel invasion in pathological stage IA lung adenocarcinoma in the 8th edition of the TNM classification. Lung Cancer. 2019;137:144–8.
Lu S, Tan KS, Kadota K, Eguchi T, Bains S, Rekhtman N, et al. Spread through Air Spaces (STAS) Is an independent predictor of recurrence and lung cancer-specific death in squamous cell carcinoma. J Thorac Oncol. 2017;12(2):223–34.
Hamanaka R, Yokose T, Sakuma Y, Tsuboi M, Ito H, Nakayama H, et al. Prognostic impact of vascular invasion and standardization of its evaluation in stage I non-small cell lung cancer. Diagn Pathol. 2015;10:17.
Jianyong Z, Yanruo H, Xiaoju T, Yiping W, Fengming L. Roles of lipid profiles in human non-small cell lung cancer. Technol Cancer Res Treat. 2021;20:15330338211041472.
Beloribi-Djefaflia S, Vasseur S, Guillaumond F. Lipid metabolic reprogramming in cancer cells. Oncogenesis. 2016;5(1):e189.
Santos CR, Schulze A. Lipid metabolism in cancer. Febs j. 2012;279(15):2610–23.
Hosokawa Y, Masaki N, Takei S, Horikawa M, Matsushita S, Sugiyama E, et al. Recurrent triple-negative breast cancer (TNBC) tissues contain a higher amount of phosphatidylcholine (32:1) than non-recurrent TNBC tissues. PLoS One. 2017;12(8):e0183724.
Tamura K, Horikawa M, Sato S, Miyake H, Setou M. Discovery of lipid biomarkers correlated with disease progression in clear cell renal cell carcinoma using desorption electrospray ionization imaging mass spectrometry. Oncotarget. 2019;10(18):1688–703.
Eggers LF, Müller J, Marella C, Scholz V, Watz H, Kugler C, et al. Lipidomes of lung cancer and tumour-free lung tissues reveal distinct molecular signatures for cancer differentiation, age, inflammation, and pulmonary emphysema. Sci Rep. 2017;7(1):11087.
Marien E, Meister M, Muley T, Fieuws S, Bordel S, Derua R, et al. Non-small cell lung cancer is characterized by dramatic changes in phospholipid profiles. Int J Cancer. 2015;137(7):1539–48.
Li T, He J, Mao X, Bi Y, Luo Z, Guo C, et al. In situ biomarker discovery and label-free molecular histopathological diagnosis of lung cancer by ambient mass spectrometry imaging. Sci Rep. 2015;5:14089.
Zhang M, He J, Li T, Hu H, Li X, Xing H, et al. Accurate Classification of Non-small Cell Lung Cancer (NSCLC) Pathology and Mapping of EGFR Mutation Spatial Distribution by Ambient Mass Spectrometry Imaging. Front Oncol. 2019;9:804.
Takanashi Y, Funai K, Sato S, Kawase A, Tao H, Takahashi Y, et al. Sphingomyelin(d35:1) as a novel predictor for lung adenocarcinoma recurrence after a radical surgery: a case-control study. BMC Cancer. 2020;20(1):800.
Takanashi Y, Funai K, Eto F, Mizuno K, Kawase A, Tao H, et al. Decreased sphingomyelin (t34:1) is a candidate predictor for lung squamous cell carcinoma recurrence after radical surgery: a case-control study. BMC Cancer. 2021;21(1):1232.
Agarwal AR, Yin F, Cadenas E. Short-term cigarette smoke exposure leads to metabolic alterations in lung alveolar cells. Am J Respir Cell Mol Biol. 2014;51(2):284–93.
Oulton M, Moores HK, Scott JE, Janigan DT, Hajela R. Effects of smoke inhalation on surfactant phospholipids and phospholipase A2 activity in the mouse lung. Am J Pathol. 1991;138(1):195–202.
Mancini NM, Béné MC, Gérard H, Chabot F, Faure G, Polu JM, et al. Early effects of short-time cigarette smoking on the human lung: a study of bronchoalveolar lavage fluids. Lung. 1993;171(5):277–91.
Liu Y, Tian Y, Cai W, Guo Y, Xue C, Wang J. DHA/EPA-Enriched Phosphatidylcholine Suppresses Tumor Growth and Metastasis via Activating Peroxisome Proliferator-Activated Receptor γ in Lewis Lung Cancer Mice. J Agric Food Chem. 2021;69(2):676–85.
Yin H, Liu Y, Yue H, Tian Y, Dong P, Xue C, et al. DHA- and EPA-Enriched Phosphatidylcholine Suppress Human Lung Carcinoma 95D Cells Metastasis via Activating the Peroxisome Proliferator-Activated Receptor γ. Nutrients. 2022;14(21):4675.
Tumours of the lung. In: WHO Classification of Tumours Editorial Board. Thoracic tumours: International Agency for Research on Cancer; 2021. p. 20-21.
Johannsen A, Susin C, Gustafsson A. Smoking and inflammation: evidence for a synergistic role in chronic disease. Periodontol 2000. 2014;64(1):111–26.
Jin C, Yuan P. Implications of lipid droplets in lung cancer: associations with drug resistance. Oncol Lett. 2020;20(3):2091–104.
Jarc E, Petan T. Lipid droplets and the management of cellular stress. Yale J Biol Med. 2019;92(3):435–52.
Herms A, Bosch M, Ariotti N, Reddy BJ, Fajardo A, Fernández-Vidal A, et al. Cell-to-cell heterogeneity in lipid droplets suggests a mechanism to reduce lipotoxicity. Curr Biol. 2013;23(15):1489–96.
Liu L, Zhang K, Sandoval H, Yamamoto S, Jaiswal M, Sanz E, et al. Glial lipid droplets and ROS induced by mitochondrial defects promote neurodegeneration. Cell. 2015;160(1–2):177–90.
Urer HN, Ahiskali R, Arda N, Batur S, Cınel L, Dekan G, et al. Interobserver agreement among histological patterns and diagnosis in lung adenocarcinomas. Turk Patoloji Derg. 2014;30(2):105–10.
Boland JM, Froemming AT, Wampfler JA, Maldonado F, Peikert T, Hyland C, et al. Adenocarcinoma in situ, minimally invasive adenocarcinoma, and invasive pulmonary adenocarcinoma–analysis of interobserver agreement, survival, radiographic characteristics, and gross pathology in 296 nodules. Hum Pathol. 2016;51:41–50.
Travis WD, Brambilla E, Rami-Porta R, Vallières E, Tsuboi M, Rusch V, et al. Visceral pleural invasion: pathologic criteria and use of elastic stains: proposal for the 7th edition of the TNM classification for lung cancer. J Thorac Oncol. 2008;3(12):1384–90.
Muranishi Y, Sato T, Ito S, Satoh J, Yoshizawa A, Tamari S, et al. The Ratios of monounsaturated to saturated phosphatidylcholines in lung adenocarcinoma microenvironment analyzed by Liquid Chromatography-Mass spectrometry and imaging Mass spectrometry. Sci Rep. 2019;9(1):8916.
Dickinson A, Saraswat M, Joenväärä S, Agarwal R, Jyllikoski D, Wilkman T, et al. Mass spectrometry-based lipidomics of oral squamous cell carcinoma tissue reveals aberrant cholesterol and glycerophospholipid metabolism - A Pilot study. Transl Oncol. 2020;13(10):100807.
Luyimbazi D, Akcakanat A, McAuliffe PF, Zhang L, Singh G, Gonzalez-Angulo AM, et al. Rapamycin regulates stearoyl CoA desaturase 1 expression in breast cancer. Mol Cancer Ther. 2010;9(10):2770–84.
Holder AM, Gonzalez-Angulo AM, Chen H, Akcakanat A, Do KA, Fraser Symmans W, et al. High stearoyl-CoA desaturase 1 expression is associated with shorter survival in breast cancer patients. Breast Cancer Res Treat. 2013;137(1):319–27.
Hamada C, Tanaka F, Ohta M, Fujimura S, Kodama K, Imaizumi M, et al. Meta-analysis of postoperative adjuvant chemotherapy with tegafur-uracil in non-small-cell lung cancer. J Clin Oncol. 2005;23(22):4999–5006.
Pignon JP, Tribodet H, Scagliotti GV, Douillard JY, Shepherd FA, Stephens RJ, et al. Lung adjuvant cisplatin evaluation: a pooled analysis by the LACE Collaborative Group. J Clin Oncol. 2008;26(21):3552–9.
Shimbo K, Kubo S, Harada Y, Oonuki T, Yokokura T, Yoshida H, et al. Automated precolumn derivatization system for analyzing physiological amino acids by liquid chromatography/mass spectrometry. Biomed Chromatogr. 2010;24(7):683–91.
Mikami H, Kimura O, Yamamoto H, Kikuchi S, Nakamura Y, Ando T, et al. A multicentre clinical validation of AminoIndex Cancer Screening (AICS). Sci Rep. 2019;9(1):13831.
We thank J. Iacona, Ph.D., from Edanz (https://jp.edanz.com/ac) for editing a draft of this manuscript.
This work used research equipment shared in MEXT Project for promoting public utilization of advanced research infrastructure (Program for supporting the construction of core facilities) Grant Number JPMXS0450200221 and supported by grants from the Smoking Research Foundation [Grant Number 108]. The funding bodies played no role in the study design, analysis, interpretation of data, or writing of the manuscript.
Ethics approval and consent to participate
This study was approved by the Ethics Committee of the Hamamatsu University School of Medicine, Hamamatsu, Japan (#22–210). Written informed consent was obtained preoperatively from patients scheduled for tissue collection.
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Takanashi, Y., Kahyo, T., Hayakawa, T. et al. Lipid biomarkers that reflect postoperative recurrence risk in lung cancer patients who smoke: a case–control study. Lipids Health Dis 22, 15 (2023). https://doi.org/10.1186/s12944-023-01778-3
- Lung cancer
- Cigarette smoking
- Recurrence risk
- Lipid biomarker
- Mass spectrometry