J Cancer 2022; 13(1):76-87. doi:10.7150/jca.59777 This issue

Research Paper

Metabolomics based plasma biomarkers for diagnosis of oral squamous cell carcinoma and oral erosive lichen planus

Xibo Li1,2, Liwei Liu3,4, Na Li5, Qingquan Jia3,4, Xiaoshuang Wang1,2, Lihua Zuo3,4, Jianglan Long6, Peng Xue7, Zhi Sun3,4 Corresponding address, Hongyu Zhao1,2 Corresponding address

1. Department of Oral Emergency, The First Affiliated Hospital of Zhengzhou University· Stomatological Hospital of Henan Province, Zhengzhou, Henan, 450052, China.
2. School and Hospital of Stomatology of Zhengzhou University, Zhengzhou, Henan, 450052, China.
3. Department of Pharmacy, the First Affiliated Hospital of Zhengzhou University, Zhengzhou, Henan, 450052, China.
4. Henan Engineering Research Center of Clinical Mass Spectrometry for Precision Medicine, Zhengzhou, Henan, 450052, China.
5. Department of Prosthodontics, The First Affiliated Hospital of Zhengzhou University· Stomatological Hospital of Henan Province, Zhengzhou, Henan, 450052, China.
6. Beijing Key Laboratory and Joint Laboratory for International Cooperation of Bio-characteristic Profiling for Evaluation of Rational Drug Use, Capital Medical University Affiliated Beijing Shijitan Hospital, Beijing, 100038, China.
7. Health Management Center, The First Affiliated Hospital of Zhengzhou University· Stomatological Hospital of Henan Province, Zhengzhou, Henan, 450052, China.

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
Citation:
Li X, Liu L, Li N, Jia Q, Wang X, Zuo L, Long J, Xue P, Sun Z, Zhao H. Metabolomics based plasma biomarkers for diagnosis of oral squamous cell carcinoma and oral erosive lichen planus. J Cancer 2022; 13(1):76-87. doi:10.7150/jca.59777. Available from https://www.jcancer.org/v13p0076.htm

File import instruction

Abstract

Graphic abstract

Backgrounds: To identify diagnostic biomarkers for differentiating oral squamous cell carcinoma (OSCC) from oral erosive lichen planus (OELP) and investigate potential biomarkers associated with malignant transformation.

Methods: In this study, 72 patients with OSCC, 75 patients with OELP subjects were recruited. Their plasma samples were analyzed by ultra-high-performance liquid chromatography quadrupole-Orbitrap high-resolution accurate mass spectrometry, (UHPLC/Q-Orbitrap HRMS). Principal component analysis, orthogonal partial least square discrimination analysis, t-test analysis and false discovery rate were used to identify different metabolites in patients with OSCC and OELP. The metabolic pathway analysis was performed by MetaboAnalyst. To further screen and identify the biomarkers of OSCC and establish a diagnostic panel, binary logistic regression analysis and receiver operating characteristic analysis were used. The data were then combined with blood samples from healthy individuals for mass spectrometry analysis to obtain biomarkers related to malignant transformation.

Results: A total of 20 kinds of endogenous metabolites were identified from plasma samples of OSCC patients and OELP patients. Metabolic pathway analysis showed that the biomarkers associated with OSCC were closely related to cholic acid metabolism and amino acid metabolism. Finally, a diagnostic panel composed of decanoylcarnitine, cysteine and cholic acid was established. This diagnostic panel had good diagnostic efficiency with the AUC=0.998. Other metabolites including uridine, taurine, glutamate, citric acid and LysoPC(18:1) were identified to be general biomarkers for malignant transformation of OELP.

Conclusion: Biomarkers based on plasma metabolomics are of great significance for the prediction of malignant transformation of OELP and early diagnosis of OSCC.

Keywords: Oral squamous cell carcinoma, Oral lichen planus, Metabolomics, Biomarker, Plasma

Introduction

Oral squamous cell carcinoma (OSCC) is one of the most prevalent tumors of head and neck and its 5-year survival rate is about 50%-60% [1]. According to the Globocan Project (http://globocan.iarc.fr/Default.aspx) data, there are about 300 000 new cases every year, and 145 000 of them will die [2]. A major cause for the high mortality of OSCC is that there is a lack of effective biomarker for early-stage diagnosis. The high mortality rate of OSCC is largely attributed to the fact that early-stage OSCC is mostly asymptomatic. In addition, the maxillofacial region has abundant blood and lymphatic vessels that facilitate invasion and metastasis of cancer cells. Therefore, most of the patients have advanced stage disease at diagnosis. Oral lichen planus (OLP) is the precancerous change of OSCC. Oral erosive lichen planus (OELP), a subtype of oral lichen planus, has a higher malignant transformation rate than the reticular type. The clinical presentations, especially oral mucosa erosions, are similar to the presentations of OSCC [3]. It is difficult to distinguish early OSCC from OELP. There is an urgent need to distinguish OELP from OSCC to better the diagnosis and treatment.

Metabolomics is one of the typical hallmarks of cancer cells and is closely related to cancer development. Cancer cells maintain rapid proliferation by reprogramming their metabolic mechanisms [4]. In general, metabolomics leads to abnormal levels of differential metabolites in blood, saliva and tissues. The changed metabolites may be potential biomarkers to distinguish malignant and benign lesions. Metabolomics is a high-throughput technique which is to measure the expression levels of small molecular compounds [5], such as lipids [6], amino acids [7], and other small molecule compounds, in body fluids or tissues [8]. Theoretically, detecting metabolic change is a feasible method for diagnosis of OSCC.

Recently, many metabolomic studies have been performed to detect metabolic changes in patients with OSCC and healthy controls. Some study groups [9, 10] successfully revealed metabolites as potential biomarkers and used these biomarkers to discriminate OSCC patients from healthy control. Salivary biomarkers have attracted attention for clinical diagnosis because of the noninvasive sampling method [11, 12]. Nevertheless, there is few studies [13] focused on searching for reliable plasma biomarkers of OSCC and OELP. Blood is more stable and contains more analytes compared with saliva [14]. Several typical tumor specific proteins in blood, such as alpha fetoprotein (AFP) [15], prostate specific antigen (PSA) [16] and carcinoembryonic antigen (CEA) [17], have been used as biomarkers for clinical diagnosis. Our team has done some studies about diagnostic biomarkers on OELP and obtained some promising results [18]. On the basis of previous research, the study aims to identify diagnostic biomarkers from metabolic reprogramming combined with mathematical model analysis, and build a reliable diagnostic panel. We also explored whether there were biomarkers that could predict the malignant transformation of OELP. The study will provide new ideas for the early diagnosis and treatment of OSCC, which has economic value and social significance.

Materials and methods

Sample Information

The study protocol was approved by the Ethics Committee of The First Affiliated Hospital of Zhengzhou University (approval No. 2020-KY-036). Written informed consent was obtained from all participants. All patients were diagnosed clinically and confirmed by pathology. Their blood samples were collected between June 6, 2019 and January 18, 2020. All the subjects were diagnosed as OSCC or OELP for the first time with no serious systemic disease. The OELP criteria for the diagnosis of OELP are the modification WHO diagnostic criteria of oral lichen planus [19]: In short, clinical criteria included presence of bilateral, more or less symmetrical lesions accompanied by erosion. Moreover, the histological criteria included a clearly defined band-like zone of cellular infiltration which is limited to the surface part of connective tissue, mainly composed of lymphocytes, in the basal cell layer, signs of 'liquefactive degeneration' and absence of epithelial dysplasia. Clinical examinations were performed by two chief physicians with more than 30-year clinical experience. Also, every patient with OELP had undergone histopathological examination. In total, 72 patients with OSCC, 75 patients with OELP were recruited in this study. To detect a 50% difference in peak areas of each ion peak with a sample size of 100 in total, we could obtain a power over 0.99. In our case, the differences in peak areas between OSCC and OELP groups are mostly more than 50%. Therefore, a total number of 147 patients is a reasonable sample size to detect the difference. The original dataset was divided randomly into training (n=98) and validation (n=49) sets in a ratio of 2:1, according to random number method. Samples were collected and stored at -80°C before UHPLC-Q-Orbitrap analysis.

Experimental equipment

Thermo Fisher centrifuge (Thermo Fisher Science, USA); VORTEX-GENIE 2 vortex oscillator (SITM, USA); BX7200HP ultrasonic cleaner (Shanghai Xinmiao medical device manufacturing Co., Ltd.); AL104 balance with accuracy of 0.0001 (Mettler Toledo Shanghai Co., Ltd., Switzerland); CF16RN centrifuge (Hitachi, Japan); Ultrapure water meter (MilliPore, USA); and -80 °C ultra-low temperature refrigerator (Thermo Fisher science, USA).

Experimental reagent

Acetonitrile (Chromatographic grade, Fisher, USA); Methanol (Chromatographic grade, Fisher, USA); Formic acid (Chromatographic grade, Shanghai Aladdin Biotechnology Co., Ltd., China); L-2-chlorophenylalanine (Bailingwei Technology Co., Ltd., China); Ketoprofen (Sigma, USA); all solutions were filtered by 0.22 μm aperture filter before used.

Sample collection

All blood samples were collected from 8:00 a.m. to 10:00 a.m. The blood was collected to the vacuum tubes containing coagulant, placed in the incubator containing ice. The sample was centrifuged at 1510×g for 10 minutes at 4°C, and the supernatant was quickly stored at -80°C until used.

Sample preparation

100 μL sample was taken out and placed in a 1.5 mL centrifuge tube after thawing. 300 μL methanol solution containing internal standards (0.05 μg/mL L-2-chlorophenylalanine and 0.5 μg/mL ketoprofen) was added. After vortex oscillation for 1 min, centrifugation was performed at 16 200 × g for 10 min at 4°C. The supernatant was aspirated to the vial for analysis.

QC sample preparation

Quality control sample (QC sample) analysis could ensure the reliability of the experimental results in the process of collecting metabolomics data of all samples. 6 QC samples were detected to monitor the pressure change before and after each injection and the shift of the main peak retention time of the total ion flow diagram. After the instrument was stable, the sample analysis started. A QC sample was inserted into every ten samples to verify the stability of the instrument. Inserted a blank sample containing only solvent after each QC sample to avoid cross-contamination.

UHPLC-Q-Orbitrap analysis

UHPLC-Q-Orbitrap system: Ultimate 3000 UHPLC (Thermo Fisher Scientific, USA), Q Exactive high resolution mass spectrometry (Thermo Fisher Scientific, USA); ACQUITY UHPLC®RBEH C18 (100 × 2.1 mm, 1.7 µm) column (Waters Company, USA).

Metabolomic analyses of plasma samples were as described previously [20]. The ultra-performance high liquid chromatography (UHPLC) system was used to separate the metabolites in plasma samples. 5 µL was extracted from each sample and injected into the ACQUITY UHPLC®RBEH C18 column (100 × 2.1 mm, 1.7 µm), and the column temperature was 40°C. The mobile phase was acetonitrile (A) with 0.1% formic acid aqueous solution (B), and gradient elution flow rate was 0.2 mL/min:0.0~1.0 min, 5%A; 1.0~9.0 min, 5%~100%A; 9.0~12.0 min, 100%A; 12.0~12.1 min, 100%~5%A; 12.1~15.0 min, 5%A.

Heated electrospray ionization (HESI) was combined with high resolution mass spectrometry to UHPLC system. The temperature of auxiliary gas was 300 °C and the flow rate was 10 arb. The ion source and capillary were 350 °C and 320 °C. The detection was performed in positive ion mode and negative ion mode with a resolution of 17 500 in full mass/DDMS2 (data dependent mass spectrometry) scanning mode. The collision energy was set at gradient from 20 eV to 60 eV. The spray voltage and sheath gas flow rate were 3.50 kV and 40 arb in positive ion mode, and 2.80 kV and 38 arb in negative ion mode. The injection sequence of all samples was random.

Data processing and statistical analysis

Data were tested for normality using a Shapiro-Wilk normality test. When the normal distribution was satisfied, an independent-samples t-test was applied. Otherwise, a non-parametric Wilcoxon test was performed. All metabolomics data were analyzed by Thermo Xcalibur™ software (Version 3.0, Thermo Scientific, USA). Specific parameters were as previously [18]. Finally, the generated data and the m /z value, retention time (RT) and peak area of each ion peak in each sample were collected. The peak areas represented the levels of metabolites. The data sets were imported into the multivariate statistical analysis software SIMCA (Version 14.0, Umetrics, Sweden). The principal component analysis (PCA) and orthogonal partial least square discriminant analysis (OPLS-DA) were performed to explore separation trend among groups. Through the establishment of the OPLS-DA model, variable importance in projection (VIP value) was obtained. Two hundred permutation tests used to evaluate whether the data was overfitted. P values were obtained using the independent-samples t-test. Also, in order to further screen the metabolites with significant difference between different groups, the false discovery rate (FDR) which was calculated by R language was conducted for metabolites with VIP value greater than 1 using SPSS 26.0 software (IBM, USA). Eventually, the metabolites were selected for identification when FDR < 0.05. The accurate m / z, ion chromatogram, retention time (RT) and other information were compared with ChemSpider and MassList database. The MS/MS data was compared with mzVault, Human Metabolome Database (HMDB, http://hmdb.ca/) and PubChem compound database. For some endogenous metabolites of which the standard substance could be obtained, the data were compared with the standard substance to determine its structure. When the data matched with the information in database, the metabolite was considered to be identified successfully. Moreover, to screen the metabolites of OELP malignant transformation, fold changes and FDR of the potential biomarkers for comparisons of OELP versus HC and OSCC versus HC were calculated. Using MetaboAnalyst (www.metaboanalyst.ca) platform, the screened differential metabolites were analyzed by thermography to show the change of the metabolites, and receiver operating characteristic (ROC) curve was drawn for each identified differential metabolite. Area under the curve (AUC) was calculated by Medcalc. A metabolic pathway network was formed according to Kyoto Encyclopedia of Genes and Genomes (KEGG, https://www.kegg.jp/kegg/pathway.html) signaling pathway database.

 Figure 1 

Flowchart. OSCC: oral squamous cell carcinoma; OELP: oral erosive lichen planus; UHPLC/Q-Orbitrap HRMS: ultra-high-performance liquid chromatography quadrupole-Orbitrap high-resolution accurate mass spectrometry.

J Cancer Image

(View in new window)

Results

Demographic baseline characteristics

The flowchart of the study was shown in Figure 1. A total of 147 subjects were enrolled in training group including 72 patients with OSCC with a mean age of 66 ± 12 (yrs.; mean ± SD), 75 patients with OELP with a mean age of 61 ± 7 (yrs.; mean ± SD). There were no significant differences in gender, age, BMI and lifestyle habits of participants among these three groups. The data of 48 OSCC patients and 50 OELP patients were used for biomarker discovery and the others for validating the effectiveness of the chosen biomarkers. The demographic baseline characteristics of these individuals were shown in Table 1.

 Table 1 

Demographic baseline characteristics of patients with OSCC, OELP and healthy controls

CharacteristicsOSCC patients (n=72)OELP patients (n=75)Healthy controls (n=47)
Age (yrs.; mean ± SD)66 ± 1261 ± 765 ± 9
Gender
Male353823
Female373724
BMI21.06±2.8121.74±2.4020.83±2.98
Habit of chewing betel nut7/725/724/72
Habit of smoking (≥1 cigarette / d for more than half a year)
Yes32/7231/7518/47
Habit of eating spicy food
Yes43/7242/7529/47
Habit of eating very hot food
Yes15/7218/7514/47
Clinical stage
I17
II21
III19
IV14
Unknown1

BMI: Body mass index; Habit of chewing betel nut: often chewing betel nut; Habit of eating spicy food: often intake spicy food; Habit of eating very hot food: often intake very high-temperature food; OSCC: oral squamous cell carcinoma; OELP: oral erosive lichen planus; SD: standard deviation; confirmed by a histopathological examination.

Primary metabolites analysis in blood samples

To gain insights into the metabolic features of OELP progressed to OSCC, UHPLC/Q-Orbitrap HRMS was performed on blood samples of 48OSCC, 50 OELP and 47 healthy controls (HC). A total of 3238 ion peaks in positive ion modes and 2663 in negative ion modes were extracted. The metabolic dataset was next analyzed by PCA. QC samples cluster tightly together indicating the good instrument stability and the reliability of the data. As shown in Figure 2A-B, the disease groups, including both the OSCC group and the OELP group, were clearly separated from HC, but there was not a sharp distinction between OSCC and OELP group. This result confirmed that it is difficult to distinguish OELP from early OSCC in clinical diagnosis, but the difference between healthy people and disease groups is very clear. Therefore, we further analyzed differential metabolites between OSCC and OELP group.

 Figure 2 

Multivariate statistical analysis revealed distinct metabolic characteristics between OSCC and OELP. HC: red pots; OELP: green pots; OSCC: blue pots; QC: yellow pots. PCA score plots among OSCC, OELP and HC in (A) positive ion mode and (B) negative ion mode built with metabolites identified by LC-MS. The abscissa and the ordinate represented the first and second principal component. PCA score plots of the OSCC and OELP group in (C) positive ion mode and (D) negative ion mode built with metabolites identified by LC-MS. The abscissa represented and the ordinate represented the first and second principal component. OPLS-DA score plots of the OSCC and OELP in (E) positive ion mode and (F) negative ion mode. The model showed a clear trend of separation between OSCC and OELP. Cross-validation plots with a permutation test repeated 200 times of the OSCC and OELP groups in (G) positive ion and (H) negative ion mode. It showed that the OPLS-DA models were not overfitting. PCA: principal component analysis; OPLS-DA: orthogonal partial least square discriminant analysis; HC: healthy control; OSCC: oral squamous cell carcinoma; OELP: oral erosive lichen planus.

J Cancer Image

(View in new window)

Screening and identifying differential metabolites

In order to further explore the unique metabolic characteristics between OSCC and OELP, PCA analysis were carried out. There was a clear trend of inter-group separation between the two groups. All samples were analyzed in both positive (Figure 2C) and negative (Figure 2D) ion modes. To better define the metabolic variations in OSCC and OELP, OPLS-DA analysis were performed. As shown in Fig. 2E-F, green dots represented patients of OELP and blue dots represented patients of OSCC. R2Y = 0.932, Q2 (cum) = 0.870 in positive ion mode and R2Y=0.905, Q2 (cum) = 0.698 in negative ion mode. The results indicated significant separation between groups. Two hundred permutation tests were further performed, with R2 = (0.0, 0.519), Q2= (0.0, -0.450) in positive ion mode and R2 = (0.0, 0.666), Q2= (0.0, -0.503) in negative ion mode, which confirmed that these OPLS-DA models were not overfitting (Figure 2G-H).

The differential metabolites were further screened by combining the P values or the fold change and VIP values of the OPLS-DA model. Volcano plots were drawn using fold change (FC) and P values. Red dots represented metabolites with P < 0.05 (- log10P > 1.30) and FC > 2.0 (log2FC > 1.0). The sites with VIP > 1 and P < 0.05 were regarded as candidate differential metabolites (Figure 3). After comparison with databases, a total of 20 endogenous metabolites between OSCC and OELP were identified. These endogenous metabolites included amino acids such as cysteine, glutamate, phenylalanine; lipids, such as lysophosphatidylcholine (LysoPC) and other small-molecule compounds. The details of these metabolites were given in Table 2. The heatmap of the differential metabolites were shown in Figure 4 which appeared the changes in metabolic signatures among OSCC and OELP. In order to better understand the relationship among metabolites, the metabolic pathway network diagram was shown in Figure 5.

 Table 2 

Metabolites different between OSCC and OELP

IDNameFormulaMolecularRTVIPPFCAUCFDR
1Uric acidC5H4N4O3168.020191.024.3530.0091.1700.6770.013
2GlutamateC5H9NO4147.045881.031.861.48E-061.2790.7803.27E-03
3Acetyl carnitineC9H17NO4203.122861.032.6414.42E-050.7150.7381.89E-04
4LeucineC6H13NO2132.086401.481.0341.95E-071.4380.8102.35E-05
5PhenylalanineC9H11NO2165.085911.953.9493.83E-071.2350.6863.25E-06
6PhenylacetylglutamineC13H16N2O4264.103883.972.3040.0031.9100.7110.017
72,4-DihydroxyacetophenoneC8H8O3152.039124.191.0475.48E-181.3930.7666.72E-13
8CysteineC3H7NO2S121.011384.471.5137.35E-1454.5850.9059.48E-12
9Acetyl phenylalanineC11H13NO3207.082114.621.0036.34E-131.8190.8880.005
104-NitrophenolC6H5NO3139.018655.281.9403.78E-050.1430.8224.16E-11
11Indoxyl sulphateC8H7NO4S213.002145.53.7503.82E-051.6100.7000.013
12OctanoylcarnitineC15H29NO4287.216645.511.6942.33E-050.6160.7631.09E-04
13GlycitinC22H22O10446.129646.081.0426.21E-060.7830.8173.53E-05
14DecanoylcarnitineC17H33NO4315.248326.262.5093.98E-060.4650.9662.39E-05
15Monobutyl phthalateC12H14O4222.081896.521.8307.04E-201.9750.8340.006
16Cholic acidC24H40O5408.281566.83.460.0093.5090.9650.012
179,10-DihydroxyoctadecanoateC18H36O4316.254368.081.1020.0020.8210.6857.40E-09
189,10-Epoxyoctadecenoic acidC18H32O3296.228009.224.2750.0040.7080.6800.003
19Oleoyl ethanolamideC20H39NO2325.304879.571.0192.16E-071.5550.7171.99E-06
205-Octadecynoic acidC18H32O2280.2330310.179.6467.30E-080.6000.8235.45E-04

OSCC, oral squamous cell carcinoma; OELP, oral erosive lichen planus; FC, fold change; RT, retention time; VIP, variable importance in projection; FC, fold change; AUC, area under curve; FDR, false discovery rate, calculation based on Benjamini and Hochberg.

 Figure 3 

The volcano plot of the OSCC vs. OELP in (A) positive ion mode and (B) negative ion mode. OSCC, oral squamous cell carcinoma; OELP, oral erosive lichen planus.

J Cancer Image

(View in new window)

 Figure 4 

Heat map of 20 differential metabolites showing significant differences between the oral squamous cell carcinoma and oral erosive lichen planus (OELP) groups. Red and yellow represented higher and lower concentrations. The abscissa represented the number of subjects. The ordinate represented differential metabolites of the two groups. Pearson's correlation-based clustering analysis was used for clustering analysis of 20 differential metabolites. OSCC, oral squamous cell carcinoma; OELP, oral erosive lichen planus.

J Cancer Image

(View in new window)

 Figure 5 

Metabolic pathway network diagram of differential metabolites in OSCC and OELP. The violet color represents downregulated metabolites, red color represents upregulated metabolites, the green color metabolites represent intermediary metabolites in the metabolic pathway of OSCC. OSCC, oral squamous cell carcinoma; OELP, oral erosive lichen planus.

J Cancer Image

(View in new window)

Pathway analysis

To further explore the underlying molecular mechanism of OSCC, the metabolic pathways of the metabolites were analyzed by MetaboAnalyst (Figure 6). The results showed that amino acid metabolism, including phenylalanine, tyrosine and tryptophan biosynthesis, D-Glutamine and D-glutamate metabolism, phenylalanine metabolism; primary bile acid biosynthesis and arginine biosynthesis were associated with OSCC. Phenylalanine, tyrosine and tryptophan biosynthesis and D-Glutamine and D-glutamate metabolism had a major impact on OSCC.

Establishment of a diagnostic panel

Binary logistic, ROC analysis and VIP were used to evaluate integrated biomarkers. ROC curve was used to assess the diagnostic performance of each metabolite. Subsequently, the metabolites of AUC > 0.900 were used to establish a diagnostic panel. The panel with reasonably high accuracy and sensitivity exhibited well-established performance. All 95% confidence intervals were given in the form: (95%CI lower, upper). Ultimately, decanoylcarnitine (Figure 7A), cysteine (Figure 7B) and cholic acid (Figure 7C) were selected to serve as a useful biomarker panel for the diagnosis of OSCC from OELP. The FC of decanoylcarnitine, cysteine and cholic acid in the observation group were 0.465, 54.585, 3.509, respectively (Table 2), the difference was statistically significant (P < 0.05). ROC analysis showed that the diagnostic capability of biomarker panel (95%CI 0.904, 0.999, P < 0.0001) (Figure 7D) was much higher than those of previous biomarkers of OSCC, including decanoylcarnitine (95%CI 0.841, 0.968, P < 0.0001), cysteine (95%CI 0.933, 0.999, P < 0.0001), and cholic acid (95%CI 0.921, 0.999, P < 0.0001).

 Figure 6 

The disturbed metabolic pathways showed various metabolic changes in OSCC and OELP group. OSCC, oral squamous cell carcinoma; OELP, oral erosive lichen planus.

J Cancer Image

(View in new window)

 Figure 7 

Receiver operating characteristic curves of multiple logistic regression models with (A) decanoylcarnitine, (B) cysteine, (C) cholic acid and (D) biomarker panel. (E) The prediction accuracy of the diagnostic panel in validation set. OSCC, oral squamous cell carcinoma; OELP, oral erosive lichen planus; QC, quality control sample; AUC, area under curves; 95%CI, 95% confidence interval.

J Cancer Image

(View in new window)

 Table 3 

Metabolites different between OSCC and HC

IDNameFormulaMolecularRT (min)VIPPFCAUCFDR
1UridineC9H12N2O6244.062271.001.0281.02E-100.7180.8828.53E-10
2Chiro-inositolC6H12O6180.055531.011.2830.0271.0930.7420.049
3LactateC3H6O390.023221.011.1580.0261.0890.6320.047
4TaurineC2H7NO3S125.006391.011.9063.35E-130.6620.9064.07E-12
5GlutamateC5H9NO4147.045011.021.9071.96E-130.6500.9122.52E-12
6Uric acidC5H4N4O3168.020171.022.8588.53E-100.8820.6200.021
7AcetylcarnitineC9H17NO4204.122851.021.6621.23E-050.6820.7742.9E-03
8Citric acidC6H8O7192.019391.273.3041.37E-080.6160.8728.67E-07
9PhenylalanineC9H11NO2165.085912.102.2773.10E-050.7840.7850.003
10OctanoylcarnitineC15H29NO4287.216615.501.0031.52E-060.6710.8017.67E-06
11Glycocholic acidC26H43NO6465.302956.211.3339.60E-052.5880.7793.27 E-04
12DecanoylcarnitineC17H33NO4315.248236.251.0339.68E-060.6040.8123.54E-05
13Cholic acidC24H40O5408.281406.802.0600.0104.7510.7790.021
14Glycoursodeoxycholic acidC26H43NO5449.307656.932.3250.0041.7630.7260.009
15N-Succinyl-2,6-diaminoheptanedioateC11H18N2O7289.106177.201.7885.92E-120.4850.8996.01E-11
1616-HETEC20H32O3320.250988.031.0168.12E-130.3930.8107.50E-06
17LysoPC(16:0)C24H50NO7P495.339518.6812.4431.81E-200.5650.9701.19E-17
18LysoPC(18:1)C26H52NO7P521.355718.8711.0904.10E-190.6630.9506.38E-19
199,10-Epoxyoctadecenoic acidC18H32O3296.228329.223.1450.0020.6600.7060.008
20Arachidonic acidC20H32O2304.2331510.035.3016.77E-090.6350.8454.45E-08
21OleamideC18H35NO281.2787210.142.0713.50E-230.0920.5024.01E-05

OSCC, oral squamous cell carcinoma; HC, healthy control; FC, fold change; RT, retention time; VIP, variable importance in projection; FC, fold change; AUC, area under curve; FDR, false discovery rate, calculation based on Benjamini and Hochberg.

 Table 4 

Metabolites different between OELP and HC

IDNameFormulaMolecularRTVIPPFCAUCFDR
1LysoPC(18:1)C26H52NO7P521.347869.221.5231.87E-170.6680.8382.94E-16
2OleamideC18H35NO281.2787210.142.1085.63E-250.0690.7603.18E-23
3Chiro-inositolC6H12O6180.055563.251.5090.0060.9080.7830.013037
4ParaxanthineC7H8N4O2180.055563.251.5008.26E-0816.4160.7422.98E-05
5CysteineC3H7NO2S121.011384.471.9862.38E-370.4180.9491.11E-08
62,4-Dihydroxybenzoic acidC7H6O4154.018436.341.4960.0070.2430.7660.32
7Arachidonic acidC20H32O2304.2332210.033.6211.98E-041.1590.8085.65E-04
8UridineC9H12N2O6244.062381.001.0346.86E-110.7940.8366.89E-04
9Indolelactic acidC11H11NO3205.066274.731.0101.14E-051.0320.8123.77E-04
10Monobutyl phthalateC12H14O4222.081896.521.2162.58E-170.9000.7368.62E-08
11Malic acidC4H6O5134.013261.071.0734.76E-220.9480.8401.01E-20
12LactateC3H6O390.023231.011.0600.0051.0680.6820.012
13GlutamateC5H9NO4147.045071.032.5504.79E-280.8710.8932.26E-26
14HistidineC6H9N3O2155.061370.942.1051.81E-1313.3440.9161.20E-12
15Uric acidC5H4N4O3168.020161.025.0459.81E-061.0880.7823.26E-05
16Citric acidCH8O7192.019391.273.9349.63E-130.8630.9385.74E-12
17N-Succinyl-2,6-diaminoheptanedioateC11H18N2O7290.106207.191.6096.79E-110.4960.8433.51E-10
18TaurineC2H7NO3S125.006371.002.0111.55E-150.8920.9071.32E-14

OELP, oral erosive lichen planus; HC, healthy control; FC, fold change; RT, retention time; VIP, variable importance in projection; FC, fold change; AUC, area under curve; FDR, false discovery rate, calculation based on Benjamini and Hochberg.

According to the ROC curve, the Youden's index was calculated. The best cut-off value was 0.664. It was used to distinguish OSCC from OELP in the verification set. Samples above the cutoff value were diagnosed as patients with OSCC and below the cutoff value were diagnosed as patients with OLP. Only one case was wrong when the diagnostic panel was used to diagnose the validation set. The results presented that the panel achieved a diagnostic accuracy of 97.9% (Figure 7E).

Discovery of malignant transformation biomarkers

In order to gain additional insight into the metabolic features when OELP progressed to OSCC, plasma metabolomic analysis was performed to distinguish metabolites in OSCC, OELP and healthy controls. An additional group of 47 healthy volunteers were recruited as healthy controls (HC), who were well matched with aspects of age, gender, and body mass index (BMI) of the patients. 47 healthy controls with a mean age of 65 ± 9 (yrs.; mean ± SD). The demographic baseline characteristics were shown in Table 1. The metabolites of OSCC and HC were listed in Table 3 and OELP and HC were listed in Table 4. Both OSCC and OELP patient groups showed decreased FC in uridine, taurine, glutamate, citric acid and LysoPC (18:1). There was a significant difference (FDR < 0.05) between disease group and HC. However, the OSCC group showed a more pronounced decreased compared to OELP. Altered metabolites could possibly contribute to cancer progression. The metabolic pathway analysis was performed respectively and some common disturbed metabolic pathways were found. D-Glutamine and D-glutamate metabolism, primary bile acid biosynthesis, alanine, aspartate and glutamate metabolism, taurine and hypotaurine metabolism and arachidonic acid metabolism may be associated with the malignant change of OELP (Figure 8).

 Figure 8 

The disturbed metabolic pathways showed various metabolic changes in (A) OELP and HC group and (B) OSCC and HC group. OSCC, oral squamous cell carcinoma; OELP, oral erosive lichen planus; HC, healthy control.

J Cancer Image

(View in new window)

Discussion

Despite the continuous advancement of diagnostic and treatment modalities in the past 20 years, the 5-year survival rate of OSCC has not substantially improved [21]. Besides, there are no specific symptoms in the early stage of the disease, and it is difficult to make a diagnosis by pathological biopsy in the early stage. These could result in a series of consequences such as pain, bleeding, infection and even death. Thus, our research group has been seeking a diagnostic method which is accurate, convenient, and less harmful. We focused on blood samples because the blood contained metabolites that could be used for successful identification of diagnostic biomarkers. The imbalance of metabolites might be associated with pathological mechanism of OSCC.

One of the central results of our study was that some biomarkers decreased more in OSCC compared to OELP. Previous studies only showed that biomarkers in OSCC or OELP by comparing healthy controls, but few studies have focused on the metabolic alterations precede malignant transformation of OELP into OSCC. In our study, we found that the level of uridine, taurine, glutamate, citric acid and LysoPC (18:1) were lower in OELP than in OSCC. Uridine was generally thought safe and harmless, but recent studies found that uridine homeostatic disorder was carcinogenic. Uridine could be used to synthesize deoxyuridine triphosphate (dUTP). While dUTP was likely to cause errors in DNA replication. DNA damage may cause or greatly increase a person's susceptibility to cancer [22, 23]. It also made a person who had a higher likelihood of OSCC. In addition, uridine synthesis originated from glutamine and glutamine could be synthesized from glutamate and ammonia [24]. Glutamate was regarded as a potential diagnostic biomarker of OSCC in a previous study [25]. Amino acids are an important unit of energy source for basic metabolic pathways in human beings. The abnormalities in amino acid metabolism may be a unique sign of OSCC [26]. There was no doubt that proliferation of OSCC cells required more energy consumed than OELP. Therefore, compensatory mechanisms promoted excessive consumption of amino acids to maintain the normal physiological function in the body, which could further aggravate metabolic disorders and exacerbate the disease. Taurine could cause tumor cell apoptosis [27], showing anti-cancer effects [28]. In our study, taurine was greatly reduced in OSCC patients, which may be related to the proliferation of OSCC cells. The proliferation of OSCC cells may in turn inhibit the anti-tumor effect of taurine. From HC to OELP to OSCC, the decrease of uridine, glutamate and taurine may be used for the mass synthesis of damaged-DNA. Cancer cells display diverse metabolic reprogramming including reductive carboxylation in the citric acid cycle to use glutamine into intracellular lipid storage [29]. Hence, citric acid and LysoPC(18:1) may be considered as the biomarker for the transformation of OELP to OSCC because progression of OSCC could exhaust more energy. Collectively, our data indicated a microenvironment that was conducive to the rapid proliferation of cancer cells was gradually formed in OSCC patients.

Previous studies on metabolomics in OSCC included healthy individuals who were regarded as controls for OSCC group [30]. Although it was a common method when identifying diagnostic biomarkers, it didn't consider potential malignant lesions which could alter metabolic features. Currently, approximately 28 000 000 of OELP occur in the world [31] and 1.1% of OLP may transform into a malignant cancer [32]. Selecting OELP patients as a control group in OSCC studies could minimize confounding factors related to benign diseases.

Another important finding in this study was a biomarker panel of OSCC and OELP. As we know, any single diagnostic biomarker had limited diagnostic accuracy (Figure 7A-C). When we combined the four biomarkers together as a “diagnostic panel”, the panel offered superior diagnostic performance and showed significantly higher sensitivity (Figure 7D). Further understanding of these biomarkers may provide more insight into the pathogenesis of OSCC and serve new therapeutic interventions. A previous study showed that the level of acylcarnitine decreased in esophageal squamous cell carcinoma [33]. Carnitine could affect metabolic mechanisms in numerous ways. Acyl-coenzyme A synthetases could catalyze the thioesterification of coenzyme A (CoA) to acyl-CoA esters while Acyl-CoA esters could be converted to acyl-carnitine esters by carnitine acyltransferases [34]. Therefore, the decrease of decanoyl carnitine may result from the inhibition of the enzymes' activity and levels by OSCC cells. In addition, the incidence of lymph node metastasis and bone metastasis were associated with acyl carnitine [35]. It was possible that decanoyl carnitine could favor the identification and demarcation of OELP and OSCC. Amino acid was an important source of energy storage, and its abnormal metabolism may be an important characterization of cancer [26]. Cysteine had been regarded as a biomarker of oral cancer in both blood and saliva [36]. Also, cysteine was the precursor for the formation of glutathione. Glutathione could regulate the redox state and immune response in the system [37]. The decrease of cysteine may be related to the inhibition of immune system and a disorder in fatty acid oxidation metabolism by OSCC. Cholic acid could activate the TGR5 receptor which can induce OSCC cell proliferation [38]. This result was in line with the increase of cholic acid production in our study. The result suggested that a microenvironment which was suitable for carcinogenesis was gradually formed in the body of OSCC patients for promoting the tumor growth.

Taken together, we used UHPLC-Q-HRMS to analyze the differential metabolites in the plasma of patients with OSCC and OELP. We found that metabolic characterization and metabolic pathways in OSCC patients were significantly different from those in OELP. These changes in endogenous metabolites and abnormal metabolic pathways may be related to the pathogenesis of OSCC, which could be used for the diagnosis of OSCC. This study provided a basis for clinical molecular diagnosis and had important significance for clinical diagnosis and treatment of OSCC. In the future, more patients, including those with cancers, precancerous lesions and healthy controls, will be recruited to further verify the clinical applicability of the biomarkers described in this study.

Conclusion

In this study, a panel of metabolites that consist of decanoylcarnitine, cysteine and cholic acid was identified for the diagnosis of OSCC. The metabolites uridine, taurine, glutamate, citric acid and LysoPC(18:1) were found to be potential biomarkers indicating malignant transformation of OELP. Biomarkers based on plasma metabolomics could be very helpful for diagnosis of OSCC.

Acknowledgements

We would like to thank all the participants, as well as all volunteers who helped us collect samples.

Funding

This study was supported by the Key Scientific Research Project of Henan Institution of Higher Education (20A350018, 21A320025 and 21A320057), Henan Province Medical Science and Technology Research Plan (SBGJ202002078 and SBGJ202002116) and Kang Meng Medical Research Foundation (No. TB204022).

Competing Interests

The authors have declared that no competing interest exists.

References

1. Gatta G, Botta L, Sánchez MJ. et al. Prognoses and improvement for head and neck cancers diagnosed in Europe in early 2000s: The EUROCARE-5 population-based study. European journal of cancer (Oxford, England: 1990). 2015;51:2130-43

2. Torre LA, Bray F, Siegel RL. et al. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65:87-108

3. Giuliani M, Troiano G, Cordaro M. et al. Rate of malignant transformation of oral lichen planus: A systematic review. Oral Dis. 2019;25:693-709

4. Peiris-Pagès M, Martinez-Outschoorn UE, Pestell RG. et al. Cancer stem cell metabolism. Breast cancer research: BCR. 2016;18:55

5. Ji EH, Cui L, Yuan X. et al. Metabolomic analysis of human oral cancer cells with adenylate kinase 2 or phosphorylate glycerol kinase 1 inhibition. Journal of Cancer. 2017;8:298-304

6. Zhan B, Wen S, Lu J. et al. Identification and causes of metabonomic difference between orthotopic and subcutaneous xenograft of pancreatic cancer. Oncotarget. 2017;8:61264-81

7. Yang XH, Zhang XX, Jing Y. et al. Amino acids signatures of distance-related surgical margins of oral squamous cell carcinoma. EBioMedicine. 2019;48:81-91

8. Khanna R, Kumar K, Roy R. A case study of primary malignancy of buccal mucosa using proton HR-MAS NMR spectroscopy on tissue specimens. J Oral Biol Craniofac Res. 2018;8:68-72

9. Ohshima M, Sugahara K, Kasahara K. et al. Metabolomic analysis of the saliva of Japanese patients with oral squamous cell carcinoma. Oncol Rep. 2017;37:2727-34

10. Li XZ, Zhang SN, Yang XY. Serum-based metabolomics characterization of patients with reticular oral lichen planus. Arch Oral Biol. 2019;99:183-9

11. Ishikawa S, Sugimoto M, Edamatsu K. et al. Discrimination of oral squamous cell carcinoma from oral lichen planus by salivary metabolomics. Oral Dis. 2020;26:35-42

12. Song X, Yang X, Narayanan R. et al. Oral squamous cell carcinoma diagnosed from saliva metabolic profiling. Proc Natl Acad Sci U S A. 2020;117:16167-73

13. Yang XY, Zhang SN, Li XZ. et al. Analysis of human serum metabolome for potential biomarkers identification of erosive oral lichen planus. Clin Chim Acta. 2017;468:46-50

14. Nunes LA, Mussavira S, Bindhu OS. Clinical and diagnostic utility of saliva as a non-invasive diagnostic fluid: a systematic review. Biochemia medica. 2015;25:177-92

15. Wu G, Wu J, Pan X. et al. Racial disparities in alpha-fetoprotein testing and alpha-fetoprotein status associated with the diagnosis and outcome of hepatocellular carcinoma patients. Cancer medicine. 2019;8:6614-23

16. Huber PM, Afzal N, Arya M. et al. Prostate Specific Antigen Criteria to Diagnose Failure of Cancer Control following Focal Therapy of Nonmetastatic Prostate Cancer Using High Intensity Focused Ultrasound. The Journal of urology. 2020;203:734-42

17. Hayakawa H, Fukasawa M, Sato T. et al. Carcinoembryonic antigen level in the pancreatic juice is effective in malignancy diagnosis and prediction of future malignant transformation of intraductal papillary mucinous neoplasm of the pancreas. Journal of gastroenterology. 2019;54:1029-37

18. Wang XS, Sun Z, Liu LW. et al. Potential Metabolic Biomarkers for Early Detection of Oral Lichen Planus, a Precancerous Lesion. Frontiers in pharmacology. 2020;11:603899

19. van der Meij EH, van der Waal I. Lack of clinicopathologic correlation in the diagnosis of oral lichen planus based on the presently available diagnostic criteria and suggestions for modifications. J Oral Pathol Med. 2003;32:507-12

20. Long J, Liu L, Jia Q. et al. Integrated biomarker for type 2 diabetes mellitus and impaired fasting glucose based on metabolomics analysis using ultra-high performance liquid chromatography quadrupole-Orbitrap high-resolution accurate mass spectrometry. Rapid communications in mass spectrometry: RCM. 2020;34:e8779

21. Ferlay J, Soerjomataram I, Dikshit R. et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136:E359-86

22. Cao Z, Ma J, Chen X. et al. Uridine homeostatic disorder leads to DNA damage and tumorigenesis. Cancer letters. 2016;372:219-25

23. Yan Y, Han X, Qing Y. et al. Inhibition of uracil DNA glycosylase sensitizes cancer cells to 5-fluorodeoxyuridine through replication fork collapse-induced DNA damage. Oncotarget. 2016;7:59299-313

24. Zhang Y, Guo S, Xie C. et al. Uridine Metabolism and Its Role in Glucose, Lipid, and Amino Acid Homeostasis. BioMed research international. 2020;2020:7091718

25. Yang XH, Jing Y, Wang S. et al. Integrated Non-targeted and Targeted Metabolomics Uncovers Amino Acid Markers of Oral Squamous Cell Carcinoma. Front Oncol. 2020;10:426

26. Tsai CK, Yeh TS, Wu RC. et al. Metabolomic alterations and chromosomal instability status in gastric cancer. World journal of gastroenterology. 2018;24:3760-9

27. He F, Ma N, Midorikawa K. et al. Anti-Cancer Mechanisms of Taurine in Human Nasopharyngeal Carcinoma Cells. Advances in experimental medicine and biology. 2019;1155:533-41

28. Baliou S, Kyriakopoulos AM, Spandidos DA. et al. Role of taurine, its haloamines and its lncRNA TUG1 in both inflammation and cancer progression. On the road to therapeutics? (Review). Int J Oncol. 2020;57:631-64

29. Hua TNM, Kim MK, Vo VTA. et al. Inhibition of oncogenic Src induces FABP4-mediated lipolysis via PPARγ activation exerting cancer growth suppression. EBioMedicine. 2019;41:134-45

30. Shahid N, Iqbal A, Siddiqui AJ. et al. Plasma metabolite profiling and chemometric analyses of tobacco snuff dippers and patients with oral cancer: Relationship between metabolic signatures. Head Neck. 2019;41:291-300

31. Olson MA, Rogers RS 3rd, Bruce AJ. Oral lichen planus. Clinics in dermatology. 2016;34:495-504

32. Aghbari SMH, Abushouk AI, Attia A. et al. Malignant transformation of oral lichen planus and oral lichenoid lesions: A meta-analysis of 20095 patient data. Oral Oncol. 2017;68:92-102

33. Xu J, Chen Y, Zhang R. et al. Global and targeted metabolomics of esophageal squamous cell carcinoma discovers potential diagnostic and therapeutic biomarkers. Molecular & cellular proteomics: MCP. 2013;12:1306-18

34. Adeva-Andany MM, Calvo-Castro I, Fernández-Fernández C. et al. Significance of l-carnitine for human health. IUBMB life. 2017;69:578-94

35. Zoni E, Minoli M, Bovet C. et al. Preoperative plasma fatty acid metabolites inform risk of prostate cancer progression and may be used for personalized patient stratification. BMC cancer. 2019;19:1216

36. Guerra EN, Acevedo AC, Leite AF. et al. Diagnostic capability of salivary biomarkers in the assessment of head and neck cancer: A systematic review and meta-analysis. Oral Oncol. 2015;51:805-18

37. Mu Y, Zhou Y, Wang Y. et al. Serum Metabolomics Study of Nonsmoking Female Patients with Non-Small Cell Lung Cancer Using Gas Chromatography-Mass Spectrometry. J Proteome Res. 2019;18:2175-84

38. Amonyingcharoen S, Suriyo T, Thiantanawat A. et al. Taurolithocholic acid promotes intrahepatic cholangiocarcinoma cell growth via muscarinic acetylcholine receptor and EGFR/ERK1/2 signaling pathway. Int J Oncol. 2015;46:2317-26

Author contact

Corresponding address Corresponding authors: Hongyu Zhao, Department of Oral Emergency, The First Affiliated Hospital of Zhengzhou University· Stomatological Hospital of Henan Province, Zhengzhou, Henan, 450052, China. E-mail: zhongyu93com; Zhi Sun, Department of Pharmacy, the First Affiliated Hospital of Zhengzhou University, Zhengzhou, Henan, 450052, China. E-mail: sunzhi2013com.


Received 2021-2-24
Accepted 2021-11-2
Published 2022-1-1