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.

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.

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.

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.

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.

Hyperglycemia magnifies bupivacaine-induced cell apoptosis triggered by mitochondria dysfunction and endoplasmic reticulum stress.

Nerve cell injury associated with apoptosis plays an important role in the development of diabetic peripheral neuropathy (DPN). However, it remains unclear whether preexisting or potential neurocyte damage induced by hyperglycemia increases sensitivity to local anesthetics. SH-SY5Y cells were pretreated with a high concentration of glucose in vitro, to imitate DPN prior to administration of bupivacaine or placebo. Cell viability and apoptosis were investigated with a CCK-8 assay and flow cytometry, respectively. In addition, mitochondrial membrane potential, reactive oxygen species (ROS), mitochondrially generated ROS, and activity of mitochondrial complexes I and III were studied to explore the molecular mechanism of bupivacaine-induced mitochondrial injury. Grp78 and caspase-12 expression were measured by qRT-PCR and Western blot, representing endoplasmic reticulum (ER) stress. Cell structure was also assessed via transmission electron microscopy. Incubation with bupivacaine decreased the activity of mitochondrial complexes I and III; increased ROS production at cell and mitochondrial levels, mitochondrial potential depolarization, and Grp78 and caspase-12 expression at both transcription and translation levels; and affected cell structure, which could be enhanced by glucose pretreatment. These findings indicate that mitochondrial dysfunction and ER stress related to ROS are involved in bupivacaine-induced apoptosis and may be enhanced by glucose administration.

Spleen stiffness and splenoportal venous flow: assessment before and after transjugular intrahepatic portosystemic shunt placement.

OBJECTIVES: To prospectively assess changes in spleen stiffness and splenoportal venous flow before and after transjugular intrahepatic portosystemic shunt (TIPS) placement. METHODS: We prospectively evaluated spleen stiffness measured by the mean shear wave velocity with acoustic radiation force impulse imaging and the splenoportal venous velocity with color Doppler sonography in 12 patients (mean age ± SD, 42.6 ± 11.0 years; range, 29-65 years) who underwent TIPS placement for portal hypertension and gastroesophageal bleeding. The mean shear wave velocity and angle-corrected splenoportal venous velocity at the main portal and splenic veins were measured 1 day before and 3 to 9 days after TIPS placement (mean interval, 6.0 ± 1.95 days; range, 4-10 days) and were compared with portal vein pressure measured during the procedure. RESULTS: There was a significant difference in portal vein pressure before and after TIPS (25.34 ± 6.21 versus 15.66 ± 6.07 mm Hg; P = .0005). After TIPS, the mean shear wave velocity decreased significantly in all 12 cases (3.50 ± 0.46 versus 3.15 ± 0.39 m/s before and after TIPS; P = .00015). The flow velocity at the main portal vein increased significantly after TIPS (22.21 ± 4.13 versus 47.25 ± 12.37 cm/s; P = .0000051). The splenic vein velocity and spleen index measured 25.57 ± 6.98 cm/s and 55.99 ± 21.27 cm(2), respectively, before TIPS and 35.72 ± 11.10 cm/s and 50.11 ± 21.12 cm(2) after TIPS (P = .0004 and .003). CONCLUSIONS: A significant decrease in the mean shear wave velocity and increase in the splenoportal venous velocity occurred with reduced portal vein pressure after TIPS placement. Hence, both parameters can be used as noninvasive quantitative markers for monitoring TIPS function after placement.

Comparison and development of machine learning for thalidomide-induced peripheral neuropathy prediction of refractory Crohn's disease in Chinese population.

BACKGROUND: Thalidomide is an effective treatment for refractory Crohn's disease (CD). However, thalidomide-induced peripheral neuropathy (TiPN), which has a large individual variation, is a major cause of treatment failure. TiPN is rarely predictable and recognized, especially in CD. It is necessary to develop a risk model to predict TiPN occurrence. AIM: To develop and compare a predictive model of TiPN using machine learning based on comprehensive clinical and genetic variables. METHODS: A retrospective cohort of 164 CD patients from January 2016 to June 2022 was used to establish the model. The National Cancer Institute Common Toxicity Criteria Sensory Scale (version 4.0) was used to assess TiPN. With 18 clinical features and 150 genetic variables, five predictive models were established and evaluated by the confusion matrix receiver operating characteristic curve (AUROC), area under the precision-recall curve (AUPRC), specificity, sensitivity (recall rate), precision, accuracy, and F1 score. RESULTS: The top-ranking five risk variables associated with TiPN were interleukin-12 rs1353248 [P = 0.0004, odds ratio (OR): 8.983, 95% confidence interval (CI): 2.497-30.90], dose (mg/d, P = 0.002), brain-derived neurotrophic factor (BDNF) rs2030324 (P = 0.001, OR: 3.164, 95%CI: 1.561-6.434), BDNF rs6265 (P = 0.001, OR: 3.150, 95%CI: 1.546-6.073) and BDNF rs11030104 (P = 0.001, OR: 3.091, 95%CI: 1.525-5.960). In the training set, gradient boosting decision tree (GBDT), extremely random trees (ET), random forest, logistic regression and extreme gradient boosting (XGBoost) obtained AUROC values > 0.90 and AUPRC > 0.87. Among these models, XGBoost and GBDT obtained the first two highest AUROC (0.90 and 1), AUPRC (0.98 and 1), accuracy (0.96 and 0.98), precision (0.90 and 0.95), F1 score (0.95 and 0.98), specificity (0.94 and 0.97), and sensitivity (1). In the validation set, XGBoost algorithm exhibited the best predictive performance with the highest specificity (0.857), accuracy (0.818), AUPRC (0.86) and AUROC (0.89). ET and GBDT obtained the highest sensitivity (1) and F1 score (0.8). Overall, compared with other state-of-the-art classifiers such as ET, GBDT and RF, XGBoost algorithm not only showed a more stable performance, but also yielded higher ROC-AUC and PRC-AUC scores, demonstrating its high accuracy in prediction of TiPN occurrence. CONCLUSION: The powerful XGBoost algorithm accurately predicts TiPN using 18 clinical features and 14 genetic variables. With the ability to identify high-risk patients using single nucleotide polymorphisms, it offers a feasible option for improving thalidomide efficacy in CD patients.

The isl2a transcription factor regulates pituitary development in zebrafish.

Background: ISL LIM homeobox 2, also known as insulin gene enhancer protein ISL-2 (ISL2), is a transcription factor gene that participates in a wide range of developmental events. However, the role of ISL2 in the hypothalamus-pituitary-thyroid axis is largely unknown. In the present study, we characterized the expression patterns of ISL2 and revealed its regulative role during embryogenesis using zebrafish. Methods: We used the CRISPR/Cas9 system to successfully establish homozygous ISL2-orthologue (isl2a and isl2b) knockout zebrafish. Moreover, we utilized these knockout zebrafish to analyze the pituitary and thyroid phenotypes in vivo. For further molecular characterization, in situ hybridization and immunofluorescence were performed. Results: The isl2a mutant zebrafish presented with thyroid hypoplasia, reduced whole-body levels of thyroid hormones, increased early mortality, gender imbalance, and morphological retardation during maturity. Additionally, thyrotropes, a pituitary cell type, was notably decreased during development. Importantly, the transcriptional levels of pituitary-thyroid axis hormones-encoding genes, such as tshba, cga, and tg, were significantly decreased in isl2a mutants. Finally, the thyroid dysplasia in isl2a mutant larvae may be attributed to a reduction in proliferation rather than changes in apoptosis. Conclusions: In summary, isl2a regulates the transcriptional levels of marker genes in hypothalamus-pituitary-thyroid axis, and isl2a knockout causing low thyroid hormone levels in zebrafish. Thus, isl2a identified by the present study, is a novel regulator for pituitary cell differentiation in zebrafish, resulting in thyroid gland hypoplasia and phenotypes of hypothyroidism.

Liver and spleen stiffness measured by acoustic radiation force impulse elastography for noninvasive assessment of liver fibrosis and esophageal varices in patients with chronic hepatitis B.

OBJECTIVES: To evaluate the performance of liver and spleen stiffness measured by acoustic radiation force impulse (ARFI) elastography for noninvasive assessment of liver fibrosis and esophageal varices in patients with chronic hepatitis B virus. METHODS: Two hundred sixty-four participants, of whom 60 were healthy volunteers (classified as stage 0), 66 were patients with chronic hepatitis B who had undergone liver biopsy, and 138 were patients with hepatitis B-related cirrhosis, were enrolled in this study. Median liver and spleen stiffness values (meters per second) from 10 successful measurements per participant were obtained. Patients with cirrhosis were examined by upper endoscopy. RESULTS: Significant linear correlations were found between liver (Spearman ρ = 0.87; P < .001) and spleen (Spearman ρ = 0.76; P < .001) stiffness and the fibrosis stage. Liver and spleen stiffness values increased as fibrosis progressed; however, overlaps in liver stiffness were detected in stages 0 and 1 and 1 and 2, and overlaps in spleen stiffness were observed in stages 0 and 1, 1 and 2, and 2 and 3. Liver stiffness cutoff values were 1.69 m/s for predicting stage 3 or greater (area under the receiver operating characteristic curve [AUROC] = 0.99) and 1.88 m/s for stage 4 (AUROC = 0.97). The spleen stiffness cutoff value was 2.72 m/s for stage 4 (AUROC = 0.96). Liver stiffness was not correlated with the varix grade, whereas a significant linear correlation (Spearman ρ = 0.65; P < .001) between spleen stiffness and the varix grade was found. The optimal spleen stiffness cutoff value for predicting varices was 3.16 m/s (AUROC = 0.83). CONCLUSIONS: Liver and spleen stiffness values measured by ARFI elastography are reliable predictors of liver fibrosis. Spleen stiffness measured by ARFI can be used as a non-invasive method for determining the presence and severity of esophageal varices; however, evidence to support a similar role for liver stiffness is lacking.

Lymphocyte infiltration and thyrocyte destruction are driven by stromal and immune cell components in Hashimoto's thyroiditis.

Hashimoto's thyroiditis (HT) is the most common autoimmune disease characterized by lymphocytic infiltration and thyrocyte destruction. Dissection of the interaction between the thyroidal stromal microenvironment and the infiltrating immune cells might lead to a better understanding of HT pathogenesis. Here we show, using single-cell RNA-sequencing, that three thyroidal stromal cell subsets, ACKR1+ endothelial cells and CCL21+ myofibroblasts and CCL21+ fibroblasts, contribute to the thyroidal tissue microenvironment in HT. These cell types occupy distinct histological locations within the thyroid gland. Our experiments suggest that they might facilitate lymphocyte trafficking from the blood to thyroid tissues, and T cell zone CCL21+ fibroblasts may also promote the formation of tertiary lymphoid organs characteristic to HT. Our study also demonstrates the presence of inflammatory macrophages and dendritic cells expressing high levels of IL-1ß in the thyroid, which may contribute to thyrocyte destruction in HT patients. Our findings thus provide a deeper insight into the cellular interactions that might prompt the pathogenesis of HT.

Detection of BRAF V600E in Fine-Needle Aspiration Samples of Thyroid Nodules by Droplet Digital PCR.

Background: BRAF exon 15 p.V600E (BRAF V600E) mutation has been established as an important molecular marker for papillary thyroid carcinoma diagnosis by ultrasound-guided fine-needle aspiration biopsy (FNAB). Sanger sequencing is the gold standard for detecting BRAF V600E mutations but fails to identify low-frequency mutations. However, droplet digital PCR (ddPCR) is a popular new method for detecting low-frequency mutations. Here, we compare the efficiency of droplet digital PCR (ddPCR) and Sanger sequencing for detection of the BRAF V600E mutation in thyroid fine-needle aspiration (FNA) samples. Methods: Thyroid fine-needle aspiration samples from 278 patients with 310 thyroid nodules were collected. Sanger sequencing and ddPCR were conducted to detect the BRAF V600E mutation. Results: The BRAF V600E mutation was found in 94 nodules (30.32%) by ddPCR and 40 nodules (12.90%) by Sanger sequencing in 310 FNA samples. A total of 119 nodules were confirmed PTC by postsurgical pathology. Among which the BRAF mutation was found in 80 (67.23%) nodules by ddPCR and 31 (26.05%) by Sanger sequencing. All nodules carrying the mutation detected by Sanger sequencing (SS+) were verified by ddPCR (ddPCR+). Also, all nodules with no mutation detected by ddPCR were interpreted as wild-type by Sanger sequencing (SS-). In addition. Almost all SS+/ddPCR + nodules (95.00%; 38/40) and SS-/ddPCR + nodules (100.00%; 54/54) displayed a BRAF mutation rate of >5% and <15%, respectively, indicating easy misdetection by Sanger sequencing when the mutation rate is between 5 and 15%. Conclusion: ddPCR has higher sensitivity than Sanger sequencing and we propose ddPCR as a supplement to Sanger sequencing in molecular testing of BRAF using FNAB samples.

Three-dimensional microscopy and image fusion reconstruction analysis of the thyroid gland during morphogenesis.

Thyroid dysgenesis (TD) is a major cause of primary congenital hypothyroidism; however, the molecular mechanism underlying this process is unclear. Current knowledge regarding the morphogenesis of the thyroid gland and vascular anomalies affecting thyroid development is limited. To monitor the early stages of thyroid gland development, we generated double transgenic zebrafish embryos Tg(tg:mCherry/flk1:EGFP). We described the volume of the thyroid from 2 days postfertilization (dpf) to 5 dpf using 3D reconstruction images. We treated zebrafish embryos with the fibroblast growth factor (FGF) inhibitor PD166866 to better understand the impact of vascular defects on thyroid development and the effects of drug administration at specific time periods on different stages of thyroid development. The 3D reconstruction data revealed that the thyroid glands underwent significant transformation at critical time points. PD166866 treatment from 48 to 72 hours postfertilization (hpf) and from 72 to 96 hpf did not cause obvious reductions in thyroid volume but did result in observable abnormalities in thyroid morphology. The treatment also affected thyroid volume from 36 to 48 hpf, thus indicating that there are time-point-specific effects of drug administration during thyroid development. Three-dimensional image reconstruction provides a comprehensive picture of thyroid anatomy and can be used to complement anatomical fluorescence information. The effects of an FGF pathway inhibitor on thyroid development were determined to be time-point-dependent.

Compelling Evidence Linking CD40 Gene With Graves' Disease in the Chinese Han Population.

Mutations in CD40 have been widely reported to be risk factors for Graves' disease (GD). The gene, along with its cognate ligand CD40L, may regulate pro-inflammatory and immune responses. Rs1883832, located at the -1 position of the Kozak sequence, is the most well-studied single nucleotide polymorphism (SNP) of CD40, and has been confirmed to predispose those with the alteration to GD, regardless of ethnicity. Our genome-wide association study (GWAS) indicated that several SNPs, including rs1883832 located within the vicinity of CD40 were associated with GD in the Han Chinese population. Aiming at identifying the most consequential SNP and its underlying pathogenic mechanism, we performed a two-stage refined study on 8,171 patients with GD and 7,906 controls, and found rs1883832 was the most significantly GD-associated SNP in the CD40 gene region (PCombined = 9.17×10-11, OR = 1.18). Through searching the cis-expression quantitative trait locus database and using quantitative RT-PCR, we further discovered that the rs1883832 genotype can influence CD40 gene transcription. Furthermore, we demonstrated that rs1883832 is a susceptibility locus for pTRAb+ GD patients. In conclusion, the current study provides robust evidence that rs1883832 can regulate CD40 gene expression and affect serum TRAb levels, which ultimately contributes to the development of GD.

Genetic Study in a Large Cohort Supported Different Pathogenesis of Graves' Disease and Hashimoto's Hypothyroidism.

CONTEXT: Hashimoto's thyroiditis (HT) and Graves' disease (GD) are the 2 main autoimmune thyroid diseases that have both similarities and differences. Determining the genetic basis that distinguishes HT from GD is key for a better understanding of the differences between these closely related diseases. OBJECTS: To identify the susceptibility genes for HT in the Chinese cohort and compare susceptibility genes between GD and HT. DESIGN: In the current study, 18 SNPs from 18 established GD risk loci were selected and then genotyped in 2682 patients with HT, 4980 patients with GD, and 3892 controls. The association analysis between HT and controls and heterogeneity analysis between HT and GD were performed on SPSS, with the logistic regression analysis adjusted for sex and age. RESULTS: We identified 11 susceptibility loci for HT in the Chinese Han population, with 4 loci, including the rs1265883 in SLAMF6 locus, rs1024161 in CTLA4, rs1521 in HLA-B, and rs5912838 in GPR174/ ITM2A at X chromosome, reaching genome-wide significance of 5 × 10-8. Five loci were reported to be associated with HT for the first time. We also identified 6 susceptibility loci with heterogeneity between GD and HT. Out of them, 4 loci were associated with GD but not with HT, including HLA-DPB1, CD40, TSHR, and TG; the association of HLA-B with GD was stronger than that with HT, but the association of SLAMF6 was reversed. CONCLUSION: Our findings suggested that the pathogenesis of HT and GD was different.

A Weighted Genetic Risk Score Using Known Susceptibility Variants to Predict Graves Disease Risk.

CONTEXT: Graves disease (GD) is a common thyroid-specific autoimmune disease and one of the most heritable diseases in the population. We present a risk-prediction model, including confirmed, known genetic variants associated with GD. DESIGN: To construct a stable-prediction model, we used known GD susceptibility single nucleotide polymorphisms (SNPs) as markers and trained and tested our model in a cohort of 4897 patients with GD and 5098 healthy controls. We weighted the contribution of each SNP to the disease to calculate the weighted genetic risk score (wGRS) for each individual. The efficiency of this model can be estimated by the area under the curve (AUC) receiver operator characteristic curve and the specificity and sensitivity of each wGRS. RESULTS: With the 20 confirmed GD risk-related SNPs, our wGRS-prediction model could predict patients with GD from the general population (AUC 0.70 [95% CI: 0.69 to 0.71]) and did especially well in predicting patients with GD with persisting thyroid-stimulating hormone receptor antibody positive [pTRAb+; AUC 0.74 (95% CI: 0.72 to 0.76)]. We also evaluated how the four pTRAb+ specific risk SNPs predicted patients with GD with pTRAb+ among all patients with GD [AUC 0.62 (95% CI: 0.61 to 0.63)]. For clinical use, we partitioned subjects in each set into different risk categories to generate the wGRS cutoff of high risk for reference. CONCLUSIONS: Our study provides an approach to predict GD risk in the general population by the calculation of the wGRS of 20 known GD susceptibility variants. The wGRS-prediction model was more stable and convenient, whereas the prediction performance was still modest.

Uterus globulin associated protein 1 (UGRP1) is a potential marker of progression of Graves' disease into hypothyroidism.

Approximately 20% of Graves' disease (GD) patients may result eventually in hypothyroidism in their natural course. Uterus globulin-associated protein 1 (UGRP1) was associated with GD in our previous study. Here we investigated the role of UGRP1 in the development of autoimmune thyroid disease (AITD). The results showed that UGRP1 was expressed in the thyrocytes of most Hashimoto's thyroiditis (HT) patients and a proportion of GD patients (293 HT and 198 GD). The pathologic features of UGRP1-positive thyrocytes resembled "Hürthle cells", and were surrounded by infiltrated leukocytes. The positivity rate of TPOAb in UGRP1-positive GD patients was much higher than that in -negative GD patients. Moreover, UGRP1 was co-expressed with Fas and HLA-DR in the thyrocytes of AITD patients. We also found IL-1ß but not Th1 or Th2 cytokines was able to upregulate the expression of UGRP1. Our findings indicated that UGRP1 may be a novel marker in thyrocytes to predict GD patients who develop hypothyroidism.

Assessment of Molecular Subtypes in Thyrotoxic Periodic Paralysis and Graves Disease Among Chinese Han Adults: A Population-Based Genome-Wide Association Study.

Importance: Thyrotoxic periodic paralysis (TPP) is a potentially lethal complication of hyperthyroidism. However, only 1 specific susceptibility locus for TPP has been identified. Additional genetic determinants should be detected so that a prediction model can be constructed. Objective: To investigate the genetic architecture of TPP and distinguish TPP from Graves disease cohorts. Design, Setting, and Participants: This population-based case-control study used a 2-stage genome-wide association study to investigate the risk loci of TPP and weighted genetic risk score to construct a TPP prediction model with data from a Chinese Han population recruited in hospitals in China from March 2003 to December 2015. The analysis was conducted from November 2014 to August 2016. Main Outcomes and Measures: Loci specifically associated with TPP risk and those shared with Graves disease and prediction model of joint effects of TPP-specific loci. Results: A total of 537 patients with TPP (mean [SD] age, 35 [11] years; 458 male) 1519 patients with Graves disease and no history of TPP (mean [SD] age, 38 [13] years; 366 male), and 3249 healthy participants (mean [SD] age, 46 [10] years; 1648 male) were recruited from the Han population by hospitals throughout China. Two new TPP-specific susceptibility loci were identified: DCHS2 on 4q31.3 (rs1352714: odds ratio [OR], 1.58; 95% CI, 1.35-1.85; P = 1.24 × 10-8) and C11orf67 on 11q14.1 (rs2186564: OR, 1.50; 95% CI, 1.29-1.74; P = 2.80 × 10-7). One previously reported specific locus was confirmed on 17q24.3 near KCNJ2 (rs312729: OR, 2.08; 95% CI, 1.83-2.38; P = 8.02 × 10-29). Meanwhile, 2 risk loci (MHC and Xq21.1) were shared by Graves disease and TPP. After 2 years of treatment, the ratio of persistent thyrotropin receptor antibody positivity was higher in patients with TPP than in patients with Graves disease and no history of TPP (OR, 3.82; 95% CI, 2.04-7.16; P = 7.05 × 10-6). The prediction model using a weighted genetic risk score and 11 candidate TPP-specific single-nucleotide polymorphisms had an area under the curve of 0.80. Conclusions and Relevance: These findings provide evidence that TPP is a novel molecular subtype of Graves disease. The newly identified loci, along with other previously reported loci, demonstrate the growing complexity of the heritable contribution to TPP pathogenesis. A complete genetic architecture will be helpful to understand the pathophysiology of TPP, and a useful prediction model could prevent the onset of TPP.

The genetic characteristics of congenital hypothyroidism in China by comprehensive screening of 21 candidate genes.

OBJECTIVE: Congenital hypothyroidism (CH), the most common neonatal metabolic disorder, is characterized by impaired neurodevelopment. Although several candidate genes have been associated with CH, comprehensive screening of causative genes has been limited. DESIGN AND METHODS: One hundred ten patients with primary CH were recruited in this study. All exons and exon-intron boundaries of 21 candidate genes for CH were analyzed by next-generation sequencing. And the inheritance pattern of causative genes was analyzed by the study of family pedigrees. RESULTS: Our results showed that 57 patients (51.82%) carried biallelic mutations (containing compound heterozygous mutations and homozygous mutations) in six genes (DUOX2, DUOXA2, DUOXA1, TG, TPO and TSHR) involved in thyroid hormone synthesis. Autosomal recessive inheritance of CH caused by mutations in DUOX2, DUOXA2, TG and TPO was confirmed by analysis of 22 family pedigrees. Notably, eight mutations in four genes (FOXE1, NKX2-1, PAX8 and HHEX) that lead to thyroid dysgenesis were identified in eight probands. These mutations were heterozygous in all cases and hypothyroidism was not observed in parents of these probands. CONCLUSIONS: Most cases of congenital hypothyroidism in China were caused by thyroid dyshormonogenesis rather than thyroid dysgenesis. This study identified previously reported causative genes for 57/110 Chinese patients and revealed DUOX2 was the most frequently mutated gene in these patients. Our study expanded the mutation spectrum of CH in Chinese patients, which was significantly different from Western countries.