This article has Open Peer Review reports available.
The influence of tumor size and environment on gene expression in commonly used human tumor lines
© Gieseg et al; licensee BioMed Central Ltd. 2004
Received: 05 March 2004
Accepted: 15 July 2004
Published: 15 July 2004
The expression profiles of solid tumor models in rodents have been only minimally studied despite their extensive use to develop anticancer agents. We have applied RNA expression profiling using Affymetrix U95A GeneChips to address fundamental biological questions about human tumor lines.
To determine whether gene expression changed significantly as a tumor increased in size, we analyzed samples from two human colon carcinoma lines (Colo205 and HCT-116) at three different sizes (200 mg, 500 mg and 1000 mg). To investigate whether gene expression was influenced by the strain of mouse, tumor samples isolated from C.B-17 SCID and Nu/Nu mice were also compared. Finally, the gene expression differences between tissue culture and in vivo samples were investigated by comparing profiles from lines grown in both environments.
Multidimensional scaling and analysis of variance demonstrated that the tumor lines were dramatically different from each other and that gene expression remained constant as the tumors increased in size. Statistical analysis revealed that 63 genes were differentially expressed due to the strain of mouse the tumor was grown in but the function of the encoded proteins did not link to any distinct biological pathways. Hierarchical clustering of tissue culture and xenograft samples demonstrated that for each individual tumor line, the in vivo and in vitro profiles were more similar to each other than any other profile. We identified 36 genes with a pattern of high expression in xenograft samples that encoded proteins involved in extracellular matrix, cell surface receptors and transcription factors. An additional 17 genes were identified with a pattern of high expression in tissue culture samples and encoded proteins involved in cell division, cell cycle and RNA production.
The environment a tumor line is grown in can have a significant effect on gene expression but tumor size has little or no effect for subcutaneously grown solid tumors. Furthermore, an individual tumor line has an RNA expression pattern that clearly defines it from other lines even when grown in different environments. This could be used as a quality control tool for preclinical oncology studies.
Preclinical animal models of human tumors represent a major tool for the selection and development of effective anticancer agents. There is a considerable number of well characterized human cancer cell lines, many of which can be grown as solid tumors (either subcutaneously or orthotopicaly) in immunodeficient mice. The ease of use and low cost of these models make them desirable for screening in vivo activity of anticancer compounds as compared to induced or transgenic rodent tumor models. Rodent models also allow pharmacodynamic and pharmacokinetic parameters to be directly measured and related to antitumor efficacy. Commonly used human cancer cell lines, such as the panel of 60 lines (NCI60) used by the Developmental Therapeutics Program at the National Cancer Institute, have been extensively studied at the molecular level in vitro and to a lesser extent in vivo. Recent advances in genomic technology allow molecular characterization of these models to an extent never before possible using RNA expression profiling, comparative genomic hybridization, proteomic and metabonomic profiling. Ross et al.  looked at the in vitro gene expression in the NCI60 panel. The gene expression pattern for many lines indicated a relationship with the tumor tissue of origin and also correlated with doubling time and drug metabolism. Virtanen et al.  have analyzed the expression patterns of 85 lung tumor samples from both clinical samples and established tumor lines. The fresh tumors clustered according to pathological subtype with many of the cell lines also clustering within the same groups. Pedersen et al.  performed a similar study with human small cell lung cancer cell lines and compared them to ressected tissue samples. Dan et al.  recently performed RNA expression profiling on 39 human cancer cell lines and related the gene expression patterns to chemosensitivity of 55 anticancer agents. Zembutsu et al.  performed a similar study but profiled 85 human cancer xenografts.
What was apparent in all these studies was that each tumor model was distinctly different and could be distinguished from related models using routine data analysis methods. Furthermore, it appears that it is possible to identify models that are no longer true to their origins. Ross et al.  identified MDA-MB-435 as a possible melanoma derived line, which is at odds with its supposed origins of a metastatic breast carcinoma . However, such identification also raises the possibility that the sample that Ross et al. obtained was not the true "MDA-MB-435." Clearly there is the opportunity to use profiling technology as an advanced quality control method that could not only identify mislabeled tumor lines but possibly genetic drift.
While these reports have tremendous value they do not address some basic questions about human tumor models that could impact the design of in vivo drug development studies. For example, most, if not all, xenografts demonstrate Gompertzian growth kinetics with continuously increasing doubling times as they grow larger . Despite this obvious change in biology, it is unknown whether xenografts alter their expression as they increase in size or whether expression differences exist for the same tumor line grown both as a solid tumor in a mouse or in tissue culture. To investigate these questions we grew human colon tumors (HCT-116 and Colo205) in Nu/Nu mice and harvested tumors at three sizes (200 mg, 500 mg, 1000 mg). Isolated RNA was analyzed on Affymetrix GeneChips. The effect of mouse strain on tumor gene expression was also investigated by comparing tumors at 500 mg grown in Nu/Nu and C.B-17 SCID mouse strains. Lastly we profiled additional tumor models to investigate the changes in gene expression that occur when a tumor line is grown in vivo or in vitro.
Xenograft and tissue culture methods
Cell lines were grown in DMEM/F12 supplemented with 10% fetal bovine serum (Invitrogen, Carlsbad, California). They were passaged in 75-cm2 tissue culture flasks in an atmosphere of 5% CO2 in air and were subcultured weekly, using 0.05% trypsin EDTA (Invitrogen, Carlsbad, California). RNA from tissue culture samples was harvested at mid-log phase. Human tumor xenografts were grown in either Nu/Nu mice or C.B-17 SCID female mice obtained from Charles River Laboratories (Wilmington, Massachusetts) between four and five weeks of age. Mice were housed five to a cage in animal rooms maintained at between 21–25°C with a 12 h alternating dark/light cycle. All animal studies were conducted under Veterinary Use Protocols approved by the Institutional Animal Care and Use Committee. Tumors were maintained by serial passage of 30 mg tumor fragments between animals, implanted subcutaneously into the right axillary region using a trocar needle aseptically. Tumors were passed when the primary had reached between 500 and 1000 mg and were never passed more than ten times. Tumor growth was followed by caliper measurements of perpendicular measures of the tumor. The weight in mg was estimated by the formula: tumor weight = a(b2)/2, where a and b are the tumor length and width respectively in mm. Tumor tissue was harvested immediately following animal sacrifice by excising the tumor and powdering it in a liquid nitrogen cooled crucible and pestle. Tumor powder was stored at -80°C until RNA isolation.
RNA was extracted using TRIZOL® reagent (Invitrogen, Carlsbad, California) according to the manufacturer's protocol. RNA integrity was monitored using denaturing agarose gel electrophoresis in 1X MOPS. Biotinylated target RNA was prepared from 15 μg of total RNA using the Affymetrix protocol. Briefly, double-stranded cDNA was prepared from the RNA template using a modified oligo-dT primer containing a 5' T7 RNA polymerase promoter sequence and the Superscript Choice System for cDNA Synthesis (Invitrogen, Carlsbad, California). Following phenol-chloroform extraction and ethanol precipitation, one-half of the cDNA reaction (0.5 – 1.0 μg) was used as the template in an in vitro transcription reaction containing T7 RNA polymerase, a mixture of unlabeled ATP, CTP, GTP, and UTP, and biotin-11-CTP and biotin-16-UTP (BioArray High Yield Kit, ENZO, Farmingdale, New York). The resulting biotinylated-cRNA "target" was purified on an affinity resin (RNeasy, Qiagen, Valencia, California) and quantified using the convention that 1 O.D. 260 nm corresponds to 40 μg/mL of RNA. Typical yields ranged from 50 – 100 μg with transcript sizes between 3.0 to 0.25 kilonucleotides as determined by denaturing gel electrophoresis. Fifteen micrograms of biotinylated cRNA was randomly fragmented to an average size of 50 nucleotides by incubating at 94°C for 35 minutes in 40 mM TRIS-acetate, pH 8.1, 100 mM potassium acetate, and 30 mM magnesium acetate. The fragmented cRNA was hybridized in a solution containing 100 mM MES, 1 M [Na+], 20 mM EDTA, 0.01% TWEEN 20, 50 pM of Control Oligonucleotide B2, 0.1 mg/mL of sonicated herring sperm DNA, and 0.5 mg/mL BSA for 16 hours at 45°C on either the Human U95A or the Mouse U74A Affymetrix GeneChips® (Affymetrix, Santa Clara, California). Each hybridization included a mixture of four bacterial biotinylated-RNA transcripts (BioB, BioC, BioD, and cre) spiked at 1.5, 5, 25, and 100 pM, respectively. The hybridization reactions were processed and scanned according to the standard Affymetrix protocols.
All arrays were global scaled to a target intensity value of 600 using the standard Affymetrix protocol. Calculation of the scaling factor, background, noise and percent present, was performed according to Affymetrix protocols using the Data Mining Tool (Affymetrix Santa Clara, California). All resulting data sets were filtered using the absolute call metric (present or absent) using Microsoft Access (Microsoft Corporation, Redmond, Washington). Genes selected had expression levels classified as present at least once in the samples selected for the particular analysis.
To determine the relationship between tumor samples harvested at different sizes the filtered RNA profiling data was analyzed with classic multidimensional scaling (MDS), implemented in R [8, 9]. MDS is an unsupervised learning technique that attempts to preserve the relationship between points from high dimensional space at lower dimensional spaces. The program R is a free integrated suite of software facilities for data manipulation, calculation and graphical display http://www.r-project.org.
Analysis of variance (ANOVA) was used to test the effects of tumor size and tumor line, using the following model: Expression of gene i ~ tumor.size + tumor.line + tumor.line*tumor.size. To test the effect due to mouse background, we used the following model: Expression of gene i ~ mouse background. RNA profiling data filtered on the absolute call metric was used for this analysis and ANOVA was implemented in R. ANOVA is a statistical linear modeling procedure that partitions the total variance into parts corresponding to various sources in the model [10, 11]. It has been used previously in microarray data analysis [12–15].
To rigorously select genes with expression differences between samples, ANOVA p-values were adjusted using multiple comparison procedures. Multiple comparison procedures are tools to adjust p-values that might be inflated as a result of performing multiple hypothesis tests. The Benjamini and Hochberg procedure controls the false discovery rate, which is the expected fraction of false discoveries in all rejected hypothesis . This procedure is less stringent than methods controlling the family wise error rate (e.g. the Bonferroni correction); hence it is more powerful.
Hierarchical clustering was performed in GeneSpring 5.0 (Silicon Genetics, Redwood City, California). The distinction calculation from Spotfire DecisionSite 6.2 (Spotfire Inc. Somerville, Massachusetts) was used to select genes differentially expressed in xenograft samples or tissue culture samples. All data from the tissue culture samples that had an in vivo pair (8 samples) were selected into one group and all data from the xenograft samples (8 samples) were selected into a second group. Genes were prefiltered using the absolute call metric by selecting genes that were present at least once in the selected samples. A distinction value score and p-value was calculated for each gene. The score (≥1) and p-value (≤0.001) was then used to select genes that were differentially expressed between xenograft samples and tissue culture samples. To functionally classify gene lists, web resources such as NCBI (National Center for Biotechnology Information, http://www.ncbi.nlm.nih.gov/) were searched and the data compiled. Further searching for gene associations in PubMed (NCBI) was also performed.
Variation in tumor xenograft gene expression due to size
We focused on two human colon carcinoma xenografts (HCT-116 and Colo205) to investigate the effects of tumor size and mouse strain on gene expression. Samples were harvested in quadruplicate at three different tumor sizes (200 mg, 500 mg and 1000 mg) for both tumor models grown in Nu/Nu mice (except for the 500 mg sample of Colo205 where five samples were harvested). These sizes were selected as they represent the range at which sensitivity to anticancer agents are traditionally tested and because most models approximate log-linear growth at these sizes. RNA expression profiling data was obtained from Affymetrix U95A GeneChips containing approximately 12600 genes. Genes present (above background) once or more across all samples were selected for further analysis (approximately 7600 genes).
To rigorously identify genes that may suggest functionally significant changes as a tumor increases in size from 200 mg to 1000 mg the ANOVA p-values were adjusted using multiple comparison procedures. Following this analysis none of the previously identified 154 genes had p-values < 0.05. This indicated that there was no change in gene expression as a human tumor xenograft increased in size from 200 mg to 1000 mg. This was true for both the HCT-116 and the Colo205 tumor lines.
Variation in tumor xenograft gene expression due to mouse strain
Genes with high expression in HCT-116 tumors grown in C.B-17 SCID mice
Affymetrix Probe Set ID
Sequence Probe Derived From
Translation/RNA binding/RNA splicing
splicing factor 3a, subunit 1, 120 kDa
dyskeratosis congenita 1, dyskerin
eukaryotic translation elongation factor 1 alpha 2
splicing factor 3a, subunit 2, 66 kDa
squamous cell carcinoma antigen recognised by T cells 3
DEAD/H box polypeptide 21
hepatocyte growth factor-regulated tyrosine kinase substrate
casein kinase 1, epsilon
hepatoma-derived growth factor
ras-related C3 botulinum toxin substrate
Small Molecule transport
solute carrier family 7, member 5
calcium channel, voltage-dependent, beta 1 subunit
solute carrier family, member 2
Protein Complex Assembly
nucleosome assembly protein 1-like 4
Protein transport and folding
chromosome 12 open reading frame 8
Unknown or other functions
signal recognition particle 72 kDa
hypothetical protein FLJ20719
KIAA0193 gene product
bromodomain containing 2
bromodomain containing 2
coagulation factor VIII-associated
DNA2 DNA replication helicase 2-like
BCL2/adenovirus E1B 19 kDa interacting protein 3-like
lysosomal-associated membrane protein 1
arsenate resistance protein ARS2
syntrophin, beta 2
dual-specificity tyrosine-(Y)-phosphorylation regulated kinase 2
likely ortholog of mouse nervous system polycomb 1
MCM5 minichromosome maintenance deficient 5, cell division cycle 46
Genes with high expression in HCT-116 tumors grown in Nu/Nu mice
Affymetrix Probe Set ID
Sequence Probe Derived From
Translation/RNA binding/RNA splicing
ribosomal protein L32
ribosomal protein L15
ribosomal protein L18a
ribosomal protein L5
ribosomal protein S15a
chromosome 20 open reading frame 14
DEAD/H box polypeptide 18 (Myc-regulated)
eukaryotic translation elongation factor 1 alpha 1
Rho guanine nucleotide exchange factor (GEF) 12
src family associated phosphoprotein 2
protein phosphatase 2, regulatory subunit B (B56), gamma isoform
GDP dissociation inhibitor 2
Unknown or Other Functions
a disintegrin and metalloproteinase domain 17
acetyl-Coenzyme A carboxylase beta
annexin A2 pseudogene 1
zinc finger protein 259
cytochrome P450, 51
heat shock 70 kDa protein 8
Homo sapiens mRNA sequence, IMAGE clone 417820
BTG family, member 3
B-cell receptor-associated protein BAP29
prothymosin, alpha (gene sequence 28)
serine (or cysteine) proteinase inhibitor, clade H, member 2
antigen identified by monoclonal antibody Ki-67
ubiquitin carrier protein
zinc finger protein 91
Comparison of lines grown in both tissue culture and solid tumor xenografts
Human tumor lines used in this study
Tissue of Origin
Tissue culture sample
Number of genes with fold change values greater than or equal to 2 with tissue culture sample as the baseline
Number of genes with ≥2 fold increase
Number of genes with ≥2 fold decrease
Total with ≥2 fold change
Another striking feature revealed by clustering analysis was that none of the tumor lines were particularly similar to each other. The tumor lines did not necessarily cluster based upon the tissue of origin from which each line was thought to be derived. However, a group of breast carcinoma lines that did cluster together (BT474, ZR-75-1, MDA-MB-453, MCF-7, SKBr-3 and MDA-MB-468), although MDA-MB-435 and MDA-MB-231 did not cluster with this group. The two lung carcinoma lines (H2009 and H125) did cluster closely, but MDA-MB-231 and HCT-116 (colon carcinoma) also clustered with this group. MDA-MB-435 and Colo205 clustered independently from all groups.
Selection and functional assessment of differential expressed genes in xenograft or tissue culture samples
Genes with high expression in xenograft samples
Affymetrix Probe Set ID
Sequence Probe Set Derived From
Extracellular Matrix/Cell Adhesion
Beta-galactoside-binding lectin, mRNA sequence
collagen, type I, alpha 2
armadillo repeat gene deletes in velocardiofacial syndrome
matrix metalloproteinase 11 (stromelysin 3)
major histocompatibility complex, class II, DO beta
O-linked N-acetylglucosamine (GlcNAc) transferase
Cell surface receptor
corticotropin releasing hormone receptor 2
glutamate receptor, metabotropic 4
5-hydroxytryptamine (serotonin) receptor 1B
CD3E antigen, epsilon polypeptide (TiT3 complex)
tumor necrosis factor receptor superfamily, member 25
zinc finger protein
pleiomorphic adenoma gene-like 1
TAF1 RNA polymerase II, TATA box binding protein-associated factor
thyroid hormone receptor interactor 8
adaptor-related protein complex 2, alpha 2 subunit
vesicle-associated membrane protein 1 (synaptobrevin 1)
Translation/RNA binding/RNA splicing
RNA binding motif protein 5
Unknown or other functions
hypothetical protein CG003
putative ankyrin-repeat containing protein
Yes-associated protein 1, 65 kDa
par-3 partitioning defective 3 homolog
SOCS box-containing WD protein SWiP-1
echinoderm microtubule associated protein like 2
serine (or cysteine) proteinase inhibitor, clade F, member 1
CDC-like kinase 1
iduronate 2-sulfatase (Hunter syndrome)
pyruvate dehydrogenase kinase, isoenzyme 1
Genes with high expression in tissue culture samples
Affymetrix ProbeSet ID
Sequence ProbeSet Derived From
serine/threonine kinase 6
tubulin, gamma 1
kinesin-like 5 (mitotic kinesin-like protein 1)
CDC45 cell division cycle 45-like (S. cerevisiae)
CSE1 chromosome segregation 1-like (yeast)
Translation/RNA binding/RNA splicing
eukaryotic translation initiation factor 3, subunit 9 eta, 116 kDa
small nuclear ribonucleoprotein polypeptide E
RNA binding motif protein 14
Cockayne syndrome 1 (classical)
polymerase (RNA) II (DNA directed) polypeptide E, 25 kDa
Unknown or other functions
immature colon carcinoma transcript 1
glutaryl-Coenzyme A dehydrogenase
ATP synthase, subunit f, isoform 2
thymidine kinase 1, soluble
The first part of our study was designed to address whether gene expression patterns in xenografts varied with tumor size. No significant changes (p-value < 0.05) were found in any genes after statistical analysis of the RNA profiling data. This was a somewhat surprising result as there is evidence to suggest that tumors show size-dependent biological variation in both a clinical setting and as implanted solid tumors in rodents. It is generally assumed clinically that the larger the tumor mass the lower the likelihood of curing the patient, regardless of the treatment [17, 18]. Pathological and genetic heterogeneity increase with increasing tumor cell number and contribute to the emergence of clinical drug resistance . As tumors grow larger their vascular surface decreases, intercapillary distance increases, interstial pressure increases and necrotic foci develop . Tumor doubling times and cell loss also increase with increasing tumor size .
Implanted animal tumor models most likely represent clinical end-stage disease and may not necessarily recapitulate clinical behavior seen during tumor development in humans. There are few reports on the molecular and biochemical changes that occur as implanted solid tumors increase in size. Massaad et al.  found that some drug metabolizing enzyme systems, including glutathione-S-transferases (GST) did change both activity and expression with increasing tumor size in the mouse colon adenocarcinoma Co38.
Preclinical tumor models show variability in response to treatment with a given anti-tumor agent despite being derived from a single tumor and being implanted into inbred mice . This has been attributed to heterogeneity of the host or the tumor cell populations [23, 24] along with tumor-to-tumor variations in perfusion and drug distribution . We found there was no significant difference in gene expression between tumors of either the same size or different size as assessed in four replicate animals. If differences in tumor population, host metabolism or perfusion exist between animals implanted with the same tumors, our data suggests they are not detectable using gene expression analysis of the whole tumor.
Our results suggest that once tumor xenografts have grown to 200 mg, physiological influences on tumor transcription (such as the interaction with host stroma, nutrient supply and oxygenation) have reached steady state and do not change appreciably as the tumor grows to 1000 mg. It is also conceivable that tumors smaller and larger than the sizes investigated in this study may show gene expression changes. When a tumor is first implanted the hypoxic environment is likely to drive the expression of pro-angiogenic cytokines that will recruit new vessels to the tumor. As a tumor becomes very large (>2 g) angiogenic factors may again become highly expressed due to the inability for the existing vasculature to adequately support the large tumor mass. We have not excluded the possibility that there may be changes in protein activities or levels, or that other tumor lines may show changes as the size increases. But since the phenomenon was independently observed in two tumor lines it suggests it may be common to many other lines. The results also imply that tumors harvested within this size range without drug treatment should be very comparable. It remains possible that selective isolation and analysis of tumor regions or subpopulations by microdissection may identify expression differences.
By comparing the same tumor line grown in different mouse strains we identified 63 genes with a significant change in gene expression with biological functions ranging from cellular signaling, RNA processing and translation. Although different genes were differentially expressed in different mouse strains, particular gene functions did not associate with one strain. Therefore, the biological significance of these changes remains unclear. It should be noted that although the changes were statistically significant, in general the magnitude of the change was small. Only 12% of the significant genes displayed a change greater or equal to 2-fold. This raises the question of whether statistical significance equates to biological significance and how relevant thresholds can be determined for analyzing expression profiling data.
In the second part of the study we compared the gene expression patterns of tumor lines grown as xenografts with the same lines grown in tissue culture. We found that each xenograft – tissue culture pair was more similar to one another than any other line. This has been previously observed with RNA expression profiling of small cell lung cancer lines grown as both xenografts and tissue culture . In our data, the one exception to this clustering pattern, ZR-75-1, exhibited aberrant growth characteristics; it had a rapid doubling time (2–3 days) and was not estrogen dependent as has been reported [26, 27].
There was a gene expression pattern consistent with growth either in a xenograft or tissue culture environment, but remarkably it consisted of only a small number of genes. We identified 36 genes with high relative expression in xenograft samples and 17 genes with high relative expression in tissue culture samples. Their biological functions were consistent with the differences between a tissue culture sample and a xenograft sample. Genes expressed at increased levels in tissue culture samples were primarily involved in cell division, cell cycle, transcription and translation and were consistent with a greater percentage of cycling cells in the tissue culture samples, which we have previously observed with flow cytometry (data not shown). Many of the genes expressed in xenograft samples encoded for proteins involved in extracellular matrix, cell adhesion, cell surface receptors transcription and translation and suggest increased interactions with an in vivo stromal environment, including the development of a 3-dimensional matrix.
There are several reports in the literature of tumor line expression profiling studies demonstrating that lines derived from a common tissue of origin group together in most cases [1, 4, 5]. We observed that most of the breast tumor lines clustered together but the other lines generally did not show clustering patterns based upon tissue of origin. However, our study evaluated fewer samples and consisted predominantly of breast cancer tumor lines with only two representatives of non-small cell lung cancer and colon cancer.
Since each tumor line was substantially different to every other line tesed it suggests that it would be possible to identify a group of genes that could be used to distinguish individual lines. Also, it has been our experience that it is possible to identify aberrant samples resulting from labeling error by expression profiling that were missed by other methods. What is less clear is whether genetic drift causes changes in RNA expression that can be identified using RNA profiling technology. Recently a report has been published showing that MCF-7 sublines were substantially different to each other both at the genetic and RNA expression levels  suggesting that variation within individual lines can be identified. Given this, it seems reasonable that RNA expression profiling could be used as a comprehensive methodology for identifying aberrant or incorrectly labeled samples. This would provide an additional quality control tool to standardize tumor models and the in vivo testing of therapeutic agents.
Our data suggest that the environment a tumor line is grown in can have a significant effect on gene expression even though tumor size has little or no effect for a subcutaneously grown solid tumor. Furthermore, an individual tumor line has an RNA expression pattern that clearly defines it from other lines even when grown in different environments such as tissue culture or in vivo. Routine RNA expression profiling of a selected set of genes could be used as a quality control tool for preclinical oncology studies and would allow the standardization of tumor models used worldwide.
We kindly acknowledge Dr Sarah Fowler for assistance with proofing this manuscript and Dr William Elliott for helpful discussion about the growth characteristics of solid tumor xenografts.
- Ross DT, Scherf U, Eisen MB, Perou CM, Rees C, Spellman P, Iyer V, Jeffrey SS, Van de RM, Waltham M, Pergamenschikov A, Lee JC, Lashkari D, Shalon D, Myers TG, Weinstein JN, Botstein D, Brown PO: Systematic variation in gene expression patterns in human cancer cell lines. Nat Genet. 2000, 24: 227-235. 10.1038/73432.View ArticlePubMedGoogle Scholar
- Virtanen C, Ishikawa Y, Honjoh D, Kimura M, Shimane M, Miyoshi T, Nomura H, Jones MH: Integrated classification of lung tumors and cell lines by expression profiling. Proc Natl Acad Sci USA. 2002, 99: 12357-12362. 10.1073/pnas.192240599.View ArticlePubMedPubMed CentralGoogle Scholar
- Pedersen N, Mortensen S, Sorensen SB, Pedersen MW, Rieneck K, Bovin LF, Poulsen HS: Transcriptional gene expression profiling of small cell lung cancer cells. Cancer Res. 2003, 63: 1943-1953.PubMedGoogle Scholar
- Dan S, Tsunoda T, Kitahara O, Yanagawa R, Zembutsu H, Katagiri T, Yamazaki K, Nakamura Y, Yamori T: An integrated database of chemosensitivity to 55 anticancer drugs and gene expression profiles of 39 human cancer cell lines. Cancer Res. 2002, 62: 1139-1147.PubMedGoogle Scholar
- Zembutsu H, Ohnishi Y, Tsunoda T, Furukawa Y, Katagiri T, Ueyama Y, Tamaoki N, Nomura T, Kitahara O, Yanagawa R, Hirata K, Nakamura Y: Genome-wide cDNA microarray screening to correlate gene expression profiles with sensitivity of 85 human cancer xenografts to anticancer drugs. Cancer Res. 2002, 62: 518-527.PubMedGoogle Scholar
- Cailleau R, Olive M, Cruciger QV: Long-term human breast carcinoma cell lines of metastatic origin: preliminary characterization. In Vitro. 1978, 14: 911-915.View ArticlePubMedGoogle Scholar
- Jackson RC: From The Growth Kinetics of Large Tumors. In The Theoretical Foundations of Cancer Chemotherapy Introduced by Computer Models. Edited by: Jackson RC. 1992, San Diego, California: Academic Press Inc, 259-298.Google Scholar
- Ihaka R, Gentleman RR: A Language for Data Analysis and Graphics. Journal of Computational and Graphical Statistics. 1996, 5: 299-314.Google Scholar
- Venables WN, Ripley BD: Modern Applied Statistics with S-PLUS. 1997, Springer, SecondView ArticleGoogle Scholar
- Fisher RD: Statistical methods for research workers. 1950, London: Oliver and BoydGoogle Scholar
- Fisher LD, van Belle G: Biostatistics: a methodology for the health sciences. 1993, New York: John Wiley & Sons IncGoogle Scholar
- Kerr MK, Martin M, Churchill GA: Analysis of variance for gene expression microarray data. J Comput Biol. 2000, 7: 819-837. 10.1089/10665270050514954.View ArticlePubMedGoogle Scholar
- Coombes KR, Highsmith WE, Krogmann TA, Baggerly KA, Stivers DN, Abruzzo LV: Identifying and quantifying sources of variation in microarray data using high-density cDNA membrane arrays. J Comput Biol. 2002, 9: 655-669. 10.1089/106652702760277372.View ArticlePubMedGoogle Scholar
- Pritchard CC, Hsu L, Delrow J, Nelson PS: Project normal: defining normal variance in mouse gene expression. Proc Natl Acad Sci USA. 2001, 98: 13266-13271. 10.1073/pnas.221465998.View ArticlePubMedPubMed CentralGoogle Scholar
- Pavlidis P, Noble WS: Analysis of strain and regional variation in gene expression in mouse brain. Genome Biol. 2001, 2: 1-15. 10.1186/gb-2001-2-10-research0042.View ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate; a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B. 1995, 57: 289-300.Google Scholar
- DeVita VT: The James Ewing lecture. The relationship between tumor mass and resistance to chemotherapy. Implications for surgical adjuvant treatment of cancer. Cancer. 1983, 51: 1209-1220.View ArticlePubMedGoogle Scholar
- Michaelson JS, Silverstein M, Wyatt J, Weber G, Moore R, Halpern E, Kopans DB, Hughes K: Predicting the survival of patients with breast carcinoma using tumor size. Cancer. 2002, 95: 713-723. 10.1002/cncr.10742.View ArticlePubMedGoogle Scholar
- Baylin SB: From Clonal selection and heterogeneity of human solid neoplasms. In Design of Models for Testing Cancer Therapeutic Agents. Edited by: Fidler IJ, White RJ. 1981, New York: Van Nostrand Reinhold, 50-63.Google Scholar
- Jain RK: Delivery of novel therapeutic agents in tumors: physiological barriers and strategies. J Natl Cancer Inst. 1989, 81: 570-576.View ArticlePubMedGoogle Scholar
- Massaad L, Chabot GG, Toussaint C, Koscielny S, Morizet J, Bissery MC, Gouyette A: Influence of tumor size on the main drug-metabolizing enzyme systems in mouse colon adenocarcinoma Co38. Cancer Chemother Pharmacol. 1994, 34: 497-502. 10.1007/s002800050179.View ArticlePubMedGoogle Scholar
- Schabel FM, Griswold DP, Corbett TH, Laster WR, Lloyd HH, Rose WC: From Variable response of advanced solid tumors of mice to treatment with anticancer drugs. In Design of Models for Testing Cancer Therapeutic Agents. Edited by: Fidler IJ, White RJ. 1981, New York: Van Nostrand Reinhold, 95-113.Google Scholar
- Nelson JA, Hokanson JA, Jenkins VK: Role of the host in the variable chemotherapeutic response of advanced Ridgway osteogenic sarcoma. Cancer Chemother Pharmacol. 1982, 9: 148-155.View ArticlePubMedGoogle Scholar
- Trope C: From Different susceptibilities of tumor cell subpopulations to cytotoxic agents. In Design of Models for Testing Cancer Therapeutic Agents. Edited by: Fidler IJ, White RJ. 1981, New York: Van Nostrand Reinhold, 64-79.Google Scholar
- Simpson-Herren L, Noker PE, Wagoner SD: Variability of tumor response to chemotherapy. II. Contribution of tumor heterogeneity. Cancer Chemother Pharmacol. 1988, 22: 131-136.View ArticlePubMedGoogle Scholar
- Gutman M, Couillard S, Roy J, Labrie F, Candas B, Labrie C: Comparison of the effects of EM-652 (SCH57068), tamoxifen, toremifene, droloxifene, idoxifene, GW-5638 and raloxifene on the growth of human ZR-75-1 breast tumors in nude mice. Int J Cancer. 2002, 99: 273-278. 10.1002/ijc.10302.View ArticlePubMedGoogle Scholar
- Cameron DA, Ritchie AA, Langdon S, Anderson TJ, Miller WR: Tamoxifen induced apoptosis in ZR-75-1 breast cancer xenografts antedates tumor regression. Breast Cancer Res Treat. 1997, 45: 99-107. 10.1023/A:1005850827825.View ArticlePubMedGoogle Scholar
- Nugoli M, Chuchana P, Vendrell J, Orsetti B, Ursule L, Nguyen C, Birnbaum D, Douzery EJ, Cohen P, Theillet C: Genetic variability in MCF-7 sublines: evidence of rapid genomic and RNA expression profile modifications. BMC Cancer. 2003, 3: 13-10.1186/1471-2407-3-13.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/4/35/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.