Pathogenic variations in MAML2 and MAMLD1 contribute to congenital hypothyroidism due to dyshormonogenesis by regulating the Notch signalling pathway.

BACKGROUND: In several countries, thyroid dyshormonogenesis is more common than thyroid dysgenesis in patients with congenital hypothyroidism (CH). However, known pathogenic genes are limited to those directly involved in hormone biosynthesis. The aetiology and pathogenesis of thyroid dyshormonogenesis remain unknown in many patients. METHODS: To identify additional candidate pathogenetic genes, we performed next-generation sequencing in 538 patients with CH and then confirmed the functions of the identified genes in vitro using HEK293T and Nthy-ori 3.1 cells, and in vivo using zebrafish and mouse model organisms. RESULTS: We identified one pathogenic MAML2 variant and two pathogenic MAMLD1 variants that downregulated canonical Notch signalling in three patients with CH. Zebrafish and mice treated with N-[N-(3,5-difluorophenacetyl)-l-alanyl]-S-phenylglycine t-butylester, a γ-secretase inhibitor exhibited clinical manifestations of hypothyroidism and thyroid dyshormonogenesis. Through organoid culture of primary mouse thyroid cells and transcriptome sequencing, we demonstrated that Notch signalling within thyroid cells directly affects thyroid hormone biosynthesis rather than follicular formation. Additionally, these three variants blocked the expression of genes associated with thyroid hormone biosynthesis, which was restored by HES1 expression. The MAML2 variant exerted a dominant-negative effect on both the canonical pathway and thyroid hormone biosynthesis. MAMLD1 also regulated hormone biosynthesis through the expression of HES3, the target gene of the non-canonical pathway. CONCLUSIONS: This study identified three mastermind-like family gene variants in CH and revealed that both canonical and non-canonical Notch signalling affected thyroid hormone biosynthesis.

Large-scale forward genetic screening of zebrafish affecting thyroid development.

The thyroid follicular cells originate from the foregut endoderm and elucidating which genes and signaling pathways regulate their development is crucial for understanding developmental disorders as well as diseases in adulthood. We exploited unique advantages of the zebrafish model to carry an ENU-based forward mutagenesis screen aiming at identifying genes involved in the development and function of the thyroid follicular cells. ENU is an excellent chemical mutagen due to its high mutation efficiency and an indiscriminate selection of genes. A total of 1606 F2 families from 36 ENU treated founders was raised and embryos from F3 generation were collected at 5dpf to perform the whole embryo in situ hybridization with a cocktail probe of thyroid marker thyroglobulin(tg), pituitary marker thyroid stimulating hormone (tshba) to determine the mutagenic phenotype. Among the 1606 F2 families, 112 F2 mutant families with normal development stages except for thyroid dysfunction were identified and divided into three different groups according to their phenotypic characteristics. Further studies of the mutants are likely to shed more insights into the molecular basis of both the thyroid development and function in the zebrafish and vertebrate.

The mutation screening in candidate genes related to thyroid dysgenesis by targeted next-generation sequencing panel in the Chinese congenital hypothyroidism.

OBJECTIVE: Congenital hypothyroidism (CH) is known to be due to thyroid dyshormonogenesis (DH), which is mostly inherited in an autosomal recessive inheritance pattern or thyroid dysgenesis (TD), whose inheritance pattern is controversial and whose molecular etiology remains poorly understood. DESIGN AND METHODS: The variants in 37 candidate genes of CH, including 25 genes related to TD, were screened by targeted exon sequencing in 205 Chinese patients whose CH cannot be explained by biallelic variants in genes related to DH. The inheritance pattern of the genes was analyzed in family trios or quartets. RESULTS: Of the 205 patients, 83 patients carried at least one variant in 19 genes related to TD, and 59 of those 83 patients harbored more than two variants in distinct candidate genes for CH. Biallelic or de novo variants in the genes related to TD in Chinese patients are rare. We also found nine probands carried only one heterozygous variant in the genes related to TD that were inherited from a euthyroid either paternal or maternal parent. These findings did not support the monogenic inheritance pattern of the genes related to TD in CH patients. Notably, in family trio or quartet analysis, of 36 patients carrying more than two variants in distinct genes, 24 patients carried these variants inherited from both their parents, which indicated that the oligogenic inheritance pattern of the genes related to TD should be considered in CH. CONCLUSIONS: Our study expanded the variant spectrum of the genes related to TD in Chinese CH patients. It is rare that CH in Chinese patients could be explained by monogenic germline variants in genes related to TD. The hypothesis of an oligogenic origin of the CH should be considered.

Tpo knockout in zebrafish partially recapitulates clinical manifestations of congenital hypothyroidism and reveals the involvement of TH in proper development of glucose homeostasis.

Congenital hypothyroidism (CH) is a highly prevalent but treatable neonatal endocrine disorder. Thyroid peroxidase (TPO) catalyzes key reactions in thyroid hormone (TH) synthesis. TPO mutations have been found to underlie approximately 5% of congenital hypothyroidism in Chinese patients with more severe phenotypes, the treatment of whom usually requires a higher dose of L-thyroxine. The Tpo gene of zebrafish has 66% homology with the human TPO gene, and synteny analysis has indicated that it is likely a human TPO ortholog. In this study, we generated a tpo-/- mutant zebrafish line through knockout of tpo with CRISPR/Cas9 and investigated the associated phenotypes. Tpo-/- mutant zebrafish displayed growth retardation; an increased number of thyroid follicular cells; and abnormal extrathyroidal phenotypes including pigmentation defects, erythema in the thoracic region, delayed scale development and failure of swim bladder secondary lobe formation. All these abnormal phenotypes were reversed by 30 nM thyroxine (T4) treatment starting at 1 month of age. Tpo-/- mutants also showed increased glucose levels during larval stages, and the increases were induced at least in part by increasing glucagon and decreasing insulin expression. Our work indicates that tpo-mutant zebrafish may serve as a human congenital hypothyroidism model for studying TPO- and TH-related disease mechanisms.

Upregulation of GBP1 in thyroid primordium is required for developmental thyroid morphogenesis.

PURPOSE: Congenital hypothyroidism (CH) is a common congenital endocrine disorder in humans. CH-related diseases such as athyreosis, thyroid ectopy, and hypoplasia are primarily caused by dysgenic thyroid development. However, the underlying molecular mechanisms remain unknown. METHODS: To identify novel CH candidate genes, 192 CH patients were enrolled, and target sequencing of 21 known CH-related genes was performed. The remaining 98 CH patients carrying no known genes were subjected to exome sequencing (ES). The functions of the identified variants were confirmed using thyroid epithelial cells in vitro and in zebrafish model organisms in vivo. RESULTS: Four pathogenic GBP1 variations from three patients were identified. In zebrafish embryos, gbp1 knockdown caused defective thyroid primordium morphogenesis and hypothyroidism. The thyroid cells were stuck together and failed to dissociate from each other to form individual follicles in gbp1-deficient embryos. Furthermore, defects were restored with wild-type human GBP1 (hGBP1) messenger RNA (mRNA) except for mutated hGBP1 (p.H150Y, p.L187P) overexpression. GBP1 promoted ß-catenin translocation into the cytosol and suppressed the formation of cellular adhesion complexes. Suppression of cell-cell adhesion restored the thyroid primordium growth defect observed in gbp1-deficient zebrafish embryos. CONCLUSION: This study provides further understanding regarding thyroid development and shows that defective cellular remodeling could cause congenital hypothyroidism.

Correlation of DUOX2 residual enzymatic activity with phenotype in congenital hypothyroidism caused by biallelic DUOX2 defects.

DUOX2 is the most frequently mutated gene in patients with congenital hypothyroidism (CH) in China. However, no reliable genotype-phenotype relationship has been found in patients with DUOX2 mutations. In this study, DUOX2 mutations were screened in 266 CH patients, and the enzymatic activity of 89 DUOX2 variants was determined in vitro. Furthermore, the DUOX2 residual activity in 76 CH patients caused by DUOX2 biallelic mutations was calculated. The thyroid-stimulating hormone (TSH) and free thyroxine (FT4) levels were found to be higher and lower in patients with DUOX2 residual activity ≤22%, respectively, compared to patients with residual enzymatic activity >22%. Moreover, we interpreted the pathogenicity of DUOX2 variants by applying the ACMG classification criteria with or without PS3/BS3 evidence. The results indicated that residual DUOX2 enzymatic activity was closely related to the clinical phenotypes of CH patients caused by DUOX2 biallelic mutations. These findings suggest that the residual enzymatic activity of 22% may be a cutoff value for estimating the severity of hypothyroidism in CH patients with biallelic DUOX2 mutations. Well-established functional studies are useful and necessary to evaluate the pathogenicity of DUOX2 variants, improving the accuracy and scope of genetic consultations.

The effect of radioiodine treatment on the characteristics of TRAb in Graves' disease.

BACKGROUND: Graves' disease (GD) is one of the most common autoimmune thyroid diseases (AITDs) in humans, and thyrotropin receptor antibody (TRAb) is a characterized autoantibody in GD. The use of radioactive iodine therapy (RAI) for GD treatment is increasing. OBJECTIVES: We studied the biological properties of TRAb and evaluated the effect of RAI therapy on TRAb in GD patients. METHODS: In total, 225 patients (22 onset GD patients without 131I therapy, 203 GD patients treated with 131I therapy) and 20 healthy individuals as normal controls were included in this study. Clinical assessments were performed, and we examined in vitro the biological properties of TRAb in the 22 onset GD patients and 20 controls as well as 84 GD patients with 131I therapy. RESULTS: Serum TRAb and thyroid peroxidase antibody (TPOAb) levels increased in the initial year of RAI treatment, and both antibodies decreased gradually after one year. After 5 years from radioiodine treatment, TRAb and TPOAb levels decreased in 88% and 65% of GD patients, respectively. The proportion of patients positive for thyroid-stimulatory antibody (TSAb) was significantly higher in the 7-12-month group, and thyroid-blocking antibody (TBAb) levels were elevated after one year in half of the patients who received 131I treatment. CONCLUSIONS: Treatment of GD patients with radioiodine increased TPOAb and TRAb (their main biological properties were TSAbs) within the first year after therapy, and the main biological properties of elevated TRAb were TBAbs after 1 year.

A five-gene panel refines differential diagnosis of thyroid nodules.

BACKGROUND: Molecular testing for oncogenic mutations in fine-needle aspiration has showed high predictive value in identifying malignant lesions from thyroid nodules with indeterminate cytology. METHODS: To figure out an efficient and economical gene panel for most medical institutions in China, we designed a five-gene panel including BRAF/NRAS/KRAS/HRAS/TERT genes and conducted a retrospective study to evaluate the role of this five-gene diagnostic panel in differential diagnosis of thyroid nodules. RESULTS: A total of 665 patients with 695 thyroid nodules were investigated in the current study. The fine-needle aspiration biopsy and surgically separated thyroid tissue specimens were harvested to test BRAF, TERT, NRAS, KRAS, and HRAS mutations. We identified 261 mutations in 665 patients, including 177 V600E mutations in BRAF. Three hundred and sixty-nine patients who underwent thyroid surgery after completion of the initial clinical and cytological evaluation were enrolled in the final analysis. The diagnostic sensitivity, specificity, and accuracy of the combination of FNAB cytology and five-gene detection were 74.7%, 93.8%, and 84.8%, respectively. BRAF V600E and five-gene panel could recognize 46.4% and 53.6% of papillary thyroid carcinoma in the patients with cytologically indeterminate nodules. CONCLUSION: The five-gene panel can effectively improve the sensitivity, negative predictive value, and accuracy of fine-needle aspiration biopsy cytology, especially in the patients with cytologically indeterminate nodules.

Fine mapping MHC associations in Graves' disease and its clinical subtypes in Han Chinese.

BACKGROUND: The classical human leucocyte antigen (HLA) genes were the most important genetic determinant for Graves' disease (GD). The aim of the study was to fine map causal variants of the HLA genes. METHODS: We applied imputation with a Pan-Asian HLA reference panel to thoroughly investigate themajor histocompatibility complex (MHC) associations with GD down to the amino acid level of classical HLA genes in 1468 patients with GD and 1490 controls of Han Chinese. RESULTS: The strongest finding across the HLA genes was the association with HLA-DPß1 position 205 (Pomnibus=2.48×10-33). HLA-DPA1*02:02 was the strongest association among the classical HLA alleles, which was in perfect linkage disequilibrium with HLA-DPα1 residue Met11 (OR=1.90, Pbinary=1.76×10-31). Applying stepwise conditional analysis, we identified amino acid position 205 in HLA-DPß1, position 66 and 99 in HLA-B and position 28 in HLA-DRß1 explain majority of the MHC association to GD risk. We further evaluated risk of two clinical subtypes of GD, namely persistent thyroid stimulating hormone receptor antibody -positive (pTRAb+) group and 'non-persistent TRAb positive' (pTRAb-) group after antithyroid drug therapy. We found that HLA-B residues Lys66-Arg69-Val76 could drive pTRAb- GD risk alone, while HLA-DPß1 position 205, HLA-B position 69 and 199 and HLA-DRß1 position 28 drive pTRAb+ GD risk. The risk heterogeneity between pTRAb+ and pTRAb- GD might be driven by HLA-DPα1 Met11. CONCLUSIONS: Four amino acid positions could account for the associations of MHC with GD in Han Chinese. These distinct HLA association patterns indicated the two subtypes have distinct molecular mechanisms of pathogenesis.

A dense mapping study of six European AITD susceptibility regions in a large Chinese Han Cohort of Graves' disease.

OBJECTIVE: We aimed to investigate the six susceptibility loci of GD identified from European population in Chinese Han population and further to estimate the genetic heterogeneity of them in stratification of our GD patients. DESIGN: Dense mapping studies based on GWAS. PATIENTS: A total of 1536 GD patients and 1516 controls in GWAS stage and 1994 GD patients and 2085 controls and 5033 GD patients and 5389 controls in two replication stages. MEASUREMENTS: Based on our previous GWAS data, independently GD-associated SNPs in each region were identified by TagSNP analysis and logistic regression analysis. The association of these SNPs was investigated in 1994 GD patients and 2085 controls, and then, the significantly associated SNPs (P < 0.05) were further genotyped in a second cohort including 5033 GD patients and 5389 controls. RESULTS: After the first replication stage, four SNPs from three regions with Pfirst  < 0.05 were further selected and genotyped in another independent cohort. The association of two SNPs with GD was confirmed in combined Chinese cohorts: rs12575636 at 11q21 (Pcombined  = 7.55 × 10-11 , OR = 1.27) and rs1881145 in TRIB2 at 2p25.1 (Pcombined  = 5.59 × 10-8 , OR = 1.14). Further study disclosed no significant difference for these SNPs between GD subsets. However, eQTL data revealed that SESN3 could be a potential susceptibility gene of GD in 11q21 region. CONCLUSIONS: Out of the six susceptibility loci of GD identified from European population, two risk loci were confirmed in a large Chinese Han population. There is variability in GD genetic susceptibility in different ethnic groups. SESN3 is a potential susceptible gene of GD in 11q21.

Multiregion sequencing reveals the intratumor heterogeneity of driver mutations in TP53-driven non-small cell lung cancer.

Intratumor heterogeneity (ITH) in non-small cell lung cancer (NSCLC) may account for resistance after a period of targeted therapies because drugs destroy only a portion of tumor cells. The recognition of ITH helps identify high-risk patients to make effective treatment decisions. However, ITH studies are confounded by interpatient heterogeneity in NSCLC and a large amount of passenger mutations. To address these issues, we recruited NSCLC patients carrying TP53 mutations and selected driver mutations within recurrently mutated genes in NSCLC. A total of 12-paired normal-tumor tissues were subjected to whole-genome/whole-exome sequencing. From these, 367 non-silent mutations were selected as driver mutations and deeply sequenced in 61 intratumoral microdissections. We identified a universal prevalence of heterogeneity in all 12 tumors, indicating branched evolution. Although TP53 mutations were observed in single biopsy of all 12 tumors, most tumors consist of both TP53 mutated and non-mutated cells in separate regions within the same tumor. This suggests the late molecular timing of the acquisition of TP53 mutations; therefore, the detection of TP53 mutations in a single biopsy may simply not reflect the early malignant potential. In addition, we identified regions of loss of heterozygosity surrounding TP53 and CDKN2A mutations in tumor 711, which also exhibited heterogeneity in different regional samples. Because the ITH of driver mutations likely has clinical consequences, further efforts are needed to limit the impact of ITH and to improve therapeutic efficiency, which will benefit NSCLC patients receiving targeted treatments.

Genome-wide association study identifies a novel susceptibility gene for serum TSH levels in Chinese populations.

Thyroid-stimulating hormone (TSH) is a sensitive indicator of thyroid function. High and low TSH levels reflect hypothyroidism and hyperthyroidism, respectively. Even within the normal range, small differences in TSH levels, on the order of 0.5-1.0 mU/l, are associated with significant differences in blood pressure, BMI, dyslipidemia, risk of atrial fibrillation and atherosclerosis. Most of the variance in TSH levels is thought to be genetically influenced. We conducted a genome-wide association study of TSH levels in 1346 Chinese Han individuals. In the replication study, we genotyped four candidate SNPs with the top association signals in an independent isolated Chinese She cohort (n = 3235). We identified a novel serum TSH susceptibility locus within XKR4 at 8q12.1 (rs2622590, Pcombined = 2.21 × 10(-10)), and we confirmed two previously reported TSH susceptibility loci near FOXE1 at 9q22.33 and near CAPZB at 1p36.13, respectively. The rs2622590_T allele at XKR4 and the rs925489_C allele near FOXE1 were correlated with low TSH levels and were found to be nominally associated to patients with papillary thyroid carcinoma (PTC) (OR = 1.41, P= 0.014 for rs2622590_T, and OR = 1.61, P= 0.030 for rs925489_C). The rs2622590 and rs925489 genotypes were also correlated with the expression levels of FOXE1 and XKR4, respectively, in PTC tissues (P = 2.41 × 10(-4) and P= 0.02). Our findings suggest that the SNPs in XKR4 and near FOXE1 are involved in the regulation of TSH levels.

Robust evidence for five new Graves' disease risk loci from a staged genome-wide association analysis.

Graves' disease (GD), characterized by autoantibodies targeting antigens specifically expressed in thyroid tissues causing hyperthyroidism, is triggered by a combination of genetic and environmental factors. However, only a few loci for GD risk were confirmed in the various ethnic groups, and additional genetic determinants have to be detected. In this study, we carried out a three-stage study in 9529 patients with GD and 9984 controls to identify new risk loci for GD and found genome-wide significant associations in the overall populations for five novel susceptibility loci: the GPR174-ITM2A at Xq21.1, C1QTNF6-RAC2 at 22q12.3-13.1, SLAMF6 at 1q23.2, ABO at 9q34.2 and an intergenic region harboring two non-coding RNAs at 14q32.2 and one previous indefinite locus, TG at 8q24.22 (Pcombined < 5 × 10(-8)). The genotypes of corresponding variants at 14q32.2 and 8q24.22 were correlated with the expression levels of C14orf64 and a TG transcript skipping exon 46, respectively. This study increased the number of GD loci with compelling evidence and indicated that non-coding RNAs might be potentially involved in the pathogenesis of GD.

Identification of BACH2 as a susceptibility gene for Graves' disease in the Chinese Han population based on a three-stage genome-wide association study.

The BACH2 gene regulates B cell differentiation and function and has been reported to be a shared susceptibility gene for several autoimmune diseases. Our previous genome-wide association study (GWAS) indicated that several single nucleotide polymorphisms (SNPs) in the BACH2 gene are associated with Graves' disease (GD) in the Chinese Han population; however, the association did not achieve genome-wide significance levels. Recently, this association of BACH2 with GD was confirmed in Caucasians in the UK population, but fine mapping in this region has not yet been reported. Here, we provide a refined analysis of a 331-kb region in the BACH2 gene, which harbors 359 SNPs, using GWAS data from 1,442 GD patients and 1,468 controls. The SNPs rs2474619 and rs9344996 were implied as the independent variants associated with GD by forward and two-locus logistic regression analysis. We genotyped eight out of 10 tagSNPs with P < 1 × 10(-3) in 3,508 GD patients and 3,209 controls, the results also showed that rs2474619 was independently associated with GD in the combined population from GWAS and the second stage (P = 1.81 × 10(-5)). The rs2474619 and rs9344996 were further genotyped in the third stage cohorts, and rs2474619 showed evidence of association with GD at genome-wide significance levels in the combined population (P = 3.28 × 10(-8), odds ratio = 1.13). The association of rs9344996 with GD can be explained by its linkage to rs2474619 in the combined population. Our study clearly demonstrated that BACH2 is a susceptibility gene for GD in the Chinese Han population and further supported rs2474619, in intron 2 of BACH2, is the best association signal with GD. However, the mechanism by which BACH2 confers increased risk of GD requires further study.

An X chromosome-wide association analysis identifies variants in GPR174 as a risk factor for Graves' disease.

BACKGROUND: Graves' disease is a female preponderant autoimmune illness and the contribution of the X chromosome to its risk has long been appreciated. However, no X-linked susceptibility loci have been indentified from recent genome-wide association studies (GWAS). METHODS: We re-examined the X chromosome data from our recent GWAS for Graves' disease by including males that were previously excluded from the X chromosome analyses. The data were analysed using logistic regression analysis including sex as a covariate, and an additive method assuming X chromosome inactivation, implemented in snpMatrix. RESULTS: A cluster of single nucleotide polymorphism (SNPs) at Xq21.1 was found showing association with genome-wide significance, among which rs3827440 was a non-synonymous SNP of GPR174 (P(logistic regression)= 9.52×10(-8); P(snpMatrix)=4.60×10(-9); OR=1.76, 95% CI 1.45 to 2.13). The association was reproduced in an independent sample collection set including 4564 Graves' disease cases and 3968 sex matched controls (combined P(logistic regression)=5.53×10(-21); combined P(snpMatrix)=4.26×10(-22); OR=1.69, 95% CI 1.53 to 1.86). Notably, GPR174 was widely expressed in immune related tissues and rs3827440 genotypes were associated with distinct mRNA levels (p=0.002). GPR174 did not show sex biased gene expression in our expression analysis. Resequencing study suggested the contribution of some rare variants in the GPR174 gene region to disease risk with a collapsing p value of 1.16×10(-3). CONCLUSIONS: The finding of an X-linked risk locus for Graves' disease expands our understanding of the role of the X chromosome in disease susceptibility.

Contactin 6, A Novel Causative Gene for Congenital Hypothyroidism, Mediates Thyroid Hormone Biosynthesis Through Notch Signaling.

Background: Congenital hypothyroidism (CH) is the most common neonatal metabolic disorder. In patients with CH in China, thyroid dyshormonogenesis is more common than thyroid dysgenesis; however, the genetic causes of CH due to thyroid dyshormonogenesis remain largely unknown. Therefore, we aimed at identifying novel candidate causative genes for CH. Methods: To identify novel CH candidate genes, a total of 599 patients with CH were enrolled and next-generation sequencing was performed. The functions of the identified variants were confirmed using HEK293T and FTC-133 cell lines in vitro and in a mouse model organism in vivo. Results: Three pathogenic contactin 6 (CNTN6) variants were identified in two patients with CH. Pedigree analysis showed that CH caused by CNTN6 variants was inherited in an autosomal recessive pattern. The CNTN6 gene was highly expressed in the thyroid in humans and mice. Cntn6 knockout mice presented with thyroid dyshormonogenesis and CH due to the decreased expression of crucial genes for thyroid hormone biosynthesis (Slc5a5, Tpo, and Duox2). All three CNTN6 variants resulted in the blocking of the release of the Notch intracellular domain, which could not translocate into the nucleus, impaired NOTCH1 transcriptional activity, and decreased expression of SLC5A5, TPO, and DUOX2. Further, we found that DTX1 was required for CNTN6 to promote thyroid hormone biosynthesis through Notch signaling. Conclusions: This study demonstrated that CNTN6 is a novel causative gene for CH through the mediation of thyroid hormone biosynthesis via Notch signaling, which provides new insights into the genetic background and mechanisms involved in CH and thyroid dyshormonogenesis.

Identification of Eukaryotic Translation Initiation Factor 4B as a Novel Candidate Gene for Congenital Hypothyroidism.

CONTEXT: Congenital hypothyroidism (CH) is the most common endocrine disorder in neonates, but its etiology is still poorly understood. OBJECTIVE: We performed whole exome sequencing to identify novel causative gene for CH and functional studies to validate its role in the occurrence of CH. METHODS: Whole exome sequencing in 98 CH patients not harboring known CH candidate genes and bioinformatic analysis were performed. Functional analysis was performed using morpholino, a synthetic short antisense oligonucleotide that contains 25 DNA bases on a methylene morpholine backbone, in zebrafish and CRISPR‒Cas9-mediated gene knockout in mice. RESULTS: Eukaryotic translation initiation factor 4B (EIF4B) was identified as the most promising candidate gene. The EIF4B gene was inherited in an autosomal recessive model, and one patient with thyroid dysgenesis carried EIF4B biallelic variants (p.S430F/p.P328L). In zebrafish, the knockdown of eif4ba/b expression caused thyroid dysgenesis and growth retardation. Thyroid hormone levels were significantly decreased in morphants compared with controls. Thyroxine treatment in morphants partially rescued growth retardation. In mice, the homozygous conceptuses of Eif4b+/- parents did not survive. Eif4b knockout embryos showed severe growth retardation, including thyroid dysgenesis and embryonic lethality before E18.5. CONCLUSION: These experimental data supported a role for EIF4B function in the pathogenesis of the hypothyroid phenotype seen in CH patients. Our work indicated that EIF4B was identified as a novel candidate gene in CH. EIF4B is essential for animal survival, but further studies are needed to validate its role in the pathogenesis of CH.

Deficiency of the HGF/Met pathway leads to thyroid dysgenesis by impeding late thyroid expansion.

The mechanisms of bifurcation, a key step in thyroid development, are largely unknown. Here we find three zebrafish lines from a forward genetic screening with similar thyroid dysgenesis phenotypes and identify a stop-gain mutation in hgfa and two missense mutations in met by positional cloning from these zebrafish lines. The elongation of the thyroid primordium along the pharyngeal midline was dramatically disrupted in these zebrafish lines carrying a mutation in hgfa or met. Further studies show that MAPK inhibitor U0126 could mimic thyroid dysgenesis in zebrafish, and the phenotypes are rescued by overexpression of constitutively active MEK or Snail, downstream molecules of the HGF/Met pathway, in thyrocytes. Moreover, HGF promotes thyrocyte migration, which is probably mediated by downregulation of E-cadherin expression. The delayed bifurcation of the thyroid primordium is also observed in thyroid-specific Met knockout mice. Together, our findings reveal that HGF/Met is indispensable for the bifurcation of the thyroid primordium during thyroid development mediated by downregulation of E-cadherin in thyrocytes via MAPK-snail pathway.

TSHR Variant Screening and Phenotype Analysis in 367 Chinese Patients With Congenital Hypothyroidism.

Background: Genetic defects in the human thyroid-stimulating hormone (TSH) receptor (TSHR) gene can cause congenital hypothyroidism (CH). However, the biological functions and comprehensive genotype-phenotype relationships for most TSHR variants associated with CH remain unexplored. We aimed to identify TSHR variants in Chinese patients with CH, analyze the functions of the variants, and explore the relationships between TSHR genotypes and clinical phenotypes. Methods: In total, 367 patients with CH were recruited for TSHR variant screening using whole-exome sequencing. The effects of the variants were evaluated by in-silico programs such as SIFT and polyphen2. Furthermore, these variants were transfected into 293T cells to detect their Gs/cyclic AMP and Gq/11 signaling activity. Results: Among the 367 patients with CH, 17 TSHR variants, including three novel variants, were identified in 45 patients, and 18 patients carried biallelic TSHR variants. In vitro experiments showed that 10 variants were associated with Gs/cyclic AMP and Gq/11 signaling pathway impairment to varying degrees. Patients with TSHR biallelic variants had lower serum TSH levels and higher free triiodothyronine and thyroxine levels at diagnosis than those with DUOX2 biallelic variants. Conclusions: We found a high frequency of TSHR variants in Chinese patients with CH (12.3%), and 4.9% of cases were caused by TSHR biallelic variants. Ten variants were identified as loss-of-function variants. The data suggest that the clinical phenotype of CH patients caused by TSHR biallelic variants is relatively mild. Our study expands the TSHR variant spectrum and provides further evidence for the elucidation of the genetic etiology of CH.

Dense mapping of IL2RA shows no association with Graves' disease in Chinese Han population.

OBJECTIVE: Associations between IL2RA and various autoimmune diseases have been reported in Caucasians. We investigated whether genetic polymorphisms at the IL2RA locus were associated with Graves' disease (GD) in the Chinese Han population. DESIGN: We performed a genome-wide association study (GWAS) in 1 536 GD patients and 1 516 controls. The 1000 Genomes Project data were adopted as references for imputation analysis. After forward and conditional logistic regressions, we found that rs11256313 was the major risk variant in the CD25/IL2RA region. Thus, we further genotyped rs11256313 in a replication cohort with 3 694 GD patients and 3 510 controls using ABI 7900HT TaqMan Real-Time PCR System. RESULTS: Nine single nucleotide polymorphisms (SNPs) in the IL2RA block were nominally associated with GD in our GWAS (0·01 < P < 0·05). After imputation analysis, 13 imputed SNPs in the IL2RA block were weakly associated with GD (P ≤ 0·05). Logistic regression analysis suggested that the imputed rs11256313 could represent the IL2RA block (P = 0·003). However, we failed to replicate the association of rs11256313 in a larger cohort (P = 0·145). A subphenotype analysis of rs11256313 on thyroid hormone receptor antibody (TRAb) and gender showed that there was no association in any of the subphenotype groups (P > 0·05). CONCLUSIONS: The results suggested that common genetic polymorphisms at IL2RA do not exert a significant genetic effect on the development of GD in the Chinese Han population. Previously reported associations between CD25/IL2RA and autoimmune diseases including GD in Caucasians again imply that heterogeneity exists in different ethnic populations.