Skip to main content

ORIGINAL RESEARCH article

Front. Neurosci., 07 December 2023
Sec. Brain Imaging Methods
This article is part of the Research Topic Diffusion-Weighted MR Imaging (DW-MRI) and Diffusion-Weighted MR Spectroscopy (DW-MRS) View all 9 articles

Practical considerations of diffusion-weighted MRS with ultra-strong diffusion gradients

  • 1The Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, MD, United States
  • 2Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, MD, United States
  • 3Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, United Kingdom
  • 4CIBM Center for Biomedical Imaging, EPFL CIBM-AIT, EPFL Lausanne, Lausanne, Switzerland
  • 5Siemens Healthcare Ltd., Camberly, United Kingdom
  • 6Department of Radiology, Universität Bern, Bern, Switzerland
  • 7Leeds Institute of Cardiovascular & Metabolic Medicine, University of Leeds, Leeds, United Kingdom
  • 8Clinical Sciences Institue, Brighton and Sussex Medical School, Brighton, United Kingdom
  • 9Center for NeuroImaging Research (CENIR), Paris Brain Institute (ICM), Pitié-Salpêtrière Hospital, Paris, France
  • 10Inserm U1127, CNRS U7225, Sorbonne Universités, Paris, France
  • 11Brain Research Imaging Centre, School Physics and Astronomy, Cardiff University, Cardiff, United Kingdom
  • 12Image Sciences Institute, University Medical Center Utrecht, Utrecht, Netherlands

Introduction: Diffusion-weighted magnetic resonance spectroscopy (DW-MRS) offers improved cellular specificity to microstructure—compared to water-based methods alone—but spatial resolution and SNR is severely reduced and slow-diffusing metabolites necessitate higher b-values to accurately characterize their diffusion properties. Ultra-strong gradients allow access to higher b-values per-unit time, higher SNR for a given b-value, and shorter diffusion times, but introduce additional challenges such as eddy-current artefacts, gradient non-uniformity, and mechanical vibrations.

Methods: In this work, we present initial DW-MRS data acquired on a 3T Siemens Connectom scanner equipped with ultra-strong (300 mT/m) gradients. We explore the practical issues associated with this manner of acquisition, the steps that may be taken to mitigate their impact on the data, and the potential benefits of ultra-strong gradients for DW-MRS. An in-house DW-PRESS sequence and data processing pipeline were developed to mitigate the impact of these confounds. The interaction of TE, b-value, and maximum gradient amplitude was investigated using simulations and pilot data, whereby maximum gradient amplitude was restricted. Furthermore, two DW-MRS voxels in grey and white matter were acquired using ultra-strong gradients and high b-values.

Results: Simulations suggest T2-based SNR gains that are experimentally confirmed. Ultra-strong gradient acquisitions exhibit similar artefact profiles to those of lower gradient amplitude, suggesting adequate performance of artefact mitigation strategies. Gradient field non-uniformity influenced ADC estimates by up to 4% when left uncorrected. ADC and Kurtosis estimates for tNAA, tCho, and tCr align with previously published literature.

Discussion: In conclusion, we successfully implemented acquisition and data processing strategies for ultra-strong gradient DW-MRS and results indicate that confounding effects of the strong gradient system can be ameliorated, while achieving shorter diffusion times and improved metabolite SNR.

1 Introduction

Diffusion-weighted magnetic resonance imaging (DW-MRI) is usually sensitized to the displacement of water, and provides a myriad of tissue microstructure metrics that aid in the study of many neuropathologies, including traumatic brain injury (Hutchinson et al., 2018), neurodegeneration (Goveas et al., 2015), and measuring treatment response in cancer therapy (Patterson et al., 2008), to name a few. However, the ubiquity of water molecules—present in both intra- and extracellular spaces—complicates the modeling of water diffusion as a measure of cellular microstructure. Magnetic resonance spectroscopy (MRS) is a non-invasive technique providing quantitative measures of metabolites and neurotransmitters which are present in the brain at millimolar concentrations. Diffusion-weighted MRS (DW-MRS) introduces diffusion gradients into MRS sequences, utilising MRS as a filter to sensitize the MR signal to different metabolites which are almost exclusively intra-cellular, with some considered predominantly glial—myo-inositol (mI) and choline compounds (tCho)—and others predominantly neuronal—N-acetyl-aspartate (NAA) and glutamate (Glu) (Choi et al., 2007). While the signal-to-noise ratio (SNR) and spatial resolution are reduced compared to conventional water-based imaging, the specificity afforded by DW-MRS greatly simplifies diffusion modeling and interpretation, and provides a valuable non-invasive window into metabolism and cellular microstructure, complimentary to water-based diffusion imaging (Ronen et al., 2014; Najac et al., 2016; Palombo et al., 2018; Ligneul et al., 2019; Genovese et al., 2021b).

The apparent diffusion coefficients (ADCs) of metabolites are at least five times smaller than those of water (Ellegood et al., 2011), which necessitates higher b-values to adequately characterize metabolite diffusion properties. A common approach is to employ a DW-STEAM sequence, whereby metabolite diffusion occurs during the mixing time (TM), with the diffusion time uncoupled from the echo time (TE). However, this comes with the caveat that STEAM generates a stimulated echo, reducing the available SNR compared to spin-echo localisation methods by a factor of two. The shorter TE of STEAM can ameliorate this, but the long diffusion times required to achieve adequate diffusion weighting may be undesired. If the goal is to probe short diffusion times and/or high b-values, then the spin-echo-based diffusion-weighted Point RESolved Spectroscopy sequence (Bottomley, 1987) (DW-PRESS) provides an alternative, offering better SNR than STEAM, without the additional TE restrictions imposed by adiabatic pulse pairing required for a diffusion-weighted semi-LASER (DW-sLASER) sequence. Larger b-values can be achieved by increasing the DW-gradient amplitude (limited by the gradient system) and/or by increasing the diffusion time. Achieving the latter with DW-PRESS can be challenging. With the diffusion time coupled to the choice of TE, the available SNR at high b-value is restricted by metabolite T2 relaxation. The introduction of ultra-strong gradient systems can mitigate this. The Siemens Connectom scanner is fitted with a gradient system capable of reaching 300 mT/m per axis. This provides larger b-values for a given TE, and access to shorter diffusion times while maintaining the required b-value range. Shorter diffusion times can reduce the variability resulting from motion artefacts and, crucially, can provide additional cell-specific microstructural properties of highly-restrictive compartments (e.g., subcellular organelles) and cellular viscosity (Setsompop et al., 2013; Palombo et al., 2017; Jones et al., 2018).

However, the introduction of ultra-strong gradients poses additional practical challenges. Specifically, eddy currents become increasingly prevalent at larger gradient amplitudes. Switched gradient fields induce eddy currents which produce time-varying magnetic fields, distorting the lineshape of MR spectra and hampering MRS modeling attempts. Moreover, eddy current correction in conventional MRS relies on acquiring a water-unsuppressed reference scan (Klose, 1990); however, in DW-MRS, the water signal is heavily attenuated at high b-values, complicating the extraction of the relevant phase information. In addition to eddy current effects, it becomes increasingly difficult to maintain uniform gradient fields on ultra-strong gradient systems. This leads to a deviation in the applied gradient field from the nominal gradient field, which is of particular relevance for diffusion studies (Mesri et al., 2020). These gradient non-uniformities will spatially modulate the b-matrix and image geometry, and must be corrected in order to obtain reliable estimates (Bammer et al., 2003). Finally, mechanical vibrations—caused by Lorentz forces generated during rapid gradient switching—can cause anomalous signal loss in diffusion experiments (Gallichan et al., 2010; Weidlich et al., 2019, 2020), and is a particular concern for strong gradient systems.

At the time of writing, only one full study has been published on DW-MRS with ultra-strong gradients in humans, which focused on utilising the hardware to measure macromolecular background profiles for MRS (Şimşek et al., 2022). Despite increasing community interest in strong gradients for diffusion encoding (Setsompop et al., 2013; Jones et al., 2018; Jenkins et al., 2020; Şimşek et al., 2020, 2021; Huang et al., 2021; Jenkins, 2021; Fan et al., 2022; Döring et al., 2023), no studies have specifically addressed the challenges of ultra-strong gradients in the context of DW-MRS. Furthermore, data processing software specific to DW-MRS data is limited, and variation between existing methods has been previously shown (Najac et al., 2022). With increased accessibility of high performance scanner gradient systems, it is crucial to develop open-source software that broaches the challenges posed by this hardware.

In this work, we present initial data acquired using a 3T Connectom—a research-only scanner equipped with a 300 mT/m gradient system. We introduce the practical issues associated with this manner of acquisition and steps that may be taken to mitigate their impact on the data. In this study, we limited the scope of our investigation specifically to the diffusion-weighted PRESS sequence, which was implemented with a flexible bipolar diffusion gradient scheme. A DW-MRS data processing pipeline was implemented and evaluated against pilot DW-PRESS data. Phantom experiments were conducted and in vivo data were acquired from a small cohort of individuals in order to validate the measurements and demonstrate the capabilities of ultra-strong gradients compared to more conventional gradient settings.

2 Materials and methods

For the DW-MRS acquisition, an in-house developed DW-PRESS sequence with bipolar diffusion gradients (Branzoli, 2015) was adapted for use on the 300 mT/m Connectom scanner. The full pulse sequence diagram is shown in Figure 1A.

FIGURE 1
www.frontiersin.org

Figure 1. (A) Pulse sequence diagram for the DW-PRESS sequence. Slice-select gradients are shown for each dimension, x (blue), y (red), z (yellow), and diffusion gradients are shown in gray. (B) The maximum achievable b-value as a function of the echo time for the DW-PRESS sequence, which is limited by the maximum achievable diffusion-gradient strengths. The curves were calculated for maximal gradient amplitudes of 40, 80, and 300 mT/m corresponding to the Siemens TRIO, Prisma, and Connectom scanner configurations, respectively. (C) For a given b-value, reduction of the echo time by going from 80 mT/m to 300 mT/m will result in SNR gain. SNR gains were estimated for tNAA, tCho and tCr using their relaxation times T2 (Ke et al., 2002; Ganji et al., 2012).

2.1 SNR simulations

We first conducted simulations to explore the relationship between TE, maximum b-value, and diffusion time for the DW-PRESS sequence, and subsequently, the potential TE-based SNR improvements possible with ultra-strong gradients. The maximum-achievable b-value is constrained by the gradient characteristics, and has a corresponding minimum-achievable TE. We investigated this relationship at three maximum gradient amplitudes: 40, 80, and 300 mT/m. The maximum achievable b-value was calculated numerically using in-house code written in Matlab. The slew rates—corresponding to each of the maximum gradient amplitudes—were 200 T/m/s, 200 T/m/s and 83 T/m/s, respectively.

For a given maximum b-value, the minimum-achievable TE will define the magnitude of T2-based signal attenuation, with short TE implying higher SNR. We investigated the expected SNR increase when moving from a Siemens Prisma gradient system (80mT/m) to a Connectom gradient system (300mT/m) for three major metabolite measures—creatine + phosphocreatine (tCr), tCho, and NAA+N-acetyl-aspartyl-glutamate (tNAA). The ratio of SNR values between the 300 mT/m and 80 mT/m gradient systems was calculated using previously-published metabolite T2 values (Ke et al., 2002; Ganji et al., 2012):

S(TE)~S0e-TE/T2,metabolite    (1)
SNR300SNR80(b)=S(TEmin,300(b))S(TEmin,80(b))=exp-TEmin,300(b)+TEmin,80(b)T2,metabolite    (2)

where, SNR300SNR80(b) is the ratio of TE-dependent SNR values corresponding to the respective gradient systems, S is the signal amplitude, S0 is the signal amplitude before T2-weighting, and T2, metabolite is the metabolite-specific T2. Note that the SNR ratio is independent of S0. The SNR ratio was calculated for a range of maximum b-values (and hence, TEs) for each metabolite. Apart from the pure T2 relaxation-driven reduction in SNR at longer TE, J-evolution contributes to metabolite dephasing and increases fitting uncertainties as well (Landheer et al., 2020); however, these effects were not included in our simple model, as we focused on the major singlet resonances.

2.2 Data acquisition

In order to evaluate measurement procedures and data processing strategies pilot data were acquired in small test group of 3 healthy participants—1 female; age 33 ± 10 years (mean ± standard deviation). Spectroscopic and structural data were collected on a research-only 3T Connectom MRI scanner, a modified 3T MAGNETOM Skyra system equipped with 300 mT/m gradients (Siemens Healthcare, Erlangen, Germany) and a 32-channel receive array head coil (Nova Medical, Wilmington, United States). The study was performed with ethics approval from the Cardiff University School of Psychology ethics review board and written informed consent was obtained from all participants. Additional Supplementary material were acquired on an isotropic diffusion phantom (NIST) (Palacios et al., 2017).

2.2.1 Structural MRI

Each in vivo MR protocol included a Magnetization Prepared Rapid Acquisition Gradient Echoes (MPRAGE) (Mugler III and Brookeman, 1990) scan. Scanning parameters for the MPRAGE sequence were: 1mm isotropic resolution, TR = 2.3s, TE = 2ms, TI = 837ms, Flip-angle = 9°, field-of-view=256 × 256.

2.2.2 DW-MRS

To reduce the impact of pulsation artefacts, all in vivo DW-MRS scans were cardiac triggered to avoid systolic pulsation using a pulse oximeter placed on the participant's forefinger. DW-MRS voxels were positioned using the T1-weighted MPRAGE.

In order to avoid cross-term effects from gradient overlap, diffusion gradients were applied along the physical single-gradient axes. To facilitate this, no rotations were applied to the voxels and the dimensions were adjusted to maximize the coverage of the relevant tissue volume in each case.

2.2.2.1 300 mT/m vs. 80 mT/m maximum gradient strength

The DW-PRESS sequence's functionality was extended to allow precise specification of the gradients applied. This allowed us to simulate different gradient configurations on the Connectom scanner by selecting appropriate minimum rise times and the maximum gradient amplitudes, Gmax. While this single-scanner experimental design limited the conclusions we could draw about the performance of the Connectom compared to a different system with lower maximum gradient strength, our approach allowed comparisons to be drawn between different gradient strengths without the confounding impact of certain experimental influences. Specifically, this approach removed the risk of different voxel placement, ensured the equivalence of scanner calibrations, and allowed us to study gradient-strength differences without introducing differences in gradient non-uniformities. Here, we considered conservative limits for the gradient characteristics [slightly below the hardware and physiological limits (Setsompop et al., 2013)] for the Connectom system, which can be potentially extended and thus the maximum-achievable b-values increased.

For the DW-MRS protocol, we selected TEs of 116ms and 74ms for Gmax = 80 mT/m (Prisma-like setting) and Gmax = 300 mT/m, respectively. Diffusion weighting was applied along a single gradient axis. The diffusion-gradient configurations provided as {δ, ϵ, a, Δ [ms]} and shown in Figure 1A were {20.7, 0.5, 9.2, 56.2} and {7.1, 3.6, 9.2, 35.2} for 80 mT/m and 300 mT/m, respectively, giving a maximum b-value of ≈30000s/mm2. Further specification of the sequence parameters are listed in Table 1. TEs and flat-top gradient durations were rounded to the nearest ms by the sequence implementation, resulting in unavoidable differences in b-values between the two conditions. However, these deviations were minimized to the greatest extent possible, and less than 1.5% in all cases. For each gradient setting, 24 transients were acquired with water suppression, along with 8 water-unsuppressed transients. A single 22 × 20 × 22 mm3 DW-MRS voxel was acquired for both gradient conditions in the grey-matter rich occipital lobe (OCC) of one participant in a single scan session.

TABLE 1
www.frontiersin.org

Table 1. Sequence parameters for the 300 and 80 mT/m gradient conditions.

2.2.2.2 Gray matter vs. white matter

To compare tissue-type specific differences, DW-MRS voxels were acquired in grey matter (GM) and white matter (WM) rich brain regions of the other two participants with a modified protocol. In one participant, a 24 × 20 × 25 mm3 voxel was placed in the OCC, centred on the mid-line and as posterior as possible without including the sagittal sinus or skull, maximising the GM fraction. In a second participant, a 27 × 16 × 16 mm3 voxel was placed in the sub-cortical white matter of the corona radiata (CR) in the left hemisphere, maximising the WM fraction. Contrary to the previous section, three diffusion weightings were applied along each of the physical gradient axes with a single maximum gradient amplitude (295 mT/m), and reduced TE of 70 ms giving a gradient configuration of {7.1, 3.6, 9.2, 31.2} and a reduced maximum b-value of 21763 s/mm2. The 8 nominal b-values were 0, 557.8, 1266.2, 2260.6, 5107.4, 9099.3, 14237, and 21763 s/mm2 along each of the three orthogonal single-gradient axes. Both DW-MRS acquisitions otherwise used the same parameters: TRmin = 2, 500 ms, spectral width = 4, 000Hz, 24 water-suppressed, and 8 water-unsuppressed transients.

2.2.2.3 Table vibration phantom scans

In order to investigate table vibrations induced by the gradient system, we acquired non-water-suppressed DW-MRS spectra in an isotropic diffusion phantom with a volume fraction of 50% polyvinylpyrrolidone (PVP, NIST). Spectra were acquired for the 300 mT/m and 80 mT/m gradient settings, with three diffusion directions, along each of the three orthogonal single-gradient physical axes. For each b-value, 4 transients were acquired, as well as 4 corresponding transients with inverted gradient polarity i.e., g = [0, ±59, ±89, ±118, ±177, ±236, ±290] mT/m. The inverted scans facilitated eddy current correction for the non-water-suppressed data, as described in the following section.

2.3 Data processing

2.3.1 MRS data processing

A DW-MRS pre-processing pipeline—conforming to MRS consensus recommendations (Near et al., 2021)—was implemented using the Matlab-based MRS toolkit, FID-A (Simpson et al., 2017)1. Relative coil phasing was applied using the water-unsuppressed b = 0 acquisition, weighting individual coil elements based on signal-to-noise ratio (Hall et al., 2014). Motion-corrupted transients were identified using a likeliness metric, comparing FIDs to the first acquisition for each respective diffusion condition (Simpson et al., 2017). Transients which varied by more than 2 standard deviations were omitted prior to averaging to reduce the impact of motion on the final results. To minimize signal losses due to phase and frequency drift, spectral registration (Near et al., 2015) was used to align individual transients for each diffusion condition, separately. The resulting spectra were then manually inspected for residual water/lipid or motion contamination, automated data quality cutoffs were used for tNAA full-width at half-maximum greater than 0.1 PPM and metabolite SNR less than 3 (Wilson et al., 2019).

Tarquin (Wilson et al., 2011) V4.3.10 was used for linear-combination modeling (LCM), with TE-specific simulated basis sets including -CrCH2 (relaxation correction basis function), alanine (Ala), aspartate (Asp), creatine (Cr), GABA, glycerophosphocholine (GPC), glucose (Glc), glutamine (Gln), glutathione (GSH), Glu, glycine (Gly), mI, lactate (Lac), NAA, N-acetylaspartylglutamate (NAAG), phosphocholine (PCh), phosphocreatine (PCr), phosphorylethanolamine (PE), scyllo-inositol (sI), and taurine (Tau). Reported metabolites are total creatine (tCr = Cr + PCr), total NAA (tNAA = NAA + NAAG), and total choline (tCho = Cho + GPC). Additional basis functions were incorporated for macromolecular and lipid resonances, with the baseline approximated by a Gaussian window function. The extracted metabolite amplitudes and Cramér-Rao lower bounds (CRLBs) were then used for diffusion modeling. Water phantom scans were processed without LCM, with the water signal amplitude quantified by taking the magnitude of the first point of the complex FID.

2.3.2 Eddy current correction

While the DW-PRESS sequence was designed to minimize eddy current artefacts using a bipolar gradient scheme (Branzoli, 2015), further eddy-current correction was required, particularly for the highest b-values. For in vivo data, we performed eddy current correction using a non-water-suppressed reference scan. The signal of the water-suppressed scan, Smet(t) can be divided by the signal of the non-water-suppressed reference Sref(t) to remove the eddy-current-related phase, ϕeddy, (Klose, 1990):

Smet(t)=|Smet|ei(ϕmet(t)+ϕeddy(t))    (3)
Sref(t)=|Sref|eiϕeddy(t)    (4)

where the phase information of the MR spectrum, ϕmet, is retained. To mitigate issues resulting from the attenuation of the water signal at higher b-values, the relevant phase correction term, ϕeddy was extracted from the reference scan using an LPSVD (Vanhamme et al., 1998), improving robustness to the noisier high-b-value water transients.

The Klose et al. method was used for all in vivo data in our study. Testing revealed satisfactory performance, even at high b-value; however, this might not be the case for all types of acquisition. To supplement this approach, we also demonstrate a second method of post-hoc eddy current correction in a nickel-doped water phantom, using gradient polarity inversion (Lin et al., 1994), whereby alternating transients of a particular diffusion gradient strength are acquired with opposite gradient polarity. This reverses the phase development due to eddy currents, such that:

Sinvert(t)=|S(t)|eiϕ(t)-iϕeddy(t)    (5)
S(t)Sinvert(t)=|S(t)|2eiϕ(t)    (6)

where, ϕ is the phase, ϕeddy is the phase specifically resulting from eddy currents, and S and Sinvert are the signals resulting from the initial and inverted gradient polarities, respectively. Multiplying the resulting time-domain signals removes phase evolution due to eddy currents.

2.3.3 Gradient non-uniformity correction

A mask representing the DW-MRS voxel was created in the anatomical image space using SPM12 (Penny et al., 2011; Edden et al., 2014). A vendor-supplied spherical harmonic description of the spatial dependence of the field generated by each gradient coil was used to assess geometric deviations of the voxel. A coil tensor, L, was computed by taking the partial derivatives of the field—normalized by the nominal gradient strength—where the elements of the coil tensor contain the spatially varying deviations for each of the gradient axes. The coil tensor was subsequently used to compute the effective b-matrix and b-value from the nominal ones, i.e., Beff=LBLT and beff = trace(Beff) (Bammer et al., 2003). Nominal b-values were estimated using in-house software written in Matlab (Mathworks, Natick, MA, USA), then subsequently, the nominal b-values were corrected for gradient non-uniformities (Bammer et al., 2003). This gave rise to a distribution of corrected b-values for each diffusion condition—rather than the single nominal value. The corrected mask was used in conjunction with the corrected b-values, to ascertain the effective b-value within the DW-MRS voxel.

2.3.4 Diffusion measures

Data were fitted in the low (b < 3, 000smm−2) and intermediate (b < 9, 500smm−2) b-value range with a mono-exponential and kurtosis diffusion representation, respectively, and over the full b-value range with an astrocylinder model (fully dispersed cylinders). All diffusion fitting was performed in Matlab 2021b using trust-region reflective optimisation. The inverse of the metabolite CRLBs was used to weight each data point, to reduce the impact of individual poorly-fit spectra on quantification.

2.3.4.1 Diffusion representation

Metabolite ADCs were estimated for each direction independently. The reported metabolite amplitudes for each b-value were modeled using two approaches. Firstly, all b-values below 3, 000 s/mm2 were fit using a mono-exponential decay:

ln(S)=ln(S0)-b·ADCe    (7)

where ln(S) is the natural logarithm of the fitted metabolite amplitude, and b is the b-value. S0 and ADCe are the non-diffusion-weighted signal amplitude and apparent diffusion coefficient (where e indicates the ADC from a mono-exponential fit), respectively. Non-mono-exponential behavior was observed beyond this b-value range, so b-values up to 9, 500 s/mm2 were fit using the diffusion kurtosis representation (Jensen et al., 2005; Yablonskiy and Sukstanskii, 2010):

ln(S)=ln(S0)-b·ADCk+16(b·ADCk)2·K    (8)

where ADCk is the apparent diffusion coefficient from the kurtosis representation and k is the kurtosis. To stabilize fitting, ADCe estimates from the mono-exponential representation were used to initiate the fit of the kurtosis representation.

Following gradient non-uniformity correction, each single nominal b-value, bi, is replaced with a distribution of corrected b-values (denoted bi,effj). To investigate the impact this had on metabolite diffusion measures, we estimated mono-exponential fits considering the distribution of b-values within the DW-MRS voxel:

[lnS(b1)lnS(b1)lnS(bN)]=[1b1,eff11b1,effM1bN,effM]·[lnS0ADCd]    (9)

where ADCd is the apparent diffusion coefficient (d indicates that it is estimated from a distribution of b-values), and the signal of the ith b-value, S(bi), is repeated across M voxels at 1 mm3 isotropic resolution (i.e., the resolution of the MPRAGE).

2.3.4.2 Diffusion modeling

Code adapted from the Multidimensional diffusion MRI (MD-dMRI) analysis framework was used to implement a single-compartment model of fully-dispersed cylinders (astrocylinders) for bipolar diffusion encoding, in agreement with DW-MRS recommendation (Nilsson et al., 2017, 2018; Ligneul et al., 2023). To rate model performance and investigate ultra-high b-values with Connectom settings (higher SNR, shorter TE, and diffusion-time) the corrected Akaike information criterion (AICc) was used. The model parameters include the signal amplitude S0 at b = 0, the free diffusivity D0 and the cylinder radius RC with boundaries 0 ≤ RC ≤ 20μm. Moreover, fitting included outlier detection to remove data-points potentially affected by motion by iteratively removing a single b-value from the fitting and identify the fit with the lowest root-mean-square-error (c.f., Supplementary Figures S3, S4). In the case of multi-directional data, the directional average of the signal was used for fitting to mitigate effects from tissue anisotropy. The fitting uncertainties of the parameters were estimated by residual bootstrapping with 250 random noise realizations (Jelescu et al., 2022).

3 Results

3.1 SNR simulations

Figures 1B, C shows the results of the theoretical investigation into the relationship between maximum-achievable b-value and minimum-achievable TE for the DW-PRESS sequence, as well as the related SNR simulation results. As anticipated, ultra-strong gradients allow for a far shorter TE for a given b-value, a benefit which becomes even more apparent as the required maximum b-value increases. Moreover, one can see from the color bar that the minimum achievable diffusion time for a given maximum b-value is strongly related to the maximum gradient strength. In terms of SNR, we found that tCr—with the shortest T2 of the metabolites considered—had the largest SNR gains. For a maximum b-value of 31 000 s/mm2, our simulations suggest a potential SNR improvement of 50%. Likewise, for tNAA and tCho, our simulations suggest an expected SNR gain of about 10% and 21% moving from 80 mT/m to 300 mT/m, respectively.

3.2 Eddy currents

Figure 2 shows a summary of sequence validation results, including examples of gradient-polarity-inversion eddy current correction (ECC) in a phantom (Figure 2A), and examples of water-reference-based ECC in vivo (Figure 2B). Although bipolar diffusion gradients were applied in the water PVP phantom, Figure 2A shows clear eddy current artifacts for both positive and negative gradient polarity configurations. However, by combining both datasets, effects from eddy currents can be largely prevented. Apart from the benefits of limiting eddy-current artifacts, combining consecutive transients with inverted gradient polarity also reduces contributions from linear gradient cross-terms. This allows for an accurate estimation of the b-value, even without taking slice-selection, and crushing gradients into account. If no gradient polarity inversion is used, one has to calculate the real b-value directly from the full gradient chronograms (Mattiello et al., 1997). For the in vivo acquisitions, where only a single gradient polarity configuration was used, the effects from eddy currents are clearly visible from Figure 2B (red line, Before ECC). This effect is more pronounced at higher b-values, where steeper gradients give rise to stronger induction currents. However, when the acquired water reference is used for ECC (black line, After ECC) non-linear phase distortions from eddy-currents can be widely prevented. Despite large differences in diffusion gradient characteristics, MRS fit residuals were comparable between the 80 mT/m and 300 mT/m acquisitions (Supplementary Figure S1). The mean ratio between the fit residual and the noise level—as reported by Tarquin—was 2.98 and 3.04 for the 300 and 80 mT/m acquisitions, respectively. The comparability of the residuals alludes to a lack of significant differences in modeling performance, despite the larger eddy current artifacts induced when ultra-strong gradients are applied.

FIGURE 2
www.frontiersin.org

Figure 2. (A) Phantom water spectra acquired with a given gradient polarity (green), inverted gradient polarity (red), and the combined, eddy-current corrected spectrum (blue) (B) Three examples of spectra before (red) and after Klose eddy-current correction (black) at 50 mT/m (top), 150 mT/m (middle), and 295 mT/m (bottom). (C) DW-MRS quantification of the isotropic diffusion phantom. Individual diffusion directions (x, green; y, blue;, z, red) are plotted on top of each other for the 300 mT/m (dashed fit line above, “x”) and 80 mT/m (dotted fit line below, “o”) acquisitions. The individual quantified ADCs are also noted in the colour corresponding to the direction acquired. Note the strong overlap between directions, and vertical offset between gradient conditions.

3.3 Gradient non-uniformities

Deviations in the voxel mask geometry and b-values due to gradient non-uniformities were corrected by taking the nonlinear spatial gradient profiles into account. While voxel deformations made little-to-no difference to the voxel volume and position (< 1% deviation), b-values were affected significantly. This is perhaps expected, as unlike spatial encoding, b-values are driven by the squared gradient amplitude. Figure 3 shows the magnitude and spatial distribution of deviations from the nominally-specified b-value due to gradient non-uniformities. While the fractional deviation from the nominal b-value is constant, the actual deviation increased with gradient strength. A relatively narrow distribution is observed about zero at the lowest non-zero b-value and a much wider, non-zero-centered distribution is seen for the highest b-value. For the smallest nominal b-value of 567 s/mm2, the mean (standard deviation) of the distributions of corrected b-values in x, y, and z were 563 (5), 567 (5), and 559 (5) s/mm2, respectively, and for the largest nominal b-value of 22,107 s/mm2–21,948 (198), 22,137 (185), 21,802 (199) s/mm2, respectively. The mean deviation—calculated by taking subtracting nominal b-value from the corrected one—were −4, 0, and −8 s/mm2 for 557.8 s/mm2, and –159, 30, and –305 s/mm2 for a nominal value of 21,763 s/mm2.

FIGURE 3
www.frontiersin.org

Figure 3. Deviation from nominally specified b-values due to gradient non-uniformity. (A) The distribution of deviation from the nominal b-value within the voxel. From dark blue to yellow, the distributions correspond to nominal b-values of 567, 2,296, 5,188, 9,243, 14,462, and 22,107 s/mm2. Higher b-values exhibit broader distributions of deviations. (B) An example of the spatial distribution of deviations from a nominal b-value of 22,107 s/mm2. The rows correspond to different gradient directions (consistent with the first panel), and the columns represent orthogonal slices. The images in the top row are overlaid with the DW-MRS voxel location used to determine the distributions shown in the left panel.

Figure 4 shows in vivo mono-exponential fitting results if b-values were not corrected for non-uniformity (blue); corrected for non-uniformity, fitting the the full in-voxel distribution (red); and corrected for non-uniformity using the median of the distribution for fitting (green). After correction, we observed no significant difference in the estimated ADCs when either using the full b-value distribution or the median of that distribution, but overestimated ADCs in the absense of non-uniformity correction (approximately 4% faster diffusion). Thus, in the subsequent analysis, the median b-value was used to inform ADC, kurtosis, and microstructural fitting.

FIGURE 4
www.frontiersin.org

Figure 4. This figure shows three different approaches to representing the b-value during mono-exponential modelling. Directly fitting using the nominally-specified b-value (blue), fitting using the median of the distribution of gradient-non-linearity-corrected b-values (green), and fitting using the full distribution of corrected b-values (red). The resulting ADCs are shown in their respective colours. The (left) shows the tNAA diffusivity, while the (right) shows tCho.

3.4 Comparison of 80 and 300 mT/m acquisitions

Figure 5A shows the individual in vivo spectra and LCM fitting results left with the Connectom (Gmax = 300mT/m) and right with the Prisma (Gmax = 80mT/m) settings. Overall spectral SNR of all acquisitions was well above minimum consensus recommendations, even for the highest b-value. The mean SNR gain for tNAA—calculated by taking the ratio of the metabolite peak amplitude to the standard deviation of the noise—was 1.16, with a 16% increase in SNR when ultra-strong gradients are used. Likewise, for the other major metabolites, we found an SNR gain of 1.36 for tCr, and 1.27 for tCho. To further elucidate the results, we expanded our SNR calculations to examine J-coupled metabolites—Glu and mI (Wyss et al., 2018)—and found similar SNR gains: a predicted 30% improvement for Glu was experimentally verified, and a 29% improvement was observed for mI, compared to the expected 31%.

FIGURE 5
www.frontiersin.org

Figure 5. Comparison of the 300 mT/m acquisition to the 80 mT/m acquisition. (A) The individual spectra (black) for all b-values are overlaid by the corresponding MRS fits for the 300 mT/m (blue) and 80 mT/m (red), respectively. The approximate b-values are noted on the left, with full details in Table 1. (B) The mono-exponential (dashed coloured line) fits for the 300 mT/m (blue) and 80 mT/m acquisition (red) for tNAA. The vertical line indicates the b-value data exclusion threshold for mono-exponential modeling.

The mean full-width half-maximum of tNAA was found to be 6.2 Hz across both acquisitions, indicating a good shimming and linewidth for LCM quantification. The tNAA CRLBs were all less than 27%—even including the highest b-values—with a mean value of 17% across all b-values. Similarly, the mean CRLB for tCr was 14%, with all-but-one value falling below 22%. CRLBs were higher for tCho with 4 acquisitions greater than 27%. Only a single transient was rejected by the motion corruption metric—for the 300 mT/m acquisition at the highest b-value. Full data quality measures are reported in Supplementary Table S2.

Figure 5B shows the tNAA diffusion decay fitted to a mono-exponential for the single diffusion direction acquired in the OCC. For the mono-exponential representation, the apparent diffusion coefficients of tNAA were ADCe, 300 = 0.168 μm2/ms, and ADCe, 80 = 0.158 μm2/ms for 300 and 80 mT/m, respectively. Results for the other metabolites—tCho and tCr—are reported in Table 2. The Kurtosis representation was not included for these data, as only a single data point fell in the intermediate b-value range of 3, 000–9, 500 s/mm2.

TABLE 2
www.frontiersin.org

Table 2. Listing of the in vivo results applying a mono-exponential and kurtosis representation for the three major metabolites.

Table 3 lists, in the upper rows, the estimated microstructural properties derived from the astrocylinder model comparing the 300 and 80 mT/m settings acquired in the OCC. The estimated AICc's are similar between 80 and 300 mT/m settings. It is interesting to note that D0 is, by average, higher for the Prisma settings, which is also apparent from the signal attenuation at low b-values in Supplementary Figure S3, and in line with the estimated ADCs. Moreover, RC tends toward zero for tNAA for the Prisma configuration.

TABLE 3
www.frontiersin.org

Table 3. Listing of the in vivo results of microstructural measures (free diffusivity D0, cylinder radius RC) estimated from an astrocylinder model.

3.5 Grey matter and white matter

Figure 6A shows the voxel positions and tissue segmentation in the OCC and CR estimated from 3D T1 MPRAGE (Penny et al., 2011; Edden et al., 2014). The OCC voxel contained predominantly grey matter (GM/WM/CSF: 75/17/8 %) and the CR voxel contained predominantly white matter (GM/WM/CSF: 15/83/2 %). Spectral quality was generally high for both voxels. The minimum SNR—as measured by the maximum point of tNAA—was 10.9 in the OCC and 8.9 in the CR at b = 0. The mean FWHM across the DW-MRS acquisitions was 5.65 Hz for the OCC voxel, and 6.31 Hz for the CR voxel, well within consensus recommended limits (Wilson et al., 2019). The tNAA CRLBs ranged from 3.4–9.7% for the OCC, and 5.1–51.9% for the CR, with difference between voxels likely driven by the lower SNR in the CR voxel, or perhaps reflects increased variability due to the proximity of the CR voxel to the ventricles, as previously reported (Genovese et al., 2021a). We also noted lower data quality for the OCC voxel in the y-direction, which exhibited lower SNR and increased CRLBs, compared to the other two directions. Regarding motion-corruption, we excluded a greater number of transients in the data sets for the GM/WM comparison than in the previous analysis. In all cases, no more than 2 transients were removed per diffusion weighting, with one exception for the CR voxel, where 3 transients were removed for a single diffusion condition. All spectra and fits are shown in Supplementary Figure S2. Figure 6B shows the fitting results for tNAA in the OCC and CR voxels. Full data quality measures are reported in Supplementary Table S2.

FIGURE 6
www.frontiersin.org

Figure 6. Results of the in vivo grey and white matter acquisitions. (A) T1-weighted images are overlaid with tissue-specific segmentation masks for grey matter (GM), white matter (WM), and CSF in the occipital (above) and white matter (below) voxels. Grey matter, white matter, and CSF voxel fractions are reported in the corresponding text. (B) The three tNAA diffusion decays for the OCC and CR voxels. Mono-exponential (dashed lines) and kurtosis (solid lines) fits are shown for the x (blue), y (green), and z (red) directions. The vertical lines indicate the b-value thresholds of data exclusion for mono-exponential and kurtosis modelling, respectively.

Metabolite ADCe for the mono-exponential representation and ADCK and K for the kurtosis representation are reported in Table 2 in the lower rows for each individual diffusion-encoding direction. The mean directional-averaged ADCes for tNAA/tCho/tCr are 0.14/0.17/0.18μm2ms−1 in the CR and 0.13/0.17/0.18μm2ms−1 in the OCC. For ADCK slightly higher values were found with 0.16/0.17/0.18μm2ms−1 in the CR and 0.17/0.14/0.20μm2ms−1 in the OCC, and corresponding Kurtosis values of 1.90/1.14/1.33 in the CR and 1.75/1.66/1.39 in the OCC.

The estimated microstructural properties for the astrocylinder model are listed in Table 3 in the lower rows for the voxels placed in the WM-rich CR and GM-rich OCC acquired with Connectom settings. The AICc's are comparable between these two brain regions and do not indicate greater model validity in either WM or GM.

3.6 Table vibration

Figure 2C shows the signal decay for the isotropic diffusion phantom acquired along the three orthogonal directions x, y, and z for the Connectom and Prisma settings. The ADC for both settings only shows a slightly faster diffusion along the z-direction, which may point to stronger table vibrations along z. However, the overall estimated ADCs are well in line with the calibration value of 0.293 μm2/ms at 22°C. Considering our results in vivo shown in Figure 6 we found a stronger signal decay for tNAA, but also tCr and tCho (not shown), in the CR when diffusion-encoding is applied along the z-direction at ultra-high b-values.

4 Discussion

In this study, we investigated the practical feasibility of DW-MRS on a Connectom MR scanner equipped with 300 mT/m gradients to measure metabolite diffusion at ultra-high b-values in the human brain. By comparing these results between clinical and ultra-high gradient amplitudes, we show the benefits of new gradient systems for DW-MRS, but also present strategies to mitigate the impact of methodological challenges.

4.1 SNR simulations

Simulated SNR benefits of applying ultra-strong gradients were mostly in line with experimental measurements in vivo. Comparing the Prisma (Gmax = 80 mT/m) and Connectom (Gmax = 300 mT/m) settings, the theoretical SNR gain of 10% for tNAA was even exceeded experimentally and a 16% improvement was achieved. A similar pattern is observed for tCho, where experimental SNR gains exceed predictions by 7 percentage points. J-coupled metabolites mI and Glu were approximately in-line with expectations. For tCr, while we still report an increase in the SNR at 300 mT/m, the magnitude of the SNR increase is not as large as predicted by simulation—36% vs. the expected 50%. This variation from expected SNR gains is attributed to inaccuracies in the assumed T2, but further verification would require acquisition of T2 relaxometry alongside the DW-MRS. Interestingly, two studies that investigated this—in humans and mice respectively—found that tCr was the only metabolite that showed a slight dependence of ADC on TE (?Mougel et al., 2022).

In general, SNR should be carefully considered in the context of diffusion imaging (Ellegood et al., 2011). The SNR dependence of MRS reconstruction and fitting methods can lead to anomalous signal decay, particularly at higher b-values, making it an important consideration for high-b-value DW-MRS. The shorter TE afforded by ultra-strong gradients improved SNR and reduced CRLBs, but strict consideration of SNR and MRS fit uncertainty is still necessary.

4.2 In vivo results

Previous work (Ronen et al., 2014) found tNAA diffusivity values of 0.076 and 0.34 μm2/ms for diffusion gradients orthogonal and parallel to the main fibre orientation of the corpus callosum, respectively. The DW-MRS voxels considered in this study were arranged to contain predominantly grey or white matter, but inevitably contained mixed fibre orientations. Our measured ADCs fell within the range reported by Ronen et al. (2014), with our lowest ADC reported in grey matter (0.096 μm2/ms) and the highest in kurtosis model of white matter (0.194 μm2/ms).

Furthermore, higher ADCs in GM compared to WM have been previously reported (Ellegood et al., 2011; Kan et al., 2012; Najac et al., 2016) and this finding is confirmed by our results for tNAA in our grey-white matter comparison, where we find a mean ADC of 0.126 μm2/ms for the OCC GM voxel, and 0.140 μm2/ms for the CR WM voxel. Kurtosis value ranges for tNAA (1.5–2.2), tCr (1.0–1.6), and—excluding one low-kurtosis fit (0.192)—tCho (1.5–1.9) fall within previously reported diffusion-time-dependent ranges for metabolites (Döring et al., 2023). Lower AICc values were found in the data averaged over multiple diffusion directions. This difference could potentially be attributed to a higher SNR or better agreement with the model, resulting from powder averaging. Furthermore, we examined the data to assess the feasibility of employing more extensive modeling strategies, specifically using a two-compartment model as described in Supplementary Figure S4 and Supplementary Table S1. However, the estimated AICc values were consistently were higher.

We report a 6% decrease in the tNAA ADC of 80 mT/m acquisition, compared to the 300 mT/m acquisition. Conversely, a significant increase was noted for the other two metabolites (tCho, tCr), in line with an even more pronounced increase in the free diffusivity D0. Although the goodness-of-fit was lower for the 80 mT/m configuration, a possible explanation could be that the longer TE reduces signal contributions from restricted compartments—such as organelles—and consequently increases ADCs and D0s. Previous studies within a TE range of 35 to 70 ms have not found a correlation between metabolite diffusion and T2 relaxation, but others at 7T and longer TE report on faster diffusion with increasing TE (?Ligneul et al., 2023). This is in line with our results for tCho and tCr indicating faster diffusion when TE increases from 74 to 116 ms, although T2s were different due to the lower field-strength and diffusion-times due to different maximum gradient amplitudes. Moreover, despite the higher SNR, shorter TE, and shorter diffusion-times achieved with the 300 mT/m settings, AICc's indicate similar model support.

It is important to note that the relatively small size of our volunteer group limits the confidence with which we can draw conclusions about our quantitative analyses. Uncertainty in the MRS modeling procedures—due to a myriad of factors, including lineshape distortion, motion, uncharacterized signals from macromolecules, and issues in baseline characterisation—will affect the measured diffusion properties necessitating a larger participant population to achieve sufficient statistical power. Here, we utilize the Cramer-Rao lower bounds, and weight the fitting accordingly in an attempt to mitigate this. However, in the future it might be beneficial to use metabolite-cycling to further reduce motion and eddy-current artifacts (Döring et al., 2018), and expanding the sample size of the study.

4.3 Eddy currents

Diffusion-weighted imaging and spectroscopy use strong gradients to achieve the desired diffusion weighting. Eddy currents generated during ramp up and down times can cause time-dependent frequency variations in the time domain and consequently distort the MR spectra. Correcting for eddy-current effects is vital in MRS—and particularly so for DW-MRS—in order to accurately reflect lineshape distortions while modeling MR spectra. Furthermore, as per Faraday's law, the magnitude of the induced eddy-current effects increase with the applied gradient strength, further compounding their relevance in the context of ultra-strong gradient DW-MRS. Thankfully, in the particular case of the Siemens Connectom scanner, the shield coil design minimizes the relative size of eddy current effects. Setsompop et al. (2013) reported that a 7.5-fold increase in gradient amplitude on the Connectom scanner corresponded to just a 2–3-fold eddy current increase, in absolute terms. Our sequence, acquisition scheme, and post processing pipeline further reduce the impact of eddy currents. Our DW-PRESS sequence utilized a bipolar diffusion gradient scheme, which has the potential to reduce eddy current effects on the acquired DW-MR spectra (Alexander et al., 1997; Branzoli, 2015). We further applied two methods of post-hoc eddy current correction.

Our in vivo DW-MRS protocols included unsuppressed water-signal, which we used to correct for eddy-current-induced phase and frequency shifts (Klose, 1990). Though this method requires acquisition of only few additional spectra, faster water diffusion may have reduced the effectiveness of this method at the highest b-values applied in this work. Therefore, we also demonstrate a second approach (Lin et al., 1994). For this, we inverted the polarity of diffusion-weighting gradients for half of metabolite spectra acquired at the highest b-values. The inversion of diffusion gradient polarities inverts the phase of the generated eddy currents, while the intrinsic signal phase and the magnitude of the diffusion-weighting are expected to remain unchanged. No additional acquisitions are required in this case. However, it should be noted that motion corruption of individual transients would necessitate the pairwise removal of the corresponding inverted diffusion condition, potentially doubling motion-related transient exclusions. Furthermore, our implementation of the method did not include the inversion of slice-selective gradients, which may explain the residual asymmetry of the lineshapes in our combined spectra.

The quality of eddy current correction of the above methods relies on the quality of the acquired water and metabolite spectra. A lot of work has been done to monitor the field perturbations in presence of diffusion gradients and their effect of the readout trajectory and the consequent deterioration in image quality in DW-MRI (Chan et al., 2014; Wilm et al., 2015). External probes could be a valuable tool in characterizing and correcting for the eddy current effects present in the FID. The existing DW-PRESS sequence would require a trigger event prior the FID acquisition and a careful temporal and spatial alignment between the FID from the spectroscopy voxel and the estimated eddy current effects from FIDs measured by the probes.

4.4 Gradient non-uniformities

Gradient coils are often designed with a limited field of view (FOV) to minimize nerve stimulation caused by rapidly changing magnetic fields over time. However, this design choice will increase gradient spatial non-uniformity, which becomes more noticeable as one moves away from the isocenter. In high-performance gradient systems like the one used here, gradient uniformity is often further compromized to achieve better performance. There is increased awareness of the importance of correcting for gradient deviations in diffusion MRI studies (Bammer et al., 2003; Guo et al., 2020; Mesri et al., 2020; Morez et al., 2021), but the specification of a gradient system is commercially sensitive information and not widely available.

Two strategies to represent the effective b-value were tested, i.e., incorporating the distribution of b-values, and representing that distribution using the median. Both strategies show similar deviations of the estimated ADC in DW-MRS, but both varied from the uncorrected nominal b-value, it is thus important to take into account the effective b-values when this information is available. Also, at lower gradient strength, and lower b-values—where the absolute signal change as a function of b-value is typically largest—gradient non-uniformities can bias quantitative estimates. In addition to b-value deviations, gradient non-uniformities will also result in distortions of the slice profile and thus voxel geometry. While the distortions we observed were small, it should be a consideration when placing voxels close to tissue boundaries, and conservative placement within the relevant tissue is advised, particularly when deviating far from the isocenter, where gradient non-uniformities are larger. Furthermore, there will likely be a small impact on chemical shift displacement error (CSDE). While CSDE is typically a linear chemical-shift-dependent translation of the effective voxel, when gradient non-uniformities are substantial, this effect is no longer a simple translation, and will result in metabolite-specific voxel deformation. This is an area that requires further study; however, we expect this effect to be minimal due to the lower amplitude of the imaging gradients, and the linear—rather than quadratic—effect that non-uniformities have on slice-selective gradients.

4.5 Table vibrations

The switching of strong diffusion gradients can cause mechanical vibrations at low frequencies within the scanner system and table. These vibrations could then be transmitted to the subject being scanned, resulting in artificial signal decay (Hiltunen et al., 2006) which artificially increases the measured ADC. While the Connectom scanner does, indeed, provide access to stronger gradients, it is not trivial how this may manifest as mechanical vibration. Interestingly, Setsompop et al. (2013) reported that the effects of acoustic noise generated by a 300 mT/m Connectom scanner were, in fact, lower than those at more conventional gradient strengths (40 mT/m). They postulated that the thicker, larger, and heavier gradient system of the Connectom scanner somewhat offset the larger amplitude of vibrational forces it generated (Setsompop et al., 2013).

Although the relatively lower directional variance of ADCs measured in the NIST phantom at the Gmax = 300 mT/m setting was surprising, it is important to note that in vivo tissue stiffness and composition is different and likely more prone to vibration-induced artifacts. Indeed, our measurements in the CR show a particularly strong signal decay when diffusion-encoding is applied along z-direction. This might originate from tissue anisotropy and major white-matter tracts pointing along z (e.g., cortico-spinal or cortico-pontine tract) when diffusion of tNAA is considered. However, tCho shows consistent signal dropouts and diffusion in the low b-value range remains unaffected, which could point to effects related to table vibrations at high gradient amplitudes. Our hypothesis is that the CR—further away from the contact point between the head and the coil—is more strongly affected by random rotations induced by table vibrations than the OCC. This would be in line with previous observations of largest vibrations on the Connectom when diffusion-encoding is applied along the z-direction (Mueller et al., 2019). However, additional dedicated investigations with different sequence parameters (voxel position, diffusion-time) are required to elucidate this further.

In general, this artefact will be hardware- and sequence-specific but can be mitigated at the hardware level through careful consideration of the coil/gradient mounting (Ogura et al., 2006; Mueller et al., 2019), or even participant positioning, as different placement of padding around the head can lead to differing vibrational coupling between the brain and RF coil housing (Gallichan et al., 2010).

4.6 Other considerations

4.6.1 Macromolecules

The macromolecular background is a potential confound to DW-MRS quantification (Ronen and Valette, 2007), with slowly-diffusing signal components persistent throughout the diffusion conditions. The impact of macromolecules is mitigated in this study as we acquired data at a longer TE, and parameterized the MM background during modeling. Experimentally acquired MM backgrounds are preferable—and can be acquired with ultra-strong gradient DW-MRS (Şimşek et al., 2022)—but are parameterization dependent, and are perhaps best-acquired at the cohort level in larger studies (Zöllner et al., 2023).

4.6.2 Other sequences

While our study focused on a DW-PRESS implementation, the methodological considerations we identified are transferable to other DW-MRS localisation methods. The gradient non-uniformities can be corrected using the same principles, with simple modifications made to account for the specific gradient profile. Similarly, both of the eddy current correction methods considered here are also applicable to other single-shot localization schemes, and DW-PRESS might actually present a particularly challenging example. PRESS localization is known to be more susceptible to chemical shift displacement error than STEAM or sLASER, and this—coupled with the strong spatiotemporal dependence of eddy currents—can adversely affect the efficacy of the eddy-current correction in a metabolite-specific manner. The T2-based SNR gains provided by ultra-strong gradients are also transferable to other sequences with some caveats. Both DW-STEAM and DW-sLASER will benefit from the T2-based SNR gains, but are less sensitive to anomalous J-modulation effects than PRESS; however, our data suggest that this has a minimal effect on our DW-PRESS data. DW-STEAM also benefits from the decoupling of diffusion time from the TE, and this can be used to circumvent some of the T2 limitations to achieve high b-values. However, when short diffusion times are required—for example, to probe shorter length scales and/or minimize motion artefacts—ultra-strong gradients can still provide benefits. Furthermore, when measuring metabolite diffusion in small anatomical structures such as thalamus, hippocampus, brainstem, or spinal cord, higher gradient amplitudes combined with DW-STEAM or DW-sLASER localization can enhance an accurate localization (with a minimal CSDE) and should be prefered over DW-PRESS 2.

Another practical consideration that we didn't directly investigate in this work is the effect of concomitant gradient fields. These are nonlinear transverse components to an applied gradient field that arise from Maxwell's equations, and are typically more prominent at lower static field strengths and higher gradient field strengths (Baron et al., 2012). Our DW-PRESS sequence implementation utilized symmetric gradient timing to mitigate the impact of this effect, but for non-symmetric gradient implementations—for example, free gradient waveforms—this effect becomes an important consideration that should be corrected in a prospective manner (Baron et al., 2012; Szczepankiewicz et al., 2019).

4.6.3 Other preprocessing steps

In this study, a retrospective outlier rejection method was used to identify and remove suspected motion-corrupted transients by comparing subsequent transients to the first. While 3 transients were removed in one such case—at the highest b-value—this did not impact the ADC estimations using the mono-exponential and kurtosis representations, which were performed using lower b-values. However, this may be adversely affect the SNR of highest b-value acquisitions in the other models. Weighting using the CRLBs mitigated the impact this had on model performance. Furthermore, if the initial transient of a given series is corrupted by motion, a high rate of outlier rejection would be observed. While this case wasn't encountered in these data, a secondary process—perhaps comparing to the median spectrum—would circumvent this issue. Our approach was conservative, and an alternative approach could be to monitor bulk motion in a prospective manner (Andrews-Shigaki et al., 2011; Saleh et al., 2020).

We have previously investigated the impact of phase and frequency correction methods on DW-MRS data (Jenkins, 2021), and the inherent SNR-dependent performance of such methods. Signal denoising techniques could mitigate this somewhat, and show promise for DW-MRS data (Mosso et al., 2022); however, care must be taken not to invalidate assumptions about noise characteristics during further modeling steps (Dziadosz et al., 2023).

4.6.4 Modeling

While we focused on the practical considerations here, alternative DW-MRS modeling procedures could be beneficial, particularly leveraging high b-values and/or the higher SNR to further disentangle compartments (i.e., cylinders and spheres, Supplementary Figure S4, Supplementary Table S1). With the current data, sufficient support was not found for a more complex model, but future work could investigate this with specifically-acquired data.

Multi-spectrum modeling of MRS data—fitting the frequency and b-value dimensions in a single step, rather than independently—is another promising avenue which improves model parsimony and the stability of diffusion measures to experimental variation (Tal, 2023), and notably, allows stable parameterization of the macromolecules and baseline. Several groups have worked to develop this methodology in recent years e.g., FitAid (Chong et al., 2011; Adalid et al., 2017), FSL-MRS (Clarke et al., 2021), and Osprey (Oeltzschner et al., 2020).

4.7 Limitations

We were able to demonstrate the SNR improvement of incorporating ultra-strong gradients. However, cardiac and nerve stimulation limits place firm lower bounds on gradient ramp time and restrict the achievable gradient amplitudes at shorter TE. As a result, at TEs below 70 ms, the maximum gradient amplitude of the Connectom scanner is not achievable within the echo time, at least with the DW-PRESS sequence and diffusion gradient configuration we used. This limits the potential benefit provided by the ultra-strong gradient system. However, in practice, the ultra-high b-values we achieved in this study might not be necessary for all applications, and sacrificing diffusion weighting in favour of the improved SNR brought by lower TE might be preferable. This would be particularly beneficial when smaller voxel sizes or shorter scan times are a necessity. Moreover, for DW-PRESS, shorter TE also reduces the diffusion time, probing smaller length scales and different aspects of cell morphology. In practice, lower b-values and shorter TE might be preferable for detecting cell morphology changes in order to reduce scan time and/or voxel size.

Our relatively small cohort of volunteers limits the interpretability of the reported ADCs. A larger study is required to elucidate the effects of diffusion time, and systematically validate microstructure across tissue types. Moreover, while our multiple-gradient-condition acquisition allowed reasonable comparisons with conventional gradient systems, a systematic repeatability study using independent gradient systems would be necessary to fully vindicate these results. Indeed, the Connectom gradient design differs from that of a lower-gradient system, affecting gradient non-uniformity even at low gradient amplitudes, and cross-scanner comparisons would automatically include such differences (Gudino and Littin, 2023).

Furthermore, while our focus in this work was in the implementation of DW-MRS with ultra-strong gradients, and circumvention of major artefacts, we primarily limited our analysis to the major metabolites.

5 Conclusion

We successfully implemented acquisition and data processing strategies for ultra-strong gradient DW-MRS. We report diffusion coefficients which conform with the existing literature. Simulated SNR gains are experimentally confirmed, and results indicate that confounding effects of the strong gradient system can be ameliorated.

Data availability statement

The datasets presented in this article are not readily available because we do not have ethical consent to make the in-vivo datasets acquired for this study publicly available. Requests to access the datasets should be directed to TaxC@cardiff.ac.uk.

Ethics statement

The studies involving humans were approved by Cardiff University School of Psychology. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

CD-J: Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing. AD: Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing. FF: Methodology, Software, Writing – review & editing. EK: Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing. LM: Investigation, Methodology, Software, Writing – review & editing. CE: Conceptualization, Methodology, Supervision, Writing – review & editing. MA: Software, Writing – review & editing. DJ: Funding acquisition, Resources, Supervision, Writing – review & editing. IR: Conceptualization, Methodology, Software, Supervision, Writing – review & editing. FB: Conceptualization, Methodology, Software, Supervision, Writing – review & editing. CT: Conceptualization, Formal analysis, Investigation, Methodology, Software, Supervision, Visualization, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was funded in whole, or in part, by a Wellcome Trust Investigator Award (096646/Z/11/Z) and a Wellcome Trust Strategic Award (104943/Z/14/Z). CT is supported by the Wellcome Trust [215944/Z/19/Z] and a Veni grant (17331) from the Dutch Research Council (NWO). AD is supported by a Swiss National Science Foundation Fellowship (SNSF #202962). FB acknowledges support from Investissements d'avenir [grant Nos. ANR-10-IAIHU-06 and ANR-11-INBS-0006]. MA and LM are supported by a Wellcome Trust Investigator Award (219536/Z/19/Z). For the purpose of open access, the author has applied a CC BY public copyright licence to any Author Accepted Manuscript version arising from this submission.

Acknowledgments

The authors would like to acknowledge the contribution of Umesh Rudrapatna whose efforts aided in the development of the DW-MRS sequence and data processing strategies.

Conflict of interest

FF was employed by company Siemens Healthcare Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2023.1258408/full#supplementary-material

Supplementary Figure S1. Figure showing the MRS fitting residuals for the 300 mT/m (left) and 80 mT/m (right) acquisitions. Residuals are plotted in red, ascending in b-value from top to bottom.

Supplementary Figure S2. This figure shows the individual processed spectra (black) and fits (red) for the OCC (left) and CR (right) voxels for all b-values. Individual diffusion directions are grouped by diffusion value, and are ordered from top to bottom: x, y, and z, respectively.

Supplementary Figure S3. Eddy current phase evolution from three different gradient amplitudes: 295 mT/m (blue), 150 mT/m (red), and 50 mT/m (yellow).

Supplementary Figure S4. Fitting results applying an astrocylinder model and the estimated microstructural measures (free diffusivity D0, cylinder radius RC). Outliers were identified by the RMSE when iteratively removing a single data point from fitting.

Supplementary Figure S5. Fitting results applying a two-compartment model and the estimated microstructural measures (free diffusivity D0, sphere radius RS, cylinder radius RC). The fraction of cylinders fC was kept fixed at 0.8 and the standard deviation was estimated from bootstrapping. Outliers were identified by the RMSE when iteratively removing a single data point from fitting.

Supplementary Table S1. Listing of the in vivo results of microstructural measures (free diffusivity D0, sphere radius RS, cylinder radius RC) estimated from a two-compartment model. The double lined column delineation shows in the upper half the comparison of the 300 and 80 mT/m settings where diffusion-encoding was applied only along the z-direction, and in the lower half the comparison of the averaged diffusion metrics over three orthogonal diffusion-directions in the white matter rich corona radiata (CR) and grey matter rich occipital lobe (OCC) using the 300 mT/m setting. The fitting results can be found in Supplementary Figure S4. *Fraction of cylinders fC was kept fixed at 0.8, c.f., text.

Supplementary Table S2. Full data quality information for the in vivo data. Each section contains entries for the Cramer-Rao lower bound (CRLB) of tNAA, tCho, and tCr, as well as the full width half maximum (FWHM), signal-to-noise ratio (SNR), points excluded by the motion-corruption metric (Excl.), and the mean frequency (Freq-Mean) and phase (Phs-Mean) shifts applied by spectral registration. Entries are ordered by increasing b-value per section.

Footnotes

1. ^Full code available https://github.com/CWDAVIESJENKINS/FIDA-Tarquin-Processing/tree/main/PipeLines/DWMRS.

2. ^At ultra-high field (≥7T) PRESS localization is severely compromised by CSD, and STEAM or sLASER localization is recommended.

References

Adalid, V., Döring, A., Kyathanahally, S. P., Bolliger, C. S., Boesch, C., and Kreis, R. (2017). Fitting interrelated datasets: metabolite diffusion and general lineshapes. Magn. Reson. Mater. Phys. Biol. Med. 30, 429–448. doi: 10.1007/s10334-017-0618-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Alexander, A. L., Tsuruda, J. S., and Parker, D. L. (1997). Elimination of eddy current artifacts in diffusion-weighted echo-planar images: the use of bipolar gradients. Magn. Reson. Med. 38, 1016–1021. doi: 10.1002/mrm.1910380623

PubMed Abstract | CrossRef Full Text | Google Scholar

Andrews-Shigaki, B. C., Armstrong, B. S., Zaitsev, M., and Ernst, T. (2011). Prospective motion correction for magnetic resonance spectroscopy using single camera retro-grate reflector optical tracking. J. Magn. Reson. Imag. 33, 498–504. doi: 10.1002/jmri.22467

PubMed Abstract | CrossRef Full Text | Google Scholar

Bammer, R., Markl, M., Barnett, A., Acar, B., Alley, M., Pelc, N., et al. (2003). Analysis and generalized correction of the effect of spatial gradient field distortions in diffusion-weighted imaging. Magn. Reson. Med. 50, 560–569. doi: 10.1002/mrm.10545

PubMed Abstract | CrossRef Full Text | Google Scholar

Baron, C., Lebel, R., Wilman, A., and Beaulieu, C. (2012). The effect of concomitant gradient fields on diffusion tensor imaging. Magn. Reson. Med. 68, 1190–1201. doi: 10.1002/mrm.24120

PubMed Abstract | CrossRef Full Text | Google Scholar

Bottomley, P. A. (1987). Spatial localization in NMR spectroscopy in vivo. Ann. NY. Acad. Sci. 508, 333–348. doi: 10.1111/j.1749-6632.1987.tb32915.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Branzoli, F. (2015).

Google Scholar

Chan, R. W., von Deuster, C., Giese, D., Stoeck, C. T., Harmer, J., Aitken, A. P., et al. (2014). Characterization and correction of eddy-current artifacts in unipolar and bipolar diffusion sequences using magnetic field monitoring. J. Magn. Reson. 244, 74–84. doi: 10.1016/j.jmr.2014.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, J.-K., Dedeoglu, A., and Jenkins, B. G. (2007). Application of MRS to mouse models of neurodegenerative illness. NMR Biomed. 20, 216–237. doi: 10.1002/nbm.1145

CrossRef Full Text | Google Scholar

Chong, D. G., Kreis, R., Bolliger, C. S., Boesch, C., and Slotboom, J. (2011). Two-dimensional linear-combination model fitting of magnetic resonance spectra to define the macromolecule baseline using fitaid, a fitting tool for arrays of interrelated datasets. Magn. Reson. Mater. Phys. Biol. Med. 24, 147–164. doi: 10.1007/s10334-011-0246-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Clarke, W. T., Stagg, C. J., and Jbabdi, S. (2021). FSL-MRS: an end-to-end spectroscopy analysis package. Magn. Reson. Med. 85, 2950–2964. doi: 10.1002/mrm.28630

PubMed Abstract | CrossRef Full Text | Google Scholar

Döring, A., Adalid, V., Boesch, C., and Kreis, R. (2018). Diffusion-weighted magnetic resonance spectroscopy boosted by simultaneously acquired water reference signals. Magn. Reson. Med. 80, 2326–2338. doi: 10.1002/mrm.27222

PubMed Abstract | CrossRef Full Text | Google Scholar

Döring, A., Rösler, F., Şimşek, K., Afzali, M., Kreis, R., Jones, D. K., et al. (2023). “Time dependent diffusion and kurtosis of human brain metabolites, in Proceedings of the Virtual Conference of ISMRM.

Google Scholar

Dziadosz, M., Rizzo, R., Kyathanahally, S. P., and Kreis, R. (2023). Denoising single MR spectra by deep learning: miracle or mirage? Magn. Reson. Med. 90, 1749–1761. doi: 10.1002/mrm.29762

PubMed Abstract | CrossRef Full Text | Google Scholar

Edden, R. A., Puts, N. A., Harris, A. D., Barker, P. B., and Evans, C. J. (2014). Gannet: a batch-processing tool for the quantitative analysis of gamma-aminobutyric acid-edited mr spectroscopy spectra. J. Magn. Reson. Imag. 40, 1445–1452. doi: 10.1002/jmri.24478

PubMed Abstract | CrossRef Full Text | Google Scholar

Ellegood, J., Hanstock, C. C., and Beaulieu, C. (2011). Considerations for measuring the fractional anisotropy of metabolites with diffusion tensor spectroscopy. NMR Biomed. 24, 270–280. doi: 10.1002/nbm.1586

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, Q., Eichner, C., Afzali, M., Mueller, L., Tax, C. M., Davids, M., et al. (2022). Mapping the human connectome using diffusion MRI at 300 mt/m gradient strength: methodological advances and scientific impact. NeuroImage 254, 118958. doi: 10.1016/j.neuroimage.2022.118958

PubMed Abstract | CrossRef Full Text | Google Scholar

Gallichan, D., Scholz, J., Bartsch, A., Behrens, T. E., Robson, M. D., and Miller, K. L. (2010). Addressing a systematic vibration artifact in diffusion-weighted MRI. Hum. Brain Mapp. 31, 193–202. doi: 10.1002/hbm.20856

PubMed Abstract | CrossRef Full Text | Google Scholar

Ganji, S. K., Banerjee, A., Patel, A. M., Zhao, Y. D., Dimitrov, I. E., Browning, J. D., et al. (2012). T2 measurement of J-coupled metabolites in the human brain at 3 T. NMR Biomed. 25, 523–529. doi: 10.1002/nbm.1767

PubMed Abstract | CrossRef Full Text | Google Scholar

Genovese, G., Marjańska, M., Auerbach, E. J., Cherif, L. Y., Ronen, I., Lehéricy, S., et al. (2021a). In vivo diffusion-weighted mrs using semi-laser in the human brain at 3 T: Methodological aspects and clinical feasibility. NMR Biomed. 34, e4206. doi: 10.1002/nbm.4206

PubMed Abstract | CrossRef Full Text | Google Scholar

Genovese, G., Palombo, M., Santin, M. D., Valette, J., Ligneul, C., Aigrot, M.-S., et al. (2021b). Inflammation-driven glial alterations in the cuprizone mouse model probed with diffusion-weighted magnetic resonance spectroscopy at 11.7 T. NMR Biomed. 34, e4480. doi: 10.1002/nbm.4480

PubMed Abstract | CrossRef Full Text | Google Scholar

Goveas, J., O'Dwyer, L., Mascalchi, M., Cosottini, M., Diciotti, S., De Santis, S., et al. (2015). Diffusion-MRI in neurodegenerative disorders. Magn. Reson. Imag. 33, 853–876. doi: 10.1016/j.mri.2015.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Gudino, N., and Littin, S. (2023). Advancements in gradient system performance for clinical and research MRI. J. Magn. Reson. Imag. 57, 57–70. doi: 10.1002/jmri.28421

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, F., de Luca, A., Parker, G., Jones, D., Viergever, M., Leemans, A., et al. (2020). The effect of gradient nonlinearities on fiber orientation estimates from spherical deconvolution of diffusion magnetic resonance imaging data. Hum. Brain Mapp. 42, 367–383. doi: 10.1002/hbm.25228

PubMed Abstract | CrossRef Full Text | Google Scholar

Hall, E. L., Stephenson, M. C., Price, D., and Morris, P. G. (2014). Methodology for improved detection of low concentration metabolites in MRS: optimised combination of signals from multi-element coil arrays. Neuroimage 86, 35–42. doi: 10.1016/j.neuroimage.2013.04.077

PubMed Abstract | CrossRef Full Text | Google Scholar

Hiltunen, J., Hari, R., Jousmki, V., Muller, K., Sepponen, R., and Joensuu, R. (2006). Quantification of mechanical vibration during diffusion tensor imaging at 3T. NeuroImage 32, 93–103. doi: 10.1016/j.neuroimage.2006.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, S. Y., Witzel, T., Keil, B., Scholz, A., Davids, M., Dietz, P., et al. (2021). Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso-and macro-connectome. Neuroimage 243, 118530. doi: 10.1016/j.neuroimage.2021.118530

PubMed Abstract | CrossRef Full Text | Google Scholar

Hutchinson, E. B., Schwerin, S. C., Avram, A. V., Juliano, S. L., and Pierpaoli, C. (2018). Diffusion MRI and the detection of alterations following traumatic brain injury. J. Neurosci. Res. 96, 612–625. doi: 10.1002/jnr.24065

PubMed Abstract | CrossRef Full Text | Google Scholar

Jelescu, I. O., de Skowronski, A., Geffroy, F., Palombo, M., and Novikov, D. S. (2022). Neurite exchange imaging (NEXI): a minimal model of diffusion in gray matter with inter-compartment water exchange. NeuroImage 256, 119277. doi: 10.1016/j.neuroimage.2022.119277

PubMed Abstract | CrossRef Full Text | Google Scholar

Jenkins, C., Kleban, E., Mueller, L., Evans, C. J., Rudrapanta, U., Jones, D., et al. (2020). “Dw-MRS with ultra-strong diffusion gradients,” in ISMRM.

Google Scholar

Jenkins, C. W. (2021). “The influence of spectral registration on diffusion-weighted magnetic resonance spectroscopy ADC estimates,” in ISMRM.

Google Scholar

Jensen, J. H., Helpern, J. A., Ramani, A., Lu, H., and Kaczynski, K. (2005). Diffusional kurtosis imaging: the quantification of non-gaussian water diffusion by means of magnetic resonance imaging. Magn. Reson. Med. 53, 1432–1440. doi: 10.1002/mrm.20508

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, D. K., Alexander, D. C., Bowtell, R., Cercignani, M., Dell'Acqua, F., McHugh, D. J., et al. (2018). Microstructural imaging of the human brain with a super-scanner: 10 key advantages of ultra-strong gradients for diffusion mri. NeuroImage 182, 8–38. doi: 10.1016/j.neuroimage.2018.05.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Kan, H. E., Techawiboonwong, A., Van Osch, M. J., Versluis, M. J., Deelchand, D. K., Henry, P.-G., et al. (2012). Differences in apparent diffusion coefficients of brain metabolites between grey and white matter in the human brain measured at 7 t. Magn. Reson. Med. 67, 1203–1209. doi: 10.1002/mrm.23129

PubMed Abstract | CrossRef Full Text | Google Scholar

Ke, Y., Cohen, B., Lowen, S., Hirashima, F., Nassar, L., and Renshaw, P. (2002). Biexponential transverse relaxation (t2) of the proton mrs creatine resonance in human brain. Magn. Reson. Med. 47, 232–238. doi: 10.1002/mrm.10063

PubMed Abstract | CrossRef Full Text | Google Scholar

Klose, U. (1990). In vivo proton spectroscopy in presence of eddy currents. Magn. Reson. Med. 14, 26–30. doi: 10.1002/mrm.1910140104

PubMed Abstract | CrossRef Full Text | Google Scholar

Landheer, K., Gajdošík, M., and Juchem, C. (2020). A semi LASER, single-voxel spectroscopic sequence with a minimal echo time of 20.1 ms in the human brain at 3 T. NMR Biomed. 33, 1–12. doi: 10.1002/nbm.4324

PubMed Abstract | CrossRef Full Text | Google Scholar

Ligneul, C., Najac, C., Döring, A., Beaulieu, C., Branzoli, F., Clarke, W. T., et al. (2023). Diffusion-weighted MR spectroscopy: consensus, recommendations and resources from acquisition to modelling. arXiv preprint arXiv:2305.10829.

PubMed Abstract | Google Scholar

Ligneul, C., Palombo, M., Hernández-Garzón, E., Carrillo-de Sauvage, M.-A., Flament, J., Hantraye, P., et al. (2019). Diffusion-weighted magnetic resonance spectroscopy enables cell-specific monitoring of astrocyte reactivity in vivo. NeuroImage 191, 457–469. doi: 10.1016/j.neuroimage.2019.02.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, C., Wendt, R. E., Evans, H. J., Rowe, R. M., Hedrick, T. D., and LeBlanc, A. D. (1994). Eddy current correction in volume-localized MR spectroscopy. J. Magn. Reson. Imaging 4, 823–827. doi: 10.1002/jmri.1880040614

PubMed Abstract | CrossRef Full Text | Google Scholar

Mattiello, J., Basser, P. J., and Le Bihan, D. (1997). The B matrix in diffusion tensor echo-planar imaging. Magn. Reson. Med. 37, 292–300. doi: 10.1002/mrm.1910370226

PubMed Abstract | CrossRef Full Text | Google Scholar

Mesri, H. Y., David, S., Viergever, M. A., and Leemans, A. (2020). The adverse effect of gradient nonlinearities on diffusion MRI: from voxels to group studies. NeuroImage 205, 116127. doi: 10.1016/j.neuroimage.2019.116127

PubMed Abstract | CrossRef Full Text | Google Scholar

Morez, J., Sijbers, J., Vanhevel, F., and Jeurissen, B. (2021). Constrained spherical deconvolution of nonspherically sampled diffusion MRI data. Hum. Brain Mapp. 42, 521–538. doi: 10.1002/hbm.25241

PubMed Abstract | CrossRef Full Text | Google Scholar

Mosso, J., Simicic, D., Şimşek, K., Kreis, R., Cudalbu, C., and Jelescu, I. O. (2022). MP-PCA denoising for diffusion mrs data: promises and pitfalls. NeuroImage 263, 119634. doi: 10.1016/j.neuroimage.2022.119634

PubMed Abstract | CrossRef Full Text | Google Scholar

Mougel, E., Malaquin, S., and Valette, J. (2022). Assessing potential correlation between t2 relaxation and diffusion of lactate in the mouse brain. Magn. Reson. Med. 88, 2277–2284. doi: 10.1002/mrm.29395

PubMed Abstract | CrossRef Full Text | Google Scholar

Mueller, L., Rudrapatna, S. U., Tax, C. M. W., Wise, R. G., and Jones, D. K. (2019). “Diffusion MRI with b = 1000 s/mm2 at te < 22 ms using single-shot spiral readout and ultra-strong gradients: Implications for microstructure imaging, in Proceedings of the 27th Annual Meeting of ISMRM.

Google Scholar

Mugler III, J. P., and Brookeman, J. R. (1990). Three-dimensional magnetization-prepared rapid gradient-echo imaging (3D mp rage). Magn. Reson. Med. 15, 152–157. doi: 10.1002/mrm.1910150117

PubMed Abstract | CrossRef Full Text | Google Scholar

Najac, C., Branzoli, F., Ronen, I., and Valette, J. (2016). Brain intracellular metabolites are freely diffusing along cell fibers in grey and white matter, as measured by diffusion-weighted mr spectroscopy in the human brain at 7 T. Brain Struct. Funct. 221, 1245–1254. doi: 10.1007/s00429-014-0968-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Najac, C., Döring, A., Clarke, W., Genovese, G., Just, N., Kreis, R., et al. (2022). “Differences in diffusion-weighted mrs processing and fitting pipelines, and their effect on tissue modeling: results from a workshop challenge, in Proceedings of the International Social Magnetic Resonance in Medicine 30.

Google Scholar

Near, J., Edden, R., Evans, C. J., Paquin, R., Harris, A., and Jezzard, P. (2015). Frequency and phase drift correction of magnetic resonance spectroscopy data by spectral registration in the time domain. Magn. Reson. Med. 73, 44–50. doi: 10.1002/mrm.25094

PubMed Abstract | CrossRef Full Text | Google Scholar

Near, J., Harris, A. D., Juchem, C., Kreis, R., Marjańska, M., Öz, G., et al. (2021). Preprocessing, analysis and quantification in single-voxel magnetic resonance spectroscopy: experts' consensus recommendations. NMR Biomed. 34, e4257. doi: 10.1002/nbm.4257

PubMed Abstract | CrossRef Full Text | Google Scholar

Nilsson, M., Lasič, S., Drobnjak, I., Topgaard, D., and Westin, C.-F. (2017). Resolution limit of cylinder diameter estimation by diffusion MRI: The impact of gradient waveform and orientation dispersion. NMR Biomed. 30, e3711. doi: 10.1002/nbm.3711

PubMed Abstract | CrossRef Full Text | Google Scholar

Nilsson, M., Szczepankiewicz, F., Lampinen, B., Ahlgren, A., de Martins, J. P., Lasic, S., et al. (2018). “An open-source framework for analysis of multidimensional diffusion MRI data implemented in MATLAB, in Proceedings of the International Social Magnetic Resonance in Medicine 5355.

Google Scholar

Oeltzschner, G., Zöllner, H. J., Hui, S. C., Mikkelsen, M., Saleh, M. G., Tapper, S., et al. (2020). Osprey: Open-source processing, reconstruction &estimation of magnetic resonance spectroscopy data. J. Neurosci. Methods 343, 108827. doi: 10.1016/j.jneumeth.2020.108827

CrossRef Full Text | Google Scholar

Ogura, A., Maeda, F., Miyai, A., Hayashi, K., and Hongoh, T. (2006). Effect of vibration caused by time-varying magnetic fields on diffusion-weighted mri. Nihon Hoshasen Gijutsu Gakkai Zasshi 62, 565–569. doi: 10.6009/jjrt.62.565

PubMed Abstract | CrossRef Full Text | Google Scholar

Palacios, E. M., Martin, A. J., Boss, M. A., Ezekiel, F., Chang, Y. S., Yuh, E. L., et al. (2017). Toward precision and reproducibility of diffusion tensor imaging: a multicenter diffusion phantom and traveling volunteer study. Am. J. Neuroradiol. 38, 537–545. doi: 10.3174/ajnr.A5025

PubMed Abstract | CrossRef Full Text | Google Scholar

Palombo, M., Ligneul, C., and Valette, J. (2017). Modeling diffusion of intracellular metabolites in the mouse brain up to very high diffusion-weighting: diffusion in long fibers (almost) accounts for non-monoexponential attenuation. Magn. Reson. Med. 77, 343–350. doi: 10.1002/mrm.26548

PubMed Abstract | CrossRef Full Text | Google Scholar

Palombo, M., Shemesh, N., Ronen, I., and Valette, J. (2018). Insights into brain microstructure from in vivo dw-mrs. Neuroimage 182, 97–116. doi: 10.1016/j.neuroimage.2017.11.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Patterson, D. M., Padhani, A. R., and Collins, D. J. (2008). Technology insight: water diffusion mria potential new biomarker of response to cancer therapy. Nat. Clin. Pract. Oncol. 5, 220–233. doi: 10.1038/ncponc1073

PubMed Abstract | CrossRef Full Text | Google Scholar

Penny, W. D., Friston, K. J., Ashburner, J. T., Kiebel, S. J., and Nichols, T. E. (2011). Statistical Parametric Mapping: The Analysis of Functional Brain Images. London: Elsevier.

Google Scholar

Ronen, I., Budde, M., Ercan, E., Annese, J., Techawiboonwong, A., and Webb, A. (2014). Microstructural organization of axons in the human corpus callosum quantified by diffusion-weighted magnetic resonance spectroscopy of n-acetylaspartate and post-mortem histology. Brain Struct. Funct. 219, 1773–1785. doi: 10.1007/s00429-013-0600-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Ronen, I., and Valette, J. (2007). Diffusion-weighted magnetic resonance spectroscopy. EMagRes. 4, 733–750. doi: 10.1002/9780470034590.emrstm1471

CrossRef Full Text | Google Scholar

Saleh, M. G., Edden, R. A., Chang, L., and Ernst, T. (2020). Motion correction in magnetic resonance spectroscopy. Magn. Reson. Med. 84, 2312–2326. doi: 10.1002/mrm.28287

PubMed Abstract | CrossRef Full Text | Google Scholar

Setsompop, K., Kimmlingen, R., Eberlein, E., Witzel, T., Cohen-Adad, J., McNab, J. A., et al. (2013). Pushing the limits of in vivo diffusion MRI for the human connectome project. Neuroimage 80, 220–233. doi: 10.1016/j.neuroimage.2013.05.078

PubMed Abstract | CrossRef Full Text | Google Scholar

Simpson, R., Devenyi, G. A., Jezzard, P., Hennessy, T. J., and Near, J. (2017). Advanced processing and simulation of mrs data using the fid appliance (fid-a) an open source, matlab-based toolkit. Magn. Reson. Med. 77, 23–33. doi: 10.1002/mrm.26091

PubMed Abstract | CrossRef Full Text | Google Scholar

Şimşek, K., Döring, A., Pampel, A., Möller, H. E., and Kreis, R. (2020). “Diffusion-weighted mrs at short te using a connectom system: non-gaussian metabolite diffusion and macromolecular signals in human brain,” in Proceedings of the Virtual Conference of ISMRM 364.

Google Scholar

Şimşek, K., Döring, A., Pampel, A., Möller, H. E., and Kreis, R. (2021). “Simultaneous spectral and bi-exponential diffusion modeling of doubly motion-corrected human brain spectra with very high b-values,” in Proceedings of the Annual Meeting of ISMRM 4263.

Google Scholar

Şimşek, K., Döring, A., Pampel, A., Möller, H. E., and Kreis, R. (2022). Macromolecular background signal and non-gaussian metabolite diffusion determined in human brain using ultra-high diffusion weighting. Magn. Reson. Med. 88, 1962–1977. doi: 10.1002/mrm.29367

PubMed Abstract | CrossRef Full Text | Google Scholar

Szczepankiewicz, F., Westin, C.-F., and Nilsson, M. (2019). Maxwell-compensated design of asymmetric gradient waveforms for tensor-valued diffusion encoding. Magn. Reson. Med. 82, 1424–1437. doi: 10.1002/mrm.27828

PubMed Abstract | CrossRef Full Text | Google Scholar

Tal, A. (2023). The future is 2d: spectral-temporal fitting of dynamic mrs data provides exponential gains in precision over conventional approaches. Magn. Reson. Med. 89, 499–507. doi: 10.1002/mrm.29456

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanhamme, L., Fierro, R. D., Van Huffel, S., and de Beer, R. (1998). Fast removal of residual water in proton spectra. J. Magn. Reson. 132, 197–203. doi: 10.1006/jmre.1998.1425

PubMed Abstract | CrossRef Full Text | Google Scholar

Weidlich, D., Honecker, J., Gmach, O., Wu, M., Burgkart, R., Ruschke, S., et al. (2019). Measuring large lipid droplet sizes by probing restricted lipid diffusion effects with diffusion-weighted mrs at 3T. Magn. Reson. Med. 81, 3427–3439. doi: 10.1002/mrm.27651

PubMed Abstract | CrossRef Full Text | Google Scholar

Weidlich, D., Zamskiy, M., Maeder, M., Ruschke, S., Marburg, S., and Karampinos, D. C. (2020). Reduction of vibration-induced signal loss by matching mechanical vibrational states: Application in high b-value diffusion-weighted mrs. Magn. Reson. Med. 84, 39–51. doi: 10.1002/mrm.28128

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilm, B. J., Nagy, Z., Barmet, C., Vannesjo, S. J., Kasper, L., Haeberlin, M., et al. (2015). Diffusion MRI with concurrent magnetic field monitoring. Magn. Reson. Med. 74, 925–933. doi: 10.1002/mrm.25827

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilson, M., Andronesi, O., Barker, P. B., Bartha, R., Bizzi, A., Bolan, P. J., et al. (2019). Methodological consensus on clinical proton mrs of the brain: review and recommendations. Magn. Reson. Med. 82, 527–550. doi: 10.1002/mrm.27742

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilson, M., Reynolds, G., Kauppinen, R. A., Arvanitis, T. N., and Peet, A. C. (2011). A constrained least-squares approach to the automated quantitation of in vivo 1h magnetic resonance spectroscopy data. Magn. Reson. Med. 65, 1–12. doi: 10.1002/mrm.22579

PubMed Abstract | CrossRef Full Text | Google Scholar

Wyss, P. O., Bianchini, C., Scheidegger, M., Giapitzakis, I. A., Hock, A., Fuchs, A., et al. (2018). In vivo estimation of transverse relaxation time constant (t2) of 17 human brain metabolites at 3T. Magn. Reson. Med. 80, 452–461. doi: 10.1002/mrm.27067

PubMed Abstract | CrossRef Full Text | Google Scholar

Yablonskiy, D. A., and Sukstanskii, A. L. (2010). Theoretical models of the diffusion weighted mr signal. NMR Biomed. 23, 661–681. doi: 10.1002/nbm.1520

PubMed Abstract | CrossRef Full Text | Google Scholar

Zöllner, H. J., Davies-Jenkins, C. W., Murali-Manohar, S., Gong, T., Hui, S. C., Song, Y., et al. (2023). Feasibility and implications of using subject-specific macromolecular spectra to model short echo time magnetic resonance spectroscopy data. NMR Biomed. 36, e4854. doi: 10.1002/nbm.4854

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: diffusion-weighted MRS, ultra-strong gradients, gradient non-uniformity, eddy currents, metabolites

Citation: Davies-Jenkins CW, Döring A, Fasano F, Kleban E, Mueller L, Evans CJ, Afzali M, Jones DK, Ronen I, Branzoli F and Tax CMW (2023) Practical considerations of diffusion-weighted MRS with ultra-strong diffusion gradients. Front. Neurosci. 17:1258408. doi: 10.3389/fnins.2023.1258408

Received: 14 July 2023; Accepted: 03 November 2023;
Published: 07 December 2023.

Edited by:

Yi Zhang, Zhejiang University, China

Reviewed by:

Min Wang, Zhejiang University, China
Yanqin Lin, Xiamen University, China

Copyright © 2023 Davies-Jenkins, Döring, Fasano, Kleban, Mueller, Evans, Afzali, Jones, Ronen, Branzoli and Tax. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Christopher W. Davies-Jenkins, cdavies9@jh.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.