PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of molcancBioMed CentralBiomed Central Web Sitesearchsubmit a manuscriptregisterthis articleMolecular CancerJournal Front Page
 
Mol Cancer. 2010; 9: 146.
Published online 2010 June 12. doi:  10.1186/1476-4598-9-146
PMCID: PMC2898705

Integrative and comparative genomics analysis of early hepatocellular carcinoma differentiated from liver regeneration in young and old

Abstract

Background

Hepatocellular carcinoma (HCC) is the third-leading cause of cancer-related deaths worldwide. It is often diagnosed at an advanced stage, and hence typically has a poor prognosis. To identify distinct molecular mechanisms for early HCC we developed a rat model of liver regeneration post-hepatectomy, as well as liver cells undergoing malignant transformation and compared them to normal liver using a microarray approach. Subsequently, we performed cross-species comparative analysis coupled with copy number alterations (CNA) of independent early human HCC microarray studies to facilitate the identification of critical regulatory modules conserved across species.

Results

We identified 35 signature genes conserved across species, and shared among different types of early human HCCs. Over 70% of signature genes were cancer-related, and more than 50% of the conserved genes were mapped to human genomic CNA regions. Functional annotation revealed genes already implicated in HCC, as well as novel genes which were not previously reported in liver tumors. A subset of differentially expressed genes was validated using quantitative RT-PCR. Concordance was also confirmed for a significant number of genes and pathways in five independent validation microarray datasets. Our results indicated alterations in a number of cancer related pathways, including p53, p38 MAPK, ERK/MAPK, PI3K/AKT, and TGF-β signaling pathways, and potential critical regulatory role of MYC, ERBB2, HNF4A, and SMAD3 for early HCC transformation.

Conclusions

The integrative analysis of transcriptional deregulation, genomic CNA and comparative cross species analysis brings new insights into the molecular profile of early hepatoma formation. This approach may lead to robust biomarkers for the detection of early human HCC.

Background

Hepatocellular carcinoma (HCC) is the fifth most common cancer type, and is the third leading cause of cancer mortality worldwide [1,2]. Recent reports show that HCC is becoming more wide-spread and has dramatically increased in North America Western Europe and Japan [2-4]. Additionally, there is an increasing incidence of the disease among younger age groups that warrants further investigation [5,6].

Recently considerable attention has been placed on global gene expression studies as well as genomic aberrations in order to understand the pathogenesis of HCC, and to look for possible early markers of detection [7-13]. Although notable successes have been achieved, there still exist significant challenges due to the heterogeneous nature of HCC (and other cancers) as well as the complexity of the molecular pathogenesis of this disease. Depending on etiological and accompanying pathological conditions, such as viral infection, cirrhosis, inflammation, fibrosis and others, the HCC signature genes identified thus far vary considerably. Additionally, the study of tumor formation in the liver is difficult due to the continuous transcriptome changes that occur during regeneration after hepatectomy [14], as well as age related gene expression changes [15-17]. Similarly, cancer progresses through a series of histopathological stages during which genetic alterations accumulate, and a natural consequence of this are the dynamic changes in gene expression patterns that occur during hepatocellular carcinogenesis. Developing animal models of HCC provide an experimental ground for dissecting the genetic and biological complexities of human cancer and contribute to our ability to identify and characterize pathogenic modifications relevant to early stages of cancer development and progression [18,19]. The previous studies have used cross-species comparative genomics approach successfully to understand the molecular pathogenesis of various cancers [20-23]. Hence, combining cross-species comparative and/or functional genomics approaches with independent datasets from human and animal models of HCC along with genomic DNA copy number alterations enhances the ability to identify robust predictive markers for HCC [23-26].

Here we present a comparative and integrative functional genomics approach to find an early marker for HCC. We developed a rat model and analyzed the transcriptomes of early HCC versus regenerated liver and normal liver in both young and old age animals using a microarray of more than 27,000 annotated genes from Celera and public repositories. We then performed cross-species comparative genomics analysis to identify genes that are conserved in rat and human early HCCs by re-analyzing independent datasets for human early HCC microarray expression profiling data [8,27] and also comparing with the Stanford HCC microarray data [7]. Finally, we performed an integrative analysis of DNA genomic copy number alterations (CNAs) and gene expression profiles (schematically outlined in Figure Figure1).1). Our findings include some genes already reported to be associated with human HCC, thus validating our approach. We also report many other novel genes which were not reported previously in liver cancers. Furthermore, we validated the high expression of eight potential biomarker genes from the blood of patients with early HCC using realtime RT-PCR. Our comparative and integrative genomics approach involving the integration of multiple high dimensional independent datasets may lead to robust biomarkers for the detection of early HCC.

Figure 1
Integrative and cross-species comparative genomics approach to identify evolutionary conserved inter-species biomarkers for early HCC differentiated from liver regeneration. Gene expression signature for early rat HCC is differentiated from liver regeneration ...

Results

Gene Expression Profiling Confirms Pathological Classification

We performed genome-wide gene expression profiling of 24 samples for early HCC, regenerated liver and normal liver of both young and old rats using Applied Biosystems Rat Genome Survey microarray which includes more than 27,000 annotated genes from Celera and public repositories. To find genes that were differentially expressed across three different "treatment" types (i.e. early HCC, regenerated, and normal), and two age groups (old and young), we performed two-factor ANOVA to look for variations due to treatment, age and their interactions. The ANOVA identified 432 and 4063 genes that were significantly modulated by treatment type and age with p < 0.01, respectively. In addition, we found 322 genes that showed a significant interaction of age and treatment effect (data not shown). The unsupervised two-dimensional hierarchical clustering as well as principal component analysis (PCA) using genes which varied significantly with the treatment effect clustered samples according to their treatment type for both old and young (Figure (Figure2A2A and and2B),2B), hence supporting the conclusion that gene expression profiles robustly reflected the histological classification.

Figure 2
Early HCC signature genes conserved across old and young rats. (A) The unsupervised two-dimensional hierarchical clustering using genes that were significantly modulated due to treatment type across all samples (p < 0.01) clustered samples based ...

Identifying HCC Specific Genes Conserved Across Old and Young

The ANOVA identified 432 genes that showed significant expression differences due to treatment in both age groups (p-value < 0.01) were subjected to a template matching algorithm (TMA) [28] to identify HCC specific genes conserved across both age groups. We identified 96 up-regulated and 38 down-regulated genes specific to HCC (template-match p-values < 0.01) (Figure (Figure2C2C and Additional file 1).

The gene ontology and functional network analyses of HCC specific genes were performed using the Ingenuity knowledge base. The biological functions assigned to the dataset are ranked according to the significance of that biological function to the dataset. The enriched functional categories and diseases include carcinogenesis, cell cycle, immune response, cell morphology, cellular development, and growth and proliferation (Figure (Figure2D).2D). The PANTHER also revealed that signal transduction (p-value = 7.17 × 10-4), proteolysis (p-value = 1.59 × 10-3), cell motility (p-value = 3.91 × 10-3), immunity and defense (p-value = 1.07 × 10-2), and cell proliferation and differentiation (p-value = 4.9 × 10-2) were among the most enriched biological processes in the HCC specific genes. The most significantly altered pathways include p53, p38 MAPK, insulin/IGF pathway-protein kinase B signaling cascade, apoptosis, interleukin and integrin signaling pathways. The gene interaction network also corroborated with the altered pathways (Figure (Figure2E2E).

Age-Dependent Differences in Early HCC Differentiated from Regeneration

We found a significant interaction of age and treatment effect in our two-factor ANOVA analysis. Indeed, the hepatic transcriptome changes with ageing as in other cancer types, and age is a potential confounding factor embedded in gene expression profiles [15,16]. Therefore, we stratified samples as young and old cohorts, and identified HCC and regeneration specific genes using one-way ANOVA in each age group separately. The ANOVA identified 925 and 408 significantly dysregulated genes (up or down) due to the different treatment types in old and young animal groups (p-value < 0.02), respectively. The hierarchical clustering in both dimensions (samples and genes), as well as the PCA clearly separated samples based on the treatment type (Additional file 2), The gene expression clustering distance between the HCC group and the other two groups (regenerated and normal) was the greatest in both age groups (Additional file 2).

Early HCC signature genes in each age group were obtained by overlapping gene lists. Each circle in the Venn diagram represents the differential expression between two treatment types. We identified 80 genes and 100 genes specific to hepatoma in young and old, respectively (Figure (Figure3A3A and Additional file 3). As seen from the heatmap of HCC specific genes, these sets of genes were exclusively up/down regulated in the HCC group only (Figure (Figure3B3B and Additional file 3). The expression of Pbsn, Lum, Adam8, Ctse, Calb3, Fbn1, Agtpbp1, Prom1, Ela1, Tnfsf13, and Ap2b1 were significantly altered in both young and old rats with early HCC. The genes A2m, Cdh13, Mas1, Slack, Cidea, and Dcn were significantly dysregulated exclusively in HCC in the young; whereas Cxcl5, Lox, Slc25a2, Rmt1, and Nid2, were specific to HCC in old rats. We also identified regeneration specific genes in old and young rats in a similar approach (data not shown).

Figure 3
Heatmap and gene interaction networks of early HCC specific genes in young. (A) Venn diagram characterizing differential gene expression between and specific to different treatment types: early rat HCC (DY), regenerated (RY), and normal (NY). The number ...

Functional Comparison of Hepatoma and Regeneration in Young and Old

Early HCC signature genes in young animals were mainly associated with cancer, cell cycle, immune response, cellular function and maintenance, development, cell adhesion-mediated signaling, proteolysis, and signal transduction, whereas genes in old animals were highly associated with cancer, cellular movement, extracellular transport and import, cell adhesion, tissue development, cell morphology, cell-to-cell signaling and interaction. On the other hand, the regeneration specific genes were mainly associated with mRNA transcription and regulation, lipid metabolism, protein modification, protein phosphorylation, cell morphology, cellular development, small molecule biochemistry, and cellular growth and proliferation (Table (Table1).1). To gain more insights into HCC pathogenesis for young and old age early HCCs, we carried out gene interaction networks of HCC specific genes in young and old (Figure 3C, D and and3E3E and Additional file 3, respectively). The interaction networks highlight the important role of p53, p38 MAPK, ERK/MAPK, PI3K/AKT signaling, NF-κB and TGF-β pathways in early rat HCC.

Table 1
Functional comparison of hepatoma and regeneration in young and old

Cross-Species Comparative Genomic Analysis

To identify how many of our early rat HCC signature genes were conserved in early human HCCs, we re-analyzed two independently performed microarray datasets for early human HCC from Wurmbach et. al. [8] and Mas et. al. [27]. The Wurmbach et. al.'s dataset composed of 19 early HCC patients and 10 normal controls, and the Mas et. al. dataset was composed of 16 cirrhotic livers with early HCC, 38 HCV associated HCC, and 19 normal livers. In addition, we also compared our signature genes with the OncodB.HCC database for 57 HCC patients from the Stanford HCC microarray data. The comparison of our rat early HCC signature genes (human orthologous) with the re-analyzed early human HCCs datasets revealed that 154 unique genes were conserved with early human HCCs (p < 0.001), and 35 of those were shared by all datasets analyzed (Table (Table2,2, Figure Figure1).1). We found that unsupervised clustering using the conserved signature gene list across species was sufficient to separate individuals in both Wurmbach's and Mas' human HCC samples as either early HCC patients or normal controls (data not shown). The gene interaction network analysis of the 154 signature genes indicated the importance of NF-κB, RAS and JNK activation in early hepatoma formation (Additional file 4). The network analysis of 35 cross-species conserved early HCC signature genes reveals the important roles of ERK/MAPK, PI3K/AKT, and TGF-β pathways (Figure (Figure4).4). It also indicated a potential critical regulatory role of MYC, ERBB2, HNF4A, and SMAD3 for malignant transformation to early HCC (Figure (Figure4B4B).

Figure 4
The gene interaction networks of early HCC potential biomarker genes that are conserved in rat early HCCs and in multiple independent human early HCCs. The network analysis of 35 early HCC signature genes indicated the activation of ERK/MAPK, PI3K/AKT ...
Table 2
List of 35 cross-species conserved early HCC signature genes with qRT-PCR and independent human/rat early HCC validations overlaid

Integrative Analysis of Transcriptional Deregulation with Genomic Copy Number Aberrations

Various studies have reported chromosomal instability at chromosomal regions associated with many cancers, including human HCC copy number (CN) status [7,9,26,29,30]. These genomic modifications, which in part are reflected in changes in DNA copy number (CN), may alter the transcriptional control mechanism, and hence impact gene expression levels [31,32]. Hence, we compiled genes located in CNA regions reported in three independent genome copy number studies of human HCC [7,9,10]. The integration of 154 early HCC signature genes with the copy number data resulted in 75 genes that mapped to human CNA chromosomal loci genes including COL1A1, CCNA2, NFATC2, F2, DCK, MMP2, GJA1, VIM, LGALS3BP, and SP100. The interaction network of those genes further corroborated with the activation of the NF-κB, p38 MAPK, AP1, and JNK pathways (Figure (Figure5).5). We found that almost 50% of the 35 cross-species conserved signature genes that are common to different types of early HCCs analyzed were mapped to genomic locations within the CNA regions (Table (Table22).

Figure 5
The interaction network analysis of 75 early HCC signature genes conserved across species and having genomic alterations. The network analysis of 75 cross-species conserved signature genes with CN alterations indicated the importance of NF-κB, ...

Independent Validation Set Analysis

As a validation of our results, we analyzed four independently performed microarray datasets for early human and rat HCCs [19,33-35] using the analysis procedure defined in the "Methods" section on the new datasets. The first validation dataset was from Chiang et.al. [33] using Affymetrix short oligo arrays. The dataset was composed of 91 HCV-related HCC tumor samples, of which 65 were very early or early stage disease, which we used in our re-analysis and comparison. The re-analyzed validation dataset showed a significant number of genes (p < 10-5) in common with our analysis results. In fact, more than 50% of our cross-species conserved genes were also differentially expressed in early HCC compared to normal controls in the validation dataset (Table (Table2).2). The significance of overlaps was calculated using hypergeometric distributional assumption [36] and p-values were adjusted using Bonferroni correction for multiple comparisons [37]. In addition, unsupervised clustering was performed using our 35- gene signature to cluster the samples from Chiang et. al. We found that using our signature gene list was sufficient to separate individuals in Chiang et. al.'s study as either early HCC patients or normal controls (Additional file 5).

Moreover, we found a significant number of genes in common with the second dataset from Boyault et. al. [34] which consisted of 57 human HCCs and five samples of pooled non-tumorus tissues, and a third gene expression dataset from Liao et. al. [35] consisting of human HCC from various stages (we used expression data for only the early stage of the disease) (Table (Table2).2). Furthermore, we obtained consistent results with the DEN-induced HCC in rats from [19,38]. The IPA functional and network analysis of all validation datasets revealed a significant number of overrepresented functional categories and pathways in common with our results. Of note, cell death, cancer, cellular development, cellular growth and proliferation, organismal development, transport, and cell cycle came up as significantly enriched categories in both the validation datasets and our analyses. The interaction networks analyses of significantly dysregulated genes in validation datasets highlight the important roles of MYC, ERK/MAPK, AKT, NF-κB and TGF-β signaling pathways. The similar findings between our results and the independent validation sets argue against random chance accounting for the observed enrichment of these functional categories and pathways.

Validation of Microarray Data for Early Rat HCC by Realtime RT-PCR

To confirm the microarray results by an independent method, we validated expression levels of six randomly selected differentially regulated genes (Pbsn, Cdh13, Lum, Nid2, Dcn, Slc22a5) in early rat HCC by realtime quantitative RT-PCR. A highly significant correlation existed between the microarray and realtime RT-PCR results (r = 0.97, p value < 0.001) (Figure (Figure6),6), thus demonstrating the reliability of our gene expression measurements. The selected genes and their interaction networks with other genes are shown (Figures 3C, D and and3E,3E, ,4A,4A, and Additional file 3).

Figure 6
Confirmation of the microarray gene expression for six randomly selected significantly regulated genes in rat early HCC by realtime qRT-PCR. Ratio of expression (fold change) for each gene in (A) early HCC in young (DY) compared to normal (NY); (B) DY ...

Validation of Potential Biomarker Genes from Whole Blood of Patients with Early HCC Using qRT-PCR

To further validate the differential expression of potential biomarker genes using realtime RT-PCR from the blood of early HCC patients and healthy control subjects (ten subjects in each group), we selected eight genes (GJA1, VIM, IGFBP3, COL1A1, SP100, MMP2, LGALS3BP, and DPP4) among early HCC gene signature with CNA (denoted with asterisk in Table Table2)2) and were differentially regulated in at least one of the independent datasets. We confirmed a statistically significant increase in the expression of these biomarker genes in early HCC patients relative to healthy control subjects (p-value < 0.05) (Table (Table2,2, Figure Figure7);7); hence demonstrating the robustness of the cross-species integrated genomics procedure.

Figure 7
Differential expression of a subset of genes was confirmed in whole blood of human early HCC subjects with qRT-PCR. The up-regulation of expression of eight genes from Table 2 was confirmed in blood of early HCC patients compared to normal controls by ...

Discussion

The present study sought to identify evolutionarily conserved inter-species biomarkers for early HCC differentiated from liver regeneration using integrative and cross-species comparative genomic approaches. The main contributions of this study are as follows: First, we developed a rat model of liver regeneration post-hepatectomy (return to quiescence), as well as liver cells undergoing malignant transformation and compared them to normal liver using a comprehensive microarray of 27,000 publicly available and Celera annotated rat genes. We included the liver regeneration in our model, as regeneration is a critical component in the surgical treatment of HCC, and frequently associated with HCC occurrence [39]. Though liver cells can regenerate, they do not typically transform and lead to HCC [40]. Therefore, an early HCC marker needs to be unique for the tumorigenic process and not overlap with the transcriptome changes that occur in regenerating or normal liver tissue. As ageing is also known to be a confounding factor embedded in gene expression profile data [15-17], we included age as a factor in our multi-factor statistical analysis, and identified age-specific differences in early HCC. Secondly, we performed cross-species comparative analysis to identify genes that are conserved in early rat HCCs and in multiple independently performed early human HCCs which would facilitate the identification of critical regulatory modules conserved across species in the expression profiles. Finally, we integrated genomic CNA data associated with the human HCCs with the transcriptomic profile, and performed validation analyses both in silico and with quantitative realtime RT-PCR (as schematically outlined in Figure Figure1).1). As CNAs have clear impact on expression levels in a variety of tumors [26,29,30], this dual strategy is very effective for interpreting the DNA and RNA level anomalies in cancer, in order to identify genes involved with tumor initiation and progression [24,26].

The validation analyses demonstrated great concordance of our results with other data sets using various microarray platforms. The ABI 1700 system has a unique approach in identifying dysregulated genes since it targets genes from both Celera and Public databases and utilizes chemiluminescently enhanced detection that is likely to determine relatively rare mRNAs. Also, our confirmatory quantitative realtime RT-PCR experiments displayed a strong correlation with the microarray results, adding to the validity of the present observations. This is in agreement with some recent studies showing a linear relationship for real-time and conventional reverse transcription and therefore validates the robustness of mRNA quantification using either microarrays or quantitative RT-PCR[41]. Hence this allowed us to identify potential biomarkers for human early HCC and to gain further insight into the mechanism of early hepatoma formation.

We performed a two-step algorithm to identify early rat HCC signature genes: In the first step, a two-way ANOVA was performed including treatment and age as well as their interactions into our statistical model and we identified genes uniquely expressed in early HCC in both young and old rats (Figure (Figure2).2). In the second step, because the interaction between age and treatment was significant, we stratified our samples as young and old cohorts, and HCC specific genes were identified using two one-way ANOVA in each age group separately [42]. Finally, the HCC specific genes in both young and old identified in the two steps were combined before performing the cross-species comparison (as detailed in "Material and Methods" section, and schematically shown in Figure Figure1).1). When comparing the early HCC group with the regeneration or normal cohorts to identify the differentially expressed genes we used a set of criteria: S/N ratio > 3 in > 50% of the samples, a p-value < 0.02 and absolute fold change > 1.8. These observations are consistent with the study of Guo et. al., in which gene lists ranked by fold change and filtered with non-stringent statistically significant tests were more reproducible across platforms than those generated through other analytical procedures [43,44]. In addition, as Ghosh et. al. discuss on combining data from multiple gene expression studies, if two studies independently discover that the same gene/protein to be differentially expressed, then the chance of error is significantly reduced [45].

The comparison of our signature genes with three different independently performed early human HCC microarray data sets revealed a significant number of early rat HCC genes (human orthologous) conserved across early human HCCs (p < 0.001). Indeed, many of those genes were related to cancer. For example, LUM, CCNA2, IGFBP3, HPX, COL1A1, SRPX, VIM, TGFBR1, DCN, MMP2, CD14, DCK, BIRC3, GJA1, LOX, SP100, PROM1 and CREB1 were known to regulate tumorigenesis, neoplasia, apoptosis, growth, differentiation and proliferation. Some of the most significantly activated canonical pathways included hepatic fibrosis/hepatic stellate cell activation (CD14, COL1A1, FGFR2, IGFBP3, MMP2, and TGFBR1), and p38 MAPK signaling pathways (CREB1, PLA2G2A, TGFBR1). The network analysis of early HCC signature genes indicated the activation of ERK/MAPK, PI3K/AKT, and TGF-β signaling pathway, and a potential critical regulatory role of MYC, ERbB-2, HNF4A, and SMAD3 for early HCC (Figures (Figures2,2, ,33 and and4).4). MAPKs are implicated in diverse cellular processes such as cell survival, differentiation, adhesion, and proliferation [46]. The gene network analysis of differentially expressed genes further confirmed the altered pathways. Moreover, it also indicated the importance of NF-κB, RAS and JNK activation in early hepatoma formation (Figure (Figure55 and Additional file 4).The role of MYC in various types of carcinogenesis has been extensively investigated [47]. Most recently, JNK1 activation[48], and increased expression of ErbB-2 were found to be associated with HCC [49]. Thus, our current findings are consistent with previously performed independent cancer studies, including those for HCC. However, the novelty of our approach is that using comparative and integrative genomics, we provide evidence for the potential central role of these genes in the earliest phase of liver malignant transformation.

Our comparative genomics analysis resulted in a 35-gene cross-species conserved signature for all types of early HCCs. Over 70% of the conserved genes were associated with cancer according to the IPA knowledgebase, including LGALS3BP[50,51], VIM[52,53], DCN[54,55], IGFBP3[56], FGFR2[57], GJA1[58], SP100[59], DPP4[60], PROM1[61], BIRC3[62], MMP2[63], and COL1A1[64,65] (Table (Table2).2). Furthermore, using literature mining tools, such as MILANO [66], we found that almost 90% of our signature genes were reported to be cancer related.

There are areas of genomic instability reported in many cancers, including HCC, and some regions commonly exhibit either deletion or increased gene dosage, leading to changes in DNA copy number (CN) [9,26,29,30]. Integrating the gene expression with the CN data reveals the chromosomal regions with concordantly altered genomic and transcriptional status in tumors [24,32,67]. Hence, focusing on differentially expressed genes with concomitant altered DNA copy number may identify novel early HCC markers of malignant transformation and progression. The presence of altered DNA CN and LOH may contribute to cancer formation [30,31,68]. Therefore, the pattern of genomic modifications in a tumor represents a structural fingerprint that may include the transcriptional control mechanisms and locally impact gene expression levels [31,32]. We identified that more than 50% of our cross-species conserved early HCC signature genes were found to be copy number dependent (Table (Table22).

We found significant expression of LGALS3BP (Lectin, galactoside binding soluble 3 binding protein) and COL1A1 located on Chromosome 17q. The LGALS3BP is a 90-kD protein, designated serum protein 90 K that was found at elevated concentrations in the serum of patients with various types of breast, lung, colorectal, ovarian, and endometrial cancer [50,51]. It is a secreted glycoprotein that binds galectins, beta1-integrins, collagens, and fibronectin, and has some relevance in cell-cell and cell-extracellular matrix adhesion [69]. Another gene which could be a potential biomarker for early HCC is dipeptidyl peptidase IV (DPP4). DPP4 is a serine protease, which plays an important role in immune regulation, signal transduction, and apoptosis. It has been shown that DPP4 may have a critical function in tumor progression in several human malignancies [60,70]. Matrix metalloproteinases (MMP) also are involved with early carcinogenic events, tumor growth, tumor invasion and metastasis [63,71,72]. Matrix metalloproteinases (MMPs) are zinc-dependent endopeptidases that cleave and degrade a wide spectrum of extracellular matrix components, and are involved with extracellular matrix remodeling during the process of tumor invasion and metastasis [72]. Alterations in MMP expression and their endogenous inhibitor (TIMP) may contribute to HCC metastasis [71-73].

It is worth mentioning that, the gene "Probasin" (Pbsn), was significantly up-regulated (fold change > 15 in both old and young early rat HCCs). The high expression of Pbsn in our rat model was also confirmed with realtime RT-PCR. Pbsn is a member of the lipocalin family and has not yet been associated with HCC in rats and has no known human ortholog. However, it has been shown that Pbsn is highly expressed in prostate and implicated in both benign prostatic hyperplasia and prostate cancer [74-76] and taste bud tumorigenesis in rats [74]. Also since the promoter of this gene exhibits strong androgen receptor-specific and tissue-specific regulation, Pbsn is proposed to be a potential candidate for targeted therapies for advanced prostate cancer [77].

We have also found significant expression of lumican (LUM) and decorin (DCN) in both early rat and human HCCs. LUM and DCN are members of a small leucine-rich proteoglycan (SLRP) family. Lumican has been shown to participate in the maintenance of tissue homeostasis and modulation of cellular functions including cell proliferation, migration, adhesion, and differentiation [78]. Decorin has been reported to have a number of functions including suppressing cancer cell growth and metastasis andacting with extracellular matrix molecules to influence cell adhesion and fibril stability [55]. DCN acts as a natural inhibitor of TGF and is considered to be a specific antagonist of EGFR [54]. In addition, the altered expression of lumican and decorin has been associated with various human cancers including breast, pancreatic, lung, ovarian, melanoma, colorectal, osteosarcoma and ductal adenocarcinoma [54,78-82].

Genes whose protein products are released into the extracellular space would be ideal tumor markers for clinical applications, as it would be possible to detect these proteins in patients' biological fluids rather than through the use of invasive biopsies. Moreover, previous studies have found that cells derived from peripheral blood could be used to assess disease-associated gene signatures [83-89]. In our study, we confirmed the high expression of eight selected candidate biomarker genes (GJA1, VIM, IGFBP3, COL1A1, SP100, MMP2, LGALS3BP, and DPP4) by using realtime RT-PCR from the blood of early HCC patients. These genes and other potential biomarker genes identified through our integrated-comparative genomics approach (listed in Table Table2)2) and their encoded proteins will be further studied in a large cohort of patients to determine if they have a role in early HCC pathology and if they could be novel early HCC biomarkers detectable in biological fluids.

Conclusions

In summary, to our knowledge, this is the first study to examine HCC differentiated from regeneration in both old and young rats, and coupled with a cross-species comparative and integrative genomics approach to identify genes that could be potential biomarkers for early human HCC. The results of our study include the depiction of refined and delineated biological pathways differentially modulated in HCC that is built around TP53, p38 MAPK, ERK/MAPK, PI3K/Akt, NF-κB, TGF-β, MYC, and ERbB-2, including their target genes that were not previously implicated with early HCC. Our cross-species comparative and integrative genomics approach which involved integration of multiple high dimensional independent datasets has led to potentially robust biomarkers for the detection of early HCC. The signature genes that we identified could be considered as "evolutionarily conserved cross-species biomarkers for early HCC with genomic copy number alterations". Further studies are needed to identify if any of the potential biomarkers identified in this study can be readily and reproducibly detected in blood, urine or other bodily fluids. This could then form the basis of a useful diagnostic test for the detection of early HCC.

Methods

Animals

Male Sprague-Dawley rats were maintained at the King Fahad National Centre for Children's Cancer and Research Animal Facility. This facility is managed in accordance with AALAS regulations. Ten young adult (5 months) and ten old adult (21 months) animals were subjected to partial hepatectomy. Actual survival rates allowed for four animals in each group to be analyzed. The re-growth of one lobe of the liver was completed within one month, by which time the liver cells again became quiescent, which was confirmed by histological analysis. In parallel, separate animals were treated with diethylnitrosoamine (200 mg/kg), which was injected intraperitoneally to induce the formation of HCC. Once the early HCC formation became apparent within 2-4 weeks, the rats were sacrificed. In the partial hepatectomized animals, one unaffected lobe and the regenerated lobe of the liver were removed independently. In the carcinogenic treated animals the tumors were carefully dissected to avoid removing normal tissue. All tissues were snap-frozen and stored at -80°C until required for RNA isolation. Small pieces of tissue were removed for formalin fixation to be used for histological examination.

Human Subjects

Twenty blood samples were collected for this study (10 early HCC and 10 from healthy controls). Histopathological classification of HCC and clinical staging of early HCC were performed according to International Working Party [90] as previously described [8]. Patients diagnosed with the early HCC and healthy controls were recruited under an institutional review board-approved project (RAC# 2060040); all subjects provided written, informed consent before entry in the study. A total of 4 ml (in two separate PaxGene tubes) of whole blood samples were collected for each individual according to manufacturer's guidelines (QIAGEN Inc., Valencia, CA, USA). The total RNA isolation was performed using PreAnalytiX - PAXgene Blood RNA System (QIAGEN Inc.) by strictly following the manual and protocols provided with the kit-system.

Microarray Hybridization

Total RNA was isolated according to standard protocols. Quality Control of RNA was done using Bioanalyzer 2100 RNA 6000 NanoAssay and RNA above RIN = 8 was included to the study (Agilent Technologies, Santa Clara, USA). Rat Genome Survey Microarray (Applied Biosystems, Foster City, CA, USA) was utilized for microarray studies. cDNA synthesis, cRNA and labeling, chemiluminescence detection, image acquisition and analysis were performed following the manufacturer's protocols, guidelines and recommendations.

Microarray Data Analysis

Images were auto-gridded and the chemiluminescent signals were quantified, then background subtracted using the Applied Biosystems 1700 Chemiluminescent Microarray Analyzer software v 1.1. For transcriptome analysis, detection thresholds were used following the manufacturer's recommendations. Detection threshold was set as S/N > 3 and quality flag < 5000. The microarray data were analyzed from 24 samples (2 samples were excluded for quality reasons). The open source R software package http://www.r-project.org and tools from the BioConductor project were used for normalization and determination of differentially expressed genes [91]. Two-factor Analysis of Variance (ANOVA) was performed to include both "treatment" (HCC, regenerated and normal, which will be referred to as treatment in the remainder of the manuscript), as well as age (old and young) factor together with feature selection algorithm (also known as template matching(TMA)) [28] to look for treatment as well as age specific variation. Significantly modulated genes specific to HCC were defined as those with ANOVA (treatment) p- value < 0.01, and TMA p-value < 0.01. Additionally, samples were stratified as young and old cohorts, and HCC specific genes identified using one-way ANOVA in each age group separately. When comparing HCC group with regenerated and normal controls to identify the differentially expressed genes specific to the HCC, we used a combination of three criteria. We considered genes that are "present" in at least half of the samples in either group. HCC specific genes were defined as those with in absolute fold change > 1.8 and p-value < 0.02. These observations are consistent with the study of Guo et. al., in which gene lists ranked by fold change and filtered with non-stringent statistically significant tests were more reproducible across platforms than those generated through other analytical procedures [43]. A validation datasets were generated from three independent human HCC studies by Chiang et. al. [33] (GSE9843), which was composed of 91 HCV-related HCC tumor samples, of which 65 were very early or early stage disease (we used in our re-analysis only very early and early HCC datasets) and Boyault et. al. [34] (E-TABM-36) which consisted of 57 human HCCs and five samples of pooled non-tumorus tissues, and from Liao et. al. [35] (GSE 6222) consisting of various stages of HCC (we used expression data for only the early stage of the disease). Furthermore, we compared our results with the results from two independent studies [19,38] with the DEN-induced HCC in rats. The raw data was analyzed by using dChip[92] and open source R/Bioconductor packages. The dChip outlier detection algorithm was used to identify outlier arrays, and probes "present" in at least 50% of the samples in either group were filtered. The data was normalized by the GC Robust Multi-array Average (GC-RMA) algorithm [93,94]. Unpaired t-tests were performed to determine significant differences in gene expression levels between patients and normal controls. The Hierarchical clustering using Pearson's correlation with average linkage clustering was performed by MeV 4.0 [95].

Information about genes participating in known biological process and pathways were derived by using DAVID Bioinformatics Resources[96], Expression Analysis Systematic Explorer (EASE)[97], and PANTHER (Protein ANalysis THrough Evolutionary Relationships) Classification Systems [98]. For each molecular function, biological process or pathway term, PANTHER calculates the number of genes identified in that category in both a list of differentially regulated genes and a reference list containing all the probe sets present on the chip and compares these results using the binomial test to determine if there are more genes than expected in the differentially regulated list [99]. Over-representation is defined by p < 0.05. Statistical analyses were performed with the MATLAB software packages (Mathworks, Natick, MA, USA), R and Bioconductor and PARTEK Genomics Suite (Partek Inc, St. Lois, MO, USA).

Functional Pathway and Network Analysis

Functional pathway, gene ontology and network analyses were executed using Ingenuity Pathways Analysis (IPA) 6.3 (Ingenuity Systems, Mountain View, CA). The differentially expressed signature gene lists for hepatoma and regeneration in different age groups were mapped to its corresponding gene object in the Ingenuity pathway knowledge base. These so-called focus genes were then used as a starting point for generating biological networks. A score was assigned to each network in the dataset to estimate the relevance of the network to the uploaded gene list. This score reflects the negative logarithm of the P that indicates the likelihood of the focus genes in a network being found together due to random chance. Using a 99% confidence level, scores of 2 were considered significant. Significances for biological functions or pathways in the signature genes for such functions or pathways compared with the ABI Rat Genome Survey Microarray as a reference set. A right-tailed Fisher's exact test was used to calculate a p-value determining the probability that the biological function (or pathway) assigned to that data set is explained by chance alone.

Cross-Species Comparative and Integrative Genomic Analysis

Human early HCC datasets from two independent studies by Mas et. al. [27] using Affymetrix HG-U133A 2.0 array, and Wumbach et. al. [8] using Affymetrix HG-U133 Plus 2.0 were re-analyzed. The raw data were analyzed using R/Bioconductor packages and Partek Genomics Suite (Partek Inc.). The data were normalized by the GC Robust Multi-array Average (GC-RMA) algorithm. Unpaired t-tests were performed to determine significant differences in gene expression levels between patients and normal controls. The cross mapping of Applied Biosystems Rat Genome Survey microarray probes were mapped to human orthologs through "AB1700 rat annotation spreadsheet" designed by Applied Biosystems on the basis of sequence identity. The transcripts present on both platforms (AB1700 and Affymetrix) were identified using Resourcerer [100]. Genes within copy number altered regions based on three independent genome CNA studies of human HCC [7,9,10] were determined using NCBI MapViewer http://www.ncbi.nlm.nih.gov/mapview, and integrated those with the gene expression profiling data (Figure (Figure11).

Realtime RT-PCR Experiments

Confirmatory realtime RT-PCR experiments were performed using the ABI 7500 Sequence Detection System (Applied Biosystems). 50 ng total RNA procured from the same microarray study samples were transcribed into cDNA using a Sensicript Kit (QIAGEN Inc., Valencia, CA, USA) under the following conditions: 25°C for 10 min, 42°C for 2 hrs, and 70°C for 15 min in a total volume of 20 μl. Six differentially expressed rat genes (Pbsn, Cdh13, Lum, Nid2, Dcn, Slc22a5) and eight human genes (GJA1, VIM, IGFBP3, COL1A1, SP100, MMP2, LGALS3BP, and DPP4) were selected and primers designed using Primer3 software. For the human samples, blood total RNA was utilized. After primer optimization, realtime PCR experiments were performed with 6 μl cDNA using Quantitech SyBr Green Kit (QIAGEN), employing GAPDH as the endogenous control gene. All reactions were conducted in triplicates and the data was analyzed using the delta delta CT method [101,102].

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

DC and NK conceived the research problem and designed the methodology. MAC developed the rat model. MG provided suggestions for the rat model. NK and AA carried out microarray experiments. DC collected and processed the data, performed statistical and bioinformatics analyses, and drafted the manuscript together with NK and BHP. JQ and MMS provided suggestions for method design and commented on the manuscript. MAC, PTO, BHP, and AAQ provided general interpretation of results. All authors read and approved the final manuscript.

Supplementary Material

Additional file 1:

Selected HCC specific genes, conserved across both age groups (old and young), and significantly modulated with respect to regenerated and normal liver.

Additional file 2:

Comparison of expression profiles of HCC and regeneration within the same age group. (A, D) Heatmap of significantly dysregulated genes due to different treatment types in young and old, respectively. (B, E) Hierarchical clustering of samples separated based on treatment type in young and old, respectively. The gene expression clustering distance between the HCC group and other two groups (regenerated and normal) was the greatest in both age groups (C, F) Principle component analysis (PCA) which contained almost 76% of the variance in the data matrix clearly separated samples based on the treatment type in young and old, respectively.

Additional file 3:

Heatmap and gene interaction networks of HCC specific genes in the old age group. (A) Venn diagram characterizing differential gene expression between and specific to different treatment types (the HCC, the regenerated and the normal). The number of HCC specific genes, 100, is circled in black. (B) Heatmap of HCC specific genes exclusively dysregulated (up/down regulated) in the HCC group only. (C-E) Functional network analysis of HCC specific genes. Top three scoring gene interaction networks (with highest relevance scores) are shown. Nodes represent genes, with their shape representing the functional class of the gene product, and edges indicate biological relationship between the nodes (see legend in Figure Figure2).2). (F) Top network functions associated with three networks shown. An IPA score of three indicates that there is 1/1000 (score = -log (p-value)) chance that the focus genes are assigned to a network randomly.

Additional file 4:

The gene interaction network analysis early HCC signature genes that are conserved in rat early HCC and in either of multiple human early HCCs (A, B) The top two scoring gene interaction networks of 154 cross-species conserved signature genes indicated the importance of NF-κB, RAS and JNK activation in early hepatoma formation. Nodes represent genes, with their shape representing the functional class of the gene product, and edges indicate biological relationship between the nodes (see legend in Figure Figure55).

Additional file 5:

The unsupervised Principle Components Analysis (PCA) was performed using our 35-gene signature to cluster samples from independent validation dataset of Chiang et. al. Our signature gene list was sufficient to separate individuals in Chiang et al.'s study as either early HCC patients or normal controls.

Acknowledgements

We thank Maqbool Ahmad for providing help in RNA isolation, and Bedri Karakas for critically reading the manuscript. We would like to thank to King Faisal Specialist Hospital and Research Center for the financial support.

References

  • El-Serag HB, Rudolph KL. Hepatocellular carcinoma: epidemiology and molecular carcinogenesis. Gastroenterology. 2007;132:2557–2576. doi: 10.1053/j.gastro.2007.04.061. [PubMed] [Cross Ref]
  • Altekruse SF, McGlynn KA, Reichman ME. Hepatocellular carcinoma incidence, mortality, and survival trends in the United States from 1975 to 2005. J Clin Oncol. 2009;27:1485–1491. doi: 10.1200/JCO.2008.20.7753. [PMC free article] [PubMed] [Cross Ref]
  • Nguyen MH, Whittemore AS, Garcia RT, Tawfeek SA, Ning J, Lam S, Wright TL, Keeffe EB. Role of ethnicity in risk for hepatocellular carcinoma in patients with chronic hepatitis C and cirrhosis. Clin Gastroenterol Hepatol. 2004;2:820–824. doi: 10.1016/S1542-3565(04)00353-2. [PubMed] [Cross Ref]
  • Farazi PA, DePinho RA. Hepatocellular carcinoma pathogenesis: from genes to environment. Nat Rev Cancer. 2006;6:674–687. doi: 10.1038/nrc1934. [PubMed] [Cross Ref]
  • El-Serag HB. Hepatocellular carcinoma: recent trends in the United States. Gastroenterology. 2004;127:S27–34. doi: 10.1053/j.gastro.2004.09.013. [PubMed] [Cross Ref]
  • El-Serag HB. Epidemiology of hepatocellular carcinoma in USA. Hepatol Res. 2007;37(Suppl 2):S88–94. doi: 10.1111/j.1872-034X.2007.00168.x. [PubMed] [Cross Ref]
  • Su WH, Chao CC, Yeh SH, Chen DS, Chen PJ, Jou YS. OncoDB.HCC: an integrated oncogenomic database of hepatocellular carcinoma revealed aberrant cancer target genes and loci. Nucleic Acids Res. 2007;35:D727–731. doi: 10.1093/nar/gkl845. [PubMed] [Cross Ref]
  • Wurmbach E, Chen YB, Khitrov G, Zhang W, Roayaie S, Schwartz M, Fiel I, Thung S, Mazzaferro V, Bruix J. Genome-wide molecular profiles of HCV-induced dysplasia and hepatocellular carcinoma. Hepatology. 2007;45:938–947. doi: 10.1002/hep.21622. [PubMed] [Cross Ref]
  • Luo JH, Ren B, Keryanov S, Tseng GC, Rao UN, Monga SP, Strom S, Demetris AJ, Nalesnik M, Yu YP. Transcriptomic and genomic analysis of human hepatocellular carcinomas and hepatoblastomas. Hepatology. 2006;44:1012–1024. doi: 10.1002/hep.21328. [PMC free article] [PubMed] [Cross Ref]
  • Woo HG, Park ES, Lee JS, Lee YH, Ishikawa T, Kim YJ, Thorgeirsson SS. Identification of potential driver genes in human liver carcinoma by genomewide screening. Cancer Res. 2009;69:4059–4066. doi: 10.1158/0008-5472.CAN-09-0164. [PMC free article] [PubMed] [Cross Ref]
  • Ye QH, Qin LX, Forgues M, He P, Kim JW, Peng AC, Simon R, Li Y, Robles AI, Chen Y. Predicting hepatitis B virus-positive metastatic hepatocellular carcinomas using gene expression profiling and supervised machine learning. Nat Med. 2003;9:416–423. doi: 10.1038/nm843. [PubMed] [Cross Ref]
  • Shackel NA, Seth D, Haber PS, Gorrell MD, McCaughan GW. The hepatic transcriptome in human liver disease. Comp Hepatol. 2006;5:6. doi: 10.1186/1476-5926-5-6. [PMC free article] [PubMed] [Cross Ref]
  • Smith MW, Yue ZN, Geiss GK, Sadovnikova NY, Carter VS, Boix L, Lazaro CA, Rosenberg GB, Bumgarner RE, Fausto N. Identification of novel tumor markers in hepatitis C virus-associated hepatocellular carcinoma. Cancer Res. 2003;63:859–864. [PubMed]
  • Nam SW, Park JY, Ramasamy A, Shevade S, Islam A, Long PM, Park CK, Park SE, Kim SY, Lee SH. Molecular changes from dysplastic nodule to hepatocellular carcinoma through gene expression profiling. Hepatology. 2005;42:809–818. doi: 10.1002/hep.20878. [PubMed] [Cross Ref]
  • Geigl JB, Langer S, Barwisch S, Pfleghaar K, Lederer G, Speicher MR. Analysis of gene expression patterns and chromosomal changes associated with aging. Cancer Res. 2004;64:8550–8557. doi: 10.1158/0008-5472.CAN-04-2151. [PubMed] [Cross Ref]
  • Yau C, Fedele V, Roydasgupta R, Fridlyand J, Hubbard A, Gray JW, Chew K, Dairkee SH, Moore DH, Schittulli F. Aging impacts transcriptomes but not genomes of hormone-dependent breast cancers. Breast Cancer Res. 2007;9:R59. doi: 10.1186/bcr1765. [PMC free article] [PubMed] [Cross Ref]
  • Thomas RP, Guigneaux M, Wood T, Evers BM. Age-associated changes in gene expression patterns in the liver. J Gastrointest Surg. 2002;6:445–453. doi: 10.1016/S1091-255X(01)00010-5. discussion 454. [PubMed] [Cross Ref]
  • Lee JS, Chu IS, Mikaelyan A, Calvisi DF, Heo J, Reddy JK, Thorgeirsson SS. Application of comparative functional genomics to identify best-fit mouse models to study human cancer. Nat Genet. 2004;36:1306–1311. doi: 10.1038/ng1481. [PubMed] [Cross Ref]
  • Perez-Carreon JI, Lopez-Garcia C, Fattel-Fazenda S, Arce-Popoca E, Aleman-Lazarini L, Hernandez-Garcia S, Le Berre V, Sokol S, Francois JM, Villa-Trevino S. Gene expression profile related to the progression of preneoplastic nodules toward hepatocellular carcinoma in rats. Neoplasia. 2006;8:373–383. doi: 10.1593/neo.05841. [PMC free article] [PubMed] [Cross Ref]
  • Paoloni M, Davis S, Lana S, Withrow S, Sangiorgi L, Picci P, Hewitt S, Triche T, Meltzer P, Khanna C. Canine tumor cross-species genomics uncovers targets linked to osteosarcoma progression. BMC Genomics. 2009;10:625. doi: 10.1186/1471-2164-10-625. [PMC free article] [PubMed] [Cross Ref]
  • Sweet-Cordero A, Mukherjee S, Subramanian A, You H, Roix JJ, Ladd-Acosta C, Mesirov J, Golub TR, Jacks T. An oncogenic KRAS2 expression signature identified by cross-species gene-expression analysis. Nat Genet. 2005;37:48–55. [PubMed]
  • Ellwood-Yen K, Graeber TG, Wongvipat J, Iruela-Arispe ML, Zhang J, Matusik R, Thomas GV, Sawyers CL. Myc-driven murine prostate cancer shares molecular features with human prostate tumors. Cancer Cell. 2003;4:223–238. doi: 10.1016/S1535-6108(03)00197-1. [PubMed] [Cross Ref]
  • Lee JS, Thorgeirsson SS. Comparative and integrative functional genomics of HCC. Oncogene. 2006;25:3801–3809. doi: 10.1038/sj.onc.1209561. [PubMed] [Cross Ref]
  • Garraway LA, Widlund HR, Rubin MA, Getz G, Berger AJ, Ramaswamy S, Beroukhim R, Milner DA, Granter SR, Du J. Integrative genomic analyses identify MITF as a lineage survival oncogene amplified in malignant melanoma. Nature. 2005;436:117–122. doi: 10.1038/nature03664. [PubMed] [Cross Ref]
  • Thorgeirsson SS, Lee JS, Grisham JW. Molecular prognostication of liver cancer: end of the beginning. J Hepatol. 2006;44:798–805. doi: 10.1016/j.jhep.2006.01.008. [PubMed] [Cross Ref]
  • Cifola I, Spinelli R, Beltrame L, Peano C, Fasoli E, Ferrero S, Bosari S, Signorini S, Rocco F, Perego R. Genome-wide screening of copy number alterations and LOH events in renal cell carcinomas and integration with gene expression profile. Mol Cancer. 2008;7:6. doi: 10.1186/1476-4598-7-6. [PMC free article] [PubMed] [Cross Ref]
  • Mas VR, Maluf DG, Archer KJ, Yanek K, Kong X, Kulik L, Freise CE, Olthoff KM, Ghobrial RM, McIver P, Fisher R. Genes involved in viral carcinogenesis and tumor initiation in hepatitis C virus-induced hepatocellular carcinoma. Mol Med. 2009;15:85–94. doi: 10.2119/molmed.2008.00110. [PMC free article] [PubMed] [Cross Ref]
  • Pavlidis P, Noble WS. Analysis of strain and regional variation in gene expression in mouse brain. Genome Biol. 2001;2:RESEARCH0042. doi: 10.1186/gb-2001-2-10-research0042. [PMC free article] [PubMed] [Cross Ref]
  • Tsafrir D, Bacolod M, Selvanayagam Z, Tsafrir I, Shia J, Zeng Z, Liu H, Krier C, Stengel RF, Barany F. Relationship of gene expression and chromosomal abnormalities in colorectal cancer. Cancer Res. 2006;66:2129–2137. doi: 10.1158/0008-5472.CAN-05-2569. [PubMed] [Cross Ref]
  • Moinzadeh P, Breuhahn K, Stutzer H, Schirmacher P. Chromosome alterations in human hepatocellular carcinomas correlate with aetiology and histological grade--results of an explorative CGH meta-analysis. Br J Cancer. 2005;92:935–941. doi: 10.1038/sj.bjc.6602448. [PMC free article] [PubMed] [Cross Ref]
  • Albertson DG, Collins C, McCormick F, Gray JW. Chromosome aberrations in solid tumors. Nat Genet. 2003;34:369–376. doi: 10.1038/ng1215. [PubMed] [Cross Ref]
  • Pollack JR, Sorlie T, Perou CM, Rees CA, Jeffrey SS, Lonning PE, Tibshirani R, Botstein D, Borresen-Dale AL, Brown PO. Microarray analysis reveals a major direct role of DNA copy number alteration in the transcriptional program of human breast tumors. Proc Natl Acad Sci USA. 2002;99:12963–12968. doi: 10.1073/pnas.162471999. [PubMed] [Cross Ref]
  • Chiang DY, Villanueva A, Hoshida Y, Peix J, Newell P, Minguez B, LeBlanc AC, Donovan DJ, Thung SN, Sole M. Focal gains of VEGFA and molecular classification of hepatocellular carcinoma. Cancer Res. 2008;68:6779–6788. doi: 10.1158/0008-5472.CAN-08-0742. [PMC free article] [PubMed] [Cross Ref]
  • Boyault S, Rickman DS, de Reynies A, Balabaud C, Rebouissou S, Jeannot E, Herault A, Saric J, Belghiti J, Franco D. Transcriptome classification of HCC is related to gene alterations and to new therapeutic targets. Hepatology. 2007;45:42–52. doi: 10.1002/hep.21467. [PubMed] [Cross Ref]
  • Liao YL, Sun YM, Chau GY, Chau YP, Lai TC, Wang JL, Horng JT, Hsiao M, Tsou AP. Identification of SOX4 target genes using phylogenetic footprinting-based prediction from expression microarrays suggests that overexpression of SOX4 potentiates metastasis in hepatocellular carcinoma. Oncogene. 2008. [PubMed]
  • Ivanova NB, Dimos JT, Schaniel C, Hackney JA, Moore KA, Lemischka IR. A stem cell molecular signature. Science. 2002;298:601–604. doi: 10.1126/science.1073823. [PubMed] [Cross Ref]
  • Dudoit SSP, Boldrick JC. Multiple hypothesis testing in microarray Experiments. Statist Scie. 2003;18:71–103. doi: 10.1214/ss/1056397487. [Cross Ref]
  • Liu YF, Zha BS, Zhang HL, Zhu XJ, Li YH, Zhu J, Guan XH, Feng ZQ, Zhang JP. Characteristic gene expression profiles in the progression from liver cirrhosis to carcinoma induced by diethylnitrosamine in a rat model. J Exp Clin Cancer Res. 2009;28:107. doi: 10.1186/1756-9966-28-107. [PMC free article] [PubMed] [Cross Ref]
  • Thorgeirsson SS, Grisham JW. Molecular pathogenesis of human hepatocellular carcinoma. Nat Genet. 2002;31:339–346. doi: 10.1038/ng0802-339. [PubMed] [Cross Ref]
  • Wanless IR. International consensus on histologic diagnosis of early hepatocellular neoplasia. Hepatol Res. 2007;37(Suppl 2):S139–141. doi: 10.1111/j.1872-034X.2007.00177.x. [PubMed] [Cross Ref]
  • Francois P, Garzoni C, Bento M, Schrenzel J. Comparison of amplification methods for transcriptomic analyses of low abundance prokaryotic RNA sources. J Microbiol Methods. 2007;68:385–391. doi: 10.1016/j.mimet.2006.09.022. [PubMed] [Cross Ref]
  • Kleinbaum DG, Kupper L.L, Morgenstern H. Epidemiological Research- Principles and Quantitative Methods. New York; Confounding; pp. 243–267.
  • Guo L, Lobenhofer EK, Wang C, Shippy R, Harris SC, Zhang L, Mei N, Chen T, Herman D, Goodsaid FM. Rat toxicogenomic study reveals analytical consistency across microarray platforms. Nat Biotechnol. 2006;24:1162–1169. doi: 10.1038/nbt1238. [PubMed] [Cross Ref]
  • Shi L, Reid LH, Jones WD, Shippy R, Warrington JA, Baker SC, Collins PJ, de Longueville F, Kawasaki ES, Lee KY. The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. Nat Biotechnol. 2006;24:1151–1161. doi: 10.1038/nbt1239. [PMC free article] [PubMed] [Cross Ref]
  • Ghosh D, Barette TR, Rhodes D, Chinnaiyan AM. Statistical issues and methods for meta-analysis of microarray data: a case study in prostate cancer. Funct Integr Genomics. 2003;3:180–188. doi: 10.1007/s10142-003-0087-5. [PubMed] [Cross Ref]
  • Panteva M, Korkaya H, Jameel S. Hepatitis viruses and the MAPK pathway: is this a survival strategy? Virus Res. 2003;92:131–140. doi: 10.1016/S0168-1702(02)00356-8. [PubMed] [Cross Ref]
  • Wu CH, Sahoo D, Arvanitis C, Bradon N, Dill DL, Felsher DW. Combined analysis of murine and human microarrays and ChIP analysis reveals genes associated with the ability of MYC to maintain tumorigenesis. PLoS Genet. 2008;4:e1000090. doi: 10.1371/journal.pgen.1000090. [PMC free article] [PubMed] [Cross Ref]
  • Chang Q, Chen J, Beezhold KJ, Castranova V, Shi X, Chen F. JNK1 activation predicts the prognostic outcome of the human hepatocellular carcinoma. Mol Cancer. 2009;8:64. doi: 10.1186/1476-4598-8-64. [PMC free article] [PubMed] [Cross Ref]
  • Liu J, Ahiekpor A, Li L, Li X, Arbuthnot P, Kew M, Feitelson MA. Increased expression of ErbB-2 in liver is associated with hepatitis B × antigen and shorter survival in patients with liver cancer. Int J Cancer. 2009;125:1894–1901. doi: 10.1002/ijc.24580. [PMC free article] [PubMed] [Cross Ref]
  • Fusco O, Querzoli P, Nenci I, Natoli C, Brakebush C, Ullrich A, Iacobelli S. 90K (MAC-2 BP) gene expression in breast cancer and evidence for the production of 90 K by peripheral-blood mononuclear cells. Int J Cancer. 1998;79:23–26. doi: 10.1002/(SICI)1097-0215(19980220)79:1<23::AID-IJC5>3.0.CO;2-Y. [PubMed] [Cross Ref]
  • Kim SJ, Lee SJ, Sung HJ, Choi IK, Choi CW, Kim BS, Kim JS, Yu W, Hwang HS, Kim IS. Increased serum 90 K and Galectin-3 expression are associated with advanced stage and a worse prognosis in diffuse large B-cell lymphomas. Acta Haematol. 2008;120:211–216. doi: 10.1159/000193223. [PubMed] [Cross Ref]
  • Vasko V, Espinosa AV, Scouten W, He H, Auer H, Liyanarachchi S, Larin A, Savchenko V, Francis GL, de la Chapelle A. Gene expression and functional evidence of epithelial-to-mesenchymal transition in papillary thyroid carcinoma invasion. Proc Natl Acad Sci USA. 2007;104:2803–2808. doi: 10.1073/pnas.0610733104. [PubMed] [Cross Ref]
  • Yamada S, Ohira M, Horie H, Ando K, Takayasu H, Suzuki Y, Sugano S, Hirata T, Goto T, Matsunaga T. Expression profiling and differential screening between hepatoblastomas and the corresponding normal livers: identification of high expression of the PLK1 oncogene as a poor-prognostic indicator of hepatoblastomas. Oncogene. 2004;23:5901–5911. doi: 10.1038/sj.onc.1207782. [PubMed] [Cross Ref]
  • Ma W, Cai S, Du J, Tan Y, Chen H, Guo Z, Hu H, Fang R, Cai S. SDF-1/54-DCN: a novel recombinant chimera with dual inhibitory effects on proliferation and chemotaxis of tumor cells. Biol Pharm Bull. 2008;31:1086–1091. doi: 10.1248/bpb.31.1086. [PubMed] [Cross Ref]
  • Bi X, Tong C, Dockendorff A, Bancroft L, Gallagher L, Guzman G, Iozzo RV, Augenlicht LH, Yang W. Genetic deficiency of decorin causes intestinal tumor formation through disruption of intestinal cell maturation. Carcinogenesis. 2008;29:1435–1440. doi: 10.1093/carcin/bgn141. [PMC free article] [PubMed] [Cross Ref]
  • Luo SM, Tan WM, Deng WX, Zhuang SM, Luo JW. Expression of albumin, IGF-1, IGFBP-3 in tumor tissues and adjacent non-tumor tissues of hepatocellular carcinoma patients with cirrhosis. World J Gastroenterol. 2005;11:4272–4276. [PubMed]
  • Sato T, Oshima T, Yoshihara K, Yamamoto N, Yamada R, Nagano Y, Fujii S, Kunisaki C, Shiozawa M, Akaike M. Overexpression of the fibroblast growth factor receptor-1 gene correlates with liver metastasis in colorectal cancer. Oncol Rep. 2009;21:211–216. [PubMed]
  • Pollmann MA, Shao Q, Laird DW, Sandig M. Connexin 43 mediated gap junctional communication enhances breast tumor cell diapedesis in culture. Breast Cancer Res. 2005;7:R522–534. doi: 10.1186/bcr1042. [PMC free article] [PubMed] [Cross Ref]
  • Bea S, Salaverria I, Armengol L, Pinyol M, Fernandez V, Hartmann EM, Jares P, Amador V, Hernandez L, Navarro A. Uniparental disomies, homozygous deletions, amplifications, and target genes in mantle cell lymphoma revealed by integrative high-resolution whole-genome profiling. Blood. 2009;113:3059–3069. doi: 10.1182/blood-2008-07-170183. [PubMed] [Cross Ref]
  • Kajiyama H, Kikkawa F, Suzuki T, Shibata K, Ino K, Mizutani S. Prolonged survival and decreased invasive activity attributable to dipeptidyl peptidase IV overexpression in ovarian carcinoma. Cancer Res. 2002;62:2753–2757. [PubMed]
  • Zhu L, Gibson P, Currle DS, Tong Y, Richardson RJ, Bayazitov IT, Poppleton H, Zakharenko S, Ellison DW, Gilbertson RJ. Prominin 1 marks intestinal stem cells that are susceptible to neoplastic transformation. Nature. 2009;457:603–607. doi: 10.1038/nature07589. [PMC free article] [PubMed] [Cross Ref]
  • Ma O, Cai WW, Zender L, Dayaram T, Shen J, Herron AJ, Lowe SW, Man TK, Lau CC, Donehower LA. MMP13, Birc2 (cIAP1), and Birc3 (cIAP2), amplified on chromosome 9, collaborate with p53 deficiency in mouse osteosarcoma progression. Cancer Res. 2009;69:2559–2567. doi: 10.1158/0008-5472.CAN-08-2929. [PMC free article] [PubMed] [Cross Ref]
  • Huang W, Yu LF, Zhong J, Wu W, Zhu JY, Jiang FX, Wu YL. Stat3 is involved in angiotensin II-induced expression of MMP2 in gastric cancer cells. Dig Dis Sci. 2009;54:2056–2062. doi: 10.1007/s10620-008-0617-z. [PubMed] [Cross Ref]
  • Ibanez de Caceres I, Dulaimi E, Hoffman AM, Al-Saleem T, Uzzo RG, Cairns P. Identification of novel target genes by an epigenetic reactivation screen of renal cancer. Cancer Res. 2006;66:5021–5028. doi: 10.1158/0008-5472.CAN-05-3365. [PubMed] [Cross Ref]
  • Chen C, Mendez E, Houck J, Fan W, Lohavanichbutr P, Doody D, Yueh B, Futran ND, Upton M, Farwell DG. Gene expression profiling identifies genes predictive of oral squamous cell carcinoma. Cancer Epidemiol Biomarkers Prev. 2008;17:2152–2162. doi: 10.1158/1055-9965.EPI-07-2893. [PMC free article] [PubMed] [Cross Ref]
  • Rubinstein R, Simon I. MILANO--custom annotation of microarray results using automatic literature searches. BMC Bioinformatics. 2005;6:12. doi: 10.1186/1471-2105-6-12. [PMC free article] [PubMed] [Cross Ref]
  • Patil MA, Chua MS, Pan KH, Lin R, Lih CJ, Cheung ST, Ho C, Li R, Fan ST, Cohen SN. An integrated data analysis approach to characterize genes highly expressed in hepatocellular carcinoma. Oncogene. 2005;24:3737–3747. doi: 10.1038/sj.onc.1208479. [PubMed] [Cross Ref]
  • Zhao X, Weir BA, LaFramboise T, Lin M, Beroukhim R, Garraway L, Beheshti J, Lee JC, Naoki K, Richards WG. Homozygous deletions and chromosome amplifications in human lung carcinomas revealed by single nucleotide polymorphism array analysis. Cancer Res. 2005;65:5561–5570. doi: 10.1158/0008-5472.CAN-04-4603. [PubMed] [Cross Ref]
  • Marchetti A, Tinari N, Buttitta F, Chella A, Angeletti CA, Sacco R, Mucilli F, Ullrich A, Iacobelli S. Expression of 90 K (Mac-2 BP) correlates with distant metastasis and predicts survival in stage I non-small cell lung cancer patients. Cancer Res. 2002;62:2535–2539. [PubMed]
  • Kholova I, Ludvikova M, Ryska A, Topolcan O, Pikner R, Pecen L, Cap J, Holubec L Jr. Diagnostic role of markers dipeptidyl peptidase IV and thyroid peroxidase in thyroid tumors. Anticancer Res. 2003;23:871–875. [PubMed]
  • Halbersztadt A, Halon A, Pajak J, Robaczynski J, Rabczynski J, St Gabrys M. [The role of matrix metalloproteinases in tumor invasion and metastasis] Ginekol Pol. 2006;77:63–71. [PubMed]
  • Samantaray S, Sharma R, Chattopadhyaya TK, Gupta SD, Ralhan R. Increased expression of MMP-2 and MMP-9 in esophageal squamous cell carcinoma. J Cancer Res Clin Oncol. 2004;130:37–44. doi: 10.1007/s00432-003-0500-4. [PubMed] [Cross Ref]
  • McKenna GJ, Chen Y, Smith RM, Meneghetti A, Ong C, McMaster R, Scudamore CH, Chung SW. A role for matrix metalloproteinases and tumor host interaction in hepatocellular carcinomas. Am J Surg. 2002;183:588–594. doi: 10.1016/S0002-9610(02)00833-4. [PubMed] [Cross Ref]
  • Asamoto M, Hokaiwado N, Cho Y-M, Shirai T. Effects of genetic background on prostate and taste bud carcinogenesis due to SV40 T antigen expression under probasin gene promoter control. Carcinogenesis. 2002;23:463–467. doi: 10.1093/carcin/23.3.463. [PubMed] [Cross Ref]
  • Yamashita S, Wakazono K, Nomoto T, Tsujino Y, Kuramoto T, Ushijima T. Expression Quantitative Trait Loci Analysis of 13 Genes in the Rat Prostate. Genetics. 2005;171:1231–1238. doi: 10.1534/genetics.104.038174. [PubMed] [Cross Ref]
  • Dillner K, Kindblom J, Flores-Morales A, Shao R, Tornell J, Norstedt G, Wennbo H. Gene Expression Analysis of Prostate Hyperplasia in Mice Overexpressing the Prolactin Gene Specifically in the Prostate. Endocrinology. 2003;144:4955–4966. doi: 10.1210/en.2003-0415. [PubMed] [Cross Ref]
  • Maffey AH, Ishibashi T, He C, Wang X, White AR, Hendy SC, Nelson CC, Rennie PS, Ausio J. Probasin promoter assembles into a strongly positioned nucleosome that permits androgen receptor binding. Mol Cell Endocrinol. 2007;268:10–19. doi: 10.1016/j.mce.2007.01.009. [PubMed] [Cross Ref]
  • Nikitovic D, Berdiaki A, Zafiropoulos A, Katonis P, Tsatsakis A, Karamanos NK, Tzanakakis GN. Lumican expression is positively correlated with the differentiation and negatively with the growth of human osteosarcoma cells. Febs J. 2008;275:350–361. [PubMed]
  • Eshchenko TY, Rykova VI, Chernakov AE, Sidorov SV, Grigorieva EV. Expression of different proteoglycans in human breast tumors. Biochemistry (Mosc) 2007;72:1016–1020. doi: 10.1134/S0006297907090143. [PubMed] [Cross Ref]
  • Ishiwata T, Cho K, Kawahara K, Yamamoto T, Fujiwara Y, Uchida E, Tajiri T, Naito Z. Role of lumican in cancer cells and adjacent stromal tissues in human pancreatic cancer. Oncol Rep. 2007;18:537–543. [PubMed]
  • Seya T, Tanaka N, Shinji S, Yokoi K, Koizumi M, Teranishi N, Yamashita K, Tajiri T, Ishiwata T, Naito Z. Lumican expression in advanced colorectal cancer with nodal metastasis correlates with poor prognosis. Oncol Rep. 2006;16:1225–1230. [PubMed]
  • Kelemen LE, Couch FJ, Ahmed S, Dunning AM, Pharoah PD, Easton DF, Fredericksen ZS, Vierkant RA, Pankratz VS, Goode EL. Genetic variation in stromal proteins decorin and lumican with breast cancer: investigations in two case-control studies. Breast Cancer Res. 2008;10:R98. doi: 10.1186/bcr2201. [PMC free article] [PubMed] [Cross Ref]
  • Yee J, Sadar MD, Sin DD, Kuzyk M, Xing L, Kondra J, McWilliams A, Man SF, Lam S. Connective tissue-activating peptide III: a novel blood biomarker for early lung cancer detection. J Clin Oncol. 2009;27:2787–2792. doi: 10.1200/JCO.2008.19.4233. [PMC free article] [PubMed] [Cross Ref]
  • Marshall KW, Mohr S, Khettabi FE, Nossova N, Chao S, Bao W, Ma J, Li XJ, Liew CC. A blood-based biomarker panel for stratifying current risk for colorectal cancer. Int J Cancer. 2010;126:1177–1186. [PubMed]
  • Lonneborg A. Biomarkers for Alzheimer disease in cerebrospinal fluid, urine, and blood. Mol Diagn Ther. 2008;12:307–320. [PubMed]
  • Aaroe J, Lindahl T, Dumeaux V, Saebo S, Tobin D, Hagen N, Skaane P, Lonneborg A, Sharma P, Borresen-Dale AL. Gene expression profiling of peripheral blood cells for early detection of breast cancer. Breast Cancer Res. 2010;12:R7. doi: 10.1186/bcr2472. [PMC free article] [PubMed] [Cross Ref]
  • Kurian SM, Le-Niculescu H, Patel SD, Bertram D, Davis J, Dike C, Yehyawi N, Lysaker P, Dustin J, Caligiuri M. Identification of blood biomarkers for psychosis using convergent functional genomics. Mol Psychiatry. 2009. in press . [PubMed]
  • Lonneborg A, Aaroe J, Dumeaux V, Borresen-Dale AL. Found in transcription: gene expression and other novel blood biomarkers for the early detection of breast cancer. Expert Rev Anticancer Ther. 2009;9:1115–1123. doi: 10.1586/era.09.31. [PubMed] [Cross Ref]
  • Borovecki F, Lovrecic L, Zhou J, Jeong H, Then F, Rosas HD, Hersch SM, Hogarth P, Bouzou B, Jensen RV, Krainc D. Genome-wide expression profiling of human blood reveals biomarkers for Huntington's disease. Proc Natl Acad Sci USA. 2005;102:11023–11028. doi: 10.1073/pnas.0504921102. [PubMed] [Cross Ref]
  • Terminology of nodular hepatocellular lesions. International Working Party. Hepatology. 1995;22:983–993. [PubMed]
  • Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80. doi: 10.1186/gb-2004-5-10-r80. [PMC free article] [PubMed] [Cross Ref]
  • Li C, Hung Wong W. Model-based analysis of oligonucleotide arrays: model validation, design issues and standard error application. Genome Biol. 2001;2:RESEARCH0032. [PMC free article] [PubMed]
  • Wu Z, Irizarry RA. Stochastic models inspired by hybridization theory for short oligonucleotide arrays. J Comput Biol. 2005;12:882–893. doi: 10.1089/cmb.2005.12.882. [PubMed] [Cross Ref]
  • Wu Z, Irizarry RA. Preprocessing of oligonucleotide array data. Nat Biotechnol. 2004;22:656–658. doi: 10.1038/nbt0604-656b. author reply 658. [PubMed] [Cross Ref]
  • Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M. TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003;34:374–378. [PubMed]
  • Sherman BT, Huang DW, Tan Q, Guo Y, Bour S, Liu D, Stephens R, Baseler MW, Lane HC, Lempicki RA. DAVID Knowledgebase: a gene-centered database integrating heterogeneous gene annotation resources to facilitate high-throughput gene functional analysis. BMC Bioinformatics. 2007;8:426. doi: 10.1186/1471-2105-8-426. [PMC free article] [PubMed] [Cross Ref]
  • Hosack DA, Dennis G Jr, Sherman BT, Lane HC, Lempicki RA. Identifying biological themes within lists of genes with EASE. Genome Biol. 2003;4:R70. doi: 10.1186/gb-2003-4-10-r70. [PMC free article] [PubMed] [Cross Ref]
  • Thomas PD, Campbell MJ, Kejariwal A, Mi H, Karlak B, Daverman R, Diemer K, Muruganujan A, Narechania A. PANTHER: a library of protein families and subfamilies indexed by function. Genome Res. 2003;13:2129–2141. doi: 10.1101/gr.772403. [PubMed] [Cross Ref]
  • Thomas PD, Kejariwal A, Guo N, Mi H, Campbell MJ, Muruganujan A, Lazareva-Ulitsky B. Applications for protein sequence-function evolution data: mRNA/protein expression analysis and coding SNP scoring tools. Nucleic Acids Res. 2006;34:W645–650. doi: 10.1093/nar/gkl229. [PMC free article] [PubMed] [Cross Ref]
  • Tsai J, Sultana R, Lee Y, Pertea G, Karamycheva S, Antonescu V, Cho J, Parvizi B, Cheung F, Quackenbush J. RESOURCERER: a database for annotating and linking microarray resources within and across species. Genome Biol. 2001;2:SOFTWARE0002. doi: 10.1186/gb-2001-2-11-software0002. [PMC free article] [PubMed] [Cross Ref]
  • Reiner A, Yekutieli D, Benjamini Y. Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics. 2003;19:368–375. doi: 10.1093/bioinformatics/btf877. [PubMed] [Cross Ref]
  • Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402–408. doi: 10.1006/meth.2001.1262. [PubMed] [Cross Ref]

Articles from Molecular Cancer are provided here courtesy of BioMed Central