Next Article in Journal
Lymph Node Involvement in Recurrent Serous Borderline Ovarian Tumors: Current Evidence, Controversies, and a Review of the Literature
Previous Article in Journal
Early Onset Colorectal Cancer in Arabs, Are We Dealing with a Distinct Disease?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 15
Issue 3
10.3390/cancers15030879
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of Novel Murine BRAFV600E-Driven Papillary Thyroid Cancer Cell Lines for Modeling of Disease Progression and Preclinical Evaluation of Therapeutics

by
Grace Purvis Branigan
1,†,
Victoria Casado-Medrano
1,†,
Alison B. O’Neill
1,
Julio C. Ricarte-Filho
1,
Nicole Massoll
2,
Madeleine Salwen
1,
Zachary Spangler
1,
Michele Scheerer
1,
Edward K. Williamson
3,
Andrew J. Bauer
1 and
Aime T. Franco
1,4,*
1
Division of Endocrinology and Diabetes, Children’s Hospital of Philadelphia, Philadelphia, PA 19104, USA
2
Department of Pathology, Winthrop P. Rockefeller Cancer Institute, University of Arkansas for Medical Sciences, Little Rock, AR 72205, USA
3
Department of Pathology and Laboratory Medicine, Children’s Hospital of Philadelphia, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA
4
Abramson Cancer Center, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Cancers 2023, 15(3), 879; https://doi.org/10.3390/cancers15030879
Submission received: 4 October 2022 / Revised: 18 November 2022 / Accepted: 24 November 2022 / Published: 31 January 2023
(This article belongs to the Section Methods and Technologies Development)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Laboratory experimental models are essential for advancements in thyroid cancer translational research. In this study, we discuss the development and characterization of six cell-based models of BRAFV600E-driven papillary thyroid cancer that closely mimic the heterogeneous pathological progression of the disease seen in patients despite sharing a common driver mutation. We anticipate that these newly developed models will allow for the discovery of novel mechanisms that determine differences in disease progression among BRAFV600E-driven thyroid cancers and facilitate the testing of therapeutic interventions.

Abstract

The Cancer Genome Atlas study in thyroid cancer exposed the genomic landscape of ~500 PTCs and revealed BRAFV600E-mutant tumors as having different prognosis, contrasting indolent cases and those with more invasive disease. Here, we describe the generation and characterization of six novel BRAFV600E-driven papillary thyroid cancer (PTC) cell lines established from a BrafV600E+/−/Pten+/−/TPO-Cre mouse model that spontaneously develop thyroid tumors. The novel cell lines were obtained from animals representing a range of developmental stages and both sexes, with the goal of establishing a heterogeneous panel of PTC cell lines sharing a common driver mutation. These cell lines recapitulate the genetics and diverse histopathological features of BRAFV600E-driven PTC, exhibiting differing degrees of growth, differentiation, and invasive potential that may help define mechanisms of pathogenesis underlying the heterogeneity present in the patient population. We demonstrate that these cell lines can be used for a variety of in vitro applications and can maintain the potential for in vivo transplantation into immunocompetent hosts. We believe that these novel cell lines will provide powerful tools for investigating the molecular basis of thyroid cancer progression and will lead to the development of more personalized diagnostic and treatment strategies for BRAFV600E-driven PTC.

1. Introduction

Thyroid cancer is the most common endocrine malignancy and the most frequently diagnosed cancer among adolescents and young adults between ages 15–29 [1,2,3]. Papillary thyroid cancer (PTC) is the most prevalent form of thyroid cancer and accounts for more than 90% of all thyroid cancer diagnoses (SEER Cancer Statistics Review, 1975–2017, National Cancer Institute). While well-differentiated (WD) PTC has an excellent prognosis, with a 5-year survival rate upwards of 90%, a small percentage of these tumors are believed to progress to poorly differentiated (PDTC) or anaplastic thyroid cancer (ATC). PDTC and ATC have a far worse prognosis, as they are often refractory to treatment and account for the majority of thyroid cancer-related mortality [4,5,6,7]. BRAFV600E is the most common driver mutation of PTC in adults and the second most common driver mutation in children [8,9,10]. The BRAFV600E mutation confers constitutive activation of the kinase, leading to sustained activation of the mitogen-activated protein kinase (MAPK) signaling pathway [11,12]. Clinical studies have associated the BRAFV600E mutation with extrathyroidal extension, lymphatic and vascular invasion, local lymph node metastasis, higher clinical stage, persistent disease, reduced avidity for radioactive iodine, and lower survival rates in adults with PTC [10,11,12,13,14,15,16,17,18]. In contrast, the BRAF mutation has not been associated with the most invasive and metastatic disease in children with PTC; therefore, more data are needed to understand the natural history of BRAF-driven PTC in this population [8,19,20,21]. Notably, genetic alterations that trigger the activation of both the MAPK and phosphoinositide 3-kinase (PI3K)/Akt pathway become increasingly prevalent as disease progresses in adults. Moreover, it has been suggested that progression from PTC to ATC may be facilitated by mutations leading to the activation of the PI3K/Akt pathway co-occurring with BRAF mutations due to the increased mutation frequency in higher grade tumors [22,23,24,25].
Most PTCs are successfully managed with standard multimodal therapy. Patients undergo surgical tumor resection and lymph node dissection for localized disease, often followed by radioactive iodine (RAI) ablation, as adjuvant therapy in patients with high risk of recurrence or metastatic disease. Patients are then placed on thyroid hormone replacement therapy, with or without long-term thyroid hormone suppression, based on risk of recurrence. In the subset of patients with advanced and progressive RAI-refractory disease, oncogene-targeted therapies represent essential advancements in care and continue to be important areas of preclinical and clinical investigation in both adult and pediatric thyroid cancer.
It remains imperative to develop tumor-specific in vitro and in vivo preclinical models that recapitulate patient disease to prospectively evaluate treatments for efficacy and to decipher molecular mechanisms that confer therapeutic resistance. Although cell lines derived from human thyroid tumors provide a powerful tool to study novel therapies, research using human cells is limited to the tissue culture environment or xenograft models in immunocompromised hosts, which fails to recapitulate the interactions between tumor cells and the immune system in the host. Further, it was discovered that many previously reported thyroid tumor cell lines were either not unique or not actually of thyroid origin, further limiting the utility of some of these models [26]. In contrast, cell lines established from mouse models of thyroid cancer can be transplanted into immunocompetent syngeneic hosts, providing researchers with the ability to study drug responses in the physiologic tumor microenvironment. However, these cell lines do not always represent the heterogeneity observed in patient samples, and many of the previously reported lines are representative of advanced and anaplastic thyroid cancers [27,28]. We have previously reported the derivation of novel follicular thyroid cancer cell lines from HrasG12V-driven models of FTC [29] and sought here to develop cell lines derived from BRAFV600E-driven PTCs.
In this study, we detail the establishment of six independent cell lines from the BrafV600E+/−/Pten+/−/TPO-Cre mouse model of PTC. We have previously described a mouse model of PTC whereby thyroid-specific expression of BrafV600E induces development of thyroid tumors with 100% penetrance by five weeks of age [30]. Additionally, we have reported a mouse model of PTC utilizing the same thyroid-specific expression of BrafV600E and homozygous Pten inactivation in mice to achieve concomitant MAPK and PI3K/Akt pathway activation. These mice develop PTCs that rapidly progress to PDTC with 100% penetrance and lethality by weaning [31]. These models recapitulate the genetics and histopathological features of PTC as well as progress to a poorly differentiated state [30,31]. Here, we describe the establishment and characterization of six independent cell lines from murine thyroid tumors expressing BrafV600E and heterozygous Pten. Three cell lines were isolated from well differentiated PTC tumors (BB19, BB57, and BB342) and the other three cell lines were isolated from higher grade tumors with more invasive characteristics (BB1845, B1865, and B1866). This variability reproduces the clinical heterogeneity within the PTC patient population and, thus, allows for evaluation of characteristic features, in a subset of cell lines, that may be critical for disease progression. These cell lines represent novel and physiologically relevant research tools that can be used to evaluate therapeutic strategies and to illuminate factors that impact the progression of disease.

2. Materials and Methods

2.1. Derivation of Murine Thyroid Tumor Cell Lines and Wild-Type Thyrocyte Cultures

BB19, BB57, BB342, and BB1845 tumor cell lines were generated from thyroid tumors of BrafV600E+/−/Pten+/−/TPO-Cre mice of a pure 129/SvImJ genetic background (Jackson Laboratory mouse strain #:002448 here on referred to as 129SvJ). B1865 and B1866 tumor cell lines were established from thyroid tumors of BrafV600E+/−/Pten+/−/TPO-Cre mice of a mixed genetic background. Thyroid tumors or WT thyroids were dissected and minced. The pieces were then resuspended in Ham’s F12 medium (Corning, Glendale, AZ, USA) supplemented with 10% fetal bovine serum (FBS, Gibco), 2 mM L-Glutamine (Gibco), Penicillin/Streptomycin/Amphotericin B (Millipore-Sigma, Burlington, MA, USA), 10 µg/mL Piperacillin sodium salt (Millipore-Sigma, Burlington, MA, USA), and 10 µg/mL Ciprofloxacin Hydrochloride (Corning, Glendale, AZ, USA). Cells were then seeded into bind tissue culture flasks (Corning, Glendale, AZ, USA), and each flask was given an additional 500 μL of FBS. They were maintained at 37 °C in 5% CO2. After 24–48 h, the media were replaced, and the cells were maintained in 10% FBS F12 media supplemented with antibiotic for the next two weeks. WT thyrocytes were only cultured for at most one week before RNA was extracted to limit outgrowth of non-thyrocyte cells in culture. Then, antibiotic-free complete Ham’s 12 media was used. To rule out stromal cell contamination, all tumor cell lines were passaged at least 5 times after seeding and then genotyped using primers specific for BrafV600E to check for complete recombination (Figure S2).

2.2. H&E and Trichrome Staining

Formalin fixed tissue was submitted to the Children’s Hospital of Philadelphia Pathology Core Laboratory for processing and embedding. All sectioning and staining were performed by the CHOP Pathology Core Laboratory. Trichrome histology slides were also scanned by the CHOP pathology core. H&E slides were evaluated by a pathologist blinded to experimental conditions.

2.3. RT-PCR Analysis

Total RNA from all cell lines was extracted using the Direct-zol RNA Miniprep (Zymo Research, Irvine, CA, USA). Equal amounts of RNA template were reverse transcribed using the Verso cDNA synthesis kit (Thermo Scientific, Waltham, MA, USA). The differential mRNA expressions of Tg, Tshr, Nkx2-1 (Ttf1), Pax8, Slc5a5 (Nis), and B2m were measured using pre-designed probes (Integrated DNA Technologies, Coralville, IA, USA) and Taqman master mix (Thermo Scientific, Waltham, MA, USA). A total of four µL of cDNA from tumor samples and independent passages of each cell line and WT thyroids were run in triplicate on the QuantStudio 3 (Thermo Fisher, Waltham, MA, USA). B2m was used as an internal control to normalize the expression levels of thyroid specific genes. Data analysis was based on the delta-delta Ct method.

2.4. Proliferation and Viability Assays

Cells were plated, 300 per well, in quadruplicate into 96-well plates. At designated time points, cell growth was assessed using the Cell Titer-Glo luminescent cell viability assay (Promega, Madison, WI, USA). Luminescence was determined using the Varioskan Lux (Thermo Scientific, Waltham, MA, USA). Day 0 represents 18 h after plating. Inhibitor studies were conducted using Dabrafenib methanesulfonate salt (D-5699, LC Laboratories, Woburn, MA, USA) and Trametinib (T-8123, LC Laboratories, Woburn, MA, USA) dissolved in DMSO. To assess the relative viability of cells in the presence of single inhibitors, Dabrafenib, Trametinib, or vehicle (DMSO) were added 18h after plating. Relative viability was assessed using the Cell Titer-Glo luminescent cell viability assay, and concentrations of Dabrafenib 10 nM and Trametinib 1 nM were selected for viability assays of the complete cell line panel.

2.5. Western Blot Analysis

To prepare total cellular protein lysates, 200,000 cells per well were seeded in six-well plates and allowed to attach for 24 h in complete medium. Cells were serum starved overnight before any specific treatment or stimulation. Then, cells were lysate in equal volumes of radioimmunoprecipitation (RIPA) buffer (50 mM Tris-HCl pH 7.4, 1% NP-40, 0.25% Na-deoxycholate, 150 mM NaCl, 1mM EDTA) supplemented with phosphatase and protease inhibitors (Thermo Scientific, Waltham, MA, USA) and 4× sample buffer (BioRad, Hercules, CA, USA) containing 200 mM DTT and boiled for 5 min. A small aliquot was taken after harvesting to measure protein concentration by BCA. To perform the SDS–polyacrylamide gel electrophoresis (PAGE), equal quantities of total cellular proteins from each sample were loaded into 10% acrylamide gels containing SDS-PAGE. Protein sizes were determined by comparison to a set of standards (Thermo Fisher, Waltham, MA, USA). To perform the immunoblot, proteins were transferred to a PVDF membrane pore size 0.45 μm (Millipore-Sigma, Burlington, MA, USA) after methanol activation. Primary antibodies (1:1000) specific to pERK (9101S), ERK (9102S), pAKT (9271S), AKT (4685S), PTEN (9188T), E-cadherin (14472s), Vimentin (5741S), and β-actin (3700S) (Cell Signaling Technologies (Danvers, MA, USA) were diluted in 5% bovine serum albumin (BSA) and 0.1% Tween-20 tris buffer (TBS-T) and incubated overnight at 4 °C. Specific binding of antibodies was detected with appropriate HRP-coupled secondary antibodies anti-rabbit (1:3000) or anti-mouse (1:5000) and visualized using enhanced chemiluminescence reagent (Thermo Fisher, Waltham, MA, USA). Protein expression levels were analyzed and compared by densitometry using iBright Analysis Software (Thermo Fisher, Waltham, MA, USA). Relative protein phosphorylation levels were obtained by normalizing the amount of phosphorylated protein to that of total protein. Original, uncropped Western blot membrane can be found in Figure S5.

2.6. Organoids

Organoids were grown using previously described methods [32]. Briefly, to analyze the organoid formation rate in the presence of inhibitors, Dabrafenib, Trametinib, or vehicle (DMSO) were added to the media 18 h after plating. Organoids were grown for 10 days before imaging. Organoid size was measured using ImageJ particle analysis tool within each 3220 × 2415 micron area per image.

2.7. Migration Assay

Migration assays were performed in 24-well plates with cell culture inserts with 8.0 μm pores (Falcon, Chicago, IL, USA); 75,000 tumor cells were seeded in the top of each insert in serum free Hams F-12 media. Complete 10% FBS Hams F12 media was added to the bottom chamber of the plate. Cells were then incubated at 37 °C for 18 h. The cells were then fixed using 4% Paraformaldehyde, and the cells were removed from the top of the chamber using cotton tipped applicators. The underside of the chambers was then stained with crystal violet. Images were taken on the EVOS M7000 imaging system (Invitrogen by Thermo Fisher Scientific).

2.8. Imaging

Phase contrast imaging was performed using the EVOS M7000 imaging system (Invitrogen by Thermo Fisher Scientific). H&E histology slides were imaged using the EVOS M7000 imaging system (Invitrogen by Thermo Fisher Scientific). Trichrome histology slides were scanned by the CHOP pathology core and images captured using Image Scope. For differential interference contrast (DIC) imaging live cells were imaged with a Zeiss Axiovert 200M platform using a 40× EC Plan Neo Fluar, 1.3 N.A., or a 63× Plan Apochromat, 1.4 N.A. objective (Carl Zeiss, Jena, Germany) equipped with an Orca ER CCD camera (Hamamatsu, Bridgewater, NJ, USA). Images were captured using Slidebook 6 imaging software (Intelligent Imaging Innovations, Denver, CO, USA) using an Orca ER CCD camera (Hamamatsu).

2.9. In Vivo Tumorigenicity Assay

All animal experiments were performed at the Children’s Hospital of Philadelphia and approved by the IACUC. 5.0 × 105 tumor cells (BB19, BB57, BB342, and BB1845) were pelleted by centrifugation, resuspended in 100 µL of a 1:1 mixture of Hams F-12 and Matrigel (Corning, Corning, NY, USA), and injected subcutaneously into the right hind flanks of 4–5-week-old male and female wild-type 129SvJ recipient mice. Tumor development was monitored weekly for 15 weeks. Mice were sacrificed due to tumor burden or loss of skin integrity as needed during that period or were collected finally at the end of the 15 weeks. Tumors were collected and fixed in 10% formalin for histological analysis.

2.10. Statistical Analysis

All data were analyzed using Prism 8 software (GraphPad, San Diego, CA, USA). Differences with p-values of ≤0.05 were considered statistically significant. All data comparing more than two groups were analyzed via ANOVA analysis with post hoc analysis (indicated in figure legends).

3. Results

3.1. Generation of Stable Tumor Cell Lines from BrafV600E Murine Thyroid Tumors

Papillary thyroid carcinomas developed in BrafV600E+/−/Pten+/−/TPO-Cre mice of a pure 129SvJ genetic background (cell lines BB19, BB57, BB342, and BB1845) or a mixed genetic background (cell lines B1865 and B1866). Tumors were collected from different sex animals and at a variety of time points from 3–24 weeks of age (Table 1). The majority of the tumor was digested to derive the cell lines with approximately 1/4 preserved for histopathologic analysis. Hematoxylin and eosin staining of the founding tumors for each cell line confirmed that the model exhibits PTC-like histologic features including papillae, nuclear grooves, psammoma bodies, fine chromatin (Figure 1A), and lymphocytic infiltration. BB57 exhibited classic PTC features such as fine chromatin, nuclear grooves, and psammoma bodies (Figure 1B), as well as invasion into the muscle (Figure 1C). In B1866, extra-capsular tumor extension was observed in association with the identified intra-thyroidal tumor (Figure 1D). Additionally, this tumor exhibits moderately differentiated PTC histology that maintains papillae and grooves but displays localized areas of necrosis (Figure 1E). Six independent cell lines were established from these six founding tumors, sharing the same genotype but derived from animals that varied in age and sex (Table 1). The lungs from the six founding mice were each examined for gross metastasis at time of dissection. No gross metastasis was observed at time of dissection, but suspicious areas were noted in the lungs of BB1845. Not metastasis could be confirmed by microscopic evaluation. Despite arising from tumors with the same driver mutation, the cells show heterogeneous morphologies among the six independent cell lines (Figure 1F). Heterogeneity in morphology and cell size is also observed within each cell line.
BrafV600E+/−/Pten+/−/TPO-Cre mice undergo Braf recombination specific to thyroid epithelium via the action of TPO-Cre. To confirm the presence of BrafV600E recombination in established cell lines, PCR analysis was performed using primers specific for Braf WT, BrafV600E recombined, and BrafV600E unrecombined alleles. As expected, each cell line exhibited one Braf WT band and one BrafV600E recombined band (Figure S1). In comparison, DNA isolated from non-target tissues from the tumor-bearing animals from which each the cell line was derived exhibited one Braf WT band and one BrafV600E unrecombined band, confirming that BrafV600E recombination was limited to thyroid epithelial cells. The absence of a BrafV600E unrecombined band in the six cell lines derived from these tumors confirms that the cell lines contained only thyroid epithelial cells, despite the morphologic heterogeneity of the cell lines.

3.2. Expression of Thyroid-Specific Genes, Proliferation Rate, and Activation of MAPK and PI3K Signaling Pathways in Novel BrafV600E Murine Cell Lines

The thyroid epithelial lineage and differentiation status of the cell lines was determined by measuring the mRNA levels of thyroid-specific genes including Tg, Tshr, Nkx2-1 (Ttf1), Pax8, and Slc5a5 (Sodium Iodide symporter: Nis) (Figure 2A). Pax8 and Ttf1 are transcription factors crucial to thyroid organogenesis, and their absence results in varying degrees of thyroid dysgenesis [33,34]. The expression of Ttf1 was similar among the six cell lines and was higher when compared to the expression in WT thyrocytes. Pax8 expression was lower than in WT thyrocytes in most cell lines, except BB19 an BB57, which had expression similar to that observed in WT thyrocytes. Pax8 and Ttf1 regulate the expression of the thyroid-specific genes Tg, Tshr, and Slc5a5 (Nis) [34]. Tg and Tshr expression levels were similar among the six cell lines but were very low compared to expression in WT thyrocytes. In contrast, Nis expression was variable between the six cell lines but was consistently much lower compared to WT thyrocytes. Like thyroid transcription factors, loss of Tg, Tshr, or Nis expression has been associated with PTC disease progression [35,36,37,38]. Loss of Nis expression is particularly relevant clinically as it is associated with disease refractory to RAI ablation therapy [39,40,41].
We next sought to evaluate the proliferation rates of the cell lines in 2D culture using the Cell Titer Glo luminescence-based viability assay (Figure 2B). The Cell Titer Glo method is based on the quantitation of ATP to determine the number of metabolically active cells. Most of the cell lines had a lag-phase of approximately 48 h. At day three, BB1845, B1866, and BB19 appeared to enter logarithmic growth, whereas BB57 entered logarithmic growth at day 4. B1865 and BB342 did not enter logarithmic growth phases within 4 days. Maximal differences were observed at day 4 clearly showing heterogeneity in proliferation rates among the six cell lines. BB1845 was the most proliferative cell line, while B1865 and BB342 had the longest lag phases.
Next, Western blot was used to assess the MAPK and PI3K pathway activation in response to acute FBS stimulation and to quantify their basal levels after growth factor deprivation (Figure 2C). The cell lines displayed heterogeneous basal levels of ERK, pERK, Akt, and pAKT, with little increase in pERK and pAKT in response to FBS stimulation (Figure 2C,D). As described previously, all of the cell lines were derived from mice harboring tumors with thyrocyte-specific heterozygous deletion of Pten. To determine whether levels of PTEN protein expression correlated with basal pAKT levels, Western blot analysis was performed for PTEN. Consistent with their elevated basal pAKT expression, BB57, BB1845, B1865, and B1866 all showed reduced or no expression of PTEN protein (Figure 2C,E). PCR analysis using primers specific for WT Pten, un-recombined mutant Pten, and recombined Pten, subsequently showed that for cell lines BB1845, B1865, and B1866, the WT copy of Pten had been lost (Figure S2). However, the WT copy of Pten was retained in BB19, BB57, and BB342 suggesting that the decreased PTEN expression seen in BB57 is due to another mechanism.
Historically, the most common method of in vitro cell culture has been the growth of cells in a two-dimensional monolayer on cell culture plastic. However, three-dimensional (3D) models of tumor growth may be more relevant as they better recapitulate the natural environment and structure of solid tumors. To characterize the growth of these novel cell lines in more physiologic conditions, we used an organoid forming assay (Figure 2D). BB342, BB1845, and B1865 formed large organoids, many with multiple lobes or protrusions from a central mass (Figure 2F,G). In contrast, BB19, BB57, and B1866 formed more numerous organoids, but they remained small in size (Figure 2F,G).

3.3. Evaluation of EMT Properties of BrafV600E Cell Lines

Heterogeneity in cell morphology can be observed among the six cell lines (Figure 1F and Figure 3A). Despite their shared thyroid epithelial cell origin, each display different morphology and adhesion to neighboring cells. Tumor cells of epithelial origin are known to be able to undergo epithelial to mesenchymal transition (EMT) to aid in their escape from the primary tumor site and travel to a distant metastatic site [42]. The process of EMT alters cell morphology and transcriptional output from the tightly adherent epithelial cells expressing epithelial markers such as E-cadherin and cytokeratin to a fusiform appearance with expression of vimentin and fibronectin, typical of mesenchymal cells [42]. As the founding tumors of our novel cell line panel displayed varying degrees of invasion, we sought to determine whether proteins associated with EMT were dysregulated in the cell lines. Western blot analysis showed that BB19, BB57, and B1866 all maintained E-cadherin expression and had little vimentin expression. In contrast, cell lines BB342, BB1845, and B1865 had low to absent expression of E-cadherin and expressed high levels of vimentin (Figure 3B–D). To determine whether loss of E-cadherin correlated with migratory capacity of the cells, a modified Boyden-chamber migration assay was performed. A greater proportion of cells migrated through the transwell in the three cell lines lacking E-cadherin protein expression (BB342, BB1845, and B1865) (Figure 3E,F). In contrast, cells lines that maintained E-cadherin protein expression (BB19, BB57, and B1866) had reduced migratory capacity (Figure 3D,E).

3.4. Growth Suppressive Effect of Pathway-Specific Inhibitors in BrafV600E-Driven Thyroid Tumor Cell Lines in 2D Monolayer and 3D Organoid Cultures

We next evaluated the efficacy of MAPK specific inhibitors on the growth of these novel cell lines. Each cell line was treated with a variety of concentrations of either Dabrafenib (BRAF inhibitor) or Trametinib (MEK1/2 inhibitor) for 3 and 6 days (Figure S3). Both of these inhibitors are being utilized in clinical trials in thyroid cancer [43] and have been utilized in the pediatric setting [44]. In all cell lines, the inhibition of BRAF or MEK1/2 alone resulted in a significant reduction of growth at 3 days (Figure 4A). There was similar efficacy between BRAF and MEK1/2 inhibition in BB19, BB57, BB342, and B1865, whereas BB1845 and B1866 seemed to show slightly more resistance to BRAF inhibition alone at 3 days (Figure 4A). By six days of treatment, BB1845, B1865, and B1866 showed reduced or no significant growth inhibition with the BRAF inhibitor, and BB1845 and B1865 similarly showed reduced or no growth inhibition with the MEK1/2 inhibitor (Figure 4B). Notably, the cell lines showing reduced response to pathway-specific inhibition at 6 days correlated with those derived from higher grade founding tumors and with loss of PTEN expression (Figure S3).
To evaluate the growth-inhibiting effects of these drugs in a more physiologic environment, we established 3D organoid cultures of each cell line and assessed the size of the organoids following treatment with Dabrafenib, Trametinib, or DMSO control for 10 days (Figure 4C). Average organoid size and total organoid area measurements recapitulated expected growth inhibition following inhibitor treatment, as demonstrated in monolayer cultures (Figure 4D,E). Interestingly, in organoid cultures, the duration of response seemed to be maintained. By 6 days post treatment, most of the cell lines showed reduced response to inhibition in 2D, whereas inhibition of spheroid size as maintained up to 10 days in 3D culture.

3.5. Development of a Syngeneic Subcutaneous Tumor Model

The subcutaneous injection of tumor cells into murine hosts has been proven to be a valuable and convenient method for studying tumorigenesis. The very rapid growth and progression of BrafV600E+/−/Pten+/−/TPO-Cre de novo tumors in the mouse thyroid gland precludes many in situ experiments. Therefore, the development of a subcutaneous model in immune competent animals would enable temporal models to study disease progression and to potentially serve as a future model for investigating targeted therapeutic interventions. The strengths of the subcutaneous tumor models coupled with the rapidly penetrant phenotype of our de novo BrafV600E+/−/Pten+/−/TPO-Cre murine model [30,31] led us to develop a syngeneic tumor model of PTC. We evaluated the tumorigenic capacity of the pure 129SvJ genetic background cell lines (BB19, BB57, BB342, and BB1845) in immunocompetent hosts. Wild-type 129SvJ mice were injected with tumor cells subcutaneously into the right hind flank at 4–5 weeks of age. The mice were monitored three times a week for 15 weeks to track the development of tumors.
Similar to our previous studies with HrasG12V-driven tumor cells, variable tumor penetrance was observed between each of the cell lines. BB19 and BB57 were the least penetrant, whereas BB342 had an intermediate penetrance and BB184-injected mice developed tumors with 100% penetrance (Table 2). The latency and progression of tumors was much slower in BB19, BB57, and BB342 where all animals survived to 15 weeks post injection. By contrast, tumors developed very rapidly in mice injected with BB1845 whereby all mice had to be collected before the end point at 15 weeks (Table 2 and Figure 5A).
We sought to determine whether subcutaneous tumors were able to recruit a tumor microenvironment similar to the observed in BrafV600E+/−/Pten+/−/TPO-Cre de novo tumors in the thyroid [31]. To corroborate that characteristic, we analyzed the recruitment of activated fibroblasts that secrete a collagen-rich stroma (Figure S4). Trichrome staining of tumor sections confirmed a collagen-rich microenvironment.
We also sought to determine whether tumors developed in these mice would also undergo spontaneous metastasis to distant sites. All mice were evaluated for distant metastases at sacrifice. Lung metastases were observed in one out of four tumor-bearing mice for the cell line BB342 and in two out of nine tumor-bearing mice for the cell line BB1845 (Table 2 and Figure 5D).

4. Discussion

Cancer cell lines harboring clinically relevant mutations derived from genetically engineered mouse models are invaluable tools for in vitro and in vivo preclinical research to drive advancements in our understanding of disease mechanisms and treatment approaches in thyroid cancer. Here, we sought to generate cell lines that would recapitulate the heterogeneity of PTC and that would allow us to examine factors that contribute to disease progression. PTC is the most prevalent form of thyroid cancer, and incidence is steadily increasing [3]. BRAFV600E is the most common driver mutation in adult PTC and the second most common mutation in pediatric PTC [5,9,10,11,12,45]. Genetic alterations leading to activation in both the MAPK and PI3K/Akt pathway commonly occur in PTC and become increasingly more prevalent as disease progresses to PDTC and ATC [22,23,24,25]. According with these genetic trends, we have previously shown that thyroid-specific expression of BrafV600E alone or combined with Pten inactivation (BrafV600E+/−/Pten−/−/TPO-Cre) leads to the development of PTC in mice [30,31]. In addition to harboring clinically relevant mutations, these models also reproduced similar histopathologic PTC features and stages of progression to PDTC.
In this study, we described six novel independent murine PTC cell lines that were derived from BrafV600E+/−/Pten+/−/TPO-Cre mouse thyroid tumors. These cell lines share the same initiating driver mutations but represent a heterogenous population regarding sex, age, and histologic differentiation status. Additionally, we shown that these cell lines can be used for a variety of in vitro applications including mechanistic studies of targeted inhibitors and formation of 3D organoid cultures, while also demonstrating successful transplantation into syngeneic hosts for in vivo studies. Altogether, our results demonstrate heterogeneity within these cell lines that recapitulate differences observed clinically in the patient population. This novel cell line panel provides an innovative research tool to identify differences among tumor cells with heterogenous responses to treatment and/or invasive behavior despite sharing the same initiating driver mutations. Previous studies have indicated that progression from PTC to PDTC or ATC may be facilitated by mutations leading to the activation of the PI3K/Akt pathway co-occurring with BRAF mutations [22,23,24,25]. Interestingly, cell lines BB1845, B1865, and B1866, which arose from higher grade tumors, have lost their wild-type copy of Pten (Figure S3). These cell lines also demonstrated decreased susceptibility to MAPK pathway inhibition compared with the three cell lines derived from tumors with well-differentiated histology. These findings support the utility of our new cell line panel in studying mechanisms of disease advancement and elucidating prognostic features that may allow the early identification of more indolent versus more invasive tumors despite a common oncogene mutation.
Epithelial-to-mesenchymal transition is associated with loss of E-cadherin expression resulting in loss of cellular adhesion and polarity that leads to an invasive cell phenotype with increased metastatic potential [42,46,47].Previous studies have described EMT in the setting of PTC characterized by over-expression of vimentin in addition to increased tumor invasion and lymph node metastases [48]. We demonstrated that our cell lines showing loss of E-cadherin expression and gain of vimentin expression had greater migratory properties when evaluated using a trans-well migration assay (Figure 5C,D). Interestingly, the thyroid cell lines maintaining the most robust E-cadherin expression formed the smallest yet most frequent organoids. This is in contrast to studies in prostate organoid cultures whereby loss of E-cadherin results in reduced organoid size and disorganized organoid formation [49]. Our in vivo findings further support that loss of E-cadherin expression and increased vimentin expression is advantageous to hematogenous metastasis, indicated by the lung metastasis observed in animals injected with BB342 and BB1845. These novel cell lines will be foundational to future studies to investigate the differential invasiveness and metastatic potential of PTC, and the role of EMT proteins in organoid formation.
The current understanding of EMT in thyroid carcinomas suggests that EMT is implicated in the progression from well-differentiated thyroid cancers to PDTC and ATC [42,50,51,52]. However, our data have mixed results to this end. One of our cell lines with well-differentiated histology, BB342, displayed a mesenchymal-like phenotype characterized by loss of E-cadherin expression, gain of vimentin expression, and increased migratory properties as well as increased metastatic potential in vivo. These findings, in the setting of the well-differentiated histology of the founding tumor, suggest that EMT can occur early in tumorigenesis and may lead to early tumor progression. Interestingly, one of our cell lines that originated from a higher-grade tumor, B1866, was shown to have maintained its epithelial characteristics despite coming from an advanced stage tumor. Together, these results highlight the complex and multifactorial nature of primary tumor progression and the closely related but distinct processes of invasion and metastasis. Furthermore, our findings demonstrate that this novel cell line panel provides a platform that will allow researchers to investigate the mechanisms by which these cell lines developed heterogeneous tumor cell phenotypes despite sharing a common initiating tumorigenic event. Using this model, investigators may gain crucial insights into early prognostic indicators of tumor progression that may supplement the identification of the underlying oncogene mutation and allow for advanced precision in therapeutic strategies for patients with newly diagnosed, recurrent, or RAI-refractory PTC.
In this study, we sought to validate a syngeneic subcutaneous tumor model. This provides a readily accessible method of exploring tumorigenesis that recapitulates the complex interactions between tumor cells, the immune system, and physiologic components of the tumor microenvironment. Previous studies have shown that immune cell recruitment, particularly macrophages, plays a role in thyroid cancer pathogenesis [53,54,55]. Further, immune modulating therapy is now being tested in clinical trials for advance thyroid cancer. These novel models will enable comprehensive investigation of thyroid oncogenesis and response to therapy, including immune modulatory therapies. We demonstrate that the derived cell lines display different latencies and penetrance of tumor formation, and differential capacity to spontaneously metastasize to the lungs. One limitation of our in vivo findings using this model is that the histology observed in the subcutaneous tumors is characteristic of advanced disease for all the cell lines evaluated, despite the heterogeneous differentiation states of the founding tumor for each cell line. However, we anticipate that these cell lines could be used in more advanced in vivo modeling systems, including orthotopic implantation directly into the thyroid of wild-type recipients. It is possible that under those circumstances the resulting tumors could have histological features that are more representative of well-differentiated PTC.
Additionally, this model could be used to evaluate the role of age at onset of disease in PTC development and progression by injecting hosts of different age equivalents to compare penetrance, disease progression, and metastasis [56,57,58,59,60,61]. While the BRAFV600E mutation is very common in adult PTC and moderately common in pediatric PTC, the mechanisms behind the disparity in clinical outcomes between these two populations is not well-understood. We believe that these cell lines model the heterogeneity seen in the population of patients with BRAF-driven PTC and provide a novel platform to investigate features critical to PTC disease progression and response to therapies. These tools provide a valuable platform for more reductionist and mechanistic approaches to investigate Braf-driven disease progression.

5. Conclusions

The establishment of six independent cell lines modeling BRAFV600E-driven papillary thyroid cancer represents an advancement in the development of much-needed validated preclinical models for thyroid cancer research. We anticipate that these cell lines will provide a powerful platform for in vitro and in vivo studies that will lead to the development of precision diagnostic and treatment strategies for patients with newly diagnosed, advanced, or treatment-refractory thyroid cancer. Ongoing work in the laboratory is seeking to identify molecular mechanisms that drive heterogeneity observed in morphology and responses to therapy of these cell lines, despite sharing identical driver mutations.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15030879/s1, Figure S1: Trichrome of original de novo tumors; Figure S2: Braf Recombined DNA gels on cell lines and non-target tissues; Figure S3: Pten Recombined DNA gel on cell lines and non-target tissues and original Pten Genotyping gels for BB1845, B1865, and B1866; Figure S4: Dose-dependent growth-suppressive effects of targeted BRAF and MEK1/2 inhibition in BrafV600E+/−/Pten+/−/TPO-Cre cell lines; Figure S5. Original, uncropped Western blot membrane.

Author Contributions

Conceptualization, G.P.B., V.C.-M., A.B.O. and A.T.F.; Data curation, G.P.B., V.C.-M., Z.S., E.K.W., A.J.B. and A.T.F.; Formal analysis, V.C.-M., A.B.O., J.C.R.-F., N.M. and A.T.F.; Funding acquisition, A.T.F.; Investigation, G.P.B., V.C.-M., A.B.O., M.S. (Madeleine Salwen), M.S. (Michele Scheerer) and E.K.W.; Methodology, G.P.B., V.C.-M., A.B.O., J.C.R.-F., E.K.W. and A.T.F.; Project administration, A.T.F.; Resources, A.T.F.; Supervision, A.T.F.; Validation, A.T.F.; Visualization, A.T.F.; Writing—original draft, G.P.B., V.C.-M., A.B.O., J.C.R.-F. and A.T.F.; Writing—review and editing, G.P.B., V.C.-M., A.B.O., J.C.R.-F., Z.S., M.S. (Michele Scheerer), A.J.B. and A.T.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by NIH R01CA214511 (A.T.F.), NICHD-NIGMS T32 Fellowship in Pediatric Clinical Pharmacology T32GM008562 (A.B.O.) and The Children’s Hospital of Philadelphia Frontier Program Grant (A.T.F. and A.J.B.).

Institutional Review Board Statement

No patient samples were used in the current study therefore no IRB approval was necessary.

Data Availability Statement

Data supporting the findings of this study are contained within the article and Supplementary Materials or are available from the corresponding author (A.T.F.) upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Siegel, R.L.; Jakubowski, C.D.; Fedewa, S.A.; Davis, A.; Azad, N.S. Colorectal Cancer in the Young: Epidemiology, Prevention, Management. Am. Soc. Clin. Oncol. Educ. Book 2020, 40, 1–14. [Google Scholar] [CrossRef] [PubMed]
  2. Qian, Z.J.; Jin, M.C.; Meister, K.D.; Megwalu, U.C. Pediatric Thyroid Cancer Incidence and Mortality Trends in the United States, 1973–2013. JAMA Otolaryngol. Neck Surg. 2019, 145, 617–623. [Google Scholar] [CrossRef] [PubMed]
  3. Lim, H.; Devesa, S.S.; Sosa, J.A.; Check, D.; Kitahara, C.M. Trends in Thyroid Cancer Incidence and Mortality in the United States, 1974–2013. JAMA 2017, 317, 1338–1348. [Google Scholar] [CrossRef]
  4. Saini, S.; Tulla, K.; Maker, A.V.; Burman, K.D.; Prabhakar, B.S. Therapeutic Advances in Anaplastic Thyroid Cancer: A Current Perspective. Mol. Cancer 2018, 17, 154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Ricarte-Filho, J.C.; Ryder, M.; Chitale, D.A.; Rivera, M.; Heguy, A.; Ladanyi, M.; Janakiraman, M.; Solit, D.; Knauf, J.A.; Tuttle, R.M.; 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. 2009, 69, 4885–4893. [Google Scholar] [CrossRef] [Green Version]
  6. Abe, I.; Lam, A.K. Anaplastic Thyroid Carcinoma: Current Issues in Genomics and Therapeutics. Curr. Oncol. Rep. 2021, 23, 31. [Google Scholar] [CrossRef]
  7. Hescheler, D.A.; Hartmann, M.J.M.; Riemann, B.; Michel, M.; Bruns, C.J.; Alakus, H.; Chiapponi, C. Anaplastic thyroid cancer: Genome-based search for new targeted therapy options. Endocr. Connect. 2022, 11, e210624. [Google Scholar] [CrossRef]
  8. Franco, A.T.; Ricarte-Filho, J.C.; Isaza, A.; Jones, Z.; Jain, N.; Mostoufi-Moab, S.; Surrey, L.; Laetsch, T.W.; Li, M.M.; DeHart, J.C.; et al. Fusion Oncogenes Are Associated with Increased Metastatic Capacity and Persistent Disease in Pediatric Thyroid Cancers. J. Clin. Oncol. 2022, 40, 1081–1090. [Google Scholar] [CrossRef]
  9. Puxeddu, E.; Moretti, S.; Elisei, R.; Romei, C.; Pascucci, R.; Martinelli, M.; Marino, C.; Avenia, N.; Rossi, E.D.; Fadda, G.; et al. BRAF(V599E) mutation is the leading genetic event in adult sporadic papillary thyroid carcinomas. J. Clin. Endocrinol. Metab. 2004, 89, 2414–2420. [Google Scholar] [CrossRef] [Green Version]
  10. Bauer, A.J. Molecular Genetics of Thyroid Cancer in Children and Adolescents. Endocrinol. Metab. Clin. N. Am. 2017, 46, 389–403. [Google Scholar] [CrossRef]
  11. Xing, M. Braf Mutation in Thyroid Cancer. Endocr. Relat. Cancer 2005, 12, 245–262. [Google Scholar] [CrossRef] [Green Version]
  12. Xing, M. Molecular Pathogenesis and Mechanisms of Thyroid Cancer. Nat. Rev. Cancer 2013, 13, 184–199. [Google Scholar] [CrossRef] [PubMed]
  13. Lee, J.H.; Lee, E.S.; Kim, Y.S. Clinicopathologic Significance of Braf V600e Mutation in Papillary Carcinomas of the Thyroid: A Meta-Analysis. Cancer 2007, 110, 38–46. [Google Scholar] [CrossRef] [PubMed]
  14. Elisei, R.; Ugolini, C.; Viola, D.; Lupi, C.; Biagini, A.; Giannini, R.; Romei, C.; Miccoli, P.; Pinchera, A.; Basolo, F. BRAF(V600E) mutation and outcome of patients with papillary thyroid carcinoma: A 15-year median follow-up study. J. Clin. Endocrinol. Metab. 2008, 93, 3943–3949. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Frasca, F.; Nucera, C.; Pellegriti, G.; Gangemi, P.; Attard, M.; Stella, M.; Loda, M.; Vella, V.; Giordano, C.; Trimarchi, F.; et al. BRAF(V600E) mutation and the biology of papillary thyroid cancer. Endocr. Relat. Cancer 2008, 15, 191–205. [Google Scholar] [CrossRef]
  16. Namba, H.; Nakashima, M.; Hayashi, T.; Hayashida, N.; Maeda, S.; Rogounovitch, T.I.; Ohtsuru, A.; Saenko, V.A.; Kanematsu, T.; Yamashita, S. Clinical implication of hot spot BRAF mutation, V599E, in papillary thyroid cancers. J. Clin. Endocrinol. Metab. 2003, 88, 4393–4397. [Google Scholar] [CrossRef] [Green Version]
  17. Kebebew, E.; Weng, J.; Bauer, J.; Ranvier, G.; Clark, O.H.; Duh, Q.Y.; Shibru, D.; Bastian, B.; Griffin, A. The Prevalence and Prognostic Value of Braf Mutation in Thyroid Cancer. Ann. Surg. 2007, 246, 466–471; discussion 70–71. [Google Scholar] [CrossRef]
  18. Liu, J.; Liu, R.; Shen, X.; Zhu, G.; Li, B.; Xing, M. The Genetic Duet of Braf V600e and Tert Promoter Mutations Robustly Predicts Loss of Radioiodine Avidity in Recurrent Papillary Thyroid Cancer. J. Nucl. Med. 2020, 61, 177–182. [Google Scholar] [CrossRef] [Green Version]
  19. Pekova, B.; Sykorova, V.; Dvorakova, S.; Vaclavikova, E.; Moravcova, J.; Katra, R.; Astl, J.; Vlcek, P.; Kodetova, D.; Vcelak, J.; et al. Ret, Ntrk, Alk, Braf, and Met Fusions in a Large Cohort of Pediatric Papillary Thyroid Carcinomas. Thyroid 2020, 30, 1771–1780. [Google Scholar] [CrossRef]
  20. Henke, L.E.; Perkins, S.M.; Pfeifer, J.D.; Ma, C.; Chen, Y.; DeWees, T.; Grigsby, P.W. Braf V600e Mutational Status in Pediatric Thyroid Cancer. Pediatr. Blood Cancer 2014, 61, 1168–1172. [Google Scholar] [CrossRef]
  21. Givens, D.J.; Buchmann, L.O.; Agarwal, A.M.; Grimmer, J.F.; Hunt, J.P. Braf V600e Does Not Predict Aggressive Features of Pediatric Papillary Thyroid Carcinoma. Laryngoscope 2014, 124, E389–E393. [Google Scholar] [CrossRef] [PubMed]
  22. Hou, P.; Liu, D.; Shan, Y.; Hu, S.; Studeman, K.; Condouris, S.; Wang, Y.; Trink, A.; El-Naggar, A.K.; Tallini, G.; et al. Genetic Alterations and Their Relationship in the Phosphatidylinositol 3-Kinase/Akt Pathway in Thyroid Cancer. Clin. Cancer Res. 2007, 13, 1161–1170. [Google Scholar] [CrossRef] [Green Version]
  23. Santarpia, L.; El-Naggar, A.K.; Cote, G.J.; Myers, J.N.; Sherman, S.I. Phosphatidylinositol 3-Kinase/Akt and Ras/Raf-Mitogen-Activated Protein Kinase Pathway Mutations in Anaplastic Thyroid Cancer. J. Clin. Endocrinol. Metab. 2008, 93, 278–284. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Song, E.; Jin, M.; Jang, A.; Jeon, M.J.; Song, D.E.; Yoo, H.J.; Kim, W.B.; Shong, Y.K.; Kim, W.G. Mutation in Genes Encoding Key Functional Groups Additively Increase Mortality in Patients with Braf(V600e)-Mutant Advanced Papillary Thyroid Carcinoma. Cancers 2021, 13, 5846. [Google Scholar] [CrossRef] [PubMed]
  25. Liu, Z.; Hou, P.; Ji, M.; Guan, H.; Studeman, K.; Jensen, K.; Vasko, V.; El-Naggar, A.K.; Xing, M. Highly Prevalent Genetic Alterations in Receptor Tyrosine Kinases and Phosphatidylinositol 3-Kinase/Akt and Mitogen-Activated Protein Kinase Pathways in Anaplastic and Follicular Thyroid Cancers. J. Clin. Endocrinol. Metab. 2008, 93, 3106–3116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Schweppe, R.E.; Klopper, J.P.; Korch, C.; Pugazhenthi, U.; Benezra, M.; Knauf, J.A.; Fagin, J.A.; Marlow, L.A.; Copland, J.A.; Smallridge, R.C.; et al. 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. 2008, 93, 4331–4341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Knauf, J.A.; Luckett, K.A.; Chen, K.Y.; Voza, F.; Socci, N.D.; Ghossein, R.; Fagin, J.A. Hgf/Met Activation Mediates Resistance to Braf Inhibition in Murine Anaplastic Thyroid Cancers. J. Clin. Investig. 2018, 128, 4086–4097. [Google Scholar] [CrossRef]
  28. Vanden Borre, P.; McFadden, D.G.; Gunda, V.; Sadow, P.M.; Varmeh, S.; Bernasconi, M.; Jacks, T.; Parangi, S. The Next Generation of Orthotopic Thyroid Cancer Models: Immunocompetent Orthotopic Mouse Models of Braf V600e-Positive Papillary and Anaplastic Thyroid Carcinoma. Thyroid 2014, 24, 705–714. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Caperton, C.O.; Jolly, L.A.; Massoll, N.; Bauer, A.J.; Franco, A.T. Development of Novel Follicular Thyroid Cancer Models Which Progress to Poorly Differentiated and Anaplastic Thyroid Cancer. Cancers 2021, 13, 1094. [Google Scholar] [CrossRef]
  30. Franco, A.T.; Malaguarnera, R.; Refetoff, S.; Liao, X.H.; Lundsmith, E.; Kimura, S.; Pritchard, C.; Marais, R.; Davies, T.F.; Weinstein, L.S.; et al. Thyrotrophin Receptor Signaling Dependence of Braf-Induced Thyroid Tumor Initiation in Mice. Proc. Natl. Acad. Sci. USA 2011, 108, 1615–1620. [Google Scholar] [CrossRef]
  31. Jolly, L.A.; Novitskiy, S.; Owens, P.; Massoll, N.; Cheng, N.; Fang, W.; Moses, H.L.; Franco, A.T. Fibroblast-Mediated Collagen Remodeling within the Tumor Microenvironment Facilitates Progression of Thyroid Cancers Driven by Brafv600e and Pten Loss. Cancer Res. 2016, 76, 1804–1813. [Google Scholar] [CrossRef] [Green Version]
  32. Vilgelm, A.E.; Bergdorf, K.; Wolf, M.; Bharti, V.; Shattuck-Brandt, R.; Blevins, A.; Jones, C.; Phifer, C.; Lee, M.; Lowe, C.; et al. Fine-Needle Aspiration-Based Patient-Derived Cancer Organoids. iScience 2020, 23, 101408. [Google Scholar] [CrossRef]
  33. Fernandez, L.P.; Lopez-Marquez, A.; Santisteban, P. Thyroid Transcription Factors in Development, Differentiation and Disease. Nat. Rev. Endocrinol. 2015, 11, 29–42. [Google Scholar] [CrossRef] [Green Version]
  34. Kimura, S. Thyroid-Specific Transcription Factors and Their Roles in Thyroid Cancer. J. Thyroid Res. 2011, 2011, 710213. [Google Scholar] [CrossRef] [Green Version]
  35. Kanjanahattakij, N.; Chayangsu, P.; Kanoksil, W.; Chontong, S.; Sriphrapradang, C. Pitfall in Immunohistochemical Staining for Thyroglobulin in Case of Thyroid Metastasis from Lung Carcinoma. Cytojournal 2015, 12, 27. [Google Scholar] [CrossRef]
  36. Liu, T.; Men, Q.; Su, X.; Chen, W.; Zou, L.; Li, Q.; Song, M.; Ouyang, D.; Chen, Y.; Li, Z.; et al. Downregulated Expression of Tshr Is Associated with Distant Metastasis in Thyroid Cancer. Oncol. Lett. 2017, 14, 7506–7512. [Google Scholar] [CrossRef]
  37. Riesco-Eizaguirre, G.; Rodriguez, I.; De la Vieja, A.; Costamagna, E.; Carrasco, N.; Nistal, M.; Santisteban, P. The Brafv600e Oncogene Induces Transforming Growth Factor Beta Secretion Leading to Sodium Iodide Symporter Repression and Increased Malignancy in Thyroid Cancer. Cancer Res. 2009, 69, 8317–8325. [Google Scholar] [CrossRef] [Green Version]
  38. Tavares, C.; Coelho, M.J.; Eloy, C.; Melo, M.; da Rocha, A.G.; Pestana, A.; Batista, R.; Ferreira, L.B.; Rios, E.; Selmi-Ruby, S.; et al. Nis Expression in Thyroid Tumors, Relation with Prognosis Clinicopathological and Molecular Features. Endocr. Connect. 2018, 7, 78–90. [Google Scholar] [CrossRef] [Green Version]
  39. Spitzweg, C.; Bible, K.C.; Hofbauer, L.C.; Morris, J.C. Advanced Radioiodine-Refractory Differentiated Thyroid Cancer: The Sodium Iodide Symporter and Other Emerging Therapeutic Targets. Lancet Diabetes Endocrinol. 2014, 2, 830–842. [Google Scholar] [CrossRef]
  40. Hou, P.; Bojdani, E.; Xing, M. Induction of Thyroid Gene Expression and Radioiodine Uptake in Thyroid Cancer Cells by Targeting Major Signaling Pathways. J. Clin. Endocrinol. Metab. 2010, 95, 820–828. [Google Scholar] [CrossRef]
  41. Fu, H.; Cheng, L.; Sa, R.; Jin, Y.; Chen, L. Combined Tazemetostat and Mapki Enhances Differentiation of Papillary Thyroid Cancer Cells Harbouring Braf(V600e) by Synergistically Decreasing Global Trimethylation of H3k27. J. Cell. Mol. Med. 2020, 24, 3336–3345. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Shakib, H.; Rajabi, S.; Dehghan, M.H.; Mashayekhi, F.J.; Safari-Alighiarloo, N.; Hedayati, M. Epithelial-to-Mesenchymal Transition in Thyroid Cancer: A Comprehensive Review. Endocrine 2019, 66, 435–455. [Google Scholar] [CrossRef]
  43. Busaidy, N.L.; Konda, B.; Wei, L.; Wirth, L.J.; Devine, C.; Daniels, G.A.; DeSouza, J.A.; Poi, M.; Seligson, N.D.; Cabanillas, M.E.; et al. Dabrafenib Versus Dabrafenib + Trametinib in ≪I≫Braf≪/I≫-Mutated Radioactive Iodine Refractory Differentiated Th205yroid Cancer: Results of a Randomized, Phase 2, Open-Label Multicenter Trial. Thyroid 2022, 32, 1184–1992. [Google Scholar] [CrossRef]
  44. Kieran, M.W.; Geoerger, B.; Dunkel, I.J.; Broniscer, A.; Hargrave, D.; Hingorani, P.; Aerts, I.; Bertozzi, A.I.; Cohen, K.J.; Hummel, T.R.; et al. A Phase I and Pharmacokinetic Study of Oral Dabrafenib in Children and Adolescent Patients with Recurrent or Refractory Braf V600 Mutation-Positive Solid Tumors. Clin. Cancer Res. 2019, 25, 7294–7302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Paulson, V.A.; Rudzinski, E.R.; Hawkins, D.S. Thyroid Cancer in the Pediatric Population. Genes 2019, 10, 723. [Google Scholar] [CrossRef] [Green Version]
  46. Yang, J.; Weinberg, R.A. Epithelial-Mesenchymal Transition: At the Crossroads of Development and Tumor Metastasis. Dev. Cell 2008, 14, 818–829. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Thiery, J.P. Epithelial-Mesenchymal Transitions in Tumour Progression. Nat. Rev. Cancer 2002, 2, 442–454. [Google Scholar] [CrossRef]
  48. Vasko, V.; Espinosa, A.V.; Scouten, W.; He, H.; Auer, H.; Liyanarachchi, S.; Larin, A.; Savchenko, V.; Francis, G.L.; de la Chapelle, A.; et al. Gene Expression and Functional Evidence of Epithelial-to-Mesenchymal Transition in Papillary Thyroid Carcinoma Invasion. Proc. Natl. Acad. Sci. USA 2007, 104, 2803–2808. [Google Scholar] [CrossRef] [Green Version]
  49. Olson, A.; Le, V.; Aldahl, J.; Yu, E.J.; Hooker, E.; He, Y.; Lee, D.H.; Kim, W.K.; Cardiff, R.D.; Geradts, J.; et al. The Comprehensive Role of E-Cadherin in Maintaining Prostatic Epithelial Integrity During Oncogenic Transformation and Tumor Progression. PLoS Genet. 2019, 15, e1008451. [Google Scholar] [CrossRef] [Green Version]
  50. Wiseman, S.M.; Griffith, O.L.; Deen, S.; Rajput, A.; Masoudi, H.; Gilks, B.; Goldstein, L.; Gown, A.; Jones, S.J. Identification of Molecular Markers Altered during Transformation of Differentiated into Anaplastic Thyroid Carcinoma. Arch. Surg. 2007, 142, 717–729; discussion 27–29. [Google Scholar] [CrossRef]
  51. Knauf, J.A.; Sartor, M.A.; Medvedovic, M.; Lundsmith, E.; Ryder, M.; Salzano, M.; Nikiforov, Y.E.; Giordano, T.J.; Ghossein, R.A.; Fagin, J.A. Progression of Braf-Induced Thyroid Cancer Is Associated with Epithelial-Mesenchymal Transition Requiring Concomitant Map Kinase and Tgfbeta Signaling. Oncogene 2011, 30, 3153–3162. [Google Scholar] [CrossRef] [Green Version]
  52. Jung, C.W.; Han, K.H.; Seol, H.; Park, S.; Koh, J.S.; Lee, S.S.; Kim, M.J.; Choi, I.J.; Myung, J.K. Expression of Cancer Stem Cell Markers and Epithelial-Mesenchymal Transition-Related Factors in Anaplastic Thyroid Carcinoma. Int. J. Clin. Exp. Pathol. 2015, 8, 560–568. [Google Scholar]
  53. Gunda, V.; Gigliotti, B.; Ndishabandi, D.; Ashry, T.; McCarthy, M.; Zhou, Z.; Amin, S.; Freeman, G.J.; Alessandrini, A.; Parangi, S. Combinations of Braf Inhibitor and Anti-Pd-1/Pd-L1 Antibody Improve Survival and Tumour Immunity in an Immunocompetent Model of Orthotopic Murine Anaplastic Thyroid Cancer. Br. J. Cancer 2018, 119, 1223–1232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Ryder, M.; Ghossein, R.A.; Ricarte-Filho, J.C.; Knauf, J.A.; Fagin, J.A. Increased Density of Tumor-Associated Macrophages Is Associated with Decreased Survival in Advanced Thyroid Cancer. Endocr. Relat. Cancer 2008, 15, 1069–1074. [Google Scholar] [CrossRef] [Green Version]
  55. Zhang, P.; Guan, H.; Yuan, S.; Cheng, H.; Zheng, J.; Zhang, Z.; Liu, Y.; Yu, Y.; Meng, Z.; Zheng, X.; et al. Targeting Myeloid Derived Suppressor Cells Reverts Immune Suppression and Sensitizes Braf-Mutant Papillary Thyroid Cancer to Mapk Inhibitors. Nat. Commun. 2022, 13, 1588. [Google Scholar] [CrossRef] [PubMed]
  56. Kreisle, R.A.; Stebler, B.A.; Ershler, W.B. Effect of Host Age on Tumor-Associated Angiogenesis in Mice. J. Natl. Cancer Inst. 1990, 82, 44–47. [Google Scholar] [CrossRef]
  57. Fane, M.; Weeraratna, A.T. How the Ageing Microenvironment Influences Tumour Progression. Nat. Rev. Cancer 2020, 20, 89–106. [Google Scholar] [CrossRef]
  58. Jackaman, C.; Radley-Crabb, H.G.; Soffe, Z.; Shavlakadze, T.; Grounds, M.D.; Nelson, D.J. Targeting Macrophages Rescues Age-Related Immune Deficiencies in C57bl/6j Geriatric Mice. Aging Cell 2013, 12, 345–357. [Google Scholar] [CrossRef]
  59. Reed, M.J.; Karres, N.; Eyman, D.; Cruz, A.; Brekken, R.A.; Plymate, S. The Effects of Aging on Tumor Growth and Angiogenesis Are Tumor-Cell Dependent. Int. J. Cancer 2007, 120, 753–760. [Google Scholar] [CrossRef]
  60. Vincent-Chong, V.K.; DeJong, H.; Rich, L.J.; Patti, A.; Merzianu, M.; Hershberger, P.A.; Seshadri, M. Impact of Age on Disease Progression and Microenvironment in Oral Cancer. J. Dent. Res. 2018, 97, 1268–1276. [Google Scholar] [CrossRef]
  61. Pettan-Brewer, C.; Morton, J.; Coil, R.; Hopkins, H.; Fatemie, S.; Ladiges, W. B16 Melanoma Tumor Growth Is Delayed in Mice in an Age-Dependent Manner. Pathobiol. Aging Age Relat. Dis. 2012, 2, 19182. [Google Scholar] [CrossRef] [PubMed]
Cancers 15 00879 g001 550
Figure 1. Histology and pathological features of BrafV600E driven murine thyroid tumors and the morphology of stable tumor cell lines derived from them. (A) Hematoxylin and eosin (H&E) staining of de novo tumors for cell lines BB19, BB57, BB342, BB1845, B1865, and B1866. Scale bar in black is 750 µM. (B) High magnification of BB57 H&E showing papillary thyroid carcinoma with fine chromatin, nuclear grooves (blue arrow) and psammoma bodies (black arrow). (C) High magnification of BB57 H&E showing papillary thyroid carcinoma (blue arrow) with invasion into the muscle (black arrow). (D) High magnification of an extra-thyroidal nodule outside of the thyroid capsule (blue arrow), adjacent to intra-thyroidal tumor (orange arrow) in B1866. (E) High magnification of B1866 H&E showing moderately differentiated papillary thyroid carcinoma, which maintains papillae and grooves but has associated necrosis (blue arrow). (F) Representative phase contrast images of the individual cell lines in culture derived from the indicated tumors. Scale bar in black is 300 µM.
Figure 1. Histology and pathological features of BrafV600E driven murine thyroid tumors and the morphology of stable tumor cell lines derived from them. (A) Hematoxylin and eosin (H&E) staining of de novo tumors for cell lines BB19, BB57, BB342, BB1845, B1865, and B1866. Scale bar in black is 750 µM. (B) High magnification of BB57 H&E showing papillary thyroid carcinoma with fine chromatin, nuclear grooves (blue arrow) and psammoma bodies (black arrow). (C) High magnification of BB57 H&E showing papillary thyroid carcinoma (blue arrow) with invasion into the muscle (black arrow). (D) High magnification of an extra-thyroidal nodule outside of the thyroid capsule (blue arrow), adjacent to intra-thyroidal tumor (orange arrow) in B1866. (E) High magnification of B1866 H&E showing moderately differentiated papillary thyroid carcinoma, which maintains papillae and grooves but has associated necrosis (blue arrow). (F) Representative phase contrast images of the individual cell lines in culture derived from the indicated tumors. Scale bar in black is 300 µM.
Cancers 15 00879 g001
Cancers 15 00879 g002 550
Figure 2. Expression of thyroid-specific genes, proliferation rate, and activation of MAPK and PI3K pathway signaling in novel BRAFV600E murine cell lines. (A) Real-time qPCR analysis of the expression levels of thyroid-specific genes in primary wild-type (WT) thyrocytes and thyroid tumor cell lines in culture. Each sample was analyzed in triplicate, graph represents n = 3. Fold change is indicated on the y-axis. (B) Growth curves quantitating the growth of the cell lines in 2D cultures using Cell-Titer Glo cell viability assay. Day 0 represents 18 hours after plating. Each cell line was plated in quadruplicate, graph represents n = 3, * p > 0.05 compared to BB19, ** p > 0.002 compared to BB19. (C) Western Blot analysis characterizing differential activation of the MAPK and PI3K pathway between cell lines and differential levels of Pten. (D) Quantitation of ratio of pERK to total ERK levels in each cell line following 10 min 10% FBS simulation after overnight growth factor withdrawal. n = minimum of 3 independent replicates. (E) Quantitation of basal pAKT levels in the 6 cell lines. n = minimum of 3 independent replicates. (F) Representative images of organoid cultures from each cell line. Top panel represents bright field images imaged at 4× magnification, bottom panel is immunofluorescence images with Phalloidin in pink and DAPI (Nuclei) in blue, imaged at 40× magnification. (G) Quantitation of average organoid size from each cell line. Average sizes from two independent replicates.
Figure 2. Expression of thyroid-specific genes, proliferation rate, and activation of MAPK and PI3K pathway signaling in novel BRAFV600E murine cell lines. (A) Real-time qPCR analysis of the expression levels of thyroid-specific genes in primary wild-type (WT) thyrocytes and thyroid tumor cell lines in culture. Each sample was analyzed in triplicate, graph represents n = 3. Fold change is indicated on the y-axis. (B) Growth curves quantitating the growth of the cell lines in 2D cultures using Cell-Titer Glo cell viability assay. Day 0 represents 18 hours after plating. Each cell line was plated in quadruplicate, graph represents n = 3, * p > 0.05 compared to BB19, ** p > 0.002 compared to BB19. (C) Western Blot analysis characterizing differential activation of the MAPK and PI3K pathway between cell lines and differential levels of Pten. (D) Quantitation of ratio of pERK to total ERK levels in each cell line following 10 min 10% FBS simulation after overnight growth factor withdrawal. n = minimum of 3 independent replicates. (E) Quantitation of basal pAKT levels in the 6 cell lines. n = minimum of 3 independent replicates. (F) Representative images of organoid cultures from each cell line. Top panel represents bright field images imaged at 4× magnification, bottom panel is immunofluorescence images with Phalloidin in pink and DAPI (Nuclei) in blue, imaged at 40× magnification. (G) Quantitation of average organoid size from each cell line. Average sizes from two independent replicates.
Cancers 15 00879 g002
Cancers 15 00879 g003 550
Figure 3. Characterization of epithelial and mesenchymal features of BrafV600E cell lines. (A) Representative differential interference contrast (DIC) images of each cell line. Top: DIC imaging of live cells at 40×. Bottom: DIC imaging of live cells at 63×. (B) Western blot analysis of levels of epithelial (E cadherin) and mesenchymal (Vimentin) proteins. (C) Quantitation of E cadherin expression in each of the six cell lines N = 3, ** p < 0.002 compared to levels in BB19. (D) Quantitation of vimentin expression in each of the six cell lines N = 3. No Statistically different comparison compared to levels in BB19. (E) Representative images of cells that have migrated through the membrane of a trans-well, imaged at 4× magnification. (F) Quantification of the cells that have migrated through a trans-well. Percent of the membrane area covered by cells that have migrated through the filter is indicated on the y-axis. N = 2 independent experiments. ** p < 0.02 compared to BB19, *** p < 0.0001 compared to BB19.
Figure 3. Characterization of epithelial and mesenchymal features of BrafV600E cell lines. (A) Representative differential interference contrast (DIC) images of each cell line. Top: DIC imaging of live cells at 40×. Bottom: DIC imaging of live cells at 63×. (B) Western blot analysis of levels of epithelial (E cadherin) and mesenchymal (Vimentin) proteins. (C) Quantitation of E cadherin expression in each of the six cell lines N = 3, ** p < 0.002 compared to levels in BB19. (D) Quantitation of vimentin expression in each of the six cell lines N = 3. No Statistically different comparison compared to levels in BB19. (E) Representative images of cells that have migrated through the membrane of a trans-well, imaged at 4× magnification. (F) Quantification of the cells that have migrated through a trans-well. Percent of the membrane area covered by cells that have migrated through the filter is indicated on the y-axis. N = 2 independent experiments. ** p < 0.02 compared to BB19, *** p < 0.0001 compared to BB19.
Cancers 15 00879 g003
Cancers 15 00879 g004 550
Figure 4. Growth-suppressive effects of targeted BRAF and MEK1/2 inhibition in BrafV600E+/−/Pten+/−/TPO-Cre cell lines in 2D monolayer and 3D organoid cultures. (A) Summary plot showing viability of novel cell lines in 2D culture following treatment with the BRAF inhibitor Dabrafenib (D), MEK1/2 inhibitor Trametinib (T) or DMSO vehicle (V) compared to control cells grown in standard culture medium. Cells were treated with selected concentrations of each inhibitor or vehicle for 3 days prior to analysis using CellTiter Glo Luminescent Cell Viability Assay. N = 3 independent experiments in triplicate. (B) Summary plot showing viability of tumor cell lines treated with inhibitors as in (A) for 6 days prior to analysis using CellTiter Glo Luminescent Cell Viability Assay. (C) Representative images of novel cell lines grown in 3D organoid culture and treated with the BRAF inhibitor Dabrafenib, MEK1/2 inhibitor Trametinib or DMSO vehicle for 10 days. Images were captured at 4× magnification. (D) Quantitation of average organoid area of novel cell lines treated as in (C).
Figure 4. Growth-suppressive effects of targeted BRAF and MEK1/2 inhibition in BrafV600E+/−/Pten+/−/TPO-Cre cell lines in 2D monolayer and 3D organoid cultures. (A) Summary plot showing viability of novel cell lines in 2D culture following treatment with the BRAF inhibitor Dabrafenib (D), MEK1/2 inhibitor Trametinib (T) or DMSO vehicle (V) compared to control cells grown in standard culture medium. Cells were treated with selected concentrations of each inhibitor or vehicle for 3 days prior to analysis using CellTiter Glo Luminescent Cell Viability Assay. N = 3 independent experiments in triplicate. (B) Summary plot showing viability of tumor cell lines treated with inhibitors as in (A) for 6 days prior to analysis using CellTiter Glo Luminescent Cell Viability Assay. (C) Representative images of novel cell lines grown in 3D organoid culture and treated with the BRAF inhibitor Dabrafenib, MEK1/2 inhibitor Trametinib or DMSO vehicle for 10 days. Images were captured at 4× magnification. (D) Quantitation of average organoid area of novel cell lines treated as in (C).
Cancers 15 00879 g004
Cancers 15 00879 g005 550
Figure 5. Cell lines develop syngeneic subcutaneous tumors. (A) Kaplan Meyer plot of tumor free survival of animals subcutaneously injected with Braf-driven tumor cell lines. (Note only 4 cell lines were on congenic 129SvJ background and therefore able to be syngeneically injected into WT 129SvJ recipients). (B) Representative H&E and images of subcutaneous tumors. Black scale bars are 750 µM. (C) Representative Trichrome images of subcutaneous tumors. Collagen is blue, nuclei are black, and muscle tissue erythrocytes and cytoplasm are red to pink. Black scale bar is 200 µM. (D) Representative H&E image of a lung metastasis. Black scale bar is 300 µM.
Figure 5. Cell lines develop syngeneic subcutaneous tumors. (A) Kaplan Meyer plot of tumor free survival of animals subcutaneously injected with Braf-driven tumor cell lines. (Note only 4 cell lines were on congenic 129SvJ background and therefore able to be syngeneically injected into WT 129SvJ recipients). (B) Representative H&E and images of subcutaneous tumors. Black scale bars are 750 µM. (C) Representative Trichrome images of subcutaneous tumors. Collagen is blue, nuclei are black, and muscle tissue erythrocytes and cytoplasm are red to pink. Black scale bar is 200 µM. (D) Representative H&E image of a lung metastasis. Black scale bar is 300 µM.
Cancers 15 00879 g005
Table
Table 1. Demographic information of the mice from which tumor cell lines were derived.
Table 1. Demographic information of the mice from which tumor cell lines were derived.
Cell LineGenotypeSexAge (Weeks) at Collection
BB19BrafV600EHet/PtenHet/TPO-Cre+Female6
BB57BrafV600EHet/PtenHet/TPO-Cre+Male4
BB342BrafV600EHet/PtenHet/TPO-Cre+Male3
BB1845BrafV600EHet/PtenHet/TPO-Cre+Female32
B1865BrafV600EHet/PtenHet/TPO-Cre+Male24
B1866BrafV600EHet/PtenHet/TPO-Cre+Male24
Table
Table 2. Summary of tumor penetrance and lung metastasis of injected cell lines.
Table 2. Summary of tumor penetrance and lung metastasis of injected cell lines.
Cell LineNumber of Animals with Tumors% Penetrance
(/N)		Number of Animals with Lung Metastasis
BB19110%
(1/10)		0
BB57220%
(2/10)		0
BB3424 40%
(4/10)		1
BB18459100%
(9/9)		2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Branigan, G.P.; Casado-Medrano, V.; O’Neill, A.B.; Ricarte-Filho, J.C.; Massoll, N.; Salwen, M.; Spangler, Z.; Scheerer, M.; Williamson, E.K.; Bauer, A.J.; et al. Development of Novel Murine BRAFV600E-Driven Papillary Thyroid Cancer Cell Lines for Modeling of Disease Progression and Preclinical Evaluation of Therapeutics. Cancers 2023, 15, 879. https://doi.org/10.3390/cancers15030879

AMA Style

Branigan GP, Casado-Medrano V, O’Neill AB, Ricarte-Filho JC, Massoll N, Salwen M, Spangler Z, Scheerer M, Williamson EK, Bauer AJ, et al. Development of Novel Murine BRAFV600E-Driven Papillary Thyroid Cancer Cell Lines for Modeling of Disease Progression and Preclinical Evaluation of Therapeutics. Cancers. 2023; 15(3):879. https://doi.org/10.3390/cancers15030879

Chicago/Turabian Style

Branigan, Grace Purvis, Victoria Casado-Medrano, Alison B. O’Neill, Julio C. Ricarte-Filho, Nicole Massoll, Madeleine Salwen, Zachary Spangler, Michele Scheerer, Edward K. Williamson, Andrew J. Bauer, and et al. 2023. "Development of Novel Murine BRAFV600E-Driven Papillary Thyroid Cancer Cell Lines for Modeling of Disease Progression and Preclinical Evaluation of Therapeutics" Cancers 15, no. 3: 879. https://doi.org/10.3390/cancers15030879

APA Style

Branigan, G. P., Casado-Medrano, V., O’Neill, A. B., Ricarte-Filho, J. C., Massoll, N., Salwen, M., Spangler, Z., Scheerer, M., Williamson, E. K., Bauer, A. J., & Franco, A. T. (2023). Development of Novel Murine BRAFV600E-Driven Papillary Thyroid Cancer Cell Lines for Modeling of Disease Progression and Preclinical Evaluation of Therapeutics. Cancers, 15(3), 879. https://doi.org/10.3390/cancers15030879

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop