Automatic hepatic tumor segmentation in intra-operative ultrasound: a supervised deep-learning approach.

Purpose: Training and evaluation of the performance of a supervised deep-learning model for the segmentation of hepatic tumors from intraoperative US (iUS) images, with the purpose of improving the accuracy of tumor margin assessment during liver surgeries and the detection of lesions during colorectal surgeries. Approach: In this retrospective study, a U-Net network was trained with the nnU-Net framework in different configurations for the segmentation of CRLM from iUS. The model was trained on B-mode intraoperative hepatic US images, hand-labeled by an expert clinician. The model was tested on an independent set of similar images. The average age of the study population was 61.9 ± 9.9 years. Ground truth for the test set was provided by a radiologist, and three extra delineation sets were used for the computation of inter-observer variability. Results: The presented model achieved a DSC of 0.84 (p=0.0037), which is comparable to the expert human raters scores. The model segmented hypoechoic and mixed lesions more accurately (DSC of 0.89 and 0.88, respectively) than hyper- and isoechoic ones (DSC of 0.70 and 0.60, respectively) only missing isoechoic or >20 mm in diameter (8% of the tumors) lesions. The inclusion of extra margins of probable tumor tissue around the lesions in the training ground truth resulted in lower DSCs of 0.75 (p=0.0022). Conclusion: The model can accurately segment hepatic tumors from iUS images and has the potential to speed up the resection margin definition during surgeries and the detection of lesion in screenings by automating iUS assessment.

Anatomically informed multi-level fiber tractography for targeted virtual dissection.

OBJECTIVES: Diffusion-weighted MRI can assist preoperative planning by reconstructing the trajectory of eloquent fiber pathways, such as the corticospinal tract (CST). However, accurate reconstruction of the full extent of the CST remains challenging with existing tractography methods. We suggest a novel tractography algorithm exploiting unused fiber orientations to produce more complete and reliable results. METHODS: Our novel approach, referred to as multi-level fiber tractography (MLFT), reconstructs fiber pathways by progressively considering previously unused fiber orientations at multiple levels of tract propagation. Anatomical priors are used to minimize the number of false-positive pathways. The MLFT method was evaluated on synthetic data and in vivo data by reconstructing the CST while compared to conventional tractography approaches. RESULTS: The radial extent of MLFT reconstructions is comparable to that of probabilistic reconstruction: [Formula: see text] for the left and [Formula: see text] for the right hemisphere according to Wilcoxon test, while achieving significantly higher topography preservation compared to probabilistic tractography: [Formula: see text]. DISCUSSION: MLFT provides a novel way to reconstruct fiber pathways by adding the capability of including branching pathways in fiber tractography. Thanks to its robustness, feasible reconstruction extent and topography preservation, our approach may assist in clinical practice as well as in virtual dissection studies.

Insights from the IronTract challenge: Optimal methods for mapping brain pathways from multi-shell diffusion MRI.

Limitations in the accuracy of brain pathways reconstructed by diffusion MRI (dMRI) tractography have received considerable attention. While the technical advances spearheaded by the Human Connectome Project (HCP) led to significant improvements in dMRI data quality, it remains unclear how these data should be analyzed to maximize tractography accuracy. Over a period of two years, we have engaged the dMRI community in the IronTract Challenge, which aims to answer this question by leveraging a unique dataset. Macaque brains that have received both tracer injections and ex vivo dMRI at high spatial and angular resolution allow a comprehensive, quantitative assessment of tractography accuracy on state-of-the-art dMRI acquisition schemes. We find that, when analysis methods are carefully optimized, the HCP scheme can achieve similar accuracy as a more time-consuming, Cartesian-grid scheme. Importantly, we show that simple pre- and post-processing strategies can improve the accuracy and robustness of many tractography methods. Finally, we find that fiber configurations that go beyond crossing (e.g., fanning, branching) are the most challenging for tractography. The IronTract Challenge remains open and we hope that it can serve as a valuable validation tool for both users and developers of dMRI analysis methods.

Tracking the Corticospinal Tract in Patients With High-Grade Glioma: Clinical Evaluation of Multi-Level Fiber Tracking and Comparison to Conventional Deterministic Approaches.

While the diagnosis of high-grade glioma (HGG) is still associated with a considerably poor prognosis, neurosurgical tumor resection provides an opportunity for prolonged survival and improved quality of life for affected patients. However, successful tumor resection is dependent on a proper surgical planning to avoid surgery-induced functional deficits whilst achieving a maximum extent of resection (EOR). With diffusion magnetic resonance imaging (MRI) providing insight into individual white matter neuroanatomy, the challenge remains to disentangle that information as correctly and as completely as possible. In particular, due to the lack of sensitivity and accuracy, the clinical value of widely used diffusion tensor imaging (DTI)-based tractography is increasingly questioned. We evaluated whether the recently developed multi-level fiber tracking (MLFT) technique can improve tractography of the corticospinal tract (CST) in patients with motor-eloquent HGGs. Forty patients with therapy-naïve HGGs (mean age: 62.6 ± 13.4 years, 57.5% males) and preoperative diffusion MRI [repetition time (TR)/echo time (TE): 5000/78 ms, voxel size: 2x2x2 mm3, one volume at b=0 s/mm2, 32 volumes at b=1000 s/mm2] underwent reconstruction of the CST of the tumor-affected and unaffected hemispheres using MLFT in addition to deterministic DTI-based and deterministic constrained spherical deconvolution (CSD)-based fiber tractography. The brain stem was used as a seeding region, with a motor cortex mask serving as a target region for MLFT and a region of interest (ROI) for the other two algorithms. Application of the MLFT method substantially improved bundle reconstruction, leading to CST bundles with higher radial extent compared to the two other algorithms (delineation of CST fanning with a wider range; median radial extent for tumor-affected vs. unaffected hemisphere - DTI: 19.46° vs. 18.99°, p=0.8931; CSD: 30.54° vs. 27.63°, p=0.0546; MLFT: 81.17° vs. 74.59°, p=0.0134). In addition, reconstructions by MLFT and CSD-based tractography nearly completely included respective bundles derived from DTI-based tractography, which was however favorable for MLFT compared to CSD-based tractography (median coverage of the DTI-based CST for affected vs. unaffected hemispheres - CSD: 68.16% vs. 77.59%, p=0.0075; MLFT: 93.09% vs. 95.49%; p=0.0046). Thus, a more complete picture of the CST in patients with motor-eloquent HGGs might be achieved based on routinely acquired diffusion MRI data using MLFT.

Virtual brain grafting: Enabling whole brain parcellation in the presence of large lesions.

Brain atlases and templates are at the heart of neuroimaging analyses, for which they facilitate multimodal registration, enable group comparisons and provide anatomical reference. However, as atlas-based approaches rely on correspondence mapping between images they perform poorly in the presence of structural pathology. Whilst several strategies exist to overcome this problem, their performance is often dependent on the type, size and homogeneity of any lesions present. We therefore propose a new solution, referred to as Virtual Brain Grafting (VBG), which is a fully-automated, open-source workflow to reliably parcellate magnetic resonance imaging (MRI) datasets in the presence of a broad spectrum of focal brain pathologies, including large, bilateral, intra- and extra-axial, heterogeneous lesions with and without mass effect. The core of the VBG approach is the generation of a lesion-free T1-weighted image, which enables further image processing operations that would otherwise fail. Here we validated our solution based on Freesurfer recon-all parcellation in a group of 10 patients with heterogeneous gliomatous lesions, and a realistic synthetic cohort of glioma patients (n = 100) derived from healthy control data and patient data. We demonstrate that VBG outperforms a non-VBG approach assessed qualitatively by expert neuroradiologists and Mann-Whitney U tests to compare corresponding parcellations (real patients U(6,6) = 33, z = 2.738, P < .010, synthetic-patients U(48,48) = 2076, z = 7.336, P < .001). Results were also quantitatively evaluated by comparing mean dice scores from the synthetic-patients using one-way ANOVA (unilateral VBG = 0.894, bilateral VBG = 0.903, and non-VBG = 0.617, P < .001). Additionally, we used linear regression to show the influence of lesion volume, lesion overlap with, and distance from the Freesurfer volumes of interest, on labeling accuracy. VBG may benefit the neuroimaging community by enabling automated state-of-the-art MRI analyses in clinical populations using methods such as FreeSurfer, CAT12, SPM, Connectome Workbench, as well as structural and functional connectomics. To fully maximize its availability, VBG is provided as open software under a Mozilla 2.0 license (https://github.com/KUL-Radneuron/KUL_VBG).

Cross-scanner and cross-protocol multi-shell diffusion MRI data harmonization: Algorithms and results.

Cross-scanner and cross-protocol variability of diffusion magnetic resonance imaging (dMRI) data are known to be major obstacles in multi-site clinical studies since they limit the ability to aggregate dMRI data and derived measures. Computational algorithms that harmonize the data and minimize such variability are critical to reliably combine datasets acquired from different scanners and/or protocols, thus improving the statistical power and sensitivity of multi-site studies. Different computational approaches have been proposed to harmonize diffusion MRI data or remove scanner-specific differences. To date, these methods have mostly been developed for or evaluated on single b-value diffusion MRI data. In this work, we present the evaluation results of 19 algorithms that are developed to harmonize the cross-scanner and cross-protocol variability of multi-shell diffusion MRI using a benchmark database. The proposed algorithms rely on various signal representation approaches and computational tools, such as rotational invariant spherical harmonics, deep neural networks and hybrid biophysical and statistical approaches. The benchmark database consists of data acquired from the same subjects on two scanners with different maximum gradient strength (80 and 300 â€‹mT/m) and with two protocols. We evaluated the performance of these algorithms for mapping multi-shell diffusion MRI data across scanners and across protocols using several state-of-the-art imaging measures. The results show that data harmonization algorithms can reduce the cross-scanner and cross-protocol variabilities to a similar level as scan-rescan variability using the same scanner and protocol. In particular, the LinearRISH algorithm based on adaptive linear mapping of rotational invariant spherical harmonics features yields the lowest variability for our data in predicting the fractional anisotropy (FA), mean diffusivity (MD), mean kurtosis (MK) and the rotationally invariant spherical harmonic (RISH) features. But other algorithms, such as DIAMOND, SHResNet, DIQT, CMResNet show further improvement in harmonizing the return-to-origin probability (RTOP). The performance of different approaches provides useful guidelines on data harmonization in future multi-site studies.

Endothelial progenitor cells and plaque burden in stented coronary artery segments: an optical coherence tomography study six months after elective PCI.

BACKGROUND: Endothelial progenitor cells (EPC) are involved in neovascularization and endothelial integrity. They might be protective in atherosclerosis. Optical coherence tomography (OCT) is a precise intracoronary imaging modality that allows assessment of subintimal plaque development. We evaluated the influence of EPC on coronary plaque burden in stable disease and implemented a novel computational plaque analysis algorithm using OCT. METHODS: Forty-three patients (69.8% males, 69.6 ± 7.7 years) were investigated by OCT during re-angiography 6 months after elective stent implantation. Different subpopulations of EPCs were identified by flow cytometry according to their co-expression of antigens (CD34+, CD133+, kinase domain receptor, KDR+). An algorithm was applied to calculate the underlying total plaque burden of the stented segments from OCT images. Plaque morphology was assessed according to international consensus in OCT imaging. RESULTS: A cumulative sub-strut plaque volume of 10.87 ± 12.7 mm3 and a sub-stent plaque area of 16.23 ± 17.0 mm2 were found within the stented vessel segments with no significant differences between different stent types. All EPC subpopulations (mean of EPC levels: CD34+/CD133+: 2.66 ± 2.0%; CD34+/KDR+: 7.50 ± 5.0%; CD34+/CD133+/KDR+: 1.12 ± 1.0%) inversely correlated with the identified underlying total plaque volume and plaque area (p ≤ 0.012). CONCLUSIONS: This novel analysis algorithm allows for the first time comprehensive quantification of coronary plaque burden by OCT and illustration as spread out vessel charts. Increased EPC levels are associated with less sub-stent coronary plaque burden which adds to previous findings of their protective role in atherosclerosis.