Reciprocity based metabolite quantification at 3T
Niklaus Zoelch1, Andreas Hock1,2, and Anke Henning1,3

1Institute for Biomedical Engineering, UZH and ETH Zurich, Zurich, Switzerland, 2Department of Psychiatry, Psychotherapy and Psychosomatics, University of Zurich, Zurich, Switzerland, 3Max Planck Institute for Biological Cybernetics, Tuebingen, Germany

Synopsis

At 1.5 T reciprocity principle based quantification strategies have been successfully used to quantify brain metabolites. But these methods all rely on the assumption that the magnitude of the RF transmission field B1+ and the reception field B1- are equal at all points in the subject. This is not true at higher field strengths and for example differences in the concentrations measured in the left and right hemisphere are observed when these methods are directly applied at higher fields. Here a further development is presented, proposing a correction for deviations of B1+ from B1- to allow concentration measurements at 3T and even higher field strength without the need of assumptions about concentrations of an internal reference. The obtained metabolite concentrations in vivo in 31 healthy volunteers highly agree with values estimated with internal water referencing, demonstrating the capabilities of this new method, which might make concentration measurements in diseased tissue more reliable.

Introduction

When the metabolite signal $$$S_{m}$$$ detected with MRS is compared to a reference signal measured in a phantom with known concentrations and relaxation rates, it is in principle possible to determine in vivo metabolite concentrations without the need of assumptions about an internal reference. However, the changes in sensitivity, i.e. the relation between the transversal magnetization and the detected signal, between the in vivo and the phantom measurement need to be carefully considered. At low field strength one can estimate the magnitude of the $$$B_1^{+}$$$-field per unit current ($$$\mid{\hat{B_1^{+}}}\mid$$$) obtained in the VOI and use it as measure of the local sensitivity1. This approach based on the reciprocity principle2 uses the assumption that the magnitude of the reception field ($$$\mid{\hat{B_1^{-}}}\mid$$$) at a given position is equal to the magnitude of the transmission field ($$$\mid{\hat{B_1^{+}}}\mid$$$). This is true at lower fields but has been demonstrated to be wrong at higher field strengths3,4. Hence at 3T $$$\mid{\hat{B_1^{+}}}\mid$$$ is not sufficient to estimate $$$\mid{\hat{B_1^{-}}}\mid$$$ at all positions in the sample and can not be used solely to correct for sensitivity changes. In this work, we present a reciprocity based reference method for quantitative MRS at 3T. We propose to use a volume selective power optimization method5 (volPO) to determine automatically a measure of $$$\mid{\hat{B_1^{+}}}\mid$$$. This shall not only simplify the measurement of the sensitivity changes but also guarantee signal detection with an optimal flip angle. To account for the deviations of $$$\mid{\hat{B_1^{+}}}\mid$$$ from $$$\mid{\hat{B_1^{-}}}\mid$$$ we introduce a correction factor Cp using reception sensitivity maps determined from contrast minimized images6. The proposed method is extensively validated.

Methods

Experiments were performed on a Philips Achieva 3T MRI scanner with a commercial T/R head coil. With the mentioned volPO the scaling factor $$$\frac{1}{DS}$$$, which is proportional to the RF pulse voltage needed to yield a given flip angle, is determined prior to each spectroscopic scan and used as a measure of $$$\mid{\hat{B_1^{+}}}\mid$$$. To investigate the stability of this measure the ratio of the water signal $$$S_{H2O}$$$ and $$$\frac{1}{DS}$$$ was measured several times at the same position in a water containing phantom. The ability to correct for the loading of the coil with $$$\frac{1}{DS}$$$ was tested with a series of measurements in a small spherical water containing phantom, while a plastic bottle containing a NaCl solution was stepwise inserted into the coil.

In 31 healthy volunteers always three PRESS-localized voxels (TE:30ms, TR:1800ms, 256 averages, 5ml) aranged in one plane were measured. One VOI was placed in the longitudinal cerebral fissure just above the corpus callosum (position M) containing mostly gray matter and the other two were shifted by 20 mm respectively to the left (positon L) and to the right (position R) into cerebral white matter. Metabolite concentrations using $$$\frac{1}{DS}$$$ were calculated as follows: $$c_{m}^{DS}=\frac{S_m}{S_{ref,c}}\frac{R_{ref,c}}{R_{m}}\frac{DS}{DS_{c}}\frac{C_{p,c}}{C_{p}}\frac{T}{T_{c}}c_{ref,c}\frac{l_{ref,c}}{l_{m}}\frac{1}{1-f_{CSF}}\qquad\text{[Eq.1]}$$ The subscript $$$c$$$ denotes quantities from the calibration measurement, were the reference compound $$$ref$$$ with a known concentration $$$c_{ref}$$$ was measured. $$$R$$$ are approximated relaxation attenuation factors [Table1], $$$T$$$ the temperature, $$$l$$$ the number of spins per molecule and $$$f_{CSF}$$$ the volume fraction of the CSF. $$$C_p$$$ is the factor to correct deviations of $$$\mid{\hat{B_1^{+}}}\mid$$$ from $$$\mid{\hat{B_1^{-}}}\mid$$$. This factor is calculated in each VOI as the ratio of the achieved $$$\mid{B_1^{+}}\mid$$$ taken from a $$$B_1$$$-map7 to the reception sensitivity $$$Rx$$$ calculated using 2D reception sensitivity maps6. Both maps were scaled with the value measured at the center of the coil.$$C_{p}=\frac{\mid{B_{1,VOI} ^{+}}\mid}{Rx_{VOI}}\qquad\text{[Eq.2]}$$Metabolite concentrations obtained with Eq.1 are compared to values using tissue water as reference8 (IWR). Total internal water concentrations obtained with Eq.1 (omitting the last term) are compared with values predicted based on the individual voxel composition:$$c^{seg}=55.126\cdot(f_{CSF}\cdot\alpha_{CSF}+f_{GM}\cdot\alpha_{GM}+f_{WM}\cdot\alpha_{WM})\qquad\text{[Eq.3]}$$

Results/Discussion

Changes in the loading of the coil can be precisely corrected with the proposed quantity $$$\frac{1}{DS}$$$ [Figure1]. During a time of 2.5 month $$$S_{H2O}$$$/$$$\frac{1}{DS}$$$ was measured 15 times in vitro with a CV of 0.78 % demonstrating a good stability of the proposed method. The need and the effect of the correction factor $$$C_p$$$ in Eq.1 is demonstrated in vitro [Figure2] and in vivo [Figure3]. The correction with $$$C_p$$$ allows to reproduce the $$$c^{seg}$$$ values at all three measured positions. Metabolite concentrations obtained via Eq.1 and IWR highly agree as demonstrated in [Figure4]. However so far only deviations of $$$\mid{\hat{B_1^{+}}}\mid$$$ from $$$\mid{\hat{B_1^{-}}}\mid$$$ along the RL and AP direction are considered. The use of a 3D reception sensitivity map might further improve the accordance of the obtained values. In conclusion, the proposed calibration strategy allows to obtain reliable metabolite concentrations in vivo without the use of additional hardware nor assumptions about an internal reference at 3T.

Acknowledgements

No acknowledgement found.

References

[1] G Helms. A precise and user-independent quantification technique for regional comparison of single volume proton mr spectroscopy of the human brain. NMR in Biomedicine, 13(7):398–406, 2000.

[2] DI Hoult. The principle of reciprocity in signal strength calculations—a mathematical guide. Concepts in Magnetic resonance, 12(4):173–187, 2000.

[3] CM Collins, QX Yang, JH Wang, X Zhang, H Liu, S Michaeli, X-H Zhu, G Adriany, JT Vaughan, P Anderson, et al. Different excitation and reception distributions with a single-loop transmit-receive surface coil near a head-sized spherical phantom at 300 mhz. Magnetic Resonance in Medicine, 47(5):1026–1028, 2002

[4] H Watanabe. Investigation of the asymmetric distributions of rf transmission and reception fields at high static field. Magnetic Resonance in Medical Sciences, 11(2):129–135, 2012.

[5] J Wang, M Qiu, QX Yang, MB Smith, and RT Constable. Measurement and correction of transmitter and receiver induced nonuniformities in vivo. Magnetic Resonance in Medicine, 53(2):408–417, 2005.

[6] MJ Versluis, HE Kan, MA Buchem, and AG Webb. Improved signal to noise in proton spectroscopy of the human calf muscle at 7 T using localized b1 calibration. Magnetic Resonance in Medicine, 63(1):207–211, 2010.

[7] LI Sacolick, F Wiesinger, I Hancu, and MW Vogel. B1 mapping by bloch-siegert shift. Magnetic Resonance in Medicine, 63(5):1315–1322, 2010.

[8] R Kreis, T Ernst, and BD Ross. Absolute quantitation of water and metabolites in the human brain. ii. metabolite concentrations. Journal of magnetic resonance, Series B, 102(1):9–19, 1993.

[9] F Träber, W Block, R Lamerichs, J Gieseke, and HH Schild. 1h metabolite relaxation times at 3.0 tesla: Measurements of t1 and t2 values in normal brain and determination of regional differences in transverse relaxation. Journal of Magnetic Resonance Imaging, 19(5):537–545, 2004.

[10] W Zaaraoui, L Fleysher, R Fleysher, S Liu, BJ Soher, and O Gonen. Human brain-structure resolved t2 relaxation times of proton metabolites at 3 tesla. Magnetic Resonance in Medicine, 57(6):983–989, 2007.

[11] GJ Stanisz, EE Odrobina, J Pun, M Escaravage, SJ Graham, MJ Bronskill, and RM Henkelman. T1, t2 relaxation and magnetization transfer in tissue at 3t. Magnetic Resonance in Medicine, 54(3):507–512, 2005.

[12] L Chen, MA Bernstein, J Huston, and S Fain. Measurements of t1 relaxation times at 3.0 t: implications for clinical mra. In Proceedings of the 9th Annual Meeting of ISMRM, Glasgow, Scotland, 2001.

[13] SK Piechnik, J Evans, LH Bary, RG Wise, and P Jezzard. Functional changes in csf volume estimated using measurement of water t2 relaxation. Magnetic Resonance in Medicine, 61(3):579–586, 2009

Figures

[Figure1] Coil loading. 1/DS scales with coil load in the same manner as the water signal SH20 over the whole range of investigated loading conditions, which was a wider range than observed in vivo. Consequently the ratios SH2O/(1/DS) are independent of coil loading with a CV of only 0.52%.

[Figure2] Correction factor Cp in vitro. (A) A map of the ratio B1/Rx (Eq.2), which is averaged over each VOI to obtain the correction factor Cp. Only with the product 1/DS *Cp (rose crosses, B) the water signal SH20 is corrected for sensitivity differences at all VOI positions (C).

[Figure3] (A) Exemplary map of B1/Rx (Eq.2) measured in vivo to obtain the correction factor Cp. (B) Internal water concentrations obtained without Cp in Eq.1 (cDS-) and with Cp (cDS) measured at the locations (R,M,L) are compared with the cseg water concentrations for 31 volunteers. (C) Values summarized in boxplots.

[Figure4] tNAA (A) and tCr (B) in vivo concentrations obtained with Eq.1 (cDS) are plotted against respective values obtained with IWR (cIWR). The scatter plots show the regression line (solid) with R2. The limits of agreement (LoA) with 95% confidence interval are given in the Bland-Altman plot on the right.

[Table1] General form of the relaxation attenuation factors Rm, the relaxation times taken from literature and relative water densities. For the internal water signal the relaxation attenuation factor was determined based on the voxel composition in each subject. For the metabolites relaxation times for mixed GM and WM were used.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
3989