A biomechanical comparison of 360° stabilizations for corpectomy and total spondylectomy: a cadaveric study in the thoracolumbar spine
© Kim et al. 2015
Received: 28 February 2015
Accepted: 15 June 2015
Published: 1 July 2015
To date, there has been no adequate biomechanical model that would allow a quantitative comparison in terms of stability/stiffness between a corpectomy with the posterior column preserved and a total spondylectomy with the posterior column sacrificed. The objective of this study was to perform a biomechanical comparison of 360° stabilizations for corpectomy and total spondylectomy, using the human thoracolumbar spine.
Five human cadaveric thoracolumbar spines (T8-L2) were tested according to the following loading protocol: axial compression, flexion, extension, lateral bending to the right and left, and axial rotation to the right and left. This loading protocol was applied three times. Each specimen was tested intact, after corpectomy, and after total spondylectomy. The relative stiffness of each motion segment was determined for each test.
There was no significant difference in stiffness after reconstruction of total spondylectomy versus corpectomy in our thoracolumbar model. Our construct consisted of an anterior cage and four-level pedicle screw instrumentation (two above and two below) and provided similar stiffness in both models. Despite the additional bone resection in a total spondylectomy versus corpectomy, the constructs did not differ biomechanically. Additionally, there was no significant difference in stiffness between the intact specimen and either reconstruction model.
A classic corpectomy, which leaves the posterior column intact, is no better in terms of stability/stiffness than a total spondylectomy carried out using a shorter cage, followed by compression using posterior instrumentation.
KeywordsCorpectomy Total spondylectomy Biomechanics Thoracolumbar spine Stiffness Stability
Spine surgeons have used combined anterior and posterior surgery for the treatment of patients with spinal fractures, tumors, and deformities. Corpectomy is generally followed by anterior reconstruction with a strut graft or a cage and posterior instrumentation with pedicle screws. This procedure is performed to provide sufficient mechanical stability. However, total spondylectomy (complete resection of the affected vertebra, including the posterior column) by means of a posterior approach is sometimes preferred, because pedicle screws and a posterior approach technique can provide secure fixation. Posterior total spondylectomy can be useful in cases where the anterior approach is considered difficult, such as in cases of severe adhesions due to previous anterior surgery, suboptimal pulmonary function, severe kyphosis, and trunk shortening [1–3].
Classic corpectomy leaves the posterior column intact, thus, offering good biomechanical stability, since it provides support at both the front and back [4–6]. On the other hand, total spondylectomy sacrifices the posterior column and leads to complete loss of spinal continuity and stability during the operation. However, there is a report on total spondylectomy using a 10-mm shorter cage, followed by compression using posterior instrumentation; this resulted in spinal shortening and stability . Ten millimeters of shortening is within the acceptable range for the spinal cord [8–10]. Despite the effectiveness of shortening, spine surgeons may question whether total spondylectomy with posterior instrumentation can provide stability comparable to combined anterior corpectomy and posterior instrumentation.
To date, there has been no adequate biomechanical model that would allow a quantitative comparison in terms of stability/stiffness between a corpectomy and a total spondylectomy. We decided to compare biomechanical stabilities using a thoracolumbar reconstructed spine (T9-L1), with anterior cage and multilevel posterior instrumentation after corpectomy (T11), and the same spine with a 10-mm shorter anterior cage and multilevel posterior instrumentation after total spondylectomy (T11).
Materials and methods
Thoracolumbar spines (T8-L2) were harvested from five fresh human cadavers (mean age, 82 years; range, 70–95 years; two females, three males). The specimens were examined grossly and radiographically to rule out malignancy or fractures that could have interfered with the results. The specimens were then frozen at −20 °C until the day before testing. After thawing to room temperature overnight, the surrounding soft tissues and muscles were dissected meticulously to preserve the osseous and ligamentous structures. The specimens were kept moist throughout the procedure.
Bone mineral density
The bone mineral density (BMD) of all specimens was measured in a water bath, using dual-energy X-ray absorptiometry (GE Lunar Prodigy, GE Medical Systems, Madison, WI, USA). Using the World Health Organization definition of osteoporosis as BMD <0.8 g/cm2, all of the vertebrae were classified as osteoporotic.
The intact specimen
Step 1 (establishing the center of rotation): The center of rotation for flexion-extension and lateral bending was established as follows. A pure compressive load of 50 N was applied through a roller bearing to the top surface of the upper end-cap containing T8. With the roller bearing positioned, the load was applied, then repositioned and applied again, until no angular rotation in the coronal or sagittal planes could be detected when a compressive force was applied to the upper end-cap. This spot on the top surface of the end-cap was deemed to be the center of rotation and remained as the reference spot in subsequent testing on that specific specimen. A cyclic compression force was then applied as conditioning (40 ± 10 N at 0.5 Hz for 15 min) to remove excess fluid from the discs and to return the disc to its predeath height .
Step 2 (compression): The specimen was then loaded in pure compression (with the load applied at the center of rotation) at a displacement rate of 0.25 cm/min. The load was applied up to a maximum of 200 N. A load-deformation curve was obtained. The specimen was then loaded two more times, and load-deformation curves were obtained each time. The three load-deformation curves were essentially identical, and only one was used to calculate the stiffness. Stiffness was calculated as the slope of the linear portion of the load-deformation curve .
Step 3 (flexion): The test was repeated with a 100 N load applied 2 cm anterior to the center of rotation, to produce a maximum bending moment of 2 Nm (i.e., 0.02 m × 100 N = 2 Nm) in flexion. The load was then released and applied twice more. The three load-deformation curves were essentially identical, and only one was used to calculate the stiffness.
Step 4 (extension): The test was repeated with a 100 N load applied 2 cm posterior to the center of rotation, to send the specimen into extension. The load was then released and applied twice more. The three load-deformation curves were essentially identical, and only one was used to calculate the stiffness.
Step 5 and 6 (lateral bending): The test was repeated with a 100 N load applied 2 cm to the right of the center of rotation, to send the specimen into right lateral bending. The load was then released and applied twice more. Three load-deformation curves were obtained for right lateral bending. The three right load-deformation curves were essentially identical, and only one was used to calculate the stiffness. The experiment was then repeated for left lateral bending, and the same criteria applied.
Step 7 and 8 (right and left axial rotation): To apply axial rotation, the specimen was first compressed to 100 N, and an axial torque was then applied in a clockwise motion (about the center of rotation) to a maximum of 10°. Three torque-angle deformation curves were obtained for right rotation. The three right torque-angle deformation curves were essentially identical, and only one was used to calculate the stiffness. The experiment was then repeated for left rotation, and the same criteria applied.
The corpectomy model
After the test on the intact specimen had been carried out (as described above), pedicle screws were inserted at T9, T10, T12, and L1 (5.5-mm diameter pedicle screws, 40-mm long at T9 and T10; 6.5-mm diameter pedicle screws, 45-mm long at T12 and L1 (Expedium System, Depuy Spine., Raynham, MA, USA)). A corpectomy at T11 was performed by incising and removing the T10-11 and T11-12 intervertebral discs and the cartilaginous endplates. The entire T11 vertebral body was removed after creating bilateral osteotomies through the T11 pedicles and resecting the anterior and posterior longitudinal ligaments at T11; the posterior elements were left intact. Vertebral body reconstruction was performed by implanting a Harms titanium mesh cage (Harms Cage, Depuy Spine) between T10 and T12, depending on the corpectomy height (Figs. 2 and 3). The diameter of each cage was 16 mm. After the corpectomy reconstruction had been assembled, the biomechanical tests (step 2 through step 8 as described above) were carried out on the specimen, and stiffness curves were obtained, as described above.
The total spondylectomy model
After the test on the corpectomy model had been carried out, the rod on one side was removed to allow access for the removal of the posterior elements. A total spondylectomy at T11 was carried out by removing the posterior elements. A cage 10 mm shorter than the one used in the corpectomy model was inserted at T11. After the “short” cage was inserted into the anterior defect, the posterior instrumentation was adjusted, using a compressor and rod-screw junctions, to compress the inserted cage (Figs. 2 and 3). With the total spondylectomy model assembled, the biomechanical tests (step 2 through step 8 as described above) were then carried out on the specimen, and stiffness curves were obtained, as described above.
Thus, each specimen was biomechanically tested three times: intact, as a corpectomy model, and as a total spondylectomy model. The order of assembling and testing the models was the same for all five specimens.
Results for stiffness are presented as the mean of determinations, with error bars representing standard deviations. Means and standard deviations were computed for the five specimens, for each measure under each experimental condition. A repeated measures analysis of variance with the three models (intact, corpectomy, and total spondylectomy) was performed, followed by post hoc pairwise comparisons based on a Bonferroni correction. Post hoc comparisons were made for the intact spine versus each of the two models and for the corpectomy model versus the total spondylectomy model. The level of significance was P < 0.05. For statistical analysis, SAS version 9.2 (Cary, NC, USA) was used.
None of the specimens fractured, and none of the pedicle screws or rods loosened during the experiments. Visual inspection revealed no major pedicle violations. There was one minor lateral perforation, with threads barely visible outside of the pedicle, but this did not compromise fixation.
Means and standard deviation of stiffness
Axial compression (N/mm)
235.0 ± 63.3
243.6 ± 45.5
240. 1 ± 60.6
149.4 ± 38.6
168.1 ± 53.6
219.9 ± 78.4
228.2 ± 85.4
230.5 ± 98.2
246.3 ± 66.7
Lateral bending (N/mm)
171.8 ± 51.3
204.8 ± 50.9
179.8 ± 26.9
Right rotation (Nmm/degree)
580 ± 270
730 ± 260
690 ± 180
Left rotation (Nmm/degree)
620 ± 310
710 ± 260
680 ± 200
No significant difference was observed between the two reconstruction models, or between the intact specimen and either of the two reconstruction models.
Both reconstruction models provided greater stiffness than that shown by the intact specimen, but the differences were not significant. The total spondylectomy model provided greater stiffness than the corpectomy model, but the difference was not significant.
The two reconstruction models provided greater stiffness than the intact specimen, but the differences were not significant. The total spondylectomy model provided greater stiffness than the corpectomy model, but the difference was not significant.
The values of stiffness measured for right and left lateral bending were averaged for each specimen, since each of our specimens showed some slight scoliotic changes (three left scoliosis, two right scoliosis). Both of the reconstruction models provided greater stiffness than that shown by the intact specimen, but the differences were not significant. The corpectomy model provided greater stiffness than the total spondylectomy model, but the difference was not significant.
Left axial rotation
The two reconstruction models provided greater stiffness than that shown by the intact specimen, but the differences were not significant. The corpectomy model provided greater stiffness than the total spondylectomy model, but the difference was not significant.
Right axial rotation
The two reconstruction models provided greater stiffness than that shown by the intact specimen, but the differences were not significant. The corpectomy model provided greater stiffness than the total spondylectomy model, but the difference was not significant.
Spinal pathologic conditions, such as fractures, tumors, and degenerative diseases, lead to spinal instability, deformity, and neural element compression [12–14]. The question is, How should corpectomy be approached and stabilized (anteriorly, posteriorly, or combined anteroposteriorly)? This has been the subject of debate for some time. Classic corpectomy in a piecemeal fashion leaves the posterior column intact, providing up to 40 % of residual segmental stability [4–6]. Various authors assessed stiffness of several stabilizations using corpectomy and total spondylectomy models; these authors demonstrated that 360° stabilization had higher primary stiffness values, compared to posterior or anterior alone. A 360° stabilization is sometimes preferred, since it stabilizes a spinal segment on both sides of the center of rotation [4, 15–17]. More recently, a combination of anterior strut graft or cage and multilevel posterior instrumentation with pedicle screws by an anterior-posterior approach is commonly performed after corpectomy for the treatment of various spinal pathological conditions, such as fractures, tumors, and deformities. However, when transferring these results to clinical application, a combination of anterior strut graft and multilevel posterior instrumentation after corpectomy requires prolonged operative time and increases the surgical risks, including enlarged surgical approaches, increased cardiopulmonary complications, and greater blood loss [12, 18, 19]. Therefore, posterior total spondylectomy can be useful in cases where the anterior approach is considered difficult, such as in cases of severe adhesions due to previous anterior surgery, suboptimal pulmonary function, severe kyphosis, trunk shortening, and poor general health status [1–3].
Total spondylectomy (complete resection of the affected vertebra, including the posterior column) leads to complete loss of spinal continuity and stability. Therefore, subsequent secure and stable spinal reconstruction is absolutely required. Some authors reported that applied compression forces on the anterior strut graft or cage using posterior instrumentation improves stability [16, 20, 21]. As noted above, Kato carried out a biomechanical study using human cadavers and recommended reconstruction using a cage 10 mm shorter in order to provide more stability after total spondylectomy . One reason for using a 10-mm shorter cage after total spondylectomy in our study is that in an actual spinal surgery, it can be very difficult to replace the removed vertebral body with a cage of the exact height. It is therefore common to use a shorter cage for the removed vertebral body. In addition, a compressive force is usually applied to the cage by means of posterior instrumentation to increase spinal stability in reconstruction. This maneuver results in slight spinal shortening (5–10 mm). A 10-mm shortening is within the safe range for the spinal cord, according to a report by Kawahara et al., who studied the safety limits and physiological effects of acute spinal shortening on the spinal cord in dogs and reported an increase in spinal cord blood flow, as well as no spinal cord injury .
In order to evaluate spinal stability after a total spondylectomy, the stiffness of the reconstructed spine and the stress generated in the instrument used for the spinal reconstruction have been investigated by a finite element method and loading test using a cadaveric spine. Experimentation using a cadaveric spine is a direct means to analyze the behavior of the implant through the displacement response of instrumented segments to applied forces, compared to intact segments. However, experimental results are sensitive to limitations caused by testing procedures and variability among specimens. A finite element model presents a great advantage, because it enables the same procedure (number of vertebrae in the segment, boundary conditions) and the same vertebral segment (geometry, mechanical characteristics) to be used for different spinal implants and for each implant to analyze the influence of a different parameter. Therefore, a validated finite element model would constitute a powerful simulation tool, for the clinician as well as for the implant designer . Nevertheless, as far as we are aware, there is no published report on direct biomechanical comparison of 360° stabilizations for a corpectomy and a total spondylectomy in a finite element method or cadaveric models. This lack of clear-cut, relevant biomechanical data has meant that choosing an anterior corpectomy and a posterior total spondylectomy model has been based mainly on empirical evidence. Therefore, we compared biomechanical stabilities, using a thoracolumbar reconstructed spine (T9-L1) with anterior cage and multilevel posterior instrumentation after corpectomy (T11), and the same spine with a 10-mm shorter anterior cage and multilevel posterior instrumentation after total spondylectomy (T11).
Posterior pedicle screw segmental instrumentation could allow more rigid fixation with three columns. The length of the posterior fixation is a major determinant for rigidity of the construct. Long posterior fixations were significantly more stable compared to intact specimens. The functional relationship between length of posterior fixation and spinal stiffness following spondylectomy and subsequent reconstruction is emphasized by the decreased range of motion in the long posterior fixation group, when compared to short fixation and an intact specimen. In most of the biomechanical studies dealing with corpectomy or total spondylectomy models, an anterior cage combined with multilevel posterior fixation (two above and two below) was able to provide more stiffness than the intact specimens. On the other hand, there is evidence to suggest that eliminating mechanical loads on healing bone when using rigid fixation may result in negative bone remodeling and net bone loss. Thus, less rigid fixation that permits a certain degree of micromotion may accelerate the time to union. Akamaru et al. made a three-dimensional finite element model of the reconstructed spine (T10-L4) following total spondylectomy at T12. They concluded that a reconstruction method with multilevel posterior instrumentation (two above and two below) and no anterior fixation should be better at allowing stress for remodeling of the bone graft inside the titanium mesh cage . Therefore, we concluded that four-level pedicle screw instrumentation (two above and two below) is an adequate method of reconstruction after spondylectomy.
There was no significant difference in stiffness after reconstruction of total spondylectomy versus corpectomy in this cadaveric model. Additionally, there was no significant difference in stiffness between the intact specimen and either of the two reconstruction models. Despite the additional bone resection performed when doing a total spondylectomy versus corpectomy, the constructs did not differ biomechanically. Some authors similarly refer to the importance of applied compression forces on the anterior strut graft to improve stability [16, 20, 21]. The stabilizing effect of such compression in reconstruction using a 10-mm shorter cage after total spondylectomy is supported by the results of our tests. This study made a direct biomechanical comparison between anterior cage and multilevel posterior instrumentation after corpectomy and anterior shorter cage (10 mm) and multilevel posterior instrumentation after total spondylectomy. We would expect that total spondylectomy using a 10-mm shorter cage would restore stiffness to a level equivalent to or greater than that of the intact spine and anterior corpectomy using a cage. In general, postoperative stability after reconstruction is accomplished by spinal fusion. Spinal internal fixation devices provide stability until spinal fusion occurs. Therefore, our own results show that each of the two reconstruction models restored stiffness to a level equivalent to that of spinal fusion in actual spinal surgery. In other words, posterior total spondylectomy using a 10-mm shorter cage can provide sufficient stability.
Several limitations to the study should be considered. 1) The loads and torques chosen to test the specimens in various configurations were low. We chose low loads and torques since there were many testing steps, and we wanted to avoid damage to the specimens before the entire sequence had been completed. 2) Our three models were tested immediately after reconstruction, so conclusions cannot be drawn about the long-term stability in a spine subject to cyclic loads applied daily. 3) The relatively small sample size may also be a limitation, although the results are very clear (no significant differences). The vagaries of postmortem collection always present problems. Most of the spines available for collection and subsequent testing are usually from older cadavers and often have multiple massive osteophytes that render the spines unsuitable for the testing described above. Most of the thoracic spines made available to us did not meet our selection criteria. 4) Thoracolumbar specimens were used in the present study, and only a T11 spondylectomy model was evaluated. Corpectomy and total spondylectomy are performed not only in the lumbar spine but also in the proximal and midthoracic spine. The biomechanical environment of the thoracic spine is quite different from that of the lumbar spine because of the stabilizing effects of the costovertebral joints and rib cage. T11 at the thoracolumbar junction is exposed to forces similar to the lumbar spine. Therefore, reconstruction stability after total spondylectomy in the thoracic spine remains unknown. Further biomechanical studies examining the comparative multisegmental mechanics of total spondylectomy in the thoracic region are required for certainty. 5) The results of the present study do not allow for a recommendation regarding the length of constructs to be used clinically, when performing either total spondylectomy or corpectomy. Depending on the bone purchase achieved by the pedicle screws, as well as the overall sagittal alignment, additional pedicle screw fixation may be needed in clinical practice.
The choice of management for spinal surgery depends on many factors, which surgeons must weigh before determining the ideal surgical approach. These in vitro biomechanical data provide helpful guidelines for planning surgical management.
A classic corpectomy, which leaves the posterior column intact, is no better in terms of stability/stiffness than a total spondylectomy, when the total spondylectomy is carried out using a shorter cage followed by compression using posterior instrumentation.
This work was supported by grant from Inje University, 2011.
- Kim KT, Park KJ, Lee JH. Osteotomy of the spine to correct the spinal deformity. Asian Spine J. 2009;3:113–23.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim SS, Cho BC, Kim JH, Lim DJ, Park JY, Lee BJ, et al. Complications of posterior vertebral resection for spinal deformity. Asian Spine J. 2012;6:257–65.PubMed CentralPubMedView ArticleGoogle Scholar
- Suk SI, Kim JH, Kim WJ, Lee SM, Chung ER, Nah KH. Posterior vertebral column resection for severe spinal deformities. Spine. 2002;27:2374–82.PubMedView ArticleGoogle Scholar
- Disch AC, Luzzati A, Melcher I, Schaser KD, Feraboli F, Schmoelz W. Three-dimensional stiffness in a thoracolumbar en-bloc spondylectomy model: a biomechanical in vitro study. Clin Biomech. 2007;22:957–64.View ArticleGoogle Scholar
- Gradl G. Combined stabilization of thoracolumbar spine fractures. Eur J Trauma. 2006;32:249–52.View ArticleGoogle Scholar
- James KS, Wenger KH, Schlegel JD, Dunn HK. Biomechanical evaluation of the stability of thoracolumbar burst fractures. Spine. 1994;19:1731–40.PubMedView ArticleGoogle Scholar
- Kato S, Murakami H, Higshino K, Okada M, Ito Z, Demura S, et al. The effect of spinal shortening after total en bloc spondylectomy: a biomechanical study in the thoracic spine. J Spinal Disord Tech. 2012;25:E183–90.PubMedView ArticleGoogle Scholar
- Tomita K, Kawahara N, Murakami H, Demura S. Total en bloc spondylectomy for spinal tumors: improvement of the technique and its associated basic background. J Orthop Sci. 2006;11:3–12.PubMed CentralPubMedView ArticleGoogle Scholar
- Kawahara N, Tomita K, Murakami H, Demura S. Total en bloc spondylectomy for spinal tumors: surgical techniques and related basic background. Orthop Clin North Am. 2009;5:47–63.View ArticleGoogle Scholar
- Kawahara N, Tomita K, Kobayashi T, Abdel-Wanis ME, Murakami H, Akamaru T. Influence of acute shortening on the spinal cord: an experimental study. Spine. 2005;30:613–20.PubMedView ArticleGoogle Scholar
- Volkman T, Horton WC, Hutton WC. Transfacet screws with lumbar interbody reconstruction: biomechanical study of motion segment stiffness. J Spinal Disord. 1996;9:425–32.PubMedView ArticleGoogle Scholar
- Bishop FS, Samuelson MM, Finn MA, Bachus KN, Brodke DS, Schmidt MH. The biomechanical contribution of varying posterior constructs following anterior thoracolumbar corpectomy and reconstruction. J Neurosurg Spine. 2010;13:234–9.PubMedView ArticleGoogle Scholar
- Brodke DS, Gollogly S, Bachus KN, Alexander MR, Nguyen BK. Anterior thoracolumbar instrumentation: stiffness and load sharing characteristics of plate and rod systems. Spine. 2003;28:1794–801.PubMedView ArticleGoogle Scholar
- Faro FD, White KK, Ahn JS, Oka RS, Mahar AT, Bawa M, et al. Biomechanical analysis of anterior instrumentation for lumbar corpectomy. Spine. 2003;28:E468.PubMedView ArticleGoogle Scholar
- Knoller SM, Meyer G, Eckhardt C, Lill CA, Schneider E, Linke B. Range of motion in reconstruction situations following corpectomy in the lumbar spine: a question for the thoracolumbar spine. Eur Spine J. 2005;30:E229–35.View ArticleGoogle Scholar
- Oda I, Cunningham BW, Abumi K, Kaneda K, McAfee PC. The stability of reconstruction methods after thoracolumbar total spondylectomy. Spine. 1999;24:1634–8.PubMedView ArticleGoogle Scholar
- Shannon FJ, DiResta GR, Ottaviano D, Castro A, Healey JH, Boland PJ. Biomechanical analysis of anterior poly-methylmethacrylate reconstruction following total spondylectomy for metastatic disease. Spine. 2004;29:2096–102.PubMedView ArticleGoogle Scholar
- Suk SI, Chung ER, Kim JH, Kim SS, Lee JS, Choi WK. Posterior vertebral column resection for severe rigid scoliosis. Spine. 2005;30:1682–7.PubMedView ArticleGoogle Scholar
- Suk SI, Chung ER, Lee SM, Lee JH, Kim SS, Kim JH. Posterior vertebral column resection in fixed lumbosacral deformity. Spine. 2005;30:E703–10.PubMedView ArticleGoogle Scholar
- Knop C, Lange U, Bastian L, Blauth M. Three-dimensional motion analysis with Synex: comparative biomechanical test series with a new vertebral body replacement for the thoracolumbar spine. Eur Spine J. 2000;9:472–85.PubMed CentralPubMedView ArticleGoogle Scholar
- Abumi K, Panjabi MM, Duranceau J. Biomechanical evaluation of spinal fixation devices. III. Stability provided by six spinal fixation devices and interbody bone graft. Spine. 1989;14:1249–55.PubMedView ArticleGoogle Scholar
- Skalli W, Robin S, Lavaste F, Dubousset J. A biomechanical analysis of short segment spinal fixation using a three-dimensional geometric and mechanical model. Spine. 1993;18:536–45.PubMedView ArticleGoogle Scholar
- Akamaru T, Kawahara N, Skamoto J, Yoshida A, Murakami H, Hato T, et al. The transmission of stress to grafted bone inside a titanium mesh cage used in anterior column reconstruction after total spondylectomy: a finite-element analysis. Spine. 2005;30:2783–7.PubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.