Fully Integrated Quantitative Multiparametric Analysis of Non-Small Cell Lung Cancer at 3-T PET/MRI: Toward One-Stop-Shop Tumor Biological Characterization at the Supervoxel Level.

INTRODUCTION: The aim of this study was to study the feasibility of a fully integrated multiparametric imaging framework to characterize non-small cell lung cancer (NSCLC) at 3-T PET/MRI. PATIENTS AND METHODS: An 18F-FDG PET/MRI multiparametric imaging framework was developed and prospectively applied to 11 biopsy-proven NSCLC patients. For each tumor, 12 parametric maps were generated, including PET full kinetic modeling, apparent diffusion coefficient, T1/T2 relaxation times, and DCE full kinetic modeling. Gaussian mixture model-based clustering was applied at the whole data set level to define supervoxels of similar multidimensional PET/MRI behaviors. Taking the multidimensional voxel behaviors as input and the supervoxel class as output, machine learning procedure was finally trained and validated voxelwise to reveal the dominant PET/MRI characteristics of these supervoxels at the whole data set and individual tumor levels. RESULTS: The Gaussian mixture model-based clustering clustering applied at the whole data set level (17,316 voxels) found 3 main multidimensional behaviors underpinned by the 12 PET/MRI quantitative parameters. Four dominant PET/MRI parameters of clinical relevance (PET: k2, k3 and DCE: ve, vp) predicted the overall supervoxel behavior with 97% of accuracy (SD, 0.7; 10-fold cross-validation). At the individual tumor level, these dimensionality-reduced supervoxel maps showed mean discrepancy of 16.7% compared with the original ones. CONCLUSIONS: One-stop-shop PET/MRI multiparametric quantitative analysis of NSCLC is clinically feasible. Both PET and MRI parameters are useful to characterize the behavior of tumors at the supervoxel level. In the era of precision medicine, the full capabilities of PET/MRI would give further insight of the characterization of NSCLC behavior, opening new avenues toward image-based personalized medicine in this field.

Simultaneous multi-slice spin- and gradient-echo dynamic susceptibility-contrast perfusion-weighted MRI of gliomas.

Although combined spin- and gradient-echo (SAGE) dynamic susceptibility-contrast (DSC) MRI can provide perfusion quantification that is sensitive to both macrovessels and microvessels while correcting for T1 -shortening effects, spatial coverage is often limited in order to maintain a high temporal resolution for DSC quantification. In this work, we combined a SAGE echo-planar imaging (EPI) sequence with simultaneous multi-slice (SMS) excitation and blipped controlled aliasing in parallel imaging (blipped CAIPI) at 3 T to achieve both high temporal resolution and whole brain coverage. Two protocols using this sequence with multi-band (MB) acceleration factors of 2 and 3 were evaluated in 20 patients with treated gliomas to determine the optimal scan parameters for clinical use. ΔR2 *(t) and ΔR2 (t) curves were derived to calculate dynamic signal-to-noise ratio (dSNR), ΔR2 *- and ΔR2 -based relative cerebral blood volume (rCBV), and mean vessel diameter (mVD) for each voxel. The resulting SAGE DSC images acquired using MB acceleration of 3 versus 2 appeared visually similar in terms of image distortion and contrast. The difference in the mean dSNR from normal-appearing white matter (NAWM) and that in the mean dSNR between NAWM and normal-appearing gray matter were not statistically significant between the two protocols. ΔR2 *- and ΔR2 -rCBV maps and mVD maps provided unique contrast and spatial heterogeneity within tumors.

AZTEK: Adaptive zero TE k-space trajectories.

PURPOSE: Because of short signal lifetimes and respiratory motion, 3D lung MRI is still challenging today. Zero-TE (ZTE) pulse sequences offer promising solutions as they overcome the issue of short T2∗ . Nevertheless, as they rely on continuous readout gradients, the trajectories they follow in k-space are not adapted to retrospective gating and inferred motion correction. THEORY AND METHODS: We propose AZTEK (adaptive ZTE k-space trajectories), a set of 3D radial trajectories featuring three tuning parameters, to adapt the acquisition to any moving organ while keeping seamless transitions between consecutive spokes. Standard ZTE and AZTEK trajectories were compared for static and moving phantom acquisitions as well as for human thoracic imaging performed on 3 volunteers (1 healthy and 2 patients with lung cancer). RESULTS: For the static phantom, we observe comparable image qualities with standard and AZTEK trajectories. For the moving phantom, spatially coherent undersampling artifacts observed on gated images with the standard trajectory are alleviated with AZTEK. The same improvement in image quality is obtained in human, so details are more delineated in the lung with the use of the adaptive trajectory. CONCLUSION: The AZTEK technique opens the possibility for 3D dynamic ZTE lung imaging with retrospective gating. It enables us to uniformly sample the k-space for any arbitrary respiratory motion gate, while preserving static image quality, improving dynamic image quality and guaranteeing continuous readout gradient transitions between spokes, which makes it appropriate to ZTE.

18F-FDG PET and DCE kinetic modeling and their correlations in primary NSCLC: first voxel-wise correlative analysis of human simultaneous [18F]FDG PET-MRI data.

OBJECTIVES: To decipher the correlations between PET and DCE kinetic parameters in non-small-cell lung cancer (NSCLC), by using voxel-wise analysis of dynamic simultaneous [18F]FDG PET-MRI. MATERIAL AND METHODS: Fourteen treatment-naïve patients with biopsy-proven NSCLC prospectively underwent a 1-h dynamic [18F]FDG thoracic PET-MRI scan including DCE. The PET and DCE data were normalized to their corresponding T1-weighted MR morphological space, and tumors were masked semi-automatically. Voxel-wise parametric maps of PET and DCE kinetic parameters were computed by fitting the dynamic PET and DCE tumor data to the Sokoloff and Extended Tofts models respectively, by using in-house developed procedures. Curve-fitting errors were assessed by computing the relative root mean square error (rRMSE) of the estimated PET and DCE signals at the voxel level. For each tumor, Spearman correlation coefficients (rs) between all the pairs of PET and DCE kinetic parameters were estimated on a voxel-wise basis, along with their respective bootstrapped 95% confidence intervals (n = 1000 iterations). RESULTS: Curve-fitting metrics provided fit errors under 20% for almost 90% of the PET voxels (median rRMSE = 10.3, interquartile ranges IQR = 8.1; 14.3), whereas 73.3% of the DCE voxels showed fit errors under 45% (median rRMSE = 31.8%, IQR = 22.4; 46.6). The PET-PET, DCE-DCE, and PET-DCE voxel-wise correlations varied according to individual tumor behaviors. Beyond this wide variability, the PET-PET and DCE-DCE correlations were mainly high (absolute rs values > 0.7), whereas the PET-DCE correlations were mainly low to moderate (absolute rs values < 0.7). Half the tumors showed a hypometabolism with low perfused/vascularized profile, a hallmark of hypoxia, and tumor aggressiveness. CONCLUSION: A dynamic "one-stop shop" procedure applied to NSCLC is technically feasible in clinical practice. PET and DCE kinetic parameters assessed simultaneously are not highly correlated in NSCLC, and these correlations showed a wide variability among tumors and patients. These results tend to suggest that PET and DCE kinetic parameters might provide complementary information. In the future, this might make PET-MRI a unique tool to characterize the individual tumor biological behavior in NSCLC.

Diffusion-weighted Imaging Voxelwise-matched Analyses of Lung Cancer at 3.0-T PET/MRI: Reverse Phase Encoding Approach for Echo-planar Imaging Distortion Correction.

Background PET/MRI has drawn increasing interest in thoracic oncology due to the simultaneous acquisition of PET and MRI data. Geometric distortions related to diffusion-weighted imaging (DWI) limit the evaluation of voxelwise multimodal analyses. Purpose To assess the effectiveness of reverse phase encoding in correcting DWI geometric distortion for multimodal PET/MRI voxelwise lung tumor analyses. Materials and Methods In this prospective study, reverse phase encoding method was implemented with 3.0-T PET/MRI to correct geometric distortions related to DWI. The method was validated in dedicated phantom and then applied to 12 consecutive patients (mean age, 66 years ± 13 [standard deviation]; 10 men) suspected of having lung cancer who underwent fluorodeoxyglucose PET/MRI between October 2018 and April 2019. The effects on DWI-related image matching and apparent diffusion coefficient (ADC) regional map computation were assessed. Consequences on multimodal PET/MRI voxelwise lung tumor analyses were evaluated. Spearman correlation coefficients (rs) between the standardized uptake value (SUV) and ADC data corrected for distortion were computed from optimal realigned DWI PET data, along with bootstrap confidence intervals. Results Phantom results showed that in highly distorted areas, correcting the distortion significantly reduced the mean error against the ground truth (-25% ± 10.6 to -18.4% ± 12.6; P < .001) and the number of voxels with more than 20% error (from 85.3% to 31.4%). In the 12 patients, the coregistration of multimodal PET/MRI tumor data was improved by using the reverse phase encoding method (0.4%-44%). In all tumors, voxelwise correlations (rs) between ADC and SUV revealed null or weak monotonic relationships (mean rs of 0.016 ± 0.24 with none above 0.5). Conclusion Reverse phase encoding is a simple-to-implement method for improved diffusion-weighted multimodal PET/MRI voxelwise-matched analyses in lung cancer. © RSNA, 2020 Online supplemental material is available for this article. See also the editorial by Colletti in this issue.

Subject-specific bone attenuation correction for brain PET/MR: can ZTE-MRI substitute CT scan accurately?

In brain PET/MR applications, accurate attenuation maps are required for accurate PET image quantification. An implemented attenuation correction (AC) method for brain imaging is the single-atlas approach that estimates an AC map from an averaged CT template. As an alternative, we propose to use a zero echo time (ZTE) pulse sequence to segment bone, air and soft tissue. A linear relationship between histogram normalized ZTE intensity and measured CT density in Hounsfield units ([Formula: see text]) in bone has been established thanks to a CT-MR database of 16 patients. Continuous AC maps were computed based on the segmented ZTE by setting a fixed linear attenuation coefficient (LAC) to air and soft tissue and by using the linear relationship to generate continuous µ values for the bone. Additionally, for the purpose of comparison, four other AC maps were generated: a ZTE derived AC map with a fixed LAC for the bone, an AC map based on the single-atlas approach as provided by the PET/MR manufacturer, a soft-tissue only AC map and, finally, the CT derived attenuation map used as the gold standard (CTAC). All these AC maps were used with different levels of smoothing for PET image reconstruction with and without time-of-flight (TOF). The subject-specific AC map generated by combining ZTE-based segmentation and linear scaling of the normalized ZTE signal into [Formula: see text] was found to be a good substitute for the measured CTAC map in brain PET/MR when used with a Gaussian smoothing kernel of [Formula: see text] corresponding to the PET scanner intrinsic resolution. As expected TOF reduces AC error regardless of the AC method. The continuous ZTE-AC performed better than the other alternative MR derived AC methods, reducing the quantification error between the MRAC corrected PET image and the reference CTAC corrected PET image.