Spin-lattice relaxation times of bound water--its determination and implications for tissue discrimination

Strategically Acquired Gradient Echo (STAGE) Imaging, part IV: Constrained Reconstruction of White Noise (CROWN) Processing as a Means to Improve Signal-to-Noise in STAGE Imaging at 3 Tesla

Increasing the signal-to-noise ratio (SNR) has always been of critical importance for magnetic resonance imaging. Although increasing field strength provides a linear increase in SNR, it is more and more costly as field strength increases. Therefore, there is a major effort today to use signal processing methods to improve SNR since it is more efficient and economical. There are a variety of methods to improve SNR such as averaging the data at the expense of imaging time, or collecting the data with a lower resolution, all of these methods, including imaging processing methods, usually come at the expense of loss of image detail or image blurring. Therefore, we developed a new mathematical approach called CROWN (Constrained Reconstruction of White Noise) to enhance SNR without loss of structural detail and without affecting scanning time. In this study, we introduced and tested the concept behind CROWN specifically for STAGE (strategically acquired gradient echo) imaging. The concept itself is presented first, followed by simulations to demonstrate its theoretical effectiveness. Then the SNR improvement on proton spin density (PSD) and R2⁎ maps was investigated using brain STAGE data acquired from 10 healthy controls (HCs) and 10 patients with Parkinson's disease (PD). For the PSD and R2* maps, the SNR and CNR between white matter and gray matter were improved by a factor of 1.87 ± 0.50 and 1.72 ± 0.88, respectively. The white matter hyperintensity lesions in PD patients were more clearly defined after CROWN processing. Using these improved maps, simulated images for any repeat time, echo time or flip angle can be created with improved SNR. The potential applications of this technology are to trade off the increased SNR for higher resolution images and/or faster imaging.

hMRI - A toolbox for quantitative MRI in neuroscience and clinical research

Neuroscience and clinical researchers are increasingly interested in quantitative magnetic resonance imaging (qMRI) due to its sensitivity to micro-structural properties of brain tissue such as axon, myelin, iron and water concentration. We introduce the hMRI-toolbox, an open-source, easy-to-use tool available on GitHub, for qMRI data handling and processing, presented together with a tutorial and example dataset. This toolbox allows the estimation of high-quality multi-parameter qMRI maps (longitudinal and effective transverse relaxation rates R1 and R2⋆, proton density PD and magnetisation transfer MT saturation) that can be used for quantitative parameter analysis and accurate delineation of subcortical brain structures. The qMRI maps generated by the toolbox are key input parameters for biophysical models designed to estimate tissue microstructure properties such as the MR g-ratio and to derive standard and novel MRI biomarkers. Thus, the current version of the toolbox is a first step towards in vivo histology using MRI (hMRI) and is being extended further in this direction. Embedded in the Statistical Parametric Mapping (SPM) framework, it benefits from the extensive range of established SPM tools for high-accuracy spatial registration and statistical inferences and can be readily combined with existing SPM toolboxes for estimating diffusion MRI parameter maps. From a user's perspective, the hMRI-toolbox is an efficient, robust and simple framework for investigating qMRI data in neuroscience and clinical research.

Quantitative T1 and proton density mapping with direct calculation of radiofrequency coil transmit and receive profiles from two-point variable flip angle data

Quantitative T1 mapping of brain tissue is frequently based on the variable flip angle (VFA) method, acquiring spoiled gradient echo (GE) datasets at different excitation angles. However, accurate T1 calculation requires a knowledge of the sensitivity profile B1 of the radiofrequency (RF) transmit coil. For an additional derivation of proton density (PD) maps, the receive coil sensitivity profile (RP) must also be known. Mapping of B1 and RP increases the experiment duration, which may be critical when investigating patients. In this work, a method is presented for the direct calculation of B1 and RP from VFA data. Thus, quantitative maps of T1 , PD, B1 and RP can be obtained from only two spoiled GE datasets. The method is based on: (1) the exploitation of the linear relationship between 1/PD and 1/T1 in brain tissue and (2) the assumption of smoothly varying B1 and RP, so that a large number of data points can be fitted across small volume elements where B1 and RP are approximately constant. The method is tested and optimized on healthy subjects. Copyright © 2016 John Wiley & Sons, Ltd.

A general linear relaxometry model of R1 using imaging data

PURPOSE

The longitudinal relaxation rate (R1 ) measured in vivo depends on the local microstructural properties of the tissue, such as macromolecular, iron, and water content. Here, we use whole brain multiparametric in vivo data and a general linear relaxometry model to describe the dependence of R1 on these components. We explore a) the validity of having a single fixed set of model coefficients for the whole brain and b) the stability of the model coefficients in a large cohort.

METHODS

Maps of magnetization transfer (MT) and effective transverse relaxation rate (R2 *) were used as surrogates for macromolecular and iron content, respectively. Spatial variations in these parameters reflected variations in underlying tissue microstructure. A linear model was applied to the whole brain, including gray/white matter and deep brain structures, to determine the global model coefficients. Synthetic R1 values were then calculated using these coefficients and compared with the measured R1 maps.

RESULTS

The model's validity was demonstrated by correspondence between the synthetic and measured R1 values and by high stability of the model coefficients across a large cohort.

CONCLUSION

A single set of global coefficients can be used to relate R1 , MT, and R2 * across the whole brain. Our population study demonstrates the robustness and stability of the model.

Contrast enrichment of spinal cord MR imaging using a ratio of T1-weighted and T2-weighted signals

PURPOSE

We aimed to assess if the T1-weighted (T1w)/T2-weighted (T2w) signal ratio could be used to improve image contrast in MR spinal cord imaging.

MATERIALS AND METHODS

T1w and T2w cervical spinal cord MR images were acquired from 23 normal subjects using 3 Tesla (T) MR scanner. In addition, a multiple sclerosis patient, and a cervical spondylotic myelopathy patient were evaluated. White matter (WM) and gray matter (GM) signal intensities were measured for each image (T1w, T2w, and T1w/T2w) for seven cervical segments in each subject to calculate the contrast. Age-related changes in signal intensity were assessed at each location (lateral column, anterior column, dorsal column, and GM) for each image. Additionally, the imaging results of two subjects with spinal diseases and the controls were numerically compared.

RESULTS

The contrast between the WM and GM in the T1w/T2w ratio image was approximately twice as much as that in the T1w and T2w images (mean ± SD = 1.8 ± 0.4). The signal intensity ratio was related to age. For both clinical patients, the signal intensities were significantly lower in the lesion areas in the ratio images.

CONCLUSION

The T1w/T2w ratio images demonstrated increased image contrast compared with T1w and T2w images alone and, reduced inter-individual signal intensity differences.

Human skin permittivity determined by millimeter wave reflection measurements

Millimeter wave reflection from the human skin was studied in the frequency range of 37-74 GHz in steps of 1 GHz. The forearm and palm data were used to model the skin with thin and thick stratum corneum (SC), respectively. To fit the reflection data, a homogeneous unilayer and three multilayer skin models were tested. Skin permittivity in the mm-wave frequency range resulted from the permittivity of cutaneous free water which was described by the Debye equation. The permittivity increment found from fitting to the experimental data was used for determination of the complex permittivity and water content of skin layers. Our approach, first tested in pure water and gelatin gels with different water contents, gave good agreement with literature data. The homogeneous skin model fitted the forearm data well. Permittivity of the forearm skin obtained with this model was close to the skin permittivity reported by others. To fit reflection from the palmar skin with a thick SC, a skin model containing at least two layers was required. Multilayer models provided better fitting to both the forearm and palmar skin reflection data. The fitting parameters obtained with different models were consistent with each other.

Three-pool model of white matter

PURPOSE

To investigate the use of a three-pool relaxation model to measure myelin, myelinated-axon, and mixed water-pool fractions in white matter (WM) during myelination.

MATERIALS AND METHODS

MRI at 1.9 Tesla, and conventional spin-echo imaging were used to acquire T1 and T2 relaxation data in 15 normal children ranging in age from 3 months to 13 years 4 months. Three equations with three unknowns were solved to calculate three water-pool fractions for each child in a frontal association-fiber area and a frontal-parietal projection-fiber area. The temporal trend of the fractions was compared with a theoretical three-pool myelination model.

RESULTS

The myelin level in the projection-fiber area rose earlier than in the association-fiber area following the standard caudal-to-rostral trend. The temporal trend of the three-pool fractions followed that predicted by the theoretical myelination model in both brain areas. The myelinated-axon and mixed pool sizes were significantly different in the two WM areas following early myelination, although their myelin pools were similar. T1 values correlated more highly with the myelinated-axon and mixed pool fractions than with the myelin pool fraction.

CONCLUSION

The three-pool relaxation model provides measurements of water-pool fractions in WM that follow values predicted during myelination.

Interregional variation of longitudinal relaxation rates in human brain at 3.0 T: relation to estimated iron and water contents

In a study of interregional variation of the longitudinal relaxation rate (R(1)) in human brain at 3 T, R(1) maps were acquired from 12 healthy adults using a multi-slice implementation of the T one by multiple readout pulses (TOMROP) sequence. Mean R(1) values were obtained from the prefrontal cortex (0.567 +/- 0.020 sec(-1)), caudate head (0.675 +/- 0.019 sec(-1)), putamen (0.749 +/- 0.023 sec(-1)), substantia nigra (0.873 +/- 0.037 sec(-1)), globus pallidus (0.960 +/- 0.034 sec(-1)), thalamus (0.822 +/- 0.027 sec(-1)), and frontal white matter (1.184 +/- 0.057 sec(-1)). For gray matter regions other than the thalamus, R(1) showed a strong correlation (r = 0.984, P < 0.0001) with estimated regional nonheme iron concentrations ([Fe]). These R(1) values also showed a strong correlation (r = 0.976, P < 0.0001) with estimates of 1/f(w) obtained from MRI relative proton density measurements, where f(w) represents tissue water content. When white matter is included in the consideration, 1/f(w) is a better predictor of R(1) than is [Fe]. An analysis based on the fast-exchange two-state model of longitudinal relaxation suggests that interregional differences in f(w) account for the majority of the variation of R(1) across gray matter regions. Magn Reson Med 45:71-79, 2001.

An absolute measurement of brain water content using magnetic resonance imaging in two focal cerebral ischemic rat models

Magnetic resonance imaging (MRI) was utilized to obtain absolute estimates of regional brain water content (W), and results were compared with those obtained with conventional wet/dry measurements. In total, 31 male Long-Evans rats were studied and divided into two groups based on the surgical procedures used to induce cerebral focal ischemia: suture (n = 18) and three-vessel ligation (TVL: n = 13) groups. Both relative spin density and T1 were extracted from the acquired MR images. After correcting for radiofrequency field inhomogeneities, T2* signal decay, and temperature effects, in vivo regional brain water content, in absolute terms, was obtained by normalizing the measured relative brain spin density of animals to that of a water phantom. A highly linear relationship between MR-estimated brain water content based on the normalized spin density and wet/dry measurements was obtained with slopes of 0.989 and 0.986 for the suture (r = 0.79) and TVL (r = 0.83) groups, respectively. Except for the normal subcortex of the TVL group (P < 0.02) and the normal hemisphere of the suture group (P < 0.003), no significant differences were observed between MR-estimated and wet/dry measurements of brain water content. In addition, a highly linear relationship between MR-measured R1 (= 1/T1) and 1/W of wet/dry measurements was obtained. However, slopes of the linear regression lines in the two groups were significantly different (P < 0.02), indicating that different R1 values were associated with the same water content depending on the model. These results show that an absolute measurement of in vivo regional brain water content can be obtained with MRI and potentially serves as a noninvasive means to monitor different therapeutic interventions for the management of brain edema subsequent to stroke and head trauma.

Absolute measurements of water content using magnetic resonance imaging: preliminary findings in an in vivo focal ischemic rat model

Using a magnetic resonance (MR) imaging method, absolute measurements of in vivo brain water content were obtained in 15 male Long Evans rats that underwent a 90-min focal cerebral ischemia. A strong linear relationship (r = 0.80) with a slope of 1 was observed when correlating MR-measured water content to that obtained with the ex vivo wet/dry measurements. This effective spin-density-based method is the first-reported in vivo absolute quantification of brain tissue water content associated with a pathophysiological state and preliminary findings suggest that a noninvasive measurement of brain water content can be obtained with MRI.

T1 and magnetization transfer at 7 Tesla in acute ischemic infarct in the rat

T1 and magnetization transfer at a field strength of 7 Tesla were used to discriminate between water accumulation and protein mobilization in tissue undergoing infarction. Twelve rats subjected to acute stroke via intralumenal suture occlusion of the middle cerebral artery, and 19 controls, were studied. In MRI studies to 6 hr post-ictus, serial data acquisition allowed the measurement of cerebral blood flow (CBF), apparent diffusion coefficient of water (ADCw), equilibrium magnetization (M0) and T1, and equilibrium magnetization and T1 under an off-resonance partial saturation of the macromolecular pool (Msat and T1sat). Using these parameters, the apparent forward transfer rate of magnetization between the free water proton pool and the macromolecular proton pool, k(fa), was calculated. Regions of interest (ROIs) were chosen using depressed areas in maps of the ADCw. T1 measurements in bovine serum albumin at 7T were not affected by the mobility of the macromolecular pool (P > 0.2), but magnetization transfer between free water and protein depended strongly on the mobility of the macromolecular pool (P < 0.001). For 6 hr after ictus, k(fa) uniformly and strongly decreased in the region of the infarct (P < 0.0001). Ratios (ischemic/non-ischemic) of parameters M0, Msat, T1, and T1sat all uniformly and strongly increased in the infarct. The ratio T1/T1sat in the region of infarction showed that a progressive accumulation of free water in the region of interest was the major (>80%) contribution to the decrease in k(fa). There also existed a small contribution due to changes at the water-macromolecular interface, possibly due to proteolysis (P = 0.005).

Use of magnetic resonance imaging for in vivo measurements of water content in human brain: method and normal values

OBJECT

The authors present a quantitative in vivo magnetic resonance (MR) imaging method and propose its use for the accurate assessment of brain water in humans.

METHODS

With this technique, a pure T1-weighted image of a selected brain slice in a patient is generated, and the image is subsequently converted to a pure water image by means of an equation derived from a tissue relaxation model. The image intensity in the resulting water map directly yields absolute measures of water expressed in grams of water per gram of tissue at a given anatomical location. The method has been validated previously in a series of phantom experiments and in an infusion model of brain edema in cats. In this report, the authors evaluate the method by using samples of tissue harvested from patients who underwent surgery for brain tumor removal and apply the technique to a series of normal volunteers, providing average regional brain water content (f(w)) values for a range of tissues. Application of the method in pathological conditions such as head trauma, tumor, and hydrocephalus allows quantification of regional or global increases in f(w) that result from edema.

CONCLUSIONS

It is now possible to obtain accurate brain water measurements with the anatomical resolution of MR imaging. This permits monitoring of the development and resolution of edema in a variety of clinical circumstances, thus enhancing understanding of the underlying pathophysiological processes.

Quantitative regional brain water measurement with magnetic resonance imaging in a focal ischemia model

Therapeutic approaches to cerebral edema require an understanding of both the magnitude and location of changes in brain water content. It is desirable to have a sensitive, accurate means of measuring brain water noninvasively so that effective therapies for cerebral edema in stroke, head trauma, and other conditions can be investigated. In this work, a three-dimensional magnetic resonance imaging technique that is able to provide both spin density and T1 simultaneously is described. This method was used to quantitate regional changes in brain water content in a rat model of focal cerebral ischemia. Brain water contents estimated from both relative spin density and relative T1 measurements made in vivo were compared with ex vivo measurements of relative tissue water content based on the wet-dry technique. Correlation coefficients of 0.95 and 0.98 were obtained between the wet-dry measurements and magnetic resonance measurements of T1 and spin density, respectively. Notably, the slope of the relationship between T1 and tissue water content changed dramatically after the injection of a paramagnetic contrast agent while precontrast and postcontrast spin density measurements remained essentially invariant. In addition, a plot of absolute spin density (obtained by normalizing spin density from agar gelatin phantoms of different water contents to the spin density of a sample of 100% water) was linearly related to wet-dry measurements with a slope of 0.99 (R2 = 0.99).

In vivo measurement of T2 distributions and water contents in normal human brain

Using a 32-echo imaging pulse sequence, T2 relaxation decay curves were acquired from five white- and six gray-matter brain structures outlined in 12 normal volunteers. The water contents of white and gray matter were 0.71 (0.01) and 0.83 (0.03) g/ml, respectively. All white-matter structures had significantly higher myelin water percentages (signal percentage with T2 between 10 and 50 ms) than all gray-matter structures. The range in geometric mean T2 of the main peak for both white and gray matter was from 70 to 86 ms. T2 distributions from the posterior internal capsules and splenium of the corpus callosum were significantly wider (width is related to water environment inhomogeneity) than those from any other white- or gray-matter structures. Thus, quantitative measurement and analysis of T2 relaxation reveals differences in brain tissue water environments not discernible on conventional MR images. These differences may make short T2 components reliable markers for normal myelin.

In vivo visualization of myelin water in brain by magnetic resonance

We exploit the intrinsic difference in magnetic resonance spin-spin relaxation time, T2, between water associated with myelin sheaths and water in other central nervous system tissue in order to measure myelin water content within any region of an image or to generate indirectly a myelin map of the brain. In normal volunteers, myelin water maps give the expected myelin distribution. In multiple sclerosis patients, lesions exhibit different myelin water contents providing insight into the demyelination process unavailable from conventional magnetic resonance images. In vivo myelin measurement has important applications in the clinical management of multiple sclerosis and other white matter diseases.

In vivo brain water determination by T1 measurements: effect of total water content, hydration fraction, and field strength

This work is concerned with the accurate quantification of brain water content under routine clinical conditions. Gelatin solutions of varying water content are first employed as a model of an edematous brain and longitudinal relaxation measurements are performed at proton Larmor frequencies of 5, 41, 63, and 100 MHz. These are followed with in vivo measurements in an experimental animal model of brain edema at 41 MHz. The results underscore the dominant role of total water content W in the relaxation process and verify the expected linearity between 1/T1 and 1/W. A scheme is presented and experimentally verified at 41 MHz for deducing the exact relationship of 1/T1 vs 1/W at any frequency. Knowledge of this relationship along with precise measurements of 1/T1 at a given field strength permits quantitative in vivo measures of brain water content to be obtained with a precision of about 0.01. It is concluded that routine, accurate, and noninvasive brain water measurements are possible by magnetic resonance imaging in a clinical environment.

Investigation of cerebral ischemia using magnetization transfer contrast (MTC) MR imaging

The effects of cerebral ischemia in rat brain were monitored as a function of time using proton MR imaging. Spin-spin relaxation time (T2), proton density, and magnetization transfer contrast (MTC) were measured by MR imaging at various time intervals during a 1-week period following the induction of ischemic damage. Ischemic injury was characterized by a maximization of both T2 value and MTC appearance at 24 hr postischemic injury. These changes were accompanied by a gradual increase in MR observable water density over the first few days of ischemia. A reduction in the magnetization exchange rate between "free" and "bound" water protons as measured by MTC imaging is at least partially responsible for the elevation in T2 values observed during ischemia, and may accompany breakdown of cellular structure.

Experimental studies for use of magnetic resonance in brain water measurements

The accurate description and quantification of altered brain water resulting from different pathologic conditions is of critical clinical importance. In this work we determined the influence of total water content, hydration fraction and magnetic field strength on observed proton relaxation rates by means of in vitro and in vitro model studies and developed a scheme for determining water content at any field strength. Equations relating T1 relaxation times and brain water content are derived. This allows a non-invasive measure of brain water to be determined in the clinical setting.

Two-site exchange revisited: a new method for extracting exchange parameters in biological systems

A new analysis is presented which links real volume fractions, relaxation rates, and intracompartmental exchange rates directly with apparent volume fractions and relaxation rates obtained from biexponential fits of transverse magnetization decay curves. The analysis differs from previous methods in that measurements from two paramagnetic doping levels are used to close the two-site exchange equations. Both the new method and one previously described by Herbst and Goldstein (HG) have been applied to paramagnetically doped whole-blood data sets. Significant differences in the calculated exchange parameters are found between the two methods. A small dependence of the intracellular relaxation rate on extracellular paramagnetic agent concentration, assumed nonexistent with the HG method, is inferred from the new analysis. The analysis was also applied to published data on perfused rat hearts, and we obtained a limited assessment of two-site exchange in this system.