Skip to main content

Attainment rate as a surrogate indicator of the intervertebral neutral zone length in lateral bending: an in vitro proof of concept study

Abstract

Background

Lumbar segmental instability is often considered to be a cause of chronic low back pain. However, defining its measurement has been largely limited to laboratory studies. These have characterised segmental stability as the intrinsic resistance of spine specimens to initial bending moments by quantifying the dynamic neutral zone. However these measurements have been impossible to obtain in vivo without invasive procedures, preventing the assessment of intervertebral stability in patients. Quantitative fluoroscopy (QF), measures the initial velocity of the attainment of intervertebral rotational motion in patients, which may to some extent be representative of the dynamic neutral zone. This study sought to explore the possible relationship between the dynamic neutral zone and intervertebral rotational attainment rate as measured with (QF) in an in vitro preparation. The purpose was to find out if further work into this concept is worth pursuing.

Method

This study used passive recumbent QF in a multi-segmental porcine model. This assessed the intrinsic intervertebral responses to a minimal coronal plane bending moment as measured with a digital force guage. Bending moments about each intervertebral joint were calculated and correlated with the rate at which global motion was attained at each intervertebral segment in the first 10° of global motion where the intervertebral joint was rotating.

Results

Unlike previous studies of single segment specimens, a neutral zone was found to exist during lateral bending. The initial attainment rates for left and right lateral flexion were comparable to previously published in vivo values for healthy controls. Substantial and highly significant levels of correlation between initial attainment rate and neutral zone were found for left (Rho = 0.75, P = 0.0002) and combined left-right bending (Rho = 0.72, P = 0.0001) and moderate ones for right alone (Rho = 0.55, P = 0.0012).

Conclusions

This study found good correlation between the initial intervertebral attainment rate and the dynamic neutral zone, thereby opening the possibility to detect segmental instability from clinical studies. However the results must be treated with caution. Further studies with multiple specimens and adding sagittal plane motion are warranted.

Background

Low back pain (LBP) is a growing problem which is responsible for major population disability [1]. In the absence of a specific pathological or neurological cause, most LBP is classified as ‘non-specific’ and is often assumed to be mechanical if the pain is made better or worse by movement or position [24]. Lumbar segmental instability is thought to be an important factor in this, but for which there is no single definition or clinically available method for detection in patients [5, 6]. However, a generally accepted definition of clinical instability is “loss of the normal pattern of spinal motion causing pain and/or neurologic dysfunction” [7]. Many laboratory studies have explored this in terms of the neutral zone (NZ), which is the size of the zone of displacement when the bending moment is minimal [8, 9]. This measure has been found to be a more sensitive motion parameter in defining the onset and progression of spinal injury than the elastic zone or range of motion [10]. However, the measurement of the NZ has traditionally been impossible to obtain in vivo without invasive procedures, preventing its use in patient assessment.

A number of studies have used quantitative fluoroscopy (QF) for the measurement of inter-vertebral motion in vivo [1115]. QF provides continuous inter-vertebral motion information in both flexion-extension and lateral flexion. Patients lie passively on a robotic passive motion platform (Fig. 1a, b) which bends them at a standardised range and velocity while fluoroscopic sequences of inter-vertebral motion are obtained for measurement using image processing codes. This patient orientation minimises muscle activity and allows the intrinsic passive holding element (disc and ligament) restraints to be characterised.

Fig. 1
figure 1

QF passive motion platform: a Swung left b Swung right

QF has been used in vivo to study lumbar intervertebral motion patients and healthy controls [15]. An early version of this technology used weight bearing cineradiography and manual image registration to measure sagittal intervertebral angular motion as trunk motion progressed and claimed to be a surrogate for the NZ [16]. Later studies using fluoroscopy described this parameter as “the slope of the IVFE curve” and “the intervertebral attainment rate” [11, 13, 17].

Although most studies have concentrated on flexion-extension motion, lateral flexion has also been linked to segmental instability [18, 19]. Furthermore, lateral flexion stability has been shown both to be affected by discectomy and altered in lower limb amputees [20, 21]. Studies by our group used the ratios of the intersegmental bending gradients in the initial 10° of standardised trunk lateral flexion to express the initial attainment rate in an attempt to obtain a more standardised NZ surrogate [15] (Fig. 2).

Fig. 2
figure 2

Porcine lumbar spine testing apparatus and motion platform seen from above

Both the initial attainment rate and the NZ are expressions of intervertebral laxity. If a relationship is found to exist between them, it would provide evidence of the criterion validity of the former and demonstrate that this in vivo assessment of intrinsic lumbar segmental robustness might be used as a relatively non-invasive diagnostic tool in patients with persistent back disability where stability is in question. This study therefore sought to explore a methodology for determining this using a multi-segmented porcine lumbar spine with segments L1 to L5. The bending moments, intervertebral motion and global motion were recorded together using QF, using the same procedures as in lateral flexion QF studies of patients.

Methods

Apparatus

A fresh 5-segment porcine lumbar spine (L1 to L5) was prepared as recommended for the biomechanical testing of vertebral specimens [22] The porcine spine is said to have an anatomy that geometrically and biomechanically resembles that of the human spine [23, 24]. The paraspinal muscles were completely excised and all ligamentous components, including the interspinous ligament were preserved [24]. The specimen was preserved wrapped in saline-soaked gauze, covered in cling film and frozen for storage. It was thawed over 12 h before testing, mounted in a horizontal testing frame with the L1 and L5 vertebrae secured by metal halos and circumferential bolts. The same robotic horizontal motion platform used to provide controlled passive motion in patients receiving quantitative fluoroscopy examinations was used for testing (Atlas Clinical Ltd.) (Fig. 2). L1 was attached to the movable segment of the platform and L5 to the fixed segment.

A digital force guage (Omega Engineering Ltd DFG35-10, range 50 N, resolution 0.05 N, sampled at 125Hz) was rigidly connected to the movable part of the motion platform holding the superior vertebral segment. The motion of a connecting rod forced the specimen through a 40° arc, as applied in patient protocols [25], simultaneously transmitting continuous force data from the rod to a laptop computer. The force data were co-ordinated with the digital time stamp output of the motion platform’s motor, which moved the specimen at a uniform velocity of 6° per second at a standardised ramp-up speed over the first second of the motion. This velocity derived from the need to replicate the image recording protocol used in patients, where trade-offs on tolerance, safety and X-ray exposure led to a consensus on these settings [25].

Data collection

Fluoroscopic sequences of left and right lateral flexion were recorded at 15 frames per second over 15 s using a Siemens Arcadis Avantic fluoroscope VC10A portable C-arm fluoroscope (CE012), whose primary beam was centred on the disc space between the L3 and L4 vertebrae of the specimen. The image field included all 5 segments in all frames so that each vertebra could be tracked and the fluoroscope incorporated automatic distortion correction. Before recording the motion, a calibration image was acquired using a radiographic ruler comprising of two metallic beads of known diameter (4.4 mm) set 100 mm apart into a plastic bar and placed adjacent to the porcine spine and perpendicular to the primary-ray beam in the image field. A single fluoroscopic image was acquired so that this could be used as a scaling factor to calculate the distances between objects in the image sequences.

As in the protocol for patient recordings, the spine was preconditioned by performing four consecutive out and return lateral flexion sequences increasing from 10° up to 40° to replicate this. Ten consecutive recordings were then made of 40° left lateral flexion sequences. The spine was then replaced in a ‘neutral’ position where the force applied by the motion platform was as close to zero as possible. The same procedure was followed for right lateral flexion, however, due to the configuration of the apparatus only a maximum of 30° was achievable for right lateral flexion.

Image analysis

Outlines of the vertebral body borders of the first image were marked using the computer’s cursor in the first of each sequence of images in a manner identical to the patient mark-up protocol. The positions each of the vertebrae in each of the fluoroscopic images were calculated using automated frame to frame registration codes written in Matlab (The Mathworks Ltd. Cambridge) producing continuous tracking of each vertebral body image throughout the sequences [12]. Trackings were verified visually by a trained operator and the means of the positions of each vertebral section were generated as an output. Average angular motion was smoothed by Tikhonov regularization to reduce inter image variation as with the analysis in living subjects [26, 27].

The changing intervertebral angles of the specimen were co-ordinated with the timing and position of the motion platform. The intervertebral angles of the specimen when the motion platform reached 10°, the moments applied at each intervertebral joint and the motion platform rotation were recorded dynamically. The positions of the point of load application/measurement and the individual joint centres were derived from the trackings of each vertebra in each image frame. Since the centres of rotation between vertebrae are not generally to be found in the joint centre and due to the elasticity of the intervertebral joint, these distances varied slightly during motion and were incorporated into the continuous calculation of moments as detailed below. Forces and moments could not be measured directly at each joint, therefore estimation of forces and moments of forces were derived from the kinematics and inertial properties of the spine by applying the process of inverse dynamics. Modelling the spine as a series of free bending rods of negligible thickness and with uniform mass distribution, an estimation of forces and moments was derived based on D’Alembert’s principle (Fig. 3). One can write the Newton-Euler equations as:

$$ {{\displaystyle \sum \mathbf{F}={m}_i\mathbf{a}}}_i\iff \left(-{\mathbf{F}}_{i-1}\right)+{\mathbf{F}}_i+{m}_i\mathbf{g}={m}_i{\mathbf{a}}_i $$
(1)
$$ {\displaystyle \sum \mathbf{M}={I}_i{\boldsymbol{\upalpha}}_i}\iff \left(-{\mathbf{M}}_{i-1}\right)+{\mathbf{r}}_{i-1}\times \left(-{\mathbf{F}}_{i-1}\right)+{\mathbf{M}}_i+{\mathbf{r}}_i\times {\mathbf{F}}_i={I}_i{\boldsymbol{\upalpha}}_i $$
(2)
Fig. 3
figure 3

Mechanical model of two successive vertebrae, modelled as having negligible thickness and uniform mass distribution. The figure shows action and reaction forces, net moments of force, and all linear and angular accelerations. Gravitational forces are ignored as they are not applicable in the plane of motion

where F app is the applied force, F i is the reaction force, r i is the distance from the segment centre of mass to F i , since the geometrical centre is considered to be the centre of mass − r i is the distance from the segment centre of mass to F app , m i is the mass of segment i, g is gravity vector, α i is the angular acceleration, Ii is the moment of inertia and × represents the vector (cross) product. Since gravity is acting perpendicular to the plane of measurement it can be ignored as in Fig. 3.

From Equations 1 and 2 one can calculate each reaction force (F i ) and joint moment (M i ) acting on each vertebra of the spine.

Initial attainment rate was calculated as the ratio of the slopes of the first 10° of platform rotation and inter-vertebral rotation over the contemporaneous outward displacement of the latter (Fig. 3). If the motion segment did not rotate by at least 2.5° over this part of the motion (being twice the inter-observer error of the measurement of rotational deformation with this method) the segment was considered stiff and the initial attainment rate was not calculated [25] (Fig. 4).

Fig. 4
figure 4

Examples of initial attainment rate calculation: Gradients of inter-vertebral and platform motion in first 10 degrees of platform motion (two intervertebral levels)

The dynamic NZ was taken to be the inter-vertebral angle at the end of the region confined by a slope of +0.05 Nm/degree [28]. Samples of the force-deformation curves for all levels and directions in the specimen were examined to confirm that this was a reasonable assumption for this experiment (Fig. 5).

Fig. 5
figure 5

Example of a force deformation curve from an L3-4 motion segment undergoing left and right lateral flexion

Statistical analysis

All data were tested for normality using the Shapiro-Wilk test. The inter-vertebral angle at 10° of platform motion, the dynamic NZs and the initial attainment rates were calculated for each intervertebral level and direction. Correlations between the dynamic NZs and the initial attainment rates in each segment were determined for the pooled data (n = 52) and for left and right separately using the Spearman rank correlation coefficient for non-normally distributed data. The cut-off for statistical significance was set at a P value of 0.05.

Results

The mean (SD) (L1-5) ranges of motion for each direction, as measured on the fluoroscopic images were: left 33.5°(1.2) and right 28.3°(0.9), which represented 83 % and 94 % of platform motion respectively. The initial attainment rates for left and right lateral flexion and the pooled data are shown in Table 1.

Table 1 Median segmental initial attainment rates for left and right lateral flexion

The levels of nonparametric correlation between initial attainment rate and dynamic NZ (Fig. 5) were substantial and highly significant for left and combined left-right and moderate for right alone [29] (Table 2).

Table 2 Correlations between initial attainment rate and dynamic NZ for pooled levels (L1-2 to L4-5)

Discussion

Main result

These results are similar to previously published in vivo values for healthy human controls [15] and suggest that there is a relationship between the initial attainment rate and the dynamic NZ. The range of upper quartiles for initial attainment rate (0.204–0.413) were comparable to the upper reference ranges found in vivo (0.290–0.429) [15]. However, initial attainment rate and the dynamic NZ are not usually perfectly coincident because they do not measure the same thing; NZ reflects resistance to a pure moment and attainment rate the inter-vertebral motion velocity compared to trunk motion. Furthermore, it is not suggested that the NZ can be calculated from the initial attainment rate, but merely that they are linked in a way that would allow the order of NZ length to be determined from a set of specimens or patients based on initial attainment rate results. In this experiment, they both appear to reflect the intrinsic restraining properties of the inter-vertebral linkages, although the differences need further explanation. In addition, the 10° cut-off used historically to define initial attainment was arbitrary. A better justified calculation may be provided by considering the subsequent work of Smit et al. [30].

Learning points as an exploratory study

Some of the motion of the frame (40°) was not transferred to the vertebral segments, as 6.5° (left) and 1.5° (right) respectively were lost. This may be due to the use of retaining bolt heads into the bone, calling for a better fixation method. This may have affected the correlations. In addition, two of the segments (L1-2 and L4-5 left) did not reach the required 2.5° required for initial attainment rate to be reported (Table 1). This is likely to be a prevailing feature of multi-segmental examinations, especially if segmental levels are not challenged. Future experimental setups should ensure that equal ranges of the motion platform are obtained.

In calculating the point of inter-vertebral motion from which initial attainment rate measurement begins, fluctuations can occur. If these are prominent, the initial attainment rate value may alter and the method chosen for smoothing to obtain an average value, as well as the ramp-up speed, could affect initial attainment rate values. An international forum on the use of QF suggested this preferred smoothing function but that these values should be kept under review [2527].

Another question might be why there was not symmetry in the measurement results. The dynamic NZs were generally of a greater order for left lateral flexion, (median left = 7.17°, median right = 4.70°) but over a higher range (range left = 5.90°, range right = 6.49°). This might be the result of lower ranges of right global bending during pre-conditioning and repeated motion and/or alternatively, greater laxity at L3-4 in left lateral flexion (representing the upper cluster in Fig. 6) as a physiological variant. Further studies using multiple specimens and symmetrical testing should clarify this.

Fig. 6
figure 6

Scatter plot of dynamic NZ (degrees) against initial attainment for left and right lateral flexion

Relevance to clinical studies

In patient and volunteer research studies, the presence of a greater volume of soft tissue between the motion frame and the segment will add noise to the calculation of the initial attainment rate. It might be expected that laxity would be associated with a greater overall range of the segment, but may also be affected by the soft tissue mass. The extent of this might be explored in vivo by comparing the initial attainment rates to the overall segmental ranges obtained using QF and to body-mass index.

An additional major challenge in passive system spine kinematics research lies in the complexity of upright motion. This adds the influence of unaccounted variations arising from muscle motor control and body segment mass. However, it also extends the scope of the kind of stability parameters that can be considered.

Suggestions for further work

The present studies were limited to lateral flexion, where in some circumstances stability may be important. However, the greatest interest in stability, especially for purposes of surgical decisions, focuses on the sagittal plane, where translation is the main kinematic measure used in estimating stability [5]. Studies of the correlation between this and initial attainment rate in the sagittal plane would further inform the use of initial attainment rate in the assessment of patients for segmental laxity.

Conclusion

The ability to measure inter-vertebral laxity in patients with QF is a step forward in the assessment of chronic back pain where mechanics is thought to be important. This study used the passive recumbent QF protocol in a multi-segmental porcine model for assessing the intrinsic intervertebral responses to a minimal bending moment. It found there to be good correlation between the initial attainment rate and the dynamic NZ, thereby opening the possibility to measure passive system inter-vertebral laxity in clinical studies. However, this was an exploratory study based on repeated measurements in a single specimen, albeit a multilevel one. Therefore, the results, although likely to be important, should be treated with caution. Further, multi-specimen in vitro studies are now warranted.

References

  1. Hoy D, March L, Brooks P, Blyth F, Woolf A, Bain C, et al. The global burden of low back pain: estimates from the Global Burden of Disease 2010 study. Ann Rheum Dis. 2014;73(6):968–74.

    Article  PubMed  Google Scholar 

  2. NICE Back pain - low (without radiculopathy) - Management. Clinic Knowledge Summaries, 2015. http://cks.nice.org.uk/back-pain-low-without-radiculopathy.

  3. Deyo RA, Dworkin SF, Amtmann D, Andersson G, Borenstein D, Carragee E, et al. Report of the National Institutes of Health Task Force on Research Standards for Chronic Low Back Pain. J Manipulative Physiol Ther. 2014;37(7):449–67.

    Article  PubMed  Google Scholar 

  4. Campbell J, Colvin LA. Management of Back Pain. Br Med J. 2013;347(Dec 5):f3148.

    Article  Google Scholar 

  5. Kanemura A, Doita M, Kasahara K, Sumi M, Kurosaka M, Iguchi T. The influence of sagittal instability factors on clinical lumbar spinal symptoms. J Spinal Disord Tech. 2009;22(7):479–85.

    Article  PubMed  Google Scholar 

  6. Steiger F, Becker H-J, Standaert CJ, Nalague F, Vader J-P, Porchet F, et al. Surgery in lumbar degenerative spondylolisthesis: indications, outcomes and complications. A systematic review. Eur Spine J. 2014;23:945–73.

    Article  CAS  PubMed  Google Scholar 

  7. Panjabi MM. Clinical spinal instability and low back pain. J Electromyogr Kinesiol. 2003;13(4):371–9.

    Article  PubMed  Google Scholar 

  8. Panjabi MM. The stabilising system of the spine - Part 2: Neutral zone and instability hypothesis. J Spinal Disord. 1992;5(4):390–7.

    Article  CAS  PubMed  Google Scholar 

  9. Crawford NR, Peles JD, Dickman CA. The spinal lax zone and neutral zone: measurement techniques and parameter comparisons. J Spinal Disord. 1998;11:416–29.

    Article  CAS  PubMed  Google Scholar 

  10. Oxland TR, Panjabi MM. The onset and progression of spinal injury: A demonstration of neutral zone sensitivity. J Biomech. 1992;25(10):1165–72.

    Article  CAS  PubMed  Google Scholar 

  11. Wong KWN, Leong JCY, Chan M-K, Lu WW. The flexion-extension profile of lumbar spine in 100 healthy volunteers. Spine. 2004;29(15):1636–41.

    Article  PubMed  Google Scholar 

  12. Breen A, Muggleton J, Mellor F. An objective spinal motion imaging assessment (OSMIA): reliability, accuracy and exposure data. BMC Musculoskelet Disord. 2006;7(1):1–10.

    Article  PubMed Central  PubMed  Google Scholar 

  13. Teyhen DS, Flynn TW, Childs JD, Kuklo TR, Rosner MK, Polly DW, et al. Fluoroscopic Video to Identify Aberrant Lumbar Motion. Spine. 2007;32(7):E220–E9.

    Article  PubMed  Google Scholar 

  14. Ahmadi A, Maroufi N, Behtash H, Zekavat H, Parnianour M. Kinematic analysis of dynamic lumbar motion in patients with lumbar segmental instability using digital videofluoroscopy. Eur Spine J. 2009;18:1677–85.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Mellor F, Muggleton JM, Bagust J, Mason W, Thomas PW, Breen AC. Mid-lumbar lateral flexion stability measured in healthy volunteers by in-vivo fluoroscopy. Spine. 2009;34(22):E811–E7.

    Article  PubMed  Google Scholar 

  16. Kanayama M, Abumi K, Kaneda K, Tadano S, Ukai T. Phase Lag of the Intersegmental Motion in Flexion-Extension of the Lumbar and Lumbosacral Spine: An In Vivo Study. Spine. 1996;21(12):1416–22.

    Article  CAS  PubMed  Google Scholar 

  17. Wong K, Luk K, Leong J, Wong S, Wong K. Continuous dynamic spinal motion analysis. Spine. 2006;31(4):414–9.

    Article  PubMed  Google Scholar 

  18. Miles M, Sullivan WE. Lateral bending at the lumbar and lumbosacral joints. Anat Rec. 1961;139(3):387–98.

    Article  Google Scholar 

  19. Kirkaldy-Willis WH, Farfan HF. Instability of the lumbar spine. Clin Orthop Relat Res. 1982;165:110–23.

    PubMed  Google Scholar 

  20. Goel VK, Nishiyama K, Weinstein J, Liu YK. Mechanical properties of lumbar spinal motion segments as affected by partial disc removal. Spine. 1986;11:1008–12.

    Article  CAS  PubMed  Google Scholar 

  21. Hendershot BD. Alterations and Asymmetries in Trunk Mechanics and Neuromuscular Control among Persons with Lower-Limb Amputation: Exploring Potential Pathways of Low Back Pain. Blacksburg: Virginia Polytechnic Institute and State University; 2012.

    Google Scholar 

  22. Wilke HJ, Wenger K, Claes L. Testing criteria for spinal implants: recommendatons for the standardization of in vitro stability testing of spinal implants. Eur Spine J. 1998;7:148–54.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Busscher I, van der Veen AJ, van Dieen JH, Kingma I, Verkerke GJ, Veldhuizen AG. In Vitro Biomechanical Characteristics of the Spine: A comparison between Human and Porcine Spinal Segments. Spine. 2010;35(2):E35–42.

    Article  PubMed  Google Scholar 

  24. Tai C, Hsieh P, Chen W, Chen L, Chen W, Lai P. Biomechanical comparison of lumbar spine instability between laminectomy and bilateral laminotomy for spinal stenosis syndrome-an experimental study in porcine model. BMC Musculoskelet Disord. 2008;9:84.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Breen AC, Teyhan DS, Mellor FE, Breen AC, Wong KWN, Deitz A. Measurement of InterVertebral Motion Using Quantitative Fluoroscopy: Report of an International Forum and Proposal for Use in the Assessment of Degenerative Disc Disease in the Lumbar Spine. Adv Orthop. 2012;2012:802350.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Eilers PH. A perfect smoother. Anal Chem. 2003;75(14):3631–6. Epub 2003/10/23.

    Article  CAS  PubMed  Google Scholar 

  27. Lubansky AS, Yeow YL, Leong Y-K, Wickramasinghe SR, Han B. A general method of computing the derivative of experimental data. AIChE J. 2006;52(1):323–32.

    Article  CAS  Google Scholar 

  28. Thompson RE, Barker TM, Pearcy MJ. Defining the Neutral Zone of sheep intervertebral joints during dynamic motions: an in vitro study. Clin Biomech. 2003;18:89–98.

    Article  Google Scholar 

  29. Landis JR, Koch GG. The Measurement of Observer Agreement for Categorical Data. Biometrics. 1977;33(1):159–74.

    Article  CAS  PubMed  Google Scholar 

  30. Smit TH, van Tunen MSLM, van der Veen AJ, Kingma I, van Dieen JH. Quantifying Intervertebral disc mechanics: a new definition of the neutral zone. BMC Musculoskelet Disord. 2011;12:38.

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Alexander C Breen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contribution

AB carried out the specimen experiments, analysed the results and drafted themanuscript. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Breen, A.C., Dupac, M. & Osborne, N. Attainment rate as a surrogate indicator of the intervertebral neutral zone length in lateral bending: an in vitro proof of concept study. Chiropr Man Therap 23, 28 (2015). https://doi.org/10.1186/s12998-015-0073-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12998-015-0073-8

Keywords