key: cord-0701989-r36uxfh4
authors: Zheng, Cong; He, Min-hui; Huang, Jian-rong; He, Yi
title: Causal Relationships Between Social Isolation and Osteoarthritis: A Mendelian Randomization Study in European Population
date: 2021-10-14
journal: Int J Gen Med
DOI: 10.2147/ijgm.s331864
sha: 8f9643303966d1e13b548b113bdd5b716be1f340
doc_id: 701989
cord_uid: r36uxfh4

OBJECTIVE: We aimed to investigate the causal relationships between social isolation and osteoarthritis. METHODS: Publicly available genome-wide association study (GWAS) summary statistics of social isolation and osteoarthritis in European population were obtained from the Neale lab Consortium and the Medical Research Council-Integrative Epidemiology Unit (MRC-IEU) consortium, respectively. Single nucleotide polymorphisms (SNPs) associated with the traits were identified by P < 5×10(−8) and linkage disequilibrium r(2) < 0.1. Three Mendelian randomization (MR) methods including the inverse-variance weighted (IVW) method, weighted median estimator, and MR-Egger regression were used to assess the potential causal effect of social isolation on osteoarthritis and the potential causal effect of osteoarthritis on social isolation. Leave-one-out analysis and test of directional horizontal pleiotropy via MR-Egger regression were performed as sensitivity analyses. RESULTS: When evaluating the causal effect of social isolation on osteoarthritis, five SNPs (rs12364432, rs13291079, rs2352075, rs4958586, rs599550) significantly associated with social isolation were studied as instruments, and social isolation was found to be causally associated with increased risk of osteoarthritis (odds ratio [OR] 1.197 (95% confidence interval (CI) 1.096–1.308) estimated by the IVW method). When evaluating the causal effect of osteoarthritis on social isolation, six SNPs (rs10405617, rs12133235, rs13107325, rs2290573, rs3771501, rs75621460) significantly associated with osteoarthritis were studied as instruments but no causal effect of osteoarthritis on social isolation was found (OR 1.104 (95% CI 0.887–1.375) estimated by the IVW method). Consistent causal relationships were observed when estimated by the weighted median estimator and MR-Egger regression. Leave-one-out analysis and test of directional horizontal pleiotropy suggested the robustness of the above findings. CONCLUSION: Social isolation is causally associated with osteoarthritis, and further work is needed to investigate the potential mechanisms.

Social relationships are well-established determinants of both mental and physical health. 1 As one of the social determinants of health, social isolation can generally be defined as a state of complete or near-complete lack of contact between an individual and society, 2 although the way it was defined and measured might vary between studies. [3] [4] [5] An increased prevalence of social isolation has been observed, especially among the elderly, which is heightened by the coronavirus disease 2019 (COVID- 19) pandemic. 6 Therefore, to understand the potential effects of social isolation on well-being is important in a post-COVID-19 world. 7, 8 So far, social isolation has been demonstrated as a risk factor for various conditions, including coronary artery disease, stroke, and cancers. [9] [10] [11] [12] As a rheumatic musculoskeletal disorder, osteoarthritis is another highly prevalent condition in the elderly which has a considerable impact on an individual patient (eg, pain and disability) and causes heavy economic burden. 13 The association between osteoarthritis and social isolation has been investigated. A prospective observational study included about 2000 elderly individuals from six European countries and found clinical osteoarthritis increased the risk of social isolation. 14 Intuitively, it is rather reasonable to speculate that osteoarthritis may have a causal effect on developing social isolation. However, given the nature of the observational study design, association does not imply causation. Considering that bidirectional associations are common between social determinants and health, 15 reverse causality cannot be ruled out between osteoarthritis and social isolation. Obviously, it is important to investigate the potential causal relationships between osteoarthritis and social isolation, since social isolation has been widely adopted as a measure against the COVID-19 pandemic.

Methodologically, a well-conducted randomized controlled trial is an ideal tool to investigate potential causal relationships, while for social determinants, such as social isolation, it would be impracticable due to ethical concern. As an increasingly used observational epidemiological method, Mendelian randomization (MR) is an analytical method that uses genetic variants as instrumental variables to assess the causal effect between exposure and outcome. 16, 17 Similar to a natural experiment of randomization, when cells divide during meiosis, genetic variants in the populations have been randomly allocated, making MR less likely to be affected by confounding or reverse causation than conventional observational studies. 17, 18 In addition, compared to the traditional single sample MR method, as one of the statistical methods for MR, the two-sample MR method provides an opportunity to substantially increase the statistical power by incorporating data from multiple sources, in which two different study samples are used to estimate the instrument-exposure and instrument-outcome associations to estimate a causal effect of the exposure on the outcome. 18, 19 This method has been successfully applied to the investigation into the causal relationships between social isolation and Alzheimer's disease, 20 while currently there are no available studies that use the MR method to study the causal relationships between social isolation and osteoarthritis. Thus, in the present study, we used the twosample MR method to explore the potential causality between social isolation and osteoarthritis.

Via the R package "TwoSampleMR" (version 0.5.6), we accessed to the publicly available genome-wide association study (GWAS) summary statistics provided by the Integrative Epidemiology Unit (IEU) GWAS database (http://gwas-api.mrcieu.ac.uk/, Medical Research Council (MRC), University of Bristol). 21 As presented in Table 1 , genetic variants associated with social isolation came from the Neale lab Consortium, which contained 58,752 cases and 273,511 controls. The total number of single nucleotide polymorphism (SNP) was 10,894,596; genetic variants associated with osteoarthritis from the MRC-IEU consortium, which contained 38,472 osteoarthritis patients, 424,461 controls, as well as 9,851,867 SNPs.

When evaluating the causal effect of social isolation on osteoarthritis, to select the SNPs associated with social isolation as instruments, we used the threshold of statistical significance as previous (P < 5 × 10 −8 , linkage disequilibrium r 2 < 0.1). 22 Similarly, when evaluating the causal effect of osteoarthritis on social isolation, we identified the SNPs associated with osteoarthritis as instruments using the same criteria (P < 5 × 10 −8 , linkage disequilibrium r 2 < 0.1).

Three MR methods were utilized to explore the causal relationships between social isolation and osteoarthritis (including the causal effect of social isolation on osteoarthritis, and the causal effect of osteoarthritis on social isolation). The inverse-variance weighted (IVW) method is the most classical statistical MR method, which uses a meta-analysis approach to combine the Wald ratio estimates of the causal effect obtained from different SNPs. 21, 23 The IVW method relies on the assumption that all the included SNPs are valid instrumental variables, 23 weighted median estimator allows up to half of the SNPs are not valid instrumental instruments. 24 We also use the MR-Egger regression to estimate the investigated causal relationships, which do not rely on non-zero mean pleiotropy. 25 To further examine the robustness of the estimated causal relationships, a leave-one-out analysis (which estimated the causal effect after excluding every single SNP) and a test of directional horizontal pleiotropy (via the intercept term of MR-Egger regression) were performed as sensitivity analyses. The R program (R Core Team (2019). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www. R-project.org/) was used to conduct all the analyses.

When evaluating the causal effect of social isolation on osteoarthritis, five SNPs were included as instruments for social isolation (rs12364432, rs13291079, rs2352075, rs4958586, rs599550). As presented in Table 2 , the variation in the associations with social isolation was 0.056% among these five SNPs. The F statistic is a commonly used indicator for evaluating the bias of weak instrumental variables. 26 In the present study, the F statistic was 37, greater than 10 to effectively avoid the bias caused by weak instrumental variables. Associations of the five SNPs with the outcome of osteoarthritis are also presented in Table 2 .

The causal relationship between social isolation and osteoarthritis estimated by using each single SNP as the instrument is presented in Figure 1 . When the multiple SNPs were used as an instrument for social isolation, social isolation was found to be causally associated with increased risk of osteoarthritis (odds ratio [OR] 1.197 per SD increase (95% confidence interval (CI) 1.096-1.308) estimated by the IVW method). Consistent causal relationships were observed when estimated by the weighted median estimator, and MR-Egger regression (Table 3 and Figure 1 ). Figure 2 presents the relation between the effect sizes of the SNPsocial isolation association (SD units) and the SNP-osteoarthritis associations (log OR) with standard error bars, in which the slopes of the lines correspond to causal estimates using each of the three different methods.

When evaluating the causal effect of osteoarthritis on social isolation, six SNPs were included as instruments for social isolation (rs10405617, rs12133235, rs13107325, rs2290573, rs3771501, rs75621460). As presented in Table 4 , the variation in the associations with social isolation was 0.062% among these six SNPs, and the F statistic was 46, which was also greater than 10. Associations of the six SNPs with the outcome of social isolation are also presented in Table 4 .

The causal relationship between osteoarthritis and social isolation estimated by using each single SNP as the instrument is presented in Figure 3 . When the multiple SNPs were Figure 1 Forest plot of the causal effect of social isolation on osteoarthritis. Black points represent the log odds ratio for osteoarthritis per standard deviation increase in social isolation, which is produced by using each single nucleotide polymorphism (SNP) selected as a separate instrument. Red points show the combined causal estimate using all SNPs together as a single instrument, using the three different Mendelian randomization methods. Horizontal line segments denote 95% confidence intervals of the estimate. used as an instrument for osteoarthritis, no causal effect of osteoarthritis on social isolation was observed (OR 1.104 per SD increase, 95% CI 0.887-1.375, estimated by the IVW method). Consistent results were found when estimated by the weighted median estimator, and MR-Egger regression ( Table 5 , Figures 3 and 4 ).

The result of leave-one-out analysis indicated that the causal estimation was not biased or driven by any included SNP when evaluating the causal effect of social isolation on osteoarthritis ( Figure 5 ). For the test of directional horizontal pleiotropy, the P value of the pleiotropy was relatively large, and the result was less likely to be affected by genetic pleiotropy (MR-Egger regression intercept= −0.003, standard error=0.002, P=0.215), suggesting the robustness of the above findings. Result of leave-one-out analysis for the evaluation of the causal effect of osteoarthritis on social isolation is presented in Figure 6 .

To the best of our knowledge, this is the first two-sample MR analysis to investigate the causal relationships between social isolation and osteoarthritis. This study Figure 2 Scatter plot of the causal effect of social isolation on osteoarthritis. The plot presents the effect sizes of the single nucleotide polymorphism (SNP)-social isolation association (x-axis, standard deviation units) and the SNP-osteoarthritis association (y-axis, log (odds ratio)) with 95% confidence intervals. The regression slopes of the lines correspond to causal estimates using the three Mendelian randomization methods. indicated that social isolation was causally associated with an increased risk of osteoarthritis. In contrast, osteoarthritis did not causally contribute to the level of social isolation. The main implication of the current study includes:

(1) it establishes a causal relationship between social isolation and osteoarthritis; (2) given the global exposure to social isolation during the COVID-19 pandemic, the established causal relationship suggests there might be an increase in developing osteoarthritis after the pandemic, and efforts are warranted either for prevention, screening, or management; (3) the findings also warrant future studies on the underlying mechanism of the causal relationship Figure 3 Forest plot of the causal effect of osteoarthritis on social isolation. Black points represent the log odds ratio for social isolation by osteoarthritis, which is produced by using each single nucleotide polymorphism (SNP) selected as a separate instrument. Red points show the combined causal estimate using all SNPs together as a single instrument, using the three different Mendelian randomization methods. Horizontal line segments denote 95% confidence intervals of the estimate. between social isolation and osteoarthritis, which may also provide insight into novel target for management of osteoarthritis.

The main strength of our study is the MR study design, which helps to establish a causal relationship between social isolation and osteoarthritis. It has been well known that social isolation is harmful to both physical and psychological health. 1, 11, 12 Previous observational studies have already established the association between social isolation and osteoarthritis. 14, 27 However, the traditional observational study designs are susceptible to interference from reverse causality and confounding factors, which limits the understanding of the relationship between social isolation and osteoarthritis. More importantly, different from general physical risk factors, the causal relationships between mental determinants and physical conditions are usually more complicated, which could be bidirectional. The development of MR methods, especially the twosample MR method, makes it easier and feasible to explore causal relationships in an observational way. To explore whether there was also a bidirectional effect between social isolation and osteoarthritis, we used the twosample MR method to evaluate the potential causal effect of social isolation on osteoarthritis, and the potential causal effect of osteoarthritis on social isolation, respectively. Such an investigation has not been reported before, which helps to fill in the knowledge gap in this topic. In addition, Figure 4 Scatter plot of the causal effect of osteoarthritis on social isolation. The plot presents the effect sizes of the single nucleotide polymorphism (SNP)-osteoarthritis association (x-axis, standard deviation units) and the SNP-social isolation association (y-axis, log (odds ratio)) with 95% confidence intervals. The regression slopes of the lines correspond to causal estimates using the three Mendelian randomization methods. 

in this study, we included large sample sizes through the two-sample MR approach, which helps to mitigate bias due to weak instrumental variables and provides the study with enough power to detect robust causality. Several possible mechanisms may explain the causal effect of social isolation on osteoarthritis, in which the active inflammation status of individuals under social isolation may play a role. Elevated levels of proinflammatory cytokines have been observed in individuals under social isolation, including tumor necrosis factor alpha (TNFα), and interleukin-6 (IL-6). 28 Meanwhile, the level of IL-6 is associated with the degree of articular degeneration in osteoarthritis. 29 Gut microbiota could be another possible pathway. The change of gut microbiota in osteoarthritis has attracted more and more attention from researchers, 30, 31 while it has been reported that psychosocial stress under social isolation leads to the activation of the hypothalamic pituitary adrenal axis and affects gut microbiota. 32 Nevertheless, further studies are warranted since our study cannot directly provide evidence about the underlying mechanisms.

There were inevitably some limitations in our study. First, all the included participants were derived from European, so further research is needed to explore whether the finding is able to be generalizable to other populations. Second, both social isolation and osteoarthritis were based on the self-report data, rather than objective measurements, which might lead to measurement error. Third, due to the limitation of GWAS summary data, we were not able to conduct further subgroup analysis or get a more prominent view of how social isolation affected the occurrence and development of osteoarthritis. Fourth, as an inherent limitation of the two-sample MR method, sample overlap is possible in our study. Due to data limitation, we are unable to validate or improve this limitation. However, the large F statistic for the instrument-exposure association in our study suggests that this concern can be minimized. 33 Figure 5 Leave-one-out analysis of the causal effect of social isolation on osteoarthritis. Each black point represents result of the inverse-variance weighted (IVW) method applied to estimate the causal effect of social isolation on osteoarthritis excluding each particular SNP from the analysis. Each red point depicts the IVW estimate using all SNPs. No single SNP is strongly driving the overall effect of social isolation on osteoarthritis in this leave-one-out sensitivity analysis. 

This study showed that social isolation is causally associated with osteoarthritis. Further work is needed to investigate the underlying mechanisms and to improve the management of both conditions.

This article does not include any studies with human participants or animals performed by any of the authors.

This work was supported by the Science and Technology Project of Guangzhou (201707010081) and Foundation of Guangdong (2017A030313606).

Social relationships and health: a flashpoint for health policy

An overview of systematic reviews on the public health consequences of social isolation and loneliness

Rapid systematic review: the impact of social isolation and loneliness on the mental health of children and adolescents in the context of COVID-19

The potential public health relevance of social isolation and loneliness: prevalence, epidemiology, and risk factors

Social isolation and loneliness among older adults in the context of COVID-19: a global challenge

Social isolation and loneliness: imperatives for health care in a post-COVID world

Each black point represents result of the inverse-variance weighted (IVW) method applied to estimate the causal effect of osteoarthritis on social isolation excluding each particular SNP from the analysis. Each red point depicts the IVW estimate using all SNPs. No single SNP is strongly driving the overall effect of osteoarthritis on social isolation in this leave-one-out sensitivity analysis

The effects of social isolation on well-being and life satisfaction during pandemic

Social isolation and mortality in US Black and White men and women

Loneliness and depression in patients with cancer during COVID-19

Social isolation and loneliness as risk factors for myocardial infarction, stroke and mortality: UK Biobank cohort study of 479 054 men and women

Associations between social relationship measures, serum brain-derived neurotrophic factor, and risk of stroke and dementia

Osteoarthritis year in review 2019: epidemiology and therapy

Association between osteoarthritis and social isolation: data from the EPOSA study

Exploring the bidirectional associations between loneliness and cognitive functioning over 10 years: the English longitudinal study of ageing

Mendelian randomization

Evaluating the potential role of pleiotropy in Mendelian randomization studies

Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians

Recent developments in Mendelian randomization studies

Social isolation, social interaction, and Alzheimer's disease: a Mendelian randomization study

The MR-Base platform supports systematic causal inference across the human phenome. Elife

The causal relationship between body mass index and the risk of osteoarthritis

A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization

Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator

Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption

Collaboration CCG. Avoiding bias from weak instruments in Mendelian randomization studies

Mental health-ill health differences in disease severity and its sociodemographic biobehavioral predictors among patients with knee osteoarthritis

Loneliness and acute stress reactivity: a systematic review of psychophysiological studies

Muscle strength, muscle balance, physical function and plasma interleukin-6 (IL-6) levels in elderly women with knee osteoarthritis (OA)

The effects of metabolic syndrome, obesity, and the gut microbiome on load-induced osteoarthritis

Gut microbiota and obesity-associated osteoarthritis

The role of the gut-brain axis in depression: endocrine, neural, and immune pathways

Efficient design for Mendelian randomization studies: subsample and 2-sample instrumental variable estimators

The authors declare that they have no competing interests.

The International Journal of General Medicine is an international, peer-reviewed open-access journal that focuses on general and internal medicine, pathogenesis, epidemiology, diagnosis, monitoring and treatment protocols. The journal is characterized by the rapid reporting of reviews, original research and clinical studies across all disease areas. The manuscript management system is completely online and includes a very quick and fair peer-review system, which is all easy to use. Visit http://www.dovepress.com/ testimonials.php to read real quotes from published authors.