- Research
- Open access
- Published:
Causal inference of sex hormone-binding globulin on venous thromboembolism: evidence from Mendelian randomisation
Thrombosis Journal volume 21, Article number: 109 (2023)
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.
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).
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.
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
Wendelboe AM, Raskob GE. Global burden of Thrombosis: epidemiologic aspects. Circ Res. 2016;118:1340–7.
Lutsey PL, Zakai NA. Epidemiology and prevention of venous thromboembolism. Nat Rev Cardiol. 2023;20:248–62.
Heit JA. Venous thromboembolism: Disease burden, outcomes and risk factors. J Thromb Haemostasis. 2010;3:1611–7.
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.
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.
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.
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.
Brenner B. Haemostatic changes in pregnancy. Thromb Res. 2004;114:409–14.
Migeon BR. The role of X inactivation and cellular mosaicism in women’s health and sex-specific Diseases. JAMA. 2006;295:1428–33.
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.
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.
Mercier Cea. A testosterone binding globulin. In Proceedings of the Second Symposium on the Steroid Hormones. Excerpta Medica International Congress Series. 1966;101:212.
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.
Joseph DR. Structure, function, and regulation of androgen-binding protein/sex hormone-binding globulin. Vitam Horm. 1994;49:197–280.
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.
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.
Simons P, Valkenburg O, Stehouwer CDA, Brouwers M. Sex hormone-binding globulin: biomarker and hepatokine? Trends Endocrinol Metab. 2021;32:544–53.
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.
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.
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.
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.
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.
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.
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.
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.
Davies NM, Holmes MV, Davey Smith G. Reading mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ. 2018;362:k601.
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.
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.
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.
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.
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.
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.
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.
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.
Bruni-Fitzgerald KR. Venous thromboembolism: an overview. J Vasc Nurs. 2015;33:95–9.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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
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
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.
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.
About this article
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
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12959-023-00553-9