Zixiang Ye1,2
Jingang Zheng1,2*- 1Department of Cardiology, Peking University China-Japan Friendship School of Clinical Medicine, Beijing, China
- 2Department of Cardiology, China-Japan Friendship Hospital, Beijing, China
Objective: ADAMTS13 plays a crucial role in several diseases. Many observational studies have reported the relationship between ADAMTS13 and some cardiovascular diseases but have drawn different conclusions, likely attributed to confounding factors lacking adjustment. Identifying the role of ADAMTS13 in cardiovascular diseases is pivotal for prevention as well as early intervention in patients with latent cardiovascular diseases. This study aims to estimate whether the level and activity of ADAMTS13 are causally associated with common cardiovascular diseases.
Methods: We applied a two-sample Mendelian randomization approach incorporating genome-wide association summary statistics to verify the causal association between ADAMTS13 level, as well as activity and cardiovascular diseases.
Results: Lower ADAMTS13 activity was causally associated with the increased risks for coronary heart diseases (b = −0.0041, se = 0.0019, p < 0.05) as well as myocardial infarction (b = −0.0048, se = 0.0022, p < 0.05). Standard inverse-variance weighted Mendelian randomization results suggested no genetic support for a causal association between ADAMTS13 level and cardiovascular diseases including coronary heart disease, myocardial infarction, atrial fibrillation, heart failure, and venous thromboembolism (p > 0.05).
Conclusion: The causal effect of lower ADAMTS13 activity on the increased odds of having cardiovascular diseases was coronary heart disease and myocardial infarction.
Introduction
A disintegrin and metalloprotease with thrombospondin type 1 repeats 13 (ADAMTS13) is one of the most well-known members in the ADAMTS family (Zhong and Khalil, 2019). It is produced and secreted by hepatic stellate cells (Haegele et al., 2018) and endothelial cells (Wang et al., 2016). ADAMTS13 plays a pivotal role in platelet tethering and platelet adhesion (Crawley et al., 2011), putting platelet into smaller and less procoagulant forms (Levy et al., 2001). The multi-domain structure helps ADAMTS13 to function through allosteric effects (Zheng et al., 2001).
Some studies have reported that increased ADAMTS13 levels are risk factors for various diseases such as diabetes (Domingueti et al., 2013), preeclampsia (Cines and Levine, 2017), liver cirrhosis (Uemura et al., 2010), and thrombotic thrombocytopenic purpura (Levy et al., 2001) because of the imbalance between ADAMTS13 and von Willebrand factor (VWF) levels in circulation. However, epidemiologic evidence shows conflicting results on the association between ADAMTS13 level and activity and the risk of cardiovascular diseases such as myocardial infarction (MI), angina pectoris, and stroke (Danesh et al., 2004; Crawley et al., 2008; Horii et al., 2008; Peyvandi et al., 2010; Wieberdink et al., 2010; Hanson et al., 2011). Furthermore, because of the lack of direct evidence, it remains unclear whether ADAMTS13 is a causal independent risk factor of cardiovascular disease.
Nevertheless, the limited evidence of direct correlation between ADAMTS13 and cardiovascular diseases is due to the lack of comparability of interdisciplinary research, weak statistical correlation, and the inability to adjust for confounding factors. To minimize potential methodological limitations such as confounding and reverse causality, the causal relationship between these features requires a more effective method for inference.
Mendelian randomization (MR) is a statistical method to detect and quantify the causality in observational epidemiological studies using genetic instrumental variables (IVs) as proxies (Emdin et al., 2017). Given an appropriate instrument, MR studies, compared with traditional observational studies, will not be disturbed by confounding factors (Smith and Ebrahim, 2003). The schematic presentation of our MR is given in Figure 1. In the past decades, one-sample MR was conducted using genetic instruments, exposure, and outcome of interest from individuals measured in the same sample. However, when measuring exposure and outcome data in different samples, MR can also be used to assess potential causality, which is called a two-sample MR (Burgess et al., 2015). The publicly available genome-wide association study (GWAS) summary data can be utilized in two-sample MR analyses. It is greatly increasing analysis scope compared with one-sample MR (Hemani et al., 2016; Zheng et al., 2017). MR analysis commonly uses single-nucleotide polymorphisms (SNPs) as individual instruments to identify genetic variants related to disease risk (Davies et al., 2018).
Figure 1. Mendelian randomization analysis schematic diagram. Our study aims to estimate the causal relationship between risk factors and outcomes using SNPs as instrumental variables. SNPs should not be associated with any confounding factors. SNP, single-nucleotide polymorphism.
This study aimed to determine the causal effect of ADAMTS13 level and activity on common cardiovascular diseases by conducting a two-sample MR analysis.
Materials and Methods
Study Design
Based on previous research, the MR design was based on IV analysis. By using genetic variation as an IV for exposure, it can strengthen inferences about the causal nature of the association between ADAMTS13 and the outcome. Two-sample MR studies were conducted to determine the causal association between ADAMTS13 level and activity and outcomes based on the available summary level data from GWAS.
In our two-sample MR analysis, the genetic IVs that were applied as IVs should have satisfied the following three assumptions: (1) the IVs must be closely associated with the exposure; (2) IVs must not be associated with any known confounding factors; and (3) the selected IVs should be independent of the outcomes, exposure, and confounders.
Data Source
Summary-level data of genetic predictors of SNPs related to ADAMTS13 levels were extracted from the GWAS study where 153 loci were identified from 3,244 healthy individuals (Ma et al., 2017). The SNPs predicting ADAMTS13 activity were acquired from de Vries’ study wherein four SNPs were identified from 5,448 European individuals (de Vries et al., 2015). Furthermore, we extracted the summary statistics for coronary heart disease (CHD) and MI from the Coronary Artery Disease Genome-Wide Replication and Meta-analysis plus the Coronary Artery Disease Genetics (CARDIoGRAMplusC4D) consortium with 1000 Genomes-based genome-wide association meta-analysis of 48 studies (CARDIoGRAMplusC4D1) (Nikpay et al., 2015). Atrial fibrillation (AF) data were contributed by Neale lab analysis of UK Biobank phenotypes with 3,818 AF cases and 333,381 controls2. Heart failure (HF), ischemic stroke, and venous thromboembolism (VTE) data were also acquired from UK Biobank phenotypes with a sample size of 361,194.
Definition of Outcomes
CHD was defined as MI, acute coronary syndrome (ACS), angina pectoris, or coronary artery stenosis > 50%, medical history, clinical diagnosis determined at the time of the study, procedures, medications, symptoms indicating angina, or self-report diagnosis. MI was defined as (a) the appearance of persistent ischemic chest pain; (b) electrocardiogram showing ischemic changes with dynamic evolution; and (c) increases in the levels of cardiac biomarkers. For AF, cases were defined as individuals with paroxysmal or permanent AF or atrial flutter. HF was defined as self-reported HF, pulmonary edema, or cardiomyopathy or an international classification of disease codes that indicate any cause of heart, ventricular failure, or cardiomyopathy. The definition of ischemic stroke was mainly the same as that listed in the Trial of Org 10,172 in Acute Stroke Treatment (TOAST) criteria. VTE referred to both a deep venous thrombosis and a pulmonary embolism.
Single-Nucleotide Polymorphism Validation and Linkage Disequilibrium Assessment
The crucial prerequisite for any MR analysis is to maintain the independence of the included instruments related to the exposure. However, we can take measures in the two-sampling MR analysis to explain any correlations about linkage disequilibrium (LD). Clustering the reference GWAS data sets of samples with similar ancestry is an effective way to ensure the independence of all instruments. R software’s clumping procedure with the “TwoSampleMR” package was performed to automatically prune SNPs with linkage dependence. The comparisons for SNP pairs in the European population were applied via LDlink to ensure that all the IVs are selected without LD with each other (Machiela and Chanock, 2015). We retained the SNPs with smaller p-value among all SNP pairs, violating the independence assumption with r2 > 0.1.
Harmonizing Single-Nucleotide Polymorphism Effects
By assigning SNPs to the same effect allele, matching can be made between data sources. We must ensure that the impact and standard error of each identified SNP on exposure and outcome corresponds to alleles of the same effect. Common sources of unexpected bias include wrong effect alleles, palindromic SNPs, and incompatible alleles. When it is a palindrome SNP, it is difficult to identify effector alleles, because the letters on the forward and reverse strands are the same. In this case, the result alleles should be flipped to match the exposed alleles, and then the effect alleles should be aligned. Palindromic SNPs were deleted to prevent unexpected biases.
Pleiotropy Assessment
We performed the MR–Egger regression to evaluate the horizontal pleiotropy in which IVs influenced outcomes through more than one biological pathway. To examine whether there is evidence of publication bias, the MR–Egger regression is often used in meta-analysis as a recognized statistical method (Egger et al., 1997). The potential pleiotropic effect across the IVs using the MR–Egger method and intercepts that deviate from the origin may provide evidence for potential pleiotropic effects in genetic IVs. Furthermore, funnel charts were used to visually inspect symmetry, which has a similar function in the meta-analysis literature, indicating whether the causal estimates of weaker variants tend to tilt in one direction and where any deviation may suggest potential pleiotropic effects (Bowden et al., 2015).
Heterogeneity Assessment
To avoid the problem that the minor allele frequency (MAF) differences among various ancestries may contribute to SNPs related to both outcomes and ancestry, we enrolled SNPs from studies where at least the majority included European individuals in our MR analysis. However, educational background, lifestyle, and other confounding factors might still exist in the European subgroup. Therefore, we performed Cochran’s Q test to estimate the heterogeneity of SNPs, which is based on the weighted sum of the squared differences between the effect of a single SNP and the aggregate effect of all SNPs.
Mendelian Randomization Estimates
The role of each exposure feature in the outcomes was evaluated via a two-sample MR. In short, we selected the SNP closely related to the specific exposure (p < 5 × 108) as the IV, and then we obtained the corresponding effect estimate from the summary statistics. After that, a two-sample MR analysis was carried out by weighting the impact of each SNP on the outcomes via weighting the impact of each SNP on each exposure. Different models were then used to aggregate these estimates to provide a comprehensive summary statistic of the effect of ADAMTS13 level and activity. Nevertheless, MR analysis would be hampered if the SNPs selected in the MR were weak. Thus, we performed a robust adjusted profile score, unbiased when there are many weak instruments and robust to systematic and idiosyncratic pleiotropy, considering the measurement error in SNP-exposure effects. The weighted-median method was used to sort the MR estimates according to their inverse variance from the smallest to largest weight. Then, the weighted-median estimate was considered as the 50% weighted percentile. If > 50% of the weight was from a valid SNP, the weighted-median method was considered to produce an unbiased estimate of the causal effect of MR based on Instrument Strength Independent of Direct Effect (InSIDE). This method not only improves the ability to detect causal effects and reduces type I errors than the MR–Egger regression but also serves as a supplement of the MR–Egger method (Bowden et al., 2016), thereby providing more reliable effect estimation in the MR analysis. All analyses were performed through the “TwoSampleMR” package in the R software. The result of the statistical threshold p < 0.05 was considered to indicate causality.
Sensitivity Analysis
We used the leave-one-out analysis to assess whether the impact of a single SNP affected the correlation disproportionately and had an obvious horizontal pleiotropic effect (Emdin et al., 2017). We left out each SNP in sequence via performing two-sample MR analysis to estimate the effect. The leave-one-out analysis eliminates SNPs in order and recalculates the MR estimate, which in the causal effect estimate helps to identify SNPs causing a dramatic change. The sensitivity of the SNP is reflected by the variation in the results before and after deleting each SNP.
Results
Single-Nucleotide Polymorphism Selection and Validation
In all, seven of 153 SNPs (rs28673647, rs28407036, rs3124753, rs4962155, rs176691, rs41297217, and rs3118665) robustly associated with ADAMTS13 level were acquired from Ma’s study based on the LD independence test (Supplementary Table 1). All seven SNPs were identified when we evaluated the possible causal relationship between outcomes including CHD, MI, AF, HF, and VTE and ADAMTS13 level, all of which showed a strong association with these outcomes. We also identified four SNPs (rs41314453, rs10456544, rs3118667, and rs139911703) that predicted ADAMTS13 activity (Supplementary Table 2). Three SNPs (rs41314453, rs10456544, and rs3118667) were robustly qualified in the causal association test between ADAMTS13 and CHD, MI, and HF, respectively. Because we could not identify any proxies for the SNP rs139911703, it was excluded. Besides, all four SNPs that we assessed exhibited a causal relation between ADAMTS13 and outcomes (AF, ischemic stroke, and VTE). The details of each SNP we used after harmonization are listed in Supplementary Tables 3–7.
In addition, we performed a Q test based on its Q parameter to assess heterogeneity (Supplementary Tables 8, 9). We used the MR–Egger regression and the appearance of the funnel plot to test whether the horizontal pleiotropy has a significant influence on the results. The MR–Egger regression showed that there was a pleiotropic effect influencing the verification of the causal relationship between ADAMTS13 level and ischemic stroke (p < 0.05) (Supplementary Figures 1, 2 and Supplementary Tables 10, 11).
Table 1 shows the main characteristics of the included SNP data and the level and activity of ADAMTS13. Supplementary Tables 8–11 show information related to horizontal pleiotropic analysis and heterogeneity testing, respectively.
Table 1. Characteristics of the SNP summary statistics for ADAMTS13 and outcomes.
Analysis by Two-Sample Mendelian Randomization
Based on the inverse-variance weighted (IVW) analysis results, the ADAMTS13 level was not causally associated with CHD (b = −0.001, se = 0.001, p > 0.05), MI (b = −0.001, se = 0.001, p > 0.05), AF (b = 0.00003, se < 0.001, p > 0.05), HF (b = −0.00003, se < 0.001, p > 0.05), and VTE (b = −0.00003, se < 0.001, p > 0.05).
In addition, the casual association of ADAMTS13 activity has also been investigated. According to the IVW analysis results, ADAMTS13 activity was causally associated with CHD (b = −0.0041, se = 0.0019, p = 0.0366) and MI (b = −0.0048, se = 0.0022, p = 0.0266). The association between ADAMTS13 activity and AF (b = 0.00006, se < 0.001, p = 0.0986) could not be demonstrated by the standard IVW or by HF (b = −0.00004, se < 0.001, p = 0.0544), ischemic stroke (b = −0.00005, se < 0.001, p = 0.1893), and VTE (b = −0.00005, se < 0.001, p = 0.1546) based on UK Biobank.
The various statistical methods of MR analysis used to estimate the causal relationship between ADAMTS13 and the outcomes are summarized in Supplementary Tables 12, 13 and shown with a graph in Supplementary Figures 3, 4. The single SNP analysis is shown in Supplementary Figures 5, 6.
Sensitivity Analyses
Figure 2 and Supplementary Figure 7 depict the sensitivity analysis of each group using the leave-one-out method. In this analysis, we excluded each SNP in turn by recalculating the MR analysis to estimate the effect, thereby aiding in the identification of SNPs that led to a dramatic change in the MR analysis. Visually, the leave-one-out analysis chart showed that none of the SNPs exhibited a significant deviation from the result. Supplementary Tables 14, 15 showed the primary data of the leave-one-out analysis.
Figure 2. Leave-one-out analysis of ADAMTS13 level. (A) Coronary heart disease, (B) myocardial infarction, (C) atrial fibrillation, (D) heart failure, and (E) venous thromboembolism.
Discussion
We detected that ADAMTS13 level was not associated with CHD, MI, AF, HF, and VTE by performing a two-sample MR analysis via GWAS statistics related to ADAMTS13 and cardiovascular diseases. Meanwhile, a lower ADAMTS13 activity was causally associated with CHD and MI, while ADAMTS13 activity was not causally related to AF, HF, ischemic stroke, and VTE. Our results are similar to those reported in some large-scale observational studies in general, supporting the assumption on the causal effect of ADAMTS13 on CHD and MI by eliminating the potential bias caused by confounding factors in a genetic method and casting doubts on studies with a different conclusion (Chion et al., 2007). This study was verified via the two-sample MR study using different IVs for the ADAMTS13 level and activity identified by genome-wide significance in a large GWAS.
Over the past decade, the role of ADAMTS13 in CHD was controversial, and whether a causal relationship exists between them remains to be determined. Both in vitro and animal studies have shown that ADAMTS13 has antithrombotic properties (Chauhan et al., 2006), and reduced ADAMTS13 level and activity can increase the risk of cardiovascular disease. One study has shown that young patients with low ADAMTS13 levels are five times more likely to develop cardiovascular disease than those with normal ADAMTS13 levels (Marenberg et al., 1994). A prospective cohort study showed that low ADAMTS13 activity is linked to increased risk of CHD in the senior group (Sonneveld et al., 2016). In addition, a genetic study showed a significant correlation between the ADAMTS13 900 V variant and cardiac death (Schettert et al., 2010). In a meta-analysis of observational studies, Maino et al. (2015) found a relationship between lower ADAMTS13 levels and increased MI risk. However, the result of a large clinical study (Chion et al., 2007) showed that the risk of MI increased with ADAMTS13 level. Nevertheless, the patients in Chion’s report were all male. In this case, sex-based differences and discrimination between ADAMTS13 level and activity should be taken into consideration. A recent MR study focused on the association between ADAMTS13 activity and ischemic heart disease (Schooling et al., 2018), suggesting that ADAMTS13’s role may be sex-specific owing to lower ADAMTS13 activity in male than female patients, which is the same as that reported in another study (Kokame et al., 2011). However, the association between ADAMTS13 and other cardiovascular diseases such as AF, HF, ischemic stroke, and VTE has been poorly reported. Our MR study showed the causal relationship between low ADAMTS13 level and CHD from a genetic perspective, as well as ADAMTS13 activity.
The level of ADAMTS13 and its activity are distinct concepts worth differentiating between. The activity of ADAMTS13 depends mainly on its concentration in circulation and its ability to interact with VWF (Ma et al., 2017), so its plasma level is not equal to its activity, although it is an active form of ADAMTS13 secreted into the circulation. The regulation of ADAMTS13 activity is mainly allosteric and depends on specific explant–substrate interactions (Muia et al., 2014). The VWF levels are determined by multiple environmental and genetic factors and vary widely. However, the changes in ADAMTS13 activity appeared to be smaller than the changes in VWF concentration in healthy populations (Kokame et al., 2011; Desch et al., 2013). A recent study reported a GWAS of ADAMTS13 level (Ma et al., 2017), while de Vries et al. (2015) uncovered an association of single-nucleotide variants related with ADAMTS13 activity. Our study verified the causal relationship between ADAMTS13 level and activity and six types of cardiovascular disease using the MR method.
The ABO gene is located near the ADAMTS13 locus, which may serve as a confounding factor in this analysis. A previous meta-analysis demonstrated that there was no significant correlation between the ABO serotype marker allele and ADAMTS13 levels (Ma et al., 2017). Furthermore, some studies also stated that there was no association between ADAMTS13 activity and ABO blood groups (Chion et al., 2007). In our study, we utilized the MR-Egger regression method to estimate horizontal pleiotropy. No significant pleiotropic effect was evident in most cases. Therefore, we did not exclude the ABO loci from the IVs for the analysis (Sabater-Lleal et al., 2019).
The correlation of ADAMTS13 with cardiovascular diseases has been gradually recognized, although its role in the pathogenesis of cardiovascular diseases is inconclusive (Chion et al., 2007). It has been hypothesized that a significant increase in the VWF/ADAMTS13 ratios in coronary blood flow plays a pathogenic role in ACSs (Pedrazzini et al., 2016). It was also reported that the total VWF and VWF/ADAMTS13 ratios increased in patients with the first ST-segment elevated MI (STEMI), while plasma ADAMTS13 levels decreased (Rutten et al., 2015). VWF is known to be a large plasma glycoprotein that binds to platelets in a variety of conditions to mediate platelet adhesion and aggregation that has a crucial and direct relationship with ADAMTS13. Specifically, ADAMTS13 highly regulates the protein size determining the VWF activity, and the degradation of VWF by ADAMTS13 can decrease platelet activation (Zheng et al., 2001). Researchers have expressed concern about whether ADAMTS13 and VWF had a co-effect on cardiovascular disease. After adjusting for risk factors such as high-density lipoprotein cholesterol (Lowe et al., 2003), Lowe et al. pointed out that both ADAMTS13 and VWF were risk factors for non-fatal severe MI, with ADAMTS13 negatively and VWF positively associated with MI risk. Although decreased ADAMTS13 activity and increased VWF level increase the risk of atherosclerosis (Andersson et al., 2012), most scholars believe that ADAMTS13 and VWF level were two independent risk factors. A prospective cohort study showed that low ADAMTS13 activity, independent of VWF, is associated with increased CHD risk in the senior group (Sonneveld et al., 2016). Furthermore, genetic analysis pointed out a genetic polymorphism of VWF (V1565L), which modulates ADAMTS13 activity. The ADAMTS13 haplotype QAGA or H4 had an independent protective effect on coronary artery disease (CAD) (Lasom et al., 2017). Thus, the role of VWF in cardiovascular disease needs further investigation.
It is logical to predict that low-level ADAMTS13 is strongly associated with ischemic stroke. ADAMTS13 plays a vital part in the ischemic stroke pathogenesis confirmed by animal and clinical studies, indicating that ADAMTS13 may be a potential therapeutic target (Bongers et al., 2006; Lee et al., 2012). Most studies have linked ADAMTS13 levels to the risk of ischemic stroke in the general population. A recent study reported an association between low ADAMTS13 levels and high stroke incidence (Lambers et al., 2013). The level of ADAMTS13 in AIS patients was significantly lower than that in normal people and patients with chronic cardiovascular disease (Denorme et al., 2017). Furthermore, in a pediatric cohort study, researchers regarded reduced ADAMTS13 activity as a risk factor for pediatric ischemic stroke (Lambers et al., 2013). However, an animal model study showed that a lack of ADAMTS13 in mice led to a pre-thrombotic state instead of stroke (Banno et al., 2006). In summary, we used the widely accepted IVs for ADAMTS13 and the summary data from GWAS reports to demonstrate that ADAMTS13 activity was not associated with ischemic stroke risk in the general population.
At present, less attention has been paid to the correlation between ADAMTS13 and other cardiovascular diseases such as AF, VTE, and HF. A recent Framingham Heart Study showed that ADAMTS13 remained significantly associated with the risk of AF after the adjustments (Ko et al., 2019). In another study, scientists found that ADAMTS13 remained an independent predictor of recurrent AF (Freynhofer et al., 2011). A clinical study pointed out that ADAMTS13 can regulate the VWF thrombogenic activity, causing thromboembolism in AF (Ammash et al., 2011). However, a recent study found that the VWF level was an independent predictor of non-valvular AF instead of the ADAMTS13 level (Wysokinski et al., 2021). Our results suggested that both ADAMTS13 level and activity are not correlated with the occurrence of AF.
Many studies have reported an imbalance between VWF and ADAMTS13 in patients with arterial thrombosis, while several studies have also suggested a pathogenic role for ADAMTS13 in VTE. To be specific, as one clinical study pointed out, ADAMTS13 activity increased in patients with VTE (Mazetto et al., 2012). A recent study demonstrated that ADAMTS13 activity was significantly associated with VTE onset and survival in COVID-19 patients (Delrue et al., 2021). Regarding HF, a study showed that chronic HF (CHF) was associated with the dysregulation of endothelial activation. As ADAMTS13 can integrate the impaired liver synthesis ability and endothelial regulation, its abnormal activity can be used as an independent predictor of CHF (Gombos et al., 2009). Our results demonstrated there was no direct causal link between ADAMTS13 levels and VTE, as well as HF, and there was also no causal effect between ADAMTS13 activity and these outcomes. However, many researchers cast doubts regarding the basic assumption in a MR analysis that if a p < 0.05, it suggests that an event is very unlikely to occur. They thought that sample sizes play an indispensable part wherein instead of stating p < 0.05 or > 0.05, the actual calculated p-values and corresponding sample sizes should be disclosed by researchers (Laber and Shedden, 2017). Although MR results cannot demonstrate the causal relationship of ADAMTS13, we believe that further investigations are required to prove this.
There are several crucial strengths in this MR analysis. First, in this study, we used summary-level data from GWAS to systematically identify the association of risk factors with several cardiovascular diseases. Second, we used high statistical powers to detect the causal association of risk factors and ADAMTS13 with cardiovascular disease by using the MR approach. The two-sample MR analysis minimized the distortion by confounding factors and reverse causality, which is a major limitation of traditional observational studies. Finally, as most of the study populations were of European descent, our results were more accurate and reasonable and reduced population stratification bias.
Our study also has some limitations. First, all genetic associations come from studies conducted on European individuals whose data are not representative on a global scale. Second, we have not yet fully understood the underlying mechanism of exposure traits affecting the genetic effects of the cardiovascular system. Third, we noted that our MR–Egger results were less compelling than the result of IVW and median-based methods. However, MR–Egger is usually only considered a sensitivity analysis to estimate the validity of MR results. Fourth, the number of genetic indicators of ADAMTS13 activity may be limited, and some crucial genetic predictors of ADAMTS13 level lack external replication such as SUTPH, which could result in false positives in the research.
In summary, our MR results showed (1) the causal association of higher ADAMTS13 activity with reduced CHD and MI risk and (2) that no causal association existed of ADAMTS13 level with five kinds of cardiovascular diseases. Our results help to identify the role of ADAMTS13, including the level and activity, in cardiovascular disease via a genetic approach. However, our findings need further validation, and potential mechanisms are further needed.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
ZY contributed to data collection, manuscript writing, graphic mapping, and data processing. JZ contributed to the research design and data proofreading. Both authors contributed to the article and approved the submitted version
Funding
This work was supported by the National Natural Science Foundation of China (91639110), National Key Clinical Specialty Construction Project (2020-QTL-009), and Beijing Natural Science Foundation (7172195).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.660989/full#supplementary-material
Supplementary Figure 1 | Funnel plot visualizing the horizontal pleiotropy of SNPs related to ADAMTS13 level. (A) Coronary heart disease, (B) myocardial infarction, (C) atrial fibrillation, (D) heart failure, and (E) venous thromboembolism.
Supplementary Figure 2 | Funnel plot visualizing the horizontal pleiotropy of SNPs related to ADAMTS13 activity. (A) Coronary heart disease, (B) myocardial infarction, (C) atrial fibrillation, (D) heart failure, (E) ischemic stroke, and (F) venous thromboembolism.
Supplementary Figure 3 | Comparison of different statistical methods for MR analysis to assess the causality between ADAMTS13 level and cardiovascular diseases. ADAMTS13 level has a causally negative correlation with coronary heart disease (A), myocardial infarction (B), heart failure (D), and venous thromboembolism (E), ADAMTS13 has a positive causal association with atrial fibrillation (C).
Supplementary Figure 4 | Comparison of different statistical methods for MR analysis to assess causality between ADAMTS13 activity and cardiovascular diseases. ADAMTS13 activity has a causally negative correlation with coronary heart disease (A) and myocardial infarction (B). No correlation was observed between ADAMTS13 activity and atrial fibrillation (C), heart failure (D), ischemic stroke (E), and venous thromboembolism (F).
Supplementary Figure 5 | Single SNP analysis of ADAMTS13 level.
Supplementary Figure 6 | Single SNP analysis of ADAMTS13 activity.
Supplementary Figure 7 | Leave-one-out analysis of ADAMTS13 activity. (A) Coronary heart disease, (B) myocardial infarction, (C) atrial fibrillation, (D) heart failure, (E) ischemic stroke, and (F) venous thromboembolism.
Footnotes
References
Ammash, N., Konik, E. A., McBane, R. D., Chen, D., Tange, J. I., Grill, D. E., et al. (2011). Left atrial blood stasis and Von Willebrand factor-ADAMTS13 homeostasis in atrial fibrillation. Arterioscler. Thromb .Vasc. Biol. 31, 2760–2766. doi: 10.1161/ATVBAHA.111.232991
Andersson, H. M., Siegerink, B., Luken, B. M., Crawley, J. T., Algra, A., Lane, D. A., et al. (2012). High VWF, low ADAMTS13, and oral contraceptives increase the risk of ischemic stroke and myocardial infarction in young women. Blood 119, 1555–1560. doi: 10.1182/blood-2011-09-380618
Banno, F., Kokame, K., Okuda, T., Honda, S., Miyata, S., Kato, H., et al. (2006). Complete deficiency in ADAMTS13 is prothrombotic, but it alone is not sufficient to cause thrombotic thrombocytopenic purpura. Blood 107, 3161–3166. doi: 10.1182/blood-2005-07-2765
Bongers, T. N., de Maat, M. P., van Goor, M. L., Bhagwanbali, V., van Vliet, H. H., Gomez Garcia, E. B., et al. (2006). High von Willebrand factor levels increase the risk of first ischemic stroke: influence of ADAMTS13, inflammation, and genetic variability. Stroke 37, 2672–2677. doi: 10.1161/01.STR.0000244767.39962.f7
Bowden, J., Davey Smith, G., and Burgess, S. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int. J. Epidemiol. 44, 512–525. doi: 10.1093/ije/dyv080
Bowden, J., Davey Smith, G., Haycock, P. C., and Burgess, S. (2016). Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet. Epidemiol. 40, 304–314. doi: 10.1002/gepi.21965
Burgess, S., Scott, R. A., Timpson, N. J., Smith, G. Davey, Thompson, S. G., and Epic-InterAct Consortium. (2015). Using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur. J. Epidemiol. 30, 543–552. doi: 10.1007/s10654-015-0011-z
Chauhan, A. K., Motto, D. G., Lamb, C. B., Bergmeier, W., Dockal, M., Plaimauer, B., et al. (2006). Systemic antithrombotic effects of ADAMTS13. J. Exp. Med. 203, 767–776. doi: 10.1084/jem.20051732
Chion, C. K., Doggen, C. J., Crawley, J. T., Lane, D. A., and Rosendaal, F. R. (2007). ADAMTS13 and von Willebrand factor and the risk of myocardial infarction in men. Blood 109, 1998–2000. doi: 10.1182/blood-2006-07-038166
Cines, D. B., and Levine, L. D. (2017). Thrombocytopenia in pregnancy. Hematology Am. Soc. Hematol. Educ. Program 2017, 144–151. doi: 10.1182/asheducation-2017.1.144
Crawley, J. T., de Groot, R., Xiang, Y., Luken, B. M., and Lane, D. A. (2011). Unraveling the scissile bond: how ADAMTS13 recognizes and cleaves von Willebrand factor. Blood 118, 3212–3221. doi: 10.1182/blood-2011-02-306597
Crawley, J. T., Lane, D. A., Woodward, M., Rumley, A., and Lowe, G. D. (2008). Evidence that high von Willebrand factor and low ADAMTS-13 levels independently increase the risk of a non-fatal heart attack. J. Thromb. Haemost. 6, 583–588. doi: 10.1111/j.1538-7836.2008.02902.x
Danesh, J., Wheeler, J. G., Hirschfield, G. M., Eda, S., Eiriksdottir, G., Rumley, A., et al. (2004). C-reactive protein and other circulating markers of inflammation in the prediction of coronary heart disease. N. Engl. J. Med. 350, 1387–1397. doi: 10.1056/NEJMoa032804
Davies, N. M., Holmes, M. V., and Davey Smith, G. (2018). Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ 362:k601. doi: 10.1136/bmj.k601
de Vries, P. S., Boender, J., Sonneveld, M. A., Rivadeneira, F., Ikram, M. A., Rottensteiner, H., et al. (2015). Genetic variants in the ADAMTS13 and SUPT3H genes are associated with ADAMTS13 activity. Blood 125, 3949–3955. doi: 10.1182/blood-2015-02-629865
Delrue, M., Siguret, V., Neuwirth, M., Joly, B., Beranger, N., Sene, D., et al. (2021). von Willebrand factor/ADAMTS13 axis and venous thromboembolism in moderate-to-severe COVID-19 patients. Br. J. Haematol. 192, 1097–1100. doi: 10.1111/bjh.17216
Denorme, F., Kraft, P., Pareyn, I., Drechsler, C., Deckmyn, H., Vanhoorelbeke, K., et al. (2017). Reduced ADAMTS13 levels in patients with acute and chronic cerebrovascular disease. PLoS One 12:e0179258. doi: 10.1371/journal.pone.0179258
Desch, K. C., Ozel, A. B., Siemieniak, D., Kalish, Y., Shavit, J. A., Thornburg, C. D., et al. (2013). Linkage analysis identifies a locus for plasma von Willebrand factor undetected by genome-wide association. Proc. Natl. Acad. Sci. U. S. A. 110, 588–593. doi: 10.1073/pnas.1219885110
Domingueti, C. P., Dusse, L. M., Carvalho, M. D., Gomes, K. B., and Fernandes, A. P. (2013). Hypercoagulability and cardiovascular disease in diabetic nephropathy. Clin. Chim. Acta 415, 279–285. doi: 10.1016/j.cca.2012.10.061
Egger, M., Davey Smith, G., Schneider, M., and Minder, C. (1997). Bias in meta-analysis detected by a simple, graphical test. BMJ 315, 629–634. doi: 10.1136/bmj.315.7109.629
Emdin, C. A., Khera, A. V., and Kathiresan, S. (2017). Mendelian Randomization. JAMA 318, 1925–1926. doi: 10.1001/jama.2017.17219
Freynhofer, M. K., Bruno, V., Jarai, R., Gruber, S., Hochtl, T., Brozovic, I., et al. (2011). Levels of von Willebrand factor and ADAMTS13 determine clinical outcome after cardioversion for atrial fibrillation. Thromb. Haemost. 105, 435–443. doi: 10.1160/TH10-09-0615
Gombos, T., Mako, V., Cervenak, L., Papassotiriou, J., Kunde, J., Harsfalvi, J., et al. (2009). Levels of von Willebrand factor antigen and von Willebrand factor cleaving protease (ADAMTS13) activity predict clinical events in chronic heart failure. Thromb. Haemost. 102, 573–580. doi: 10.1160/TH09-01-0036
Haegele, S., Fuxsteiner, J., Pereyra, D., Koeditz, C., Rumpf, B., Schuetz, C., et al. (2018). Elevated ADAMTS13 Activity is Associated with Poor Postoperative Outcome in Patients Undergoing Liver Resection. Sci. Rep. 8:16823. doi: 10.1038/s41598-018-34794-w
Hanson, E., Jood, K., Karlsson, S., Nilsson, S., Blomstrand, C., and Jern, C. (2011). Plasma levels of von Willebrand factor in the etiologic subtypes of ischemic stroke. J. Thromb. Haemost. 9, 275–281. doi: 10.1111/j.1538-7836.2010.04134.x
Hemani, G., Zheng, J., Wade, K. H., Laurin, C., Elsworth, B., Burgess, S., et al. (2016). MR-Base: a platform for systematic causal inference across the phenome using billions of genetic associations. Preprint.
Horii, M., Uemura, S., Uemura, M., Matsumoto, M., Ishizashi, H., Imagawa, K., et al. (2008). Acute myocardial infarction as a systemic prothrombotic condition evidenced by increased von Willebrand factor protein over ADAMTS13 activity in coronary and systemic circulation. Heart Vessels 23, 301–307. doi: 10.1007/s00380-008-1053-x
Ko, D., Benson, M. D., Ngo, D., Yang, Q., Larson, M. G., Wang, T. J., et al. (2019). Proteomics Profiling and Risk of New-Onset Atrial Fibrillation: framingham Heart Study. J. Am. Heart Assoc. 8:e010976. doi: 10.1161/JAHA.118.010976
Kokame, K., Sakata, T., Kokubo, Y., and Miyata, T. (2011). von Willebrand factor-to-ADAMTS13 ratio increases with age in a Japanese population. J. Thromb. Haemost. 9, 1426–1428. doi: 10.1111/j.1538-7836.2011.04333.x
Laber, E. B., and Shedden, K. (2017). Statistical Significance and the Dichotomization of Evidence: the Relevance of the ASA Statement on Statistical Significance and p-values for Statisticians. J. Am. Stat. Assoc. 112, 902–904. doi: 10.1080/01621459.2017.1311265
Lambers, M., Goldenberg, N. A., Kenet, G., Kirkham, F. J., Manner, D., Bernard, T., et al. (2013). Role of reduced ADAMTS13 in arterial ischemic stroke: a pediatric cohort study. Ann. Neurol. 73, 58–64. doi: 10.1002/ana.23735
Lasom, S., Komanasin, N., Settasatian, N., Settasatian, C., Kukongviriyapan, U., Intharapetch, P., et al. (2017). Protective effect of a disintegrin and metalloproteinase with a thrombospondin type 1 motif, member 13 haplotype on coronary artery disease. Blood Coagul. Fibrinolysis 28, 286–294. doi: 10.1097/MBC.0000000000000594
Lee, M., Rodansky, E. S., Smith, J. K., and Rodgers, G. M. (2012). ADAMTS13 promotes angiogenesis and modulates VEGF-induced angiogenesis. Microvasc. Res. 84, 109–115. doi: 10.1016/j.mvr.2012.05.004
Levy, G. G., Nichols, W. C., Lian, E. C., Foroud, T., McClintick, J. N., McGee, B. M., et al. (2001). Mutations in a member of the ADAMTS gene family cause thrombotic thrombocytopenic purpura. Nature 413, 488–494. doi: 10.1038/35097008
Lowe, G. D., Woodward, M., Rumley, A., Morrison, C. E., and Nieuwenhuizen, W. (2003). Associations of plasma fibrinogen assays, C-reactive protein and interleukin-6 with previous myocardial infarction. J. Thromb. Haemost. 1, 2312–2316. doi: 10.1046/j.1538-7836.2003.00467.x
Ma, Q., Jacobi, P. M., Emmer, B. T., Kretz, C. A., Ozel, A. B., McGee, B., et al. (2017). Genetic variants in ADAMTS13 as well as smoking are major determinants of plasma ADAMTS13 levels. Blood Adv. 1, 1037–1046. doi: 10.1182/bloodadvances.2017005629
Machiela, M. J., and Chanock, S. J. (2015). LDlink: a web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants. Bioinformatics 31, 3555–3557. doi: 10.1093/bioinformatics/btv402
Maino, A., Siegerink, B., Lotta, L. A., Crawley, J. T., le Cessie, S., Leebeek, F. W., et al. (2015). Plasma ADAMTS-13 levels and the risk of myocardial infarction: an individual patient data meta-analysis. J. Thromb. Haemost. 13, 1396–1404. doi: 10.1111/jth.13032
Marenberg, M. E., Risch, N., Berkman, L. F., Floderus, B., and de Faire, U. (1994). Genetic susceptibility to death from coronary heart disease in a study of twins. N. Engl. J. Med. 330, 1041–1046. doi: 10.1056/NEJM199404143301503
Mazetto, B. M., Orsi, F. L., Barnabe, A., De Paula, E. V., Flores-Nascimento, M. C., and Annichino-Bizzacchi, J. M. (2012). Increased ADAMTS13 activity in patients with venous thromboembolism. Thromb. Res. 130, 889–893. doi: 10.1016/j.thromres.2012.09.009
Muia, J., Zhu, J., Gupta, G., Haberichter, S. L., Friedman, K. D., Feys, H. B., et al. (2014). Allosteric activation of ADAMTS13 by von Willebrand factor. Proc. Natl. Acad. Sci. U. S. A. 111, 18584–18589. doi: 10.1073/pnas.1413282112
Nikpay, M., Goel, A., Won, H. H., Hall, L. M., Willenborg, C., Kanoni, S., et al. (2015). A comprehensive 1,000 Genomes-based genome-wide association meta-analysis of coronary artery disease. Nat. Genet. 47, 1121–1130. doi: 10.1038/ng.3396
Pedrazzini, G., Biasco, L., Sulzer, I., Anesini, A., Moccetti, T., Kremer Hovinga, J. A., et al. (2016). Acquired intracoronary ADAMTS13 deficiency and VWF retention at sites of critical coronary stenosis in patients with STEMI. Blood 127, 2934–2936. doi: 10.1182/blood-2015-12-688010
Peyvandi, F., Hollestelle, M. J., Palla, R., Merlini, P. A., Feys, H. B., Vanhoorelbeke, K., et al. (2010). Active platelet-binding conformation of plasma von Willebrand factor in young women with acute myocardial infarction. J. Thromb. Haemost. 8, 1653–1656. doi: 10.1111/j.1538-7836.2010.03896.x
Rutten, B., Maseri, A., Cianflone, D., Laricchia, A., Cristell, N. A., Durante, A., et al. (2015). Plasma levels of active Von Willebrand factor are increased in patients with first ST-segment elevation myocardial infarction: a multicenter and multiethnic study. Eur. Heart J. Acute Cardiovasc. Care 4, 64–74. doi: 10.1177/2048872614534388
Sabater-Lleal, M., Huffman, J. E., de Vries, P. S., Marten, J., Mastrangelo, M. A., Song, C., et al. (2019). Genome-Wide Association Transethnic Meta-Analyses Identifies Novel Associations Regulating Coagulation Factor VIII and von Willebrand Factor Plasma Levels. Circulation 139, 620–635. doi: 10.1161/CIRCULATIONAHA.118.034532
Schettert, I. T., Pereira, A. C., Lopes, N. H., Hueb, W. A., and Krieger, J. E. (2010). Association between ADAMTS13 polymorphisms and risk of cardiovascular events in chronic coronary disease. Thromb. Res. 125, 61–66. doi: 10.1016/j.thromres.2009.03.008
Schooling, C. M., Luo, S., and Johnson, G. (2018). ADAMTS-13 activity and ischemic heart disease: a Mendelian randomization study. J. Thromb. Haemost. 16, 2270–2275. doi: 10.1111/jth.14267
Smith, G. D., and Ebrahim, S. (2003). ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease?. Int. J. Epidemiol. 32, 1–22. doi: 10.1093/ije/dyg070
Sonneveld, M. A., Kavousi, M., Ikram, M. A., Hofman, A., Rueda Ochoa, O. L., Turecek, P. L., et al. (2016). Low ADAMTS-13 activity and the risk of coronary heart disease - a prospective cohort study: the Rotterdam Study. J. Thromb. Haemost. 14, 2114–2120. doi: 10.1111/jth.13479
Uemura, M., Fujimura, Y., Ko, S., Matsumoto, M., Nakajima, Y., and Fukui, H. (2010). Pivotal role of ADAMTS13 function in liver diseases. Int. J. Hematol. 91, 20–29. doi: 10.1007/s12185-009-0481-4
Wang, A., Duan, Q., Wu, J., Liu, X., and Sun, Z. (2016). The expression of ADAMTS13 in human microvascular endothelial cells. Blood Coagul. Fibrinolysis 27, 464–466. doi: 10.1097/MBC.0000000000000405
Wieberdink, R. G., van Schie, M. C., Koudstaal, P. J., Hofman, A., Witteman, J. C., de Maat, M. P., et al. (2010). High von Willebrand factor levels increase the risk of stroke: the Rotterdam study. Stroke 41, 2151–2156. doi: 10.1161/STROKEAHA.110.586289
Wysokinski, W. E., Melduni, R. M., Ammash, N. M., Vlazny, D. T., Konik, E., Saadiq, R. A., et al. (2021). Von Willebrand Factor and ADAMTS13 as Predictors of Adverse Outcomes in Patients With Nonvalvular Atrial Fibrillation. CJC Open 3, 318–326. doi: 10.1016/j.cjco.2020.10.018
Zheng, J., Baird, D., Borges, M. C., Bowden, J., Hemani, G., Haycock, P., et al. (2017). Recent Developments in Mendelian Randomization Studies. Curr. Epidemiol. Rep. 4, 330–345. doi: 10.1007/s40471-017-0128-6
Zheng, X., Chung, D., Takayama, T. K., Majerus, E. M., Sadler, J. E., and Fujikawa, K. (2001). Structure of von Willebrand factor-cleaving protease (ADAMTS13), a metalloprotease involved in thrombotic thrombocytopenic purpura. J. Biol. Chem. 276, 41059–41063. doi: 10.1074/jbc.C100515200
Keywords: ADAMTS13, cardiovascular disease, coronary heart disease, myocardial infarction, Mendelian randomization
Citation: Ye Z and Zheng J (2021) Verification of the Role of ADAMTS13 in the Cardiovascular Disease Using Two-Sample Mendelian Randomization. Front. Genet. 12:660989. doi: 10.3389/fgene.2021.660989
Received: 03 February 2021; Accepted: 25 May 2021;
Published: 01 July 2021.
Edited by:
Alain Stepanian, Assistance Publique Hopitaux de Paris, FranceReviewed by:
Stella Trompet, Leiden University Medical Center, NetherlandsAndrew Johnson, National Heart, Lung, and Blood Institute (NHLBI), United States
Copyright © 2021 Ye and Zheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jingang Zheng, bWRqaW5nYW5nemhlbmdAMTYzLmNvbQ==