Mutation screening of eight genes and comparison of the clinical data in a Chinese cohort with congenital hypothyroidism

Clinical efficacy of multigene panels in the management of congenital hypothyroidism with gland in situ

Congenital hypothyroidism (CHT) is a diverse condition with various genetic etiologies. This study aimed to investigate the utility of next-generation sequencing (NGS) analysis in guiding treatment decisions and predicting prognosis for CHT patients with gland in situ (GIS). A retrospective analysis was conducted on 33 CHT patients with GIS who underwent NGS analysis at a single institution between 2018 and 2023. Patients were classified as having permanent (PCH), transient congenital hypothyroidism, or ambiguous congenital hypothyroidism (ACH) CHT based on their response to levothyroxine discontinuation at 3 years of age. Among the 33 patients, genetic variants were identified in 26, with the most prevalent variants found in DUOX2 (26.92%), TSHR (30.77%), TG (19.35%), and DUOXA2 (19.23%). Patients with high initial thyroid-stimulating hormone levels (>50 mIU/L) and low free thyroxine levels (<0.89 ng/dL) at diagnosis tended to have compound heterozygous or homozygous variants in DUOX2, DUOXA2, and TG, and were more likely to develop PCH. In contrast, patients with heterozygous variants in these genes often exhibited ACH. TSHR variants were associated with diverse clinical manifestations, ranging from PCH to ACH, and were more common in patients with initial thyroid-stimulating hormone levels <50 mIU/L. The study highlights the potential utility of NGS analysis in predicting the clinical course and guiding treatment decisions for CHT patients with GIS. Genetic analysis may aid in determining the appropriate duration of levothyroxine therapy and monitoring strategies, particularly in cases where traditional clinical indicators are inconclusive.

GLIS3 expression in the thyroid gland in relation to TSH signaling and regulation of gene expression

Loss of GLI-Similar 3 (GLIS3) function in mice and humans causes congenital hypothyroidism (CH). In this study, we demonstrate that GLIS3 protein is first detectable at E15.5 of murine thyroid development, a time at which GLIS3 target genes, such as Slc5a5 (Nis), become expressed. This, together with observations showing that ubiquitous Glis3KO mice do not display major changes in prenatal thyroid gland morphology, indicated that CH in Glis3KO mice is due to dyshormonogenesis rather than thyroid dysgenesis. Analysis of GLIS3 in postnatal thyroid suggested a link between GLIS3 protein expression and blood TSH levels. This was supported by data showing that treatment with TSH, cAMP, or adenylyl cyclase activators or expression of constitutively active PKA enhanced GLIS3 protein stability and transcriptional activity, indicating that GLIS3 activity is regulated at least in part by TSH/TSHR-mediated activation of PKA. The TSH-dependent increase in GLIS3 transcriptional activity would be critical for the induction of GLIS3 target gene expression, including several thyroid hormone (TH) biosynthetic genes, in thyroid follicular cells of mice fed a low iodine diet (LID) when blood TSH levels are highly elevated. Like TH biosynthetic genes, the expression of cell cycle genes is suppressed in ubiquitous Glis3KO mice fed a LID; however, in thyroid-specific Glis3 knockout mice, the expression of cell cycle genes was not repressed, in contrast to TH biosynthetic genes. This indicated that the inhibition of cell cycle genes in ubiquitous Glis3KO mice is dependent on changes in gene expression in GLIS3 target tissues other than the thyroid.

Results of neonatal screening for congenital hypothyroidism and hyperphenylalaninemia in Zhejiang province from 1999 to 2022

OBJECTIVES

To analyze the results of neonatal screening for congenital hypothyroidism (CH) and hyperphenylalaninemia (HPA) in Zhejiang province from 1999 to 2022.

METHODS

A total of 11 922 318 newborns were screened from September 1999 and December 2022 in Zhejiang province. The blood thyroid stimulating hormone (TSH) levels were measured by a fluorescence method and blood phenylalanine (Phe) levels were measured by fluorescence method or tandem mass spectrometry. TSH≥9 μIU/mL was considered positive for CH, while Phe>120 μmol/L and/or Phe/Tyr ratio>2.0 were considered positive for HPA. The positive newborns in screening were recalled, and the gene variations were detected by high-throughput sequencing and MassARRAY tests.

RESULTS

The overall neonatal screening rate during 1999-2022 was 89.41% (11 922 318/13 333 929) and the screening rate was increased from 6.46% in 1999 to 100.0% in 2022. A total of 8924 cases of CH were diagnosed among screened newborns with an incidence rate of 1/1336. A total of 563 cases of HPA were diagnosed, including 508 cases of classic phenylketonuria (cPKU) and 55 cases of tetrahydrobiopterin deficiency (BH4D), with an incidence rate of 1/21 176. Ninety-seven out of 8924 cases of CH underwent genetic analysis. Gene mutations were detected in 9 CH related genes, the highest frequency mutations were found in DUOX2 gene (69.0%) with c.3329G>A (p.R1110Q) (18.2%) and c.1588A>T (p.K530X) (17.3%) as the hotspot mutations. There were 81 PAH gene variants detected in a total of 250 cases of cPKU, and c728G>A (p.R243Q) (24.4%), c.721C>T (p.R241C) (15.0%) were the hotspot mutations. Meanwhile 7 novel variants in PAH gene were detected: c.107C>A (p.S36*), c.137G>T (p.G46V), c.148A>G(p.K50E), c.285C>T (p.I95I), c.843-10delTTCC, exon4-7del and c.1066-2A>G. There were 12 PTS gene variants detected in 36 cases of BH4D, and c.259C>T (p.P87S) (31.9%) was the hotspot mutation.

CONCLUSIONS

The incident of CH has increased from 1999 to 2022 in Zhejiang province, and it is higher than that of national and global levels; while the incidence of HPA is similar to the national average. DUOX2 gene variation is the most common in CH patients; c.728G>A (p.R243Q) is the hotspot mutation in cPKU patients, while c.259C>T (p.P87S) is the hotspot mutation in BH4D patients.

Role of GLIS3 in thyroid development and in the regulation of gene expression in thyroid specific Glis3KO mice

Loss of GLI-Similar 3 (GLIS3) function in mice and humans causes congenital hypothyroidism (CH). In this study, we demonstrate that GLIS3 protein is first detectable at E15.5 of murine thyroid development, a time when GLIS3 target genes, such as Slc5a5 (Nis), become also expressed. We further show that Glis3KO mice do not display any major changes in prenatal thyroid gland morphology indicating that CH in Glis3KO mice is due to dyshormonogenesis rather than thyroid dysgenesis. Analysis of thyroid-specific Glis3 knockout (Glis3-Pax8Cre) mice fed either a normal or low-iodine diet (ND or LID) revealed that, in contrast to ubiquitous Glis3KO mice, thyroid follicular cell proliferation and the expression of cell cycle genes were not repressed suggesting that the inhibition of thyroid follicular cell proliferation in ubiquitous Glis3KO mice is related to loss of GLIS3 function in other cell types. However, the expression of several thyroid hormone biosynthesis-, extracellular matrix (ECM)-, and inflammation-related genes was still suppressed in Glis3-Pax8Cre mice particularly under conditions of high blood levels of thyroid stimulating hormone (TSH). We further demonstrate that treatment with TSH, protein kinase A (PKA) or adenylyl cyclase activators or expression of constitutively active PKA enhances GLIS3 protein and activity, suggesting that GLIS3 transcriptional activity is regulated in part by TSH/TSHR-mediated activation of the PKA pathway. This mechanism of regulation provides an explanation for the dramatic increase in GLIS3 protein expression and the subsequent induction of GLIS3 target genes, including several thyroid hormone biosynthetic genes, in thyroid follicular cells of mice fed a LID.

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.