Clinical Study
Stress in lumbar intervertebral discs during distraction: a cadaveric study
The Spine Journal 8 (2008) 982–990
a) Department of Physical Medicine and Rehabilitation, Mayo Clinic, 200 1st St. SW, Rochester, MN 55905, USA b) Biomechanics Laboratory, Division of Orthopedic Research, Mayo Clinic, 200 1st St. SW, Rochester, MN 55905, USA c) Northwestern Health Sciences University, 2501 W 84th St., Bloomington, MN 55431, USA
Received 12 December 2006; accepted 31 July 2007
AUTHORS: Ralph E. Gay, MD, DC *a,b; Brice Ilharreborde, MD *b; Kristin D. Zhao, MS *b; Lawrence J. Berglund, BS *b; Gert Bronfort, DC, PhD *c; Kai-Nan An, PhD *b
abstract
BACKGROUND CONTEXT: The intervertebral disc is a common source of low back pain (LBP). Prospective studies suggest that treatments that intermittently distract the disc might be beneficial for chronic LBP. Although the potential exists for distraction therapies to affect the disc biomechanically, their effect on intradiscal stress is debated.
PURPOSE: To determine if distraction alone, distraction combined with flexion, or distraction combined with extension can reduce nucleus pulposus pressure and posterior annulus compressive stress in cadaveric lumbar discs compared with simulated standing or lying.
STUDY DESIGN: Laboratory study using single cadaveric motion segments.
OUTCOME MEASURES: Strain gauge measures of nucleus pulposus pressure and compressive stress in the anterior and posterior annulus fibrosus.
METHODS: Intradiscal stress profilometry was performed on 15 motion segments during 5 simulated conditions: standing, lying, and 3 distracted conditions. Disc degeneration was graded by inspection from 1 (normal) to 4 (severe degeneration).
RESULTS: All distraction conditions markedly reduced nucleus pressure compared with either simulated standing or lying. There was no difference between distraction with flexion and distraction with extension in regard to posterior annulus compressive stress. Discs with little or no degeneration appeared to distribute compressive stress differently than those with moderate or severe degeneration.
CONCLUSIONS: Distraction appears to predictably reduce nucleus pulposus pressure. The effect of distraction therapy on the distribution of compressive stress may be dependent in part on the health of the disc. 2008 Elsevier Inc. All rights reserved.
Introduction
Low back pain (LBP) is a ubiquitous problem in developed countries. The cost of LBP to the United States economy is estimated to be more than 100 billion dollars annually [1,2]. The relationship between disc degeneration and back pain is incompletely understood. Disc degeneration is a progressive process that results in biomechanical compromise of the motion segment. Nucleus pulposus pressure decreases in proportion to the degree of degeneration in persons with chronic LBP [3]. The tensile modulus and Poisson’s ratio of the annulus fibrosus are likewise reduced [4]. As a result, annulus fibrosus fibers fail at lower loads leading to further degeneration [5] and abnormal spinal motion [6–8]. Although the course of disc degeneration cannot be predictably altered, many investigators are seeking ways to enhance disc physiology and retard or reverse degeneration.
Many treatments using traction (axial distraction) have been devised in an attempt to relieve LBP by affecting the disc and nerve roots. A meta-analysis of the traction literature concluded that, as a group, there was no evidence that traction therapies were beneficial for LBP [9].
Nonetheless, some randomized trials have suggested that chronic LBP might be relieved by traction methods [10- 13] and these treatments continue to be used in practice. The most commonly used methods are intermittent axial traction (which includes various proprietary devices issued under various trade names) and distraction manipulation. Distraction manipulation combines axial distraction with intermittent off-axis moments, usually flexion or extension. It is different than typical spinal manipulative therapy which uses a high velocity impulse during treatment. It is commonly used by chiropractors [14] as well as physical therapists and osteopathic and medical physicians.
Several mechanisms have been proposed to explain how distraction therapies might affect the disc. These include reducing nucleus pulposus pressure, changing the position of the nucleus relative to the posterior annulus, reducing posterior annulus stress, and changing the disc-nerve interface [15-17]. Although both axial distraction and distraction manipulation may temporarily reduce nucleus pulposus pressure [18,19], their effect on the distribution of stress in the disc is unknown.
The objective of this study was to determine the effect of distraction therapies (axial distraction, distraction with flexion, and distraction with extension) on vertical (compressive) and horizontal stress in anterior annulus, posterior annulus, and nucleus pulposus regions of the disc. We used the technique of intradiscal stress profilometry to estimate the stress in human cadaver motion segments under five conditions [20,21]. We hypothesized that all three forms of distraction would significantly reduce nucleus pulposus stress compared with axial loads simulating standing or sitting. We also hypothesized that distraction with flexion would reduce posterior disc stress more than axial distraction or distraction with extension. Finally, we sought to determine if degenerative discs were affected by distraction differently than relatively healthy ones.
Materials and methods
Specimens
Ten fresh, frozen (-20°C) cadaveric lumbar spines (LlSl) with mean age of 66.4 years (SD 13.8 years, range 40 to 82) were chosen for testing. These spines were screened for HIV/AIDS, Hepatitis B and C, tuberculosis, and Creutzfeldt-Jakob disease. Prospective specimens were imaged with anterior-posterior and lateral radiographs; those with severe osteoporosis, posttraumatic deformity, bone pathology, or significant anatomical anomaly were excluded. Spines with diffuse (multilevel), severe degenerative changes that might make stress profilometry testing difficult were also excluded.
Cadavers were thawed overnight in a refrigerator and Ll through S3 was removed en bloc. The iliotransverse ligaments were sacrificed in this process. Nonligamentous soft tissues were then removed leaving intact the lumbar vertebral bodies and all ligamentous structures including anterior longitudinal ligament, posterior longitudinal ligament, interspinous ligaments, intertransverse ligaments, and facet joint capsules. Each specimen was divided into either two or three separate vertebra-disc-vertebra units, yielding 25 motion segments. K-wires were placed in the vertebral bodies and facet joints and each motion segment was potted in circular acrylic fixtures using polymethylmethacrylate. A custom jig kept the superior and inferior fixtures parallel to each other and to the plane of the disc. Specimens were kept moist with saline-soaked toweling during preparation and testing.
Preliminary testing
The degree of distraction force and flexion-extension moments necessary to simulate these treatments in isolated motion segments was unknown. Studies have reported the amount of vertebral displacement occurring during traction [22-24] and distraction manipulation [16] in vivo or in cadavers. Therefore, we conducted preliminary tests with four randomly chosen lumbar motion segments to estimate the forces needed to produce similar displacements. Distraction of 90 N produced an average increase in posterior disc height of 1.8 mm. Adding a pure moment of 5 Nm to the distracted motion segments produced an average angular displacement of 4.3° inflexion and 4.5° in extension. These displacements were similar to those previously reported to occur during distraction. An axial load of 500 N was used to simulate the load on the lumbar spine during quiet standing and 300 N to simulate lying (nonweightbearing) [25].
The transducer (Model OrthoAR; Medical Measurements, Inc., Hackensack, NJ, USA) was previously shown to accurately measure positive hydrostatic pressure up to 2 MPa [26] but the linearity of measurements in the negative range had not been reported. Because negative values might be encountered during distraction, negative pressure values (from a custom calibration chamber) were plotted against the transducer output. The response was linear to -30 kPa (-225 mmHg) with R2 =0.9996. This range includes the negative pressures reported to occur during distraction therapies [ 18, 19].
Intradiscal stress profilometry technique
The technique of intradiscal stress profilometry was performed as described by McNally and Adams [20]. Measurements were obtained with a high-pressure strain gauge transducer mounted on a blunt 1.3 mmx 15 cm needle. By pulling the transducer through the disc at a constant rate, a “stress profile” is produced. Stress in the normal nucleus pulposus is isotropic (equal in all directions) but stress in the annulus is typically anisotropic. Therefore, the transducer was oriented to measure both the vertical and horizontal stress values in each specimen.
The Biomechanical testing set-up
The potted motion segments were attached to a custom testing device that could apply pure bending moments and axial compression or distraction simultaneously (Fig. 1). The lower vertebra was centered on a 6 degrees of freedom load cell (JR3; Woodland, CA, USA) and maintained in a neutral (0 moment) position with respect to the global co-ordinate system. Compression and distraction loads were applied to the upper vertebra using pneumatic actuators. Pure moments in flexion or extension were applied with a pulley apparatus fixed to the upper acrylic fixture with force supplied by pneumatic actuators. Angular displace-ment of the upper and lower fixtures (relative to the trans-verse or X axis) was measured with miniature tilt sensors with a resolution of 0.03 over their 20 range (Model CXTLA02; Crossbow Technology Inc., San Jose, CA, USA). The transducer was extracted by a stepper motor/pulley system that pulled a cable attached to the needle hub at 2 mm/second. LabVIEW software (National Instru-ments Inc., Austin, TX, USA) was used for data acquisition and to control transducer extraction. Data was collected at 30 Hz. The transducer was calibrated using a custom pres-sure chamber and a known amount of positive and negative pressure before tests.

Biomechanical testing
The 21 remaining motion segments were tested in the same manner. A preload of 300-N compression was applied for 30 minutes to expel excess fluid [20]. A 1.3 mm spinal needle with stylet (ground to a point) was then introduced into the anterior disc and advanced in the mid-sagittal plane through the posterior annulus under fluoroscopic guidance (Fig. 2). This created a track for the transducer midway between the vertebral end plates. The guide needle was removed and the blunt transducer needle with transducer inserted and oriented to measure the vertical component of stress. The first condition was then applied to the motion segment. The cable from the needle hub to a stepper motor/pulley was properly aligned and the transducer was with-drawn at 2 mm/second. The needle was then reinserted to measure the horizontal component of stress and again ex-tracted. Each of five test conditions were applied in a con-stant order: 1) axial compression 300 N (simulation of nonweightbearing or lying) [25], 2) axial compression 500 N (simulation of relaxed standing) [25], 3) axial dis-traction 90 N (simulation of axial distraction in neutral or traction), 4) axial distraction 90 N and extension 5 Nm (simulation of extension-distraction), and 5) axial distrac-tion 90 N and flexion 5 Nm (simulation of flexion-distrac-tion). There was at least 1 minute between conditions to allow for viscoelastic recovery.

Grading of disc degeneration
After testing, each disc was sectioned in sagittal and coronal planes and graded by two observers (an orthopedic spine surgeon and a rehabilitation physician) as normal (grade 1), mild, moderate, or severe (grades 2, 3, or 4, respectively) according to the scale of Adams et al. [5]. In the case of a disagreement between observers, a third observer (an orthopedic spine surgeon) determined the final grade. All graders were blinded to results of individual motion segment tests.
Data reduction and analysis
The relative stress values in the posterior, middle, and anterior disc regions were examined by partitioning the data into thirds. Because these regions could best be identified on profiles collected during compressive loading, each 500-N stress profile was reviewed to ensure that the middle third of the data was consistent with the hydrostatic region, which represented the functional nucleus pulposus [20]. The anterior and posterior thirds of the data (excluding the outermost data points with a precipitous drop in stress) were taken to represent the anterior and posterior disc regions (annulus fibrosus). Vertical and horizontal data were analyzed separately for each test condition in each motion segment. Peak vertical stress values were calculated for the anterior and posterior regions by averaging the single highest point value with the point values before and after it (an average of 3 point values).
The effect of the five conditions on regional vertical and horizontal stress values was examined using repeated measures ANOVA. When global F-tests were significant (p!.05), pairwise comparisons (contrasts) of nucleus stress (pressure) were made between the axial compression (500 N and 300 N) and each of the three distracted conditions. Because of the limited number of motion segments, the degenerative grades were collapsed into low degeneration (grades 1 and 2) and high degeneration (grades 3 and 4) groups. The effect of test condition and degeneration were evaluated using two-way ANOVA with repeated measures; generalized estimating equations were used to account for the correlation of the data within motion segments. After this, analysis using one-way ANOVA for repeated measures was performed by degenerative group. Finally, the distribution of vertical stress among the anterior, nucleus, and posterior disc regions was qualitatively examined in each of the five conditions. Analyses were carried out with SAS (SAS Institute Inc., Cary, NC). An a level of .05 (two-tailed) was used for all tests.
Results
Three motion segments from a single spine (L1–2, L3–4, and L5–S1) were excluded as a result of unexpected pathology found upon grading. Two more L5–S1 motion segments could not be tested because of difficulty in obtaining stable potting, so the one remaining L5–S1 segment was excluded. Data from the remaining 15 motion segments (9 lumbar spines) were analyzed. Distribution of disc levels was L1–2 (2), L2–3 (5), L3–4 (3), and L4–5 (5). The effect of disc level was not formally examined because of the small sample size and the risk of type II error; nonetheless, ANOVA indicated no large differences between disc levels suggesting that pooling of the levels was appropriate.
Distribution of degenerative grades was grade 1 (3), grade 2 (5), grade 3 (4), and grade 4 (3). This resulted in 8 in the low degeneration group and 7 in the high degeneration group. Only one cadaver was female. Fig. 3 shows a representative set of vertical stress profiles for five conditions recorded from a single disc with mild (grade 2) degeneration. These profiles are representative of the raw data collected.

Regional vertical and horizontal stress values
Table 1 shows the regional vertical and horizontal stress values for the five conditions for all specimens combined, low degeneration discs (n58), and high degeneration discs (n57). The regional (mean) vertical stress values for all specimens combined during each of the five conditions are shown graphically in Fig. 4. The vertical and horizontal stress values in each disc region are compared in Table 2. Vertical and horizontal values were statistically different (paired t tests) only in the anterior disc region and only for some conditions. Vertical and horizontal peak values in the posterior and anterior disc regions are also included in Table 2 but no statistical comparison was made as the peak values within a region did not always coincide with the same point value position.



Effect of distraction on disc stress measures Vertical and horizontal stress values in the nucleus pulposus were nearly the same (Table 1), suggesting that the measures from the nucleus were consistent with nucleus pressure. Nucleus pressure, posterior vertical stress, and anterior vertical stress were all significantly decreased in the three distracted conditions as compared with either 300 N or 500 N compression (pairwise post hoc contrasts, p!.001 for each comparison). This was also true when both low-and high-degeneration groups were analyzed separately with one exception. Comparison of anterior vertical stress in high degeneration discs between 300-N compression and flexion-distraction was not significant (p5.76). Axial distraction (without flexion or extension) yielded the lowest mean nucleus pressure. Compared with 300-N compression (simulated lying), nucleus pressure decreased 99% with axial distraction, 73% with extension-distraction, and 65% with flexion-distraction. Statistical analysis of the differences between disc regions was not carried out. Effect of flexion-distraction and extension-distraction on stress distribution The mean vertical stress values in each disc region of low degeneration discs (grades 1 and 2) during flexiondistraction and extension-distraction are shown in Fig. 5.

The highest mean value for both conditions was in the nucleus. There was little difference between the two conditions in any disc region. The same data for high degeneration discs (grades 3 and 4) are shown in Fig. 6. A formal statistical comparison was not made because of the small numbers in each group. Inspection of Fig. 6 suggests no statistical difference with flexion and distraction with extension. Yet, vertical stress appeared to be distributed differently in these conditions when the trend from anterior to posterior was considered. During extension-distraction of high degeneration discs, vertical stress was greater in the posterior and nucleus regions and least in the anterior region. A very different pattern was seen during flexion-distraction. The vertical stress appears to decrease from anterior to posterior suggesting a gradient.

Discussion
In this experiment, all three distraction conditions temporarily reduced nucleus pressure compared with simulated standing and lying. The largest effect was observed during axial distraction without flexion or extension which reduced pressure to near zero. Degenerated discs responded differently than relatively normal discs; they had greater temporary net reductions in nucleus pressure. Although not examined quantitatively, the distribution of stress among disc regions in normal or minimally degenerated discs (grade 1 or 2) was similar in flexion-distraction and extension- distraction. This could be a result of the nucleus being pressurized and efficiently distributing the stress. In discs with higher amounts of degeneration (grades 3 and 4), the nucleus had much less pressure when extension or flexion was introduced indicating that stress distribution may have been dependent on the moment applied to the segment. Flexion-distraction resulted in compressive stress being temporarily qualitatively lower in the posterior region compared with the nucleus and anterior regions. Conversely, extension-distraction of degenerated discs yielded similar vertical stress in all three regions.
Nucleus pulposus pressure has been used to calculate axial loads on the spine [25,27]. This is appropriate because the normal nucleus acts as a fluid with the stress being hydrostatic or isotropic (equal in all directions). As such, it is a scalar quantity that can be measured with strain gauge technology. Quantifying stress in the annulus is more problematic. Annular stress is not isotropic but anisotropic with different vertical and horizontal components [5]. Pressure and stress have the same SI unit of measure (Pascal). Although strain gauge transducers have been used to estimate stress in the annulus, it is debatable exactly what the measurements represent. Rao et al. interpreted the output from strain gauges placed in the annulus (to detect vertical stress) to be ‘‘intradiscal pressure in the axial direction’’ [28] despite the fact that pressure is nondirectional. McMillan et al. attempted to determine the validity of strain gauge transducer measures in the annulus and found the output of their transducer to be linearly proportional to the vertical force applied to the disc. They reasoned that the output was also proportional to the compressive stress perpendicular to the transducer membrane [21]. Interestingly, they found that the same calibration coefficient was applicable to liquids, nucleus pulposus, and all but the outer 2 to 4 mm of the annulus fibrosus. Although we recorded both horizontal and vertical stress components in this experiment, we were primarily interested in the ability of distraction to ‘‘unload’’ the disc, that is, reduce the vertical or compressive stress. Therefore, we have referred to the vertical measures as vertical stress. Vertical stress measures in the nucleus were essentially the same as the horizontal measures, and therefore were interpreted as nucleus pressure.
The normal lumbar nucleus is displaced anteriorly by extension and posteriorly by flexion when lying [29,30] but changes in nucleus pressure and position in degenerated discs are not as predictable [29,31,32] and degenerated discs have been noted to bulge posteriorly with extension [30,33]. Our findings are consistent with reports that degenerated discs may respond differently from healthy discs to flexion and extension [30,31,33] and extend that observation to include flexion and extension combined with distraction. The qualitative differences we observed in stress distribution between relatively healthy and degenerated discs might be because of the degenerated discs being unable to generate or maintain nucleus pressure. They may also be explained in part by anatomy. When the motion segment is extended, the facet joints contact each other and the center of rotation moves posteriorly toward the facets, causing the anterior disc space to widen. This effectively shields the posterior disc from further compression [32]. Conversely, flexion-distraction of degenerated discs may result in anterior compression and an anterior shift of the center of rotation. This appears to produce a stress distribution with the least compressive force in the posterior annulus. These observations suggest that the normal response of lumbar discs to flexion and extension is dependent to some extent on the health of the disc.
The primary mechanical theory underlying the use of distraction therapies for disc herniation is that they reduce nucleus pressure and pull peripheral nucleus tissue toward the center of the disc [34–36]. Distraction has been shown to produce temporary negative pressure in the nucleus of living patients [18]. Nucleus pressure in the present experiment became negative during axial distraction in 4 of 8 low degeneration discs but in only 1 of 7 high degeneration discs. Gudavalli et al. [19], recorded negative pressures during flexion-distraction in a whole cadaver model but we did not observe that in this study. This may have been a result of violation of the annular ‘‘seal’’ with the transducer, but that is unlikely considering the instruments used by Gudavalli et al. were similar to the ones we used. Other possible explanations include dissimilar forces used during flexion-distraction or the difference between whole cadaver and single motion segment models. Gudavalli et al. used intermittently applied, short-duration forces and continuous measurement. We measured pressures 1 to 2 minutes after the force was applied which might also explain this difference.
This study has several weaknesses that should be considered. First, a cadaver model may not accurately represent the response of the disc to loading in vivo. At this time there is no safe and acceptable method of obtaining similar in vivo measurements in humans. The age of tissue donors was generally older than persons presenting with discogenic back pain. The effects of freezing and thawing lumbar spine tissues is not thought to significantly affect the physical properties of human spine specimens [37]. Yet, dehydration and prolonged exposure to room temperatures are known to affect their material properties. The specimens in this experiment were kept moist [38] and the exposure to room temperature minimized. Our results were not likely affected by soft-tissue changes because of exposure. Second, the method we used to simulate treatments is most consistent with intermittent traction and lasting 1 to 2 minutes. It may not reflect the exact time course of stress change during shorter treatments such as distraction manipulation. Third, although the output of the transducer we used has been shown to be proportional to the applied compressive stress (perpendicular to the sensing element), it may not provide a highly accurate measure of compressive stress. Nonetheless, it provides a reasonable measure of stress change within specimens [21]. Fourth, we excluded all L5–S1 motion segments from our data. The L5–S1 segment has different ligamentous anatomy and slightly different kinematics than the other lumbar segments. Further, it can be difficult to secure and test. We did encounter difficulties with potting and as a result elected to exclude the single L5–S1 motion segment with usable data from analysis. Fifth, the results must be considered carefully in light of the small sample size and risk of error. Yet, the study was designed as a repeated measures study to maximize the power.
Our findings provide insight into the mechanical effects of distraction therapies but they do not establish a mechanism by which distraction might benefit those with back pain or sciatica because of disc injury. It is possible that the motion or change in stress results in mechanobiological events that lead to pain relief or promote disc health [39,40]. Studies using both animal and in vitro models have demonstrated that mechanical stress may play a role in the regulation of both degradative and anabolic processes in discs [41–43]. Kroeber et al. [42] using a rabbit model found that degenerated discs (created by compression) treated with distraction had restoration of disc height and histological evidence of regeneration. Although the method of producing degeneration in that model can be questioned, the results provide preliminary evidence that distraction might potentially have a beneficial affect on disc physiology. Distraction might also reduce local stress peaks in the annulus fibrosus which are thought to produce LBP [44]. Further studies are needed to establish a clear clinical benefit of distraction therapies. Additionally, studies are needed to examine the relationship between stress distribution and clinical markers of disc biology such as the degree of nucleus hydration [45].
Acknowledgments
The authors would like to thank Jonathan Bridges for technical assistance.
References
[1] Frymoyer JW, Durett CL. The economics of spinal disorders. In: Frymoyer JW, editor. 2nd ed., The adult spine: principles and practice, Vol 1 Philadelphia: Lippincott-Raven, 1997:144–9.
[2] Katz JN. Lumbar disc disorders and low-back pain: socioeconomic factors and consequences. J Bone Joint Surg Am 2006;88(Suppl 2):21–4.
[3] Sato K, Kikuchi S, Yonezawa T. In vivo intradiscal pressure measurement in healthy individuals and in patients with ongoing back problems. Spine 1999;24:2468–74.
[4] Acaroglu ER, Iatridis JC, Setton LA, Foster RJ, Mow VC, Weidenbaum M. Degeneration and aging affect the tensile behavior of human lumbar anulus fibrosus. Spine 1995;20:2690–701.
[5] Adams MA, McNally DS, Dolan P. ‘Stress’ distributions inside intervertebral discs. The effects of age and degeneration. J Bone Joint Surg Br 1996;78:965–72.
[6] Fujiwara A, Tamai K, An HS, et al. The relationship between disc degeneration, facet joint osteoarthritis, and stability of the degenerative lumbar spine. J Spinal Disord 2000;13:444–50.
[7] Okawa A, Shinomiya K, Komori H, Muneta T, Arai Y, Nakai O. Dynamic motion study of the whole lumbar spine by videofluoroscopy. Spine 1998;23:1743–9.
[8] Takayanagi K, Takahashi K, Yamagata M, Moriya H, Kitahara H, Tamaki T. Using cineradiography for continuous dynamic-motion analysis of the lumbar spine. Spine 2001;26:1858–65.
[9] Clarke JA, van Tulder MW, Blomberg SE, de Vet HC, van der Heijden GJ, Bronfort G. Traction for low-back pain with or without sciatica. Cochrane Database Syst Rev 2005; CD003010.
[10] Gudavalli MR, Cambron JA, McGregor M, et al. A randomized clinical trial and subgroup analysis to compare flexion-distraction with active exercise for chronic low back pain. Eur Spine J 2006;15: 1070–82.
[11] Ljunggren AE, Weber H, Larsen S. Autotraction versus manual traction in patients with prolapsed lumbar intervertebral discs. Scand J Rehabil Med 1984;16:117–24.
[12] Cambron JA, Gudavalli MR, Hedeker D, et al. One-year follow-up of a randomized clinical trial comparing flexion distraction with an exercise program for chronic low-back pain. J Altern Complement Med 2006;12:659–68.
[13] Sherry E, Kitchener P, Smart R. A prospective randomized controlled study of VAX-D and TENS for the treatment of chronic low back pain. Neurol Res 2001;23:780–4.
[14] Christensen MG, Kerkhoff D, Kollasch MW. Job analysis of chiropractic. Greeley, Colorado: National Board of Chiropractic Examiners, 2000.
[15] Cox JM, Feller J, Cox-Cid J. Distraction chiropractic adjusting: clinical application and outcomes of 1,000 cases. Top Clin Chiropr 1996;3:45–59.
[16] Gudavalli MR. Biomechanics research on flexion-distraction procedure. In: Cox JM, editor. Low back pain: mechanisms, diagnosis and treatment. 6th ed.. Philadelphia: Lippincott Williams & Wilkins, 1998:263–8.
[17] Knutsson E, Skoglund CR, Natchev E. Changes in voluntary muscle strength, somatosensory transmission and skin temperature concomitant with pain relief during autotraction in patients with lumbar and sacral root lesions. Pain 1988;33:173–9.
[18] Ramos G, Martin W. Effects of vertebral axial decompression on intradiscal pressure. J Neurosurg 1994;81:350–3.
[19] Gudavalli MR, Cox JM, Baker JA, Cramer GD, Patwardhan AG. Intervertebral disc pressure changes during a chiropractic procedure. Adv Bioeng 1997;36:215–6.
[20] McNally DS, Adams MA. Internal intervertebral disc mechanics as revealed by stress profilometry. Spine 1992;17:66–73.
[21] McMillan DW, McNally DS, Garbutt G, Adams MA. Stress distributions inside intervertebral discs: the validity of experimental ‘‘stress profilometry’’. Proc Inst Mech Eng [H] 1996;210:81–7.
[22] Twomey LT. Sustained lumbar traction. An experimental study of long spine segments. Spine 1985;10:146–9.
[23] Awad EA. Effects of pelvic traction on the intervertebral disc, a radiological study. Arch Phys Med Rehabil 1988;69:785.
[24] Mathews JA. The effects of spinal traction. Physiotherapy 1972;58: 64–6.
[25] Nachemson AL. Disc pressure measurements. Spine 1981;6:93– 7.
[26] McNally DS, Adams MA, Goodship AE. Development and validation of a new transducer for intradiscal pressure measurement. J Biomed Eng 1992;14:495–8.
[27] Wilke HJ, Neef P, Caimi M, Hoogland T, Claes LE. New in vivo measurements of pressures in the intervertebral disc in daily life. Spine 1999;24:755–62.
[28] Rao RD, Wang M, Singhal P, McGrady LM, Rao S. Intradiscal pressure and kinematic behavior of lumbar spine after bilateral laminotomy and laminectomy. Spine J 2002;2:320–6.
[29] Fennell AJ, Jones AP, Hukins DW. Migration of the nucleus pulposus within the intervertebral disc during flexion and extension of the spine. Spine 1996;21:2753–7.
[30] Beattie PF, Brooks WM, Rothstein JM, et al. Effect of lordosis on the position of the nucleus pulposus in supine subjects: a study using magnetic resonance imaging. Spine 1994;19:2096–102.
[31] Merriam WF, Quinnell RC, Stockdale HR, Willis DS. The effect of postural changes on the inferred pressures within the nucleus pulposus during lumbar discography. Spine 1984;9:405–8.
[32] Adams MA, May S, Freeman BJ, Morrison HP, Dolan P. Effects of backward bending on lumbar intervertebral discs. Relevance to physical therapy treatments for low back pain. Spine 2000;25: 431–7.
[33] Schnebel BE, Simmons JW, Chowning J, Davidson R. A digitizing technique for the study of movement of intradiscal dye in response to flexion and extension of the lumbar spine. Spine 1988;13:309–12.
[34] Cox JM, Hazen LJ, Mungovan M. Distraction manipulation reduction of an L5-S1 disk herniation. J Manipulative Physiol Ther 1993;16: 342–6
[35] Quellette JP. Low back pain: an orthopedic medicine approach. Can Fam Physician 1987;33:693–4.
[36] Cyriax J. 8th ed.In: Textbook of orthopaedic medicine, Vol 1 London: Bailliere Tindall, 1982.
[37] Panjabi MM, Krag M, Summers D, Videman T. Biomechanical timetolerance of fresh cadaveric human spine specimens. J Orthop Res 1985;3:292–300.
[38] Pflaster DS, KragMH, Johnson CC, Haugh LD, PopeMH. Effect of test environment on intervertebral disc hydration. Spine 1997;22:133–9.
[39] Iatridis JC, MacLean JJ, Roughley PJ, Alini M. Effects of mechanical loading on intervertebral disc metabolism in vivo. J Bone Joint Surg Am 2006;88(Suppl 2):41–6.
[40] MacLean JJ, Lee CR, Grad S, Ito K, Alini M, Iatridis JC. Effects of immobilization and dynamic compression on intervertebral disc cell gene expression in vivo. Spine 2003;28:973–81.
[41] MacLean JJ, Lee CR, Alini M, Iatridis JC. The effects of short-term load duration on anabolic and catabolic gene expression in the rat tail intervertebral disc. J Orthop Res 2005;23:1120–7.
[42] Kroeber M, Unglaub F, Guehring T, et al. Effects of controlled dynamic disc distraction on degenerated intervertebral discs: an in vivo study on the rabbit lumbar spine model. Spine 2005;30: 181–7.
[43] Walsh AJ, Lotz JC. Biological response of the intervertebral disc to dynamic loading. J Biomech 2004;37:329–37.
[44] McNally DS, Shackleford IM, Goodship AE, Mulholland RC. In vivo stress measurement can predict pain on discography. Spine 1996;21: 2580–7.
[45] Han SM, Lee SY, Cho MH, Lee JK. Disc hydration measured by magnetic resonance imaging in relation to its compressive stiffness in rat models. Proc Inst Mech Eng [H] 2001;215:497–501