Abstract
BACKGROUND AND PURPOSE: TSE-based T2-weighted imaging of the spine has long scan times. This work proposes a fast imaging protocol using variable refocusing flip angles, optimized for blurring and specific absorption rate.
MATERIALS AND METHODS: A variable refocusing flip angle echo-train was optimized for the spine to improve the point spread function and minimize the specific absorption rate, yielding images with improved spatial resolution and SNR compared with the constant flip angle sequence. Data were acquired from 51 patients (35 lumbar, 16 whole-spine) using conventional TSE and the proposed sequence, with a single-shot variant for whole-spine. Noninferiority analysis was performed to evaluate the efficiency of the proposed technique.
RESULTS: The proposed multishot sequence resulted in a 2× shorter scan time with a >1.5× lower specific absorption rate. The variable flip angle sequence was noninferior to the conventional TSE (P < .025) for all image-quality and clinical criteria except signal-to-noise ratio for the lumbar spine protocol. However, mean image scores for the TSE-variable refocusing flip angle were ≥4.3 for all criteria, and concordance analysis showed high agreement (>90%) with the TSE, indicating clinical equivalence. The single-shot sequence resulted in 4× shorter whole-spine scans, and image scores were ≥4.4 for all criteria, attesting to its clinical utility.
CONCLUSIONS: We present a fast T2-weighted spine protocol using variable refocusing flip angles, including a single-shot variant. The sequences have better point spread function behavior than their constant flip angle counterparts and, being faster, should be less sensitive to patient motion, often seen in the longer TSE scans.
ABBREVIATIONS:
- ETL
- echo-train length
- SAR
- specific absorption rate
- VFA
- variable refocusing flip angle
T2-weighted MR imaging protocols are routinely used in the clinic for spine imaging. T2WI sequences such as 2D fast spin-echo or turbo spin-echo have proved useful for the detection and diagnosis of osteomyelitis, stenosis, nerve root compression, disc characterization, and other pathology.1,2 A typical spine protocol includes sequences for T2WI, T1WI, and short-tau inversion recovery fat-suppressed T2 imaging, all with FSE/TSE readouts. While TSE/FSE with multiecho readout reduces scan times compared with conventional spin-echo sequences, the scan times are still on the order of 3–4 minutes, and these double or triple when performing combination spine examinations (eg, thoracolumbar). This scan time is uncomfortable for patients with painful spine conditions and makes the imaging more vulnerable to motion artifacts. While approaches using parallel imaging3 have been explored to reduce scan times, the use of higher acceleration factors is limited by image-quality degradation arising from poor coil coverage. While the scan time of a multislice TSE sequence could also be reduced by increasing the echo-train length (ETL), the corresponding increase in the specific absorption rate (SAR) results in an increased minimum TR, which offsets any reduction in scan time. In addition, increased ETL exacerbates image blurring from T2 decay during the long echo-train.
Single-shot sequences such as single-shot fast spin-echo and HASTE have been used for lumbar spine myelography,4 cervical spine motion imaging,5 and functional imaging of the spine.6 However, these sequences are not preferred for routine T2WI because they have severe image blurring and SAR limitations arising from the long echo-trains, often addressed by scanning at a low spatial resolution.7
Various approaches have been proposed to minimize radiofrequency power deposition in TSE-based sequences.8⇓–10 These approaches optimize the refocusing flip angle train to improve the imaging point spread function, leading to reduced T2 blurring and an improved signal-to-noise ratio and SAR performance. While the variable refocusing flip angle (VFA) technique has been successfully used for 3D TSE-based brain,11,12 knee,13,14 and spine15 imaging, its use in 2D imaging has been rather limited. Recently, VFA along with 2D single-shot TSE was used in abdominal and pelvic imaging16 with increased sharpness and significant reduction in scan times due to the reduced SAR compared with conventional single-shot TSE.
In this work, we propose a fast T2-weighted spine imaging protocol using a multishot turbo spin-echo sequence with variable refocusing flip angles (TSE-VFA) optimized for spatial resolution, SNR, and SAR. We also explore the use of a single-shot variant (HASTE-VFA) as a fast T2WI screening sequence for whole-spine protocols, in which the use of conventional TSE is prohibitively long.
Materials and Methods
Pulse Sequence and Parameter Optimization
The vendor-supplied TSE and HASTE sequences were modified to incorporate the refocusing flip angle modulation scheme proposed by Busse et al,10 which parameterizes the refocusing flip angle train by using 4 control angles (αstart, αmin, αcent, αend). Following the approach of Loening et al,16 flip angles were optimized to maximize the SNR for the desired effective TE and minimize the SAR, with minimal spatial blurring. Numeric simulations were performed using Matlab (MathWorks, Natick, Massachusetts) to narrow the range of “optimal” refocusing flip angle values. For each possible value of the control angle vector (αstart, αmin, αcent, αend), T2 signal evolution from the turbo spin-echo pulse sequence was simulated using the Echo Phase Graph algorithm, assuming T1 and T2 values for the spine at 3T (T1 = 1060 ms, T2 = 69 ms), echo-train length = 56, and effective TE = 105 ms. The imaging point spread function was computed using the simulated signal. The peak and full width at half maximum value of the point spread function were measured because these are indicative of the SNR and resolution, respectively. Relative SAR was computed as the sum of squared refocusing flip angles divided by the echo-train length: , where αi is flip angle of the ith refocusing pulse.
The optimal flip angles were chosen by using a 2-step search algorithm: 1) The search space was reduced by choosing control angles that maximize the point spread function on the basis of a preset threshold while maintaining the SNR at the effective TE. 2) Control angles that yield minimum SAR were chosen to generate the tailored refocusing flip angle train. This approach ensures adequate resolution and SNR for the VFA sequence compared with a constant flip angle train for the same echo-train length. To minimize motion-related signal dephasing at very low values of αmin17, we chose a value of αmin = 45°. Due to the nonconvex nature of the search space, the resultant flip angle parameters are one of many possible optimal solutions.
Phantom and Volunteer Scanning
TSE-VFA and HASTE-VFA sequences were implemented and tested on a 3T Siemens Skyra scanner (Siemens, Erlangen, Germany). Phantom experiments were performed to study the blurring behavior and SNR of the proposed technique based on flip angle parameters chosen from simulations. Two nickel-doped agarose gel phantoms with T2 = 45 ms, T1 = 1000 ms and T2 = 80 ms, T1 = 1300 ms were imaged using the following sequences: a conventional TSE sequence at ETL = 21, TSE at ETL = 56, and TSE-VFA at ETL = 56. Data were acquired at an effective TE = 100 ms for all 3 sequences. The mean signal in the 2 phantoms was measured by drawing ROIs inside the object and in the background noise. The acquisition SNR was computed as . In addition, SNR efficiency was calculated as . The use of variable refocusing flip angles can alter TSE image contrast. To ensure comparable contrast behavior of TSE-VFA, we also measured relative contrast between the 2 phantoms: . The SAR computed by the scanner was also recorded for the 3 sequences.
Because a complete sweep of the flip angle parameter space is not practical in vivo, the optimal flip angles obtained from simulations were validated using data acquired from volunteers. Sagittal lumbar spine data were acquired from 5 volunteers after informed consent, using the conventional TSE and TSE-VFA sequences with the parameters listed in Table 1.
Clinical Imaging
The optimized TSE-VFA sequence was added to the lumbar and cervical spine protocols as an addition to the routine T2WI TSE sequences. The HASTE-VFA sequence was added to the whole-spine protocol (covering the cervical, thoracic, and lumbar spine), which currently omits the T2WI TSE sequence due to long scan times (∼9 minutes). This omission was possible due to the lower spatial resolution typically used in whole-spine screening protocols, which makes single-shot VFA with an adequate SNR viable. The effective TE was set to 108 ms for all the sequences, and images were acquired in the sagittal plane in accordance with the standard clinical protocol. The scan parameters for conventional TSE, TSE-VFA, and HASTE-VFA are shown in Table 1. As part of a prospective study, data were acquired from 51 patients (35 lumbar spine, 16 whole-spine) after written informed consent and in compliance with the institutional review board.
Image-Quality Assessment
Performance of the proposed VFA sequences was assessed independently by 2 neuroradiologists (2 and 12 years of experience) in a blinded fashion for both the volunteer and clinical data. To avoid recall bias, images from the 2 sequences were scored at least 2 weeks apart with random assignment of one or the other sequence for each patient in each reading session. Images were graded on a scale of 1–5 (1, nondiagnostic images; 2, severely limited; 3, limited; four, 1 or 2 suboptimal attributes but still diagnostic; and 5, optimal image quality for diagnosis). Images were assessed on the basis of the following criteria: edge sharpness, motion, artifacts, and noise. Clinical utility of the sequences was assessed by quantifying the ability of the reader to interpret facet joints, endplates, nerve roots, spinal cord, and discs, also on a scale of 1–5. All levels of the spine were evaluated in terms of the image-quality metrics and the degree to which each structure could be interpreted. The lowest score was then applied for each of the anatomic structures.
Quantitative Analysis
An ROI analysis was performed on the volunteer and clinical lumbar spine data to quantitatively measure the SNR and relative contrast of the 2 sequences. ROIs were placed in the L3 vertebral body and the L3–L4 disc to measure the mean signal intensity. The noise SD was measured from the background air region. The vertebral body SNR was calculated as , and the SNR efficiency was defined as SNR per unit scan time computed as . The relative contrast between the disc and the vertebral body was calculated to quantify image contrast, defined as . A 2-tailed t test was performed to compare the mean SNR, relative contrast, and SAR from the 2 sequences for the 35 clinical scans.
Statistical Analysis
Noninferiority Analysis.
Statistical analysis was performed using R statistical and computing software, Version 3.4.1 (http://www.r-project.org/).18 A noninferiority analysis between TSE-VFA and conventional TSE was performed using a 1-sided Wilcoxon signed rank test19 for the healthy volunteers and clinical patients. In this analysis, we sought to show the noninferiority of the proposed TSE-VFA compared with the conventional sequence—that is, the proposed technique does not significantly differ from the conventional method by testing the following hypothesis: H0:MVFA − MTSE ≤ −Δ; Ha:MVFA − MTSE > −Δ, where MTSE and MVFA are the median scores of the 2 sequences respectively, H0 is the null hypothesis and Ha is the alternate hypothesis. The clinically acceptable range of differences is defined using the noninferiority margin Δ. In this context, the null hypothesis states that the median difference in image quality scores between TSE and TSE-VFA is greater than a predetermined noninferiority margin Δ. The null hypothesis implies differences in the performance between the 2 sequences, and by rejecting the null at the significance level (α), we show noninferiority in the performance of TSE-VFA compared with TSE. The analysis was performed with a noninferiority margin of Δ = 0.5 and significance level α = .025. The noninferiority margin Δ = 0.5 was chosen on the basis of preliminary data using recommendations from Ahn et al19 for noninferiority studies in radiology. To reduce type I error from multiple comparisons, we adjusted P values from the hypothesis test with the Bonferroni correction factor of 9. A P value < α results in the rejection of the null hypothesis and implies noninferiority of the proposed sequence compared with conventional TSE.
Agreement Analysis
To compare the clinical utility of the 2 sequences, we performed an agreement analysis20,21 using TSE as the non-reference standard for the 5 clinical criteria considered (ie, facet joints, endplates, nerve roots, spinal cord, and discs). A McNemar test of proportions was performed for each of the clinical criteria on the basis of the assumption that a score of ≥4 indicates diagnostic image quality. In accordance with FDA statistical guidelines,21 an overall percentage agreement and the positive percentage agreement were also computed to characterize the diagnostic agreement of the proposed TSE-VFA sequence with TSE.
The Cohen κ statistic was initially used to assess interobserver reliability. However, it was observed that the statistic yielded low values for grading criteria that have a high percentage of agreement and are skewed in the score distribution. To overcome this so-called κ paradox,22 we instead used a weighted Gwet AC1 statistic to measure interobserver variability because it has been shown23 to be more robust to skewed distributions. For HASTE-VFA, images were graded by the 2 observers and the Gwet AC1 was computed.
Results
Figure 1 shows surface plots of the peak and full width at half maximum of the point spread function and SAR computed from the simulated TSE signal evolution as a function of αmin and αcent. For a fixed echo-train length and effective TE, a lower αmin and higher αcent results in better peak signal and resolution, with a very small increase in SAR. Figure 1 also illustrates a VFA scheme and compares the simulated T2 signal decay curves for the constant and proposed VFA echo-trains. From the point spread function plots in Fig 1F, the improvement in spatial resolution with VFA echo-trains can be clearly observed at longer ETLs, albeit with a slight reduction in signal compared with a short ETL constant flip angle sequence.
T2WI generated using the agarose gel phantoms is shown in Fig 2A. To visualize the spatial resolution differences between the acquisitions, a line profile through the 2 phantoms is shown in Fig 2B. The conventional TSE sequence with a constant flip angle exhibits ringing artifacts at long echo-train lengths as observed in the acquired images and the line plots. Table 2 shows the SNR and relative contrast values for the 3 different sequences along with the SAR value computed by the scanner. The proposed TSE-VFA has a 1.5× reduction in SAR compared with the short ETL conventional TSE sequence and 2.3× reduction compared with the long ETL variant. The proposed sequence can generate images with comparable contrast in a 2× shorter scan time. While prolongation of the echo-train and insufficient magnetization recovery lead to a slight loss in SNR, the reduction in scan time results in slightly improved SNR efficiency.
To verify the improvement in point spread function behavior in vivo, we acquired lumbar spine data from 5 healthy volunteers using both TSE and TSE-VFA with an ETL of 56. The conventional TSE sequence with a shorter ETL of 21 was also acquired as a reference. As seen in Fig 3, TSE-VFA has considerably reduced blurring compared with TSE with ETL = 56 and comparable image quality with TSE with ETL = 21. The TSE-VFA also had a 1.6× lower SAR, resulting in scan time reduction from 3.9 to 1.7 minutes. With a long ETL, conventional TSE has poor image quality and higher SAR, resulting in longer scan times (3 min 12 sec) compared with TSE-VFA.
Qualitative assessment of images from 5 healthy volunteers by neuroradiologists (Table 3) showed that TSE-VFA had comparable mean scores for all the different criteria except SNR. Table 4 shows mean values from the ROI analysis on the volunteer images. The proposed sequence generates T2WI with contrast comparable with that of conventional TSE at a lower SAR. However, TSE-VFA had slightly lower measured SNR, consistent with the qualitative scores and the phantom experiments.
Figure 4 shows representative examples of lumbar spine T2WI acquired from 2 subjects using TSE and TSE-VFA sequences. Figure 4A illustrates the presence of motion artifacts on a patient diagnosed with multifocal osseous metastases during the 3 min 24 sec TSE scan. On the other hand, TSE-VFA (1 min 30 sec scan time) generated images with comparable image quality and reduced motion artifacts. Figure 4C, -D shows images from a subject with multilevel degenerative disc disease. The shorter TSE-VFA sequence (D) has resolution comparable with that of the TSE but received a SNR score of 4, while TSE received a score of 5. Note a slight reduction in CSF signal due to insufficient recovery of the magnetization in TSE-VFA.
Table 5 summarizes the mean scores from the 2 observers for TSE and TSE-VFA for lumbar examinations. The median difference in scores between the 2 sequences for lumbar scans was significantly (P < .025) less than the noninferiority margin (Δ = 0.5) for all criteria except SNR, implying noninferiority of TSE-VFA compared with conventional TSE. Note that our null hypothesis states that median difference in the image-quality scores between TSE and TSE-VFA is greater than Δ. The interobserver reliability AC1 scores between readers listed in Table 5 indicate a strong agreement between the 2 readers.
The clinical utility of the TSE-VFA spine protocols is further affirmed by the high positive percentage agreement (>80%) and overall agreement values shown in Table 6 for the lumbar spine. In addition, the McNemar test of proportions showed no significant difference (P > .05) between the conventional and the proposed sequences for the clinical criteria.
The SNR and relative contrast for the lumbar spine were measured for 35 clinical subjects as shown in Table 4. There was a difference in the mean SNR and relative contrast between clinical subjects and healthy volunteers for both TSE and TSE-VFA, presumably due to the hyperintense vertebral body signal caused by edema. There was no significant difference (P = .32) in the vertebrae-disc relative contrast between TSE-VFA and TSE. The proposed sequence had a significantly lower (P < .001) SAR compared with TSE. Note that the mean SAR for TSE was higher in clinical subjects than in healthy volunteers due to higher variability in patient weight. Furthermore, the TSE sequence had a significantly (P = .006) higher SNR than TSE-VFA.
Figure 5 illustrates the utility of the fast HASTE-VFA sequence in a whole-spine protocol acquired in <2 minutes. While HASTE-VFA has a slightly lower resolution (1.25 × 1.0 mm compared with TSE with 1.18 × 0.82 mm), it is about 4 times faster than the TSE, making it viable in a clinical setting and potentially more robust to patient motion. HASTE-VFA (Fig 5B) clearly depicts pathology while exhibiting reduced motion artifacts in the cervical spine section compared with the 9-minute STIR-TSE scan (Fig 5A).
Mean scores for whole-spine scans along with the Gwet AC1 for interobserver agreement are shown in Table 7. There is excellent interobserver agreement in the scores except for SNR, which shows a moderate concordance. The image scores follow a trend similar to that of the lumbar TSE-VFA scores.
Discussion
We have demonstrated the use of variable refocusing flip angle TSE as a fast adjunct for routine spine imaging. Flip angle modulation reduced the blurring and SAR typically associated with long-ETL TSE sequences. The TSE-VFA sequence was roughly 2.3× faster than the conventional TSE sequence (1 min 28 sec versus 3 min 25 sec). The refocusing flip angles of the TSE-VFA sequence were optimized to maximize the spatial resolution and minimize SAR. Image scores were ≥4 for all the image quality and clinical criteria used and for all the clinical criteria considered; TSE-VFA was noninferior to TSE, attesting to the equivalence of the 2 sequences. While we have shown the clinical viability of TSE-VFA for sagittal T2WI, this technique can also be used as a fast alternative to the conventional TSE to generate T1WI and STIR contrasts in axial and sagittal planes. A spine protocol encompassing TSE-VFA for T1, T2, and STIR contrasts would reduce overall scan times.
We have also demonstrated the use of variable refocusing flip angle single-shot sequences as a valuable addition for whole-spine imaging protocols (4× shorter scan time), whereas TSE sequences are prohibitively long. The HASTE-VFA sequence was added to the whole-spine protocol as a fast T2WI scan (1 min 54 sec), and its diagnostic value was confirmed by the excellent image-quality scores (clinical criteria had an average score of 4.71, and image quality metrics showed artifact behavior and SNR comparable with those in the lumbar spine). The HASTE-VFA protocol had a slightly lower resolution than the TSE-VFA (1.25 × 1 mm versus 0.81 × 0.73 mm) to minimize image blurring and maintain SNR, and this limits its use for high-resolution lumbar spine imaging. We are exploring the use of optimal radiofrequency pulses to further reduce echo-spacing and enable higher resolution single-shot T2 imaging.
While the mean motion score of TSE-VFA was noninferior to TSE for the clinical subjects (Table 5), presumably due to the relative lack of patient motion in lumbar spine imaging, we observed a trend toward significance when we eliminated 9 subjects with motion scores of 5 on both sequences (4.43 for TSE-VFA versus 4.04 for TSE, P < .1). This suggests that TSE-VFA could be a good alternative to conventional TSE in patients who cannot tolerate long scans or for cervical spine imaging, in which motion artifacts are common from swallowing.
One possible limitation of the proposed VFA technique is the slight reduction in CSF signal due to insufficient magnetization recovery. This needs to be studied further in a larger clinical study, even though our study showed no detrimental effects from this phenomenon and subtle changes in CSF signal in the spine were adequately visualized in the VFA sequences due to the inherent contrast between CSF and other structures such as the spinal cord and intervertebral discs.
In TSE-VFA sequences, spatial resolution is limited by the available SNR, resulting in slightly lower SNRs for the same resolution as conventional TSE. This feature is indicated both by the slightly lower qualitative SNR scores compared with the TSE sequence (4.36 versus 4.84) and quantitative SNR measurements (45.18 versus 55.45). However, SNR efficiency is slightly higher for TSE-VFA (37.27 versus 29.93) due to the reduced scan time. Preliminary studies on volunteers indicate that the use of anterior coil arrays, atypical in spine imaging, in addition to the posterior coil array would improve SNR and enable comparable spatial resolution between TSE-VFA and TSE. We are currently exploring the feasibility of this clinically.
In this study, the TSE-VFA protocol parameters obtained from numeric simulations have been validated using phantoms and healthy volunteers and then tested on clinical patients. In the patient data, we did see our VFA protocol perform comparably with the conventional TSE, indicating that we are not far off from the volunteer and phantom optimal values. However, an analysis on a larger cohort of clinical patients is necessary to study the optimality of the proposed protocol and its clinical diagnostic utility.
VFA sequences have also been shown17 to be sensitive to motion, a feature determined by the choice of the flip angle parameters. While low values of αmin can enable longer echo-trains, they can cause signal dephasing in the presence of substantial motion. VFA-based sequences are also sensitive to B1 field homogeneities24 because B1 variations cause changes in the refocusing flip angle that could lead to SNR loss in very large patients. We used αmin values of 45°–50°, which were a good compromise between blurring and motion sensitivity.
Conclusions
We have developed fast T2WI spine protocols using variable refocusing flip angles, including a single-shot variant. VFA-based sequences have better point spread function behavior than their constant flip angle counterparts and are less sensitive to patient motion, often observed in longer TSE scans. These features enable the use of VFA sequences as an adjunct or replacement in clinical spine protocols, especially for pediatric 3T spine imaging where TSE scans are highly SAR limited.
Acknowledgments
The authors thank Dr Srinivasan Vedantham for helping with the statistical analysis and interpretation of the results.
Footnotes
Disclosures: Manojkumar Saranathan—UNRELATED: Consultancy: ABC Medical Education.
Paper previously presented, in part, at: Annual Meeting and Exhibition of the International Society for Magnetic Resonance in Medicine, April 22–28, 2017; Honolulu, Hawaii.
References
- Received September 22, 2017.
- Accepted after revision May 14, 2018.
- © 2018 by American Journal of Neuroradiology