New insights into the hypothalamic-pituitary-thyroid axis: a transcriptome- and proteome-wide association study.

Introduction: Thyroid hormones have systemic effects on the human body and play a key role in the development and function of virtually all tissues. They are regulated via the hypothalamic-pituitary-thyroid (HPT) axis and have a heritable component. Using genetic information, we applied tissue-specific transcriptome-wide association studies (TWAS) and plasma proteome-wide association studies (PWAS) to elucidate gene products related to thyrotropin (TSH) and free thyroxine (FT4) levels. Results: TWAS identified 297 and 113 transcripts associated with TSH and FT4 levels, respectively (25 shared), including transcripts not identified by genome-wide association studies (GWAS) of these traits, demonstrating the increased power of this approach. Testing for genetic colocalization revealed a shared genetic basis of 158 transcripts with TSH and 45 transcripts with FT4, including independent, FT4-associated genetic signals within the CAPZB locus that were differentially associated with CAPZB expression in different tissues. PWAS identified 18 and ten proteins associated with TSH and FT4, respectively (HEXIM1 and QSOX2 with both). Among these, the cognate genes of five TSH- and 7 FT4-associated proteins mapped outside significant GWAS loci. Colocalization was observed for five plasma proteins each with TSH and FT4. There were ten TSH and one FT4-related gene(s) significant in both TWAS and PWAS. Of these, ANXA5 expression and plasma annexin A5 levels were inversely associated with TSH (PWAS: P = 1.18 × 10-13, TWAS: P = 7.61 × 10-12 (whole blood), P = 6.40 × 10-13 (hypothalamus), P = 1.57 × 10-15 (pituitary), P = 4.27 × 10-15 (thyroid)), supported by colocalizations. Conclusion: Our analyses revealed new thyroid function-associated genes and prioritized candidates in known GWAS loci, contributing to a better understanding of transcriptional regulation and protein levels relevant to thyroid function.

The role of thyroid function in borderline personality disorder and schizophrenia: a Mendelian Randomisation study.

BACKGROUND: Genome-wide association studies have reported a genetic overlap between borderline personality disorder (BPD) and schizophrenia (SCZ). Epidemiologically, the direction and causality of the association between thyroid function and risk of BPD and SCZ are unclear. We aim to test whether genetically predicted variations in TSH and FT4 levels or hypothyroidism are associated with the risk of BPD and SCZ. METHODS: We employed Mendelian Randomisation (MR) analyses using genetic instruments associated with TSH and FT4 levels as well as hypothyroidism to examine the effects of genetically predicted thyroid function on BPD and SCZ risk. Bidirectional MR analyses were employed to investigate a potential reverse causal association. RESULTS: Genetically predicted higher FT4 was not associated with the risk of BPD (OR: 1.18; P = 0.60, IVW) or the risk of SCZ (OR: 0.93; P = 0.19, IVW). Genetically predicted higher TSH was not associated with the risk of BPD (OR: 1.11; P = 0.51, IVW) or SCZ (OR: 0.98, P = 0.55, IVW). Genetically predicted hypothyroidism was not associated with BPD or SCZ. We found no evidence for a reverse causal effect between BPD or SCZ on thyroid function. CONCLUSIONS: We report evidence for a null association between genetically predicted FT4, TSH or hypothyroidism with BPD or SCZ risk. There was no evidence for reverse causality.

Towards Personalized TSH Reference Ranges: A Genetic and Population-Based Approach in Three Independent Cohorts.

Background: Serum thyroid-stimulating hormone (TSH) measurement is the diagnostic cornerstone for primary thyroid dysfunction. There is high inter-individual but limited intra-individual variation in TSH concentrations, largely due to genetic factors. The currently used wide population-based reference intervals may lead to inappropriate management decisions. Methods: A polygenic score (PGS) including 59 genetic variants was used to calculate genetically determined TSH reference ranges in a thyroid disease-free cohort (n = 6,834). Its effect on reclassification of diagnoses was investigated when compared to using population-based reference ranges. Next, results were validated in a second independent population-based thyroid disease-free cohort (n = 3,800). Potential clinical implications were assessed in a third independent population-based cohort including individuals without thyroid disease (n = 26,321) as well as individuals on levothyroxine (LT4) treatment (n = 1,132). Results: PGS was a much stronger predictor of individual TSH concentrations than FT4 (total variance in TSH concentrations explained 9.2-11.1% vs. 2.4-2.7%, respectively) or any other nongenetic factor (total variance in TSH concentrations explained 0.2-1.8%). Genetically determined TSH reference ranges differed significantly between PGS quartiles in all cohorts, while the differences in FT4 concentrations were absent or only minor. Up to 24.7-30.1% of individuals, previously classified as having subclinical hypo- and hyperthyroidism when using population-based TSH reference ranges, were reclassified as euthyroid when genetically determined TSH reference ranges were applied. Individuals in the higher PGS quartiles had a higher probability of being prescribed LT4 treatment compared to individuals from the lower PGS quartiles (3.3% in Q1 vs. 5.2% in Q4, Pfor trend =1.7 × 10-8). Conclusions: Individual genetic profiles have the potential to personalize TSH reference ranges, with large effects on reclassification of diagnosis and LT4 prescriptions. As the currently used PGS can only predict approximately 10% of inter-individual variation in TSH concentrations, it should be further improved when more genetic variants determining TSH concentrations are identified in future studies.

Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications.

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Novel (sulfated) thyroid hormone transporters in the solute carrier 22 family.

Objective: Thyroid hormone (TH) transport represents a critical first step in governing intracellular TH regulation. It is still unknown whether the full repertoire of TH transporters has been identified. Members of the solute carrier (SLC) 22 family have substrates in common with the known TH transporters of the organic anion-transporting peptide family. Therefore, we screened the SLC22 family for TH transporters. Methods: Uptake of 1 nM of iodothyronines or sulfated iodothyronines in COS1 cells expressing SLC22 proteins was performed. Results: We first tested 25 mouse (m) SLC22 proteins for TH uptake and found that the majority of the organic anion transporter (OAT) clade were capable of 3,3',5-triiodothyronine and/or thyroxine (T4) transport. Based on phylogenetic tree analysis of the mouse and human (h) SLC22 family, we selected eight hSLC22s that grouped with the newly identified mouse TH transporters. Of these, four tested positive for uptake of one or more substrates, particularly hSLC22A11 showed robust (3-fold over control) uptake of T4. Uptake of sulfated iodothyronines was strongly (up to 17-fold) induced by some SLC22s, most notably SLC22A8, hSLC22A9, mSLC22A27 and mSLC22A29. Finally, the zebrafish orthologues of SLC22A6/8 drOatx and drSlc22a6l also transported almost all (sulfated) iodothyronines tested. The OAT inhibitors lesinurad and probenecid inhibited most SLC22 proteins. Conclusions: Our results demonstrated that members of the OAT clade of the SLC22 family constitute a novel, evolutionary conserved group of transporters for (sulfated) iodothyronines. Future studies should reveal the relevance of these transporters in TH homeostasis and physiology.

The Effects of Common Genetic Variation in 96 Genes Involved in Thyroid Hormone Regulation on TSH and FT4 Concentrations.

OBJECTIVE: While most of the variation in thyroid function is determined by genetic factors, single nucleotide polymorphisms (SNPs) identified via genome-wide association analyses have only explained ~5% to 9% of this variance so far. Most SNPs were in or nearby genes with no known role in thyroid hormone (TH) regulation. Therefore, we performed a large-scale candidate gene study investigating the effect of common genetic variation in established TH regulating genes on serum thyrotropin [thyroid-stimulating hormone (TSH)] and thyroxine (FT4) concentrations. METHODS: SNPs in or within 10 kb of 96 TH regulating genes were included (30 031 TSH SNPs, and 29 962 FT4 SNPs). Associations were studied in 54 288 individuals from the ThyroidOmics Consortium. Linkage disequilibrium-based clumping was used to identify independently associated SNPs. SNP-based explained variances were calculated using SumHer software. RESULTS: We identified 23 novel TSH-associated SNPs in predominantly hypothalamic-pituitary-thyroid axis genes and 25 novel FT4-associated SNPs in mainly peripheral metabolism and transport genes. Genome-wide SNP variation explained ~21% (SD 1.7) of the total variation in both TSH and FT4 concentrations, whereas SNPs in the 96 TH regulating genes explained 1.9% to 2.6% (SD 0.4). CONCLUSION: Here we report the largest candidate gene analysis on thyroid function, resulting in a substantial increase in the number of genetic variants determining TSH and FT4 concentrations. Interestingly, these candidate gene SNPs explain only a minor part of the variation in TSH and FT4 concentrations, which substantiates the need for large genetic studies including common and rare variants to unravel novel, yet unknown, pathways in TH regulation.

Thyroid Function and Mood Disorders: A Mendelian Randomization Study.

Background: Observational studies suggest that even minor variations in thyroid function are associated with the risk of mood disorders, including major depressive disorder (MDD) and bipolar disorder (BD). However, it is unknown whether these associations are causal or not. We used a Mendelian randomization (MR) approach to investigate causal effects of minor variations in thyrotropin (TSH) and free thyroxine (fT4) levels on MDD and BD risk. Materials and Methods: We performed two-sample MR analyses using data from the largest publicly available genome-wide association studies on normal-range TSH (n = 54,288) and fT4 (n = 49,269) levels, MDD (170,756 cases, 329,443 controls) and BD (20,352 cases, 31,358 controls). Secondary MR analyses investigated the effects of TSH and fT4 levels on specific MDD and BD subtypes. Reverse MR was also performed to assess the effects of MDD and BD on TSH and fT4 levels. Results: There were no associations between genetically predicted TSH and fT4 levels and MDD risk, nor MDD subtypes and minor depressive symptoms. A one standard deviation increase in fT4 levels was nominally associated with an 11% decrease in the overall BD risk (odds ratio [OR] = 0.89, 95% confidence interval [CI] = 0.80-0.98, p = 0.022) and a 13% decrease in the BD type 1 risk (OR = 0.87, CI = 0.75-1.00, p = 0.047). In the reverse direction, genetic predisposition to MDD and BD was not associated with TSH nor fT4 levels. Conclusions: Variations in normal-range TSH and fT4 levels have no effects on the risk of MDD and its subtypes, and neither on minor depressive symptoms. This indicates that depressive symptoms should not be attributed to minor variations in thyroid function. Borderline associations with BD and BD type 1 risks suggest that further clinical studies should investigate the effect of thyroid hormone treatment in BD.

Effects of Thyroid Function on Hemostasis, Coagulation, and Fibrinolysis: A Mendelian Randomization Study.

Background: Untreated hypothyroidism is associated with acquired von Willebrand syndrome, and hyperthyroidism is associated with increased thrombosis risk. However, the causal effects of thyroid function on hemostasis, coagulation, and fibrinolysis are unknown. Methods: In a two-sample Mendelian randomization (MR) study with genome-wide association variants, we assessed causality of genetically predicted hypothyroidism (N = 134,641), normal-range thyrotropin (TSH; N = 54,288) and free thyroxine (fT4) (N = 49,269), hyperthyroidism (N = 51,823), and thyroid peroxidase antibody positivity (N = 25,821) on coagulation (activated partial thromboplastin time, von Willebrand factor [VWF], factor VIII [FVIII], prothrombin time, factor VII, fibrinogen) and fibrinolysis (D-dimer, tissue plasminogen activator [TPA], plasminogen activator inhibitor-1) from the CHARGE Hemostasis Consortium (N = 2583-120,246). Inverse-variance-weighted random effects were the main MR analysis followed by sensitivity analyses. Two-sided p < 0.05 was nominally significant, and p < 0.0011[ = 0.05/(5 exposures × 9 outcomes)] was Bonferroni significant for the main MR analysis. Results: Genetically increased TSH was associated with decreased VWF [ß(SE) = -0.020(0.006), p = 0.001] and with decreased fibrinogen [ß(SE) = -0.008(0.002), p = 0.001]. Genetically increased fT4 was associated with increased VWF [ß(SE) = 0.028(0.011), p = 0.012]. Genetically predicted hyperthyroidism was associated with increased VWF [ß(SE) = 0.012(0.004), p = 0.006] and increased FVIII [ß(SE) = 0.013(0.005), p = 0.007]. Genetically predicted hypothyroidism and hyperthyroidism were associated with decreased TPA [ß(SE) = -0.009(0.024), p = 0.024] and increased TPA [ß(SE) = 0.022(0.008), p = 0.008], respectively. MR sensitivity analyses showed similar direction but lower precision. Other coagulation and fibrinolytic factors were inconclusive. Conclusions: In the largest genetic studies currently available, genetically increased TSH and fT4 may be associated with decreased and increased synthesis of VWF, respectively. Since Bonferroni correction may be too conservative given the correlation between the analyzed traits, we cannot reject nominal associations of thyroid traits with coagulation or fibrinolytic factors.