ORIGINAL RESEARCH article

Front. Oncol., 15 October 2019

Sec. Gastrointestinal Cancers

Volume 9 - 2019 | https://doi.org/10.3389/fonc.2019.01069

Metabolomic Profiles for HBV Related Hepatocellular Carcinoma Including Alpha-Fetoproteins Positive and Negative Subtypes

  • 1. Beijing You'an Hospital, Capital Medical University, Beijing, China

  • 2. NHC Key Laboratory of Systems Biology of Pathogens, Institute of Pathogen Biology, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China

Abstract

Background: Hepatocellular carcinoma (HCC) is very common globally prevalent cancer. Due to its poor clinical prognosis, increasing the diagnostic rate of HCC is urgently needed. Herein, we validate discovered metabolomic biomarkers to distinguish Hepatitis B virus (HBV)-related HCC, including alpha-fetoprotein (AFP) negative (AFP–) and positive (AFP+) individuals.

Methods: We recruited 130 HCC subjects (independent case-control, randomized clinical cohorts) to our study. We separated the subjects randomly into two panels: (1) 58 individuals for the discovery panel; and (2) 72 individuals for the validation panel. For each panel, gender and age-matched hepatitis B group (HBG) and healthy group were included as controls. Plasma samples were collected for metabolic profiling by liquid chromatography—mass spectrometry—based metabolomics assays. We applied both non-targeted metabolomics analyses and targeted metabolomics analyses. Significantly changed metabolites (SCMs) were identified. The power of SCMs to discriminate HCC and HBG or healthy group was determined by receiver operating characteristic curve (ROC) analysis.

Results: Ten SCMs were selected form the discovery panel, and further verified in the validation panel. ROC analyses indicated that 1 SCMs (LysoPC (24:0)) could discriminate HCC from HBG (AUC = 0.765). Further, 8 SCMs including (LysoPC (17:0), LysoPC (20:4(8Z,11Z,14Z,17Z)), LysoPC (22:0), LysoPC (24:0), PE (P-16:0/22:4(7Z,10Z,13Z,16Z)), SM (d18:1/22:1(13Z)), Creatinine, and L-Isoleucine) displayed a heightened ability to discriminate between HCC and healthy controls (AUC were more than 0.800). Most of these SCMs were important in lipid metabolism.

Conclusions: LysoPC (24:0) could distinguished HCC from HBG, and 8 SCMs distinguished HCC from healthy controls. LysoPC and other metabolites have the potential to serve as non-invasive biomarkers for HBV related AFP– and AFP+ HCC.

Introduction

Hepatocellular carcinoma (HCC) is the third leading cause of cancer mortality in the world. About 278.07/100000 new cases, and 167.89/100000 deaths were reported each year in China (1). Hepatitis B virus (HBV) is the most common cause of HCC, accounting for more than 50% of total mortality (24). HCC diagnosed at a late, advanced stage of progression, and usually with ineffective therapies being available, which contributes to the high mortality rate in HCC (5). In addition, a lack of reliable biomarkers to diagnose HCC remains a significant challenge.

Alpha-Fetoprotein (AFP), which was first identified in the 1970s, has often been used as a blood biomarker for HCC (6). However, the sensitivity of AFP is limited to 65% for clinical HCC diagnoses and <40% for preclinical predictions (79). The sensitivity and specificity of AFP are not fully applicable to clinical use at this time. For AFP negative (AFP–) HCC, liver biopsy is the current gold standard in clinical diagnosis; however, as an invasive diagnostic approach, it suffers from many limitations, and a non-invasive diagnostic tool is required. Therefore, it is necessary to develop reliable biomarkers to further strengthen the goal of increasing the diagnosis rate in HCC.

Metabolomics represents a rapidly emerging approach to study small molecules (i.e., molecular weight <1,800 Da) that define the metabolic status of a biological system. Such small molecules play important roles in biological systems. Thus, metabolomic profiling has been used to explore the discovery of useful early biomarkers in many human diseases (10). Metabolomics has the potential to precisely measure the metabolite complement in living systems, and can capture the global dynamic responses to changes that from both endogenous and exogenous factors. Metabolomics clearly provides us with many potential applications and advantages for research into complex systems, and a useful platform to detect and study dysregulated metabolic pathways in the field of HCC genesis (11).

It has been reported that a serological panel of metabolite biomarkers exhibits good diagnostic performance in the early detection of HCC as compared detecting the levels of AFP (12). Through a metabolomics approach, Wang et al. described the identification of 13 potential biomarkers and corresponding pathways that were significantly apparent in patients with HCC (10). Of the candidate biomarkers, glycochenodeoxycholic acid was proposed as a potential indicator for HCC diagnosis (13). Although, there are many studies striving to identify new biomarkers for HCC by metabolomics approaches (1013), the majority of these studies have not included AFP–HCC as a study group. In this context, we have recruited gender- and age-matched AFP negative (AFP–) HCC and compared this group with AFP positive (AFP+) HCC, with the over-arching aim of identifying metabolomic biomarkers that are characteristic of HCC.

Moreover, to exclude the effects of HCC therapy, we chose HCC that never received HCC therapy when recruited. Our study quantified metabolomic profiles in plasma samples amongst the different groups, and to discover diagnostic biomarkers for HCC, a special metabolic pathway is also included.

Materials and Methods

Clinical Design and Sample Collection

One hundred and thirty individuals participated in this study at Beijing You'an Hospital of the Capital Medical University from 2014 to 2017, in a case-control, randomized clinical cohort study. Diagnosis of HCC and chronic hepatitis B were based on the EASL Clinical Practice Guidelines (14) and the AASLD practice guidelines (15). Entry criteria included: (1) being diagnosed as HCC by liver biopsy or CT/MRI; (2) individuals with positive HBV surface antigen (HBsAg). Exclusion criteria included: (1) individuals undergoing treatment; for example, interventional surgery, radiofrequency ablation, liver transplantation, TACE, liver cancer resection, and any other treatments for HCC; (2) other hepatitis virus infections; (3) combined presentation of metabolic related disease, including diabetes or cardiovascular diseases; (4) combined presentation of autoimmune-related disease; (5) individuals with a recurrent tumor.

Next, the level of AFP (Roche, Cobas e601, Germany) in the AFP–HCC group was normal (AFP ≤ 7 ng/ml), and the level of AFP in the AFP+HCC group was AFP > 20 ng/ml. Based on the criteria, HCC individuals were separated into a discovery panel and a validation panel according to their age (40–50 years for the discovery panel and 50–70 years for the validation panel). The gender- and age-matched hepatitis B group (HBG) and healthy groups were recruited as controls for this study and were randomly separated into the discovery and validation panels.

All individuals were subjected to the written informed consent requirement, and the local ethics committee of Beijing You'an hospital, Capital Medical University approved this study. The plasma samples were collected immediately after being recruited, following which, the levels of HBsAg, hepatitis B virus surface antibody, hepatitis B virus e antigen, hepatitis B virus e antibody, hepatitis B virus core antibody, HBV DNA, liver function and blood routine analyses were measured. All plasma samples were stored at −80°C until use.

Plasma Sample Preparation

Plasma samples were thawed on ice at 4°C for 30–60 min, then vortexed for a few seconds and centrifuged at 12,000 rpm for 10 min at 4°C. Aliquots of 100 μl plasma were placed into labeled microcentrifuge tubes for analysis. Pooled quality control samples (QCs) consisted of small aliquots of each biological sample being pooled, which were mixed together to monitor the stability and reproducibility of the results.

Screening Metabolite Profiles of Different Groups

The metabolites in the discovery panel including 58 individuals (Figure 1), were screened using the Thermo Scientific™ Dionex™ UltiMate™ 3000 Rapid Separation LC (RSLC) system. A reversed phase C18 column and hydrophilic interaction liquid chromatography column (HILIC) were used for UHPLC separations (positive and negative mode ions combined).

Figure 1

RP Separation for Lipids

For the C18 separation, mobile phase A was acetonitrile/water (60/40) and mobile phase B was isopropanol/acetonitrile (90/10). Both mobile phase A and B contained 0.1% formic acid and 10 mmol/L ammonium acetate. Gradient conditions for reversed phase C18 separation is shown in Table S1. The column was HSS T3 (2.1 × 100 mm, 1.8 μm, water) that was operated at 45°C. The flow rate was 300 μL/min and the injection volume was 1 μL.

HILIC Separation

For the HILIC separation, mobile phase A was acetonitrile and mobile phase B was water; both mobile phase A and B contained 0.1% formic acid and 10 mmol/L ammonium acetate. The column was a BEH amide column (2.1 × 100 mm, 1.7 μm, water) that was operated at 40°C. The flow rate was 300 μL/min and the injection volume was 1 μL. The gradient conditions for HILIC separation is shown in Table S1.

Data Acquisition

After the chromatographic separations, using a Thermo Scientific™ Q Exactive™ hybrid quadrupole Orbitrap mass spectrometer equipped with a HESI-II probe, we could perform data acquisition. The positive and negative HESI-II spray voltages were 3.7 and 3.5 kV, respectively. The heated capillary temperature was set to 320°C. The sheath gas pressure was 30 psi. The auxiliary gas setting was 10 psi, and the heated vaporizer temperature was 300°C. Both the sheath gas and the auxiliary gas were nitrogen. The collision gas was also nitrogen at a pressure of 1.5 mTorr. The parameters of the full mass scan were as follows: a resolution of 70,000, an auto gain control target that was under 1 × 106, a maximum isolation time of 50 ms, and an m/z range of 50–1,500. The LC-MS system was controlled using Xcalibur 2.2 SP1.48 software (Thermo Fisher Scientific), and data were collected and processed with the same software. All data that was obtained using both positive and negative ion modes were processed using the Progenesis QI data analysis software (Non-linear Dynamics, Newcastle, UK). The ranges of automatic peak picking for the C18 and HILIC assays were between 1 and 19 min and between 1 and 12 min, respectively. Next, the adduct ions of each “feature” (m/z, tR) were deconvoluted, and these features were identified in the human metabolome database (HMDB) and lipid maps.

Metabolomic Data Analysis

The raw data were screened by correcting individual bias using QC and blank data sets. The screened data was subjected to Principal Component Analysis (PCA), Orthogonal signal correction Partial Least Square Discrimination Analysis (OPLS-DA), Variable Importance in Projection (VIP), and coefficients vs. VIP spots with the SIMCA 14.1 software program (Umetrics AB, Umea, Sweden).

Verification of the Metabolite Profiles in the HCC Group

To verify the metabolites in the HCC group, we used an independent cohort with 72 individuals as a validation panel (Figure 1). All samples were subjected to UHPLC separation by the Thermo Scientific™ Dionex™ UltiMate™ 3000 Rapid Separation LC (RSLC) system. The gradient conditions for the C18 column and the HILIC column was the same as the above.

Metabolite Enrichment

Metabo-Analyst version 4.0 was used for pathway enrichment analysis. The software was obtained from http://www.metaboanalyst.ca/faces/ModuleView.xhtml.

Statistical Analysis

GraphPad Prism software (version 6.0, San Diego, California, USA) was used for statistical analysis. Continuous variables presented as mean ± standard deviation (SD). We used the Mann-Whitney U-test and Fisher's exact test to compare the statistical significance between different groups, respectively. The level of statistical significance was set at p < 0.05 two-sided for all tests. The area under the receiver-operating characteristic (ROC) curve (AUC) was calculated to evaluate the classification performance.

Results

Clinical Characteristics of the Subjects

Overall, 30 individuals with HBV related AFP+HCC (median 53.8 years, 25 males and 5 females), and 40 individuals with HBV related AFP–HCC (median 53.4 years, 34 males and 6 females) were recruited. Thirty HBG and 30 healthy individuals were recruited as controls. Among them, 58 individuals in the discovery panel, and 72 individuals in the validation panel were recruited (Figure 1, Table 1).

Table 1

GroupHCCControlp*
AFP+HCCAFP–HCCHBGHealthy
58 SUBJECTS
Number15131614
Male/female13/211/215/112/20.701
Age46.6 ± 3.22543.85 ± 5.36742.31 ± 3.97944.79 ± 2.6940.090
Ethnicity0.183
  Han14121413
  Hui1020
  Mongols0101
BCLC stage
  002
  A15
  B54
  C82
  D10
Tumor size
  ≥3 cm94
   <3 cm27
  Missing42
AFP24,042 ± 33,9523.162 ± 1.2373.122 ± 1.1830.002
WBC6.774 ± 2.5684.673 ± 1.5765.154 ± 1.9210.400
HGB142.5 ± 14.42139.5 ± 19.7155.2 ± 16.510.015
PLT150.7 ± 76.18165.5 ± 73.62163.3 ± 80.260.861
ALT45.49 ± 21.5140.25 ± 31.0141.09 ± 27.2224.49 ± 6.5770.087
AST53.76 ± 23.0937.46 ± 13.0230.13 ± 11.6630.29 ± 5.2070.001
TBIL21.66 ± 12.218.57 ± 7.96413.84 ± 4.4550.041
DBIL7.233 ± 5.6884.846 ± 2.82.756 ± 1.1830.002
TP68.09 ± 6.45265.58 ± 4.21271.44 ± 5.8310.014
ALB38.27 ± 4.38840.41 ± 3.76345.04 ± 4.9030.000
CREA64.01 ± 15.3363.52 ± 7.64965.72 ± 200.252
CHOL4.698 ± 1.6423.813 ± 0.72594.319 ± 1.0930.592
r-GT215.7 ± 233.650.95 ± 35.9745.65 ± 28.940.057
ALP127.4 ± 53.1687.18 ± 30.8373.29 ± 19.560.027
PT12.27 ± 1.50712.19 ± 1.24911.03 ± 1.2290.003
PT%89.47 ± 15.8489.92 ± 14.44106 ± 18.120.003
PTINR1.094 ± 0.12971.086 ± 0.10670.9867 ± 0.10560.001
APTT32.51 ± 5.31333.07 ± 2.95233.18 ± 3.340.804
HBV DNA >20,000 IU/mL3/15 (20.00%)3/13 (23.08%)4/16 (25.00%)0.199
72 SUBJECTS
Number15271416
Male/female12/323/410/413/30.561
Age59.71 ± 6.9657.23 ± 5.9055.86 ± 4.9956.4 ± 4.6870.247
Ethnicity0.422
  Han14261315
  Man0110
  Mongols1001
BCLC stage
  000
  A26
  B411
  C35
  D65
Tumor size
  ≥3 cm713
   <3 cm514
  Missing30
AFP2,185 ± 3,3403.64 ± 1.7613.202 ± 2.6830.003
WBC5.639 ± 1.9786.089 ± 2.5564.751 ± 1.890.087
HGB138.2 ± 15.13140.3 ± 18.49148.4 ± 9.5570.064
PLT140.6 ± 69.4151.2 ± 63.34139.2 ± 50.40.892
ALT38.56 ± 20.238.61 ± 23.6347.62 ± 65.6319.65 ± 16.070.439
AST42.2 ± 16.2934.3 ± 11.4736.55 ± 34.3929.91 ± 10.020.052
TBIL19.2 ± 9.419.91 ± 11.6615.92 ± 6.1940.305
DBIL5.612 ± 2.7346.277 ± 8.0843.257 ± 1.460.023
TP64.24 ± 6.53166.91 ± 5.7869.12 ± 11.380.046
ALB36.99 ± 4.79140.48 ± 3.80442.91 ± 4.8380.009
CREA67.69 ± 14.2266.67 ± 11.8268.98 ± 18.63> 0.999
CHOL3.929 ± 0.83633.857 ± 0.77684.384 ± 0.99110.070
r-GT99.65 ± 100.668.1 ± 48.0342.22 ± 27.40.028
ALP88.01 ± 43.579.35 ± 28.1584.6 ± 23.320.584
PT12.33 ± 1.3811.71 ± 1.1911.52 ± 0.890.389
PT%88.41 ± 16.6595.53 ± 15.796.96 ± 12.160.475
PTINR1.098 ± 0.11981.046 ± 0.10521.02 ± 0.080.261
APTT32.28 ± 3.32232.77 ± 5.48733.33 ± 3.280.254
HBV DNA >20,000 IU/mL2/15 (13.33%)3/27 (11.11%)1/14 (7.14%)0.199

Clinical characteristics of subjects belong to HCC and Control group.

*

Compared between HCC and control groups. Bold and italic represent significant differences between the HCC and control groups. HCC, Hepatocellular carcinoma; BCLC, Barcelona Clinic Liver Cancer stage; AFP, alpha-fetoprotein; HBG, Hepatitis B group; WBC, White blood cell; HGB, Hemoglobin; PLT, Platelets; ALT, Alanine aminotransferase; AST, Aspartate aminotransferase; TBIL, Total bilirubin; DBIL, Direct bilirubin; TP, total protein; ALB, albumin; CREA, Creatinine; CHOL, Cholesterol; r-GT, r-glutamyl transpeptadase; ALP, Alkaline phosphatase; PT, Prothrombin time; PT%, Prothrombin time percent; PTINR, Prothrombin time international normalized ration; APTT, Activated partial thromboplatin time; HBV DNA, Hepatitis B virus deoxyribonucleic acid.

Fifty eight individuals including AFP+HCC (n = 15), AFP–HCC (n = 13), HBG (n = 16), and healthy controls (n = 14) were analyzed for metabolic profiles (Table 1). Clinical characteristics were compared between the AFP+HCC, AFP–HCC, HBG, and healthy control groups (Table 1), and the mean age were 46.6 ± 3.225, 43.85 ± 5.367, 42.31 ± 3.979, and 44.79 ± 2.694, respectively. The level of AFP was <7 ng/ml in the AFP–HCC and HBG groups.

An independent cohort that included 72 individuals was recruited for validation (Table 1), included AFP+HCC (n = 15), AFP–HCC (n = 27), HBG (n = 14), and healthy controls (n = 16), and the mean age were 59.71 ± 6.96, 57.23 ± 5.90, 55.86 ± 4.99, and 56.4 ± 4.687, respectively. The level of AFP was <7 ng/ml in the AFP–HCC and HBG groups.

There was no significant difference found in terms of the age and gender among the HCC groups, and the control groups (including HBG and healthy groups) in those two panels (p > 0.05). The clinical features that included assay of direct bilirubin (DBIL), total protein (TP), and albumin (ALB) were significantly different when comparing the HCC and control groups (p < 0.05), in both the discovery panel and the validation panel (Table 1).

The clinical characteristics of the discovery panel and the validation panel among the AFP+HCC, AFP–HCC, HBG, and healthy groups were all compared. With the exception of age, a significantly different outcome was identified for the discovery panel and validation panel (p < 0.05), and most of the others that were studied were statistically non-significant (Table 2).

Table 2

GroupAFP+HCCAFPHCCHBGHealthy
Discovery panelValidation panelp*Discovery panelValidation panelp*Discovery panelValidation panelp*Discovery panelValidation panelp*
Number1515132716141416
Male/female13/212/30.85211/223/40.94815/110/40.31312/213/30.854
Age46.6 ± 3.22559.71 ± 6.960.00043.85 ± 5.36757.23 ± 5.900.00042.31 ± 3.97955.86 ± 4.990.00044.79 ± 2.69456.4 ± 4.6870.000
Ethnicity1.00.7890.4950.391
  Han1414122614131315
  Hui10002000
  Man00010100
  Mongols01100011
BCLC stage1.0000.147
  00020
  A1256
  B54411
  C8325
  D1605
Tumor size1.0000.301
  ≥3 cm97413
   <3 cm25714
  Missing4320
AFP24,042 ± 33,9522,185 ± 3,3400.0643.162 ± 1.2373.64 ± 1.7610.4253.122 ± 1.1833.202 ± 2.6830.120
WBC6.774 ± 2.5685.639 ± 1.9780.2024.673 ± 1.5766.089 ± 2.5560.0305.154 ± 1.9214.751 ± 1.890.580
HGB142.5 ± 14.42138.2 ± 15.130.502139.5 ± 19.7140.3 ± 18.490.927155.2 ± 16.51148.4 ± 9.5570.154
PLT150.7 ± 76.18140.6 ± 69.40.794165.5 ± 73.62151.2 ± 63.340.488163.3 ± 80.26139.2 ± 50.40.473
ALT45.49 ± 21.5138.56 ± 20.20.31340.25 ± 31.0138.61 ± 23.630.59241.09 ± 27.2247.62 ± 65.630.58024.49 ± 6.57719.65 ± 16.070.003
AST53.76 ± 23.0942.2 ± 16.290.23037.46 ± 13.0234.3 ± 11.470.44030.13 ± 11.6636.55 ± 34.390.98430.29 ± 5.20729.91 ± 10.020.142
TBIL21.66 ± 12.219.2 ± 9.40.73718.57 ± 7.96419.91 ± 11.660.96913.84 ± 4.45515.92 ± 6.1940.448
DBIL7.233 ± 5.6885.612 ± 2.7340.6284.846 ± 2.86.277 ± 8.0840.9692.756 ± 1.1833.257 ± 1.460.313
TP68.09 ± 6.45264.24 ± 6.5310.24665.58 ± 4.21266.91 ± 5.780.47271.44 ± 5.83169.12 ± 11.380.637
ALB38.27 ± 4.38836.99 ± 4.7910.62840.41 ± 3.76340.48 ± 3.8040.92745.04 ± 4.90342.91 ± 4.8380.142
CREA64.01 ± 15.3367.69 ± 14.220.41163.52 ± 7.64966.67 ± 11.820.41065.72 ± 2068.98 ± 18.630.984
CHOL4.698 ± 1.6423.929 ± 0.83630.3503.813 ± 0.72593.857 ± 0.77680.6484.319 ± 1.0934.384 ± 0.99110.780
r-GT215.7 ± 233.699.65 ± 100.60.31350.95 ± 35.9768.1 ± 48.030.18445.65 ± 28.9442.22 ± 27.40.561
ALP127.4 ± 53.1688.01 ± 43.50.14287.18 ± 30.8379.35 ± 28.150.86573.29 ± 19.5684.6 ± 23.320.354
PT12.27 ± 1.50712.33 ± 1.380.76612.19 ± 1.24911.71 ± 1.190.30011.03 ± 1.22911.52 ± 0.890.384
PT%89.47 ± 15.8488.41 ± 16.650.76689.92 ± 14.4495.53 ± 15.70.366106 ± 18.1296.96 ± 12.160.650
PTINR1.094 ± 0.12971.098 ± 0.11980.7371.086 ± 0.10671.046 ± 0.10520.3130.9867 ± 0.10561.02 ± 0.080.384
APTT32.51 ± 5.31332.28 ± 3.3220.57633.07 ± 2.95232.77 ± 5.4870.45633.18 ± 3.3433.33 ± 3.280.837
HBV DNA >20,000 IU/mL3/15 (20.00%)2/15 (13.33%)0.7103/13 (23.08%)3/27 (11.11%)0.5054/16 (25.00%)1/14 (7.14%)0.423

Clinical characteristics of discovery panel and validation panel.

*

Compared between discovery panel and validation panel. Bold and italic represent significant difference between discovery panel and validation panel. HCC, Hepatocellular carcinoma; BCLC, Barcelona Clinic Liver Cancer stage; AFP, alpha-fetoprotein; HBG, Hepatitis B group; WBC, White blood cell; HGB, Hemoglobin; PLT, Platelets; ALT, Alanine aminotransferase; AST, Aspartate aminotransferase; TBIL, Total bilirubin; DBIL, Direct bilirubin; TP, total protein; ALB, albumin; CREA, Creatinine; CHOL, Cholesterol; r-GT, r-glutamyl transpeptadase; ALP, Alkaline phosphatase; PT, Prothrombin time; PT%, Prothrombin time percent; PTINR, Prothrombin time international normalized ration; APTT, Activated partial thromboplatin time; HBV DNA, Hepatitis B virus deoxyribonucleic acid.

PCA and OPLS-DA

Liquid chromatography was applied to separate lipid and polarized components from the plasma samples. The metabolic profiles of the plasma samples in the AFP+HCC, AFP–HCC, HBG, and healthy groups were displayed by PCA (Figures 2A,B). PCA analysis suggested distinct clusters of samples from different groups in both the C18 column (R2X = 0.481, Q2 = 0.392; Figure 2A) and the HILIC columns (R2X = 0.468, Q2 = 0.272; Figure 2B). OPLS-DA was further used to analyze the metabolic differences among the AFP+HCC with healthy groups (Figure 2C for C18 column, Figure 2D for HILIC columns), AFP–HCC with healthy groups (Figure 2E for C18 column, Figure 2F for HILIC columns), HBG with healthy groups (Figure 2G for C18 column, Figure 2H for HILIC columns). All OPLS-DA models between different groups indicated adequate classification.

Figure 2

Metabolic Profiles in Different Groups Compared With the Healthy Group

Significantly changed metabolites (SCM; defined by based VIP > 1 and p < 0.05) were selected for subsequent chemical structure identification. The metabolic profiles of the AFP+HCC, AFP–HCC, HBG were substantially different from that of the healthy control group, both in the C18 column (Figure 3A and Table S2) and the HILIC column (Figure 3B and Table S3). In addition, a substantial difference was observed between the AFP+HCC, AFP–HCC, and HBG groups both in the C18 column (Table S2) and the HILIC column (Table S3).

Figure 3

Metabolic Profiles in the HCC Group as Compared With HBG

We then selected common metabolites in the C18 column (Figure 3C) and the HILIC column (Figure 3D) from the AFP+HCC and the AFP–HCC groups, respectively, and as compared with HBG. There were 13 SCMs for the C18 column (Table S4) and 7 SCMs for the HILIC column (Table S5) identified for both the AFP+HCC and the AFP–HCC groups. In addition, there were 13 SCMs for the C18 column and 6 out of 7 SCMs for the HILIC column (with the exception of creatine) with significant difference when compared HCC with healthy control groups (p < 0.05; Figure 4). Otherwise, comparing HCC group with HBG group, we found 7 out of 13 SCMs for C18 column (LysoPC (24:0), LysoPC (17:0), PC (20:4(8Z,11Z,14Z,17Z)/20:4(5Z,8Z,11Z,14Z)),PC (20:4(5Z,8Z,11Z,14Z)/18:2(9Z,12Z)),PC (18:0/22:5(7Z,10Z,13Z,16Z,19Z)),PC (20:4(5Z,8Z,11Z,14Z)/15:0), PC (18:0/22:6(4Z,7Z,10Z,13Z,16Z,19Z))), and 2 out of 7 SCMs for HILIC (L-asparagine and ornithine) with significant difference (p < 0.05; Figure 4). Moreover, there were no significant difference between the AFP+HCC and the AFP–HCC groups for these SCMs (13 SCMs for the C18 column and 7 SCMs for the HILIC column, data not shown).

Figure 4

Validation of the Potential Metabolic Biomarkers

We validated these metabolites (13 SCMs for the C18 column and 7 SCMs for the HILIC column) in a validation panel. Targeted metabolomic analysis was performed and the concentration of SCMs was measured. Compared with the healthy control group, 7 SCMs for the C18 column and 5 SCMs for the HILIC column were identified (p < 0.05; Figure 5). In addition, when compared with HBG group, there were 5 SCMs for the C18 column and none SCMs for the HILIC column were identified (p < 0.05; Figure 5). Moreover, there were also no significant difference when comparing the AFP+HCC group and the AFP–HCC group for these SCMs (13 SCMs for the C18 column and 7 SCMs for the HILIC column, data not shown).

Figure 5

Comparing the tendency of these metabolites (13 SCMs for the C18 and 7 SCMs for the HILIC columns) in the training and validation panels, we found that 10 SCMs had the same metabolic tendency in both panels, and all of them were significantly decreased in the HCC group as compared with the healthy controls. Then, compared with the HBG group, only LysoPC (24:0) was significantly decreased in the HCC groups for both panels.

ROC analyses, which calculate the area under the curve (AUC), were applied for 10 SCMs. Six SCMs for the C18 column, including LysoPC (17:0) (Figure 6A), LysoPC (20:4(8Z,11Z,14Z,17Z)) (Figure 6B), LysoPC (22:0) (Figure 6C), LysoPC (24:0) (Figure 6D), PE (P-16:0/22:4(7Z,10Z,13Z,16Z)) (Figure 6E), SM (d18:1/22:1(13Z)) (Figure 6F) had a diagnostic value in distinguishing the HCC group from the healthy group. The AUCs were 0.845, 0.865,0.863, 0.829, 0.821, and 0.851, respectively. Two SCMs for the HILIC column including creatinine (Figure 6G), and L-isoleucine (Figure 6H) had diagnostic value in distinguishing the HCC group from the healthy group, with AUC values of 0.934 and 0.824, respectively. Finally, to detect the diagnostic value of LysoPC (24:0) (Figure 6I) between the HCC and HBG groups, the ROC analyses showed that the AUC was 0.765. And none of the SCMs for the HILIC column had diagnostic value in distinguishing HCC from HBG.

Figure 6

The Metabolic Pathway of HCC

MetaboAnalyst version 4.0 was used to identity the related pathways of those SCMs. Three metabolic pathways (with larger circles in the Figure 7A) were significantly enriched, including Glycerophospholipid metabolism (circle 19), Arginine and proline metabolism (circle 17), and Glycine, serine, and threonine metabolism (circle 14). Six of eight SCMs were correlated with lipid metabolism, including LysoPC (22:0), LysoPC (24:0), LysoPC (17:0), LysoPC (20:4(8Z,11Z,14Z,17Z)), PE (P-16:0/22:4(7Z,10Z,13Z,16Z)), and SM (d18:1/22:1(13 Z)). Additionally, all of them were significantly decreased in the HCC group (Figure 5).

Figure 7

To better understand the mechanism of these SCMs, the reported signaling pathways were involved in cell growth, cell proliferation, and cell cycle in HCC and the potential metabolic signaling pathways that were identified in our work were schematically depicted in Figure 7B. Importantly, those identified signaling pathways indicated metabolic dysregulation during the process of HCC tumorigenesis.

Discussion

In order to improve the sensitivity and specificity of diagnostic relevant metabolomic biomarkers for HCC patients, we collected plasma samples from HBV infected HCC. This included data sets for gender- and age-matched AFP–HCC and AFP+HCC subjects, as well as HCC groups, and HBG and healthy individuals as control groups, so that we could perform a comprehensive metabolomic profiling analysis. In results we found 20 common SCMs in both the AFP–HCC and the AFP+HCC groups. Considering that metabolomics analyses might be affected by age, we divided the discovery panel and validation panel to individuals that were aged 40–50 and 50–70 years, respectively. Then, we further validated these 20 SCMs in 72 subjects, and found that there were 10 SCMs displaying the same tendency in both the discovery panel and the validation panel.

Next, we analyzed the diagnostic value of the markers among HCC, HBG, and healthy groups. The ROC curve analyses indicated that 8 SCMs had a high discriminatory ability when comparing HCC and healthy controls, indicating that they might serve as potential biomarkers. Compared to HBG, the abundance of LysoPC (24:0) in the HCC group was significantly decreased in both the discovery and validation panel. The AUC value for LysoPC (24:0) reached a value of 0.765, which suggested that it can discriminate HCC and HBG individuals (Figure 6I).

The functional cluster analysis showed that 6 out of 8 SCMs, including LysoPC (17:0) (Figure 6A), LysoPC (20:4(8Z,11Z,14Z,17Z)) (Figure 6B), LysoPC (22:0) (Figure 6C), LysoPC (24:0) (Figure 6D), PE (P-16:0/22:4(7Z,10Z,13Z,16Z)) (Figure 6E), and SM (d18:1/22:1(13Z)) (Figure 6F), whose molecular formulas are C25H52NO7P, C28H50NO7P, C30H62NO7P, C32H66NO7P, C43H78NO7P, and C45H89N2O6P, were correlated with lipid metabolism. Among them, LysoPC (22:0) belongs to members of the lysophosphatidylcholines (LPCs), which are formed by hydrolysis of phosphatidylcholines (PCs) by the enzyme phospholipase A2 (16). The other SCMs were PE (P-16:0/22:4(7Z, 10Z, 13Z, 16Z)) and SM (d18:1/22:1(13Z)), which are essential components of membrane structure, signal transduction and lipoprotein metabolism.

L-isoleucine has a molecular formula of C6H13NO2 and is an essential amino acid that is involved in protein synthesis in the human body, and is an aromatic amino acid. The metabolites included L-isoleucine, and provided a diagnostic model that could discriminate HCC and normal control very well (17).

The pathway analysis showed that the different molecules played important roles in activated tumor cell function and tumorigenesis, and included protein and lipid metabolism, and especially sphingolipid metabolism (18). Our ROC analysis showed that lipid metabolism was significantly different among the HCC, HBG, and healthy control groups. Prior studies considered that the development of HCC can cause a hypoxic microenvironment; thus, tumor cells need lipid metabolism to change the metabolic mode and adapt to the change (19, 20). Our observations further support this notion, and increased cellular lipid metabolism and increased consumption resulted in decreased expression levels of lipid metabolite in the plasma. According to the metabolic character of the lipid, activation of lipid metabolism in cells, and especially that of LysoPE and LysoPC, can promote synthesis of human phosphatidyl-ethanolamine-binding protein 4 (hPEBP4), which is an anti-apoptotic protein, and can associate with HBV to activate the PI3K/Akt/mTOR signaling pathway to promote tumorigenesis, development, and metastasis (2123).

Nutritional and functional decline (cachexia) is a common phenomenon seen in cancer patients, and even the loss of skeletal muscle is tightly connected to poor performance status, and poor response to treatment and a poor prognostic outcome for those patients (24, 25). In our results, creatinine and creatine in HCC were down-regulated as compared with the healthy control group. In addition, it has been reported that it might imply the loss of muscle in HCC patients, which might then contribute to the observed high morbidity and mortality seen in HCC (24).

As the main pathogenic agent of HCC genesis, many studies have found that AFP can mediate HBx protein activation, and promote the development of HCC (26), with an accompanied suppression of host protective immunity (27). Recently, research studies have found that AFP high and AFP low subjects display differential proteomic profiles, and that the poor prognosis seen in HCC is associated with cholesterol homeostasis (28). According to the clinical data and pathway analysis, our results also showed that HBV DNA levels were higher in the AFP+HCC group than that found in the AFP–HCC group. In addition, both HBV and AFP can biologically collaborate to promote HCC development. At the same time, AFP possesses a variety of biological functions. The clinical data showed that the level of DBIL in the AFP+HCC group was higher than that found in the AFP–HCC group (Table 1, p < 0.05). According to the above analysis, different AFP expression levels might have a differential metabolic profile.

In this study, we found some biomarkers might have utility in diagnosing HBV-associated HCC. Several studies have also proposed biomarkers for each HCC stage, from prevention to progression (29). And in our results, we also found that some metabolomic biomarkers were correlated with tumor stage or tumor size (data not shown). Thus, future studies will strive to investigate those biomarkers in larger subject cohorts, including different stages of HCC. We will also seek to explore the different metabolic biomarkers that might distinguish AFP–HCC and AFP+HCC.

The employment of more than one biomarker could provide a dynamic and powerful approach to advance our understanding of the spectrum of HCC, with diverse applications in clinical epidemiology, laboratory screening, and physician confirmed diagnosis and prognosis. Our results indicated that fatty acid and lipid metabolism might promote HCC genesis, including the interplay of glycerophospholipid metabolism and some select amino acid metabolic activity; however, these ideas require formal confirmation. On the other hand, the limitation of this study, is that there were several factors, which included environmental or behavioral factors (e.g., diet) that might have formally played an important role in the concentrations of metabolic products in an individual patient (30)—a consideration that was not further explored in this particular study.

Statements

Data availability statement

The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.

Ethics statement

The studies involving human participants were reviewed and approved by the local ethics committee of Beijing You'an hospital, of the Capital Medical University. The patients/participants provided their written informed consent to participate in this study.

Author contributions

JS, TZ, and YonghongZ carried out the studies, participated in collecting data, and drafted the manuscript. YananZ performed the statistical analysis. LQ participated in its design. KL, HS, and YanZ helped to draft the manuscript. All authors read and approved the final manuscript.

Funding

The work was supported by CAMS Innovation Fund for Medical Sciences (CIFMS, Grant No. 2016-I2M-1-018), National S & T Major Project for Infectious Diseases Control (Grant No. 2017ZX10201101001008), the Beijing Municipal Administration of Hospitals Ascent Plan (Grant No. DFL20181701), and the Biomarkers of infection related diseases Beijing key laboratory (Grant No. BZ0373).

Acknowledgments

The authors thank all participants for providing blood samples that were used in this study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2019.01069/full#supplementary-material

References

  • 1.

    ChenWSunKZhengRZengHZhangSXiaCet al. Cancer incidence and mortality in China, 2014. Chin J Cancer Res. (2018) 30:112. 10.21147/j.issn.1000-9604.2018.01.01

  • 2.

    El-SeragHB. Epidemiology of viral hepatitis and hepatocellular carcinoma. Gastroenterology. (2012) 142:126473.e1261. 10.1053/j.gastro.2011.12.061

  • 3.

    RingelhanMO'ConnorTProtzerUHeikenwalderM. The direct and indirect roles of HBV in liver cancer: prospective markers for HCC screening and potential therapeutic targets. J Pathol. (2015) 235:35567. 10.1002/path.4434

  • 4.

    ParkinDM. The global health burden of infection-associated cancers in the year 2002. Int J Cancer. (2006) 118:303044. 10.1002/ijc.21731

  • 5.

    HongTPGowPJFinkMDevARobertsSKNicollAet al. Surveillance improves survival of patients with hepatocellular carcinoma: a prospective population-based study. Med J Aust. (2018) 209:34854. 10.5694/mja18.00373

  • 6.

    SauzayCPetitABourgeoisAMBarbareJCChauffertBGalmicheAet al. Alpha-foetoprotein (AFP): a multi-purpose marker in hepatocellular carcinoma. Clin Chim Acta. (2016) 463:3944. 10.1016/j.cca.2016.10.006

  • 7.

    ManiniMASangiovanniAFornariFPiscagliaFBiolatoMFanigliuloLet al. Clinical and economical impact of 2010 AASLD guidelines for the diagnosis of hepatocellular carcinoma. J Hepatol. (2014) 60:9951001. 10.1016/j.jhep.2014.01.006

  • 8.

    MarreroJALokAS. Newer markers for hepatocellular carcinoma. Gastroenterology. (2004) 127:S1139. 10.1053/j.gastro.2004.09.024

  • 9.

    MeguroMMizuguchiTNishidateTOkitaKIshiiMOtaSet al. Prognostic roles of preoperative alpha-fetoprotein and des-gamma-carboxy prothrombin in hepatocellular carcinoma patients. World J Gastroenterol. (2015) 21:493345. 10.3748/wjg.v21.i16.4933

  • 10.

    WangXZhangASunH. Power of metabolomics in diagnosis and biomarker discovery of hepatocellular carcinoma. Hepatology. (2013) 57:20727. 10.1002/hep.26130

  • 11.

    RessomHWXiaoJFTuliLVargheseRSZhouBTsaiTHet al. Utilization of metabolomics to identify serum biomarkers for hepatocellular carcinoma in patients with liver cirrhosis. Anal Chim Acta. (2012) 743:90100. 10.1016/j.aca.2012.07.013

  • 12.

    LuoPYinPHuaRTanYLiZQiuGet al. A Large-scale, multicenter serum metabolite biomarker identification study for the early detection of hepatocellular carcinoma. Hepatology. (2018) 67:66275. 10.1002/hep.29561

  • 13.

    ChenTXGWangXFanJQiuYZhengXQiXet al. Serum and urine metabolite profiling reveals potential biomarkers of human hepatocellular carcinoma. Mol Cell Proteomics. (2011) 10:M110.004945. 10.1074/mcp.M110.004945

  • 14.

    LiverEAftSot. EASL clinical practice guidelines: management of hepatocellular carcinoma. J Hepatol. (2018) 69:182236. 10.1016/j.jhep.2018.03.019

  • 15.

    TerraultNALokASFMcMahonBJChangKMHwangJPJonasMMet al. Update on prevention, diagnosis, and treatment of chronic hepatitis B: AASLD 2018 hepatitis B guidance. Hepatology. (2018) 67:156099. 10.1002/hep.29800

  • 16.

    DuHWangLLiuBWangJSuHZhangTet al. Analysis of the metabolic characteristics of serum samples in patients with multiple myeloma. Front Pharmacol. (2018) 9:884. 10.3389/fphar.2018.00884

  • 17.

    XueRLinZDengCDongLLiuTWangJet al. A serum metabolomic investigation on hepatocellular carcinoma patients by chemical derivatization followed by gas chromatography/mass spectrometry. Rapid Commun Mass Spectrom. (2008) 22:30618. 10.1002/rcm.3708

  • 18.

    GuriYColombiMDazertEHindupurSKRoszikJMoesSet al. mTORC2 promotes tumorigenesis via lipid synthesis. Cancer Cell. (2017) 32:80723.e812. 10.1016/j.ccell.2017.11.011

  • 19.

    BalabanSLeeLSSchreuderMHoyAJ. Obesity and cancer progression: is there a role of fatty acid metabolism?Biomed Res Int. (2015) 2015:274585. 10.1155/2015/274585

  • 20.

    Vander HeidenMGCantleyLCThompsonCB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. (2009) 324:102933. 10.1126/science.1160809

  • 21.

    QiuJXiaoJHanCLiNShenXJiangHet al. Potentiation of tumor necrosis factor-alpha-induced tumor cell apoptosis by a small molecule inhibitor for anti-apoptotic protein hPEBP4. J Biol Chem. (2010) 285:122417. 10.1074/jbc.M110.111898

  • 22.

    Zijian WuBLXueminZhengHuijingHouYingLi. Role of the PEBP4 protein in the development and metastasis of gastric cancer. Oncotarget. (2017) 8:1817784. 10.18632/oncotarget.15255

  • 23.

    YuGHuangBChenGMiY. Phosphatidylethanolamine-binding protein 4 promotes lung cancer cells proliferation and invasion via PI3K/Akt/mTOR axis. J Thorac Dis. (2015) 7:180616. 10.3978/j.issn.2072-1439.2015.10.17

  • 24.

    McMillanDC. Systemic inflammation, nutritional status and survival in patients with cancer. Curr Opin Clin Nutr Metab Care. (2009) 12:2236. 10.1097/MCO.0b013e32832a7902

  • 25.

    FearonKCVossACHusteadDSCancer Cachexia Study Group. Definition of cancer cachexia_ effect of weight loss, reduced food intake, and systemic inflammation on functional status and prognosis. Am J Clin Nutr. (2006) 83:134550. 10.1093/ajcn/83.6.1345

  • 26.

    ZhangCChenXLiuHLiHJiangWHouWet al. Alpha fetoprotein mediates HBx induced carcinogenesis in the hepatocyte cytoplasm. Int J Cancer. (2015) 137:181829. 10.1002/ijc.29548

  • 27.

    WangXWangQ. Alpha-fetoprotein and hepatocellular carcinoma immunity. Can J Gastroenterol Hepatol. (2018) 2018:9049252. 10.1155/2018/9049252

  • 28.

    JiangYSunAZhaoYYingWSunHYangXet al. Proteomics identifies new therapeutic targets of early-stage hepatocellular carcinoma. Nature. (2019) 567:25761. 10.1038/s41586-019-0987-8

  • 29.

    Juarez-HernandezEMotola-KubaDChavez-TapiaNCUribeMBarbero BecerraV. Biomarkers in hepatocellular carcinoma: an overview. Expert Rev Gastroenterol Hepatol. (2017) 11:54958. 10.1080/17474124.2017.1311785

  • 30.

    SenguptaSParikhND. Biomarker development for hepatocellular carcinoma early detection: current and future perspectives. Hepat Oncol. (2017) 4:11122. 10.2217/hep-2017-0019

Summary

Keywords

hepatocellular carcinomas, metabolomics, alpha-fetoproteins, biomarkers, HBV

Citation

Sun J, Zhao Y, Qin L, Li K, Zhao Y, Sun H, Zhang T and Zhang Y (2019) Metabolomic Profiles for HBV Related Hepatocellular Carcinoma Including Alpha-Fetoproteins Positive and Negative Subtypes. Front. Oncol. 9:1069. doi: 10.3389/fonc.2019.01069

Received

05 July 2019

Accepted

30 September 2019

Published

15 October 2019

Volume

9 - 2019

Edited by

Marina Pajic, Garvan Institute of Medical Research, Australia

Reviewed by

Sean Porazinski, The Kinghorn Cancer Centre, Australia; Ashleigh Louise Parkin, Garvan Institute of Medical Research, Australia

Updates

Copyright

*Correspondence: Ting Zhang Yonghong Zhang

This article was submitted to Gastrointestinal Cancers, a section of the journal Frontiers in Oncology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics