Bone microarchitecture has been shown to provide useful information regarding the evaluation of skeleton quality with an added value to areal bone mineral density, which can be used for the di-agnosis of several bone diseases. Bone mineral density estimated from dual-energy x-ray absorp-tiometry (DXA) has shown to be a limited tool to identify patients’ risk stratification and therapy delivery. Magnetic resonance imaging (MRI) has been proposed as another technique to assess bone quality and fracture risk by evaluating the bone structure and microarchitecture.
A large number of studies have demonstrated the substantial burden of bone disorders worldwide [1][2][3]. Considered as the second greatest cause of disability [1], musculoskeletal pathologies account for 6.8% of total disability worldwide [2]. Bone pathologies are usually affecting the bones solid phase, which is composed of both cortical and cancellous/trabecular types of bone. Bone alterations commonly include cortical shell thinning, increased porosity of both cortical and trabecular bone phases [4][5], and reduced density, volume, and regenerative power. These bone modifications generally account for a reduced resistivity and flexibility eventually leading to an increased risk of fragility fractures accompanied by long-term disabilities. Recent studies have shown that people over the age of 50 with a high risk of osteoporotic fractures represented more than 150 million people worldwide with 137 million women [6]. This number is expected to exceed 300 million by 2040 [6]. Fragility fractures lead to more than half a million hospitalizations each year in North America alone, with an annual direct cost, which has been estimated to be $17 billion dollars in 2005. This cost is expected to rise by almost 50% by 2025 [7]. Overall, the early identification of bone fragility risk is a major health issue [8]. In the clinical context, bone disorders are usually assessed using dual-energy X-ray absorptiometry (DXA), which is able to assess the bone mineral density (BMD). The BMD score is then compared to a reference range of values calculated in healthy (25–35 years old) volunteers taking into account sex and ethnicity. Accordingly, a score (T-score) is generated indicating how far, in terms of SD (standard deviation), the measured BMD is from the reference values. A T-score between −1 and −2.5 indicates a low bone mass or osteopenia while a value lower than −2.5 is indicative of osteoporosis. The corresponding method has good sensitivity (around 88% for both men and post-menopausal women), but the specificity is poor (around 41% for post-menopausal women and 55% for men) [9] resulting in a low clinical diagnostic accuracy (70%) [10]. In addition, DXA measurements do not take into consideration microarchitectural alterations, which have also been recognized as part of the structural picture in osteoporosis. Of interest, bone microarchitecture can be assessed using quantitative computed tomography (qCT) [11][12]. Given that both DXA and qCT are both radiative imaging techniques, non-radiative alternatives would be of great interest. Over the last decades, magnetic resonance imaging (MRI) [13][14][15] has been indicated as a non-ionizing and non-invasive technique.
Using MRI, a large number of studies have attempted to assess bone microarchitecture in bone disorders and more particularly in osteoporosis [16][17][18][19]. The corresponding studies have been conducted at different magnetic field strengths, using different Radio Frequency coils and pulse sequences. Although, the results were compelling, the sensitivity of the corresponding microarchitecture metrics for diagnostic purposes and the assessment of the disease severity is still a matter of debate.
On the basis of a comparative survey of MRI, computed tomography, and DXA-based metrics, we intended to address the issues related to the diagnostic potential of the corresponding metrics and their capacity to predict disease severity. The final section will be devoted to potential perspectives offered by magnetic resonance spectroscopy (MRS) and chemical shift encoding (CSE-MRI), solid-state MRI, and quantitative susceptibility mapping (QSM).
Bone is a multiphase material composed of a solid phase and a viscoelastic component. The solid phase is considered as hierarchical, anisotropic, and heterogeneous and is composed of 65% of inorganic matrix (mostly calcium hydroxyapatite crystals) and 35% of organic matrix (type I collagen, proteoglycans, and bound water) [20]. While the inorganic matrix is characterized by a high rigidity, a high resistivity, and an elastic behavior, the organic matrix is deformable thereby providing the tissue with tensile strength. Due to the combination of these two materials, bone tissue is simultaneously deformable and rigid [21]. The solid phase creates a shell for the bone marrow, which is the viscoelastic component. The bone marrow on the other hand has a double function. It provides nutriments to the solid phase allowing higher regenerative rate and is able, due to its viscoelastic properties, to spread the dynamics of an impulsive action, reducing the risk of fractures due to impacts [22]. Bone tissue is composed of both trabecular and cortical bone phases. Cortical bone covers the whole surface of the bone. It is compact, dense, and characterized by overlapped and parallel lamellae, which provide a large resistivity [20]. Trabecular bone is the inner compartment of bone tissue. It is composed of 25% of bone and 75% of marrow [23]. At the microstructural level, trabecular bone appears as a complex 3D network of interconnected trabeculae rods and plates responsible for tissue resistance to loading forces. The bone inner architecture is an important contributor to bone strength independent of bone mass [20]. It is characterized by a high porosity so that trabecular bone is lighter and less dense than cortical bone. In fact, cortical bone mainly works in compression while trabecular bone principally works in flexion and torsion reaching a higher area under the stress–strain curve [23].
Bone is actually a dynamic porous structure and this porosity can change as a result of pathological processes but also as an adaptive response to mechanical or physiological stimuli. This change in both cortical and trabecular bone porosity can strongly affect the corresponding mechanical properties [23].
A non-invasive alternative to DXA and qCT could be MRI. Over the last two decades, a large number of studies have intended to assess bone microstructure using MRI. The initial investigations have been performed using T1-weighted spin echo sequences characterized by short TR (<1200 ms) and short TE (<25 ms) in distal radius and calcaneus [16][24][25][26]. Due to technical advances, tibiae [17][27][28], spine [24][29], and proximal femur [18][30][31][32] have been investigated. MRI of trabecular microstructure can be obtained by imaging the marrow phase inside the bone segment, which appears as a hyperintense signal in conventional MR images. Using higher field MRI, i.e., 3T one can expect an increased signal to noise ratio (SNR), which can be translated either in a reduced acquisition time or an increased image resolution. Over the last decades, due to the higher availability of high-field (HF) MRI scanners, a large number of studies have been dedicated to the MRI assessment of osteoporosis [17][18][30][33][28][32][31]. Very recently, clinical FDA and CE-approved ultra-high field (i.e., 7T UHF) MRI scanners with announced MSK applications have become available. Their clinical availability is still poor and the coming results will be of utmost importance to decide about the future of UHF MRI for clinical purposes.
Using MRI, the most common extrapolated features are the bone volume fraction (BVF), the trabecular thickness (Tb.Th), spacing (Tb.Sp), and number (Tb.N) [18][33].
A signal to noise ratio (SNR) of 10 has been reported as the minimum value for the investigation of bone microarchitecture [34]. The scan time considered acceptable for clinical examination has to range between 10 and 15 min. As a result the minimum voxel size, which has been obtained at 1.5T was between 0.135 and 0.250 mm while the slice thickness was between 0.3 and 1.5 mm. One has to keep in mind that SNR would be higher for superficial anatomical sites (radius or calcaneus compared to deeper anatomical sites, e.g., proximal femur) leading to higher resolution or shorter acquisition time. Moreover, SNR can be increased at higher field strengths and/or using multichannel coils [34][35][36][37].
MRI pulse sequences such as gradient recalled echo (GRE) and spin echo (SE) have also been tested at different field strengths [17][32][38]. It has been shown that SE sequences were less susceptible to partial volume effects as compared to GRE sequences and that GRE were more sensitive to trabecular broadening than SE. These results indicate that SE sequences would provide more accurate results regarding trabecular characteristics [17][38]. However, the use of these pulse sequences might be problematic using ultra-high field (UHF) MRI considering power-deposition issues.
A list of the main literature references, scanned regions, sequences, and principal MRI setup parameters is reported in Table 1.
Table 1. List of the main magnetic resonance imaging (MRI) parameters and sequences.
Anatomical Site |
Clinical History |
Specimen /Patient |
Acq. Time |
Sl. Thickness [mm] [mm] |
Pix. Size [mm] |
FOV [mm] |
Sequence |
Main Field |
N° |
Reference |
distal radii |
type 2 diabetes |
patient |
12 min 9 s |
1 |
0.195 × 0.195 |
100 × 100 |
FSE |
1T |
[39] |
Pritchard et al. |
calcaneus |
osteoporotic hip fractures |
patient |
15 min 15 s |
0.5 |
0.195 × 0.195 |
100 × 100 |
GE |
1.5T |
[26] |
Link et al. |
distal radii |
healthy |
patient |
16 min 25 s |
0.5 |
0.156 × 0.156 |
80 × 45 |
3D FLASE |
1.5T |
[36] |
Techawiboonwong et al. |
distal radii |
healthy |
patient |
3 min 15 s |
0.5 |
0.156 × 0.156 |
80 × 45 |
3D SSFP |
1.5T |
[36] |
Techawiboonwong et al. |
distal radii |
NA |
specimen |
15 min |
0.3 |
0.156 × 0.156 |
80 |
GE |
1.5T |
[13] |
Majumdar et al. |
lumbar spine |
osteoporotic |
patient |
16 min |
0.7 |
0.156 × 0.156 |
80 × 80 |
GE |
1.5T |
[24] |
Majumdar et al. |
distal radii |
hip fractures |
patient |
NA |
0.5 |
0.156 × 0.156 |
80 × 80 |
GE |
1.5T |
[16] |
Majumdar et al. |
distal radii |
NA |
specimen |
58 min (1) 16 min (2) |
0.3 (1) 0.9 (2) |
0.153 × 0.153 |
49×78 |
SE |
1.5T |
[40] |
Link et al. |
prox. femur |
NA |
specimen |
74 min (1) 27 min (2) |
0.3 (1) 0.9 (2) |
0.195 × 0.195 |
75 × 100 |
SE |
1.5T |
[41] |
Link et al. |
prox. femur |
healthy |
patient |
6 min 12 s |
1.5 |
0.234 × 0.234 |
NA |
3D FIESTA |
1.5T |
[32] |
Krug et al. |
distal tibiae |
NA |
specimen |
40 min |
0.16 |
0.160 × 0.160 |
70 × 63 |
3D FLASE |
1.5T |
[42] |
Rajapakse et al. |
lumbar spine |
NA |
specimen |
15 min 23 s |
0.41 |
0.137 × 0.137 |
70 × 64 × 13 |
3D FLASE |
1.5T |
[29] |
Rajapakse et al. |
distal radii(1) distal tibiae(2) |
osteopenic and osteoporotic |
patient |
12 min (1) 16 min (2) |
0.4 |
0.137 × 0.137 |
70 × 40(1) 70 × 50(2) |
3D FLASE |
1.5T |
[25] |
Ladinsky et al. |
distal femur |
cerebral palsy (children) |
patient |
9 min 52 s |
0.7 |
0.175 × 0.175 |
90 |
3D fast GE |
1.5T |
[43] |
Modlesky et al. |
distal radii(1) distal tibi.ae(2) |
osteoporotic |
patient |
12 min (1) 16 min (2) |
0.41 |
0.137 × 0.137 |
70 × 40 × 13 (1) 70 × 50 × 13 (2) |
3D FLASE |
1.5T |
[44] |
Rajapakse et al. |
prox. femur |
NA |
specimen |
16 min 55 s |
1.1 |
0.21 × 0.21 |
120 |
TSE |
3T |
[45] |
Soldati et al. |
prox. femur |
healthy |
patient |
12 min 43 s |
1.5 |
0.234 × 0.235 |
NA |
3D FIESTA |
3T |
[32] |
Krug et al. |
distal radii, distal tibiae |
NA |
specimen |
< 10 min |
0.5 |
0.156 × 0.156 |
NA |
GE |
3T |
[38] |
Krug et al. |
distal radii, distal tibiae |
NA |
specimen |
< 10 min |
0.5 |
0.156 × 0.156 |
NA |
GRE |
3T |
[38] |
Krug et al. |
distal radii, distal tibiae |
NA |
specimen |
< 10 min |
0.5 |
0.156 × 0.156 |
NA |
SE |
3T |
[38] |
Krug et al. |
distal tibiae |
osteoporotic |
patient |
15 min |
0.41 |
0.137 × 0.137 |
70 × 64 × 13 |
3D FLASE |
3T |
[28] |
Zhang et al. |
prox. femur |
fragility fractured |
patient |
25 min 30 s |
1.5 |
0.234 × 0.234 |
120 |
FLASH |
3T |
[30] |
Chang et al. |
prox. femur |
long-term glucocorticoid |
patient |
15 min 18 s |
1.5 |
0.234 × 0.234 |
100 |
FLASH |
3T |
[31] |
Chang et al. |
distal radii |
HR+ breast cancer |
patient |
7 min |
0.34 |
0.170 × 0.170 |
65 |
GE |
3T |
[46] |
Baum et al. |
distal femur |
osteoarthritis |
patient |
9 min 18 s |
1 |
0.180 × 0.180 |
100 |
3D B-FFE |
3T |
[47] |
Liu et al. |
prox. tibia |
osteoarthritis |
patient |
3 min |
2.8 |
0.230 × 0.240 |
120 × 123 |
SE |
3T |
[48] |
MacKey et al. |
prox. tibia, distal femur |
osteoarthritis |
patient |
NA |
1 |
0.195 × 0.195 |
100 |
FIESTA-c |
3T |
[49] |
Chiba et al. |
prox. tibia, distal femur |
osteoarthritis |
patient |
NA |
1 |
0.195 × 0.195 |
160 |
SPGR |
3T |
[49] |
Chiba et al. |
distal tibiae |
NA |
specimen |
7 min |
0.41 |
0.137 × 0.137 |
70 × 53 × 13 |
3D FLASE |
3T |
[19] |
Rajapakse et al. |
prox. femur |
NA |
specimen |
16 min 45 s |
1.5 |
0.13 × 0.13 |
130 |
TSE |
7T |
[50] |
Soldati et al. |
prox. femur |
NA |
specimen |
37 min 36 s |
0.5 |
0.170 × 0.170 |
140 × 140 |
GRE |
7T |
[33] |
Guenoun et al. |
distal tibiae |
healthy |
patient |
19 min 10 s |
0.5 |
0.156 × 0.156 |
NA |
SE |
7T |
[17] |
Krug et al. |
distal tibiae |
healthy |
patient |
18 min 25 s |
0.5 |
0.156 × 0.157 |
NA |
FP |
7T |
[17] |
Krug et al. |
vertebrae (1 axial, 2 sagittal) |
NA |
specimen |
34 min (1) 51 min (2) |
0.4 (1) 0.5 (2) |
0.170 × 0.170 |
140 × 140 |
GRE |
7T |
[51] |
Guenoun et al. |
distal femur |
fragility fractured |
patient |
7 min 9 s |
1 |
0.234 × 0.234 |
120 |
FLASH |
7T |
[18] |
Chang et al. |
femurs, tibiae, vertebrae |
NA |
specimen |
120 min |
0.05 |
0.05 × 0.05 |
6.4 × 6.4 × 25.6 |
SE |
9.4T |
[52] |
Rajapakse et al. |
In the majority of MRI literature, the morphological parameters that are reported are BVF, Tb.Th, Tb.Sp, and Tb.N [32][38][39]. In addition, some groups have proposed some other features such as an erosion index, trabecular rod- and plate-like structures, trabecular plate-to-rod ratio, trabecular isolation, and fractal lacunarity [18][40].
These microarchitectural parameters have been generated from the post-processing of both 2D and 3D images. The corresponding analyses were performed in binarized images or in original grey level intensities. All these approaches have tried to take into account partial volume effects occurring given the poor resolution of MRI as compared to the trabeculae dimension [33][41][42]. So far, no standard reference has been suggested.
Studies performed at different MRI field strength in postmenopausal woman with fragility fractures have illustrated microstructural alterations (reduced BVF and increased Tb.Sp) whereas DXA T-scores were unchanged. In a study conducted in distal radii at 1.5T, Kijowsky et al. showed that post-menopausal woman had a slightly lower (−9%) bone volume fraction and a higher erosion index (+17%) compared to controls [43]. Krug et al. in a study conducted on the proximal femurs of six healthy males and females using both 1.5T and 3T MRI showed good correlation (r up to 0.86) between structural parameters obtained from the two different field strengths. However, they reported that bone structure of the proximal femur was substantially better depicted at 3T than 1.5T [32]. Microstructure alterations have been reported in a large variety of cases including chronic kidney disease (CKD) [11][53], HIV-infection [44], glucocorticoid-induced osteoporosis [31], or disuse osteoporosis [45].
In a 3T MRI study conducted in distal tibiae of 20 patients with CKD, Ruderman et al. reported trabecular deterioration together with reduced cortical thickness [53]. Moreover, a study conducted on 30 patients affected by end stage renal disease (ESRD) it has been shown that Tb.N, Tb.Th, and whole bone stiffness were significantly lower (p < 0.01) is ESRD compared to controls [46]. A similar study conducted on distal tibiae of 11 kidney transplant recipient patients have high-lightened post-transplant deterioration in trabecular bone quality [47]. In a study conducted in proximal femurs at 3T, glucocorticoid treated patients had a largely reduced (−50.3%) Tb.N, trabecular plate-to-rod ratio (−20.1%), and a largely increased (+191%) Tb.Sp [31]. Patients with a disuse osteoporosis displayed similar anomalies for BVF (−30%), Tb.N (−21%), Tb.Th (−12%), and Tb.Sp (+48%) [45]. Chang et al. [18] further supported and extended these results in a study conducted in distal femur at 7T. In 31 subjects with fragility fractures, they reported a lower BVF (–3%), Tb.N (–6%), and erosion index (–6%). Moreover, in a 7T MRI study conducted in the distal radius of 24 women, Griffin et al. reported a trabecular bone microarchitecture gradient with an overall higher quality (+123% BVF, +16% Tb.N) distally (epiphysis) than proximally (diaphysis) [48].
Ultra-high field MRI can provide images with a smaller pixel size (0.156 mm × 0.156 mm) as compared to the resolution achieved at lower field strength (0.234 mm × 0.234 mm for example at 3T [17][32]). In a dual 3T-7T study conducted in distal tibiae of 10 healthy volunteers, Krug et al. reported that metrics computed at higher field strength were different than those quantified from 3T MR images. More specifically, UHF measurements illustrated increased BVF (+22%) and Tb.Th (+25%) whereas Tb.Sp (−21%) and Tb.N (−4%) were both decreased [88]. These results suggest a higher discriminative power of UHF MRI for trabecular features.
In a study conducted in 32 postmenopausal women, Kang et al. showed a good correlation between DXA-based BMD and MRI T2 and T2 * in calcaneus (r = −0.8, p < 0.001) and spine (r = −0.53, p = 0.002) [50]. Similar results have been reported for the femoral neck [65,90] with a good correlation (r = 0.74, p < 0.001) between DXA-based BMD and T2 * values [52]. T2 * relaxation time illustrates the susceptibility differences between trabecular and bone marrow leading to signal loss due to magnetic field inhomogeneities. MRI-derived T2 * has been shown to correlate with DXA results in several anatomical areas such as calcaneus, distal radius, and Ward’s area in the femoral neck [54][55]. Based on T2 * measurements, Schmeel et al. reported a significant difference between benign and malignant neoplastic vertebral compression fractures (VCFs). A 72% diagnostic accuracy was computed [56]. Furthermore, a strong negative correlation was found between the pelvic bone marrow adipose tissue (BMAT) calculated in 56 healthy women using MRI and the corresponding DXA-based BMD (r = −0.646, p < 0.001) [57]. The negative correlation indicates that patients with decreased bone mineral density are characterized by an increased fat content in bone marrow [57][58][59].
Based on highly resolved MR images (0.150–0.300 mm in-plane pixel size), Chang et al. showed a lack of significant correlation between DXA-computed BMD T-scores and MRI computed microarchitectural parameters in the femoral neck in both controls and glucocorticoid-treated patients [31][60]. Similar results were also reported more recently in subchondral tibiae [61], proximal femurs [33][42], vertebrae [62], and on patients affected by diabetes [63][64]. Guenoun et al. reported that the combination of BVF and BMD was able to improve the prediction of the failure stress (from r2 = 0.384 for BMD alone to r2 = 0.414). All the presented results suggest that although density and structure metrics illustrate bone quality, microarchitectural parameters provide additional information regarding skeletal fragility.
Results from the literature showed that image resolution is a key parameter for the assessment of bone microarchitecture. Importantly, a distinction must be made between in-plane and through-planes resolution. For specific oriented plane (mostly perpendicular to the trabecular), an in-plane MRI pixel size in the same order of magnitude than Tb.Th dimension is enough to measure morphological parameters similar to those extrapolated using gold standard method and so both ex vivo (µCT) [42][65] and in vivo (HR-pQCT) [38]. If one intends to assess bone microstructure using small isovolumetric voxels (0.15 mm), close to the actual thickness of the trabeculae, with an acceptable SNR, acquisition times would exceed the in vivo acceptable duration. One can increase the SNR and reduce the acquisition time with an increased slice thickness while keeping the plane pixel size constant. Accordingly, the radius morphological parameters computed from similar in-plane pixel sizes and different slice thicknesses (0.156 mm × 0.156 mm × 0.3 mm [13], 0.156 mm × 0.156 mm × 0.5 mm [16][66], 0.156 mm × 0.156 mm × 0.7 mm [24], and 0.153 mm × 0.153 mm × 0.9 mm [67]) were comparable. In fact, the bone inner microarchitecture appeared to be a mixture of oriented plates- and rod-like structures. The parallel trabecular plates structures are separated by bone marrow and are perpendicular to the coronal plane [68]. On that basis, increasing the in-plane pixel size should provide more accurate results independently of the slice thickness. As reported by Mulder et al., the calculated volume of ellipsoid at high resolution (0.1 mm × 0.1 mm) is independent from the anisotropy factor but related to the orientation [69].
Different studies performed in distal radii at 1.5T, using similar in-plane pixel size and using different slice thicknesses above 0.3 mm, reported comparable morphological results [13][16][43]. However, in a study conducted by Majumdar et al. in 39 distal radii specimens acquired using 1.5T MRI and contact radiograph, 0.9-mm thick MR images performed better than those obtained from 0.3-mm images. This was explained with the significantly higher SNR (18.2 in 0.9-mm thick images and 9.3 in 0.3-mm sections) [67]. Similar results were obtained in vivo in distal radii scanned at 1.5T (0.156 mm × 0.156 mm × 0.5 mm) with an acceptable SNR around 10 [36]. Moreover, wrists and distal tibiae scanned in patients using 1.5T with pixel sizes in the same range of trabecular thickness (0.156 mm × 0.156 mm × 0.410 mm [25] and 0.137 mm × 0.137 mm × 0.410 mm [70]) reporting lower acquisition time for wrist (12 min) than for tibiae (16 min) and good image quality in both anatomical regions. In a second study conducted by Majumdar et al., 31 cadaveric proximal femurs were scanned at 1.5T with an in-plane pixel size of 0.195 × 0.195 and comparing two different slice thicknesses (0.9 and 0.3 mm). The SNR achieved was 25.2 and 13.8 for the larger and smaller slice thickness respectively. The corresponding acquisition times were very long (27:19 and 73:14 min), i.e., much longer than what could be accepted in clinics [71].
The knee articulation has also been assessed in the study of Rajapakse et al., 17 distal tibiae specimens were scanned at 3T (0.137 mm × 0.137 mm × 0.410 mm) in 7 min [19]. These results where extended in vivo by Zhang et al., in the distal tibiae of 20 postmenopausal women with osteoporosis. The scanning time using 3T MRI (0.137 mm × 0.137 mm × 0.410 mm) was less than 15 min [28]. Krug et al. further confirmed these results in a study comparing 3T MRI (0.156 mm × 0.156 mm × 0.5 mm) and X-ray based techniques both ex vivo (5 tibiae and 3 radii) and in vivo (5 radii and 6 tibiae). While the scanning time was less than 10 min, correlations were reported between both methods and so for the whole set of parameters, i.e., BVF (r = 0.83) and Tb.Sp (r = 0.7) [72]. Liu et al. also reported 3T MR images (pixel size 0.180 mm × 0.180 mm, acq. time 9:18 min) of 92 distal femurs divided in three groups (without osteoarthritis, mild osteoarthritis, and severe osteoarthritis) reporting progressively lower BVF and higher erosion index from healthy patients to those affected by severe osteoarthritis [73], extending previous results [74][75][76].
The technical advantages of moving from 1.5T to 3T or 7T MR scanners were clearly visible in the acquisition of deeper anatomical sites keeping the spatial pixel size in the same order of the trabecular thickness, the acquisition time (acq. Time), and the SNR (>10) being clinically compatible. On that basis, 7T MR scanners have been tested mostly for the acquisition of distal and proximal femur, which represent a clinical important fracture site and one of the most invalidating [77].
In a comparative study conducted in vivo in proximal femur at 1.5 and 3T, Krug et al., reported as expected a 1.6 time-SNR increase together with a corresponding contrast-to-noise ratio (CNR) increase at higher magnetic field. While the 3T images clearly showed the trabecular bone structure, the image resolution did not allow a proper trabecular morphological analysis [32]. In a more recent study in the knee joint of 16 healthy volunteers scanned at 1.5T (0.6 mm × 0.6 mm × 0.6 mm, acq. time 7:15 min) and 3T (0.5 mm × 0.5 mm × 0.5 mm, acq. time 6:51 min), Abdulaal et al. reported significantly higher SNR (p < 0.05) allowing a better trabecular characterization at 3T than 1.5T [78]. Moreover, 3T MRI could be used to successfully scan radii with an in-plane pixel size comparable to the trabecular thickness and an acquisition time (10 min) lower than what commonly needed at 1.5T [66][79]. Jarraya et al., on a study conducted in 50 distal radii scanned at both 3T (0.2 mm × 0.2 mm × 2.0 mm, acq. time 4:29 min) and 7T MR (0.125 mm × 0.125 mm × 2.0 mm, acq. time 3:16 min), reported a statistical significant difference of horizontal and fractal dimensions between patients with chronic wrist disease and controls [80]. A similar comparative analysis has been performed between 3T and 7T MRI (0.156 mm × 0.156 mm × 0.5 mm, acq. time lower than 10 min) and HR-pQCT. Krug et al. showed that tibial trabecular structures were over-represented at higher field strength. Due to susceptibility-induced broadening smaller trabeculae normally not visible due to partial volume effects may be emphasized at 7T [49]. Moreover, using UHF MRI (0.234 mm × 0.234 mm × 1.0 mm, acq. time 7 min), Chang et al. reported that microarchitectural parameters could discriminate between patients and controls and could detect bone deterioration in women with fragility fractures for whom BMD was normal [18]. In addition to the effects of magnetic field strength, Krug et al. also assessed the potential differences between GRE and TSE sequences at 7T. SNR was slightly higher for GRE sequences (13.2 vs. 11.9) while the bone marrow signal was more homogeneous using TSE sequences. This large homogeneity is related to a reduced susceptibility-induced broadening of the trabeculae so that the morphological analysis showed decreased BVF (−13%) and Tb.Th (−23%). These values were closer to those reported using the HR-pQCT reference method [49]. Furthermore, in a study conducted in three cadaveric proximal femurs scanned at 7T (0.130 mm × 0.130 mm × 1.5 mm, acq. time 16 min) and using µCT, Soldati et al. reported no statistical difference between the methods and so for the whole set of morphological parameters [42]. These preliminary results strongly suggest that UHF MRI could be of interest for the in vivo assessment of bone microarchitecture particularly for the deep anatomical regions.
Validation of the bone morphological parameters derived from the high-resolution MR images has usually been performed through the comparison with X-ray based techniques (qCT, HRpQCT, and μCT).
Ex vivo studies have been performed in different body parts. However, due to the samples size (<5 cm3) and the commonly used preparation protocols (replacement of marrow), they remain poorly representative of the in vivo conditions [13][39][42][67][81]. One of the first studies validating MR bone structure measurements was performed by Hipp et al. in cubic bovine trabecular bone from several anatomical sites using optical and micro-MRI methods. BVF and Tb.N were linearly related (r2 = 0.81 and r2 = 0.53 respectively) and did not differ statistically (p = 0.96 and p = 0.17) [39][81]. These results were confirmed and extended in human specimens by Majumdar et al., in a study conducted in 7 cubic specimens of trabecular bone extracted from cadaveric radii scanned at 1.5T (0.156 mm × 0.156 mm × 0.3 mm) and using μCT (0.018 mm isovolumetric). The results showed a good correlation for the whole set of metrics with BVF and Tb.Th performing the best (r = 0.77 and 0.87 respectively) and Tb.Sp and Tb.N the worst (r = 0.53 and 0.6 respectively). However a significative statistical difference (p > 0.01) was reported for all the calculated features [13]. MRI images with an in plane pixel-size lower than the smallest trabecular thickness order (0.1 mm) are not easily reachable. On that basis, one cannot expect to fully characterize it. Moreover, these findings were further extended in a larger study conducted in 39 distal radius specimens scanned at 1.5T MRI (0.152 mm × 0.152 mm × 0.9 mm) and using contact radiography (0.05 mm isovolumetric). The results showed a significant correlation (r > 0.61) between bone microstructure parameters derived from both methods with Tb.Sp and BVF providing the highest correlations (r = 0.69 and p = 0.75 respectively) [67]. More recently, Rajakapse et al. conducted a study in 13 cylindrical specimens (7 proximal femurs, 3 proximal tibiae, and 3 third lumbar vertebrae) extracted from 7 human donors and computed microarchitectural parameters using 9.4T micro-MRI (0.050 mm isovolumetric) and μCT (0.021 mm isovolumetric). Architectural parameters were found to highly correlate between these two modalities with a slope close to unity (r2 ranging from 0.78 to 0.97) [82]. In a more recent study conducted in three cadaveric entire proximal femurs evaluating the trabecular morphology using 7T MRI (0.13 mm × 0.13 mm × 1.5 mm) and comparing the results with those acquired using μCT (0.051 mm isovolumetric) (Figure 1), Soldati et al. showed a good intraclass correlation coefficient for all the parameters (ICC > 0.54) between 7T and μCT [42] illustrating that bone morphological metrics of human specimens can be properly assessed using MRI. Moreover, due to the comparison between MR images and gold standard high-resolution CT images, it has been shown that trabecular features derived from images with a similar pixel size provide statistically comparable results. However, when assessing bone trabeculae using MRI, partial volume effects will occur and will affect image segmentation and trabeculae quantification.
Figure 1. Comparison between MRI and CT. (first row) MR images of in vivo distal tibia acquired using gradient echo sequence at 7T MRI (a) (0.156 mm × 0.156 mm × 0.5 mm) and 3T MRI (b) (0.156 mm × 0.156 mm × 0.5 mm), and compared with high-resolution peripheral computed tomography (HR-pQCT) (c) (0.082 mm3) (reproduced from J. of Mag. Res. Im. 27:854–859 (2008)). (second row) MR images of cadaveric proximal femur acquired using turbo spin echo sequence at 7T MRI (d) (0.13 mm × 0.13 mm × 1.5 mm) and 3T MRI (e) (0.21 mm × 0.21 mm × 1.1 mm), and compared with µCT (f) (0.051 mm3). Note that using MRI, the trabecular bone appears black and bone marrow delivers the bright signal whereas for HR-pQCT and µCT the trabecular bone is shown bright. Additionally, note that the trabecular network is clearly more enhanced at 7T compared to 3T.
The MRI potential for the bone microstructure has also been assessed in vivo in anatomical regions more affected by osteoporosis, i.e., tibiae and radii, vertebrae [24][83][84], distal [18][73][74][75][85], and proximal femurs [30][32][31]. Microarchitectural parameters extrapolated from 3T MRI (0.156 mm × 0.156 mm × 0.5 mm) and compared to HR-pQCT of tibiae and radii of 11 healthy volunteers showed good correlation for BVF (r = 0.83) and Tb.Sp (r = 0.7) in tibiae and good correlation for all the microarchitecture parameters investigated in radii (r = 0.65, 0.95, 0.83, and 0.63 for BVF, Tb.N, Tb.Sp, and Tb.Th respectively) [38]. Kazakia et al. extended these results in a study conducted in tibiae and radii of 52 postmenopausal scanned at 3T MRI (0.156 mm × 0.156 mm × 0.5 mm) and using HR-pQCT. A significant correlation between MRI and HR-pQCT has been reported for Tb.N (r2 = 0.52) and Tb.Sp (r2 = 0.54–0.60) with no statistical difference for these two parameters. Poor correlations were reported for BVF and Tb.Th (r2 = 0.18–0.34) [86]. Similar results were also reported by Folkesson et al., in a study conducted in 52 postmenopausal women scanned at 3T (0.156 mm × 0.156 mm × 0.5 mm) and using HR-pQCT in both tibiae and radii. All the structural parameters derived from MRI were highly correlated to those obtained from HR-pQCT (Tb.N was equal to 0.68 and 0.73 and Tb.Sp was equal to 0.77 and 0.67 for tibiae and radii respectively) with the exception of BVF and Tb.Th for which correlations were less significant (BVF was equal to 0.61 and 0.39 and Tb.Th was equal to 0.43 and 0.32 for tibiae and radii respectively) [79]. Furthermore, Krug et al. confirmed and extended these results in a study conducted in distal tibiae of 10 healthy volunteers scanned at 3T and 7T (0.156 mm × 0.156 mm × 0.5 mm for both techniques). The results showed that microarchitectural parameters extracted from HR-pQCT images had higher correlation with those extracted from 7T MR images (r equal to 0.73 for BVF, 0.69 for Tb.N, 0.89 for Tb.Sp, and 0.13 for Tb.Th) as compared to 3T MR images (r = 0.83, 0.49, 0.67, and 0.15 for BVF, Tb.N, Tb.Sp, and Tb.N respectively) (Figure 1). Interestingly, the corresponding absolute values did only differ by 0.6% for 7T and 3% for 3T [49]. All the findings reported above indicate good correlations for Tb.Sp and Tb.N between MRI and HR-pQCT. In contrast, this was not the case for BVF and Tb.Th. The limited resolution in MRI leads to partial volume effects responsible for the exclusion of the smallest trabeculae, while susceptibility artifacts enhance the remaining trabeculae leading to an overestimation of Tb.Th. This double effect seems limited when using UHF MRI. Indeed, good correlations were found between MRI and HR-pQCT metrics although a poor correlation was still existing for Tb.Th.
This entry is adapted from the peer-reviewed paper 10.3390/ijms22052509