Original research

Ultrasound imaging measures of vertebral bony landmark distances are weakly to moderately correlated with intervertebral disc height as assessed by MRI

Abstract

Objectives To assess the validity and reliability of ultrasound-derived interbony landmark distances as a proxy for MRI-derived intervertebral disc (IVD) height.

Methods This is a cross-sectional criterion validity study. Twelve college-aged participants without current low back pain completed both MRI and ultrasound imaging of the lumbar spine in a prone position. Single-segment and multisegment distances between the spinous and mammillary processes at the lumbar segments (L2/L3, L3/L4, L4/L5) were measured twice using ultrasound and analysed digitally. Sagittal slices of the lumbar spine were taken via T1-weighted MRI and IVD height, and the overall distance between IVDs L2/L3 and L4/L5 was imaged once and measured twice.

Results There was moderate correlation between multilevel-based measurements (overall distance between L2 and L5, r=0.677, p=0.016) and the average across three levels (r=0.596, p=0.041) when using the spinous processes as bony landmarks. Single-segment measures were not significantly correlated (all: p>0.092). Accuracy and precision were better for the overall MRI-derived distance between the three IVDs from L2 and L5 MRI and the distance measured between the spinous processes L2–L5. There was excellent reliability within multiple measurements at each location, with intraclass correlation coefficient, ICC(3,1), ranging from 0.93 to 0.99 (95% CI 0.82 to 0.99) for ultrasound and from 0.98 to 0.99 (95% CI 0.92 to 0.99) for MRI.

Conclusion Findings do not support the use of ultrasound imaging for estimating single-segment IVD height, yet it may be used to measure the change in distance over time with a certain degree of precision based on its excellent reliability.

Key messages

What is already known on this topic

  • Changes in lumbar spine height are related to water fluctuations which drastically change with ageing and disease processes.

  • Ultrasound imaging has been used to measure lumbar spine height, but the results have not been validated against data obtained from MRI, or the data were not collected with the participant in a prone position.

What this study adds

  • There is some validity in using ultrasound-derived multisegment but not single-segment distances to estimate related intervertebral disc height.

How this study might affect research, practice or policy

  • The assessment of spinal height changes over time could become an additional tool for clinicians who use point-of-care ultrasound imaging.

Introduction

The spine height of a fully grown individual decreases with load-bearing and increases with unloading.1 For example, during space flight, a spinal height increase of 2 inches is common,2 which is largely due to water fluctuations occurring in the intervertebral disc (IVD), and places the astronaut at an increased risk of IVD trauma.3 Measuring IVD height and assessing changes thereof can provide valuable information on the hydration status of the IVD.2 Decreases in IVD height occur with ageing and disease processes.4 Together with other radiological signs, such as the presence of osteophytes and endplate sclerosis, decreased IVD height is considered a sign of disc degeneration.5 In addition, loss of disc height has been correlated with low back pain,6 and therefore knowledge of this phenomenon could be of value for the treating clinician. Point-of-care ultrasound imaging is becoming more popular as a form of examination. It is used with ‘the intent of clarifying uncertain clinical examination findings to enhance the quality and effectiveness of a physical therapy intervention’.7 The results of the systematic review published in 20198 endorse the addition of point-of-care ultrasound imaging to the assessment and treatment process.

Various methods have been employed to study immediate changes in IVD height, such as MRI, stadiometry, X-ray and ultrasound imaging.1 9–12 MRI and X-rays can measure IVD height directly, while stadiometry and ultrasound imaging estimate IVD height indirectly.13 MRI is considered the reference standard when imaging the spine,14 yet is often not available or obtainable for different reasons. X-ray, while also yielding high-resolution images, exposes patients to potentially harmful ionising radiation.15 Ultrasound imaging does not provide exposure to ionising radiation and has the additional advantages of being financially cost-effective and portable. However, the method yields comparatively lower-resolution images and may be less accurate. Nevertheless, ultrasound imaging has been used to measure spinal segment height in vivo1 2 16–18 and in vitro.13 Distances between different bony landmarks, including lumbar transverse processes,17 mammillary processes,1 anterior vertebral bodies2 and spinous processes,16 18 have been used. Still, the results were either not validated against data obtained from MRI,1 2 17 18 or the data were not collected with the participant in prone position.1 2 16 Ultrasound imaging of the lower back cannot be easily done with the patient in supine position due to inaccessibility of the spine and associated probe placement issues.

The specific aims of this study were to (1) assess the validity of estimating IVD height from ultrasound images by measuring single-segment and multisegment distances between spinous processes and mammillary processes and compare them with single-segment and multisegment MRI-derived IVD height, respectively; and (2) assess intrarater test–retest reliability and intrarater reliability of measuring distances between posterior bony landmarks using ultrasound imaging and IVD heights on MRI scans. We hypothesised that intermammillary distances are better proxies for IVD heights than interspinous processes distances.

Methods

This is a cross-sectional criterion validity study. Participants were recruited from the university community and by word of mouth. We included individuals aged 18 years and older without current, that is, within the previous week, low back pain. Informed consent was obtained before data collection, and the rights of the participants were protected. Each participant was informed that their data would be submitted for publication and that their confidentiality would be protected. The exclusion criteria were (1) inability to lie on stomach for 30 min; (2) history of spinal surgery; (3) history of traumatic injury to the spine; (4) known scoliosis for which prior medical consultation was sought; (5) known claustrophobia; and (6) absolute contraindications to MRI.19

Imaging

Ultrasound and MRI were performed in a prone lying position to directly compare the measurements and keep the spine in a consistent position across imaging modalities. Due to the ferromagnetic characteristics of the ultrasound machine, the participant had to be moved out of the MRI room and repositioned for the ultrasound imaging procedure. A folded-up small pillow was placed under the hips to decrease lordosis of the lumbar spine. A 3-inch high rolled-up towel also supported the forehead, and the arms were resting at the side.

Ultrasound imaging

Ultrasound images (Logiq e, Probe, GE Healthcare, Chicago, Illinois, USA) were captured by placing the probe (12 MHz linear) in longitudinal orientation on the participant’s skin with a water-based gel as a coupling agent between the probe and the skin. The frequency ranged from 8 MHz to 12 MHz as needed to enhance imaging of the bony structure. The spinous processes were first visualised. Separate images were taken of each of the following pairs of spinous processes: L2 and L3, L3 and L4, as well as L4 and L5. Keeping the probe in longitudinal orientation, it was moved about 2 cm laterally off midline to locate the posterior-most hyperechoic points, typically indicating the mammillary processes of L5–L2. Using the LogiqView mode, a panoramic image was captured. The probe was placed distally with the right edge of the sacrum in the field of view, the orientation marker on the probe pointed cephalad. Slowly the probe was glided cranially along with the mammillary processes until L2 was visualised. This was repeated for the mammillary processes on the contralateral side. Images of the spinous processes and intermammillary processes were recorded for processing after the assessment. This imaging process was performed twice for each side, with the probe briefly removed from the skin between each scan without repositioning the participant. No specific breathing instructions were given to the participant for the scan because ultrasound imaging is less prone to artefacts with restful breathing compared with MRI.20

The same examiner, who has 1 year of experience with this imaging modality, collected all images and performed the postimaging measurements using the OsiriX (Pixmeo, Geneva, Switzerland) software program in a blinded fashion. The distances between the mammillary processes were measured by first drawing straight vertical lines through the apex (most posterior aspect) of each mammillary process (L2/L3, L3/L4 and L4/L5). Horizontal lines were drawn between each vertical line (see figure 1). Interspinous process distance was measured by drawing vertical lines through the apex of the spinous processes with a horizontal connecting line between the two (see figure 2). In rare cases, there were two apices on one lumbar spinous process. In these cases, we chose the most posterior one.

Figure 1
Figure 1

Straight vertical lines were drawn through the mammillary processes, signifying the apex of each facet joint (L2/L3, L3/L4 and L4/L5). Horizontal lines were then drawn between neighbouring vertical lines and measured. *Indicates the apex of the hyperechoic bone cortex, the location of the mammillary processes.

Figure 2
Figure 2

Vertical lines were drawn through the apices of the spinous processes of vertebrae L2, L3, L4 and L5. Horizontal lines were then drawn between neighbouring vertical lines and measured.

Magnetic resonance imaging

MRI was performed using the 3-Tesla Siemens MAGNETOM Tim Trio (Siemens Healthineers, Erlangen Germany) scanner with a four-channel flexible coil. The participant was placed prone on the imaging table. The inferior border of the 51.2 cm × 22 cm coil was positioned caudally on the posterior superior iliac spines and placed in a cranial direction along the long axis of the spine. The participant was instructed to breathe shallowly during the scan. A total of 36 sagittal slices were acquired using a fast, low-angle shot sequence using the following parameters: field of view 192×173×119 mm, repetition time 13 ms, echo time 2 ms, slices 36, flip angle 10°, resolution 1.5×1.5 mm on the plane, 3 mm through the plane and no fat suppression. The scan took approximately 38 s. After obtaining images in Digital Imaging and Communications in Medicine (DICOM) format, they were viewed on OsiriX (Pixmeo). A midsagittal plane view was used and the IVDs of interest (L2/L3, L3/L4 and L4/L5) were manually segmented from the two-dimensional T1-weighted structural scan images. Lines connecting the superior anterior and inferior anterior borders and the superior posterior and inferior posterior borders were drawn, representing the anterior and posterior height of the IVD. A horizontal line was drawn connecting the anterior and posterior lines. An additional horizontal line was drawn to mark the halfway point of the full IVD width. We then drew a line perpendicular to the vertebral body at the halfway mark of the IVD.21 An additional vertical line was drawn from the superior midpoint at IVD segment L2/L3 to the inferior midpoint at IVD segment L4/L5 to capture the multisegment distance that includes three IVDs (figure 3).

Figure 3
Figure 3

Anterior/posterior and mid-IVD height at segments L2/L3, L3/L4 and L4/L5 as well as the distance between the superior border of L2/L3 to the inferior border of L4/L5 were measured. The latter was done by drawing a straight line between the cephalad midpoint of IVD at segment L2/L3 and the caudal midpoint of IVD at segment L4/L5. IVD, intervertebral disc.

Statistical analyses were conducted using Stata V.16 statistical software. For assessing measurement precision using four ultrasound examinations per participant, it is recommended that at least nine participants are measured22; hence, our study was capable of detecting reasonable estimates of reliability. The strength and direction of association between MRI and ultrasound-derived IVD height were assessed by Pearson correlation coefficient (weak: 0.10–0.39; moderate: 0.40–0.69; strong: >0.7).23 Univariate linear regression was used to quantify the proportion of variance of MRI-derived IVD height explained by ultrasound-derived methods. Bland-Altman plots were generated to assess agreement between MRI and ultrasound imaging methods and to graphically depict the ultrasound–MRI difference as a function of the mean for each subject.24 Analyses used average IVD height at segments L2/L3, L3/L4 and L4/L5 (ie, the average of the anterior, mid and posterior height), the average from across the three segments, and the total distance from the superior border of L2/L3 to the inferior border of L4/L5 (L2–L5) as obtained from the MRI. An alpha level of 0.05 was adopted for statistical tests. Calculation of intraclass correlation coefficient (ICC) was performed using SPSS V.28 for consistency with random subjects and fixed raters, mixed model ICC (model3,k). Intrarater test–retest reliability, that is, the same rater measures the first and second image, of the ultrasound and intrarater reliability, that is, the rater measures the same image twice, of MRI measurements were obtained by computing a two-way mixed-model ICC(3,1) with absolute agreement. To identify testing error we calculated the SE of the measurement (SEm) and per cent SEm (%SEm) and the minimum detectable difference at 95% (MDD at 95%) using the following equations: SEm=SD (Sq rt 1–ricc); 95% CI SEm=muscle mean±(1.96×SEm); MDD=SEm×1.96×√2.

Patient and public involvement

Patients and/or the public were not involved in the design, or conduct, or reporting or dissemination plans of this research.

Results

Twelve asymptomatic participants (8 men: mean±SD age 24.1±1.7 years, height 182.0±7.9 cm, weight 84.6±23.6 kg; 4 women: age 23.5±2.2 years, height 171.0±6.3 cm, weight 63.0±8.0 kg) completed both MRI and ultrasound imaging of the lumbar spine within several minutes of each other.

The average distances for both ultrasound and MRI imaging are presented in table 1.

Table 1
|
Ultrasound and MRI measurements

Comparison between methods

The total distance between the L2 and L5 spinous processes (ultrasound) moderately correlated with MRI-derived distance from the superior border of IVD at the L2/L3 level to the inferior border of IVD at the L4/L5 level (r=0.677, p=0.016). IVD height, when measured indirectly via spinous process distance using ultrasound imaging, moderately correlated with MRI-derived IVD height when averaged across the three IVDs examined (r=0.596, p=0.041), yet not at each individual IVD (all, p≥0.092; table 2).

Table 2
|
Correlations between MRI and ultrasound methods

Linear regression revealed that the distance between the L2 and L5 spinous processes explained approximately 46% (R2=0.458) of the variance in MRI-derived distance from the superior border of L2/L3 IVD to the inferior border of L4/L5 IVD. Moreover, the average interspinous distance explained approximately 36% (R2=0.355) of the variance in average MRI IVD height. These values were greater than the variance explained by mammillary-based methods (all: <2.7%).

Bland-Altman analysis revealed a mean bias of −0.159 cm and wide limits of agreement (LoA) ranging from −1.288 to 0.970 for the overall distance between L2 and L5 measured by MRI and ultrasound imaging using spinous processes as bony reference (online supplemental material figure 1). The mean bias for the averages (average IVD height derived from levels L2/L3, L3/L4 and L4/L5 and average distances between spinous processes L2, L3, L4 and L5) was −2.123 cm, and the LoA were narrower, ranging from −2.518 cm to −1.729 cm.

Within-method repeatability

The methods of measuring the distance between the spinous processes and the mammillary processes and IVD height were reliable for ultrasound imaging and MRI, respectively (table 3).

Table 3
|
Distances between spinal segments and reliability indicators

Discussion

The results of our study indicate that there is some validity in using ultrasound-derived multisegment bony distances (between mammillary processes and between spinous processes) as a measure of multisegment IVD height. Single-segment measurements, however, showed little validity. More specifically, we observed significant moderate correlations (r=0.60–0.68) for the multivertebral segment-based measurements (L2–L5) and the average across these three segments, while single-segment measurements were not significantly correlated with their respective MRI-derived IVD height. This suggests that averaging or combining distances likely helps account for variations between individual vertebral segments.

Our hypothesis that intermammillary distances are better proxies for IVD heights than interspinous processes distances was not confirmed. Considering that any distance of the same landmark measured between adjacent vertebrae includes a part or the whole height of the vertebral body, we recognise that this problem is the same for any bony landmark. However, measuring the distance between two different bony landmarks in neighbouring vertebrae, for example the mammillary process of the inferior vertebra to the transverse process of the superior vertebra, could potentially lessen this problem.

Bland-Altman

To further understand the relationship between ultrasound-derived and MRI-derived methods we used to estimate lumbar IVD height, we also employed Bland-Altman plots. The precision of the measurements (ie, the narrower the LoA, the more precise the measurements) and degree of agreement between the two methods can be gleaned from the Bland-Altman plots. The LoA in our study ranged from −1.288 cm to 0.970 cm for MRI-derived overall distance between three IVDs (L2–L5) compared with the ultrasound-derived distance between spinous processes L2–L5. The LoA for the average MRI-derived disc height and the average ultrasound-derived spinous process distance were smaller (range: −2.518 cm to −1.729 cm) and should therefore be considered more precise. The former LoA were not narrow enough to conclude that the method of indirectly measuring IVD height was precise. The latter indicated that the ultrasound imaging measurements were within 0.4 cm of the MRI-derived values. We suggest that the 0.4 cm difference between ultrasound imaging measurements and MRI-derived IVD height could be explained by the respective anatomical starting and stopping points of the two measurement locations: the MRI measurements started at the inferior border of the vertebral body of L2 and ended at the superior border of L5, while the ultrasound measurement began at the apex of the L2 spinous process, which corresponds more closely to the middle of the body of L2 and ended at the apex of L5 spinous process, which again corresponds more closely to the middle of the vertebral body. Thus, based on the differences in the measured structures, it is not surprising that the distances measured are dissimilar, although consistently so.

The degree of agreement between the ultrasound and MRI measures can be extrapolated from the bias shown in the Bland-Altman plot (online supplemental material figure 1; only the significant correlations are shown). The x-axes on the Bland-Altman plots represent the total (online supplemental material figure 1A) and average (online supplemental material figure 1B) distances between L2 and L5 as measured by ultrasound (spinous processes) and by MRI (three IVDs and two vertebral bodies) (distance=(MRI distance+ultrasound distance)/2) for each participant. The y-axes on the Bland-Altman plots represent the absolute differences (total and average distance difference, respectively) between the two values (distance difference=MRI distance–ultrasound distance) for each participant. The mean bias for the data representing the overall distance between segments L2 and L5 was −0.159 cm, while it was −2.123 cm for the average distance. Since the mean biases are smaller than 0, it indicates that the ultrasound-derived measures are consistently larger than the MRI measures.

The positive correlation between IVD height and the discrepancy between ultrasound and MRI measurements seen in both Bland-Altman graphs indicate that the discrepancy between MRI and ultrasound imaging values becomes larger with increasing multisegment distance. This could be a function of increasing IVD or/and increasing vertebral body heights. Regardless, clinically this implies that ultrasound measurements tend to underestimate the overall distance (and overestimate to a lesser degree the average distance) between three IVDs in people with shorter lumbar spines and overestimate it in people with longer lumbar spines.

Linear regression

The linear regression analysis sheds further light on the relationship between the ultrasound-derived and MRI-derived methods we used to estimate lumbar IVD height. Linear regression revealed that the ultrasound-derived multisegment distance between the L2 and L5 spinous processes accounted for about 46% of the variance of the data obtained via MRI. This finding supports a relationship between the two measurement methods and modalities, although a moderate one. However, we believe that this result is more a function of the correlation between lumbar vertebral body height when measured with ultrasound imaging and compared with our reference standard obtained via MRI. This is based on the fact that, on average, vertebral bodies make up more than 60% of the distance from segments L2–L5.25

Reliability and repeatability

The second aim of this study was to assess intrarater test–retest reliability and intrarater reliability. Our methods, ultrasound imaging (captured twice, measured once) and MRI (captured once, measured twice), showed excellent reliability, with ICC between 0.93 and 1.00. Interestingly, the ICC for the single-segment measurements of intermammillary distances was lower than the multisegment measurements; this was reversed for measurements using interspinous process distances (both using ultrasound imaging).

Minimum detectable difference

The MDD is a useful clinical variable that indicates the amount of change needed to have a significant difference in a given population. Our results for MDD at 95% (table 3) suggest that the ultrasound imaging protocol can detect a single-segment distance change of 2–3 mm, depending on the lumbar segment, when using the mammillary processes as bony landmarks. When using the spinous processes as bony landmarks, this value is even smaller, ranging between 0.7 mm and 1.5 mm. We submit that these values are small enough to suggest this measuring tool is reliable. To confirm this proposition we also assessed the possible measuring error and relative measuring error, as indicated by the %SEm. The values for our %SEm range between 2.6 and 2.9. This indicates that we can expect a measuring error of up to 2.9% or 0.7 mm on a measured distance of 2.7 cm (eg, between the L4/L5 mammillary processes), which is considered to be good.26 Our results lead us to conclude that our ultrasound imaging protocol is reliable in making multisegment measurements, although not to estimate IVD height. Of note, our results apply to asymptomatic individuals rather than those with low back pain, so generalisation should be made with caution.

Limitations

Despite our effort to design a carefully controlled study, it also had several potential limitations. Our sample size is relatively small but larger than the sample size of other studies.16 18 27 Our population sample included only healthy young adults who have a relatively small likelihood of having IVD or vertebral body pathology. In the case of presence of corner osteophytes and Schmorl’s nodes, the MRI-derived IVD height measurement protocol would have to be adjusted. We only assessed the validity and reliability of segments L2/L3, L3/L4 and L4/L5. We did not consider any segments above and below. While ultrasound imaging cannot be performed in a supine position, it is the preferred position for MRI. The lordotic angle will differ between these two positions and thus affect the anterior and posterior IVD dimensions (as measured on MRI) and the interbony distances (as measured via ultrasound imaging). Lastly, selection of MRI images may not correspond to image selection in ultrasound imaging. This is most pertinent to our interspinous process distance measurements since the MRI-derived and ultrasound-derived images should be collected in the midsagittal plane.

Conclusions

This study demonstrates that there is some validity in using ultrasound-derived multisegment bony distances as a measure of multisegment IVD height. However, single-segment intermammillary and interspinous process distances are not good estimates of the actual respective single-segment lumbar IVD height. Nevertheless, because the ultrasound-derived data demonstrate excellent intrarater test–retest reliability, these measurements could be investigated for their ability to assess changes in lumbar spine height over time. This is especially likely for the multisegment measurements using either measurement method because their measurement repeatability is excellent as indicated by their high ICC, low SEm, low %SEm and low MDD. The assessment of spinal height changes over time could therefore become an additional tool for clinicians who use point-of-care ultrasound imaging.