- Open Access
Identification of frequent cytogenetic aberrations in hepatocellular carcinoma using gene-expression microarray data
Genome Biology volume 3, Article number: research0075.1 (2002)
Hepatocellular carcinoma (HCC) is a leading cause of death worldwide. Frequent cytogenetic abnormalities that occur in HCC suggest that tumor-modifying genes (oncogenes or tumor suppressors) may be driving selection for amplification or deletion of these particular genetic regions. In many cases, however, the gene(s) that drive the selection are unknown. Although techniques such as comparative genomic hybridization (CGH) have traditionally been used to identify cytogenetic aberrations, it might also be possible to identify them indirectly from gene-expression studies. A technique we have called comparative genomic microarray analysis (CGMA) predicts regions of cytogenetic change by searching for regional gene-expression biases. CGMA was applied to HCC gene-expression profiles to identify regions of frequent cytogenetic change and to identify genes whose expression is misregulated within these regions.
Using CGMA, 104 HCC gene-expression microarray profiles were analyzed. CGMA identified 13 regions of frequent cytogenetic change in the HCC samples. Ten of these regions have been detected in previous CGH studies (+lq, -4q, +6p, -8p, +8q, -13q, -16q, -17p, +17q, +20q). CGMA identified three additional regions that have not been previously identified by CGH (+5q, +12q, +19p). Genes located in regions of frequent cytogenetic change were examined for changed expression in the HCC samples.
Our results suggest that CGMA predictions using gene-expression microarray datasets are a practical alternative to CGH profiling. In addition, CGMA might be useful for identifying candidate genes within cytogenetically abnormal regions.
Aneuploidy is a common feature of cancer. Genetic alterations such as amplification, deletion, translocation and rearrangement could result in either gain-of-function or loss-of-function mutations in genes that modulate aspects of cell proliferation, differentiation, motility and survival. Whereas cytogenetic profiling techniques, such as comparative genomic hybridization (CGH) , have been useful in finding genetic abnormalities, other experimental approaches are frequently used to identify which specific gene(s) drive selection for the genetic aberration and contribute most to tumor progression. Common gene identification techniques include determining if a candidate gene contains a sequence mutation and/or determining if the candidate gene or gene product is abnormally expressed. As mutation analysis and protein expression studies are time-consuming, increasingly high-throughput gene-expression profiling is being used to identify abnormally expressed genes within a region of cytogenetic change [2,3,4,5,6].
Recently, several groups have observed that chromosomal changes can lead to regional biases in gene-expression values both in yeast (Saccharomyces cerevisiae) and in human tumors and tumor-derived cell lines [2,3,7,8]. These studies suggest that a fraction of gene-expression values (15-25%) are regulated in concordance with gene dosage. A computational technique termed comparative genomic microarray analysis (CGMA) has previously been used to identify regions of allelic imbalance indirectly from gene-expression profiles of human tumors . CGMA predicts chromosomal amplifications and deletions by organizing gene-expression data by genomic mapping location and scanning for regions that contain a statistically significant number of gene-expression values that change in the same relative direction. In this study, we apply CGMA analysis to a large hepatocellular carcinoma microarray dataset to demonstrate its validity as an alternative to CGH and to identify candidate genes in regions of frequent cytogenetic change.
Primary liver cancer in adults is the sixth most common form of cancer and the fourth leading cause of death from cancer worldwide [9,10]. Through the examination of hepatitis B virus (HBV) - and hepatitis C virus (HCV)-induced tumors, two landmark CGH studies have suggested that a subset of cytogenetic changes frequently occurs in HCC [11,12]. These include frequent gain of chromosomes 1q, 6p, 8q, 17q and 20q and frequent loss of chromosomes 1p, 4q, 6q, 8p, 13, 16 and 17p [11,12]. In particular, gain of chromosomes 1q and 8q has been associated with the early development of HCC , whereas loss of chromosome 4q has been linked to increased aggressiveness of established tumors . To determine whether gene-expression data could be used to identify cytogenetic changes accurately, we applied CGMA to a microarray dataset of HCC tumors and compared the CGMA predictions to existing CGH data. For HCC, CGMA was able to predict nearly all chromosomal aberrations identified previously by CGH. In addition, from the gene-expression data we also identified a set of genes whose expression values change most within the regions of cytogenetic change. These genes may represent candidate genes whose expression changes drive selection for chromosomal gains or losses.
CGMA predictions of cytogenetic changes
Normalized, log-transformed gene-expression data from 104 HCC gene-expression arrays  were obtained from the Stanford Microarray Database . As CGMA analysis yields more intuitive predictions if the tumor expression data is compared to a normal tissue reference, the original HCC gene-expression data was mathematically transformed such that the pooled cell-line reference was replaced by a normal tissue reference (, see Materials and methods). Using this transformation, gene-expression values from the tumor sample are compared to corresponding values from non-cancerous tissue. Genomic regions that contain a disproportionate number of genes that change expression in the same relative direction (that is, show a gene-expression bias) may indicate an underlying chromosomal gain or loss (Figure 1a) [2,3,7,8]. Chromosomal regions that contained a gene-expression bias with at least 95% confidence (a sign test z-statistic of at least 1.96, see Materials and methods) were identified for all 104 HCC expression profiles (Figure 1b). In addition, genomic regions that contained significant gene-expression biases in at least 35% of non-replicate samples were identified (Figure 1c). A 35% threshold was chosen because in previous CGMA profiling experiments this threshold yielded the highest CGMA to CGH agreement ( and data not shown). CGMA predicted frequent gains for chromosome 1q (gained in 72% of tumor samples), 6p (56%), 8q (49%), 17q (46%), 20q (46%), 5q (42%), 19P (37%) and 12q (35%). Frequent chromosomal losses were predicted for chromosome 4q (lost in 66% of tumor samples), 17p (48%), 13 (39%), 16 (37%), and 8p (35%). To determine if CGMA predictions were consistent with other cytogenetic profiling studies, the CGMA data were compared with data from two of the largest CGH profiling studies (67 and 50 samples, respectively) using HCC tumors [11,12]. Of the 13 regions detected by CGMA, 10 (77%) were also implicated by CGH (Figure 1c). CGMA also detected three gained regions, chromosomes 5q, 12q and 19p, which were not implicated in the CGH analysis. CGMA failed to discover two regions of loss detected by CGH - chromosomes 1p and 6q. It is noteworthy that these particular losses were not identified in both CGH studies. These data suggest that CGMA predictions produce results very similar to traditional CGH profiling studies.
Comparison to previous HCC studies
To date there have been at least 20 reports on the application of CGH to HCC [11,12,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]. To determine whether the differences between the CGMA predictions and the two large CGH studies were similar to the experimental variation observed between different HCC CGH studies, predictions produced by CGMA were compared to 13 different HCC CGH profiling studies (Figure 2). CGMA produced 10 of 13 (77%) predictions that matched a consensus chromosomal aberration profile. On average, each CGH study matched the consensus profile in 78 ± 14% of the chromosomal regions analyzed. Therefore, the variation in CGMA results was similar to the variations between independent CGH studies. These results suggest that CGMA profiling is able to predict regions of frequent chromosomal imbalance in HCC as well as CGH profiling. As only 4 of the 104 (4%) samples analyzed scored positive for HCV infection, we could not use this dataset to detect significant cytogenetic differences between HCV-infected verses HBV-infected individuals (data not shown).
CGMA predictions of multifocal tumors
Included in the set of HCC gene-expression profiles were several cases in which multiple tumor nodules were removed from the same patient. In some cases the nodules had related gene-expression profiles (patients HK63, HK64, HK66) whereas in other cases tumors from the same patient had distinctive profiles (HK65, HK67, HK85) . In particular, in patient HK67 the gene-expression profile from nodule HK67.1 was distinct from the expression profiles from nodules HK67.2 and HK67.3 . Array CGH was used previously to determine the cytogenetic profiles of the tumor nodules from this patient . Array CGH identified common cytogenetic abnormalities in patient HK67's tumors, including loss of chromosome 15q, an unusual gain of chromosome 19q, and loss of the centromeric region of chromosome 22. However, additional cytogenetic changes were found in tumors HK67.2 and HK67.3 that were not present in HK67.1. The cytogenetic profiling data coupled with the observation that HK67.2 and HK67.3 were both smaller in size and had an increased mitotic index, suggested that HK67.1 was the primary tumor and HK67.2 and HK67.3 were divergent HK67.1 subclones. To determine whether CGMA predictions agree with this monoclonal origin hypothesis, CGMA profiles of patient HK67's tumor nodules were isolated and organized by hierarchical clustering (Figure 3). CGMA also detected the common chromosome 19q gain and chromosome 15q loss in the HK67 tumors. CGMA did not identify a common loss on chromosome 22, however, CGMA identified other genetic aberrations (+8q, -16q and -19q) consistently found in the HK67 tumors. CGMA also identified additional aberrations (+2q, +5q, +12q) present in H67.2 and H67.3 that were not found in H67.1. Taken together, the CGMA data supports the hypothesis that the H67.2 and HK67.3 tumor nodules probably arose from HK67.1, but that additional distinct cytogenetic events had occurred in these nodules during tumor progression. In contrast, tumor nodules from patient HK64 have very similar gene-expression profiles and very similar cytogenetic profiles as predicted by CGMA, suggesting that these tumors have common origins and these nodules have not diverged significantly from the original lineage. In addition, tumor nodules from patient HK85 showed distinctive expression profiles and distinct HBV integration sites . Similarly, the tumors from patient HK85 also show distinct CGMA-predicted cytogenetic profiles, reflecting the independent transforming mechanism (Figure 3).
Identification of candidate genes in regions of frequent cytogenetic change
Frequent cytogenetic abnormalities suggest that tumor-modifying genes (oncogenes or tumor suppressors) may be driving selection for the amplification or deletion of these particular genetic regions [6,11,12,34,35,36]. An advantage of using CGMA profiling rather than traditional molecular genetic profiling is that access to gene-expression data is inherent in the analysis. CGMA allows cytogenetic analysis and the candidate gene approach to be performed with the same dataset. For example, the c-myc oncogene has been postulated to drive selection for frequent chromosome 8q amplification. Though c-myc is located on a region that both CGMA and CGH identify as frequently gained, c-myc's expression is increased more than twofold in less than 6% of the samples. In fact, in 52% of the HCC tissue samples, c-myc's expression is downregulated (Table 1). This implies that increased c-myc expression is not driving the selection for the amplification of chromosomal region 8q in these samples. In the small region of chromosome 8q presented in Table 1, two other genes (for squalene monooxygenase and pro2000) do show increased expression in a majority of HCC samples. Consistent with previous reports examining gene-expression levels in regions of cytogenetic change, expression levels for a large percentage of genes in this amplified region remain unchanged [3,5,8].
The set of genes that are consistently misregulated in regions of frequent cytogenetic change as predicted by CGMA are shown in Table 2. Platelet-derived growth factor receptor alpha is consistently downregulated in a region of frequent cytogenetic loss and this suggests that loss of a member of this receptor gene family is important in HCC progression. It has previously been reported that a transcript (PRLTS) with sequence similarity to the extracellular domain of platelet-derived growth factor receptor may also be a tumor suppressor for HCC . In addition, consistently increased expression of the pituitary tumor transforming gene 1 oncogene (PTTG) is observed in these samples (Table 2). PTTG maps to chromosome 5q, a region that was identified as frequently changed by CGMA, but not identified in the majority of CGH profiling studies. PTTP overexpression in NIH 3T3 cells induces these cells to form tumors when injected into nude mice. Overexpression of this gene may result from frequent chromosomal amplification and may participate in HCC tumor progression.
CGMA prediction software
To assist in identifying regions of unidirectional gene-expression bias, we have constructed a web-based program that processes two-color gene-expression data and identifies genomic regions that contain gene-expression biases. The input for this program is a simple tab-delimited gene-expression matrix file consisting of columns for the probe sequence identifier, probe name, and gene-expression ratios. Because different microarray technologies use different identifiers to describe the microarray probe, the program translates probe sequence identifiers (ids) such as GenBank accession numbers and UniGene cluster ids to Ensembl transcript ids using precompiled sequence comparisons. After data analysis, a summary table is displayed showing chromosomal regions that show significant (α ≤ 0.05) unidirectional gene-expression bias highlighted in either red or green, indicating either increased or decreased expression biases, respectively. The program can also send several output files to the user via e-mail. These files include a summary report that contains the z-statistic for each chromosomal region (positive for upregulated regions and negative for downregulated regions) and a list of genes located in regions of frequent cytogenetic change. The program is available at .
In this study we have used gene-expression profiling data to predict cytogenetic changes that frequently occur in HCC. Two landmark CGH analyses identified 12 different regions of frequent imbalance. However, one study found 8 regions and the other study found 11 [11,12]. Five of these 12 regions were not found in both experiments. CGMA successfully identified 10 of 12 regions previously distinguished by CGH. CGMA also detected three regions that have not been implicated by these CGH studies. On average however, 22% of genomic regions indentified in a particular HCC CGH study are not constantly identified in other studies. Therefore, the three inconsistent CGMA predictions (3 of 13; 23%) are comparable to the inconsistencies between independent CGH studies for HCC.
Three additional regions were identified by CGMA that were not identified by CGH. While these CGMA-predicted regions were near the 35% cutoff for detection, they could represent other regions of allelic imbalance yet to be detected by CGH. It is also possible that biological mechanisms other than cytogenetic change could influence expression in large genomic regions and produce regional gene-expression biases. Additional molecular genetic work will be required to resolve these differences.
If CGH data are not available for a particular cancer type, but gene-expression profiling data are, then CGMA could allow rapid prediction of the cytogenetic abnormalities that frequently occur within that cancer type. Moreover, in instances where gene-expression profiling reveals previously unrecognized cancer subtypes, CGMA could determine whether cytogenetic differences are responsible for these different subtypes. In cancer types where traditional cytogenetic profiling has already been carried out, CGMA predictions could serve to confirm existing cytogenetic profiling data and be used further to examine candidate genes whose expression changes most within a region of frequent cytogenetic change. In this way CGMA can be combined with the candidate gene approach to identify genes that are directly involved in tumor progression.
CGMA can be used to indicate chromosomal imbalances by detecting chromosomal regions that contain a disproportionate number of gene-expression values that change in the same relative direction. This analysis provides good evidence that CGMA is a practical alternative to CGH cytogenetic profiling when gene-expression profiling data is available.
Materials and methods
Normalization and filtering
Normalized, log-transformed gene-expression data for 104 HCC samples and 76 corresponding non-cancerous liver gene-expression profiles  were obtained from the Stanford Microarray Database . Genes that were present in at least 75% of samples (10,037 genes) were used for further analysis. In this study, both the tumor samples and normal tissue samples were compared to a pooled cell-line reference . To allow comparison of tumor gene-expression values to gene-expression values from surrounding non-cancerous tissue, new gene expression ratios, tumor verse normal (T/N), were estimated. To create the new ratios, log-transformed non-cancerous tissue ratios (N/U) were subtracted from the log-transformed HCC tissue ratios (T/U) for each gene such that log2(T/N) = log2(T/U) - log2(N/U). If an HCC sample did not have a corresponding non-cancerous sample, the global mean of the non-cancerous tissue gene-expression ratios were used.
To identify regional gene-expression biases, gene-expression values that map within a given chromosomal arm were collected and a sign test for a one-sample mean/median was used to determine whether a significant upward or downward bias was present in the expression values. An exception was made for chromosomes 13-16, 21 and 22. These chromosomes are more telocentric and therefore only their q-arms were tested for expression biases. Sequence comparisons were used to map microarray probe sequences (the sequences that are placed on the microarray) to predicted Ensembl transcripts (Ensembl version 6.28) . Included in the Ensembl transcript annotations are chromosomal mapping locations at base-pair resolution. Redundancy introduced by replicate probes on the array and/or multiple probes mapping to the same gene were eliminated by averaging expression values that map to identical transcripts. Of the filtered set of 10,037 genes, 6,274 genes (63%) were unique and had associated genomic mapping information.
A sign test for the one-sample mean/median was used to determine whether a significant number of genes that map to a given chromosomal region change in a unidirectional manner. The algorithm scores a gene as up (+) or down (-) regulated if the magnitude of the expression value change is at least 1.8-fold. The sign test computes the probability, in the form of a z-statistic, of finding x upregulated genes out of n genes that change in a given genomic region. For simplicity, the z-statistic is computed using the normal approximation to the binomial distribution such that z = (2x - n)/sqrt(n). Genomic regions that contained less than 15 changed gene-expression values were excluded from further analysis. On average, 160 gene-expression values were located to each genomic region. The sign test z-statistic can be converted to a significance value (α) based on the two-tailed z-statistic (zα/2) critical values. For example, if z = 1.96, then α = 0.05; if z = 2.58 then α = 0.01, and so on.
Kallioniemi A, Kallioniemi OP, Sudar D, Rutovitz D, Gray JW, Waldman F, Pinkel D: Comparative genomic hybridization for molecular cytogenetic analysis of solid tumors. Science. 1992, 258: 818-821.
Virtaneva K, Wright FA, Tanner SM, Yuan B, Lemon WJ, Caligiuri MA, Bloomfield CD, de La Chapelle A, Krahe R: Expression profiling reveals fundamental biological differences in acute myeloid leukemia with isolated trisomy 8 and normal cytogenetics. Proc Natl Acad Sci USA. 2001, 98: 1124-1129. 10.1073/pnas.98.3.1124.
Phillips JL, Hayward SW, Wang Y, Vasselli J, Pavlovich C, Padilla-Nash H, Pezullo JR, Ghadimi BM, Grossfeld GD, Rivera A, et al: The consequences of chromosomal aneuploidy on gene expression profiles in a cell line model for prostate carcinogenesis. Cancer Res. 2001, 61: 8143-8149.
Monni O, Barlund M, Mousses S, Kononen J, Sauter G, Heiskanen M, Paavola P, Avela K, Chen Y, Bittner ML, Kallioniemi A: Comprehensive copy number and gene expression profiling of the 17q23 amplicon in human breast cancer. Proc Natl Acad Sci USA. 2001, 98: 5711-5716. 10.1073/pnas.091582298.
Platzer P, Upender MB, Wilson K, Willis J, Lutterbaugh J, Nosrati A, Willson JK, Mack D, Ried T, Markowitz S: Silence of chromosomal amplifications in colon cancer. Cancer Res. 2002, 62: 1134-1138.
Xu XR, Huang J, Xu ZG, Qian BZ, Zhu ZD, Yan Q, Cai T, Zhang X, Xiao HS, Qu J, et al: Insight into hepatocellular carcinogenesis at transcriptome level by comparing gene expression profiles of hepatocellular carcinoma with those of corresponding non-cancerous liver. Proc Natl Acad Sci USA. 2001, 98: 15089-15094. 10.1073/pnas.241522398.
Hughes TR, Roberts CJ, Dai H, Jones AR, Meyer MR, Slade D, Burchard J, Dow S, Ward TR, Kidd MJ, et al: Widespread aneuploidy revealed by DNA microarray expression profiling. Nat Genet. 2000, 25: 333-337. 10.1038/77116.
Haddad R, Furge KA, Miller J, Schoumans J, Haab B, Teh B, Barr L, Webb C: Genomic profiling and cDNA microarray analysis of human colon adenocarcinoma and associated peritoneal metastasis reveals consistent cytogenetic and transcriptional aberrations associated with progression of multiple metastases. Appl Genomics Proteomics. 2002, 1: 123-134.
Pisani P, Parkin DM, Ferlay J: Estimates of the worldwide mortality from eighteen major cancers in 1985. Implications for prevention and projections of future burden. Int J Cancer. 1993, 55: 891-903.
Parkin DM, Pisani P, Ferlay J: Estimates of the worldwide incidence of eighteen major cancers in 1985. Int J Cancer. 1993, 54: 594-606.
Marchio A, Meddeb M, Pineau P, Danglot G, Tiollais P, Bernheim A, Dejean A: Recurrent chromosomal abnormalities in hepatocellular carcinoma detected by comparative genomic hybridization. Genes Chromosomes Cancer. 1997, 18: 59-65. 10.1002/(SICI)1098-2264(199701)18:1<59::AID-GCC7>3.0.CO;2-0.
Wong N, Lai P, Lee SW, Fan S, Pang E, Liew CT, Sheng Z, Lau JW, Johnson PJ: Assessment of genetic changes in hepatocellular carcinoma by comparative genomic hybridization analysis: relationship to disease stage, tumor size, and cirrhosis. Am J Pathol. 1999, 154: 37-43.
Chen X, Cheung ST, So S, Fan ST, Barry C, Higgins J, Lai K, Dudoit S, Ng I, Rijn M, et al: Gene expression patterns in human liver cancer. Mol Biol Cell. 2002, 13: 1929-1939. 10.1091/mbc.02-02-0023..
Sherlock G, Hernandez-Boussard T, Kasarskis A, Binkley G, Matese JC, Dwight SS, Kaloper M, Weng S, Jin H, Ball CA, et al: The Stanford Microarray Database. Nucleic Acids Res. 2001, 29: 152-155. 10.1093/nar/29.1.152.
Wong N, Lai P, Pang E, Fung LF, Sheng Z, Wong V, Wang W, Hayashi Y, Perlman E, Yuna S, et al: Genomic aberrations in human hepatocellular carcinomas of differing etiologies. Clin Cancer Res. 2000, 6: 4000-4009.
Kusano N, Shiraishi K, Kubo K, Oga A, Okita K, Sasaki K: Genetic aberrations detected by comparative genomic hybridization in hepatocellular carcinomas: their relationship to clinicopathological features. Hepatology. 1999, 29: 1858-18562.
Marchio A, Pineau P, Meddeb M, Terris B, Tiollais P, Bernheim A, Dejean A: Distinct chromosomal abnormality pattern in primary liver cancer of non-B, non-C patients. Oncogene. 2000, 19: 3733-3738. 10.1038/sj.onc.1203713.
Guan XY, Fang Y, Sham JS, Kwong DL, Zhang Y, Liang Q, Li H, Zhou H, Trent JM: Recurrent chromosome alterations in hepatocellular carcinoma detected by comparative genomic hybridization. Genes Chromosomes Cancer. 2000, 29: 110-116. 10.1002/1098-2264(2000)9999:9999<::AID-GCC1022>3.0.CO;2-V.
Zondervan PE, Wink J, Alers JC, Ijzermans JN, Schalm SW, de Man RA, van Dekken H: Molecular cytogenetic evaluation of virus-associated and non-viral hepatocellular carcinoma: analysis of 26 carcinomas and 12 concurrent dysplasias. J Pathol. 2000, 192: 207-215. 10.1002/1096-9896(2000)9999:9999<::AID-PATH690>3.0.CO;2-#.
Tornillo L, Carafa V, Richter J, Sauter G, Moch H, Minola E, Gambacorta M, Bianchi L, Vecchione R, Terracciano LM: Marked genetic similarities between hepatitis B virus-positive and hepatitis C virus-positive hepatocellular carcinomas. J Pathol. 2000, 192: 307-312. 10.1002/1096-9896(2000)9999:9999<::AID-PATH706>3.0.CO;2-O.
Collonge-Rame MA, Bresson-Hadni S, Koch S, Carbillet JP, Blagosklonova O, Mantion G, Miguet JP, Heyd B, Bresson JL: Pattern of chromosomal imbalances in non-B virus related hepatocellular carcinoma detected by comparative genomic hybridization. Cancer Genet Cytogenet. 2001, 127: 49-52. 10.1016/S0165-4608(00)00421-0.
Niketeghad F, Decker HJ, Caselmann WH, Lund P, Geissler F, Dienes HP, Schirmacher P: Frequent genomic imbalances suggest commonly altered tumour genes in human hepatocarcinogenesis. Br J Cancer. 2001, 85: 697-704. 10.1054/bjoc.2001.1963.
Koo SH, Ihm CH, Kwon KC, Park JW, Kim JM, Kong G: Genetic alterations in hepatocellular carcinoma and intrahepatic cholangiocarcinoma. Cancer Genet Cytogenet. 2001, 130: 22-28. 10.1016/S0165-4608(01)00460-5.
Kitay-Cohen Y, Amiel A, Ashur Y, Fejgin MD, Herishanu Y, Afanasyev F, Bomstein Y, Lishner M: Analysis of chromosomal aberrations in large hepatocellular carcinomas by comparative genomic hybridization. Cancer Genet Cytogenet. 2001, 131: 60-64. 10.1016/S0165-4608(01)00492-7.
Chang J, Kim NG, Piao Z, Park C, Park KS, Paik YK, Lee WJ, Kim BR, Kim H: Assessment of chromosomal losses and gains in hepatocellular carcinoma. Cancer Lett. 2002, 182: 193-202. 10.1016/S0304-3835(02)00083-6.
Balsara BR, Pei J, De Rienzo A, Simon D, Tosolini A, Lu YY, Shen FM, Fan X, Lin WY, Buetow KH, et al: Human hepatocellular carcinoma is characterized by a highly consistent pattern of genomic imbalances, including frequent loss of 16q23.1-24.1. Genes Chromosomes Cancer. 2001, 30: 245-253. 10.1002/1098-2264(2000)9999:9999<::AID-GCC1083>3.0.CO;2-M.
Zimonjic DB, Keck CL, Thorgeirsson SS, Popescu NC: Novel recurrent genetic imbalances in human hepatocellular carcinoma cell lines identified by comparative genomic hybridization. Hepatology. 1999, 29: 1208-1214.
Harada T, Shiraishi K, Kusano N, Umayahara K, Kondoh S, Okita K, Sasaki K: Evaluation of the reliability of chromosomal imbalances detected by combined use of universal DNA amplification and comparative genomic hybridization. Jpn J Cancer Res. 2000, 91: 1119-1125.
Rao UN, Gollin SM, Beaves S, Cieply K, Nalesnik M, Michalopoulos GK: Comparative genomic hybridization of hepatocellular carcinoma: correlation with fluorescence in situ hybridization in paraffin-embedded tissue. Mol Diagn. 2001, 6: 27-37. 10.1054/modi.2001.22021.
Takeo S, Arai H, Kusano N, Harada T, Furuya T, Kawauchi S, Oga A, Hirano T, Yoshida T, Okita K, Sasaki K: Examination of oncogene amplification by genomic DNA microarray in hepatocellular carcinomas: comparison with comparative genomic hybridization analysis. Cancer Genet Cytogenet. 2001, 130: 127-132. 10.1016/S0165-4608(01)00479-4.
Kusano N, Okita K, Shirahashi H, Harada T, Shiraishi K, Oga A, Kawauchi S, Furuya T, Sasaki K: Chromosomal imbalances detected by comparative genomic hybridization are associated with outcome of patients with hepatocellular carcinoma. Cancer. 2002, 94: 746-751. 10.1002/cncr.10254.
Yasui K, Arii S, Zhao C, Imoto I, Ueda M, Nagai H, Emi M, Inazawa J: TFDP1, CUL4A, and CDC16 identified as targets for amplification at 13q34 in hepatocellular carcinomas. Hepatology. 2002, 35: 1476-1484.
Chen X, Cheung ST, So S, Fan ST, Barry C, Higgins J, Lai KM, Ji J, Dudoit S, Ng IO, et al: Gene expression patterns in human liver cancers. Mol Biol Cell. 2002, 13: 1929-1939. 10.1091/mbc.02-02-0023..
Wong N, Lam WC, Lai PB, Pang E, Lau WY, Johnson PJ: Hypomethylation of chromosome 1 heterochromatin DNA correlates with q-arm copy gain in human hepatocellular carcinoma. Am J Pathol. 2001, 159: 465-471.
Fujiwara Y, Ohata H, Kuroki T, Koyama K, Tsuchiya E, Monden M, Nakamura Y: Isolation of a candidate tumor suppressor gene on chromosome 8p21.3-p22 that is homologous to an extracellular domain of the PDGF receptor beta gene. Oncogene. 1995, 10: 891-895.
Kuroki T, Fujiwara Y, Nakamori S, Imaoka S, Kanematsu T, Nakamura Y: Evidence for the presence of two tumour-suppressor genes for hepatocellular carcinoma on chromosome 13q. Br J Cancer. 1995, 72: 383-385.
CGMA - comparative analysis of microarray data. [http://www.vai.org/vari/downloads/development/clam-0.2.pl]
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.
Funding was through the generosity of the Van Andel Research Institute and Michigan Center for Biological Information (MCBI). We specially thank Ramsi Haddad for helpful discussion in the preparation of this manuscript.
About this article
Cite this article
Crawley, J.J., Furge, K.A. Identification of frequent cytogenetic aberrations in hepatocellular carcinoma using gene-expression microarray data. Genome Biol 3, research0075.1 (2002) doi:10.1186/gb-2002-3-12-research0075
- Comparative Genomic Hybridization
- Pituitary Tumor Transforming Gene
- Comparative Genomic Hybridization Analysis
- Cytogenetic Change
- Comparative Genomic Hybridization Data