Skip to main content

The pattern of gene copy number alteration (CNAs) in hepatocellular carcinoma: an in silico analysis

Abstract

Background

Hepatocellular carcinoma (HCC) is the most common type of liver cancer that occurs predominantly in patients with previous liver conditions. In the absence of an ideal screening modality, HCC is usually diagnosed at an advanced stage. Recent studies show that loss or gain of genomic materials can activate the oncogenes or inactivate the tumor suppressor genes to predispose cells toward carcinogenesis. Here, we evaluated both the copy number alteration (CNA) and RNA sequencing data of 361 HCC samples in order to locate the frequently altered chromosomal regions and identify the affected genes.

Results

Our data show that the chr1q and chr8p are two hotspot regions for genomic amplifications and deletions respectively. Among the amplified genes, YY1AP1 (chr1q22) possessed the largest correlation between CNA and gene expression. Moreover, it showed a positive correlation between CNA and tumor grade. Regarding deleted genes, CHMP7 (chr8p21.3) possessed the largest correlation between CNA and gene expression. Protein products of both genes interact with other cellular proteins to carry out various functional roles. These include ASH1L, ZNF496, YY1, ZMYM4, CHMP4A, CHMP5, CHMP2A and CHMP3, some of which are well-known cancer-related genes.

Conclusions

Our in-silico analysis demonstrates the importance of copy number alterations in the pathology of HCC. These findings open a door for future studies that evaluate our results by performing additional experiments.

Background

Hepatocellular carcinoma (HCC) is the most common type of liver cancer that occurs predominantly in patients with chronic liver disease and cirrhosis. The well-established risk factors of HCC are hepatitis B virus infection, alcoholism and metabolic disorders [1, 2]. Although liver lesions are usually detectable on computed tomography (CT), many HCC tumors are asymptomatic. As a result, HCC is usually diagnosed after development of end-stage clinical symptoms [3]. Despite the numerous researches, no effective screening modality exists for HCC. This necessitates identifying the key genetic factors and new potential therapeutic targets. However, this is a complicated process because of the heterogeneous nature of HCC [4].

Recent findings demonstrate that Copy number alteration (CNA) can lead to activation of oncogenes and inactivation of tumor suppressor genes in various cancers, which in turn predispose the cells toward malignancy. These genetic alterations include various chromosomal losses or gains that can cause the genome copy number of the affected cells to deviate from normal diploid state [5, 6]. Subsequently, this affects the genome stability and drives the cells toward tumorigenesis [7]. Chromosomal aberrations are associated with poor prognosis and function (age of onset, metastatic state, drug resistant and liver failure phenotypes) [8]. Therefore, it is helpful to identify the hotspot chromosomal region at which these genetic alterations are concentrated as well as the list of the genes that are affected by these alterations [9].

Recent advances in functional genomics along with the invention of comparative genomic hybridization (CGH) and microarray-based techniques permit the scientists to characterize the cytogenetic signatures of the malignant cells. The emergence of bioinformatic tools adds another level of proficiency that is needed for translating high-throughput data to improve cancer research, speeding up the discovery of hundreds of new potential targets, and finally, facilitating the developments of new prognostic, diagnostic and therapies approaches [10, 11]. Here, by means of an in silico analysis, we have evaluated the patterns of CNA in liver HCC samples to identify hotspot chromosomal regions that are altered during the course of disease. Our research provides the preliminary support for further laboratory experiments.

Results

Chromosomal cytobands chr1q and chr8q were strongly associated with gene amplification in HCC samples

Of 24,776 examined genes, 1024 (4.1%) were amplified in HCC samples (linear copy number values cut-off ≥ 0.5). These genes were mapped on chr1q (86%; n = 880), chr8q (13%; n = 131) and chr1p (1.3%; n = 13) (Fig. 1). Figure 2 shows the frequency with which different cytobands are observed for all the 1024 amplified genes. The top 34 genes with the high amplification scores were ADAM15, FLAD1, ZBTB7B, PYGO2, DCST2, LENEP, SHC1, PBXIP1, PMVK, DCST1, KCNN3, SLC50A1, DPM3, SCAMP3, FAM189B, MTX1, KRTCAP2, GBA, GBAP1, MUC1, THBS3, TRIM46, ASH1L, RUSC1, CLK2, HCN3, FDPS, PKLR, YY1AP1, MSTO1, GON4L, RIT1, SYT11 and KIAA0907 (Table 1). Majority of these 34 amplified genes are located on chr1q22 (67.6%; n = 23) while the remaining are on chr1q21.3 (32.4%; n = 11). Amplification of 13 genes were associated with tumor grade (Table 2). However, we observed no association between gene amplification and other clinicopathologic parameters, such as tumor size, stage and metastasis. Although all the 34 selected genes were amplified in HCC samples, we only observed a strong correlation between the CNA of oncogene YY1AP1 (Yin Yang-1 Associated Protein 1) and its expression (72% of studied samples; r > 0.6). PANTHER classified the protein products of these 34 genes in 15 cellular processes (Table 1), with many of them being well-recognized signaling pathways. These pathways include Wnt signaling pathway, angiogenesis signaling pathway, Ras signaling pathways, EGF signaling pathway, FGF signaling pathway, interleukin signaling pathway, and PDGF signaling pathway. We also noticed the involvement of some other pathways in HCC that are less emphasized in literature, if not completely ignored. Around 12.5% of these proteins are associated with cholesterol biosynthesis. Other pathways include flavin signaling pathway, glycolysis pathway, inflammation mediated by chemokine and cytokine signaling pathway, integrin signaling pathway, pyruvate metabolism pathway, and synaptic vesicle trafficking.

Fig. 1
figure1

Human chromosome ideograms showing regions of copy number variation (deletion and/or amplification) in a high resolution (244 K) whole genome hybridization (aCGH) analysis of 361 subjects with HCC. Deletions are shown in red and amplification in green. Data was extracted from cBioPortal dataset (https://www.cbioportal.org) through Bioconductor package cgdsr and the ideogram was generated by Bioconductor package IdeoViz

Fig. 2
figure2

Presentation of hotspot regions for all amplified genes on chr 1 and chr 8. 361 HCC samples posessing both expression data as well as CNV data were used to draw this plot. Data was extracted from cBioPortal dataset (https://www.cbioportal.org)

Table 1 The role of 34 amplified genes and 17 deleted genes with the highest scores according to PANTHER (http://pantherdb.org/)
Table 2 The association of copy number variations with clinicopathologic parameters of 361 HCC samples. Only candidate genes with significant association have been presented in the table

Analysis of protein interaction networks predicts 21 protein–protein interactions involving YY1AP1 and other proteins either directly or through intermediary. The top five interacting partners with scores greater than 0.6 are ASH1L, ZNF496 (Zinc Finger Protein 496), YY1 (Yin Yang 1), ZMYM4 (Zinc Finger MYM-Type Containing 4) and ACAD9 (Acyl-CoA dehydrogenase family member 9) (Fig. 3). Functional analysis shows that some of these genes belong to critical cellular components such as Nucleolar Remodeling Complexes (NoRC) that negatively regulates rRNA expression, and transcriptional machinery (Table 3).

Fig. 3
figure3

The top protein intractions with YY1AP1 protein, obtained from STRING version 11 (r > 0.5; cyan = from curated databases; pink = experimentally determined; olive = textmining; dark azure = gene co-occurrence; black = co-expression)

Table 3 The functional role of proteins interacting with YY1AP1 and CHMP7 according to g:profiler

MicroRNA target prediction also identified two miRNAs—hsa-miR-375 and hsa-miR-222-3p—that target YY1AP1, thereby they are able to regulate the concentration of YY1AP1 protein in the cell (Table 4).

Table 4 Human miRNAs target human YY1AP1 and CHMP7. Data were extracted from mirWalk dataset (http://mirtarbase.cuhk.edu.cn/php/search.php)

Cytoband chr8p was strongly associated with gene deletion in HCC samples

Of 24,776 examined genes, 172 (0.69%) were lost in samples (n = 361) (linear copy number values cut-off ≥ 0.5). Majority of these genes are located on chr8p21-23 (Fig. 1). We selected the top 17 genes with the high deletion scores (Table 1)—TNFRSF10B, RHOBTB2, PEBP4, TNFRSF10C, CHMP7, TNFRSF10A, ENTPD4, EGR3, BIN3, PDLIM2, R3HCC1, LOXL2, STC1, PIWIL2, SLC25A37, TNFRSF10D and CSMD1—for further analysis. Among these 17 genes, 16 are on chr8p21.3 (94%) and one on chr8p23.2 (6%) (Fig. 4). These genes encode proteins that function as either transcription factors or signaling/cytoskeletal proteins and take part in several processes including EGF, FGF, apoptosis, and P53 signaling pathways (Table 1). We didn’t observe any strong association between deletions of these 17 genes with tumor metastasis or grade. Although all the mentioned genes were downregulated in studied cancerous samples, only CHMP7 (Charged Multivesicular Body Protein 7) gene showed a moderate correlation between CNA and expression values (r = 0.5) in 70% of studied samples.

Fig. 4
figure4

The frequency of cytogenetic bands involved in all deleted genes on chromosme 8. 361 HCC samples posessing both expression data as well as CNV data were used to draw this plot. Data were extracted from cBioPortal dataset (https://www.cbioportal.org)

STRING database predicts direct protein–protein interactions between the CHMP7 protein and ten other proteins, among which the top five interactions with scores greater than 0.9 are with CHMP4A, CHMP5, CHMP2A, CHMP3 and ENSG00000249884 (RNF103-CHMP3 gene) (Fig. 5). According to functional analysis, some of these genes belong to vital cellular pathways including endosome transport, cytokinesis, vacuolar transport and viral assembly. (Table 3).

Fig. 5
figure5

The top protein interactions with CHMP7, obtained from STRING version 11 (r > 0.9; cyan = from curated databases; pink = experimentally determined; olive = textmining; dark azure = gene co-occurrence; black = co-expression)

MicroRNA target prediction also recognized four miRNAs hsa-miR-26b-5p, hsa-miR-505-5p, hsa-miR-484 and hsa-miR-15b-5p that target CHMP7 to influence its concentration in the cell (Table 4).

Discussion

We identified two hotspot regions chr1q and chr8p as locations that carry most of the chromosomal gains and losses respectively. This finding is in agreement with results from earlier studies. Comparative genomic hybridization on hepatocellular carcinoma samples shows involvement of both of these chromosomal regions, among several other regions [12]. Takafumi Nishimura et al. also noticed that one of the most common features of HCC is the recurrent chromosomal gain at chr1q [13]. Besides, Chen and coworkers discovered that chr1q and chr8p are among the frequently occurring gains and losses in HCC respectively [14]. Other independent studies have also highlighted the involvement of chr8p cytoband in HCC [15, 16]. Intriguingly, Yosuke Kishimoto et al. claim that HCC may develop from cirrhotic cells carrying chr8p loss [17].

We were also able to identify list of the genes that are influenced by these copy number alterations. Majority of these amplified and deleted genes are involved in various pathways, several of which are repeatedly reported as critical routs in HCC [18, 19]. Whereas amplified genes were located on either chr1q, chr1p or chr8q, deleted genes were located on chr8p. We noticed that 13 out of top 34 amplified genes are associated with tumor grade. Tae-Min Kim et al. have also found that gain in 8q11.21-q24.3 is associated with the high tumor grade [5]. On the other hand, Xin‐Yuan Guan et al. have reported an association between gain of 1q and 8q and stage of hepatocellular carcinoma [20] and Shirley M-H Sy and colleagues suggested a higher incidence of regional 1q21–q22 in association with the advanced stage tumors [21]. Despite the reported association of these chromosomal regions with tumor stage in general, we didn’t find any significant association between amplification of selected 34 genes and tumor stage.

Combing the RNA-seq and CNA data allowed us to identify subset of these genes that are simultaneously affected by both copy number alterations and gene expression pattern. Around 70% of the samples with gain of YY1AP1 or loss of the CHMP7 were also upregulated or downregulated respectively. YY1AP1 protein is a component of the INO80 chromatin remodeling complex, which is responsible for transcriptional regulation, DNA repair, and replication [22]. YY1AP1 acts as a critical oncoprotein in HCC by promoting the cell proliferation; maintaining stem cell features and epigenetically regulating of transcriptional networks. Moreover, silencing YY1AP1 eliminates the oncogenic properties of HCC cells by altering the chromatin structure and triggers apoptosis in vitro [4]. Sela et al. have also noticed that YY1AP1 possess exons that are intronic within normal tissues and are recognized as exons only within cancerous tissues. These exons are known as Alu exons [23] which are derived from Alu elements: the most common transposable element in the primate genomes [24]. This would explain the significance of YY1AP1 as a gene with both altered copy number and expression values.

Interestingly of five top proteins that interact with YY1AP1, four—ASH1L, ZNF496, YY1 and ZMYM4—are either transcription factors or transcriptional activators, while ACAD9 is involved in the beta-oxidation of fatty acyl-CoA in mitochondria. ASH1L is a histone methyltransferase that catalyzes the dimethylation and trimethylation of H3K36 and therefore is a transcriptional activator [25]. Silencing ASH1L through RNAi increases the sensitivity of HepG2 cells to sorafenib, suggesting that ASH1L may play a role in drug resistance in HCC cells [26]. In lung cancer, ASH1L also stimulates migration of cancer cells through Cdk5/p35 pathway [27]. YY1 is a famous oncogene widely involved in the development of many cancers including hepatocellular carcinogenesis. YY1AP1 acts as a coactivator of YY1 transcription factor [4]. YY1 is a critical component of epigenetic regulatory networks that dictates progression of hepatocellular carcinoma and suppresses the differentiation of hepatocytes [28, 29]. Both ZMYM4 (Zinc Finger MYM-Type Containing 4) and ZNF496 (Zinc finger protein 496) are identifiable from serum samples of patients with HCC [30]. Unfortunately, we don’t know anything about their function apart from the fact that they are zinc finger proteins overexpressed in multiple cancers [31, 32].

Adding to the complexity of interactions, two microRNAs can target YY1AP1 with opposite consequences. While hsa-miR-375 generally has a tumor suppressor role [33], hsa-miR-222-3p shows oncogenic properties [34].

Considering all the aspects, YY1AP1 is a key player in HCC with complicated transcriptional and epigenetic roles, exerting its effect either directly or through binding partners. This scheme is consistent with our finding that YY1AP1 shows significant changes in both copy number and gene expression data.

CHMP7 is a key component of ESCRT-III machinery that is involved in membrane abscission during cytokinesis, endosomal sorting, plasma membrane repair and nuclear envelope reforming [35]. Although, there are no publication regarding the importance of CHMP7 gene in HCC, activation of ESCRT-III machinery leads to membrane repair by shedding damaged parts of cell membranes. ESCRT-III is associated with resistance to cell death, especially when treated with anticancer agents [36].

In the same way as CHMP7, interacting proteins CHMP4A, CHMP5, CHMP2A, CHMP3 are components of mammalian ESCRT-III endosomal sorting complex while ENSG00000249884 is the result of readthrough transcription of CHMP3 and the adjacent gene RNF103. Apart from CHMP4A, the rest of these proteins show changes in their expression pattern in HCC [37, 38]. Our functional analysis suggests that CHMP7 and its interacting proteins promotes mitotic cytokinetic, midbody abscission and multivesicular body organization. Interestingly, some of these functional roles involve viral infection, which might explain the significance of ESCRT-III complex in progression of HCC. Infection by hepatitis B and C viruses is a major cause of hepatocellular carcinoma [39, 40]. Both of these viruses are dependent to ESCRT-III complex for their production, maturation and releases [41,42,43]. In the TCGA HCC dataset that we studied, 46% of total patients were in stage I alone, with less than 0.1% of them in late stage IV, suggesting that down-regulation of ESCRT-III is happening mostly at early stages of HCC but its significance is unclear.

We also found that four microRNAs target CHMP7: hsa-miR-26b-5p, hsa-miR-505-5p, hsa-miR-484 and hsa-miR-15b-5p. hsa-miR-26b-5p is involved in hepatitis B virus mediated HCC [44]; hsa-miR-484 shows both oncogenic and tumor-suppressor properties depending on the interacting partners [45, 46] and miR-15b-5p is a potential tumor suppressor [47, 48]. Altogether, the complex network of interactions suggests that CHMP7 could play important role(s) in mediating the HCC.

Conclusions

We combined RNA-seq and copy number data for hepatocellular carcinoma. Copy number alteration mainly affected the chromosomal regions chr1q and chr8p, mostly influencing 34 genes with gains and 17 genes with losses. These genes are linked to various pathways that are famous for their role in HCC. Among these, YY1AP1 and CHMP7 genes demonstrated significant changes in both copy number and gene expression data. Considering these findings, we suggest that candidate regions chr1q and chr8p in HCC could be subjects of further researches, with a major emphasis on the role of two genes YY1AP1 and CHMP7.

Methods

As chromosomal CNAs may affect the level of gene expression, the CNA and RNA-seq data were analyzed in parallel. The data for Liver Hepatocellular Carcinoma (TCGA, Provisional) possessing 24,777 genes and 361 samples were downloaded from cBioPortal [49, 50] through the “cgdsr” package [51] and analyzed in R v3.5. The CNA data were filtered by considering the linear copy number values cut-off ≥ 0.5 and cut-off ≤ -0.5 as thresholds for gene amplification and deletion respectively. The filtered results were further examined to identify the genes with significant correlation between their CNA and expression values. As in strong positive correlation, the linear correlation coefficient (r) is close to + 1, the results were filtered based on the r > 0.7.

The association of CNV variants with pathological tumor grade and stage was examined in R. The class type of selected proteins and their involvement in various pathways were estimated using the PANTHER classification system [52]. STRINGDB [53] was employed to identify the proteins that have functional and physical interactions with genes that possess both significant CNA and expression changes. Involvement of these proteins as groups in various biological processes were assessed by g:Profiler [54]. MiRWalk [55] was used to predict the presence of seed sequence in the selected genes. Two databases, one for oncogenes (http://ongene.bioinfo-minzhao.org/) and one for tumor suppressor genes (https://bioinfo.uth.edu/TSGene/), helped to identify if any of the selected genes are listed among oncogenes/tumor suppressor genes.

Availability of data and materials

The datasets used and/or analyzed during the current study are available on cbioportal.org and are also accessible through cgdsr R package.

References

  1. 1.

    Tahmasebi Birgani M, Carloni V. Tumor microenvironment, a paradigm in hepatocellular carcinoma progression and therapy. Int J Mol Sci. 2017;18(2):405.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  2. 2.

    Mazzocca A, Birgani MT, Sabba C, Carloni V. Tetraspanin-enriched microdomains and hepatocellular carcinoma progression. Cancer Lett. 2014;351(1):23–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  3. 3.

    Patel N, Yopp AC, Singal AG. Diagnostic delays are common among patients with hepatocellular carcinoma. J Natl Compr Cancer Netw. 2015;13(5):543–9.

    Article  Google Scholar 

  4. 4.

    Zhao X, Parpart S, Takai A, Roessler S, Budhu A, Yu Z, et al. Integrative genomics identifies YY1AP1 as an oncogenic driver in EpCAM(+) AFP(+) hepatocellular carcinoma. Oncogene. 2015;34(39):5095–104.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Kim TM, Yim SH, Shin SH, Xu HD, Jung YC, Park CK, et al. Clinical implication of recurrent copy number alterations in hepatocellular carcinoma and putative oncogenes in recurrent gains on 1q. Int J Cancer. 2008;123(12):2808–15.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    Albertson DG, Collins C, McCormick F, Gray JW. Chromosome aberrations in solid tumors. Nat Genet. 2003;34(4):369–76.

    CAS  PubMed  Article  Google Scholar 

  7. 7.

    Vargas-Rondon N, Villegas VE, Rondon-Lagos M. The role of chromosomal instability in cancer and therapeutic responses. Cancers (Basel). 2017;10(1):4.

    Article  CAS  Google Scholar 

  8. 8.

    Sansregret L, Swanton C. The role of aneuploidy in cancer evolution. Cold Spring Harb Perspect Med. 2017;7(1):1–18.

    Article  CAS  Google Scholar 

  9. 9.

    Shlien A, Malkin D. Copy number variations and cancer. Genome Med. 2009;1(6):62.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  10. 10.

    Wu D, Rice CM, Wang X. Cancer bioinformatics: a new approach to systems clinical medicine. BMC Bioinform. 2012;13:71.

    Article  Google Scholar 

  11. 11.

    Monni O, Hyman E, Mousses S, Barlund M, Kallioniemi A, Kallioniemi OP. From chromosomal alterations to target genes for therapy: integrating cytogenetic and functional genomic views of the breast cancer genome. Semin Cancer Biol. 2001;11(5):395–401.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    Rao UN, Gollin SM, Beaves S, Cieply K, Nalesnik M, Michalopoulos GK. Comparative genomic hybridization of hepatocellular carcinoma: correlation with fluorescence in situ hybridization in paraffin-embedded tissue. Mol Diagn. 2001;6(1):27–37.

    CAS  PubMed  Article  Google Scholar 

  13. 13.

    Nishimura T, Nishida N, Itoh T, Komeda T, Fukuda Y, Ikai I, et al. Discrete breakpoint mapping and shortest region of overlap of chromosome arm 1q gain and 1p loss in human hepatocellular carcinoma detected by semiquantitative microsatellite analysis. Genes Chromosom Cancer. 2005;42(1):34–43.

    CAS  PubMed  Article  Google Scholar 

  14. 14.

    Chen YJ, Yeh SH, Chen JT, Wu CC, Hsu MT, Tsai SF, et al. Chromosomal changes and clonality relationship between primary and recurrent hepatocellular carcinoma. Gastroenterology. 2000;119(2):431–40.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Roessler S, Long EL, Budhu A, Chen Y, Zhao X, Ji J, et al. Integrative genomic identification of genes on 8p associated with hepatocellular carcinoma progression and patient survival. Gastroenterology. 2012;142(4):957–66.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Qin LX, Tang ZY, Sham JS, Ma ZC, Ye SL, Zhou XD, et al. The association of chromosome 8p deletion and tumor metastasis in human hepatocellular carcinoma. Cancer Res. 1999;59(22):5662–5.

    CAS  PubMed  Google Scholar 

  17. 17.

    Kishimoto Y, Shiota G, Wada K, Kitano M, Nakamoto K, Kamisaki Y, et al. Frequent loss in chromosome 8p loci in liver cirrhosis accompanying hepatocellular carcinoma. J Cancer Res Clin Oncol. 1996;122(10):585–9.

    CAS  PubMed  Article  Google Scholar 

  18. 18.

    Dimri M, Satyanarayana A. Molecular signaling pathways and therapeutic targets in hepatocellular carcinoma. Cancers. 2020;12(2):491.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  19. 19.

    Whittaker S, Marais R, Zhu A. The role of signaling pathways in the development and treatment of hepatocellular carcinoma. Oncogene. 2010;29(36):4989–5005.

    CAS  PubMed  Article  Google Scholar 

  20. 20.

    Guan XY, Fang Y, Sham JS, Kwong DL, Zhang Y, Liang Q, et al. Recurrent chromosome alterations in hepatocellular carcinoma detected by comparative genomic hybridization. Genes Chromosom Cancer. 2000;29(2):110–6.

    CAS  PubMed  Article  Google Scholar 

  21. 21.

    Sy SM, Wong N, Lai PB, To K-F, Johnson PJ. Regional over-representations on chromosomes 1q, 3q and 7q in the progression of hepatitis B virus-related hepatocellular carcinoma. Mod Pathol. 2005;18(5):686–92.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Guo DC, Duan XY, Regalado ES, Mellor-Crummey L, Kwartler CS, Kim D, et al. Loss-of-function mutations in YY1AP1 lead to grange syndrome and a fibromuscular dysplasia-like vascular disease. Am J Hum Genet. 2017;100(1):21–30.

    CAS  PubMed  Article  Google Scholar 

  23. 23.

    Sela N, Mersch B, Gal-Mark N, Lev-Maor G, Hotz-Wagenblatt A, Ast G. Comparative analysis of transposed element insertion within human and mouse genomes reveals Alu’s unique role in shaping the human transcriptome. Genome Biol. 2007;8(6):1–19.

    Article  CAS  Google Scholar 

  24. 24.

    Ule J. Alu elements: at the crossroads between disease and evolution. Biochem Soc Trans. 2013;41(6):1532–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  25. 25.

    Miyazaki H, Higashimoto K, Yada Y, Endo TA, Sharif J, Komori T, et al. Ash1l methylates Lys36 of histone H3 independently of transcriptional elongation to counteract polycomb silencing. PLoS Genet. 2013;9(11):e1003897.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  26. 26.

    Li G-M, Wang Y-G, Pan Q, Wang J, Fan J-G, Sun C. RNAi screening with shRNAs against histone methylation-related genes reveals determinants of sorafenib sensitivity in hepatocellular carcinoma cells. Int J Clin Exp Pathol. 2014;7(3):1085.

    PubMed  PubMed Central  Google Scholar 

  27. 27.

    Demelash A, Rudrabhatla P, Pant HC, Wang X, Amin ND, McWhite CD, et al. Achaete-scute homologue-1 (ASH1) stimulates migration of lung cancer cells through Cdk5/p35 pathway. Mol Biol Cell. 2012;23(15):2856–66.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Zhang S, Jiang T, Feng L, Sun J, Lu H, Wang Q, et al. Yin Yang-1 suppresses differentiation of hepatocellular carcinoma cells through the downregulation of CCAAT/enhancer-binding protein alpha. J Mol Med. 2012;90(9):1069–77.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Rong Z, Wang Z, Wang X, Qin C, Geng W. Molecular interplay between linc01134 and YY1 dictates hepatocellular carcinoma progression. J Exp Clin Cancer Res. 2020;39:1–15.

    Article  Google Scholar 

  30. 30.

    Gao H-J, Chen Y-J, Zuo D, Xiao M-M, Li Y, Guo H, et al. Quantitative proteomic analysis for high-throughput screening of differential glycoproteins in hepatocellular carcinoma serum. Cancer Biol Med. 2015;12(3):246.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Cibis H, Biyanee A, Dörner W, Mootz HD, Klempnauer K-H. Characterization of the zinc finger proteins ZMYM2 and ZMYM4 as novel B-MYB binding proteins. Sci Rep. 2020;10(1):1–14.

    Article  CAS  Google Scholar 

  32. 32.

    Revill K, Wang T, Lachenmayer A, Kojima K, Harrington A, Li J, et al. Genome-wide methylation analysis and epigenetic unmasking identify tumor suppressor genes in hepatocellular carcinoma. Gastroenterology. 2013;145(6):1424–35.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  33. 33.

    Chang Y, Yan W, He X, Zhang L, Li C, Huang H, et al. miR-375 inhibits autophagy and reduces viability of hepatocellular carcinoma cells under hypoxic conditions. Gastroenterology. 2012;143(1):177–87.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  34. 34.

    Wang X, Liao X, Huang K, Zeng X, Liu Z, Zhou X, et al. Clustered microRNAs hsa-miR-221-3p/hsa-miR-222-3p and their targeted genes might be prognostic predictors for hepatocellular carcinoma. J Cancer. 2019;10(11):2520–33.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Olmos Y, Perdrix-Rosell A, Carlton JG. Membrane binding by CHMP7 coordinates ESCRT-III-dependent nuclear envelope reformation. Curr Biol. 2016;26(19):2635–41.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Liu J, Kang R, Tang D. ESCRT-III-mediated membrane repair in cell death and tumor resistance. Berlin: Nature Publishing Group; 2020.

    Google Scholar 

  37. 37.

    Liu GM, Zeng HD, Zhang CY, Xu JW. Key genes associated with diabetes mellitus and hepatocellular carcinoma. Pathology. 2019;215(11):152510.

    CAS  Google Scholar 

  38. 38.

    Li Y, Chen G, Yan Y, Fan Q. CASC15 promotes epithelial to mesenchymal transition and facilitates malignancy of hepatocellular carcinoma cells by increasing TWIST1 gene expression via miR-33a-5p sponging. Eur J Pharmacol. 2019;860:172.

    Google Scholar 

  39. 39.

    Chan SL, Wong V, Qin S, Chan H. Infection and cancer: the case of hepatitis B. J Clin Oncol. 2016;34(1):83–90.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  40. 40.

    Axley P, Ahmed Z, Ravi S, Singal AK. Hepatitis C virus and hepatocellular carcinoma: a narrative review. J Clin Transl Hepatol. 2018;6(1):79.

    PubMed  Google Scholar 

  41. 41.

    Ariumi Y, Kuroki M, Maki M, Ikeda M, Dansako H, Wakita T, et al. The ESCRT system is required for hepatitis C virus production. PLoS ONE. 2011;6(1):e14517.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Corless L, Crump CM, Griffin SD, Harris M. Vps4 and the ESCRT-III complex are required for the release of infectious hepatitis C virus particles. J Gen Virol. 2010;91(2):362–72.

    CAS  PubMed  Article  Google Scholar 

  43. 43.

    Lambert C, Döring T, Prange R. Hepatitis B virus maturation is sensitive to functional inhibition of ESCRT-III, Vps4, and γ2-adaptin. J Virol. 2007;81(17):9050–60.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Cui XW, Qian ZL, Li C, Cui SC. Identification of miRNA and mRNA expression profiles by PCR microarray in hepatitis B virus-associated hepatocellular carcinoma. Mol Med Rep. 2018;18(6):5123–32.

    CAS  PubMed  Google Scholar 

  45. 45.

    Qiu L, Huang Y, Li Z, Dong X, Chen G, Xu H, et al. Circular RNA profiling identifies circ ADAMTS 13 as a miR-484 sponge which suppresses cell proliferation in hepatocellular carcinoma. Mol Oncol. 2019;13(2):441–55.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    Yi R, Feng J, Yang S, Huang X, Liao Y, Hu Z, et al. miR-484/MAP2/c-Myc-positive regulatory loop in glioma promotes tumor-initiating properties through ERK1/2 signaling. J Mol Hist. 2018;49(2):209–18.

    CAS  Article  Google Scholar 

  47. 47.

    Yang Y, Hou N, Wang X, Wang L, Se C, He K, et al. miR-15b-5p induces endoplasmic reticulum stress and apoptosis in human hepatocellular carcinoma, both in vitro and in vivo, by suppressing Rab1A. Oncotarget. 2015;6(18):16227.

    PubMed  PubMed Central  Article  Google Scholar 

  48. 48.

    Li H, Zhang J, Lee M-J, Yu G-R, Han X, Kim D-G. OIP5, a target of miR-15b-5p, regulates hepatocellular carcinoma growth and metastasis through the AKT/mTORC1 and β-catenin signaling pathways. Oncotarget. 2017;8(11):18129.

    PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. In.: AACR; 2012.

  50. 50.

    Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  51. 51.

    Jacobsen A, Luna A. cgdsr: R-Based API for Accessing the MSKCC Cancer Genomics Data Server (CGDS), 2015. R package version. 2015;1(5). https://CRAN.R-project.org/package=cgdsr.

  52. 52.

    Thomas PD, Campbell MJ, Kejariwal A, Mi H, Karlak B, Daverman R, et al. PANTHER: a library of protein families and subfamilies indexed by function. Genome Res. 2003;13(9):2129–41.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  53. 53.

    Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.

    CAS  Article  Google Scholar 

  54. 54.

    Raudvere U, Kolberg L, Kuzmin I, Arak T, Adler P, Peterson H, et al. g: Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. 2019;47(W1):W191–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  55. 55.

    Dweep H, Gretz N. miRWalk2.0: a comprehensive atlas of microRNA-target interactions. Nat Methods. 2015;12(8):697–697.

    CAS  PubMed  Article  Google Scholar 

Download references

Acknowledgements

All the authors are thankful to Prof.Vinicio Carloni at university of Florence for generous assistance in manuscript revision and clinical analysis of data.

Funding

The present study was approved by Ethical Committee of Ahvaz Jundishapur University of Medical Sciences research affair, Ahvaz, Iran (Grant Number: CMRC-9613 and code of ethics: IR.AJUMS.REC.1396.367) and experimentally performed in Department of Medical Genetics.

Author information

Affiliations

Authors

Contributions

Study conception and design: MT-B; Acquisition of data: AS; Analysis and interpretation of data: AS, MT-B; Drafting of manuscript: HA, ZM and JM-A; Critical revision: VC; AS. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Maryam Tahmasebi-Birgani.

Ethics declarations

Ethics approval and consent to participate

The present study was approved by Ethical Committee of Ahvaz Jundishapur University of Medical Sciences research affair, Ahvaz, Iran (Grant number: CMRC-9613 and code of ethics: IR.AJUMS.REC.1396.367) and experimentally performed in Department of Medical Genetics.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher's Note

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.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Shahrisa, A., Tahmasebi-Birgani, M., Ansari, H. et al. The pattern of gene copy number alteration (CNAs) in hepatocellular carcinoma: an in silico analysis. Mol Cytogenet 14, 33 (2021). https://doi.org/10.1186/s13039-021-00553-2

Download citation

Keywords

  • Hepatocellular carcinoma
  • Copy number
  • RNA dysregulation
  • chr8p
  • chr1q