Gene networks and transcriptional regulators associated with liver cancer development and progression
BMC Medical Genomics volume 14, Article number: 41 (2021)
Treatment options for hepatocellular carcinoma (HCC) are limited, and overall survival is poor. Despite the high frequency of this malignoma, its basic disease mechanisms are poorly understood. Therefore, the aim of this study was to use different methodological approaches and combine the results to improve our knowledge on the development and progression of HCC.
Twenty-three HCC samples were characterized by histological, morphometric and cytogenetic analyses, as well as comparative genomic hybridization (aCGH) and genome-wide gene expression followed by a bioinformatic search for potential transcriptional regulators and master regulatory molecules of gene networks.
Histological evaluation revealed low, intermediate and high-grade HCCs, and gene expression analysis split them into two main sets: GE1-HCC and GE2-HCC, with a low and high proliferation gene expression signature, respectively. Array-based comparative genomic hybridization demonstrated a high level of chromosomal instability, with recurrent chromosomal gains of 1q, 6p, 7q, 8q, 11q, 17q, 19p/q and 20q in both HCC groups and losses of 1p, 4q, 6q, 13q and 18q characteristic for GE2-HCC. Gene expression and bioinformatics analyses revealed that different genes and gene regulatory networks underlie the distinct biological features observed in GE1-HCC and GE2-HCC. Besides previously reported dysregulated genes, the current study identified new candidate genes with a putative role in liver cancer, e.g. C1orf35, PAFAH1B3, ZNF219 and others.
Analysis of our findings, in accordance with the available published data, argues in favour of the notion that the activated E2F1 signalling pathway, which can be responsible for both inappropriate cell proliferation and initial chromosomal instability, plays a pivotal role in HCC development and progression. A dedifferentiation switch that manifests in exaggerated gene expression changes might be due to turning on transcriptional co-regulators with broad impact on gene expression, e.g. POU2F1 (OCT1) and NFY, as a response to accumulating cell stress during malignant development. Our findings point towards the necessity of different approaches for the treatment of HCC forms with low and high proliferation signatures and provide new candidates for developing appropriate HCC therapies.
Hepatocellular carcinoma (HCC) is the sixth most prevalent cancer and the second most lethal tumour worldwide; its incidence continues to rise . Prognosis and treatment options for HCC are poor and mainly dependent on tumour stage  and histological grade, which reflects tumour biology . High-grade HCC can develop from low-grade HCC in a course of several months in 75% of patients . Notably, cancer tissues with different histological grades can be frequently found in individual HCC nodules , a phenomenon that provides a ‘nodule-in-nodule appearance’. It may be assumed by these histological findings that there is a stepwise development from low- to high-grade HCC based on molecular mechanisms that are still unknown. Likewise, a stepwise increase in chromosomal instability and aneuploidy occurs during progression from well differentiated to moderately or low differentiated HCC . Chromosomal instability predicts drug resistance and poor prognosis in multiple cancer types, but the causes of chromosomal instability in HCC and other tumour types are still poorly understood [7, 8]. By gene expression profiling, HCC may be divided into two main groups: one with upregulation of genes responsible for cell proliferation and anti-apoptosis and poor outcome, and the other with the opposite findings [9, 10].
Although a number of studies have been performed using individual analytical methods, a combined approach for the investigation of HCC at histological, genomic and gene expression levels with a bioinformatics search for transcriptional regulators and master regulators of signal transduction pathways is rare. Therefore, the aim of our study was to gain a closer insight into the network of histological alterations, chromosomal aberrations and gene expression profiles and apply bioinformatics tools to find possible key players in the process of malignant transformation. For this purpose, we analysed 23 HCC samples of different histological grades and adjacent non-tumourous (NT) liver tissues with regard to the above-mentioned aspects. With this approach, we identified different genes and gene regulatory networks linked to distinct biological features observed in HCCs with low and high proliferative profiles. These findings provide additional candidate genes for the development of novel, subtype-dependent options for HCC therapy. Based on our results and comprehensive analysis of available published data, we also hypothesise key molecular events that underlie HCC onset and progression.
The samples for this study comprised tissue from partial hepatectomy of 15 HCC patients treated at the Department of Abdominal Surgery of the Inselspital (Bern, Switzerland) in 2002–2008 (clinicopathological data in Additional file 1). Twenty-three tumour samples were collected from these surgical specimens, as well as 17 NT samples adjacent to the carcinoma nodules. Tissue samples were formalin fixed and paraffin embedded following standard procedures. Edmonson–Steiner grading  was applied using haematoxylin and eosin (H&E) stained slides. For molecular analyses, tumours and paired NT tissues were macro-dissected after defining the regions of interest by an experienced pathologist to ensure a content of tumour cells of at least 75%.
Gene expression analyses and data evaluation
Total RNA was isolated using the Qiagen AllPrep DNA/RNA FFPE kit strictly according to the manufacturer’s protocol. DNA was checked for contamination with qPCR using Brilliant SYBRGreen QPCR Master Mix (Stratagene) and qPCR normalisation primers for single-copy regions of non-coding genomic DNA from the SideStep II Cell Lysis Analysis Kit (Agilent). From the kit, Primer Set 1 (233 bp) and Primer Set 2 (244 bp) were used, and the PCR cycling program was 95 °C for 10 min followed by 40 cycles of 95 °C for 30 s, 60 °C for 30 s, and 72 °C for 30 s. PCR reactions were run on a Versant kPCR real-time PCR instrument (Siemens) with MxPro™ QPCR Software for the Mx3000P/3005P QPCR System (Stratagene). Analysis of RNA quality was performed using an Agilent RNA 6000 Nano Assay on an Agilent 2100 Bioanalyzer. RNA integrity numbers for the isolated RNAs were 2.1–2.5.
Further procedures strictly followed the Gene Expression FFPE Workflow v. 2.0.1 protocol developed by Agilent Technologies. Briefly, a complementary DNA (cDNA) library was generated from the total RNA using the TransPlex™ Whole Transcriptome Amplification Kit (Sigma). cDNA was labelled with Cy3 using the Genomic DNA ULS Labeling Kit (Agilent) and purified with the Agilent KREApure column. Next, 500 ng of Cy3-cDNA was hybridized to one-colour Agilent SurePrint G3 Gene Expression v2 8 × 60 K Microarray (G4851B). The array contained 50,559 probes covering over 40,000 transcripts, including lncRNAs and transcripts of unknown coding potential (TUCPs). Hybridization, washing, and scanning procedures were performed exactly as described in the One-Colour Microarray-Based Gene Expression Analysis Protocol, v. 6.6 (Agilent). Washed arrays were immediately scanned on the Agilent C scanner and Agilent Feature Extraction software v. 220.127.116.11 was used to extract signal intensities from images.
Microarray data evaluation was performed using the United States Food and Drug Administration’s genomic tool ArrayTrack with embedded PCA and HCA . Raw gene expression values were normalized using the 75th percentile of the overall signal value, as recommended by the manufacturer, and log2 transformed prior to further statistical and bioinformatic analyses. Genes with significant differences in their expressions were identified by one-way ANOVA in multiple-group testing, or by unpaired Welch’s t test in pair-wise comparisons. For multiple probe sets for one gene, the median value was calculated.
For functional analysis of differently expressed genes, we employed Pathway Enrichment Analysis implemented in ArrayTrack which is based on canonical pathways supplied by the KEGG database. Fisher’s exact test was applied to identify biological pathways with statistically significant enrichment in regulated genes, the cut off being P ≤ 0.05.
Bioinformatics search for transcriptional regulators and master regulators of gene networks
To identify potential transcriptional regulators and master regulator molecules, we applied promoter analysis and upstream network analysis of differentially expressed genes provided by the GeneXplain platform (https://genexplain.com/) with integrated MATCH™ tool , TRANSFAC®  and TRANSPATH®  databases (release 2.4.1).
Gene promoters were defined as sequences from − 1000 to + 100 relative to the transcription start sites (Ensemble database). Promoter sequences of co-regulated genes were searched for potential transcription factor binding sites (TFBSs) using MATCH™ tool and TRANSFAC® library of positional weight matrices. The frequency of putative TFBSs in the genes differentially expressed in tumours (Yes set) was compared to the corresponding value for the default set of 300 human housekeeping genes of which expression is kept unchanged (No set). Transcription factors with overrepresented binding sites in the Yes set versus No set (ratio > 1) were considered as potential regulators.
Master regulator molecules in the signal transduction pathways upstream of the dysregulated genes were identified using the curated database of pathways and protein interactions TRANSPATH® with default parameters. Molecules with the lowest Rank Sum values were considered as the most likely candidates.
Real-time polymerase chain reaction (qPCR)
For validation of gene expression data obtained with microarrays, we employed qPCR for selected genes. Equal amounts of cDNAs from 12 NTs, 8 GE1-HCCs and 12 GE2-HCCs were pooled into the corresponding groups, and 10–20 ng of each cDNA pool were used to perform three independent qPCR experiments, each experiment in triplicate. We used a Versant kPCR real-time PCR instrument (Siemens) with the MxPro™ QPCR Software for Mx3000P/3005P QPCR System (Stratagene). PCR assays were performed for MKI67, HJURP, KPNA2, XPOT and TSNARE1 (internal control) with predesigned TaqMan® Gene Expression Assays and the TaqMan® Gene Expression Master Mix (Applied Biosystems) according to the manufacturer’s protocol. Assays for non-coding RNAs H19, HNF-AS1, SNAR-A3, SNAR-B2, SNAR-H were performed with Brilliant II SYBR® Green QPCR Low Rox Master Mix (Agilent) according to the manufacturer’s ‘Recommended Protocol with Two-Step Cycling (All Targets)’ with subsequent ‘Dissociation Program for All Targets’. The SNAR gene primers were derived from Parrott et al. . We designed H19, HNF-AS1 and TSNARE1 (an internal control) primers with the NCBI Primer-BLAST tool (https://www.ncbi.nlm.nih.gov/tools/primer-blast/) and purchased them from Biomers (Ulm, Germany). The sequences are as follows: SNAR-A3, F 5′-AGCCATTGTGGCTCAGGC-3′, R 5′-TTTTCCGACCCATGTGGACC-3′; SNAR-B2, F 5′-GCCATTGTGGCTCCGGC-3′, R 5′-AACCCATGTGGACCAGGTTG-3′; SNAR-H, F 5′-CCACTGTGGCTCCGGC-3′, R 5′-AATGTGGACCAGGTTGGCCT-3′; H19, F 5′-GACGTGACAAGCAGGACATGA-3′, R 5′-TAAGGTGTTCAGGAAGGCCG-3′; HNF1A-AS1, F 5′-ACTAAAATTCGGGCGAGGCA-3′, R 5′-GACTGGCTGAAGGGACACTC-3′; and TSNARE1, F 5′-GAAGAAAATTGCAGAAAAGTCCAGA-3′, R 5′-GTCACTCCCGTTAAAGACCTTC-3′. All PCR products were checked with an agarose gel. We used the comparative threshold cycle method (ΔΔCt) to quantify the relative amount of product transcripts. ∆Ct values = ((Ct target gene)—(Ct reference gene)) from three independent experiments were evaluated with t test to identify the statistically significant differences between the means of the sample groups. We calculated the FC for each gene as 2−∆∆Ct values, where ΔΔCt = ΔCtHCC − ΔCtNT.
Array-based comparative genomic hybridisation (aCGH)
aCGH was performed exactly as previously described ; the detailed protocol is also freely accessible in the ArrayExpress database under E-MTAB-8886. Cy5-labelled genomic DNA from tissue samples and Cy3-labeled human male genomic DNA (Promega, Germany) as a reference were hybridised to Agilent G4450A arrays. We normalised, visualised and analysed data with the Agilent Genomic Workbench v7.0 software using ADM-2 algorithm. We set cut-offs for chromosomal amplifications/deletions to the minimum number of probes in aberration interval ≥ 3 and average absolute log2 ratio ≥ 0.2, which corresponds to the presence of single copy aberration in about one third of analysed cell population.
Nuclear size measurement
Nuclear size measurement was done using Leica Application Suite, v.3.7.0 on 2-µM thick H&E-stained tissue sections in the same tumour and non-tumour regions used for molecular analyses. At least 100 cells were measured.
Fluorescence in situ hybridisation (FISH).
FISH was performed on 4 µm thick tissue sections using locus-specific probes for centromeres, CEP1, CEP3, CEP7, CEP8, and CEP17 (all supplied by Abbott, Wiesbaden, Germany) according to the supplier’s instructions. An Axio Imager microscope (Zeiss, Jena, Germany) was used for evaluation of fluorescent signals. A total of 100 nuclei were evaluated in each sample for calculation of the mean value.
Gene expression analysis divided HCC into two distinct groups: GE1 and GE2
Analysis of 23 HCC and 17 NT samples (Additional file 1) yielded 458 differently expressed probes with a mean absolute fold change (FC) ≥ 2 and false discovery rate (FDR) ≤ 0.05. These probes represented 321 unique genes: 131 were upregulated and 190 were downregulated in tumours (Additional file 2). Principal component analysis (PCA) (Fig. 1a) and hierarchical clustering analysis (HCA) (Fig. 1b) based on the gene expression of 458 probes divided all samples into two main groups of 13 HCC samples and another group of 10 HCC and 17 NT samples. The full image for HCA analysis with gene expression values is shown in Additional file 3.
Gene expression patterns of 10 HCC samples (GE1-HCC) were rather similar to the NT tissues. In contrast, gene expression pattern of 13 HCCs (GE2-HCC) clearly differed from GE1-HCC and adjacent NT tissues. After analyses, one tumour sample, HCC-17, was re-classified as a cholangiocellular carcinoma (CCC). The expression profile of this sample was very close to GE1-HCC.
To identify genes specifically deregulated in the two HCC subtypes, we performed three pair-wise t tests—GE1-HCC versus NT, GE2-HCC versus NT and GE2-HCC vs GE1-HCC—and multiple comparison of GE1-HCC, GE2-HCC and NT groups using one-way analysis of variance (ANOVA). The cut-offs for a significant change in gene expression for GE2-HCC versus NT and GE2-HCC versus GE1-HCC were FDR ≤ 0.05 and mean absolute FC ≥ 2. For the GE1-HCC versus NT comparison, we applied a less stringent cut-off, namely P ≤ 0.01, instead of a more stringent FDR threshold to recover enough differentially expressed genes for further analyses. All genes identified in the three pairwise t tests showed a significant change in a multiple-group testing (P ≤ 0.05), with the exception of two genes. Additional files 4–6 list the details for each gene.
In general, there was a strong perturbation in gene expression with regard to the number of regulated genes as well as expression fold changes in the GE2-HCC group. In contrast, in GE1-HCC there were only moderate gene expression changes and far fewer differentially regulated genes compared with surrounding NT tissue. In GE1-HCC, there were 209 differentially expressed probes with a mean absolute FC in the range 2 to 4.9, despite applying the less stringent P value cut-off. These probes represented 141 unique genes: 26 upregulated and 115 downregulated. In GE2-HCC, 1254 probes were differentially regulated in comparison with NT adjacent tissues, with an absolute FC from 2 to 31.7. They represented 859 unique genes: 395 upregulated and 464 downregulated. Analysis of the differentially expressed genes in tumours showed that a part of them were significantly regulated in both groups of tumours (Fig. 2a).
The direction of the common gene expression changes, i.e. activation or repression, was concordant in both tumour groups, with the exception of long non-coding RNA H19, which was repressed in GE1-HCC and increased in GE2-HCC. Differentially expressed genes in GE1-HCC and GE2-HCC were further used for analyses of functional pathways, gene networks and gene promoters. Pathway enrichment analysis (Fisher’s exact test) revealed five pathways specifically deregulated in GE1-HCC (Fig. 2b). In contrast, 26 pathways were affected in GE2-HCC. The dysregulated pathways in GE1-HCC were linked to immune response and NAD metabolism, both reported to be involved in cancer [18, 19]. The prominent feature of GE2-HCC was strong perturbation in the cell cycle, in particular DNA replication and the p53 signalling pathway. Additionally, cytochrome P450–catalysed drug metabolism, retinoic acid metabolism, fatty acid metabolism, cysteine and methionine metabolism and the metabolism of several other amino acids were specifically affected in GE2-HCC. Pathways related to cancer and immune diseases were also significantly changed in GE2-HCC. There were only two commonly regulated pathways in both HCC subgroups: complement and coagulation cascades and arachidonic acid metabolism. These pathways are involved in immune defence and inflammation and reportedly have a significant role in the development and progression of numerous malignancies [20, 21]. In general, the functional analysis of genes regulated in tumours and associated pathways indicated that GE2-HCC tumours exhibited a high proliferative capacity and strongly perturbed cell metabolism. They likely represent a much more aggressive phenotype compared with GE1-HCC tumours.
Array-based comparative genomic hybridisation (aCGH) revealed deletions to be typical for GE2-HCC
Numerous chromosomal aberrations were identified, including both gains and losses, which were shared by 30–81.8% of the HCC samples; some of these chromosomal aberrations were also found in the adjacent non-tumorous tissue (Table 1, Additional files 7–8). There were no frequently occurring deletions in NT samples, but some gains, e.g., 7q, 8q, 11p, and 19p, were observed in 50%-69.2% of non-tumorous tissues, derived from both GE1 and GE2 groups.
When analysing the GE1-HCC and GE2-HCC subgroups separately, obvious differences in their genomic patterns are observed (Table 2, Additional files 9–10). Most gains were found repeatedly in both the GE1-HCC and GE2-HCC subgroups. In contrast, losses were typical for the GE2-HCC subgroup, in which they occurred in 40–80% of the samples. Some loci, e.g., 1p, 16q, and 17p, which demonstrated recurrent gains in GE1 samples showed recurrent losses in the GE2-HCC subgroup. In the GE1-HCC subgroup, the only repeatedly identified loss was at 8p, which was shared by 30–50% of the samples.
Some amplified chromosomal loci encode genes with tumour-promoting functions that were strongly upregulated in tumours, as denoted by gene expression analysis (Table 3).
Nuclear size was significantly increased in tumours, most notably in GE2-HCC
Nuclear size, which is an important diagnostic criterion for tumour grading, was significantly larger in all HCC samples when compared to the adjacent NT tissues (Fig. 3). The median nuclear size in NT tissues was 7.44 µM, which falls in the range reported for the average normal hepatocyte nucleus of 7–8 µM . Meanwhile, the median nuclear size values in the GE1-HCC and GE2-HCC subgroups fell above this range, 8.49 µM and 10.50 µM, respectively. The interquartile ranges were 0.4 µM, 1.4 µM, and 1.8 µM for the NT, GE1-HCC and GE2-HCC groups, respectively, indicating a much more variable nuclear size in the HCC groups.
GE2-HCC showed high-level aneuploidy
Fluorescence in situ hybridisation (FISH) for the five selected chromosomes revealed a mean number of signals in GE1-HCC close to NT samples, with the highest value of 2.73 (Table 4). In contrast, GE2-HCC, which typically also had a higher histological grade and enlarged nuclei compared with GE1-HCC, were highly aneuploid, showing a variable increased number of chromosomes (up to 5.44).
Notably, histological grading and molecular group assignment did not match completely. Whereas all grade 1 HCC samples were assigned to GE1-HCC, based on gene expression profiling, and all grade 3 or 4 HCC samples were allotted to GE2-HCC, grade 2 tumours were in both molecular groups. Furthermore, GE2-HCC included more samples with a hepatitis background, although there was no pivotal difference between the two groups. Overall, 77% of GE2-HCC and 56% of GE1-HCC samples (1 sample had unknown status) were derived from patients with positive hepatitis serology.
Validation of microarray gene expression changes by real-time polymerase chain reaction (qPCR)
Microarray gene expression data were confirmed by qPCR analyses for selected genes and non-coding RNAs that were overexpressed in the GE2-HCC subgroup. (Table 5). There was good concordance between qPCR and microarray data for most genes; in particular, the strong up regulation of these genes in GE2 HCC was confirmed. In a number of cases, the fold changes in gene expression obtained via the qPCR method were much larger than those identified via microarray analyses, possibly due to the higher sensitivity of the qPCR assay, especially in cases of genes expressed at low levels. For example, in the case of both MKI67 and HJURP, microarray analysis revealed significant increases in expression in the GE2-HCC subgroup only, whereas qPCR analysis detected upregulation in both of the HCC subgroups, although a much higher relative level was detected in the GE2-HCC subgroup as well.
Master regulator molecules of gene networks in GE1-HCC and GE2-HCC
Master regulator molecules in the signal transduction pathways, at a distance of up to 10 steps upstream of the dysregulated genes, were identified for up- and downregulated genes in the GE1-HCC and GE2-HCC subgroups (Table 6, Additional files 11–14) using the curated database of pathways and protein interactions, TRANSPATH®, integrated in the GeneXplain platform. The top three master regulators identified for the 26 genes upregulated in the GE1-HCC subgroup included two cell-surface proteins involved in pro-survival signaling in the liver and E2F1, a key regulator of the cell cycle. Moderate and strong overexpression of E2F1 was observed in the GE1-HCC and GE2-HCC subgroups, respectively. The top three master regulator molecules identified for the 395 upregulated genes in the GE2-HCC subgroup comprised essential regulators of progression through the mitotic phase of the cell cycle, and overexpression of the genes coding these proteins was specifically observed in the GE2-HCC subgroup.
The top regulatory molecules identified in the networks of genes downregulated in the GE1-HCC (115 input genes) and GE2-HCC (464 input genes) were related to reduced anti-tumour immune response and inflammation.
Transcriptional regulators of gene networks in GE1-HCC and GE2-HCC
A bioinformatic search for potential transcription factor binding sites (TFBSs) using the MATCH™ tool and TRANSFAC® library of positional weight matrices integrated in the GeneXplain platform identified a number of transcription factors potentially involved in the coordinated regulation of gene expression in the GE1-HCC and GE2-HCC subgroups (Additional files 15–18). The majority of the top potential transcriptional regulators identified have already been found to play roles in various cancers, including HCC (Tables 7 and 8), thus demonstrating the validity of our in silico findings.
The results in GE1-HCC (Table 7, Additional file 15) focus at first on androgen receptor signalling (AR) (Yes/No ratio = 14.1), whereas its role in GE2-HCC seems to be low, as indicated by a Yes/No ratio of 1.1 (Additional file 16). Furthermore, expression data showed AR downregulation (mean FC = -3.22) in advanced GE2-HCC (Additional file 5).
Upregulated genes in GE1-HCC also showed enrichment in TFBSs for progesterone (PGR). Unlike AR, PGR binding sites enrichment was found also in GE2-HCC (Table 8). We also identified potential transcriptional factors with known oncogenic roles in HCC, e.g. ZEB1, or in other cancer types, e.g. HES1, NR6A1, PATZ1, PAX4 and MTF1, in GE1-HCC. A novel finding was the transcription factor ZNF219; its binding site showed a high Yes/No ratio of 4.4 and was found in promoters of more than half of upregulated genes in GE1-HCC.
In GE2-HCC, top potential transcriptional activators included transcription factors with reported tumour-promoting activity in HCC and/or other cancer types, e.g., NKX6-1, POU2F1/2, cJUN, E2F1, HSF2, SRF, TFCP2 and NFY (Table 8).
Some transcription factors seem to have an important role in activation of gene expression in both GE1-HCC and GE2-HCC. These were general transcription factors (GTF2A1/2) with a role in general enhancement of transcription and PGR, as mentioned above. A few transcription factors with high rank in one group (Tables 7 and 8) might also play a less significant role in the other group, where they were ranked lower (Additional files 15 and 16). This phenomenon relates to ZEB1 and PATZ1, which binding sites were highly presented in genes upregulated in GE1-HCC but less so in GE2-HCC. By contrast, the rank of E2F1, GABP and TFCP2 in regulation of gene networks in liver carcinomas strongly increased from GE1-HCC to GE2-HCC.
As described by Edmondson and Steiner in 1954, HCC includes a heterogeneous group of carcinomas . Morphological dedifferentiation is accompanied by increased genomic instability, in particular chromosomal instability, and a worse prognosis [6,7,8]. In this study, we demonstrated that molecular tumour typing based on gene expression profiling correlated to histological grade and chromosomal instability—defining two main types of HCC—and may be used for detection of molecular markers to improve HCC therapy options.
DNA level/chromosomal instability
Previously reported HCC chromosomal changes include gains of 1q, 6p, 8q, 17q and 20q and losses of 1p, 4q, 6q, 8p and 13q. Furthermore, 4q and 13 q deletions have been repeatedly correlated to HCC dedifferentiation [17, 60,61,62]. We also observed these characteristic aberrations in our set of hepatic malignancies. Amplification of genomic material could contribute to ectopic expression of genes with essential roles in cell growth, survival, and proliferation (Table 3). Notably, while gains of genetic material were observed in both the GE1-HCC and GE2-HCC subgroups, losses were observed mostly in the GE2-HCC subgroup (i.e., high-grade tumours). We assume that the loss of genetic material observed in GE2-HCCs is a result of severely dysregulated control of mitotic mechanisms and chromosomal missegregation resulting from progressive malignant derailment. Interestingly, recurrent gains at 7q, 8q, 11p, 19p, and 19q were also observed in the adjacent NT tissues from both HCC sample groups and most likely represent initial cancer-predisposing events.
Similar to the results in our previous study of childhood HCCs  we observed recurrent gain of chromosome 19 in GE1-HCC (100% frequency) and GE2-HCC (60% frequency) as well as in 50% of adjacent NT tissues. This finding may indicate a particular role for this chromosomal aberration in the development of liver carcinogenesis.
Overall, gene expression analysis revealed two main HCC groups: GE1-HCC and GE2-HCC, with GE1-HCC closer to NT than to GE2-HCC—the latter typically having a higher histological grade and higher level of chromosomal instability (as described above). There were a number of dysregulated genes in GE1-HCC and GE2-HCC, suggesting their essential role in the initiation and development of liver carcinogenesis. By contrast, other genes that were only dysregulated in GE2-HCC indicated their specific roles in cancer progression.
Genes upregulated in GE1-HCC and GE2-HCC
The genes significantly dysregulated in GE1-HCC and GE2-HCC revealed changes in the same direction; most of the dysregulation was more pronounced in GE2-HCC. The only exception was the long non-coding RNA (lncRNA) H19, which is frequently reported as dysregulated in cancer including HCC, sometimes up- and sometimes downregulated [63, 64]. In our samples, H19 was highly expressed in NT, downregulated in GE1-HCC and upregulated in GE2-HCC. These findings support the notion about the essential but context-specific role of H19 during all steps of liver carcinogenesis. Most of the genes upregulated in both GE1-HCC and GE-2 HCC have been reported to exert an oncogenic action. For example, the lncRNA HNF1A-AS1 is functionally involved in various carcinomas including HCC [65,66,67,68,69]. Another example is MDK (encodes a secreted growth factor): its upregulation correlated with a worse prognosis and chemotherapy resistance in diverse malignant tumours including HCC [70,71,72]. MDK-targeted molecular therapy approaches are under investigation . Other genes with increased expression in GE2-HCC compared with GE1-HCC included transcription factor E2F1, with a crucial role in cell cycle regulation, BIRC5 (survivin), an inhibitor of apoptosis, and helicase RECQL4, essential for maintaining genomic stability. The overexpression of these genes is critical for the development of various malignancies (including HCC) [47, 74,75,76,77,78,79]. Moreover, there have been reported direct interactions of E2F1 and retinoblastoma protein (RB) with the promoter of RECQL4 in prostate cancer cells  and RECQL4 with BIRC5 in breast cancer cells . These data suggest a ↑E2F1 → ↑RECQL4 → ↑BIRC5 axis in oncogenic derailment. A new finding in our study was upregulation of C1orf35, a novel potential oncogene . To our knowledge, our study is the first to show its overexpression in solid tumours. Notably, we found the genomic loci coding for H19 (11p15.5), RECQL4 (8q24.3), E2F1 (20q11.2), BIRC5 (17q25.3) and C1orf35 (1q42.13) to be frequently amplified in both GE1-HCC and GE2-HCC. Furthermore, loci 11p15.5 and 8q24.3 were frequently amplified in NT tissues. These findings suggest the role of gene amplification in the overexpression of these genes.
Genes downregulated in GE1-HCC and GE2-HCC
The genes repressed in both GE1-HCC and GE2-HCC included those with supposed tumour-inhibiting activity, e.g. CXCL14, CCBE1, IGFBP1, RND3, BMP10 and COLEC10, which encode proteins involved in extracellular matrix remodelling, migration and the immune response. Their decreased expression in various cancers has been reported [81,82,83,84,85,86].
Genes upregulated in GE1-HCC
A few genes specifically upregulated in GE1-HCC have been found to promote carcinogenesis of various cancer types when overexpressed. Some examples are CACNA1H, DBP, B3GnT8, ASAP3 and CYHR1, with CACNA1H and DBP linked to chemoresistance [87, 88]. One exception is DNM3, a gene that encodes a protein with growth-inhibiting function in HCC cells .
Genes upregulated in GE2-HCC: cell proliferation
Most of the genes strongly upregulated only in the GE2-HCC subgroup (Additional files 5 and 6) encode proteins involved in the cell cycle and its regulation. Many of these genes have previously been linked to advanced stages of different cancer types including HCC. These genes enclose KIFC1 [90, 91], GPC3 , PEG10 [93,94,95], and UBE2C [96,97,98].
In GE2-HCC, there was marked upregulation (up to 17 fold) in the expression of several small non-coding RNAs—SNAR-A3, SNAR-B2, SNAR-H, SNAR-I, SNAR-D and SNAR-G2—located on 19q13.33, a genomic region that showed recurrent gain in HCC and NT samples. An in vivo study showed that the SNAR family members are associated with NF90, a double-stranded RNA (dsRNA) binding protein implicated in transcriptional and translational control  and suggested to be involved in viral replication, e.g. hepatitis C virus . ILF3, which encodes the NF90 protein, was also upregulated in GE2-HCC and located on the amplified 19p13.2 region. A recent study reported ectopic expression of some SNAR family members in HCC patients. This dysregulation promoted invasion and migration in HCC cells .
Genes upregulated in GE2-HCC: nuclear size and transport
The distinct feature of tumour cells is an enlarged nucleus, although the molecular mechanisms that regulate nuclear size in normal and cancer cells remains unknown. We found that the expression of two major nuclear transport receptor genes, importin alpha 1—also known as karyopherin alpha 2 (KPNA2) —and exportin (XPOT), were significantly increased specifically in GE2-HCC. These changes could lead to profound dysregulation of nuclear-cytoplasmic transport . KPNA2 is the most abundant nuclear import receptor for proteins. An excess of this protein directly increased the nuclear size in a Xenopus model . It contributes to enhanced nuclear entry of structural proteins necessary for building up the nuclear envelope. The signalling pathways are also dependent on the nuclear transport system, namely import of activated transcription factors and signalling proteins and export of target gene products. Thus, KPNA2 overexpression might strongly deregulate gene transcription associated with increased proliferation and survival in GE2 carcinomas. Increased KPNA2 expression has been found in HCC patients, and its overexpression correlated with the expression of essential mitotic proteins CCNB2 and CDK1 in HCC cells . These findings are consistent with our study. Similar results have been reported for exportin-t (encoded by XPOT), a nuclear export receptor responsible for most export of mature tRNA from the nucleus to cytoplasm. The enhanced export of mature tRNA could hasten protein synthesis in the cytoplasm. A study revealed XPOT upregulation in HCC tumour tissues, poor outcome for patients with the high expression of the protein and association of XPOT expression with expression of molecules that regulate the cell cycle (i.e. several cyclins and cyclin-dependent kinases) in HCC cells .
Genes upregulated in GE2-HCC: chromosomal instability
GE2-HCC showed a high level of chromosomal instability, with large genomic losses and aneuploidy. This dysfunction might be linked to overexpression of genes involved in DNA repair, e.g. HJURP and PARP1, and in chromosome segregation and spindle checkpoint function, e.g. CDC20, ANAPC1, ANAPC11, TTK, BUB1B, MAD2L1, KIFC1, PTTG1 and UBE2C. HJURP encodes a protein involved in the repair of double-stranded DNA breaks. Its overexpression leads to mitotic defects in human cells . HJURP has already been shown to be an independent prognostic marker in HCC [107, 108]. PARP1 is involved in microhomology-mediated end joining (MMEJ), an error-prone DNA repair pathway. PARP1 is aberrantly expressed in different cancers, where its overexpression rather than underexpression is associated with genomic instability, likely due to accumulation of errors through increased MMEJ . Higher PARP1 expression has been correlated to a higher stage in HCC . PARP1 inhibitors are already in clinical use to suppress various tumours . We also observed high expression of some other genes involved in the spindle checkpoint mechanism and chromosome segregation, e.g. CDC20, MAD2L1 and TTK. Their overexpression, as opposed to mutation or underexpression, has been found in cancer cells with genomic instability and mitotic spindle defects , and has also been reported in HCC, where it was linked to low patient survival [112, 113].
Master regulators in GE1-HCC
Master regulators identified for upregulated gene networks in GE1-HCC included membrane-type metalloproteinases (MT-MMPs), TGM2 and E2F1. TGM2 is a cell-surface protein involved in cell binding to fibronectin. Fibronectin is a reactive extracellular matrix component implicated in cell adhesion, migration and signalling in the liver . Abnormal fibronectin accumulation in the liver due to sustained injury or disease has been suggested to elicit aberrant cell signalling that promotes further tissue damage and contributes to tumour genesis [114, 115]. E2F1 plays a critical role in the control of cellular proliferation and can activate genes that promote cell proliferation or apoptosis. Pro-survival signalling might abrogate the pro-apoptotic capacity of E2F1 and thus favour E2F1-mediated cell proliferation. MMP15 (associated with the plasma membrane) and TGM2, besides their involvement in invasion and metastasis, can convey cell signalling to enhance cell survival [116, 117]. Therefore, MT-MMPs, TGM2 and E2F1 might be conjoined by pathological processes in the liver to trigger inappropriate cell proliferation and tumour development. In our study, there was a nearly twofold increase in E2F1 in GE1-HCC. We discuss its pivotal role in the process of liver transformation below.
Transcriptional regulators in GE1-HCC
A bioinformatics search identified AR as an important transcription factor involved in upregulation of genes in GE1-HCC, with BIRC5 (an apoptosis inhibitor) and C1orf35 (a transcript with transforming activity) as its possible target genes. Androgen-dependent and androgen-independent AR activation has been reported in hepatocarcinogenesis . However, anti-androgenic treatment did not benefit male patients with advanced HCC . Consistent with this observation, our analysis showed that in GE2-HCC, AR gene expression and its role in gene transcription (Yes/No ratio in regulated genes) were significantly reduced compared with GE1-HCC. Information about PGR role in HCC is rare . The results of our study support the notion that PGR may be implicated in HCC carcinogenesis at early and advanced stages of the disease. Understanding the roles of AR and PGR in HCC is important because AR and PGR inhibitors are already in clinical use [24, 27]. A number of other potential transcriptional regulators identified in GE1-HCC—e.g. ZEB1, HES1, NR6A1, PATZ1, PAX4 and MTF1—have a reported oncogenic role in cancer types other than HCC [28,29,30,31,32,33,34,35]. Our results implicate them in hepatocarcinogenesis too. Transcription factor ZNF219, whose potential binding sites are found in several genes with an obvious role in carcinogenesis, i.e. E2F1, MDK, REQL4 and C1orf35, has emerged as a new potential transcriptional regulator in GE1-HCC.
Master regulators in GE2-HCC
The most significant master regulators identified for upregulated gene networks in GE2-HCC are important regulators of the cell cycle, namely the G2/M transition and mitotic spindle formation. The kinase PKMYT1 inhibits the G2/M transition, allowing time for DNA damage recovery. It is overexpressed rather than underexpressed in HCC, and exerts oncogenic activity in HCC cells . Targeted inhibition of PKMYT1 activity would shorten the time between checkpoint abrogation and mitotic entry and, therefore, should preferably damage cancer cells with an already defective G1 checkpoint mechanism. Efforts have been made to develop small molecule inhibitors of PKMYT1 as a promising anti-cancer therapy approach . CCNB2 overexpression hyperactivates the centrosomal kinases AURKA and PLK1, leading to accelerated centrosome separation, lagging chromosomes and aneuploidy . The APC/C complex bound to the regulatory protein CDC20 initiates anaphase entry for chromosome separation, and CDC20 is a critical molecule in the spindle assembly checkpoint mechanism. The ubiquitin-conjugating enzyme UBE2C collaborates with the APC/C–CDC20 complex in the removal of mitotic regulators and thereby contributes to cell cycle progression. Increased CDC20 expression has been associated with defective spindle formation and progression of multiple tumours, including HCC . UBE2C overexpression also leads to chromosomal missegregation and tumour formation . Nath et al.  linked UBE2C overexpression and chromosomal instability in cancer cells to excess E2F1, which recruits the CDC20–APC/C to upregulate UBE2C transcription. This phenomenon highlights the role of deregulated RB–E2F1 pathway in premature anaphase, chromosomal abnormalities and aneuploidy beyond its well-documented role in G1/S progression. All aforementioned genes, i.e. CCNB2, PKMYT, some components of the APC/C complex (e.g. ANAPC1 and ANAPC11), CDC20, UBE2C, E2F1 and many others involved in progression and control of G2/M phase, were strongly upregulated in GE2-HCC. They might represent candidates for targeted therapy in advanced HCCs with high-proliferation signature.
Transcriptional regulators in GE2-HCC
Promoter analysis of genes upregulated in GE2-HCC suggests the activation of transcription factor NKX6.1. A few upregulated genes with potential TFBSs for NKX6-1 in their promoters are critical genes essential for cell proliferation (CDK1 and two replication-dependent histones), growth (KPNA2 and XPOT) and motility (ROBO1) and might mediate the oncogenic effects of NKX6-1. NKX6-1 overexpression in HCC is associated with progressive features and an unfavourable prognosis .
Transcription factor POU2F1 might be pro-oncogenic in multiple contexts. It has been implicated in cell proliferation, immune modulation, oxidative and cytotoxic stress resistance, metabolic regulation, stem cell identity and a poised transcription state [38, 125]. It has been hypothesised to integrate multiple signal inputs through diverse posttranslational modifications and interaction with various partners to direct changes in gene transcription in general, by opening chromatin, and in cell-specific contexts. In prostate cancer, POU2F1 is a co-regulator of AR, leading to AR hypersensitivity and driving androgen-independent cancer progression . A similar phenomenon might occur in liver carcinogenesis. Accumulation of DNA damage and cell stress during the disease progression might activate POU2F1, leading, in turn, to changes of gene transcription induced by other transcription factors through modifying their interaction with ligands and/or binding to DNA. Enhanced POU2F1 activity has been described in different epithelial malignancies as a result of focal amplification, increased mRNA, augmented protein level and enhanced activity . We found amplification of the 1q24.2 locus, which encodes POU2F1, in more than 80% of all HCCs, but its transcriptional level was not changed. These data suggest the involvement of other mechanisms in the enhancement of its transcriptional activity. POU2F1 has already emerged as a prognostic marker and a potential therapeutic target in several types of cancer including HCC . The fact that POU2F1 is important for the cell response to genotoxic and oxidative stress, but is not critical under standard conditions, makes its inhibition an attractive anti-cancer therapeutic approach. One tactic is targeting the POU2F1-binding DNA sequence by pyrrole-imidazole polyamides, which has been successfully implemented in mice .
The second highest transcriptional regulator in GE2-HCC was the proto-oncogene JUN, overexpression of which has been shown to reprogramme oestrogen receptor signalling in breast cancer cells and confer them resistance to tamoxifen . JUN overexpression in HCC has been linked to sorafenib resistance .
Another significant transcription factor was E2F1, the expression of which was also strongly increased in GE2-HCC. As discussed in the section about master regulator molecules, the E2F family plays a crucial role in controlling the cell cycle and the action of tumour suppressor proteins . When overexpressed, it can be responsible for promoting the cell cycle as well as genomic instability . Sustained activation of E2F1 and its target genes might be induced by impaired regulatory mechanisms and aberrant signalling due to chronic liver disorders or injury or by DNA tumour virus oncoproteins, e.g. the hepatitis B virus X protein known to target the Rb–E2F1 pathway . Subtle yet significant gains in E2F1 and E2F3 copy number rather than mutations have been found in advanced HCC, and a direct, cell-autonomous role for E2f activators in driving HCC in mice has been demonstrated . In our HCC cohort, the genetic locus 20q11.2, which encodes E2F1, was recurrently amplified in tumours, E2F1 expression was moderately and strongly upregulated in GE1-HCC and GE2-HCC, respectively.
To achieve growth regulation, E2F1 requires partners. One of these is NFY . More than a half of the genes upregulated in GE2-HCC with a high proliferation signature were enriched in CCAAT box and NFY binding motif for CCAAT-binding transcription factor NFY. This finding agrees with studies of other cancers where CCAAT boxes were over-represented and bound to NFY in variety of upregulated genes . NFY is an essential regulator of cell cycle progression. It activates transcription of numerous genes that regulate cell cycle, including E2F1. It also has a direct non-transcriptional role in the overall efficiency of DNA replication. NFY increases cell type–specific gene expression by promoting chromatin accessibility for cell type–specific master transcription factors . Therefore, NFY activation might contribute to a broad augmentation of gene expression changes, first of all, E2F1 overexpression, associated with cancer progression in GE2-HCC. Considerable efforts have being made to find compounds that specifically inhibit NFY activity in cancer cells .
The transcription factor GABP is a negative prognostic biomarker in HCC and inhibits HCC cell migration and invasion . On the other hand, GABP strongly binds and activates a mutant promoter of the telomerase reverse transcriptase gene (TERT), reactivating its expression and thereby increasing the replicative potential of tumour cells . Two specific non-coding mutations in the TERT promoter required for the activation of TERT transcription by GABP occur with high frequency in aggressive cancer types. In one study, these mutations were found in 47% of HCCs . In our study, the promoter analysis identified GABP as potential important transcriptional regulator in GE2-HCC, but we did not find a significant difference in TERT transcription among GE1-HCC, GE2-HCC and NT samples.
Two additional transcriptional regulators upregulated in GE2-HCC were SRF and TFCP2. Increased SRF and TFCP2 levels accelerate tumour cell migration and invasion and are linked to cancer progression [49, 55] and to high-grade HCC [50, 56]. TFCP2 is involved in regulation of cell proliferation, invasion, angiogenesis, metastasis and chemoresistance in HCC . Small molecular inhibitors of TCPF2 have emerged as promising potent and effective therapeutics for HCC .
Candidate molecules for targeted anti-tumour therapy in GE1-HCC and GE2-HCC
As described above, the GE1-HCC and GE2-HCC subgroups were determined based on specific changes in histological, DNA and RNA levels that reflect their different biology. As a result, we uncovered many potentially valuable targets for anti-cancer therapy, as already discussed earlier and listed in Table 9.
GE gene expression, HCC hepatocellular carcinoma, MR master regulator, TR transcriptional regulator
Our findings point to a demand for differentiated approaches for treatment of HCC with low and high proliferation. Whereas E2F1, MDK, BIRC5, IL10, PGR, C1orf35 and TFCP2 (among others) seem to represent promising targets for intervention in all HCCs, therapy directed against KIFC1, PAFAH1B3, PKMYT1, PEG10, PARP1, POU2F1, NFY or TNF might be effective specifically in advanced GE2-HCC. Conversely, anti-androgen therapy, anti-ZNF219 and SOSC3-peptidomimetics may have a therapeutic potential for HCCs with low proliferation signature.
The findings from this study, in accordance with abundant previously published data, argue for a pivotal role of dysregulation of the E2F1 pathway in liver carcinogenesis. This dysfunction, due to diverse pathological processes in the liver, is capable of initiating both inappropriate cell proliferation and chromosomal instability.
A dedifferentiation switch manifested by strong propagation of gene expression changes and genomic instability might be linked to turning on transcriptional co-regulators, e.g. POU2F1 (OCT1) and NFY, as a response to accumulating cell stress during malignant development (Fig. 4).
Activation of transcriptional co-activators might explain, at least in part, the insensitivity of advanced HCC to anti-androgen treatment.
Our findings highlight the demand for different approaches to treat HCC entities that exhibit low or high proliferation signatures and provide new, promising candidates for developing targeted HCC therapy.
Availability of data and materials
The transcriptomic dataset generated during the current study and the detailed protocols used have been deposited in the ArrayExpress database at EMBL-EBI (www.ebi.ac.uk/arrayexpress) under accession number E-MTAB-8887. The genomic dataset generated during the current study and the detailed protocols used are deposited in the ArrayExpress database at EMBL-EBI (www.ebi.ac.uk/arrayexpress) under accession number E-MTAB-8886.
Gene expression sub-group 1/2
Non-tumourous liver tissue
Array-based comparative genomic hybridisation
Principal component analysis
Hierarchical clustering analysis
False discovery rate
Transcription factor binding site
Villanueva A. Hepatocellular carcinoma. N Engl J Med. 2019;380:1450–62.
Liu C, Duan LG, Lu WS, Yan LN, Xiao GQ, Jiang L, et al. Prognosis evaluation in patients with hepatocellular carcinoma after hepatectomy: comparison of BCLC, TNM and Hangzhou criteria staging systems. PLoS ONE. 2014;9:e103228.
Zhou L, Rui JA, Zhou WX, Bin WS, Chen SG, Qu Q. Edmondson-Steiner grade: a crucial predictor of recurrence and survival in hepatocellular carcinoma without microvascular invasio. Pathol Res Pract. 2017;213:824–30.
Sugihara S, Nakashima O, Kojiro M, Majima Y, Tanaka M, Tanikawa K. The morphologic transition in hepatocellular carcinoma. A comparison of the individual histologic features disclosed by ultrasound-guided fine-needle biopsy with those of autopsy. Cancer. 1992;70:1488–92.
Kenmochi K, Sugihara S, Kojiro M. Relationship of histologic grade of hepatocellular carcinoma (HCC) to tumor size, and demonstration of tumor cells of multiple different grades in single small HCC. Liver. 1987;7:18–26.
Wilkens L, Flemming P, Gebel M, Bleck J, Terkamp C, Wingen L, et al. Induction of aneuploidy by increasing chromosomal instability during dedifferentiation of hepatocellular carcinoma. Proc Natl Acad Sci USA. 2004;101:1309–14.
Rao CV, Asch AS, Yamada HY. Frequently mutated genes/pathways and genomic instability as prevention targets in liver cancer. Carcinogenesis. 2017;38:2–11.
McClelland SE. Role of chromosomal instability in cancer progression. Endocr Relat Cancer. 2017;24:T23-31.
Lee JS, Chu IS, Heo J, Calvisi DF, Sun Z, Roskams T, et al. Classification and prediction of survival in hepatocellular carcinoma by gene expression profiling. Hepatology. 2004;40:667–76.
Ke K, Chen G, Cai Z, Huang Y, Zhao B, Wang Y, et al. Evaluation and prediction of hepatocellular carcinoma prognosis based on molecular classification. Cancer Manag Res. 2018;10:5291–302.
Edmondson HA, Steiner PE. Primary carcinoma of the liver. A study of 100 cases among 48,900 necropsies. Cancer. 1954;7:462–503.
Fang H, Harris SC, Su Z, Chen M, Qian F, Shi L, et al. ArrayTrack: an FDA and public genomic tool. Methods Mol Biol. 2009;563:379–98.
Kel AE, Gössling E, Reuter I, Cheremushkin E, Kel-Margoulis OV, Wingender E. MATCH: a tool for searching transcription factor binding sites in DNA sequences. Nucleic Acids Res. 2003;31:3576–9.
Matys V, Kel-Margoulis OV, Fricke E, Liebich I, Land S, Barre-Dirrie A, et al. TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes. Nucleic Acids Res. 2006;34:D108–10.
Krull M, Pistor S, Voss N, Kel A, Reuter I, Kronenberg D, et al. TRANSPATH: an information resource for storing and visualizing signaling pathways and their pathological aberrations. Nucleic Acids Res. 2006;34:D546–51.
Parrott AM, Tsai M, Batchu P, Ryan K, Ozer HL, Tian B, et al. The evolution and expression of the snaR family of small non-coding RNAs. Nucleic Acids Res. 2011;39:1485–500.
Tan L, Meier T, Kuhlmann M, Xie F, Baier C, Zhu Z, et al. Distinct set of chromosomal aberrations in childhood hepatocellular carcinoma is correlated to hepatitis B virus infection. Cancer Genet. 2016;209:87–96.
Gonzalez H, Hagerling C, Werb Z. Roles of the immune system in cancer: From tumor initiation to metastatic progression. Genes Dev. 2018;32:1267–84.
Yaku K, Okabe K, Hikosaka K, Nakagawa T. NAD metabolism in cancer therapeutics. Front Oncol. 2018;8:622.
Afshar-Kharghan V. The role of the complement system in cancer. J Clin Investig. 2017;127:780–9.
Yang P, Cartwright CA, Li J, Wen S, Prokhorova IN, Shureiqi I, et al. Arachidonic acid metabolism in human prostate cancer. Int J Oncol. 2012;41:1495–503.
Gallik S. Exercise 8. Cell fractionation II: isolation and visualization of liver cell nuclei using density gradient centrifugation and fluorescence microscopy. Cell Biol OLM 4.0. http://stevegallik.org/cellbiologyolm_Ex07_P02.html. Accessed 10 Dec 2020.
Antonarakis ES. AR signaling in humanmalignancies: Prostate cancer and beyond. Cancers. 2018;10:22.
Kanda T, Takahashi K, Nakamura M, Nakamoto S, Wu S, Haga Y, et al. Androgen receptor could be a potential therapeutic target in patients with advanced hepatocellular carcinoma. Cancers. 2017;9:43.
Hieb AR, Halsey WA, Betterton MD, Perkins TT, Kugel JF, Goodrich JA. TFIIA changes the conformation of the DNA in TBP/TATA complexes and increases their kinetic stability. J Mol Biol. 2007;372:619–32.
Check JH. The role of progesterone and the progesterone receptor in cancer. Expert Rev Endocrinol Metab. 2017;12:187–97.
Yeh YT, Chang CW, Wei RJ, Wang SN. Progesterone and related compounds in hepatocellular carcinoma: basic and clinical aspects. Biomed Res Int. 2013;2013:1–9.
Li X, Cao Y, Li M, Jin F. Upregulation of HES1 promotes cell proliferation and invasion in breast cancer as a prognosis marker and therapy target via the AKT pathway and EMT process. J Cancer. 2018;9:757–66.
Sánchez-Tilló E, Siles L, de Barrios O, Cuatrecasas M, Vaquero EC, Castells A, et al. Expanding roles of ZEB factors in tumorigenesis and tumor progression. Am J Cancer Res. 2011;1:897–912.
Zhou Y-M, Cao L, Li B, Zhang R-X, Sui C-J, Yin Z-F, et al. Clinicopathological significance of ZEB1 protein in patients with hepatocellular carcinoma. Ann Surg Oncol. 2012;19:1700–6.
Qin Y, Yu J, Zhang M, Qin F, Lan X. ZEB1 promotes tumorigenesis and metastasis in hepatocellular carcinoma by regulating the expression of vimentin. Mol Med Rep. 2019;19:2297–306.
Mathieu R, Evrard B, Fromont G, Rioux-Leclercq N, Godet J, Cathelineau X, et al. Expression screening of cancer/testis genes in prostate cancer identifies NR6A1 as a novel marker of disease progression and aggressiveness. Prostate. 2013;73:1103–14.
Keskin N, Deniz E, Eryilmaz J, Un M, Batur T, Ersahin T, et al. PATZ1 is a DNA damage-responsive transcription factor that inhibits p53 function. Mol Cell Biol. 2015;35:1741–53.
Li CG, Eccles MR. PAX genes in cancer; friends or foes? Front Genet. 2012;3:6.
Ji L, Zhao G, Zhang P, Huo W, Dong P, Watari H, et al. Knockout of MTF1 inhibits the epithelial to mesenchymal transition in ovarian cancer cells. J Cancer. 2018;9:4578–85.
Tseng I-C, Yeh MM, Yang C-Y, Jeng Y-M. NKX6-1 is a novel immunohistochemical marker for pancreatic and duodenal neuroendocrine tumors. Am J Surg Pathol. 2015;39:850–7.
Huang LL, Zhang Y, Zhang JX, He LJ, Lai YR, Liao YJ, et al. Overexpression of NKX61 is closely associated with progressive features and predicts unfavorable prognosis in human primary hepatocellular carcinoma. Tumor Biol. 2015;36:4405–15.
Vázquez-Arreguín K, Tantin D. The Oct1 transcription factor and epithelial malignancies: old protein learns new tricks. Biochim Biophys Acta Gene Regul Mech. 2016;1859:792–804.
Zhu HY, Cao GY, Wang SP, Chen Y, Liu GD, Gao YJ, et al. POU2F1 promotes growth and metastasis of hepatocellular carcinoma through the FAT1 signaling pathway. Am J Cancer Res. 2017;7:1665–79.
Wang SM, Tie J, Wang WL, Hu SJ, Yin JP, Yi XF, et al. POU2F2-oriented network promotes human gastric cancer metastasis. Gut. 2016;65:1427–38.
He H, Sinha I, Fan R, Haldosen LA, Yan F, Zhao C, et al. C-Jun/AP-1 overexpression reprograms ERα signaling related to tamoxifen response in ERα-positive breast cancer. Oncogene. 2018;37:2586–600.
Chen W, Xiao W, Zhang K, Yin X, Lai J, Liang L, et al. Activation of c-Jun predicts a poor response to sorafenib in hepatocellular carcinoma: preliminary clinical evidence. Sci Rep. 2016;6:22976.
Lin X, Li Y, Jin W, Han F, Lu S, Yu W, et al. LHX3 is an early stage and radiosensitivity prognostic biomarker in lung adenocarcinoma. Oncol Rep. 2017;38:1482–90.
Yanai H, Negishi H, Taniguchi T. The IRF family of transcription factors inception, impact and implications in oncogenesis. Oncoimmunology. 2012;1:1376–86.
Li P, Du Q, Cao Z, Guo Z, Evankovich J, Yan W, et al. Interferon-gamma induces autophagy with growth inhibition and cell death in human hepatocellular carcinoma (HCC) cells through interferon-regulatory factor-1 (IRF-1). Cancer Lett. 2012;314:213–22.
Denechaud PD, Fajas L, Giralt A. E2F1, a novel regulator of metabolism. Front Endocrinol. 2017;8:311.
Kent LN, Bae S, Tsai SY, Tang X, Srivastava A, Koivisto C, et al. Dosage-dependent copy number gains in E2f1 and E2f3 drive hepatocellular carcinoma. J Clin Investig. 2017;127:830–42.
Björk JK, Åkerfelt M, Joutsen J, Puustinen MC, Cheng F, Sistonen L, et al. Heat-shock factor 2 is a suppressor of prostate cancer invasion. Oncogene. 2016;35:1770–84.
Prencipe M, Fabre A, Murphy TB, Vargyas E, O’Neill A, Bjartell A, et al. Role of serum response factor expression in prostate cancer biochemical recurrence. Prostate. 2018;78:724–30.
Park MY, Kim KR, Park HS, Park B-H, Choi HN, Jang KY, et al. Expression of the serum response factor in hepatocellular carcinoma: implications for epithelial-mesenchymal transition. Int J Oncol. 2007;31:1309–15.
Bell RJA, Rube HT, Kreig A, Mancini A, Fouse SD, Nagarajan RP, et al. The transcription factor GABP selectively binds and activates the mutant TERT promoter in cancer. Science. 2015;348:1036–9.
Liu R, Zhang T, Zhu G, Xing M. Regulation of mutant TERT by BRAF V600E/MAP kinase pathway through FOS/GABP in human cancer. Nat Commun. 2018;9:579.
Zhang S, Zhang K, Zheng X, Jin J, Feng M, Liu P. GABPA predicts prognosis and inhibits metastasis of hepatocellular carcinoma. BMC Cancer. 2017;17:380.
Hu Q, Zhang B, Chen R, Fu C, Jun A, Fu X, et al. ZFHX3 is indispensable for ERβ to inhibit cell proliferation via MYC downregulation in prostate cancer cells. Oncogenesis. 2019;8:28.
Kotarba G, Krzywinska E, Grabowska AI, Taracha A, Wilanowski T. TFCP2/TFCP2L1/UBP1 transcription factors in cancer. Cancer Lett. 2018;420:72–9.
Yoo BK, Emdad L, Gredler R, Fuller C, Dumur CI, Jones KH, et al. Transcription factor late SV40 factor (LSF) functions as an oncogene in hepatocellular carcinoma. Proc Natl Acad Sci USA. 2010;107:8357–62.
Gurtner A, Manni I, Piaggio G. NF-Y in cancer: Impact on cell transformation of a gene essential for proliferation. Biochim Biophys Acta Gene Regul Mech. 2017;1860:604–16.
Cao B, Zhao Y, Zhang Z, Li H, Xing J, Guo S, et al. Gene regulatory network construction identified NFYA as a diffuse subtype-specific prognostic factor in gastric cancer. Int J Oncol. 2018;53:1857–68.
Mori Y, Olaru AV, Cheng Y, Agarwal R, Yang J, Luvsanjav D, et al. Novel candidate colorectal cancer biomarkers identified by methylation microarray-based scanning. Endocr Relat Cancer. 2011;18:465–78.
Hertz S, Rothämel T, Skawran B, Giere C, Steinemann D, Flemming P, et al. Losses of chromosome arms 4q, 8p, 13q and gain of 8q are correlated with increasing chromosomal instability in hepatocellular carcinoma. Pathobiology. 2008;75:312–22.
Moinzadeh P, Breuhahn K, Stützer H, Schirmacher P. Chromosome alterations in human hepatocellular carcinomas correlate with aetiology and histological grade—results of an explorative CGH meta-analysis. Br J Cancer. 2005;92:935–41.
Steinemann D, Skawran B, Becker T, Tauscher M, Weigmann A, Wingen L, et al. Assessment of differentiation and progression of hepatic tumors using array-based comparative genomic hybridization. Clin Gastroenterol Hepatol. 2006;4:1283–91.
Raveh E, Matouk IJ, Gilon M, Hochberg A. The H19 Long non-coding RNA in cancer initiation, progression and metastasis—a proposed unifying theory. Mol Cancer. 2015;14:184.
Ardelt MA, Pachmayr J. The long non-coding RNA H19 - a new player in hepatocellular carcinoma. Cell Stress. 2017;1:4–6.
Yang X, Song JH, Cheng Y, Wu W, Bhagat T, Yu Y, et al. Long non-coding RNA HNF1A-AS1 regulates proliferation and migration in oesophageal adenocarcinoma cells. Gut. 2014;63:881–90.
Wang C, Mou L, Chai HX, Wang F, Yin YZ, Zhang XY. Long non-coding RNA HNF1A-AS1 promotes hepatocellular carcinoma cell proliferation by repressing NKD1 and P21 expression. Biomed Pharmacother. 2017;89:926–32.
Fang C, Qiu S, Sun F, Li W, Wang Z, Yue B, et al. Long non-coding RNA HNF1A-AS1 mediated repression of miR-34a/SIRT1/p53 feedback loop promotes the metastatic progression of colon cancer by functioning as a competing endogenous RNA. Cancer Lett. 2017;410:50–62.
Wu Y, Liu H, Shi X, Yao Y, Yang W, Song Y. The long non-coding RNA HNF1A-AS1 regulates proliferation and metastasis in lung adenocarcinoma. Oncotarget. 2015;6:9160–72.
Liu Z, Wei X, Zhang A, Li C, Bai J, Dong J. Long non-coding RNA HNF1A-AS1 functioned as an oncogene and autophagy promoter in hepatocellular carcinoma through sponging hsa-miR-30b-5p. Biochem Biophys Res Commun. 2016;473:1268–75.
Ma M-C, Chen Y-J, Chiu T-J, Lan J, Liu C-T, Chen Y-C, et al. Positive expression of Midkine predicts early recurrence and poor prognosis of initially resectable combined hepatocellular cholangiocarcinoma. BMC Cancer. 2018;18:227.
Kishida S, Kadomatsu K. Involvement of midkine in neuroblastoma tumourigenesis. Br J Pharmacol. 2014;171:896–904.
Zhu W-W, Guo J-J, Guo L, Jia H-L, Zhu M, Zhang J-B, et al. Evaluation of midkine as a diagnostic serum biomarker in hepatocellular carcinoma. Clin Cancer Res. 2013;19:3944–54.
Muramatsu T, Kadomatsu K. Midkine: an emerging target of drug development for treatment of multiple diseases. Br J Pharmacol. 2014;171:811–3.
Kent LN, Leone G. The broken cycle: E2F dysfunction in cancer. Nat Rev Cancer. 2019;19:326–38.
Li D, Hu C, Li H. Survivin as a novel target protein for reducing the proliferation of cancer cells (review). Biomed Rep. 2018;8:399–406.
Su C. Survivin in survival of hepatocellular carcinoma. Cancer Lett. 2016;379:184–90.
Su Y, Meador JA, Calaf GM, Proietti De-Santis L, Zhao Y, Bohr VA, et al. Human RecQL4 helicase plays critical roles in prostate carcinogenesis. Cancer Res. 2010;70:9207–17.
Fang H, Nie L, Chi Z, Liu J, Guo D, Lu X, et al. RecQL4 helicase amplification is involved in human breast tumorigenesis. PLoS ONE. 2013;8:e69600.
Li J, Jin J, Liao M, Dang W, Chen X, Wu Y, et al. Upregulation of RECQL4 expression predicts poor prognosis in hepatocellular carcinoma. Oncol Lett. 2018;15:4248–54.
Luo SQ, Xiong DH, Li J, Li G, Wang Y, Zhang JM, et al. C1orf35 contributes to tumorigenesis by activating c-MYC transcription in multiple myeloma. Oncogene. 2020;39:3354–66.
Gu XL, Ou ZL, Lin FJ, Yang XL, Luo JM, Shen ZZ, et al. Expression of CXCL14 and its anticancer role in breast cancer. Breast Cancer Res Treat. 2012;135:725–35.
Mesci A, Huang X, Taeb S, Jahangiri S, Kim Y, Fokas E, et al. Targeting of CCBE1 by miR-330-3p in human breast cancer promotes metastasis. Br J Cancer. 2017;116:1350–7.
Dai B, Ruan B, Wu J, Wang J, Shang R, Sun W, et al. Insulin-like growth factor binding protein-1 inhibits cancer cell invasion and is associated with poor prognosis in hepatocellular carcinoma. Int J Clin Exp Pathol. 2014;7:5645–54.
Paysan L, Piquet L, Saltel F, Moreau V. Rnd3 in cancer: a review of the evidence for tumor promoter or suppressor. Mol Cancer Res. 2016;14:1033–44.
Ye L, Bokobza S, Li J, Moazzam M, Chen J, Mansel RE, et al. Bone morphogenetic protein-10 (BMP-10) inhibits aggressiveness of breast cancer cells and correlates with poor prognosis in breast cancer. Cancer Sci. 2010;101:2137–44.
Zhang B, Wu H. Decreased expression of COLEC10 predicts poor overall survival in patients with hepatocellular carcinoma. Cancer Manag Res. 2018;10:2369–75.
Fornaro L, Vivaldi C, Lin D, Xue H, Falcone A, Wang Y, et al. Prognostic relevance of a T-type calcium channels gene signature in solid tumours: a correlation ready for clinical validation. PLoS ONE. 2017;12:e0182818.
Huang Y-F, Wu Y-H, Cheng W-F, Peng S-L, Shen W-L, Chou C-Y. Vitamin D-binding protein enhances epithelial ovarian cancer progression by regulating the insulin-like growth factor-1/Akt pathway and vitamin d receptor transcription. Clin Cancer Res. 2018;24:3217–28.
Gu C, Yao J, Sun P. Dynamin 3 suppresses growth and induces apoptosis of hepatocellular carcinoma cells by activating inducible nitric oxide synthase production. Oncol Lett. 2017;13:4776–84.
Xiao Y-X, Yang W-X. KIFC1: a promising chemotherapy target for cancer treatment? Oncotarget. 2016;7:48656–70.
Fu X, Zhu Y, Zheng B, Zou Y, Wang C, Wu P, et al. KIFC1, a novel potential prognostic factor and therapeutic target in hepatocellular carcinoma. Int J Oncol. 2018;52:1912–22.
Zhou F, Shang W, Yu X, Tian J. Glypican-3: A promising biomarker for hepatocellular carcinoma diagnosis and treatment. Med Res Rev. 2018;38:741–67.
Deng X, Hu Y, Ding Q, Han R, Guo Q, Qin J, et al. PEG10 plays a crucial role in human lung cancer proliferation, progression, prognosis and metastasis. Oncol Rep. 2014;32:2159–67.
Bang H, Ha SY, Hwang SH, Park CK. Expression of PEG10 is associated with poor survival and tumor recurrence in hepatocellular carcinoma. Cancer Res Treat. 2015;47:844–52.
Xie T, Pan S, Zheng H, Luo Z, Tembo KM, Jamal M, et al. PEG10 as an oncogene: Expression regulatory mechanisms and role in tumor progression. Cancer Cell Int. 2018;18:112.
Xie C, Powell C, Yao M, Wu J, Dong Q. Ubiquitin-conjugating enzyme E2C: a potential cancer biomarker. Int J Biochem Cell Biol. 2014;47:113–7.
Wei ZI, Liu YI, Qiao SH, Li XU, Li QI, Zhao JI, et al. Identification of the potential therapeutic target gene ube2c in human hepatocellular carcinoma: an investigation based on geo and tcga databases. Oncol Lett. 2019;17:5409–18.
Cacciola NA, Calabrese C, Malapelle U, Pellino G, De Stefano A, Sepe R, et al. UbcH10 expression can predict prognosis and sensitivity to the antineoplastic treatment for colorectal cancer patients. Mol Carcinog. 2016;55:793–807.
Parrott AM, Mathews MB. Novel rapidly evolving hominid RNAs bind nuclear factor 90 and display tissue-restricted distribution. Nucleic Acids Res. 2007;35:6249–58.
Patiño C, Haenni AL, Urcuqui-Inchima S. NF90 isoforms, a new family of cellular proteins involved in viral replication? Biochimie. 2015;108:20–4.
Shi Z, Wei D, Wu H, Ge J, Lei X, Guo Z, et al. Long non-coding RNA snaR is involved in the metastasis of liver cancer possibly through TGF-β1. Oncol Lett. 2019;17:5565–71.
Stelma T, Chi A, Van Der Watt PJ, Verrico A, Lavia P, Leaner VD. Targeting nuclear transporters in cancer: diagnostic, prognostic and therapeutic potential. IUBMB Life. 2016;68:268–80.
Levy DL, Heald R. Nuclear Size Is Regulated by Importin α and Ntf2 in Xenopus. Cell. 2010;143:288–98.
Gao CL, Wang GW, Yang GQ, Yang H, Zhuang L. Karyopherin subunit-α2 expression accelerates cell cycle progression by upregulating CCNB2 and CDK1 in hepatocellular carcinoma. Oncol Lett. 2018;15:2815–20.
Lin J, Hou Y, Huang S, Wang Z, Sun C, Wang Z, et al. Exportin-T promotes tumor proliferation and invasion in hepatocellular carcinoma. Mol Carcinog. 2019;58:293–304.
Mishra PK, Au WC, Choy JS, Kuich PH, Baker RE, Foltz DR, et al. Misregulation of Scm3p/HJURP causes chromosome instability in saccharomyces cerevisiae and human cells. PLoS Genet. 2011;7:e1002303.
Hu B, Wang Q, Wang Y, Chen J, Li P, Han M. Holliday junction–recognizing protein promotes cell proliferation and correlates with unfavorable clinical outcome of hepatocellular carcinoma. Onco Targets Ther. 2017;10:2601–7.
Chen T, Huang H, Zhou Y, Geng L, Shen T, Yin S, et al. HJURP promotes hepatocellular carcinoma proliferation by destabilizing p21 via the MAPK/ERK1/2 and AKT/GSK3β signaling pathways. J Exp Clin Cancer Res. 2018;37:193.
Wang L, Liang C, Li F, Guan D, Wu X, Fu X, et al. PARP1 in carcinomas and PARP1 inhibitors as antineoplastic drugs. Int J Mol Sci. 2017;18:2111.
Li J, Dou D, Li P, Luo W, Lv W, Zhang C, et al. PARP-I serves as a novel molecular marker for hepatocellular carcinoma in a Southern Chinese Zhuang population. Tumor Biol. 2017;39:1–8.
Yuan B, Xu Y, Woo JH, Wang Y, Bae YK, Yoon DS, et al. Increased expression of mitotic checkpoint genes in breast cancer cells with chromosomal instability. Clin Cancer Res. 2006;12:405–10.
Fan G, Tu Y, Chen C, Sun H, Wan C, Cai X. DNA methylation biomarkers for hepatocellular carcinoma. Cancer Cell Int. 2018;18:140.
Miao R, Wu Y, Zhang H, Zhou H, Sun X, Csizmadia E, et al. Utility of the dual-specificity protein kinase TTK as a therapeutic target for intrahepatic spread of liver cancer. Sci Rep. 2016;6:33121.
Aziz-Seible RS, Casey CA. Fibronectin: functional character and role in alcoholic liver disease. World J Gastroenterol. 2011;17:2482–99.
Torbenson M, Wang J, Choti M, Ashfaq R, Maitra A, Wilentz RE, et al. Hepatocellular carcinomas show abnormal expression of fibronectin protein. Mod Pathol. 2002;15:826–30.
Abraham R, Schäfer J, Rothe M, Bange J, Knyazev P, Ullrich A. Identification of MMP-15 as an anti-apoptotic factor in cancer cells. J Biol Chem. 2005;280:34123–32.
Boroughs LK, Antonyak MA, Cerione RA. A novel mechanism by which tissue transglutaminase activates signaling events that promote cell survival. J Biol Chem. 2014;289:10115–25.
Trinchet JC. Randomized trial of leuprorelin and flutamide in male patients with hepatocellular carcinoma treated with tamoxifen. Hepatology. 2004;40:1361–9.
Liu L, Wu J, Wang S, Luo X, Du Y, Huang D, et al. PKMYT1 promoted the growth and motility of hepatocellular carcinoma cells by activating beta-catenin/TCF signaling. Exp Cell Res. 2017;358:209–16.
Schmidt M, Rohe A, Platzer C, Najjar A, Erdmann F, Sippl W. Regulation of G2/M transition by inhibition of WEE1 and PKMYT1 Kinases. Molecules. 2017;22:2045.
Nam H-J, van Deursen JM. Cyclin B2 and p53 control proper timing of centrosome separation. Nat Cell Biol. 2014;16:538–49.
Li J, Gao J-Z, Du J-L, Huang Z-X, Wei L-X. Increased CDC20 expression is associated with development and progression of hepatocellular carcinoma. Int J Oncol. 2014;45:1547–55.
Van Ree JH, Jeganathan KB, Malureanu L, Van Deursen JM. Overexpression of the E2 ubiquitin-conjugating enzyme UbcH10 causes chromosome missegregation and tumor formation. J Cell Biol. 2010;188:83–100.
Nath S, Chowdhury A, Dey S, Roychoudhury A, Ganguly A, Bhattacharyya D, et al. Deregulation of Rb-E2F1 axis causes chromosomal instability by engaging the transactivation function of Cdc20—anaphase-promoting complex/cyclosome. Mol Cell Biol. 2015;35:356–69.
Pance A. Oct-1, to go or not to go? That is the PolII question. Biochim Biophys Acta Gene Regul Mech. 2016;1859:820–4.
Obinata D, Takayama K, Fujiwara K, Suzuki T, Tsutsumi S, Fukuda N, et al. Targeting Oct1 genomic function inhibits androgen receptor signaling and castration-resistant prostate cancer growth. Oncogene. 2016;35:6350–8.
Wang W-H, Hullinger RL, Andrisani OM. Hepatitis B virus X protein via the p38MAPK pathway induces E2F1 release and ATR kinase activation mediating p53 apoptosis. J Biol Chem. 2008;283:25455–67.
Van Ginkel PR, Hsiao KM, Schjerven H, Farnham PJ. E2F-mediated growth regulation requires transcription factor cooperation. J Biol Chem. 1997;272:18367–74.
Oldfield AJ, Yang P, Conway AE, Cinghu S, Freudenberg JM, Yellaboina S, et al. Histone-fold domain protein NF-Y promotes chromatin accessibility for cell type-specific master transcription factors. Mol Cell. 2014;55:708–22.
Quaas A, Oldopp T, Tharun L, Klingenfeld C, Krech T, Sauter G, et al. Frequency of TERT promoter mutations in primary tumors of the liver. Virchows Arch. 2014;465:673–7.
Rajasekaran D, Siddiq A, Willoughby JLS, Biagi JM, Christadore LM, Yunes SA, et al. Small molecule inhibitors of Late SV40 Factor (LSF) abrogate hepatocellular carcinoma (HCC): Evaluation using an endogenous HCC model. Oncotarget. 2015;6:26266–77.
Hayatigolkhatmi K, Padroni G, Su W, Fang L, Gómez-Castañeda E, Hsieh YC, et al. Investigation of a minor groove-binding polyamide targeted to E2F1 transcription factor in chronic myeloid leukaemia (CML) cells. Blood Cells Mol Dis. 2018;69:119–22.
Garg H, Suri P, Gupta JC, Talwar GP, Dubey S. Survivin: a unique target for tumor therapy. Cancer Cell Int. 2016;16:49.
Sekino Y, Oue N, Koike Y, Shigematsu Y, Sakamoto N, Sentani K, et al. KIFC1 Inhibitor CW069 induces apoptosis and reverses resistance to docetaxel in prostate cancer. J Clin Med. 2019;8:225.
La Manna S, Lee E, Ouzounova M, Di Natale C, Novellino E, Merlino A, et al. Mimetics of suppressor of cytokine signaling 3: novel potential therapeutics in triple breast cancer. Int J Cancer. 2018;143:2177–86.
Oft M. IL-10: master switch from tumor-promoting inflammation to antitumor immunity. Cancer Immunol Res. 2014;2:194–9.
Wang M, Wu M, Yang T. The synergistic effect of sorafenib and TNF-α inhibitor on hepatocellular carcinoma. EBioMedicine. 2019;40:11–2.
We thank Dr. Olga Kel-Margoulis and Dr. Alexander Kel from the GeneXplain team for excellent support in the use of this bioinformatics platform.
Ethics approval and consent to participate
The study was conducted in compliance with ethical standards. The biological material and health-related data were studied anonymously. The use of archived tissue samples for retrospective analysis was approved by the Ethics Committee of the Canton of Bern (KEK 2014/200) which provided a waiver of the need for informed consent.
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.
. Clinicopathological data of the patients with liver carcinoma used in the study.
. Differentially expressed probe sets in hepatocellular carcinoma (HCC) versus non-tumourous (NT) tissues. Unpaired Welch’s t-test was performed using gene expression microarray data from 23 HCC and 17 NT samples. The probe sets differentially expressed in the HCC group (n = 23) versus the NT group (n = 17) showing a mean absolute fold change (FC) ≥ 2, with a false discovery rate (FDR) ≤ 0.05, were selected and used for hierarchical cluster analysis and principal component analysis.
. Two subtypes of hepatocellular carcinoma (HCC), with distinctive gene expression patterns, were defined by hierarchical cluster analysis. Gene expression values from 23 HCC samples of various histological grades and 17 surrounding non-tumourous (NT) tissues were analysed. For the analysis, 458 probe sets differently expressed in the HCC samples versus the NT tissues with a mean absolute ratio ≥ 2 and false discovery rate ≤ 0.05 (Additional file 2) were used. The scale in the top right shows the colour codes representing gene expression values (Log2) – green for low expression and red for high expression.
. List of differentially expressed probe sets in GE1-hepatocellular carcinoma (HCC) versus non-tumourous (NT) tissues. Unpaired Welch’s t-test was performed using gene expression microarray data from GE1-HCC (n=10) and NT (n=17) samples. The probe sets demonstrating expression changes in the GE1-HCC group versus the NT group with a mean absolute fold change (FC) ≥ 2 and P ≤ 0.01 were selected and used for further bioinformatic analyses.
. List of differentially expressed probe sets in 13 GE2-hepatocellular carcinoma (HCC) versus 17 non-tumourous (NT) tissues. Unpaired Welch’s t-test was performed using gene expression microarray data from GE2-HCC (n=13) and NT (n=17) samples. The probe sets demonstrating expression change in the GE2-HCC group versus the NT group with a mean absolute fold change (FC) ≥ 2 and false discovery rate (FDR) ≤ 0.05 were selected and used for further bioinformatics analyses.
. List of differentially expressed probe sets in 13 GE2-hepatocellular carcinoma (HCC) versus 10 GE1-HCC samples. Unpaired Welch’s t-test was performed using gene expression microarray data from GE2-HCC (n=13) and GE1-HCC (n=10) samples. The probe sets that expression change in the GE2-HCC group versus the GE1-HCC group showed mean absolute fold change (FC) ≥ 2 and false discovery rate (FDR) ≤ 0.05 were selected.
. Percentage of aberrations shared by multiple samples in the set of hepatocarcinoma (HCC) tissues. Genomic DNA isolated from HCC samples (n=22) was analysed using high resolution, microarray-based, comparative genomic hybridisation (aCGH). Human male genomic DNA (Promega) was used as a reference.
. Percentage of aberrations shared by multiple samples in the set of non-tumorous (NT) tissues. Genomic DNA isolated from NT samples (n=13) was analysed using high resolution, microarray-based, comparative genomic hybridisation (aCGH). Human male genomic DNA (Promega) was used as a reference.
. Percentage of aberrations shared by multiple samples in the set of hepatocarcinoma (HCC) of GE1 group. Genomic DNA isolated from GE1-HCC samples (n=10) was analysed using high resolution, microarray-based, comparative genomic hybridisation (aCGH). Human male genomic DNA (Promega) was used as a reference.
. Percentage of aberrations shared by multiple samples in the set of hepatocarcinoma (HCC) of GE2 group. Genomic DNA isolated from GE2-HCC samples (n=12) was analysed using high resolution, microarray-based, comparative genomic hybridisation (aCGH). Human male genomic DNA (Promega) was used as a reference.
. Master regulator molecules of genes upregulated in GE1-hepatocellular carcinoma (HCC). Significantly up regulated genes in GE1-HCC (n=10) versus NT (n=17) samples with mean fold change (FC) ≥ 2 and P ≤ 0.01 were identified (Additional file 4) and used as an input gene list for bioinformatics analyses. Master regulator molecules in the signal transduction pathways, at a distance of up to 10 steps upstream of the dysregulated genes, were identified using the curated database of pathways and protein interactions TRANSPATH®.
. Master regulator molecules of genes upregulated in GE2-hepatocellular carcinoma (HCC). Significantly up regulated genes in GE2-HCC (n=13) versus NT (n=17) samples with mean fold change (FC) ≥ 2 and false discovery rate (FDR) ≤ 0.05 were identified (Additional file 5) and used as an input gene list for bioinformatics analyses. Master regulator molecules in the signal transduction pathways, at a distance of up to 10 steps upstream of the dysregulated genes, were identified using the curated database of pathways and protein interactions TRANSPATH®.
. Master regulator molecules of genes downregulated in GE1-hepatocellular carcinoma (HCC). Significantly down regulated genes in GE1-HCC (n=10) versus NT (n=17) samples with mean fold change (FC) ≤ -2 and P ≤ 0.01 were identified (Additional file 4) and used as an input gene list for bioinformatics analyses. Master regulator molecules in the signal transduction pathways, at a distance of up to 10 steps upstream of the dysregulated genes, were identified using the curated database of pathways and protein interactions TRANSPATH®.
. Master regulator molecules of genes downregulated in GE2-hepatocellular carcinoma (HCC). Significantly down regulated genes in GE2-HCC (n=13) versus NT (n=17) samples with mean fold change (FC) ≤ -2 and false discovery rate (FDR) ≤ 0.05 were identified (Additional file 5) and used as an input gene list for bioinformatics analyses. Master regulator molecules in the signal transduction pathways, at a distance of up to 10 steps upstream of the dysregulated genes, were identified using the curated database of pathways and protein interactions TRANSPATH®.
. Lists of overrepresented transcription factor binding sites (TFBSs) in genes upregulated in GE1-hepatocellular carcinoma (HCC). Significantly upregulated genes in GE1-HCC (n=10) versus NT (n=17) samples with mean fold change (FC) ≥ 2 and P ≤ 0.01 (Additional file 4) were used as an input gene list for bioinformatics analyses. Promoter sequences of the co-regulated in tumours genes (“Yes” set) and the control gene set (“No” set) were searched for potential transcription factor binding sites (TFBSs) using the MATCHTM tool and TRANSFAC® library of positional weight matrices. The frequencies of TFBSs in the “Yes” gene set and “No” gene set were identified, and the ratio of “Yes” versus “No” was calculated. Yes/No ratio > 1 indicates overrepresentation of TFBSs in the co-regulated genes and suggests potential involvement of the corresponding transcription factor in regulation of the observed expression changes.
. Lists of overrepresented transcription factor binding sites (TFBSs) in genes upregulated in GE2-hepatocellular carcinoma (HCC). Significantly up regulated genes in GE2-HCC (n=13) versus NT (n=17) samples with mean fold change (FC) ≥ 2 and false discovery rate (FDR) ≤ 0.05 (Additional file 5) were used as an input gene list for bioinformatics analyses. Promoter sequences of the co-regulated in tumours genes (“Yes” set) and the control gene set (“No” set) were searched for potential transcription factor binding sites (TFBSs) using the MATCHTM tool and TRANSFAC® library of positional weight matrices. The frequencies of TFBSs in the “Yes” gene set and “No” gene set were identified, and the ratio of “Yes” versus “No” was calculated. Yes/No ratio > 1 indicates overrepresentation of TFBSs in the co-regulated genes and suggests potential involvement of the corresponding transcription factor in regulation of the observed expression changes.
. Lists of overrepresented transcription factor binding sites (TFBSs) in genes downregulated in GE1-hepatocellular carcinoma (HCC). Significantly down regulated genes in GE1-HCC (n=10) versus NT (n=17) samples with mean fold change (FC) ≤ -2 and P ≤ 0.01 (Additional file 4) were used as an input gene list for bioinformatics analyses. Promoter sequences of the co-regulated in tumours genes (“Yes” set) and the control gene set (“No” set) were searched for potential transcription factor binding sites (TFBSs) using the MATCHTM tool and TRANSFAC® library of positional weight matrices. The frequencies of TFBSs in the “Yes” gene set and “No” gene set were identified, and the ratio of “Yes” versus “No” was calculated. Yes/No ratio > 1 indicates overrepresentation of TFBSs in the co-regulated genes and suggests potential involvement of the corresponding transcription factor in regulation of the observed expression changes.
. Lists of overrepresented transcription factor binding sites (TFBSs) in genes downregulated in GE2-hepatocellular carcinoma (HCC). Significantly down regulated genes in GE2-HCC (n=13) versus NT (n=17) samples with mean fold change (FC) ≤ -2 and false discovery rate (FDR) ≤ 0.05 (Additional file 5) were used as an input gene list for bioinformatics analyses. Promoter sequences of the co-regulated in tumours genes (“Yes” set) and the control gene set (“No” set) were searched for potential transcription factor binding sites (TFBSs) using the MATCHTM tool and TRANSFAC® library of positional weight matrices. The frequencies of TFBSs in the “Yes” gene set and “No” gene set were identified, and the ratio of “Yes” versus “No” was calculated. Yes/No ratio > 1 indicates overrepresentation of TFBSs in the co-regulated genes and suggests potential involvement of the corresponding transcription factor in regulation of the observed expression changes.
About this article
Cite this article
Meier, T., Timm, M., Montani, M. et al. Gene networks and transcriptional regulators associated with liver cancer development and progression. BMC Med Genomics 14, 41 (2021). https://doi.org/10.1186/s12920-021-00883-5