- Open Access
Elucidating the role of missense SNP of protein kinase C epsilon in HCV-induced hepatocellular carcinoma
BMC Cancer volume 23, Article number: 147 (2023)
The protein kinase C (PKC) family of serine/threonine kinases contains more than ten isozymes that are involved in multiple signaling pathways, including cell cycle regulation and carcinogenesis. The PKCε isozyme is an oncogene known to be upregulated in various signaling pathways involved in hepatitis C virus (HCV)-induced hepatocellular carcinoma (HCC). However, there is no known association of missense SNPs in PKCε with this disease, which can be a potential biomarker for early diagnosis and treatment. This research reveals a novel missense SNP in PKCε that is associated with HCV-induced HCC in the Pakistani population.
The PKCε SNP with amino acid substitution of E14K was chosen for wet lab analysis. Tetra ARMS-PCR was employed for the identification of high-risk SNP in PKCε of HCV-induced HCC patients. Liver function testing was also performed for comparison between the liver condition of the HCC patient and control group, and the viral load of HCC patient samples was evaluated to determine any alteration in the viral infectivity between different genotypes of the selected high-risk PKCε variant SNP.
Frequency distribution of the homozygous GG genotype was found to be highest among HCV-induced HCC patients and was also found to be significantly associated with disease development and progression. The p values of comparative data obtained for the other two genotypes, heterozygous AG and homozygous AA, of the SNP also showed the significance of the data for these alleles. Still, their odds ratio and relative risk analysis did not indicate their association with HCV-induced HCC.
The distribution of a genotype GG of PKCε has been found in HCV- induced HCC patients. Therefore, these PKCε SNP have the potential to be biomarkers for HCV-induced HCC. Further investigation using a larger sample size would provide additional insight into these initial data and open a new avenue for a better prognosis of this disease.
Protein kinase C (PKC) is a large family of serine-threonine kinases which are expressed in various types of cells. More than ten types of PKC isozymes have been discovered so far which have been categorized into three sub-classes, the classical cPKCs (α, β1, β2, γ), novel nPKCs (δ, ε, η), and atypical aPKCs (λ, ζ) . PKCε is a calcium independent and phorbol ester /Diacyl glycerol (DAG) dependent isoform in the novel PKC subclass. It is activated and compartmentalized to specific cellular sites depending on the secondary messengers and the adaptor proteins activated in response to the extracellular stimuli [2, 3]. Various receptors like tyrosine kinase receptors, tyrosine kinase-coupled receptors and G-protein-coupled receptors can lead to PKCε induction , which in turn is involved in various physiological processes such as immune system activation, inflammatory responses, exocrine, endocrine and nervous system stimulation [5, 6]. PKCε has been positively associated with several kinds of cancers and plays major roles in several tumorigenic pathways [7,8,9,10,11]. Various studies have shown that PKCε directly phosphorylates AKT at S473  and SMAD3 at S213 , activating their respective pathways resulting in uncontrolled cell survival and proliferation. Evidence indicates the role of PKCε in modulating the expression of AFP, which has significance as a diagnostic and follow-up marker of HCC . PKCε has also shown to facilitate hepatic steatosis and impairment of hepatic insulin signaling, thus, promoting insulin resistance which is a major factor in liver disease and development of HCC in patients with chronic HCV infection [15,16,17]. Besides that, PKCε activates ATF2 signaling which is involved in CAP2 expression that is induced in conditions of ER stress and promotes epithelial to mesenchymal transformation (EMT), an important hallmark of cancer, via Rac1 and ERK in HCC .
Most of the data available so far have only indicated the carcinogenic role of PKCε at the functional level. However, genotype association studies for other PKC isoforms with cancer susceptibility have been conducted, indicating the pathogenic effect of certain SNPs [19, 20]. Recently, the pathogenic SNPs in PKCε were also predicted through bioinformatics tools that indicated its missence SNPs pathogenicity and potential impact on the structure and function of PKCε protein and non-coding variants impact on the gene expression of the protein [21, 22].
HCC is the most frequently occurring type of liver cancer accounting for ~ 85% of primary liver cancers. The risk factors associated with it include environmental exposures and hepatitis B and C viruses [23, 24]. Nearly 25% of HCC cases worldwide can be attributed to chronic HCV infection in patients. Chronic HCV infection results in development of inflammation-induced lesions in the liver which causes steatohepatitis leading to fibrosis. Over the span of many years, these events progress to cirrhosis or HCC in some patients [25, 26]. Studies at molecular levels have shown the involvement of PKCs in the progression of HCC , and research demonstrates that PKCε is engaged in various mechanisms leading to the progression of HCC [28, 29]. However, there are no known studies investigating the association of PRKCE SNPs with risk of HCC.
Based on evidence that nsSNPs in other PKC isoforms are correlated with the risk of cancers, we examined the association of a non-synonymous SNP (snSNP) in PRKCE, i.e. rs1553369874, in HCV-induced HCC in Pakistani population. Previously, nsSNP rs1553369874 was predicted as deleterious through several bioinformatics tools . The study also indicated this SNP impact on the PKCε stability, flexibility, structure, and function. Therefore, this nsSNP was selected to probe its association with HCV-induced HCC.
Project approval was received from the Institution Review Board (IRB) of ASAB, NUST. Only patients with HCV-induced HCC were included in this research and written consent was taken from them prior to sample acquisition. HCV negative HCC cases were ignored as they were not the part of the hypothesis of this study. The control group taken was of healthy people without any known serious health or liver condition. Furthermore, a total of 63 patients belong to stage 1&2 and 37 from stage 3&4. None of the patients received disease treatment. Blood samples were acquired from patients of both genders whose ages were between 20 and 60 years. One hundred samples from HCC patients and an equal number of healthy control samples were collected for carrying out the procedures for identification and association of selected PKCε SNP with HCV-induced HCC.
Genomic DNA was extracted from blood using a phenol-chloroform extraction method . This two day procedure utilizes four types of buffer solutions (namely A,B,C, and D), SDS, and proteinase K. Solution A was composed of 0.32 M sucrose, 5mM magnesium chloride, Triton X-100 and Tris (pH 7.5). Solution B consisted of 400mM sodium chloride, 2mM EDTA, and Tris (pH 7.5) while solution C solely consisted of phenol. Solution D was a mixture of isoamylalcohol and chloroform. 750 µl of blood was processed with solution A, B, 20% sodium dodecyl sulfate (SDS) and proteinase K on the first day followed by overnight incubation at 37 °C. On the second day the processing was completed using solution C, D, and DNA pellet was precipitated, washed, and dried using sodium acetate, chilled iso-propanol and ethanol in multiple steps. DNA was stored at 4 °C. The presence of DNA was confirmed with gel electrophoresis through 1% agarose at 120 volts. The quantification of DNA was performed on NanoDrop™ 8000 Spectrophotometer (Thermofisher scientific). Nanodrop value near 1.8 for 280/260 ratio was considered pure for extracted genomic content.
For the analysis of PKCε gene using tetra ARMS-PCR, four primers were designed using primer 1 software. Forward inner (A allele) CTTCTTAAGATCAAAATCTTCA, forward outer CACAAGGTGTAGGGAGTGT, reverse inner (G allele) CTTCAAGCTCACGGCATC, and reverse outer GCTGTTGGTCTTCTGCTT. The tetra ARMS PCR reaction mixture contained 2 µl PCR buffer, 4.5 µl MgCl2, 2 µl dNTPs, 0.4 µl Taq Polymerase, 1 µl of each of the four primers, 1.5 µl the DNA sample and 10.6 µl PCR water. The PCR reaction entailed initial denaturation for 5 min at 92 °C, 35 cycles for 30 s each with denaturation 95 °C, annealing 50 °C, and extension at 72 °C along with the final extension at 72 °C for 7 min. PCR products were visualized by 2% (w/v) agarose gel electrophoresis.
Measurement of ALT levels
Alanine aminotransferase (ALT) levels of both control and disease samples were measured to detect the change in liver functions of the HCC patients. An ALT assay kit by Merck (Darmstadt,
Germany) was used for this purpose. Blood serum from each of the samples was taken for ALT measurement. The chemical reaction between serum and kit reagents produced lactate as the end product and its absorbance was measured at 340 nm using a (UV-visible) spectrophotometer (Microlab 300 – Semi-automated chemistry analyzer). The concentration of ALT was calculated according to the formula, where kit factor is 1746 for any absorbance using 340 nm wavelength at room temperature:
Normal ALT values lie below 40–41 IUL− 1.
Viral RNA extraction
The viral load was measured in HCV-induced HCC blood samples for comparison and association of SNPs with the severity of disease. For that, RNA of hepatitis C virus was extracted from the blood serum of the patients using Favorgen kit. The extracted RNA was dissolved in RNase free water and stored at -20 °C for qRT-PCR.
Quantitative reverse transcriptase polymerase chain reaction (qRT-PCR)
Extraction of viral genome was performed from the blood plasma of patients through Instant virus RNA kit CE IVD ªAnalytik Jena and qRT-PCR employing TaqMan probes (TaqMan™ Fast Virus 1-Step Master Mix) was used for quantifying the viral RNA in real time. PCR reaction mix concentration and protocol for reaction was used as per the manufacturer’s instructions. Two standards with known concentrations were used to determine the upper and lower limit of the concentrations obtained from the samples. The CT (cyclic-threshold) value was determined for the samples; those samples with lower CT value contain more viral load than samples with higher CT value.
GraphPad prism 9.0 (GraphPad, California, USA) and Microsoft Excel (Microsoft cooperation 2016) were employed for statistical analysis and making graphs. Frequency distribution of genotypes, odds ratio, relative risk along with the genotypes association with the disease was found out through Fisher Exact test. P-values below 0.05 were taken as significant.
The data from tetra ARMS PCR was collected and statistical tools were applied for exploring its association with HCV induced HCC, which is presented in Tables 1 and 2. The comparitive genotype analysis between the control and patient sample data of all alleles of the selected locus of PKCε gene showed significance (P < 0.05) of the acquied data for all 3 genotypes of that locus. The homozygous wild genotype GG demonstrated high association with HCC with an odds ratio 4.19 and relative risk 2.1 (P < 0.0001). The other two genotypes were calculated to have odds ratio and relative risk less than 1, suggesting they are not associated with HCC development and progression (Table 1).
Gender based comparison of the genotype data of control and patient samples demonstrated significant difference in the data of all three genotypes in males and GG genotype in females (P < 0.05). Concerning the GG genotype, the data depicted a slight difference in males and females, with the odds ratio and relative risk of GG association with HCC higher in males than in females. In males the odds ratio of GG is 6.824 and relative risk is 2.737, whereas in females the odds ratio is 2.794 and relative risk is 1.711. The P-value of other two genotypes AG in both males and females, and AA in males show their protective significance in HCV-induced HCC (Table 2). The results of genotype data obtained from tetra ARMS PCR was also analyzed with respect to the stages of HCC. The result of this analysis shows significance of the data for all genotypes in all stages of HCC except for genotype AA and AG in stages 3 and 4. The data shows that there is a high odds ratio and relative risk of GG genotype in all 4 stages of HCC (Table 3).
ALT levels analysis
ALT level test is an auxillary test that is performed to check the deviation in the normal liver function. It was performed for only functional analysis of the liver and has no correlation with the PKCε gene expression. The mean levels of ALT in control versus HCV-induced HCC patients is illustrated in Fig. 1. The ALT values of the control group (37.5 ± 10.2) compared to those of HCV-induced HCC patients (106.84 ± 52.5) indicate a considerable difference between the groups.
Viral load analysis
The viral load measured by qRT-PCR was plotted in the graph with respect to the three genotypes of the selected locus of PKCε as shown in Fig. 2. This graph shows clear association of the homozygous wild GG genotype with high abundance of viral load in HCV patients.
The presence of high-risk SNPs in PKCε represent a potential cause for its oncogenic properties observed in several types of cancers. From the mechanistic point of view, the role of PKCε in endothelial to mesenchymal transformation (EMT) may contribute to the invasive, migratory and metastatic properties of the cancer cells . Its progressively elevated levels have been observed in the hyperplastic cells of rat liver in comparison to the normal liver cells which might have been associated with their malignant development . However, the underlying cause can be a genetic polymorphism such as that present in several other genes which have been proved to be associated with development of disease and its progression. The presence of two novel SNPs in PKCε has been previously observed to be associated with the risk of occurrence of prostate cancer in the Han population . Similarly, SNPs in five different genes involved in nucleotide excision repair pathways has been associated with the risk of development of HCC . Therefore, we investigated a potentially damaging missense SNP of PKCε and examined the frequency of variants in HCV-induced HCC samples versus control samples.
Previously, a thorough in silico study to analyze the genetic variations of PKCε and their effects on protein structure, folding and function was carried out. Numerous in silico tools and software were employed for this purpose. The study aimed to identify potential disease-causing variants of PKCε gene for which the damaging non-synonymous variants were focused. The structure prediction of the protein with those variants along with further analysis exhibited the impact of variation on the dynamic behavior of protein. The delineation of phenotypic impact on PKCε due to the presence of variations helped in selecting a high-risk SNP that can further impact the functional properties of the protein .
The high-risk SNP was selected after screening through well-known and frequently used tools for this purpose [34, 35]. After careful screening the potential high-risk SNP of PKCε rs1553369874 was chosen for exploration of its association with HCV-induced HCC.
To see the damage inflicted on the liver of patients in comparison to the healthy individual, the ALT level of patient and control samples were measured. The results of this test provided a clear indication about the extent of liver damage in patients. Previous evidence suggests that ALT levels increase in patients with chronic liver damage and are significantly associated with the chronic HCV infection leading to HCC development .
Apart from that, elevated viral load with the GG genotype in HCV-induced HCC patients was found. This could be an indication of the increased susceptibility of idividuals with this genotype to hepatitis C viral pathogenicity. Previously, SNPs have been associated with increased viral load and disease progression in individuals with chronic viral infections [37, 38].
The statistical analysis of data obtained from ARMS-PCR after amplification of the selected SNP (rs1553369874) in 100 controls and 100 HCV-induced HCC patient’s samples pointed out that the homozygous wild genotype GG has significant association with HCV-induced HCC with 2.1 relative risk and 4.19 odds ratio. The analysis of this data with respect to gender also indicated the association of the same genotype in both males and females, although its relative risk in males and odds ratio is slightly greater than in females. This difference could also be due to the relatively greater incidence rate of liver cancer, more specifically HCC in males than in females. The male to female ratio of incidence and mortality of HCC exceeds 2.5 for both [39, 40].
The damage sustained by the livers of HCV-induced HCC patients was confirmed though the measurement of ALT levels in them. The results from tetra ARMS PCR suggest the wild-type GG variant to be significantly associated with the disease progression. Viral load of this GG genotype was also very high among all three genotypes of HCV-induced HCC patients, indicating the increased susceptibility to disease in patients with this genotype. Thus, this SNP can be further explored for a potential biomarker for HCC prognosis by in vitro and in vivo studies. Investigation of a larger sample size along with the transcriptome and expression analysis of PKCε containing this SNP will provide better insight into its role in HCV-induced HCC occurrence.
All data generated or analyzed during this study are included in this article.
Schaefer M. Protein Kinase C 2007.
Webb BL, Hirst SJ, Giembycz MA. Protein kinase C isoenzymes: a review of their structure, regulation and role in regulating airways smooth muscle tone and mitogenesis. Br J Pharmacol. 2000;130(7):1433–52.
Newton AC. Protein kinase C: perfectly balanced. Crit Rev Biochem Mol Biol. 2018;53(2):208–30.
Akita Y. Protein kinase C-ε (PKC-ε): its unique structure and function. J Biochem. 2002;132(6):847–52.
Aksoy E, Goldman M, Willems F. Protein kinase C epsilon: a new target to control inflammation and immune-mediated disorders. Int J Biochem Cell Biol. 2004;36(2):183–8.
Obis T, et al. The novel protein kinase C epsilon isoform modulates acetylcholine release in the rat neuromuscular junction. Mol Brain. 2015;8(1):1–16.
Kang J-H. Protein kinase C (PKC) isozymes and cancer New Journal of Science, 2014. 2014.
Jain K, Basu A. The multifunctional protein kinase C-ε in cancer development and progression. Cancers. 2014;6(2):860–78.
Gorin MA, Pan Q. Protein kinase Cε: an oncogene and emerging tumor biomarker. Mol Cancer. 2009;8(1):1–8.
Khan K et al. Unravelling Structure, Localization, and Genetic Crosstalk of KLF3 in Human Breast Cancer BioMed research international, 2020. 2020.
Safi S, et al. Predicting 3D structure, Cross Talks, and Prognostic significance of KLF9 in Cervical Cancer. Front Oncol. 2021;11:797007–7.
Zhang J, et al. Functional proteomic analysis of a three-tier PKCε-Akt-eNOS signaling module in cardiac protection. Am J Physiol Heart Circ Physiol. 2005;288(2):H954–61.
Xu W, et al. Crosstalk of protein kinase C ε with Smad2/3 promotes tumor cell proliferation in prostate cancer cells by enhancing aerobic glycolysis. Cell Mol Life Sci. 2018;75(24):4583–98.
Doumit R, et al. The role of protein kinase C isoforms in telomerase activity and alphafetoprotein secretion by the hepatocellular carcinoma cell lines. Cancer Res. 2016;76(14Supplement):940–0.
Liu LF, Purushotham A, Belury M. Activation of peroxisome proliferator-activated receptor (PPAR) α by conjugated linoleic acid is modulated by protein kinase C ε. Wiley Online Library; 2006.
Hung C-H, et al. Insulin resistance is associated with hepatocellular carcinoma in chronic hepatitis C infection. World J gastroenterology: WJG. 2010;16(18):2265.
Mu W, et al. Potential nexus of non-alcoholic fatty liver disease and type 2 diabetes mellitus: insulin resistance between hepatic and peripheral tissues. Front Pharmacol. 2019;9:1566.
Yoon S, Shin B, Woo HG. Endoplasmic reticulum stress induces CAP2 expression promoting epithelial-mesenchymal transition in Liver Cancer cells. Mol Cells. 2021;44(8):569.
Parzefall T et al. PRKCA overexpression is frequent in young oral tongue squamous cell carcinoma patients and is associated with poor prognosis Cancers, 2021;13(9):p. 2082.
Liu X, et al. Genetic variants of the peroxisome proliferator-activated receptor (PPAR) signaling pathway genes and risk of pancreatic cancer. Mol Carcinog. 2020;59(8):930–9.
Khan K, et al. Influence of PRKCE non-synonymous variants on protein dynamics and functionality. Human Molecular Genetics; 2022.
Khan K, et al. PRKCE non-coding variants influence on transcription as well as translation of its gene. RNA Biol. 2022;19(1):1115–29.
Srivatanakul P, Sriplung H, Deerasamee S. Epidemiology of liver cancer: an overview. Asian Pac J Cancer Prev. 2004;5(2):118–25.
Kim E, Viatour P. Hepatocellular carcinoma: old friends and new tricks. Exp Mol Med. 2020;52(12):1898–907.
Vescovo T, et al. Molecular mechanisms of hepatitis C virus–induced hepatocellular carcinoma. Clin Microbiol Infect. 2016;22(10):853–61.
Khatun M, Ray R, Ray RB. Hepatitis C virus associated hepatocellular carcinoma. Adv Cancer Res. 2021;149:103–42.
Lu H-C, et al. Expression of protein kinase C family in human hepatocellular carcinoma. Pathol Oncol Res. 2010;16(3):385–91.
Yan R, Niu C-Y, Tian Y. Roles of Autophagy and protein kinase C-epsilon in lipid metabolism of nonalcoholic fatty liver cell models. Arch Med Res. 2018;49(6):381–90.
Pascut D, et al. HCV Proteins modulate the host cell miRNA expression contributing to Hepatitis C pathogenesis and Hepatocellular Carcinoma Development. Cancers. 2021;13(10):2485.
Ghatak S, Muthukumaran RB, Nachimuthu SK. A simple method of genomic DNA extraction from human samples for PCR-RFLP analysis. J Biomol Tech. 2013;24(4):224–31.
Perletti G, et al. ϵPKC Acts like a marker of Progressive Malignancy in Rat Liver, but fails to enhance tumorigenesis in rat hepatoma cells in culture. Biochem Biophys Res Commun. 1996;221(3):688–91.
Li Q, et al. Two novel PRKCI polymorphisms and prostate cancer risk in an Eastern Chinese Han population. Mol Carcinog. 2015;54(8):632–41.
Wang B, et al. The association of six polymorphisms of five genes involved in three steps of nucleotide excision repair pathways with hepatocellular cancer risk. Oncotarget. 2016;7(15):20357.
Kumar CV et al. Computational analysis reveals the association of threonine 118 methionine mutation in PMP22 resulting in CMT-1A Advances in bioinformatics, 2014. 2014.
Suybeng V, et al. Comparison of pathogenicity prediction tools on somatic variants. J Mol Diagn. 2020;22(12):1383–92.
Hayashi J, et al. A relationship between the evolution of hepatitis C virus variants, liver damage, and hepatocellular carcinoma in patients with hepatitis C viremia. J Infect Dis. 2000;181(5):1523–7.
Anokhin VV et al. Previously unidentified single nucleotide polymorphisms in HIV/AIDS cases associate with clinical parameters and disease progression BioMed research international, 2016. 2016.
Zhong Y-W, et al. IL28B SNP rs12979860 is the critical predictor for sustained viral response in chinese children aged 1 to 6 years with chronic hepatitis C. Int J Biol Sci. 2016;12(11):1357.
Singal AG, Lampertico P, Nahon P. Epidemiology and surveillance for hepatocellular carcinoma: new trends. J Hepatol. 2020;72(2):250–61.
Hafeez S, et al. Trends in Cancer Prevalence in Punjab, Pakistan: a systematic study from 2010 to 2016. J Bioresource Manage. 2020;7(2):8.
The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project no. (IFKSURG- 1582).
.The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project no. (IFKSURG- 1582).
Ethics approval and consent to participate
The experimental protocol for the use of human blood samples was approved (Ref: No: IRB No. 04-2019-03/06) by the ethical committee of Combined Military Hospital and ASAB, NUST. Informed consent was obtained from all subjects and/or their legal guardian(s). All the steps/ methods were performed in accordance with the relevant guidelines and regulations (Declaration of Helsinki).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Rehman, A., Shabbir, M., Badshah, Y. et al. Elucidating the role of missense SNP of protein kinase C epsilon in HCV-induced hepatocellular carcinoma. BMC Cancer 23, 147 (2023). https://doi.org/10.1186/s12885-023-10618-7
- Protein kinase C (PKC)
- Hepatocellular carcinoma (HCC)
- Single nucleotide polymorphism (SNP)
- Hepatitis C virus (HCV)