Model-Based High-Definition Dynamic Contrast Enhanced MRI for Concurrent Estimation of Perfusion and Microvascular Permeability

Implementation and evaluation of a dynamic contrast-enhanced MR perfusion protocol for glioblastoma using a 0.35 T MRI-Linac system

PURPOSE

MRI-linear accelerator (MRI-Linac) systems allow for daily tracking of MRI changes during radiotherapy (RT). Since one common MRI-Linac operates at 0.35 T, there are efforts towards developing protocols at that field strength. In this study we demonstrate the implementation of a post-contrast 3DT1-weighted (3D-T1w) and dynamic contrast-enhancement (DCE) protocol to assess glioblastoma response to RT using a 0.35 T MRI-Linac.

METHODS AND MATERIALS

The protocol implemented was used to acquire 3D-T1w and DCE data from a flow phantom and two patients with glioblastoma (a responder and a non-responder) who underwent RT on a 0.35 T MRI-Linac. The detection of post-contrast-enhanced volumes was evaluated by comparing the 3DT1w images from the 0.35 T MRI-Linac to images obtained using a 3 T scanner. The DCE data were tested temporally and spatially using data from a flow phantom and patients. Ktrans maps were derived from DCE at three time points (a week before treatment-Pre RT, four weeks through treatment-Mid RT, and three weeks after treatment-Post RT) and were validated with patients' treatment outcomes.

RESULTS

The 3D-T1w contrast-enhancement volumes were visually and volumetrically similar between 0.35 T MRI-Linac and 3 T. DCE images showed temporal stability, and associated Ktrans maps were consistent with patient response to treatment. On average, Ktrans values showed a 54 % decrease and 8.6 % increase for a responder and non-responder respectively when Pre RT and Mid RT images were compared.

CONCLUSION

Our findings support the feasibility of obtaining post-contrast 3D-T1w and DCE data from patients with glioblastoma using a 0.35 T MRI-Linac system.

Model-Based Needle Identification Using Image Analysis and Needle Library Matching for Ultrasound-Guided Kidney Biopsy: A Feasibility Study

OBJECTIVE

The aim of the work described here was to determine the feasibility of using a novel biopsy needle detection technique that achieves high sensitivity and specificity in a trade-off of resolution, detectability and depth of imaging.

METHODS

The proposed needle detection method consists of a model-based image analysis, temporal needle projection and needle library matching: (i) Image analysis was formulated under the signal decomposition framework; (ii) temporal projection converted the time-resolved needle dynamics into a single image of the desired needle; and (iii) the enhanced needle structure was spatially refined by matching a long, straight linear object in the needle library. The efficacy was examined with respect to different needle visibility.

RESULTS

Our method effectively eliminated confounding effects of the background tissue artifacts more robustly than conventional methods, thus improving needle visibility even with the low contrast between the needle and tissue. The improvement in needle structure further resulted in an improvement in estimation performance for the trajectory angle and tip position.

CONCLUSION

Our three-step needle detection method can reliably detect needle position without the need for external devices, increasing the needle conspicuity and reducing motion sensitivity.

Joint Optimization of k-t Sampling Pattern and Reconstruction of DCE MRI for Pharmacokinetic Parameter Estimation

This work proposes to develop and evaluate a deep learning framework that jointly optimizes k-t sampling patterns and reconstruction for head and neck dynamic contrast-enhanced (DCE) MRI aiming to reduce bias and uncertainty of pharmacokinetic (PK) parameter estimation. 2D Cartesian phase encoding k-space subsampling patterns for a 3D spoiled gradient recalled echo (SPGR) sequence along a time course of DCE MRI were jointly optimized in a deep learning-based dynamic MRI reconstruction network by a loss function concerning both reconstruction image quality and PK parameter estimation accuracy. During training, temporal k-space data sharing scheme was optimized as well. The proposed method was trained and tested by multi-coil complex digital reference objects of DCE images (mcDROs). The PK parameters estimated by the proposed method were compared with two published iterative DCE MRI reconstruction schemes using normalized root mean squared errors (NRMSEs) and Bland-Altman analysis at temporal resolutions of [Formula: see text] = 2s, 3s, 4s, and 5s, which correspond to undersampling rates of R = 50, 34, 25, and 20. The proposed method achieved low PK parameter NRMSEs at all four temporal resolutions compared with the benchmark methods on testing mcDROs. The Bland-Altman plots demonstrated that the proposed method reduced PK parameter estimation bias and uncertainty in tumor regions at temporal resolution of 2s. The proposed method also showed robustness to contrast arrival timing variations across patients. This work provides a potential way to increase PK parameter estimation accuracy and precision, and thus facilitate the clinical translation of DCE MRI.

Joint Reconstruction of Vascular Structure and Function Maps in Dynamic Contrast Enhanced MRI Using Vascular Heterogeneity Priors

This work introduces a novel, joint reconstruction of vascular structure and microvascular function maps directly from highly undersampled data in k - t space using vascular heterogeneity priors for high-definition, dynamic contrast-enhanced (DCE) MRI. In DCE MRI, arteries and veins are characterized by rapid, high uptake and wash-out of contrast agents (CA). On the other hand, depending on CA uptake and wash-out signal patterns, capillary tissues can be categorized into highly perfused, moderately perfused, and necrotic regions. Given the above considerations, macrovascular maps are generated as a prior to differentiate penalties on arteries relative to capillary tissues during image reconstruction. Furthermore, as a microvascular prior, contrast dynamics in capillary regions are represented in a low dimensional space using a finite number of basic vectors that reflect actual tissue-specific signal patterns. Both vascular structure and microvascular function maps are jointly estimated by solving a constrained optimization problem in which the above vascular heterogeneity priors are represented by spatially weighted nonnegative matrix factorization. Retrospective and prospective experiments are performed to validate the effectiveness of the proposed method in generating well-defined vascular structure and microvascular function maps for patients with brain tumor at high reduction factors.

A deep matrix factorization framework for identifying underlying tissue-specific patterns of DCE-MRI: applications for molecular subtype classification in breast cancer

Objective.Breast cancer is heterogeneous in that different angiogenesis and blood flow characteristics could be present within a tumor. The pixel kinetics of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) can assume several distinct signal patterns related to specific tissue characteristics. Identification of the latent, tissue-specific dynamic patterns of intratumor heterogeneity can shed light on the biological mechanisms underlying the heterogeneity of tumors.Approach.To mine this information, we propose a deep matrix factorization-based dynamic decomposition (DMFDE) model specifically designed according to DCE-MRI characteristics. The time-series imaging data were decomposed into tissue-specific dynamic patterns and their corresponding proportion maps. The image pixel matrix and the reference matrix of population-level kinetics obtained by clustering the dynamic signals were used as the inputs. Two multilayer neural network branches were designed to collaboratively project the input matrix into a latent dynamic pattern and a dynamic proportion matrix, which was justified using simulated data. Clinical implications of DMFDE were assessed by radiomics analysis of proportion maps obtained from the tumor/parenchyma region for classifying the luminal A subtype.Main results.The decomposition performance of DMFDE was evaluated by the root mean square error and was shown to be better than that of the conventional convex analysis of mixtures (CAM) method. The predictive model withK = 3, 4, and 5 dynamic proportion maps generated AUC values of 0.780, 0.786 and 0.790, respectively, in distinguishing between luminal A and nonluminal A tumors, which are better than the CAM method (AUC = 0.726). The combination of statistical features from images with different proportion maps has the highest prediction value (AUC = 0.813), which is significantly higher than that based on CAM.Conclusion.This proposed method identified the latent dynamic patterns associated with different molecular subtypes, and radiomics analysis based on the pixel compositions of the uncovered dynamic patterns was able to determine molecular subtypes of breast cancer.

Assessing the reproducibility of high temporal and spatial resolution dynamic contrast-enhanced magnetic resonance imaging in patients with gliomas

Temporal and spatial resolution of dynamic contrast-enhanced MR imaging (DCE-MRI) is critical to reproducibility, and the reproducibility of high-resolution (HR) DCE-MRI was evaluated. Thirty consecutive patients suspected to have brain tumors were prospectively enrolled with written informed consent. All patients underwent both HR-DCE (voxel size, 1.1 × 1.1 × 1.1 mm³; scan interval, 1.6 s) and conventional DCE (C-DCE; voxel size, 1.25 × 1.25 × 3.0 mm³; scan interval, 4.0 s) MRI. Regions of interests (ROIs) for enhancing lesions were segmented twice in each patient with glioblastoma (n = 7) to calculate DCE parameters (K^trans, V_p, and V_e). Intraclass correlation coefficients (ICCs) of DCE parameters were obtained. In patients with gliomas (n = 25), arterial input functions (AIFs) and DCE parameters derived from T2 hyperintense lesions were obtained, and DCE parameters were compared according to WHO grades. ICCs of HR-DCE parameters were good to excellent (0.84–0.95), and ICCs of C-DCE parameters were moderate to excellent (0.66–0.96). Maximal signal intensity and wash-in slope of AIFs from HR-DCE MRI were significantly greater than those from C-DCE MRI (31.85 vs. 7.09 and 2.14 vs. 0.63; p < 0.001). Both 95^th percentile K^trans and V_e from HR-DCE and C-DCE MRI could differentiate grade 4 from grade 2 and 3 gliomas (p < 0.05). In conclusion, HR-DCE parameters generally showed better reproducibility than C-DCE parameters, and HR-DCE MRI provided better quality of AIFs.