- Research article
- Open access
- Published:
Assessing breast cancer cell lines as tumour models by comparison of mRNA expression profiles
Breast Cancer Research volume 17, Article number: 114 (2015)
Abstract
Introduction
Breast cancer researchers use cell lines to model myriad phenomena ranging from DNA repair to cancer stem cell phenotypes. Though appropriate, and even requisite, for many studies, the suitability of cell lines as tumour models has come into question owing to possibilities of tissue culture artefacts and clonal selection. These issues are compounded by the inability of cancer cells grown in isolation to fully model the in situ tumour environment, which also contains a plethora of non-tumour cell types. It is thus important to understand similarities and differences between cancer cell lines and the tumours that they represent so that the optimal tumour models can be chosen to answer specific research questions.
Methods
In the present study, we compared the RNA-sequencing transcriptomes of a collection of breast cancer cell lines to transcriptomes obtained from hundreds of tumours using The Cancer Genome Atlas. Tumour purity was accounted for by analysis of stromal and immune scores using the ESTIMATE algorithm so that differences likely resulting from non-tumour cells could be accounted for.
Results
We found the transcriptional characteristics of breast cancer cell lines to mirror those of the tumours. We identified basal and luminal cell lines that are most transcriptionally similar to their respective breast tumours. Our comparison of expression profiles revealed pronounced differences between breast cancer cell lines and tumours, which could largely be attributed to the absence of stromal and immune components in cell culture. A focus on the Wnt pathway revealed the transcriptional downregulation or absence of several secreted Wnt antagonists in culture. Gene set enrichment analysis suggests that cancer cell lines have enhanced proliferation and glycolysis independent of stromal and immune contributions compared with breast cancer cells in situ.
Conclusions
This study demonstrates that many of the differences between breast cancer cell lines and tumours are due to the absence of stromal and immune components in vitro. Hence, extra precautions should be taken when modelling extracellular proteins in vitro. The specific differences discovered emphasize the importance of choosing an appropriate model for each research question.
Introduction
Since the establishment of the HeLa cell line in 1951, cell lines have been an integral part of cancer research, and their use has tremendously advanced understanding of molecular cancer biology [1]. However, the suitability of these models has come into question, as many in vitro phenomena are challenging to replicate in vivo. Interpreting the potential clinical significance of discoveries made using cell lines requires an understanding of the extent to which these cell lines represent in vivo tumours.
Since the first breast cancer cell line, BT-20, was established in 1958 [2], various other immortalized primary tumour cell lines have been established at exceptionally poor efficiencies [3, 4]. This low efficiency has often been attributed to slow growth rates of tumour cells in culture as compared with associated stromal cells, such as fibroblasts [5]. To overcome this issue, most established breast cancer lines have been derived from pleural effusions, which provide an abundance of dissociated, aggressive tumour cells with very few contaminating cell types. The pattern of growth of these tumour cells is characterized by a slow initial proliferation, followed by exponential expansion of a few cells, suggestive of clonal selection for cells that are particularly proliferative and amenable to culture [6–8].
Another caveat of cell culture is the loss of the in vivo microenvironment (changes summarized in [9]). During the derivation process, tumour cells are removed from a very complex, partially hypoxic three-dimensional microenvironment; maintained in nutrient media supplemented with a surplus of growth factors, including glucose; and passaged indefinitely at relatively high atmospheric oxygen levels. In such a drastically altered microenvironment, it would not be surprising if cell lines differed substantially from the tumours they were established to represent.
Genomic and transcriptional differences between cancer cell lines and tumour samples have been investigated in several studies [10–13]. For example, in gliomas, it was shown that expression profiles of tumour cell primary cultures were much closer to profiles obtained from clinically resected tumours than to profiles of immortalized cancer cell lines [14]. In breast cancer, clustering based on expression profiles has elucidated the many clinically relevant subtypes in cell lines and tumours (summarized in [15]) [16–20]. However, modern RNA-sequencing (RNA-seq) data have not yet been used to directly compare the expression profiles of breast cancer cell lines with breast tumours. As well, in vitro signatures are the combined effect of adaptation to cell culture and selection for specific cellular subtypes. Dissecting out the influence of either of these two phenomena has remained a substantial obstacle in any cell line–tumour transcriptional comparison.
Recent transcriptional profiling of a collection of breast cancer cell lines [21] and hundreds of tumours from The Cancer Genome Atlas (TCGA) [19] has enabled a direct mRNA comparison of cell lines and tumours. In this study, we focus on RNA-seq transcriptional profiles mined from TCGA and the Gene Expression Omnibus (GEO) series [GEO:GSE48213] and investigate the strengths and weaknesses of cell lines as in vitro breast cancer models. In addition, we seek to identify the breast cancer cell lines that are most transcriptionally representative of their respective tumour subtype. Importantly, we are able to correlate most of the highly differentially expressed genes to tumour stromal or immune signatures, highlighting the importance of considering the entire niche in cancer modelling. Finally, we summarize relevant breast cancer cell line genomic alterations. In our study, we used RNA-seq data to broaden the dynamic range of transcript detection and extend earlier efforts by including more cell lines and by considering and quantifying stromal and immune cell contributions to help elucidate the origin of detected differences.
Methods
Datasets
Level 3 TCGA RNAseqV2 gene expression data were obtained from the TCGA Data Portal [22] in August 2014. RNA-seq expression data were retrieved in September 2014 for 50 luminal and basal breast cancer cell lines profiled in the GEO database [GEO:GSE48213] [21]. Oestrogen receptor (ER) status and subtype data were available in the original publication for cell lines and were accessed via the UCSC Cancer Genomics Browser (RNASeqV2 defined [23]) for tumours in September 2014. Breast cancer cell line copy number information and mutation data for 1651 genes were retrieved from the cBio Cancer Genomics Portal [24] for the Cancer Cell Line Encyclopaedia (CCLE) in March 2015.
Data preparation
Relative abundance (in transcripts per million [TPM]) was calculated for 975 breast tumours by multiplying the scaled estimate data by 106, and for 50 breast cell lines by converting fragments per kilobase of exons per million mapped reads to TPM. To avoid infinite values in log calculations, a value of 1 was added to all TPM values before log2 transformation. Values for the genes that were available in both datasets (16,282 coding genes in total) were combined for further analysis.
Gene expression profiling analysis
The top 5000 genes by variance across the combined dataset were chosen for principal component analysis as well as for hierarchical clustering using 1 − c (where c is Pearson’s correlation coefficient) as the distance and Ward’s agglomeration method (ward.D2). The 5000 most variable genes were also used to compute the Pearson’s correlation coefficient of all the cell line–tumour pairs in a subtype-specific manner. The cell lines were ranked based on their average correlation with all tumours of their respective subtype. Significant differences in relative transcript abundances between cell lines and tumours were calculated with Welch’s t test, and p values were corrected for multiple testing using the Benjamini-Hochberg method. Enrichment for functionally related genes between the two datasets was tested using Generally Applicable Gene-set Enrichment (GAGE v2.12.3; Bioconductor: [25]) with Kyoto Encyclopedia of Genes and Genomes [26] gene sets with fewer than 200 items.
Tumour purity
Stromal and immune scores were defined for tumours by ESTIMATE scores (Estimation of STromal and Immune cells in MAlignant Tumour tissues using Expression data) using RNASeqV2 data as previously described [27], and accessed in October 2014 via the bioinformatics portal at the Department of Bioinformatics and Computational Biology, University of Texas MD Anderson Cancer Center [28]. Pearson’s correlation coefficient was used to calculate the association of specific genes with stromal and immune signatures. To decrease hits of transcripts more likely due to tumour purity issues, transcripts that correlated with stromal or immune scores (|r|>0.2) were filtered from the list of differentially expressed genes, and a new ranked list was generated. Stromal and immune correlations were calculated for each gene set by averaging the stromal and immune Pearson’s correlation coefficients of the essential genes (as determined using the GAGE package).
Cell line genomic analysis
The CCLE [29] investigators examined the mutational status of 1651 genes by hybrid capture sequencing and genome-wide copy number analysis. In our genomics summary, we considered all breast cancer cell lines that were available in both the CCLE and at GEO accession number [GEO:GSE48213]. The fraction of the genome altered represents the fraction of the genome that has a log2 copy number value above 0.2 or below −0.2. Selected mutational events were considered if they were found to be significantly altered from or in associated healthy tissue in the original TCGA study [19]. Copy number status was investigated for significantly mutated genes in the TCGA study that also displayed frequent copy number amplifications or deletions.
PubMed citation analysis
The number of PubMed abstracts that mentioned 1 of the 50 breast cancer cell lines was determined as an estimator of frequency of use in laboratories. Hits were determined using the PubMed search function [30] on 18 March 2015. Several punctuation alternatives were used for the cell line names. For the cell lines LY2 and MX1, searches were conducted with the term ‘cells’ to eliminate the inclusion of abstracts that mentioned the LY2 and MX1 genes.
Statistical analysis
We conducted all analyses and visualizations in the RStudio programming environment (v0.98.501; [31]). The R/Bioconductor packages ggplot2, plyr, gplots, ggdendro and GAGE were used as appropriate.
Results
Comparison of cell lines and tumour expression profiles
To evaluate the transcriptional fidelity of breast cancer cell lines to tumours, we compared the mean expression values of 16,282 coding genes in oestrogen receptor–positive (ER+) and oestrogen receptor–negative (ER−) cell lines to ER+ and ER− tumours. The mean expression values of cell lines and tumours were similar, though the mean expression values of ER− cell lines and tumours are more highly correlated (r = 0.90) than ER+ cell lines and tumours (r = 0.88) (Fig. 1a, ER+; Fig. 1b, ER−). However, closer inspection revealed a point of interest: Almost all outliers were genes with high expression in tumours and low expression in cell lines.
We further explored the relationship of cell lines and tumours by conducting principal component analysis (Fig. 1c) and found four clusters clearly divided based on sample group (cell line or tumour; principal component 1 [PC1]) and ER status (principal component 2 [PC2]). One of the main differences between cell lines and tumours is the absence of certain cellular components (e.g., stromal and immune cells). Given that many of the outliers in Fig. 1a, b were genes that had higher expression in tumours than in cell lines, and PC1 in Fig. 1c was largely responsible for the distance between tumours and cell lines, we correlated PC1 with stromal and immune scores in tumours as determined by using the ESTIMATE paradigm [27]. We found that stromal scores strongly positively correlated (r = 0.74) with PC1 (Fig. 1d). Thus, the loss of the stromal component likely has significant repercussions in vitro.
In the principal component analysis, we observed that ER− cell lines clustered closer to their respective tumours than ER+ cell lines, indicating again that ER− cell lines may be more representative of their tumour counterparts than ER+ cell lines. Expression-based, unsupervised hierarchal clustering also revealed this trend. Although cell lines cluster apart from all tumours, they cluster closer to the largely ER−/basal subtype division of tumours than to the largely ER+/luminal subtype division (Additional file 1: Figure S1).
Top differentially expressed genes are genes correlated with stromal and immune scores
We found that the top 1 % of genes differentially expressed in cell culture were all genes that had lower or undetectable expression in culture compared with tumours (Fig. 2). To determine the contribution of stromal and immune cellular compartments to this observation, we correlated the expression of all genes with stromal and immune scores in tumours. We found that 134 of the top 163 differentially expressed genes were highly correlated with stromal or immune scores in tumours (r > 0.5). Representative correlations are shown in Fig. 2c, e, g and Additional file 2.
Recently, Winslow et al. investigated different breast cancer cellular compartments using laser capture microdissection [32]. When we examined the genes that they found to be upregulated in breast cancer stromal cells compared with malignant cells, we found that 99 % of them were significantly downregulated in breast cancer cell lines and that their average correlation with the stromal score was 0.65 (data not shown). Taken together, this supports the theory that the downregulation of many genes in cell culture is likely a consequence of losing stromal and immune cellular compartments.
Wnt antagonists are underrepresented in cell culture
Because previous studies have suggested that Wnt pathway components are provided by stromal cells in certain situations [33, 34], we investigated the expression of various Wnt pathway members in the datasets. We determined that numerous putative Wnt antagonist transcripts were underrepresented in cell culture and highly correlated with stromal scores (Fig. 3a). However, apart from WNT2, all other Wnt agonists (Fig. 3b) and receptors (Fig. 3c) did not display this pattern. This provides evidence that the stromal compartment of tumours provides a unique and non-redundant role in tumours that cannot be modelled accurately using cancer cell monoculture.
Top differentially expressed genes not correlated with stromal and immune scores
It is expected that cell line in vitro signatures are a combined result of selection for the malignant subtype of cells and in vitro adaptation. Though the expression level changes connected to stromal and immune contributions seem to be the most pronounced, we were interested in investigating subtler distinctions that are more likely a result of in vitro adaptations. In an attempt to overcome the contributions from stromal and immune cell lineages, we identified a new top 1 % of differentially expressed genes after removing the genes whose expression correlated with stromal or immune scores (|r| > 0.2) (Fig. 4 and Additional file 3). This list is more likely to reflect changes in cancer cells induced by the cell-culturing process.
Gene set enrichment analysis reveals the enrichment of proliferative gene sets in cell culture
Gene set enrichment analysis revealed 41 upregulated and 35 downregulated gene sets in cell lines compared with tumours (Fig. 5). Cell lines were enriched for gene sets associated with proliferation and metabolism, whereas gene sets associated with extracellular matrix interactions were underrepresented, akin to a pattern that was previously observed when tumours and cell lines were compared [35]. Stromal and immune correlations were calculated for each pathway’s essential genes and averaged for the entire pathway to provide a measure of the contribution of stromal and immune components to resultant gene set perturbation (Table 1). We determined that downregulated gene sets were quite strongly correlated with stromal and immune scores. Specifically, genes that allow tumour cells to interact with extracellular components, such as syndecan 2 and CD36, were lower in cell lines than in tumours. However, upregulated gene sets were not as strongly correlated with stromal and immune scores, indicating that their perturbation is likely a result of the differences in the tumour cells themselves.
Ranking of cell lines by transcriptional similarity to their tumour counterparts
To assess the transcriptional suitability of individual cell lines as tumour models, we calculated the correlation coefficients of the top 5000 most variable genes in all subtype-specific cell line–tumour pairs and ranked the cell lines based on their average correlation coefficient (Table 2, Fig. 6). Although this evaluation is not fully comprehensive of all potential genomic and epigenomic differences, it does provide a reasonable guide for choosing cell lines that are most transcriptionally representative of their respective tumour subtype. Ranking the breast cancer cell lines based on correlation leads to a spread of the cell lines from most representative (Table 2, top) to least representative (Table 2, bottom). Keeping with the trend previously observed, the highest ranked basal cell line (HCC70; r = 0.58) was more strongly correlated with respective tumours than the highest ranked luminal cell lines (BT483; r = 0.52). It is also reassuring to note that two of the most extensively published luminal cell lines, T47D and MCF7, are ranked fourth and fifth, respectively, of the 27 luminal lines that were evaluated. However, the top ranked luminal and basal cell lines (luminal: BT483, ZR7530 and 600MPE; basal: HCC70, MX1 and HCC3153) are infrequently used as breast cancer models and account for only 0.4 % of publications on this cell line panel.
We went further and investigated mutation status and copy number alterations of breast cancer cell lines profiled by the CCLE. With our transcriptional correlation ranking, we created a summary of all of these events (Fig. 6), hoping that it could help inform breast cancer cell line choice. The frequency of these somatic mutational events in cell lines mirrors the frequency found in tumours, with a few notable differences: TP53, PTEN, NF1 and PTPRD were mutated at significantly higher frequencies in cell lines as compared with tumours (p < 0.05 by binomial test) (Additional file 4).
Discussion
This study is the first transcriptional comparison of cancer cell lines and tumours to methodologically account for the contributions of tumour stromal and immune cellular components. We demonstrate, for the first time to our knowledge, using RNA-seq data, that breast cancer cell lines generally represent breast tumours, with notable exceptions. First, many extracellular proteins thought to be lost in breast cancers may actually be supplied in situ by the stroma. Second, many genes associated with proliferation and metabolism are highly expressed in culture. Hence, whereas certain aspects of breast cancer biology can be studied using breast cancer cell lines alone, others (in particular those involving factors in the extracellular space) should include additional relevant cell types.
This study revealed that, in general, basal/ER− cell lines were more representative of their respective tumours than luminal/ER+ cell lines. In addition, 60 % of cell lines in this study were ER−, as compared with only 23 % of the primary tumours (p < 0.0001 by two-tailed Fisher’s exact test), an overrepresentation of the ER− status in cell lines, which has been observed previously [1]. The reason for this discrepancy remains unknown. However, it may be due to the fact that most cell lines were obtained from metastatic tumours and pleural effusions and thus represent the most aggressive variants that could be adapted to culture (a trend previously reported in renal cancer [36]). We would expect this phenomenon to be especially pronounced for the ER+/luminal subtype, which is characteristically a less aggressive subtype of breast cancer. Additionally, as cells are grown in culture, the epithelial phenotype is lost in favour of more mesenchymal traits, a type of in vitro epithelial–mesenchymal transition which would result in greater transcriptional distance between the more epithelial ER+/luminal cell lines and the respective tumours [1].
Despite the transcriptional differences between cell lines and tumours, we were nonetheless interested in determining the most transcriptionally representative breast cancer cell lines. In our analysis, we found that the correlation coefficients of individual breast cancer cell lines versus tumours varied from 0.41 to 0.58. This was remarkably similar to the range of 0.43–0.60 that was observed in an analysis of ovarian cell lines and tumours [11]. Interestingly, the top correlation of any individual cell line could be exceeded by a fictional cell line composed of the averages of all cell line gene expression values (luminal, 0.52 for BT483 vs. 0.62; basal, 0.58 for HCC70 vs. 0.60) (data not shown). This points to the importance of including multiple cell lines in any analysis to ensure that any observed phenomenon is not a product of a single outlier.
A fundamental limitation of cell culture models is that the environment created by culture conditions is markedly different from the breast cancer microenvironment [9]. The loss of stromal and immune cells in culture is one major drawback of monoculture models. Emerging evidence supports the notion that tumour stromal cells play exceptionally important roles in tumour initiation, progression and metastasis [37–39]. In fact, studies have shown that depletion of fibroblast activation protein–expressing stromal cells leads to suppression of primary tumour growth and metastasis [40]. Our research indicates that loss of the stromal and immune components is the principal transcriptional difference between cell lines and tumours. It also suggests that the stroma has a unique and significant role that often is not accounted for in in vitro studies. For example, several studies have looked at the expression levels and functional roles of various Wnt antagonists (e.g., secreted frizzled-related proteins [SFRPs]) in cell culture, and researchers have drawn conclusions about their absence and mechanisms of action in this context [41–44]. However, given that we found the expression levels of various SFRPs to be high in tumours and strongly correlated with stromal scores, we should recognize that looking at these proteins in tumour cell monoculture may not be appropriate. In fact, given their roles as matricellular proteins, it would not be surprising if their effects in vivo are quite different than those observed in vitro.
In broader investigations using gene set enrichment analysis, we observed an enrichment in cell line proliferative and metabolic gene sets, similar to those reported in other studies [45–47]. The upregulation of these gene sets could be due to two phenomena: (1) malignant cellular adaptation/selection or (2) genes more highly expressed in the malignant cells are upregulated in cell lines as a result of the enrichment of this cell subtype in culture. If the latter is true, we would expect a negative correlation with stromal/tumour purity score. For one of the gene sets, DNA replication, we observed such a negative correlation with stromal score (r = −0.27). Thus, the expansion of malignant cells in cell culture likely plays a role in the upregulation of this gene set. However, none of the other upregulated proliferative/metabolic gene sets display this correlation. This suggests, on the one hand, that either the derivation process or the continuous culturing of cell lines selects for a highly proliferative subset of cells. On the other hand, many of the underrepresented gene sets were matrix- or immune-related and tightly correlated with stromal or immune scores, once again indicating that loss of the stromal and immune compartments has pronounced consequences in transcriptional programs observed in cell culture.
Conclusions
Important efforts are being made to systematically compare tumours and cell lines using DNA mutation, copy number and gene expression data from a diverse spectrum of tumour types [5, 9–11, 13, 35]. In this study, we focused on breast cancer expression data and sought to identify major transcriptional differences between cell lines and tumours while accounting for variation resulting from stromal and immune components. We determined that basal cell lines are transcriptionally better models of their respective tumours than luminal cell lines. We ranked cell lines based on their transcriptional similarity to tumour samples and recommend that cell line choices be informed by this summary. We have also pointed out situations where cell line monoculture may not be the best tumour model. Fortunately, there exist many other tumour models (e.g., patient-derived xenografts, co-cultures and three-dimensional systems) that may more appropriately represent these situations. Knowing in which contexts cell lines have high or low fidelity to tumours can help direct tumour model choice, optimizing the clinical relevance of future research efforts.
Abbreviations
- CCLE:
-
Cancer Cell Line Encyclopaedia
- ER:
-
Oestrogen receptor
- ESTIMATE:
-
Estimation of STromal and Immune cells in MAlignant Tumours using Expression data
- GAGE:
-
Generally applicable gene set enrichment
- GEO:
-
Gene Expression Omnibus
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- PC1:
-
Principal component 1
- PC2:
-
Principal component 2
- RNA-seq:
-
RNA sequencing
- SFRP:
-
Secreted frizzled-related protein
- TCGA:
-
The Cancer Genome Atlas
- TPM:
-
transcripts per million
References
Lacroix M, Leclercq G. Relevance of breast cancer cell lines as models for breast tumours: an update. Breast Cancer Res Treat. 2004;83:249–89.
Lasfargues EY, Ozzello L. Cultivation of human breast carcinomas. J Natl Cancer Inst. 1958;21:1131–47.
Gazdar AF, Kurvari V, Virmani A, Gollahon L, Sakaguchi M, Westerfield M, et al. Characterization of paired tumor and non-tumor cell lines established from patients with breast cancer. Int J Cancer. 1998;78:766–74.
Amadori D, Bertoni L, Flamigni A, Savini S, De Giovanni C, Casanova S, et al. Establishment and characterization of a new cell line from primary human breast carcinoma. Breast Cancer Res Treat. 1993;28:251–60.
Burdall SE, Hanby AM, Lansdown MR, Speirs V. Breast cancer cell lines: friend or foe? Breast Cancer Res. 2003;5:89–95.
Giard DJ, Aaronson SA, Todaro GJ, Arnstein P, Kersey JH, Dosik H, et al. In vitro cultivation of human tumors: establishment of cell lines derived from a series of solid tumors. J Natl Cancer Inst. 1973;51:1417–23.
Engel LW, Young NA, Tralka TS, Lippman ME, O’Brien SJ, Joyce MJ. Establishment and characterization of three new continuous cell lines derived from human breast carcinomas. Cancer Res. 1978;38:3352–64.
Cailleau R, Young R, Olivé M, Reeves WJ. Breast tumor cell lines from pleural effusions. J Natl Cancer Inst. 1974;53:661–74.
van Staveren WCG, Solís DYW, Hébrant A, Detours V, Dumont JE, Maenhaut C. Human cancer cell lines: experimental models for cancer cells in situ? For cancer stem cells? Biochim Biophys Acta. 2009;1795:92–103.
Neve RM, Chin K, Fridlyand J, Yeh J, Baehner FL, Fevr T, et al. A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes. Cancer Cell. 2006;10:515–27.
Domcke S, Sinha R, Levine DA, Sander C, Schultz N. Evaluating cell lines as tumour models by comparison of genomic profiles. Nat Commun. 2013;4:2126.
Wistuba II, Behrens C, Milchgrub S, Syed S, Ahmadian M, Virmani AK, et al. Comparison of features of human breast cancer cell lines and their corresponding tumors. Clin Cancer Res. 1998;4:2931–8.
Dairkee SH, Ji Y, Ben Y, Moore DH, Meng Z, Jeffrey SS. A molecular “signature” of primary breast cancer cultures; patterns resembling tumor tissue. BMC Genomics. 2004;5:47.
Lee J, Kotliarova S, Kotliarov Y, Li A, Su Q, Donin NM, et al. Tumor stem cells derived from glioblastomas cultured in bFGF and EGF more closely mirror the phenotype and genotype of primary tumors than do serum-cultured cell lines. Cancer Cell. 2006;9:391–403.
Holliday DL, Speirs V. Choosing the right cell line for breast cancer research. Breast Cancer Res. 2011;13:215.
Perou CM, Sørlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, et al. Molecular portraits of human breast tumours. Nature. 2000;406:747–52.
Ross DT, Scherf U, Eisen MB, Perou CM, Rees C, Spellman P, et al. Systematic variation in gene expression patterns in human cancer cell lines. Nat Genet. 2000;24:227–35.
Charafe-Jauffret E, Ginestier C, Monville F, Finetti P, Adélaïde J, Cervera N, et al. Gene expression profiling of breast cell lines identifies potential new basal markers. Oncogene. 2006;25:2273–84.
Cancer Genome Atlas Network. Comprehensive molecular portraits of human breast tumours. Nature. 2012;490:61–70.
Kao J, Salari K, Bocanegra M, Choi YL, Girard L, Gandhi J, et al. Molecular profiling of breast cancer cell lines defines relevant tumor models and provides a resource for cancer gene discovery. PLoS One. 2009;4:e6146.
Daemen A, Griffith OL, Heiser LM, Wang NJ, Enache OM, Sanborn Z, et al. Modeling precision treatment of breast cancer. Genome Biol. 2013;14:R110. A published erratum appears in. Genome Biol. 2015;16:95.
TCGA Data Portal. https://tcga-data.nci.nih.gov/tcga/.
USCS Cancer Genomics Browser. https://genome-cancer.ucsc.edu/.
cBioPortal for Cancer Genomics. http://www.cbioportal.org/. Accessed 17 July 2015.
GAGE v2.12.3; Bioconductor. http://www.bioconductor.org/packages/release/bioc/html/gage.html.
Kyoto Encyclopedia of Genes and Genomes. http://www.genome.jp/kegg/.
Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.
Bioinformatics portal at the Department of Bioinformatics and Computational Biology, University of Texas MD Anderson Cancer Center. http://bioinformatics.mdanderson.org/main/Main_Page.
Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, et al. The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012;483:603–7.
National Library of Medicine, National Center for Biotechnology Information. PubMed. http://www.ncbi.nlm.nih.gov/pubmed/. Accessed 17 July 2015.
RStudio programming environment. https://www.rstudio.com/.
Winslow S, Leandersson K, Edsjö A, Larsson C. Prognostic stromal gene signatures in breast cancer. Breast Cancer Res. 2015;17:23.
Kabiri Z, Greicius G, Madan B, Biechele S, Zhong Z, Zaribafzadeh H, et al. Stroma provides an intestinal stem cell niche in the absence of epithelial Wnts. Development. 2014;141:2206–15.
Pongracz J, Hare K, Harman B, Anderson G, Jenkinson EJ. Thymic epithelial cells provide WNT signals to developing thymocytes. Eur J Immunol. 2003;33:1949–56.
Ertel A, Verghese A, Byers SW, Ochs M, Tozeren A. Pathway-specific differences between tumor cell lines and normal and tumor tissue cells. Mol Cancer. 2006;5:55.
Ebert T, Bander NH, Finstad CL, Ramsawak RD, Old LJ. Establishment and characterization of human renal cancer and normal kidney cell lines. Cancer Res. 1990;50:5531–6.
Mao Y, Keller ET, Garfield DH, Shen K, Wang J. Stromal cells in tumor microenvironment and breast cancer. Cancer Metastasis Rev. 2013;32:303–15.
Östman A, Augsten M. Cancer-associated fibroblasts and tumor growth – bystanders turning into key players. Curr Opin Genet Dev. 2009;19:67–73.
Pietras K, Östman A. Hallmarks of cancer: interactions with the tumor stroma. Exp Cell Res. 2010;316:1324–31.
Loeffler M, Krüger JA, Niethammer AG, Reisfeld RA. Targeting tumor-associated fibroblasts improves cancer chemotherapy by increasing intratumoral drug uptake. J Clin Invest. 2006;116:1955–62.
Suzuki H, Toyota M, Caraway H, Gabrielson E, Ohmura T, Fujikane T, et al. Frequent epigenetic inactivation of Wnt antagonist genes in breast cancer. Br J Cancer. 2008;98:1147–56.
Suzuki H, Gabrielson E, Chen W, Anbazhagan R, van Engeland M, Weijenberg MP, et al. A genomic screen for genes upregulated by demethylation and histone deacetylase inhibition in human colorectal cancer. Nat Genet. 2002;31:141–9.
Nojima M, Suzuki H, Toyota M, Watanabe Y, Maruyama R, Sasaki S, et al. Frequent epigenetic inactivation of SFRP genes and constitutive activation of Wnt signaling in gastric cancer. Oncogene. 2007;26:4699–713.
Lee J, Yoon YS, Chung JH. Epigenetic silencing of the WNT antagonist DICKKOPF-1 in cervical cancer cell lines. Gynecol Oncol. 2008;109:270–4.
Perou CM, Jeffrey SS, Van de Rijn M, Rees CA, Eisen MB, Ross DT, et al. Distinctive gene expression patterns in human mammary epithelial cells and breast cancers. Proc Natl Acad Sci U S A. 1999;96:9212–7.
Birgersdotter A, Sandberg R, Ernberg I. Gene expression perturbation in vitro—a growing case for three-dimensional (3D) culture systems. Semin Cancer Biol. 2005;15:405–12.
Sandberg R, Ernberg I. The molecular portrait of in vitro growth by meta-analysis of gene-expression profiles. Genome Biol. 2005;6:R65.
Acknowledgements
This work was supported by an Alberta Innovates – Health Solutions Translational Health Chair in cancer and a Canadian Breast Cancer Foundation operating grant (to LMP). LMP was the recipient of the Peter-Lougheed Premier New Investigator Award from the Canadian Institutes of Health Research. SDF is a recipient of an Ontario Graduate Scholarship, and KMV is a Vanier Scholar.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare they have no competing interests.
Authors’ contributions
KMV, SDF and LMP conceived and designed the analysis. KMV prepared the data and performed all data analyses. All authors participated in interpreting the results. KMV wrote the manuscript. SDF and LMP participated in revising the manuscript. All authors read and approved the final manuscript.
Additional files
Additional file 1: Figure S1.
Unsupervised hierarchal clustering of 50 breast cancer cell lines and 1025 TCGA breast cancer tumour samples shows that cell lines cluster apart from tumour samples. Hierarchical clustering on the 5000 most variable genes was performed using 1 − c (where c is Pearson’s correlation coefficient) as the distance and Ward’s agglomeration method. Though cell lines cluster apart from tumours, basal cell lines cluster closer to their respective tumours than luminal cell lines do. Cell line clustering followed previously observed subtype divisions. (PDF 175 kb)
Additional file 2:
Top 1 % of differentially expressed genes between tumours and cell lines (with adjusted p values by t test). (XLSX 26 kb)
Additional file 3:
Top 1 % of differentially expressed genes between tumours and cell lines after filtering out genes correlated with stromal and immune scores in tumours (absolute tumour stromal and immune scores, Pearson correlation coefficient <0.2, adjusted p values by t test). (XLSX 29 kb)
Additional file 4:
Comparison of the frequency of select mutational events in breast cancer tumours versus cell lines. Differences in mutational frequencies of selected genes in TCGA BRCA tumours (n = 503) and CCLE breast cancer cell lines (n = 31). (XLSX 11 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Vincent, K.M., Findlay, S.D. & Postovit, L.M. Assessing breast cancer cell lines as tumour models by comparison of mRNA expression profiles. Breast Cancer Res 17, 114 (2015). https://doi.org/10.1186/s13058-015-0613-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13058-015-0613-0