Assessment of (31)P relaxation times in the human calf muscle: a comparison between 3 T and 7 T in vivo

Robust dual-angle T 1 $$ {T}_1 $$ measurement in magnetization transfer spectroscopy by time-optimal control

Magnetization transfer spectroscopy relies heavily on the robust determination ofT 1 $$ {T}_1 $$ relaxation times of nuclei participating in metabolic exchange. Challenges arise due to the use of surface RF coils for transmission (highB 1 + $$ {B}_1^{+} $$ variation) and the broad resonance band of most X nuclei. These challenges are particularly pronounced when fastT 1 $$ {T}_1 $$ mapping methods, such as the dual-angle method, are employed. Consequently, in this work, we develop resonance offset andB 1 + $$ {B}_1^{+} $$ robust excitation RF pulses for 31P magnetization transfer spectroscopy at 7T through ensemble-based time-optimal control. In our approach, we introduce a cost functional for designing robust pulses, incorporating the full Bloch equations as constraints, which are solved using symmetric operator splitting techniques. The optimal control design of the RF pulses developed demonstrates improved accuracy, desired phase properties, and reduced RF power when applied to dual-angleT 1 $$ {T}_1 $$ mapping, thereby improving the precision of exchange-rate measurements, as demonstrated in a preclinical in vivo study quantifying brain creatine kinase activity.

Fast in vivo assay of creatine kinase activity in the human brain by 31 P magnetic resonance fingerprinting

A new and efficient magnetisation transfer 31 P magnetic resonance fingerprinting (MT-31 P-MRF) approach is introduced to measure the creatine kinase metabolic ratek CK between phosphocreatine (PCr) and adenosine triphosphate (ATP) in human brain. The MRF framework is extended to overcome challenges in conventional 31 P measurement methods in the human brain, enabling reduced acquisition time and specific absorption rate (SAR). To address the challenge of creating and matching large multiparametric dictionaries in an MRF scheme, a nested iteration interpolation method (NIIM) is introduced. As the number of parameters to estimate increases, the size of the dictionary grows exponentially. NIIM can reduce the computational load by breaking dictionary matching into subsolutions of linear computational order. MT-31 P-MRF combined with NIIM providesT 1 PCr ,T 1 ATP andk CK estimates in good agreement with those obtained by the exchange kinetics by band inversion transfer (EBIT) method and literature values. Furthermore, the test-retest reproducibility results showed that MT-31 P-MRF achieves a similar or better coefficient of variation (<12%) forT 1 ATP andk CK measurements in 4 min 15 s, than EBIT with 17 min 4 s scan time, enabling a fourfold reduction in scan time. We conclude that MT-31 P-MRF in combination with NIIM is a fast, accurate, and reproducible approach for in vivok CK assays in the human brain, which enables the potential to investigate energy metabolism in a clinical setting.

In vivo phosphorus magnetic resonance spectroscopic imaging of the whole human liver at 7 T using a phosphorus whole-body transmit coil and 16-channel receive array: Repeatability and effects of principal component analysis-based denoising

Quantitative three-dimensional (3D) imaging of phosphorus (³¹ P) metabolites is potentially a promising technique with which to assess the progression of liver disease and monitor therapy response. However, ³¹ P magnetic resonance spectroscopy has a low sensitivity and commonly used ³¹ P surface coils do not provide full coverage of the liver. This study aimed to overcome these limitations by using a ³¹ P whole-body transmit coil in combination with a 16-channel ³¹ P receive array at 7 T. Using this setup, we determined the repeatability of whole-liver ³¹ P magnetic resonance spectroscopic imaging (³¹ P MRSI) in healthy subjects and assessed the effects of principal component analysis (PCA)-based denoising on the repeatability parameters. In addition, spatial variations of ³¹ P metabolites within the liver were analyzed. 3D ³¹ P MRSI data of the liver were acquired with a nominal voxel size of 20 mm isotropic in 10 healthy volunteers twice on the same day. Data were reconstructed without denoising, and with PCA-based denoising before or after channel combination. From the test-retest data, repeatability parameters for metabolite level quantification were determined for 12 ³¹ P metabolite signals. On average, ³¹ P MR spectra from 100 ± 25 voxels in the liver were analyzed. Only voxels with contamination from skeletal muscle or the gall bladder were excluded and no voxels were discarded based on (low) signal-to-noise ratio (SNR). Repeatability for most quantified ³¹ P metabolite levels in the liver was good to excellent, with an intrasubject variability below 10%. PCA-based denoising increased the SNR ~ 3-fold, but did not improve the repeatability for mean liver ³¹ P metabolite quantification with the fitting constraints used. Significant spatial heterogeneity of various ³¹ P metabolite levels within the liver was observed, with marked differences for the phosphomonoester and phosphodiester metabolites between the left and right lobe. In conclusion, using a ³¹ P whole-body transmit coil in combination with a 16-channel ³¹ P receive array at 7 T allowed ³¹ P MRSI acquisitions with full liver coverage and good to excellent repeatability.

Apparent diffusion coefficients of 31P metabolites in the human calf muscle at 7 T

PURPOSE

In this study, we aimed to measure the apparent diffusion coefficients (ADCs) of major phosphorous metabolites in the human calf muscle at 7 T with a diffusion-weighted (DW)-STEAM sequence.

METHODS

A DW-STEAM sequence with bipolar gradients was implemented at 7 T, and DW MR spectra were acquired in three orthogonal directions in the human calf muscle of six healthy volunteers (TE/TM/TR = 15 ms/750 ms/5 s) at three b-values (0, 800, and 1200 s/mm²). Frequency and phase alignments were applied prior to spectral averaging. Averaged DW MR spectra were analyzed with LCModel, and ADCs of ³¹P metabolites were estimated.

RESULTS

Four metabolites (phosphocreatine (PCr), adenosine triphosphate (ATP), inorganic phosphate (Pi) and glycerol phosphorylcholine (GPC)) were quantified at all b-values with mean CRLBs below 10%. The ADC values of PCr, ATP, Pi, and GPC were (0.24 ± 0.02, 0.15 ± 0.04, 0.43 ± 0.14, 0.40 ± 0.09) × 10^-3 mm²/s, respectively.

CONCLUSION

The ADCs of four ³¹P metabolites were successfully measured in the human calf muscle at 7 T, among which those of ATP, Pi and GPC were reported for the first time in humans. This study paves the way to investigate ³¹P metabolite diffusion properties in health and disease on the clinical MR scanner.

Phosphorus transversal relaxation times and metabolite concentrations in the human brain at 9.4 T

A method to estimate phosphorus (³¹ P) transversal relaxation times (T₂ s) of coupled spin systems is demonstrated. Additionally, intracellular and extracellular pH and relaxation-corrected metabolite concentrations are reported. Echo time (TE) series of ³¹ P metabolite spectra were acquired using stimulated echo acquisition mode (STEAM) localization. Spectra were fitted using LCModel with accurately modeled Versatile Simulation, Pulses and Analysis (VeSPA) basis sets accounting for J-evolution of the coupled spin systems. T₂ s were estimated by fitting a single exponential two-parameter model across the TE series. Fitted inorganic phosphate frequencies were used to calculate pH, and estimated relaxation times were used to determine the relaxation-corrected brain metabolite concentrations on an assumption of 3 mM γ-ATP. The method was demonstrated in healthy human brain at a field strength of 9.4 T. T₂ times of ATP and nicotinamide adenine dinucleotide (NAD) were shortest between 8 and 20 ms, followed by T₂ s of inorganic phosphate between 25 and 50 ms, and phosphocreatine with a T₂ of 100 ms. Phosphomonoesters and phosphodiesters had the longest T₂ s of about 130 ms. The measured T₂ s are comparable with literature values and fit in a decreasing trend with increasing field strengths. Calculated pHs and metabolite concentrations are also comparable with literature values.

Evaluation of Acute Supplementation With the Ketone Ester (R)-3-Hydroxybutyl-(R)-3-Hydroxybutyrate (deltaG) in Healthy Volunteers by Cardiac and Skeletal Muscle 31P Magnetic Resonance Spectroscopy

In this acute intervention study, we investigated the potential benefit of ketone supplementation in humans by studying cardiac phosphocreatine to adenosine-triphosphate ratios (PCr/ATP) and skeletal muscle PCr recovery using phosphorus magnetic resonance spectroscopy (³¹P-MRS) before and after ingestion of a ketone ester drink. We recruited 28 healthy individuals: 12 aged 23–70 years for cardiac ³¹P-MRS, and 16 aged 60–75 years for skeletal muscle ³¹P-MRS. Baseline and post-intervention resting cardiac and dynamic skeletal muscle ³¹P-MRS scans were performed in one visit, where 25 g of the ketone monoester, deltaG^®, was administered after the baseline scan. Administration was timed so that post-intervention ³¹P-MRS would take place 30 min after deltaG^® ingestion. The deltaG^® ketone drink was well-tolerated by all participants. In participants who provided blood samples, post-intervention blood glucose, lactate and non-esterified fatty acid concentrations decreased significantly (−28.8%, p ≪ 0.001; −28.2%, p = 0.02; and −49.1%, p ≪ 0.001, respectively), while levels of the ketone body D-beta-hydroxybutyrate significantly increased from mean (standard deviation) 0.7 (0.3) to 4.0 (1.1) mmol/L after 30 min (p ≪ 0.001). There were no significant changes in cardiac PCr/ATP or skeletal muscle metabolic parameters between baseline and post-intervention. Acute ketone supplementation caused mild ketosis in blood, with drops in glucose, lactate, and free fatty acids; however, such changes were not associated with changes in ³¹P-MRS measures in the heart or in skeletal muscle. Future work may focus on the effect of longer-term ketone supplementation on tissue energetics in groups with compromised mitochondrial function.

Assessing the Feasibility of Dynamic 31P Spectroscopy for Metabolic Studies with a 1.0T Extremity Scanner

Objective: The feasibility of conducting in vivo non-localized ³¹P Magnetic Resonance Spectroscopy (MRS) with a 1.0T extremity scanner and the potential to increase accessibility of this important diagnostic tool for low cost applications is revisited. Methods: This work presents a custom transmit-only quadrature birdcage, four-element receive coil array, and spectrometer interfaced to a commercial ONI 1.0T magnet for enabling multi-channel, non-1H frequency capabilities. A custom, magnetic resonance compatible plantar flexion-extension exercise device was also developed to enable exercise protocols. The coils were assessed with bench measurements and ³¹P phantom studies before an in vivo demonstration. Results: In pulse and acquire spectroscopy of a phantom, the array was found to improve the signal-to-noise ratio (SNR) by a factor of 1.31 and reduce the linewidth by 13.9% when compared to a large loop coil of the same overall size. In vivo testing results show that two averages and a four second repetition time for a temporal resolution of eight seconds was sufficient to obtain phosphocreatine recovery values and baseline pH levels aligned with expected literature values. Conclusion: Initial in vivo human skeletal muscle ³¹P MRS allowed successful monitoring of metabolic changes during an 18-minute exercise protocol. Significance: Adding an array coil and multinuclear capability to a commercial low-cost 1.0T extremity scanner enabled the observation of characteristic ³¹P metabolic information, such as the phosphocreatine recovery rate and underlying baseline pH.

7 Tesla and Beyond: Advanced Methods and Clinical Applications in Magnetic Resonance Imaging

ABSTRACT

Ultrahigh magnetic fields offer significantly higher signal-to-noise ratio, and several magnetic resonance applications additionally benefit from a higher contrast-to-noise ratio, with static magnetic field strengths of B0 ≥ 7 T currently being referred to as ultrahigh fields (UHFs). The advantages of UHF can be used to resolve structures more precisely or to visualize physiological/pathophysiological effects that would be difficult or even impossible to detect at lower field strengths. However, with these advantages also come challenges, such as inhomogeneities applying standard radiofrequency excitation techniques, higher energy deposition in the human body, and enhanced B0 field inhomogeneities. The advantages but also the challenges of UHF as well as promising advanced methodological developments and clinical applications that particularly benefit from UHF are discussed in this review article.

Influence of Surface Ligand Molecular Structure on Phospholipid Membrane Disruption by Cationic Nanoparticles

Cationic nanoparticles are known to interact with biological membranes and often cause serious membrane damage. Therefore, it is important to understand the molecular mechanism for such interactions and the factors that impact the degree of membrane damage. Previously, we have demonstrated that spatial distribution of molecular charge at cationic nanoparticle surfaces plays an important role in determining the cellular uptake and membrane damage of these nanoparticles. In this work, using diamond nanoparticles (DNPs) functionalized with five different amine-based surface ligands and small phospholipid unilamellar vesicles (SUVs), we further investigate how chemical features and conformational flexibility of surface ligands impact nanoparticle/membrane interactions. ³¹P-NMR T₂ relaxation measurements quantify the mobility changes in lipid dynamics upon exposing the SUVs to functional DNPs, and coarse-grained molecular dynamics simulations further elucidate molecular details for the different modes of DNP-SUV interactions depending on the surface ligands. Collectively, our results show that the length of the hydrophobic segment and conformational flexibility of surface ligands are two key factors that dictate the degree of membrane damage by the DNP, while the amount of surface charge alone is not predictive of the strength of interaction.

3D 31 P MR spectroscopic imaging of the human brain at 3 T with a 31 P receive array: An assessment of 1 H decoupling, T1 relaxation times, 1 H-31 P nuclear Overhauser effects and NAD. NMR IN BIOMEDICINE 2021;34:e4169. [PMID: 31518036 PMCID: PMC8244063 DOI: 10.1002/nbm.4169]

³¹ P MR spectroscopic imaging (MRSI) is a versatile technique to study phospholipid precursors and energy metabolism in the healthy and diseased human brain. However, mainly due to its low sensitivity, ³¹ P MRSI is currently limited to research purposes. To obtain 3D ³¹ P MRSI spectra with improved signal-to-noise ratio on clinical 3 T MR systems, we used a coil combination consisting of a dual-tuned birdcage transmit coil and a ³¹ P eight-channel phased-array receive insert. To further increase resolution and sensitivity we applied WALTZ4 ¹ H decoupling and continuous wave nuclear Overhauser effect (NOE) enhancement and acquired high-quality MRSI spectra with nominal voxel volumes of ~ 17.6 cm³ (effective voxel volume ~ 51 cm³ ) in a clinically relevant measurement time of ~ 13 minutes, without exceeding SAR limits. Steady-state NOE enhancements ranged from 15 ± 9% (γ-ATP) and 33 ± 3% (phosphocreatine) to 48 ± 11% (phosphoethanolamine). Because of these improvements, we resolved and detected all ³¹ P signals of metabolites that have also been reported for ultrahigh field strengths, including resonances for NAD⁺ , NADH and extracellular inorganic phosphate. T₁ times of extracellular inorganic phosphate were longer than for intracellular inorganic phosphate (3.8 ± 1.4s vs 1.8 ± 0.65 seconds). A comparison of measured T₁ relaxation times and NOE enhancements at 3 T with published values between 1.5 and 9.4 T indicates that T₁ relaxation of ³¹ P metabolite spins in the human brain is dominated by dipolar relaxation for this field strength range. Even although intrinsic sensitivity is higher at ultrahigh fields, we demonstrate that at a clinical field strength of 3 T, similar ³¹ P MRSI information content can be obtained using a sophisticated coil design combined with ¹ H decoupling and NOE enhancement.

Phosphorus Kβ X-ray emission spectroscopy detects non-covalent interactions of phosphate biomolecules in situ

Phosphorus is ubiquitous in biochemistry, being found in the phosphate groups of nucleic acids and the energy-transferring system of adenine nucleotides (e.g. ATP). Kβ X-ray emission spectroscopy (XES) of phosphorus has been largely unexplored, with no previous applications to biomolecules. Here, the potential of P Kβ XES to study phosphate-containing biomolecules, including ATP and NADPH, is evaluated, as is the application of the technique to aqueous solution samples. P Kβ spectra offer a detailed picture of phosphate valence electronic structure, reporting on subtle non-covalent effects, such as hydrogen bonding and ionic interactions, that are key to enzymatic catalysis. Spectral features are interpreted using density functional theory (DFT) calculations, and potential applications to the study of biological energy conversion are highlighted.

Phosphorus X-ray emission spectroscopy probes non-covalent interactions and electronic structure of phosphate biomolecules in both solid and solution samples.

Interval-induced metabolic perturbation determines tissue fluid shifts into skeletal muscle

Intense interval exercise has proven to be as effective as traditional endurance exercise in improving maximal oxygen uptake. Shared by these two exercise regimes is an acute reduction in plasma volume, which is a suggested stimulus behind exercise-induced increases in blood volume and maximal oxygen uptake. This study aimed to link exercise-induced metabolic perturbation with volume shifts into skeletal muscle tissue. Ten healthy subjects (mean age 33 ± 8 years, 5 males and 5 females) performed three 30 s all-out sprints on a cycle ergometer. Upon cessation of exercise magnetic resonance imaging, ³¹ Phosphorus magnetic resonance spectroscopy and blood samples were used to measure changes in muscle volume, intramuscular energy metabolites and plasma volume. Compared to pre-exercise, muscle volume increased from 1147.1 ± 35.6 ml to 1283.3 ± 11.0 ml 8 min post-exercise. At 30 min post-exercise, muscle volume was still higher than pre-exercise (1147.1 ± 35.6 vs. 1222.2 ± 6.8 ml). Plasma volume decreased by 16 ± 3% immediately post-exercise and recovered back to - 5 ± 6% after 30 min. Principal component analysis of exercise performance, muscle and plasma volume changes as well as changes in intramuscular energy metabolites showed generally strong correlations between metabolic and physiological variables. The strongest predictor for the volume shifts of muscle and plasma was the magnitude of glucose-6-phosphate accumulation post-exercise. Interval training leads to large metabolic and hemodynamic perturbations with accumulation of glucose-6-phosphate as a possible key event in the fluid flux between the vascular compartment and muscle tissue.

Magnetic resonance spectroscopy reveals mitochondrial dysfunction in amyotrophic lateral sclerosis

Mitochondrial dysfunction is postulated to be central to amyotrophic lateral sclerosis (ALS) pathophysiology. Evidence comes primarily from disease models and conclusive data to support bioenergetic dysfunction in vivo in patients is currently lacking. This study is the first to assess mitochondrial dysfunction in brain and muscle in individuals living with ALS using 31P-magnetic resonance spectroscopy (MRS), the modality of choice to assess energy metabolism in vivo. We recruited 20 patients and 10 healthy age and gender-matched control subjects in this cross-sectional clinico-radiological study. 31P-MRS was acquired from cerebral motor regions and from tibialis anterior during rest and exercise. Bioenergetic parameter estimates were derived including: ATP, phosphocreatine, inorganic phosphate, adenosine diphosphate, Gibbs free energy of ATP hydrolysis (ΔGATP), phosphomonoesters, phosphodiesters, pH, free magnesium concentration, and muscle dynamic recovery constants. Linear regression was used to test for associations between brain data and clinical parameters (revised amyotrophic functional rating scale, slow vital capacity, and upper motor neuron score) and between muscle data and clinico-neurophysiological measures (motor unit number and size indices, force of contraction, and speed of walking). Evidence for primary dysfunction of mitochondrial oxidative phosphorylation was detected in the brainstem where ΔGATP and phosphocreatine were reduced. Alterations were also detected in skeletal muscle in patients where resting inorganic phosphate, pH, and phosphomonoesters were increased, whereas resting ΔGATP, magnesium, and dynamic phosphocreatine to inorganic phosphate recovery were decreased. Phosphocreatine in brainstem correlated with respiratory dysfunction and disability; in muscle, energy metabolites correlated with motor unit number index, muscle power, and speed of walking. This study provides in vivo evidence for bioenergetic dysfunction in ALS in brain and skeletal muscle, which appears clinically and electrophysiologically relevant. 31P-MRS represents a promising technique to assess the pathophysiology of mitochondrial function in vivo in ALS and a potential tool for future clinical trials targeting bioenergetic dysfunction.

PCA denoising and Wiener deconvolution of 31 P 3D CSI data to enhance effective SNR and improve point spread function

Purpose

This study evaluates the performance of 2 processing methods, that is, principal component analysis‐based denoising and Wiener deconvolution, to enhance the quality of phosphorus 3D chemical shift imaging data.

Methods

Principal component analysis‐based denoising increases the SNR while maintaining spectral information. Wiener deconvolution reduces the FWHM of the voxel point spread function, which is increased by Hamming filtering or Hamming‐weighted acquisition. The proposed methods are evaluated using simulated and in vivo 3D phosphorus chemical shift imaging data by 1) visual inspection of the spatial signal distribution; 2) SNR calculation of the PCr peak; and 3) fitting of metabolite basis functions.

Results

With the optimal order of processing steps, we show that the effective SNR of in vivo phosphorus 3D chemical shift imaging data can be increased. In simulations, we show we can preserve phosphorus‐containing metabolite peaks that had an SNR < 1 before denoising. Furthermore, using Wiener deconvolution, we were able to reduce the FWHM of the voxel point spread function with only partially reintroducing Gibb‐ringing artifacts while maintaining the SNR. After data processing, fitting of the phosphorus‐containing metabolite signals improved.

Conclusion

In this study, we have shown that principal component analysis‐based denoising in combination with regularized Wiener deconvolution allows increasing the effective spectral SNR of in vivo phosphorus 3D chemical shift imaging data, with reduction of the FWHM of the voxel point spread function. Processing increased the effective SNR by at least threefold compared to Hamming weighted acquired data and minimized voxel bleeding. With these methods, fitting of metabolite amplitudes became more robust with decreased fitting residuals.

Musculoskeletal MR Imaging Applications at Ultra-High (7T) Field Strength

Regulatory approval of ultrahigh field (UHF) MR imaging scanners for clinical use has opened new opportunities for musculoskeletal imaging applications. UHF MR imaging has unique advantages in terms of signal-to-noise ratio, contrast-to-noise ratio, spectral resolution, and multinuclear applications, thus providing unique information not available at lower field strengths. But UHF also comes with a set of technical challenges that are yet to be resolved and may not be suitable for all imaging applications. This review focuses on the latest research in musculoskeletal MR imaging applications at UHF including morphologic imaging, T₂, T₂∗, and T_1ρ mapping, chemical exchange saturation transfer, sodium imaging, and phosphorus spectroscopy imaging applications.

Multi-parametric MR in Becker muscular dystrophy patients

Quantitative MRI and MRS of muscle are increasingly being used to measure individual pathophysiological processes in Becker muscular dystrophy (BMD). In particular, muscle fat fraction was shown to be highly associated with functional tests in BMD. However, the muscle strength per unit of contractile cross-sectional area is lower in patients with BMD compared with healthy controls. This suggests that the quality of the non-fat-replaced (NFR) muscle tissue is lower than in healthy controls. Consequently, a measure that reflects changes in muscle tissue itself is needed. Here, we explore the potential of water T₂ relaxation times, diffusion parameters and phosphorus metabolic indices as early disease markers in patients with BMD. For this purpose, we examined these measures in fat-replaced (FR) and NFR lower leg muscles in patients with BMD and compared these values with those in healthy controls. Quantitative proton MRI (three-point Dixon, multi-spin-echo and diffusion-weighted spin-echo echo planar imaging) and 2D chemical shift imaging ³¹ P MRS data were acquired in 24 patients with BMD (age 18.8-66.2 years) and 13 healthy controls (age 21.3-63.6 years). Muscle fat fractions, phosphorus metabolic indices, and averages and standard deviations (SDs) of the water T₂ relaxation times and diffusion tensor imaging (DTI) parameters were assessed in six individual leg muscles. Phosphodiester levels were increased in the NFR and FR tibialis anterior, FR peroneus and FR gastrocnemius lateralis muscles. No clear pattern was visible for the other metabolic indices. Increased T₂ SD was found in the majority of FR muscles compared with NFR and healthy control muscles. No differences in average water T₂ relaxation times or DTI indices were found between groups. Overall, our results indicate that primarily muscles that are further along in the disease process showed increases in T₂ heterogeneity and changes in some metabolic indices. No clear differences were found for the DTI indices between groups.

Evidence of Mitochondrial Dysfunction in Fibromyalgia: Deviating Muscle Energy Metabolism Detected Using Microdialysis and Magnetic Resonance

In fibromyalgia (FM) muscle metabolism, studies are sparse and conflicting associations have been found between muscle metabolism and pain aspects. This study compared alterations in metabolic substances and blood flow in erector spinae and trapezius of FM patients and healthy controls. FM patients (n = 33) and healthy controls (n = 31) underwent a clinical examination that included pressure pain thresholds and physical tests, completion of a health questionnaire, participation in microdialysis investigations of the etrapezius and erector spinae muscles, and also underwent phosphorus-31 magnetic resonance spectroscopy of the erector spinae muscle. At the baseline, FM had significantly higher levels of pyruvate in both muscles. Significantly lower concentrations of phosphocreatine (PCr) and nucleotide triphosphate (mainly adenosine triphosphate) in erector spinae were found in FM. Blood flow in erector spinae was significantly lower in FM. Significant associations between metabolic variables and pain aspects (pain intensity and pressure pain threshold PPT) were found in FM. Our results suggest that FM has mitochondrial dysfunction, although it is unclear whether inactivity, obesity, aging, and pain are causes of, the results of, or coincidental to the mitochondrial dysfunction. The significant regressions of pain intensity and PPT in FM agree with other studies reporting associations between peripheral biological factors and pain aspects.

Effects of contrast agents on relaxation properties of 31P metabolites

PURPOSE

Phosphorous MR spectroscopy (31P-MRS) forms a powerful, non-invasive research tool to quantify the energetics of the heart in diverse patient populations. 31P-MRS is frequently applied alongside other radiological examinations, many of which use various contrast agents that shorten relaxation times of water in conventional proton MR, for a better characterisation of cardiac function, or following prior computed tomography (CT). It is, however, unknown whether these agents confound 31P-MRS signals, for example, 2,3-diphosphoglycerate (2,3-DPG).

METHODS

In this work, we quantitatively assess the impact of non-ionic, low osmolar iodinated CT contrast agent (iopamidol/Niopam), gadolinium chelates (linear gadopentetic acid dimeglumine/Magnevist and macrocyclic gadoterate meglumine/Dotarem) and superparamagnetic iron oxide nanoparticles (ferumoxytol/Feraheme) on the nuclear T1 and T2 of 31P metabolites (ie, 2,3-DPG), and 1H in water in live human blood and saline phantoms at 11.7 T.

RESULTS

Addition of all contrast agents led to significant shortening of all relaxation times in both 1H and 31P saline phantoms. On the contrary, the T1 relaxation time of 2,3-DPG in blood was significantly shortened only by Magnevist (P = .03). Similarly, the only contrast agent that influenced the T2 relaxation times of 2,3-DPG in blood samples was ferumoxytol (P = .02).

CONCLUSION

Our results show that, unlike conventional proton MR, phosphorus MRS is unconfounded in patients who have had prior CT with contrast, not all gadolinium-based contrast agents influence 31P-MRS data in vivo, and that ferumoxytol is a promising contrast agent for the reduction in 31P-MRS blood-pool signal.

Quantifying the effect of dobutamine stress on myocardial Pi and pH in healthy volunteers: A 31 P MRS study at 7T

Purpose

Phosphorus spectroscopy (³¹P‐MRS) is a proven method to probe cardiac energetics. Studies typically report the phosphocreatine (PCr) to adenosine triphosphate (ATP) ratio. We focus on another ³¹P signal: inorganic phosphate (Pi), whose chemical shift allows computation of myocardial pH, with Pi/PCr providing additional insight into cardiac energetics. Pi is often obscured by signals from blood 2,3‐diphosphoglycerate (2,3‐DPG). We introduce a method to quantify Pi in 14 min without hindrance from 2,3‐DPG.

Methods

Using a ³¹P stimulated echo acquisition mode (STEAM) sequence at 7 Tesla that inherently suppresses signal from 2,3‐DPG, the Pi peak was cleanly resolved. Resting state UTE‐chemical shift imaging (PCr/ATP) and STEAM ³¹P‐MRS (Pi/PCr, pH) were undertaken in 23 healthy controls; pH and Pi/PCr were subsequently recorded during dobutamine infusion.

Results

We achieved a clean Pi signal both at rest and stress with good 2,3‐DPG suppression. Repeatability coefficient (8 subjects) for Pi/PCr was 0.036 and 0.12 for pH. We report myocardial Pi/PCr and pH at rest and during catecholamine stress in healthy controls. Pi/PCr was maintained during stress (0.098 ± 0.031 [rest] vs. 0.098 ± 0.031 [stress] P = .95); similarly, pH did not change (7.09 ± 0.07 [rest] vs. 7.08 ± 0.11 [stress] P = .81). Feasibility for patient studies was subsequently successfully demonstrated in a patient with cardiomyopathy.

Conclusion

We introduced a method that can resolve Pi using 7 Tesla STEAM ³¹P‐MRS. We demonstrate the stability of Pi/PCr and myocardial pH in volunteers at rest and during catecholamine stress. This protocol is feasible in patients and potentially of use for studying pathological myocardial energetics.

Multinuclear MRS at 7T Uncovers Exercise Driven Differences in Skeletal Muscle Energy Metabolism Between Young and Seniors

Purpose: Aging is associated with changes in muscle energy metabolism. Proton (¹H) and phosphorous (³¹P) magnetic resonance spectroscopy (MRS) has been successfully applied for non-invasive investigation of skeletal muscle metabolism. The aim of this study was to detect differences in adenosine triphosphate (ATP) production in the aging muscle by ³¹P-MRS and to identify potential changes associated with buffer capacity of muscle carnosine by ¹H-MRS.

Methods: Fifteen young and nineteen elderly volunteers were examined. ¹H and ³¹P-MRS spectra were acquired at high field (7T). The investigation included carnosine quantification using ¹H-MRS and resting and dynamic ³¹P-MRS, both including saturation transfer measurements of phosphocreatine (PCr), and inorganic phosphate (Pi)-to-ATP metabolic fluxes.

Results: Elderly volunteers had higher time constant of PCr recovery (τ_PCr) in comparison to the young volunteers. Exercise was connected with significant decrease in PCr-to-ATP flux in both groups. Moreover, PCr-to-ATP flux was significantly higher in young compared to elderly both at rest and during exercise. Similarly, an increment of Pi-to-ATP flux with exercise was found in both groups but the intergroup difference was only observed during exercise. Elderly had lower muscle carnosine concentration and lower postexercise pH. A strong increase in phosphomonoester (PME) concentration was observed with exercise in elderly, and a faster Pi:PCr kinetics was found in young volunteers compared to elderly during the recovery period.

Conclusion: Observations of a massive increment of PME concentration together with high Pi-to-ATP flux during exercise in seniors refer to decreased ability of the muscle to meet the metabolic requirements of exercise and thus a limited ability of seniors to effectively support the exercise load.

31 P magnetic resonance spectroscopy in skeletal muscle: Experts' consensus recommendations

Skeletal muscle phosphorus-31 31 P MRS is the oldest MRS methodology to be applied to in vivo metabolic research. The technical requirements of 31 P MRS in skeletal muscle depend on the research question, and to assess those questions requires understanding both the relevant muscle physiology, and how 31 P MRS methods can probe it. Here we consider basic signal-acquisition parameters related to radio frequency excitation, TR, TE, spectral resolution, shim and localisation. We make specific recommendations for studies of resting and exercising muscle, including magnetisation transfer, and for data processing. We summarise the metabolic information that can be quantitatively assessed with 31 P MRS, either measured directly or derived by calculations that depend on particular metabolic models, and we give advice on potential problems of interpretation. We give expected values and tolerable ranges for some measured quantities, and minimum requirements for reporting acquisition parameters and experimental results in publications. Reliable examination depends on a reproducible setup, standardised preconditioning of the subject, and careful control of potential difficulties, and we summarise some important considerations and potential confounders. Our recommendations include the quantification and standardisation of contraction intensity, and how best to account for heterogeneous muscle recruitment. We highlight some pitfalls in the assessment of mitochondrial function by analysis of phosphocreatine (PCr) recovery kinetics. Finally, we outline how complementary techniques (near-infrared spectroscopy, arterial spin labelling, BOLD and various other MRI and 1 H MRS measurements) can help in the physiological/metabolic interpretation of 31 P MRS studies by providing information about blood flow and oxygen delivery/utilisation. Our recommendations will assist in achieving the fullest possible reliable picture of muscle physiology and pathophysiology.

Feasibility of Using a 1T Extremity Scanner with a Four-Element Array to Detect 31P in the Human Calf

The feasibility of conducting in vivo non-localized skeletal muscle ³¹P Magnetic Resonance Spectroscopy (MRS) with a low-cost extremity 1 Tesla magnet is demonstrated. We designed and built a transmit-only quadrature birdcage, four-element receive coil array, and employed a home-built spectrometer interfaced with a commercial ONI 1.0T magnet. In phantom comparison tests with a large loop coil of comparable size, the array was found to improve the SNR by a factor of 1.8 and the linewidth from 0.72 ppm to 0.45 ppm. Phantom and in vivo testing results show only 6 averages with a 4 second repetition time are required to obtain quantifiable ³¹P spectra. Initial in vivo human skeletal muscle ³¹P spectra successfully allowed for peak characterization. A low-cost approach to MRS could enable more widespread use of this tool in clinical diagnosis and in vivo metabolic research.

Low SAR 31 P (multi-echo) spectroscopic imaging using an integrated whole-body transmit coil at 7T

Phosphorus (³¹ P) MRSI provides opportunities to monitor potential biomarkers. However, current applications of ³¹ P MRS are generally restricted to relatively small volumes as small coils are used. Conventional surface coils require high energy adiabatic RF pulses to achieve flip angle homogeneity, leading to high specific absorption rates (SARs), and occupy space within the MRI bore. A birdcage coil behind the bore cover can potentially reduce the SAR constraints massively by use of conventional amplitude modulated pulses without sacrificing patient space. Here, we demonstrate that the integrated ³¹ P birdcage coil setup with a high power RF amplifier at 7 T allows for low flip angle excitations with short repetition time (T_R ) for fast 3D chemical shift imaging (CSI) and 3D T₁ -weighted CSI as well as high flip angle multi-refocusing pulses, enabling multi-echo CSI that can measure metabolite T₂ , over a large field of view in the body. B₁⁺ calibration showed a variation of only 30% in maximum B₁ in four volunteers. High signal-to-noise ratio (SNR) MRSI was obtained in the gluteal muscle using two fast in vivo 3D spectroscopic imaging protocols, with low and high flip angles, and with multi-echo MRSI without exceeding SAR levels. In addition, full liver MRSI was achieved within SAR constraints. The integrated ³¹ P body coil allowed for fast spectroscopic imaging and successful implementation of the multi-echo method in the body at 7 T. Moreover, no additional enclosing hardware was needed for ³¹ P excitation, paving the way to include larger subjects and more space for receiver arrays. The increase in possible number of RF excitations per scan time, due to the improved B₁⁺ homogeneity and low SAR, allows SNR to be exchanged for spatial resolution in CSI and/or T₁ weighting by simply manipulating T_R and/or flip angle to detect and quantify ratios from different molecular species.

Stimulated echo double diffusion encoded imaging of closed pores: Influence and removal of unbalanced terms

Nuclear magnetic resonance (NMR) diffusion pore imaging has been proposed to study the shape of arbitrary closed pores filled with an NMR-detectable medium by use of nonclassical diffusion encoding schemes. Potential applications can be found in biomedical imaging and porous media research. When studying non-point-symmetric pores, NMR signals with nonvanishing imaginary parts arise containing the pore shape information, which is lost for classical diffusion encoding schemes. Key limitations are the required high magnetic field gradient amplitudes and T2 relaxation while approaching the diffusion long-time limit. To benefit from the slower T1 decay, we demonstrate the feasibility of diffusion pore imaging with stimulated echoes using Monte Carlo simulations and experiments with hyperpolarized xenon-129 gas in well-defined geometries and show that the necessary complex-valued signals can be acquired. Analytical derivation of the stimulated echo double diffusion encoded signal was performed to investigate the effect of the additionally arising undesired terms on the complex phase information. These terms correspond to signals arising for spin-echo sequences with unbalanced gradients. For most possible applications, the unbalanced terms can be neglected. If non-negligible, selection of the appropriate signal component using a phase cycling scheme was demonstrated experimentally. Using stimulated echoes may be a step towards application of diffusion pore imaging to larger pores with gradient amplitudes available today in preclinical systems.

Shortening of apparent transverse relaxation time of inorganic phosphate as a breast cancer biomarker

Phosphorus MRS offers a non-invasive tool for monitoring cell energy and phospholipid metabolism and can be of additional value in diagnosing cancer and monitoring cancer therapy. In this study, we determined the transverse relaxation times of a number of phosphorous metabolites in a group of breast cancer patients by adiabatic multi-echo spectroscopic imaging at 7 T. The transverse relaxation times of phosphoethanolamine, phosphocholine, inorganic phosphate (P_i ), glycerophosphocholine and glycerophosphatidylcholine were 184 ± 8 ms, 203 ± 17 ms, 87 ± 8 ms, 240 ± 56 ms and 20 ± 10 ms, respectively. The transverse relaxation time of P_i in breast cancer tissue was less than half that of healthy fibroglandular tissue. This effect is most likely caused by an up-regulation of glycolysis in breast cancer tissue that leads to interaction of P_i with the GAPDH enzyme, which forms part of the reversible pathway of exchange of P_i with gamma-adenosine tri-phosphate, thus shortening its apparent transverse relaxation time. As healthy breast tissue shows very little glycolytic activity, the apparent T₂ shortening of P_i due to malignant transformation could possibly be used as a biomarker for cancer.

Measuring inorganic phosphate and intracellular pH in the healthy and hypertrophic cardiomyopathy hearts by in vivo 7T 31P-cardiovascular magnetic resonance spectroscopy

BACKGROUND

Cardiovascular phosphorus MR spectroscopy (31P-CMRS) is a powerful tool for probing energetics in the human heart, through quantification of phosphocreatine (PCr) to adenosine triphosphate (ATP) ratio. In principle, 31P-CMRS can also measure cardiac intracellular pH (pHi) and the free energy of ATP hydrolysis (ΔGATP). However, these require determination of the inorganic phosphate (Pi) signal frequency and amplitude that are currently not robustly accessible because blood signals often obscure the Pi resonance. Typical cardiac 31P-CMRS protocols use low (e.g. 30°) flip-angles and short repetition time (TR) to maximise signal-to-noise ratio (SNR) within hardware limits. Unfortunately, this causes saturation of Pi with negligible saturation of the flowing blood pool. We aimed to show that an adiabatic 90° excitation, long-TR, 7T 31P-CMRS protocol will reverse this balance, allowing robust cardiac pHi measurements in healthy subjects and patients with hypertrophic cardiomyopathy (HCM).

METHODS

The cardiac Pi T1 was first measured by the dual TR technique in seven healthy subjects. Next, ten healthy subjects and three HCM patients were scanned with 7T 31P-MRS using long (6 s) TR protocol and adiabatic excitation. Spectra were fitted for cardiac metabolites including Pi.

RESULTS

The measured Pi T1 was 5.0 ± 0.3 s in myocardium and 6.4 ± 0.6 s in skeletal muscle. Myocardial pH was 7.12 ± 0.04 and Pi/PCr ratio was 0.11 ± 0.02. The coefficients of repeatability were 0.052 for pH and 0.027 for Pi/PCr quantification. The pH in HCM patients did not differ (p = 0.508) from volunteers. However, Pi/PCr was higher (0.24 ± 0.09 vs. 0.11 ± 0.02; p = 0.001); Pi/ATP was higher (0.44 ± 0.14 vs. 0.24 ± 0.05; p = 0.002); and PCr/ATP was lower (1.78 ± 0.07 vs. 2.10 ± 0.20; p = 0.020), in HCM patients, which is in agreement with previous reports.

CONCLUSION

A 7T 31P-CMRS protocol with adiabatic 90° excitation and long (6 s) TR gives sufficient SNR for Pi and low enough blood signal to permit robust quantification of cardiac Pi and hence pHi. Pi was detectable in every subject scanned for this study, both in healthy subjects and HCM patients. Cardiac pHi was unchanged in HCM patients, but both Pi/PCr and Pi/ATP increased that indicate an energetic impairment in HCM. This work provides a robust technique to quantify cardiac Pi and pHi.

Modular 31 P wideband inversion transfer for integrative analysis of adenosine triphosphate metabolism, T1 relaxation and molecular dynamics in skeletal muscle at 7T

PURPOSE

For efficient and integrative analysis of de novo adenosine triphosphate (ATP) synthesis, creatine-kinase-mediated ATP synthesis, T₁ relaxation time, and ATP molecular motion dynamics in human skeletal muscle at rest.

METHODS

Four inversion-transfer modules differing in center inversion frequency were combined to generate amplified magnetization transfer (MT) effects in targeted MT pathways, including Pi ↔ γ-ATP, PCr ↔ γ-ATP, and ³¹ P_γ(α)ATP ↔ ³¹ P_βATP . MT effects from both forward and reverse exchange kinetic pathways were acquired to reduce potential bias and confounding factors in integrated data analysis.

RESULTS

Kinetic data collected using 4 wideband inversion modules (8 minutes each) yielded the forward exchange rate constants, k_{PCr →γ ATP} = 0.31 ± 0.05 s^-1 and k_{Pi →γ ATP} = 0.064 ± 0.012 s^-1 , and the reverse exchange rate constants, k_γATP→Pi = 0.034 ± 0.006 s^-1 and k_γATP→PCr = 1.37 ± 0.22 s^-1 , respectively. The cross-relaxation rate constant, σ_{γ(α) ↔ βATP} was -0.20 ± 0.03 s^-1 , corresponding to ATP rotational correlation time τ_c of 0.8 ± 0.1 × 10^-7 seconds. The intrinsic T₁ relaxation times were Pi (9.2 ± 1.4 seconds), PCr (6.2 ± 0.4 seconds), γ-ATP (1.8 ± 0.1 seconds), α-ATP (1.4 ± 0.1 seconds), and β-ATP (1.1 ± 0.1 seconds). Muscle ATP T₁ values were found to be significantly longer than those previously measured in the brain using a similar method.

CONCLUSION

A combination of multiple inversion transfer modules provides a comprehensive and integrated analysis of ATP metabolism and molecular motion dynamics. This relatively fast technique could be potentially useful for studying metabolic disorders in skeletal muscle.

MR Imaging of the Musculoskeletal System Using Ultrahigh Field (7T) MR Imaging

MR imaging is an indispensable instrument for the diagnosis of musculoskeletal diseases. In vivo MR imaging at 7T offers many advantages, including increased signal-to-noise ratio, higher spatial resolution, improved spectral resolution for spectroscopy, improved sensitivity for X-nucleus imaging, and decreased image acquisition times. There are also however technical challenges of imaging at a higher field strength compared with 1.5 and 3T MR imaging systems. We discuss the many potential opportunities as well as the challenges presented by 7T MR imaging systems and highlight recent developments in in vivo research imaging of musculoskeletal applications in general and cartilage, skeletal muscle, and bone in particular.

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.

In vivo bone 31 P relaxation times and their implications on mineral quantification

PURPOSE

The intersubject variations in bone phosphorus-31 (³¹ P) T₁ and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> , as well as the implications on in vivo ³¹ P MRI-based bone mineral quantification, were investigated at 3T field strength.

METHODS

A technique that isolates the bone signal from the composite in vivo ³¹ P spectrum was first evaluated via simulation and experiments ex vivo and subsequently applied to measure the T₁ of bone ³¹ P collectively with a spectroscopic saturation recovery sequence in a group of healthy subjects aged 26 to 76 years. <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> was derived from the bone signal linewidth. The density of bone ³¹ P was derived for all subjects from ³¹ P zero TE images acquired in the same scan session using the measured relaxation times. Test-retest experiments were also performed to evaluate repeatability of this in vivo MRI-based bone mineral quantification protocol.

RESULTS

The T₁ obtained in vivo using the proposed spectral separation method combined with saturation recovery sequence is 38.4 ± 1.5 s for the subjects studied. Average ³¹ P density found was 6.40 ± 0.58 mol/L (corresponding to 1072 ± 98 mg/cm³ mineral density), in good agreement with an earlier study in specimens from donors of similar age range. Neither the relaxation times (P = 0.18 for T₁ , P = 0.99 for <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> ) nor ³¹ P density (P = 0.55) were found to correlate with subject age. Average coefficients of variation for the repeat study were 1.5%, 2.6%, and 4.4% for bone ³¹ P T₁ , <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> , and density, respectively.

CONCLUSION

Neither ³¹ P T₁ nor <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> varies significantly in healthy adults across a 50-year age range, therefore obviating the need for subject-specific measurements.

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.

Using 31P-MRI of hydroxyapatite for bone attenuation correction in PET-MRI: proof of concept in the rodent brain

BACKGROUND

The correction of γ-photon attenuation in PET-MRI remains a critical issue, especially for bone attenuation. This problem is of great importance for brain studies due to the density of the skull. Current techniques for skull attenuation correction (AC) provide indirect estimates of cortical bone density, leading to inaccurate estimates of brain activity. The purpose of this study was to develop an alternate method for bone attenuation correction based on NMR. The proposed approach relies on the detection of hydroxyapatite crystals by zero echo time (ZTE) MRI of 31P, providing individual and quantitative assessment of bone density. This work presents a proof of concept of this approach. The first step of the method is a calibration experiment to determine the conversion relationship between the 31P signal and the linear attenuation coefficient μ. Then 31P-ZTE was performed in vivo in rodent to estimate the μ-map of the skull. 18F-FDG PET data were acquired in the same animal and reconstructed with three different AC methods: 31P-based AC, AC neglecting the bone and the gold standard, CT-based AC, used to comparison for the other two methods.

RESULTS

The calibration experiment provided a conversion factor of 31P signal into μ. In vivo 31P-ZTE made it possible to acquire 3D images of the rat skull. Brain PET images showed underestimation of 18F activity in peripheral regions close to the skull when AC neglected the bone (as compared with CT-based AC). The use of 31P-derived μ-map for AC leads to increased peripheral activity, and therefore a global overestimation of brain 18F activity.

CONCLUSIONS

In vivo 31P-ZTE MRI of hydroxyapatite provides μ-map of the skull, which can be used for attenuation correction of 18F-FDG PET images. This study is limited by several intrinsic biases associated with the size of the rat brain, which are unlikely to affect human data on a clinical PET-MRI system.

31 P magnetic resonance fingerprinting for rapid quantification of creatine kinase reaction rate in vivo

The purpose of this work was to develop a 31 P spectroscopic magnetic resonance fingerprinting (MRF) method for fast quantification of the chemical exchange rate between phosphocreatine (PCr) and adenosine triphosphate (ATP) via creatine kinase (CK). A 31 P MRF sequence (CK-MRF) was developed to quantify the forward rate constant of ATP synthesis via CK ( kfCK), the T1 relaxation time of PCr ( T1PCr), and the PCr-to-ATP concentration ratio ( MRPCr). The CK-MRF sequence used a balanced steady-state free precession (bSSFP)-type excitation with ramped flip angles and a unique saturation scheme sensitive to the exchange between PCr and γATP. Parameter estimation was accomplished by matching the acquired signals to a dictionary generated using the Bloch-McConnell equation. Simulation studies were performed to examine the susceptibility of the CK-MRF method to several potential error sources. The accuracy of nonlocalized CK-MRF measurements before and after an ischemia-reperfusion (IR) protocol was compared with the magnetization transfer (MT-MRS) method in rat hindlimb at 9.4 T (n = 14). The reproducibility of CK-MRF was also assessed by comparing CK-MRF measurements with both MT-MRS (n = 17) and four angle saturation transfer (FAST) (n = 7). Simulation results showed that CK-MRF quantification of kfCK was robust, with less than 5% error in the presence of model inaccuracies including dictionary resolution, metabolite T2 values, inorganic phosphate metabolism, and B1 miscalibration. Estimation of kfCK by CK-MRF (0.38 ± 0.02 s-1 at baseline and 0.42 ± 0.03 s-1 post-IR) showed strong agreement with MT-MRS (0.39 ± 0.03 s-1 at baseline and 0.44 ± 0.04 s-1 post-IR). kfCK estimation was also similar between CK-MRF and FAST (0.38 ± 0.02 s-1 for CK-MRF and 0.38 ± 0.11 s-1 for FAST). The coefficient of variation from 20 s CK-MRF quantification of kfCK was 42% of that by 150 s MT-MRS acquisition and was 12% of that by 20 s FAST acquisition. This study demonstrates the potential of a 31 P spectroscopic MRF framework for rapid, accurate and reproducible quantification of chemical exchange rate of CK in vivo.

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.

Spatially localized phosphorous metabolism of skeletal muscle in Duchenne muscular dystrophy patients: 24-month follow-up

Objectives

To assess the changes in phosphodiester (PDE)-levels, detected by ³¹P magnetic resonance spectroscopy (MRS), over 24-months to determine the potential of PDE as marker for muscle tissue changes in Duchenne Muscular Dystrophy (DMD) patients.

Methods

Spatially resolved phosphorous datasets were acquired in the right lower leg of 18 DMD patients (range: 5–15.4 years) and 12 age-matched healthy controls (range: 5–14 years) at three time-points (baseline, 12-months, and 24-months) using a 7T MR-System (Philips Achieva). 3-point Dixon images were acquired at 3T (Philips Ingenia) to determine muscle fat fraction. Analyses were done for six muscles that represent different stages of muscle wasting. Differences between groups and time-points were assessed with non-parametric tests with correction for multiple comparisons. Coefficient of variance (CV) were determined for PDE in four healthy adult volunteers in high and low signal-to-noise ratio (SNR) datasets.

Results

PDE-levels were significantly higher (two-fold) in DMD patients compared to controls in all analyzed muscles at almost every time point and did not change over the study period. Fat fraction was significantly elevated in all muscles at all time points compared to healthy controls, and increased significantly over time, except in the tibialis posterior muscle. The mean within subject CV for PDE-levels was 4.3% in datasets with high SNR (>10:1) and 5.7% in datasets with low SNR.

Discussion and conclusion

The stable two-fold increase in PDE-levels found in DMD patients in muscles with different levels of muscle wasting over 2-year time, including DMD patients as young as 5.5 years-old, suggests that PDE-levels may increase very rapidly early in the disease process and remain elevated thereafter. The low CV values in high and low SNR datasets show that PDE-levels can be accurately and reproducibly quantified in all conditions. Our data confirms the great potential of PDE as a marker for muscle tissue changes in DMD patients.

Large improvement of RF transmission efficiency and reception sensitivity for human in vivo31P MRS imaging using ultrahigh dielectric constant materials at 7T

In vivo³¹P MRS provides a unique and important imaging tool for studying high-energy phosphate metabolism and bioenergetics noninvasively. However, compared to ¹H MRS, ³¹P MRS with a relatively low gyromagnetic ratio (γ) has a lower and limited sensitivity even at ultrahigh field. The proof of concept has been recently demonstrated that the use of high dielectric constant (HDC) materials between RF coil and object sample could increase MRI signal and reduce required RF transmission power for reaching the same RF pulse flip angle in the region of interest. For low-γ MRS applications operated at relatively lower frequency, however, it demands the dielectric materials with a much higher permittivity for achieving optimal performance. We conducted a ³¹P MRS imaging study using ultra-HDC (uHDC; with a relative permittivity of ~1200) material blocks incorporated with an RF volume coil at ultrahigh field of 7.0T. The experimental results from phantom and human calf muscle demonstrate that the uHDC technique significantly enhanced RF magnetic transmit field (B₁⁺) and reception field (B₁^-) and the gain could reach up to two folds in the tissue near the uHDC blocks. The overall results indicate that the incorporation of the uHDC materials having an appropriate permittivity value with a RF coil can significantly increase detection sensitivity and reduces RF transmission power for X-nuclei MRS applications at ultrahigh field. The uHDC technology could provide an efficient, cost-effective engineering solution for achieving high detection sensitivity and concurrently minimizing tissue heating concern for human MRS and MRI applications.

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.

Interleaved multivoxel 31 P MR spectroscopy

PURPOSE

Separate measurements are required when investigating multiple exercising muscles with singlevoxel-localized dynamic 31 P-MRS. With multivoxel spectroscopy, 31 P-MRS time-series spectra are acquired from multiple independent regions during one exercise-recovery experiment with the same time resolution as for singlevoxel measurements.

METHODS

Multiple independently selected volumes were localized using temporally interleaved semi-LASER excitations at 7T. Signal loss caused by mutual saturation from shared excitation or refocusing slices was quantified at partial and full overlap, and potential contamination was investigated in phantom measurements. During an exercise-recovery experiment both gastrocnemius medialis and soleus of two healthy volunteers were measured using multivoxel acquisitions with a total TR of 6 s, while avoiding overlap of excitation slices.

RESULTS

Signal reduction by shared adiabatic refocusing slices selected 1 s after the preceding voxel was between 10% (full overlap) and 20% (half overlap), in a phantom measurement. In vivo data were acquired from both muscles within the same exercise experiment, with 13-18% signal reduction. Spectra show phosphocreatine, inorganic phosphate, adenosine-triposphate, phosphomonoesters, and phosphodiesters.

CONCLUSION

Signal decrease was relatively low compared to the 2-fold increase in information. The approach could help to improve the understanding in metabolic research and is applicable to other organs and nuclei. Magn Reson Med 77:921-927, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Phosphodiester content measured in human liver by in vivo 31 P MR spectroscopy at 7 tesla

Purpose

Phosphorus (³¹P) metabolites are emerging liver disease biomarkers. Of particular interest are phosphomonoester and phosphodiester (PDE) “peaks” that comprise multiple overlapping resonances in ³¹P spectra. This study investigates the effect of improved spectral resolution at 7 Tesla (T) on quantifying hepatic metabolites in cirrhosis.

Methods

Five volunteers were scanned to determine metabolite T₁s. Ten volunteers and 11 patients with liver cirrhosis were scanned at 7T. Liver spectra were acquired in 28 min using a 16‐channel ³¹P array and 3D chemical shift imaging. Concentrations were calculated using γ‐adenosine‐triphosphate (γ‐ATP) = 2.65 mmol/L wet tissue.

Results

T₁ means ± standard deviations: phosphatidylcholine 1.05 ± 0.28 s, nicotinamide‐adenine‐dinucleotide (NAD⁺) 2.0 ± 1.0 s, uridine‐diphosphoglucose (UDPG) 3.3 ± 1.4 s. Concentrations in healthy volunteers: α‐ATP 2.74 ± 0.11 mmol/L wet tissue, inorganic phosphate 2.23 ± 0.20 mmol/L wet tissue, glycerophosphocholine 2.34 ± 0.46 mmol/L wet tissue, glycerophosphoethanolamine 1.50 ± 0.28 mmol/L wet tissue, phosphocholine 1.06 ± 0.16 mmol/L wet tissue, phosphoethanolamine 0.77 ± 0.14 mmol/L wet tissue, NAD⁺ 2.37 ± 0.14 mmol/L wet tissue, UDPG 2.00 ± 0.22 mmol/L wet tissue, phosphatidylcholine 1.38 ± 0.31 mmol/L wet tissue. Inorganic phosphate and phosphatidylcholine concentrations were significantly lower in patients; glycerophosphoethanolamine concentrations were significantly higher (P < 0.05).

Conclusion

We report human in vivo hepatic T₁s for phosphatidylcholine, NAD⁺, and UDPG for the first time at 7T. Our protocol allows high signal‐to‐noise, repeatable measurement of metabolite concentrations in human liver. The splitting of PDE into its constituent peaks at 7T may allow more insight into changes in metabolism. Magn Reson Med 78:2095–2105, 2017. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

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.

Elevated phosphodiester and T2 levels can be measured in the absence of fat infiltration in Duchenne muscular dystrophy patients

Quantitative MRI and MRS are increasingly important as non-invasive outcome measures in therapy development for Duchenne muscular dystrophy (DMD). Many studies have focussed on individual measures such as fat fraction and metabolite levels in relation to age and functionality, but much less attention has been given to how these indices relate to each other. Here, we assessed spatially resolved metabolic changes in leg muscles of DMD patients, and classified muscles according to the degree of fat replacement compared with healthy controls. Quantitative MRI (three-point Dixon and multi-spin echo without fat suppression and a tri-exponential fit) and 2D-CSI ³¹ P MRS scans were obtained from 18 DMD patients and 12 healthy controls using a 3 T and a 7 T MR scanner. Metabolite levels, T₂ values and fat fraction were individually assessed for five lower leg muscles. In muscles with extensive fat replacement, phosphodiester over adenosine triphosphate (PDE/ATP), inorganic phosphate over phosphocreatine, intracellular tissue pH and T₂ were significantly increased compared with healthy controls. In contrast, in muscles without extensive fat replacement, only PDE/ATP and T₂ values were significantly elevated. Overall, our results show that PDE levels and T₂ values increase prior to the occurrence of fat replacement and remain elevated in later stages of the disease. This suggests that these individual measures could not only function as early markers for muscle damage but also reflect potentially reversible pathology in the more advanced stages.

Adiabatic excitation for 31 P MR spectroscopy in the human heart at 7 T: A feasibility study

Purpose

Phosphorus magnetic resonance spectroscopy (³¹P‐MRS) provides a unique tool for assessing cardiac energy metabolism, often quantified using the phosphocreatine (PCr)/adenosine triphosphate (ATP) ratio. Surface coils are typically used for excitation for ³¹P‐MRS, but they create an inhomogeneous excitation field across the myocardium, producing undesirable, spatially varying partial saturation. Therefore, we implemented adiabatic excitation in a 3D chemical shift imaging (CSI) sequence for cardiac ³¹P‐MRS at 7 Tesla (T).

Methods

We optimized an adiabatic half passage pulse with bandwidth sufficient to excite PCr and γ‐ATP together. In addition, the CSI sequence was modified to allow interleaved excitation of PCr and γ‐ATP, then 2,3‐DPG, to enable PCr/ATP determination with blood correction. Nine volunteers were scanned at 2 transmit voltages to confirm that measured PCr/ATP was independent of B1+ (i.e. over the adiabatic threshold). Six septal voxels were evaluated for each volunteer.

Results

Phantom experiments showed that adiabatic excitation can be reached at the depth of the heart using our pulse. The mean evaluated cardiac PCr/ATP ratio from all 9 volunteers corrected for blood signal was 2.14 ± 0.16. Comparing the two acquisitions with different voltages resulted in a minimal mean difference of −0.005.

Conclusion

Adiabatic excitation is possible in the human heart at 7 T, and gives consistent PCr/ATP ratios. Magn Reson Med 78:1667–1673, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

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

Phosphorus MRSI (³¹ P-MRSI) using a spiral-trajectory readout at 7 T was developed for high temporal resolution mapping of the mitochondrial capacity of exercising human skeletal muscle. The sensitivity and localization accuracy of the method was investigated in phantoms. In vivo performance was assessed in 12 volunteers, who performed a plantar flexion exercise inside a whole-body 7 T MR scanner using an MR-compatible ergometer and a surface coil. In five volunteers the knee was flexed (~60°) to shift the major workload from the gastrocnemii to the soleus muscle. Spiral-encoded MRSI provided 16-25 times faster mapping with a better point spread function than elliptical phase-encoded MRSI with the same matrix size. The inevitable trade-off for the increased temporal resolution was a reduced signal-to-noise ratio, but this was acceptable. The phosphocreatine (PCr) depletion caused by exercise at 0° knee angulation was significantly higher in both gastrocnemii than in the soleus (i.e. 64.8 ± 19.6% and 65.9 ± 23.6% in gastrocnemius lateralis and medialis versus 15.3 ± 8.4% in the soleus). Spiral-encoded ³¹ P-MRSI is a powerful tool for dynamic mapping of exercising muscle oxidative metabolism, including localized assessment of PCr concentrations, pH and maximal oxidative flux with high temporal and spatial resolution.

Efficient 31 P band inversion transfer approach for measuring creatine kinase activity, ATP synthesis, and molecular dynamics in the human brain at 7 T

PURPOSE

To develop an efficient ³¹ P magnetic resonance spectroscopy (MRS) method for measuring creatine kinase (CK) activity, adenosine triphosphate (ATP) synthesis, and motion dynamics in the human brain at 7 Tesla (T).

METHODS

Three band inversion modules differing in center frequency were used to induce magnetization transfer (MT) effect in three exchange pathways: (i) CK-mediated reaction PCr → γ-ATP; (ii) de novo ATP synthesis Pi → γ-ATP; and (iii) ATP intramolecular ³¹ P-³¹ P cross-relaxation γ-(α-) ↔ β-ATP. The resultant MT data were analyzed using a 5-pool model in the format of magnetization matrix according to Bloch-McConnell-Solomon formalism.

RESULTS

With a repetition time (TR) of 4 s, the scan time for each module was approximately 8 min. The rate constants were k_{PCr → γATP} 0.38 ± 0.02 s^-1 , k_{Pi → γATP} 0.19 ± 0.02 s^-1 , and σ_{γ(α) ↔ βATP} 0.19 ± 0.04 s^-1 , corresponding to ATP rotation correlation time τ_c (0.8 ± 0.2) ·10^-7 s. The T₁ relaxation times were Pi 7.26 ± 1.76 s, PCr 5.99 ± 0.58 s, γ-ATP 0.98 ± 0.07 s, α-ATP 0.95 ± 0.04 s, and β-ATP 0.68 ± 0.03 s.

CONCLUSION

Short-TR band inversion modules provide a time-efficient way of measuring brain ATP metabolism and could be useful in studying metabolic disorders in brain diseases. Magn Reson Med 78:1657-1666, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Clinical applications at ultrahigh field (7 T). Where does it make the difference?

Presently, three major MR vendors provide commercial 7-T units for clinical research under ethical permission, with the number of operating 7-T systems having increased to over 50. This rapid increase indicates the growing interest in ultrahigh-field MRI because of improved clinical results with regard to morphological as well as functional and metabolic capabilities. As the signal-to-noise ratio scales linearly with the field strength (B0 ) of the scanner, the most obvious application at 7 T is to obtain higher spatial resolution in the brain, musculoskeletal system and breast. Of specific clinical interest for neuro-applications is the cerebral cortex at 7 T, for the detection of changes in cortical structure as a sign of early dementia, as well as for the visualization of cortical microinfarcts and cortical plaques in multiple sclerosis. In the imaging of the hippocampus, even subfields of the internal hippocampal anatomy and pathology can be visualized with excellent resolution. The dynamic and static blood oxygenation level-dependent contrast increases linearly with the field strength, which significantly improves the pre-surgical evaluation of eloquent areas before tumor removal. Using susceptibility-weighted imaging, the plaque-vessel relationship and iron accumulation in multiple sclerosis can be visualized for the first time. Multi-nuclear clinical applications, such as sodium imaging for the evaluation of repair tissue quality after cartilage transplantation and (31) P spectroscopy for the differentiation between non-alcoholic benign liver disease and potentially progressive steatohepatitis, are only possible at ultrahigh fields. Although neuro- and musculoskeletal imaging have already demonstrated the clinical superiority of ultrahigh fields, whole-body clinical applications at 7 T are still limited, mainly because of the lack of suitable coils. The purpose of this article was therefore to review the clinical studies that have been performed thus far at 7 T, compared with 3 T, as well as those studies performed at 7 T that cannot be routinely performed at 3 T. Copyright © 2015 John Wiley & Sons, Ltd.

A simple approach to evaluate the kinetic rate constant for ATP synthesis in resting human skeletal muscle at 7 T

Inversion transfer (IT) is a well-established technique with multiple attractive features for analysis of kinetics. However, its application in measurement of ATP synthesis rate in vivo has lagged behind the more common saturation transfer (ST) techniques. One well-recognized issue with IT is the complexity of data analysis in comparison with much simpler analysis by ST. This complexity arises, in part, because the γ-ATP spin is involved in multiple chemical reactions and magnetization exchanges, whereas Pi is involved in a single reaction, Pi → γ-ATP. By considering the reactions involving γ-ATP only as a lumped constant, the rate constant for the reaction of physiological interest, kPi→γATP , can be determined. Here, we present a new IT data analysis method to evaluate kPi→γATP using data collected from resting human skeletal muscle at 7 T. The method is based on the basic Bloch-McConnell equation, which relates kPi→γATP to m˙Pi, the rate of Pi magnetization change. The kPi→γATP value is accessed from m˙Pi data by more familiar linear correlation approaches. For a group of human subjects (n = 15), the kPi→γATP value derived for resting calf muscle was 0.066 ± 0.017 s(-1) , in agreement with literature-reported values. In this study we also explored possible time-saving strategies to speed up data acquisition for kPi→γATP evaluation using simulations. The analysis indicates that it is feasible to carry out a (31) P IT experiment in about 10 min or less at 7 T with reasonable outcome in kPi→γATP variance for measurement of ATP synthesis in resting human skeletal muscle. We believe that this new IT data analysis approach will facilitate the wide acceptance of IT to evaluate ATP synthesis rate in vivo. Copyright © 2015 John Wiley & Sons, Ltd.