Research Article | Open Access
Molecular Signatures of Recurrent Hepatocellular Carcinoma Secondary to Hepatitis C Virus following Liver Transplantation
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.
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 . 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) . Unfortunately, the recurrence of HCC is a major cause of mortality in surgically treated patients . 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 . RNA derived from FFPE specimens is partially degraded, resulting from the variability in tissue handling/processing, tissue sources, and RNA extraction methods . However, advanced RNA isolation techniques have been developed . 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.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) . 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 . Similar studies by Abramovitiz  and Ton  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 . 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 . 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 . 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 . 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.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).
|Body mass index; **alpha-fetoprotein level.|
***Time specimens removed in surgical procedure until being fixed in formalin.
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.
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).
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 , was downregulated in recurrent HCC tumors.
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)).
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.
| = 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).
FC = fold change (based on our expression data).|
Recurrent HCC following surgical treatment or liver transplantation continues to be a serious health problem . 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 . 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 ; 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 . 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 . The increased expression of RIOK3 has been observed in metastatic head and neck cancers compared with nonrecurrent tumors . 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 , 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.
|HCV:||Hepatitis C virus|
|BMI:||Body mass index|
|IRB:||Institutional Review Board|
|PCRS:||Predicting cancer recurrence score|
|DASL:||cDNA-mediated annealing, selection, extension, and ligation|
|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|
|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|
|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|
|CRIP2:||Cysteine-rich protein 2|
|DBH:||Dopamine beta-hydroxylase (dopamine beta-monooxygenase)|
|DLG4:||Discs, large homolog 4 (Drosophila)|
|DYNLRB2:||Dynein, light chain, roadblock-type 2|
|ESR1:||Estrogen receptor 1|
|FMO3:||Flavin containing monooxygenase 3|
|FMO4:||Flavin containing monooxygenase 4|
|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|
|KLF11:||Kruppel-like factor 11|
|MASP1:||Mannan-binding lectin serine peptidase 1 (C4/C2 activating component of Ra-reactive factor)|
|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|
Conflict of Interests
No conflict of interests was disclosed.
- 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.
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.