Causality of genetically determined metabolites and metabolic pathways on osteoarthritis: a two-sample mendelian randomization study
Journal of Translational Medicine volume 21, Article number: 357 (2023)
Osteoarthritis (OA) is one of the most prevalent musculoskeletal diseases and is the leading cause of pain and disability in the aged population. However, the underlying biological mechanism has not been fully understood. This study aims to reveal the causal effect of circulation metabolites on OA susceptibility.
A two-sample Mendelian Randomization (MR) analysis was performed to estimate the causality of GDMs on OA. A genome-wide association study (GWAS) of 486 metabolites was used as the exposure, whereas 8 different OA phenotypes, including any-site OA (All OA), knee and/or hip OA (knee/hip OA), knee OA, hip OA, spine OA, finger and/or thumb OA (hand OA), finger OA, thumb OA, were set the outcomes. Inverse-variance weighted (IVW) was used for calculating causal estimates. Methods including weight mode, weight median, MR-egger, and MR-PRESSO were used for the sensitive analysis. Furthermore, metabolic pathway analysis was performed via the web-based Metaconflict 4.0. All statistical analyses were performed in R software.
In this MR analysis, a total of 235 causative associations between metabolites and different OA phenotypes were observed. After false discovery rate (FDR) correction and sensitive analysis, 9 robust causative associations between 7 metabolites (e.g., arginine, kynurenine, and isovalerylcarnitine) and 5 OA phenotypes were finally identified. Additionally, eleven significant metabolic pathways in 4 OA phenotypes were identified by metabolic pathway analysis.
The finding of our study suggested that identified metabolites and metabolic pathways can be considered useful circulating metabolic biomarkers for OA screening and prevention in clinical practice, and can also serve as candidate molecules for future mechanism exploration and drug target selection.
Osteoarthritis (OA) is the most prevalent degenerative musculoskeletal disorder and is the leading cause of progressive pain and chronic disability in the aged population. OA can affect various synovial joints, with knee joints most commonly involved, followed by the lumbar spine, cervical spine, hands, ankle, and hip joints . With a growing global population and increasing life expectancy, OA has become a major public health problem due to its high prevalence and disability rate. From the recent Global Burden of Disease estimation , approximately 344 million people worldwide suffer from osteoarthritis, with a 114% increase in prevalence and a 115% increase in the population living with disability since 1990. Although the pathological process of OA formation varies in different types of joints, it is generally characterized by progressive degradation of the articular cartilage along with secondary episodic synovitis and bone remodeling . Although considerable efforts have been undertaken to understand the nature of OA, the mechanism and risk factors of OA are still elusive.
In recent years, there are growing evidence that metabolic dysregulation has a close association with the development of OA. Both in vivo and in vitro studies indicated potential links between OA and metabolic disorders including hypertension, diabetes mellitus (or insulin resistance), dyslipidemia, and obesity [4,5,6,7,8]. There are also studies suggesting that they may share common pathological processes or pathways [4, 5]. Currently, OMIC technologies, including genomics and metabolomics, have been introduced into the investigation of underlying pathophysiological mechanisms and potential therapeutic strategies for human diseases. Several recent metabolomic studies have reported numerous circulating biomarkers including amino acids, carbohydrates, and lipids in both human and animal models [9,10,11,12]. However, due to sample size limitations and confounding factors, the causal effect of blood metabolites on OA still cannot be confirmed.
With the development of high-throughput technologies, measuring hundreds of circulating metabolites and performing genotyping in large-scale populations in parallel are now capable . Genome-wide association studies (GWAS) can offer molecular insights into the complex interplay between environmental and genetic factors in the pathogenesis of diseases. Additionally, a considerable number of single nucleotide polymorphisms (SNPs) have been identified with strong associations with serum metabolites. However, there is still a great barrier to translating these genetic findings into biological mechanisms of OA development, which calls for deep analysis to reveal the causal interaction of serum metabolites on OA susceptibility.
Mendelian randomization (MR) analysis is a novel and powerful epidemiological tool that uses genetic variants as unconfounded instrumental variables to investigate the causal relationships among exposures and clinical outcomes of diseases. Exploiting the fact that genotypes are determined at conception and not generally susceptible to confounders, MR analysis can provide unbiased estimates. Recently, GWAS were extended to metabolic phenotypes which generated an atlas of genetically determined metabolites (GDMs) . Inspired by the given dataset of GDMs, we hereby employed systematic MR methods to (1) assess the causal effects of human serum metabolites on OA in both weight-bearing and non-weight-bearing joints; (2) identify common metabolites that had causal effects on multiple OA phenotype; and (3) identify metabolic pathways that might contribute to the development of the OA.
Material and methods
We systematically assessed the causal association between human circulating metabolites and the risk of OA using a two-sample MR design. The performance of a convincing MR study should comply with 3 fundamental assumptions: (1) the genetic instruments are supposed to have direct associations with exposure (i.e., metabolites in this study); (2) the genetic instruments are supposed to be unrelated with the outcome (i.e., OA in this study) and independent of any known or unknown confounding factors; (3) the effects of IVs on the outcomes are solely mediated by the exposures of interest. Genetic information for metabolites and osteoarthritis was obtained from independent GWAS datasets separately to avoid sample overlap. The overview of this MR study was presented in Fig. 1.
GWAS data for human serum metabolites
The genome-wide association summary datasets involving 486 metabolites were obtained from the study by Shin et al. . This is the currently most comprehensive analysis of human metabolites, and the full summary statistics of which were publicly available via Metabolomics GWAS Server (http://metabolomics.helmholtz-muenchen.de/gwas/). A total of 7824 adult individuals from two European cohorts (TwinsUK and KORA cohorts) and approximately 2.1 million SNPs were included in this GWAS analysis. Among the 486 metabolites, 309 were known metabolites that can be assigned to 8 broad metabolic groups (amino acids, carbohydrates, cofactors and vitamins, energy, lipids, nucleotides, peptides, and xenobiotic metabolism) as defined in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (Additional file 1: Table S2). Another 177 were unknown metabolites whose chemical identity had not been conclusively determined.
GWAS data for OA
The summary datasets of OA were derived from a GWAS meta-analysis of 826,690 individuals (177517 OA patients and 649173 controls) from 9 different populations, this study identified 100 independently associated risk variants in 11 osteoarthritis phenotypes and is the most comprehensive GWAS analysis available on osteoarthritis . After excluding surgical phenotypes that were not relevant to the purpose of this study, we selected 8 OA phenotypes for further analysis, which were any-site OA (All OA), knee and/or hip OA (knee/hip OA), knee OA, hip OA, spine OA, finger and/or thumb OA (hand OA), finger OA, thumb OA. The characteristics of the summary datasets for OA were shown in Table 1 The definition of OA in the datasets included self-reported OA, clinical diagnosis, and the tenth edition of the International Classification of Diseases (ICD-10). According to the function of these OA-involved joints, they can be divided into weight-bearing joints (knee, hip, and spine) and non-weight-bearing joints (finger, thumb). Useful information for each SNP (e.g., effect size, standard error, effect allele, and P value) was retained for further analysis.
To ensure the rationality of our analysis, we performed quality control of SNPs by removing non-dual allele SNPs, all SNPs with strand ambiguous alleles, SNPs without rsID, SNPs with duplicate rsID or base pair positions, SNPs not in the 1000 Genomes Project phase 3, SNPs with base pair positions or alleles that do not match those in 1000 Genomes Project phase 3, SNPs with interpolation information < 0.9 and all SNPs on chromosome X and Y.
Selection of instrumental variables
The selection of instrumental variables (IV) in this MR analysis was based on the 3 fundamental assumptions. Firstly, for each metabolite, we extracted SNPs with association thresholds at P < 1 × 10–5. Secondly, independent variants were identified using a clumping procedure implemented in R software, in which a linkage-disequilibrium threshold of R2 < 0.001 within a 500 kilobase (kb) distance in the European 1000 Genomes Project Phase 3 reference panel was set. Finally, to quantitatively verify whether the selected SNPs were strong instruments, we calculated the proportion of phenotypic variation explained (PVE) and the F statistic of instruments for each metabolite. Typically, a threshold of F > 10 was suggested for further MR analysis.
The causal associations between metabolites and OA for this MR analysis were mainly estimated using a standard inverse variance weighted (IVW) method. When the instrumental variables satisfy all 3 assumptions, the IVW method can provide consistent estimates of the causal effect of the exposure and is considered to be the strongest MR method. However, if some instruments contradict the IV assumptions, the analysis may give incorrect results. We have therefore performed the following sensitivity analysis: (1) Q-test for IVW and MR-Egger was used to detect potential violations of the assumption by the heterogeneity of the association between individual IVs; (2) MR-Egger was applied to estimate horizontal pleiotropy according to its intercept, ensuring that genetic variation was independently associated with exposure and outcome; (3) we applied the additional analyses of MR methods with different modeling assumptions and strengths (weighted median and weighted mode) to increase the stability and robustness of the results; (4) we applied MR-PRESSO to detect outliers and correct for horizontal pleiotropy; (5) we applied individual SNP analysis and leave-one-out analysis to assess the likelihood of associations observed by individual SNP driver.
Metabolic pathway analysis
Metabolic pathways were analyzed via the web-based Metaconflict 4.0. (https://www.metaboanalyst.ca/) . The functional enrichment analyses and pathway analyses module were used to identify potential metabolite groups or pathways that may be relevant to the biological processes of OA. Two libraries, the Small Molecule Pathway Database (SMPDB) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database were used in this study, and the significance level for pathway analysis was set at 0.10.
Statistical analyses were performed in R3.5.3 software, and MR analyses were performed using the MendelianRandomization package. MR-PRESSO was performed using the MRPRESSO package. False discovery rate (FDR) correction was used to control for false positives in multiple testing. A statistically significant association was considered if the estimated causal effect of a given metabolite had a FDR < 0.05.
Selection of IVs
The number of selected IVs for 486 metabolites ranged from 3 to 631, with a median number of 15 (Additional file 1: Table S4). These generated IVs could explain 0.0156–3.327% of the variance of their respective metabotypes. Importantly, the minimum F statistics for the validity test were all above 10 (ranging from 17.63 to 21.96) (Additional file 1: Table S4), indicating that weak instrumental bias is unlikely to occur .
Causal effects of metabolites on 8 OA phenotypes
For a better interpretation of metabolic changes, we excluded 177 unknown metabolites while including 309 with known structures and functions. Using these IVs, we estimated the causal association between these 309 metabolites and 8 OA phenotypes and identified a total of 235 suggestive associations (p < 0.05, corresponding to 140 unique metabolites) (Fig. 2). Using the IVW method, 10 causal associations with multiple-testing corrected significance (FDR < 0.05) could be observed, involving eight metabolites, which included four metabolites from the amino acid pathways, two from the lipid metabolism pathways, one from the peptide pathways and one from the xenobiotic pathways (Fig. 3). They were as follows: finger OA, 1-linoleoylglycerophosphocholine (odds ratio [OR] = 0.210, 95% confidence intervals CI 0.102–0.432, FDR < 0.001, p = 2.28 × 10–5); hand OA, X-11423–O-sulfo-L-tyrosine (OR = 0.440, 95% CI 0.293–0.663, FDR < 0.001, p = 8.41 × 10–5); hand OA, isovalerylcarnitine (OR = 1.510,95% CI 0.293–0.663,FDR < 0.001, p = 3.17 × 10–4); hand OA, 1-linoleoylglycerophosphocholine (OR = 0.357, 95% CI 0.222–0.574, FDR < 0.001, p = 2.07 × 10–5); hip OA, arginine (OR = 0.559, 95% CI 0.408–0.765, FDR < 0.001, p = 2.84 × 10–4); hip OA, ADpSGEGDFXAEGGGVR* (OR = 0.655, 95% CI 0.552–0.822, FDR < 0.001, p = 2.60 × 10–4); hip OA, 4-acetaminophen sulfate (OR = 1.023, 95% CI 1.011 = 1.036, FDR < 0.001, p = 2.50 × 10–4); knee-hip OA, kynurenine (OR = 1.438, 95% CI 1.196–1.730, FDR < 0.001, p = 1.12 × 10–4); thumb OA, taurocholate (OR = 1.352, 95% CI 1.151–1.589, FDR < 0.001, p = 2.39 × 10–4); thumb OA, X-11423–O-sulfo-L-tyrosine (OR = 0.320, 95% CI 0.186–0.553, FDR < 0.001, p = 4.32 × 10–5).
Additionally, we observed that some OA phenotypes have shared causal metabolites. For example, 1-linoleoylglycerophosphocholine had significant causal associations with both finger OA and hand OA, as well as a suggestive causal association with thumb OA (OR = 0.4398, 95% CI 0.221–0.873, p = 0.01889). X-11423-O-sulfo-L-tyrosine had causal associations with both hand OA, thumb OA, and finger OA (OR = 0.400, 95% CI 0.223–0.718, p = 0.00213). Besides knee-hip OA, kynurenine had been found to have suggestive associations with knee OA (OR = 1.435, 95% CI 1.150–1.789, p = 0.00133), hip OA (OR = 1.450, 95% CI 1.082–2.064, p = 0.00699), all OA (OR = 1.194, 95% CI 1.037–1.376, p = 0.01368) and spine OA (OR = 1.494, 95% CI 1.083–2.064, p = 0.01465). Interestingly, the associations between shared metabolites with OA phenotypes showed a trend in the differential distribution between weight-bearing and non-weight-bearing joints.
Sensitivity analyses were performed to avoid the horizontal pleiotropy for MR estimate. Figure 3 shows the results of the sensitivity analyses for 10 pairs of metabolites and OA phenotypes with significant causal associations. Generally, the causal association was robust when statistical significances (p < 0.05) were observed in two additional MR tests, typically the Weighted Median test and the MR-PRESSO test. Nine of the 10 pairs of associations were considered to be robust, which respectively were 1-linoleoylglycerophosphocholine on hand OA and finger OA, X-11423–O-sulfo-L-tyrosine on hand OA and thumb OA, isovalerylcarnitine on hand OA, arginine on hip OA, ADpSGEGDFXAEGGGVR* on hip OA, kynurenine on knee-hip OA, taurocholate on thumb OA (Fig. 3). The association of 4-acetaminophen sulfate on hip OA showed nonsignificant in the Weighted Median test (p = 0.060), but it still can be considered as a potential causal association as the significance be observed in IVW method and MR-PRESSO test (p = 1.48 × 10–4). We further screened out possible horizontal pleiotropy in all associations by MR-Egger’s intercept term and MR-PRESSO's global test (Additional file 1: Table S3). Additionally, scatter (Fig. 4) and funnel plots (Fig. 5) ruled out the possibility of potential outliers and horizontal pleiotropy for all identified metabolites. The full results of sensitive and pleiotropy analysis were listed in Additional file 1: Table S4.
Metabolic pathway analysis
The metabolic pathway analysis identified 11 significant metabolic pathways in 4 of the included OA phenotypes (Fig. 6). Our results show that the “arginine biosynthesis” (p = 5.58 × 10–6), “alanine, aspartate and glutamate metabolism” (p = 1.05 × 10–4), “D-Glutamine and D-glutamate metabolism” (p = 0.00131), “glyoxylate and dicarboxylate metabolism” (p = 0.00316) and “arginine and proline metabolism” (p = 0.00518) pathways were found to be associated with the pathogenetic process of thumb OA, whereas “valine, leucine and isoleucine biosynthesis” pathway was considered to be associated with hip OA. We also found some OA phenotypes may share common metabolic pathways, such as the “caffeine metabolism” pathway for all OA (p = 5.65 × 10–4), knee-hip OA (p = 5.50 × 10–5), and hip OA (p = 0.00243), and “aminoacyl-tRNA biosynthesis” pathway for thumb OA (p = 8.91 × 10–4) and hip OA (p = 3.46 × 10–4) (Table 2).
In this two-sample MR study, we found 235 causal relationships involving 140 metabolites and 8 OA phenotypes, nine of which were observed with multiple-testing corrected significance and considered to be robust associations as they passed all sensitivity analyses. Additionally, eleven significant metabolic pathways involved in 4 OA phenotypes were detected (Fig. 7). To our knowledge, this is the first MR study to systematically appraise the causal role of human blood metabolites in the issue of OA. Our study offers new insights to reveal the role of gene-environment interactions in the pathogenesis of osteoarthritis and to provide potential inspiration for further precision treatments.
In the recent decade, OA has been increasingly recognized as a metabolism-related disease, not only because it can be found in combination with a variety of metabolic disorders, but also because OA-related metabolites and metabolic pathways have been continuously discovered in metabolomic studies [4, 5, 19, 20]. Blood, synovial fluid, cartilage, and subchondral bone are typical sample sources for metabolomic identification . Among them, blood has been considered a good source because it contains numerous detectable metabolites and can easily be obtained in large sample sizes, thus contributing to the screening of circulation markers for OA risk [22, 23]. Metabolomic studies using plasma/serum have identified altered metabolic profiles in patients with OA, with the most commonly reported metabolites being amino acids, taurine, and phospholipids [19, 20]. Our study confirmed the presence of OA-specific metabolic profiles and further identified some key metabolites and metabolic pathways that causally contribute to the pathogenesis of OA.
Our study identified four metabolites from the amino acid metabolism pathways been causally associated with the pathogenic process of OA. Amino acids are not only the components to build peptides and proteins but also precursors for many small molecules (e.g., nitric oxide, dopamine, 5-hydroxytryptamine, polyamines, and glutathione) with irreplaceable physiological functions . Studies have established that amino acids serve as nutrients, immune modulators, and oxidative regulators in joint development and diseases . Among these amino acid-related metabolites in OA, the role of arginine and kynurenine has been increasingly focused on.
Arginine is a semi-essential amino acid involved in the urea cycle and arginine/proline metabolism . In humans, the metabolism of arginine is mainly via the arginase pathway to produce urea and ornithine; and the nitric oxide synthase (NOS) pathway to produce nitric oxide (NO) and citrulline . The anti-inflammatory and antioxidant effects of arginine have been demonstrated in various tissues and cells, and are suggested to act mainly through NO . Existing studies suggested that NO and its redox derivatives may play a protective role in the joints by inhibiting multiple inflammatory pathways and alleviating immune cell infiltrations [27, 28]. Metabolomic analysis indicated the presence of arginine depletion in patients with OA and suggested that it may be due to overactive arginine catabolism . Werdyani et al.  identified arginine deficit as one of 3 clinical endotypes of OA, with the other two being muscle weakness and low inflammatory OA. A multi-omic analysis on injury-induced OA animal models also observed an enrichment of arginine metabolism-related genes, suggesting that the arginine metabolic pathway was an important modulated metabolic pathway in OA . Choi et al.  observed that the gene encoding Arg-II, an arginine metabolizing enzyme, was specifically upregulated in cartilage samples from both human OA patients and mouse models. They also found that adenovirus-mediated overexpression of Arg-II in mouse joint tissue led to OA pathogenesis and that genetic ablation of Arg2 (Arg2-/-) in mice eliminated the pathological manifestations of OA. Li et al.  found that exogenous arginine supplementation helped to alleviate lipid peroxidation and inflammatory responses in osteoblasts-osteoarthritis cells, suggesting a novel therapeutic target for OA treatment. Our results further support these findings and highlight the importance of arginine as a protective factor in the progression of hip OA.
As a main degradation product of tryptophan, kynurenine has been shown to have direct damaging effects in a variety of tissues [32,33,34]. As an essential amino acid for protein synthesis, tryptophan undergoes extensive and complex metabolism along several pathways, resulting in many biologically active molecules acting in different organs through various action mechanisms . Metabolism of tryptophan in the joint mainly involves the kynurenine pathway and the 5-hydroxytryptophan pathway, with the kynurenine pathway being the main pathway in the pathology of OA . Lögters et al.  found that kynurenine could inhibit chondrocyte (ATDC5) cell proliferation in a dose-dependent way. In OA rat models, Wang et al.  found that activation of the kynurenine-aryl hydrocarbon receptor axis can impair chondrogenesis and chondroprotection in human umbilical cord-derived mesenchymal stromal cells. Based on the given evidence and our findings, kynurenine can be considered an important circulation biomarker and a key therapeutic target for OA.
Carnitine and acylcarnitine are involved in fatty acid β-oxidation and BCAAs (including valine, leucine, and isoleucine) metabolism, and have been demonstrated to have immunomodulatory functions in joints development and diseases . Studies showed that L-carnitine stimulated the proliferation of human primary chondrocytes and may alleviate symptoms of knee pain [40, 41]. Whereas, altered concentrations of multiple acylcarnitines, such as acetylcarnitine, hexanoylcarnitine, and butanoylcarnitine, were observed in OA patients by metabolomic studies [12, 42]. Our MR analysis identified isovalerylcarnitine, another acylcarnitine, to be causally linked to hand OA. Isovalerylcarnitine is a carnitine substrate of the isovaleryl-CoA dehydrogenase, an enzyme that is involved in the degradation of leucine and fatty acids. However, the role of isovalerylcarnitine in OA development has been rarely reported.
Another finding of this study is that we detected 1-linoleoylglycerophosphocholine from the lysolipid pathway had a protective effect on hand OA and finger OA. Dysregulated lipid metabolism has been demonstrated in previous studies to be present in OA as an important pathophysiological feature of the disease [10, 43]. Phospholipids are important components of synovial fluid and help to lubricate joints. Accumulating experimental evidence suggested that disturbed phosphatidylcholine/lysophosphatidylcholine metabolism was a characteristic metabolic alteration in OA and could be used as candidate circulation markers and drug targets [10, 44, 45]. Additionally, phosphatidylcholine/lysophosphatidylcholine metabolism was found to be a shared metabolic pathway between OA and metabolic disorders including diabetes mellitus, partly explaining the concomitant development of OA with these diseases [19, 46]. Known as a lysophosphatidylcholine, although the role of linoleoylglycerophosphocholine in OA has not been reported in the previous literature, studies have confirmed its potential as a biomarker of insulin resistance and impaired glucose metabolism . Given the fact that hand OA, especially the erosive OA subtypes, is closely associated with diabetes , 1-lysophatidylcholine could be a promising direction for further studies.
Additionally, we identified X-11423–O-sulfo-L-tyrosine from the phenylalanine and tyrosine metabolism pathway, taurocholate belonging to bile acid metabolism, ADpSGEGDFXAEGGGVR* (a fibrinogen cleavage peptide) and 4-acetaminophen sulfate (a paracetamol related drug metabolite) were causal factors of OA risk. However, mechanisms involved in these metabolites affecting OA have not been fully understood, which warrants further experimental exploration.
In this study, we also identified metabolic pathways that are causally associated with the development of OA, some of which are well documented in experimental studies for their role in the pathogenesis of OA. As mentioned above, arginine-related pathways (including arginine biosynthesis and arginine and proline metabolism) and the valine, leucine and isoleucine biosynthesis pathways have been demonstrated to be involved in the pathogenesis of OA [29, 31]. Glutamine metabolism is involved in biosynthesis and redox reactions. Glutamine metabolism is associated with biosynthesis and redox reactions and has been proven to have a role in regulating cartilage inflammation and cartilage repair [49, 50].
Importantly, our results showed that some OA phenotypes have shared metabolites and pathways. Moreover, there was a visible trend of differential distribution of these shared metabolites among OA in different types of joints. For example, the risky effects of kynurenine were more noticeable in weight-bearing joints including the spine, knee, and hip. Whereas, the causal associations of taurocholate and X-11423–O-sulfo-L-tyrosine with OA were mainly seen in non-weight-bearing joints, such as hand, finger, and thumb. Epidemiological studies have demonstrated differences in the prevalence, risk factors and susceptible populations of OA in weight-bearing and non-weight-bearing joints [1, 51]. Joint type-specific genetic variants have also been observed in GWAS analysis, indicating differences in pathogenic mechanisms between weight-bearing and non-weight-bearing joints .
The present study has several strengths. First, the most significant strength of this study is the large scale of genetic variables we covered to analyze the relationship between blood metabolites and different OA phenotypes. Specifically, excluding the ones that have not been identified, a total of 486 metabolites were covered in this study. Meanwhile, the genome-wide data source for the OA genetic variables included 826,690 individuals from 9 populations with 11 different OA phenotypes that can be broadly classified into OA of the weight-bearing joints and non-weight-bearing joints. Benefiting from these two GWAS datasets, we have made our study a relatively comprehensive and systematic analysis of the metabolic profile related to OA development. Second, using the MR design, our study largely avoided reverse causality and residual confounding factors. As extensive sensitivity analyses ruled out the possibility of variable polymorphisms, inferences on causal relationships between metabolites and OA risk in this study were considered to be robust.
There are still some limitations in our study. Firstly, this MR study was a blood metabolome-based analysis. Although blood has been suggested as a good source of sample for metabolite data, as the direct pathological site of osteoarthritis is within the joint cavity, further studies to analyze metabolite changes in synovial fluid or cartilage would help identify more promising biomarkers and drug targets. Secondly, the data for metabolites were mainly from European populations, which limited the cross-ethnic extrapolation of our findings. Thirdly, this study covered a relatively comprehensive spectrum of metabolites, however, the functions and mechanisms of some metabolites in the disease were not fully understood, which limited our interpretation of our findings from this MR analysis.
In conclusion, this is the first systematic MR analysis to evaluate the causal relationship between serum metabolites and different OA phenotypes using genome-wide data, providing preliminary evidence on the impact of circulating metabolic disturbance on the risk of OA. By IVW and multiple sensitive analysis, nine robust causal associations between 7 metabolites and 5 OA phenotypes were finally identified. Eleven significant metabolic pathways in 4 OA phenotypes were identified by metabolic pathway analysis. The finding of our study suggested that these metabolites can consider useful circulating metabolic biomarkers for OA screening and prevention in clinical practice, and can also serve as candidate molecules for future mechanism exploration and drug target selection.
Availability of data and materials
Publicly available datasets were analyzed in this study. This data can be found here: http://metabolomics.helmholtz-muenchen.de/gwas/. All the data generated by the MR analysis is in the included Supplementary Material.
Cho HJ, Morey V, Kang JY, Kim KW, Kim TK. Prevalence and risk factors of spine, shoulder, hand, hip, and knee osteoarthritis in community-dwelling koreans older than age 65 years. Clin Orthop Relat Res. 2015;473(10):3307–14.
Cieza A, Causey K, Kamenov K, Hanson SW, Chatterji S, Vos T. Global estimates of the need for rehabilitation based on the global burden of disease study 2019: a systematic analysis for the global burden of disease study 2019. Lancet. 2021;396(10267):2006–17.
Sacitharan PK. Ageing and osteoarthritis. Subcell Biochem. 2019;91:123–59.
Gandhi R, Woo KM, Zywiel MG, Rampersaud YR. Metabolic syndrome increases the prevalence of spine osteoarthritis. Orthop Surg. 2014;6(1):23–7.
Kluzek S, Newton JL, Arden NK. Is osteoarthritis a metabolic disorder? Br Med Bull. 2015;115(1):111–21.
Godziuk K, Prado CM, Woodhouse LJ, Forhan M. Prevalence of sarcopenic obesity in adults with end-stage knee osteoarthritis. Osteoarthr Cartil. 2019;27(12):1735–45.
Kendzerska T, King LK, Lipscombe L, Croxford R, Stanaitis I, Hawker GA. The impact of hip and knee osteoarthritis on the subsequent risk of incident diabetes: a population-based cohort study. Diabetologia. 2018;61(11):2290–9.
Liang X, Chou OHI, Cheung CL, Cheung BMY. Is hypertension associated with arthritis? The United States national health and nutrition examination survey 1999–2018. Ann Med. 2022;54(1):1767–75.
Southan J, McHugh E, Walker H, Ismail HM. Metabolic signature of articular cartilage following mechanical injury: an integrated transcriptomics and metabolomics analysis. Front Mol Biosci. 2020;7: 592905.
Tootsi K, Vilba K, Martson A, Kals J, Paapstel K, Zilmer M. Metabolomic signature of amino acids biogenic amines and lipids in blood serum of patients with severe osteoarthritis. Metabolites. 2020;10(8):323.
Pertusa C, Mifsut D, Morales JM, Tarin JJ, Cano A, Monleon D, Garcia-Perez MA. Metabolomic analysis of seere osteoarthritis in a Spanish population of women compared to healthy and osteoporotic subjects. Metabolites. 2022;12(8):677.
Werdyani S, Liu M, Zhang H, Sun G, Furey A, Randell EW, Rahman P, Zhai G. Endotypes of primary osteoarthritis identified by plasma metabolomics analysis. Rheumatology. 2021;60(6):2735–44.
Kettunen J, Demirkan A, Wurtz P, Draisma HH, Haller T, Rawal R, Vaarhorst A, Kangas AJ, Lyytikainen LP, Pirinen M, et al. Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA. Nat Commun. 2016;7:11122.
Shin SY, Fauman EB, Petersen AK, Krumsiek J, Santos R, Huang J, Arnold M, Erte I, Forgetta V, Yang TP, et al. An atlas of genetic influences on human blood metabolites. Nat Genet. 2014;46(6):543–50.
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40:D109-114.
Boer CG, Hatzikotoulas K, Southam L, Stefansdottir L, Zhang Y, Coutinho de Almeida R, Wu TT, Zheng J, Hartley A, Teder-Laving M, et al. Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations. Cell. 2021;184(24):6003–5.
Chong J, Xia J. 2104. Using MetaboAnalyst 4 0 for Metabolomics Data Analysis, Interpretation and Integration with Other Omics Data. Methods Mol Biol. 11 337–360.
Burgess S, Small DS, Thompson SG. A review of instrumental variable estimators for mendelian randomization. Stat Methods Med Res. 2017;26(5):2333–55.
Zhai G. Alteration of metabolic pathways in osteoarthritis. Metabolites. 2019;9(1):11.
Showiheen SAA, Sun AR, Wu X, Crawford R, Xiao Y, Wellard RM, Prasadam I. Application of metabolomics to osteoarthritis: from basic science to the clinical approach. Curr Rheumatol Rep. 2019;21(6):26.
Li JT, Zeng N, Yan ZP, Liao T, Ni GX. A review of applications of metabolomics in osteoarthritis. Clin Rheumatol. 2021;40(7):2569–79.
Maher AD, Coles C, White J, Bateman JF, Fuller ES, Burkhardt D, Little CB, Cake M, Read R, McDonagh MB, et al. 1H NMR spectroscopy of serum reveals unique metabolic fingerprints associated with subtypes of surgically induced osteoarthritis in sheep. J Proteome Res. 2012;11(8):4261–8.
Costello CA, Hu T, Liu M, Zhang W, Furey A, Fan Z, Rahman P, Randell EW, Zhai G. Metabolomics signature for non-responders to total joint replacement surgery in primary osteoarthritis patients: the newfoundland osteoarthritis study. J Orthop Res. 2020;38(4):793–802.
Hou Y, Yin Y, Wu G. Dietary essentiality of “nutritionally non-essential amino acids” for animals and humans. Exp Biol Med. 2015;240(8):997–1007.
Li Y, Xiao W, Luo W, Zeng C, Deng Z, Ren W, Wu G, Lei G. Alterations of amino acid metabolism in osteoarthritis: its implications for nutrition and health. Amino Acids. 2016;48(4):907–14.
Morris SM Jr. Arginine metabolism: boundaries of our knowledge. J Nutr. 2007;137(6 Suppl 2):1602S-1609S.
Marti ILAA, Reith W. Arginine-dependent immune responses. Cell Mol Life Sci. 2021;78(13):5303–24.
Zhai G. Clinical relevance of biochemical and metabolic changes in osteoarthritis. Adv Clin Chem. 2021;101:95–120.
Zhang W, Sun G, Likhodii S, Liu M, Aref-Eshghi E, Harper PE, Martin G, Furey A, Green R, Randell E, et al. Metabolomic analysis of human plasma reveals that arginine is depleted in knee osteoarthritis patients. Osteoarthrs Cart. 2016;24(5):827–34.
Choi WS, Yang JI, Kim W, Kim HE, Kim SK, Won Y, Son YO, Chun CH, Chun JS. Critical role for arginase II in osteoarthritis pathogenesis. Ann Rheum Dis. 2019;78(3):421–8.
Li J, Zhang Z, Huang X. l-Arginine and allopurinol supplementation attenuates inflammatory mediators in human osteoblasts-osteoarthritis cells. Int J Biol Macromol. 2018;118(Pt A):716–21.
Del’Arco AE, Argolo DS, Guillemin G, Costa MFD, Costa SL, Pinheiro AM. Neurological infection, kynurenine pathway, and parasitic infection by neospora caninum. Front Immunol. 2021;12: 714248.
Eryavuz Onmaz D, Sivrikaya A, Isik K, Abusoglu S, Albayrak Gezer I, Humeyra Yerlikaya F, Abusoglu G, Unlu A, Tezcan D. Altered kynurenine pathway metabolism in patients with ankylosing spondylitis. Int Immunopharmacol. 2021;99: 108018.
Meier TB, Savitz J. The kynurenine pathway in traumatic brain injury: implications for psychiatric outcomes. Biol Psychiatry. 2022;91(5):449–58.
Modoux M, Rolhion N, Mani S, Sokol H. Tryptophan metabolism as a pharmacological target. Trends Pharmacol Sci. 2021;42(1):60–73.
Igari T, Obara K, Ono S, Toba Y. Tryptophan metabolism in the joint diseases. Acta Vitaminol Enzymol. 1975;29(1–6):194–7.
Logters TT, Laryea MD, Jager M, Schadel-Hopfner M, Windolf J, Flohe S, Altrichter J, Scholz M, Paunel-Gorgulu AN. Kynurenine inhibits chondrocyte proliferation and is increased in synovial fluid of patients with septic arthritis. J Orthop Res. 2010;28(11):1490–6.
Wang X, Zhao Y, Li S, Wang Y, Jia C, Yang X, Li S, Zhang B, Wei W, Chang Y. Activation of the kynurenine-aryl hydrocarbon receptor axis impairs the chondrogenic and chondroprotective effects of human umbilical cord-derived mesenchymal stromal cells in osteoarthritis rats. Hum Cell. 2022;36(1):163.
Ferrara F, Bertelli A, Falchi M. Evaluation of carnitine, acetylcarnitine and isovalerylcarnitine on immune function and apoptosis. Drugs Exp Clin Res. 2005;31(3):109–14.
Stoppoloni D, Politi L, Dalla Vedova P, Messano M, Koverech A, Scandurra R, Scotto d’Abusco A. L-carnitine enhances extracellular matrix synthesis in human primary chondrocytes. Rheumatol Int. 2013;33(9):2399–403.
Malek Mahdavi A, Mahdavi R, Kolahi S. Effects of l-carnitine supplementation on serum inflammatory factors and matrix metalloproteinase enzymes in females with knee osteoarthritis: a randomized, double-blind, placebo-controlled pilot study. J Am Coll Nutr. 2016;35(7):597–603.
Zhang W, Likhodii S, Zhang Y, Aref-Eshghi E, Harper PE, Randell E, Green R, Martin G, Furey A, Sun G, et al. Classification of osteoarthritis phenotypes by metabolomics analysis. BMJ Open. 2014;4(11): e006286.
Pousinis P, Gowler PRW, Burston JJ, Ortori CA, Chapman V, Barrett DA. Lipidomic identification of plasma lipids associated with pain behaviour and pathology in a mouse model of osteoarthritis. Metabolomics. 2020;16(3):32.
Zhai G, Pelletier JP, Liu M, Randell EW, Rahman P, Martel-Pelletier J. Serum lysophosphatidylcholines to phosphatidylcholines ratio is associated with symptomatic responders to symptomatic drugs in knee osteoarthritis patients. Arthritis Res Ther. 2019;21(1):224.
Carlson AK, Rawle RA, Adams E, Greenwood MC, Bothner B, June RK. Application of global metabolomic profiling of synovial fluid for osteoarthritis biomarkers. Biochem Biophys Res Commun. 2018;499(2):182–8.
Zhai G. The role of metabolomics in precision medicine of osteoarthritis: how far are we? Osteoarthr Cartil Open. 2021;3(4): 100170.
Perez-Matos MC, Morales-Alvarez MC, Toloza FJK, Ricardo-Silgado ML, Mantilla-Rivas JO, Pinzon-Cortes JA, Perez-Mayorga M, Jimenez E, Guevara E, Mendivil CO. The Phospholipid linoleoylglycerophosphocholine as a biomarker of directly measured insulin resistance. Diabetes Metab J. 2017;41(6):466–73.
Plotz B, Bomfim F, Sohail MA, Samuels J. Current epidemiology and risk factors for the development of hand osteoarthritis. Curr Rheumatol Rep. 2021;23(8):61.
Wang C, Ying J, Niu X, Li X, Patti GJ, Shen J, O’Keefe RJ. Deletion of Glut1 in early postnatal cartilage reprograms chondrocytes toward enhanced glutamine oxidation. Bone Res. 2021;9(1):38.
Arra M, Swarnkar G, Adapala NS, Naqvi SK, Cai L, Rai MF, Singamaneni S, Mbalaviele G, Brophy R, Abu-Amer Y. Glutamine metabolism modulates chondrocyte inflammatory response. Elife. 2022;11:e80725.
Andrianakos AA, Kontelis LK, Karamitsos DG, Aslanidis SI, Georgountzos AI, Kaziolas GO, Pantelidou KV, Vafiadou EV, Dantis PC, Group, ES. Prevalence of symptomatic knee, hand, and hip osteoarthritis in Greece the ESORDIG study. J Rheumatol. 2006;33(12):2507–13.
We thank all the consortium studies for making the summary association statistics data publicly available.
This research was supported by grants from the National Natural Science Foundation of China [81902235, 82102631, 82102616] and the Shanghai Sailing Program [19YF1447600].
Ethics approval and consent to participate
This article contains human participants collected from several studies performed by previous studies. All participants gave informed consent in all the corresponding original studies, as described in the Methods. Our study is based on large-scale GWAS datasets, and not individual-level data. Hence, ethical approval was not applicable.
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1: Information of the instrumental variables an SNPs identified by Mendelian randomization estimates. Table S2: Raw data of the results from IVW analysis between metabolites and 8 osteoarthritis phenotypes. Table S3: The results of MR analysis for 10 causal associations with multiple test correction significance. Table S4: Raw data of the results of sensitivity analysis between all metabolites and osteoarthritis phenotypes. Table S5: Raw data of the results for the metabolic pathway analysis.
About this article
Cite this article
Gu, Y., Jin, Q., Hu, J. et al. Causality of genetically determined metabolites and metabolic pathways on osteoarthritis: a two-sample mendelian randomization study. J Transl Med 21, 357 (2023). https://doi.org/10.1186/s12967-023-04165-9
- Genetically determined metabolites
- Mendelian randomization