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.

Whole-brain deuterium metabolic imaging via concentric ring trajectory readout enables assessment of regional variations in neuronal glucose metabolism.

Deuterium metabolic imaging (DMI) is an emerging magnetic resonance technique, for non-invasive mapping of human brain glucose metabolism following oral or intravenous administration of deuterium-labeled glucose. Regional differences in glucose metabolism can be observed in various brain pathologies, such as Alzheimer's disease, cancer, epilepsy or schizophrenia, but the achievable spatial resolution of conventional phase-encoded DMI methods is limited due to prolonged acquisition times rendering submilliliter isotropic spatial resolution for dynamic whole brain DMI not feasible. The purpose of this study was to implement non-Cartesian spatial-spectral sampling schemes for whole-brain 2H FID-MR Spectroscopic Imaging to assess time-resolved metabolic maps with sufficient spatial resolution to reliably detect metabolic differences between healthy gray and white matter regions. Results were compared with lower-resolution DMI maps, conventionally acquired within the same session. Six healthy volunteers (4 m/2 f) were scanned for ~90 min after administration of 0.8 g/kg oral [6,6']-2H glucose. Time-resolved whole brain 2H FID-DMI maps of glucose (Glc) and glutamate + glutamine (Glx) were acquired with 0.75 and 2 mL isotropic spatial resolution using density-weighted concentric ring trajectory (CRT) and conventional phase encoding (PE) readout, respectively, at 7 T. To minimize the effect of decreased signal-to-noise ratios associated with smaller voxels, low-rank denoising of the spatiotemporal data was performed during reconstruction. Sixty-three minutes after oral tracer uptake three-dimensional (3D) CRT-DMI maps featured 19% higher (p = .006) deuterium-labeled Glc concentrations in GM (1.98 ± 0.43 mM) compared with WM (1.66 ± 0.36 mM) dominated regions, across all volunteers. Similarly, 48% higher (p = .01) 2H-Glx concentrations were observed in GM (2.21 ± 0.44 mM) compared with WM (1.49 ± 0.20 mM). Low-resolution PE-DMI maps acquired 70 min after tracer uptake featured smaller regional differences between GM- and WM-dominated areas for 2H-Glc concentrations with 2.00 ± 0.35 mM and 1.71 ± 0.31 mM, respectively (+16%; p = .045), while no regional differences were observed for 2H-Glx concentrations. In this study, we successfully implemented 3D FID-MRSI with fast CRT encoding for dynamic whole-brain DMI at 7 T with 2.5-fold increased spatial resolution compared with conventional whole-brain phase encoded (PE) DMI to visualize regional metabolic differences. The faster metabolic activity represented by 48% higher Glx concentrations was observed in GM- compared with WM-dominated regions, which could not be reproduced using whole-brain DMI with the low spatial resolution protocol. Improved assessment of regional pathologic alterations using a fully non-invasive imaging method is of high clinical relevance and could push DMI one step toward clinical applications.

Universal dynamic fitting of magnetic resonance spectroscopy.

PURPOSE: Dynamic (2D) MRS is a collection of techniques where acquisitions of spectra are repeated under varying experimental or physiological conditions. Dynamic MRS comprises a rich set of contrasts, including diffusion-weighted, relaxation-weighted, functional, edited, or hyperpolarized spectroscopy, leading to quantitative insights into multiple physiological or microstructural processes. Conventional approaches to dynamic MRS analysis ignore the shared information between spectra, and instead proceed by independently fitting noisy individual spectra before modeling temporal changes in the parameters. Here, we propose a universal dynamic MRS toolbox which allows simultaneous fitting of dynamic spectra of arbitrary type. METHODS: A simple user-interface allows information to be shared and precisely modeled across spectra to make inferences on both spectral and dynamic processes. We demonstrate and thoroughly evaluate our approach in three types of dynamic MRS techniques. Simulations of functional and edited MRS are used to demonstrate the advantages of dynamic fitting. RESULTS: Analysis of synthetic functional 1H-MRS data shows a marked decrease in parameter uncertainty as predicted by prior work. Analysis with our tool replicates the results of two previously published studies using the original in vivo functional and diffusion-weighted data. Finally, joint spectral fitting with diffusion orientation models is demonstrated in synthetic data. CONCLUSION: A toolbox for generalized and universal fitting of dynamic, interrelated MR spectra has been released and validated. The toolbox is shared as a fully open-source software with comprehensive documentation, example data, and tutorials.

Diffusion-weighted MR spectroscopy: Consensus, recommendations, and resources from acquisition to modeling.

Brain cell structure and function reflect neurodevelopment, plasticity, and aging; and changes can help flag pathological processes such as neurodegeneration and neuroinflammation. Accurate and quantitative methods to noninvasively disentangle cellular structural features are needed and are a substantial focus of brain research. Diffusion-weighted MRS (dMRS) gives access to diffusion properties of endogenous intracellular brain metabolites that are preferentially located inside specific brain cell populations. Despite its great potential, dMRS remains a challenging technique on all levels: from the data acquisition to the analysis, quantification, modeling, and interpretation of results. These challenges were the motivation behind the organization of the Lorentz Center workshop on "Best Practices & Tools for Diffusion MR Spectroscopy" held in Leiden, the Netherlands, in September 2021. During the workshop, the dMRS community established a set of recommendations to execute robust dMRS studies. This paper provides a description of the steps needed for acquiring, processing, fitting, and modeling dMRS data, and provides links to useful resources.

Cerebellar GABA Change during Visuomotor Adaptation Relates to Adaptation Performance and Cerebellar Network Connectivity: A Magnetic Resonance Spectroscopic Imaging Study.

Motor adaptation is crucial for performing accurate movements in a changing environment and relies on the cerebellum. Although cerebellar involvement has been well characterized, the neurochemical changes in the cerebellum underpinning human motor adaptation remain unknown. We used a novel magnetic resonance spectroscopic imaging (MRSI) technique to measure changes in the inhibitory neurotransmitter GABA in the human cerebellum during visuomotor adaptation. Participants (n = 17, six female) used their right hand to adapt to a rotated cursor in the scanner, compared with a control task requiring no adaptation. We spatially resolved adaptation-driven GABA changes at the cerebellar nuclei and cerebellar cortex in the left and the right cerebellar hemisphere independently and found that simple right-hand movements increase GABA in the right cerebellar nuclei and decreases GABA in the left. When isolating adaptation-driven GABA changes, we found that GABA in the left cerebellar nuclei and the right cerebellar nuclei diverged, although GABA change from baseline at the right cerebellar nuclei was not different from zero at the group level. Early adaptation-driven GABA fluctuations in the right cerebellar nuclei correlated with adaptation performance. Participants showing greater GABA decrease adapted better, suggesting early GABA change is behaviorally relevant. Early GABA change also correlated with functional connectivity change in a cerebellar network. Participants showing greater decreases in GABA showed greater strength increases in cerebellar network connectivity. Results were specific to GABA, to adaptation, and to the cerebellar network. This study provides first evidence for plastic changes in cerebellar neurochemistry during motor adaptation. Characterizing these naturally occurring neurochemical changes may provide a basis for developing therapeutic interventions to facilitate human motor adaptation.SIGNIFICANCE STATEMENT Despite motor adaptation being fundamental to maintaining accurate movements, its neurochemical basis remains poorly understood, perhaps because measuring neurochemicals in the human cerebellum is technically challenging. Using a novel magnetic resonance spectroscopic imaging method, this study provides evidence for GABA changes in the left compared with the right cerebellar nuclei driven by both simple movement and motor adaptation. Although right cerebellar GABA changes were not significantly different from zero at the group level, the adaptation-driven GABA fluctuations in the right cerebellar nuclei correlated with adaptation performance and with functional connectivity change in a cerebellar network. These results show the first evidence for plastic changes in cerebellar neurochemistry during a cerebellar learning task. This provides the basis for developing therapeutic interventions that facilitate these naturally occurring changes to amplify cerebellar-dependent learning.

GABA and Glutamate in hMT+ Link to Individual Differences in Residual Visual Function After Occipital Stroke.

BACKGROUND: Damage to the primary visual cortex following an occipital stroke causes loss of conscious vision in the contralateral hemifield. Yet, some patients retain the ability to detect moving visual stimuli within their blind field. The present study asked whether such individual differences in blind field perception following loss of primary visual cortex could be explained by the concentration of neurotransmitters γ-aminobutyric acid (GABA) and glutamate or activity of the visual motion processing, human middle temporal complex (hMT+). METHODS: We used magnetic resonance imaging in 19 patients with chronic occipital stroke to measure the concentration of neurotransmitters GABA and glutamate (proton magnetic resonance spectroscopy) and functional activity in hMT+ (functional magnetic resonance imaging). We also tested each participant on a 2-interval forced choice detection task using high-contrast, moving Gabor patches. We then measured and assessed the strength of relationships between participants' residual vision in their blind field and in vivo neurotransmitter concentrations, as well as visually evoked functional magnetic resonance imaging activity in their hMT+. Levels of GABA and glutamate were also measured in a sensorimotor region, which served as a control. RESULTS: Magnetic resonance spectroscopy-derived GABA and glutamate concentrations in hMT+ (but not sensorimotor cortex) strongly predicted blind-field visual detection abilities. Performance was inversely related to levels of both inhibitory and excitatory neurotransmitters in hMT+ but, surprisingly, did not correlate with visually evoked blood oxygenation level-dependent signal change in this motion-sensitive region. CONCLUSIONS: Levels of GABA and glutamate in hMT+ appear to provide superior information about motion detection capabilities inside perimetrically defined blind fields compared to blood oxygenation level-dependent signal changes-in essence, serving as biomarkers for the quality of residual visual processing in the blind-field. Whether they also reflect a potential for successful rehabilitation of visual function remains to be determined.

Event-related functional magnetic resonance spectroscopy.

Proton-Magnetic Resonance Spectroscopy (MRS) is a non-invasive brain imaging technique used to measure the concentration of different neurochemicals. "Single-voxel" MRS data is typically acquired across several minutes, before individual transients are averaged through time to give a measurement of neurochemical concentrations. However, this approach is not sensitive to more rapid temporal dynamics of neurochemicals, including those that reflect functional changes in neural computation relevant to perception, cognition, motor control and ultimately behaviour. In this review we discuss recent advances in functional MRS (fMRS) that now allow us to obtain event-related measures of neurochemicals. Event-related fMRS involves presenting different experimental conditions as a series of trials that are intermixed. Critically, this approach allows spectra to be acquired at a time resolution in the order of seconds. Here we provide a comprehensive user guide for event-related task designs, choice of MRS sequence, analysis pipelines, and appropriate interpretation of event-related fMRS data. We raise various technical considerations by examining protocols used to quantify dynamic changes in GABA, the primary inhibitory neurotransmitter in the brain. Overall, we propose that although more data is needed, event-related fMRS can be used to measure dynamic changes in neurochemicals at a temporal resolution relevant to computations that support human cognition and behaviour.

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 (31 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 31 P-MRSI could make cardiac stress protocols or creatine kinase rate measurements (which involve repeated scans) more tolerable for patients without sacrificing spatial resolution.

Continuous Automated Analysis Workflow for MRS Studies.

Magnetic resonance spectroscopy (MRS) can non-invasively measure levels of endogenous metabolites in living tissue and is of great interest to neuroscience and clinical research. To this day, MRS data analysis workflows differ substantially between groups, frequently requiring many manual steps to be performed on individual datasets, e.g., data renaming/sorting, manual execution of analysis scripts, and manual assessment of success/failure. Manual analysis practices are a substantial barrier to wider uptake of MRS. They also increase the likelihood of human error and prevent deployment of MRS at large scale. Here, we demonstrate an end-to-end workflow for fully automated data uptake, processing, and quality review.The proposed continuous automated MRS analysis workflow integrates several recent innovations in MRS data and file storage conventions. They are efficiently deployed by a directory monitoring service that automatically triggers the following steps upon arrival of a new raw MRS dataset in a project folder: (1) conversion from proprietary manufacturer file formats into the universal format NIfTI-MRS; (2) consistent file system organization according to the data accumulation logic standard BIDS-MRS; (3) executing a command-line executable of our open-source end-to-end analysis software Osprey; (4) e-mail delivery of a quality control summary report for all analysis steps.The automated architecture successfully completed for a demonstration dataset. The only manual step required was to copy a raw data folder into a monitored directory.Continuous automated analysis of MRS data can reduce the burden of manual data analysis and quality control, particularly for non-expert users and multi-center or large-scale studies and offers considerable economic advantages.

Uncertainty in denoising of MRSI using low-rank methods.

PURPOSE: Low-rank denoising of MRSI data results in an apparent increase in spectral SNR. However, it is not clear if this translates to a lower uncertainty in metabolite concentrations after spectroscopic fitting. Estimation of the true uncertainty after denoising is desirable for downstream analysis in spectroscopy. In this work, the uncertainty reduction from low-rank denoising methods based on spatiotemporal separability and linear predictability in MRSI are assessed. A new method for estimating metabolite concentration uncertainty after denoising is proposed. Automatic rank threshold selection methods are also assessed in simulated low SNR regimes. METHODS: Assessment of denoising methods is conducted using Monte Carlo simulation of proton MRSI data and by reproducibility of repeated in vivo acquisitions in 5 subjects. RESULTS: In simulated and in vivo data, spatiotemporal based denoising is shown to reduce the concentration uncertainty, but linear prediction denoising increases uncertainty. Uncertainty estimates provided by fitting algorithms after denoising consistently underestimate actual metabolite uncertainty. However, the proposed uncertainty estimation, based on an analytical expression for entry-wise variance after denoising, is more accurate. It is also shown automated rank threshold selection using Marchenko-Pastur distribution can bias the data in low SNR conditions. An alternative soft-thresholding function is proposed. CONCLUSION: Low-rank denoising methods based on spatiotemporal separability do reduce uncertainty in MRS(I) data. However, thorough assessment is needed as assessment by SNR measured from residual baseline noise is insufficient given the presence of non-uniform variance. It is also important to select the right rank thresholding method in low SNR cases.

NIfTI-MRS: A standard data format for magnetic resonance spectroscopy.

PURPOSE: Multiple data formats in the MRS community currently hinder data sharing and integration. NIfTI-MRS is proposed as a standard spectroscopy data format, implemented as an extension to the Neuroimaging informatics technology initiative (NIfTI) format. This standardized format can facilitate data sharing and algorithm development as well as ease integration of MRS analysis alongside other imaging modalities. METHODS: A file format using the NIfTI header extension framework incorporates essential spectroscopic metadata and additional encoding dimensions. A detailed description of the specification is provided. An open-source command-line conversion program is implemented to convert single-voxel and spectroscopic imaging data to NIfTI-MRS. Visualization of data in NIfTI-MRS is provided by development of a dedicated plugin for FSLeyes, the FMRIB Software Library (FSL) image viewer. RESULTS: Online documentation and 10 example datasets in the proposed format are provided. Code examples of NIfTI-MRS readers are implemented in common programming languages. Conversion software, spec2nii, currently converts 14 formats where data is stored in image-space to NIfTI-MRS, including Digital Imaging and Communications in Medicine (DICOM) and vendor proprietary formats. CONCLUSION: NIfTI-MRS aims to solve issues arising from multiple data formats being used in the MRS community. Through a single conversion point, processing and analysis of MRS data are simplified, thereby lowering the barrier to use of MRS. Furthermore, it can serve as the basis for open data sharing, collaboration, and interoperability of analysis programs. Greater standardization and harmonization become possible. By aligning with the dominant format in neuroimaging, NIfTI-MRS enables the use of mature tools present in the imaging community, demonstrated in this work by using a dedicated imaging tool, FSLeyes, for visualization.

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.

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.

Myocardial Energetics in Obesity: Enhanced ATP Delivery Through Creatine Kinase With Blunted Stress Response.

BACKGROUND: Obesity is strongly associated with exercise intolerance and the development of heart failure. Whereas myocardial energetics and diastolic function are impaired in obesity, systolic function is usually preserved. This suggests that the rate of ATP delivery is maintained, but this has never been explored in human obesity. We hypothesized that ATP transfer rate through creatine kinase (CK) (kfCKrest) would be increased, compensating for depleted energy stores (phosphocreatine/ATP), but potentially limiting greater ATP delivery during increased workload. We hypothesized that these changes would normalize with weight loss. METHODS: We recruited 80 volunteers (35 controls [body mass index 24±3 kg/m2], 45 obese [body mass index 35±5 kg/m2]) without coexisting cardiovascular disease. Participants underwent body composition analysis, magnetic resonance imaging of abdominal, liver, and myocardial fat content, left ventricular function, and 31P magnetic resonance spectroscopy to assess phosphocreatine/ATP and CK kinetics, at rest and during dobutamine stress. Obese volunteers were assigned to a dietary weight loss intervention, before reexamination. RESULTS: At rest, although myocardial phosphocreatine/ATP was 14% lower in obesity (1.9±0.3 versus 2.2±0.2, P<0.001), kfCkrest was 33% higher (0.23±0.07 s-1 versus 0.16±0.08 s-1, P=0.002), yielding no difference in overall resting ATP delivery (obese 2.5±0.9 µmol·g-1·s-1 versus control 2.2±1.1 µmol·g-1·s-1, P=0.232). In controls, increasing cardiac workload led to an increase in both kfCK (+86%, P<0.001) and ATP delivery (+80%, P<0.001). However, in obesity, similar stress led to no significant increase in either kfCK (P=0.117) or ATP delivery (P=0.608). This was accompanied by reduced systolic augmentation (absolute increase in left ventricular ejection fraction, obese +16±7% versus control +21±4%, P=0.031). Successful weight loss (-11±5% body weight) was associated with improvement of these energetic changes such that there was no significant difference in comparison with controls. CONCLUSIONS: In the obese resting heart, the myocardial CK reaction rate is increased, maintaining ATP delivery despite reduced phosphocreatine/ATP. During increased workload, although the nonobese heart increases ATP delivery through CK, the obese heart does not; this is associated with reduced systolic augmentation and exercise tolerance. Weight loss reverses these energetic changes. This highlights myocardial energy delivery through CK as a potential therapeutic target to improve symptoms in obesity-related heart disease, and a fascinating modifiable pathway involved in the progression to heart failure, as well.

Cardiac Energetics in Patients With Aortic Stenosis and Preserved Versus Reduced Ejection Fraction.

BACKGROUND: Why some but not all patients with severe aortic stenosis (SevAS) develop otherwise unexplained reduced systolic function is unclear. We investigate the hypothesis that reduced creatine kinase (CK) capacity and flux is associated with this transition. METHODS: We recruited 102 participants to 5 groups: moderate aortic stenosis (ModAS) (n=13), SevAS, left ventricular (LV) ejection fraction ≥55% (SevAS-preserved ejection fraction, n=37), SevAS, LV ejection fraction <55% (SevAS-reduced ejection fraction, n=15), healthy volunteers with nonhypertrophied hearts with normal systolic function (normal healthy volunteer, n=30), and patients with nonhypertrophied, non-pressure-loaded hearts with normal systolic function undergoing cardiac surgery and donating LV biopsy (non-pressure-loaded heart biopsy, n=7). All underwent cardiac magnetic resonance imaging and 31P magnetic resonance spectroscopy for myocardial energetics. LV biopsies (AS and non-pressure-loaded heart biopsy) were analyzed for CK total activity, CK isoforms, citrate synthase activity, and total creatine. Mitochondria-sarcomere diffusion distances were calculated by using serial block-face scanning electron microscopy. RESULTS: In the absence of failure, CK flux was lower in the presence of AS (by 32%, P=0.04), driven primarily by reduction in phosphocreatine/ATP (by 17%, P<0.001), with CK kf unchanged (P=0.46). Although lowest in the SevAS-reduced ejection fraction group, CK flux was not different from the SevAS-preserved ejection fraction group (P>0.99). Accompanying the fall in CK flux, total CK and citrate synthase activities and the absolute activities of mitochondrial-type CK and CK-MM isoforms were also lower (P<0.02, all analyses). Median mitochondria-sarcomere diffusion distances correlated well with CK total activity (r=0.86, P=0.003). CONCLUSIONS: Total CK capacity is reduced in SevAS, with median values lowest in those with systolic failure, consistent with reduced energy supply reserve. Despite this, in vivo magnetic resonance spectroscopy measures of resting CK flux suggest that ATP delivery is reduced earlier, at the moderate AS stage, where LV function remains preserved. These findings show that significant energetic impairment is already established in moderate AS and suggest that a fall in CK flux is not by itself a necessary cause of transition to systolic failure. However, because ATP demands increase with AS severity, this could increase susceptibility to systolic failure. As such, targeting CK capacity and flux may be a therapeutic strategy to prevent and treat systolic failure in AS.

Frequency drift in MR spectroscopy at 3T.

PURPOSE: Heating of gradient coils and passive shim components is a common cause of instability in the B0 field, especially when gradient intensive sequences are used. The aim of the study was to set a benchmark for typical drift encountered during MR spectroscopy (MRS) to assess the need for real-time field-frequency locking on MRI scanners by comparing field drift data from a large number of sites. METHOD: A standardized protocol was developed for 80 participating sites using 99 3T MR scanners from 3 major vendors. Phantom water signals were acquired before and after an EPI sequence. The protocol consisted of: minimal preparatory imaging; a short pre-fMRI PRESS; a ten-minute fMRI acquisition; and a long post-fMRI PRESS acquisition. Both pre- and post-fMRI PRESS were non-water suppressed. Real-time frequency stabilization/adjustment was switched off when appropriate. Sixty scanners repeated the protocol for a second dataset. In addition, a three-hour post-fMRI MRS acquisition was performed at one site to observe change of gradient temperature and drift rate. Spectral analysis was performed using MATLAB. Frequency drift in pre-fMRI PRESS data were compared with the first 5:20 minutes and the full 30:00 minutes of data after fMRI. Median (interquartile range) drifts were measured and showed in violin plot. Paired t-tests were performed to compare frequency drift pre- and post-fMRI. A simulated in vivo spectrum was generated using FID-A to visualize the effect of the observed frequency drifts. The simulated spectrum was convolved with the frequency trace for the most extreme cases. Impacts of frequency drifts on NAA and GABA were also simulated as a function of linear drift. Data from the repeated protocol were compared with the corresponding first dataset using Pearson's and intraclass correlation coefficients (ICC). RESULTS: Of the data collected from 99 scanners, 4 were excluded due to various reasons. Thus, data from 95 scanners were ultimately analyzed. For the first 5:20 min (64 transients), median (interquartile range) drift was 0.44 (1.29) Hz before fMRI and 0.83 (1.29) Hz after. This increased to 3.15 (4.02) Hz for the full 30 min (360 transients) run. Average drift rates were 0.29 Hz/min before fMRI and 0.43 Hz/min after. Paired t-tests indicated that drift increased after fMRI, as expected (p < 0.05). Simulated spectra convolved with the frequency drift showed that the intensity of the NAA singlet was reduced by up to 26%, 44 % and 18% for GE, Philips and Siemens scanners after fMRI, respectively. ICCs indicated good agreement between datasets acquired on separate days. The single site long acquisition showed drift rate was reduced to 0.03 Hz/min approximately three hours after fMRI. DISCUSSION: This study analyzed frequency drift data from 95 3T MRI scanners. Median levels of drift were relatively low (5-min average under 1 Hz), but the most extreme cases suffered from higher levels of drift. The extent of drift varied across scanners which both linear and nonlinear drifts were observed.

Comparison of Assays for Therapeutic Monitoring of Infliximab and Adalimumab in Patients With Inflammatory Bowel Diseases.

Comparison data regarding anti-tumor necrosis factor drug concentrations in inflammatory bowel disease (IBD) between the enzyme-linked immunosorbent assay (ELISA) and the homogenous mobility shift assay (HMSA) are scarce.1-3 As decisions in clinical practice depend on the thresholds that define a therapeutic drug concentration, it is important to determine if this varies based on the type of assay used for therapeutic drug monitoring.4 We recently showed a discrepancy between a commercially available ELISA and the HMSA for both infliximab and adalimumab concentrations in patients with IBD.5 Based on the results of the study, Prometheus Laboratories (San Diego, CA) initiated a comprehensive review of their HMSA assays and found that there was an upward drift for both infliximab (from December 2017 to May 2019) and adalimumab (from August 2017 to May 2019), including when our study was performed. Prometheus Laboratories corrected the errant values and reported the revised drug concentrations to physicians (Supplementary Methods). We aimed to compare the corrected infliximab and adalimumab concentrations with the original ELISA values.

FSL-MRS: An end-to-end spectroscopy analysis package.

PURPOSE: We introduce FSL-MRS, an end-to-end, modular, open-source MRS analysis toolbox. It provides spectroscopic data conversion, preprocessing, spectral simulation, fitting, quantitation, and visualization. METHODS: The FSL-MRS package is modular. Its programs operate on data in a standard format (Neuroimaging Informatics Technology Initiative [NIfTI]) capable of storing single-voxel and multivoxel spectroscopy, including spatial orientation information. The FSL-MRS toolbox includes tools for preprocessing of raw spectroscopy data, including coil combination, frequency and phase alignment, and filtering. A density matrix simulation program is supplied for generation of basis spectra from simple text-based descriptions of pulse sequences. Fitting is based on linear combination of basis spectra and implements Markov chain Monte Carlo optimization for the estimation of the full posterior distribution of metabolite concentrations. Validation of the fitting is carried out on independently created simulated data, phantom data, and three in vivo human data sets (257 single-voxel spectroscopy and 8 MRSI data sets) at 3 T and 7 T. Interactive HTML reports are automatically generated by processing and fitting stages of the toolbox. The FSL-MRS package can be used on the command line or interactively in the Python language. RESULTS: Validation of the fitting shows low error in simulation (median error of 11.9%) and in phantom (3.4%). Average correlation between a third-party toolbox (LCModel) and FSL-MRS was high (0.53-0.81) in all three in vivo data sets. CONCLUSION: The FSL-MRS toolbox is designed to be flexible and extensible to new forms of spectroscopic acquisitions. Custom fitting models can be specified within the framework for dynamic or multivoxel spectroscopy. It is available as part of the FMRIB Software Library.

Multi-site harmonization of 7 tesla MRI neuroimaging protocols.

Increasing numbers of 7 T (7 T) magnetic resonance imaging (MRI) scanners are in research and clinical use. 7 T MRI can increase the scanning speed, spatial resolution and contrast-to-noise-ratio of many neuroimaging protocols, but technical challenges in implementation have been addressed in a variety of ways across sites. In order to facilitate multi-centre studies and ensure consistency of findings across sites, it is desirable that 7 T MRI sites implement common high-quality neuroimaging protocols that can accommodate different scanner models and software versions. With the installation of several new 7 T MRI scanners in the United Kingdom, the UK7T Network was established with an aim to create a set of harmonized structural and functional neuroimaging sequences and protocols. The Network currently includes five sites, which use three different scanner platforms, provided by two different vendors. Here we describe the harmonization of functional and anatomical imaging protocols across the three different scanner models, detailing the necessary changes to pulse sequences and reconstruction methods. The harmonized sequences are fully described, along with implementation details. Example datasets acquired from the same subject on all Network scanners are made available. Based on these data, an evaluation of the harmonization is provided. In addition, the implementation and validation of a common system calibration process is described.

Multi-centre, multi-vendor reproducibility of 7T QSM and R2* in the human brain: Results from the UK7T study.

INTRODUCTION: We present the reliability of ultra-high field T2* MRI at 7T, as part of the UK7T Network's "Travelling Heads" study. T2*-weighted MRI images can be processed to produce quantitative susceptibility maps (QSM) and R2* maps. These reflect iron and myelin concentrations, which are altered in many pathophysiological processes. The relaxation parameters of human brain tissue are such that R2* mapping and QSM show particularly strong gains in contrast-to-noise ratio at ultra-high field (7T) vs clinical field strengths (1.5-3T). We aimed to determine the inter-subject and inter-site reproducibility of QSM and R2* mapping at 7T, in readiness for future multi-site clinical studies. METHODS: Ten healthy volunteers were scanned with harmonised single- and multi-echo T2*-weighted gradient echo pulse sequences. Participants were scanned five times at each "home" site and once at each of four other sites. The five sites had 1× Philips, 2× Siemens Magnetom, and 2× Siemens Terra scanners. QSM and R2* maps were computed with the Multi-Scale Dipole Inversion (MSDI) algorithm (https://github.com/fil-physics/Publication-Code). Results were assessed in relevant subcortical and cortical regions of interest (ROIs) defined manually or by the MNI152 standard space. RESULTS AND DISCUSSION: Mean susceptibility (χ) and R2* values agreed broadly with literature values in all ROIs. The inter-site within-subject standard deviation was 0.001-0.005 ppm (χ) and 0.0005-0.001 ms-1 (R2*). For χ this is 2.1-4.8 fold better than 3T reports, and 1.1-3.4 fold better for R2*. The median ICC from within- and cross-site R2* data was 0.98 and 0.91, respectively. Multi-echo QSM had greater variability vs single-echo QSM especially in areas with large B0 inhomogeneity such as the inferior frontal cortex. Across sites, R2* values were more consistent than QSM in subcortical structures due to differences in B0-shimming. On a between-subject level, our measured χ and R2* cross-site variance is comparable to within-site variance in the literature, suggesting that it is reasonable to pool data across sites using our harmonised protocol. CONCLUSION: The harmonized UK7T protocol and pipeline delivers on average a 3-fold improvement in the coefficient of reproducibility for QSM and R2* at 7T compared to previous reports of multi-site reproducibility at 3T. These protocols are ready for use in multi-site clinical studies at 7T.