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

Sequence building block for magnetic resonance spectroscopy on Siemens VE-series scanners

We present a sequence building block (SBB) that embeds magnetic resonance spectroscopy (MRS) into another sequence on the Siemens VE platform without any custom hardware. This enables dynamic studies such as functional MRS (fMRS), dynamic shimming and frequency correction, and acquisition of navigator images for motion correction. The SBB supports nonlocalised spectroscopy (free induction decay), STimulated Echo Acquisition Mode single voxel spectroscopy, and 1D, 2D and 3D phase-encoded chemical shift imaging. It can embed 1H or X-nuclear MRS into a 1H sequence; and 1H-MRS into an X-nuclear sequence. We demonstrate integration into the vendor's gradient-recalled echo sequence. We acquire test data in phantoms with three coils (31P/1H, 13C/1H and 2H/1H) and in two volunteers on a 7-T Terra MRI scanner. Fifteen lines of code are required to insert the SBB into a sequence. Spectra and images are acquired successfully in all cases in phantoms, and in human abdomen and calf muscle. Phantom comparison of signal-to-noise ratio and linewidth showed that the SBB has negligible effects on image and spectral quality, except that it sometimes produces a nuclear Overhauser effect (NOE) signal enhancement for multinuclear applications in line with conventional 1H NOE pulses. Our new SBB embeds MRS into a host imaging or spectroscopy sequence in 15 lines of code. It allows homonuclear and heteronuclear interleaving. The package is available through the standard C2P procedure. We hope this will lower the barrier for entry to studies applying dynamic fMRS and for online motion correction and B0-shim updating.

Measurement of metabolite levels and treatment-induced changes in hepatic metastases of gastro-esophageal cancer using 7-T phosphorus magnetic resonance spectroscopic imaging

Methods for early treatment response evaluation to systemic therapy of liver metastases are lacking. Tumor tissue often exhibits an increased ratio of phosphomonoesters to phosphodiesters (PME/PDE), which can be noninvasively measured by phosphorus magnetic resonance spectroscopy (31P MRS), and may be a marker for early therapy response assessment in liver metastases. However, with commonly used 31P surface coils for liver 31P MRS, the liver is not fully covered, and metastases may be missed. The objective of this study was to demonstrate the feasibility of 31P MRS imaging (31P MRSI) with full liver coverage to assess 31P metabolite levels and chemotherapy-induced changes in liver metastases of gastro-esophageal cancer, using a 31P whole-body birdcage transmit coil in combination with a 31P body receive array at 7 T. 3D 31P MRSI data were acquired in two patients with hepatic metastases of esophageal cancer, before the start of chemotherapy and after 2 (and 9 in patient 2) weeks of chemotherapy. 3D 31P MRSI acquisitions were performed using an integrated 31P whole-body transmit coil in combination with a 16-channel body receive array at 7 T, with a field of view covering the full abdomen and a nominal voxel size of 20-mm isotropic. From the 31P MRSI data, 12 31P metabolite signals were quantified. Prior to chemotherapy initiation, both PMEs, that is, phosphocholine (PC) and phosphoethanolamine (PE), were significantly higher in all metastases compared with the levels previously determined in the liver of healthy volunteers. After 2 weeks of chemotherapy, PC and PE levels remained high or even increased further, resulting in increased PME/PDE ratios compared with healthy liver tissue, in correspondence with the clinical assessment of progressive disease after 2 months of chemotherapy. The suggested approach may present a viable tool for early therapy (non)response assessment of tumor metabolism in patients with liver metastases.

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.

Repeatability of deuterium metabolic imaging of healthy volunteers at 3 T

BACKGROUND

Magnetic resonance (MR) imaging of deuterated glucose, termed deuterium metabolic imaging (DMI), is emerging as a biomarker of pathway-specific glucose metabolism in tumors. DMI is being studied as a useful marker of treatment response in a scan-rescan scenario. This study aims to evaluate the repeatability of brain DMI.

METHODS

A repeatability study was performed in healthy volunteers from December 2022 to March 2023. The participants consumed 75 g of [6,6'-2H2]glucose. The delivery of 2H-glucose to the brain and its conversion to 2H-glutamine + glutamate, 2H-lactate, and 2H-water DMI was imaged at baseline and at 30, 70, and 120 min. DMI was performed using MR spectroscopic imaging on a 3-T system equipped with a 1H/2H-tuned head coil. Coefficients of variation (CoV) were computed for estimation of repeatability and between-subject variability. In a set of exploratory analyses, the variability effects of region, processing, and normalization were estimated.

RESULTS

Six male participants were recruited, aged 34 ± 6.5 years (mean ± standard deviation). There was 42 ± 2.7 days between sessions. Whole-brain levels of glutamine + glutamate, lactate, and glucose increased to 3.22 ± 0.4 mM, 1.55 ± 0.3 mM, and 3 ± 0.7 mM, respectively. The best signal-to-noise ratio and repeatability was obtained at the 120-min timepoint. Here, the within-subject whole-brain CoVs were -10% for all metabolites, while the between-subject CoVs were -20%.

CONCLUSIONS

DMI of glucose and its downstream metabolites is feasible and repeatable on a clinical 3 T system.

TRIAL REGISTRATION

ClinicalTrials.gov, NCT05402566 , registered the 25th of May 2022.

RELEVANCE STATEMENT

Brain deuterium metabolic imaging of healthy volunteers is repeatable and feasible at clinical field strengths, enabling the study of shifts in tumor metabolism associated with treatment response.

KEY POINTS

• Deuterium metabolic imaging is an emerging tumor biomarker with unknown repeatability.  • The repeatability of deuterium metabolic imaging is on par with FDG-PET.  • The study of deuterium metabolic imaging in clinical populations is feasible.

Linewidth-related bias in modelled concentration estimates from GABA-edited 1H-MRS

J-difference-edited MRS is widely used to study GABA in the human brain. Editing for low-concentration target molecules (such as GABA) typically exhibits lower signal-to-noise ratio (SNR) than conventional non-edited MRS, varying with acquisition region, volume and duration. Moreover, spectral lineshape may be influenced by age-, pathology-, or brain-region-specific effects of metabolite T2, or by task-related blood-oxygen level dependent (BOLD) changes in functional MRS contexts. Differences in both SNR and lineshape may have systematic effects on concentration estimates derived from spectral modelling. The present study characterises the impact of lineshape and SNR on GABA+ estimates from different modelling algorithms: FSL-MRS, Gannet, LCModel, Osprey, spant and Tarquin. Publicly available multi-site GABA-edited data (222 healthy subjects from 20 sites; conventional MEGA-PRESS editing; TE = 68 ms) were pre-processed with a standardised pipeline, then filtered to apply controlled levels of Lorentzian and Gaussian linebroadening and SNR reduction. Increased Lorentzian linewidth was associated with a 2-5% decrease in GABA+ estimates per Hz, observed consistently (albeit to varying degrees) across datasets and most algorithms. Weaker, often opposing effects were observed for Gaussian linebroadening. Variations are likely caused by differing baseline parametrization and lineshape constraints between models. Effects of linewidth on other metabolites (e.g., Glx and tCr) varied, suggesting that a linewidth confound may persist after scaling to an internal reference. These findings indicate a potentially significant confound for studies where linewidth may differ systematically between groups or experimental conditions, e.g. due to T2 differences between brain regions, age, or pathology, or varying T2* due to BOLD-related changes. We conclude that linewidth effects need to be rigorously considered during experimental design and data processing, for example by incorporating linewidth into statistical analysis of modelling outcomes or development of appropriate lineshape matching algorithms.

Multi-channel magnetic resonance spectroscopy graphical user interface (McMRSGUI)

This work introduces an open-sourced graphical user interface (GUI) software enabling the combination of multi-channel magnetic resonance spectroscopy data with different literature-based methods for the improvement of the quality and reliability of combined spectra. The multi-channel magnetic resonance spectroscopy graphical user interface (McMRSGUI) is a MATLAB-based spectroscopy processing GUI equipped to load multi-channel MRS data, pre-process, combine, and export combined data for evaluation with open-source quantification software (jMRUI). A literature-based, decision-tree process was incorporated into the combination type selection to serve as a guide to minimize spectral distortion in selecting between weighting methods. Multi-channel, simulated spectra were combined with the different combination techniques and evaluated for spectral distortion to validate the code. The incorporation of the combination methods into a single processing software enables multi-channel magnetic resonance spectroscopy (MRS) data to be combined and compared for improved spectral quality with little user knowledge of combination techniques. Through the spectral peak distortion simulation of the combination methods, combined signal-to-noise ratio (SNR) values from the literature were verified. The spectral peak distortion simulation provides a secondary tool for researchers to estimate the spectral SNR levels when spectral distortion could occur and use this knowledge to further guide the selection of their combination technique. The McMRSGUI provides a software toolkit for evaluating multi-channel MRS data and their combination. Simulations evaluating spectral distortion at different noise levels were performed for each combination method to validate the GUI and demonstrate a method for researchers to assess the combined SNR levels at which they could be introducing spectral distortion.

Compartment-based reconstruction of acquisition-weighted 31P cardiac MRSI reduces sensitivity to cardiac motion and scan planning

Motivation: 31P magnetic resonance spectroscopic imaging (31P MRSI) is a powerful technique for investigating the metabolic effects of treatments for heart failure in vivo, allowing a better understanding of their mechanism of action in patient cohorts. Unfortunately, cardiac 31P MRSI is fundamentally limited by low SNR, which leads to compromises in acquisition, such as no cardiac or respiratory gating or low spatial resolution, in order to achieve reasonable scan times. Spectroscopy with linear algebra modeling (SLAM) reconstruction may be able to address these challenges and therefore improve repeatability by incorporating a segmented localizer into the reconstruction. Methods: Six healthy volunteers were scanned twice in a test-retest procedure to allow quantification of repeatability. Each scan consisted of anatomical localizers and two acquisition-weighted (AW) 31P MRSI acquisitions, which were acquired with and without cardiac gating. Five patients with heart failure with a preserved ejection fraction were then scanned with the same 31P MRSI sequence without cardiac gating. All 31P MRSI datasets were reconstructed with both conventional Fourier transform (FT)-based reconstruction and SLAM reconstruction, which were compared statistically. The effect of shifting the 31P MRSI acquisition field of view was also investigated. Results: In the healthy volunteer cohort, the spectral fit of the SLAM reconstructions had significantly improved Cramer-Rao lower bounds (CRLBs) compared to the FT-based reconstruction of non-cardiac gated data, as well as improved coefficients of variability and repeatability. The SLAM reconstruction found a significant difference in the PCr/ATP ratio between the healthy volunteer and patient cohorts, which the FT-based reconstruction did not find. Furthermore, the SLAM reconstruction was less influenced by the placement of the field of view (FOV) of the 31P MRSI acquisition in post hoc analysis. Discussion: The experimental benefits of the SLAM reconstruction for AW data were demonstrated by the improvements in fit confidence and repeatability seen in the healthy volunteer cohort and post hoc FOV analysis. The benefit of SLAM reconstruction of AW data for clinical studies was then illustrated by the patient cohort, which suggested improved sensitivity to clinically significant changes in the PCr/ATP ratio.

Quo Vadis Hyperpolarized 13C MRI?

Over the last two decades, hyperpolarized 13C MRI has gained significance in both preclinical and clinical studies, hereby relying on technologies like PHIP-SAH (ParaHydrogen-Induced Polarization-Side Arm Hydrogenation), SABRE (Signal Amplification by Reversible Exchange), and dDNP (dissolution Dynamic Nuclear Polarization), with dDNP being applied in humans. A clinical dDNP polarizer has enabled studies across 24 sites, despite challenges like high cost and slow polarization. Parahydrogen-based techniques like SABRE and PHIP offer faster, more cost-efficient alternatives but require molecule-specific optimization. The focus has been on imaging metabolism of hyperpolarized probes, which requires long T1, high polarization and rapid contrast generation. Efforts to establish novel probes, improve acquisition techniques and enhance data analysis methods including artificial intelligence are ongoing. Potential clinical value of hyperpolarized 13C MRI was demonstrated primarily for treatment response assessment in oncology, but also in cardiology, nephrology, hepatology and CNS characterization. In this review on biomedical hyperpolarized 13C MRI, we summarize important and recent advances in polarization techniques, probe development, acquisition and analysis methods as well as clinical trials. Starting from those we try to sketch a trajectory where the field of biomedical hyperpolarized 13C MRI might go.

Role of Cardiac Energetics in Aortic Stenosis Disease Progression: Identifying the High-risk Metabolic Phenotype

BACKGROUND

Severe aortic stenosis (AS) is associated with left ventricular (LV) hypertrophy and cardiac metabolic alterations with evidence of steatosis and impaired myocardial energetics. Despite this common phenotype, there is an unexplained and wide individual heterogeneity in the degree of hypertrophy and progression to myocardial fibrosis and heart failure. We sought to determine whether the cardiac metabolic state may underpin this variability.

METHODS

We recruited 74 asymptomatic participants with AS and 13 healthy volunteers. Cardiac energetics were measured using phosphorus spectroscopy to define the myocardial phosphocreatine to adenosine triphosphate ratio. Myocardial lipid content was determined using proton spectroscopy. Cardiac function was assessed by cardiovascular magnetic resonance cine imaging.

RESULTS

Phosphocreatine/adenosine triphosphate was reduced early and significantly across the LV wall thickness quartiles (Q2, 1.50 [1.21-1.71] versus Q1, 1.64 [1.53-1.94]) with a progressive decline with increasing disease severity (Q4, 1.48 [1.18-1.70]; P=0.02). Myocardial triglyceride content levels were overall higher in all the quartiles with a significant increase seen across the AV pressure gradient quartiles (Q2, 1.36 [0.86-1.98] versus Q1, 1.03 [0.81-1.56]; P=0.034). While all AS groups had evidence of subclinical LV dysfunction with impaired strain parameters, impaired systolic longitudinal strain was related to the degree of energetic impairment (r=0.219; P=0.03). Phosphocreatine/adenosine triphosphate was not only an independent predictor of LV wall thickness (r=-0.20; P=0.04) but also strongly associated with myocardial fibrosis (r=-0.24; P=0.03), suggesting that metabolic changes play a role in disease progression. The metabolic and functional parameters showed comparable results when graded by clinical severity of AS.

CONCLUSIONS

A gradient of myocardial energetic deficit and steatosis exists across the spectrum of hypertrophied AS hearts, and these metabolic changes precede irreversible LV remodeling and subclinical dysfunction. As such, cardiac metabolism may play an important and potentially causal role in disease progression.

Assessment of the sensitivity of 2 H MR spectroscopy measurements of [2,3-2 H2 ]fumarate metabolism for detecting tumor cell death

Imaging the metabolism of [2,3-2 H2 ]fumarate to produce malate can be used to detect tumor cell death post-treatment. Here, we assess the sensitivity of the technique for detecting cell death by lowering the concentration of injected [2,3-2 H2 ]fumarate and by varying the extent of tumor cell death through changes in drug concentration. Mice were implanted subcutaneously with human triple negative breast cancer cells (MDA-MB-231) and injected with 0.1, 0.3, and 0.5 g/kg [2,3-2 H2 ]fumarate before and after treatment with a multivalent TRAlL-R2 agonist (MEDI3039) at 0.1, 0.4, and 0.8 mg/kg. Tumor conversion of [2,3-2 H2 ]fumarate to [2,3-2 H2 ]malate was assessed from a series of 13 spatially localized 2 H MR spectra acquired over 65 min using a pulse-acquire sequence with a 2-ms BIR4 adiabatic excitation pulse. Tumors were then excised and stained for histopathological markers of cell death: cleaved caspase 3 (CC3) and DNA damage (terminal deoxynucleotidyl transferase dUTP nick end labeling [TUNEL]). The rate of malate production and the malate/fumarate ratio plateaued at tumor fumarate concentrations of 2 mM, which were obtained with injected [2,3-2 H2 ]fumarate concentrations of 0.3 g/kg and above. Tumor malate concentration and the malate/fumarate ratio increased linearly with the extent of cell death determined histologically. At an injected [2,3-2 H2 ]fumarate concentration of 0.3 g/kg, 20% CC3 staining corresponded to a malate concentration of 0.62 mM and a malate/fumarate ratio of 0.21. Extrapolation indicated that there would be no detectable malate at 0% CC3 staining. The use of low and nontoxic fumarate concentrations and the production of [2,3-2 H2 ]malate at concentrations that are within the range that can be detected clinically suggest this technique could translate to the clinic.

Retained Metabolic Flexibility of the Failing Human Heart

BACKGROUND

The failing heart is traditionally described as metabolically inflexible and oxygen starved, causing energetic deficit and contractile dysfunction. Current metabolic modulator therapies aim to increase glucose oxidation to increase oxygen efficiency of adenosine triphosphate production, with mixed results.

METHODS

To investigate metabolic flexibility and oxygen delivery in the failing heart, 20 patients with nonischemic heart failure with reduced ejection fraction (left ventricular ejection fraction 34.9±9.1) underwent separate infusions of insulin+glucose infusion (I+G) or Intralipid infusion. We used cardiovascular magnetic resonance to assess cardiac function and measured energetics using phosphorus-31 magnetic resonance spectroscopy. To investigate the effects of these infusions on cardiac substrate use, function, and myocardial oxygen uptake (MVo2), invasive arteriovenous sampling and pressure-volume loops were performed (n=9).

RESULTS

At rest, we found that the heart had considerable metabolic flexibility. During I+G, cardiac glucose uptake and oxidation were predominant (70±14% total energy substrate for adenosine triphosphate production versus 17±16% for Intralipid; P=0.002); however, no change in cardiac function was seen relative to basal conditions. In contrast, during Intralipid infusion, cardiac long-chain fatty acid (LCFA) delivery, uptake, LCFA acylcarnitine production, and fatty acid oxidation were all increased (LCFA 73±17% of total substrate versus 19±26% total during I+G; P=0.009). Myocardial energetics were better with Intralipid compared with I+G (phosphocreatine/adenosine triphosphate 1.86±0.25 versus 2.01±0.33; P=0.02), and systolic and diastolic function were improved (LVEF 34.9±9.1 baseline, 33.7±8.2 I+G, 39.9±9.3 Intralipid; P<0.001). During increased cardiac workload, LCFA uptake and oxidation were again increased during both infusions. There was no evidence of systolic dysfunction or lactate efflux at 65% maximal heart rate, suggesting that a metabolic switch to fat did not cause clinically meaningful ischemic metabolism.

CONCLUSIONS

Our findings show that even in nonischemic heart failure with reduced ejection fraction with severely impaired systolic function, significant cardiac metabolic flexibility is retained, including the ability to alter substrate use to match both arterial supply and changes in workload. Increasing LCFA uptake and oxidation is associated with improved myocardial energetics and contractility. Together, these findings challenge aspects of the rationale underlying existing metabolic therapies for heart failure and suggest that strategies promoting fatty acid oxidation may form the basis for future therapies.

Assessment of Cardiac Energy Metabolism, Function, and Physiology in Patients With Heart Failure Taking Empagliflozin: The Randomized, Controlled EMPA-VISION Trial

BACKGROUND

Sodium-glucose co-transporter 2 inhibitors (SGLT2i) have emerged as a paramount treatment for patients with heart failure (HF), irrespective of underlying reduced or preserved ejection fraction. However, a definite cardiac mechanism of action remains elusive. Derangements in myocardial energy metabolism are detectable in all HF phenotypes, and it was proposed that SGLT2i may improve energy production. The authors aimed to investigate whether treatment with empagliflozin leads to changes in myocardial energetics, serum metabolomics, and cardiorespiratory fitness.

METHODS

EMPA-VISION (Assessment of Cardiac Energy Metabolism, Function and Physiology in Patients With Heart Failure Taking Empagliflozin) is a prospective, randomized, double-blind, placebo-controlled, mechanistic trial that enrolled 72 symptomatic patients with chronic HF with reduced ejection fraction (HFrEF; n=36; left ventricular ejection fraction ≤40%; New York Heart Association class ≥II; NT-proBNP [N-terminal pro-B-type natriuretic peptide] ≥125 pg/mL) and HF with preserved ejection fraction (HFpEF; n=36; left ventricular ejection fraction ≥50%; New York Heart Association class ≥II; NT-proBNP ≥125 pg/mL). Patients were stratified into respective cohorts (HFrEF versus HFpEF) and randomly assigned to empagliflozin (10 mg; n=35: 17 HFrEF and 18 HFpEF) or placebo (n=37: 19 HFrEF and 18 HFpEF) once daily for 12 weeks. The primary end point was a change in the cardiac phosphocreatine:ATP ratio (PCr/ATP) from baseline to week 12, determined by phosphorus magnetic resonance spectroscopy at rest and during peak dobutamine stress (65% of age-maximum heart rate). Mass spectrometry on a targeted set of 19 metabolites was performed at baseline and after treatment. Other exploratory end points were investigated.

RESULTS

Empagliflozin treatment did not change cardiac energetics (ie, PCr/ATP) at rest in HFrEF (adjusted mean treatment difference [empagliflozin - placebo], -0.25 [95% CI, -0.58 to 0.09]; P=0.14) or HFpEF (adjusted mean treatment difference, -0.16 [95% CI, -0.60 to 0.29]; P=0.47]. Likewise, there were no changes in PCr/ATP during dobutamine stress in HFrEF (adjusted mean treatment difference, -0.13 [95% CI, -0.35 to 0.09]; P=0.23) or HFpEF (adjusted mean treatment difference, -0.22 [95% CI, -0.66 to 0.23]; P=0.32). No changes in serum metabolomics or levels of circulating ketone bodies were observed.

CONCLUSIONS

In patients with either HFrEF or HFpEF, treatment with 10 mg of empagliflozin once daily for 12 weeks did not improve cardiac energetics or change circulating serum metabolites associated with energy metabolism when compared with placebo. Based on our results, it is unlikely that enhancing cardiac energy metabolism mediates the beneficial effects of SGLT2i in HF.

REGISTRATION

URL: https://www.

CLINICALTRIALS

gov; Unique identifier: NCT03332212.

Deuterium echo-planar spectroscopic imaging (EPSI) in the human liver in vivo at 7 T

PURPOSE

To demonstrate the feasibility of deuterium echo-planar spectroscopic imaging (EPSI) to accelerate 3D deuterium metabolic imaging in the human liver at 7 T.

METHODS

A deuterium EPSI sequence, featuring a Hamming-weighted k-space acquisition pattern for the phase-encoding directions, was implemented. Three-dimensional deuterium EPSI and conventional MRSI were performed on a water/acetone phantom and in vivo in the human liver at natural abundance. Moreover, in vivo deuterium EPSI measurements were acquired after oral administration of deuterated glucose. The effect of acquisition time on SNR was evaluated by retrospectively reducing the number of averages.

RESULTS

The SNR of natural abundance deuterated water signal in deuterium EPSI was 6.5% and 5.9% lower than that of MRSI in the phantom and in vivo experiments, respectively. In return, the acquisition time of in vivo EPSI data could be reduced retrospectively to 2 min, beyond the minimal acquisition time of conventional MRSI (of 20 min in this case), while still leaving sufficient SNR. Three-dimensional deuterium EPSI, after administration of deuterated glucose, enabled monitoring of hepatic glucose dynamics with full liver coverage, a spatial resolution of 20 mm isotropic, and a temporal resolution of 9 min 50 s, which could retrospectively be shortened to 2 min.

CONCLUSION

In this work, we demonstrate the feasibility of accelerated 3D deuterium metabolic imaging of the human liver using deuterium EPSI. The acceleration obtained with EPSI can be used to increase temporal and/or spatial resolution, which will be valuable to study tissue metabolism of deuterated compounds over time.

Efficacy and tolerability of an endogenous metabolic modulator (AXA1125) in fatigue-predominant long COVID: a single-centre, double-blind, randomised controlled phase 2a pilot study

Background

'Long COVID' describes persistent symptoms, commonly fatigue, lasting beyond 12 weeks following SARS-CoV-2 infection. Potential causes include reduced mitochondrial function and cellular bioenergetics. AXA1125 has previously increased β-oxidation and improved bioenergetics in preclinical models along with certain clinical conditions, and therefore may reduce fatigue associated with Long COVID. We aimed to assess the efficacy, safety and tolerability of AXA1125 in Long COVID.

Methods

Patients with fatigue-dominant Long COVID were recruited in this single-centre, double-blind, randomised controlled phase 2a pilot study completed in the UK. Patients were randomly assigned (1:1) using an Interactive Response Technology to receive either AXA1125 or matching placebo in a clinical-based setting. Each dose (33.9 g) of AXA1125 or placebo was administered orally in a liquid suspension twice daily for four weeks with a two-week follow-up period. The primary endpoint was the mean change from baseline to day 28 in the phosphocreatine (PCr) recovery rate following moderate exercise, assessed by ³¹P-magnetic resonance spectroscopy (MRS). All patients were included in the intention to treat analysis. This trial was registered at ClinicalTrials.gov, NCT05152849.

Findings

Between December 15th 2021, and May 23th 2022, 60 participants were screened, and 41 participants were randomised and included in the final analysis. Changes in skeletal muscle phosphocreatine recovery time constant (τ_PCr) and 6-min walk test (6MWT) did not significantly differ between treatment (n = 21) and placebo group (n = 20). However, treatment with AXA1125 was associated with significantly reduced day 28 Chalder Fatigue Questionnaire [CFQ-11] fatigue score when compared with placebo (least squares mean difference [LSMD] -4.30, 95% confidence interval (95% CI) -7.14, -1.47; P = 0.0039). Eleven (52.4%, AXA1125) and four (20.0%, placebo) patients reported treatment-emergent adverse events; none were serious or led to treatment discontinuation.

Interpretation

Although treatment with AXA1125 did not improve the primary endpoint (τ_PCr-measure of mitochondrial respiration), when compared to placebo, there were significant improvements in fatigue-based symptoms among patients living with Long COVID following a four-week treatment period. Further multicentre studies are needed to validate our findings in a larger cohort of patients with fatigue-dominant Long COVID.

Funding

Axcella Therapeutics.

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.

Compartment-based reconstruction of 3D acquisition-weighted 31 P cardiac magnetic resonance spectroscopic imaging at 7 T: A reproducibility study

Even at 7 T, cardiac 31 P magnetic resonance spectroscopic imaging (MRSI) is fundamentally limited by low signal-to-noise ratio (SNR), leading to long scan times and poor temporal and spatial resolutions. Compartment-based reconstruction algorithms such as magnetic resonance spectroscopy with linear algebraic modeling (SLAM) and spectral localization by imaging (SLIM) may improve SNR or reduce scan time without changes to acquisition. Here, we compare the repeatability and SNR performance of these compartment-based methods, applied to three different acquisition schemes at 7 T. Twelve healthy volunteers were scanned twice. Each scan session consisted of a 6.5-min 3D acquisition-weighted (AW) cardiac 31 P phase encode-based MRSI acquisition and two 6.5-min truncated k-space acquisitions with increased averaging (4 × 4 × 4 central k-space phase encodes and fractional SLAM [fSLAM] optimized k-space phase encodes). Spectra were reconstructed using (i) AW Fourier reconstruction; (ii) AW SLAM; (iii) AW SLIM; (iv) 4 × 4 × 4 SLAM; (v) 4 × 4 × 4 SLIM; and (vi) fSLAM acquisition-reconstruction combinations. The phosphocreatine-to-adenosine triphosphate (PCr/ATP) ratio, the PCr SNR, and spatial response functions were computed, in addition to coefficients of reproducibility and variability. Using the compartment-based reconstruction algorithms with the AW 31 P acquisition resulted in a significant increase in SNR compared with previously published Fourier-based MRSI reconstruction methods while maintaining the measured PCr/ATP ratio and improving interscan reproducibility. The alternative acquisition strategies with truncated k-space performed no better than the common AW approach. Compartment-based spectroscopy approaches provide an attractive reconstruction method for cardiac 31 P spectroscopy at 7 T, improving reproducibility and SNR without the need for a dedicated k-space sampling strategy.

Detection and alterations of acetylcarnitine (AC) in human liver by 1 H MRS at 3T after supplementation with l-carnitine

PURPOSE

Acetylcarnitine can be assessed in vivo using proton MRS (1 H-MRS) with long TEs and this has been previously applied successfully in muscle. The aim of this study was to evaluate a 1 H-MRS technique for liver acetylcarnitine quantification in healthy humans before and after l-carnitine supplementation.

METHOD

Baseline acetylcarnitine levels were quantified using a STEAM sequence with prolonged TE in 15 healthy adults. Using STEAM with four different TEs was evaluated in phantoms. To assess reproducibility of the measurements, five of the participants had repeated 1 H-MRS without receiving l-carnitine supplementation. To determine if liver acetylcarnitine could be changed after l-carnitine supplementation, acetylcarnitine was quantified 2 h after intravenous l-carnitine supplementation (50 mg/kg body weight) in the other 10 participants. Hepatic lipids were also quantified from the 1 H-MRS spectra.

RESULTS

There was good separation between the acetylcarnitine and fat in the phantoms using TE = 100 ms. Hepatic acetylcarnitine levels were reproducible (coefficient of reproducibility = 0.049%) and there was a significant (p < 0.001) increase in the relative abundance after a single supplementation of l-carnitine. Hepatic allylic, methyl, and methylene peaks were not altered by l-carnitine supplementation in healthy volunteers.

CONCLUSION

Our results demonstrate that our 1 H-MRS technique could be used to measure acetylcarnitine in the liver and detect changes following intravenous supplementation in healthy adults despite the presence of lipids. Our techniques should be explored further in the study of fatty liver disease, where acetylcarnitine is suggested to be altered due to hepatic inflexibilities.

Deuterium body array for the simultaneous measurement of hepatic and renal glucose metabolism and gastric emptying with dynamic 3D deuterium metabolic imaging at 7 T

Deuterium metabolic imaging (DMI) is a novel noninvasive method to assess tissue metabolism and organ (patho)physiology in vivo using deuterated substrates, such as [6,6'-² H₂ ]-glucose. The liver and kidneys play a central role in whole-body glucose homeostasis, and in type 2 diabetes, both hepatic and renal glucose metabolism are dysregulated. Diabetes is also associated with gastric emptying abnormalities. In this study, we developed a four-channel ² H transmit/receive body array coil for DMI in the human abdomen at 7 T and assessed its performance. In addition, the feasibility of simultaneously measuring gastric emptying, and hepatic and renal glucose uptake and metabolism with dynamic 3D DMI upon administration of deuterated glucose, was investigated. Simulated and measured B₁ ⁺ patterns were in good agreement. The intrasession variability of the natural abundance deuterated water signal in the liver and right kidney, measured in nine healthy volunteers, was 5.6% ± 0.9% and 4.9% ± 0.7%, respectively. Dynamic 3D DMI scans with oral administration of [6,6'-² H₂ ]-glucose showed similar kinetics of deuterated glucose appearance and disappearance in the liver and kidney. The measured gastric emptying half time was 80 ± 10 min, which is in good agreement with scintigraphy measurements. In conclusion, DMI with oral administration of [6,6'-² H₂ ]-glucose enables simultaneous assessment of gastric emptying and liver and kidney glucose uptake and metabolism. When applied in patients with diabetes, this approach may advance our understanding of the interplay between disturbances in liver and kidney glucose uptake and metabolism and gastric emptying, at a detail that cannot be achieved by any other method.

Three-dimensional, 2.5-minute, 7T phosphorus magnetic resonance spectroscopic imaging of the human heart using concentric rings

A three-dimensional (3D), density-weighted, concentric rings trajectory (CRT) magnetic resonance spectroscopic imaging (MRSI) sequence is implemented for cardiac phosphorus (³¹ P)-MRS at 7 T. The point-by-point k-space sampling of traditional phase-encoded chemical shift imaging (CSI) sequences severely restricts the minimum scan time at higher spatial resolutions. Our proposed CRT sequence implements a stack of concentric rings, with a variable number of rings and planes spaced to optimise the density of k-space weighting. This creates flexibility in acquisition time, allowing acquisitions substantially faster than traditional phase-encoded CSI sequences, while retaining high signal-to-noise ratio (SNR). We first characterise the SNR and point-spread function of the CRT sequence in phantoms. We then evaluate it at five different acquisition times and spatial resolutions in the hearts of five healthy participants at 7 T. These different sequence durations are compared with existing published 3D acquisition-weighted CSI sequences with matched acquisition times and spatial resolutions. To minimise the effect of noise on the short acquisitions, low-rank denoising of the spatiotemporal data was also performed after acquisition. The proposed sequence measures 3D localised phosphocreatine to adenosine triphosphate (PCr/ATP) ratios of the human myocardium in 2.5 min, 2.6 times faster than the minimum scan time for acquisition-weighted phase-encoded CSI. Alternatively, in the same scan time, a 1.7-times smaller nominal voxel volume can be achieved. Low-rank denoising reduced the variance of measured PCr/ATP ratios by 11% across all protocols. The faster acquisitions permitted by 7-T CRT ³¹ P-MRSI could make cardiac stress protocols or creatine kinase rate measurements (which involve repeated scans) more tolerable for patients without sacrificing spatial resolution.

Assessment of Bioenergetic Deficits in Patients With Parkinson Disease and Progressive Supranuclear Palsy Using 31P-MRSI

BACKGROUND AND OBJECTIVE

Bioenergetic disturbance, mainly caused by mitochondrial dysfunction, is an established pathophysiologic phenomenon in neurodegenerative movement disorders. The in vivo assessment of brain energy metabolism by 31phosphorus magnetic resonance spectroscopy imaging could provide pathophysiologic insights and serve in the differential diagnosis of parkinsonian disorders. In this study, we investigated such aspects of the underlying pathophysiology in patients with idiopathic Parkinson disease (PwPD) and progressive supranuclear palsy (PwPSP).

METHODS

In total, 30 PwPD, 16 PwPSP, and 25 healthy control subjects (HCs) underwent a clinical examination, structural magnetic resonance imaging, and 31phosphorus magnetic resonance spectroscopy imaging of the forebrain and basal ganglia in a cross-sectional study.

RESULTS

High-energy phosphate metabolites were remarkably decreased in PwPD, particularly in the basal ganglia (-42% compared with HCs and -43% compared with PwPSP, p < 0.0001). This result was not confounded by morphometric brain differences. By contrast, PwPSP had normal levels of high-energy energy metabolites. Thus, the combination of morphometric and metabolic neuroimaging was able to discriminate PwPD from PwPSP with an accuracy of up to 0.93 [95%-CI: 0.91-0.94].

DISCUSSION

Our study shows that mitochondrial dysfunction and bioenergetic depletion contribute to idiopathic Parkinson disease pathophysiology but not to progressive supranuclear palsy. Combined morphometric and metabolic imaging could serve as an accompanying diagnostic biomarker in the neuroimaging-guided differential diagnosis of these parkinsonian disorders.

CLASSIFICATION OF EVIDENCE

This study provides Class III evidence that 31phosphorus magnetic resonance spectroscopy imaging combined with morphometric MRI can differentiate PwPD from PwPSP.

Maternal Cardiac Changes in Women With Obesity and Gestational Diabetes Mellitus

OBJECTIVE

We investigated if women with gestational diabetes mellitus (GDM) in the third trimester of pregnancy exhibit adverse cardiac alterations in myocardial energetics, function, or tissue characteristics.

RESEARCH DESIGN AND METHODS

Thirty-eight healthy, pregnant women and 30 women with GDM were recruited. Participants underwent phosphorus MRS and cardiovascular magnetic resonance for assessment of myocardial energetics (phosphocreatine [PCr] to ATP ratio), tissue characteristics, biventricular volumes and ejection fractions, left ventricular (LV) mass, global longitudinal shortening (GLS), and mitral in-flow E-wave to A-wave ratio.

RESULTS

Participants were matched for age, gestational age, and ethnicity. The following data are reported as mean ± SD. The women with GDM had higher BMI (27 ± 4 vs. 33 ± 5 kg/m2; P = 0.0001) and systolic (115 ± 11 vs. 121 ± 13 mmHg; P = 0.04) and diastolic (72 ± 7 vs. 76 ± 9 mmHg; P = 0.04) blood pressures. There was no difference in N-terminal pro-brain natriuretic peptide concentrations between the groups. The women with GDM had lower myocardial PCr to ATP ratio (2.2 ± 0.3 vs. 1.9 ± 0.4; P < 0.0001), accompanied by lower LV end-diastolic volumes (76 ± 12 vs. 67 ± 11 mL/m2; P = 0.002) and higher LV mass (90 ± 13 vs. 103 ± 18 g; P = 0.001). Although ventricular ejection fractions were similar, the GLS was reduced in women with GDM (-20% ± 3% vs. -18% ± 3%; P = 0.008).

CONCLUSIONS

Despite no prior diagnosis of diabetes, women with obesity and GDM manifest impaired myocardial contractility and higher LV mass, associated with reductions in myocardial energetics in late pregnancy compared with lean women with healthy pregnancy. These findings may aid our understanding of the long-term cardiovascular risks associated with GDM.

Deuterium MRSI of tumor cell death in vivo following oral delivery of 2 H-labeled fumarate

PURPOSE

There is an unmet clinical need for direct and sensitive methods to detect cell death in vivo, especially with regard to monitoring tumor treatment response. We have shown previously that tumor cell death can be detected in vivo from 2 H MRS and MRSI measurements of increased [2,3-2 H2 ]malate production following intravenous injection of [2,3-2 H2 ]fumarate. We show here that cell death can be detected with similar sensitivity following oral administration of the 2 H-labeled fumarate.

METHODS

Mice with subcutaneously implanted EL4 tumors were fasted for 1 h before administration (200 μl) of [2,3-2 H2 ]fumarate (2 g/kg bodyweight) via oral gavage without anesthesia. The animals were then anesthetized, and after 30 min, tumor conversion of [2,3-2 H2 ]fumarate to [2,3-2 H2 ]malate was assessed from a series of 13 2 H spectra acquired over a period of 65 min. The 2 H spectra and 2 H spectroscopic images were acquired using a surface coil before and at 48 h after treatment with a chemotherapeutic drug (etoposide, 67 mg/kg).

RESULTS

The malate/fumarate signal ratio increased from 0.022 ± 0.03 before drug treatment to 0.12 ± 0.04 following treatment (p = 0.023, n = 4). Labeled malate was undetectable in spectroscopic images acquired before treatment and increased in the tumor area following treatment. The increase in the malate/fumarate signal ratio was similar to that observed previously following intravenous administration of labeled fumarate.

CONCLUSION

Orally administered [2,3-2 H2 ]fumarate can be used to detect tumor cell death noninvasively following treatment with a sensitivity that is similar to that obtained with intravenous administration.

Imaging Glioblastoma Response to Radiotherapy Using 2H Magnetic Resonance Spectroscopy Measurements of Fumarate Metabolism

Early detection of tumor cell death in glioblastoma following treatment with chemoradiation has the potential to distinguish between true disease progression and pseudoprogression. Tumor cell death can be detected noninvasively in vivo by imaging the production of [2,3-2H2]malate from [2,3-2H2]fumarate using 2H magnetic resonance (MR) spectroscopic imaging. We show here that 2H MR spectroscopy and spectroscopic imaging measurements of [2,3-2H2]fumarate metabolism can detect tumor cell death in orthotopically implanted glioblastoma models within 48 hours following the completion of chemoradiation. Following the injection of [2,3-2H2]fumarate into tumor-bearing mice, production of [2,3-2H2]malate was measured in a human cell line-derived model and in radiosensitive and radioresistant patient-derived models of glioblastoma that were treated with temozolomide followed by targeted fractionated irradiation. The increase in the [2,3-2H2]malate/[2,3-2H2]fumarate signal ratio posttreatment, which correlated with histologic assessment of cell death, was a more sensitive indicator of treatment response than diffusion-weighted and contrast agent-enhanced 1H MRI measurements, which have been used clinically to detect responses of glioblastoma to chemoradiation. Overall, early detection of glioblastoma cell death using 2H MRI of malate production from fumarate could help improve the clinical evaluation of response to chemoradiation.

SIGNIFICANCE

2H magnetic resonance imaging of labeled fumarate metabolism can detect early evidence of tumor cell death following chemoradiation, meeting a clinical need to reliably detect treatment response in glioblastoma.

Deuterium metabolic imaging and hyperpolarized 13C-MRI of the normal human brain at clinical field strength reveals differential cerebral metabolism

Deuterium metabolic imaging (DMI) and hyperpolarized 13C-pyruvate MRI (13C-HPMRI) are two emerging methods for non-invasive and non-ionizing imaging of tissue metabolism. Imaging cerebral metabolism has potential applications in cancer, neurodegeneration, multiple sclerosis, traumatic brain injury, stroke, and inborn errors of metabolism. Here we directly compare these two non-invasive methods at 3 T for the first time in humans and show how they simultaneously probe both oxidative and non-oxidative metabolism. DMI was undertaken 1-2 h after oral administration of [6,6'-2H2]glucose, and 13C-MRI was performed immediately following intravenous injection of hyperpolarized [1-13C]pyruvate in ten and nine normal volunteers within each arm respectively. DMI was used to generate maps of deuterium-labelled water, glucose, lactate, and glutamate/glutamine (Glx) and the spectral separation demonstrated that DMI is feasible at 3 T. 13C-HPMRI generated maps of hyperpolarized carbon-13 labelled pyruvate, lactate, and bicarbonate. The ratio of 13C-lactate/13C-bicarbonate (mean 3.7 ± 1.2) acquired with 13C-HPMRI was higher than the equivalent 2H-lactate/2H-Glx ratio (mean 0.18 ± 0.09) acquired using DMI. These differences can be explained by the route of administering each probe, the timing of imaging after ingestion or injection, as well as the biological differences in cerebral uptake and cellular physiology between the two molecules. The results demonstrate these two metabolic imaging methods provide different yet complementary readouts of oxidative and reductive metabolism within a clinically feasible timescale. Furthermore, as DMI was undertaken at a clinical field strength within a ten-minute scan time, it demonstrates its potential as a routine clinical tool in the future.

Coexistent Diabetes Is Associated With the Presence of Adverse Phenotypic Features in Patients With Hypertrophic Cardiomyopathy

OBJECTIVE

Type 2 diabetes mellitus (T2DM) is associated with worsened clinical outcomes in hypertrophic cardiomyopathy (HCM) patients. We sought to investigate whether HCM patients with T2DM comorbidity exhibit adverse cardiac alterations in myocardial energetics, function, perfusion, or tissue characteristics.

RESEARCH DESIGN AND METHODS

A total of 55 participants with concomitant HCM and T2DM (HCM-DM) (n = 20) or isolated HCM (n = 20) and healthy volunteers (HV) (n = 15) underwent 31P-MRS and cardiovascular MRI. The HCM groups were matched for HCM phenotype.

RESULTS

Mean ± SD European Society of Cardiology sudden cardiac death risk scores were comparable between the HCM groups (HCM 2.2 ± 1.5%, HCM-DM 1.9 ± 1.2%; P = not significant), and sarcomeric mutations were equally common. HCM-DM patients had the highest median NT-proBNP levels (HV 42 ng/L [interquartile range 35-66], HCM 298 ng/L [157-837], HCM-DM 726 ng/L [213-8,695]; P < 0.0001). Left ventricular (LV) ejection fraction, mass, and wall thickness were similar between the HCM groups. HCM-DM patients displayed a greater degree of fibrosis burden with higher scar percentage and lower global longitudinal strain compared with HCM patients. PCr/ATP (the relative concentrations of phosphocreatine and ATP) was significantly lower in the HCM-DM group than in both HCM and HV (HV 2.17 ± 0.49, HCM 1.93 ± 0.38, HCM-DM 1.54 ± 0.27; P = 0.002). In a similar pattern, stress myocardial blood flow was significantly lower in the HCM-DM group than in both HCM and HV (HV 2.06 ± 0.42 mL/min/g, HCM 1.74 ± 0.44 mL/min/g, HCM-DM 1.39 ± 0.42 mL/min/g; P = 0.002).

CONCLUSIONS

We show for the first time that HCM-DM patients display greater reductions in myocardial energetics, perfusion, and contractile function and higher myocardial scar burden and serum NT-proBNP levels compared with patients with isolated HCM despite similar LV mass and wall thickness and presence of sarcomeric mutations. These adverse phenotypic features may be important components of the adverse clinical manifestation attributable to a combined presence of HCM and T2DM.

Comparison of seven modelling algorithms for γ-aminobutyric acid-edited proton magnetic resonance spectroscopy

Edited MRS sequences are widely used for studying γ-aminobutyric acid (GABA) in the human brain. Several algorithms are available for modelling these data, deriving metabolite concentration estimates through peak fitting or a linear combination of basis spectra. The present study compares seven such algorithms, using data obtained in a large multisite study. GABA-edited (GABA+, TE = 68 ms MEGA-PRESS) data from 222 subjects at 20 sites were processed via a standardised pipeline, before modelling with FSL-MRS, Gannet, AMARES, QUEST, LCModel, Osprey and Tarquin, using standardised vendor-specific basis sets (for GE, Philips and Siemens) where appropriate. After referencing metabolite estimates (to water or creatine), systematic differences in scale were observed between datasets acquired on different vendors' hardware, presenting across algorithms. Scale differences across algorithms were also observed. Using the correlation between metabolite estimates and voxel tissue fraction as a benchmark, most algorithms were found to be similarly effective in detecting differences in GABA+. An interclass correlation across all algorithms showed single-rater consistency for GABA+ estimates of around 0.38, indicating moderate agreement. Upon inclusion of a basis set component explicitly modelling the macromolecule signal underlying the observed 3.0 ppm GABA peaks, single-rater consistency improved to 0.44. Correlation between discrete pairs of algorithms varied, and was concerningly weak in some cases. Our findings highlight the need for consensus on appropriate modelling parameters across different algorithms, and for detailed reporting of the parameters adopted in individual studies to ensure reproducibility and meaningful comparison of outcomes between different studies.

Increased cardiac Pi/PCr in the diabetic heart observed using phosphorus magnetic resonance spectroscopy at 7T

Phosphorus magnetic resonance spectroscopy (³¹P-MRS) has previously demonstrated decreased energy reserves in the form of phosphocreatine to adenosine-tri-phosphate ratio (PCr/ATP) in the hearts of patients with type 2 diabetes (T2DM). Recent ³¹P-MRS techniques using 7T systems, e.g. long mixing time stimulated echo acquisition mode (STEAM), allow deeper insight into cardiac metabolism through assessment of inorganic phosphate (Pi) content and myocardial pH, which play pivotal roles in energy production in the heart. Therefore, we aimed to further explore the cardiac metabolic phenotype in T2DM using STEAM at 7T. Seventeen patients with T2DM and twenty-three healthy controls were recruited and their cardiac PCr/ATP, Pi/PCr and pH were assessed at 7T. Diastolic function of all patients with T2DM was assessed using echocardiography to investigate the relationship between diastolic dysfunction and cardiac metabolism. Mirroring the decreased PCr/ATP (1.70±0.31 vs. 2.07±0.39; p<0.01), the cardiac Pi/PCr was increased (0.13±0.07 vs. 0.10±0.03; p = 0.02) in T2DM patients in comparison to healthy controls. Myocardial pH was not significantly different between the groups (7.14±0.12 vs. 7.10±0.12; p = 0.31). There was a negative correlation between PCr/ATP and diastolic function (R² = 0.33; p = 0.02) in T2DM. No correlation was observed between diastolic function and Pi/PCr and (R² = 0.16; p = 0.21). In addition, we did not observe any correlation between cardiac PCr/ATP and Pi/PCr (p = 0.19). Using STEAM ³¹P-MRS at 7T we have for the first time explored Pi/PCr in the diabetic human heart and found it increased when compared to healthy controls. The lack of correlation between measured PCr/ATP and Pi/PCr suggests that independent mechanisms might contribute to these perturbations.

Hyperpolarized 13 C MRI Reveals Large Changes in Pyruvate Metabolism During Digestion in Snakes

Purpose

Methods

Results

Conclusion

The NADPARK study: A randomized phase I trial of nicotinamide riboside supplementation in Parkinson's disease

We conducted a double-blinded phase I clinical trial to establish whether nicotinamide adenine dinucleotide (NAD) replenishment therapy, via oral intake of nicotinamide riboside (NR), is safe, augments cerebral NAD levels, and impacts cerebral metabolism in Parkinson's disease (PD). Thirty newly diagnosed, treatment-naive patients received 1,000 mg NR or placebo for 30 days. NR treatment was well tolerated and led to a significant, but variable, increase in cerebral NAD levels-measured by ³¹phosphorous magnetic resonance spectroscopy-and related metabolites in the cerebrospinal fluid. NR recipients showing increased brain NAD levels exhibited altered cerebral metabolism, measured by ¹⁸fluoro-deoxyglucose positron emission tomography, and this was associated with mild clinical improvement. NR augmented the NAD metabolome and induced transcriptional upregulation of processes related to mitochondrial, lysosomal, and proteasomal function in blood cells and/or skeletal muscle. Furthermore, NR decreased the levels of inflammatory cytokines in serum and cerebrospinal fluid. Our findings nominate NR as a potential neuroprotective therapy for PD, warranting further investigation in larger trials.

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.

Myocardial Energy Response to Glyceryl Trinitrate: Physiology Revisited

Objective: Although intravenous nitrates are commonly used in clinical medicine, they have been shown to increase myocardial oxygen consumption and inhibit complex IV of the electron transport chain. As such we sought to measure whether myocardial energetics were impaired during glyceryl trinitrate (GTN) infusion. Methods: 10 healthy volunteers underwent cardiac magnetic resonance imaging to assess cardiac function and 31phosphorus magnetic resonance spectroscopy to measure Phosphocreatine/ATP (PCr/ATP) ratio and creatine kinase forward rate constant (CK kf ) before and during an intravenous infusion of GTN. Results: During GTN infusion, mean arterial pressure (78 ± 7 vs. 65 ± 6 mmHg, p < 0.001), left ventricular (LV) stroke work (7,708 ± 2,782 vs. 6,071 ± 2,660 ml mmHg, p < 0.001), and rate pressure product (7,214 ± 1,051 vs. 6,929 ± 976 mmHg bpm, p = 0.06) all fell. LV ejection fraction increased (61 ± 3 vs. 66 ± 4%, p < 0.001), with cardiac output remaining constant (6.2 ± 1.5 vs. 6.5 ± 1.4 l/min, p = 0.37). Myocardial PCr/ATP fell during GTN infusion (2.17 ± 0.2 vs. 1.99 ± 0.22, p = 0.03) with an increase in both CK kf (0.16 ± 0.07 vs. 0.25 ± 0.1 s-1, p = 0.006) and CK flux (1.8 ± 0.8 vs. 2.6 ± 1.1 μmol/g/s, p = 0.03). Conclusion: During GTN infusion, despite reduced LV stroke work and maintained cardiac output, there was a 44% increase in myocardial ATP delivery through CK. As PCr/ATP fell, this increase in ATP demand coincided with GTN-induced impairment of mitochondrial oxidative phosphorylation. Overall, this suggests that while GTN reduces cardiac work, it does so at the expense of increasing ATP demand beyond the capacity to increase ATP production.

Deuterium MRS of early treatment-induced changes in tumour lactate in vitro

Elevated production of lactate is a key characteristic of aberrant tumour cell metabolism and can be non-invasively measured as an early marker of tumour response using deuterium (² H) MRS. Following treatment, changes in the ² H-labelled lactate signal could identify tumour cell death or impaired metabolic function, which precede morphological changes conventionally used to assess tumour response. In this work, the association between apoptotic cell death, extracellular lactate concentration, and early treatment-induced changes in the ² H-labelled lactate signal was established in an in vitro tumour model. Experiments were conducted at 7 T on acute myeloid leukaemia (AML) cells, which had been treated with 10 μg/mL of the chemotherapeutic agent cisplatin. At 24 and 48 h after cisplatin treatment the cells were supplied with 20 mM of [6,6'-² H₂ ]glucose and scanned over 2 h using a two-dimensional ² H MR spectroscopic imaging sequence. The resulting signals from ² H-labelled glucose, lactate, and water were quantified using a spectral fitting algorithm implemented on the Oxford Spectroscopy Analysis MATLAB toolbox. After scanning, the cells were processed for histological stains (terminal deoxynucleotidyl transferase UTP nick end labelling and haematoxylin and eosin) to assess apoptotic area fraction and cell morphology respectively, while a colorimetric assay was used to measure extracellular lactate concentrations in the supernatant. Significantly lower levels of ² H-labelled lactate were observed in the 48 h treated cells compared with the untreated and 24 h treated cells, and these changes were significantly correlated with an increase in apoptotic fraction and a decrease in extracellular lactate. By establishing the biological processes associated with treatment-induced changes in the ² H-labelled lactate signal, these findings suggest that ² H MRS of lactate may be valuable in evaluating early tumour response.

Energetic Basis for Exercise-Induced Pulmonary Congestion in Heart Failure With Preserved Ejection Fraction

BACKGROUND

Transient pulmonary congestion during exercise is emerging as an important determinant of reduced exercise capacity in heart failure with preserved ejection fraction (HFpEF). We sought to determine whether an abnormal cardiac energetic state underpins this process.

METHODS

We recruited patients across the spectrum of diastolic dysfunction and HFpEF (controls, n=11; type 2 diabetes, n=9; HFpEF, n=14; and severe diastolic dysfunction attributable to cardiac amyloidosis, n=9). Cardiac energetics were measured using phosphorus spectroscopy to define the myocardial phosphocreatine to ATP ratio. Cardiac function was assessed by cardiovascular magnetic resonance cine imaging and echocardiography and lung water using magnetic resonance proton density mapping. Studies were performed at rest and during submaximal exercise using a magnetic resonance imaging ergometer.

RESULTS

Paralleling the stepwise decline in diastolic function across the groups (E/e' ratio; P<0.001) was an increase in NT-proBNP (N-terminal pro-brain natriuretic peptide; P<0.001) and a reduction in phosphocreatine/ATP ratio (control, 2.15 [2.09, 2.29]; type 2 diabetes, 1.71 [1.61, 1.91]; HFpEF, 1.66 [1.44, 1.89]; cardiac amyloidosis, 1.30 [1.16, 1.53]; P<0.001). During 20-W exercise, lower left ventricular diastolic filling rates (r=0.58; P<0.001), lower left ventricular diastolic reserve (r=0.55; P<0.001), left atrial dilatation (r=-0.52; P<0.001), lower right ventricular contractile reserve (right ventricular ejection fraction change, r=0.57; P<0.001), and right atrial dilation (r=-0.71; P<0.001) were all linked to lower phosphocreatine/ATP ratio. Along with these changes, pulmonary proton density mapping revealed transient pulmonary congestion in patients with HFpEF (+4.4% [0.5, 6.4]; P=0.002) and cardiac amyloidosis (+6.4% [3.3, 10.0]; P=0.004), which was not seen in healthy controls (-0.1% [-1.9, 2.1]; P=0.89) or type 2 diabetes without HFpEF (+0.8% [-1.7, 1.9]; P=0.82). The development of exercise-induced pulmonary congestion was associated with lower phosphocreatine/ATP ratio (r=-0.43; P=0.004).

CONCLUSIONS

A gradient of myocardial energetic deficit exists across the spectrum of HFpEF. Even at low workload, this energetic deficit is related to markedly abnormal exercise responses in all 4 cardiac chambers, which is associated with detectable pulmonary congestion. The findings support an energetic basis for transient pulmonary congestion in HFpEF.

Obesity modifies the energetic phenotype of dilated cardiomyopathy

AIMS

We sought to determine if myocardial energetics could distinguish obesity cardiomyopathy as a distinct entity from dilated cardiomyopathy.

METHODS AND RESULTS

Sixteen normal weight participants with dilated cardiomyopathy (DCMNW), and 27 with DCM and obesity (DCMOB), were compared to 26 normal weight controls (CTLNW). All underwent cardiac magnetic resonance imaging and 31P spectroscopy to assess function and energetics. Nineteen DCMOB underwent repeat assessment after a dietary weight loss intervention. Adenosine triphosphate (ATP) delivery through creatine kinase (CK flux) was 55% lower in DCMNW than in CTLNW (P = 0.004), correlating with left ventricular ejection fraction (LVEF, r = 0.4, P = 0.015). In contrast, despite similar LVEF (DCMOB 41 ± 7%, DCMNW 38 ± 6%, P = 0.14), CK flux was two-fold higher in DCMOB (P < 0.001), due to higher rate through CK [median kf 0.21 (0.14) vs. 0.11 (0.12) s-1, P = 0.002]. During increased workload, the CTLNW heart increased CK flux by 97% (P < 0.001). In contrast, CK flux was unchanged in DCMNW and fell in DCMOB (by >50%, P < 0.001). Intentional weight loss was associated with positive left ventricular remodelling, with reduced left ventricular end-diastolic volume (by 8%, P < 0.001) and a change in LVEF (40 ± 9% vs. 45 ± 10%, P = 0.002). This occurred alongside a fall in ATP delivery rate with weight loss (by 7%, P = 0.049).

CONCLUSIONS

In normal weight, DCM is associated with reduced resting ATP delivery. In obese DCM, ATP demand through CK is greater, suggesting reduced efficiency of energy utilization. Dietary weight loss is associated with significant improvement in myocardial contractility, and a fall in ATP delivery, suggesting improved metabolic efficiency. This highlights distinct energetic pathways in obesity cardiomyopathy, which are both different from dilated cardiomyopathy, and may be reversible with weight loss.

Diabetic mitochondria are resistant to palmitoyl CoA inhibition of respiration, which is detrimental during ischemia

The bioactive lipid intermediate palmitoyl CoA (PCoA) can inhibit mitochondrial ADP/ATP transport, though the physiological relevance of this regulation remains unclear. We questioned whether myocardial ischemia provides a pathological setting in which PCoA regulation of ADP/ATP transport would be beneficial, and secondly, whether the chronically elevated lipid content within the diabetic heart could make mitochondria less sensitive to the effects of PCoA. PCoA acutely decreased ADP‐stimulated state 3 respiration and increased the apparent K_m for ADP twofold. The half maximal inhibitory concentration (IC₅₀) of PCoA in control mitochondria was 22 µM. This inhibitory effect of PCoA on respiration was blunted in diabetic mitochondria, with no significant difference in the K_m for ADP in the presence of PCoA, and an increase in the IC₅₀ to 32 µM PCoA. The competitive inhibition by PCoA was localised to the phosphorylation apparatus, particularly the ADP/ATP carrier (AAC). During ischemia, the AAC imports ATP into the mitochondria, where it is hydrolysed by reversal of the ATP synthase, regenerating the membrane potential. Addition of PCoA dose‐dependently prevented this wasteful ATP hydrolysis for membrane repolarisation during ischemia, however, this beneficial effect was blunted in diabetic mitochondria. Finally, using ³¹P‐magnetic resonance spectroscopy we demonstrated that diabetic hearts lose ATP more rapidly during ischemia, with a threefold higher ATP decay rate compared with control hearts. In conclusion, PCoA plays a role in protecting mitochondrial energetics during ischemia, by preventing wasteful ATP hydrolysis. However, this beneficial effect is blunted in diabetes, contributing to the impaired energy metabolism seen during myocardial ischemia in the diabetic heart.

Noninvasive assessment of steatosis and viability of cold-stored human liver grafts by MRI

Purpose

A shortage of suitable donor livers is driving increased use of higher risk livers for transplantation. However, current biomarkers are not sensitive and specific enough to predict posttransplant liver function. This is limiting the expansion of the donor pool. Therefore, better noninvasive tests are required to determine which livers will function following implantation and hence can be safely transplanted. This study assesses the temperature sensitivity of proton density fat fraction and relaxometry parameters and examines their potential for assessment of liver function ex vivo.

Methods

Six ex vivo human livers were scanned during static cold storage following normothermic machine perfusion. Proton density fat fraction, T₁, T₂, and T2* were measured repeatedly during cooling on ice. Temperature corrections were derived from these measurements for the parameters that showed significant variation with temperature.

Results

Strong linear temperature sensitivities were observed for proton density fat fraction (R² = 0.61, P < .001) and T₁ (R² = 0.78, P < .001). Temperature correction according to a linear model reduced the coefficient of repeatability in these measurements by 41% and 36%, respectively. No temperature dependence was observed in T₂ or T2* measurements. Comparing livers deemed functional and nonfunctional during normothermic machine perfusion by hemodynamic and biochemical criteria, T₁ differed significantly: 516 ± 50 ms for functional versus 679 ± 60 ms for non-functional, P = .02.

Conclusion

Temperature correction is essential for robust measurement of proton density fat fraction and T₁ in cold-stored human livers. These parameters may provide a noninvasive measure of viability for transplantation.

Water-suppression cycling 3-T cardiac 1 H-MRS detects altered creatine and choline in patients with aortic or mitral stenosis

Cardiac proton spectroscopy (1 H-MRS) is widely used to quantify lipids. Other metabolites (e.g. creatine and choline) are clinically relevant but more challenging to quantify because of their low concentrations (approximately 10 mmol/L) and because of cardiac motion. To quantify cardiac creatine and choline, we added water-suppression cycling (WSC) to two single-voxel spectroscopy sequences (STEAM and PRESS). WSC introduces controlled residual water signals that alternate between positive and negative phases from transient to transient, enabling robust phase and frequency correction. Moreover, a particular weighted sum of transients eliminates residual water signals without baseline distortion. We compared WSC and the vendor's standard 'WET' water suppression in phantoms. Next, we tested repeatability in 10 volunteers (seven males, three females; age 29.3 ± 4.0 years; body mass index [BMI] 23.7 ± 4.1 kg/m2 ). Fat fraction, creatine concentration and choline concentration when quantified by STEAM-WET were 0.30% ± 0.11%, 29.6 ± 7.0 μmol/g and 7.9 ± 6.7 μmol/g, respectively; and when quantified by PRESS-WSC they were 0.30% ± 0.15%, 31.5 ± 3.1 μmol/g and 8.3 ± 4.4 μmol/g, respectively. Compared with STEAM-WET, PRESS-WSC gave spectra whose fitting quality expressed by Cramér-Rao lower bounds improved by 26% for creatine and 32% for choline. Repeatability of metabolite concentration measurements improved by 72% for creatine and 40% for choline. We also compared STEAM-WET and PRESS-WSC in 13 patients with severe symptomatic aortic or mitral stenosis indicated for valve replacement surgery (10 males, three females; age 75.9 ± 6.3 years; BMI 27.4 ± 4.3 kg/m2 ). Spectra were of analysable quality in eight patients for STEAM-WET, and in nine for PRESS-WSC. We observed comparable lipid concentrations with those in healthy volunteers, significantly reduced creatine concentrations, and a trend towards decreased choline concentrations. We conclude that PRESS-WSC offers improved performance and reproducibility for the quantification of cardiac lipids, creatine and choline concentrations in healthy volunteers at 3 T. It also offers improved performance compared with STEAM-WET for detecting altered creatine and choline concentrations in patients with valve disease.

Liver cT1 decreases following direct-acting antiviral therapy in patients with chronic hepatitis C virus

PURPOSE

Direct-acting antiviral therapies (DAAs) for treatment of chronic hepatitis C virus (HCV) have excellent rates of viral eradication, but their effect on regression of liver fibrosis is unclear. The primary aim was to use magnetic resonance imaging (MRI) and spectroscopy (MRS) to evaluate changes in liver fibrosis, liver fat and liver iron content (LIC) in patients with chronic HCV following treatment with DAAs.

METHODS

In this prospective study, 15 patients with chronic HCV due to start treatment with DAAs and with transient elastography (TE) > 8 kPa were recruited consecutively. Patients underwent MRI and MRS at baseline (before treatment), and at 24 weeks and 48 weeks after the end of treatment (EoT) for the measurement of liver cT1 (fibroinflammation), liver fat and T2* (LIC).

RESULTS

All patients achieved a sustained virological response. Liver cT1 showed significant decreases from baseline to 24 weeks post EoT (876 vs 806 ms, p = 0.002, n = 15), baseline to 48 weeks post EoT (876 vs 788 ms, p = 0.0002, n = 13) and 24 weeks post EoT to 48 weeks post EoT (806 vs 788 ms, p = 0.016, n = 13). Between baseline and 48 weeks EoT significant reduction in liver fat (5.17% vs 2.65%, p = 0.027) and an increase in reported LIC (0.913 vs 0.950 mg/g, p = 0.021) was observed.

CONCLUSION

Liver cT1 decreases in patients with chronic HCV undergoing successful DAA treatment. The relatively fast reduction in cT1 suggests a reduction in inflammation rather than regression of fibrosis.

Investigating the effect of trigger delay on cardiac 31P MRS signals

The heart’s geometry and its metabolic activity vary over the cardiac cycle. The effect of these fluctuations on phosphorus (³¹P) magnetic resonance spectroscopy (MRS) data quality and metabolite ratios was investigated. 12 healthy volunteers were measured using a 7 T MR scanner and a cardiac ³¹P-¹H loop coil. ³¹P chemical shift imaging data were acquired untriggered and at four different times during the cardiac cycle using acoustic triggering. Signals of adenosine-triphosphate (ATP), phosphocreatine (PCr), inorganic phosphate (Pi) and 2,3-diphosphoglycerate (2,3-DPG) and their fit quality as Cramér-Rao lower bounds (CRLB) were quantified including corrections for contamination by ³¹P signals from blood, flip angle, saturation and total acquisition time. The myocardial filling factor was estimated from cine short axis views. The corrected signals of PCr and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma$$\end{document}γ-ATP were higher during end-systole and lower during diastasis than in untriggered acquisitions (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P<0.05$$\end{document}P<0.05). Signal intensities of untriggered scans were between those with triggering to end-systole and diastasis. Fit quality of PCr and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma$$\end{document}γ-ATP peaks was best during end-systole when blood contamination of ATP and Pi signals was lowest. While metabolite ratios and pH remained stable over the cardiac cycle, signal amplitudes correlated strongly with myocardial voxel filling. Triggering of cardiac ³¹P MRS acquisitions improves signal amplitudes and fit quality if the trigger delay is set to end-systole. We conclude that triggering to end-systole is superior to triggering to diastasis.

Schulte R, Laustsen C. Comprehensive Literature Review of Hyperpolarized Carbon-13 MRI: The Road to Clinical Application

This review provides a comprehensive assessment of the development of hyperpolarized (HP) carbon-13 metabolic MRI from the early days to the present with a focus on clinical applications. The status and upcoming challenges of translating HP carbon-13 into clinical application are reviewed, along with the complexity, technical advancements, and future directions. The road to clinical application is discussed regarding clinical needs and technological advancements, highlighting the most recent successes of metabolic imaging with hyperpolarized carbon-13 MRI. Given the current state of hyperpolarized carbon-13 MRI, the conclusion of this review is that the workflow for hyperpolarized carbon-13 MRI is the limiting factor.

MR spectroscopy in pediatric neuroradiology

Magnetic resonance spectroscopy (MRS), being able to identify and measure some brain components (metabolites) in pathologic lesions and in normal-appearing tissue, offers a valuable additional diagnostic tool to assess several pediatric neurological diseases. In this review we will illustrate the basic principles and clinical applications of brain proton (H¹; hydrogen) MRS (H¹MRS), by now the only MRS method widely available in clinical practice. Performing H¹MRS in the brain is inherently less complicated than in other tissues (e.g., liver, muscle), in which spectra are heavily affected by magnetic field inhomogeneities, respiration artifacts, and dominating signals from the surrounding adipose tissues. H¹MRS in pediatric neuroradiology has some advantages over acquisitions in adults (lack of motion due to children sedation and lack of brain iron deposition allow optimal results), but it requires a deep knowledge of pediatric pathologies and familiarity with the developmental changes in spectral patterns, particularly occurring in the first two years of life. Examples from our database, obtained mainly from a 1.5 Tesla clinical scanner in a time span of 15 years, will demonstrate the efficacy of H¹MRS in the diagnosis of a wide range of selected pediatric pathologies, like brain tumors, infections, neonatal hypoxic-ischemic encephalopathy, metabolic and white matter disorders.

Left atrial 4D flow cardiovascular magnetic resonance: a reproducibility study in sinus rhythm and atrial fibrillation

BACKGROUND

Four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) allows sophisticated quantification of left atrial (LA) blood flow, and could yield novel biomarkers of propensity for intra-cardiac thrombus formation and embolic stroke. As reproducibility is critically important to diagnostic performance, we systematically investigated technical and temporal variation of LA 4D flow in atrial fibrillation (AF) and sinus rhythm (SR).

METHODS

Eighty-six subjects (SR, n = 64; AF, n = 22) with wide-ranging stroke risk (CHA2DS2VASc 0-6) underwent LA 4D flow assessment of peak and mean velocity, vorticity, vortex volume, and stasis. Eighty-five (99%) underwent a second acquisition within the same session, and 74 (86%) also returned at 30 (27-35) days for an interval scan. We assessed variability attributable to manual contouring (intra- and inter-observer), and subject repositioning and reacquisition of data, both within the same session (same-day scan-rescan), and over time (interval scan). Within-subject coefficients of variation (CV) and bootstrapped 95% CIs were calculated and compared.

RESULTS

Same-day scan-rescan CVs were 6% for peak velocity, 5% for mean velocity, 7% for vorticity, 9% for vortex volume, and 10% for stasis, and were similar between SR and AF subjects (all p > 0.05). Interval-scan variability was similar to same-day scan-rescan variability for peak velocity, vorticity, and vortex volume (all p > 0.05), and higher for stasis and mean velocity (interval scan CVs of 14% and 8%, respectively, both p < 0.05). Longitudinal changes in heart rate and blood pressure at the interval scan in the same subjects were associated with significantly higher variability for LA stasis (p = 0.024), but not for the remaining flow parameters (all p > 0.05). SR subjects showed significantly greater interval-scan variability than AF patients for mean velocity, vortex volume, and stasis (all p < 0.05), but not peak velocity or vorticity (both p > 0.05).

CONCLUSIONS

LA peak velocity and vorticity are the most reproducible and temporally stable novel LA 4D flow biomarkers, and are robust to changes in heart rate, blood pressure, and differences in heart rhythm.

Rapid, B 1 -insensitive, dual-band quasi-adiabatic saturation transfer with optimal control for complete quantification of myocardial ATP flux

PURPOSE

Phosphorus saturation-transfer experiments can quantify metabolic fluxes noninvasively. Typically, the forward flux through the creatine kinase reaction is investigated by observing the decrease in phosphocreatine (PCr) after saturation of γ-ATP. The quantification of total ATP utilization is currently underexplored, as it requires simultaneous saturation of inorganic phosphate ( P i ) and PCr. This is challenging, as currently available saturation pulses reduce the already-low γ-ATP signal present.

METHODS

Using a hybrid optimal-control and Shinnar-Le Roux method, a quasi-adiabatic RF pulse was designed for the dual saturation of PCr and P i to enable determination of total ATP utilization. The pulses were evaluated in Bloch equation simulations, compared with a conventional hard-cosine DANTE saturation sequence, before being applied to perfused rat hearts at 11.7 T.

RESULTS

The quasi-adiabatic pulse was insensitive to a >2.5-fold variation in B 1 , producing equivalent saturation with a 53% reduction in delivered pulse power and a 33-fold reduction in spillover at the minimum effective B 1 . This enabled the complete quantification of the synthesis and degradation fluxes for ATP in 30-45 minutes in the perfused rat heart. While the net synthesis flux (4.24 ± 0.8 mM/s, SEM) was not significantly different from degradation flux (6.88 ± 2 mM/s, P = .06) and both measures are consistent with prior work, nonlinear error analysis highlights uncertainties in the P_i -to-ATP measurement that may explain a trend suggesting a possible imbalance.

CONCLUSIONS

This work demonstrates a novel quasi-adiabatic dual-saturation RF pulse with significantly improved performance that can be used to measure ATP turnover in the heart in vivo.

Nicotinic acid receptor agonists impair myocardial contractility by energy starvation

Nicotinic acid receptor agonists have previously been shown to cause acute reductions in cardiac contractility. We sought to uncover the changes in cardiac metabolism underlying these alterations in function. In nine humans, we recorded cardiac energetics and function before and after a single oral dose of nicotinic acid using cardiac MRI to demonstrate contractile function and Phosphorus-31 (31 P) magnetic resonance spectroscopy to demonstrate myocardial energetics. Nicotinic Acid 400 mg lowered ejection fraction by 4% (64 ± 8% to 60 ± 7%, P = .03), and was accompanied by a fall in phosphocreatine/ATP ratio by 0.4 (2.2 ± 0.4 to 1.8 ± 0.1, P = .04). In four groups of eight Wistar rats, we used pyruvate dehydrogenase (PDH) flux studies to demonstrate changes in carbohydrate metabolism induced by the nicotinic acid receptor agonist, Acipimox, using hyperpolarized Carbon-13 (13 C) magnetic resonance spectroscopy. In rats which had been starved overnight, Acipimox caused a fall in ejection fraction by 7.8% (67.5 ± 8.9 to 60 ± 3.1, P = .03) and a nearly threefold rise in flux through PDH (from 0.182 ± 0.114 to 0.486 ± 0.139, P = .002), though this rise did not match pyruvate dehydrogenase flux observed in rats fed carbohydrate rich chow (0.726 ± 0.201). In fed rats, Acipimox decreased pyruvate dehydrogenase flux (to 0.512 ± 0.13, P = .04). Concentration of plasma insulin fell by two-thirds in fed rats administered Acipimox (from 1695 ± 891 ng/L to 550 ± 222 ng/L, P = .005) in spite of glucose concentrations remaining the same. In conclusion, we demonstrate that nicotinic acid receptor agonists impair cardiac contractility associated with a decline in cardiac energetics and show that the mechanism is likely a combination of reduced fatty acid availability and a failure to upregulate carbohydrate metabolism, essentially starving the heart of fuel.

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.

Rescue of myocardial energetic dysfunction in diabetes through the correction of mitochondrial hyperacetylation by honokiol

Cardiac energetic dysfunction has been reported in patients with type 2 diabetes (T2D) and is an independent predictor of mortality. Identification of the mechanisms driving mitochondrial dysfunction, and therapeutic strategies to rescue these modifications, will improve myocardial energetics in T2D. We demonstrate using ³¹P-magnetic resonance spectroscopy (³¹P-MRS) that decreased cardiac ATP and phosphocreatine (PCr) concentrations occurred before contractile dysfunction or a reduction in PCr/ATP ratio in T2D. Real-time mitochondrial ATP synthesis rates and state 3 respiration rates were similarly depressed in T2D, implicating dysfunctional mitochondrial energy production. Driving this energetic dysfunction in T2D was an increase in mitochondrial protein acetylation, and increased ex vivo acetylation was shown to proportionally decrease mitochondrial respiration rates. Treating T2D rats in vivo with the mitochondrial deacetylase SIRT3 activator honokiol reversed the hyperacetylation of mitochondrial proteins and restored mitochondrial respiration rates to control levels. Using ¹³C-hyperpolarized MRS, respiration with different substrates, and enzyme assays, we localized this improvement to increased glutamate dehydrogenase activity. Finally, honokiol treatment increased ATP and PCr concentrations and increased total ATP synthesis flux in the T2D heart. In conclusion, hyperacetylation drives energetic dysfunction in T2D, and reversing acetylation with the SIRT3 activator honokiol rescued myocardial and mitochondrial energetics in T2D.

Pharmacologically targeting mitochondrial acetylation provides a mechanism to rescue impaired myocardial energy generation in type 2 diabetes.

Osprey: Open-source processing, reconstruction & estimation of magnetic resonance spectroscopy data

BACKGROUND

Processing and quantitative analysis of magnetic resonance spectroscopy (MRS) data are far from standardized and require interfacing with third-party software. Here, we present Osprey, a fully integrated open-source data analysis pipeline for MRS data, with seamless integration of pre-processing, linear-combination modelling, quantification, and data visualization.

NEW METHOD

Osprey loads multiple common MRS data formats, performs phased-array coil combination, frequency-and phase-correction of individual transients, signal averaging and Fourier transformation. Linear combination modelling of the processed spectrum is carried out using simulated basis sets and a spline baseline. The MRS voxel is coregistered to an anatomical image, which is segmented for tissue correction and quantification is performed based upon modelling parameters and tissue segmentation. The results of each analysis step are visualized in the Osprey GUI. The analysis pipeline is demonstrated in 12 PRESS, 11 MEGA-PRESS, and 8 HERMES datasets acquired in healthy subjects.

RESULTS

Osprey successfully loads, processes, models, and quantifies MRS data acquired with a variety of conventional and spectral editing techniques.

COMPARISON WITH EXISTING METHOD(S)

Osprey is the first MRS software to combine uniform pre-processing, linear-combination modelling, tissue correction and quantification into a coherent ecosystem. Compared to existing compiled, often closed-source modelling software, Osprey's open-source code philosophy allows researchers to integrate state-of-the-art data processing and modelling routines, and potentially converge towards standardization of analysis.

CONCLUSIONS

Osprey combines robust, peer-reviewed data processing methods into a modular workflow that is easily augmented by community developers, allowing the rapid implementation of new methods.

Use of cardiac magnetic resonance to detect changes in metabolism in heart failure

The heart has a massive adenosine triphosphate (ATP) requirement, produced from the oxidation of metabolic substrates such as fat and glucose. Magnetic resonance spectroscopy offers a unique opportunity to probe this biochemistry: ³¹Phosphorus spectroscopy can demonstrate the production of ATP and quantify levels of the transport molecule phosphocreatine while 13Carbon spectroscopy can demonstrate the metabolic fates of glucose in real time. These techniques allow the metabolic deficits in heart failure to be interrogated and can be a potential future clinical tool.