Skip to main content

Study protocol for patient response to spinal manipulation – a prospective observational clinical trial on physiological and patient-centered outcomes in patients with chronic low back pain



Low back pain (LBP) is a major health issue due to its high prevalence rate and socioeconomic cost. While spinal manipulation (SM) is recommended for LBP treatment by recently published clinical guidelines, the underlying therapeutic mechanisms remain unclear. Spinal stiffness is routinely examined and used in clinical decisions for SM delivery. It has also been explored as a predictor for clinical improvement. Flexion-relaxation phenomenon has been demonstrated to distinguish between LBP and healthy populations. The primary objective of the current study is to collect preliminary estimates of variability and effect size for the associations of these two physiological measures with patient-centered outcomes in chronic LBP patients. Additionally biomechanical characteristics of SM delivery are collected with the intention to explore the potential dose–response relationship between SM and LBP improvement.


This is a prospective, observational study applying side-lying, high velocity, low amplitude SM as treatment for patients with LBP over a course of 6 weeks. Approximately 80 participants will be enrolled if they present with chronic LBP of 1, 2 or 3 in Quebec Task Force Classification for spinal disorders, a Roland-Morris Disability Questionnaire (RMDQ) score ≥ 6, and persistent LBP ≥ 2 with a maximum ≥ 4 using numerical rating scale. Patient-centered outcomes include LBP using visual analog scale, RMDQ, and PROMIS-29. Lumbar spine stiffness is assessed using palpation, a hand-held instrumented device, and an automated device. Flexion-relaxation is assessed using surface electromyography at the third level of the lumbar spine. Biomechanical characteristics of SM are assessed using a self-reported, itemized description system, as well as advanced kinetic measures that will be applied to estimate forces and moments at the lumbar segment level targeted by SM.


Beside alterations in material properties of the passive components of the spine, increased neuromuscular activity may also contribute to a stiffened spine. Examining changes in both spinal stiffness and flexion-relaxation along the course of the treatment provides an opportunity to understand if the therapeutic effect of SM is associated with its action on active and/or passive components of the spine.

Trial registration

NCT01670292 on

Peer Review reports


Low back pain (LBP) is a major public health problem due to its high prevalence and staggering socioeconomic impact[15]. The pathophysiology of chronic LBP in particular is not well understood, with approximately 90% of cases categorized as either idiopathic or non-specific[6, 7]. Conservative approaches for LBP treatment include medication, exercise, and manual therapy[8]. Spinal manipulation (SM), a form of manual therapy delivered by physical therapists, osteopathic physicians and most commonly by doctors of chiropractic (DCs), has been recommended by recently published clinical practice guidelines as an effective treatment option for LBP[8, 9]. Systematic reviews of clinical trials demonstrate that SM has therapeutic effects comparable to other non-invasive treatment methods such as physical therapy and core muscle exercise[1013].

In the U.S., at least 8% of the population seeks care from DCs annually, representing approximately 190 million patient visits[14]. Among patients reporting back or neck pain, 20% seek chiropractic care[15]. Additionally, patients are highly satisfied with chiropractic care[16, 17]. In spite of its relatively widespread use by the public, the underlying therapeutic mechanisms of SM are largely unknown. Thus, it is important to investigate physiological measures that may serve both as markers of LBP severity and clinical response to SM. Two such measures are spinal stiffness and the flexion-relaxation phenomenon.

Manual therapy practitioners routinely assess spinal stiffness (e.g. the perceived spinal resistance to manually applied force) as one component of the clinical evaluation[18]. Combined with information from the overall clinical picture, DCs may elect to employ SM targeted toward regions of higher resistance with the intention to improve mobility, which in turn is believed to contribute to improved clinical outcomes[19, 20].

Studies show increased spinal stiffness using instrumented stiffness measurements in patients with LBP compared to those who are asymptomatic[21, 22]. One recent study by Fritz et al.[23] suggests that spinal stiffness assessed using an automated device may be a valuable predictor of clinical outcome for patients with LBP[23]. In particular, the authors found an immediate decrease in spinal stiffness following the first treatment and lower baseline spinal stiffness were associated with better Oswestry disability index scores after 1 week of SM treatment. However, there are conflicting reports regarding whether or not SM can reduce spinal stiffness[2328].

The flexion-relaxation phenomenon (FRP) provides information regarding the nature of neuromuscular functioning in patients with LBP[29]. It is characterized by relaxation of lumbar paraspinal muscles during full torso flexion while standing. The muscular relaxation is likely associated with the principle of energy conservation as passive tissues are thought to support and stabilize the spine at the end range position without the need for active muscle contraction[30]. The phenomenon can readily be seen by visual inspection of a plot of electromyographic (EMG) data taken from the lumbar paraspinal musculature while a participant stands erect, bends forward as far as he/she can, maintains the fully flexed position briefly, and then returns to the upright position. For asymptomatic individuals, the plot typically shows very little muscular activity while the participant stands erect, an increase in activity while bending forward due to eccentric paraspinal muscle contraction, very little activity while fully flexed, and a greater activity increase while arising due to concentric muscular contraction. Several ways of quantifying the FRP have been reported in the literature, although no one method has emerged as universally accepted[31, 32]. In general, comparisons have been made between 1) the highest muscle activity recorded while the participant is bending forward vs. the mean value while fully flexed, and 2) the highest muscle activity recorded while extending to an upright position vs. the mean value while fully flexed.

The phenomenon is most clearly seen in patients without LBP[33]. The absence of or a reduced FRP, indicating paraspinal muscle activation instead of relaxation during full torso flexion, is thought to represent neuromuscular system dysfunction, and is typical in patients with LBP[29]. Based on a review of the literature, there are indications that FRP may be a valuable objective clinical tool to aid in the diagnosis and treatment of patients with LBP[34].

Given the above, the primary objective (i.e. Specific Aim 1) of the current study is to collect preliminary estimates of variability and effect size for the associations of these two physiological measures with patient-centered outcomes in patients with chronic LBP. Study findings will aid us in planning properly-powered clinical trials that examine if baseline lumbar spine stiffness and flexion-relaxation 1) are associated with baseline patient-centered outcomes; 2) predict improvement in patient-centered outcomes after treatment with SM; and 3) improve along the course of treatment.

Secondary objectives (i.e. Specific Aim 2) will explore the therapeutic mechanism of SM from a dose–response perspective. We will obtain preliminary estimates of variability and effect size to determine if differences in spinal manipulation delivery, as estimated by thrust contact force and spinal segment load, are related to patient-centered outcomes. Clinicians use a variety of SM procedures but most commonly contact the patient with their hands while delivering quick dynamic loads at specific locations of the spine with a duration range of 100–500 milliseconds[3539]. These quickly applied dynamic forces are known as high velocity, low amplitude SM (HVLA-SM). Because considerable variation exists in HVLA-SM in terms of rate of loading, pre-load, peak force, and duration of loading, quantification and standardization of delivery may be important for maximizing therapeutic effects[38, 4042]. In the current study, DC self-reported descriptions of HVLA-SM and lumbar spine segmental load estimation methods have been developed and implemented to explore this area. Additionally, dynamic spinal stiffness during SM is explored by measuring apparent mass and driving-point impedance[43, 44].


Study population

Ethics approval was obtained from the Palmer College of Chiropractic Institution Review Board (IRB) with the Assurance Number - X2011X141. Participants are eligible if they are between 21–65 years old, meeting the classification of 1, 2, or 3 in Quebec Task Force (QTF) classification for Spinal Disorders[45, 46], present with LBP for more than 12 weeks, have a Roland Morris Disability Questionnaire (RMDQ) ≥ 6 at phone screen, and a self-reported average pain within the past 24 hours ≥ 2 at phone screen and baseline visits 1 and 2 (BL1 and BL2), and ≥ 4 at phone screen or BL1 on an 11 point numerical rating scale (NRS). A more detailed description of the inclusion/exclusion criteria is summarized in Table 1. Potential participants are recruited using a multimedia approach of newspaper ads and direct mail promotions that have been used in our previous LBP studies[4749]. Additionally press releases, media interviews, a local newspaper website, community flyers, roadside signs, Craigslist, and other online classified ads are used. All advertisement methods involved in the recruitment process have been approved by the Palmer College of Chiropractic IRB.

Table 1 Study inclusion and exclusion criteria

Study design

The current study is a prospective observational trial to evaluate the delivery of HVLA-SM using patient-centered outcomes and physiological measures over a course of 6 weeks. A total of 80 participants are projected to be recruited over a period of 20 months (September 2012 through May 2014, roughly 4 participants per month). Primary outcomes include patient-centered outcomes of back pain and function and physiological measures. Physiological measures consist of posterior-anterior stiffness of the lumbar spine using three assessment methods and flexion-relaxation. The assessment for physiological measures is implemented in the following sequence: lumbar spine stiffness measurements using hand palpation, a hand-held instrumented stiffness device, an automated stiffness device, and flexion-relaxation. The sequence was chosen to improve visit efficiency and reduce participant burden. Secondary outcomes include the descriptive record of SM delivered in each treatment visit (TV) and lumbar spine segmental load during SM delivery. Demographic and clinical characteristics of the participants are collected during BL1.

Patient visit protocol

Participant screening begins with an initial eligibility determination via phone screen. Interested and eligible participants are scheduled for an in-person BL1 consisting of an informed consent process and comprehensive clinical examination. Following the BL1 a panel renders eligibility decisions for criteria requiring clinical decision-making in a process called case review, as described in more details below. Final eligibility determination and enrollment for participants still eligible occurs during BL2. The treatment phase consists of 2 TVs per week over an approximate 6 week period. The BL2 (pre enrollment) and TV1 (post enrollment) occur on the same day. Patient-centered outcomes are collected at BL1, TV6, and TV13. The physiological assessments are performed immediately before and after SM delivery during TV1, TV5, and TV12 (6 sets of assessments in total). During these same TVs, full kinetic measures during SM delivery are captured for lumbar spine segmental load analysis. SM thrust profiles (i.e. simple kinetics including contact force and displacement and acceleration at DC’s thrusting hand) are assessed separately at TV3 and TV9. An exit interview occurs at the 13th TV following the completion of treatment and all study related activities. The detailed visit-by-visit activities are summarized in Table 2.

Table 2 Summary of clinical activities in the visit-by-visit basis

Eligibility determination

Interested individuals contact our center by phone or with a return pre-stamped postcard. Study personnel administer a short computer-assisted telephone interview to screen for provisional eligibility and, if eligible and still interested, schedule a BL1 visit. During a BL1visit, a study coordinator reviews the informed consent document, study flow chart, and specific visit activities with the participant. A short video describing study procedures is also viewed. The participant has opportunity to read the informed consent document and ask questions. Those who wish to further participate sign the written informed consent document. Vital signs, height, and weight are then measured by a study coordinator and the participant completes several research forms including patient-centered outcomes, demographics, Beck Depression Inventory, a substance assessment questionnaire, medication use, and a health history. Study coordinators review forms for incomplete data and answer computer-based questions programmed with eligibility criteria.

A DC reviews health-related research documents, noting where further information is needed to determine eligibility. The DC then conducts a focused low back diagnostic examination by first obtaining a focused LBP history. Lumbar spine X-rays and/or urinalysis may be obtained to assist in diagnosis and provide information to render a safety determination. If additional laboratory procedures or diagnostic tests are required to evaluate the participant’s LBP or health status, he or she is excluded and referred to an appropriate healthcare provider. With participant consent, other health records may be requested and reviewed to determine eligibility. Following the BL1, participants are scheduled for the BL2.

The clinician who performed the BL1 examination presents findings to a Case Review panel consisting of research clinicians, study coordinators, and at least one investigator. The Case Review panel renders eligibility determination for criteria requiring subjective and clinical decision-making. Case Review provides a formally structured eligibility determination process including at least 1 person from each of 3 roles: the consenting coordinator, research clinician, and investigator to: 1) facilitate consistent interpretation of eligibility criteria; 2) ensure participant safety; and 3) mitigate selection bias. This research team has used a similar case review process to determine eligibility in previous clinical trials[47, 49]. Following eligibility determination, a web module programmed with explicit exclusion criteria is completed for each participant. Eligible participants proceed to the BL2. Participants no longer eligible receive a phone call from the examining clinician who informs them of the determination, a summary of the exam findings and appropriate clinical recommendations.

Spinal manipulation intervention

A standardized side lying, HVLA-SM is the treatment procedure utilized. Participants lie in a lateral recumbent position with the lower leg straightened and the superior or free leg flexed at both the hip and knee and adducted across midline. The DC, standing and facing the participant, stabilizes the superior shoulder or upper arm with one hand (i.e. non-thrusting hand) as the participant’s forearms rest across the chest or abdomen. The participant’s free leg is stabilized by the DC’s thigh or lower leg (i.e. contacting leg, the same side as the thrusting hand). The manipulative load is delivered by a quick, and short controlled movement of the DC’s shoulder, arm, and hand in combination with a slight body drop. The areas where the thrusting hand contacts the participant include lumbo-pelvic tissues over or adjacent to vertebral mammillary and spinous processes, between and slightly lateral to the posterior superior iliac spine of the ilium, ischial tuberosity, and the sacrum. The manipulative load is also transmitted to participants through the clinician’s stabilizing thigh/leg, resulting in a twisting force directed to the pelvis. In the current study the thrust is only delivered using the palmar aspect of the hand. Other side lying HVLA-SM procedures utilizing other manual contacts, such as pulling on a lumbar spinous process with fingers, are not utilized. The HVLA-SM is delivered by a team of trained and experienced clinicians who have at least 15 years of clinical experience using HVLA-SM for treating patients with LBP. To enable the assessment of segmental load at the lumbosacral region during the TV1, TV5, and TV12, the patient’s upper body is restrained by an external strap system while the clinician places the non-thrusting hand on a rest located near the participant’s free shoulder[50, 51].

Primary outcome measures

Patient-centered outcomes

The primary patient-centered outcomes include the Visual Analog Scale (VAS) for LBP, the 24-item RMDQ, and the Patient Reported Outcomes Measurement Information System (PROMIS)-29 questionnaire. The RMDQ has been validated in previous studies in patients with LBP and is more responsive to change over time than most other functional status measures[5254]. The VAS is used to evaluate the worst, least, and average LBP over the past 24 hours using a 100 millimeter horizontal scale (0 = no pain; 100 = worst imaginable pain). Current pain is also monitored with the VAS at the beginning and end of each TV. The VAS has excellent metric properties, is easy to administer and score, and has received much use in LBP research[11]. The comprehensive LBP status is evaluated using the PROMIS-29 questionnaire, which is a collection of questions measuring physical function, anxiety, depression, fatigue, sleep disturbance, satisfaction with social role, pain interference and pain intensity. It aims to provide clinicians and researchers access to efficient, precise, valid, and responsive adult- and child-reported measures of health and well-being[55]. Biomechanical assessors and treating clinicians are blinded to patient-centered outcomes data.

Posterior-anterior stiffness of the lumbar spine

Lumbar spine stiffness is assessed sequentially with three approaches using hand palpation, a hand-held instrumented device, and an automated device, respectively. The hand-held instrument approach is the primary measure as it has been validated by our study team and demonstrated good reliability (0.79)[56, 57]. The same set-up allows us to implement hand palpation in a manner that mimics spinal stiffness assessments in clinical practice and will allow us to explore the association between stiffness as measured by palpation versus the hand-held instrument. Because the hand palpation approach was thought to be less intrusive, it is tested before the hand-held approach in the trial. As we finalized the protocol for study launch, a study using an automated stiffness assessment method with high within- and between-day reliabilities (0.99 and 0.98) was published and suggests that spinal stiffness may have important clinical implications[23]. We contacted the authors who loaned us the equipment for use as an added stiffness measure. However, due to the study timeline, and technical and training requirements, the automated stiffness procedure was not implemented until the clinical trial enrolled 20 participants.During stiffness assessments, participants lie prone while study clinicians mark the skin with a sterile surgical marking pen at the following locations: posterior superior iliac spines, spinous processes the seventh cervical vertebra (C7), the eleventh thoracic vertebra (T11) to the fifth lumbar vertebra (L5), interspinous spaces from T10 to the first level of sacrum (S1), and the most concave lumbar segment while lying prone (Figure 1A). During the hand palpation assessment, the participants lie face down on a custom-made chiropractic treatment table with force plates (Model 4060-NC, Bertec Inc., Columbus, OH) embedded under the cushions that support the thoracic and pelvic regions and acquire the palpatory force transmitted through the participant. The clinician applies a gentle anteriorly directed force consistent with the level typically used during examination to the vertebral levels L1 to L5. An infrared smart marker (Optotrak 3020/Certus hybrid system, Northern Digital Inc., Waterloo, Ontario, Canada) is placed on the clinician’s palpatory hand to simultaneously acquire displacement. The sampling rate for the Optotrak system is set at 100 Hz while the sampling rate for the force plates is set at 1000 Hz. No specific instruction is given to participants in terms of their breathing during the test. Figure 1B demonstrates the setup for the palpatory stiffness measurement.The hand-held stiffness device consists of a force transducer (Model # LC201-50, Omegadyne Inc., Sunbury, OH) and infrared smart markers for accurate force and displacement measurement. The assessment is performed with the participant lying prone on the treatment table. The examiner assesses tolerance of participants by pressing gradually up to a maximum of 80 N with the stiffness device in an anterior direction over the pre-marked spinous process of each lumbar vertebra from L1-L5. Participants inform the examiner if the pressure causes discomfort. When the maximum pressure of 80 N is reached an audible tone sounds as feedback to the examiner to stop compressing and withdraw the device. If pressure from the stiffness device causes discomfort, the test is reapplied more slowly and the participant indicates the level of pressure causing discomfort. If the tolerance test does not cause discomfort with 80 N of pressure at any vertebral level, the subsequent stiffness tests are conducted with a maximum pressure of 80 N. If the participant indicates discomfort during the tolerance test, the subsequent stiffness test is conducted at all vertebral levels using a maximum force value 10 N lower than the lowest value causing discomfort at any vertebral level during the tolerance test. If tolerance to the hand-held stiffness device is less than 60 N at any vertebral level, the hand-held stiffness test, as well as the automated stiffness test, is not performed. Both force and displacement are recorded using a Motion Monitor data acquisition workstation (Innovative Sport Training, Chicago, IL) at a sampling rate of 1000 Hz. The stiffness measures are performed by pressing anteriorly over each lumbar spinous process for 5 cycles at a rate of approximately 1 repetition per second. The participants are instructed to make a deep breath in, a deep breath out, and then inhale half way and hold their breath during which the test is performed. The average stiffness obtained over the last 4 cycles is used for analysis. Figure 1C demonstrates the setup for the hand-held stiffness measurement.

Figure 1
figure 1

Illustration of spinal stiffness assessment tests. Skin marks over bony landmarks of the T11-L5 spinous processes and posterior superior iliac spine (A), palpatory stiffness assessment (B), stiffness assessment with a hand-held device (C); and stiffness assessment with an automated device (D).

The automated stiffness test consists of 1) a height-adjustable treatment table equipped with safety switches to allow either the participant or examiner to immediately lower the table thereby withdrawing from the automated stiffness device; 2) an automated stiffness device that extends and applies a designated force; and 3) a computer system to program device movement based on displacement/force data sent from the device in real-time using a custom-written program in LabVIEW with a sampling rate at 200Hz (Version 8, National Instruments Corp., Austin, TX)[23, 58, 59]. The automated stiffness device is positioned over the spinous process nearest the most concave point of the lumbar curve while in the prone position, typically corresponding with the L3 spinous process. When the lumbar spine is not concave in shape (e.g. flat or convex lumbar spine) or the L2 or L4 level appears to be at similar height as L3, the automated stiffness device is placed over the L3 spinous process. Participants are instructed to lay prone, inhale then exhale normally, and hold their breath at the end of exhalation for approximately 10 seconds during which the automated stiffness test is performed. During the test, the actuator on the automated device advances at a rate of 2.5 mm per second. After reaching a preload force of 5 N on the target spinous process, the actuator holds the preload for one second before advancing further and applying a maximum load of 60 N. The maximum load is maintained for one second before the actuator retracts. To decrease participant’s burden in holding the breath, the initial position of the actuator is placed as close as possible to the target spinous process as long as the participants can breathe normally without pushing against the actuator. The procedure is conducted three times during each assessment with a testing trial performed prior to data collection to orient the participant to the procedure and to determine safety. The average stiffness of the three trials is used for analysis. Figure 1D demonstrates the setup for the automated stiffness measurement.

Data reduction: For all three stiffness measurement tests, force and displacement data are acquired simultaneously. Data acquired from the hand palpation and hand-held stiffness tests are exported into a custom-written MATLAB program (MathWorks Inc., Natick, MA) and the graphs of displacement against force are plotted. The user identifies the points corresponding to 10 N and 60 N and the software computes the slope of the curve within the region using linear regression. For the automated stiffness test, the LabVIEW program is used to compute stiffness using the same algorithm as that of the MATLAB program. It is worth noting that spinal stiffness obtained in vivo in the current study is different from the inter-segmental stiffness typically obtained in specimen studies. To distinguish between them, the term global stiffness was suggested to describe the stiffness measurement over a region of the spine containing several segments[60]. Consultation with DCs performing treatments in this study revealed that relative perceived stiffness as opposed to the absolute stiffness of a given spinal segment provides the most useful clinical information. Therefore, to be consistent with the clinical interpretation utilized by study DCs the global stiffness variation (GSV) is calculated as:

GSV=max g s 1 , g s 2 , g s 3 , g s 4 , g s 5 -min g s 1 , g s 2 , g s 3 , g s 4 , g s 5

where gsi is the global stiffness measured at the lumbar spine and i = 1,2,…,5.

To account for the between-subject variation in the global stiffness measurement (e.g. a participant may have an overall stiffer or more flexible lumbar spine), the normalized GSV (nGSV) is calculated as:

nGSV= GSV mean g s 1 , g s 2 , g s 3 , g s 4 , g s 5

Both GSV and nGSV are obtained from the palpatory and hand-held stiffness tests and will be used in data analysis.

Flexion-relaxation phenomenon

Surface EMG measurements of lumbar paraspinal muscle activity for FRP analysis in this study were initially taken with a system using Ag/AgCl electrodes built into a plastic enclosure along with a solid state preamplifier (Therapeutics Unlimited model 544, Iowa City, IA – no longer in business). Motion data were initially taken with an electronic goniometer (SG150B, Biometrics Ltd, Newport, UK) that was synchronized with the EMG system. After the first 30 participants, both measurement systems were replaced with a wireless EMG system with integrated accelerometers to provide motion data (Delsys Inc. Trigno Wireless System, Natick, MA). The major difference between the two EMG systems is that the original system had fixed-distance electrodes, whereas the current Delsys system uses two separate electrodes. The EMG and motion data are recorded using the Motion Monitor workstation with a sampling rate of 1000 Hz.FRP measurements take place pre and post treatment during TV1, TV5, and TV12. The skin over the lumbar paraspinal muscles is prepared by vigorously rubbing the skin with an alcohol wipe. The participant, seated on a stool, leans forward while disposable, self-adhesive, pre-gelled surface EMG Ag/AgCl snap electrodes (Disposable ECG electrode white foam type, 36 mm in diameter, AMBU Inc., Copenhagen, Denmark) are placed on the skin over the mid-lumbar paraspinal muscles (Figure 2). Specifically the electrodes are placed with one above and one below the L3 level on each side with approximately 1 cm vertical distance between the edges of the electrodes in a semi-flexed position such that when the participant sits erect, the electrodes do not touch. Because the electrodes are removed to enable HVLA-SM, the outline of each electrode is drawn with a skin marking pen so new ones can be placed in the same location for the post treatment measurements.

Figure 2
figure 2

Illustration of EMG electrode placement (A) and example EMG and motion signals (B). Root-mean-square (RMS) EMG data and motion data are scaled according to the maximum of RMS EMG data (B). B contains 3 cycles of flexion (B1), holding in position while fully flexed (B2), and extension (B3).

Participants are instructed to place their feet shoulder width apart with their arms hanging loosely at their side. Keeping their knees straight, they bend forward as far as they can and hold that position for 3 seconds, then return to the upright standing position. Participants perform this movement once as practice and to allow research personnel to ensure proper functioning of the system. The participant then performs the procedure 3 times in a continuous manner while the EMG and motion data are collected.

Data reduction: In order to quantify the degree to which the FRP is present, two different forms of a flexion-relaxation ratio (FRR) are used[31, 32]. One is calculated by dividing the maximum root mean square (RMS) EMG activity level during flexion (while bending forward) by the lowest mean EMG activity as measured over a one second interval during the fully flexed phase. Another FRR is similarly calculated by dividing the maximum RMS EMG activity level during extension (while returning to standing erect) by the same minimum. The beginning and end of the fully flexed phase for each cycle are determined from the plot of the motion data. A macro written in Visual Basic for Applications, within Microsoft Excel (Microsoft Corp., Redmond, WA), is used to identify the specific regions of the data plots. These are regions of full flexion, the regions in which the trunk is flexing forward, and the regions in which the trunk is extending while returning to the upright position. The macro allows the user to adjust the limits of the various regions of the plot in order to exclude occasional spurious spikes in the RMS EMG data.

In all, 12 FRRs are calculated each time the participant performs this test. A set of 3 FRRs is based on the maximum RMS EMG activity during flexion – 1 FRR for each of the 3 cycles from the left paraspinal muscles. Another set of 3 is made from the right paraspinal muscles. Similarly, 2 more sets of 3 are made: one based on the maximum RMS EMG activity levels during extension from the left and the other similarly from the right paraspinal muscles. The FRRs for the 3 cycles of each situation are averaged which yields 4 mean FRRs each time a participant performs this activity: a mean FRR from flexion and one from extension for the left and right side. Once the regions have been identified the macro performs the necessary calculations.

Secondary outcome measures

DC self-reported record of SM delivery over the course of treatment

DCs report all HVLA-SM delivered using a record log as illustrated in Figure 3. This information will allow us to better identify SM delivery patterns from the clinical perspective.

Figure 3
figure 3

Self-reported record for high velocity, low amplitude spinal manipulation. Level: the segment over which the manual contact occurs (L1-5, sacrum, and sacroiliac joint); Side Up: side of body contacted by clinician providing trunk twist motion direction; Contact Point: manual contact of the thrusting hand; SP: spinous process; MP: mammillary process; DW: interspinous space 1 cm lateral to spinous processes; PSIS: posterior superior iliac spine; Ischial Tubes: Ischial tuberosity; Fossa: midpoint between ischial tuberosity and PSIS; Base: Medial aspect of the superior sacrum; Apex: inferior sacrum; Ala: sacral ala or lateral superior sacrum; Direction of thrust: direction of thrust applied by the thrusting hand; A: anterior; S: superior; I: inferior; M: medial; L: lateral.

Contact force at the clinician’s thrusting hand

The measurement of contact force characteristics occurs as a part of the dynamic stiffness assessment (see more details below). The data reduction program used for segmental load is applied to extract variables from the contact force as measured under the clinician’s thrusting hand. These variables include the force exerted during the preload phase before the HVLA thrust, the loading rate, and peak force.

Segmental load

The method for estimating lumbar spine segmental load proposed by Triano and Schultz[50] and Triano et al.[51] was adapted for use in the current study with two modifications. Rather than using a rigid chest panel to limit the participant’s upper body movement, we stabilize the upper shoulder with a shoulder strap crossed behind the participant (Figure 4). The strap restricts the predominant motion of the participants receiving HVLA-SM, which tends to include whole-body movement in the cephalic direction and a rolling movement of the lower trunk toward the clinician. The strap design attempts to mimic the clinician’s stabilization hand (non-thrusting hand). The strap is then secured to the force plate right under the participant’s thoracic cage such that force and moment occurring during the HVLA-SM procedure are monitored. The clinician’s non-thrusting hand is placed on a moveable rest located near but not touching the participants shoulder to approximate the typical position attained during the procedure and to help the clinician maintain balance (Figure 4). The location of the lumbar segment targeted by the contact hand is traced with three Optotrak smart marker rigid bodies during the HVLA-SM procedure instead of the one-time digitization during the preload phase used by Triano et al.[51]. As a result, it is possible to compensate for the effect that lumbar spine movement may have on segmental load.

Figure 4
figure 4

Illustration of segmental load assessment during side-lying high-velocity spinal manipulation.

The method described in Triano and Schultz[50] and Triano et al.[51], however, does not consider the effect of lumbar segment orientation on segmental load. As a result, only the magnitude of the segmental force and moment will be reported. To estimate three-dimensional segmental load, the location of the joint center and its orientation of the target lumbar segment need to be determined. We adopted the method developed by Splittstoesser[61] to fulfill this task. The procedure involves the creation of individualized lumbar spine models in the standing posture using skin marks over key spinal landmarks[62]. Additionally four Optotrak smart marker rigid bodies are placed at C7, T7, T10, and S1 to estimate the orientation of the target segment using the individualized lumbar spine model (Figure 4)[61]. The limitation of Splittstoesser’s method is that it is only designed to estimate sagittal plane motion. It does not account for lumbar spine axial rotation or lateral flexion, which may also occur during HVLA-SM. Thus, investigation of the effects of segment orientation on segmental load estimation in this study is exploratory.

Data reduction: A custom-written MATLAB program is applied for segmental load data reduction. Variables extracted from the segmental load data include the force and moment during the preload phase of the HVLA-SM, loading rate, and peak force and moment.

Dynamic stiffness

To investigate dynamic spinal stiffness (i.e. the displacement, acceleration and force characteristics in time), a mixed approach is applied to measure: the movement of the doctor’s thrusting hand using a triaxial accelerometer (Model 356A17, PCB Piezotronics Inc., Depew, NY) with a sampling rate at 1000 Hz; the Optotrak motion capture system with a sampling rate at 100 Hz; and the arithmetic summation of the thrusting force over the area at the base of the palm using a thin pressure pad (Pliance system, Novel Electronics, St Paul, MN) with a sampling rate of 50 Hz. Our previous study demonstrated that a sampling rate at 50 Hz is sufficient for acquiring contact force during HVLA-SM[36]. Two characteristics of dynamic spinal stiffness are computed using the force/displacement/acceleration-time profiles: 1) apparent mass (M) that takes into account the force applied (F) and the resulting acceleration (a) of the contact interface (M = F/a); and 2) driving-point mechanical impedance (Z) that takes into account the force applied (F) and the resulting velocity (v) of the contact interface (Z = F/v). The dynamic spinal stiffness measures are obtained during TV3 and TV9.

Data collection, management and quality control

A customized submission, tracking and reporting web-based system was developed for the study. It is comprised of multiple sub-modules integrated into one comprehensive study management web application that includes sub-modules for computer assisted telephone interviewing, participant eligibility checking, participant tracking and report generation. The ASP. NET web application elements were programmed in C# and Structured Query Language (SQL) using Microsoft Visual Studio 2010 (Microsoft, Redmond, WA, USA) and Microsoft SQL Server Management Studio. User-friendly data entry interfaces were programmed with appropriate participant flow restrictions, validation schemes and skip patterns. The system uses a Project/Users Permissions System to control project personnel access to web modules. The web system is password-protected and uses a Microsoft SQL Server database platform to store all data.

Information is collected at every stage of recruitment, and throughout treatment and assessment so that the patient flow can be reported according to the Consolidated Standards of Reporting Trials guidelines[63]. We collect recruitment source, total number of responses per recruitment source, the resolution of these responses (ineligible, refused or enrolled), the number of withdrawals and reason for withdrawal, and the number of participants completing the study. For each enrolled participant, we track compliance to the treatment protocol and the assessment data that were collected at each visit.

Participant self-report questionnaires and clinician-reported SM delivery records are paper-based forms with unique participant identification numbers. Study coordinators have oversight for all paper data collection forms, log each completed form into a form tracking interface of the web system and submit data forms for key-entry weekly. These forms are double key-entered by trained data entry clerks in an MS Windows program using range and validation checks to improve accuracy and are stored in locked filing cabinets. The electronic biomechanical data are stored on a password protected network file server. Data reduction is completed within 2 weeks of data collection, transferred to the data manager in Microsoft Excel workbooks, and then uploaded to Microsoft SQL Server and stored with all other project-related data.

Data management and quality control of all data are performed using SQL views, stored procedures and real-time, web-based reports. Automated reports are viewed by the data manager to determine if quality improvement actions must occur, such as improved documentation, protocol revisions or personnel retraining. Final project datasets of combined web, paper and biomechanical data are assembled by transferring data from Microsoft SQL Server to SAS System for Windows (Release 9.3, SAS Institute Inc., Cary, NC, USA). Back-up tapes of the network drive are produced nightly.

Protection of human subjects and assessment of safety

Protection of human subjects

The study protocol was approved by the Palmer College of Chiropractic IRB.

Data and Safety Monitoring Committee (DSMC)

This trial is being monitored by an independent DSMC comprised of epidemiologists with expertise in LBP and CAM clinical trials, a biostatistician and a doctor of chiropractic. The DSMC’s role is to provide scientific and ethical oversight by evaluating the following data provided to them on a semi-annual basis: participant recruitment, accrual, and retention; adverse events (AEs); protocol deviations; data monitoring; and participant characteristics. To further monitor participant safety, serious AEs are reported to the chair of the DSMC within 5 business days. The DSMC makes recommendations to the Principal Investigator and the Office of Clinical and Regulatory Affairs at the National Center for Complementary and Alternative Medicine regarding continuation, termination, or other modifications to the study.

Adverse events

An IRB-approved AE grading and reporting protocol defines the process by which AEs are monitored and categorized for this study. This protocol also outlines when and how participant safety concerns are reported to the IRB, DSMC and funding agency. An AE is any untoward medical occurrence that may present itself during the conduct of the study and which may or may not have a causal relationship with the study procedures[64]. A serious AE is defined as any adverse event resulting in any of the following outcomes: death, a life-threatening adverse experience, hospitalization or prolongation of existing hospitalization, a persistent or significant disability/incapacity, or a congenital anomaly/birth defect resulting from a pregnancy. Participants are asked about adverse experiences at each TV and instructed to contact study personnel if they experience significant pain, discomfort or distress that they believe may be associated with treatment.

Statistical analysis

SAS System for Windows and R ( will be used for data analysis. Descriptive statistics of participant baseline characteristics and all measurements for each of the 3 assessment times will be calculated. The 3 different methods of collecting spinal stiffness will be compared with intraclass correlation coefficients. The general approach to data analysis will be to first summarize the variations in the physiological and kinetic measures using hierarchical linear regression models and then to use that information to determine the most appropriate statistical methods, including linear mixed regression models[65] and conditional linear mixed regression models[65, 66].

For Specific Aim 1, the variations in physiological measures will be summarized using a hierarchical linear regression model with a random effect for participant. If there is a non-significant random effect for participant, then this estimated random effect will be considered a sufficient summary of the physiological measures over time and linear mixed regression models will be used to examine the association between that random effect and the patient-centered outcome variables. If there is a significant random effect for participant, then there is variation in the measures over time. In this case, a linear mixed regression model with time-varying covariates will be fit to assess the association between the physiological measurements at a particular time point and patient-centered outcomes at the subsequent visit(s).

For Specific Aim 2, the variations in kinetic measures will be summarized using a hierarchical linear regression model with random effects for clinician and for participant nested within clinician. Intraclass correlation coefficients will be computed to establish if there is substantial consistency within clinician within participants and within clinician across participants and if there is sufficient variation between clinicians to allow evaluation of how variations in kinetic measures can predict patient-centered outcomes. If the random effect for participant within clinician is non-significant, or the intraclass correlation coefficient is 0.7 or higher, then the estimated random effect for clinician will be considered a sufficient summary of the delivery of care by an individual clinician. In this case, the association between this random effect and patient-centered outcomes will be evaluated using a linear mixed regression model. If there is not sufficient consistency within clinician between patients, but there is consistency within a clinician within participants, it will indicate that clinicians vary their care by participant but provide a consistent level of care across the course of treatment. In this case, the linear mixed regression model will be used to estimate the association between random effects specific to the clinician/participant combination and patient-centered outcomes. If consistency within clinicians and within participants cannot be established, this indicates substantial variation in the care provided across time by clinicians within a participant. In this case, we will fit a linear mixed regression model with time-varying covariates to assess the association between the kinetic measurements at a particular time point and patient-centered outcomes at the subsequent visit(s).

In the analysis for both specific aims 1 and 2, the complexity introduced by modeling with time-varying covariates may be better served with conditional linear mixed models. Regardless of what model is chosen, alternative models will be fit to assess the sensitivity of the statistical model to any assumptions about consistency.

Sample size

A minimum of 80 participants will be enrolled in the current study. The sample size should be sufficient to obtain preliminary estimates of variability and effect size through our planned data analysis, taking into consideration the possibility of dropouts and technical issues that may occur given the rigor and complexity of this study protocol.


The primary objective of the current study is to examine the predictive value of spinal stiffness and FRP for patient-centered outcomes in those with chronic LBP and treated with HVLA-SM. However, it also allows us to examine two domains that may contribute to the therapeutic effects of SM (i.e. neurophysiological and biomechanical processes). Analyzing both spinal stiffness and FRP measures in LBP patients undergoing SM treatment may help elucidate the contributions of passive (skeletal, ligamentous and disc structures) and active (muscles) structures during and following HVLA-SM, thus providing evidence for suspected therapeutic mechanisms[67].We hypothesize that changes in lumbar spine stiffness may be in part attributed to material changes in spinal tissues and altered paraspinal muscle activity induced by SM. Figure 5 illustrates two potential mechanistic pathways (biomechanical and neurophysiological) that are shared as putative mechanisms of LBP, whereby SM may lead to altered spinal stiffness and flexion-relaxation phenomenon in LBP patients. Knowledge gained from this study will further our understanding of the neuro-biomechanical mechanisms influenced by SM and may provide objective assessment measures for use in future clinical trials.

Figure 5
figure 5

Putative effects of low back pain (A) and spinal manipulation (B).



Adverse events


Baseline Visit 1


Baseline Visit 2


T1-T12, L1-L5, and S1-S3: Segments of cervical, thoracic, and lumbar spines, and sacrum


Doctors of chiropractic


Data and safety monitoring committee




Flexion-relaxation phenomenon


Flexion-relaxation ratio


Global stiffness variation


High velocity low amplitude spinal manipulation


Institution Review Board


Low back pain


National Institute of Health


Numerical rating scale


Patient Reported Outcomes Measurement Information System


Quebec Task Force (classification for Spinal Disorders)


Roland Morris Disability Questionnaire


Spinal manipulation




Structured Query Language


Treatment visit


Visual analog scale


  1. Dagenais S, Caro J, Haldeman S: A systematic review of low back pain cost of illness studies in the United States and internationally. Spine J. 2008, 8: 8-20. 10.1016/j.spinee.2007.10.005.

    Article  PubMed  Google Scholar 

  2. Luo X, Pietrobon R, Sun SX, Liu GG, Hey L: Estimates and patterns of direct health care expenditures among individuals with back pain in the United States. Spine. 2004, 29: 79-86. 10.1097/01.BRS.0000105527.13866.0F.

    Article  PubMed  Google Scholar 

  3. Stewart WF, Ricci JA, Chee E, Morganstein D: Lost productive work time costs from health conditions in the United States: results from the American Productivity Audit. J Occup Environ Med. 2003, 45: 1234-1246. 10.1097/01.jom.0000099999.27348.78.

    Article  PubMed  Google Scholar 

  4. Vos T, Flaxman AD, Naghavi M, Lozano R, Michaud C, Ezzati M, Shibuya K, Salomon JA, Abdalla S, Aboyans V, Abraham J, Ackerman I, Aggarwal R, Ahn SY, Ali MK, Alvarado M, Anderson HR, Anderson LM, Andrews KG, Atkinson C, Baddour LM, Bahalim AN, Barker-Collo S, Barrero LH, Bartels DH, Basanez MG, Baxter A, Bell ML, Benjamin EJ, Bennett D: Years lived with disability (YLDs) for 1160 sequelae of 289 diseases and injuries 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012, 380: 2163-2196. 10.1016/S0140-6736(12)61729-2.

    Article  PubMed  Google Scholar 

  5. Walker BF: The prevalence of low back pain: a systematic review of the literature from 1966 to 1998. J Spinal Disord. 2000, 13: 205-217. 10.1097/00002517-200006000-00003.

    Article  CAS  PubMed  Google Scholar 

  6. Nachemson AL, Jonsson E: Neck and Back Pain: The scientific evidence of causes, diagnosis and treatment. 2000, Philadelphia: Lippincott Williams & Wilkins

    Google Scholar 

  7. Pransky G, Borkan JM, Young AE, Cherkin DC: Are we making progress?: the tenth international forum for primary care research on low back pain. Spine. 2011, 36: 1608-1614. 10.1097/BRS.0b013e3181f6114e.

    Article  PubMed  Google Scholar 

  8. Chou R, Qaseem A, Snow V, Casey D, Cross JT, Shekelle P, Owens DK: Diagnosis and treatment of low back pain: a joint clinical practice guideline from the American College of Physicians and the American Pain Society. Ann Intern Med. 2007, 147: 478-491. 10.7326/0003-4819-147-7-200710020-00006.

    Article  PubMed  Google Scholar 

  9. Koes BW, Van TM, Lin CW, Macedo LG, McAuley J, McAuley J, Maher C: An updated overview of clinical guidelines for the management of non-specific low back pain in primary care. Eur Spine J. 2010, 19: 2075-2094. 10.1007/s00586-010-1502-y.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Bronfort G, Haas M, Evans RL, Bouter LM: Efficacy of spinal manipulation and mobilization for low back pain and neck pain: a systematic review and best evidence synthesis. Spine J. 2004, 4: 335-356. 10.1016/j.spinee.2003.06.002.

    Article  PubMed  Google Scholar 

  11. Goertz CM, Pohlman KA, Vining RV, Brantingham JW, Long CR: Patient-centered outcomes of high-velocity, low-amplitude spinal manipulation for low back pain: A systematic review. J Electromyogr Kinesiol. 2012, 22 (5): 670-691. 10.1016/j.jelekin.2012.03.006.

    Article  CAS  PubMed  Google Scholar 

  12. Lawrence DJ, Meeker W, Branson R, Bronfort G, Cates JR, Haas M, Haneline M, Micozzi M, Updyke W, Mootz R, Triano JJ, Hawk C: Chiropractic management of low back pain and low back-related leg complaints: a literature synthesis. J Manipulative Physiol Ther. 2008, 31: 659-674. 10.1016/j.jmpt.2008.10.007.

    Article  PubMed  Google Scholar 

  13. Rubinstein SM, van MM, Assendelft WJ, de Boer MR, van Tulder MW: Spinal manipulative therapy for chronic low-back pain: an update of a Cochrane review. Spine. 2011, 36: E825-E846. 10.1097/BRS.0b013e3182197fe1.

    Article  PubMed  Google Scholar 

  14. Barnes PM, Bloom B, Nahin RL: Complementary and alternative medicine use among adults and children: United States. Natl Health Stat Report. 2007, 2008: 1-23.

    Google Scholar 

  15. Wolsko PM, Eisenberg DM, Davis RB, Kessler R, Phillips RS: Patterns and perceptions of care for treatment of back and neck pain: results of a national survey. Spine. 2003, 28: 292-297.

    PubMed  Google Scholar 

  16. Carey TS, Garrett J, Jackman A, McLaughlin C, Fryer J, Smucker DR: The outcomes and costs of care for acute low back pain among patients seen by primary care practitioners, chiropractors, and orthopedic surgeons. N Engl J Med. 1995, 333: 913-917. 10.1056/NEJM199510053331406.

    Article  CAS  PubMed  Google Scholar 

  17. Cherkin DC, MacCornack FA: Patient evaluations of low back pain care from family physicians and chiropractors. West J Med. 1989, 150: 351-355.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. Petty NJ, Maher C, Latimer J: Manual examination of accessory movements-seeking R1. Man Ther. 2002, 7: 39-43. 10.1054/math.2001.0431.

    Article  CAS  PubMed  Google Scholar 

  19. Maitland GD, Banks K, English K, Hengeveld E: Maitland's vertebral manipulation. 2005, Oxford, UK: Butterworth-Heinemann, 7

    Google Scholar 

  20. Murtagh JE, Kenna CJ: Back pain and spinal manipulation: a practical guide. 1997, Oxford, UK: Butterworth-Heinemann, 2

    Google Scholar 

  21. Brodeur RR, DelRe L: Stiffness of the thoracolumbar spine for subjects with and without low back pain. J Neuromusculoskeletal Syst. 1999, 7: 127-133.

    Google Scholar 

  22. Latimer J, Lee M, Adams R, Moran CM: An investigation of the relationship between low back pain and lumbar posteroanterior stiffness. J Manipulative Physiol Ther. 1996, 19: 587-591.

    CAS  PubMed  Google Scholar 

  23. Fritz JM, Koppenhaver SL, Kawchuk GN, Teyhen DS, Hebert JJ, Childs JD: Preliminary investigation of the mechanisms underlying the effects of manipulation: exploration of a multivariate model including spinal stiffness, multifidus recruitment, and clinical findings. Spine. 2011, 36: 1772-1781. 10.1097/BRS.0b013e318216337d.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Allison G, Edmonston S, Kiviniemi K, Lanigan H, Simonsen AV, Walcher S: Influence of standardized mobilization on the posteroanterior stiffness of the lumbar spine in asymptomatic subjects. Physiother Res Int. 2001, 6: 145-156. 10.1002/pri.223.

    Article  CAS  PubMed  Google Scholar 

  25. Ferreira ML, Ferreira PH, Latimer J, Herbert RD, Maher C, Refshauge K: Relationship between spinal stiffness and outcome in patients with chronic low back pain. Man Ther. 2009, 14: 61-67. 10.1016/j.math.2007.09.013.

    Article  PubMed  Google Scholar 

  26. Goodsell M, Lee M, Latimer J: Short-term effects of lumbar posteroanterior mobilization in individuals with low-back pain. J Manipulative Physiol Ther. 2000, 23: 332-342. 10.1016/S0161-4754(00)90208-2.

    Article  CAS  PubMed  Google Scholar 

  27. Snodgrass SJ, Haskins R, Rivett DA: A structured review of spinal stiffness as a kinesiological outcome of manipulation: its measurement and utility in diagnosis, prognosis and treatment decision-making. J Electromyogr Kinesiol. 2012, 22: 708-723. 10.1016/j.jelekin.2012.04.015.

    Article  PubMed  Google Scholar 

  28. Stamos-Papastamos N, Petty NJ, Williams JM: Changes in bending stiffness and lumbar spine range of movement following lumbar mobilization and manipulation. J Manipulative Physiol Ther. 2011, 34: 46-53. 10.1016/j.jmpt.2010.11.006.

    Article  PubMed  Google Scholar 

  29. Mayer TG, Neblett R, Brede E, Gatchel RJ: The quantified lumbar flexion-relaxation phenomenon is a useful measurement of improvement in a functional restoration program. Spine. 2009, 34: 2458-2465. 10.1097/BRS.0b013e3181b20070.

    Article  PubMed  Google Scholar 

  30. McGill SM, Jones K, Bennett G, Bishop PJ: Passive stiffness of the human neck in flexion, extension and lateral bending. Clin Biomech. 1994, 9: 193-198. 10.1016/0268-0033(94)90021-3.

    Article  CAS  Google Scholar 

  31. Neblett R, Mayer TG, Gatchel RJ, Keeley J, Proctor T, Anagnostis C: Quantifying the lumbar flexion-relaxation phenomenon: theory, normative data, and clinical applications. Spine. 2003, 28: 1435-1446.

    PubMed  Google Scholar 

  32. Schinkel-Ivy A, Nairn BC, Drake JD: Evaluation of methods for the quantification of the flexion-relaxation phenomenon in the lumbar erector spinae muscles. J Manipulative Physiol Ther. 2013, 36: 349-358. 10.1016/j.jmpt.2013.05.017.

    Article  PubMed  Google Scholar 

  33. Geisser ME, Ranavaya M, Haig AJ, Roth RS, Zucker R, Ambroz C, Caruso M: A meta-analytic review of surface electromyography among persons with low back pain and normal, healthy controls. J Pain. 2005, 6: 711-726. 10.1016/j.jpain.2005.06.008.

    Article  PubMed  Google Scholar 

  34. Colloca CJ, Hinrichs RN: The biomechanical and clinical significance of the lumbar erector spinae flexion-relaxation phenomenon: a review of literature. J Manipulative Physiol Ther. 2005, 28: 623-631. 10.1016/j.jmpt.2005.08.005.

    Article  PubMed  Google Scholar 

  35. Christensen MG, Kollasch MW, Hyland JK: Practice Analysis of Chiropractic 2010: A project report, survey analysis, and summary of chiropractic practice in the United States. 2010, Greeley, Colorado: National Board of Chiropractic Examiners

    Google Scholar 

  36. Gudavalli MR, DeVocht J, Tayh A, Xia T: Effect of sampling rates on the quantification of forces, durations, and rates of loading of simulated side posture high-velocity, low-amplitude lumbar spine manipulation. J Manipulative Physiol Ther. 2013, 36: 261-266. 10.1016/j.jmpt.2013.05.010.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Herzog W, Conway PJ, Kawchuk GN, Zhang Y, Hasler EM: Forces exerted during spinal manipulative therapy. Spine. 1993, 18: 1206-1212. 10.1097/00007632-199307000-00014.

    Article  CAS  PubMed  Google Scholar 

  38. Triano JJ: The Mechanics of Spinal Manipulation. Clinical Biomechanics of Spinal Manipulation. Edited by: Herzog W, Herzog W. 2000, New York: Churchill Livingstone, 92-190.

    Google Scholar 

  39. Van Zoest GG, Gosselin G: Three-dimensionality of direct contact forces in chiropractic spinal manipulative therapy. J Manipulative Physiol Ther. 2003, 26: 549-556. 10.1016/j.jmpt.2003.08.001.

    Article  PubMed  Google Scholar 

  40. Brennan PC, Cramer GD, Kirstukas SJ, Cullum ME: Basic science research in chiropractic: the state of the art and recommendations for a research agenda. J Manipulative Physiol Ther. 1997, 20: 150-168.

    CAS  PubMed  Google Scholar 

  41. Cramer G, Budgell B, Henderson C, Khalsa P, Pickar J: Basic science research related to chiropractic spinal adjusting: the state of the art and recommendations revisited. J Manipulative Physiol Ther. 2006, 29: 726-761. 10.1016/j.jmpt.2006.09.003.

    Article  PubMed  Google Scholar 

  42. Herzog W: The Mechanical Neuromuscular, and Physiologic Effects Produced by Spinal Manipulation. Clinical Biomechanics of Spinal Manipulation. Edited by: Herzong W. 2000, New York: Churchill Livingstone, 191-207.

    Google Scholar 

  43. Wilder DG, Woodworth BB, Frymoyer JW, Pope MH: Vibration and the human spine. Spine. 1982, 7: 243-254. 10.1097/00007632-198205000-00008.

    Article  CAS  PubMed  Google Scholar 

  44. Mansfield NJ: Impedance methods (apparent mass, driving point mechanical impedance and absorbed power) for assessment of the biomechanical response of the seated person to whole-body vibration. Ind Health. 2005, 43: 378-389. 10.2486/indhealth.43.378.

    Article  PubMed  Google Scholar 

  45. Spitzer WO, LeBlanc FE, Dupuis M: Scientific approach to the assessment and management of activity-related spinal disorders. A monograph for clinicians. Report of the Quebec Task Force on Spinal Disorders. Spine. 1987, 12: S1-S59. 10.1097/00007632-198701000-00001.

    Article  Google Scholar 

  46. Spitzer WO, Skovron ML, Salmi LR, Cassidy JD, Duranceau J, Suissa S, Zeiss E: Scientific monograph of the Quebec Task Force on Whiplash-Associated Disorders: redefining "whiplash" and its management. Spine. 1995, 20: 1S-73S. 10.1097/00007632-199501000-00001.

    Article  CAS  PubMed  Google Scholar 

  47. Goertz CM, Salsbury SA, Vining RD, Long CR, Andresen AA, Jones ME, Lyons KJ, Hondras MA, Killinger LZ, Wolinsky FD, Wallace RB: Collaborative Care for Older Adults with low back pain by family medicine physicians and doctors of chiropractic (COCOA): study protocol for a randomized controlled trial. Trials. 2013, 14: 18-10.1186/1745-6215-14-18.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Hondras MA, Long CR, Haan AG, Spencer LB, Meeker WC: Recruitment and enrollment for the simultaneous conduct of 2 randomized controlled trials for patients with subacute and chronic low back pain at a CAM research center. J Altern Complement Med. 2008, 14: 983-992. 10.1089/acm.2008.0066.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Wilder DG, Vining RD, Pohlman KA, Meeker WC, Xia T, DeVocht JW, Gudavalli RM, Long CR, Owens EF, Goertz CM: Effect of spinal manipulation on sensorimotor functions in back pain patients: study protocol for a randomised controlled trial. Trials. 2011, 12: 161-10.1186/1745-6215-12-161.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Triano J, Schultz AB: Loads transmitted during lumbosacral spinal manipulative therapy. Spine. 1997, 22: 1955-1964. 10.1097/00007632-199709010-00003.

    Article  CAS  PubMed  Google Scholar 

  51. Triano JJ, Rogers CM, Combs S, Potts D, Sorrels K: Developing skilled performance of lumbar spine manipulation. J Manipulative Physiol Ther. 2002, 25: 353-361. 10.1067/mmt.2002.126132.

    Article  PubMed  Google Scholar 

  52. Deyo RA: Measuring the functional status of patients with low back pain. Arch Phys Med Rehabil. 1988, 69: 1044-1053.

    CAS  PubMed  Google Scholar 

  53. Riddle DL, Stratford PW, Binkley JM: Sensitivity to change of the Roland-Morris Back Pain Questionnaire: part 2. Phys Ther. 1998, 78: 1197-1207.

    CAS  PubMed  Google Scholar 

  54. Stratford PW, Binkley JM, Riddle DL: Development and initial validation of the back pain functional scale. Spine. 2000, 25: 2095-20102. 10.1097/00007632-200008150-00015.

    Article  CAS  PubMed  Google Scholar 

  55. Patient Reported Outcomes Measurement Information System (PROMIS). []

  56. Owens EF, DeVocht JW, Wilder DG, Gudavalli MR, Meeker WC: The reliability of a posterior-to-anterior spinal stiffness measuring system in a population of patients with low back pain. J Manipulative Physiol Ther. 2007, 30: 116-123. 10.1016/j.jmpt.2006.12.006.

    Article  PubMed  Google Scholar 

  57. Owens EF, DeVocht JW, Gudavalli MR, Wilder DG, Meeker WC: Comparison of posteroanterior spinal stiffness measures to clinical and demographic findings at baseline in patients enrolled in a clinical study of spinal manipulation for low back pain. J Manipulative Physiol Ther. 2007, 30: 493-500. 10.1016/j.jmpt.2007.07.009.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Kawchuk GN, Fauvel OR: Sources of variation in spinal indentation testing: indentation site relocation, intraabdominal pressure, subject movement, muscular response, and stiffness estimation. J Manipulative Physiol Ther. 2001, 24: 84-91. 10.1067/mmt.2001.112566.

    Article  CAS  PubMed  Google Scholar 

  59. Wong AY, Kawchuk G, Parent E, Prasad N: Within- and between-day reliability of spinal stiffness measurements obtained using a computer controlled mechanical indenter in individuals with and without low back pain. Man Ther. 2013, 18: 395-402. 10.1016/j.math.2013.02.003.

    Article  PubMed  Google Scholar 

  60. Stanton TR, Kawchuk GN: Reliability of assisted indentation in measuring lumbar spinal stiffness. Man Ther. 2009, 14: 197-205. 10.1016/j.math.2008.01.011.

    Article  PubMed  Google Scholar 

  61. Splittstoesser RE: Prediction of lumbar motion segment angles using trunk angle and anthropometry. Proceedings of the human factors and ergonomics society 50th annual meeting: 16–20 October, 2006. 2006, San Francisco: SAGE Publications: Human Factors and Ergonomics Society, 1279-1283.

    Google Scholar 

  62. Sicard C, Gagnon M: A geometric model of the lumbar spine in the sagittal plane. Spine. 1993, 18: 646-658. 10.1097/00007632-199304000-00018.

    Article  CAS  PubMed  Google Scholar 

  63. Schulz KF, Altman DG, Moher D: CONSORT 2010 statement: updated guidelines for reporting parallel group randomised trials. BMJ. 2010, 340: c332-10.1136/bmj.c332.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Guidance on Reporting Adverse Events to Institutional Review Boards for NIH-supported Multicenter Clinical Trials. 1999, Bethesda, MD: National Institutes of Health, []

  65. Verbeke G, Molenberghs G: Linear Mixed Models for Longitudinal Data. 2009, New York, NY: Springer, 1

    Google Scholar 

  66. Verbeke G, Spiessens B, Lesaffre E: Conditional Linear Mixed Models. Am Stat. 2001, 55: 25-34. 10.1198/000313001300339905.

    Article  Google Scholar 

  67. Bialosky JE, Bishop MD, Price DD, Robinson ME, George SZ: The mechanisms of manual therapy in the treatment of musculoskeletal pain: a comprehensive model. Man Ther. 2009, 14: 531-538. 10.1016/j.math.2008.09.001.

    Article  PubMed  Google Scholar 

Pre-publication history

Download references


The current study is in part funded by National Center for Complementary and Alternative Medicine, National Institute of Health, Grant Number 5U19AT004663 and conducted in a facility constructed with support from Research Facilities Improvement Grant Number C06 RR15433 from the National Center for Research Resources, National Institutes of Health. The study was supervised by an external advisory committee that provided their expertise in support of a variety of facets of the study from study design to technical details. The study was made possible with Mr. Dean Macken’s service and insight in fabricating necessary hardware for conducting physiological and kinetic measures. We would like to thank Dr. Julie Fritz, Department of Physical Therapy, University of Utah for loaning the automated spinal stiffness assessment equipment. We would like to thank the clinical team Drs. Paige Morgenthal, Michael Seidman, James Boysen, Amy Mathis, William Alexander, Anna Walden, Nicole Homb, Elissa Twist, Julie Hartman, Janice Kane, Ms Therese Devlin, and Ms Jennifer Toenjes for their contribution in conducting the clinical trial, and Ms Katie Hoyt for her support in recruitment. We would also like to thank Mr. Lance Corber and Mr. Gregory Boyer for their support in data and safety management and monitoring. Finally we would like to thank Mr. Varun Hariwan and Dr. Cosmin Berceanu for their assistance in data collection and reduction. All photographic images used in the current manuscript were taken from research personnel for illustration purpose and written permission was obtained for using their images.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Ting Xia.

Additional information

Competing interests

None of the authors have any competing interests arising from this research.

Authors’ contributions

CMG is the principal investigator of the parent grant of the current study. TX and DGW are the study co-leaders. TX, RMG, JWD, DGW, RDV, and CMG participated in the conception of the study. CMG, CRL, TX, DGW, RMG, JWD, RDV, GNK, and KAP participated in the design of the trial. TX, RMG, JWD, DGW, GNK, and RDV led efforts in implementing individual components of the trial. CRL developed the data analysis plan. All authors participated in the writing, and have read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xia, T., Wilder, D.G., Gudavalli, M.R. et al. Study protocol for patient response to spinal manipulation – a prospective observational clinical trial on physiological and patient-centered outcomes in patients with chronic low back pain. BMC Complement Altern Med 14, 292 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: