Front. Bioeng. Biotechnol. Frontiers in Bioengineering and Biotechnology Front. Bioeng. Biotechnol. 2296-4185 Frontiers Media S.A. 10.3389/fbioe.2020.00159 Bioengineering and Biotechnology Original Research Spinal Compressive Forces in Adolescent Idiopathic Scoliosis With and Without Carrying Loads: A Musculoskeletal Modeling Study Schmid Stefan 1 2 3 * Burkhart Katelyn A. 1 2 Allaire Brett T. 1 Grindle Daniel 1 4 Bassani Tito 5 Galbusera Fabio 5 Anderson Dennis E. 1 2 1Center for Advanced Orthopaedic Studies, Beth Israel Deaconess Medical Center, Boston, MA, United States 2Department of Orthopaedic Surgery, Harvard Medical School, Boston, MA, United States 3Spinal Movement Biomechanics Group, Division of Physiotherapy, Department of Health Professions, Bern University of Applied Sciences, Bern, Switzerland 4Division of Engineering Mechanics, Department of Biomedical Engineering and Mechanics, Virginia Polytechnic Institute and State University, Blacksburg, VA, United States 5Laboratory of Biological Structures Mechanics (LABS), IRCCS Istituto Ortopedico Galeazzi, Milan, Italy

Edited by: Alexandros E. Tsouknidas, University of Western Macedonia, Greece

Reviewed by: Rizwan Arshad, Royal Military College of Canada, Canada; Mohammad Nikkhoo, Islamic Azad University, Iran

 *Correspondence: Stefan Schmid, stefanschmid79@gmail.com

This article was submitted to Biomechanics, a section of the journal Frontiers in Bioengineering and Biotechnology

 03 03 2020 2020 8 159 16 12 2019 17 02 2020 Copyright © 2020 Schmid, Burkhart, Allaire, Grindle, Bassani, Galbusera and Anderson. 2020 Schmid, Burkhart, Allaire, Grindle, Bassani, Galbusera and Anderson

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

The pathomechanisms of curve progression in adolescent idiopathic scoliosis (AIS) remain poorly understood and biomechanical data are limited. A deeper insight into spinal loading could provide valuable information toward the improvement of current treatment strategies. This work therefore aimed at using subject-specific musculoskeletal full-body models of patients with AIS to predict segmental compressive forces around the curve apex and to investigate how these forces are affected by simulated load carrying. Models were created based on spatially calibrated biplanar radiographic images from 24 patients with mild to moderate AIS and validated by comparing predictions of paravertebral muscle activity with reported values from in vivo studies. Spinal compressive forces were predicted during unloaded upright standing as well as standing with external loads of 10, 15, and 20% of body weight (BW) applied to the scapulae to simulate carrying a backpack in the regular way on the back as well as in front of the body and over the shoulder on the concave and convex sides of the scoliotic curve. The predicted muscle activities around the curve apex were higher on the convex side for the erector spinae (ES) and multifidi (MF) muscles, which was comparable to the EMG-based in vivo measurements from the literature. In terms of spinal loading, the implementation of spinal deformity resulted in a 10% increase of compressive force at the curve apex during unloaded upright standing. Apical compressive forces further increased by 50–62% for a simulated 10% BW load and by 77–94% and 103–128% for 15% and 20% BW loads, respectively. Moreover, load-dependent compressive force increases were the lowest in the regular backpack and the highest in the frontpack and convex conditions, with concave side-carrying forces in between. The predictions indicated increased segmental compressive forces during unloaded upright standing, which could be ascribed to the scoliotic deformation. When carrying loads, compressive forces further increased depending on the carrying mode and the weight of the load. These results can be used as a basis for further studies investigating segmental loading in AIS patients during functional activities. Models can thereby be created using the same approach as proposed in this study.

 OpenSim inverse dynamics validation spinal loading AIS muscle asymmetry subject-specific modeling backpack carrying Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung10.13039/501100001711 Introduction

Adolescent idiopathic scoliosis (AIS) is a complex three-dimensional deformity of the spine, which affects as many as 4 out of 100 adolescents and occurs in early puberty (Cheng et al., 2015). Among congenital, neuromuscular, and other types of scoliosis, AIS is by far the most common form and is characterized by a poorly understood etiology and pathogenesis (de Seze and Cugy, 2012; Schlosser et al., 2014). Once the diagnosis of AIS is made, adolescents are generally treated conservatively using scoliosis-specific exercises (SSEs) and other forms of physiotherapy in order to minimize curve progression (Romano et al., 2012). For 10% of the initially diagnosed AIS patients, however, scoliosis exceeds a lateral deformation of 20–30° and brace treatment is indicated (Negrini et al., 2015). And finally, in about one fourth of the patients treated with exercises and braces, the progression of deformation cannot be contained and surgical intervention is required (Dolan and Weinstein, 2007).

In order to avoid bracing, surgery, and any associated health problems (e.g., chronic back pain), stopping curve progression in an early stage by means of SSEs is highly important. However, the effects of SSEs are not evident and further research is needed to clearly define the best types of SSEs as well as the frequency and intensity with which they should be administered (Romano et al., 2012). Since the pathomechanics of the AIS spine are not well understood, SSE protocols that affect spinal loading in a targeted, and scientifically sound manner cannot be developed. The literature currently lacks any studies reporting on spinal loading during functional activities or exercises in AIS patients.

Furthermore, it has been suggested that schoolbag carrying might play a role in the progression of scoliotic deformity and contribute to the development of back pain in AIS patients (Chow et al., 2006; Sahli et al., 2013). Sahli et al. (2013) recommended limiting backpack loads to 10% of body weight (BW) and carrying the load equally over both shoulders or over the shoulder on the concave rather than the convex side. However, these statements should be considered with caution since none of these studies investigated the effects of load carrying on spinal loading or trunk muscle forces.

Due to recent advancements in radiography-based geometric 3D reconstruction (Bassani et al., 2017) and musculoskeletal modeling (Bruno et al., 2015; Schmid et al., 2019), such parameters can be studies non-invasively and do not require invasive procedures such as intradiscal pressure or implant-based vertebral load measurements. The aim of this study was twofold: (1) To create subject-specific musculoskeletal full-body models of patients with mild to moderate AIS and validate predicted muscle activities with EMG data available in literature, and (2) to predict segmental compressive forces around the curve apex and how forces are affected by load carrying conditions.

 Materials and Methods Development of Subject-Specific Models Base Models

The base models for this study were created using our previously validated OpenSim-based musculoskeletal full-body models for children and adolescents aged 6–18 years (Schmid et al., 2019). They include a fully articulated thoracolumbar spine with a rib cage and are age- and gender-adjusted for sagittal spinal alignment as well as segmental inertial properties and maximum trunk muscle force capacity.

We enhanced the models with non-linear stiffness properties for flexion-extension and lateral-bending motions in all segments between T1/2 and L5/S1 using values from a recent meta-regression analysis over 45 studies involving experiments on adult cadaveric spines (Zhang et al., 2020). The properties were implemented using standard built-in linear bushing elements (expression-based bushing forces), which create reaction moments based on the rotational displacements of the adjacent vertebrae (Meng et al., 2015; Senteler et al., 2016). To avoid stiffness-related reaction moments in the neutral position of the spinal segments (i.e., position of the spinal segments when standing in an upright position), the bushing frames were oriented accordingly. Passive moments for segmental axial rotation were modeled using reserve coordinate actuators.

 Implementing Spinal Deformity

Using our enhanced base models, we created subject-specific models for 24 patients with mild to moderate AIS (Table 1). Spinal deformity was thereby implemented based on existing simultaneously captured and spatially calibrated anterior–posterior and lateral radiographic images (EOS Imaging, France) that were acquired within a previous study conducted at the IRCCS Istituto Ortopedico Galeazzi in Milan, Italy (Bassani et al., 2017). The protocol for this study was approved by the local ethics commission and patient assent and parental permission to use the anonymized radiological data were given by signing a written informed consent form.

Demographics of the patients with adolescent idiopathic scoliosis (AIS) from which the biplanar radiographic images were used to create the musculoskeletal models for the current study.

Patient Sex Age (years) Height (cm) Mass (kg) Curve type1 Cobb (°) Convexity
P1 Female 14 157 47 4B 23.9 Left
P2 Female 13 168 46 5 18.1 Left
P3 Female 12 140 35 1B 14.0 Right
P4 Female 14 165 50 3B 24.0 Right
P5 Female 15 168 48 5 24.5 Left
P6 Female 16 172 54 1A 25.5 Right
P7 Female 14 160 50 5 23.4 Right
P8 Male 11 134 39 5 26.5 Left
P9 Female 12 160 45 5 22.1 Left
P10 Female 15 165 45 5 22.4 Right
P11 Male 14 180 50 5 16.3 Right
P12 Female 13 155 44 5 18.3 Left
P13 Male 15 175 52 1A 28.3 Left
P14 Female 15 170 60 5 19.2 Left
P15 Female 15 157 48 5 27.7 Right
P16 Female 9 133 29 5 18.8 Right
P17 Female 11 144 35 5 18.3 Right
P18 Female 14 165 42 5 27.7 Left
P19 Female 13 167 60 5 21.1 Left
P20 Female 14 167 48 5 27.9 Right
P21 Female 16 169 63 1A 25.0 Right
P22 Female 15 168 51 1A 13.6 Right
P23 Female 17 170 53 5 24.7 Right
P24 Female 14 160 45 5 18.3 Left
Average (SD) 13.8 (1.8) 161.2 (12.3) 47.5 (8.0) – 22.1 (4.4) –
1Curve type classification according to Lenke et al. (2001).

Three-dimensional position and orientation of each vertebra from T1 to L5 was extracted from the radiographs using a custom MATLAB script (Bassani et al., 2017). In brief, this script contains a graphical user interface (GUI), which allows for the manual identification of nine characteristic landmarks per vertebra, i.e., upper and lower vertebral corners in the sagittal and frontal planes as well as the location of spinous process in the frontal plane. Based on these landmarks, sagittal and frontal vertebral orientations were calculated as the average of the slopes of the lines connecting upper and lower vertebral corners. The axial vertebral orientation was obtained by transforming the landmark which identifies the spinous process in 3D space through evaluating a fitted referential anatomic mesh model for the vertebra under assessment (Bassani et al., 2017). To account for individual spine height, we determined the height of each vertebral body using the geometric centers of the proximal and distal intervertebral disc spaces in the sagittal plane (centroids of the lower and upper corners of the proximal and distal vertebrae, respectively) and the vertebral tilt angle in the frontal plane (Figure 1).

Determination of vertebral body height (HTotal) from radiographic images in two steps: (1) Extracting the distance between the centroids of the upper and lower intervertebral disc spaces in the sagittal plane (HSagittal), and (2) correct it for scoliosis-dependent projection errors using the vertebral tilt angle measured in the frontal plane.

Subject-specific models were created in four steps: (1) scaling the base model from the corresponding age- and gender-group by body height and body mass, (2) implementing 3D spinal deformity by adjusting the orientation of the vertebral bodies T1 to L5 in the flexion/extension, lateral-flexion and axial rotation directions, (3) scaling intersegmental joint distances by vertebral body height, and (4) re-adjusting 3D orientation of the lumped head and neck segment as well as arms and ribs to a neutral position (Figure 2).

Model creation in four steps: (1) scaling of base model by body height and body mass, (2) implementing 3D spinal deformity by adjusting the orientation of the vertebral bodies T1 to L5, (3) scaling intersegmental joint distances by vertebral body height, and (4) re-adjusting 3D orientation of the lumped head and neck segment as well as arms and ribs to a neutral position.

 Evaluation of Muscle Geometry

It was previously reported that AIS is associated with side-to-side asymmetries in erector spinae (ES) and multifidi (MF) muscle geometry (Zoabli et al., 2007; Zapata et al., 2015). To ensure appropriate handling of the muscle geometry in our models, we therefore estimated bilateral CSAs of the modeled ES and MF muscles for each thoracic and lumbar vertebral mid-plane by summing the CSAs of the individual fascicles [calculated by dividing the maximum force generating capacity of the respective fascicle by an assumed uniform maximal muscle stress (MMS) of 100 N/cm2] crossing the respective mid-plane to compute an equivalent muscle group CSA at that level (Bruno et al., 2015). We refer to this procedure as a “virtual CT scan” of the model, which enables the comparison of model muscle geometry with conventional medical imaging studies. For the comparison with the literature, an asymmetry ratio was calculated by dividing the CSA of the muscles on the convex side by the CSA of the muscles on the concave side. In accordance with the in vivo studies, this ratio was calculated for ES at the curve apex and for MF at the levels T8, L1, and L4 as well as at the curve apex.

 Simulations

All simulations were carried out using OpenSim 3.3 (Delp et al., 2007) and MATLAB R2019a (MathWorks Inc., Natick, MA, United States). Models were solved using an inverse dynamics based static optimization with a cost function that minimized the sum of squared muscle activation (Herzog, 1987). Due to uncertainties of how the scoliotic deformation might affect muscle physiology, we solved the models without considering force-length relationships.

Model Validation

To comply with the best practice guidelines for verification and validation of musculoskeletal models (Hicks et al., 2015), we aimed at providing a reasonable validation for our models by comparing model predictions to in vivo values available from the literature. Unfortunately, the literature lacks in vivo studies reporting on functional spinal loading (e.g., segmental compressive forces or intradiscal pressure in upright standing conditions) in AIS patients. For this reason, we conducted simulations of trunk muscle activity in standing and prone positions to compare them with the results of three reported in vivo studies using surface electromyography (EMG) (Cheung et al., 2005; Kwok et al., 2015; Stetkarova et al., 2016). For each of the simulations, we selected the AIS models that matched best the respective in vivo study population in terms of curve location.

To evaluate the accuracy of our models to predict ES muscle activity, we placed the models in a neutral upright standing position and compared the convex to concave ratios of the average activation levels of the ES muscle fascicles in the lumbar and thoracic regions as well as at the curve apex and upper and lower curve limits (two levels above and below the apex) to surface EMG-based in vivo measurements in AIS patients with main thoracic and thoracolumbar curves (Kwok et al., 2015) as well as patients with non-progressive AIS (Cheung et al., 2005), respectively. Accuracy of MF muscle activity predictions was evaluated by placing the models in a prone position and comparing the convex to concave ratio of the average activation levels of the MF muscle fascicles at the curve apex to needle EMG-based in vivo measurements in AIS patients with main thoracic curves (Stetkarova et al., 2016).

Muscle fascicles were thereby selected based on the surface electrode placement and needle electrode insertion locations described in the respective in vivo studies. To simulate the floor and the table where the models were “standing or lying on,” we used residual point actuators with maximum activation at a force of 10 kN, which was shown to be large enough to provide the required support with minimal expenses in the static optimization (Schmid et al., 2019). Model predictions and in vivo measurements were compared qualitatively.

 Prediction of Compressive Forces

In order to investigate spinal compressive forces in upright standing AIS patients with and without carrying loads, we conducted simulations in five different conditions (Table 2).

Conditions for simulating spinal compressive forces during upright standing with and without carrying loads in patients with adolescent idiopathic scoliosis (AIS).

Condition Models External load
Unloaded AIS and undeformed1 No external load applied
Backpack Only AIS External loads of 10, 15, and 20% of body weight (BW) applied 20° dorsally angled to the lateral third of the upper edge of the scapulae (equally distributed between sides) to simulate a regular backpack
Frontpack Only AIS External loads of 10, 15, and 20% of BW applied 25° ventrally angled to the lateral third of the upper edge of the scapulae (equally distributed between sides) to simulate a backpack carried in front of the body
Sidepack concave Only AIS External loads of 10, 15, and 20% of BW applied 5° dorsally and 10° laterally angled to the lateral third of the upper edge of the scapula on the concave side to simulate a backpack carried unilaterally over the shoulder on the concave side
Sidepack convex Only AIS External loads of 10, 15, and 20% of BW applied 5° dorsally and 10° laterally angled to the lateral third of the upper edge of the scapula on the convex side to simulate a backpack carried unilaterally over the shoulder on the convex side
The angles for the external loads were derived from standardized photographs of one of the investigators carrying a backpack in the different modes. 1Models before implementation of spinal deformity.

After solving the models, joint reaction analysis was carried out to calculate the axial compressive forces acting on the spinal segments at the curve apex as well as one and two levels above and below. In addition to the absolute force magnitudes, compressive forces in AIS patients in the unloaded condition (1) were expressed as a percentage of the forces derived from the undeformed models, whereas compressive forces in AIS patients in the loaded conditions (2–5) were expressed as a percentage of the unloaded condition. To provide a coherent overview of these percentages, data were presented using violin plots with superimposed boxplots and individual values. Wilcoxon signed rank tests with an alpha-level set to 5% were conducted to test for differences from 100%.

 Results Muscle Geometry

The predicted mean CSA ratio for the ES muscle at the apex was 1.05 SD 0.08, which compares to the ratio of 1.02 calculated from MRI-derived ES muscle volumes reported in the literature (Zoabli et al., 2007; Figure 3). For the MF muscle, predicted mean CSA ratios for the levels T8, L1, and L4 were 0.95 SD 0.21, 0.92 SD 0.09, and 1.01 SD 0.03, respectively. This corresponded well to the values reported in the literature for T8 (0.96), L1 (0.95), and L4 (0.98) (Zapata et al., 2015). Predicted mean CSA ratio for the MF muscle at the curve apex was 0.89 SD 0.11, however, no in vivo values were available for comparison.

Convex to concave ratios for cross-sectional areas (CSA) of the erector spinae (ES) muscle at the curve apex (left) as well as the multifidi (MF) muscles at the levels T8, L1, L4, and the curve apex (right).

 Validation Studies

The predicted mean ES muscle activity ratios during upright standing indicated higher activity on the convex side of the muscle for thoracic (thoracic region: 1.29 SD 0.71; lumbar region: 1.18 SD 0.78) and thoracolumbar curves (thoracic region: 1.05 SD 0.19; lumbar region: 1.12 SD 0.31) (Figure 4). Furthermore, ES muscle activity within the curve also indicated higher convex activity at the apex (1.72 SD 1.06) as well as the upper (1.37 SD 0.78) and lower curve ends (1.02 SD 0.83). These ratios compare reasonably well to the literature for the thoracic and lumbar portions of the ES muscle in patients with main thoracic curves (1.87 SD 1.7 and 1.7 SD 0.85, respectively) as well as at the curve apex (2.1 SD 1.38) and the lower curve end (0.96 SD 0.32) (Cheung et al., 2005; Kwok et al., 2015). For MF muscle activity, models predictions indicated higher activity at the apex on the convex side of the curve (1.21 SD 0.50), which goes along with the in vivo measured ratio of 1.38 (Stetkarova et al., 2016).

Convex to concave ratios for erector spinae (ES) muscle activity in the thoracic and lumbar regions (left) as well as at the apex and the upper and lower ends of the curve (middle) during upright standing, and for multifidi (MF) muscle activity at the curve apex (right) in prone position.

 Spinal Compressive Forces

The implementation of spinal deformity resulted in higher median compressive forces within the scoliotic curve, with forces constantly increasing from two levels above the apex (103% IQR 13%, p = 0.092), to the apex (110% IQR 13%, p < 0.001), to two levels below the apex (118% IQR 16%, p < 0.001) (Figure 5).

Spinal compressive forces at the curve apex as well as two levels above and below in patients with adolescent idiopathic scoliosis (AIS) expressed as percentages of the forces in the same individuals but without spinal deformity (% of undeformed). Data are presented in form of violin plots with superimposed boxplots and individual values.

All loaded median axial compressive forces were significantly different from 100% (i.e., from an unloaded condition) at a level of p < 0.001. When carrying a load which corresponded to 10% of BW, median axial compressive forces increased on average by 50% for the regular backpack, 62% for the backpack carried in front and 60% and 54% for the backpack carried over the shoulder on the convex and concave sides, respectively (Figure 6). In the regular backpack as well as both sidepack conditions, compressive force increased the most above the apex and the least below the apex, whereas in the frontpack condition, compressive force increased about equally on all spinal levels.

Spinal compressive forces at the curve apex as well as two levels above and below in patients with adolescent idiopathic scoliosis (AIS) when carrying a backpack with a weight corresponding to 10% of body weight in the regular way (top left), in front of the body (top right) as well as over the shoulder on the concave (bottom left), and convex sides (bottom right) of the scoliotic curve. Forces are expressed as percentages of unloaded upright standing and presented in form of violin plots with superimposed boxplots and individual values.

When applying loads corresponding to 15 and 20% of BW, median axial compressive forces increased on average by 77 and 103% for the backpack, 93 and 125% for the frontpack, 94 and 128% for the convex sidepack, and 85 and 116% for the concave sidepack conditions, respectively (Figures 7, 8). The load distribution pattern within the spinal levels of the scoliotic curve remained similar as described for the 10% of BW load.

Spinal compressive forces at the curve apex as well as two levels above and below in patients with adolescent idiopathic scoliosis (AIS) when carrying a backpack with a weight corresponding to 15% of body weight in the regular way (top left), in front of the body (top right) as well as over the shoulder on the concave (bottom left), and convex sides (bottom right) of the scoliotic curve. Forces are expressed as percentages of unloaded upright standing and presented in form of violin plots with superimposed boxplots and individual values.

Spinal compressive forces at the curve apex as well as two levels above and below in patients with adolescent idiopathic scoliosis (AIS) when carrying a backpack with a weight corresponding to 20% of body weight in the regular way (top left), in front of the body (top right) as well as over the shoulder on the concave (bottom left), and convex sides (bottom right) of the scoliotic curve. Forces are expressed as percentages of unloaded upright standing and presented in form of violin plots with superimposed boxplots and individual values.

A complete set of the absolute force magnitude and relative force percentage values can be found in the Supplementary Material.

 Discussion

We created 24 subject-specific musculoskeletal full-body models from biplanar radiographic images of patients with mild to moderate AIS and validated these models by comparing predictions of paravertebral muscle activity with reported values from in vivo studies. Moreover, we predicted apical spinal loads with and without simulated load carrying, i.e., carrying a backpack in the regular way, carrying a backpack in front of the body and carrying a backpack over the shoulder on the concave and convex sides of the scoliotic curve.

The evaluation of muscle geometry indicated that the implementation of spinal deformity resulted in side-to-side asymmetries, which agreed with reports in the literature. The validation studies showed higher convex ES and MF muscle activity around the curve apex, which was comparable to the EMG-based in vivo measurements from the literature. Measurements of overall thoracic and lumbar ES muscle activity agreed well for thoracic but less thoracolumbar curves. In terms of spinal loading, the implementation of spinal deformity resulted in a 10% increase of compressive force at the curve apex during unloaded upright standing. Apical compressive forces further increased by 50–62% for a simulated 10% BW load and by 77–94% and 103–128% for 15 and 20% BW loads, respectively. Moreover, load-dependent compressive force increases were the lowest in the regular backpack and the highest in the frontpack and convex conditions, with concave side-carrying forces in between.

Even though the evaluation of muscle geometry estimated from our models indicated larger ES muscle CSA on the convex side of the scoliotic deformation, about one third of our models showed in fact larger ES muscle CSA on the concave side. This agrees with the study of Zoabli et al. (2007), reporting that even though there was an overall tendency for larger ES muscle volume on the convex side, some patients also presented larger volumes on the concave side. For the MF muscle at the curve apex, a larger CSA on the concave side was found in about 90% of the models, whereas the remaining 10% showed a larger CSA on the convex side. Despite these uncertainties, however, the muscle geometry estimated from our AIS models was comparable to what was measured in vivo. This raises the question of whether the previously reported muscle volume/thickness asymmetries in AIS patients represent an actual change in muscle, or is just a result of measurement of muscle with different spinal geometry.

The partial disagreements between our model predictions and the ES muscle activity ratios derived from the results reported by Kwok et al. (2015) might be related to differences in curve characteristics (i.e., location of the apex and severity of deformation) between the respective patient populations. Especially when evaluating overall ES muscle activity in the thoracic and lumbar regions with the same electrode placement for all patients, different curve characteristics could have a significant effect on muscle activation. Furthermore, muscle fiber redistribution with a higher proportion of type I fibers on the convex side of the AIS curve (Gonyea et al., 1985; Meier et al., 1997; Mannion et al., 1998; Stetkarova et al., 2016) might be another contributing factor, since muscle fiber type seems to have an influence on the EMG signal (Poosapadi Arjunan et al., 2016). However, due to the lack of appropriate data, the consideration of fiber distribution change in our current models would be associated with too many assumptions. Finally, it is not known whether AIS has an influence on the force-length-relationship of the paravertebral muscles, which could also have an influence on EMG activity. Based on our validation studies, however, we consider the muscle activation patterns predicted by our models comparable to the patterns reported in the in vivo studies.

This is the first study using inverse dynamics-based musculoskeletal full-body modeling to investigate the immediate effect of AIS-related spinal deformity on axial compressive forces within the scoliotic curve. The results suggested that spinal deformity causes an overall increase in compressive forces, with forces increasing the most at the lower and the least at the upper end of the curve. When considering the different curve types, it appears that the average increase in segmental loading was not significantly related to the spinal level of the curve apex or the curve severity, i.e., the Cobb angles (Pearson correlation: r = 0.28, p = 0.193 and r = 0.22, p = 0.302, respectively). It should be considered, however, that all the patients in this study had mild to moderate AIS, providing a relatively small range of Cobb angles. The relationship between compressive force increase and curve severity would probably be more pronounced with larger range of Cobb angles in the dataset. It should also be noted that a small fraction of the predictions resulted in compressive force decreases. When looking at these cases, especially those predicting compressive forces of less than 90% of the undeformed, it appears that these patients tended to have particularly flat thoracic sagittal profiles (i.e., <25° of thoracic kyphosis), which might have resulted in a significant reduction of muscular effort (specifically the ES muscle) and therefore in reduced compressive forces.

The simulation of load carrying in AIS patients indicated compressive force increments that were dependent on the carrying mode as well as the weight of the load. Carrying the load in front of the body resulted in considerably higher compressive forces compared to regular backpack carrying. This is not surprising since the front carrying mode would be assumed to cause higher muscular effort to prevent increased flexion. Interestingly, carrying the backpack on the concave side yielded compressive forces that were slightly higher than the ones for regular backpack carrying, but lower compared to carrying the load in front. This indicates that shifting the load from the back to the concave side does result in increased compressive forces – most likely due to the higher muscular effort on the contralateral side – but not as much as carrying the load symmetrically in front. In addition, carrying the load on the convex side caused compressive forces that were comparable to the front carrying mode, most likely because the load was acting more directly on the spinal segments. The different patterns of compressive force increments between the frontpack condition (constant over all levels) and the other carrying conditions (decrease in force increments from two levels above to two levels below the apex) might be related to the direction of the applied external load in the sagittal plane. In the frontpack condition, the external load was directed anteriorly, which resulted in increased paraspinal muscle activity. In the other carrying conditions, however, the external load is directed posteriorly, causing an increased activity predominantly of the abdominal muscles. Considering that most of the abdominal muscles are only indirectly connected to the spine, i.e., over the rib cage, it seems plausible that this affected spinal loading differently than in the frontpack condition. In any way, these results do not allow any conclusions on whether carrying a load in front, on either side or regularly on the back is advantageous to minimize curve progression or prevent complications such as joint degeneration or back pain. In fact, it is possible that carrying the load on the convex side might put the patients at a higher risk for back pain, but at the same time slow down curve progression by modulating the segmental load in a way that vertebral growth is positively affected due to the Hueter-Volkmann law, where bone growth is slowed with compression and accelerated with distraction (de Seze and Cugy, 2012). Future studies should therefore address this issue by using a combination of motion capture-driven musculoskeletal and finite element models.

This study has some important limitations that should be discussed. First of all, due to a lack of appropriate data for healthy children and adolescents as well as patients with AIS, the passive segmental stiffness properties that were implemented in the current base models were derived from healthy adult cadaveric spines. This issue could be addressed in the future by conducting clinical studies in AIS patients using approaches such as the intraoperative determination of load-displacement behavior proposed by Reutlinger et al. (2012). However, the lack of appropriate stiffness properties did not affect the current predictions since all simulations were conducted in a neutral position, i.e., with the spinal segments assumed to be in the neutral zone (Smit et al., 2011). Stiffness-related reaction moments would only have occurred with induced segmental rotations. When using the models for future investigations involving simulations beyond the neutral position of the spine, on the other hand, these limitations will have to be considered. Secondly, we did only consider compressive forces in this study, but forces in other directions (i.e., anterior–posterior and medial–lateral shear forces) might also be strongly affected by scoliotic deformities. Furthermore, the patient populations of the in vivo studies used for the validation of our models did not exactly match the population from which they were created. It is therefore advised that future studies investigating spinal loading in AIS include EMG measurements of the paraspinal muscles for more specific validations of the respective models. Lastly, the current simulations were not based on real-life kinematics, i.e., they were not driven by motion capture data. Especially for the load carrying investigations, it can be assumed that real subjects would have slightly adapted their posture based on the applied load, such as previously observed for regular backpack carrying in healthy young adults (Neuschwander et al., 2010).

In conclusion, this study used validated subject-specific OpenSim-based musculoskeletal full-body models to provide an insight into spinal loading in patients with AIS with and without carrying loads. The predictions indicated increased segmental compressive forces of about 10% around the curve apex during unloaded upright standing. When carrying loads, compressive forces further increased depending on the carrying mode and the weight of the load. These results can be used as a basis for further studies investigating segmental loading in AIS patients during functional activities. Models can thereby be created using the same approach as proposed in this study.

 Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

 Author Contributions

SS substantially contributed to the conception and design of the study, created the subject-specific models and performed the simulations, analyzed the data, interpreted the results, and wrote the manuscript. KB, BA, and DG provided substantial assistance to the creation of the models and the simulations. TB and FG provided the EOS images and MATLAB scripts for the extraction of the spinal deformity parameters. DA substantially contributed to the conception and design of the study, and provided substantial assistance to the creation of the models and the simulations. All authors critically revised the manuscript and approved the version to be published.

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Funding. This work was funded by the Swiss National Science Foundation (SNSF, grant no. 178427) and the National Center for Simulation in Rehabilitation Research (NCSRR, sub-award of NIH grant no. 5P2CHD065690).

The authors thank Dr. Bram Verhofste for assistance in the clinical interpretation of the radiographs.

 Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2020.00159/full#supplementary-material

 References Bassani T. Ottardi C. Costa F. Brayda-Bruno M. Wilke H. J. Galbusera F. (2017). Semiautomated 3D spine reconstruction from biplanar radiographic images: prediction of intervertebral loading in scoliotic subjects. Front. Bioeng. Biotechnol. 5:1. 10.3389/fbioe.2017.00001 28164082 Bruno A. G. Bouxsein M. L. Anderson D. E. (2015). Development and validation of a musculoskeletal model of the fully articulated thoracolumbar spine and rib cage. J. Biomech. Eng. 137:081003. 10.1115/1.4030408 25901907 Cheng J. C. Castelein R. M. Chu W. C. Danielsson A. J. Dobbs M. B. Grivas T. B. (2015). Adolescent idiopathic scoliosis. Nat. Rev. Dis. Primers 1:15030. 10.1038/nrdp.2015.30 27188385 Cheung J. Halbertsma J. P. Veldhuizen A. G. Sluiter W. J. Maurits N. M. Cool J. C. (2005). A preliminary study on electromyographic analysis of the paraspinal musculature in idiopathic scoliosis. Eur. Spine J. 14 130–137. 10.1007/s00586-004-0780-7 15368104 Chow D. H. Kwok M. L. Cheng J. C. Lao M. L. Holmes A. D. Au-Yang A. (2006). The effect of backpack weight on the standing posture and balance of schoolgirls with adolescent idiopathic scoliosis and normal controls. Gait. Posture 24 173–181. 10.1016/j.gaitpost.2005.08.007 16198112 de Seze M. Cugy E. (2012). Pathogenesis of idiopathic scoliosis: a review. Ann. Phys. Rehabil. Med. 55 128–138. 10.1016/j.rehab.2012.01.003 22321868 Delp S. L. Anderson F. C. Arnold A. S. Loan P. Habib A. John C. T. (2007). OpenSim: open-source software to create and analyze dynamic simulations of movement. IEEE Trans. Biomed. Eng. 54 1940–1950. 10.1109/tbme.2007.901024 18018689 Dolan L. A. Weinstein S. L. (2007). Surgical rates after observation and bracing for adolescent idiopathic scoliosis: an evidence-based review. Spine 32 S91–S100. 17728687 Gonyea W. J. Moore-Woodard C. Moseley B. Hollmann M. Wenger D. (1985). An evaluation of muscle pathology in idiopathic scoliosis. J. Pediatr. Orthop. 5 323–329. 10.1097/01241398-198505000-00013 3998135 Herzog W. (1987). Individual muscle force estimations using a non-linear optimal design. J. Neurosci. Methods 21 167–179. 10.1016/0165-0270(87)90114-2 3682873 Hicks J. L. Uchida T. K. Seth A. Rajagopal A. Delp S. L. (2015). Is my model good enough? Best practices for verification and validation of musculoskeletal models and simulations of movement. J. Biomech. Eng. 137:020905. 10.1115/1.4029304 25474098 Kwok G. Yip J. Cheung M. C. Yick K. L. (2015). Evaluation of myoelectric activity of paraspinal muscles in adolescents with idiopathic scoliosis during habitual standing and sitting. Biomed. Res. Int. 2015 958450. 10.1155/2015/958450 26583151 Lenke L. G. Betz R. R. Harms J. Bridwell K. H. Clements D. H. Lowe T. G. (2001). Adolescent idiopathic scoliosis: a new classification to determine extent of spinal arthrodesis. J. Bone Joint Surg. Am. 83 1169–1181. 10.2106/00004623-200108000-00006 11507125 Mannion A. F. Meier M. Grob D. Muntener M. (1998). Paraspinal muscle fibre type alterations associated with scoliosis: an old problem revisited with new evidence. Eur. Spine J. 7 289–293. 10.1007/s005860050077 9765036 Meier M. P. Klein M. P. Krebs D. Grob D. Muntener M. (1997). Fiber transformations in multifidus muscle of young patients with idiopathic scoliosis. Spine 22 2357–2364. 10.1097/00007632-199710150-00008 9355216 Meng X. Bruno A. G. Cheng B. Wang W. Bouxsein M. L. Anderson D. E. (2015). incorporating six degree-of-freedom intervertebral joint stiffness in a lumbar spine musculoskeletal model-method and performance in flexed postures. J. Biomech. Eng. 137:101008. 10.1115/1.4031417 26299207 Negrini S. Minozzi S. Bettany-Saltikov J. Chockalingam N. Grivas T. B. Kotwicki T. (2015). Braces for idiopathic scoliosis in adolescents. Cochrane Database Syst. Rev. 6 Cd006850. Neuschwander T. B. Cutrone J. Macias B. R. Cutrone S. Murthy G. Chambers H. (2010). The effect of backpacks on the lumbar spine in children: a standing magnetic resonance imaging study. Spine 35 83–88. 10.1097/brs.0b013e3181b21a5d 20023607 Poosapadi Arjunan S. Kumar D. K. Wheeler K. Shimada H. Siddiqi A. (2016). Effect of number of motor units and muscle fibre type on surface electromyogram. Med. Biol. Eng. Comput. 54 575–582. 10.1007/s11517-015-1344-1 26223565 Reutlinger C. Hasler C. Scheffler K. Buchler P. (2012). Intraoperative determination of the load-displacement behavior of scoliotic spinal motion segments: preliminary clinical results. Eur. Spine J. 21(Suppl. 6), S860–S867. 22315034 Romano M. Minozzi S. Bettany-Saltikov J. Zaina F. Chockalingam N. Kotwicki T. (2012). Exercises for adolescent idiopathic scoliosis. Cochrane Database Syst. Rev. 8 Cd007837. Sahli S. Rebai H. Ghroubi S. Yahia A. Guermazi M. Elleuch M. H. (2013). The effects of backpack load and carrying method on the balance of adolescent idiopathic scoliosis subjects. Spine J. 13 1835–1842. 10.1016/j.spinee.2013.06.023 24095102 Schlosser T. P. Van Der Heijden G. J. Versteeg A. L. Castelein R. M. (2014). How ‘idiopathic’ is adolescent idiopathic scoliosis? A systematic review on associated abnormalities. PLoS One 9:e97461. 10.1371/journal.pone.0097461 24820478 Schmid S. Burkhart K. A. Allaire B. T. Grindle D. Anderson D. E. (2019). Musculoskeletal full-body models including a detailed thoracolumbar spine for children and adolescents aged 6-18 years. J. Biomech. 10.1016/j.jbiomech.2019.07.049 31471110 Senteler M. Weisse B. Rothenfluh D. A. Snedeker J. G. (2016). Intervertebral reaction force prediction using an enhanced assembly of OpenSim models. Comput. Methods Biomech. Biomed. Engin. 19 538–548. 10.1080/10255842.2015.1043906 26031341 Smit T. H. Van Tunen M. S. Van Der Veen A. J. Kingma I. Van Dieen J. H. (2011). Quantifying intervertebral disc mechanics: a new definition of the neutral zone. BMC Musculoskelet. Disord. 12:38. 10.1186/1471-2474-12-38 21299900 Stetkarova I. Zamecnik J. Bocek V. Vasko P. Brabec K. Krbec M. (2016). Electrophysiological and histological changes of paraspinal muscles in adolescent idiopathic scoliosis. Eur Spine J. 25 3146–3153. 10.1007/s00586-016-4628-8 27246349 Zapata K. A. Wang-Price S. S. Sucato D. J. Dempsey-Robertson M. (2015). Ultrasonographic measurements of paraspinal muscle thickness in adolescent idiopathic scoliosis: a comparison and reliability study. Pediatr. Phys. Ther. 27 119–125. 10.1097/PEP.0000000000000131 25695194 Zhang C. Mannen E. M. Sis H. L. Cadel E. S. Wong B. M. Wang W. (2020). Moment-rotation behavior of intervertebral joints in flexion-extension, lateral bending, and axial rotation at all levels of the human spine: a structured review and meta-regression analysis. J. Biomech. 100 109579. 10.1016/j.jbiomech.2019.109579 31911050 Zoabli G. Mathieu P. A. Aubin C. E. (2007). Back muscles biometry in adolescent idiopathic scoliosis. Spine J. 7 338–344. 10.1016/j.spinee.2006.04.001 17482118