An Update in Qualitative Imaging of Bone Using Ultrashort Echo Time Magnetic Resonance
An Update in Qualitative Imaging of Bone Using Ultrashort Echo Time Magnetic Resonance
An Update in Qualitative Imaging of Bone Using Ultrashort Echo Time Magnetic Resonance
Surgery, University of California, San Diego, San Diego, CA, United States, 3 Research Service, Veterans Affairs San Diego
Healthcare System, San Diego, CA, United States
Water in cortical and trabecular bone exist in different states signal when using conventional clinical pulse sequences with
and at various locations (7, 9). In healthy bone, the main portion echo times (TEs) of a several milliseconds or longer (24, 25). The
of water exists in “bound” form to HA crystals and to the lack of direct signal originating from bone impairs the ability of
collagenous matrix (6, 8, 10). The remaining water volume in conventional MRI sequences to provide any qualitative or
bone resides in pores ranging from sub-microns to hundreds of quantitative bone assessments. It should be noted that MRI has
micrometers in size (7, 8). Bound water indicates the bone mineral been used in the past to measure bone microstructure via
and collagenous matrix, while pore water indicates the porosity of indirect visualization of the dark regions (bone) in high-
bone (11, 12). Average bound water T2* is approximately 300 µs resolution conventional acquisitions, however, this approach is
whereas pore water T2* is longer than 1,000 µs and can reach up to limited to distal bone sites and is very motion-sensitive (9, 26).
several milliseconds (10, 13–15). Collagen protons have extremely Recently, new MRI techniques such as ultrashort echo time
short T2*s, on the order of several microseconds (10). (UTE)-MRI have been developed for direct bone imaging and
Bone mineral density (BMD) measurement has been the associated quantitative measurements (5, 6, 16, 22).
standard bone evaluation method in clinics performed using Qualitative bone imaging can be achieved using conventional,
x-ray-based techniques including dual-energy x-ray absorptiometry ultrashort echo time (UTE), adiabatic inversion recovery UTE (IR-
(DEXA) and quantitative computed tomography (QCT) (2, 16, UTE), dual-inversion recovery UTE (Dual-IR-UTE), double-
17). DEXA-based measurement of BMD is non-reliable due to inversion recovery UTE (Double-IR-UTE), UTE with rescaled
very low resolution and its 2D nature. BMD as a predictive echo subtraction (UTE-RS), Fat suppression UTE, Water- and fat-
clinical measurement is quite limited in its representation of suppressed proton projection imaging (WASPI), and zero echo
bone microstructure and, consequently, of bone fragility, time (ZTE) sequences. The contrast mechanisms as well as the
functionality, and fracture risk (18–21). However, these non- advantages and disadvantages of each technique are discussed in
mineral components may describe the bone microstructural and detail. A brief comparison between discussed MRI techniques is
biomechanical properties independently from BMD. Although, presented in Table 1. This review will be an update to our
QCT enables the measurement of bone microstructure in previously published review paper in 2013 (6). As UTE-MRI
addition to BMD however it comes with a high radiation dose. bone imaging field is experiencing fast growth, it is believed that
Employing magnetic resonance imaging (MRI) for bone revisiting this review topic would be benedictional to bone
evaluation has been increasingly reported in the literature. imaging society.
MRI-based techniques for bone evaluation avoid the potential
harm associated with x-ray-based imaging techniques (5, 6, 16,
22). MRI-based bone evaluation can also provide valuable UTE AND IR-UTE PULSE SEQUENCES
evaluation of the surrounding soft tissues including tendons AND THEIR CONTRAST MECHANISM
(23) and muscles, advantages that are not available in x-ray-
based techniques. Bone has a short apparent transverse Both UTE and IR-UTE sequences have been developed for
relaxation time (T2*) and is typically visualized with a void imaging of cortical and trabecular bone as described in the
Conventional FSE (27) Very Low Water in large High (reverse Partial cortical Insensitive Insensitive High High
pores contrast) bone
Conventional STE (28) Very Low Water in large High (reverse Partial cortical Sensitive Insensitive High High
pores contrast) bone
Basic UTE (6, 8, 10, 13, 15, 29) High Bound and Low Cortical bone Insensitive Insensitive Relatively high High
pore water
IR-UTE (25, 30–39) High Bound water High Cortical and Insensitive Insensitive Low High
trabecular bone
Dual-IR UTE (33, 40, 41) High Bound water High Cortical and Sensitive Insensitive Low Moderate
trabecular bone
Double-IR UTE (42) High Bound water High Cortical and Insensitive Insensitive Low High
trabecular bone
Fat suppression UTE (43–45) High Bound and Moderate Cortical and Sensitive Insensitive Moderate Moderate
pore water trabecular bone
UTE-RS (46, 47) High Bound and High Cortical bone Sensitive Insensitive Moderate Moderate
pore water
WASPI (48–51) High Bound water High Cortical and Sensitive Sensitive Relatively low Moderate
trabecular bone
ZTE (52–56) Moderate (low Bound and Low Cortical bone Insensitive Insensitive Relatively high High
flip angle) pore water
following sections. Representative 3D UTE and IR-UTE sequences UTE (IR-UTE) sequences have been developed for selective
are shown in Figure 1. The basic 3D UTE sequence (Figure 1A) imaging of collagen-bound water (Figure 1C). With the IR-UTE
employs a short radiofrequency (RF) rectangular pulse (duration = sequence, a Silver-Hoult adiabatic inversion pulse (duration = 8.64
26–52 µs) for signal excitation followed by 3D radial ramp ms) is used to invert the pore water longitudinal magnetization.
sampling with minimal nominal TEs of 8 to 50 µs depending on The longitudinal magnetization of bound water with a very short
the hardware. For 3D UTE Cones sequence, sampling is T2* cannot be inverted, but can be largely saturated, by the
performed over a k-space divided into multiple cone shapes adiabatic IR pulse due to major relaxation of bound water
with twisted radial trajectories along each cone. The Cones magnetization during the adiabatic inversion pulse. The UTE
trajectories are more time-efficient than radial trajectories in acquisition is initiated to selectively detect signal from bound
covering 3D k-space (46), and resolve the limitations associated water only after an inversion time (TI) when the inverted pore
with 2D UTE sequences, namely, sensitivity to eddy currents (57). water magnetization approaches the null point (Figure 1D).
Furthermore, the 3D UTE Cones sequence allows anisotropic
fields of view and spatial resolution, resulting in reduced scan
times (58–60). For 2D UTE imaging the rectangular pulse is CONVENTIONAL MR FOR BONE IMAGING
replaced with a half RF pulse for slice selective excitation. Both
bound water and pore water can be detected using the basic UTE Conventional MRI sequences generally visualize the bone tissue
sequences (Figure 1B) (14). Adiabatic inversion recovery prepared with a signal void surrounded by bright signal from adjacent soft
A
B
C D
E F
FIGURE 1 | (A) The 3D UTE and (C) IR-UTE sequences, as well as the contrast mechanisms for imaging of (B) total water and (D) collagen-bound water. The UTE
sequence employs a short rectangular pulse for signal excitation followed by 3D radial ramp sampling with a minimal nominal TE of 8 µs. The schematic time point
for UTE acquisition at TE of 8 µs is shown in B. The IR-UTE sequence employs an adiabatic inversion pulse to invert and null the pore water magnetization. The
collagen-bound water magnetization is not inverted and is detected by the subsequent UTE data acquisition. For 2D UTE and IR-UTE imaging the rectangular pulse
is replaced with a half RF pulse for slice selective excitation. (E) Dual-IR and (F) Double-IR UTE sequences combined with the contrast mechanism for bone imaging
(bound water). The Dual-IR UTE sequence utilizes two adiabatic IR pulses with center frequencies of 0 and −440 Hz with narrow bandwidth of around 500 Hz (fat
frequency offset at 3T) to invert long T2 water and fat, respectively. The Double-IR UTE sequence utilizes two identical adiabatic IR pulses with a center frequency of
−220 Hz with broad bandwidth of no less than 1000 Hz (half of the fat frequency offset at 3T) to invert both long T2 water and fat. Bone signal is largely saturated by
each adiabatic IR pulse due to its fast T2* relaxation during the inversion period. At the optimized TI1 and TI2 in both Dual- and Double-IR UTE sequences, both long
T2 water and fat can be suppressed simultaneously, and the bone signals can be acquired with multiple UTE spokes This figure was partially presented before by
Ma et al. (42). Reprinting permission is granted through Rightslink system. The figure is modified for presentation purposes. Minor modifications were performed for
presentation purposes. RF, radio frequency pulse; TE, echo time; Gx,y,z, gradient magnetic field in X, Y, Z directions; DAW, data acquisition window; Mz, longitudinal
magnetization; Mxy, transverse magnetization; CW, collagen-bound water; PW, pore water; IR, inversion recovery; TI, inversion time.
tissues. Pore water may compose up to a quarter of the bone volume Collagen protons have extremely short T2*s (T2* on the order of
and possesses a short T2* and relatively long T2 (up to 100 ms) (10, several microseconds) and are not detectable, even when using
13–15). Therefore, some conventional fast spin echo (FSE) (27) and UTE MRI sequences. The effective TEs in UTE are significantly
short echo time (STE) c sequences have the potential to image pore longer than the T2*s of collagen backbone protons (≈ 10 µs)
water in cortical bone, even though clinical gradient recalled echo because of the use of a relatively long RF excitation pulse as well as
(GRE) sequences are not capable of pore water imaging, as shown time-consuming ramp sampling (8, 10, 15, 29). UTE MRI imaging
in Figure 2 (Figures 2A, B) and Figure 3A. Consequently, FSE and results in high in vivo signal of bone, but it is still lower than the
STE techniques might be useful in qualitative imaging of highly signal from surrounding soft tissues, which in turn leads to
porous cortical bone sites, which can be found more in elderly relatively low contrast as demonstrated in Figure 3B. Improving
cohorts. Nevertheless, qualitative imaging of cortical bone with low the UTE contrast in bone imaging requires fat suppression in
porosity requires more advanced MR imaging techniques. trabecular bone sites while fat and soft tissue suppression in cortical
bone sites. Moreover, pore water suppression can improve the
cortical bone UTE contrast by imaging only the bound water.
UTE MR BONE IMAGING
UTE MRI sequences with nominal TEs around tens of IR-UTE MR IMAGING OF CORTICAL BONE
microseconds or less can detect signals from both bound water
(T2* ≈ 300 µs) and pore water (T2* > 1,000 µs) pools in cortical Adiabatic inversion recovery (IR) preparation pulses have been
bone, as demonstrated in Figure 2 (Figures 2C, D) (6, 13). suggested in different studies to suppress long T2 tissue components
FIGURE 2 | (A) Conventional fast spin echo (FSE), (B) gradient echo (GRE), (C) 3D ultrashort echo time (3D UTE), (D) 3D adiabatic inversion recovery UTE (3D IR-
UTE) MRI images in axial plane of a piece of human cortical bone harvested from the anterior tibial midshaft. FSE and GRE only detect signal from pore water, likely
in Haversian canals (indicated with yellow arrow). UTE MRI results in high signal from all sites of the bone specimen. IR-UTE detects signal from bone with much
higher contrast than UTE MRI can. his figure was previously presented by Du et al. (13). Reprinting permission is granted through Rightslink system. The figure is
modified for presentation purposes. Minor modifications were performed for presentation purposes.
FIGURE 3 | (A) Conventional gradient echo (GRE), (B) ultrashort echo time (UTE), and (C) adiabatic inversion recovery UTE (IR-UTE) image in axial plane of the
lower leg. GRE results in void signal in bone. UTE MRI results in high signal for bone, but low contrast. IR-UTE results in higher contrast between cortical bone and
surrounding soft tissue, similar to that can be seen in CT. This figure was previously presented by Jerban et al. (61). Reprinting permission is granted through
Rightslink system. The figure is modified for presentation purposes. Minor modifications were performed for presentation purposes. 3D UTE sequence parameters:
FOV = 14×14×14 cm3; voxel size = 0.7×0.7×5 mm3, scan time » 3.2 mins. 3D IR-UTE sequence parameters: FOV = 14×14×14 cm3; voxel size = 0.7×0.7×5 mm3,
TR/TI = 300/110 ms, scan time » 3.5 mins.
and increase the contrast in UTE MR imaging of cortical bone The IR-UTE-based qualitative bone imaging sequence has
(Figure 1). The bound water component in cortical bone can be been applied to different bone sites in vivo such as tibia and fibula
selectively imaged with both 2D and 3D IR-UTE sequences. In both (36–38), radius (36), hip (34), and shoulder (39); however, the
cases, a relatively long adiabatic IR pulse (e.g., Silver-Hoult pulse, contrast and image quality depend on the anatomical location,
8.64 ms in duration) is employed to invert the longitudinal coil quality and size, B0 and B1 homogeneity, bone thickness,
magnetizations of first, long T2 water (e.g., free water image resolution, and slice thickness. These factors affect the
components of cortical bone, and water in muscle) and second, proper visualization of the bone structure, signal to noise ratio
fat in cortical bone and bone marrow (25, 30, 32–34). However, the (SNR), suppression uniformity, and image artifacts. Bone sites in
simultaneous suppression of pore water, fat, and muscle is highly axial skeleton or deeply located in the body such as spine and hip
challenging due to the significant differences between MR are more difficult to be scanned compared with tibia and radius
properties of these three proton pools. Thus, IR-UTE based in peripheral bone sites. Specifically, sophisticated thin bone
imaging of bound water may coexist with some level of pore structure requires high resolution imaging which is highly
water, fat, and muscle contamination. The 2D or 3D UTE data sensitive to motion. Larger required coils result in greater B1
acquisition starts at an inversion time (TI) designed to allow the inhomogeneity and non-uniform suppression. Figure 4
inverted free water and fat longitudinal magnetizations to be at or illustrates coronal images of the hip and femoral head of a
close to the null point (Figure 1) (30). Figure 2 shows a comparison healthy young volunteer using 2D FSE and 3D IR-UTE Cones
between conventional MRI (FSE and GRE), UTE, and IR-UTE sequences (31).
sequences performed on a piece of anterior tibial cortex imaged in
axial plane. The conventional FSE sequence produces higher signal
in cortical bone compared with the GRE sequence, though both
techniques only detect signal from pore water residing in Haversian DUAL-IR-UTE IMAGING OF CORTICAL
canals. The UTE sequence can detect free water in the pores (high BONE
signal with fine structure) and bound water (uniform background
signal). The signal of the fine structure will disappear using the IR- Dual-adiabatic inversion recovery (Dual-IR) pulses can also be
UTE sequence which suppresses the free water signal while the employed to invert and null signals from long T2 water and fat,
uniform background signal from bound water remains. respectively (33). Dual-IR followed by selective UTE imaging can
Figure 3 shows tibial and fibular midshaft bone in a healthy detect bound water in cortical bone and provide qualitative bone
young volunteer imaged with GRE, UTE and IR-UTE sequences. evaluations. In this approach, two successive long adiabatic
In UTE MRI, bone shows as high signal compared with inversion pulses are employed to invert the longitudinal
conventional MRI technique, but as low signal compared with magnetization of long T2 water and long T2 fat, respectively (33,
surrounding tissues such as marrow fat and muscle. The IR-UTE 40, 41). The pulse sequence diagram for Dual-IR-UTE is similar to
technique greatly suppresses signals from both the marrow fat IR-UTE sequence hovers with an additional inversion recovery
and muscle, resulting in visualization of cortical bone with a CT- pulse (Figure 1). As a result of a significant transverse relaxation of
like contrast that can be used for qualitative bone evaluation. cortical bone magnetization with short T2 during the long adiabatic
Adiabatic inversion recovery (AIR) UTE (AIR-UTE) (35, 36) inversion process, its longitudinal magnetization cannot be inverted
is an alternative abbreviation used for abovementioned IR-UTE (62). The UTE MRI data acquisition begins after the first delay time
technique by other groups. AIR-UTE is used to visualize bound (TI1), which is required for the nulling of the inverted long T2 water
water in cortical bone and provide a qualitative image of cortical magnetization, and the second delay time (TI2), which is required
bone structure. for the nulling of the inverted fat magnetization. Figure 5 shows
FIGURE 4 | In vivo imaging of the hip of a 24‐year‐old female volunteer with a (A) clinical 2D T2‐weighted FSE and (B) 3D IR‐UTE Cones sequences. Soft tissues
are well‐suppressed in the 3D IR‐UTE Cones image, while they are bright in the clinical T2‐FSE images. This figure was previously presented by Ma et al. (31).
Reprinting permission is granted through Rightslink system. The figure is modified for presentation purposes. Minor modifications were performed for presentation
purposes. 3D IR-UTE sequence parameters: FOV = 38 cm3 × 38 cm3 × 20 cm3; voxel size = 2.4 mm3 × 2.4 mm3 × 5 mm3, TR/TI = 150/64 ms, and scan time ≈
9.5 min.
representative images of the left distal tibia of a volunteer using IR-UTE utilizes two identical adiabatic inversion pulses. The center
clinical 2D GRE, UTE and Dual-IR-UTE techniques. Cortical bone frequency of these identical pulses are located at the water peak
is visualized with a void signal using the 2D GRE pulse sequence, while their spectral widths are broad enough to cover both water
however with a poor contrast using the basic UTE pulse sequence and fat frequencies (42). These two adiabatic inversion pulses are
because of the significant signals acquired from surrounding muscle applied in sequence with two different inversion times (TI1 and
and fat. The Dual-IR-UTE sequence suppresses long T2 water TI2) in order to invert and null the longitudinal magnetizations of
signals, including those from muscle and free water in bone, as long T2 muscle and fat. Similar to Dual-IR-UTE technique, due to
well as from marrow fat, and displays bound water in cortical bone the significant transverse relaxation of the cortical bone
with high signal and contrast. Dual-IR UTE can be applied to magnetization during the long adiabatic inversion pluses, its
different body regions. However, this technique only works well for longitudinal magnetization cannot be inverted. The pulse
simultaneous water and fat suppression when the frequency offset sequence diagram for Double-IR-UTE is similar to IR-UTE
due to the B0 inhomogeneity is less than half of the bandwidth of sequence hovers with an additional inversion recovery pulse
the used adiabatic IR pulse. (Figure 1). Figure 6 shows representative images of the left distal
tibia of a volunteer using clinical 3D GRE, UTE and Double-IR-
UTE techniques (42). The Double-IR-UTE demonstrated cortical
bone with high signal and contrast as the signal from long T2 water
DOUBLE-INVERSION RECOVERY UTE and fat were suppressed robustly. Notably, bone was visualized with
(DOUBLE-IR-UTE) a signal void using the 3D GRE sequence, and with a poor contrast
using the conventional UTE sequence.
Double inversion recovery (Double-IR) UTE sequence also All the IR-UTE-based techniques (IR-UTE, Dual-IR-UTE,
employs two adiabatic inversion pulses to invert and null signals Double-IR-UTE, and AIR-UTE) appear to provide a uniform
from long T2 tissues (42). Despite Dual-IR-UTE sequence, Double- suppression of long T2 water and fat signals, which in turn
FIGURE 5 | The tibial midshaft of a healthy young volunteer imaged with (A) GRE, (B) UTE and dual adiabatic inversion recovery UTE (Dual-IR-UTE) sequences. The
Dual-IR-UTE image (C) selectively suppresses signal from fat and muscle, and which creates high contrast for cortical bone. This figure was previously presented by
Du et al. (33). Reprinting permission is granted through Rightslink system. The figure is modified for presentation purposes. Minor modifications were performed for
presentation purposes. 2D UTE sequence parameters: FOV = 10 cm2 × 10 cm2; voxel size = 0.2 cm2 × 0.2 cm2 and scan time ≈ 3 min. 3D Dual-IR-UTE sequence
parameters: FOV = 10 cm2 × 10 cm2; voxel size = 0.2×0.2, TR/TI1/TI2 = 300/140/110 ms, and scan time ≈ 3 min.
FIGURE 6 | In vivo imaging of the lower leg of young healthy volunteer using (A) clinical 3D GRE sequence (no signal in cortical bone), (B) the conventional 3D UTE
sequence (detects signal in cortical bone signal but with low contrast), and (C) the 3D Double-IR-UTE sequence which shows simultaneous suppression of muscle
and marrow fat, providing higher contrast between cortical bone and surrounding soft tissue. This figure was previously presented by Ma et al. (42). Reprinting
permission is granted through Rightslink system. The figure is modified for presentation purposes. Minor modifications were performed for presentation purposes. 3D
UTE sequence parameters: FOV = 14 cm3 × 14 cm3 × 12 cm3; voxel size = 0.55 mm3 × 0.55 mm3 × 6 mm3, and scan time ≈ 1.3 min. 3D Double-IR-UTE
sequence parameters: FOV = 14×14×12 cm3; voxel size = 1.1 mm3 × 1.1 mm3 × 6 mm3, TR/TI1/TI2 = 200/100/45 ms, and scan time ≈ 2.9 min.
provides a good qualitative image of the targeted bone. This is chemical shifts and average signal oscillation observed in the
because adiabatic IR pulses are relatively insensitive to B1 and B0 multi-echo MRI in T2 fitting analyses (63). Fat suppression
inhomogeneities (41, 62). techniques can be used to remove fat signal contamination in
To the best of our knowledge, there is no detailed comparison bone assessment and to improve the bone contrast in UTE-MRI.
study for these techniques. Based on our experience, IR-UTE Chemical shift fat saturation (FatSat), soft-hard water excitation,
sequence is more time-efficient than both Dual- and Double-IR and single point Dixon methods have been employed to suppress
UTE sequences. IR- and Double-IR UTE sequences are much less fat in UTE bone imaging (43, 44). FatSat is widely used in clinical
sensitive to the B0 inhomogeneity than Dual-IR UTE sequence. MR sequences, however, it is not suitable for bone imaging due to
Dual-IR sequence provides a better SNR compared with Double- the strong signal saturation of the wide spectrum band of bone. The
IR UTE sequence. It would be interesting to perform a future study soft-hard pulse has been proposed to overcome the signal
to compare all these techniques regarding the performance of attenuation effect via utilizing a low power soft-pulse for fat
SNR, CNR, and the corresponding scan efficiency. The excitation in the opposite direction of the following hard pulse
comparisons presented in Table 1 have been prepared based on (43). Figure 7 shows a comparison between tibial and fibular bone
the practical experience of the authors. image contrast in vivo obtained using basic UTE, FatSat, and the
soft-hard water excitation techniques (43). The fat signal could be
suppressed well using both the soft-hard pulse and the FatSat
FAT SUPPRESSION UTE module in the three illustrated slices. However, the cortical bone
signal (indicated by yellow arrow) were much better preserved in
Human bone exists in combination with bone marrow which the soft-hard excitation images (Figures 7D–F) (43). Single-point
possess high percentage of fat. The fat presence results in Dixon method is a postprocessing method to separate water and fat
FIGURE 7 | In vivo UTE Cones imaging of tibial midshaft of young healthy volunteer using (A–C) a single hard pulse excitation (basic UTE), (D–F) the soft-hard water
excitation pulse, and (G-I) the conventional FatSat module. UTE images are presented in three representative slices. Fat was well suppressed by both the soft-hard
pulse and the FatSat module. The cortical bone and coil elements (indicated by yellow arrows in D–F) were much better preserved in the soft-hard excitation images
(D–F) compared with FatSat images (G–I). This figure was previously presented by Ma et al. (43). The reprinting permission is granted through Rightslink system.
The figure is modified for presentation purposes. Minor modifications were performed for presentation purposes. 3D UTE sequence parameters: FOV = 12 cm3 × 12
cm3 × 16 cm3; voxel size = 0.63 mm3 × 0.63 mm3 × 2 mm3, and scan time ≈ 3.4 min.
signals from a dual-echo UTE acquisition (44). The calculated water WATER- AND FAT-SUPPRESSED
and fat maps can then be used to suppress fat in the UTE image. PROTON PROJECTION IMAGING (WASPI)
OF CORTICAL BONE
UTE WITH RESCALED ECHO Water- and fat-suppressed proton projection MRI (WASPI) is
SUBTRACTION (UTE-RS) another MRI sequence developed for selective imaging of bone
water bound to the organic matrix (48–50). In this technique,
UTE with rescaled echo subtraction (UTE-RS) can provide two long-duration, yet low-power, rectangular RF pulses are used
qualitative imaging of the cortical bone (46, 47). In UTE-RS, the to selectively saturate signals from long T2 water and fat. Since
free induction decay (FID) image is scaled down so that signals from bound water has a short T2*, it will remain largely unsaturated
muscle and fat become lower than those from the second echo (46). and provide qualitative imaging of bound water in bone. Figure
In the subtraction image, signals from muscle and fat are negative, 9 shows bone WASPI images in the wrist joint of a healthy
whereas those from short-T2 species are positive, separating them volunteer in transverse, coronal and sagittal planes (51).
from air which has a signal intensity fluctuating around zero (46).
The UTE-RS technique can be efficient in creating high positive
contrast for short T2 species such as cortical bone. Regular unscaled ZTE BONE IMAGING
echo subtraction may reduce bone contrast in this situation.
Performing UTE-RS in trabecular bone sites would be challenging An alternative approach for bone imaging is the zero echo time
due to the fat presence in bone marrow and signal oscillation. (ZTE) sequence, which employs a short rectangular excitation
Rescaled echo subtraction method has been also used with ZTE pulse during the fully ramped up readout gradient, followed by
technique (64). Figure 8 shows the tibial cortex of a healthy young fast radial sampling (52, 53). Alternatively, frequency-modulated
volunteer imaged in the coronal plane using the UTE-RS technique. pulse with interleaved transmit-receive operation, can also be
FIGURE 8 | Dual-echo 3D UTE imaging of the tibial cortex of a healthy young volunteer in coronal plane with (A) TE = 8 ms and (B) TE = 2.2 ms. (C) Subtraction of
the second echo from the first echo image rescaled down by a factor of 0.4. This figure was previously presented by Du et al. (46). The reprinting permission is
granted through Rightslink system. The figure is modified for presentation purposes. Minor modifications were performed for presentation purposes. 3D UTE
sequence parameters: FOV = 24 × 24 cm3 × 24 cm3; voxel size = 0.8 mm3 × 0.8 mm3× 0.8 mm3, TEs = 0.008/2.2 ms, and scan time ≈ 5.5 min.
FIGURE 9 | In vivo water- and fat-suppressed proton projection imaging (WASPI) of the wrist joint of a healthy volunteer in (A) transverse, (B) coronal and (C)
sagittal slices from the 3D WASPI image dataset. This figure was previously presented by Wu et al. (51). Reprinting permission is granted through Rightslink system.
The figure is modified for presentation purposes. Minor modifications were performed for presentation purposes. 3D WASPI sequence parameters: FOV = 12 cm3 ×
12 cm3 × 12 cm3; voxel size = 0.9 mm3 × 0.9 mm3 × 0.9 mm3, and scan time ≈ 18 min.
used for bone imaging which is known as sweep imaging with filled using a Cartesian single-point imaging technique, which is
Fourier transformation (SWIFT) (54). Eliminating the rapid known as pointwise encoding time reduction with radial
gradient switching between TR intervals results in the decreased acquisition (PETRA) (55). Figure 10 shows representative
acoustic noises during scans and in the reduced eddy current postprocessed ZTE and conventional FSE images of the
artifacts (5). These sequences have the potential to image both shoulder of a young symptomatic patient. CT-like ZTE-based
bound water and pore water in bone. These sequences suffer from image (Figure 11A) was obtained after bias-correction and
a gap of data at the center of k-space as a result of a dead time inverse-logarithmic rescaling (52). The qualitative ZTE bone
caused by the finite RF pulse, transmit-receive switching, and image clearly visualizes the bone fragmentation, which is not
digital bandpass filtering (53). The missing data in the ZTE clear in FSE image (52). To improve the image contrast ZTE
technique can be compensated via oversampled acquisition and techniques has been also used with inversion recovery preparation
mathematical reconstruction (53). The ZTE data gap can be also which suppress long-T2 signal (56).
FIGURE 10 | (A) Zero echo time (ZTE)-based and (B) proton-density-weighted FSE MR imaging of the shoulder of a young symptomatic patient. The CT-like ZTE-
based image (A) is obtained after bias-correction and inverse-logarithmic rescaling (52). Qualitative bone imaging by ZTE shows the bone fragment (indicated with
arrow), which is not clear in FSE image. This figure was previously presented by Breighner et al. (52). Reprinting permission is granted from Radiology journal (RSNA).
The figure is modified for presentation purposes. Minor modifications were performed for presentation purposes. 3D ZTE sequence parameters: FOV = 28 cm3 ×
28 cm3 × 28 cm3; voxel size = 0.87 mm3 × 0.87 mm3 × 1.5 mm3, and scan time ≈ 5 min.
FIGURE 11 | 3D IR-UTE-Cones imaging of (A) the femur/spine and (B) the shoulder sample in the coronal plane, demonstrating high contrast imaging of cortical
and trabecular bone at various sites in the body using a clinical whole-body 3T scanner. 3D IR-UTE sequence parameters for femur/spine: FOV = 45×45×20.8 cm3;
voxel size = 2.5×2.5×4 mm3, TR/TI = 183/78 ms, and scan time » 10 mins. 3D IR-UTE sequence parameters for shoulder: FOV = 20× 20 × 10 cm3; voxel size =
0.8×0.8×2 mm3, TR/TI = 140/61 ms, and scan time » 6 mins. This figure is original and based on data from (31). A version of this figure has been presented before
in ISMRM 2019 conference (poster 3754).
REFERENCES 4. Marks Jr. SC, Odgren PR. Chapter 1 - Structure and Development of the
Skeleton. In: LGR John, P Bilezikian, A Lawrence, G Raisz, A Gideon, A
1. NIH Consensus Development Panel on Osteoporosis Prevention, Diagnosis, Rodan, JP Bilezikian, GA Rodan, editors. Principles of Bone Biology, 2nd ed.
and Therapy. Osteoporosis Prevention, Diagnosis, and Therapy. JAMA. San Diego: Academic Press (2002). p. 3–15. doi: 10.1016/B978-012098652-
(2001) 285(6):785–95. doi: 10.1001/jama.285.6.785 1.50103-7
2. MacDonald HM, Nishiyama KK, Kang J, Hanley DA, Boyd SK. Age-related 5. Chang EY, Du J, Chung CB. UTE imaging in the musculoskeletal system. J Magn
patterns of trabecular and cortical bone loss differ between sexes and skeletal Reson Imaging (2015) 41:870–83. doi: 10.1002/jmri.24713
sites: A population-based HR-pQCT study. J Bone Miner Res (2011) 26:50–62. 6. Du J, Bydder GM. Qualitative and quantitative ultrashort-TE MRI of cortical
doi: 10.1002/jbmr.171 bone. NMR BioMed (2013) 26:489–506. doi: 10.1002/nbm.2906
3. Ritchie RO, Buehler MJ, Hansma P. Plasticity and toughness in bone. Phys 7. Cowin SC. Bone poroelasticity. J Biomech (1999) 32:217–38. doi: 10.1016/
Today (2009) 62:41–7. doi: 10.1063/1.3156332 S0021-9290(98)00161-4
8. Wang X, Ni Q. Determination of cortical bone porosity and pore size 28. Akbari A, Abbasi-Rad S, Rad HS. T1 correlates age: A short-TE MR
distribution using a low field pulsed NMR approach. J Orthop Res (2003) relaxometry study in vivo on human cortical bone free water at 1.5T. Bone
21:312–9. doi: 10.1016/S0736-0266(02)00157-2 (2016) 83:17–22. doi: 10.1016/j.bone.2015.10.006
9. Wehrli FW, Song HK, Saha PK, Wright AC. Quantitative MRI for the 29. Diaz E, Chung CB, Bae WC, Statum S, Znamirowski R, Bydder GM, et al.
assessment of bone structure and function. NMR BioMed (2006) 19:731–64. Ultrashort echo time spectroscopic imaging (UTESI): an efficient method for
doi: 10.1002/nbm quantifying bound and free water. NMR BioMed (2012) 25:161–8.
10. Nyman JS, Ni Q, Nicolella DP, Wang X. Measurements of mobile and bound doi: 10.1002/nbm.1728
water by nuclear magnetic resonance correlate with mechanical properties of 30. Du J, Bydder M, Takahashi AM, Chung CB. Two-dimensional ultrashort echo
bone. Bone (2008) 42:193–9. doi: 10.1016/j.bone.2007.09.049 time imaging using a spiral trajectory. Magn Reson Imaging (2008) 26:304–12.
11. Horch RA, Gochberg DF, Nyman JS, Does MD. Non-invasive predictors of doi: 10.1016/j.mri.2007.08.005
human cortical bone mechanical properties: T2-Discriminated 1H NMR 31. Ma YJ, Chen Y, Li L, Cai Z, Wei Z, Jerban S, et al. Trabecular bone imaging
compared with high resolution X-ray. PloS One (2011) 6:1–5. doi: 10.1371/ using a 3D adiabatic inversion recovery prepared ultrashort TE Cones
journal.pone.0016359 sequence at 3T. Magn Reson Med (2020) 83:1640–51. doi: 10.1002/mrm.28027
12. Bae WC, Chen PC, Chung CB, Masuda K, D’Lima D, Du J. Quantitative 32. Du J, Takahashi AM, Chung CB. Ultrashort TE spectroscopic imaging
ultrashort echo time (UTE) MRI of human cortical bone: Correlation with (UTESI): Application to the imaging of short T2 relaxation tissues in the
porosity and biomechanical properties. J Bone Miner Res (2012) 27:848–57. musculoskeletal system. J Magn Reson Imaging (2009) 29:412–21.
doi: 10.1002/jbmr.1535 doi: 10.1002/jmri.21465
13. Du J, Hermida JC, Diaz E, Corbeil J, Znamirowski R, D’Lima DD, et al. 33. Du J, Takahashi AM, Bae WC, Chung CB, Bydder GM. Dual inversion
Assessment of cortical bone with clinical and ultrashort echo time sequences. recovery, ultrashort echo time (DIR UTE) imaging: Creating high contrast for
Magn Reson Med (2013) 70:697–704. doi: 10.1002/mrm.24497 short-T2species. Magn Reson Med (2010) 63:447–55. doi: 10.1002/mrm.22257
14. Biswas R, Bae WC, Diaz E, Masuda K, Chung CB, Bydder GM, et al. 34. Nazaran A, Carl M, Ma Y, Jerban S, Lu X, Du J, et al. Three-dimensional
Ultrashort echo time (UTE) imaging with bi-component analysis: Bound adiabatic inversion recovery prepared ultrashort echo time cones (3D IR-
and free water evaluation of bovine cortical bone subject to sequential drying. UTE-Cones) imaging of cortical bone in the hip. Magn Reson Imaging (2017)
Bone (2012) 50:749–55. doi: 10.1016/j.bone.2011.11.029 44:60–4. doi: 10.1016/j.mri.2017.07.012
15. Horch RA, Nyman JS, Gochberg DF, Dortch RD, Does MD. Characterization 35. Nyman JS, Gorochow LE, Adam Horch R, Uppuganti S, Zein-Sabatto A,
of 1H NMR signal in human cortical bone for magnetic resonance imaging. Manhard MK. Does MD. Partial removal of pore and loosely bound water by
Magn Reson Med (2010) 64:680–7. doi: 10.1002/mrm.22459 low-energy drying decreases cortical bone toughness in young and old donors.
16. Manhard MK, Nyman JS, Does MD. Advances in imaging approaches to J Mech Behav BioMed Mater (2013) 22:136–45. doi: 10.1016/
fracture risk evaluation. Transl Res (2017) 181:1–14. doi: 10.1016/ j.jmbbm.2012.08.013
j.trsl.2016.09.006 36. Manhard MK, Horch RA, Gochberg DF, Nyman JS, Does MD. In Vivo
17. Moser E, Rejnmark L, Sikjaer T, Mosekilde L, Amstrup AK, Jakobsen NFB. Quantitative MR Imaging of Bound and Pore Water in cortical bone.
Association between bone indices assessed by DXA, HR-pQCT and QCT Radiology (2015) 277:221–30. doi: 10.1148/radiol.2015140336
scans in post-menopausal women. J Bone Miner Metab (2015) 34:638–45. 37. Jerban S, Ma Y, Li L, Jang H, Wan L, Guo T, et al. Volumetric Mapping of
doi: 10.1007/s00774-015-0708-9 Bound and Pore Water as well as Collagen Protons in Cortical Bone Using 3D
18. Cummings SR. Are patients with hip fractures more osteoporotic? Review of Ultrashort Echo Time Cones MR Imaging Techniques. Bone (2019) 127:120–
the evidence. Am J Med. (1985) 78(3):487–94. doi: 10.1016/0002-9343(85) 8. doi: 10.1016/j.bone.2019.05.038
90343-2. 38. Zhao X, Song HK, Seifert AC, Li C, Wehrli FW. Feasibility of assessing bone
19. Marshall D, Johnell O, Wedel H. Meta-analysis of how well measures of bone matrix and mineral properties in vivo by combined solidstate 1H and 31P
mineral density predict occurrence of osteoporotic fractures. Br Med J (1996). MRI. PloS One (2017) 12:1–16. doi: 10.1371/journal.pone.0173995
312(7041):1254–9. doi: 10.1136/bmj.312.7041.1254 39. Ma Yj, West J, Nazaran A, Cheng X, Hoenecke H, Du J, et al. Feasibility of
20. Faulkner KG. Bone matters: are density increases necessary to reduce fracture using an inversion-recovery ultrashort echo time (UTE) sequence for
risk? J Bone Miner Res. (2000) 15(2):183–7. doi: 10.1359/jbmr.2000.15.2.183 quantification of glenoid bone loss. Skeletal Radiol (2018) 47:973–80.
21. Schuit SCE, Van Der Klift M, Weel AEAM, De Laet CEDH, Burger H, Seeman doi: 10.1007/s00256-018-2898-4
E, et al. Fracture incidence and association with bone mineral density in 40. Du J, Carl M, Bae WC, Statum S, Chang EY, Bydder GM, et al. Dual inversion
elderly men and women: The Rotterdam Study. Bone (2004) 34(1):195–202. recovery ultrashort echo time (DIR-UTE) imaging and quantification of the
doi: 10.1016/j.bone.2003.10.001 zone of calcified cartilage (ZCC). Osteoarthr Cartil (2013) 21(1):77–85.
22. Wehrli FW. Magnetic resonance of calcified tissues. J Magn Reson (2013) doi: 10.1016/j.joca.2012.09.009
229:35–48. doi: 10.1016/j.jmr.2012.12.011 41. Garwood M, DelaBarre L. The return of the frequency sweep: Designing
23. Jerban S, Ma Y, Namiranian B, Ashir A, Shirazian H, Zhao W, et al. Age- adiabatic pulses for contemporary NMR. J Magn Reson (2001) 153(2):155–77.
related decrease in collagen proton fraction in tibial tendons estimated by doi: 10.1006/jmre.2001.2340
magnetization transfer modeling of ultrashort echo time magnetic resonance 42. Ma YJ, Zhu Y, Lu X, Carl M, Chang EY, Du J. Short T 2 imaging using a 3D
imaging (UTE-MRI). Sci Rep (2019) 9(1):17974. doi: 10.1038/s41598-019- double adiabatic inversion recovery prepared ultrashort echo time cones (3D
54559-3 DIR-UTE-Cones) sequence. Magn Reson Med (2017) 00:1–9. doi: 10.1002/
24. Robson MD, Gatehouse PD, Bydder M, Bydder GM. Magnetic Resonance: An mrm.26908
Introduction to Ultrashort TE (UTE) Imaging. J Comput Assist Tomogr (2003) 43. Ma Y, Jerban S, Jang H, Chang EY, Du J. Fat suppression for ultrashort echo
27:825–46. doi: 10.1097/00004728-200311000-00001 time imaging using a novel soft-hard composite radiofrequency pulse. Magn
25. Reichert ILH, Robson MD, Gatehouse PD, He T, Chappell KE, Holmes J, et al. Reson Med (2019) 82:2178–87. doi: 10.1002/mrm.27885
Magnetic resonance imaging of cortical bone with ultrashort TE pulse 44. Jang H, Carl M, Ma Y, Jerban S, Guo T, Zhao W, et al. Fat Suppression for
sequences. Magn Reson Imaging (2005) 23:611–8. doi: 10.1016/ Ultrashort Echo Time Imaging Using a Single Point Dixon Method. NMR
j.mri.2005.02.017 BioMed (2019) 32:e4069. doi: 10.1002/nbm.4069
26. Rajapakse CS, Magland JF, Wald MJ, Liu XS, Zhang XH, Guo XE, et al. 45. Wurnig MC, Calcagni M, Kenkel D, Vich M, Weiger M, Andreisek G, et al.
Computational Biomechanics of the Distal Tibia from High- Resolution MR Characterization of trabecular bone density with ultra-short echo-time MRI at
and Micro-CT Images. Bone (2010) 47:556–63. doi: 10.1016/ 1.5, 3.0 and 7.0 T - comparison with micro-computed tomography. NMR
j.bone.2010.05.039 BioMed (2014) 27:1159–66. doi: 10.1002/nbm.3169
27. Bae WC, Patil S, Biswas R, Li S, Chang EY, Statum S, et al. Magnetic resonance 46. Du J, Bydder M, Takahashi AM, Carl M, Chung CB, Bydder GM. Short T2
imaging assessed cortical porosity is highly correlated with mCT porosity. contrast with three-dimensional ultrashort echo time imaging. Magn Reson
Bone (2014) 66:56–61. doi: 10.1016/j.bone.2014.06.004 Imaging (2011) 29:470–82. doi: 10.1016/j.mri.2010.11.003
47. Johnson EM, Vyas U, Ghanouni P, Pauly KB, Pauly JM. Improved cortical 61. Jerban S, Ma Y, Nazaran A, Dorthe EW, Cory E, Carl M, et al. Detecting stress
bone specificity in UTE MR Imaging. Magn Reson Med (2017) 77:684–95. injury (fatigue fracture) in fibular cortical bone using quantitative ultrashort
doi: 10.1002/mrm.26160 echo time-magnetization transfer (UTE-MT): An ex vivo study. NMR Biomed
48. Wu Y, Ackerman JL, Chesler DA, Graham L, Wang Y, Glimcher MJ. Density (2018) 31:e3994. doi: 10.1002/nbm.3994
of organic matrix of native mineralized bone measured by water- and fat- 62. Larson PEZ, Conolly SM, Pauly JM, Nishimura DG. Using adiabatic inversion
suppressed proton projection MRI. Magn Reson Med (2003) 50:59–68. pulses for long-T2 suppression in ultrashort echo time (UTE) imaging. Magn
doi: 10.1002/mrm.10512 Reson Med (2007) 58:952–61. doi: 10.1002/mrm.21341
49. Cao H, Ackerman JL, Hrovat MI, Graham L, Glimcher MJ, Wu Y. 63. Lu X, Jerban S, Wan L, Ma Y, Jang H, Le N, et al. Three Dimensional
Quantitative bone matrix density measurement by water- and fat- Ultrashort Echo Time Imaging with Tri-component Analysis for Human
suppressed proton projection MRI (WASPI) with polymer calibration Cortical Bone. Magn Reson Med (2019) 82:348–55. doi: 10.1002/mrm.27718
phantoms. Magn Reson Med (2008) 60(6):1433–43. doi: 10.1002/ 64. Lee HM, Weiger M, Giehr C, Froidevaux R, Brunner DO, Rösler MB, et al.
mrm.21771 Long-T2-suppressed zero echo time imaging with weighted echo subtraction
50. Cao H, Nazarian A, Ackerman JL, Snyder BD, Rosenberg AE, Nazarian RM, and gradient error correction. Magn Reson Med (2020) 83:412–26.
et al. Quantitative 31P NMR spectroscopy and 1H MRI measurements of bone doi: 10.1002/mrm.27925
mineral and matrix density differentiate metabolic bone diseases in rat 65. Wehrli FW. Structural and functional assessment of trabecular and cortical
models. Bone (2010) 46(6):1582–90. doi: 10.1016/j.bone.2010.02.020 bone by micro magnetic resonance imaging. J Magn Reson Imaging (2007)
51. Wu Y, Hrovat MI, Ackerman JL, Reese TG, Cao H, Ecklund K, et al. Bone 25:390–409. doi: 10.1002/jmri.20807
matrix imaged in vivo by water- and fat-suppressed proton projection MRI 66. Robson MD, Gatehouse PD, Bydder GM, Neubauer S. Human Imaging of
(WASPI) of animal and human subjects. J Magn Reson Imaging (2010) 31 Phosphorus in Cortical and Trabecular Bone In Vivo. Magn Reson Med
(4):954–63. doi: 10.1002/jmri.22130 (2004) 51:888–92. doi: 10.1002/mrm.20055
52. Breighner RE, Endo Y, Konin GP, Gulotta LV, Koff MF, Potter HG. Zero echo 67. Wu Y, Reese TG, Cao H, Hrovat MI, Toddes SP, Lemdiasov RA, et al. Bone
time imaging of the shoulder: Enhanced osseous detail by using MR imaging. mineral imaged in vivo by 31P solid state MRI of human wrists. J Magn Reson
Radiology (2018) 286:960–6. doi: 10.1148/radiol.2017170906 Imaging (2011) 34:623–33. doi: 10.1002/jmri.22637
53. Weiger M, Pruessmann KP. MRI with Zero Echo Time. eMagRes, (2012) 1: 68. Zhao X, Song HK, Wehrli FW. In vivo bone 31P relaxation times and their
311–22. doi: 10.1002/9780470034590.emrstm1292 implications on mineral quantification. Magn Reson Med (2018) 80:2514–24.
54. Garwood M, Idiyatullin D, Corum CA, Chamberlain R, Moeller S, Kobayashi N, doi: 10.1002/mrm.27230
et al. Capturing Signals from Fast-relaxing Spins with Frequency-Swept MRI: 69. Wu Y, Dai G, Ackerman JL, Hrovat MI, Glimcher MJ, Snyder BD, et al.
SWIFT. EMagRes (2012) 1:322–32. doi: 10.1002/9780470034590.emrstm1259 Water- and fat-suppressed proton projection MRI (WASPI) of rat femur
55. Grodzki DM, Jakob PM, Heismann B. Ultrashort echo time imaging using bone. Magn Reson Med (2007) 57:554–67. doi: 10.1002/mrm.21174
pointwise encoding time reduction with radial acquisition (PETRA). Magn 70. Weiger M, Stampanoni M, Pruessmann KP. Direct depiction of bone
Reson Med (2011) 67:510–8. doi: 10.1002/mrm.23017 microstructure using MRI with zero echo time. Bone (2013) 54:44–7.
56. Weiger M, Wu M, Wurnig MC, Kenkel D, Boss A, Andreisek G, et al. ZTE doi: 10.1016/j.bone.2013.01.027
imaging with long-T2 suppression. NMR BioMed (2015) 28:247–54. 71. Carl M, Bydder GM, Du J. UTE imaging with simultaneous water and fat
doi: 10.1002/nbm.3246 signal suppression using a time-efficient multispoke inversion recovery pulse
57. Wansapura JP, Daniel BL, Pauly J, Butts K. Temperature mapping of frozen sequence. Magn Reson Med (2016) 76:577–82. doi: 10.1002/mrm.25823
tissue using eddy current compensated half excitation RF pulses. Magn Reson 72. Tannús A, Garwood M. Adiabatic pulses. NMR Biomed (1997) 10:423–34.
Med (2001) 46:985–92. doi: 10.1002/mrm.1285 doi: 10.1002/(SICI)1099-1492(199712)10:8<423::AID-NBM488>3.0.CO;2-X
58. Hong W, He Q, Fan S, Carl M, Shao H, Chen J, et al. Imaging and
quantification of iron-oxide nanoparticles (IONP) using MP-RAGE and Conflict of Interest: The authors declare that the research was conducted in the
UTE based sequences. Magn Reson Med (2016) 78:226–32. doi: 10.1002/ absence of any commercial or financial relationships that could be construed as a
mrm.26371 potential conflict of interest.
59. Chen J, Carl M, Ma Y, Shao H, Lu X, Chen B, et al. Fast volumetric imaging of
bound and pore water in cortical bone using three-dimensional ultrashort-TE Copyright © 2020 Jerban, Chang, Ma, Jang, Chang and Du. This is an open-access
(UTE) and inversion recovery UTE sequences. NMR BioMed (2016) 29:1373– article distributed under the terms of the Creative Commons Attribution License
80. doi: 10.1002/nbm.3579 (CC BY). The use, distribution or reproduction in other forums is permitted, provided
60. Ma Y, Carl M, Shao H, Tadros AS, Chang EY, Du J. Three-dimensional the original author(s) and the copyright owner(s) are credited and that the original
ultrashort echo time cones T 1r (3D UTE-cones-T 1r ) imaging. NMR publication in this journal is cited, in accordance with accepted academic practice. No
BioMed (2017) 30:e3709. doi: 10.1002/nbm.3709 use, distribution or reproduction is permitted which does not comply with these terms.