MR susceptibility imaging

Quantitative Susceptibility Mapping Using Structural Feature Based Collaborative Reconstruction (SFCR) in the Human Brain

The reconstruction of MR quantitative susceptibility mapping (QSM) from local phase measurements is an ill posed inverse problem and different regularization strategies incorporating a priori information extracted from magnitude and phase images have been proposed. However, the anatomy observed in magnitude and phase images does not always coincide spatially with that in susceptibility maps, which could give erroneous estimation in the reconstructed susceptibility map. In this paper, we develop a structural feature based collaborative reconstruction (SFCR) method for QSM including both magnitude and susceptibility based information. The SFCR algorithm is composed of two consecutive steps corresponding to complementary reconstruction models, each with a structural feature based l 1 norm constraint and a voxel fidelity based l 2 norm constraint, which allows both the structure edges and tiny features to be recovered, whereas the noise and artifacts could be reduced. In the M-step, the initial susceptibility map is reconstructed by employing a k -space based compressed sensing model incorporating magnitude prior. In the S-step, the susceptibility map is fitted in spatial domain using weighted constraints derived from the initial susceptibility map from the M-step. Simulations and in vivo human experiments at 7T MRI show that the SFCR method provides high quality susceptibility maps with improved RMSE and MSSIM. Finally, the susceptibility values of deep gray matter are analyzed in multiple head positions, with the supine position most approximate to the gold standard COSMOS result.

Imaging whole-brain cytoarchitecture of mouse with MRI-based quantitative susceptibility mapping

The proper microstructural arrangement of complex neural structures is essential for establishing the functional circuitry of the brain. We present an MRI method to resolve tissue microstructure and infer brain cytoarchitecture by mapping the magnetic susceptibility in the brain at high resolution. This is possible because of the heterogeneous magnetic susceptibility created by varying concentrations of lipids, proteins and irons from the cell membrane to cytoplasm. We demonstrate magnetic susceptibility maps at a nominal resolution of 10-μm isotropic, approaching the average cell size of a mouse brain. The maps reveal many detailed structures including the retina cell layers, olfactory sensory neurons, barrel cortex, cortical layers, axonal fibers in white and gray matter. Olfactory glomerulus density is calculated and structural connectivity is traced in the optic nerve, striatal neurons, and brainstem nerves. The method is robust and can be readily applied on MRI scanners at or above 7T.

Advances in MRI-based computational neuroanatomy: from morphometry to in-vivo histology

PURPOSE OF REVIEW

Current computational neuroanatomy based on MRI focuses on morphological measures of the brain. We present recent methodological developments in quantitative MRI (qMRI) that provide standardized measures of the brain, which go beyond morphology. We show how biophysical modelling of qMRI data can provide quantitative histological measures of brain tissue, leading to the emerging field of in-vivo histology using MRI (hMRI).

RECENT FINDINGS

qMRI has greatly improved the sensitivity and specificity of computational neuroanatomy studies. qMRI metrics can also be used as direct indicators of the mechanisms driving observed morphological findings. For hMRI, biophysical models of the MRI signal are being developed to directly access histological information such as cortical myelination, axonal diameters or axonal g-ratio in white matter. Emerging results indicate promising prospects for the combined study of brain microstructure and function.

SUMMARY

Non-invasive brain tissue characterization using qMRI or hMRI has significant implications for both research and clinics. Both approaches improve comparability across sites and time points, facilitating multicentre/longitudinal studies and standardized diagnostics. hMRI is expected to shed new light on the relationship between brain microstructure, function and behaviour, both in health and disease, and become an indispensable addition to computational neuroanatomy.

Quantitative [Fe]MRI of PSMA-targeted SPIONs specifically discriminates among prostate tumor cell types based on their PSMA expression levels

We report the development, experimental verification, and application of a general theory called [Fe]MRI (pronounced fem-ree) for the non-invasive, quantitative molecular magnetic resonance imaging (MRI) of added magnetic nanoparticles or other magnetic contrast agents in biological tissues and other sites. [Fe]MRI can easily be implemented on any MRI instrument, requiring only measurements of the background nuclear magnetic relaxation times (T₁, T₂) of the tissue of interest, injection of the magnetic particles, and the subsequent acquisition of a pair of T₁-weighted and T₂-weighted images. These images, converted into contrast images, are subtracted to yield a contrast difference image proportional to the absolute nanoparticle, iron concentration, ([Fe]) image. [Fe]MRI was validated with the samples of superparamagnetic iron oxide nanoparticles (SPIONs) both in agarose gels and bound to human prostate tumor cells. The [Fe]MRI measurement of the binding of anti-prostate specific membrane antigen (PSMA) conjugated SPIONs to PSMA-positive LNCaP and PSMA-negative DU145 cells in vitro allowed a facile discrimination among prostate tumor cell types based on their PSMA expression level. The low [Fe] detection limit of ~2 μM for SPIONs allows sensitive MRI of added iron at concentrations considerably below the US Food and Drug Administration’s human iron dosage guidelines (<90 μM, 5 mg/kg).

Foundations of MRI phase imaging and processing for Quantitative Susceptibility Mapping (QSM)

Quantitative Susceptibility Mapping (QSM) is a novel MRI based technique that relies on estimates of the magnetic field distribution in the tissue under examination. Several sophisticated data processing steps are required to extract the magnetic field distribution from raw MRI phase measurements. The objective of this review article is to provide a general overview and to discuss several underlying assumptions and limitations of the pre-processing steps that need to be applied to MRI phase data before the final field-to-source inversion can be performed. Beginning with the fundamental relation between MRI signal and tissue magnetic susceptibility this review covers the reconstruction of magnetic field maps from multi-channel phase images, background field correction, and provides an overview of state of the art QSM solution strategies.

Isolated cortical vein thrombosis: case series

OBJECT

Isolated cortical vein thrombosis (ICVT) accounts for less than 1% of all cerebral infarctions. ICVT may cause irreversible parenchymal damage, rendering early and accurate diagnosis critical. This case series and literature review presents the clinical and radiological findings in 7 patients with ICVT, and highlights risk factors and imaging modalities that may be most beneficial in rendering an accurate and timely diagnosis.

METHODS

Patients with CT and MRI findings consistent with ICVT examined between January 2011 and June 2014 were included in this retrospective review.

RESULTS

Seven patients (5 females, 2 males), ranging in age from 11 months to 34 years, met the inclusion criteria. The most common clinical presentations were headaches (n = 4) and seizures (n = 3). The most common comorbidities noted in these patients were hypercoagulable states (n = 4) and intracranial hypotension (n = 3). Five patients had intraparenchymal involvement. CT suggested the correct diagnosis in 4 patients, and MRI confirmed the diagnosis in all 7 patients. All patients who received anticoagulation therapy (n = 5) experienced complete resolution of their symptoms.

CONCLUSIONS

The majority of these patients were adult females, consistent with published data. Seizures and headaches were the most common presenting symptoms. Hypercoagulable state and intracranial hypotension, both known risk factors for thrombosis, were the most commonly noted ICVT risk factors. Intraparenchymal involvement was prevalent in nearly all ICVT cases and presented as vasogenic edema, early intraparenchymal hemorrhage, or hemorrhagic venous infarction. Susceptibility-weighted imaging was the most sensitive imaging technique in diagnosing ICVT.

Quantitative susceptibility mapping (QSM): Decoding MRI data for a tissue magnetic biomarker

In MRI, the main magnetic field polarizes the electron cloud of a molecule, generating a chemical shift for observer protons within the molecule and a magnetic susceptibility inhomogeneity field for observer protons outside the molecule. The number of water protons surrounding a molecule for detecting its magnetic susceptibility is vastly greater than the number of protons within the molecule for detecting its chemical shift. However, the study of tissue magnetic susceptibility has been hindered by poor molecular specificities of hitherto used methods based on MRI signal phase and T2* contrast, which depend convolutedly on surrounding susceptibility sources. Deconvolution of the MRI signal phase can determine tissue susceptibility but is challenged by the lack of MRI signal in the background and by the zeroes in the dipole kernel. Recently, physically meaningful regularizations, including the Bayesian approach, have been developed to enable accurate quantitative susceptibility mapping (QSM) for studying iron distribution, metabolic oxygen consumption, blood degradation, calcification, demyelination, and other pathophysiological susceptibility changes, as well as contrast agent biodistribution in MRI. This paper attempts to summarize the basic physical concepts and essential algorithmic steps in QSM, to describe clinical and technical issues under active development, and to provide references, codes, and testing data for readers interested in QSM.

A method for estimating and removing streaking artifacts in quantitative susceptibility mapping

Quantitative susceptibility mapping (QSM) is a novel MRI method for quantifying tissue magnetic property. In the brain, it reflects the molecular composition and microstructure of the local tissue. However, susceptibility maps reconstructed from single-orientation data still suffer from streaking artifacts which obscure structural details and small lesions. We propose and have developed a general method for estimating streaking artifacts and subtracting them from susceptibility maps. Specifically, this method uses a sparse linear equation and least-squares (LSQR)-algorithm-based method to derive an initial estimation of magnetic susceptibility, a fast quantitative susceptibility mapping method to estimate the susceptibility boundaries, and an iterative approach to estimate the susceptibility artifact from ill-conditioned k-space regions only. With a fixed set of parameters for the initial susceptibility estimation and subsequent streaking artifact estimation and removal, the method provides an unbiased estimate of tissue susceptibility with negligible streaking artifacts, as compared to multi-orientation QSM reconstruction. This method allows for improved delineation of white matter lesions in patients with multiple sclerosis and small structures of the human brain with excellent anatomical details. The proposed methodology can be extended to other existing QSM algorithms.

The role of iron in brain ageing and neurodegenerative disorders

SUMMARY

In the CNS, iron in several proteins is involved in many important processes such as oxygen transportation, oxidative phosphorylation, myelin production, and the synthesis and metabolism of neurotransmitters. Abnormal iron homoeostasis can induce cellular damage through hydroxyl radical production, which can cause the oxidation and modification of lipids, proteins, carbohydrates, and DNA. During ageing, different iron complexes accumulate in brain regions associated with motor and cognitive impairment. In various neurodegenerative diseases, such as Alzheimer's disease and Parkinson's disease, changes in iron homoeostasis result in altered cellular iron distribution and accumulation. MRI can often identify these changes, thus providing a potential diagnostic biomarker of neurodegenerative diseases. An important avenue to reduce iron accumulation is the use of iron chelators that are able to cross the blood-brain barrier, penetrate cells, and reduce excessive iron accumulation, thereby affording neuroprotection.

Phase-corrected bipolar gradients in multi-echo gradient-echo sequences for quantitative susceptibility mapping

OBJECTIVE

Large echo spacing of unipolar readout gradients in current multi-echo gradient-echo (GRE) sequences for mapping fields in quantitative susceptibility mapping (QSM) can be reduced using bipolar readout gradients thereby improving acquisition efficiency.

MATERIALS AND METHODS

Phase discrepancies between odd and even echoes in the bipolar readout gradients caused by non-ideal gradient behaviors were measured, modeled as polynomials in space and corrected for accordingly in field mapping. The bipolar approach for multi-echo GRE field mapping was compared with the unipolar approach for QSM.

RESULTS

The odd-even-echo phase discrepancies were approximately constant along the phase encoding direction and linear along the readout and slice-selection directions. A simple linear phase correction in all three spatial directions was shown to enable accurate QSM of the human brain using a bipolar multi-echo GRE sequence. Bipolar multi-echo acquisition provides QSM in good quantitative agreement with unipolar acquisition while also reducing noise.

CONCLUSION

With a linear phase correction between odd-even echoes, bipolar readout gradients can be used in multi-echo GRE sequences for QSM.

The role of iron in brain ageing and neurodegenerative disorders

In the CNS, iron in several proteins is involved in many important processes such as oxygen transportation, oxidative phosphorylation, myelin production, and the synthesis and metabolism of neurotransmitters. Abnormal iron homoeostasis can induce cellular damage through hydroxyl radical production, which can cause the oxidation and modification of lipids, proteins, carbohydrates, and DNA. During ageing, different iron complexes accumulate in brain regions associated with motor and cognitive impairment. In various neurodegenerative diseases, such as Alzheimer's disease and Parkinson's disease, changes in iron homoeostasis result in altered cellular iron distribution and accumulation. MRI can often identify these changes, thus providing a potential diagnostic biomarker of neurodegenerative diseases. An important avenue to reduce iron accumulation is the use of iron chelators that are able to cross the blood-brain barrier, penetrate cells, and reduce excessive iron accumulation, thereby affording neuroprotection.

Localizing the human primary auditory cortex in vivo using structural MRI

Currently there are no routine methods to delineate the primary auditory cortex (PAC) of humans in vivo. Due to the large differences in the location of the PAC between subjects, labels derived from post-mortem brains may be inaccurate when applied to different samples of in vivo brains. Recent magnetic resonance (MR) imaging studies suggested that MR-tissue properties can be used to define the location of the PAC region in vivo. The basis for such an approach is that the PAC region is more strongly myelinated than the secondary areas. We developed a fully automatic method to identify the PAC in conventional anatomical data using a combination of two complementary MR contrasts, i.e., T1 and T2, at 3T with 0.7mm isotropic resolution. Our algorithm maps the anatomical MR data to reconstructed cortical surfaces and uses a classification approach to create an artificial contrast that is highly sensitive to the effects of an increased myelination of the cortex. Consistent with the location of the PAC defined in post-mortem brains, we found a compact region on the medial two thirds of Heschl's gyrus in both hemispheres of all 39 subjects. With further improvements in signal-to-noise ratio of the anatomical data and manual correction of segmentation errors, the results suggest that the primary auditory cortex can be defined in the living brain of single subjects.

Ultrafast in vivo diffusion imaging of stroke at 21.1 T by spatiotemporal encoding

PURPOSE

This study quantifies in vivo ischemic stroke brain injuries in rats using ultrahigh-field single-scan MRI methods to assess variations in apparent diffusion coefficients (ADCs).

METHODS

Magnitude and diffusion-weighted spatiotemporally encoded imaging sequences were implemented on a 21.1 T imaging system, and compared with spin-echo and echo-planar imaging diffusion-weighted imaging strategies. ADC maps were calculated and used to evaluate the sequences according to the statistical comparisons of the ipsilateral and contralateral ADC measurements at 24, 48, and 72 h poststroke.

RESULTS

Susceptibility artifacts resulting from normative anatomy and pathological stroke conditions were particularly intense at 21.1 T. These artifacts strongly distorted single-shot diffusion-weighted echo-planar imaging experiments, but were reduced in four-segment interleaved echo-planar imaging acquisitions. By contrast, nonsegmented diffusion-weighted spatiotemporally encoded images were largely immune to field-dependent artifacts. Effects of stroke were apparent in both magnitude images and ADC maps of all sequences. When stroke recovery was followed by ADC variations, spatiotemporally encoded, echo-planar imaging, and spin-echo acquisitions revealed statistically significant increase in ADCs.

CONCLUSIONS

Consideration of experiment duration, image quality, and mapped ADC values provided by spatiotemporally encoded demonstrates that this single-shot acquisition is a method of choice for high-throughput, ultrahigh-field in vivo stroke quantification.

Susceptibility tensor imaging of the kidney and its microstructural underpinnings

PURPOSE

The purpose of this study was to determine whether susceptibility tensor imaging (STI) could overcome limitations of current techniques to detect tubules throughout the kidney.

METHODS

Normal mouse kidneys (n = 4) were imaged at 9.4T using a three-dimensional gradient multi-echo sequence (55-micron isotropic resolution). Phase images from 12 orientations were obtained to compute the susceptibility tensor. Diffusion tensor imaging (DTI) with 12 encoding directions was compared with STI. Tractography was performed to visualize and track the course of tubules with DTI and STI. Confocal microscopy was used to identify which tubular segments of the nephron were detected by DTI and STI.

RESULTS

Diffusion anisotropy was limited to the inner medulla of the kidney. DTI did not find a significant number of coherent tubular tracks in the outer medulla or cortex. With STI, we found strong susceptibility anisotropy and many tracks in the inner and outer medulla and in limited areas of the cortex.

CONCLUSION

STI was able to track tubules throughout the kidney, whereas DTI was limited to the inner medulla. STI provides a novel contrast mechanism related to local tubule microstructure and may offer a powerful method to study the nephron.

Magnetic resonance imaging at ultrahigh fields

Since the introduction of 4 T human systems in three academic laboratories circa 1990, rapid progress in imaging and spectroscopy studies in humans at 4 T and animal model systems at 9.4 T have led to the introduction of 7 T and higher magnetic fields for human investigation at about the turn of the century. Work conducted on these platforms has demonstrated the existence of significant advantages in SNR and biological information content at these ultrahigh fields, as well as the presence of numerous challenges. Primary difference from lower fields is the deviation from the near field regime; at the frequencies corresponding to hydrogen resonance conditions at ultrahigh fields, the RF is characterized by attenuated traveling waves in the human body, which leads to image nonuniformities for a given sample-coil configuration because of interferences. These nonuniformities were considered detrimental to the progress of imaging at high field strengths. However, they are advantageous for parallel imaging for signal reception and parallel transmission, two critical technologies that account, to a large extend, for the success of ultrahigh fields. With these technologies, and improvements in instrumentation and imaging methods, ultrahigh fields have provided unprecedented gains in imaging of brain function and anatomy, and started to make inroads into investigation of the human torso and extremities. As extensive as they are, these gains still constitute a prelude to what is to come given the increasingly larger effort committed to ultrahigh field research and development of ever better instrumentation and techniques.

Integrated Laplacian-based phase unwrapping and background phase removal for quantitative susceptibility mapping

Quantitative susceptibility mapping (QSM) is a recently developed MRI technique that provides a quantitative measure of tissue magnetic susceptibility. To compute tissue magnetic susceptibilities based on gradient echoes, QSM requires reliable unwrapping of the measured phase images and removal of contributions caused by background susceptibilities. Typically, the two steps are performed separately. Here, we present a method that simultaneously performs phase unwrapping and HARmonic (background) PhasE REmovaL using the LAplacian operator (HARPERELLA). Both numerical simulations and in vivo human brain images show that HARPERELLA effectively removes both phase wraps and background phase, whilst preserving all low spatial frequency components originating from brain tissues. When compared with other QSM phase preprocessing techniques, such as path-based phase unwrapping followed by background phase removal, HARPERELLA preserves the tissue phase signal in gray matter, white matter and cerebrospinal fluid with excellent robustness, providing a convenient and accurate solution for QSM. The proposed algorithm is provided, together with QSM and susceptibility tensor imaging (STI) tools, in a shared software package named 'STI Suite'.

Integrated Laplacian-based phase unwrapping and background phase removal for quantitative susceptibility mapping

Quantitative susceptibility mapping (QSM) is a recently developed MRI technique that provides a quantitative measure of tissue magnetic susceptibility. To compute tissue magnetic susceptibilities based on gradient echoes, QSM requires reliable unwrapping of the measured phase images and removal of contributions caused by background susceptibilities. Typically, the two steps are performed separately. Here, we present a method that simultaneously performs phase unwrapping and HARmonic (background) PhasE REmovaL using the LAplacian operator (HARPERELLA). Both numerical simulations and in vivo human brain images show that HARPERELLA effectively removes both phase wraps and background phase, whilst preserving all low spatial frequency components originating from brain tissues. When compared with other QSM phase preprocessing techniques, such as path-based phase unwrapping followed by background phase removal, HARPERELLA preserves the tissue phase signal in gray matter, white matter and cerebrospinal fluid with excellent robustness, providing a convenient and accurate solution for QSM. The proposed algorithm is provided, together with QSM and susceptibility tensor imaging (STI) tools, in a shared software package named 'STI Suite'.

Differential developmental trajectories of magnetic susceptibility in human brain gray and white matter over the lifespan

As indicated by several recent studies, magnetic susceptibility of the brain is influenced mainly by myelin in the white matter and by iron deposits in the deep nuclei. Myelination and iron deposition in the brain evolve both spatially and temporally. This evolution reflects an important characteristic of normal brain development and ageing. In this study, we assessed the changes of regional susceptibility in the human brain in vivo by examining the developmental and ageing process from 1 to 83 years of age. The evolution of magnetic susceptibility over this lifespan was found to display differential trajectories between the gray and the white matter. In both cortical and subcortical white matter, an initial decrease followed by a subsequent increase in magnetic susceptibility was observed, which could be fitted by a Poisson curve. In the gray matter, including the cortical gray matter and the iron-rich deep nuclei, magnetic susceptibility displayed a monotonic increase that can be described by an exponential growth. The rate of change varied according to functional and anatomical regions of the brain. For the brain nuclei, the age-related changes of susceptibility were in good agreement with the findings from R2* measurement. Our results suggest that magnetic susceptibility may provide valuable information regarding the spatial and temporal patterns of brain myelination and iron deposition during brain maturation and ageing.

Imaging of iron

Magnetic resonance imaging (MRI) enables a noninvasive in vivo quantification of iron in various organs. Several techniques have been developed that detect signal alterations derived mainly from the magnetic properties of ferritin and hemosiderin, the major iron storage compounds. High magnetic susceptibility of ferritin shortens the transversal relaxation time of nearby water protons and thus induces a focal signal extinction of iron-rich areas in T2-weighted (T2w) MRI. T2w tissue contrast is additionally influenced by other factors such as water content, myelin density, and the presence of other metals. Therefore, more specific methods are needed with higher specificity to iron. These in vivo techniques can be divided into three groups: relaxometry, magnetic field correlation imaging and phase-based contrast covering susceptibility-weighted imaging, and quantitative susceptibility mapping. The differential diagnosis of various neurological disorders is aided by characteristic patterns of iron depositions. Reliable estimates of cerebral tissue iron concentration are equally important in studying physiological age-related as well as pathological conditions in neurodegenerative, neuroinflammatory, and vascular diseases. In the future, monitoring changes in iron storage and content may serve as sensitive biomarker for diagnosis as well as treatment monitoring.