{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,4,8]],"date-time":"2026-04-08T04:29:45Z","timestamp":1775622585701,"version":"3.50.1"},"reference-count":46,"publisher":"Frontiers Media SA","license":[{"start":{"date-parts":[[2026,1,29]],"date-time":"2026-01-29T00:00:00Z","timestamp":1769644800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["frontiersin.org"],"crossmark-restriction":true},"short-container-title":["Front. Endocrinol."],"abstract":"<jats:p>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.<\/jats:p>","DOI":"10.3389\/fendo.2026.1748458","type":"journal-article","created":{"date-parts":[[2026,1,29]],"date-time":"2026-01-29T05:22:55Z","timestamp":1769664175000},"update-policy":"https:\/\/doi.org\/10.3389\/crossmark-policy","source":"Crossref","is-referenced-by-count":1,"title":["Single-cell analysis reveals the cellular and transcriptional heterogeneity of the normal pediatric thyroid"],"prefix":"10.3389","volume":"17","author":[{"given":"Erin R.","family":"Reichenberger","sequence":"first","affiliation":[{"name":"Department of Biomedical and Health Informatics, Children\u2019s Hospital of Philadelphia","place":["Philadelphia, PA, United States"]}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Nicholas E.","family":"Bambach","sequence":"additional","affiliation":[{"name":"Division of Endocrinology and Diabetes, Children\u2019s Hospital of Philadelphia","place":["Philadelphia, PA, United States"]}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Zachary","family":"Spangler","sequence":"additional","affiliation":[{"name":"Division of Endocrinology and Diabetes, Children\u2019s Hospital of Philadelphia","place":["Philadelphia, PA, United States"]}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Julio C.","family":"Ricarte-Filho","sequence":"additional","affiliation":[{"name":"Division of Endocrinology and Diabetes, Children\u2019s Hospital of Philadelphia","place":["Philadelphia, PA, United States"]}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Kyle","family":"Hinkle","sequence":"additional","affiliation":[{"name":"Division of Endocrinology and Diabetes, Children\u2019s Hospital of Philadelphia","place":["Philadelphia, PA, United States"]}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Amber","family":"Isaza","sequence":"additional","affiliation":[{"name":"Division of Endocrinology and Diabetes, Children\u2019s Hospital of Philadelphia","place":["Philadelphia, PA, United States"]}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Tricia","family":"Bhatti","sequence":"additional","affiliation":[{"name":"Department of Pathology and Laboratory Medicine, Children\u2019s Hospital of Philadelphia, University of Pennsylvania","place":["Philadelphia, PA, United States"]}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Andrew J.","family":"Bauer","sequence":"additional","affiliation":[{"name":"Division of Endocrinology and Diabetes, Children\u2019s Hospital of Philadelphia","place":["Philadelphia, PA, United States"]}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Aime T.","family":"Franco","sequence":"additional","affiliation":[{"name":"Division of Endocrinology and Diabetes, Children\u2019s Hospital of Philadelphia","place":["Philadelphia, PA, United States"]},{"name":"Abramson Cancer Center, University of Pennsylvania","place":["Philadelphia, PA, United States"]}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"1965","published-online":{"date-parts":[[2026,1,29]]},"reference":[{"key":"B1","doi-asserted-by":"publisher","DOI":"10.1038\/nrendo.2018.18","article-title":"Global epidemiology of hyperthyroidism and hypothyroidism","volume":"14","author":"Taylor","year":"2018","journal-title":"Nat Rev Endocrinol"},{"key":"B2","doi-asserted-by":"publisher","DOI":"10.1001\/jamapediatrics.2016.0486","article-title":"Thyroid disorders in children and adolescents: A review","volume":"170","author":"Hanley","year":"2016","journal-title":"JAMA Pediatr"},{"key":"B3","doi-asserted-by":"publisher","first-page":"e50612","DOI":"10.15252\/embr.202050612","article-title":"Single-cell transcriptome analysis reveals thyrocyte diversity in the zebrafish thyroid gland","volume":"21","author":"Gillotay","year":"2020","journal-title":"EMBO Rep"},{"key":"B4","doi-asserted-by":"publisher","DOI":"10.1210\/endocr\/bqad029","article-title":"Single cell analysis of human thyroid reveals the transcriptional signatures of aging","volume":"164","author":"Hong","year":"2023","journal-title":"Endocrinology"},{"key":"B5","doi-asserted-by":"publisher","first-page":"451","DOI":"10.4172\/2155-9899.1000451","article-title":"Immune suppression mediated by myeloid and lymphoid derived immune cells in the tumor microenvironment facilitates progression of thyroid cancers driven by hras(G12v) and pten loss","volume":"7","author":"Jolly","year":"2016","journal-title":"J Clin Cell Immunol"},{"key":"B6","doi-asserted-by":"publisher","DOI":"10.1158\/0008-5472.CAN-15-2351","article-title":"Fibroblast-mediated collagen remodeling within the tumor microenvironment facilitates progression of thyroid cancers driven by brafv600e and pten loss","volume":"76","author":"Jolly","year":"2016","journal-title":"Cancer Res"},{"key":"B7","doi-asserted-by":"publisher","DOI":"10.3389\/fendo.2022.1019072","article-title":"Integrated analysis of fibroblasts molecular features in papillary thyroid cancer combining single-cell and bulk rna sequencing technology","volume":"13","author":"Li","year":"2022","journal-title":"Front Endocrinol (Lausanne)"},{"key":"B8","doi-asserted-by":"publisher","DOI":"10.1172\/jci169653","article-title":"Anaplastic transformation in thyroid cancer revealed by single-cell transcriptomics","volume":"133","author":"Lu","year":"2023","journal-title":"J Clin Invest"},{"key":"B9","doi-asserted-by":"publisher","first-page":"6058","DOI":"10.1038\/s41467-021-26343-3","article-title":"Single-cell transcriptomic analysis of the tumor ecosystems underlying initiation and progression of papillary thyroid carcinoma","volume":"12","author":"Pu","year":"2021","journal-title":"Nat Commun"},{"key":"B10","article-title":"Percentage of epithelial cells, fibroblasts and endothelial cells in the dog thyroid","volume":"180","author":"Dow","year":"1986","journal-title":"Comptes Rendus Des Seances la Societe Biologie ses Filiales"},{"key":"B11","doi-asserted-by":"publisher","DOI":"10.1002\/dvdy.20653","article-title":"The developing mouse thyroid: embryonic vessel contacts and parenchymal growth pattern during specification, budding, migration, and lobulation","volume":"235","author":"Fagman","year":"2006","journal-title":"Dev Dynamics"},{"key":"B12","doi-asserted-by":"publisher","first-page":"299","DOI":"10.1002\/aja.1000480203","article-title":"Innervation of the thyroid gland. Ii. Origin and course of the thyroid nerves in the dog","volume":"48","author":"Nonidez","year":"1931","journal-title":"Am J Anat"},{"key":"B13","doi-asserted-by":"publisher","first-page":"51","DOI":"10.1016\/j.exphem.2018.09.004","article-title":"The human cell atlas bone marrow single-cell interactive web portal","volume":"68","author":"Hay","year":"2018","journal-title":"Exp Hematol"},{"key":"B14","doi-asserted-by":"publisher","DOI":"10.1038\/s41586-018-0414-6","article-title":"Rna velocity of single cells","volume":"560","author":"La Manno","year":"2018","journal-title":"Nature"},{"key":"B15","doi-asserted-by":"publisher","first-page":"495","DOI":"10.1038\/nbt.3192","article-title":"Spatial reconstruction of single-cell gene expression data","volume":"33","author":"Satija","year":"2015","journal-title":"Nat Biotechnol"},{"key":"B16","doi-asserted-by":"publisher","DOI":"10.3389\/fimmu.2018.02425","article-title":"A single-cell sequencing guide for immunologists","volume":"9","author":"See","year":"2018","journal-title":"Front Immunol"},{"key":"B17","doi-asserted-by":"publisher","first-page":"1888","DOI":"10.1016\/j.cell.2019.05.031","article-title":"Comprehensive integration of single-cell data","volume":"177","author":"Stuart","year":"2019","journal-title":"Cell"},{"key":"B18","doi-asserted-by":"publisher","DOI":"10.1038\/s41590-018-0276-y","article-title":"Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage","volume":"20","author":"Aran","year":"2019","journal-title":"Nat Immunol"},{"key":"B19","doi-asserted-by":"publisher","DOI":"10.1038\/nbt.4096","article-title":"Integrating single-cell transcriptomic data across different conditions, technologies, and species","volume":"36","author":"Butler","year":"2018","journal-title":"Nat Biotechnol"},{"key":"B20","doi-asserted-by":"publisher","first-page":"3573","DOI":"10.1016\/j.cell.2021.04.048","article-title":"Integrated analysis of multimodal single-cell data","volume":"184","author":"Hao","year":"2021","journal-title":"Cell"},{"key":"B21","doi-asserted-by":"publisher","first-page":"1688","DOI":"10.1038\/s41467-022-29212-9","article-title":"Comprehensive generation, visualization, and reporting of quality control metrics for single-cell rna sequencing data","volume":"13","author":"Hong","year":"2022","journal-title":"Nat Commun"},{"key":"B22","doi-asserted-by":"publisher","first-page":"741","DOI":"10.12688\/f1000research.14966.1","article-title":"Isee: interactive summarizedexperiment explorer","volume":"7","author":"Rue-Albrecht","year":"2018","journal-title":"F1000Res"},{"key":"B23","doi-asserted-by":"publisher","first-page":"15","DOI":"10.1186\/s13059-017-1382-0","article-title":"Scanpy: large-scale single-cell gene expression data analysis","volume":"19","author":"Wolf","year":"2018","journal-title":"Genome Biol"},{"key":"B24","doi-asserted-by":"publisher","first-page":"e1006245","DOI":"10.1371\/journal.pcbi.1006245","article-title":"Exploring the single-cell rna-seq analysis landscape with the scrna-tools database","volume":"14","author":"Zappia","year":"2018","journal-title":"PloS Comput Biol"},{"key":"B25","doi-asserted-by":"publisher","first-page":"e1007794","DOI":"10.1371\/journal.pcbi.1007794","article-title":"Scedar: A scalable python package for single-cell rna-seq exploratory data analysis","volume":"16","author":"Zhang","year":"2020","journal-title":"PloS Comput Biol"},{"key":"B26","doi-asserted-by":"publisher","DOI":"10.1038\/ncomms14049","article-title":"Massively parallel digital transcriptional profiling of single cells","volume":"8","author":"Zheng","year":"2017","journal-title":"Nat Commun"},{"key":"B27","doi-asserted-by":"publisher","DOI":"10.1101\/2025.05.14.654073","article-title":"Swans: A highly configurable analysis pipeline for single-cell and single-nucleus rna-sequencing data","author":"Beigel","year":"2025","journal-title":"bioRxiv"},{"key":"B28","doi-asserted-by":"publisher","DOI":"10.1093\/gigascience\/giaa151","article-title":"Soupx removes ambient rna contamination from droplet-based single-cell rna sequencing data","volume":"9","author":"Young","year":"2020","journal-title":"GigaScience"},{"key":"B29","doi-asserted-by":"publisher","first-page":"329","DOI":"10.1016\/j.cels.2019.03.003","article-title":"Doubletfinder: doublet detection in single-cell rna sequencing data using artificial nearest neighbors","volume":"8","author":"McGinnis","year":"2019","journal-title":"Cell Syst"},{"key":"B30","doi-asserted-by":"publisher","DOI":"10.1093\/gigascience\/giy083","article-title":"Clustering trees: A visualization for evaluating clusterings at multiple resolutions","volume":"7","author":"Zappia","year":"2018","journal-title":"GigaScience"},{"key":"B31","unstructured":"R: A language and environment for statistical computing\n          \n          \n          \n          2013"},{"key":"B32","unstructured":"Fast gene set enrichment analysis\n          \n          \n            \n              Korotkevich G\n              SV\n            \n            \n              Sergushichev\n              A\n            \n          \n          \n          2024"},{"key":"B33","doi-asserted-by":"publisher","DOI":"10.1073\/pnas.0506580102","article-title":"Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles","volume":"102","author":"Subramanian","year":"2005","journal-title":"Proc Natl Acad Sci"},{"key":"B34","doi-asserted-by":"publisher","DOI":"10.1016\/j.cels.2015.12.004","article-title":"The molecular signatures database hallmark gene set collection","volume":"1","author":"Liberzon","year":"2015","journal-title":"Cell Syst"},{"key":"B35","doi-asserted-by":"publisher","DOI":"10.1093\/bioinformatics\/btr260","article-title":"Molecular signatures database (Msigdb) 3.0","volume":"27","author":"Liberzon","year":"2011","journal-title":"Bioinformatics"},{"key":"B36","doi-asserted-by":"publisher","first-page":"581","DOI":"10.1210\/er.2003-0027","article-title":"Vascular endothelial growth factor: basic science and clinical progress","volume":"25","author":"Ferrara","year":"2004","journal-title":"Endocr Rev"},{"key":"B37","doi-asserted-by":"publisher","DOI":"10.3389\/fimmu.2019.00518","article-title":"The interplay between lymphatic vessels and chemokines","volume":"10","author":"Farnsworth","year":"2019","journal-title":"Front Immunol"},{"key":"B38","doi-asserted-by":"publisher","DOI":"10.1161\/CIRCULATIONAHA.109.910703","article-title":"Targeting lymphatic vessel activation and ccl21 production by vascular endothelial growth factor receptor-3 inhibition has novel immunomodulatory and antiarteriosclerotic effects in cardiac allografts","volume":"121","author":"Nyk\u00e4nen","year":"2010","journal-title":"Circulation"},{"key":"B39","doi-asserted-by":"publisher","first-page":"217","DOI":"10.1186\/s12864-019-5551-2","article-title":"Tadkb: family classification and a knowledge base of topologically associating domains","volume":"20","author":"Liu","year":"2019","journal-title":"BMC Genomics"},{"key":"B40","doi-asserted-by":"publisher","DOI":"10.3389\/fmolb.2021.793493","article-title":"Mir100hg regulates cald1 gene expression by targeting mir-142-5p to affect the progression of bladder cancer cells in vitro, as revealed by transcriptome sequencing","volume":"8","author":"Zhang","year":"2021","journal-title":"Front Mol Biosci"},{"key":"B41","doi-asserted-by":"publisher","DOI":"10.1016\/j.cell.2014.09.050","article-title":"Integrated genomic characterization of papillary thyroid carcinoma","volume":"159","year":"2014","journal-title":"Cell"},{"key":"B42","doi-asserted-by":"publisher","first-page":"32","DOI":"10.1186\/s13578-023-00979-8","article-title":"Glis3 regulates transcription of thyroid hormone biosynthetic genes in coordination with other thyroid transcription factors","volume":"13","author":"Kang","year":"2023","journal-title":"Cell Biosci"},{"key":"B43","unstructured":"2024"},{"key":"B44","doi-asserted-by":"publisher","DOI":"10.1016\/j.semcancer.2013.02.001","article-title":"Immunomonitoring and prognostic relevance of neutrophils in clinical trials","volume":"23","author":"Donskov","year":"2013","journal-title":"Semin Cancer Biol"},{"key":"B45","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1007\/978-3-319-22542-5_1","article-title":"Thyroid C-cell biology and oncogenic transformation","volume":"204","author":"Cote","year":"2015","journal-title":"Recent Results Cancer Res"},{"key":"B46","doi-asserted-by":"publisher","first-page":"216","DOI":"10.1038\/s41420-023-01514-x","article-title":"The lncrna rmst is drastically downregulated in anaplastic thyroid carcinomas where exerts a tumor suppressor activity impairing epithelial-mesenchymal transition and stemness","volume":"9","author":"De Martino","year":"2023","journal-title":"Cell Death Discov"}],"container-title":["Frontiers in Endocrinology"],"original-title":[],"link":[{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fendo.2026.1748458\/full","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,2,24]],"date-time":"2026-02-24T09:16:23Z","timestamp":1771924583000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fendo.2026.1748458\/full"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,1,29]]},"references-count":46,"alternative-id":["10.3389\/fendo.2026.1748458"],"URL":"https:\/\/doi.org\/10.3389\/fendo.2026.1748458","relation":{},"ISSN":["1664-2392"],"issn-type":[{"value":"1664-2392","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,1,29]]},"article-number":"1748458"}}