Osteoarthritis, coronary artery disease, and myocardial infarction: A mendelian randomization study

Background: Observational studies indicate that osteoarthritis (OA) and coronary artery disease (CAD), as well as myocardial infarction (MI), are often diagnosed as comorbid diseases. We performed a bidirectional Mendelian randomization (MR) study to demonstrate whether there is a causal relationship between OA, CAD, and MI.

Methods: We extracted single nucleotide polymorphisms (SNPs) related to OA in the Genetics of Osteoarthritis (GO) Consortium as instrumental variables to assess whether OA is associated with CAD and MI in the CARDIoGRAMplusC4D 1,000 Genomes genome-wide association study (GWAS). In the reverse MR, we used CAD-associated and MI-associated SNPs to the GWAS of OA to analyze their causality. These GWASs included 766,690 individuals of OA, 184,305 individuals of CAD, and 166,065 individuals of MI. MR was conducted using several methods, including the inverse variance weighted (IVW) method, the weighted median method, the MR-Egger method, and the MR-Pleiotropy RESidual Sum and Outlier (MR-PRESSO) method.

Results: The forward causal effect of OA on CAD and MI was not observed. In reverse analysis, no causal effect was discovered for CAD on the risk of OA. Notably, we observed a causal association between MI and total OA [IVW odds ratio (OR) = 0.95, 95% CI = 0.93, 0.98, P = 4E−04] and spine OA (IVW OR = 0.92, 95% CI = 0.88, 0.97, P = 0.001) but a null association between MI and knee OA, hip OA, hand OA, and thumb OA.

Conclusion: This MR study identifies a potentially protective effect of genetically predicted MI on total and spine OA risks.

Introduction

Cardiovascular diseases (CVDs), including coronary artery disease (CAD) and myocardial infarction (MI), and musculoskeletal disorders, including osteoarthritis (OA), are the main leading causes of mortality and disability worldwide (1, 2). They are often diagnosed as comorbid diseases and impair the quality of life significantly, especially among older adults (3). The possible association between CVDs and OA may be due to systemic inflammation. Since CAD, MI, and OA share common risk factors, such as aging and obesity (4–6), the association between the two prevailing conditions was explored increasingly in epidemiological studies. Jonsson et al. reported that hand OA in older women was associated with coronary atherosclerosis in the population-based multidisciplinary study of aging in the older adult population of Reykjavik (7). The presence of OA significantly increased the risk of new CAD in older adults according to the findings from the Progetto Veneto Anziano Study Cohort after 4.4 years of follow-up (8). The results from a systematic review and meta-analysis of population-based studies demonstrated that OA was related to a 31% increased risk of MI in the general population (9). However, no significantly higher risk of MI was detected in people with OA compared with the non-OA group from another systematic review and meta-analysis (10). The exact relationships between OA with CAD and MI remain unclear due to the contradictory evidence.

On the contrary, conclusions derived from observational studies were potentially biased due to residual confounding and reverse causality (11). Moreover, the cause-effect relationship between two diseases cannot be demonstrated by randomized controlled trials (RCTs) because of ethical and practical reasons (12). Mendelian randomization (MR) is a novel method that follows the law of independent assortment and uses genetic variants as instrumental variables (IVs) to assess the causal effects of exposure on the outcome (13, 14). Since the genotype of an individual is determined at conception and cannot be changed, this method largely avoids the reverse causality between the genetic phenotype and the associated outcome.

In the present study, we attempted to verify whether OA was causally associated with CAD and MI and to investigate whether CAD and/or MI were causally correlated with OA by leveraging summary genome-wide association studies (GWASs) data. A two-sample bidirectional MR framework was used to disentangle the causality and the direction of this association.

Methods

Study population

We chose genetic variants and extracted summary statistics of OA and its phenotypes from the largest genome-wide association study (GWAS) on the Genetics of Osteoarthritis (GO) Consortium participants of mainly European ancestry (177,517 cases and 649,173 controls) (15). OA was defined by the GO based on either self-reported status, the hospital diagnosed, the International Classification of Diseases-10 (ICD-10) codes, or radiography defined by the TREAT-OA Consortium (16). The study population is given in Table 1.

TABLE 1

Table 1. Baseline characteristics of the study population.

Genetic architectures of CAD were publicly available from a GWAS meta-analysis contributed to by the CARDIoGRAMplusC4D 1,000 Genomes, involving ~184,305 participants of mainly European ancestry (77%) (17). Summary statistics data of MI were obtained from the same GWAS, including 166,065 individuals, and used as a subgroup analysis. In the GWAS, CAD diagnosis included MI, acute coronary syndrome, chronic stable angina, or coronary stenosis >50%.

Ethical approval was not required in the present study as this was a secondary analysis of existing published data.

Study design

The design of the bidirectional MR study is shown in Figure 1. Briefly, the causal effects of OA on CAD and MI were first estimated. Then, the causal effects of CAD and MI on OA were analyzed. The genetic variants used as IVs must meet the following three stringent assumptions: (1) genetic variants are strongly associated with the exposures; (2) genetic variants are not associated with any modifiable confounders; and (3) the selected genetic variants are independent of any pathway that is involved in the outcome, except the exposure pathway (18). The genetic variants satisfying the second and third assumptions' mean independence from pleiotropy.

FIGURE 1

Figure 1. The design of bidirectional MR study. The solid paths are significant; dashed paths should not exist in the MR study. SNP, single nucleotide polymorphism; CAD, coronary artery disease; MI, myocardial infarction; OA, osteoarthritis; IVW, inverse variance weighted; WM, weighted median; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier.

Single nucleotide polymorphism selection

Single nucleotide polymorphisms (SNPs) identified as IVs must be strongly associated with the exposure. Thus, the selected p-value was at a genome-wide significant level (p < 5 × 10⁻⁸). Besides, we performed a linkage disequilibrium (LD) test on each SNP identified as IVs to prune the SNPs for pairs with r² > 0.01 so as to ensure independence among the SNPs.

We searched all of the SNPs associated with exposure in the PhenoScanner database to identify whether they were associated with any modifiable confounders and outcomes (19). These satisfied the last two assumptions of IVs, that is, the genetic variants involved in the outcome only through the exposure pathway. We estimated the causal effects of the SNPs on outcomes. Proxy SNPs with strong linkage disequilibrium (r² > 0.8) were used to replace SNPs not available in the outcome GWAS. If no proxy could be obtained, the SNP was excluded from our MR study. Finally, R² and F-statistics were calculated to assess the strength of IVs (20). The details of instrument SNPs are given in Supplementary Tables 1–4.

Mendelian randomization analysis

We conducted the MR analysis using several methods, including the inverse variance weighted (IVW) method, the weighted median (WM) method, the MR-Egger method, and the MR-Pleiotropy RESidual Sum and Outlier (MR-PRESSO) method. The IVW was performed as our primary MR analysis methodology to estimate the overall effect of genetically predicted exposure on outcome. The WM, MR-Egger, and MR-PRESSO methods were employed as additional sensitivity analyses to evaluate the two-sample MR assumptions (Figure 1). Specifically, the WM method provides an accurate estimate even when 50% of the genetic variants violate the core MR assumptions (21). The MR-Egger method may provide a correct estimate of causality, though no genetic variants satisfy the core MR assumptions (22). Besides, we introduced the MR-PRESSO method, which is a method for the detection and correction of outliers in IVW linear regression (23). Moreover, we used Cochran's Q statistics as a measure of heterogeneity (P-values for Cochran's Q < 0.05 suggest the existence of heterogeneity). We adopted the MR-Egger method to analyze the existence of directional pleiotropy, expressed by the size of the intercept (P-values for the MR-Egger method intercept < 0.05 suggest the existence of directional pleiotropy) (24). The “leave-one-out analysis” by removal of every single SNP at turn could assess the reliability of the results.

All of the analyses above were performed with the “TwoSampleMR” and “MR-PRESSO” packages in R software. The effect estimates of exposure on outcome were presented as odds ratios (ORs) and 95% CIs. Given the issue of multiple testing, the main results possessed statistical significance at P-values < 0.002 (0.05/24) after a Bonferroni correction.

Results

Causal effect of osteoarthritis and its phenotypes on coronary artery disease and myocardial infarction

After evaluating all the OA-associated SNPs through the PhenoScanner database and matching them with the summary data of SNP outcome (CAD or MI), we discovered that eight of them exhibited association with confounding factors (such as body mass index and/or waist-to-hip ratio), and one SNP associated with knee OA and one SNP associated with spine OA were not found in the outcome GWAS. Lastly, 16, 18, 36, 7, and 4 genetic variants with LD independence were taken as IVs for the total, knee, hip, hand, and thumb OA, respectively (Supplementary Table 1). F-statistics for IVs of OA were above the threshold of 10, suggesting that the IVs were strong instruments and reduced the bias of IV estimates.

As seen in Table 2, using the IVW method, total OA (OR = 0.93, 95% CI = 0.65, 1.32, P = 0.675) and each of its subtypes were not related to the risk of CAD. These findings were similar to the other MR estimates. Significant heterogeneity between the selected IVs of the total OA and the knee OA and CAD was observed (P < 0.05), but no heterogeneity was detected for the hip, hand, and thumb OA. The MR-Egger analysis did not suggest any directional pleiotropy for the IVs.

TABLE 2

Table 2. MR results of the association between OA and CAD.

Additionally, total OA (IVW OR = 0.95, 95% CI = 0.70, 1.29, P = 0.739) and its subtypes had no causal effect on the risk of MI, as seen in Table 3. Similar findings were obtained using the other MR estimate methods. No heterogeneity was detected by Cochran's Q-test for OA and its phenotypes, except for the total and knee OA. The MR-Egger analysis did not demonstrate any directional pleiotropy for the IVs.

TABLE 3

Table 3. MR results of the association between OA and MI.

Causal effect of coronary artery disease and myocardial infarction on osteoarthritis and its phenotypes

After calculating F-statistics and searching the PhenoScanner database, three SNPs associated with confounding factors were removed. Among the rest of the SNPs, three of them were not available in the outcome GWAS. All of the unavailable SNPs were unable to be replaced by their proxy SNPs and were removed. Last, a total of 35 LD-independent genetic variants were taken as IVs for CAD (Supplementary Table 3). Similarly, 20 LD-independent genetic variants were taken as IVs for MI.

As seen in Table 4, using the IVW method, we discovered that CAD was not associated with total OA (OR = 0.96, 95% CI = 0.94, 0.99, P = 0.008) and each of its subtypes. The results of the MR-Egger, WM, and MR-PRESSO models were consistent with IVW. Significant heterogeneity was detected by Cochran's Q-test for the total, knee, and hip OA (P < 0.05). Instead, Cochran's Q-test did not present heterogeneity for the spine, hand, and thumb OA. Moreover, the MR-Egger regression result did not reveal any bias by directional pleiotropy.

TABLE 4

Table 4. MR results of the association between CAD and OA.

As seen in Table 5, we observed that genetically predicted MI was associated with the reduced risk of total OA (IVW OR = 0.95, 95% CI = 0.93, 0.98, P = 4E−04) and spine OA (IVW OR = 0.92, 95% CI = 0.88, 0.97, P = 0.001) rather than with the knee OA, hip OA, hand OA, and thumb OA after correction for multiple comparisons. Notably, the results of the MR-Egger, WM, and MR-PRESSO models did not reflect a significant association but gave similar effect size estimates (Supplementary Figure 1). No heterogeneity was detected by Cochran's Q-test for OA and its phenotypes, except for hip OA (P < 0.05). The MR-Egger analysis did not suggest any directional pleiotropy for the IVs. The results of funnel plots and leave-one-out sensitivity analysis demonstrated that no single SNP has a significant effect on the pooled results, verifying the stability of our results, as given in Supplementary Figures 2, 3.

TABLE 5

Table 5. MR results of the association between MI and OA.

Discussion

To the best of our knowledge, this is the first two-sample MR study to investigate the bidirectional causal association between OA and CAD as well as MI. The results suggested no causality of OA on CAD and MI. In reverse analyses, genetically increased odds of MI have a causal effect on a lower risk of total and spine OA, while a null association between MI and knee OA, hip OA, hand OA, and thumb OA is provided. No causal effect was discovered for CAD on the risk of OA.

The association of total- and site-specific OA with CAD and MI is still inconclusive. Some cross-sectional studies observed higher CVD risk in patients with OA compared to controls (25, 26). A case–control study indicated that OA was significantly associated with CAD (27). In a cohort study, OA increased the risk of CVD (28). However, a prospective population-based cohort study reported that participants with OA were not at an increased risk of CVD (29). In the present MR study, we revealed a null causality of OA on CAD and MI using primary MR estimate methods (IVW method). Cochran's Q-test suggested heterogeneity in the MR (total OA to CAD and MI MR and knee OA to CAD and MI MR). The MR-Egger regression result did not present any bias by directional pleiotropy. However, all of the sensitivity analyses, including the WM, MR-Egger, and MR-PRESSO methods, supported that OA was not causally associated with CAD and MI.

Our findings of a protective relationship between MI and OA are surprising. However, the same association was also observed in observational studies. A cross-sectional study reported an inverse association between CVD and OA after adjustment for age, gender, and CVD risk factors (30). Another case–control study indicated that cardiovascular events were slightly associated with reduced hand OA risk (7). However, these observational associations cannot provide adequate evidence of causality since they are influenced by confounding factors and reverse causality. Our MR study explored the causality from the association of a selected exposure, predicted by genetic variants, with corresponding outcomes, which can avoid these biases. The MR result is in favor of the causal association between MI and total (P = 4E−04) and spine OA (P = 0.001), using the IVW method. Potential mechanisms for a protective effect of MI on OA involve inhibition of autophagy and inhibition of inflammation. After MI, cardiomyocytes appear as oxidative stress and autophagic cell death under the mediation of myocardin and p53. To inhibit the induction of autophagy, long non-coding RNA (lncRNA) cardiac autophagy inhibitory factor (CAIF) was correspondingly overexpressed (31). Meanwhile, CAIF overexpression mediates the downregulation of miR-1246 and inhibits the secretion of interleukin-6 (IL-6) from CHON-001 cells. Therefore, CAIF would inhibit inflammation through the CAIF/miR-1246/IL-6 pathway to improve OA (32). However, the inhibition of autophagy and the inhibition of inflammation may be the result of disease compensation, and, therefore, may be related to disease progression rather than disease initiation. A further longitudinal study of the MI progression rate in OA coexisting with MI may confirm this notion. We did not find causal associations between genetically predicted MI and the knee, hip, hand, or thumb OA, revealing the joint site's specificity. The site specificity of the causal relationship may be caused by the association of spondyloarthropathy with aortitis and biomechanical differences.

This study has some advantages. First, our study is the first MR study to identify the bidirectional causal relationship between multiple OA phenotypes and CAD, as well as MI, in the largest OA GWAS to date. Second, the MR method is used to infer the causal relationship between two diseases. This is not feasible in randomized controlled trials owing to an ethical problem to let patients with one disease not receive treatment to observe the development of another disease. Third, using summary data from the GWAS dataset in two-sample MR increase statistical power, especially in testing the causality of binary disease outcomes (33). Furthermore, given the possible genetic interaction between CVD and OA, we performed bidirectional MR analysis and observed negative evidence in only one direction.

However, this study has several limitations. First, with summary data from a large GWAS dataset, it is unlikely to be able to perform an analysis on the relevant stratifying risk factors, such as the duration of the disease, severity and treatment undertaken, and comorbidities. An independent cohort with individual-level data is warranted to replicate the findings later. Besides, participants with CAD and MI were not screened for OA at baseline and vice versa. Many exposed datasets with outcomes would lead to exaggerated causal estimates in the MR analysis. However, this is an inevitable limitation of two-sample MR analysis in the absence of individual-level data. Finally, not all the study participants are of European ancestry, and the variability of allele frequencies between populations may affect the results.

In conclusion, this MR study does not observe the causal effect of OA on CAD and MI or the causal effect of CAD on OA. Notably, we found a causal association between MI and total and spine OA but a null association between MI and knee OA, hip OA, hand OA, and thumb OA.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author/s.

Ethics statement

Ethical review and approval was not required for this study in accordance with the local legislation and institutional requirements.

Author contributions

HX mainly designed and performed analysis, verified data, and wrote the manuscript. YLin, HJ, and YLi performed experiments and analysis. MJ supervised the entire project. All authors have read and approved the final manuscript and provided critical feedback on intellectual content.

Acknowledgments

The authors thank all the participants and researchers for their participation in this MR study.

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/fcvm.2022.892742/full#supplementary-material

References

1. 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:2006–17. doi: 10.1016/S0140-6736(20)32340-0

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Safiri S, Kolahi AA, Cross M, Hill C, Smith E, Carson-Chahhoud K, et al. Prevalence, deaths, and disability-adjusted life years due to musculoskeletal disorders for 195 countries and territories 1990-2017. Arth Rheumatol. (2021) 73:702–14. doi: 10.1002/art.41571

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Fernandes GS, Valdes AM. Cardiovascular disease and osteoarthritis: common pathways and patient outcomes. Eur J Clin Invest. (2015) 45:405–14. doi: 10.1111/eci.12413

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Hunter DJ, Bierma-Zeinstra S. Osteoarthritis. Lancet. (2019) 393:1745–59. doi: 10.1016/S0140-6736(19)30417-9

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Messerli FH, Hofstetter L, Rimoldi SF, Rexhaj E, Bangalore S. Risk factor variability and cardiovascular outcome: jacc review topic of the week. J Am Coll Cardiol. (2019) 73:2596–603. doi: 10.1016/j.jacc.2019.02.063

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Wong ND. Epidemiological studies of chd and the evolution of preventive cardiology. Nat Rev Cardiol. (2014) 11:276–89. doi: 10.1038/nrcardio.2014.26

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Jonsson H, Helgadottir GP, Aspelund T, Eiriksdottir G, Sigurdsson S, Ingvarsson T, et al. Hand osteoarthritis in older women is associated with carotid and coronary atherosclerosis: the ages reykjavik study. Ann Rheum Dis. (2009) 68:1696–700. doi: 10.1136/ard.2008.096289

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Veronese N, Trevisan C, De Rui M, Bolzetta F, Maggi S, Zambon S, et al. Association of osteoarthritis with increased risk of cardiovascular diseases in the elderly: findings from the progetto veneto anziano study cohort. Arthr Rheumatol. (2016) 68:1136–44. doi: 10.1002/art.39564

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Schieir O, Tosevski C, Glazier RH, Hogg-Johnson S, Badley EM. Incident myocardial infarction associated with major types of arthritis in the general population: a systematic review and meta-analysis. Ann Rheum Dis. (2017) 76:1396–404. doi: 10.1136/annrheumdis-2016-210275

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Hall AJ, Stubbs B, Mamas MA, Myint PK, Smith TO. Association between osteoarthritis and cardiovascular disease: systematic review and meta-analysis. Eur J Prev Cardiol. (2016) 23:938–46. doi: 10.1177/2047487315610663

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Smith GD, Ebrahim S. 'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. (2003) 32:1–22. doi: 10.1093/ije/dyg070

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Gianicolo EAL, Eichler M, Muensterer O, Strauch K, Blettner M. Methods for evaluating causality in observational studies. Dtsch Arztebl Int. (2020) 116:101–7. doi: 10.3238/arztebl.2020.0101

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Davies NM, Holmes MV, Davey Smith G. Reading mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ. (2018) 362:k601. doi: 10.1136/bmj.k601

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Human Mol Genet. (2014) 23:R89–98. doi: 10.1093/hmg/ddu328

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Boer CG, Hatzikotoulas K, Southam L, Stefánsdóttir L, Zhang Y, Coutinho de Almeida R, et al. Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations. Cell. (2021) 184:4784–818.e17. doi: 10.1016/j.cell.2021.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Kerkhof HJ, Meulenbelt I, Akune T, Arden NK, Aromaa A, Bierma-Zeinstra SM, et al. Recommendations for standardization and phenotype definitions in genetic studies of osteoarthritis: the treat-oa consortium. Osteoarthr Cartil. (2011) 19:254–64. doi: 10.1016/j.joca.2010.10.027

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Nikpay M, Goel A, Won HH, Hall LM, Willenborg C, Kanoni S, et al. A comprehensive 1,000 genomes-based genome-wide association meta-analysis of coronary artery disease. Nat Genet. (2015) 47:1121–30. doi: 10.1038/ng.3396

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Zheng J, Baird D, Borges MC, Bowden J, Hemani G, Haycock P, et al. Recent developments in mendelian randomization studies. Curr Epidemiol Rep. (2017) 4:330–45. doi: 10.1007/s40471-017-0128-6

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Kamat MA, Blackshaw JA, Young R, Surendran P, Burgess S, Danesh J, et al. Phenoscanner V2: an expanded tool for searching human genotype-phenotype associations. Bioinformatics. (2019) 35:4851–3. doi: 10.1093/bioinformatics/btz469

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Pierce BL, Ahsan H, Vanderweele TJ. Power and instrument strength requirements for mendelian randomization studies using multiple genetic variants. Int J Epidemiol. (2011) 40:740–52. doi: 10.1093/ije/dyq151

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genetic Epidemiol. (2016) 40:304–14. doi: 10.1002/gepi.21965

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Bowden J, Davey Smith G, 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

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Verbanck M, Chen CY, Neale B, 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

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Burgess S, Thompson SG. Interpreting findings from mendelian randomization using the mr-egger method. Eur J Epidemiol. (2017) 32:377–89. doi: 10.1007/s10654-017-0255-x

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Singh G, Miller JD, Lee FH, Pettitt D, Russell MW. Prevalence of cardiovascular disease risk factors among us adults with self-reported osteoarthritis: data from the third national health and nutrition examination survey. Am J Managed Care. (2002) 8(Suppl. 15):S383–91.

PubMed Abstract | Google Scholar

26. Calvet J, Orellana C, Larrosa M, Navarro N, Chillarón JJ, Pedro-Botet J, et al. High prevalence of cardiovascular co-morbidities in patients with symptomatic knee or hand osteoarthritis. Scand J Rheumatol. (2016) 45:41–4. doi: 10.3109/03009742.2015.1054875

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Kadam UT, Jordan K, Croft PR. Clinical comorbidity in patients with osteoarthritis: a case-control study of general practice consulters in england and wales. Ann Rheum Dis. (2004) 63:408–14. doi: 10.1136/ard.2003.007526

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Kendzerska T, Jüni P, King LK, Croxford R, Stanaitis I, Hawker GA. The longitudinal relationship between hand, hip and knee osteoarthritis and cardiovascular events: a population-based cohort study. Osteoarthr Cartil. (2017) 25:1771–80. doi: 10.1016/j.joca.2017.07.024

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Hoeven TA, Leening MJ, Bindels PJ, Castaño-Betancourt M, van Meurs JB, Franco OH, et al. Disability and not osteoarthritis predicts cardiovascular disease: a prospective population-based cohort study. Ann Rheum Dis. (2015) 74:752–6. doi: 10.1136/annrheumdis-2013-204388

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Nielen MM, van Sijl AM, Peters MJ, Verheij RA, Schellevis FG, Nurmohamed MT. Cardiovascular disease prevalence in patients with inflammatory arthritis, diabetes mellitus and osteoarthritis: a cross-sectional study in primary care. BMC Muscul Disord. (2012) 13:150. doi: 10.1186/1471-2474-13-150

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Liu CY, Zhang YH, Li RB, Zhou LY, An T, Zhang RC, et al. Lncrna caif inhibits autophagy and attenuates myocardial infarction by blocking P53-mediated myocardin transcription. Nat Commun. (2018) 9:29. doi: 10.1038/s41467-017-02280-y

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Qi K, Lin R, Xue C, Liu T, Wang Y, Zhang Y, et al. Long non-coding rna (lncrna) caif is downregulated in osteoarthritis and inhibits lps-induced interleukin 6 (Il-6) upregulation by downregulation of mir-1246. Med Sci Monitor. (2019) 25:8019–24. doi: 10.12659/MSM.917135

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Lawlor DA. Commentary: two-sample mendelian randomization: opportunities and challenges. Int J Epidemiol. (2016) 45:908–15. doi: 10.1093/ije/dyw127

PubMed Abstract | CrossRef Full Text | Google Scholar