Xin Xi1,2†
Zhi Ling3†
Cong Wang3Chunya Gu4Xuqiang Zhan1,2Haixin Yu5Siqi Lu1,2
Tsung-Yuan Tsai3,6*
Yan Yu1,2*Liming Cheng1,2- 1Department of Spine Surgery, Tongji Hospital, School of Medicine, Tongji University, Shanghai, China
- 2Key Laboratory of Spine and Spinal Cord Injury Repair and Regeneration, Ministry of Education, Department of Spine Surgery, Tongji Hospital, School of Medicine, Tongji University, Shanghai, China
- 3School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China
- 4Department of Spinal Rehabilitation, Tongji Hospital, School of Medicine, Tongji University, Shanghai, China
- 5Department of Orthopedic Surgery, Yangpu Hospital, School of Medicine, Tongji University, Shanghai, China
- 6TAOiMAGE Medical Technologies Corporation, Shanghai, China
Traditional optical motion capture (OMC) with retroreflective markers is commonly used to measure joint kinematics but was also reported with unavoidable soft tissue artifacts (STAs) when quantifying the motion of the spine. Additionally, the patterns of the STA on the lumbar spine remain unclear. This study aimed to 1) quantify the in vivo STAs of the human lower back in three-dimensional directions during weight-bearing forward–backward bending and 2) determine the effects of the STAs on the calculated flexion angles between the upper and lower lumbar spines and adjacent vertebrae by comparing the skin marker (SM)- and virtual bone marker (VM)-based measurements. Six healthy volunteers were imaged using a biplanar radiographic system, and thirteen skin markers were mounted on every volunteer’s lower back while performing weight-bearing forward–backward bending. The STAs in the anterior/posterior (AP), medial/lateral (ML), and proximal/distal (PD) directions were investigated. The flexion angles between the upper and lower lumbar segments and adjacent intervertebral segments (L2–L5) throughout the cycle were calculated. For all the participants, STAs continuously increased in the AP direction and exhibited a reciprocal trend in the PD direction. During flexion, the STA at the lower lumbar region (L4–L5: 13.5 ± 6.5 mm) was significantly higher than that at the upper lumbar (L1–L3: 4.0 ± 1.5 mm) in the PD direction (p < 0.01). During extension, the lower lumbar (L4–L5: 2.7 ± 0.7 mm) exhibited significantly less STAs than that exhibited by the upper lumbar region (L1–L3: 6.1 ± 3.3 mm) (p < 0.05). The STA at the spinous process was significantly lower than that on both sides in the AP direction (p < 0.05). The present results on STAs, based on dual fluoroscopic measurements in healthy adult subjects, presented an anatomical direction, marker location, and anatomic segment dependency, which might help describe and quantify STAs for the lumbar spine kinematics and thus help develop location- and direction-specific weighting factors for use in global optimization algorithms aimed at minimizing the effects of STAs on the calculation of lumbar joint kinematics in the future.
1 Introduction
Low back pain (LBP) is one of the most common musculoskeletal disorders and a leading cause of disability, globally creating a substantial personal, community, and financial burden (Airaksinen et al., 2006; Delitto et al., 2012; Itz et al., 2013; Maher et al., 2017). Approximately 85%–95% of LBP cases have no identifiable cause or pathology and are, therefore, classified as non-specific LBP (NS-LBP). Hoy et al. (2012), Maher et al. (2017), and Wallbott (1989) characterized the way patients with NS-LBP perform daily activities and defined movement quality (MQ) as “the way in which human movement is executed with respect to the dimensions of time and space.” Observation, analysis, and the influence of whole-body movements on the MQ are the key elements of LBP management (Laird et al., 2012; Hodges et al., 2013).
The common final manifestation of LBP is a change in spine kinematics (Liguori and ACSM, 2021). Accurately quantifying three-dimensional (3D) joint kinematics during in vivo lumbar motion is essential for a detailed understanding of joint function, investigating pathologies, and assessing the success of therapies (Aiyangar et al., 2014). 3D motion analysis, such as optical motion capture (OMC), using skin markers is the most common method for measuring spine joint kinematics in vivo (Wang et al., 2016; Huang et al., 2021). Unfortunately, when skin-marker trajectories are used in human motion analysis, the skin-mounted markers move over the underlying bone, generating the so-called soft tissue artifact (STA), which makes the estimation of the instantaneous skeletal pose difficult (Leardini et al., 2005). Motions of soft tissue covering the spine may cause large errors (Mo¨rl and Blickhan, 2006). STA represents one of the main limitations in obtaining accurate and reliable skeletal kinematics from motion capture (Cereatti et al., 2017).
Many previous studies have demonstrated that the extent of STAs, caused by a combination of skin stretching and sliding, muscle contraction, gravity, and inertia, is unique to each specific body segment (Akbarshahi et al., 2010), the physical characteristics of individuals (Holden et al., 1997), marker locations (Schwartz et al., 2004), and the nature of the performed movement task (Fuller et al., 1997). Several studies have validated the STAs of skin markers on different body regions (Wang et al., 2016; Camomilla and Bonci, 2020; Metcalf et al., 2020). However, validation studies for lumbar spine segments are rare (Zemp et al., 2014, Mo¨rl and Blickhan, 2006). Moreover, most studies investigated lumbar STAs in a seated upright position (Mo¨rl and Blickhan, 2006) or static sitting position (Zemp et al., 2014), and a few have been validated against global spinal shape.
It is also important to quantify the propagation of STAs to estimate lumbar joint kinematics. Various studies have used intrusive techniques to quantify STAs; for example, intra-cortical pins (Fuller et al., 1997; Benoit et al., 2006; Dal Maso et al., 2016), external fixators (Cappozzo et al., 1996), and percutaneous skeletal trackers (Manal et al., 2000) have been used to quantify joint motion in vivo. Unfortunately, these devices restrict the movement of the participant and alter the normal, unimpeded sliding of the soft tissues relative to the underlying bone. Moreover, these techniques are invasive; therefore, their use is limited. To overcome these problems, non-invasive methods, such as magnetic resonance imaging (MRI) (Sangeux et al., 2006; Zemp et al., 2014) and X-ray fluoroscopy (Garling et al., 2007), have been used to quantify joint motion in vivo. However, the methods proposed in these studies have several limitations: 1) they capture several terminal static poses in different directions rather than measuring a continuous dynamic motion; 2) they use a skin marker set and compare the measurements collected from 2D images to quantify STAs at the lower back; and 3) they investigate only a single motor task (Garling et al., 2007; Zemp et al., 2014). The dual fluoroscopy imaging system (DFIS) has been recognized as an accurate and validated in vivo evaluation technique for determining vertebra translations and orientations with the accuracy of 0.3 mm and 0.70° (Wang et al., 2008). Recent studies have demonstrated that the motion tracking technique based on dual fluoroscopic imaging is a reliable method for determining intervertebral motion and adjacent segment kinematics in various functional spine motions (Wawrose et al., 2020; Zhou et al., 2020). Therefore, in this study, we used a DFIS combined with a validated 3D-to-2D registration technique (Wang et al., 2020; Zhou et al., 2021) to quantify in vivo intervertebral and upper-lower lumbar segment kinematics in asymptomatic subjects during forward–backward bending.
The present study aimed to describe patterns and magnitudes of STAs in vivo in three different directions on the full region of the human lumbar spine during weight-bearing forward–backward bending. Vertebral movement through the cycle was measured using the DFIS. It also aimed to determine the effects of STAs on the calculated flexion angles between the upper and lower lumbar spine and adjacent intervertebral levels for L2–L5 by comparing the skin marker (SM)- and virtual bone marker (VM)-based measurements. We hypothesized that the in vivo patterns of skin markers in the lumbar region would be direction- and segment-dependent patterns and that STAs would significantly affect the kinematic variables in human spinal motion investigation.
2 Materials and methods
2.1 Participant preparation
The present study was conducted with the approval of the institutional review board at Tongji Hospital, Shanghai, China (Protocol Number: 2021-011-SK), and followed the guidelines of the Helsinki Declaration (2013) (World Medical Association, 2013). The recruitment was conducted in the personal and workplace environment of the investigators. Ten asymptomatic adults without prior diagnosed spinal disorders were recruited. All the participants provided written informed consent before any personal or health-related data were collected. The inclusion criteria were as follows: age between 18 and 75 years, ability to perform the required functional tasks, and sufficient understanding of Chinese. Individuals were excluded in the case of any history of LBP in the past 6 months; injuries or operations on the spine, hip, knee, or ankle; and any comorbidities or circumstances (e.g., pregnancy) that could limit the forward–backward bending capabilities. Additionally, yoga practitioners, dancers, gymnasts, and those who had physical therapy within the last three months were ineligible owing to potential bias toward bending activities.
All the participants underwent a lumbar spine computed tomography (CT) scan (United Imaging, uCT760, voltage 120 kV, resolution 512 × 512, layer spacing 1.0 mm) in a supine, relaxed position. The CT images were imported into the 3D visualization and modeling software, Amira 6.7 (Thermo Fisher Scientific, Rockford, IL, United States), to reconstruct 3D bone models of the lumbar spine. All images were processed using Digital Imaging and Communications in Medicine (DICOM) and bitmap file formats. Four participants were excluded because of poor biplane image quality of the lumbar spine due to abdominal obesity and difficulties in later registration. Eventually, six pain-free female adults (age: 26.17 ± 3.55 years, body mass index (BMI): 19.76 ± 1.36 kg/m2) were included in this study.
2.2 Experiment procedure
A trained physiotherapist (7 years of experience in motion analysis) equipped all the participants through manual palpation with the subject standing upright. Based on our lumbar marker set (Gombatto et al., 2017; Papi et al., 2019), 13 custom reflective spherical markers (diameter: 14 mm) with a lead core (diameter: 3 mm) were used. The central marker was always positioned on the L1, L3, and L5 spinous processes (Markers 1, 3, and 5), which were identified following the palpation guidelines (Field and Hutchinson, 2012). Six lateral markers were placed 5 cm on either side from the spinous process of the first (Markers L1 and R1), third (Markers L3 and R3), and fifth (Markers L5 and R5) lumbar vertebrae separately. Four markers were attached 2.5 cm to each side of the spinous process of the second (Markers L2 and R2) and fourth (Markers L4 and R4) lumbar vertebrae (Papi et al., 2019). The primary goal of this marker set (Table1) was to constrain the lumbar segment (L1–L5) so that the flexion angles between the upper and lower lumbar spine and adjacent vertebrae could be obtained via kinematic calculation (Figure 1).
TABLE 1. Marker placement and abbreviations.
FIGURE 1. (A) Schematic of marker placement and spine anatomical frames of reference (a1); joint coordinate system axes of rotation for the upper and lower lumbar segment considered (a2), (B) lumbar DR image (b1) and retroreflective marker schematic (b2).
Quasi-static dual-digital radiography (DR) images (TAOiMAGE, Shanghai, China; image resolution 2804 × 2804 pixels) were initially captured during a static standing trial with the subject’s feet shoulder-width apart for two seconds. Next, each participant performed weight-bearing forward–backward bending from maximum extension to approximately 45° of flexion (measured using a protractor) to keep the participants within view of the system (Figure 2). In order to reduce the overlapping of skin markers and various internal bony structures in dual fluoroscopic lumbar spine images and enhance the efficiency and accuracy of registration and recognition, the participants stood in an oblique position of 45°. Before the formal collection, we pasted marked lines at different heights of the participants’ lower limbs, instructed the participants to place their palms in front of their lower limbs, and moved down to the designated marking line in an orderly manner. While reaching the designated marked position, they were asked to cross their hands behind their heads and keep still. Participants were trained to keep both knees straight throughout the process. After the preparation above, anteroposterior and lateral biplanar images of the lumbar vertebrae were recorded synchronously at different bending angles. Participants practiced the activities before the actual testing to get familiar with the tasks. The testing was repeated if the participants violated the task instructions, resulting in non-valid trials. During the experiment, all the participants were required to wear lead protective clothing in non-collection areas to reduce unnecessary radiation. During testing, the participant was exposed to approximately seven pairs of fluoroscopic projections. The entire image collection process took around 10 min to complete. According to the product manual provided by the manufacturer, the effective dose under our typical testing conditions (78 Kv, 40 mA) is 0.16 mSv for each paired image of the lumbar or less than 1.20 mSv per test, approximately 6% of the average dose limit for occupational exposure in a year.
FIGURE 2. Experimental setup of the dual fluoroscopic system for capturing lumbar spine and marker positions in vivo; virtual reproduction of the dual fluoroscopic system and vertebral positions.
2.3 In vivo kinematics measurements
Two multi-segmental kinematic models of the lumbar spine were used to derive lumbar spine kinematics based on virtual bone and skin markers. At each vertebral level and lumbar segment, anatomical coordinate systems were defined based on the recommendations of the International Society of Biomechanics (Wu et al., 2002).
1) At each vertebral level, the origin was in the middle of the upper and lower endplates. The y-axis (proximal/distal, PD) was defined as the line passing through the centers of the upper and lower vertebral endplates. The z-axis (medial/lateral, ML) was defined as the line parallel to a line joining similar landmarks on the bases of the right and left pedicles, pointing to the right. The x-axis (anterior/posterior, AP) was perpendicular to the y- and z-axes and pointed anteriorly.
2) The upper lumbar segment was defined by its origin at L3, a vertical axis from the spinous process of L3 to L1 (PD, +y), an AP axis (+x) orthogonal to the plane formed by the PD axis, the line connecting the left and right transverse processes of L3, and a horizontal axis (ML, +z) cross-product of AP and PD axes. The lower lumbar segment was defined by its origin at L5, a vertical axis from L5 to L3 (PD, +y), an AP axis orthogonal to the plane formed by the PD axis, the line connecting the left and right transverse processes of L5, and a horizontal axis cross-product of the AP and PD axes (Gombatto et al., 2017) [Figure 1A (a2)].
2.4 Data processing and analysis
The STA-free trajectories of the markers (i.e., “virtual bone markers”) were derived from the kinematics of the bones. 2D DR images and 3D models were imported into a registration program (MATLAB, R2018a MathWorks, Natick, MA, United States). The 3D positions of each lumbar vertebra and marker were adjusted until they matched the corresponding outlines on the DR images. The 3D positions of the skin markers relative to the underlying vertebra in each frame were defined. This technique was validated by the roentgen stereophotogrammetric analysis technique as the gold standard for achieving submillimeter accuracy in our preliminary study (Li et al., 2009). Bone landmark positions from CT measurements were transformed into the biplanar DR coordinate system and combined with a skin marker to determine the flexion angles of the upper lumbar relative to the lower lumbar segment. In order to evaluate the higher derivatives of the motion parameter based on the two systems, flexion angles between the adjacent lumbar vertebra (L2–L3, L3–L4, and L4–L5) were also examined. Lumbar joint angles were calculated using the relevant segment poses following a z-x-y cardanic rotation sequence (Robertson et al., 2004; Lees et al., 2010), which corresponded to flexion+/extension−, right bend+/left bend−, and right twist+/left twist−. The effects of STAs on lumbar kinematics were quantified throughout the motion using the differences between the SM-determined and VM-determined kinematics.
Each local coordinate system was located at the center of gravity of the respective vertebral bodies. The global coordinates of the markers were converted into the corresponding local vertebral coordinate system. STA was calculated by subtracting the displacement of the skin marker under the coordinate system of each vertebra from the corresponding marker in the static standing position. The overall flexion angle of the lumbar spine was described as the angle between the fitting line through the center of all vertebrae in each frame and the static standing position. The STA at each tested flexion angle of each subject was normalized to one bending cycle, from full extension (0% of normalized cycle) to maximal flexion (100% of normalized cycle), through data interpolation. The flexion angle ranged from −15° to 45° with an interval of 7.5°. The mean values and range of the STAs in the AP, ML, and PD directions among all the participants were investigated. The mean STAs at the lower lumbar (Markers R5, 5, L5, R4, and L4) and upper lumbar regions (Markers L3, 3, R3, L2, R2, L1, 1, and R1) among all the participants during extension and flexion were calculated. Meanwhile, the STAs at the spinous processes (Markers 1, 3, and 5) and both sides (Markers L1, R1, L2, R2, L3, R3, L4, R4, L5, and R5) were calculated.
2.5 Statistical analysis
Continuous variables are presented as means and standard errors. Normal distribution was verified using the Shapiro Wilk test and Q-Q-plot inspection for discrete parameters. A paired t-test with a significance level of 0.05 was used to compare 1) the STAs at the lower and upper lumbar regions, 2) the STAs at the spinous processes and both sides, and 3) the SM- and VM-determined flexion angles among all the participants at each tested frame. The correlation between the flexion angles based on the two measurements throughout the entire cycle was estimated using the Spearman correlation coefficient.
3 Results
3.1 Lumbar soft tissue artifact patterns and magnitudes in three directions during forward–backward bending
3.1.1 Anterior–posterior direction
The quantitative plots of STAs during the bending cycle demonstrated that the amount of movement of skin markers relative to the underlying bi-planar DR-measured bone STAs at each lumbar level showed a similar increasing trend from full extension to flexion in the AP direction (Figure 3). All the markers were anteriorly shifted during the entire cycle. The smallest STA through the bending cycle occurred at marker spinous processes of L1, ranging from −2.3 ± 3.3 mm at full extension to 4.7 ± 1.9 mm at maximal flexion. The largest STAs occurred at the marker right side of the spinous process of L3, ranging from −6.9 ± 4.2 mm to 17.8 ± 5.9 mm (Table 2).
TABLE 2. Direction-related (AP, PD, and ML) mean marker artifact (mean) with their standard deviations (SDs) of all the lumbar skin markers in the flexed and extended positions.
FIGURE 3. STA of the thirteen skin markers from L1 to L5 level during the bending cycle (0%, 25%, and 100% represent extension, neutral, and flexion, respectively). The solid line represents the mean STA for all the subjects, and the shaded-area error bands represent the standard deviation of the STA. The vertical dotted line represents the neutral position.
3.1.2 Proximal–distal direction
In the PD direction, the STA exhibited an approximately reciprocal pattern during the lumbar extension and flexion (Figure 3). The STA underwent a sharp increase during flexion at the L5 level (Markers L5, 5, and R5 shifted 16.3 ± 5.2 mm, 14.8 ± 3.8 mm, and 19.1 ± 5.2 mm proximal from the static position to maximal flexion, respectively), followed by the L4 level (Markers L4 and R4 shifted 9.5 ± 4.3 mm and 8.0 ± 6.5 mm proximal from the static position to maximal flexion, respectively), and slightly increased at the L1-L3 level (maximal STA < 6.5 mm). At the extension phase, the STA in the PD direction slightly decreased at the L3–L5 level (maximal STA < 4.9 mm), and the greatest decrease occurred at the L1 level (Markers L1, 1, and R1 shifted 11.4 ± 5.6 mm, 8.9 ± 4.0 mm, and 8.9 ± 6.0 mm distal from full extension to the static position, respectively), followed by L2 (Markers L2 and R2 shifted 5.5 ± 4.1 mm and 4.5 ± 3.5 mm distal from full extension to the static position, respectively) (Table 2). The positions of the markers did not move either to the right or left during the bending motion.
3.1.3 Medial–lateral direction
Throughout the bending cycle, the changing trend of the STA in the ML direction was stable relative to the AP and PD. In the extension phase, the largest STA occurred at Marker R1 (−3.3 ± 5.5 mm), and the smallest STA occurred at Marker R3 (−0.3 ± 2.1 mm). During flexion, the smallest STA occurred at Marker R4 (0.1 ± 2.9 mm), and the largest STA occurred at Marker L1 (−4.0 ± 5.3 mm) (Table 2).
In conclusion, as the gross spinal angle increased, the STA became prominent in the AP and PD directions. Specifically, the skin markers moved from posterior to anterior with respect to the underlying vertebra through the bending cycle in AP, whereas, in the PD direction, it shifted from proximal to distal during extension and from distal to proximal during flexion (Figure 3).
3.2 Quantification and characteristics of lumbar soft tissue artifact in upper and lower lumbar regions during forward–backward bending
The results showed that the range of STAs during flexion was 13.0 ± 3.4, 7.7 ± 5.6, and −0.8 ± 1.8 mm in the AP, PD, and ML directions, respectively, and −5.2 ± 2.0, 4.8 ± 3.1, and −0.7 ± 2.0 mm in the AP, PD, ML directions, respectively, during extension. During flexion, the STA at the lower lumbar (L4, L5: 13.5 ± 6.5 mm) was significantly higher than that at the upper lumbar regions (L1–L3: 4.0 ± 1.5 mm) in the PD direction (p < 0.01). During extension, the STA at the lower lumbar (L4–L5: 2.7 ± 0.7 mm) was significantly less than that at the upper lumbar region (L1–L3: 6.1 ± 3.3 mm) in the PD direction (p < 0.05). The STA at the spinous process was significantly lower than that on both sides in the AP direction during extension and flexion (p = 0.02) (Table 3).
TABLE 3. Segment-related (upper lumbar, lower lumbar) mean marker artifact (mean) with their standard deviations (SD) of the lumbar skin markers in the flexed and extended positions and the 95% confidence interval (95% CI) of difference.
3.3 Effects of soft tissue artifact on the calculation of flexion angles between the upper and lower lumbar spine and adjacent intervertebral levels
Significant effects of STAs on the calculated vertebral flexion angles were observed. VM measurement demonstrated a consistent increase in intervertebral flexion with overall body flexion, with an intervertebral range of motion (ROM) of 2.70 ± 4.23°, 2.30 ± 2.74°, and 3.06 ± 3.95° during extension and 8.62 ± 5.78°, 13.57 ± 3.36°, and 15.39 ± 2.61° during flexion for L2–L3, L3–L4, and L4–L5 levels, respectively. The overall patterns of joint angles measured by SM over the bending cycle were opposite to VM for the L2–L3 (r = −0.95, p < 0.01) and L4–L5 levels (r = −0.62, p = 0.08), but similar for the L3–L4 level (r = 0.98, p < 0.01) (Figure 4). However, the amplitude of variation of the flexion angle at the L3–L4 level was overestimated by SM measurements, with L3 more extended during the 0%–50% cycle and more flexed during the 50%–100% cycle. Significant differences in the joint angles among all the participants occurred during the 0%–25%, 75%–100%; 0%–37.5%, 100%, and 0%–50% cycles for L2–L3, L3–L4, and L4–L5 levels, respectively (p < 0.05). The zero kinematic errors of L2/L3, L3/L4, and L4/L5 under the two measurement methods are at 46%, 59%, and 72%, respectively (Figure 4), indicating that the skin movement characteristics of the upper and lower lumbar segments were different. These discrepancies can be explained by the skin deformation mechanism caused by complex lumbar muscle deformation. Second, it is speculated that the kinematic error of zero for the two measurements at 50% of the bending cycle occurs near the L3 level, which may be related to its anatomical and biomechanical characteristics because L3 is the apex of lumbar lordosis and the center of lumbar motion.
FIGURE 4. [(A) ①–④] In vivo kinematics measured by two measurements at adjacent vertebral levels (L2–L3, L3–L4, L4–L5) and upper and lower lumbar levels (upper to lower). The solid line represents the mean values for all the subjects measured by DR. The dotted line represents the mean values for all the subjects measured by skin markers. Statistical significance between the two measurements is marked by the solid orange line along the x-axis of top-row graphs. [(B) ①–④] Rotational error at adjacent vertebral levels (L2–L3, L3–L4, L4–L5) and upper and lower lumbar levels (upper to lower). The shaded area represents the standard deviation of the values.
The flexion angles between the upper and lower lumbar levels from full extension to maximal flexion demonstrated a similar increasing trend in SM and VM measurements (Figure 4). The VM- and SM-derived flexion angles were 4.99 ± 6.57° and 9.44 ± 11.84° during extension and 23.53 ± 4.66° and 29.00 ± 7.87° during flexion, respectively. A significant difference in the flexion angle was observed among all the participants during the entire cycle (p < 0.05). A linear relationship was found between the flexion angles during the bending cycle and both measurements (r = 1, p < 0.01). The error in the measured flexion angle decreased from the extension to the flexion phase.
4 Discussion
The most important finding of the present study is that the extent of STAs in the lumbar spine differed by anatomical direction, marker location, lumbar segment, and bending phases when the participants performed forward–backward bending under the weight-bearing condition. Specifically, the STA was prominent in the AP and PD directions. As the flexion angle increased, the STA continuously increased in the AP direction and exhibited a reciprocal trend in the PD direction during extension and flexion. The STA in the lower lumbar region was significantly smaller than that in the upper lumbar region during extension and larger during flexion in the PD direction. The STA at the spinous process was significantly smaller than that on either side. These findings were consistent with those of the previous studies that described the factors that influence the STA extent (Fuller et al., 1997; Holden et al., 1997; Schwartz et al., 2004). These specific patterns and magnitudes of STAs in the lumbar spine might help clinicians and physical therapists assess kinematic errors and devise methods to minimize its effects on the calculation of lumbar joint kinematics.
STA is particularly evident in the analysis of the lumbar spine because the markers are close to each other, and small relative movements among them may lead to large errors in the assessment of vertebral position and orientation (Ciavarro et al., 2006). Previous studies showed that some general features might be identified. For instance, for tasks performed by non-obese volunteers and motor tasks that do not involve high accelerations or impacts, the largest proportion of the STA-affecting markers located on a given body segment is closely correlated to the relevant joint movement (Cappozzo et al., 1996; Bonci et al., 2014; Camomilla et al., 2015). Evidence showed that forward bending (Corkery et al., 2014; van Wingerden et al., 2008), lifting (Hemming et al., 2018; Pranata et al., 2018), and walking (Henchoz et al., 2015; Lamoth et al., 2006) are key activities in assessing subjects with NS-LBP, also corresponding to the most frequently reported treatment goals of patients with LBP on the Patient-Specific Functional Scale (Stratford, 1995). In this study, we chose forward–backward bending, which is relatively simple in execution and acquisition, as the daily activity for evaluating MQ in healthy participants. To ensure the accuracy and validity of the registration and final quantification data, we finally selected six participants with a relatively low BMI and good image quality.
To the best of our knowledge, only a few studies have reported STA results in the spinous process (Mo¨rl and Blickhan, 2006; Zemp et al., 2014). They showed that the largest STA was 27.4 mm for marker spinous process L5 in the flexed position, and the largest absolute value of STA was 18.3 mm for marker spinous process L3 in the extended position (Zemp et al., 2014). It is higher than our results for the corresponding marker location and task position (L3: 7.6 mm, L5: 20.1 mm). This may be because our experiment was performed under the weight-bearing condition, which is different from those in previous studies using magnetic resonance imaging MRI or CT machines, which limit the positions of the subject to only passive posture (Mo¨rl and Blickhan, 2006; Zemp et al., 2014). In fact, increasing the load significantly affects the movement of the lumbar joints (Wen et al., 2022), and flexion of the spine under active conditions is the most common mechanical cause of lumbar injuries (Shahbazi Moheb Seraj et al., 2018). During extension, the STA in the lower lumbar region was significantly smaller than that in the upper lumbar region in the PD direction, whereas, in the flexion phase, the STA gradually increased from proximal to distal. The maximum STA reached approximately 20 mm at the L5 level, followed by approximately 10 mm at the L4 level, and fell below 10 mm from L1–L3. A high degree of lumbar STA variability was observed during lumbar extension and flexion. The characteristics of STAs at each specific lumbar level may help quantitatively describe the lumbar segment-related kinematic errors caused by STAs, thus providing a better understanding of the normal and pathological biomechanics of the lumbar spine during functional activities.
Various marker placement methods have already been used to track the movement of the human lumbar spine, which is mainly separated into upper and lower segments, to study the kinematic characteristics of LBP (Gombatto et al., 2015; Hemming et al., 2018; Papi et al., 2019). Our study is the first to investigate the in vivo 3D STA of markers on either side of the spinous process. The overall direction of the marker shift was similar to that reported in a previous study (Zemp et al., 2014). We demonstrated that the 3D STA at the spinous process was significantly smaller than that on either side, which means that markers at anatomical landmarks would be least affected by lumbar motion. This is consistent with not only the STA’s markers’ location-related characteristics but also the theoretical conjecture. Because the marker at the spinous process is the closest to the osseous structure of the internal spinous process, there are relatively few soft tissues, such as muscle and fat, between the two, which ultimately results in the least amount of surface-skin deformation caused by internal-muscle deformation.
The present study demonstrated a linear correlation between the flexion angle measured using skin markers and dual DR focused on the upper and lower lumbar segments. This result supports previous findings that the marker set is reliable in assessing lumbar kinematics regardless of the STAs (Seerden et al., 2016). On the contrary, intervertebral rotations measured by the skin markers in our study were far from the intrinsic values measured by fluoroscopy for the L2–L3 and L4–L5 levels. Although a high correlation between the flexion angles at the L3–L4 level (r = 0.98, p < 0.01) derived by the two methods was found, the irregular STA of markers on different lumbar levels in the AP direction together led to a significant overestimation of the flexion angles. The consistency of marker positions and vertebral motion at the L3–L4 level was aligned with some previous studies (Mo¨rl and Blickhan, 2006; Hashemirad et al., 2013); however, the errors should be considered when performing marker-based motion tracking on the lumbar spine. Moreover, given that the overall distribution of the STAs in the lumbar area is lacking, the effectiveness of our markers in estimating the joint angles needs further verification.
This study had several limitations. First, only young women with relatively low BMI were included in the experiment. This experiment was designed to maintain high image quality for accurate in vivo vertebral-motion tracking with thinner fat on the back of the participants. However, large individual differences still exist because of the heterogeneity of each person. While the STA and its effects on spinal kinematics may be affected by the subject’s BMI, age, and sex (Tsai et al., 2011), participants of a larger age and male sex could be involved in the future. Second, to improve the accuracy of STAs and reduce unnecessary radiation, we focused on the flexion-extension direction for acquisition and analysis; therefore, the distribution characteristics and quantitative results of lumbar spine STAs in other degrees of freedom were lacking. However, learning the displacement of the marker set in the full region of the lumbar forward–backward bending could provide a useful reference for estimating 6DOF lumbar motions using SMs. Third, the quasi-static images captured in the current study presented the STA during different bending phases but still could not represent the instantaneous dynamic positions of the marker relative to the vertebra. However, “quasi-static” indicates a dynamic effect, and quasi-static load indicates that the load is applied slowly such that the structure deforms significantly slowly (at a considerably low strain rate); therefore, the inertial force is significantly small and can be ignored. Last but not least, the experimental collection was carried out asynchronously under the motion capture and DFIS, and the errors in OMC measurement differ from the SM measured radiographically, which might limit the utility of this study in future research. However, the quantitative results and patterns of STAs on the lumbar spine might help develop factors for use in global optimization algorithms aimed at minimizing the effects of STAs on the building of specific lumbar kinematic models and calculation of lumbar joint kinematics in the future.
In summary, the extent of STAs in the lumbar spine differed with respect to anatomical direction, marker location, lumbar segment, and flexion phases when the participants performed forward–backward bending under the weight-bearing condition. The STA was prominent in the AP and PD directions during the forward–backward bending motion. As the flexion angle increased, the STA continuously increased in the AP direction and exhibited a reciprocal trend in the PD direction during extension and flexion. The STA in the lower lumbar region was significantly smaller than that in the upper lumbar region during extension and larger during flexion in the PD direction. Skin markers at the spinous process were least affected by the motion of the lumbar spine. The STA characteristics at different lumbar vertebral levels may help clinicians and physical therapists assess kinematic errors and devise methods to minimize its effects on the calculation of lumbar joint kinematics. Credible and valid lumbar spine kinematic results are important for future studies aiming to provide personalized clinical management for LBP.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the institutional review board at Tongji Hospital, Shanghai, China. The patients/participants provided their written informed consent to participate in this study.
Author contributions
XX designed the research methods, performed the experiments, analyzed the data, and drafted the manuscript. ZL contributed to acquisition, interpretation of data, and critical review of the manuscript. CW contributed to critical review of the manuscript. CG performed the experiments. XZ, HY, and SL participated in the experiments. T-YT and YY designed the research, interpreted data, and revised the manuscript. LC revised the manuscript. All authors meet the International Committee of Medical Journal Editors (ICMJE) criteria for authorship for this article and are responsible for the integrity of the work as a whole. All authors have read and approved the final submitted manuscript.
Funding
This study was funded by the Clinical Research Plan of SHDC (SHDC2020CR 2035B). The funder had no involvement in the study design, data analysis, and interpretation of the results.
Conflict of interest
Author T-YT was employed by TaoImage Medical Technologies Corporation.
The remaining 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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Airaksinen, O., Brox, J. I., Cedraschi, C., Hildebrandt, J., Klaber-Moffett, J., Kovacs, F., et al. (2006). Chapter 4. European guidelines for the management of chronic nonspecific low back pain. Eur. Spine J. 15 (2), S192–S300. doi:10.1007/s00586-006-1072-1
Aiyangar, A. K., Zheng, L., Tashman, S., Anderst, W. J., and Zhang, X. (2014). Capturing three-dimensional in vivo lumbar intervertebral joint kinematics using dynamic stereo-X-ray imaging. J. Biomech. Eng. 136, 011004. doi:10.1115/1.4025793
Akbarshahi, M., Schache, A. G., Fernandez, J. W., Baker, R., Banks, S., and Pandy, M. G. (2010). Non-invasive assessment of soft-tissue artifact and its effect on knee joint kinematics during functional activity. J. Biomech. 43, 1292–1301. doi:10.1016/j.jbiomech.2010.01.002
Benoit, D. L., Ramsey, D. K., Lamontagne, M., Xu, L., Wretenberg, P., and Renström, P. (2006). Effect of skin movement artifact on knee kinematics during gait and cutting motions measured in vivo. Gait Posture 24, 152–164. doi:10.1016/j.gaitpost.2005.04.012
Bonci, T., Camomilla, V., Dumas, R., Chèze, L., and Cappozzo, A. (2014). A soft tissue artefact model driven by proximal and distal joint kinematics. J. Biomech. 47, 2354–2361. doi:10.1016/j.jbiomech.2014.04.029
Camomilla, V., and Bonci, T. (2020). A joint kinematics driven model of the pelvic soft tissue artefact. J. Biomech. 111, 109998. doi:10.1016/j.jbiomech.2020.109998
Camomilla, V., Bonci, T., Dumas, R., Chèze, L., and Cappozzo, A. (2015). A model of the soft tissue artefact rigid component. J. Biomech. 48, 1752–1759. doi:10.1016/j.jbiomech.2015.05.007
Cappozzo, A., Catani, F., Leardini, A., Benedetti, M. G., and Croce, U. D. (1996). Position and orientation in space of bones during movement: Experimental artefacts. Clin. Biomech. (Bristol, Avon) 11, 90–100. doi:10.1016/0268-0033(95)00046-1
Cereatti, A., Bonci, T., Akbarshahi, M., Aminian, K., Barré, A., Begon, M., et al. (2017). Standardization proposal of soft tissue artefact description for data sharing in human motion measurements. J. Biomech. 62, 5–13. doi:10.1016/j.jbiomech.2017.02.004
Ciavarro, G. L., Santambrogio, G. C., and Andreoni, G. (2006). Assessment of soft tissue artifacts for flexion/extension movements of the lumbar spine. J. Biomech. 39, S172–S173. doi:10.1016/s0021-9290(06)83605-5
Corkery, M. B., O’Rourke, B., Viola, S., Yen, S.-C., Rigby, J., Singer, K., et al. (2014). An exploratory examination of the association between altered lumbar motor control, joint mobility and low back pain in athletes. Asian J. Sports Med. 5, e24283. doi:10.5812/asjsm.24283
Dal Maso, F., Blache, Y., Raison, M., Lundberg, A., and Begon, M. (2016). Glenohumeral joint kinematics measured by intracortical pins, reflective markers, and computed tomography: A novel technique to assess acromiohumeral distance. J. Electromyogr. Kinesiol. 29, 4–11. doi:10.1016/j.jelekin.2015.07.008
Delitto, A., George, S. Z., Van Dillen, L., Whitman, J. M., Sowa, G., Shekelle, P., et al. (2012). Low back pain. J. Orthop. Sports Phys. Ther. 42, A1–A57. doi:10.2519/jospt.2012.42.4.A1
Field, D., and Hutchinson, J. O. (2012). Field’s anatomy, palpation & surface markings. 5th ed. London, England: Churchill Livingstone.
Fuller, J., Liu, L.-J., Murphy, M. C., and Mann, R. W. (1997). A comparison of lower-extremity skeletal kinematics measured using skin- and pin-mounted markers. Hum. Mov. Sci. 16, 219–242. doi:10.1016/s0167-9457(96)00053-x
Garling, E. H., Kaptein, B. L., Mertens, B., Barendregt, W., Veeger, H. E. J., Nelissen, R. G. H. H., et al. (2007). Soft-tissue artefact assessment during step-up using fluoroscopy and skin-mounted markers. J. Biomech. 40 (1), S18–S24. doi:10.1016/j.jbiomech.2007.03.003
Gombatto, S. P., Brock, T., DeLork, A., Jones, G., Madden, E., and Rinere, C. (2015). Lumbar spine kinematics during walking in people with and people without low back pain. Gait Posture 42, 539–544. doi:10.1016/j.gaitpost.2015.08.010
Gombatto, S. P., D’Arpa, N., Landerholm, S., Mateo, C., O’Connor, R., Tokunaga, J., et al. (2017). Differences in kinematics of the lumbar spine and lower extremities between people with and without low back pain during the down phase of a pick-up task, an observational study. Musculoskelet. Sci. Pract. 28, 25–31. doi:10.1016/j.msksp.2016.12.017
Hashemirad, F., Hatef, B., Jaberzadeh, S., and Ale Agha, N. (2013). Validity and reliability of skin markers for measurement of intersegmental mobility at L2-3 and L3-4 during lateral bending in healthy individuals: A fluoroscopy study. J. Bodyw. Mov. Ther. 17, 46–52. doi:10.1016/j.jbmt.2012.04.010
Hemming, R., Sheeran, L., van Deursen, R., and Sparkes, V. (2018). Non-specific chronic low back pain: Differences in spinal kinematics in subgroups during functional tasks. Eur. spine J. 27 (1), 163–170. doi:10.1007/s00586-017-5217-1
Henchoz, Y., Soldini, N., Peyrot, N., and Malatesta, D. (2015). Energetics and mechanics of walking in patients with chronic low back pain and healthy matched controls. Eur. J. Appl. Physiol. 115, 2433–2443. doi:10.1007/s00421-015-3227-4
Hodges, P. W., Van Dillen, L. R., Mcgill, S., Brumange, S., and Hides, J. A. (2013). “State-of-the-art approach to clinical rehabilitation of low back and pelvic pain,” in The rehabilitation of back pain: State of the art and science. Editors P. W. Hodges, J. Cholewicki, and J. H. Van Dieën (Edinburgh: Elsevier).
Holden, J. P., Orsini, J. A., Siegel, K. L., Kepple, T. M., Gerber, L. H., and Stanhope, S. J. (1997). Surface movement errors in shank kinematics and knee kinetics during gait. Gait Posture 5, 217–227. doi:10.1016/s0966-6362(96)01088-0
Hoy, D., Bain, C., Williams, G., March, L., Brooks, P., Blyth, F., et al. (2012). A systematic review of the global prevalence of low back pain. Arthritis Rheum. 64, 2028–2037. doi:10.1002/art.34347
Huang, X., Lin, D., Liang, Z., Deng, Y., He, Z., Wang, M., et al. (2021). Mechanical parameters and trajectory of two Chinese cervical manipulations compared by a motion capture system. Front. Bioeng. Biotechnol. 9, 714292. doi:10.3389/fbioe.2021.714292
Itz, C. J., Geurts, J. W., van Kleef, M., and Nelemans, P. (2013). Clinical course of non-specific low back pain: A systematic review of prospective cohort studies set in primary care: Clinical course of non-specific low back pain. Eur. J. Pain 17, 5–15. doi:10.1002/j.1532-2149.2012.00170.x
Laird, R. A., Kent, P., and Keating, J. L. (2012). Modifying patterns of movement in people with low back pain -does it help? A systematic review. BMC Musculoskelet. Disord. 13, 169. doi:10.1186/1471-2474-13-169
Lamoth, C. J. C., Meijer, O. G., Daffertshofer, A., Wuisman, P. I. J. M., and Beek, P. J. (2006). Effects of chronic low back pain on trunk coordination and back muscle activity during walking: Changes in motor control. Eur. Spine J. 15, 23–40. doi:10.1007/s00586-004-0825-y
Leardini, A., Chiari, L., Della Croce, U., and Cappozzo, A. (2005). Human movement analysis using stereophotogrammetry. Gait Posture 21, 212–225. doi:10.1016/j.gaitpost.2004.05.002
Lees, A., Barton, G., and Robinson, M. (2010). The influence of Cardan rotation sequence on angular orientation data for the lower limb in the soccer kick. J. Sports Sci. 28, 445–450. doi:10.1080/02640410903540352
Li, G., Wang, S., Passias, P., Xia, Q., Li, G., and Wood, K. (2009). Segmental in vivo vertebral motion during functional human lumbar spine activities. Eur. Spine J. 18, 1013–1021. doi:10.1007/s00586-009-0936-6
Liguori, G.American College of Sports Medicine (ACSM) (2021). ACSM’s guidelines for exercise testing and prescription. 11th ed. Baltimore, MD: Wolters Kluwer Health.
Maher, C., Underwood, M., and Buchbinder, R. (2017). Non-specific low back pain. Lancet 389, 736–747. doi:10.1016/s0140-6736(16)30970-9
Manal, K., McClay, I., Stanhope, S., Richards, J., and Galinat, B. (2000). Comparison of surface mounted markers and attachment methods in estimating tibial rotations during walking: An in vivo study. Gait Posture 11, 38–45. doi:10.1016/s0966-6362(99)00042-9
Metcalf, C. D., Phillips, C., Forrester, A., Glodowski, J., Simpson, K., Everitt, C., et al. (2020). Quantifying soft tissue artefacts and imaging variability in motion capture of the fingers. Ann. Biomed. Eng. 48, 1551–1561. doi:10.1007/s10439-020-02476-2
Mörl, F., and Blickhan, R. (2006). Three-dimensional relation of skin markers to lumbar vertebrae of healthy subjects in different postures measured by open MRI. Eur. Spine J. 15, 742–751. doi:10.1007/s00586-005-0960-0
Papi, E., Bull, A. M. J., and McGregor, A. H. (2019). Spinal segments do not move together predictably during daily activities. Gait Posture 67, 277–283. doi:10.1016/j.gaitpost.2018.10.031
Pranata, A., Perraton, L., El-Ansary, D., Clark, R., Mentiplay, B., Fortin, K., et al. (2018). Trunk and lower limb coordination during lifting in people with and without chronic low back pain. J. Biomech. 71, 257–263. doi:10.1016/j.jbiomech.2018.02.016
Robertson, D. G. E., Hamill, J., Caldwell, G. E., Kamen, G., and Whittlesey, S. (2004). Research methods in biomechanics. Champaign, IL: Human Kinetics.
Sangeux, M., Marin, F., Charleux, F., Dürselen, L., and Ho Ba Tho, M. C. (2006). Quantification of the 3D relative movement of external marker sets vs. bones based on magnetic resonance imaging. Clin. Biomech. (Bristol, Avon) 21, 984–991. doi:10.1016/j.clinbiomech.2006.05.006
Schwartz, M. H., Trost, J. P., and Wervey, R. A. (2004). Measurement and management of errors in quantitative gait data. Gait Posture 20, 196–203. doi:10.1016/j.gaitpost.2003.09.011
Seerden, Stefan, Dankaerts, Wim, and Vanwanseele, Benedicte (2016). Reliability of a multi-segmental kinematic model to analyse 3D movement in the lumbar spine. Gait Posture. 49 53. doi:10.1016/j.gaitpost.2016.07.117
Shahbazi Moheb Seraj, M., Sarrafzadeh, J., Maroufi, N., Ebrahimi Takamjani, I., Ahmadi, A., and Negahban, H. (2018). The ratio of lumbar to hip motion during the trunk flexion in patients with mechanical chronic low back pain according to O’sullivan classification system: A cross-sectional study. Arch. Bone Jt. Surg. 6, 560–569.
Stratford, P. (1995). Assessing disability and change on individual patients: A report of a patient specific measure. Physiother. Can. 47, 258–263. doi:10.3138/ptc.47.4.258
Tsai, T.-Y., Lu, T.-W., Kuo, M.-Y., and Lin, C.-C. (2011). Effects of soft tissue artifacts on the calculated kinematics and kinetics of the knee during stair-ascent. J. Biomech. 44 (6), 1182–1188. doi:10.1016/j.jbiomech.2011.01.009
van Wingerden, J.-P., Vleeming, A., and Ronchetti, I. (2008). Differences in standing and forward bending in women with chronic low back or pelvic girdle pain: Indications for physical compensation strategies. Spine 33, E334–E341. doi:10.1097/brs.0b013e318170fcf6
Wallbott, H. G. (1989). “Movement quality changes in psychopathological disorders,” in Normalities and abnormalities in human movement. Editor B. Kirkcaldy (Basel: Karger), 128–146. doi:10.1159/000416532
Wang, H., Zhou, C., Yu, Y., Wang, C., Tsai, T. Y., Han, C., et al. (2020). Quantifying the ranges of relative motions of the intervertebral discs and facet joints in the normal cervical spine. J. Biomech. 112, 110023. doi:10.1016/j.jbiomech.2020.110023
Wang, J., Lui, Z., Qian, Z., and Ren, L. (2016). Soft tissue artifact evaluation of the cervical spine in motion patterns of flexion and lateral bending: A preliminary study. PeerJ 4, e1893. doi:10.7717/peerj.1893
Wang, S., Passias, P., Li, G., Li, G., and Wood, K. (2008). Measurement of vertebral kinematics using noninvasive image matching method-validation and application. Spine 33, E355–E361. doi:10.1097/brs.0b013e3181715295
Wawrose, R. A., LeVasseur, C. M., Byrapogu, V. K., Dombrowski, M. E., Donaldson, W. F., Shaw, J. D., et al. (2020). In vivo changes in adjacent segment kinematics after lumbar decompression and fusion. J. Biomechanics 26, 109515. doi:10.1016/j.jbiomech.2019.109515
Wen, W., Xu, H., Zhang, Z., Kou, B., Sun, Q., and Miao, J. (2022). The effect of various weight-bearing activities on the motion of lumbar facet joints in vivo. J. Orthop. Surg. Res. 17, 114. doi:10.1186/s13018-022-03016-5
World Medical Association (2013). World medical association declaration of Helsinki: Ethical principles for medical research involving human subjects: Ethical principles for medical research involving human subjects. JAMA.310, 2191–2194. doi:10.1001/jama.2013.281053
Wu, G., Siegler, S., Allard, P., Kirtley, C., Leardini, A., Rosenbaum, D., et al. (2002). ISB recommendation on definitions of joint coordinate system of various joints for the reporting of human joint motion—part I: Ankle, hip, and spine. J. Biomech. 35, 543–548. doi:10.1016/s0021-9290(01)00222-6
Zemp, R., List, R., Gülay, T., Elsig, J. P., Naxera, J., Taylor, W. R., et al. (2014). Soft tissue artefacts of the human back: Comparison of the sagittal curvature of the spine measured using skin markers and an open upright MRI. PLoS ONE 9, e95426. doi:10.1371/journal.pone.0095426
Zhou, C., Li, G., Wang, C., Wang, H., Yu, Y., Tsai, T. Y., et al. (2021). In vivo intervertebral kinematics and disc deformations of the human cervical spine during walking. Med. Eng. Phys. 87, 63–72. doi:10.1016/j.medengphy.2020.11.010
Keywords: lumbar spine, forward–backward bending, in vivo kinematics, soft tissue artifacts, dual fluoroscopy
Citation: Xi X, Ling Z, Wang C, Gu C, Zhan X, Yu H, Lu S, Tsai T-Y, Yu Y and Cheng L (2022) Lumbar segment-dependent soft tissue artifacts of skin markers during in vivo weight-bearing forward–Backward bending. Front. Bioeng. Biotechnol. 10:960063. doi: 10.3389/fbioe.2022.960063
Received: 02 June 2022; Accepted: 18 July 2022;
Published: 17 August 2022.
Edited by:
Cheng-Feng Lin, National Cheng Kung University, TaiwanReviewed by:
Alexander Breen, Bournemouth University, United KingdomAlex Kernkamp, Leiden University Medical Center, Netherlands
Copyright © 2022 Xi, Ling, Wang, Gu, Zhan, Yu, Lu, Tsai, Yu and Cheng. 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.
*Correspondence: Tsung-Yuan Tsai, tytsai@sjtu.edu.cn; Yan Yu, yyu15@tongji.edu.cn
†These authors share first authorship