Quantitative transport mapping (QTM) of the kidney with an approximate microvascular network

Quantitative non-contrast perfusion MRI in the body using arterial spin labeling

Arterial spin labeling (ASL) is a non-invasive magnetic resonance imaging (MRI) method that enables the assessment and the quantification of perfusion without the need for an exogenous contrast agent. ASL was originally developed in the early 1990s to measure cerebral blood flow. The utility of ASL has since then broadened to encompass various organ systems, offering insights into physiological and pathological states. In this review article, we present a synopsis of ASL for quantitative non-contrast perfusion MRI, as a contribution to the special issue titled "Quantitative MRI-how to make it work in the body?" The article begins with an introduction to ASL principles, followed by different labeling strategies, such as pulsed, continuous, pseudo-continuous, and velocity-selective approaches, and their role in perfusion quantification. We proceed to address the technical challenges associated with ASL in the body and outline some of the innovative approaches devised to surmount these issues. Subsequently, we summarize potential clinical applications, challenges, and state-of-the-art ASL methods to quantify perfusion in some of the highly perfused organs in the thorax (lungs), abdomen (kidneys, liver, pancreas), and pelvis (placenta) of the human body. The article concludes by discussing future directions for successful translation of quantitative ASL in body imaging.

Quantitative transport mapping of multi-delay arterial spin labeling MRI detects early blood perfusion alterations in Alzheimer's disease

BACKGROUND

Quantitative transport mapping (QTM) of blood velocity, based on the transport equation has been demonstrated higher accuracy and sensitivity of perfusion quantification than the traditional Kety's method-based cerebral blood flow (CBF). This study aimed to investigate the associations between QTM velocity and cognitive function in Alzheimer's disease (AD) using multiple post-labeling delay arterial spin labeling (ASL) MRI.

METHODS

A total of 128 subjects (21 normal controls (NC), 80 patients with mild cognitive impairment (MCI), and 27 AD) were recruited prospectively. All participants underwent MRI examination and neuropsychological evaluation. QTM velocity and traditional CBF maps were computed from multiple delay ASL. Regional quantitative perfusion measurements were performed and compared to study group differences. We tested the hypothesis that cognition declines with reduced cerebral blood perfusion with consideration of age and gender effects.

RESULTS

In cortical gray matter (GM) and the hippocampus, QTM velocity and CBF showed decreased values in the AD group compared to NC and MCI groups; QTM velocity, but not CBF, showed a significant difference between MCI and NC groups. QTM velocity and CBF showed values decreasing with age; QTM velocity, but not CBF, showed a significant gender difference between male and female. QTM velocity and CBF in the hippocampus were positively correlated with cognition, including global cognition, memory, executive function, and language function.

CONCLUSION

This study demonstrated an increased sensitivity of QTM velocity as compared with the traditional Kety's method-based CBF. Specifically, we observed only in QTM velocity, reduced perfusion velocity in GM and the hippocampus in MCI compared with NC. Both QTM velocity and CBF demonstrated a reduction in AD vs. controls. Decreased QTM velocity and CBF in the hippocampus were correlated with poor cognitive measures. These findings suggest QTM velocity as potential biomarker for early AD blood perfusion alterations and it could provide an avenue for early intervention of AD.

Identifiability of spatiotemporal tissue perfusion models

Objective.Standard models for perfusion quantification in DCE-MRI produce a bias by treating voxels as isolated systems. Spatiotemporal models can remove this bias, but it is unknown whether they are fundamentally identifiable. The aim of this study is to investigate this question in silico using one-dimensional toy systems with a one-compartment blood flow model and a two-compartment perfusion model.Approach.For each of the two models, identifiability is explored theoretically and in-silico for three systems. Concentrations over space and time are simulated by forward propagation. Different levels of noise and temporal undersampling are added to investigate sensitivity to measurement error. Model parameters are fitted using a standard gradient descent algorithm, applied iteratively with a stepwise increasing time window. Model fitting is repeated with different initial values to probe uniqueness of the solution. Reconstruction accuracy is quantified for each parameter by comparison to the ground truth.Main results.Theoretical analysis shows that flows and volume fractions are only identifiable up to a constant, and that this degeneracy can be removed by proper choice of parameters. Simulations show that in all cases, the tissue concentrations can be reconstructed accurately. The one-compartment model shows accurate reconstruction of blood velocities and arterial input functions, independent of the initial values and robust to measurement error. The two-compartmental perfusion model was not fully identifiable, showing good reconstruction of arterial velocities and input functions, but multiple valid solutions for the perfusion parameters and venous velocities, and a strong sensitivity to measurement error in these parameters.Significance.These results support the use of one-compartment spatiotemporal flow models, but two-compartment perfusion models were not sufficiently identifiable. Future studies should investigate whether this degeneracy is resolved in more realistic 2D and 3D systems, by adding physically justified constraints, or by optimizing experimental parameters such as injection duration or temporal resolution.

Quantitative transport mapping of multi-delay arterial spin labeling MRI detects early blood perfusion alteration in Alzheimer's disease

Background

Quantitative transport mapping (QTM) of blood velocity, based on the transport equation has been demonstrated higher accuracy and sensitivity of perfusion quantification than the traditional Kety's method-based blood flow (Kety flow). This study aimed to investigate the associations between QTM velocity and cognitive function in Alzheimer's disease (AD) using multiple post-labeling delay arterial spin labeling (ASL) MRI.

Methods

A total of 128 subjects (21 normal controls (NC), 80 patients with mild cognitive impairment (MCI), and 27 AD) were recruited prospectively. All participants underwent MRI examination and neuropsychological evaluation. QTM velocity and traditional Kety flow maps were computed from multiple delay ASL. Regional quantitative perfusion measurements were performed and compared to study group differences. We tested the hypothesis that cognition declines with reduced cerebral blood flow with consideration of age and gender effects.

Results

In cortical gray matter (GM) and the hippocampus, QTM velocity and Kety flow showed decreased values in AD group compared to NC and MCI groups; QTM velocity, but not Kety flow, showed a significant difference between MCI and NC groups. QTM velocity and Kety flow showed values decreasing with age; QTM velocity, but not Kety flow, showed a significant gender difference between male and female. QTM velocity and Kety flow in the hippocampus were positively correlated with cognition, including global cognition, memory, executive function, and language function.

Conclusion

This study demonstrated an increased sensitivity of QTM velocity as compared with the traditional Kety flow. Specifically, we observed only in QTM velocity, reduced perfusion velocity in GM and the hippocampus in MCI compared with NC. Both QTM velocity and Kety flow demonstrated reduction in AD vs controls. Decreased QTM velocity and Kety flow in the hippocampus were correlated with cognitive measures. These findings suggest QTM velocity as an improved biomarker for early AD blood flow alterations.

Current status in spatiotemporal analysis of contrast-based perfusion MRI

In perfusion MRI, image voxels form a spatially organized network of systems, all exchanging indicator with their immediate neighbors. Yet the current paradigm for perfusion MRI analysis treats all voxels or regions-of-interest as isolated systems supplied by a single global source. This simplification not only leads to long-recognized systematic errors but also fails to leverage the embedded spatial structure within the data. Since the early 2000s, a variety of models and implementations have been proposed to analyze systems with between-voxel interactions. In general, this leads to large and connected numerical inverse problems that are intractible with conventional computational methods. With recent advances in machine learning, however, these approaches are becoming practically feasible, opening up the way for a paradigm shift in the approach to perfusion MRI. This paper seeks to review the work in spatiotemporal modelling of perfusion MRI using a coherent, harmonized nomenclature and notation, with clear physical definitions and assumptions. The aim is to introduce clarity in the state-of-the-art of this promising new approach to perfusion MRI, and help to identify gaps of knowledge and priorities for future research.

MR Perfusion Imaging for Kidney Disease

Renal perfusion reflects overall function of a kidney. As an important indicator of kidney diseases, renal perfusion can be noninvasively measured by multiple methods of MR imaging, such as dynamic contrast-enhanced MR imaging, intravoxel incoherent motion analysis, and arterial spin labeling method. In this article we introduce the principle of the methods, review their recent technical improvements, and then focus on summarizing recent applications of the methods in assessing various renal diseases. By this review, we demonstrate the capability and clinical potential of the imaging methods, with the hope of accelerating their adoption to clinical practice.

Advances in Image Processing for Epileptogenic Zone Detection with MRI

Focal epilepsy is a common and severe neurologic disorder. Neuroimaging aims to identify the epileptogenic zone (EZ), preferably as a macroscopic structural lesion. For approximately a third of patients with chronic drug-resistant focal epilepsy, the EZ cannot be precisely identified using standard 3.0-T MRI. This may be due to either the EZ being undetectable at imaging or the seizure activity being caused by a physiologic abnormality rather than a structural lesion. Computational image processing has recently been shown to aid radiologic assessments and increase the success rate of uncovering suspicious regions by enhancing their visual conspicuity. While structural image analysis is at the forefront of EZ detection, physiologic image analysis has also been shown to provide valuable information about EZ location. This narrative review summarizes and explains the current state-of-the-art computational approaches for image analysis and presents their potential for EZ detection. Current limitations of the methods and possible future directions to augment EZ detection are discussed.

Cerebral oxygen extraction fraction declines with ventricular enlargement in patients with normal pressure hydrocephalus

OBJECTIVE

Normal pressure hydrocephalus (NPH) is a neurodegenerative disease that is potentially reversible by shunt surgery in approximately 60% of patients. Imaging may provide a means to investigate brain tissue viability and oxygen metabolism in NPH patients.

METHODS

Oxygen extraction fraction (OEF) mapping was generated from 3D multi-echo gradient echo MRI (mGRE) data using QQ-CCTV algorithm and cerebral blood flow (CBF) using 3D arterial spin labeling (ASL) MRI data, thereby calculating the cerebral metabolic rate of oxygen (CMRO2 = CBF × OEF × [H]a) in 16 NPH patients. Regression analyses using cortical gray matter and deep gray matter regions were conducted with age, gender, CSF stroke volume and normalized ventricular volume as independent variables.

RESULTS

OEF showed significant negative correlations with normalized brain ventricular volumes in the whole brain (p = 0.004, q = 0.01), cortical gray matter (p = 0.004, q = 0.01), caudate (p = 0.02, q = 0.04), and pallidum (p = 0.03, q = 0.04), but no significant correlation with CSF stroke volume (q > 0.05). There was no significant finding with CBF or CMRO2.

CONCLUSION

In NPH patients, low OEF in several regions was significantly correlated with large ventricular volumes, indicating decreased tissue oxygen metabolism with increased NPH severity. OEF mapping may provide a functional understanding of neurodegeneration in NPH and may improve monitoring of disease course and treatment outcomes.

Fluid Mechanics Approach to Perfusion Quantification: Vasculature Computational Fluid Dynamics Simulation, Quantitative Transport Mapping (QTM) Analysis of Dynamics Contrast Enhanced MRI, and Application in Nonalcoholic Fatty Liver Disease Classification

OBJECTIVE

We quantify liver perfusion using quantitative transport mapping (QTM) method that is free of arterial input function (AIF). QTM method is validated in a vasculature computational fluid dynamics (CFD) simulation and is applied for processing dynamic contrast enhanced (DCE) MRI images in differentiating liver with nonalcoholic fatty liver disease (NAFLD) from healthy controls using pathology reference in a preclinical rabbit model.

METHODS

QTM method was validated on a liver perfusion simulation based on fluid dynamics using a rat liver vasculature model and the mass transport equation. In the NAFLD grading task, DCE MRI images of 7 adult rabbits with methionine choline-deficient diet-induced nonalcoholic steatohepatitis (NASH), 8 adult rabbits with simple steatosis (SS) were acquired and processed using QTM method and dual-input two compartment Kety's method respectively. Statistical analysis was performed on six perfusion parameters: velocity magnitude | u | derived from QTM, liver arterial blood flow LBFa, liver venous blood flow LBFv, permeability Ktrans, blood volume Vp and extravascular space volume Ve averaged in liver ROI.

RESULTS

In the simulation, QTM method successfully reconstructed blood flow, reduced error by 48% compared to Kety's method. In the preclinical study, only QTM |u| showed significant difference between high grade NAFLD group and low grade NAFLD group.

CONCLUSION

QTM postprocesses DCE-MRI automatically through deconvolution in space and time to solve the inverse problem of the transport equation. Comparing with Kety's method, QTM method showed higher accuracy and better differentiation in NAFLD classification task.

SIGNIFICANCE

We propose to apply QTM method in liver DCE MRI perfusion quantification.

Prediction of Lung Shunt Fraction for Yttrium-90 Treatment of Hepatic Tumors Using Dynamic Contrast Enhanced MRI with Quantitative Perfusion Processing

There is no noninvasive method to estimate lung shunting fraction (LSF) in patients with liver tumors undergoing Yttrium-90 (Y90) therapy. We propose to predict LSF from noninvasive dynamic contrast enhanced (DCE) MRI using perfusion quantification. Two perfusion quantification methods were used to process DCE MRI in 25 liver tumor patients: Kety's tracer kinetic modeling with a delay-fitted global arterial input function (AIF) and quantitative transport mapping (QTM) based on the inversion of transport equation using spatial deconvolution without AIF. LSF was measured on SPECT following Tc-99m macroaggregated albumin (MAA) administration via hepatic arterial catheter. The patient cohort was partitioned into a low-risk group (LSF ≤ 10%) and a high-risk group (LSF > 10%). Results: In this patient cohort, LSF was positively correlated with QTM velocity |u| (r = 0.61, F = 14.0363, p = 0.0021), and no significant correlation was observed with Kety's parameters, tumor volume, patient age and gender. Between the low LSF and high LSF groups, there was a significant difference for QTM |u| (0.0760 ± 0.0440 vs. 0.1822 ± 0.1225 mm/s, p = 0.0011), and Kety's Ktrans (0.0401 ± 0.0360 vs 0.1198 ± 0.3048, p = 0.0471) and Ve (0.0900 ± 0.0307 vs. 0.1495 ± 0.0485, p = 0.0114). The area under the curve (AUC) for distinguishing between low LSF and high LSF was 0.87 for |u|, 0.80 for Ve and 0.74 for Ktrans. Noninvasive prediction of LSF is feasible from DCE MRI with QTM velocity postprocessing.

Quantitative transport mapping (QTM) for differentiating benign and malignant breast lesion: Comparison with traditional kinetics modeling and semi-quantitative enhancement curve characteristics

PURPOSE

To test the feasibility of using quantitative transport mapping (QTM) method, which is based on the inversion of transport equation using spatial deconvolution without any arterial input function, for automatically postprocessing dynamic contrast enhanced MRI (DCE-MRI) to differentiate malignant and benign breast tumors.

MATERIALS AND METHODS

Breast DCE-MRI data with biopsy confirmed malignant (n = 13) and benign tumors (n = 13) was used to assess QTM velocity (|u|) and diffusion coefficient (D), volume transfer constant (K^trans), volume fraction of extravascular extracellular space (V_e) from kinetics method, and traditional enhancement curve characteristics (ECC: amplitude A, wash-in rate α, wash-out rate β). A Mann-Whitney U test and receiver operating characteristic curve (ROC) analysis were performed to assess the diagnostic performance of these parameters for distinguishing between benign and malignant tumors.

RESULTS

Between malignant and benign tumors, there was a significant difference in |u| and K^trans, (p = 0.0066, 0.0274, respectively), but not in D, V_e, A, α and β (p = 0.1119, 0.2382, 0.4418,0.2592 and 0.9591, respectively). ROC area-under-the-curve was 0.82, 0.75 (95% confidence level 0.60-0.95, 0.51-0.90) for |u| and K^trans, respectively.

CONCLUSION

QTM postprocesses DCE-MRI automatically through deconvolution in space and time to solve the inverse problem of the transport equation. Comparing with traditional kinetics method and ECC, QTM method showed better diagnostic accuracy in differentiating benign from malignant breast tumors in this study.

DCE-MRI quantitative transport mapping for noninvasively detecting hypoxia inducible factor-1α, epidermal growth factor receptor overexpression, and Ki-67 in nasopharyngeal carcinoma patients

BACKGROUND

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) has the potential to noninvasively detect expression of hypoxia inducible factor-1-alpha (HIF-1α), epidermal growth factor receptor (EGFR), and Ki-67 in nasopharyngeal carcinoma (NPC) by quantitatively measuring tumor blood flow, vascularity, and permeability.

PURPOSE

We aim to explore the utility of DCE-MRI in detecting HIF-1α, EGFR, and Ki-67 expression levels using traditional Kety's/Tofts' modeling and quantitative transport mapping (QTM).

MATERIALS AND METHODS

Eighty-nine NPC patients underwent DCE-MRI before treatment were enrolled. DCE-MRI was processed to generate the following kinetic parameters: |u| and D from the QTM model, tumor blood flow (TBF) from Kety's model, and K^trans, V_e, and K_ep from Tofts' model. Pretreatment levels of HIF-1α, EGFR, and Ki-67 were assessed by immunohistochemistry and classified into low and high expression groups.

RESULTS

|u| (p < 0.001) and TBF (p = 0.015) values were significantly higher in the HIF-1α high-expression group compared to low-expression group. Only K^trans (p = 0.016) was significantly higher in the EGFR high-expression group. Only |u| (p < 0.001) values were significantly higher in the Ki-67 high-expression group compared to low-expression group. Multiple linear regression analyses showed that |u| independently correlated with HIF-1α and Ki-67 expression, and K^trans independently correlated with EGFR. The areas under the ROC curves of |u| for HIF-1α and Ki-67, and K^trans for EGFR were 0.83, 0.74, and 0.70, respectively.

CONCLUSION

|u| and K^trans derived from DCE-MRI may be considered as noninvasive imaging markers for detecting hypoxia and proliferation in NPC patients.

Quantitative transport mapping (QTM) of the kidney with an approximate microvascular network

Purpose

Proof‐of‐concept study of mapping renal blood flow vector field according to the inverse solution to a mass transport model of time resolved tracer‐labeled MRI data.

Theory and Methods

To determine tissue perfusion according to the underlying physics of spatiotemporal tracer concentration variation, the mass transport equation is integrated over a voxel with an approximate microvascular network for fitting time‐resolved tracer imaging data. The inverse solution to the voxelized transport equation provides the blood flow vector field, which is referred to as quantitative transport mapping (QTM). A numerical microvascular network modeling the kidney with computational fluid dynamics reference was used to verify the accuracy of QTM and the current Kety’s method that uses a global arterial input function. Multiple post‐label delay arterial spin labeling (ASL) of the kidney on seven subjects was used to assess QTM in vivo feasibility.

Results

Against the ground truth in the numerical model, the error in flow estimated by QTM (18.6%) was smaller than that in Kety’s method (45.7%, 2.5‐fold reduction). The in vivo kidney perfusion quantification by QTM (cortex: 443 ± 58 mL/100 g/min and medulla: 190 ± 90 mL/100 g/min) was in the range of that by Kety’s method (482 ± 51 mL/100 g/min in the cortex and 242 ± 73 mL/100 g/min in the medulla), and QTM provided better flow homogeneity in the cortex region.

Conclusions

QTM flow velocity mapping is feasible from multi‐delay ASL MRI data based on inverting the transport equation. In a numerical simulation, QTM with deconvolution in space and time provided more accurate perfusion quantification than Kety’s method with deconvolution in time only.