Acessibilidade / Reportar erro

Sagittal balance and intervertebral disc composition in patients with low back pain

Abstract

The aim of this study was to verify the relationship between quantitative T2 relaxation measurements of lumbar intervertebral discs (IVDs) and spinopelvic parameters in patients with chronic low back pain. The study was approved by the Clinical Hospital of the Ribeirao Preto Medical School (USP) Ethics Committee, and written consent was obtained from all patients. A total of 455 IVDs from 91 consecutive patients with chronic low back pain were included in this prospective study. All subjects were assessed using the Oswestry Disability Index and visual analogue scale questionnaires and were confirmed to have no other spine diseases except disc degeneration. Spinopelvic parameters including the pelvic incidence (PI), pelvic tilt (PT), sacral slope (SS), sagittal vertical axis (SVA), global tilt (GT), T1 pelvic angle (TPA), lumbar lordosis (LL), thoracic kyphosis (TK), pelvic incidence minus lumbar lordosis mismatch (PI-LL), and lack of lumbar lordosis (LLL) were measured. The study group was categorized according to the Roussouly classification. Sagittal T2 maps were acquired to extract the IVD relaxation times, and the complete manual segmentation of IVDs at all levels was performed using Display® software. Lumbar IVD T2 relaxation times showed significant correlation with PT (P<0.01), GT (P<0.01), TPA (P<0.01), PI-LL (P=0.01), and LLL (P=0.01). No difference was noted between Roussouly subtypes regarding T2 relaxation times at any disc level. Data from questionnaires showed no correlation with T2 relaxation times. Global tilt and T1 pelvic angle were correlated with IVD composition changes (T2 relaxometry). There was no correlation between clinical symptoms and IVD T2 relaxation times.

Low back pain; Intervertebral disc degeneration; Magnetic resonance imaging; Spine; Radiography


Introduction

Degeneration of the intervertebral discs (IVD) has a high prevalence and is associated with low back pain (11. Jarvik JG, Deyo RA. Diagnostic evaluation of low back pain with emphasis on imaging. Ann Intern Med 2002; 137: 586-597, doi: 10.7326/0003-4819-137-7-200210010-00010.
https://doi.org/10.7326/0003-4819-137-7-...
-2. Brinjikji W, Diehn FE, Jarvik JG, Carr CM, Kallmes DF, Murad MH, et al. MRI findings of disc degeneration are more prevalent in adults with low back pain than in asymptomatic controls: a systematic review and meta-analysis. AJNR Am J Neuroradiol 2015; 36: 2394-2399, doi: 10.3174/ajnr.A4498.
https://doi.org/10.3174/ajnr.A4498...
33. Hancock M, Maher C, Macaskill P, Latimer J, Kos W, Pik J. MRI findings are more common in selected patients with acute low back pain than controls? Eur Spine J 2012; 21: 240-246, doi: 10.1007/s00586-011-1955-7.
https://doi.org/10.1007/s00586-011-1955-...
). The evaluation of disc degeneration based on the Pfirrmann classification, which uses the T2-weighted magnetic resonance sequence, although useful, is subjective and limited in the assessment of early onset of IVD degeneration (44. Pfirrmann CW, Metzdorf A, Elfering A, Hodler J, Boos N. Effect of aging and degeneration on disc volume and shape: a quantitative study in asymptomatic volunteers. J Orthop Res 2006; 24: 1086-1094, doi: 10.1002/jor.20113.
https://doi.org/10.1002/jor.20113...
). The T2 relaxometry MRI technique is, otherwise, capable of quantifying and providing information on early biochemical changes in the IVD degeneration (55. Weidenbaum M, Foster RJ, Best BA, Saed-Nejad F, Nickoloff E, Newhouse J, et al. Correlating magnetic resonance imaging with the biochemical content of the normal human intervertebral disc. J Orthop Res 1992; 10: 552-561, doi: 10.1002/jor.1100100410.
https://doi.org/10.1002/jor.1100100410...
-6. Perry J, Haughton V, Anderson PA, Wu Y, Fine J, Mistretta C. The value of T2 relaxation times to characterize lumbar intervertebral disks: preliminary results. AJNR Am J Neuroradiol 2006; 27: 337-342. 7. Chiu EJ, Newitt DC, Segal MR, Hu SS, Lotz JC, Majumdar S. Magnetic resonance imaging measurement of relaxation and water diffusion in the human lumbar intervertebral disc under compression in vitro. Spine (Phila Pa 1976) 2001; 26: E437-E444, doi: 10.1097/00007632-200110010-00017.
https://doi.org/10.1097/00007632-2001100...
8. Blumenkrantz G, Zuo J, Li X, Kornak J, Link TM, Majumdar S. In vivo 3.0-tesla magnetic resonance T1rho and T2 relaxation mapping in subjects with intervertebral disc degeneration and clinical symptoms. Magn Reson Med 2010; 63: 1193-1200, doi: 10.1002/mrm.22362.
https://doi.org/10.1002/mrm.22362...
99. Marinelli NL, Haughton VM, Anderson PA. T2 relaxation times correlated with stage of lumbar intervertebral disk degeneration and patient age. AJNR Am J Neuroradiol 2010; 31: 1278-1282, doi: 10.3174/ajnr.A2080.
https://doi.org/10.3174/ajnr.A2080...
). Recent studies have established the relationship between T2 relaxometry and semiquantitative Pfirrmann classification, and the decrease in T2 relaxation times correlated significantly with increased disc degeneration (1010. Welsch GH, Trattnig S, Paternostro-Sluga T, Bohndorf K, Goed S, Stelzeneder D, et al. Parametric T2 and T2* mapping techniques to visualize intervertebral disc degeneration in patients with low back pain: initial results on the clinical use of 3.0 Tesla MRI. Skeletal Radiol 2011; 40: 543-551, doi: 10.1007/s00256-010-1036-8.
https://doi.org/10.1007/s00256-010-1036-...
-11. Trattnig S, Stelzeneder D, Goed S, Reissegger M, Mamisch TC, Paternostro-Sluga T, et al. Lumbar intervertebral disc abnormalities: comparison of quantitative T2 mapping with conventional MR at 3.0 T. Eur Radiol 2010; 20: 2715-2722. 1212. Detiger SE, Holewijn RM, Hoogendoorn RJ, van Royen BJ, Helder MN, Berger FH, et al. MRI T2* mapping correlates with biochemistry and histology in intervertebral disc degeneration in a large animal model. Eur Spine J 2015; 24: 1935-1943, doi: 10.1007/s00586-014-3498-1.
https://doi.org/10.1007/s00586-014-3498-...
). In addition, the relationship between sagittal balance and degeneration of the IVD using the semiquantitative Pfirrmann was demonstrated (1313. Torrie PA, McKay G, Byrne R, Morris SA, Harding IJ. The influence of lumbar spinal subtype on lumbar intervertebral disc degeneration in young and middle-aged adults. Spine Deform 2015; 3: 172-179, doi: 10.1016/j.jspd.2014.08.006.
https://doi.org/10.1016/j.jspd.2014.08.0...
-14. Menezes-Reis R, Bonugli GP, Dalto VF, da Silva Herrero CF, Defino HL, Nogueira-Barbosa MH. Association between lumbar spine sagittal alignment and L4-L5 disc degeneration among asymptomatic young adults. Spine (Phila Pa 1976) 2016; 41: E1081-E1087, doi: 10.1097/BRS.0000000000001568.
https://doi.org/10.1097/BRS.000000000000...
15. Mardare M, Oprea M, Popa I, Zazgyva A, Niculescu M, Poenaru DV. Sagittal balance parameters correlate with spinal conformational type and MRI changes in lumbar degenerative disc disease: results of a retrospective study. Eur J Orthop Surg Traumatol 2016; 26: 735-743, doi: 10.1007/s00590-016-1842-3.
https://doi.org/10.1007/s00590-016-1842-...
1616. Ogon I, Takashima H, Morita T, Oshigiri T, Terashima Y, Yoshimoto M, et al. Association between spinopelvic alignment and lumbar intervertebral disc degeneration quantified with magnetic resonance imaging t2 mapping in patients with chronic low back pain. Spine Surg Relat Res 2019; 4: 135-141, doi: 10.22603/ssrr.2019-0051.
https://doi.org/10.22603/ssrr.2019-0051...
). However, we did not identify studies in the literature that evaluated the potential correlation between sagittal alignment and the biochemical composition of the IVD evaluated by quantitative magnetic resonance techniques. Given that previous studies revealed a negative correlation between T2 values and disc degeneration (66. Perry J, Haughton V, Anderson PA, Wu Y, Fine J, Mistretta C. The value of T2 relaxation times to characterize lumbar intervertebral disks: preliminary results. AJNR Am J Neuroradiol 2006; 27: 337-342.-7. Chiu EJ, Newitt DC, Segal MR, Hu SS, Lotz JC, Majumdar S. Magnetic resonance imaging measurement of relaxation and water diffusion in the human lumbar intervertebral disc under compression in vitro. Spine (Phila Pa 1976) 2001; 26: E437-E444, doi: 10.1097/00007632-200110010-00017.
https://doi.org/10.1097/00007632-2001100...
88. Blumenkrantz G, Zuo J, Li X, Kornak J, Link TM, Majumdar S. In vivo 3.0-tesla magnetic resonance T1rho and T2 relaxation mapping in subjects with intervertebral disc degeneration and clinical symptoms. Magn Reson Med 2010; 63: 1193-1200, doi: 10.1002/mrm.22362.
https://doi.org/10.1002/mrm.22362...
), we expected that patients with sagittal imbalance assessed through spinopelvic parameters would have more severe disc degeneration and consequently lower T2 values. Our hypothesis was that spinopelvic parameters indicating a sagittal imbalance are related with disc degeneration quantified by the T2 relaxometry technique. The purpose of this study was to verify the correlation between quantitative T2 relaxation measurements of IVD with spinopelvic parameters and clinical symptoms in patients with chronic low back pain.

Material and Methods

Patients

Ninety-one consecutive patients with chronic lumbar pain were included in this cross-sectional and prospective study. There were 56 female patients (mean age 53.5 years; range 23-76 years) and 35 male patients (mean age 53.6 years; range 19-73 years). Inclusion criteria were age greater than or equal to 18 years and chronic low back pain (visual analogue scale ≥3) with recurrent episodes in the last 6 months.

Exclusion criteria consisted of scoliosis or spondylolisthesis, vertebral fracture, neoplasia, postoperative status, disabling hip disease, infection, or bone metabolic diseases. Of the 107 evaluated patients, 23 were excluded for the following reasons: 10 individuals with spondylolisthesis, 5 individuals with scoliosis, and 1 individual with vertebral fracture. All individuals were assessed by the Oswestry disability index (ODI) and visual analogue scale (VAS). The study was approved by the Clinical Hospital of the Ribeirao Preto Medical School (USP) ethics committee, and all research was performed in accordance with the current regulations. Informed consent was obtained from all participants of this study.

Spinopelvic parameters evaluation

Each patient underwent lateral panoramic radiography acquired with a CR Long Length Vertical Imaging System (Kodak Direct View, Carestream Health, USA). Subjects stood with arms supported on a stand, shoulders at 30° of flexion, and elbows slightly flexed according to previous literature (1717. Marks M, Stanford C, Newton P. Which lateral radiographic positioning technique provides the most reliable and functional representation of a patient's sagittal balance? Spine (Phila Pa 1976) 2009; 34: 949-954, doi: 10.1097/BRS.0b013e318199650a.
https://doi.org/10.1097/BRS.0b013e318199...
) to minimize possible postural compensations.

Surgimap® software (Nemaris Inc., USA, version 2.2.9.6) was used to measure spinopelvic parameters and vertebral curvature angles. The following parameters were evaluated (Figure 1): pelvic incidence (PI), pelvic-tilt (PT), sacral slope (SS), sagittal vertical axis (SVA), global tilt (GT), T1 pelvic angle (TPA), lumbar lordosis (LL), thoracic kyphosis (TK), pelvic incidence minus lumbar lordosis mismatch (PI-LL), and lack of lumbar lordosis (LLL). We used the formula proposed by Schwab and co-workers (1818. Schwab F, Patel A, Ungar B, Farcy JP, Lafage V. Adult spinal deformity-postoperative standing imbalance: how much can you tolerate? An overview of key parameters in assessing alignment and planning corrective surgery. Spine (Phila Pa 1976) 2010; 35: 2224-2231, doi: 10.1097/BRS.0b013e3181ee6bd4.
https://doi.org/10.1097/BRS.0b013e3181ee...
) to determine ideal LL, where LL = PI + 9. Therefore, using measured LL and ideal theoretical lordosis, we can establish the LLL of each patient. The contours of the femoral heads were marked, and lines were drawn adjacent to the sacral cartilage, at the upper plateau of the L1 vertebra, the lower plateau of T12, the upper plateau of T4, the upper plateau of T1, and the center of C7 and C2 vertebral bodies. From these markings, the software automatically calculated the spinopelvic parameters and the vertebral curvatures.

Figure 1
Panoramic radiograph of a 65-year-old male patient illustrating the spinopelvic parameter measurements: A, sacral slope (SS); B, pelvic incidence (PI); C, pelvic tilt (PT); D, thoracic kyphosis (TK); E, lumbar lordosis (LL); F, sagittal vertical axis (SVA); G, global tilt (GT); H, T1 pelvic angle (TPA).

Disc composition evaluation

We evaluated lumbar IVD composition using a 1.5T MRI scanner (Achieva; Philips Healthcare, Best, The Netherlands) with a 16-channel spinal coil. To avoid possible physiologic daily variations in disc water content, MRIs were acquired in the afternoon. All images of the lumbar spine were acquired in the sagittal plane using the following parameters: field of view=22×22 cm, thickness=4 mm, number of slices=16, and matrix=220×217. A 2D fast spin-echo T2-weighted sequence, echo time (TE)=120 ms, and repetition time (TR)=3900 ms were used for anatomical reference and IVD segmentation. We used a T2 relaxometry multi-echo sequence as follows: TE=20/40/60/80/100/120/140/160 ms and TR=3900 ms.

We generated T2 relaxometry maps using MINC tools and Display® software (McConell Brain Imaging Centre, Canada) to extract the IVD relaxation times. T2 maps were computed on a pixel-by-pixel basis using an exponential decay model, where S0 is the equilibrium magnetization signal and S(TE) is the signal acquired with the echo time (TE):

S(TE) = S) * exp TET2

Using T2-weighted images, the complete volume of each lumbar IVD was manually segmented (Figure 2). The segmentation encompassed both the nucleus pulposus and annulus fibrosus using all the sagittal images in which the five discs were identified, being careful not to include the subchondral bone.

Figure 2
Image obtained from a 25-year-old woman with low back pain. Graphical representation of segmentation, encompassing the nucleus pulposus and annulus fibrosus of each lumbar intervertebral disc, with overlapping relaxometry maps. The schematic color scale represents the T2 relaxation times in milliseconds (ms).

Image evaluation

The first examiner, a radiologist with four years of experience, was responsible for spinopelvic measurements and IVD segmentations. The second examiner, a researcher with four years of experience in spinal imaging who was blind to the measurements of the first examiner performed the same measurements to assess interobserver reliability. Participants were categorized into study groups by the two examiners according to the Roussouly classification for the 4 postural subtypes (1919. Roussouly P, Gollogly S, Berthonnaud E, Dimnet J. Classification of the normal variation in the sagittal alignment of the human lumbar spine and pelvis in the standing position. Spine (Phila Pa 1976) 2005; 30: 346-353, doi: 10.1097/01.brs.0000152379.54463.65.
https://doi.org/10.1097/01.brs.000015237...
). When disagreement occurred, the opinion of a third examiner, a senior researcher and radiologist with more than 17 years of experience in musculoskeletal radiology, was considered. The distribution of the Roussouly subgroups was as follows: Type 1, n=13; Type 2, n=23; Type 3, n=44; and Type 4, n=11.

Statistical analysis

Student's t-test was used to compare the composition of the IVD according to gender. We performed a mixed effects linear regression between the spinopelvic parameters PI, PT, SS, SVA, GT, TPA, LL, TK, PI-LL, LLL and the T2 relaxation times, controlling for the effect of age. Roussouly subtypes and T2 relaxation times were compared with analysis of variance (ANOVA). The correlation between questionnaire data and T2 relaxation time was analyzed using the Spearman correlation coefficient (ρ). The interobserver agreement was analyzed by intraclass correlation coefficient with a 95% confidence interval, considering zero to 0.19 as poor agreement, 0.20-0.39 as fair agreement, 0.40-0.59 as moderate agreement, 0.60-0.79 as substantial agreement, and 0.80-1.00 as almost perfect agreement. All statistical analysis was performed with a statistical software (SAS 9.2; SAS Institute, USA). For this study, a significance level of 5% was adopted. The statistical power (1-β error probability) of our sample (91 subjects) was 0.86 with an α error of 5%.

Results

Table 1 presents the values of T2 relaxation times of the five IVD levels as means and standard deviations. The L2-L3 and L3-L4 discs presented the highest values for T2 relaxation times, and these values decreased as they approached the upper (L1-L2) and lower (L4-L5 and L5-S1) ends.

Table 1
T2 relaxation times of patients with chronic low back pain for each disc level.

The mean T2 relaxation time was T2=91.3±8.24 ms for men and 91.83±9.41 ms for women. No statistically significant differences in T2 relaxometry values were noted between genders at any disc level.

The spinopelvic parameters of our study group are provided in Table 2. Age was negatively correlated with T2 relaxation time in all segments (L1-L2: R: -0.49, P<0.01; L2-L3: R: -0.48, P<0.01; L3-L4: -0.26, P<0.01; L4-L5: -0.28, P<0.01; L5-S1: -0.20, P<0.01). A negative relationship was noted between T2 relaxometry values and spinopelvic parameters PT (beta=-0.26, P<0.01), GT (beta=0.25, P=0.01), TPA (beta=-0.31, P<0.01), PI-LL (beta=-0.18, P=0.01), and LLL (beta=-0.18, P=0.01) (Table 3).

Table 2
Spinopelvic parameters of patients with chronic low back pain.
Table 3
Correlation between spinopelvic parameters and T2 relaxation times of lumbar discs, controlling for the effect of age.

No difference was noted between Roussouly subtypes regarding T2 relaxation times at any disc level.

The mean ODI was 23.86±9.41 points (minimum: 5, maximum: 40, median: 25) and the mean VAS value was 7.48±1.88 (minimum: 2, maximum: 10, median: 8). No statistically significant correlation was noted between T2 relaxation values and questionnaire data at any disc level.

The ICC for T2 relaxometry values demonstrated almost perfect inter-observer agreement for manual disc segmentations (ICC: 0.96) and for spinopelvic parameters (ICC: 0.85 to 0.96) (Table 4).

Table 4
Intraclass correlation coefficients (ICC) and 95% confidence intervals (CI) of manual disc segmentation and spinopelvic parameter measurements.

Discussion

This study assessed the correlation between quantitative T2 relaxation measurements of lumbar intervertebral discs and spinopelvic parameters reports using a consecutive series of 91 patients. We observed a significant negative relationship between T2 relaxometry values and the spinopelvic parameters PT, GT, TPA, PI-LL, and LLL. Roussouly subtype did not correlate with T2 relaxation time. We did not observe a significant correlation between mean T2 relaxometry values at the IVDs and questionnaire data.

Patients with lumbar degenerative disease are characterized by anterior sagittal imbalance, loss of LL, and increased PT (1212. Detiger SE, Holewijn RM, Hoogendoorn RJ, van Royen BJ, Helder MN, Berger FH, et al. MRI T2* mapping correlates with biochemistry and histology in intervertebral disc degeneration in a large animal model. Eur Spine J 2015; 24: 1935-1943, doi: 10.1007/s00586-014-3498-1.
https://doi.org/10.1007/s00586-014-3498-...
,1515. Mardare M, Oprea M, Popa I, Zazgyva A, Niculescu M, Poenaru DV. Sagittal balance parameters correlate with spinal conformational type and MRI changes in lumbar degenerative disc disease: results of a retrospective study. Eur J Orthop Surg Traumatol 2016; 26: 735-743, doi: 10.1007/s00590-016-1842-3.
https://doi.org/10.1007/s00590-016-1842-...
,2020. Barrey C, Jund J, Noseda O, Roussouly P. Sagittal balance of the pelvis-spine complex and lumbar degenerative diseases. A comparative study about 85 cases. Eur Spine J 2007; 16: 1459-1467, doi: 10.1007/s00586-006-0294-6.
https://doi.org/10.1007/s00586-006-0294-...
,2121. Savarese LG, Menezes-Reis R, Bonugli GP, Herrero CFPD, Defino HLA, Nogueira-Barbosa MH. Spinopelvic sagittal balance: what does the radiologist need to know? Radiol Bras 2020; 53: 175-184, doi: 10.1590/0100-3984.2019.0048.
https://doi.org/10.1590/0100-3984.2019.0...
). Ogon et al. (1616. Ogon I, Takashima H, Morita T, Oshigiri T, Terashima Y, Yoshimoto M, et al. Association between spinopelvic alignment and lumbar intervertebral disc degeneration quantified with magnetic resonance imaging t2 mapping in patients with chronic low back pain. Spine Surg Relat Res 2019; 4: 135-141, doi: 10.22603/ssrr.2019-0051.
https://doi.org/10.22603/ssrr.2019-0051...
) found an association between anterior annulus fibrosus degeneration in all lumbar disc levels and hypolordosis of the lumbar spine, anterior translation of the trunk and posterior inclination of the pelvis in chronic low back pain using quantitative MRI. SVA and PT are spinopelvic parameters that indicate the severity of adult spine deformity, but some points should be considered.

First, the SVA measurement can be reduced by postural compensatory mechanisms, such as pelvic retroversion. Therefore, a high PT may “hide” a larger spinal deformity when only the SVA is considered. SVA and PT are related and the magnitude of one affects the other. In patients who may have developed pelvic retroversion (high PT) to compensate for an underlying spinal malalignment and maintain their head over the pelvis (low SVA), the evaluation of the sagittal alignment by SVA alone would not detect their sagittal malalignment. Lafage et al. suggested that PT should be considered together with SVA to detect patients with spinal deformity in the sagittal plane without high SVA due to pelvic compensation (2222. Lafage V, Schwab F, Patel A, Hawkinson N, Farcy JP. Pelvic tilt and truncal inclination: two key radiographic parameters in the setting of adults with spinal deformity. Spine (Phila Pa 1976) 2009; 34: E599-E606, doi: 10.1097/BRS.0b013e3181aad219.
https://doi.org/10.1097/BRS.0b013e3181aa...
). Of note, a successful realignment plan should not only restore the spinopelvic relationship but also “zero out” compensatory mechanisms, which drain energy and affect patients’ quality of life. Our results without relationship between SVA and T2 relaxometry reinforces that SVA should not be assessed alone in the evaluation of sagittal plane. In contrast, TPA and GT have advantages in assessing global alignment given that they account for pelvic retroversion and trunk anteversion and are not affected by postural or radiographic calibrations (2323. Protopsaltis T, Schwab F, Bronsard N, Smith JS, Klineberg E, Mundis G, et al. TheT1 pelvic angle, a novel radiographic measure of global sagittal deformity, accounts for both spinal inclination and pelvic tilt and correlates with health-related quality of life. J Bone Joint Surg Am 2014; 96: 1631-1640, doi: 10.2106/JBJS.M.01459.
https://doi.org/10.2106/JBJS.M.01459...
,2424. Obeid I, Hauger O, Aunoble S, Bourghli A, Pellet N, Vital JM. Global analysis of sagittal spinal alignment in major deformities: correlation between lack of lumbar lordosis and flexion of the knee. Eur Spine J 2011; 20: 681-685, doi: 10.1007/s00586-011-1936-x.
https://doi.org/10.1007/s00586-011-1936-...
). Furthermore, they are strongly correlated with SVA, PT, and PI-LL (2525. Banno T, Togawa D, Arima H, Hasegawa T, Yamato Y, Kobayashi S, et al. The cohort study for the determination of reference values for spinopelvic parameters (T1 pelvic angle and global tilt) in elderly volunteers. Eur Spine J 2016; 25: 3687-3693, doi: 10.1007/s00586-016-4411-x.
https://doi.org/10.1007/s00586-016-4411-...
) and showed relation with disc degeneration quantitatively assessed by T2 relaxometry.

The purpose of sagittal misalignment correction surgeries is to achieve a good global sagittal alignment and restore the LL, considering the morphology of the pelvis.

We observed no relationship between LL and disc degeneration quantified by T2 relaxometry. However, we identified the individual's LLL calculated by the difference between the observed lordosis and the ideal lordosis estimated by the Schwab formula. We also identified a relationship between the PI-LL, a valuable tool in the intraoperative planning (2626. Legaye J, Duval-Beaupàre G, Hecquet J, Marty C. Pelvic incidence: a fundamental pelvic parameter for three-dimensional regulation of spinal sagittal curves. Eur Spine J 1998; 7: 99-103, doi: 10.1007/s005860050038.
https://doi.org/10.1007/s005860050038...
,2727. Schwab F, Lafage V, Boyce R, Skalli W, Farcy JP. Gravity line analysis in adult volunteers: age-related correlation with spinal parameters, pelvic parameters, and foot position. Spine (Phila Pa 1976) 2006; 31: E959-E967, doi: 10.1097/01.brs.0000248126.96737.0f.
https://doi.org/10.1097/01.brs.000024812...
), and T2 relaxometry values. These findings are expected when we consider the relationship between pelvic morphology, represented by PI, and LL. As postulated, LL has a strong correlation with PI; a high PI is accompanied by a high LL and a low PI is accompanied by a low LL (2828. Morvan G, Mathieu P, Vuillemin V, Guerini H, Bossard P, Zeitoun F, et al. Standardized way for imaging of the sagittal spinal balance. Eur Spine J 2011; 20: 602-608, doi: 10.1007/s00586-011-1927-y.
https://doi.org/10.1007/s00586-011-1927-...
). Thus, when two patients have the same LL and different PIs, the LL can be normal or reduced, depending on the PI value, and the calculation of the LLL provides this information. Therefore, individuals with a greater mismatch between pelvic incidence and lumbar lordosis (PI-LL) and higher LLL estimated by the difference between the expected lordosis and the observed lordosis, exhibited lower T2 relaxometry values in our study, indicating discs with more severe degeneration.

Consistent with other studies (2929. Niu G, Yang J, Wang R, Dang S, Wu EX, Guo Y. MR imaging assessment of lumbar intervertebral disk degeneration and age-related changes: apparent diffusion coefficient versus T2 quantitation. AJNR Am J Neuroradiol 2011; 32: 1617-1623, doi: 10.3174/ajnr.A2556.
https://doi.org/10.3174/ajnr.A2556...
-30. Wu N, Liu H, Chen J, Zhao L, Zuo W, Ming Y, et al. Comparison of apparent diffusion coefficient and T2 relaxation time variation patterns in assessment of age and disc level related intervertebral disc changes. PLoS One 2013; 8: e69052, doi: 10.1371/journal.pone.0069052.
https://doi.org/10.1371/journal.pone.006...
3131. Menezes-Reis R, Salmon CE, Bonugli GP, Mazoroski D, Tamashiro MH, Savarese LG, et al. Lumbar intervertebral discs T2 relaxometry and T1ρ relaxometry correlation with age in asymptomatic young adults. Quant Imaging Med Surg 2016; 6: 402-412, doi: 10.21037/qims.2016.08.01.
https://doi.org/10.21037/qims.2016.08.01...
), we observed a negative linear correlation between age and T2 values. As older individuals tend to have more dehydrated discs and therefore lower T2 relaxometry values, we used a mixed effects statistical model that included age as a fixed effect. By fixing the age, we were able to evaluate the influence of spinopelvic parameters on disc degeneration without the bias of age. With that, we identified a significant negative correlation between the parameters PT, GT, TPA, PI-LL, LLL, and T2 values.

Patients with different Roussouly subtypes did not present significant differences in T2 relaxation times of IVD in our study. Of note, the classification of Roussouly et al. (3232. Roussouly P, Pinheiro-Franco JL. Biomechanical analysis of the spino-pelvic organization and adaptation in pathology. Eur Spine J 2011; 20: 609-618, doi: 10.1007/s00586-011-1928-x.
https://doi.org/10.1007/s00586-011-1928-...
) was originally described in 160 asymptomatic individuals aged between 18 and 48 years. This classification system has its limitations given that individuals with low SS also tend to have low LL if they remain sagittally balanced in relation to the lumbar segment of the spine. Torrie et al. observed no significant correlation between postural Roussouly subtypes and semiquantitative Pfirrmann grading of disc degeneration (1313. Torrie PA, McKay G, Byrne R, Morris SA, Harding IJ. The influence of lumbar spinal subtype on lumbar intervertebral disc degeneration in young and middle-aged adults. Spine Deform 2015; 3: 172-179, doi: 10.1016/j.jspd.2014.08.006.
https://doi.org/10.1016/j.jspd.2014.08.0...
). Consistent with this finding, our results reinforced that there is no correlation between Roussouly subtypes and disc degeneration in patients with chronic low back pain. A possible explanation is the fact that this classification is based on the SS and LL parameters, which exhibited no correlation with the T2 relaxometry values in our study.

We observed higher T2 relaxation values in the L2-L3 and L3-L4 discs and these values decreased in the upper (L1-L2) and lower (L4-L5 and L5-S1) lumbar spine extremities. It is reasonable to speculate that the L1-L2, L4-L5, and L5-S1 levels offer greater movement in the sagittal plane compared with other levels, given that L1-L2 is closer to the thoracolumbar junction. In addition, L4-L5 and L5-S1 support the largest amount of load due to gravitational force.

The T2 values found in the IVD were very close to those reported in the study by Niu et al. (2929. Niu G, Yang J, Wang R, Dang S, Wu EX, Guo Y. MR imaging assessment of lumbar intervertebral disk degeneration and age-related changes: apparent diffusion coefficient versus T2 quantitation. AJNR Am J Neuroradiol 2011; 32: 1617-1623, doi: 10.3174/ajnr.A2556.
https://doi.org/10.3174/ajnr.A2556...
), who used ROIs that covered the nucleus pulposus and annulus fibrosus in the most central sagittal slice, and the study by Menezes-Reis et al. (3333. Menezes-Reis R, Salmon CE, Carvalho CS, Bonugli GP, Chung CB, Nogueira-Barbosa MH. T1ρ and T2 mapping of the intervertebral disk: comparison of different methods of segmentation. AJNR Am J Neuroradiol 2015; 36: 606-611, doi: 10.3174/ajnr.A4125.
https://doi.org/10.3174/ajnr.A4125...
), who explored the segmentation of the vertebral disc throughout its extension. Most of the studies in the literature preferred to use the segmentation of NP and AF individually (1010. Welsch GH, Trattnig S, Paternostro-Sluga T, Bohndorf K, Goed S, Stelzeneder D, et al. Parametric T2 and T2* mapping techniques to visualize intervertebral disc degeneration in patients with low back pain: initial results on the clinical use of 3.0 Tesla MRI. Skeletal Radiol 2011; 40: 543-551, doi: 10.1007/s00256-010-1036-8.
https://doi.org/10.1007/s00256-010-1036-...
,3030. Wu N, Liu H, Chen J, Zhao L, Zuo W, Ming Y, et al. Comparison of apparent diffusion coefficient and T2 relaxation time variation patterns in assessment of age and disc level related intervertebral disc changes. PLoS One 2013; 8: e69052, doi: 10.1371/journal.pone.0069052.
https://doi.org/10.1371/journal.pone.006...
,3434. Stelzeneder D, Welsch GH, Kovács BK, Goed S, Paternostro-Sluga T, Vlychou M, et al. Quantitative T2 evaluation at 3.0T compared to morphological grading of the lumbar intervertebral disc: a standardized evaluation approach in patients with low back pain. Eur J Radiol 2012; 81: 324-330, doi: 10.1016/j.ejrad.2010.12.093.
https://doi.org/10.1016/j.ejrad.2010.12....
,3535. Wang YX, Zhao F, Griffith JF, Mok GS, Leung JC, Ahuja AT, et al. T1rho and T2 relaxation times for lumbar disc degeneration: an in vivo comparative study at 3.0-Tesla MRI. Eur Radiol 2013; 23: 228-234, doi: 10.1007/s00330-012-2591-2.
https://doi.org/10.1007/s00330-012-2591-...
). Moreover, for practical reasons, most previous studies used small geometric ROIs for the evaluation of degenerated discs (1111. Trattnig S, Stelzeneder D, Goed S, Reissegger M, Mamisch TC, Paternostro-Sluga T, et al. Lumbar intervertebral disc abnormalities: comparison of quantitative T2 mapping with conventional MR at 3.0 T. Eur Radiol 2010; 20: 2715-2722.,3434. Stelzeneder D, Welsch GH, Kovács BK, Goed S, Paternostro-Sluga T, Vlychou M, et al. Quantitative T2 evaluation at 3.0T compared to morphological grading of the lumbar intervertebral disc: a standardized evaluation approach in patients with low back pain. Eur J Radiol 2012; 81: 324-330, doi: 10.1016/j.ejrad.2010.12.093.
https://doi.org/10.1016/j.ejrad.2010.12....
). Intervertebral discs classified as Pfirrmann I and II can be better distinguished when they are segmented separately. As degeneration progresses to Pfirrmann III and IV degrees, this distinction becomes difficult or impossible (3636. Pfirrmann CW, Metzdorf A, Zanetti M, Hodler J, Boos N. Magnetic resonance classification of lumbar intervertebral disc degeneration. Spine (Phila Pa 1976) 2001; 26: 1873-1878, doi: 10.1097/00007632-200109010-00011.
https://doi.org/10.1097/00007632-2001090...
). Thus, a relative advantage considered in the choice of the segmentation method of this study was that the segmentation of the entire area or the complete volume of the disc allows a safer comparison between discs with different degrees of degeneration.

A study performed by Menezes-Reis et al. (3737. Menezes-Reis R, Garrido Salmon CE, Bonugli GP, Mazoroski D, Savarese LG, Herrero CFPS, et al. Association between spinal alignment and biochemical composition of lumbar intervertebral discs assessed by quantitative magnetic resonance imaging. Quant Imaging Med Surg 2021; 11: 2428-2441, doi: 10.21037/qims-20-40.
https://doi.org/10.21037/qims-20-40...
) with asymptomatic volunteers did not find an association between global spinopelvic parameters and IVD composition. Because their sample consisted of asymptomatic volunteers and the spine posture in these subjects could be considered within normal limits, a lack of association could be explained by the lack of individuals without sagittal plane imbalance.

There were some limitations in our study. First, our study used a cross-sectional design, and thus it is impossible to infer cause and effect relationships between the variables. Moreover, the study group consisted only of individuals with chronic low back pain, and asymptomatic individuals were not evaluated. Although the radiographic technique did not allow direct evaluation of lower limb compensation, systematic control was performed during the acquisition to prevent small compensatory mechanisms.

In conclusion, several spinopelvic parameters were correlated with disc degeneration quantitatively assessed using the T2 relaxometry technique. Given that this study involved the pathophysiology of mechanical association with disc composition, we believe that longitudinal studies with long follow-up are necessary to evaluate the potential future use of T2 relaxometry in clinical practice. To our knowledge, this is the first study correlating spinopelvic parameters GT and TPA with disc degeneration evaluated by T2 relaxometry.

Acknowledgments

This study was supported by National Council of Technological and Scientific Development (CNPq, Grant No. 473688/2012-3), Coordination of Superior Level Staff Improvement (CAPES, Scholarship Grant No. 31806584875 and Finance Code 001), Financing of Studies and Projects of Brazil (FINEP, Ref. 0184/07), and Foundation of Support to Teaching, Research and Assistance (FAEPA 2044/2017 and 728/2002).

References

  • 1
    Jarvik JG, Deyo RA. Diagnostic evaluation of low back pain with emphasis on imaging. Ann Intern Med 2002; 137: 586-597, doi: 10.7326/0003-4819-137-7-200210010-00010.
    » https://doi.org/10.7326/0003-4819-137-7-200210010-00010
  • 2
    Brinjikji W, Diehn FE, Jarvik JG, Carr CM, Kallmes DF, Murad MH, et al. MRI findings of disc degeneration are more prevalent in adults with low back pain than in asymptomatic controls: a systematic review and meta-analysis. AJNR Am J Neuroradiol 2015; 36: 2394-2399, doi: 10.3174/ajnr.A4498.
    » https://doi.org/10.3174/ajnr.A4498
  • 3
    Hancock M, Maher C, Macaskill P, Latimer J, Kos W, Pik J. MRI findings are more common in selected patients with acute low back pain than controls? Eur Spine J 2012; 21: 240-246, doi: 10.1007/s00586-011-1955-7.
    » https://doi.org/10.1007/s00586-011-1955-7
  • 4
    Pfirrmann CW, Metzdorf A, Elfering A, Hodler J, Boos N. Effect of aging and degeneration on disc volume and shape: a quantitative study in asymptomatic volunteers. J Orthop Res 2006; 24: 1086-1094, doi: 10.1002/jor.20113.
    » https://doi.org/10.1002/jor.20113
  • 5
    Weidenbaum M, Foster RJ, Best BA, Saed-Nejad F, Nickoloff E, Newhouse J, et al. Correlating magnetic resonance imaging with the biochemical content of the normal human intervertebral disc. J Orthop Res 1992; 10: 552-561, doi: 10.1002/jor.1100100410.
    » https://doi.org/10.1002/jor.1100100410
  • 6
    Perry J, Haughton V, Anderson PA, Wu Y, Fine J, Mistretta C. The value of T2 relaxation times to characterize lumbar intervertebral disks: preliminary results. AJNR Am J Neuroradiol 2006; 27: 337-342.
  • 7
    Chiu EJ, Newitt DC, Segal MR, Hu SS, Lotz JC, Majumdar S. Magnetic resonance imaging measurement of relaxation and water diffusion in the human lumbar intervertebral disc under compression in vitro Spine (Phila Pa 1976) 2001; 26: E437-E444, doi: 10.1097/00007632-200110010-00017.
    » https://doi.org/10.1097/00007632-200110010-00017
  • 8
    Blumenkrantz G, Zuo J, Li X, Kornak J, Link TM, Majumdar S. In vivo 3.0-tesla magnetic resonance T1rho and T2 relaxation mapping in subjects with intervertebral disc degeneration and clinical symptoms. Magn Reson Med 2010; 63: 1193-1200, doi: 10.1002/mrm.22362.
    » https://doi.org/10.1002/mrm.22362
  • 9
    Marinelli NL, Haughton VM, Anderson PA. T2 relaxation times correlated with stage of lumbar intervertebral disk degeneration and patient age. AJNR Am J Neuroradiol 2010; 31: 1278-1282, doi: 10.3174/ajnr.A2080.
    » https://doi.org/10.3174/ajnr.A2080
  • 10
    Welsch GH, Trattnig S, Paternostro-Sluga T, Bohndorf K, Goed S, Stelzeneder D, et al. Parametric T2 and T2* mapping techniques to visualize intervertebral disc degeneration in patients with low back pain: initial results on the clinical use of 3.0 Tesla MRI. Skeletal Radiol 2011; 40: 543-551, doi: 10.1007/s00256-010-1036-8.
    » https://doi.org/10.1007/s00256-010-1036-8
  • 11
    Trattnig S, Stelzeneder D, Goed S, Reissegger M, Mamisch TC, Paternostro-Sluga T, et al. Lumbar intervertebral disc abnormalities: comparison of quantitative T2 mapping with conventional MR at 3.0 T. Eur Radiol 2010; 20: 2715-2722.
  • 12
    Detiger SE, Holewijn RM, Hoogendoorn RJ, van Royen BJ, Helder MN, Berger FH, et al. MRI T2* mapping correlates with biochemistry and histology in intervertebral disc degeneration in a large animal model. Eur Spine J 2015; 24: 1935-1943, doi: 10.1007/s00586-014-3498-1.
    » https://doi.org/10.1007/s00586-014-3498-1
  • 13
    Torrie PA, McKay G, Byrne R, Morris SA, Harding IJ. The influence of lumbar spinal subtype on lumbar intervertebral disc degeneration in young and middle-aged adults. Spine Deform 2015; 3: 172-179, doi: 10.1016/j.jspd.2014.08.006.
    » https://doi.org/10.1016/j.jspd.2014.08.006
  • 14
    Menezes-Reis R, Bonugli GP, Dalto VF, da Silva Herrero CF, Defino HL, Nogueira-Barbosa MH. Association between lumbar spine sagittal alignment and L4-L5 disc degeneration among asymptomatic young adults. Spine (Phila Pa 1976) 2016; 41: E1081-E1087, doi: 10.1097/BRS.0000000000001568.
    » https://doi.org/10.1097/BRS.0000000000001568
  • 15
    Mardare M, Oprea M, Popa I, Zazgyva A, Niculescu M, Poenaru DV. Sagittal balance parameters correlate with spinal conformational type and MRI changes in lumbar degenerative disc disease: results of a retrospective study. Eur J Orthop Surg Traumatol 2016; 26: 735-743, doi: 10.1007/s00590-016-1842-3.
    » https://doi.org/10.1007/s00590-016-1842-3
  • 16
    Ogon I, Takashima H, Morita T, Oshigiri T, Terashima Y, Yoshimoto M, et al. Association between spinopelvic alignment and lumbar intervertebral disc degeneration quantified with magnetic resonance imaging t2 mapping in patients with chronic low back pain. Spine Surg Relat Res 2019; 4: 135-141, doi: 10.22603/ssrr.2019-0051.
    » https://doi.org/10.22603/ssrr.2019-0051
  • 17
    Marks M, Stanford C, Newton P. Which lateral radiographic positioning technique provides the most reliable and functional representation of a patient's sagittal balance? Spine (Phila Pa 1976) 2009; 34: 949-954, doi: 10.1097/BRS.0b013e318199650a.
    » https://doi.org/10.1097/BRS.0b013e318199650a
  • 18
    Schwab F, Patel A, Ungar B, Farcy JP, Lafage V. Adult spinal deformity-postoperative standing imbalance: how much can you tolerate? An overview of key parameters in assessing alignment and planning corrective surgery. Spine (Phila Pa 1976) 2010; 35: 2224-2231, doi: 10.1097/BRS.0b013e3181ee6bd4.
    » https://doi.org/10.1097/BRS.0b013e3181ee6bd4
  • 19
    Roussouly P, Gollogly S, Berthonnaud E, Dimnet J. Classification of the normal variation in the sagittal alignment of the human lumbar spine and pelvis in the standing position. Spine (Phila Pa 1976) 2005; 30: 346-353, doi: 10.1097/01.brs.0000152379.54463.65.
    » https://doi.org/10.1097/01.brs.0000152379.54463.65
  • 20
    Barrey C, Jund J, Noseda O, Roussouly P. Sagittal balance of the pelvis-spine complex and lumbar degenerative diseases. A comparative study about 85 cases. Eur Spine J 2007; 16: 1459-1467, doi: 10.1007/s00586-006-0294-6.
    » https://doi.org/10.1007/s00586-006-0294-6
  • 21
    Savarese LG, Menezes-Reis R, Bonugli GP, Herrero CFPD, Defino HLA, Nogueira-Barbosa MH. Spinopelvic sagittal balance: what does the radiologist need to know? Radiol Bras 2020; 53: 175-184, doi: 10.1590/0100-3984.2019.0048.
    » https://doi.org/10.1590/0100-3984.2019.0048
  • 22
    Lafage V, Schwab F, Patel A, Hawkinson N, Farcy JP. Pelvic tilt and truncal inclination: two key radiographic parameters in the setting of adults with spinal deformity. Spine (Phila Pa 1976) 2009; 34: E599-E606, doi: 10.1097/BRS.0b013e3181aad219.
    » https://doi.org/10.1097/BRS.0b013e3181aad219
  • 23
    Protopsaltis T, Schwab F, Bronsard N, Smith JS, Klineberg E, Mundis G, et al. TheT1 pelvic angle, a novel radiographic measure of global sagittal deformity, accounts for both spinal inclination and pelvic tilt and correlates with health-related quality of life. J Bone Joint Surg Am 2014; 96: 1631-1640, doi: 10.2106/JBJS.M.01459.
    » https://doi.org/10.2106/JBJS.M.01459
  • 24
    Obeid I, Hauger O, Aunoble S, Bourghli A, Pellet N, Vital JM. Global analysis of sagittal spinal alignment in major deformities: correlation between lack of lumbar lordosis and flexion of the knee. Eur Spine J 2011; 20: 681-685, doi: 10.1007/s00586-011-1936-x.
    » https://doi.org/10.1007/s00586-011-1936-x
  • 25
    Banno T, Togawa D, Arima H, Hasegawa T, Yamato Y, Kobayashi S, et al. The cohort study for the determination of reference values for spinopelvic parameters (T1 pelvic angle and global tilt) in elderly volunteers. Eur Spine J 2016; 25: 3687-3693, doi: 10.1007/s00586-016-4411-x.
    » https://doi.org/10.1007/s00586-016-4411-x
  • 26
    Legaye J, Duval-Beaupàre G, Hecquet J, Marty C. Pelvic incidence: a fundamental pelvic parameter for three-dimensional regulation of spinal sagittal curves. Eur Spine J 1998; 7: 99-103, doi: 10.1007/s005860050038.
    » https://doi.org/10.1007/s005860050038
  • 27
    Schwab F, Lafage V, Boyce R, Skalli W, Farcy JP. Gravity line analysis in adult volunteers: age-related correlation with spinal parameters, pelvic parameters, and foot position. Spine (Phila Pa 1976) 2006; 31: E959-E967, doi: 10.1097/01.brs.0000248126.96737.0f.
    » https://doi.org/10.1097/01.brs.0000248126.96737.0f
  • 28
    Morvan G, Mathieu P, Vuillemin V, Guerini H, Bossard P, Zeitoun F, et al. Standardized way for imaging of the sagittal spinal balance. Eur Spine J 2011; 20: 602-608, doi: 10.1007/s00586-011-1927-y.
    » https://doi.org/10.1007/s00586-011-1927-y
  • 29
    Niu G, Yang J, Wang R, Dang S, Wu EX, Guo Y. MR imaging assessment of lumbar intervertebral disk degeneration and age-related changes: apparent diffusion coefficient versus T2 quantitation. AJNR Am J Neuroradiol 2011; 32: 1617-1623, doi: 10.3174/ajnr.A2556.
    » https://doi.org/10.3174/ajnr.A2556
  • 30
    Wu N, Liu H, Chen J, Zhao L, Zuo W, Ming Y, et al. Comparison of apparent diffusion coefficient and T2 relaxation time variation patterns in assessment of age and disc level related intervertebral disc changes. PLoS One 2013; 8: e69052, doi: 10.1371/journal.pone.0069052.
    » https://doi.org/10.1371/journal.pone.0069052
  • 31
    Menezes-Reis R, Salmon CE, Bonugli GP, Mazoroski D, Tamashiro MH, Savarese LG, et al. Lumbar intervertebral discs T2 relaxometry and T1ρ relaxometry correlation with age in asymptomatic young adults. Quant Imaging Med Surg 2016; 6: 402-412, doi: 10.21037/qims.2016.08.01.
    » https://doi.org/10.21037/qims.2016.08.01
  • 32
    Roussouly P, Pinheiro-Franco JL. Biomechanical analysis of the spino-pelvic organization and adaptation in pathology. Eur Spine J 2011; 20: 609-618, doi: 10.1007/s00586-011-1928-x.
    » https://doi.org/10.1007/s00586-011-1928-x
  • 33
    Menezes-Reis R, Salmon CE, Carvalho CS, Bonugli GP, Chung CB, Nogueira-Barbosa MH. T1ρ and T2 mapping of the intervertebral disk: comparison of different methods of segmentation. AJNR Am J Neuroradiol 2015; 36: 606-611, doi: 10.3174/ajnr.A4125.
    » https://doi.org/10.3174/ajnr.A4125
  • 34
    Stelzeneder D, Welsch GH, Kovács BK, Goed S, Paternostro-Sluga T, Vlychou M, et al. Quantitative T2 evaluation at 3.0T compared to morphological grading of the lumbar intervertebral disc: a standardized evaluation approach in patients with low back pain. Eur J Radiol 2012; 81: 324-330, doi: 10.1016/j.ejrad.2010.12.093.
    » https://doi.org/10.1016/j.ejrad.2010.12.093
  • 35
    Wang YX, Zhao F, Griffith JF, Mok GS, Leung JC, Ahuja AT, et al. T1rho and T2 relaxation times for lumbar disc degeneration: an in vivo comparative study at 3.0-Tesla MRI. Eur Radiol 2013; 23: 228-234, doi: 10.1007/s00330-012-2591-2.
    » https://doi.org/10.1007/s00330-012-2591-2
  • 36
    Pfirrmann CW, Metzdorf A, Zanetti M, Hodler J, Boos N. Magnetic resonance classification of lumbar intervertebral disc degeneration. Spine (Phila Pa 1976) 2001; 26: 1873-1878, doi: 10.1097/00007632-200109010-00011.
    » https://doi.org/10.1097/00007632-200109010-00011
  • 37
    Menezes-Reis R, Garrido Salmon CE, Bonugli GP, Mazoroski D, Savarese LG, Herrero CFPS, et al. Association between spinal alignment and biochemical composition of lumbar intervertebral discs assessed by quantitative magnetic resonance imaging. Quant Imaging Med Surg 2021; 11: 2428-2441, doi: 10.21037/qims-20-40.
    » https://doi.org/10.21037/qims-20-40

Publication Dates

  • Publication in this collection
    11 Nov 2022
  • Date of issue
    2022

History

  • Received
    7 Apr 2022
  • Accepted
    17 Oct 2022
Associação Brasileira de Divulgação Científica Av. Bandeirantes, 3900, 14049-900 Ribeirão Preto SP Brazil, Tel. / Fax: +55 16 3315-9120 - Ribeirão Preto - SP - Brazil
E-mail: bjournal@terra.com.br