Thyroid cancer cell lines: an overview

Human thyroid cancer cell lines are the most used models for thyroid cancer studies. They must be used with detailed knowledge of their characteristics. These in vitro cell lines originate from differentiated and dedifferentiated in vivo human thyroid tumors. However, it has been shown that mRNA expression profiles of these cell lines were closer to dedifferentiated in vivo thyroid tumors (anaplastic thyroid carcinoma, ATC) than to differentiated ones. Here an overview of the knowledge of these models was made. The mutational status of six human thyroid cancer cell lines (WRO, FTC133, BCPAP, TPC1, K1, and 8505C) was in line with previously reported findings for 10 genes frequently mutated in thyroid cancer. However, the presence of a BRAF mutation (T1799A: V600E) in WRO questions the use of this cell line as a model for follicular thyroid carcinoma (FTC). Next, to investigate the biological meaning of the modulated mRNAs in these cells, a pathway analysis on previously obtained mRNA profiles was performed on five cell lines. In five cell lines, the MHC class II pathway was down-regulated and in four of them, ribosome biosynthesis and translation pathways were up-regulated. mRNA expression profiles of the cell lines were also compared to those of the different types of thyroid cancers. Three datasets originating from different microarray platforms and derived from distinct laboratories were used. This meta-analysis showed a significant higher correlation between the profiles of the thyroid cancer cell lines and ATC, than to differentiated thyroid tumors (i.e., PTC or FTC) specifically for DNA replication. This already observed higher correlation was obtained here with an increased number of in vivo tumors and using different platforms. In summary, this would suggest that some papillary thyroid carcinoma or follicular thyroid carcinoma (PTC or FTC) cell lines (i.e., TPC-1) might have partially lost their original DNA synthesis/replication regulation mechanisms during their in vitro cell adaptation/evolution.

Introduction

The experimental study of human cancers uses in vitro and in vivo models. Among the various possible experimental models, human cancer cell lines are frequently used (van Staveren et al., 2009). They have retained hallmarks of cancer cells; they are pure, genetically identical, easily propagated and can be genetically manipulated. A cell line originates from a tissue and is obtained by selection of the most rapidly proliferating and resistant cells in monolayer during passages.

Results obtained on a cell line are sometimes directly extrapolated for in vivo cancers which produced this cell line (Yeung et al., 2007; Wang et al., 2008). However, the representativity of the cancer cell line may be distorted by a cross contamination of one cell line by another (Ribeiro et al., 2008; Schweppe et al., 2008), an in vitro evolution of the cell line (van Staveren et al., 2009), a strong genomic instability due to the number of passages or a risk of infection (Harlin and Gajewski, 2008). Thus, a systematic verification of the status of these cell lines is important (van Staveren et al., 2007; Ribeiro et al., 2008; Schweppe et al., 2008).

Thyroid cancer is the most frequent endocrine cancer (Kondo et al., 2006; Sipos and Mazzaferri, 2010). There are various types of thyroid carcinomas, the carcinomas from thyrocyte are largely the most frequent: papillary thyroid carcinoma (PTC), follicular thyroid carcinoma (FTC) anaplastic thyroid carcinoma (ATC) and an intermediate form between PTC/FTC and ATC, the poorly differentiated thyroid cancer. Each type is characterized by a set of mutations leading to increased cellular proliferation and dedifferentiation (Catalano et al., 2010).

PTC is the most frequent type of human thyroid carcinoma (Kondo et al., 2006). The genetic alterations most often found in PTC are BRAF point mutations, accounting for 40–60% of the cases, and RET/PTC rearrangements which are present in about 20% of the cases. The most frequent BRAF mutation occurs in the serine/threonine kinase domain (V600E) and leads to the constitutive kinase activity of the protein (Xing, 2010). This mutation could be a key mutation for the treatment and the diagnosis of the most aggressive PTC (Nucera et al., 2010; Xing, 2010). RET/PTC rearrangements are the result of a fusion between the 3′ end of a receptor of the tyrosine kinase family (RET) and the 5′ end of a gene constitutively expressed in thyrocytes. The most frequent rearrangements of this type are RET/PTC1 and RET/PTC3. The fusion results in the constitutive activation of the truncated tyrosine kinase portion of RET by autophosphorylation due to the dimerization domain of the heterologous gene (Catalano et al., 2010).

The genetic changes most often found in FTC are RAS point mutations (approximately 45% of the cases) and PAX8/PPARγ rearrangements (approximately 35% of the cases). Mutations of the RAS gene activate the mitogenic MAPK and the PI3K pathways. PAX8/PPARγ rearrangements are the result of a fusion between the 5′ end of the PAX8 transcription factor and the 3′ end of PPARγ, a member of the nuclear hormone receptor superfamily, constitutively expressed in thyrocytes. Different rearrangements between the 2 genes have been described, and sometimes reported within a single tumor (Lacroix et al., 2004). PAX8/PPARγ chimeric proteins have a dominant negative effect on the wild-type PPARγ. In addition to these, alterations in the tumor suppressor gene PTEN (10% of the cases) and the PI3KCA oncogenes (10% of the cases) have been described in this tumor type.

ATC is the most dedifferentiated, aggressive thyroid cancer (Kondo et al., 2006). The genetic alterations most often observed are the ones already described in PTC and FTC with the exception of RET/PTC and PAX8/PPARγ rearrangements, but in addition also include mutations in TP53 and β-catenin genes.

For each type of human thyroid carcinoma, derived cell lines have been generated. In this study, we used six different human cancer cell lines derived from different types of thyroid cancers. The WRO (Estour et al., 1989) and FTC133 (Goretzki et al., 1990) cell lines are derived from FTC; the BCPAP (Fabien et al., 1994) commonly known as a PTC cell line is derived from a poorly differentiated PTC; TPC1 (Tanaka et al., 1987) and K1 (Challeton et al., 1997) cell lines are derived, respectively, from PTC and from a metastasis of a well-differentiated PTC (Ribeiro et al., 2008); the 8505C (Ito et al., 1993) cell line from ATC. These six cell lines are among those mostly used for thyroid cancer research. The mutational status of each of these cell lines was investigated and compared to the literature. Next, pathways were further investigated by using our previously obtained data on mRNA expression patterns of these cell lines that had been compared to primary cultured normal human thyrocytes (van Staveren et al., 2007). Previously we showed that mRNA expression profiles from different human thyroid tumor cell lines, including cell lines of the present study, evolved in vitro into similar phenotypes with mRNA expression profiles closer to undifferentiated in vivo thyroid tumors (ATC) than to differentiated thyroid cancers (PTC and FTC) (van Staveren et al., 2007). To investigate this further, a correlation analysis between mRNA expression profiles from cell lines and mRNA expression profiles from each in vivo thyroid tumor type was performed by using three datasets generated in two different laboratories that were derived from different microarray platforms.

Materials and Methods

Cell Culture

Cell lines were obtained from various laboratories. The WRO, FTC133, BCPAP, and 8505C cell lines were received from Prof. G. Brabant (Medizinische Hochschule, Hannover, Germany); the TPC1 cell line from Prof. M. Mareel (University of Ghent, Belgium); the K1 cell line from Dr. Zaruhi Poghosyan (Cardiff University, Scholl of Medecine, Cardiff, UK). The BCPAP, TPC1 and 8505C cell lines were cultured in RPMI 1640 (w L-glutamate) completed by foetal bovine serum (10%), penicillin/streptomycin (2%) and amphotericine B (1%); the FTC133, WRO and K1 cell lines were cultured in Dulbecco's Modified Eagle Medium, Nutrient mixture F-12 (1:1, by volume) completed by foetal bovine serum (10%), penicillin/streptomycin (2%), amphotericine B (1%). All culture reagents were purchased from Gibco (Paisely, UK).

RNA Extraction

After removal of the culture medium, 2 ml of TRIzol Reagent (Ambion, Austin, USA) was added to the cells grown in 10 cm dishes. Total RNA was extracted according to the manufacturer's instructions, followed by purification on miRNeasy columns (QIAGEN, Hilden, Germany). RNA was spectrophotometrically quantified, and its integrity was verified by automated gel electrophoresis (Experion, Bio-Rad, Hercules, USA).

RT-PCR

After a DNase treatment with DNase I Amplification Grade kit (Invitrogen, Carlsbad, USA), 1 μg of total RNA was used for reverse transcription using hexamers [(3.6 μg/μl) (Roche, Basel, Switzerland)] and reverse transcriptase (Superscript II RNase H Reverse Transcriptase kit, Invitrogen). The PCR reactions were performed with the recombining Taq DNA polymerase kit (Invitrogen). Each PCR reaction was performed in the presence of 5 μl 10X PCR buffer, 1.5 μl MgCl₂ (50 mM), 1 μl dNTP mix (10 mM, Invitrogen), 1 μl of the forward and the reverse primer (10 μM each), 0.4 μl Taq DNA polymerase (5U/μl), 2 μl DNA from the RT reaction and 38.1 μl of water. Primers and PCR conditions are detailed in Table 1.

TABLE 1

Table 1. Targeted exons, sequences of the primers used for PCR amplification, amplicon length and PCR conditions for each investigated gene.

Sequencing

PCR products were purified with the QIAquick PCR purification kit (QIAGEN) according to the manufacturer's instructions. Sequencing was performed with the BigDye Terminator V3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, USA) with the sequencer ABI PRISM 3130 (Applied Biosystems) and the genetic analysis program 3130-XI. Sequences were analyzed with BLAST (Basic Local Alignment Search Tool program) from NCBI (http://blast.ncbi.nlm.nih.gov/Blast.cgi).

Pre-Processing of the Different Datasets

A merged file was created which contained three different datasets derived from three different microarray platforms. The first dataset consisted of five cell lines hybridized on IRIBHM custom slides, using primary cultured normal thyrocytes as a reference (www.ulb.ac.be//medecine/iribhm/microarray/data/) (van Staveren et al., 2007). The second dataset was composed of 11 ATC, 49 PTC, and 45 normal thyroid tissues hybridized on the Affymetrix HG U133 Plus 2.0 platform (GSE33630). The third dataset contained 4 ATC, 51 PTC, 13 FTC, and 4 normal thyroid tissues hybridized on the Affymetrix HG U133A platform (Giordano et al., 2005) (GSE27155). Each platform has been validated in previous studies [Giordano et al., 2005; van Staveren et al., 2006; Hébrant et al., 2012; Dom et al. (submitted)]. For both independent Affymetrix platforms, mRNA expression profiles of normal tissues were pooled as a reference and an mRNA expression ratio between each tumor and the normal reference was calculated for each gene. For the comparison of cell lines with in vivo thyroid tumors, only genes were included that was expressed in both cell lines and primary thyrocytes. Finally, the merged file contained the mRNA expressions of 2392 genes, expressed in log₂, present on the three different microarray platforms.

DAVID Pathway Analysis

In order to investigate the biological meaning of the modulated mRNAs, Database for Annotation, Visualization and Integrated Discovery (DAVID pathway analysis) (Huang da et al., 2009) was used. For the cell line analyses, only genes showing a regulation of at least |0.58 | in log₂, i.e., 1.5-fold or more were included. To calculate the background, all the genes present on the IRIBHM custom slides were used. For the comparison between ATC and cell lines, only significant commonly regulated genes were considered. The background was calculated using all the genes present in the merged file. For both analyses, a pathway was selected when its False Discovery Rate (FDR) was <5%. The analysis was performed separately for up-regulated and down-regulated genes.

Correlation of mRNA Expression Between in vitro Human Thyroid Cancer Cell Lines and in vivo Human Thyroid Tumors

Spearman's correlation coefficients were computed using R (version 2.14.1). Based on the merged file, correlation coefficients between gene expression profiles of ATC and thyroid cancer cell lines were compared to correlation coefficients between gene expression profiles of PTC or FTC and thyroid cancer cell lines, using a Mann–Whitney U-test. This assessed whether ATC and cell lines had more similar mRNA expression profiles than FTC or PTC and cell lines. Comparisons of those coefficients were made between samples hybridized on the same microarray platform. 2D-multidimensional scaling was plotted based on those Spearman correlation coefficients using the R function isoMDS from the package MASS v.7.3–18.

Analysis of the Genes Commonly and Oppositely Modulated in ATC in vivo and Cancer Cell Lines

Scores and q-values were computed based on a slightly modified version of Rank Products (Breitling et al., 2004) to find genes regulated in ATC compared to normal tissues and in cell lines compared to primary cultured normal thyrocytes. Rank Products is a non-parametric method allowing to study differential expression and to integrate the information of different datasets regardless of the potential technical differences, e.g., from different laboratories or different microarray platforms. It has been shown to outperform SAM (Zang et al., 2007) and is more powerful with fewer samples than t-test based statistics. The genes with a q−value ≤ 0.05 were selected from ATC and cell lines. Genes that were regulated in the same direction in ATC and cell lines were kept and are shown in a heatmap using R library: gplots v.2.10.1 (http://cran.r-project.org/web/packages/gplots/index.html). Similarly, genes that had significant regulations in opposite directions between ATC and cell lines were further analyzed.

Results

Mutational Status

The mutational status was investigated for known thyroid oncogenes and tumor suppressor genes in six human thyroid cancer cell lines. For each cell line, the presence of point mutations in BRAF, NRAS, HRAS, KRAS, TP53, PI3KCA, PTEN genes was explored as well as the presence of RET/PTC1, RET/PTC3, and PAX8/PPARγ rearrangements. Therefore, primers were designed to target the point mutations and the rearrangements most frequently found in thyroid carcinomas (Kondo et al., 2006) and each mutation was verified by sequencing both DNA strands. The results are summarized in Table 2. Neither NRAS, HRAS, KRAS point mutations nor RET/PTC3 or PAX8/PPARγ gene rearrangements were detected in any of the six cell lines. As reported earlier, FTC133 cells showed PTEN and TP53 mutations, BCPAP and 8505C cells showed BRAF and TP53 mutations, and TPC1 cells harbored the RET/PTC1 gene rearrangement. WRO cells presented the heterozygous BRAF T1799A (V600E) mutation and K1 the heterozygous PIK3CA G1624A (E542K) mutation, in addition to the BRAF mutation (Figure 1).

TABLE 2

Table 2. Mutational status of six human thyroid cancer cell lines derived from three types of thyroid carcinomas.

FIGURE 1

Figure 1. (A) Sequence of BRAF in WRO cells. The arrow shows the T1799A mutation. (B) Sequence of PI3KCA in K1 cells. The arrow indicates the G1624A mutation.

Pathway Analysis of the Cell Lines

To define the modulated pathways in thyroid cancer cell lines, the genes from previously published mRNA expression profiles from five thyroid cancer cell lines (FTC133, BCPAP, 8505C, TPC1, WRO) were further analysed (van Staveren et al., 2007). These data reported the differential mRNA expression between each cell line and a pool of human primary cultured normal thyrocytes. DAVID pathway analyses were performed on each individual cell line (Supplementary information A). Four of the five cell lines, i.e., WRO, BCPAP, TPC1, and 8505C, showed commonly regulated pathways. An upregulation of ribosome biosynthesis and translation was detected in these four cell lines, and in addition, expression of genes involved in the cell cycle was increased in WRO and BCPAP cells. Expression of genes involved in DNA replication was up-regulated in both BCPAP and 8505C cells. Furthermore, in all cell lines, expression of genes involved in MHC class II biosynthesis was down-regulated and in four cell lines (not in 8505C) the expression of genes involved in negative regulation of cell death/apoptosis was also down-regulated. In both 8505C and TPC1 cells, expression of genes involved in the immune response was decreased and for the latter, a decrease in the protease inhibitor pathway was also observed. In contrast, the FTC133 cell line showed some differences in the regulated pathways. Although the downregulation of the MHC class II and the negative regulation of cell death/apoptosis were present, no increase of expression of genes involved in ribosome biosynthesis and translation were observed. Instead, FTC133 cells showed a strong upregulation of global RNA processing.

Correlations between mRNA Expression Profiles of in vitro Thyroid Cancer Cell Lines and mRNA Expression Profiles of in vivo Thyroid Tumors

Two different datasets of mRNA expression profiles from in vivo human thyroid tumors were correlated to the mRNA profiles of the WRO, FTC133, BCPAP, TPC1, and 8505C cell lines. One dataset was composed of 11 ATC and 49 PTC hybridized on Affymetrix HG U133 Plus 2.0 arrays. The other dataset consisted of 4 ATC, 13 FTC, and 51 PTC hybridized on Affymetrix HG U133A arrays. For both independent datasets, the Spearman correlation coefficient was the highest between cell lines and ATC compared to the cell lines and PTC or FTC (Figures 2A,B). Furthermore, there was a significant difference between the correlations of ATC versus cell lines compared to those calculated with the PTC or FTC profiles (Figures 2A,B).

FIGURE 2

Figure 2. Boxplots of correlations between cell lines and in vivo thyroid tumors. (A and B) Correlations were calculated based on 2,392 mRNA expressions comparing human thyroid cancer cell lines (CL) with different types of in vivo human thyroid tumors by using three datasets. The first was composed of five thyroid cancer cell lines (WRO, FTC133, BCPAP, TPC1, 8505C) hybridized on IRIBHM custom slides. The second was composed of 11 ATC and 49 PTC hybridized on Affymetrix HG U133 Plus 2.0 (A). The third was composed of 4 ATC, 13 FTC, and 51 PTC hybridized on Affymetrix HG U133A (B). P-values indicate whether ATC and cell lines had more similar profiles than FTC or PTC and cell lines. P-values were significant for both independent datasets. (C) 2D-multidimensional scaling representation of the mRNA expression profiles of cell lines, ATC, PTC, and FTC based on 2392 mRNA expressions.

These results are in line with previous results which showed by 2D-multidimensional scaling analysis the highest similarity between cell lines and ATC, compared to cell lines with PTC or FTC (van Staveren et al., 2007). Thus, this already observed relation could be confirmed in the present study using more tumor samples that had been hybridized on different microarray platforms. As shown in Figure 2C, although a platform effect was visible, in both independent datasets ATC were closer to cell lines.

Analysis of the Genes Commonly and Oppositely Modulated in ATC in vivo and Cancer Cell Lines

To further investigate the similarity between ATC and cell lines, their mRNA expression profiles, obtained on different microarray platforms, were compared. Rank Products and SAM were used to find genes significantly regulated in the same direction or in the opposite direction comparing ATC with cell lines. Sixty-one genes were found to be up-regulated in both ATC and cell lines and eighty-four genes were found to be commonly down-regulated (Figure 3). Eighteen genes were regulated in an opposite direction between ATC and cell lines.

FIGURE 3

Figure 3. Heatmap of 145 mRNA expressions in thyroid cancer cell lines (CL) and in human thyroid tumors. These genes were the most commonly up and down-regulated genes in ATC and in the cell lines. The level of modulation of these genes in PTC and FTC are also shown. Profiles of the five cell lines are underlined in black. Profiles of the samples hybridized on Affymetrix HG U133 Plus 2.0 platform and on Affymetrix HG U133A platform are underlined in blue and violet respectively. Genes are indicated by their official gene symbols.

To gain more insight into those regulations, a DAVID pathway analysis was performed (Supplementary information B). Some common up-regulated genes in ATC and cell lines are involved in DNA replication and DNA metabolic processes, e.g., TOP2A, MCM7, CHAF1A, or HMGB2. Most of the 18 genes that were regulated in an opposite direction were down-regulated in cell lines but up-regulated in ATC. They are involved in MHC class II pathway or in cell adhesion, e.g., HLA-DPB1, HLA-DQB1, HLA-DRA, ICAM1, CCL2, or COL6A3.

Discussion

In this work, the aim was to further investigate the characteristics for different human thyroid cancer cell lines. First, genetic alterations in these cell lines were confirmed. Second, regulated pathways in these cell lines were further explored. Third, correlations, based on gene expressions, between in vitro human thyroid cell lines and in vivo human thyroid tumors were analysed. Finally, the similarity between cell lines and ATC was further explored on the gene level.

Mutational Status

We have investigated the mutational status of 10 commonly mutated genes in thyroid cancer by using six different cell lines. STR analyses on these cell lines (van Staveren et al., 2007 and Supplementary information C) showed that they were identical to those published by Schweppe and coauthors (Schweppe et al., 2008). No differences with reported literature data were observed for the mutational status of these genes in BCPAP, FTC133, TPC1, and 8505C cell lines. However, a difference could be observed for the status of BRAF in WRO and PI3KCA in K1 cells.

For WRO cells described as a FTC cell line, the BRAF T1799A (V600E) mutation was detected although these cells were originally described as BRAF wild-type (Namba et al., 2003; Xu et al., 2003) and are still frequently used as such today (Palona et al., 2006; Liu et al., 2007; Burrows et al., 2010; Hou et al., 2010). The presence of a BRAF mutation in these cells was nevertheless reported by Schweppe and coworkers (Schweppe et al., 2008) and they demonstrated that the cell line was not contaminated by another cell line harboring this mutation. It has been suggested that two different WRO cell lines were distributed. The WRO cells from Schweppe's laboratory has the same DNA profile (STR analysis) than ours (van Staveren et al., 2007; Schweppe et al., 2008), but is different from that of the European Collection for Cell Culture. The presence of this mutation in WRO cells was unexpected because this cell line has been derived from a FTC, and the V600E BRAF mutation is described as a specific characteristic for PTC. This might suggest that, in spite of its initial diagnosis of follicular cancer, the tumor at the origin of the creation of the WRO cell line would have rather been a papillary type. The presence of a mixture of BRAF mutated and non-mutated cells in the original tumor sample cannot be excluded. However, it may also be proposed that this cell line acquired this mutation in vitro, during its successive divisions in an artificial environment of selection (van Staveren et al., 2007, 2009). These observations question the representativity and the utility of this cell line as a model for human FTC and show the necessity of a systematic check of the cell lines. Furthermore, the V600E BRAF mutation is not specific of thyroid cancer. The origin of the WRO cell line as a thyroid cancer cell line is not completely certified.

For the K1 cell line, a mutation in the PI3KCA gene has been reported previously. By sequencing, the presence of a G1624A substitution in codon 542 (E542K) was confirmed in accordance with previously published data (Garcia-Rostan et al., 2005; Meireles et al., 2007). In another publication, Ricarte-Filho et al. (Ricarte-Filho et al., 2009) reported the presence of a E545K (G1633A) mutation in these cells. We could observe a glutamic acid at position 545. Thus, the G1633A mutation might potentially exist but was never described as a new mutation by Ricarte–Filho et al.

Pathway Analysis of the Cell Lines

Pathway analysis of mRNA expression profiles in WRO, FTC133, BCPAP, TPC1, and 8505C cell lines was performed. In general, the cell lines abrogate gene expressions which are not necessary for their propagation in in vitro cultures and express genes that promote in growth and survival. The common loss of MHC2 biosynthesis and, in some cell lines, of the immune response, can be explained because they would be not necessary for a cell in culture. The common downregulation of the negative regulation of the cell death/apoptosis pathway might be explained by an ineffective feedback resulting from an excess of proliferation.

Except for FTC133 cells, a common overexpression of ribosomal biosynthesis and translation pathways were observed. This may reflect the increase in proliferation of the cancer cells induced by gene mutations. However, some recent studies suggest that the up-regulation of the ribosomal biosynthesis may induce cancer by a downregulation of the cell tumor suppressor potential (Montanaro et al., 2012). For our data, the two explanations are possible. However, the growth rate of the FTC133 cell line is higher than for the other cell lines [Dom et al. (submitted)]. For this cell line, an overexpression of the global RNA processing including alternative splicing was observed. This might allow the strong increase in proliferation of this cell line.

Correlations between mRNA Expression Profiles of in vitro Thyroid Cancer Cell Lines and mRNA Expression Profiles of in vivo Thyroid Tumors

The correlations of the mRNA expression profiles between cell lines and ATC were significantly higher than those compared to other in vivo differentiated thyroid tumors (PTC, FTC). By using data derived from different microarray platforms and from two different laboratories, thereby including a larger number of in vivo tumors, the current data confirm our previous finding that mRNA expression profiles from human thyroid tumor cell lines with different origins are the closest to undifferentiated in vivo tumors (van Staveren et al., 2007). This is explained by an in vitro evolution of the cell lines that should be taken into account when extrapolating results obtained from these cells. An in vitro cell line is the result of the selection of the most rapidly proliferating and resistant cells during numerous passages. Cell lines contain one cell type, whereas in vivo tumors are complex mixtures of various cell types. Therefore, it might be argued that the comparisons between thyroid cancer cell lines and in vivo thyroid tumors could be affected by the presence of different cell types such as for instance lymphocytes, endothelial cells. However, in this study only genes were compared that were expressed in thyrocytes, i.e., were expressed in both thyroid cancer cell lines and in primary thyrocytes. This should reduce the impact of the potential presence of different cell types but, unfortunately, it cannot assure the complete elimination of a background due to commonly expressed genes in the different cell types.

Analysis of the Genes Commonly and Oppositely Modulated in ATC in vivo and Cancer Cell Lines

Genes commonly up-regulated in ATC and cell lines were involved in DNA replication and DNA metabolic processes. The correlation, based on mRNA expression profiles, between cell lines and ATC thus seems to be due to the strong proliferation capacity of these two types of cells. Indeed, ATC is one of the most proliferating human tumors and human thyroid cancer cell lines have a strong proliferation capacity. In our data this property was linked to ribosomal biosynthesis and translation pathways for four cell lines. As expected, the genes found regulated in an opposite direction between ATC and cell lines might reflect the differences between in vivo tumors and in vitro cultures. They are involved in the MHC class II pathway or in cell adhesion. As mentioned above in the pathway analysis of the cell lines, they abrogate gene expressions which are not necessary for their survival and growth in in vitro cultures. The increase of those genes in ATC, and the implication of this increase to growth and survival of the tumor in vivo, would be interesting to explore. However, there are a number of factors which could also explain this increase. For instance, the potential presence of immune cells, which express MHC class II in the tumor or the activation of the MHC class II pathway in these cells by an inflammatory process, might explain an increase in these genes in ATC. In line with this, thyrocytes indeed can express MHC class II after an inflammatory stimulation and inflammation is often associated with tumor development. Furthermore, the increase of expression might also be due to a combination of the factors.

Conclusions

In this work, the aim was to further investigate the characteristics of human thyroid cancer cell lines derived from different origins. First, genetic alterations in six cell lines were in line with literature data and furthermore, the BRAF mutation in WRO cells questions this cell line as a model for FTC, indeed this mutation does not occur in FTC. Second, in five cell lines the MHC class II pathway was down-regulated and in four of them, ribosome biosynthesis and translation pathways were up-regulated. Third, by using microarray data from ATC, PTC, and FTC, we have preliminary found that mRNA expression profiles of human thyroid cancer cell lines, originally derived from ATC, poorly differentiated PTC, FTC, or PTC, were closer to ATC. Specifically, we found genes commonly up-regulated in ATC involved in DNA replication which were also up-regulated in the thyroid cancer cell lines, in accordance with their higher proliferation rate in vitro. In summary, these results would suggest that some human PTC or FTC derived cell lines (i.e., TPC-1) might have partially lost their original DNA synthesis/replication regulation mechanisms during their in vitro cell adaptation/evolution. Also, further comparisons using microarray data set from poorly differentiated thyroid cancers are needed and could be a near future goal in our laboratory.

During our study, the origin of the WRO cell line as thyroid cancer cell line has been questioned. However, the pathway analysis and the comparisons of expression profiles showed similar results for WRO cell line and the other thyroid cancer cell lines. To address the origin of the WRO cell line, more complete analyses (including cell lines from different cancer types) could be a near future goal in our laboratory.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by the “Fonds de la Recherche Scientifique Médicale” (FRSM), Fondation contre le Cancer, Télévie, F.R.I.A., Fondation Van Buuren and Fondation Rose et Jean Hoguet. Authors want to thank Professor J. E Dumont for his conceptual support. Authors want to thank Bernadette Bournonville and Chantal Degraef for her technical support.

Supplementary Material

The Supplementary Material for this article can be found online at http://www.frontiersin.org/Cancer_Endocrinology/10.3389/fendo.2012.00133/abstract

Supplementary information A. DAVID pathway analysis of the genes modulated in each human thyroid cancer cell line. Up, analysis of pathways using all the upregulated genes (>0.58 log2) in the cell line. Down: analysis of pathways using all the downregulated genes (<−0.58 log2) in the cell line.

Supplementary information B. DAVID pathway analysis of the genes commonly and oppositely modulated in human thyroid cancer cell lines and ATC. Up: analysis of pathways using the 61 up-regulated genes that were upregulated in both. Down: analysis of pathways using the 84 genes down-regulated in both ATC and cell lines. Modulation: analysis of pathways using the 18 genes that were regulated in an opposite direction comparing cell lines with ATC.

Supplementary information C. DNA profiling by STR analysis of the K1 cell line.

References

Breitling, R., Armengaud, P., Amtmann, A., and Herzyk, P. (2004). Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Lett. 573, 83–92.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Burrows, N., Resch, J., Cowen, R. L., von, W. R., Hoang-Vu, C., West, C. M., et al. (2010). Expression of hypoxia-inducible factor 1 alpha in thyroid carcinomas. Endocr. Relat. Cancer 17, 61–72.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Catalano, M. G., Poli, R., Pugliese, M., Fortunati, N., and Boccuzzi, G. (2010). Emerging molecular therapies of advanced thyroid cancer. Mol. Aspects Med. 31, 215–226.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Challeton, C., Branea, F., Schlumberger, M., Gaillard, N., de Vathaire, F., Badie, C., et al. (1997). Characterization and radiosensitivity at high or low dose rate of four cell lines derived from human thyroid tumors. Int. J. Radiat. Oncol. Biol. Phys. 37, 163–169.

Pubmed Abstract | Pubmed Full Text

Estour, B., Van Herle, A. J., Juillard, G. J., Totanes, T. L., Sparkes, R. S., Giuliano, A. E., et al. (1989). Characterization of a human follicular thyroid carcinoma cell line (UCLA RO 82 W-1). Virchows Arch. B Cell Pathol. Incl. Mol. Pathol. 57, 167–174.

Pubmed Abstract | Pubmed Full Text

Fabien, N., Fusco, A., Santoro, M., Barbier, Y., Dubois, P. M., and Paulin, C. (1994). Description of a human papillary thyroid carcinoma cell line. Morphologic study and expression of tumoral markers. Cancer 73, 2206–2212.

Pubmed Abstract | Pubmed Full Text

Garcia-Rostan, G., Costa, A. M., Pereira-Castro, I., Salvatore, G., Hernandez, R., Hermsem, M. J., et al. (2005). Mutation of the PIK3CA gene in anaplastic thyroid cancer. Cancer Res. 65, 10199–10207.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Giordano, T. J., Kuick, R., Thomas, D. G., Misek, D. E., Vinco, M., Sanders, D., et al. (2005). Molecular classification of papillary thyroid carcinoma: distinct BRAF, RAS, and RET/PTC mutation-specific gene expression profiles discovered by DNA microarray analysis. Oncogene 24, 6646–6656.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Goretzki, P. E., Frilling, A., Simon, D., and Roeher, H. D. (1990). Growth regulation of normal thyroids and thyroid tumors in man. Recent Results Cancer Res. 118, 48–63.

Pubmed Abstract | Pubmed Full Text

Harlin, H., and Gajewski, T. F. (2008). Diagnosis and treatment of mycoplasma-contaminated cell cultures. Curr. Protoc. Cytom. 43, A.3C.1–A.3C.7.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hébrant, A., Dom, G., Dewaele, M., Andry, G., Trésallet, C., Leteurtre, E., et al. (2012). mRNA expression in papillary and anaplastic thyroid carcinoma: molecular anatomy of a killing switch. PLoS ONE 7:e37807. doi: 10.1371/journal.pone.0037807

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hou, P., Bojdani, E., and Xing, M. (2010). Induction of thyroid gene expression and radioiodine uptake in thyroid cancer cells by targeting major signaling pathways. J. Clin. Endocrinol. Metab. 95, 820–828.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Huang da, W., Sherman, B. T., and Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ito, T., Seyama, T., Iwamoto, K. S., Hayashi, T., Mizuno, T., Tsuyama, N., et al. (1993). In vitro irradiation is able to cause RET oncogene rearrangement. Cancer Res. 53, 2940–2943.

Pubmed Abstract | Pubmed Full Text

Kondo, T., Ezzat, S., and Asa, S. L. (2006). Pathogenetic mechanisms in thyroid follicular-cell neoplasia. Nat. Rev. Cancer 6, 292–306.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lacroix, L., Mian, C., Barrier, T., Talbot, M., Caillou, B., Schlumberger, M., et al. (2004). PAX8 and peroxisome proliferator-activated receptor gamma 1 gene expression status in benign and malignant thyroid tissues. Eur. J. Endocrinol. 151, 367–374.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Liu, D., Liu, Z., Jiang, D., Dackiw, A. P., and Xing, M. (2007). Inhibitory effects of the mitogen-activated protein kinase kinase inhibitor CI-1040 on the proliferation and tumor growth of thyroid cancer cells with BRAF or RAS mutations. J. Clin. Endocrinol. Metab. 92, 4686–4695.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Meireles, A. M., Preto, A., Rocha, A. S., Rebocho, A. P., Maximo, V., Pereira-Castro, I., et al. (2007). Molecular and genotypic characterization of human thyroid follicular cell carcinoma-derived cell lines. Thyroid 17, 707–715.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Montanaro, L., Trere, D., and Derenzini, M. (2012). Changes in ribosome biogenesis may induce cancer by down-regulating the cell tumor suppressor potential. Biochim. Biophys. Acta 1825, 101–110.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Namba, H., Nakashima, M., Hayashi, T., Hayashida, N., Maeda, S., Rogounovitch, T. I., et al. (2003). Clinical implication of hot spot BRAF mutation, V599E, in papillary thyroid cancers. J. Clin. Endocrinol. Metab. 88, 4393–4397.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Nucera, C., Porrello, A., Antonello, Z. A., Mekel, M., Nehs, M. A., Giordano, T. J., et al. (2010). B-Raf(V600E) and thrombospondin-1 promote thyroid cancer progression. Proc. Natl. Acad. Sci. U.S.A. 107, 10649–10654.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Palona, I., Namba, H., Mitsutake, N., Starenki, D., Podtcheko, A., Sedliarou, I., et al. (2006). BRAFV600E promotes invasiveness of thyroid cancer cells through nuclear factor kappaB activation. Endocrinology 147, 5699–5707.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ribeiro, F. R., Meireles, A. M., Rocha, A. S., and Teixeira, M. R. (2008). Conventional and molecular cytogenetics of human non-medullary thyroid carcinoma: characterization of eight cell line models and review of the literature on clinical samples. BMC Cancer 8:371. doi: 10.1186/1471-2407-8-371

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ricarte-Filho, J. C., Ryder, M., Chitale, D. A., Rivera, M., Heguy, A., Ladanyi, M., et al. (2009). Mutational profile of advanced primary and metastatic radioactive iodine-refractory thyroid cancers reveals distinct pathogenetic roles for BRAF, PIK3CA, and AKT1. Cancer Res. 69, 4885–4893.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Schweppe, R. E., Klopper, J. P., Korch, C., Pugazhenthi, U., Benezra, M., Knauf, J. A., et al. (2008). Deoxyribonucleic acid profiling analysis of 40 human thyroid cancer cell lines reveals cross-contamination resulting in cell line redundancy and misidentification. J. Clin. Endocrinol. Metab. 93, 4331–4341.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sipos, J. A., and Mazzaferri, E. L. (2010). Thyroid cancer epidemiology and prognostic variables. Clin. Oncol. (R. Coll. Radiol.) 22, 395–404.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Tanaka, J., Ogura, T., Sato, H., and Hatano, M. (1987). Establishment and biological characterization of an in vitro human cytomegalovirus latency model. Virology 161, 62–72.

Pubmed Abstract | Pubmed Full Text

van Staveren, W. C., Solis, D. W., Delys, L., Duprez, L., Andry, G., Franc, B., et al. (2007). Human thyroid tumor cell lines derived from different tumor types present a common dedifferentiated phenotype. Cancer Res. 67, 8113–8120.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

van Staveren, W. C., Solis, D. W., Delys, L., Venet, D., Cappello, M., Andry, G., et al. (2006). Gene expression in human thyrocytes and autonomous adenomas reveals suppression of negative feedbacks in tumorigenesis. Proc. Natl. Acad. Sci. U.S.A. 103, 413–418.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

van Staveren, W. C., Solis, D. Y., Hebrant, A., Detours, V., Dumont, J. E., and Maenhaut, C. (2009). Human cancer cell lines: experimental models for cancer cells in situ? For cancer stem cells? Biochim. Biophys. Acta 1795, 92–103.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wang, H., Mannava, S., Grachtchouk, V., Zhuang, D., Soengas, M. S., Gudkov, A. V., et al. (2008). c-Myc depletion inhibits proliferation of human tumor cells at various stages of the cell cycle. Oncogene 27, 1905–1915.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Xing, M. (2010). Prognostic utility of BRAF mutation in papillary thyroid cancer. Mol. Cell Endocrinol. 321, 86–93.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Xu, X., Quiros, R. M., Gattuso, P., Ain, K. B., and Prinz, R. A. (2003). High prevalence of BRAF gene mutation in papillary thyroid carcinomas and thyroid tumor cell lines. Cancer Res. 63, 4561–4567.

Pubmed Abstract | Pubmed Full Text

Yeung, S. C., She, M., Yang, H., Pan, J., Sun, L., and Chaplin, D. (2007). Combination chemotherapy including combretastatin A4 phosphate and paclitaxel is effective against anaplastic thyroid cancer in a nude mouse xenograft model. J. Clin. Endocrinol. Metab. 92, 2902–2909.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Zang, S., Guo, R., Zhang, L., and Lu, Y. (2007). Integration of statistical inference methods and a novel control measure to improve sensitivity and specificity of data analysis in expression profiling studies. J. Biomed. Inform. 40, 552–560.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text