Characterization of the diffusion coefficient of blood

SpinDoctor-IVIM: A virtual imaging framework for intravoxel incoherent motion MRI

Intravoxel incoherent motion (IVIM) imaging is increasingly recognised as an important tool in clinical MRI, where tissue perfusion and diffusion information can aid disease diagnosis, monitoring of patient recovery, and treatment outcome assessment. Currently, the discovery of biomarkers based on IVIM imaging, similar to other medical imaging modalities, is dependent on long preclinical and clinical validation pathways to link observable markers derived from images with the underlying pathophysiological mechanisms. To speed up this process, virtual IVIM imaging is proposed. This approach provides an efficient virtual imaging tool to design, evaluate, and optimise novel approaches for IVIM imaging. In this work, virtual IVIM imaging is developed through a new finite element solver, SpinDoctor-IVIM, which extends SpinDoctor, a diffusion MRI simulation toolbox. SpinDoctor-IVIM simulates IVIM imaging signals by solving the generalised Bloch-Torrey partial differential equation. The input velocity to SpinDoctor-IVIM is computed using HemeLB, an established Lattice Boltzmann blood flow simulator. Contrary to previous approaches, SpinDoctor-IVIM accounts for volumetric microvasculature during blood flow simulations, incorporates diffusion phenomena in the intravascular space, and accounts for the permeability between the intravascular and extravascular spaces. The above-mentioned features of the proposed framework are illustrated with simulations on a realistic microvasculature model.

Spatial profiling of in vivo diffusion-weighted MRI parameters in the healthy human kidney

OBJECTIVE

Diffusion-weighted MRI is a technique that can infer microstructural and microcirculatory features from biological tissue, with particular application to renal tissue. There is extensive literature on diffusion tensor imaging (DTI) of anisotropy in the renal medulla, intravoxel incoherent motion (IVIM) measurements separating microstructural from microcirculation effects, and combinations of the two. However, interpretation of these features and adaptation of more specific models remains an ongoing challenge. One input to this process is a whole organ distillation of corticomedullary contrast of diffusion metrics, as has been explored for other renal biomarkers.

MATERIALS AND METHODS

In this work, we probe the spatial dependence of diffusion MRI metrics with concentrically layered segmentation in 11 healthy kidneys at 3 T. The metrics include those from DTI, IVIM, a combined approach titled "REnal Flow and Microstructure AnisotroPy (REFMAP)", and a multiply encoded model titled "FC-IVIM" providing estimates of fluid velocity and branching length.

RESULTS

Fractional anisotropy decreased from the inner kidney to the outer kidney with the strongest layer correlation in both parenchyma (including cortex and medulla) and medulla with Spearman correlation coefficients and p-values (r, p) of (0.42, <0.001) and (0.37, <0.001), respectively. Also, dynamic parameters derived from the three models significantly decreased with a high correlation from the inner to the outer parenchyma or medulla with (r, p) ranges of (0.46-0.55, <0.001).

CONCLUSIONS

These spatial trends might find implications for indirect assessments of kidney physiology and microstructure using diffusion MRI.

Influence of Magnetic Field Strength on Intravoxel Incoherent Motion Parameters in Diffusion MRI of the Calf

Background: The purpose of this study was to investigate the dependence of Intravoxel Incoherent Motion (IVIM) parameters measured in the human calf on B0. Methods: Diffusion-weighted image data of eight healthy volunteers were acquired using five b-values (0-600 s/mm2) at rest and after muscle activation at 0.55 and 7 T. The musculus gastrocnemius mediale (GM, activated) was assessed. The perfusion fraction f and diffusion coefficient D were determined using segmented fits. The dependence on field strength was assessed using Student's t-test for paired samples and the Wilcoxon signed-rank test. A biophysical model built on the three non-exchanging compartments of muscle, venous blood, and arterial blood was used to interpret the data using literature relaxation times. Results: The measured perfusion fraction of the GM was significantly lower at 7 T, both for the baseline measurement and after muscle activation. For 0.55 and 7 T, the mean f values were 7.59% and 3.63% at rest, and 14.03% and 6.92% after activation, respectively. The biophysical model estimations for the mean proton-density-weighted perfusion fraction were 3.37% and 6.50% for the non-activated and activated states, respectively. Conclusions: B0 may have a significant effect on the measured IVIM parameters. The blood relaxation times suggest that 7 T IVIM may be arterial-weighted whereas 0.55 T IVIM may exhibit an approximately equal weighting of arterial and venous blood.

Impact of velocity- and acceleration-compensated encodings on signal dropout and black-blood state in diffusion-weighted magnetic resonance liver imaging at clinical TEs

PURPOSE

The study aims to develop easy-to-implement concomitant field-compensated gradient waveforms with varying velocity-weighting (M1) and acceleration-weighting (M2) levels and to evaluate their efficacy in correcting signal dropouts and preserving the black-blood state in liver diffusion-weighted imaging. Additionally, we seek to determine an optimal degree of compensation that minimizes signal dropouts while maintaining blood signal suppression.

METHODS

Numerically optimized gradient waveforms were adapted using a novel method that allows for the simultaneous tuning of M1- and M2-weighting by changing only one timing variable. Seven healthy volunteers underwent diffusion-weighted magnetic resonance imaging (DWI) with five diffusion encoding schemes (monopolar, velocity-compensated (M1 = 0), acceleration-compensated (M1 = M2 = 0), 84%-M1-M2-compensated, 67%-M1-M2-compensated) at b-values of 50 and 800 s/mm2 at a constant echo time of 70 ms. Signal dropout correction and apparent diffusion coefficients (ADCs) were quantified using regions of interest in the left and right liver lobe. The blood appearance was evaluated using two five-point Likert scales.

RESULTS

Signal dropout was more pronounced in the left lobe (19%-42% less signal than in the right lobe with monopolar scheme) and best corrected by acceleration-compensation (8%-10% less signal than in the right lobe). The black-blood state was best with monopolar encodings and decreased significantly (p < 0.001) with velocity- and/or acceleration-compensation. The partially M1-M2-compensated encoding schemes could restore the black-blood state again. Strongest ADC bias occurred for monopolar encodings (difference between left/right lobe of 0.41 μm2/ms for monopolar vs. < 0.12 μm2/ms for the other encodings).

CONCLUSION

All of the diffusion encodings used in this study demonstrated suitability for routine DWI application. The results indicate that a perfect value for the level of M1-M2-compensation does not exist. However, among the examined encodings, the 84%-M1-M2-compensated encodings provided a suitable tradeoff.

Nuclear magnetic resonance footprint of Wharton Jelly mesenchymal stem cells death mechanisms and distinctive in‐cell biophysical properties in vitro

The importance of the biophysical characterization of mesenchymal stem cells (MSCs) was recently pointed out for supporting the development of MSC‐based therapies. Among others, tracking MSCs in vivo and a quantitative characterization of their regenerative impact by nuclear magnetic resonance (NMR) demands a full description of MSCs’ MR properties. In the work, Wharton Jelly MSCs are characterized in a low magnetic field (LF) in vitro by using different approaches. They encompass various settings: MSCs cultured in a Petri dish and cell suspensions; experiments‐ 1D‐T₁, 1D‐T₂, 1D diffusion, 2D T₁‐T₂ and D‐T₂; devices‐ with a bore aperture and single‐sided one. Complex NMR analysis with the aid of random walk simulations allows the determination of MSCs T₁ and T₂ relaxation times, cells and nuclei sizes, self‐diffusion coefficients of the nucleus and cytoplasm. In addition, the influence of a single layer of cells on the effective diffusion coefficient of water is detected with the application of a single‐sided NMR device. It also enables the identification of apoptotic and necrotic cell death and changed diffusional properties of cells suspension caused by compressing forces induced by the subsequent cell layers. The study delivers MSCs‐specific MR parameters that may help tracking MSCs in vivo.

Non-Invasive Blood Flow Speed Measurement Using Optics

Non-invasive measurement of the arterial blood speed gives important health information such as cardio output and blood supplies to vital organs. The magnitude and change in arterial blood speed are key indicators of the health conditions and development and progression of diseases. We demonstrated a simple technique to directly measure the blood flow speed in main arteries based on the diffused light model. The concept is demonstrated with a phantom that uses intralipid hydrogel to model the biological tissue and an embedded glass tube with flowing human blood to model the blood vessel. The correlation function of the measured photocurrent was used to find the electrical field correlation function via the Siegert relation. We have shown that the characteristic decorrelation rate (i.e., the inverse of the decoherent time) is linearly proportional to the blood speed and independent of the tube diameter. This striking property can be explained by an approximate analytic solution for the diffused light equation in the regime where the convective flow is the dominating factor for decorrelation. As a result, we have demonstrated a non-invasive method of measuring arterial blood speed without any prior knowledge or assumption about the geometric or mechanic properties of the blood vessels.

Estimation of the Morphofunctional Status of the Brain in Hypertensive Wistar Rats Using Diffusion-Weighted MRI

Morphofunctional changes of the brain tissues of Wistar rats were studied based on the development of a multifactor cardiovasorenal model of arterial hypertension using MRI. An increase of the signal on the diffusion brain maps was recorded in 3 months, which indicated fluid accumulation in the intra- and extracellular space of the brain tissue. The data characterize the development of the pathogenetic mechanism of the hypervolemic variant of experimental arterial hypertension. The development of endothelial dysfunction in the brain vessels was manifested by predominance of abnormal constrictor reactions. In 6 months after arterial hypertension simulation, structural changes in the brain developed, such as leukoareosis, cystic encephalomalacia with dilated cerebrospinal fluid spaces and limited blood supply to brain tissue in the basins of the large cerebral arteries.

Effect of compression garments on muscle perfusion in delayed-onset muscle soreness: A quantitative analysis using intravoxel incoherent motion MR perfusion imaging

The aim of this prospective cohort study was to evaluate the effect of compression garments under resting conditions and after the induction of delayed-onset muscle soreness (DOMS) by MR perfusion imaging using intravoxel incoherent motion (IVIM). Magnetic resonance imaging of both lower legs of 16 volunteers was performed before and after standardized eccentric exercises that induced DOMS. A compression garment (21-22 mmHg) was worn during and for 6 h after exercise on one randomly selected leg. IVIM MR imaging, represented as total muscle perfusion D*f, perfusion fraction f and tissue diffusivity D, were compared between baseline and directly, 30 min, 6 h and 48 h after exhausting exercise with and without compression. Creatine kinase levels and T2-weighted images were acquired at baseline and after 48 h. DOMS was induced in the medial head of the gastrocnemius muscle (MGM) in all volunteers. Compression garments did not show any significant effect on IVIM perfusion parameters at any time point in the MGM or the tibialis anterior muscle (p > 0.05). Microvascular perfusion in the MGM increased significantly in both the compressed and noncompressed leg between baseline measurements and those taken directly after and 30 min after the exercise: the relative median f increased by 31.5% and 24.7% in the compressed and noncompressed leg, respectively, directly after the exercise compared with the baseline value. No significant change in tissue perfusion occurred 48 h after the induction of DOMS compared with baseline. It was concluded that compression garments (21-22 mmHg) do not alter microvascular muscle perfusion at rest, nor do they have any significant effect during the regeneration phase of DOMS. In DOMS, only a short-term effect of increased muscle perfusion (30 min after exercise) was observed, with normalization occurring during regeneration after 6-48 h. The normalization of perfusion independently of compression after 6 h may have implications for diagnostic and therapeutic strategies and for the better understanding of pathophysiological pathways in DOMS.

MR cell size imaging with temporal diffusion spectroscopy

Cytological features such as cell size and intracellular morphology provide fundamental information on cell status and hence may provide specific information on changes that arise within biological tissues. Such information is usually obtained by invasive biopsy in current clinical practice, which suffers several well-known disadvantages. Recently, novel MRI methods such as IMPULSED (imaging microstructural parameters using limited spectrally edited diffusion) have been developed for direct measurements of mean cell size non-invasively. The IMPULSED protocol is based on using temporal diffusion spectroscopy (TDS) to combine measurements of water diffusion over a wide range of diffusion times to probe cellular microstructure over varying length scales. IMPULSED has been shown to provide rapid, robust, and reliable mapping of mean cell size and is suitable for clinical imaging. More recently, cell size distributions have also been derived by appropriate analyses of data acquired with IMPULSED or similar sequences, which thus provides MRI-cytometry. This review summarizes the basic principles, practical implementations, validations, and example applications of MR cell size imaging based on TDS and demonstrates how cytometric information can be used in various applications. In addition, the limitations and potential future directions of MR cytometry are identified including the diagnosis of nonalcoholic steatohepatitis of the liver and the assessment of treatment response of cancers.

3.0 T MRI IVIM-DWI for predicting the efficacy of neoadjuvant chemoradiation for locally advanced rectal cancer

Purpose

The purpose of this study was to determine the diagnostic performance of intravoxel incoherent motion (IVIM) on assessing response to neoadjuvant chemoradiation (nCRT) in patients with Locally Advanced Rectal Cancer (LARC).

Methods

50 patients with rectal cancer who underwent magnetic resonance (MR) imaging before and after nCRT, the values of pre-nCRT and post-nCRT IVIM-DWI parameters apparent diffusion coefficient (ADC), diffusion coefficient (D), false diffusion coefficient (D*), and perfusion fraction (f), together with the percentage changes (∆% parametric value) induced by nCRT were calculated. According to the patient's response to nCRT, the patients were divided into pathological complete response (pCR) and non-pCR groups, Good Response (GR) group and Poor Response (PR) group, and the above values were compared between different groups. Univariate and multiple logistic regression analysis were done to investigate the relation between different parameters and patient nCRT. Draw ROC curve according to sensitivity and specificity, and compare its diagnostic efficacy.

Results

There were no significant differences in the baseline data of 50 patients. After nCRT, the ADC and D values for LARC increased significantly (all p < 0.05). The pCR group (n = 9) had higher preD*, pref, postD*, ∆%ADC and ∆%D values than the non-pCR group (n = 41) (all p < 0.05). The GR group (n = 17) exhibited higher post D, ∆%ADC and ∆%D values than the PR group (n = 33) (all p < 0.05). From the results of Logistic regression analysis found that ∆%ADC and ∆%D were significantly correlated with patients' response to nCRT. Based on ROC analysis, ∆%D had a higher area under the curve value than ∆%ADC (p = 0.009) in discriminating the pCR from non-pCR groups.

Conclusions

IVIM-DWI technology may be helpful in identifying the pCR and GR patients to nCRT for LARC.

The haemodynamics of the human placenta in utero

We have used magnetic resonance imaging (MRI) to provide important new insights into the function of the human placenta in utero. We have measured slow net flow and high net oxygenation in the placenta in vivo, which are consistent with efficient delivery of oxygen from mother to fetus. Our experimental evidence substantiates previous hypotheses on the effects of spiral artery remodelling in utero and also indicates rapid venous drainage from the placenta, which is important because this outflow has been largely neglected in the past. Furthermore, beyond Braxton Hicks contractions, which involve the entire uterus, we have identified a new physiological phenomenon, the ‘utero-placental pump’, by which the placenta and underlying uterine wall contract independently of the rest of the uterus, expelling maternal blood from the intervillous space.

MRI provides important new insights into the function of the human placenta, revealing slow net flow and high, uniform oxygenation in healthy pregnancies, detecting changes that will lead to compromised oxygen delivery to the fetus in preeclampsia, and identifying a new physiological phenomenon, the ‘utero-placental pump’.

Mapping hepatocyte size in vivo using temporal diffusion spectroscopy MRI

PURPOSE

The goal of this study is to implement a noninvasive method for in vivo mapping of hepatocyte size. This method will have a broad range of clinical and preclinical applications, as pathological changes in hepatocyte sizes are relevant for the accurate diagnosis and assessments of treatment response of liver diseases.

METHODS

Building on the concepts of temporal diffusion spectroscopy in MRI, a clinically feasible imaging protocol named IMPULSED (Imaging Microstructural Parameters Using Limited Spectrally Edited Diffusion) has been developed, which is able to report measurements of cell sizes noninvasively. This protocol acquires a selected set of diffusion imaging data and fits them to a model of water compartments in tissues to derive robust estimates of the cellular structures that restrict free diffusion. Here, we adapt and further develop this approach to measure hepatocyte sizes in vivo. We validated IMPULSED in livers of mice and rats and implemented it to image healthy human subjects using a clinical 3T MRI scanner.

RESULTS

The IMPULSED-derived mean hepatocyte sizes for rats and mice are about 15-20 µm and agree well with histological findings. Maps of mean hepatocyte size for humans can be achieved in less than 15 minutes, a clinically feasible scan time.

CONCLUSION

Our results suggest that this method has potential to overcome major limitations of liver biopsy and provide noninvasive mapping of hepatocyte sizes in clinical applications.

Optimal acquisition scheme for flow-compensated intravoxel incoherent motion diffusion-weighted imaging in the abdomen: An accurate and precise clinically feasible protocol

PURPOSE

Flow-compensated (FC) diffusion-weighted MRI (DWI) for intravoxel-incoherent motion (IVIM) modeling allows for a more detailed description of tissue microvasculature than conventional IVIM. The long acquisition time of current FC-IVIM protocols, however, has prohibited clinical application. Therefore, we developed an optimized abdominal FC-IVIM acquisition with a clinically feasible scan time.

METHODS

Precision and accuracy of the FC-IVIM parameters were assessed by fitting the FC-IVIM model to signal decay curves, simulated for different acquisition schemes. Diffusion-weighted acquisitions were added subsequently to the protocol, where we chose the combination of b-value, diffusion time and gradient profile (FC or bipolar) that resulted in the largest improvement to its accuracy and precision. The resulting two optimized FC-IVIM protocols with 25 and 50 acquisitions (FC-IVIMopt25 and FC-IVIMopt50 ), together with a complementary acquisition consisting of 50 diffusion-weighting (FC-IVIMcomp ), were acquired in repeated abdominal free-breathing FC-IVIM imaging of seven healthy volunteers. Intersession and intrasession within-subject coefficient of variation of the FC-IVIM parameters were compared for the liver, spleen, and kidneys.

RESULTS

Simulations showed that the performance of FC-IVIM improved in tissue with larger perfusion fraction and signal-to-noise ratio. The scan time of the FC-IVIMopt25 and FC-IVIMopt50 protocols were 8 and 16 min. The best in vivo performance was seen in FC-IVIMopt50 . The intersession within-subject coefficients of variation of FC-IVIMopt50 were 11.6%, 16.3%, 65.5%, and 36.0% for FC-IVIM model parameters diffusivity, perfusion fraction, characteristic time and blood flow velocity, respectively.

CONCLUSIONS

We have optimized the FC-IVIM protocol, allowing for clinically feasible scan times (8-16 min).

Performing IVIM-DWI using the multifunctional nanosystem for the evaluation of the antitumor microcirculation changes

OBJECTIVES

There is a controversy about the D* and f values of intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) for mid- and long-term efficacy monitoring of tumor blood perfusion. To monitor the antitumor efficacy of the F/A-PLGA@DOX/SPIO nanosystem via IVIM-DWI and to explore the value of parameters pseudo-diffusion (D*) and fraction of pseudo-diffusion (f) for evaluating therapeutic effect in non-small cell lung cancer (NSCLC).

MATERIALS AND METHODS

Thirty-six A549 tumor-bearing mice were divided randomly into three groups (each n = 12). Group 1 (G1) was injected with saline (the control group). Group 2 (G2) and group 3(G3) were injected with DOX and F/A-PLGA@DOX/SPIO, respectively. Each group underwent IVIM-DWI scanning at baseline and 3, 14, 21, and 28 days after treatment. D* and f values were derived using GE AW 4.5 post-processing station. All mice were sacrificed for pathological examination.

RESULTS

The D* value of all three groups showed an upward trend, with the highest increase in G1 and the lowest in G3. Conversely, the f value of all groups trended to decrease within 7 days, of which G3 showed the most significant decline. Immunohistochemical staining revealed that vascular endothelial growth factor (VEGF)-positive staining rate and the microvessel density (MVD) of the tumors in G3 were significantly lower than those of the other groups (P < 0.05). The D* and f values were significantly and positively correlated to CD31 (r = 0.654, P < 0.001; r = 0.712, P < 0.001) and VEGF (r = 0.694, P < 0.001; r = 0.664, P < 0.001).

CONCLUSION

IVIM-DWI-derived parameters D* and f are valuable indicators for the evaluation of the antitumor microcirculation changes of multifunctional nanosystem.

Dynamics of driftless preconcentration using ion concentration polarization leveraged by convection and diffusion

Over the past several decades, separation and preconcentration methods of (bio)molecules have been actively developed for various biomedical and chemical processes such as disease diagnostics, point of care test and environmental monitoring. Among the great developments of the electrokinetic method in a micro/nanofluidic platform is the ion concentration polarization (ICP) phenomenon, in which a target molecule is accumulated near a permselective nanoporous membrane under an applied electric field. ICP method has been actively studied due to its easy implementation and high preconcentration/separation efficiency. However, the dynamic behavior of preconcentrated analytes has not yet been fully studied, especially driftless migration, where the applied electric field is orthogonal to the direction of the drift migration. Here, we demonstrate anomalous shapes of preconcentrated analytes (either plug or dumbbell shape) and the morphologies were analytically modeled by the leverage of convection and diffusion migration. This model was experimentally verified with various lengths of DNA and the limiting cases (convection-free environment in paper-based microfluidic device and extremely low diffusivity of red blood cells) were also shown to confirm the model. Thus, this study not only provides an insight into the fundamental electrokinetic dynamics of molecules in an ICP platform but also plays a guiding role for the design of a nanofluidic preconcentrator for a lab on a chip application.

Evidence of the diffusion time dependence of intravoxel incoherent motion in the brain

PURPOSE

To investigate the diffusion time (T_D ) dependence of intravoxel incoherent motion (IVIM) signals in the brain.

METHODS

A 3-compartment IVIM model was proposed to characterize 2 types of microcirculatory flows in addition to tissue water in the brain: flows that cross multiple vascular segments (pseudo-diffusive) and flows that stay in 1 segment (ballistic) within T_D . The model was first evaluated using simulated flow signals. Experimentally, flow-compensated (FC) pulsed-gradient spin-echo (PGSE) and oscillating-gradient spin-echo (OGSE) sequences were tested using a flow phantom and then used to examine IVIM signals in the mouse brain with T_D ranging from ~2.5 ms to 40 ms on an 11.7T scanner.

RESULTS

By fitting the model to simulated flow signals, we demonstrated the T_D dependency of the estimated fraction of pseudo-diffusive flow and the pseudo-diffusion coefficient (D*), which were dictated by the characteristic timescale of microcirculatory flow (τ). Flow phantom experiments validated that the OGSE and FC-PGSE sequences were not susceptible to the change in flow velocity. In vivo mouse brain data showed that both the estimated fraction of pseudo-diffusive flow and D* increased significantly as T_D increased.

CONCLUSION

We demonstrated that IVIM signals measured in the brain are T_D -dependent, potentially because more microcirculatory flows approach the pseudo-diffusive limit as T_D increases with respect to τ. Measuring the T_D dependency of IVIM signals may provide additional information on microvascular flows in the brain.

Motion-Induced Signal Loss in In Vivo Cardiac Diffusion-Weighted Imaging

LEVEL OF EVIDENCE

5 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2020;51:319-320.