Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Causal relationship between COVID-19 and chronic pain: A mendelian randomization study

  • Yuchao Fan ,

    Roles Conceptualization, Data curation, Formal analysis, Writing – original draft, Writing – review & editing

    yuchaofan_pain@126.com

    ‡ YF and XL are joint co-first authors on this work.

    Affiliation Department of Anesthesiology, Sichuan Clinical Research Center for Cancer, Sichuan Cancer Hospital & Institute, Sichuan Cancer Center, Affiliated Cancer Hospital of University of Electronic Science and Technology of China, Chengdu, China

  • Xiao Liang

    Roles Writing – original draft

    ‡ YF and XL are joint co-first authors on this work.

    Affiliation Department of Anesthesiology, West China Hospital, Sichuan University, Chengdu, Sichuan Province, China

Abstract

Objective

COVID-19 is a highly transmissible disease that can result in long-term symptoms, including chronic pain. However, the mechanisms behind the persistence of long-COVID pain are not yet fully elucidated, highlighting the need for further research to establish causality. Mendelian randomization (MR), a statistical technique for determining a causal relationship between exposure and outcome, has been employed in this study to investigate the association between COVID-19 and chronic pain.

Material and methods

The IVW, MR Egger, and weighted median methods were employed. Heterogeneity was evaluated using Cochran’s Q statistic. MR Egger intercept and MR-PRESSO tests were performed to detect pleiotropy. The Bonferroni method was employed for the correction of multiple testing. R software was used for all statistical analyses.

Result

Based on the IVW method, hospitalized COVID-19 patients exhibit a higher risk of experiencing lower leg joint pain compared to the normal population. Meanwhile, the associations between COVID-19 hospitalization and back pain, headache, and pain all over the body were suggestive. Additionally, COVID-19 patients requiring hospitalization were found to have a suggestive higher risk of experiencing neck or shoulder pain and pain all over the body compared to those who did not require hospitalization. Patients with severe respiratory-confirmed COVID-19 showed a suggestive increased risk of experiencing pain all over the body compared to the normal population.

Conclusion

Our study highlights the link between COVID-19 severity and pain in different body regions, with implications for targeted interventions to reduce COVID-19 induced chronic pain burden.

Introduction

COVID-19, a highly infectious respiratory disease caused by the SARS-CoV-2 virus, has emerged as a global pandemic since late 2019, claiming more than 6.5 million confirmed deaths and over 600 million confirmed cases as of current records [1]. Transmission of the virus occurs predominantly via respiratory droplets emitted by infected individuals, as well as through contact with contaminated surfaces, followed by contact with the face or mouth [2]. Symptoms of COVID-19 include fever, cough, fatigue, loss of taste or smell, muscle aches, and shortness of breath, with some patients also experiencing gastrointestinal symptoms, ranging from mild to severe [3, 4]. Older adults individuals and those with underlying health conditions, such as heart disease, diabetes, and obesity, face the highest risk of severe illness and mortality [57].

Although most people recover within a few weeks after contracting COVID-19, a subset of individuals experience persistent symptoms and complications that can last for weeks or even months after the initial infection [8, 9]. This condition, known as long- COVID, is a prolonged and distressing health issue that affects a significant proportion of COVID-19 survivors [10]. Pain, in its diverse forms such as headache, muscle pain, joint pain, abdominal pain, and neuropathic pain, is a common manifestation of long- COVID [11]. The type, severity, and duration of pain can greatly vary among individuals, making management and treatment challenging [12, 13]. Although the exact mechanisms of pain in long- COVID remain incompletely understood, they may involve a combination of immune dysregulation, neurosensitivity, and inflammation [14]. Physical and emotional stress, sleep disturbances and reduced physical function may also exacerbate the onset of pain [15, 16].

The debilitating pain experienced by individuals with long-COVID has a profound impact on their quality of life, often leading to depression, anxiety, and other mental health problems [1, 9, 17]. This burden extends beyond the individual and imposes significant costs on society, including increased healthcare expenses, decreased productivity, and diminished quality of life for both patients and their families [18, 19]. Managing the pain associated with long-COVID is crucial to mitigate the physical and emotional toll it exacts. Therefore, focused monitoring and specialized management are essential to prevent and treat pain resulting from COVID-19, ultimately reducing the incidence of long-COVID.

The majority of studies on pain development following COVID-19 infection are observational, with a dearth of high-quality research to identify the specific COVID-19 patient cohorts at greater risk of long-term pain. As such, a comprehensive understanding of the factors that contribute to pain after COVID-19 is yet to be established. While various theories exist to explain the development of persistent pain after COVID-19, more rigorous investigations are necessary to establish a causal relationship between COVID-19 and long-term pain.

Mendelian randomization (MR) is a statistical method employed in epidemiology and genetics to establish causal relationships between exposures and outcomes. This technique employs genetic variants linked to the exposure of interest as IVs to estimate the causal impact of exposure on outcome [20]. MR studies offer various benefits over observational and prospective studies [2124]. Firstly, they provide robust evidence for causality by mitigating the influence of confounding factors on genetic variants. Consequently, MR studies can distinguish between potentially causal and non-causal associations. Secondly, MR studies are less prone to reverse causality bias, commonly encountered in observational studies. Thirdly, MR studies enable researchers to explore potential causal pathways for the association between exposure and outcome, providing insight into the biological mechanisms by which exposure influences outcome and identifying novel targets for intervention. Finally, MR studies can harness data from large-scale genetic consortia, offering high statistical power to detect causal effects. This study employs the MR technique to investigate the causal link between COVID-19 and pain, a potent tool to inform the development of interventions and treatments aimed at reducing the incidence of long-COVID pain.

Material and methods

Study design

This study employed a two-sample MR approach utilizing summary statistics from genome-wide association studies (GWAS) to investigate the potential causal association between different conditions with COVID-19 and bodily pain in various regions. The investigation aimed to explore the relationship between different conditions with COVID-19 as the exposure and bodily pain as the outcome, while circumventing the need for ethical approval through the use of publicly accessible data. To satisfy the IVs assumptions in MR, the study utilized genetic variation to estimate the causal effect of the exposure on the outcome. The essential criteria for genetic variation to fulfill the IVs assumptions in this study are as follows Fig 1): (1) The genetic variant must exhibit an association with the exposure; (2) the genetic variant must not exhibit an association with any confounding factors that may influence the association between the exposure and outcome; and (3) the genetic variant should not directly influence the outcome, except through its association with the exposure.

thumbnail
Fig 1. MR study design examining the causal effect of different conditions with COVID-19 on bodily pain in various regions, requiring a genetic variant that meets three criteria: (1) The genetic variant must exhibit an association with the exposure; (2) the genetic variant must not exhibit an association with any confounding factors that may influence the association between the exposure and outcome; and (3) the genetic variant should not directly influence the outcome, except through its association with the exposure.

https://doi.org/10.1371/journal.pone.0295982.g001

Data source

This study relied on genetic associations derived from independent GWAS datasets featuring the same ancestral population, in order to account for confounding factors. The COVID-19 GWAS dataset comprised four separate datasets sourced from European populations, which reflected varying levels of COVID-19 prevalence. These datasets included COVID-19 (infection vs. normal population), COVID-19 (hospitalized vs. normal population), COVID-19 (hospitalized vs. not hospitalized), and COVID-19 (very severe respiratory confirmed vs. normal population). The 13 pain-related datasets used in this study encompassed Pain in joint, Pain in joint (Lower leg), Low back pain, Low back pain (Lumbar region), Pain in limb (Lower leg), Back pain, Facial pain, Headache, Hip pain, Knee pain, Neck or shoulder pain, Stomach or abdominal pain, and Pain all over the body. These datasets reflected the prevalence of pain in the study population. Although pain is experienced throughout the body and some of the datasets depict overlapping pain sites, they complement each other by either identifying distinct pain sites or localizing to more specific areas. For example, while back pain and low back pain may overlap in certain areas, there are also areas of back pain that are not represented in low back pain. On the other hand, Low back pain (Lumbar region) corresponds to a more specific area of low back pain. A detailed overview of the datasets featured in this study is provided in Table 1.

Selection of instrumental variables

In this investigation, the selection of Instrumental variables (IVs) was conducted meticulously, adhering to strict criteria. Single-nucleotide polymorphisms (SNPs) were deemed valid IVs if they exhibited significant genome-wide association with the exposure (P < 1e-5) [25]. Moreover, SNPs with a minor allele frequency (MAF) greater than 0.01 in the outcome and a linkage disequilibrium (LD) r2 of less than 0.001 within a 10,000 kb distance were selected as IVs. SNPs linked with confounding factors or outcomes according to the Phenoscanner database (http://www.phenoscanner.medschl.cam.ac.uk/) were excluded from the study. To assess the strength of the IVs, we calculated the proportion of variance accounted for by each SNP, and the F-statistic. A value of less than 10 for the F-statistic indicated the selected genetic variants were weak IVs, which could result in biased outcomes [26]. Thus, we cautiously interpreted the findings. The determination of the attributable fraction of variance ascribed to individual SNPs was computed through the formula: R2 = 2 × β2 × EAF × (1 − EAF)/(2 × β2 × EAF × (1 − EAF) + 2 × SE2 × N × EAF × (1 − EAF)). Simultaneously, the F-statistic was ascertained using the following equation: F = (N − k − 1)/k ×R2 /(1 − R2), wherein ’N’ denotes the count of samples subjected to the GWAS, ’k’ signifies the number of IVs, and ’R2’ characterizes the degree to which IVs expound upon the exposure under investigation.

MR analysis

To explore the potential causal link between different conditions with COVID-19 and pain in distinct anatomical regions, we employed three distinct MR methodologies. The primary approach was the Inverse Variance Weighting (IVW) method, while we employed MR-Egger and weighted median methods as complementary measures. Significant substantive outcomes necessitated concurrence among the results derived from MR-Egger, the weighted median, and IVW methodologies in terms of directional implications. In the absence of such alignment, the observed significance would remain predominantly nominal in nature.

To mitigate the challenge associated with conducting numerous statistical tests, we implemented a stringent correction method based on the Bonferroni principle, setting the significance threshold at p < 0.0125 (0.05 divided by 4). P-values falling within the interval of 0.0125 to 0.05 were deemed indicative of preliminary indications of potential causal relationships, warranting subsequent validation and confirmation. The effect estimates were presented as odds ratios (ORs) with corresponding 95% confidence intervals (CIs).

We employed Cochran’s Q statistic, using MR-Egger and IVW methods, to evaluate heterogeneity [27, 28]. A p-value greater than 0.05 indicated the absence of heterogeneity. To assess the influence of individual SNPs on the causal relationship between exposure and outcome, we conducted a "leave-one-out" analysis. In situations where heterogeneity existed, we utilized a random effects model to estimate the causal association. Moreover, we performed MR-Egger intercept and MR-PRESSO tests to detect possible pleiotropy of IVs [29]. A p-value greater than 0.05 was regarded as evidence of the absence of pleiotropy. The statistical analyses for all MR procedures were carried out using the R software (version 4.2.2) and R Package "TwoSampleMR" and "MRPRESSO."

Results

The aim of our study was to examine the causal relationship between various COVID-19 conditions and bodily pain in different regions, using a two-sample MR analysis approach. This method allowed us to investigate whether COVID-19 conditions were a potential causal factor for experiencing pain in different parts of the body. In adherence to our stringent criteria for SNP selection, a defined set of SNPs has been meticulously curated to function as IVs in the context of our MR study. Detailed information pertaining to these chosen SNPs, which includes Number of SNPs extracted from exposure, Number of SNPs extracted from outcome, Number of SNPs after harmonising with outcome, SNPs removed due to harmonization, and Filtering for SNPs significantly associated with exposure in outcome data, has been meticulously documented in S1 Table for comprehensive reference. No SNP were removed due to MAF.

The primary method employed to assess causal effect was the IVW method, which revealed that although COVID-19 patients do not exhibit a correlation with pain risk overall, certain cases of COVID-19 are associated with an elevated risk of pain (S2 Table, Fig 2B). Specifically, we found that patients hospitalized with COVID-19 were at a higher risk of developing various types of pain compared to the normal population. The odds ratios (ORs) for Pain in joint (Lower leg), Back pain, Headache, Knee pain, Neck or shoulder pain, and Pain all over the body were 1.001104 (95% CI: 1.000484–1.001725; IVW P value = 0.000487), 1.003429 (95% CI: 1.000014–1.006855; IVW P value = 0.049055), 1.003588 (95% CI: 1.000320–1.006867; IVW P value = 0.031366), 1.004784 (95% CI: 1.001635–1.007942; IVW P value = 0.002877), 1.002835 (95% CI: 1.000078–1.005599; IVW P value = 0.043864), and 1.001085 (95% CI: 1.000178–1.001994; IVW P value = 0.019048), respectively (Figs 3A–3D, 4A–4D, Table 2). These results indicate that COVID-19 patients who require hospitalization are at greater risk of experiencing pain in different parts of the body than those in the general population. However, owing to the incongruity in the directional outcomes between IVW and MR-Egger analyses, the increased risk of knee pain and neck or shoulder pain among patients hospitalized with COVID-19 compared to the normal population was merely nominal (Table 2). Furthermore, considering the Bonferroni-corrected threshold set at 0.0125 for multiple testing, it is noteworthy that patients hospitalized with COVID-19 exhibit a significantly increased risk of experiencing Pain in joint (Lower leg) compared to the normal population. However, for causal effects related to back pain, headache, and pain all over the body, the evidence remains suggestive, necessitating further investigation for confirmation (Table 2, Fig 2B).

thumbnail
Fig 2. Forest plot of the potential causal association between different conditions with COVID-19 and bodily pain in various regions.

(A) COVID-19 (infection vs. normal population), (B) COVID-19 (hospitalized vs. normal population), (C) COVID-19 (hospitalized vs. not hospitalized), and (D) COVID-19 (very severe respiratory confirmed vs. normal population). The different conditions with COVID-19 has been identified as a potential risk factor in the pain occurrence, as indicated by its highlighted representation in Red and Bolded font. The symbol ’a’ signifies that in the context of MR-inverse-variance weighted methods, the direction of results deviates from the outcomes observed in other methodologies, indicating a mere nominal association. The symbol ’b’ signifies a inverse-variance weighted p-value less than 0.05 but greater than 0.0125, indicating that there is a suggestive association between the condition with COVID-19 and pain.

https://doi.org/10.1371/journal.pone.0295982.g002

thumbnail
Fig 3. Scatter plot displaying individual MR estimates of the effect of SNPs on COVID-19 and pain risk, including regression lines for different MR methods.

COVID-19 (hospitalized vs. normal population) on (A) Pain in joint(Lower leg), (B) Back pain, (C) Headache, (D) Pain all over the body; COVID-19 (hospitalized vs. not hospitalized) on (E) Neck or shoulder pain, (F) Pain all over the body; (G) COVID-19 (very severe respiratory confirmed vs. normal population) on Pain all over the body. IV, instrumental variant; IVW, inverse variance weighted; MR, Mendelian randomization; SNP, single-nucleotide polymorphism.

https://doi.org/10.1371/journal.pone.0295982.g003

thumbnail
Fig 4. Forest plot illustrating the MR effect size for different conditions with COVID-19 on pain risk for each of SNPs.

COVID-19 (hospitalized vs. normal population) on (A) Pain in joint(Lower leg), (B) Back pain, (C) Headache, (D) Pain all over the body; COVID-19 (hospitalized vs. not hospitalized) on (E) Neck or shoulder pain, (F) Pain all over the body; (G) COVID-19 (very severe respiratory confirmed vs. normal population) on Pain all over the body. MR, Mendelian randomization; SNP, single-nucleotide polymorphism.

https://doi.org/10.1371/journal.pone.0295982.g004

thumbnail
Table 2. MR estimates of associations between COVID-19 (hospitalized vs. normal population) and bodily pain in various regions across different methods.

https://doi.org/10.1371/journal.pone.0295982.t002

Additionally, COVID-19 patients who required hospitalization were found to be at a higher risk of experiencing Neck or shoulder pain (OR: 1.002855; 95% CI: 1.000450–1.005266; IVW P value = 0.019966) and Pain all over the body (OR: 1.000819; 95% CI: 1.000037–1.001601; IVW P value = 0.040095) compared to those who did not require hospitalization (Figs 2C, 3E, 3F, 4E and 4F, S3 Table). Furthermore, patients with severe respiratory confirmed COVID-19 had an increased risk for pain all over the body as compared to the general population, with an OR of 1.000638 (95% CI: 1.000016–1.001260; IVW P value = 0.044438) (Figs 2D, 3G and 4G, S4 Table). Consistent estimates were also obtained from the MR Egger and weighted median analyses (S3 and S4 Tables). Nevertheless, due to the P value exceeding 0.0125, all these associations remain in the realm of being suggestive, demanding further empirical validation through additional research endeavors.

To assess the robustness of our results, we conducted a series of sensitivity analyses, including Cochran’s Q test, MR-Egger intercept test, and MR-PRESSO global test. While heterogeneity was observed in the Q test analysis for COVID-19 (hospitalized vs. normal population) on Headache (Table 2), no heterogeneity was detected in other analyses. Notably, the MR estimates were not invalidated by this heterogeneity as random-effect IVW was used in this study, which balanced the pooled heterogeneity. Furthermore, while there were potential outliers in some of the IVs, such as COVID-19 (hospitalized vs. normal population) on Back pain, Headache, Neck or shoulder pain, and Pain all over the body, and COVID-19 (hospitalized vs. not hospitalized) on Neck or shoulder pain and Pain all over the body, as well as COVID-19 (very severe respiratory confirmed vs. normal population) on Pain all over the body, further MR-Egger intercept test and MR-PRESSO analysis did not identify any significant outliers (Table 2, S2S4 Tables). Therefore, there was insufficient evidence to suggest the presence of horizontal pleiotropy in the association between COVID-19 conditions and pain. Additionally, we ensured the eligibility of the SNPs included in the analysis, as they did not show any association with pain or related confounders in Phenoscanner and had F-values greater than 10, indicating the absence of weak instruments (S5 Table). Additionally, the leave-one-out analysis revealed that none of the individual independent IVs exerted an exaggerated influence on the causal effect (Fig 5A–5G). These sensitivity analyses provide further support for the robustness and validity of our findings.

thumbnail
Fig 5. Leave-one-out sensitivity analysis displaying the MR analysis for different conditions with COVID-19 on pain risk and the effect of excluding individual SNPs.

COVID-19 (hospitalized vs. normal population) on (A) Pain in joint(Lower leg), (B) Back pain, (C) Headache, (D) Pain all over the body; COVID-19 (hospitalized vs. not hospitalized) on (E) Neck or shoulder pain, (F) Pain all over the body; (G) COVID-19 (very severe respiratory confirmed vs. normal population) on Pain all over the body. MR, Mendelian randomization; SNP, single-nucleotide polymorphism.

https://doi.org/10.1371/journal.pone.0295982.g005

Discussion

Our study investigated the potential causal relationship between b different COVID-19 prevalence conditions and pain in various regions of the body. A two-sample MR analysis was used to investigate this relationship. IVW methods were used as the primary method to assess causal effects, revealing that hospitalised COVID-19 patients were at higher risk of developing Pain in joint (Lower leg) compared to the normal population. However, due to the stringent Bonferroni correction for multiple testing, the causal effects indicating an elevated risk of experiencing back pain, headache, and pain all over the body among hospitalized COVID-19 patients when compared to the general population were suggestive in nature. In addition, patients with COVID-19 who required hospitalisation were at a suggestive higher risk of developing neck or shoulder pain and pain all over the body compared to those who did not require hospitalisation. Patients with a severe respiratory confirmation of COVID-19 also had a suggestive increased risk of pain all over the body compared to the general population. Consistent direction was obtained from MR Egger and weighted median analyses, providing further support for the robustness and validity of the study results. To ensure the validity of the results, a series of sensitivity analyses were conducted. While heterogeneity was found in the Q-test analysis of the COVID-19 regarding headache, no heterogeneity was found in the other analyses. The MR estimates were not invalidated by this heterogeneity as the random effects IVW model was used in this study. In addition, the MR-Egger intercept test and MR-PRESSO analysis did not reveal any significant outliers, indicating the absence of pleiotropy.

Our findings suggest that the risk of pain in the overall COVID-19 population is not higher than that in the general population. This may be due to the prevalence of asymptomatic or minimally symptomatic cases of SARS-CoV-2 [30], which could dilute the overall effect of COVID-19. However, our results indicate a possible association between pain and COVID-19 in specific prevalence groups. For instance, hospitalised COVID-19 patients are at a higher risk of developing various types of pain, such as back pain and headache. According to Istvan et al [31]., the frequency of web search terms associated with these pain increased after the COVID-19 outbreak, which is consistent with our findings.

Headache is a common neurological symptom of long-COVID syndrome [32]. Cross-sectional studies reveal that long-COVID headaches are mainly bilateral, frontal or periocular, and compressive [33, 34]. Tension-like headache and migraine-like headache are the most common phenotypes, and refractory headache in the acute phase of infection is a risk factor associated with long-COVID headache [35]. Our results indicate that hospitalised COVID-19 patients are at a suggestive higher risk of developing headaches compared to the general population, suggesting that the risk of prolonged pain may increase with the severity of COVID-19. These findings are supported by evidence suggesting that acute infections can trigger chronic pain, and the severity and duration of the infection must be sufficient to disrupt normal activity [36, 37].

In addition, our results show that hospitalized COVID-19 patients may exhibit a heightened risk of experiencing pain all over the body compared to the general population. Similarly, a comparison between hospitalized and non-hospitalized COVID-19 patients also revealed an elevated risk of pain all over the body in the former group. Furthermore, individuals afflicted with very severe respiratory confirmed COVID-19 may also display an increased likelihood of developing pain all over the body when compared to the general population. This suggests that the presence of body pain may be a prevalent symptom of long-COVID manifestations in severely ill COVID-19 patients. Notably, hospitalized COVID-19 patients appear to be at an elevated risk of experiencing joint pain, particularly in the lower leg, compared to the general population. Additional exposure outcome datasets suggest that knee pain may also be increased, indicating that the knee may be a significant site of attack for long-COVID pain. Several longitudinal studies in Turkey, France, and Italy followed COVID-19 patients for six months after hospital discharge, and approximately 60% of patients continued to experience at least one associated symptom, most commonly fatigue, myalgia, and arthralgia. Arthralgia and myalgia were widespread throughout the body in the majority of patients [38]. A study by Sansin et al. found that myalgia, arthralgia, and back pain were prevalent in patients with COVID, with the knee being one of the primary sites of arthralgia [39]. Furthermore, there was a significant correlation between arthralgia and disease severity. There have been reports of early-onset arthritis in a few patients after severe acute infection, with this complication more frequently occurring in men and more likely to manifest in the lower leg [40, 41].

Chronic pain caused by COVID-19 can have a significant impact on a person’s quality of life, making it imperative to investigate this phenomenon. Such investigations may lead to a better understanding of the pathophysiology of COVID-19, which, in turn, could inform the development of new treatments and interventions. Furthermore, an examination of chronic pain may aid healthcare professionals in identifying patients who are at a higher risk of developing this condition. In our study, we found that hospitalized and severely respiratory confirmed COVID-19 patients were among those at high risk. Early identification of those at risk could prompt timely intervention and management, potentially preventing chronic pain from developing [42]. One potential mechanism through which COVID-19 may contribute to chronic pain is the development of long-COVID. long-COVID is a condition in which individuals experience symptoms of COVID-19 long after they have recovered from the acute phase of the disease [43]. Some individuals with long-COVID report persistent pain or discomfort months after their initial diagnosis of COVID-19. The pain they report is typically diffuse and can affect multiple areas of the body, including the back, neck, and joints [44], which is consistent with our results. Although the relationship between COVID-19 and chronic pain is not yet fully understood, there is evidence to suggest that the mechanisms involved may be multifactorial, involving changes in the immune system, central nervous system, and psychological factors.

Activation of the immune system is a potential mechanism underlying the development of chronic pain. Inflammation is a natural response of the immune system to injury or infection; however, excessive or chronic inflammation can result in tissue damage and organ dysfunction. Inflammation is a hallmark of COVID-19 and plays a pivotal role in the disease’s pathogenesis [14, 30]. The cytokine storm, a state of excessive immune response that results in widespread inflammation throughout the body, including muscles, joints, nerves, and the brain, is a characteristic feature of COVID-19 [45]. In COVID-19, the immune system overreacts, leading to the release of large quantities of pro-inflammatory cytokines such as interleukin-6, tumour necrosis factor-alpha, and interferon-gamma [14]. These cytokines can activate other immune cells and recruit them to the infection site, resulting in further inflammation. The sensitization of nerve fibres in the affected region by the release of chemicals from immune cells can lead to an increased response to stimuli, resulting in pain. Inflammation has the potential to activate glial cells within the central nervous system, resulting in increased sensitivity of nerves and intensified pain [46]. Futhermore, multiple studies have indicated that the SARS-CoV-2 virus is capable of infiltrating the central nervous system, potentially resulting in pain [47]. The entry of SARS-CoV-2 into human cells is facilitated by the binding of the spike protein to the angiotensin-converting enzyme 2 (ACE2) receptor, which is distributed across various tissues and organs. This interaction is a crucial pathophysiological mechanism underlying the development of COVID-19 [48]. The binding of SARS-CoV-2 to ACE2 and the consequent disruption of the renin-angiotensin system have been demonstrated to be implicated in endothelial and adipose dysfunction, as well as the initiation of neuronal sensitisation. Such mechanisms may lead to joint pain and osteoarthritis-like changes in bones [38]. Meanwhile, SARS-CoV-2 may also conceivably infect ACE2-expressing cells in the human spinal cord dorsal horn, leading to pain [49, 50]. The development of chronic pain in long-COVID may also involve psychological factors. The stress and anxiety brought about by the pandemic, coupled with social isolation during lockdowns, can trigger a physiological response that worsens pain [15, 16].

In this study, we employed the MR approach to examine the potential association between chronic pain and COVID-19. This method offers several advantages over retrospective and randomised controlled trials [51]. MR studies utilise genetic markers as IVs, which are randomly allocated at the time of conception and remain unaffected by lifestyle or environmental factors [21]. Therefore, the MR approach can significantly reduce the bias that is typically present in retrospective and randomised controlled trials [20]. Moreover, the MR method can utilise large-scale genetic data to improve statistical power, while the sample sizes required for retrospective and randomised controlled trials can often be prohibitively large [24]. Additionally, MR studies are often more cost-effective than randomised controlled or retrospective studies, as genetic data can be obtained from pre-existing databases or biobanks, and the data collection process is less expensive than that of conducting a randomised controlled trial [24, 52].

Our study has several limitations that warrant consideration. Firstly, we conducted a two-sample MR analysis, assuming the genetic tools used to represent COVID-19 exposure were valid and not influenced by polymorphism. While we performed sensitivity analyses to assess the validity of the results, unmeasured confounders could have influenced the observed associations [24, 51]. Secondly, the study’s focus on COVID-19 patients from Western populations limits the generalizability of the findings to other populations with distinct genetic and environmental factors. Thirdly, we did not examine the duration of COVID-19 symptoms and its potential impact on the risk of developing long-term pain. Future studies should investigate the temporal relationship between COVID-19 exposure and pain outcomes. Moreover, while meticulous corrections for multiple testing have been diligently applied, our investigation has illuminated noteworthy associations between specific COVID-19 conditions and chronic pain. Nevertheless, it is essential to acknowledge that the effect sizes derived from the MR analysis are of a relatively modest magnitude. Consequently, we advise circumspection when interpreting these particular conclusions, emphasizing the imperative need for subsequent comprehensive research endeavors and inquiries to corroborate and substantiate these findings. Finally, the availability of data for certain pain outcomes restricted the scope of our analysis. Future studies must investigate additional pain sites and different levels of pain intensity. Confirming the associations observed in this study will require larger epidemiological investigations.

Conclusion

In conclusion, the detrimental impact of chronic pain on an individual’s quality of life underscores the importance of studying its prevalence in the context of COVID-19. Our research sheds light on the potential association between varying degrees of COVID-19 severity and the manifestation of pain across different regions of the body. These findings hold significant implications for clinical management of COVID-19 patients. By identifying those at high risk of developing chronic pain, targeted interventions can be developed to mitigate the onset and progression of this debilitating condition, thereby reducing the burden on individuals, healthcare systems, and society at large.

Supporting information

S1 Table. Detailed information pertaining to SNPs in the two-sample MR analysis.

https://doi.org/10.1371/journal.pone.0295982.s001

(DOCX)

S2 Table. MR estimates of associations between COVID-19 (infection vs. normal population) and bodily pain in various regions across different methods.

https://doi.org/10.1371/journal.pone.0295982.s002

(DOCX)

S3 Table. MR estimates of associations between COVID-19 (hospitalized vs. not hospitalized) and bodily pain in various regions across different methods.

https://doi.org/10.1371/journal.pone.0295982.s003

(DOCX)

S4 Table. MR estimates of associations between COVID-19 (very severe respiratory confirmed vs. normal population) and bodily pain in various regions across different methods.

https://doi.org/10.1371/journal.pone.0295982.s004

(DOCX)

S5 Table. Details of SNPs extracted from exposure and F of them.

https://doi.org/10.1371/journal.pone.0295982.s005

(DOCX)

References

  1. 1. Mizrahi B, Sudry T, Flaks-Manov N, Yehezkelli Y, Kalkstein N, et al. (2023) long-COVID outcomes at one year after mild SARS-CoV-2 infection: nationwide cohort study. BMJ 380: e072529. pmid:36631153
  2. 2. Li J, Lai S, Gao GF, Shi W (2021) The emergence, genomic diversity and global spread of SARS-CoV-2. Nature 600: 408–418. pmid:34880490
  3. 3. Struyf T, Deeks JJ, Dinnes J, Takwoingi Y, Davenport C, et al. (2020) Signs and symptoms to determine if a patient presenting in primary care or hospital outpatient settings has COVID-19 disease. Cochrane Database Syst Rev 7: CD013665. pmid:32633856
  4. 4. Chow EJ, Schwartz NG, Tobolowsky FA, Zacks R, Huntington-Frazier M, et al. (2020) Symptom Screening at Illness Onset of Health Care Personnel With SARS-CoV-2 Infection in King County, Washington. JAMA 323: 2087–2089. pmid:32301962
  5. 5. Gentile S, Strollo F, Ceriello A (2020) COVID-19 infection in Italian people with diabetes: Lessons learned for our future (an experience to be used). Diabetes Res Clin Pract 162: 108137. pmid:32259611
  6. 6. Hu X, Hu C, Yang Y, Chen J, Zhong P, et al. (2020) Clinical characteristics and risk factors for severity of COVID-19 outside Wuhan: a double-center retrospective cohort study of 213 cases in Hunan, China. Ther Adv Respir Dis 14: 1753466620963035. pmid:33138694
  7. 7. Yang J, Hu J, Zhu C (2021) Obesity aggravates COVID-19: A systematic review and meta-analysis. J Med Virol 93: 257–261. pmid:32603481
  8. 8. Callard F, Perego E (2021) How and why patients made long-COVID. Soc Sci Med 268: 113426. pmid:33199035
  9. 9. Herman E, Shih E, Cheng A (2022) long-COVID: Rapid Evidence Review. Am Fam Physician 106: 523–532. pmid:36379497
  10. 10. Bowyer R, Huggins C, Toms R, Shaw RJ, Hou B, et al. (2023) Characterising patterns of COVID-19 and long-COVID symptoms: evidence from nine UK longitudinal studies. Eur J Epidemiol 38: 199–210. pmid:36680646
  11. 11. Oktay Gültekin E, Gültekin O, Coskun A, Aksak T (2022) Antibody response three months after SARS-CoV-2 infection. J Med Virol 94: 4712–4718. pmid:35655391
  12. 12. Puntillo F, Giglio M, Brienza N, Viswanath O, Urits I, et al. (2020) Impact of COVID-19 pandemic on chronic pain management: Looking for the best way to deliver care. Best Pract Res Clin Anaesthesiol 34: 529–537. pmid:33004164
  13. 13. Clauw DJ, Häuser W, Cohen SP, Fitzcharles MA (2020) Considering the potential for an increase in chronic pain after the COVID-19 pandemic. Pain 161: 1694–1697. pmid:32701829
  14. 14. Chen G, Wu D, Guo W, Cao Y, Huang D, et al. (2020) Clinical and immunological features of severe and moderate coronavirus disease 2019. J Clin Invest 130: 2620–2629. pmid:32217835
  15. 15. Chaturvedi SK (2020) Health anxiety, health-related life events, and somatization during COVID-19 pandemic can increase chronic pain. Pain 161: 2652. pmid:33065698
  16. 16. Lee CC, Wu HJ, Chiou WB (2022) Reminders of COVID-19 social distancing can intensify physical pain. Br J Soc Psychol 61: 587–598. pmid:34510481
  17. 17. Petersen MS, Kristiansen MF, Hanusson KD, Danielsen ME, Á Steig B, et al. (2021) long-COVID in the Faroe Islands: A Longitudinal Study Among Nonhospitalized Patients. Clin Infect Dis 73: e4058–4058e4063. pmid:33252665
  18. 18. Jacobs LG, Gourna Paleoudis E, Lesky-Di Bari D, Nyirenda T, Friedman T, et al. (2020) Persistence of symptoms and quality of life at 35 days after hospitalization for COVID-19 infection. PLoS One 15: e0243882. pmid:33306721
  19. 19. Carfì A, Bernabei R, Landi F, Gemelli Against COVID-19 Post-Acute Care Study Group (2020) Persistent Symptoms in Patients After Acute COVID-19. JAMA 324: 603–605.
  20. 20. Davey Smith G, Hemani G (2014) Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet 23: R89–98. pmid:25064373
  21. 21. Skrivankova VW, Richmond RC, Woolf B, Yarmolinsky J, Davies NM, et al. (2021) Strengthening the Reporting of Observational Studies in Epidemiology Using Mendelian Randomization: The STROBE-MR Statement. JAMA 326: 1614–1621. pmid:34698778
  22. 22. Sanderson E (2021) Multivariable Mendelian Randomization and Mediation. Cold Spring Harb Perspect Med 11. pmid:32341063
  23. 23. Sekula P, Del Greco M F, Pattaro C, Köttgen A (2016) Mendelian Randomization as an Approach to Assess Causality Using Observational Data. J Am Soc Nephrol 27: 3253–3265. pmid:27486138
  24. 24. Richmond RC, Davey Smith G (2022) Mendelian Randomization: Concepts and Scope. Cold Spring Harb Perspect Med 12. pmid:34426474
  25. 25. Sanna S, van Zuydam NR, Mahajan A, Kurilshikov A, Vich Vila A, et al. (2019) Causal relationships among the gut microbiome, short-chain fatty acids and metabolic diseases. Nat Genet 51: 600–605. pmid:30778224
  26. 26. Burgess S, Thompson SG, CRP CHD Genetics Collaboration (2011) Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol 40: 755–764.
  27. 27. Bowden J, Davey Smith G, Burgess S (2015) Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol 44: 512–525. pmid:26050253
  28. 28. Burgess S, Thompson SG (2017) Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol 32: 377–389. pmid:28527048
  29. 29. Verbanck M, Chen CY, Neale B, Do R (2018) Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet 50: 693–698. pmid:29686387
  30. 30. Sharma A, Ahmad Farouk I, Lal SK (2021) COVID-19: A Review on the Novel Coronavirus Disease Evolution, Transmission, Detection, Control and Prevention. Viruses 13. pmid:33572857
  31. 31. Szilagyi IS, Ullrich T, Lang-Illievich K, Klivinyi C, Schittek GA, et al. (2021) Google Trends for Pain Search Terms in the World’s Most Populated Regions Before and After the First Recorded COVID-19 Case: Infodemiological Study. J Med Internet Res 23: e27214.
  32. 32. Sudre CH, Murray B, Varsavsky T, Graham MS, Penfold RS, et al. (2021) Attributes and predictors of long-COVID. Nat Med 27: 626–631. pmid:33692530
  33. 33. Caronna E, Ballvé A, Llauradó A, Gallardo VJ, Ariton DM, et al. (2020) Headache: A striking prodromal and persistent symptom, predictive of COVID-19 clinical evolution. Cephalalgia 40: 1410–1421. pmid:33146036
  34. 34. Caronna E, Alpuente A, Torres-Ferrus M, Pozo-Rosich P (2021) Toward a better understanding of persistent headache after mild COVID-19: Three migraine-like yet distinct scenarios. Headache 61: 1277–1280. pmid:34363619
  35. 35. Planchuelo-Gómez Á, Trigo J, de Luis-García R, Guerrero ÁL, Porta-Etessam J, García-Azorín D (2020) Deep Phenotyping of Headache in Hospitalized COVID-19 Patients via Principal Component Analysis. Front Neurol 11: 583870. pmid:33391152
  36. 36. Hickie I, Davenport T, Wakefield D, Vollmer-Conna U, Cameron B, et al. (2006) Post-infective and chronic fatigue syndromes precipitated by viral and non-viral pathogens: prospective cohort study. BMJ 333: 575. pmid:16950834
  37. 37. Wessely S, Chalder T, Hirsch S, Pawlikowska T, Wallace P, Wright DJ (1995) Postinfectious fatigue: prospective cohort study in primary care. Lancet 345: 1333–1338. pmid:7752755
  38. 38. Lauwers M, Au M, Yuan S, Wen C (2022) COVID-19 in Joint Ageing and Osteoarthritis: Current Status and Perspectives. Int J Mol Sci 23. pmid:35054906
  39. 39. Tuzun S, Keles A, Okutan D, Yildiran T, Palamar D (2021) Assessment of musculoskeletal pain, fatigue and grip strength in hospitalized patients with COVID-19. Eur J Phys Rehabil Med 57: 653–662. pmid:33393277
  40. 40. Gasparotto M, Framba V, Piovella C, Doria A, Iaccarino L (2021) Post-COVID-19 arthritis: a case report and literature review. Clin Rheumatol 40: 3357–3362. pmid:33587197
  41. 41. Kocyigit BF, Akyol A (2021) Reactive arthritis after COVID-19: a case-based review. Rheumatol Int 41: 2031–2039. pmid:34550429
  42. 42. Shanthanna H, Nelson AM, Kissoon N, Narouze S (2022) The COVID-19 pandemic and its consequences for chronic pain: a narrative review. Anaesthesia 77: 1039–1050. pmid:35848380
  43. 43. Crook H, Raza S, Nowell J, Young M, Edison P (2021) long-COVID-mechanisms, risk factors, and management. BMJ 374: n1648. pmid:34312178
  44. 44. Aiyegbusi OL, Hughes SE, Turner G, Rivera SC, McMullan C, et al. (2021) Symptoms, complications and management of long-COVID: a review. J R Soc Med 114: 428–442. pmid:34265229
  45. 45. Nalbandian A, Sehgal K, Gupta A, Madhavan MV, McGroder C, et al. (2021) Post-acute COVID-19 syndrome. Nat Med 27: 601–615. pmid:33753937
  46. 46. Ji RR, Chamessian A, Zhang YQ (2016) Pain regulation by non-neuronal cells and inflammation. Science 354: 572–577. pmid:27811267
  47. 47. Wu Y, Xu X, Chen Z, Duan J, Hashimoto K, et al. (2020) Nervous system involvement after infection with COVID-19 and other coronaviruses. Brain Behav Immun 87: 18–22. pmid:32240762
  48. 48. Jackson CB, Farzan M, Chen B, Choe H (2022) Mechanisms of SARS-CoV-2 entry into cells. Nat Rev Mol Cell Biol 23: 3–20. pmid:34611326
  49. 49. Yamagata R, Nemoto W, Nakagawasai O, Takahashi K, Tan-No K (2020) Downregulation of spinal angiotensin converting enzyme 2 is involved in neuropathic pain associated with type 2 diabetes mellitus in mice. Biochem Pharmacol 174: 113825. pmid:31987854
  50. 50. Nemoto W, Nakagawasai O, Yaoita F, Kanno S, Yomogida S, et al. (2013) Angiotensin II produces nociceptive behavior through spinal AT1 receptor-mediated p38 mitogen-activated protein kinase activation in mice. Mol Pain 9: 38. pmid:23898828
  51. 51. Emdin CA, Khera AV, Kathiresan S (2017) Mendelian Randomization. JAMA 318: 1925–1926. pmid:29164242
  52. 52. Zheng J, Baird D, Borges MC, Bowden J, Hemani G, et al. (2017) Recent Developments in Mendelian Randomization Studies. Curr Epidemiol Rep 4: 330–345. pmid:29226067