We implemented and optimized an Off-Resonance Saturation (ORS) pulse sequence on a 7T- preclinical scanner, with the aim of quantifying iron present in brain phantoms doped with ferritin-bound and magnetite-bound iron. These mineralized iron forms are known to be involved in ageing, and also in toxic cellular pathways in neurodegenerative diseases. We show that the ORS method offers the possibility to quantify, and possibly differentiate, iron ions in brains affected by neurodegenerative diseases which are characterized by disturbed iron homeostasis. Additionally, our results are discussed in the light of susceptibility values of the same minerals.
While iron dis-homeostasis is known to occur in the brain of patients with neurodegenerative diseases, including Alzheimer’s disease (AD), only recently has attention been paid to the presence of different forms of mineralized iron1,2,3. Experimental evidence suggests that magnetite is present in the core of amyloid plaques, where it may contribute to initiate/progress the amyloid cascade and induce toxicity4. Another crucial form of iron, in terms of MRI contrast and biological relevance, is ferritin: a nanoglobular protein responsible for iron-storage. Iron and ferritin contents are altered across the cortices, in AD5,6.
Gradient-Echo images are routinely employed to detect iron overload in tissue7 and to obtain Quantitative Susceptibility Mapping (QSM) maps, which have been used to indirectly assess iron ex-vivo and in-vivo8,9. However, these methods do not offer speciation of tissue iron.
Here we develop and optimize an Off-Resonance Saturation (ORS) pulse sequence to obtain absolute concentrations of iron contained in ferritin and magnetite nanoparticles. We discuss our results and compare them with QSM values.
1.5%-agarose brain phantoms, one containing different concentrations of 20nm-magnetite nanoparticles, a second containing different concentrations of horse spleen ferritin, and a third one containing both species, were fabricated. Iron concentrations were chosen to cover the range of physiological and pathological conditions. The ORS preparation module consisted of multiple saturation sech-pulses, with 90° phase cycling in-between pulses, followed by a spoiler gradient in the slice-selection direction. The acquisition consisted of a post-excitation SSFP sequence, with four segments. An additional 90° phase cycling in between segments was implemented to improve spoiling. The Z-spectrum was inspected to assess the presence of magnetization transfer (MT) effects, which may confound the analysis. Afterward, the contrast curve was fitted to the following expression:
$$C=S_0 \frac{N}{\pi^2}\left[ \arctan \left( \frac{\alpha}{B_{eq}[Fe]}(\omega_0 + \Delta \omega_0/2)\right)-\arctan \left( \frac{\alpha}{B_{eq}[Fe]}(\omega_0 -\Delta \omega_0/2) \right) \right]+off$$
Where S0 is the magnitude of the unsaturated image, N is the number of ORS pulses, α a numerical factor10, ω0 and Δω the frequency and bandwidth of the ORS pulse respectively, [Fe] the iron concentration and off a constant. These last two are the fitting parameters. The equatorial field, Beq, was fixed to 167 mT for magnetite10 and 1.2 mT for ferritin11. Δω was fixed to 500 Hz, and ω0 was varied. Multi-echo gradient-Echo images were acquired on the magnetite-doped phantom, for comparative analysis. STI Suite toolbox was used for QSM reconstruction12.
Fig.1 shows positive contrast on the first two samples. For these, the fit to the equation above returned the nominal [Fe], within the error bar (Fig. 2). For the mixed sample, a weighted fit was used to extract [Fe] bound to ferritin and magnetite, respectively. The fit returned the concentration of magnetite-iron with a good accuracy, while it underestimated the ferritin-iron fraction (Fig. 3).
QSM maps were used to extract mean susceptibility values (χ) from the magnetite sample. A correlative slope of 0.65 ± 0.10 ppb/µg iron/ml was found (Fig. 4), in agreement with the calculated value of 0.5 ppb/µg iron/ml13.
Figure 1. Unsaturated and ORS images (saturated at ω0=300Hz and subtracted from the unsaturated one) for magnetite and ferritin-doped phantoms. The iron concentration, in mM, is written in the sample compartments. Top row: Acquisition parameters for the magnetite-phantom were TE/TR=1.634/3.268 ms and FA=10°, voxel size=0.25x0.25x1.5 mm3. ORS images were obtained by varying ω0 in a [-700, 700] Hz range. Three ORS pulses were used. Bottom row: Acquisition parameters for the ferritin-phantom were as above except for TE/TR=20/40 ms, voxel size=0.25x0.25x1mm3. Five ORS pulses were used.