Skip to main content

Causal inference of sex hormone-binding globulin on venous thromboembolism: evidence from Mendelian randomisation

Abstract

Background

Previous cohort studies have shown that exogenous sex hormone use, such as testosterone replacement therapy and oestrogen-containing contraceptives, can increase the risk of venous thromboembolism (VTE). However, the relationship between endogenous sex hormone levels and VTE remains unclear. The goal of the present study was to explore the causal roles of endogenous sex hormones, including hormone-binding globulin (SHBG), bioactive testosterone (BT), and total testosterone (TT), in VTE and its two subgroups, deep vein thrombosis (DVT) and pulmonary embolism (PE).

Methods

We used a genome-wide association study of sex hormones as exposure data and Finnish VTE data as the outcome. Inverse variance weighting, MR-Egger, and weighted median were used for two-sample Mendelian randomisation (MR). Sensitivity analyses included MR-Egger, MR-PRESSO, Cochrane Q test, MR Steiger, leave-one-out analysis, and funnel plot, combined with multivariate MR and replicated MR analyses using larger VTE data from the global biobank meta-analysis initiative. Linkage disequilibrium score regression (LDSC) was used to determine genetic associations and estimate sample overlap.

Results

Our findings genetically predicted that an increase in serum SHBG levels by one standard deviation (SD) caused 25% higher odds for VTE (OR: 1.25, 95% CI: 1.01−1.55) and 58% higher odds for PE (OR: 1.58, 95% CI: 1.20−2.08). LDSC supported the genetic correlation between these two traits and replicated analyses confirm SHBG’s genetic effect on VTE in both sexes (OR: 1.46, 95% CI: 1.20−1.78) and in females (OR: 1.49, 95% CI: 1.17−1.91). In addition, an increase in serum TT levels by one SD caused 32% higher odds for VTE (OR: 1.32, 95% CI: 1.08−1.62) and 31% higher odds for DVT (OR: 1.31, 95% CI: 1.01−1.69); however, LDSC and replicated analyses did not find a genetic correlation between TT and VTE or its subtypes. No significant correlation was observed between BT and all three outcome traits.

Conclusion

Our study provides evidence that elevated serum SHBG levels, as predicted by genetics, increase VTE risk. However, the causal effect of testosterone levels on VTE requires further investigation.

Background

Venous thromboembolism (VTE), which includes interrelated diseases such as deep vein thrombosis (DVT) and pulmonary embolism (PE), is among the top five most common vascular diseases in many regions [1]. In Western populations, the lifetime incidence of VTE is approximately one in 12 individuals. The survival rate after VTE diagnosis is much lower than expected; approximately 20% of individuals die within 1 year of VTE diagnosis, from VTE or other cardiovascular diseases, neoplastic conditions, or respiratory system disorders; and complications tend to manifest frequently in survivors [2,3,4]. The pathogenesis of VTE is complex and involves interactions among genetic, environmental, and lifestyle factors. Identifying underlying causal factors is important for developing effective preventive and therapeutic strategies.

Exogenous sex hormones such as hormone replacement therapy and oral oestrogen-containing contraceptives increase the risk of VTE in women [3, 5]. Additionally, some studies have shown that testosterone therapy may increase the risk of VTE in men [6, 7]. Regarding endogenous sex hormones, substantial evidence showed that endogenous oestrogen can reduce thrombosis by regulating coagulation and inflammation and inhibiting the vascular injury response [8, 9]. However, limited studies are available on the impact of endogenous testosterone and sex hormone-binding globulin (SHBG) on the risk of VTE. A prospective study incorporating 3,051 postmenopausal women and 3,925 middle-aged to older men without hormone replacement therapy (HRT) suggested that endogenous testosterone was not associated with VTE risk in a median follow-up of 17.6 years [10]. An analogous prospective study incorporated 4658 women and 4673 men and found that endogenous testosterone was not associated with VTE and its two subtypes during a follow-up period of 21 years [11]. However, the effects of SHBG on VTE risk remain unclear. Since its identification by Mercier et al. in 1966 [12], SHBG has been regarded as a glycoprotein synthesised in the liver which exhibits a strong binding affinity to circulating sex hormones, including 5α-dihydrotestosterone, testosterone, and 17β-oestradiol protein, and plays a crucial role in regulating their bioavailability at the target location [13,14,15]. However, recent studies have yielded new perspectives on the role of SHBG, suggesting that it can serve as both a biomarker and a potential drug candidate for multiple diseases, independent of sex hormones [16,17,18]. Odlind et al. first proposed SHBG as a standalone prognosticator of VTE risk following the administration of hormonal contraceptives [19], and subsequent studies supported this view [20, 21]. A case-control study conducted by Luuk et al. suggested that the risk of VTE shows a dose-dependent response to serum SHBG levels in women aged ≤ 45 years (369 cases and 296 controls) without exposure to exogenous hormones or pregnancy [22]. In addition, Raps et al. suggested that SHBG can be a risk marker of VTE for contraceptive preparations based on the positive correlation between SHBG and the normalized activate protein C (APC) sensitivity ratio and VTE risk in contraceptive users [20]. However, some studies suggested that SHBG has no significant effect on VTE. The aforementioned prospective study found no significant correlation between endogenous SHBG levels and VTE [10]. Frank et al. suggested that SHBG levels are not a useful risk marker for VTE in women using oral contraceptives [23]. Therefore, it remains unclear whether endogenous testosterone and SHBG levels are associated with VTE risk. Although some observational studies have small sample sizes, their results are contradictory and may be confounded by multiple confounding factors as well as reverse causalities. Moreover, the causal relationships between SHBG, endogenous testosterone, and VTE remain unclear. Hence, there is a need for a rigorous investigation of the causal effects of SHBG and testosterone levels on VTE risk.

Mendelian randomisation (MR) is a potent statistical method that establishes causal links between exposures and outcomes using single nucleotide polymorphisms (SNPs) as instrumental variables (IVs) [24]. In contrast to observational studies, MR can provide more robust evidence of causality by circumventing the intrinsic limitations of these studies, such as the risk of confounding, reverse causation, and measurement errors, by leveraging randomly assigned genetic variants during gametogenesis [25]. Although randomised controlled trials provide the best clinical evidence compared with MR research, they have limitations such as large investment in manpower and material resources, long follow-up time, and difficulty in implementation [25]. Therefore, the objective of this study was to further explore the potential causal effects of SHBG and testosterone on VTE, including its two subtypes, PE and DVT, using two-sample MR within the largest genetic population samples of SHBG and testosterone.

Methods

Study design

We conducted a two-sample MR study to assess whether there is a causal effect of genetically predicted SHBG, total testosterone (TT), and bioactive testosterone (BT) on VTE and its subtypes PE and DVT. MR analysis utilizes genetic variation as a source of natural randomisation, and SNPs, which represent worldwide human genetic diversity, were chosen as IVs that are not subject to environmental or behavioural selection bias. There were three fundamental hypotheses of the MR analysis: (1) IVs are highly correlated with the exposure variable; (2) IVs remain unaffected by any confounding factors; and (3) IVs only impact the outcome variable through the exposure variable. This includes relevance, independence, and exclusion restrictions [26]. A diagram of the two-sample MR study is shown in Fig. 1A. Second, MR sensitivity analyses were performed using MR-Egger, MR-PRESSO, Cochrane Q test, MR-Steiger, leave-one-out, funnel plot, and Multivariable Mendelian randomisation MR (MVMR). Then, replicated MR analyses using larger VTE data from the Global Biobank Meta-analysis Initiative (GBMI) were performed to verify the results of the primary MR analyses. Finally, linkage disequilibrium score regression (LDSC) was used to determine genetic associations and estimate sample overlap. A diagram of the study design is shown in Fig. 1B.

Fig. 1
figure 1

Diagram of the two-sample MR study (A) and the study design (B). MR: Mendelian randomisation; IVs: Instrumental variables; SNPs: Single nucleotide polymorphisms; SHBG: Sex hormone-binding globulin; TT: Total testosterone; BT: Bioactive testosterone; VTE: Venous thromboembolism; PE: Pulmonary embolism; DVT: Deep vein thrombosis; LDSC: Linkage disequilibrium score regression; MVMR: Multivariable Mendelian randomisation

Exposure data source

The exposure data for SHBG, TT, and BT were procured from publicly available genome-wide association studies (GWAS) provided by Ruth et al. [16] using data from UK Biobank, a national resource that has been described extensively elsewhere (https://www.ukbiobank.ac.uk) [27]. Individuals aged 40 to 69 years (men and women), who were enrolled in the National Health Service and residing within a maximum distance of 25 miles from any of the 22 designated study assessment centres, were eligible for participation between 2006 and 2010 [16]. Serum SHBG and TT levels were quantified in nmol/L using a two-step sandwich immunoassay analysis and one-step competitive analysis, respectively, in 425,097 participants of European ancestry. BT (N = 382,988) was calculated as nmol/L from albumin and TT levels [16]. Additionally, given that previous studies have shown relationships between body mass index (BMI) and SHBG as well as between BMI and VTE, we used BMI-adjusted SHBG data provided by the authors to mitigate the influence of BMI as a potential confounder on the study’s outcomes.

Outcome data source

The outcome data on VTE were sourced from the R8 version of the FinnGen GWAS database (https://r8.finngen.fi, date of retrieval: 2023-03-25), which contains electronic health record data on over five million Finnish individuals, including hospital discharge diagnoses and prescribed medication information. The phenotypic codes for VTE, DVT, and PE were ‘I9_VTE’ (17,048 cases, 325,451 controls), ‘I9_PHLETHROMBDVTLOW’ (8.077 cases, 295,014 controls), and ‘I9_PULMEMB’ (8,170 cases, 333,487 controls), respectively. The FinnGen GWAS database incorporated a VTE population aged 10 to 100 years. The median ages at the first event were 60.20, 57.80, and 62.96 years for both sexes, women, and men, respectively, in VTE cases; 58.86, 56.99, and 60.84 years, respectively, in DVT cases; and 65.50, 64.28, and 66.76 years, respectively, in PE cases (https://risteys.finregistry.fi/). The diagnoses of these VTE events were defined using nationwide registries of hospital discharge and cause of death data, using the International Classification of Diseases (ICD) revisions 8, 9, and 10 (https://risteys.finregistry.fi/) [28]. Replication was performed using another larger VTE dataset with 915,868 participants of European ancestry provided by the Global Biobank Meta-analysis Initiative (GBMI, https://www.globalbiobankmeta.org/, date of retrieval: 2023-04-15). To initially explore the causal relationship between SHBG, testosterone, and VTE in specific sexes, we further used the female VTE dataset from the GBMI as the outcome for verification, which contained 454,926 female participants, mostly of European ancestry (92.4%) [29]. The disease endpoints in the GBMI were defined following the ICD-9 or ICD-10 codes. The average age of all participants in the GBMI was 39.91–69.86 years, and the standard deviation of age was 7.9–25 years; however, the GBMI did not provide age span information for the VTE population specifically [29].

Selection of instrumental variables

To select SNPs that could serve as valid genetically predicted IVs for SHBG and testosterone, we initially established a rigorous significance threshold of P < 5 × 10− 8. Second, to avoid linkage disequilibrium, we grouped SNPs (r2 ≤ 0.001 and kb = 10,000). Subsequently, SNPs with an F-statistic < 10 were excluded. The F-statistic for each SNP was computed using the formula F = (N-2) × R2 / (1-R2), where N represents the sample size, and R2 represents the proportion of variance in the exposure variable explained by each SNP. For palindromic SNPs, researchers have attempted to deduce positive-strand alleles by utilizing allele frequencies, removing SNPs with intermediate allele frequencies because their orientation concerning exposure and outcome in GWASs cannot be determined with certainty. The resulting SNPs were used as IVs in the MR analysis and the overall R2 of the IVs for each trait pair was calculated to assess the efficacy of the IVs.

Mendelian randomisation analysis

The analyses were performed using R software version 4.2.3, using the ‘Two-Sample MR’ and ‘MR-PRESSO’ packages [30]. Multiple MR methods have been used to examine causal relationships, including inverse variance weighting (IVW) [31], MR-Egger [32], weighted median [33], weighted mode, simple mode, and MR-pleiotropy residual sum and outlier (MR-PRESSO) [34]. Each method assumes different criteria for instrumental variable (IV) validity; however, the IVW method is widely regarded as the most robust approach. Essentially, IVW is a meta-analysis technique that integrates causal estimates from each IV through weighted aggregation to obtain the final estimate of the causal effect. Therefore, in this study, IVW was used as the principal methodology to ascertain the causal associations between exposure and outcome. The remaining methods were utilized as ancillary approaches to complement the main analysis or offer supplementary insights. If MR-PRESSO identified any outlier SNPs, they were initially excluded, followed by a reassessment of the remaining IVs. Scatter plots were constructed to show the causal estimates of the five MR methods. Single SNP analyses were performed and forest and funnel plots of the IVW and MR-Egger methods were used to describe the results. Odds ratios (ORs) and 95% confidence intervals (CI) were used to quantify their influence on VTE risk. A significance level of P < 0.05 was deemed indicative of a potential association. The Benjamini-Hochberg (BH) procedure was employed to control the false discovery rate for multiple comparisons by adjusting the P-values. To rule out reverse causality, we conduct an MR-Steiger directionality test.

Sensitivity analysis

For sensitivity analysis, we first employed two other MR methods, MR-Egger and median MR, which are known for their increased robustness against pleiotropy. Directionally concordant outcomes with IVW strengthened the validity of our causal inferences. We also used the MR-PRESSO test, which identifies and corrects for pleiotropy by detecting and removing outlier SNPs that violate the assumption of no horizontal pleiotropy [34]. Finally, we assessed pleiotropy by conducting a leave-one-out analysis wherein we systematically excluded each SNP and assessed its impact on the causal estimate. Cochrane’s Q-statistic, derived from the IVW method, was used to assess the heterogeneity among the estimates obtained from individual SNPs. Moreover, LDSCs have been used to assess genome-wide genetic associations between exposure and outcomes. Intercepts from the LDSC also indicated a potential sample overlap between the two GWAS. Finally, the two-sample MR results were validated by using an independent outcome dataset with a larger sample size. MVMR was performed to rule out bias due to the correlation between the sex hormone levels.

Results

Selection of instrumental variables

IVs for SHBG, TT, and BT were selected based on established quality control criteria after clumping and excluding SNPs with an F-statistic of less than 10, suggesting that the presence of a weak IV bias is unlikely to exert a significant effect. After harmonization of the alleles and effects between the exposure and outcome in primary analyses, there were 148,148,147 SNPs used for causal analyses of SHBG with VTE, PE, and DVT (R2 = 0.062); 138,140,140 SNPs for TT with VTE, PE, and DVT (R2 = 0.027); and 86,89,91 SNPs for BT with VTE, PE, and DVT (R2 = 0.019). In the replicated analyses, there were 147 SNPs used for causal analysis of SHBG with VTE (R2 = 0.060), 149 SNPs for TT with VTE (R2 = 0.031), 93 SNPs for BT with VTE (R2 = 0.019), 88 SNPs for SHBG in women (SHBGw) with VTE in women (VTEw) (R2 = 0.105), 193 SNPs for TT in women (TTw) with VTEw (R2 = 0.075), and 127 SNPs for BT in women (BTw) with VTEw (R2 = 0.052). The detailed characteristics of these IVs are shown in Additional File 1 (Supplementary Table S1).

Causal estimates of genetically predicted SHBG, TT and BT with VTE

We conducted MR analysis to investigate the causal effects of SHBG, TT, and BT on VTE, PE, and DVT. As a result, nine trait pairs were analysed, four of which exhibited statistical differences before BH correction (P-value in both IVW and MR-PRESSO < 0.05). An increase in serum SHBG levels by one standard deviation (SD) caused 25% higher odds of VTE (OR: 1.25, 95%, CI: 1.01−1.55, P(IVW) = 0.0391, P(PRESSO) = 0.0409), and 58% higher odds for PE (OR: 1.58, 95%, CI: 1.20−2.08, P(IVW) = 0.0011, P(PRESSO) = 0.0014). An increase in serum total testosterone levels by one SD caused 32% higher odds for VTE (OR: 1.32, 95%, CI: 1.08−1.62, P(IVW) = 0.0068, P(PRESSO) = 0.0075), and 31% higher odds for DVT (OR, 1.31; 95%, CI: 1.01−1.69, P(IVW) = 0.0400, P(PRESSO), 0.0483). Two of the four pairs (TT-VTE and SHBG-PE) exhibited significant differences after BH-correct. The MR Steiger test of directionality showed ‘TRUE’ for all nine trait pairs. All the results are shown in Fig. 2 and Additional File 1 (Supplementary Table S2). Forest plots are shown in Fig. 2 and Additional File 2 (Supplementary Fig. 1). Scatter and funnel plots of the four positive trait pairs are shown in Fig. 3, and those of the negative trait pairs are shown in Additional File 2 (Supplementary Fig. 2 and Supplementary Fig. 3, respectively).

Fig. 2
figure 2

Causal relationships between sex hormones and VTE in primary two-sample MR analyses. The effect of a one-standard deviation (SD) increase in sex hormone levels on the odds ratio (OR) of VTE is represented by the OR and its corresponding 95% confidence interval (CI). MR: Mendelian randomisation; Nsnp: Number of SNPs participated in the analysis; R2: The proportion of variance elucidated by instrumental variables (IVs) in the exposure variable in each trait pair; Pa: P-value in IVW after multiple comparison BH correction; P(PRESSO): P-value after outlier-corrected in MR-PRESSO; Action = 3: Correct causal direction in MR Steiger test of directionality; SHBG: Sex hormone-binding globulin; TT: Total testosterone; BT: Bioactive testosterone; VTE: Venous thromboembolism; PE: Pulmonary embolism; DVT: Deep vein thrombosis

Fig. 3
figure 3

Scatter plots and funnel plots of the primary two-sample MR analyses (positive trait pairs). (A-D) In the scatter plots, the horizontal axis denotes the impact of instrumental variables (IVs) on the sex hormones, and the vertical axis represents the effect of IVs on VTE. Each black dots represents an individual SNP and the vertical and horizontal lines represent its corresponding 95% confidence interval (CI). The slope of the line represents the estimated causal effect of the various MR methods. (E-H) The horizontal axis of the funnel plot represents the estimated effect of each SNP on the exposure variable (βIV), while the vertical axis of 1/SEIV reflects the precision or uncertainty of these estimates. MR: Mendelian randomisation; SHBG: Sex hormone-binding globulin; TT: Total testosterone; VTE: Venous thromboembolism; PE: Pulmonary embolism; DVT: Deep vein thrombosis

MR sensitivity analysis

When assessing horizontal pleiotropy, the different MR methods (IVW, MR-Egger, and median MR) showed directionally consistent results, as shown in Fig. 2 and Additional File 1 (Supplementary Table S2), and all MR-Egger regression intercepts remained at zero without deviation, as shown in Table 1 and Additional File 1 (Supplementary Table S2). Moreover, the IVs exhibited no signs of horizontal pleiotropy, as indicated by intercept P-values exceeding 0.05, as shown in Table 1 and Additional File 1 (Supplementary Table S2). The leave-one-out analysis confirmed that no particular IVs drove any causal relationships, as shown in the leave-one-out plot (Fig. 4 and Additional File 2, Supplementary Fig. 4). In the LDSC results, three trait pairs showed significant genetic correlations (P-LDSC ≤ 0.05, Rg ≥ 0.1): SHBG-VTE, SHBG-PE, and SHBG-DVT. None of the nine trait pairs showed significant sample overlap in LDSC (P ≥ 0.05). As shown in Table 2. Replicated MR analyses were performed to verify the causal effect of SHBG, TT, and BT on VTE using data from both sexes and females from the GBMI Database. An increase in serum SHBG levels by one SD caused 46% higher odds for VTE (OR: 1.46, 95%, CI: 1.20−1.78, P(IVW) = 0.0002, P(PRESSO) = 0.0002) in both sexes and 49% higher odds for VTE (OR: 1.49, 95%, CI: 1.17−1.91, P(IVW) = 0.0014, P(PRESSO) = 0.0011) in females specifically, as shown in Fig. 5 and Additional File 1 (Supplementary Table S2), which further confirmed the results of the two-sample MR. MVMR analysed the causality between sex hormones as potentially relevant exposure variables; VTE, PE, and DVT from the Finnish database; and VTE from the GBMI database as outcomes. With outcome data from the Finnish database, the results showed a significant causal effect of SHBG on VTE (P = 0.0258) and PE (P = 0.0028) but not on DVT, which is consistent with the primary two-sample MR analyses, with the difference being that BT also showed a significant causal effect on PE (P = 0.0445). Using outcome data from the GBMI database, the results showed a significant causal effect of SHBG on VTE in both sexes (P = 0.0053) and females (P = 0.0000), which is also consistent with the replicative two-sample MR analyses. Additionally, female data showed a significant causal effect of BT (P = 0.0003) and TT (P = 0.0012) on VTE, as shown in Fig. 6 and Additional File 1 (Supplementary Table S3). Cochran’s Q test was performed to assess heterogeneity, and the P-values were below 0.05 (Table 1 and Additional File 1, Supplementary Table S2), indicating heterogeneity among the IVs and supporting the adoption of a random-effects model in the IVW analyses for these instances.

Table 1 Pleiotropy and Heterogeneity Tests for the Primary Two-sample Mendelian Randomisation (MR) Analyses
Fig. 4
figure 4

Leave-one-out plots of primary (A-D) and replicated (E, F) two-sample MR analyses (positive trait pairs). Each SNP is sequentially excluded to evaluate its impact on the causal estimate. MR: Mendelian randomisation; SHBG: Sex hormone-binding globulin; TT: Total testosterone; VTE: Venous thromboembolism; PE: Pulmonary embolism; DVT: Deep vein thrombosis; SHBGw: Sex hormone-binding globulin in women; VTEw: Venous thromboembolism in women

Table 2 Genetic correlation and genetic overlap estimated by Linkage Disequilibrium Score regression (LDSC)
Fig. 5
figure 5

Causal relationships between sex hormones and VTE in replicated two-sample MR analyses. The effect of a one-standard deviation (SD) increase in sex hormone levels on the odds ratio (OR) of VTE is represented by the OR and its corresponding 95% confidence interval (CI). MR: Mendelian randomisation; Nsnp: Number of SNPs participated in the analysis; R2: The proportion of variance elucidated by instrumental variables (IVs) in the exposure variable in each trait pair; Pa: P-value in IVW after multiple comparison BH correction; P(PRESSO): P-value after outlier-corrected in MR-PRESSO; Action = 3: Correct causal direction in MR Steiger test of directionality; SHBG: Sex hormone-binding globulin; TT: Total testosterone; BT: Bioactive testosterone; SHBGw: Sex hormone-binding globulin in women; TTw: Total testosterone in women; BTw: Bioactive testosterone in women; VTE: Venous thromboembolism; VTEw: Venous thromboembolism in women

Fig. 6
figure 6

Causal relationships between sex hormones and VTE in multivariable Mendelian randomisation (MVMR). MR: Mendelian randomisation; Nsnp: Number of SNPs participated in analyses; SHBG: Sex hormone-binding globulin; TT: Total testosterone; BT: Bioactive testosterone; SHBGw: Sex hormone-binding globulin in women; TTw: Total testosterone in women; BTw: Bioactive testosterone in women; VTE: Venous thromboembolism; VTEgbmi: Venous thromboembolism from GBMI database in replicative MR analysis; VTEw: Venous thromboembolism in women; PE: Pulmonary embolism; DVT: Deep vein thrombosis

Discussion

In this study, we explored the genetic associations and causality of SHBG, TT, and BT with VTE, PE, and DVT using a combination of two-sample MR, MRPRESSO, and LDSC. In the primary MR analyses, we found that elevated serum SHBG levels increased the risk of VTE and PE, whereas elevated serum TT levels increased the risk of VTE and DVT. Furthermore, the LDSC and MR sensitivity analyses supported the causal effect of SHBG on VTE. However, the causal effects of testosterone levels have not been verified.

Sex hormone-binding globulin (SHBG) and venous thromboembolism (VTE)

Previous studies have shown that the use of exogenous oestrogen, androgens, and endogenous oestrogen increases the risk of VTE. However, the association between endogenous androgens such as SHBG and VTE remains controversial.

Numerous investigations have substantiated the independent prognostic value of SHBG in assessing VTE risk subsequent to the use of hormonal contraceptives [19,20,21]. A recent case-control study involving women aged ≤ 45 years revealed that the occurrence risk of VTE exhibits a dose-response association with SHBG levels, indicating an OR of up to 2.0 for individuals in the highest quartile of SHBG levels compared to those in the lowest quartile, suggesting the pathophysiological effects of SHBG VTE in young women not using contraceptives [22]. In this study, primary IVW before BH correction and MR-PRESSO suggested a causal effect of SHBG on VTE. The replicated MR analyses in the largest VTE database currently knownand MVMR also supported this result, combined with the results of LDSC, which supports the genetic correlation between SHBG and VTE. We suggest that SHBG is causally associated with VTE, consistent with the conclusions of previous studies. Simultaneously, owing to the differences in sex hormones and SHBG between genders, we used the largest GWAS data available on female VTE to confirm the causal relationship between SHBG and VTE in females. Unfortunately, there are currently no available male VTE data that do not significantly overlap with our exposure data (VTE data for both males and females can be obtained in UKBB, but our previous LDSC analysis suggests a high overlap with exposure data). In addition, it must be emphasized that in the original GWAS study of the exposure data, the authors analysed and pointed out that the genetic structure of SHBG levels is highly consistent between males and females [16]; therefore, the results of the MR analysis using population-wide SHBG genetic IVs are credible, we concluded that a causal relationship exists between SHBG and VTE. Among the subtypes of VTE, SHBG was causally associated with PE in the primary MR and MVMR analyses, which is consistent with the genetic correlation of LDSC. However, we recognize that DVT and PE are considered two clinical presentations of VTE, with approximately 50–70% of PE patients concurrently experiencing DVT and 50% of DVT patients having an associated asymptomatic PE noted on lung scans [35]. The conclusion that SHBG is causally associated with PE but not DVT requires further validation. The transition from DVT to PE involves intricate physiological regulation, and several risk factors which predict the occurrence of PE in DVT patients, such as hypertension, diabetes, long lying state, glucocorticoid therapy, and D-dimer levels, have been reported [36, 37]. Therefore, our findings may allow us to speculate that SHBG could potentially be associated with a higher risk of PE development in DVT patients within the context of VTE formation. However, clinical research on the high-risk factors and underlying mechanisms for the occurrence of PE in DVT patients are limited. Our speculation in this regard necessitates further validation through subsequent research endeavours.

However, the mechanisms underlying the regulation of VTE by SHBG remain unclear. Among women using hormonal contraceptives, investigators have considered SHBG a marker of the ‘oestrogenicity’ of contraceptives, and oestrogen is an established risk factor for VTE [19]. Some researchers have also hypothesized that hepatically metabolized hormonal contraceptives may influence the production of SHBG and coagulation factors, ultimately contributing to VTE [20]. In addition, SHBG can be regarded as a potent amplifier of steroid activity, and it is the main transporter of hydrophobic androgens in hydrophilic blood, but perhaps less so for oestrogens, which can also undergo glucuronidation and sulfation to facilitate their transportation to target organs [38]. Therefore, their effects may also be via testosterone [39].Furthermore, the interaction between SHBG and megalin might act as an additional autocrine-controlled mechanism, actively transporting steroids into cells and potentially modulating hormone signalling pathways [38]. It is also believed that SHBG may be implicated in the modulation of chronic inflammation, cellular proliferation, and lipid metabolism, and these processes can potentially impact the development of VTE independent of sex hormones [17, 40]. Regardless of the underlying mechanism, SHBG may be a potential clinical biomarker of VTE. Follow-up studies should focus on exploring the pathophysiological mechanisms of SHBG in VTE to look forward to identifying novel therapeutic targets.

Testosterone and venous thromboembolism (VTE)

Multiple observational studies have indicated a potential elevation in short-term VTE risk among men receiving testosterone therapy [6, 7, 41]. A recent case-crossover study demonstrated that testosterone therapy was associated with an elevated risk of VTE within 12 months of follow-up in men with and without hypogonadism, with an OR of 2.32 (95% CI, 1.97−2.74) and 2.02 (95% CI, 1.47−2.77), respectively [42]. The mechanism underlying this effect might be attributed to the elevation of Hct levels in men following testosterone treatment, which is significantly associated with VTE [7, 43]. The current observational findings do not support an association between endogenous testosterone levels and the risk of VTE in men or women [10, 11, 22, 44]. However, a recent MR study observed a positive association between endogenous testosterone, genetically predicted by mutations in the JMJD1C gene region, and thromboembolism (VTE, arterial embolism, and thrombosis) (OR: 2.09, 95% CI: 1.27−3.46) in men, whereas no such association was found in women, and no association was observed between endogenous testosterone, genetically predicted by mutations in the SHBG gene region, and thromboembolism [45]. However, similar to this study, the correlation between testosterone and VTE failed to obtain consistent results in validation studies. Coupled with the small sample size (3225 men), we believe that their conclusions need to be further verified. In this study, the results of both IVW and MR-PRESSO suggested a causal relationship between TT but not BT and VTE. However, the replicated MR analyses in both sexes and females did not support this result, nor did the MVMR data in either sex or LDSC. Additionally, in the original GWAS of exposure data, the authors highlighted that the genetic contribution to the variability in circulating testosterone levels varies considerably between males and females, with several variants having genome-wide significance for testosterone in opposite directions, thus emphasising that sex-disaggregated data are best used in genetic association analyses of testosterone [16]. Additionally, in our MVMR analyses of female data, TT and BT had a significant causal effect on VTE. Therefore, we cannot completely exclude the correlation between testosterone levels and VTE when conducting sex-specific analyses. In conclusion, although this study had a larger sample size than previous MR studies, given the negative results of the LDSC and replicated analyses and the lack of better available sex-disaggregated data, we considered whether there was a causal association between testosterone and VTE needs the support of further large-sample sex-disaggregated MR analysis results.

Limitations and prospects

We conducted the first MR study to confirm the genetically predicted causal relationship between SHBG and VTE; however, this study had a few limitations. First, because the three assumptions of MR studies cannot be precisely satisfied, corresponding bias cannot be eliminated. In terms of instrument relevance, we adopted a stricter P-value (P < 5 × 10− 8), calculated the F value of each SNP, and deleted SNPs with F ≤ 10. The R2 of the IVs used in all MR analyses ranged from 1.9 to 10.5%, thereby reducing the bias introduced by weak IVs. Horizontal pleiotropy was the main reason for the bias in instrument validity. We used MR-PRESSO to exclude outliers that may have pleiotropic effects and then performed a subsequent MR analysis. Although IVW was utilized as the principal methodology to ascertain causal associations, we ensured the consistency of direction using various MR analysis methods. In addition, the leave-one-out analyses showed no significant horizontal pleiotropy, and the MVMR analyses showed robust results, consistent with the main MR analysis of the overall results, which were mainly used to rule out bias owing to the correlation between sex hormones. Additionally, we performed replicated MR analyses using an independent outcome database and LDSC to confirm the genetic correlation and reliability of primary MR analyses. Reverse causality was ruled out using MR Steiger. Second, heterogeneity tests revealed some heterogeneity in this study, which may reflect the different responses of different populations or environments to research effects. In response, a random-effects model was applied for IVW analyses. Third, this study was conducted in a European population; therefore, generalising the results to other races requires further verification. Fourth, sex hormones differ between sexes, and the genetic IVs of testosterone also show large differences. However, owing to the lack of sex-specific outcome data that did not overlap with our exposure data, our sex-disaggregated data were only for women. Therefore, the correlation between testosterone levels and VTE was used as the reference in this study. However, because the genetic IVs of SHBG described by the original author were not significantly different between the sexes, the results of this study on SHBG remain credible. Further, large-sample sex-disaggregated MR analyses of sex hormones, especially testosterone, in VTE are necessary to determine their causal relationship. Therefore, SHBG may be a potential clinical biomarker for VTE, and follow-up studies should focus on exploring the pathophysiological mechanisms of SHBG in VTE to look forward to developing novel therapeutic targets.

Conclusions

Our MR study suggests that elevated serum SHBG levels predicted by genetics may increase the risk of VTE; however, the causal relationship between testosterone and VTE requires further research. SHBG may play a pathological role in VTE and can be used as a predictor of lifetime risk and a novel therapeutic target for VTE. However, further studies of this mechanism are required.

Data Availability

The exposure data supporting the conclusions of this article are openly accessible and available in the UK Biobank (https://www.nealelab.is/uk-biobank). The outcome data in primary MR analyses are openly accessible and available in R8 release of the FinnGen GWAS results (https://r8.finngen.fi). The outcome data in replicated MR analyses are openly accessible and available in GBMI database(https://www.globalbiobankmeta.org/).

Abbreviations

BH:

Benjamini-Hochberg

BMI:

Body mass index

BT:

Bioactive testosterone

CI:

Confidence interval

DVT:

Deep vein thrombosis

GBMI:

Global Biobank Meta-analysis Initiative

GWAS:

Genome-wide association study

ICD:

International Classification of Diseases

IVs:

Instrumental variables

IVW:

Inverse variance weighting

LDSC:

Linkage disequilibrium score regression

MR:

Mendelian randomisation

MR-PRESSO:

MR-pleiotropy residual sum and outlier

MVMR:

Multivariable Mendelian randomization

OR:

Odds ratio

R2 :

Percentage of the variation explained by the IVs

PE:

Pulmonary embolism

SD:

Standard deviation

SHBG:

Sex hormone-binding globulin

SNP:

Single nucleotide polymorphisms

TT:

Total testosterone

VTE:

Venous thromboembolism

References

  1. Wendelboe AM, Raskob GE. Global burden of Thrombosis: epidemiologic aspects. Circ Res. 2016;118:1340–7.

    CAS  PubMed  Google Scholar 

  2. Lutsey PL, Zakai NA. Epidemiology and prevention of venous thromboembolism. Nat Rev Cardiol. 2023;20:248–62.

    PubMed  Google Scholar 

  3. Heit JA. Venous thromboembolism: Disease burden, outcomes and risk factors. J Thromb Haemostasis. 2010;3:1611–7.

    Google Scholar 

  4. Pastori D, Cormaci VM, Marucci S, Franchino G, Del Sole F, Capozza A et al. A Comprehensive Review of Risk factors for venous thromboembolism: from epidemiology to pathophysiology. Int J Mol Sci. 2023;24.

  5. Ge SQ, Tao X, Cai LS, Deng XY, Hwang MF, Wang CL. Associations of hormonal contraceptives and infertility medications on the risk of venous thromboembolism, ischemic Stroke, and Cardiovascular Disease in women. J Investig Med. 2019;67:729–35.

    PubMed  Google Scholar 

  6. Martinez C, Suissa S, Rietbrock S, Katholing A, Freedman B, Cohen AT, et al. Testosterone treatment and risk of venous thromboembolism: population based case-control study. BMJ. 2016;355:i5968.

    PubMed  PubMed Central  Google Scholar 

  7. Ory J, Nackeeran S, Balaji NC, Hare JM, Ramasamy AR. Secondary Polycythemia in men receiving testosterone therapy increases risk of major adverse Cardiovascular events and venous thromboembolism in the First Year of Therapy. J Urol. 2022;207:1295–301.

    PubMed  Google Scholar 

  8. Brenner B. Haemostatic changes in pregnancy. Thromb Res. 2004;114:409–14.

    CAS  PubMed  Google Scholar 

  9. Migeon BR. The role of X inactivation and cellular mosaicism in women’s health and sex-specific Diseases. JAMA. 2006;295:1428–33.

    CAS  PubMed  Google Scholar 

  10. Roetker NS, MacLehose RF, Hoogeveen RC, Ballantyne CM, Basu S, Cushman M, et al. Prospective study of endogenous hormones and incidence of venous thromboembolism: the Atherosclerosis risk in communities Study. Thromb Haemost. 2018;118:1940–50.

    PubMed  PubMed Central  Google Scholar 

  11. Holmegard HN, Nordestgaard BG, Schnohr P, Tybjaerg-Hansen A, Benn M. Endogenous sex hormones and risk of venous thromboembolism in women and men. J Thromb Haemost. 2014;12:297–305.

    CAS  PubMed  Google Scholar 

  12. Mercier Cea. A testosterone binding globulin. In Proceedings of the Second Symposium on the Steroid Hormones. Excerpta Medica International Congress Series. 1966;101:212.

  13. Dunn JF, Nisula BC, Rodbard D. Transport of steroid hormones: binding of 21 endogenous steroids to both testosterone-binding globulin and corticosteroid-binding globulin in human plasma. J Clin Endocrinol Metab. 1981;53:58–68.

    CAS  PubMed  Google Scholar 

  14. Joseph DR. Structure, function, and regulation of androgen-binding protein/sex hormone-binding globulin. Vitam Horm. 1994;49:197–280.

    CAS  PubMed  Google Scholar 

  15. Goldman AL, Bhasin S, Wu FCW, Krishna M, Matsumoto AM, Jasuja R. A reappraisal of Testosterone’s binding in circulation: physiological and clinical implications. Endocr Rev. 2017;38:302–24.

    PubMed  PubMed Central  Google Scholar 

  16. Ruth KS, Day FR, Tyrrell J, Thompson DJ, Wood AR, Mahajan A, et al. Using human genetics to understand the Disease impacts of testosterone in men and women. Nat Med. 2020;26:252–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Simons P, Valkenburg O, Stehouwer CDA, Brouwers M. Sex hormone-binding globulin: biomarker and hepatokine? Trends Endocrinol Metab. 2021;32:544–53.

    CAS  PubMed  Google Scholar 

  18. Bourebaba N, Ngo T, Śmieszek A, Bourebaba L, Marycz K. Sex hormone binding globulin as a potential drug candidate for liver-related metabolic disorders treatment. Biomed Pharmacother. 2022;153:113261.

    CAS  PubMed  Google Scholar 

  19. Odlind V, Milsom I, Persson I, Victor A. Can changes in sex hormone binding globulin predict the risk of venous thromboembolism with combined oral contraceptive pills? Acta Obstet Gynecol Scand. 2002;81:482–90.

    PubMed  Google Scholar 

  20. Raps M, Helmerhorst F, Fleischer K, Thomassen S, Rosendaal F, Rosing J, et al. Sex hormone-binding globulin as a marker for the thrombotic risk of hormonal contraceptives. J Thromb Haemost. 2012;10:992–7.

    CAS  PubMed  Google Scholar 

  21. van Rooijen M, Silveira A, Hamsten A, Bremme K. Sex hormone-binding globulin-A surrogate marker for the prothrombotic effects of combined oral contraceptives. Am J Obstet Gynecol. 2004;190:332–7.

    PubMed  Google Scholar 

  22. Scheres LJJ, van Hylckama Vlieg A, Ballieux B, Fauser B, Rosendaal FR, Middeldorp S, et al. Endogenous sex hormones and risk of venous thromboembolism in young women. J Thromb Haemost. 2019;17:1297–304.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Stanczyk FZ, Grimes DA. Sex hormone-binding globulin: not a surrogate marker for venous thromboembolism in women using oral contraceptives. Contraception. 2008;78:201–3.

    CAS  PubMed  Google Scholar 

  24. Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG. Using published data in mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. 2015;30:543–52.

    PubMed  PubMed Central  Google Scholar 

  25. Hartwig FP, Davies NM, Hemani G, Davey Smith G. Two-sample mendelian randomization: avoiding the downsides of a powerful, widely applicable but potentially fallible technique. Int J Epidemiol. 2016;45:1717–26.

    PubMed  Google Scholar 

  26. Davies NM, Holmes MV, Davey Smith G. Reading mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ. 2018;362:k601.

    PubMed  PubMed Central  Google Scholar 

  27. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562:203–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Kurki MI, Karjalainen J, Palta P, Sipilä TP, Kristiansson K, Donner KM, et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature. 2023;613:508–18.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Zhou W, Kanai M, Wu KH, Rasheed H, Tsuo K, Hirbo JB, et al. Global Biobank Meta-analysis Initiative: powering genetic discovery across human Disease. Cell Genom. 2022;2:100192.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7:e34408.

    PubMed  PubMed Central  Google Scholar 

  31. Lawlor DA, Harbord RM, Sterne JA, Timpson N, Davey Smith G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. 2008;27:1133–63.

    PubMed  Google Scholar 

  32. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44:512–25.

    PubMed  PubMed Central  Google Scholar 

  33. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some Invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40:304–14.

    PubMed  PubMed Central  Google Scholar 

  34. Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from mendelian randomization between complex traits and Diseases. Nat Genet. 2018;50:693–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Bruni-Fitzgerald KR. Venous thromboembolism: an overview. J Vasc Nurs. 2015;33:95–9.

    PubMed  Google Scholar 

  36. Dong HY, Tong MS, Wang J, Liu Y, Tao GY, Petersen RH, et al. Risk factors for Pulmonary Embolism in Lung cancer patients with lower limb deep venous Thrombosis: a case-control study. Transl Lung Cancer Res. 2023;12:1539–48.

    PubMed  PubMed Central  Google Scholar 

  37. Jin ZY, Li CM, Zheng K, Qu H, Yang WT, Wen JH et al. Nomogram for predicting Pulmonary Embolism in gynecologic inpatients with isolated distal deep venous Thrombosis. Int J Gynaecol Obstet. 2023.

  38. Thaler MA, Seifert-Klauss V, Luppa PB. The biomarker sex hormone-binding globulin - from established applications to emerging trends in clinical medicine. Best Pract Res Clin Endocrinol Metab. 2015;29:749–60.

    CAS  PubMed  Google Scholar 

  39. Ramachandran S, Hackett GI, Strange RC. Sex hormone binding globulin: a review of its interactions with testosterone and age, and its impact on Mortality in Men with type 2 Diabetes. Sex Med Rev. 2019;7:669–78.

    PubMed  Google Scholar 

  40. Simó R, Sáez-López C, Barbosa-Desongles A, Hernández C, Selva DM. Novel insights in SHBG regulation and clinical implications. Trends Endocrinol Metab. 2015;26:376–83.

    PubMed  Google Scholar 

  41. Sharma R, Oni OA, Chen G, Sharma M, Dawn B, Sharma R, et al. Association between Testosterone Replacement Therapy and the incidence of DVT and Pulmonary Embolism: a retrospective cohort study of the Veterans Administration Database. Chest. 2016;150:563–71.

    PubMed  Google Scholar 

  42. Walker RF, Zakai NA, MacLehose RF, Cowan LT, Adam TJ, Alonso A, et al. Association of Testosterone Therapy with risk of venous thromboembolism among men with and without hypogonadism. JAMA Intern Med. 2020;180:190–7.

    CAS  PubMed  Google Scholar 

  43. Rezende SM, Lijfering WM, Rosendaal FR, Cannegieter SC. Hematologic variables and venous Thrombosis: red cell distribution width and blood monocyte count are associated with an increased risk. Haematologica. 2014;99:194–200.

    PubMed  PubMed Central  Google Scholar 

  44. Svartberg J, Braekkan SK, Laughlin GA, Hansen JB. Endogenous sex hormone levels in men are not associated with risk of venous thromboembolism: the Tromso study. Eur J Endocrinol. 2009;160:833–8.

    CAS  PubMed  Google Scholar 

  45. Luo S, Au Yeung SL, Zhao JV, Burgess S, Schooling CM. Association of genetically predicted testosterone with thromboembolism, Heart Failure, and Myocardial Infarction: mendelian randomisation study in UK Biobank. BMJ. 2019;364:l476.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to extend our gratitude to the participants and investigators associated with the UK Biobank data, FinnGen study, and GBMI database. due to its open and accessible data, it is convenient for us to further conduct various bioinformatics analyses such as MR.

Funding

National Key R&D Program of China, Grant/Award Number: 2018YFC1311505; Gansu Provincial Clinical Research Center for Cardiovascular Diseases, Grant/Award Number:18JR2FA005; Department of Science and Technology of Gansu Province, Grant/Award Number: 20JR5RA357; National Administration of Traditional Chinese Medicine “refractory hypertension”, Grant/Award Number: Medical Administrative Department of the National Administration of Traditional Chinese Medicine[2018] NO.3; The First Hospital of Lanzhou University, Grant/Award Number: ldyyyn2021-90.

Author information

Authors and Affiliations

Authors

Contributions

ZZ: Writing-Review, Editing, Supervision. MB: Writing-Review, Supervision. SPW: Research Design, Data Analysis, Visualization, Writing-Draft. YXW: Methodology, Visualization, Writing-Review. YP: Writing-Review, Editing. PL and DZ: Data Collection and Acquisition, Editing. BPZ: Data Collection, Visualization, Acquisition. PYZ: Methodology, Editing. All authors contributed to the article and approved the final manuscript.

Corresponding author

Correspondence to Zheng Zhang.

Ethics declarations

Ethics approval and consent to participate

The research conducted in this study solely utilized publicly available online data and did not involve the collection of original data. The ethical approval for each of the studies can be traced back to the primary publications. All studies adhered to the ethical principles outlined in the Declaration of Helsinki, and informed consent was obtained from all participants.

Consent for publication

Not applicable.

Competing interests

The authors affirm that there are no conflicts of interest to disclose.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, S., Wang, Y., Bai, M. et al. Causal inference of sex hormone-binding globulin on venous thromboembolism: evidence from Mendelian randomisation. Thrombosis J 21, 109 (2023). https://doi.org/10.1186/s12959-023-00553-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12959-023-00553-9

Keywords