- 1Department of Neurosurgery, The First Affiliated Hospital of Harbin Medical University, Harbin, Heilongjiang, China
- 2School of Life Science, Northeast Agricultural University, Harbin, Heilongjiang, China
- 3NHC Key Laboratory of Cell Transplantation, Harbin, Heilongjiang, China
Background: Hemorrhagic stroke (HS), a leading cause of death and disability worldwide, has not been clarified in terms of the underlying biomolecular mechanisms of its development. Circulating metabolites have been closely associated with HS in recent years. Therefore, we explored the causal association between circulating metabolomes and HS using Mendelian randomization (MR) analysis and identified the molecular mechanisms of effects.
Methods: We assessed the causal relationship between circulating serum metabolites (CSMs) and HS using a bidirectional two-sample MR method supplemented with five ways: weighted median, MR Egger, simple mode, weighted mode, and MR-PRESSO. The Cochran Q-test, MR-Egger intercept test, and MR-PRESSO served for the sensitivity analyses. The Steiger test and reverse MR were used to estimate reverse causality. Metabolic pathway analyses were performed using MetaboAnalyst 5.0, and genetic effects were assessed by linkage disequilibrium score regression. Significant metabolites were further synthesized using meta-analysis, and we used multivariate MR to correct for common confounders.
Results: We finally recognized four metabolites, biliverdin (OR 0.62, 95% CI 0.40–0.96, PMVMR = 0.030), linoleate (18. 2n6) (OR 0.20, 95% CI 0.08–0.54, PMVMR = 0.001),1-eicosadienoylglycerophosphocholine* (OR 2.21, 95% CI 1.02–4.76, PMVMR = 0.044),7-alpha-hydroxy-3 -oxo-4-cholestenoate (7-Hoca) (OR 0.27, 95% CI 0.09–0.77, PMVMR = 0.015) with significant causal relation to HS.
Conclusion: We demonstrated significant causal associations between circulating serum metabolites and hemorrhagic stroke. Monitoring, diagnosis, and treatment of hemorrhagic stroke by serum metabolites might be a valuable approach.
1 Introduction
In recent studies, stroke has been considered the second leading cause of death (6.6 million persons) and disability [143 million Disability Life Years (DALYs) loss] worldwide. In the last three decades, the global incidence of stroke has increased by 70%, the prevalence of stroke has increased by 85%, the mortality rate has risen by 43%, and DALYs attributable to stroke have risen by 32%, amongst which hemorrhagic stroke (HS) is a major contributor (1, 2), there is an urgent need to explore the underlying molecular mechanisms of hemorrhagic stroke. HS is the deadliest form of stroke and includes subtypes of intracerebral hemorrhage (ICH) and subarachnoid hemorrhage (SAH) (3). Currently, neuroinflammation, oxidative stress, chronic damage to the vascular wall, cell death, and hemodynamic alterations have been discovered to be tightly related to the development of hemorrhagic stroke (4, 5). Thus, searching for the molecules that drive these microscopic mechanisms is particularly significant.
Circulating serum metabolites are a series of metabolic substrates and products closely associated with human health (6). Circulating serum metabolites (CSMs) drive various cellular functions, such as energy production and storage, signal transduction, and apoptosis (7). They can act in the development of a variety of diseases, such as chronic kidney disease, oncological diseases, and myocardial infarction (8–10). As for the nervous system, lipid metabolites with a function of triggering nerve repair were discovered (11), and a group of circulating serum metabolites consisting of serine, isoleucine, betaine, PC (5:0/5:0), and LysoPE (18:2) were judged to be excellent biomarkers for the diagnostic and progression process of acute ischemic stroke (AIS) (12). Also, a meta-analysis from seven prospective cohorts investigating the connections between serum or circulating metabolites and stroke risk found that the amino acid histidine, the glycolysis-related metabolite pyruvate, the acute-phase response marker glycoprotein acetyls, and several lipoprotein subfractions were linked to stroke risk. All stroke events, as well as ischemic and hemorrhagic events, were included in the analysis (13); likewise, after stroke onset can reduce human blood plasma metabolites such as branched-chain amino acids (valine, leucine, and isoleucine) that correlate with poor neurological outcomes (14). These studies reveal an inextricable link between CSMs and stroke.
Mendelian randomization (MR) is an application of instrumental variable analysis. It utilizes genetic variation to ascertain whether an observed association between risk factors and outcome corresponds to the causal effect. The principle refers to Mendel’s second law about the independent segregation of genetic alleles when DNA is passed from parents to offspring at gamete formation (15). Since MR analysis generates a random distribution of genetic variants through natural arbitrary classification, which is immutable, this method not only lowers the risk of confounding but also decreases the danger of reverse causality bias and is not influenced by disease state (16). Meanwhile, randomized controlled trials (RCTs) as the gold standard of causality have unavoidable problems such as substantial financial investment, weak patient compliance, and the ethics of randomized treatment allocation. Further, in situations where RCTs are unavailable, MR is extensively employed to infer a causal connection between exposure and outcome, with equally compelling results.
Based on the background mentioned above, we utilized MR analysis to detect the causal relationship between the circulating metabolome and hemorrhagic stroke to search for the metabolic pathways and mechanisms of effect behind the causality.
2 Materials and methods
2.1 Ethics statement and study design
This study did not require additional ethical approval since we used publicly available data, which had already received approval from the appropriate ethical and institutional review boards.
In this study, we evaluated the causal relationship between human blood metabolites and the risk of hemorrhagic stroke using a two-sample bidirectional MR design. A scientific MR study should comply with three hypotheses: (1) Instrumental variables (IVs) are strongly associated with the exposure of interest; (2) IVs must be independent of confounders; and (3) IVs are not related to the outcome, but only affect the outcome through exposure. And follow the Strengthening the Reporting of Observational Studies in Epidemiology Using Mendelian Randomization guidelines (STROBE-MR) (Supplementary Table S1) (17).
Notably, we each used four different GWAS summary data to perform a comprehensive MR analysis with hemorrhagic stroke as the outcome. False discovery rate correction for p-value, meta-analysis, and multivariate Mendelian randomization were performed to increase the reliability and stability of the results, linkage disequilibrium score regression (LDSC) was used by us to assess the genetic correlation between the screened exposures and the outcome, and similarly, we also explored the potential mechanisms of circulating serum metabolites for hemorrhagic stroke by metabolite pathway analyses, as outlined in the present study in Figure 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1376889/fnut-11-1376889-HTML/image_m/fnut-11-1376889-g001.jpg)
Figure 1. Overall study overview flowchart. SNP, single nucleotide polymorphism; MR, Mendelian randomization; ICH, intracerebral hemorrhage; SAH, subarachnoid hemorrhage; MR-PRESSO, MR pleiotropy residual sum and outlier.
2.2 Data source for circulating serum metabolites
We downloaded summary-type GWAS data for human serum metabolites from the Metabolomics GWAS Server.1 Notably, this is the most comprehensive report on the genetic loci of blood metabolites, which successfully screened out 486 metabolites with genetic influences on human serum metabolites by Shin et al. (18). Specifically, the study included 7,824 Europeans, including 1768 from the KORA F4 study in Germany and 6,056 from the United Kingdom Twin Study. This study analyzed 529 metabolites in plasma or serum from 7,824 adult individuals from two European population studies via liquid chromatography and gas chromatography separation-coupled tandem mass spectrometry. Four hundred and eighty-six serum metabolites were obtained after rigorous quality control. Among the 486 metabolites, 177 were defined as unknown due to poorly defined chemical properties. Another 309 metabolites were chemically authenticated and allocated to eight broad metabolic groups, including amino acid, carbohydrate, cofactors and vitamin, energy, lipid, nucleotide, peptide, and xenobiotic metabolism, as documented in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. We record additional information in Supplementary Table S2.
2.3 Data source for hemorrhagic stroke
Hemorrhagic stroke includes two subtypes, intracerebral hemorrhage (ICH) and subarachnoid hemorrhage (SAH). GWAS summary data for hemorrhagic stroke were obtained from the FinnGen and European Bioinformatics Institute (EBI) databases. For more detail, for ICH and SAH summary data, we chose the R9, R7, and R5 versions of the FinnGen database as well as the registry numbers downloaded from the GWAS Catalog for the data sets GCST90018870 and GCST90018923 datasets, the latter being the result of the meta-analysis of FinnGen R3 and UK Biobank by Saori Sakaue et al. (19). We show the specific features of the above datasets in Supplementary Table S2. More information can be acquired by accessing (https://www.finngen.fi/fi) and (https://www.ebi.ac.uk/gwas/).
2.4 Selection of instrumental variables
We selected instrumental variables to satisfy the three main assumptions mentioned above, and the first step was to screen SNPs for subsequent MR analyses by an association threshold of p < 1 × 10−5 to maximize the amount of genetic variance explained by genetic predictors (20). Second, SNPs were clustered in the European 1,000 Genomes Project Phase III reference panel using the R software with a linkage disequilibrium threshold r2 < 0.01 within 500 kilobases (kb), a condition that has been widely used in previous studies (21, 22). Finally, at the same time, the F-statistic was used as a reliable measure to assess the tool’s robustness. Ultimately, the F-statistic serves as a robust and reliable metric for evaluating the tool. We used the following formula to calculate:
Where R2 denotes explained variance, n denotes sample size, and k denotes the number of selected IVs.
2.5 Bidirectional MR analysis and sensitivity analysis
For univariable MR analysis, given that the random-effects inverse variance weighted (IVW) estimates were derived from a pooled analysis of Wald ratios across all genetic variants and the premise that IVW could provide the most accurate assessment of causal effects based on the assumption that there was no horizontal pleiotropy across all the included SNPs (23). We chose a p < 0.05 for the results of IVW as a preliminary assessment of a causal relationship between circulating serum metabolites (CSMs) and hemorrhagic stroke (HS). We also used five methods, Weighted median, MR Egger, Simple mode, Weighted mode, and MR-PRESSO, as complementary analyses to the IVW method, although the Weighted median method produces unbiased estimates even when as much as 50% of the data come from invalid instruments (24). The Weighted mode method is reliable when most individual instrument causal effect estimates come from valid instruments, even if some are considered invalid (25). Simple mode represents an unweighted empirical density function for causal estimation (26). In addition, the MR-Egger method is a valuable tool for estimating causal effects through the slope coefficients of Egger regressions, which helps to identify and address potential slight study bias (27). When the condition of IVW method p < 0.05, even though the p value of the five methods of supplementary analysis does not satisfy all <0.05, but the β value of all methods show the same direction effect, the result is recognized as positive circulating serum metabolites (28), which are included in the next stage of analysis. Meanwhile, the MR-PRESSO method results in a p < 0.05, which can further increase the accuracy and stability of positive results. As for the reverse MR analysis, we chose HS for exposure, significant and potentially causal CSM as the outcome, and applied the identical SNP selection conditions of the forward MR analysis by which to test the directionality of causality.
For the MR results of positive primary screening in the above steps, we corrected the p value of false discovery rate (FDR) according to the different types of HS and the different classifications of CSM to derive the corresponding PFDR. When the p < 0.05 and PFDR < 0.1, we considered that there was a significant causative relationship between the exposure and outcome (29). Moreover, we thought of a potential causal relationship between exposure and outcome when p < 0.05 but PFDR ≥ 0.1.
For sensitivity analyses of significant and potential causality, we employed three methods, Cochran’s Q test, MR-Egger intercept test, and MR-PRESSO, to identify horizontal pleiotropy and heterogeneity and to mitigate their effects by removing outliers. Cochran’s Q-test was used to detect heterogeneity, while the MR-Egger intercept test and MR-PRESSO (Global test) were used to detect horizontal pleiotropy of selected SNP, and the MR analysis was repeated after excluding these pleiotropic SNP. p > 0.05 indicates the absence of heterogeneity or pleiotropy (27, 30, 31). Additionally, to ensure unbiased causal estimation, we performed a leave-one-out analysis, which assesses whether the results are affected by the severity of a single SNP by discarding each SNP in turn and then performing MR analysis (32).
2.6 Linkage disequilibrium score and directionality tests
Although SNPs associated with HS were excluded during screening for IV, SNPs that can mediate the inheritance of HS still exist. And MR analyses may violate the causal effect in the presence of a genetic correlation between exposure and outcome (33). LDSC regression observes the genetic contribution of complex traits and characteristics by estimating the strength of association between SNPs and traits. Therefore, we utilized the LDSC to examine the genetic correlation between positive circulating serum metabolites and HS to ascertain that the genetic concordance of exposure and outcome did not confound the inter-causal effects. We also performed a Steiger test to eliminate bias due to reverse causality and validate causality’s directionality (34), providing a more informative interpretation of the reverse MR results.
2.7 Metabolic pathway analysis
We explored the possible pathways and mechanisms of action behind circulating serum metabolites with significant and potentially causal effects on hemorrhagic stroke using MetaboAnalyst 5.02 (35).
2.8 Meta-analysis and multivariate MR analysis
Aiming to make the final screening of circulating serum metabolites causally associated with hemorrhagic stroke comprehensive, precise, and constant, ICH and SAH in hemorrhagic stroke, we used four different GWAS data each. Then, a meta-analysis of the results of the IVW model was carried out using the R software package, followed by the screening of the corresponding PMeta, with PMeta > 0.05 removed. Meanwhile, we applied I2 to perform the heterogeneity test of meta-results (36).
To adequately adjust for confounders, circulating serum metabolites after meta-analysis selection were further included as confounders (hypertension, obesity, alcohol consumption) in our MR analysis for multivariate MR analysis (MVMR). IVW, weighted median, and MR-Egger regression were also analyzed in the MVMR analysis. Besides, Egger-intercept and Cochran’s Q tests were evaluated to assess the multiplicity and heterogeneity of the results (37).
2.9 Analysis software and packages
All the statistical analyses and data visualizations were performed using R software version 4.3.1. Bidirectional and multivariate MR analyses were implemented with the packages “TwoSampleMR” (26), “MR-PRESSO” (31), and “MendelianRandomization” (38); meta-analysis was done using the “meta” (39) package.
3 Results
3.1 Instrumental variables selection in causal analyses
Following rigorous instrumental variable selection rules and procedures, we performed MR analyses on 309 known circulating serum metabolites, acquiring IVs ranging from 3 to 485 SNPs. All selected SNPs had F values exceeding 10, showing instrumental solid efficacy.
3.2 Causal assessment of circulating metabolome leading to hemorrhagic stroke
Detailed results of the IVs and each method selected for positive MR analysis of 309 CSMs of known structure and function with HS are shown in their entirety in Supplementary Tables S3, S4, and the visualization of the results is available in Figures 2, 3. We judged that CSMs with IVW model p < 0.05 were positive for the presence of a causal effect, provided that the β values of the five analytical methods were in the same direction, and PFDR < 0.1 were considered to have a significant causal effect.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1376889/fnut-11-1376889-HTML/image_m/fnut-11-1376889-g002.jpg)
Figure 2. Causal relationships between circulating serum metabolites and intracerebral hemorrhage based on different data versions; respectively (A) FinnGen R9 (B) FinnGen R7 (C) FinnGen R5 (D) EBI database. From outside to inside, the β values of the five methods of inverse variance weighted (IVW), weighted median, MR egger, simple mode, weighted mode, the p values of IVW model, the PFDR, and the metabolite species grouping are shown sequentially.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1376889/fnut-11-1376889-HTML/image_m/fnut-11-1376889-g003.jpg)
Figure 3. Causal relationships between circulating serum metabolites and subarachnoid hemorrhage based on different data versions; respectively (A) FinnGen R9 (B) FinnGen R7 (C) FinnGen R5 (D) EBI database. From outside to inside, the β values of the five methods of inverse variance weighted (IVW), weighted median, MR egger, simple mode, weighted mode, the p values of IVW model, the PFDR, and the metabolite species grouping are shown sequentially.
As for the ICH results, following the above selection criteria, we identified a total of 39 positive CSMs in seven significant categories for the four ICH GWAS data, classified according to circulating serum metabolites, of which 18 were in Lipid, 9 in Amino acid, 3 in Peptide, 3 in Cofactors and vitamins, 3 in Xenobiotics, 1 in Carbohydrate, and 2 in Nucleotide. Corresponding β, OR, and p values can be searched in Supplementary Table S4 and Figure 2.
For the SAH results, with the same screening standard, we recognized 34 positive CSMs in six major categories for the four SAH GWAS data, according to the circulating serum metabolite classification, of which 18 belong to Lipid, 8 to Amino acid, 4 to Peptide, 2 to Carbohydrate, 1 to Cofactors and vitamins and 1 to Nucleotide. Matching β, OR, and p values are shown in Supplementary Table S2 and Figure 3.
Of the above 73 positive circulating serum metabolites after rectification (Figure 4), 7 of these (PFDR < 0.1) were determined to be significant causative CSMs, ICH was the outcome for 4, namely pyroglutamylglycine in Peptide (OR 2.13, 95% CI 1.28–3.53, p = 0.00348, PFDR = 0.0939), biliverdin in Cofactors and vitamins (OR 0.59, 95% CI 0.40–0.86, p = 0.00679, PFDR = 0.0883), linoleate (18:2n6) in Lipid (OR 0.11, 95% CI 0.04–0.33, p = 0.0000966, PFDR = 0.0118) and eicosenoate (20:1n9 or 11) (OR 0.26, 95% CI 0.13–0.52, p = 0.000141, PFDR = 0.00862) in Lipid, and 3 for SAH as an ending, for gamma-glutamylmethionine* in Peptide(OR 3.25, 95% CI 1.51–7.00, p = 0.00258, PFDR = 0.0697), 1-eicosadienoylglycerophosphocholine* in Lipid (OR 3.22, 95% CI 1.72–6.04, p = 0.000255, PFDR = 0.0311) and 7-alpha-hydroxy-3-oxo-4-cholestenoate (7-Hoca) in Lipid (OR 0.13, 95% CI 0.04–0.42, p = 0.000605, PFDR = 0.0738). The above results are presented in Table 1 and Figures 2–4, 5A.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1376889/fnut-11-1376889-HTML/image_m/fnut-11-1376889-g004.jpg)
Figure 4. Circulating serum metabolites with causal effects on outcome in different versions of hemorrhagic stroke. Demonstrate the magnitude of ORs and p-values based on the IVW model. OR, odds ratio; IVW, inverse-variance weighted.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1376889/fnut-11-1376889-HTML/image_m/fnut-11-1376889-t001.jpg)
Table 1. Results of MR estimation of circulating serum metabolites significantly and causally associated with hemorrhagic stroke.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1376889/fnut-11-1376889-HTML/image_m/fnut-11-1376889-g005.jpg)
Figure 5. (A) Forest plot of five MR analysis methods for circulating serum metabolites significantly and causally associated with hemorrhagic stroke. (B) Significant metabolic pathways enriched based on different versions of hemorrhagic stroke outcome. MetRatio = Metabolites enriched/Total number of metabolites in the metabolic pathway.
3.3 Sensitivity tests in positive causal analysis
We did multiplicity and heterogeneity tests and leave-one-out analyses in the sensitivity analyses. MR-Egger regression intercepts were not significantly different from zero (all Pegger-intercept > 0.05) except for gamma-glutamylmethionine (Pegger-intercept = 0.046) for Peptide in FinnGen R7 SAH, in which the resulting causality was not plausible since there was horizontal multiplicity in this serum substance. The Cochran’s Q test for heterogeneity p-values were all more than 0.05; after removing all outliers of X-14208-phenylalanylserine from Peptide in EBI SAH MR-PRESSO Global test P- values were more than 0.05, without significant heterogeneity and pleiotropy (Supplementary Tables S5, S6). Moreover, leave-one-out analysis demonstrated that no particular SNPs driving the correlation between CSMs and HS existed (Supplementary Table S7).
3.4 Genetic association assessment and directionality test for positive causality
Genetic correlation through linkage disequilibrium (LDSC), when assessing causality using circulating serum metabolites as exposures, indicated that 3-methylhistidine (Rg = 1.199, SERg = 0.543, PRg = 0.027) as the outcome of FinnGen R7 ICH, ursodeoxycholate (Rg = 1.026, SERg = 0.475, PRg = 0.031) as the outcome of FinnGen R5 ICH, gamma-glutamylvaline (Rg = 0.580, SERg = 0.285, PRg = 0.042) as the outcome of FinnGen R9 SAH, guanosine (Rg = 1.441, SERg = 0.725, PRg = 0.047) as the outcome of FinnGen R5 ICH, These four metabolites were genetically correlated with the outcome meaning that the results of the MR analysis mixed and coexisted with the genetic component, so the initially gained causality between them was not sufficiently robust, as well as some of the metabolites for which the genetic correlation was calculated to be negative, which we denoted as NA (Supplementary Table S8). Also, the directionality check Steiger test results implied that none of the reverse causality effects influenced the forward causality (Supplementary Table S5).
3.5 Alterations in the circulating metabolome after hemorrhagic stroke
Reverse MR analyses involving circulating serum metabolites as outcome to estimate causality included FinnGen R7 ICH as exposure to mannose (β 0.026, 95% CI 0.008–0.044, p = 0.004, PFDR = 0.065), EBI ICH as exposure to taurochenodeoxycholate (β 0.063, 95% CI 0.016–0.109, p = 0.008, PFDR = 0.067) and tetradecanedioate with EBI ICH as the exposure (β 0.047, 95% CI 0.013–0.080, p = 0.007, PFDR = 0.110) were considered to have a causal effect, whereby the levels of the mentioned three metabolites increased after the occurrence of ICH (Supplementary Tables S9, S10). Also, we did not find heterogeneity and pleiotropy, nor were SNPs significantly correlated with the results in leave-one-out analysis (Supplementary Tables S11, S12).
3.6 Effect pathway analysis of the circulating metabolome
By performing metabolic pathway analyses on GWAS data from 8 ICH and SAH as outcomes probing for causality, we characterized altogether 10 possible pathways of effect (p < 0.05), “Biosynthesis of unsaturated fatty acids” (p = 0.00088), “Linoleic acid metabolism” (p = 0.02873), “Ascorbate and aldarate metabolism” (p = 0.04562), “Porphyrin and chlorophyll metabolism” (p = 0.00212), and “Histidine metabolism” (p = 0.04069) might be the pathway mechanisms by which circulating serum metabolites contribute to the progression of ICH, while “Arginine biosynthesis” (p = 0.02687), “Glycine, serine and threonine metabolism” (p = 0.00625), “Aminoacyl-tRNA biosynthesis” (p = 0.01301), “Pantothenate and CoA biosynthesis “(p = 0.04818) may be related to the SAH disease process, “Glycerophospholipid metabolism” (p = 0.01705 for EBI ICH, p = 0.00742 for FinnGen R7 SAH, p = 0.00306 for FinnGen R9 SAH) is a probable metabolite action pathway common to the development of ICH and SAH (Supplementary Table S13 and Figure 5B).
3.7 Meta-analyses based on multiple versions of hemorrhagic stroke
Aiming to further refine the exploration and evaluation of causal effects, we conducted a meta-analysis of the IVW model for each of the seven CSMs with significant causal associations with HS based on the results of the MR analyses of the four outcome data we selected, revealing that the results of biliverdin (OR 0.72, 95% CI 0.60–0.86, Prandom effects < 0.01, I2 = 0%), pyroglutamylglycine (OR 1.83, 95% CI 1.24–2.72, Prandom effects < 0.01, I2 = 40%), linoleate (18:2n6) (OR 0.23, 95% CI 0.12–0.46, Prandom effects < 0.01, I2 = 4%) with a constant causal relationship between ICH, 1-eicosadienoylglycerophosphocholine* (OR 2.13, 95% CI 1.07–4.25, Prandom effects = 0.03, I2 = 77%) (I2 suggests that there is heterogeneity in this result, which needs to be interpreted with caution), gamma- glutamylmethionine* (OR 3.10, 95% CI 1.86–5.17, Prandom effects < 0.01, I2 = 0%), 7-alpha-hydroxy-3-oxo-4-cholestenoate (7-Hoca) (OR 0.28, 95% CI 0.14–0.55, Prandom effects < 0.01, I2 = 0%) stablely influenced the progression of SAH, while eicosenoate (20:1n9 or 11) (OR 0.59, 95% CI 0.30–1.17, Prandom effects = 0.13, I2 = 67%) with ICH causality was not plausible, thus it will be excluded from the subsequent analysis (Figures 6A,B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1376889/fnut-11-1376889-HTML/image_m/fnut-11-1376889-g006.jpg)
Figure 6. Meta-analysis of circulating serum metabolites with significant causality based on four versions of the outcome IVW model. Of these, (A) with ICH as the outcome (B) with SAH as the outcome. (C) De-confounded multivariate MR analysis incorporating alcohol consumption, obesity, and essential hypertension.
3.8 Removing confounding effects with multivariate MR analysis
To ensure that the eventually obtained CSMs directly influenced HS rather than through confounders, we incorporated alcohol consumption, hypertension, and obesity as confounders, using the respective types of HS data as the outcome and performing multivariate MR analyses considering them separately and collectively. The results of the separate analyses are presented in Supplementary Table S14, and the results of the collective studies are shown in Table 2 and Figure 6C. In summary, four CSMs were identified that directly contributed to HS, of which biliverdin belonging to Cofactors and vitamins and linoleate (18:2n6) belonging to lipid were potential protective factors for ICH, while 1-eicosadienoylglycerophosphocholine* belonging to lipid and 7- alpha-hydroxy-3-oxo-4-cholestenoate (7-Hoca) belonging to lipid, the former being a possible risk factor for SAH and the latter a potential protective factor for SAH (Figure 7). We similarly conducted the Cochran’s Q test for MVMR analysis (Multivariable IVW and Multivariable MR-Egger) and MR-Egger’s intercept without potential heterogeneity and multiplicity (Supplementary Table S15).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1376889/fnut-11-1376889-HTML/image_m/fnut-11-1376889-t002.jpg)
Table 2. Results of multivariable MR estimation of circulating serum metabolites and hemorrhagic stroke.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1376889/fnut-11-1376889-HTML/image_m/fnut-11-1376889-g007.jpg)
Figure 7. The analytical procedure and detailed information (category, OR values and outcome of action) of the 4 main metabolites were eventually obtained.
4 Discussion
For this study, we selected to explore the causal associations between 309 known CSMs and HS by using eight different large GWAS datasets on hemorrhagic stroke. We performed several bidirectional two-sample MR analyses, combined with robustness checks of the results, after correction for p-values, Pyroglutamylglycine (Peptide), biliverdin (Cofactors and vitamins), linoleate (18:2n6) (Lipid), eicosenoate (20:1n9 or 11) (Lipid), gamma- glutamylmethionine*(Peptide), 1-eicosadienoylglycerophosphocholine* (Lipid), 7-alpha-hydroxy-3-oxo-4-cholestenoate (7-Hoca) (Lipid), these 7 CSMs were recognized as having a significant causal link with HS. Subsequently, we made use of meta-analysis to thoroughly investigate the causal relationships of these CSMs with distinct versions of hemorrhagic stroke. We found that the relevance of eicosenoate (20:1n9 or 11) (Lipid) to the outcome became insignificant. Moreover, we conducted multivariate MR analysis to remove potential confounders. Ultimately, we identified four CSMs that had significant causality with HS. Specifically, biliverdin (Cofactors and vitamins) and linoleate (18:2n6) (lipid) can serve as protective factors for ICH, and increased levels of 1-eicosadienoylglycerophosphocholine* (lipid) raise the risk of SAH. In contrast, high levels of 7-alpha-hydroxy-3-oxo-4-cholestenoate (7-Hoca) (lipid) lower the risk of SAH. Besides, we also appraised the genetic correlation through LDSC from CSMs and HS to clarify the coexistence of inherited effects, strengthening the reliability of our findings. Finally, metabolic pathway analysis was applied to find ten potential action pathways binding the two, including a shared path (glycerophospholipid metabolism) that CSMs act on both ICH and SAH. In recent similar studies, Wang et al. (40) focused on the effect of selected serum metabolites produced by gut flora on stroke, and our attention is centered on the role of the more comprehensive circulating metabolome on the mechanisms of hemorrhagic stroke. Zhang et al. conducted a broad and general exploration of the causal relationship between serum metabolites and various subtypes of stroke, including hemorrhagic and ischemic stroke (41). While existing studies lack in-depth studies on each of these subtypes, our research focusing on hemorrhagic stroke fills this gap. We used multiple versions rather than a single version of the ending data to comprehensively analyze the causality, then combined the results of the individual versions utilizing meta-analysis, followed by the MVMR method to remove the interference of common confounders associated with HS, and the interlocking analyses led us to obtain a wealthy and robust causal relationship between CSMs and HS. We believe that the circulating serum metabolites discovered in this comprehensive study can offer assistance in the prevention, treatment, and deeper mechanistic investigation of hemorrhagic stroke.
Recently, a large cohort study that explored the association of L-alpha glycerylphosphorylcholine (α-GPC) with the risk of subsequent stroke at 10 years noted a higher hazard of hemorrhagic stroke in α-GPC users compared to non-users, after adjusting for traditional cerebrovascular risk factors (42). This indicates that exogenous metabolites can influence the occurrence of HS in the form of altered intake in a dose–response manner. In a meta-analysis of 26 prospective cohort studies and 12 randomized controlled trials, fish and long chain omega 3 fatty acids intake were found to be moderately negatively related to cerebrovascular risk (43). Fish can supply a wealth of diverse nutrients which, upon entering the body, make their transformation into a variety of metabolites entering the circulatory system and thereby intervening in the occurrence of cerebrovascular events, suggesting that HS is also capable of being interfered with by variations in the categories and concentrations of circulating metabolites due to food intake. As for endogenous circulating metabolites, current findings demonstrate that greater circulating concentrations of vitamin D correlate with a lower incidence of cerebrovascular disease, whereas higher levels of circulating calcium relate to an increased risk of cerebrovascular disease (44). Conversely, the post-HS CSMs also altered correspondingly. Reduced plasma concentrations of L-arginine in patients early after intracerebral hemorrhage were identified as an independent risk factor for poor outcomes (45). Meanwhile, following subarachnoid hemorrhage, nitric oxide metabolites nitrite and nitrate levels are also markedly elevated (46–48). The aforementioned discoveries confirm that CSMs and HS share a reciprocal impact. During our study, we initially identified 73 CSMs that were causally associated with HS; combined with reverse MR, we found similar mannose (carbohydrate), taurochenodeoxycholate (lipid), and tetradecanedioate (lipid), three metabolites that are capable of interfering with ICH. Mannose and taurochenodeoxycholate serve as potential protective factors for ICH, and tetradecanedioate is a potentially hazardous element in ICH. Following the occurrence of ICH, the level of all three rises. Mannose synthesizes glycoproteins and participates in immunomodulation as a widely distributed monosaccharide in body fluids and tissues (49). Hence, it is upregulated after the onset of ICH as a timely energy supplier and immunomodulator. Taurochenodeoxycholate has been linked to the inhibition of cell death (50) and, in turn, may act as a neuroprotective factor that is activated upon the rise of stress following the development of ICH. Regarding tetradecanedioate, the pathway mediating FA β-oxidation, peroxisomal FA β-oxidation or FA α-oxidation may damage vascular components contributing to ICH (51), Whereas, once ICH occurs, it might remain persistently elevated owing to its complex biological functions for example in vivo transporter protein assessment of biomarkers (52).
We further adjusted for common confounders of hemorrhagic stroke via multivariate MR analysis, which ultimately identified four metabolites having significant causal effects. To begin with, biliverdin (cofactors and vitamins) was characterized as an ICH protective element. It acts as a by-product of heme degradation, and numerous studies have confirmed that its biological functions are tightly linked to anti-inflammatory, anti-apoptotic, and anti-oxidative stress (53). Moreover, it also inhibits the expression of the pro-inflammatory cytokines interleukin 1β, tumor necrosis factor α, and interleukin 6 to achieve an anti-inflammatory role (54). Studies have shown that the anti-oxidative stress effect of biliverdin (cofactors and vitamins) can be accomplished by scavenging superoxide (55); meanwhile, biliverdin protects vascular tissue from vessel damage by reducing c-Jun NH2 terminal kinase activation and preventing endothelial cell apoptosis (56); additionally, it regulates the Nrf2/A20/eEF1A2 axis to suppress cellular death and thereby attenuates cerebral ischemia–reperfusion injury (57). However, the main mechanisms involved in the pathogenesis of ICH include damage to the cerebral vascular wall, vitellosis, and lipid deposition (4); in the latest study, oxidative stress erythrocyte-associated erythrocyte-brain endothelial interactions that can induce microglia activation in vivo also lead to cerebral hemorrhage (58). In summary, biliverdin probably blocks the occurrence and progression of ICH by decreasing vascular inflammation and cellular pyroptosis caused by the release of cellular inflammatory factors, preventing vascular endothelial apoptosis and oxidative stress in erythrocytes. Linoleate (18:2n6) (lipid) is likewise recognized as a protection factor for ICH by us, and it is an n-6 polyunsaturated fatty acid essential for average growth and development. It is found in low concentrations within the brain (<2% of total fatty acids). It constitutes a necessary precursor to benign fatty acids and arachidonic acid in the brain, which has important implications for neurodevelopment and the regulation of pain and inflammatory signaling in peripheral tissues (59). The bio function of linoleate (18:2n6) is more sophisticated. There is a belief that long-term consumption of a low-level linoleate diet might protect the brain from inflammation. It has been found that a dietary structure in which linoleate is completely deprived for an extended period leads to a rapid loss of the beneficial components of ceruloplasmin in the rat brain and, conversely, an accumulation of anti-inflammatory lipids in the rat brain following low intake of linoleate. In contrast, excessive intake promotes neuroinflammation in the rats (60, 61). The genotype is closely related to how linoleate exerts inflammatory and metabolic responses in the human body (62). Briefly, linoleate (18:2n6) may act as an ICH protective factor through an anti-inflammatory mechanism of action, but long-term low-level intake is required to fulfill this conclusion. When ingested in excess over time, it could become a risk factor for ICH by enhancing the inflammatory response. Therefore, more in-depth mechanisms need to be investigated further.
SAH usually occurs in association with aneurysm rupture, and the mechanisms behind it are primarily hemodynamic stress and vascular wall injury with inhibition of its repair (5). Apart from the significant discovery of ICH and circulating serum metabolites, we have also identified circulating serum metabolites relevant to SAH. Our analysis showed that 1-eicosadienoylglycerophosphocholine* (lipid) promoted SAH. 1-eicosadienoylglycerophosphocholine*, alias lysophosphatidylcholine (20:2), a cleavage product of phosphatidylcholine, and in addition, lysophosphatidylcholine (20:2) (LPC) plays a critical role in promoting the progression of atherosclerosis and other cardiovascular diseases by affecting endothelial cells, vascular smooth muscle cells, and arteries (63). LPC is recognized as a significant component of oxidatively damaged low-density lipoprotein (oxLDL), which induces migration of lymphocytes and macrophages, increases the production of pro-inflammatory cytokines, and can aggregate inflammation, trigger oxidative stress, and promote apoptosis. On the other hand, it activates NLRP3 and NLRC4 inflammatory vesicles in microglia and astrocytes. It also synergizes with Procaspase-1 to induce activation of the ROS promoter CYP1B1 and a robust inflammatory response in human arterial endothelial cells (64, 65). Furthermore, lysophosphatidylcholine (20:2), besides inducing oxidative stress in human endothelial cells through NOX5-mediated increases in intracellular calcium (66), also activated voltage-dependent calcium channels to increase calcium influx and enhanced 5-HT-induced contraction of vascular smooth muscle cells in umbilical arteries (67). These micro-mechanisms provide more evidence that LPC can result in cerebral neurovascular inflammatory injury, oxidative stress, and rapid hemodynamic changes, consistent with our results. Hence, elevated levels of 1-eicosadienoylglycerophosphocholine* are likely to be the mechanism behind the development of SAH. 7-Hoca is a naturally existing cholesterol metabolite in human blood, produced by the metabolism of 7α-hydroxy-4-cholestene-3-one and 27-hydroxycholesterol in the brain. It is also an intermediate metabolite in the synthesis of bile acids (BA) (68). A significant rise in the concentration of 7-Hoca after SAH was observed (69); while the conversion of 27-hydroxycholesterol to 7-Hoca is an essential mechanism for the elimination of preoxidized sterols of the brain, 27-hydroxycholesterol has been implicated in the protection of both brain and cognition after injury (70). From this, the increase in 7-Hoca after SAH may indicate accelerated metabolism of 27-hydroxycholesterol to mediate cerebral protection. BA, a downstream product of 7-Hoca, could lower plasma triglycerides by inhibiting hepatic SREBP-1c expression or modulating glucose-induced adipogenesis. Meanwhile, BA’s activation of the farnesoid X receptor (FXR) improved dyslipidemia in mice (71). So, 7-Hoca might be a metabolite acting as a cerebroprotective agent. Therefore, we consider that 7-Hoca exerts a potential protective function in SAH primarily via abolishing 27-hydroxycholesterol and facilitating BA as a downstream product, thereby mediating lipid regulation. Yet more direct mechanisms of action require more thorough exploration.
Although we have done complete and systematic research, some limitations remain. Firstly, the sources we chose for the dataset are all European, with results that are not guaranteed in terms of population generalizability. Secondly, to obtain a sufficient number of SNPs for MR analysis of the exposure data, we slightly liberalized the filtering criteria; this may have made the IVs less effective, although it is common practice among other studies, And the fact that the IVs we got are all strong instrumental variables and the directionality test has no wrong causal direction also strengthens the credibility of our results in some extent. Lastly, we have examined many metabolites of recognized functional structure, yet many unknown circulating serum metabolites could not be researched. Nevertheless, we went through a series of rigorous and comprehensive MR analyses, which led to the identification of a reliable causal relationship, which can provide a high reference value for carrying out more in-depth interaction research between CSMs and HS.
5 Conclusion
In conclusion, we comprehensively and systematically evaluated the causality among circulating serum metabolites and hemorrhagic stroke utilizing MR analysis. Four circulating serum metabolites with statistically significant and robust causal effects with hemorrhagic stroke were eventually identified. Of these, biliverdin and linoleate (18:2n6) are capable of decreasing the risk of ICH, 1-eicosadienoylglycerophosphocholine* is a hazard factor for SAH and 7-Hoca is a protective element for SAH, the mediation of inflammation, oxidative stress, apoptosis, lipid homeostasis and hemodynamics are the likely mechanisms for the action of these metabolites. Additionally, another ten notable metabolic pathways have been characterized. These results suggest that these metabolites might be considered biomarkers for HS prevention and monitoring and also provide some references and assistance for future research on the selection of circulating metabolites and the exploration of the mechanism of preventive and curative targets.
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 authors.
Author contributions
YW: Conceptualization, Investigation, Visualization, Writing – original draft, Writing – review & editing. YS: Data curation, Software, Visualization, Writing – original draft, Writing – review & editing. QL: Resources, Software, Writing – review & editing. HX: Data curation, Investigation, Resources, Writing – review & editing. AG: Data curation, Investigation, Methodology, Writing – original draft. KL: Conceptualization, Investigation, Software, Writing – review & editing. YR: Data curation, Methodology, Resources, Writing – review & editing. SG: Methodology, Resources, Writing – review & editing. HL: Funding acquisition, Supervision, Writing – review & editing. XZ: Funding acquisition, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundations of China (nos. 82271340, 82071368) and the Outstanding Youth Funding of Harbin Medical University (no. HYD2020JQ0016).
Acknowledgments
We thank the researchers and volunteers who participated in the study, FinnGen, EBI, and all the joint studies for making their data publicly available.
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2024.1376889/full#supplementary-material
Footnotes
References
1. Campbell, BCV, and Khatri, P. Stroke. Lancet (London, England). (2020) 396:129–42. doi: 10.1016/S0140-6736(20)31179-X
2. Owolabi, MO, Thrift, AG, Mahal, A, Ishida, M, Martins, S, Johnson, WD, et al. Primary stroke prevention worldwide: translating evidence into action. Lancet Public Health. (2022) 7:e74–85. doi: 10.1016/S2468-2667(21)00230-9
3. Ohashi, SN, DeLong, JH, Kozberg, MG, Mazur-Hart, DJ, van Veluw, SJ, Alkayed, NJ, et al. Role of inflammatory processes in hemorrhagic stroke. Stroke. (2023) 54:605–19. doi: 10.1161/STROKEAHA.122.037155
4. Sheth, KN. Spontaneous intracerebral hemorrhage. N Engl J Med. (2022) 387:1589–96. doi: 10.1056/NEJMra2201449
5. Claassen, J, and Park, S. Spontaneous subarachnoid haemorrhage. Lancet (London, England). (2022) 400:846–62. doi: 10.1016/S0140-6736(22)00938-2
6. Diener, C, Dai, CL, Wilmanski, T, Baloni, P, Smith, B, Rappaport, N, et al. Genome-microbiome interplay provides insight into the determinants of the human blood metabolome. Nat Metab. (2022) 4:1560–72. doi: 10.1038/s42255-022-00670-1
7. Johnson, CH, Ivanisevic, J, and Siuzdak, G. Metabolomics: beyond biomarkers and towards mechanisms. Nat Rev Mol Cell Biol. (2016) 17:451–9. doi: 10.1038/nrm.2016.25
8. Hu, JR, Coresh, J, Inker, LA, Levey, AS, Zheng, Z, Rebholz, CM, et al. Serum metabolites are associated with all-cause mortality in chronic kidney disease. Kidney Int. (2018) 94:381–9. doi: 10.1016/j.kint.2018.03.008
9. Nathan, JA. Metabolite-driven antitumor immunity. Science (New York, NY). (2022) 377:1488–9. doi: 10.1126/science.ade3697
10. Floegel, A, Kühn, T, Sookthai, D, Johnson, T, Prehn, C, Rolle-Kampczyk, U, et al. Serum metabolites and risk of myocardial infarction and ischemic stroke: a targeted metabolomic approach in two German prospective cohorts. Eur J Epidemiol. (2018) 33:55–66. doi: 10.1007/s10654-017-0333-0
11. Crunkhorn, S. Lipid metabolite triggers neural repair. Nat Rev Drug Discov. (2023) 22:787. doi: 10.1038/d41573-023-00142-5
12. Liu, P, Li, R, Antonov, AA, Wang, L, Li, W, Hua, Y, et al. Discovery of metabolite biomarkers for acute ischemic stroke progression. J Proteome Res. (2017) 16:773–9. doi: 10.1021/acs.jproteome.6b00779
13. Vojinovic, D, Kalaoja, M, Trompet, S, Fischer, K, Shipley, MJ, Li, S, et al. Association of circulating metabolites in plasma or serum and risk of stroke: meta-analysis from seven prospective cohorts. Neurology. (2020) 96:e1110–23. doi: 10.1212/WNL.0000000000011236
14. Kimberly, WT, Wang, Y, Pham, L, Furie, KL, and Gerszten, RE. Metabolite profiling identifies a branched chain amino acid signature in acute cardioembolic stroke. Stroke. (2013) 44:1389–95. doi: 10.1161/STROKEAHA.111.000397
15. Larsson, SC, Butterworth, AS, and Burgess, S. Mendelian randomization for cardiovascular diseases: principles and applications. Eur Heart J. (2023) 44:4913–24. doi: 10.1093/eurheartj/ehad736
16. Emdin, CA, Khera, AV, and Kathiresan, S. Mendelian randomization. JAMA. (2017) 318:1925–6. doi: 10.1001/jama.2017.17219
17. Skrivankova, VW, Richmond, RC, Woolf, BAR, Yarmolinsky, J, Davies, NM, Swanson, SA, et al. Strengthening the reporting of observational studies in epidemiology using Mendelian randomization: the STROBE-MR statement. JAMA. (2021) 326:1614–21. doi: 10.1001/jama.2021.18236
18. Shin, S-Y, Fauman, EB, Petersen, A-K, Krumsiek, J, Santos, R, Huang, J, et al. An atlas of genetic influences on human blood metabolites. Nat Genet. (2014) 46:543–50. doi: 10.1038/ng.2982
19. Sakaue, S, Kanai, M, Tanigawa, Y, Karjalainen, J, Kurki, M, Koshiba, S, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet. (2021) 53:1415–24. doi: 10.1038/s41588-021-00931-x
20. Sanna, S, van Zuydam, NR, Mahajan, A, Kurilshikov, A, Vich Vila, A, Võsa, U, et al. Causal relationships among the gut microbiome, short-chain fatty acids and metabolic diseases. Nat Genet. (2019) 51:600–5. doi: 10.1038/s41588-019-0350-x
21. Long, Y, Tang, L, Zhou, Y, Zhao, S, and Zhu, H. Causal relationship between gut microbiota and cancers: a two-sample Mendelian randomisation study. BMC Med. (2023) 21:66. doi: 10.1186/s12916-023-02761-6
22. Gu, Y, Jin, Q, Hu, J, Wang, X, Yu, W, Wang, Z, et al. Causality of genetically determined metabolites and metabolic pathways on osteoarthritis: a two-sample mendelian randomization study. J Transl Med. (2023) 21:357. doi: 10.1186/s12967-023-04165-9
23. Pierce, BL, and Burgess, S. Efficient design for Mendelian randomization studies: subsample and 2-sample instrumental variable estimators. Am J Epidemiol. (2013) 178:1177–84. doi: 10.1093/aje/kwt084
24. Bowden, J, Davey Smith, G, Haycock, PC, and Burgess, S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. (2016) 40:304–14. doi: 10.1002/gepi.21965
25. Hartwig, FP, Davey Smith, G, and Bowden, J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. (2017) 46:1985–98. doi: 10.1093/ije/dyx102
26. 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:7. doi: 10.7554/eLife.34408
27. Bowden, J, Davey Smith, G, and Burgess, S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. (2015) 44:512–25. doi: 10.1093/ije/dyv080
28. Chen, X, Kong, J, Diao, X, Cai, J, Zheng, J, Xie, W, et al. Depression and prostate cancer risk: a Mendelian randomization study. Cancer Med. (2020) 9:9160–7. doi: 10.1002/cam4.3493
29. Storey, JD, and Tibshirani, R. Statistical significance for genomewide studies. Proc Natl Acad Sci USA. (2003) 100:9440–5. doi: 10.1073/pnas.1530509100
30. Cohen, JF, Chalumeau, M, Cohen, R, Korevaar, DA, Khoshnood, B, and Bossuyt, PM. Cochran’s Q test was useful to assess heterogeneity in likelihood ratios in studies of diagnostic accuracy. J Clin Epidemiol. (2015) 68:299–306. doi: 10.1016/j.jclinepi.2014.09.005
31. Verbanck, M, Chen, CY, Neale, B, and 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. doi: 10.1038/s41588-018-0099-7
32. Ye, X, Liu, B, Bai, Y, Cao, Y, Lin, S, Lyu, L, et al. Genetic evidence strengthens the bidirectional connection between gut microbiota and periodontitis: insights from a two-sample Mendelian randomization study. J Transl Med. (2023) 21:674. doi: 10.1186/s12967-023-04559-9
33. O’Connor, LJ, and Price, AL. Distinguishing genetic correlation from causation across 52 diseases and complex traits. Nat Genet. (2018) 50:1728–34. doi: 10.1038/s41588-018-0255-0
34. Hindy, G, Engström, G, Larsson, SC, Traylor, M, Markus, HS, Melander, O, et al. Role of blood lipids in the development of ischemic stroke and its subtypes: a Mendelian randomization study. Stroke. (2018) 49:820–7. doi: 10.1161/STROKEAHA.117.019653
35. Xia, J, Psychogios, N, Young, N, and Wishart, DS. MetaboAnalyst: a web server for metabolomic data analysis and interpretation. Nucleic Acids Res. (2009) 37:W652–60. doi: 10.1093/nar/gkp356
36. Arya, S, Schwartz, TA, and Ghaferi, AA. Practical guide to Meta-analysis. JAMA Surg. (2020) 155:430–1. doi: 10.1001/jamasurg.2019.4523
37. D’Urso, S, Arumugam, P, Weider, T, Hwang, LD, Bond, TA, Kemp, JP, et al. Mendelian randomization analysis of factors related to ovulation and reproductive function and endometrial cancer risk. BMC Med. (2022) 20:419. doi: 10.1186/s12916-022-02585-w
38. Yavorska, OO, and Burgess, S. MendelianRandomization: an R package for performing Mendelian randomization analyses using summarized data. Int J Epidemiol. (2017) 46:1734–9. doi: 10.1093/ije/dyx034
39. Bertrand, A, Kostine, M, Barnetche, T, Truchetet, ME, and Schaeverbeke, T. Immune related adverse events associated with anti-CTLA-4 antibodies: systematic review and meta-analysis. BMC Med. (2015) 13:211. doi: 10.1186/s12916-015-0455-8
40. Wang, Q, Dai, H, Hou, T, Hou, Y, Wang, T, Lin, H, et al. Dissecting causal relationships between gut microbiota, blood metabolites, and stroke: a Mendelian randomization study. Journal of stroke. (2023) 25:350–60. doi: 10.5853/jos.2023.00381
41. Zhang, T, Cao, Y, Zhao, J, Yao, J, and Liu, G. Assessing the causal effect of genetically predicted metabolites and metabolic pathways on stroke. J Transl Med. (2023) 21:822. doi: 10.1186/s12967-023-04677-4
42. Lee, G, Choi, S, Chang, J, Choi, D, Son, JS, Kim, K, et al. Association of L-α Glycerylphosphorylcholine with subsequent stroke risk after 10 years. JAMA Netw Open. (2021) 4:e2136008. doi: 10.1001/jamanetworkopen.2021.36008
43. Chowdhury, R, Stevens, S, Gorman, D, Pan, A, Warnakula, S, Chowdhury, S, et al. Association between fish consumption, long chain omega 3 fatty acids, and risk of cerebrovascular disease: systematic review and meta-analysis. BMJ. (2012) 345:e6698. doi: 10.1136/bmj.e6698
44. Chowdhury, R, Stevens, S, Ward, H, Chowdhury, S, Sajjad, A, and Franco, OH. Circulating vitamin D, calcium and risk of cerebrovascular disease: a systematic review and meta-analysis. Eur J Epidemiol. (2012) 27:581–91. doi: 10.1007/s10654-012-9729-z
45. Mader, MM, Böger, R, Appel, D, Schwedhelm, E, Haddad, M, Mohme, M, et al. Intrathecal and systemic alterations of L-arginine metabolism in patients after intracerebral hemorrhage. J Cereb. Blood Flow Metab. (2021) 41:1964–77. doi: 10.1177/0271678X20983216
46. Suzuki, M, Asahara, H, Endo, S, Inada, K, Doi, M, Kuroda, K, et al. Increased levels of nitrite/nitrate in the cerebrospinal fluid of patients with subarachnoid hemorrhage. Neurosurg Rev. (1999) 22:96–8. doi: 10.1007/s101430050038
47. Sadamitsu, D, Kuroda, Y, Nagamitsu, T, Tsuruta, R, Inoue, T, Ueda, T, et al. Cerebrospinal fluid and plasma concentrations of nitric oxide metabolites in postoperative patients with subarachnoid hemorrhage. Crit Care Med. (2001) 29:77–9. doi: 10.1097/00003246-200101000-00018
48. Suzuki, Y, Osuka, K, Noda, A, Tanazawa, T, Takayasu, M, Shibuya, M, et al. Nitric oxide metabolites in the cisternal cerebral spinal fluid of patients with subarachnoid hemorrhage. Neurosurgery. (1997) 41:807–12. doi: 10.1097/00006123-199710000-00008
49. Wei, Z, Huang, L, Cui, L, and Zhu, X. Mannose: good player and assister in pharmacotherapy. Biomed. Pharmacother. (2020) 129:110420. doi: 10.1016/j.biopha.2020.110420
50. Zhao, Q, Dai, MY, Huang, RY, Duan, JY, Zhang, T, Bao, WM, et al. Parabacteroides distasonis ameliorates hepatic fibrosis potentially via modulating intestinal bile acid metabolism and hepatocyte pyroptosis in male mice. Nat Commun. (2023) 14:1829. doi: 10.1038/s41467-023-37459-z
51. Wanders, RJ, Komen, J, and Kemp, S. Fatty acid omega-oxidation as a rescue pathway for fatty acid oxidation disorders in humans. FEBS J. (2011) 278:182–94. doi: 10.1111/j.1742-4658.2010.07947.x
52. Müller, F, Sharma, A, König, J, and Fromm, MF. Biomarkers for in vivo assessment of transporter function. Pharmacol Rev. (2018) 70:246–77. doi: 10.1124/pr.116.013326
53. Ayer, A, Zarjou, A, Agarwal, A, and Stocker, R. Heme Oxygenases in cardiovascular health and disease. Physiol Rev. (2016) 96:1449–508. doi: 10.1152/physrev.00003.2016
54. Fondevila, C, Shen, XD, Tsuchiyashi, S, Yamashita, K, Csizmadia, E, Lassman, C, et al. Biliverdin therapy protects rat livers from ischemia and reperfusion injury. Hepatology (Baltimore, Md). (2004) 40:1333–41. doi: 10.1002/hep.20480
55. Vasavda, C, Kothari, R, Malla, AP, Tokhunts, R, Lin, A, Ji, M, et al. Bilirubin links Heme metabolism to neuroprotection by scavenging superoxide. Cell Chem. Biol. (2019) 26:1450–60.e7. doi: 10.1016/j.chembiol.2019.07.006
56. Nakao, A, Murase, N, Ho, C, Toyokawa, H, Billiar, TR, and Kanno, S. Biliverdin administration prevents the formation of intimal hyperplasia induced by vascular injury. Circulation. (2005) 112:587–91. doi: 10.1161/CIRCULATIONAHA.104.509778
57. Bai, W, Huo, S, Zhou, G, Li, J, Yang, Y, and Shao, J. Biliverdin modulates the Nrf2/A20/eEF1A2 axis to alleviate cerebral ischemia-reperfusion injury by inhibiting pyroptosis. Biomed Pharmacother. (2023) 165:115057. doi: 10.1016/j.biopha.2023.115057
58. Zhang, H, Sumbria, RK, Chang, R, Sun, J, Cribbs, DH, Holmes, TC, et al. Erythrocyte-brain endothelial interactions induce microglial responses and cerebral microhemorrhages in vivo. J Neuroinflammation. (2023) 20:265. doi: 10.1186/s12974-023-02932-5
59. Taha, AY. Linoleic acid-good or bad for the brain? NPJ Sci Food. (2020) 4:1. doi: 10.1038/s41538-019-0061-9
60. Taha, AY, Blanchard, HC, Cheon, Y, Ramadan, E, Chen, M, Chang, L, et al. Dietary linoleic acid lowering reduces lipopolysaccharide-induced increase in brain arachidonic acid metabolism. Mol Neurobiol. (2017) 54:4303–15. doi: 10.1007/s12035-016-9968-1
61. Lin, LE, Chen, CT, Hildebrand, KD, Liu, Z, Hopperton, KE, and Bazinet, RP. Chronic dietary n-6 PUFA deprivation leads to conservation of arachidonic acid and more rapid loss of DHA in rat brain phospholipids. J Lipid Res. (2015) 56:390–402. doi: 10.1194/jlr.M055590
62. Lankinen, MA, Fauland, A, Shimizu, BI, Ågren, J, Wheelock, CE, Laakso, M, et al. Inflammatory response to dietary linoleic acid depends on FADS1 genotype. Am J Clin Nutr. (2019) 109:165–75. doi: 10.1093/ajcn/nqy287
63. Liu, P, Zhu, W, Chen, C, Yan, B, Zhu, L, Chen, X, et al. The mechanisms of lysophosphatidylcholine in the development of diseases. Life Sci. (2020) 247:117443. doi: 10.1016/j.lfs.2020.117443
64. Freeman, L, Guo, H, David, CN, Brickey, WJ, Jha, S, and Ting, JP. NLR members NLRC4 and NLRP3 mediate sterile inflammasome activation in microglia and astrocytes. J Exp Med. (2017) 214:1351–70. doi: 10.1084/jem.20150237
65. Lu, Y, Nanayakkara, G, Sun, Y, Liu, L, Xu, K, Ct, D, et al. Procaspase-1 patrolled to the nucleus of proatherogenic lipid LPC-activated human aortic endothelial cells induces ROS promoter CYP1B1 and strong inflammation. Redox Biol. (2021) 47:102142. doi: 10.1016/j.redox.2021.102142
66. da Silva, JF, Alves, JV, Silva-Neto, JA, Costa, RM, Neves, KB, Alves-Lopes, R, et al. Lysophosphatidylcholine induces oxidative stress in human endothelial cells via NOX5 activation - implications in atherosclerosis. Clin Sci (Lond). (2021) 135:1845–58. doi: 10.1042/CS20210468
67. Okatani, Y, Wakatsuki, A, Watanabe, K, Taniguchi, K, and Fukaya, T. Melatonin counteracts potentiation by lysophosphatidylcholine of serotonin-induced vasoconstriction in human umbilical artery: relation to calcium influx. J Pineal Res. (2001) 30:116–22. doi: 10.1034/j.1600-079X.2001.300207.x
68. Saeed, AA, Edström, E, Pikuleva, I, Eggertsen, G, and Björkhem, I. On the importance of albumin binding for the flux of 7α-hydroxy-3-oxo-4-cholestenoic acid in the brain. J Lipid Res. (2017) 58:455–9. doi: 10.1194/jlr.P073403
69. Nagata, K, Seyama, Y, and Shimizu, T. Changes in the level of 7 alpha-hydroxy-3-oxo-4-cholestenoic acid in cerebrospinal fluid after subarachnoid hemorrhage. Neurol Med Chir. (1995) 35:294–7. doi: 10.2176/nmc.35.294
70. Meaney, S, Heverin, M, Panzenboeck, U, Ekström, L, Axelsson, M, Andersson, U, et al. Novel route for elimination of brain oxysterols across the blood-brain barrier: conversion into 7alpha-hydroxy-3-oxo-4-cholestenoic acid. J Lipid Res. (2007) 48:944–51. doi: 10.1194/jlr.M600529-JLR200
Keywords: causal association, circulating metabolome, hemorrhagic stroke, linkage disequilibrium score regression, Mendelian randomization
Citation: Wang Y, Shen Y, Li Q, Xu H, Gao A, Li K, Rong Y, Gao S, Liang H and Zhang X (2024) Exploring the causal association between genetically determined circulating metabolome and hemorrhagic stroke. Front. Nutr. 11:1376889. doi: 10.3389/fnut.2024.1376889
Edited by:
Jagannath Misra, Indiana University – Purdue University Indianapolis, United StatesReviewed by:
Georgia Damoraki, National and Kapodistrian University of Athens, GreeceDeng-Feng Zhang, Kunming Institute of Zoology, Chinese Academy of Sciences (CAS), China
Marcus Scotti, Federal University of Paraíba, Brazil
Copyright © 2024 Wang, Shen, Li, Xu, Gao, Li, Rong, Gao, Liang and Zhang. 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: Hongsheng Liang, lianghongsheng@hrbmu.edu.cn; Xiangtong Zhang, zgxgtg@sina.com
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work and share last authorship