Journal of Transplantation

Volume 2013 (2013), Article ID 878297, 14 pages

http://dx.doi.org/10.1155/2013/878297

Research Article

Molecular Signatures of Recurrent Hepatocellular Carcinoma Secondary to Hepatitis C Virus following Liver Transplantation

Trina Das,¹ Deborah L. Diamond,² Matthew Yeh,¹ Sajida Hassan,³ Janine T. Bryan,⁴ Jorge D. Reyes,¹ and James D. Perkins¹

¹Division of Transplantation, Department of Surgery, University of Washington School of Medicine, Seattle, WA 98195-6410, USA

²Department of Pharmaceutics, University of Washington, Seattle, WA 98121-1023, USA

³Department of Laboratory Medicine, University of Washington Medical Center, Seattle, WA 98104-2499, USA

⁴Department of Microbiology, University of Washington, Seattle, WA 98195-7735, USA

Received 3 July 2013; Accepted 25 September 2013

Academic Editor: P. Burra

Copyright © 2013 Trina Das et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Chronic hepatitis C virus (HCV) induced hepatocellular carcinoma (HCC) is a primary indication for liver transplantation (LT). In western countries, the estimated rate of HCC recurrence following LT is between 15% and 20% and is a major cause of mortality. Currently, there is no standard method to treat patients who are at high risk for HCC recurrence. The aim of this study was to investigate the molecular signatures underlying HCC recurrence that may lead to future studies on gene regulation contributing to new therapeutic options. Two groups of patients were selected, one including patients with HCV who developed HCC recurrence (HCC-R) ≤3 years from LT and the second group including patients with HCV who did not have recurrent HCC (HCC-NR). Microarray analysis containing more than 29,000 known genes was performed on formalin-fixed-paraffin-embedded (FFPE) liver tissue from explanted livers. Gene expression profiling revealed 194 differentially regulated genes between the two groups. These genes belonged to cellular networks including cell cycle G1/S checkpoint regulators, RAN signaling, chronic myeloid leukemia signaling, molecular mechanisms of cancer, FXR/RXR activation and hepatic cholestasis. A subset of molecular signatures associated with HCC recurrence was found. The expression levels of these genes were validated by quantitative PCR analysis.

1. Introduction

Hepatitis C virus (HCV) infection is the most common cause of hepatocellular carcinoma (HCC) in the USA, Europe, and Japan, accounting for 47%–49%, 56%, and 75% of cases, respectively [1, 2]. HCC causes >600,000 deaths annually worldwide and is the most common primary liver cancer [3]. Definitive treatment for HCC is surgical resection when possible or liver transplantation for patients with end-stage liver disease and liver tumors. In USA, transplant guidelines stipulate that eligibility for liver transplantation is determined by the patient’s liver tumor(s) meeting the Milan criteria (a single tumor ≤5 cm in diameter or up to 3 tumors with individual diameters ≤3 cm and no macrovascular invasion) [4]. Unfortunately, the recurrence of HCC is a major cause of mortality in surgically treated patients [5]. There is no standard therapy for patients who are at high risk for HCC recurrence. Hence, a better understanding of the molecular mechanisms involved in the recurrence of HCC post LT is necessary to develop an efficient surveillance protocol and seek new potential therapies.

Gene expression profiling is best performed on fresh or frozen tissue to lessen the degradation of RNA. However, collecting and storing this tissue are burdensome and costly, and obtaining Institutional Review Board (IRB) approval for longitudinal studies can be unsuccessful. Sample sizes are typically small, and collection requires years of time. In contrast, all clinical pathology laboratories have huge storage files containing formalin-fixed paraffin-embedded (FFPE) tissue. Formalin fixation has been used for several decades for preserving tissue. FFPE tissue would supply a huge resource for genomic studies to utilize previously collected samples with long-term clinical records. However, the formalin and storage can negatively influence the integrity of RNA extracted from the FFPE tissues [6]. RNA derived from FFPE specimens is partially degraded, resulting from the variability in tissue handling/processing, tissue sources, and RNA extraction methods [7]. However, advanced RNA isolation techniques have been developed [8]. Likewise, advanced technologies using a range of microarray platforms allow improved microarray gene expression profiling using FFPE tissues [9, 10]. This present study uses the advanced technologies of the complementary DNA-mediated annealing, selection, extension, and ligation (DASL) assay (Illumina, Inc., San Diego, CA, USA) to analyze FFPE stored liver tissue to look for potential molecular signatures of recurrent HCC.

Several gene-profiling studies have identified molecular signatures associated with HCC metastatic potential [11–13]. However, few common genes were identified in these studies [11–13]. This may have been due to the use of differing microarrays and gene selection procedures (i.e., algorithms), types of tissue sample studied, patients’ clinical variability, the use of HCCs of different viral origins, that is, hepatitis B or C virus (HBV or HCV), which are major causes of human HCC, and to the small number of specimens studied [14, 15]. Most of these studies were performed on biopsy liver tissue samples. Hoshida et al. demonstrated the feasibility of genome-wide expression profiling of FFPE tissues using the DASL assay only for 6000 transcriptionally informative genes [16, 17]. This present feasibility study uses the method of transcriptome profiling of FFPE tissues to analyze 29,000 transcripts (DASL assay) for the clinical outcome of molecular signatures of HCC recurrence secondary to HCV chronic infection.

To help control for the variability encountered in other genomic studies, the present study performed microarray analysis utilizing FFPE stored liver tissue among a select group of patients, mostly male Caucasians with an average age of 55 years. Genes were successfully identified that are differentially expressed in tissue from post-LT patients with tumor recurrence and HCV compared with tissue from patients with no tumor recurrence and HCV after LT. This project has provided a subset of potential molecular signatures for future studies to explore their therapeutic and possible early marker potential for patients with HCC-R. This feasibility study has contributed information towards the better utilization of FFPE stored tissue and future large-scale multicenter expression analysis studies.

2. Methods

2.1. Patient Population

After IRB approval from the University of Washington, all patients receiving liver transplants in the University of Washington program from January 1, 2000 until July 1, 2007 were retrospectively screened. Each patient was categorized for the presence of HCC in the explanted liver. All patients with HCC were followed until death or for at least 3 years until July 1, 2010. Recipient factors recorded included age, gender, race, body mass index (BMI), etiology of liver disease, pretransplant alpha-fetoprotein (AFP) level, any prior treatment for HCV, any treatment for HCC prior to transplantation, date of HCC recurrence after liver transplantation, and date of death. Tumor characteristics recorded at the time of explant included HCC histology, size of largest tumor, number of tumors, presence of tumor in both lobes of liver, presence of macroinvasion of HCC, and calculated predicting cancer recurrence score (PCRS) [5]. Two groups, patients with HCC recurrence and patients without HCC recurrence following liver transplantation, were identified and matched for the presence of HCV and the same gender for further genomic study.

2.2. Tissue Preparation

All explanted livers were stored in refrigerators at 4°C until being sliced into 1 cm sections to evaluate for the presence of HCC. HCC tissue was fixed in formalin and embedded in paraffin. Any characteristic tumor nodule was stained for diagnosis. In preparation for RNA extraction, the FFPE serial sections were cut into sections with a thickness of 10 μm. The time period from storing the FFPE fixed tissue until preparation for RNA extraction was recorded. To differentiate nontumor tissue (adjacent cirrhotic tissue) and tumor tissue, hematoxylin and eosin (H+E) stained slides were microscopically evaluated by a pathologist. The tumor tissue was dissected from the nontumor tissue and sent for RNA extraction.

2.3. RNA Extraction, Purification, and Quality Assessment

Following deparaffinization, the tumor tissue was scraped from the slides into tubes and the RNA extracted using a modified protocol based on the Qiagen RNeasy FFPE kit (Qiagen, Valencia, CA, USA). Agilent 2100 bioanalysis (Agilent Technologies, Santa Clara, CA, USA) was performed to view the RNA quality. The cDNA synthesis was performed using the Qiagen Quantitect Reverse Transcription kit (Qiagen SpA, Milano, Italy). TaqMan quantitative polymerase chain reaction (qPCR) for the RPL-13A gene was performed using 200 ng of input RNA and 40 cycles of amplification as a quality assessment of the RNA; this is the RNA level quality control (QC) standard optimized by Illumina [16]. Similar studies by Abramovitiz [18] and Ton [19] had shown that a threshold value of ≤29 would signify sufficient RNA quality to give reproducible results in the cDNA-mediated annealing, selection, extension, and ligation (DASL) assay; this threshold was used in our study.

2.4. Gene Expression Profiling

The whole genome DASL assay developed by Illumina is highly specialized to detect 1.3- to 2-fold changes in intact and partially degraded RNA from FFPE samples [10]. As previously described and validated by Hoshida et al., DASL assay (6000 genes) allows highly reproducible expression profiles with RNA derived from FFPE samples from HCC resection specimens [17]. The Illumina WG-DASL HumanRef-12_V4 Expression BeadChip microarrays, which quantify approximately 29000 transcripts, were used for our study. In the DASL assay, total RNA was first converted into cDNA in a reverse transcription reaction using biotinylated primers and resulted in fluorescent-labeled PCR products that were annealed to the BeadChips. The fluorescence of hybridized transcripts was captured by laser confocal microscopy scanning using the Illumina BeadArray Reader system [20]. The DASL assay was performed by the Fred Hutchinson Cancer Research Center (FHCRC) Genomics Resource.

2.5. Data Analysis

The raw DASL array data was exported via genomeStudio version 1.0 (Illumina, Inc.) to JMP Genomics (SAS Institute, Inc., Cary, NC, USA) for quantile normalization, principal component analysis (PCA), and log transformation. The data analysis QC standard was performed on the sample and chip level to ensure uniformly high quality microarray data. Initially, the quality of each BeadArray was assessed by examining the percent present calls (defined as the percent of bead types having a detection call value <0.05), as well as plots of signal intensities for housekeeping, cy3_hyb, low_stringency_hyb, labeling, and biotin control bead types. In addition, log 2 transformation was performed, followed by quantile normalization and variance stabilizing transformation. For differential expression analyses, we followed the quintile with the greatest variance as previously described by Mittempergher and colleagues [21]. Genes differentially expressed between HCC-R and HCC-NR tumors were selected based on significance criteria of a false discovery rate less than 0.05 using analysis of variance (ANOVA) (≤0.05, with a fold change of 2). Heat maps were generated based on gene expression patterns through -means clustering and viewed using CLUSTER 3.0 and TreeView 1.45 (software at http://eisenlab.org/), respectively. The threshold applied to our statistical tests was influenced by the large patient-to-patient variability in gene expression, although heterogeneity in expression is expected in human samples. Survival analysis was conducted using Kaplan-Meier curves and compared using the log-rank test.

2.6. Functional Data Analysis

Network analysis of the differentially expressed genes was performed using Ingenuity Pathway Analysis software (IPA) (Ingenuity Systems, Inc., Redwood City, CA, USA). The top canonical pathways that represent differentially regulated genes in the tumor tissues from the HCC recurrence group were evaluated. The top-scoring network of interactions is presented for the concurrent downregulated and the concurrent upregulated gene sets. This software analyzes molecular data in the context of known biological response and regulatory networks as well as other higher-order response pathways. Ingenuity functional analysis identified biological functions and/or diseases that were most significantly enriched and generated values to determine the probability that each biological function assigned to that data set was due to chance alone.

2.7. Quantitative PCR (qPCR) Analysis

Quantitative PCR was performed on selected genes. The mRNA expression levels were quantified using a ViiA 7 real-time PCR System (Applied Biosystems, Life Technologies Corporation, Carlsbad, CA, USA) according to the manufacturer’s instructions. In brief, total RNA was extracted from tumor tissues. First-strand cDNA synthesis was performed using the Quantitect Reverse Transcription kit (Qiagen SpA, Milano, Italy). Quantitative PCR for respective custom primer-probe sets (Applied Biosystems) was performed using 9.0 ng of input RNA and 40 cycles of amplification using the ViiA 7 real-time PCR system (Applied Biosystems). RPL-13A was used as the endogenous control for mRNA levels. Each experiment was run in duplicate, including RPL-13A as the endogenous control and repeated 3 times. Relative quantification was performed using the delta method relative to RPL-13A as an internal control. The values were calculated with Student’s -test.

3. Results

3.1. Clinicopathological Characteristics of the Patient Population

Of the 768 patients undergoing transplantation at the University of Washington from January 1, 2000 until July 1, 2007, 93 patients (12%) had HCC due to HCV. Of these, 11 male patients developed recurrent tumors within 3 years following transplantation. From this group, 10 male patients with recurrence (HCC-R) were matched with second cohort group of 20 patients in whom HCC did not recur (HCC-NR) for HCV status, age and BMI. Following review of the quality control (QC) metrics (see Section 3.2), only 8 patients from each group were considered for further study. The clinical characteristics of the HCC-R group () and the HCC-NR () group revealed all patients to be male, and most of them were Caucasian. There were no significant differences between the two groups with respect to age, BMI, prior HCV or HCC treatment, number of tumors, and presence of tumors in both lobes (Table 1). The HCC-R group showed a higher AFP level, ( mg/dL; ) when compared with the HCC-NR group ( mg/dL). Likewise, the presence of macroinvasion and the PCRS was trending to be higher in the HCC recurrence group. The presence of poorly differentiated tumor was significantly higher () in the HCC-R group (75%) versus the HCC-NR group (0%). Kaplan-Meier survival curves revealed that the HCC-R group had significantly () lower 3-year survival (50%) compared with the HCC-NR group (100%) (Figure 1).

Table 1: Patient characteristics: clinical characteristics of the 16 patients from tumor recurrence and tumor nonrecurrence groups (mean ± SD or proportion).

Study group
HCC-R group
HCC-NR group
P-value

Number
8
8

Age, years

0.1

Male gender
100%
100%
1

Race

    Caucasian
75%
100%
0.2

    Asian
12.5%
0%

    Alaska native
12.5%
0%

BMI*

0.5

AFP level, mg/dL**

0.07

Cirrhosis, postnecrotic-type C
100%
100%
1

Prior HCV treatment
75%
50%
0.3

Prior HCC treatment
37.5%
0%
0.2

Hours to process***

0.06

Largest diameter Tumor, cm

0.1

Number of tumors

    1
25%
50%
0.2

    2
25%
12.5%

    (3–5)
0%
37.5%

    ≥6
50%
0.0%

Tumor in both lobes
37.5%
12.5%
0.6

Histology

    Poorly differentiated
75%
0%
<0.01

    Moderately differentiated
25%
75%

    Well differentiated
0%
25%

Macroinvasion
50%
0%
0.08

Body mass index; **alpha-fetoprotein level.

***Time specimens removed in surgical procedure until being fixed in formalin.

Figure 1: Kaplan-Meier survival curves comparing the high mortality rate in the tumor recurrence group with the mortality rate in the nonrecurrence group.

3.2. Quality Control Measure

The present genomic study began with a sample size of 30 (HCC-R () and HCC-NR ()). Following performance of QC metrics, 29 samples out of 30 (HCC-R 10/10 and HCC-NR 19/20) passed RNA level QC ( ≤ 29) and met the criteria for the WG-DASL array. Of these 29 samples following performance of the WG-DASL array, 16 (55%) samples (HCC-R: 8/10 (80%) and HCC-NR: 8/20 (40%)) met the data analysis QC standard and were selected for further comparative gene expression analysis (Figure 2). Overall, 14 (47%) of the samples were excluded from the final analysis.

Figure 2: Experimental strategy for comparative transcriptome analysis of HCV-associated HCC with or without recurrence, using FFPE tissue to determine molecular signatures of HCC recurrence.

3.3. Transcriptome Analysis Defines Differentially Regulated Genes, and Major Canonical Pathways Are Associated with Recurrence of HCC

The comparative transcriptome analyses aimed at identifying molecular signatures representative of HCC recurrence were carried out as described in Figure 2. Hierarchical clustering by gene expression segregated all HCC-R from HCC-NR. A total of 194 genes were identified to be differentially expressed, with 151 genes upregulated and 43 downregulated in HCC-R (Figure 3).

Figure 3: Hierarchical clustering for 194 genes differentially expressed in HCC-R (FDR corrected ≤0.05). For visual comparison, genes differentially expressed in HCC-R and HCC-NR were clustered by the TreeView program. The red color represents genes upregulated, and the green color represents genes downregulated. Saturation was set at ±2-fold change.

The biological significance of the altered gene expression pattern described above was investigated by classifying the associated gene within the context of biologically relevant functions using IPA. HCC with recurrence exhibited enrichment of upregulated genes mapping to signaling or disease pathways associated with cell cycle regulators (CDKN2B, E2F2, E2F5, GNL3, HDAC2, MDM2, MYC, and PA2G4), including the genes that encode the proteins that control molecular mechanisms of cancer (FANCD2, FZD3, PLCB1, and PMAIP1). (; Table 2). Gene-encoding proteins implied that nucleo-cytoplasmic transport was also overexpressed (KPNA2, KPNB1, RANBP1, and RCC1). The presence of downregulated pathways reflects the concurrent downregulation of directly related genes categorized in nicotine degradation, complement system, hepatic cholestasis, and catecholamine biosynthesis (; Table 2). FXR/RXR activation, which is associated with hepatoprotection [53], was downregulated in recurrent HCC tumors.

Table 2: Top canonical pathways: significant molecular pathways regulated in HCC-R.

Up-regulated in HCC-R
P value
Genes

Cell cycle: G1/S checkpoint regulation

CDKN2B, E2F2, E2F5, GNL3, HDAC2, MDM2, MYC, PA2G4

RAN signaling

KPNA2, KPNB1, RANBP1, RCC1

Chronic myeloid leukemia signaling

E2F2, E2F5, HDAC2, MDM2, MYC, PA2G4

Molecular mechanisms of cancer

CDKN2B, E2F2, E2F5, FANCD2, FZD3, MDM2, MYC, PA2G4, PLCB1, PMAIP1

Cyclins and cell cycle regulation

CDKN2B, E2F2, E2F5, HDAC2, PA2G4

Downregulated in HCC-R
-value
Genes

Nicotine degradation II

FMO3, FMO4, INMT

FXR/RXR activation

FETUB, SLC10A1, SLC22A7

Complement system

C7, MASP1

Hepatic cholestasis

GCGR, SLC10A1, SLC22A7

Catecholamine biosynthesis

DBH

3.4. Network Analysis Defines Top-Scoring Upregulated Genes Associated with HCC Recurrence

Interestingly, when network analysis was performed, which allows exploration of the biological relationship between any two genes, it was seen that the largest number of genes with higher expression in HCC-R tumor tissue was dominated by major transcriptional regulators. These genes are associated with cellular malignancies by controlling cell cycle progression, cell growth and proliferation, cell-to-cell signaling, and cell survival and death. Additionally, a novel set of highly significant genes associated with tumor recurrence was identified (Figure 4(a)).

(a)

(b)

Figure 4: Network analysis of directly related genes which are upregulated and downregulated in HCC-R tumor tissues compared with HCC-NR tumor tissues. (a) The upregulated network includes genes dominated by regulators of cell cycle progression, cell growth and proliferation, DNA recombination, and cell-to-cell signaling genes. These are represented as respective shapes based on functional property in HCC-R tumor tissues. (b) The downregulated network includes genes reflecting the concurrent downregulation of genes related to innate immune response, cell-to-cell signaling, cell morphology, and cellular metabolism genes. These are represented as respective shapes based on their functional properties in HCC-R tumor tissues. The gray outlined nodules without expression are linker genes that are not altered in HCC-R but are statistically enriched for interaction with the altered genes.

3.4.1. Activation of the Major Transcriptional Regulators

The expression data revealed an increased (≥3-fold) expression of transcriptional regulators in HCC-R tumor tissues (MYC (3.72), CTNNB1 (3.19), and MDM2 (4.27)). The respective -score for MYC is 2.080 (), as obtained from the IPA transcription factor (TF) analysis, which predicts the activation of TFs based on the expression levels of their known targets. -catenin (CTNNB1), another major transcriptional regulator of specific oncogenes, was also overexpressed in the HCC-R tissue. Many of the downstream targets of the transcription factors MYC and CTNNB1 were upregulated in a fashion consistent with their increased abundance.

3.4.2. Genes Associated with Cellular Malignancies

A number of key genes related to cellular malignancies were upregulated in HCC-R tumors, specifically, HMGA1 (4.16), SPP1 (3.90), GNL3 (3.86), and PPARG (3.84). Two major kinases, including the cyclin-dependent kinase inhibitor 2A (CDKN2A) and cyclin-dependent kinase inhibitor 2B (CDKN2B), showed significantly increased expression (fold changes 3.52 and 3.99, resp.) and are known to be associated with abnormal cell growth. A highly significant up-regulation of a large number of genes () implicated cell death and survival. The network included genes with fold change ≥4: RIPK2 (4.35) and NCKAP1 (4.02). Additionally, MDK (3.99), NME1 (4.82), and PA2G4 (4.21) showed increased expression in the tumor tissues from the patients with recurrent HCC.

3.4.3. Novel Highly Significant Genes Associated with Tumor Recurrence

A set of genes including PRPF38A, RIOK3, QSER1, PSMC3IP, ATAD3B, MGC12982, and C20ORF27 were significantly overexpressed in HCC-R tumor tissue with fold changes ≥4 ( = 0.001 to 0.0001). Most of these genes are known to be involved in various regulatory mechanisms such as cell cycle regulation and DNA replication and recombination mechanism (Table 3); however, the mechanism of action is not well characterized. QSER1, MGC12982 and C20ORF27 are completely uncharacterized. Genes such as MCM7 (4.78), DFFA (5.24), and PRPF38A (4.28) are engaged in gene regulation while E2F5 (6.60), RPS6KA3 (4.77), and YWHAZ (5.18) are associated with cell proliferation. Anti-apoptotic genes such as RFFL (5.32), EIF3H (4.27), and HDAC2 (4.05) were overexpressed. RIOK3 (5.00; ) and RCC1 (5.32; ) are the top significant genes associated with cytoskeletal architecture and the cellular transport mechanism.

Table 3: Highly significant upregulated genes in HCC-R tumors with known and unknown functions (references are based on PUBMED search).

ILMN Gene ID no.
Genes
P
FC
Known cellular function (Ref.)
Role in HCC (Ref.)
HCC-R secondary to HCV (Ref.)

1667213
DNA fragmentation factor, 45 kDa, alpha polypeptide (DFFA)
0.0011
5.24
(1) Tumor suppressor gene [22].

(2) Protein encoded by DFFA gene known as inhibitor of caspase-activated DNase (ICAD), involved in DNA fragmentation during apoptosis [23].
Unknown
Unknown

1664111
Ring finger and FYVE-like domain containing E3 ubiquitin protein ligase (RFFL)
0.0005
4.72
(1) RFFL is a member of caspase 8/10-associated RING proteins (CARPs) are a recently described family of protein ubiquitin ligases that interact with and negatively regulate death receptor-mediated apoptosis [24].
Unknown
Unknown

1675626
PRP38 pre-mRNA processing factor 38 (yeast) domain containing (PRPF38A)
0.0005
4.28
(1) PRP38 pre-mRNA processing factor 38 (yeast) domain containing A protein involved in pre-mRNA splicing [25].
Unknown
Unknown

1767747
Histone deacetylase 2 (HDAC2)

4.05
(1) This gene product belongs to the histone deacetylase family known to play an important role in transcriptional regulation, cell cycle progression and developmental events.

(2) Role in antiapoptosis in tumor cells [26].
[27]
Unknown

2147863
Glutamine and serine rich 1 (QSER1)

3.98
Unknown
Unknown
[28]

2354211
Regulator of chromosome condensation 1 (RCC1)
0.0002
5.32
(1) Regulator of chromosomal condensation [29].
Unknown
Unknown

2404135
RIO kinase 3 (yeast) (RIOK3)
<0.0001
5.00
(1) The specific function of this gene has not yet been determined (NCBI).

(2) Recent study identified Rio kinase 3 (RIOK3) as an amplified gene that alters cytoskeletal architecture as well as promotes pancreatic ductal cell migration and invasion [30].
Unknown

1770822
Ribosomal protein S6 kinase, 90 kDa, polypeptide 3 (RPS6KA3)
0.0001
4.77
(1) Known to regulating cell growth and differentiation [31].

(2) Known to play a role in the cellular antiviral response [32].
[33]
Unknown

1663195
Minichromosome maintenance complex component 7 (MCM7)
0.0031
4.78
(1) Respective protein involved in regulating DNA replication [34].

(2) MCM7 overexpression and amplification and several human malignancies [35].
[36]
Unknown

1801928
Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide (YWHAZ)
0.0007
5.18
(1) The encoded protein interacts with IRS1 protein, suggesting a role in regulating insulin sensitivity [37].

(2) Known to promote malignancies in many cell types [38].
[39]
Unknown

2057981
Zinc finger, C2HC-type containing 1A (FAM164A)
0.0008
4.90
Unknown
Unknown
Unknown

1782551
E2F transcription factor 5, p130-binding (E2F5)
<0.0001
6.60
(1) A transcription factor regulates cell proliferation and differentiation [40].

(2) Over-expressed in several human malignancies [41].
[37]
Unknown

3307266
PSMC3 interacting protein (PSMC3IP)
0.0006
5.23
(1) This gene encodes a protein that functions in meiotic recombination [42].
Unknown
Unknown

2131936
ATPase family, AAA domain containing 3B (ATAD3B)
0.001
4.55
(1) ATAD3B is mitochondrial membrane proteins that contribute to the stabilization of large mitochondrial DNA (mtDNA) protein complexes called nucleoids [43].
Unknown
Unknown

1683660
Eukaryotic translation initiation factor 3, subunit H (EIF3H)
0.0016
4.27
(1) At the cellular level, EIF3H overexpression increases proliferation, growth and survival [44].

(2) Inhibit induction of apoptosis [44].
[45]
Unknown

3309468
Hypothetical protein MGC12982 (MGC12982)
0.0006
4.09
Unknown
Unknown
Unknown

1697363
Chromosome 20 open reading frame 27 (C20ORF27)
0.0003
4.15
Unknown
Unknown
Unknown

= value; FC = fold change.

3.5. Network Analysis Defines Top-Scoring Downregulated Genes Associated with HCC Recurrence

Underexpressed genes were categorized in regulating innate immune response, cell-to-cell signaling and interaction, and the inflammatory response (Figure 4(b)). Two major transcriptional regulators, the hepatocyte nuclear factor 4, alpha (HNF4A) and ubiquitin C (UBC) genes, had no predictive expression value in HCC-R tumor tissues.

3.5.1. Inactivation of the Major Transcriptional Regulators

The downregulated network includes two major transcriptional regulators, HNF4A and UBC, with no predictive expression values. However, most of the downstream targets of the transcription factors HNF4A and UBC were downregulated in a fashion consistent with their decreased abundance.

3.5.2. Genes Associated with Innate Immunity, Cell-to-Cell Signaling and Interaction

The genes that displayed the most dramatic downregulation were MASP1 (−5.11), C7 (−5.06), DBH (−4.66), and FBLN5 (−4.48). However, such genes as IFI27 (−4.38), SLC22A7 (−3.97), IFIT1 (−3.91), TGFB3 (−3.71), and IFN- induced (−3.72) also had stronger suppression at the transcript level in HCC-R tumor tissues.

3.5.3. Genes Associated with Inflammatory Response

The inflammatory response-related genes (CCL14 (−4.52), LEP (−3.73) and PTGDS (−3.84)) showed decreased expression in the recurrence tumor tissues.

3.5.4. Other Significant Genes

Several genes, SLC10A1, GCGR, FMO3, and INMT, were previously known to be downregulated in HCV-induced HCC [28, 54]. These genes showed similar expression changes with the fold change ≥3.5. Additionally, a set of genes, including AFM, DYNLRB2, FDX1, and SHBG, were significantly downregulated in HCC-R tumor tissues with fold changes ≥4 involved in various small molecule biochemistry and molecular transport mechanisms. These genes were previously known to be suppressed in HCC [55–57].

3.6. Validation of Gene Expression in HCC-R and HCC-NR Tumor Tissues

Total RNA from the same two groups of patients, the HCC-R group and the HCC-NR group, was used for qPCR analysis to validate microarray expression data. The following genes were selected for validation from a set of novel highly significant upregulated genes from our expression dataset (Table 3): DFFA, RIOK3, E2F5, EIF3H, YWHAZ, QSER1, RPS6KA3, PRPF38A, MCM7, and C20ORF27. Additionally, we selected a few genes (CTNNB1, PPARG, HIF1A, HMGA1, MYC, and CDKN2A with values <0.05; Table 4) that are a hallmark of HCC in the literature [47, 49, 50, 58, 59]. The expression levels determined by qPCR were comparable to the microarray data in HCC-R versus HCC-NR tissue (Figures 5 and 6).

Table 4: Gene list known to be deregulated in HCC selected from our microarray data for cross validation as positive control (references are based on PUBMED search).

Genes
FC
Role in HCC (Ref.)

v-myc myelocytomatosis viral oncogene homolog (avian) (MYC)
3.72
[46]

Catenin (cadherin-associated protein), beta 1, 88 kDa (CTNNB1)
3.19
[47]

Hypoxia inducible factor 1, alpha subunit (basic helix-loop-helix transcription

factor) (HIF1A)

3.06
[48]

High mobility group AT-hook 1 (HMGA1)
4.16
[49]

Peroxisome proliferator-activated receptor gamma (PPARG)
3.84
[50]

Cyclin-dependent kinase inhibitor 2A (CDKN2A)
4.52
[51, 52]

FC = fold change (based on our expression data).

(a)

(b)

Figure 5: Validation of microarray data by qPCR. The values were calculated with Student’s -test (*).

(a)

(b)

Figure 6: Validation of microarray data by qPCR. The values were calculated with Student’s -test (*** < 0.05; ** < 0.07; * < 0.25).

4. Discussion

Recurrent HCC following surgical treatment or liver transplantation continues to be a serious health problem [60]. New therapeutic methods need to be developed. The aim of genomic analysis is to enable development of these new treatment options. Presently, the molecular mechanisms involved in HCV-infected individuals who develop HCC recurrence are largely unknown.

To find these molecular mechanisms, several studies using many samples will need to be conducted. This present study has revealed FFPE tissue to be a good source of such study material; the use of FFPE tissue will greatly expand the number of samples available for study.

The present study revealed both known and previously unknown molecular patterns associated with HCC recurrence; however, this is a proof of concept study on a small subset of HCC patients. All patients were chronically infected with HCV, most were Caucasian, and the average age was 55 years. Tumor characteristics varied in terms of histology, vascular invasion, tumor grade, and number and size of tumors. The small sample size and wide variety of tumor characteristics do not allow for any strong conclusions to be drawn from this study; the highly differentiated levels of gene expression are interesting to review, however.

In this study, the use of FFPE-stored tissues for whole transcriptome analysis provided both previously known and novel insights into HCV-associated HCC with or without recurrence. A set of 194 genes differentially expressed was observed in two groups of patients with HCV and either experiencing or not experiencing HCC recurrence. The present study identified and quantified 6 genes that were highly overexpressed in our HCC-R samples (CTNNB1, PPARG, HIF1A, HMGA1, MYC, and CDKN2A) that other investigators have determined to be hallmarks of HCC [46, 47, 49, 50]. Additionally, the present study identified, for the first time, a set of highly significant upregulated genes in the HCC-R tumor tissues. These genes are involved with gene regulation (DFFA, MCM7, and PRPF38A), cell proliferation (E2F5, RPS6KA3, and YWHAZ) cytoskeletal architecture (RIOK3), anti-apoptosis (EIF3H, RFFL, and HDAC2) and uncharacterized functions (QSER1, MGC12982, and C20ORF27). The increased expression of anti-apoptotic genes, including EIF3H, RFFL and HDAC2, may favorably assist uncontrolled cell growth and proliferation to enhance tumor growth in patients with HCC-R. Previous studies have shown that RFFL and HDAC2 are known to also promote tumor formation by inhibiting the apoptosis process [26, 27, 61], whereas EIF3H is known to increase cell proliferation, growth and survival and inhibit the apoptosis process [44]. We observed an increased expression of the YWHAZ gene, which is known as a potential metastasis factor with the antiapoptotic property of promoting cellular malignancy [38, 62]. Overall data suggests that HCC recurrence is associated with increased apoptosis inhibition, and uncontrolled cell growth leads to a tumor metastasis microenvironment that favors tumor recurrence in patients with HCC-R. Of note is that a major hallmark of an aggressive HCC is its ability to metastasize [63]; the recurrence of HCC supports the metastatic phenomenon. The observed increased expression of RCC1 is one of the most important members of the RAN signaling pathway, which is involved in the nucleocytoplasmic transport of macromolecules [64, 65]. Recent findings have shown that silencing RAN expression could induce more apoptosis in cancer cells, and therefore is a promising cancer therapeutic target [66]. This suggests a novel link between the elevated RAN signaling pathway in HCC recurrence and a potentially important role for nucleo-cytoplasmic transport mechanisms of RCC1 during HCC progression. The increased expression of RIOK3 is known to alter the cytoskeletal architecture, as well as promoting pancreatic ductal cell migration and invasion [30]. The increased expression of RIOK3 has been observed in metastatic head and neck cancers compared with nonrecurrent tumors [67]. These observations raise an interesting prospect that similarly, increased expression of RIOK3 may contribute to cytoskeletal architecture alteration to influence cell migration and tumor invasion in HCC patients having tumor recurrence. Thus, the findings of the present study further point toward new avenues of research aimed at evaluating the impact of anti-apoptosis, cytoskeletal architecture alteration, and RAN signaling on HCC recurrence.

A limitation of this study was the use of only 50% of the FFPE tissue with microarray analysis and the low number of tissues used for the final analysis. Further investigation with larger cohort is warranted. Comparing specific subgroups of different liver diseases, races/ethnicities, ages, and tumor characteristics could reveal clinical implications that could potentially aid in patient selection for liver transplantation. Future studies with recent advanced technology such Next Generation RNA Sequencing (RNA-seq) might offer greater potential for the use of FFPE samples, with a tremendous increase in the number of samples to study. RNA-seq, a recently developed approach to transcriptome profiling that uses deep-sequencing technologies [68], could offer a greater opportunity to use FFPE samples [8, 69] to bring gene expression results into the clinical treatment of HCC.

In conclusion, this pilot expression profiling study using FFPE tissue has shown that stored FFPE tissue is a vital resource and has identified molecular patterns for HCC-R tumor tissue consistent with prior studies. We also identified a set of genes not previously reported to be associated with HCC-R. All of these genes may be potential targets for future therapeutic interventions.

Abbreviations

HCV: | Hepatitis C virus
HCC: | Hepatocellular carcinoma
BMI: | Body mass index
FFPE: | Formalin-fixed paraffin-embedded
IRB: | Institutional Review Board
AFP: | Alpha-fetoprotein
PCRS: | Predicting cancer recurrence score
DASL: | cDNA-mediated annealing, selection, extension, and ligation
QC: | Quality control
IPA: | Ingenuity pathway analysis
PCA: | Principal component analysis
ATAD3A/ATAD3B: | ATPase family, AAA domain containing 3A
DASL: | cDNA-mediated annealing, selection, extension, and ligation
CDK4: | Cyclin-dependent kinase 4
CDKN2A: | Cyclin-dependent kinase inhibitor 2A
CDKN2B: | Cyclin-dependent kinase inhibitor 2B (p15, inhibits CDK4)
DFFA: | DNA fragmentation factor, 45 kDa, alpha polypeptide
E2F2: | E2F transcription factor 2
E2F5: | E2F transcription factor 5, p130-binding
EIF3H: | Eukaryotic translation initiation factor 3, subunit H
ENO1: | Enolase 1, (alpha)
EZH2: | Enhancer of zeste homolog 2 (Drosophila)
FDR: | False discovery rate
G6PC3: | Glucose 6 phosphatase, catalytic, 3
GNL3: | Guanine nucleotide binding protein-like 3 (nucleolar)
HCC-NR: | Hepatocellular carcinoma nonrecurrence
HCC-R: | Hepatocellular carcinoma recurrence
HDAC2: | Histone deacetylase 2
IL7: | Interleukin 7
LT: | Liver transplantation
MCM7: | Minichromosome maintenance complex component 7
MDK: | Midkine (neurite growth-promoting factor 2)
MDM2: | Mdm2, p53 E3 ubiquitin protein ligase homolog (mouse)
MYC: | v-myc myelocytomatosis viral oncogene homolog (avian)
NCKAP1: | NCK-associated protein 1
NME1: | NME/NM23 nucleoside diphosphate kinase 1
PA2G4: | Proliferation-associated 2G4, 38 kDa
POLA2: | Polymerase (DNA directed), alpha 2, accessory subunit
PPARG: | Peroxisome proliferator-activated receptor gamma
PRPF38A: | PRP38 pre-mRNA processing factor 38 (yeast) domain containing A
PSMC3IP: | PSMC3 interacting protein
qPCR: | Quantitative real-time polymerase chain reaction
QSER1: | Glutamine and serine rich 1
RCC1: | Regulator of chromosome condensation 1
RFFL: | Ring finger and FYVE-like domain containing E3 ubiquitin protein ligase
RIOK3: | RIO kinase 3 (yeast)
RIPK2: | Receptor-interacting serine-threonine kinase 2
RPS6KA3: | Ribosomal protein S6 kinase, 90 kDa, polypeptide 3
RPSA: | Ribosomal protein SA
RUVBL2: | RuvB-like 2 (E. coli)
SPP1: | Secreted phosphoprotein 1
SSBP1: | Single-stranded DNA binding protein 1, mitochondrial
TP53: | Tumor protein p53
YWHAZ: | Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide
ZC2HC1A: | Zinc finger, C2HC-type containing 1A
ADHFE1: | Alcohol dehydrogenase, iron containing, 1
AFM: | Afamin
AHR: | Aryl hydrocarbon receptor
ANO1: | Anoctamin 1, calcium activated chloride channel
APBB3: | Amyloid beta (A4) precursor protein-binding, family B, member 3
APP: | Amyloid beta (A4) precursor protein
BAG6: | BCL2-associated athanogene 6
BDH1: | 3-Hydroxybutyrate dehydrogenase, type 1
BRCA1: | Breast cancer 1, early onset
C7: | Complement component 7
CCL14: | Chemokine (C-C motif) ligand 14
CEBPB: | CCAAT/enhancer binding protein (C/EBP), beta
CLU: | Clusterin
CRIP2: | Cysteine-rich protein 2
DBH: | Dopamine beta-hydroxylase (dopamine beta-monooxygenase)
DLG4: | Discs, large homolog 4 (Drosophila)
DPT: | Dermatopontin
DYNLRB2: | Dynein, light chain, roadblock-type 2
ESR1: | Estrogen receptor 1
FBLN5: | Fibulin 5
FDX1: | Ferredoxin 1
FETUB: | Fetuin B
FMO3: | Flavin containing monooxygenase 3
FMO4: | Flavin containing monooxygenase 4
GCGR: | Glucagon receptor
GSK3B: | Glycogen synthase kinase 3 beta
HNF4A: | Hepatocyte nuclear factor 4, alpha
HSD11B1: | Hydroxysteroid (11-beta) dehydrogenase 1
IFIT1: | Interferon-induced protein with tetratricopeptide repeats 1
INMT: | Indolethylamine N-methyltransferase
KLF11: | Kruppel-like factor 11
LEP: | Leptin
LIPC: | Lipase, hepatic
MASP1: | Mannan-binding lectin serine peptidase 1 (C4/C2 activating component of Ra-reactive factor)
NLGN2: | Neuroligin 2
OLFML1: | Olfactomedin-like 1
PBLD: | Phenazine biosynthesis-like protein domain containing
PGRMC1: | Progesterone receptor membrane component 1
PTGDS: | Prostaglandin D2 synthase 21 kDa (brain)
SHBG: | Sex hormone-binding globulin
SLC10A1: | Solute carrier family 10 (sodium/bile acid cotransporter family), member 1
SLC22A7: | Solute carrier family 22 (organic anion transporter), member 7
TGFB3: | Transforming growth factor, beta 3
TGM2: | Transglutaminase 2 (C polypeptide, protein-glutamine-gamma-glutamyltransferase)
TK1: | Thymidine kinase 1, soluble
TRAF6: | TNF receptor-associated factor 6, E3 ubiquitin protein ligase
UBC: | Ubiquitin C.

Conflict of Interests

No conflict of interests was disclosed.

References

A. M. di Bisceglie, A. C. Lyra, M. Schwartz et al., “Hepatitis C-related hepatocellular carcinoma in the United States: influence of ethnic status,” American Journal of Gastroenterology, vol. 98, no. 9, pp. 2060–2063, 2003.

K. Kiyosawa, T. Umemura, T. Ichijo et al., “Hepatocellular carcinoma: recent trends in Japan,” Gastroenterology, vol. 127, no. 5, supplement 1, pp. S17–S26, 2004.

J. M. Llovet, A. Burroughs, and J. Bruix, “Hepatocellular carcinoma,” The Lancet, vol. 362, no. 9399, pp. 1907–1917, 2003.

V. Mazzaferro, E. Regalia, R. Doci et al., “Liver transplantation for the treatment of small hepatocellular carcinomas in patients with cirrhosis,” The New England Journal of Medicine, vol. 334, no. 11, pp. 693–699, 1996.

E. Y. Chan, A. M. Larson, O. K. Fix et al., “Identifying risk for recurrent hepatocellular carcinoma after liver transplantation: implications for surveillance studies and new adjuvant therapies,” Liver Transplantation, vol. 14, no. 7, pp. 956–965, 2008.

S. von Ahlfen, A. Missel, K. Bendrat, and M. Schlumpberger, “Determinants of RNA quality from FFPE samples,” PLoS One, vol. 2, no. 12, Article ID e1261, 2007.

S. M. Farragher, A. Tanney, R. D. Kennedy, and D. Paul Harkin, “RNA expression analysis from formalin fixed paraffin embedded tissues,” Histochemistry and Cell Biology, vol. 130, no. 3, pp. 435–445, 2008.

L. Turner, J. D. Heath, and N. Kurn, “Gene expression profiling of RNA extracted from FFPE tissues: NuGEN technologies' whole-transcriptome amplification system,” Methods in Molecular Biology, vol. 724, pp. 269–280, 2011.

R. A. Coudry, S. I. Meireles, R. Stoyanova et al., “Successful application of microarray technology to microdissected formalin-fixed, paraffin-embedded tissue,” Journal of Molecular Diagnostics, vol. 9, no. 1, pp. 70–79, 2007.

C. April, B. Klotzle, T. Royce et al., “Whole-genome gene expression profiling of formalin-fixed, paraffin-embedded tissue samples,” PLoS One, vol. 4, no. 12, Article ID e8162, 2009.

N. Iizuka, M. Oka, H. Yamada-Okabe et al., “Oligonucleotide microarray for prediction of early intrahepatic recurrence of hepatocellular carcinoma after curative resection,” The Lancet, vol. 361, no. 9361, pp. 923–929, 2003.

Y. Kurokawa, R. Matoba, I. Takemasa et al., “Molecular-based prediction of early recurrence in hepatocellular carcinoma,” Journal of Hepatology, vol. 41, no. 2, pp. 284–291, 2004.

V. R. Mas, R. A. Fisher, K. J. Archer et al., “Genes associated with progression and recurrence of hepatocellular carcinoma in hepatitis C patients waiting and undergoing liver transplantation: preliminary results,” Transplantation, vol. 83, no. 7, pp. 973–981, 2007.

C. L. Lai, V. Ratziu, M.-F. Yuen, and T. Poynard, “Viral hepatitis B,” The Lancet, vol. 362, no. 9401, pp. 2089–2094, 2003.

T. Poynard, M.-F. Yuen, V. Ratziu, and C. Lung Lai, “Viral hepatitis C,” The Lancet, vol. 362, no. 9401, pp. 2095–2100, 2003.

J.-B. Fan, J. M. Yeakley, M. Bibikova et al., “A versatile assay for high-throughput gene expression profiling on universisal array matrices,” Genome Research, vol. 14, no. 5, pp. 878–885, 2004.

Y. Hoshida, A. Villanueva, M. Kobayashi, et al., “Gene expression in fixed tissues and outcome in hepatocellular carcinoma,” The New England Journal of Medicine, vol. 359, no. 19, pp. 1995–2004, 2008.

M. Abramovitz, M. Ordanic-Kodani, Y. Wang et al., “Optimization of RNA extraction from FFPE tissues for expression profiling in the DASL assay,” BioTechniques, vol. 44, no. 3, pp. 417–423, 2008.

C. C. Ton, N. Vartanian, X. Chai et al., “Gene expression array testing of FFPE archival breast tumor samples: an optimized protocol for WG-DASL sample preparation,” Breast Cancer Research and Treatment, vol. 125, no. 3, pp. 879–883, 2011.

M. Ravo, M. Mutarelli, L. Ferraro et al., “Quantitative expression profiling of highly degraded RNA from formalin-fixed, paraffin-embedded breast tumor biopsies by oligonucleotide microarrays,” Laboratory Investigation, vol. 88, no. 4, pp. 430–440, 2008.

L. Mittempergher, J. J. de Ronde, M. Nieuwland et al., “Gene expression profiles from formalin fixed paraffin embedded breast cancer tissue are largely comparable to fresh frozen matched tissue,” PLoS One, vol. 6, no. 2, Article ID e17163, 2011.

Y. Z. Chen, E. Soeda, H. W. Yang et al., “Homozygous deletion in a neuroblastoma cell line defined by a high-density STS map spanning human chromosome band Ip36,” Genes Chromosomes and Cancer, vol. 31, no. 4, pp. 326–332, 2001.

X. Liu, H. Zou, C. Slaughter, and X. Wang, “DFF, a heterodimeric protein that functions downstream of caspase-3 to trigger DNA fragmentation during apoptosis,” Cell, vol. 89, no. 2, pp. 175–184, 1997.

W. Yang, L. M. Rozan, E. R. McDonald III et al., “CARPs are ubiquitin ligases that promote MDM2-independent p53 and phospho-p53ser20 degradation,” Journal of Biological Chemistry, vol. 282, no. 5, pp. 3273–3281, 2007.

A. Vinayagam, U. Stelzl, R. Foulle et al., “A directed protein interaction network for investigating intracellular signal transduction,” Science Signaling, vol. 4, no. 189, p. rs8, 2011.

W.-W. Lei, K.-H. Zhang, X.-C. Pan et al., “Histone deacetylase 1 and 2 differentially regulate apoptosis by opposing effects on extracellular signal-regulated kinase 1/2,” Cell Death and Disease, vol. 1, no. 5, article e44, 2010.

J. H. Noh, K. H. Jung, J. K. Kim et al., “Aberrant regulation of HDAC2 mediates proliferation of hepatocellular carcinoma cells by deregulating expression of G1/S cell cycle proteins,” PLoS One, vol. 6, no. 11, Article ID e28103, 2011.

E. Wurmbach, Y.-B. Chen, G. Khitrov et al., “Genome-wide molecular profiles of HCV-induced dysplasia and hepatocellular carcinoma,” Hepatology, vol. 45, no. 4, pp. 938–947, 2007.

M. E. Nemergut, C. A. Mizzen, T. Stukenberg, C. D. Allis, and I. G. Macara, “Chromatin docking and exchange activity enhancement of RCC1 by histones H2A and H2B,” Science, vol. 292, no. 5521, pp. 1540–1543, 2001.

A. C. Kimmelman, A. F. Hezel, A. J. Aguirre, et al., “Genomic alterations link Rho family of GTPases to the highly invasive phenotype of pancreas cancer,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 49, pp. 19372–19377, 2008.

J. Xing, D. D. Ginty, and M. E. Greenberg, “Coupling of the RAS-MAPK pathway to gene activation by RSK2, a growth factor-regulated CREB kinase,” Science, vol. 273, no. 5277, pp. 959–963, 1996.

S. Kakugawa, M. Shimojima, H. Goto et al., “Mitogen-activated protein kinase-activated kinase RSK2 plays a role in innate immune responses to influenza virus infection,” Journal of Virology, vol. 83, no. 6, pp. 2510–2517, 2009.

C. Guichard, G. Amaddeo, S. Imbeaud et al., “Integrated analysis of somatic mutations and focal copy-number changes identifies key genes and pathways in hepatocellular carcinoma,” Nature Genetics, vol. 44, no. 6, pp. 694–698, 2012.

L. Homesley, M. Lei, Y. Kawasaki, S. Sawyer, T. Christensen, and B. K. Tye, “Mcm10 and the MCM2-7 complex interact to initiate DNA synthesis and to release replication factors from origins,” Genes and Development, vol. 14, no. 8, pp. 913–926, 2000.

K. A. Honeycutt, Z. Chen, M. I. Koster et al., “Deregulated minichromosomal maintenance protein MCM7 contributes to oncogene driven tumorigenesis,” Oncogene, vol. 25, no. 29, pp. 4027–4032, 2006.

J. Liu, S.-K. Yao, and F. Yin, “Expression and significance of minichromosome maintenance protein7 in hepatocellular carcinoma,” Zhonghua Gan Zang Bing Za Zhi, vol. 16, no. 1, pp. 55–56, 2008.

R. W. Furlanetto, B. R. Dey, W. Lopaczynski, and S. P. Nissley, “14-3-3 proteins interact with the insulin-like growth factor receptor but not the insulin receptor,” Biochemical Journal, vol. 327, part 3, pp. 765–771, 1997.

M. Niemantsverdriet, K. Wagner, M. Visser, and C. Backendorf, “Cellular functions of 14-3-3ζ in apoptosis and cell adhesion emphasize its oncogenic character,” Oncogene, vol. 27, no. 9, pp. 1315–1319, 2008.

J. E. Choi, W. Hur, C. K. Jung et al., “Silencing of 14-3-3ζ over-expression in hepatocellular carcinoma inhibits tumor growth and enhances chemosensitivity to cis-diammined dichloridoplatium,” Cancer Letters, vol. 303, no. 2, pp. 99–107, 2011.

C. Sardet, M. Vidal, D. Cobrinik et al., “E2F-4 and E2F-5, two members of the E2F family, are expressed in the early phases of the cell cycle,” Proceedings of the National Academy of Sciences of the United States of America, vol. 92, no. 6, pp. 2403–2407, 1995.

N. Kothandaraman, V. B. Bajic, P. N. K. Brendan et al., “E2F5 status significantly improves malignancy diagnosis of epithelial ovarian cancer,” BMC Cancer, vol. 10, article 64, 2010.

G. V. Petukhova, R. J. Pezza, F. Vanevski, M. Ploquin, J.-Y. Masson, and R. D. Camerini-Otero, “The Hop2 and Mnd1 proteins act in concert with Rad51 and Dmc1 in meiotic recombination,” Nature Structural and Molecular Biology, vol. 12, no. 5, pp. 449–453, 2005.

J. He, C.-C. Mao, A. Reyes et al., “The AAA+ protein ATAD3 has displacement loop binding properties and is involved in mitochondrial nucleoid organization,” Journal of Cell Biology, vol. 176, no. 2, pp. 141–146, 2007.

L. Zhang, Z. Smit-McBride, X. Pan, J. Rheinhardt, and J. W. B. Hershey, “An oncogenic role for the phosphorylated h-subunit of human translation initiation factor eIF3,” Journal of Biological Chemistry, vol. 283, no. 35, pp. 24047–24060, 2008.

H. Okamoto, K. Yasui, C. Zhao, S. Arii, and J. Inazawa, “PTK2 and EIF3S3 genes may be amplification targets at 8q23-q24 and are associated with large hepatocellular carcinomas,” Hepatology, vol. 38, no. 5, pp. 1242–1249, 2003.

M. Saegusa, Y. Takano, H. Kishimoto, G. Wakabayashi, K. Nohga, and M. Okudaira, “Comparative analysis of p53 and c-myc expression and cell proliferation in human hepatocellular carcinomas—an enhanced immunohistochemical approach,” Journal of Cancer Research and Clinical Oncology, vol. 119, no. 12, pp. 737–744, 1993.

Y. Hoshida, S. M. B. Nijman, M. Kobayashi et al., “Integrative transcriptome analysis reveals common molecular subclasses of human hepatocellular carcinoma,” Cancer Research, vol. 69, no. 18, pp. 7385–7392, 2009.

F. Simon, M. Bockhorn, C. Praha et al., “Deregulation of HIF1-alpha and hypoxia-regulated pathways in hepatocellular carcinoma and corresponding non-malignant liver tissue-influence of a modulated host stroma on the prognosis of HCC,” Langenbeck's Archives of Surgery, vol. 395, no. 4, pp. 395–405, 2010.

Z.-G. Chang, L.-Y. Yang, W. Wang et al., “Determination of high mobility group A1 (HMGA1) expression in hepatocellular carcinoma: a potential prognostic marker,” Digestive Diseases and Sciences, vol. 50, no. 10, pp. 1764–1770, 2005.

B. Shen, E. S. H. Chu, G. Zhao et al., “PPARgamma inhibits hepatocellular carcinoma metastases in vitro and in mice,” British Journal of Cancer, vol. 106, no. 9, pp. 1486–1494, 2012.

A. Csepregi, M. P. A. Ebert, C. Röcken et al., “Promoter methylation of CDKN2A and lack of p16 expression characterize patients with hepatocellular carcinoma,” BMC Cancer, vol. 10, Article ID 3171471, 2010.

Y. Matsuda, T. Ichida, J. Matsuzawa, K. Sugimura, and H. Asakura, “p16(INK4) is inactivated by extensive CpG methylation in human hepatocellular carcinoma,” Gastroenterology, vol. 116, no. 2, pp. 394–400, 1999.

Y. Liu, J. Binz, M. J. Numerick et al., “Hepatoprotection by the farnesoid X receptor agonist GW4064 in rat models of intra- and extrahepatic cholestasis,” Journal of Clinical Investigation, vol. 112, no. 11, pp. 1678–1687, 2003.

V. R. Mas, D. G. Maluf, K. J. Archer et al., “Genes involved in viral carcinogenesis and tumor initiation in hepatitis C virus-induced hepatocellular carcinoma,” Molecular Medicine, vol. 15, no. 3-4, pp. 85–94, 2009.

D. L. Diamond, A. L. Krasnoselsky, K. E. Burnum, et al., “Proteome and computational analyses reveal new insights into the mechanisms of hepatitis C virus-mediated liver disease posttransplantation,” Hepatology, vol. 56, no. 1, pp. 28–38, 2012.

K. Furuta, S. Sato, T. Yamauchi, and S. Kakumu, “Changes in intrahepatic gene expression profiles from chronic hepatitis to hepatocellular carcinoma in patients with hepatitis C virus infection,” Hepatology Research, vol. 38, no. 7, pp. 673–682, 2008.

J. Jiang, L. Yu, X. Huang et al., “Identification of two novel human dynein light chain genes, DNLC2A and DNLC2B, and their expression changes in hepatocellular carcinoma tissues from 68 Chinese patients,” Gene, vol. 281, no. 1-2, pp. 103–113, 2001.

S. Boyault, D. S. Rickman, A. de Reyniès et al., “Transcriptome classification of HCC is related to gene alterations and to new therapeutic targets,” Hepatology, vol. 45, no. 1, pp. 42–52, 2007.

C. Schlaeger, T. Longerich, C. Schiller et al., “Etiology-dependent molecular mechanisms in human hepatocarcinogenesis,” Hepatology, vol. 47, no. 2, pp. 511–520, 2008.

M. W. Welker, W. O. Bechstein, S. Zeuzem, and J. Trojan, “Recurrent hepatocellular carcinoma after liver transplantation—an emerging clinical challenge,” Transplant International, vol. 26, no. 2, pp. 109–118, 2013.

E. R. McDonald III and W. S. El-Deiry, “Suppression of caspase-8- and -10-associated RING proteins results in sensitization to death ligands and inhibition of tumor cell growth,” Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 16, pp. 6170–6175, 2004.

C. H. Chen, S. M. Chuang, M. F. Yang, J. W. Liao, S. L. Yu, and J. J. Chen, “A novel function of YWHAZ/beta-catenin axis in promoting epithelial-mesenchymal transition and lung cancer metastasis,” Molecular Cancer Research, vol. 10, no. 10, pp. 1319–1331, 2012.

A. Budhu, M. Forgues, Q.-H. Ye et al., “Prediction of venous metastases, recurrence, and prognosis in hepatocellular carcinoma based on a unique immune response signature of the liver microenvironment,” Cancer Cell, vol. 10, no. 2, pp. 99–111, 2006.

T. Seki, N. Hayashi, and T. Nishimoto, “RCC1 in the Ran pathway,” Journal of Biochemistry, vol. 120, no. 2, pp. 207–214, 1996.

T. Tachibana, N. Imamoto, H. Seino, T. Nishimoto, and Y. Yoneda, “Loss of RCC1 leads to suppression of nuclear protein import in living cells,” Journal of Biological Chemistry, vol. 269, no. 40, pp. 24542–24545, 1994.

H.-F. Yuen, K.-K. Chan, C. Grills et al., “Ran is a potential therapeutic target for cancer cells with molecular changes associated with activation of the PI3K/Akt/mTORC1 and Ras/MEK/ERK pathways,” Clinical Cancer Research, vol. 18, no. 2, pp. 380–391, 2012.

U. Giri, C. L. Ashorn, L. Ramdas et al., “Molecular signatures associated with clinical outcome in patients with high-risk head-and-neck squamous cell carcinoma treated by surgery and radiation,” International Journal of Radiation Oncology Biology Physics, vol. 64, no. 3, pp. 670–677, 2006.

A. Mortazavi, B. A. Williams, K. McCue, L. Schaeffer, and B. Wold, “Mapping and quantifying mammalian transcriptomes by RNA-Seq,” Nature Methods, vol. 5, no. 7, pp. 621–628, 2008.

J. D. Morlan, K. Qu, and D. V. Sinicropi, “Selective depletion of rRNA enables whole transcriptome profiling of archival fixed tissue,” PLoS One, vol. 7, no. 8, Article ID e42882, 2012.