GWAS of thyroid stimulating hormone highlights pleiotropic effects and inverse association with thyroid cancer

No causal associations of genetically predicted birth weight and life course BMI with thyroid function and diseases

OBJECTIVE

Observational studies have suggested associations of birth weight, childhood BMI, and adulthood BMI with thyroid function or diseases. However, the causal relationships remain unclear due to residual confounding inherent in conventional epidemiological studies.

METHODS

We performed a two-sample Mendelian randomization (MR) study to investigate causal relationships of genetically predicted birth weight, childhood BMI, and adulthood BMI with a range of clinically relevant thyroid outcomes. Additionally, we conducted a reverse MR analysis on adulthood BMI. Data on exposures and outcomes were obtained from large-scale genome-wide association study meta-analyses predominantly composed of individuals of European ancestry.

RESULTS

The MR analysis revealed no evidence of causal associations of birth weight or BMI at different life stages with thyrotropin (TSH) levels, hypothyroidism, hyperthyroidism, autoimmune thyroid disorders, or thyroid cancer. Contrarily, thyroid cancer demonstrated a significant causal relationship with increased adulthood BMI (β = 0.010, 95% CI: 0.006-0.015; p = 5.21 × 10-6).

CONCLUSIONS

Our comprehensive MR did not find causal links of birth weight, childhood BMI, or adulthood BMI with thyroid diseases but provided evidence that thyroid cancer may play a role in weight gain. Our research findings offer valuable insights into the intricate relationship between body weight and thyroid health throughout an individual's life.

Periodontitis and thyroid function: A bidirectional Mendelian randomization study

BACKGROUND AND OBJECTIVE

Previous studies suggest interaction between periodontitis and thyroid function, while the causality has not yet been established. We applied the Mendelian randomization (MR) method to assess bidirectional causal association between periodontitis and thyroid-related traits, including free thyroxine (FT4), thyroid stimulating hormone (TSH), hypothyroidism, hyperthyroidism and autoimmune thyroid disease (AITD).

METHODS

Genetic instruments were extracted from large-scale genome-wide association studies on normal-range FT4 (N = 49 269) and TSH (N = 54 288) levels, TSH in full range (N = 119 715); hypothyroidism (discovery/replication cohorts: N = 53 423/334 316), hyperthyroidism (discovery/replication cohorts: N = 51 823/257 552), AITD (N = 755 406) and periodontitis (N = 45 563). Here, the inverse variance weighted (IVW) method was applied as the primary analysis, and robustness of results were assessed by several pleiotropic-robust methods. Results were adjusted for Bonferroni correction thresholds with significant p < .004 (0.05/13) and suggestive p between .004 and .05.

RESULTS

The IVW analysis revealed a suggestively causal linkage between genetic predisposition to periodontitis and the increased risk of hypothyroidism (discovery cohort: odds ratio [OR] = 1.24, 95% confidence interval [CI] = 1.05-1.46, p = .012; replication cohort: OR = 1.06, 95% CI = 1.01-1.11, p = .011). No evidence was found for supporting the impact of periodontitis on hyperthyroidism and AITD risks (associated p ≥ .209), as well as thyroid-related traits on periodontitis risk (associated p ≥ .105). These findings were robust and consistent through sensitivity analysis with other MR models.

CONCLUSION

This bidirectional MR reveals periodontitis should not be attributed to variations in thyroid function but it has potential causal effect on hypothyroidism risk, which provides a better understanding of the relationship between periodontitis and thyroid function, and potential evidence for the clinical intervention of hypothyroidism. Further investigations are warranted to elucidate the nature and underlying mechanisms of this finding.

A polygenic risk score for alcohol-associated cirrhosis among heavy drinkers with European ancestry

BACKGROUND

Polygenic Risk Scores (PRS) based on results from genome-wide association studies offer the prospect of risk stratification for many common and complex diseases. We developed a PRS for alcohol-associated cirrhosis by comparing single-nucleotide polymorphisms among patients with alcohol-associated cirrhosis (ALC) versus drinkers who did not have evidence of liver fibrosis/cirrhosis.

METHODS

Using a data-driven approach, a PRS for ALC was generated using a meta-genome-wide association study of ALC (N=4305) and an independent cohort of heavy drinkers with ALC and without significant liver disease (N=3037). It was validated in 2 additional independent cohorts from the UK Biobank with diagnosed ALC (N=467) and high-risk drinking controls (N=8981) and participants in the Indiana Biobank Liver cohort with alcohol-associated liver disease (N=121) and controls without liver disease (N=3239).

RESULTS

A 20-single-nucleotide polymorphisms PRS for ALC (PRSALC) was generated that stratified risk for ALC comparing the top and bottom deciles of PRS in the 2 validation cohorts (ORs: 2.83 [95% CI: 1.82 -4.39] in UK Biobank; 4.40 [1.56 -12.44] in Indiana Biobank Liver cohort). Furthermore, PRSALC improved the prediction of ALC risk when added to the models of clinically known predictors of ALC risk. It also stratified the risk for metabolic dysfunction -associated steatotic liver disease -cirrhosis (3.94 [2.23 -6.95]) in the Indiana Biobank Liver cohort -based exploratory analysis.

CONCLUSIONS

PRSALC incorporates 20 single-nucleotide polymorphisms, predicts increased risk for ALC, and improves risk stratification for ALC compared with the models that only include clinical risk factors. This new score has the potential for early detection of heavy drinking patients who are at high risk for ALC.

Functional variants in a TTTG microsatellite on 15q26.1 cause familial nonautoimmune thyroid abnormalities

Insufficient thyroid hormone production in newborns is referred to as congenital hypothyroidism. Multinodular goiter (MNG), characterized by an enlarged thyroid gland with multiple nodules, is usually seen in adults and is recognized as a separate disorder from congenital hypothyroidism. Here we performed a linkage analysis of a family with both nongoitrous congenital hypothyroidism and MNG and identified a signal at 15q26.1. Follow-up analyses with whole-genome sequencing and genetic screening in congenital hypothyroidism and MNG cohorts showed that changes in a noncoding TTTG microsatellite on 15q26.1 were frequently observed in congenital hypothyroidism (137 in 989) and MNG (3 in 33) compared with controls (3 in 38,722). Characterization of the noncoding variants with epigenomic data and in vitro experiments suggested that the microsatellite is located in a thyroid-specific transcriptional repressor, and its activity is disrupted by the variants. Collectively, we presented genetic evidence linking nongoitrous congenital hypothyroidism and MNG, providing unique insights into thyroid abnormalities.

Multidimensional data analysis revealed thyroiditis-associated TCF19 SNP rs2073724 as a highly ranked protective variant in thyroid cancer

BACKGROUND

Thyroid cancer represents the most prevalent malignant endocrine tumour, with rising incidence worldwide and high mortality rates among patients exhibiting dedifferentiation and metastasis. Effective biomarkers and therapeutic interventions are warranted in aggressive thyroid malignancies. The transcription factor 19 (TCF19) gene has been implicated in conferring a malignant phenotype in cancers. However, its contribution to thyroid neoplasms remains unclear.

RESULTS

In this study, we performed genome-wide and phenome-wide association studies to identify a potential causal relationship between TCF19 and thyroid cancer. Our analyses revealed significant associations between TCF19 and various autoimmune diseases and human cancers, including cervical cancer and autoimmune thyroiditis, with a particularly robust signal for the deleterious missense variation rs2073724 that is associated with thyroid function, hypothyroidism, and autoimmunity. Furthermore, functional assays and transcriptional profiling in thyroid cancer cells demonstrated that TCF19 regulates important biological processes, especially inflammatory and immune responses. We demonstrated that TCF19 could promote the progression of thyroid cancer in vitro and in vivo and the C>T variant of rs2073724 disrupted TCF19 protein binding to target gene promoters and their expression, thus reversing the effect of TCF19 protein.

CONCLUSIONS

Taken together, these findings implicate TCF19 as a promising therapeutic target in aggressive thyroid malignancies and designate rs2073724 as a causal biomarker warranting further investigation in thyroid cancer.

Investigating the Bidirectional Association of Rheumatoid Arthritis and Thyroid Function: A Methodologic Assessment of Mendelian Randomization

OBJECTIVE

Rheumatoid arthritis (RA) and thyroid dysfunction are frequently observed in the same patient. However, whether they co-occur or exhibit a causal relationship remains uncertain. We aimed to systematically investigate the causal relationship between RA and thyroid function using a large sample and advanced methods.

METHODS

Bidirectional two-sample Mendelian randomization (MR) analysis was performed based on RA and six thyroid function trait data sets from the European population. The robustness of the results was demonstrated using multiple MR methods and a series of sensitivity analyses. Multivariable MR using Bayesian model averaging (MR-BMA) was performed to adjust for possible competing risk factors. A sensitivity data set, which included data from patients with seropositive RA and controls, was used to repeat the analyses. Furthermore, enrichment analysis was employed to discover the underlying mechanism between RA and thyroid functions.

RESULTS

A significantly positive causal effect was identified for RA on autoimmune thyroid disease (AITD) as well as for AITD on RA (P < 0.001). Further sensitivity analyses showed consistent causal estimates from a variety of MR methods. After removing the outliers, MR-BMA results showed that RA and AITD were independent risk factors in their bidirectional causality, even in the presence of other competing risk factors (adjusted P < 0.05). Enrichment analysis showed immune cell activation and immune response play crucial roles in them.

CONCLUSION

Our results illustrate the significant bidirectional causal effect of RA and AITD, which holds even in multiple competing risk factors. Clinical screening for thyroid dysfunction in patients with RA deserves further attention, and vice versa.

Impaired thyroid hormone sensitivity increases the risk of papillary thyroid cancer and cervical lymph node metastasis

BACKGROUND

The association of thyroid hormone sensitivity with papillary thyroid carcinoma (PTC) is unclear. This study investigated the relationship between the thyroid hormone sensitivity indices and the risk of PTC and the influence of thyroid hormone sensitivity on the aggressive clinicopathologic features of PTC.

METHODS

This retrospective study recruited 1225 PTC patients and 369 patients with benign nodules undergoing surgery in Zhongshan Hospital in 2020. The thyroid hormone sensitivity indices were thyroid feedback quantile-based index (TFQI), TSH index (TSHI) and thyrotropin thyroxine resistance index (TT4RI). We employed logistic regression models to explore the correlation between the thyroid hormone sensitivity indices and the risk of PTC and its cervical lymph node metastasis (LNM).

RESULTS

PTC patients had significantly higher levels of TSH, TFQI, TSHI and TT4RI compared to the patients with benign nodules, but thyroid hormone levels did not differ significantly between the two groups. Logistic regression analysis revealed that the higher levels of TFQI, TSHI, and TT4RI were associated with an increased risk of PTC after adjustment for multiple risk factors (TFQI: OR = 1.92, 95% CI: 1.39-2.65, P < 0.001; TSHI: OR = 2.33, 95% CI:1.67-3.26, P < 0.001; TT4RI: OR = 2.41, 95% CI:1.73-3.36, P < 0.001). In addition, patients with decreased thyroid hormone sensitivity had a higher risk of cervical LNM in multiple logistic regression analysis (TFQI: OR = 1.38, 95% CI:1.03-1.86, P = 0.03; TSHI: OR = 1.37, 95% CI:1.02-1.84, P = 0.04; TT4RI: OR = 1.41, 95% CI:1.05-1.89, P = 0.02).

CONCLUSION

Impaired sensitivity to thyroid hormone was associated with an increased risk of PTC, and it is also associated with a higher risk of cervical LNM in PTC patients.

Causal relationships between gut microbiota and hypothyroidism: a Mendelian randomization study

BACKGROUND

Previous studies have shown that the gut microbiota plays an important role in the maintenance of thyroid homeostasis. We aimed to evaluate the causal relationships between gut microbiota and hypothyroidism.

METHODS

Summary statistics for 211 gut microbiota taxa were obtained from the largest available genome-wide association study (GWAS) meta-analysis conducted by the MiBioGen consortium. Summary statistics for hypothyroidism were obtained from two distinct sources: the FinnGen consortium R9 release data (40,926 cases and 274,069 controls) and the UK Biobank data (22,687 cases and 440,246 controls), respectively. A two-sample Mendelian randomization (MR) design was employed, and thorough sensitivity analyses were carried out to ensure the reliability of the results.

RESULTS

Based on the FinnGen consortium, we found increased levels of Intestinimonas (OR = 1.09; 95%CI = 1.02-1.16; P = 0.01) and Ruminiclostridium5 (OR = 1.11; 95%CI = 1.02-1.22; P = 0.02) may be associated with a higher risk of hypothyroidism, while increased levels of Butyrivibrio (OR = 0.95; 95%CI = 0.92-0.99; P = 0.02), Eggerthella (OR = 0.93; 95%CI = 0.88-0.98; P = 0.01), Lachnospiraceae UCG008 (OR = 0.92; 95%CI = 0.85-0.99; P = 0.02), Ruminococcaceae UCG011 (OR = 0.95; 95%CI = 0.90-0.99; P = 0.02), and Actinobacteria (OR = 0.88; 95%CI = 0.80-0.97; P = 0.01) may be associated with a lower risk. According to the UK Biobank data, Eggerthella and Ruminiclostridium5 remain causally associated with hypothyroidism. The sensitivity analysis demonstrates consistent results without evidence of heterogeneity or pleiotropy.

CONCLUSION

This study highlights the impact of specific gut microbiota on hypothyroidism. Strategies to change composition of gut microbiota may hold promise as potential interventions.

Causal relationship between hypothyroidism and temporomandibular disorders: evidence from complementary genetic methods

BACKGROUND

The role of thyroid health in temporomandibular disorders (TMDs) has been emphasized in observational studies. However, whether the causation exists is unclear, and controversy remains about which specific disorder, such as hypothyroidism or hyperthyroidism, is destructive in TMDs. This study aims to investigate the overall and specific causal effects of various thyroid conditions on TMDs.

METHODS

Mendelian randomization (MR) studies were performed using genetic instruments for thyrotropin (TSH, N = 119,715), free thyroxine (fT4, N = 49,269), hypothyroidism (N = 410,141), hyperthyroidism (N = 460,499), and TMDs (N = 211,023). We assessed the overall effect of each thyroid factor via inverse-variance weighted (IVW), weighted median, and MR-Egger methods, and performed extensive sensitivity analyses. Additionally, multivariable MR was conducted to evaluate the direct or indirect effects of hypothyroidism on TMDs whilst accounting for TSH, fT4 and hyperthyroidism, and vice versa.

RESULTS

Univariable MR analyses revealed a causal effect of hypothyroidism on an increased risk of TMDs (IVW OR: 1.12, 95% CI: 1.05-1.20, p = 0.001). No significant association between genetically predicted hyperthyroidism, TSH, or fT4 and TMDs. In the multivariable MR analyses, the effects of hypothyroidism on TMDs occurrence remained significant even after adjSusting for TSH, fT4 and hyperthyroidism (multivariable IVW OR: 1.10, 95% CI: 1.03-1.17, p = 0.006). No pleiotropy and heterogeneity were detected in the analyses (p > 0.05).

CONCLUSIONS

Hypothyroidism might causally increase the risk of TMDs through a direct pathway, highlighting the critical role of managing thyroid health in the prevention of TMDs. Clinicians should give heightened attention to patients with hypothyroidism when seeking medical advice for temporomandibular discomfort. However, caution is warranted due to the potential confounders, pleiotropy, and selection bias in the MR study.

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.

The Rotterdam Study. Design update and major findings between 2020 and 2024

The Rotterdam Study is a population-based cohort study, started in 1990 in the district of Ommoord in the city of Rotterdam, the Netherlands, with the aim to describe the prevalence and incidence, unravel the etiology, and identify targets for prediction, prevention or intervention of multifactorial diseases in mid-life and elderly. The study currently includes 17,931 participants (overall response rate 65%), aged 40 years and over, who are examined in-person every 3 to 5 years in a dedicated research facility, and who are followed-up continuously through automated linkage with health care providers, both regionally and nationally. Research within the Rotterdam Study is carried out along two axes. First, research lines are oriented around diseases and clinical conditions, which are reflective of medical specializations. Second, cross-cutting research lines transverse these clinical demarcations allowing for inter- and multidisciplinary research. These research lines generally reflect subdomains within epidemiology. This paper describes recent methodological updates and main findings from each of these research lines. Also, future perspective for coming years highlighted.

Quantitative Detection of Thyroid-Stimulating Hormone in Patient Samples with a Nanomechanical Single-Antibody Spectro-Immunoassay

Functional disorders of the thyroid remain a global challenge and have profound impacts on human health. Serving as the barometer for thyroid function, thyroid-stimulating hormone (TSH) is considered the single most useful test of thyroid function. However, the prevailing TSH immunoassays rely on two types of antibodies in a sandwich format. The requirement of repeated incubation and washing further complicates the issue, making it unable to meet the requirements of the shifting public health landscape that demands rapid, sensitive, and low-cost TSH tests. Herein, a systematic study is performed to investigate the clinical translational potential of a single antibody-based biosensing platform for the TSH test. The biosensing platform leverages Raman spectral variations induced by the interaction between a TSH antigen and a Raman molecule-conjugated TSH antibody. In conjunction with machine learning, it allows TSH concentrations in various patient samples to be predicted with high accuracy and precision, which is robust against substrate-to-substrate, intra-substrate, and day-to-day variations. It is envisioned that the simplicity and generalizability of this single-antibody immunoassay coupled with the demonstrated performance in patient samples pave the way for it to be widely applied in clinical settings for low-cost detection of hormones, other molecular biomarkers, DNA, RNA, and pathogens.

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.

The Influence of Autoimmune Thyroid Diseases on Viral Pneumonia Development, Including COVID-19: A Two-Sample Mendelian Randomization Study

The association between thyroid function and viral pneumonia has undergone extensive examination, yet the presence of a causal link remains uncertain. The objective of this paper was to employ Two-Sample Mendelian Randomization (MR) analysis to investigate the connections between three thyroid diseases and thyroid hormone indicators with viral pneumonia and COVID-19. We obtained summary statistics datasets from seven genome-wide association studies (GWASs). The primary method used for estimating relationships was inverse-variance weighting (IVW). In addition, we employed weighted median, weighted mode, MR-Egger, and MR-PRESSO as supplementary analytical tools. Sensitivity analyses encompassed Cochran's Q test, MR-Egger intercept test, and MR-PRESSO. Our study revealed significant causal relationships between having a genetic predisposition to autoimmune thyroid disease (AITD) and an increased susceptibility to viral pneumonia (odds ratio [OR]: 1.096; 95% confidence interval [CI]: 1.022-1.176). Moreover, it demonstrated a heightened susceptibility and severity of COVID-19 (OR for COVID-19 susceptibility, COVID-19 hospitalization, and COVID-19 critical illness, with 95% CIs of 1.016, 1.001-1.032; 1.058, 1.003-1.116; 1.045, 1.010-1.081). However, no statistically significant associations were found between TSH, FT4, subclinical hypo- or hyperthyroidism, and the risk of viral pneumonia incidence, or the susceptibility and severity of COVID-19 (all p > 0.05). This study establishes a cause-and-effect relationship between AITD and the development of viral pneumonia, as well as the susceptibility and severity of COVID-19.

Mechanisms of thyrotropin receptor-mediated phenotype variability deciphered by gene mutations and M453T-knockin model

The clinical spectrum of thyrotropin receptor-mediated (TSHR-mediated) diseases varies from loss-of-function mutations causing congenital hypothyroidism to constitutively active mutations (CAMs) leading to nonautoimmune hyperthyroidism (NAH). Variation at the TSHR locus has also been associated with altered lipid and bone metabolism and autoimmune thyroid diseases. However, the extrathyroidal roles of TSHR and the mechanisms underlying phenotypic variability among TSHR-mediated diseases remain unclear. Here we identified and characterized TSHR variants and factors involved in phenotypic variability in different patient cohorts, the FinnGen database, and a mouse model. TSHR CAMs were found in all 16 patients with NAH, with 1 CAM in an unexpected location in the extracellular leucine-rich repeat domain (p.S237N) and another in the transmembrane domain (p.I640V) in 2 families with distinct hyperthyroid phenotypes. In addition, screening of the FinnGen database revealed rare functional variants as well as distinct common noncoding TSHR SNPs significantly associated with thyroid phenotypes, but there was no other significant association between TSHR variants and more than 2,000 nonthyroid disease endpoints. Finally, our TSHR M453T-knockin model revealed that the phenotype was dependent on the mutation's signaling properties and was ameliorated by increased iodine intake. In summary, our data show that TSHR-mediated disease risk can be modified by variants at the TSHR locus both inside and outside the coding region as well as by altered TSHR-signaling and dietary iodine, supporting the need for personalized treatment strategies.

A Detailed Histologic and Molecular Assessment of the Diffuse Sclerosing Variant of Papillary Thyroid Carcinoma

Diffuse sclerosing variant papillary thyroid carcinoma (DS-PTC) is characterized clinically by a predilection for children and young adults, bulky neck nodes, and pulmonary metastases. Previous studies have suggested infrequent BRAFV600E mutation but common RET gene rearrangements. Using strict criteria, we studied 43 DS-PTCs (1.9% of unselected PTCs in our unit). Seventy-nine percent harbored pathogenic gene rearrangements involving RET, NTRK3, NTRK1, ALK, or BRAF; with the remainder driven by BRAFV600E mutations. All 10 pediatric cases were all gene rearranged (P = .02). Compared with BRAFV600E-mutated tumors, gene rearrangement was characterized by psammoma bodies involving the entire lobe (P = .038), follicular predominant or mixed follicular architecture (P = .003), pulmonary metastases (24% vs none, P = .04), and absent classical, so-called "BRAF-like" atypia (P = .014). There was no correlation between the presence of gene rearrangement and recurrence-free survival. Features associated with persistent/recurrent disease included pediatric population (P = .030), gene-rearranged tumors (P = .020), microscopic extrathyroidal extension (P = .009), metastases at presentation (P = .007), and stage II disease (P = .015). We conclude that DS-PTC represents 1.9% of papillary thyroid carcinomas and that actionable gene rearrangements are extremely common in DS-PTC. DS-PTC can be divided into 2 distinct molecular subtypes and all BRAFV600E-negative tumors (1.5% of papillary thyroid carcinomas) are driven by potentially actionable oncogenic fusions.

Thyroid Function and Metabolic Syndrome: A Two-Sample Bidirectional Mendelian Randomization Study

CONTEXT

Thyroid function has been associated with metabolic syndrome (MetS) in a number of observational studies but the direction of effects and the exact causal mechanism of this relationship is still unknown.

OBJECTIVE

To examine genetically predicted effects of thyroid function on MetS risk and its components, and vice versa, using large-scale summary genetic association data.

METHODS

We performed a two-sample bidirectional Mendelian randomization (MR) study using summary statistics from the most comprehensive genome-wide association studies (GWAS) of thyroid-stimulating hormone (TSH, n = 119 715), free thyroxine (fT4, n = 49 269), MetS (n = 291 107), and components of MetS: waist circumference (n = 462 166), fasting blood glucose (n = 281 416), hypertension (n = 463 010), triglycerides (TG, n = 441 016) and high-density lipoprotein cholesterol (HDL-C, n = 403 943). We chose the multiplicative random effects inverse variance weighted (IVW) method as the main analysis. Sensitivity analysis included weighted median and mode analysis, as well as MR-Egger and Causal Analysis Using Summary Effect estimates (CAUSE).

RESULTS

Our results suggest that higher fT4 levels lower the risk of developing MetS (OR = 0.96, P = .037). Genetically predicted fT4 was also positively associated with HDL-C (β = 0.02, P = .008), while genetically predicted TSH was positively associated with TG (β = 0.01, P = .044). These effects were consistent across different MR analyses and confirmed with the CAUSE analysis. In the reverse direction MR analysis, genetically predicted HDL-C was negatively associated with TSH (β = -0.03, P = .046) in the main IVW analysis.

CONCLUSION

Our study suggests that variations in normal-range thyroid function are causally associated with the diagnosis of MetS and with lipid profile, while in the reverse direction, HDL-C has a plausible causal effect on reference-range TSH levels.

Association Between Educational Attainment and Thyroid Function: Results From Mendelian Randomization and the NHANES Study

CONTEXT

Many observational studies have reported on the association between educational attainment (EA) and thyroid function, but the causal relationship remains unclear.

OBJECTIVE

We aimed to obtain causal effects of EA on thyroid function and to quantify the mediating effects of modifiable risk factors.

METHODS

Two-sample mendelian randomization (MR) was performed by using summary statistics from large genome-wide association studies (GWAS) to assess the effect of EA on thyroid function, including hypothyroidism, hyperthyroidism, thyrotropin (TSH), and free thyroxine (FT4). A multivariable analysis was conducted to assess the mediating role of smoking and help to explain the association between EA and thyroid function. Similar analysis was further performed using data from the National Health and Nutrition Examination Survey (NHANES) 1999 to 2002.

RESULTS

In MR analysis, EA was causally associated with TSH (β = .046; 95% CI, 0.015-0.077; P = 4.00 × 10-3), rather than hypothyroidism, hyperthyroidism, and FT4. Importantly, smoking could serve as a mediator in the association between EA and TSH, in which the mediating proportion was estimated to be 10.38%. After adjusting for smoking in the multivariable MR analysis, the β value of EA on TSH was attenuated to 0.030 (95% CI, 0.016-0.045; P = 9.32 × 10-3). Multivariable logistic regression model in NHANES suggested a dose-response relationship between TSH (quartile [Q]4 vs Q1: odds ratio = 1.33; 95% CI, 1.05-1.68; P for trend = .023) and EA. Smoking, systolic blood pressure, and body mass index partially mediated the association between EA and TSH, with the proportion of the mediation effects being 43.82%, 12.28%, and 6.81%, respectively.

CONCLUSION

There is a potentially causal association between EA and TSH, which could be mediated by several risk factors, such as smoking.

Role of ST6GAL1 in Thyroid Cancers: Insights from Tissue Analysis and Genomic Datasets

Thyroid cancer is the predominant endocrine-related malignancy. ST6 β-galactoside α2,6-sialyltransferase 1 (ST6GAL1) has been studied in various types of cancers; however, the expression and function of ST6GAL1 in thyroid cancer has not been investigated so far. Previously, we conducted two genome-wide association studies and have identified the association of the ST6GAL1 gene with plasma thyroglobulin (Tg) levels. Since Tg levels are altered in thyroid pathologies, in the current study, we wanted to evaluate the expression of ST6GAL1 in thyroid cancer tissues. We performed an immunohistochemical analysis using human thyroid tissue from 89 patients and analyzed ST6GAL1 protein expression in papillary thyroid cancer (including follicular variant and microcarcinoma) and follicular thyroid cancer in comparison to normal thyroid tissue. Additionally, ST6GAL1 mRNA levels from The Cancer Genome Atlas (TCGA, n = 572) and the Genotype-Tissue Expression (GTEx) project (n = 279) were examined. The immunohistochemical analysis revealed higher ST6GAL1 protein expression in all thyroid tumors compared to normal thyroid tissue. TCGA data revealed increased ST6GAL1 mRNA levels in both primary and metastatic tumors versus controls. Notably, the follicular variant of papillary thyroid cancer exhibited significantly higher ST6GAL1 mRNA levels than classic papillary thyroid cancer. High ST6GAL1 mRNA levels significantly correlated with lymph node metastasis status, clinical stage, and reduced survival rate. ST6GAL1 emerges as a potential cancer-associated glycosyltransferase in thyroid malignancies, offering valuable insights into its diagnostic and prognostic significance.

Genetically predicted thyroid function and risk of colorectal cancer: a bidirectional Mendelian randomization study

BACKGROUND

Observational studies have reported an association between thyroid function and colorectal cancer (CRC), with conflicting results. Elucidating the causal relationship between thyroid function and CRC facilitates the development of new preventive strategies to reduce CRC incidence.

METHOD

We applied a two-sample Mendelian randomization (MR) method to evaluate the causal relationship between five thyroid-related indexes, including hyperthyroidism, hypothyroidism, thyroid stimulating hormone (TSH), free thyroxine (FT4) and basal metabolic rate (BMR), and CRC. Genome-wide association study statistics for thyroid-related phenotypes were obtained from the ThyroidOmics consortium, and summary statistics for genetic associations with CRC were obtained from the FinnGen consortium. We set a series of criteria to screen single nucleotide polymorphisms (SNPs) as instrumental variables and then performed bidirectional MR analysis, stratified analysis and extensive sensitivity analysis. Multiplicative random-effects inverse variance weighted was the primary analysis method, supplemented by weighted median and MR-Egger.

RESULT

We identified 12 SNPs for hyperthyroidism, 10 SNPs for hypothyroidism, 41 SNPs for TSH, 18 SNPs for FT4, and 556 SNPs for BMR. Genetically predicted hyperthyroidism, hypothyroidism, TSH, and FT4 were not associated with CRC risk (all P > 0.05). Sensitivity analysis revealed no heterogeneity or pleiotropy. Genetically predicted BMR was significantly associated with increased CRC risk after removing outlier (OR = 1.30, P = 0.0029). Stratified analysis showed that BMR was significantly associated with colon cancer (OR = 1.33, P = 0.0074) but not rectal cancer. In the reverse analysis, there was no evidence of an effect of CRC on thyroid function (all P > 0.05).

CONCLUSION

Our bidirectional MR analysis provides new insights into the relationship between thyroid function and CRC. CRC prevention may benefit from enhanced screening of high BMR populations.

Association Between Serum Thyroid Measurements and Hyperhomocysteinemia in Euthyroid Subjects: A Retrospective Cross-Sectional Study

Objective

The objective of this study was to investigate the association between serum thyroid measurements and homocysteine (HCy) in euthyroid participants.

Methods

This retrospective study was based on Hospital Information Systems. After excluding participants with thyroid dysfunction and those who had recently taken medications that affected serum HCy, 775 participants were enrolled. We compared the serum thyroid function measurements of patients with or without hyperhomocysteinemia (HHCy) and analyzed the effect of thyroid indicators on HHCy prevalence and HCy levels. Multivariate regression analysis was utilized to analyze the association of thyroid-stimulating hormone (TSH) and thyroid peroxidase (TPOAb) with HCy.

Results

The serum TSH level (2.10 ± 1.06 mIU/L) of HHCy patients (n = 98) was significantly higher than controls (n = 677) (1.65 ± 0.90 mIU/L) (p < 0.05), as was the positive rate of TPOAb (19.4% vs 10.0%, p < 0.05). The serum HCy levels in subjects with TSH within the highest quartile were significantly higher than those in the lowest quartile (13.49 ± 7.78 vs 9.81 ± 3.59 μmol/L, p < 0.05). HCy was also significantly higher in TPOAb-positive patients than in negative subjects (14.06 ± 8.89 vs 11.48 ± 5.47 μmol/L, p < 0.05). Among the TSH quartiles, the prevalence of HHCy showed a similar significant upward trend to that described above. The prevalence of HHCy was also significantly higher in TPOAb-positive patients. The results of multivariate regression analysis suggested that both TSH elevation and TPOAb positivity were independent risk factors for HCy elevation and HHCy prevalence. However, we found no definitive association between linear increases in TPOAb titers and HCy concentrations or HHCy prevalence.

Conclusion

Patients with HHCy had significantly higher TSH levels and positive rates of TPOAb. Elevated TSH and positive TPOAb levels were independent risk factors for elevated HCy concentrations and HHCy risk.

Genome-wide association study of thyroid-stimulating hormone highlights new genes, pathways and associations with thyroid disease

Thyroid hormones play a critical role in regulation of multiple physiological functions and thyroid dysfunction is associated with substantial morbidity. Here, we use electronic health records to undertake a genome-wide association study of thyroid-stimulating hormone (TSH) levels, with a total sample size of 247,107. We identify 158 novel genetic associations, more than doubling the number of known associations with TSH, and implicate 112 putative causal genes, of which 76 are not previously implicated. A polygenic score for TSH is associated with TSH levels in African, South Asian, East Asian, Middle Eastern and admixed American ancestries, and associated with hypothyroidism and other thyroid disease in South Asians. In Europeans, the TSH polygenic score is associated with thyroid disease, including thyroid cancer and age-of-onset of hypothyroidism and hyperthyroidism. We develop pathway-specific genetic risk scores for TSH levels and use these in phenome-wide association studies to identify potential consequences of pathway perturbation. Together, these findings demonstrate the potential utility of genetic associations to inform future therapeutics and risk prediction for thyroid diseases.

Causal link between thyroid function and schizophrenia: a two-sample Mendelian randomization study

Schizophrenia is a chronic psychiatric disorder with inconsistent behavioral and cognitive abnormalities with profound effects on the individual and the society. Individuals with schizophrenia have altered thyroid function, but results from observational studies are conflicting. To date, it remains unclear whether and in which direction there is a causal relationship between thyroid function and schizophrenia. To investigate causal paths, a bidirectional two-sample Mendelian randomization (MR) study was conducted using summary statistics from genome-wide association studies including up to 330,132 Europeans. Thyroid function was described by the normal-range thyroid-stimulating hormone (TSH) and free thyroxine levels as well as an increased and decreased TSH status. The iterative radial inverse-variance weighted approach with modified second order weights was used as the main method. Based on a discovery and replication sample for schizophrenia, pooled effect estimates were derived using a fixed-effect meta-analysis. Robustness of results was assessed using both a range of pleiotropy robust methods and a network analysis that clustered genetic instruments potentially responsible for horizontal pleiotropy. Genetic liability for hypothyroidism was inversely associated with schizophrenia ([Formula: see text]; 95% CI: (-0.10; -0.02); [Formula: see text]). No notable associations were observed between other thyroid parameters and schizophrenia. Furthermore, no associations could be detected in the reverse direction. Our results suggest that an elevated level of TSH reduce the risk for schizophrenia. The role of thyroid function and the hypothalamic-pituitary-thyroid axis in the development of schizophrenia should be subject of further research.

Whole Exome-Wide Association Identifies Rare Variants in GALNT9 Associated with Middle Eastern Papillary Thyroid Carcinoma Risk

Papillary thyroid carcinoma (PTC) is the commonest thyroid cancer. The majority of inherited causes of PTC remain elusive. However, understanding the genetic underpinnings and origins remains a challenging endeavor. An exome-wide association study was performed to identify rare germline variants in coding regions associated with PTC risk in the Middle Eastern population. By analyzing exome-sequencing data from 249 PTC patients (cases) and 1395 individuals without any known cancer (controls), GALNT9 emerged as being strongly associated with rare inactivating variants (RIVs) (4/249 cases vs. 1/1395 controls, OR = 22.75, p = 5.09 × 10-5). Furthermore, three genes, TRIM40, ARHGAP23, and SOX4, were enriched for rare damaging variants (RDVs) at the exome-wide threshold (p < 2.5 × 10-6). An additional seven genes (VARS1, ZBED9, PRRC2A, VWA7, TRIM31, TRIM40, and COL8A2) were associated with a Middle Eastern PTC risk based on the sequence kernel association test (SKAT). This study underscores the potential of GALNT9 and other implicated genes in PTC predisposition, illuminating the need for large collaborations and innovative approaches to understand the genetic heterogeneity of PTC predisposition.

Genetic determinants of thyroid function in children

OBJECTIVE

Genome-wide association studies in adults have identified 42 loci associated with thyroid stimulating hormone (TSH) and 21 loci associated with free thyroxine (FT4) concentrations. While biologically plausible, age-dependent effects have not been assessed. We aimed to study the association of previously identified genetic determinants of TSH and FT4 with TSH and FT4 concentrations in newborns and (pre)school children.

METHODS

We selected participants from three population-based prospective cohorts with data on genetic variants and thyroid function: Generation R (N = 2169 children, mean age 6 years; N = 2388 neonates, the Netherlands), the Avon Longitudinal Study of Parents and Children (ALSPAC; N = 3382, age 7.5 years, United Kingdom), and the Brisbane Longitudinal Twin Study (BLTS; N = 1680, age 12.1 years, Australia). The association of single nucleotide polymorphisms (SNPs) with TSH and FT4 concentrations was studied with multivariable linear regression models. Weighted polygenic risk scores (PRSs) were defined to combine SNP effects.

RESULTS

In childhood, 30/60 SNPs were associated with TSH and 11/31 SNPs with FT4 after multiple testing correction. The effect sizes for AADAT, GLIS3, TM4SF4, and VEGFA were notably larger than in adults. The TSH PRS explained 5.3%-8.4% of the variability in TSH concentrations; the FT4 PRS explained 1.5%-4.2% of the variability in FT4 concentrations. Five TSH SNPs and no FT4 SNPs were associated with thyroid function in neonates.

CONCLUSIONS

The effects of many known thyroid function SNPs are already apparent in childhood and some might be notably larger in children as compared to adults. These findings provide new knowledge about genetic regulation of thyroid function in early life.

Mendelian randomization study of thyroid function and anti-Müllerian hormone levels

Objective

Although previous studies have reported an association between thyroid function and anti-Müllerian hormone (AMH) levels, which is considered a reliable marker of ovarian reserve, the causal relationship between them remains uncertain. This study aims to investigate whether thyrotropin (TSH), free thyroxine (fT4), hypo- and hyperthyroidism are causally linked to AMH levels.

Methods

We obtained summary statistics from three sources: the ThyroidOmics Consortium (N = 54,288), HUNT + MGI + ThyroidOmics meta-analysis (N = 119,715), and the most recent AMH genome-wide association meta-analysis (N = 7,049). Two-sample MR analyses were conducted using instrumental variables representing TSH and fT4 levels within the normal range. Additionally, we conducted secondary analyses to explore the effects of hypo- and hyperthyroidism. Subgroup analyses for TSH were also performed.

Results

MR analyses did not show any causality relationship between thyroid function and AMH levels, using normal range TSH, normal range fT4, subclinical hypothyroidism, subclinical hyperthyroidism and overt hypothyroidism as exposure, respectively. In addition, neither full range TSH nor TSH with individuals <50 years old was causally associated with AMH levels. MR sensitivity analyses guaranteed the robustness of all MR results, except for the association between fT4 and AMH in the no-DIO1+DIO2 group.

Conclusion

Our findings suggest that there was no causal association between genetically predicted thyroid function and AMH levels in the European population.

Genome-wide association studies: utility and limitations for research in physiology

Physiological systems are subject to interindividual variation encoded by genetics. Genome-wide association studies (GWAS) operate by surveying thousands of genetic variants from a substantial number of individuals and assessing their association to a trait of interest, be it a physiological variable, a molecular phenotype (e.g. gene expression), or even a disease or condition. Through a myriad of methods, GWAS downstream analyses then explore the functional consequences of each variant and attempt to ascertain a causal relationship to the phenotype of interest, as well as to delve into its links to other traits. This type of investigation allows mechanistic insights into physiological functions, pathological disturbances and shared biological processes between traits (i.e. pleiotropy). An exciting example is the discovery of a new thyroid hormone transporter (SLC17A4) and hormone metabolising enzyme (AADAT) from a GWAS on free thyroxine levels. Therefore, GWAS have substantially contributed with insights into physiology and have been shown to be useful in unveiling the genetic control underlying complex traits and pathological conditions; they will continue to do so with global collaborations and advances in genotyping technology. Finally, the increasing number of trans-ancestry GWAS and initiatives to include ancestry diversity in genomics will boost the power for discoveries, making them also applicable to non-European populations.

Advances in Antioxidant Applications for Combating 131I Side Effects in Thyroid Cancer Treatment

Thyroid cancer is the most common endocrine cancer, and its prevalence has been increasing for decades. Approx. 95% of differentiated thyroid carcinomas are treated using ¹³¹iodine (¹³¹I), a radionuclide with a half-life of 8 days, to achieve optimal thyroid residual ablation following thyroidectomy. However, while ¹³¹I is highly enriched in eliminating thyroid tissue, it can also retain and damage other body parts (salivary glands, liver, etc.) without selectivity, and even trigger salivary gland dysfunction, secondary cancer, and other side effects. A significant amount of data suggests that the primary mechanism for these side effects is the excessive production of reactive oxygen species, causing a severe imbalance of oxidant/antioxidant in the cellular components, resulting in secondary DNA damage and abnormal vascular permeability. Antioxidants are substances that are capable of binding free radicals and reducing or preventing the oxidation of the substrate in a significant way. These compounds can help prevent damage caused by free radicals, which can attack lipids, protein amino acids, polyunsaturated fatty acids, and double bonds of DNA bases. Based on this, the rational utilization of the free radical scavenging function of antioxidants to maximize a reduction in ¹³¹I side effects is a promising medical strategy. This review provides an overview of the side effects of ¹³¹I, the mechanisms by which ¹³¹I causes oxidative stress-mediated damage, and the potential of natural and synthetic antioxidants in ameliorating the side effects of ¹³¹I. Finally, the disadvantages of the clinical application of antioxidants and their improving strategies are prospected. Clinicians and nursing staff can use this information to alleviate ¹³¹I side effects in the future, both effectively and reasonably.

Association Between Systemic Lupus Erythematosus and Primary Hypothyroidism: Evidence from Complementary Genetic Methods

INTRODUCTION

Systemic lupus erythematosus (SLE) and hypothyroidism often coexist in observational studies; however, the causal relationship between them remains controversial.

METHODS

Complementary genetic approaches, including genetic correlation, Mendelian randomization (MR), and colocalization analysis, were conducted to assess the potential causal association between SLE and primary hypothyroidism using summary statistics from large-scale genome-wide association studies. The association between SLE and thyroid-stimulating hormone (TSH) was further analyzed to help interpret the findings. In addition, findings were verified using a validation data set, as well as through different MR methods with different model assumptions.

RESULTS

The linkage disequilibrium score regression revealed a shared genetic structure between SLE and primary hypothyroidism, with the significant genetic correlation estimated to be 0.2488 (P = 6.00 × 10-4). MR analysis with the inverse variance weighted method demonstrated a bidirectional causal relationship between SLE and primary hypothyroidism. The odds ratio (OR) of SLE on primary hypothyroidism was 1.037 (95% CI, 1.013-1.061; P = 2.00 × 10-3) and that of primary hypothyroidism on SLE was 1.359 (95% CI, 1.217-1.520; P < 0.001). The OR of SLE on TSH was 1.007 (95% CI, 1.001-1.013; P = 0.032). However, TSH was not causally associated with SLE (P = 0.152). Similar results were found using different MR methods. In addition, colocalization analysis suggested that shared causal variants existed between SLE and primary hypothyroidism. The results of the validation analysis indicated a bidirectional causal relationship between SLE and primary hypothyroidism, as well as shared loci.

CONCLUSION

In summary, a bidirectional causal relationship between SLE and primary hypothyroidism was observed with complementary genetic approaches.

Migraine and thyroid dysfunction: Co-occurrence, shared genes and biological mechanisms

BACKGROUND AND PURPOSE

Migraine and thyroid dysfunction, particularly hypothyroidism, are common medical conditions and are known to have high heritability. Thyroid function measures, thyroid stimulating hormone (TSH) and free thyroxine (fT4), are also known to be genetically influenced. Although observational epidemiological studies report an increased co-occurrence of migraine and thyroid dysfunction, a clear and combined interpretation of the findings is currently lacking. A narrative review is provided of the epidemiological and genetic association evidence linking migraine, hypothyroidism, hyperthyroidism and thyroid hormones TSH and fT4.

METHODS

An extensive literature search was conducted in the PubMed database for epidemiological, candidate gene and genome-wide association studies using the terms migraine, headache, thyroid hormones, TSH, fT4, thyroid function, hypothyroidism and hyperthyroidism.

RESULTS

Epidemiological studies suggest a bidirectional relationship between migraine and thyroid dysfunction. However, the nature of the relationship remains unclear, with some studies suggesting migraine increases the risk for thyroid dysfunction whilst other studies suggest the reverse. Early candidate gene studies have provided nominal evidence for MTHFR and APOE, whilst more recently genome-wide association studies have provided robust evidence for THADA and ITPK1 being associated with both migraine and thyroid dysfunction.

CONCLUSIONS

These genetic associations improve our understanding of the genetic relationship between migraine and thyroid dysfunction, provide an opportunity to develop biomarkers to identify migraine patients most likely to benefit from thyroid hormone therapy, and indicate that further cross-trait genetic studies have excellent potential to provide biological insight into their relationship and inform clinical interventions.

Variants in TPO rs2048722, PTCSC2 rs925489 and SEMA4G rs4919510 affect thyroid carcinoma susceptibility risk

BACKGROUND

Thyroid carcinoma (THCA) is a malignant endocrine tumor all around the world, which is influenced by genetic and environmental factors.

OBJECTIVE

To explore the association between TPO rs2048722, PTCSC2 rs925489, SEMA4G rs4919510 polymorphisms and THCA susceptibility in Chinese population.

METHODS

We recruited 365 THCA patients and 498 normal controls for the study. Logistic regression analysis was used to evaluate the association between TPO rs2048722, PTCSC2 rs925489, SEMA4G rs4919510 polymorphisms and THCA susceptibility. MDR was used to assess the genetic interactions among the three SNPs.

RESULTS

Overall analysis demonstrated that rs925489 of PTCSC2 was evidently associated with increased risk of THCA in multiple genetic models (OR = 1.59, 95%CI = 1.12-2.24, p = 0.009). The results of stratified analysis illustrated that rs2048722 of TPO can significantly increase the THCA susceptibility of participants less than or equal to 44 years old and smokers. Similarly, rs925489 of PTCSC2 obviously improved the risk of THCA among participants older than 44 years, males, smokers and drinkers. However, rs4919510 of SEMA4G has a protective effect on the development of THCA among participants with less than or equal to 44 years old and non-drinkers. Interestingly, there was a strong genetic interaction among the three SNPs in the occurrence of THCA risk.

CONCLUSION

TPO rs2048722, PTCSC2 rs925489 and SEMA4G rs4919510 polymorphisms were evidently associated with the risk of THCA in the Chinese population, which was affected by age, gender, smoking and drinking consumption.

Shared genetics and causal relationships between migraine and thyroid function traits

BACKGROUND

Epidemiological studies have reported a comorbid relationship between migraine and thyroid dysfunction.

METHODS

We investigated the genetic relationship between migraine and thyroid function traits using genome-wide association study (GWAS) data.

RESULTS

We found a significant genetic correlation (r_g) with migraine for hypothyroidism (r_g = 0.0608), secondary hypothyroidism (r_g = 0.195), free thyroxine (fT4) (r_g = 0.0772), and hyperthyroidism (r_g = -0.1046), but not thyroid stimulating hormone (TSH). Pairwise GWAS analysis revealed two shared loci with TSH and 11 shared loci with fT4. Cross-trait GWAS meta-analysis of migraine identified novel genome-wide significant loci: 17 with hypothyroidism, one with hyperthyroidism, five with secondary hypothyroidism, eight with TSH, and 15 with fT4. Of the genes at these loci, six (RERE, TGFB2, APLF, SLC9B1, SGTB, BTBD16; migraine + hypothyroidism), three (GADD45A, PFDN1, RSPH6A; migraine + TSH), and three (SSBP3, BRD3, TEF; migraine + fT4) were significant in our gene-based analysis (p_{Fisher's combined P-value} < 2.04 × 10^-6). In addition, causal analyses suggested a negative causal relationship between migraine and hyperthyroidism (p = 8.90 × 10^-3) and a positive causal relationship between migraine and secondary hypothyroidism (p = 1.30 × 10^-3).

CONCLUSION

These findings provide strong evidence for genetic correlation and suggest complex causal relationships between migraine and thyroid traits.

Modification of TSH-related genetic effects by indicators of socioeconomic position

OBJECTIVE

Thyroid-stimulating hormone (TSH) is influenced by genetic and environmental factors such as socioeconomic position (SEP). However, interactions between TSH-related genetic factors and indicators of SEP have not been investigated to date. The aim of the study was to determine whether education and income as SEP indicators may interact with TSH-related genetic effect allele sum scores (GESTSH_2013 and GESTSH_2020) based on two different GWAS meta-analyses that affect TSH values in a population-based study.

METHODS

In 4085 participants of the Heinz Nixdorf Recall Study associations between SEP indicators, GESTSH and TSH were quantified using sex- and age-adjusted linear regression models. Interactions between SEP indicators and GESTSH were assessed by GESTSH × SEP interaction terms, single reference joint effects and calculating genetic effects stratified by SEP group.

RESULTS

Participants within the highest education group showed the strongest genetic effect with on average 1.109-fold (95% CI: 1.067-1.155) higher TSH values per GESTSH_2013 SD, while in the lowest education group, the genetic effect was less strong (1.061-fold (95% CI: 1.022-1.103)). In linear regression models including interaction terms, some weak indication for a positive GESTSH_2013 by education interaction was observed showing an interaction effect size estimate of 1.005 (95% CI: 1.000-1.010) per year of education and GESTSH_2013 SD. No indication for interaction was observed for using income as SEP indicator. Using the GESTSH_2020, similar results were observed.

CONCLUSION

Our results gave some indication that education may affect the expression of TSH-related genetic effects. Stronger genetic effects in high-education groups may be explained by environmental factors that have an impact on gene expression and are more prevalent in high SEP groups.

Genome-wide association studies for thyroid physiology and diseases

Genetic factors are involved in the etiology of most diseases, but prior to 2000, the methods for identifying such factors were very limited. Genome-wide association study (GWAS), developed in the 2000s, is an analytical method that can be applied to most diseases, including endocrine disorders. GWAS has provided a wealth of information on disease risks and the molecular pathogenesis of many human diseases. This review summarizes key findings from GWAS for thyroid physiology and diseases, and illustrates how GWAS is a powerful research tool to elucidate the molecular mechanisms of the diseases.

Low thyroid function is associated with an increased risk of advanced fibrosis in patients with metabolic dysfunction-associated fatty liver disease

AIMS

Observational studies showed that low thyroid function may perturb liver function. We aimed to evaluate the association of low thyroid function with both metabolic dysfunction-associated fatty liver disease (MAFLD) and advanced hepatic fibrosis.

METHODS

Participants who underwent abdominal ultrasonography and thyroid function test in a Chinese hospital from 2015 to 2021were enrolled. Fibrosis-4 index (FIB-4) > 2.67 and/or non-alcoholic fatty liver disease fibrosis score (NFS) > 0.676 were used to define advanced fibrosis. Descriptive analyses were performed to characterize the epidemiology of MAFLD according to levels of thyroid-stimulating hormone (TSH). The logistic regression model was applied to estimate the association of low thyroid function with MAFLD and advanced fibrosis.

RESULTS

A total of 19,946 participants (52.78% males, mean age: 47.31 years, 27.55% MAFLD) were included, among which 14,789 were strict-normal thyroid function, 4,328 were low-normal thyroid function, 829 were subclinical hypothyroidism. TSH levels were significantly higher in MAFLD patients with a FIB-4 > 2.67 and /or NFS > 0.676 than their counterparts. The logistic regression model adjusted for age and sex showed that low-normal thyroid function increased the risk of MAFLD (odds ratio [OR] = 1.09; 95% confidence interval [CI] 1.01-1.18). Multivariable regression model adjusted for age, sex, body mass index, type 2 diabetes, and hypertension showed low-normal thyroid function increased the risk of advanced fibrosis in patients with MAFLD (FIB-4 > 2.67: OR = 1.41, 95% CI 1.02-1.93; NFS > 0.676: OR = 1.72, 95% CI 1.08-2.72).

CONCLUSION

Elevated TSH concentrations are associated with advanced hepatic fibrosis, even in the euthyroid state.

Cross-Trait Genetic Analyses Indicate Pleiotropy and Complex Causal Relationships between Headache and Thyroid Function Traits

Epidemiological studies have reported a comorbid relationship between headache and thyroid traits; however, little is known about the shared genetics and causality that contributes to this association. We investigated the genetic overlap and associations between headache and thyroid function traits using genome-wide association study (GWAS) data. We found a significant genetic correlation (rg) with headache and hypothyroidism (rg = 0.09, p = 2.00 × 10−4), free thyroxine (fT4) (rg = 0.08, p = 5.50 × 10−3), and hyperthyroidism (rg = −0.14, p = 1.80 × 10−3), a near significant genetic correlation with secondary hypothyroidism (rg = 0.20, p = 5.24 × 10−2), but not with thyroid stimulating hormone (TSH). Pairwise-GWAS analysis revealed six, 14, four and five shared (pleiotropic) loci with headache and hypothyroidism, hyperthyroidism, secondary hypothyroidism, and fT4, respectively. Cross-trait GWAS meta-analysis identified novel genome-wide significant loci for headache: five with hypothyroidism, three with secondary hypothyroidism, 12 with TSH, and nine with fT4. Of the genes at these loci, six (FAF1, TMX2-CTNND1, AARSD1, PLCD3, ZNF652, and C20orf203; headache-TSH) and six (HMGB1P45, RPL30P1, ZNF462, TMX2-CTNND1, ITPK1, SECISBP2L; headache-fT4) were significant in our gene-based analysis (pFisher’s combined p-value < 2.09 × 10−6). Our causal analysis suggested a positive causal relationship between headache and secondary hypothyroidism (p = 3.64 × 10−4). The results also suggest a positive causal relationship between hypothyroidism and headache (p = 2.45 × 10−3) and a negative causal relationship between hyperthyroidism and headache (p = 1.16 × 10−13). These findings suggest a strong evidence base for a genetic correlation and complex causal relationships between headache and thyroid traits.

Genetically predicted selenium concentrations and thyroid function: A two-sample Mendelian randomization study

OBJECTIVE

The impact of selenium (Se) on human thyroid function remains unclear, with inconsistent results from recent epidemiological studies. Moreover, the observed associations are prone to bias due to potential confounding and reverse causation. Mendelian randomization (MR) analysis facilitates the large minimization of biases produced by environmental and lifestyle influences, providing unconfounded estimates of causal effects using instrumental variables. We aim to examine the association between Se concentrations and human thyroid function using a two-sample MR analysis.

DESIGN AND METHODS

Genetic instruments for Se concentrations, including toenail and blood (TAB) and blood Se concentrations, were identified from a genome-wide association study (GWAS) of blood Se (n = 5477) and toenail Se levels (n = 4162). GWAS summary statistics on thyroid phenotypes were downloaded from the ThyroidOmics consortium, including thyroid-stimulating hormone (TSH) (n = 54,288), free thyroxin (FT4) (n = 49,269), hypo (n = 53,423), and hyperthyroidism (n = 51,823). The MR study was conducted using the inverse-variance weighted (IVW) method, supplemented with the weighted median and the mode-based method.

RESULTS

Genetically determined TAB Se was negatively associated with FT4 (β = -.067; 95% confidence interval [CI] = -0.106, -0.028; p = 0.001) using the IVW analyses, as well in the additional analyses using the weighted median and weighted-mode methods. No evidence in heterogeneity, pleiotropy or outlier single-nucleotide polymorphisms was detected (all p > 0.05). Suggestive casual association between increased genetically determined TAB Se concentrations and decreased hypothyroidism risk was found by the IVW method (odds ratio [OR] = 0.847; 95% CI = 0.728, 0.985; p = 0.031). The causal effect of TAB Se on FT4 was observed in women (β = -.076; 95% CI = -0.129, -0.024; p = 0.004). However, the influence of genetically determined higher Se concentrations on TSH levels and hyperthyroidism revealed insignificance in the primary and sensitivity analyses.

CONCLUSIONS

The present MR study indicated that high Se concentration enable the decreasing of FT4 levels, and the effects of Se concentrations on FT4 remain sex-specific.

Meta-analysis fine-mapping is often miscalibrated at single-variant resolution

Meta-analysis is pervasively used to combine multiple genome-wide association studies (GWASs). Fine-mapping of meta-analysis studies is typically performed as in a single-cohort study. Here, we first demonstrate that heterogeneity (e.g., of sample size, phenotyping, imputation) hurts calibration of meta-analysis fine-mapping. We propose a summary statistics-based quality-control (QC) method, suspicious loci analysis of meta-analysis summary statistics (SLALOM), that identifies suspicious loci for meta-analysis fine-mapping by detecting outliers in association statistics. We validate SLALOM in simulations and the GWAS Catalog. Applying SLALOM to 14 meta-analyses from the Global Biobank Meta-analysis Initiative (GBMI), we find that 67% of loci show suspicious patterns that call into question fine-mapping accuracy. These predicted suspicious loci are significantly depleted for having nonsynonymous variants as lead variant (2.7×; Fisher's exact p = 7.3 × 10^-4). We find limited evidence of fine-mapping improvement in the GBMI meta-analyses compared with individual biobanks. We urge extreme caution when interpreting fine-mapping results from meta-analysis of heterogeneous cohorts.

COVID-19 and Thyroid Function: A Bi-Directional Two-Sample Mendelian Randomization Study

Background: Thyroid dysfunction has been observed among some patients with coronavirus disease (COVID-19). It is unclear whether severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection (or its severity) leads to the development of thyroid dysfunction, or vice versa. In this study, we examined the bi-directional causal relationship between host genetic liability to three COVID-19 phenotypes (including SARS-CoV-2 infection, hospitalized and severe COVID-19) and three thyroid dysfunction traits (including hyperthyroidism, hypothyroidism, and autoimmune thyroid disease [AITD]) and three continuous traits of thyroid hormones (including thyrotropin [TSH] and free thyroxine [fT4] within reference range, and TSH in full range). Methods: Summary statistics from the largest available meta-analyses of human genome-wide association studies were retrieved for the following variables: SARS-CoV-2 infection (n = 1,348,701), COVID-19 hospitalization (n = 1,557,411), severe COVID-19 (n = 1,059,456), hyperthyroidism (n = 51,823), hypothyroidism (n = 53,423), AITD (n = 755,406), TSH within reference range (n = 54,288), fT4 within reference range (n = 49,269), and TSH in full range (n = 119,715). Using a two-sample Mendelian randomization (MR) approach, the inverse-variance weighted (IVW) method was adopted as the main MR analysis. Weighted median, contamination mixture, MR-Egger, and MR pleiotropy residual sum and outlier (MR-PRESSO) methods were applied as sensitivity analyses. Results: Host genetic susceptibility to SARS-CoV-2 infection was causally associated with hypothyroidism in the main IVW analysis (per doubling in prevalence of SARS-CoV-2 infection, odds ratio [OR] = 1.335; 95% confidence interval [CI]: 1.167-1.526; p = 2.4 × 10^-5, surpassing the Bonferroni multiple-testing threshold). Similar causal estimates were observed in the sensitivity analyses (weighted median: OR = 1.296; CI: 1.066-1.575; p = 9 × 10^-3; contamination mixture: OR = 1.356; CI: 1.095-1.818; p = 0.013; MR-Egger: OR = 1.712; CI: 1.202-2.439; p = 2.92 × 10^-3, and MR-PRESSO: OR = 1.335; CI: 1.156-1.542; p = 5.73 × 10^-4). Host genetic liability to hospitalized or severe COVID-19 was not associated with thyroid dysfunction or thyroid hormone levels. In the reverse direction, there was no evidence to suggest that genetic predisposition to thyroid dysfunction or genetically determined thyroid hormone levels altered the risk of the COVID-19 outcomes. Conclusions: This bi-directional MR study supports that host response to SARS-CoV-2 viral infection plays a role in the causal association with increased risk of hypothyroidism. Long-term follow-up studies are needed to confirm the expected increased hypothyroidism risk.

Thyroid Stimulating Hormone Levels Are Associated With Genetically Predicted Nonalcoholic Fatty Liver Disease

CONTEXT

Findings from observational studies indicate an association of thyroid hormone levels with the risk of nonalcoholic fatty liver disease (NAFLD); however, conflicting results remain and reverse causality may be a possibility.

OBJECTIVE

This study aimed to evaluate the associations between NAFLD and both plasma thyroxine (T4) and thyroid stimulating hormone (TSH) at the phenotypic and genetic levels.

METHODS

We included 14 797 participants, aged 20 to 74 years who had undergone abdominal ultrasonography during the Third National Health and Nutrition Examination Survey (NHANES III). Multivariable logistic regression analyses were used to examine the observational associations of TSH and T4 with NAFLD. Mediation analyses were performed to study whether the relationship between NAFLD and TSH levels was mediated via potential confounders. A bidirectional, two-sample Mendelian randomization (MR) analysis was used to determine the potential causal relationship.

RESULTS

Multivariable logistic regression model suggested a "dose-response" relationship between TSH (Q4 vs Q1: OR = 1.29; 95% CI, 1.10-1.52; Ptrend = 0.001) and NAFLD. BMI and ALT partially mediated the association between TSH and NAFLD, while the proportion of the mediation effects of BMI and ALT were 39.1% and 22.3%, respectively. In MR analyses, the inverse-variance weighted method was selected as primary method and suggested a putative causal effect of NAFLD on serum TSH levels (OR = 1.022; 95% CI, 1.002-1.043). The result was further validated in the sensitivity analyses.

CONCLUSION

Circulating TSH levels were associated with the risk of NAFLD. MR analysis suggested a putative causal effect of NAFLD on TSH levels.

New Horizons: the value of UK Biobank to research on endocrine and metabolic disorders

UK Biobank is an intensively characterized prospective study of 500 000 men and women, aged 40 to 69 years when recruited, between 2006 and 2010, from the general population of the United Kingdom. Established as an open-access resource for researchers worldwide to perform health research that is in the public interest, UK Biobank has collected (and continues to collect) a vast amount of data on genetic, physiological, lifestyle, and environmental factors, with prolonged follow-up of heath conditions through linkage to administrative electronic health records. The study has already demonstrated its unique value in enabling research into the determinants of common endocrine and metabolic diseases. The importance of UK Biobank, heralded as a flagship project for UK health research, will only increase over time as the number of incident disease events accrue, and the study is enhanced with additional data from blood assays (such as whole-genome sequencing, metabolomics, and proteomics), wearable technologies (including physical activity and cardiac monitors), and body imaging (magnetic resonance imaging and dual-energy X-ray absorptiometry). This unique research resource is likely to transform our understanding of the causes, diagnosis, and treatment of many endocrine and metabolic disorders.

Targeted Next-Generation Sequencing of Congenital Hypothyroidism-Causative Genes Reveals Unexpected Thyroglobulin Gene Variants in Patients with Iodide Transport Defect

Congenital iodide transport defect is an uncommon autosomal recessive disorder caused by loss-of-function variants in the sodium iodide symporter (NIS)-coding SLC5A5 gene and leading to dyshormonogenic congenital hypothyroidism. Here, we conducted a targeted next-generation sequencing assessment of congenital hypothyroidism-causative genes in a cohort of nine unrelated pediatric patients suspected of having a congenital iodide transport defect based on the absence of 99mTc-pertechnetate accumulation in a eutopic thyroid gland. Although, unexpectedly, we could not detect pathogenic SLC5A5 gene variants, we identified two novel compound heterozygous TG gene variants (p.Q29* and c.177-2A>C), three novel heterozygous TG gene variants (p.F1542Vfs*20, p.Y2563C, and p.S523P), and a novel heterozygous DUOX2 gene variant (p.E1496Dfs*51). Splicing minigene reporter-based in vitro assays revealed that the variant c.177-2A>C affected normal TG pre-mRNA splicing, leading to the frameshift variant p.T59Sfs*17. The frameshift TG variants p.T59Sfs*17 and p.F1542Vfs*20, but not the DUOX2 variant p.E1496Dfs*51, were predicted to undergo nonsense-mediated decay. Moreover, functional in vitro expression assays revealed that the variant p.Y2563C reduced the secretion of the TG protein. Our investigation revealed unexpected findings regarding the genetics of congenital iodide transport defects, supporting the existence of yet to be discovered mechanisms involved in thyroid hormonogenesis.

Causal Association of Thyroid Signaling with C-Reactive Protein: A Bidirectional Mendelian Randomization

Methods

Based on the latest genome-wide association study summary data, bidirectional two-sample Mendelian randomization (MR) was employed to detect the causal relationship and effect direction between TSH, fT4, and CRP. Furthermore, in view of obesity being an important risk factor of CVD, obesity trait waist-hip ratio (WHR) and body mass index (BMI) were treated as the research objects in MR analyses for exploring the causal effects of TSH and fT4 on them, respectively.

Results

Genetically increased CRP was associated with increased TSH (β = −0.02, P = 0.011) and with increased fT4 (β = 0.043, P = 0.001), respectively, but there was no evidence that TSH or fT4 could affect CRP. In further analyses, genetically increased TSH was associated with decreased WHR (β = −0.02, P = 3.99e − 4). Genetically increased WHR was associated with decreased fT4 (β = −0.081, P = 0.002). Genetically increased BMI was associated with increased TSH (β = 0.03, P = 0.028) and with decreased fT4 (β = −0.078, P = 1.05e − 4). Causal associations of WHR and BMI with thyroid signaling were not supported by weighted median analysis in sensitivity analyses.

Conclusion

TSH and fT4 were increased due to the higher genetically predicted CRP. WHR was decreased due to the higher genetically predicted TSH. These findings will provide reference for the prevention and treatment of inflammation and metabolic syndrome.

Identifying Causes of Fracture Beyond Bone Mineral Density: Evidence From Human Genetics

New therapies may help to prevent osteoporotic fractures other than through increasing bone mineral density (BMD). Because fracture risk has an important genetic component, we aim to identify loci increasing fracture risk that do not decrease BMD, using a recently-proposed structural equation model adapted to remove genetic influences of BMD on fracture risk. We used summary statistics of the largest genome-wide association studies (GWASs) for BMD and for fracture in these analyses. We next estimated the genetic correlation between the non-BMD or BMD-related genetic effects and other clinical risk factors for fracture. Last, based on white British participants in the UK Biobank, we conducted genetic risk score analyses to assess whether the aggregated genetic effects conferred increased major osteoporotic fracture risk. We found that only three loci affecting fracture risk exhibited genetic effects not mediated by BMD: SOST, CPED1-WNT16, and RSPO3, while these three loci simultaneously conferred BMD-related effects. No strong genetic associations between non-BMD or BMD-related effects and 16 clinical risk factors were observed. However, non-BMD effects might be genetic correlated with hip bone size. In the UK Biobank, a 1 standard deviation (1-SD) increase in the non-BMD genetic risk score conferred an odds ratio of 1.17 for incident major osteoporotic fracture, compared to 1.29 by a BMD-related genetic risk score. Our study suggests that the majority of common genetic predisposition toward fracture risk acts upon BMD. Although non-BMD genetic effects may exist, they are not strongly correlated with most traditional clinical risk factors. Risk loci harboring non-BMD genetic effects may influence other perspectives of bone quality, or confer effects that existing GWASs fail to capture, but they demonstrate weaker impact on fracture risk than BMD-related genetic effects. These findings suggest that most successful drug development programs for osteoporosis should focus on pathways identified through BMD-associated loci. © 2022 American Society for Bone and Mineral Research (ASBMR).

Thyroid Hormones Enhance the Growth of Estrogen Receptor-Positive Breast Cancers

BACKGROUND

Thyroid disorders have been associated with breast cancer. In fact, breast cancer is the most common secondary malignancy in female patients with thyroid cancer (1). Moreover, hyperthyroidism is associated with an 11% increased risk of breast cancer in women (2). Importantly, up to 30% of patients with breast cancer are treated with thyroid hormone replacement therapy (THRT) for overt or subclinical hypothyroidism (3). These observations, coupled with the preclinical data showing growth stimulatory effects of thyroid hormones (THs) in various cancer models (4), formed the rationale for the study by Wahdan-Alaswad et al. aimed at investigating the role of THRT on the outcome of patients with nonmetastatic breast cancer (3).

METHODS

The authors conducted an observational study analyzing the association between THRT, disease-free survival (DFS), and disease-specific survival (DSS) in two cohorts of patients with nonmetastatic breast cancer. The first cohort consistent of 820 patients followed for a median of 10 years and treated for breast cancer between 1962 and1993, with THRT implemented in 69 patients. The second cohort included 160 patients treated more recently (between 2006 and 2009) and followed for a median of 8.8 years, with 50 patients exposed to THRT. The data on the age, tumor size, presence or absence of steroid (estrogen and/or progesterone) receptors (SR+/SR-), and treatment regimen were incorporated in the multivariate model analyzing the association between DFS/DSS and THRT at 5 and 10 years. To better understand the results of the observational cohort study, the authors performed functional in vitro and in vivo experiments to investigate the molecular mechanisms underlying TH effects on breast cancer cells and to test the interactions between estrogen receptors (ERs) and TH receptors (THRs).

RESULTS

In patients with SR+ breast cancer, THRT was associated with a significantly increased risk of recurrence (DFS RR, 2.9; P<0.001) and death (DSS RR, 3.4; P<0.001), independent of age, tumor size and grade, while THRT in patients with SR- breast cancer was not associated with worse outcomes. Moreover, patients with SR+ breast cancer undergoing therapy with aromatase inhibitor combined with THRT were characterized by a shorter DFS (P<0.042) and a higher 10-year recurrence rate of 14%, as compared with 2% in patient treated with the aromatase inhibitor alone.The functional in vitro and in vivo studies revealed growth stimulatory effects of monotherapy with TH or estrogens that were further potentiated with combination therapy in ER+ breast cancer cell lines and mice xenografts. The RNA-Seq analysis revealed that combination therapy was associated with a significant activation of the cell cycle, mismatch repair, homologous recombination, and DNA replication signaling, as well as induced thyroid-specific genes and estrogen-mediated signatures. These effects were abrogated by the knockdown or inhibition of ER and/or THRa, suggesting that cross-talk and nuclear colocalization of ERs and THRs are major drivers of pro-oncogenic signaling in the ER+ breast cancer model.

CONCLUSIONS

The study reveals clinically significant associations between THRT and worse outcomes in patients with nonmetastatic SR+ breast cancer that are likely driven by interactions between the nuclear ERs and THRs, leading to upregulation of pro-oncogenic signaling. These results suggest that overuse of THRT in patients with hypothyroidism and concurrent breast cancer should be avoided.

An Atlas of Genetic Correlations and Genetically Informed Associations Linking Psychiatric and Immune-Related Phenotypes

Importance

Certain psychiatric and immune-related disorders are reciprocal risk factors. However, the nature of these associations is unclear.

Objective

To characterize the pleiotropy between psychiatric and immune-related traits, as well as risk factors of hypothesized relevance.

Design, Setting, and Participants

This genetic association study was conducted from July 10, 2020, to January 15, 2022. Analyses used genome-wide association (GWA) statistics related to 14 psychiatric traits; 13 immune-related phenotypes, ie, allergic, autoimmune, and inflammatory disorders; and 15 risk factors related to health-related behaviors, social determinants of health, and stress response. Genetically correlated psychiatric-immune pairs were assessed using 2-sample mendelian randomization (MR) with sensitivity analyses and multivariable adjustment for genetic associations of third variables. False discovery rate correction (Q value < .05) was applied for each analysis.

Exposures

Genetic associations.

Main Outcomes and Measures

Genetic correlations and MR association estimates with SEs and P values. A data-driven approach was used that did not test a priori planned hypotheses.

Results

A total of 44 genetically correlated psychiatric-immune pairs were identified, including 31 positive correlations (most consistently involving asthma, Crohn disease, hypothyroidism, and ulcerative colitis) and 13 negative correlations (most consistently involving allergic rhinitis and type 1 diabetes). Correlations with third variables were especially strong for psychiatric phenotypes. MR identified 7 associations of psychiatric phenotypes on immune-related phenotypes that were robust to multivariable adjustment, including the positive association of (1) the psychiatric cross-disorder phenotype with asthma (odds ratio [OR], 1.04; 95% CI, 1.02-1.06), Crohn disease (OR, 1.09; 95% CI, 1.05-1.14), and ulcerative colitis (OR, 1.09; 95% CI, 1.05-1.14); (2) major depression with asthma (OR, 1.25; 95% CI, 1.13-1.37); (3) schizophrenia with Crohn disease (OR, 1.12; 95% CI, 1.05-1.18) and ulcerative colitis (OR, 1.14; 95% CI, 1.07-1.21); and a negative association of risk tolerance with allergic rhinitis (OR, 0.77; 95% CI, 0.67-0.92).

Conclusions and Relevance

Results of this genetic association study suggest that genetic liability for psychiatric disorders was associated with liability for several immune disorders, suggesting that vertical pleiotropy related to behavioral traits (or correlated third variables) contributes to clinical associations observed in population-scale data.

Appraising the Effect of Potential Risk Factors on Thyroid Cancer: A Mendelian Randomization Study

CONTEXT

Various risk factors have been associated with the risk of thyroid cancer in observational studies. However, the causality of the risk factors is not clear given the susceptibility of confounding and reverse causation.

OBJECTIVE

A 2-sample Mendelian randomization approach was used to estimate the effect of potential risk factors on thyroid cancer risk.

METHODS

Genetic instruments to proxy 55 risk factors were identified by genome-wide association studies (GWAS). Associations of these genetic variants with thyroid cancer risk were estimated in GWAS of the FinnGen Study (989 cases and 217 803 controls). A Bonferroni-corrected threshold of P = 9.09 × 10-4 was considered significant, and P < 0.05 was considered to be suggestive of an association.

RESULTS

Telomere length was significantly associated with increased thyroid cancer risk after correction for multiple testing (OR 4.68; 95% CI, 2.35-9.31; P = 1.12 × 10-5). Suggestive associations with increased risk were noted for waist-to-hip ratio (OR 1.85; 95% CI, 1.02-3.35; P = 0.042) and diastolic blood pressure (OR 1.60; 95% CI, 1.08-2.38; P = 0.019). Suggestive associations were noted between hemoglobin A1c (HbA1c) (OR 0.20; 95% CI, 0.05-0.82; P = 0.025) and decreased risk of thyroid cancer. Risk of thyroid cancer was not associated with sex hormones and reproduction, developmental and growth, lipids, diet and lifestyle, or inflammatory factors (All P > 0.05).

CONCLUSION

Our study identified several potential targets for primary prevention of thyroid cancer, including central obesity, diastolic blood pressure, HbA1c, and telomere length, which should inform public health policy.

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.

Analysis of rare genetic variation underlying cardiometabolic diseases and traits among 200,000 individuals in the UK Biobank

Cardiometabolic diseases are the leading cause of death worldwide. Despite a known genetic component, our understanding of these diseases remains incomplete. Here, we analyzed the contribution of rare variants to 57 diseases and 26 cardiometabolic traits, using data from 200,337 UK Biobank participants with whole-exome sequencing. We identified 57 gene-based associations, with broad replication of novel signals in Geisinger MyCode. There was a striking risk associated with mutations in known Mendelian disease genes, including MYBPC3, LDLR, GCK, PKD1 and TTN. Many genes showed independent convergence of rare and common variant evidence, including an association between GIGYF1 and type 2 diabetes. We identified several large effect associations for height and 18 unique genes associated with blood lipid or glucose levels. Finally, we found that between 1.0% and 2.4% of participants carried rare potentially pathogenic variants for cardiometabolic disorders. These findings may facilitate studies aimed at therapeutics and screening of these common disorders.