Thyroid Hormone Promotes Fetal Neurogenesis

Abstract

Maternal low thyroxine (T4) serum levels during the first trimester of pregnancy correlate with cerebral cortex volume and mental development of the progeny, but why neural cells during early fetal brain development are vulnerable to maternal T4 levels remains unknown. In this study, using iPSCs obtained from a boy with a loss-of-function mutation in MCT8-a transporter previously identified as critical for thyroid hormone uptake and action in neural cells-we demonstrate that thyroid hormones induce transcriptional changes that promote the progression of human neural precursor cells along the dorsal projection trajectory. Consistent with these findings, single-cell, spatial, and bulk transcriptomics from MCT8-deficient cerebral organoids and cultures of human neural precursor cells underscore the necessity for optimal thyroid hormone levels for these cells to differentiate into neurons. The controlled intracellular activation of T4 signaling occurs through the transient expression of the enzyme type 2 deiodinase, which converts T4 into its active form, T3, alongside the coordinated expression of thyroid hormone nuclear receptors. The intracellular activation of T4 in NPCs results in transcriptional changes important for their division mode and cell cycle progression. Thus, T4 is essential for fetal neurogenesis, highlighting the importance of adequate treatment for mothers with hypothyroidism.

Keywords: DIO2; MCT8; T3 signaling; brain; deiodinase; thyroid.

Figure 1.

A. Schematic representation of generating D50 MCT8-COs and obtaining single-cell suspension for sc-RNA seq analysis. B. UMAP plot showing the cell clusters identified by principal component analysis in the indicated COs C. UMAP plots showing the distribution of markers for NPCs and neurons. D. Dot plot showing the relative expression levels of the gene markers used to identify each cell population. E. UMAP plot showing the cell types identified by manual curation in control COs. F. Dot plot showing the relative expression levels of SLC16A2, THRB, and THRA in each cell population of the dorsal projection trajectory. G. Interpretation of sequential gene expression changes during dorsal projection trajectory progression. H. UMAP plot showing the cell types identified by manual curation in MCT8-COs. I. Histogram of the relative number of cells in control and MCT8-COs. J. Volcano plots showing the distribution of differentially expressed genes in MCT8-deficient vs. control COs. iPSCs, induced pluripotent stem cells; COs, cortical organoids; NPCs, neural precursor cells; IPCs, intermediate precursor cells; DPT, dorsal projection trajectory; UMAP, Uniform Manifold Approximation and Projection. The χ² test for multiple comparisons and pairwise cell proportions. **P < 0.01; ***P < 0.0001. Differentially expressed genes thresholds: p-value < 0.05, and Average Log₂ Fold-Change of 0.26.

Figure 2.

A. Schematic representation of obtaining spatial transcriptomic data. B. Image showing the cell segmentation on the indicated COs, based on transcript localization data. Scale bar: 100 μm. C. Heatmap showing the relative expression levels of the genes used to identify the indicated neural cell types. D. Heatmap of the differentially expressed genes in the indicated neural cells between control vs. MCT8-COs. E. Same as in B, except the cells are classified into the indicated types. F. Distribution of distances between the indicated cells and groups. G. Histogram of the distribution of cell densities. Considering a radius of 100 μm, if only one cell was in contact with another cell, it was considered as sparse density. If one cell was in contact with five or more cells, it was considered a dense density. H. Violin plots show the expression levels of the indicated genes, considering the indicated cellular densities. I. Ligand-receptor interaction between the indicated neural cells in Control and MCT8-deficient COs. Differentially expressed genes threshold: p-value < 0.05.

Figure 3.

A. Schematic representation of generating NPCs from iPSCs and the subsequent differentiation of the NPCs into neurons. B. Brightfield and confocal fluorescence images showing iPSCs-derived NPCs stained for SOX2 (magenta), NESTIN (yellow), MCT8 (magenta), and Dapi (blue; nuclear). C. Relative mRNA levels of the indicated genes in Control and MCT8 NPCs after 1, 6, and 12 days of neurodifferentiation. D. Brightfield images showing NPCs after six days of neurodifferentiation. E. Principal component plot illustrating differences between control and MCT8-deficient NPCs. F. Volcano plots showing the distribution of differentially expressed genes in control vs. MCT8-NPCs; each point represents the average of 5 control and 5 MCT8-deficient samples of pooled NPCs for each transcript. G. Heatmap depicting the top 20 differentially expressed genes related to neurogenesis between control vs. MCT8-NPCs identified by bulk-RNA seq. H. Venn comparison of differentially expressed genes belonging to the gene set related to neurogenesis between control vs. MCT8-NPCs identified by bulk-RNA seq. (blue) and sc-RNA seq analysis (green) and between control vs. MCT8-IPCs identified by sc-RNA seq analysis (purple); common differentially expressed genes were identified (grey box). I. Brightfield images of control and MCT8-NPCs after twelve days of neurodifferentiation; TUJ1 staining in green and Dapi (blue; nuclear). J. The upper two panels are MCT8 staining in red, TUJ1 in green, and Dapi (blue); the lower panels are RBFOX3 staining in green, NEUROD1 in red, and Dapi in blue. K. MCT8-NPCs after being treated with 60nM T3 during the twelve days of neurodifferentiation. TUJ1 staining is green, and SOX2 is red. L. SOX2 staining in red and Dapi in blue on the indicated cells and treatments. Scale bars: B: 50 μm; D, I, J: 100 μm. Differentially expressed genes thresholds: p-value < 0.05, and Average Log₂ Fold-Change of 1.5 in the Partek Flow platform. Expression values are mean ± SD of n = 3–6 RNA samples, each of them consisting of 2 pooled 6-well plates of NPCs from either control or MCT8-NPCs; Two-tailed Student’s test for comparing D2 deiodination and relative mRNA expression between D1, D6 and D12 of neurodifferentiation; *P < 0.05, **P < 0.01, ***P < 0.001.

Figure 4.

A. Schematic representation of treating iPSCs-derived NPCs with radioactive T4^I125 to measure T3^I125 production. B. Representative chromatograms of the medium after control and MCT8-NPCs were incubated with T4^I125 for 24 hours. C. Quantitation of the DIO2 deiodination in control and MCT8 NPCs; n = 5 DIO2 assays. D. Volcano plots showing the distribution of differentially expressed genes in control + 1nM T4 vs. control NPCs; each point represents the average of five control + 1nM T4 and five control samples of pooled NPCs for each transcript. E. Interpretation of the findings in A-D. F. Schematic of the generation and timing of COs generation, starting with iPSCs to a culture of embryoid bodies, followed by neural induction, neuroepithelial bud expansion, and maturation. G. Quantitation of DIO2 deiodination in control COs during their first 20 days in culture. n = 4 DIO2 assays per timepoint, each consisting of 4 pooled COs from control COs. H. Relative SOX2 mRNA levels in control COs during their first 20 days in culture. Expression values are mean ± SD of n = 3–6 RNA samples, each of them consisting of 4 pooled COs from control COs; I. Schematic representation of the experiment: COs are treated with 1nM T4 from D7 to D50 and then dissociated into a single-cell suspension for sc-RNA seq. J. UMAP plot showing the cell types identified. K. Histogram of the relative number of cells in T4-COs and control COs. L. Histograms of the relative number of cells in clusters of NPCs. The identification number of each cell cluster is indicated at the bottom right corner of each rectangle. M. Volcano plots showing the distribution of differentially expressed genes in T4-CO vs. control COs. H. Gene set enrichment analysis reveals gene ontology terms enriched in T4_TX COs. O. Histograms of the relative number of cells undergoing the indicated cell cycle phase in clusters one and nine of control and NPCs. Two-tailed Student’s test for comparing DIO2 deiodination in iPSC-derived NPCs; ***P < 0.001. Differentially expressed genes thresholds: p-value < 0.05, and Average Log₂ Fold-Change of 0.26.

Figure 5.

A. Schematic representation of the protocol to isolate mNPCs, propagate them in culture, and differentiate into neurons. B, C. Brightfield images showing representative neurospheres containing mNPCs (B) and the mNPCs cultured in collagen-coated plasticware (C). D-F. Confocal images showing mNPCs expressing Nestin (green) and Sox2 (red) (D), and Mct8 (red) (E). The insets in E depict two dividing mNPCs that exhibited higher intensity of Mct8 immunofluorescence. F. Quantitation of the DIO2 deiodination in mNPCs; n = 6 DIO2 assays. G. Brightfield images showing representative mNPCs after two days of neurodifferentiation. H-J. Relative mRNA levels of the indicated genes in mNPCs under the indicated conditions; n = 4. K. Brightfield images showing representative mNPCs after four days of neurodifferentiation. Note the neuronal process extension. L. Tuj1+ cells under the indicated conditions and the quantitation of the percentage of Tuj1+ cells. Values are the mean ± SD of 5 replicates. Scale bars: B: 300 μm, C-L: 25 μm. Two-tailed Student’s test for comparing DIO2 deiodination in iPSC-derived NPCs; **P < 0.01; ***P < 0.001.