Using flow relaxography to elucidate flow relaxivity

Magnetic resonance multitasking for motion-resolved quantitative cardiovascular imaging

Quantitative cardiovascular magnetic resonance (CMR) imaging can be used to characterize fibrosis, oedema, ischaemia, inflammation and other disease conditions. However, the need to reduce artefacts arising from body motion through a combination of electrocardiography (ECG) control, respiration control, and contrast-weighting selection makes CMR exams lengthy. Here, we show that physiological motions and other dynamic processes can be conceptualized as multiple time dimensions that can be resolved via low-rank tensor imaging, allowing for motion-resolved quantitative imaging with up to four time dimensions. This continuous-acquisition approach, which we name cardiovascular MR multitasking, captures — rather than avoids — motion, relaxation and other dynamics to efficiently perform quantitative CMR without the use of ECG triggering or breath holds. We demonstrate that CMR multitasking allows for T₁ mapping, T₁-T₂ mapping and time-resolved T₁ mapping of myocardial perfusion without ECG information and/or in free-breathing conditions. CMR multitasking may provide a foundation for the development of setup-free CMR imaging for the quantitative evaluation of cardiovascular health.

Intracellular water preexchange lifetime in neurons and astrocytes

PURPOSE

To determine the intracellular water preexchange lifetime, τ_i , the "average residence time" of water, in the intracellular milieu of neurons and astrocytes. The preexchange lifetime is important for modeling a variety of MR data sets, including relaxation, diffusion-sensitive, and dynamic contrast-enhanced data sets.

METHODS

Herein, τ_i in neurons and astrocytes is determined in a microbead-adherent, cultured cell system. In concert with thin-slice selection, rapid flow of extracellular media suppresses extracellular signal, allowing determination of the transcytolemmal-exchange-dominated, intracellular T₁ . With this knowledge, and that of the intracellular T₁ in the absence of exchange, τ_i can be derived.

RESULTS

Under normal culture conditions, τ_i for neurons is 0.75 ± 0.05 s versus 0.57 ± 0.03 s for astrocytes. Both neuronal and astrocytic τ_i s decrease within 30 min after the onset of oxygen-glucose deprivation, with the astrocytic τ_i showing a substantially greater decrease than the neuronal τ_i .

CONCLUSIONS

Given an approximate intra- to extracellular volume ratio of 4:1 in the brain, these data imply that, under normal physiological conditions, an MR experimental characteristic time of less than 0.012 s is required for a nonexchanging, two-compartment (intra- and extracellular) model to be valid for MR studies. This characteristic time shortens significantly (i.e., 0.004 s) under injury conditions. Magn Reson Med 79:1616-1627, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Diffusion effects on longitudinal relaxation in poorly mixed compartments

Diffusion of spins between physical or virtual, communicating compartments having different states of longitudinal magnetization leads to diffusion-driven longitudinal relaxation. Herein, in two model systems, the effects of diffusion-driven longitudinal relaxation are explored experimentally and analyzed quantitatively. In the first case, longitudinal relaxation in a single slice of a water phantom is monitored spectroscopically as a function of slice thickness. In the second case, mimicking vascular flow/diffusion effects, longitudinal relaxation is monitored in a two-compartment, semi-permeable fiber phantom. In both cases, apparent longitudinal relaxation, though clearly multi-exponential, is well-modeled as bi-exponential.

Three-dimensional T1 mapping of the mouse heart using variable flip angle steady-state MR imaging

Cardiac MR T(1) mapping is a promising quantitative imaging tool for the diagnosis and evaluation of cardiomyopathy. Here, we present a new preclinical cardiac MRI method enabling three-dimensional T(1) mapping of the mouse heart. The method is based on a variable flip angle analysis of steady-state MR imaging data. A retrospectively triggered three-dimensional FLASH (fast low-angle shot) sequence (3D IntraGate) enables a constant repetition time and maintains steady-state conditions. 3D T(1) mapping of the complete mouse heart could be achieved in 20 min. High-quality, bright-blood T(1) maps were obtained with homogeneous T(1) values (1764 ± 172 ms) throughout the myocardium. The repeatability coefficient of R(1) (1/T(1) ) in a specific region of the mouse heart was between 0.14 and 0.20 s(-1) , depending on the number of flip angles. The feasibility for detecting regional differences in ΔR(1) was shown with pre- and post-contrast T(1) mapping in mice with surgically induced myocardial infarction, for which ΔR(1) values up to 0.83 s(-1) were found in the infarct zone. The sequence was also investigated in black-blood mode, which, interestingly, showed a strong decrease in the apparent mean T(1) of healthy myocardium (905 ± 110 ms). This study shows that 3D T(1) mapping in the mouse heart is feasible and can be used to monitor regional changes in myocardial T(1), particularly in relation to pathology and in contrast-enhanced experiments to estimate local concentrations of (targeted) contrast agent.

Quantification of slow flow using FAIR

Phase contrast (PC)-based MRI methods are considered to be the most accurate approach for spatially resolved flow quantification, but the measurement of very slow velocities requires signal detection at long echo times and the application of strong field gradients. On the other hand, measurements based on time-of-flight or inflow effects can be conducted at short echo times and without flow-encoding gradients. A method for imaging flow at velocities of the order of 0.1 mm/s is presented and validated here. It consists of measuring the apparent spin-lattice relation rate (R(1)*) of the flowing fluid using magnetization preparation by alternating slice-selective and nonselective inversion pulses (FAIR or flow-sensitive alternating inversion recovery) and a fast gradient-echo detection sequence. This method is appropriate for the quantitative imaging of slow flow at low Reynolds numbers in fluids where the T(2) values are too short to allow sensitive flow measurements by phase contrast-based methods.

Intracellular water-specific MR of microbead-adherent cells: the HeLa cell intracellular water exchange lifetime

Quantitative characterization of the intracellular water (1)H MR signal from cultured cells will provide critical biophysical insight into the MR signal from tissues in vivo. Microbeads provide a robust immobilization substrate for the many mammalian cell lines that adhere to surfaces and also provide sufficient cell density for observation of the intracellular water MR signal. However, selective observation of the intracellular water MR signal from perfused, microbead-adherent mammalian cells requires highly effective suppression of the extracellular water MR signal. We describe how high-velocity perfusion of microbead-adherent cells results in short apparent (1)H MR longitudinal and transverse relaxation times for the extracellular water in a thin slice selected orthogonal to the direction of flow. When combined with a spin-echo pulse sequence, this phenomenon provides highly effective suppression of the extracellular water MR signal. This new method is exploited here to quantify the kinetics of water exchange from the intracellular to extracellular spaces of HeLa cells. The time constant describing water exchange from intracellular to extracellular spaces, also known as the exchange lifetime for intracellular water, is 119 +/- 14 ms.

Magnetic field and tissue dependencies of human brain longitudinal1H2O relaxation in vivo

Brain water proton (1H2O) longitudinal relaxation time constants (T1) were obtained from three healthy individuals at magnetic field strengths (B0) of 0.2 Tesla (T), 1.0T, 1.5T, 4.0T, and 7.0T. A 5-mm midventricular axial slice was sampled using a modified Look-Locker technique with 1.5 mm in-plane resolution, and 32 time points post-adiabatic inversion. The results confirmed that for most brain tissues, T1 values increased by more than a factor of 3 between 0.2T and 7T, and over this range were well fitted by T1 (s)=0.583(B0)0.382, T1(s)=0.857(B0)0.376, and T1(s)=1.35(B0)0.340 for white matter (WM), internal GM, and blood 1H2O, respectively. The ventricular cerebrospinal fluid (CSF) 1H2O T1 value did not change with B0, and its average value (standard deviation (SD)) across subjects and magnetic fields was 4.3 (+/-0.2) s. The tissue 1/T1 values at each field were well correlated with the macromolecular mass fraction, and to a lesser extent tissue iron content. The field-dependent increases in 1H2O T1 values more than offset the well-known decrease in typical MRI contrast reagent (CR) relaxivity, and simulations predict that this leads to lower CR concentration detection thresholds with increased magnetic field.

Evidence for shutter-speed variation in CR bolus-tracking studies of human pathology

The standard pharmacokinetic model for the analysis of MRI contrast reagent (CR) bolus-tracking (B-T) data assumes that the mean intracellular water molecule lifetime (tau(i)) is effectively zero. This assertion is inconsistent with a considerable body of physiological measurements. Furthermore, theory and simulation show the B-T time-course shape to be very sensitive to the tau(i) magnitude in the physiological range (hundreds of milliseconds to several seconds). Consequently, this standard model aspect can cause significant underestimations (factors of 2 or 3) of the two parameters usually determined: K(trans), the vascular wall CR transfer rate constant, and v(e), the CR distribution volume (the extracellular, extravascular space fraction). Analyses of animal model data confirmed two predicted behaviors indicative of this standard model inadequacy: (1) a specific temporal pattern for the mismatch between the best-fitted curve and data; and (2) an inverse dependence of the curve's K(trans) and v(e) magnitudes on the CR dose. These parameters should be CR dose-independent. The most parsimonious analysis allowing for realistic tau(i) values is the 'shutter-speed' model. Its application to the experimental animal data essentially eliminated the two standard model signature inadequacies. This paper reports the first survey for the extent of this 'shutter-speed effect' in human data. Retrospective analyses are made of clinical data chosen from a range of pathology (the active multiple sclerosis lesion, the invasive ductal carcinoma breast tumor, and osteosarcoma in the leg) that provides a wide variation, particularly of K(trans). The signature temporal mismatch of the standard model is observed in all cases, and is essentially eliminated by use of the shutter-speed model. Pixel-by-pixel maps show that parameter values from the shutter-speed analysis are increased by more than a factor of 3 for some lesion regions. This endows the lesions with very high contrast, and reveals heterogeneities that are often not seen in the standard model maps. Normal muscle regions in the leg allow validation of the shutter-speed model K(trans), v(e), and tau(i) magnitudes, by comparison with results of previous careful rat leg studies not possible for human subjects.

p38 mitogen-activated protein kinase modulates expression of tumor necrosis factor-related apoptosis-inducing ligand induced by interferon-gamma in fetal brain astrocytes

This study describes the involvement of the p38 mitogen-activated protein kinase (MAPK) during interferon-gamma (IFN-gamma) signaling in fetal brain astrocytes. In some pathological conditions of brain, p38 MAPK transduces stress-related signals, increases expression of proinflammatory cytokines, and induces cellular damage or apoptosis. In astrocytes, the tumor necrosis factor (TNF)-related apoptosis-inducing ligand (TRAIL) expression level was increased by IFN-gamma. AG490, a JAK inhibitor, blocked TRAIL expression induced by IFN-gamma. SB203580, a specific p38alpha and p38beta2 MAPK inhibitor, decreased the TRAIL expression induced by IFN-gamma. The phosphorylation of the Ser727 site of STAT1, but not the Tyr701 site, was inhibited by SB203580. These results suggest that p38 MAPK modulates STAT1 phosphorylation in IFN-gamma signaling in fetal brain astrocytes.

Parsimonious basis selection in exponential spectral analysis

Sums of decaying real exponentials (SDREs) are frequently used in models of time-varying processes. First-order compartmental models are widely employed to describe mass transit in chemical and biological systems. In these models the evolution of compartment concentration versus time is represented as the convolution of an input function with an SDRE. In exponential spectral analysis (ESA) the nonlinear problem of estimating the SDRE rate constants is replaced by the linear estimation of the coefficients of a preselected set of exponential basis functions (EBFs). This work addresses the problem of selecting the number of EBFs and the rate constant of each basis element. Basis dimension is established via model selection, in which approximation error and parameter redundancy are the criteria. The latter is estimated via simulation of the fitted model over multiple noise realizations. A constrained Cramér-Rao lower bound is derived for ESA parameters. The resulting parsimonious ESA algorithm (PESA) ameliorates the inherent problem of non-uniqueness in ESA parameters. Consequently, sets of time series may be compared in a statistically meaningfully way in terms of physically or physiologically significant parameters. PESA is applied to compare the retention of two radiotracers in the artificially perfused rabbit heart.

A model of blood-brain barrier permeability to water: accounting for blood inflow and longitudinal relaxation effects

A noninvasive technique for measuring the permeability of the blood-brain barrier (BBB) to water could help to evaluate changes in the functional integrity of the BBB that occur in different pathologies, such as multiple sclerosis or growth of brain tumor. Recently, Schwarzbauer et al. (Magn Reson Med 1997;37:769-777) proposed an MR method to measure this permeability based on the T(1) reductions induced by injecting various doses of paramagnetic contrast agent. However, this method may be difficult to implement in a clinical environment. Described here is a two-point technique, in which a spatially selective inversion is used to measure T(1) prior to and after injection of an intravascular contrast agent. Measurements made in the rat brain are compared to numerical simulations generated with a physiological model that accounts for blood flow and includes two different blood volumes: nonexchanging and exchanging blood volumes. Our results suggest that BBB permeability could be evaluated from the change in T(1) caused by the vascular contrast agent. This technique might provide an approach for monitoring changes in BBB permeability to water in clinical studies.

Two-dimensional time correlation relaxometry of skeletal muscle in vivo at 3 Tesla

A hybrid two-dimensional relaxometry (2DR) sequence was used to simultaneously measure both the spin-spin (R2) and spin-lattice relaxation rates (R1) of skeletal muscle in vivo. The 2DR sequence involved a 180 degrees inversion pulse followed by a variable delay time (30 values from 40 to 7000 ms); a projection presaturation (PP) scheme to localize a 16-ml cylindrical voxel; and a CPMG sequence (950 even echoes, effective echo spacing = 1.2 ms, equilibrium time = 12 s). The 2DR data were collected at 3.0 Tesla from the flexor digitorum profundus of eight healthy males, 26 +/- 2 years old. Analysis was performed with a 2D version of the non-negative least-squares algorithm and a one-way ANOVA. All subjects exhibited at least three spin-groups (R2 < 200 s(-1)), designated B, C, and D, with R2 values of 42.7, 26.5, and 8.1 s(-1), and fractional volumes of 52, 35, and 7%, respectively. The R1 values of B and C were similar, congruent with0.7 s(-1), but different from that of D (P < 0.001), which had an R1 of 1.0 s(-1). The results suggest that exchange between B and C ranges from 0.7-16.2 s(-1), while exchange between either of these spin-groups with D is slower. If the data are interpreted with a compartment model, in which spin-groups with short and long R2 values are attributed to extra- and intracellular fluid, respectively, the exchange of water across the cell membrane in living skeletal muscle is slow or intermediate relative to both R1 and R2.

Determination of the MRI contrast agent concentration time course in vivo following bolus injection: effect of equilibrium transcytolemmal water exchange

For bolus-tracking studies, it is commonly assumed that CR concentration bears a linear relationship with the measured (usually longitudinal) (1)H(2)O relaxation rate constant, R*(1) identical with(T(1) *)(-1). This requires that equilibrium transcytolemmal water exchange be in the fast exchange limit (FXL). However, though systems remain in fast exchange, the FXL will not usually obtain. Here, the consequences are considered: 1) the measurement of R(1) * itself can be affected, 2) the resultant non-linear [CR]-dependence causes significant error by assuming FXL, 3) the thermodynamic [CR] (based on the space in which CR is actually distributed) can be determined, 4) transcytolemmal water permeability may be estimated, and 5) the pharmacokinetic parameters can be factored. For a 30-sec, 0.17 mmol/kg dose of GdDTPA(2-), the FXL assumption underestimates the [CR] maximum in rat thigh muscle by a factor of almost two. Similar results are obtained for a rat brain GS-9L gliosarcoma tumor model.

A model system for perfusion quantification using FAIR

Flow-sensitive experiments (FAIR) have been performed on a tube-flow phantom in order to validate quantitative perfusion measurements on humans. A straight-forward correspondence between perfusion and bulk-flow is found. It is shown that the flow phantom model only holds when the slice profiles of the involved RF pulses are taken into account. A small flow-independent off-set may be present in the data. The off-set is explained by the model. Based on the correspondence between the phantom and the in vivo models, it is shown that the lowest flow values that could be measured in the phantom correspond to perfusion values lower than the cortical perfusion in the brain. Thus, the experimental accuracy and the computational methods for quantitative perfusion measurements in vivo can be validated by a tube-flow phantom.

PURR-TURBO: a novel pulse sequence for longitudinal relaxographic imaging

A novel pulse sequence based on a segmented phase-encoding scheme for measuring the longitudinal relaxation time (T(1)) value of an NMR signal in an imaging context is introduced. This pulse sequence is a hybrid version of the "single-shot" and "one-shot" inversion recovery (IR) snapshot-FLASH methods. These are also known as "multipoint" IR techniques. The new sequence presented here collects multiple k-space rows at each time point during magnetization recovery, rather than the entire dataset (as in a "single-shot" method) or just one row (as in a "one-shot" method). Thus, it reduces the scanning time without significant sacrifice of the small sampling time advantage of the one-shot IR snapshot-FLASH methods. Furthermore, this approach does not require a high-performance gradient system. Here, we demonstrate that a single slice human brain (1)H(2)O T(1) map with a nominal in-plane resolution of less than (1 mm)2can be obtained at 4 T in about 4 min. Published 2000 Wiley-Liss, Inc.

Equilibrium transcytolemmal water-exchange kinetics in skeletal muscle in vivo

It is commonly assumed that equilibrium transcytolemmal water exchange in tissue is sufficiently frequent as to be fast on any NMR time scale achievable with an extracellular contrast agent (CR) in vivo. A survey of literature values for cell membrane diffusional permeability coefficients (P) and cell sizes suggests that this should not really be so. To evaluate this issue experimentally, we used a programmed intravenous CR infusion protocol for the rat with several rate plateaus, each of which achieved an increased steady-state concentration of GdDTPA(2-) in the blood plasma. Interleaved rigorous measurements of (1)H(2)O inversion recoveries were made from arterial blood and from a region of homogeneous thigh muscle tissue throughout the CR infusion. We made careful relaxographic analyses for the blood and muscle (1)H(2)O longitudinal relaxation times. The combined data from several animals were evaluated with a two-site model for equilibrium transcytolemmal water exchange. An excellent fitting was achieved, with parameters that agreed very well with the relevant physiological properties available in the literature. The fraction of water in the extracellular space, 0.11, is quite consistent with published values, as well as with reported tissue CR concentrations when one accounts for the restriction of CR to this space. The derived average lifetime for a water molecule in the thigh muscle sarcoplasm, 1.1 +/- 0.4 sec, implies a sarcolemmal P of 13 x 10(-4) cm/sec, which is well within the range of literature values determined in vitro. Moreover, we find that because of the exchange, the (1)H(2)O longitudinal relaxation rate constant exhibits a decided nonlinear dependence on the tissue or thermodynamic (extracellular) concentration of GdDTPA(2-). The muscle system departs the fast-exchange limit at a [CR] value of <100 micromol/L. This has significant implications for the quantitative use of CRs as MRI tracers. Magn Reson Med 42:467-478, 1999. Published 1999 Wiley-Liss, Inc.

Intimate combination of low- and high-resolution image data: I. Real-space PET and (1)H(2)O MRI, PETAMRI

Two different types of (co-registered) images of the same slice of tissue will generally have different spatial resolutions. The judicious pixel-by-pixel combination of their data can be accomplished to yield a single image exhibiting properties of both. Here, axial (18)FDG PET and (1)H(2)O MR images of the human brain are used as the low- and high-resolution members of the pair. A color scale is necessary in order to provide for separate intensity parameters from the two image types. However, not all color scales can accommodate this separability. The HSV color model allows one to choose a color scale in which the intensity of the low-resolution image type is coded as hue, while that of the high-resolution type is coded as value, a reasonably independent parameter. Furthermore, the high-resolution image must have high contrast and be quantitative in the same sense as the low-resolution image almost always is. Here, relaxographic MR images (naturally segmented quantitative (1)H(2)O spin-density components) are used. Their essentially complete contrast serves to effect an apparent editing function when encoded as the value of the color scale. Thus, the combination of (18)FDG PET images with gray-matter (GM) relaxographic (1)H(2)O images produces visually "GM-edited" (18)FDG PETAMR (positron emission tomography and magnetic resonance) images. These exhibit the high sensitivity to tracer amounts characteristic of PET along with the high spatial resolution of (1)H(2)O MRI. At the same time, however, they retain the complete quantitative measures of each of their basis images. Magn Reson Med 42:345-360, 1999. Published 1999 Wiley-Liss, Inc.