Dynamic 31 P-MRSI using spiral spectroscopic imaging can map mitochondrial capacity in muscles of the human calf during plantar flexion exercise at 7 T

3D ultra-short echo time 31P-MRSI with rosette k-space pattern: Feasibility and comparison with conventional weighted CSI

Phosphorus-31 magnetic resonance spectroscopic imaging (31P-MRSI) provides valuable non-invasive in vivo information on tissue metabolism but is burdened by poor sensitivity and prolonged scan duration. Ultra-short echo time (UTE) acquisitions minimize signal loss when probing signals with relatively short spin-spin relaxation time (T2), while also preventing first-order dephasing. Here, a three-dimensional (3D) UTE sequence with a rosette k-space trajectory is applied to 31P-MRSI at 3T. Conventional chemical shift imaging (CSI) employs highly regular Cartesian k-space sampling, susceptible to substantial artifacts when accelerated via undersampling. In contrast, this novel sequence's "petal-like" pattern offers incoherent sampling more suitable for compressed sensing (CS). These results showcase the competitive performance of UTE rosette 31P-MRSI against conventional weighted CSI with simulation, phantom, and in vivo leg muscle comparisons.

Reproducibility of 3D MRSI for imaging human brain glucose metabolism using direct (2H) and indirect (1H) detection of deuterium labeled compounds at 7T and clinical 3T

INTRODUCTION

Deuterium metabolic imaging (DMI) and quantitative exchange label turnover (QELT) are novel MR spectroscopy techniques for non-invasive imaging of human brain glucose and neurotransmitter metabolism with high clinical potential. Following oral or intravenous administration of non-ionizing [6,6'-2H2]-glucose, its uptake and synthesis of downstream metabolites can be mapped via direct or indirect detection of deuterium resonances using 2H MRSI (DMI) and 1H MRSI (QELT), respectively. The purpose of this study was to compare the dynamics of spatially resolved brain glucose metabolism, i.e., estimated concentration enrichment of deuterium labeled Glx (glutamate+glutamine) and Glc (glucose) acquired repeatedly in the same cohort of subjects using DMI at 7T and QELT at clinical 3T.

METHODS

Five volunteers (4 m/1f) were scanned in repeated sessions for 60 min after overnight fasting and 0.8 g/kg oral [6,6'-2H2]-glucose administration using time-resolved 3D 2H FID-MRSI with elliptical phase encoding at 7T and 3D 1H FID-MRSI with a non-Cartesian concentric ring trajectory readout at clinical 3T.

RESULTS

One hour after oral tracer administration regionally averaged deuterium labeled Glx4 concentrations and the dynamics were not significantly different over all participants between 7T 2H DMI and 3T 1H QELT data for GM (1.29±0.15 vs. 1.38±0.26 mM, p=0.65 & 21±3 vs. 26±3 µM/min, p=0.22) and WM (1.10±0.13 vs. 0.91±0.24 mM, p=0.34 & 19±2 vs. 17±3 µM/min, p=0.48). Also, the observed time constants of dynamic Glc6 data in GM (24±14 vs. 19±7 min, p=0.65) and WM (28±19 vs. 18±9 min, p=0.43) dominated regions showed no significant differences. Between individual 2H and 1H data points a weak to moderate negative correlation was observed for Glx4 concentrations in GM (r=-0.52, p<0.001), and WM (r=-0.3, p<0.001) dominated regions, while a strong negative correlation was observed for Glc6 data GM (r=-0.61, p<0.001) and WM (r=-0.70, p<0.001).

CONCLUSION

This study demonstrates that indirect detection of deuterium labeled compounds using 1H QELT MRSI at widely available clinical 3T without additional hardware is able to reproduce absolute concentration estimates of downstream glucose metabolites and the dynamics of glucose uptake compared to 2H DMI data acquired at 7T. This suggests significant potential for widespread application in clinical settings especially in environments with limited access to ultra-high field scanners and dedicated RF hardware.

Chemical exchange saturation transfer imaging of creatine, phosphocreatine, and protein arginine residue in tissues

Chemical exchange saturation transfer (CEST) MRI has become a promising technique to assay target proteins and metabolites through their exchangeable protons, noninvasively. The ubiquity of creatine (Cr) and phosphocreatine (PCr) due to their pivotal roles in energy homeostasis through the creatine phosphate pathway has made them prime targets for CEST in the diagnosis and monitoring of disease pathologies, particularly in tissues heavily dependent on the maintenance of rich energy reserves. Guanidinium CEST from protein arginine residues (i.e. arginine CEST) can also provide information about the protein profile in tissue. However, numerous obfuscating factors stand as obstacles to the specificity of arginine, Cr, and PCr imaging through CEST, such as semisolid magnetization transfer, fast chemical exchanges such as primary amines, and the effects of nuclear Overhauser enhancement from aromatic and amide protons. In this review, the specific exchange properties of protein arginine residues, Cr, and PCr, along with their validation, are discussed, including the considerations necessary to target and tune their signal effects through CEST imaging. Additionally, strategies that have been employed to enhance the specificity of these exchanges in CEST imaging are described, along with how they have opened up possible applications of protein arginine residues, Cr and PCr CEST imaging in the study and diagnosis of pathology. A clear understanding of the capabilities and caveats of using CEST to image these vital metabolites and mitigation strategies is crucial to expanding the possibilities of this promising technology.

Reproducibility of 3D MRSI for imaging human brain glucose metabolism using direct ( 2 H) and indirect ( 1 H) detection of deuterium labeled compounds at 7T and clinical 3T

Introduction

Deuterium metabolic imaging (DMI) and quantitative exchange label turnover (QELT) are novel MR spectroscopy techniques for non-invasive imaging of human brain glucose and neurotransmitter metabolism with high clinical potential. Following oral or intravenous administration of non-ionizing [6,6'- ² H ₂ ]-glucose, its uptake and synthesis of downstream metabolites can be mapped via direct or indirect detection of deuterium resonances using ² H MRSI (DMI) and ¹ H MRSI (QELT), respectively. The purpose of this study was to compare the dynamics of spatially resolved brain glucose metabolism, i.e., estimated concentration enrichment of deuterium labeled Glx (glutamate+glutamine) and Glc (glucose) acquired repeatedly in the same cohort of subjects using DMI at 7T and QELT at clinical 3T.

Methods

Five volunteers (4m/1f) were scanned in repeated sessions for 60 min after overnight fasting and 0.8g/kg oral [6,6'- ² H ₂ ]-glucose administration using time-resolved 3D ² H FID-MRSI with elliptical phase encoding at 7T and 3D ¹ H FID-MRSI with a non-Cartesian concentric ring trajectory readout at clinical 3T.

Results

One hour after oral tracer administration regionally averaged deuterium labeled Glx ₄ concentrations and the dynamics were not significantly different over all participants between 7T ² H DMI and 3T ¹ H QELT data for GM (1.29±0.15 vs. 1.38±0.26 mM, p=0.65 & 21±3 vs. 26±3 µM/min, p=0.22) and WM (1.10±0.13 vs. 0.91±0.24 mM, p=0.34 & 19±2 vs. 17±3 µM/min, p=0.48). Also, the observed time constants of dynamic Glc ₆ data in GM (24±14 vs. 19±7 min, p=0.65) and WM (28±19 vs. 18±9 min, p=0.43) dominated regions showed no significant differences. Between individual ² H and ¹ H data points a weak to moderate negative correlation was observed for Glx ₄ concentrations in GM (r=-0.52, p<0.001), and WM (r=-0.3, p<0.001) dominated regions, while a strong negative correlation was observed for Glc ₆ data GM (r=- 0.61, p<0.001) and WM (r=-0.70, p<0.001).

Conclusion

This study demonstrates that indirect detection of deuterium labeled compounds using ¹ H QELT MRSI at widely available clinical 3T without additional hardware is able to reproduce absolute concentration estimates of downstream glucose metabolites and the dynamics of glucose uptake compared to ² H DMI data acquired at 7T. This suggests significant potential for widespread application in clinical settings especially in environments with limited access to ultra-high field scanners and dedicated RF hardware.

MR-Nucleomics: The study of pathological cellular processes with multinuclear magnetic resonance spectroscopy and imaging in vivo

Clinical medicine has experienced a rapid development in recent decades, during which therapies targeting specific cellular signaling pathways, or specific cell surface receptors, have been increasingly adopted. While these developments in clinical medicine call for improved precision in diagnosis and treatment monitoring, modern medical imaging methods are restricted mainly to anatomical imaging, lagging behind the requirements of precision medicine. Although positron emission tomography and single photon emission computed tomography have been used clinically for studies of metabolism, their applications have been limited by the exposure risk to ionizing radiation, the subsequent limitation in repeated and longitudinal studies, and the incapability in assessing downstream metabolism. Magnetic resonance spectroscopy (MRS) or spectroscopic imaging (MRSI) are, in theory, capable of assessing molecular activities in vivo, although they are often limited by sensitivity. Here, we review some recent developments in MRS and MRSI of multiple nuclei that have potential as molecular imaging tools in the clinic.

In vivo magnetic resonance <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msup><mml:mrow/> <mml:mrow><mml:mn>31</mml:mn></mml:mrow> </mml:msup> </mml:mrow> </mml:math> P-Spectral Analysis With Neural Networks: 31P-SPAWNN

PURPOSE

We have introduced an artificial intelligence framework, 31P-SPAWNN, in order to fully analyze phosphorus-31 ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msup><mml:mrow/> <mml:mrow><mml:mn>31</mml:mn></mml:mrow> </mml:msup> </mml:mrow> <mml:annotation>$$ {}^{31} $$</mml:annotation></mml:semantics> </mml:math> P) magnetic resonance spectra. The flexibility and speed of the technique rival traditional least-square fitting methods, with the performance of the two approaches, are compared in this work.

THEORY AND METHODS

Convolutional neural network architectures have been proposed for the analysis and quantification of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msup><mml:mrow/> <mml:mrow><mml:mn>31</mml:mn></mml:mrow> </mml:msup> </mml:mrow> <mml:annotation>$$ {}^{31} $$</mml:annotation></mml:semantics> </mml:math> P-spectroscopy. The generation of training and test data using a fully parameterized model is presented herein. In vivo unlocalized free induction decay and three-dimensional <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msup><mml:mrow/> <mml:mrow><mml:mn>31</mml:mn></mml:mrow> </mml:msup> </mml:mrow> <mml:annotation>$$ {}^{31} $$</mml:annotation></mml:semantics> </mml:math> P-magnetic resonance spectroscopy imaging data were acquired from healthy volunteers before being quantified using either 31P-SPAWNN or traditional least-square fitting techniques.

RESULTS

The presented experiment has demonstrated both the reliability and accuracy of 31P-SPAWNN for estimating metabolite concentrations and spectral parameters. Simulated test data showed improved quantification using 31P-SPAWNN compared with LCModel. In vivo data analysis revealed higher accuracy at low signal-to-noise ratio using 31P-SPAWNN, yet with equivalent precision. Processing time using 31P-SPAWNN can be further shortened up to two orders of magnitude.

CONCLUSION

The accuracy, reliability, and computational speed of the method open new perspectives for integrating these applications in a clinical setting.

Emerging MRI techniques for molecular and functional phenotyping of the diseased heart

Recent advances in cardiac MRI (CMR) capabilities have truly transformed its potential for deep phenotyping of the diseased heart. Long known for its unparalleled soft tissue contrast and excellent depiction of three-dimensional (3D) structure, CMR now boasts a range of unique capabilities for probing disease at the tissue and molecular level. We can look beyond coronary vessel blockages and detect vessel disease not visible on a structural level. We can assess if early fibrotic tissue is being laid down in between viable cardiac muscle cells. We can measure deformation of the heart wall to determine early presentation of stiffening. We can even assess how cardiomyocytes are utilizing energy, where abnormalities are often precursors to overt structural and functional deficits. Finally, with artificial intelligence gaining traction due to the high computing power available today, deep learning has proven itself a viable contender with traditional acceleration techniques for real-time CMR. In this review, we will survey five key emerging MRI techniques that have the potential to transform the CMR clinic and permit early detection and intervention. The emerging areas are: (1) imaging microvascular dysfunction, (2) imaging fibrosis, (3) imaging strain, (4) imaging early metabolic changes, and (5) deep learning for acceleration. Through a concerted effort to develop and translate these areas into the CMR clinic, we are committing ourselves to actualizing early diagnostics for the most intractable heart disease phenotypes.

Interleaved and simultaneous multi-nuclear magnetic resonance in vivo. Review of principles, applications and potential

Magnetic resonance signals from different nuclei can be excited or received at the same time,rendering simultaneous or rapidly interleaved multi-nuclear acquisitions feasible. The advan-tages are a reduction of total scan time compared to sequential multi-nuclear acquisitions or that additional information from heteronuclear data is obtained at thesame time and anatomical position. Information content can be qualitatively increased by delivering a more comprehensive MR-based picture of a transient state (such as an exercise bout). Also, combiningnon-proton MR acquisitions with ¹ Hinformation (e.g., dynamic shim updates and motion correction) can be used to improve data quality during long scans and benefits image coregistration. This work reviews the literature on interleaved and simultaneous multi-nuclear MRI and MRS in vivo. Prominent use cases for this methodology in clinical and research applications are brain and muscle, but studies have also been carried out in other targets, including the lung, knee, breast and heart. Simultaneous multi-nuclear measurements in the liver and kidney have also been performed, but exclusively in rodents. In this review, a consistent nomenclature is proposed, to help clarify the terminology used for this principle throughout the literature on in-vivo MR. An overview covers the basic principles, the technical requirements on the MR scanner and the implementations realised either by MR system vendors or research groups, from the early days until today. Considerations regarding the multi-tuned RF coils required and heteronuclear polarisation interactions are briefly discussed, and fields for future in-vivo applications for interleaved multi-nuclear MR pulse sequences are identified.

Achieving high-resolution 1H-MRSI of the human brain with compressed-sensing and low-rank reconstruction at 7 Tesla

Low sensitivity MR techniques such as magnetic resonance spectroscopic imaging (MRSI) greatly benefit from the gain in signal-to-noise provided by ultra-high field MR. High-resolution and whole-slab brain MRSI remains however very challenging due to lengthy acquisition, low signal, lipid contamination and field inhomogeneity. In this study, we propose an acquisition-reconstruction scheme that combines 1H free-induction-decay (FID)-MRSI sequence, short TR acquisition, compressed sensing acceleration and low-rank modeling with total-generalized-variation constraint to achieve metabolite imaging in two and three dimensions at 7 Tesla. The resulting images and volumes reveal highly detailed distributions that are specific to each metabolite and follow the underlying brain anatomy. The MRSI method was validated in a high-resolution phantom containing fine metabolite structures, and in five healthy volunteers. This new application of compressed sensing acceleration paves the way for high-resolution MRSI in clinical setting with acquisition times of 5 min for 2D MRSI at 2.5 mm and of 20 min for 3D MRSI at 3.3 mm isotropic.

Accelerated MR spectroscopic imaging-a review of current and emerging techniques

Over more than 30 years in vivo MR spectroscopic imaging (MRSI) has undergone an enormous evolution from theoretical concepts in the early 1980s to the robust imaging technique that it is today. The development of both fast and efficient sampling and reconstruction techniques has played a fundamental role in this process. State-of-the-art MRSI has grown from a slow purely phase-encoded acquisition technique to a method that today combines the benefits of different acceleration techniques. These include shortening of repetition times, spatial-spectral encoding, undersampling of k-space and time domain, and use of spatial-spectral prior knowledge in the reconstruction. In this way in vivo MRSI has considerably advanced in terms of spatial coverage, spatial resolution, acquisition speed, artifact suppression, number of detectable metabolites and quantification precision. Acceleration not only has been the enabling factor in high-resolution whole-brain 1 H-MRSI, but today is also common in non-proton MRSI (31 P, 2 H and 13 C) and applied in many different organs. In this process, MRSI techniques had to constantly adapt, but have also benefitted from the significant increase of magnetic field strength boosting the signal-to-noise ratio along with high gradient fidelity and high-density receive arrays. In combination with recent trends in image reconstruction and much improved computation power, these advances led to a number of novel developments with respect to MRSI acceleration. Today MRSI allows for non-invasive and non-ionizing mapping of the spatial distribution of various metabolites' tissue concentrations in animals or humans, is applied for clinical diagnostics and has been established as an important tool for neuro-scientific and metabolism research. This review highlights the developments of the last five years and puts them into the context of earlier MRSI acceleration techniques. In addition to 1 H-MRSI it also includes other relevant nuclei and is not limited to certain body regions or specific applications.

Atlas-based GABA mapping with 3D MEGA-MRSI: Cross-correlation to single-voxel MRS

The purpose of this work is to develop and validate a new atlas-based metabolite quantification pipeline for edited magnetic resonance spectroscopic imaging (MEGA-MRSI) that enables group comparisons of brain structure-specific GABA levels. By using brain structure masks segmented from high-resolution MPRAGE images and coregistering these to MEGA-LASER 3D MRSI data, an automated regional quantification of neurochemical levels is demonstrated for the example of the thalamus. Thalamic gamma-aminobutyric acid + coedited macromolecules (GABA+) levels from 21 healthy subjects scanned at 3 T were cross-validated both against a single-voxel MEGA-PRESS acquisition in the same subjects and same scan sessions, as well as alternative MRSI processing techniques (ROI approach, four-voxel approach) using Pearson correlation analysis. In addition, reproducibility was compared across the MRSI processing techniques in test-retest data from 14 subjects. The atlas-based approach showed a significant correlation with SV MEGA-PRESS (correlation coefficient r [GABA+] = 0.63, P < 0.0001). However, the actual values for GABA+, NAA, tCr, GABA+/tCr and tNAA/tCr obtained from the atlas-based approach showed an offset to SV MEGA-PRESS levels, likely due to the fact that on average the thalamus mask used for the atlas-based approach only occupied 30% of the SVS volume, ie, somewhat different anatomies were sampled. Furthermore, the new atlas-based approach showed highly reproducible GABA+/tCr values with a low median coefficient of variance of 6.3%. In conclusion, the atlas-based metabolite quantification approach enables a more brain structure-specific comparison of GABA+ and other neurochemical levels across populations, even when using an MRSI technique with only cm-level resolution. This approach was successfully cross-validated against the typically used SVS technique as well as other different MRSI analysis methods, indicating the robustness of this quantification approach.

Utilizing Dynamic Phosphorous-31 Magnetic Resonance Spectroscopy for the Early Detection of Acute Compartment Syndrome: A Pilot Study on Rats

INTRODUCTION

Disasters, including terrorism and earthquakes, are significant threats to people and may lead to many people requiring rescue. The longer the rescue takes, the higher the chances of an individual contracting acute compartment syndrome (ACS). ACS is fatal if diagnosed too late, and early diagnosis and treatment are essential.

OBJECTIVE

To assess the ability of dynamic phosphorus magnetic resonance spectroscopy (31P-MRS) in the early detection of muscular damage in ACS.

MATERIALS AND METHODS

Six ACS model rats were used for serial 31P-MRS scanning (9.4 Tesla). Skeletal muscle metabolism, represented by the levels of phosphocreatine (PCr), inorganic phosphate (Pi), and adenosine triphosphate (ATP), was assessed. The PCr/(Pi + PCr) ratio, which decreases with ischemia, was compared with simultaneously sampled plasma creatine phosphokinase (CPK), a muscle damage marker.

RESULTS

The PCr/(Pi + PCr) ratio significantly decreased after inducing ischemia (from 0.86 ± 0.10 to 0.18 ± 0.06; p < 0.05), while CPK did not change significantly (from 89 ± 29.46 to 241.50 ± 113.28; p > 0.05). The intracellular and arterial pH index decreased over time, revealing significant differences at 120 min post-ischemia (from 7.09 ± 0.01 to 6.43 ± 0.13, and from 7.47 ± 0.03 to 7.39 ± 0.04, respectively). In the reperfusion state, the spectra and pH did not return to the original values.

CONCLUSIONS

The dynamic 31P-MRS technique can rapidly detect changes in muscle bioenergetics. This technique is a promising non-invasive method for determining early muscular damage in ACS.

High-Resolution Dynamic 31P-MR Spectroscopic Imaging for Mapping Mitochondrial Function

OBJECTIVE

To enable non-invasive dynamic metabolic mapping in rodent model studies of mitochondrial function using 31P-MR spectroscopic imaging (MRSI).

METHODS

We developed a novel method for high-resolution dynamic 31P-MRSI. The method synergistically integrates physics-based models of spectral structures, biochemical modeling of molecular dynamics, and subspace learning to capture spatiospectral variations. Fast data acquisition was achieved using rapid spiral trajectories and sparse sampling of (k, t, T)-space; image reconstruction was accomplished using a low-rank tensor-based framework.

RESULTS

The proposed method provided high-resolution dynamic metabolic mapping in rat hindlimb at spatial and temporal resolutions of 4[Formula: see text]2 mm3 and 1.28 s, respectively. This allowed for in vivo mapping of the time-constant of phosphocreatine resynthesis, a well established index of mitochondrial oxidative capacity. Multiple rounds of in vivo experiments were performed to demonstrate reproducibility, and in vitro experiments were used to validate the accuracy of the estimated metabolite maps.

CONCLUSIONS

A new model-based method is proposed to achieve high-resolution dynamic 31P-MRSI. The proposed method's ability to delineate metabolic heterogeneity was demonstrated in rat hindlimb.

SIGNIFICANCE

Abnormal mitochondrial metabolism is a key cellular dysfunction in many prevalent diseases such as diabetes and heart disease; however, current understanding of mitochondrial function is mostly gained from studies on isolated mitochondria under nonphysiological conditions. The proposed method has the potential to open new avenues of research by allowing in vivo and longitudinal studies of mitochondrial dysfunction in disease development and progression.

In vivo imaging of phosphocreatine with artificial neural networks

Phosphocreatine (PCr) plays a vital role in neuron and myocyte energy homeostasis. Currently, there are no routine diagnostic tests to noninvasively map PCr distribution with clinically relevant spatial resolution and scan time. Here, we demonstrate that artificial neural network-based chemical exchange saturation transfer (ANNCEST) can be used to rapidly quantify PCr concentration with robust immunity to commonly seen MRI interferences. High-quality PCr mapping of human skeletal muscle, as well as the information of exchange rate, magnetic field and radio-frequency transmission inhomogeneities, can be obtained within 1.5 min on a 3 T standard MRI scanner using ANNCEST. For further validation, we apply ANNCEST to measure the PCr concentrations in exercised skeletal muscle. The ANNCEST outcomes strongly correlate with those from 31P magnetic resonance spectroscopy (R = 0.813, p < 0.001, t test). These results suggest that ANNCEST has potential as a cost-effective and widely available method for measuring PCr and diagnosing related diseases.

Interleaved 31 P MRS/1 H ASL for analysis of metabolic and functional heterogeneity along human lower leg muscles at 7T

PURPOSE

MR offers the unique possibility to noninvasively investigate cellular energy metabolism via 31P MRS, while blood perfusion, which provides oxygen and substrates to the tissue, is accessible by arterial spin labeling (ASL) 1H MRI. Because metabolic and hemodynamic parameters are linked, it would be desirable to study them simultaneously. A 3D-resolved method is presented that allows such measurements with high spatiotemporal resolution and has the potential to discern differences along an exercising muscle.

METHODS

Multi-voxel localized ³¹ P MRS was temporally interleaved with multi-slice pASL 1H MRI. Phosphorus spectra were collected from two adjacent positions in gastrocnemius medialis (GM) during rest, submaximal plantar flexion exercise and recovery, while perfusion and T 2 * -weighted axial images were acquired at the same time. Seventeen healthy volunteers (9 f / 8 m) were studied at 7 T.

RESULTS

An increase of postexercise perfusion and T 2 * -weighted signal in GM positively correlated with end-exercise PCr depletion and pH drop. At proximal positions functional and metabolic activity was higher than distally, that is, perfusion increase and peak T 2 * -weighted signal, end-exercise PCr depletion, end-exercise pH, and PCr recovery time constant were significantly different. An NOE-induced SNR increase of approximately 20 % (P < .001), at rest, was found in interleaved ³¹ P spectra, when comparing to ³¹ P-only acquisitions.

CONCLUSIONS

A technique for fast, simultaneous imaging of muscle functional heterogeneity in ASL, T 2 * and acquisition of time-resolved ³¹ P  MRS data is presented. These single exercise recovery experiments can be used to investigate local variations during disease progression in patients suffering from vascular or muscular diseases.

Magnetic Resonance Imaging of the Musculoskeletal System at 7T: Morphological Imaging and Beyond

In 2017, a whole-body 7T magnetic resonance imaging (MRI) device was given regulatory approval for clinical use in both the EU and United States for neuro and musculoskeletal applications. As 7 Tesla allows for higher signal-to-noise , which results in higher resolution images than those obtained on lower-field-strength scanners, it has attracted considerable attention from the musculoskeletal field, as evidenced by the increasing number of publications in the last decade. Besides morphological imaging, the quantitative MR methods, such as T₂, T₂∗, T_1ρ mapping, sodium imaging, chemical-exchange saturation transfer, and spectroscopy, substantially benefit from ultrahigh field scanning. In this review, we provide technical considerations for the individual techniques and an overview of (mostly) clinical applications for the assessment of cartilage, tendon, meniscus, and muscle. The first part of the review is dedicated to morphological applications at 7T, and the second part describes the most recent developments in quantitative MRI at 7T.

MR compatible ergometers for dynamic 31P MRS

Magnetic Resonance (MR) compatible ergometers are specialized ergometers used inside the MR scanners for the characterization of tissue metabolism changes during physical stress. They are most commonly used for dynamic phosphorous magnetic resonance spectroscopy (31P MRS), but can also be used for lactate production measurements, perfusion studies using arterial spin labelling or muscle oxygenation measurements by blood oxygen dependent contrast sequences. We will primarily discuss the importance of ergometers in the context of dynamic 31P MRS. Dynamic 31P MRS can monitor muscle fatigue and energy reserve during muscle contractions as well as the dynamics of recuperation of skeletal muscle tissue during the following recovery through signal changes of phosphocreatine (PCr), inorganic phosphate and adenosine triphosphate (ATP). Based on the measured data it is possible to calculate intracellular pH, metabolic flux of ATP through creatine-kinase reaction, anaerobic glycolysis and oxidative phosphorylation and other metabolic parameters as mitochondrial capacity. This review primarily focuses on describing various technical designs of MR compatible ergometers for dynamic 31P MRS that must be constructed with respect to the presence of magnetic field. It is also expected that the construction of ergometers will be easy for the handling and well accepted by examined subjects.

Comparison of compressed sensing reconstruction algorithms for 31P magnetic resonance spectroscopic imaging

Phosphorus MR spectroscopy and spectroscopic imaging (³¹P-MRS/MRSI) provide information about energy metabolism, membrane degradation and pH in vivo. In spite of their proven utility, ³¹P-MRS/MRSI are not often used primarily because of the challenges imposed by the low sensitivity and low concentration of metabolites leading to low signal to noise ratio (SNR), coarse spatial resolution and prolonged acquisition time. More recently there has been considerable interest in compressed sensing as an acceleration method for MR signal acquisition. This approach takes advantage of the intrinsic sparsity of the spectral data. In this work, we present a ³¹P-MRSI sequence that combines a flyback EPSI trajectory and compressed sensing, and we compared two different reconstruction methods, L1 norm minimization and low rank Hankel matrix completion. Our phantom results showed good preservation of spectral quality for both ×2.0 and ×3.0 acceleration factors, using both CS reconstruction methods. However, in vivo ³¹P-MRS brain data showed the low rank reconstruction approach was most suitable. Overall, this study shows the feasibility of combining a flyback EPSI trajectory and compressed sensing in the acquisition of ³¹P-MRSI as well as the better suitability of a low rank reconstruction approach.

Dynamic 31 P spectroscopic imaging of skeletal muscles combining flyback echo-planar spectroscopic imaging and compressed sensing

PURPOSE

Dynamic phosphorus MR spectroscopic imaging (³¹ P-MRSI) experiments require temporal resolution on the order of seconds to concurrently assess different muscle groups. A highly accelerated pulse sequence combining flyback echo-planar spectroscopic imaging (EPSI) and compressed sensing was developed and tested in a phantom and healthy humans during an exercise-recovery challenge of the lower leg muscles, using a clinical 3T MRI.

METHODS

A flyback EPSI readout designed to achieve 2.25 × 2.25 cm² resolution over a 18 × 18 cm² field of view (i.e., 8 × 8 matrix) was combined with compressed sensing through the inclusion of pseudorandom gradient blips to sub-sample the ky-kt dimensions by a factor of 2.7×, achieving a temporal resolution of 9 s. The sequence was first tested in a phantom to assess performance compared to fully sampled EPSI (fidEPSI) and phase encoded chemical shift imaging (fidCSI). Then, tests were performed in 11 healthy volunteers during an exercise-recovery challenge of the lower leg muscles. Voxels containing tissue from different muscle groups were evaluated measuring percentage phosphocreatine (%PCr) depletion, time constant of PCr recovery (τPCr) and intracellular pH at rest and following exercise.

RESULTS

The sequence was capable to track the dynamic PCr response of multiple muscles simultaneously. No statistical differences were found in the metabolite ratio, pH or linewidth when compared with fidEPSI and fidCSI in the phantom study. Dynamic experiments showed differences in PCr depletion when comparing soleus with gastrocnemius muscles. Intracellular pH, τPCr and %PCr decrease were consistent with reported values.

CONCLUSION

Highly accelerated ³¹ P-MRSI combining flyback EPSI and compressed sensing is capable of assessing concurrent energy metabolism in multiple muscle groups using a clinical 3T MR system.

Pros and cons of ultra-high-field MRI/MRS for human application

Magnetic resonance imaging and spectroscopic techniques are widely used in humans both for clinical diagnostic applications and in basic research areas such as cognitive neuroimaging. In recent years, new human MR systems have become available operating at static magnetic fields of 7 T or higher (≥300 MHz proton frequency). Imaging human-sized objects at such high frequencies presents several challenges including non-uniform radiofrequency fields, enhanced susceptibility artifacts, and higher radiofrequency energy deposition in the tissue. On the other side of the scale are gains in signal-to-noise or contrast-to-noise ratio that allow finer structures to be visualized and smaller physiological effects to be detected. This review presents an overview of some of the latest methodological developments in human ultra-high field MRI/MRS as well as associated clinical and scientific applications. Emphasis is given to techniques that particularly benefit from the changing physical characteristics at high magnetic fields, including susceptibility-weighted imaging and phase-contrast techniques, imaging with X-nuclei, MR spectroscopy, CEST imaging, as well as functional MRI. In addition, more general methodological developments such as parallel transmission and motion correction will be discussed that are required to leverage the full potential of higher magnetic fields, and an overview of relevant physiological considerations of human high magnetic field exposure is provided.

Dynamic multivoxel-localized 31 P MRS during plantar flexion exercise with variable knee angle

Exercise studies investigating the metabolic response of calf muscles using ³¹ P MRS are usually performed with a single knee angle. However, during natural movement, the distribution of workload between the main contributors to force, gastrocnemius and soleus is influenced by the knee angle. Hence, it is of interest to measure the respective metabolic response of these muscles to exercise as a function of knee angle using localized spectroscopy. Time-resolved multivoxel ³¹ P MRS at 7 T was performed simultaneously in gastrocnemius medialis and soleus during rest, plantar flexion exercise and recovery in 12 healthy volunteers. This experiment was conducted with four different knee angles. PCr depletions correlated negatively with knee angle in gastrocnemius medialis, decreasing from 79±14 % (extended leg) to 35±23 %(∼40°), and positively in soleus, increasing from 20±21 % to 36±25 %; differences were significant. Linear correlations were found between knee angle and end-exercise PCr depletions in gastrocnemius medialis (R² =0.8) and soleus (R² =0.53). PCr recovery times and end-exercise pH changes that correlated with PCr depletion were consistent with the literature in gastrocnemius medialis and differences between knee angles were significant. These effects were less pronounced in soleus and not significant for comparable PCr depletions. Maximum oxidative capacity calculated for all knee angles was in excellent agreement with the literature and showed no significant changes between different knee angles. In conclusion, these findings confirm that plantar flexion exercise with a straight leg is a suitable paradigm, when data are acquired from gastrocnemius only (using either localized MRS or small surface coils), and that activation of soleus requires the knee to be flexed. The present study comprises a systematic investigation of the effects of the knee angle on metabolic parameters, measured with dynamic multivoxel ³¹ P MRS during muscle exercise and recovery, and the findings should be used in future study design.

Oxidative capacity varies along the length of healthy human tibialis anterior

KEY POINTS

During exercise skeletal muscles use the energy buffer phosphocreatine. The post-exercise recovery of phosphocreatine is a measure of the oxidative capacity of muscles and is traditionally assessed by ³¹ P magnetic resonance spectroscopy of a large tissue region, assuming homogeneous energy metabolism. To test this assumption, we collected spatially resolved spectra along the length of human tibialis anterior using a home-built array of ³¹ P detection coils, and observed a striking gradient in the recovery rate of phosphocreatine, decreasing along the proximo-distal axis of the muscle. A similar gradient along this muscle was observed in signal changes recorded by ¹ H muscle functional MRI. These findings identify intra-muscular variation in the physiology of muscles in action and highlight the importance of localized sampling for any methodology investigating oxidative metabolism of this, and potentially other muscles.

ABSTRACT

The rate of phosphocreatine (PCr) recovery (k_PCr ) after exercise, characterizing muscle oxidative capacity, is traditionally assessed with unlocalized ³¹ P magnetic resonance spectroscopy (MRS) using a single surface coil. However, because of intramuscular variation in fibre type and oxygen supply, k_PCr may be non-uniform within muscles. We tested this along the length of the tibialis anterior (TA) muscle in 10 male volunteers. For this purpose, we employed a 3T MR system with a ³¹ P/¹ H volume transmit coil combined with a home-built ³¹ P phased-array receive probe, consisting of five coil elements covering the TA muscle length. Mono-exponential k_PCr was determined for all coil elements after 40 s of submaximal isometric dorsiflexion (SUBMAX) and incremental exercise to exhaustion (EXH). In addition, muscle functional MRI (¹ H mfMRI) was performed using the volume coil after another 40 s of SUBMAX. A strong gradient in k_PCr was observed along the TA (P < 0.001), being two times higher proximally vs. distally during SUBMAX and EXH. Statistical analysis showed that this gradient cannot be explained by pH variations. A similar gradient was seen in the slope of the initial post-exercise ¹ H mfMRI signal change, which was higher proximally than distally in both the TA and the extensor digitorum longus (P < 0.001) and strongly correlated with k_PCr . The pronounced differences along the TA in functional oxidative capacity identify regional variation in the physiological demand of this muscle during everyday activities and have implications for the bio-energetic assessment of interventions to modify its performance and of neuromuscular disorders involving the TA.

Phosphorus magnetic resonance spectroscopic imaging using flyback echo planar readout trajectories

OBJECT

To present and evaluate a fast phosphorus magnetic resonance spectroscopic imaging (MRSI) sequence using echo planar spectroscopic imaging with flyback readout gradient trajectories.

MATERIALS AND METHODS

Waveforms were designed and implemented using a 3 Tesla MRI system. ³¹P spectra were acquired with 2 × 2 cm² and 3 × 3 cm² resolution over a 20- and 21-cm field of view and spectral bandwidths up to 1923 Hz. The sequence was first tested using a 20-cm-diameter phosphate phantom, and subsequent in vivo tests were performed on healthy human calf muscles and brains from five volunteers.

RESULTS

Flyback EPSI achieved 10× and 7× reductions in acquisition time, with 68.0 ± 1.2 and 69.8 ± 2.2% signal-to-noise ratio (SNR) per unit of time efficiency (theoretical SNR efficiency was 74.5 and 76.4%) for the in vivo experiments, compared to conventional phase-encoded MRSI for the 2 × 2 cm² and 3 × 3 cm² resolution waveforms, respectively. Statistical analysis showed no difference in the quantification of most metabolites. Time savings and SNR comparisons were consistent across phantom, leg and brain experiments.

CONCLUSION

EPSI using flyback readout trajectories was found to be a reliable alternative for acquiring ³¹P-MRSI data in a shorter acquisition time.

Density-weighted concentric rings k-space trajectory for 1 H magnetic resonance spectroscopic imaging at 7 T

It has been shown that density-weighted (DW) k-space sampling with spiral and conventional phase encoding trajectories reduces spatial side lobes in magnetic resonance spectroscopic imaging (MRSI). In this study, we propose a new concentric ring trajectory (CRT) for DW-MRSI that samples k-space with a density that is proportional to a spatial, isotropic Hanning window. The properties of two different DW-CRTs were compared against a radially equidistant (RE) CRT and an echo-planar spectroscopic imaging (EPSI) trajectory in simulations, phantoms and in vivo experiments. These experiments, conducted at 7 T with a fixed nominal voxel size and matched acquisition times, revealed that the two DW-CRT designs improved the shape of the spatial response function by suppressing side lobes, also resulting in improved signal-to-noise ratio (SNR). High-quality spectra were acquired for all trajectories from a specific region of interest in the motor cortex with an in-plane resolution of 7.5 × 7.5 mm2 in 8 min 3 s. Due to hardware limitations, high-spatial-resolution spectra with an in-plane resolution of 5 × 5 mm2 and an acquisition time of 12 min 48 s were acquired only for the RE and one of the DW-CRT trajectories and not for EPSI. For all phantom and in vivo experiments, DW-CRTs resulted in the highest SNR. The achieved in vivo spectral quality of the DW-CRT method allowed for reliable metabolic mapping of eight metabolites including N-acetylaspartylglutamate, γ-aminobutyric acid and glutathione with Cramér-Rao lower bounds below 50%, using an LCModel analysis. Finally, high-quality metabolic mapping of a whole brain slice using DW-CRT was achieved with a high in-plane resolution of 5 × 5 mm2 in a healthy subject. These findings demonstrate that our DW-CRT MRSI technique can perform robustly on MRI systems and within a clinically feasible acquisition time.

Assessing tissue metabolism by phosphorous-31 magnetic resonance spectroscopy and imaging: a methodology review

Many human diseases are caused by an imbalance between energy production and demand. Magnetic resonance spectroscopy (MRS) and magnetic resonance imaging (MRI) provide the unique opportunity for in vivo assessment of several fundamental events in tissue metabolism without the use of ionizing radiation. Of particular interest, phosphate metabolites that are involved in ATP generation and utilization can be quantified noninvasively by phosphorous-31 (31P) MRS/MRI. Furthermore, 31P magnetization transfer (MT) techniques allow in vivo measurement of metabolic fluxes via creatine kinase (CK) and ATP synthase. However, a major impediment for the clinical applications of 31P-MRS/MRI is the prohibitively long acquisition time and/or the low spatial resolution that are necessary to achieve adequate signal-to-noise ratio. In this review, current 31P-MRS/MRI techniques used in basic science and clinical research are presented. Recent advances in the development of fast 31P-MRS/MRI methods are also discussed.

OXSA: An open-source magnetic resonance spectroscopy analysis toolbox in MATLAB

In vivo magnetic resonance spectroscopy provides insight into metabolism in the human body. New acquisition protocols are often proposed to improve the quality or efficiency of data collection. Processing pipelines must also be developed to use these data optimally. Current fitting software is either targeted at general spectroscopy fitting, or for specific protocols. We therefore introduce the MATLAB-based OXford Spectroscopy Analysis (OXSA) toolbox to allow researchers to rapidly develop their own customised processing pipelines. The toolbox aims to simplify development by: being easy to install and use; seamlessly importing Siemens Digital Imaging and Communications in Medicine (DICOM) standard data; allowing visualisation of spectroscopy data; offering a robust fitting routine; flexibly specifying prior knowledge when fitting; and allowing batch processing of spectra. This article demonstrates how each of these criteria have been fulfilled, and gives technical details about the implementation in MATLAB. The code is freely available to download from https://github.com/oxsatoolbox/oxsa.

13C/31P MRS Metabolic Biomarkers of Disease Progression and Response to AAV Delivery of hGAA in a Mouse Model of Pompe Disease

The development of therapeutic clinical trials for glycogen storage disorders, including Pompe disease, has called for non-invasive and objective biomarkers. Glycogen accumulation can be measured in vivo with ¹³C MRS. However, clinical implementation remains challenging due to low signal-to-noise. On the other hand, the buildup of glycolytic intermediates may be detected with ³¹P MRS. We sought to identify new biomarkers of disease progression in muscle using ¹³C/³¹P MRS and ¹H HR-MAS in a mouse model of Pompe disease (Gaa^−/−). We evaluated the sensitivity of these MR biomarkers in vivo after treatment using an adeno-associated virus vector 2/9 encoding hGAA driven by the desmin promotor. ³¹P MRS showed significantly elevated phosphomonoesters (PMEs) in Gaa^−/− compared to control at 2 (0.06 ± 0.02 versus 0.03 ± 0.01; p = 0.003), 6, 12, and 18 months of age. Correlative ¹H HR-MAS measures in intact gastrocnemius muscles revealed high glucose-6-phosphate (G-6-P). After intramuscular AAV injections, glycogen, PME, and G-6-P were decreased within normal range. The changes in PME levels likely partly resulted from changes in G-6-P, one of the overlapping phosphomonoesters in the ³¹P MR spectra in vivo. Because ³¹P MRS is inherently more sensitive than ¹³C MRS, PME levels have greater potential as a clinical biomarker and should be considered as a complementary approach for future studies in Pompe patients.

Dynamic phosphocreatine imaging with unlocalized pH assessment of the human lower leg muscle following exercise at 3T

PURPOSE

To develop a high temporal resolution imaging method that measures muscle-specific phosphocreatine (PCr) resynthesis time constant (τ_PCr ) and pH changes in muscles of the lower leg following exercise on a clinical 3T MRI scanner.

METHODS

We developed a frequency-selective 3D non-Cartesian FLORET sequence to measure PCr with 17-mm nominal isotropic resolution (28 mm actual resolution) and 6-s temporal resolution to capture dynamic metabolic muscle activity. The sequence was designed to additionally collect inorganic phosphate spectra for pH quantification, which were localized using sensitivity profiles of individual coil elements. Nineteen healthy volunteers were scanned while performing a plantar flexion exercise on an in-house developed ergometer. Data were acquired with a dual-tuned multichannel coil array that enabled phosphorus imaging and proton localization for muscle segmentation.

RESULTS

After a 90-s plantar flexion exercise at 0.66 Hz with resistance set to 40% of the maximum voluntary contraction, τ_PCr was estimated at 22.9 ± 8.8 s (mean ± standard deviation) with statistical coefficient of determination r²  = 0.89 ± 0.05. The corresponding pH values after exercise were in the range of 6.9-7.1 in the gastrocnemius muscle.

CONCLUSION

The developed technique allows measurement of muscle-specific PCr resynthesis kinetics and pH changes following exercise, with a temporal resolution and accuracy comparable to that of single voxel ³¹ P-MRS sequences. Magn Reson Med 79:974-980, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

High-resolution dynamic 31 P-MRSI using a low-rank tensor model

PURPOSE

To develop a rapid ³¹ P-MRSI method with high spatiospectral resolution using low-rank tensor-based data acquisition and image reconstruction.

METHODS

The multidimensional image function of ³¹ P-MRSI is represented by a low-rank tensor to capture the spatial-spectral-temporal correlations of data. A hybrid data acquisition scheme is used for sparse sampling, which consists of a set of "training" data with limited k-space coverage to capture the subspace structure of the image function, and a set of sparsely sampled "imaging" data for high-resolution image reconstruction. An explicit subspace pursuit approach is used for image reconstruction, which estimates the bases of the subspace from the "training" data and then reconstructs a high-resolution image function from the "imaging" data.

RESULTS

We have validated the feasibility of the proposed method using phantom and in vivo studies on a 3T whole-body scanner and a 9.4T preclinical scanner. The proposed method produced high-resolution static ³¹ P-MRSI images (i.e., 6.9 × 6.9 × 10 mm³ nominal resolution in a 15-min acquisition at 3T) and high-resolution, high-frame-rate dynamic ³¹ P-MRSI images (i.e., 1.5 × 1.5 × 1.6 mm³ nominal resolution, 30 s/frame at 9.4T).

CONCLUSIONS

Dynamic spatiospectral variations of ³¹ P-MRSI signals can be efficiently represented by a low-rank tensor. Exploiting this mathematical structure for data acquisition and image reconstruction can lead to fast ³¹ P-MRSI with high resolution, frame-rate, and SNR. Magn Reson Med 78:419-428, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

In-vivo31P-MRS of skeletal muscle and liver: A way for non-invasive assessment of their metabolism

In addition to direct assessment of high energy phosphorus containing metabolite content within tissues, phosphorus magnetic resonance spectroscopy (³¹P-MRS) provides options to measure phospholipid metabolites and cellular pH, as well as the kinetics of chemical reactions of energy metabolism in vivo. Even though the great potential of ³¹P-MR was recognized over 30 years ago, modern MR systems, as well as new, dedicated hardware and measurement techniques provide further opportunities for research of human biochemistry. This paper presents a methodological overview of the ³¹P-MR techniques that can be used for basic, physiological, or clinical research of human skeletal muscle and liver in vivo. Practical issues of ³¹P-MRS experiments and examples of potential applications are also provided. As signal localization is essential for liver ³¹P-MRS and is important for dynamic muscle examinations as well, typical localization strategies for ³¹P-MR are also described.