Biomechanical analysis and modeling of different vertebral growth patterns in adolescent idiopathic scoliosis and healthy subjects
© Shi et al; licensee BioMed Central Ltd. 2011
Received: 11 March 2011
Accepted: 23 May 2011
Published: 23 May 2011
The etiology of AIS remains unclear, thus various hypotheses concerning its pathomechanism have been proposed. To date, biomechanical modeling has not been used to thoroughly study the influence of the abnormal growth profile (i.e., the growth rate of the vertebral body during the growth period) on the pathomechanism of curve progression in AIS. This study investigated the hypothesis that AIS progression is associated with the abnormal growth profiles of the anterior column of the spine.
A finite element model of the spinal column including growth dynamics was utilized. The initial geometric models were constructed from the bi-planar radiographs of a normal subject. Based on this model, five other geometric models were generated to emulate different coronal and sagittal curves. The detailed modeling integrated vertebral body growth plates and growth modulation spinal biomechanics. Ten years of spinal growth was simulated using AIS and normal growth profiles. Sequential measures of spinal alignments were compared.
(1) Given the initial lateral deformity, the AIS growth profile induced a significant Cobb angle increase, which was roughly between three to five times larger compared to measures utilizing a normal growth profile. (2) Lateral deformities were absent in the models containing no initial coronal curvature. (3) The presence of a smaller kyphosis did not produce an increase lateral deformity on its own. (4) Significant reduction of the kyphosis was found in simulation results of AIS but not when using the growth profile of normal subjects.
Results from this analysis suggest that accelerated growth profiles may encourage supplementary scoliotic progression and, thus, may pose as a progressive risk factor.
Keywordsfinite element model growth profile of the vertebral body adolescent idiopathic scoliosis bone growth modulation scoliosis pathomechanism
Adolescent idiopathic scoliosis (AIS) is a 3D spinal deformity with unknown etiology . Often, spinal column overgrowth during the peripubertal period is observed in AIS patients [2, 3]. Correspondingly, others reported scoliotic spines to be longer than control subjects (particularly in the thoracic segments) , progression of scoliotic spinal deformity occurs during the adolescent growth spurt [5–7], and curve progression is correlated with the rapid spinal growth period . Adolescents with the most common type of thoracic scoliosis were also found to be taller, leaner, and with hypokyphotic thoracic spines when compared to normal subjects [9, 10]. In particular, the anterior spinal column was found to have relative overgrowth in AIS over normal subjects . MRI studies have further confirmed the presence of longer vertebral column lengths both in AIS with thoracic or thoracolumbar curves without any corresponding changes in spinal cord length [12, 13].
Many studies have reported significant differences in the pattern of growth and growth velocity between AIS and normal adolescents [9, 10, 14]. The mean age and the magnitude of peak sitting height growth velocity were also found to differ significantly between girls that finally progressed to scoliosis and those that did not . Hägglund et al. observed above average height in scoliotic girls two years before the onset of the pubertal growth spurt . In addition, radiographs of 274 AIS patients between the age of 6.5~18.5 compared to 212 age-matched controls demonstrated an early start and later cessation of the pubertal spinal growth spurt in AIS patients . Stokes also documented a different growth profile in AIS patients compared to controls .
Based on the Hueter-Volkmann law for bone growth modulation, the "vicious cycle" qualitatively explained the mechanism of scoliotic progression in an iterative manner: the asymmetrical stress distribution leads to asymmetrical growth, which in turn causes the vertebral wedging and contributes to the spinal deformity . Stokes quantitatively modeled the effect of loading asymmetry in scoliotic spines on the rate of scoliotic progression to confirm the plausibility of the "vicious cycle" principle . Plaats et al. and Azegami et al. simulated the 'buckling' effect on the progression of scoliosis and showed that, on its own, buckling will not initiate scoliosis [17, 18].
Finite element modeling (FEM) is an effective and objective technique that allows the direct investigation of variables of interest and can be used to test different pathomechanical hypotheses [19–22]. Villemure et al. tested the contribution of different pathogenesis hypotheses related to initial asymmetrical loads in scoliotic progression . Huynh et al. demonstrated that the asymmetry of pedicle growth rate alone will contribute neither to the initiation nor the progression of the scoliotic deformity . Driscoll et al. tested the influences of concave-convex biases on the progression of scoliotic curves using a FEM integrating the anterior spine and a detailed representation of growth physiology and dynamics , and found that concave-convex biases are potential factors that influence the progression of scoliotic curves. Until now, proper biomechanical modeling has not been used to study in depth the influence of the abnormal growth profile on the pathomechanism of curve progression in AIS.
The purpose of this study is to explore the hypothesis that the progression of AIS curve deformity, during the peripubertal period, may result from abnormal differential growth profiles of the vertebral column in AIS when compared to normal adolescent controls.
Finite Element Model
Material properties of finite element model
Young's Modulus (MPa)
Newly Formed Bone
The spine models loading was set according to the data as reported by Schultz i.e., the spine loading increases from 14% to 57% of the body weight along the spine from T1 to L5, with the 2.6% body weight increase between succeeding vertebrae. Two distinct spinal loading techniques were programmed. Under the gravitational approach, the loading direction was maintained axially to simulate forces which, when coupled with the selected boundary conditions (T1 restricted in the transverse plane and L5 limited in all degrees of freedom) provided appropriate spinal stability. Another loading technique, the "follower load", was alternatively simulated in a fashion that the resultant forces from cumulative loads on each vertebra was maintained tangential to the curvature of the spine in the sagittal plane and remained axial in the coronal plane. This follower load type proved to allow improved spinal stability under ex-vivo spinal loading  and previously utilized in spinal finite element analyses . Therefore, the vector direction gravitational loading approach remained constant throughout the iterative process described below whereas the direction follower load method regulated itself in order to maintain sagittal tangential loading under each numerical cycle. Both these approaches were explored separately. The model was built and computed using ANSYS 11.0 finite element package (ANSYS Inc., USA).
Model of spinal growth
In order to represent the global result of the growth mechanism in each individual vertebra, the local deformation of the elements of the newly formed bone layer within the growth plate was simulated as in previous finite element analyses of scoliotic spines . An iteration of vertebral body growth consists of four consecutive steps, namely applying forces, measuring stresses in each element of the growth plates, calculating growth, and updating the geometry. The Cobb angles, kyphosis angle, and the lordosis angle were output after each iteration.
The validity of the developed modeling platform to comply with scoliotic progression was explored using patient data. Three patients were selected with different curve types: Lenke type-1A, Lenke type-2A, and Lenke type-3C, with no significant alteration in sagittal spinal alignment i.e., kyphotic curves between 20 and 35 with less than 5 degrees modification over time. These patients previously underwent an annual radiographic follow up of 3, 2, and 2 years respectively. The formerly described simulation methods were performed utilizing regular adolescent growth rates (G m ) of 0.8 mm/year and 1.1 mm/year in thoracic and lumbar spines respectively . Starting from initial patient curves, the model was constructed and its ability to corroborate with patient data was deemed successful if curve patterns were replicated within 5 degrees for the Cobb angles.
Simulation of Different Vertebral Column Growth Patterns
Several steps were undertaken to ensure that the model corroborated with reality while behaving in a robust fashion. The growth algorithm utilized in this analysis was acquired from in-vivo experimentation  while its application utilizing finite element analysis to explore progressive scoliotic spines has previously been demonstrated [21, 23]. The longitudinal stresses measured in the intervertebral disc L5 showed agreement with in-vivo measurements . Furthermore, in order to explore the influence of the adopted numerical assumptions on the results, several sensitivity analyses were performed. These additional simulations explored the influence of the selected growth constant (β = 0.4 to 0.6 MPa-1, a range of plausible physiologic values ), the loading configuration (gravitational and follower type spinal loading), and the magnitude of the growth velocities (G m = ± 15% of values reported in Figure 3).
However, when an initial coronal deformity was present, both AIS and normal growth profiles resulted in increased lateral deformity. This progressive trend was significantly amplified when using the AIS growth profile. More specifically, after ten years of simulated growth cases 3, 4, 5, and 6 underwent Cobb angle increases with magnitudes respectively measured at 5.0, 4.3, 5.2, and 3.4 times larger using AIS growth profile compared to the normal one.
In addition to the augmented increase in the lateral deformity under the AIS growth profile, another finding suggest that the presence of an initial lateral deformity encourages a decrease in kyphosis under AIS growth profile while no such trends were observed with the normal growth profile. This observation held true despite the variety of explored initial kyphosis angles. Conversely, no significant changes were found in the lordosis angle of all the six cases.
The sensitivity analysis of β (0.4 to 0.6 MPa-1) magnified spinal Cobb angle amplitudes as β was increased. The final angle measures (after 10 years of simulated growth) were roughly doubled under a β of 0.6 compared to 0.4. However, relative differences between final measures of AIS and healthy growth rates remained as described above. That is, the relative increase in coronal Cobb (CobbAIS growth/Cobbnormal growth after 10 years of simulated growth), for cases 3 to 6, varied lightly between 1.85 to 2.61 and 1.79 to 2.86 when using a beta of 0.4 and 0.6 respectively. Therefore, it was consistently maintained that the AIS growth rate significantly encouraged additional scoliotic progression. Analysis of spinal loading (gravity force or sagittal plane follower load) also proved not to significantly alter the tendency of AIS growth to promote progression. To elaborate, the average Cobb angle increase initiated by the use of the AIS growth profile was between 1.85 to 2.61 and 1.50 to 1.88 when adopting gravity and sagittal follower loads respectively. Analysis of growth velocities (G m = ± 15%) altered magnitude of measured Cobb angles, to a lesser extent than factor β, and again, the relative comparisons were not significantly changed. More specifically, even when a 15% decrease in AIS growth rates was coupled with a 15% increase in normal growth velocities, final coronal Cobb angle related to AIS growth remained 1.8 times larger than that of the normal growth.
This study is the first to explore the influence of different vertebral growth patterns on AIS progression using state-of-the-art biomechanical modeling techniques. This research utilizes the "vicious cycle" notion of scoliotic progression under different spinal growth rates (accelerated AIS and normal). Simulation results of the finite element models, which were reaffirmed via sensitivity analyses, suggest that when an initial deformity is present, a faster AIS growth profile significantly encourages scoliotic progression in the coronal plane and decreases kyphosis in the sagittal plane. This result is consistent with the observations made between the height and angle velocities in AIS patients  and agrees with the tendency for scoliotic patients to adopt a reduced kyphosis [9, 10]. Result from Case 2 also suggests that the presence of a small kyphosis angle cannot lead to scoliosis exclusive of an initial coronal deformity. Therefore, results suggest that the abnormal growth pattern of the anterior spine may play a secondary instead of a primary role in the development of AIS.
To investigate if the spinal deformity was incurred by the axial rotate on of each vertebral body, axial rotation angle of every vertebral body, with respect to a fixed global reference plane, was analyzed. Minor axial rotation (less than 5 degrees) was measured. This observation is consistent with the simulation results reported in the literature [21, 38]. Therefore, it is perceivable that other mechanisms are involved in transverse plane deformations that include vertebral rotation, axial torsion, and rib hump.
A potential limitation of this analysis resides with the mere modeling of anterior spine growth. However, in a preceding study , it was shown that pedicle growth rate asymmetry (neuro-central growth plate) was neither able to independently generate a scoliosis nor to act in conjunction with other deformations to initiate scoliotic spinal curves. In addition, the cylindrical shape used to model the vertebral bodies may influence local stress distributions over the growth plates. However, this study seeks to draw comparative conclusions utilizing identical platforms while altering a single variable of interest (i.e. the growth rate) and, for the purpose of this analysis, such a factor was not deemed important. Moreover, assumptions that were considered to have important influence on the reported results were explored under complementing sensitivity analyses. Therefore, although the implemented numerical approach contains simplifications, sensitivity analyses suggest that adopted numerical techniques do not interfere with relative conclusion reported herein. Finally, as always in biomechanics, finite element modeling is a technique for simulating a mechanism of interest performed under logical assumptions rather than completely reconstructing reality. Therefore, results and conclusions of this study should be interpreted within the prescribed conditions. These modeling simplifications are not able to fully account for the functional limitations of the posterior elements and the coupling between loads in the different directions. However, resulting contact forces on facet joints might be more important in loading modes such as torsion, flexion/extension, and lateral bending as compared to compression, which may modify transmitted spinal loads. Thus, although not explored, these contact forces might potentially play a role in the scoliosis deformation process.
In the growth modulation equation, G m represents the growth rate (result of growth profile), β is the sensitivity factor and may be linked with the biological influences, and σ reflects mechanical factors (asymmetrical stress distribution). This is the first study that isolated and quantified the impact of growth profile (G m ) on scoliotic progression, showing that the augmented Gm, combined with an initial coronal deformity, will lead to progression of coronal deformity as well as decreasing the kyphosis angle in the spine.
The magnitude of the adopted parameters used in the analyses may have influenced the simulations results and therefore the conclusion. More specifically, although the sensitivity factor β was held constant between the explored cases, β may vary with the patient age . Moreover, AIS patients' progression is mostly concentrated immediately prior to puberty, which may be related to circulatory hormones  (i.e., estrogen and melatonin). Such notions feed the speculation that, due to a variation in biological factors, the sensitivity factor β may be influenced by the disturbed growth plate mechanotransduction. For this reason, a sensitivity analysis exploring the influence of this variable was performed. The outcome of such study was encouraging and demonstrated that the identified association between scoliotic progression and increased growth velocity was robust. Although the current modeling settings were determined according to existing literature, one must recognize the possibility of inter-patient variability. Nevertheless, the developed finite element platform effectively allowed for the isolation of growth rate influence to implicitly explore its impact on the progressive profiles of scoliotic patients.
This study presented the biomechanical comparison of accelerated AIS and normal vertebral body growth patterns on scoliotic progression using finite element modeling. Result of this analysis suggests that amplified AIS growth velocity could indeed lead to the supplementary progression of scoliosis and thus pose as a progressive risk factor. Whether the documentation of patient growth profiles has any clinical predictive value for curve progression deserves further investigations.
adolescent idiopathic scoliosis
finite element model
Funded by the Natural Sciences and Engineering Research Council of Canada and the Canada Research Chair Program.
- Lowe TG, Edgar M, Margulies JY, Miller NH, Raso VJ, Reinker KA, Rivard CH: Etiology of Idiopathic Scoliosis: Current Trends in Research. J Bone Joint Surg Am. 2000, 82-A (8): 1157-1168.PubMed
- Goldberg CJ, Dowling FE, Fogarty EE: Adolescent idiopathic scoliosis: is rising growth rate the triggering factor in progression?. European Spine Journal. 1993, 2 (1): 29-36. 10.1007/BF00301052.PubMedView Article
- Sanders JO, Browne RH, Cooney TE, Finegold DN, McConnell SJ, Margraf SA: Correlates of the Peak Height Velocity in Girls With Idiopathic Scoliosis. Spine. 2006, 31 (20): 2289-2295. 10.1097/01.brs.0000236844.41595.26.PubMedView Article
- Bjure J, Grimby G, Nachemson A: Correction of body height in predicting spirometric values in scoliotic patients. Scand J Clin Lab Invest. 1968, 21: 191-192.PubMedView Article
- Hefti FL, McMaster MJ: The effect of the adolescent growth spurt on early posterior spinal fusion in infantile and juvenile idiopathic scoliosis. J Bone Joint Surg Br. 1983, 65: 247-254.PubMed
- Little DG, Song KM, Katz D, Herring JA: Relationship of peak height velocity to other maturity indicators in idiopathic scoliosis in girls. J Bone Joint SurgAm. 2000, 82: 685-693.
- Song KM, Little DG: Peak height velocity as a maturity indicator for males with idiopathic scoliosis. J Pediatr Orthop. 2000, 20: 286-288.
- Ylikoski M: Spinal growth and progression of adolescent idiopathic scoliosis. European Spine Journal. 1993, 1 (4): 236-239. 10.1007/BF00298366.PubMedView Article
- Nissinen M, Heliovaara M, Seitsamo J, Poussa M: Trunk asymmetry, posture, growth, and risk of scoliosis. A three-year followup of Finnish prepubertal school children. Spine. 1993, 18: 8-13. 10.1097/00007632-199301000-00002.PubMedView Article
- Skogland LB, Miller JA: The length and proportions of the thoracolumbar spine in children with idiopathic scoliosis. Acta Orthop. Scandinavica. 1981, 52: 177-185. 10.3109/17453678108991778.View Article
- Guo X, Chau WW, Chan YL, Cheng JCY: Relative anterior spinal overgrowth in adolescent idiopathic scoliosis. The Journal of Bone and Joint Surgery. 2003, 85-B: 1026-1031.View Article
- Chu WC, Lam WW, Chan YL, Ng BK, Lam TP, Lee KM, Guo X, Cheng JC: Relative shortening and functional tethering of spinal cord in adolescent idiopathic scoliosis?: study with multiplanar reformat magnetic resonance imaging and somatosensory evoked potential. Spine. 2006, 31: E19-25. 10.1097/01.brs.0000193892.20764.51.PubMedView Article
- Chu WC, Man GC, Lam WW, Yeung BH, Chau WW, Ng BK, Lam TP, Lee KM, Cheng JC: Morphological and functional electrophysiological evidence of relative spinal cord tethering in adolescent idiopathic scoliosis. Spine. 2008, 33: 673-680. 10.1097/BRS.0b013e318166aa58.PubMedView Article
- Hägglund G, Karlberg J, Willner S: Growth in girls with adolescent idiopathic scoliosis. Spine. 1992, 17 (1): 108-111. 10.1097/00007632-199201000-00016.PubMedView Article
- Stokes IA: Analysis and simulation of progressive adolescent scoliosis by biomechanical growth modulation. Eur Spine J. 2007, 16 (10): 1621-1628. 10.1007/s00586-007-0442-7.PubMedPubMed CentralView Article
- Stokes IA, Spence H, Aronsson DD, Kilmer N: Mechanical modulation of vertebral body growth: implications for scoliosis progression. Spine. 1996, 21 (10): 1162-1167. 10.1097/00007632-199605150-00007.PubMedView Article
- Azegami H, Murachi S, Kitoh J, Ishida Y, Kawakami N, Makino M: Etiology of idiopathic scoliosis: Computational study. Clin Orthop Relat Res. Dec. 1998, 229-236. 357
- Plaats A, Veldhuizen AG, Verkerke GJ: Numerical simulation of asymmetrically altered growth as initiation mechanism of scoliosis. Ann Biomed Eng. 2007, 35 (7): 1206-1215. 10.1007/s10439-007-9256-3.PubMedPubMed CentralView Article
- Carrier J, Aubin CE, Villemure I, Labelle H: Biomechanical modelling of growth modulation following rib shortening or lengthening in adolescent idiopathic scoliosis. Med Biol Eng Comput. 2004, 42: 541-548. 10.1007/BF02350997.PubMedView Article
- Huynh AM, Aubin CE, Rajwani T, Bagnall KM, Villemure I: Pedicle growth asymmetry as a cause of adolescent idiopathic scoliosis: a biomechanical study. Eur Spine J. 2007, 16: 523-529. 10.1007/s00586-006-0235-4.PubMedPubMed CentralView Article
- Villemure I, Aubin CE, Dansereau J, Labelle H: Simulation of progressive deformities in adolescent idiopathic scoliosis using a biomechanical model integrating vertebral growth modulation. J Biomech Eng. 2002, 124: 784-790. 10.1115/1.1516198.PubMedView Article
- Villemure I, Aubin CE, Dansereau J, Labelle H: Biomechanical simulations of the spine deformation process in adolescent idiopathic scoliosis from different pathogenesis hypotheses. Eur Spine J. 2004, 13: 83-90. 10.1007/s00586-003-0565-4.PubMedPubMed CentralView Article
- Driscoll M, Aubin CE, Moreau A, Villemure I, Parent S: The role of spinal concave-convex biases in the progression of idiopathic scoliosis. Eur Spine J. 2009, 18 (2): 180-187. 10.1007/s00586-008-0862-z.PubMedPubMed CentralView Article
- Delorme S, Petit Y, de Guise JA, Labelle H, Aubin CE, Dansereau J: Assessment of the 3D reconstruction and high-resolution geometrical modeling of the human skeletal trunk from 2D radiographic images. IEEE Trans Biomed Eng. 2003, 50: 989-998. 10.1109/TBME.2003.814525.PubMedView Article
- Sylvestre PL, Villemure I, Aubin CE: Finite element modeling of the growth plate in a detailed spine model. Med Biol Eng Comput. 2007, 45: 977-988. 10.1007/s11517-007-0220-z.PubMedView Article
- Price J, Oyajobi B, Russell R: The cell biology of bone growth. Eur J Clin Nutr. 1994, 48: 131-149.
- Schultz A, Andersson GB, Ortengren R, Bjork R, Nordin M: Analysis and quantitative myoelectric measurements of loads on the lumbar spine when holding weights in standing postures. Spine. 1982, 7: 390-397. 10.1097/00007632-198207000-00009.PubMedView Article
- Patwardhan A, Havey R: A follower load increases the load-carrying capacity of the lumbar spine in compression. Spine. 1999, 24 (10): 1003-1009. 10.1097/00007632-199905150-00014.PubMedView Article
- Shirazi A, Parpianopour M: Load bearing and stress analysis of the human spine under a novel wrapping compression loading. 2000, 15: 718-725.
- Dime'glio A, Bonnel F: Le rachis en croissance scoliose, taille assise et puberté. 1990, Springer-Verlag, Paris, New York
- Stokes IA: Mechanical effects on skeletal growth. J Musculoskel Neuron Interact. 2002, 2: 277-280.
- Stokes I, Aronsson D, Dimock A, Cortright V, Beck S: Endochondral growth in growth plates of three species at two anatomical locations modulated by mechanical compression and tension. J Orthop Res. 2006, 24 (6): 1327-1333. 10.1002/jor.20189.PubMedPubMed CentralView Article
- Stokes IAF, Windisch L: Vertebral Height Growth Predominates Over Intervertebral Disc Height Growth in Adolescents With Scoliosis. Spine. 2006, 31 (14): 1600-1604. 10.1097/01.brs.0000222008.15750.1f.PubMedPubMed CentralView Article
- Stokes I, Aronsson D, Dimock A, Cortright , Beck S: Endochondral growth in growth plates of three species at two anatomical locations modulated by mechanical compression and tension. J of Orthopaedic Research. 2006, 10: 1327-33.View Article
- Wilke H, Neef P, Caimi M, Hoogland T, Claes L: New In vivo measurements pressures intervertebral disc in daily life. Spine. 1999, 24 (8): 755-762. 10.1097/00007632-199904150-00005.PubMedView Article
- Stokes IA: Mechanical effects on skeletal growth. J Musculoskel Neuro Interact. 2002, 2: 277-280.
- Escalada F, Marco E, Duarte E, Ma Muniesa J, Boza R, Tejero M, Cáceres E: Assessment of angle velocity in girls with adolescent idiopathic scoliosis. Scoliosis. 2009, 4: 20-10.1186/1748-7161-4-20.PubMedPubMed CentralView Article
- Stokes IA, Gookin DM, Reid S, Hazard RG: Effects of axis placement on measurement of isokinetic flexion and extension torque in the lumbar spine. J Spinal Disord. 1990, 3 (2): 114-118.PubMed
- Serrat MA, Lovejoy CO, King D: Age- and site-specific decline in insulin-like growth factor-I receptor expression is correlated with differential growth plate activity in the mouse hindlimb. Anat Rec (Hoboken). 2007, 290 (4): 375-381.View Article
- Leboeuf D, Letellier K, Alos N, Edery P, Moldovan F: Do estrogens impact adolescent idiopathic scoliosis?. Trends Endocrinol Metab. 2009, 20 (4): 147-152. 10.1016/j.tem.2008.12.004.PubMedView Article
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.