Next Article in Journal
Inflammaging and Blood Pressure Profiles in Late Life: The Screening for CKD among Older People across Europe (SCOPE) Study
Previous Article in Journal
Thrombotic Events during Lenvatinib Treatment: A Single Institution Experience
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Factors Causing Thyroid Dyshormonogenesis as the Major Etiologies for Primary Congenital Hypothyroidism: Clinical and Genetic Characterization of 33 Patients

1
Department of Pediatrics, The First Affiliated Hospital, Jinan University, Guangzhou 510630, China
2
Department of Pediatrics, Huizhou No. 2 Women’s and Children’s Healthcare Hospital, Huizhou 516000, China
3
Neonatal Screening Center, Dongguan Maternal and Child Healthcare Hospital, Dongguan 523125, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
J. Clin. Med. 2022, 11(24), 7313; https://doi.org/10.3390/jcm11247313
Submission received: 15 November 2022 / Revised: 1 December 2022 / Accepted: 7 December 2022 / Published: 9 December 2022
(This article belongs to the Section Endocrinology & Metabolism)

Abstract

:
Background and aims: Although the significance of primary congenital hypothyroidism (CH) is supported by an increasing amount of evidence, the clinical and genetic characteristics of this condition are still poorly understood. This study aimed to explore the underlying genetic etiologies in a cohort of primary CH patients. Subjects and Methods: The clinical data of 33 patients with primary CH were collected and analyzed via a cross-sectional study. Genetic analysis was performed by high-throughput sequencing and Sanger verification, and the pathogenicity of the novel missense variants was predicted using a variety of comprehensive bioinformatic tools. Results: Among the 33 patients, 22 (22/33, 66.7%) harbored pathogenic variants in the causative genes of thyroid dysgenesis or dyshormonogenesis, with DUOX2 (15/33, 45.5%) topping the list, followed by TG, TPO, DUOXA2 and PAX8. Four novel genetic variants were detected, including a pathogenic frameshift and three likely pathogenic missense variants. Positive neonatal screening for TSH, neonatal jaundice and abnormal thyroid morphology were the main positive findings among all cases. Although 31 of the total 33 CH patients exhibited normal anthropometric and social performance, the other 2 had poor prognosis in this study. Conclusions: This study reported 33 new CH patients bearing four novel genetic variants, which enriched the variant spectrum of CH genes. In this cohort, genetic factors causing thyroid dyshormonogenesis were the main etiologies of CH development. Most patients exhibited a favorable prognosis; however, systematic management remains a challenge in achieving improved clinical outcomes for CH patients.

1. Introduction

Congenital hypothyroidism (CH) is estimated to occur in approximately 1 in 3000–4000 live births worldwide, and timely diagnosis is crucial, considering the essential role of thyroid hormones in central nervous system maturation [1,2,3]. For years, the most prevalent cause of CH was iodine deficiency, but today, the significance of defects in the thyroid gland itself, which is defined as primary CH, is supported by an increasing amount of evidence [4]. Due to the advent of next-generation sequencing (NGS), more than 24 genes have been identified to be linked with primary CH [5,6], which is divided into thyroid dysgenesis (TD) and thyroid dyshormonogenesis (TDH) [7]. The former, which exhibits abnormal thyroid gland organogenesis, is associated with the genes TSHR, TUBB1, PAX8, NKX2-1, NKX2-5, GLIS3 CDC8A and NTN1 [3,7,8,9], while the latter, which is characterized by defects in a structurally normal thyroid gland to produce normal quantities of thyroid hormone, is suspected to be caused by the genes TG, TPO, SLC26A4, SLC26A7, SLC5A5, DUOX1, DUOX2 or DUOXA2 [3,10,11]. The inheritance model of TDH is generally autosomal recessive; however, the inheritance mode of TD remains debated, with several mechanisms that need to be considered, such as monogenic, polygenic and epigenetic inheritance [12].
The typical presentations of CH include myxedematous face, macroglossia, umbilical hernia, hypotonia, large fontanels, a distended abdomen, decreased activity and increased sleep, constipation, feeding difficulty and prolonged neonatal jaundice [5,13]. However, these clinical manifestations are often subtle at birth, and even negative in some cases with elevated thyroid-stimulating hormone (TSH) levels but normal thyroid hormone levels, making diagnosis of etiology and genetic counseling greatly challenging [14]. In recent years, next-generation sequencing (NGS) has facilitated the identification of the etiology of genetic conditions, and it is now possible to analyze the genetic etiologies of CH patients simultaneously and systematically [7,15].
The incidence of CH varies by race and ethnicity. It has been reported at 1:1016 in Asian people, 1:1815 in Caucasian people and 1:1902 in Black people [16]. China is a vast country with a huge population, and a TSH-based neonatal screening for CH revealed an incidence of approximately 1 in 2000–2500 live births [10,17,18]. However, the clinical and genetic characterization of this condition remain poorly understood in our country. The present study aimed to explore the underlying genetic etiologies in a cohort of primary CH patients.

2. Materials and Methods

2.1. Subjects and Ethical Approval

A total of 33 patients from November 2006 to July 2021, including three pairs of siblings (patients 1 and 2, 3 and 4 and 11 and 12), with primary CH were recruited as the research subjects. We recruited patients via a TSH-based neonatal screening program taking place between 72 h and 7 days after their births, with blood samples collected from the heel, and the TSH levels were measured by time-resolved fluorescence assay. All diagnoses were confirmed by at least one recall test(s) of the serum TSH, free thyroxine (FT3) and FT4, which were analyzed by electrochemistry immunoassay. CH was diagnosed based on TSH elevation (TSH ≥ 10 mIU/L) with low FT4 levels (FT4, with the reference range 12–27 pmol/L), and thereafter treated with levothyroxine (L-T4) at an initial dosage of 10–15 μg/kg per day, which was adjusted according to the serum levels of TSH and FT4 [14,19,20,21,22]. Patients with known causes of transient CH, including maternal use of antithyroid drugs, premature birth gestation (<37 weeks), history of iodine deficiency or iodine excess during pregnancy and patients who refused to undergo genetic testing or who were lost to follow-up were excluded. All the research subjects in this study met the CH diagnostic criteria [19,20,21].
The current study was conducted in accordance with the Declaration of Helsinki, and approved by the Committee for Medical Ethics, The First Affiliated Hospital of Jinan University (KY-2019-048 approved on 28 November 2019 and KY-2021-001 on 18 January 2021) and written informed consent was obtained from the parents of the patients.

2.2. Clinical Evaluation

Clinical data including main complaints, history, physical and laboratory examination, treatment and outcomes were collected via a comprehensive review of the patients’ medical records on serial clinic follow-up. Notably, all CH patients underwent a detailed thyroid ultrasonography examination in the longitudinal and axial planes by thyroid imaging specialists who had at least 4 years of experience in pediatric ultrasonography. Moreover, developmental quotient (DQ) evaluation was performed during the period of treatment by board-certified physiatrists. The scales used for DQ evaluation included Griffiths Mental Development Scales (GMDSs) and Denver Developmental Screening Tests (DDSTs).

2.3. Targeted NGS Analysis

Genomic DNA was extracted from peripheral blood leukocytes using a DNA extraction kit (Simen) according to the manufacturer’s protocol to create customized NimbleGen SeqCap probes (Roche NimbleGen, Madison, WI, USA) containing target genes associated with thyroid diseases, such as DUOX2, DUOXA2, TG, TSHR, TPO and PAX8. All exons and the flanking nucleotide sequences of the targeted genes were captured using a liquid-phase capture kit (Mygenostics, Beijing, China), and NGS was performed on the Illumina Hiseq 2000 sequencer (Illumina, San Diego, CA, USA). The overall sequencing coverage of the target regions was >95%, and the average depth of sequencing was not less than 200×. Sequence variant nomenclature was in agreement with current guidelines of the Human Genome Variation Society (http://www.hgvs.org/rec.html, accessed on 26 July 2022).

2.4. Sanger Sequencing

Sanger sequencing was performed to validate the relevant variants detected on NGS. DNA fragments were amplified by polymerase chain reaction (PCR) in 100μL of reaction system with the following temperature conditions: pre-denaturation at 98 °C for 1 min; 35 cycles of denaturation at 98 °C for 25 s, annealing at 65 °C for 30 s and elongation at 72 °C for 30 s; and a final extension at 72 °C for 5 min. The variants were identified by comparison with the reference sequence on the National Center Biotechnology Information (NCBI) website (https://www.ncbi.nlm.nih.gov/,accessed on 26 July 2022).
The sequencing results were aligned with the known gene sequence, which is available at Ensembl Genome Browser (www.ensembl.org, accessed on 26 July 2022), using the DNAman software (version 5.2.2; Lynnon Biosoft Corporation, San Ramon, CA, USA) and analyzed using the Chromas software (version 2.6.6; Technelysium Pty, Ltd., South Brisbane, Australia). The allele frequency of the identified novel variant was investigated in five population databases, including the 1000 Genomes Project (https://www.internationalgenome.org/, accessed on 26 July 2022), the Genome Aggregation Database (gnomAD, http://www.gnomad-sg.org/, accessed on 26 July 2022), the Exome Sequencing Project (https://esp.gs.washington.edu/drupal/, accessed on 26 July 2022), the Exome Aggregation Consortium (https://exac.broadinstitute.org, accessed on 26 July 2022) and the Human Gene Mutation Database (https://www.hgmd.cf.ac.uk/ac/index.php, accessed on 26 July 2022). Variants were filtered for allele frequencies < 1%. All the variants were classified according to the American College of Medical Genetics and Genomics (ACMG) guidelines [23].

2.5. Bioinformatic Analysis

A total of five prediction programs were used to predict the pathogenicity of the novel missense variants. PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2/, accessed on 26 July 2022) analysis identifies a variant as ‘probably damaging’ if the probability is >0.85 and ‘possibly damaging’ if the probability is >0.15 [24]. SIFT (http://sift.jcvi.org/www/SIFT_chr_coords_submit.html, accessed on 26 July 2022) identifies the variant as being ‘deleterious’ if the SIFT score is <0.05 [25]. MutationTaster (http://mutationtaster.org/MutationTaster/index.html, accessed on 26 July 2022) prediction identifies probabilities ~1 as a ‘high security’ of the prediction [26]. PROVEN (http://provean.jcvi.org/index.php, accessed on 26 July 2022) analysis identifies a variant as ‘deleterious’ if the score is equal to or less than −2.5, and ‘neutral’ if the score is above −2.5 [27]. MutationAssessor (http://mutationassessor.org/, accessed on 26 July 2022) scores a mutation by global and subfamily specific conservation patterns as low, medium or high [28].
Comparative alignment of the amino acid sequences of DUOX2, TG and PAX8 orthologs, collected from 151, 114 and 170 species, respectively, from the Ensembl Genome Browse (https://asia.ensembl.org/index.html, accessed on 26 July 2022), was performed. The 151 DUOX2 orthologs were classified into five taxonomy subgroups: primates (24 species), rodents and lagomorphs (33 species), other mammals (47 species), other vertebrates (44 species) and other species (3 species). The 114 TG orthologs were also classified into five taxonomy subgroups: primates (22 species), rodents and lagomorphs (30 species), other mammals (27 species), other vertebrates (33 species) and other species (2 species). Meanwhile, the 170 PAX8 orthologs were classified into four taxonomy subgroups: primates (21 species), rodents and lagomorphs (37 species), other mammals (45 species), other vertebrates (67 species) and other species (0 species).
Moreover, the three-dimensional (3D) structural models of DUOX2, TG and PAX8 protein were built using the online software SWISS-MODEL automated protein modeling server (https://swissmodel.expasy.org/, accessed on 26 July 2022). The tertiary structures of the wild-type and mutant protein were compared using SWISSPdb Viewer 4.1.0 (http://www.expasy.org/spdbv, accessed on 26 July 2022) to evaluate the effect of the variant on the structure of the molecule.

3. Results

3.1. Results of Genetic Testing

Among the 33 CH patients studied, 22 were found to have genetic causes, with a total rate of 66.7% (22/33). The genes causing TD (PAX8) and TDH (DUOX2, DUOXA2, TG and TPO) accounted for 3.0% (1/33) and 63.6% (21/33) of total cases, respectively (Table 1, Figure 1). DUOX2 was the most common CH-causing gene (15/33, 45.5%), followed by TG (2/33, 6.1%), TPO (2/33, 6.1%), DUOXA2 (2/33, 6.1%) and PAX8 (1/33, 3.0%). The genetic etiologies remained unclear in the other 11 (11/33, 33.3%) CH patients (Supplementary Table S1).
Of note, c.3904T > C and c.2018dupA in DUOX2 as well as c.7404G > C in TG were detected in trans with another reported pathogenic mutation, besides a heterozygous variant, c.307T > A, in the PAX8 gene, in the relevant genes (Figure 2). To the best of our knowledge, these four variants, including three missenses and one duplication, have not been previously reported in any official references in the database Pubmed; however, the Exome Aggregation Consortium database has recorded the variant c.7404G > C with an allele frequency of 0.004‰, and gnomAD cited the variants c.3904T > C and c.2018dupA with an allele frequency of 0.00398‰ and 0.00398‰. The four novel variants are not currently included in the 1000 Genomes Project, the Exome Sequencing Project or the Human Gene Mutation Database.

3.2. Bioinformatic Findings

In ployphen-2 analysis, the missense variants c.3904T > C (p.Trp1302Arg) and c.7404G > C (p.Lys2468Asn) were predicted to be ‘probably damaging’, while the c.307T > A (p.Phe103Ile) variant was predicated to be ‘benign’. The three variants were all predicted to be ‘disease-causing’ on Mutation Taster and ‘deleterious’ on PROVEAN and SIFT. Moreover, MutationAssessor yielded a result of ‘high’ for the first variant, ‘medium’ for the second and ‘neutral’ for the last. The results are shown in Supplementary Table S2.
To predict the pathogenicity of the three novel missense variants, the affected amino acids in a total of 151, 114 and 170 peptides homologous to human DUOX2, TG and PAX8, respectively, were comparatively aligned. The amino acids p.Trp1302 in DUOX2 and p.Phe103 in PAX8 were consistent in all homologous peptides (Supplementary Tables S3 and S4). The amino acid p.Lys2468 in TG was relatively conserved in 22/22 of the primates including humans (Supplementary Table S5), 30/30 of rodents and lagomorphs, 26/27 of the other mammals, 25/33 of the other vertebrates and 2/2 of the other species.
Moreover, in SWISS-MODEL analysis (Supplementary Figure S1), the three novel missense variants exhibited alterations in the length and/or number of the hydrogen bonds between specific amino acid residues in the DUOX2, TG and PAX8 proteins, thus distorting the molecular structure of the relevant protein.

3.3. Clinical Presentations

Among the 33 CH patients, 16 (16/33, 48.5%) presented with neonatal jaundice, 2 (2/33, 6.1%) with increased thyroid gland sizes and 4 (4/33, 12.1%) with thyroid aplasia, while another 2 (2/33, 6.1%) presented with thyroid ectopy in ultrasonography analysis (Table 2).
Oral L-T4 was given to 29 (29/33, 87.9%) CH patients as early as the neonatal period. However, in the other four (4/33, 12.1%) cases (patients 11, 21, 23 and 25), L-T4 supplementation was initiated at ages 36, 36, 46 and 90 days, respectively, due to incompliance or inaccessibility to a timely recall test, treatment and regular clinic follow-up (Table 3 and Supplementary Figure S2).
As a result of L-T4 uptake, the thyroid function indices in most cases (31/33, 93.9%) generally improved within eight weeks (Table 3 and Supplementary Figure S2). Although elevated TSH levels were observed occasionally in 24 patients (24/33, 72.7%), they all ultimately improved or recovered in response to subsequent adjustment of the L-T4 dosage (Supplementary Figure S2). Then, the L-T4 dosages were gradually reduced, and clinic follow-up was performed in 29 CH patients. Of note, in the other four affected individuals (patients 12, 13, 29 and 31), L-T4 treatment was eventually stopped at ages 4.5, 3.5, 0.1 and 4.3 years and thereafter their TSH levels were observed to be normal.
Regarding the clinical outcomes (n = 33) at the latest follow-up, at the age of 47 (17.5, 71.5) months, 28 patients exhibited normal anthropometric and social performance with a normal developmental quotient (DQ) (Supplementary Table S6). Although DQ results are not available for patients 3, 4 and 26, no abnormal anthropometric or social performance was observed by their parents or school teachers. However, patient 21 received a low DQ on evaluation with retardation of personal–social skills, fine motor skills and language development, while patient 25 exhibited microcephaly along with a significantly poor learning ability and school performance due to incompliance/inaccessibility to the recall test, treatment and regular follow-up (Table 2 and Supplementary Table S6).

4. Discussion

This study reported 33 new primary CH patients and identified four novel genetic variants. The molecular findings suggest that genetic factors constituted the main (61.1%, 22/33) CH etiologies in this cohort. According to the ACMG standards [23], the novel DUOX2 variant c.2018dupA (p.Val674Glyfs*39) was classified as pathogenic; since frameshift disrupted gene function (PVS1), the novel variant appeared in trans with known pathogenic variants (PM3), and the clinical phenotypes of the patient were specific for CH (PP4). Additionally, the novel DUOX2 variant c.3904T > C (p.Trp1302Arg), TG variant c.7404G > C (p.Lys2468Asn) and PAX8 variant c.307T > A (p.Phe103Ile) were all classified as likely pathogenic. In detail, they were absent or found at extremely low frequencies in the 1000 Genomes Project, Exome Sequencing Project, gnomAD and Exome Aggregation Consortium (PM2). The novel variants appeared in trans with known pathogenic variants (PM3), and a variety of bioinformatic tools suggested that the three novel missense variants may be detrimental to the functionality of the protein (PP3). Finally, the clinical phenotype of the four patients were quite specific for CH (PP4). It was well known that, when combined with other evidence of the disease in question, a variant classified as likely pathogenic has sufficient evidence for a health-care provider to use the molecular testing information in clinical decision making [23]. In combination with the clinical and laboratory findings, the four novel genetic variants enriched the variant spectrum of CH-associated genes and provided important molecular evidence for the diagnosis of etiology and genetic counseling for CH patients and their families.
The DUOX2- and TG-encoded proteins play important roles in the process of thyroid hormonogenesis. DUOX2 encodes a dual oxidase to generate hydrogen peroxide, which is needed by thyroid peroxidase for the incorporation of iodine into thyroglobulin [29]. The amino acid residue p.Trp1302 is located in the flavine adenosine dinucleotide (FAD)-binding domain of the DUOX2 protein (https://www.uniprot.org/, accessed on 26 July 2022), and the novel DUOX2 variant c.3904T > C (Trp1302Arg) distorted the domain structure and hence might impair its FAD-binding ability. Another novel DUOX2 variant c.2018dupA gives rise to a truncated protein (p.Val674Glyfs*39) lacking the C-terminal FAD-binding domain and two EF-hand motifs (https://www.uniprot.org/, accessed on 26 July 2022), which inevitably affect the dual oxidase activity of DUOX2 needed to generate hydrogen peroxide. Additionally, TG-encoded thyroglobulin functions as the substrate for the synthesis of thyroid hormones [5]. The novel TG variant c.7404G > C (p.Lys2468Asn) affects the Cholinesterase-like (ChEL) domain of the TG protein (https://www.uniprot.org/, accessed on 26 July 2022). Since a variety of ChEL domain mutants impair upstream TG folding and eliminate TG secretion from the apical membrane of the follicular cells to iodination sites in the follicle [30,31,32], this variant might cause CH via a similar mechanism by changing the TG molecule structure.
The gene PAX8 encodes the PAX8 protein, which is involved in the development of thyroid follicular cells and plays an essential role in the regulation of thyroid development, functional differentiation and organogenesis in human embryos [33,34], acting as a transcription factor characterized by the presence of a DNA-binding paired domain at the N-terminus [35]. The PAX8 novel variant c.307T > A affects the amino acid residue p.Phe103, which corresponds to its DNA-binding domain (https://www.uniprot.org/, accessed on 26 July 2022) and causes structural distortion of the PAX8 molecule via changing hydrogen bonds (Supplementary Figure S1). It has been reported that PAX8 variants involved in this domain usually cause reduced TG and TPO promoter transcriptional activation ability in transient transfection studies, due to an impaired DNA-binding ability; thus, the differentiation of thyroid cells could be influenced [33,34]. Overall, it is not surprising that the individuals harboring these four novel genetic variants presented with positive neonatal screening for TSH, neonatal jaundice and/or abnormal thyroid morphology, as shown in Table 2.
In the 22 CH patients with positively detected genetic causes, only one patient was associated with the TD-related gene PAX8, while 21 were associated with the TDH-related genes DUOX2, DUOXA2, TG and TPO (Table 1), suggesting that TDH is more common than TD as the major mechanism for CH development. This finding is consistent with previous publications [6,9,36,37]. Among the TDH-associated genes positively detected in this CH cohort, as shown in Table 1 and Figure 1, DUOX2 topped the list, followed by DUOXA2, TG and TPO. Of note, DUOX2 variants were frequently detected in CH patients of Asian origin, including those from China, Japan and Korea [10,36,38,39]. The mutated DUOX2 allele frequencies in the primary CH patients from the Chinese cities of Changzhou, Shanghai, Guangxi, Qingdao and Guangzhou were as high as 37.1%, 32.8%, 29.0%, 21.2% and 31.3%, respectively [40,41,42,43,44]. On the other hand, TD has been reported as the most common cause of CH [5,16,45,46], and the involved genes include TSHR, PAX8, NKX2-1, NKX2-5, GLIS3 and NTN1 [3,7,8,9]. Considering the possible sampling bias in this study, data for a larger sampling size are still in need to evaluate the significance of TDH and TD in CH development.
Additionally, there were still 11 patients with typical CH manifestations but unclear underlying genetic etiologies. This finding might have arisen to due to the following reasons. First, these patients might carry hidden genetic variants such as large deletions/insertions, variants in the noncoding regions affecting transcription factors or copy number variations (CNVs) [47]. Second, there might be more relevant CH-causing genes which have not been discovered yet [5]. Third, other modes of inheritance including epigenetics and multigenic involvement could not be excluded [9]. Actually, patient 27 in this study had oligogenic variants in the genes DUOX2 and TG, and the possibility of an oligogenic mechanism, a specific multigenic manner, for his CH development could not be completely ruled out.
It is well-known that oral L-T4 is the main treatment for primary CH patients. However, this treatment does not need to continue for the whole life duration of patients, especially in those with transient but not permanent CH [6,14,41,48]. The penetrance was 100% in the 22 CH-affected children in this study, and even in such cases, L-T4 withdrawal was achieved in patients 12, 13, 29 and 31 (Table 3). This phenomenon might be explained by the fact that the demand for thyroid hormone is five to seven times higher in the neonatal period than in the adult stage, and gradually reduces after the infantile period [41]. As a possible explanation, functional compensation might exist in CH-affected subjects. For example, patient 12 in this study exhibited biallelic DUOX2 variants (Table 1), and DUOX1 can compensate for the DUOX2 mutation in thyroid cells to produce H2O2. The insufficient supply of H2O2 using DUOX1 alone in the neonatal and infantile periods might become sufficient as thyroid hormone requirements decrease with age [41,48].
Once diagnosed and treated in time, primary CH patients usually demonstrate benign prognosis. Unfortunately, due to incompliance or inaccessibility to a timely recall test, treatment and regular clinic follow-up, two patients in this study exhibited intelligence impairment, and similar unfavorable prognoses have been observed in previous studies [42,48,49,50]. These findings strongly indicate that, although general neonatal screening has been performed for years, the systematic management of CH patients by timely recall tests, diagnosis, treatment and regular clinic follow-up remains an issue today.
This study reported 33 new CH patients bearing four genetic variants, which enriched the variant spectrum of CH genes, and provided genetic evidence supporting the diagnosis of etiology and genetic counseling for the affected individuals. In this cohort, genetic factors causing thyroid dyshormonogenesis were the main etiologies of CH development. Most patients exhibited favorable clinical outcomes, and in some transient cases, the L-T4 treatment was temporary. However, the systematic management of this condition remains a significant challenge.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jcm11247313/s1, Figure S1: DUOX2, TG and PAX8 molecular alterations caused by the novel missense variants c.3904T > C(p.Trp1302Arg), c.7404G > C (p.Lys2468Asn) and c.307T > A(p.Phe103Ile), respectively. The figure illustrates the ball-and-stick models of the human DUOX2 (A and B), TG (C and D) and PAX8 protein (E and F). The solid lines in white, red and blue represented carbon, oxygen and nitrogen atoms, respectively, while the dashed lines in green, hydrogen bonds. The figures A, C and E illustrated the wild-type models of DUOX2, TG and PAX8, respectively, while B, D and F, the mutant models of DUOX2, TG and PAX8, respectively. The three-dimensional (3D) structural models were built up using the online software SWISS-MODEL automated protein modeling server (https://swissmodel.expasy.org/, accessed on 26 July 2022). The tertiary structures of the wild type and mutant protein were compared by SWISSPdb Viewer 4.1.0 (http://www.expasy.org/spdbv, accessed on 26 July 2022); Figure S2: Laboratory indices over time in the patients. The arrow (↓) indicated the ages when the treatment was started. y: year; m: month; FT3: free thyroxine 3; FT4: free thyroxine 4; TSH: thyroid stimulating hormone; Test values: FT3 and FT4 in pmol/L; TSH in mIU/L; The data of Patient 3 before 3 years of age was lost; Table S1: Genetic findings of the 11 CH patients with unclear etiology; Table S2: Pathogenicity prediction of the three novel missense variants; Table S3: Comparative alignment of the primate peptides homologous for human DUOX2; Table S4: Comparative alignment of the primate peptides homologous for human PAX8; Table S5: Comparative alignment of the primate peptides homologous for human TG; Table S6: DQ evaluation in 33 CH patients on clinic follow up.

Author Contributions

Conceptualization, R.L. and Y.-Z.S.; methodology, R.L., J.-L.T., X.-L.H. and Y.-Z.S.; validation, R.L., J.-L.T., X.-L.H. and Y.-Z.S.; investigation, J.-L.T., X.-L.H. and Y.-Z.S.; data curation, R.L., J.-L.T., X.-L.H. and Y.-Z.S.; writing—original draft preparation, R.L.; writing, review and editing, R.L. and Y.-Z.S.; supervision, Y.-Z.S.; project administration, Y.-Z.S.; funding acquisition, Y.-Z.S.; All authors have read and agreed to the published version of the manuscript.

Funding

This paper was supported financially by the Science and Technology Plan Project of Guangzhou city (No. 202201020088) as well as Clinical Frontier Technology Program of the First Affiliated Hospital of Jinan University, China (No. JNU1AF-CFTP-2022-a01228).

Institutional Review Board Statement

This study was approved by the Committee for Medical Ethics, The First Affiliated Hospital of Jinan University (KY-2019-048 approved on 28 November 2019 and KY-2021-001 on 18 January 2021).

Informed Consent Statement

Written informed consent was obtained from the parents of the patients.

Data Availability Statement

The datasets presented in this study are available upon request to corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Deladoey, J.; Ruel, J.; Giguere, Y.; Van Vliet, G. Is the incidence of congenital hypothyroidism really increasing? A 20-year retrospective population-based study in Quebec. J. Clin. Endocrinol. Metab. 2011, 96, 2422–2429. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Fisher, D.A. Second International Conference on Neonatal thyroid screening: Progress report. J. Pediatr. 1983, 1025, 653–654. [Google Scholar] [CrossRef] [PubMed]
  3. Peters, C.; van Trotsenburg AS, P.; Schoenmakers, N. Diagnosis of Endocrine Disease: Congenital hypothyroidism: Update and perspectives. Eur. J. Endocrinol. 2018, 179, R297–R317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Szinnai, G. Clinical genetics of congenital hypothyroidism. Endocr. Dev. 2014, 26, 60–78. [Google Scholar]
  5. Kostopoulou, E.; Miliordos, K.; Spiliotis, B. Genetics of primary congenital hypothyroidism-a review. Hormones 2021, 20, 225–236. [Google Scholar] [CrossRef]
  6. Yamaguchi, T.; Nakamura, A.; Nakayama, K.; Hishimura, N.; Morikawa, S.; Ishizu, K.; Tajima, T. Targeted Next-Generation Sequencing for Congenital Hypothyroidism with Positive Neonatal TSH Screening. J. Clin. Endocrinol. Metab. 2020, 105, dgaa308. [Google Scholar] [CrossRef]
  7. Persani, L.; Rurale, G.; de Filippis, T.; Galazzi, E.; Muzza, M.; Fugazzola, L. Genetics and management of congenital hypothyroidism. Best Pract. Res. Clin. Endocrinol. Metab. 2018, 32, 387–396. [Google Scholar] [CrossRef]
  8. Cherella, C.E.; Wassner, A.J. Congenital hypothyroidism: Insights into pathogenesis and treatment. Int. J. Pediatr. Endocrinol. 2017, 2017, 11. [Google Scholar] [CrossRef]
  9. Sun, F.; Zhang, J.; Yang, C.; Gao, G.; Zhu, W.; Han, B.; Zhang, L.; Wan, Y.; Ye, X.; Ma, Y.; et al. The Genetic Characteristics of Congenital Hypothyroidism in China by Comprehensive Screening of 21 Candidate Genes. Eur. J. Endocrinol. 2018, 178, 623–633. [Google Scholar] [CrossRef] [Green Version]
  10. Tan, M.; Huang, Y.; Jiang, X.; Li, P.; Tang, C.; Jia, X.; Chen, Q.; Chen, W.; Sheng, H.; Feng, Y.; et al. The Prevalence, Clinical, and Molecular Characteristics of Congenital Hypothyroidism Caused by DUOX2 Mutations: A Population-Based Cohort Study in Guangzhou. Horm. Metab. Res. 2016, 48, 581–588. [Google Scholar] [CrossRef] [Green Version]
  11. Niu, D.M.; Hwang, B.; Chu, Y.K.; Liao, C.J.; Wang, P.L.; Lin, C.Y. High prevalence of a novel mutation (2268 insT) of the thyroid peroxidase gene in Taiwanese patients with total iodide organification defect, and evidence for a founder effect. J. Clin. Endocrinol. Metab. 2002, 87, 4208–4212. [Google Scholar] [CrossRef] [PubMed]
  12. Stoupa, A.; Kariyawasam, D.; Carre, A.; Polak, M. Update of Thyroid Developmental Genes. Endocrinol. Metab. Clin. N. Am. 2016, 45, 243–254. [Google Scholar] [CrossRef] [PubMed]
  13. Cao, X.; Jiang, X.; Dou, Z.; Rakeman, M.A.; Zhang, M.; O’Donnell, K.; Ma, T.; Amette, K.; DeLong, N.; DeLong, G.R. Timing of vulnerability of the brain to iodine deficiency in endemic cretinism. N. Engl. J. Med. 1994, 331, 1739–1744. [Google Scholar] [CrossRef]
  14. Fu, C.; Luo, S.; Zhang, S.; Wang, J.; Zheng, H.; Yang, Q.; Xie, B.; Hu, X.; Fan, X.; Luo, J.; et al. Next-generation sequencing analysis of DUOX2 in 192 Chinese subclinical congenital hypothyroidism (SCH) and CH patients. Clin. Chim. Acta 2016, 458, 30–34. [Google Scholar] [CrossRef]
  15. Lu, J.T.; Campeau, P.M.; Lee, B.H. Genotype–Phenotype Correlation—Promiscuity in the Era of Next-Generation Sequencing. N. Engl. J. Med. 2014, 371, 593–596. [Google Scholar] [CrossRef]
  16. Rastogi, M.V.; LaFranchi, S.H. Congenital hypothyroidism. Orphanet J. Rare Dis. 2010, 5, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Fu, C.; Luo, S.; Li, Y.; Li, Q.; Hu, X.; Li, M.; Zhang, Y.; Su, J.; Hu, X.; Chen, Y.; et al. The incidence of congenital hypothyroidism (CH) in Guangxi, China and the predictors of permanent and transient CH. Endocr. Connect. 2017, 6, 926–934. [Google Scholar] [CrossRef] [Green Version]
  18. Zhao, D.H.; Shen, Y.; Gong, J.M.; Meng, Y.; Su, L.; Zhang, X. Newborn screening for congenital hypothyroidism in Henan province, China. Clin. Chim. Acta 2016, 452, 58–60. [Google Scholar] [CrossRef]
  19. Wassner, A.J. Congenital hypothyroidism. Clin. Perinatol. 2018, 45, 1–18. [Google Scholar] [CrossRef]
  20. Wassner, A.J. Pediatric Hypothyroidism: Diagnosis and Treatment. Paediatr. Drugs 2017, 19, 291–301. [Google Scholar] [CrossRef]
  21. van Trotsenburg, P.; Stoupa, A.; Leger, J.; Rohrer, T.; Peters, C.; Fugazzola, L.; Cassio, A.; Heinrichs, C.; Beauloye, V.; Pohlenz, J.; et al. Congenital Hypothyroidism: A 2020–2021 Consensus Guidelines Update-An ENDO-European Reference Network Initiative Endorsed by the European Society for Pediatric Endocrinology and the European Society for Endocrinology. Thyroid 2021, 31, 387–419. [Google Scholar] [CrossRef] [PubMed]
  22. Group of Endocrinology and Inborn Metabolic Diseases. Consensus of diagnosis and treatment of congenital hypothyroidism. Chin. J. Pediatr. 2011, 49, 421–424. [Google Scholar]
  23. Richards, S.; Aziz, N.; Bale, S.; Bick, D.; Das, S.; Gastier-Foster, J.; Grody, W.W.; Hegde, M.; Lyon, E.; Spector, E.; et al. Standards and guidelines for the interpretation of sequence variants: A joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet. Med. 2015, 17, 405–424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Adzhubei, I.A.; Schmidt, S.; Peshkin, L.; Ramensky, V.E.; Gerasimova, A.; Bork, P.; Kondrashov, A.S.; Sunyaev, S.R. A method and server for predicting damaging missense mutations. Nat. Methods 2010, 7, 248–249. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kumar, P.; Henikoff, S.; Ng, P.C. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat. Protoc. 2009, 4, 1073–1081. [Google Scholar] [CrossRef] [PubMed]
  26. Schwarz, J.M.; Cooper, D.N.; Schuelke, M.; Seelow, D. MutationTaster2: Mutation prediction for the deep-sequencing age. Nat. Methods 2014, 11, 361–362. [Google Scholar] [CrossRef]
  27. Choi, Y.; Chan, A.P. PROVEAN web server: A tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics 2015, 31, 2745–2747. [Google Scholar] [CrossRef] [Green Version]
  28. Bacalhau, M.; Pratas, J.; Simoes, M.; Mendes, C.; Ribeiro, C.; Santos, M.J.; Diogo, L.; Macario, M.C.; Grazina, M. In silico analysis for predicting pathogenicity of five unclassified mitochondrial DNA mutations associated with mitochondrial cytopathies’ phenotypes. Eur. J. Med. Genet. 2017, 60, 172–177. [Google Scholar] [CrossRef]
  29. Wang, F.; Lu, K.; Yang, Z.; Zhang, S.; Lu, W.; Zhang, L.; Liu, S.; Yan, S. Genotypes and phenotypes of congenital goitre and hypothyroidism caused by mutations in dual oxidase 2 genes. Clin. Endocrinol. 2014, 81, 452–457. [Google Scholar] [CrossRef]
  30. Bruellman, R.J.; Watanabe, Y.; Ebrhim, R.S.; Creech, M.K.; Abdullah, M.A.; Dumitrescu, A.M.; Refetoff, S.; Weiss, R.E. Increased Prevalence of TG and TPO Mutations in Sudanese Children with Congenital Hypothyroidism. J. Clin. Endocrinol. Metab. 2020, 105, 1564–1572. [Google Scholar] [CrossRef]
  31. Lee, J.; Di Jeso, B.; Arvan, P. Maturation of thyroglobulin protein region I. J. Biol. Chem. 2011, 286, 33045–33052. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Vono-Toniolo, J.; Rivolta, C.M.; Targovnik, H.M.; Medeiros-Neto, G.; Kopp, P. Naturally occurring mutations in the thyroglobulin gene. Thyroid 2005, 15, 1021–1033. [Google Scholar] [CrossRef] [PubMed]
  33. Qian, F.; Li, G.Y.; Wu, X.J.; Jia, Q.; Lyu, G.T.; Wang, M.L.; Wang, J. Novel non-synonymous mutations of PAX8 in a cohort of Chinese with congenital hypothyroidism. Chin. Med. J. 2019, 132, 1322–1327. [Google Scholar] [CrossRef] [PubMed]
  34. Magliano, M.P.; di Lauro, R.D.; Zannini, M. Pax8 has a key role in thyroid cell differentiation. Proc. Natl. Acad. Sci. USA 2000, 97, 13144–13149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Iwahashi-Odano, M.; Fujisawa, Y.; Ogata, T.; Nakashima, S.; Muramatsu, M.; Narumi, S. Identification and functional characterization of a novel PAX8 mutation (p.His39Pro) causing familial thyroid hypoplasia. Clin. Pediatr. Endocrinol. 2020, 29, 173–178. [Google Scholar] [CrossRef] [PubMed]
  36. Sorapipatcharoen, K.; Tim-Aroon, T.; Mahachoklertwattana, P.; Chantratita, W.; Iemwimangsa, N.; Sensorn, I.; Panthan, B.; Jiaranai, P.; Noojarern, S.; Khlairit, P.; et al. DUOX2 variants are a frequent cause of congenital primary hypothyroidism in Thai patients. Endocr. Connect. 2020, 9, 1121–1134. [Google Scholar] [CrossRef]
  37. Lof, C.; Patyra, K.; Kuulasmaa, T.; Vangipurapu, J.; Undeutsch, H.; Jaeschke, H.; Pajunen, T.; Kero, A.; Krude, H.; Biebermann, H.; et al. Detection of Novel Gene Variants Associated with Congenital Hypothyroidism in a Finnish Patient Cohort. Thyroid 2016, 26, 1215–1224. [Google Scholar] [CrossRef]
  38. Narumi, S.; Muroya, K.; Asakura, Y.; Aachi, M.; Hasegawa, T. Molecular basis of thyroid dyshormonogenesis: Genetic screening in population-based Japanese patients. J. Clin. Endocrinol. Metab. 2011, 96, E1838–E1842. [Google Scholar] [CrossRef]
  39. Park, K.J.; Park, H.K.; Kim, Y.J.; Lee, K.R.; Park, J.H.; Park, J.H.; Park, H.D.; Lee, S.Y.; Kim, J.W. DUOX2 Mutations Are Frequently Associated with Congenital Hypothyroidism in the Korean Population. Ann. Lab. Med. 2016, 36, 145–153. [Google Scholar] [CrossRef] [Green Version]
  40. Long, W.; Lu, G.; Zhou, W.; Yang, Y.; Zhang, B.; Zhou, H.; Jiang, L.; Yu, B. Targeted next-generation sequencing of thirteen causative genes in Chinese patients with congenital hypothyroidism. Endocr. J. 2018, 65, 1019–1028. [Google Scholar] [CrossRef] [Green Version]
  41. Zheng, Z.; Yang, L.; Sun, C.; Wu, J.; Luo, F.; Zhou, W.; Lu, W. Genotype and phenotype correlation in a cohort of Chinese congenital hypothyroidism patients with DUOX2 mutations. Ann. Transl. Med. 2020, 8, 1649. [Google Scholar] [CrossRef] [PubMed]
  42. Fu, C.; Zhang, S.; Su, J.; Luo, S.; Zheng, H.; Wang, J.; Qin, H.; Chen, Y.; Shen, Y.; Hu, X.; et al. Mutation screening of DUOX2 in Chinese patients with congenital hypothyroidism. J. Endocrinol. Investig. 2015, 38, 1219–1224. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, F.; Zang, Y.; Li, M.; Liu, W.; Wang, Y.; Yu, X.; Li, H.; Wang, F.; Liu, S. DUOX2 and DUOXA2 Variants Confer Susceptibility to Thyroid Dysgenesis and Gland-in-situ With Congenital Hypothyroidism. Front. Endocrinol. 2020, 11, 237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Yu, B.; Long, W.; Yang, Y.; Wang, Y.; Jiang, L.; Cai, Z.; Wang, H. Newborn Screening and Molecular Profile of Congenital Hypothyroidism in a Chinese Population. Front. Genet. 2018, 9, 509. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Narumi, S.; Muroya, K.; Asakura, Y.; Adachi, M.; Hasegawa, T. Transcription factor mutations and congenital hypothyroidism: Systematic genetic screening of a population-based cohort of Japanese patients. J. Clin. Endocrinol. Metab. 2010, 95, 1981–1985. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Nilsson, M.; Fagman, H. Mechanisms of thyroid development and dysgenesis: An analysis based on developmental stages and concurrent embryonic anatomy. Curr. Top. Dev. Biol. 2013, 106, 123–170. [Google Scholar]
  47. De Marco, G.; Agretti, P.; Montanelli, L.; Di Cosmo, C.; Bagattini, B.; De Servi, M.; Ferrarini, E.; Dimida, A.; Freitas Ferreira, A.C.; Molinaro, A.; et al. Identification and functional analysis of novel dual oxidase 2 (DUOX2) mutations in children with congenital or subclinical hypothyroidism. J. Clin. Endocrinol. Metab. 2011, 96, E1335–E1339. [Google Scholar] [CrossRef] [Green Version]
  48. Maruo, Y.; Nagasaki, K.; Matsui, K.; Mimura, Y.; Mori, A.; Fukami, M.; Takeuchi, Y. Natural course of congenital hypothyroidism by dual oxidase 2 mutations from the neonatal period through puberty. Eur. J. Endocrinol. 2016, 174, 453–463. [Google Scholar] [CrossRef] [Green Version]
  49. Aminzadeh, M. Higher prevalence of permanent congenital hypothyroidism in the Southwest of Iran mostly caused by dyshormonogenesis: A five-year follow-up study. Arch. Endocrinol. Metab. 2018, 62, 602–608. [Google Scholar] [CrossRef] [Green Version]
  50. Cangul, H.; Aycan, Z.; Kendall, M.; Bas, V.N.; Saglam, Y.; Barrett, T.G.; Maher, E.R. A truncating DUOX2 mutation (R434X) causes severe congenital hypothyroidism. J. Pediatr. Endocrinol. Metab. 2014, 27, 323–327. [Google Scholar] [CrossRef]
Figure 1. Genetic analytic findings in the 33 CH patients.
Figure 1. Genetic analytic findings in the 33 CH patients.
Jcm 11 07313 g001
Figure 2. Sanger sequencing results of DUOX2, TG and PAX8 gene in Patients 1, 2, 3, 4, 20 and 22. (A) Patients 1 and 2 were compound heterozygous for the DUOX2 variants c.3904T > C (p.Trp1302Arg) and c.1588A > T (p.Lys530*), while their father was a carrier of c.3904T > C and their mother was a carrier of c.1588A > T. (B) Patients 3 and 4 were compound heterozygous for the TG variants c.4736_4746delTGATCTTCGACinsA (p.Val1579Glufs) and c.7404G > C (p.Lys2468Asn) from their father and mother, respectively. (C) Patient 20 was compound heterozygous for the DUOX2 variants c.3329G > A (p.Arg1110Gln) and c.2018dupA (p.Val674Glyfs*39) from father and mother, respectively. (D) Patient 22 was a CH patient with PAX8 variant c.307T > A (p.Phe103Ile) from her father, while her mother was normal.
Figure 2. Sanger sequencing results of DUOX2, TG and PAX8 gene in Patients 1, 2, 3, 4, 20 and 22. (A) Patients 1 and 2 were compound heterozygous for the DUOX2 variants c.3904T > C (p.Trp1302Arg) and c.1588A > T (p.Lys530*), while their father was a carrier of c.3904T > C and their mother was a carrier of c.1588A > T. (B) Patients 3 and 4 were compound heterozygous for the TG variants c.4736_4746delTGATCTTCGACinsA (p.Val1579Glufs) and c.7404G > C (p.Lys2468Asn) from their father and mother, respectively. (C) Patient 20 was compound heterozygous for the DUOX2 variants c.3329G > A (p.Arg1110Gln) and c.2018dupA (p.Val674Glyfs*39) from father and mother, respectively. (D) Patient 22 was a CH patient with PAX8 variant c.307T > A (p.Phe103Ile) from her father, while her mother was normal.
Jcm 11 07313 g002
Table 1. Genotyping results of the 22 CH patients with clear genetic etiologies.
Table 1. Genotyping results of the 22 CH patients with clear genetic etiologies.
PatientsMOI Causative
Genes
The Allele Frequency in gnomAD Genotypes
Index Patient Father Mother
1AR DUOX2 0.00398‰/0.0119‰ c.3904T > C (p.Trp1302Arg)/c.1588A > T (p.Lys530*) c.3904T > C (p.Trp1302Arg)/- c.1588A > T (p.Lys530*)/-
2AR DUOX2 0.00398‰/0.0119‰ c.3904T > C (p.Trp1302Arg)/c.1588A > T (p.Lys530*) c.3904T > C(p.Trp1302Arg)/- c.1588A > T (p.Lys530*)/-
3AR TG NA/0.00398‰ c.4736_4746delTGATCTTCGACinsA (p.Val1579Glufs)/c.7404G > C (p.Lys2468Asn) c.4736_4746delTGATCTTCGACinsA (p.Val1579Glufs)/- c.7404G > C (p.Lys2468Asn)/-
4AR TG NA/0.00398‰ c.4736_4746delTGATCTTCGACinsA (p.Val1579Glufs)/c.7404G > C (p.Lys2468Asn) c.4736_4746delTGATCTTCGACinsA (p.Val1579Glufs)/- c.7404G > C (p.Lys2468Asn)/-
5AR TPO 0.117‰ c.2268dupT (p.Glu757*)/c c.2268dupT (p.Glu757*) c.2268dupT (p.Glu757*)/- c.2268dupT (p.Glu757*)/-
6AR DUOX2 0.598‰; 0.343‰/
0.417‰
[c.4027C > T (p.Leu1343Phe); c.2048G > T (p.Arg683Leu)]/
c.2654G > T (p.Arg885Leu)
[c.4027C > T (p.Leu1343Phe); c.2048G > T (p.Arg683Leu)]/- c.2654G > T (p.Arg885Leu)/-
7AR DUOX2 0.417‰/0.113‰ c.2654G > T (p.Arg885Leu)/c.2654G > A (p.Arg885Gln) c.2654G > T (p.Arg885Leu)/- c.2654G > A (p.Arg885Gln)/-
8AR DUOX2 0.343‰; 0.598‰/
0.191‰; NA
[c.2048G > T (p.Arg683Leu); c.4027C > T (p.Leu1343Phe/
[c.3329G > A (p.Arg1110Gln); c.1094A > G (p.Gln365Arg)]
[c.2048G > T (p.Arg683Leu); c.4027C > T (p.Leu1343Phe)]/- [c.3329G > A (p.Arg1110Gln); c.1094A > G (p.Gln365Arg)]/-
9AR DUOX2 0.0119‰/
0.598‰; 0.343‰
c.1588A > T (p.Lys530*)/
[c.4027C > T (p.Leu1343Phe); c.2048G > T (p.Arg683Leu)]
c.1588A > T (p.Lys530*)/- [c.4027C > T (p.Leu1343Phe); c.2048G > T (p.Arg683Leu)]/-
10AR DUOX2 0.191‰/0.0191‰ c.3329G > A (p.Arg1110Gln)/c.1588A > T (p.Lys530*) c.3329G > A (p.Arg1110Gln)/- c.1588A > T (p.Lys530*)/-
11AR DUOX2 0.0119‰/
0.343‰; 0.598‰
c.1588A > T (p.Lys530*)/
[c.2048G > T (p.Arg683Leu); c.4027C > T (p.Leu1343Phe)]
c.1588A > T (p.Lys530*)/- [c.2048G > T (p.Arg683Leu); c.4027C > T (p.Leu1343Phe)]/-
12AR DUOX2 0.019‰/
0.343‰; 0.598‰
c.1588A > T (p.Lys530*)/
[c.2048G > T (p.Arg683Leu); c.4027C > T (p.Leu1343Phe)]
c.1588A > T (p.Lys530*)/- [c.2048G > T (p.Arg683Leu); c.4027C > T (p.Leu1343Phe)]/-
13AR DUOXA2 0.231‰/0.143‰ c.413dupA (p.Tyr138*)/c.738C > G (p.Tyr246*) c.413dupA (p.Tyr138*)/- c.738C > G (p.Tyr246*)/-
14AR DUOX2 0.417‰/0.113‰ c.2654G > T (p.Arg885Leu)/c.3693 + 1G > T c.2654G > T (p.Arg885Leu)/- c.3693 + 1G > T/-
15AR DUOXA2 0.231‰/0.00402‰ c.413dupA (p.Tyr138*)/c.515dupA (p. Tyr173Valfs*57) c.413dupA (p.Tyr138*)/- c.515dupA (p. Tyr173Valfs*57)/-
16AR DUOX2 0.0707‰/0.113‰ c.2104_2106delGGA (p.Gly702del)/c.3693 + 1G > Tc.2104_2106delGGA (p.Gly702del)/- c.3693 + 1G > T/-
17AR DUOX2 0.0119‰/0.0679‰ c.1588A > T(p.Lys530*)/c.1946C > A(p.Ala649Glu) c.1588A > T(p.Lys530*)/- c.1946C > A (p.Ala649Glu)/-
18AR DUOX2 0.191‰/0.0119‰ c.3329G > A(p.Arg1110Gln)/c.1588A > T(p.Lys530*) c.3329G > A(p.Arg1110Gln)/- c.1588A > T (p.Lys530*)/-
19AR DUOX2 0.0119‰/0.0278‰ c.1588A > T (p.Lys530*)/c.3478_3480del (p.Leu1160del) c.1588A > T (p.Lys530*)/- c.3478_3480del (p.Leu1160del)/-
20AR DUOX2 0.191‰/0.00398‰ c.3329G > A (p.Arg1110Gln)/c.2018dupA (p.Val674Glyfs*39) c.3329G > A (p.Arg1110Gln)/- c.2018dupA (p.Val674Glyfs*39)/-
21AR TPO 0.117‰/NA c.2268dupT (p.Glu757*)/c.1804C > T (p.Arg602Cys) c.2268dupT (p.Glu757*)/- c.1804C > T (p.Arg602Cys)/-
22AD PAX8 NA/- c.307T > A (p.Phe103Ile)/- c.307T > A (p.Phe103Ile)/- -/-
MOI: Mode of inheritance; AD: Autosomal dominant; AR: Autosomal recessive; TPO: Thyroid peroxidase; TG: Thyroglobulin; DUOX2: Dual oxidase 2; DUOXA2: dual oxidase maturation factor 2; PAX8: Paired Box Gene 8 Mutations; NA: Frequency of the variant was not available in gnomAD.
Table 2. Clinical information of the 33 CH patients.
Table 2. Clinical information of the 33 CH patients.
PatientsSexGA(Weeks)BW(g)Delivery ModeChief ComplaintFamily HistoryThyroidMorphology/Age at Latest Examination (y)OtherPresentationsClinical Outcomes/Assessment Ages (y)
1Female40+22700Vaginal deliveryAbnormal NSNNormal/3.0NoneNormal DQ/3.0
2Male39+32580Vaginal deliveryAbnormal NSHis sister had CHNormal/0.1NoneNormal DQ/0.5
3Female414500CaesareanAbnormal NSNIncreased/9.3MacrosomiaNormal school performance/10.0
4Male403500CaesareanAbnormal NSHis sister had CHNormal/2.7NoneNormal school performance/7.3
5Male40+42900Vaginal deliveryAbnormal NSMother hadHashimoto’s thyroiditisNormal/3.3NoneNormal DQ/6.7
6Male403800CaesareanAbnormal NSNNormal/3.8NoneNormal DQ/3.7
7Male382450Vaginal deliveryAbnormal NSNNormal/2.3Neonatal jaundiceNormal DQ/4.3
8Male40+13000Vaginal deliveryAbnormal NSNNormal/0.1Neonatal jaundiceNormal DQ/0.8
9Female42+33000CaesareanAbnormal NSNNormal/0.1NoneNormal DQ/0.9
10Female39+13300CaesareanAbnormal NSNIncreased/0.1Neonatal jaundiceNormal DQ/2.0
11Female382400Vaginal deliveryAbnormal NSNNormal/3.7Neonatal jaundiceNormal DQ/4.4
12Male403250Vaginal deliveryAbnormal NSHis sister had CHNormal/2.6NoneNormal DQ/6.3
13Female382970Vaginal deliveryAbnormal NSHer brother had CHNormal/9.2NoneNormal DQ/4.8
14Male37+52590CaesareanAbnormal NSNNormal/2.5One of the monozygotic twins, Full termNormal DQ/3.9
15Male38+53150Vaginal deliveryAbnormal NSNNormal/1.5Neonatal jaundiceNormal DQ/3.0
16Female393500CaesareanAbnormal NSMother hadHashimoto’s thyroiditisNormal/3.5Neonatal jaundiceNormal DQ/4.4
17Female39+13400Vaginal deliveryAbnormal NSNNormal/0.3Neonatal jaundiceNormal DQ/4.0
18Male382500CaesareanAbnormal NSNNormal/3.5NoneNormal DQ/3.5
19Female393400Vaginal deliveryAbnormal NSNNormal/0.1NoneNormal DQ/0.8
20Female39+13100Vaginal deliveryAbnormal NSHer sister had CHNormal/0.1Neonatal jaundiceNormal DQ/1.2
21Male413500Vaginal deliveryAbnormal NSNNormal/0.7Neonatal jaundiceLow DQ/4.5
22Female372740Vaginal deliveryAbnormal NSHer father had a history of hypothyroidismNormal/0.1NoneNormal DQ/2.2
23Female39+42780Vaginal deliveryAbnormal NSNNot doneNeonatal jaundiceAtrial septal defect, Pierre-Robin syndromeNormal DQ/0.9
24Female40+43000Vaginal deliveryAbnormal NSAunt had hyperthyroidismThyroid agenesis/0.1NoneNormal DQ/2.5
25Male39+62700Vaginal deliveryAbnormal NSMother had hyperthyroidism in late pregnancyNot doneNoneMicrocephaly and poor learning ability and school performance/6.3
26Female37+3900CaesareanAbnormal NSNNot doneOne of the monozygotic twins, Full term, Small for gestational ageNormal developmental milestones/0.8
27Female393100Vaginal deliveryAbnormal NSNThyroid agenesis/0.1Neonatal jaundiceNormal DQ/6.1
28Female40+62900Vaginal deliveryAbnormal NSNThyroid ectopy/4.7Neonatal jaundiceNormal DQ/3.0
29Male383450CaesareanAbnormal NSMother had hyperthyroidism in pregnancyNormal/0.1Neonatal jaundiceNormal DQ/0.3
30Female38+53330CaesareanAbnormal NSNThyroid ectopy/0.1NoneNormal DQ/2.8
31Male403370Vaginal deliveryAbnormal NSMother had hypothyroidismNormal/0.8Neonatal jaundiceNormal DQ/2.6
32Female413350CaesareanAbnormal NSMother had hyperthyroidismThyroid agenesis/2.4Neonatal jaundiceNormal DQ/1.6
33Female413870CaesareanAbnormal NSNThyroid agenesis/1.1Neonatal jaundiceNormal DQ/1.1
NS: newborn screening; N: no family history of thyroid diseases; GA: gestational age; BW: Birth weight; y: years; CH: congenital hypothyroidism; DQ: Developmental quotient.
Table 3. Thyroid function indices over time in the 33 CH patients.
Table 3. Thyroid function indices over time in the 33 CH patients.
PatientsAge When L-T4 Initiated
(days)
Diagnostic EvaluationInitial L-T4 Dosage
(µg/kg/day)
Days of Age at Normalization of Thyroid Function (Days)Maintenance L-T4 Dosage (µg/kg/day)/
the Latest Ages (Years)
FT3
(3.5–10 pmol/L)
FT4
(12–27 pmol/L)
T3
(1.3–3.8 nmol/L)
T4
(79–192 nmol/L)
TSH
(0.27–4.2 mIU/L)
1 212.15<5.15-->1008.3562.8/2.9
2 201.766.631.7630.80>10011.7373.7/0.5
3 21---->10011.1-1.8/10.0
4 223.655.131.7846.75>1009.43(y)3.5/7.3
5 27--0.31->15010.7573.2/6.7
6 215.2011.801.9087.8222.428.5901.0/3.7
7 225.707.252.3944.9489.4612.5371.6/4.3
8 277.136.063.2435.94>1009.4462.0/0.8
9 172.89<5.151.0624.14>10010.7311.9/0.9
10 204.89<5.151.8126.60>10015341.9/2.0
11 366.339.01--68.544.9971.0/8.5
12 176.9610.17--46.416.8600/4.5
13 171.622.96-->15010.4380/3.5
14 145.687.34--10012.5302.2/3.5
15 251.891.80--1008.5452.2/4.3
16 173.912.83--1009.4341.5/4.4
17 143.794.89--1009.8320.4/4.0
18 275.076.05--99.7711.4451.8/4.0
19 143.705.79--10013.0393.7/0.8
20161.861.54--10010.7373.4/1.6
21 360.591.29--1008.3603.1/4.9
22 172.932.96--1007.8301.0/2.2
23 465.875.93--145.979.8703.9/0.9
24 214.38.071.6359.25>1007.4713.3/2.5
25 90---->1005.0-1.7/6.3
26 123.6810.831.3392.2542.41101913.8/0.2
27 144.059.40-->1507.0902.1/8.3
28 22-4.031.1-223.17.6501.6/8.3
29 154.5611.841.4157.2764.8912.5210/0.1
30 133.706.69--10013.0251.5/3.7
31 141.392.19--10013.0450/4.3
32 111.482.19--1008.61204.7/2.9
33182.192.71--56.57.7303.1/2.1
T3: thyroxine 3; T4: thyroxine 4; FT3: free thyroxine 3; FT4: free thyroxine 4; TSH: thyroid-stimulating hormone; y: years; -: data lost or not tested.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, R.; Tian, J.-L.; Huang, X.-L.; Song, Y.-Z. Genetic Factors Causing Thyroid Dyshormonogenesis as the Major Etiologies for Primary Congenital Hypothyroidism: Clinical and Genetic Characterization of 33 Patients. J. Clin. Med. 2022, 11, 7313. https://doi.org/10.3390/jcm11247313

AMA Style

Liu R, Tian J-L, Huang X-L, Song Y-Z. Genetic Factors Causing Thyroid Dyshormonogenesis as the Major Etiologies for Primary Congenital Hypothyroidism: Clinical and Genetic Characterization of 33 Patients. Journal of Clinical Medicine. 2022; 11(24):7313. https://doi.org/10.3390/jcm11247313

Chicago/Turabian Style

Liu, Rui, Jing-Li Tian, Xiao-Ling Huang, and Yuan-Zong Song. 2022. "Genetic Factors Causing Thyroid Dyshormonogenesis as the Major Etiologies for Primary Congenital Hypothyroidism: Clinical and Genetic Characterization of 33 Patients" Journal of Clinical Medicine 11, no. 24: 7313. https://doi.org/10.3390/jcm11247313

APA Style

Liu, R., Tian, J. -L., Huang, X. -L., & Song, Y. -Z. (2022). Genetic Factors Causing Thyroid Dyshormonogenesis as the Major Etiologies for Primary Congenital Hypothyroidism: Clinical and Genetic Characterization of 33 Patients. Journal of Clinical Medicine, 11(24), 7313. https://doi.org/10.3390/jcm11247313

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop