Sex- and age-specific genetic analysis of chronic back pain

Gene discovery and pleiotropic architecture of Chronic Pain in a Genome-wide Association Study of >1.2 million Individuals

Chronic pain is highly prevalent worldwide, and genome-wide association studies (GWAS) have identified a growing number of chronic pain loci. To further elucidate its genetic architecture, we leveraged data from 1,235,695 European ancestry individuals across three biobanks. In a meta-analytic GWAS, we identified 343 independent loci for chronic pain, 92 of which were new. Sex-specific meta-analyses revealed 115 independent loci (12 of which were new) for males (N = 583,066) and 12 loci (two of which were new) for females (N = 241,266). Multi-omics gene prioritization analyses highlighted 490 genes associated with chronic pain through their effects on brain- and blood-specific regulation. Loci associated with increased risk for chronic pain were also associated with increased risk for multiple other traits, with Mendelian randomization analyses showing that chronic pain was causally associated with psychiatric disorders, substance use disorders, and C-reactive protein levels. Chronic pain variants also exhibited pleiotropic associations with cortical area brain structures. This study expands our knowledge of the genetics of chronic pain and its pathogenesis, highlighting the importance of its pleiotropy with multiple disorders and elucidating its multi-omic pathophysiology.

Characterizing Gabapentinoid Use Among Patients With Isolated Low Back Pain

STUDY DESIGN

Retrospective cohort study.

OBJECTIVES

Gabapentinoids have been increasingly studied as a non-narcotic option for neuropathic and postoperative pain. However, there is evidence suggesting that off-label use of these medications for the treatment of isolated LBP is not effective. The purpose of the current study was to evaluate prescription patterns for gabapentinoids among patients with isolated low back pain.

METHODS

Adult patients with LBP were abstracted from the 2019 M91Ortho PearlDiver dataset using International Classification of Diseases (ICD-10) code M54.5. Patients were excluded if they had a diagnosis of neurologic symptoms, history of spinal surgery, spinal fracture, or conditions for which gabapentinoids are FDA approved. Gabapentinoid and narcotic prescriptions within one year of LBP diagnosis were identified. Patient characteristics and prescriber specialty were extracted from the dataset and predictors of gabapentinoid prescriptions were determined using univariate and multivariate analyses.

RESULTS

Among the 1,158,875 isolated LBP patients, gabapentinoids were prescribed for 11%, narcotics for 8%, and both for 3%. The most common prescriber specialties included: primary care physicians (45%), nurse practitioners (15%), and pain management physicians (5%). Independent predictors of gabapentinoid prescriptions included: female sex, region of the country, and insurance type (P-values <.001).

CONCLUSIONS

Of nearly 1.2 million patients with isolated LBP, 14.2% were prescribed gabapentinoids within one year of their LBP diagnosis. Patient characteristics predictive of receiving gabapentinoids were predominantly non-clinical factors.

Multi-ancestry meta-analysis of genome-wide association studies discovers 67 new loci associated with chronic back pain

This multi-ancestry meta-analysis of genome-wide association studies (GWAS) investigated the genetic factors underlying chronic back pain (CBP) in a sample from the Million Veteran Program comprised of 553,601 Veterans of African (19.2%), European (72.6%), and Hispanic (8.2%) ancestry. The results revealed novel (N = 67) and known (N = 20) genome-wide significant loci associated with CBP, with 43 independent variants replicating in a non-overlapping contemporary meta-GWAS of the spinal pain dorsalgia phenotype. The most significant novel variant was rs12533005 (chr7:114416000, p = 1.61 × 10-20, OR = 0.96 (95% CI: 0.95-0.97), EA = C, EAF = 0.39), in an intron of the FOXP2 gene. In silico functional characterization revealed enrichment in brain and pituitary tissues. Mendelian randomization analysis of 62 variants for CBP-MVP revealed 48 with causal links to dorsalgia. Notably, four genes (INPP5B, DRD2, HTT, SLC30A6) associated with these variants are targets of existing drugs. Our findings more than double the number of previously reported genetic predictors across all spinal pain phenotypes.

Gene Expression Correlates with Disability and Pain Intensity in Patients with Chronic Low Back Pain and Modic Changes in a Sex-Specific Manner

Chronic low back pain (cLBP) lacks clear physiological explanations, and the treatment options are of limited effect. We aimed to elucidate the underlying biology of cLBP in a subgroup of patients with Modic changes type I (suggestive of inflammatory vertebral bone marrow lesions) by correlating gene expression in blood with patient-reported outcomes on disability and pain intensity and explore sex differences. Patients were included from the placebo group of a clinical study on patients with cLBP and Modic changes. Blood was collected at the time of inclusion, after three months, and after one year, and gene expression was measured at all time points by high-throughput RNA sequencing. The patients reported disability using the Roland-Morris Disability Questionnaire, and pain intensity was assessed as a mean of three scores on a 0-10 numeric rating scale: current LBP, worst LBP within the last two weeks, and mean LBP within the last two weeks. The gene expression profiles were then correlated to the reported outcomes. Changes in gene expression over time correlated significantly with changes in both disability and pain. The findings showed distinct patterns in men and women, with negligible overlap in correlated genes between the sexes. The genes involved were enriched in immunological pathways, particularly T cell receptor complex and immune responses related to neutrophils. Several of the genes harbour polymorphisms that previously have been found to be associated with chronic pain. Taken together, our results indicate gender differences in the underlying biology of disability and pain intensity in patients with low back pain.

Trends of Korean medicine service utilization for lumbar disc herniation and spinal stenosis: A 10-year analysis of the 2010 to 2019 data

This retrospective, cross-sectional, and descriptive study aimed to analyze the trend of utilizing traditional Korean medicine services in patients with lumbar disc herniation (LDH) and/or lumbar spinal stenosis (LSS). In this study, based on the national patient sample data provided by Health Insurance Review and Assessment Service (HIRA), the trend of Korean medicine service utilization was investigated, including the following information: demographic characteristics of the patients, the total expenditure, number of claim statements per category, medical care expenditure per category, and routes of visiting traditional Korean medicine institutions. The study population comprised patients who visited Korean medicine institutions at least once from January 2010 to December 2019, with LSS and LDH as the primary diagnosis. LDH patients who used traditional Korean medicine services for treatment increased by about 1.36 times. LDH and LSS patients under 45 years of age were more likely to be males, but women accounted for a higher percentage among those over 45 years of age. Overall, women accounted for a slightly higher percentage than their counterparts for both diseases. From details of treatments received that were extracted from the claims data, acupuncture treatment accounted for the highest percentage for both disorders. Moreover, 50.7% of the patients who visited Korean medicine institutions to treat the two diseases also visited conventional Western medicine institutions. These patients, who were diagnosed with their condition at a Korean medicine institution, visited a conventional institution and then returned; the conventional institutions were primarily used for examination (40.5%). Increased utilization of traditional Korean medicine services was confirmed among patients with LDH and/or LSS; in particular, a sharp increase was noted among patients with LSS. The results of this study will be useful as basic research data for clinicians, researchers, and policy makers.

A systematic review of genome-wide association studies for pain, nociception, neuropathy, and pain treatment responses

ABSTRACT

Pain is the leading cause of disability worldwide, imposing an enormous burden on personal health and society. Pain is a multifactorial and multidimensional problem. Currently, there is (some) evidence that genetic factors could partially explain individual susceptibility to pain and interpersonal differences in pain treatment response. To better understand the underlying genetic mechanisms of pain, we systematically reviewed and summarized genome-wide association studies (GWASes) investigating the associations between genetic variants and pain/pain-related phenotypes in humans. We reviewed 57 full-text articles and identified 30 loci reported in more than 1 study. To check whether genes described in this review are associated with (other) pain phenotypes, we searched 2 pain genetic databases, Human Pain Genetics Database and Mouse Pain Genetics Database. Six GWAS-identified genes/loci were also reported in those databases, mainly involved in neurological functions and inflammation. These findings demonstrate an important contribution of genetic factors to the risk of pain and pain-related phenotypes. However, replication studies with consistent phenotype definitions and sufficient statistical power are required to validate these pain-associated genes further. Our review also highlights the need for bioinformatic tools to elucidate the function of identified genes/loci. We believe that a better understanding of the genetic background of pain will shed light on the underlying biological mechanisms of pain and benefit patients by improving the clinical management of pain.

Big data, big consortia, and pain: UK Biobank, PAINSTORM, and DOLORisk

Chronic pain (CP) is a common and often debilitating disorder that has major social and economic impacts. A subset of patients develop CP that significantly interferes with their activities of daily living and requires a high level of healthcare support. The challenge for treating physicians is in preventing the onset of refractory CP or effectively managing existing pain. To be able to do this, it is necessary to understand the risk factors, both genetic and environmental, for the onset of CP and response to treatment, as well as the pathogenesis of the disorder, which is highly heterogenous. However, studies of CP, particularly pain with neuropathic characteristics, have been hindered by a lack of consensus on phenotyping and data collection, making comparisons difficult. Furthermore, existing cohorts have suffered from small sample sizes meaning that analyses, especially genome-wide association studies, are insufficiently powered. The key to overcoming these issues is through the creation of large consortia such as DOLORisk and PAINSTORM and biorepositories, such as UK Biobank, where a common approach can be taken to CP phenotyping, which allows harmonisation across different cohorts and in turn increased study power. This review describes the approach that was used for studying neuropathic pain in DOLORisk and how this has informed current projects such as PAINSTORM, the rephenotyping of UK Biobank, and other endeavours. Moreover, an overview is provided of the outputs from these studies and the lessons learnt for future projects.

Genotype-by-environment interactions in chronic back pain

BACKGROUND CONTEXT

Chronic back pain (CBP) is a common debilitating condition with substantial societal impact. While understanding genotype-by-environment (GxE) interactions may be crucial to achieving the goals of personalized medicine, there are few large-scale studies investigating this topic for CBP. None of them systematically explore multiple CBP risk factors.

PURPOSE

To estimate the extent to which genetic effects on CBP are modified by known demographic and clinical risk factors.

RESEARCH DESIGN

Case-control study, genome-wide GxE interaction study.

PATIENT SAMPLE

Data on up to 331,610 unrelated participants (57,881 CBP cases and 273,729 controls) from the UK Biobank cohort were used. UK Biobank is a prospective cohort with collected deep genetic and phenotypic data on approximately 500,000 individuals across the UK.

OUTCOME MEASURES

Self-reported chronic back pain.

METHODS

We applied a whole-genome approach to estimate the proportion of phenotypic variance explained by interactions between genotype and 12 known risk factors. We also analyzed if effects of common single-nucleotide polymorphisms on CBP are changed in presence of known risk factors.

RESULTS

The results indicate a modest, if any, modification of genetic effects by examined risk factors in CBP. Our estimates suggest that detecting such weak effects would require a sample size of millions of individuals.

CONCLUSIONS

The GxE interactions with examined common risk factors for CBP are either weak or absent. Interactions of such magnitude are unlikely to have the potential to inform and influence treatment strategies. Risk estimation models may use common genetic variation and the considered risk factors as independent predictors, without accounting for GxE.

Development and Replication of a Genome-Wide Polygenic Risk Score for Chronic Back Pain

Chronic back pain (CBP) is a complex heritable trait and a major cause of disability worldwide. We developed and validated a genome-wide polygenic risk score (PRS) for CBP using a large-scale GWAS based on UK Biobank participants of European ancestry (N = 265,000). The PRS showed poor overall predictive ability (AUC = 0.56 and OR = 1.24 per SD, 95% CI: 1.22-1.26), but individuals from the 99th percentile of PRS distribution had a nearly two-fold increased risk of CBP (OR = 1.82, 95% CI: 1.60-2.06). We validated the PRS on an independent TwinsUK sample, obtaining a similar magnitude of effect. The PRS was significantly associated with various ICD-10 and OPCS-4 diagnostic codes, including chronic ischemic heart disease (OR = 1.1, p-value = 4.8 × 10^-15), obesity, metabolism-related traits, spine disorders, disc degeneration, and arthritis-related disorders. PRS and environment interaction analysis with twelve known CBP risk factors revealed no significant results, suggesting that the magnitude of G × E interactions with studied factors is small. The limited predictive ability of the PRS that we developed is likely explained by the complexity, heterogeneity, and polygenicity of CBP, for which sample sizes of a few hundred thousand are insufficient to estimate small genetic effects robustly.

A Shared Genetic Signature for Common Chronic Pain Conditions and its Impact on Biopsychosocial Traits

The multiple comorbidities & dimensions of chronic pain present a formidable challenge in disentangling its aetiology. Here, we performed genome-wide association studies of 8 chronic pain types using UK Biobank data (N =4,037-79,089 cases; N = 239,125 controls), followed by bivariate linkage disequilibrium-score regression and latent causal variable analyses to determine (respectively) their genetic correlations and genetic causal proportion (GCP) parameters with 1,492 other complex traits. We report evidence of a shared genetic signature across chronic pain types as their genetic correlations and GCP directions were broadly consistent across an array of biopsychosocial traits. Across 5,942 significant genetic correlations, 570 trait pairs could be explained by a causal association (|GCP| >0.6; 5% false discovery rate), including 82 traits affected by pain while 410 contributed to an increased risk of chronic pain (cf. 78 with a decreased risk) such as certain somatic pathologies (eg, musculoskeletal), psychiatric traits (eg, depression), socioeconomic factors (eg, occupation) and medical comorbidities (eg, cardiovascular disease). This data-driven phenome-wide association analysis has demonstrated a novel and efficient strategy for identifying genetically supported risk & protective traits to enhance the design of interventional trials targeting underlying causal factors and accelerate the development of more effective treatments with broader clinical utility. PERSPECTIVE: Through large-scale phenome-wide association analyses of >1,400 biopsychosocial traits, this article provides evidence for a shared genetic signature across 8 common chronic pain types. It lays the foundation for further translational studies focused on identifying causal genetic variants and pathophysiological pathways to develop novel diagnostic & therapeutic technologies and strategies.

Challenges and opportunities for omics-based precision medicine in chronic low back pain

PURPOSE

Chronic low back pain (cLBP) is a common health condition worldwide and a leading cause of disability with an estimated lifetime prevalence of 80-90% in industrialized countries. However, we have had limited success in treating cLBP likely due to its non-specific heterogeneous nature that goes beyond detectable anatomical changes. We propose that omics technologies as precision medicine tools are well suited to provide insight into its pathophysiology and provide diagnostic markers and therapeutic targets. Therefore, in this review, we explore the current state of omics technologies in the diagnosis and classification of cLBP. We identify factors that may serve as markers to differentiate between acute and chronic cases of low back pain (LBP). Finally, we also discuss some challenges that must be overcome to successfully apply precision medicine to the diagnosis and treatment of cLBP.

METHODS

A literature search for the current applications of omics technologies to chronic low back pain was performed using the following search terms- "back pain," "low back pain," "proteomics," "transcriptomics", "epigenomics," "genomics," "omics." We reviewed molecular markers identified from 35 studies which hold promise in providing information regarding molecular insights into cLBP.

RESULTS

GWAS studies have found evidence for the role of single nucleotide polymorphisms (SNPs) associated with pain pathways in individuals with cLBP. Epigenomic modifications in patients with cLBP have been found to be enriched among genes involved in immune signaling and inflammation. Transcriptomics profiles of patients with cLBP show multiple lines of evidence for the role of inflammation in cLBP. The glycomics profiles of patients with cLBP are similar to those of patients with inflammatory conditions. Proteomics and microbiomics show promise but have limited studies currently.

CONCLUSION

Omics technologies have identified associations between inflammatory and pain pathways in the pathophysiology of cLBP. However, in order to integrate information across the range of studies, it is important for the field to identify and adopt standardized definitions of cLBP and control patients. Additionally, most papers have applied a single omics method to a sampling of cLBP patients which have yielded limited insight into the pathophysiology of cLBP. Therefore, we recommend a multi-omics approach applied to large global consortia for advancing subphenotyping and better management of cLBP, via improved identification of diagnostic markers and therapeutic targets.

Age and Sex in Back Pain Intensity—Retrospective Study of Conservatory vs. Surgical Discopathy Treatment

(1) Background: The frequency of back pain diagnosis and treatment has markedly increased in recent years. Back pain may be caused by many factors and discopathy is one of them. The aim of the study was to assess the impact of age and sex on back pain intensity in rehabilitated patients with discopathy treated conservatively and surgically; (2) Methods: The study included 137 patients: 96 undergoing conservative therapy and 41 after back surgery due to discopathy. VAS and the Laitinen scale were used for pain assessment. All patients underwent a multidirectional rehabilitation program at the Department of Rehabilitation and Physical Medicine of the Medical University in Łódź; (3) Results: No statistically significant effect of age and sex was observed on the level of pain intensity on VAS and the Laitinen scale; (4) Conclusions: Age and sex do not seem to affect back pain intensity in rehabilitated patients treated conservatively and surgically for discopathy. The problem requires further research on a larger group of patients.

Brain-specific genes contribute to chronic but not to acute back pain

Introduction

Back pain is the leading cause of disability worldwide. Although most back pain cases are acute, 20% of acute pain patients experience chronic back pain symptoms. It is unclear whether acute pain and chronic pain have similar or distinct underlying genetic mechanisms.

Objectives

To characterize the molecular and cellular pathways contributing to acute and chronic pain states.

Methods

Cross-sectional observational genome-wide association study.

Results

A total of 375,158 individuals from the UK Biobank cohort were included in the discovery of genome-wide association study. Of those, 70,633 (19%) and 32,209 (9%) individuals met the definition of chronic and acute back pain, respectively. A total of 355 single nucleotide polymorphism grouped into 13 loci reached the genome-wide significance threshold (5x10^-8) for chronic back pain, but none for acute. Of these, 7 loci were replicated in the Nord-Trøndelag Health Study (HUNT) cohort (19,760 chronic low back pain cases and 28,674 pain-free controls). Single nucleotide polymorphism heritability was 4.6% (P=1.4x10^-78) for chronic back pain and 0.81% (P=1.4x10-8) for acute back pain. Similar differences in heritability estimates between acute and chronic back pain were found in the HUNT cohort: 3.4% (P=0.0011) and 0.6% (P=0.851), respectively. Pathway analyses, tissue-specific heritability enrichment analyses, and epigenetic characterization suggest a substantial genetic contribution to chronic but not acute back pain from the loci predominantly expressed in the central nervous system.

Conclusion

Chronic back pain is substantially more heritable than acute back pain. This heritability is mostly attributed to genes expressed in the brain.

Comorbidity and familial aggregation of back/neck pain in the NIMH Family Study of Affective Spectrum Disorders

OBJECTIVE

Back pain is associated with substantial Global Burden of Disease and is highly comorbid with mood and anxiety symptoms and syndromes. However, mechanisms underlying this association have not been well-elucidated. Here we apply data from the NIMH Family Study of Affective Spectrum Disorders to investigate the comorbidity, familial aggregation, and cross-aggregation of back/neck pain with mood disorder subtypes.

METHODS

The sample includes 519 probands and 560 interviewed first-degree relatives. Lifetime DSM-IV Bipolar I, Bipolar II, and Major Depressive Disorder [MDD] were derived from semi-structured diagnostic interviews. Lifetime history of back or neck pain and its age of onset were self-reported retrospectively. Familial aggregation and cross-aggregation were estimated via mixed effects models in probands and interviewed first-degree relatives, while heritability and co-heritability (endophenotypic ranking value [ERV]) were estimated using full pedigrees.

RESULTS

Over 45% of participants endorsed a history of back/neck pain. Back/neck pain was familial (adjusted odds ratio [aOR] 1.5, p = 0.04; h2 = 0.24, p = 0.009). Back/neck pain in probands was associated with MDD in relatives (aOR 1.5, p = 0.04; ERV = 0.17, p = 0.024), but not with bipolar disorder. Onset of back/neck pain occurred earlier in those with bipolar disorder compared to controls.

CONCLUSION

Findings suggest common familial risk factors underlying back/neck pain with MDD, whereas there was within-individual comorbidity of bipolar with back/neck pain. Future studies that identify common factors that lead to either back/neck pain or MDD can inform prevention and interventions.

A protocol for recruiting and analyzing the disease-oriented Russian disc degeneration study (RuDDS) biobank for functional omics studies of lumbar disc degeneration

Lumbar intervertebral disc degeneration (DD) disease is one of the main risk factors for low back pain and a leading cause of population absenteeism and disability worldwide. Despite a variety of biological studies, lumbar DD is not yet fully understood, partially because there are only few studies that use systematic and integrative approaches. This urges the need for studies that integrate different omics (including genomics and transcriptomics) measured on samples within a single cohort. This protocol describes a disease-oriented Russian disc degeneration study (RuDDS) biobank recruitment and analyses aimed to facilitate further omics studies of lumbar DD integrating genomic, transcriptomic and glycomic data. A total of 1,100 participants aged over 18 with available lumbar MRI scans, medical histories and biological material (whole blood, plasma and intervertebral disc tissue samples from surgically treated patients) will be enrolled during the three-year period from two Russian clinical centers. Whole blood, plasma and disc tissue specimens will be used for genotyping with genome-wide SNP-arrays, glycome profiling and RNA sequencing, respectively. Omics data will be further used for a genome-wide association study of lumbar DD with in silico functional annotation, analysis of plasma glycome and lumbar DD disease interactions and transcriptomic data analysis including an investigation of differential expression patterns associated with lumbar DD disease. Statistical tests applied in each of the analyses will meet the standard criteria specific to the attributed study field. In a long term, the results of the study will expand fundamental knowledge about lumbar DD development and contribute to the elaboration of novel personalized approaches for disease prediction and therapy. Additionally to the lumbar disc degeneration study, a RuDDS cohort could be used for other genetic studies, as it will have unique omics data.

Trial registration numberNCT04600544.

Acute inflammatory response via neutrophil activation protects against the development of chronic pain

The transition from acute to chronic pain is critically important but not well understood. Here, we investigated the pathophysiological mechanisms underlying the transition from acute to chronic low back pain (LBP) and performed transcriptome-wide analysis in peripheral immune cells of 98 participants with acute LBP, followed for 3 months. Transcriptomic changes were compared between patients whose LBP was resolved at 3 months with those whose LBP persisted. We found thousands of dynamic transcriptional changes over 3 months in LBP participants with resolved pain but none in those with persistent pain. Transient neutrophil-driven up-regulation of inflammatory responses was protective against the transition to chronic pain. In mouse pain assays, early treatment with a steroid or nonsteroidal anti-inflammatory drug (NSAID) also led to prolonged pain despite being analgesic in the short term; such a prolongation was not observed with other analgesics. Depletion of neutrophils delayed resolution of pain in mice, whereas peripheral injection of neutrophils themselves, or S100A8/A9 proteins normally released by neutrophils, prevented the development of long-lasting pain induced by an anti-inflammatory drug. Analysis of pain trajectories of human subjects reporting acute back pain in the UK Biobank identified elevated risk of pain persistence for subjects taking NSAIDs. Thus, despite analgesic efficacy at early time points, the management of acute inflammation may be counterproductive for long-term outcomes of LBP sufferers.

Derivation of clinical prediction rules for identifying patients with non-acute low back pain who respond best to a lumbar stabilization exercise program at post-treatment and six-month follow-up

Low back pain (LBP) remains one of the most common and incapacitating health conditions worldwide. Clinical guidelines recommend exercise programs after the acute phase, but clinical effects are modest when assessed at a population level. Research needs to determine who is likely to benefit from specific exercise interventions, based on clinical presentation. This study aimed to derive clinical prediction rules (CPRs) for treatment success, using a lumbar stabilization exercise program (LSEP), at the end of treatment and at six-month follow-up. The eight-week LSEP, including clinical sessions and home exercises, was completed by 110 participants with non-acute LBP, with 100 retained at the six-month follow-up. Physical (lumbar segmental instability, motor control impairments, posture and range of motion, trunk muscle endurance and physical performance tests) and psychological (related to fear-avoidance and home-exercise adherence) measures were collected at a baseline clinical exam. Multivariate logistic regression models were used to predict clinical success, as defined by ≥50% decrease in the Oswestry Disability Index. CPRs were derived for success at program completion (T8) and six-month follow-up (T34), negotiating between predictive ability and clinical usability. The chosen CPRs contained four (T8) and three (T34) clinical tests, all theoretically related to spinal instability, making these CPRs specific to the treatment provided (LSEP). The chosen CPRs provided a positive likelihood ratio of 17.9 (T8) and 8.2 (T34), when two or more tests were positive. When applying these CPRs, the probability of treatment success rose from 49% to 96% at T8 and from 53% to 92% at T34. These results support the further development of these CPRs by proceeding to the validation stage.

Investigating the Causal Relationship Between Physical Activity and Chronic Back Pain: A Bidirectional Two-Sample Mendelian Randomization Study

Background: Recent observational studies have reported a negative association between physical activity and chronic back pain (CBP), but the causality of the association remains unknown. We introduce bidirectional Mendelian randomization (MR) to assess potential causal inference between physical activity and CBP. Materials and Methods: This two-sample MR used independent genetic variants associated with physical activity and CBP as genetic instruments from large genome-wide association studies (GWASs). The effects of both directions (physical activity to CBP and CBP to physical activity) were examined. Inverse variance-weighted meta-analysis and alternate methods (weighted median and MR-Egger) were used to combine the MR estimates of the genetic instruments. Multiple sensitivity analyses were conducted to examine the robustness of the results. Results: The MR set parallel GWAS cohorts, among which, those involved in the primary analysis were comprised of 337,234 participants for physical activity and 158,025 participants (29,531 cases) for CBP. No evidence of a causal relationship was found in the direction of physical activity to CBP [odds ratio (OR), 0.98; 95% CI, 0.85-1.13; p = 0.81]. In contrast, a negative causal relationship in the direction of CBP to physical activity was detected (β = -0.07; 95% CI, -0.12 to -0.01; p = 0.02), implying a reduction in moderate-vigorous physical activity (approximately 146 MET-minutes/week) for participants with CBP relative to controls. Conclusion: The negative relationship between physical activity and CBP is probably derived from the reduced physical activity of patients experiencing CBP rather than the protective effect of physical activity on CBP.