To investigate the repeatability and reproducibility of radiomic features extracted from MR images and provide a workflow to identify robust features. Methods: T 2-weighted images of a pelvic phantom were acquired on three scanners of two manufacturers and two magnetic field strengths. The repeatability and reproducibility of features were assessed by the intraclass correlation coefficient and the concordance correlation coefficient, respectively, and by the within-subject coefficient of variation, considering repeated acquisitions with and without phantom repositioning, and with different scanner and acquisition parameters. The features showing intraclass correlation coefficient or concordance correlation coefficient >0.9 were selected, and their dependence on shape information (Spearman's ρ > 0.8) analyzed. They were classified for their ability to distinguish textures, after shuffling voxel intensities of images. Results: From 944 two-dimensional features, 79.9% to 96.4% showed excellent repeatability in fixed position across all scanners. A much lower range (11.2% to 85.4%) was obtained after phantom repositioning. Three-dimensional extraction did not improve repeatability performance. Excellent reproducibility between scanners was observed in 4.6% to 15.6% of the features, at fixed imaging parameters. In addition, 82.4% to 94.9% of the features showed excellent agreement when extracted from images acquired with echo times 5 ms apart, but decreased with increasing 1714 | BIANCHINI et Al.

Purpose Mesorectal lymph node staging plays an important role in treatment decision making. Here, we explore the benefit of higher‐order diffusion MRI models accounting for non‐Gaussian diffusion effects to classify mesorectal lymph nodes both 1) ex vivo at ultrahigh field correlated with histology and 2) in vivo in a clinical scanner upon patient staging. Methods The preclinical investigation included 54 mesorectal lymph nodes, which were scanned at 16.4 T with an extensive diffusion MRI acquisition. Eight diffusion models were compared in terms of goodness of fit, lymph node classification ability, and histology correlation. In the clinical part of this study, 10 rectal cancer patients were scanned with diffusion MRI at 1.5 T, and 72 lymph nodes were analyzed with Apparent Diffusion Coefficient (ADC), Intravoxel Incoherent Motion (IVIM), Kurtosis, and IVIM‐Kurtosis. Results Compartment models including restricted and anisotropic diffusion improved the preclinical data fit, as well as the lymph node classification, compared to standard ADC. The comparison with histology revealed only moderate correlations, and the highest values were observed between diffusion anisotropy metrics and cell area fraction. In the clinical study, the diffusivity from IVIM‐Kurtosis was the only metric showing significant differences between benign (0.80 ± 0.30 μm2/ms) and malignant (1.02 ± 0.41 μm2/ms, P = .03) nodes. IVIM‐Kurtosis also yielded the largest area under the receiver operating characteristic curve (0.73) and significantly improved the node differentiation when added to the standard visual analysis by experts based on T2‐weighted imaging. Conclusion Higher‐order diffusion MRI models perform better than standard ADC and may be of added value for mesorectal lymph node classification in rectal cancer patients.

Noninvasive characterization of lymph node involvement in cancer is an enduring onerous challenge. In rectal cancer, pathologic lymph node status constitutes the most important determinant of local recurrence and overall survival, and patients with involved lymph nodes may benefit from preoperative chemo and/or radiotherapy. However, knowledge of lymph node status before surgery is currently hampered by limited imaging accuracy. Here, we introduce Susceptibility-Perturbation MRI (SPI) as a novel source of contrast to map malignant infiltration into mesorectal lymph nodes. SPI involves multigradient echo (MGE) signal decays presenting a nonmonoexponential nature, which we show is sensitive to the underlying microstructure via susceptibility perturbations. Using numerical simulations, we predicted that the large cell morphology and the high cellularity of tumor within affected mesorectal lymph nodes would induce signature SPI decays. We validated this prediction in mesorectal lymph nodes excised from total mesorectal excision specimens of patients with rectal cancer using ultrahigh field (16.4 T) MRI. SPI signals distinguished benign from malignant nodal tissue, both qualitatively and quantitatively, and our histologic analyses confirmed cellularity and cell size were the likely underlying sources for the differences observed. SPI was then adapted to a clinical 1.5 T scanner, added to patients' staging protocol, and compared with conventional assessment by two expert radiologists. Nonmonoexponential decays, similar to those observed in the ex vivo study, were demonstrated, and SPI classified lymph nodes more accurately than standard highresolution T 2 -weighted imaging assessment. These findings suggest this simple, yet highly informative, method can improve rectal cancer patient selection for neoadjuvant therapy.Significance: These findings introduce an MRI methodology tailored to detect magnetic susceptibility perturbations induced by subtle alterations in tissue microstructure.