Practical considerations in T1 mapping of prostate for dynamic contrast enhancement pharmacokinetic analyses

Radiology of fibrosis part III: genitourinary system

Fibrosis is a pathological process involving the abnormal deposition of connective tissue, resulting from improper tissue repair in response to sustained injury caused by hypoxia, infection, or physical damage. It can impact any organ, leading to their dysfunction and eventual failure. Additionally, tissue fibrosis plays an important role in carcinogenesis and the progression of cancer.Early and accurate diagnosis of organ fibrosis, coupled with regular surveillance, is essential for timely disease-modifying interventions, ultimately reducing mortality and enhancing quality of life. While extensive research has already been carried out on the topics of aberrant wound healing and fibrogenesis, we lack a thorough understanding of how their relationship reveals itself through modern imaging techniques.This paper focuses on fibrosis of the genito-urinary system, detailing relevant imaging technologies used for its detection and exploring future directions.

Comparison of DCE-MRI parametric mapping using MP2RAGE and variable flip angle T1 mapping

Quantitative dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) measures the rate of transfer of contrast agent from the vascular space to the tissue space by fitting signal-time data to pharmacokinetic models. However, these models are very sensitive to errors in T₁ mapping. Accurate T₁ mapping is necessary for high quality quantitative DCE-MRI studies. This study compares magnetization prepared rapid (two) gradient echo sequence (MP2RAGE) T₁-mapping accuracy to the conventional variable flip angle (VFA) approach, and also determines the effect of the new T₁-mapping method on the K^trans parameter. VFA and MP2RAGE T₁ values were compared to the gold standard inverse recovery (IR) method in phantom over manually drawn ROIs. In vivo, ROIs were manually drawn over prostate and prostatic lesions. Average T₁ values over ROIs were compared and K^trans maps for each method were calculated via the extended Tofts model. VFA-T₁ maps overestimated T₁ values by up to 50% compared to gold standard IR T₁ values in phantom. MP2RAGE differed by up to 9%. MP2RAGE-T₁ and K^trans values were significantly different from VFA values over prostatic lesions (p < 0.05). K^trans was consistently underestimated using VFA compared to MP2RAGE (p < 0.05). MP2RAGE T₁ maps are shown to be more accurate, leading to more reliable pharmacokinetic modeling. This can potentially lead to better lesion characterization and improve clinical outcomes.

Effect of magnetic resonance imaging pre-processing on the performance of model-based prostate tumor probability mapping

Objective. Multi-parametric magnetic resonance imaging (mpMRI) has become an important tool for the detection of prostate cancer in the past two decades. Despite the high sensitivity of MRI for tissue characterization, it often suffers from a lack of specificity. Several well-established pre-processing tools are publicly available for improving image quality and removing both intra- and inter-patient variability in order to increase the diagnostic accuracy of MRI. To date, most of these pre-processing tools have largely been assessed individually. In this study we present a systematic evaluation of a multi-step mpMRI pre-processing pipeline to automate tumor localization within the prostate using a previously trained model.Approach. The study was conducted on 31 treatment-naïve prostate cancer patients with a PI-RADS-v2 compliant mpMRI examination. Multiple methods were compared for each pre-processing step: (1) bias field correction, (2) normalization, and (3) deformable multi-modal registration. Optimal parameter values were estimated for each step on the basis of relevant individual metrics. Tumor localization was then carried out via a model-based approach that takes both mpMRI and prior clinical knowledge features as input. A sequential optimization approach was adopted for determining the optimal parameters and techniques in each step of the pipeline.Main results. The application of bias field correction alone increased the accuracy of tumor localization (area under the curve (AUC) = 0.77;p-value = 0.004) over unprocessed data (AUC = 0.74). Adding normalization to the pre-processing pipeline further improved diagnostic accuracy of the model to an AUC of 0.85 (p-value = 0.000 12). Multi-modal registration of apparent diffusion coefficient images to T2-weighted images improved the alignment of tumor locations in all but one patient, resulting in a slight decrease in accuracy (AUC = 0.84;p-value = 0.30).Significance. Overall, our findings suggest that the combined effect of multiple pre-processing steps with optimal values has the ability to improve the quantitative classification of prostate cancer using mpMRI. Clinical trials: NCT03378856 and NCT03367702.

Evaluation of the differentiation of benign and malignant breast lesions using synthetic relaxometry and the Kaiser score

Objective

To investigate whether there is added value of quantitative parameters from synthetic magnetic resonance imaging (SyMRI) as a complement to the Kaiser score (KS) to differentiate benign and malignant breast lesions.

Materials and methods

In this single-institution study, 122 patients who underwent breast MRI from March 2020 to May 2021 were retrospectively analyzed. SyMRI and dynamic contrast-enhanced MRI were performed using a 3.0-T system. Two experienced radiologists independently assigned the KS and measured the quantitative values of T1 relaxation time (T1), T2 relaxation time (T2), and proton density (PD) from SyMRI. Pathology was regarded as the gold standard. The diagnostic values were compared using the appropriate statistical tests.

Results

There were 122 lesions (86 malignant and 36 benign) in 122 women. The T1 value was identified as the only independent factor for the differentiation of malignant and benign lesions. The diagnostic accuracy of incorporating the T1 into the KS protocol (T1+KS) was 95.1% and 92.1% for all lesions (ALL) and The American College of Radiology (ACR) Breast Imaging Reporting and Data System (BI-RADS) category 4 lesions, respectively, which was significantly higher than that of either T1 (ALL: 82.8%, P = 0.0001; BI-RADS 4: 78.9%, P = 0.002) or KS (ALL: 90.2%, P = 0.031; BI-RADS 4: 84.2%, P = 0.031) alone. The sensitivity and specificity of T1+KS were also higher than those of the T1 or KS alone. The combined diagnosis could have avoided another 15.6% biopsies compared with using KS alone.

Conclusions

Incorporating T1 into the KS protocol improved both the sensitivity and specificity to differentiate benign and malignant breast lesions, thus avoiding unnecessary invasive procedures.

Flow blockage disrupts cilia-driven fluid transport in the epileptic brain

A carpet of ependymal motile cilia lines the brain ventricular system, forming a network of flow channels and barriers that pattern cerebrospinal fluid (CSF) flow at the surface. This CSF transport system is evolutionary conserved, but its physiological function remains unknown. Here we investigated its potential role in epilepsy with studies focused on CDKL5 deficiency disorder (CDD), a neurodevelopmental disorder with early-onset epilepsy refractory to seizure medications and the most common cause of infant epilepsy. CDKL5 is a highly conserved X-linked gene suggesting its function in regulating cilia length and motion in the green alga Chlamydomonas might have implication in the etiology of CDD. Examination of the structure and function of airway motile cilia revealed both the CDD patients and the Cdkl5 knockout mice exhibit cilia lengthening and abnormal cilia motion. Similar defects were observed for brain ventricular cilia in the Cdkl5 knockout mice. Mapping ependymal cilia generated flow in the ventral third ventricle (v3V), a brain region with important physiological functions showed altered patterning of flow. Tracing of cilia-mediated inflow into v3V with fluorescent dye revealed the appearance of a flow barrier at the inlet of v3V in Cdkl5 knockout mice. Analysis of mice with a mutation in another epilepsy-associated kinase, Yes1, showed the same disturbance of cilia motion and flow patterning. The flow barrier was also observed in the Foxj1^± and FOXJ1Cre^ERT:Cdkl5^y/fl mice, confirming the contribution of ventricular cilia to the flow disturbances. Importantly, mice exhibiting altered cilia-driven flow also showed increased susceptibility to anesthesia-induced seizure-like activity. The cilia-driven flow disturbance arises from altered cilia beating orientation with the disrupted polarity of the cilia anchoring rootlet meshwork. Together these findings indicate motile cilia disturbances have an essential role in CDD-associated seizures and beyond, suggesting cilia regulating kinases may be a therapeutic target for medication-resistant epilepsy.

Augmented Reality-Assisted CT-Guided Puncture: A Phantom Study

Purpose

To investigate the feasibility of a novel augmented reality system for CT-guided liver interventions and to compare it with free-hand interventions in a phantom setting.

Methods and materials

A newly developed augmented reality interface was used, with projection of CT-imaging in multiplanar reconstruction and live rendering of the needle position, a bull`s eye view of the needle trajectory and a visualization of the distance to the target. Punctures were performed on a custom-made abdominal phantom by three interventional radiologists with different levels of expertise. Time and needle placement accuracy were measured. Two-tailed Wilcoxon signed rank test (p < 0.05) was performed to evaluate intraparticipant difference.

Results

Intraparticipant puncture times were significantly shorter for each operator in the augmented reality condition (< 0.001 for the resident, < 0.001 for the junior staff member and 0.027 for the senior staff member). The junior staff member had an improvement in accuracy of 1 mm using augmented reality (p 0.026); the other two participants showed no significant improvement regarding accuracy.

Conclusion

In this small series, it appears that the novel augmented reality system may improve the speed of CT-guided punctures in the phantom model compared to the free-hand procedure while maintaining a similar accuracy.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00270-022-03195-y.

Automated high-speed 3D imaging of organoid cultures with multi-scale phenotypic quantification

Current imaging approaches limit the ability to perform multi-scale characterization of three-dimensional (3D) organotypic cultures (organoids) in large numbers. Here, we present an automated multi-scale 3D imaging platform synergizing high-density organoid cultures with rapid and live 3D single-objective light-sheet imaging. It is composed of disposable microfabricated organoid culture chips, termed JeWells, with embedded optical components and a laser beam-steering unit coupled to a commercial inverted microscope. It permits streamlining organoid culture and high-content 3D imaging on a single user-friendly instrument with minimal manipulations and a throughput of 300 organoids per hour. We demonstrate that the large number of 3D stacks that can be collected via our platform allows training deep learning-based algorithms to quantify morphogenetic organizations of organoids at multi-scales, ranging from the subcellular scale to the whole organoid level. We validated the versatility and robustness of our approach on intestine, hepatic, neuroectoderm organoids and oncospheres.

Tracking the progress of inflammation with PET/MRI in a canine model of myocardial infarction

BACKGROUND

Following myocardial infarction, tissue undergoes pathophysiological changes involving inflammation and scar tissue formation. However, little is known about the pathophysiology and prognostic significance of any corresponding changes in remote myocardium. The aim of this study was to investigate the potential application of a combined constant infusion of 18F-FDG and Gd-DTPA to quantitate inflammation and extracellular volume (ECV) from 3 to 40 days after myocardial infarction.

METHODS

Eight canine subjects were imaged at multiple time points following induction of an MI with a 60-minute concurrent constant infusion of Gd-DTPA and 18F-FDG using a hybrid PET/MRI scanner.

RESULTS

There was a significant increase in ECV in remote myocardium on day 14 post-MI (P = .034) and day 21 (P = .021) compared to the baseline. ECV was significantly elevated in the infarcted myocardium compared to remote myocardium at all time points post-MI (days 3, 7, 14, 21, and 40) (P < .001) while glucose uptake was also increased within the infarct on days 3, 7, 14, and 21 but not 40.

CONCLUSIONS

The significant increase in ECV in remote tissue may be due to an ongoing inflammatory process in the early weeks post-infarct.

Repeatability and robustness of MP-GRASP T1 mapping

PURPOSE

To demonstrate the repeatability of fast 3D T1 mapping using Magnetization-Prepared Golden-angle RAdial Sparse Parallel (MP-GRASP) MRI and its robustness to variation of imaging parameters including flip angle and spatial resolution in phantoms and the brain.

THEORY AND METHODS

Multiple imaging experiments were performed to (1) assess the robustness of MP-GRASP T1 mapping to B1 inhomogeneity using a single tube phantom filled with uniform MnCl2 liquid; (2) compare the repeatability of T1 mapping between MP-GRASP and inversion recovery-based spin-echo (IR-SE; over 12 scans), using a commercial T1MES phantom; (3) evaluate the longitudinal variation of T1 estimation using MP-GRASP with varying imaging parameters, including spatial resolution, flip angle, TR/TE, and acceleration rate, using the T1MES phantom (106 scans performed over a period of 12 months); and (4) evaluate the variation of T1 estimation using MP-GRASP with varying imaging parameters in the brain (24 scans in a single visit). In addition, the accuracy of MP-GRASP T1 mapping was also validated against IR-SE by performing linear correlation and calculating the Lin's concordance correlation coefficient (CCC).

RESULTS

MP-GRASP demonstrates good robustness to B1 inhomogeneity, with intra-slice variability below 1% in the single tube phantom experiment. The longitudinal variability is good both in the phantom (below 2.5%) and in the brain (below 2%) with varying imaging parameters. The T1 values estimated from MP-GRASP are accurate compared to that from the IR-SE imaging (R2 = 0.997, Lin's CCC = 0.996).

CONCLUSION

MP-GRASP shows excellent repeatability of T1 estimation over time, and it is also robust to variation of different imaging parameters evaluated in this study.

A non-invasive, automated diagnosis of Menière's disease using radiomics and machine learning on conventional magnetic resonance imaging: A multicentric, case-controlled feasibility study

Purpose

This study investigated the feasibility of a new image analysis technique (radiomics) on conventional MRI for the computer-aided diagnosis of Menière’s disease.

Materials and methods

A retrospective, multicentric diagnostic case–control study was performed. This study included 120 patients with unilateral or bilateral Menière’s disease and 140 controls from four centers in the Netherlands and Belgium. Multiple radiomic features were extracted from conventional MRI scans and used to train a machine learning-based, multi-layer perceptron classification model to distinguish patients with Menière’s disease from controls. The primary outcomes were accuracy, sensitivity, specificity, positive predictive value, and negative predictive value of the classification model.

Results

The classification accuracy of the machine learning model on the test set was 82%, with a sensitivity of 83%, and a specificity of 82%. The positive and negative predictive values were 71%, and 90%, respectively.

Conclusion

The multi-layer perceptron classification model yielded a precise, high-diagnostic performance in identifying patients with Menière’s disease based on radiomic features extracted from conventional T2-weighted MRI scans. In the future, radiomics might serve as a fast and noninvasive decision support system, next to clinical evaluation in the diagnosis of Menière’s disease.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11547-021-01425-w.

Determining the longitudinal accuracy and reproducibility of T1 and T2 in a 3T MRI scanner

Purpose

To determine baseline accuracy and reproducibility of T₁ and T₂ relaxation times over 12 months on a dedicated radiotherapy MRI scanner.

Methods

An International Society of Magnetic Resonance in Medicine/National Institute of Standards and Technology (ISMRM/NIST) System Phantom was scanned monthly on a 3T MRI scanner for 1 year. T₁ was measured using inversion recovery (T₁‐IR) and variable flip angle (T₁‐VFA) sequences and T₂ was measured using a multi‐echo spin echo (T₂‐SE) sequence. For each vial in the phantom, accuracy errors (%bias) were determined by the relative differences in measured T₁ and T₂ times compared to reference values. Reproducibility was measured by the coefficient of variation (CV) of T₁ and T₂ measurements across monthly scans. Accuracy and reproducibility were mainly assessed on vials with relaxation times expected to be in physiological ranges at 3T.

Results

A strong linear correlation between measured and reference relaxation times was found for all sequences tested (R² > 0.997). Baseline bias (and CV[%]) for T₁‐IR, T₁‐VFA and T₂‐SE sequences were +2.0% (2.1), +6.5% (4.2), and +8.5% (1.9), respectively.

Conclusions

The accuracy and reproducibility of T₁ and T₂ on the scanner were considered sufficient for the sequences tested. No longitudinal trends of variation were deduced, suggesting less frequent measurements are required following the establishment of baselines.

Rapid high-resolution volumetric T1 mapping using a highly accelerated stack-of-stars Look Locker technique

PURPOSE

To develop a fast volumetric T1 mapping technique.

MATERIALS AND METHODS

A stack-of-stars (SOS) Look Locker technique based on the acquisition of undersampled radial data (>30× relative to Nyquist) and an efficient multi-slab excitation scheme is presented. A principal-component based reconstruction is used to reconstruct T1 maps. Computer simulations were performed to determine the best choice of partitions per slab and degree of undersampling. The technique was validated in phantoms against reference T1 values measured with a 2D Cartesian inversion-recovery spin-echo technique. The SOS Look Locker technique was tested in brain (n = 4) and prostate (n = 5). Brain T1 mapping was carried out with and without kz acceleration and results between the two approaches were compared. Prostate T1 mapping was compared to standard techniques. A reproducibility study was conducted in brain and prostate. Statistical analyses were performed using linear regression and Bland Altman analysis.

RESULTS

Phantom T1 values showed excellent correlations between SOS Look Locker and the inversion-recovery spin-echo reference (r2 = 0.9965; p < 0.0001) and between SOS Look Locker with slab-selective and non-slab selective inversion pulses (r2 = 0.9999; p < 0.0001). In vivo results showed that full brain T1 mapping (1 mm3) with kz acceleration is achieved in 4 min 21 s. Full prostate T1 mapping (0.9 × 0.9 × 4 mm3) is achieved in 2 min 43 s. T1 values for brain and prostate were in agreement with literature values. A reproducibility study showed coefficients of variation in the range of 0.18-0.2% (brain) and 0.15-0.18% (prostate).

CONCLUSION

A rapid volumetric T1 mapping technique was developed. The technique enables high-resolution T1 mapping with adequate anatomical coverage in a clinically acceptable time.

Sparse precontrast T1 mapping for high-resolution whole-brain DCE-MRI

PURPOSE

To develop and evaluate an efficient precontrast T₁ mapping technique suitable for quantitative high-resolution whole-brain dynamic contrast-enhanced-magnetic resonance imaging (DCE-MRI).

METHODS

Variable flip angle (VFA) T₁ mapping was considered that provides 1 × 1 × 2 mm³ resolution to match a recent high-resolution whole-brain DCE-MRI protocol. Seven FAs were logarithmically spaced from 1.5° to 15°. T₁ and M₀ maps were estimated using model-based reconstruction. This approach was evaluated using an anatomically realistic brain tumor digital reference object (DRO) with noise-mimicking 3T neuroimaging and fully sampled data acquired from one healthy volunteer. Methods were also applied on fourfold prospectively undersampled VFA data from 13 patients with high-grade gliomas.

RESULTS

T₁ -mapping precision decreased with undersampling factor R, althoughwhereas bias remained small before a critical R. In the noiseless DRO, T₁ bias was <25 ms in white matter (WM) and <11 ms in brain tumor (BT). T₁ standard deviation (SD) was <119.5 ms in WM (coefficient of variation [COV] ~11.0%) and <253.2 ms in BT (COV ~12.7%). In the noisy DRO, T₁ bias was <50 ms in WM and <30 ms in BT. For R ≤ 10, T₁ SD was <107.1 ms in WM (COV ~9.9%) and <240.9 ms in BT (COV ~12.1%). In the healthy subject, T₁ bias was <30 ms for R ≤ 16. At R = 4, T₁ SD was 171.4 ms (COV ~13.0%). In the prospective brain tumor study, T₁ values were consistent with literature values in WM and BT.

CONCLUSION

High-resolution whole-brain VFA T₁ mapping is feasible with sparse sampling, supporting its use for quantitative DCE-MRI.

Arterial spin labelling as a gadolinium-free alternative in the detection of prostate cancer

PURPOSE

To determine the capability of Gadolinium-free arterial spin labelling (ASL) sequences as novel, contrast-free, non-invasive alternative perfusion imaging method to differentiate prostate cancer (PCA) from benign prostate tissue compared to conventional DCE MRI.

METHODS

Thirty men with histologically confirmed PCA were included in this prospectively enrolled single center cohort study. All patients received multiparametric MRI (T2, DWI, DCE) at 3 T with additional ASL of the PCA lesion. Primary endpoint was differentiability of PCA versus benign prostate tissue by signal intensities (SI) and contrast ratios (CR) in ASL in comparison to DCE. For DCE also Signal-Enhancement-Ratio (SER) of native and early contrast enhancement SI was assessed. Secondary objectives were differences regarding PCA localisation in peripheral (PZ) or transition zone (TZ) and PCA detection.

RESULTS

In both, ASL and DCE, average SI of PCA differed significantly from SI in benign tissue in the TZ and PZ (p < 0,01, respectively). ASL had significantly higher CR discerning PCA and benign tissue in PZ and TZ (PZ = 5.19; TZ = 6.45) compared to DCE SI (PZ = 1.61; TZ = 1.43) and DCE SER (PZ = 1.59; TZ = 1.43) (p < 0.01, respectively). In subjective evaluation, PCA could be detected in ASL in 28 patients, compared to 29 in DCE.

CONCLUSION

ASL had significantly higher CR differentiating PCA from benign tissue in PZ and TZ compared to DCE. Visual detection of PCA does not differ significantly between the two sequences. As perfusion gadolinium-based contrast media is seen more critical in the last few years, ASL seems to be a promising alternative to DCE in PCA detection.

MRI Mapping of Renal T1: Basic Concept

In renal MRI, measurement of the T₁ relaxation time of water molecules may provide a valuable biomarker for a variety of pathological conditions. Due to its sensitivity to the tissue microenvironment, T₁ has gained substantial interest for noninvasive imaging of renal pathology, including inflammation and fibrosis. In this chapter, we will discuss the basic concept of T₁ mapping and different T₁ measurement techniques and we will provide an overview of emerging preclinical applications of T₁ for imaging of kidney disease.This chapter is based upon work from the COST Action PARENCHIMA, a community-driven network funded by the European Cooperation in Science and Technology (COST) program of the European Union, which aims to improve the reproducibility and standardization of renal MRI biomarkers. This introduction chapter is complemented by two separate chapters describing the experimental procedure and data analysis.

Quantitative MRI: Defining repeatability, reproducibility and accuracy for prostate cancer imaging biomarker development

INTRODUCTION

Quantitative MRI (qMRI) parameters have been increasingly used to develop predictive models to accurately monitor treatment response in prostate cancer after radiotherapy. To reliably detect changes in signal due to treatment response, predictive models require qMRI parameters with high repeatability and reproducibility. The purpose of this study was to measure qMRI parameter uncertainties in both commercial and in-house developed phantoms to guide the development of robust predictive models for monitoring treatment response.

MATERIALS AND METHODS

ADC, T1, and R2* values were acquired across three 3 T scanners with a prostate-specific qMRI protocol using the NIST/ISMRM system phantom, RSNA/NIST diffusion phantom, and an in-house phantom. A B1 field map was acquired to correct for flip angle inhomogeneity in T1 maps. All sequences were repeated in each scan to assess within-session repeatability. Weekly scans were acquired on one scanner for three months with the in-house phantom. Between-session repeatability was measured with test-retest scans 6-months apart on all scanners with all phantoms. Accuracy, defined as percentage deviation from reference value for ADC and T1, was evaluated using the system and diffusion phantoms. Repeatability and reproducibility coefficients of variation (%CV) were calculated for all qMRI parameters on all phantoms.

RESULTS

Overall, repeatability CV of ADC was <2.40%, reproducibility CV was <3.98%, and accuracy ranged between -8.0% to 2.7% across all scanners. Applying B1 correction on T1 measurements significantly improved the repeatability and reproducibility (p<0.05) but increased error in accuracy (p<0.001). Repeatability and reproducibility of R2* was <4.5% and <7.3% respectively in the system phantom across all scanners.

CONCLUSION

Repeatability, reproducibility, and accuracy in qMRI parameters from a prostate-specific protocol was estimated using both commercial and in-house phantoms. Results from this work will be used to identify robust qMRI parameters for use in the development of predictive models to longitudinally monitor treatment response for prostate cancer in current and future clinical trials.

Characterization of 3-Dimensional Printing and Casting Materials for use in Magnetic Resonance Imaging Phantoms at 3 T

Imaging phantoms are used to calibrate and validate the performance of magnetic resonance imaging (MRI) systems. Many new materials have been developed for additive manufacturing (three-dimensional [3D] printing) processes that may be useful in the direct printing or casting of dimensionally accurate, anatomically accurate, patient-specific, and/or biomimetic MRI phantoms. The T1, T2, and T2* spin relaxation times of polymer samples were tested to discover materials for use as tissue mimics and structures in MRI phantoms. This study included a cohort of polymer compounds that was tested in cured form. The cohort consisted of 101 standardized polymer samples fabricated from: two-part silicones and polyurethanes used in commercial casting processes; one-part optically cured polyurethanes used in 3D printing; and fused deposition thermoplastics used in 3D printing. The testing was performed at 3 T using inversion recovery, spin echo, and gradient echo sequences for T1, T2, and T2*, respectively. T1, T2, and T2* values were plotted with error bars to allow the reader to assess how well a polymer matches a tissue for a specific application. A correlation was performed between T1, T2, T2* values and material density, elongation, tensile strength, and hardness. Two silicones, SI_XP-643 and SI_P-45, may be usable mimics for reported liver values; one silicone, SI_XP-643, may be a useful mimic for muscle; one silicone, SI_XP-738, may be a useful mimic for white matter; and four silicones, SI_P-15, SI_GI-1000, SI_GI-1040, and SI_GI-1110, may be usable mimics for spinal cord. Elongation correlated to T2 (p = 0.0007), tensile strength correlated to T1 (p = 0.002), T2 (p = 0.0003), and T2* (p = 0.003). The 80 samples not providing measurable signal with T1, T2, T2* relaxation values too short to measure with the standard sequences, may be useful for MRI-invisible fixturing and medical devices at 3 T.

Virtual measurements of paracellular permeability and chronic inflammation via color coded pixel-wise T1 mapping

To assess whether quantitative T₁ relaxometry can measure permeability, chronic inflammation and mural thickening of mouse bladder wall. Adult female C57BL6 mice unexposed to radiation (controls) or 40 wk postirradiation of 10 Gy were scanned at 9.4 T before and after instillation (0.1 mL) of aqueous, novel contrast mixture (NCM) containing 4 mM gadobutrol and 5 mM ferumoxytol. Rapid acquisition with refocused echo (RARE) sequence was used with variable repetition times (TR). Pixel-wise maps of T₁ relaxation times for the segmented bladder wall layers were generated from voxel-wise, nonlinear least square data fitting of TR-dependent signal intensity acquired with TR array of 0.4-10 s followed by the histology of harvested bladder. Significant differences between precontrast and postcontrast T₁ (ΔT₁) were noted in urothelium and lamina propria of both groups but only in detrusor of irradiated group (P < 0.001; 2-way ANOVA). Nearly twofold higher gadobutrol permeability (550 ± 73 vs. 294 ± 160 μM; P < 0.01) derived as per 1/ΔT₁ = r1. [C] in urothelium of irradiated group. Inflammation and bladder wall thickening (0.75 ± 0. vs. 0.44 ± 0.08 mm; P < 0.001) predicted by MRI was subsequently confirmed by histology and altered expression of CD45 and zonula occludens-1 (ZO-1) relative to controls. NCM enhanced MRI relies on the retention of large molecular weight ferumoxytol in lumen for negative contrast, while permeation of the non-ionic, small molecular weight gadobutrol through ZO-1 generates positive contrast in bladder wall for virtual measurement of paracellular permeability and assessment of chronic inflammation in thin and distensible bladder wall, which is also defined by its variable shape and location within pelvis.

MRI-visible liquid crystal thermometer

PURPOSE

MRI parameters, such as T1 , T2 , and ADC, of tissue-mimicking materials in MRI phantoms can exhibit temperature dependence, and bore temperatures can vary over a 10°C range across different MRI systems. If this variation is not accurately corrected for, the quantitative nature of reference or phantom measurements is irrelevant. Available thermometers require opening the phantoms to probe the temperature, which can introduce contaminants that may affect the stability and accuracy of the phantom. An integrated, MRI-visible thermometer that can be read using typical imaging protocols is needed.

THEORY AND METHODS

An MRI-compatible thermometer was designed using liquid crystals (LCs) that exhibit rapid transitions between the LC cholesteric state and isotropic state in the room temperature range spanning 17°C to 23°C in 1.0°C increments. The LC thermometer was assessed visually and using superconducting quantum interference device magnetometry, NMR, and MRI techniques.

RESULTS

The signal generated from the LC thermometer was visible with spin-echo and gradient-echo MRI images. The LC state transition temperatures were visually referenced to a National Institute of Standards and Technology-traceable thermometer, and these LC state transitions were confirmed using superconducting quantum interference device magnetometry and NMR.

CONCLUSIONS

The LC MR-visible thermometer had measurable changes in relative signal with temperature, which were invariant to a variety of imaging sequences used.

Prostate magnetic resonance imaging technique

Multiparametric magnetic resonance (MR) imaging of the prostate is an excellent tool to detect clinically significant prostate cancer, and it has widely been incorporated into clinical practice due to its excellent tissue contrast and image resolution. The aims of this article are to describe the prostate MR imaging technique for detection of clinically significant prostate cancer according to PI-RADS v2.1, as well as alternative sequences and basic aspects of patient preparation and MR imaging artifact avoidance.

Blood-brain barrier disruption and delivery of irinotecan in a rat model using a clinical transcranial MRI-guided focused ultrasound system

We investigated controlled blood-brain barrier (BBB) disruption using a low-frequency clinical transcranial MRI-guided focused ultrasound (TcMRgFUS) device and evaluated enhanced delivery of irinotecan chemotherapy to the brain and a rat glioma model. Animals received three weekly sessions of FUS, FUS and 10 mg/kg irinotecan, or irinotecan alone. In each session, four volumetric sonications targeted 36 locations in one hemisphere. With feedback control based on recordings of acoustic emissions, 98% of the sonication targets (1045/1071) reached a pre-defined level of acoustic emission, while the probability of wideband emission (a signature for inertial cavitation) was than 1%. BBB disruption, evaluated by mapping the R1 relaxation rate after administration of an MRI contrast agent, was significantly higher in the sonicated hemisphere (P < 0.01). Histological evaluation found minimal tissue effects. Irinotecan concentrations in the brain were significantly higher (P < 0.001) with BBB disruption, but SN-38 was only detected in <50% of the samples and only with an excessive irinotecan dose. Irinotecan with BBB disruption did not impede tumor growth or increase survival. Overall these results demonstrate safe and controlled BBB disruption with a low-frequency clinical TcMRgFUS device. While irinotecan delivery to the brain was not neurotoxic, it did not improve outcomes in the F98 glioma model.

Modified acquisition strategy for reduced motion artifact in super resolution T 2 FSE multislice MRI: Application to prostate

PURPOSE

To reduce slice-to-slice motion effects in multislice T 2 -weighted fast-spin-echo ( T 2 FSE) imaging, manifest as "scalloping" in reformats, by modification of the acquisition strategy and to show applicability in prostate MRI.

METHODS

T 2 FSE images of contiguous or overlapping slices are typically acquired using multiple passes in which each pass is comprised of multiple slices with slice-to-slice gaps. Combination of slices from all passes provides the desired sampling. For enhancement of through-plane resolution with super resolution or for reformatting into other orientations, subtle ≈1 mm motion between passes can cause objectionable "scalloping" artifact. Here we address this by subdivision of each pass into multiple segments. Interleaving of segments from the multiple passes causes all slices to be acquired over substantially the same time, reducing pass-to-pass motion effects. This was implemented in acquiring 78 overlapped T 2 FSE axial slices and studied in phantoms and in 14 prostate MRI patients. Super-resolution axial images and sagittal reformats from the original and new segmented acquisitions were evaluated by 3 uroradiologists.

RESULTS

For all criteria of sagittal reformats, the segmented acquisition was statistically superior to the original. For all sharpness criteria of axial images, although the trend preferred the original acquisition, the difference was not significant. For artifact in axial images, the segmented acquisition was significantly superior.

CONCLUSIONS

For prostate MRI the new segmented acquisition significantly reduces the scalloping motion artifact that can be present in reformats due to long time lags between the acquisition of adjacent or overlapped slices while retaining image sharpness in the acquired axial slices.

Evaluation of T1 relaxation time in prostate cancer and benign prostate tissue using a Modified Look-Locker inversion recovery sequence

Purpose of this study was to evaluate the diagnostic performance of T1 relaxation time (T1) for differentiating prostate cancer (PCa) from benign tissue as well as high- from low-grade PCa. Twenty-three patients with suspicion for PCa were included in this prospective study. 3 T MRI including a Modified Look-Locker inversion recovery sequence was acquired. Subsequent targeted and systematic prostate biopsy served as a reference standard. T1 and apparent diffusion coefficient (ADC) value in PCa and reference regions without malignancy as well as high- and low-grade PCa were compared using the Mann-Whitney U test. The performance of T1, ADC value, and a combination of both to differentiate PCa and reference regions was assessed by receiver operating characteristic (ROC) analysis. T1 and ADC value were lower in PCa compared to reference regions in the peripheral and transition zone (p < 0.001). ROC analysis revealed high AUCs for T1 (0.92; 95%-CI, 0.87-0.98) and ADC value (0.97; 95%-CI, 0.94 to 1.0) when differentiating PCa and reference regions. A combination of T1 and ADC value yielded an even higher AUC. The difference was statistically significant comparing it to the AUC for ADC value alone (p = 0.02). No significant differences were found between high- and low-grade PCa for T1 (p = 0.31) and ADC value (p = 0.8). T1 relaxation time differs significantly between PCa and benign prostate tissue with lower T1 in PCa. It could represent an imaging biomarker for PCa.

Selection of Fitting Model and Arterial Input Function for Repeatability in Dynamic Contrast-Enhanced Prostate MRI

RATIONALE AND OBJECTIVES

Analysis of dynamic contrast-enhanced (DCE) magnetic resonance imaging is notable for the variability of calculated parameters. The purpose of this study was to evaluate the level of measurement variability and error/variability due to modeling in DCE magnetic resonance imaging parameters.

MATERIALS AND METHODS

Two prostate DCE scans were performed on 11 treatment-naïve patients with suspected or confirmed prostate peripheral zone cancer within an interval of less than two weeks. Tumor-suspicious and normal-appearing regions of interest (ROI) in the prostate peripheral zone were segmented. Different Tofts-Kety based models and different arterial input functions, with and without bolus arrival time (BAT) correction, were used to extract pharmacokinetic parameters. The percent repeatability coefficient (%RC) of fitted model parameters K_trans, v_e, and k_ep was calculated. Paired t-tests comparing parameters in tumor-suspicious ROIs and in normal-appearing tissue evaluated each parameter's sensitivity to pathology.

RESULTS

Although goodness-of-fit criteria favored the four-parameter extended Tofts-Kety model with the BAT correction included, the simplest two-parameter Tofts-Kety model overall yielded the best repeatability scores. The best %RC in the tumor-suspicious ROI was 63% for k_ep, 28% for v_e, and 83% for K_trans . The best p values for discrimination between tissues were p <10^-5 for k_ep and K_trans, and p = 0.11 for v_e. Addition of the BAT correction to the models did not improve repeatability.

CONCLUSION

The parameter k_ep, using an arterial input functions directly measured from blood signals, was more repeatable than K_trans. Both K_trans and k_ep values were highly discriminatory between healthy and diseased tissues in all cases. The parameter v_e had high repeatability but could not distinguish the two tissue types.

Acoustic feedback enables safe and reliable carboplatin delivery across the blood-brain barrier with a clinical focused ultrasound system and improves survival in a rat glioma model

The blood-brain barrier (BBB) restricts delivery of most chemotherapy agents to brain tumors. Here, we investigated a clinical focused ultrasound (FUS) device to disrupt the BBB in rats and enhance carboplatin delivery to the brain using the F98 glioma model.

Methods: In each rat, 2-3 volumetric sonications (5 ms bursts at 1.1 Hz for 75s) targeted 18-27 locations in one hemisphere. Sonication was combined with Definity microbubbles (10 µl/kg) and followed by intravenous carboplatin (50 mg/kg). Closed-loop feedback control was performed based on acoustic emissions analysis.

Results: Safety and reliability were established in healthy rats after three sessions with carboplatin; BBB disruption was induced in every target without significant damage evident in MRI or histology. In tumor-bearing rats, concentrations of MRI contrast agent (Gadavist) were 1.7 and 3.3 times higher in the tumor center and margin, respectively, than non-sonicated tumors (P<0.001). Tissue-to-plasma ratios of intact carboplatin concentrations were increased by 7.3 and 2.9 times in brain and tumor respectively, at one hour after FUS and 4.2 and 2.4 times at four hours. Tumor volume doubling time in rats receiving FUS and carboplatin increased by 96% and 126% compared to rats that received carboplatin alone and non-sonicated controls, respectively (P<0.05); corresponding increases in median survival were 48% and 66% (P<0.01).

Conclusion: Overall, this work demonstrates that actively-controlled BBB disruption with a clinical device can enhance carboplatin delivery without neurotoxicity at level that reduces tumor growth and improves survival in an aggressive and infiltrative rat glioma model.

Assessment of histotripsy-induced liquefaction with diagnostic ultrasound and magnetic resonance imaging in vitro and ex vivo

Histotripsy is a therapeutic ultrasound modality under development to liquefy tissue mechanically via bubble clouds. Image guidance of histotripsy requires both quantification of the bubble cloud activity and accurate delineation of the treatment zone. In this study, magnetic resonance (MR) and diagnostic ultrasound imaging were combined to assess histotripsy treatment in vitro and ex vivo. Mechanically ablative histotripsy pulses were applied to agarose phantoms or porcine livers. Bubble cloud emissions were monitored with passive cavitation imaging (PCI), and hyperechogenicity via plane wave imaging. Changes in the medium structure due to bubble activity were assessed with diagnostic ultrasound using conventional B-mode imaging and T ₁-, T ₂-, and diffusion-weighted MR images acquired at 3 Tesla. Liquefaction zones were correlated with diagnostic ultrasound and MR imaging via receiver operating characteristic (ROC) analysis and Dice similarity coefficient (DSC) analysis. Diagnostic ultrasound indicated strong bubble activity for all samples. Histotripsy-induced changes in sample structure were evident on conventional B-mode and T ₂-weighted images for all samples, and were dependent on the sample type for T ₁- and diffusion-weighted imaging. The greatest changes observed on conventional B-mode or MR imaging relative to baseline in the samples did not necessarily indicate the regions of strongest bubble activity. Areas under the ROC curve for predicting phantom or liver liquefaction were significantly greater than 0.5 for PCI power, plane wave and conventional B-mode grayscale, T ₁, T ₂, and ADC. The acoustic power mapped via PCI provided a better prediction of liquefaction than assessment of the liquefaction zone via conventional B-mode or MR imaging for all samples. The DSC values for T ₂-weighted images were greater than those derived from conventional B-mode images. These results indicate diagnostic ultrasound and MR imaging provide complimentary sets of information, demonstrating that multimodal imaging is useful for assessment of histotripsy liquefaction.

Radiofrequency transmit calibration: A multi-center evaluation of vendor-provided radiofrequency transmit mapping methods

PURPOSE

To determine the accuracy and test-retest repeatability of fast radiofrequency (RF) transmit measurement approaches used in Dynamic Contrast Enhanced Magnetic Resonance Imaging (DCE-MRI). Spatial variation in the transmitted RF field introduces bias and increased variance in quantitative DCE-MRI metrics including tracer kinetic parameter maps. If unaccounted for, these errors can dominate all other sources of bias and variance. The amount and pattern of variation depend on scanner-specific hardware and software.

METHODS

Human tissue mimicking torso and brain phantoms were constructed. RF transmit maps were measured and compared across eight different commercial scanners, from three major vendors, and three clinical sites. Vendor-recommended rapid methods for RF mapping were compared to a slower reference method. Imaging was repeated at all sites after 2 months. Ranges and magnitude of RF inhomogeneity were compared scanner-wise at two time points. Limits of Agreement of vendor-recommended methods and double-angle reference method were assessed.

RESULTS

At 3 T, B₁ ⁺ inhomogeneity spans across 35% in the head and 120% in the torso. Fast vendor provided methods are within 30% agreement with the reference double angle method for both the head and the torso phantom.

CONCLUSIONS

If unaccounted for, B₁ ⁺ inhomogeneity can severely impact tracer-kinetic parameter estimation. Depending on the scanner, fast vendor provided B₁ ⁺ mapping sequences allow unbiased and reproducible measurements of B₁ ⁺ inhomogeneity to correct for this source of bias.

MRI-guided transurethral insonation of silica-shell phase-shift emulsions in the prostate with an advanced navigation platform

PURPOSE

In this study, the efficacy of transurethral prostate ablation in the presence of silica-shell ultrasound-triggered phase-shift emulsions (sUPEs) doped with MR contrast was evaluated. The influence of sUPEs on MR imaging assessment of the ablation zone was also investigated.

METHODS

sUPEs were doped with a magnetic resonance (MR) contrast agent, Gd2 O3 , to assess ultrasound transition. Injections of saline (sham), saline and sUPEs alone, and saline and sUPEs with Optison microbubbles were performed under guidance of a prototype interventional MRI navigation platform in a healthy canine prostate. Treatment arms were evaluated for differences in lesion size, T1  contrast, and temperature. In addition, non-perfused areas (NPAs) on dynamic contrast-enhanced (DCE) MRI, 55°C isotherms, and areas of 240 cumulative equivalent minutes at 43°C (CEM43 ) dose or greater computed from MR thermometry were measured and correlated with ablated areas indicated by histology.

RESULTS

For treatment arms including sUPEs, the computed correlation coefficients between the histological ablation zone and the NPA, 55°C isotherm, and 240 CEM43 area ranged from 0.96-0.99, 0.98-0.99, and 0.91-0.99, respectively. In the absence of sUPEs, the computed correlation coefficients between the histological ablation zone and the NPA, 55°C isotherm, and 240 CEM43 area were 0.69, 0.54, and 0.50, respectively. Across all treatment arms, the areas of thermal tissue damage and NPAs were not significantly different (P = 0.47). Areas denoted by 55°C isotherms and 240 CEM43 dose boundaries were significantly larger than the areas of thermal damage, again for all treatment arms (P = 0.009 and 0.003, respectively). No significant differences in lesion size, T1 contrast, or temperature were observed between any of the treatment arms (P > 0.0167). Lesions exhibiting thermal fixation on histological analysis were present in six of nine insonations involving sUPE injections and one of five insonations involving saline sham injections. Significantly larger areas (P = 0.002), higher temperatures (P = 0.004), and more frequent ring patterns of restricted diffusion on ex vivo diffusion-weighted imaging (P = 0.005) were apparent in lesions with thermal fixation.

CONCLUSIONS

T1 contrast suggesting sUPE transition was not evident in sUPE treatment arms. The use of MR imaging metrics to predict prostate ablation was not diminished by the presence of sUPEs. Lesions generated in the presence of sUPEs exhibited more frequent thermal fixation, though there were no significant changes in the ablation areas when comparing arms with and without sUPEs. Thermal fixation corresponded to some qualitative imaging features.

Impact of contrast agent injection duration on dynamic contrast-enhanced MRI quantification in prostate cancer

The volume transfer constant K^trans , which describes the leakage of contrast agent (CA) from vasculature into tissue, is the most commonly reported quantitative parameter for dynamic contrast-enhanced (DCE-) MRI. However, the variation in reported K^trans values between studies from different institutes is large. One of the primary sources of uncertainty is quantification of the arterial input function (AIF). The aim of this study is to determine the influence of the CA injection duration on the AIF and tracer kinetic analysis (TKA) parameters (i.e. K^trans , k_ep and v_e ). Thirty-one patients with prostate cancer received two DCE-MRI examinations with an injection duration of 5 s in the first examination and a prolonged injection duration in the second examination, varying between 7.5 s and 30 s. The DCE examination was carried out on a 3.0 T MRI scanner using a transversal T₁ -weighted 3D spoiled gradient echo sequence (300 s duration, dynamic scan time of 2.5 s). Data of 29 of the 31 were further analysed. AIFs were determined from the phase signal in the left and right femoral arteries. K^trans , k_ep and v_e were estimated with the standard Tofts model for regions of healthy peripheral zone and tumour tissue. We observed a significantly smaller peak height and increased width in the AIF for injection durations of 15 s and longer. However, we did not find significant differences in K^trans , k_ep or v_e for the studied injection durations. The study demonstrates that the TKA parameters K^trans , k_ep and v_e , measured in the prostate, do not show a significant change as a function of injection duration.

Magnetic resonance imaging and spectroscopy evidence of efficacy for adrenal and gonadal hormone replacement therapy in anorexia nervosa

PURPOSE

Dehydroepiandrosterone (DHEA)+estrogen/progestin therapy for adolescent girls with anorexia nervosa (AN) has the potential to arrest bone loss. The primary aim of this study was to test the effects of DHEA+estrogen/progestin therapy in adolescent girls with AN on bone marrow in the distal femur using magnetic resonance imaging (MRI) and spectroscopy.

METHODS

Seventy adolescent girls with AN were enrolled in a double blind, randomized, placebo-controlled trial at two urban hospital-based programs.

INTERVENTION

Seventy-six girls were randomly assigned to receive 12months of either oral micronized DHEA or placebo. DHEA was administered with conjugated equine estrogens (0.3mg daily) for 3months, then an oral contraceptive (20μg ethinyl estradiol/ 0.1mg levonorgestrel) for 9months. The primary outcome measure was bone marrow fat by MRI and magnetic resonance spectroscopy (MRS).

RESULTS

T2 of the water resonance dropped significantly less in the active vs. placebo group over 12months at both the medial and lateral distal femur (p=0.02). Body mass index (BMI) was a significant effect modifier for T1 and for T2 of unsaturated (T2_unsat) and saturated fat (T2_sat) in the lateral distal femur. Positive effects of the treatment of DHEA+estrogen/progestin were seen primarily for girls above a BMI of about 18kg/m².

CONCLUSIONS

These findings suggest treatment with oral DHEA+estrogen/progestin arrests the age- and disease-related changes in marrow fat composition in the lateral distal femur reported previously in this population.

Accuracy, repeatability, and interplatform reproducibility of T1 quantification methods used for DCE-MRI: Results from a multicenter phantom study

PURPOSE

To determine the in vitro accuracy, test-retest repeatability, and interplatform reproducibility of T₁ quantification protocols used for dynamic contrast-enhanced MRI at 1.5 and 3 T.

METHODS

A T₁ phantom with 14 samples was imaged at eight centers with a common inversion-recovery spin-echo (IR-SE) protocol and a variable flip angle (VFA) protocol using seven flip angles, as well as site-specific protocols (VFA with different flip angles, variable repetition time, proton density, and Look-Locker inversion recovery). Factors influencing the accuracy (deviation from reference NMR T₁ measurements) and repeatability were assessed using general linear mixed models. Interplatform reproducibility was assessed using coefficients of variation.

RESULTS

For the common IR-SE protocol, accuracy (median error across platforms = 1.4-5.5%) was influenced predominantly by T₁ sample (P < 10^-6 ), whereas test-retest repeatability (median error = 0.2-8.3%) was influenced by the scanner (P < 10^-6 ). For the common VFA protocol, accuracy (median error = 5.7-32.2%) was influenced by field strength (P = 0.006), whereas repeatability (median error = 0.7-25.8%) was influenced by the scanner (P < 0.0001). Interplatform reproducibility with the common VFA was lower at 3 T than 1.5 T (P = 0.004), and lower than that of the common IR-SE protocol (coefficient of variation 1.5T: VFA/IR-SE = 11.13%/8.21%, P = 0.028; 3 T: VFA/IR-SE = 22.87%/5.46%, P = 0.001). Among the site-specific protocols, Look-Locker inversion recovery and VFA (2-3 flip angles) protocols showed the best accuracy and repeatability (errors < 15%).

CONCLUSIONS

The VFA protocols with 2 to 3 flip angles optimized for different applications achieved acceptable balance of extensive spatial coverage, accuracy, and repeatability in T₁ quantification (errors < 15%). Further optimization in terms of flip-angle choice for each tissue application, and the use of B₁ correction, are needed to improve the robustness of VFA protocols for T₁ mapping. Magn Reson Med 79:2564-2575, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Dynamic contrast-enhanced magnetic resonance imaging of prostate cancer: A review of current methods and applications

In many areas of oncology, dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) has proven to be a clinically useful, non-invasive functional imaging technique to quantify tumor vasculature and tumor perfusion characteristics. Tumor angiogenesis is an essential process for tumor growth, proliferation, and metastasis. Malignant lesions demonstrate rapid extravasation of contrast from the intravascular space to the capillary bed due to leaky capillaries associated with tumor neovascularity. DCE-MRI has the potential to provide information regarding blood flow, areas of hypoperfusion, and variations in endothelial permeability and microvessel density to aid treatment selection, enable frequent monitoring during treatment and assess response to targeted therapy following treatment. This review will discuss the current status of DCE-MRI in cancer imaging, with a focus on its use in imaging prostate malignancies as well as weaknesses that limit its widespread clinical use. The latest techniques for quantification of DCE-MRI parameters will be reviewed and compared.

Robust and efficient pharmacokinetic parameter non-linear least squares estimation for dynamic contrast enhanced MRI of the prostate

PURPOSE

To describe an efficient numerical optimization technique using non-linear least squares to estimate perfusion parameters for the Tofts and extended Tofts models from dynamic contrast enhanced (DCE) MRI data and apply the technique to prostate cancer.

METHODS

Parameters were estimated by fitting the two Tofts-based perfusion models to the acquired data via non-linear least squares. We apply Variable Projection (VP) to convert the fitting problem from a multi-dimensional to a one-dimensional line search to improve computational efficiency and robustness. Using simulation and DCE-MRI studies in twenty patients with suspected prostate cancer, the VP-based solver was compared against the traditional Levenberg-Marquardt (LM) strategy for accuracy, noise amplification, robustness to converge, and computation time.

RESULTS

The simulation demonstrated that VP and LM were both accurate in that the medians closely matched assumed values across typical signal to noise ratio (SNR) levels for both Tofts models. VP and LM showed similar noise sensitivity. Studies using the patient data showed that the VP method reliably converged and matched results from LM with approximate 3× and 2× reductions in computation time for the standard (two-parameter) and extended (three-parameter) Tofts models. While LM failed to converge in 14% of the patient data, VP converged in the ideal 100%.

CONCLUSION

The VP-based method for non-linear least squares estimation of perfusion parameters for prostate MRI is equivalent in accuracy and robustness to noise, while being more reliably (100%) convergent and computationally about 3× (TM) and 2× (ETM) faster than the LM-based method.

Bone marrow fat content in 70 adolescent girls with anorexia nervosa: Magnetic resonance imaging and magnetic resonance spectroscopy assessment

BACKGROUND

Adolescents and women with anorexia nervosa have increased bone marrow fat and decreased bone formation, at least in part due to hormonal changes leading to preferential stem cell differentiation to adipocytes over osteoblasts.

OBJECTIVE

The purpose of this study was to evaluate marrow fat content and correlate with age and disease severity using knee MRI with T1 relaxometry (T1-R) and MR spectroscopy (MRS) in 70 adolescents with anorexia nervosa.

MATERIALS AND METHODS

We enrolled 70 girls with anorexia nervosa who underwent 3-T knee MRI with coronal T1-W images, T1-R and single-voxel proton MRS at 30 and 60 ms TE. Metaphyses were scored visually on the T1-W images for red marrow. Visual T1 score, T1 relaxometry values, MRS lipid indices and fat fractions were analyzed by regression on age, body mass index (BMI) and bone mineral density (BMD) as disease severity markers. MRS measures included unsaturated fat index, T2 water, unsaturated and saturated fat fractions.

RESULTS

All red marrow measures declined significantly with age. T1-R values were associated negatively with BMI and BMD for girls ≤16 years (P=0.03 and P=0.002, respectively) and positively for those≥17 years (P=0.05 and P=0.003, respectively). MRS identified a strong inverse association between T2 water and saturated fat fraction from 60 ms TE data (r=-0.85, P<0.0001). There was no association between unsaturated fat index and BMI or BMD.

CONCLUSIONS

The association between T1 and BMI and BMD among older girls suggests more marrow fat in those with severe anorexia nervosa. In contrast, the physiological association between marrow fat content and age remained dominant in younger patients. The strong association between T2 water and saturated fat may relate to the restricted mobility of water with increasing marrow fat.

T2-weighted MRI-derived textural features reflect prostate cancer aggressiveness: preliminary results

PURPOSE

To evaluate the diagnostic relevance of T2-weighted (T2W) MRI-derived textural features relative to quantitative physiological parameters derived from diffusion-weighted (DW) and dynamic contrast-enhanced (DCE) MRI in Gleason score (GS) 3+4 and 4+3 prostate cancers.

MATERIALS AND METHODS

3T multiparametric-MRI was performed on 23 prostate cancer patients prior to prostatectomy. Textural features [angular second moment (ASM), contrast, correlation, entropy], apparent diffusion coefficient (ADC), and DCE pharmacokinetic parameters (K^trans and V_e) were calculated from index tumours delineated on the T2W, DW, and DCE images, respectively. The association between the textural features and prostatectomy GS and the MRI-derived parameters, and the utility of the parameters in differentiating between GS 3+4 and 4+3 prostate cancers were assessed statistically.

RESULTS

ASM and entropy correlated significantly (p < 0.05) with both GS and median ADC. Contrast correlated moderately with median ADC. The textural features correlated insignificantly with K^trans and V_e. GS 4+3 cancers had significantly lower ASM and higher entropy than 3+4 cancers, but insignificant differences in median ADC, K^trans, and V_e. The combined texture-MRI parameters yielded higher classification accuracy (91%) than the individual parameter sets.

CONCLUSION

T2W MRI-derived textural features could serve as potential diagnostic markers, sensitive to the pathological differences in prostate cancers.

KEY POINTS

• T2W MRI-derived textural features correlate significantly with Gleason score and ADC. • T2W MRI-derived textural features differentiate Gleason score 3+4 from 4+3 cancers. • T2W image textural features could augment tumour characterization.

Dynamic contrast-enhanced case-control analysis in 3T MRI of prostate cancer can help to characterize tumor aggressiveness

PURPOSE

The aim of this work is to establish normality and tumor tissue ranges for perfusion parameters from dynamic contrast-enhanced (DCE) MR of the peripheral prostate at 3T and to compare the diagnostic performance of quantitative and semi-quantitative parameters.

MATERIALS AND METHODS

Thirty-six patients with prostate carcinomas (18 Gleason-6, 15 Gleason-7, and 3 Gleason-8) and 33 healthy subjects were included. Image analysis workflow comprised four steps: manual segmentation of whole prostate and lesions, series registration, voxelwise T1 mapping and calculation of pharmacokinetic and semi-quantitative parameters.

RESULTS

K^trans, v_e, upslope and AUC60 showed statistically significant differences between healthy peripheral areas and tumors. Curve type showed no association with healthy/tumor peripheral areas (chi-square=0.702). Areas under the ROC curves were 0.64 (95% CI: 0.54-0.75), 0.70 (0.60-0.80), 0.62 (0.51-0.72) and 0.63 (0.52-0.74) for K^trans, v_e, upslope and AUC60, respectively. The optimal cutoff values were: K^trans=0.21min^-1 (sensitivity=0.61, specificity=0.64), v_e=0.36 (0.63, 0.71), upslope=0.59 (0.59, 0.59) and AUC60=2.4 (0.63, 0.64). Significant differences were found between Gleason scores 6 and 7 for normalized K^trans, upslope and AUC60, with good diagnostic accuracy (area under ROC curve 0.80, 95% CI: 0.60-1.00).

CONCLUSION

Quantitative (K^trans and v_e) and semi-quantitative (upslope and AUC60) perfusion parameters showed significant differences between tumors and control areas in the peripheral prostate. Normalized K^trans, upslope and AUC60 values might characterize tumor aggressiveness.

Computer-aided Detection of Prostate Cancer with MRI: Technology and Applications

One in six men will develop prostate cancer in his lifetime. Early detection and accurate diagnosis of the disease can improve cancer survival and reduce treatment costs. Recently, imaging of prostate cancer has greatly advanced since the introduction of multiparametric magnetic resonance imaging (mp-MRI). Mp-MRI consists of T2-weighted sequences combined with functional sequences including dynamic contrast-enhanced MRI, diffusion-weighted MRI, and magnetic resonance spectroscopy imaging. Because of the big data and variations in imaging sequences, detection can be affected by multiple factors such as observer variability and visibility and complexity of the lesions. To improve quantitative assessment of the disease, various computer-aided detection systems have been designed to help radiologists in their clinical practice. This review paper presents an overview of literatures on computer-aided detection of prostate cancer with mp-MRI, which include the technology and its applications. The aim of the survey is threefold: an introduction for those new to the field, an overview for those working in the field, and a reference for those searching for literature on a specific application.

A novel alternative to classify tissues from T 1 and T 2 relaxation times for prostate MRI

OBJECTIVE

To segment and classify the different attenuation regions from MRI at the pelvis level using the T 1 and T 2 relaxation times and anatomical knowledge as a first step towards the creation of PET/MR attenuation maps.

MATERIALS AND METHODS

Relaxation times were calculated by fitting the pixel-wise intensities of acquired T 1- and T 2-weighted images from eight men with inversion-recovery and multi-echo multi-slice spin-echo sequences. A decision binary tree based on relaxation times was implemented to segment and classify fat, muscle, prostate, and air (within the body). Connected component analysis and an anatomical knowledge-based procedure were implemented to localize the background and bone.

RESULTS

Relaxation times at 3 T are reported for fat (T 1 = 385 ms, T 2 = 121 ms), muscle (T 1 = 1295 ms, T 2 = 40 ms), and prostate (T 1 = 1700 ms, T 2 = 80 ms). The relaxation times allowed the segmentation-classification of fat, prostate, muscle, and air, and combined with anatomical knowledge, they allowed classification of bone. The good segmentation-classification of prostate [mean Dice similarity score (mDSC) = 0.70] suggests a viable implementation in oncology and that of fat (mDSC = 0.99), muscle (mDSC = 0.99), and bone (mDSCs = 0.78) advocates for its implementation in PET/MR attenuation correction.

CONCLUSION

Our method allows the segmentation and classification of the attenuation-relevant structures required for the generation of the attenuation map of PET/MR systems in prostate imaging: air, background, bone, fat, muscle, and prostate.

Practical aspects of prostate MRI: hardware and software considerations, protocols, and patient preparation

The use of multiparametric MRI scans for the evaluation of men with prostate cancer has increased dramatically and is likely to continue expanding as new developments come to practice. However, it has not yet gained the same level of acceptance of other imaging tests. Partly, this is because of the use of suboptimal protocols, lack of standardization, and inadequate patient preparation. In this manuscript, we describe several practical aspects of prostate MRI that may facilitate the implementation of new prostate imaging programs or the expansion of existing ones.

Bolus arrival time and its effect on tissue characterization with dynamic contrast-enhanced magnetic resonance imaging

Matching the bolus arrival time (BAT) of the arterial input function (AIF) and tissue residue function (TRF) is necessary for accurate pharmacokinetic (PK) modeling of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). We investigated the sensitivity of volume transfer constant ([Formula: see text]) and extravascular extracellular volume fraction ([Formula: see text]) to BAT and compared the results of four automatic BAT measurement methods in characterization of prostate and breast cancers. Variation in delay between AIF and TRF resulted in a monotonous change trend of [Formula: see text] and [Formula: see text] values. The results of automatic BAT estimators for clinical data were all comparable except for one BAT estimation method. Our results indicate that inaccuracies in BAT measurement can lead to variability among DCE-MRI PK model parameters, diminish the quality of model fit, and produce fewer valid voxels in a region of interest. Although the selection of the BAT method did not affect the direction of change in the treatment assessment cohort, we suggest that BAT measurement methods must be used consistently in the course of longitudinal studies to control measurement variability.

Feasibility and reproducibility of BOLD and TOLD measurements in the liver with oxygen and carbogen gas challenge in healthy volunteers and patients with hepatocellular carcinoma

PURPOSE

To quantify baseline relaxation rates R2* and R1 in the abdomen, their changes after respiratory challenges, and their reproducibility in healthy volunteers and patients with hepatocellular carcinoma (HCC) at 1.5T and 3.0T.

MATERIALS AND METHODS

R2* measurements were acquired in the liver in 8 volunteers and 27 patients with 34 HCCs using multiecho T2* at baseline and after respiratory challenges with 100% oxygen (O2 ) and carbogen (CB = 95%O2 /5%CO2 ). R1 was measured at 1.5T in one volunteer and 21 patients with 23 HCCs. Test-retest coefficient of variation (CV) was assessed in 10 subjects. Intra- and interobserver variability of R2* and R1 measurements was assessed in 12 and 10 patients, respectively. Parameters for HCC, liver, and muscle were compared between baseline and after gas challenges.

RESULTS

We observed that R2* and R1 imaging of HCCs with O2 and CB is feasible and reproducible (test-retest CV R2*<15%/R1 <5%; intra- and interobserver intraclass correlation coefficient R2*>0.88/R1 >0.7 and CV R2*<7%/R1 <3% at 1.5T). R2* measurements were observed to be less reproducible at 3.0T (CV<35%). There was a statistically significant decrease in R2* values in HCC before and after O2 (P = 0.02) and increase in R1 after O2 (P = 0.004). CB had no significant effect (P R2* = 0.47/R1  = 0.278).

CONCLUSION

R2* measurements in HCC and liver parenchyma are more reproducible at 1.5T than at 3.0T, and with O2 than with CB challenge. We observed a decrease in R2* and an increase in R1 of HCCs from baseline in response to O2 challenge, as expected with increased tissue and blood oxygenation.

Quantitative pharmacokinetic analysis of prostate cancer DCE-MRI at 3T: comparison of two arterial input functions on cancer detection with digitized whole mount histopathological validation

Accurate pharmacokinetic (PK) modeling of dynamic contrast enhanced MRI (DCE-MRI) in prostate cancer (PCa) requires knowledge of the concentration time course of the contrast agent in the feeding vasculature, the so-called arterial input function (AIF). The purpose of this study was to compare AIF choice in differentiating peripheral zone PCa from non-neoplastic prostatic tissue (NNPT), using PK analysis of high temporal resolution prostate DCE-MRI data and whole-mount pathology (WMP) validation. This prospective study was performed in 30 patients who underwent multiparametric endorectal prostate MRI at 3.0T and WMP validation. PCa foci were annotated on WMP slides and MR images using 3D Slicer. Foci ≥0.5cm(3) were contoured as tumor regions of interest (TROIs) on subtraction DCE (early-arterial - pre-contrast) images. PK analyses of TROI and NNPT data were performed using automatic AIF (aAIF) and model AIF (mAIF) methods. A paired t-test compared mean and 90th percentile (p90) PK parameters obtained with the two AIF approaches. Receiver operating characteristic (ROC) analysis determined diagnostic accuracy (DA) of PK parameters. Logistic regression determined correlation between PK parameters and histopathology. Mean TROI and NNPT PK parameters were higher using aAIF vs. mAIF (p<0.05). There was no significant difference in DA between AIF methods: highest for p90 volume transfer constant (K(trans)) (aAIF differences in the area under the ROC curve (Az) = 0.827; mAIF Az=0.93). Tumor cell density correlated with aAIF K(trans) (p=0.03). Our results indicate that DCE-MRI using both AIF methods is excellent in discriminating PCa from NNPT. If quantitative DCE-MRI is to be used as a biomarker in PCa, the same AIF method should be used consistently throughout the study.

The role of pathology correlation approach in prostate cancer index lesion detection and quantitative analysis with multiparametric MRI

Rationale and Objectives

Development of imaging biomarkers often relies on their correlation with histopathology. Our aim was to compare two approaches for correlating pathology to multiparametric magnetic resonance (MR) imaging (mpMRI) for localization and quantitative assessment of prostate cancer (PCa) index tumor using whole mount (WM) pathology (WMP) as the reference.

Materials and Methods

Patients (N = 30) underwent mpMRI that included diffusion-weighted imaging and dynamic contrast-enhanced (DCE) MRI at 3 T before radical prostatectomy (RP). RP specimens were processed using WM technique (WMP) and findings summarized in a standard surgical pathology report (SPR). Histology index tumor volumes (HTVs) were compared to MR tumor volumes (MRTVs) using two approaches for index lesion identification on mpMRI using annotated WMP slides as the reference (WMP) and using routine SPR as the reference. Consistency of index tumor localization, tumor volume, and mean values of the derived quantitative parameters (mean apparent diffusion coefficient [ADC], K^trans, and v_e) were compared.

Results

Index lesions from 16 of 30 patients met the selection criteria. There was WMP/SRP agreement in index tumor in 13 of 16 patients. ADC-based MRTVs were larger (P < .05) than DCE-based MRTVs. ADC MRTVs were smaller than HTV (P < .005). There was a strong correlation between HTV and MRTV (Pearson r > 0.8; P < .05). No significant differences were observed in the mean values of K^trans and ADC between the WMP and SPR.

Conclusions

WMP correlation is superior to SPR for accurate localization of all index lesions. The use of WMP is however not required to distinguish significant differences of mean values of quantitative MRI parameters within tumor volume.