GLIS3 regulates transcription of thyroid hormone biosynthetic genes in coordination with other thyroid transcription factors

BACKGROUND: Loss of the transcription factor GLI-Similar 3 (GLIS3) function causes congenital hypothyroidism (CH) in both humans and mice due to decreased expression of several thyroid hormone (TH) biosynthetic genes in thyroid follicular cells. Whether and to what extent, GLIS3 regulates thyroid gene transcription in coordination with other thyroid transcriptional factors (TFs), such as PAX8, NKX2.1 and FOXE1, is poorly understood. METHODS: PAX8, NKX2.1, and FOXE1 ChIP-Seq analysis with mouse thyroid glands and rat thyrocyte PCCl3 cells was performed and compared to that of GLIS3 to analyze the co-regulation of gene transcription in thyroid follicular cells by these TFs. RESULTS: Analysis of the PAX8, NKX2.1, and FOXE1 cistromes identified extensive overlaps between these TF binding loci and those of GLIS3 indicating that GLIS3 shares many of the same regulatory regions with PAX8, NKX2.1, and FOXE1, particularly in genes associated with TH biosynthesis, induced by thyroid stimulating hormone (TSH), and suppressed in Glis3KO thyroid glands, including Slc5a5 (Nis), Slc26a4, Cdh16, and Adm2. ChIP-QPCR analysis showed that loss of GLIS3 did not significantly affect PAX8 or NKX2.1 binding and did not cause major alterations in H3K4me3 and H3K27me3 epigenetic signals. CONCLUSIONS: Our study indicates that GLIS3 regulates transcription of TH biosynthetic and TSH-inducible genes in thyroid follicular cells in coordination with PAX8, NKX2.1, and FOXE1 by binding within the same regulatory hub. GLIS3 does not cause major changes in chromatin structure at these common regulatory regions. GLIS3 may induce transcriptional activation by enhancing the interaction of these regulatory regions with other enhancers and/or RNA Polymerase II (Pol II) complexes. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s13578-023-00979-8.