# Single-cell analysis reveals the cellular and transcriptional heterogeneity of the normal pediatric thyroid

**Authors:** Erin R. Reichenberger, Nicholas E. Bambach, Zachary Spangler, Julio C. Ricarte-Filho, Kyle Hinkle, Amber Isaza, Tricia Bhatti, Andrew J. Bauer, Aime T. Franco

PMC · DOI: 10.3389/fendo.2026.1748458 · Frontiers in Endocrinology · 2026-01-29

## TL;DR

This study uses single-cell RNA sequencing to explore the diversity of cell types and functions in the normal pediatric thyroid.

## Contribution

The first comprehensive analysis of cellular and transcriptional heterogeneity in the normal pediatric thyroid without malignancy.

## Key findings

- Seven major cell types were identified, including thyrocytes, endothelial cells, and immune cells.
- Thyrocytes showed significant heterogeneity, with subsets linked to hormone synthesis and transport.
- An immune subpopulation was detected in one sample, highlighting transcriptional differences without visible histological changes.

## Abstract

To enhance the understanding of cellular heterogeneity within the pediatric thyroid, single-nuclei RNA sequencing was used to recover 38,069 non-pathogenic cells from thyroid tissue of three female pediatric/young adult patients. The recovered cells were analyzed using the SWANS (Single Entity Workflow ANalysiS pipeline (version 1.0). Analysis revealed seven major cell types: thyrocytes, endothelial cells, fibroblasts, C cells, T cells, B cells, and myeloid cells. Thyrocytes were the most prominent and heterogeneous cell type. Initially, two dominant thyrocyte subsets were identified based on transcriptional activity, which were subsequently subdivided into seven subclusters. Differentially expressed genes within each cluster support distinct cellular functions, including a metabolically active subset which may be involved in hormone synthesis and a subset involved in the transport of thyroid hormone into circulation. We identified an immune subpopulation originating predominantly from a single sample that was histologically and morphologically similar to the other two samples. This supports that transcriptional changes can be detected and used to identify populations of cells, even in the absence of histologically observable changes. This characterization represents the first comprehensive portraiture of pediatric thyroid gland cells and the first description of normal patient thyrocyte and stromal cell heterogeneity in the absence of adjacent malignancy.

## Linked entities

- **Species:** Homo sapiens (taxon 9606)

## Full-text entities

- **Genes:** MS4A1 (membrane spanning 4-domains A1) [NCBI Gene 931] {aka B1, Bp35, CD20, CVID5, FMC7, LEU-16}, LYVE1 (lymphatic vessel endothelial hyaluronan receptor 1) [NCBI Gene 10894] {aka CRSBP-1, HAR, LYVE-1, XLKD1}, COX1 (cytochrome c oxidase subunit I) [NCBI Gene 4512] {aka COI, MTCO1}, SRGN (serglycin) [NCBI Gene 5552] {aka PPG, PRG, PRG1}, FRY (FRY microtubule binding protein) [NCBI Gene 10129] {aka 13CDNA73, 214K23.2, C13orf14, CG003, bA207N4.2, bA37E23.1}, NKAIN2 (sodium/potassium transporting ATPase interacting 2) [NCBI Gene 154215] {aka FAM77B, NKAIP2, TCBA, TCBA1}, COX8A (cytochrome c oxidase subunit 8A) [NCBI Gene 1351] {aka COX, COX8, COX8-2, COX8L, MC4DN15, VIII}, CD74 (CD74 molecule) [NCBI Gene 972] {aka CLIP, DHLAG, HLADG, II, Ia-GAMMA, p33}, RPS4X (ribosomal protein S4 X-linked) [NCBI Gene 6191] {aka CCG2, DXS306, RPS4, S4, SCAR, SCR10}, COL8A1 (collagen type VIII alpha 1 chain) [NCBI Gene 1295] {aka C3orf7}, COX4I1 (cytochrome c oxidase subunit 4I1) [NCBI Gene 1327] {aka COX IV-1, COX4, COX4-1, COXIV, COXIV-1, MC4DN16}, NDUFS5 (NADH:ubiquinone oxidoreductase subunit S5) [NCBI Gene 4725] {aka CI-15k, CI15K}, RPS21 (ribosomal protein S21) [NCBI Gene 6227] {aka HLDF, S21, eS21}, EPB41L3 (erythrocyte membrane protein band 4.1 like 3) [NCBI Gene 23136] {aka 4.1B, DAL-1, DAL1}, STK32A (serine/threonine kinase 32A) [NCBI Gene 202374] {aka YANK1}, EEF1G (eukaryotic translation elongation factor 1 gamma) [NCBI Gene 1937] {aka EF1G, GIG35}, BANK1 (B cell scaffold protein with ankyrin repeats 1) [NCBI Gene 55024] {aka BANK}, ND2 (NADH dehydrogenase subunit 2) [NCBI Gene 4536] {aka MTND2}, CYTB (cytochrome b) [NCBI Gene 4519] {aka MTCYB}, MAP2 (microtubule associated protein 2) [NCBI Gene 4133] {aka MAP-2, MAP2A, MAP2B, MAP2C}, SLIT2 (slit guidance ligand 2) [NCBI Gene 9353] {aka SLIL3, Slit-2}, NKX2-1 (NK2 homeobox 1) [NCBI Gene 7080] {aka BCH, BHC, NK-2, NKX2.1, NKX2A, NMTC1}, EEF1A1 (eukaryotic translation elongation factor 1 alpha 1) [NCBI Gene 1915] {aka CCS-3, CCS3, EE1A1, EEF-1, EEF1A, EF-Tu}, RASSF8 (Ras association domain family member 8) [NCBI Gene 11228] {aka C12orf2, HOJ1}, NDUFC2 (NADH:ubiquinone oxidoreductase subunit C2) [NCBI Gene 4718] {aka B14.5b, CI-B14.5b, HLC-1, MC1DN36, NADHDH2}, SLA (Src like adaptor) [NCBI Gene 6503] {aka SLA1, SLAP}, FOXE1 (forkhead box E1) [NCBI Gene 2304] {aka BAMLAZ, FKHL15, FOXE2, HFKH4, HFKL5, NMTC4}, PAX8 (paired box 8) [NCBI Gene 7849] {aka PAX-8}, CALD1 (caldesmon 1) [NCBI Gene 800] {aka CDM, H-CAD, HCAD, L-CAD, LCAD, NAG22}, NDUFA3 (NADH:ubiquinone oxidoreductase subunit A3) [NCBI Gene 4696] {aka B9, CI-B9}, EEF2 (eukaryotic translation elongation factor 2) [NCBI Gene 1938] {aka EEF-2, EF-2, EF2, SCA26}, COL12A1 (collagen type XII alpha 1 chain) [NCBI Gene 1303] {aka BA209D8.1, BTHLM2, COL12A1L, DJ234P15.1, EDSMYP, UCMD2}, RPL29 (ribosomal protein L29) [NCBI Gene 6159] {aka HIP, HUMRPL29, L29, RPL29P10, RPL29_3_370, eL29}, RAB6A (RAB6A, member RAS oncogene family) [NCBI Gene 5870] {aka RAB6}, CDH2 (cadherin 2) [NCBI Gene 1000] {aka ACOGS, ADHD8, ARVD14, CD325, CDHN, CDw325}, DIO1 (iodothyronine deiodinase 1) [NCBI Gene 1733] {aka 5DI, THMA2, TXDI1}, RAB31 (RAB31, member RAS oncogene family) [NCBI Gene 11031] {aka Rab22B}, TGFB1 (transforming growth factor beta 1) [NCBI Gene 7040] {aka CAEND1, CED, DPD1, IBDIMDE, LAP, TGF-beta1}, TAGLN (transgelin) [NCBI Gene 6876] {aka SM22, SM22-alpha, SMCC, TAGLN1, TGLN, WS3-10}, COL3A1 (collagen type III alpha 1 chain) [NCBI Gene 1281] {aka EDS4A, EDSVASC, PMGEDSV}, CLDN1 (claudin 1) [NCBI Gene 9076] {aka CLD1, ILVASC, SEMP1}, MEG3 (maternally expressed 3) [NCBI Gene 55384] {aka FP504, GTL2, LINC00023, Lnc-DLK1-35, NCRNA00023, PRO0518}, SKAP1 (src kinase associated phosphoprotein 1) [NCBI Gene 8631] {aka HEL-S-81p, SCAP1, SKAP55}, DOCK4 (dedicator of cytokinesis 4) [NCBI Gene 9732], SEMA6D (semaphorin 6D) [NCBI Gene 80031], ATP5MJ (ATP synthase membrane subunit j) [NCBI Gene 9556] {aka 6.8PL, ATP5MPL, C14orf2, MLQ, MP68, PLPM}, THRB (thyroid hormone receptor beta) [NCBI Gene 7068] {aka C-ERBA-2, C-ERBA-BETA, ERBA2, GRTH, NR1A2, PRTH}, RPL3 (ribosomal protein L3) [NCBI Gene 6122] {aka ASC-1, L3, TARBP-B, uL3}, GRIK2 (glutamate ionotropic receptor kainate type subunit 2) [NCBI Gene 2898] {aka EAA4, GLR6, GLUK6, GLUR6, GluK2, MRT6}, RPL35 (ribosomal protein L35) [NCBI Gene 11224] {aka DBA19, L35, uL29}, PDE10A (phosphodiesterase 10A) [NCBI Gene 10846] {aka ADSD2, HSPDE10A, IOLOD, PDE10A19}, NDUFB7 (NADH:ubiquinone oxidoreductase subunit B7) [NCBI Gene 4713] {aka B18, CI-B18, MC1DN39}, RPLP2 (ribosomal protein lateral stalk subunit P2) [NCBI Gene 6181] {aka D11S2243E, LP2, P2, RPP2}, COX3 (cytochrome c oxidase subunit III) [NCBI Gene 4514] {aka COIII, MTCO3}, RAB11FIP1 (RAB11 family interacting protein 1) [NCBI Gene 80223] {aka NOEL1A, RCP, rab11-FIP1}, ND1 (NADH dehydrogenase subunit 1) [NCBI Gene 4535] {aka MTND1}, GLIS3 (GLIS family zinc finger 3) [NCBI Gene 169792] {aka NDH, ZNF515}, NDUFB11 (NADH:ubiquinone oxidoreductase subunit B11) [NCBI Gene 54539] {aka CI-ESSS, ESSS, MC1DN30, NP17.3, Np15, P17.3}, TPO (thyroid peroxidase) [NCBI Gene 7173] {aka MSA, TDH2A, TPX}, DUOX1 (dual oxidase 1) [NCBI Gene 53905] {aka LNOX1, NOXEF1, THOX1}
- **Diseases:** Thyroid dysfunction (MESH:D013959), TB (MESH:D014390), Thyroid Cancer (MESH:D013964), Hypothyroidism (MESH:D007037), thyroid tumorigenesis (MESH:D063646), hyperthyroidism (MESH:D006980), thyroid (MESH:D013966), Hashimoto's thyroiditis (MESH:D050031), Graves' disease (MESH:D006111), fibrosis (MESH:D005355), inflammation (MESH:D007249), malignancy (MESH:D009369), ATC (MESH:D065646)
- **Chemicals:** ATP (MESH:D000255), Iodine (MESH:D007455), T3 (MESH:D014284), T4 (MESH:D013974), H2O2 (MESH:D006861), Debris (-), tyrosine (MESH:D014443), cyclic nucleotide (MESH:D009712), cyclic AMP (MESH:D000242), iodide (MESH:D007454)
- **Species:** Danio rerio (leopard danio, species) [taxon 7955], Homo sapiens (human, species) [taxon 9606], Mus musculus (house mouse, species) [taxon 10090]
- **Mutations:** brafv600e
- **Cell lines:** S2 — Drosophila melanogaster (Fruit fly), Spontaneously immortalized cell line (CVCL_Z232)

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12917210/full.md

## Figures

5 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12917210/full.md

## References

46 references — full list in the complete paper: https://tomesphere.com/paper/PMC12917210/full.md

---
Source: https://tomesphere.com/paper/PMC12917210