Single-cell transcriptome analysis reveals cell-cell communication and thyrocyte diversity in the zebrafish thyroid gland

The thyroid gland regulates growth and metabolism via production of thyroid hormone in follicles composed of thyrocytes. So far, thyrocytes have been assumed to be a homogenous population. To uncover genetic heterogeneity in the thyrocyte population, and molecularly characterize the non-thyrocyte cells surrounding the follicle, we developed a single-cell transcriptome atlas of the zebrafish thyroid gland. The 6249-cell atlas includes profiles of thyrocytes, blood vessels, lymphatic vessels, immune cells and fibroblasts. Further, the thyrocytes could be split into two sub-populations with unique transcriptional signature, including differential expression of the transcription factorpax2a. To validate thyrocyte heterogeneity, we generated a CRISPR/Cas9-basedpax2aknock-in line, which demonstrated specificpax2aexpression in the thyrocytes. However, a population ofpax2a-low mature thyrocytes interspersed within individual follicles could be distinguished, corroborating heterogeneity within the thyrocyte population. Our results identify and validate transcriptional differences within the nominally homogenous thyrocyte population.