Abstract

Objective. Observational studies show the correlation between thyroid dysfunction and risk of venous thromboembolism. However, the causal effects remain uncertain. Our study was conducted to evaluate whether thyroid function and dysfunction were causally linked to the risk of venous thromboembolism. Methods. Publicly available summary data of thyrotropin (TSH) and free thyroxine (FT4), hypothyroidism, and hyperthyroidism were obtained from the ThyroidOmics Consortium and the UK Biobank. With single nucleotide polymorphisms (SNPs) as instrumental variables, the casual effects of genetically predicted TSH and FT4 and hypo- and hyperthyroidism on venous thromboembolism outcome were estimated through Mendelian randomization analysis methods (inverse variance weighted (IVW), MR-Egger, weighted median, simple mode, and weighted mode). Cochran’s Q test was performed to evaluate the heterogeneity and horizontal pleiotropy. Results. Our study selected 15 FT4-, 36 TSH-, 3 hyperthyroidism-, and 79 hypothyroidism-associated SNPs as instrumental variables. The IVW analysis results showed that the odds ratio of venous thromboembolism for hyperthyroidism was 1.124 (95% confidence interval: 1.019-1.240; ), demonstrating the casual effect of hyperthyroidism not FT4, TSH, and hypothyroidism on venous thromboembolism. No heterogeneity or horizontal pleiotropy was observed according to Cochran’s Q test. Conclusion. Our Mendelian randomization analysis supports the causal effect of hypothyroidism on risk of venous thromboembolism. There is no evidence that genetically predicted TSH, FT4, and hypothyroidism have casual effects on venous thromboembolism. Future studies should be conducted to elucidate the underlying pathophysiological mechanisms.

1. Introduction

Venous thromboembolism, including deep vein thrombosis and pulmonary embolism, is a chronic disease that affects approximately 10 million individuals per year globally [1, 2]. The annual incidence of acute venous thromboembolism is 1-2 cases per 1000 people and is four times higher in high-income countries compared with low-income countries [3]. The diagnostic methods of venous thromboembolism consist of sequential examinations combining clinical probabilities, D-dimer test, and imaging assessments [4]. Venous thromboembolism therapy is aimed at preventing thrombus extension and embolism, cardiopulmonary failure, deaths, recurrence, and long-term complications [5]. Direct oral anticoagulants and thrombin inhibitor are the major treatment strategies in the management of venous thromboembolism [6]. Venous thromboembolism is a multifactorial illness induced by the interaction of multiple predisposing factors [4]. Observational studies have showed the correlation between thyroid dysfunction and risk of venous thromboembolism. A prospective cohort study showed that hyperthyroidism not hypothyroidism was associated with lower risk of recurrent venous thromboembolism among elderly patients [7]. A meta-analysis of cohort studies also showed that hyperthyroidism was a risk factor of venous thromboembolism [8]. However, in a retrospective study, increased risk of venous thromboembolism was observed among patients with hypothyroidism not hyperthyroidism [9]. Nevertheless, it remains unclear whether the above observed associations are causal, because observational studies could have selection bias, residual confounding, and reverse causality.

Mendelian randomization is a research method that uses genetic variation as an instrumental variable to establish a model to test the causal relationship between variable exposure factors and diseases [10, 11]. Single nucleotide polymorphism (SNP) loci discovered by genome-wide association study (GWAS) based on large samples provide a large number of genetic variation tool variables for Mendelian randomization analysis, which can effectively reduce the bias of causality estimates [1214]. Mendelian randomization has been widely used in venous thromboembolism studies, revealing the intrinsic link between various phenotypes and venous thromboembolism [1517]. For instance, a Mendelian randomization study showed that obesity is casually correlated with venous thromboembolism [18]. Another study found that taller height acts as a risk factor of venous thromboembolism [19]. Based on previous research, this study collected currently available GWAS data on thyroid function and dysfunction. The causal relationships of free thyroxine (FT4), thyrotropin (TSH), hyperthyroidism, and hypothyroidism with venous thromboembolism were assessed using two-sample Mendelian randomization, which might provide a basis for prevention and treatment of venous thromboembolism.

2. Materials and Methods

2.1. Summary Data and Study Population

Our study acquired publicly available genetic summary data from two large GWAS cohorts (the ThyroidOmics Consortium [20] and the UK Biobank [21]) (Table 1). GWAS summary data on FT4 and TSH were accessed from the ThyroidOmics Consortium, composed of 72,167 venous thromboembolism samples. Meanwhile, GWAS summary data on hyperthyroidism and hypothyroidism were obtained from the UK Biobank, comprising 337,159 venous thromboembolism samples.

2.2. Selection of Instrumental Variables

The screening criteria of FT4-, TSH-, hyperthyroidism-, and hypothyroidism-associated SNPs as instrumental variables were as follows: (i) ; (ii) linkage disequilibrium was removed ( and ).

2.3. Two-Sample Mendelian Randomization Estimates

Five two-sample Mendelian randomization analysis methods containing inverse variance weighted (IVW), MR-Egger, weighted median, simple mode, and weighted mode were applied for estimating the casual effects of FT4, TSH, hyperthyroidism, and hypothyroidism on risk of venous thromboembolism. Among the five approaches, IVW was the main Mendelian randomization analysis to evaluate the causal effects because this approach is stable and accurate when directional pleiotropy is absent. All results were displayed as odds ratio (OR) and 95% confidence interval (CI), and values < 0.05 were considered statistically significant. In addition, the results were visualized into forest and scatter plots of the relationships of FT4-, TSH-, hyperthyroidism-, and hypothyroidism-associated SNPs with risk of venous thromboembolism.

2.4. Heterogeneity and Pleiotropy Test

Mendelian randomization analysis exists heterogeneity because of the differences in platforms, experimental conditions, inclusion populations, and SNPs. Cochran’s Q test of IVW analysis was conducted to detect heterogeneity. values < 0.05 indicated significant heterogeneity.

2.5. Statistical Analysis

All statistical analysis was implemented utilizing Two-Sample MR package in R (version 3.6.2) [22].

3. Results

3.1. Selection of FT4-, TSH-, Hyperthyroidism-, and Hypothyroidism-Associated SNPs as IVs

This study was conducted to investigate the causal relationships of FT4, TSH, hyperthyroidism, and hypothyroidism with venous thromboembolism via adopting two-sample Mendelian randomization. Based on the ThyroidOmics Consortium, we identified SNPs associated with thyroid function (FT4 and TSH). According to , LD , and , we identified 15 FT4-associated SNPs as IVs, including rs10119187, rs10739496, rs10818937, rs10946313, rs11039355, rs113107469, rs17185536, rs2235544, rs225014, rs4149056, rs4842131, rs4954192, rs56069042, rs6785807, and rs9356988 (Table 2). Meanwhile, 36 TSH-associated SNPs were used as IVs, as follows: rs1042673, rs1045476, rs1079418, rs10814915, rs10917469, rs10957494, rs11159482, rs11255790, rs1157994, rs11732089, rs1203944, rs12284404, rs1265091, rs12893151, rs13015993, rs13329353, rs1663070, rs17020122, rs17477923, rs17767491, rs2127387, rs2439301, rs28502438, rs30227, rs334725, rs398745, rs4445669, rs4804413, rs4933466, rs59381142, rs7329958, rs8015085, rs9298749, rs9381266, rs9497965, and rs963384 (Table 3). Additionally, on the basis of UK Biobank, with the same criteria, we determined 3 hyperthyroidism-associated SNPs (rs2160215, rs28360997, and rs3087243; Table 4) as well as 79 hypothyroidism-associated SNPs as IVs (Table 5). As depicted in Figures 1(a)1(d), forest plots displayed the estimates of each FT4-, TSH-, hyperthyroidism-, and hypothyroidism-associated SNP on venous thromboembolism.

3.2. Two-Sample Mendelian Randomization Analysis for Evaluating Causal Effects of FT4, TSH, Hyperthyroidism, and Hypothyroidism on Venous Thromboembolism

Five different two-sample Mendelian randomization analysis approaches (IVW, MR-Egger, weighted median, sample mode, and weighted mode) were conducted to assess the casual effects of FT4, TSH, hyperthyroidism, and hypothyroidism on risk of venous thromboembolism, as shown in Figures 2(a)2(d). For the IVW analysis results, the OR of venous thromboembolism for hyperthyroidism was 1.124 (95% CI: 1.019-1.240; ), without statistical significance for the MR-Egger, weighted median, sample mode, and weighted mode analysis results (Table 6). Meanwhile, no statistical significance from the five Mendelian randomization analysis approaches was found for FT4, TSH, and hypothyroidism on venous thromboembolism. The above data demonstrated that hyperthyroidism was a risk factor of venous thromboembolism.

3.3. Heterogeneity and Pleiotropy Test

Through Cochran’s Q test of IVW, we evaluated the heterogeneity and horizontal pleiotropy. A significant heterogeneity was found among FT4-associated SNPs on venous thromboembolism (). Moreover, there was no heterogeneity among TSH- (), hyperthyroidism- (), or hypothyroidism-associated SNPs () on venous thromboembolism (Table 7). Further analysis showed that FT4- (), TSH- (), hyperthyroidism- (), and hypothyroidism- () associated SNPs had no horizontal pleiotropy between FT4, TSH, hyperthyroidism, and hypothyroidism and risk of venous thromboembolism. As shown in funnel plots, when a single FT4-, TSH-, hyperthyroidism-, or hypothyroidism-associated SNP was used as an IV, the points representing the causality were relatively symmetrically distributed, suggesting that the cause was unlikely to be influenced by potential bias (Figures 3(a)3(d)).

4. Discussion

Randomized controlled trials are the gold standard for causal inference in clinical research but are expensive and time-consuming and often difficult to implement due to ethical factors and subject constraints [23, 24]. Observational studies are relatively easy to implement and are widely applied for initial determination of etiology [25, 26]. Nonetheless, when there are reverse causal effects and potential confounding factors, the results of observational studies are often difficult to be recognized [27]. Mendelian randomization method provides an effective way to solve the above problems, using genetic variation as an instrumental variable effectively to avoid confounding factors caused by the environment [28]. Furthermore, compared with the immediate results obtained from randomized controlled trials, the exposure factors obtained from a genetic point of view are often accompanied for life and can even be passed on to the next generation [29]. Based on above evidence, we carried out the two-sample Mendelian randomization study to analyze the casual effects of FT4, TSH, hyperthyroidism, and hypothyroidism on venous thromboembolism.

The IVW method applies the ratio method to calculate causality estimates for individual instrumental variables and aggregates each estimate to perform a weighted linear regression to obtain an overall estimate [30]. However, this method cannot calculate the bias caused by pleiotropy, so an additional Cochran’s Q test is needed to evaluate pleiotropy [31]. Our IVW analysis results showed that the ORR of venous thromboembolism for hyperthyroidism was 1.124 (95% CI: 1.019-1.240; ), demonstrating the casual effect of hyperthyroidism not FT4, TSH, and hypothyroidism on venous thromboembolism. Additionally, there was no heterogeneity or horizontal pleiotropy in accordance with Cochran’s Q test. Genome-wide association analysis uncovers that most variations in thyroid function are genetically determined [20]. In a nationwide population-based cohort study, hypothyroidism was in relation to higher risk of venous thromboembolism [32]. Decreased venous thromboembolism risk of hypothyroid patients treated with thyroxine replacement therapy was not investigated. However, our Mendelian randomization analysis showed that there was no causal effect of thromboembolism on venous thromboembolism. Another cohort study showed that TSH levels within the normal range were not associated with venous thromboembolism risk, but abnormal TSH levels were associated with moderate risk of venous thromboembolism [33]. Case-control studies have shown the significant association between serum TSH levels and increased risk of venous thromboembolism [34]. Consistently, TSH concentration was an independent risk factor of venous thromboembolism regardless of thyroid function [35]. Our causal analysis demonstrated that TSH was not casually associated with venous thromboembolism.

However, there are several limitations in our study. Firstly, the sample sets included in this study are all from the European region. Whether the research conclusions can be extended to other populations still needs further research to be proven. Secondly, there may be population sharing of sample sets from Europe, leading to overuse of genetic information.

5. Conclusion

Altogether, in the two-sample Mendelian randomization study with genome-wide association variants, our evidence demonstrated that hypothyroidism was casually associated with venous thromboembolism outcome. However, there was no evidence that TSH, FT4, and hypothyroidism had casual effects on venous thromboembolism. Despite this, further studies are required to validate the biological mechanisms underlying the casual association.

Abbreviations

SNP:Single nucleotide polymorphism
GWAS:Genome-wide association study
FT4:Free thyroxine
TSH:Thyrotropin
IVW:Inverse variance weighted
OR:Odds ratio
CI:Confidence interval.

Data Availability

The datasets analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.