THE EFFECTS OF 12 WEEKS OF CHIROPRACTIC SPINAL ADJUSTMENTS ON PHYSIOLOGICAL BIOMARKERS IN ADULTS
 
   

The Effects of 12 Weeks of Chiropractic
Spinal Adjustments on Physiological Biomarkers
in Adults: A Pragmatic Randomized Controlled Trial

This section is compiled by Frank M. Painter, D.C.
Send all comments or additions to:
   Frankp@chiro.org
 
   

FROM:   PLoS One 2025 (Dec 11); 20 (12): e0338730 ~ FULL TEXT

  OPEN ACCESS   


Imran Amjad • Imran Khan Niazi • Nitika Kumari • Usman Ghani • Usman Rashid
Felipe Coutinho Kullmann Duarte • Federico Fortuna • Silvia Iglesias
Diego Gonzalez • Alex Sumich • Bibiana Fabre • Kelly Holt • Heidi Haavik

Centre for Chiropractic Research,
New Zealand College of Chiropractic,
Auckland, New Zealand.


Background:   Longer-term effects of chiropractic care on neuroplasticity, stress, and immune biomarkers remain unclear.

Objective:   This study evaluates the effects of chiropractic care on physiological biomarkers, including

brain-derived neurotrophic factor (BDNF),
cortisol (saliva, blood, hair), and
inflammatory cytokines
    interleukin-6 (IL-6),
    tumor necrosis factor-alpha (TNF-α),
    interferon-gamma (IFN-γ),
    C-reactive protein (CRP),
    B-lymphocytes (CD19),
    T-helper cells (CD4),
    cytotoxic T cells (CD8), and
    natural killer cells (CD56)

in subclinical spinal pain patients.

Methods:   Parallel-group, pragmatic randomized controlled trial conducted at the Rehabilitation Center of Railway General Hospital, Rawalpindi, Pakistan. Intervention: 12 weeks; follow-up: 16 weeks (May-December 2022). Participants with subclinical spinal pain were randomly assigned by using simple lottery method to either 12 weeks of chiropractic or sham care. We aimed to recruit up to 150 participants over three months; however, given the pragmatic nature of the trial and logistical constraints, including the availability of chiropractors, the final number enrolled was determined by how many eligible participants could be recruited during this time. Adults aged 20-60 years with subclinical spinal pain (n = 106 randomized; 88 completed 12-week measures; 73 completed 16-week follow-up). Among those who finished 12 weeks: chiropractic, 26 males/15 females, mean age 37.49 ± 12.39 years; sham, 24 males/23 females, mean age 26.85 ± 7.13 years.

The primary outcome blood BDNF and secondary outcome, including saliva, blood and hair cortisol, IL-6, TNF-α, IFN-γ, CRP, CD19, CD4, CD8, and CD56 levels were measured at baseline, after 12 weeks of intervention, and at a 16-week follow-up. Linear and linear mixed-effects regression models were used to assess the effect of care and time on biological measures.

Results:   Significant between-group differences were observed after 12 weeks of intervention, with higher salivary cortisol 5 ± 2 [0, 10], p = 0.045 and blood BDNF150 ± 60 (40, 270), p = 0.009 and IL-6 1.0 ± 0.3 [0.5, 1.5], p < 0.001 levels in the chiropractic care group.

At the 16-week follow-up, blood cortisol -9 ± 4 [-17, -1], p = 0.024, IFN-γ - 22 ± 7 [-35, -9], and TNF-α -2 ± 1 [-5, 0], p = 0.028 levels increased in the sham group.

Within-group comparisons showed a non-significant 10 ± 20 [-20, 50], p = 0.439 reduction in hair cortisol levels in the chiropractic group at 12 weeks, along with increased levels of blood cortisol, BDNF, CD8, CD4, IL-6, and CD19.

Conclusion:   12 weeks of Chiropractic care modulates biomarkers linked to neuroplasticity, inflammation, and stress. Increases in brain-derived neurotrophic factor and interleukin-6 suggest enhanced neuroplasticity and inflammatory responses, while decreases in tumor necrosis factor-alpha indicate a regulatory effect on systemic inflammation. These findings support the notion that chiropractic care modulates physiological systemic biomarkers, which may underscore its benefits on clinical outcomes.

Trial registration:   ClinicalTrials.gov NCT05369156



From the FULL TEXT Article:

Introduction

Chiropractors specialise in musculoskeletal health, primarily emphasizing the function and disorders of the spine, including spinal pain. During chiropractic care, chiropractors employ various conservative interventions with emphasis on spinal manipulative therapy (SMT), typically administered over multiple sessions. [1, 2] Based on previous research, it has been suggested that vertebral subluxations impact central neural function due to a breakdown in central segmental motor control (CSMC problems) that results in altered afferent feedback from muscle spindles in the paraspinal muscles at subluxated levels of the spine. [1, 3]

This, in turn, is thought to result in maladaptive neural plastic changes in the central nervous system that result in altered sensorimotor control. The chiropractic adjustments are shown to alter brain activity, including sensorimotor, prefrontal and cerebellar function, and improve the ability of peripheral muscles to produce force. [1, 4–14] Whilst maladaptive central activity might thus be improved following chiropractic care (CC) [1, 3], longer-term effects, and indeed their precise mechanisms, remain unclear.

Although chiropractic care is widely used and its clinical benefits on pain and function are recognized, little is known about its sustained effects on physiological biomarkers. understanding these effects is important to establish the clinical relevance and mechanistic basis of chiropractic interventions.

For example, one feasible mechanism could be through the immune and stress systems; that is, intervention could impact protective mechanisms (e.g., brain-derived Neurotrophic Factor, BDNF) and inflammation, which would be expected to have enduring effects on brain function.

There is also likely a bidirectional relationship between the function of key brain areas and inflammation, which could be manipulated through CC. For example, the prefrontal cortex (PFC) is centrally involved in multimodal sensory integration. Thus, alteration in PFC activity following chiropractic care is associated with

improved joint position sense error [15],
cortical processing [9, 16],
reflex excitability [17],
reaction time [16],
cortical sensorimotor integration [9, 10],
motor control [10] and
lower limb muscle strength. [18]

Through its involvement in regulating hypothalamic function and autonomic activity, the PFC also modulates endocrine and immune systems. [19–30] On the other hand, PFC is vulnerable to the damaging effects of stress and chronically raised inflammation. [31] Brain-derived neurotrophic factor (BDNF) protects against hypothalamic–pituitary–adrenal (HPA) dysregulation [32], as well as stress- and inflammation-induced neuronal damage. [33] A recent review [34] suggests chiropractic adjustments modulate immune mediators [35, 36], at least within the same day [36–40]; although again, longer-term effects are unclear. These findings indicate that chiropractic care can impact neuroplasticity, stress, and immune function through multiple pathways. however, evidence to date is mostly derived from acute or single-session studies in healthy adults, which limits our understanding of whether these changes are sustained or clinically meaningful in patient populations.

Immune status can be assessed by measuring the concentration of several inflammation mediators, such as tumor necrosis factor-alpha (TNF-α), C-reactive protein (CRP), interleukin (IL)-1 and IL-6, which are essential cytokines produced by monocytes and macrophages. Also, IFN-gamma activates the cells involved in the elimination of bacteria, viruses, fungi, and tumor cells – e.g., monocytes, macrophages, cytotoxic T-lymphocytes (CD3), natural killer cells (CD56) – and promotes the generation of antigen-specific IgG2 by B-lymphocytes (CD19). The ratio of T-helper cells (CD4) and Cytotoxic T cells (CD8) is an indicator of long-term activation of the immune system. [41] Moreover, increased levels of glucocorticoids (e.g., cortisol) can alter cytokine levels (such as increasing IL-6 pro-inflammatory cytokines), subsequently altering inflammation levels throughout the body. [42]

Cortisol can be measured from hair, blood, and saliva samples. Blood and salivary cortisol provide a snapshot of acute stress, reflecting current circulating levels. In contrast, hair cortisol is considered a marker of chronic stress, as it provides information on cortisol levels over weeks to months. [43] Recent studies [44–48] have highlighted the role of BDNF, cortisol, and cytokines in chronic pain, stress regulation, and healthy aging. however, very few have examined these biomarkers in the context of chiropractic care, particularly using longer-term follow-up designs. this study was therefore undertaken to address this gap and provide contemporary evidence.

Mechanical loading of various body tissues has been shown to modulate the release of pro- and anti-inflammatory cytokines, which are crucial for inflammation and tissue repair processes. [49–53] High-load cyclic compressive massage, in a pre-clinical study, increased macrophage infiltration and elevated muscle inflammation in the healthy tibialis anterior muscle. [54] Macrophages also play a significant role in the release of inflammatory cytokines, contributing to inflammation, intercellular signalling, and tissue repair. [55] Lower loading magnitudes of compressive massage stimulate the M2 subpopulation of macrophages, which are involved in tissue repair and regeneration, as opposed to the M1 subpopulation, which is associated with the secretion of pro-inflammatory cytokines. [56]

In another recent study, higher compared to lower SMT force raised plasma pro-inflammatory and dual-role cytokines in healthy, young adults. [57] Taken together, these data suggest that pro- and anti-inflammatory markers are sensitive to manual therapy and that CC might well influence these markers. However, most prior studies have been limited to single sessions or acute outcomes, few have investigated immune cell subset changes (CD4, CD8, CD19, CD56), and almost none have linked biomarker changes to clinical outcomes in patients with spinal pain. these gaps underscore the need for longitudinal, well-controlled studies.

Demonstrating longer-term effects of chiropractic care (CC) on immune and neuroplasticity markers would represent an important advance in understanding the physiological mechanisms underlying its clinical effects. We hypothesized that 12 weeks of chiropractic care would modulate biomarkers related to stress (cortisol), neuroplasticity (BDNF), and inflammation (IL-6, TNF-α, IFN-γ, CRP, and immune cell profiles) compared to sham care in adults with subclinical spinal pain.

To test this hypothesis, we conducted a pragmatic randomized controlled trial (RCT) examining the effects of CC on a range of blood, saliva, and hair-derived physiological biomarkers over a 16-week period. This study is clinically significant because it offers new mechanistic insights for researchers, provides clinicians with biomarker-informed evidence to guide practice, and demonstrates to patients that chiropractic care may improve neurophysiological resilience, immune balance, and stress regulation beyond symptom relief.



Materials and methods

      Design and setting

A parallel-group, pragmatic RCT was conducted at the Rehabilitation Center of Railway General Hospital, Rawalpindi, Pakistan (from May to December 2022). The Ethical Review Committee of Riphah International University, Pakistan, approved the study (Riphah/RCRS/REC/01288), which was registered with the National Institutes of Health ClinicalTrials.gov clinical trial registry (NCT05369156; first registration on 05 May 2022). The data collection period was from 15 May to 1 November 2022. All participants provided written informed consent before enrollment in the study.

      Study participants

Participants were included if they were aged between 20 and 60 years and had subclinical spinal pain. Subclinical spinal pain refers to recurrent spinal aches, pain, or stiffness, for which the person has not yet sought treatment. [58] The participants were pain-free at the time of their baseline assessment, i.e., less than or equal to 3/10 on a pain-based visual analogue scale (VAS). People with subclinical pain were recruited because subclinical pain participants are more likely to have altered neurophysiological and biomechanical function than completely pain-free participants. [58]

Participants were excluded from the study if any of the following conditions applied:

  • no evidence of spinal dysfunction is present,

  • they are in current pain (above 3/10 on VAS),

  • they had sought previous treatment for their spinal issues,

  • they were unable to perform the assessment procedures due to contra-indications or movement limitations,

  • were diagnosed with immune dysfunction and/or utilized a prescribed immunosuppressive medication.

Participants were also excluded if they had uncontrolled asthma, had nasal polyps, used of an intranasal steroid spray one month or less before the study, were HIV-positive, had abnormal laboratory values (i.e., any change from reference values in the last 12 months) for complete blood count or hormonal changes, were pregnant or breastfeeding women, had a history of drug abuse, was unable or unwilling to comply with the study protocol (i.e., relating to blood draws), were participating in another research study during the time of data collection, had any diagnosed comorbidity or concomitant disease (diabetes, hypertension or any neurological condition), had donated blood within last month, had allergies to yeast or yeast-derived products, or had chronic sinusitis and/ or recent (within the last six weeks) episode of acute sinusitis. Written consent was obtained from all volunteers before they participated in the study.

      Sample size

The sample size was calculated based on 80% power and a 95% significance level. A total of 150 participants was sufficient to detect an effect size of 0.25 in the tests of multimodal data included as primary outcome measures in this study. [59] We planned to allocate 75 participants to the chiropractic care (CC) group and 75 to the control group. Recruitment was open for a maximum of three months, during which we aimed to enrol as many eligible participants as possible, to reach 150 participants. Recruitment was limited to a fixed three-month period due to the availability of chiropractors to deliver the intervention.

      Interventions

The interventions were 12 weeks of chiropractic and 12 weeks of sham chiropractic care. Participants were allowed to receive usual care if required as part of the best ethical practice to ensure symptom relief, avoid withholding standard care, and improve the acceptability of participation.

      Chiropractic intervention.

In the chiropractic group, a New Zealand registered chiropractor checked participants for central segmental motor control (CSMC) problems and delivered chiropractic adjustments to the appropriate segments during the experimental sessions.. Participants were checked by the chiropractor approximately three times per week for 12 weeks, with interventions delivered pragmatically to reflect routine clinical practice, allowing adjustments to be tailored based on individual clinical assessments. At each experimental session the chiropractor assessed the participant’s spine for indicators of CSMC problems. These indicators are routinely used by chiropractors to determine the appropriate segments of the spine and pelvis that require chiropractic adjustments and, when used in combination, have been shown to be reliable indicators for determining CSMC problems. [60, 61]

The indicators used include static and motion palpation to assess for spinal tenderness, restricted intersegmental motion, increased muscle tension, and reduced joint play. Once the site of the CSMC problem has been identified the chiropractor delivers the adjustments either with a manual, high-velocity, low amplitude thrust (HVLA) or with an instrument-assisted thrust. [62] Typically, multiple sites of the spine were adjusted at each experimental session, as deemed appropriate by the chiropractor based on the spinal assessment. Each chiropractic visit lasted approximately 15–20 minutes.

      Sham (Control) chiropractic intervention.

Research involving physical or manual interventions can face challenges with the blinding of participants. [63, 64] Whilst it is difficult to fully blind participants in a study of this nature, it is possible to control contextual factors to enhance the chances of successful blinding. This study was conducted in Pakistan where chiropractic care is relatively unknown. [36] In a recent survey of university students in Lahore, Pakistan, including pharmacy students, more than two-thirds of respondents were unaware that chiropractic care involved spinal manipulation and that it is used as an intervention for low back pain. [65]

The limited understanding of chiropractic care presents a unique opportunity to investigate its effects while maximizing the chances of successful participant blinding. By carefully controlling contextual factors in this study design, we can more effectively isolate and analyze the specific contributions of chiropractic treatment beyond those of placebo and contextual effects. The control group, in this case, received a sham chiropractic intervention.

This involves participants attending experimental sessions at the same frequency and duration as the intervention group. At these sessions the same chiropractor performed a spinal assessment but no chiropractic adjustments are performed. Instead, the chiropractor would position the patient as they normally would to deliver the adjustment, but with no thrust on the spine. Alternatively, the chiropractor would use the adjusting instrument on a low setting to produce the clicking sound, but with no contact on the participants spine.

      Outcome measures

All outcome measures were assessed at baseline (before the first chiropractic session), 12 weeks (post-intervention), and at 16 weeks (to assess retention effects). Although multiple measures were collected as part of the broader study, this paper specifically focuses on blood-based, saliva and Hair markers.

The primary outcome measure was the blood BDNF levels. Secondary outcome measures included the B-Lymphocytes (CD19), CD4, CD8, IL-6, TNF-alpha, C-reactive protein (CRP), interferon-gamma, Natural killer cells (CD56) and blood, hair, and saliva cortisol.

Following the 16-week (follow-up session), participants in both groups were asked to indicate whether they thought they had received active chiropractic care using a yes or no response to test the effectiveness of participant blinding.

Blood samples were collected from participants via venipuncture into sterile blood separation anticoagulant tubes. A total of 5 mL of blood was drawn from each participant. After collection, the samples were allowed to clot at room temperature for 30 minutes and were subsequently centrifuged at 3000 x g for 10 minutes to separate the blood. The blood was divided into two microcentrifuge tubes and immediately stored at –80°C until analysis.

Each participant was provided with a labelled plastic container for saliva collection and instructed to collect their saliva early in the morning, before 7 o’clock. These containers were then stored at –80°C. Each sample underwent one freeze-thaw cycle before analysis to minimize the degradation of biomarkers. Hair samples were collected by selecting a small section of hair from the vertex area, cutting it with scissors, and placing it in aluminium foil. The samples were then stored at room temperature.

Blood concentrations of biomarkers (explained earlier) and cortisol concentration from saliva were evaluated using commercially available enzyme-linked immunosorbent assay (ELISA) kits (abcam, Elabscience, ELK Biotechnology) following the manufacturer’s instructions. The optical density (OD) was measured using a microplate reader (RaytoRC 2100 C, China) at a wavelength of 450 nm. The samples were analyzed at Health Plus Lab, Rawalpindi, Pakistan. Calibration curves were generated for each assay, and sample concentrations were calculated using a four-parameter logistic (4PL) curve-fitting model provided by the ELISA kit software.

Assay reliability was confirmed by intra- and inter-assay coefficients of variation consistently <10–15% across all kits. Each plate included standardised controls and calibration curves for quantification. Validity was further supported by manufacturer data, which confirmed high specificity with no significant cross-reactivity (<5%) and reproducibility across various biological matrices.

Analytical sensitivity (LOD/LLOQ) varied by biomarker, ranging from

0.19 ng/mL for CD4,
0.23 ng/mL for CRP, and
0.118 ng/mL for CD19 to
18.75 pg/mL for BDNF.

Detection ranges spanned several orders of magnitude, for example:

Cortisol 6.25–400 ng/Ml;
IFN-γ 15.63–1000 pg/mL;
TNF-α 7.81–500 pg/mL;
CD8 1.25–80 ng/mL;
IL-6 1.56–100 pg/mL
CD56 1.56–100 ng/mL;
CD19 0.32–20 ng/mL

Recovery rates for serum and plasma samples were typically within ~90–110%, in line with immunoassay validation guidelines, supporting the accuracy of quantification.

Hair samples were collected using scissors from the posterior vertex, as close to the scalp as possible, to evaluate the hair cortisol concentration. Each sample was stored in a paper envelope at room temperature. Cortisol was extracted and quantified using an automated chemiluminescent method (Immulite 2000 autoanalyzer, Siemens, LA, USA) following the protocol reported by Gonzalez et al. Hair cortisol concentration was expressed in pg/mg of hair. [66]

Potential harms or adverse events were investigated by asking participants, during scheduled intervention visits, about any injuries or perceived adverse effects of care that may have occurred during the trial, such as muscle soreness or pain. [67]

      Randomization and blinding

Randomization was carried out following the baseline assessment by a simple lottery method. All participants and the outcome assessors were blinded to group allocation. An investigator undertook and collated all outcome measurements and remained blinded to participant allocation throughout the study, the outcome assessor was not involved in administering the interventions.

The statistician who analyzed the data was also blinded to group allocation, as all recorded data were anonymized and coded before being provided for analysis. The chiropractors providing chiropractic care could not be blinded to group allocation.

      Statistical analysis

The primary null hypothesis for the analysis was that there was no difference in blood BDNF levels between the chiropractic and sham groups. Pre-specified secondary null hypotheses stated that there were no differences between the two groups in CD19, CD4, CD8, IL-6, TNF-alpha, CRP, interferon-gamma, CD56 and blood and saliva cortisol.

Linear and linear mixed effects regression models were constructed to analyze the pre- to post-intervention change scores of the outcome measures. Multiple outcome measures were included in a single model if they were correlated and had similar error variance. This multivariate analysis of covariance using linear mixed effects modelling is statistically superior to separate models for each outcome as it results in better parameter estimation. [68] The models regressed the change scores on the pre-intervention values, group, and time in the case of outcomes with follow-up, outcome, and interactions between these fixed factors. The within-participant correlations arising from repeated measures were accounted for by including participant-level random intercepts in the models.

Where a model included multiple outcomes, these random intercepts were estimated for each outcome. The linear regression model was only fitted to the hair cortisol as it was only measured at pre- and post-intervention time points. Linear mixed effects regression was used for all the other outcomes. The assumptions of normality and homogeneity of error variance were evaluated with QQ-plots and errors-versus-fitted values plots. A uniform dispersion of the residuals around the zero line for all fitted values was the expected norm.

Any deviations from this norm, such as a fan-out or fan-in pattern, were taken as an indication of poor model fit to the data. In such cases, a separate model was fitted for each outcome, rather than pooling multiple outcomes in a single model.

The model parameters are used to estimate between and within-group differences at successive time points. These mean differences in change scores were presented with their standard errors, 95% confidence intervals and hypothesis test statistics (t-values, degrees of freedom, p-values) for 0 change at 5% type-I error rate.

To test these hypotheses, data collected during the study were collated in an Excel (Microsoft Corp., Redmond, WA) spreadsheet with groups relabeled for blinding. Statistical analyses were conducted in R version 4.0.0 using the packages lme4. [69–71] A detailed report of the statistical analysis is available in the supplementary file (S1 File).

Mean differences for between-group and within-group differences are reported with standard errors and 95% confidence intervals that were obtained from the models. P-values for testing the primary and secondary null hypotheses were obtained with z or t-tests based on these estimates. Cohen’s d corresponding to the mean differences is also reported. The statistical significance level was set at 0.05. Šidák’s correction was applied to multiple comparisons from the same model to control the family-wise error rate for correlated outcomes and timepoints.

Due to the difference in the number of participants willing to give samples for different tests after group allocation, a per-protocol analysis was performed. Nonetheless, to evaluate the sensitivity of the statistically significant results to missing data, the primary models were refit with multiple imputation for missing values. [72] Statistically significant results that remained statistically significant in the sensitivity analysis were marked with an obelisk ().



Results

Figure 1

Table 1

Figure 2

Table 2

Table 3
page 12

One hundred forty-four volunteers were recruited and checked for eligibility during a three-month recruitment period. 106 participants met the inclusion criteria and agreed to participate. 88 completed the post-12-week measure(n = 41 in chiro group,n = 47 in sham group), and 73 completed the 16-week (n = 32 in chiro group,n = 41 in sham group) follow-up session (Figure 1, CONSORT diagram). A few participants (n = 5) reported mild transient muscle soreness and stiffness in the back muscles following the first intervention visit in the interventional group, resolving spontaneously within 24 hours without intervention; otherwise, no adverse events or reports of harm were received during the trial. Participants’ clinical characteristics are given in Table 1, who finished 12 weeks post-intervention.

      Participant blinding

After the 16-week intervention period, all participants who completed the trial were asked if they thought they had received active chiropractic care. 39 of the 41 participants in the chiro group and 43 of the 47 participants in the sham group believed they had received active chiropractic care. This suggests that participant blinding was successful, with 94% of participants across the two groups believing they had received active chiropractic care, with no between-group differences present.

Figure 2 shows line graphs representing group differences at 12 and 16 weeks. Table 2 shows t-statistics for between-group comparisons and 12 and 16 weeks. Table 3 shows t-statistics for within-subjects effects.

      Blood BDNF concentration

There was a significant effect of Group for BDNF levels (Chiro>Sham) at 12 weeks (mean difference = 150 ± 60; 95% CI: 40–270; p = 0.009), but this effect failed to reach significance at 16 weeks (120 ± 60; 95% CI: 0–250; p = 0.059).

Within-group comparison shows a significant increase in blood BDNF levels in the chiropractic care group after 12 weeks (80 ± 40; 95% CI: 0–170; p = 0.050), compared to the control group. At the 16-week follow-up, there is no significant difference between the chiropractic care group and the control group.

      Blood, hair and saliva cortisol

No significant between-group difference for serum cortisol at 12 weeks (–2 ± 4; 95% CI: –9–6; p = 0.627). At 16 weeks, the chiropractic group showed significantly lower serum cortisol (–9 ± 4; 95% CI: –17 to –1; p = 0.024).

For hair cortisol no significant difference at 12 weeks (10 ± 20; 95% CI: –20–50; p = 0.439). A within-group trend toward reduction was observed in the chiropractic group at 12 weeks (20 ± 10; 95% CI: 0–50; p = 0.070). At 12 weeks, salivary cortisol was higher in the chiropractic group (5 ± 2; 95% CI: 0–10; p = 0.045), but no significant effect was found at 16 weeks (1 ± 3; 95% CI: –5–6; p = 0.775).

      Blood CD19, CD4, CD56, and CD8

There were no significant effects of Group for CD19, CD4, CD56, or CD8 at 12 weeks or 16 weeks. However, there was a significant increase in CD19, CD4, and CD8 in the chiropractic care group at 12 weeks, which was not seen in the sham group. But CD56 showed a significant increase in the chiropractic care group. There was no significant effect of Group or Time for CD19, CD4, CD56, or CD8 at 16 weeks.

      Blood CRP, IL-6, INF-γ and TNF-α concentration

No significant effects of Group in CRP were seen at either 12 (1 ± 1; 95% CI: –1–3; p = 0.371) or 16 weeks (0 ± 1; 95% CI: –3–2; p = 0.675). However, there was a significant increase between baseline and 12 weeks, with no further change between 12 and 16 weeks. The chiropractic care group showed significantly higher IL-6 than Sham at 12 weeks(1.0 ± 0.3; 95% CI: 0.5 to 1.5; p < 0.001), but not at 16 weeks (0.4 ± 0.3; 95% CI: –0.2 to 0.9; p = 0.170).

No effect on blood INF-γ levels at 12 weeks (4 ± 6; 95% CI: –8–17; p = 0.465), but significantly lower in the chiropractic group at 16 weeks (–22 ± 7; 95% CI: –35 to –9; p = 0.001). This appeared to be due to a significant increase in the Sham group between 12 and 16 weeks. There was a significant lower TNF-alpha at 12 (–2 ± 1; 95% CI: –4–0; p = 0.023) and 16 weeks (–2 ± 1; 95% CI: –5–0; p = 0.028) in the chiropractic group.



Discussion

This is the first investigation of the longer-term effects of Chiropractic care on biological markers of neuroplasticity, stress-response and immune system function. In line with hypotheses, present findings demonstrate Chiropractic care affects BDNF levels (higher in the chiropractic care group), blood TNF-α levels (lower in the chiropractic care group) and saliva cortisol levels (higher in the chiropractic care group) after 12 weeks. Additionally, at the 16-week follow-up, there were significant differences in blood cortisol levels (higher in the sham group), blood INF-γ levels (higher in the Sham group), and blood TNF-α levels (higher in the Sham group).

Within-group comparisons showed a significant reduction in hair cortisol levels in the chiropractic care group after 12 weeks, while blood cortisol, saliva cortisol, blood BDNF, CD8, CD4, INF-γ, TNF-α, IL-6, and CD19 levels increased significantly in the chiropractic care group after 12 weeks. There were no significant differences between groups or within groups for blood CD8, CD4, CD56, or CRP levels at both time points.

The relatively small effect sizes and wide confidence intervals for some biomarkers (e.g., cortisol, IFN-γ) suggest that while trends were observed, the intervention’s effectiveness may be modest, and results should be interpreted with caution.

This variability highlights the need for larger trials with sufficient power. Although we anticipated that changes in immune and stress biomarkers may be partly mediated by altered central neural function, direct measures of brain activity were not reported in this study. Our previous studies have shown changes in cortical excitability and prefrontal activation following CC [9, 13], providing initial evidence that neuroimmune modulation may cause such effects. In future, we plan to integrate neurophysiological and biomarker measures to formally test this mediation hypothesis using structural equation or multilevel mediation models.

The increase in BDNF levels observed in the chiropractic group suggests enhanced neuroplasticity compared to the control group, indicating that chiropractic intervention may stimulate neuroplastic activity in the central nervous system (CNS). [73] BDNF, part of the neurotrophin family along with nerve growth factor (NGF), is crucial for the initiation and maintenance of CNS adaptations. [74] For instance, previous evidence studying the effects of spinal manipulation on growth factors observed that spinal manipulation, commonly used in chiropractic practice, stretches spinal muscles and stimulates the mechano-growth factor (MGF), a splice product of the insulin-like growth factor (IGF-1) gene. MGF is involved in muscle and neuronal growth and repair. [75] Furthermore, non-noxious mechanical skin stimulation has been reported to stimulate NGF release from cortical neurons, which plays a key role in neuronal survival and function. [76] These findings underscore the potential for chiropractic interventions to influence neurotrophic factors and promote neuroplasticity, thereby enhancing neuronal health and function. [77]

The chiropractic intervention also resulted in higher serum levels of IL-6 after 12 weeks. Other inflammatory cytokines, such as TNF-α, were significantly higher in the sham group compared to the chiropractic group at both 12 and 16 weeks, and IFN-γ was higher in the sham group at 16 weeks. Although IL-6 is well known to be released in large amounts during infection, autoimmunity and cancer, it is also released by skeletal muscles during exercise. [78] When muscles are used intensely, they release IL-6, in which case it is considered a myokine, where it plays a role in local regulation of the inflammatory process and skeletal muscle regeneration. [78] Thus, the changes in inflammatory markers found in the current study may be explained in the following manner. It has been documented that deep paraspinal muscles undergo numerous maladaptive plastic processes over time following an injury, including muscle atrophy, fatty infiltration, changes in fiber type and morphological changes of the muscle spindles. [79–83]

Thus, the chiropractic intervention in this study, which consisted of repeated HVLA adjustments of such dysfunctional spinal segments, may be considered akin to a kind of exercise for these deep paraspinal muscles. The HVLA intervention used as part of the chiropractic intervention rapidly stretches these deep muscles and is known to produce ongoing increased intervertebral range of motion at that segment for a time after the adjustment. [84] Thus, the increase in IL-6 found in this study may reflect an increased release of this myokine in response to the repeated HVLA adjustment provided over the 12 weeks.

As the adjustments stopped at 12 weeks, the release of IL-6 would naturally be reduced by 16 weeks; hence, we found no significant difference in IL-6 levels at 16 weeks. IFN-γ and TNF-α are produced by immune cells (like macrophages and T cells) and non-immune cells (such as fibroblasts and endothelial cells). [57] They are known to play pro-inflammatory and signalling roles, while others, like IL-4 and IL-10, have anti-inflammatory effects. Interestingly, these pro-inflammatory cytokines were not increased in the chiropractic group at either the 12-week mark or the 16-week mark. However, they were significantly higher in the sham group, suggesting the sham group, which did not receive CC, was steadily getting worse (more inflamed) over the study period. This finding supports the notion that the IL-6 changes seen in the chiropractic group at 12 weeks most likely reflect skeletal muscle regeneration.

The cortisol changes observed in this study were interesting. In the current study, the saliva cortisol levels were significantly increased in the chiropractic group at 12 weeks, yet the serum cortisol levels were significantly decreased at 16 weeks. Cortisol generally rises following acute injury to modulate and limit the initial pro-inflammatory response. [85] In the context of CC, the increase in saliva cortisol at 12 weeks could be a response to the mechanical stress and strain associated with the HVLA chiropractic adjustments. [57, 75–77] Interestingly, the overall serum cortisol levels at 16 weeks ended up lower than at baseline, suggesting that overall inflammation levels were significantly reduced 4 weeks after the 12-week intervention had stopped. This finding is in line with the IFN-γ and TNF-α immune markers as well.

IFN-γ plays a role in activating the immune system in response to damaged cells. As there was no increase in IFN-γ in the chiropractic group, while this significantly increased in the sham group, this suggests that chiropractic care does not damage cells and is more likely to be linked to exercise for the deep paraspinal muscles surrounding the vertebral subluxations.

The elevated cortisol levels observed after the chiropractic intervention may be due to an influence of the neuroimmunomodulation pathway and the influence of the HPA axis. [86, 87] Increased cortisol levels could suppress systemic inflammatory agents like TNF-α, helping fine-tune immune and inflammatory activity during tissue repair. [86] This complex interplay highlights the body’s regulatory mechanisms in balancing inflammation and immune response following chiropractic care.

After 12 weeks of intervention, there were no significant changes in CRP. When IL-6 is released as a pro-inflammatory marker, it regulates the acute release of CRP from the liver. Yet the current study did not find an increase in CRP following the chiropractic intervention. Thus, following 12 weeks of chiropractic care, it is more likely that the increase seen in the IL-6 may be due to its release as myokines in the deep paraspinal muscles following the repeated HVLA adjustments.

This notion is further supported by the finding of TNF-α levels being significantly less in the chiropractic group at both the post-12-week intervention period and at the post-16-week time point. This reduction in TNF-α levels is also consistent with previous studies showing that chiropractic care, including SMT, can lead to decreased TNF-α levels. [88–91] TNF-α plays a crucial role in initiating and propagating acute inflammation. Its reduction could result from a negative feedback mechanism in response to the initial acute inflammatory reaction, aiming to modulate overall systemic inflammation, which is largely mediated by TNF-α. [92]

      Strengths and limitations

Statistically significant between-group differences were observed in BDNF, IL-6, TNF, and salivary cortisol levels at the 12-week mark with a sample size of 88 participants. These results suggest that the study was adequately powered for these specific outcome measures at the primary endpoint. The study was designed pragmatically to recruit up to 150 participants; however, the final enrollment (n = 106 randomized) reflected real-world feasibility constraints, including the availability of chiropractors within a fixed three-month recruitment window.

However, we acknowledge that this sample size may have been insufficient to detect between-group differences in other outcome measures. Consequently, the possibility of type II errors cannot be ruled out. As an exploratory study, we employed a wide range of outcome measures and performed multiple comparisons without applying adjustments to p-values. While this approach increases the likelihood of type I errors, it is considered appropriate when investigating novel areas of research, such as the present study. [93]

Nevertheless, this remains a limitation that should be addressed in future studies. Moreover, the relationship and implications of the expression of these physiological markers should be evaluated in relation to clinical outcomes. In the current study, we also have some data on clinical outcomes, and we plan to conduct a secondary analysis on this.

One strength of the current study is the use of a broad array of inflammatory markers, which provided a more comprehensive understanding of the effects of chiropractic care on physiological markers. This approach will allow for a more detailed and specific evaluation of the mechanisms behind the effectiveness of chiropractic interventions. Future large-scale RCTs should target more specific mediators of inflammation and immune responses influenced by cortical activity. It may be beneficial to explore specific biochemical cascades or pathways, employing more refined and sophisticated outcome measurement techniques, and potentially even examining changes at the level of messenger RNA (mRNA) expression. The sample size estimates from this study can serve as a guide for powering these future trials. Moreover, future research should incorporate longer intervention periods and extended follow-up durations, potentially up to six months, with a more pragmatic approach to participant care.

Another strength of the current study is the success in maintaining participant blinding. Over 94% of participants, regardless of their group allocation, believed they had undergone an active chiropractic intervention, demonstrating that blinding was adequately achieved. This is particularly challenging in trials involving manual interventions [64, 94], yet it reinforces the conclusion that observed between-group differences were not merely a result of contextual or placebo effects.



Conclusion

In conclusion, this study demonstrates that 12 weeks of chiropractic care was associated with increased blood brain-derived neurotrophic factor (BDNF) and IL-6 levels and decreased TNF-α levels at 12 weeks, along with reduced blood cortisol and IFN-γ levels at 16 weeks compared with sham care. These findings suggest that chiropractic care may modulate neurotrophic, stress, and inflammatory pathways. While the effects were modest and some confidence intervals were wide, the results provide preliminary evidence for physiological changes beyond musculoskeletal outcomes. Larger and longer-term trials are needed to confirm these findings and to clarify their clinical relevance.


Supplementary Material

S1 File Detailed statistical analysis report.

CONSORT 2025 Checklist CONSORT 2025 checklist item description


Acknowledgments

The authors thank all the patients and the staff from the Rehabilitation Centre of Railway General Hospital, Rawalpindi, Pakistan, who contributed to this study. We would also like to acknowledge Dr Abdul Rahaman from Balmoral Chiropractic Clinic (Auckland, New Zealand) for his valuable clinical collaboration and support during this project.


Abbreviations

BDNF = Brain-Derived Neurotrophic Factor
CNS = Central Nervous System
CSMC = Central Segmental Motor Control
CC = Chiropractic Care
CRP = C-Reactive Protein
CD8 = Cytotoxic T cells
ELISA = Enzyme-Linked Immunosorbent Assay
IGF-1 = Insulin-Like Growth Factor
IL = Interleukin
CD19 = Lymphocytes
MGF = Mechano-Growth Factor
CD56 = Natural Killer Cells
NGF = Nerve Growth Factor
PFC = Prefrontal Cortex
RCT = Randomized Controlled Trial
SMT = Spinal Manipulative Therapy
CD4 = T-helper cells
TNF-α = Tumor Necrosis Factor-Alpha
VAS = Visual Analogue Scale.


Funding Statement

This study was funded by the Centre for Chiropractic Research Supporters Program at the New Zealand College of Chiropractic, the Australian Spinal Research Foundation (ASRF), Australia, the NCMIC Foundation, Iowa, USA, and the Hamblin Trust, New Zealand. Funding was awarded to Heidi Haavik (ASRF, NCMIC Foundation) and Imran Khan Niazi (Hamblin Trust, ASRF). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.



References:

  1. Haavik H, Murphy B.
    The Role of Spinal Manipulation in Addressing Disordered
    Sensorimotor Integration and Altered Motor Control

    J Electromyogr Kinesiol. 2012 (Oct); 22 (5): 768–776

  2. Henderson CN.
    The Basis for Spinal Manipulation:
    Chiropractic Perspective of Indications and Theory

    J Electromyography and Kinesiology 2012 (Oct); 22 (5): 632–642

  3. Haavik Taylor H, Holt K, Murphy B.
    Exploring the Neuromodulatory Effects of the
    Vertebral Subluxation and Chiropractic Care

    Chiropractic J Australia 2010 (Mar); 40 (1): 37–44

  4. Christiansen TL, Niazi IK, Holt K, Nedergaard RW, et al.
    The Effects of a Single Session of Spinal Manipulation
    on Strength and Cortical Drive in Athletes

    Eur J Appl Physiol 2018 (Apr); 118 (4): 737-749

  5. Daligadu J, Haavik H, Yielder PC, Baarbe J, Murphy B.
    Alterations in Cortical and Cerebellar Motor
    Processing in Subclinical Neck Pain Patients
    Following Spinal Manipulation

    J Manipulative Physiol Ther. 2013 (Oct); 36 (8): 527–537

  6. Haavik H, Niazi IK, Holt K, Murphy B.
    Effects of 12 Weeks of Chiropractic Care on Central
    Integration of Dual Somatosensory Input in
    Chronic Pain Patients: A Preliminary Study

    J Manipulative Physiol Ther. 2017 (Mar); 40 (3): 127–138

  7. Haavik H, Niazi IK, Jochumsen M, Sherwin D.
    Impact of Spinal Manipulation on Cortical Drive
    to Upper and Lower Limb Muscles

    Brain Sci. 2017 (Jan); 7 (1): 2

  8. Haavik H, Niazi IK, Jochumsen M, Ugincius P, et al.
    Chiropractic spinal manipulation alters TMS induced
    I-wave excitability and shortens the
    cortical silent period.
    J Electromyogr Kinesiol. 2018;42:24–35.
    doi: 10.1016/j.jelekin.2018.06.010

  9. Haavik-Taylor H, Murphy B.
    Cervical Spine Manipulation Alters Sensorimotor
    Integration: A Somatosensory Evoked Potential Study

    Clin Neurophysiol. 2007 (Feb); 118 (2): 391–402

  10. Taylor HH, Murphy B.
    Altered sensorimotor integration with
    cervical spine manipulation.
    J Manipulative Physiol Ther. 2008;31(2):115–26.
    doi: 10.1016/j.jmpt.2007.12.011

  11. Taylor HH, Murphy B.
    Altered Central Integration of Dual Somatosensory
    Input After Cervical Spine Manipulation

    J Manipulative Physiol Ther. 2010 (Mar); 33 (3): 178–188

  12. Holt K, Niazi IK, Nedergaard RW, Duehr J, Amjad I, et al.
    The effects of a single session of chiropractic care
    on strength, cortical drive, and spinal
    excitability in stroke patients.
    Sci Rep. 2019;9(1):2673.
    doi: 10.1038/s41598-019-39577-5

  13. Lelic D, Niazi IK, Holt K, Jochumsen M, Dremstrup K, Yielder P, et al.
    Manipulation of Dysfunctional Spinal Joints Affects
    Sensorimotor Integration in the Prefrontal Cortex:
    A Brain Source Localization Study

    Neural Plast. 2016 (Mar 7); 2016: 3704964

  14. Niazi IK, Türker KS, Flavel S, Kinget M, Duehr J, Haavik H.
    Changes in H-reflex and V-waves
    Following Spinal Manipulation

    Exp Brain Res. 2015 (Apr); 233 (4): 1165–1173

  15. Haavik H, Murphy B.
    Subclinical Neck Pain and the Effects of Cervical
    Manipulation on Elbow Joint Position Sense

    J Manipulative Physiol Ther. 2011 (Feb); 34 (2): 88–97

  16. Kelly DD, Murphy BA, Backhouse DP.
    Use of a mental rotation reaction-time paradigm to
    measure the effects of upper cervical adjustments
    on cortical processing: a pilot study.
    J Manipulative Physiol Ther. 2000;23(4):246–51.
    doi: 10.1067/mmt.2000.106099

  17. Herzog W, Scheele D, Conway PJ.
    Electromyographic Responses of Back and Limb Muscles
    Associated With Spinal Manipulative Therapy

    Spine (Phila Pa 1976) 1999 (Jan 15); 24 (2): 146–152

  18. Hillermann B, Gomes AN, Korporaal C, Jackson D.
    A pilot study comparing the effects of spinal
    manipulative therapy with those of extra-spinal
    manipulative therapy on quadriceps muscle strength.
    J Manipulative Physiol Ther. 2006;29(2):145–9.
    doi: 10.1016/j.jmpt.2005.12.003

  19. Critchley HD, Harrison NA.
    Visceral influences on brain and behavior.
    Neuron. 2013;77(4):624–38.
    doi: 10.1016/j.neuron.2013.02.008

  20. Machado S, Cunha M, Velasques B, Minc D, Teixeira S, et al.
    Sensorimotor integration: basic concepts, abnormalities
    related to movement disorders and sensorimotor
    training-induced cortical reorganization.
    Rev Neurol. 2010;51(7):427–36.

  21. Thayer JF.
    Vagal tone and the inflammatory reflex.
    Cleve Clin J Med. 2009;76 Suppl 2:S23-6.
    doi: 10.3949/ccjm.76.s2.05

  22. Ahern GL, Sollers JJ, Lane RD, Labiner DM, Herring AM, et al.
    Heart rate and heart rate variability changes in
    the intracarotid sodium amobarbital test.
    Epilepsia. 2001;42(7):912–21.
    doi: 10.1046/j.1528-1157.2001.042007912.x

  23. Berthoud HR, Neuhuber WL.
    Functional and chemical anatomy of the afferent vagal system.
    Autonomic Neuroscience. 2000;85(1–3):1–17.

  24. Thayer JF, Sternberg EM.
    Neural aspects of immunomodulation:
    focus on the vagus nerve.
    Brain Behav Immun. 2010;24(8):1223–8.
    doi: 10.1016/j.bbi.2010.07.247

  25. Bankenahally R, Krovvidi H.
    Autonomic nervous system: anatomy, physiology, and
    relevance in anaesthesia and critical care medicine.
    BJA Education. 2016;16(11):381–7.
    doi: 10.1093/bjaed/mkw011

  26. Kulchynskyi AB, Kyjenko VM, Zukow W, Popovych IL.
    Causal neuro-immune relationships at patients with
    chronic pyelonephritis and cholecystitis.
    correlations between parameters EEG,
    HRV and white blood cell count.
    Open Medicine. 2017;12(1):201–13.

  27. Ohira H, Matsunaga M, Osumi T, Fukuyama S, et al.
    Vagal nerve activity as a moderator of
    brain-immune relationships.
    J Neuroimmunol. 2013;260(1–2):28–36.
    doi: 10.1016/j.jneuroim.2013.04.011

  28. McCraty R, Shaffer F.
    Heart Rate Variability: New Perspectives on
    Physiological Mechanisms, Assessment of
    Self-regulatory Capacity, and Health risk.
    Glob Adv Health Med. 2015;4(1):46–61.
    doi: 10.7453/gahmj.2014.073

  29. Hänsel A, von Känel R.
    The ventro-medial prefrontal cortex: a major link
    between the autonomic nervous system, regulation
    of emotion, and stress reactivity?
    Biopsychosoc Med. 2008;2:21.
    doi: 10.1186/1751-0759-2-21

  30. Diorio D, Viau V, Meaney MJ.
    The role of the medial prefrontal cortex
    (cingulate gyrus) in the regulation of
    hypothalamic-pituitary-adrenal
    responses to stress.
    J Neurosci. 1993;13(9):3839–47.
    doi: 10.1523/JNEUROSCI.13-09-03839.1993

  31. Joyce MKP, Uchendu S, Arnsten AFT.
    Stress and Inflammation Target Dorsolateral Prefrontal
    Cortex Function: Neural Mechanisms Underlying
    Weakened Cognitive Control.
    Biol Psychiatry. 2025;97(4):359–71.
    doi: 10.1016/j.biopsych.2024.06.016

  32. Correia AS, Cardoso A, Vale N.
    BDNF Unveiled: Exploring Its Role in Major
    Depression Disorder Serotonergic Imbalance
    and Associated Stress Conditions.
    Pharmaceutics. 2023;15(8):2081.
    doi: 10.3390/pharmaceutics15082081

  33. Carballedo A, Morris D, Zill P, Fahey C, Reinhold E, et al.
    Brain-derived neurotrophic factor Val66Met
    polymorphism and early life adversity
    affect hippocampal volume.
    Am J Med Genet B Neuropsychiatr Genet. 2013;162B(2):183–90.
    doi: 10.1002/ajmg.b.32130

  34. Haavik H, Niazi IK, Kumari N, Amjad I, Duehr J, Holt K.
    The Potential Mechanisms of High-Velocity,
    Low-Amplitude, Controlled Vertebral Thrusts
    on Neuroimmune Function: A Narrative Review

    Medicina (Kaunas) 2021 (Jun); 57 (6): 536

  35. Kovanur-Sampath K, Mani R, Cotter J, Gisselman AS, Tumilty S.
    Changes in biochemical markers following spinal
    manipulation-a systematic review
    and meta-analysis.
    Musculoskelet Sci Pract. 2017;29:120–31.
    doi: 10.1016/j.msksp.2017.04.004

  36. Plaza-Manzano G, Molina-Ortega F, Lomas-Vega R.
    Changes in Biochemical Markers of Pain Perception
    and Stress Response After Spinal Manipulation

    J Orthop Sports Phys Ther. 2014 (Apr); 44 (4): 231–239

  37. Teodorczyk-Injeyan JA, Injeyan HS, McGregor M.
    Enhancement of In Vitro Interleukin-2 Production
    in Normal Subjects Following a Single Spinal
    Manipulative Treatment

    Chiropractic & Osteopathy 2008 (May 28); 16: 5

  38. Teodorczyk-Injeyan JA, McGregor M, Ruegg R, Injeyan HS.
    Interleukin 2-regulated In Vitro Antibody Production
    Following A Single Spinal Manipulative Treatment
    in Normal Subjects

    Chiropractic & Osteopathy 2010 (Sep 8); 18: 26

  39. Molina-Ortega F, Lomas-Vega R, Hita-Contreras F, et al.
    Immediate effects of spinal manipulation on
    nitric oxide, substance P and pain perception.
    Man Ther. 2014;19(5):411–7.
    doi: 10.1016/j.math.2014.02.007

  40. Brennan PC, Kokjohn K, Kaltinger CJ, Lohr GE, et al.
    Enhanced Phagocytic Cell Respiratory Burst Induced
    by Spinal Manipulation: Potential Role of Substance P

    J Manipulative Physiol Ther 1991 (Sep); 14 (7): 399–408

  41. Calder PC.
    Immunological parameters: what do they mean?
    J Nutr. 2007;137(3 Suppl 2):773S-80S.
    doi: 10.1093/jn/137.3.773S

  42. Herkenham M, Kigar SL.
    Contributions of the adaptive immune system to
    mood regulation: Mechanisms and pathways
    of neuroimmune interactions.
    Prog Neuropsychopharmacol Biol Psychiatry. 2017;79(Pt A):49–57.
    doi: 10.1016/j.pnpbp.2016.09.003

  43. El-Farhan N, Rees DA, Evans C.
    Measuring cortisol in serum, urine and saliva -
    are our assays good enough?
    Ann Clin Biochem. 2017;54(3):308–22.
    doi: 10.1177/0004563216687335

  44. Sims S-K, Wilken-Resman B, Smith CJ, Mitchell A, McGonegal L.
    Brain-Derived Neurotrophic Factor and Nerve
    Growth Factor Therapeutics for Brain Injury:
    The Current Translational Challenges
    in Preclinical and Clinical Research.
    Neural Plast. 2022;2022:3889300.
    doi: 10.1155/2022/3889300

  45. Xiong H-Y, Hendrix J, Schabrun S, Wyns A, Campenhout JV, et al.
    The Role of the Brain-Derived Neurotrophic Factor
    in Chronic Pain: Links to Central Sensitization
    and Neuroinflammation.
    Biomolecules. 2024;14(1):71.
    doi: 10.3390/biom14010071

  46. Lukkahatai N, Ong IL, Benjasirisan C, Saligan LN.
    Brain-Derived Neurotrophic Factor (BDNF) as a Marker
    of Physical Exercise or Activity Effectiveness
    in Fatigue, Pain, Depression, and Sleep
    Disturbances: A Scoping Review.
    Biomedicines. 2025;13(2):332.
    doi: 10.3390/biomedicines13020332

  47. Enrico VT, Anderst W, Bell KM, Coelho JP, Darwin J, et al.
    Plasma Pro- and Anti-Inflammatory Cytokines in an
    Observational Chronic Low Back Pain Cohort.
    JOR Spine. 2025;8(3):e70095.
    doi: 10.1002/jsp2.70095

  48. Brown MC, Kosinski AS, Fillipo R, Howell G, et al.
    Patterns and trajectories of peripheral inflammatory
    cytokines, immune tolerance, and lymphocyte
    differentiation predict transition from
    acute to chronic low back pain in
    a sex- and age-specific manner.
    Pain. 2025.
    doi: 10.1097/j.pain.0000000000003811

  49. Lima CR, Martins DF, Reed WR.
    Physiological Responses Induced by Manual Therapy
    in Animal Models: A Scoping Review

    Frontiers in Neuroscience 2020 (May 8); 14: 430

  50. d’Agostino MC, Craig K, Tibalt E, Respizzi S.
    Shock wave as biological therapeutic tool: From
    mechanical stimulation to recovery and healing,
    through mechanotransduction.
    Int J Surg. 2015;24(Pt B):147–53.
    doi: 10.1016/j.ijsu.2015.11.030

  51. Wong VW, Paterno J, Sorkin M, Glotzbach JP, Levi K, et al.
    Mechanical force prolongs acute inflammation via
    T-cell-dependent pathways during scar formation.
    FASEB J. 2011;25(12):4498–510.
    doi: 10.1096/fj.10-178087

  52. Waters-Banker C, Dupont-Versteegden EE, Kitzman PH.
    Investigating the mechanisms of massage efficacy:
    the role of mechanical immunomodulation.
    J Athl Train. 2014;49(2):266–73.
    doi: 10.4085/1062-6050-49.2.25

  53. Huang C, Holfeld J, Schaden W, Orgill D, Ogawa R.
    Mechanotherapy: revisiting physical therapy and
    recruiting mechanobiology for a
    new era in medicine.
    Trends Mol Med. 2013;19(9):555–64.
    doi: 10.1016/j.molmed.2013.05.005

  54. Waters-Banker C, Butterfield TA, Dupont-Versteegden EE.
    Immunomodulatory effects of massage on nonperturbed
    skeletal muscle in rats.
    J Appl Physiol (1985). 2014;116(2):164–75.
    doi: 10.1152/japplphysiol.00573.2013

  55. Butterfield TA, Best TM, Merrick MA.
    The dual roles of neutrophils and macrophages in
    inflammation: a critical balance between
    tissue damage and repair.
    J Athl Train. 2006;41(4):457–65.

  56. Butterfield TA, Zhao Y, Agarwal S, Haq F, Best TM.
    Cyclic compressive loading facilitates recovery
    after eccentric exercise.
    Med Sci Sports Exerc. 2008;40(7):1289–96.
    doi: 10.1249/MSS.0b013e31816c4e12

  57. Duarte FCK, Funabashi M, Starmer D, Partata WA, et al.
    Effects of Distinct Force Magnitude of Spinal
    Manipulative Therapy on Blood Biomarkers of
    Inflammation: A Proof of Principle Study
    in Healthy Young Adults.
    J Manipulative Physiol Ther. 2022;45(1):20–32.
    doi: 10.1016/j.jmpt.2022.03.012

  58. Daligadu J, Haavik H, Yielder PC, Baarbe J, Murphy B.
    Alterations in Cortical and Cerebellar Motor
    Processing in Subclinical Neck Pain Patients
    Following Spinal Manipulation

    J Manipulative Physiol Ther. 2013 (Oct); 36 (8): 527–537

  59. Brandie D, Mitchell L, Greig L, Parkinson E, Brown V, Morrissey D.
    What are small, medium and large effect sizes for
    exercise treatments of tendinopathy?
    A systematic review and meta-analysis.

  60. Holt K, Russell D, Cooperstein R, Young M, Sherson M, Haavik H.
    Interexaminer Reliability of a Multidimensional Battery
    of Tests Used to Assess for Vertebral Subluxations

    Chiropractic Journal of Australia 2016; 46 (1): 100–117

  61. Triano JJ, Budgell B, Bagnulo A, Roffey B, Bergmann T, et al.
    Review Of Methods Used By Chiropractors To Determine
    The Site For Applying Manipulation

    Chiropractic & Manual Therapies 2013 (Oct 21); 21 (1): 36

  62. Cooperstein R, Gleberzon B.
    Technique Systems in Chiropractic.
    Churchill Livingstone; page needed. 2004.

  63. Hancock MJ, Maher CG, Latimer J, McAuley JH.
    Selecting an appropriate placebo for a trial
    of spinal manipulative therapy.
    Aust J Physiother. 2006;52(2):135–8.
    doi: 10.1016/s0004-9514(06)70049-6

  64. Rosner AL.
    Evidence-based Medicine:
    Revisiting the Pyramid of Priorities

    J Bodywork & Movement Therapies 2012 (Jan); 16 (1): 42–49

  65. Ashraf M, Saeed H, Saleem Z, Rathore HA, Rasool F, et al.
    A cross-sectional assessment of knowledge, attitudes
    and self-perceived effectiveness of complementary
    and alternative medicine among pharmacy and
    non-pharmacy university students.
    BMC Complement Altern Med. 2019;19(1):95.
    doi: 10.1186/s12906-019-2503-y

  66. Gonzalez D, Jacobsen D, Ibar C, Pavan C, Monti J, et al.
    Hair Cortisol Measurement by an Automated Method.
    Sci Rep. 2019;9(1):8213.
    doi: 10.1038/s41598-019-44693-3

  67. Pohlman KA, Funabashi M, O’Beirne M, Cassidy JD, et al.
    What’s the harm? Results of an active surveillance
    adverse event reporting system for chiropractors
    and physiotherapists.
    PLoS One. 2024;19(8):e0309069.
    doi: 10.1371/journal.pone.0309069

  68. Brown VA.
    An introduction to linear mixed-effects modeling in
    R. Advances in Methods and Practices in
    Psychological Science. 2021;4(1):2515245920960351.

  69. Bates D, Mächler M, Bolker B, Walker S.
    Fitting linear mixed-effects models using lme4.
    arXiv preprint arXiv:14065823. 2014.

  70. Koller M.
    robustlmm: An R Package for Robust Estimation of
    Linear Mixed-Effects Models.
    J Stat Soft. 2016;75(6).
    doi: 10.18637/jss.v075.i06

  71. Team RC.
    R: A Language and Environment for Statistical Computing.
    2020.

  72. Van Buuren S, Groothuis-Oudshoorn K.
    mice: Multivariate Imputation by Chained Equations inR.
    J Stat Soft. 2011;45(3).
    doi: 10.18637/jss.v045.i03

  73. Kowianski P, Lietzau G, Czuba E, Waskow M.
    BDNF: A Key Factor with Multipotent Impact on
    Brain Signaling and Synaptic Plasticity.
    Cell Mol Neurobiol. 2018;38(3):579–93.
    doi: 10.1007/s10571-017-0510-4

  74. Thoenen H.
    Neurotrophins and neuronal plasticity.
    Science. 1995;270(5236):593–8.
    doi: 10.1126/science.270.5236.593

  75. Johnson IP.
    Hypothesis: upregulation of a muscle-specific
    isoform of insulin-like growth factor-1
    (IGF-1) by spinal manipulation.
    Med Hypotheses. 2008;71(5):715–21.
    doi: 10.1016/j.mehy.2008.06.038

  76. Hotta H, Watanabe N, Piché M, Hara S, Yokawa T, Uchida S.
    Non-noxious skin stimulation activates the nucleus
    basalis of Meynert and promotes NGF secretion in
    the parietal cortex via nicotinic ACh receptors.
    J Physiol Sci. 2014;64(4):253–60.
    doi: 10.1007/s12576-014-0313-z

  77. Maltese PE, Michelini S, Baronio M, Bertelli M.
    Molecular foundations of chiropractic therapy.
    Acta Biomed. 2019;90(10-S):93–102.
    doi: 10.23750/abm.v90i10-S.8768

  78. Belizário JE, Fontes-Oliveira CC, Borges JP.
    Skeletal muscle wasting and renewal:
    a pivotal role of myokine IL-6.
    Springerplus. 2016;5:619.
    doi: 10.1186/s40064-016-2197-2

  79. Mandelli F, Nüesch C, Zhang Y, Halbeisen F, et al.
    Assessing Fatty Infiltration of Paraspinal Muscles
    in Patients With Lumbar Spinal Stenosis:
    Goutallier Classification and
    Quantitative MRI Measurements.
    Front Neurol. 2021;12:656487.
    doi: 10.3389/fneur.2021.656487

  80. Hodges PW, Danneels L.
    Changes in Structure and Function of the Back
    Muscles in Low Back Pain: Different Time
    Points, Observations, and Mechanisms.
    J Orthop Sports Phys Ther. 2019;49(6):464–76.
    doi: 10.2519/jospt.2019.8827

  81. He K, Head J, Mouchtouris N, Hines K, Shea P, Schmidt R, et al.
    The Implications of Paraspinal Muscle Atrophy in
    Low Back Pain, Thoracolumbar Pathology, and
    Clinical Outcomes After Spine Surgery:
    A Review of the Literature.
    Global Spine J. 2020;10(5):657–66.
    doi: 10.1177/2192568219879087

  82. Gupta P, Sareen A, Kumar N, Boruah T, Ghasi RG, Singh DK.
    Morphological Changes in Paraspinal Muscles on
    MR Imaging in Chronic Low Back Pain Patients
    with Unilateral Lumbar Disc Herniation:
    A Cross-sectional Study.
    Journal of Clinical & Diagnostic Research. 2023;17(12).

  83. Haavik H, Kumari N, Holt K, Niazi IK, Amjad I, et al.
    The Contemporary Model of Vertebral Column Joint Dysfunction
    and Impact of High-velocity, Low-amplitude Controlled
    Vertebral Thrusts on Neuromuscular Function

    European J Applied Physiology 2021 (Oct); 121 (10): 2675–2720

  84. Anderst WJ, Gale T, LeVasseur C, Raj S, Gongaware K, Schneider M.
    Intervertebral Kinematics of the Cervical Spine
    Before, During, and After High-velocity
    Low-amplitude Manipulation

    Spine J 2018 (Dec); 18 (12): 2333–2342

  85. Vukelic S, Stojadinovic O, Pastar I, Rabach M, et al.
    Cortisol synthesis in epidermis is induced by
    IL-1 and tissue injury.
    J Biol Chem. 2011;286(12):10265–75.
    doi: 10.1074/jbc.M110.188268

  86. Adams S, Wuescher LM, Worth R, Yildirim-Ayan E.
    Mechano-Immunomodulation: Mechanoresponsive Changes
    in Macrophage Activity and Polarization.
    Ann Biomed Eng. 2019;47(11):2213–31.
    doi: 10.1007/s10439-019-02302-4

  87. Watkins AD.
    Hierarchical cortical control of
    neuroimmunomodulatory pathways.
    Neuropathol Appl Neurobiol. 1994;20(5):423–31.
    doi: 10.1111/j.1365-2990.1994.tb00992.x

  88. Teodorczyk-Injeyan JA, Injeyan HS, Ruegg R.
    Spinal Manipulative Therapy Reduces Inflammatory Cytokines
    but Not Substance P Production in Normal Subjects

    J Manipulative Physiol Ther 2006 (Jan); 29 (1): 14–21

  89. Licciardone JC, Kearns CM, Hodge LM, Bergamini MVW.
    Associations of cytokine concentrations with key
    osteopathic lesions and clinical outcomes in
    patients with nonspecific chronic low back pain:
    results from the OSTEOPATHIC Trial.
    J Am Osteopath Assoc. 2012;112(9):596–605.
    doi: 10.7556/jaoa.2012.112.9.596

  90. Licciardone JC, Kearns CM, Hodge LM, Minotti DE.
    Osteopathic manual treatment in patients with diabetes
    mellitus and comorbid chronic low back pain:
    subgroup results from the OSTEOPATHIC Trial.
    J Am Osteopath Assoc. 2013;113(6):468–78.

  91. Teodorczyk-Injeyan JA, Triano JJ, Gringmuth R, DeGraauw C.
    Effects of Spinal Manipulative Therapy on Inflammatory
    Mediators in Patients with Non-specific Low Back Pain:
    A Non-randomized Controlled Clinical Trial

    Chiropractic & Manual Therapies 2021 (Jan 8); 29 (1): 3

  92. Jang D-I, Lee A-H, Shin H-Y, Song H-R, Park J-H, Kang T-B, et al.
    The Role of Tumor Necrosis Factor Alpha (TNF-α)
    in Autoimmune Disease and Current TNF-α
    Inhibitors in Therapeutics.
    Int J Mol Sci. 2021;22(5):2719.
    doi: 10.3390/ijms22052719

  93. Feise RJ.
    Do multiple outcome measures require
    p-value adjustment?
    BMC Medical Research Methodology. 2002;2:1–4.

  94. Chaibi A, Šaltyte Benth J, Bjørn Russell M.
    Validation of Placebo in a Manual Therapy
    Randomized Controlled Trial

    Sci Rep. 2015 (Jul 6); 5: 11774

Return to SUBLUXATION

Since 1-09-2026

© 1995–2026 ~ The Chiropractic Resource Organization ~ All Rights Reserved