| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Biology, Pathobiology and Genetics |
1 Department of Mathematics and Computer Science, Eindhoven University of Technology, Eindhoven, the Netherlands and Departments of 2 Medical Oncology, 3 Clinical Chemistry, and 4 Pathology, VU University Medical Center, Amsterdam, the Netherlands
Requests for reprints: Bauke Ylstra, Department of Pathology-Microarray Core Facility, VU University Medical Center, Medical Faculty, Room B356, Van der Boechorststraat 7-9, NL-1081 BT Amsterdam, the Netherlands. Phone: 31-204448299; Fax: 31-204448318; E-mail: b.ylstra{at}vumc.nl.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
| Materials and Methods |
|---|
|
|
|---|
Expression microarrays. A total of four microarray expression experiments were done (Table 1). Arrays were done essentially according to Bergman et al. (4). For the expression arrays, the mouse oligo library version 1.0 (Compugen/Sigma-Aldrich Chemie B.V., Zwijndrecht, the Netherlands) was spotted in duplicate using the SA 72 (Perkin-Elmer, Zaventem, Belgium) with Telechem SMP3 pins (TeleChem International, Inc., Sunnyvale, CA), containing 7,524 oligonucleotides of 65 bp in length representing 7,230 separate genes. RNA isolations and labeling, as well as scanning and image analysis, were done as previously described (4). All data were of high quality: Background signals were very low compared with the foreground signals. Hence, neither spot filtering nor background subtractions were done. Oligo nucleotide sequences (accessible at http://www.ensembl.org/) were verified for the 13 identified genes by Blast onto Mus musculus build 33.1. All 13 sequences matched their designated gene ID. However, RPL13 only had homology with the distal 32 bases of the oligonucleotide sequence. No other significant similarities were found with the first 28 bases; thus, no artifact can be expected from this maldesigned oligonucleotide sequence. Furthermore, the sequence for the spotted oligonucleotide representing HBB-B1 (AF071431) was AAACCCCCTTTCCTGATTTTGCCTGTGAACAATGGTTAATTGTTCCCAAGAGAGCATCTGTCAGT, whereas the actual chromosomal sequence at 7E1 for HBB-B1 has 2 bp differences in the regions printed in bold, which should have been GCTCT according to Ensembl. No significant homologies with other mouse genes were observed.
|
The Bayesian hierarchical model implemented is made of several levels that mimic some of the crucial steps in the microarray experiment, such as hybridization and slide preparation. Errors are propagated through the model, resulting in realistic log ratio estimates. One of the consequences is that the more inaccurately measured spots are down-weighted in the analysis (see Supplementary data). The parametric model provides an accurate compromise between the k-fold rule and t test. The k-fold rule is a biological rather than statistical criterion that preferentially selects for most differentially expressed genes. The standard t test is limited inherently considering the small number of arrays and focuses on statistical significance. This compromise is the posterior probability of a gene being k-fold expressed. Such a probability depends on both the intensity of the mean differential gene expression and the amount of agreement between measurements of the genes on the three arrays of the solid tumors. Alongside, the 2-fold adjusted expression ratios (AER2) across the three arrays are computed, which is the amount by which the log 2 ratio exceeds 1 [= log 2(2)] divided by the SD.
Real-time LightCycler PCR. Seven genes identified by the microarray experiments were verified by real-time PCR with a LightCycler 1.0. (Roche Diagnostics, Almere, the Netherlands) according to the protocols described by Bergman et al. (4). Primers for all seven murine genes RNA were based on the sequence of the gene (http://www.ncbi.nlm.nih.gov/entrez/) designed with the program Primer3 (http://frodo.wi.mit.edu/); OLFR683 forward primer: 5'-GATCAAAGCAGAGGGAGCTG, reverse primer: 5'-AAGGTTCCGTATTCCCTGCT; TRIM21 forward: 5'-CCATGGTGGAGCCTATGAGT, reverse: 5'-GGTGAAGCTTCTCTCCATGC; POLM forward: 5'-CCCGAGTCAACTCAGCTTTC, reverse: 5'-CTGCACAACACCTCACTGCT; COX4I1 forward: 5'-TTCTACTTCGGTGTGCCTTC, reverse: 5'-GCGAAGCTCTCGTTAAACTG; RPL13 forward: 5'-TACTGAAGCCCCACTTCCAC, reverse: 5'-CGGACCTTGGTGTGGTATCT; KRT2-8 forward: 5'-ATCGAGATCACCACCTACCG and reverse: 5'-TGAAGCCAGGGCTAGTGAGT. Primer sequences for RRM1 and the reference gene ß-actin were described by Bergman et al. (4). Expression levels were quantified on five different tumor passages relative to ß-actin (4). Table 1 shows the median log 2 expression ratios of the gene in 26G relative to the expression in 26A with SD, thereby eliminating ß-actin in the equation, such that it is directly comparable with the M:log 2 ratio G/A of the arrays.
Oligo comparative genomic hybridization microarray. DNA from tumor and normal liver samples were isolated using the Wizard Genomic DNA purification kit according to the protocol of the manufacturer (Promega Benelux BV, Leiden, the Netherlands). Labeling and hybridization procedures for the oligo array CGH were done as previously described (8) with the following modifications: mouse oligo library version 2.0 (Compugen/Sigma-Aldrich Chemie B.V., Zwijndrecht, the Netherlands) containing 21,997 oligonucleotides (65 bp) representing 21,587 exon regions, mapped to the University of California Santa Cruz Mouse genome 2003 freeze, was spotted on the arrays. Prehybridization was omitted and Cot-1 concentrations during the hybridization were reduced to 100 µg. Hybridizations were done using a Hybstation12 (Perkin-Elmer). CGH arrays were scanned using a laser scanner (Agilent Technologies, Amstelveen, the Netherlands) and analyzed using Bluefuse Software v.2.0 (Bluegnome Ltd., Cambridge, United Kingdom). Images show fused values; values with confidence higher then 0.35 with the overall Cy3 and Cy5 channels normalized to a log 2 ratio of 0. DNA was isolated from tumors 26A and 26G. The log 2 ratio for tumor 26G was slightly noisier than for tumor 26A, which is generally a result of the DNA quality. Therefore, DNA from a second tumor 26G was isolated and hybridized. The results of this second hybridization were of the same quality as the first 26G hybridization and completely confirmed the results. Because bacterial artificial chromosome CGH arrays are often spotted in triplicate (9), we present the Bluefuse confidence-based average of values of the two spots on the separate arrays for 26G. Original data files for all three arrays were uploaded in MIAME format for expression arrays at GEO (http://www.ncbi.nlm.nih.gov/geo/; 26G accession nos. GSM44665 and GSM44666, and 26A accession no. GSM44664). For interpretation and visualization purposes, smoothing was done by version 2 of aCGH smooth (10), with
set to 2.0.
| Results |
|---|
|
|
|---|
|
|
|
Three of 13 genes with differential expression in resistant colon tumors are involved in DNA replication. To identify a common functional denominator or pathway, we collected data for all 13 genes, such as gene description, gene ontology, oligo sequence, and chromosomal location (Table 2).
One of the highly up-regulated genes identified by our statistical approach, RRM1, was verified by both real-time PCR (Table 2) and Western blot analysis (4). For RRM1, DNTT, POLM, and KRT2-8, our statistical approach not only showed high log 2 ratios and high posterior probabilities of >2-fold differential expression, but also high intensity levels. These intensity levels indicate high RNA expression levels especially because the hybridization efficiency of the oligo probes are largely normalized in their design. According to their gene ontology, 3 of 13 genes are involved in DNA replication, which includes RRM1 in addition to POLM and DNTT (Table 2).
Five of 13 genes identified by expression array analysis are located in a 3Mb region on chromosome 7E1. Five of the 13 genes were found to reside in a 3 Mb region of chromosome 7E1. Moreover, two genes that are both down-regulated are located on a 2.5 Mb region of chromosome 8. The close proximity of the genes and their respective up- or down-regulation suggest chromosomal copy number changes as a consequence of the gemcitabine treatment. To test this hypothesis, we did oligo array CGH. Array CGH detected high-level amplifications at 7E1 spanning a 3.1 Mb chromosomal region (Fig. 2). This high-level amplification encompasses all five genes (Fig. 2C). No additional genes in this region were detected with our expression array analysis. The reduced RNA expression of the two genes at 8E1, COX4I1, and RPL13 coincide with a decrease in copy number of the entire chromosome 8 relative to the sensitive tumor 26A. In addition to the anticipated chromosomal changes, a decrease in copy number of chromosome 19 is observed in tumor 26G relative to 26A.
| Discussion |
|---|
|
|
|---|
7,000 genes were identified by expression microarray analyses that are highly regulated in gemcitabine-resistant versus nonresistant tumors using a Bayesian parametric model. Seven of the 13 genes were verified by LC-PCR, confirming our expression array platform and data analysis. The identification of 8 of these 13 genes could straightforwardly be verified and explained by literature, gene ontology, or chromosomal copy number changes as measured by oligo array CGH. The gene RRM1 was previously identified to be involved in drug resistance by several independent research groups in both cell lines and solid human tumors (2, 3, 11). Gemcitabine is a fluorinated deoxycytidine analogue and thereby interferes with DNA replication (1). Our expression array analysis identified DNTT and POLM as highly up-regulated, highly expressed, with high posterior probabilities and the same gene ontology as RRM1:DNA replication. This suggests that they could act in the same pathway to convey resistance to the tumors. Four additional genes, HBB-BH1, HBB-B1, OLFR683, and TRIM21, are located at the same chromosomal location, 7E1. HBB-BH1, OLFR683, and TRIM21 are also up-regulated along with RRM1. By oligo array CGH, a high and specific amplification of 7E1, encompassing all five genes, was identified. Gene function or ontology of OLFR683 or the hemoglobin genes does not give a logical explanation for their identification by our expression array analysis. The up-regulation of OLFR683, HBB-BH1, and TRIM21, however, can be explained by chromosomal location rather than by gene function.
The particular region on 7E1 has been studied in great detail because it does not only contain an olfactor receptor cluster, but also the ß-globin gene cluster (12). ß-globin gene expression is highly specific and molecularly one of the best-characterized regions of the mouse and human genomes (13). The genes are arranged along the chromosome in order of their expression during development (14), such that the expression of embryonic and adult ß-globin are mutually exclusive. HBB-BH1 is an embryonic ß-globin and is highly up-regulated in the resistant tumor probably as a consequence of the high-level amplification. The expression of HBB-B1, which is located downstream of HBB-BH1, and one of the adult ß-globins is down-regulated, despite its high-level amplification. In addition, one of the adult
-globin genes at chromosome 11 is also down-regulated in the resistant tumor. Thus, it is tempting to speculate that, as a consequence of the high up-regulation of HBB-BH1, both HBB-B1 and HBA-A1 RNA expression is down-regulated. Thus, HBB-BH1 is then a consequence of the amplification along with RRM1 and the down-regulation of HBB-B1 and HBA-A1 a consequence of a consequence.
Interestingly, the 7E1 region corresponds to the imprinted 11p15.5 region in humans (15). Loss or gain of imprinting is a feature of many tumors where changes in dosage compensation leads to changes in the expression levels of the genes involved with subsequent selective growth advantage of these tumor cells. To our knowledge, the relation between gemcitabine and 11p15 has only been described for nonsmall cell lung cancer (2). In the latter case, deletions are assumed to be involved. Our data show that a similar relation exists in colon tumor cells but involves amplification leading to overexpression of RRM1, HBB-B1, OLFR683, and TRIM21. Down-regulation of two other genes, COX4I1 and RPL13, could also be explained by our array CGH analysis because down-regulation of these genes coincides with a deletion of chromosome 8. The genes identified by expression array analysis thus biologically validate that both our oligo array expression platform as well our parametric approach are robust.
Chromosomal instability is generally studied by DNA replication inhibitors such as Aphidicolin (16). It is well established that chromosomal instability is directly associated with DNA repair or delayed DNA replication (17). Here, we show in a solid tumor that, indeed, the nucleoside analogue gemcitabine is capable to convey drug resistance to tumors by chromosomal instability. The role of nucleotide availability in DNA replication and chromosomal stability has recently been shown by Debatisse et al. (18). Acquired drug resistance leading to chromosomal copy number was also observed in cell lines for other nucleoside analogues (19). Our studies, therefore, make it conceivable that this group of nucleoside analogues can convey chromosomal instability as a direct consequence of their interference in DNA replication.
| Acknowledgments |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
We thank Paul van den Ijssel, Paul P. Eijk, and Marianne Tijssen (VU University Medical Center Microarray Facility) for helpful assistance; Joris Vermeesch (Leuven, Belgium) for critically reading the manuscript before submission and Nir Gamliel (Compugen USA, Inc.) and Tony Cox (Sanger Institute, Cambridge, United Kingdom) for making the oligo sequences publicly available and mapping them onto the human genome in Ensembl.
| Footnotes |
|---|
Received 3/ 4/05. Revised 9/ 1/05. Accepted 9/13/05.
| References |
|---|
|
|
|---|
, ß,
,
. J Clin Invest 2005;115:14703.[CrossRef][Medline]
- and ß-globin loci. Genes Dev 1999;13:1124.This article has been cited by other articles:
![]() |
P. Tooker, W.-C. Yen, S.-C. Ng, A. Negro-Vilar, and T. W. Hermann Bexarotene (LGD1069, Targretin), a Selective Retinoid X Receptor Agonist, Prevents and Reverses Gemcitabine Resistance in NSCLC Cells by Modulating Gene Amplification Cancer Res., May 1, 2007; 67(9): 4425 - 4433. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Ylstra, P. van den IJssel, B. Carvalho, R. H. Brakenhoff, and G. A. Meijer BAC to the future! or oligonucleotides: a perspective for micro array comparative genomic hybridization (array CGH) Nucleic Acids Res., January 26, 2006; 34(2): 445 - 450. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. van den IJssel, M. Tijssen, S.-F. Chin, P. Eijk, B. Carvalho, E. Hopmans, H. Holstege, D. K. Bangarusamy, J. Jonkers, G. A. Meijer, et al. Human and mouse oligonucleotide-based array CGH Nucleic Acids Res., December 16, 2005; 33(22): e192 - e192. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |