Open Access
Open Peer Review

This article has Open Peer Review reports available.

How does Open Peer Review work?

Gene expression profile analysis of human hepatocellular carcinoma using SAGE and LongSAGE

  • Hui Dong1, 2, 5,
  • Xijin Ge3,
  • Yan Shen4,
  • Linlei Chen6,
  • Yalin Kong7,
  • Hongyi Zhang7,
  • Xiaobo Man2,
  • Liang Tang2,
  • Hong Yuan6,
  • Hongyang Wang2,
  • Guoping Zhao1, 4, 5Email author and
  • Weirong Jin4, 5Email author
Contributed equally
BMC Medical Genomics20092:5

DOI: 10.1186/1755-8794-2-5

Received: 09 October 2008

Accepted: 26 January 2009

Published: 26 January 2009

Abstract

Background

Hepatocellular carcinoma (HCC) is one of the most common cancers worldwide and the second cancer killer in China. The initiation and malignant transformation of cancer result from accumulation of genetic changes in the sequences or expression level of cancer-related genes. It is of particular importance to determine gene expression profiles of cancers on a global scale. SAGE and LongSAGE have been developed for this purpose.

Methods

We performed SAGE in normal liver and HCC samples as well as the liver cancer cell line HepG2. Meanwhile, the same HCC sample was simultaneously analyzed using LongSAGE. Computational analysis was carried out to identify differentially expressed genes between normal liver and HCC which were further validated by real-time quantitative RT-PCR.

Results

Approximately 50,000 tags were sequenced for each of the four libraries. Analysis of the technical replicates of HCC indicated that excluding the low abundance tags, the reproducibility of SAGE data is high (R = 0.97). Compared with the gene expression profile of normal liver, 224 genes related to biosynthesis, cell proliferation, signal transduction, cellular metabolism and transport were identified to be differentially expressed in HCC. Overexpression of some transcripts selected from SAGE data was validated by real-time quantitative RT-PCR. Interestingly, sarcoglycan-ε (SGCE) and paternally expressed gene (PEG10) which is a pair of close neighboring genes on chromosome 7q21, showed similar enhanced expression patterns in HCC, implicating that a common mechanism of deregulation may be shared by these two genes.

Conclusion

Our study depicted the expression profile of HCC on a genome-wide scale without the restriction of annotation databases, and provided novel candidate genes that might be related to HCC.

Background

Being a worldwide malignant liver tumor, hepatocellular carcinoma (HCC) ranks the fifth in frequency among common human solid tumors and the fourth leading cause of cancer-related death [1, 2]. The majority of HCC cases occur in Asia and Sub-Saharan Africa, but incidence has been increasing in Western Europe and the United States in recent years [3, 4]. In China, HCC is now the second cancer killer [5]. Numerous studies have been carried out in an effort to elucidate molecular mechanism of hepatocarcinogenesis, metastasis and/or prognosis. Multiple genes have been reported to be involved in the development of HCC. For instance, mutations of p53, β-catenin and AXIN1 [610], activation of oncogenes such as c-myc, c-met, c-jun, N-ras and nuclear factor κB [1116], and up-regulation of a set of genes including GPC3, LCN2, and IGFBP-1 were observed in a certain number of HCC cases [1719]. These studies contributed greatly to our understanding of HCC in terms of individual molecules.

It is well known that the initiation and malignant transformation of cancer result from accumulation of genetic changes in the sequences or expression level of cancer-related genes. Thus it is of particular importance to determine gene expression profiles of cancers on a global scale. Several techniques have been developed for this purpose. Serial Analysis of Gene Expression (SAGE) allows quantitative measurement of gene expression profile through sequencing of short tags. Compared with microarrays, the expression profiles obtained by SAGE is exploratory in nature as it is not restricted to available annotations. It has been widely used in cancer studies since it was first introduced by Velculescu et al. in 1995 [20]. Saha et al. further developed LongSAGE which detects 21-bp tags instead of the original 14-bp tags, thus making it feasible to directly map the LongSAGE tags to genomic sequence data [21].

In the present study, we employed SAGE to comprehensively analyze gene expression profiles of normal human liver, HCC and the liver cancer cell line HepG2. The same HCC sample was further analyzed simultaneously using LongSAGE. To the best of our knowledge, there is no SAGE or LongSAGE expression profile of HCC available in the public domain. Our study provided a wealth of information on HCC expression profile, and comparative analysis of datasets of HCC vs normal liver led to the identification of a series of differentially expressed genes. The following real-time quantitative RT-PCR analysis further identified novel candidate genes potentially involved in hepatocarcinogenesis of HCC.

Methods

Cell culture and tissue samples

Human HCC cell line HepG2 was cultured in Dulbecco's modified Eagle's medium (Invitrogen, Carlsbad, CA) with 10% fetal bovine serum under 5% CO2 in a humidified incubator at 37°C. Cells were harvested at 80–90% confluence. Normal liver tissue was obtained from a patient who underwent hepatectomy because of hepatic hemangioma. Clinical HCC samples and corresponding non-tumorous liver tissues were derived from HCC patients and kept frozen in liquid nitrogen immediately after separation. All samples were collected with informed consent in accordance with the standards of Institutional Human Subjects Protection Review Board.

SAGE and LongSAGE

Total RNAs were extracted from HepG2 and liver tissues using TRIzol reagent (Invitrogen, Carlsbad, CA) and then treated with DNaseI (Rhoch Diagnostics, Almere, The Netherlands) according to the manufacture's protocol. Polyadenylated mRNAs were then isolated using Oligotex Direct mRNA Midi Kit (QIAGEN). Three SAGE libraries were constructed starting with 200 ng polyA+ mRNA from HepG2, normal liver tissue and HCC tissue respectively according to SAGE protocol [20, 22]. LongSAGE was performed with the same amount of polyA+ mRNA from the same HCC tissue used in SAGE following the procedures described previously [21]. Sequencing of SAGE and LongSAGE tags were carried out using BigDye Terminator Cycle Sequencing Kits and ABI 3700 DNA Sequencers (Applied Biosystems, Foster City, CA). The sequence and frequency of 14-bp tags from SAGE and 21-bp tags from LongSAGE were extracted from the raw sequence files of concatenated di-tags using SAGE2000 analysis software version 4.5 (kindly provided by Dr. K. Kinzler, Johns Hopkins University School of Medicine).

Computational Analysis

To exclude sequencing errors, only tags detected at least twice in the four SAGE libraries (three SAGE libraries and one LongSAGE library) were included for further analysis. SAGEmap reliable mapping (http://www.ncbi.nlm.nih.gov/SAGE, UniGene Build #182) [23] was used to establish tag-to-gene assignments. Comparison analysis between SAGE libraries was performed using the IDEG6 software http://telethon.bio.unipd.it/bioinfo/IDEG6/[24]. Differentially expressed tags were selected according to the cut-off threshold (p value < 0.05, fold change > 3) calculated according to pairwise comparison algorithm [25]. The Cluster and Treeview programs were used to generate the average linkage hierarchical clustering and visualize changes of gene expression [26]. Functional classification of genes was performed using DAVID program http://david.abcc.ncifcrf.gov/ obtained from NIAID (National Institute of Allergy and Infectious Disease).

Real-time quantitative RT-PCR

Total RNAs were extracted from HCC samples and corresponding nontumorous liver tissues using TRIzol reagent (Invitrogen) followed by treating with DNaseI according to the manufacture's protocol. 2 μg of total RNA was reverse transcribed in a total volume of 20 μl containing 200 units of SuperScript II RNase H- Reverse Transcriptase (Invitrogen), 50 mM Tris-HCl (pH8.3), 75 mM KCl, 3 mM MgCl2, 10 mM dithiothreitol, 500 μM dNTPs each, 500 ng oligo(dT)23 primer and 40 units of RNaseOUT at 42°C for 50 min, followed by inactivating at 70°C for 15 min.

For real-time quantitative RT-PCR, 1 μl of the first strand cDNA and 5 μl of 2× SYBR Green PCR mix (Applied Biosystems, Foster City, CA) was added to a final volume of 10 μl. The thermal cycles were performed on LightCycler (Roche) following conditions at 95°C for 30 s, 40 cycles at 95°C for 5 s, 60°C for 5 s and 72°C for 30 s. Each reaction was performed in triplicate. The threshold cycle (Ct) value was used to calculate the expression ratio of target genes to internal control gene (GAPDH) with the formula 2Ct(GAPDH)-Ct(target gene). Two-side student's t test was applied in the statistical analysis of quantitative RT-PCR data.

All the sequences of differently expressed transcripts were obtained from GenBank http://www.ncbi.nlm.nih.gov/Genbank/index.html. PCR Primers were designed using Primer 3.0 http://frodo.wi.mit.edu/. One set of primers were also designed for GAPDH which is considered an internal control for RT-PCR. The sequences of primers for each transcript were listed in Table 1.
Table 1

Sequences of primers used in RT-PCR

Unigene

GenBank

Accession No.

Product

Size

Primers

GAPDH

NM_002046

106 bp

F:ATGGGTGTGAACCATGAGAAG

R:AGTTGTCATGGATGACCTTGG

CCL20

BC020698

200 bp

F:GCGCAAATCCAAAACAGACT

R:CAAGTCCAGTGAGGCACAAA

S100P

BC006819

116 bp

F:TACCAGGCTTCCTGCAGAGT

R:AGCCACGAACACGATGAACT

PEG10

NM_015068

111 bp

F:TCCTGTCTTCGCAGAGGAGT

R:TTCACTTCTGTGGGGATGGA

SGCE

NM_003919

220 bp

F:ATGCAAACACCAGACATCCA

R:TCTGATGTGGCAAGTTCTGC

XAGE-1

AF251237

159 bp

F:TCTGCAAGAGCTGCATCAGT

R:AGCTTGCGTTGTTTCAGCTT

COL4A1

NM_001845

228 bp

F:ACGGGGGAAAACATAAGACC

R:TGGCGCACTTCTAAACTCCT

ZFYVE1

NM_178441

213 bp

F:AACTGCTACGAAGCCAGGAA

R:GTTGTGGCAGTGGAGGATTT

ZNF83

BC050407

224 bp

F:TGGGAAGGTCTTCGGTCTAA

R:CGGCATGAAATCTCTGATGA

TNPO2

NM_013433

213 bp

F:GGCGTTGTGCAGGACTTTAT

R:GGCAGTCTCCATGATCACCT

TM4SF1

NM_014220

169 bp

F:AGGGCCAGTACCTTCTGGAT

R:AAGCCACATATGCCTCCAAG

Results

Summary of SAGE data

A total of 56,984, 63,800 and 47,149 tags were detected from the normal liver, HCC and HepG2 SAGE libraries respectively, and 51,632 tags were obtained from the HCC LongSAGE library. Approximately 31%~34% of the total tags of each SAGE library were unique, while unique tags accounted for 45% of total tags of the HCC LongSAGE library (Table 2). When using SAGEmap to annotate these tags, we found that about 50% of unique tags in a SAGE library matched to specific genes (single match), 25~28% to multiple genes (multiple match) and 22~25% were novel tags that do not match to any known genes. The pattern was quite different in the LongSAGE library, with 43% of unique tags being single matched tags, 3.5% multiple matched tags and 53% novel tags. The result indicated that by extending the tag length from 14 to 21 bp, the specificity of SAGE tags was substantially improved. More than 70% of total tags in each library were single copy tags, while tags with high copy numbers (> 5 copies) account for less than 10% of the total tags (Table 3). This result is consistent with the fact that the majority of human genes are expressed at low levels [27].
Table 2

Tag counts of SAGE and LongSAGE libraries

 

Total

Unique

(%)

Single matched

(%)

Multiple matched

(%)

Novel

(%)

Normal Liver

56,984

17,719

(31.1)

8,823

(49.8)

4,452

(25.1)

4,444

(25.1)

HepG2

47,149

15,827

(33.6)

7,956

(50.2)

4,427

(28.0)

3,444

(21.8)

HCC

63,800

20,313

(31.8)

10,048

(49.4)

5,375

(26.5)

4,890

(24.1)

HCC LongSAGE

51,632

23,093

(44.7)

9,973

(43.2)

801

(3.5)

12,319

(53.3)

Table 3

Frequency distributions of unique tags from the SAGE and LongSAGE libraries

Tag copy No.

Normal Liver

(%)

HepG2

(%)

HCC

(%)

HCC LongSAGE

(%)

1

12,842

(72.5)

11,223

(71.0)

14,477

(71.3)

18,171

(78.7)

2

2,123

2,015

2,495

2,205

3

872

847

1,002

947

4

447

429

558

469

5–10

864

837

1,083

823

11–20

291

248

364

256

21–99

231

166

268

180

100 and >

49

62

66

42

Total

17,719

15,827

20,313

23,093

Correlation analysis of SAGE data

As one set of SAGE data of normal liver is available in the public Gene Expression Omnibus database (http://www.ncbi.nlm.nih.gov/geo/, accession number GSM785), we compared our data of normal liver with it to evaluate the correlation between these two independent libraries. The result indicated that the data from public database correlated well with our data, with Pearson's correlation coefficient 0.74 (Figure 1A).
https://static-content.springer.com/image/art%3A10.1186%2F1755-8794-2-5/MediaObjects/12920_2008_Article_69_Fig1_HTML.jpg
Figure 1

Correlation analysis of SAGE and LongSAGE libraries. (A) Correlation between normal liver SAGE data from public database ftp://ftp.ncbi.nlm.nih.gov/pub/sage/obsolete/seq/ and our own normal liver data. Pearson's correlation coefficient 0.74. (B) Correlation between SAGE and LongSAGE data of HCC. Pearson's correlation coefficient 0.97. (C) Correlation between HCC and HepG2 SAGE data. Pearson's correlation coefficient 0.55. (D) Correlation between HCC and our own normal liver SAGE data. Pearson's correlation coefficient 0.70.

We also compared the SAGE and LongSAGE data which were obtained from the same HCC tissue sample of the same patient. We extracted the 14-bp SAGE tags from the 21-bp LongSAGE tags. The converted SAGE library could be considered as a technical replicate. Totally 6,815 overlapping tags were detected between these two libraries. Thus among the 20,313 unique tags in the original SAGE library, only 33.5% of them are detected by the second library derived from LongSAGE tags. Most of the missed tags are low abundance tags for which detection by random sampling is a stochastic process. Another reason might be possible sequencing errors that resulted in many singleton tags. However, analysis of these overlapping tags showed remarkable consistency between these two libraries with Pearson's correlation coefficient 0.97 (Figure 1B). This result reassured us the reproducibility of SAGE method on measuring gene expression.

We further compared the data of HCC with that of HepG2 to find out how much the gene expression profile has changed in cell line comparing with the HCC tissue. The difference between HCC and HepG2 was very notable, with Pearson's correlation coefficient 0.55, even greater than difference between HCC and normal liver (Pearson's correlation coefficient 0.70) (Figure 1C and 1D). This result suggested that the gene expression profile of cell line has changed too much to be regarded as the reference of its corresponding tissue.

Identification and functional analysis of genes differentially expressed between normal liver and HCC

In order to identify genes differentially expressed between normal liver and HCC, we used the IDEG6 software to compare the normal liver library with both the HCC library and the second HCC library derived from LongSAGE. Only genes consistently up-regulated or down-regulated in both comparisons were considered to be differentially expressed. Totally 224 genes were identified to be differentially expressed, with the significance p < 0.05 and fold change > 3 (Additional file 1). Figure 2 shows the expression patterns of these altered genes in normal liver, HCC and HepG2. The top 20 up-regulated and down-regulated genes ranked by fold change were listed in Table 4 and Table 5 respectively. Consistent with previous results [18, 28], our data also indicated that GPC3 gene was significantly up-regulated in HCC with fold change of 69. Among the 224 differentially expressed genes, the expression level of 104 genes was significantly higher in HCC than in normal liver (p < 0.05), while 120 genes was significantly lower in HCC (p < 0.05). In addition to these 224 annotated tags, our list of differentially expressed tags also includes 99 tags that could be matched to multiple genes, and 109 tags that were differentially expressed but could not be mapped to any known gene (data not shown). Some of these novel tags are extremely highly expressed in HCC. For example, the tag TAAGTTTGGG is detected 24 and 17 times in two HCC libraries but is not detected at all in the other two libraries. It will be of interest to further study these tags. In the current study, however, we will focus on the annotated tags.
https://static-content.springer.com/image/art%3A10.1186%2F1755-8794-2-5/MediaObjects/12920_2008_Article_69_Fig2_HTML.jpg
Figure 2

Hierarchical clustering for genes differentially expressed in HCC. For visual comparison, genes differentially expressed in HCC and normal liver were clustered by Treeview program. The expression pattern of these genes in HepG2 was also shown. The red color represents genes up-regulated, and the green color represents genes down-regulated. 1: normal liver, 2: normal liver from public database, 3: HCC SAGE, 4: HCC LongSAGE, 5: HepG2.

Table 4

Top 20 genes up-regulated in HCC

Tag

UniGene

Symbol

Name

Function

Fold Change

ACCCGCCGGG

Hs.335106

ZFYVE1

Zinc finger, FYVE domain containing 1

Transport

81

GATTTCTTTG

Hs.435036

GPC3

Glypican 3

extracellular matrix

69

GCCAAGAATC

Hs.307720

DTNB

Dystrobrevin, beta

Calcium and Zinc ion binding

56

CTAACTAGTT

Hs.473583

NSEP1

Nuclease sensitive element binding protein 1

Response to external stimulus

50

GCTTGAATAA

Hs.116724

AKR1B10

Aldo-keto reductase family 1, member B10 (aldose reductase)

Steroid metabolism

47

GTGACCACGG

Hs.436980

GRIN2C

Glutamate receptor, ionotropic, N-methyl D-aspartate 2C

Ion transport

39

CTCATCCTAC

Hs.416049

TNPO2

Transportin 2 (importin 3, karyopherin beta 2b)

Protein transport

36

GATGTATGTG

Hs.461085

 

EST, strongly similar to NP_775842.1

Unclassified

36

ACCAGCACTC

Hs.436911

AMBP

Alpha-1-microglobulin/bikunin precursor

Signal Transduction

32

CCGACGGGCG

Hs.198161

PLA2G4B

Phospholipase A2, group IVB (cytosolic)

Signal Transduction

32

GGTCAGTCGG

Hs.409352

FLJ20701

Hypothetical protein FLJ20701

Unclassified

29

ATGTAAAAAA

Hs.429608

C5orf18

Chromosome 5 open reading frame 18

Unclassified

28

GACCCACCAT

Hs.436911

AMBP

Alpha-1-microglobulin/bikunin precursor

Signal Transduction

28

CAAGTTTGCT

Hs.439552

EEF1A1

Eukaryotic translation elongation factor 1 alpha 1

Signal Transduction

26

TTGTAAACTT

Hs.509226

FKBP3

FK506 binding protein 3, 25kDa

Cellular metabolism

25

GAAGGTGATC

Hs.112208

XAGE-1

X antigen family, member 1

Cellular metabolism

25

GGGACGAGTG

Hs.351316

TM4SF1

Transmembrane 4 L six family member 1

Unclassified

21

GAGGGTGGCG

Hs.473847

SH3BGR

SH3 domain binding glutamic acid-rich protein

Cellular metabolism

18

TAAATACAGT

Hs.221497

PRO0149

PRO0149 protein

Unclassified

18

GTGGATGGAC

Hs.6418

GPR175

G protein-coupled receptor 175

Organogenesis

17

Table 5

Top 20 genes down-regulated in HCC

Tag

UniGene

Symbol

Name

Function

Fold Change

GGCAGAGCCT

Hs.1955

SAA2

Serum amyloid A2

Response to external stimulus

56

TTAACTTTAT

Hs.368626

RTN1

Reticulon 1

Signal Transduction

56

TAAGCCCCGC

Hs.2899

HPD

4-hydroxyphenylpyruvate dioxygenase

Cellular metabolism

49

GACACCGAGG

Hs.333383

FCN3

Ficolin (collagen/fibrinogen domain containing) 3 (Hakata antigen)

Transport

44

TTGGCTAGAC

Hs.100914

Cep192

Centrosomal protein 192 kDa

Unclassified

40

TTCCAGAGGC

Hs.8821

HAMP

Hepcidin antimicrobial peptide

Cellular metabolism

38

GATCCCAACT

Hs.418241

MT2A

Metallothionein 2A

Metal Ion Binding

37

AAGGACGCCG

Hs.117367

SLC22A1

Solute carrier family 22 (organic cation transporter), member 1

Transport

37

AACTCACAGA

Hs.168718

AFM

Afamin

Transport

36

AGAATAAGAA

Hs.418167

ALB

Albumin

Regulation of body fluid (Blood coagulation)

26

GCGGCCCCCC

Hs.839

IGFALS

Insulin-like growth factor binding protein, acid labile subunit

Signal Transduction

24

TACAGCCTGT

Hs.406184

FGFR1OP2

FGFR1 oncogene partner 2

Unclassified

24

CTTGGGTTTT

Hs.355888

PLCB2

Phospholipase C, beta 2

Lipid metabolism

24

CACTTCAAGG

Hs.521903

LY6E

Lymphocyte antigen 6 complex, locus E

Signal Transduction

23

CCACCCCGAA

Hs.35052

TEGT

Testis enhanced gene transcript (BAX inhibitor 1)

Cell Death

22

CCATTACCTC

Hs.134958

RODH-4

Retinol dehydrogenase 16 (all-trans and 13-cis)

Lipid metabolism

21

AGTCTGGCCT

Hs.416707

ABCA4

ATP-binding cassette, sub-family A (ABC1), member 4

Response to external stimulus

21

CCAGCAAGAG

Hs.11900

ZGPAT

Zinc finger, CCCH-type with G patch domain

Unclassified

20

AGAACCTTCC

Hs.181244

HLA-A

Major histocompatibility complex, class I, A

Antigene presentation

19

CCCTTCTTTG

Hs.356368

HAO2

Hydroxyacid oxidase 2 (long chain)

Cellular metabolism

19

Functional analysis of the 224 differently expressed genes was performed using DAVID program. The 104 up-regulated genes were grouped into eight functional clusters and one unclassified cluster (Additional file 2). These genes activated in HCC were involved in biological processes of biosynthesis, cell proliferation, signal transduction, transport, response to external stimulus, and cellular metabolism et al. Genes related to biosynthesis include 10 ribosomal proteins. Consistent with our observation, increased expression pattern of ribosomal proteins in HCC were demonstrated by previous studies [29, 30]. Genes participate in cell proliferation and signal transduction processes, such as MCM7 and IGFBP1 were also reported to be up-regulated in HCC by other groups [31, 32].

The 120 down-regulated genes were assembled into seven functional Gene Ontology (level 3) categories including blood coagulation, cell death, cellular metabolism, transport, signal transduction et al., and one unclassified group (Additional file 3). We also tested the over-representation of GO terms of all levels in this list. The most significant cluster is 8 genes belongs to a subcategory related to acute inflammatory response with p < 0.0005 after Benjamini correction of multiple testing. This cluster includes complement factor I (CFI) and several genes (C1S, C1QA, and C8B) encoding subcomponents of complement component 1 and 8. Another significantly enhanced cluster includes 6 genes related to blood coagulation (p < 0.03). These results indicated the loss of liver function in HCC.

Taken together, functional classification of genes differentially expressed in HCC revealed dysregulated pathways which may contribute to the disease. Similar to our results, microarray studies carried out by other groups also indicated genes involved in protein biosynthesis and cell signaling were up-regulated in HCC, while genes expressed at lower level in HCC included complement proteins and clotting factors [33, 34]. Thus, data obtained from SAGE and microarray methods could cross-validate each other and reveal some consistent events in HCC despite its heterogeneity.

Validation of differentially expressed genes

To validate the differently expressed genes identified by SAGE, real-time quantitative RT-PCR analysis was performed in 20 pairs of HCC and adjacent nontumorous liver tissues. Totally 10 genes from SAGE data were selected for examination. CCL20 and S100P were significantly up-regulated in HepG2 but not in HCC data, while ZFYVE1, PEG10, SGCE, XAGE-1, COL4A1, ZNF83, TNPO2 and TM4SF1 were up-regulated in HCC data. SAGE tag abundance of these genes was listed in Table 6. SGCE was not included in the list of 224 differently expressed genes which were derived only from single matched tags, because the tag matched to SGCE was a multiple matched tag. Notably, SGCE locates on chromosome 7q21 together with PEG10 in a head-to-head manner, separated by only 115 base pairs between the 5' ends of these two genes. Among these genes, CCL20 and PEG10 have been demonstrated to be up-regulated in HCC by other groups [3537], and were used as positive controls in our study. XAGE-1, COL4A1, S100P and TM4SF1 were observed to show elevated expression in various cancers, but the expression level of these genes in HCC has not been reported previously [3843]. The possible association existing between ZNF83, TNPO2, ZFYVE1 and cancer was discussed in the present study for the first time.
Table 6

SAGE tag abundance of selected genes

  

Tag Counts (Tag Per Million)

Tag

Unigene

Normal liver

HCC

HCC-Long

HepG2

GAGGGTTTAG

CCL20

25

23

5

844

TACCTCTGAT

S100P

25

23

5

1091

GAAGGTGATC

XAGE-1

5

224

286

5

GACCGCAGGA

COL4A1

5

188

146

5

ACCCGCCGGG

ZFYVE1

5

553

1107

203

ATTTTGTCGT

ZNF83

5

188

122

5

CTCATCCTAC

TNPO2

5

315

427

5

GAAGTTATAA

TM4SF1

25

169

685

499

GAAGTTATAA

PEG10

5

133

216

128

TTGGCAGTAT

SGCE; DDX43

25

242

263

30

By real-time quantitative RT-PCR analysis, increased expression was confirmed in 8 of these 10 genes. As shown in Figure 3, CCL20, S100P, PEG10, SGCE, XAGE-1, COL4A1, ZNF83, and TM4SF1 was significantly up-regulated in HCC compared to the corresponding nontumorous liver (p < 0.05). Expression level of PEG10 and SGCE was further examined in more samples (totally 32 pairs of HCC and corresponding nontumorous liver). Similar expression patterns of these two genes were observed in HCC, with Pearson's correlation coefficient 0.71 (Figure 4). No significant differences in expression level of TNPO2 and ZFYVE1 was detected between HCC and adjacent nontumorous liver, reflecting the heterogeneous feature of HCC. These results further validated the quantitative data of SAGE profiles and provided novel candidate genes which may help to illustrate the molecular mechanism of hepatocarcinogenesis.
https://static-content.springer.com/image/art%3A10.1186%2F1755-8794-2-5/MediaObjects/12920_2008_Article_69_Fig3_HTML.jpg
Figure 3

Real-time quantitative RT-PCR validation of SAGE data. Expression level of candidate genes identified to be up-regulated in HCC or HepG2 by SAGE data was examined in 20 pairs of HCC and corresponding nontumorous liver tissues by quantitative RT-PCR. Increased expression level in HCC was observed in CCL20, COL4A1, XAGE-1, PEG10, SGCE, S100P, TM4SF1 and ZNF83 genes (p < 0.05). Red bars represent expression level of HCC tissues, blue bars represent nontumorous liver tissues, and horizontal bars represent SD values.

https://static-content.springer.com/image/art%3A10.1186%2F1755-8794-2-5/MediaObjects/12920_2008_Article_69_Fig4_HTML.jpg
Figure 4

Similar expression patterns of PEG10 and SGCE in HCC patients. Expression level of PEG10 and SGCE was detected by real-time quantitative RT-PCR in 32 pairs of HCC and corresponding nontumorous liver tissues. The fold changes between HCC and nontumorous liver were represented by the height of bars (logarithmic scale). Red bars represent PEG10, and blue bars represent SGCE. These two genes showed similar expression patterns in 23 out of 32 patients, with Pearson's correlation coefficient 0.71.

Discussion

In this study, we applied SAGE to analyze gene expression profiles of normal liver, HepG2 and HCC on the genome scale. LongSAGE was also carried out to analyze the same HCC sample. From our results we found that the proportion of unique tags was higher in the HCC LongSAGE library (45%) than in the SAGE library (32%), suggesting that LongSAGE is definitely more efficient in detecting unique tags and tag-to-gene mapping than SAGE. Among the 20,313 unique tags in the HCC SAGE library, 6,815 (33.5%) of them were detected by LongSAGE tags with high correlation (Pearson's correlation 0.97). This indicated that the data of SAGE is repeatable even by LongSAGE. Thus it is reasonable that Pearson's correlation observed between our normal liver data and those in the public database was as high as 0.74 although they were obtained from totally independent experiments. We could draw a conclusion from these results that comparison among SAGE data produced by different labs is practical and reliable because of the high repeatability of SAGE technique.

It is not to our surprise that HepG2, a cultured hepatocellular carcinoma cell line, exhibited a different gene expression pattern compared to the pattern of HCC. Within passage cells, a lot of changes might have occurred on the genome scale which would subsequently affect the gene expression level. Relatively low correlation (Pearson's correlation coefficient 0.55) was observed between HepG2 and HCC in our SAGE data set. This result suggested that the expression profile of cultured cells should be cautiously used as reference of the corresponding tissue. However, analysis of HepG2 data still could provide some useful clues in identifying genes differentially expressed between normal liver and HCC. For example, increased expression of both CCL20 and S100P genes were detected in HCC by quantitative RT-PCR, although they were found to be significantly up-regulated only in HepG2 but not in HCC SAGE data. Consistent with our result, recent studies carried out by Rubie et al indicated that CCL20 was activated in HCC and supposed to be involved in hepatocarcinogenesis [35, 36]. Up-regulation of S100P was reported in breast cancer, prostate cancer and early-stage non-small cell lung cancer [3840], and it was also suggested to be a key factor in the aggressiveness of pancreatic cancer and promote cancer growth, survival and invasion [41]. Our results indicated that S100P was up-regulated in HCC, suggesting a possible role of S100P may play in liver cancer.

Enhanced expression of COL4A1, TM4SF1, XAGE-1, ZNF83, PEG10 and SGCE in HCC was detected by our SAGE data, and further validated by real-time quantitative RT-PCR. Potential roles of COL4A1, TM4SF1, and XAGE-1 in various cancers have been demonstrated including oral squamous cell carcinomas, gastric cancer, breast cancer and lung cancer et al, but not liver cancer [4244]. Our results gave a clue to the probability that these genes may also play a role in HCC pathogenesis. The maternally imprinted gene PEG10 was recently identified as a potential biomarker for HCC diagnosis, and genomic gain accounts for the major cause of its up-regulation in HCC [37]. SGCE was also a maternally imprinted gene, locating on chromosome 7q21 as close as 115 base pairs apart from PEG10 in a head-to-head manner. Due to the original SAGE tag TTGGCAGTAT matched to SGCE and DDX43 (DEAD box polypeptide 43) gene simultaneously, it was not included in our later analysis of differentially expressed genes in HCC which were derived only from single matched SAGE tags. However, PEG10 and SGCE were demonstrated to be up-regulated in a parallel way in B-cell chronic lymphocytic leukemia patients [45]. We examined the expression of these two genes in HCC and found that PEG10 and SGCE show similar expression patterns with close correlation (Pearson's correlation coefficient 0.71), pointing towards the possibility that these two genes may share common mechanism of deregulation in HCC. Whether PEG10 and SGCE are co-regulated via chromosome amplification at the 7q21 locus or their common promoter region is under investigating by our further study.

Conclusion

In summary, this study depicted the expression profile of HCC using both SAGE and LongSAGE techniques, which would provide valuable clues for understanding of molecular mechanism of HCC. Further investigations are needed for the candidate genes suggested by this study, especially for PEG10 and SGCE genes locating together on chromosome 7q21, and for many of the tags highly expressed in HCC but could not be matched to any known genes.

Notes

Abbreviations

HCC: 

hepatocellular carcinoma

SAGE: 

serial analysis of gene expression

Declarations

Acknowledgements

This study was supported by Shanghai Commission for Science and Technology (07ZR14083), the Scientific Research Programs Foundation during the Eleventh Five-Year Plan of PLA (06Q44), the Chinese High-Tech Research and Development Program (863) (2006AA020704, 2007AA021702) and Shanghai Rising-Star Program (07QB14015). Xijin Ge was supported in part by Agricultural Experiment Station of the South Dakota State University and Susan G. Komen Breast Cancer Foundation.

Authors’ Affiliations

(1)
Department of Microbiology and Microbial Engineering, School of Life Sciences, Fudan University
(2)
International Cooperation Laboratory on Signal Transduction, Eastern Hepatobiliary Surgery Institute, Second Military Medical University
(3)
Department of Mathematics and Statistics, South Dakota State University
(4)
National Engineering Center for Biochip at Shanghai
(5)
Chinese National Human Genome Center at Shanghai
(6)
Center for Clinical Pharmacology, Third Xiangya Hospital, Central South University
(7)
Department of Hepatobiliary Surgery, General Hospital of Air Force PLA

References

  1. Parkin DM, Pisani P, Ferlay J: Estimates of the worldwide incidence of 25 major cancers in 1990. Int J Cancer. 1999, 80: 827-841. 10.1002/(SICI)1097-0215(19990315)80:6<827::AID-IJC6>3.0.CO;2-P.View ArticlePubMedGoogle Scholar
  2. Parkin DM, Bray F, Ferlay J, Pisani P: Estimating the world cancer burden: Globocan 2000. Int J Cancer. 2001, 94: 153-156. 10.1002/ijc.1440.View ArticlePubMedGoogle Scholar
  3. Taylor-Robinson SD, Foster GR, Arora S, Hargreaves S, Thomas HC: Increase in primary liver cancer in the UK, 1979–94. Lancet. 1997, 350: 1142-1143. 10.1016/S0140-6736(05)63789-0.View ArticlePubMedGoogle Scholar
  4. Deuffic S, Poynard T, Buffat L, Valleron AJ: Trends in primary liver cancer. Lancet. 1998, 351: 214-215. 10.1016/S0140-6736(05)78179-4.View ArticlePubMedGoogle Scholar
  5. Yang L, Parkin DM, Ferlay J, Li L, Chen Y: Estimates of cancer incidence in China for 2000 and projections for 2005. Cancer Epidemiol Biomarkers Prev. 2005, 14: 243-250. 10.1158/1055-9965.EPI-04-0680.View ArticlePubMedGoogle Scholar
  6. Rashid A, Wang JS, Qian GS, Lu BX, Hamilton SR, Groopman JD: Genetic alterations in hepatocellular carcinomas: association between loss of chromosome 4q and p53 gene mutations. Br J Cancer. 1999, 80: 59-66. 10.1038/sj.bjc.6690321.View ArticlePubMedPubMed CentralGoogle Scholar
  7. Teramoto T, Satonaka K, Kitazawa S, Fujimori T, Hayashi K, Maeda S: p53 gene abnormalities are closely related to hepatoviral infections and occur at a late stage of hepatocarcinogenesis. Cancer Res. 1994, 54: 231-235.PubMedGoogle Scholar
  8. de La Coste A, Romagnolo B, Billuart P, Renard CA, Buendia MA, Soubrane O, Fabre M, Chelly J, Beldjord C, Kahn A, Perret C: Somatic mutations of the beta-catenin gene are frequent in mouse and human hepatocellular carcinomas. Proc Natl Acad Sci USA. 1998, 95: 8847-8851. 10.1073/pnas.95.15.8847.View ArticlePubMedGoogle Scholar
  9. Miyoshi Y, Iwao K, Nagasawa Y, Aihara T, Sasaki Y, Imaoka S, Murata M, Shimano T, Nakamura Y: Activation of the beta-catenin gene in primary hepatocellular carcinomas by somatic alterations involving exon 3. Cancer Res. 1998, 58: 2524-2527.PubMedGoogle Scholar
  10. Seidensticker MJ, Behrens J: Biochemical interactions in the wnt pathway. Biochim Biophys Acta. 2000, 1495: 168-182. 10.1016/S0167-4889(99)00158-5.View ArticlePubMedGoogle Scholar
  11. Abou-Elella A, Gramlich T, Fritsch C, Gansler T: c-myc amplification in hepatocellular carcinoma predicts unfavorable prognosis. Mod Pathol. 1996, 9: 95-98.PubMedGoogle Scholar
  12. Zhang XK, Huang DP, Qiu DK, Chiu JF: The expression of c-myc and c-N-ras in human cirrhotic livers, hepatocellular carcinomas and liver tissue surrounding the tumors. Oncogene. 1990, 5: 909-914.PubMedGoogle Scholar
  13. Suzuki K, Hayashi N, Yamada Y, Yoshihara H, Miyamoto Y, Ito Y, Ito T, Katayama K, Sasaki Y, Ito A: Expression of the c-met protooncogene in human hepatocellular carcinoma. Hepatology. 1994, 20: 1231-1236.View ArticlePubMedGoogle Scholar
  14. Twu JS, Lai MY, Chen DS, Robinson WS: Activation of protooncogene c-jun by the X protein of hepatitis B virus. Virology. 1993, 192: 346-350. 10.1006/viro.1993.1041.View ArticlePubMedGoogle Scholar
  15. Takada S, Koike K: Activated N-ras gene was found in human hepatoma tissue but only in a small fraction of the tumor cells. Oncogene. 1989, 4: 189-193.PubMedGoogle Scholar
  16. Hildt E, Saher G, Bruss V, Hofschneider PH: The hepatitis B virus large surface protein (LHBs) is a transcriptional activator. Virology. 1996, 225: 235-239. 10.1006/viro.1996.0594.View ArticlePubMedGoogle Scholar
  17. Yamashita T, Kaneko S, Hashimoto S, Sato T, Nagai S, Toyoda N, Suzuki T, Kobayashi K, Matsushima K: Serial analysis of gene expression in chronic hepatitis C and hepatocellular carcinoma. Biochem Biophys Res Commun. 2001, 282: 647-654. 10.1006/bbrc.2001.4610.View ArticlePubMedGoogle Scholar
  18. Patil MA, Chua MS, Pan KH, Lin R, Lih CJ, Cheung ST, Ho C, Li R, Fan ST, Cohen SN, Chen X, So S: An integrated data analysis approach to characterize genes highly expressed in hepatocellular carcinoma. Oncogene. 2005, 24: 3737-3747. 10.1038/sj.onc.1208479.View ArticlePubMedGoogle Scholar
  19. Kondoh N, Wakatsuki T, Ryo A, Hada A, Aihara T, Horiuchi S, Goseki N, Matsubara O, Takenaka K, Shichita M, Tanaka K, Shuda M, Yamamoto M: Identification and characterization of genes associated with human hepatocellular carcinogenesis. Cancer Res. 1999, 59: 4990-4996.PubMedGoogle Scholar
  20. Velculescu VE, Zhang L, Vogelstein B, Kinzler KW: Serial analysis of gene expression. Science. 1995, 270: 484-487. 10.1126/science.270.5235.484.View ArticlePubMedGoogle Scholar
  21. Saha S, Sparks AB, Rago C, Akmaev V, Wang CJ, Vogelstein B, Kinzler KW, Velculescu VE: Using the transcriptome to annotate the genome. Nat Biotechnol. 2002, 20: 508-512. 10.1038/nbt0502-508.View ArticlePubMedGoogle Scholar
  22. Velculescu VE, Zhang L, Zhou W, Vogelstein J, Basrai MA, Bassett DE, Hieter P, Vogelstein B, Kinzler KW: Characterization of the yeast transcriptome. Cell. 1997, 88: 243-251. 10.1016/S0092-8674(00)81845-0.View ArticlePubMedGoogle Scholar
  23. Lash AE, Tolstoshev CM, Wagner L, Schuler GD, Strausberg RL, Riggins GJ, Altschul SF: SAGEmap: a public gene expression resource. Genome Res. 2000, 10: 1051-1060. 10.1101/gr.10.7.1051.View ArticlePubMedPubMed CentralGoogle Scholar
  24. Romualdi C, Bortoluzzi S, D'Alessi F, Danieli GA: IDEG6: a web tool for detection of differentially expressed genes in multiple tag sampling experiments. Physiol Genomics. 2003, 12: 159-162.View ArticlePubMedGoogle Scholar
  25. Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7: 986-995.PubMedGoogle Scholar
  26. Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.View ArticlePubMedPubMed CentralGoogle Scholar
  27. Bishop JO, Morton JG, Rosbash M, Richardson M: Three abundance classes in HeLa cell messenger RNA. Nature. 1974, 250: 199-204. 10.1038/250199a0.View ArticlePubMedGoogle Scholar
  28. Man XB, Tang L, Zhang BH, Li SJ, Qiu XH, Wu MC, Wang HY: Upregulation of Glypican-3 expression in hepatocellular carcinoma but downregulation in cholangiocarcinoma indicates its differential diagnosis value in primary liver cancers. Liver Int. 2005, 25: 962-966. 10.1111/j.1478-3231.2005.01100.x.View ArticlePubMedGoogle Scholar
  29. Kondoh N, Shuda M, Tanaka K, Wakatsuki T, Hada A, Yamamoto M: Enhanced expression of S8, L12, L23a, L27 and L30 ribosomal protein mRNAs in human hepatocellular carcinoma. Anticancer Res. 2001, 21: 2429-2433.PubMedGoogle Scholar
  30. Shuda M, Kondoh N, Tanaka K, Ryo A, Wakatsuki T, Hada A, Goseki N, Igari T, Hatsuse K, Aihara T, Horiuchi S, Shichita M, Yamamoto N, Yamamoto M: Enhanced expression of translation factor mRNAs in hepatocellular carcinoma. Anticancer Res. 2000, 20: 2489-2494.PubMedGoogle Scholar
  31. Iizuka N, Tsunedomi R, Tamesa T, Okada T, Sakamoto K, Hamaguchi T, Yamada-Okabe H, Miyamoto T, Uchimura S, Hamamoto Y, Oka M: Involvement of c-myc-regulated genes in hepatocellular carcinoma related to genotype-C hepatitis B virus. J Cancer Res Clin Oncol. 2006, 132: 473-481. 10.1007/s00432-006-0094-8.View ArticlePubMedGoogle Scholar
  32. Kondoh N, Wakatsuki T, Ryo A, Hada A, Aihara T, Horiuchi S, Goseki N, Matsubara O, Takenaka K, Shichita M, Tanaka K, Shuda M, Yamamoto M: Identification and characterization of genes associated with human hepatocellular carcinogenesis. Cancer Res. 1999, 59: 4990-4996.PubMedGoogle Scholar
  33. Chen X, Cheung ST, So S, Fan ST, Barry C, Higgins J, Lai KM, Ji J, Dudoit S, Ng IO, Rijn Van De M, Botstein D, Brown PO: Gene expression patterns in human liver cancers. Mol Biol Cell. 2002, 13: 1929-1939. 10.1091/mbc.02-02-0023..View ArticlePubMedPubMed CentralGoogle Scholar
  34. Wurmbach E, Chen YB, Khitrov G, Zhang W, Roayaie S, Schwartz M, Fiel I, Thung S, Mazzaferro V, Bruix J, Bottinger E, Friedman S, Waxman S, Llovet JM: Genome-wide molecular profiles of HCV-induced dysplasia and hepatocellular carcinoma. Hepatology. 2007, 45: 938-947. 10.1002/hep.21622.View ArticlePubMedGoogle Scholar
  35. Rubie C, Frick VO, Wagner M, Rau B, Weber C, Kruse B, Kempf K, Tilton B, König J, Schilling M: Enhanced expression and clinical significance of CC-chemokine MIP-3 alpha in hepatocellular carcinoma. Scand J Immunol. 2006, 63: 468-477. 10.1111/j.1365-3083.2006.001766.x.View ArticlePubMedGoogle Scholar
  36. Rubie C, Frick VO, Wagner M, Weber C, Kruse B, Kempf K, König J, Rau B, Schilling M: Chemokine expression in hepatocellular carcinoma versus colorectal liver metastases. World J Gastroenterol. 2006, 12: 6627-6633.View ArticlePubMedPubMed CentralGoogle Scholar
  37. Ip WK, Lai PB, Wong NL, Sy SM, Beheshti B, Squire JA, Wong N: Identification of PEG10 as a progression related biomarker for hepatocellular carcinoma. Cancer Lett. 2007, 250: 284-291. 10.1016/j.canlet.2006.10.012.View ArticlePubMedGoogle Scholar
  38. Guerreiro Da Silva ID, Hu YF, Russo IH, Ao X, Salicioni AM, Yang X, Russo J: S100P calcium-binding protein overexpression is associated with immortalization of human breast epithelial cells in vitro and early stages of breast cancer development in vivo. Int J Oncol. 2000, 16: 231-240.PubMedGoogle Scholar
  39. Hammacher A, Thompson EW, Williams ED: Interleukin-6 is a potent inducer of S100P, which is up-regulated in androgen-refractory and metastatic prostate cancer. Int J Biochem Cell Biol. 2005, 37: 442-450. 10.1016/j.biocel.2004.07.011.View ArticlePubMedGoogle Scholar
  40. Diederichs S, Bulk E, Steffen B, Ji P, Tickenbrock L, Lang K, Zänker KS, Metzger R, Schneider PM, Gerke V, Thomas M, Berdel WE, Serve H, Müller-Tidow C: S100 family members and trypsinogens are predictors of distant metastasis and survival in early-stage non-small cell lung cancer. Cancer Res. 2004, 64: 5564-5569. 10.1158/0008-5472.CAN-04-2004.View ArticlePubMedGoogle Scholar
  41. Arumugam T, Simeone DM, Van Golen K, Logsdon CD: S100P promotes pancreatic cancer growth, survival, and invasion. Clin Cancer Res. 2005, 11: 5356-5364. 10.1158/1078-0432.CCR-05-0092.View ArticlePubMedGoogle Scholar
  42. Suhr ML, Dysvik B, Bruland O, Warnakulasuriya S, Amaratunga AN, Jonassen I, Vasstrand EN, Ibrahim SO: Gene expression profile of oral squamous cell carcinomas from Sri Lankan betel quid users. Oncol Rep. 2007, 18: 1061-1075.PubMedGoogle Scholar
  43. Kaneko R, Tsuji N, Kamagata C, Endoh T, Nakamura M, Kobayashi D, Yagihashi A, Watanabe N: Amount of expression of the tumor-associated antigen L6 gene and transmembrane 4 superfamily member 5 gene in gastric cancers and gastric mucosa. Am J Gastroenterol. 2001, 96: 3457-3458. 10.1111/j.1572-0241.2001.05355.x.View ArticlePubMedGoogle Scholar
  44. Egland KA, Kumar V, Duray P, Pastan I: Characterization of overlapping XAGE-1 transcripts encoding a cancer testis antigen expressed in lung, breast, and other types of cancers. Mol Cancer Ther. 2002, 1: 441-450.PubMedGoogle Scholar
  45. Kainz B, Shehata M, Bilban M, Kienle D, Heintel D, Krömer-Holzinger E, Le T, Kröber A, Heller G, Schwarzinger I, Demirtas D, Chott A, Döhner H, Zöchbauer-Müller S, Fonatsch C, Zielinski C, Stilgenbauer S, Gaiger A, Wagner O, Jäger U: Overexpression of the paternally expressed gene 10 (PEG10) from the imprinted locus on chromosome 7q21 in high-risk B-cell chronic lymphocytic leukemia. Int J Cancer. 2007, 121: 1984-1993. 10.1002/ijc.22929.View ArticlePubMedGoogle Scholar
  46. Pre-publication history

    1. The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/2/5/prepub

Copyright

© Dong et al; licensee BioMed Central Ltd. 2009

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Advertisement