Sodium MRI of human articular cartilage of the wrist: a feasibility study on a clinical 3T MRI scanner

Sodium MRI of the Lumbar Intervertebral Discs of the Human Spine: An Ex Vivo Study

BACKGROUND

Lower back pain affects 75%-85% of people at some point in their lives. The detection of biochemical changes with sodium (23Na) MRI has potential to enable an earlier and more accurate diagnosis.

PURPOSE

To measure 23Na relaxation times and apparent tissue sodium concentration (aTSC) in ex-vivo intervertebral discs (IVDs), and to investigate the relationship between aTSC and histological Thompson grade.

STUDY TYPE

Ex-vivo.

SPECIMEN

Thirty IVDs from the lumbar spines of 11 human body donors (4 female, 7 male, mean age 86 ± 8 years).

FIELD STRENGTH/SEQUENCE

3 T; density-adapted 3D radial sequence (DA-3D-RAD).

ASSESSMENT

IVD 23Na longitudinal (T1), short and long transverse (T2s* and T2l*) relaxation times and the proportion of the short transverse relaxation (ps) were calculated for one IVD per spine sample (11 IVDs). Furthermore, aTSCs were calculated for all IVDs. The degradation of the IVDs was assessed via histological Thompson grading.

STATISTICAL TESTS

A Kendall Tau correlation (τ) test was performed between the aTSCs and the Thompson grades. The significance level was set to P < 0.05.

RESULTS

Mean 23Na relaxation parameters of a subset of 11 IVDs were T1 = 9.8 ± 1.3 msec, T2s* = 0.7 ± 0.1 msec, T2l* = 7.3 ± 1.1 msec, and ps = 32.7 ± 4.0%. A total of 30 IVDs were examined, of which 3 had Thompson grade 1, 4 had grade 2, 5 had grade 3, 5 had grade 4, and 13 had grade 5. The aTSC decreased with increasing degradation, being 274.6 ± 18.9 mM for Thompson grade 1 and 190.5 ± 29.5 mM for Thompson grade 5. The correlation between whole IVD aTSC and Thompson grade was significant and strongly negative (τ = -0.56).

DATA CONCLUSION

This study showed a significant correlation between aTSC and degenerative IVD changes. Consequently, aTSC has potential to be useful as an indicator of degenerative spinal changes.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

Recent technical developments and clinical research applications of sodium (23Na) MRI

Sodium is an essential ion that plays a central role in many physiological processes including the transmembrane electrochemical gradient and the maintenance of the body's homeostasis. Due to the crucial role of sodium in the human body, the sodium nucleus is a promising candidate for non-invasively assessing (patho-)physiological changes. Almost 10 years ago, Madelin et al. provided a comprehensive review of methods and applications of sodium (23Na) MRI (Madelin et al., 2014) [1]. More recent review articles have focused mainly on specific applications of 23Na MRI. For example, several articles covered 23Na MRI applications for diseases such as osteoarthritis (Zbyn et al., 2016, Zaric et al., 2020) [2,3], multiple sclerosis (Petracca et al., 2016, Huhn et al., 2019) [4,5] and brain tumors (Schepkin, 2016) [6], or for imaging certain organs such as the kidneys (Zollner et al., 2016) [7], the brain (Shah et al., 2016, Thulborn et al., 2018) [8,9], and the heart (Bottomley, 2016) [10]. Other articles have reviewed technical developments such as radiofrequency (RF) coils for 23Na MRI (Wiggins et al., 2016, Bangerter et al., 2016) [11,12], pulse sequences (Konstandin et al., 2014) [13], image reconstruction methods (Chen et al., 2021) [14], and interleaved/simultaneous imaging techniques (Lopez Kolkovsky et al., 2022) [15]. In addition, 23Na MRI topics have been covered in review articles with broader topics such as multinuclear MRI or ultra-high-field MRI (Niesporek et al., 2019, Hu et al., 2019, Ladd et al., 2018) [16-18]. During the past decade, various research groups have continued working on technical improvements to sodium MRI and have investigated its potential to serve as a diagnostic and prognostic tool. Clinical research applications of 23Na MRI have covered a broad spectrum of diseases, mainly focusing on the brain, cartilage, and skeletal muscle (see Fig. 1). In this article, we aim to provide a comprehensive summary of methodological and hardware developments, as well as a review of various clinical research applications of sodium (23Na) MRI in the last decade (i.e., published from the beginning of 2013 to the end of 2022).

Sodium Dynamics in the Cellular Environment

Sodium ions are essential for the functions of biological cells, and they are maintained at the balance between intra- and extracellular environments. The quantitative assessment of intra- and extracellular sodium as well as its dynamics can provide crucial physiological information on a living system. ²³Na nuclear magnetic resonance (NMR) is a powerful and noninvasive technique to probe the local environment and dynamics of sodium ions. However, due to the complex relaxation behavior of the quadrupolar nucleus in the intermediate-motion regime and because of the heterogeneous compartments and diverse molecular interactions in the cellular environment, the understanding of the ²³Na NMR signal in biological systems is still at the early stage. In this work, we characterize the relaxation and diffusion of sodium ions in the solutions of proteins and polysaccharides, as well as in the in vitro samples of living cells. The multi-exponential behavior of ²³Na transverse relaxation has been analyzed according to the relaxation theory to derive the crucial information related to the ionic dynamics and molecular binding in the solutions. The bi-compartment model of transverse relaxation and diffusion measurements can corroborate each other to quantify the fractions of intra- and extracellular sodium. We show that ²³Na relaxation and diffusion can be used to monitor the viability of human cells, which offers versatile NMR metrics for in vivo studies.

Evaluation of Sodium Relaxation Times and Concentrations in the Achilles Tendon Using MRI

Sodium magnetic resonance imaging (MRI) can be used to evaluate the change in the proteoglycan content in Achilles tendons (ATs) of patients with different AT pathologies by measuring the ²³Na signal-to-noise ratio (SNR). As ²³Na SNR alone is difficult to compare between different studies, because of the high influence of hardware configurations and sequence settings on the SNR, we further set out to measure the apparent tissue sodium content (aTSC) in the AT as a better comparable parameter. Ten healthy controls and one patient with tendinopathy in the AT were examined using a clinical 3 Tesla (T) MRI scanner in conjunction with a dual tuned ¹H/²³Na surface coil to measure ²³Na SNR and aTSC in their ATs. ²³Na T₁ and T₂* of the AT were also measured for three controls to correct for different relaxation behavior. The results were as follows: ²³Na SNR = 11.7 ± 2.2, aTSC = 82.2 ± 13.9 mM, ²³Na T₁ = 20.4 ± 2.4 ms, ²³Na T_2s* = 1.4 ± 0.4 ms, and ²³Na T_2l* = 13.9 ± 0.8 ms for the whole AT of healthy controls with significant regional differences. These are the first reported aTSCs and ²³Na relaxation times for the AT using sodium MRI and may serve for future comparability in different studies regarding examinations of diseased ATs with sodium MRI.

Chemical Exchange Saturation Transfer for Lactate-Weighted Imaging at 3 T MRI: Comprehensive In Silico, In Vitro, In Situ, and In Vivo Evaluations

Based on in silico, in vitro, in situ, and in vivo evaluations, this study aims to establish and optimize the chemical exchange saturation transfer (CEST) imaging of lactate (Lactate-CEST—LATEST). To this end, we optimized LATEST sequences using Bloch−McConnell simulations for optimal detection of lactate with a clinical 3 T MRI scanner. The optimized sequences were used to image variable lactate concentrations in vitro (using phantom measurements), in situ (using nine human cadaveric lower leg specimens), and in vivo (using four healthy volunteers after exertional exercise) that were then statistically analyzed using the non-parametric Friedman test and Kendall Tau-b rank correlation. Within the simulated Bloch−McConnell equations framework, the magnetization transfer ratio asymmetry (MTRasym) value was quantified as 0.4% in the lactate-specific range of 0.5−1 ppm, both in vitro and in situ, and served as the imaging surrogate of the lactate level. In situ, significant differences (p < 0.001) and strong correlations (τ = 0.67) were observed between the MTRasym values and standardized intra-muscular lactate concentrations. In vivo, a temporary increase in the MTRasym values was detected after exertional exercise. In this bench-to-bedside comprehensive feasibility study, different lactate concentrations were detected using an optimized LATEST imaging protocol in vitro, in situ, and in vivo at 3 T, which prospectively paves the way towards non-invasive quantification and monitoring of lactate levels across a broad spectrum of diseases.

Fast, Accurate, and Robust T2 Mapping of Articular Cartilage by Neural Networks

For T2 mapping, the underlying mono-exponential signal decay is traditionally quantified by non-linear Least-Squares Estimation (LSE) curve fitting, which is prone to outliers and computationally expensive. This study aimed to validate a fully connected neural network (NN) to estimate T2 relaxation times and to assess its performance versus LSE fitting methods. To this end, the NN was trained and tested in silico on a synthetic dataset of 75 million signal decays. Its quantification error was comparatively evaluated against three LSE methods, i.e., traditional methods without any modification, with an offset, and one with noise correction. Following in-situ acquisition of T2 maps in seven human cadaveric knee joint specimens at high and low signal-to-noise ratios, the NN and LSE methods were used to estimate the T2 relaxation times of the manually segmented patellofemoral cartilage. In-silico modeling at low signal-to-noise ratio indicated significantly lower quantification error for the NN (by medians of 6−33%) than for the LSE methods (p < 0.001). These results were confirmed by the in-situ measurements (medians of 10−35%). T2 quantification by the NN took only 4 s, which was faster than the LSE methods (28−43 s). In conclusion, NNs provide fast, accurate, and robust quantification of T2 relaxation times.

Quantification of Sodium Relaxation Times and Concentrations as Surrogates of Proteoglycan Content of Patellar CARTILAGE at 3T MRI

Sodium MRI has the potential to depict cartilage health accurately, but synovial fluid can influence the estimation of sodium parameters of cartilage. Therefore, this study aimed to reduce the impact of synovial fluid to render the quantitative compositional analyses of cartilage tissue technically more robust. Two dedicated protocols were applied for determining sodium T1 and T2* relaxation times. For each protocol, data were acquired from 10 healthy volunteers and one patient with patellar cartilage damage. Data recorded with multiple repetition times for T1 measurement and multi-echo data acquired with an additional inversion recovery pulse for T2* measurement were analysed using biexponential models to differentiate longitudinal relaxation components of cartilage (T1,car) and synovial fluid (T1,syn), and short (T2s*) from long (T2l*) transversal relaxation components. Sodium relaxation times and concentration estimates in patellar cartilage were successfully determined: T1,car = 14.5 ± 0.7 ms; T1,syn = 37.9 ± 2.9 ms; c(T1-protocol) = 200 ± 48 mmol/L; T2s* = 0.4 ± 0.1 ms; T2l* = 12.6 ± 0.7 ms; c(T2*-protocol) = 215 ± 44 mmol/L for healthy volunteers. In conclusion, a robust determination of sodium relaxation times is possible at a clinical field strength of 3T to quantify sodium concentrations, which might be a valuable tool to determine cartilage health.

In-Situ Cartilage Functionality Assessment Based on Advanced MRI Techniques and Precise Compartmental Knee Joint Loading through Varus and Valgus Stress

Stress MRI brings together mechanical loading and MRI in the functional assessment of cartilage and meniscus, yet lacks basic scientific validation. This study assessed the response-to-loading patterns of cartilage and meniscus incurred by standardized compartmental varus and valgus loading of the human knee joint. Eight human cadaveric knee joints underwent imaging by morphologic (i.e., proton density-weighted fat-saturated and 3D water-selective) and quantitative (i.e., T1ρ and T2 mapping) sequences, both unloaded and loaded to 73.5 N, 147.1 N, and 220.6 N of compartmental pressurization. After manual segmentation of cartilage and meniscus, morphometric measures and T2 and T1ρ relaxation times were quantified. CT-based analysis of joint alignment and histologic and biomechanical tissue measures served as references. Under loading, we observed significant decreases in cartilage thickness (p < 0.001 (repeated measures ANOVA)) and T1ρ relaxation times (p = 0.001; medial meniscus, lateral tibia; (Friedman test)), significant increases in T2 relaxation times (p ≤ 0.004; medial femur, lateral tibia; (Friedman test)), and adaptive joint motion. In conclusion, varus and valgus stress MRI induces meaningful changes in cartilage and meniscus secondary to compartmental loading that may be assessed by cartilage morphometric measures as well as T2 and T1ρ mapping as imaging surrogates of tissue functionality.

Longitudinal T2 Mapping and Texture Feature Analysis in the Detection and Monitoring of Experimental Post-Traumatic Cartilage Degeneration

Background: Traumatic cartilage injuries predispose articulating joints to focal cartilage defects and, eventually, posttraumatic osteoarthritis. Current clinical-standard imaging modalities such as morphologic MRI fail to reliably detect cartilage trauma and to monitor associated posttraumatic degenerative changes with oftentimes severe prognostic implications. Quantitative MRI techniques such as T2 mapping are promising in detecting and monitoring such changes yet lack sufficient validation in controlled basic research contexts. Material and Methods: 35 macroscopically intact cartilage samples obtained from total joint replacements were exposed to standardized injurious impaction with low (0.49 J, n = 14) or high (0.98 J, n = 14) energy levels and imaged before and immediately, 24 h, and 72 h after impaction by T2 mapping. Contrast, homogeneity, energy, and variance were quantified as features of texture on each T2 map. Unimpacted controls (n = 7) and histologic assessment served as reference. Results: As a function of impaction energy and time, absolute T2 values, contrast, and variance were significantly increased, while homogeneity and energy were significantly decreased. Conclusion: T2 mapping and texture feature analysis are sensitive diagnostic means to detect and monitor traumatic impaction injuries of cartilage and associated posttraumatic degenerative changes and may be used to assess cartilage after trauma to identify “cartilage at risk”.

Differentiating rheumatoid and psoriatic arthritis: a systematic analysis of high-resolution magnetic resonance imaging features-preliminary findings

BACKGROUND

Because of overlapping phenotypical presentations, the diagnostic differentiation of rheumatoid arthritis (RA) and psoriatic arthritis (PsA) remains challenging. Thus, this study aimed to examine the diagnostic value of distinct imaging features obtained by high-resolution 3-T MRI for the diagnostic differentiation.

MATERIALS AND METHODS

Seventeen patients with PsA and 28 patients with RA were imaged at high resolution using 3-T MRI scanners and a dedicated 16-channel hand coil. All images were analyzed according to the outcome measures in rheumatology clinical trials' (OMERACT) RAMRIS (Rheumatoid Arthritis Magnetic Resonance Imaging Score) and PsAMRIS (Psoriatic Arthritis Magnetic Resonance Imaging Score) for the presence and intensity of synovitis, flexor tenosynovitis, bone edema, bone erosion, periarticular inflammation, bone proliferation, and joint space narrowing. Next, odds ratios (OR) were calculated to determine the strength of the associations between these imaging features, demographic characteristics, and the outcome RA vs. PsA.

RESULTS

PsA could be differentiated from RA by extracapsular inflammatory changes (PsAMRIS sub-score "periarticular inflammation"), with low odds for the presence of RA (OR of 0.06, p < 0.01) at all metacarpophalangeal (MCP) joints. A prediction model informed by the items that were strongest associated with the presence of RA or PsA demonstrated excellent differentiating capability with an area under the curve of 98.1%.

CONCLUSION

High-resolution imaging is beneficial for the identification of relevant imaging features that may assist the clinical differentiation of inflammatory conditions of the hand. At the MCP level, extracapsular inflammatory changes were strongly associated with PsA and may consequently allow the imaging differentiation of PsA and RA.

Aggrecan, the Primary Weight-Bearing Cartilage Proteoglycan, Has Context-Dependent, Cell-Directive Properties in Embryonic Development and Neurogenesis: Aggrecan Glycan Side Chain Modifications Convey Interactive Biodiversity

This review examines aggrecan's roles in developmental embryonic tissues, in tissues undergoing morphogenetic transition and in mature weight-bearing tissues. Aggrecan is a remarkably versatile and capable proteoglycan (PG) with diverse tissue context-dependent functional attributes beyond its established role as a weight-bearing PG. The aggrecan core protein provides a template which can be variably decorated with a number of glycosaminoglycan (GAG) side chains including keratan sulphate (KS), human natural killer trisaccharide (HNK-1) and chondroitin sulphate (CS). These convey unique tissue-specific functional properties in water imbibition, space-filling, matrix stabilisation or embryonic cellular regulation. Aggrecan also interacts with morphogens and growth factors directing tissue morphogenesis, remodelling and metaplasia. HNK-1 aggrecan glycoforms direct neural crest cell migration in embryonic development and is neuroprotective in perineuronal nets in the brain. The ability of the aggrecan core protein to assemble CS and KS chains at high density equips cartilage aggrecan with its well-known water-imbibing and weight-bearing properties. The importance of specific arrangements of GAG chains on aggrecan in all its forms is also a primary morphogenetic functional determinant providing aggrecan with unique tissue context dependent regulatory properties. The versatility displayed by aggrecan in biodiverse contexts is a function of its GAG side chains.