| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Division of Endocrinology, Diabetes and Metabolism (R.E.S., J.P.K., U.P., B.R.H.), Department of Medicine and University of Colorado Cancer Center DNA Sequencing and Analysis Core (C.K.), University of Colorado Cancer Center, University of Colorado Denver, Aurora, Colorado 80045; Human Oncology and Pathogenesis Program (M.B., J.A.K., J.A.F.), Memorial Sloan-Kettering Cancer Center, New York, New York, 10021; and Division of Endocrinology (R.C.S.), Internal Medicine Department, Department of Cancer Biology (L.A.M., J.A.C.), Mayo Clinic Comprehensive Cancer Center, Jacksonville, Florida 32224
Address all correspondence and requests for reprints to: Rebecca E. Schweppe, Ph.D., Division of Endocrinology, Diabetes and Metabolism, University of Colorado Denver, Aurora, Colorado 80045. E-mail: Rebecca.Schweppe{at}ucdenver.edu.
| Abstract |
|---|
|
|
|---|
Objective: We evaluated 40 reported thyroid cancer-derived cell lines using short tandem repeat and single nucleotide polymorphism array analysis.
Results: Only 23 of 40 cell lines tested have unique genetic profiles. The following groups of cell lines are likely derivatives of the same cell line: BHP5-16, BHP17-10, BHP14-9, and NPA87; BHP2-7, BHP10-3, BHP7-13, and TPC1; KAT5, KAT10, KAT4, KAT7, KAT50, KAK1, ARO81-1, and MRO87-1; and K1 and K2. The unique cell lines include BCPAP, KTC1, TT2609-C02, FTC133, ML1, WRO82-1, 8505C, SW1736, Cal-62, T235, T238, Uhth-104, ACT-1, HTh74, KAT18, TTA1, FRO81-2, HTh7, C643, BHT101, and KTC-2. The misidentified cell lines included the DRO90-1, which matched the melanoma-derived cell line, A-375. The ARO81-1 and its derivatives matched the HT-29 colon cancer cell line, and the NPA87 and its derivatives matched the M14/MDA-MB-435S melanoma cell line. TTF-1 and Pax-8 mRNA levels were determined in the unique cell lines.
Conclusions: Many of these human cell lines have been widely used in the thyroid cancer field for the past 20 yr and are not only redundant, but not of thyroid origin. These results emphasize the importance of cell line integrity, and provide the short tandem repeat profiles for a panel of thyroid cancer cell lines that can be used as a reference for comparison of cell lines from other laboratories.
| Introduction |
|---|
|
|
|---|
A comprehensive analysis of commonly used human thyroid cancer cell lines is currently lacking, despite three recent reports that have identified redundant cell lines (11, 12, 13). Meireles et al. (11) reported that the papillary thyroid cancer (PTC)-derived FB2 and TPC1 cell lines, which both harbor the RET/PTC1 rearrangement, are genetically identical, and van Staveren (12) and Rocha (13) et al. reported redundancy among the KAK-1, KAT10, and KAT4 cell lines, which all share a point mutation in E-cadherin, a rare event in thyroid cancer. We used short tandem repeat (STR) profiling and single nucleotide polymorphism (SNP) array analysis to systematically evaluate the integrity of the thyroid cancer cell lines in our laboratories (3). STR profiling was originally developed for forensic applications (14) and has been established as an international reference standard for genetically profiling cell lines (3). All major cell line repositories use STR profiling to authenticate their cell lines, thus increasing the power of this approach in the form of searchable STR profile databases. In this study we genetically analyzed 40 thyroid cancer cell lines. We found significant redundancy among the cell lines and misidentification of some cell lines that are likely not of thyroid cancer origin.
| Materials and Methods |
|---|
|
|
|---|
Cell lines, as well as year of first publication, and reported mutation status are listed in Table 1
, and were obtained from the primary source or the American Type Culture Collection (ATCC) with the exception of the following. BCPAP cells were kindly provided by Dr. M. Santoro (Medical School, University "Federico II " of Naples, Naples, Italy). K1 and K2 cell lines were provided by Dr. Wynford-Thomas (Cardiff University, Cardiff, UK). SW1736, C643, HTh7, and HTh74 cells were provided by Dr. K. Ain (University of Kentucky, Lexington, KY) with permission from Dr. N.-E. Heldin (University Hospital, Uppsala, Sweden), and the TPC1 cells were kindly provided by Dr. S. Jhiang (Ohio State University, Columbus, OH). Although we were unable to obtain these cell lines from the source, these cell lines matched their respective STR profiles shown in previously published studies (11, 15), and/or in the German Collection of Microorganisms and Cell Cultures (DSMZ) and European Collection for Cell Cultures (ECACC) databases at all loci tested in common (Table 2
). In most of these cases, we were unable to obtain cell lines from the source due to a lack of current stocks. All cell lines were grown at 37 C with 5% CO2 in a humidified environment (media used for individual cell lines are available upon request). To avoid cross-contamination, cell lines were cultured separately in a sterile tissue culture hood, which was UV irradiated between each cell line, and a separate bottle of media was used for each cell line.
|
|
Genomic DNA was prepared separately for each cell line in a clean hood using DNeasy Tissue kit (QIAGEN, Inc., Valencia, CA). To avoid cross-contamination, Pipetman were UV irradiated between each sample, and filtered pipette tips were used throughout the DNA preparation. The earliest or original passage of each cell line received in our laboratories was analyzed in each case. The Applied Biosystems Profiler Plus kit (ABI no. 4303326; Foster City, CA) was used with modifications. Genomic DNA (
1 ng) was used as template, and cycling conditions were those suggested by the manufacturer, and multiple aliquots of diluted PCRs were electrophoresed in an ABI 3730 genetic analyzer using the POP7 polymer and 50 cm. The data were collected using ABIs Data Collection software version 3.0, and the sizes of the PCR products for the different STR loci were determined using GeneMapper 4.0. The identities of the products were determined by comparing the product sizes with those of the different alleles present in the standard allelic ladder in the Profiler Plus kit. A no DNA template control was included in all stages of the STR analysis to exclude contamination. The profiles were compared visually or by using a searchable database (Korch, C. and J. West, personal communication) containing in-house data and data downloaded from the ATCC web site (www.atcc.org/common/cultures/str.cfm). The ECACC and DSMZ databases were manually searched with the STR profiles.
SNP arrays
The genetic identity of some of the commonly used thyroid cancer cell lines was further characterized by SNP genotyping, using Affymetrix 250K SNP arrays (Santa Clara, CA). Each sample was restricted with NspI, labeled, and hybridized to the Affymetrix GeneChip Human Mapping 250K SNP array, according to the manufacturers instructions. Genotyping calls were generated using GType 1.1.
BRAF mutation analysis
Exon 15 of BRAF was amplified from genomic DNA (120 ng) using standard touchdown PCR methods with the following primers: BRAF-forward, 5'-TCATAATGCTTGCTCTGATAGGA-3'; and BRAF-reverse, 5'-GGCCAAAAATTTAATCAGTGGA-3'. Residual primers and deoxynucleotide triphosphates were removed before sequencing using the ExoSAP-IT kit (USB Corp., Cleveland, OH). DNA samples were sequenced using the BRAF-For primer by the University of Colorado Cancer Center DNA Sequencing & Analysis Core (http://loki.uchsc.edu).
RET/PTC1 analysis
Total RNA was isolated from cells/tissues using RNeasy Mini Kit (QIAGEN) and reverse transcribed using Moloney Murine leukemia virus reverse transcriptase (Life Technologies, Inc., Gaithersburg, MD). The primers used for RET/PTC1 amplification were: RET/PTC1 (F) 5'-ATTGTCATCTCGCCGTTC-3' and RET/PTC1 (R) 5'-CTTTCAGCATCTTCACGG-3'. The final PCR mixture contained 200 nM each of forward and reverse primers, 1x SYBR PCR mix (Applied Biosystems), and 5 µl cDNA template. Real-time PCR was performed with an ABI Prism 7700 sequence detector (Applied Biosystems). The thermal cycling conditions were: 2 min at 50 C, 10 min at 95 C, followed by 40 cycles of 30 sec at 95 C, 30 sec at 45 C, and 1 min at 72 C. Dissociation curves were recorded after each run, and the amplified products were visualized by 2% agarose gel electrophoresis.
Analysis of Pax-8 and TTF-1 by real-time quantitative RT-PCR (qRT-PCR)
The mRNAs for human Pax-8 and TTF-1 were measured by real-time qRT-PCR using ABI Prism 7700 Sequence detector, as previously described (16). The following primer/probe sequences were used: Pax-8 (F) 5'-GGACTACAAACGCCAGAACCC-3', Pax-8 (R) 5'-GGGCACAGTGTCATTGTCACA-3', PAX8 (Probe) 6FAM- ACCATGTTTGCCTGGGAGATCCGA-TAMRA, TTF-1 (F) 5'-ACGACTCCGTTCTCAGTGTCTG-3', TTF-1 (R) 5'-GCCCTCCA TGCCCACTTT-3', TTF1 (probe) 6FAM- CATCTTGAGTCCCCTGGAGGAAAGCTACA-TAMRA. The TaqMan probes were purchased from ABI, 5'-labeled with 6-carboxyfluorescein (FAM) and 3'-labeled with 6-carboxy-tetramethylrhodamine (TAMRA). Thermal cycling conditions were as follows, and reactions were performed in duplicate. Reverse transcription was performed at 48 C for 30 min, followed by activation of AmpliTaq Gold at 95 C for 10 min, followed by 40 cycles of amplification at 95 C for 15 sec and 60 C for 1 min. Each PCR run included the standard curve (10-fold serially diluted human normal thyroid RNA), test samples, no-template, and no reverse transcriptase controls. The standard curve was used to calculate the relative amounts of targets in test samples. Quantities of Pax-8 and TTF1 in test samples were normalized to the corresponding 18s rRNA (PE ABI, P/N 4308310). Genomic DNA was isolated from normal and papillary thyroid carcinoma tissue (Cooperative Human Tissue Network, National Cancer Institute, Bethesda, MD), or anaplastic thyroid carcinoma (ATC) tissue (Mayo Clinic, Jacksonville, FL). Tissue sample no. 9N was used for the standard curve.
| Results |
|---|
|
|
|---|
DNA profiling using STRs was performed on the 40 cell lines using 10 different genetic markers, including nine STR loci (D3S1358, vWA, FGA, D8S1179, D21S11, D18S51, D5S818, D13S317, and D7S820) and amelogenin. For these loci the probability of two samples selected at random having an identical profile is less than one in 1010. As shown in Table 2
, the PTC cell lines BHP5-16, BHP17-10, and BHP14-9 were identical to each other at all 10 loci. The BHP5-16, BHP17-10, and BHP14-9 were also identical to the PTC-derived NPA87 cell line at nine loci. The STR profile of the NPA87 cell line from the Haugen and Fagin laboratories was identical. BHP5-16, BHP17-10, and BHP14-9 are heterozygous at the D7S820 locus (8- and 10-repeat alleles), whereas the NPA87 cells have lost the 10-repeat allele (Table 2
). These results indicate that the BHP5-16, BHP17-10, and BHP14-9 cell lines are not only identical to each other, but also to the NPA87 cell line.
The ATCC database was searched with BHP5-16, BHP17-10, BHP14-9, and NPA87 profiles, and identified the MDA-MB-435S cell line as a match at all five loci tested in common (Table 2
; vWA, D5S818, D13S317, D7S820, plus amelogenin). Of note, the MDA-MB-435 breast cancer cell line has recently been identified as a derivative of the M14 melanoma cell line (17, 18). Because the original M14 cell line is not available through the ATCC, we obtained and analyzed the MDA-MB-435S genomic DNA from the ATCC and found it to have an identical STR profile to the aforementioned cancer cell lines at all loci except the D18S51, D5S818, and D3S1358 loci, where the MDA-MB-435S exhibited a loss of heterozygosity (LOH) (Table 2
). These profiles are consistent with the STR profile of an original stock of the MDA-MB-435 cell line from MD Anderson Cancer Center in addition to a cell stock from the Lombardi Comprehensive Cancer Center (Table 2
and Ref. 18). The M14 melanoma cell line was established in 1976 from a metastatic amelanotic lesion (19), predating the reported establishment of the NPA87 cell line and its derivatives. Because NPA87 is a commonly used thyroid cancer cell line, we further examined this using 250K SNP arrays. This showed that NPA87 is identical to BHP5-16, BHP17-10, BHP14-9, and MDA-MB-435S/M14 at 96.1% of the polymorphic sites, confirming the genetic identity of these cell lines.
The MDA-MB-435/M14 cell line has been reported to express a heterozygous BRAF V600E mutation (20, 21). In contrast, the NPA87 cells have been reported to express both a homo/hemizygous BRAF mutation (22, 23) and a heterozygous BRAF mutation (24). The NPA87 cell line and the BHP5-16, BHP17-10, and BHP14-9 cell lines in our laboratories (B.R.H. and J.A.F.) express a homo/hemizygous BRAF V600E mutation (22) (data not shown). We analyzed the BRAF V600E mutational status in the MDA-MB-435S cells from the ATCC and found that the MDA-MB-435S cell lines express a heterozygous BRAF V600E mutation (data not shown). This difference is likely due to a LOH that occurred during culturing (25).
The thyroid cancer derived KAT5, KAT10, KAT4, KAT7, KAT50, KAK1, ARO81-1, and MRO87-1 cell lines are identical to each other and to the colon cancer-derived HT-29 cell line
The PTC-derived KAT5 and KAT10 cells, the ATC-derived KAT4 and ARO81-1, and the follicular thyroid carcinoma (FTC)-derived MRO87-1 cells were analyzed by STR, and found to be genetically similar and likely identical. The KAT10 cells differed at the D18S51 locus, which in addition to a 13-repeat allele, also retained a 12-repeat allele (Table 2
). These results are consistent with recent studies, which found the KAT4, KAT10, and FTC-derived KAK1 cells to be identical by STR profiling (12, 13). In addition, the KAT4 STR profile shown here is identical to that published recently by Lee et al. (15), except at the D8S1179 locus there was a 16-repeat allele in addition to the 10-repeat allele that we observed here (Table 2
). These differences are likely due to clonal drift, and indicate that the KAT5, KAT10, and KAT4 cells are identical. We also found that the KAT5, KAT10, and KAT4 cell lines were not only identical to each other, but also to the MRO87-1 and ARO81-1 cell lines (Table 2
). The ARO81-1 cell line from the Haugen and Fagin laboratories showed identical STR profiles. By 250K SNP array analysis, we found that the KAT7 and KAK1 cell lines from the Fagin laboratory are identical to the ARO81-1 cell lines at 97.1% of the sites, and the KAT50 and ARO81-1 cell lines share identity at 97.4% of the sites, proving that they are derivatives from the same line.
By searching the ATCC database, we discovered that the STR profiles of these cell lines matched that of the HT-29 colon cancer cell line at those loci tested in common (Table 2
). Therefore, we obtained the HT-29 cell line from the ATCC and found that the STR profile of the HT-29 cell line is identical to that of the ARO81-1, MRO87-1, KAT5, KAT10, and KAT4 except at the D21S11 locus, where the HT-29 cell line retained the 29-repeat allele. The HT-29 cells also differed from the KAT10 cells at the D18S51 locus as discussed previously. By SNP 250K array genotyping, ARO81-1 and HT-29 were identical at 93.1% of the sites, confirming that they are genetically identical.
Consistent with previous studies, we found that the KAT5, KAT10, and KAT4 cells express a heterozygous BRAF mutation (24, 26, 27) (data not shown). Previous studies regarding BRAF mutational status in ARO81-1 cells have shown that these cells are homo/hemizygous (24) or heterozygous (23, 26, 27). The ARO81-1 cells in our laboratories (B.R.H. and J.A.F.) express a heterozygous BRAF mutation, consistent with that reported by two other laboratories (23, 27), and consistent with the ARO81-1 cell line matching the KAT5, KAT10, KAT4, KAK1, and KAT7 cell lines. The BRAF V600E mutation is also prevalent in colon cancer (28), and we found that the HT-29 colon cancer cell line from the ATCC expressed a heterozygous BRAF mutation, consistent with a previous report (29) and the genetic identity of these cell lines. The HT-29 colon cancer cell line was established in 1964 (30), before the establishment of the ARO81-1 cell line, which was established in 1981.
The ATC cell line DRO90-1 is genetically identical to the A-375 melanoma-derived cell line
The ATC-derived DRO90-1 cell line was profiled by STR analysis (Table 2
). This cell line was unique among our panel of thyroid cancer cell lines; however, the DRO90-1 matched the STR profile of the A-375 melanoma cell line (ATCC database) at all five overlapping loci. Therefore, we obtained the earliest available passage of the A-375 cell line (passage 10 from 1972), and we found that these two cell lines were identical at all 10 loci (Table 2
). We further confirmed that the A-375 and DRO90-1 cell lines both express a homo/hemizygous BRAF mutation and lack TTF-1 and Pax-8 expression (data not shown). These results indicate that the DRO90-1 cell line (established in 1990) is genetically identical to the A-375 melanoma cell line by STR profiling.
The PTC-derived BHP2-7, BHP10-3, BHP7-13, and TPC1 cells are genetically identical
Table 2
shows that the BHP2-7, BHP10-3, and BHP7-13 cells are not only identical to each other by STR analysis, but also to the TPC1 cell line, which was the original cell line isolated (31). We have confirmed the presence of the RET/PTC1 rearrangement in the BHP2-7, BHP10-3, BHP7-13, and TPC1 cell lines by SYBR green RT-PCR (Table 3
and data not shown). Expression levels of the thyroid-specific transcription factors, Pax-8 and TTF-1, were evaluated by qRT-PCR, and high levels of Pax-8 mRNA were observed by qRT-PCR (Table 3
; 407 pg mRNA/ng rRNA), further supporting the TPC1 cell line as a cell line derived from thyroid cancer. TTF-1 expression was not detected by qRT-PCR (Table 3
), consistent with previous studies (11).
|
The PTC-derived K1 and K2 cells lines have identical STR profiles as shown in Table 2
. The STR profiles obtained for the K1 and K2 cells in our study match the K1 STR profile provided by the ECACC (personal communication) for those loci tested in common (Table 2
). We have confirmed the presence of the BRAF V600E mutation in addition to high levels of Pax-8 mRNA by qRT-PCR (Table 3
; 627 pg mRNA/ng rRNA), which suggests that this cell line is likely of thyroid origin. TTF-1 mRNA was not detected, consistent with previous studies (11).
The PTC-derived cell lines, BCPAP and KTC-1 are unique
The PTC-derived BCPAP and KTC-1 cell line STR profiles are shown in Table 2
. The BCPAP cell line does not match any of the cell lines tested in this study or any other cell line listed in the ATCC, DSMZ, or ECACC databases. The BCPAP profile reported here is identical to that reported by the DSMZ for those loci tested in common (Table 2
) and also that reported by van Staveren et al. (12). Consistent with previous studies (11, 12), the BCPAP cells also express high levels of Pax-8 and TTF-1 mRNA, as determined by qRT-PCR (Table 3
). We have also confirmed the presence of the BRAF V600E mutation (homo/hemizygous) (Table 3
) (11). Likewise, the KTC-1 cells express high levels of Pax-8 and TTF-1, and express a heterozygous BRAF mutation, in agreement with previous studies (23, 32). The KTC-1 cell line does not match any other cell line listed in the ATCC, DSMZ, or ECACC databases.
The FTC-derived TT2609-CO2, FTC133, ML1, and WRO82-1 cell lines are unique
The STR profiles of the TT2609-CO2, FTC133, ML1, and WRO82-1 cell lines are shown in Table 2
. The TT2609-CO2 and ML1 STR profiles are consistent with that reported in the DSMZ database, and the FTC133 STR profile is consistent with that of the ECACC (personal communication) (Table 2
). The WRO82-1 STR profile from the Haugen laboratory is in agreement with a previously published profile (12), but not with that of the ECACC (personal communication) (data not shown). We found the TT2609-CO2, ML1, and FTC133 cells express wild-type (WT) BRAF, consistent with the follicular origin of these cells (Table 3
). However, the WRO82-1 cell line from the Haugen laboratory expresses mutant BRAF V600E. The WRO82-1 cell line from the Haugen laboratory expresses moderate levels of Pax-8 (126 pg mRNA/ng rRNA) and undetectable TTF-1 (Table 3
). In contrast, the TT2609-CO2, FTC133, and ML1 cell lines express high levels of TTF-1 (205-687 pg mRNA/ng rRNA), and low to undetectable levels of Pax-8 (0-76 pg mRNA/ng rRNA; Table 3
).
The ATC-derived cell lines 8505C, SW1736, Cal-62, T235, T238, Uhth-104, ACT-1, HTh74, KAT18, TTA1, FRO81-2, HTh7, C643, BHT101, and KTC-2 are unique
The STR profiles of the 8505C, SW1736, Cal-62, T235, T238, Uhth-104, ACT-1, HTh74, KAT18, TTA1, FRO81-2, HTh7, C643, BHT101, and KTC-2 cell lines are shown in Table 2
. These cells are unique among our panel of cell lines. The STR profiles of the 8505C, SW1736, Uhth-104, HTh74, HTh7, and C643 (Table 2
) are in agreement with previously published profiles, with the exception of the SW1736 cell line, which showed a difference at the D3S1358 locus, the Uhth-104 cell line, which differed at the D18S51 locus, and the C643 cell line, which differed at the D18S51 locus (12, 15). The STR profiles of the 8505C, Cal-62, and BHT101 cell lines are also consistent with their respective profiles in the DSMZ database. The STR profiles of the HTh7, HTh74, SW1736, and C643 cell lines, which were tested independently by the Fagin and Haugen laboratories, are in agreement with each other, with the exception of minor differences due to LOH. The Haugen and Fagin laboratory KAT18 STR profiles are identical. BRAF mutational status of these cell lines is in agreement with previous reports (if known), with exception of the FRO81-2 cell line from the Haugen laboratory, which express WT BRAF (Table 3
). The SW1736, Cal-62, T235, T238, Uhth-104, and ACT-1 cell lines expressed moderate (80-160 pg mRNA/ng rRNA) to high (230-658 pg mRNA/ng rRNA) levels of Pax-8; and the 8505C, Uhth-104, and ACT-1 cell lines expressed high (580 pg mRNA/ng rRNA) to moderate (79-114 pg mRNA/ng rRNA) levels of TTF-1 (Table 3
). The remaining cell lines expressed low to undetectable levels of Pax-8 and TTF-1.
| Discussion |
|---|
|
|
|---|
The expression levels of Pax-8 and TTF-1 in the unique cell lines reported here are in agreement with previous studies analyzing the expression of these genes in many of our cell lines, including the TPC1, BCPAP, and K1 cells (11). In other cases, including the C643, 8505C, HTh74, and WRO82-1 cell lines, we observed low to undetectable levels of Pax-8 and/or TTF-1, whereas previous studies reported expression of one or both of these transcription factors (11, 12). The reason for these differences is unclear but is likely due to different cell culture or RT-PCR conditions, as discussed below. Because the expression of tissue-specific genes is often reduced or lost in undifferentiated cancer, we also analyzed the expression of these genes in normal vs. thyroid tumor tissues using qRT-PCR. Interestingly, decreased levels of Pax-8 and TTF-1 mRNA are only observed in a subset of the tumor tissue samples compared with matched normal samples, indicating that these transcription factors are not necessarily down-regulated in papillary thyroid carcinoma (Table 3
). Pax-8 mRNA levels in the PTC-derived cell lines are similar to that observed in the PTC and ATC tissue samples, and are lower in the FTC- and ATC-derived cell lines. In contrast, TTF-1 expression is lower in all of the cell lines, compared with normal, PTC, and ATC tissue samples (Table 3
).
The results of our BRAF mutational analysis are also consistent with previous studies, with the exception of the FRO81-2 and WRO82-1 cell lines from the Haugen laboratory, which were found to express WT and mutant BRAF, respectively (Table 3
). The Haugen laboratory WRO82-1 STR profile matches that previously published by van Staveren et al. (12), but not that of the ECACC (personal communication) or the Fagin laboratory (data not shown). This suggests that two different WRO82-1 and FRO81-2 cell lines have been distributed based on STR profiling and BRAF mutation. Both cell lines were obtained directly from the source by each laboratory (B.R.H. and J.A.F.).
For the cell lines we believe to be genetically identical, we and others have reported distinct responses which can seem confusing. For example, we have shown that the BHP5-16, BHP17-10, and BHP14-9 cell lines are more sensitive to retinoid and rexinoid treatment than the NPA87 cell line (16), even though we have shown that they are genetically identical by STR profiling. We believe that these differences in response are due to clonal drift, which can occur when cells are maintained in culture for long periods (33, 34, 35). In contrast, we and others have found that many of the sublines identified in this study, including the BHP2-7, BHP7-13, and BHP10-3 (16, 36), and the KAK1, KAT5, KAT7, KAT10, and ARO81-1 (37, 38) indeed respond as expected for identical cell lines, suggesting that long-term culture does not necessarily result in phenotypical drift of all cell lines.
Cell line cross-contamination and misidentification are long-standing problems, and redundancy among thyroid cancer cell lines has been previously reported (KAT5, KAT10, KAT4, KAK1; and FB2, TPC1) (11, 12, 13). After initial submission of this manuscript, the DSMZ (www.dsmz.de/) reported that the ONCO-DG-1 PTC-derived cell line (39) is actually a derivative of the OVCAR-3 ovarian cancer cell line. We have also confirmed that the ONCO-DG-1 and OVCAR-3 cell lines are identical by STR profiling (data not shown).
STR profiling is commonly used to determine genetic identity; however, STR results can be complicated by clonal drift that occurs in culture. Indeed, we observed minor differences in cell lines we believe to be identical. For example, the NPA87 cell line differed from the BHP5-16, BHP17-10, BHP14-9, and MDA-MB-435S cell lines at the D7S820 locus, where the BHP cell lines potentially gained a 10-repeat allele or alternatively, the NPA87 cell line lost the 10-repeat allele (Table 2
). These cell lines also differed at the D3S1358, D18S51, and D5S818 loci by a loss or gain of an allele (Table 2
). Previous studies evaluating the STR profile of two different stocks of the MDA-MB-435 cell line from the MD Anderson Cancer Center and the Lombardi Cancer Center Tissue Shared Resource also exhibit a small degree of clonal drift, which is comparable to our results (Table 2
and Ref. 18). Despite these differences in STR profiles due to a loss or gain in alleles, in no case did we find differences in the number of allele repeats in the cell lines we believe to be identical. Finally, none of the differences in what we believe to be identical cell lines reduced the similarity between cell line profiles to less than the 80% cutoff proposed by Masters et al. (3). Furthermore, by 250K SNP array analysis, the NPA87 cell line and its derivates matched the MDA-MB-435S/M14 cell line at 96.1% of the polymorphic sites, providing further evidence that these cell lines are identical.
For the cell lines that we believe to be identical, we found that these cell lines harbor identical BRAF mutations, with the exception of the NPA87, BHP5-16, BHP17-10, and BHP14-9 cell lines, which are homo/hemizygous (23) (data not shown), whereas the MDA-MB-435S/M14 cell line is BRAF heterozygous (20, 21) (data not shown). The precise reason for this difference is not known but could be due to LOH through loss of the WT allele or to uniparental disomy, either of which can occur during long-term culturing, as discussed previously, and has been shown for the BRAF WT allele at the T1799A site in other cell lines (40, 41). We also discovered that the Haugen laboratory MRO87-1 cell line expresses a heterozygous BRAF mutation, which is not consistent with the FTC-origin of this cell line but agrees with the genetic identity with the ARO81-1, KAT5, KAT10, KAT4, and HT-29 cell lines by STR profiling (Table 2
). Previous studies have shown that the MRO87-1 cell line expresses either WT BRAF or mutant BRAF V600E (24, 37), suggesting that there are two different MRO87-1 cell lines in use. Again, this cell line was obtained directly from the laboratory that generated MRO87-1.
For the cell lines we believe to be of nonthyroid origin, previous studies have shown expression of thyroid specific genes, including Pax-8 and TTF-1 in the KAT-4, KAT-10, and KAK-1 cell lines (12, 42, 43). In addition, increased expression of thyroid-specific genes after pharmacological manipulation has been shown in some of the cell lines that we believe are of nonthyroid origin, including NPA87, DRO90-1, and ARO81-1 (37, 44, 45). In most cases expression of thyroid-specific genes was detected using PCR-based assays, and it is conceivable that trace amounts of these gene products may be expressed under certain conditions in cells that are not of thyroid origin. Alternatively, RT-PCR based methods can result in the amplification of nonspecific products. Indeed, we detected moderate levels of TTF-1 mRNA in some of the cell lines we believe to be of nonthyroid origin, including the ARO81-1 cell line, using SYBR green RT-PCR (data not shown). However, using this approach we noted variable melting temperatures between different cell lines (75-82C, data not shown) and, therefore, used TaqMan qRT-PCR to measure more accurately TTF-1 and Pax-8 mRNA levels (Table 3
). Using this approach we were able to specifically measure levels of TTF-1 and Pax-8 mRNA in our cell lines and did not detect expression of either of these transcription factors in the cell lines we believe to be of nonthyroid origin, including the ARO81-1 and NPA87 cells (data not shown). Lack of expression of thyroid-specific genes has been attributed to a dedifferentiated phenotype. Given the results of the current study, it is possible that the low expression of thyroid-specific genes is due to these cell lines being of nonthyroid origin.
Although the STR profiling approach used in this study represents an accurate and powerful approach for determining cell line identity, STR profiling does not provide any information regarding tissue of origin. Because original patient material from which most of these misidentified cell lines were derived is not available, we cannot definitely identify the tissue of origin for any of these cell lines based on STR or SNP profiling. However, given the time line of establishment and lack of Pax-8 and TTF-1 expression, we believe that the NPA87, ARO81-1, MRO87-1, and DRO90-1 cell lines and their respective sublines were established by cross-contamination at the source. The "90," "81," and "87" in the names DRO90-1, ARO81-1, MRO87-1, and NPA87 refer to the year the cell lines were established (Canlapan, D., University of California Los Angeles, Los Angeles, CA, personal communication). Given that the first publication of the A-375 cell line was in 1973 (46), we believe that the DRO90-1 is a subline of A-375 that resulted from cell line contamination.
Likewise, the M14 melanoma cell line was established in 1976 (19), and the MDA-MB-435 misidentification with the M14 cell line occurred before 1982 (17), both of which predate the establishment of the NPA87 cell line. Of note, the M14 cell line was also established at University of California Los Angeles (originally called UCLA-SO-M14) and was used by the originating investigators of the NPA87 cell line before the establishment of the NPA87 cell line (47). Thus, it is likely that the establishment of the NPA87 cell line is also due to cross- contamination.
The HT-29 colon cancer cell line was established in 1964 (30) and deposited into the ATCC in 1981 (48, 49). Previous studies have shown that the 1981 HT-29 stock deposited in the ATCC is genetically identical to an original HT-29 stock from 1967 (49), indicating that the ARO81-1 cell line and the MRO87-1 cell line from the Haugen laboratory were likely established by cross-contamination. The NPA87, ARO81-1, DRO90-1, and MRO87-1 cell lines and their respective sublines have been widely used as thyroid cancer models for the past 20 yr.
In conclusion, we have identified significant redundancy and misidentification among 20 of 40 thyroid cancer cell lines tested (50%). We found that Pax-8 and/or TTF-1 is expressed in the majority (15 of 23) of the unique cell lines. For cell lines that do not express either thyroid-specific gene, all are of anaplastic origin, suggesting that the loss of these genes could be due to the dedifferentiated phenotype. Based on the studies presented here, we believe the cell lines listed in Table 3
represent viable models to study thyroid cancer and at the time of this report, are unique.
| Acknowledgments |
|---|
| Footnotes |
|---|
Disclosure Statement: The authors have nothing to disclose.
First Published Online August 19, 2008
Abbreviations: ATC, Anaplastic thyroid carcinoma; ATCC, American Type Culture Collection; DSMZ, German Collection of Microorganisms and Cell Cultures; ECACC, European Collection for Cell Cultures; FAM, 6-carboxyfluorescein; FTC, follicular thyroid carcinoma; LOH, loss of heterozygosity; PTC, papillary thyroid cancer; qRT-PCR, quantitative RT-PCR; SNP, single nucleotide polymorphism; STR, short tandem repeat; TAMRA, 6-carboxy-tetramethylrhodamine; WT, wild type.
Received May 22, 2008.
Accepted August 12, 2008.
| References |
|---|
|
|
|---|
receptors in human and rat thyroid cells and regulation of the receptors by thyrotropin. Endocrinology 125:1783–1788
activation induces cell cycle arrest via the p53-independent pathway in human anaplastic thyroid cancer cells. Jpn J Cancer Res 93:1358–1365[CrossRef]This article has been cited by other articles:
![]() |
E. Diallo-Krou, J. Yu, L. A. Colby, K. Inoki, J. E. Wilkinson, D. G. Thomas, T. J. Giordano, and R. J. Koenig Paired Box Gene 8-Peroxisome Proliferator-Activated Receptor-{gamma} Fusion Protein and Loss of Phosphatase and Tensin Homolog Synergistically Cause Thyroid Hyperplasia in Transgenic Mice Endocrinology, November 1, 2009; 150(11): 5181 - 5190. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Jin, T. Jiang, D. M. Rosen, B. D. Nelkin, and D. W. Ball Dual Inhibition of Mitogen-Activated Protein Kinase Kinase and Mammalian Target of Rapamycin in Differentiated and Anaplastic Thyroid Cancer J. Clin. Endocrinol. Metab., October 1, 2009; 94(10): 4107 - 4112. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. J. Robinson-White, I. Bossis, H.-P. Hsiao, M. Nesterova, W. W. Leitner, and C. A. Stratakis 8-Cl-Adenosine Inhibits Proliferation and Causes Apoptosis in B-Lymphocytes via Protein Kinase A-Dependent and Independent Effects: Implications for Treatment of Carney Complex-Associated Tumors J. Clin. Endocrinol. Metab., October 1, 2009; 94(10): 4061 - 4069. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Liu, P. Hou, Z. Liu, G. Wu, and M. Xing Genetic Alterations in the Phosphoinositide 3-Kinase/Akt Signaling Pathway Confer Sensitivity of Thyroid Cancer Cells to Therapeutic Targeting of Akt and Mammalian Target of Rapamycin Cancer Res., September 15, 2009; 69(18): 7311 - 7319. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. Prante, S. Maschauer, V. Fremont, J. Reinfelder, R. Stoehr, M. Szkudlinski, B. Weintraub, A. Hartmann, and T. Kuwert Regulation of Uptake of 18F-FDG by a Follicular Human Thyroid Cancer Cell Line with Mutation-Activated K-Ras J. Nucl. Med., August 1, 2009; 50(8): 1364 - 1370. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. DEEN, O. L. GRIFFITH, H. MASOUDI, A. GOWN, S. J.M. JONES, and S. M. WISEMAN Anaplastic Thyroid Carcinoma Exhibits Intratumoral Molecular Homogeneity for a Therapeutic Target Panel Anticancer Res, July 1, 2009; 29(7): 2437 - 2444. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Rivas, B. Mellstrom, B. Torres, G. Cali, A. M. Ferrara, D. Terracciano, M. Zannini, G. Morreale de Escobar, and J. R. Naranjo The DREAM Protein Is Associated with Thyroid Enlargement and Nodular Development Mol. Endocrinol., June 1, 2009; 23(6): 862 - 870. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. C. Ricarte-Filho, M. Ryder, D. A. Chitale, M. Rivera, A. Heguy, M. Ladanyi, M. Janakiraman, D. Solit, J. A. Knauf, R. M. Tuttle, et al. Mutational Profile of Advanced Primary and Metastatic Radioactive Iodine-Refractory Thyroid Cancers Reveals Distinct Pathogenetic Roles for BRAF, PIK3CA, and AKT1 Cancer Res., June 1, 2009; 69(11): 4885 - 4893. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. E. Schweppe, A. A. Kerege, J. D. French, V. Sharma, R. L. Grzywa, and B. R. Haugen Inhibition of Src with AZD0530 Reveals the Src-Focal Adhesion Kinase Complex as a Novel Therapeutic Target in Papillary and Anaplastic Thyroid Cancer J. Clin. Endocrinol. Metab., June 1, 2009; 94(6): 2199 - 2203. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Nucera, J. Eeckhoute, S. Finn, J. S. Carroll, A. H. Ligon, C. Priolo, G. Fadda, M. Toner, O. Sheils, M. Attard, et al. FOXA1 Is a Potential Oncogene in Anaplastic Thyroid Carcinoma Clin. Cancer Res., June 1, 2009; 15(11): 3680 - 3689. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. C Smallridge, L. A Marlow, and J. A Copland Anaplastic thyroid cancer: molecular pathogenesis and emerging therapies Endocr. Relat. Cancer, March 1, 2009; 16(1): 17 - 44. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. D. Castellone, V. De Falco, D. M. Rao, R. Bellelli, M. Muthu, F. Basolo, A. Fusco, J. S. Gutkind, and M. Santoro The {beta}-Catenin Axis Integrates Multiple Signals Downstream from RET/Papillary Thyroid Carcinoma Leading to Cell Proliferation Cancer Res., March 1, 2009; 69(5): 1867 - 1876. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. A. Marlow, L. A. Reynolds, A. S. Cleland, S. J. Cooper, M. L. Gumz, S. Kurakata, K. Fujiwara, Y. Zhang, T. Sebo, C. Grant, et al. Reactivation of Suppressed RhoB is a Critical Step for the Inhibition of Anaplastic Thyroid Cancer Growth Cancer Res., February 15, 2009; 69(4): 1536 - 1544. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Endocrinology | Endocrine Reviews | J. Clin. End. & Metab. |
| Molecular Endocrinology | Recent Prog. Horm. Res. | All Endocrine Journals |