Use of proteomics to identify mechanisms of hepatocellular carcinoma with the CYP2D6*10 polymorphism and identification of ANGPTL6 as a new diagnostic and prognostic biomarker
Journal of Translational Medicine volume 19, Article number: 359 (2021)
Although an association between the cytochrome P4502D6 (CYP2D6) *10 (100C>T) polymorphism and hepatocellular carcinoma (HCC) is known, the mechanism remains unclear. Here we aimed to explore mechanisms of CYP2D6*10 (100C>T) polymorphism conferring to HCC, and screen markers for HCC.
Label-free global proteome profiling with 34 normal livers and peritumor tissue from 61 HCC patients was performed, and angiopoietin-like protein-6 (ANGPTL6) was evaluated in 2 liver samples validation cohorts and 2 blood specimens validation cohorts.
We found a significantly decreased frequency of TT in HCC patients which reduced HCC susceptibility by 69.2% and was accompanied by lowered enzymatic activity for CYP2D6. Proteomic analysis revealed 1342 differentially expressed proteins (DEPs) that were associated with HCC and 88 DEPs were identified as 100 TT-related proteins, likely underlying the susceptibility to HCC. Twenty-two upregulated DEPs and 66 downregulated DEPs were mainly related to lipid metabolism and the extracellular matrix, respectively. High ANGPTL6 was associated with a higher risk to HCC and worse prognosis. ANGPTL6 was both an independent risk factor and an independent prognostic factor for HCC and exhibited strong potential for predicting HCC occurrence, with comparable AUC values and higher sensitivity compared with alpha-fetoprotein.
The TT genotype-associated decreased risk of HCC appears to be related to lowered CYP2D6 activity and altered protein expression in the tumor microenvironment, and ANGPTL6 is a promising new diagnostic and prognostic biomarker for HCC. Our findings reveal new mechanistic insights for polymorphisms related to HCC risk and provide avenues for screening for HCC.
Numerous studies have reported the important role of gene polymorphisms in the development of numerous diseases [1,2,3,4,5]. However, there is still a lack of underlying mechanisms for these polymorphisms. Indeed, phenotypic changes related to diseases caused by genotypic changes involve numerous proteins and complex processes, making it difficult to determine mechanisms. However, currently available mass-spectrometry-based proteomics, which had evolved as the preferred method for the analysis of complex proteomes, provided unprecedented opportunities to elucidate the mechanisms by which gene polymorphisms lead to cancer.
Hepatocellular carcinoma (HCC), which is characterized by high morbidity and mortality rates, is the most common type of liver cancer and accounts for 90% of cases . Hepatocarcinogenesis is well known to be associated with numerous genetic polymorphisms that are involved in its complex and multi-factorial development and contribute to malignant transformation [7,8,9,10]. To date, no effective intervention methods for HCC are based on genetics. Therefore, discovery of new genetic polymorphisms in HCC and exploration of their mechanisms might provide new approaches for identifying novel therapeutic targets and designing new or personalized treatment for HCC.
Cytochrome P450 2D6 (CYP2D6), is one of the most investigated CYPs in relation to its highly complicated polymorphisms, with 144 single nucleotide polymorphisms (SNPs) reported , resulting in remarkable differences in the activity of this enzyme [12, 13]. In recent years CYP2D6 has gained increasing attention on those polymorphisms that are related disease susceptibility, due to altered metabolism of its metabolically abundant substrates including approximately 25% of therapeutic drugs (about 223 drugs) and known and potential carcinogens such as nitrosamine, nicotine, and cotinine found in tobacco . An increasing evidence suggest that CYP2D6 polymorphisms are related to susceptibility to numerous diseases, including autoimmune diseases , Alzheimer's disease  and Parkinson's disease . It has been demonstrated that the frequency of CYP2D6-inactivating mutations is significantly lower among liver cancer patients than in healthy controls and cirrhotic subjects who did not develop liver cancer . Moreover, CYP2D6 poor metabolizer genotypes were reported significantly more frequent in healthy subjects and carriers than in hepatitis/cirrhosis and HCC patients . Although association between CYP2D6 gene polymorphisms and disease risk had been investigated in several studies, it is still unclear how CYP2D6 polymorphisms confer susceptibility to HCC.
Our previous study systematically assessed changes in 10 major CYP450 enzymes in HCC patients  and found that differences in CYP activity between physiological and pathological states were associated with increased risk for hepatofibrosis  and hepatocarcinogenesis [1, 22, 23]. CYP2D6*10 (100C>T) has been reported to be most closely associated with significantly altered metabolic activity and susceptibility to HCC, both in our and other studies [24,25,26]. However, the underlying mechanism by which this 100C>T mutation impacts hepatocarcinogenesis is not established. Thus, to address this problem, we performed label-free global proteome profiling with 34 normal livers and peritumors of 61 HCC patients to elucidate changes in the tumor microenvironment (TME) related to the CYP2D6*10 polymorphism. In addition, ANGPTL6 as a potential diagnosis and prognostic biomarker of HCC was evaluated in the above liver tissues and validated in the GEO database, another set of liver samples, and serum samples of 23 normal healthy controls and 42 HCC patients.
The study included a training cohort and four independent validation cohorts. The workflow of our study was summarized in Fig. 1. In the training cohort 105 healthy subjects and 102 HCC patients who had undergone hepatic surgery were included. Human and rat liver tissues, paraffin specimens, human serum samples and GEO database (GSE142987) constituted the validation cohorts.
Human liver samples
Normal liver samples were from hepatic haemangioma patients with normal liver function and histologic appearance. Peritumor liver specimens were obtained about 2 cm distant from the tumor tissues from HCC patients without tumor radiotherapy, chemotherapy or targeted drug therapy before surgery. All liver specimens were stored in liquid nitrogen within 30 min of resection and until use. The data were recorded via an inpatient case record, including gender, age, laboratory tests such as aspartate amino transferase (AST), alanine aminotransferase (ALT), gamma-glutamyl transferase (GGT), alkaline phosphatase (ALP), platelets (PLT), alpha‐fetoprotein (AFP), bilirubin and so on. All 207 liver tissue samples were collected at The First Affiliated Hospital of Zhengzhou University, Henan Provincial People's Hospital and Henan Provincial Tumor Hospital (Henan, China); demographic characteristics are described in our previous research . Paraffin samples of human livers included 6 normal livers, 6 peritumor tissues and tumor cores of HCC patients derived from the Second Affiliated Hospital of Zhengzhou University (Henan, China).
Serum samples were drawn at 6:00–8:00 AM after an overnight fast from a validation cohort that included 27 healthy controls and 42 HCC patients from The First Affiliated Hospital of Zhengzhou University. A single venous blood sample (3–5 mL) was taken from each participant under aseptic conditions. Serum was separated by centrifugation and stored at − 20 °C, avoiding repeated freeze–thaw cycles prior to analysis.
Genotypes of CYP2D6*10 (100C>T)
Genomic DNA was isolated from human liver tissues using a genomic DNA purification kit (QIAGEN Translational Medicine Co., Ltd, China). Genotyping of CYP2D6 100C>T of all participants was determined by Polymerase Chain Reaction (PCR) sequencing according to our previous study .
Determination of enzyme activity of CYP2D6
Enzymic activity of CYP2D6 in human liver microtomes (HLMs) of 105 healthy subjects and 102 HCC patients was measured by metabolism of dextromethorphan (National Institute for the Food and Drug Control), a probe drug for CYP2D6. HLMs were prepared by hypothermal differential centrifugation. Briefly, incubation mixtures contained HLMs (0.2 mg/mL protein), 100 mM phosphate buffer (pH 7.4), different concentrations of substrate and 1 mM NADPH. The mixture was pre-incubated for 5 min at 37 °C. Optimal incubation time was 20 min for dextromethorphan O-demethylation. For biotransformation, 0.625–960 μM dextromethorphan was added. Reactions were terminated by adding 20 μL ice-cold acetonitrile; substrate metabolite was determined by high performance liquid chromatography-fluorescence detection (HPLC-FLD). The values of the Michaelis constant (Km) and Vmax were determined by nonlinear regression analysis using GraphPad Prism 8.03 software. The intrinsic clearance (CLint) was calculated from the ratio of Vmax to Km.
Proteomics and data analysis
Proteomic analysis of 34 normal livers and peritumors of 61 HCC patients was performed by the State Key Laboratory of Proteomics, Beijing Proteome Research Center, (Beijing, China). Briefly, proteins were prepared by means of in-solution or in-gel digestion and subjected to liquid chromatography tandem mass spectrometry analysis (Q-Exactive HF LC–MS/MS). Raw data from mass spectrometry were analyzed against the human UniProt protein sequence database (version 20140922, 20,193 sequences) with MaxQuant software31 (version 220.127.116.11). The false discovery rate (FDR) for peptides was set to < 1%. The maximum mass error of the primary ion and product ion scans for the initial mass spectrometry were set to 20 ppm and 0.5 Da, respectively. The corrected error of the precursor ion was set to 5 ppm. A 6947 × 95 protein expression matrix was obtained from the 3.2.4 MaxQuant result file. Intensity-based absolute protein quantification (iBAQ) based on peak intensity was used to determine protein expression levels. The R/Bioconductor package limma v.3.24.15 34 was used to identify differentially expressed proteins (DEPs) between normal liver tissues and peritumor tissues or between CC and TT genotype groups. Proteins detected in < 50% of samples were excluded, and missing values were replaced with half the minimal value for each protein. Further details of the procedure for digesting proteins into peptides and proteomic analysis are described elsewhere .
Functional enrichment analysis
Gene ontology (GO) analysis was performed to determine the potential biological functions of ANGPTL6 and were carried out using Enrich GO function in the R package “clusterProfiler”. GO analysis was according to the threshold of P < 0.05 and q < 0.05.
Validation of Angiopoietin-like protein-6 (ANGPTL6) in liver tissues
Expression of ANGPTL6 in livers from 6 normal samples and 6 tumor-extrinsic tissues of HCC patients was determined by western blotting (WB) with primary anti-ANGPTL6 antibody (GTX118035, diluted at 1:1000). Further validation of ANGPTL6 in paraffin samples including 6 human normal livers and 6 tumor-extrinsic tissues of HCC patients was performed by immunohistochemistry (IHC) using anti-ANGPTL6 (GTX118035) antibody at 1:500 dilution. Antigen Retrieval: Trilogy™ (EDTA based, pH 8.0) buffer, 15 min.
Allograft transplantation model
Sprague Dawley (SD) rats (180–220 g), purchased from ZhengZhou University, were anesthetized with 300 mg/kg chloral hydrate by intraperitoneal injection and then placed in a supine position on the operating table. To expose the left lobe of the liver, a 1.5 cm longitudinal incision was made below the xiphoid process. Walker256 cells (1.5 × 106) in 50 μL RPMI-1640 medium (GIBCO RPMI 1640, HCC model group, n = 7) or 50 μL serum-free RPMI-1640 medium (Normal control group, n = 6) were slowly injected into the left liver lobe using a microinjector at a 30-degree angle. After allowing 3 min for the cells to disperse into the liver tissue, the skin was sutured and the surgical incision disinfected. The rats were kept in a Specific Pathogen Free (SPF) warm incubator until they had recovered from the anesthesia and then returned to animal room.
Validation of ANGPTL6 in serum
The serum level of ANGPTL6 in 23 healthy controls and 42 HCC patients was determined by ELISA. A human ANGPTL6 enzyme-linked immunosorbent assay (ELISA) Kit (Catalog No. EK19979) was used to measure ANGPTL6 levels according to the manufacturer’s instructions. The analytical (linear) measurement range was 0.312–20 ng/mL with sensitivity of 0.077 ng/mL. The intra-assay coefficient of variation (CV) and inter-assay CV of ANGPTL6 were less than 10% and 15% respectively. In addition, quality controls were included in all measurements with the results within the expected range.
Validation of ANGPTL6 in the GEO database
Serum mRNA levels of ANGPTL6 in 30 healthy controls and 35 liver cancers were obtained from the GEO database (GSE142987).
Statistical analysis was performed using SPSS 23.0 software. Graphs were plotted by R software (Version 4.0.3) and Graphpad Prism 8.3.0 software. The data were summarized as mean ± SE. Differences between two independent groups were tested with the Wilcoxon Mann–Whitney test. Receiver operating characteristics (ROC) curves were constructed to assess sensitivity, specificity, and respective areas under the curves (AUCs). A two-tailed P value < 0.05 was regarded as statistically significant.
Decreasing susceptibility to HCC is skewed in 100TT HCC patients and associated with metabolic activity of CYP2D6
A Chi-square test revealed that the study population is in Hardy–Weinberg equilibrium for the CYP2D6 gene (P > 0.05), which suggest that the genotype frequency distribution had a group representation (Additional file 1: Table S1). Significantly decreased frequency of TT was discovered in HCC patients compared with healthy subjects (Fig. 2a, 16.7% vs 46.7%, p < 0.01), as evidenced by the frequency of CC, CT and TT (6.5%, 46.5% and 16.7% respectively) in HCC patients and healthy subjects (31.4%, 21.9% and 46.7%, respectively) (Additional file 2: Table S2).
CYP2D6*10 polymorphisms were significantly associated with susceptibility to HCC, as the risk of HCC in individuals carrying the TT genotype decreased by 69.2% [(Odds Ratio (OR) = 0.308, 95% confidence interval (CI) 0.147–0.644, P = 0.002)] when compared with CC genotype carriers (Additional file 2: Table S2), suggesting that the TT genotype is a protective factor for HCC. The above findings were confirmed by the overall survival (OS) of HCC patients. A Kaplan–Meier survival curve revealed that patients with the TT genotype had longer OS than those with the CC genotype (median OS 702 vs 344 days, P = 0.365, Fig. 2b). The survival difference did not reach statistical significance, possibly due to the small number of cases. This conclusion was further substantiated by the biochemical data obtained in HCC patients. As shown in Table 1, platelets (PLT) and the albumin/globulin ratio (A/G) were higher in patients with the TT genotype (P = 0.022 and P = 0.020, respectively), while total bilirubin (TBIL), direct bilirubin (DBIL) and prothrombin time (PT) were lower when compared with individuals carrying the CC genotype (P = 0.009, P = 0.043, P = 0.043, respectively). There was no significant difference in distribution of gender and age between the two groups. Taken together, these results suggest that homozygous TT mutations may play an important role in the decreased susceptibility to HCC.
To uncover the underlying mechanism, we next assessed the effects of CYP2D6*10 polymorphisms on enzyme activity. As measured by dextromethorphan O-demethylation metabolism in human liver microsomes, the intrinsic clearance (Clint) of CYP2D6 increased significantly in HCC patients compared with normal controls (Median: 6.86 μL/min/mg vs 3.49 μL/min/mg, P < 0.01, Fig. 2c). What’s more, lower metabolic activity was observed among TT genotype subjects when compared with CC genotype individuals both in healthy samples (Median: 1.775 μL/min/mg vs 5.079 μL/min/mg, P < 0.001) and in the HCC group (Median: 1.792 μL/min/mg vs 6.228 μL/min/mg, P = 0.045, Fig. 2d). Taken together, these results suggest that a significantly reduced frequency of TT genotype was associated with decreased HCC susceptibility, at least partially, probably by decreasing the metabolic activity of CYP2D6.
Proteomic mechanisms of CYP2D6*10 polymorphisms that confer susceptibility to HCC
To investigate the mechanism by which the CYP2D6*10 mutation reduces HCC susceptibility, we performed proteomic analysis of liver tissues from 34 normal subjects and 61 HCC patients to identify HCC related proteins. A total of 1342 DEPs were identified. We then focused on 18 CC genotypes and 10 TT genotypes derived from the above HCC patients to identify DEPs by using a filter criterion of |log2-fold change (FC)|≥ 0.20 and P value < 0.05. As shown in the volcano plot (Fig. 3a) and heat-map (Fig. 3b), 88 statistically significant DEPs were identified, including 22 upregulated and 66 downregulated proteins (Additional file 3: Table S3).
The gene ontology (GO) term enrichment biological process analysis revealed that upregulated DEPs were significantly enriched in lipid transport, lipid localization, sterol metabolic process, lipid storage and cholesterol biosynthetic processes (Fig. 3c). Downregulated DEPs were mostly enriched in extracellular matrix organization, extracellular structure organization, and cell − substrate adhesion categories (Fig. 3d). Patients with the TT genotype had a lower risk of HCC which might be related to an altered tumor-extrinsic microenvironment. Subsequently, we investigated the differential expression of these CYP2D6*10-related DEPs in normal subjects and HCC patients. Surprisingly, an opposite trend further confirmed that the peritumoral microenvironment of these two genotypes was different. As shown in Fig. 3e, most of the up-regulated proteins (68.2%, 15 out of 22) in the TT genotype were lowly expressed, while the vast majority of downregulated proteins (97.0%, 64 out of 66) in this genotype were highly expressed in HCC (Additional file 4: Table S4). The significantly upregulated top nine DEPs—EFEMP1, CYP3A7, LOXL1, EFEMP2, MFI2, CHST4, CCL21, NPNT, LTBP4 and ANGPTL6, and downregulated top five DEPs—MIF2, ACE, SLC9B2, GRAMD4 and ACACB in the peritumoral microenvironment of HCC are shown in Fig. 3f and g. These results clearly demonstrate that patients with the TT genotype had lower risk of HCC that might result, at least in part, from an altered tumor-extrinsic microenvironment, perhaps by upregulating proteins related to lipid metabolism and downregulating proteins related to the extracellular matrix.
Screening biomarkers for HCC based on proteomics
A receiver operating characteristic (ROC) curve on the 100C >T-related DEPs revealed that 30 of the 88 proteins showed high ability to distinguish HCC patients from healthy subjects with a high AUC (> 0.70) (Additional file 5: Table S5). The top five proteins with the AUC > 0.80 were ANGPTL6 (AUC 0.900, sensitivity 91.8%, specificity 81.1%), LTBP1 (AUC 0.830, sensitivity 77.0%, specificity 88.2%), EFEMP2 (AUC 0.814, sensitivity 82.0%, specificity 73.5%), EFEMP1 (AUC 0.813, sensitivity 78.7%, specificity 79.4%) and KRT7 (AUC 0.804, sensitivity 77.0%, specificity 79.4%). The largest AUC was for ANGPTL6, which was comparable to serum AFP with an AUC of 0.900 (Fig. 4a and b). Interestingly, the detection rate of AFP in normal liver tissues and peritumoral tissues of HCC was 2.94% (1/34) and 18.03% (11/61), respectively, which was much lower than those of the top five proteins. What’s more, the combination of ANGPTL6 and AFP increased the diagnostic accuracy to a great extent when compared with ANGPTL6 or AFP alone, yielding an AUC of 0.963 and a sensitivity of 90.2% for HCC (Fig. 4c). Collectively, the above results show that in the parallel study the top five proteins, particularly ANGPTL6, exhibited a significantly high diagnostic performance whether alone or in combination with AFP in discriminating HCC from the healthy controls.
Among the 88 newly identified CYP2D6*10 related DEPs, we focused on ANGPTL6, a secretory protein, due to its high expression in peritumoral tissues and highest AUC for predicting HCC. In line with the proteomic results, ANGPTL6 was higher in peritumor specimens than normal livers both in human liver samples and in a rat HCC model (P < 0.01, respectively, Fig. 5a). Meanwhile, paraffin specimen analysis further confirmed these results and demonstrated almost no expression of ANGPTL6 in tumor cores (Fig. 5b).
High ANGPTL6 is a risk factor for HCC and associated with poor prognosis
To further elucidate the diagnostic potential of ANGPTL6, univariate and multivariate logistic regression analysis was performed on the training set. As shown in the forest maps (Fig. 6a), the univariate logistic regression analysis revealed that a high level of ANGPTL6 was associated with high risk for HCC [Odds ratio (OR): 3.417, 95% confidence interval (CI) 2.124–5.498; P < 0.001]. Multivariate logistic analysis showed that ANGPTL6 was an independent risk factor for HCC (OR: 2.327, 95%CI 1.191–4.547, P = 0.013, Fig. 6b). Collectively, these data strongly suggest that up-regulation of ANGPTL6 in peritumoral tissues of HCC might contribute to the pathogenesis of HCC.
Kaplan–Meier survival analysis indicated that higher ANGPTL6 in patients with HCC was associated with worse overall survival (OS), suggesting an important role of ANGPTL6 in HCC progression (Fig. 6c). Moreover, compared with the low expression group, patients with higher ANGPTL6 had lower PTL (P = 0.048), higher alanine aminotransferase (ALT, P = 0.030) and aspartate aminotransferase (AST, P = 0.032 Fig. 6d). Interestingly, older patients had a significantly higher level of ANGPTL6 than young patients (Additional file 6: Figure S1A). Serum levels of gamma-glutamyl transpeptidase (GGT) and PT showed an increasing trend in the high ANGPTL6 groups (Additional file 6: Figure S1B and S1C). What’s more, ANGPTL6 seemed to be higher in CC genotype patients than TT group (Additional file 6: Figure S1D). Univariate analysis revealed that high expression of ANGPTL6 showed a significant correlation with poor OS [hazard ratio (HR): 1.781, 95% confidence interval (CI) 1.009–3.144; P = 0.047, Fig. 6e]. Multivariate Cox analysis showed that ANGPTL6 was an independent prognostic factor for HCC (OR: 2.336, 95%CI 1.186–4.600, P = 0.014, Fig. 6f). Taken together, these data strongly suggest that high expression of ANGPTL6 is an independent prognostic indicator that correlates with poor prognosis.
Validation ANGPTL6 as diagnostic marker in serum
Serum ANGPTL6 and AFP from another validation cohort with 27 normal healthy volunteers and 42 HCC patients was determined to explore whether the serum level of ANGPTL6 can be used as an indicator for the diagnosis of HCC. There was a significant increase in serum ANGPTL6 and AFP levels in patients with HCC before surgery compared with healthy controls (Fig. 7a). Similar results were corroborated by the validation set obtained from the GEO database (GSE142987, Fig. 7b) in 30 healthy control and 35 liver cancers for mRNA expression levels. Surprisingly, the mean concentration of plasma ANGPTL6 before surgery was 3.51 (SD 0.83) ng/mL, and values dropped to 2.02 (SD 0.88) ng/mL (P < 0.001) after surgery; however, AFP was still high in patients seven days after undergoing hepatic surgery or receiving interventional therapy (Fig. 7a), suggesting that ANGPTL6 can monitor the curative effect dynamically.
We further assessed the accuracy of serum ANGPTL6 levels in the diagnosis of HCC. ROC curves showed that serum ANGPTL6 (AUC 0.826, sensitivity 92.6%, specificity 69.6%) had a comparable diagnostic value with AFP (AUC 0.807, sensitivity 72.7%, specificity 81.2%, Fig. 7c). The combination of ANGPTL6 with AFP provided better diagnostic performance of preclinical HCC before clinical diagnosis (AUC 0.852, sensitivity 72.7%, specificity 87.5%, Fig. 7c). Similar results were corroborated by the validation set obtained from the GEO database (GSE142987, AUC, 0.688 vs 0.708, Fig. 7d). These results conclusively demonstrate that serum ANGPTL6 is an indicator for the diagnosis of HCC.
In the present study, we demonstrated that CYP2D6 polymorphisms alter the susceptibility to HCC by decreasing the enzyme’s activity. Furthermore, based on proteomic analysis, we demonstrated a significant difference in the TME between CC- and TT-genotype patients that decreases the susceptibility to HCC in TT carriers, at least partly, by upregulating proteins related to lipid metabolism and downregulating proteins related to the extracellular matrix. Our current findings also reveal ANGPTL6 as a promising diagnostic and prognostic biomarker for HCC.
Relationships between genetic polymorphisms and cancer risk have been widely studied [28, 29]. However, most studies focused on polymorphisms and tumor susceptibility on the basis of epidemiologic studies, with the underlying mechanisms undiscovered. Recent advances in high-throughput proteomic technology provide expanded approaches into proteomic analyses with breadth and throughput not previously possible [30, 31]. Here, proteomic analysis was carried out to uncover the mechanisms of susceptibility to HCC underlying a CYP2D6 polymorphism. Our findings provide a unique insight into possible early intervention for HCC, especially in patients with cirrhosis.
In this study we found a significantly decreased frequency of 100TT-related genotypes with decreased susceptibility to HCC, consistent with previous studies that reported that CYP2D6 poor metabolizer (PM) genotypes or inactivating alleles were significantly more frequent in healthy subjects and carriers than in patients with hepatitis/cirrhosis and HCC [18, 19]. We speculated that the decrease in HCC susceptibility in TT genotypes may be partly explained by decreased CYP2D6 activity, due to its important role in activation of procarcinogens such as nitrosamine, nicotine, cotinine [14, 32].
Subsequently, based on proteomics, we identified a total of 1342 DEPs related to HCC, among which, 88 DEPs constitute a different TME between CC and TT and susceptibility to HCC. These findings could be explained by upregulated and downregulated DEPs that were mainly involved in lipid metabolism and the extracellular matrix, respectively, two biologic processes closely related with HCC development [33, 34]. Moreover, most of downregulated DEPs, which are involved in extracellular matrix organization, extracellular structure organization, and cell − substrate adhesion were upregulated in HCC. These findings together further validated that the TT variant was a protective factor in the occurrence of HCC and confirmed that an altered TME played a key role in polymorphisms related to susceptibility to HCC.
One of the main reasons for the high mortality rate of HCC is the lack of sensitive and reliable early diagnostic markers for screening HCC at early stage. Alpha fetoprotein (AFP) is the only currently available blood test commonly regarded as a characteristic biomarker for screening, detection and surveillance of HCC . Unfortunately, the clinical utility is limited because of low sensitivity. Therefore, it is urgent to find novel effective predictive biomarkers with early diagnostic value for HCC. ANGPTL-6, also called angiopoietin-related growth factor protein 5 (ARP5), initially known as a circulating proangiogenic factor mainly secreted from the liver [36,37,38], represents an attractive therapeutic target for anti-tumor therapy . In the current study, ANGPLT6 is dramatically increased both in peritumor tissues and serum and enriched in the extracellular matrix of HCC, suggesting an important role in the initiation and development of HCC. Upregulated ANGPTL6 expression resulted in the proliferation of undifferentiated glioblastoma cells , and knockdown of ANGPLT6 in vivo significantly inhibited AFP-producing gastric cancer by regulating endothelial cell migration and tube formation . The diagnostic accuracy of ANGPTL6 was comparable to AFP; the AUC values for the two biomarkers were 0.900 and 0.902, respectively, while the sensitivity of ANGPTL6 was markedly superior to AFP (91.8% vs 73.8%). Regrettably, the number of serum samples in our study was too small to allow stratified analyses. In particular, some AFP data were missing, possibly leading to bias in the results, such as combination of ANGTPL6 and AFP not showing significant advantages over a single biomarker. A larger sample size is needed to confirm the clinical application value of ANGPTL6.
Hematological results showed high ANGPTL6 levels were associated with lower PTL and higher ALT and AST levels, and conferred a shorter OS. Univariate and multivariate Cox proportional hazards regression analysis revealed that high expression of ANGPTL6 was an independent prognostic indicator that correlated with poor prognosis. As a consequence, our current findings highlight ANGPTL6 as a promising diagnostic and prognostic biomarker for HCC.
In summary, our findings support the mechanism by which CYP2D6*10 (100C>T) polymorphisms confer susceptibility to HCC by altering CYP2D6 activity and TME, which contributes to hepatocarcinogenesis. We further highlight ANGPTL6 as a promising diagnostic and prognostic biomarker for HCC. Our findings open new mechanistic insights for polymorphisms related HCC risk and provide avenues for clinical screening for HCC.
Availability of data materials
All the data generated and analyzed during this study are included in the manuscript and the additional materials.
Area under the curve
- Clint :
Differentially expressed proteins
Single nucleotide polymorphisms
Receiver operator characteristics curve
Khera AV, Chaffin M, Aragam KG, Haas ME, Roselli C, Choi SH, et al. Genome-wide polygenic scores for common diseases identify individuals with risk equivalent to monogenic mutations. Nat Genet. 2018;50:1219–24. https://doi.org/10.1038/s41588-018-0183-z.
Whibley C, Pharoah PD, Hollstein M. p53 polymorphisms: cancer implications. Nat Rev Cancer. 2009;9:95–107. https://doi.org/10.1038/nrc2584.
Schmiedel BJ, Singh D, Madrigal A, Valdovino-Gonzalez AG, White BM, Zapardiel-Gonzalo J, et al. Impact of genetic polymorphisms on human immune cell gene expression. Cell. 2018;175:1701-15.e16. https://doi.org/10.1016/j.cell.2018.10.022.
Aleksandrova E, Vlaykova T, Ananiev J, Gulubova M. Association of IL-12Bpro polymorphism with tumor-infiltrating dendritic cells in colorectal cancer. J BUON. 2017;22:888–93.
Dopazo J, Amadoz A, Bleda M, Garcia-Alonso L, Alemán A, García-García F, et al. 267 Spanish exomes reveal population-specific differences in disease-related genetic variation. Mol Biol Evol. 2016;33:1205–18. https://doi.org/10.1093/molbev/msw005.
Laursen L. A preventable cancer. Nature. 2014;516:S2-3. https://doi.org/10.1038/516S2a.
Nahon P, Zucman-Rossi J. Single nucleotide polymorphisms and risk of hepatocellular carcinoma in cirrhosis. J Hepatol. 2012;57:663–74. https://doi.org/10.1016/j.jhep.2012.02.035.
Thomas DL, Thio CL, Martin MP, Qi Y, Ge D et al. Genetic variation in IL28B and spontaneous clearance of hepatitis C virus. Nature. 2009;461:798–801. doi: https://doi.org/10.1038/nature08463
Bressac B, Kew M, Wands J, Ozturk M. Selective G to T mutations of p53 gene in hepatocellular carcinoma from southern Africa. Nature. 1991;350:429–31. https://doi.org/10.1038/350429a0.
Niu ZS, Niu XJ, Wang WH. Genetic alterations in hepatocellular carcinoma: an update. World J Gastroenterol. 2016;22:9069–95. https://doi.org/10.3748/wjg.v22.i41.9069.
Preissner SC, Hoffmann MF, Preissner R, Dunkel M, Gewiess A, Preissner S. Polymorphic cytochrome P450 enzymes (CYPs) and their role in personalized therapy. PLoS ONE. 2013;8: e82562. https://doi.org/10.1371/journal.pone.0082562.
Neafsey P, Ginsberg G, Hattis D, Sonawane B. Genetic polymorphism in cytochrome P450 2D6 (CYP2D6): Population distribution of CYP2D6 activity. J Toxicol Environ Health B Crit Rev. 2009;12:334–61. https://doi.org/10.1080/10937400903158342.
Zanger UM, Schwab M. Cytochrome P450 enzymes in drug metabolism: regulation of gene expression, enzyme activities, and impact of genetic variation. Pharmacol Ther. 2013;138:103–41. https://doi.org/10.1016/j.pharmthera.2012.12.007.
Zhou SF. Polymorphism of human cytochrome P450 2D6 and its clinical significance: Part I. Clin Pharmacokinet. 2009;48:689–723. https://doi.org/10.2165/11318030-000000000-00000.
Lee YH, Bae SC. Association between Functional CYP2D6 polymorphisms and susceptibility to autoimmune diseases: a meta-analysis. Immunol Invest. 2017;46:109–22. https://doi.org/10.1080/08820139.2016.1226898.
Lu Y, Qin X, Li S, Zhang X, He Y, Peng Q, et al. Quantitative assessment of CYP2D6 polymorphisms and risk of Alzheimer’s disease: a meta-analysis. J Neurol Sci. 2014;343:15–22. https://doi.org/10.1016/j.jns.2014.05.033.
Lu Y, Peng Q, Zeng Z, Wang J, Deng Y, Xie L, et al. CYP2D6 phenotypes and Parkinson’s disease risk: a meta-analysis. J Neurol Sci. 2014;336:161–8. https://doi.org/10.1016/j.jns.2013.10.030.
Agúndez JA, Ledesma MC, Benítez J, Ladero JM, Rodríguez-Lescure A, Díaz-Rubio E, et al. CYP2D6 genes and risk of liver cancer. Lancet. 1995;345:830–1. https://doi.org/10.1016/s0140-6736(95)92965-7.
Silvestri L, Sonzogni L, De Silvestri A, Gritti C, Foti L, Zavaglia C, et al. CYP enzyme polymorphisms and susceptibility to HCV-related chronic liver disease and liver cancer. Int J Cancer. 2003;104:310–7. https://doi.org/10.1002/ijc.10937.
Zhang HF, Wang HH, Gao N, Wei JY, Tian X, Zhao Y, et al. Physiological content and intrinsic activities of 10 cytochrome P450 isoforms in human normal liver microsomes. J Pharmacol Exp Ther. 2016;358:83–93. https://doi.org/10.1124/jpet.116.233635.
Thomson S. A system for quantitative studies on interactions between Friend leukemia virus and hemopoietic cells. Proc Soc Exp Biol Med. 1969;130:227–32. https://doi.org/10.3181/00379727-130-33527.
Gao J, Wang Z, Wang GJ, Gao N, Li J, Zhang YF, et al. From hepatofibrosis to hepatocarcinogenesis: Higher cytochrome P450 2E1 activity is a potential risk factor. Mol Carcinog. 2018;57:1371–82. https://doi.org/10.1002/mc.22851.
Gao J, Wang Z, Wang GJ, Zhang HX, Gao N, Wang J, et al. Higher CYP2E1 activity correlates with hepatocarcinogenesis induced by diethylnitrosamine. J Pharmacol Exp Ther. 2018;365:398–407. https://doi.org/10.1124/jpet.117.245555.
Zhou J, Wen Q, Li SF, Zhang YF, Gao N, Tian X, et al. Significant change of cytochrome P450s activities in patients with hepatocellular carcinoma. Oncotarget. 2016;7:50612–23. https://doi.org/10.18632/oncotarget.9437.
Gao J, Tian X, Zhou J, Cui MZ, Zhang HF, Gao N, et al. From genotype to phenotype: cytochrome P450 2D6-mediated drug clearance in humans. Mol Pharm. 2017;14:649–57. https://doi.org/10.1021/acs.molpharmaceut.6b00920.
Mochizuki J, Murakami S, Sanjo A, Takagi I, Akizuki S, Ohnishi A. Genetic polymorphisms of cytochrome P450 in patients with hepatitis C virus-associated hepatocellular carcinoma. J Gastroenterol Hepatol. 2005;20:1191–7. https://doi.org/10.1111/j.1440-1746.2005.03808.x.
Jiang Y, Sun A, Zhao Y, Ying W, Sun H, Yang X, et al. Proteomics identifies new therapeutic targets of early-stage hepatocellular carcinoma. Nature. 2019;567:257–61. https://doi.org/10.1038/s41586-019-0987-8.
Turajlic S, Sottoriva A, Graham T, Swanton C. Resolving genetic heterogeneity in cancer. Nat Rev Genet. 2019;20:404–16. https://doi.org/10.1038/s41576-019-0114-6.
Bao Y, Prescott J, Yuan C, Zhang M, Kraft P, Babic A, et al. Leucocyte telomere length, genetic variants at the TERT gene region and risk of pancreatic cancer. Gut. 2017;66:1116–22. https://doi.org/10.1136/gutjnl-2016-312510.
Zhang Z, Wu S, Stenoien DL, Paša-Tolić L. High-throughput proteomics. Annu Rev Anal Chem. 2014;7:427–54. https://doi.org/10.1146/annurev-anchem-071213-020216.
Karczewski KJ, Snyder MP. Integrative omics for health and disease. Nat Rev Genet. 2018;19:299–310. https://doi.org/10.1038/nrg.2018.4.
Nebert DW, McKinnon RA, Puga A. Human drug-metabolizing enzyme polymorphisms: effects on risk of toxicity and cancer. DNA Cell Biol. 1996;15:273–80. https://doi.org/10.1089/dna.1996.15.273.
Ma C, Kesarwala AH, Eggert T, Medina-Echeverz J, Kleiner DE, Jin P, et al. NAFLD causes selective CD4(+) T lymphocyte loss and promotes hepatocarcinogenesis. Nature. 2016;531:253–7. https://doi.org/10.1038/nature16969.
Affo S, Yu LX, Schwabe RF. The role of cancer-associated fibroblasts and fibrosis in liver cancer. Annu Rev Pathol. 2017;12:153–86. https://doi.org/10.1146/annurev-pathol-052016-100322.
Chen T, Dai X, Dai J, Ding C, Zhang Z, Lin Z, et al. AFP promotes HCC progression by suppressing the HuR-mediated Fas/FADD apoptotic pathway. Cell Death Dis. 2020;11:822. https://doi.org/10.1038/s41419-020-03030-7.
Kim I, Kim HG, Kim H, Kim HH, Park SK, Uhm CS, et al. Hepatic expression, synthesis and secretion of a novel fibrinogen/angiopoietin-related protein that prevents endothelial-cell apoptosis. Biochem J. 2000;346(Pt 3):603–10.
Oike Y, Ito Y, Maekawa H, Morisada T, Kubota Y, Akao M, et al. Angiopoietin-related growth factor (AGF) promotes angiogenesis. Blood. 2004;103:3760–5. https://doi.org/10.1182/blood-2003-04-1272.
Urano T, Ito Y, Akao M, Sawa T, Miyata K, Tabata M, et al. Angiopoietin-related growth factor enhances blood flow via activation of the ERK1/2-eNOS-NO pathway in a mouse hind-limb ischemia model. Arterioscler Thromb Vasc Biol. 2008;28:827–34. https://doi.org/10.1161/ATVBAHA.107.149674.
Chen E, Tang C, Peng K, Cheng X, Wei Y, Liu T. ANGPTL6-mediated angiogenesis promotes alpha fetoprotein-producing gastric cancer progression. Pathol Res Pract. 2019;215: 152454. https://doi.org/10.1016/j.prp.2019.152454.
Cui JG, Zhao Y, Sethi P, Li YY, Mahta A, Culicchia F, et al. Micro-RNA-128 (miRNA-128) down-regulation in glioblastoma targets ARP5 (ANGPTL6), Bmi-1 and E2F–3a, key regulators of brain cell proliferation. J Neurooncol. 2010;98:297–304. https://doi.org/10.1007/s11060-009-0077-0.
We thank the department of liver and gallbladder surgery, Pathology and Laboratory of The First and Second Affiliated Hospital of Zhengzhou University, Henan Provincial People's Hospital and Henan Provincial Tumor Hospital (Henan, China).
This work was supported by National Natural Science Foundation of China (NSFC) (Nos. 81473279, 81673507, 81872931 and 82073930); Zhengzhou Major Scientific and Technological Innovation Projects (No. 2020CXZX0076).
Ethics approval and consent to participate
This study was carried out abiding by the principles established by the Declaration of Helsinki and approved by the Ethical Committee of the Zhengzhou University (Henan, China). All participants were provided with an informed consent that they signed prior to the experiments.
Consent for publication
The authors declare no potential conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Hardy–Weinberg Equilibrium (HWE).
The genotype frequency of CYP2D6*10 (100 C>T) in healthy volunteers (n = 105) and HCC patients (n = 96).
Compared with CC groups, the 22 upregulated and 66 downregulated DEPs in TT groups based on proteomics analysis.
Compared with healthy subjects, the 17 upregulated and 71 downregulated DEPs in HCC groups based on proteomics analysis.
DEPs with the area under the ROC curve (AUC) > 7.0.
ANGPTL6 has clinical significance. a High ANGPTL6 expression was more common in the elderly group. Serum levels of GGT (b) and PT (c) show an increasing trend in the high ANGPTL6 groups. d ANGPTL6 in TT and CC group. HCC, hepatocellular carcinoma. GGT, gamma-glutamyl transpeptidase. PT, prothrombin time.
About this article
Cite this article
Hu, G., Gao, F., Wang, G. et al. Use of proteomics to identify mechanisms of hepatocellular carcinoma with the CYP2D6*10 polymorphism and identification of ANGPTL6 as a new diagnostic and prognostic biomarker. J Transl Med 19, 359 (2021). https://doi.org/10.1186/s12967-021-03038-3
- Hepatocellular carcinoma