Table of Contents Author Guidelines Submit a Manuscript
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

1Division of Transplantation, Department of Surgery, University of Washington School of Medicine, Seattle, WA 98195-6410, USA
2Department of Pharmaceutics, University of Washington, Seattle, WA 98121-1023, USA
3Department of Laboratory Medicine, University of Washington Medical Center, Seattle, WA 98104-2499, USA
4Department 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 [1113]. However, few common genes were identified in these studies [1113]. 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).

tab1
Table 1: Patient characteristics: clinical characteristics of the 16 patients from tumor recurrence and tumor nonrecurrence groups (mean ± SD or proportion).
878297.fig.001
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.

878297.fig.002
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).

878297.fig.003
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.

tab2
Table 2: Top canonical pathways: significant molecular pathways regulated in HCC-R.
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)).

fig4
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.

tab3
Table 3: Highly significant upregulated genes in HCC-R tumors with known and unknown functions (references are based on PUBMED search).
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 [5557].

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).

tab4
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).
fig5
Figure 5: Validation of microarray data by qPCR. The values were calculated with Student’s -test (* ).
fig6
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

  1. 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. View at Google Scholar · View at Scopus
  2. 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. View at Publisher · View at Google Scholar · View at Scopus
  3. J. M. Llovet, A. Burroughs, and J. Bruix, “Hepatocellular carcinoma,” The Lancet, vol. 362, no. 9399, pp. 1907–1917, 2003. View at Publisher · View at Google Scholar · View at Scopus
  4. 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. View at Publisher · View at Google Scholar · View at Scopus
  5. 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. View at Publisher · View at Google Scholar · View at Scopus
  6. 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. View at Publisher · View at Google Scholar · View at Scopus
  7. 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. View at Publisher · View at Google Scholar · View at Scopus
  8. 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. View at Publisher · View at Google Scholar · View at Scopus
  9. 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. View at Publisher · View at Google Scholar · View at Scopus
  10. 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. View at Publisher · View at Google Scholar · View at Scopus
  11. 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. View at Publisher · View at Google Scholar · View at Scopus
  12. 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. View at Publisher · View at Google Scholar
  13. 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. View at Publisher · View at Google Scholar · View at Scopus
  14. C. L. Lai, V. Ratziu, M.-F. Yuen, and T. Poynard, “Viral hepatitis B,” The Lancet, vol. 362, no. 9401, pp. 2089–2094, 2003. View at Publisher · View at Google Scholar · View at Scopus
  15. T. Poynard, M.-F. Yuen, V. Ratziu, and C. Lung Lai, “Viral hepatitis C,” The Lancet, vol. 362, no. 9401, pp. 2095–2100, 2003. View at Publisher · View at Google Scholar · View at Scopus
  16. 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. View at Publisher · View at Google Scholar · View at Scopus
  17. 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. View at Google Scholar
  18. 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. View at Publisher · View at Google Scholar · View at Scopus
  19. 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. View at Publisher · View at Google Scholar · View at Scopus
  20. 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. View at Publisher · View at Google Scholar · View at Scopus
  21. 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. View at Publisher · View at Google Scholar · View at Scopus
  22. 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. View at Publisher · View at Google Scholar · View at Scopus
  23. 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. View at Google Scholar · View at Scopus
  24. 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. View at Publisher · View at Google Scholar · View at Scopus
  25. 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. View at Publisher · View at Google Scholar · View at Scopus
  26. 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. View at Publisher · View at Google Scholar · View at Scopus
  27. 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. View at Publisher · View at Google Scholar · View at Scopus
  28. 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. View at Publisher · View at Google Scholar · View at Scopus
  29. 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. View at Publisher · View at Google Scholar · View at Scopus
  30. 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. View at Publisher · View at Google Scholar
  31. 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. View at Google Scholar · View at Scopus
  32. 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. View at Publisher · View at Google Scholar · View at Scopus
  33. 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. View at Publisher · View at Google Scholar · View at Scopus
  34. 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. View at Google Scholar · View at Scopus
  35. 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. View at Publisher · View at Google Scholar · View at Scopus
  36. 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. View at Google Scholar · View at Scopus
  37. 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. View at Google Scholar · View at Scopus
  38. 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. View at Publisher · View at Google Scholar · View at Scopus
  39. 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. View at Publisher · View at Google Scholar · View at Scopus
  40. 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. View at Publisher · View at Google Scholar · View at Scopus
  41. 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. View at Publisher · View at Google Scholar · View at Scopus
  42. 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. View at Publisher · View at Google Scholar · View at Scopus
  43. 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. View at Publisher · View at Google Scholar · View at Scopus
  44. 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. View at Publisher · View at Google Scholar · View at Scopus
  45. 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. View at Publisher · View at Google Scholar · View at Scopus
  46. 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. View at Google Scholar · View at Scopus
  47. 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. View at Publisher · View at Google Scholar · View at Scopus
  48. 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. View at Publisher · View at Google Scholar · View at Scopus
  49. 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. View at Publisher · View at Google Scholar · View at Scopus
  50. 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. View at Publisher · View at Google Scholar · View at Scopus
  51. 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. View at Publisher · View at Google Scholar · View at Scopus
  52. 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. View at Publisher · View at Google Scholar · View at Scopus
  53. 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. View at Publisher · View at Google Scholar · View at Scopus
  54. 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. View at Publisher · View at Google Scholar · View at Scopus
  55. 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. View at Publisher · View at Google Scholar
  56. 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. View at Publisher · View at Google Scholar · View at Scopus
  57. 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. View at Publisher · View at Google Scholar · View at Scopus
  58. 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. View at Publisher · View at Google Scholar · View at Scopus
  59. C. Schlaeger, T. Longerich, C. Schiller et al., “Etiology-dependent molecular mechanisms in human hepatocarcinogenesis,” Hepatology, vol. 47, no. 2, pp. 511–520, 2008. View at Publisher · View at Google Scholar · View at Scopus
  60. 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. View at Publisher · View at Google Scholar
  61. 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. View at Publisher · View at Google Scholar · View at Scopus
  62. 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. View at Publisher · View at Google Scholar
  63. 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. View at Publisher · View at Google Scholar · View at Scopus
  64. T. Seki, N. Hayashi, and T. Nishimoto, “RCC1 in the Ran pathway,” Journal of Biochemistry, vol. 120, no. 2, pp. 207–214, 1996. View at Google Scholar · View at Scopus
  65. 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. View at Google Scholar · View at Scopus
  66. 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. View at Publisher · View at Google Scholar · View at Scopus
  67. 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. View at Publisher · View at Google Scholar · View at Scopus
  68. 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. View at Publisher · View at Google Scholar · View at Scopus
  69. 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. View at Google Scholar