- Research article
- Open Access
Gene expressions and copy numbers associated with metastatic phenotypes of uterine cervical cancer
© Lyng et al; licensee BioMed Central Ltd. 2006
- Received: 06 June 2006
- Accepted: 20 October 2006
- Published: 20 October 2006
A better understanding of the development of metastatic disease and the identification of molecular markers for cancer spread would be useful for the design of improved treatment strategies. This study was conducted to identify gene expressions associated with metastatic phenotypes of locally advanced cervical carcinomas and investigate whether gains or losses of these genes could play a role in regulation of the transcripts. Gene expressions and copy number changes were determined in primary tumors from 29 patients with and 19 without diagnosed lymph node metastases by use of cDNA and genomic microarray techniques, respectively.
Thirty-one genes that differed in expression between the node positive and negative tumors were identified. Expressions of eight of these genes (MRPL11, CKS2, PDK2, MRPS23, MSN, TBX3, KLF3, LSM3) correlated with progression free survival in univariate analysis and were therefore more strongly associated with metastatic phenotypes than the others. Immunohistochemistry data of CKS2 and MSN showed similar relationships to survival. The prognostic genes clustered into two groups, suggesting two major metastatic phenotypes. One group was associated with rapid proliferation, oxidative phosphorylation, invasiveness, and tumor size (MRPS23, MRPL11, CKS2, LSM3, TBX3, MSN) and another with hypoxia tolerance, anaerobic metabolism, and high lactate content (PDK2, KLF3). Multivariate analysis identified tumor volume and PDK2 expression as independent prognostic variables. Gene copy number changes of the differentially expressed genes were not frequent, but correlated with the expression level for seven genes, including MRPS23, MSN, and LSM3.
Gene expressions associated with known metastatic phenotypes of cervical cancers were identified. Our findings may indicate molecular mechanisms underlying development of these phenotypes and be useful as markers of cancer spread. Gains or losses of the genes may be involved in development of the metastatic phenotypes in some cases, but other mechanisms for transcriptional regulation are probably important in the majority of tumors.
- Cervical Cancer
- Negative Tumor
- Copy Number Change
- Metastatic Phenotype
- Hypoxia Tolerance
Lymph node involvement is the first indication of spreading and a strong prognostic factor for epithelial cancers . A better understanding of the development of metastatic tumor phenotypes and the identification of molecular markers for lymphatic spread would be useful in design of improved treatment strategies . Cervical carcinomas have been studied extensively during the last years in the search for biological characteristics that are associated with lymph node involvement. Tumor volume is among the strongest prognostic factors. Severe hypoxia, high level of lactate, high proliferation rate, increased angiogenesis, high interstitial fluid pressure, and low apoptotic activity have also been associated with poor prognosis [3–9], suggesting that these are metastasis promoting phenotypes. The molecular biology behind these phenotypes has not been clarified, but the proliferation proteins EGFR and ERBB2, the anti-apoptosis proteins cIAP1 and BCL2, and the glucose transporter GLUT1 may be involved [10–14]. A more comprehensive characterization of the tumors is, however, needed to achieve a complete understanding of how metastatic phenotypes develop in cervical cancers.
Gene expression microarrays are useful for discovery of new genes that are regulated in metastatic tumors [2, 15]. Panels of genes associated with lymph node metastasis have been identified for several cancer types by using this technique [2, 16], suggesting that the gene expression program is altered in metastatic as compared to nonmetastatic tumors. Genes related to treatment outcome or lymph node metastasis have also been identified in microarray studies of cervical cancers [17–19]. Less than 20 patients were included in these studies, making it difficult to draw firm conclusions from the results. Moreover, protein expressions or copy numbers of the genes were not addressed.
In the present work we have used microarrays to identify gene expressions associated with metastatic phenotypes of cervical cancers and to investigate whether gains or losses of these genes could play a role in the transcriptional regulation and phenotype development. Locally advanced primary tumors of squamous cell origin, all receiving curative radiotherapy, were included. Treatment outcome differs considerably among these patients, emphasizing the need for identifying and exploring risk factors. We report 31 genes that differed in expression between lymph node positive and negative tumors, as diagnosed from magnetic resonance (MR) images. These genes provided a basis for further analyses and enabled us to generate hypotheses of gene functions in metastatic tumors. The importance of gains and losses for the expression of these genes were explored by including genomic microarray data in the analyses. The gene data were further related to progression free survival, since this end point is a stronger indicator of the metastatic capacity of the tumors than the diagnosed metastatic status per se. Unsupervised hierarchical clustering was performed to identify coregulated genes that were probably associated with the same phenotypes.
Differentially expressed genes in node positive versus node negative tumors
Genes differentially expressed in node positive versus node negative cervical tumors
Upregulated in node positive tumors
DnaJ (Hsp(40) homolog, subfamily C, member 9
Mitochondrial ribosomal protein S23
Structural component of mitochondrial ribosome
Von Willebrand factor
Lysosomal proteinase inhibition
Pyruvate dehydrogenase kinase, isoenzyme 2
CDC28 protein kinase regulatory subunit 2
Cell cycle regulation
In multiple clusters
Hypothetical protein MGC14151
Hypothetical protein LOC285378
Mitochondrial ribosomal protein L11
Structural component of mitochondrial ribosome
Brain-specific angiogenesis inhibitor 3
RNA terminal phosphate cyclase-like 1
Maturation of 18S RNA
NIMA (never in mitosis gene a) – related kinase 1
DNA damage repair
Hypothetical protein DKFZp586I1420
Downregulated in node positive tumors
Erythrocyte membrane protein band 4.1 like 4B
Structural constituent of cytoskeleton
Degradation of extracellular matrix components
Muscleblind-like 2 (Drosophila)
Hypothetical protein FLJ13291
Krüppel-like factor 3 (basic)
Glucose metabolism, growth, apoptosis, angiogenesis
In multiple clusters
ERO1-like (S. cerevisiae)
Oxidative protein folding
Cell-cell and -extracellular matrix interaction
T-box 3 (ulnar mammary syndrome)
Regulation of growth and developmental processes
Growth regulation, apoptosis, calcium signaling
In multiple clusters
Cell-cell and -extracellular matrix interaction
Myocyte enhancer factor 2A
Growth factor related transcription, differentiation
Gene copy numbers
Progression free survival
Cox regression analysis of progression free survival
Model containing clinical variables alone
1.01 – 1.02
1.01 – 1.02
Number of pathological lymph nodes
1.1 – 1.7
0.64 – 4.3
Model containing gene variables alone
2.3 – 14
1.9 – 13
1.3 – 7.9
1.4 – 10
1.1 – 6.6
0.15 – 0.90
0.09 – 0.67
1.1 – 6.9
0.10 – 0.92
0.15 – 0.93
0.97 – 7.4
Model containing clinical and gene variables
1.01 – 1.02
1.3 – 8.3
Univariate analysis of the gene copy number changes, categorised as gain, loss, or no change, showed a significant relationship to progression free survival for MRPS23 (p = 0.0002), PDK2 (p = 0.004), and MSN (p = 0.016). Although the aberrations were not frequent, occurring in 5 (MRPS23), 7 (PDK2), and 14 (MSN) tumors, they seemed to have a pronounced influence on the survival probability when they occurred (Figure 4B). The relationship for MRPS23 and PDK2 was similar, except for two additional losses of PDK2, reflecting their colocation on chromosome 17q. The three tumors with gain of PDK2 had a high PDK2 expression, however, there was no general correlation between expression and aberration for this gene (Figure 2B). The results in Figure 4B for PDK2 could therefore be a consequence of its colocalization with MRPS23. The results of MRPS23 and MSN, on the other hand, were consistent with our observations that the aberrations correlated with the expression levels of these genes (Figure 2B) and that the expression levels showed prognostic significance (Figure 4A).
In univariate analysis of the protein expressions the immunostaining score of the tumor cells was included as continuous data, without considering the score in stroma cells. The data showed a significant relationship to progression free survival for CKS2 (p = 0.034) and MSN (p = 0.025), but not for CSTA, MEF2A, and HK2. The largest difference in survival was achieved when using a cut off for the immunostaining score of 4 and 9 for CKS2 and MSN, respectively (Figure 4C). The protein data were therefore in concordance with the gene expression results, showing that high expression of CKS2 and low expression of MSN were associated with poor survival (Figure 4A).
To ensure that the expression data of each phenotype separated the patients into groups with different survival probability, unsupervised clustering of the patients was performed. The combined expression data (log2ratios) of each group of coregulated genes were used in two separate analyses. Clustering based on MSN, TBX3, LSM3, CKS2, MRPL1, and MRPS23 expression identified four tumor groups, for which patients of group 1 and 4, with high expression of LSM3, CKS2, MRPL11, and MRPS23 and, in general, low expression of MSN and TBX3 compared to group 2 and 3, had the lowest survival probability (Figure 5B). Seven of the eight tumors with gain of MRPS23 or loss of MSN were in these groups. A similar analysis based on KLF3 and PDK2 expression, separated two major tumor groups, for which patients in the group with high PDK2 and low KLF3 expression had the lowest survival probability (Figure 5C). All tumors with gain of PDK2 were in this group. Each of the two phenotypes indicated by the coregulated genes in Figure 5A were therefore associated with poor progression free survival probability and therefore with metastasis development.
Genes that differed in expression between node positive and negative cervical tumors and therefore may be related to metastatic phenotypes, were identified in our study. Our data on protein expressions and gene copy numbers provided information on the cell type expressing the genes and the regulation mechanisms involved. The frequent copy number changes, especially on chromosome 1q, 3q, 3p, and 5p, were consistent with previous reports [20, 21]. Copy number changes of the differentially expressed genes were, however, less common. Such changes do probably not play a role in development of the metastatic phenotypes in the majority of tumors, for which other transcriptional regulation mechanisms seem to be important. It should be noted that no general conclusion about the role of gene copy number changes in development of the metastatic phenotypes could be drawn from our study, since only selected genes were considered. Gains or losses of other genes may be important and even influence the expression of the genes addressed here.
The protein data were not correlated with the gene expressions of two of the five proteins investigated, HK2 and MEF2A. Cross-reaction of the antibodies used for immunohistochemistry to other proteins may explain this apparent discrepancy. Hence, although the MEF2A antibody used was recommended for this protein, the producer states that cross-reactions to MEF2C and MEF2D may occur to a lesser extent. The protein data may also be less representative of the entire tumor than the gene expressions, since they were derived from a single biopsy whereas several samples were used in the microarray analyses. Moreover, post-transcriptional control of the protein levels is a likely explanation of these results as well .
We identified two independent groups of genes with prognostic significance, suggesting the existence of at least two major metastatic phenotypes of the locally advanced stages of cervical carcinomas. None of the prognostic genes have previously been associated with metastasis in this tumor type. Genes, such as EGFR, ERBB2, BCL2, cIAP, and GLUT1, which have shown correlations to survival in protein studies [10–14], were, however, not identified here. None of these were differentially expressed between the node positive and negative tumors with the cut off used in our study and therefore not considered in the further analyses. A separate analysis showed that EGFR expression correlated with survival (data not shown), in concordance with previous reports . The other genes may be regulated post-transcriptionally or be prognostic in larger studies. Our strategy was therefore not suited to find all prognostic genes, but ensured that the ones identified were truly associated with metastatic phenotypes. Based on the current knowledge of gene function it was possible to relate the prognostic genes to biological characteristics that have been associated with metastasis development in cervical cancers and thereby propose molecular mechanisms underlying these characteristics.
MRPS23, MRPL11, CKS2, TBX3, LSM3, and MSN constituted the largest group of coregulated and prognostic genes. The former five genes probably play a role in tumor growth, since their expression showed a significant relationship to tumor volume. MRPS23 and MRPL11 are structural components of the mitochondria. Upregulation of these genes suggests increased mitochondrial activity, energy production through oxidative phosphorylation, and a high oxygen consumption rate. CKS2 is involved in cell cycle control, and its activation has been associated with high proliferation of lymphoma cells and incidence of metastasis in colon carcinomas [23, 24]. LSM3 and TBX3 may also participate in cell cycle control; upregulation of LSM3 promotes proliferation of pancreatic cancer cells, and TBX3 may interact with the cell cycle protein CDKN2A [25, 26]. Our data are therefore consistent with increased mitochondrial activity and cellular proliferation in metastatic tumors, probably leading to rapid growth and large volumes. These phenotypes have previously been associated with poor prognosis of cervical cancer , and our findings point to genes that are associated with and might participate in their development. Furthermore, when gains and losses of MRPS23 and LSM3 occur, they may influence the gene expressions and possibly the development of metastatic disease.
Clustering of MSN together with MRPS23, MRPL11, CKS2, TBX, and LSM suggests a role of this gene in cell proliferation, consistent with previous studies . Loss of MSN function may also cause cell depolarization, increased motility, and invasiveness, and therefore be directly involved in several steps during spreading . MSN repression may therefore actively promote metastasis development in some rapidly growing tumors. Our data further suggest that gene loss leads to MSN downregulation and is associated with a metastatic phenotype. Loss of heterozygosity studies have shown increased imbalance on the X chromosome in lymph node metastases as compared to the primary cervical tumors, in agreement with this hypothesis .
The other group of coregulated and prognostic genes, PDK2 and KLF3, is probably associated with hypoxia tolerance. Activation of PDK2 directs glycolytically derived pyruvate towards anaerobic metabolism and lactate output rather than oxidative phosphorylation through repression of the pyruvate dehydrogenase (PDH) complex . Repression of KLF3 may indicate increased glucose metabolism under hypoxia, although other functions involving growth, apoptosis, and angiogenesis, have also been proposed [30, 31]. Our data on PDK2 and KLF3 therefore suggest activation of strategies to conserve energy and, hence, survive under hypoxia. Hypoxia and high lactate content are prognostic factors for cervical cancer [3, 5, 8], consistent with our hypothesis. Moreover, repression of PDH subunit alpha 1 has previously been associated with treatment outcome of this disease , and similar survival strategies are probably used in other tumor types, showing activation of lactate dehydrogenases [32, 33]. Our data showed that tumors with concurrent gain of PDK2 and MRPS23 had high expression of MRPS23 and PDK2, large volumes above the median value, a very short survival time and, therefore, an extremely high risk for metastasis development, suggesting the combined characteristics of rapid proliferation, high oxygen consumption, and hypoxia tolerance.
The regulation of the other differentially expressed genes also contributes to the characteristics of metastatic disease. Many of them were coregulated with the prognostic ones, and could be markers of a high proliferation activity or regulated in response to hypoxia. Repression of NTN4 and HYAL1 may promote cell migration and invasive growth, whereas BAI3, VWF, and EPB41L4B probably participate in angiogenesis, attachment of tumor cells to endothelial surfaces, or reflect vascular structures in the tumors [34–36]. Moreover, regulation of NEK1, CSTA, ANXA4, and DDOST indicates activation of DNA damage repair (NEK1) and resistance to apoptosis (CSTA and possibly ANXA4 and DDOST) [37–39], whereas transcriptional downregulation of HK2 may be a result of glucose deprivation . The roles of MBNL2, RCL1, MGC14151, ERO1L, DNAJC9, PLAC2, and MEF2A are more unclear, but may involve regulation of proteins like the insulin receptor (MBNL2)  or proteins participating in development (MEF2A). Elevated VWF plasma levels have been found in node positive cervical cancer patients , consistent with our results. Moreover, increased CSTA protein expression has been correlated with poor prognosis of breast carcinoma, whereas repression of ERO1L has been shown to increase the recurrence probability of pulmonary adenocarcinoma [42, 43], suggesting that these genes are related to metastasis development also in other tumor types.
We have identified genes associated with major metastatic phenotypes of cervical cancers that probably are related to rapid proliferation and hypoxia. Copy number changes of the genes may be involved in development of these phenotypes in some cases, but other mechanisms for transcriptional regulation are probably important in the majority of tumors. The gene expressions may be useful markers of cancer spread, since they showed a stronger relationship to progression free survival than the metastatic status and number of pathological lymph nodes and added prognostic information to tumor volume. Although our findings need to be confirmed in larger studies with more patients, they provide a useful basis for further investigations to reveal the mechanisms underlying development of these phenotypes.
Patients, lymph node involvement, and treatment
Forty-eight patients with primary squamous cell carcinoma of the uterine cervix, diagnosed in the period 2001 – 2004 at Health Enterprise Rikshospitalet-Radiumhospitalet, were included. Tumor stage (The Fédération Internationale des Gynaecologistes et Obstetristes) was 2a (1 patient), 2b (29), 3b (14), and 4a (4), and tumor grade was 2 or 3. Lymph nodes were evaluated at the time of diagnosis with MR imaging, using axial fast spin echo T2-weighted images. A pelvic phased array coil was used for nodes below the promontory. The body coil was used for retroperitoneal nodes from the promontory to the upper pole of the kidneys. Lymph nodes were classified as pathologic when the minimal axial diameter exceeded 10 mm in oval nodes or 8 mm in round nodes. Using these criteria for metastasis in pelvic lymph nodes in prostatic and urinary bladder carcinoma, the sensitivity and specificity were 75% and 98%, respectively . Some tumors may therefore falsely be classified as node negative. Pretreatment tumor volume, determined from the MR images as π/6·x·y·z, where x is largest diameter in the sagittal plane through the cervix length axis and y and z are diameters in a plane orthogonal to this axis, ranged from 17 – 336 cm3 (median 46 cm3).
All patients received external irradiation and brachytherapy with curative intent. External irradiation, 50 Gy to tumor and parametria and 45 Gy to the rest of the pelvic region, was delivered by use of a linear accelerator in 25 fractions five times per week. A 192Ir stepping source was used for endocavitary brachytherapy, employing 21Gy in five fractions to point A. Adjuvant cisplatin (40 mg/m2) was given weekly in maximum 6 courses during the period of external radiation. Twenty-two patients completed all cisplatin courses, whereas the others had dose reduction or delay due to toxicity problems. The follow up included clinical examinations every 3rd month for the first 2 years, twice a year the next three years, and thereafter once a year. MR imaging of pelvis and retroperitoneum and X-ray of thorax were performed when symptoms of recurrent disease were seen. Progression free survival, defined as the time between diagnosis and the first event of locoregional and/or distant relapse or cancer related death, was used as end point. Four patients, who died due to cerebral hemorrhage or treatment related complications, were censored. All other deaths were cancer related. Observation time, calculated for the patients that were still alive, ranged from 15 – 51 months (median 36 months). The study was approved by the regional committee of medical research ethics in southern Norway (REK no. S-01129), and written informed-consent was achieved from all patients.
One – four biopsies (median 3 biopsies), approximately 5 × 5 × 5 mm in size, were used in the microarray analyses of each tumor, minimizing confounding effects caused by intratumor heterogeneity in gene expressions and copy numbers . The biopsies were taken from different locations of the tumor at the time of diagnosis, immediately snap-frozen in liquid nitrogen, and stored at -80°C. All biopsies had more than 50% tumor cells in hematoxylin and eosin stained sections, derived from the central part of the specimen. Median tumor cell fraction was 70% (range 50–90%) for both the node positive and negative tumor groups. Total RNA was isolated by use of Trizol reagent (Invitrogen, Carlsbad, CA) followed by double precipitation with isopropanol and final precipitation with 5 M lithium chloride . Genomic DNA was isolated according to a standard protocol, including proteinase K, phenol, chloroform, and isoamylalcohol . A separate biopsy was fixed in neutral 4% buffered formalin and used for immunohistochemistry.
Gene expression microarray analysis
Gene expressions were studied using array slides produced at the Microarray Facility at Health Enterprise Rikshospitalet-Radiumhospitalet, containing 15000 cDNA clones. Labeled cDNA was synthesized from 20 μg total RNA by anchored oligo(dT)-primed reverse transcription, using SuperScript II reverse transcriptase (Invitrogen) and Fairplay labeling kit (Stratagene, La Jolla, CA) in the presence of either Cy3-dUTP or Cy5-dUTP (Amersham Pharmacia Biotech Inch., Piscataway, NJ). The labeled cDNA of each tumor was co-hybridized with that of a reference sample (Universal Human Reference RNA, Stratagene) to the array slides overnight at 65°C, using an automated hybridization station (GeneTAC, Genomic Solutions/Perkin Elmer, Boston, MA). All hybridizations were performed twice in a dye-swap design. Scanning and image analysis were performed with an Agilent scanner (Agilent Technologies Inc., Palo Alto, CA) and the GenePix 4.1 image analysis software (Axon Instruments Inc., Union City, CA), respectively. Data preprosessing included correction of saturated intensities , filtering of weak and bad spots, and lowess normalization. Our protocol for the cDNA microarray experiments has revealed reliable results in agreement with northern analyses in previous studies on cell lines .
Quantitative real time PCR
For validation of the microarray data we performed qRT PCR of four genes (HK2, KLF3, MRPS23, PDK2) in all 48 tumors and two genes (CSTA, DDOST) in twelve tumors. Pre-designed, gene-specific TaqMan probe and primer sets (Applied Biosystems, Foster City, CA) and a 7500 Sequence Detector (Perkin-Elmer) were used (Additional file 4). Ten ng cDNA, which was synthesized from total RNA by use of Superscript II transcriptase (Invitrogen), was employed. Conditions for amplification were one cycle of 95°C in 10 min followed by 40 cycles of 95°C in 15 sec and 60°C in 1 min. The reactions were carried out in triplicate in a 25 μl reaction volume and a 96-well plate format. Gene expression was calculated using the delta-delta method, were the B2M gene served as an internal control and the reference sample used in the microarray experiments served as a calibrator .
Genomic microarray analysis
Gains and losses of the genes identified in the cDNA microarray analysis were assessed using genomic array slides produced by the Microarray Facility at Health Enterprise Rikshospitalet-Radiumhospitalet, containing 5000 unique BAC clones (RPCI-11 library, Wellcome Trust Sanger Institute, Cambridge, UK) that covered the whole genome with a resolution of approximately 1 Mbp. The BAC clone localizations were obtained from Ensembl . The clone identifications corresponding to the genes of interest are listed in Additional file 2. Forty-seven of the 48 tumors were included, for which genomic DNA was available. One half μg of digested and ethanol-purified DNA and normal female reference DNA were labeled by a random primer reaction with Cy3-dUTP and Cy5-dUTP (Amersham Pharmacia Biotech Inch.), respectively, and co-hybridized to the array slides in 48 hours at 37°C by use of an automated hybridization station (GeneTAC, Genomic Solutions/Perkin Elmer) and our established protocol . Scanning, image analysis, and normalization were performed as described previously, using an Agilent scanner (Agilent Technologies Inc.), the GenePix 4.1 image analysis software (Axon Instruments Inc.), and the M-CGH data preprosessing software, respectively . Gene copy number changes were determined from the log2ratios by using a cut off value specific for each tumor, taking into account the tumor cell fraction, the DNA ploidy, and the empirically determined dynamics of the experiments (Figure 2A). This procedure provides results in agreement with fluorescence-in situ-hybridization data [50, 51]. The data were classified as gain, loss, or no change in copy number relative to the modal DNA content.
Commercially available antibodies for CKS2, CSTA, HK2, MEF2A, and MSN were used. Paraffin tissue sections of all tumors were stained using the Dako EnVision™ + System, Peroxidase (DAB) (K4007 and K4011, Dako Corporation, CA) and Dakoautostainer. Monoclonal antibodies MSN (clone 38/87, Sigma, MI), CSTA (clone WR-23/2/3/3, Sigma), CKS2 (clone 1F7G5, Zymed Laboratories Inc., CA), and polyclonal antibodies HK2 (AB3279, Chemicon International, CA) and MEF2A (C-21, Santa Cruz Biotechnology, CA) were used. Detailed information on the immunohistochemical staining has been published previously . All series included positive and negative controls that gave satisfactory results. Negative controls included substitution of the monoclonal antibody with mouse myeloma protein of the same subclass and concentration as the monoclonal antibody and substitution of the polyclonal primary antibody with normal rabbit IgG of the same concentration as the polyclonal antibody. Immunostaining was scored on a 3-tiered scale for both intensity (absent/weak, 1; moderate, 2; strong, 3) and extent of staining (percentage of positive tumor cells < 10%, 1; 10–50%, 2; > 50%, 3). The scoring results of intensity and extent were multiplied to give a composite score ranging from 1 to 9 for each tumor.
The SAM (Significance Analysis of Microarrays) program  with a false discovery rate of 10% was applied on the log2ratios to find the genes that differed most in expression between the node positive and negative tumors. This method does not take into account the increased information and certainty in the data achieved from dye-swap designs with two experiments on each tumor. The analysis was therefore performed on the raw data sets. The differentially expressed genes were selected and the average expression ratios were calculated from the log2ratios of the two dye-swap experiments for further use.
The data of the differentially expressed genes were included in Cox univariate and multivariate analysis of progression free survival, using the SPSS software. When the variables are semi continuous, the covariate values may be any uniform function of these. We tested three functions for the gene expression data: the log2ratio, the ratio, and the categorized form, and genes showing significance (p < 0.05) for one of these functions were identified as major metastasis genes. In the analysis of categorized data, the tumors were divided into three groups; one group of 12 tumors with the highest expression ratios, another group of 24 tumors with intermediate ratios, and a third group of 12 tumors with the lowest ratios. The survival probability of the three groups were compared in Kaplan Meier plots to ensure that there was a monotonous increase or decrease in survival among the groups. The categorized form was chosen as approximately the group of 12 tumors versus the others. Unsupervised hierarchical clustering of the gene expression data was performed based on the log2ratios, using a program developed for this purpose .
Array express accession
The raw data and the processed data file from the cDNA microarray platform have been deposited to the Array Express repository (E-TABM-146).
We thank Mai Thi Phuong Nguyen, Esther Hasper and Ann Larsen at Department of Pathology for technical assistance. The study was supported by The National Programme for Research in Functional Genomics (FUGE) in the Research Council of Norway.
- Eccles SA: Cell biology of lymphatic metastasis. The potential role of c-erbB oncogene signalling. Recent Results Cancer Res. 2000, 157: 41-54.PubMedView ArticleGoogle Scholar
- Roepman P, Wessels LF, Kettelarij N, Kemmeren P, Miles AJ, Lijnzaad P, Tilanus MG, Koole R, Hordijk GJ, van der Vliet PC, Reinders MJ, Slootweg PJ, Holstege FC: An expression profile for diagnosis of lymph node metastases from primary head and neck squamous cell carcinomas. Nat Genet. 2005, 37 (2): 182-186. 10.1038/ng1502.PubMedView ArticleGoogle Scholar
- Fyles A, Milosevic M, Hedley D, Pintilie M, Levin W, Manchul L, Hill RP: Tumor hypoxia has independent predictor impact only in patients with node-negative cervix cancer. J Clin Oncol. 2002, 20 (3): 680-687. 10.1200/JCO.20.3.680.PubMedView ArticleGoogle Scholar
- Hockel M, Schlenger K, Hockel S, Vaupel P: Hypoxic cervical cancers with low apoptotic index are highly aggressive. Cancer Res. 1999, 59 (18): 4525-4528.PubMedGoogle Scholar
- Lyng H, Sundfor K, Trope C, Rofstad EK: Disease control of uterine cervical cancer: relationships to tumor oxygen tension, vascular density, cell density, and frequency of mitosis and apoptosis measured before treatment and during radiotherapy. Clin Cancer Res. 2000, 6 (3): 1104-1112.PubMedGoogle Scholar
- Milosevic M, Fyles A, Hedley D, Pintilie M, Levin W, Manchul L, Hill R: Interstitial fluid pressure predicts survival in patients with cervix cancer independent of clinical prognostic factors and tumor oxygen measurements. Cancer Res. 2001, 61 (17): 6400-6405.PubMedGoogle Scholar
- Tsang RW, Wong CS, Fyles AW, Levin W, Manchul LA, Milosevic M, Chapman W, Li YQ, Pintilie M: Tumour proliferation and apoptosis in human uterine cervix carcinoma II: correlations with clinical outcome. Radiother Oncol. 1999, 50 (1): 93-101. 10.1016/S0167-8140(98)00119-4.PubMedView ArticleGoogle Scholar
- Walenta S, Wetterling M, Lehrke M, Schwickert G, Sundfor 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 (4): 916-921.PubMedGoogle Scholar
- West CM, Cooper RA, Loncaster JA, Wilks DP, Bromley M: Tumor vascularity: a histological measure of angiogenesis and hypoxia. Cancer Res. 2001, 61 (7): 2907-2910.PubMedGoogle Scholar
- Airley R, Loncaster J, Davidson S, Bromley M, Roberts S, Patterson A, Hunter R, Stratford I, West C: Glucose transporter glut-1 expression correlates with tumor hypoxia and predicts metastasis-free survival in advanced carcinoma of the cervix. Clin Cancer Res. 2001, 7 (4): 928-34.PubMedGoogle Scholar
- Crawford RA, Caldwell C, Iles RK, Lowe D, Shepherd JH, Chard T: Prognostic significance of the bcl-2 apoptotic family of proteins in primary and recurrent cervical cancer. Br J Cancer. 1998, 78 (2): 210-214.PubMedPubMed CentralView ArticleGoogle Scholar
- Imoto I, Tsuda H, Hirasawa A, Miura M, Sakamoto M, Hirohashi S, Inazawa J: Expression of cIAP1, a target for 11q22 amplification, correlates with resistance of cervical cancers to radiotherapy. Cancer Res. 2002, 62 (17): 4860-4866.PubMedGoogle Scholar
- Kim YT, Park SW, Kim JW: Correlation between expression of EGFR and the prognosis of patients with cervical carcinoma. Gynecol Oncol. 2002, 87 (1): 84-89. 10.1006/gyno.2002.6803.PubMedView ArticleGoogle Scholar
- Nakano T, Oka K, Ishikawa A, Morita S: Correlation of cervical carcinoma c-erb B-2 oncogene with cell proliferation parameters in patients treated with radiation therapy for cervical carcinoma. Cancer. 1997, 79 (3): 513-520. 10.1002/(SICI)1097-0142(19970201)79:3<513::AID-CNCR13>3.0.CO;2-6.PubMedView ArticleGoogle Scholar
- Ramaswamy S, Ross KN, Lander ES, Golub TR: A molecular signature of metastasis in primary solid tumors. Nat Genet. 2003, 33 (1): 49-54. 10.1038/ng1060.PubMedView ArticleGoogle Scholar
- Bidus MA, Risinger JI, Chandramouli GV, Dainty LA, Litzi TJ, Berchuck A, Barrett JC, Maxwell GL: Prediction of lymph node metastasis in patients with endometrioid endometrial cancer using expression microarray. Clin Cancer Res. 2006, 12 (1): 83-88. 10.1158/1078-0432.CCR-05-0835.PubMedView ArticleGoogle Scholar
- Harima Y, Togashi A, Horikoshi K, Imamura M, Sougawa M, Sawada S, Tsunoda T, Nakamura Y, Katagiri T: Prediction of outcome of advanced cervical cancer to thermoradiotherapy according to expression profiles of 35 genes selected by cDNA microarray analysis. Int J Radiat Oncol Biol Phys. 2004, 60 (1): 237-248. 10.1016/j.ijrobp.2004.02.047.PubMedView ArticleGoogle Scholar
- Wong YF, Selvanayagam ZE, Wei N, Porter J, Vittal R, Hu R, Lin Y, Liao J, Shih JW, Cheung TH, Lo KW, Yim SF, Yip SK, Ngong DT, Siu N, Chan LK, Chan CS, Kong T, Kutlina E, McKinnon RD, Denhardt DT, Chin KV, Chung TK: Expression genomics of cervical cancer: molecular classification and prediction of radiotherapy response by DNA microarray. Clin Cancer Res. 2003, 9 (15): 5486-5492.PubMedGoogle Scholar
- Grigsby PW, Watson M, Powell MA, Zhang Z, Rader JS: Gene expression patterns in advanced human cervical cancer. Int J Gynecol Cancer. 2006, 16 (2): 562-567. 10.1111/j.1525-1438.2006.00389.x.PubMedView ArticleGoogle Scholar
- Heselmeyer K, Macville M, Schröck E, Blegen H, Hellström A-C, Shah K, Auer G, Ried T: Advanced-stage cervical carcinomas are defined by a recurrent pattern of chromosomal aberrations reveling high genetic instability and a consistent gain of chromoksome arm 3q. Genes Chromosomes Cancer. 1997, 19: 233-240. 10.1002/(SICI)1098-2264(199708)19:4<233::AID-GCC5>3.0.CO;2-Y.PubMedView ArticleGoogle Scholar
- Lyng H, Beigi M, Svendsrud DH, Brustugun OT, Stokke T, Kristensen GB, Sundfor K, Skjonsberg A, De Angelis PM: Intratumor chromosomal heterogeneity in advanced carcinomas of the uterine cervix. Int J Cancer. 2004, 111 (3): 358-366. 10.1002/ijc.20258.PubMedView ArticleGoogle Scholar
- Pedersen PL, Mathupala S, Rempel A, Geschwind JF, Ko YH: Mitochondrial bound type II hexokinase: a key player in the growth and survival of many cancers and an ideal prospect for therapeutic intervention. Biochim Biophys Acta. 2002, 1555 (1–3): 14-20.PubMedView ArticleGoogle Scholar
- Urbanowicz-Kachnowicz I, Baghdassarian N, Nakache C, Gracia D, Mekki Y, Bryon PA, Ffrench M: ckshs expression is linked to cell proliferation in normal and malignant human lymphoid cells. Int J Cancer. 1999, 82 (1): 98-104. 10.1002/(SICI)1097-0215(19990702)82:1<98::AID-IJC17>3.0.CO;2-A.PubMedView ArticleGoogle Scholar
- Li M, Lin YM, Hasegawa S, Shimokawa T, Murata K, Kameyama M, Ishikawa O, Katagiri T, Tsunoda T, Nakamura Y, Furukawa Y: Genes associated with liver metastasis of colon cancer, identified by genome-wide cDNA microarray. Int J Oncol. 2004, 24 (2): 305-312.PubMedGoogle Scholar
- Fraser MM, Watson PM, Fraig MM, Kelley JR, Nelson PS, Boylan AM, Cole DJ, Watson DK: CaSm-mediated cellular transformation is associated with altered gene expression and messenger RNA stability. Cancer Res. 2005, 65 (14): 6228-6236. 10.1158/0008-5472.CAN-05-0650.PubMedView ArticleGoogle Scholar
- Packham EA, Brook JD: T-box genes in human disorders. Hum Mol Genet. 2003, 12 (Spec No 1): R37-44. 10.1093/hmg/ddg077.PubMedView ArticleGoogle Scholar
- McClatchey AI: Merlin and ERM proteins: unappreciated roles in cancer development?. Nat Rev Cancer. 2003, 3 (11): 877-883. 10.1038/nrc1213.PubMedView ArticleGoogle Scholar
- Kersemaekers AM, van de Vijver MJ, Kenter GG, Fleuren GJ: Genetic alterations during the progression of squamous cell carcinomas of the uterine cervix. Genes Chromosomes Cancer. 1999, 26 (4): 346-354. 10.1002/(SICI)1098-2264(199912)26:4<346::AID-GCC9>3.0.CO;2-D.PubMedView ArticleGoogle Scholar
- Sugden MC, Holness MJ: Recent advances in mechanisms regulating glucose oxidation at the level of the pyruvate dehydrogenase complex by PDKs. Am J Physiol Endocrinol Metab. 2003, 284 (5): E855-E862.PubMedView ArticleGoogle Scholar
- Black AR, Black JD, Azizkhan-Clifford J: Sp1 and kruppel-like factor family of transcription factors in cell growth regulation and cancer. J Cell Physiol. 2001, 188 (2): 143-160. 10.1002/jcp.1111.PubMedView ArticleGoogle Scholar
- Hasan NM, MacDonald MJ: Sp/Kruppel-like transcription factors are essential for the expression of mitochondrial glycerol phosphate dehydrogenase promoter B. Gene. 2002, 296 (1–2): 221-234. 10.1016/S0378-1119(02)00865-X.PubMedView ArticleGoogle Scholar
- Ha CS, Kong JS, McLaughlin P, Tucker SL, Fayad LE, Hess MA, Wilder RB, Cabanillas F, Cox JD: Stage III follicular lymphoma: long-term follow-up and patterns of failure. Int J Radiat Oncol Biol Phys. 2003, 57 (3): 748-754. 10.1016/S0360-3016(03)00641-2.PubMedView ArticleGoogle Scholar
- Koukourakis MI, Giatromanolaki A, Sivridis E, Bougioukas G, Didilis V, Gatter KC, Harris AL: Lactate dehydrogenase-5 (LDH-5) overexpression in non-small-cell lung cancer tissues is linked to tumour hypoxia, angiogenic factor production and poor prognosis. Br J Cancer. 2003, 89 (5): 877-885. 10.1038/sj.bjc.6601205.PubMedPubMed CentralView ArticleGoogle Scholar
- Bertrand P, Courel MN, Maingonnat C, Jardin F, Tilly H, Bastard C: Expression of HYAL2 mRNA, hyaluronan and hyaluronidase in B-cell non-Hodgkin lymphoma: relationship with tumor aggressiveness. Int J Cancer. 2005, 113 (2): 207-212. 10.1002/ijc.20562.PubMedView ArticleGoogle Scholar
- Kennedy TE: Cellular mechanisms of netrin function: long-range and short-range actions. Biochem Cell Biol. 2000, 78 (5): 569-575. 10.1139/bcb-78-5-569.PubMedView ArticleGoogle Scholar
- Schwartsmann G, Damin DC, Roisemberg I: Malignant disease and von Willebrand factor. Lancet Oncol. 2001, 2 (12): 716-717. 10.1016/S1470-2045(01)00580-0.PubMedView ArticleGoogle Scholar
- Hayes MJ, Moss SE: Annexins and disease. Biochem Biophys Res Commun. 2004, 322 (4): 1166-1170. 10.1016/j.bbrc.2004.07.124.PubMedView ArticleGoogle Scholar
- Kelleher DJ, Gilmore R: DAD1, the defender against apoptotic cell death, is a subunit of the mammalian oligosaccharyltransferase. Proc Natl Acad Sci USA. 1997, 94 (10): 4994-4999. 10.1073/pnas.94.10.4994.PubMedPubMed CentralView ArticleGoogle Scholar
- Vigneswaran N, Wu J, Nagaraj N, Adler-Storthz K, Zacharias W: Differential susceptibility of metastatic and primary oral cancer cells to TRAIL-induced apoptosis. Int J Oncol. 2005, 26 (1): 103-112.PubMedGoogle Scholar
- Ho TH, Charlet B, Poulos MG, Singh G, Swanson MS, Cooper TA: Muscleblind proteins regulate alternative splicing. EMBO J. 2004, 23 (15): 3103-3112. 10.1038/sj.emboj.7600300.PubMedPubMed CentralView ArticleGoogle Scholar
- Gadducci A, Baicchi U, Marrai R, Facchini V, del Bravo B, Fosella PV, Fioretti P: Pretreatment plasma levels of fibrinopeptide-A (FPA), D-dimer (DD), and von Willebrand factor (vWF) in patients with operable cervical cancer: influence of surgical-pathological stage, tumor size, histologic type, and lymph node status. Gynecol Oncol. 1993, 49 (3): 354-358. 10.1006/gyno.1993.1139.PubMedView ArticleGoogle Scholar
- Kuopio T, Kankaanranta A, Jalava P, Kronqvist P, Kotkansalo T, Weber E, Collan Y: Cysteine proteinase inhibitor cystatin A in breast cancer. Cancer Res. 1998, 58 (3): 432-436.PubMedGoogle Scholar
- Endoh H, Tomida S, Yatabe Y, Konishi H, Osada H, Tajima K, Kuwano H, Takahashi T, Mitsudomi T: Prognostic model of pulmonary adenocarcinoma by expression profiling of eight genes as determined by quantitative real-time reverse transcriptase polymerase chain reaction. J Clin Oncol. 2004, 22 (5): 811-819. 10.1200/JCO.2004.04.109.PubMedView ArticleGoogle Scholar
- Jager GJ, Barentsz JO, Oosterhof GO, Witjes JA, Ruijs SJ: Pelvic adenopathy in prostatic and urinary bladder carcinoma: MR imaging with a three-dimensional TI-weighted magnetization-prepared-rapid gradient-echo sequence. AJR Am J Roentgenol. 1996, 167 (6): 1503-1507.PubMedView ArticleGoogle Scholar
- Lyng H, Landsverk KS, Kristiansen E, DeAngelis PM, Ree AH, Myklebost O, Hovig E, Stokke T: Response of malignant B lymphocytes to ionizing radiation: gene expression and genotype. Int J Cancer. 2005, 115 (6): 935-942. 10.1002/ijc.20962.PubMedView ArticleGoogle Scholar
- De Angelis PM, Clausen OP, Schjolberg A, Stokke T: Chromosomal gains and losses in primary colorectal carcinomas detected by CGH and their associations with tumour DNA ploidy, genotypes and phenotypes. Br J Cancer. 1999, 80 (3–4): 526-535. 10.1038/sj.bjc.6690388.PubMedView ArticleGoogle Scholar
- Lyng H, Badiee A, Svendsrud DH, Hovig E, Myklebost O, Stokke T: Profound influence of microarray scanner characteristics on gene expression ratios: analysis and procedure for correction. BMC Genomics. 2004, 5 (1): 10-10.1186/1471-2164-5-10.PubMedPubMed CentralView ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.PubMedView ArticleGoogle Scholar
- Ensembl. [http://www.ensembl.org/Homo_sapiens/searchview]
- Galteland E, Sivertsen EA, Svendsrud DH, Smedshammer L, Kresse SH, Meza-Zepeda LA, Myklebost O, Suo Z, Mu D, DeAngelis PM, Stokke T: Translocation t(14;18) and gain of chromosome 18/BCL2: effects on BCL2 expression and apoptosis in B-cell non-Hodgkin's lymphomas. Leukemia. 2005, 19 (12): 2313-2323. 10.1038/sj.leu.2403954.PubMedView ArticleGoogle Scholar
- Sivertsen EA, Galteland E, Mu D, Holte H, Meza-Zepeda L, Myklebost O, Patzke S, Smeland EB, Stokke T: Gain of chromosome 6p is an infrequent cause of increased PIM1 expression in B-cell non-Hodgkin's lymphomas. Leukemia. 2006, 20 (3): 539-542. 10.1038/sj.leu.2404094.PubMedView ArticleGoogle Scholar
- Wang Y, Kringen P, Kristensen GB, Holm R, Baekelandt MM, Olivier M, Skomedal H, Hainaut P, Trope CG, Abeler VM, Nesland JM, Borresen-Dale AL, Helland A: Effect of the codon 72 polymorphism (c.215G>C, p.Arg72Pro) in combination with somatic sequence variants in the TP53 gene on survival in patients with advanced ovarian carcinoma. Hum Mutat. 2004, 24 (1): 21-34. 10.1002/humu.20055.PubMedView ArticleGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.PubMedPubMed CentralView ArticleGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.PubMedPubMed CentralView ArticleGoogle Scholar
- NCBI UniGene database. [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/UniGene]
- SOURCE database. [http://source.stanford.edu]
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.