Relationship between cathepsins and cardiovascular diseases: a Mendelian randomized study.

Background: Cardiovascular diseases (CVDs) are the leading age-related disorders worldwide, with their prevalence increasing annually. Cathepsins are protein-degrading enzymes essential for processes such as intracellular protein breakdown, apoptosis, and immune responses. Recent studies suggest a potential link between cathepsins and CVDs, yet the exact causal relationship remains to be elucidated. To address this, we propose using Mendelian randomization (MR) to explore the causal relationships between cathepsins and CVDs. Methods: We obtained single nucleotide polymorphism (SNP) data for cathepsins from the INTERVAL study, a publicly accessible genome-wide association study (GWAS) dataset. Outcome SNP data were sourced from seven distinct GWAS datasets, ensuring a comprehensive analysis across multiple cardiovascular outcomes. For MR analysis, we primarily employed the inverse variance weighted (IVW) method, known for its efficiency when all SNPs are valid instruments. This was supplemented by the weighted median and MR-Egger methods to provide robustness against potential violations of MR assumptions, such as pleiotropy. The IVW method offers precision and efficiency, the weighted median method adds robustness against invalid instruments, and the MR-Egger method helps identify and correct for pleiotropic biases. Cochran's Q test was utilized to assess heterogeneity, and sensitivity analyses were conducted using MR-PRESSO and the leave-one-out approach. Results: The strength of the associations between exposure and outcome was measured using odds ratios (ORs), and results were presented with 95% confidence intervals (CIs). The cathepsin E increases the risk of myocardial infarction (MI) (OR = 1.053%, 95% CI: 1.007-1.101, p = 0.024) and ischemic stroke (IS) (OR = 1.06%, 95% CI: 1.019-1.103, p = 0.004). Conversely, cathepsin L2 decreases the risk of chronic heart failure (CHF) (OR = 0.922%, 95% CI: 0.859-0.99, p = 0.025) and atrial fibrillation (AF) (OR = 0.956%, 95% CI: 0.918-0.996, p = 0.033). Cathepsin O was associated with an increased risk of IS (OR = 1.054%, 95% CI: 1.008-1.102, p = 0.021) and AF (OR = 1.058%, 95% CI: 1.02-1.098, p = 0.002). Conclusion: Our MR analysis reveals that cathepsin E is a risk factor for MI and IS, cathepsin L2 offers protective effects against CHF and AF, and cathepsin O increases the risk for IS and AF.

Li Q ，Zhou Z ，Xu T ，Gao X ，Lou Y ，Chen Z ，Zhang M ，Fang Q ，Tan J ，Huang J ... -  《Frontiers in Pharmacology》

Genetic insights into the role of cathepsins in cardiovascular diseases: a Mendelian randomization study.

This study aimed to explore the causal relationships between cathepsins and cardiovascular diseases (CVDs) by Mendelian randomization (MR) analysis. Single nucleotide polymorphisms (SNPs) associated with nine cathepsin types (cathepsins B, E, F, G, H, O, S, L2, and Z) were obtained from the INTERVAL study (3301 individuals). CVDs data were acquired from the UK Biobank (coronary atherosclerosis: 14 334 cases, 346 860 controls) and a genome-wide association study (GWAS) (myocardial infarction: 20 917 cases, 440 906 controls; myocarditis: 633 cases, 427 278 controls; chronic heart failure: 14 262 cases, 471 898 controls; angina pectoris: 30 025 cases, 440 906 controls; stable angina pectoris: 17 894 cases, 325 132 controls; unstable angina pectoris: 9481 cases, 446 987 controls; pericarditis: 1795 cases, 453 370 controls). Inverse variance weighted (IVW), MR-Egger, weighted median methods were adopted to conduct univariable MR (UVMR), reverse MR, multivariable MR (MVMR) analyses to estimate causality. The UVMR analyses demonstrated significant causal relationships between higher cathepsin E levels and increased risk of coronary atherosclerosis [IVW: P = 0.0051, odds ratio (OR) = 1.0033, 95% confidence interval (CI) = 1.0010-1.0056] and myocardial infarction (IVW: P = 0.0097, OR = 1.0553, 95% CI = 1.0131-1.0993), while elevated cathepsin L2 levels were causally related to reduced risk of myocarditis (IVW: P = 0.0120, OR = 0.6895, 95% CI = 0.5158-0.9216) and chronic heart failure (IVW: P = 0.0134, OR = 0.9316, 95% CI = 0.8807-0.9854). Reverse MR analyses revealed that myocardial infarction increased cathepsin O levels (IVW: P = 0.0400, OR = 1.0708, 95% CI = 1.0031-1.1431). MVMR analyses treating nine cathepsins together revealed that the positive causality between cathepsin E levels and coronary atherosclerosis risk (IVW: P = 0.0390, OR = 1.0030, 95% CI = 1.0000-1.0060), and the protective effect of cathepsin L2 levels on myocarditis (IVW: P = 0.0030, OR = 0.6610, 95% CI = 0.5031-0.8676) and chronic heart failure (IVW: P = 0.0090, OR = 0.9259, 95% CI = 0.8737-0.9812) remained, as higher cathepsin O levels were found to be causally related to increased risks of myocarditis (IVW: P = 0.0030, OR = 1.6145, 95% CI = 1.1829-2.2034) and chronic heart failure (IVW: P = 0.0300, OR = 1.0779, 95% CI = 1.0070-1.1537). The study highlights the causalities of cathepsin E, L2, and O on CVDs, offering insights into their roles in cardiovascular biomarkers and therapeutic targets development. Further research is required to apply these genetic findings clinically.

Zeng R ，Zhou Z ，Liao W ，Guo B ... -  《ESC Heart Failure》

Identification of atrial fibrillation-related genes through transcriptome data analysis and Mendelian randomization.

Atrial fibrillation (AF) is a common persistent arrhythmia characterized by rapid and chaotic atrial electrical activity, potentially leading to severe complications such as thromboembolism, heart failure, and stroke, significantly affecting patient quality of life and safety. As the global population ages, the prevalence of AF is on the rise, placing considerable strains on individuals and healthcare systems. This study utilizes bioinformatics and Mendelian Randomization (MR) to analyze transcriptome data and genome-wide association study (GWAS) summary statistics, aiming to identify biomarkers causally associated with AF and explore their potential pathogenic pathways. We obtained AF microarray datasets GSE41177 and GSE79768 from the Gene Expression Omnibus (GEO) database, merged them, and corrected for batch effects to pinpoint differentially expressed genes (DEGs). We gathered exposure data from expression quantitative trait loci (eQTL) and outcome data from AF GWAS through the IEU Open GWAS database. We employed inverse variance weighting (IVW), MR-Egger, weighted median, and weighted model approaches for MR analysis to assess exposure-outcome causality. IVW was the primary method, supplemented by other techniques. The robustness of our results was evaluated using Cochran's Q test, MR-Egger intercept, MR-PRESSO, and leave-one-out sensitivity analysis. A "Veen" diagram visualized the overlap of DEGs with significant eQTL genes from MR analysis, referred to as common genes (CGs). Additional analyses, including Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and immune cell infiltration studies, were conducted on these intersecting genes to reveal their roles in AF pathogenesis. The combined dataset revealed 355 differentially expressed genes (DEGs), with 228 showing significant upregulation and 127 downregulated. Mendelian randomization (MR) analysis identified that the autocrine motility factor receptor (AMFR) [IVW: OR = 0.977; 95% CI, 0.956-0.998; P = 0.030], leucine aminopeptidase 3 (LAP3) [IVW: OR = 0.967; 95% CI, 0.934-0.997; P = 0.048], Rab acceptor 1 (RABAC1) [IVW: OR = 0.928; 95% CI, 0.875-0.985; P = 0.015], and tryptase beta 2 (TPSB2) [IVW: OR = 0.971; 95% CI, 0.943-0.999; P = 0.049] are associated with a reduced risk of atrial fibrillation (AF). Conversely, GTPase-activating SH3 domain-binding protein 2 (G3BP2) [IVW: OR = 1.030; 95% CI, 1.004-1.056; P = 0.024], integrin subunit beta 2 (ITGB2) [IVW: OR = 1.050; 95% CI, 1.017-1.084; P = 0.003], glutaminyl-peptide cyclotransferase (QPCT) [IVW: OR = 1.080; 95% CI, 1.010-0.997; P = 1.154], and tripartite motif containing 22 (TRIM22) [IVW: OR = 1.048; 95% CI, 1.003-1.095; P = 0.035] are positively associated with AF risk. Sensitivity analyses indicated a lack of heterogeneity or horizontal pleiotropy (P > 0.05), and leave-one-out analysis did not reveal any single nucleotide polymorphisms (SNPs) impacting the MR results significantly. GO and KEGG analyses showed that CG is involved in processes such as protein polyubiquitination, neutrophil degranulation, specific and tertiary granule formation, protein-macromolecule adaptor activity, molecular adaptor activity, and the SREBP signaling pathway, all significantly enriched. The analysis of immune cell infiltration demonstrated associations of CG with various immune cells, including plasma cells, CD8T cells, resting memory CD4T cells, regulatory T cells (Tregs), gamma delta T cells, activated NK cells, activated mast cells, and neutrophils. By integrating bioinformatics and MR approaches, genes such as AMFR, G3BP2, ITGB2, LAP3, QPCT, RABAC1, TPSB2, and TRIM22 are identified as causally linked to AF, enhancing our understanding of its molecular foundations. This strategy may facilitate the development of more precise biomarkers and therapeutic targets for AF diagnosis and treatment.

Zhang Y ，Lian Q ，Nie Y ，Zhao W ... -  《Frontiers in Cardiovascular Medicine》

[Genetic Causation Analysis of Hyperandrogenemia Testing Indicators and Preeclampsia].

Some epidemiological studies have shown that pregnant women who develop preeclampsia (PE) have elevated levels of testosterone in their maternal plasma compared to women with normal blood pressure during pregnancy, revealing a potential association between hyperandrogenism in women and PE. To explore the causal relationship between hyperandrogenism and PE, this study selected total testosterone (TT), bioavailable testosterone (BIOT), and sex hormone binding globulin (SHBG) as exposure factors and PE and chronic hypertension with superimposed PE as disease outcomes. Two-sample Mendelian randomization (MR) analyses were used to genetically dissect the causal relationships between the three exposure factors (TT, BIOT, and SHBG) and the outcomes of PE and chronic hypertension with superimposed PE. Two independent genome-wide association study (GWAS) databases were used for the two-sample MR analysis. In the GWAS data of female participants from the UK Biobank cohort, single nucleotide polymorphisms (SNPs) associated with TT, BIOT, and SHBG were analyzed, involving 230454, 188507, and 188908 samples, respectively. GWAS data on PE and chronic hypertension with superimposed PE from the Finnish database were used to calculate SNP, involving 3556 PE cases and 114735 controls, as well as 38 cases of chronic hypertension with superimposed PE and 114735 controls. To meet the assumptions of instrumental relevance and independence in MR analysis, SNPs associated with exposure were identified at the genome-wide level (P<5.0×10-8), and those in linkage disequilibrium interference were excluded based on clustering thresholds of R 2<0.001 and an allele distance greater than 10000 kb. Known confounding factors, including previous PE, chronic kidney disease, chronic hypertension, diabetes, systemic lupus erythematosus, or antiphospholipid syndrome, were also identified and the relevant SNPs were removed. Finally, we extracted the outcome data based on the exposure-related SNPs in the outcome GWAS, integrating exposure and outcome data, and removing palindromic sequences. Five genetic causal analysis methods, including inverse variance-weighted method (IVW), MR-Egger regression, weighted median method, simple mode method, and weighted mode method, were used to infer causal relationships. In the IVW, it was assumed that the selected SNPs satisfied the three assumptions and provided the most ideal estimate of the effect. IVW was consequently used as the primary analysis method in this study. Considering the potential heterogeneity among the instrumental variables, random-effects IVW was used for MR analysis. The results were interpreted using odds ratios (OR) and the corresponding 95% confidence interval (CI) to explain the impact of exposure factors on PE and chronic hypertension with superimposed PE. If the CI did not include 1 and had a P value less than 0.05, the difference was considered statistically significant. Sensitivity analysis was conducted to assess heterogeneity and pleiotropy. Heterogeneity was examined using Cochran's Q test, and pleiotropy was assessed using MR-Egger intercept analysis. Additionally, leave-one-out analysis was conducted to examine whether individual SNPs were driving the causal associations. To further validate the findings, MR analyses were performed using the same methods and outcome variables, but with different exposure factors, including waist-to-hip ratio adjusted for BMI (WHRadjBMI) and 25-hydroxyvitamin D levels, with MR results for WHRadjBMI and PE serving as the positive controls and MR results for 25-hydroxyvitamin D levels and PE as the negative controls. According to the criteria for selecting genetic instrumental variables, 186, 127, and 262 SNPs were identified as genetic instrumental variables significantly associated with testosterone indicators TT, BIOT, and SHBG. MR analysis did not find a causal relationship between the TT, BIOT, and SHBG levels and the risk of developing PE and chronic hypertension with superimposed PE. The IVW method predicted that genetically predicted TT (OR [95% CI]=1.018 [0.897-1.156], P=0.78), BIOT (OR [95% CI]=1.11 [0.874-1.408], P=0.392), and SHBG (OR [95% CI]=0.855 [0.659-1.109], P=0.239) were not associated with PE. Similarly, genetically predicted TT (OR [95% CI]=1.222 [0.548-2.722], P=0.624), BIOT (OR [95% CI]=1.066 [0.242-4.695], P=0.933), and SHBG (OR [95% CI]=0.529 [0.119-2.343], P=0.402) were not significantly associated with chronic hypertension with superimposed PE. Additionally, MR analysis using the MR-Egger method, weighted median method, simple mode method, and weighted mode method yielded consistent results, indicating no significant causal relationship between elevated testosterone levels and PE or chronic hypertension with superimposed PE. Heterogeneity was observed for SHBG in the analysis with PE (Cochran's Q test, P=0.01), and pleiotropy was detected for BIOT in the analysis with PE (MR-Egger intercept analysis, P=0.014), suggesting that the instrumental variables did not affect PE through BIOT. Other instrumental variables did not show significant heterogeneity or pleiotropy. Leave-one-out analysis confirmed that the results of the MR analysis were not driven by individual instrumental variables. Consistent with previous MR studies, the results of the control MR analyses using WHRadjBMI and 25-hydroxyvitamin D levels supported the accuracy of the MR analysis approach and the methods used in this study. The MR analysis results suggest that current genetic evidence does not support a causal relationship between TT, BIOT, and SHBG levels and the development of PE and chronic hypertension with superimposed PE. This study suggests that elevated testosterone may be a risk factor for PE but not a direct cause.

Lin C ，Chen J ，Zhao X 《-》

Systemic lupus erythematosus and the risk of cardiovascular diseases: A two-sample Mendelian randomization study.

Previous observational studies have suggested that the causal role of systemic lupus erythematosus (SLE) in the risk of cardiovascular diseases (CVDs) remained inconsistent. In this study, we aimed to investigate the causal relationship between SLE and CVDs by two-sample Mendelian randomization (MR) analysis. Genetic instruments for SLE were obtained from a public genome-wide association study (GWAS) with 4,036 patients with SLE and 6,959 controls. Summary statistical data for CVDs, including coronary artery disease (CAD), myocardial infarction (MI), atrial fibrillation (AF), ischemic stroke (IS), and its subtypes, were identified from other available GWAS meta-analyses. The inverse-variance weighted (IVW) method was used as the primary method to estimate the causal effect. The simple- and weighted-median method, MR-Egger method, and MR pleiotropy residual sum and outlier (MR-PRESSO) were provided as a supplement to the IVW method. Besides, we performed sensitivity analyses, including Cochran's Q test, MR-Egger intercept test, and leave-one-out analysis, to evaluate the robustness of the results. A total of 15 single-nucleotide polymorphisms (SNPs) were identified after excluding linkage disequilibrium (LD) and potential confounding factors. According to the IVW results, our MR study indicated that genetically predicted SLE was not causally connected with the risk of CVDs [CAD: odds ratio (OR) = 1.005, 95% confidence interval (CI) = 0.986-1.024, p-value = 0.619; MI: OR = 1.002, 95% CI = 0.982-1.023, p-value = 0.854; AF: OR = 0.998, 95% CI = 0.982-1.014, p-value = 0.795; IS: OR = 1.006, 95% CI = 0.984-1.028, p-value = 0.621; cardioembolic stroke (CES): OR = 0.992, 95% CI = 0.949-1.036, p-value = 0.707; small vessel stroke (SVS): OR = 1.014, 95% CI = 0.964-1.067, p-value = 0.589; large artery stroke (LAS): OR = 1.030, 95% CI = 0.968-1.096, p-value = 0.352]. Analogical findings could be observed in supplementary MR methods. Sensitivity analyses suggested that the causal estimates were robust. Our two-sample MR analysis provided no evidence that genetically determined SLE was causally associated with the risk of CVDs.

Huang S ，Huang F ，Mei C ，Tian F ，Fan Y ，Bao J ... -  《Frontiers in Cardiovascular Medicine》