Skip to main content
NIHPA Author Manuscripts logoLink to NIHPA Author Manuscripts
. Author manuscript; available in PMC: 2016 Nov 5.
Published in final edited form as: Cell Stem Cell. 2015 Oct 22;17(5):527–542. doi: 10.1016/j.stem.2015.09.004

Regeneration of thyroid function by transplantation of differentiated pluripotent stem cells

Anita A Kurmann 1,2,*, Maria Serra 2,3,*, Finn Hawkins 2,3, Scott A Rankin 4, Munemasa Mori 2,3, Inna Astapova 1, Soumya Ullas 6, Sui Lin 5, Melanie Bilodeau 7, Janet Rossant 7,8, Jyh C Jean 2,3, Laertis Ikonomou 2,3, Robin R Deterding 9, John M Shannon 5, Aaron M Zorn 4, Anthony N Hollenberg 1,#, Darrell N Kotton 2,3,#
PMCID: PMC4666682  NIHMSID: NIHMS732811  PMID: 26593959

Abstract

Differentiation of functional thyroid epithelia from pluripotent stem cells (PSCs) holds the potential for application in regenerative medicine. However, progress toward this goal is hampered by incomplete understanding of the signaling pathways needed for directed differentiation without forced over-expression of exogenous transgenes. Here we use mouse PSCs to identify key conserved roles for BMP and FGF signaling in regulating thyroid lineage specification from foregut endoderm in mouse and Xenopus. Thyroid progenitors derived from mouse PSCs can be matured into thyroid follicular organoids that provide functional secretion of thyroid hormones in vivo and rescue hypothyroid mice after transplantation. Moreover, by stimulating the same pathways we were also able to derive human thyroid progenitors from normal and disease-specific iPSCs generated from patients with hypothyroidism resulting from NKX2-1 haploinsufficiency. Our studies have therefore uncovered the regulatory mechanisms that underlie early thyroid organogenesis and provide a significant step toward cell-based regenerative therapy for hypothyroidism.

Introduction

Recent progress in the differentiation of pluripotent stem cells (PSCs) in vitro has allowed the derivation of desired cell lineages and in some instances their in vitro self-assembly into 3D structures, referred to as organoids (Lancaster and Knoblich, 2014; McCracken et al., 2014). Transplantation and in vivo function of engineered cells has been less successful, typically due to poor engraftment. The derivation of PSC-derived endocrine tissues, such as pancreatic islets or thyroid follicles, provides particularly attractive opportunities to achieve in vivo function without orthotopic transplantation, since engraftment of these hormone-secreting tissues in any location with access to circulating blood would potentially achieve function and even clinical rescue of hormone deficiency. Indeed, the differentiation of PSCs into pancreatic islet-like cells has produced cells able to secrete insulin in vivo following transplantation (Pagliuca et al., 2014). In marked contrast, thyroid epithelial cells displaying in vivo functional potential have only been generated from PSCs through forced over-expression of transcription factors (Antonica et al., 2012; Ma et al., 2013, 2015). The directed differentiation of PSCs in vitro into thyroid epithelial cells using growth factor supplemented medias previously has only resulted in immature cells that fail to express the full genetic program necessary for functional thyroid hormone biosynthesis (Arufe et al., 2006; Arufe et al., 2009; Jiang et al., 2010; Longmire et al., 2012; Ma et al., 2009).

The primary hurdle preventing the successful differentiation of PSCs into mature thyroid cells has been a lack of knowledge of the pathways that regulate early thyroid embryonic development. “Directed differentiation” of PSCs utilizes sequential exposure of undifferentiated PSCs to a series of growth factor-supplemented medias designed to recapitulate the sequence of developmental milestones that normally occurs during embryonic differentiation. This approach has been employed to successfully produce a wide variety of non-thyroid lineages from PSCs (Murry and Keller, 2008). The derivation of thyroid epithelial cells via directed differentiation remains a compelling goal given the known capacity of thyroid epithelia to self-organize and form follicular structures in vitro (Hilfer et al., 1968; Mallette and Anthony, 1966). Moreover, recent work has demonstrated that thyroid cells generated from embryonic stem cells (ESCs) through forced over-expression of transcription factors also formed follicles in vitro or after transplantation in vivo (Antonica et al., 2012; Ma et al., 2015).

Here we identify combinatorial bone morphogenetic protein (BMP) and fibroblast growth factor (FGF) signaling as necessary to induce thyroid fate in developing anterior foregut endoderm in multiple species from amphibians to humans, and we utilize these pathways to induce thyroid fate in PSC-derived endodermal precursors via the technique of directed differentiation. Thyroid progenitors derived in this way can be sorted to purity using lineage-selective fluorochrome reporters for further culture expansion and subsequent maturation. During maturation in 3D culture the resulting cells form thyroid follicular organoids with an organized monolayered epithelium consistent with thyroid follicles, express genes required for hormone biosynthesis, and can function in vivo following transplantation into hypothyroid mice. Thus we have developed an in vitro system able to reveal the basic developmental mechanisms and gene programs of thyroid cells at early stages of embryonic development and we have produced an inexhaustible source of sortable cells with structural and functional thyroid follicular capacity.

Results

Nkx2-1+ thyroid cells are specified from ESC-derived endodermal precursors upon stimulation with exogenous BMP4 and FGF2

We previously developed a directed differentiation protocol to produce ESC-derived definitive endodermal cells competent to form thyroid and lung epithelial lineages (Longmire et al., 2012). We employed an ESC line (Nkx2-1GFP) carrying a knock-in GFP reporter cDNA targeted to the locus encoding the homeodomain-containing transcription factor Nkx2-1, an essential regulator of thyroid and lung development. Within endoderm, Nkx2-1 is only expressed in lung or thyroid epithelia and is the first known protein induced upon endodermal specification to lung or thyroid fates. Nkx2-1+ cells induced from ESC-derived endoderm using our published (Longmire et al., 2012) 6 factor cocktail (WFKBE+2; Figure. 1A) expressed thyroid epithelial-specific genes, such as thyroglobulin (Tg). However, these cells did not display full thyroid maturation as the genes required for iodine metabolism and thyroid hormone biosynthesis were not robustly expressed (ie the sodium iodine symporter (Nis) and thyroid peroxidase (Tpo)); (Figure S1A and data not shown).

Figure 1. Exogenous BMP4 and FGF2 are sufficient to specify Nkx2-1+ thyroid progenitors from ESC-derived endoderm.

Figure 1

(A) Schematic comparing protocols for directed differentiation of mouse embryonic stem cells (ESCs) into Nkx2-1+ cells. (B) Representative sort gates used to purify Nkx2-1GFP+ cells on Day 14, showing efficiency of specification in each media. (C) Schematic of wildtype (WT) vs targeted Nkx2-1 alleles in the Nkx2-1GFP and Nkx2-1mCherry reporter ESC lines. Boxes represent exons. WT = wild type allele, UTR = untranslated region, IRES = internal ribosomal entry site, pA = polyadenylation site. (D) Representative sort gate on Day 14 used to sort Nkx2-1mCherry+ cell specified with WB+2. (E) mRNA expression by real time RT-PCR comparing the Nkx2-1GFP and Nkx2-1mCherry ESC lines. Bars indicate fold change in gene expression over Day 0 ESCs±SEM 2(−ΔΔCT); n=2 independent experiments). See also Figure S1 and Tables S1 and S2.

Microarray analysis of the global mRNA expression profiles of these ESC-derived Nkx2-1GFP+ endodermal progenitors (Figure 1A and Tables S1 and S2), revealed multiple active signaling pathways including BMP, Wnt, EGF, and FGF signaling. Hence, to identify the essential factors for thyroid cell development we performed sequential withdrawal of each individual factor from the WFKBE+2 cocktail. We found BMP4 and FGF2 to be indispensible for efficient Nkx2-1+ endodermal induction. An induction media composed of BMP4, FGF2 and Wnt3a led to the maximum efficiency of Nkx2-1+ cell specification by Day 14 and sorted Nkx2-1GFP+ cells at this time point, either with or without Wnt3a, were competent to subsequently express thyroid differentiation markers, Tg and Tsh-r (thyroid stimulating hormone receptor), and maturation markers, Nis and Tpo (Figure 1B and Figure S1A). Therefore, BMP4 and FGF2 are sufficient to specify endodermal cells in vitro towards Nkx2-1+ endodermal progenitor cells with thyroid potential.

Nkx2-1 haploinsufficiency affects the differentiation of thyroid epithelial cells, but not the efficiency of Nkx2-1+ endodermal lineage specification

Because the Nkx2-1GFP knock-in reporter in our engineered ESCs replaced one Nkx2-1 allele potentially creating haploinsufficiency (Figure 1C), we sought to test the effects of Nkx2-1 gene dosage on thyroid differentiation in vitro. We compared the Nkx2-1GFP ESC line to a non-haploinsufficient Nkx2-1mCherry ESC line where the mCherry fluorophore is knocked-in 3’ to the coding sequence (Bilodeau et al., 2014) (Figure 1C). The two lines were differentiated in parallel and Nkx2-1GFP+ or Nkx2-1mCherry+ cells were sorted to purity on Day 14 (Figure 1D), at which time the expression of Nkx2-1 mRNA was approximately 50% lower in the haploinsufficient Nkx2-1GFP+ cells than in the non-haploinsufficient Nkx2-1mCherry+ cells (Figure 1E). The expression of both early (Pax8) as well as all mature thyroid markers analyzed (Tg, Tsh-r, Nis, and Tpo) in the outgrowth of sorted Nkx2-1+ cells was diminished in Nkx2-1 haploinsufficient cells (Figure 1E). This was confirmed when the Nkx2-1GFP line was compared to its parental syngeneic, non-haploinsufficient clone (Figure S1C). In contrast we did not detect any difference in the rates of lineage specification of Nkx2-1+ endodermal cells on Day 14, comparing haploinsufficient to normal ESC lines. The percentages and cell numbers of Nkx2-1+ cells specified by Day 14 was similar in the Nkx2-1GFP line and its parental clone (Figure S1B) as well as in the Nkx2-1mCherry line (data not shown). Thus, Nkx2-1 haploinsufficiency does not affect specification of Nkx2-1+ endodermal cells, but does affect subsequent maturation of thyroid epithelial cells. This is supported by observations in mice and humans who have Nkx2-1 haploinsufficiency and hypothyroidism (Krude et al., 2002).

The expression of Pax8 distinguishes a subpopulation within Nkx2-1+ endodermal cells with thyroid potential

Having established a protocol able to produce putative thyroid follicular cells, next we sought to determine whether thyroid and lung lineage specification occurred in distinct endodermal Nkx2-1+ precursors as is thought to occur in vivo. The Pax8 transcription factor is co-expressed with Nkx2-1 only in thyroid epithelial cells, beginning at mouse E8.5 and human E20 in the foregut endoderm within the region of the prospective thyroid primordium (Fagman and Nilsson, 2010; Trueba et al., 2005). In contrast Pax8 is not expressed at any developmental stage in lung epithelial cells. Hence, we developed a bifluorescent reporter system to track thyroid and non-thyroid fates in Nkx2-1+ cells by breeding mice with the GFP reporter targeted to the Nkx2-1 locus (Nkx2-1GFP (Longmire et al., 2012), the Cre recombinase targeted to the Pax8 locus (Pax8Cre (Bouchard et al., 2004)), and a conditional Rosa26 tdTomato reporter (lox-stop-lox-tdTomato). As expected, these mice (hereafter Nkx2-1GFP;Pax8tdTomato trace) co-expressed both GFP and tdTomato reporters only in the thyroid epithelium, whereas GFP without tdTomato was expressed throughout the developing respiratory epithelium (Figures 2A, 2B). As previously reported, GFP was also expressed within developing ectoderm in the ventral forebrain (Longmire et al., 2012), and tdTomato without GFP was expressed in tissues known to express Pax8 including the kidney, hindbrain, midbrain, and inner ear (Figure S2A and as published (Bouchard et al., 2004)). We derived iPSC clones carrying Nkx2-1GFP;Pax8tdTomato trace reporters by reprogramming mouse embryonic fibroblasts (MEFs) isolated from the bifluorescent mouse (Figures 2C, 2D and Figure S2B) and differentiated 6 lines in vitro using our Wnt3a, BMP4, and FGF2 Nkx2-1 induction protocol (Figures 2E, 2F). By Day 14 all clones exhibited induction of the Nkx2-1GFP reporter with kinetics and efficiency similar to our previous report (Longmire et al., 2012) and Pax8tdTomato trace+ cells represented 4.9% (+/−4.7) of all Nkx2-1+ cells. Based on GFP and tdTomato expression, we sorted four distinguishable populations on Day 14 (Figure 2F) and further differentiated each population until Day 25 (Figure 2E and Figure S2C). Analysis of cells at the time of sorting (Day 14) as well as until the full 25 days of differentiation of 3 clones revealed that only the Day14 sorted GFP+ tdTomato+ population gave rise to cells expressing the constellation of 4 transcription factors (Parlato et al., 2004) that together uniquely identifies thyroid epithelia (Nkx2-1, Pax8, Foxe1, and Hhex) and the differentiation marker Tg (Figure 2E). In contrast, this population was depleted of cells competent to express differentiated lung marker genes, Sftpc and Scgb1a1 (data not shown). In separate experiments using ESC lines that lack Pax8 reporters, we used FACS analyses of intracellular protein to verify that Pax8 protein is expressed in approximately 5% of Nkx2-1+ ESC-derived cells (data not shown). We conclude that, as in the developing embryo, the thyroid lineage is specified in a distinct Nkx2-1+ endodermal population.

Figure 2. Thyroid lineage specification in developing embryos and in iPSC-derived endodermal precursors is restricted to cells co-expressing Nkx2-1 and Pax8.

Figure 2

(A) Representative E14.5 mouse embryo carrying Nkx2-1GFP (forebrain and lung) and Pax8tdTomato trace (midbrain, hindbrain, inner ear, and kidney) reporters (B) Dissection of E14.5 mouse embryo anterior foregut showing the developing lung and thyroid (arrow), with Nkx2-1GFP and Pax8tdTomato trace signal overlapping (merge) in the thyroid but not the lung. Note: spherical Nkx2-1+/Pax8-ultimobranchial bodies adjacent to each of the two Nkx2-1+/Pax8+ thyroid lobes. (C) Schematic of induced pluripotent stem cell (iPSC) derivation from MEFs by lentiviral-mediated transduction of the STEMCCA vector. (D) Representative image of the resulting iPSCs, stained for the pluripotency marker Alkaline Phosphatase. (E) Comparison of the four different sorted populations from (F). Cells were purified on Day 14 using the gate indicated and replated for differentiation until Day 25 (Fold change mRNA expression on D25 over D0 by RT-qPCR; 2(−ΔΔCT)). Bars indicate average fold change in gene expression over ESCs±SEM (n=3 independent clones). (F) Representative FACS plot on Day 14 of differentiation of the iPSC line from (C–D). See also Figure S2.

BMP and FGF signaling are necessary and sufficient for thyroid lineage specification from mouse and Xenopus anterior foregut endoderm or ESC-derived definitive endoderm

Having demonstrated that Pax8+/Nkx2-1+ thyroid progenitors were induced from PSC-derived endoderm by addition of Wnt3a, BMP4, and FGF2 (WB+2), we next sought to determine whether each pathway was necessary for thyroid lineage specification. In keeping with prior publications suggesting Wnt is not required for thyroid specification from endoderm in vivo (Goss et al., 2009), we found Nkx2-1GFP+ cells (Figure S1A) or Nkx2-1mCherry+ cells (Figure S3A) derived with BMP4 and FGF2 induction media (Day 6–14) in the absence of Wnt3a displayed undiminished subsequent expression of thyroid markers after outgrowth. While either BMP4 or FGF2 alone induced low numbers of Nkx2-1+ cells with little thyroid competence (Figures 3A, 3B), combinatorial use of BMP4 and FGF2 induced the highest percentage of thyroid competent Nkx2-1+ progenitors, as evidenced by subsequent robust expression of all early and mature thyroid markers including Pax8, Tg, Tsh-r, Tpo, or Nis (Figure 3B and Figure S3A). Consistent with these results, withdrawal of Wnt3a and/or addition of Wnt inhibitors Dkk1 or XAV-939 had no significant impact on the efficiency of lineage specification of Nkx2-1GFP+;Pax8tdTomato trace+ iPSC-derived cells by Day 14, whereas there was a significant decrease in the percent yield of GFP+/tdTomato+ cells when either BMP or FGF2 was removed from the Nkx2-1 induction media (data not shown).

Figure 3. BMP and FGF signaling are necessary for thyroid lineage specification of ESC-derived definitive endoderm.

Figure 3

(A) Representative gates used to sort Nkx2-1mCherry+ cells specified with only FGF2, only BMP4 or FGF2+BMP4 on Day 14. (B) Fold change of mRNA expression on D26 over D0 by RT-qPCR; 2(−ΔΔCT). Cells were purified on Day 14 using the gate indicated in (A) and replated for differentiation until Day 26. (C) Quantitation by flow cytometry of the % of Nkx2-1mCherry+ cells specified on Day 14 using different specification media (Day 6 to 14) containing 4 µM of Dorsomorphin, 20 µM of PD98059 (PD), 10 µM of SB203580 (SB) or 25 µM of LY294002 (LY). Error bars indicate mean±SD (n=3 biological replicates). (D) mRNA expression on Day 26 of the Nkx2-1mCherry ESC line, comparing the specification medias used in Figure 3C. Bars indicate average fold change in gene expression [2(−ΔΔCT)] over D0 ESCs±SD (n=3 biological replicates). 2-way ANOVA (a and c) and 1-way ANOVA (b) *p≤0.05, ****p≤0.0001. See Figure S3.

To determine which downstream BMP and FGF signaling cascades are required for induction of thyroid fate, we supplemented the BMP4 + FGF2 media with specific inhibitors of either SMAD-dependent BMP signaling (Dorsomorphin), MEK1/2-dependent FGF or BMP signaling (PD98059; hereafter PD), PI3 kinase (PI3K)-dependent FGF signaling (LY294002; hereafter LY), or p38-MAPK-dependent BMP signaling (SB203580). We found reduced numbers and percentages of Nkx2-1+ cells in conditions supplemented with Dorsomophin or LY (Figure 3C). Furthermore, when the rare Nkx2-1+ cells that were induced in the presence of these inhibitors were sorted to purity on Day 14 and cultured further to allow completion of the protocol (without inhibitors), the progeny of these Nkx2-1+ cells showed reduced capacity to express mature thyroid markers, such as Tg and Tpo (Figure 3D).

Taken together, our findings in the ESC/iPSC in vitro model system suggested that canonical Wnt signaling may be dispensable for thyroid lineage specification while combinatorial SMAD-dependent BMP signaling together with FGF signaling are required. To test this hypothesis in primary cells of the developing embryo, we employed both the murine and Xenopus models of embryonic foregut organogenesis (Figure 4 and S4). First, in developing Xenopus embryos we found evidence in support of active BMP and FGF signaling in the region of the thyroid primordium by immunostaining embryos for nuclear phospho-SMAD1/5/8, nuclear phospho-ERK1/2, and phospho-AKT ser473 proteins during early thyroid specification (Figure S4A). These phosphorylated signaling effectors were present in the foregut epithelium prior to and during the evagination of Nkx2-1+ cells of the developing thyroid anlage as well as in the surrounding mesenchyme (Figure S4A; stages NF20 and NF33).

Figure 4. ESC models predict the evolutionarily conserved pathways that are necessary and sufficient for thyroid specification in mouse and Xenopus embryos.

Figure 4

(A) BMP signaling blockade abrogates thyroid bud formation in mouse: whole foreguts from E8.0 embryos (6-8 ss) were cultured for 2-3 days in the presence of the BMP antagonist DMH-1. Co-localization of Nkx2-1+ and Pax8+ shows a discrete thyroid bud formed in control cultures (upper left panel). No thyroid bud was apparent when BMP signaling was inhibited (lower left panel). Reduced pSmad1/5 content in the presence of, DMH-1 (right panel; Western blot). (B) In-vivo pharmacological loss-of-function using antagonists of type I BMP receptor (DMH-1) or FGF receptor (SU5402) activity in Xenopus embryos. Whole embryos were cultured in the presence of the antagonists from stage NF13-20 and assayed for the indicated genes at stage NF35. The number of embryos with the displayed phenotype is indicated. (C) FGF and BMP signaling is sufficient to induce thyroid gene expression in dissected Xenopus foregut endoderm explants. Explants were dissected at stage NF15, cultured until NF35, and assayed for the indicated genes. (D) Schematic of inhibition of BMP4 or FGF2 signaling blocking thyroid specification. See also Figure S4.

Next, to assess whether FGF and BMP signaling are required for thyroid specification in vivo, we incubated developing mouse as well as Xenopus embryos in inhibitors of BMP or FGF signaling (Figure 4 and S4). Developing mouse foreguts were isolated by dissection at 6–8 ss (~E8.0) prior to detectable Nkx2-1 expression in the thyroid field and incubated for 2 or 3 days with the BMP inhibitor, DMH-1. DMH-1 caused a marked reduction in phosphorylation of SMAD1/5 (Figure 4A, right panel Western blot) and blocked induction of both Nkx2-1 and Pax8 in the region of the mouse endodermal thyroid primordium (Figure 4A, left panel). Similarly, we incubated developing Xenopus embryos in inhibitors of BMP signaling (DMH-1 or an injected dominant negative BMPR) or FGF signaling (SU5402, PD161570, or an injected dominant negative FGFR), starting just after gastrulation (stage NF13) until stage NF20 (6–7 somite stage; ss). The inhibitors were then removed and embryos allowed to develop until stage NF34 (36 ss); a time by which thyroid and lung lineages are normally specified (Shifley et al., 2012). In situ hybridization for markers of pharyngeal endoderm and thyroid lineage specification foxe1, nkx2-1, pax2, and hhex, revealed that inhibition of either BMP or FGF signaling abrogated induction of all four thyroid markers in the vast majority of Xenopus embryos (Figure 4B and Figure S4B). Furthermore, inhibition of PI3K-AKT with LY also prevented thyroid specification (Figure S4B). Immunostaining of pSmad1/5/8, pERK and pAKT confirmed the efficacy of the inhibitors (Figure S4C). In contrast, when Xenopus embryos were incubated from stage NF13 through NF35 in inhibitors of canonical Wnt signaling (XAV939), RA signaling (BMS493), or VEGF signaling (KRN633), we observed normal nkx2-1 induction in the thyroid primordium (Figure S4B), indicating that Wnt, RA, and VEGF signaling at these developmental stages are dispensable for thyroid specification.

To assess the stage-dependence of these signaling requirements we varied the timing of BMP and FGF loss of function during foregut endoderm development. We observed that early inhibition of BMP or FGF signaling beginning at stage NF13 (analogous to mouse E7.5) blocked induction of nkx2-1 (Figure S4B), whereas inhibition beginning later (at stage NF20; Figure S4D) did not, suggesting that the requirement for BMP and FGF signaling in thyroid lineage specification is restricted to a narrow developmental window between stages NF13-20.

Since our mouse ESC model had predicted that FGF2 and BMP4 were sufficient to induce thyroid lineage specification, we next asked whether exogenous FGF2 and BMP4 were sufficient to induce thyroid development in Xenopus foregut endoderm (Figure 4C). Foregut explants were micro-dissected at stage NF15, prior to thyroid specification and the mesoderm was removed. The foregut endoderm explants were then cultured until stage NF35 either without growth factors or with a combination of FGF2 and BMP4. In situ hybridization revealed that only explants incubated with FGF2 and BMP4 expressed nkx2-1, pax2, and hhex (Figure 4C). We did not detect expression of sftpc in explants from sibling embryos (data not shown) suggesting that the nkx2-1 expression was thyroid and not respiratory epithelium. Taken together these results from Xenopus and mouse embryo models extended our observations made in differentiating mouse ESCs and iPSCs, confirming that FGF and BMP signaling are evolutionarily conserved pathways required for the specification of thyroid fate from developing endoderm both in vitro and in vivo (Figure 4D).

Thyroid stimulating hormone and 3D culture promotes ESC-derived thyroid follicular maturation and organoid formation

Having interrogated the signals required for the induction of thyroid fate, next we focused on augmenting the maturation state of the thyroid epithelial progenitors generated from PSCs, using the Nkx2-1mCherry ESCs. In contrast to lineage specification and early development, the expression of thyroid genes necessary for iodine metabolism, Nis and Tpo, is associated with later gland maturation (Figure 5A) and has been shown in vivo to require TSH receptor activation (Postiglione et al., 2002). Hence we tested the effect of TSH at various developmental stages of our in vitro protocol, either on induction of Nkx2-1+ progenitors (Day 9–12) or during outgrowth of sorted Nkx2-1+ cells (Day 14–22). We found addition of TSH prior to lineage specification (Day 9–12) had no detectable effect on the efficiency of thyroid lineage specification or the competence of thyroid cells to subsequently differentiate (data not shown). In contrast, addition of TSH after lineage specification (Day 14–22) resulted in increased Nis, Tsh-r, and Tpo expression (Figure S5A), but no significant change in the expression of lineage markers Nkx2-1 and Pax8 (data not shown). An additional 3 days of culture maturation in the presence of TSH was accompanied by further augmentation of Nis and Tsh-r expression (Figure S5A).

Figure 5. In vitro development of ESC derived thyroid follicular cells by directed differentiation.

Figure 5

(A) H&E staining of sectioned mouse thyroid tissue (left panel) and schematic of a thyroid follicular cell (right panel). (B) Effect of time in culture on thyroid marker gene expression comparing D26 vs D30. (C–D) Comparison of 3D vs 2D culture on thyroid gene expression. (B–D) Bars indicate average fold change in gene expression by real time RT-PCR over D0 ESCs±SD (n=3 biological replicates, post-natal mouse thyroid tissue n=1). (E–G) Histology on Day 30 of outgrowth and differentiation of Nkx2-1mCherry+ cells sorted on Day 12 as H&E stained paraffin sections (E). Immunofluorescence microscopy of Day 30 ESC-derived follicular-like structures after immunostaining for Nkx2-1 and Tg (F), Nkx2-1 and Pax8 (G). Nuclei are counterstained with DAPI (F–G). Scale bars 100 µm (E) and 10 µm (F, G). Student’s t-test *p≤0.05, **p≤0.01, ***p≤0.001, ****p≤0.0001. See also Figure S5.

We previously published augmented thyroid gene expression in ESC-derived Nkx2-1+ cells when using a thyroid media containing IGF-1, insulin, transferrin, and selenium (Longmire et al., 2012). Therefore, we included these additives, together with TSH to our prior base media of FGF2 + FGF10. This “thyroid outgrowth media” was used for replating the sorted Nkx2-1+ cells on Day 12 (Figure S5B). For the final maturation stages after expansion of sorted cells in “thyroid outgrowth media” we sought to further optimize expression of Nis and Tpo. Since we had observed induction of these markers in response to either DCI+K (Figure 1E) or in response to 1–100 mU/ml of TSH (data not shown), we tested the effect of combinations of each of these factors (Figure 5C; Day 22–26). We found that withdrawal of either cAMP or dexamethasone from the DCI+K cocktail adversely affected maturation (Figure S5C), however substitution of TSH for cAMP in the DCI+K media retained maturation and resulted in the most robust expression of Nis (Figure S5C). Thus, the combination of dexamethasone and TSH was included in the final 4 days of maturation (Day 26–30), resulting in the thyroid directed differentiation protocol summarized in Figure S5B and achieving the derivation from ESCs of a population of cells expressing thyroid marker transcripts at levels equal to (Nkx2-1, Hhex, Foxe1, Tpo) or within 10-fold of (Pax8, Tg, Tshr, Nis) post-natal murine thyroid control tissue (Figure 5B).

To facilitate thyroid follicular structural formation during outgrowth of ESC-derived Nkx2-1+ cells, based on prior reports of Matrigel as a favorable thyroid culture substrate (Martin et al., 1993), we tested the effect of 3D culture in Matrigel (Days 12–30; Figures 5C, 5D). Compared to 2D culture, we found culture of purified Nkx2-1mCherry+ cells in 3D Matrigel conditions resulted in higher expression of thyroid genes by qPCR (Nkx2-1, Pax8, Tg, Tsh-r, Tpo, and Nis; Figure 5C), increased epithelial gene expression (increased E-Cadherin and EpCam expression and reduced mesenchymal markers, Snail1, Twist, and Col1a1; Figure 5D), and triggered the formation of follicular-like clusters of cells (Figure 5E). By Day 30 we observed 130+/−17 follicular clusters emerging from each 65,500 Nkx2-1mCherry+ cells plated on Day 12 in a well of a 24-well plate. By Day 30 in culture, 50% of the progeny of the Nkx2-1mCherry+ sorted cells expressed Pax8 protein by FACS, 46% expressed thyroglobulin, and 96% continued to express Nkx2-1 by either microscopy or FACS for the mCherry reporter (data not shown). Most importantly, phenotyping of the resulting clusters or organoids, revealed them to have the 3D structure and molecular phenotype typical of thyroid follicles, including formation of organized monolayered Nkx2-1+ Pax8+ epithelia surrounding a lumen filled with thyroglobulin (Figures 5F, 5G). These organoids could be further passaged and maintained in these culture conditions with stable or increased expression of Nkx2-1, Pax8, Hhex, Foxe1, Tg, Nis, Tshr, and Tpo until at least Day 52 (Figure S5D). In addition, after incubation in iodinated media for 48 hours the organoids displayed the functional capacity to organify iodine and produce small amounts of T4 hormone bound to TG in vitro at different time points during the differentiation protocol, beginning on Day 30 (Figure S5E).

In vivo function of purified ESC-derived thyroid follicular organoids following transplantation into hypothyroid mouse recipients

The defining functional characteristic of thyroid follicular epithelial cells is their capacity to synthesize and secrete thyroid hormones in vivo. Hence, we sought to test whether ESC-derived Nkx2-1mCherry+ thyroid progenitors sorted on Day 12 of differentiation and further matured into organoids until Day 30 could function in vivo after transplantation, and moreover whether these cells could rescue hypothyroid mice. Syngeneic mouse recipients underwent radioiodine ablation of their native thyroid tissue after 8 days of a low iodine diet in order to induce severe hypothyroidism, evident as significantly reduced circulating plasma T4 and T3 levels and elevated circulating TSH levels (Figures 6A, 6B). Four weeks later, we prepared 1 positive control mouse group (n=5) that had not undergone radioiodine ablation, 1 negative control group that received sham surgery (hereafter “sham”; n=4), and 2 experimental recipient groups that received either undifferentiated (Day 0 ESCs; hereafter “D0-ESCs” n=5) or differentiated ESCs that had been sorted on Day 12 based on Nkx2-1mCherry+ expression and further differentiated until Day 30 in our complete 3D protocol (hereafter D30-ESCs; n=10). Following ESC transplantation with 2.5–3 × 106 cells beneath the left kidney capsule of each mouse recipient we followed all recipient groups for at least 8 weeks post transplantation monitoring circulating T4, T3, and TSH levels, and we employed both MRI and SPECT/CT imaging to screen for in vivo graft growth and function (Figure 6 and Figure S6).

Figure 6. In vivo function of transplanted ESC-derived putative thyroid follicles in hypothyroid mice.

Figure 6

(A) Experiment overview for generation of hypothyroid mice and rescue with in vitro-derived thyroid follicular cells. (B) Recovery of circulating total T4, T3, and TSH 8 weeks after transplantation of thyroid follicular-like cells derived from ESCs (D30-ESCs) but not undifferentiated (D0) ESCs. Data indicate each individual mouse plasma level (dots, squares, and triangles) as well as mean ± SEM for each group. 1-way ANOVA **p≤0.01, ****p≤0.0001. Over the time course shown 4 mice in the D30-ESCs group died. (C) Circulating TSH plotted versus T4 for individual recipients 8 weeks after surgery. (D) 8 weeks after transplantation MRI shows a mass (red arrow) on the left kidney (left panel); SPECT co-registered with MRI shows overlay with the mass (red arrow: right panel). (E) Tc99M SPECT/CT shows signal in the anatomical region of the native thyroid gland (green arrow) in a positive control mouse (left panel) and a hot spot (red arrow) on the left kidney in a D30-ESCs recipient (right panel) 8 weeks after surgery but no signal of the native thyroid tissue in the ablated mouse (right panel). Tc99M uptake in bladder (*), stomach (**), and submandibular gland (black triangle). (F–I) Kidney tissue sections 8 weeks after transplantation of D30-ESCs. (F) H&E staining of transplanted cells in the host kidney expressing nuclear Nkx2-1, nuclear Pax8, and mCherry. (G–H) Immunostaining for Nkx2-1 and Tg (G), or Nkx2-1 and T4 (H). PCNA immunohistochemistry of D30 ESC-derived cells (I). Nuclei are counterstained with DAPI (G, H) or H&E (I). Scale bars 100 µm (F, I) and 10 µm (G, H). See also Figure S6.

Starting 2–4 weeks post transplantation, recipients of D30-ESCs displayed augmented circulating T4 and T3 levels that continued to increase over time until they approached normal levels by 8 weeks after transplantation (Figure 6B and Figure S6A, S6C). Reconstitution of T4 levels was accompanied by amelioration of the elevated plasma TSH levels (Figures 6B, 6C and Figure S6B), further emphasizing that physiologic rescue was occurring. In a repeat experiment some recipients followed for longer periods were observed to display circulating T4 reconstitution at normal levels for at least 29 weeks after transplantation (data not shown). In contrast, we observed no reconstitution of circulating thyroid hormone levels in either sham transplanted controls or recipients of D0-ESCs, excluding the possibility that either residual undifferentiated cells or re-growth of native thyroid tissue was responsible for the return of circulating T4 or the dampening of the TSH response.

To further evaluate the function of the transplanted cells we incorporated imaging with sodium pertechtenate (Tc99M) and SPECT/CT co-registered with MRI imaging. Importantly the uptake of Tc99M into thyroid tissue is mediated by Nis allowing its visualization to be a surrogate for iodine uptake (Zuckier et al., 2004). In recipients of D30-ESCs we observed Tc99M uptake occurring in the region of the kidney capsule grafts (Figures 6D, 6E and Figure S6E). We also observed low levels of Tc99M uptake in the large masses that derived from kidney capsule transplantation of D0-ESCs (Figure S6D), but in no D0-ESC recipient was this signal accompanied by detectable return of circulating T4 levels. As expected, Tc99M uptake was not detected in the anatomical region of the native thyroid tissue of any ablated recipient (Figure 6E and Figure S6E), further suggesting that return of thyroid function in D30-ESC recipients was not due to regrowth of native thyroid tissue.

To confirm the presence of engrafted D30-ESCs in vivo we evaluated the histology of the renal lesions that had been visualized by MRI and SPECT (Figures 6F-6I). As expected, recipients of D0-ESCs developed large teratomas (Figure S6D) without evidence of organized follicles and only rare, if any, cells expressed nuclear Nkx2-1 protein (data not shown). In contrast, no teratomas were observed in recipients of D30-ESCs; instead these recipients displayed persistent grafted cells beneath their kidney capsules organized into many follicular-like structures per high power field that were morphologically consistent with thyroid epithelial follicles (Figure 6F). These structures consisted of mCherry+ epithelial cells expressing nuclear Nkx2-1 and Pax8 proteins by immunostaining, surrounding central lumens filled with colloid-like material that contained Tg and T4 (Figures 6F-6H). A subset of the grafted cells expressed the proliferation marker, PCNA (Figure 6I). We did not detect expression of pro-Sftpc in any grafted cell suggesting that Nkx2-1 expression resulted from thyroid and not respiratory epithelium (Figure S6F).

To assess whether the engrafted D30-ESCs functioned in a regulated fashion we exposed additional D30-ESC transplanted mice, negative sham controls, as well as positive normal controls (n=5 per group) to a high dose of exogenous T3 expected to fully suppress endogenous TSH and hence T4 production. Indeed in all groups administration of T3 resulted in TSH suppression and abrogated T4 secretion (data not shown). To further interrogate TSH responsiveness (Figure S6G), we then injected each mouse with exogenous bovine TSH and three hours later observed an increase in circulating T4 in positive control mice as well as in recipients of D30-ESCs, but not in sham controls, consistent with the regulation of engrafted cells by TSH.

These results indicate that ESC-derived thyroid follicular epithelial cells prepared by directed differentiation in serum-free, feeder-free conditions exhibit in vivo functional potential, including the capacity to rescue hypothyroid mice who otherwise lack native functional thyroid tissue.

Derivation of human thyroid progenitors from ESCs and hypothyroid patient-specific iPSCs

Finally, we sought to adapt our protocol to direct the differentiation of human ESCs and iPSCs into thyroid progenitors. Having observed that thyroid specification was not affected by Nkx2-1 haploinsufficiency in mouse PSCs (Figure 1 and Figure S1), we hypothesized that BMP4 and FGF2 might also generate human thyroid progenitors in normal as well as in haploinsufficient lines. Hence we obtained dermal fibroblasts from three hypothyroid children previously diagnosed with the Brain-Lung-Thyroid syndrome arising from three different respective NKX2-1 coding sequence mutations (Figure S7A, S7B) that were predicted to cause NKX2-1 haploinsufficiency (Hamvas et al., 2013). Using the STEMCCA reprogramming system (Somers et al., 2010) we generated iPSC lines from each child’s fibroblasts (hereafter T1, T3, and T4; Figures 7A, 7B, supp. Figure 7C–7E). Next, we derived definitive endoderm from each iPSC line as well as from control human ESCs (RUES2) and control iPSCs (BU3 and iPS17; Figure 7C and supplemental methods). We employed dual inhibition of TGFβ and BMP signaling as originally described by Green et al. (Green et al., 2011) in order to mimic anterior foregut patterning of each ESCs/iPSCs-derived endodermal population and then tested the capacity of varying concentrations of BMP4 and FGF2 to specify putative thyroid cells. Using the control RUES2 line, we varied the doses of BMP4 and FGF2 to identify a dosing combination that resulted in induction of clusters of NKX2-1+/PAX8+ co-expressing cells (Figure 7D, and Figure S7F). This same dose induced NKX2-1+/PAX8+ cell clusters from human iPSC lines BU3 and iPS17, as well as from all three patient-specific iPSC lines (Figure 7E left panel and Figure S7G). BMP4 or FGF2 alone induced rare, if any, NKX2-1+/PAX8+ co-expressing cells in the T1, T3, T4, and RUES2 lines (data not shown). Furthermore, in BU3, iPS17, and T3 iPSC-derived endodermal cells subjected in parallel to published lung inducing media (Huang et al., 2014) that contains BMP4 (with Chir, FGF10, FGF7, and RA; “CFKBRA”) but lacks FGF2, we observed only induction of NKX2-1+ cells without co-expression of PAX8, suggesting these cells were respiratory rather than thyroid (Figure 7E right panel). As in our mouse ESC/iPSC experiments, NKX2-1+/PAX8+, putative thyroid progenitors induced after 10 days of exposure to BMP4 and FGF2 expressed low to undetectable levels of TPO and TSH-R transcripts (data not shown). However, after further culture in the same thyroid differentiation and maturation conditions used for mouse ESCs these thyroid differentiation markers were upregulated and NKX2-1, PAX8, and TG expression was maintained for at least 42 days in culture (Figure 7F). As expected, without the use of a human NKX2-1 reporter for purification of human thyroid cells from these heterogeneous endodermal cultures, these thyroid markers were detected at significantly lower levels than adult human thyroid control tissue biopsies (data not shown). Taken together these findings support the capacity of combinatorial BMP and FGF signaling to induce thyroid lineage specification from developing endoderm across species, from Xenopus to mice to humans, and enable the generation of patient-specific thyroid progenitors from individuals with genetic or congenital hypothyroidism.

Figure 7. Derivation of human thyroid progenitors from human ESCs/iPSCs and hypothyroid patient-specific iPSCs.

Figure 7

(A) Reprogramming strategy to generate iPSCs lines from a cohort of three children with Brain-Lung-Thyroid Syndrome. Dermal fibroblasts were reprogrammed with the hSTEMCCA-loxP lentivirus and clones were screened for a single integrated STEMCCA copy, which was then excised by transient Cre recombinase. (B) Karyotype of patient T3 (clone T3-2Cr3) indicates the known disease-causing deletion of the NKX2-1 locus (red box): interstitial deletion on chromosome 14 between bands q11.2 and q13. (C) Schematic of directed differentiation protocol of human iPSCs/ESCs into thyroid progenitors. (D) Immunofluorescence microscopy on Day 19 of differentiation of RUES2-derived NKX2-1+/PAX8+ clusters. Overview (left panel, scale bar 500 µm) and zoom (right panel, scale bars 50 µm). (E) Immunostaining for NKX2-1 and PAX8 on Day 17 of differentiation in multiple human iPSCs lines in either thyroid or lung specification media. Nuclei are counterstained with DAPI. Scale bars 50 µm. (F) mRNA expression of thyroid marker genes. Bars indicate average fold change in expression over undifferentiated iPS17 mean±SD (n=3 biological replicates). Student’s t-test *p≤0.05, **p≤0.01. See also Figure S7.

Discussion

Two applications of in vitro model systems, deriving from ESCs and iPSCs, are frequently touted: 1) the engineering of cell-based therapies able to durably ameliorate disease or regenerate tissue after a single transplantation treatment, and 2) the production of an inexhaustible supply of purified cells of defined lineage at identifiable developmental stages for basic mechanistic studies. Here we demonstrate the application of ESCs and iPSCs towards both goals. We have utilized these cells to produce functional, transplantable thyroid follicular organoids able to produce circulating thyroid hormone in vivo at levels that rescue hypothyroid mice, and we have produced thyroid precursor cells at an early developmental stage revealing that combinatorial BMP and FGF signaling together are required for thyroid lineage specification from multi-potent definitive endoderm.

While ESC/iPSC model systems frequently have been used to confirm or validate developmental signaling pathways identified in vivo or in explants, we utilized ESCs and iPSCs to discern a novel role for combinatorial BMP and FGF signaling in thyroid lineage specification, and argue against a role for any individual signaling pathway in isolation. Our findings agree with previous work suggesting a role for FGF signaling in early thyroid development. Our finding of nuclear pERK1/2 protein staining in endodermal Nkx2-1+ thyroid progenitors in Xenopus embryos is consistent with prior demonstrations of active FGF signaling in the developing Xenopus ventral foregut endoderm (Shifley et al., 2012). Moreover, a role for FGF signaling in inducing thyroid cell fate has been previously demonstrated in mouse endodermal explant studies, where addition of exogenous FGF2 prior to thyroid specification induced Nkx2-1 and Tg (Serls et al., 2005). In zebrafish models inhibition of FGF prior to specification leads to absence of a thyroid primordium (Wendl et al., 2007) and mouse genetic models establish that FGF signaling is required for normal early thyroid development, with the adjacent cardiac mesoderm likely being the source of this signal (Celli et al., 1998; Fagman et al., 2007; Fagman and Nilsson, 2010; Kameda et al., 2009; Lania et al., 2009; Vitelli et al., 2002).

In contrast only two prior studies, featuring Chordin mutant and Twisted mutant mice with perturbed BMP signaling, have reported thyroid hypoplasia, raising the possibility of a role for BMP signaling in thyroid development (Bachiller et al., 2003; Petryk et al., 2004). These prior reports could not discern whether the thyroid defects were primary effects due to inhibition of lineage specification, reduced proliferation, or secondary effects due to perturbed cardiac field development. We used our PSC model system, confirmatory Xenopus embryos, and mouse developing foregut endoderm explant models to interrogate the period of thyroid lineage specification from endoderm. We observed that the BMP signal transducers phospho-SMAD1/5/8 are present in the nuclei of endodermal thyroid progenitors at the time of lineage specification, and SMAD-dependent BMP signaling is necessary for induction of the Nkx2-1+ thyroid program. The precise mechanisms by which BMP signaling induces thyroid fate warrants further study to determine whether these are similar to those that induce endodermal Nkx2-1+ lung progenitors

Identification of the sequence of developmental signals that promote thyroid development from endoderm allowed us to perform the generation of functional thyroid follicular-like epithelial cells by the technique of “directed differentiation” of PSCs. In accordance with prior observations from mouse genetic models (Postiglione et al., 2002), we observed TSH signaling to be dispensable for early thyroid specification; however, TSH promoted the maturation of the genetic program necessary for hormone synthesis. As suggested by observations made in patients with hypothyroidism due to NKX2-1 haploinsufficiency (Krude et al., 2002), we also found thyroid differentiation subsequent to lineage specification to be profoundly sensitive to Nkx2-1 haploinsufficiency in our in vitro mouse model system. Furthermore, we were able to apply our protocol to successfully generate NKX2-1+/PAX8+ human thyroid progenitors in iPSCs generated from patients with Brain-Lung-Thyroid disease presenting with hypothyroidism. However, these double positive clusters were rare events possibly due to the NKX2-1 haploinsufficiency that causes the syndrome.

Since the defining feature of thyroid follicular epithelial cells is functional hormone biosynthesis, our transplantation studies provided an important functional assessment of the thyroid organoids we generated by directed differentiation. After transplantation beneath the kidney capsules of hypothyroid mouse recipients, we found evidence for durable, functional engraftment of the putative organoids which after 8 weeks in vivo: 1. retained recognizable structure; 2. continued to express the molecular phenotype of thyroid follicular epithelia; 3. produced circulating T4 and T3; and 4. were regulated by TSH. To prove these grafted organoids derived from Nkx2-1+ endodermal precursors an important aspect of our work involved the sorting to purity of only Nkx2-1mCherry+ endodermal cells prior to transplantation. This sorting algorithm should also help to minimize the chance of teratomas developing in vivo from any residual undifferentiated cells. Indeed, no teratomas were observed arising from the sorted, transplanted cells in any recipient in our studies.

In summary we present a novel in vitro model system based on the directed differentiation of PSCs into thyroid progenitors and mature follicular epithelial organoids for basic developmental studies and in vivo functional thyroid tissue regeneration. This system recapitulates the developmental milestones of early thyroid development, revealing previously unknown mechanisms of thyroid organogenesis. Successful functional transplantation of the resulting cells also suggests the feasibility of a potential future cell-based therapy for hypothyroidism.

Experimental Procedures

Experimental methods are further detailed in the supplement and in downloadable protocols posted online at www.bumc.bu.edu/stem cells.

Animal maintenance

All studies involving mice carrying GFP and tdTomato reporters were approved by the Institutional Animal Care and Use Committee of Boston University School of Medicine. All murine transplantation experiments were approved by Beth Israel Deaconess Medical Center Institutional Animal Care and Use Committee.

Mouse ESC and iPSC lines

Derivation of the Nkx2-1GFP and Nkx2-1mCherry knock-in reporter ESC lines by homologous recombination was previously published and is further detailed in the supplement (Bilodeau et al., 2014; Longmire et al., 2012). The Nkx2-1GFP;Pax8tdTomato trace iPSC reporter line was generated by transduction of MEFs with the dox-inducible Tet-STEMCCA reprogramming lentiviral vector. All ESCs and iPSCs were maintained on fibroblast feeder layers using serum-containing media and LIF (ESGRO Chemicon).

Mouse ESC/iPSC directed differentiation

Cell differentiation was performed in complete serum-free differentiation media (cSFDM, see Extended Experimental Procedures) supplemented with the indicated factors. Briefly, definitive endoderm was induced in embryoid body suspension cultures with 50 ng/ml Activin A (R&D) as previously published (Longmire et al., 2012). Next, 100 ng/ml mNoggin (R&D) and 10 µM SB431542 (Sigma) was applied for 24 hours to generate anterior foregut endoderm. Expression of Nkx2-1 was induced with specification media containing 100 ng/ml mWnt3a, 10 ng/ml mKGF, 10 ng/ml hFGF10, 10 ng/ml mBMP4, 20 ng/ml hEGF, and 250 ng/ml mFGF2 (all from R&D). After Nkx2-1 induction, cells were sorted on Day 12–14 by flow cytometry based on expression of GFP, mCherry, or tdTomato reporters and then replated for either 2D or 3D culture outgrowth as indicated. For 2D outgrowth plated cells were grown in gelatin-coated dishes for 8 more days in cSFDM with 250 ng/ml mFGF2, 100 ng/ml hFGF10 and 100 ng/ml Heparin Sodium Salt. On Day 22, the media was switched to DCI+K maturation media (see Extended Experimental Procedures). For 3D outgrowth cells were plated in pure growth factor reduced Matrigel drops (Corning) in cSFDM supplemented with 250 ng/ml mFGF2, 100 ng/ml hFGF10, 50 ng/ml mIGF-1, 25 ng/ml hEGF (all from R&D), 100 ng/ml Heparin Sodium Salt, 10 µg/ml insulin (Sigma), and 1 mU/ml bTSH (Los Angeles Biomedical Research Institute). On Day 26 media was switched to thyroid maturation media (Ham’s F12 supplemented with 15 mM HEPES, 0.8 mM CaCl2, 100 ng/ml Heparin Sodium Salt, 0.25% BSA, 50 ng/ml mIGF-1, 5 µg/ml insulin, 5 µg/ml ITS, 25 ng/ml hEGF, 50 nM Dexamethasone, and 100 mU/ml bTSH until the harvest day indicated.

Human iPSC derivation

Procurement of all human specimens was approved by the Institutional Review Boards of University of Colorado and Boston University. Three hypothyroid individuals diagnosed with Brain-Lung-Thyroid Syndrome arising from mutations in NKX2-1 were identified and their dermal fibroblasts were reprogrammed with the STEMCCA reprogramming system to generate iPSC lines T1, T3, and T4 (Somers et al., 2010). Control human PSC lines included: BU3 (peripheral blood-derived normal iPSC line; from the CReM of Boston University and Boston Medical Center), iPS17 (fibroblast-derived iPSC) (Crane et al., 2015); (gift of Dr. Brian Davis, University of Texas Health Sciences), and RUES2 human ESCs (gift of Dr. Ali Brivanlou, Rockefeller University). Differentiation of human ESC/iPSC into thyroid cells is detailed in the supplement.

Statistical Analysis

Data are presented as sample means and standard deviations (SD) or standard errors of the mean (SEM), as indicated in the text and in each figure legend. Sample numbers are also indicated in each case. Differences between groups were analyzed by unpaired two-tailed Student’s t-test, 1-way analysis of variance (ANOVA), or 2-way ANOVA with Tukey's multiple-comparison post hoc test as stated in the text or figure legends; p<0.05 was used to indicate significant differences between groups.

Supplementary Material

Highlights.

  • BMP4 and FGF2 are necessary and sufficient to specify the thyroid lineage

  • Pathways regulating thyroid specification are conserved across species

  • Transplantation of ESC-derived thyroid cells regenerates in vivo function

  • Patient-specific thyroid progenitors can be differentiated from human iPSCs

Acknowledgements

We dedicate this work to the memory of our co-first author, Dr. Anita Kurmann, who died in a tragic bicycle accident when this manuscript was in the final stages of formatting. Dr. Kurmann was an accomplished endocrine and general surgeon from Switzerland who sought basic science post-doctoral research fellowship training in Boston in the co-mentored setting of the Hollenberg and Kotton laboratories. She was intelligent, well-read, kind, humble, and was tirelessly committed to her patients, her thyroid research, her family and her colleagues who miss her dearly.

The authors wish to thank the members of the Kotton and Hollenberg labs for insightful discussions. We thank Ron Gluck and Elizabeth Brody-Gluck for project ignition, and Dr. Meinrad Busslinger for kindly sharing the Pax8cre mouse. We thank Brian R. Tilton and Patrick Autissier of the Boston University Flow Cytometry Core Facility for technical assistance, and Amulya Iyer of the Kotton lab and Marianne James of the iPS core at the CReM for maintenance and characterization of patient-specific iPSCs. A.A.K. was supported by a grant from the Swiss National Science Foundation (SNSF) PBBS P3-146612, L.I. by R01 HL111574 and an award co-sponsored by the chILD Foundation/American Thoracic Society, M.B. by a postdoctoral fellowships from CIHR/Canadian Lung Association/GSK partnership and from the FRQS, A.M.Z. by R01 HL114898, and D.N.K. by R01 HL095993, R01 HL108678, R01 HL122442, and U01HL110-967.

Footnotes

Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

The authors have no conflict of interest to disclose

Supplemental information

Supplemental information includes Extended Experimental Procedures, seven supplementary figures, and two supplementary tables.

Author contributions

A.A.K., M.S., A.N.H. and D.N.K. conceived the work, designed experiments, and wrote the manuscript. S.A.R. and A.M.Z. performed Xenopus experiments and wrote the manuscript. J.M.S. and S.L. performed mouse foregut explant experiments and wrote the manuscript. M.B. and J.R. designed the Nkx2-1mCherry mouse model. S.U. performed the Tc99M SPECT/CT, and MRI imaging experiments. J.C.J., R.R.D., and F.H. procured human cells and developed the human ESC/iPSC model. M.M., F.H., I.A., and L.I. provided critical intellectual input. All authors reviewed, edited and approved the manuscript.

References

  1. Antonica F, Kasprzyk DF, Opitz R, Iacovino M, Liao XH, Dumitrescu AM, Refetoff S, Peremans K, Manto M, Kyba M, et al. Generation of functional thyroid from embryonic stem cells. Nature. 2012;491:66–71. doi: 10.1038/nature11525. [DOI] [PMC free article] [PubMed] [Google Scholar]
  2. Arufe MC, Lu M, Kubo A, Keller G, Davies TF, Lin RY. Directed differentiation of mouse embryonic stem cells into thyroid follicular cells. Endocrinology. 2006;147:3007–3015. doi: 10.1210/en.2005-1239. [DOI] [PMC free article] [PubMed] [Google Scholar]
  3. Arufe MC, Lu M, Lin RY. Differentiation of murine embryonic stem cells to thyrocytes requires insulin and insulin-like growth factor-1. Biochem Biophys Res Commun. 2009;381:264–270. doi: 10.1016/j.bbrc.2009.02.035. [DOI] [PMC free article] [PubMed] [Google Scholar]
  4. Bachiller D, Klingensmith J, Shneyder N, Tran U, Anderson R, Rossant J, De Robertis EM. The role of chordin/Bmp signals in mammalian pharyngeal development and DiGeorge syndrome. Development. 2003;130:3567–3578. doi: 10.1242/dev.00581. [DOI] [PubMed] [Google Scholar]
  5. Bilodeau M, Shojaie S, Ackerley C, Post M, Rossant J. Identification of a proximal progenitor population from murine fetal lungs with clonogenic and multilineage differentiation potential. Stem cell reports. 2014;3:634–649. doi: 10.1016/j.stemcr.2014.07.010. [DOI] [PMC free article] [PubMed] [Google Scholar]
  6. Bouchard M, Souabni A, Busslinger M. Tissue-specific expression of cre recombinase from the Pax8 locus. Genesis. 2004;38:105–109. doi: 10.1002/gene.20008. [DOI] [PubMed] [Google Scholar]
  7. Celli G, LaRochelle WJ, Mackem S, Sharp R, Merlino G. Soluble dominant-negative receptor uncovers essential roles for fibroblast growth factors in multi-organ induction and patterning. EMBO J. 1998;17:1642–1655. doi: 10.1093/emboj/17.6.1642. [DOI] [PMC free article] [PubMed] [Google Scholar]
  8. Crane AM, Kramer P, Bui JH, Chung WJ, Li XS, Gonzalez-Garay ML, Hawkins F, Liao W, Mora D, Choi S, et al. Targeted Correction and Restored Function of the CFTR Gene in Cystic Fibrosis Induced Pluripotent Stem Cells. Stem cell reports. 2015;4:569–577. doi: 10.1016/j.stemcr.2015.02.005. [DOI] [PMC free article] [PubMed] [Google Scholar]
  9. Fagman H, Liao J, Westerlund J, Andersson L, Morrow BE, Nilsson M. The 22q11 deletion syndrome candidate gene Tbx1 determines thyroid size and positioning. Hum Mol Genet. 2007;16:276–285. doi: 10.1093/hmg/ddl455. [DOI] [PubMed] [Google Scholar]
  10. Fagman H, Nilsson M. Morphogenesis of the thyroid gland. Molecular and cellular endocrinology. 2010;323:35–54. doi: 10.1016/j.mce.2009.12.008. [DOI] [PubMed] [Google Scholar]
  11. Goss AM, Tian Y, Tsukiyama T, Cohen ED, Zhou D, Lu MM, Yamaguchi TP, Morrisey EE. Wnt2/2b and beta-catenin signaling are necessary and sufficient to specify lung progenitors in the foregut. Dev Cell. 2009;17:290–298. doi: 10.1016/j.devcel.2009.06.005. [DOI] [PMC free article] [PubMed] [Google Scholar]
  12. Green MD, Chen A, Nostro MC, d’Souza SL, Schaniel C, Lemischka IR, Gouon-Evans V, Keller G, Snoeck HW. Generation of anterior foregut endoderm from human embryonic and induced pluripotent stem cells. Nat Biotechnol. 2011;29:267–272. doi: 10.1038/nbt.1788. [DOI] [PMC free article] [PubMed] [Google Scholar]
  13. Hamvas A, Deterding RR, Wert SE, White FV, Dishop MK, Alfano DN, Halbower AC, Planer B, Stephan MJ, Uchida DA, et al. Heterogeneous Pulmonary Phenotypes Associated With Mutations in the Thyroid Transcription Factor Gene NKX2-1. Chest. 2013;144:794–804. doi: 10.1378/chest.12-2502. [DOI] [PMC free article] [PubMed] [Google Scholar]
  14. Hilfer SR, Iszard LB, Hilfer EK. Follicle formation in the embryonic chick thyroid. II. Reorganization after dissociation. Zeitschrift fur Zellforschung und mikroskopische Anatomie. 1968;92:256–269. doi: 10.1007/BF00335651. [DOI] [PubMed] [Google Scholar]
  15. Huang SX, Islam MN, O’Neill J, Hu Z, Yang YG, Chen YW, Mumau M, Green MD, Vunjak-Novakovic G, Bhattacharya J, et al. Efficient generation of lung and airway epithelial cells from human pluripotent stem cells. Nat Biotechnol. 2014;32:84–91. doi: 10.1038/nbt.2754. [DOI] [PMC free article] [PubMed] [Google Scholar]
  16. Jiang N, Hu Y, Liu X, Wu Y, Zhang H, Chen G, Liang J, Lu X, Liu S. Differentiation of E14 mouse embryonic stem cells into thyrocytes in vitro. Thyroid. 2010;20:77–84. doi: 10.1089/thy.2008.0291. [DOI] [PubMed] [Google Scholar]
  17. Kameda Y, Ito M, Nishimaki T, Gotoh N. FRS2alpha is required for the separation, migration, and survival of pharyngeal-endoderm derived organs including thyroid, ultimobranchial body, parathyroid, and thymus. Dev Dyn. 2009;238:503–513. doi: 10.1002/dvdy.21867. [DOI] [PubMed] [Google Scholar]
  18. Krude H, Schutz B, Biebermann H, von Moers A, Schnabel D, Neitzel H, Tonnies H, Weise D, Lafferty A, Schwarz S, et al. Choreoathetosis, hypothyroidism, and pulmonary alterations due to human NKX2-1 haploinsufficiency. J Clin Invest. 2002;109:475–480. doi: 10.1172/JCI14341. [DOI] [PMC free article] [PubMed] [Google Scholar]
  19. Lancaster MA, Knoblich JA. Organogenesis in a dish: modeling development and disease using organoid technologies. Science. 2014;345:1247125. doi: 10.1126/science.1247125. [DOI] [PubMed] [Google Scholar]
  20. Lania G, Zhang Z, Huynh T, Caprio C, Moon AM, Vitelli F, Baldini A. Early thyroid development requires a Tbx1-Fgf8 pathway. Dev Biol. 2009;328:109–117. doi: 10.1016/j.ydbio.2009.01.014. [DOI] [PMC free article] [PubMed] [Google Scholar]
  21. Longmire TA, Ikonomou L, Hawkins F, Christodoulou C, Cao Y, Jean JC, Kwok LW, Mou H, Rajagopal J, Shen SS, et al. Efficient derivation of purified lung and thyroid progenitors from embryonic stem cells. Cell Stem Cell. 2012;10:398–411. doi: 10.1016/j.stem.2012.01.019. [DOI] [PMC free article] [PubMed] [Google Scholar]
  22. Ma R, Latif R, Davies TF. Thyrotropin-independent induction of thyroid endoderm from embryonic stem cells by activin A. Endocrinology. 2009;150:1970–1975. doi: 10.1210/en.2008-1374. [DOI] [PMC free article] [PubMed] [Google Scholar]
  23. Ma R, Latif R, Davies TF. Thyroid follicle formation and thyroglobulin expression in multipotent endodermal stem cells. Thyroid. 2013;23:385–391. doi: 10.1089/thy.2012.0644. [DOI] [PMC free article] [PubMed] [Google Scholar]
  24. Ma R, Latif R, Davies TF. Human embryonic stem cells form functional thyroid follicles. Thyroid. 2015 doi: 10.1089/thy.2014.0537. [DOI] [PMC free article] [PubMed] [Google Scholar]
  25. Mallette JM, Anthony A. Growth in culture of trypsin dissociated thyroid cells from adult rats. Experimental cell research. 1966;41:642–651. doi: 10.1016/s0014-4827(66)80115-5. [DOI] [PubMed] [Google Scholar]
  26. Martin A, Valentine M, Unger P, Lichtenstein C, Schwartz AE, Friedman EW, Shultz LD, Davies TF. Preservation of functioning human thyroid organoids in the scid mouse: 1. System characterization. The Journal of clinical endocrinology and metabolism. 1993;77:305–310. doi: 10.1210/jcem.77.2.8345031. [DOI] [PubMed] [Google Scholar]
  27. McCracken KW, Cata EM, Crawford CM, Sinagoga KL, Schumacher M, Rockich BE, Tsai YH, Mayhew CN, Spence JR, Zavros Y, et al. Modelling human development and disease in pluripotent stem-cell-derived gastric organoids. Nature. 2014 doi: 10.1038/nature13863. [DOI] [PMC free article] [PubMed] [Google Scholar]
  28. Murry CE, Keller G. Differentiation of embryonic stem cells to clinically relevant populations: lessons from embryonic development. Cell. 2008;132:661–680. doi: 10.1016/j.cell.2008.02.008. [DOI] [PubMed] [Google Scholar]
  29. Pagliuca FW, Millman JR, Gurtler M, Segel M, Van Dervort A, Ryu JH, Peterson QP, Greiner D, Melton DA. Generation of functional human pancreatic beta cells in vitro. Cell. 2014;159:428–439. doi: 10.1016/j.cell.2014.09.040. [DOI] [PMC free article] [PubMed] [Google Scholar]
  30. Parlato R, Rosica A, Rodriguez-Mallon A, Affuso A, Postiglione MP, Arra C, Mansouri A, Kimura S, Di LR, De FM. An integrated regulatory network controlling survival and migration in thyroid organogenesis. DevBiol. 2004;276:464–475. doi: 10.1016/j.ydbio.2004.08.048. [DOI] [PubMed] [Google Scholar]
  31. Petryk A, Anderson RM, Jarcho MP, Leaf I, Carlson CS, Klingensmith J, Shawlot W, O’Connor MB. The mammalian twisted gastrulation gene functions in foregut and craniofacial development. Dev Biol. 2004;267:374–386. doi: 10.1016/j.ydbio.2003.11.015. [DOI] [PubMed] [Google Scholar]
  32. Postiglione MP, Parlato R, Rodriguez-Mallon A, Rosica A, Mithbaokar P, Maresca M, Marians RC, Davies TF, Zannini MS, De Felice M, et al. Role of the thyroid-stimulating hormone receptor signaling in development and differentiation of the thyroid gland. Proc Natl Acad Sci U S A. 2002;99:15462–15467. doi: 10.1073/pnas.242328999. [DOI] [PMC free article] [PubMed] [Google Scholar]
  33. Serls AE, Doherty S, Parvatiyar P, Wells JM, Deutsch GH. Different thresholds of fibroblast growth factors pattern the ventral foregut into liver and lung. Development. 2005;132:35–47. doi: 10.1242/dev.01570. [DOI] [PubMed] [Google Scholar]
  34. Shifley ET, Kenny AP, Rankin SA, Zorn AM. Prolonged FGF signaling is necessary for lung and liver induction in Xenopus. BMC Dev Biol. 2012;12:27. doi: 10.1186/1471-213X-12-27. [DOI] [PMC free article] [PubMed] [Google Scholar]
  35. Somers A, Jean JC, Sommer CA, Omari A, Ford CC, Mills JA, Ying L, Sommer AG, Jean JM, Smith BW, et al. Generation of Transgene-Free Lung Disease-Specific Human iPS Cells Using a Single Excisable Lentiviral Stem Cell Cassette. Stem Cells. 2010;28:1728–1740. doi: 10.1002/stem.495. [DOI] [PMC free article] [PubMed] [Google Scholar]
  36. Trueba SS, Auge J, Mattei G, Etchevers H, Martinovic J, Czernichow P, Vekemans M, Polak M, Attie-Bitach T. PAX8, TITF1, and FOXE1 gene expression patterns during human development: new insights into human thyroid development and thyroid dysgenesis-associated malformations. The Journal of clinical endocrinology and metabolism. 2005;90:455–462. doi: 10.1210/jc.2004-1358. [DOI] [PubMed] [Google Scholar]
  37. Vitelli F, Taddei I, Morishima M, Meyers EN, Lindsay EA, Baldini A. A genetic link between Tbx1 and fibroblast growth factor signaling. Development. 2002;129:4605–4611. doi: 10.1242/dev.129.19.4605. [DOI] [PubMed] [Google Scholar]
  38. Wendl T, Adzic D, Schoenebeck JJ, Scholpp S, Brand M, Yelon D, Rohr KB. Early developmental specification of the thyroid gland depends on han-expressing surrounding tissue and on FGF signals. Development. 2007;134:2871–2879. doi: 10.1242/dev.02872. [DOI] [PubMed] [Google Scholar]
  39. Zuckier LS, Dohan O, Li Y, Chang CJ, Carrasco N, Dadachova E. Kinetics of perrhenate uptake and comparative biodistribution of perrhenate, pertechnetate, and iodide by NaI symporter-expressing tissues in vivo. Journal of nuclear medicine : official publication, Society of Nuclear Medicine. 2004;45:500–507. [PubMed] [Google Scholar]

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

RESOURCES