PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 45913.
Published online 2017 April 5. doi:  10.1038/srep45913
PMCID: PMC5381110

Quantitative proteomics by SWATH-MS reveals sophisticated metabolic reprogramming in hepatocellular carcinoma tissues

Abstract

Hepatocellular carcinoma (HCC) is one of the most common cancers worldwide, and understanding its molecular pathogenesis is pivotal to managing this disease. Sequential window acquisition of all theoretical mass spectra (SWATH-MS) is an optimal proteomic strategy to seek crucial proteins involved in HCC development and progression. In this study, a quantitative proteomic study of tumour and adjacent non-tumour liver tissues was performed using a SWATH-MS strategy. In total, 4,216 proteins were reliably quantified, and 338 were differentially expressed, with 191 proteins up-regulated and 147 down-regulated in HCC tissues compared with adjacent non-tumourous tissues. Functional analysis revealed distinct pathway enrichment of up- and down-regulated proteins. The most significantly down-regulated proteins were involved in metabolic pathways. Notably, our study revealed sophisticated metabolic reprogramming in HCC, including alteration of the pentose phosphate pathway; serine, glycine and sarcosine biosynthesis/metabolism; glycolysis; gluconeogenesis; fatty acid biosynthesis; and fatty acid β-oxidation. Twenty-seven metabolic enzymes, including PCK2, PDH and G6PD, were significantly changed in this study. To our knowledge, this study presents the most complete view of tissue-specific metabolic reprogramming in HCC, identifying hundreds of differentially expressed proteins, which together form a rich resource for novel drug targets or diagnostic biomarker discovery.

Liver cancer is one of the most common malignant cancers in the world, with more than 850,000 new cases worldwide annually1. This neoplasm is currently the second leading cause of cancer-related death globally, and the incidence is increasing2. Among all primary liver cancers, hepatocellular carcinoma (HCC) is the most common neoplasm, accounting for approximately 90% of all cases1,3,4,5,6,7,8. Hepatitis B virus (HBV) infection, hepatitis C virus (HCV) infection, alcohol abuse and intake of aflatoxin B1 are the main factors contributing to HCC1,3,4,5,6,7. In China, HCC has been ranked as the second most frequent fatal cancer since the 1990s9, and the majority of HCCs in China are caused by HBV infection10,11.

Currently, surgical resection and liver transplantation are considered the best treatment options for early-stage HCC and are curative therapies for approximately 30% to 40% of early-stage patients3,12. Due to the asymptomatic features of HCC at early stages, patients are often diagnosed at very advanced stages. Thus, there is an urgent need to find key carcinogenesis-associated molecules for HCC diagnosis and treatment.

Mass spectrometry (MS)-based proteomic analysis of human clinical tissues is a powerful tool to investigate cancer biomarkers and therapeutic targets13. Numerous clinical studies of HCC have been reported over the past decade using various quantitative techniques14,15,16,17,18,19,20, including SILAC (stable isotope labelling by amino acids in cell culture), iTRAQ (isobaric tags for relative and absolute quantification) and CDIT (culture-derived isotope tags) labelling techniques as well as label-free proteomics approaches based on quantification by ion intensity or spectral counting. Label-free approaches are relatively cheap compared to labelling approaches; when labelling reagents are not required, high-throughput and sensitive analyses in a mass spectrometer are possible. Quantitative studies of HCC using spectral counting and ion intensities have also been reported19,20. SWATH-MS (sequential window acquisition of all theoretical mass spectra) is an emerging label-free quantification approach that combines a highly specific data independent acquisition (DIA) method with a novel targeted data extraction strategy to mine the resulting fragment ion data sets. SWATH-MS has been widely used to compare protein expression and modify alterations21,22,23,24. To our knowledge, no SWATH-MS approach has been used to study HCC proteomics until now.

In this study, we compared the protein expression of tumourous (HCC) and adjacent non-tumourous (non-HCC) tissues from 14 HBV-associated HCC patients using a SWATH-MS technique to identify new HCC biomarkers and potential therapeutic target candidates. In total, 338 differential proteins were quantified, and most down-regulated proteins were involved in metabolism. Sophisticated reprogramming of cell metabolic pathways was revealed. These observations are essential to elucidate the mechanisms underlying the occurrence and progression of HCC and contribute to the discovery of new candidates for early HCC diagnosis.

Results

Differentially expressed proteins quantified by SWATH-MS analysis in HCC tissues

The experimental scheme of the present study is shown in Fig. 1. HCC and non-HCC liver tissue samples were compared by SWATH-MS to identify differentially expressed proteins that can be used as biomarkers for HCC diagnosis or in HCC development and progression. To avoid individual differences and detect true HCC-related proteins, samples were analysed by equal pooling of two or three tissues from both groups to determine a quantitative expression ratio between HCC and non-HCC liver tissue groups based on total ion intensity normalization. Five biological replicates were analysed, and 14 pairs of tissue samples were used in total (Supplementary Table S1).

Figure 1
Quantitative proteomic workflow of human HCC and adjacent non-tumourous liver tissues analysed using a SWATH-MS approach.

The targeted identification of peptides in SWATH-MS datasets requires a priori generation of a spectral library that includes essential coordinates for each targeted peptide, such as precursor ion masses, fragment ion masses, fragment ion intensities and retention times21. For each biological replicate, a spectral library was generated with a traditional data-dependent acquisition (DDA) mass spectrometry technique as described in the Methods section. Five libraries were obtained in total. On average, the spectral libraries contained approximately 2,491 distinct protein groups, and 26,586 peptides were identified with greater than 99% confidence and passed the global false discovery rate (FDR) from fit analysis using a critical FDR of 1% (Supplementary Table S2, Supplementary Figure S1). Taken together, these findings indicated that the experimentally generated spectral libraries contained only high-confidence proteins.

Following generation of the spectral library, the identification and quantification of HCC and non-HCC proteins were performed using a SWATH-MS approach as described by Gill et al.21, with modifications. Proteins were digested by trypsin, and the peptides were separated using a gradient of 120 min on a reverse-phase nanoLC instrument. SWATH data from six injections for each biological replicate were submitted in unison to PeakView software (Version 1.2, AB Sciex) for targeted data extraction, which resulted in the quantitative export of 2,122, 2,681, 1,860, 2,115 and 2,346 unique proteins for the analysis of five biological replicates (Fig. 2a). In total, 4,216 proteins were quantified in at least one biological replicate, and 1,903 proteins were quantified in at least three biological replicates, which accounted for 45% of all quantified proteins (Fig. 2b and Supplementary Table S3).

Figure 2
Quantitative proteomic analysis between HCC and non-HCC tissues.

Relative protein quantification was analysed using MarkerView (Version 1.2.1, AB Sciex) and R (Version 3.3.1, the R foundation), as described in the Methods section. Normality distribution of each technical replicates was performed after log2 transformation for the peak intensities of all the MS measurements prior to further analysis and the histograms were shown in Supplementary Figure S2. In total, four criteria were used to filter out differential proteins. First, the Shapiro-Wilk test was used to test normality for each protein within one biological replicate, and only proteins that met normality were used for further analysis. Second, Welch’s t-test and Benjamini-Hochberg multiple test correction were performed, and an adjusted p_value < 0.05 and fold change (FC)  1.5 or FC  1/1.5 was considered statistically significant. With these two criteria, 865, 884, 550, 715 and 760 differentially expressed proteins were obtained for the five biological replicates (Fig. 2a). Third, up- or down- regulated proteins were detected in at least three biological replicates are shown in the Venn diagram (Fig. 2c). Fourth, the average ratio of up- or down-regulated proteins had to meet the 1.5-FC cutoff requirements. With the above mentioned four criteria, 191 up-regulated and 147 down-regulated proteins were obtained in total and are shown in the heatmap (Fig. 2d). Supplementary Table S4 and Table S5 present the differentially expressed proteins for the adjusted p_value < 0.05 and FC ≥1.5 or FC  1/1.5 in five biological replicates and the final analysis results.

To evaluate our SWATH-MS data, comprehensive proteomics-based quantitative expression data obtained from five previously published studies of HCC and adjacent/normal tissues were used to perform side-by-side comparisons15,17,18,20,25. The differentially expressed proteins quantified in these five studies were 151, 573, 71, 267 and 648. The comparison revealed that 49 (42), 155 (150), 28 (26), 68 (65) and 70 (62) differentially expressed proteins detected in our data were present in the five previous studies; the trends of the major proteins were consistent with published data (numbers in brackets represent proteins with the same change trend) (Supplementary Table S6). In total, 214 (199) proteins in our dataset were detected in previous studies, which accounted for 63.31% (58.88%) of the 338 differential proteins. Among the 214 proteins, 99 proteins were detected in ≥3 studies, including our data (Supplementary Table S6). Among these 99 frequently identified common proteins, 41 proteins were up-regulated, and 58 were down-regulated. Among the 338 differentially expressed proteins, nine were validated by western blotting in previous studies, and five were validated by immunohistochemical (IHC) methods. Seven of the nine western blotting-validated proteins (OTC, PEBP1, CPS1, BHMT, CLIC1, PPA1 and APEX1) and all five IHC-tested proteins (OTC, BHMT, CLIC1, PPA1 and APEX1) were detected in ≥3 data sets. All five proteins (OTC, BHMT, CLIC1, PPA1 and APEX1) were tested using both methods as HCC potential biomarkers. Our SWATH data were consistent with the western blotting- and IHC-validated results (Supplementary Table S6).

To decipher whether the differentially expressed proteins were detectable in plasma, we searched for these proteins in the Human Plasma Proteome database (HPPD)26. Approximately 85.50% of the differentially expressed proteins (289/338) appeared in this database, suggesting that they had relatively strong potential to be secreted into the blood (Supplementary Table S7). Among these 289 proteins, the numbers of up- and down-regulated proteins were 169 and 120, respectively. Furthermore, 94 proteins were detected in ≥3 studies and were present in HPPD; among these, 38 proteins were up-regulated (Supplementary Table S7). These proteins were potential biomarker candidates for HCC; some may be tested using multiple reaction monitoring (MRM), as described by Hou et al.27. The above results show the reliability of our SWATH-MS proteomic results and indicate that these differentially expressed proteins may be useful to delineate HCC properties and screen HCC biomarker candidates.

GO and KEGG pathway enrichment analysis

To investigate the function of these differentially expressed proteins, GO and KEGG pathway analyses of up-regulated and down-regulated proteins were performed separately by DAVID (Version 6.8, LHRI & DAVID Bioinformatics)28,29. This method easily determined the characteristics of up- or down-regulated proteins. Liver proteins downloaded from the human proteome map (http://humanproteomemap.org) were used as the background dataset for enrichment analysis30. For up-regulated proteins, 34 significant enrichments were identified using GO analysis (p < 0.05, p_values were corrected using the Benjamini-Hochberg procedure). These were classified into three GO categories, including biological processes (BP, 5), molecular functions (MF, 6) and cellular components (CC, 23) (Table 1). In BP, those items significantly participated in cell-cell adhesion (p = 7.12 × 10−5), mRNA splicing via the spliceosome (p = 4.60 × 10−3), SRP-dependent cotranslational protein targeting to membrane (p = 6.97 × 10−3) and translation initiation (p = 9.07 × 10−3). The most significant terms for MF and CC were poly(A) RNA binding (p = 5.60 × 10−20) and extracellular exosome (p = 1.30 × 10−15), respectively. For the down-regulated proteins, 66 significant enrichments were obtained, including 29 BP, 28 MF and 9 CC (Table 1). The five most significant terms of BP were related to metabolic processes, including oxidation-reduction processes (p = 4.93 × 10−21), xenobiotic metabolic processes (p = 1.03 × 10−8), metabolic processes (p = 1.95 × 10−8), drug metabolic processes (p = 3.32 × 10−8) and epoxygenase P450 pathway (p = 7.15 × 10−7). For MF and CC, the most significant terms were oxidoreductase activity (p = 1.03 × 10−16) and mitochondrial matrix (p = 1.69 × 10−18), respectively.

Table 1
GO and KEGG pathway analyses of up-regulated and down-regulated proteins by DAVID.

To locate the key pathways implicated in HCC development and progression, KEGG pathway enrichment analysis was performed for 338 differentially expressed proteins. As for GO analysis, up- and down-regulated proteins were analysed separately, and the background of the enrichment used liver proteins from the human proteome map30. Pathway analysis showed that the up-regulated proteins were significantly enriched in term of spliceosome (p = 4.24 × 10−2). For down-regulated proteins, 37 terms were enriched, and the most significant terms were metabolic pathways (p = 3.36 × 10−42), glycine, serine and threonine metabolism (p = 9.03 × 10−13), retinol metabolism (p = 1.65 × 10−12), drug metabolism - cytochrome P450 (p = 6.72 × 10−12) and biosynthesis of antibiotics (p = 7.05 × 10−11). The all five significant GO terms of the BP for up-regulated proteins and the top 10 most significant GO terms of the BP and KEGG pathways for down-regulated proteins are shown in Fig. 3(a–c). Those of MF and CC are in Supplementary Figure S3(a–d). All enriched terms are shown in Supplementary Table S8.

Figure 3
GO and KEGG pathway enrichment of 338 significantly regulated proteins according to DAVID functional annotation.

Western blot validation for nine selected proteins in clinical HCC tissues

Nine candidate proteins were selected for validation by western blot analysis using six sample pairs of HCC and non-HCC liver tissues that differed from those used in the proteomics studies (Supplementary Table S1 patient ID 15–20). The candidate proteins were selected based on either dramatic fold change or involvement in key metabolic pathways. Five candidate proteins were up-regulated, namely FBXO2, ACSL4, PLIN2, PKM2 and GFPT1. The four down-regulated proteins were CYP1A2, FTCD, UGT2B7 and PCK2. Here, the analysis showed differential expression of all the candidates in HCC tissues compared with non-HCC tissues. FBOX2, ACSL4 and PLIN2 showed strong expression in all six tumour samples but weak or no expression in non-tumour tissues. PKM2 and GFPT1 showed generally high expression levels in five HCC tissues compared with the control group. For all down-regulated proteins, low expression was detected in HCC tissues compared with non-HCC tissues (Fig. 4). The representative extracted ion chromatogram (XIC) comparisons in these nine proteins are shown in Supplementary Figure S4. Overall, the western blot analysis results for all nine proteins were consistent with the proteomics data, which indicated that our proteomics data were highly reliable and that some proteins are worthy of further investigation.

Figure 4
Validation of nine selected proteins in clinical HCC tissues by western blotting.

Discussion

The purpose of this study was to characterize proteomic changes in HCC tissues, provide potential protein candidates for biomarker discovery and suggest molecular mechanisms of HCC development and progression. Although much proteomic research has been performed, the biological mechanisms of HCC development and progression are still unclear.

Metabolic reprogramming is a hallmark of cancer31. Cancer cells can increase the amount of glucose and glutamine to satisfy energy needs and macromolecular synthesis demands. Therefore, understanding the metabolism of tumours remains an intense study topic with important therapeutic potential32.

Using the newly developed SWATH-MS technique, we quantified more than 4,000 proteins, and 338 proteins were differentially expressed in HCC. Sophisticated metabolic reprogramming was revealed as depicted in Fig. 5, including the following major aspects.

Figure 5
Sophisticated metabolic reprogramming in HCC.

First, the oxidative pentose phosphate pathway (PPP) was up-regulated in HCC. PPP is the first branch pathway of glycolysis. In PPP, glucose-6-phosphate becomes partially oxidized to generate NADPH and ribose-5-phosphate. PPP is frequently elevated in tumourigenesis. In our study, two key enzymes—the rate-limiting enzyme glucose-6-phophate dehydrogenase (G6PD) and transaldolase (TALDO)—were over-expressed in HCC. The over-expression of G6PD and TALDO was detected in previous HCC references18,25,33,34.

Second, the serine, glycine and sarcosine biosynthesis/metabolism pathways were significantly down-regulated in HCC. Serine biosynthesis is a key metabolic pathway for cell proliferation, contributing carbon to many anabolic products, such as protein, glutathione, nucleotide and phospholipid biosynthesis35,36. Phosphoglycerate dehydrogenase (PHGDH) and phosphoserine aminotransferase (PSAT) are two key enzymes of serine biosynthesis37,38. In our study, both PHGDH and PSAT were down-regulated in HCC tissues. The two proteins were quantified in the HCC proteomics study mentioned above, with the same trends. Overall, our data indicate that the serine biosynthesis pathway is down-regulated in HCC.

In line with PHGDH and PSAT, key enzyme of serine/glycine metabolism, serine hydroxymethyl-transferase 1 (SHMT1) in the cytosol was down-regulated. SHMT catalyses the formation of glycine from serine, thereby generating 5,10-methylene-tetrahydrofolate (5,10-MTHF)39. Glycine-N-methyl transferase (GNMT), the enzyme that generates sarcosine from glycine, showed dramatically low expression. In addition to GNMT, sarcosine levels were regulated by sarcosine dehydrogenase (SARDH), the enzyme that converts sarcosine back to glycine and dimethylglycine dehydrogenase (DMGDH), which generates sarcosine from dimethylglycine40. Both SARDH and DMGDH were down-regulated in HCC tissues. In addition, betaine-homocysteine S-methyltransferase 1 (BHMT1) converts betaine and homocysteine into dimethylglycine and methionine. In our study, BHMT1 was also down-regulated.

In addition to generating sarcosine, the glycine cleavage system is a catabolic mechanism for glycine. Three key glycine catabolism enzymes—glycine decarboxylase (GLDC), glycine cleavage system H protein (GCSH) and glycine N-acyltransferase (GLYAT)—were also decreased in HCC tissues. GLDC is a key component of the highly conserved glycine cleavage system in amino acid metabolism, which catalyses the breakdown of glycine to form CO2, NH3 and 5,10-MTHF to fuel one-carbon metabolism41. GLDC and GCSH are two members of the glycine cleavage system, and their down-regulation indicated the significantly low expression of the system. In total, 10 serine, glycine and sarcosine metabolic enzymes were down-regulated in our study. The schematic pathway and changes in these enzymes are shown in Fig. 6.

Figure 6
Down-regulation of serine, glycine and sarcosine biosynthesis/metabolism pathways in HCC.

PHGDH is the first enzyme in serine biosynthesis from a glycolic intermediate, is frequently amplified in breast cancer and is required for the growth of PHGDH-amplified cells in vitro37. The over-expression of GLDC in the tumour-initiating cells (TIC) of non-small cell lung cancer (NSCLC) was also reported to induce dramatic changes in glycolysis and glycine metabolism42. Sarcosine was identified as a differentially expressed metabolite that was highly elevated during prostate cancer progression to metastasis and is a potentially important metabolic intermediate of cancer cell invasion and aggression40. These studies showed that the up-regulation of serine/glycine metabolism may contribute to the pathogenesis of several human cancers and may provide novel targets for improving anticancer therapies.

Overall, our observation was the opposite of earlier results, which may be due to different cancer types. GLDC, together with GCSH and eight other enzymes related to glycine metabolism, mainly catabolism, were down-regulated in HCC. This indicated that glycine catabolism plays different roles in HCC other than NSCLC, and the function of glycine metabolism in HCC is worthy of further investigation. The down-regulation of the glycine catabolism pathway agreed well with the increased glycine in HCC serum compared with normal control samples in a recent metabolomics study43.

Third, fatty acid synthesis (FAS) was up-regulated in HCC. Fatty acid de novo synthesis is required for membrane synthesis and therefore for cell growth and proliferation. It has been revealed that fatty acid synthesis requires several key enzymes, including pyruvate dehydrogenase (PDH) complex, citrate synthase (CS), ATP-citrate synthase (ACLY), acetyl-CoA carboxylase (ACC) and fatty acid synthase (FASN)44. Three of these, ACLY, ACC and FASN, were frequently up-regulated in transformed cells45. In our study, one subunit of PDH, CS, ACLY and FASN were all over-expressed in HCC. The over-expression of these four enzymes, especially the dramatic increase in ACLY and FASN, indicated the up-regulation of fatty acid biosynthesis in HCC. An increased capacity for producing lipids de novo facilitates the formation of lipid bilayers but also enables the cell to adapt to oxidative stress46. We speculate that similar mechanisms exist in HCC for fatty acid synthesis as in other cancers.

Moreover, fatty acid oxidation (FAO, also known as β-oxidation) was down-regulated in HCC. In most cases, the growth and survival of cancer cells is limited by levels of cytosolic NADPH. The production of FAO-derived cytosolic NADPH by cancer cells is critical to counteract oxidative stress47. In our study, 6 enzymes of FAO were quantified. Levels of all of these enzymes were decreased, including ACADS and ACADSB (which catalyse acyl-CoA dehydrogenase activity); ECHS1, ECHD2 and ECHD3 (which catalyse enoyl-CoA hydratase activity); and ACAA2 (which catalyse acetyl-CoA C-acyltransferase activity). The schematic diagrams of FAO and all altered enzymes are shown in Supplementary Figure S5. Evidence indicates that some tumours, including prostate tumours, leukaemia and large B-cell lymphoma, utilize FAO as their main energy supply for proliferation and survival48. The down-regulation of FAO in HCC indicated that FAO plays different roles in HCC, unlike other tumours.

In addition to the above mentioned pathways, glycolysis and gluconeogenesis are key metabolic pathways in human cells. Conversion of glucose to lactate is a major pathway of glucose metabolism in cancer cells even in the presence of O2; this is known as the Warburg effect or aerobic glycolysis49. Pyruvate kinase (PK) catalyses the final rate-limiting reaction in glycolysis by transferring high-energy phosphate from phosphoenolpyruvate (PEP) to ADP to produce ATP and pyruvate50. PKM2 is the dominant M isoform in most adult tissues and is the major PK in proliferating and cancer cells50. Elevated PKM2 expression has been demonstrated in various human tumours, including lung, breast, prostate, blood, cervix, kidney, bladder and colon, when compared to matched normal tissues51. PKM2 promotes tumourigenesis by regulating the Warburg effect, as reported in a previous study50. In our proteomic data, PKM2 was over-expressed in HCC tissue, and the expression of PKM2 was also validated by western blotting. Thus, our results indicate that the change in PKM2 in HCC is consistent with the Warburg effect.

As shown in Fig. 5, pyruvate can be transformed into two different products, acetyl-CoA and oxaloacetate, by PDH and pyruvate carboxylase (PC). PDHA1, a pyruvate dehydrogenase E1 component subunit alpha, was over-expressed, while PC was down-regulated in our proteomic data set. The down-regulation of PC in HCC has been reported20. Thus, our data indicated that the increase in pyruvate was directed to acetyl-CoA for the TCA cycle, which may further be used for fatty acid synthesis.

Phosphoenolpyruvate carboxykinase (PCK) is the rate-limiting enzyme of gluconeogenesis in the liver and kidney. There are two isoforms of PCK: cytosolic PCK (PCK1) and a mitochondrial isoform of PCK (PCK2). PCK1, but not PCK2, was reported to be over-expressed in colorectal cancer and promotes tumour growth by increasing glucose and glutamine metabolism52. PCK2 expression was reported to be elevated in non-small-cell lung carcinoma (NSCLC) and is regulated by glucose and required for in vivo tumour growth53. The down-regulation of PCK2 in HCC has been documented in a previous study17. In our proteomic data set, both PCK1 and PCK2 were lowly expressed in HCC. The down-regulation of PCK2 was validated using western blot analysis and was consistent with the proteomic data. Thus, the down-regulation of PCK was a specific characteristic of HCC, different from colorectal cancer and NSCLC in the literature. This observation is in agreement with the view that cancer cells display metabolic flexibility, and not all features are universal53. Taken together, these data demonstrate roles for PKM2, PC and PCK2 in HCC, which links metabolic flux and anabolic pathways to cancer cell proliferation.

In summary, we compared the proteomic profile of HCC tumour tissues (n = 14) with paired adjacent non-tumourous liver tissues (n = 14) using a SWATH-MS quantification strategy. To obtain reliable results, five biological replicates were performed, and strict filtration criteria were applied. In total, 338 differentially expressed proteins were obtained, and these proteins were involved in various key pathways. Compared with previous studies, more than 50% of proteins were quantified by others with the same change trend. Approximately 85% of the differentially expressed proteins showed relatively strong potential to be secreted into the blood. It is noteworthy that sophisticated metabolic reprogramming in HCC was revealed by proteomic data, including the up-regulation of glycolysis, the pentose phosphate pathway and fatty acid biosynthesis and the down-regulation of gluconeogenesis; serine, glycine and sarcosine metabolism; and fatty acid-oxidation. This reprogramming showed the distinct metabolic flexibility of HCC. In total, 27 key metabolic enzymes were quantified in our study, including PCK2, PDH and G6PD, which are important for cancer development and progression. In addition to metabolic reprogramming, spliceosome pathway was significantly up-regulated. The results indicate that the SWATH-MS strategy is effective for identifying crucial proteins involved in HCC development and progression. Differentially expressed proteins form a rich resource for diagnostic biomarkers or therapeutic drug target discovery. Furthermore, our findings into HCC-specific metabolic reprogramming may provide new insights into understanding HCC biology and direct optimal treatment.

Materials and Methods

Patients and Clinical Specimens

The following methods were carried out in accordance with approved guidelines. Access to human tissues complied with the guidelines of the Ethics Committee at Beijing 302 Hospital, and informed consent was obtained from all patients for collecting the specimens. All experimental protocols were approved by the National Center of Biomedical Analysis. HCC tissues and adjacent non-tumourous liver tissue counterparts used for this study were collected from 20 HBV-associated HCC patients who underwent hepatectomy at Beijing 302 Hospital between 2013 and 2014 (14 males and 6 females). The tumourous and non-tumourous liver tissues were determined by the experienced pathologist and pathological sections. None of these patients received anti-neoplastic therapy prior to surgery. All patients used for the quantitative proteomics study and western blotting verification had liver cirrhosis. Patients details are shown in Supplementary Table S1. Tissues were immediately snap-frozen in liquid nitrogen after surgical resection and stored at −80 °C until use.

Reagents and Materials

Hepes, NaCl, Urea and thiourea were purchased from Sigma-Aldrich (St. Louis, MO, USA). Cocktail was purchased from Roche Diagnostics (Indianapolis, IN, USA). Trypsin were purchased from Promega (Madison, WI, USA). Protein Assay Dye Reagent Concentrate was purchased from Bio-Rad (Hercules, CA, USA). Dithiothreitol (DTT) was purchased from Amresco (Solon, OH, USA). Iodoacetamide (IAA) was purchased from Acros Organics (Morris Plains, NJ, USA). Acetonitrile (ACN, HPLC grade) and methanol were purchased from Fisher Scientific (Fair Lawn, NJ, USA). Formic acid (FA) was provided by Fluka (Milwaukee, WI, USA). SuperSignal West Pico was purchased from Thermo Scientific (Rockford, IL, USA). All water used in the experiments was purified using a Milli-Q system (Millipore, Billerica, MA, USA).

Protein Extractions and Trypsin Digestion

For each extraction, ~0.2 g of tissue was ground into powder in liquid nitrogen with a pre-cooled mortar and pestle. Samples were homogenized on ice in 1 ml of lysis buffer (50 mM Hepes, 6 M urea, 2 M thiourea and 1 × protease inhibitor cocktail). The samples were lysed by sonication for 1 min on ice (pulse on 3 s, pulse off 10 s). After centrifugation at 14,000 × g for 30 min at 4 °C, the supernatant was collected, and the protein concentration was determined by the Bradford method54 (Bradford Protein Assay, Bio-Rad). Twenty μg proteins from each sample were digested with trypsin using filter-aided sample preparation (FASP) as previously described55,56. After digestion, the peptides were dried in a vacuum for MS analysis.

Mass Spectrometry Analysis

MS analysis was performed using an AB Sciex 5600+ TripleTOF mass spectrometer (Concord, Ontario, Canada) interfaced to an EkspertTM NanoLC 425 system (Dublin, CA) as previously described21,22,23,24. For library construction, peptides were trapped on a NanoLC pre-column (Chromxp C18-LC-3 μm, size 0.35 × 0.5 mm, Eksigent), eluted onto an analytical column (C18-CL-120, size 0.075 × 150 mm, Eksigent) and separated by a 120-min gradient from 5 to 35% Buffer B (Buffer A: 2% ACN, 98% H2O, Buffer B: 98% ACN, 2% H2O, 0.1% FA) at a flow rate of 300 nL/min. Full-scan MS was performed in positive ion mode with a nano-ion spray voltage of 2.3 kv from 350 to 1500 (m/z), with up to 50 precursors selected for MS/MS (m/z 100–1500). The selection criteria for parent ions included an intensity greater than 150 counts/s, a charge state from +2 to +5, a mass tolerance of 50 mDa and dynamic exclusion for 15 s. Ions were fragmented in the collision cell using rolling collision energy.

In SWATH™ acquisition, the parameters were essentially the same as those described by Gillet et al.21. With the same chromatographic conditions used in the DDA run described above and a variable isolation window obtained using variable window package software (including 1 Da for the window overlap), a set of 55 overlapping windows was constructed, covering the precursor mass range of 400–1250 Da. The collision energy for each window was determined based on the appropriate collision energy for a 2+ ion centred in the window with a spread of 15 eV. The high-sensitivity mode was used, allowing accurate extraction of the fragment ion masses.

Generating the Reference Spectral Library

For each bio-replicate analysis, three DDA injections were performed to increase protein coverage as described by Hou et al.27. All three mass spectrometry files were searched in unison using ProteinPilot software (Version 4.2, AB Sciex) with the Paragon algorithm as described by Haverland et al. with minor modifications22. Samples were input as unlabelled samples with the following parameters: iodoacetamide-cysteine alkylation, digestion with trypsin and no special factors. The search was conducted using a thorough identification effort and the human UniProt database (April 2013 release)57.

SWATH-MS Data Analysis

Spectral alignment and targeted data extraction of DIA samples were performed with the SWATH Processing Micro App in Peakview (Version 1.2, AB Sciex) using the reference spectral library generated above as described by Haverland et al. with modifications22. There was a reference library for each bio-replicate comparison. Six DIA raw files in one comparison group were loaded in unison using an extraction window of 15 min and the following parameters: 5 peptides, 8 transitions and peptide confidence of >99%, including shared peptides and XIC width set at 50 ppm. After data processing, three distinct files were exported for further quantitation analysis. The processed mrkvw files containing protein information from PeakView were loaded into MarkerView (Version 1.2.1, AB Sciex) for normalization of protein intensity (peak area) for all runs using the built-in total ion intensity sum plug-in. Log2 transformation was performed prior to further statistical analysis. We plot the histogram to check the normality distribution of each technical replicate. Differential analysis was performed using R (Version 3.3.1, the R foundation). In each biological replicate, normality tests of protein expression were performed for HCC/non-HCC groups using the Shapiro-Wilk normality test. Proteins whose expression values met the normality in both groups were retained for Welch’s t-test and the Benjamini-Hochberg multiple test correction. Mean values of protein expression were used for calculation of fold change (FC). Proteins with adjusted p < 0.05 and FC  1.5 or FC  1/1.5 in at least three biological replicates and average FC  1.5 or FC  1/1.5 were regarded as differentially expressed proteins in this study.

Functional Analysis

Protein IDs were converted to gene names using the UniProt Retrieve/ID mapping tool57. Heatmaps for expression of selected proteins in the five groups were created using the R (Version 3.3.1, the R foundation) heatmap package. The DAVID webserver (Version 6.8, LHRI & DAVID Bioinformatics) was employed for Gene Ontology enrichment analysis and KEGG pathway analysis28,29.

Western Blotting Validation

Prior to western blotting, the protein concentration of the samples was determined by the Bradford assay method54. Equal amounts of 20 μg protein per sample were separated by SDS-PAGE on an 8–12% polyacrylamide gel. Proteins were subsequently transferred to PVDF membranes (Immun-Blot PVDF, Bio-Rad, Hercules, CA, USA), and membranes were blocked with 5% (w/v) skim milk in Tris-buffered saline with 0.1% Tween 20 (TBS-T) for 1 h at room temperature. The following antibodies were used: ACSL4 (1:1000), ADPR (1:1000), CYP1A2 (1:1000), FTCD (1:1000), FBXO2 (1:300) and UGT2B7 (1:100) from Proteintech (Rosemont, IL, USA); PKM2 (1:1000) and GFPT1 (1:1000) from Cell Signaling (Danvers, MA, USA); PCK2 (1:1000) from Abcam (Cambridge, MA, USA). These primary antibodies were diluted in 5% (w/v) skim milk in TBS-T and incubated with membranes overnight at 4 °C. After washing for five minutes in TBS-T three times, horseradish peroxidase-labelled secondary antibodies (Jackson ImmunoResearch, West Grove, PA, USA) were used for detection for 1 h at room temperature. Visualization of the immunoreactive proteins was accomplished using enhanced chemiluminescence (SuperSignal West Pico, Thermo, Rockford, IL, USA) followed by exposure to X-ray film (XBT, Carestream, Xiamen, Fujian, China).

Additional Information

Accession codes: All the raw data and meta data were deposited in a public repository iprox database (www.iprox.org) with the accession IPX0000859000, which is based on the National Center for Protein Sciences, Beijing (the PHOENIX) Center).

How to cite this article: Gao, Y. et al. Quantitative proteomics by SWATH-MS reveals sophisticated metabolic reprogramming in hepatocellular carcinoma tissues. Sci. Rep. 7, 45913; doi: 10.1038/srep45913 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supporting Information:
Supplementary Dataset 1:
Supplementary Dataset 2:
Supplementary Dataset 3:
Supplementary Dataset 4:
Supplementary Dataset 5:
Supplementary Dataset 6:
Supplementary Dataset 7:
Supplementary Dataset 8:

Acknowledgments

This work was supported by the National Natural Science Foundation of China (81472273), the Frontier Technology of Beijing Municipal Science and Technology Commission (Z141100000214003) and the National Key Research and Development Programme (2015BAK45B01 and 2016YFC1202400). We would like to thank Dr. Xiaotao Duan (Beijing Institute of Pharmacology and Toxicology) for criticism and help with the manuscript.

Footnotes

The authors declare no competing financial interests.

Author Contributions H.W. and K.H. supervised the study. H.W. and Y.G. conceived and designed the experiments. Y.G., X.W. and Z.S. performed the experiments. D.Y. conducted the sample collection. Z.L., X.Y. and Y.G. analysed the data and prepared the figures. X.W., D.Y., Z.S., F.L., Y.Z., H.W., P.L. and K.H. provided materials and analysis tools. D.Y., Y.G., H.W. and X.Z. provided scientific discussion and guidance. Y.G., H.W. and K.H. wrote the paper. All authors reviewed the manuscript.

References

  • Torre L. A. et al. . Global cancer statistics, 2012. CA: a cancer journal for clinicians 65, 87–108 (2015). [PubMed]
  • Mortality G. B. D. & Causes of Death C. Global, regional, and national age-sex specific all-cause and cause-specific mortality for 240 causes of death, 1990-2013: a systematic analysis for the Global Burden of Disease Study 2013. Lancet 385, 117–171 (2015). [PMC free article] [PubMed]
  • European Association For The Study Of The Liver, European Organisation For Research and Treatment Of Cancer. EASL-EORTC clinical practice guidelines: management of hepatocellular carcinoma. J Hepatol 56, 908–943 (2012). [PubMed]
  • Liu J. & Fan D. Hepatitis B in China. Lancet 369, 1582–1583 (2007). [PubMed]
  • Mohd Hanafiah K., Groeger J., Flaxman A. D. & Wiersma S. T. Global epidemiology of hepatitis C virus infection: new estimates of age-specific antibody to HCV seroprevalence. Hepatology 57, 1333–1342 (2013). [PubMed]
  • Omer R. E. et al. . Population-attributable risk of dietary aflatoxins and hepatitis B virus infection with respect to hepatocellular carcinoma. Nutr Cancer 48, 15–21 (2004). [PubMed]
  • Laursen L. A preventable cancer. Nature 516, S2–3 (2014). [PubMed]
  • Ferlay J. et al. . Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. International journal of cancer 136, E359–386 (2015). [PubMed]
  • Tang Z. Y. Hepatocellular carcinoma–cause, treatment and metastasis. World J Gastroenterol 7, 445–454 (2001). [PMC free article] [PubMed]
  • Song P. et al. . Screening for and surveillance of high-risk patients with HBV-related chronic liver disease: promoting the early detection of hepatocellular carcinoma in China. Biosci Trends 7, 1–6 (2013). [PubMed]
  • de Martel C., Maucort-Boulch D., Plummer M. & Franceschi S. World-wide relative contribution of hepatitis B and C viruses in hepatocellular carcinoma. Hepatology 62, 1190–1200 (2015). [PMC free article] [PubMed]
  • Llovet J. M. et al. . Hepatocellular carcinoma. Nature reviews. Disease primers 2, 16018, doi: (2016).10.1038/nrdp.2016.18 [PubMed] [Cross Ref]
  • Liu Z., Ma Y., Yang J. & Qin H. Upregulated and downregulated proteins in hepatocellular carcinoma: a systematic review of proteomic profiling studies. Omics: a journal of integrative biology 15, 61–71 (2011). [PubMed]
  • Orimo T. et al. . Proteomic profiling reveals the prognostic value of adenomatous polyposis coli-end-binding protein 1 in hepatocellular carcinoma. Hepatology 48, 1851–1863 (2008). [PubMed]
  • Sun W. et al. . Proteome analysis of hepatocellular carcinoma by two-dimensional difference gel electrophoresis: novel protein markers in hepatocellular carcinoma tissues. Molecular & cellular proteomics: MCP 6, 1798–1808 (2007). [PubMed]
  • Chen N. et al. . Quantitative proteome analysis of HCC cell lines with different metastatic potentials by SILAC. Proteomics 8, 5108–5118 (2008). [PubMed]
  • Chaerkady R. et al. . A quantitative proteomic approach for identification of potential biomarkers in hepatocellular carcinoma. Journal of proteome research 7, 4289–4298 (2008). [PMC free article] [PubMed]
  • Li C. et al. . Quantitative proteomics reveal up-regulated protein expression of the SET complex associated with hepatocellular carcinoma. Journal of proteome research 11, 871–885 (2012). [PubMed]
  • Lee S. et al. . Enhanced peptide quantification using spectral count clustering and cluster abundance. BMC Bioinformatics 12, 423; doi: (2011).10.1186/1471-2105-12-423 [PMC free article] [PubMed] [Cross Ref]
  • Megger D. A. et al. . Proteomic differences between hepatocellular carcinoma and nontumorous liver tissue investigated by a combined gel-based and label-free quantitative proteomics study. Molecular & cellular proteomics: MCP 12, 2006–2020 (2013). [PMC free article] [PubMed]
  • Gillet L. C. et al. . Targeted data extraction of the MS/MS spectra generated by data-independent acquisition: a new concept for consistent and accurate proteome analysis. Molecular & cellular proteomics: MCP 11, O111 016717, doi: (2012).10.1074/mcp.O111.016717 [PMC free article] [PubMed] [Cross Ref]
  • Haverland N. A., Fox H. S. & Ciborowski P. Quantitative proteomics by SWATH-MS reveals altered expression of nucleic acid binding and regulatory proteins in HIV-1-infected macrophages. Journal of proteome research 13, 2109–2119 (2014). [PMC free article] [PubMed]
  • Liu Y. et al. . Glycoproteomic analysis of prostate cancer tissues by SWATH mass spectrometry discovers N-acylethanolamine acid amidase and protein tyrosine kinase 7 as signatures for tumor aggressiveness. Molecular & cellular proteomics: MCP 13, 1753–1768 (2014). [PMC free article] [PubMed]
  • Sidoli S. et al. . Sequential Window Acquisition of all Theoretical Mass Spectra (SWATH) Analysis for Characterization and Quantification of Histone Post-translational Modifications. Molecular & cellular proteomics: MCP 14, 2420–2428 (2015). [PMC free article] [PubMed]
  • Xu B. et al. . Large-scale proteome quantification of hepatocellular carcinoma tissues by a three-dimensional liquid chromatography strategy integrated with sample preparation. Journal of proteome research 13, 3645–3654 (2014). [PubMed]
  • Nanjappa V. et al. . Plasma Proteome Database as a resource for proteomics research: 2014 update. Nucleic Acids Res 42, D959–965 (2014). [PMC free article] [PubMed]
  • Hou G. et al. . Biomarker Discovery and Verification of Esophageal Squamous Cell Carcinoma Using Integration of SWATH/MRM. Journal of proteome research 14, 3793–3803 (2015). [PubMed]
  • Huang da W., Sherman B. T. & Lempicki R. A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature protocols 4, 44–57 (2009). [PubMed]
  • Huang da W., Sherman B. T. & Lempicki R. A. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res 37, 1–13 (2009). [PMC free article] [PubMed]
  • Kim M. S. et al. . A draft map of the human proteome. Nature 509, 575–581 (2014). [PMC free article] [PubMed]
  • Pavlova N. N. & Thompson C. B. The Emerging Hallmarks of Cancer Metabolism. Cell Metab 23, 27–47 (2016). [PMC free article] [PubMed]
  • Kaushik A. K. et al. . Inhibition of the hexosamine biosynthetic pathway promotes castration-resistant prostate cancer. Nat Commun 7, 11612; doi: (2016).10.1038/ncomms11612 [PMC free article] [PubMed] [Cross Ref]
  • Liu B. et al. . Hepatitis B virus stimulates G6PD expression through HBx-mediated Nrf2 activation. Cell Death Dis 6, e1980; doi: (2015).10.1038/cddis.2015.322 [PMC free article] [PubMed] [Cross Ref]
  • Wang C. et al. . Identification of transaldolase as a novel serum biomarker for hepatocellular carcinoma metastasis using xenografted mouse model and clinic samples. Cancer Lett 313, 154–166 (2011). [PubMed]
  • Locasale J. W. Serine, glycine and one-carbon units: cancer metabolism in full circle. Nat Rev Cancer 13, 572–583 (2013). [PMC free article] [PubMed]
  • Vander Heiden M. G. et al. . Metabolic pathway alterations that support cell proliferation. Cold Spring Harb Symp Quant Biol 76, 325–334 (2011). [PubMed]
  • Locasale J. W. et al. . Phosphoglycerate dehydrogenase diverts glycolytic flux and contributes to oncogenesis. Nat Genet 43, 869–874 (2011). [PMC free article] [PubMed]
  • Yang M. & Vousden K. H. Serine and one-carbon metabolism in cancer. Nature Reviews Cancer 16, 650–662 (2016). [PubMed]
  • de Koning T. J. et al. . L-serine in disease and development. Biochem J 371, 653–661 (2003). [PubMed]
  • Sreekumar A. et al. . Metabolomic profiles delineate potential role for sarcosine in prostate cancer progression. Nature 457, 910–914 (2009). [PMC free article] [PubMed]
  • Kume A. et al. . The glycine cleavage system. Molecular cloning of the chicken and human glycine decarboxylase cDNAs and some characteristics involved in the deduced protein structures. J Biol Chem 266, 3323–3329 (1991). [PubMed]
  • Zhang W. C. et al. . Glycine decarboxylase activity drives non-small cell lung cancer tumor-initiating cells and tumorigenesis. Cell 148, 259–272 (2012). [PubMed]
  • Gao R. et al. . Serum Metabolomics to Identify the Liver Disease-Specific Biomarkers for the Progression of Hepatitis to Hepatocellular Carcinoma. Sci Rep 5, 18175; doi: (2015).10.1038/srep18175 [PMC free article] [PubMed] [Cross Ref]
  • Ward P. S. & Thompson C. B. Metabolic reprogramming: a cancer hallmark even warburg did not anticipate. Cancer cell 21, 297–308 (2012). [PMC free article] [PubMed]
  • Flavin R., Peluso S., Nguyen P. L. & Loda M. Fatty acid synthase as a potential therapeutic target in cancer. Future Oncol 6, 551–562 (2010). [PMC free article] [PubMed]
  • Rysman E. et al. . De novo lipogenesis protects cancer cells from free radicals and chemotherapeutics by promoting membrane lipid saturation. Cancer research 70, 8117–8126 (2010). [PubMed]
  • Carracedo A., Cantley L. C. & Pandolfi P. P. Cancer metabolism: fatty acid oxidation in the limelight. Nat Rev Cancer 13, 227–232 (2013). [PMC free article] [PubMed]
  • Qu Q., Zeng F., Liu X., Wang Q. J. & Deng F. Fatty acid oxidation and carnitine palmitoyltransferase I: emerging therapeutic targets in cancer. Cell Death Dis 7, e2226; doi: (2016).10.1038/cddis.2016.132 [PMC free article] [PubMed] [Cross Ref]
  • Warburg O. On the origin of cancer cells. Science 123, 309–314 (1956). [PubMed]
  • Wong N., Ojo D., Yan J. & Tang D. PKM2 contributes to cancer metabolism. Cancer Lett 356, 184–191 (2015). [PubMed]
  • Bluemlein K. et al. . No evidence for a shift in pyruvate kinase PKM1 to PKM2 expression during tumorigenesis. Oncotarget 2, 393–400 (2011). [PMC free article] [PubMed]
  • Montal E. D. et al. . PEPCK Coordinates the Regulation of Central Carbon Metabolism to Promote Cancer Cell Growth. Mol Cell 60, 571–583 (2015). [PMC free article] [PubMed]
  • Vincent E. E. et al. . Mitochondrial Phosphoenolpyruvate Carboxykinase Regulates Metabolic Adaptation and Enables Glucose-Independent Tumor Growth. Mol Cell 60, 195–207 (2015). [PubMed]
  • Bradford M. M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem 72, 248–254 (1976). [PubMed]
  • Wisniewski J. R., Zougman A., Nagaraj N. & Mann M. Universal sample preparation method for proteome analysis. Nature methods 6, 359–362 (2009). [PubMed]
  • Aradska J. et al. . Gel-free mass spectrometry analysis of Drosophila melanogaster heads. Proteomics 15, 3356–3360 (2015). [PubMed]
  • UniProt C. UniProt: a hub for protein information. Nucleic Acids Res 43, D204–212 (2015). [PMC free article] [PubMed]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group