- Open Access
Characteristic gene expression profiles in the progression from liver cirrhosis to carcinoma induced by diethylnitrosamine in a rat model
Journal of Experimental & Clinical Cancer Research volume 28, Article number: 107 (2009)
Liver cancr is a heterogeneous disease in terms of etiology, biologic and clinical behavior. Very little is known about how many genes concur at the molecular level of tumor development, progression and aggressiveness. To explore the key genes involved in the development of liver cancer, we established a rat model induced by diethylnitrosamine to investigate the gene expression profiles of liver tissues during the transition to cirrhosis and carcinoma.
A rat model of liver cancer induced by diethylnitrosamine was established. The cirrhotic tissue, the dysplasia nodules, the early cancerous nodules and the cancerous nodules from the rats with lung metastasis were chosen to compare with liver tissue of normal rats to investigate the differential expression genes between them. Affymetrix GeneChip Rat 230 2.0 arrays were used throughout. The real-time quantity PCR was used to verify the expression of some differential expression genes in tissues.
The pathological changes that occurred in the livers of diethylnitrosamine-treated rats included non-specific injury, fibrosis and cirrhosis, dysplastic nodules, early cancerous nodules and metastasis. There are 349 upregulated and 345 downregulated genes sharing among the above chosen tissues when compared with liver tissue of normal rats. The deregulated genes play various roles in diverse processes such as metabolism, transport, cell proliferation, apoptosis, cell adhesion, angiogenesis and so on. Among which, 41 upregulated and 27 downregulated genes are associated with inflammatory response, immune response and oxidative stress. Twenty-four genes associated with glutathione metabolism majorly participating oxidative stress were deregulated in the development of liver cancer. There were 19 members belong to CYP450 family downregulated, except CYP2C40 upregulated.
In this study, we provide the global gene expression profiles during the development and progression of liver cancer in rats. The data obtained from the gene expression profiles will allow us to acquire insights into the molecular mechanisms of hepatocarcinogenesis and identify specific genes (or gene products) that can be used for early molecular diagnosis, risk analysis, prognosis prediction, and development of new therapies.
HCC is a heterogeneous disease in terms of etiology, biologic and clinical behavior. Meanwhile, hepatocarcinogenesis is a long-term, multistep process associated with changes in gene expression profiles. In the last several years, there have been important gains in our understanding of the pathogenesis of HCC and our appreciation of the critical oncogenic and tumor suppressor pathways involved in hepatocarcinogenesis [1–5]. Despite this, current knowledge about the molecular pathogenesis of HCC is a result of investigations of fully developed HCC. Very little is known about how many genes concur at the molecular level of tumor development, progression and aggressiveness.
Molecular profiling has been successfully used to identify candidate genes for HCC in human and animal model systems. Although many approaches (including genome-scale studies) provide insights into some of the stages in human tumorigenesis, a sequential analysis of the development of tumors in humans is very difficult. Most of them have not given us the gene expression profiles that could point to those genes that play key roles during the whole carcinogenetic process from initiation to metastasis. Animal models of carcinogenesis have permitted the examination of the stages of neoplastic development in considerable detail.
In this study, we established the rat model of liver cancer induced by DEN to explore the processes of initiation and progression of HCC. HCC develops commonly, but not exclusively, in a setting of chronic liver cell injury, which leads to inflammation, hepatocyte regeneration, liver matrix remodeling, fibrosis, and ultimately, cirrhosis[7, 8]. The histological changes in DEN-induced liver cancer in rats are similar to those seen in human HCC. We think the similar phenotype might be based on similar gene expression profiles. Affymetrix GeneChip Rat 230 2.0 arrays were used to analyze gene expression profiles of liver tissues from control and DEN-treated rats during the process from cirrhosis to metastasis. This allowed us to obtain an almost complete picture of the early genetic alterations that are directly or indirectly involved in the development of HCC. We supposed that the genes expression profiles deregulated during the process from liver cirrhosis to carcinoma and metastasis play key roles in the hepatocarcinogenesis. The data obtained from the gene expression profiles will allow us to acquire insights into the molecular mechanisms of hepatocarcinogenesis and identify specific genes (or gene products) that can be used for early molecular diagnosis, risk analysis, prognosis prediction, and development of new therapies.
Materials and methods
Animals and treatments
Male Wistar rats weighing 120–150 g at the beginning of the experiments were obtained from SLAC Laboratory Animal Co. Ltd. (Shanghai). The animals were acclimatized to standard laboratory conditions (temperature 22–25°C, relative humidity 50–60%, and 12 hour photoperiods (lights on 07:00–19:00 h)) and were housed in stainless steel wire-mesh cages (five rats per cage). During the entire period of study, the rats were supplied with a semi-purified basal diet (Shanghai) and water ad lib. All experiments followed the Guide for the Care and Use of Laboratory Animals.
Briefly, ninety Wistar rats were randomly divided into two groups: the control and DEN group (DEN, Sigma Chemical Co. St Louis, MO; CAS 56-23-5; purity > 98%). After one week on a basal diet, rats belonging to the DEN group underwent intragastric administration of 1% aqueous solution of DEN (70 mg/kg) once a week, consecutively for 14 weeks. Animals that belonged to the control group received distilled water. There were ten rats in the control group. Daily food and water intakes were noted and the body weights of the animals from each group were recorded every second day. The rats were sacrificed with 25% (g/ml) urethane anesthesia (6 ml/kg) by intraperitoneal injection and sacrificed at different time points. At the end of the 2nd, 4th and 6th week after DEN-treatment, 3, 3, and 4 rats were sacrificed respectively at these time periods. From the end of the 8th week until the end of the 20th week after DEN-treatment, ten rats were sacrificed respectively every two weeks. Meanwhile one age-matched control rat was sacrificed. All the age-matched rats from the control group and rats of treatment groups were sacrificed on the same day.
Sample collection and histopathological analyses
Animals were chosen sequentially for necropsy. The entire liver was observed grossly and weighted. The selected liver tissues were observed for gross changes, divided into pieces of about 0.1 g, snap-frozen directly in liquid nitrogen and stored at -80°C prior to RNA isolation for microarray analysis. The remaining livers were preserved in 10% phosphate-buffered formalin. The liver tissue fixed in neutral formalin was embedded in paraffin, sectioned, and stained with hematoxylin and eosin (H&E). Histopathologic examinations of the liver sections were conducted by a pathologist and peer-reviewed.
Frozen liver sections were ground in a Mixer Mill mm 200 (Retsch GmbH and Co. KG, Haan, Germany) using pre-cooled stainless steel balls. Total RNAs were isolated from livers with Trizol Reagent (Invitrogen, CA) using manufacturer recommended procedures. The ratio of the optical densities from RNA samples measured at 260 and 280 nm was used to evaluate nucleic acid purity, and total RNA concentrations were determined by the absorbance at 260 nm. The quality of total RNA was estimated based on the integrity of 28S and 18S rRNA. RNA was separated using 1% agarose gel electrophoresis. Good RNA quality was indicated by 28S rRNA banding having twice the intensity of the 18S rRNA, without significant smearing of the rRNA bands. Samples of total RNA from livers of rats from the same time points were pooled for subsequent use in the GeneChip analysis. Prior to GeneChip analysis, the pooled total RNA samples were purified using the RNeasy Total RNA Mini Kit (Qiagen, Valencia, CA) according to manufacturer's instructions.
Oligo microarray hybridization
Biotin-labeled cRNA samples were used for hybridization of Affymetrix GeneChip Rat 230 2.0 arrays. The arrays were prepared according to the protocol supplied with the GeneChip Sample Cleanup module (P/N 900371, Affymetrix Inc., Santa Clara, CA). Briefly, 5 μg total RNA was used for cDNA synthesis with the SuperScript Choice System (Invitrogen Life Technologies, Carlsbad, CA) employing a T7-(d7)24 primer. After spin column purification, biotin-labeled cRNA was synthesized from the cDNA using the ENZO RNA Transcript Labeling Kit (Affymetrix Inc.). Spin column-purified cRNA was quality controlled using an Agilent 2100 Bioanalyzer and spectrophotometrically quantified. The cRNA (15 μg) was then fragmented in buffer supplied with the Cleanup Module and hybridized for 16 h at 45°C (Affymetrix Genechip Hybridization Oven 640). The microarrays were washed and stained with streptavidin-phycoerythrin (SAPE, Molecular Probes) on the Affymetrix Fluidics Station 450, including an amplification step according to the manufacturer's instructions. Fluorescent images were read using the Gene Array Scanner 3000. The raw data image files (DAT) were converted into RPT files using Affymetrix Microarray Suite (MAS) 5.0. In RPT files, the scan data from the 36 pixels per oligo set were averaged. Pooled RNAs of cirrhotic tissue, the dysplasia nodules, the early cancerous nodules and the cancerous nodules from the rats with lung metastases at the 12th, 14th, 16th, and 20th week after DEN-treatment and from normal control rat were processed, hybridized and analyzed individually.
The chips were scanned in an Agilent ChipScanner to detect hybridization signals. Average target intensity was set at 500 arbitrary units. Each array was assessed for quality and stability by examining replicated copies of the same gene at different locations on the array. To ensure the quality of the cRNA samples and of the Affymetrix GeneChips, quality control experiments were performed using test chips, and the same cRNA sample used in both the test chip and GeneChip.
Microarray quality control
With GeneChip Operating Software (GCOS) v1.2, dark and white spots, gradients and distortions were detected and corrected using the RPT file data. The GeneChip Rat Genome 230 2.0 Array provides the entire transcribed rat genome on a single array and enables scientists to obtain the most comprehensive view of the transcribed rat genome in order to make accurate biological conclusions. The Affymetrix Rat Genome 230 2.0 microarrays contain 31,000 probe sets corresponding to about 24,000 annotated rat genes and 6693 expressed sequence tags (ESTs). Each probe set is represented by 11–20 pairs of 25 mer oligonucleotides. Each probe pair consists of a perfect match oligo (PM) complementary to the cRNA target sequence and a mismatched oligo (MM). Using the MAS 5.0 statistical algorithms implemented in the Quality Controller software, the intensities of all 11–20 probe pairs were condensed to one intensity value per probe set associated with a statistical detection p value calculated from the intensity differences of the PM and corresponding MM oligos. This p value indicates how reliably a transcript is detected. Transcripts with p < 0.04 were designated present, whereas those with a p > 0.06 were designated absent. Transcripts with 0.04 < p < 0.06 were designated marginal, whose reliability were doubted and need to be verified by methods with higher sensitivity. After condensing (which also included overall microarray background correction) the microarray was scaled to an average signal intensity of 100, after excluding the highest and lowest 2% of the data. GeneChip Rat Genome 230 2.0 microarrays include a set of rat maintenance genes to facilitate the normalization and scaling of array experiments. These probe sets serve as a tool to normalize or scale data prior to performing data comparison. These normalization genes show consistent levels of expression over defined sample sets.
Microarray data analysis
The microarray data were analyzed using the microarray suite 5.0 software. First, using the present genes, those significantly deregulated between the DEN-treated and control groups were selected using a two-sample t-test with a p cutoff value of 0.001 in combination with n-fold regulation/ratio of means. This requires that the mean intensity of a gene has to be at least n times different between the treated and control sample group to be included. The cutoff was set at 2 times.
Secondly, genes designated present in treated samples but absent in controls, or vice versa, were determined, as these could be genes induced from or repressed to background expression levels, respectively, after treatment. From these genes, those discriminating between treated and control samples were again selected with a two-sample t-test (p < 0.001), combined with the requirement of an at least two-fold difference of the mean intensities for a given gene.
Scatter plot, gene tree
Scatter plots were used to visually examine the expressional level of genes between the control and DEN-exposed groups. Hierarchical dendrograms were drawn with the Cluster (2.0). It was created by clustering the genes according to their expression in response to the carcinogenic agent. Genes sharing similar expression profiles tended to be clustered together, and the location of a branch containing the genes can be considered a measure of how similar the gene expression was. Genes were selected for the construction of gene tree if the expression of the gene was two-fold greater or less in the treatments, relative to that in the corresponding control. The horizontal axis shows the clustering of the genes according to their expression across treatments; while the vertical axis showed the clustering according to their expression profile in the treatment.
The genechip probe array system only allows comparison of one treatment hybridizing with the probe set. In a comparison analysis, two samples were hybridized to two genechip probe arrays of the same type, they were compared against each other in order to detect and quantify changes in gene expression. One genechip was for baseline (control) and the other was for the experiment (treatment). Two sets of algorithms were generated and they were used to generate change significance and change quantity metrics for every probe set using Microarray Suite (MAS) version 5.0 (Affymetrix, CA). The change algorithm generated a Change p value and an associated fold-change value. The second algorithm gave a quantitative estimate of the change in gene expression in the form of Signal Log Ratio. In the present study, the level of gene expression can be regarded as increased if its Change p-value was less than 0.002 and the gene expression would be considered to be decreased if its Change p-value was greater than 0.997. This method has been used by other investigators. Fold change could be calculated with the following formula: fold change = 2(signal log ratio).
Validation of differential expression of genes by real-time RT-PCR
The differential expression of selected genes was further validated by real-time PCR with SYBR green-based detection (ABI) using gene-specific primer pairs that were run on an ABI 7000 fluorescent sequence detection system (Perkin-Elmer, Foster City, CA). It is known that only the cDNA of each target gene can be the template, owing to the fact all the primer sequences span one or more introns by design (Gene Runner). The rat housekeeping gene β-actin was used as the control. Quantitative values were obtained from the cycle number (Ct value) at which the increase in fluorescent signal (associated with exponential growth of PCR products) starts to be picked up by the laser detector of the detection system. Results, expressed as N-fold differences in target gene expression between the liver tissues of DEN-treated and normal rats and termed 'Ntarget' were determined using the formula: Ntarget = 2ΔCtsample (while ΔCtsample = ΔCtDEN - ΔCtNormal), where the ΔCtDEN and ΔCtnormal values of the sample were determined by subtracting the Ct value of the target gene from the average Ct value of the β-actin gene.
The histological changes of livers of the DEN-treated rats can be divided into three stages. Initially, from the 2nd to 8th week, non-specific injury occurred such as cellular swelling, fatty changes, necrosis, inflammatory infiltration and hepatocyte regeneration. On the 10th to the 14th week, significant liver fibrosis occurred. At the 10th week, the livers showed an quantitative increase in connective tissue, and encapsulation of regenerative nodules, while at the end of the 12th week, nodular cirrhosis could be seen macroscopically. At the 14th week, gray-white nodules, 3 mm to 5 mm in diameter, could be distinguished from the surrounding reddish brown cirrhosis nodules in the livers of 2/10 rats. These were histologically diagnosed as dysplastic nodules. From the 16th to the 20th week the number of nodules increased significantly. At the 16th week, nodules, 5 mm to 1.5 cm in diameter, could be distinguished in the livers of 8/10 rats, while at the 18th and the 20th week, gray-white nodules were present in the livers of all 20 rats. In addition, by the 20th week, abdominal cavity and lung metastases were observed in 2/10 rats. (Figure 1, 2)
Isolation of RNA, microarray hybridization, scanning and signal detection
Liver tissues were taken from the control and the DEN-treated rats. Total RNA was isolated from theses samples and used to prepare cRNA probes for hybridization with Affymetrix GeneChip Rat 230 2.0 arrays (Figure 3). The hybridized microarrays were then scanned and the signals acquired (Figure 4). At the 12th week, liver cirrhosis occurred in 10 of 10 rats, so we took the pooled cirrhotic tissues from the 10 rats for the microarrays. At the 14th week, dysplastic nodules occurred only in the livers of 2/10 rats, so we took the pooled dysplastic nodules from the two rats for the microarrays. At the 16th week, early tumor nodules occurred in the liver of 8/10 rats, so we took the pooled tumor nodules from the eight rats for the microarrays. At the 20th week, tumor nodules occurred in all of the ten rats(10/10), but lung metastasis only occurred in the two of them, so we took the pooled tumor nodules in the liver from the two rats with lung metastasis for the microarrays. We used the pooled liver tissues from the control rats killed at the 12th, 14th, 16th and the 20th week for the microarrays. The decision to pool the mRNA from the rat livers was made in order to obtain a representative analysis of gene expression changes across more than one animal.
Analysis of the differential expression genes
The differential expression genes of cirrhotic tissue, dysplastic nodules, early tumors nodules and tumor nodules from rats with lung metastasis compared with the tissue from normal rats were screened and to determine the upregulated and downregulated DEGs. The results are shown in Table 1.
The upregulated and the downregulated DEGs of cirrhosis tissue, dysplastic nodules, early tumor nodules and tumor nodules with lung metastasis vs the normal livers respectively were grouped into functional categories using gene ontology (GO) terms (Table 2). The GO terms such as metabolism, transport, cellular proliferation, apoptosis, adhesion, angiogenesis, etc. were chosen. Meanwhile, some other genes were associated with oxidative stress, immune response and inflammatory response.
The histological changes during the hepatocarnogenesis in DEN-treated rat models were similar to those seen in humans, including non-specific damage, fibrosis, cirrhosis, dysplastic nodules, early tumorous nodules, progression and metastasis, which appeared to be sequential events. The processes of chronic inflammation, fibrosis and cirrhosis are closely related to liver cancer, while cirrhosis was considered as the precancerous lesions. Therefore, the co-expression of deregulated genes among these four stages might suggest they play key roles in the development of hepatocellular carcinoma. Among upregulated DEGs sharing from cirrhosis to metastasis, there were 246 known genes, 39 translocation loci, 51 inferred genes and 13 unkown genes; while among downregulated DEGs sharing from cirrhosis to metastasis, there were 215 known genes, 48 translocation loci, 63 inferred genes and 19 unkown genes (see additional file 1).
Cellular proliferation, apoptosis, adhesion, migration and agiogenesis all play important roles in carcinogenesis. Extracellular matrix (ECM) is important for the growth and progression of hepatocarcinogenesis (see additional file 2 and Tables 3). There are 105 upregulated and 51 downregulated DEGs with the above functions.
The rat models of liver cancer induced by DEN occurred following chronic injury, regenetation, fiborsis and cirrhosis. Elements of the inflammatory response, immune response and oxidative stress were also involved in the process of hepatocarcinogenesis. Tables 4 and 5 show that the expression of 40 such genes was upregulated and the expression of 27 genes was downregulated.
Hierarchy cluster analysis of DEGs
The DEGs commonly present from stages of cirrhosis to metastasis were clustered by Cluster software according to the values of log2 ratio(Figure 5A,B,C). The upregulated (red) and the downregulated (green) DEGs had the same alterated tendency during the process from liver cirrhosis to metastasis. Furthermore, the DEGs involved in the metabolism of glucose, lipids and alcohol and so on (Figure 6A), DEGs associated with the metabolism of glutathione (Figure 6B) and DEGs of members belong to the CYPs family were listed (Figure 6C).
Validation of differential expression of genes by real-time RT-PCR
The DEGs detected through Affymetrix genechip analysis were confirmed in the selected tissue of DEN-treated and control rats by real-time RT-PCR, as shown in Figure 7. TWEAKR, ANXA2, CTGF were chosen from the upregulated DEGs, and EGFR, KDR, CXCL12 were chosen from the downregulated DEGs. The primer sequences for each gene were listed in Table 6. The quality and specificity of the amplified products were confirmed by visualization on a 2% agarose gel. The results confirmed the validity of the Affymetrix genechip results. The lower the ΔCt value of the target gene, the more mRNA content of the target gene there is in the tissue. The Ct value of β-actin obtained from DEN-treated and normal tissue was almost identical.
Some major drawbacks limit the usefulness of studies regarding molecular pathogenesis in HCC. First, few studies analyze the genetic and genomic alterations that emerge at different time points during the entire progressive process of the disease. Second, the limited size of the studies is often a factor that undermines the capability to provide consistent genomic data. Animal models of hepatocarcinogenesis summarize the primal biology of liver tumorigenesis and have provided reliable data for understanding the cellular development of HCC in humans[1, 10, 11]. In the present study, the pathologic changes of livers in rats treated by DEN included non-specific injuries, regeneration and repair, fibrosis, and cirrhosis, dysplastic nodules, early tumorous nodules, advanced tumorous nodules and metastasis foci, resembling the process of human hepatocarcinogenesis. DEGs obtained by compare normal rats with DEN-treated animals at stages from cirrhosis to metastasis allowed us to screen for upregulated and downregulated gene expressional profiles. The number of DEGs at each stage was large and the information obtained was powerful. We were thus able to visualize the complicated process of hepatocarcinogenesis at the genomic level.
The annotated information of the DEGs show that extensive and diverse biological processes and molecular functions are involved in hepatocaricnogenesis. Most of the DEGs are involved in metabolism and transport, indicating that significant alterations occurred in the process of metabolism and transport during the developmnet of HCC. For example, tumor cells always perform anaerobic glycolysis, even when there is an adequate oxygen supply[12, 13], partly a result of alterations in the profile of enzymes associated with glycolysis. In this study, the gene expression level of lactate dehydrogenase B increased from the cirrhosis phase to the metastasis phase. Evidence shows that some genetic changes promoting tumor growth influences glucose energy metabolism directly[14, 15]. Many intermediate products from glycolysis are used to synthesize proteins, nucleic acids and lipids by tumor cells, providing the essential materials for the growth and hyperplasia of tumor cells. For aggressive tumors, increased glycolysis and metabolism alterations often occurred. The microenvironment acidosis provided by the conversion of pyruvic acid to lactic acid promotes invasion and metastasis of tumor cells [16–18]. Therapy targeting glycolysis could preferentially kill tumor cells with metabolic alterations and there is experimental evidence that some types of tumors can be treated by inhibiting the glycolytic and phosphopentose pathways. In this study, some DEGs associated with metabolisms of glucose were shown in Figure 6A.
Fat metabolism have significant changes in the process of tumorigenesis, e.g. a high fat diet was related to the development of many tumors . Enhanced fat synthesis in tumor cells could not only support the increased membrane synthesis and energy metabolism, but also higher level of fatty acid synthetase provides the base for interpretation the relation between the fat metabolism and the capacity of hyperplasia and metastasis of tumor cells. Stearoyl-CoA desaturase (SCD), which have four known isomers, takes part in regulating lipid synthesis. SCD2 plays key roles in the early development and survival of embryos in mice, whose expressional levels in the livers of wild mice embryos and newborn mice were higher than that of adult mice. Inhibition of lipid synthesis caused by the depletion of SCD2 was related to the decreased expression level of peroxisome proliferator-activated receptor gamma (PPAR-γ). Fatty acid binding proteins (FABPs) are proteins that could bind to fatty acid and other lipids reversibly. Researchers found expression of FABP5, coding epidermal fatty acid binding protein (E-FABP-GenBank Accession), upregulated in primary tongue carcinomas. FABP4, as a bridge between the inflammation and other metabolism syndromes, could not only transport the nuclear receptor PPAR-γ from cytoplasm to nucleus but also cause increased transcript activation of it. In this study, the expressional levels of SCD2, FABP4 and FABP5 increased during the process from cirrhosis to metastasis in rat model, suggesting that an alteration of the fat metabolism occurred in hepatocarcinogenesis of rat model. Other DEGs associated with fatty metabolisms were shown in Figure 6A.
In the present study, some enzymes related to the glutathione (GSH) metabolism were found to be significantly altered. For example, the expressional level of Gstm3 (glutathione S-transferase, mu type 3) decreased in all stages of hepatocarcinogenesis, while the expression levels of of enzymes increased, which including Glul (Glutamate-ammonia ligase), Gclc (Glutamate-cysteine ligase, catalytic subunit), GPX2 (Glutathione peroxidase 2), GPX3 (Glutathioneperoxidase 3), GSR (Glutathione reductase), Yc2 (Glutathione S-transferase Yc2 subunit), Gstm5 (Glutathione S-transferase, mu 5), Gstp1 (Glutathione-S-transferase, pi 1) and GSS (Glutathione synthetase). Some studies reported that GSH and the associated enzymes were considered to promot the tumor transformation from dysplastic nodules and take part in the development and progression of hepatocarcinomas[26, 27]. In gastroenteric tumors, the increase of GSH and the associated enzymes exhibit activation, indicating the existence of ROS and anti-oxidation defense related to GSH metabolism. There were five binding sites for β-catenin/TCF at the promoter region of GPX2, indicating that GPX2 might take part in the corresponding signal pathways. Thus previous research and our data indicate that genes related to oxidative stress and GSH metabolism play important roles in the process of progression from dysplastic nodules to tumor. The expressional level of GSH increased in tissue of HCC and the active hyperplasia liver cells[30, 31]. Research has shown that DNA oxidative injury is increased in human HCC [32, 33].
Many other enzymes associated with metabolism are involved in the defense and stress reaction, such as oxidative stress. For example, AKR1B7 (aldo-keto reductase family 1, member 7) takes part in the detoxification of oxides, such as aldehyde. During the detoxification of aldehyde, the expressional level of AKR1B7 mRNA increased. There are five binding sites with NF-κB at the 5' upstream region of the AKR1B7 gene, and oxidative stress upregulates the expression of AKR1B7 mediated by NF-κB[34, 35]. The expression level of aldehyde dehydrogenase ALDH3A1 (Aldehyde dehydrogenase 3A1) also increased after oxidative stress. In the present study, the expression levels of AKR1B7, AKR1B8 and ALDH3A1 were up-regulated at all stages of hepatocarcinogenesis. In the tumor cells, reactive oxigen species (ROS) was produced through the oxidative stress. ROS as signal molecules mediate various reactions relating to growth, such as angiogenesis. ROS in endothelial cells is mainly from NADPH oxidation enzymes, consisting of Nox1, Nox2, Nox4, Nox5, p22 (phox), p47 (phox) and Rac1 (small G-protein Rac1). NADPH oxidative enzymes were activated by different factors including VEGF, angiopoietin-1, hypoxia and ischemia. Furthermore, ROS has been shown to be involved in spontaneous phosphorylation. Nox4 mediated growth factors, such as anti-apoptosis of IGF, is partly due to the ROS produced by NADPH oxidative enzymes inhibiting the key protein tyrosine phosphatases(PTPs), then continually causing JAK2 kinase phosphorylation which resists the apoptosis reaction. In this study, However, the gene expression level of NADPH oxidative enzymes decreased in the livers of our rat model at all stages of hepatocarcinogenesis. The mechanism is unclear.
Cytochrome P450s (CYPs) are key enzymes in tumorigenesis, taking part in the activation and inactivation of chemotherapeutic agents in tumor tissues. The expression level of CYP1B1 in breast carcinomas was up-regulated significantly, providing a new therapy target and phenotype biomarker. The significant increase in CYP2E1 correlates with invasiveness and is a potential prognosis factor[39, 40]. Other studies have shown that the expression of CYP could influence the synthesis of arachidonic acid derivatives, thus altering the various downstream signal pathways, which was thought to be the prelude of carcinogenesis. At present, it is not possible to always predict the efficacy of a chemotherapeutic agent or its degree of toxicity. Analysis of CYPs expressional levels in tumor cells may allow prognosis decisions and therapy predictions. In this study, only the expression level of CYP2C40 increased at all stages of hepatocarcinogenesis in rat models, while the remaining CYPs decreased (Figure 6C). Clearly, further investigation is needed to determine the role(s) of CYPs in hepatocarcinogenesis.
In addition to the deregulated expression of metabolism associated genes, we noticed that among the DEGs in the hepatocarcinogenesis of rat models, some known tumor-associated genes, such as Rb1 and Myc, showed deregulated expression occurring at all the stages of hepatocarcinogenesis. Their persisting activation or deactivation could induce the tumor phenotype, as well as play roles at the later stage of progression and metastasis. Meanwhile, some known metastasis-associated genes are found deregulated at the promotion stage of tumor development. For example, the expression level of Ndrg2 and Hrasls3 (HRAS like suppressor 3) decreased at all stages compared to the normal livers, while the expression level of Nme1 (expressed in non-metastatic cells 1) increased. Generally, it was thought that genes involved in the development of carcinoma activation participated at the early stage, while genes participating in the metastasis were activated at the latter stage of tumor progression. In opposition to the traditional model, Bernards and Weinberg proposed that the metastatic ability of tumor cells occurred at the early stage of tumor development. Some oncogenes such as Ras and Src assigned the tumor cells with the metastatic phenotype [44–46].
As we known, the important characteristic of malignant tumor cells is the capability of invading the vicinity, forming metastasis foci at the remote organ, overcoming the host's control over the microenvironment[47, 48]. The malignant transformation of liver cells occurred on the basis of chronic injury, regeneration and cirrhosis. The liver cancer cells could synthesize ECM components and the ECM surrounding liver cancer cells was found to be different from that of stroma in the normal organ [49–51]. Integrin and laminin are the major components of ECM. The interaction between integrin and laminin is closely related to the signal transduction, providing survival signals for the cells, mediating the liver cancer cells formation of pseudopodia, and adherence with laminin, which are imperative if a liver cancer cell is to migrate and invade [52–55]. In the process of hepatocarcinogenesis in this rat model, the deregulated expression of many ECM associated genes plays important roles in the hepatocarcinogenesis, e.g. Itga6, Lamc1, Col1a1 and Spp1, etc. (Table 2, 3 and additonal file 2). The differential expression profile of ECM associated genes in time course and space is very important to cellular proliferation and migration.
It is worthy of note that the expression of the annexin family members and calcium-binding proteins are deregulated in this rat model of liver cancer. In our study, the expressional level of Annexin A1, A2, A3, A5 and A7 increased compared with the normal liver tissue. Annexins consist of a conserved protein family. Annexin A2 is closely associated with cell division regulation and tumor growth, and is deregulated in many tumors[56, 57]. Two Annexin A2 molecules bind to the long chains of p11/S100A10 dimers through its N-terminals, form the isotetramer, regulating the reactions of Annexin A2 and membranes and actin in cortical areas, and the distribution of recirculating endosomes. In addition, S100A10 and Annexin A2 form isodimers, prompting the invasion and metastasis of the tumor by activating plasminogen. In the present study, the expression level of S100a10, S100a11, S100a6, S100a8 and S100a9 increased from cirrhosis to metastatic process when compared with the normal liver. S100A8/A9 form the compounds that play a role in inducing apoptosis in tumor cells. S100A8/A9 at low concentrations prompts growth activity, the phosphorylation of MAPK pathway and NF-κB is activated in cells after S100A8/A9 treatment.
The majority of HCCs slowly unfold against a background of chronic hepatitis and cirrhosis, which can be considered as preneoplastic conditions of the liver. Chronic hepatitis is characterized by persistent inflammation, cytokine and oxidative stress-mediated hepatocyte death and active proliferation of residual hepatocytes to replace the lost parenchyma[1, 60]. During the process of hepatocarcinogenesis in rat models, chronic inflammation precedes cirrhosis. Epidemiology studies showed that chronic inflammation increased the risk of tumors, and the microenvironment of tumorigenesis resembles the reaction of inflammation to injury in many ways. In the tumor microenvironment, the chemotactic factors and receptors mediated angiogenesis, recruited cells, prompting cellular survival and proliferation. On the other hand, oxidative stress occurred in inflammatory processes. The inflammatory cells and tumor cells both produce free radicals and soluble factors such as arachidonic acid, cytokines and chemotactic factors, seubsequently producing reactive oxygen. All these factors strongly recruit the inflammatory cells to produce cytokines, which promotes a vicious cycle. The intermediate products of active oxygen oxidize DNA directly or interfere with DNA repair. These oxides activate protein, carbohydrate and lipids quickly, the derived products interfere with inter- and intracellular homeostasis, favoring DNA mutation. Thus, the chronic inflammation prompts the malignant transformation of cells. Chronic inflammation also favors angiogenesis. In the present study, many DEGs are related to inflammation reaction, immune reaction and stress. (Tables 4 and 5)
In conclusion, analysis of gene expression profile is a useful tool to provide new clues and produce new research targets in the hepatocarcinogenesis field. In the present study, a rat model of liver cancer was established. We have listed the deregulated expression genes in the process from cirrhosis to liver cancer in the DEN-treated rat model. As indicated in the literature, this model shows that cirrhosis is a precancerous lesion of the liver. Although we only discuss some parts of the great quantity of information in this study, much unknown information remains. Functional analysis of these genes revealed discrete expression clusters, including cell proliferation, protein synthesis, and hepatocyte-specific functions. We still need to discern the key genes playing core roles in the promotion and progression of liver cancer. In this article, we focused our attention on the global molecular events that occurred in DEN-treated rats (and probably represent the earliest ones that start the multistep process of hepatocarcinogenesis). Additional information may be mined from this and similar studies to provide clues to many areas including the very important search for diagnostic markers, therapy targets and prognosis prediction markers.
differential expression genes.
Thorgeirsson SS, Grisham JW: Molecular pathogenesis of human hepatocellular carcinoma. Nat Genet. 2002, 31: 339-346. 10.1038/ng0802-339.
Nagai H, Pineau P, Tiollais P, Buendia MA, Dejean A: Comprehensive allelotyping of human hepatocellular carcinoma. Oncogene. 1997, 14: 2927-2933. 10.1038/sj.onc.1201136.
Lee JS, Grisham JW, Thorgeirsson SS: Comparative functional genomics for identifying models of human cancer. Carcinogenesis. 2005, 26: 1013-1020. 10.1093/carcin/bgi030.
Zender L, Xue W, Zuber J, Semighini CP, Krasnitz A, Ma B, Zender P, Kubicka S, Luk JM, Schirmacher P, McCombie WR, Wigler M, Hicks J, Hannon GJ, Powers S, Lowe SW: An oncogenomics-based in vivo RNAi screen identifies tumor suppressors in liver cancer. Cell. 2008, 13: 852-864. 10.1016/j.cell.2008.09.061.
Yachida S, Sakamoto M, Imaida K, Yokohira M, Saoo K, Okano K, Wakabayashi H, Maeta H, Suzuki Y: p27 is overexpressed in very early stages of hepatocarcinogenesis. Cancer Sci. 2008, 99: 2152-2159. 10.1111/j.1349-7006.2008.00923.x.
Jiang Sheng Ze, Fang Gang Shi, Gao Yi, Wang Shuang, Chen Feng Jian: Dynamic changes of metrix metalloproteinases in liver tissue during the development of diethylinitrosamine-induced rat heptocarcinoma. World Chin J Digestol. 2001, 9 (7): 759-762.
Umeda T, Hino O: Molecular aspects of human hepatocarcinogenesis mediated by inflammation: From hypercarcinogenic state to normo- or hypo carcinogenic state. Oncology. 2002, 62 (Suppl): 138-142. 10.1159/000048274.
Mantovani A: Cancer: Inflammation by remote control. Nature. 2005, 435: 752-753. 10.1038/435752a.
Farazi PA, DePinho RA: Hepatocellular carcinoma pathogenesis: from genes to environment. Nat Rev Cancer. 2006, 6: 674-687. 10.1038/nrc1934.
Lee JS, Chu IS, Heo J, Calvisi DF, Sun Z, Roskams T, Durnez A, Demetris AJ, Thorgeirsson SS: Classification and prediction of survival in hepatocellular carcinoma by gene expression profiling. Hepatology. 2004, 40: 667-676. 10.1002/hep.20375.
Feitelson MA, Pan J, Lian Z: Early molecular and genetic determinants of primary liver malignancy. Surg Clin North Am. 2004, 84: 339-354. 10.1016/S0039-6109(03)00226-3.
Garber K: Energy boost. the Warburg Effect returns in a new theory of cancer. J Natl Cancer Inst. 2004, 96: 1805-1806.
Chen H, Yue JX, Yang SH, Ding H, Zhao RW, Zhang S: Overexpression of transketolase-like gene 1 is associated with cell proliferation in uterine cervix cancer. J Exp Clin Cancer Res. 2009, 28: 43-10.1186/1756-9966-28-43.
Pelicano H, Martin DS, Xu RH, Huang P: Glycolysis inhibition for anticancer treatment. Oncogene. 2006, 25: 4633-4646. 10.1038/sj.onc.1209597.
Wittig R, Coy JF: The Role of Glucose Metabolism and Glucose-Associated Signalling in Cancer. Perspectives in Medicinal Chemistry. 2007, 1: 64-82.
Gatenby RA, Gillies RJ: Why do cancers have high aerobic glycolysis?. Nat Rev Cancer. 2004, 4: 891-899. 10.1038/nrc1478.
Stern R, Shuster S, Neudecker BA, Formby B: Lactate stimulates fibroblast expression of hyaluronan and CD44: the Warburg effect revisited. Exp Cell Res. 2002, 276: 24-31. 10.1006/excr.2002.5508.
Walenta S, Wetterling M, Lehrke M, Schwickert G, Sundfør K, Rofstad EK, Mueller-Klieser W: High lactate levels predict likelihood of metastases, tumor recurrence, and restricted patient survival in human cervical cancers. Cancer Res. 2000, 60: 916-921.
Philips BJ, Dhir R, Hutzley J, Sen M, Kelavkar UP: Polyunsaturated fatty acid metabolizing 15-Lipoxygenase-1 (15-LO-1) expression in normal and tumorigenic human bladder tissues. Appl Immunohistochem Mol Morphol. 2008, 16: 159-164. 10.1097/PAI.0b013e31805baa41.
Young CD, Anderson SM: Sugar and fat – that's where it's at: metabolic changes in tumors. Breast Cancer Res. 2008, 10: 202-10.1186/bcr1852.
Miyazaki M, Dobrzyn A, Elias PM, Ntambi JM: Stearoyl-CoA desaturase-2 gene expression is required for lipid synthesis during early skin and liver development. Proc Natl Acad Sci USA. 2005, 102: 12501-12506. 10.1073/pnas.0503132102.
Christianson JL, Nicoloro S, Straubhaar J, Czech MP: Stearoyl-CoA desaturase 2 is required for peroxisome proliferator-activated receptor gamma expression and adipogenesis in cultured 3T3-L1 cells. J Biol Chem. 2008, 283: 2906-1916. 10.1074/jbc.M705656200.
Uma RS, Naresh KN, D'Cruz AK, Mulherkar R, Borges AM: Metastasis of squamous cell carcinoma of the oral tongue is associated with down-regulation of epidermal fatty acid binding protein (E-FABP). Oral Oncol. 2007, 43: 27-32. 10.1016/j.oraloncology.2005.12.024.
Ordovas JM: Identification of a functional polymorphism at the adipose fatty acid binding protein gene (FABP4) and demonstration of its association with cardiovascular disease: a path to follow. Nutr Rev. 2007, 65: 130-134. 10.1111/j.1753-4887.2007.tb00290.x.
Gillilan RE, Ayers SD, Noy N: Structural basis for activation of fatty acid-binding protein 4. J Mol Biol. 2007, 372: 1246-1260. 10.1016/j.jmb.2007.07.040.
Hendrich S, Campbell HA, Pitot HC: Quantitative stereological evaluation of four histochemical markers of altered foci in multistage hepatocarcinogenesis in the rat. Carcinogenesis. 1987, 8: 1245-1250. 10.1093/carcin/8.9.1245.
Higashi K, Hiai H, Higashi T, Muramatsu M: Regulatory mechanism of glutathione S-transferase P-form during chemical hepatocarcinogenesis: old wine in a new bottle. Cancer Lett. 2004, 209: 155-163. 10.1016/j.canlet.2004.01.003.
Scibior D, Skrzycki M, Podsiad M, Czeczot H: Glutathione level and glutathione-dependent enzyme activities in blood serum of patients with gastrointestinal tract tumors. Clin Biochem. 2008, 41: 852-858. 10.1016/j.clinbiochem.2008.03.005.
Kipp A, Banning A, Brigelius-Flohé R: Activation of the glutathione peroxidase 2 (GPx2) promoter by beta-catenin. Biol Chem. 2007, 388: 1027-1033. 10.1515/BC.2007.137.
Lu SC: Regulation of hepatic glutathione synthesis: current concepts and controversies. FASEB J. 1999, 13: 1169-1183.
Huang ZZ, Chen C, Zeng Z, Yang H, Oh J, Chen L, Lu SC: Mechanism and significance of increased glutathione level in human hepatocellular carcinoma and liver regeneration. FASEB J. 2001, 15: 19-21.
Schwarz KB, Kew M, Klein A, Abrams RA, Sitzmann J, Jones L, Sharma S, Britton RS, Di Bisceglie AM, Groopman J: Increased hepatic oxidative DNA damage in patients with hepatocellular carcinoma. Dig Dis Sci. 2001, 46: 2173-2178. 10.1023/A:1011958814371.
Jungst C, Cheng B, Gehrke R, Schmitz V, Nischalke HD, Ramakers J, Schramel P, Schirmacher P, Sauerbruch T, Caselmann WH: Oxidative damage is increased in human liver tissue adjacent to hepatocellular carcinoma. Hepatology. 2004, 39: 1663-1672. 10.1002/hep.20241.
Baumann C, Davies B, Peters M, Kaufmann-Reiche U, Lessl M, Theuring F: AKR1B7 (mouse vas deferens protein) is dispensable for mouse development and reproductive success. Reproduction. 2007, 134: 97-109. 10.1530/REP-07-0022.
Jia G, Takahashi R, Zhang Z, Tsuji Y, Sone H: Aldo-keto reductase 1 family B7 is the gene induced in response to oxidative stress in the livers of Long-Evans Cinnamon rats. Int J Oncol. 2006, 29: 829-838.
Ushio-Fukai M, Nakamura Y: Reactive oxygen species and angiogenesis: NADPH oxidase as target for cancer therapy. Cancer Lett. 2008, 266: 37-52. 10.1016/j.canlet.2008.02.044.
Lee JK, Edderkaoui M, Truong P, Ohno I, Jang KT, Berti A, Pandol SJ, Gukovskaya AS: NADPH oxidase promotes pancreatic cancer cell survival via inhibiting JAK2 dephosphorylation by tyrosine phosphatases. Gastroenterology. 2007, 133: 1637-1648. 10.1053/j.gastro.2007.08.022.
Rodriguez-Antona C, Ingelman-Sundberg M: Cytochrome P450 pharmacogenetics and cancer. Oncogene. 2006, 25: 1679-1691. 10.1038/sj.onc.1209377.
McFadyen MC, Murray GI: Cytochrome P450 1B1: a novel anticancer therapeutic target. Future Oncol. 2005, 1: 259-263. 10.1517/147966184.108.40.2069.
Vaclavikova R, Hubackova M, Stribrna-Sarmanova J, Kodet R, Mrhalova M, Novotny J, Gut I, Soucek P: RNA expression of cytochrome P450 in breast cancer patients. Anticancer Res. 2007, 27: 4443-4450.
Nebert DW, Dalton TP: The role of cytochrome P450 enzymes in endogenous signalling pathways and environmental carcinogenesis. Nat Rev Cancer. 2006, 6: 947-960. 10.1038/nrc2015.
Poste G, Fidler IJ: The pathogenesis of cancer metastasis. Nature. 1980, 283: 139-146. 10.1038/283139a0.
Bernards R, Weinberg RA: A progression puzzle. Nature. 2002, 418: 823-10.1038/418823a.
Chambers AF, Tuck AB: Ras-responsive genes and tumor metastasis. Crit Rev Oncog. 1993, 4: 95-114.
Yu D, Wang SS, Dulski KM, Tsai CM, Nicolson GL, Hung MC: c-erbB-2/neu overexpression enhances metastatic potential of human lung cancer cells by induction of metastasis-associated properties. Cancer Res. 1994, 54: 3260-3266.
Myoui A, Nishimura R, Williams PJ, Hiraga T, Tamura D, Michigami T, Mundy GR, Yoneda T: C-SRC tyrosine kinase activity is associated with tumor colonization in bone and lung in an animal model of human breast cancer metastasis. Cancer Res. 2003, 63: 5028-5033.
Cairns RA, Khokha R, Hill RP: Molecular mechanisms of tumor invasion and metastasis: an integrated view. Curr Mol Med. 2003, 3: 659-671. 10.2174/1566524033479447.
Coussens LM, Werb Z: Inflammation and cancer. Nature. 2002, 420: 860-867. 10.1038/nature01322.
Grigioni WF, Garbisa S, D'Errico A, Baccarini P, Stetler-Stevenson WG, Liotta LA, Mancini AM: Evaluation of hepatocellular carcinoma aggressiveness by a panel of extracellular matrix antigens. Am J Pathol. 1991, 138: 647-654.
Torimura T, Ueno T, Inuzuka S, Kin M, Ohira H, Kimura Y, Majima Y, Sata M, Abe H, Tanikawa K: The extracellular matrix in hepatocellular carcinoma shows different localization patterns depending on the differentiation and the histological pattern of tumors: immunohistochemical analysis. J Hepatol. 1994, 21: 37-46. 10.1016/S0168-8278(94)80134-7.
Bissell DM: Chronic liver injury, TGF-beta, and cancer. Exp Mol Med. 2001, 33: 179-190.
Carloni V, Romanelli RG, Mercurio AM, Pinzani M, Laffi G, Cotrozzi G, Gentilini P: Knockout of alpha6beta1-integrin expression reverses the transformed phenotype of hepatocarcinoma cells. Gastroenterology. 1998, 115: 433-442. 10.1016/S0016-5085(98)70210-0.
Yang C, Zeisberg M, Lively JC, Nyberg P, Afdhal N, Kalluri R: Integrin alpha1beta1 and alpha2beta1 are the key regulators of hepatocarcinoma cell invasion across the fibrotic matrix microenvironment. Cancer Res. 2003, 63: 8312-8317.
Giannelli G, Bergamini C, Fransvea E, Marinosci F, Quaranta V, Antonaci S: Human hepatocellular carcinoma (HCC) cells require both alpha3beta1 integrin and matrix metalloproteinases activity for migration and invasion. Lab Invest. 2001, 81: 613-627.
Fu BH, Wu ZZ, Dong C: Integrin beta1 mediates hepatocellular carcinoma cells chemotaxis to laminin. Hepatobiliary Pancreat Dis Int. 2004, 3: 548-551.
Brichory FM, Misek DE, Yim AM, Krause MC, Giordano TJ, Beer DG, Hanash SM: An immune response manifested by the common occurrence of Annexin I and Annexin II autoantibodies and high circulating levels of IL-6 in lung cancer. Proc Natl Acad Sci USA. 2001, 98: 9824-9829. 10.1073/pnas.171320598.
Emoto K, Yamada Y, Sawada H, Fujimoto H, Ueno M, Takayama T, Kamada K, Naito A, Hirao S, Nakajima Y: Annexin II overexpression correlates with stromal tenascin-C overexpression: a prognostic marker in colorectal carcinoma. Cancer. 2001, 92: 1419-1426. 10.1002/1097-0142(20010915)92:6<1419::AID-CNCR1465>3.0.CO;2-J.
Morel E, Gruenberg J: The p11/S100A10 light chain of annexin A2 is dispensable for annexin A2 association to endosomes and functions in endosomal transport. PLoS ONE. 2007, 2: e1118-10.1371/journal.pone.0001118.
Ito Y, Arai K, Nozawa R, Yoshida H, Higashiyama T, Takamura Y, Miya A, Kobayashi K, Kuma K, Miyauchi A: S100A10 expression in thyroid neoplasms originating from the follicular epithelium: contribution to the aggressive characteristic of anaplastic carcinoma. Anticancer Res. 2007, 27: 2679-2783.
Coleman WB: Mechanisms of human hepatocarcinogenesis. Curr Mol Med. 2003, 3: 573-588. 10.2174/1566524033479546.
Coussens LM, Werb Z: Inflammation and cancer. Nature. 2002, 420: 860-867. 10.1038/nature01322.
Slaga TJ, Lichti U, Hennings H, Elgjo K, Yuspa SH: Effects of tumor promoters and steroidal anti-inflammatory agents on skin of newborn mice in vivo and in vitro. J Natl Cancer Inst. 1978, 60: 425-431.
Jackson JR, Seed MP, Kircher CH, Willoughby DA, Winkler JD: The codependence of angiogenesis and chronic inflammation. FASEB J. 1997, 11: 457-465.
This study was supported by a grant from the Science Foundation of the Ministry of Health of China (No. wkj2004 -2-12). The authors thanks for technician Yuhua Li for assistance with preparation of tissue slices. We would like to thank Shanghai Biochip Co., Ltd. And the National Engineering Center for Biochip at Shanghai for the operation of the Affymetrix genechips.
The authors declare that they have no competing interests.
YFL wrote the manuscript. BSZ performed the validation of genes. HLZ and XJZ established the animal model. YHL prepared the tissue slides. JZ helped write the manuscript. JPZ, ZQF and XHG participated in the design of the study and helped to draft the manuscript. All authors read and approved the final manuscript.
and Jian-Ping Zhang contributed equally to this work.
Electronic supplementary material
Additional file 2: The up-regulated DEGs sharing from cirrhosis to metastasis sorted out by the following GO function. for the screened DEGs sharing from stage of liver cirrhosis to metastasis sorted out by the GO words: angiogenesis, apoptosis, cell adhesion, cell migration, cell proliferation and extracellular matrix. (PDF 46 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.