Controlling glycolysis to generate characteristic volatile organic compounds of lung cancer cells

Characteristic volatile organic compounds (VOCs) are anticipated to be used for the identification of lung cancer cells. However, to date, consistent biomarkers of VOCs in lung cancer cells have not been obtained through direct comparison between cancer and healthy groups. In this study, we regulated the glycolysis, a common metabolic process in cancer cells, and employed solid phase microextraction gas chromatography mass spectrometry (SPME–GC–MS) combined with untargeted analysis to identify the characteristic VOCs shared by cancer cells. The VOCs released by three types of lung cancer cells (A549, PC-9, NCI-H460) and one normal lung epithelial cell (BEAS-2B) were detected using SPME–GC–MS, both in their resting state and after treatment with glycolysis inhibitors (2-Deoxy-d-glucose, 2-DG/3-Bromopyruvic acid, 3-BrPA). Untargeted analysis methods were employed to compare the VOC profiles between each type of cancer cell and normal cells before and after glycolysis regulation. Our findings revealed that compared to normal cells, the three types of lung cancer cells exhibited three common differential VOCs in their resting state: ethyl propionate, acetoin, and 3-decen-5-one. Furthermore, under glycolysis control, a single common differential VOC—acetoin was identified. Notably, acetoin levels increased by 2.60–3.29-fold in all three lung cancer cell lines upon the application of glycolysis inhibitors while remaining relatively stable in normal cells. To further elucidate the formation mechanism of acetoin, we investigated its production by blocking glutaminolysis. This interdisciplinary approach combining metabolic biochemistry with MS analysis through interventional synthetic VOCs holds great potential for revolutionizing the identification of lung cancer cells and paving the way for novel cytological examination techniques.

within the headspace of H226 cells.Similarly, when comparing two adenocarcinoma cell lines A549 and HCC827, A549 released more benzaldehyde whereas HCC827 showed higher toluene emissions 9 .
Different lung cancer cells produce distinct VOCs, which can be attributed to variations in genetic mutations among different cell lines 10 .For instance, Peled et al. employed the SPME-GC-MS technique and identified benzaldehyde, triethylamine/decanal, decanal and triethylamine as characteristic VOCs for V-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog (KRAS) mutation, epidermal growth factor receptor (EGFR) mutation, fusion of the EML4 gene to the ALK gene (EML4-ALK) mutation and non-mutant groups respectively 11 .Furthermore, Davies et al. using the thermal desorption-gas chromatography-mass spectrometry (TD-GC-MS) technique observed an increase in 2-methyl-pentane while a decrease in 2,2,3-trimethyl-pentane emissions from the human bronchial epithelial cell line 3KT upon TP53 gene knockout 12 .
Despite the presence of diverse gene mutations in lung cancer cells, these cells predominantly rely on glycolysis, a common biochemical process known as the Warburg effect, for their energy supply 13 .Exploiting this phenomenon, positron emission tomography-computed tomography (PET-CT) has been employed for clinical cancer diagnosis 14 .Therefore, this study proposes a novel approach to induce the production of VOCs with common features in lung cancer cells by regulating the glycolysis process.To this end, three types of lung cancer cells (A549, PC-9, NCI-H460) and one normal lung epithelial cell (BEAS-2B) were selected.SPME-GC-MS technique was used to analyze the VOCs released by these cells in both resting and glycolysis inhibitor-treated states (2-Deoxy-d-glucose, 2-DG/3-Bromopyruvic acid, 3-BrPA).Furthermore, potential biochemical sources of characteristic VOCs were investigated.The experimental process and schematic diagram are presented in Fig. S1.This study aims to induce specific VOCs production in lung cancer cells through manipulation of the glycolysis process with the objective of utilizing these VOCs as the potential biomarker for distinguishing cancer cells from normal cells in vitro.This approach may provide a novel method for developing cytological examination techniques.

Cell lines
The lung cancer cell lines A549, PC-9, NCI-H460, and the normal lung epithelial cell line BEAS-2B used in this experiment were purchased from the American Type Culture Collection (ATCC) and cryopreserved in our laboratory's liquid nitrogen tank at a temperature of − 192 °C.All cells underwent Short Tandem Repeat (STR) identification and mycoplasma detection.

Cell culture
A549, PC-9, NCI-H460, and BEAS-2B cells were cultured in RPMI-1640 medium supplemented with 10% FBS and 1% penicillin/streptomycin at 37 °C under a 5% CO 2 atmosphere in an incubator.The culture medium was refreshed every other day.To minimize the potential influence of plastic products on cell volatiles, the cells were transferred to T-25 cm 3 glass culture flask 15 at the logarithmic growth stage.Each glass flask was inoculated with 8 × 10 5 cells, and the fresh culture medium containing different concentrations of inhibitors was replaced the next day.3-6 parallel samples were set for each group of experiments.To prevent cross-contamination of VOCs released by different cell types within the incubator, parafilm membrane (p-m996; Bemis, USA) was immediately applied after completing the aforementioned procedures.The performance of parafilm membrane has been validated in our previous studies and shown not to interfere with our experimental results 15 .After culturing for 24 h at 37 °C under a 5% CO 2 atmosphere in an incubator, detection was performed using SPME-GC-MS.

Cell viability evaluation
The cell viability was assessed using the CCK8 assay according to the manufacturer's instructions in order to determine the optimal concentration of inhibitors.Take 2-DG as an example, cells were seeded at a density of 5 × 10 3 cells per well in a 96-well plate and incubated for 24 h.Subsequently, the cells were treated with different final concentrations (0, 5, 10, 20, 30, 40 mmol/L) of 2-DG for another 24 h.After removing the supernatant, each well was supplemented with 10 μL of CCK8 reagent.Following a further incubation period of 1.5 h with CCK8 reagent, absorbance at a wavelength of 450 nm was measured using Varioskan Flash microplate reader (Thermo Fisher Scientific, Rockford, IL, USA).The results are presented as percentage (%) cell viability relative to control conditions.The optimum concentration of 3-BrPA was determined using the same methodology.Low concentrations of PA that did not affect cell viability were employed in this study.

SPME process and GC-MS condition
SPME, as a pre-treatment method for the extraction of volatiles, enables the combination of sampling and preconcentration in a single step, thereby simplifying the extraction process 16,17 .The selection of SPME fiber type, extraction time, and extraction temperature are shown in Figs.S2, S3.Similar with our previous research, we selected an SPME fiber (65 μm PDMS/DVB) based on its ability to enrich VOCs and extraction efficiency 15,18 .Prior to each SPME procedure, the extracted fiber was subjected to a 30-min pre-treatment at the 200 °C injection port of GC-MS to eliminate potential residual substances that could interfere with subsequent analysis of cell volatiles.Subsequently, the SPME fiber was introduced into the headspace of the cell culture flask for VOCs extraction at 37 °C incubator for 20 min.Throughout the experiment, only authorized personnel were allowed access to maintain air composition stability within the laboratory environment.SPME-GC-MS detection and analysis were carried out on the gas chromatography triple quadrupole mass spectrometer (TSQ Quantum XLS, Thermo Fisher, USA).The column had a length of 30 m, an inner diameter of 0.32 mm, and a film thickness of 1.8 μm (Thermo, Bellefonte, USA), with the fixed-phase materials consisting of 6% cyanopropylphenyl and 94% polysiloxane.Initially, the column temperature was set at 50 °C for 0.5 min before being increased at a rate of 10 °C/min to reach 180 °C; it was then maintained for an additional 2 min before further elevation to 200 °C at a rate of 15 °C/min and subsequent maintenance for another 5 min.The flow rate of carrier gas helium (99.999%) was 1.5 mL/min with splitless injection mode.The MS conditions were consistent with our previous work 15 : electron impact energy set at 70 eV; mass scanning range from 45 to 200 amu; and utilization of the first quadruple rod as ion filter.Perfluorotributamine, MASS Grade (SynQuest Laboratories USA) was utilized as a corrective fluid to calibrate mass spectrometer.The sensitivity was checked prior to each test and its consistent performance is depicted in Fig. S4a.The relative intensity fluctuation is 7.39%, which proves the stability of the instrument.

Statistical analysis
After GC-MS detection, the Xcalibur 2.2 (Thermo, USA) software was utilized to export all sample chromatographic peak areas.The selected chromatographic peaks were confirmed by NIST database to correct the shift of chromatographic peaks caused by different experiments.Chromatographic peaks with areas exceeding 1 × 10 5 were selected and labeled as VOC1-VOC135 in ascending order of their retention time.Subsequently, orthogonal partial least squares discriminant analysis (OPLS-DA), the Mann-Whitney U test, and fold change (FC: ratio of lung cancer cell chromatographic peak area value to normal cell chromatographic peak value) were applied to analyze the 135 chromatographic peak areas of three types of lung cancer cells and one normal cell under resting state and glycolysis inhibition conditions.A substance corresponding to a chromatographic peak was considered significantly variable if it met the following criteria: variable importance projection (VIP) > 1, p-value < 0.01 (the p-value of the Mann-Whitney U test, p < 0.01, Benjamini-Hochberg corrected), and FC > 2. The qualitative analysis proceeded as follows: initially, potential substance candidates were selected based on a Reverse Search Index (RSI) value exceeding 800 from the NIST database.Subsequently, the retention index (RI) for certain VOC was determined using the experimental retention times.0][21] , the most probable VOCs were ultimately identified from the substance candidates provided by the NIST database.Finally, a portion of the VOCs were definitively identified through the use of chemical standards.Table S1 shows the specific information and qualitative parameters of differential VOCs.The OPLS-DA analysis was carried out using SIMCA P14.1 while the Mann-Whitney U test was performed using IBM SPSS Statistics 25.

Determination of the concentration range of inhibitors
In order to determine the appropriate concentration range of a glycolysis inhibitor, we investigated the relationship between the concentration of 2-DG and cell survival.The experimental results are presented in Fig. 1.Following treatment with varying concentrations of 2-DG, the survival rate of four different cell types (A549, PC-9, NCI-H460, BEAS-2B) exhibited a dose-dependent decrease with increasing inhibitor concentration.Specifically, as the concentration of 2-DG increased from 0 to 40 mmol/L, the survival rate decreased from 100 to 40% for A549 cells, from 100 to 63% for PC-9 cells, from 100 to 47% for NCI-H460 cells, and from 100 to 47% for BEAS-2B cells.Notably, BEAS-2B cells demonstrated higher sensitivity towards low concentrations of 2-DG compared to other cell types; conversely, A549 cells displayed greater sensitivity towards high concentrations of 2-DG.
To ensure that the characteristic VOCs are from living cells, according to the principle that the cell survival rate is nearly 80%, we finally chose the concentration of the experimental inhibitor as 2-DG of 10 mmol/L.

Differences in VOCs between lung cancer cells and normal cells in resting state
In the resting state, the headspace VOCs of A549, PC-9, NCI-H460, and BEAS-2B were analyzed using the SPME-GC-MS technique.OPLS-DA was conducted utilizing the VOCs released by three distinct types of lung cancer cells, as well as normal cells during resting state.The resulting score plots are presented in Fig. 2a-c, where it is evident that BEAS-2B cells exhibit a clear separation from the three lung cancer cell types.Notably, these separations exhibit high goodness of fitting (R 2 Y) and prediction (Q 2 ), indicating the robustness of the model.To further validate the efficacy of the OPLS-DA approach, 200 permutation tests were performed, and the results are depicted in Fig. 2d-f.The R 2 and Q 2 values obtained from the actual model exceed all corresponding values in the permutation models, unequivocally demonstrating the validity and statistical significance of the OPLS-DA analysis conducted in this study.
Following OPLS-DA, Mann-Whitney U test, and FC analysis, the results revealed distinct VOC profiles between each type of lung cancer cell and normal cells.Initially, based on the U-test analysis and FC a Volcano   plot (Fig. S5a-c) was generated to visualize the changes in VOCs between each lung cell lines and normal cell line during the resting state.Subsequently, with reference to the OPLS-DA VIP > 1 and NIST database RSI > 800 differential VOCs were summarized in Table 1.As depicted in Venn diagram of difference VOCs (Fig. S6a), three common substances were found: VOC19 (ethyl propionate), VOC23 (acetoin), and VOC80 (3-decen-5one).Ethyl propionate and 3-decen-5-one exhibited lower levels in the headspace of each lung cancer cell type compared to normal cells; conversely, acetoin showed higher levels.In conclusion, these three aforementioned VOCs may serve as potential biomarkers for distinguishing lung cancer cells from normal cells.The Volcano plot (Fig. S5d-f) serves as a visualization tool to illustrate the alterations in VOCs between each lung cancer cell and normal cells after 2-DG regulation.Table 2 presents the differential VOCs information.Furthermore, the Venn diagram (Fig. S6b) provides a concise representation of the overlap among the three sets of differential VOCs, revealing that only one substance, VOC23 (acetoin), was consistently detected across all three groups.All three types of lung cancer cells exhibited higher levels of acetoin in the headspace compared to normal cells.

Identification of characteristic VOC for lung cancer cells under glycolysis regulation
It is worth emphasizing that in comparison to the resting state, the glycolysis inhibition by 2-DG significantly augmented the disparity in acetoin levels in the headspace of lung cancer cells (A549, PC-9, NCI-H460) and normal cells (BEAS-2B).The FC values increased from 5.67, 6.29, 8.10 times to 12.93, 18.74, and 20.56 times respectively (Fig. 4).To further validate the association between increased acetoin levels and glycolysis inhibition, we conducted a repetitive experiment using another glycolysis inhibitor (3-BrPA).Initially, an appropriate concentration range of 100 μmol/L was selected for 3-BrPA (Fig. S7).Subsequently, SPME-GC-MS detection and statistical analysis were performed on the headspace VOCs of A549, PC-9, NCI-H460, and BEAS-2B under glycolysis regulation with 3-BrPA treatment.Relevant results are presented in Fig. S8, Fig. S5h-j, Table S2, and Fig. S6c.Importantly, acetoin was still detected among all three groups of differential VOCs.The FC values for lung cancer cells (A549, PC-9, NCI-H460) and normal cells (BEAS-2B) increased from 5.67, 6.29, 8.10 times to 13.85, 16.76, 21.97 times, respectively (Fig. 4).Overall, glycolysis inhibition can significantly amplify the disparity in acetoin levels between lung cancer cells and normal cells.This suggests that acetoin could serve as the characteristic VOC for distinguishing between lung cancer cells and normal cells under glycolysis regulation.We also noticed that Feinberg et al. 22 mentioned that the substance with m/z of 49 experienced the most significant reduction in A549 cells after the action of inhibitor 3-BrPA, which speculated to be methanethiol.In our experiment, no characteristic VOC with m/z of 49 was detected, including methanethiol, and the reason for the two differences is not clear.

Biochemical sources of characteristic VOC for lung cancer cells under glycolysis regulation
Several studies have observed the presence of acetoin [23][24][25] in the breath of individuals with lung cancer, Sanni et al. 25 have identified that this compound may originated from the glucose metabolism of oral bacteria.To our knowledge, there is currently no available research on the biochemical origins of acetoin in human cells.However, it has been confirmed that three distinct pathways exist for the production of acetoin in bacteria and yeast, which are schematically represented in Fig. S9 [26][27][28] .Pyruvate serves as the starting substrate to produce acetoin based on these three pathways.Glycolysis is a process by which glucose is broken down into pyruvate in the cytoplasm, releasing energy.When glycolysis is inhibited, mitochondria increase their utilization of glutamine to produce organic molecules such as pyruvate and acetyl CoA to maintain energy and biosynthesis needs 29 , as depicted in Fig. 5.
Therefore, we hypothesize that the glycolysis inhibition leads to an upregulation of the glutamine degradation pathway, resulting in a compensatory increase in acetoin production, as illustrated in Fig. 6a.Consequently, it is reasonable to consider whether concurrent restriction of the glutamine breakdown pathway and glycolysis inhibition would lead to a reduction in acetoin (Fig. 6b).
To investigate the potential of glycolysis inhibition in enhancing glutaminolysis for pyruvate production in lung cancer cells, we used 2-DG/3-BrPA to inhibit glycolysis while utilizing PA to block glutaminolysis 30 .As depicted in Fig. 7, compared to the resting state, treatment with 2-DG/3-BrPA significantly increased acetoin levels in the headspace of all three types of lung cancer cells.However, there was no significant alteration observed when treated with PA alone.Moreover, concurrent administration of 2-DG/3-BrPA and PA effectively suppressed the upward trend of acetoin release.These findings further validate our hypothesis (Fig. 6b).It is evident that the elevation in acetoin is associated with compensatory enhancement of glutaminolysis.

Conclusions
To address the challenge of VOC replication in lung cancer cells, we focus on glycolysis, a common metabolic process in cancer cells.Initially, we regulated the metabolic processes of lung cancer cells (A549, PC-9, NCI-H460) and normal lung epithelial cell (BEAS-2B) using the glycolysis inhibitors (2-DG and 3-BrPA).We then analyzed the cell headspace VOCs before and after glycolysis regulation using SPME-GC-MS technique combined with untargeted analysis methods.Our results showed that application of glycolysis inhibitors significantly increased acetoin release by all three types of lung cancer cells but not by normal lung epithelial cells.This observation further highlighted the disparity between lung cancer cells and normal cells under chemical intervention.Additionally, we investigated the mechanism of acetoin production by examining pharmacological effects of glycolysis inhibitors.Our interdisciplinary approach combining biochemical metabolism and mass spectrometry analysis can differentiate between lung cancer and normal cells through intervening in VOC synthesis.These findings may have significant implications for cytological examination and diagnosis.
While this study provides valuable insights into VOCs in lung cancer cells, further research is needed to understand variations in VOCs among cell lines with varying glycolysis levels and to assess its broader applicability.Future work will focus on co-culturing normal and lung cancer cells and detecting VOCs in human tumor tissues under controlled glycolysis conditions, with the goal of translating this research into clinical applications.

After the glycolysis inhibition by 2 -
DG, headspace VOCs of A549, PC-9, NCI-H460, and BEAS-2B cells were analyzed using SPME-GC-MS technique and subjected to statistical analysis.OPLS-DA scores plots and corresponding 200 permutation tests are presented in Fig. 3.All three OPLS-DA models have high goodness of fitting and goodness of prediction.Additionally, the simulated values of 200 replacement tests of the corresponding models are greater than all their original values, indicating the effectiveness of the OPLS-DA model.

Figure 4 .
Figure 4. Differences in the peak area of acetoin in the headspace before and after glycolysis inhibition between lung cancer cells and normal cells.

Figure 6 .
Figure 6.(a) Mechanism underlying the upregulation of acetoin resulting from glycolysis inhibition.(b) Mechanism explaining the downregulation of acetoin due to combined inhibition of glycolysis and blockade of glutaminolysis.The thickness of the arrows represents pathway strength or weakness.

Figure 7 .
Figure 7. Alteration in the peak area of acetoin within the headspace of lung cancer cells when 2-DG/3-BrPA inhibiting glycolysis while PA blocking glutaminolysis.