Plasma protein and venous thromboembolism: prospective cohort and mendelian randomisation analyses

Exploring genetic associations in systemic lupus erythematosus through Mendelian randomization: implications for novel biomarkers and therapeutic targets

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease characterized by a significant health burden. There is an essential need for novel biomarkers and therapeutic targets to improve diagnosis and management. Mendelian randomization (MR) was applied to explore causal links between SLE and various biomarkers like immune cells, metabolites, and inflammatory cytokines using multiple databases. Initially, biomarkers significantly associated with SLE were identified. Bidirectional MR helped clarify these relationships, and a two-step mediation MR examined their effects on SLE risk. Intersection analysis was used to identify biomarkers with consistent effects across datasets. Four biomarkers were identified as having significant associations with SLE risk: 1-palmitoyl-2-arachidonoyl-GPI levels [odds ratio (OR), 1.379; 95% confidence interval (CI), 1.180 to 1.613; FDR, 0.046], IL-17A levels (OR, 2.197; 95% CI, 1.412 to 3.418; FDR, 0.044), N-acetyl-aspartyl-glutamate (NAAG) levels (OR, 0.882; 95% CI, 0.831 to 0.936; FDR, 0.030), and ribitol levels (OR, 0.743; 95% CI, 0.644 to 0.857; FDR, 0.012). Bidirectional MR showed an inverse effect of NAAG on IL-17A levels (OR, 0.978; 95% CI, 0.962 to 0.994; p = 0.006). Mediation analysis indicated that NAAG influenced SLE risk both directly (beta =  - 0.108) and indirectly through IL-17A (beta =  - 0.018), highlighting the potential mediating role of IL-17A. After expanding the significance criteria to p < 0.05, intersection analysis across multiple datasets revealed 29 biomarkers with consistent beta directions, including 19 potential risk factors (beta > 0) and 10 protective factors (beta < 0) for SLE. This research has revealed significant genetic associations with SLE and demonstrated that IL-17A mediates the relationship between NAAG levels and SLE risk, highlighting potential new targets for personalized therapeutic interventions. Key Points • This study employs MR to identify significant genetic associations between various biomarkers and SLE, providing novel insights into potential biomarkers and therapeutic targets. • Four key biomarkers were identified as significantly associated with SLE risk: 1-palmitoyl-2-arachidonoyl-GPI, IL-17A, N-acetyl-aspartyl-glutamate (NAAG), and ribitol. • The findings suggest that NAAG levels have a protective effect against SLE, partly mediated through IL-17A, indicating a complex interplay between these biomarkers in the pathogenesis of SLE. • Intersectional analysis across multiple datasets revealed 29 biomarkers with consistent effects on SLE risk, highlighting new directions for future research and potential personalized therapeutic strategies.

Causal role of blood metabolites in HER-positive and HER-negative breast cancer: a Mendelian randomization (MR) study

BACKGROUND

Previous studies provide evidence that in vivo metabolites are associated with breast cancer (BC). However, the causal relationship between blood metabolites and BC remains unclear.

METHOD

Comprehensive two-sample Mendelian randomization analysis was conducted to determine the causal association between 1400 publicly available genetic data on metabolic factors and human epidermal growth factor receptor positive (HER+) BC or HER- BC in this study.

RESULT

Epiandrosterone sulfate levels (OR = 1.07, 95% CI = 1.02 ~ 1.10, p = 0.0013), 5alpha-androstan-3beta,17beta-diol monosulfate (2) levels (OR = 1.07, 95% CI = 1.03 ~ 1.12, p = 0.0012), glycohyocholate levels (OR = 0.85, 95% CI = 0.77 ~ 0.93, p = 0.0007) and etiocholanolone glucuronide levels (OR = 1.12, 95% CI = 1.05 ~ 1.20, p = 0.0013) were causally correlated with HER+ BC. 5 metabolites were causally correlated with HER- BC: Vanillic acid glycine levels (OR = 1.14, 95% CI = 1.06 ~ 1.22, p = 0.0003), Thyroxine levels (OR = 1.26, 95% CI = 1.11 ~ 1.44, p = 0.0004), 1-palmitoyl-2-linoleoyl-GPI (16:0/18:2) levels (OR = 0.86, 95% CI = 0.79 ~ 0.94, p = 0.0010), N-acetylphenylalanine levels (OR = 1.12, 95% CI = 1.05 ~ 1.19, p = 0.0007) and Glucose-to-mannose ratio (OR = 1.15, 95% CI = 1.06 ~ 1.24, p = 0.0008). Two common causally related metabolites were identified: Gamma-glutamyl glutamate and X-12849 levels.

CONCLUSIONS

Our study has respectively demonstrated the connection between blood metabolites and HER+ or HER- BC by genetic means, thereby offering opportunities for therapeutic targets.

Stress-Induced Hypercoagulability: Insights from Epidemiological and Mechanistic Studies, and Clinical Integration

By integrating findings from comprehensive reviews, meta-analyses, and cutting-edge genetic studies, this article illuminates the significance of stress-induced hypercoagulability in clinical medicine. In particular, the findings from numerous prospective cohort studies indicate that stress and hemostatic factors of a hypercoagulable state are associated with increased incident risk and poor prognosis for atherosclerotic cardiovascular disease and venous thromboembolism. Mendelian randomization studies suggest that these associations are partially causal. The review synthesizes extensive research on the link between acute and chronic stress and hypercoagulability, outlining a potential pathway from stress to thrombosis risk. Consistent with the allostatic load concept, acute stress-induced hypercoagulability, initially adaptive, can turn maladaptive under chronic stress or excessive acute stress, leading to arterial or venous thrombotic events. Individuals with predisposing factors, including atherosclerosis, thrombophilia, or immobilization, may exhibit an increased risk of thrombotic disease during stress. Contextual sociodemographic characteristics, the stress experience, and coping resources additionally modulate the extent of stress-induced hypercoagulability. Research into the neuroendocrine, cellular, and molecular bases reveals how stress influences platelet activation coagulation and fibrinolysis. The activation of the sympathetic nervous system and the hypothalamic-pituitary-adrenal axis, along with vagal withdrawal, and the effects of catecholamines, cortisol, and vasopressin, are the central mechanisms involved. Hemoconcentration, inflammation, endothelial dysfunction, and thrombopoiesis additionally contribute to stress-induced hypercoagulability. Further research is needed to prove a causal link between chronic stress and hypercoagulability. This includes exploring its implications for the prevention and management of thrombotic diseases in stressed individuals, with a focus on developing effective psychosocial and pharmacological interventions.

Venous thromboembolic disease genetics: from variants to function

Venous thromboembolic disease (VTE) is a prevalent and potentially life-threatening vascular disease, including both deep vein thrombosis and pulmonary embolism. This review will focus on recent insights into the heritable factors that influence an individual's risk for VTE. Here, we will explore not only the discovery of new genetic risk variants but also the importance of functional characterization of these variants. These genome-wide studies should lead to a better understanding of the biological role of genes inside and outside of the canonical coagulation system in thrombus formation and lead to an improved ability to predict an individual's risk of VTE. Further understanding of the molecular mechanisms altered by genetic variation in VTE risk will be accelerated by further human genome sequencing efforts and the use of functional genetic screens.

Targeted proteomics involved in cardiovascular health and heart rate variability in children with overweight/obesity

BACKGROUND

Children with overweight/obesity often exhibit alterations in their plasma protein profiles and reduced heart rate variability (HRV). Plasma proteomics is at the forefront of identifying biomarkers for various clinical conditions. We aimed to examine the association between plasma-targeted proteomics involved in cardiovascular health and resting vagal-related HRV parameters in children with overweight/obesity.

METHODS

Forty-four children with overweight/obesity (10.2 ± 1.1 years old; 52% boys) participated in the study. Olink's technology was used to quantify 92 proteins involved in cardiovascular health. HRV was measured using a heart rate monitor (Polar RS800CX). Four resting vagal-related HRV parameters were derived in time- and frequency-domain.

RESULTS

Eight proteins (KIM1, IgG Fc receptor II-b, IDUA, BOC, IL1RL2, TNFRSF11A, VSIG2, and TF) were associated with at least one out of the four vagal-related HRV parameters (β values ranging from -0.188 to 0.288; all p < .05), while KIM1, IDUA, and BOC associated with ≥ three vagal-related HRV parameters. Multiple hypothesis testing corrections did not reach statistical significance (false discovery rate [FDR >0.05]).

CONCLUSION

Plasma-targeted proteomics suggested novel biomarkers for resting vagal-related HRV parameters in children with overweight/obesity. Future studies using larger cohorts and longitudinal designs should confirm our findings and their potential clinical implications.

Genetic assessment of efficacy and safety profiles of coagulation cascade proteins identifies Factors II and XI as actionable anticoagulant targets

Aims

Anticoagulants are routinely used by millions of patients worldwide to prevent blood clots. Yet, problems with anticoagulant therapy remain, including a persistent and cumulative bleeding risk in patients undergoing prolonged anticoagulation. New safer anticoagulant targets are needed.

Methods and results

To prioritize anticoagulant targets with the strongest efficacy [venous thromboembolism (VTE) prevention] and safety (low bleeding risk) profiles, we performed two-sample Mendelian randomization and genetic colocalization. We leveraged three large-scale plasma protein data sets (deCODE as discovery data set and Fenland and Atherosclerosis Risk in Communities as replication data sets] and one liver gene expression data set (Institut Universitaire de Cardiologie et de Pneumologie de Québec bariatric biobank) to evaluate evidence for a causal effect of 26 coagulation cascade proteins on VTE from a new genome-wide association meta-analysis of 44 232 VTE cases and 847 152 controls, stroke subtypes, bleeding outcomes, and parental lifespan as an overall measure of efficacy/safety ratio. A 1 SD genetically predicted reduction in F2 blood levels was associated with lower risk of VTE [odds ratio (OR) = 0.44, 95% confidence interval (CI) = 0.38-0.51, P = 2.6e-28] and cardioembolic stroke risk (OR = 0.55, 95% CI = 0.39-0.76, P = 4.2e-04) but not with bleeding (OR = 1.13, 95% CI = 0.93-1.36, P = 2.2e-01). Genetically predicted F11 reduction was associated with lower risk of VTE (OR = 0.61, 95% CI = 0.58-0.64, P = 4.1e-85) and cardioembolic stroke (OR = 0.77, 95% CI = 0.69-0.86, P = 4.1e-06) but not with bleeding (OR = 1.01, 95% CI = 0.95-1.08, P = 7.5e-01). These Mendelian randomization associations were concordant across the three blood protein data sets and the hepatic gene expression data set as well as colocalization analyses.

Conclusion

These results provide strong genetic evidence that F2 and F11 may represent safe and efficacious therapeutic targets to prevent VTE and cardioembolic strokes without substantially increasing bleeding risk.

Perioperative Treatment with Rivaroxaban and Dabigatran on Changes of Coagulation and Platelet Activation Biomarkers following Left Atrial Appendage Closure

Insufficient data exist regarding the investigation of the impact of novel oral anticoagulants (NOACs) on coagulation activation biomarkers in the context of left atrial appendage closure (LAAC) and device-related thrombosis (DRT). The study was designed to investigate the changes and presence of coagulation activation biomarkers between different antithrombotic strategies following LAAC. A total of 120 nonvalvular atrial fibrillation patients intolerant of long-term anticoagulants, who underwent successful WATCHMAN closure implantation, were enrolled (rivaroxaban, n = 82; dabigatran, n = 38). Blood samples were obtained from left atrium (LA) and left atrial appendage (LAA) during the operation and fasting blood samples on the same day of LAAC and 45 days after discharge. The biochemical indicators, thrombin-antithrombin complex (TAT), soluble P-selectin (sP-selectin), von Willebrand factor (vWF), and CD40 ligand (CD40L), were measured by enzyme-linked immunosorbent assay. The primary endpoints of this study were the efficacy and safety characteristics of different antithrombotic strategies, including DRT incidence, stroke or transient ischemic attack, systemic embolism, and clinical major and nonmajor bleeding complications during the follow-up of 180 days. The results revealed that TAT, vWF, sP-selectin, and CD40L levels in vein were significantly reduced by 2.4% (p = 0.043), 5.0% (p < 0.001), 8.7% (p < 0.001), and 2.5% (p = 0.043) from their baseline levels after rivaroxaban treatment. Conversely, no significant changes were detected in the dabigatran group. Furthermore, the plasma levels of platelet activation biomarkers (CD40L and sP-selectin) in both LA and LAA groups were significantly lower after anticoagulation with rivaroxaban, as compared to dabigatran treatment (CD40L: 554.62 ± 155.54 vs. 445.02 ± 130.04 for LA p = 0.0013, 578.51 ± 156.28 vs. 480.13 ± 164.37 for LAA p = 0.0052; sP-selectin: 2849.07 ± 846.69 vs. 2225.54 ± 799.96 for LA p = 0.0105, 2915.52 ± 1402.40 vs. 2203.41 ± 1061.67 for LAA p = 0.0022). Notably, the present study suggests that rivaroxaban may be more effective in the prevention of DRT for patients undergoing LAAC.

Proteomic insights into modifiable risk of venous thromboembolism and cardiovascular comorbidities

BACKGROUND

Venous thromboembolism (VTE) has been associated with several modifiable factors (MFs) and cardiovascular comorbidities. However, the mechanisms are largely unknown.

OBJECTIVES

We aimed to decipher proteomic pathways underlying the associations of VTE with MFs and cardiovascular comorbidities.

METHODS

A 2-stage network Mendelian randomization analysis was conducted to explore the associations between 15 MFs, 1151 blood proteins, and VTE using data from a genome-wide meta-analysis including 81 190 cases of VTE. We used protein data from 35 559 individuals as the discovery analysis, and from 2 independent studies including 10 708 and 54 219 participants as the replication analyses. Based on the identified proteins, we assessed the druggability and examined the cardiovascular pleiotropy.

RESULTS

The network Mendelian randomization analyses identified 10 MF-VTE, 86 MF-protein, and 34 protein-VTE associations. These associations were overall consistent in the replication analyses. Thirty-eight pathways with directionally consistent direct and indirect effects in the MF-protein-VTE pathway were identified. Low-density lipoprotein receptor-related protein 12 (LRP12: 34.3%-58.1%) and coagulation factor (F)XI (20.6%-39.6%) mediated most of the associations between 3 obesity indicators and VTE. Likewise, coagulation FXI mediated most of the smoking-VTE association (40%; 95% CI, 20%-60%) and insomnia-VTE association (27%; 95% CI, 5%-49%). Many VTE-associated proteins were highly druggable for thrombotic conditions. Five proteins (interleukin-6 receptor subunit alpha, LRP12, prothrombin, angiopoietin-1, and low-density lipoprotein receptor-related protein 4) were associated with VTE and its cardiovascular comorbidities.

CONCLUSION

This study suggests that coagulation FXI, a druggable target, is an important mediator of the associations of obesity, smoking, and insomnia with VTE risk.

Proteome-wide mendelian randomization identifies causal plasma proteins in venous thromboembolism development

Genome-wide association studies (GWAS) have identified numerous risk loci for venous thromboembolism (VTE), but it is challenging to decipher the underlying mechanisms. We employed an integrative analytical pipeline to transform genetic associations to identify novel plasma proteins for VTE. Proteome-wide association studies (PWAS) were determined by functional summary-based imputation leveraging data from a genome-wide association analysis (14,429 VTE patients, 267,037 controls), blood proteomes (1348 cases), followed by Mendelian randomization, Bayesian colocalization, protein-protein interaction, and pathway enrichment analysis. Twenty genetically regulated circulating protein abundances (F2, F11, ABO, PLCG2, LRP4, PLEK, KLKB1, PROC, KNG1, THBS2, SERPINA1, RARRES2, CEL, GP6, SERPINE2, SERPINA10, OBP2B, EFEMP1, F5, and MSR1) were associated with VTE. Of these 13 proteins demonstrated Mendelian randomized correlations. Six proteins (F2, F11, PLEK, SERPINA1, RARRES2, and SERPINE2) had strong support in colocalization analysis. Utilizing multidimensional data, this study suggests PLEK, SERPINA1, and SERPINE2 as compelling proteins that may provide key hints for future research and possible diagnostic and therapeutic targets for VTE.

Dynamic observation of circRNA and mRNA profiles in a rat model of deep vein thrombosis

The goal of the present study was to identify different transcriptome expression profiles involved in the pathogenesis of deep vein thrombosis (DVT), and to illustrate the diagnostic and therapeutic potential of circular RNAs (circRNAs) and mRNAs in DVT progression. A Sprague-Dawley rat model of DVT was successfully established through the stenosis method and samples were sequenced at four time points (1, 6 and 12 h, and 3 days after ligation) to observe the dynamic changes in circRNAs and mRNAs during DVT progression. RNA sequencing was used to analyze the circRNA and mRNA expression profiles, and associated functions and pathways, in the blood of DVT rats at the four time points. In addition, Short Time Series Expression Miner (STEM) analysis was performed to explore temporal gene expression. Differential expression of 1,680, 4,018, 3,724, and 3,036 circRNAs, and 400, 1,176, 373, and 573 mRNAs was observed in the 1, 6 and 12 h, and 3-day groups, respectively, compared with the sham group (fold change >2.0 or <-2.0, P<0.05). Functional enrichment analysis indicated that differentially expressed mRNAs were associated with the following terms: Immune response, cell activation, blood stasis facilitated organelle, extracellular membrane-bounded organelle, and blood microparticle, oxygen transporter activity. STEM analysis indicated that the expression of 366 circRNAs in circRNA profile 45 and 270 mRNAs in mRNA profile 45 was consistent with thrombus progression. Enrichment analysis was performed on mRNA profile 45. The main Gene Ontology annotations were chromosome segregation, mitotic sister chromatid segregation, cell cycle process, and ligand-dependent nuclear receptor transcription coactivator activity. Pathway enrichment analysis identified the platelet-associated pathway, immune-associated pathway, and inflammation-relation pathway. According to the enriched platelet-associated pathways, four mRNAs and ten candidate circRNAs were selected for reverse transcription-quantitative PCR verification. The expression of nine of the ten circRNAs and all four mRNAs was consistent with the sequencing results. In summary, differentially expressed circRNAs and mRNAs are dynamically involved in DVT development. Dysregulated transcriptome profiles and the corresponding functions and pathways may provide mechanistic insights into DVT diagnosis and treatment.

Deciphering the genetic architecture of atrial fibrillation offers insights into disease prediction, pathophysiology and downstream sequelae

Aims

The study aimed to discover novel genetic loci for atrial fibrillation (AF), explore the shared genetic etiologies between AF and other cardiovascular and cardiometabolic traits, and uncover AF pathogenesis using Mendelian randomization analysis.

Methods and results

We conducted a genome-wide association study meta-analysis including 109,787 AF cases and 1,165,920 controls of European ancestry and identified 215 loci, among which 91 were novel. We performed Genomic Structural Equation Modeling analysis between AF and four cardiovascular comorbidities (coronary artery disease, ischemic stroke, heart failure, and vneous thromboembolism) and found 189 loci shared across these diseases as well as a universal genetic locus shared by atherosclerotic outcomes (i.e., rs1537373 near CDKN2B). Three genetic loci (rs10740129 near JMJD1C, rs2370982 near NRXN3, and rs9931494 near FTO) were associated with AF and cardiometabolic traits. A polygenic risk score derived from this genome-wide meta-analysis was associated with AF risk (odds ratio 2.36, 95% confidence interval 2.31-2.41 per standard deviation increase) in the UK biobank. This score, combined with age, sex, and basic clinical features, predicted AF risk (AUC 0.784, 95% CI 0.781-0.787) in Europeans. Phenome-wide association analysis of the polygenic risk score identified many AF-related comorbidities of the circulatory, endocrine, and respiratory systems. Phenome-wide and multi-omic Mendelian randomization analyses identified associations of blood lipids and pressure, diabetes, insomnia, obesity, short sleep, and smoking, 27 blood proteins, one gut microbe (genus.Catenibacterium), and 11 blood metabolites with risk to AF.

Conclusions

This genome-wide association study and trans-omic Mendelian randomization analysis provides insights into disease risk prediction, pathophysiology and downstream sequelae.